Lemaitre等向硬化弹塑性损伤耦合本构模型积分算法及程序实现EI北大核心CSCD

涉及复杂材料弹塑性损伤问题数值计算研究时,不仅需要选择恰当预测损伤和破坏的本构模型,还需要有效和稳健的本构积分算法。首先,阐述了在热力学和连续介质力学框架下建立弹塑性损伤本构模型的基本步骤;其次,基于Lemaitre等向硬化弹塑性损伤耦合本构模型、相应的本构积分算法-完全隐式返回映射算法(Fully Return Mapping Algorithm)和一致切线模量,采用C++语言在Visual 6.0环境下编制有限元本构求解程序,在塑性损伤修正步中求解返回映射方程时,选取一种简单的形式,只需迭代求解一个标量非线性方程,计算效率较高。最后,通过缺口圆棒数值算例初步验证了程序的正确性,并编制接口程序对计算结果进行可视化。研究结果表明积分算法的有效性及程序的正确性,Lemaitre等向硬化弹塑性损伤耦合本构模型能够较好地模拟韧性材料的破坏发展过程,可以求解类似的有限元边界值问题,为考虑损伤特性的韧性材料结构研究和设计奠定基础。     

基于Hoek-Brown准则的岩体弹塑性损伤模型及其应力回映算法研究

为了克服一般弹塑性损伤模型不能反映岩体结构、岩块强度、应力状态的影响以及非线性破坏特征等问题,该文基于广义的Hoek-Brown(HB)屈服准则,考虑损伤引起的刚度退化和塑性导致的流动两种破坏机制的耦合作用,同时引入修正有效应力原理来考虑孔隙水压力的作用,建立了岩体弹塑性损伤本构模型,给出了损伤变量定义及演化方程。针对该模型在数值求解过程中存在的奇异点问题,从主应力空间推导了弹塑性损伤模型的完全隐式返回映射求解算法,包括弹性预测、塑性修正和损伤修正三个步骤。通过ABAQUS软件的用户子程序接口Umat,实现了弹塑性损伤模型的数值求解过程。采用单轴、三轴压缩试验和隧道算例对模型算法进行验证和分析,结果表明,所建立的HB损伤本构模型能够很好地描述岩体材料的力学特性,在实际岩体工程的损伤模拟中效果令人满意,计算结果对工程有指导意义。     

完全隐式返回映射算法对岩土地基问题的求解

针对岩土工程中的复杂力学问题,在弹塑性力学理论框架和非线性有限元理论基础上,采用非关联等向硬化Drucker-Prager模型的完全隐式积分算法—返回映射算法(Return Mapping Algorithm)编制了有限元求解程序。该算法可以避免预测应力漂移屈服面的现象,对准静态变形条件下的本构方程可以获得准确解,在迭代中使用Newton-Raphson法获得近似平方的收敛速率,具有较高的精确性和稳定性。对岩土工程中的地基问题进行求解,计算得出位移、应力等结果,模拟了塑性区随载荷步增加的演化过程,对地基极限承载力进行了解析解和数值解的对比。结果表明了算法的优越性、程序的正确性和实用性。     

基于弹塑性损伤本构模型的复合材料层合板破坏荷载预测

在连续损伤力学和塑性力学框架内,建立一个同时考虑塑性效应和损伤累积导致材料属性退化的复合材料弹塑性损伤本构模型。基于最近点投影回映算法,开发本构模型的应变驱动隐式积分算法以更新应力及与解答相关的状态变量,并推导与所开发算法相应的数值一致性切线刚度矩阵,保证有限元分析采用NewtonRaphson迭代法解答非线性问题的计算效率。采用断裂带模型对已开发的本构模型软化段进行规则化,以减轻有限元分析结果的网格相关性问题。对损伤变量进行粘滞规则化,并推导出相应的粘滞规则化数值一致性切线刚度张量,解决了在有限元隐式计算程序中采用含应变软化段本构关系的数值分析由于计算困难而提前终止的问题。开发包含数值积分算法的用户材料子程序UMAT,并嵌于有限元程序Abaqus v6.14中。通过对力学行为展现显著塑性效应的AS4/3501-6V型开口复合材料层合板的渐进失效分析,验证本文提出的材料本构模型的有效性。结果显示,预测结果与已报道的试验结果吻合良好,并且预测精度高于其他已有弹性损伤模型。表明已建立的弹塑性损伤本构模型能够准确预测力学行为,展现显著塑性效应的复合材料层合板的破坏荷载,为其构件和结构设计提供一种有效的分析方法。     

弹塑性接触有限元混合法及其在齿轮传动中的应用

弹塑性接触问题属于一种表面非线性和材料非线性耦合的双非线性问题。本文采用有限元混合法及Newton-Raphson迭代求解。通过有限元离散,建立相互接触物体的接触边界连续方程,求出接触边界上的力学参量,再耦合推导出弹塑性接触有限元列式,编制了相应的Fortran程序,并在通过实例考核以后分析了重点工程重载齿轮的弹塑性行为。数值计算结果表明,本文所提供的方法和程序是通用而有效的。计算结果令人满意。     

RPC200塑性与各向异性损伤耦合本构建模研究

针对活性粉末混凝土 (RPC200:Reactive Powder Concrete 200MPa)的非线性行为,将连续损伤力学和塑性力学结合起来,建立了RPC200的弹塑性各向异性损伤耦合本构。在有效应力空间内,采用非关联流动法则和非均匀等向强化法则建立了RPC200的塑性本构。基于损伤能释放率建立了损伤准则,并由正交法则给出损伤演化法则,同时采用两个不同的损伤硬化法则来描述拉、压性能不同的损伤硬化。证明了该文所建模型与热动力学基本原理的相容性。最后,给出了RPC弹塑性各向异性损伤耦合本构模型数值计算流程,并以简单例证证明了所建模型的预测能力。     

饱和软粘土地基的损伤模型与震陷计算

基于各向同性弹塑性损伤和Prevost模型的基本理论,把弹塑性等向硬化、运动硬化和各向同性损伤结合起来,推导了循环荷载作用下不排水饱和软粘土弹塑性动力损伤本构模型.考虑到地震作用下土体应力的不规则性,对循环三轴试验获得的粘土残余应变经验计算公式进行了修正,最后将该残余应变引入到损伤演化方程中.通过对地基的震陷计算并与不考虑损伤的模型进行对比,结果表明,该模型能合理地考虑屈服面内应力循环对地基残余塑性应变的贡献.     

有限变形下材料的弹塑性损伤行为及结构响应分析

该文基于一种弹塑性损伤本构模型,分析了材料在有限变形下的力学行为及结构力学响应。将模型建立的本构方程进行线性化处理,嵌入到有限元计算软件ABAQUS的子程序中,实现模型在结构计算中的应用。针对铝合金材料(6061-T6),建立了两种不同的结构模型进行数值模拟。模拟结果与实验结果对比表明,模型能较好地描述材料的变形过程及结构的力学行为。最后基于模型的理论基础,分析了材料在超出弹性范围时的复杂变形行为。     

单向纤维增强聚合物复合材料压缩渐进破坏单向纤维增强聚合物复合材料压缩渐进破坏

复合材料结构在承压时破坏如何演化,是其强度破坏分析的基础和核心任务。本文提出了基于连续介质损伤力学（CDM）的单向纤维增强聚合物复合材料压缩破坏渐进损伤分析（PDA）模型。建模中考虑了材料非线性行为、失效判断及损伤演化中材料性能退化等基本问题,分别对应于拉压不对称弹塑性本构关系、Puck准则、LaRC05准则及考虑破坏面方向的刚度退化方法。该模型通过用户材料子程序接口VUMAT引入到有限元软件ABAQUS中实现了有限元求解。对文献中提供的纵向、横向及偏轴压缩案例进行了数值计算并与试验数据对比。数值分析结果与试验数据吻合较好,证明了该方法的合理性和有效性,对开展多向层合板压缩破坏分析富有参考价值。      

纤维增强复合材料本构模型研究进展

本文对纤维增强复合材料的本构模型最新研究成果进行简要阐述和归纳,并结合实例进行了相关数值模拟计算。文中阐述的复合材料为纤维增强树脂复合材料、纤维增强金属基复合材料和纤维编织复合材料。目前,关于纤维增强金属基复合材料的研究较少,大多数研究集中在纤维增强树脂复合材料和纤维编织复合材料。许多学者在经典弹性本构模型和连续介质力学基础上,推导出了一些弹塑性、粘弹性和考虑损伤的非线性本构模型。通过实例计算结果表明,考虑损伤导致刚度退化的复合材料弹塑性本构模型计算得到的拉伸应力应变曲线与测试结果基本一致。本文通过对复合材料本构模型的最新研究成果的归纳和数值模拟实例计算,为后续研究工作提供借鉴,推动其在实际工程中的应用。     

A combined elastoplastic damage model for progressive failure analysis of composite materials and structures

The paper is concerned with the development and verification of a combined elastoplastic damage model for the progressive failure analysis of composite materials and structures. The model accounts for the irreversible strains caused by plasticity effects and material properties degradation due to the damage initiation and development. The strain-driven implicit integration procedure is developed using equations of continuum damage mechanics, plasticity theory and includes the return mapping algorithm. A tangent operator consistent with the integration procedure is derived to ensure a computational efficiency of the Newton–Raphson method in the finite element analysis. The algorithm is implemented in Abaqus as a user-defined subroutine. The efficiency of the constitutive model and computational procedure is demonstrated using the analysis of the progressive failure of composite laminates containing through holes and subjected to in-plane uniaxial tensile loading. It has been shown that the predicted results agree well with the experimental data reported in the literature.     

Finite element implementation of non‐linear elastoplastic constitutive laws using local and global explicit algorithms with automatic error control

Implicit stress integration algorithms have been demonstrated to provide a robust formulation for finite element analyses in computational mechanics, but are difficult and impractical to apply to increasingly complex non‐linear constitutive laws. This paper discusses the performance of fully explicit local and global algorithms with automatic error control used to integrate general non‐linear constitutive laws into a non‐linear finite element computer code. The local explicit stress integration procedure falls under the category of return mapping algorithm with standard operator split and does not require the determination of initial yield or the use of any form of stress adjustment to prevent drift from the yield surface. The global equations are solved using an explicit load stepping with automatic error control algorithm in which the convergence criterion is used to compute automatically the coarse load increment size. The proposed numerical procedure is illustrated here through the implementation of a set of elastoplastic constitutive relations including isotropic and kinematic hardening as well as small strain hysteretic non‐linearity. A series of numerical simulations confirm the robustness, accuracy and efficiency of the algorithms at the local and global level. Published in 2001 by John Wiley & Sons, Ltd.     

Combined continuum damage‐embedded discontinuity model for explicit dynamic fracture analyses of quasi‐brittle materials

In this paper, a novel constitutive model combining continuum damage with embedded discontinuity is developed for explicit dynamic analyses of quasi‐brittle failure phenomena. The model is capable of describing the rate‐dependent behavior in dynamics and the three phases in failure of quasi‐brittle materials. The first phase is always linear elastic, followed by the second phase corresponding to fracture‐process zone creation, represented with rate‐dependent continuum damage with isotropic hardening formulated by utilizing consistency approach. The third and final phase, involving nonlinear softening, is formulated by using an embedded displacement discontinuity model with constant displacement jumps both in normal and tangential directions. The proposed model is capable of describing the rate‐dependent ductile to brittle transition typical of cohesive materials (e.g., rocks and ice). The model is implemented in the finite element setting by using the CST elements. The displacement jump vector is solved for implicitly at the local (finite element) level along with a viscoplastic return mapping algorithm, whereas the global equations of motion are solved with explicit time‐stepping scheme. The model performance is illustrated by several numerical simulations, including both material point and structural tests. The final validation example concerns the dynamic Brazilian disc test on rock material under plane stress assumption. Copyright © 2014 John Wiley & Sons, Ltd.     

On the numerical integration of a class of pressure-dependent plasticity models

An unconditionally stable algorithm for the numerical integration of elastoplastic pressure-dependent constitutive relations is analysed in detail in this paper. The application of the method to plane stress problems, in which the out-of-plane strain component is not defined kinematically, is discussed. The tangent moduli resulting from this integration algorithm are obtained by consistent linearization of the elastoplastic constitutive equations. The algorithm is applied to Gurson's constitutive model, some one-dimensional problems are solved, and comparisons with exact solutions are made. The paper closes with a numerical study of the necking of an axi-symmetric specimen using Gurson's plasticity model to describe the constitutive behaviour of the material.     

An analysis of a new class of integration algorithms for elastoplastic constitutive relations

An accuracy analysis of a new class of integration algorithms for finite deformation elastoplastic constitutive relations recently proposed by the authors, is carried out in this paper. For simplicity, attention is confined to infinitesimal deformations. The integration rules under consideration fall within the category of return mapping algorithms and follow in a straightforward manner from the theory of operator splitting applied to elastoplastic constitutive relations. General rate-independent and rate-dependent behaviour, with plastic hardening or softening, associated or non-associated flow rules and nonlinear elastic response can be efficiently treated within the present framework. Isoerror maps are presented which demonstrate the good accuracy properties of the algorithm even for strain increments much larger than the characteristic strains at yielding.     

Large strain constitutive modelling and computation for isotropic,creep elastoplastic damage solids

A three-dimensional fully coupled creep elastoplastic damage model at finite strain for isotropic non-linear material is developed. The model is based on the thermodynamics of an irreversible process and the internal state variable theory. A hyperelastic form of stress–strain constitutive relation in conjunction with the multiplicative decomposition of the deformation gradient into elastic and inelastic parts is employed. The pressure-dependent plasticity with strain hardening and the damage model with two damage internal variables are particularly considered. The rounding of stress–strain curves appearing in cycling loading is reproduced by introduction of the creep mechanism into the model. A numerical integration procedure for the coupled constitutive equations with three hierarchical phases is proposed. A consistent tangent matrix with consideration of the fully coupled effects at finite strain is derived. Numerical examples are tested to demonstrate the capability and performance of the present model at large strain.     

On the numerical prediction of the ductile fracture in metal forming

In this paper, fully coupled constitutive equations accounting for combined isotropic and kinematic hardening as well as the isotropic ductile damage are implemented into the general purpose finite element code for metal forming simulation. The associated numerical aspects concerning both the local integration of the coupled constitutive equations as well as the global (equilibrium) integration schemes are presented. Various 2D and 3D examples are given in order to show the capability of the proposed numerical methodology to predict the ductile fracture initiation and growth during metal forming processes.