ABAQUS混凝土塑性损伤模型

1、4.5.2混凝土和其它准脆性材料的塑性损伤模型这部分介绍的是ABAQUS提供分析混凝土和其它准脆性材料的混凝土塑性损伤模型。ABAQUS材料库中也包括分析混凝的其它模型如基于弥散裂纹方法的土本构模型。他们分 别是在 ABAQUS/Sta ndardAn in elastic co nstitutive model for con crete,” Section 4.5.1 ,中的弥散裂纹模型和在 ABAQUS/Explicit, A cracking model for concrete and other brittle materials, ” Section 4.5.3中的脆性开裂模型。

2、混凝土塑性损伤模型主要是用来为分析混凝土结构在循环和动力荷载作用下的提供一个普遍 分析模型。该模型也适用于其它准脆性材料如岩石、砂浆和陶瓷的分析；本节将以混凝土的 力学行为来演示本模型的一些特点。在较低的围压下混凝土表现出脆性性质，主要的失效机 制是拉力作用下的开裂失效和压力作用下的压碎。当围压足够大能够阻止裂纹开裂时脆性就 不太明显了。这种情况下混凝土失效主要表现为微孔洞结构的聚集和坍塌，从而导致混凝土 的宏观力学性质表现得像具有强化性质的延性材料那样。本节介绍的塑性损伤模型并不能有效模拟混凝土在高围压作用下的力学行为。而只能模拟混 凝土和其它脆性材料在与中等围压条件（围压通常小于单轴抗压强

3、度的四分之一或五分之一） 下不可逆损伤有关的一些特性。这些特性在宏观上表现如下：单拉和单压强度不同，单压强度是单拉强度的 10倍甚至更多；受拉软化，而受压在软化前存在强化；在循环荷载（压）下存在刚度恢复；率敏感性，尤其是强度随应变率增加而有较大的提高。概论混凝土非粘性塑性损伤模型的基本要点介绍如下: 应变率分解对率无关的模型附加假定应变率是可以如下分解的：总=占+評，是总应变率，是应变率的弹性部分，是应变率的塑性部分。应力应变关系应力应变关系为下列弹性标量损伤关系：/ = （1 一:一 曰 = D衬；（sr 一 b）其中门 是材料的初始（无损）刚度，I门铠-疋是有损刚度，是刚度退化变量其值在0

4、 （无损）到1 （完全失效）之间变化，与失效机制（开裂和压碎）相关的损伤导致了弹性 刚度的退化。在标量损伤理论框架内，刚度退化是各向同性的，它可由单个标量d来描述。按照传统连续介质力学观点，有效应力可定义如下：cr =：Cauchy应力通过标量退化变量（d ）转化为有效应力（7 = I 1 对如任何一个给定的材料截面，因子代表承力的有效面积占总截面积的比重（总截面积剪除受损面积）。在无损时d=0，有效应力等于cauchy应力。然而，当损伤发生后，有效 应力比cauchy应力更能代表实际情况，因为损伤后截面承力的是有效无损的面积。因此， 可以很方便的用有效应力来建立塑性相关公式。正如后面将要谈论

5、的那样，退化变量的演化 是由一组硬化参数和有效应力控制的：即盘盘春.硬化变量受拉和受压的损伤状态由两个独立的硬化变量和描述，他们分别代表受拉和受压时的等效塑性应变。硬化参数的演化由下式给出（下文将进一步讨论）：1混凝土的微裂纹和压碎由不断增大的硬化变量来描述。这些硬化变量控制着屈服面和弹性刚 度退化。他们也与产生新裂纹面所要消耗的断裂能有密切的关系。屈服函数屈服函数广在有效应力空间内代表一个空间曲面，它决定了失效或损伤的状态。屈服函数，至于本粘性无关的塑性损伤模型其屈服函数的具体形式稍后详细介绍。F（孔評） 0.流动法则根据流动法则，塑性流动由塑性势 G来确定，形式为：式中 为非负的流动因子，

