挡土墙主动土压力塑性临界深度的解析解

墙后塑性区的临界深度问题一直没有得到很好解决,传统计算公式仅适用于一些特殊情况。基于极限平衡理论,视墙后填土为服从Mohr-Coulomb屈服准则的理想弹塑性材料,假定塑性区的一族滑移线为直线即平面滑裂面,提出弹性覆盖层取代传统的张拉裂缝,建立了较为完善的滑楔分析模型,采用极限平衡法推导了在一般情况下的塑性临界压力、临界深度以及塑性区可能最大深度的解析解。计算结果表明,塑性临界深度的解析解与目前文献采用迭代计算的结果完全一致,传统计算公式是该解析解的一个特例。     

Response of a circular opening in a friable low-permeability medium to temperature and pore pressure changes

A general analysis using an incremental elastic, perfectly plastic constitutive stress–strain relationship for poroelastoplastic materials is presented to simulate an opening in a low-permeability friable porous medium under non-isothermal conditions. Analytical solutions are obtained for the stresses and strains around a 2-D plane strain circular borehole. An expansion potential is introduced by combining the strains induced by temperature and pore pressure changes. Steady-state pressures and temperatures are considered, and a non-associated plastic flow rule is applied to calculate plastic strains. Focusing on stress distribution near a circular opening, the classic solutions for those stresses under dry and isothermal conditions are used to compare with the newly derived solution. The general poroelastoplastic effect and the thermal effect on sand production and borehole stability are addressed. We suggest that the knowledge of stress history is critical to achieve adequate solutions for displacement and stress in friable media such as clays, shales and oil sands.     

Large strain similarity solution for a spherical or circular opening excavated in elastic‐perfectly plastic media

This paper deals with the unloading problem of a spherical or circular opening excavated in elastic‐perfectly plastic media with a nonassociated Mohr–Coulomb yield criterion. A large strain similarity solution, using incremental velocity approach, is presented by replacing partial differential equations from stress equilibrium, constitutive law, consistency condition, and displacement equation with first‐order ordinary differential equations. The classical Runge–Kutta method is used to solve the first‐order ordinary differential equations. Comparisons among small and large strain solutions are made using some data sets of soil and rock. The results show that the displacements by large strain similarity solution are smaller than those by exact small strain solution and somewhat larger than those by large strain solution using total strain approach. Copyright © 2014 John Wiley & Sons, Ltd.     

Relationship between small and large strain solutions for general cavity expansion problems in elasto-plastic soils

This paper presents a closed-form relationship between small and finite strain cavity expansion solutions. Its derivation is based on the non-linearly elastic–perfectly plastic cylindrical (or spherical) problem considering a general Mohr’s criterion and constant plastic dilatancy. It is shown, however, that it is sufficiently accurate for general expansion problems not obeying plane-strain rotationally (or spherically) symmetric conditions and involving strain-hardening/softening constitutive behaviour. Therefore, this relationship quantifies the error stemming from the computational assumption of small deformations and provides a simple and efficient way of accounting for geometric non-linearity based entirely on conventional computational methods: ‘self-correction’ of small strain analyses results.     

Elasto‐plastic seepage‐induced stresses due to tunneling

When tunneling is carried out beneath the groundwater table, hydraulic boundary is altered, resulting in seepage entering into the tunnel. The development of flow into the tunnel induces seepage stresses in the ground and the lining is subjected to additional loads. This can often cause fine particles to move, which clog the filter resulting in the long‐term hydraulic deterioration of the drainage system. However, the effect of seepage force is generally not considered in the analysis of tunnel. While several elastic solutions have been proposed by assuming seepage in an elastic medium, stress solutions have not been considered for the seepage force in a porous elasto‐plastic medium. This paper documents a study that investigates the stress behavior, caused by seepage, of a tunnel in an elasto‐plastic ground and its effects on the tunnel and ground. New elasto‐plastic solutions that adopt the Mohr–Coulomb failure criterion are proposed for a circular tunnel under radial flow conditions. A simple solution based on the hydraulic gradient obtained from a numerical parametric study is also proposed for practical use. It should be noted that the simple equation is useful for acquiring additional insight into a problem on a tunnel under drainage, because only a minimal computational effort is needed and considerable economic benefits can be gained by using it in the preliminary stage of tunnel design. The proposed equations were partly validated by numerical analysis, and their applicability is illustrated and discussed using an example problem. Comments on the tunnel analysis are also provided. Copyright © 2010 John Wiley & Sons, Ltd.     

A non‐linear coupled finite element–boundary element model for the prediction of vibrations due to vibratory and impact pile driving

