Question #702751 on Yade changed: https://answers.launchpad.net/yade/+question/702751
Status: Needs information => Open Roxana Saghafian Larijani gave more information on the question: Please find the MWE below: from yade import ymport from yade import pack from yade import utils, plot,wrapper fr = 0.38 rho = 2500 En = 0.7 Et = 0.7 poi=0.21 yoM=63e6 O.dt = 0.000005 r = 0.01 #0.002381 ##definig material mat=O.materials.append( ViscElCapMat(frictionAngle=fr, density=rho, poisson=poi, young= yoM,en=En, et=Et) ) #defining the spheres sp=pack.SpherePack() sp.makeCloud((0.6,0.6,0.05),(1,1,0.15),rMean=r) sp.toSimulation(material=mat) Nprtcl=len(O.bodies) Drum=geom.facetCylinder((0.8,0.8,0.1),radius=0.5,height=0.2,orientation=Quaternion(Vector3(0,0,1),(pi/2.0)), wallMask=4) walls = O.bodies.append(Drum) ##periodic boundary O.periodic = True O.cell.setBox(5,5,0.2) #### ##engine O.engines = [ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb(),Bo1_Wall_Aabb()],allowBiggerThanPeriod = True), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom()], [Ip2_ViscElCapMat_ViscElCapMat_ViscElCapPhys()], [Law2_ScGeom_ViscElCapPhys_Basic()], ), NewtonIntegrator(gravity=[0, -9.8, 0]), RotationEngine(ids=walls,rotationAxis=[0,0,1],rotateAroundZero=True,zeroPoint=[0.8,0.8,0.1],angularVelocity=1.68) ] from yade import qt qt.View() O.run() >What is the result with Law2_ScGeom_FrictPhys_CundallStrack? I checked it and the problem does not get solved by changing to Cundall. >But firstly I would guess the problem is due to the fact "If I remember correctly, the large bodies are not "periodically copied", so the spheres from "next" cells do not interact with the wall." discussed in [2], #4 I think this is the problem! Because it is like the particles pass through the drum wall and does not interact with it! Regards, Roxana Saghafian ________________________________ From: nore...@launchpad.net <nore...@launchpad.net> on behalf of Jan Stránský <question702...@answers.launchpad.net> Sent: Thursday, August 18, 2022 9:20 AM To: Saghafian Larijani, Roxana (UT-ET) <r.saghafianlarij...@utwente.nl> Subject: Re: [Question #702751]: periodic boundary in the direction of the length of a rotating drum (cylinder ) Your question #702751 on Yade changed: https://answers.launchpad.net/yade/+question/702751 Status: Open => Needs information Jan Stránský requested more information: > However, I still have a problem. please, provide a MWE as discussed above. > Might it be because I am using (Law2_ScGeom_ViscElCapPhys_Basic()) as it is mentioned in the question 429604? What is the result with Law2_ScGeom_FrictPhys_CundallStrack? But firstly I would guess the problem is due to the fact "If I remember correctly, the large bodies are not "periodically copied", so the spheres from "next" cells do not interact with the wall." discussed in [2], #4 If it is this case or not would be clear with the MWE. Cheers Jan [2] https://answers.launchpad.net/yade/+question/701266 -- To answer this request for more information, you can either reply to this email or enter your reply at the following page: https://answers.launchpad.net/yade/+question/702751 You received this question notification because you asked the question. -- 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@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp