Three‐dimensional finite element analyses of passive pile behaviour

Piles may be subjected to lateral soil pressures as a result of lateral soil movements from nearby construction‐related activities such as embankment construction or excavation operations. Three‐dimensional finite element analyses have been carried out to investigate the response of a single pile when subjected to lateral soil movements. The pile and the soil were modelled using 20‐node quadrilateral brick elements with reduced integration. For compatibility between the soil–pile interface elements, 27‐node quadrilateral brick elements with reduced integration were used to model the soil around the pile adjacent to the soil–pile interface. A Mohr–Coulomb elastic–plastic constitutive model with large‐strain mode was assumed for the soil. The analyses indicate that the behaviour of the pile was significantly influenced by the pile flexibility, the magnitude of soil movement, the pile head boundary conditions, the shape of the soil movement profile and the thickness of the moving soil mass. Reasonable agreement is found between some existing published solutions and those developed herein. Copyright © 2005 John Wiley & Sons, Ltd.     

Three‐dimensional finite element analysis of lined tunnels

This paper describes finite element procedures that have been developed to model the ground movements that occur when a shallow tunnel is installed in a clay soil. This study is part of a wider project concerned with the development of new methods to predict the likely extent of damage to surface structures caused by nearby shallow tunnelling. This particular paper, however, is concerned only with the numerical model of tunnel installation. The structural liner is an important component of this tunnel installation model; two different ways of modelling the liner (based on continuum elements and shell elements) are discussed in the paper. A test problem consisting of the installation of a lined tunnel in an elastic continuum is used to investigate the merits of these different approaches. When continuum elements are used to model the liner, the numerical results agree well with an analytical solution to the problem. When shell elements are used to model the liner, however, the results were found to be significantly influenced by the particular formulation adopted for the shell elements. Example analyses, involving incremental tunnel construction in a clay soil where the soil is modelled using a kinematic hardening plasticity model, are described. These analyses confirm that a thin layer of continuum elements may be used, satisfactorily, to model tunnel linings in a soil–structure interaction analysis of this sort. Copyright © 2001 John Wiley & Sons, Ltd.     

Three‐dimensional finite element modeling of tire/ground interaction

Tire/ground interaction has been an important issue in terramechanics, transportation and pavement engineering. Characterization of tire/ground interaction has been majorly investigated based on empirical approaches and field tests. So far very few dynamic tire/ground interaction models have been presented. This paper presents a three‐dimensional finite element model for tire/ground interaction. The rubber material is modeled as nearly incompressible finite strain hyperelasticity instead of being modeled as a rigid wheel. The tire is also modeled as bias type and steel rebars are embedded within the tire structure. The tire model is combined with the ground model to form a tire/ground interaction model. Both single tire/ground interaction and four‐tire vehicle/pavement interaction models are presented, which allow us to investigate a lot of issues easily and accurately, such as compaction, pavement response and pavement damage. Numerous simulations are carried out to demonstrate that the dynamic tire/terrain interaction model can be used to predict the dynamic ground response due to moving vehicle. Different tire rolling conditions can be easily incorporated into the tire/ground interaction model, which further substantiates the broad application of the model in transportation and pavement engineering. Copyright © 2011 John Wiley & Sons, Ltd.     

Three‐dimensional responses of buried corrugated pipes and ANN‐based method for predicting pipe deflections

This study investigated localized responses, such as circumferential stresses, on corrugation and pipe deflections. Also, this study examined the effect of corrugation geometry on the overall and localized response of corrugated pipes with refined three‐dimensional modeling of the entire soil–pipe interaction system, including corrugation. To investigate the availability of the traditional two‐dimensional method, the results from the three‐dimensional finite element method (FEM) were compared with those from the two‐dimensional FEM. The soil–pipe modeling techniques of this study were verified by comparing the FEM results by Utah State University and analytical results. An artificial neural network (ANN)‐based model to predict vertical deflections of buried corrugated pipes was developed to overcome the shortcomings of existing methods and obtain results that are close to the level of accuracy of FEM results. In order to train an ANN, analyses on a large amount of data were executed with various standardized pipe geometries and burial depths regulated by the Korea Highway Corporation using the two‐dimensional FEM verified in this study. The widely used back propagation algorithm was adopted. The ANN‐based model developed in this study was shown to be an effective tool by comparing the results with test data and sensitivity analyses were executed based on the data from the developed ANN. Copyright © 2010 John Wiley & Sons, Ltd.     

Three‐dimensional nonlinear finite element analysis of pile groups in saturated porous media using a new transmitting boundary

