摘要:根据碾压混凝土拱坝的结构特性与碾压混凝土材料的破坏特征,以sidoroff各向异性损伤理论为基础,建立了碾压混凝土正交异性损伤本构模型;基于正交异性损伤模型的损伤传递方式,导出了碾压混凝土材料的损伤演化方程。以碾压混凝土在复杂应力下的试验资料,对碾压混凝土正交异性损伤本构模型进行了验证。
关键词:碾压混凝土 混凝土 正交异性损伤 本构模型
高碾压混凝土拱坝结构设计面临的一个主要问题是如何有效地消散坝体在建造期间与运行期间的温度应力[1]。目前,工程师倾向于在坝体中设置少量的伸缩缝(包括横缝与诱导缝),通过对带缝碾压混凝土拱坝进行数值仿真计算和物理模型试验,分析坝体的受力规律与破坏机制,以期有组织地解决温度应力问题[2~4]。由于碾压混凝土材料的配合比选择、组分比例和成型方法与常规混凝土有很大的不同,其温度特性、变形特性和强度特性也与常规混凝土有较大的不同[5]。因此对于碾压混凝土材料,有必要研究其受荷后的破坏机理,建立合理的本构模型,为碾压混凝土拱坝进行数值计算服务。
1 碾压混凝土材料损伤本构模型
碾压混凝土损伤开裂具有明显的正交异性特征。碾压混凝土的单向拉伸与压缩试验破坏特征表明,拉伸时混凝土开裂方向与拉力垂直,其损伤与拉力方向相同;压缩时混凝土开裂方向与压力平行,其损伤与压力方向垂直。当混凝土处于复杂应力状态时,不同的应力比对应力应变之间的函数关系发生显著影响。
为了分析弹性材料的各向异性损伤,sidoroff[6]等人提出了能量等效假设,认为受损材料的弹性余能和无损材料的弹性余能在形式上相同。只需将其中的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´]=[r]t[e*(d)][r] |
(3) |
[r]为局部坐标与整体坐标之间的转换矩阵[7],
|
|
其中:li,mi,ni(i=1,2,3)分别为各主应力方向在总体坐标系中的方向余弦。
2 碾压混凝土的正交异性损伤演化方程
大体积混凝土材料在初始状态认为是各向同性的线弹性体,随着承载后材料的损伤,表现出明显的正交异性特征。大量的碾压混凝土的单向拉伸与压缩试验破坏特征表明,拉伸时混凝土开裂方向与拉力垂直,其损伤是与拉力方向相同,可称为“直接损伤”;压缩时混凝土开裂方向与压力平行,其损伤是与压力方向垂直,可称为“传递损伤”。当混凝土处于复杂应力状态时,不同的应力比对应力-应变之间的函数关系发生显著影响。在大连理工大学土木系结构室所做的双轴应力状态下碾压混凝土的力学特性研究试验中,双轴拉压状态下,一个方向压应力的增大可明显的降低另一方向的抗拉强度;而在双轴拉状态下,一个方向的拉应力对另一方向的抗拉强度影响很小[8]。基于以上的试验研究结果,可假定:(1)碾压混凝土材料在拉应变方向发生损伤同时,与之正交方向的压应变对拉应变方向的损伤有影响;(2)压应变方向损伤不受其它正交方向应变状态的影响。(3)受拉方向对其它正交方向的损伤无影响。
2.1 损伤“传递”模型 首先将下文中所用的符号变量集中说明如下:ω为损伤变量值;a为无损状态时的横截面积;ω′为裂缝扩展后损伤变量值;
为受损后有效承载面积;a′0为初始损伤面积;a′0i为初始损伤时材料内每一小裂纹或缺陷等效圆形区域面积;d为圆形区域a′01的直径;w为裂缝扩展后的宽度。a′为裂纹扩展后损伤区域面积;k、k′ 为材料参数;a′i为裂缝扩展后材料内每一小裂纹或缺陷等效圆形区域面积;
|
材料在某一方向的劣化可认为是由于材料的微缺陷与内部裂纹导致的有效承载面积减小,损伤变量ω=(a- |
假定材料在z方向受到压应变,考虑压应变对损伤区域a′0i的作用。按照线弹性断裂力学分析,在一个含有椭圆裂纹的无限介质中,当受到单向压缩时,在平行于加载方向的裂纹尖端将产生切向拉应力,使裂纹扩展[3]。这里,按照应变作用方式描述裂纹的扩展,如图2所示,在z方向压应变作用下,损伤面张开。当裂纹张开达到一定的宽度,裂纹便发生新的扩展,使x方向有效承载面积减小,造成新的损伤。
|
图1 x方向的材料微缺陷与内部裂纹模型 |
图2 压应变作用下的裂纹张开示意 |
|
设某一微损伤区域a′0i的直径为d,在损伤面张开时,裂缝张开宽度达到一定的值裂缝wc将发生扩展:这里wc为临界状态裂缝宽度。对于碾压混凝土材料,可参照常规混凝土裂缝宽度w与裂缝长度d为线性关系[10]:w=kd,k为材料参数。 下面由z方向压应变确定裂缝宽度w.直径为d的微损伤区域a′0i在z应变εz作用正直径变为d(1-εz).由图3几何关系,可算得裂缝宽度w. |
图3 裂缝宽度与裂缝长度关系 |
|
|
(4) |
|
当w>wc时,裂缝长度扩展为d′,而wc/w=d/d,所以,再由式(4)得 |
|
|
|
(5) |
微损伤区域面积由a′0i扩展为a′i
|
|
则在x方向上材料的总损伤面积a′为
|
|
(6) |
新的损伤变量值为ω′
|
|
(7) |
由于损伤导致有效承载面积的减小,有效应力随之升高。定义有效应力张量为
|
|
(8) |
|
lemaitre[11]于1971年提出了等效应变假设,认为受损材料的变形可以只通过有效应力来体现。即损伤材料的本构关系可以采用无损形式,只要将其中的应力σij替换为有效应力 |
|
d=1- |
(9) |
因此d可用损伤变量ω来定义,
|
d=ω |
(10) |
z方向压应变对x方向损伤的作用表示为,
|
d=apd0=-εz(2+εz)/ k′/d0 (εz<0) |
(11) |
|
ap=-εz(2+εz) k′(εz<0) |
(12) |
式中:d0为未考虑传递损伤时x方向的损伤度。
同样,另一方向(y方向)也为压应变作用时,也有
|
a′p=-εy(2+εy) k′(εy<0) |
(13) |
此时x方向的损伤为
|
d=apa′pd0 (εz<0,εy<0) |
(14) |
显然,公式(12、13)还必须有,当ap<1,a′p<1时,令ap=1,a′p=1.2.2 损伤演化模型 在损伤演化模型中,选择mazars一维损伤演化模型[12]描述受拉方向x的“初始损伤”d0.在单轴受拉状态时mazars一维损伤演化方程,
|
{ |
d0=0 |
(0≤ε≤εp) |
(15a) |
|
d0=1-/εp(1-at)/ε-at/exp[bl(ε-εp)](15a) |
(ε>εp) |
在单轴受压状态时mazars一维损伤演化方程,
|
{ |
d0=0 |
(0≥ε≥εp) |
(15b) |
|
d0=1-[εp(1-ac)/ε]-[acexp[bc(ε-εp)] |
(ε<εp) |
所以,考虑损伤传递效应的损伤演化方程定义为,在受拉方向,
|
{ |
d=0 |
(0≤ε≤εp) |
|
(16a) |
|
| ||||
|
|
在受压方向,
|
{ |
|
(16b) |
式中:at bt ac bc为损伤阈值;当ap<1,a′p<1时,令ap=1,a′p=1,公式(16a、16b)退化为公式(15a、15b).