This paper presents a non‐linear coupled finite element–boundary element approach for the prediction of free field vibrations due to vibratory and impact pile driving. Both the non‐linear constitutive behavior of the soil in the vicinity of the pile and the dynamic interaction between the pile and the soil are accounted for. A subdomain approach is used, defining a generalized structure consisting of the pile and a bounded region of soil around the pile, and an unbounded exterior linear soil domain. The soil around the pile may exhibit non‐linear constitutive behavior and is modelled with a time‐domain finite element method. The dynamic stiffness matrix of the exterior unbounded soil domain is calculated using a boundary element formulation in the frequency domain based on a limited number of modes defined on the interface between the generalized structure and the unbounded soil. The soil–structure interaction forces are evaluated as a convolution of the displacement history and the soil flexibility matrices, which are obtained by an inverse Fourier transformation from the frequency to the time domain. This results in a hybrid frequency–time domain formulation of the non‐linear dynamic soil–structure interaction problem, which is solved in the time domain using Newmark's time integration method; the interaction force time history is evaluated using the θ‐scheme in order to obtain stable solutions. The proposed hybrid formulation is validated for linear problems of vibratory and impact pile driving, showing very good agreement with the results obtained with a frequency‐domain solution. Linear predictions, however, overestimate the free field peak particle velocities as observed in reported field experiments during vibratory and impact pile driving at comparable levels of the transferred energy. This is mainly due to energy dissipation related to plastic deformations in the soil around the pile. Ground vibrations due to vibratory and impact pile driving are, therefore, also computed with a non‐linear model where the soil is modelled as an isotropic elastic, perfectly plastic solid, which yields according to the Drucker–Prager failure criterion. This results in lower predicted free field vibrations with respect to linear predictions, which are also in much better agreement with experimental results recorded during vibratory and impact pile driving. Copyright © 2008 John Wiley & Sons, Ltd.     

A new closed-form solution for circular openings modeled by the Unified Strength Theory and radius-dependent Young’s modulus

A new closed-form solution is presented for the stress and displacement distribution surrounding circular openings with finite external radii that are subject to uniform internal and external pressures under plane strain conditions. The specific solution for a deep circular tunnel in an infinite rock mass is also provided. It is assumed that the rock mass is elastic–brittle–plastic and governed by the Unified Strength Theory (UST). In the plastic zone, the radius-dependent Young’s modulus (RDM) model and a non-associated linear flow rule were adopted to establish the radial displacement solution. The new closed-form solution obtained in this paper is a series of results rather than one specific solution; hence, it is suitable for a wide range of rock masses and engineering structures. The traditional solutions, which are based on the Mohr–Coulomb failure criterion and the Generalized Twin Shear Stress yield criterion, can be categorized as special cases of this proposed solution. This new solution agrees reasonably well with the results of a borehole collapse test, a secondary development numerical simulation and an additional closed-form solution using the generalized non-linear Hoek–Brown failure criterion. Parametric studies were conducted to investigate the effects of intermediate principal stress, RDM and dilatancy on the results. It is shown herein that the effects of intermediate principal stress and dilatancy are significant; the RDM model is recommended as the optimum approach for calculating radial displacement and support pressure.     

Analysis of pressuremeter geometry effects in clay using critical state models

The conventional interpretation methods of pressuremeter testing effectively approximate pressuremeter membranes as infinitely long. As a result, the effects of the two‐dimensional geometry of pressuremeters are ignored, leading to an overestimation of soil shear strength by pressuremeter testing, as demonstrated in several previous studies. This paper presents results of a numerical study of two‐dimensional geometry effects on self‐boring pressuremeter tests in undrained clay. The results are obtained using critical state soil models with an effective stress formulation. This is in contrast to most (if not all) existing studies on pressuremeter geometry effects, which were based on perfectly plastic soil models (e.g. Yu (Cavity expansion theory and its application to the analysis of pressuremeters. DPhil Thesis, The University of Oxford, 1990), Yeung and Carter (Proc. 3rd Int. Symp. on Pressuremeters, 1990), and Houlsby and Carter (Géotechnique, 1993; 43 (4):567–576)). The present study suggests that the overestimation of soil strength due to the neglect of finite pressuremeter length is significantly affected by the soil model used in the calculations. It is found that for clays with a high overconsolidation ratio (OCR) the strength overestimation predicted using critical state soil models could be considerably smaller than that predicted using perfectly plastic soil models. The main conclusion of this numerical study is that care must be exercised before directly applying any numerically determined pressuremeter geometry correction factors in practice. Copyright © 2005 John Wiley & Sons, Ltd.     

Ground response curves for rock masses exhibiting strain‐softening behaviour

