yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #06569
[Question #212185]: About gravity deposition
New question #212185 on Yade:
https://answers.launchpad.net/yade/+question/212185
Hi,
I have two questions about generating a package using gravity deposition:
1. How can I get a particle package with prescribed after gravity deposition?
I have tried using the following script, but there will be a large compression due to gravity deposition in the package.
##SCRIPT
from yade import pack
num_spheres = 30000
rr = 0.2
compFricDegree = 30
finalFricDegree = 30
damp=0.7
stabilityThreshold=0.01
young=5e6
mn,mx=Vector3(0,0,0),Vector3(20,80,50)
thick = 0.01
O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=radians(compFricDegree),density=2600,label='spheres'))
O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=0,density=0,label='walls'))
walls=utils.aabbWalls([mn,mx],thickness=thick,material='walls')
wallIds=O.bodies.append(walls)
sp=pack.SpherePack()
sp.makeCloud(mn,mx,rRelFuzz=rr,num=num_spheres,porosity=0.35)
O.bodies.append([utils.sphere(center,rad,material='spheres') for center,rad in sp])
O.dt=.5*utils.PWaveTimeStep()
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
),
GravityEngine(gravity=(0,0,-9.81)),
NewtonIntegrator(damping=0.4),
]
O.run(2000, True)
while 1:
O.run(500, True)
unb=unbalancedForce()
print 'unbalanced force:',unb
if unb<0.1:
break
yade.qt.Controller(), yade.qt.View()
How should I modify my script if I want to:
(1) generate a sample which can fill the whole box; and
(2) has the prescribed porosity, i.e. 0.35
2. How can I get the stress the particles apply on the boundaries after the gravity deposition finishes?
Many thanks.
Nait
--
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.