# Analysis of the Three-Dimensional Dynamic Problems by Using a New Numerical Method

Authors: Yao Rong; Yang Sun; LiQing Zhu; Xiao Xiao
Journal: Advances in Civil Engineering
 (2021)

---

## Abstract

The problems of the consolidation of saturated soil under dynamic loading are very complex. At present, numerical methods are widely used in the research. However, some traditional methods, such as the finite element method, involve more degrees of freedom, resulting in low computational efficiency. In this paper, the scaled boundary element method (SBFEM) is used to analyze the displacement and pore pressure response of saturated soil due to consolidation under dynamic load. The partial differential equations of linear problems are transformed into ordinary differential equations and solved along the radial direction. The coefficients in the equations are determined by approximate finite elements on the circumference. As a semianalytical method, the application of scaled boundary element method in soil-structure interaction is extended. Dealing with complex structures and structural nonlinearity, it can simulate two-phase saturated soil-structure dynamic interaction in infinite and finite domain, which has an important engineering practical value. Through the research, some conclusions are obtained. The dimension of the analytical problem can be reduced by one dimension if only the boundary surface is discretized. The SBFEM can automatically satisfy the radiation conditions at infinite distances. The 3D scaled boundary finite element equation for dynamic consolidation of saturated soils is not only accurate in finite element sense but also convenient in mathematical processing.


---

## Body

## 1. Introduction

Dynamic consolidation analysis of foundation is one of the most concerned problems, such as seismic design of large underground caverns, high arch dams, super high-rise buildings, and other large-scale industrial and their response characteristics under mechanical vibration (Chen and He [1], Lu et al. [2], and Jiang and Liang [3]).The free field formed by dynamic loads not only affects the dynamic response of long-span structures but also has feedback effect on the propagation of vibration waves by different components. Therefore, it is not reasonable to study the structure-constrained system as an isolator. It is necessary to consider the dynamic interaction between structure and foundation soil. For foundation soils in semi-infinite space, many literatures regard foundation soils as single-phase materials. When the foundation soil with high water content is liquefied by earthquake load, the pore pressure will increase with time, which will have a greater impact on the safety of the structure (Du and Wang [4] and Xia et al. [5]). Therefore, considering that the foundation soil is saturated soil with two-phase coupling of water and soil, the analysis of dynamic consolidation can better reflect the engineering practice.For infinite and semi-infinite foundation soils, the main numerical analysis methods are finite element method (FEM), boundary element method, and finite element-boundary element coupling method. The simulation of semi-infinite foundation soil is much more complicated compared to the finite element method mainly because of the large discrete range of foundation soil. In order to obtain enough calculation accuracy, the free degrees of the system increased inevitably. Although the computing power has been greatly improved, the computational units needed for some special cases, such as the accurate analysis of 3D stress wave propagation in anisotropic soils, are still difficult to achieve. Generally, the finite element discrete analysis is used to study foundations around the structure. Free or fixed boundary conditions are used to simulate infinite foundation soil. At the same time, in order to eliminate false reflection on artificial boundary, a transmission boundary model is established, for example, viscous and viscoelastic boundary (Zhang [6]), superposition boundary (Smith [7]), paraxial boundary (Engquist and Majda [8]), transient transmission boundary (Liao et al. [9] and Liao et al. [10]), and multidirectional and bi-asymptotic multidirectional transmission boundary (Keys [11] and Wolf and Song [12]). The main problems of finite element method based on the infinite boundary model, when it is used to analyze the dynamic consolidation of infinite foundation soil, are the inadequate accuracy of low-order boundary and the poor stability of high-order boundary. It does not have the accuracy in the sense of finite element. In other words, the numerical solution can converge to the solution when the discrete mesh is infinitely small.The wave attenuation factor is introduced into the shape function to describe the propagation of vibration waves along infinite distances for the infinite element method (Zhao and Valliappan [13] and Khaliliet al [14]), which is also a numerical method for infinite domain problems. However, the accuracy depends largely on the selection of radial shape function. When the order of the radial shape function is low, the infinite element must be placed far away for higher accuracy, which will lead to the expansion of the solution scale of the system. The higher order radial shape functions usually leads to ill-conditioned equations. In addition, similar to the local transmission boundary, the infinite element method is not accurate enough in the finite element sense. In dealing with these problems of infinite and semi-infinite domains, the boundary element method could satisfy the radiation conditions very well at infinity, and there is no reflection problem of artificial boundary (Bonnet [15] and Xu et al. [16]). However, the basic solutions needed to form the integral equation of boundary element method are generally complex, especially for anisotropic materials, and there is not even a basic analytical solution. On the contrary, the treatment of singular integral or even hypersingular integral of boundary integral equation is very difficult, which limits the further application of boundary element method. The FEM-BEM coupling method has been widely used in the analysis of structure-foundation interaction (Yazdchiet al. [17] and Du [18]). However, it is difficult to coordinate the coupled deformation between FEM and BEM for the interface between structure and foundation.Scaled boundary finite element method (SBFEM) is a boundary element method based on finite element method. Compared with the boundary element method, it does not need to solve the fundamental solution, so it can effectively deal with the problem of anisotropic media which is particularly complex and satisfies certain conditions. The scaled boundary finite element method (SBFEM) successfully satisfies Sommer field radiation condition by choosing appropriate similar centers. Therefore, the waves emitted by the source can only dissipate to infinity in the form of de-wave, but not return from infinity in the infinite domain. At present, SBFEM has been applied to the analysis of wave problems in time domain and frequency domain and to the solution of boundary dynamic stiffness matrix of infinite foundation (Wolf and Song [19, 20]). Previous results show that the SBFEM is very accurate and effective in dealing with most problems of infinite medium, anisotropic medium, and inhomogeneous change of materials. Considering that the foundation soil is a saturated soil with two-phase coupling of water and soil, the dynamic consolidation of 3D infinite saturated soil is analyzed by the SBFEM, which has not been reported in previous literatures.


