Hi (Anton, Chiara?) Who is working on dynamic problems? Could not you explain why it occurs (see attached script)? And do test the configuration with own laws?
-- Best regards, Sergei Dorofeenko, Cand. Phys.-Mat. Sci. Institute of Problems of Chemical Physics, Chernogolovka, Moscow region, Russia. e-mail: [email protected]
# encoding: utf-8 from yade import utils,geom Rsph=0.003 density=2400 frictionAngle=radians(35) tc=0.0001 en=0.3 es=1.0 params=utils.getViscoelasticFromSpheresInteraction(tc,en,es) mat=O.materials.append(ViscElMat(density=density,frictionAngle=frictionAngle,**params)) O.bodies.append(utils.facet( [ (-3*Rsph,-0.05,-0.05), (-3*Rsph,-0.05,0.05), (-3*Rsph,0.05,0) ], material=mat)) O.bodies.append([ utils.sphere( (-0.1*Rsph,Rsph,0),Rsph,material=mat,fixed=True), utils.sphere( (-0.0060010254015300003,0.0042204209332099996,0.0009),Rsph,material=mat) ]) O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom(), Ig2_Facet_Sphere_ScGeom()], [Ip2_ViscElMat_ViscElMat_ViscElPhys()], [Law2_ScGeom_ViscElPhys_Basic()], ), GravityEngine(gravity=[0,-9.81,0]), NewtonIntegrator(damping=0), ] O.dt=.01*tc O.saveTmp() from yade import qt qt.View()
_______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