The dynamic behaviour of pile groups subjected to an earthquake base shaking is analysed. An analysis is formulated in the time domain and the effects of material nonlinearity of soil, pile–soil–pile kinematic interaction and the superstructure–foundation inertial interaction on seismic response are investigated. Prediction of response of pile group–soil system during a large earthquake requires consideration of various aspects such as the nonlinear and elasto‐plastic behaviour of soil, pore water pressure generation in soil, radiation of energy away from the pile, etc. A fully explicit dynamic finite element scheme is developed for saturated porous media, based on the extension of the original formulation by Biot having solid displacement (u) and relative fluid displacement (w) as primary variables (u–w formulation). All linear relative fluid acceleration terms are included in this formulation. A new three‐dimensional transmitting boundary that was developed in cartesian co‐ordinate system for dynamic response analysis of fluid‐saturated porous media is implemented to avoid wave reflections towards the structure. In contrast to traditional methods, this boundary is able to absorb surface waves as well as body waves. The pile–soil interaction problem is analysed and it is shown that the results from the fully coupled procedure, using the advanced transmitting boundary, compare reasonably well with centrifuge data. Copyright © 2007 John Wiley & Sons, Ltd.     

Implementation of the finite element method in the three‐dimensional discontinuous deformation analysis (3D‐DDA)

A modified three‐dimensional discontinuous deformation analysis (3D‐DDA) method is derived using four‐noded tetrahedral elements to improve the accuracy of current 3D‐DDA algorithm in practical applications. The analysis program for the modified 3D‐DDA method is developed in a C++ environment and its accuracy is illustrated through comparisons with several analytical solutions that are available for selected problems. The predicted solutions for these problems using the modified 3D‐DDA approach all show satisfactory agreement with the corresponding analytical results. Results presented in this paper demonstrate that the modified 3D‐DDA method with discontinuous modeling capabilities offers a useful computational tool to determine stresses and deformations in practical problems involving fissured elastic media with reasonable accuracy. Copyright © 2008 John Wiley & Sons, Ltd.     

Three‐dimensional finite element study of arching behavior in slope/drilled shafts system

Two‐dimensional slope stability analysis for a slope with a row of drilled shafts needs a mechanism to take into account the three‐dimensional effect of the soil arching due to the spaced drilled shafts on slope. To gain a better understanding of the arching mechanisms in a slope with evenly spaced drilled shafts socketed into a stable stratum (or a rock layer), the three‐dimensional finite element modelling technique was used for a comprehensive parametric study, where the nonlinear and plastic nature of the soil and the elastic behavior of the drilled shafts as well as the interface frictions were modelled. Various factors were varied in the parametric study to include (1) the rigidity of the drilled shafts as influenced by its diameter, modulus of elasticity, and total length; (2) shafts spacing and location on the slope; (3) the material properties of rock and the socket length of shaft; and (4) the soil movement and strength parameters. Evidences of soil arching and reduction in the stresses and displacements through the load transfer mechanisms due to the presence of the drilled shafts were elucidated through the finite element method (FEM) simulation results. Design charts based on regression analysis of FEM simulation results were created to obtain a numerical value of the load transfer factor for the arching mechanism provided by the drilled shafts on the slope. Observations of the arching behavior learned from the FEM simulations provide an insight into the behavior of drilled shafts stabilized slope. Copyright © 2009 John Wiley & Sons, Ltd.     

群桩与筏板的竖向动力相互作用研究

对沿深度离散形成一系列均质的土层,应用传递矩阵法推导出单桩的动力阻抗。借助推求的桩-桩动力相互作用因子,研究承台板与群桩的竖向动力相互作用,通常假设承台板是刚性的,并充分考虑了承台板具有一定的刚度、可以变形的柔性承台板,提出承台板-群桩动力相互作用模型,然后推出其动力相互作用方程。通过算例说明了承台板的刚度对承台板的动力作用是有影响的。研究结果表明,振动频率越低,柔性承台板的振幅与刚性承台板的振幅相差越大;随着频率的增大,二者差别也减小,并逐渐趋于相等。     

Three‐dimensional finite elements for the analysis of soil contamination using a multiple‐porosity approach

