yade-users team mailing list archive

[sebbah <question677994@xxxxxxxxxxxxxxxxxxxxx>]

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

hi,
i have wrote a program of bouncing ball on the soil and i want to save the height of sphere in each time step in a file 'height.txt'
the problem is that the script does not  calculate the heights and the  time step .
my code is:

#!/usr/bin/python
# coding: utf-8

from numpy import*
from yade import plot
 
""" add the particule to simulation """
O.bodies.append([utils.sphere((0,0,2),.5)])

""" add a wall to simulation """
O.bodies.append([utils.wall(-1,axis=2)])


O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Wall_Aabb()]),
   InteractionLoop(
       [Ig2_Sphere_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom()], # collision geometry
       [Ip2_FrictMat_FrictMat_FrictPhys()], # collision "physics"
       [Law2_ScGeom_FrictPhys_CundallStrack()] # contact law -- apply forces
   ),
   # Apply gravity force to particles. damping: numerical dissipation of energy.
   NewtonIntegrator(gravity=(0,0,-9.81),damping=0.1)
 ]

#fixing step time
O.dt=.5e-4*utils.PWaveTimeStep()

# save the simulation
O.saveTmp()

""" open file"""
fil=open('height.txt', 'w') 
fil.write('%e %e\n' % (t/T, height))

""" print the sphere's height """
def height_sphr():
    h = [O.bodies.state.pos[2]]
    fil.write('%d, %.4f\n' % (O.iter,h))

fil.close()


Any ideas why it doesn't work?

Thanks in advance,

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