← Back to team overview

yade-users team mailing list archive

[Question #248700]: Law2_ScGeom6D_CohFrictPhys_CohesionMoment

 

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

Hi,
Im trying to solve this code,
#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)
O.bodies.append([
   # fixed: particle's position in space will not change (support)
   utils.sphere(center=(0,0,0),radius=.5,fixed=True),
   # this particles is free, subject to dynamics
   utils.sphere((0,0,2),.5)
])

# 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()

O.saveTmp()

To understand better the Law2_ScGeom6D_CohFrictPhys_CohesionMoment, but the error [1] always appears:

[1]ERROR /tmp/buildd/yadedaily-1.05.0-28-g5e71fa5~precise/pkg/common/GravityEngines.cpp:19 action: GravityEngine is deprecated, consider using Newton::gravity instead (unless gravitational energy has to be tracked - not implemented with the newton attribute).
terminate called after throwing an instance of 'std::runtime_error'
  what():  Undefined or ambiguous IPhys dispatch for types FrictMat and FrictMat.

What am I doing wrong

Thanks

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