DLOAD, element deletion, new active load surface

[atomicra]

Dear all,

I am actually engaged in the decision of geological problems. So we have a some model with the soil disk, there is a circular hole in the centre, witch has initial pressure at the first time(a modelling of the soil injection with the polymer solution). How can I recalculate by the subroutine DLOAD a normal pressure inside of the circular hole in the soil model if a nearest FEs were failed. Failure criteria is given by Mohr Coulomb criterion (

https://en.wikipedia.org/wiki/Mohr–Coulomb_theory).

This is predictable behaviour :

http://files.engineering.com/getfil...-aa63-4b61-bb80-5b65df5c7f7c&file=predict.png

This is unpredictable behaviour :

http://files.engineering.com/getfil...b12-442a-bb32-d718de27dcfc&file=unpredict.png

This is my question:

http://files.engineering.com/getfil...-4d96-87c5-2bee8d7a70c5&file=Soil_Failure.pdf

This is my subroutine:

http://files.engineering.com/getfil...-467e-9869-0c910bbfa140&file=sand_v9si-02.for

This is my model:

http://files.engineering.com/getfil...164-4ad8-832f-52dba4feb314&file=Geo-si-07.inp

Many thanks for any help