Analytical layer‐element method for non‐axisymmetric consolidation of multilayered soils

A numerically efficient and stable method is developed to analyze Biot's consolidation of multilayered soils subjected to non‐axisymmetric loading in arbitrary depth. By the application of a Laplace–Hankel transform and a Fourier expansion, the governing equations are solved analytically. Then, the analytical layer‐element (i.e. a symmetric stiffness matrix) describing the relationship between generalized displacements and stresses of a layer is exactly derived in the transformed domain. Considering the continuity conditions between adjacent layers, the global stiffness matrix of multilayered soils is obtained by assembling the inter‐related layer‐elements. Once the solution in the Laplace–Hankel transformed domain that satisfies the boundary conditions has been obtained, the actual solution can be derived by the inversion of the Laplace–Hankel transform. Finally, numerical examples are presented to verify the theory and to study the influence of the layered soil properties and time history on the consolidation behavior. Copyright © 2011 John Wiley & Sons, Ltd.     

Consolidation analysis of saturated multi‐layered soils with anisotropic permeability caused by a point sink

This paper presents the analytical layer‐element method to analyze the consolidation of saturated multi‐layered soils caused by a point sink by considering the anisotropy of permeability. Starting from the governing equations of the problem, the solutions of displacements and stresses for a single soil layer are obtained in the Laplace–Hankel transformed domain. Then, the analytical layer‐element method is utilized to further derive the solutions for the saturated multi‐layered soils in the transformed domain by combining with the boundary conditions of the soil system and continuity conditions between adjacent layers. The actual solutions in the physical domain can be acquired by the inversion of Laplace–Hankel transform. Numerical results are carried out to show the accuracy and stability of the proposed method and evaluate the influence of sink depth and anisotropic permeability on excess pore pressure and surface settlement. Copyright © 2011 John Wiley & Sons, Ltd.     

层状地基上弹性地基梁的有限元-边界元耦合分析

将地基视为多层各向同性弹性体,对Euler-Bernoulli梁进行有限单元离散分析,对地基-梁接触面采用边界积分法求解,根据地基-梁接触面的竖向位移协调和光滑接触条件,应用有限元与边界元耦合的方法推导出各向同性成层弹性地基上的Euler-Bernoulli梁的半数值半解析解。基于文中理论,编制了相应的程序,通过与现有文献对比,验证了文中理论的正确性,对比分析了分层地基与等效均匀地基两种模型上的弹性地基梁。研究结果表明,分层地基与等效均匀地基两种模型上弹性地基梁性状差异较大,实际弹性地基梁计算中应采用分层地基模型。     

层状地基中轴向受荷单桩的边界元法分析

以层状地基内部作用一竖向集中力时的广义Mindlin解作为边界单元法的基本解,对层状地基中的轴向受荷单桩进行了分析,对基本解的奇异性处理方法进行了改进。考虑了桩的可压缩性和长径比对桩-土荷载传递规律和沉降特性的影响,编制了计算程序,并进行了数值分析和计算。结果表明,该方法具有较快的计算速度和良好的计算精度。     

Deformation analysis of tunnel excavation below existing pipelines in multi‐layered soils based on displacement controlled coupling numerical method

The effect of tunneling on surrounding environments, especially on existing buried pipelines is a problem that engineers designing and practicing in urban geotechnical environments encounter more frequently than in the past. However, previous studies are usually based on the assumption that the soil is homogeneous. How to reflect soil stratification is the main focus for the problem of tunneling in multi‐layered soils. A displacement controlled coupling numerical method is presented for the displacement analysis of tunnel excavation below existing pipelines in multi‐layered soils. On the basis of the layered soil model, to consider the soil nonhomogeneous characteristic, the finite element method and boundary element method are coupled to simulate the deformation of existing pipelines induced by tunneling. The solutions indicate that good agreements are obtained between the proposed coupling numerical method and the commercial software. The accuracy of the proposed numerical method is better than the two stages method based on the existing closed‐form solutions. Moreover, the results discussed in this paper show that the error obtained by the previous method of weighted average on the basis of homogeneous half space converted from layered soils is not negligible for the obvious difference of elastic parameters among successive layers. Copyright © 2012 John Wiley & Sons, Ltd.     

Influence of a point sink on a pile group in saturated multilayered soils with anisotropic permeability

