New question #657560 on Yade: https://answers.launchpad.net/yade/+question/657560
Hello Jan, Bruno and other Yade mates :) Referring to my previous question (#657126), my further task now is to find the force experienced by a hammer (trying spherical first) as it penetrates through a granular media (spherical grains packed in a cylinder here) and plot this force on hammer along with the displacement of the hammer. Currently I have the following script and can get the simulation running but I am trying to proceed further. Please assist. from yade import pack,plot pred=pack.inCylinder((0,0,0),(0,0,0.25),0.10) sp=pack.randomDensePack(pred,radius=0.01,spheresInCell=300) O.bodies.append(sp) yade.qt.Controller() idSteel=O.materials.append(FrictMat(young=210e9,poisson=.25,frictionAngle=.8,label="steel")) hammer=sphere((0,0,0.30),0.04,material=idSteel) hammerID=O.bodies.append(hammer) hammer.state.blockedDOFs='z' calm() hammer.state.vel=Vector3(0,0,-0.5) O.engines=[ForceResetter(),InsertionSortCollider([Bo1_Sphere_Aabb()]),InteractionLoop([Ig2_Sphere_Sphere_ScGeom()],[Ip2_FrictMat_FrictMat_FrictPhys()],[Law2_ScGeom_FrictPhys_CundallStrack()]),NewtonIntegrator(),PyRunner(iterPeriod=100,command="print hammer.state.displ()")] O.dt=0.5*PWaveTimeStep() O.saveTmp() -- 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
