Question #661926 on Yade changed:
https://answers.launchpad.net/yade/+question/661926
Marcus Moravia gave more information on the question:
Hello again,
I have done some scripting checks. Although the source code contemplates
homothetic displacements (Shop_02.cpp), using only two clumps with the
same geometry in the tutorial script, I was not able to observe this
behavior (the code results different overlapping values between the
initial condition and after the growth process of the spheres).
Moreover, in the simplified script below that consider only the
growParticles function (i.e., function used in the
TriaxialStressController function), the radius of the spheres that make
up the clumps seems to be updated to a value that follows:
new_radius=(multiplier^number_of_clumps)*initial_radius.
Thank you for your help.
### Simplified Script ###
O.materials.append(FrictMat(poisson=0.5,density=2600,label='sphereMat'))
from yade import pack
nElem = 7
c=pack.SpherePack([((0,0,0),.003),((.006,0,0),.003)])
sp=pack.SpherePack()
print 'Generated # of
clumps:',sp.makeClumpCloud((0,0,0),(10,10,10),[c],num=nElem,periodic=False,seed=1)
sp.toSimulation(material='sphereMat')
for c in O.bodies:
if c.isClump:
O.bodies.updateClumpProperties(discretization=250)
growParticles(multiplier=2)
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
),
NewtonIntegrator(damping=.4,gravity=(0,0,-9.81)),
# PyRunner(command='volumeClumps()',iterPeriod=50,label='volume')
]
O.dt=.5*PWaveTimeStep()
O.saveTmp()
O.step()
yade.qt.View()
--
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
