ABAQUS混凝土损伤塑性模型损伤因子

1、混凝土损伤因子的定义BY lizhenxian271损伤因子的定义损伤理论最早是1958年Kachanov提出来用于研究金属徐变的。所谓损伤，是指在各种 加载条件下，材料内凝聚力的进展性减弱，并导致体积单元破坏的现象，是受载材料由于微 缺陷(微裂纹和微孔洞)的产生和发展而引起的逐步劣化。损伤一般被作为一种“劣化因素” 而结合到弹性、塑性和粘塑性介质中去。由于损伤的发展和材料结构的某种不可逆变化，因而不同的学者采用了不同的损伤定义。一般来说，按使用的基准可将损伤分为：(1)微观基准量1,空隙的数目、长度、面积、体积；2空隙的形状、排列、由取向所决定的有效面积。(2)宏观基准量1、弹性常数、屈服应

2、力、拉伸强度、延伸率。2、密 度、电阻、超声波波速、声发射。对于第一类基准量，不能直接与宏观力学量建立物性关系，所以用它来定义损伤变量的 时候，需要对它做出一定的宏观尺度下的统计处理 (如平均、求和等)。对于第二类基准量，一般总是采用那些对损伤过程比较敏感， 在实验室里易于测量的量， 作为损伤变量的依据。由于微裂纹和微孔洞的存在，微缺陷所导致的微应力集中以及缺陷的相互作用，有效承 载面积由A减小为A'。如假定这些微裂纹和微孔洞在空间各个方向均匀分布，A与法向无关,这时可定义各向同性损伤变量 D为D= ( A- A' / A事实上，微缺陷的取向、分布及演化与受载方向密切相关，因此

3、材料损伤实际上是各向 异性的。为描述损伤的各向异性，可采用张量形式来定义。损伤表征了材损伤是一个非负的 因子，同时由于这一力学性能的不可逆性，必然有也0 dt2有效应力定义Cauchy有效应力张量A/A /(1 D)一般情况下，存在于物体内的损伤(微裂纹、空洞)是有方向性的。当损伤变量与受力面 法向相关时，是为各向异性损伤；当损伤变量与法向无关时，为各向异性损伤。这时的损伤变 量是一标量。3等效性假设损伤演化方程推导一般使用两种等效性假设，一种是应变等效性假设，另一种是能量等效性假设。采用能量等效性假设可以避免采用应变等效假设而使得各向异性损伤模型中的有 效弹性矩阵不对称的问题.以下对两种假设

