Question #248700 on Yade changed:
https://answers.launchpad.net/yade/+question/248700
Status: Open => Answered
Luc Sibille proposed the following answer:
Hi Luis,
There is a problem with the definition of the material properties (I
imagine the default material is not a CohFricMat material but rather a
FricMat?), by defining a CohFricMat material as introduced below in your
script it works:
#exercicio de Law
# basic simulation showing sphere falling ball gravity,
# bouncing against another sphere representing the support
# DATA COMPONENTS
# add 2 particles to the simulation
# they the default material (utils.defaultMat)
matSph=CohFrictMat(
young=170e6,
poisson=0.2,
frictionAngle=radians(17),
isCohesive=False,
alphaKr=1.8,
alphaKtw=0,
etaRoll=1.8,
momentRotationLaw=True,
density=2600,
label='spheres',
)
O.materials.append(matSph)
O.bodies.append([
# fixed: particle's position in space will not change (support)
utils.sphere(center=(0,0,0),radius=.5,fixed=True,material='spheres'),
# this particles is free, subject to dynamics
utils.sphere((0,0,2),.5,material='spheres')
])
# FUNCTIONAL COMPONENTS
# simulation loop -- see presentation for the explanation
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom6D()],
[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(label="cohesiveIp")],
[Law2_ScGeom6D_CohFrictPhys_CohesionMoment(label="cohesiveLaw")]
),
# apply gravity force to particles
GravityEngine(gravity=(0,0,-9.81)),
# damping: numerical dissipation of energy
NewtonIntegrator(damping=0.1)
]
O.dt=.5e-3*utils.PWaveTimeStep()
Best,
Luc
--
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
