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()
_______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
