LIGGGHTS® - User Forum

Hi

I am trying to develop a new version of fix_pour by combining the "regionExempts" function of pour/dev with pour. I have created a new pour called pour/exempt and copied accross all the relavent bits of code. I have also made the necessary links to these other files as well: mech_param_gran.cpp, region_block.h (and cylinder), fix_gravity.h, pair_gran.cpp and style_fix.cpp.

Hi all,

I am using compute pair/gran/local and the corresponding dump to look at the particle-particle contact data (and heat flow), but also need to look at particle-wall interactions. As far as I can see, these are not included in these computes/dumps?

I think this would be a very useful addition - I guess it would require assigning each wall an ID etc. In particular I need to identify particles in contact with the wall, and sum their heat fluxes etc.

Dear friends,
Could some body please let me know how can I apply rotation and translation (at the same time) to a moving mesh imported from a stl file. Any response will be truly appreciated.
Regards,
Carlos

Hi all,

I'm trying to set up a sim with a region of fixed temperature particles at the top and bottom, and particles that have variable temperature in between. The docs suggest I do not include adiabatic particles in the heat/gran fix. So I've set up a lattice of (4200) particles and defined regions as:

region upper block INF INF INF INF ${upper_limit} INF units box
region lower block INF INF INF INF INF ${lower_limit} units box
region middle union 2 upper lower side out

group thermal_group region middle
3000 atoms in group thermal_group

Hi everybody,
I would like to start the shear test using LIGGGHTS.
Could you help me in doing it.
I was using velocity command for the control of temperature. So far no success. I came to know about the shear test from the following message.
Many thanks in advance.
Samantaray

Tue, 11/30/2010 - 17:18 — tshan

Dear all,

I have already done some triaxial tests and simple shear tests, but I want to know the contact normal among the particles and then realise them in paraview to further generate the force-chain.

Hi Chris,
Could you give hints about the calculation of neighbors in LIGGGHTS?
I am using the following commands in the programme.
lattice fcc 0.0013
neighbor 0.00125 bin

It is creating
Total # of neighbors = 48418
Ave neighs/atom = 24.209
Neighbor list builds = 27

Hi

i need create 2 atoms with different youngsModulus parametr.

i try this:

fix m1 all property/global youngsModulus peratomtype 5.e7 5.e6

create_atoms 1 single 0.0 0.0 10.0
set atom 1 density 2000 diameter 2.0

create_atoms 2 single 5 0.0 10.0
set atom 2 density 2000 diameter 2.0

So i get that ERROR: Invalid atom type in create_atoms command
What's wrong?

Thanks
jikra

Hi, all

I just got an error with the 'fix' when running the LIGGGHTS examples:

fix m1 all property/global youngsModulus peratomtype 5.e6
ERROR: Invalid fix style

It seems the property/global is a LIGGGHTS style for the 'fix'.
The LAMMPS examples were done well, without any error, perhaps the LIGGGHTS was not invoked properly??

Any help is appreciated.

Hi everyone,

I'm trying to set up a simulation with 3 spheres clumped together, exactly as seen in LIGGGHTS Features > Non-spherical particles > Example - pouring a continuous stream of non-spherical particles, where can be seen particles made of 3 spheres glued together.

I'm exploring some possibilities to do a modification of this study http://www.nt.ntnu.no/users/skoge/prost/proceedings/ecce6_sep07/upload/1... by simulating a packed bed with non-spherical particle as my PhD study. My idea is starting with 3 spheres glued and then move on with different particle geometries.

Hi all,

I have so far been running LIGGGHTS with more-or-less "large" particles (mm/cm), but now want to move towards the size that interests me, around ~1 µm or so. However, if I try to use fix particletemplate/sphere and fix particledistribution/discrete to generate small particles I get an error which seems to arise from a hard-coded limit:

if(volume_expect < 1e-12) error->all("Fix template/sphere: Volume expectancy too small");

in fix_template_sphere.cpp.