Info on traceEnergy and O.trackEnergy
Hi,
I am currently developing a "traceEnergy" feature in the Law2_SCG_
#Question 1:
To calculate the plastic dissipation due to friction for a sliding contact, I see in the Law2_ScGeom_
Real dissip=
if (traceEnergy) plasticDissipation += dissip;
else if(dissip>0) scene->
The same exists in other constitutive laws, e.g. in ElasticContactL
Why do we only check the sign of "dissip" when we are about to add it in the O.energy module (i.e. traceEnergy=False), but not when we are about to add it in the OpenMPAccumulator "plasticDissipa
#Question 2:
In the Law2_SCG_
Also, when I add this dissipation to O.energy, should I check its sign, like we do for "dissip" above, and do you think it would be more convenient to reset the value in every timestep, or monitor the cumulative value (i.e. reset=false), like we do for the dissipation due to friction?
#Question 3:
In the Law2_SCG_
phys->normalForce -= phys->normalVis
When I calculate the elastic potential energy, should I consider the reduced potential energy or the originally elastic one? i.e.
Reduced potential energy:
scene-
Elastic potential energy :
scene-
Apologies in advance for not opening three different questions; I just feel these should be in one thread, as they concern a very specific part of the source code.
Many thanks,
Vasileios
[1] https:/
[2] https:/
[3] https:/
Question information
- Language:
- English Edit question
- Status:
- Expired
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply: