yade-users team mailing list archive

[Marcus Moravia <question661926@xxxxxxxxxxxxxxxxxxxxx>]

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.