LIGGGHTS® - User Forum

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

Delete particles

Submitted by RobertG on Mon, 11/10/2014 - 16:21

Hello,
im modeling an continous mixer and have two different particle types.
I continously inserte particles at two different location and the mixed up particles leaving the mixer at the same place.
Now the problem. Because of the continously use, there are build more and more particles.
Is it possible to delete all particles, which are at the same area (exit of the mixer).
I have tried it with the delate_atome command and wasn't very successfull.

If anybody can help me. I would be very grateful.

Best regards
RobertG

Restitution coefficient of rigid bodys

Submitted by Nafets on Tue, 11/04/2014 - 09:37

Hi all,
I generate a rigid body shaped like a box with moltemplate and simulate a droptest. The problem is that the body jump's after contact with the mesh to a height that is not in relation to the physical law. k = square root (h'/h)
Restitution coefficient: 0,5
Previous height: 0,5 m
Height after first contact: 0,27756295 m
According to the formula the body should just jump to a height of 0,125m.

Particle moving through .stl-walls

Submitted by RobertG on Fri, 10/31/2014 - 12:33

Hello again,
I'm working with a program, based on the tutorial Continuous Blending Mixer.
My Problem is, that at the insertion area, the particle start to "explode" at a certain time.
That means, some of the particles are going through the wall of the .stl-file.
The starting velocity of the particles is high enough to avoid, that some of the old ones are still in the area, when the new ones are created.
I'm shure, that the area they are created is not too small and there is plenty of space between the the area they are created and the wall.

Is the shear modulus equation in the doc wrong?

Submitted by xuji2010 on Fri, 10/31/2014 - 09:47

Hi all,

In the doc of DEM simulation, the effective shear modulus is:
1 2(2+v1)(1-v1) 2(2+v2)(1-v2)
-- = ------------------ + -------------------
G Y1 Y2

but, the in the global_properties.cpp, line 217:
matrix->data[i][j] = 1./(2.*(2.-vi)*(1.+vi)/Yi+2.*(2.-vj)*(1.+vj)/Yj);

The above two don't match, which one is wrong?

Appreciate any help in advance!

Problem with two or more force fields in a single input script

Submitted by b.saprativ on Wed, 10/29/2014 - 07:04

I am trying to incorporate two separate force fields for two separate block regions in a single input script as I need to giving force in two separate regions

I am using the following script for one block

#zone 1
region spray block -0.0295 0.0555 0.2 0.24 -0.016 0.014 units box
fix push all addforce 0.0 0.1 0.0 region spray

I tried to add the same set of commands for the other for the other zone but its giving error

#zone 2
region spray block -0.01 0.05 0.3 0.4 -0.06 0.01 units box
fix push all addforce 0.0 0.1 0.0 region spray

Pages

Subscribe to RSS - LIGGGHTS® - User Forum