Simulation of arch dam – foundation interaction with a new friction interface element

A new model for three-dimensional non-linear contact problems with irreversible friction is presented here for the interaction between the rock foundation and an arch dam structure. Based on the finite element method and load incremental theory, a constraint contact element with displacements and contact stresses as node parameters is developed. In this approach, four contact conditions are considered, i.e. fixed, slip, free and mixed. This model can simulate frictional slippage, decoupling and re-bonding of two bodies initially mating at a common interface or with any initial gaps. Furthermore boundary conditions for this element are discussed and treatment measures proposed. This method is shown to be effective and to have the advantage of being easily implemented into standard finite element programs. Solutions are obtained for a centrally loaded, simply supported composite beam and for an end-loaded elastica with initial gaps in regional contact with a rigid surface. The results obtained for these examples are compared to the plane stress solutions by contact friction analysis. As an application example, Quanshui arch dam located in Ruyuan County of Guangdong Province in southern China is simulated with the new element.     

Finite element analyses of soil plug response

Open-ended pipe piles are often used in offshore foundations. The response of the soil plug inside a pipe pile is poorly understood, and only limited work has been performed to quantify the response under the different loading conditions relevant to offshore platforms. This paper describes numerical analyses that have been carried out in order to assess the end-bearing capacity of the soil plug under loading conditions which range from undrained to fully drained. The soil plug has been modelled as either elastic, elastic–perfectly-plastic or elastoplastic. The soil–pile interface, an important aspect of the problem, has been examined critically. Comparison with experimental data from model test at laboratory scale indicates that the load–deformation behaviour of the soil plug is modelled well using an elastoplastic model for the soil plug, and an elastic–perfectly-plastic joint element to model the soil–pile interface. The finite element analyses show that, under typical loading conditions, adequate end bearing may be mobilized by the soil plug, largely by high effective stresses in the bottom 3–5 diameters of the soil plug.     

A multidirectional semi‐analytical method for analysis of laterally loaded pile groups in multi‐layered elastic strata

A semi‐analytical method for calculating the response of single piles and pile groups subjected to lateral loading is developed in this paper. Displacements anywhere in the soil domain are tied to the displacements of the piles through decay functions. The principle of virtual work and the calculus of variations are used to derive the governing differential equations that describe the response of the piles and soil. The eigenvalue method and the finite difference technique are used to solve the system of coupled differential equations for the piles and soil, respectively. The proposed method takes into account the soil surface displacement along and perpendicular to the loading direction and produces displacement fields that are very close to those produced by the finite element method but at lower computational effort. Compared with the previous method that considered only the soil displacement along the loading direction, accounting for the multi‐directional soil displacement field produces responses for the piles and soil that are closer to those approximated by the finite element method. Copyright © 2016 John Wiley & Sons, Ltd.     

An embedded bond-slip model for finite element modelling of soil–nail interaction

Soil nailing has been widely used as a reinforcing technique to retain excavations and stabilise slopes. Proper assessment of the interaction between the nails and the surrounding soil is central to safe and economical design of the composite reinforced soil structure. In this note, a new interface model, denoted as “embedded bond-slip model”, is proposed to model the soil–nail interaction numerically in a simplified manner. Combining the key features of the embedded element technique and the conventional interface element method, the proposed plane–strain interface model has the advantages that no special considerations have to be given to the arrangement of the finite element mesh for the soil nails, and that possible tangential slippage along the interface can be modelled. The formulation also allows pore water flow across the soil nails to be incorporated into the analysis. The proposed model has been implemented into a finite element code and verified by simple element tests under different uni-direction loading conditions. Using the proposed interface model, back analyses of a field test involving a soil-nailed cut slope subjected to a rise in groundwater table have been conducted. This note presents the details of the embedded bond-slip model and the numerical results which demonstrate that the proposed model is capable of simulating soil–nail interaction conveniently and realistically.     

ELASTIC ANALYSIS OF BURIED STRUCTURES SUBJECT TO THREE-DIMENSIONAL SURFACE LOADING

A method is presented which may be used to compute the displacements, strains and moments (both in-plane and transverse) in buried structures such as pipelines and culverts subjected to longitudinal bending. This type of bending can occur if a surface loading such as a vehicular loading or an embankment loading is applied to the soil above the pipe or culvert. Fourier transforms are used to reduce the three-dimensional problem to one involving only two spatial directions, thereby reducing the data preparation and computation time. Conventional finite element analysis is used to approximate the field quantities in the transformed two-dimensional plane. Two Fourier integral element types have been developed which have many applications in geotechnical engineering.     

