DLOAD, element deletion, new active load surface

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. 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

This is my question: Question.pdf

Many thanks for any help