← Back to team overview

yade-users team mailing list archive

[Question #294657]: Cpm model problem in Yade 1.20.0-151-0ecfcd0~trusty

 

New question #294657 on Yade:
https://answers.launchpad.net/yade/+question/294657

In the name of god
Hi all
3 days ago(28.5.2016) i installed new version of yadedaily in a sysmtem (Yade 1.20.0-151-0ecfcd0~trusty). When runnug my code i underestood there is a problem with Cpm model.
to check it i wrote a simple script (3 bouncing sphere). I was wondered by result; simulation results of the same script on diffrenet versions, was different. in 'Yade 1.20.0-126-e3271e8~trusty' , the spheres bounces naturally whereas in 'Yade 1.20.0-151-0ecfcd0~trusty' version the spheres path fixed one!!!! I also checked another script in different versions and the results differed!
What is the reason?

My script:
#####################################
mat=CpmMat(young=1e6,neverDamage=True,sigmaT=1e20,epsCrackOnset=1e20,relDuctility=10)
O.materials.append(mat)
O.bodies.append([
# fixed: particle's position in space will not change (support)
sphere(center=(0,0,0),radius=.5,fixed=True),
# this particles is free, subject to dynamics
sphere((0,0,2),.5),sphere((0,0,5),.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(gravity=(0,0,-10),damping=0),
]


O.dt = 0
O.step()
is2aabb.aabbEnlargeFactor = 1.
ss2d3dg.interactionDetectionFactor = 1.


mat.young = 1e6
O.dt = .000001*PWaveTimeStep()


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