New question #404127 on Yade:
https://answers.launchpad.net/yade/+question/404127
Hi,
I simulated a steel ball impacted on a particle assembly at a high speed, but
my code got an error, that is the steel ball not move. My codes is as follows,
please help me to solve this problem. Many thanks!
from yade import pack
##### define materials for particles and mill #####
rock=O.materials.append(CohFrictMat(young=5.98e7,poisson=0.3,density=2678,frictionAngle=0.5,normalCohesion=1.2e10,shearCohesion=1.2e9,momentRotationLaw=True))
steel=O.materials.append(FrictMat(young=3.06e11,poisson=0.29,density=7861,frictionAngle=0.545))
##### create partciles assembly #####
pred=pack.inAlignedBox((-0.1,-0.1,-0.1),(0.1,0.1,0.1))
spheres=pack.randomDensePack(pred,radius=0.005,rRelFuzz=0.3,spheresInCell=1500,material=rock,color=(1,0,1))
O.bodies.append(spheres)
##### create a rigid ball #####
ball=O.bodies.append(sphere(center=(0,0,0.5),radius=0.1,color=(0,1,1),material=steel))
ball.blockedDOFs='xyXYZ'
ball.state.vel=Vector3(0,0,-1) ##### set velocity along axis-Z
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom6D(interactionDetectionFactor=1.1)],
[Ip2_FrictMat_FrictMat_FrictPhys(),Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True)],
[Law2_ScGeom_FrictPhys_CundallStrack(),Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
),
NewtonIntegrator(damping=0.3,gravity=(0,0,-9.81))
]
from yade import qt
qt.View()
qt.Controller()
#qt.Renderer().shape=0
#qt.Renderer().intrPhys=1
O.dt=0.001*PWaveTimeStep()
--
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