4、进行简要的介绍。(1)应变等效性假设1971年Lematire提出，损伤单元在应力 作用下的应变响应与无损单元在定义的有效应 力作用下的应变响应相同。在外力作用下受损材料的本构关系可采用无损时的形式，只要 把其中的Cauchy应力简单地换成有效应力即可。在一维线弹性问题中，如以表示损伤弹性应变则E。E0(1 D) E由此可得E0(1 D)(2)能量等效性假设Sidiroff的能量等价原理，应力作用在受损材料产生的弹性余能与作用在无损材料产生的 弹性余能在形式上相同，只要将应力改为等效应力，或将弹性模量改为损伤时的等效弹性模量 即可。无损伤材料弹性余能:Woe-o2E0等效有损伤材料弹性余能:(

5、2于是得EdEo（1 D）2,则进一步可以得到Eo(1D)24单轴情况下损伤演化方程的介绍因为abaqus中用到的损伤塑性模型，在帮助文件中并没有给给出如何定义损伤。如果用 户没有自定义损伤因子，充其量是带强度硬化的塑性模型。且在abaqus中用户需要输入的只是单轴下的，相应的损伤因子与开裂应变（或开裂位移）文中单指拉伸强化。混凝土受拉时，主要表现为脆性，具有较小的不可逆变形，因此工程上常把它视为弹性 材料。从1980年开始，各国学者用损伤理论分析混凝土受载后的力学状态，提出了各种损伤 模型，并首先应用于研究材料受拉的情况。建立损伤模型可以用能量的方法，也可以用几何 的方法，而最简单又实用的是

6、用半实验半理论的方法。下面介绍一些关于混凝土材料主要的 损伤模型。（1）经典损伤理论的混凝土损伤变量计算方法：D=1-Es/Eo其中Es为应力应变曲线上任一点的割线模量。（经过分析，该公式为弹性损伤模型，计算的损伤变量偏大，不适合ABAQUS塑性损伤计算，输入后会报错）(2) Loland 模型该模型认为，在应力值以前（f ）,裂纹仅在体元中萌生和扩展，且保持在一个很小的限度内；在应力峰值以后（fu ）裂纹主要在破坏区内不稳定扩展、开裂。E1D。(1 D) D D° C1式中E为初始弹性模量。2. fu时，混凝土损伤材料损伤演化方程为E f1 D°(1D)DfC2( f)由

7、边界条件:0、D|1可以定出常数、C1、C2：CiD°f、C2Df在应力达到峰值应力前，Loland模型假定有效应力和应变之间为线性关系，与试验结果 比较吻合。而在峰值应变之后，有效应力为常数，并由此得到损伤变量与应变为线性关系, 与实际情况不符，是一种近似的模拟。(3) Mazars模型混凝土损伤变量计算方法:D 1 f(1 A) Aexp B( f)式中:A,B为材料常数，由实验确定，对一般混凝土材料0.7<A<1, 104<B<105, fe为对应于 f的应变。(计算公式较复杂，但有实验数据说话，可信度高，可以参考)对于一般混凝土单轴压缩时公式类推其中1

8、<A<1.5, 103<B<2X 103。abaqus5基于规范推荐的应力一应变曲线推导的单轴损伤演化方程损伤演化方程的推导由于选取不一样的假设前提而得到不同的方程。往往在应用 算混凝土损伤时，需要输入参数多，在其帮助文件中又没有给出具体损伤的定义，有时输入 数据容易出错 (During the conversion from crushing to plastic strain abaqus found negative and/or decreasing values of plastic strain. Verify that the degradation da

9、ta under *concrete compression这时可以根据其对应的damage is correct。在工程应用中并没有现成的混凝土应力应变曲线,相应的抗拉强度ft、抗压强度fc和弹性模量E。，近似的简化。在达到极限应力时假设其应力 应变曲线为直线，此阶段没有损伤，在极限应力峰值后采用规范给出的应力应变曲线，采用 能量等效原理得出abaqus输入的数据。f, y / f)1.0压缩应力应变曲线y / fc上式是通过混凝土结构设计规范(GB 50010-2002)»规范简化而来的，能满足大体积混凝土结构的基本需求。混凝土单轴受拉的应力-应变曲线方程可按下列公式确定(在计算

10、中前半部分认为线弹性，损伤只发生在峰值后)xt(x 1)1.7 x同理混凝土单轴受拉的应力-应变曲线方程可按下列公式确定:xy ；72d(x 1) x其中 t 0.312 ft1 ；d(x 1)2 x ,d 0.157 fc0.785 0.905采用上文中提到的能量等效原理可得出:单轴受拉损伤方程:_1.7_t(x 1) x单轴受压损伤方程:其它各种损伤可参考相关损伤力学书籍例子："MATERIAL, NAME=a0_baduan*ELASTIC2.31E10,0.2"density2400"expansion0.826e-5"Concrete Dama

11、ged Plasticity35., 0.1, 1.16,0.67,0."CONCRETE COMPRESSION HARDENING24.3E+6,0.000020.78E+6, 0.678E-0318.32E+6, 0.105E-0216.11E+6, 0.141E-0210.32E+6, 0.271E-026.45E+6, 0.445E-022.63E+6, 0.104E-011.26E+6, 0.210E-010.82E+6, 0.315E-01"CONCRETE COMPRESSION DAMAGE0.0000, 0.00000.2450,0.678E-30.34

12、37,0.105E-20.4242,0.141E-20.6238,0.271E-20.7571,0.445E-20.8959,0.104E-10.9492,0.210E-10.9664,0.315E-1"CONCRETE TENSION STIFFENING 1.780E+6, 0.0000 1.351E+6,0.725E-40.996E+6,0.150E-30.536E+6,0.401E-30.389E+6,0.638E-30.313E+6,0.873E-30.266E+6,0.111E-20.233E+6,0.134E-20.121E+6,0.346E-2"CONCRETE TENSION DAMAGE 0.0000, 0.0000 0.3317,0.725E