The behavior of a pile group is solved using the finite element method, and the fundamental solution of saturated multilayered soils with anisotropic permeability is obtained by the analytical layer element method. Based on the supposition of no slip occurring at the pile‐soil interface, the governing equations of the interaction between the pile group and the soils due to a point sink are established in the Laplace‐Hankel transformed domain by considering the pile‐soil compatibility condition. Numerical results are presented to study the effect of point sink pumping, the properties of soils, and the geometries of piles on the behavior of the pile group.     

Iteratively coupled solution strategies for a four‐field mixed finite element method for poroelasticity

In this paper, we consider numerical algorithms for modeling of the time‐dependent coupling between the fluid flow and deformation in elastic porous media. Here, we employ a four‐field formulation which uses the total stress, displacement, flux, and pressure as its primary variables and satisfies Darcy's law and linear elasticity in mixed weak form. We present four different iteratively coupled methods, known as drained, undrained, fixed‐strain, and fixed‐stress splits, in which the diffusion operator is separated from the elasticity operator and the two subproblems are solved in a staggered way while ensuring convergence of the solution at each time step. A‐priori convergence results for each iterative coupling which differs from those found when using a traditional two‐field or three‐field formulation are presented. We also present some numerical results to support the convergence estimates and to show the accuracy and efficiency of the algorithms. Copyright © 2016 John Wiley & Sons, Ltd.     

A FIC‐based stabilized mixed finite element method with equal order interpolation for solid–pore fluid interaction problems

A new mixed displacement‐pressure element for solving solid–pore fluid interaction problems is presented. In the resulting coupled system of equations, the balance of momentum equation remains unaltered, while the mass balance equation for the pore fluid is stabilized with the inclusion of higher‐order terms multiplied by arbitrary dimensions in space, following the finite calculus (FIC) procedure. The stabilized FIC‐FEM formulation can be applied to any kind of interpolation for the displacements and the pressure, but in this work, we have used linear elements of equal order interpolation for both set of unknowns. Examples in 2D and 3D are presented to illustrate the accuracy of the stabilized formulation for solid–pore fluid interaction problems. Copyright © 2016 John Wiley & Sons, Ltd.     

Three-dimensional analysis of embankments on soft soils incorporating vertical drains by finite element method

Three-dimensional behaviour of an embankment on soft soils incorporating vertical drains is analysed by a numerical model based on the finite element method. The model, which incorporates the Biot’s consolidation theory (coupled formulation of the flow and equilibrium equations) and constitutive relations simulated by the p–q–θ critical state model, is applied on both the embankment with vertical drains (three-dimensional analysis) and the same problem but without vertical drains (two-dimensional analysis). Special emphasis is given to the analysis, during and after the construction period, of the excess pore pressures, settlements, horizontal displacements and stress levels.     

A three‐phase thermo‐hydro‐mechanical finite element model for freezing soils

Artificial ground freezing (AGF) is a commonly used technique in geotechnical engineering for ground improvement such as ground water control and temporary excavation support during tunnel construction in soft soils. The main potential problem connected with this technique is that it may produce heave and settlement at the ground surface, which may cause damage to the surface infrastructure. Additionally, the freezing process and the energy needed to obtain a stable frozen ground may be significantly influenced by seepage flow. Evidently, safe design and execution of AGF require a reliable prediction of the coupled thermo‐hydro‐mechanical behavior of freezing soils. With the theory of poromechanics, a three‐phase finite element soil model is proposed, considering solid particles, liquid water, and crystal ice as separate phases and mixture temperature, liquid pressure, and solid displacement as the primary field variables. In addition to the volume expansion of water transforming into ice, the contribution of the micro‐cryo‐suction mechanism to the frost heave phenomenon is described in the model using the theory of premelting dynamics. Through fundamental physical laws and corresponding state relations, the model captures various couplings among the phase transition, the liquid transport within the pore space, and the accompanying mechanical deformation. The verification and validation of the model are accomplished by means of selected analyses. An application example is related to AGF during tunnel excavation, investigating the influence of seepage flow on the freezing process and the time required to establish a closed supporting frozen arch. Copyright © 2013 John Wiley & Sons, Ltd.     

A new development of the scaled boundary finite element method for wave motion in layered half-space

The scaled boundary finite element method (SBFEM) developed by Wolf and Song has shown certain parallels to the finite element method (FEM) and boundary element method (BEM). Because of its semi-analytical nature, SBFEM is particularly suitable for the analysis of wave propagation in unbounded domains. This paper makes a certain modification of the standard SBFEM. A new idea of scaling surface instead of a scaling center is introduced to formulate the governing SBFE equations for the analysis of wave propagation in multilayered half-space, which leads to simplifying the modeling and saving considerably the computational effort. In addition, by employing the proposed approach, some problems encountered in engineering practice, which are difficult to deal with by the conventional SBFEM, for example, 3D foundation impedance on half-space with irregular geographical features, can be effectively solved. The proposed approach also helps to simplify the solution of shell structures. Numerical examples are provided to validate the accuracy and efficiency of the proposed approach.     

饱和多孔介质中水力-力学-传质耦合过程的混合有限元法

对饱和多孔介质提出了一个含溶混污染物输运（传质）过程的混合元方法,其中污染物输运过程数学模型包含了对流、机械逸散、分子弥散和吸附等机制。固相位移、应变和有效应力,孔隙水压力、压力空间梯度和Darcy速度,污染物浓度、浓度空间梯度和浓度流量在单元内均为独立变量分别插值。基于胡海昌-Washizu三变量广义变分原理,结合可以滤掉虚假振荡的特征线方法,推导出饱和土中水力-力学-传质耦合问题控制方程的单元弱形式,并导出了混合元计算公式。数值模拟证明了所提出的方法可以提供与传统4点积分方案同样精度,同时能够提高计算效率。     

A high‐frequency open boundary for transient seepage analyses of semi‐infinite layers by extending the scaled boundary finite element method

A high‐frequency open boundary has been developed for the transient seepage analyses of semi‐infinite layers with a constant depth. The scaled boundary finite element equation of pore water pressure is formulated first in the frequency domain. With the eigenvalue problem, the equation can be decoupled into modal equations whose modal dynamic permeability equation can be determined. The continued fraction technique is adopted to formulate the continued fraction solution in the frequency domain. All constants in the solution are determined recursively at the high‐frequency limit. By introducing auxiliary variables and the continued fraction solution to the relationship between the prescribed seepage flow and the pore water pressure in the frequency domain, the open boundary condition is obtained. After transformed to the time domain, the open boundary condition is expressed as a system of fractional differential equations. No convolution integral is required. The accuracy of the analysis results increases with the increasing orders of continued fraction. Copyright © 2015 John Wiley & Sons, Ltd.     

Vertical harmonic vibration of piled raft foundations in layered soils

This paper develops a method to analyze the piled raft foundation under vertical harmonic load. This method takes into account the interactions among the piles, soil, and raft. The responses of the piles and raft are formulated as a series of equations in a suitable way and that of layered soils is simulated with the use of the analytical layer‐element method. Then, according to the equilibrium and continuity conditions at the piles–soil–raft interface, solutions for the piled raft systems are obtained and further demonstrated to be correct through comparing with the existing results. Finally, some examples are given to investigate the influence of the raft, pile length‐diameter ratio, and layering on the response of the piled raft foundations. Copyright © 2017 John Wiley & Sons, Ltd.     

Semi‐analytical elastostatic analysis of unbounded two‐dimensional domains

Unbounded plane stress and plane strain domains subjected to static loading undergo infinite displacements, even when the zero displacement boundary condition at infinity is enforced. However, the stress and strain fields are well behaved, and are of practical interest. This causes significant difficulty when analysis is attempted using displacement‐based numerical methods, such as the finite‐element method. To circumvent this difficulty problems of this nature are often changed subtly before analysis to limit the displacements to finite values. Such a process is unsatisfactory, as it distorts the solution in some way, and may lead to a stiffness matrix that is nearly singular. In this paper, the semi‐analytical scaled boundary finite‐element method is extended to permit the analysis of such problems without requiring any modification of the problem itself. This is possible because the governing differential equations are solved analytically in the radial direction. The displacement solutions so obtained include an infinite component, but relative motion between any two points in the unbounded domain can be computed accurately. No small arbitrary constants are introduced, no arbitrary truncation of the domain is performed, and no ill‐conditioned matrices are inverted. Copyright © 2002 John Wiley & Sons, Ltd.     

Finite element modelling of thick plates on two‐parameter elastic foundation

