ESyS-Particle

Particle packing with equal radius

Asked by dennis

Hi to all,

 I'm just wondering it is possible to construct a block of cohesionless
particles with equal radius in Esys-particle. If it is possible then how it
is done? The RandomBoxPacker will just generate particles with Rmin to Rmax. On the other hand,
in the InteractionGroup routine, when scaling is set to "True" does
this mean that all the particles in that group will have the same stiffness regardless of the particle radius?

Sincerely,
Dennis

Question information

Language:
English Edit question
Status:
Answered
For:
ESyS-Particle Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
SteffenAbe (s-abe) said : 		#1

Hi Dennis,

just setting Rmin=Rmax in the RandomBoxPacker should work IIRC. At least it does in gengeo which uses a similar implementation of the particle packer.

Steffen

Revision history for this message
dennis (dca1101) said : 		#2

Hi Steffen,

 Thank you for your reply. It does work in the RandomBoxPacker routine. I tried the example
in the gengeo, "gengeo_example1.py", using the InsertGenerator2D routine, when Rmin=Rmax, it wont work. No
particles were inserted in this case. Is there a work around in this? or I miss some parameters to change make it
work.

Sincerely,

Dennis

Revision history for this message
Dion Weatherley (d-weatherley) said : 		#3

Hi Dennis and Steffen,

I tried gengeo_example1.py with Rmin=Rmax and confirmed that it does not work (no particles inserted). However it does work using the 3D example (gengeo_example2.py). Looks like you may have triggered a bug in GenGeo. We'll look into this in more detail.

Another option for in-simulation geometry creation is to use the CubicBoxPacker in place of RandomBoxPacker. CubicBoxPacker generates a face-centred cubic array of equal-sized particles. The API docsfor it are here:
http://esys.esscc.uq.edu.au/esys-particle_python_doc/snapshot/pythonapi/html/esys.lsm.geometry.GeometryPy.CubicBoxPacker-class.html

I'm not sure whether/how well this works in 2D but worth a try...

Cheers,

Dion

Revision history for this message
Toketo (tomcarlosimborio) said : 		#4

Hmmm. I tried to have the minRadius and maxRadius = 1.0 in the gengeo_example2.py and it wouldn't work.

Can you help with this problem?

Provide an answer of your own, or ask dennis for more information if necessary.

To post a message you must log in.