## 2. Governing Equation

Taking 3D saturated soil areaV as the research object in Cartesian rectangular coordinates, the cross section is shown in Figure 1. If the field V is a finite field, the scale curvilinear coordinate system shown in Figure 2 is established for a unit Ve in the field, whose bottom surface is Se and the cone surface is Ae. For convenient of expression, the scaled center point C coincides with the origin point O of the Cartesian coordinate system. When the SBFEM is used, elements similar to the FEM are discretized on the Se boundary, which is on the bottom of the domain. η and ζ are the toroidal local coordinates of the elements on the Se boundary in the surface coordinate system, and ξ is the radial coordinates perpendicular to the boundary. η,ζ=±1 represents the edge of bottom Se. Radial coordinate 0≤ξ≤1. When ξ = 0 or 1, it represents the scale center point and the bottom of the finite field, respectively. The bottom s of scaled center point and finite field V are represented, respectively. When the position of the scaled center is determined, a one-to-one conversion relationship is formed between the Cartesian rectangular coordinate system and the local scaled boundary coordinate system, which is called the scaled boundary transformation.Figure 1
Cutaway view of three-dimensional saturated soil: (a) finite region and (b) unbounded media.
(a)(b)Figure 2
Transformation of global Cartesian coordinates to scaled boundary coordinates.For the 3D infinite saturated soil shown in Figure1(b), the scaled center point is selected outside the infinite domain. Similar to the finite field, the vertex of the element cone is the scaled center point C, and the bottom represents the interface (as shown in Figure 2). Then, the radial coordinates ξ≥1 and ξ = 1 correspond to the bottom of the infinite field. The whole infinite field is formed by assembling all the cones.In the local coordinate system of the scaled boundary, any point on the bottomSe of ξ = 1 can be expressed by n node coordinates of the discrete element on the bottom as follows:(1)xη,ζ=Nη,ζx=Nx,yη,ζ=Nη,ζy=Ny,zη,ζ=Nη,ζz=Nz,where the shape function is defined as Nη,ζ=N1η,ζN2η,ζN3η,ζ,….In the scaled coordinates, the relationship between pointr^=x^i+y^j+z^k in the computational element domain Ve and bottom Se boundary point r=xi+yj+zk can be expressed as follows:(2)x^ξ,η,ζ=ξxη,ζ,y^ξ,η,ζ=ξyη,ζ,z^ξ,η,ζ=ξzη,ζ.Then,(3)∂∂x^∂∂y^∂∂z^=J^−1∂∂ξ∂∂η∂∂ζ,where Jacobian matrix J^η,ζ=x^ξ,y^ξ,z^ξ,x^η,y^η,z^η,x^ζ,y^ζ,z^ζ,.According to Biot’s theory, the governing equations of consolidation of saturated soils include the dynamic equilibrium equation of soil skeleton and the continuous equation of fluid. In the Cartesian coordinate systemx^,y^,z^, the dynamic equilibrium equation expressed by the displacement (u=ux^,y^,z^=uxuyuzT) of soil skeleton in frequency domain is as follows:(4)LTσ′+mpf+ρω2u+ρb=0,where pf is pore water pressure, ω is frequency, and ρ=ρs+1−ϕρf. ρs and ρf are soil skeleton and fluid density, ϕ is porosity, b=bxbybzT, σ′=σx′σy′σz′τyzτxzτxyT and σ=σ′+mpf are unit volume force, effective stress, and total stress, respectively. m=111000T. Differential operator L is defined as(5)L=∂∂x^000∂∂z^∂∂y^0∂∂y^0∂∂z^0∂∂x^00∂∂y^∂∂y^∂∂x^0T.The continuity equation of fluid motion in saturated soils in frequency domain is as follows:(6)−kpf,ii+ω2kρfui,j+kρfbi,j+iωui,j+iωϕpfKf=0,where k and Kf are the permeability coefficient of soil and the bulk modulus of fluid, respectively.For isotropic saturated soils, effective stress, strain, and displacement vector satisfy Hooke's law:(7)σ′=Dε=DLu,where [D] is the elastic matrix of saturated soil material.According to equation (3), in the scaled curvilinear coordinate system ξ,η,ζ shown in Figure 2, the differential operator [L] is expressed as(8)L=b1∂∂ξ+1ξb2∂∂η+b3∂∂ζ,b1=gξJnxξ000nzξnyξ0nyξ0nzξ0nxξ00nzξnyξnxξ0T,b2=gηJnxη000nzηnyη0nyη0nzη0nxη00nzηnyηnxη0T,b3=gζJnxζ000nzζnyζ0nyζ0nzζ0nxζ00nzζnyζnxζ0T,where nxξ, nyξ, nzξ, nxη, nyη, nzη, nxζ,nyζ, and nzζ are the components of the outer surface unit normal (nξ, nη and nζ) of the bottom Se shown in Figure 2, respectively.Equation (8) shows that b1, b2, and b3 are independent of the radial coordinates ξ and satisfy the following formula:(9)Jb2η+Jb3ζ=−2Jb1.In the Cartesian coordinatesx^,y^,z^, the tangential force of a surface with a normal direction of n=nxi+nyj+nzk is(10)τ=τxτyτz=nx000nzny0ny0nz0nx00nznynx0σ.From equations (3), (8), and (10), it can be seen that the tangential force on the surface η,ζ, ζ,ξ, and ξ,η of any given value ξ, η, and ζ is as follows:(11)tξ=Jgξb1Tσ,tη=Jgηb2Tσ,tζ=Jgζb3Tσ.For the computational element domainVe shown in Figure 2, the displacement and pore pressure of discrete nodes are expressed by a shape function similar to the coordinates of nodes on the bottom boundary Se of ξ=1. Then, on a surface Sξ with any radial coordinate ξ, displacement and pore pressure can be expressed as(12)u=uξ,η,ζ=Nuη,ζuξ,pf=pfξ,η,ζ=Npη,ζpfξ.From equations (7), (8), and (12), it can be seen that the effective stress in the scaled curvilinear coordinates is as follows:(13)σ′=Db1Tuξ+1ξb2Tuη+b3Tuζ=DBu1uξξ+1ξBu2uξ,(14)σξ′=DBu1uξξξ+1ξBu2uξξ−1ξ2Bu2uξ,whereBu1T=b1Nu and Bu2=b2Nu,η+b3Nu,ζ. Therefore, it can be concluded that Bu1 and Bu2 are only functions of circumferential coordinates, independent of radial coordinates η,ζ.By using equations (3), (8), and (12), the differential equations (4) and (6) of saturated soil consolidation in Cartesian coordinate system x^,y^,z^ can be transformed into the scaled curvilinear coordinates, and it can be expressed as follows:(15)b1Tσξ′+mpf,ξ+1ξb2Tση′+mpf,η+b3Tσζ′+mpf,ζ+ρω2u+ρb=0,(16)−kmTb1∂∂ξb1T∂∂ξ+1ξb2T∂∂η+b3T∂∂ζ+1ξb2T∂∂ηb1T∂∂ξ+1ξb2T∂∂η+b3T∂∂ζ+b3T∂∂ζb1T∂∂ξ+1ξb2T∂∂η+b3T∂∂ζmpf+ω2kρfmTb1Tuξ+1ξb2Tuη+b3Tuζ+kρfmTb1Tbξ+1ξb2Tbη+b3Tbζ+iωmTb1Tuξ+1ξb2Tuη+b3Tuζ+iωϕpfKf=0.


## 3. Application of Weighted Residual Method

Galerkin weighted residual method is adopted in computing unit domainVe. For convenience of expression, the superscript e of the presentation unit is omitted. Interpolated shape functions of the same form as the displacement of discrete nodes are chosen for the weight function, such as wu=wuξ,η,ζ=Nuη,ζwuξ. For the dynamic equilibrium equation (12) of saturated consolidation soil, it can be obtained by the Galerkin method in the computational domain:(17)∫VwTb1Tσξ′+mpf,ξdV+∫V1ξwTb2Tση′+mpf,η+b3Tσζ′+mpf,ζdV+ρω2∫VwTudV+ρ∫VwTbdV=0.For the second item in the above equation, by using the partial integral and substituting equations (9) and (11), we can obtain(18)I=∫01ξ∫ΓξwTtζgζdη+tηgηdζ−∫Sξ−2wTb1T+wηTb2T+wζTb3Tσ′+mpfJdηdζdξ.By substituting equation (18) into equation (17), the following equation can be obtained:(19)∫01ξ2∫SξwTb1Tσξ′+mpf,ξJdηdζ+ξ∫ΓξwTtζgζdη+tηgηdζ−ξ∫Sξ−2wTb1T+wηTb2T+wζTb3Tσ′+mpfJdηdζ+ρω2ξ2∫SξwTuJdηdζ+ρξ2∫SξwTbJdηdζdξ=0.Then, we substitutedwT=wξTNuT into the upper equation, and if equality is established for the upper equation in the whole integral domain, then(20)wTξ2∫SξNuTb1Tσξ′Jdηdζ+wTξ2∫SξNuTb1TmNppfξξJdηdζ+wTξ∫Sξ2NuTb1T−NuηTb2T−NuζTb3Tσ′Jdηdζ+wTξ∫Sξ2NuTb1T−NuηTb2T−NuζTb3TmNppfξJdηdζ+wTξ∫ΓξNuTtζgζdη+tηgηdζ+wTω2ξ2∫SξNuTρNuuξJdηdζ+wTρξ2∫SξNuTbJdηdζ=0.After multiplyingwT−1 at both sides of the upper equation and substituting Bu1 and Bu2, we can obtain(21)ξ2∫SξBu1Tσξ′Jdηdζ+ξ2∫SξBu1TmNppfξξJdηdζ+ξ∫Sξ2Bu1T−Bu2Tσ′Jdηdζ+ξ∫Sξ2Bu1T−Bu2TmNppfξJdηdζ+ω2ξ2∫SξNuTρNuuξJdηdζ+ξFtξ+ξ2Fbξ=0,where Ftξ=∫ΓξNuη,ζTtζgζdη+tηgηdζ and Fbξ=ρ∫SξNuη,ζTbJdηdζ。.By using equations (13), (14), and (21), the following results can be obtained:(22)ξ2∫SξBu1TDBu1Jdηdζuξξξ+ξ2∫SξBu1TmNpJdηdζpfξξ+ξ∫SξBu1TDBu2Jdηdζuξξ+ξ∫Sξ2Bu1T−Bu2TDBu1Jdηdζuξξ−∫SξBu2TDBu2Jdηdζuξ+ω2ξ2∫SξNuTρNuJdηdζuξ+ξ∫Sξ2Bu1T−Bu2TmNpJdηdζpfξ+ξFtξ+ξ2Fbξ=0.For the continuity equation of fluid in the consolidation equation of saturated soil, the weight function is expressed by the interpolation function which is the same as the pore pressure. Galerkin's method is used in the calculation domain. Similar to the above deduction, the following can be obtained from equation (16):(23)∫V−kwTmTb1b1Tmpf,ξξ+1ξb1b2Tmpf,ξη+b1b3Tmpf,ξζ−1ξ2b1b2Tmpf,η+b1b3Tmpf,ζ+1ξb2b1ηTmpf,ξ+b1Tmpf,ξη+1ξb3b1ζTmpf,ξ+b1Tmpf,ξζ+1ξ2b2b2ηTmpf,η+b2Tmpf,ηη+1ξ2b2b3ηTmpf,ζ+b3Tmpf,ηζ+1ξ2b3b2ζTmpf,η+b2Tmpf,ηζ+1ξ2b3b3ζTmpf,ζ+b3Tmpf,ζζ+iω+ω2kρfwTmTb1Tu,ξ+1ξb2Tuη+b3Tuζ+kρfwTmTb1Tbξ+1ξb2Tbη+b3Tbζ+iωϕwTpfKfdV=0.By using the partial integral and substitutingwpT=wpξTNpT into the upper equation, the following results can be obtained:(24)−ξ2wT∫SξBp1TkBp1Jdηdζpfξξξ+ξwT∫SξBp1TkBp2Jdηdζpfξξ+ξwT∫ΓξBp2ATkBp1Jdζ+∫ΓξBp2BTkBp1Jdηpfξξ+2ξwT∫SξBp1TkBp1Jdηdζpfξξ−ξwT∫SξBp2TkBp1Jdηdζpfξξ+wT∫ΓξBp2ATkBp2Jdζ+∫ΓξBp2BTkBp2Jdηpfξ−wT∫SξBp2TkBp2Jdηdζpfξ+ξ2iω+ω2kρf∫SξNpTmTBu1Jdηdζuξξ+ξiω+ω2kρf∫ΓξBp2ATNuJdζ+Bp2BTNuJdηuξ−ξiω+ω2kρf∫SξBp2TNuJdηdζuξ+2ξiω+ω2kρf∫SξNpTmTBu1Jdηdζuξ+ξ2kρf∫SξBp1TbξJdηdζ+2ξkρf∫SξBp1TbJdηdζ+ξ2iω∫SξwTϕKfNpJdηdζpfξ=0,where Bp1=b1TmNp, Bp2AT=NpTmTb2, Bp2=b2TmNpη+b3TmNpζ, and Bp2BT=NpTmTb3.The dynamic consolidation equation of saturated soil considering pore water pressure can be obtained by equations (22) and (24). The scaled boundary finite element equation expressed by pore water pressure and displacement is as follows:(25)E0ξ2uξξξ+E3ξ2pfξξ+2E3−E4ξpfξ+2E0−E1+E1Tξuξξ−E2uξ+ω2M0ξ2uξ+ξ2Fbξ+ξFtξ=0,(26)E5ξ2pfξξξ+2E5−E6+E6T+Ft1ξpfξξ+−E7+Ft2pfξ−iωM1ξ2pfξ+Ft1ξpfξξ+−E7+Ft2×pfξ−iωM1ξ2pfξ−iω+ω2kρfE3Tξ2uξξ−iω+ω2kρf2E3T−E8+Ft3ξuξ−kρfξ2Fb1ξ+2ξFb2ξ−ξFb3ξ+ξFb4ξ=0,where(27)M0=∫SξNuTρNuJdηdζ,M1=∫SξNpTϕKfNpJdηdζ,E0=∫SξBu1TDBu1Jdηdζ,E1=∫SξBu1TDBu1Jdηdζ,E2=∫SξBu2TDBu2Jdηdζ,E3=∫SξBu1TmNpJdηdζ,E4=∫SξBu2TmNpJdηdζ,E5=∫SξBp1TkBp1Jdηdζ,E6=∫SξBp2TkBp1Jdηdζ,E7=∫SξBp2TkBp2Jdηdζ,E8=∫SξBp2TNuJdηdζ,Fb3ξ=∫SξBp2TbJdηdζ,Fb1ξ=∫SξBp1TbξJdηdζ,Fb2ξ=∫SξBp1TbJdηdζ,Fb4ξ=∫ΓξBp2ATbJdζ+Bp2BTbJdη,Ft1=∫ΓξBp2ATkBp1Jdζ+∫ΓξBp2BTkBp1Jdη,Ft2=∫ΓξBp2ATkBp2Jdζ+∫ΓξBp2BTkBp2Jdη,Ft3=∫ΓξBp2ATNuJdζ+Bp2BTNuJdη.The displacement and pore pressure vectors are expressed by generalized displacement vectorsU=uξpfξ.The second-order differential equations can be obtained by equations (25) and (26):(28)ξ2E000E5Uξξξ−ξ2iω+ω2kρfE3T00−E3−ξ2E0−E1+E1T002E5−E6+E6T+Ft1Uξξ+ξ2ω2M000−iωM1−ξiω+ω2kρf2E3T−E8+Ft300E4−2E3−E200E7−Ft2Uξ=−ξ2Fbξ−kρfFb1ξ−ξFtξ−kρf2Fb2ξ−Fb3ξ+Fb4ξ.


## 4. Numerical Analysis and Discussion

The displacement and pore pressure responses of 3D half-space surface subjected to vertical concentrated harmonic loads are investigated. In the calculation, the saturated soil layer is clay layer, and the relevant parameters in [16] are used. Some soil parameters adopted in the present study are obtained from field investigation and laboratory test, such as ρs=2000kg/m3, ρf=1000kg/m3, ϕ=0.4, Kf=2.0×109N/m2, υ=0.4, k=1.0×10−9m/s, and μ=1.0×107Pa.Using the method presented in this paper, the 3D half-space is divided into two parts. A hemisphere with a radius of 2R is a structure and the rest is a semi-infinite space. According to symmetry, 1/4 half space is taken for analysis, as shown in Figure3. Two parts of 3D half-space are analyzed by scaled boundary finite element method. The structure half-space interface is discretized by three 8-node two-dimensional elements, and the structure part is discretized by eighteen 8-node two-dimensional elements.Figure 3
SBFEM discretization of 1/4 unbounded half-space.In order to verify the accuracy of this method, ifρf=0 and n=0, the two-phase medium degenerates into a single-phase medium. By using the method presented in this paper, the solution of the degenerate 3D linear elastic dynamic Boussinesq problem can be obtained. Figure 4 shows the relationship between vertical displacement u˜z=uzrμ/P and dimensionless frequency ω˜=ρ/μωr, where r is the horizontal distance between concentrated load P and surface calculation point. It can be seen from Figure 4(a) that the shape of the curve for the present solution is very similar to the solution in the literatures (Banerjee and Mamoon [21] and Kobayashi and Nishimura [22]). The maximum value of the solution in the present solution is 0.092, and the minimum value is 0.03, which is consistent with the solution in the literature. And, the difference between them is in the range of 0.002. It can be seen from Figure 4(b) that the solution in present solution has a decreasing trend, and the curve shape is relatively smooth. The maximum value decreased from −0.017 to −0.092. By comparing the present solution in this paper with that in the literature, it can be concluded that the fluctuation of the solution in present solution is relatively small. The numerical results in this paper are in good agreement with those in literature (Banerjee and Mamoon [21] and Kobayashi and Nishimura [23]), which verifies the accuracy of the proposed method.Figure 4
Compared present solutions with the reference’s results: (a) real part of vertical displacement and (b) imaginary part of vertical displacement.
(a)(b)Figure5 shows the relationship between vertical displacement u˜z=uzrμ/P of saturated and elastic soils in 3D half space and dimensionless frequency ω˜=ρ/μωr. The results show that the dynamic responses of saturated soil and elastic soil are different at different frequencies. It can be seen from the figure that, with the increase of ω˜, the real part of the vertical displacement u˜z decreases, while the imaginary part of the vertical displacement u˜z appears at an extreme value near ω˜ = 1.3.Figure 5
The vertical displacement for the 3D half-space varied with the frequency: (a) real part of vertical displacement and (b) imaginary part of vertical displacement.
(a)(b)Figure6 shows the relationship between pore pressure p˜f at point B on the structure-soil interface (2R away from the surface) and dimensionless frequency ω˜ at different permeability coefficients k=1.0×10−7m/s, 1.0×10−9m/s, and 1.0×10−11m/s. It can be seen from the figure that the pore pressure amplitude increases with the increase of permeability coefficient of saturated soil.Figure 6
The pore pressure for the interface between the half spec and structure varied with the frequency: (a) real part of vertical displacement and (b) imaginary part of vertical displacement.
(a)(b)


