LIGGGHTS® - User Forum

LIGGGHTS® related topics can be discussed here: discussion about models, installation, feature requests and general discussion

Error libmpi_cxx.so.0

Submitted by katterfe on Thu, 01/20/2011 - 18:28

Hey everybody,

after installing Liggghts on Fedora 14 the following error message occurs when I try to start liggghts:

liggghts: error while loading shared libraries: libmpi_cxx.so.0: cannot open shared object file: No such file or directory

I have installed the newest version of openmpi and the mentioned file can be found in /usr/lib/openmpi/lib

I attach a screenshot.

Thanks for your help!
Andre

Question about the addfroce command and the size of the box

JF's picture
Submitted by JF on Tue, 01/18/2011 - 17:31

Hi all,

I have two questions:

Do you know how should I use the command (fix addforce) to push a group of spheres ?
Because this command only works during an energy minimization, but I don't understand the minimize command (it is a specific command for atomic simulation).

During a tensile test on a granular sample, it is possible that some spheres are outside the box initially defined. Do I need to change the size of the box during the simulation or there will be no errors ?

Regards
JF

The rotation viscosity

Submitted by tshan on Tue, 01/18/2011 - 07:56

Hi, all

I am thinking and finding if there is a rotation viscosity(damping) or rotation constraint command in LIGGGHTS, which may allow realistic description of granular flow behavior.

I have found a fix momentum command which allow me to zero angular momentum every several steps, but it is not what I want.

Any response from you will help! Thanks a lot!

Best,
Tong

Specify region- or group- dependant coefficients of friction

Submitted by alexander.polson on Fri, 01/14/2011 - 23:34

Hi Everybody

I wanted to set the friction to different values, depending on where the particle is in the model, e.g. in a hopper, and in the process below the hopper.

My approach was to:
- define a region, REG_hopper
- associate a group of particles with this region
- set the friction first for "all", and then for the group using "fix property/global"

The steps below shows only the applicable commands:
region REG_hopper block x_lo x_hi ... z_lo z_hi units box
...
group GRP_hopper region REG_hopper
...

Create txt files from Pizza.py

JF's picture
Submitted by JF on Thu, 01/13/2011 - 18:16

Hello everyone,

As it is my first question of the year 2011, I take this opportunity to wish all LIGGGHTS' users a Happy New Year.

About my question, I would like to know if it exists a Pizza's function to export the data from the custom dump to txt files ?

Each line of txt file will be the numerical data dumped for one atom. I need these txt files to perform posttreatment with Matlab.

Regards
JF

Illegal if command

Submitted by tshan on Wed, 01/12/2011 - 08:18

Dear Christoph,

There's a bug (Illegal if command) in early version of lammps, and it has been fixed in current lammps. However, the bug is still in liggghts. I used the patch file to liggghts that applied in lammps and it didn't work.

So could you help me to fix this? The patch file of lammps has been attached.

Many thanks!

Best,
Tong

Error unfixing a wall/gran/hertz/history with a mesh/gran

diego.peinado's picture
Submitted by diego.peinado on Thu, 12/30/2010 - 18:04

Hello Christoph.
I was playing around with liggghts and I tried to first:
fix container all mesh/gran silo.stl 1 1.0 0 0 0 0 0 0
fix container_wall all wall/gran/hertz/history 1 0 mesh/gran 1 container

Then after some running steps I tried to:

unfix container_wall
run 1000
But I have a segmentation fault in the first iteration of the run 1000 command.
I've debugged a bit, and found that the fix is deleted, but when running the deleted fix is invoked again in Verlet::run()
if (n_pre_force) modify->pre_force(vflag);

different diameters of beads

raguelmoon's picture
Submitted by raguelmoon on Thu, 12/30/2010 - 14:54

Hi,
I want to pour two different diameters particles in simulation box, I used following script but couldn't put different diameters but just single, print-out is also given below:
....
.....

fix pts1 all particletemplate/sphere 786 atom_type 1 density constant 2460 radius constant 0.0020
fix pts2 all particletemplate/sphere 786 atom_type 1 density constant 2460 radius constant 0.003
fix pdd1 all particledistribution/discrete 6778 1 pts1 0.5 pts2 0.5
.....
......

region bc cylinder z 0 0 0.088 0.251 0.4 units box #bc block 0.0 0.239 0.0 0.239 1.2 1.5 units box

Pages

Subscribe to RSS - LIGGGHTS® - User Forum