New question #695791 on Yade:
https://answers.launchpad.net/yade/+question/695791
Dear All,
I have a problem in simulating spherical particles' gravity deposition in 2D. I
use state.blockedDOFs="yXZ" to fix all the particles in y plane. But I found
that
1) all the particels fall together (i.e. free falling and no particles
collision),
2) the wall disappears (does not prevent particles go through),
3) no packing occurs.
I modified the tutorial example 02-gravity-deposition.py into my input file. It
reads:
from yade import pack, plot
O.bodies.append(geom.facetBox((500,100,100),(500,100,100),wallMask=31))
sp=pack.SpherePack()
sp.makeCloud((0,100,0),(1000,100,200),rMean=15,rRelFuzz=.3)
sp.toSimulation()
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
InteractionLoop(
# handle sphere+sphere and facet+sphere collisions
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
),
NewtonIntegrator(gravity=(0,0,-9.81),damping=0.4),
PyRunner(command='checkUnbalanced()',realPeriod=2),
PyRunner(command='addPlotData()',iterPeriod=100)
]
O.dt=.05*PWaveTimeStep()
#Block the "flat" axis (Y), and only allow rotations around it.
for i in O.bodies:
i.state.blockedDOFs="yXZ"
O.trackEnergy=True
def checkUnbalanced():
if unbalancedForce()<.005:
O.pause()
plot.saveDataTxt('bbb.txt.bz2')
# collect history of data which will be plotted
def addPlotData():
plot.addData(i=O.iter,unbalanced=unbalancedForce(),**O.energy)
plot.plots={'i':('unbalanced',None,O.energy.keys)}
plot.plot()
O.saveTmp()
Is there any problme in my input file ?
Many thanks for your help.
Best regards,
Man Yick
--
You received this question notification because your team yade-users is
an answer contact for Yade.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
