Question #280376 on Yade changed: https://answers.launchpad.net/yade/+question/280376
mohsen posted a new comment: @#6 Jérôme I added a comment on [1] @ #7 Thanks Jan. For a betters understanding of the potato script, i run it with various 'intRadious' and the results were as expected. I also try this script with three sphere; but something strange happened and particles did not move because of gravity. can you help me why? My script: mat=CpmMat(young=1e6,neverDamage=True,sigmaT=1e20,epsCrackOnset=1e20,relDuctility=10) O.bodies.append([ # fixed: particle's position in space will not change (support) sphere(center=(0,0,0),radius=.5,fixed=True), # this two particles are free, subject to dynamics sphere((0,0,2),.5),sphere((0,0,20),.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(damping=0), ] # creates initial cohesive bonds with young=1e6 and reset intRadius O.dt = 0 O.step() is2aabb.aabbEnlargeFactor = 1. ss2d3dg.interactionDetectionFactor = 1. # change young such that new interactions will have higher young mat.young = 1e6 O.dt = .1*PWaveTimeStep() ################################################### [1]: https://answers.launchpad.net/yade/+question/266828 -- 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 : [email protected] Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp

