Modeling Surface Waves on a Partially Saturated Poroelastic Half-Space

Defining relations for a partially saturated porous Biot medium, written in the variables of displacements of the skeleton and pore pressures of the fillers, are considered. The initial system of partial differential equations includes five functions (a displacement vector and two pore pressures). The model of the material corresponds to a three-component medium. A system of equations in partial derivatives and boundary-value conditions are written in Laplace transform for time variable and in direct time with initial conditions. The boundary-value problem is analyzed using the method of boundary integral equations, their solutions being sought with the boundary-element method. The numerical scheme is based on using the Green-Bettie-Somigliana formula. Quadratic interpolation polynoms are taken as form functions in describing the boundary of the body. Unknown boundary fields are sought through nodal values in interpolation nodes. The element-by-element numerical integration uses Gauss method and an adaptive integration algorithm. The boundary-element schemes are constructed, based on the consistent approximation of the boundary functions and the collocation method. The solution of the formulated system of linear algebraic equations is sought using the block-type Gauss method. The boundary integral equation method in combination with the technique of searching a boundary-element solution is oriented at a dynamic problem of an isotropic homogeneous partially saturated poroelastic half-space. The time-stepping method for numerical inverting the Laplace transform is used to obtain the solution in the time domain.