LIGGGHTS® - Developer Forum

Topics related to developing with LIGGGHTS® can be discussed here: discussion about implementation details, C++, MPI and debugging tools

Calculating force in a bonded simulation

Submitted by felipe.merino on Tue, 05/09/2017 - 22:05

Hello people. I'm working on a simulation of a tensile test with granular bonds. I'm currently using liggghts 3.5.0 with the bondspackage uploaded by richtie. The thing is that I'm trying to get a tension-deformation diagram for the case, but the force output by the command compute/reduce is just too low. The bonds are working as intended, as they break when they should break (the simulation uses a velocity in one end, while the other end is freezed, so the deformation is directly proportional to the time elapsed), but as i said, the force is orders of magnitude wrong.

Asymmetric stress tensor calculation for non-spherical particles

Submitted by vinaym on Mon, 04/24/2017 - 16:23

Hello all,

I want to export all nine components (averaged) for particle (asymmetric) stress tensor with separate kinetic and static contributions. LIGGGHTS at the moment exports 6 components (symmetric) stress tensor.

has anyone implemented this already? Clues on implementation are welcome.

With a quick look it seems fix_ave_euler and compute_stress_atom is the right place to start. Can anyone please confirm?

Thanks and kind regards,
Vinay

Order of force computes

Submitted by mschramm on Tue, 04/18/2017 - 21:49

Hello,
Could someone guide me on where I should look for where each force compute is done? (Contact, bond, other...)
I am working on the gran/bond model and would like to turn off contact forces between bonded particles (since they are part of a singular geometry).
I tried using the grep command to attempt to find where it may be but had no luck.

Thank you!

individual material properties

Submitted by jordantcox on Thu, 04/13/2017 - 10:37

Heya,

I am also working on trying to have individual material properties for each discreet element. I am, for example. working on spheres that would each have a different youngs modulus in a packed bed. I think this would be best accomplished through a restart file, where I run a simulation, and then modify the restart file (like I am doing right now to provide an individual radius based on an outside algorithm) and then to read in that restart file.

Radiation model can be added for packed bed

Submitted by wuhaoinet on Fri, 04/07/2017 - 08:11

I work for the particle radiation in dense granular system. A useful radiation model has been tested in LIGGGHTS, and it is also can be applied in CFD-DEM.
I am not permitted to open my code. But it is very easy to understand the numerical model from my paper:
1.Wu, H., Gui, N., Yang, X., Tu, J. and Jiang, S., 2016. Effect of scale on the modeling of radiation heat transfer in packed pebble beds. International Journal of Heat and Mass Transfer, 101, pp.562-569.(based on LIGGGHTS)

[SOLVED] Index handling Particle Bond Model (Potyondy) via cohesion

Submitted by Baric on Wed, 03/08/2017 - 13:25

Hello everybody,

I am currently writing a rigid bond contact model (as from "A bonded-particle model for rock" Potyondy and Cundall 2004) using the contact style cohesion. I have everything working, but I came across an odd behavior:

Pages

Subscribe to RSS - LIGGGHTS® - Developer Forum