A three‐dimensional finite‐element model of contaminant migration in fissured clays or contaminated sand which includes multiple sources of non‐equilibrium processes is proposed. The conceptual framework can accommodate a regular network of fissures in 1D, 2D or 3D and immobile solutions in the macro‐pores of aggregated topsoils, as well as non‐equilibrium sorption. A Galerkin weighted‐residual statement for the three‐dimensional form of the equations in the Laplace domain is formulated. Equations are discretized using linear and quadratic prism elements. The system of algebraic equations is solved in the Laplace domain and solution is inverted to the time domain numerically. The model is validated and its scope is illustrated through the analysis of three problems: a waste repository deeply buried in fissured clay, a storage tank leaking into sand and a sanitary landfill leaching into fissured clay over a sand aquifer. Copyright © 2005 John Wiley & Sons, Ltd.     

Two‐dimensional finite element analysis of gravitational and lateral‐driven deformation in sedimentary basins

The paper deals with the modeling of some aspects, such as the formulation of constitutive equations for sediment material or finite element approach for basin analysis, related to mechanical compaction in sedimentary basins. In addition to compaction due to gravity forces and pore‐pressure dissipation, particular emphasis is given to the study of deformation induced by tectonic sequences. The numerical model relies upon the implementation of a comprehensive constitutive model for the sediment material formulated within the framework of finite poroplasticity. The theoretical model accounts for both hydromechanical and elasticity–plasticity coupling due to the effects of irreversible large strains. From the numerical viewpoint, a finite element procedure specifically devised for dealing with sedimentary basins as open systems allows to simulate within a two‐dimensional setting the process of sediment accretion or erosion. Several basin simulations are presented. The main objective is to analyze the behavior of a sedimentary basin during the different phases of its life cycle: accretion phase, pore‐pressure dissipation phase and compressive/extensional tectonic motions. Copyright © 2013 John Wiley & Sons, Ltd.     

水平荷载下导管架平台桩基础的非线性有限元分析

导管架平台桩基础的控制荷载主要为风荷载、波浪荷载、地震荷载等水平荷载,为研究水平荷载下导管架平台桩基础的承载特性,采用非线性有限元分析方法对水平荷载下桩-土之间的相互作用进行研究,提出了有效模拟桩基水平承载特性的有限元模型,分析了模型桩的刚度、直径、土质参数中水平土压力系数、剪胀角对桩基承载特性的影响及水平荷载下群桩承载特性,并将有限元计算结果与API规范及模型试验结果进行对比。研究结果表明,非线性有限元分析方法分析水平荷载下桩-土相互作用是可行的,计算结果可为导管架平台的桩基设计提供参考。     

浅基础沉降可靠性的摄动随机有限元法

将二阶非线性摄动随机有限元理论和土性参数随机场的离散技术相结合,可以处理岩土工程中由多种不确定因素带来的结果的不确定性。与一阶摄动随机有限元相比,二阶摄动随机有限元具有较高的计算精度。笔者探讨了这一方法,考虑土性分布参数和荷载变异性,用二阶摄动随机有限元法进行了浅地基沉降的概率分析。     

Cauchy problem of three‐dimensional critical slip surfaces of slopes

Using the continuum‐based finite element approach in conjunction with the strength reduction technique, one can easily find out the factor of safety of a slope with rather high precision, but will encounter some obstacles in the accurate location of the corresponding three‐dimensional critical slip surface (CSS). On the basis of the Mohr–Coulomb failure criterion and the stress field in the limit equilibrium state of the slope, it is deduced that the three‐dimensional CSS is the solution of the Cauchy problem of a quasi‐linear partial differential equation (PDE) of first order. By means of the method of characteristics for the problem, the three‐dimensional CSS that may take any shape can be determined without the specification of its geometry. Copyright © 2010 John Wiley & Sons, Ltd.     

地道桥结构与土相互作用的有限元分析

借助有限元软件ANSYS来分析地道桥结构与土的相互作用,采用三维实体建模方法模拟空间范围内的地道桥结构及其周围土体.分析有无列车荷载等两种不同工况下地道桥结构的位移、应力分布规律以及由其引起的地表沉降和土体的应力变化规律.通过有限元分析可以求出结构和土体的位移、应力峰值及其空间位置,为地道桥结构设计提供可靠的分析依据.     

A modulus‐multiplier approach for non‐linear analysis of laterally loaded pile groups

A modulus‐multiplier approach, which applies a reduction factor to the modulus of single pile p–y curves to account for the group effect, is presented for analysing the response of each individual pile in a laterally loaded pile group with any geometric arrangement based on non‐linear pile–soil–pile interaction. The pile–soil–pile interaction is conducted using a 3D non‐linear finite element approach. The interaction effect between piles under various loading directions is investigated in this paper. Group effects can be neglected at a pile spacing of 9 times the pile diameter for piles along the direction of the lateral load and at a pile spacing of 6 times the pile diameter for piles normal to the direction of loading. The modulus multipliers for a pair of piles are developed as a function of pile spacing for departure angle of 0, 90, and 180sup>/sup> with respect to the loading direction. The procedure proposed for computing the response of any individual pile within a pile group is verified using two well‐documented full‐scale pile load tests. Copyright © 2006 John Wiley & Sons, Ltd.     

