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.