钢筋混凝土简支梁结构损伤破坏的有限元分析.docxVIP

钢筋混凝土简支梁结构损伤破坏的有限元分析 在外部负荷和环境的作用下，由于微裂纹、微孔等材料和结构的缺陷，材料和结构的质量损失是众所周知的。混凝土材料是由水泥、粗、薄材料和各种混合材料组成的多相复合材料。在形成之初，有许多微裂纹和微孔。在外部负荷的作用下，这些微裂纹继续扩展和积累，并形成了一些宏观裂纹。随着负荷的增加，宏观单元格的聚集和扩展最终导致材料失败。在附近的破坏阶段，整个结构具有明显的非线性特征。这种非线性行为主要是由于混凝土材料的破坏引起的。 材料损伤的概念最早由Kachanov在1958年研究金属蠕变破坏时提出,它是以材料有效承载面积的缩减来定义材料的损伤度,并由此引出有效应力的概念.1971年,Lemaitre提出了著名的应变等效假设,损伤材料的应变本构方程与无损材料相同,只要将Cauchy应力用有效应力代替即可.基于Kachanov型损伤理论,国内外不少学者建立了多种不同的混凝土损伤本构模型[5,6,7,8,9,10,11,12,13].对于Kachanov型损伤理论缺陷的分析讨论见文献.另外,由于有效应力定义的随意性与不对称性,将Kachanov的单轴损伤模型向多轴和各向异性情况推广一直没有取得成功.文献直接从不可逆热力学基本定律出发,建立了弹性损伤问题的普遍理论,解决了将Kachanov损伤模型向各向异性推广时所遇见的理论困难.混凝土材料在受拉与受压状态下损伤行为明显不同,这是混凝土材料损伤破坏的固有特性.文献中的研究对象是拱结构,拱的受力特点主要是轴向受压.而对于梁的弯曲损伤破坏问题,试验观测表明,受拉区与受压区的损伤破坏情况明显不同,这表明分析梁的损伤问题时必须考虑受拉与受压状态的区别.因此,本研究引入一个拉压修正系数γ,对文献中损伤本构模型进行改进,使它适用于梁的弯曲损伤问题. 1 应变张量的确定 对于各向同性弹性损伤问题,根据热力学第二定律,损伤本构方程与损伤应变能释放率Y的表达式为: σij=2μεijΜμ(D)+λεkkδijΜλ(D)；(1)Y=-12λΜ′λ(D)(εkk)2-μΜ′μ(D)εijεij.(2) 式中:σij,εij分别为Cauchy应力张量与应变张量;λ,μ为无损材料的Lame弹性常数;D是各向同性标量损伤变量;Y为标量损伤对偶变量;δij为Kronecker张量;Mλ(D),Mμ(D)是两个损伤效应函数,分别表示损伤对材料两个弹性常数的影响;撇号表示求导. 材料的损伤演化与其损伤耗散功率密切相关,可设损伤速率˙D正比于Y的m次方,损伤演化方程的一般形式是: ˙D=AYm=A[-12λΜ′λ(D)(εkk)2-μΜ′μ(D)εijεij]m.(3) 式中:A与m为与材料损伤演化性质有关的参数,为正值. 式(1)～(3)的级数形式为: δij=2μεij[1-Ν∑n=1β(n)Dn]+λεkkδij[1-Ν∑n=1α(n)Dn]；(4)Y=12λΝ∑n=1α(n)nDn-1(εkk)2+μΝ∑n=1β(n)nDn-1εijεij；(5)˙D=A[μεijεijΝ∑n=1nβ(n)Dn-1+12λ(εkk)2Ν∑n=1nα(n)Dn-1]m.(6) 式中:α(n),β(n)为无因次非负材料损伤参数. 考虑到损伤变量不是无穷小量,而损伤速率应与材料当前损伤状态有关,故取式(4)～(6)的二阶近似形式(即N=2).取α(n)=β(n),这意味着Poisson比不随损伤而变化.实际上,剪切应变与体积应变是材料的两种基本变形形式,前者反映材料单元的形状变化,后者反映材料单元的体积变化.将应变张量分解成体积应变与剪切应变两部分考虑,可得 σij=(1-α1D-α2D2)(2μεij+λεkkδij)；(7)Y=(α1+2α2D)[32με2eq+(3μ+92λ)(εm)2]m；(8)˙D=AYm=A(α1+2α2D)m[32με2eq+(3μ+92λ)(εm)2]m.(9) 式中:εm=εkk/3;εeq=√23eijeij;eij为应变偏量;εeq为等效应变;εm为体积应变,反映了材料单元体积变化对损伤的影响;εeq为von Mises等效应变,反映了材料单元形状变化对损伤的影响;α1与α2为二个非负材料参数,且α1+α2=1.当D=1时,材料完全失效,材料的有效弹性常数应降为0.因此混凝土损伤本构模型含有三个独立的材料损伤参数,即α1(或α2),A与m.对于不同的材料,这些材料参数应有不同的数值,须由试验结果来确定,这属于结构分析的反问题,即模型参数识别问题,其解的存在性与唯一性问题是尚未得到解决的一个理论难题.本研究是在相当大的范围内对三个损伤参数进行打靶计算,对于每一种可能的组合,计算出破坏荷载下的损伤场,选取最接近试验破坏形式的组合,确定出三个参数的具体数值. 混凝土材料的性能