[GeoMechanicsApplication] Add documentation for compression cap for material models #14062
Conversation
WPK4FEM
left a comment
WPK4FEM
left a comment
There was a problem hiding this comment.
Hi Mohamed,
This is well on its way and in some places already more than asked.
Please find my remarks.
Regards, Wijtze Pieter
applications/GeoMechanicsApplication/custom_constitutive/README_cap.md
Outdated
Show resolved
Hide resolved
applications/GeoMechanicsApplication/custom_constitutive/README_cap.md
Outdated
Show resolved
Hide resolved
applications/GeoMechanicsApplication/custom_constitutive/README_cap.md
Outdated
Show resolved
Hide resolved
applications/GeoMechanicsApplication/custom_constitutive/README_cap.md
Outdated
Show resolved
Hide resolved
applications/GeoMechanicsApplication/custom_constitutive/README_cap.md
Outdated
Show resolved
Hide resolved
applications/GeoMechanicsApplication/custom_constitutive/README_cap.md
Outdated
Show resolved
Hide resolved
applications/GeoMechanicsApplication/custom_constitutive/README_cap.md
Outdated
Show resolved
Hide resolved
|
|
||
| ```math | ||
| p_c = p_{c0} + H \epsilon^p | ||
| ``` | ||
| where: | ||
| - $`p_{c0}`$ = the initial cap position | ||
| - $`H`$ = the hardening modulus | ||
| - $`\epsilon^p`$ = the plastic volumetric strain | ||
|
|
||
| This produces compaction hardening: as the material densifies, it supports higher compressive stresses. |
There was a problem hiding this comment.
Please leave out hardening until the moment that we really address it.
WPK4FEM
left a comment
WPK4FEM
left a comment
There was a problem hiding this comment.
Hi Mohamed,
Thank you for this documentation, I think its very close and I hope that taking up the remarks is now a small effort. Please proceed with the implementation of the cap yield surface and the derivatives of the flow function like was done when MC was started ( just F, dG/d sigma and unit tests for that as stand alone thing )
Regards, Wijtze Pieter
|
|
||
| ### Mohr–Coulomb yield surface | ||
|
|
||
| In the $`(\sigma, \tau)`$ stress space, the Mohr-Coulmb yiel surface is expressed as: |
There was a problem hiding this comment.
| In the $`(\sigma, \tau)`$ stress space, the Mohr-Coulmb yiel surface is expressed as: | |
| In the $`(\sigma, \tau)`$ stress space, the Mohr-Coulomb yield surface is expressed as: |
applications/GeoMechanicsApplication/custom_constitutive/README_cap.md
Outdated
Show resolved
Hide resolved
applications/GeoMechanicsApplication/custom_constitutive/README_cap.md
Outdated
Show resolved
Hide resolved
applications/GeoMechanicsApplication/custom_constitutive/README_cap.md
Outdated
Show resolved
Hide resolved
applications/GeoMechanicsApplication/custom_constitutive/README_cap.md
Outdated
Show resolved
Hide resolved
applications/GeoMechanicsApplication/custom_constitutive/README_cap.md
Outdated
Show resolved
Hide resolved
| Here, we need to convert the compression cap yield surface from $(p, q)$ coordinates to $(\sigma, \tau)$ coordinates. The conversion is to be followed ... | ||
|
|
||
| ### Compression cap in $(\siga, \tau)$ coordinates | ||
| ... | ||
|
|
There was a problem hiding this comment.
This can be left out for the time being, until we solved that challenge.
applications/GeoMechanicsApplication/custom_constitutive/README_cap.md
Outdated
Show resolved
Hide resolved
WPK4FEM
left a comment
WPK4FEM
left a comment
There was a problem hiding this comment.
Hi Mohamed,
Thank you for the modifications, 1 more spelling mistake found and 2 small remarks. With that it should be good to go.
Wijtze Pieter
📝 Description
In order to start with implementation of the compression cap, it is more efficient to document the mathematical formulations, algorithms and the way how to implement. This issue focuses on documentation of the compression cap.
At the moment, the formulation is based on p and q. Later we need to find a way to find the intersection between the yield surfaces.
🆕 Changelog