yade-users team mailing list archive

[Václav Šmilauer <eudoxos@xxxxxxxx>]

> 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.