6、塑性势也是定义在有效应力空间里的。 其具体形式稍后介绍。由于 使用的是非相关联流动法则，所以刚度矩阵将会是非对称的。巨二D/ ：仏-护） 可尸（爪尹） 0t評=皿7評）評*（42 -1）小结：总之，塑性损伤本构模型的混凝土弹塑性损伤是在有效应力空间和硬化变量来描述的岳=Dj :仏一 ept) W 0t評=h(7 0)(4.5.2 1) 刃二茫G(卩式中 和F满足Kuhn-Tucker条件叭:Cauchy是由刚度退化变量川厅.評:;和有效应力按下式cr = (1 - d)a.(4.5.2 2)计算得到的。从等式4.5.2-1可以看出，弹塑性关系与刚度退化是非耦合的。式4.5.2-2的优点在于他能

7、方便计算机数值计算。此处总结的非粘性塑性损伤模型可以很轻易地进行拓展就 能考虑粘塑性影响了，只要允许有效应力超出屈服面然后对其归一化就可以了。损伤和刚度退化硬化变量，的演化规律可以很方便的先通过考虑单轴情况在推广到多轴情况来确定(但 实际上从单轴到多轴的推广往往并不容易的，译者认为)dc =或(刖(0 dc rrsjsioL.这一节里面我们约定一是正数，它代表的是单压时的应力值，即*=:匸-，正如在图4.5.2-1 中显示的那样，当从应力-应变曲线的应变软化段卸载时，可以发现卸载的响应是退化了的， 也就是说材料的弹性模量看起来变小了(损伤了)。弹性刚度的损伤在拉压试验中表现是大 不相同的。但在

8、拉压两种情况中，随着塑性变形的增加损伤效果都是越来越明显的。混凝土 的损伤响应由两个独立的单轴损伤变量和，控制，他们是塑性应变、温度和其它行变量的 函数。(1.5.2 5)必=必(犁出J/ (0fA 1). 心=鸣(刖6 齐)* (0 dc 1).图 4.5.2 -,单轴刚度退化变量是等效塑性应变的非减函数，他们的取值范围在0 （无损伤）到1 （完全损伤）之间。如果表示材料的初始弹性刚度，那么在单轴拉压下的应力 -应变关系分别为 叭=（1-阳巴侶-說h% = （1 rfc）o （8 聖）在单轴加载条件下，裂纹是沿着与应力垂直方向发展的。裂纹的成核和扩展就造成了界面有 效承载面积的减小，因此就导

9、致了有效应力的增加。在单轴压是这种承载面积减小的效果还 要稍好一点，因为开始是裂纹基本上是平行于应力方向扩展的，但是当压碎发展到比较厉害 时有效承载面积也将显著地减小。那么有效单轴内聚力 -和形式如下cr-fTf,=有效单轴内聚力决定了屈服（破坏）面的大小。 单轴循环加载在单轴循环加载条件下，刚度退化机制比较复杂，它设计到预先存在裂纹的开闭问题和裂纹 间的相互作用问题。试验观察发现，但循环加载的应力符号变号是反向加载的刚度有所恢复。 这种刚度恢复也称之为“单边效应”它是混凝土循环加载的一个显著特点。特别是当应力有 拉变为压是，效应很明显，这时压应力是的受拉形成的裂纹闭合从而是受压刚度得到恢复。

10、 混凝土塑性损伤模型假定弹性模量按标量减小变量退化E = （1 拧）E也是材料的初始（无损）模量。这个关系式在拉压曲线中都是成立的，刚度减小变量d是应力状态和单轴损伤变量，和，的 函数，在单轴循环条件下 ABAQUS假定下式成立：（1 = 1 S 內.铝 0 1.（4.5.2 G）式中 和 应力状态的函数，弓I入他们是为了反应由于反向加载是刚度恢复效应，他们定义为：= 1 0 1,0 叫 00 if 厅i 0= 1 - we（ I - 厅H ）：Wii） =其中，权系数 和这里假定为材料参数，他们分别控制应力反向是的刚度恢复能力。举例来说，考虑图4.5.2 荷载有拉变成压的情况。假定材料没有初

11、始预损伤，也就是及 ， 那么此时有（1 d） = （1 = （1 （1 Mr（ 1 拉应力）时一 1，正如预计的那样 。反之压应力 C ）时,.0如果:那么：$=；：，材料恢复到受压无损状态C二忌，反之，若-时，材料没有刚度恢复。当在0-1之间取值时表示刚度只能部分恢复。图4.5.2 受压刚度恢复参数 效应的示意图单轴循环加载时的等效塑性演化方程也可以进行推广如下:(4-5.2 7)它在单拉或单压就退化为方程 4.5.2-4的形式。多轴情况有必要把硬化变量的演化规律推广到多轴情况下，在 假定有效塑性应变率可由下式计算得到：Lee and Fen ves ( 1998)的工作基础上,max(4.

12、5.2 8)A. .j式中 和分别是塑性应变率张量的最大和最小主值。叫&)些匚討吵:0 ) 1是拉压应力权重系数，若有效应力张量三个主值全是正时为 1，反之为0。Macauley运算 定义为：。单轴加载情况下方程4.5.2-8退化为单轴定义式4.5.2-4和谀zpf -p/4.5.2-7，因为此时单拉时，单压时。若果对塑性应变率张量的主值进行排序如: H 一丄 二七、爲一那么多轴普通应力条件下 等效塑性盈利率演化可以写成一下矩阵形式：()0弹性刚度退化 混凝土塑性损伤模型认为混凝土的弹性刚度退化时各向同性的，且可以用一个单标量写成如 下形式：= （1 -0 / 1.（4A2-9）式中的刚度退化

13、标量变量d必须与单轴单调加载时的响应一致，同时还要能够反应在循环加 载退化机制带来的复杂性。对普通多轴加载情况ABAQUS假定，（1 - J） = （1 -旳心）（1 一 餅讣 0 叭、聲 1,（4.5.2 10）形式上与单轴相同，只是现在通过应力权重系数将它推广到多轴情况了：= 1 -* 口厂（疗）：0 冲 1.%=1 h-c（ 1 f厅I）：0 tre 1.显然，很容易验证方程4.5.2-10的标量退化式与单轴加载时是一致的很多准脆性材料（混凝土）的试验表明，当拉应力换到压应力时由于裂纹闭合受压刚度将会 恢复。但是另一方面，当受压是的微裂纹压碎时，由受压换到受拉时的受拉刚度将不会恢复。 鉴

14、于此，ABAQUS默认条件下，假定八-二及即只有受压刚度恢复而没有受拉刚度恢复。图4.5.2-3就是默认条件下的一个应力循环的曲线图图4.523默认条件下（宀-门，.）单轴应力循环曲线图（拉-压-拉）屈服条件本模型的屈服条件基于Lubliner等人（1989）建议的屈服函数，它综合了 Lee and Fenves （1998）的修正以考虑拉压不同时强度的不同演化规律。用有效应力表达时的屈服函数为：）（帀也叨+汙（評）（抒吋?;）一 了（一厅他期）厅*（詐）W（4A2-11）式中和是无量纲材料参数时有效静水压力，是Mises等效应力，恳一.鬥卜沐是有效应力张量的偏量部分，而是厂的代数最大主值，函

15、数形式如下孔尹）=e丄-打）口 _-（1 +“）：西（齐）亠式中和分别为有效拉压内聚力。在双轴受压时，卓賈g;了方程4.5.2-11就退化为Drucker-Prage屈服条件，材料系数可由单轴受压强度和双轴受压强度比值给出:47M） 一 cO（1 = 2疗to 仃cO般材性试验给出的单双受压强度比值在1.10 -1.16之间，那么取值在0.08 -0.12之间（Lubliner et al., 1989 ）系数只在三维受压时才出现在公式中，它可以通过比较沿拉压子午线的强度比值得到。根据定义拉子午线是满足主应力空间中：九小一耳 U - a、的轨迹线，而压子午线是满足的轨迹线。其中，和是应力主值。

16、显然易求得，沿拉压子午线其表达式为： UJ.=雋一二，= V/ -o当扎赵盅时，响应的屈服准则为：(TM)(石十 1 )(7 (十 3tk )/ = ( 1 亍十i丿q(:+3口)卩=(i 一 u)疗厂(CM)令，为静水压力，那么就有；。事实上大多数试验也并没有” _ 3(1 - Ac)证明.是变化的，因此就可求出。对于混凝土来说一般取- I，那么1 = 3”。当时，沿拉压子午线的屈服函数就简化为：(W $+1)Q (+ 3ti)j) (1 ct)听 (TM)+ 1丿 十3u)0= (1 -芥 (CM)K = 1+；同理令灯汕，那么 y在偏片面上典型的屈服面见图4.5.2-4，图4.5.2-5

17、是平面应力时的屈服面。图4.5.2-4 :对应于不同的值在片平面内的屈服面。(C.M.)|-s3图4.5.2-5平面应力时的屈服面uniaxial tensionbiaxial compression1 _(q- 3ap + (je?2) =uniaxial compressionbiaxial tension1 _B_4(q - 3a p + pa.11-a流动法则本模型取的是非关联流动法则：百=人 uc塑性势G取为Drucker-Prager双曲函数的形式G (“yj2 + 沪戶tan 映(4.5,2 12)式中 是p -q面内高围压时的膨胀角，一是单轴抗拉强度，是势函数偏心率，它描述势函

18、 数向其渐近线逼近的速度(当偏心率趋于零时，流动势函数趋于直线)。流动势函数的连续 光滑性保证了流动方向的唯一性。当围压很高时流动势函数渐近于线性Drucker-Prager势函数，且与静水轴的交角是 90 度。在 Models for granular or polymer behavior, ” Section 442 , 中对这个势函数有详细的讨论。因为采用了非关联流动法则，刚度矩阵将会出现非对称。粘塑性归一化在隐式分析程序里，当材料模型出现软化或刚度退化是往往难收敛。有些收敛困难可以通过 对模型的粘塑性归一化来解决。本模型可用粘塑性归一化，因而就允许有效应力超出屈服面。根据Duvaut

19、-Lions归一化粘塑性应变定义为：逊=丄(占一理)-(4.5.2 13)式中 是粘性参数表征粘塑性系统随时间的松弛，是非粘性backbone model的计算塑性应变。同理，在粘塑性体系里，粘性刚度退化变量为。 1 .dr =(川一),(4,5,2 14)d是非粘性backbone model的刚度退化变量，那么粘塑性模型的应力-应变关系就为：(r= (1 -心)Dj : (e -曙).(4.5,2-15)当粘塑性系统的解就趋近于非粘性情况。其中t代表时间。当用较小的粘性系数时就可以对其进行粘塑性归一化，通常能够改善模型在软化段内的收敛速度，而不会影响计算 结果。模型数值计算模型采用后退欧拉

20、法进行计算，该方法在 ABAQUS里对塑性计算用的很多。在平衡迭代中 米用了与积分算子相容的雅克比矩阵。办公室卫生管理制度一、主要内容与适用范围1 本制度规定了办公室卫生管理的工作内容和要求及检查与考核。2 此管理制度适用于本公司所有办公室卫生的管理二、定义1 公共区域：包括办公室走道、会议室、卫生间，每天由行政文员进行清扫；2个人区域：包括个人办公桌及办公区域由各部门工作人员每天自行清扫。1. 公共区域环境卫生应做到以下几点：1) 保持公共区域及个人区域地面干净清洁、无污物、污水、浮土，无死角。2) 保持门窗干净、无尘土、玻璃清洁、透明 3) 保持墙壁清洁，表面无灰尘、污迹。4)保持挂件、画框及其他装饰品表面干净整洁。5) 保持卫生间、洗手池内无污垢，经常保持清洁，毛巾放在固定(或隐蔽)的地方。