Modelling of fibre–cohesive soil mixtures

A new constitutive model for fibre-reinforced cohesive soil is proposed. The model combines a Cam-Clay like bounding surface model with an elastic–plastic one-dimensional fibrous element model. A “smearing procedure”, which can consider any spatial distribution of fibre orientation, is employed to transform discrete tensile forces developed in the fibres into stresses for the composite material. The fibre stress contribution is bounded by both degradation of soil–fibre bonding due to pull-out mechanism and tensile strength of the fibres. Eventual occurrence of fibre breakage is also considered. The model performances are analysed for both consolidation and shearing loading modes, and qualitative comparison is performed with experimental data available in the literature. For consolidation loading, tensile stresses are not developed in the fibres and thus the fibre effect is rather limited. For drained shear loading, addition of fibres can result in a consistent shear strength increase. The beneficial effect of fibres seems to be controlled by two parameters: the fibre tensile stiffness and the fibre/soil strain ratio that accounts for any possible slippage or shear deformation at the fibre/soil matrix interface. For undrained shear loading, the strengthening effect of the fibres appears to be counteracted by the increase in pore water pressure, induced by the additional confining contribution of the fibres. In agreement with published experimental data, the model suggests also that the moisture content is a key factor governing fibre effectiveness for undrained shearing. Finally, analysis of the model predicted critical states for fibre-reinforced cohesive soil is provided.     

Prediction of elastic settlement of rectangular raft foundation — a coupled FE–BE approach

A numerical method of analysis is proposed for computation of the elastic settlement of raft foundations using a FEM–BEM coupling technique. The structural model adopted for the raft is based on an isoparametric plate bending finite element and the raft–soil interface is idealized by boundary elements. Mindlin's half-space solution is used as a fundamental solution to find the soil flexibility matrix and consequently the soil stiffness matrix. Transformation of boundary element matrices are carried out to make it compatible for coupling with plate stiffness matrix obtained from the finite element method. This method is very efficient and attractive in the sense that it can be used for rafts of any geometry in terms of thickness as well as shape and loading. Depth of embedment of the raft can also be considered in the analysis. Copyright © 1999 John Wiley & Sons, Ltd.     

青藏铁路多年冻土区桥梁桩基础地震响应的试验研究与数值分析

基于对负温条件下青藏铁路高温不稳定多年冻土区桥梁桩基础的缩尺模型振动台试验,明确了动荷载作用下模型桩-土界面和两桩中间土体存在温度升高响应。在此基础上,考虑天然状态和地震荷载作用下地温升高两种条件,运用动力有限元方法,对青藏铁路清水河特大桥桩基础进行了地震响应分析。计算表明,在高温状态下(-1℃以上),多年冻土的地基及桩基础在地震荷载作用下的动态响应对温度的升高异常敏感;在50年超载概率2%的青藏人工波作用下,因桩基础的相对位移增大了地基与桩基础间的滑移、脱离现象,特别是高温状态下地基与桩基础间的变形出现了明显的不稳定滑移,影响了整个桩基础的稳定性。     

Particle‐scale effects on global axial and torsional interface shear behavior

An extensive literature on the shear behavior of continuum–particulate interfaces has been developed during the last four decades. However, relatively limited work regarding the behavior of interfaces under different loading conditions has been published. This paper presents a discrete element modeling study, along with comparisons from experimental data, of interface behavior under axial and torsional drained loading conditions. Detailed studies allow for links between micro‐scale particle behavior and observed global response to be developed and for the latter to be evaluated in light of particle–particle and particle–continuum interactions. The results of this study indicate that axial and torsional interface shear induce inherently different loading conditions, as shown by the different failure envelopes, stress paths, and induced soil volume changes and deformations. Furthermore, the results presented in this paper indicate that particle‐level mechanisms, such as particle rotations and contact slippage, play different roles in axial and torsional shear. Coordination number, polar histograms, particle displacements, particle rotations, and local void ratio measurements provide further insights into the fabric evolution, loading conditions, and failure mechanisms induced by these two shear modes. This study expands the current understanding of interface behavior and discusses potential improvements to geotechnical systems that leverage the characteristics of different imposed loading conditions. Copyright © 2016 John Wiley & Sons, Ltd.     

现浇混凝土薄壁管桩的荷载传递机理

