Analysis of undrained cavity expansion in elasto‐plastic soils with non‐linear elasticity

A large strain analysis of undrained expansion of a spherical/cylindrical cavity in a soil modelled as non‐linear elastic modified Cam clay material is presented. The stress–strain response of the soil is assumed to obey non‐linear elasticity until yielding. A power‐law characteristic or a hyperbolic stress–strain curve is used to describe the gradual reduction of soil stiffness with shear strain. It is assumed that, after yielding, the elasto‐plastic behaviour of the soil can be described by the modified Cam clay model. Based on a closed‐form stress–strain response in undrained condition, a numerical solution is obtained with the aid of simple numerical integration technique. The results show that the stresses and the pore pressure in the soil around an expanded cavity are significantly affected by the non‐linear elasticity, especially if the soil is overconsolidated. The difference between large strain and small strain solutions in the elastic zone is not significant. The stresses and the pore pressure at the cavity wall can be expressed as an approximate closed‐form solution. Copyright © 2001 John Wiley & Sons, Ltd.     

重复荷载作用下粉性路基土累积塑性变形研究

为研究不同含水率、压实度和应力水平下路基压实粉性路基土累积塑性变形影响因素与变化规律,开展了系列动三轴试验。根据累积塑性应变随加载次数的变化规律,获得了不同含水率、压实度下粉土的临界动应力。研究表明,粉性路基土临界动应力随压实度的提高而增大,随含水率的增大而减小。为避免路面结构的破坏,路基土应力状态应控制在临界应力界限范围内。鉴于此,针对路基土临界应力范围内的试验数据,建立了路基土累积塑性变形预估模型。该模型考虑了应力和加载次数的影响。回归分析表明,该模型具有较高的决定系数,即模型具有较高的合理性与可靠性。模型的建立为基于力学经验法的路基累积塑性变形计算提供了参数。最后以典型水泥混凝土路面为例,获得了不同轴载、轴型作用下粉土路基累积塑性变形,为进一步研究路基累积塑性变形对路面结构的影响提供了思路。     

Response evaluation for horizontally loaded fixed‐head pile groups using 3‐D non‐linear analysis

The response of laterally loaded pile foundations may be significantly important in the design of structures for such loads. A static horizontal pile load test is able to provide a load–deflection curve for a single free‐head pile, which significantly differs from that of a free‐ or fixed‐head pile group, depending on the particular group configuration. The aim of this paper is to evaluate the influence of the interaction between the piles of a group fixed in a rigid pile cap on both the lateral load capacity and the stiffness of the group. For this purpose, a parametric three‐dimensional non‐linear numerical analysis was carried out for different arrangements of pile groups. The response of the pile groups is compared to that of the single pile. The influence of the number of piles, the spacing and the deflection level to the group response is discussed. Furthermore, the contribution of the piles constituting the group to the total group resistance is examined. Finally, a relationship is proposed allowing a reasonable prediction of the response of fixed‐head pile groups at least for similar soil profile conditions. Copyright © 2005 John Wiley & Sons, Ltd.     

A bounding surface elasto‐viscoplastic constitutive model for non‐isothermal cyclic analysis of asphaltic materials

An elasto‐viscoplastic constitutive model for asphaltic materials is presented within the context of bounding surface plasticity theory, taking into account the effects of the stress state, void binder degree of saturation, temperature and strain rate on the material behaviour. A stress state dependent non‐linear elasticity model is introduced to represent time‐independent recoverable portion of the deformation. The consistent visco‐plasticity framework is utilised to capture the rate‐dependent, non‐recoverable strain components. The material parameters introduced in the model are identified, and their determination from conventional laboratory tests is discussed. The capability of the model to reproduce experimentally observed response of asphaltic materials is demonstrated through numerical simulations of several laboratory test data from the literature. Copyright © 2016 John Wiley & Sons, Ltd.     

On finite volume method implementation of poro‐elasto‐plasticity soil model

Accurate prediction of the interactions between the nonlinear soil skeleton and the pore fluid under loading plays a vital role in many geotechnical applications. It is therefore important to develop a numerical method that can effectively capture this nonlinear soil‐pore fluid coupling effect. This paper presents the implementation of a new finite volume method code of poro‐elasto‐plasticity soil model. The model is formulated on the basis of Biot's consolidation theory and combined with a perfect plasticity Mohr‐Coulomb constitutive relation. The governing equation system is discretized in a segregated manner, namely, those conventional linear and uncoupled terms are treated implicitly, while those nonlinear and coupled terms are treated explicitly by using any available values from previous time or iteration step. The implicit–explicit discretization leads to a linearized and decoupled algebraic system, which is solved using the fixed‐point iteration method. Upon the convergence of the iterative method, fully nonlinear coupled solutions are obtained. Also explored in this paper is the special way of treating traction boundary in finite volume method compared with FEM. Finally, three numerical test cases are simulated to verify the implementation procedure. It is shown in the simulation results that the implemented solver is capable of and efficient at predicting reasonable soil responses with pore pressure coupling under different loading situations. Copyright © 2015 John Wiley & Sons, Ltd.     

