[GeoMechanicsApplication] Capacity checks based on Mohr Coulomb in stress strain utilities #12258
Conversation
…ts, after messy merge on old branch.
WPK4FEM
changed the title
applications/GeoMechanicsApplication/custom_utilities/README.md
Outdated
Show resolved
Hide resolved
There was a problem hiding this comment.
I'm very happy with this set of tools that helps to produce output to interpret results of the Mohr-Coulomb model. Thanks for all the hard work that you have put into this, I really appreciate that. I have just a few minor remarks which I hope you will find useful.
|
|
||
| ### Green Lagrange strain tensor | ||
|
|
||
| With current configuration $x$ and reference configuration $X$, the deformation gradient $F = \frac{x}{X}$ and unit tensor I: |
There was a problem hiding this comment.
Nitpicking:
| With current configuration $x$ and reference configuration $X$, the deformation gradient $F = \frac{x}{X}$ and unit tensor I: | |
| With current configuration $x$ and reference configuration $X$, the deformation gradient $F = \frac{x}{X}$ and unit tensor $I$: |
|
|
||
| which brings: | ||
|
|
||
| $$\bar{\theta}_s = \frac{1}{3} asin( - \frac{27}{2} \frac{(\sigma_1 - p)(\sigma_2 - p)(\sigma_3 - p)}{q^3})$$ |
There was a problem hiding this comment.
Shouldn't we use $\arcsin$ rather than $asin$, to change
| double StressStrainUtilities::CalculateDenominator(const Vector& rStressVector, double Phi) | ||
| { | ||
| const double lode_angle = CalculateLodeAngle(rStressVector); | ||
| return sqrt(3.) * std::cos(lode_angle) - std::sin(lode_angle) * std::sin(Phi); |
There was a problem hiding this comment.
One minor thing:
| return sqrt(3.) * std::cos(lode_angle) - std::sin(lode_angle) * std::sin(Phi); | |
| return std::sqrt(3.) * std::cos(lode_angle) - std::sin(lode_angle) * std::sin(Phi); |
| KRATOS_EXPECT_DOUBLE_EQ(0.0, StressStrainUtilities::CalculateMohrCoulombShearCapacity( | ||
| stress_vector, cohesion, friction_angle)); |
There was a problem hiding this comment.
Is this the formatting produced by clang-format?
There was a problem hiding this comment.
Reran the formatting.
| stress_vector, cohesion, friction_angle)); | ||
| } | ||
|
|
||
| KRATOS_TEST_CASE_IN_SUITE(CheckCalculateMCShearCapacityShearOnly, KratosGeoMechanicsFastSuite) |
There was a problem hiding this comment.
Perhaps it is more readable and consistent to use MohrCoulomb rather than MC? If yes, then the same comment applies to the next two test cases.
📝 Description
Added Lode angle and Mohr Coulomb based capacity checks to geo stress/strain utilities
🆕 Changelog