yade-users team mailing list archive

[yakir <question255199@xxxxxxxxxxxxxxxxxxxxx>]

New question #255199 on Yade:
https://answers.launchpad.net/yade/+question/255199

Hello all,

I'm trying to create a pack of high adhesion balls. so i created a simple script to test the CohFrictMat engine but somehow the blue ball refuse to stick.
what do i need to do?



#simple script

from yade import ymport, pack, plot 
import operator 
rMean=0.0025

O.materials.append(FrictMat(young=210e10,poisson=.7,frictionAngle=.1,label="steel")) # material for the box and the wedge

O.materials.append(CohFrictMat(
young=50000,poisson=1,density=1074,frictionAngle=radians(0),
isCohesive=True,normalCohesion=10000000,shearCohesion=10000000,
etaRoll=0.1,momentRotationLaw=True,alphaKr=0.1,alphaKtw=0.1,label="sphereMat"))

O.bodies.append(utils.sphere(center=(2e-2,4e-2+(2*rMean),0e-2),radius=rMean,material="sphereMat",fixed=False,color=(1,0,0)))
O.bodies.append(utils.sphere(center=(2e-2-(2*rMean),4e-2,0e-2),radius=rMean,material="sphereMat",fixed=False,color=(0,0,1)))
O.bodies.append(utils.sphere(center=(2e-2,4e-2,0e-2),radius=rMean,material="sphereMat",fixed=True,color=(0,1,0)))

O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
   InteractionLoop(
   [Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Facet_Sphere_ScGeom()],
   [Ip2_FrictMat_FrictMat_FrictPhys(),
   Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,label="cohesiveIp")],
   [Law2_ScGeom_FrictPhys_CundallStrack(), 
   Law2_ScGeom6D_CohFrictPhys_CohesionMoment(useIncrementalForm=True,label='cohesiveLaw')]
   ),
   NewtonIntegrator(damping=0.8,gravity=[0,0,-9.81]), 
]
O.dt=.0005*utils.PWaveTimeStep() 

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.