Thread Previous • Date Previous • Date Next • Thread Next |
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.(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.
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
Thread Previous • Date Previous • Date Next • Thread Next |