一种适用于饱和黏土循环动力分析边界面塑性模型的隐式积分算法

修正了传统隐式回映算法,建立了适用于饱和黏土循环动力分析的边界面塑性模型的完全隐式积分格式。该模型基于无弹性域概念和临界状态理论,采用各向同性、运动硬化准则、旋转的边界面,并引入表征土体结构损伤和重塑程度的损伤变量以反映循环载荷作用下饱和黏土的各向异性、刚度、强度软化及塑性变形累积等特征。针对等压固结 和偏压固结 的饱和高岭黏土的不排水三轴试验进行模拟,采用不同的应变增量步长进行计算,并与试验数据对比,结果表明,修正隐式回映算法应用于该类边界面模型的合理性、积分格式的精确性和稳定性;另外,结合有限元软件自动时间步长的增量迭代解法,对饱和黏土应力控制的不排水动三轴试验进行预测,结果表明,修正的适用于该边界面的塑性模型隐式回映算法可以得到比较合理的数值分析结果,能够反映饱和黏土的循环刚度的退化和强度的弱化等动力特性。     

On integration of a cyclic soil plasticity model

Performance of three classes of explicit and implicit time‐stepping integrators is assessed for a cyclic plasticity constitutive model for sands. The model is representative of an important class of cyclic plasticity models for soils and includes both isotropic and nonlinear kinematic hardening. The implicit algorithm is based on the closest point projection method and the explicit algorithm follows a cutting‐plane integration procedure. A sub‐stepping technique was also implemented. The performance of these algorithms is assessed through a series of numerical simulations ranging from simulations of laboratory tests (such as triaxial and bi‐axial compression, direct shear, and cyclic triaxial tests) to the analysis of a typical boundary value problem of geotechnical earthquake engineering. These simulations show that the closest point projection algorithm remains stable and accurate for relatively large strain increments and for cases where the mean effective stress in a soil element reaches very small values leading to a liquefaction state. It is also shown that while the cutting plane (CP) and sub‐stepping (SS) algorithms provide high efficiency and good accuracy for small to medium size strain increments, their accuracy and efficiency deteriorate faster than the closest point projection method for large strain increments. The CP and SS algorithms also face convergence difficulties in the liquefaction analysis when the soil approaches very small mean effective stresses. Copyright © 2001 John Wiley & Sons, Ltd.     

粗粒土边界面塑性模型及其积分算法

粗粒土的强度、变形特性对土石坝、边坡和路基等工程的安全性与稳定性有着至关重要的影响。针对粗粒土在复杂应力状态下的强度和变形特性,在边界面塑性理论和临界状态理论框架下,通过引入状态参数和动态临界状态线建立了粗粒土状态相关边界面塑性模型。模型不仅能够模拟粗粒土的应变硬化和体积收缩行为,还能描述应变软化和体积膨胀特性。基于ABAQUS的二次开发平台,结合带误差控制的改进欧拉积分算法编写了边界面塑性模型的UMAT子程序。通过设置不同的应变增量步和误差容许值,对改进欧拉积分算法的精确性和收敛性进行了分析。对不同密实状态和围压下粗粒土三轴排水剪切试验进行了模拟,验证了带误差控制的改进欧拉积分算法应用于粗粒土边界面塑性模型的合理性,为进一步工程应用奠定了基础。     

A virtual surface concept for nested yield surface plasticity

A return mapping algorithm based on the concept of virtual surfaces is presented for stress integration in the numerical solution of boundary value problems using an elastoplastic multiple yield surface model. Both kinematic hardening and isotropic softening yield surfaces are addressed. With the implementation of virtual surfaces, the active yield surface is located before a stress update is implemented. The return mapping algorithm consists of only two steps. The first step consists of performing the intermediate update of the trial stress corresponding to the translation of all yield surfaces inside the active yield surface. Little calculation is involved in this step. The second step consists of performing the conventional explicit one-step radial return mapping of the updated trial stress on the corresponding active yield surface. The algorithm efficiency resides in the virtual surface based return mapping which preserves the elastic predictor representation of the updated trial stress, and limits iterations to at most two steps.     