现浇薄壁管桩水平受力特性三维弹塑性有限元分析

利用三维弹塑性有限元模型,对水平荷载作用下现浇薄壁管桩(简称PCC桩)受力特性进行数值模拟和分析,并与试验结果进行比较,验证了模型的合理性。利用该模型对桩长、桩径以及桩身相对刚度等影响桩受力特性的主要因素进行分析,结果表明：桩径一定的情况下,桩长对桩的破坏形式有很大的影响,弹性桩与刚性桩临界点的长径比为8左右,且弹性桩的受力主要集中在桩的上部,约为0.3倍桩长;桩径对侧向位移和弯矩都有很大影响,桩径越大,桩侧向位移越小,而最大弯矩值越大,且最大弯矩点下移;在一定程度上,提高桩体强度,有利于提高桩的抗水平变形能力,但当增加到一定程度时,效果并不十分明显。     

Homogenization framework for three‐dimensional elastoplastic finite element analysis of a grouted pipe‐roofing reinforcement method for tunnelling

This paper deals with the grouted pipe‐roofing reinforcement method that is used in the construction of tunnels through weak grounds. This system consists on installing, prior to the excavation of a length of tunnel, an array of pipes forming a kind of ‘umbrella’ above the area to be excavated. In some cases, these pipes are later used to inject grout to strengthen the ground and ‘connect’ the pipes. This system has proven to be very efficient in reducing tunnel convergence and water inflow when tunnelling through weak grounds. However, due to the geometrical and mechanical complexity of the problem, existing finite element frameworks are inappropriate to simulate tunnelling using this method. In this paper, a mathematical framework based on a homogenization technique to simulate ‘grouted pipe‐roofing reinforced ground’ and its implementation into a 3‐D finite element programme that can consider stage construction situations are presented. The constitutive model developed allows considering the main design parameters of the problem and only requires geometrical and mechanical properties of the constituents. Additionally, the use of a homogenization approach implies that the generation of the finite element mesh can be easily produced and that re‐meshing is not required as basic geometrical parameters such as the orientation of the pipes are changed. The model developed is used to simulate tunnelling with the grouted pipe‐roofing reinforcement method. From the analyses, the effects of the main design parameters on the elastic and the elastoplastic analyses are considered. Copyright © 2004 John Wiley & Sons, Ltd.     

Discrete element method simulations of the seismic response of shallow foundations including soil‐foundation‐structure interaction

A novel three‐dimensional particle‐based technique utilizing the discrete element method is proposed to analyze the seismic response of soil‐foundation‐structure systems. The proposed approach is employed to investigate the response of a single‐degree‐of‐freedom structure on a square spread footing founded on a dry granular deposit. The soil is idealized as a collection of spherical particles using discrete element method. The spread footing is modeled as a rigid block composed of clumped particles, and its motion is described by the resultant forces and moments acting upon it. The structure is modeled as a column made of particles that are either clumped to idealize a rigid structure or bonded to simulate a flexible structure of prescribed stiffness. Analysis is done in a fully coupled scheme in time domain while taking into account the effects of soil nonlinear behavior, the possible separation between foundation base and soil caused by rocking, the possible sliding of the footing, and the dynamic soil‐foundation interaction as well as the dynamic characteristics of the superstructure. High fidelity computational simulations comprising about half a million particles were conducted to examine the ability of the proposed technique to model the response of soil‐foundation‐structure systems. The computational approach is able to capture essential dynamic response patterns. The cyclic moment–rotation relationships at the base center point of the footing showed degradation of rotational stiffness by increasing the level of strain. Permanent deformations under the foundation continued to accumulate with the increase in number of loading cycles. Copyright © 2011 John Wiley & Sons, Ltd.     

某特深基坑支护的非线性三维有限元分析

讨论了润扬长江大桥北锚碇基础施工开挖中50 m特深基坑的三维非线性有限元分析,介绍了工程的有关情况,讨论了其计算模型和参数取值,给出主要计算结果并对有关影响因素进行了分析,在此基础上对该基坑支护的设计作出评价,最后结合施工结果讨论了计算分析的可靠性。