yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #03227
Re: Granular ratchetting explained
Thanks. Everything is as expected then, always good to hear.
For motion engines, I prefer not to use them as we can simply assign
velocity directly. Soon, you will not have to modify anything in Newton,
since this behaviour of non-dynamic bodies is currently being implemented.
Test scripts without Newton is not a good thing IMHO, you are putting
yourself in a special case wich will never happen in routine
simulations. It can hide existing problems or create problems that would
not exist with Newton.
Bruno
chiara modenese a écrit :
If you still want to compare ratcheting in
Law2_ScGeom_FrictPhys_Basic, a quick hack is to change
"true"->"false" at line ElasticContactLaw.cpp:77
Please tell me the result. :)
Hi Bruno,
I revised a little bit the script (I attached it again) I sent you
concerning RotationEngine (which applies the rolling) and now it works
well. I preferred to maintain the use of Rotation and StepDisplacer
engines, I did not have time to modify NewtonIntegrator and anyway
both solutions should work.
I tested it using Law2_ScGeom_FrictPhys_Basic. Avoiding granular
ratcheting (which is the default-fixed option in the code atm) no
shear force is present at the end of the loop (just run the script as
it is). Instead, I tried the case in which we want to consider
granular ratcheting (changing directly the code as you suggest above)
and a net shear force remains at the end (in the case I was running, I
have got a value of 0.183), although the ball comes back exactly at
the same initial position. Meanwhile, I have also printed out the
values of the branches applied to the rotations in ScGeom, and yes
they differ, in the case with ratcheting, between the first and second
rotation. They are exactly the same avoiding ratcheting. Maybe I try
also with other contact laws and I let you know.
Chiara
------------------------------------------------------------------------
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@xxxxxxxxxxxxxxxxxxx
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
--
_______________
Bruno Chareyre
Associate Professor
Grenoble INP
Lab. 3SR
BP 53 - 38041, Grenoble cedex 9 - France
Tél : 33 4 56 52 86 21
Fax : 33 4 76 82 70 43
________________
Follow ups
References