连续损伤理论和其在混凝土上的应用.doc

连续损伤理论在混凝土上的应用作者：J.Mazars ，G Pijaudier —Cabot，学生们，美国土木工程师学会摘要：本文提出了一种对于不同混凝土模型的观点，该模型是基于连续损伤理论，在实验室de Mtecanique et Technologie（卡尚，法国）中推导出公式的每个公式都是建立在物理观测的基础上，由不可逆过程热力学框架而来受次生各向异性，延伸性，和方向性的影响，如封闭裂缝的讨论，并提出足够的损伤模型最后，数值代入的实现提供一个良好的破坏过程的描述，以及对混凝土和钢筋混凝土结构的行为准确的预测简介：许多材料，如混凝土，岩土，木材和复合材料的破坏是由于微裂纹的开展和闭合在结构分析中这种被称为损伤的现象最常是被视为应变软化为了建立这个破坏过程的模型，人们成功的提出了各种不同类型的本构关系，包括本构模型内蕴时间塑性理论（Bazant 1986年），塑性断裂理论（Dougill 1983年，Dragon and Mroz 1979年）;总应变模型（Gerstle et al.1980; Kostovos 1980年），随屈服极限减小的塑性理论（Wastiels1980年），以及最近的微平面模型（Bazant 1980年; Pande and Sharma 1982年）。

Kachanov（1958年）在1958年提出的徐变相关问题，使得最近连续损伤力学得以应用于渐进破坏的描述[金属和复合材料的静态失效（Dufailly 1980年，Ladeveze 1986，Lemaitre and Chaboche 1978年）; 材料的疲劳和徐变（Leckie 1978）]在20世纪80年代初，证实了损伤力学可以准确的模拟混凝土应变软化反应（Krajcinovic 1983年，Ladeveze 1983年，Lemaitre and Mazars 1982年）和可以用不可逆过程热力学为框架编写相应的本构关系的公式（Lemaitre and Chaboche 1985年）考虑到材料是作为一组变量和热力学势所描述的系统，本构关系系统得随损伤运动学条件而派生然而我们仍然应做出适当的势能和损伤变量的选择（标量，张量，等）各种逐渐复杂的模型的提出以及混凝土及钢筋混凝土构件的数值实现的提出其中热力学方法的优势，由可由我们选择的容许势能组成;为了说明我们的描述，我们把自己限定在由Laboratoire de Mecanique et Technologie实验室中，由J.L. Clement,F. Collombet，C. LaBorderie，A. Zaborski共同组成的一个研究小组成果中。

损伤模式在开始我们的分析之前，我们有必要回顾下混凝土反应的几个主要的方面，这将引导我们在理论公式的推导中做出适当的选择此阶段混凝土可视为由三个成分制成的一种复合材料：水泥基质（微孔材料），骨料，连接基质和骨料的过渡区（Maso 1982年）在这个区域中，水合混凝土的结晶是具有高度方向性的（由于管壁效应）它也是复合材料中最多孔的部分，因此，也是其最薄弱的区域微观损伤机制已经可以由不同的技术观测到：X -射线（Slate and Oleski 1963年），显微镜（Dhir and Sangha，1974年），或声发射（Terrien 1980）这些调查也已被一个旨在更好地理解破坏过程（例如，Maso 1982年）的模型所完成它建立了：（1）损伤出现在阈值后，主要是在位于过渡区和水泥基质;（2）不同的损伤模式存在于应力状态和应力史之间的连系中这两种类型的损伤是可以区分的：（1）水泥基质的微孔结构的破损是一个类型的损伤，这种损伤由在材料上的静水压力造成，同时可能导致聚合;（2）微裂纹的扩展最常位于水泥基质当荷载扩大，I型的开裂占主导地位，但裂缝也可能根据加载史扩展成II或III型断裂先端的摩擦也可能会带来额外的延展性。

图1（a）总结了这些针对作用在材料上应力状态的观测数据只有裂缝具有高度方向性并且于静水压力无关才会导致各项异性的存在当由于裂纹的闭合发生所产生的荷载符号的相反时，固有存在的裂缝，初始刚度的恢复都是可以测量的这种定向的现象称为“单边效应”，它是在梁受到循环加载时观察得到的（1987年Mazars和LaBorderie）上图1（b）和（c）所示的是损伤的宏观影响图中单轴的反应存在于拉伸和压缩记录下的形状和最大应力是不同的，提出不同的损伤运动学，但最初的线性弹性行为仍然存在损伤的增长会产生一个减少的卸载重载刚度和增加的永久应变下面将介绍不同损伤参数，它们是描述被视为连续材料的反应的变化损伤增长的方程将任意推出为最大的契合实验的数据然而，每个公式只限于说明在宏观效应的损伤，即损伤变量的类型和热力学势能的形式图1（a）损伤方式及混凝土的性能；混凝土在（b）压缩；和（c）拉伸情况下的实验性性能这里所用的方法与已在文献中提出不同的塑性能的塑性理论相似理论公式在恒温下混凝土可以用弹性应变张量 ，损伤D来描述，有效塑性应变记为 可能被定义为：上式中 是塑性应变张量率并且表明张量乘积由两个因子共同制约：损伤D的数学定义在这方面不须太精密，假设总应变率 的弹性和塑性变形是分区的， 。

