碾压混凝土的正交异性损伤本构模型研究(1)程

摘要：根据碾压混凝土拱坝的结构特性与碾压混凝土材料的破坏特征，以sidoroff各向异性损伤理论为基础，建立了碾压混凝土正交异性损伤本构模型；基于正交异性损伤模型的损伤传递方式，导出了碾压混凝土材料的损伤演化方程。以碾压混凝土在复杂应力下的试验资料，对碾压混凝土正交异性损伤本构模型进行了验证。

关键词：碾压混凝土 混凝土 正交异性损伤 本构模型

　　高碾压混凝土拱坝结构设计面临的一个主要问题是如何有效地消散坝体在建造期间与运行期间的温度应力^[1]。目前，工程师倾向于在坝体中设置少量的伸缩缝(包括横缝与诱导缝)，通过对带缝碾压混凝土拱坝进行数值仿真计算和物理模型试验，分析坝体的受力规律与破坏机制，以期有组织地解决温度应力问题^[2～4]。由于碾压混凝土材料的配合比选择、组分比例和成型方法与常规混凝土有很大的不同，其温度特性、变形特性和强度特性也与常规混凝土有较大的不同^[5]。因此对于碾压混凝土材料，有必要研究其受荷后的破坏机理，建立合理的本构模型，为碾压混凝土拱坝进行数值计算服务。

1 碾压混凝土材料损伤本构模型

　　碾压混凝土损伤开裂具有明显的正交异性特征。碾压混凝土的单向拉伸与压缩试验破坏特征表明，拉伸时混凝土开裂方向与拉力垂直，其损伤与拉力方向相同；压缩时混凝土开裂方向与压力平行，其损伤与压力方向垂直。当混凝土处于复杂应力状态时，不同的应力比对应力应变之间的函数关系发生显著影响。

　　为了分析弹性材料的各向异性损伤，sidoroff^[6]等人提出了能量等效假设，认为受损材料的弹性余能和无损材料的弹性余能在形式上相同。只需将其中的cauchy应力σ换为等效应力。在此基础上，提出了一种各向异性损伤模型

基于式(1)，建立碾压混凝土材料的应力-应变关系为(在主轴系内)：

这里λ为拉梅系数，λ=e(t)／(1+μ)(1-2μ)；g为剪切模量，g=e(t)/2(1+μ)；e(t)料的弹性模量，是时间的函数；μ为泊松比d_i 为i主应变方向的损伤度。

　　将主轴坐标系下的本构关系转换到总体坐标系下，得到总体坐标系下的应力-应变关系矩阵

［r］为局部坐标与整体坐标之间的转换矩阵^［7］，

И其中：l_i,m_i,n_i(i=1,2,3)分别为各主应力方向在总体坐标系中的方向余弦。

2 碾压混凝土的正交异性损伤演化方程

　　大体积混凝土材料在初始状态认为是各向同性的线弹性体，随着承载后材料的损伤，表现出明显的正交异性特征。大量的碾压混凝土的单向拉伸与压缩试验破坏特征表明，拉伸时混凝土开裂方向与拉力垂直，其损伤是与拉力方向相同，可称为“直接损伤”；压缩时混凝土开裂方向与压力平行，其损伤是与压力方向垂直，可称为“传递损伤”。当混凝土处于复杂应力状态时，不同的应力比对应力-应变之间的函数关系发生显著影响。在大连理工大学土木系结构室所做的双轴应力状态下碾压混凝土的力学特性研究试验中，双轴拉压状态下，一个方向压应力的增大可明显的降低另一方向的抗拉强度；而在双轴拉状态下，一个方向的拉应力对另一方向的抗拉强度影响很小^［8］。基于以上的试验研究结果，可假定：(1)碾压混凝土材料在拉应变方向发生损伤同时，与之正交方向的压应变对拉应变方向的损伤有影响；(2)压应变方向损伤不受其它正交方向应变状态的影响。(3)受拉方向对其它正交方向的损伤无影响。

2.1 损伤“传递”模型　首先将下文中所用的符号变量集中说明如下：ω为损伤变量值；a为无损状态时的横截面积；ω′为裂缝扩展后损伤变量值；为受损后有效承载面积；a′₀为初始损伤面积；a′_0i为初始损伤时材料内每一小裂纹或缺陷等效圆形区域面积；d为圆形区域a′₀₁的直径；w为裂缝扩展后的宽度。a′为裂纹扩展后损伤区域面积；k、k′ 为材料参数；a′_i为裂缝扩展后材料内每一小裂纹或缺陷等效圆形区域面积；

　　假定材料在z方向受到压应变，考虑压应变对损伤区域a′_0i的作用。按照线弹性断裂力学分析，在一个含有椭圆裂纹的无限介质中，当受到单向压缩时，在平行于加载方向的裂纹尖端将产生切向拉应力，使裂纹扩展^[3]。这里，按照应变作用方式描述裂纹的扩展，如图2所示，在z方向压应变作用下，损伤面张开。当裂纹张开达到一定的宽度，裂纹便发生新的扩展，使x方向有效承载面积减小，造成新的损伤。

微损伤区域面积由a′_0i扩展为a′_i

г蛟趚方向上材料的总损伤面积a′为

新的损伤变量值为ω′

由于损伤导致有效承载面积的减小，有效应力随之升高。定义有效应力张量为

因此d可用损伤变量ω来定义，

z方向压应变对x方向损伤的作用表示为，

式中：d₀为未考虑传递损伤时x方向的损伤度。

　　同样，另一方向(y方向)也为压应变作用时，也有

此时x方向的损伤为

显然，公式(12、13)还必须有，当a_p<1,a′_p<1时，令a_p＝1,a′_p＝1.2.2 损伤演化模型 在损伤演化模型中，选择mazars一维损伤演化模型^[12]描述受拉方向x的“初始损伤”d₀.在单轴受拉状态时mazars一维损伤演化方程，

在单轴受压状态时mazars一维损伤演化方程，

所以，考虑损伤传递效应的损伤演化方程定义为，在受拉方向，И

在受压方向，

　　式中：a_t b_t a_c b_c为损伤阈值；当a_p<1，a′_p<1时，令ap=1，a′p=1，公式(16a、16b)退化为公式(15a、15b).

共2页: 1