公路路基地基承载力的离心模型试验与分析

公路路基由于基底的柔性和荷载分布为梯形,路基地基的破坏模式和确定方法与刚性基础地基有所差异,对公路路基地基承载力的研究具有重要的理论意义和实用价值。通过柔性条形基础、梯形路基以及模拟沟谷地形的柔性条形基础等3组离心模型试验,对公路路基地基的破坏机制和极限承载力进行了研究,并采用关联流动的Mohr-Coulomb内切圆屈服准则对试验结果进行了数值模拟和对比分析,研究结果表明,公路路基地基由于基底的柔性和较高的离心加速度,很难产生理想的整体剪切破坏形式;与刚柔性条形基础相比,梯形路基作用下的基底中心点沉降更大;沟谷地形对柔性条形基础基底中心点的沉降有一定的影响,但并不明显影响地基极限承载力。     

Discrete element method analysis of non‐coaxial flow under rotational shear

This study focuses on non‐coaxial flow behavior of cohesionless soil undergoing cyclic rotational shear, with a special interest in the effects of particle‐scale characteristics. To this end, we perform a series of 2D discrete element simulations with various particle shapes, inter‐particle coefficient of friction, initial density, and stress ratios. The validity and efficacy of the numerical model is established by systematically comparing numerical simulation results with existing laboratory testing results. Such comparison shows that the numerical simulations are capable of capturing mechanical behavior observed in laboratory testing under rotational shear. We further demonstrate and quantify a strong yet simple relationship between the deviatoric part of the normalized strain increment and the non‐coaxial angle, denoted by and ψ, respectively. This quantitative correlation between ψ and is independent of applied stress ratio, initial and current void ratio, and the number of cycles applied, but dependent on the principal stress orientation and particle‐scale characteristics. At the same , specimens with higher inter‐particle friction angle or smaller particle aspect ratio show greater non‐coaxial angles. A simple model is able to fit this ψ‐ relationship well, which provides a useful relationship that can be exploited in developing constitutive models for rotational shearing. Copyright © 2014 John Wiley & Sons, Ltd.     

Bifurcation and instability modelling by a multimechanism elasto‐plastic model

The bifurcation and instability conditions in geomechanics are closely related to the elasto‐plastic behaviour. In this paper the potential of a multimechanism elasto‐plastic model to predict various modes of failure is examined. First, a brief overview for the essential aspects of the constitutive model and the development of the elasto‐plastic constitutive matrix for this model are presented. Then, numerical simulations of different drained and undrained paths in the axisymmetric and plane‐strain conditions for the Hostun sand are illustrated. These examples confirm the capacity of the model to reproduce instability and strain localization phenomena. The obtained response is in agreement with experimental observations, theoretical developments and numerical analyses existing in the literature. Copyright © 2007 John Wiley & Sons, Ltd.     

Scaled boundary finite‐element analysis of a non‐homogeneous axisymmetric domain subjected to general loading

The scaled boundary finite‐element method is derived for elastostatic problems involving an axisymmetric domain subjected to a general load, using a Fourier series to model the variation of displacement in the circumferential direction of the cylindrical co‐ordinate system. The method is particularly well suited to modelling unbounded problems, and the formulation allows a power‐law variation of Young's modulus with depth. The efficiency and accuracy of the method is demonstrated through a study showing the convergence of the computed solutions to analytical solutions for the vertical, horizontal, moment and torsion loading of a rigid circular footing on the surface of a homogeneous elastic half‐space. Computed solutions for the vertical and moment loading of a smooth rigid circular footing on a non‐homogeneous half‐space are compared to analytical ones, demonstrating the method's ability to accurately model a variation of Young's modulus with depth. Copyright © 2003 John Wiley & Sons, Ltd.     

Upper bound limit analysis using linear finite elements and non‐linear programming

A new method for computing rigorous upper bounds on the limit loads for one‐, two‐ and three‐dimensional continua is described. The formulation is based on linear finite elements, permits kinematically admissible velocity discontinuities at all interelement boundaries, and furnishes a kinematically admissible velocity field by solving a non‐linear programming problem. In the latter, the objective function corresponds to the dissipated power (which is minimized) and the unknowns are subject to linear equality constraints as well as linear and non‐linear inequality constraints. Provided the yield surface is convex, the optimization problem generated by the upper bound method is also convex and can be solved efficiently by applying a two‐stage, quasi‐Newton scheme to the corresponding Kuhn–Tucker optimality conditions. A key advantage of this strategy is that its iteration count is largely independent of the mesh size. Since the formulation permits non‐linear constraints on the unknowns, no linearization of the yield surface is necessary and the modelling of three‐dimensional geometries presents no special difficulties. The utility of the proposed upper bound method is illustrated by applying it to a number of two‐ and three‐dimensional boundary value problems. For a variety of two‐dimensional cases, the new scheme is up to two orders of magnitude faster than an equivalent linear programming scheme which uses yield surface linearization. Copyright © 2001 John Wiley & Sons, Ltd.     

