Total elastic shear energy in Hertz-Mindlin model

Hi Jiangtao,
I used the Hertz-Mindlin model for some time, but I never looked at the energies.
As I read the source code, in shearEnergy not only the current increment is saved, but the sum of all increments (see below). What exactly would you like to compute?
Regards,
Bettina

Yade-1.20.0, HertzMindlin.cpp, Lines 423 - 430
/************************/
/* SHEAR ELASTIC ENERGY */
/************************/
// NOTE: shear elastic energy calculation must come after the MC criterion, otherwise displacements and forces are not updated
if (calcEnergy) {
 shearEnergy += (us_elastic-prevUs_el).dot((shearElastic+prev_FsElastic)/2.); // NOTE: no additional energy if we perform sliding since us_elastic and prevUs_el will hold the same value (in fact us_elastic is only keeping the elastic part). We work out the area of the trapezium.
}

Dear Bettina,
Thanks for your kindly reply. My question before is not so clear and let me restate my problem. If we want to know the energy flows during a specific stage. However, calcEnergy is false before the stage and the shearEnergy will start from zero. But we know that the elastic shear energy is not zero.

Best regards,
Jiangtao

Hi Jiangtao,

sorry, but it is still not clear to me, what exactly is your question. Would you like to calculate the change in shear energy between two moments in time? You write that you set calcEnergy=False in the beginning, changing it to True at some later time. Wouldn’t it be easier to set it True from the beginning?
Maybe you could write a minimal example with only two spheres, or so, and post the code. Then, you can describe, what is calculated in shearEnergy and what you would like to calculate.

Best regards,
Bettina

This question was expired because it remained in the 'Open' state without activity for the last 15 days.