Implementation of a elastic-plastic rolling resistance model

Hi everyone,

I want to implement the following rolling resistance torque into simulation.


J. Ai, J.F. Chen, J.M. Rotter, J.Y. Ooi. Assessment of rolling resistance models in discrete element simulations. Powder Technology, 2011, 206: 269-282. Redirecting

I have read and edited the file called cala_force_dem_mod.f. The expressions such as contact forces and torques in this file may not directly related to time and time step. However, the equations here involve T(t), T(t+∆t) and ∆t. Which files could I go to and edit? Are T(t) and T(t+∆t) definated in the same file?

And, is OMEGA_NEW(LL, : ) in cala_force_dem_mod.f means the angular velocity of particle LL? In which file is OMEGA_NEW(LL, : ) defined or calculated?

Thank you very much.