yade-users team mailing list archive

[Rodolfo França de Lima <question647311@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi! I am new user in yade and I would like to know if is possible build a chart from sphere velocity of  script below.  Thanks for helping.

from yade import pack, plot, qt

isopor=CohFrictMat(density=23)
idisopor=O.materials.append(isopor)

O.bodies.append([
   sphere(center=(0,2,100),radius=.5,material=idisopor,color=(1,1,0)),
])

h=0.0
#tampa
p1=(-15,15,h)
p2=(-15,-15,h)
p3=(15,-15,h)
p4=(15,15,h)

f1=utils.facet([p1,p2,p3],wire=False)
f2=utils.facet([p1,p3,p4],wire=False)

O.bodies.append(f1)
O.bodies.append(f2)
O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
   InteractionLoop(
      # handle sphere+sphere and facet+sphere collisions
      [Ig2_Sphere_Sphere_L3Geom(),Ig2_Facet_Sphere_L3Geom()],
      [Ip2_FrictMat_FrictMat_FrictPhys()],
      [Law2_L3Geom_FrictPhys_ElPerfPl()]
   ),
   DragEngine(ids=[0],Cd=0.45, Rho=1.225),
   NewtonIntegrator(damping=0.05,gravity=(0,0,-9.81)),
   
]
O.dt=.05*utils.PWaveTimeStep()

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