振动沉模现浇混凝土薄壁管桩（简称PCC桩）技术是河海大学自主开发研制的用于地基加固处理的新技术,它是一种适合于软土地区的新型高效优质桩型。通过有限元方法研究了PCC桩的荷载传递机理,考虑了桩土接触面的滑移、脱开和土的非线性。分析结果表明,荷载-沉降曲线数值模拟结果与现场试验相当接近,并揭示了内外侧摩阻力的分布、端阻力的发挥、土塞的作用等一些目前无法由现场试验得到的结果。     

局部加载条件下二维土质边坡稳定性的数值分析

借助Hill的材料应变局部化原理和土的弹塑性本构模型,应用声学张量示踪边坡失稳滑移线发生和发展,提出了局部加载条件下的确定边坡失稳滑移线的方法。在弹塑性有限元分析过程中,应用边坡材料的塑性一致性条件或滑移线两侧的应力连续条件,即应变局部化条件,检测每一个高斯点的声学张量行列式的值,当其等于或接近0时,标出滑移萌生的位置、方向,模拟边坡滑移线萌生和发展过程。     

Coupled theory of mixtures for clayey soils

In this work, elasto-plastic coupled equations are formulated in order to describe the time-dependent deformation of saturated cohesive soils (two-phase state). Formulation of these equations is based on the principle of virtual work and the theory of mixtures for inelastic porous media. The theory of mixtures for a linear elastic porous skeleton was first developed by Biot (Theory of elasticity and consolidation for a porous anisotropic solid, Journal of Applied Physics, 1955, 26, 188–185). An extension of Biot's theory into a nonlinear inelastic media was performed by Prevost (Mechanics of continuous porous media, International Journal of Engineering Science, 1980, 18, 787–800). The saturated soil is considered as a mixture of two deformable media, the solid grains and the water. Each medium is regarded as a continuum and follows its own motion. The flow of pore-water through the voids is assumed to follow Darcy's law. The coupled equations are developed for large deformations with finite strains in an updated Lagrangian reference frame. The coupled behavior of the two-phase materials (soil-water state) is implemented in a finite element program. A modified Cam-clay model is adopted and implemented in the finite element program in order to describe the plastic behavior of clayey soils. Penetration of a piezocone penetrometer in soil is numerically simulated and implemented into a finite element program. The piezocone penetrometer is assumed to be infinitely stiff. The continuous penetration of the cone is simulated by applying an incremental vertical movement of the cone tip boundary. Results of the finite element numerical simulation are compared with experimental measurements conducted at Louisiana State University using the calibration chamber. The numerical simulation is carried out for two cases. In the first case, the interface friction between the soil and the piezocone penetrometer is neglected. In the second case, interface friction is assumed between the soil and the piezocone. The results of the numerical simulations are compared with experimental laboratory measurements.     

Nonlinear three-dimensional analysis of reinforced concrete piles subjected to horizontal loading

A three-dimensional finite element approach is proposed to predict the response of reinforced concrete piles to horizontal loading. This approach allows the nonlinear effects arising from the soil–pile interaction to be properly accounted for. In particular, the occurrence of plastic strains in the soil, concrete cracking and steel yielding in the pile as well as the occurrence of slip and gap at the soil–pile interface are reliably simulated using appropriate constitutive models. Another advantage of the present method is that it requires few material parameters as input data. In addition, these parameters can be readily obtained from conventional geotechnical and structural tests. The proposed approach is used to analyse the results from some loading tests documented in the literature concerning a large-diameter pile and a large-section rectangular pile (barrette) embedded in sandy soils. The theoretical results from these analyses are found to be in fairly good agreement with the experimental measurements available from the loading tests. Remarks of practical interest on the response of the structures considered to horizontal loading are also made.     

Analytical solutions of linear finite‐ and small‐strain one‐dimensional consolidation

Analytical solutions are presented for linear finite‐strain one‐dimensional consolidation of initially unconsolidated soil layers with surcharge loading for both one‐ and two‐way drainage. These solutions complement earlier solutions for initially unconsolidated soil layers without surcharge and initially normally consolidated soil layers with surcharge. Small‐strain solutions for the consolidation of initially unconsolidated soil layers with surcharge loading are also presented, and the relationship between the earlier solutions for initially unconsolidated soil without surcharge and the corresponding small‐strain solutions, which was not addressed in the earlier work, is clarified. The new solutions for initially unconsolidated soil with surcharge loading can be applied to the analysis of low stress consolidation tests and to the partial validation of numerical solutions of non‐linear finite‐strain consolidation. They also clarify a formerly perplexing aspect of finite‐strain solution charts first noted in numerical solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

