(If an interaction is broken, then it is a different issue, i.e. that
dissipated energy of broken interactions should be saved somehow, since
otherwise we lose the number).
the calculation of plasticDissipation is incremental. If the contact
breaks, then nothing is lost, you only stop incrementing it.
Good point. Something is missed each time a contact is lost, since this
contact dissipated something on [t,t+dt] even if it is lost at time t+dt.
This can be negligible in dense packing (triaxial) but significant in
bouncing spheres.
It doesn't explain the excessive growth of plasticE though...
I'm wondering if the problem is not just the time integration scheme :
if a bouncing is described in 2-3 steps, the approximation of
derivatives is horrible, and it might well create energy artificialy.
Perhaps Cundall introduced global damping for a reason...
Bruno
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