考虑饱和黏土初始各向异性的循环边界面塑性模型

天然土体的初始各向异性通常可对其后继循环特性产生显著影响。现有考虑循环载荷作用的土体弹塑性模型,往往采用类似修正剑桥模型的椭圆形屈服面,已有研究表明,该椭圆形屈服面因其拉伸弹性区域偏大,针对天然K0固结状态的土体,其计算精度较差。基于新近提出的广义各向同性硬化准则,在边界面方程中引入初始各向异性张量,并采用空间滑动面破坏准则（SMP）的变换应力法,建立了能考虑饱和黏土初始各向异性的循环边界面塑性模型。分别针对等压和偏压固结的饱和黏土静、动三轴试验进行模拟,结果表明,该模型能合理反映土体的初始各向异性及其后继循环动力特性。     

Anisotropic bounding-surface plasticity model for the cyclic shakedown and degradation of saturated clay

The cyclic behaviours of embedded offshore structures under different cyclic loading levels are related to the cyclic shakedown and degradation of the surrounding soils. In the present study, a damage-dependent bounding-surface model based on a newly proposed hardening rule was developed to predict the cyclic shakedown and degradation of saturated clay and the effect of the initial anisotropic stress state. By extending the Masing’s rule to the bounding-surface plasticity theory, the stress reversal point is taken as the generalised homological centre of the bounding surface. With movement of the generalised homological centre, at lower stress amplitudes, the cyclic process ends at a steady state, and cyclic shakedown is reached. At higher stress amplitudes, a damage parameter related to the accumulated deviatoric plastic strain is incorporated into the form of the bounding surface, which is hence able to contract to model degradations in stiffness and strength. To take into account the effects of initial anisotropic conditions on the cyclic behaviour of soils, an initial anisotropic tensor is introduced in the bounding surface. The developed model is validated through undrained isotropic and anisotropic cyclic triaxial tests in normally consolidated and overconsolidated saturated clay under both one-way and two-way loadings. Both cyclic shakedown and degradation are well reproduced by the model, as is the anisotropy effect induced by the initial anisotropic consolidation process.     

纯摩擦型岩土介质本构积分算法及其在ABAQUS中开发应用

提出了基于Euler向后积分的自适应子增量本构积分算法,推导了相应的一致切线模量矩阵;通过引入伪屈服函数(塑性势函数),提出了对屈服面角点应力区进行两个方向应力投射的本构积分算法,使超出屈服面的试应力收敛到角点;推导了两个投射方向的一致切线模量矩阵;采用赋小值方法解决0应力屈服的问题。用上述方法编制了基于D-P准则的理想弹塑性模型ABAQUS-UMAT子程序,并进行了算例验证。     

Explicit stress integration of complex soil models

In this paper, two complex critical‐state models are implemented in a displacement finite element code. The two models are used for structured clays and sands, and are characterized by multiple yield surfaces, plastic yielding within the yield surface, and complex kinematic and isotropic hardening laws. The consistent tangent operators—which lead to a quadratic convergence when used in a fully implicit algorithm—are difficult to derive or may even not exist. The stress integration scheme used in this paper is based on the explicit Euler method with automatic substepping and error control. This scheme employs the classical elastoplastic stiffness matrix and requires only the first derivatives of the yield function and plastic potential. This explicit scheme is used to integrate the two complex critical‐state models—the sub/super‐loading surfaces model (SSLSM) and the kinematic hardening structure model (KHSM). Various boundary‐value problems are then analysed. The results for the two models are compared with each other, as well with those from standard Cam‐clay models. Accuracy and efficiency of the scheme used for the complex models are also investigated. Copyright © 2005 John Wiley & Sons, Ltd.     

Kinematic hardening of soft clay in simple shear

