ABAQUS_混凝土损伤塑性模型_损伤因子

混凝土损伤因子的定义

BY lizhenxian27

1 损伤因子的定义

损伤理论最早是1958年Kachanov提出来用于研究金属徐变的。所谓损伤，是指在各种加载条件下，材料内凝聚力的进展性减弱，并导致体积单元破坏的现象，是受载材料由于微缺陷(微裂纹和微孔洞)的产生和发展而引起的逐步劣化。损伤一般被作为一种“劣化因素”而结合到弹性、塑性和粘塑性介质中去。

由于损伤的发展和材料结构的某种不可逆变化，因而不同的学者采用了不同的损伤定义。一般来说，按使用的基准可将损伤分为:

(1) 微观基准量

1,空隙的数目、长度、面积、体积;

2空隙的形状、排列、由取向所决定的有效面积。

(2) 宏观基准量

1、弹性常数、屈服应力、拉伸强度、延伸率。

2、密度、电阻、超声波波速、声发射。

对于第一类基准量，不能直接与宏观力学量建立物性关系，所以用它来定义损伤变量的时候，需要对它做出一定的宏观尺度下的统计处理(如平均、求和等)。

对于第二类基准量，一般总是采用那些对损伤过程比较敏感，在实验室里易于测量的量，作为损伤变量的依据。

由于微裂纹和微孔洞的存在，微缺陷所导致的微应力集中以及缺陷的相互作用，有效承

载面积由

A 减小为A ’。如假定这些微裂纹和微孔洞在空间各个方向均匀分布，A ’与法向无关，这时可定义各向同性损伤变量D 为

D= ( A- A ’ )/ A

事实上，微缺陷的取向、分布及演化与受载方向密切相关，因此材料损伤实际上是各向异性的。为描述损伤的各向异性，可采用张量形式来定义。损伤表征了材损伤是一个非负的因子，同时由于这一力学性能的不可逆性，必然有

0dD

dt

≥ 2有效应力

定义Cauchy 有效应力张量'σ

''//(1)A A D σσσ==-

一般情况下，存在于物体内的损伤(微裂纹、空洞)是有方向性的。当损伤变量与受力面法向相关时，是为各向异性损伤;当损伤变量与法向无关时，为各向异性损伤。这时的损伤变量是一标量。

3等效性假设

损伤演化方程推导一般使用两种等效性假设，一种是应变等效性假设，另一种是能量等效性假设。采用能量等效性假设可以避免采用应变等效假设而使得各向异性损伤模型中的有效弹性矩阵不对称的问题.以下对两种假设进行简要的介绍。

(1) 应变等效性假设

1971年 Lematire 提出，损伤单元在应力σ作用下的应变响应与无损单元在定义的有效应力'σ作用下的应变响应相同。在外力作用下受损材料的本构关系可采用无损时的形式，只要

把其中的

Cauchy 应力简单地换成有效应力即可。在一维线弹性问题中，如以ε表示损伤弹性应变则

'

'

0(1)

E E D E σσ

σ

ε==

=

-

由此可得0(1)E D σε=- (2)能量等效性假设

Sidiroff 的能量等价原理，应力作用在受损材料产生的弹性余能与作用在无损材料产生的弹性余能在形式上相同,只要将应力改为等效应力，或将弹性模量改为损伤时的等效弹性模量即可。

无损伤材料弹性余能：

2

2e W E σ=

等效有损伤材料弹性余能：

'2

2e d

d

W

E σ=

于是得20(1)d E E D =-，则进一步可以得到

20(1)E D σε=-

4单轴情况下损伤演化方程的介绍

因为abaqus 中用到的损伤塑性模型，在帮助文件中并没有给给出如何定义损伤。如果用户没有自定义损伤因子，充其量是带强度硬化的塑性模型。且在abaqus 中用户需要输入的只是单轴下的，相应的损伤因子与开裂应变（或开裂位移）文中单指拉伸强化。

混凝土受拉时，主要表现为脆性，具有较小的不可逆变形，因此工程上常把它视为弹性材料。从1980年开始，各国学者用损伤理论分析混凝土受载后的力学状态，提出了各种损伤

模型，并首先应用于研究材料受拉的情况。建立损伤模型可以用能量的方法，也可以用几何的方法，而最简单又实用的是用半实验半理论的方法。下面介绍一些关于混凝土材料主要的损伤模型。

（1）经典损伤理论的混凝土损伤变量计算方法：

D=1-E s /E 0

其中E s 为应力应变曲线上任一点的割线模量。

（经过分析，该公式为弹性损伤模型，计算的损伤变量偏大，不适合ABAQUS 塑性损伤计算，输入后会报错）

（2）Loland 模型

该模型认为，在应力值以前（f

ε

ε≤），裂纹仅在体元中萌生和扩展，且保持在一个很

小的限度内；在应力峰值以后（f u εεε≤≤）裂纹主要在破坏区内不稳定扩展、开裂。

1．f εε≤时，混凝土损伤材料损伤演化方程为

'

1E D εσ=- '(1)D σσ=- 01D D C β

ε=+

式中E 为初始弹性模量。 2．f

u εεε≤≤时，混凝土损伤材料损伤演化方程为

'

1f E D εσ=

-

'

(1)D σσ=- 2()f f D D C εε=+- 由边界条件：|f

f

εεσσ==、|0f d d εεσ

ε

==、|1u D εε==可以定出常数β、1C 、2C ： f

f f E σβεσ=- 、0

111f D C εβ

-=+ 、21f

f u D C εε-=

在应力达到峰值应力前，Loland 模型假定有效应力和应变之间为线性关系，与试验结果比较吻合。而在峰值应变之后，有效应力为常数，并由此得到损伤变量与应变为线性关系，与实际情况不符，是一种近似的模拟。