yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #03376
Strange behavior of periodic cell
Hello,
I'm trying to do some simulations using periodic cell. In one case, I prescribed some extension of the periodic cell, but the particles in the simulation behave really strange - the dispacement of particles si continuously increasing.
So if someone ran the attached script and said what is wrong in it, I would be very grateful :-)
cheers
Jan
from yade import pack, log, qt
# create material
O.materials.append(CpmMat(young=1e7,epsCrackOnset=20,relDuctility=3,sigmaT=3e30,frictionAngle=.5,G_over_E=.2))
O.periodic=True
# create 3x3x3 particles
a=1.; r=a/6.
for i in range(3):
for j in range(3):
for k in range(3):
O.bodies.append(utils.sphere([(2*i-1)*r,(2*j-1)*r,(2*k-1)*r],r))
O.cell.refSize=(a,a,a)
O.engines=[
BoundDispatcher([Bo1_Sphere_Aabb()]),
InsertionSortCollider(),
InteractionDispatchers(
[Ig2_Sphere_Sphere_Dem3DofGeom()],
[Ip2_CpmMat_CpmMat_CpmPhys()],[Law2_Dem3DofGeom_CpmPhys_Cpm()]),
NewtonIntegrator()
]
# compression of periodic cell
O.cell.trsf+=Matrix3(-.1,0,0, 0,0,0, 0,0,0)
v=qt.View()
v.viewDir=Vector3(0,0,1)
v.center()
O.run()
Follow ups