A literature review has shown that there exist adequate techniques to obtain ground reaction curves for tunnels excavated in elastic‐brittle and perfectly plastic materials. However, for strain‐softening materials it seems that the problem has not been sufficiently analysed. In this paper, a one‐dimensional numerical solution to obtain the ground reaction curve (GRC) for circular tunnels excavated in strain‐softening materials is presented. The problem is formulated in a very general form and leads to a system of ordinary differential equations. By adequately defining a fictitious ‘time’ variable and re‐scaling some variables the problem is converted into an initial value one, which can be solved numerically by a Runge–Kutta–Fehlberg method, which is implemented in MATLAB environment. The method has been developed for various common particular behaviour models including Tresca, Mohr–Coulomb and Hoek–Brown failure criteria, in all cases with non‐associative flow rules and two‐segment piecewise linear functions related to a principal strain‐dependent plastic parameter to model the transition between peak and residual failure criteria. Some particular examples for the different failure criteria have been run, which agree well with closed‐form solutions—if existing—or with FDM‐based code results. Parametric studies and specific charts are created to highlight the influence of different parameters. The proposed methodology intends to be a wider and general numerical basis where standard and newly featured behaviour modes focusing on obtaining GRC for tunnels excavated in strain‐softening materials can be implemented. This way of solving such problems has proved to be more efficient and less time consuming than using FEM‐ or FDM‐based numerical 2D codes. Copyright © 2003 John Wiley & Sons, Ltd.     

Probabilistic yielding and cyclic behavior of geomaterials

In this paper, the novel concept of probabilistic yielding is used for 1‐D cyclic simulation of the constitutive behavior of geomaterials. Fokker–Planck–Kolmogorov equation‐based probabilistic elastic–plastic constitutive framework is applied for obtaining the complete probabilistic (probability density function) material response. Both perfectly plastic and hardening‐type material models are considered. It is shown that when uncertainties in material parameters are taken into consideration, even the simple, elastic‐perfectly plastic model captures some of the important features of geomaterial behavior, for example, modulus reduction with cyclic strain, which, deterministically, is only possible with more advanced constitutive models. Furthermore, it is also shown that the use of isotropic and kinematic hardening rules does not significantly improve the probabilistic material response. Copyright © 2010 John Wiley & Sons, Ltd.     

Cavity expansion in sands under drained loading conditions

Solutions for the expansion of cylindrical and spherical cavities in sands are presented. The sand is modelled using recently proposed critical-state models in which the values of the friction and dilation angles depend on the deformation history. Similarity solutions are obtained which enable the limit pressure to be calculated as a function of the initial conditions. Comparisons with existing perfectly plastic theories are made and consequences for the interpretation of cone penetrometer measurements are indicated.     

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.     

Spurious bifurcation modes for cohesive solids under uniform stress

Singularities leading to the calculation of spurious velocity fields have been observed in the tangent stiffness equations assembled during the finite-element analysis of square blocks of purely cohesive perfectly plastic material. These occur because of weaknesses in conventional implementations of the plastic flow rule. Both twisting and swaying modes occur which interfere with the calculation of uniform velocity fields for blocks responding to uniform applied stress. Various techniques can be used to eliminate the problem. The phenomenon does not eventuate when non-uniform stress fields are acting, so does not affect the solution of practical plasticity problems.     

Use of sparse iterative methods for the resolution of three-dimensional soil/structure interaction problems

In this paper we present a study of the performance of sparse iterative solvers regarding the resolution of three-dimensional and non-linear problems encountered in soil/structure interaction. It is composed of two parts. In the first one, we present briefly iterative methods and preconditioners used in this study, then we analyse their performance on three soil/structure interaction problems: a shallow foundation under a vertical loading, a single pile subjected to a lateral loading and the construction of a lined tunnel in a soft soil. Tests are performed assuming an elastic–perfectly plastic constitutive law for the soil material with a non-associated Mohr–Coulomb flow rule. Copyright © 1999 John Wiley & Sons, Ltd.     

Strain‐softening analysis of a spherical cavity

This paper studies the excavation of a spherical cavity subjected to hydrostatic initial stresses in the infinite homogeneous and isotropic rock mass with strain‐softening Mohr–Coulomb (M‐C) and Hoek–Brown (H‐B) behaviors. Numerical solutions of the spherical cavity are obtained and the application to determining stress–strain curve of strain‐softening M‐C and H‐B rock mass is studied. A closed‐form solution for the elastic–brittle–plastic medium is introduced first, and then a numerical procedure that simplifies the strain‐softening process into a series of brittle–plastic ones is presented. The approach is validated against the facts that the strain‐softening process evolves into a brittle–plastic one when the softening slope is very steep, whereas it evolves into an elasto‐plastic one when the softening slope approaches zero. Numerical solutions for the prediction of displacements and stresses around the spherical cavity in the strain‐softening M‐C and H‐B rock mass are presented. On the basis of the analysis of the spherical cavity in strain‐softening rock mass, the stress–strain relationship at an infinitesimal cube around the cavity is obtained and discussed with different evolution laws for the strength parameters considered. Copyright © 2011 John Wiley & Sons, Ltd.     

