plastic energy dissipation
| Question information | |
|---|---|
| Author | zhanglingran130 |
| Date | 16 Aug 2013 |
Zhanglingran130 is having issues with plastic energy dissipation calculation in Yade's source code for Ubuntu. Jan suggested modifications to the source code including adding bendingElastEnergy and twistElastEnergy to Law2_ScGeom6D_CohFrictPhys_CohesionMoment. The modified source codes were recompiled and tested in simple 2D tests, but when combined with Law2_ScGeom_FrictPhys_CundallStrack, bendingElastEnergy was found to be 'nan' or 'inf'. Bruno suggested filtering the interactions and adding elasticEnergy() to Law2_ScGeom_FrictPhys_CundallStrack. Zhanglingran130 is advised to filter the interactions using an efficient method that avoids dynamic_casts. Vaclav suggested adding an 'int index' to each Indexable and initializing it once. Zhanglingran130 wants to obtain local plastic energy dissipation values by filtering some interactions in Yade's source code, but has not succeeded yet due to the complexity of the issue. He hopes for a solution that reduces computation time.
Link to full discussion: https://answers.launchpad.net/yade/+question/234129