水利工程论文-碾压混凝土的正交异性损伤本构模型研究.doc

水利工程论文-碾压混凝土的正交异性损伤本构模型研究摘要：根据碾压混凝土拱坝的结构特性与碾压混凝土材料的破坏特征，以Sidoroff各向异性损伤理论为基础，建立了碾压混凝土正交异性损伤本构模型；基于正交异性损伤模型的损伤传递方式，导出了碾压混凝土材料的损伤演化方程。以碾压混凝土在复杂应力下的试验资料，对碾压混凝土正交异性损伤本构模型进行了验证。关键词：碾压混凝土混凝土正交异性损伤本构模型高碾压混凝土拱坝结构设计面临的一个主要问题是如何有效地消散坝体在建造期间与运行期间的温度应力1。目前，工程师倾向于在坝体中设置少量的伸缩缝(包括横缝与诱导缝)，通过对带缝碾压混凝土拱坝进行数值仿真计算和物理模型试验，分析坝体的受力规律与破坏机制，以期有组织地解决温度应力问题24。由于碾压混凝土材料的配合比选择、组分比例和成型方法与常规混凝土有很大的不同，其温度特性、变形特性和强度特性也与常规混凝土有较大的不同5。因此对于碾压混凝土材料，有必要研究其受荷后的破坏机理，建立合理的本构模型，为碾压混凝土拱坝进行数值计算服务。1碾压混凝土材料损伤本构模型碾压混凝土损伤开裂具有明显的正交异性特征。碾压混凝土的单向拉伸与压缩试验破坏特征表明，拉伸时混凝土开裂方向与拉力垂直，其损伤与拉力方向相同；压缩时混凝土开裂方向与压力平行，其损伤与压力方向垂直。当混凝土处于复杂应力状态时，不同的应力比对应力应变之间的函数关系发生显著影响。为了分析弹性材料的各向异性损伤，Sidoroff6等人提出了能量等效假设，认为受损材料的弹性余能和无损材料的弹性余能在形式上相同。只需将其中的Cauchy应力换为等效应力。在此基础上，提出了一种各向异性损伤模型(1)(2)基于式(1)，建立碾压混凝土材料的应力-应变关系为(在主轴系内)：E*(D)=(+2D)(1-D1）2（1-D1）（1-D2）（1-D1）（1-D3）对（+2G）（1-D2）2（1-D2）（1-D3）称（+2G）（1-D3）2这里为拉梅系数，=E(t)(1+)(1-2)；G为剪切模量，G=E(t)/2(1+)；E(t)料的弹性模量，是时间的函数；为泊松比Di为i主应变方向的损伤度。将主轴坐标系下的本构关系转换到总体坐标系下，得到总体坐标系下的应力-应变关系矩阵D=RTE*(D)R(3)R为局部坐标与整体坐标之间的转换矩阵7，li,mi,ni(I=1,2,3)分别为各主应力方向在总体坐标系中的方向余弦。2碾压混凝土的正交异性损伤演化方程大体积混凝土材料在初始状态认为是各向同性的线弹性体，随着承载后材料的损伤，表现出明显的正交异性特征。大量的碾压混凝土的单向拉伸与压缩试验破坏特征表明，拉伸时混凝土开裂方向与拉力垂直，其损伤是与拉力方向相同，可称为“直接损伤”；压缩时混凝土开裂方向与压力平行，其损伤是与压力方向垂直，可称为“传递损伤”。当混凝土处于复杂应力状态时，不同的应力比对应力-应变之间的函数关系发生显著影响。在大连理工大学土木系结构室所做的双轴应力状态下碾压混凝土的力学特性研究试验中，双轴拉压状态下，一个方向压应力的增大可明显的降低另一方向的抗拉强度；而在双轴拉状态下，一个方向的拉应力对另一方向的抗拉强度影响很小8。基于以上的试验研究结果，可假定：(1)碾压混凝土材料在拉应变方向发生损伤同时，与之正交方向的压应变对拉应变方向的损伤有影响；(2)压应变方向损伤不受其它正交方向应变状态的影响。(3)受拉方向对其它正交方向的损伤无影响。2.1损伤“传递”模型首先将下文中所用的符号变量集中说明如下：为损伤变量值；A为无损状态时的横截面积；为裂缝扩展后损伤变量值；为受损后有效承载面积；A0为初始损伤面积；A0i为初始损伤时材料内每一小裂纹或缺陷等效圆形区域面积；d为圆形区域A01的直径；w为裂缝扩展后的宽度。A为裂纹扩展后损伤区域面积；K、K为材料参数；Ai为裂缝扩展后材料内每一小裂纹或缺陷等效圆形区域面积；材料在某一方向的劣化可认为是由于材料的微缺陷与内部裂纹导致的有效承载面积减小，损伤变量=(A-)/A；A为无损状态时的横截面积；为受损后有效承载面积；A0为初始损伤面积，A0A-.这里，材料的微缺陷与内部裂纹描述为图1所示模型：材料中存在许多的微小裂纹与微小缺陷，这些微小裂纹与微小缺陷对某一方向X的损伤可等效为许多的微小圆形裂开区域A0i，使得这一向上的有效承载面积减小。材料在X方向的初始损伤面积A0A0i假定材料在Z方向受到压应变，考虑压应变对损伤区域A0i的作用。按照线弹性断裂力学分析，在一个含有椭圆裂纹的无限介质中，当受到单向压缩时，在平行于加载方向的裂纹尖端将产生切向拉应力，使裂纹扩展3。这里，按照应变作用方式描述裂纹的扩展，如图2所示，在Z方向压应变作用下，损伤面张开。当裂纹张开达到一定的宽度，裂纹便发生新的扩展，使X方向有效承载面积减小，造成新的损伤。图1X图2压应变作用下的裂纹张开示意设某一微损伤区域A0i的直径为d，在损伤面张开时，裂缝张开宽度达到一定的值裂缝wc将发生扩展：这里wc为临界状态裂缝宽度。对于碾压混凝土材料，可参照常规混凝土裂缝宽度w与裂缝长度d为线性关系10：wKd，K为材料参数。下面由Z方向压应变确定裂缝宽度w.直径为d的微损伤区域A0i在Z应变z作用正直径变为d(1-z).由图3几何关系，可算得裂缝宽度w.图3裂缝宽度与裂缝长度关系(4)当wwc时，裂缝长度扩展为d，而wc/wd/d，所以，再由式(4)得(5)微损伤区域面积由A0i扩展为AiX方向上材料的总损伤面积A为(6)新的损伤变量值为(7)由于损伤导致有效承载面积的减小，有效应力随之升高。定义有效应力张量为(8)Lemaitre11于1971年提出了等效应变假设，认为受损材料的变形可以只通过有效应力来体现。即损伤材料的本构关系可以采用无损形式，只要将其中的应力ij替换为有效应力。损伤体现为把无损时的弹性模量E减小为损伤后的弹性模量(1-)E.以割线模量的变化定义损伤D，表示为D1/E(9)因此D可用损伤变量来定义，D(10)z方向压应变对x方向损伤的作用表示为，DApD0-z(2+z)/K/D0(z0)(11)Ap-z(2+z)K(z0)(12)式中：D0为未考虑传递损伤时x方向的损伤度。同样，另一方向(y方向)也为压应变作用时，也有Ap-y(2+y)K(y0)(13)此时x方向的损伤为D=ApApD0(z0,y0)(14)显然，公式(12、13)还必须有，当Ap1,App)在单轴受压状态时Mazars一维损伤演化方程，D00(0p)(15b)D01-p(1-Ac)/-AcexpBc(-p)(p)D=0(0p)(16a)在受压方向，(16b)式中：AtBtAcBc为损伤阈值；当Ap1，Ap1时，令Ap=1，Ap=1，公式(16a、16b)退化为公式(15a、15b).2.3参数标定损伤演变方程(16a、16b)的未需要试验所得到的应力、应变资料来确定。材料常数At、Bt、Ac、Bc采用选择初值试算方法通过碾压混凝土单轴拉和单轴压状态应力-应变关系曲线确定。选定碾压混凝土材料常数At=1.105，Bt=5100，Ac=0.74，Bc=400，对应的应力应变关系曲线如图4、图5.图4单轴拉状态碾压混凝土应力应变关系曲线图5单轴压态碾压混凝土应力应变关系曲线DZ)按公式(16a)确定各试件的参数K，采用反复试算的方法计算得各试件的K值如表2.由于混凝土材料试验结果有较大的随机性和模糊性，确定K不应简单地取各试件的均值。参数K的频度统计规律见图5，对于碾压混凝土可选择K=1.510-31.810-3。取K=1.810-3由公式(16a)算得的各试件在达到极限强度时的拉应变方向损伤度D(如表2).用Mazars损伤模型同正交异性损伤模型进行比较。在三维情况下Mazars损伤演化方程12为：表2试件的材料参数K与损伤度D、D*编号K(10-3)DD*编号K(10-3)DD*ACC411.845.592.578ACT31.419.617.666ACC442.259.838.667ACT41.977.822.748ACC452.6031.0.730ACT51.426.500.624ACC462.481.983.713ACT61.504.539.646ACC212.8131.0.850ACT321.515.548.