New question #212480 on Yade: https://answers.launchpad.net/yade/+question/212480
Hi, I have encountered something strange when conducting gravity deposition: some particles will go outside the box. Can you give some suggestions for solving this problem? Thanks. Nait My script is as following: ####SCRIPT#### from yade import pack rm = 0.8 rr = 0.6 compFricDegree = 30 finalFricDegree = 30 stabilityThreshold=0.01 young=5e6 mn,mx=Vector3(0,0,0),Vector3(20,80,200) mxp = Vector3(20,80,150) mnp = Vector3(20,80,90) thick = 0.01 targetporosity = 0.5 O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=radians(compFricDegree),density=2600,label='spheres')) O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=0,density=0,label='walls')) walls=utils.aabbWalls([mn,mx],thickness=thick,material='walls') wallIds=O.bodies.append(walls) sp=pack.SpherePack() sp.makeCloud(mn,mxp,rMean=rm,rRelFuzz=rr,porosity=targetporosity) O.bodies.append([utils.sphere(center,rad,material='spheres') for center,rad in sp]) O.dt=.5*utils.PWaveTimeStep() # initial timestep, to not explode right away O.usesTimeStepper=True O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()]), InteractionLoop( # handle sphere+sphere and facet+sphere collisions [Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()], [Ip2_FrictMat_FrictMat_FrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack()] ), GravityEngine(gravity=(0,0,-9.81)), NewtonIntegrator(damping=0.6), ] O.run(2000, True) while 1: O.run(500, True) unb=unbalancedForce() print 'unbalanced force:',unb if unb<0.01: break yade.qt.Controller(), yade.qt.View() ####END OF SCRIPT#### -- You received this question notification because you are a member of yade-users, which is an answer contact for Yade. _______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp