非饱和黏土状态相关本构模型的数值实现

基于双重塑性机制,即剪切滑移和加载湿陷机制,及毛细滞回和变形间的耦合作用,提出了一个适用于非饱和黏土的状态相关本构模型。建立了该模型的隐式积分算法,并处理了双重固相屈服面交点处的应力积分问题。考虑了吸力与应力间的关系,推导了本构关系的一致性切线模量。最后,利用该模型预测了一组非饱和黏土的三轴试验,以此反映模型的描述能力;对比了不同应变步长下应变控制试验的计算结果,并以此验证算法的收敛性和准确性;利用有限元程序平台U_DYSAC2计算了二维情况下非饱和土地区地基固结问题,从而验证所得一致性切线模量的有效性。     

一种新的动态回弹模量有限元实现方法

美国在2004年通过NCHRP 1-37A项目提出了用于路基土评估的新的动态回弹模量模型。该模型由于考虑了体应力和偏应力的影响,已得到广泛应用。然而,目前对于该模型的有限元实现均采用基于等效动三轴试验的等效切线模量进行材料点的局部迭代。为克服该实现方法在复杂应力状态下因模型简化带来的偏差,基于广义虎克定律推导了该回弹模量模型的精确一致切线刚度矩阵。通过编写的ABAQUS用户自定义材料子程序（UMAT）对轴压和围压加载情形进行有限元模拟,发现新推导的有限元实现方法较现有基于等效三轴试验的简化方法具有更高的计算效率和精度。最后对典型沥青路面结构的分析表明,为提高复杂应力状态下路基响应计算精度,在路面结构分析中采用基于当前应力状态的一致切线刚度矩阵代替等效切线刚度矩阵是十分必要的。     

弹塑性无厚度接触面单元的本构积分算法及验证

对基于理想弹塑性理论框架、屈服准则为Mohr-Coulomb准则、采用非关联流动法则的无厚度接触面单元的本构积分算法进行了探讨,引入非关联的伪屈服函数和伪势函数,提出了将超出屈服面、处于角点应力区的试应力双向返回到屈服面的本构积分算法。据此编制了ABAQUS用户单元子程序,进行了算例验证。结果表明,提出的算法可以较好地实现土与结构物共同作用的有限元数值模拟。     

基于半隐式积分方案统一强度理论ANSYS二次开发

目前,岩土工程中常用的塑性屈服准则多为局部不可导的分段函数,在隐式有限元求解中均采用显式向前欧拉法更新应力和Newton-Raphson迭代法求解方程。但显式向前欧拉法容易导致应力偏离屈服面和计算不收敛,而Newton-Raphson迭代法需频繁求导,对于分段函数而言这些方法均不利于求解。为此,提出了采用半隐式向后欧拉法和免导数的Steffensen迭代法联合更新应力与一致切线模量的方案。根据上述方案,利用ANSYS提供的User Mat模块编写了基于统一强度理论与拉破坏复合屈服准则的理想塑性模型。用自定义的本构模型模拟深埋圆形巷道弹塑性开挖过程,数值模拟结果与理论解很好地吻合,验证了模型代码和所提出方案的正确性和可行性。采用半隐式向后欧拉法与Steffensen迭代法相结合的手段,可简化应力与一致切线模量的求解,该方法避免了应力偏离屈服面和求解流动函数偏导,算法简便,易于推广应用。     

亚塑性模型不同积分算法的数值实现

亚塑性模型是以Jaumann应力率张量及变形率张量描述的一种率型本构关系,本构关系在非线性有限元分析计算中具有关键作用,解决用应变增量求解应力响应的问题需要一个时间积分过程。针对亚塑性本构模型发展了自适应隐式和显式两种不同的积分算法,给出了误差控制的方法,同时推导了自适应隐式积分算法所需的一致切向模量,并采用了两个不同的单元,利用ABAQUS平台比较了两种积分算法的数值模拟结果。为了实现从ABAQUS/Standard到ABAQUS/Explicit的过渡,开发了UMAT-VUMAT接口,从而可以使已有的UMAT子程序用于大变形动力问题分析。算例分析证明了研究结果的正确性。     

一种高温下混凝土化学塑性-损伤耦合本构模型

提出了一个用于模拟高温下混凝土化学塑性-损伤耦合本构行为的数值模型。发展了一个化学塑性-损伤耦合分析的一致性应力返回映射算法。为了保证对于全局守恒方程Newton迭代过程的2阶收敛率,推导并形成了一致性切线模量矩阵。数值算例显示了文中发展的化学塑性-损伤耦合本构模型在模拟高温下混凝土中复杂破坏过程的能力和有效性。     

非饱和土水-力本构模型及其隐式积分算法

在已有工作基础上建立了水力-力学耦合的非饱和土本构模型,在硬化方程中考虑饱和度的影响,同时在土水特征曲线中考虑了塑性体变的影响,从而使模型可以反映非饱和土中的毛细现象与土中弹塑性变形现象的耦合行为。采用隐式积分方法,建立了非饱和土耦合模型的数值模型,并推导了得到了水力-力学耦合的非饱和土的一致切线模量。利用该算法编制了本构模型计算的子程序,使其能向外输出切线刚度矩阵,用于有限元计算。为了验证该算法和程序的正确性,用所编制程序对不同路径下的土体行为进行了预测。通过预测结果与试验结果相对比,表明程序预测结果与试验数据相吻合,模型可以较好地模拟土体的水力-力学耦合行为特性。     

基于ABAQUS的颗粒材料亚塑性模型的数值实现及探讨

颗粒材料的亚塑性模型以Jaumann应力率张量及变形率张量描述本构关系,基于Cauchy应力的Jaumann速率及变形率给出了有限元的切线刚度矩阵,由此指出在严格意义上基于ABAQUS二次开发的亚塑性模型数值实现必需借助UEL接口。为了简化程序开发,文中建议了与有限元相结合的两种近似切线模量矩阵,即基于矩阵广义逆一致性切向模量矩阵与近似切向模量矩阵,从而形成了通过UMAT接口实现亚塑性模型的数值方案。由此降低了程序开发的难度,同时也可借助ABAQUS的后处理功能,提高了工作效率。数值算例表明了所开发程序的正确性以及所建议方案的可行性。     

基于C2阶连续函数的广义Hoek-Brown强度准则屈服面与塑性势面棱角圆化方法

Hoek-Brown强度准则的屈服面与塑性势面在棱角处导数无定义,具有数值奇异性,采用圆角进行光滑过渡只能满足一阶导数连续性,而二阶导数仍然无定义,使得棱边上一致切线模量矩阵无法正确计算,导致有限元总体平衡方程组Newton-Raphson隐式迭代二阶收敛性丧失。提出基于C2阶连续函数的广义Hoek-Brown准则屈服面与塑性势面棱角圆化方法,使得棱角处函数曲面二阶连续可导,棱边上一致切线模量矩阵可精确计算。基于ABAQUS数值开发平台,采用FORTRAN语言编制Hoek-Brown准则理想弹塑性UMAT用户子程序,通过数值算例验证所提方法的正确性。     

基于修正Mohr-Coulomb准则的弹塑性本构模型及其数值实施

针对Mohr-Coulomb准则高估岩土体抗拉性能的局限性,建立考虑最大拉应力准则的修正Mohr-Coulomb模型;系统地论述隐式本构积分算法的主要内容,推导相应的一致性刚度矩阵。以ABAQUS软件为平台,采用向后欧拉隐式应力积分算法编制了UMAT本构程序,对单轴拉伸试验和三轴压缩试验进行数值模拟,对比分析ABAQUS自带模型和自编模型的优劣,结果表明编写的修正Mohr-Coulomb模型能够有效地反映岩土介质的抗拉性能,弥补了ABAQUS自带模型的不足。     

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.     

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.     

