Extended precise integration method for consolidation of transversely isotropic poroelastic layered media

A new numerical method for calculating the consolidation behavior of the stratified, transversely isotropic and poroelastic material is presented by combining the extended precise integration algorithm with the integral transformation techniques. Starting with the governing partial differential equations of a saturated medium with transversely isotropic skeleton and compressible fluid constituents, an ordinary differential matrix equation is deduced with the aid of a Laplace–Hankel transform. An extended precise integration method for internal loading situations is proposed to solve the ordinary differential matrix equation in the transformed domain, and the actual solution is recovered by a numerical inverse transformation. Numerical examples are also provided to prove the feasibility of this method.