This paper is intended to give some information about how to build a model necessary for bending analysis of rectangular and circular plates resting on a two‐parameter elastic foundation, subjected to combined loading and permitting various types of boundary conditions. The formulation of the problem takes into account the shear deformation of the plate and the surrounding interaction effect outside the plate. The numerical model based on an 18‐node zero‐thickness isoparametric interface element interacting with a thick Reissner–Mindlin plate element with three degrees of freedom at each of the nine nodes, which enforce C⁰ continuity requirements for the displacements and rotations of the midsurface, is proposed. Stiffness matrices of a special interface element are superimposed on the global stiffness matrix to represent the stiffening elastic foundation under and beyond the plate. Some numerical examples are given to illustrate the advantages of the method presented. Copyright © 2001 John Wiley & Sons, Ltd.     

A superposition scheme to obtain fundamental boundary element solutions in multi‐layered elastic media

A superposition scheme is proposed to obtain a fundamental solution for boundary elements in multi‐layered elastic media. A three‐layered elastic region is obtained by superposing two sets of bonded half‐planes and subtracting one infinite plane. Therefore, the solution for an element in a layered media can be expressed in terms of bonded half‐plane solutions and an infinite‐plane solution. The major advantages of this superposition scheme are: (1) it is unnecessary to introduce elements at the interface, (2) it can be extended to higher‐order element, and (3) it may be applicable to three dimensions easily. The accuracy and performance of the developed model is illustrated by two examples. For the problem of a pressurized two‐dimensional crack within a three‐layered system, the comparison with other numerical results shows the model is quite accurate and efficient. The model is also used for a study of a practical two‐dimensional mining problem in South Africa, i.e. stoping through a dyke with material properties different from the host rock. Copyright © 2000 John Wiley & Sons, Ltd.     

Solute transfer during consolidation based on a solid‐fluid‐solute coupling model

The migration of contaminant through soil is usually modeled using the advection‐dispersion equation and assumes that the porous media is stationary without introducing a constitutive equation to represent soil structure. Consequently, time‐dependent deformation induced by soil consolidation or physical remediation is not considered, despite the need to consider these variables during planning for the remediation of contaminated ground, the prediction of contaminated groundwater movement, and the design of engineered landfills. This study focuses on the numerical modeling of solute transfer during consolidation as a first step to resolve some of these issues. We combine a coupling theory‐based mass conservation law for soil‐fluid‐solute phases with finite element modeling to simulate solute transfer during deformation and groundwater convection. We also assessed the sensitivity of solute transfer to the initial boundary conditions. The modeling shows the migration of solute toward the ground surface as a result of ground settlement and the dissipation of excess pore water pressure. The form of solute transport is dependent on the ground conditions, including factors such as the loading schedule, contamination depth, and water content. The results indicate that an understanding of the interaction between coupling phases is essential in predicting solute transfer in ground deformation and could provide an appropriate approach to ground management for soil remediation.     

Analysis of pile‐raft foundations under complex loads in layered soils

Considering there is hardly any concerted effort to analyze the pile‐raft foundations under complex loads (combined with vertical loads, horizontal loads and moments), an analysis method is proposed in this paper to estimate the responses of pile‐raft foundations which are subjected to vertical loads, horizontal loads and moments in layered soils based on solutions for stresses and displacements in layered elastic half space. Pile to pile, pile to soil surface, soil surface to pile and soil surface to soil surface interactions are key ingredients for calculating the responses of pile‐raft foundations accurately. Those interactions are fully taken into account to estimate the responses of pile‐raft foundations subject to vertical loads, horizontal loads and moments in layered soils. The constraints of the raft on vertical movements, horizontal movements and rotations of the piles as well as the constraints of the raft on vertical movements and horizontal movements of the soils are considered to reflect the coupled effect on the raft. The method is verified through comparisons with the published methods and FEM. Then, the method is adopted to investigate the influence of soil stratigraphy on pile responses. The study shows that it is necessary to consider the soil non‐homogeneity when estimating the responses of pile‐raft foundations in layered soils, especially when estimating the horizontal responses of pile‐raft foundations. The horizontal loads and the moments have a significant impact on vertical responses of piles in pile‐raft foundations, while vertical loads have little influence on horizontal responses of piles in pile‐raft foundations in the cases of small deformations. The proposed method can provide a simple and useful tool for engineering design. Copyright © 2013 John Wiley & Sons, Ltd.