基于Hardin曲线的土体边界面本构模型在ADINA软件中的实现

为了提出一种适合于岩土地震数值模拟的土体本构模型,基于土体动应力-应变关系的Hardin曲线及其在非等幅往返荷载下的Pyke修正,采用von Mises准则在偏应力平面上构造边界面,以反向加载点和当前应力点的连线在边界面上投影的比例作为硬化参数,推导了塑性硬化模量并给出该边界面本构的具体增量表述。在有限元软件ADINA中通过自定义材料的二次开发实现了该本构模型,并利用动三轴试验对该本构模型进行了验证。数值模拟与试验结果的对比表明,本构模型能如实反映土体的应力-应变关系。针对实际工程场地的地震反应,应用边界面本构模型在ADINA中进行了二维数值模拟,与SHAKE91的计算结果进行了对比,说明了该本构模型应用于岩土地震工程问题的合理性。     

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

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

Implicit integration of elastoplastic relations with reference to soils

The stress solution in plasticity with an associated or a non-associated flow rule is considered. Upon fully implicit integration of the relevant constitutive relations the stress is obtained as the projection of the (fictitious) elastic stress onto the yield surface for given values of the hardening softening variables. This projection is defined, for a general non-associated flow rule, in adjusted complementary elastic energy, which becomes exactly the complementary energy when the flow rule is associated. Isotropic elasticity and mean-stress dependent isotropic yield criteria (pertinent to soil) are considered and the implications of a certain class of non-associated flow rules are evaluated. This class relates to dilatant (or contractant) materials and involves non-associated plastic volume change. The corresponding stress solutions are shown to be strongly influenced by Poisson's ratio and by the dilation angle.     

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.     

Implicit integration algorithm for Hoek-Brown elastic-plastic model

A realistic strength criterion often used to describe the yielding behaviour of a jointed rock mass at a continuum level is the well-known Hoek and Brown criterion. This paper is concerned with a 3-D stress generalization of the Hoek-Brown failure criterion by means of an elliptical functional which leads to a smooth deviatoric trace in the stress space. For its incorporation into a finite element analysis involving plasticity calculations, the formulation of an implicit stress integration algorithm is presented. The key computational methodology alludes to the notion of consistent tangent modulus and implicit return mapping schemes (radial and closest point return) for stress integration in a finite element analysis. Within the context of non-linear elastoplastic analysis, it is found that formulation of such consistent modulus and success into achieving numerical efficiency are closely intertwined. Indeed, the procedure results into accurate and rapid convergence of the displacement finite element scheme during the search for equilibrium. This means that considerable savings in computational time can be achieved for large scale problems. Numerical examples which focus on the Hoek-Brown plasticity model are presented in order to fully appreciate the robustness of the algorithm, and hence the viability of such method to practical problems.     

Finite element formulation and algorithms for unsaturated soils. Part I: Theory

This paper presents a complete finite‐element treatment for unsaturated soil problems. A new formulation of general constitutive equations for unsaturated soils is first presented. In the incremental stress–strain equations, the suction or the pore water pressure is treated as a strain variable instead of a stress variable. The global governing equations are derived in terms of displacement and pore water pressure. The discretized governing equations are then solved using an adaptive time‐stepping scheme which automatically adjusts the time‐step size so that the integration error in the displacements and pore pressures lies close to a specified tolerance. The non‐linearity caused by suction‐dependent plastic yielding, suction‐dependent degree of saturation, and saturation‐dependent permeability is treated in a similar way to the elastoplasticity. An explicit stress integration scheme is used to solve the constitutive stress–strain equations at the Gauss point level. The elastoplastic stiffness matrix in the Euler solution is evaluated using the suction as well as the stresses and hardening parameters at the start of the subincrement, while the elastoplastic matrix in the modified Euler solution is evaluated using the suction at the end of the subincrement. In addition, when applying subincrementation, the same rate is applied to all strain components including the suction. Copyright © 2003 John Wiley & Sons, Ltd.