Implicit and explicit integration schemes in the anisotropic bounding surface plasticity model for cyclic behaviours of saturated clay

Two integration algorithms, namely the implicit return mapping and explicit sub-stepping schemes, are adopted in the anisotropic bounding surface plasticity model for cyclic behaviours of saturated clay and are implemented into finite element code. The model is a representative of a series of bounding surface models that have typical characteristics, including isotropic and kinematic hardening rules and a rotational bounding surface to capture complex but important cyclic behaviours of soils, such as cyclic shakedown and degradation. However, there is no explicit current yield surface in the model to which the conventional implicit algorithm returns the stress state back or the sub-stepping integration corrects the drift of the stress state. Hence, necessary modifications have been made for both of the integration schemes. First, the image stress point is mapped or corrected to the bounding surface instead of mapping back or correcting the stress state to the yield surface. Second, the unloading–loading criterion is checked to determine the image stress point rather than checking the yield criterion after giving the trial stress state in a conventional way. Comparative studies on the accuracy, stability and efficiency of the two integration schemes are conducted not only at the element level but also in solving boundary value problems of monotonic and cyclic bearing behaviours of rigid footings on saturated clay. For smaller strain increments, there is no significant difference in the accuracy between the two integration schemes, but the explicit integration shows a higher efficiency and accuracy. For relatively larger increments, the implicit return mapping algorithm presents good accuracy and more robustness, while the sub-stepping algorithm shows deteriorating accuracy and suffers the convergence problem. With the tolerance used in the present model, the bearing capacity of the rigid footing predicted by the return mapping algorithm is closer to the available analytical and numerical solutions, while the bearing capacity predicted by the sub-stepping algorithm shows a marginal increase.     

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

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

The critical state behaviour of barodesy compared with the Matsuoka–Nakai failure criterion

Barodesy is a new approach to constitutive modelling of soil. It is based on Goldscheider's principles and maps stretching directions onto corresponding stress directions with the help of a simple exponential function. This mapping also determines a critical state surface in principal stress space. The article investigates this surface and relates it to the well‐known Matsuoka–Nakai failure criterion. It turns out that the difference between these two surfaces is negligible for practical applications. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling of rutting of two flexible pavements with the shakedown theory and the finite element method

This paper presents a finite element program, for the modelling of rutting of flexible pavements. In its present version, the program incorporates a permanent deformation model for unbound granular materials based on the concept of the shakedown theory developed by Zarka for metallic structures under cyclic loadings and has been used to estimate the permanent deformations of unbound granular materials (UGM) subjected to traffic loading. The calculation is performed in two steps: the first step consists in modelling the resilient behaviour of the pavement in 3D, using non-linear elastic models, to determine the stress field in the pavement. Then stress paths are derived and used to calculate the permanent deformations and the displacements, using a Drucker–Prager yield surface. An application to the prediction of the permanent deformations of experimental pavements with an unbound granular base, tested on the LCPC pavement testing facility is presented.     

A METHOD TO CORRECT YIELD SURFACE DRIFT IN SOIL PLASTICITY UNDER MIXED CONTROL AND EXPLICIT INTEGRATION

When applying an explicit integration algorithm in e.g. soil plasticity, the predicted stress point at the end of an elastoplastic increment of loading might not be situated on the updated current yield surface. This so-called yield surface drift could generally be held under control by using small integration steps. Another possibility, when circumstances might demand larger steps, is to adopt a drift correction method. In this paper, a drift correction method for mixed control in soil plasticity, under drained as well as undrained conditions, is proposed. By simulating triaxial tests in a Constitutive Driver, the capability and efficiency of this correction method, under different choices of implementation, have been analysed. It was concluded that the proposed drift correction method, for quite marginal additional computational cost, was able to correct successfully for yield surface drift giving results in close agreement to those obtained with a very large number of integration steps. © 1997 by John Wiley & Sons, Ltd.     

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

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

三维弹塑性模型在路堤软基固结分析中应用

针对基于SMP准则改进的剑桥模型(MCC-SMP),采用弹性预测和塑性修正两阶段的回映算法,给出了应力更新算法,根据其算法编制了非线性有限元子程序,进行了三轴压缩、三轴伸长及平面应变等单元试验问题的数值模拟,通过数值分析和试验结果的对比,验证了子程序的精度;采用该本构模型,对澳大利亚某试验段路堤软土地基固结进行了三维有限元分析,重点分析和对比了孔压、沉降及水平位移的变化规律。     

A fast algorithm for finding the first intersection with a non-convex yield surface

A major task in the numerical modelling of soils using complex elastoplastic material models is stress updating. This paper proposes a fast and robust numerical algorithm for locating the first intersection between a non-convex yield surface and an elastic trial stress path. The intersection problem is cast into a problem of finding the smallest positive root of a nonlinear function. Such a function may have multiple roots within the interval of interest. The method is based on the modified Steffensen method, with important modifications to address the issues arising from the non-convexity. Numerical examples demonstrate that the proposed M² Steffensen method is indeed computationally efficient and robust.     

A middle surface concept (MSC) model for saturated sands in general stress space

