ABSTRACT

At each step, K" and df " are formed, and (39) is solved using a frontal method. Total displacements and pore pressures are then accumulated via.

u J + du" (40a) J 147, = up-1 + du", (40b)

Finally, stresses and strains are determined, and the geometry is updated for the next time step. As mentioned above, any equilibrium or yield surface violations are imposed as corrections to the subsequent step.