Computation of Contact energy with .getEnergy("contact")
Hello,
Using the python interface, I am studying the fragmentation of a brittle bar using cohesive elements with the cohesive linear law that implements contact using a penalty coefficient ("penalty" in the material input file) to model contact between fragments. During my simulation, I compute the contact energy by using Ec = model.getEnergy("contact") as found in the Akantu documentation. Using a non zero penalty coefficient should consider contact and thus result in a non zero contact energy but whatever value I use for "penalty", the contact energy is always zero resulting in the total energy not being conserved and I do not understand why ?
Thank you in advance for your help