An elastoplastic constitutive model is proposed for saturated sands in general stress space using the middle surface concept (MSC). In MSC, different features of stress–strain response of a material are divided into different pseudo‐yield surfaces. The true‐yield surface representing the true response is established by using various links between the yield surfaces. In this MSC sand model, several well‐known features of sand response are represented by three different pseudo‐yield surfaces, which are developed in a simple and straightforward way. These features include the critical state behaviour, the effects of state parameter, unloading and reloading plastic deformation, the influence of fabric anisotropy, and phase transformation line related behaviour. Finally, the model predictions and test results are compared for two different types of sands under a variety of loading conditions and good comparisons are obtained. The application of MSC to saturated sand modelling shows the versatility of MSC as a general concept for modelling stress–strain response of materials. Copyright © 2006 John Wiley & Sons, Ltd.     

基于Drucker-Prager准则的岩石弹塑性损伤本构模型研究

大多数岩石材料软化本构模型在硬化函数中引入塑性内变量来表示材料的硬化/软化性质,但并不能反映岩石微裂隙损伤对材料力学性能的影响及单轴拉伸和压缩所表现的初始屈服强度f0与屈服极限fu的差异。基于D-P准则同时考虑塑性软化及损伤软化,建立岩石类材料的弹塑性本构关系及其数值算法。塑性屈服函数采用Borja等的应力张量的硬化/软化函数,反映塑性内变量及应力状态对硬化函数的影响;由于岩石损伤软化是微裂隙扩展所导致的体积膨胀引起的,因此,提出用体积应变表征岩石损伤变量的演化,并用回映隐式积分算法编制了岩石的弹塑性损伤本构程序。对单轴压缩及拉伸荷载作用下的岩石材料试验进行数值模拟,结果表明,所提出的岩石弹塑性损伤本构模型可以较好地符合岩石材料的力学特性。     

广义塑性力学中的屈服面与应力-应变关系

详细讨论了广义塑性力学中屈服面和塑性势面的对应关系以及岩土材料的三类屈服面(即体积屈服面与 q 方向上及方向上的剪切屈服面)的基本特征, 尤其是提出了能考虑剪胀与剪缩的体积屈服面和应力 Lode 角θσ方向的剪切屈服面。指出在广义塑性力学中不必采用硬化定律, 就能得出塑性应变增量与应力增量的关系, 给出了求弹塑性矩阵的方法。     

单轴压缩非贯通节理岩石尖端塑性区及弹塑性断裂模型研究

基于Drucker-Prager（下简称D-P）准则,建立压缩载荷作用下的非贯通节理岩石的弹塑性断裂模型。针对节理岩石小范围屈服翼裂纹尖端塑性区,推导了D-P屈服准则的纯I、纯II及I、II复合型3种翼裂纹无量纲塑性区径长函数,并与Mises准则的塑性区进行对比;结果表明,D-P准则的I型和复合型塑性区较Mises屈服准则的塑性区大,且其II型及I、II复合型塑性区在翼裂纹上下表面不连续。进一步,引入断裂软化因子以表征节理岩石裂隙断裂扩展后的断裂软化规律,考虑非贯通节理岩石复合型断裂软化,是由于节理尖端翼裂纹应变能密度超过最小应变能密度导致其成核扩展引起的,提出用应变能密度的指数函数形式表征断裂软化变量的演化;塑性屈服函数采用Borja等的应力张量3个不变量的硬化/软化函数,反映塑性内变量及应力状态对硬化函数的影响;建立节理岩石的弹塑性断裂本构关系及其数值算法,并用回映隐式积分算法编制了弹塑性断裂模型的程序。以单轴压缩下非贯通节理岩石为例,分析岩石断裂韧度、节理摩擦系数和节理倾角等参数的影响,结果表明,所提出的弹塑性断裂模型与数值和试验结果比较吻合。     

Monica—A computer algorithm for solving boundary value problems using the double-hardening constitutive law Monot: I. Algorithm development

The principal features involved in the implementation of the double-hardening constitutive law Monot¹ into a general purpose computer algorithm called MONICA^2,3 are described. These include details of general program structure and of steps taken to overcome problems such as computer storage, computer run-time, algorithm stability and problems associated with the stress–strain singularity which exists at the intersection of the two yield surfaces.     

An approach for providing quasi‐convexity to yield functions and a generalized implicit integration scheme for isotropic constitutive models based on 2 unknowns

Yield and plastic potential surfaces are often affected by problems related to convexity. One such problem is encountered when the yield surface that bounds the elastic domain is itself convex; however, convexity is lost when the surface expands to pass through stress points outside the current elastic domain. In this paper, a technique is proposed, which effectively corrects this problem by providing linear homothetic expansion with respect to the centre of the yield surface. A very compact implicit integration scheme is also presented, which is of general applicability for isotropic constitutive models, provided that their yield and plastic potential functions are based on a separate mathematical definition of the meridional and deviatoric sections and that stress invariants are adopted as mechanical quantities. The elastic predictor‐plastic corrector algorithm is based on the solution of a system of 2 equations in 2 unknowns only. This further reduces to a single equation and unknown in the case of yield and plastic potential surfaces with a linear meridional section. The effectiveness of the proposed convexification technique and the efficiency and stability of the integration scheme are investigated by running numerical analyses of a notoriously demanding boundary value problem.     

