ABAQUS混凝土塑性损伤模型概述-孙庆昭.pdf

Chongqing Architecture0引言钢筋混凝土结构在土木工程中的应用非常广泛。对于常规的建筑结构，多止步于弹性阶段的结构设计，并没有过多地考察其非线性特性。但是随着建筑结构设计要求的提高、经济投资的加大，越来越多的超限结构需要进行弹性设计以外的弹塑性时程分析验算，以获得结构体系在地震往复作用下的反应。近年来，随着有限元理论和电子计算机的飞速发展，大大促进了钢筋混凝土结构的非线性分析。ABAQUS作为具有强大非线性计算能力的通用有限元分析软件在土木工程结构计算中得到了广泛应用。在ABAQUS中，为混凝土材料定义了一种材料模型塑性损伤模型，它可以模拟混凝土材料的拉裂和压碎等力学现象。ABAQUS中的塑性损伤模型是在Lubliner1和Lee and Fenves2模型的基础上建立的。该模型使用各向同性损伤弹性结合各向同性拉伸和压缩塑性的模式来表示混凝土的非弹性行为，是一个基于塑性的连续介质损伤模型。同时，它基于各向相同破坏的假设，可用于单向加载、循环加载及动态加载等情况，同时考虑了由于拉压塑性应变导致的弹性刚度的退化以及循环荷载作用下的刚度恢复，具有较好的收敛性3。1 ABAQUS塑性损伤模型理论ABAQUS要求用户填入的是混凝土单轴的拉伸和压缩应力-应变曲线上的一些特征点数据，它可以直接被调用而用于ABAQUS中的一维梁单元，对于二维壳单元和三维实体单元，程序会自动运用广义胡克定律（弹性）和塑性位势理论（弹塑性）确定二维和三维应力状态下的空间应力-应变关系。混凝土材料单轴拉伸和压缩应力-应变曲线（如图1）显示：应力-应变曲线存在下降段，应力超过弹性范围以后卸载时存在刚度退化现象，这与混凝土材料内部存在微小裂缝（损伤）有关。混凝土塑性损伤模型的受拉应力-应变数据在ABAQUS中以t-0/ckot的形式输入。其中0/ckot称为开裂应变，其被定义为总应变减去材料无损伤的弹性应变，0/ckot=t-elot，elot=tE0。当拉应力超过弹性极限t0后，将按(1-dt)E0的退化刚度卸载，dt称为受拉损伤因子，它是受拉等效塑性应变0/plt的非减函数，其取值范围在0（无损伤）到1（完全损伤）之间。受拉损伤数据以dt-0/ckt的形式输入到程序中，程序自动按下式将开裂应变转化为等效塑性应变0/plt：0/plt=0/ckt-dt1-dttE0（1）混凝土塑性损伤模型的受压应力-应变数据在ABAQUS中以c-0/cm的形式输入。其中0/mc称为受压非弹性应变，其被定ABAQUS混凝土塑性损伤模型概述孙庆昭（中机中联工程有限公司，重庆400039）摘要：该文对大型非线性有限元分析软件ABAQUS中的混凝土塑性损伤模型进行了详细的介绍。着重讨论了混凝土塑性损伤模型的基本理论和损伤因子基于混凝土结构设计规范（GB50010-2010）4附录C应力-应变曲线的计算方法。同时也对ABAQUS中塑性损伤材料模型的滞回规则、屈服准则和流动法则作了简要介绍。关键词：塑性损伤；ABAQUS；损伤因子；滞回规则；屈服准则；流动法则；混凝土中图分类号：TU313文献标识码：A文章编号：1671-9107（2014）11-0070-03收稿日期：2014-10-23作者简介：孙庆昭（1989-），男，河南虞城人，本科，主要从事建筑设计工作。oooo2014.NO.1170建筑结构doi：10.3969j.issn.1671-9107.2014.11.070图1混凝土单轴拉伸和压缩应力-应变曲线（上图：拉伸；下图：压缩）ttoccucoEoEoEo(1-dt)Eo(1-dc)EopltcktoltoltplclncolcolocEoc0ooAn Overview of Plastic Damage Model of ABAQUS ConcreteAbstract:Concrete plastic damage model in nonlinear finite element analysis software ABAQUS is introduced in detail.The basic theories of concreteplastic damage model and the calculation methods of the damage factors based on stress-strain curve in appendix C of Code for Concrete Structures De-sign（GB50010-2010）4are discussed in particular.Also,the hysteresis rules，yield criterion and flow rule for plastic damage material model inABAQUS are introduced for reference.Key words:plastic damage;ABAQUS;damage factor;hysteresis rules;yield criterion;flow rule;concrete重庆建筑2014.N0.义为总应变减去材料无损伤的弹性应变0/mc=c-eloc，eloc=cE0。当压应力超过弹性极限c o后，将按(1-dc)E0的退化刚度卸载，dc称为受压损伤因子，它是受压等效塑性应变0/plc的非减函数，其取值范围在0（无损伤）到1（完全损伤）之间。受压损伤数据以dc-0/mc的形式输入到程序中，程序自动按下式将开裂应变转化为等效塑性应变0/plc：0/plc=0/mc-dc1-dccE0（2）2损伤因子的取值ABAQUS需要用户输入的材料信息包括混凝土材料的单轴的受压、受拉应力-应变关系曲线以及受压、受拉时的损伤因子-非弹性应变曲线（仅输入曲线上的一组离散点，程序会自动根据所输入的离散点进行拟合，因此离散点的选择至关重要）。材料本构曲线和损伤因子曲线一般通过实验得到，当缺少实验数据时，可以结合混凝土结构设计规范（GB50010-2010）4附录C提供的混凝土应力-应变曲线利用能量等价原理计算得到所需参数。规范4给出的混凝土单轴应力-应变曲线如图2所示，其中ft,r、fc,r分别为混凝土单轴抗拉、抗压强度代表值；t,r、c,r分别为与单轴抗拉强度代表值和单轴抗压强度代表值对应的混凝土峰值应变。混凝土单轴受拉的应力应变曲线公式如下：=(1-dt)Ec（3）当x燮1时，dt=1-t(1.2-0.2x5)（4）当x1时，dt=1-t1-t(x-1)1.7+x（5）其中，x=t,r（6）t=ft,rEct,r（7）t为混凝土单轴受拉应力-应变曲线下降段的参数值，按规范4表C.2.3取用。令y=ft,r，结合式（3）至式（7），混凝土单轴受拉应力-应变曲线可以用如下公式表示：当x燮1时，y=1.2x-0.2x6（8）当x1时，y=xt(x-1)1.7+x（9）Sidiroff的能量等价原理认为，应力作用在受损材料产生的弹性余能与作用在无损材料产生的弹性余能在形式上是相同的5，只要将应力改为等效应力，或将弹性模量改为损伤时的等效弹性模量即可。无损伤材料弹性余能：We0=22Ed（10）等效有损伤材料弹性余能：Wed=22E0（11）=1-d为有效应力，于是得Ed=E0(1-d)2，进一步可得：=E0(1-d)2（12）将式（12）归一化得：y=(1-d)2tx（13）将式（8）、式（9）分别与式（13）对应，即可解出：当x燮1时，d=1-t(1.2-0.2x5)姨（14）当x1时，d=1-tt(x-1)1.7+x姨（15）同理，可得单轴受压损伤因子计算公式如下：当x燮1时，d=1-cnn-1+xn姨（16）当x1时，d=1-cc(x-1)2+x姨（17）其中，x=c,r（18）c=fc,rEcc,r（19）n=Ecc,rEcc,r-fc,r（20）3滞回规则ABAQUS中混凝土塑性损伤模型的滞回准则由损伤因子d和刚度恢复系数共同决定。刚度恢复系数分为受拉刚度恢复系数t和受压刚度恢复系数c，分别表示混凝土应力-应变曲线从受压区过渡到受拉区和从受拉区过渡到受压区时，弹性模量的变化程度。c=1表示混凝土从受拉区进入受压区时弹性模量可完全恢复至上一次受压卸载时的弹性模量，c=0表示混凝土弹性模量不能恢复，如图3所示。t=1表示混凝土从受压ooo1 1第13卷 总第133期71ABAQUS混凝土塑性损伤模型概述建筑结构图2混凝土单轴应力-应变曲线c.rfc.r0.5fc.sc.ul.r0fl.roooChongqing Architecture责任编辑：孙苏区进入受拉区时弹性模量可完全恢复至上一次受拉卸载时的弹性模量，t=0表示混凝土弹性模量不能恢复6。很多准脆性材料（混凝土）的试验表明，当从拉应力变到压应力时，由于裂纹闭合，受压刚度将会恢复。但是另一方面，当受压的微裂纹压碎时，由受压变为受拉时的受拉刚度将不会恢复。鉴于此，ABAQUS默认条件下，假定t=0及c=1，即只有受压刚度恢复而没有受拉刚度恢复3。图4为默认条件下的一个应力循环的曲线图。4屈服准则在复杂应力状态下，物体内某一点开始产生塑性变形时，应力必须满足一定的条件，它就是复杂应力状态下的屈服准则。一般来说，它应是6个应力分量的函数，可表示为：F(x,y,z,xy,yz,zx)=C（21）式中，C为与材料有关的常数，F为屈服函数7。ABAQUS中混凝土塑性损伤模型采用的屈服准则是基于Lubliner等人（1989）1建议的屈服函数，它综合了Lee andFenves（1998）2的修正以考虑拉压不同时强度的不同演化规律。用有效应力表达时的屈服函数为：F=11-(q-3p+0/plo)-)-c(0/plc)=0（22）式中和是无量纲材料参数，p=-13I是有效静水压力，q=32(SS姨)是Mises等效应力，S=+pI是有 效应 力 张量 的偏 量部 分，而max是的代数最大主值，函数(0/plo)=c(0/plc)t(0/plt)(1-)-(1+)，式中t和c分别为有效拉压内聚力。5流动法则材料在初始弹性范围内，应力与应变之间存在着一一对应的关系，即广义胡克定律。进入塑性状态以后，一般来说，不再存在着应力与应变之间的一一对应关系，只能建立应力增量与应变增量之间的关系。这种用增量形式表示的材料本构关系，称为增量理论或流动法则7。流动法则是关于塑性应变增量方向的假定准则，经典塑性位势理论假定塑性应变增量方向满足正交流动规则，即塑性应变增量&pl的方向与塑性势函数的梯度方向或塑性势面的外法线方向一致。其数学表达式为：&pl=&坠Q()坠（23）式中&是一个非负的标量因子，表示塑性应变增量的大小；Q为塑性势函数。如果假定塑性势函数等于加载函数或屈服函数，即Q=F，这样塑性势函数与加载函数或屈服函数相关联，称为相关联的流动法则；如果QF，则&pl的方向与加载面或屈服面不正交，但仍与塑性势面正交，这时称为不相关联的流动法则，这样塑性势函数的选取就非唯一。ABAQUS中混凝土塑性损伤模型采用的是非关联流动法则，其塑性势的数学表达式取为Drucker-Prager双曲面函数的形式：G=缀t0tan)2+q2姨-ptan（24）式中为膨胀角，t0是单轴抗拉强度，缀是势函数偏心率。6结语本文详细介绍了ABAQUS混凝土塑性损伤模型的基本理论，基于能量等效原理推导了与现行混凝土结构设计规范（GB50010-2010）4附录C提供的混凝土受压和受拉应力-应变关系曲线对应的损伤因子计算公式。并对复杂应力状态下的弹塑性有限元分析中较难理解的滞回规则、屈服准则以及流动法则进行了简单阐述，有助于工程界对ABAQUS软件和弹塑性有限元分析理论的理解与掌握。参考文献：1Lubliner J，Oliver J，Oller S，et al.A Plastic-damagemodel for concreteJ.International Journal