In a separate paper, the authors have proposed a normalized, non-degrading form of the shear stress–shear strain relationship for undrained, cyclic simple shear of soft clay. This relationship is described in the present paper, and it is seen to include a single fatigue parameter—the mean effective stress. Application of the relationship therefore requires knowledge of the history of the mean effective stress during any loading history. The present paper proposes an effective stress path model which may be used for prediction of this history. The model is developed within the framework of bounding surface kinematic and isotropic hardening plasticity. It incorporates an isotropic hardening bounding surface, and a kinematic hardening yield surface, in which the elastic region vanishes, and so the yield surface reduces to the stress point. The normalized shear stress–shear strain relationship, developed on the basis of Iwan's model, is used to establish the shape of the cap of the bounding surface. A new translation rule is also incorporated in the model, allowing improved prediction of stress path development within the bounding surface during regular or irregular cyclic loading. Use of the proposed model to simulate the behaviour of soft clay in laboratory undrained cyclic simple shear tests shows excellent qualitative agreement, with most of the major features of the actual behaviour being predicted.     

Automatic substepping integration of the subloading tij model with stress path dependent hardening

This paper discusses the numerical integration of the subloading tij model. This is an elastoplastic model with stress path dependent hardening, which can predict the behaviour of normally consolidated clays or loose sands, as well as of over-consolidated clays or dense sands, with a small number of material parameters. Three features distinguish the subloading tij model from the conventional ones: (a) the use of a modified stress space given by tensor t_ij; (b) the split of the plastic strain increments in two components leading to a stress path dependent hardening; and (c) the use of two yield surfaces (subloading yield surface and normal yield surface). This last feature is based on the concept of sub-yielding stress states and adds an extra internal strain-like hardening variable, related to the relative density state, which demands its own evolution law. The three characteristics above greatly improve the prediction capabilities of the model, with respect to those of the well-known Cam clay model, at the cost of only two additional parameters. Nonetheless, the numerical integration of the constitutive equations of subloading tij model is a bit challenging, mainly due to the stress path dependent hardening. In order to integrate the equations of subloading tij model in the same way as for any conventional model, the authors reformulated its equations in a simpler and direct manner. Here, these equations are integrated using multi-step explicit schemes, such as modified-Euler and Runge–Kutta–Dormand–Price, with automatic error control. Simple forward-Euler scheme is also used for the sake of comparison. The results show that the modified-Euler scheme is more accurate as well as faster than the other schemes analysed over a wide range of error tolerance. Besides, the automatic feature of these schemes is a great convenience for the users of numerical codes.     

An elastoplastic model with combined isotropic–kinematic hardening to predict the cyclic behavior of stiff clays

This paper presents a kinematic hardening model for describing some important features of natural stiff clays under cyclic loading conditions, such as closed hysteretic loops, smooth transition from the elastic behavior to the elastoplastic one and changes of the compression slope with loading/unloading loops. The model includes two yield surfaces, an inner surface and a bounding surface. A non-associated flow rule and a kinematic hardening law are proposed for the inner surface. The adopted hardening law enables the plastic modulus to vary smoothly when the kinematic yield surface approaches the bounding surface and ensures at the same time the non-intersection of the two yield surfaces. Furthermore, the first loading, unloading, and reloading stages are treated differently by applying distinct hardening parameters. The main feature of the model is that its constitutive equations can be simply formulated based on the consistency condition for the inner yield surface based on the proposed kinematic hardening law; thereby, this model can be easily implemented in a finite element code using a classic stress integration scheme as for the modified Cam Clay model. The simulation results on the Boom Clay, natural stiff clay, have revealed the relevance of the model: a good agreement has been obtained between simulations and the experimental results from the tests with different stress paths under cyclic loading conditions. In particular, the model can satisfactorily describe the complex case of oedometric conditions where the deviator stress is positive upon loading (compression) but can become negative upon unloading (extension).     

考虑循环载荷下饱和黏土软化的损伤边界面模型研究

