Strange behavior: sand pile settles in a squared shape ?
Hi all,
When simulating the settlement of a relatively dense) packing on a plane surface, I obtain an equilibrium state which seems a bit awkward: from above, the particles tend to spread in a squared shape instead of a circular shape (as I would have expected). I tested different contact laws as well as different initial packing shapes and it seems to occur systematically (the squared shape is more or less evident depending on the case). The initial packing is built based on randomDensePack(): could this be the reason?
Please find below a MWE if you want to reproduce the "problem" (jut click and play for ~10 000 iterations when the GUI opens).
Any feedback would be appreciated.
Luc
---
# -*- coding: utf-8 -*-
from yade import pack
#### assembly
O.materials.
L=2.
D=L/2.
sphereRad=
pred=pack.
#pred=pack.
O.bodies.
dim=utils.
xinf=dim[0][0]
xsup=dim[1][0]
X=xsup-xinf
yinf=dim[0][1]
ysup=dim[1][1]
Y=ysup-yinf
zinf=dim[0][2]
zsup=dim[1][2]
Z=zsup-zinf
#### get spheres dimensions (pre-processing)
R=0
numSpheres=0.
for o in O.bodies:
if isinstance(
o.shape.
numSpheres+=1
R+=o.
Rmean=R/numSpheres
print('nb of spheres=
#### floor
O.materials.
O.bodies.
##### engines
O.engines=[
ForceResetter()
,InsertionS
,Interactio
[Ig2_
[Ip2_
[Law2_
)
,GlobalStif
,NewtonInte
]
#### open yade GI
from yade import qt
v=qt.Controller()
v=qt.View()
#### Allows to reload the simulation
O.saveTmp()
print('press play now!)')
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Luc Scholtès
- Solved:
- Last query:
- Last reply: