Dear all,
Currently I am engaged in solving a geological problem. So we have some soil model in the form of disk, there is a circular hole in the center, which has the initial pressure at the reference time (a modelling of the soil injection with the polymer compound)
How can I recalculate the normal pressure inside the circular hole of the soil model using the DLOAD subroutine if the neighbouring FEs were failed.
This is my question: Question.pdf
Failure criteria is given by Mohr Coulomb criterion (https://en.wikipedia.org/wiki/Mohr%E2%80%93Coulomb_theory).
This is predictable behaviour : predict.png
This is unpredictable behaviour : unpredict.png
This is my subroutine: sand_v9si-02.for_.txt
This is my model: Geo-si-07.inp_.txt
Many thanks for any help
| Attachment | Size |
|---|---|
| predict.png | 23.53 KB |
| unpredict.png | 24.3 KB |
| sand_v9si-02.for_.txt | 8.47 KB |
| Geo-si-07.inp_.txt | 312.04 KB |
| Question.pdf | 636.3 KB |