每个平衡状态是由一个热力学势 ，D， 的函数 数值区分（ 是材料的质量密度）一般能使 满足热力学第一准则的是二次形式比能（1978年Lemaitre和Chaboche）在Kachanov和Lemaitre的解释之后，我们认为只有材料的弹性性能受到损伤的影响因此， 可表示为： （1）其中 是损伤和弹性应变函数，和 是有效塑性应变的函数应力张量 ，损伤能量释放速率Y，和有效应力 都是用比能来定义的： （2） 永久变形和损伤是不可逆的过程，导致机械能转换成热量，在表面上生成据Clausius Duhem不等式，能量耗散率 必须保持是正的：我们在由于损伤 和塑性 引起的能量耗散率的表达上可以这么区分：一个满足Clausius Duhem不等式的充分条件可以是 和 由于在模型中塑性的引入与经典的发展非常相似（Ladeveze 1983）。

我们将集中我们的注意力在把损伤引入到弹性本构关系中我们将采用弹性势能 ：是一个四阶对称张量即割线刚度矩阵（公式2） 是损伤D的函数此时我们认为公式中 的选择可以简化为考虑损伤刚度矩阵把5式代入2和4式中得到：由于 ，损伤能量释放率Y是一个正定二次型，即损伤增加时刚度下降满足Clausius Duhem不等式的充分条件是满足这个条件的损伤增量将由表面加载方程 来决定，其中K0是损伤初始临界值这应力状态的函数的唯一性是由将f选为应变的函数而不是应力的函数来确保的（两个应变张量可能与相同的应力状态有关）考虑到加载条件，损伤演变可以定义为：F（e）是一个由实验确定的应变正函数我们可能会注意到，式11中的条件设置与塑性下的加载条件相似在下面的章节中，我们将详细介绍不同的函数f和F表达式然而，这种方法最重要的假设是 的函数 的选择， 是建立在第一段中相关的宏观观测的基础上显然函数 可以从微观力学研究中引入，这构成下一步，它目前还正在研究中它旨在更好的理解如混凝土的非均质材料的微观和宏观行为。

目前，在例子中使用的势能的数学表达式的灵感来自于这样的研究（Ladeveze 1983年）每个选择对应到不同的近似单轴刚度和压缩分布，在Ladeveze（1983）和Pijaudier-Cabot（1985）中有更广泛的详述标量损伤模型 (Mazars 1984)在这种模式下，材料应该是弹性的行为，并保持各向同性我们称 为材料初始刚度矩阵，D为损伤真应力概念（1978年Lemaitre和Chaboche）导出了弹性能的下列形式：应力 和损伤能量释放率是由式3计算得到的，耗散率可以从式4中获得：损伤标量D的范围从原始材料的0到代表均匀应变条件下的破坏（零应力）的1假设延伸是裂纹扩展的原因，我们使用的表面荷载的灵感来自圣维南最大主应变准则，图2（a）中已给出其公式为：与 等效的应变： 硬化软化参数K（D）考虑到保留以前对材料的加载史，充分利用等效应变 最初，K（D）是K0的临界值在图1（a）中定义的A损伤的方式，这种表面荷载的表达方式阻碍这种模式的正确性三种不同损伤模型的反应：（a）一个标量损伤变量D的模型；（b）两个标量损伤变量Dt和Dc的模型;及（c）永久变形的各向异性模型拉伸或压缩的反应是由两个类型损伤Dt和Dc（Mazars 1986）的耦合的下列规律来描述。

Dt和Dc分别对应单轴拉伸和单轴压缩测量的损伤从图、1（b）和1（c）所描绘的反应，其发展是由式12中两个函数Ft和Fc组成的F的综合形式所给出的总损伤D是Dt和Dc的加权总和在式11中 和 的权值是一个应变状态的函数我们称 和 为只分别出现在正和负的主应力的张量，应变张量 ， 定义为：和 的权重由以下表达式（Mazars 1984）定义为：如果 ，则Hi=1，否则，Hi= 0在演化规律中参数K0，At，Bt，Ac和Bc是由圆柱体压缩试验和梁弯曲试验独立确定由于使用两个不同的运动损伤取决于主应力的符号，系数 和 定义了每个类型的损伤一般荷载的作用从式14可以验证，在单轴拉伸中， = 1， =0，D = Dt，在压缩中反之单侧损伤模型(Ladeveze 1983; Mazars 1985)我们可能可以有效的区分损伤是由于拉伸还是压缩，而不是由式11定义的损伤运动学的平均设置当材料受到循环荷载，先前的公式不能获得在压力逆转（1987年Mazars和LaBorderie）期间观测到的刚度恢复由于损伤不能削弱（Clausius Duhem不等式）两个独立的标量，所以我们用损伤Dt和Dc。

根据应力的符号，明显的损伤将会是正压力Dt或负压力Dc如果荷载是复杂的，损伤可能是Dt和Dc的综合我们将应力张量分解成一个正的 和负的 两部分该材料是假定保持弹性，势能 是：其中 拉伸和压缩由应力的符号区分解出式2，在式16中找到相应的本构关系：其中I是个体张量在Clausius Duhem不等式中，两个损伤的能量释放率关系到每一个损伤标量的表现：一个满足式3的充分条件是两个损伤率Dt和Dc一直是正的我们建议使用每个变量的不同损伤加载面，用能量释放率表示为：Kt和Kc是与式.10中K（D）类似的硬化软化参数每个损伤标量的演化规律是以下能量释放率的函数，Dt= Ft（Yt）和Dc= Fc（Yc）图2（b）显示在应力空间的范围的相应形状，该应力空间中材料的反应最初是弹性的它是从每个表面获得的弹性域的交集与以前的标量模型图.2（a）相似，考虑了拉伸和压缩的不对称性相比标量模型图. 2（a）循环反应有很大的不同，尽管应。