← Back to team overview

yade-users team mailing list archive

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

 

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
Hi Mohsen,
I did some changes recently, will have a closer look
cheers
Jan


2016-05-31 15:37 GMT+02:00 mohsen <question294657@xxxxxxxxxxxxxxxxxxxxx>:

> 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.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users
> Post to     : yade-users@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-users
> More help   : https://help.launchpad.net/ListHelp
>

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