Following the standard procedure of the finite element method, the domain W is first divided into subdomains (elements), then the displacement are pore water pressure fields within an element are expressed in terms of a finite number of nodal values and interpolation (shape) functions as equations (9.a, 9.b).
In what follows, the applications of the weighted residual method to the partially saturated dynamic u-p formulation of (1) is demonstrated.