Boundary integral formulation and two-dimensional fundamental solutions for dynamic behavior analysis of unsaturated soils

In this paper the coupled equations governing the dynamic behavior of unsaturated soils are derived based on the poromechanics theory within the framework of the suction-based mathematical model presented by Gatmiri (1997) [Gatmiri B. Analysis of fully coupled behavior of unsaturated porous medium under stress, suction and temperature gradient. Final report of CERMES-EDF, 1997] and Gatmiri et al. (1998) [Gatmiri B, Delage P, Cerrolaza M, UDAM: a powerful finite element software for the analysis of unsaturated porous media. Adv Eng Software 1998; 29(1): 29–43]. In this formulation, the solid skeleton displacements, water pressure and air pressure are presumed to be independent variables. The Boundary Integral formulations as well as fundamental solutions for such a dynamic u–pw–pa theory are presented in this paper for the first time. The boundary integral equations are derived via the use of the weighted residuals method in a way that permits an easy discretization and implementation in a Boundary Element code. Also, the associated two dimensional (2D) fundamental solutions for such deformable porous medium with linear elastic behavior are derived in Laplace transform domain using the method of Hörmander. Finally, some numerical results are presented to show the accuracy of the proposed solutions. The derived results are verified analytically by comparison with the previously introduced corresponding fundamental solutions in elastodynamic limiting case.

Research highlights
► Development of a set of coupled equations governing the dynamic behavior of unsaturated soil.
► The Boundary Integral formulations for such a dynamic model are derived.
► The associated 2D fundamental solutions for such deformable porous medium are derived in Laplace transform domain.
► The singular behavior of the fundamental solutions is studied in order to be able to determine the unknown boundary data.
► The derived results are verified analytically by comparison with the corresponding fundamental solutions in elastodynamic limiting case.