yade-users team mailing list archive

[<_honzik@xxxxxxxxxx>]

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()