路基荷载下PCC刚性桩复合地基沉降简化计算

将路堤荷载作用下的刚性桩复合地基变形,简化为只有竖向变形,水平变形可以忽略,然后,应用土体的平衡、几何、物理方程,建立沉降计算的控制方程。桩土界面的接触条件按位移协调或产生滑移变形两种情况分别建立,采用有限差分法离散求解,可得桩、桩周土各自的沉降变形、桩土界面的摩擦力分布。与有限元分析和现场实测结果的对比表明,该简化方法具有较高的准确性,可以方便地应用到工程实际。     

Axial and lateral response of pile groups embedded in nonhomogeneous soils

A numerical method of analysis based on elasticity theory is presented for the analysis of axially and laterally loaded pile groups embedded in nonhomogeneous soils. The problem is decomposed into two systems, namely the group piles acted upon by external applied loads and pile–soil interaction forces, and a layered soil continuum acted upon by a system of pile–soil interaction forces at the imaginary positions of the piles. The group piles are discretized into discrete elements while the nonhomogeneous soil behaviour is determined from an economically viable finite element procedure. The load–deformation relationship of the pile group system is then determined by considering the equilibrium of the pile–soil interaction forces, and the compatibility of the pile and soil displacements. The influence of soil nonlinearity can be studied by limiting the soil forces at the pile–soil interface, and redistributing the ‘excess forces’ by an ‘initial stress’ process popular in elasto-plastic finite element analysis. The solutions from this approach are compared with some available published solutions for single piles and pile groups in homogeneous and nonhomogeneous soils. A limited number of field tests on pile groups are studied, and show that, in general, the computed response compares favourably with the field measurements.     

Negative skin friction on pile groups

A numerical procedure is presented for the downdrag analysis of group piles which penetrate a consolidating upper soil layer to socket into a firm bearing stratum of finite stiffness. The settlement of the consolidating upper soil layer under a surcharge load is estimated using Terzaghi's one-dimensional consolidation theory. Parametric solutions are presented to show the influence of various parameters on the performance of the socketed pile groups in terms of the development of the induced downdrag forces and associated pile head settlements. In general, pile–soil–pile interaction has the beneficial effect of reducing the downdrag forces and settlements of the group piles when compared to the corresponding single pile values, provided that the soil settlements are not so large as to cause full slippage at the interface in all the piles. Reasonable agreement is obtained between the theoretical and experimental results for pile groups subjected to negative skin friction.     

Axisymmetric finite element analysis of pile loading tests

In this paper a typical soil–structure interaction problem is considered, the case of a vertical pile installed in sand and submitted to an axial compression loading. Results from two full scale pile tests are analysed and the tests are reproduced by numerical simulations via finite elements method (FEM). The choice of the mechanical parameters for the soil and the sand–pile interface and the modelling approach are first described. A new numerical strategy is outlined to account for pile installation effects due to jacking and driving via FEM. The proposed approach is based on the application of existing empirical correlations available for the quantification of residual radial and shear stresses along the pile shaft as well as residual pressures around the pile base after the installation. This approach is proposed as an alternative to more complex methods based on the numerical modelling of the pile penetration problem. The role of the constitutive modelling of the interface is also discussed. Finally, comparative analyses of pile loading tests using FEM are provided and the comparisons between numerical and experimental results are presented and discussed.     

Effects of frictional forces acting on sidewalls of buried box culverts

The effects of frictional forces acting on the sidewalls of buried box culverts are presented as determined with finite element method (FEM) and detailed soil modelling. The possibility of reducing earth pressure on deeply buried concrete box culverts by the imperfect trench installation (ITI) method has been contemplated during the last several decades. There have been limited research results published primarily regarding the qualitative aspect of load reduction in ITIs. It was found during the course of this study that significant frictional forces develop along the sidewalls of box culverts and adjacent sidefills in ITIs. Current American Association of State Highway and Transportation Officials provisions do not consider these frictional forces, but they cannot be neglected in ITIs, as their effect is dominant. An optimum geometry for the soft zone in ITIs is presented to maximize earth load reductions. The soil–structure interaction at the box culvert–soil interface was found to have a significant effect on total earth pressure acting on the bottom slab. Predictor equations for earth load reduction rates were formulated for ITIs incorporating the optimum soft zone geometry based on the FEM. Copyright © 2007 John Wiley & Sons, Ltd.