yade-users team mailing list archive

[mohsen <question280376@xxxxxxxxxxxxxxxxxxxxx>]

Question #280376 on Yade changed:
https://answers.launchpad.net/yade/+question/280376

mohsen posted a new comment:
@#6
Jérôme I added a comment on [1]

@ #7
Thanks Jan.
For a betters understanding of the potato script, i run it with various  'intRadious' and the results were as expected. I also try this script with three sphere; but something strange happened and particles did not move because of gravity. can you help me why?
 My script:
mat=CpmMat(young=1e6,neverDamage=True,sigmaT=1e20,epsCrackOnset=1e20,relDuctility=10)

O.bodies.append([
# fixed: particle's position in space will not change (support)
sphere(center=(0,0,0),radius=.5,fixed=True),
# this two particles are free, subject to dynamics
sphere((0,0,2),.5),sphere((0,0,20),.5)
])

intRadius = 1

O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=intRadius,label='is2aabb')]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=intRadius,label='ss2d3dg')],
[Ip2_CpmMat_CpmMat_CpmPhys()],
[Law2_ScGeom_CpmPhys_Cpm()]
),
NewtonIntegrator(damping=0),
]

# creates initial cohesive bonds with young=1e6 and reset intRadius
O.dt = 0
O.step()
is2aabb.aabbEnlargeFactor = 1.
ss2d3dg.interactionDetectionFactor = 1.

# change young such that new interactions will have higher young
mat.young = 1e6
O.dt = .1*PWaveTimeStep()
###################################################
[1]:  https://answers.launchpad.net/yade/+question/266828

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.