yade-users team mailing list archive

[Luc OGER <question668315@xxxxxxxxxxxxxxxxxxxxx>]

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

Dear all,

I am doing a simple and slow tilting process of a box filled with monosize spheres. For that, I am increasing by step the angle of the gravity and calculated after that the unbalancedForce.
In real experiments the system will move by successive avalanches depending of the friction coefficient between grains.
In my simulation it is not the case as I guess the number of timestep before doing the calculation of the new unbalancedforce minium is too small : I am using 5000 !
my question is how to define or calculate or adjust this nimimal value depending of the parameter of my expriments: number of grains, friciton coefficient, and so on...

thanks for some clue for my problem.

details of the program:
def checkUnbalanced():
 global i_pas, pas_box, step0,gravity_y,gravity_z,theta_max,converg_min
   # at the very start, unbalanced force can be low as there is only few contacts, but it does not mean the packing is stable
 if O.iter<(step0+5000): return 
 if utils.unbalancedForce()<converg_min:
        i_pas = i_pas+ 1
        gravity_y = 9.81*sin( pas_box*i_pas*3.14/180.0)
        gravity_z = 9.81*cos( pas_box*i_pas*3.14/180.0)
        Newton_integrator.gravity=Vector3(0,-gravity_y,-gravity_z)
        step0 = O.iter
......



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