Poro‐elasto‐plastic response of an unconsolidated formation confined with stiff seal rocks under radial injection

In this article, we evaluate geomechanics of fluid injection from a fully penetrating vertical well into an unconsolidated formation confined with stiff seal rocks. The coupled behavior of an isotropic, homogeneous sand layer is studied under injection pressures that are high enough to induce plasticity yet not fracturing. Propagation of the significant influence zone surrounding the injection borehole, quantified by the extent of the plastic domain in the elasto‐plastic model, is examined for the first time. First, a new fully coupled axisymmetric numerical model is developed. A comprehensive assessment is performed on pore pressures, stresses/strains, and failure planes during the entire transient period of an injection cycle. Results anticipate existence of five distinctive zones in terms of plasticity state: liquefaction at wellbore; two inner plastic domains surrounding the wellbore, where failure occurs along two planes and major principal stress is in vertical direction; remaining of the plastic domain, where formation fails along one plane and major principal stress is in radial direction; and a non‐plastic region. Failure mechanism at the wellbore is found to be shear followed by liquefaction. Next, a novel methodology is proposed based on which new weakly coupled poro‐elasto‐plastic analytical solutions are derived for all three stress/strain components. Unlike previous studies, extension of the plastic zone is obtained as a function of injection pressure, incorporating plasticity effects on the subsequent elastic domain. Solutions, proven to be a good approximation of numerical simulations, offer a huge advantage as the run time of coupled numerical simulations is considerably long. Copyright © 2016 John Wiley & Sons, Ltd.     

Preconditioned IDR(s) iterative solver for non‐symmetric linear system associated with FEM analysis of shallow foundation

Non‐associated flow rule is essential when the popular Mohr–Coulomb model is used to model nonlinear behavior of soil. The global tangent stiffness matrix in nonlinear finite element analysis becomes non‐symmetric when this non‐associated flow rule is applied. Efficient solution of this large‐scale non‐symmetric linear system is of practical importance. The standard Krylov solver for a non‐symmetric solver is Bi‐CGSTAB. The Induced Dimension Reduction [IDR(s)] solver was proposed in the scientific computing literature relatively recently. Numerical studies of a drained strip footing problem on homogenous soil layer show that IDR(s = 6) is more efficient than Bi‐CGSTAB when the preconditioner is the incomplete factorization with zero fill‐in of global stiffness matrix K_ep (ILU(0)‐K_ep). Iteration time is reduced by 40% by using IDR(s = 6) with ILU(0)‐K_ep. To further reduce computational cost, the global stiffness matrix K_ep is divided into two parts. The first part is the linear elastic stiffness matrix K_e, which is formed only once at the beginning of solution step. The second part is a low‐rank matrix Δ, which is re‐formed at each Newton–Raphson iteration. Numerical studies show that IDR(s = 6) with this ILU(0)‐K_e preconditioner is more time effective than IDR(s = 6) with ILU(0)‐K_ep when the percentage of yielded Gauss points in the mesh is less than 15%. The total computation time is reduced by 60% when all the recommended optimizing methods are used. Copyright © 2013 John Wiley & Sons, Ltd.     

Linear stability of one‐dimensional non‐Darcy flow in broken rocks

Considering the effect of non‐Darcy flow, the perturbation theory and normal mode method are introduced to analyze the linear stability of one‐dimensional non‐Darcy flow of gases in broken rocks. A stability criterion for linear systems is obtained theoretically when the steady states of pressure and velocity fields are perturbed, and the effects of the physical parameters on the linear governing system are studied theoretically and numerically. It is pointed out that the deviation coefficient from Darcy's law plays an important role in the governing system; the increasing absolute value of deviation coefficient from Darcy's law stabilizes the system, and the numerical results are shown graphically. Copyright © 2015 John Wiley & Sons, Ltd.     

General coupling extended multiscale FEM for elasto‐plastic consolidation analysis of heterogeneous saturated porous media