## 5. Conclusion

In this paper, the scaled boundary element method for single-phase elastic media is extended to the dynamic consolidation of 3D saturated soils with water-soil two-phase coupling. The application of scaled boundary element method in soil-structure interaction is extended. From the theoretical analysis in this paper, it can be seen that(1)
Not the same as single-phase media, the displacement and stress matrices are not only formed in the scaled boundary finite element equation of saturated soil consolidation but also the pore water pressure effect is coupled in the scaled boundary finite element matrix. According to Biot's theory, the governing differential equation of saturated soil consolidation is solved in Cartesian rectangular coordinate system, and the proportional boundary finite element equation in the form of pore pressure and displacement is obtained.(2)
The scaled boundary of dynamic consolidation of saturated soil also has unique advantages. The dimension of the analytical problem can be reduced by one dimension if only the boundary surface is discretized. In addition, similar to the boundary element method, it can automatically satisfy the radiation conditions at infinite distances. However, the method in this paper does not need to solve the fundamental solution and does not involve singular integrals.(3)
As a semi analytical method, the 3D scaled boundary finite element equation for dynamic consolidation of saturated soils is not only accurate in finite element sense but also convenient in mathematical processing. Especially, in dealing with complex structures and structural nonlinearity, it can simulate two-phase saturated soil-structure dynamic interaction in infinite and finite domain, which has important engineering practical value.



---

*Source: 5555575-2021-05-05.xml*