A two-mechanism soil–structure interface model for three-dimensional cyclic loading

A three-dimensional (3D) soil–structure interface model is proposed within the two-mechanism constitutive theory and bounding surface theory originally established for soils. The proposed model has two main characteristics: first, the model is formulated based on two different and superposed deformation mechanisms. The first mechanism is due to the stress ratio increment, and the second is due to the normal stress increment. Each mechanism induces a shear strain component and a normal strain component. The proposed model can be reduced to the conventional single-mechanism interface model. Second, the plastic modulus and stress dilatancy are defined using the bounding surface theory. The plastic flow rule under cyclic loading is modified and assumed to be dependent on both the stress state of the mapping point and the stress reversal loading direction. The proposed model was validated against the available 3D interface tests and was found to satisfactorily reflect the salient features of the interfaces under monotonic and cyclic loading paths with different normal boundaries. The problem in which the elastic normal stiffness in conventional single-mechanism interface models is often underestimated to enhance the simulation performance under varying normal stress conditions is solved by incorporating the second mechanism. And the effect of the second mechanism on the modeling behavior is discussed. The modified plastic flow direction accurately simulates the 3D cyclic shear response, and the difference between the model simulation and test result increases with the number of cycles by use of the plastic flow direction defined in conventional bounding surface theory.     

A finite strain closed‐form solution for the elastoplastic ground response curve in tunnelling

The ground response to tunnel excavation is usually described in terms of the characteristic line of the ground (also called ‘ground response curve’, GRC), which relates the support pressure to the displacement of the tunnel wall. Under heavily squeezing conditions, very large convergences may take place, sometimes exceeding 10–20% of the excavated tunnel radius, whereas most of the existing formulations for the GRC are based on the infinitesimal deformation theory. This paper presents an exact closed‐form analytical solution for the ground response around cylindrical and spherical openings unloaded from isotropic and uniform stress states, incorporating finite deformations and linearly elastic‐perfectly plastic rock behaviour obeying the Mohr–Coulomb failure criterion with a non‐associated flow rule. Additionally, the influence of out‐of‐plane stress in the case of cylindrical cavities under plane‐strain conditions is examined. The solution is presented in the form of dimensionless design charts covering the practically relevant parameter range. Finally, an application example is included with reference to a section of the Gotthard Base tunnel crossing heavily squeezing ground. The expressions derived can be used for preliminary convergence assessments and as valuable benchmarks for finite strain numerical analyses. Copyright © 2014 John Wiley & Sons, Ltd.     

基于西原模型的圆形隧道黏弹-黏塑性解析解

基于西原模型,假设黏塑性体的偏应变张量的一阶导数与瞬态偏应力张量和稳态偏应力张量之差成正比,得到围岩黏塑性区的本构方程。采用拉普拉斯变换与逆变换,推导了圆形隧道黏弹-黏塑性解析解。当 时,该解退化成线弹性本构模型的解答;当 时,该解退化成理想弹塑性本构模型的解答。通过工程实例,分析了围岩位移场、应力场和黏塑性区半径随时间的变化规律。当支护力保持不变时,围岩不同位置位移、围岩黏塑性区半径将随时间增长而持续增大并趋于稳定;围岩黏弹-黏塑性特征对径向应力和黏弹性区切向应力影响较小,对黏塑性区切向应力影响较大,越靠近洞壁处,切向应力随时间变化越剧烈。此外,不同支护力作用下洞壁处的切向应力在支护初期均较大,因此应采用及时支护的策略;考虑到围岩黏弹-黏塑性特征对支护力的影响,建议采取让压支护技术以保证围岩和衬砌的稳定性。     

Solutions for Displacement and Stress in Strain-Softening Surrounding Rock Incorporating the Effects of Hydraulic–Mechanical Coupling and Rockbolts Effectiveness

This study focuses on the stress and displacement of a circular opening that is excavated in a strain-softening rock mass incorporating the effects of hydraulic–mechanical coupling and rockbolts effectiveness. It follows the generalized Hoek–Brown failure criterion. Moreover, an improved numerical approach and stepwise procedure are proposed. This approach considers the deterioration of the strength, deformation, and dilation angle and the variation of elastic strain in the plastic region considering the effect of the hydraulic–mechanical coupling and the rockbolts effectiveness. The presented solutions were validated by FLAC results. Several examples are conducted to demonstrate the validity and accuracy of the proposed solution through MATLAB programming. Parametric studies are also conducted to highlight the influences of hydraulic–mechanical coupling and rockbolts effectiveness on stress and displacement. Results show that stress and displacement, incorporating the effects of hydraulic–mechanical coupling and rockbolts effectiveness, are between those when hydraulic–mechanical coupling or rockbolts effectiveness is considered separately. However, this theory needs more verification from practical engineering.