yade-users team mailing list archive

[Wenbin Fei <question269513@xxxxxxxxxxxxxxxxxxxxx>]

Question #269513 on Yade changed:
https://answers.launchpad.net/yade/+question/269513

Wenbin Fei posted a new comment:
Hi Jose,

I am so interested in your work.

However, I tried to run your scripts, it still cannot work normally.

Ther scripts I ran is as follows:

##################################
from yade import polyhedra_utils
import random

polyMat = PolyhedraMat(density=3000,young=1e10,poisson=.5,frictionAngle=atan(0.8))
frictMat = FrictMat(density=8000,young=1e9,poisson=.5,frictionAngle=atan(0.8))

O.materials.append((polyMat,frictMat))

t=polyhedra_utils.polyhedra(polyMat,(2,2,2),v=[])
t.state.pos = (0,0,5)
O.bodies.append(t)

O.bodies.append(geom.facetBox((0,0,-2),(5,5,10)))

O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Polyhedra_Aabb(),Bo1_Facet_Aabb()],verletDist=0.5),   
   InteractionLoop(
      [Ig2_Facet_Polyhedra_PolyhedraGeom()], # [1]
      [Ip2_PolyhedraMat_PolyhedraMat_PolyhedraPhys()], # [2]
      [Law2_PolyhedraGeom_PolyhedraPhys_Volumetric()], # [3]
),
    NewtonIntegrator(damping=0.5,gravity=(0,0,-9.81))
]

O.dt=0.025*polyhedra_utils.PWaveTimeStep()


from yade import qt
v=qt.View()
v.axes=True

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.