LIGGGHTS® - User Forum

Hi,

I'm looking to calculate the dissipated energy for colliding spheres. In the documentation, the computeDissipatedEnergy option for a granular model references a fix style to use later that throws up an error when I tried to use it - so I am unable to access these values. In a previous forum post I read it seems that this is a premium feature (I am only using the public version). My question is, is there a way in LIGGGHTS-public to efficiently track the energy dissipation?

Hello !

I've been using LIGGGHTS for around a year now mainly with spherical particles. I want to create a '2D' model using cylindrical particles, it would be a 3D plan symmetrical model. To do so I created superquadric particles with N1=8 and N2=2. But whenever I try to use the command 'write_data' to save state the program crashes. I get the following error :

Hi
I am currently using the 'LIGGGHTS flexible-fibers' version. As far as I can see there is option for elastic bonds. Is there a way to introduce elasto-plasticity in the bonds.

Thanks

Hi everyone,
I am a new user to LIGGGHTS and I am trying to simulate a packing including both non-cohesive and cohesive particles and the interaction between them. I already tried to simulate the in.cohesion example in LIGGGHTS/ Tutorial_public but when I tried to add a larger size particle (non-cohesive) together with the existing cohesive particles, I saw no cohesive behaviour between them. Is there anybody has any ideas on how to insert a packing of non-cohesive and cohesive particles? I really really appreciate any helps.

I have been able to use the granular model effectively to simulate monodisperse and binary mixtures of spherical particles. I especially like that I can use read_data to set the initial position/velocity/etc of each individual particle. I now would like to simulate mixtures of particles with more complex shapes - e.g. a mixture of cylindrical particles.

Hello all,

I am trying to simulate the isotropic compression of a spherical wall/bounding box along the negative 'r' direction in spherical coordinates. Can anyone direct me towards the approach to edit mesh_module_stress_servo.cpp(/h) and/or fix_deform.cpp(/h) so that these modules work in spherical coordinates? Thank you very much!

Hi
I am using servo wall to maintain constant force on wall. And I have some doubts.

1. The formula referred in documentation for PID controller
output = kp * error + ki * errorsum + kd * errorchange
I want to print or save the output, error, errorsum and errorchange values. How can I do this in LIGGGHTS.

2. Is there any thumb rule for selecting kp, ki, and kd values.

3. How can I change a variable value in every step.

Hi
I am new to LIGGGHTS and modeling an impact problem. As shown in attached figure, impact is happening between a web and fixed wall. I want to compute the total impact force with time (sum of contact force for all contact points). Please note that the contact points change as the simulation progresses. Please help me by suggesting a suitable command for that.

Thanks

Hi

it's regarding storing the variables at particular timestep,
Then need to use those variables in the next loop

The script is like this

#####################################################
start of the program
>> Scripts intitiates
>> define varaibles
>> variables are changing at every timestep

>> ran for ten timestep
run 10
unfix the all the fixes
>> save the varables at 10

Hi,

i m trying to run my first simple simulation but an error comes

ERROR: Mesh error: dimensions too small - use different unit system (../multi_node_mesh_parallel_I.h:584)

I actually used SI unit and stl. file in ASCII file.

I can post some codes of mein and look forward to get some advices:

#Definition of boundaries
variable xmin equal -0.0102
variable xmax equal 0.0102

variable ymin equal -0.042
variable ymax equal 0.042

variable zmin equal -0.042
variable zmax equal 0.042