研究表明,循环载荷作用下饱和黏土将发生软化,其机制主要有两个：一是孔压的积累;二是土体原有结构的不断损伤和新结构的不断重塑。针对上述机制,基于广义各向同性硬化准则建立了考虑饱和黏土循环软化的损伤单面模型。该模型在有效应力空间中引入损伤变量,表征土体结构的损伤和重塑程度,在连续的循环加载下,损伤不断累积,边界面则随着损伤的累积不断收缩,以模拟饱和黏土刚度和强度的软化;以应力反向点作为边界面的广义各向同性硬化中心和映射法则的映射中心,灵活地选择塑性模量的插值公式以模拟塑性变形和孔压的累积以及应力-应变的滞回特性。应用该模型对不排水循环三轴试验进行模拟,并且考查了循环周次、循环应力水平和固结历史对饱和黏土循环软化特性的影响,并与相关试验比较,验证了模型的有效性。     

岩石弹塑性损伤本构模型建立及在隧道工程中的应用

在实际隧道施工过程中,隧道开挖引起地下岩体应力重分布使得围岩的微裂纹扩展损伤,并伴随有塑性流动变形。在地下水环境中对于孔隙和微裂隙围岩介质受到应力作用时,在内部将产生高孔隙水压力影响岩石的力学性质,也改变了围岩的破坏模式。为了研究损伤引起的刚度退化和塑性导致的流动两种破坏机制的耦合作用,从弹塑性力学和损伤理论的角度出发,同时引入修正有效应力原理来考虑孔隙水压力的作用,建立基于Drucker-Prager屈服准则的弹塑性损伤本构模型;针对该本构模型推导了孔隙水压力作用下弹塑性损伤本构模型的数值积分算法-隐式返回映射算法,分别对预测应力返回到屈服面的光滑圆锥面或尖点奇异处两种可能的情况给出了详细的描述,隐式返回映射算法具有稳定性和准确性的特点;大多数弹塑性损伤模型中涉及参数多且不易确定的问题,采用反分析方法获得损伤参数,解决了损伤参数不易确定的难题;采用面向对象的编程方法,使用C++语言编制了弹塑性损伤本构求解程序,并对所建立的弹塑性损伤模型和所编程序进行了试验和数值两个方面的验证;最后将其在吉林抚松隧道工程中进行应用,模拟了塑性区和损伤区的发展变化。研究结果表明：所建立的弹塑性损伤本构模型能够较好地描述岩石的力学性能、塑性和损伤变化趋势,所编程序能够进行实际工程问题的模拟,对现场施工给予一定的指导。     

Numerical integration and FE implementation of a hypoplastic constitutive model

Hypoplastic constitutive equation based on nonlinear tensor functions possesses a failure surface but no yield surface. In this paper, we consider the numerical integration and FE implementation of a simple hypoplastic constitutive equation. The accuracy of several integration methods, including implicit and explicit methods, is examined by performing a set of triaxial compression tests. Adaptive explicit schemes show the best performance. In addition, the stress drift away from the failure surface is corrected with a predictor-corrector scheme, which is verified by two boundary value problems, i.e. rigid footing tests and slope stability.     

基于上下负荷面的弹黏塑性本构模型及应力积分算法实现

为了描述天然软土的时间相依以及结构性特征,提出了一种能考虑土体超固结和结构性的实用弹黏塑性本构模型。它以Asaoka和Hashiguchi的上下负荷面作为某一应变速率下的参考屈服面,按照相对过应力的基本思路,新引入了两个能通过不同应变速率三轴压缩试验测定的率敏性参数 和 ,建立了以当前应力、黏塑性应变以及黏塑性应变速率为状态变量的动屈服准则函数,并给出了基于Newton-Raphson迭代的应力积分算法,且成功地将其嵌入到大型有限元软件ABAQUS中。最终通过数值算例来验证模型的正确性以及应力积分算法的可靠性。结果表明：该模型能同时描述土体的率敏性、蠕变以及结构性特征,模型参数物理意义明确、易懂可测,预测结果与试验数据吻合良好,可用于复杂边值问题的有限元计算。     

Development of an implicit material point method for geotechnical applications