基于极限平衡有限元的岩质高边坡安全系数长期预测分析

引入基于Dijkstra算法的极限平衡有限元方法到岩质边坡长期稳定性研究中,得到边坡长期安全系数,并将该方法应用于糯扎渡工程坝肩高边坡,在弹塑性下比较刚体极限平衡和基于Dijkstra极限平衡有限元方法计算结果,认为该算法可信。基于糯扎渡水电站坝肩高边坡有限元流变分析中各时步应力场结果,采用基于Dijkstra极限平衡有限元方法得到各时步潜在滑动面和安全系数。由此判定该边坡长期稳定性,并给出安全系数随时间变化趋势和长期安全系数。     

Hoek-Brown准则的主应力回映算法及其二次开发

在弹塑性模型数值实现的过程中,需要进行应力更新的回映算法。针对三维应力空间回映算法在奇异点收敛性方面的不足,提出主应力空间的回映算法,讨论了算法实现过程的应力空间转化问题,分析了应力更新过程中确定回映区域的方法,建立了相应的一致性刚度矩阵。基于大型有限元软件ABAQUS提供的用户材料子程序UMAT接口,利用 Fortran 编程语言,实现了Hoek-Brown准则主应力空间的隐式积分算法,利用开发的模型,进行了岩石常规三轴压缩试验的数值模拟,通过与ABAQUS内部的Mohr-Coulomb准则计算结果的对比,验证了模型和程序的可行性和准确性     

土的三重屈服面应力-应变模型与SMP准则的结合

土的三重屈服面应力应变模型在国内有一定的影响力,它把土的塑性应变分成3部分,每一部分对应一个屈服面,分别为压缩屈服面、剪切屈服面和剪胀屈服面,从而在一定程度上反应了土的基本特性。但它在由三轴压缩应力状态向其它应力状态转化时,塑性系数需要重新确定,比较繁琐。而SMP准则能合理反映土的破坏特性,在国际上有一定的影响力。借用应力变换三维化方法,通过把土的三重屈服面应力应变模型和SMP准则相结合,使得原有的模型在不做任何假设的条件下,采用统一的塑性系数,由三轴压缩应力状态简单地转化到一般应力状态,并且能够考虑应变的分叉特性,能合理地预测已有的试验数据。     

山区谷底沉陷预测模型及其参数反演

我国西部地区地势复杂,沟壑纵横,地下开采极易导致边坡失稳,引发采动滑坡。在地下采动沉降与滑坡体挤压上升的叠加影响下,谷底区域地表沉降值明显小于类似地质采矿条件下的平原地区。为准确预测山区谷底区域地表沉降值,基于简支梁的弹性变形理论,并借助概率密度函数建立了山区谷底区域地表沉陷预计修正模型,明确模型参数物理意义及其取值方法。依据修正模型,以实测值和预测值之差平方和最小为原则构建适应值函数,基于模拟退火粒子群算法提出新的模型参数反演方法,借助MATLAB语言编制了相应的参数反演程序。最后将研究成果应用于山西某矿,得到谷底区域预测结果中误差为73 mm,与实测值基本一致,取得了较好的工程实践效果。     

An improved implicit numerical integration of a non‐associated,three‐invariant cap plasticity model with mixed isotropic–kinematic hardening for geomaterials

The Sandia GeoModel is a continuum elastoplastic constitutive model that captures many features of the mechanical response for geological materials over a wide range of porosities and strain rates. Among the specific features incorporated into the formulation are a smooth compression cap, isotropic/kinematic hardening, nonlinear pressure dependence, strength differential effect, and rate sensitivity. This study attempts to provide enhancements regarding computational tractability, domain of applicability, and robustness of the model. A new functional form is presented for the yield and plastic potential functions. This reformulation renders a more accurate, robust, and efficient model as it eliminates spurious solutions attributed to the original form. In addition, we achieve a high‐performance implementation, because the local iterative method is allowed to recast residual vectors with a uniform dimensionality. The model is also furnished with a smooth, elliptical tension cap to account for the tensile yielding. Moreover, an efficient algorithm is introduced, which decreases the computational cost by differentiating the updated shear yield surface from the cap surface based on the trial relative stress state. Finally, various numerical examples including a large‐scale boundary value problem are presented to demonstrate the fidelity of the modified model and to analyze its numerical performance. Copyright © 2015 John Wiley & Sons, Ltd.     

A unified approach to the implicit integration of standard,non‐standard and viscous plasticity models

It is shown how modern concepts to integrate the elasto‐plastic rate equations of standard plasticity via an implicit algorithm can be generalized to plasticity without an explicitly defined yield surface and to overstress‐type models of viscoplasticity, where the stress point can be located outside the loading surface. For completeness, a tangent operator is derived that is consistent with the update algorithm. Copyright © 2002 John Wiley & Sons, Ltd.