LIGGGHTS® - User Forum

Hello Christoph,

I wonder which commands can complete the reduction of particle diameter over time in a specified region just as you have gived example lammps_particleReduction_.wmv in this web. I think it is maybe lammps's some command or tips, Could you kindly give me help or instruction?

Thank you.

best regards,

leexz

Hi,

When I use
pair_style peri/lps for the pair style, I get this error:
ERROR: Invalid atom style

I wonder if it is possible to use this pair style with LIGGGHTS or not.

Thanks,

Ozgur Keles

Hi LIGGGHTS Users,

we try to create a cone-style insert region (just like a hopper) with the pour/dev/packing command. The input script looks like this:

#particle insertion region
region ins_cone cone z 0.04 0.0 0.01 0.065 0.09 0.16 units box
region ins_cyl cylinder z 0.04 0.0 0.065 0.08 0.16 units box
region no_ins_region intersect 2 ins_cyl ins_cone side out units box

...

#region and insertion
fix ins all pour/dev/packing 1 distributiontemplate pdd vol ${alphastart} 10 overlapcheck no region ins_cyl regionexempts 1 no_ins_region

Dear LIGGGHTS users!

Attached you find a (not very sophisticated) Python script to produce a data file as input for a new simulation generated from a specified time step of a dump file of another calculation. Maybe this is useful for somebody as a starting point for similar problems.
The dump/custom command in the first simulation includes "id type type x y z ... radius mass".

Christian

PS: You have to rename the file extension to .py

Hi everybody,

I am interested in a simulation of a big sphere containing a many smaller particles. I tried this simulation by adding particles of a defined radius on a regular lattice. The problem is that this reguar structure works a a kind of crystalline system with slip planes when I try to deform it.

However I'd like to try another approach by adding the particles without a lattice with the "pour command". Is there a way to use this command with a sphere as the input region (and not a cube or a cylinder)?

regards,
Martin Koester

Hi all,

For the gran/hooke pair style the documentation says that the equation for the elastic constant is:

k_n = 15/16* ...

In the code the elastic constant is calculated as follows:

kn=(16./15.)*sqrt(reff)* ...

So which one is right, the code or the documentation?

Thanks
Yansan

Hi LIGGGHTS users,
In new fix pour/dev, there is no pouring after i continue my simulation with restart file. Can anybody solve my problem?

Thanks in advance...
Ram

Hi,
Could some body please, let me know what causes the following error message:
*****************************************************************
Setting up run ...
Memory usage per processor = 10.1914 Mbytes
Step Atoms KinEng 1 Volume
1 2000 0.0047170791 0 0.0015
......
.........
37000 2000 0.0094914152 0 0.0015
38000 2000 0.0094929527 2.9568818e-05 0.0015
39000 2000 0.0097540535 0.00030106194 0.0015
40000 2000 0.0101535 0.00085611282 0.0015

Hi all,

Christoph,

Is it possible to define multiple materials in a simulation, and if so, do you have any example scripts that do this? I haven't been able to figure this out yet. For example, I want to be able to assign different friction properties for particle-particle and particle-wall interactions.

Thanks,
Chris