> Put a "Real plasticWork" in the functor. Compute the energy dissipated > at one contact on time increment dt, and include a "plasticWork+= ..." > in the if(plasticityCondition) bracket of the functor. > Don't ask me how to define plastic work at contact with an > elasto-plastic Hertz-based law... ;)
I think since the formulation is incremental, you can use |Δε|*|σ| for energy dissipated in plastic slip by Δε at plastic stress σ...? Generally, there is no unified function for dissipation. You could, though, sum kinetic energy of particles, potential energy (if there is some potential field), and subtract cummulative external work (boundary conditions). That should give you pretty good image of system energy evolution, including damping, plasticity and dissipation by numerical entropy ;-) Cheers, v. _______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
