钢筋混凝土结构有限元分析的力学模型研究2006防灾减灾及防护工程.docVIP

钢筋混凝土结构有限元分析的力学模型研究 2006 防灾减灾及防护工程 曹明 摘要：碾压混凝土的正交异性损伤本构模型研究根据碾压混凝土拱坝的结构特性与碾压混凝土材料的破坏特征，以Sidoroff各向异性损伤理论为基础，建立了碾压混凝土正交异性损伤本构模型；基于正交异性损伤模型的损伤传递方式，导出了碾压混凝土材料的损伤演化方程。以碾压混凝土在复杂应力下的试验资料，对碾压混凝土正交异性损伤本构模型进行了验证。 关键词：碾压混凝土 混凝土 正交异性损伤 本构模型 碾压混凝土的正交异性损伤本构模型研究问题目前计算分析规律方法研究 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-)/A；A为无损状态时的横截面积；为受损后有效承载面积；A′0为初始损伤面积，A′0＝A-.这里，材料的微缺陷与内部裂纹描述为图1所示模型：材料中存在许多的微小裂纹与微小缺陷，这些微小裂纹与微小缺陷对某一方向X的损伤可等效为许多的微小圆形裂开区域A′0i，使得这一向上的有效承载面积减小。材料在X方向的初始损伤面积A′0＝∑A′0i。 　　假定材料在Z方向受到压应变，考虑压应变对损伤区域A′0i的作用。按照线弹性断裂力学分析，在一个含有椭圆裂纹的无限介质中，当受到单向压缩时，在平行于加载方向的裂纹尖端将产生切向拉应力，使裂纹扩展[3]。这里，按照应变作用方式描述裂纹的扩展，如图2所示，在Z方向压应变作用下，损伤面张开。当裂纹张开达到一定的宽度，裂纹便发生新的扩展，使X方向有效承载面积减小，造成新的损伤。 图1 X方向的材料微缺陷与内部裂纹模型 图2 压应变作用下的裂纹张开