《ABAQUS混凝土塑性损伤模型》由会员分享，可在线阅读，更多相关《ABAQUS混凝土塑性损伤模型（11页珍藏版）》请在金锄头文库上搜索。
1、4.52混凝土和其它准脆性材料的塑性损伤模型这部分介绍的是ABAQUS提供分析混凝土和其它准脆性材料的混凝土塑性损伤模型。 ABAQUS材料库中也包括分析混凝的其它模型如基于弥散裂纹方法的土本构模型。他们分 另 H是在 ABAQUS/Standard “An in elastic con stitutive model for con crete,” Secti on 4.5.1,中 的弥散裂纹模型和在 ABAQUS/Explicit, “A cracking model for concrete and other brittle materials,” Section 4.5.3中的脆性开
2、裂模型。混凝土塑性损伤模型主要是用来为分析混凝土结构在循环和动力荷载作用下的提供一个普遍 分析模型。该模型也适用于其它准脆性材料如岩石、砂浆和陶瓷的分析；本节将以混凝土的 力学行为来演示本模型的一些特点。在较低的围压下混凝土表现出脆性性质，主要的失效机 制是拉力作用下的开裂失效和压力作用下的压碎。当围压足够大能够阻止裂纹开裂时脆性就 不太明显了。这种情况下混凝土失效主要表现为微孔洞结构的聚集和坍塌，从而导致混凝土 的宏观力学性质表现得像具有强化性质的延性材料那样。本节介绍的塑性损伤模型并不能有效模拟混凝土在高围压作用下的力学行为。而只能模拟混 凝土和其它脆性材料在与中等围压条件（围压通常小于单
3、轴抗压强度的四分之一或五分之一） 下不可逆损伤有关的一些特性。这些特性在宏观上表现如下：单拉和单压强度不同，单压强度是单拉强度的10倍甚至更多；受拉软化，而受压在软化前存在强化；在循环荷载（压）下存在刚度恢复；率敏感性，尤其是强度随应变率增加而有较大的提高。概论 混凝土非粘性塑性损伤模型的基本要点介绍如下：应变率分解对率无关的模型附加假定应变率是可以如下分解的：合二亠尹:是总应变率，-是应变率的弹性部分，：是应变率的塑性部分。应力应变关系应力应变关系为下列弹性标量损伤关系：护一 （1 仲:W -日）一 D创:（ -护）.其中门是材料的初始（无损）刚度，是有损刚度，是刚度退化变量其值在0 （无损
4、）到1 （完全失效）之间变化，与失效机制（开裂和压碎）相关的损伤导致了弹性 刚度的退化。在标量损伤理论框架内，刚度退化是各向同性的，它可由单个标量d来描述。 按照传统连续介质力学观点，有效应力可定义如下：D；/ :（-Cauchy应力通过标量退化变量（d）转化为有效应力挣=机对如任何一个给定的材料截面，因子I代表承力的有效面积占总截面积的比重（总截面积剪除受损面积）。在无损时d=0，有效应力等于cauchy应力。然而，当损伤发生后，有效 应力比cauchy应力更能代表实际情况，因为损伤后截面承力的是有效无损的面积。因此， 可以很方便的用有效应力来建立塑性相关公式。正如后面将要谈论的那样，退化变
5、量的演化 是由一组硬化参数和有效应力控制的：即d=d（7 出）.硬化变量受拉和受压的损伤状态由两个独立的硬化变量和描述，他们分别代表受拉和受压时的等 效塑性应变。硬化参数的演化由下式给出（下文将进一步讨论）：混凝土的微裂纹和压碎由不断增大的硬化变量来描述。这些硬化变量控制着屈服面和弹性刚 度退化。他们也与产生新裂纹面所要消耗的断裂能有密切的关系。屈服函数屈服函数厂在有效应力空间内代表一个空间曲面，它决定了失效或损伤的状态。屈服函数，至于本粘性无关的塑性损伤模型其屈服函数的具体形式稍后详细介绍。流动法则八啟型根据流动法则，塑性流动由塑性势G来确定，形式为：式中;为非负的流动因子，塑性势也是定义在
6、有效应力空间里的。其具体形式稍后介绍。由于 使用的是非相关联流动法则，所以刚度矩阵将会是非对称的。評評）旳（4.5.2. 1）讥、DG（fr）如小结：总之，塑性损伤本构模型的混凝土弹塑性损伤是在有效应力空间和硬化变量来描述的Eg 評 衍|F伉評兰0.泸二h（讥評）言叭（15.2- 1）说、dG（/T）-式中;.和F满足Kuhn-Tucker条件：二，11 Cauchy是由刚度退化变量和有效应力按下式w = （1 尙乳（4,5.2 2计算得到的。从等式4.5.2-1可以看出，弹塑性关系与刚度退化是非耦合的。式4.5.2-2的优 点在于他能方便计算机数值计算。此处总结的非粘性塑性损伤模型可以很轻易
7、地进行拓展就 能考虑粘塑性影响了，只要允许有效应力超出屈服面然后对其归一化就可以了。损伤和刚度退化硬化变量，的演化规律可以很方便的先通过考虑单轴情况在推广到多轴情况来确定（但 实际上从单轴到多轴的推广往往并不容易的，译者认为）(4.5.2-5) =皿用曲、 1）. 心=畀（喫血加0）r/cl）.单轴情况演化: 首先假定单轴应力-应变关系可以通过下式转化成应力-塑性应变关系:.pi二对_ *上切打式中下表t c分别代表拉压。和，是拉压时的等效塑性应变率，和評=人兰出是拉压等型塑性应变，卩是温度,）是其它预定义常变量。在单 轴拉压情况下有效塑性应变率为：= e? . in uniaxial ten
8、sion and；. .帖必丄）=需.in uniaxial comprrssion.这一节里面我们约定一是正数，它代表的是单压时的应力值，即J-。正如在图4.5.2-1 中显示的那样，当从应力-应变曲线的应变软化段卸载时，可以发现卸载的响应是退化了的， 也就是说材料的弹性模量看起来变小了（损伤了）。弹性刚度的损伤在拉压试验中表现是大 不相同的。但在拉压两种情况中，随着塑性变形的增加损伤效果都是越来越明显的。混凝土 的损伤响应由两个独立的单轴损伤变量：和汀，控制，他们是塑性应变、温度和其它行变量的 函数。(4.5.2-5) = W&站 （o is心二心（型J汰（）!）-图4.521混凝土单轴拉
9、和压应力-应变曲线单轴刚度退化变量是等效塑性应变的非减函数，他们的取值范围在0 （无损伤）到1 （完全损 伤）之间。如果表示材料的初始弹性刚度，那么在单轴拉压下的应力-应变关系分别为 囚二门如E。怎帑）,o古=（1 必）场匕-聖），在单轴加载条件下，裂纹是沿着与应力垂直方向发展的。裂纹的成核和扩展就造成了界面有 效承载面积的减小，因此就导致了有效应力的增加。在单轴压是这种承载面积减小的效果还 要稍好一点，因为开始是裂纹基本上是平行于应力方向扩展的，但是当压碎发展到比较厉害 时有效承载面积也将显著地减小。那么有效单轴内聚力和形式如下有效单轴内聚力决定了屈服（破坏）面的大小。单轴循环加载在单轴循环
10、加载条件下，刚度退化机制比较复杂，它设计到预先存在裂纹的开闭问题和裂纹 间的相互作用问题。试验观察发现，但循环加载的应力符号变号是反向加载的刚度有所恢复。 这种刚度恢复也称之为“单边效应”它是混凝土循环加载的一个显著特点。特别是当应力有 拉变为压是，效应很明显，这时压应力是的受拉形成的裂纹闭合从而是受压刚度得到恢复。 混凝土塑性损伤模型假定弹性模量按标量减小变量退化-!是材料的初始（无损）模量。这个关系式在拉压曲线中都是成立的，刚度减小变量d是应力状态和单轴损伤变量和川的 函数，在单轴循环条件下ABAQUS假定下式成立：.i：l 出）=1 旳也）. 0 矢* 1,（4,Ct.2 G）式中-和、
11、应力状态的函数，引入他们是为了反应由于反向加载是刚度恢复效应，他们定义为: 旳=1 一 禎吓（61）：0 1.0 wr 00 if rtLJ o仏=I - Wa(l - r*(an): 其中，*权系数-和：这里假定为材料参数，他们分别控制应力反向是的刚度恢复能力。举例来说, 考虑图4.5.2-2荷载有拉变成压的情况。假定材料没有初始预损伤，也就是= 及-, 那么此时有（1- = （! -f） = （1- （ -拉应力（门）时，正如预计的那样-。反之压应力（.）时, = 二.。如果一那么 ，材料恢复到受压无损状态，反之，若* 时，一 =，材料没有刚度恢复。当：在0-1之间取值时表示刚度只能部分恢
12、复。图4.5.2-2受压刚度恢复参数：效应的示意图TtoEod-dt)E0单轴循环加载时的等效塑性演化方程也可以进行推广如下:V -宀血c j r c 1 I丄謬二_(1 隹监(4.5,2-7)它在单拉或单压就退化为方程4.5.2-4的形式。多轴情况有必要把硬化变量的演化规律推广到多轴情况下，在Lee and Fenves (1998)的工作基础上, 假定有效塑性应变率可由下式计算得到：(4,5.2-8)式中：“；-和；，:分别是塑性应变率张量，的最大和最小主值。是拉压应力权重系数，若有效应力张量三个主值全是正时为1，反之为0。Macauley运算 定义为：；门一 ;1 - o单轴加载情况下方
13、程4.5.2-8退化为单轴定义式4.5.2-4和汶-i一4.5.2-7，因为此时单拉时.，单压时。若果对塑性应变率张量的主值进piqj哼吧pt行排序如：-，那么多轴普通应力条件下等效塑性盈利率演化可以写成一下矩阵形式：-(r)0)3 q Lc200-( 一 F0)，弹性刚度退化混凝土塑性损伤模型认为混凝土的弹性刚度退化时各向同性的，且可以用一个单标量写成如 下形式：- (1 -旳 Df: 0d4.5.2 10)形式上与单轴相同，只是现在通过应力权重系数将它推广到多轴情况了：趴=1 () L8e = 1 We(l 一0 ir- I.显然，很容易验证方程4.5.2-10的标量退化式与单轴加载时是一
14、致的。很多准脆性材料(混凝土)的试验表明，当拉应力换到压应力时由于裂纹闭合受压刚度将会 恢复。但是另一方面，当受压是的微裂纹压碎时，由受压换到受拉时的受拉刚度将不会恢复。 鉴于此，ABAQUS默认条件下，假定-及一-即只有受压刚度恢复而没有受拉刚度 恢复。图4.5.2-3就是默认条件下的一个应力循环的曲线图 图4.5.2-3默认条件下3 -匚，一】)单轴应力循环曲线图(拉-压-拉)屈服条件本模型的屈服条件基于Lubliner等人（1989）建议的屈服函数，它综合了 Lee and Fenves（1998）的修正以考虑拉压不同时强度的不同演化规律。用有效应力表达时的屈服函数为：F回評） （可一（厅 了（一厅映/）厅&（岂）W（）（15,2- U）式中g和一是无量纲材料参数-1- 丁P =矿：I时有效静水压力，是Mises等效应力，是有效应力张量的偏量部分，而一是戸的代数最大主值，函数 形式如下式中和分别为有效拉压内聚力。在双轴受压时,11方程4.5.2-11就退化为
