The Boundary-Element Approach to Modeling the Dynamics of Poroelastic Bodies

The present paper is dedicated to dynamic behavior of poroelastic solids. Biot's model of poroelastic media with four base functions is employed in order to describe wave propagation process, base functions are skeleton displacements and pore pressure of the fluid filler. In order to study the boundary-value problem boundary integral equations (BIE) method is applied, and to find their solutions boundary element method (BEM) for obtaining numerical solutions. The solution of the original problem is constructed in Laplace transforms, with the subsequent application of the algorithm for numerical inversion. The numerical scheme is based on the Green-Betty-Somigliana formula. To introduce BE-discretization, we consider the regularized boundary-integral equation. The collocation method is applied. As a result, systems of linear algebraic equations will be formed and can be solved with the parallel calculations usage. Modified Durbin's algorithm of numerical inversion of Laplace transform is applied to perform solution in time domain. A problem of the three-dimensional poroelastic prismatic solid clamped at one end, and subjected to uniaxial and uniform impact loading and a problem of poroelastic cube with cavity subjected to a normal internal pressure are considered.