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