An implicit material point method (MPM), a variant of the finite element method (FEM), is presented in this paper. The key feature of MPM is that the spatial discretisation uses a set of material points, which are allowed to move freely through the background mesh. All history-dependent variables are tracked on the material points and these material points are used as integration points similar to the Gaussian points. A mapping and re-mapping algorithm is employed, to allow the state variables and other information to be mapped back and forth between the material points and background mesh nodes during an analysis. In contrast to an explicit time integration scheme utilised by most researchers, an implicit time integration scheme has been utilised here. The advantages of such an approach are twofold: firstly, it addresses the limitation of the time step size inherent in explicit integration schemes, thereby potentially saving significant computational costs for certain types of problems; secondly, it enables an improved algorithm accuracy, which is important for some constitutive behaviours, such as elasto-plasticity. The main purpose of this paper is to provide a unified MPM framework, in which both quasi-static and dynamic analyses can be solved, and to demonstrate the model behaviour. The implementation closely follows standard FEM approaches, where possible, to allow easy conversion of other FEM codes. Newton’s method is used to solve the equation of motion for both cases, while the formation of the mass matrix and the required updating of the kinematic variables are unique to the dynamic analysis. Comparisons with an Updated Lagrangian FEM and an explicit MPM code are made with respect to the algorithmic accuracy and time step size in a couple of representative examples, which helps to illustrate the relative performance and advantages of the implicit MPM. A geotechnical application is then considered, illustrating the capabilities of the proposed method when applied in the geotechnical field.     

Explicit integration of bounding surface model for the analysis of earthquake soil liquefaction

This paper presents a new plasticity model developed for the simulation of monotonic and cyclic loading of non‐cohesive soils and its implementation to the commercial finite‐difference code FLAC, using its User‐Defined‐Model (UDM) capability. The new model incorporates the framework of Critical State Soil Mechanics, while it relies upon bounding surface plasticity with a vanished elastic region to simulate the non‐linear soil response. Stress integration of constitutive relations is performed using a recently proposed explicit scheme with automatic error control and substepping, which so far has been employed in the literature only for constitutive models aiming at monotonic loading. The overall accuracy of this scheme is evaluated at element level by simulating cyclic loading along complex stress paths and by using iso‐error maps for paths involving change of the Lode angle. The performance of the new constitutive model and its stress integration scheme in complex boundary value problems involving earthquake‐induced liquefaction is evaluated, in terms of accuracy and computational cost, via a number of parametric analyses inspired by the successful simulation of the VELACS centrifuge Model Test No. 2 studying the lateral spreading response of a liquefied sand layer. Copyright © 2009 John Wiley & Sons, Ltd.     

Implicit numerical integration for a kinematic hardening soil plasticity model

Soil models based on kinematic hardening together with elements of bounding surface plasticity, provide a means of introducing some memory of recent history and stiffness variation in the predicted response of soils. Such models provide an improvement on simple elasto‐plastic models in describing soil behaviour under non‐monotonic loading. Routine use of such models requires robust numerical integration schemes. Explicit integration of highly non‐linear models requires extremely small steps in order to guarantee convergence. Here, a fully implicit scheme is presented for a simple kinematic hardening extension of the Cam clay soil model. The algorithm is based on the operator split methodology and the implicit Euler backward integration scheme is proposed to integrate the rate form of the constitutive relations. This algorithm maintains a quadratic rate of asymptotic convergence when used with a Newton–Raphson iterative procedure. Various strain‐driven axisymmetric triaxial paths are simulated in order to demonstrate the efficiency and good performance of the proposed algorithm. Copyright © 2001 John Wiley & Sons, Ltd.     

A simple procedure to improve the explicit integration of Cam-Clay models

A simple procedure is proposed to reduce the error in the calculation of the stresses when simulating stress paths applying explicit stress integration schemes in Cam-Clay models with an associated flow rule. The procedure consists of defining the value of the preconsolidation stress (the hardening parameter of a Cam-Clay model) so that the drift of the yield surface is zero at each calculation step. The results demonstrate the efficiency of the proposed procedure, especially along softening paths. The method can easily be implemented in existing explicit integration modules, which facilitates its practical application.