An Potential Error in source code ComputeLinearElasticPFFractureStress.C #27663
Unanswered
DanielChou0916
asked this question in
Q&A Modules: Solid mechanics
Replies: 1 comment
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
-
Dear MOOSE community
I'm currently doing research on phase field fracture model, and when I take a look as the source code "ComputeLinearElasticPFFractureStress.C"
I found that in "computeStrainSpectral" function, there is an potential error.
In the source file:
https://github.com/idaholab/moose/blob/next/modules/solid_mechanics/src/materials/ComputeLinearElasticPFFractureStress.C
Line 106 is given as
F_neg = -lambda * etrneg * etrneg / 2.0 + mu * nval;
However, the negative sign on the first term
-lambda * etrneg * etrneg / 2.0
does not make sense to me.In the official website, the compressive part of elastic energy is given as shwon in following figure:
https://mooseframework.inl.gov/source/materials/ComputeLinearElasticPFFractureStress.html
I also checked several papers to verify this portion, such as
https://www.sciencedirect.com/science/article/pii/S0020768324001744?ref=cra_js_challenge&fr=RR-1 (page 3, equation 11)
And
https://materialstheory.springeropen.com/articles/10.1186/s41313-021-00037-1 (page 5, equation 16).
Could someone verify if this negative sign is indeed an error? If so, correcting it might be important for ensuring the accuracy of the energy calculations in the phase field fracture models.
Appreciate for your time and patience.
Sincerely
Beta Was this translation helpful? Give feedback.
All reactions