This paper presents a general coupling extended multiscale FEM (GCEMs) for solving the coupling problem of elasto‐plastic consolidation of heterogeneous saturated porous media. In the GCEMs, the numerical multiscale base functions for the solid skeleton and fluid phase of the coupling system are all constructed on the basis of the equivalent stiffness matrix of the unit cell, which not only contain the interaction between the solid and fluid phases but also consider the time effect. Furthermore, in order to improve the computational accuracy for two‐dimensional problems, a multi‐node coarse element strategy for the GCEMs is proposed, and a two‐scale iteration algorithm for the elasto‐plastic consolidation analysis is developed. Some one‐dimensional and two‐dimensional homogeneous and heterogeneous numerical examples are carried out to validate the proposed method through the comparison with the coupling multiscale FEM and standard FEM. Numerical results show that the newly developed GCEMs can almost preserve the same convergent property as the standard FEM and also possesses the advantages of high computational efficiency. In addition, the GCEMs can be easily applied to other coupling multifield and multiphase transient problems. Copyright © 2014 John Wiley & Sons, Ltd.     

A finite strain solution for the elastoplastic ground response curve in tunnelling: rocks with non‐linear failure envelopes

This paper generalizes the finite strain Coulomb solution of Vrakas and Anagnostou (Int J Numer Anal Meth Geomech 2014; 38(11): 1131–1148) for the classic tunnel mechanics problem of the ground response curve to elastoplastic grounds satisfying a non‐linear Mohr's failure criterion. A linear (Coulomb‐type) plastic potential function is used, leading to a non‐associated flow law, and edge plastic flow is considered in the plastic zone. The solution for a general non‐linear Mohr's failure criterion is semi‐analytical in that it requires the evaluation of definite integrals. In the special case of the Hoek–Brown criterion, however, these integrals are calculated analytically, resulting in a rigorous closed‐form series solution. The applicability of the derived solution is illustrated through the example of the Yacambú‐Quibor tunnel, where very large deformations were observed when crossing of weak graphitic phyllites. Copyright © 2016 John Wiley & Sons, Ltd.     

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.     

Displacement‐controlled method and its applications to material non‐linearity

For the analysis of non‐linear problems, the displacement‐controlled method (DCM) has a more extensive application scope and more powerful abilities than the load‐controlled method (LCM). However, difficulties of the DCM's procedure not amenable to most finite element implementations of the conventional LCM have restricted its applications in geomechanics. By means of Sherman–Morrison's theorem, the solution of DCM is improved. The improved procedure is characterized by high efficiency, good numerical stability and a programme structure similar to LCM. Two aspects of applications of DCM are illustrated. The first application is to compute the response of a structure under a given load level like the conventional finite element analysis. The second application is to trace the equilibrium path of a structure under a given load distribution type. A simple but effective algorithm is presented for automatically adjusting the step length in tracing the equilibrium path. Examples illustrate that the proposed procedures are suited for modelling complicated non‐linear problems in geomechanics. Copyright © 2005 John Wiley & Sons, Ltd.     

Non‐linear consolidation analysis of soil with variable compressibility and permeability under cyclic loadings

In this paper, a simple semi‐analytical method has been developed to solve the one‐dimensional non‐linear consolidation problems by considering the changes of compressibility and permeability of the soil layer, subjected to complicated time‐dependent cyclic loadings at the ground surface. The solution presented here takes into account e ～ lg k_v and e ～ lg σ′ linear responses. With c_k the slope of the e ～ lg k_v line and c_c as the slope of the e ～ lg σ′ line, the identified parameter c_c/c_k is found to control the rate of consolidation. Using the solutions obtained, some diagrams are prepared and the relevant behaviours of one‐dimensional non‐linear consolidation of saturated soft soil under cyclic loadings are discussed. The method in this paper does not require any special data; conventional oedometer data can be used. Therefore, the method is particularly efficient and convenient for engineering practice. Copyright © 2006 John Wiley & Sons, Ltd.     

An elasto‐plastic three phase model for partially saturated soil for the finite element simulation of compressed air support in tunnelling

This paper presents a fully coupled finite element formulation for partially saturated soil as a triphasic porous material, which has been developed for the simulation of shield tunnelling with heading face support using compressed air. While for many numerical simulations in geotechnics use of a two‐phase soil model is sufficient, the simulation of compressed air support demands the use of a three‐phase model with the consideration of air as a separate phase. A multiphase model for soft soils is developed, in which the individual constituents of the soil—the soil skeleton, the fluid and the gaseous phase—and their interactions are considered. The triphasic model is formulated within the framework of the theory of porous media, based upon balance equations and constitutive relations for the soil constituents and their mixture. An elasto‐plastic, cam–clay type model is extended to partially saturated soil conditions by incorporating capillary pressure according to the Barcelona basic model. The hydraulic properties of the soil are described via DARCY 's law and the soil–water characteristic curve after VAN GENUCHTEN . Water is modelled as an incompressible and air as a compressible phase. The model is validated by means of selected benchmark problems. The applicability of the model to geotechnical problems is demonstrated by results from the simulation of a compressed air intervention in shield tunnelling. Copyright © 2009 John Wiley & Sons, Ltd.