Re: [Yade-users] [Question #293342]: Problem with simulation
Question #293342 on Yade changed: https://answers.launchpad.net/yade/+question/293342 Status: Open => Expired Launchpad Janitor expired the question: This question was expired because it remained in the 'Open' state without activity for the last 15 days. -- 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
Re: [Yade-users] [Question #293342]: Problem with simulation
Question #293342 on Yade changed: https://answers.launchpad.net/yade/+question/293342 Status: Answered => Open Alex is still having a problem: Hello Jérôme and Klauss, Thank you for the answers. They have helped a lot. I'm sorry for the delay in replying and for the title... The next time I will be more specific and more original. ;-) I rewrote my code and I continue with the same problem, that is, when I run the simulation the cylinder goes inside the spheres and that's not what I want. Once again, I ask if someone could help in this problem. Below is my new code. import sys from yade import ymport from yade import qt from numpy import linspace material definition def sphereMat(): return CohFrictMat(young=8e5,poisson=0.3,density=4e3, frictionAngle=radians(30),normalCohesion=4e4,shearCohesion=1e5,momentRotationLaw=False,label='sphereMat') def gridNodeMat(): return CohFrictMat(young=8e5,poisson=0.3,density=4e3,frictionAngle=radians(30), normalCohesion=1e5,shearCohesion=1e5,momentRotationLaw=True,label='gridNodeMat') def gridCoMat(): return CohFrictMat(young=8e5,poisson=0.3,density=4e3, frictionAngle=radians(30), normalCohesion=1e5,shearCohesion=1e5,momentRotationLaw=True,label='gridCoMat') end material definition ### O.engines=[ ForceResetter(), InsertionSortCollider([ Bo1_GridConnection_Aabb(), Bo1_Sphere_Aabb() ]), InteractionLoop( Geometric interactions [ Ig2_GridNode_GridNode_GridNodeGeom6D(), Ig2_Sphere_GridConnection_ScGridCoGeom() ], [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,setCohesionOnNewContacts=True), ], [Law2_ScGeom6D_CohFrictPhys_CohesionMoment(label='law'), Law2_ScGridCoGeom_CohFrictPhys_CundallStrack() ] ), Motion equation NewtonIntegrator(damping=0.4,gravity=[0,0,0]), PyRunner(command='main()',iterPeriod=10) ] O.dt=1e-7 rCyl=0.025 nL= 2 L=0.3 Create all nodes first : nodesIds=[] for i in linspace(0,L,nL): nodesIds.append( O.bodies.append( gridNode([0.3,i + 0.4,2.02],rCyl,wire=False,fixed=False,material=gridNodeMat) ) ) Now create connection between the nodes for i,j in zip( nodesIds[:-1], nodesIds[1:]): O.bodies.append( gridConnection(i,j,rCyl, material=gridCoMat ) ) read a sphere's package PACKING = '10K_P1' O.bodies.append(ymport.text(PACKING+'.spheres',scale=1.,shift=Vector3(0,0,0),material=sphereMat)) qt.View(); def main(): global Fn,Ft O.bodies[0].dynamic=False O.bodies[0].state.vel[2]=-500.0 O.bodies[1].state.vel[2]=-500.0 O.bodies[2].state.vel[2]=-500.0 -- 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
Re: [Yade-users] [Question #293342]: Problem with simulation
Question #293342 on Yade changed: https://answers.launchpad.net/yade/+question/293342 Status: Open => Answered Jérôme Duriez proposed the following answer: Hi, According to the name of Ig2_Sphere_ChainedCylinder_CylScGeom6D(), the interaction between Sphere and chainedCylinder has a "geom" of CylScGeom6D type Thus, you need to include in your simulations a "Law2" that can handle such interactions, wich is not the case here (none of your Law2 functors include CylScGeom6D in its name). I'd advice you to give a close look to https://yade- dem.org/doc/user.html#functors-choice, instead of using brute force and stuffing as many functors as possible (without the right one, by lack of luck) into your simulation ;-) ;-) PS: the title of your question is not so original -- 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
[Yade-users] [Question #293342]: Problem with simulation
New question #293342 on Yade: https://answers.launchpad.net/yade/+question/293342 Hello friends, I created the simulation below and I had a problem. My objective is to make the cylinder collide with spheres. But, I can't do this. When I run the simulation the cylinder goes inside the spheres. I ask if someone could help in this problem. import sys from yade import ymport, utils , plot import math Material microproperties -> SHALE! intR=1.22955 DENS=2659.6453 YOUNG=40e9 FRICT=25 ALPHA=0.25 TENS=10e5 COH=10e6 aRes=1e-3 damp=0.2 jointFRICT=25 jointDIL=15 material definition def sphereMat(): return JCFpmMat(type=1,density=DENS,young=YOUNG,poisson = ALPHA,frictionAngle=radians(FRICT), tensileStrength=TENS,cohesion=COH,jointNormalStiffness=YOUNG/10.,jointShearStiffness=0.2*YOUNG/10., jointTensileStrength=0.,jointCohesion=0.,jointFrictionAngle=radians(jointFRICT),jointDilationAngle=radians(jointDIL)) O.engines=[ ForceResetter(), InsertionSortCollider([ Bo1_ChainedCylinder_Aabb(), Bo1_Sphere_Aabb() ]), InteractionLoop( # Geomeetric interactions [Ig2_ChainedCylinder_ChainedCylinder_ScGeom6D(), Ig2_Sphere_ChainedCylinder_CylScGeom6D(), Ig2_Sphere_Sphere_ScGeom6D(), Ig2_Sphere_Sphere_L3Geom(),Ig2_Facet_Sphere_L3Geom() ], [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,setCohesionOnNewContacts=True), Ip2_FrictMat_FrictMat_FrictPhys() ], [Law2_ScGeom6D_CohFrictPhys_CohesionMoment(label='law'), Law2_ScGeom_FrictPhys_CundallStrack(), Law2_L3Geom_FrictPhys_ElPerfPl(), ] ), ## Motion equation NewtonIntegrator(damping=0.4,gravity=[0,0,0]), PyRunner(command='main()',iterPeriod=10) ] O.dt=1e-7 # create a cylinder cil = chainedCylinder(begin=Vector3(0.3,0.4,2.02), end=Vector3(0.3,0.6,2.02),radius=0.025,color=Vector3(0.6,0.5,0.5)) # read a sphere's package PACKING = '10K_P1' O.bodies.append(ymport.text(PACKING+'.spheres',scale=1.,shift=Vector3(0,0,0),material=sphereMat)) yade.qt.View(); def main(): global Fn,Ft O.bodies[0].dynamic=False O.bodies[0].state.vel[2]=-500.0 -- 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
