LIGGGHTS® - User Forum

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

How to simulate an oedometer test?

Submitted by faat1107 on Mon, 02/27/2012 - 12:03

I am a beginner of liggghts and need to make a simulation of oedometer test which is a one-dimensional test and the clay sample is immersed in water, the top wall and bottom wall allow the passage of water during consolidation. The followings are my questions :

1. Should there be two types of particles , one is the clay particle and the other one is water particle?
2. How can I create the wall which allow the passage of water particles but not for clay particles?
3. Could anyone send me example of the script of simulation of oedometer test?

Very Low efficiency of simulation + warning of Maxmimum number of particle-tri neighbors...

Submitted by venes520 on Sun, 02/26/2012 - 16:37

Hi guys,

I am now still testing the packing of particles and I used the packing example.
Differently, I used the pair_style of gran/hooke/history/stiffness 1 0 and tried to used slightly smaller particle size (I tried 3mm but failed), I also import mesh from stl files generated by Gmsh.
My problem is:

1. the Pour/dev function actually insert much less particles into the region, (less than 100), but the region is big enough to contain a lot more particles.

Unexpected Crash when integration (maybe)

Submitted by venes520 on Fri, 02/17/2012 - 15:59

Hi guys,

Base on the example of in.packing, I am trying to do a packing similarly.
Different from the example, I use pair style of gran/hooke/history/stiffness instead. I don't think it would be problematic but the system crash down when it came to the nve integration,
Here is the output of the screen:

LIGGGHTS 1.4.4 based on lammps-10Mar10
Created orthogonal box = (-1 -1 -1) to (2 2 2)
1 by 1 by 1 processor grid
0 atoms in group upper
0 atoms in group lower
0 atoms in group total

How to move an ideal wall ?

Submitted by venes520 on Thu, 02/16/2012 - 10:50

I created ideal walls and want to use the fix wall/gran/hertz/history/stiffness command and the syntax of this command is:

fix ID group-ID wall/gran/(*) friction_damp_flag cohesionflag wallstyle args keyword values ...

I don't really understand the group-ID here, does it mean the ID of the group that I want to create walls to bound? or the ID for the group of the walls themselves? How to "group" these walls? they don't consist of particles.

Gravity and Archimedes

Submitted by tdl on Wed, 02/15/2012 - 13:36

Noob question.
I would like to simulate particles falling down in a container filled with water. For now I do not need to couple with CFD, but I would like to able to tell the simulation that once the particles reach a given region, they are subject to an Archimedean lift.
But I am not sure how I can do it? Any simple example that I may have missed?
Thanks!

How to define a container & void ratio

Submitted by venes520 on Wed, 02/15/2012 - 08:03

Hi guys,

I am trying to simulate a Direct Shear Box Test for soil.

I want to create two containers (two granular boxes, upper & lower) to contain the soil assemble and do the shearing.

To create the containers, I am thinking of wall/gran, however, I don't know how to group this walls together so that they become a container and move together, be able to read the force on it.

Please help me out if you have any thoughts or did something similar.

Remarks on Benchmark 3

Submitted by SamWong on Wed, 02/15/2012 - 07:02

Hi there,

Can someone comment on one of the results in Benchmark 3 where the speed of LIGGGHTS is slower than the commercial code for the continuous blending mixer case? Why is it slower in this test case?

Also, which commercial code is it? Can you please email/message me privately? I'm doing some benchmark myself too.

Thanks,
Sam

invalid atom style

Submitted by neelu09 on Tue, 02/14/2012 - 20:09

hi

I am not able to use atom style ellipsoid for my simulation. I have done work using liggghts with granular atom style and had no problems. But when i try to use ellipsoid atom style, it gives an error saying "invalid atom style". Same goes for colloid and dipole atom styles as well. Here's my input file
atom_style ellipsoid
boundary f f f
newton off
atom_modify sort 0 0

communicate single vel yes
units si

region reg cylinder y 0 0 0.0225 0 .08 units box

create_box 2 reg

neighbor 0.001 bin
neigh_modify delay 0

shape 1 1.0 1.0 1.0
shape 2 1.0 4.0 1.0

Pages

Subscribe to RSS - LIGGGHTS® - User Forum