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.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
