多裂纹扩展分析的边界元方法Ξ.PDFVIP

Vol. 23 No. 4 　　 　　　　　　 　 　　　　　　　　　 　 第23卷 第4期 固 体 力 学 学 报 2002年 　12月 ACTA MECHANICA SOLIDA SINICA December 2002 多裂纹扩展分析的边界元方法 秦 　飞 (北京工业大学机械电子工程学院 ,北京 ,100022) 岑章志 　　杜庆华 (清华大学工程力学系 ,北京 ,100084) 摘 　要 　采用边界元数值模拟和即时等效材料常数计算相结合的办法 ,只需模拟一个裂纹 的扩展情况便可预测裂纹体的整体响应. 针对准脆性材料的特点 ,采用粘性裂纹模型模拟裂纹 开裂行为 ;采用二次裂纹扩展量作为增量控制变量 ,避免了软化及失稳分析中用力或位移作控 制变量时遇到的困难. 针对二次裂纹扩展路径未知情况 ,给出了预测和修正裂纹开裂界面的迭 代技术. 分析计算了含多个规则分布裂纹石膏板受压时的响应 ,并与实验结果进行了比较. 结果 表明该文方法的可行性与有效性. 关键词 　边界元 ,等效材料常数 ,多裂纹 ,扩展 ,粘性裂纹模型 1 　引言 岩石 、混凝土类材料由于在形成或制造过程中受到机械或环境载荷作用 ,在它们被使用 时已含有裂纹等缺陷. 这些裂纹可能是随机分布 ,也可能规则分布. 本文中称含这样裂纹的 固体为裂纹体. 当裂纹体受压时 ,初始裂纹面间产生摩擦滑动导致二次裂纹在各个初始裂纹 的两端萌生并随载荷增加而扩展. 二次裂纹的扩展方向一般平行于最大受压方向 ,但会受其 它裂纹和试件边界影响. 二次裂纹是张开裂纹 ,是岩石类材料“容胀效应”的主要原因. 二次 裂纹扩展的方向选择性导致裂纹体响应的各向异性[1 ] . ( 为了更深入理解裂纹体的力学行为 ,较合理的研究方法应该从分析裂纹行为 闭合、摩 ) 擦滑动、扩展等 开始 ,预测裂纹体整体响应. 这样 ,裂纹扩展与整体响应之间的关系便一目 了然. 本文即采用这样的研究路线. 但首先必须解决两个关键问题 :一是裂纹扩展的数值模 拟 ;二是如何同时模拟数目较多的裂纹. 数值计算是模拟裂纹扩展的有效方法. 和有限元方法相比 ,边界元在模拟裂纹问题时有 明显优势[2 ,3 ] ,尤其对裂纹扩展问题. 文献[4 ]采用边界元分区算法 ,用二次裂纹扩展量作为 控制变量成功分析了混凝土构件中裂纹扩展导致的软化现象 ,并采用粘性裂纹模型描述混 凝土材料中裂纹行为特征. 但在他们的分析中 ,只有两个裂纹. 当裂纹数目较多时 , 由于沿裂 纹面分区非常困难 ,使得一般分区算法不再适用. 为解决多裂纹问题 ,作者提出了新的分析方法. 该方法将边界元分区算法与等效材料常 数方法相结合 ,利用等效材料常数计算避开了分析多裂纹的困难 , 同时利用边界元分区算法 模拟裂纹扩展问题的优势对单个裂纹进行扩展模拟. 具体做法是 ,选取裂纹体中一代表裂 收到第 1 稿 收到修改稿. © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 　　　　　　　　　　　　　　　　　固体力学学报 　　　　　　　　　　　　　　　2002 年 第 23 卷 4·32 · 纹 ,采用边界元分区算法对其扩展过程进行模拟 , 同时将其它裂纹的影响通过改变材料弹性