外文翻译--对塔式起重机臂三维裂纹的有限元断裂分析 中文版

第 38 卷增刊 上海交通大学学报 Vol.38.sup 2004 年 6 月 JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY Jun.2004 对塔式起重机臂三维裂纹的有限元断裂分析 曹宗杰 ,匡振邦和李昌平 中国上海交通大学工程力学系， 200240 中国 长春航空空军大学机械工程系 ， 130022 ,,L 摘要 ：本文介绍了在三维裂纹前端位移域的局部解和构造三维裂纹问题的一种 新方法 ，和提出了一种新方法用于计算三维裂纹问题的应力强度因子。 用 现有的方法， 对 结构强度与纵向闭合裂纹轴向张力进行了分析。裂纹扩展是可能的，因为轴压塔式起重机主臂的 关节附近有环拉伸应力。应力和应力强度因子 用于 塔式起重机主臂关节计算。探讨了焊接残余应力和应力强度因子的影响。 关键词：应力强度因子；残余应力奇异有限元；裂纹 1 引言 电力系统塔式起重机 是 塔式起重机类型 中 最重的 的一个 。一些灾难性事故发生时塔式起重机 都是 在使用工程中的。磁测量，封闭式纵向裂缝在塔式起重机主臂的表面 经 常发现。裂缝影响结构的强度，因此 对 纵向闭合裂纹塔 式起重 机主臂的强度必须进行分析。因为裂纹尖端 ， 一般方法很难准确地确定裂纹如何影响塔式起重机结构强度。许多研究人员在裂纹尖端的奇异领域提出了各种模拟和处理16（ ）方 法 。每一种方法都有其优点和缺点存在。他们中的大多数 是 需要密集的网格计算的有限元方法。一些更为精确的元素，如杂交 /混合有限元，有更精确的结果，但 是它 很难控制结果的稳定性。因此需要在工程上应用一种特殊的算法。 奇异有限元求解应力强度因子裂纹问题 时有 特殊 的 优势。在本文中的三维裂纹的前端位移域解决方案，介绍了一种构建二维裂纹弹性问题和三维裂纹问题的新的奇异单元， 并 提出了一种新方法用于计算应力强度因子 的 三维裂纹问题。此方法不仅可以节省内存和计算时间，也增加了对涉及到 应力强度因子有限域问题的计算解的精度。与现有的方法 相比 ， 对 纵向闭合裂纹 的 结构强度轴向张力进行了分析。 对 塔起重机的主臂关节应力和应力强度因子 进行 计算 。 探讨了焊接残余应力和应力强度因子的影响。 第 38 卷增刊 上海交通大学学报 Vol.38.sup 2004 年 6 月 JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY Jun.2004 2 奇异有限元方法 由于应力裂纹奇异性的存在，常规有限元应力域精确的表达强度是非常困难的。因此，任意裂纹形状的三维物体，局部位移场，应力场和应变场可以表示为如下形式： l l lSK （ 1） l l lBK（ 2） l l lMK （ 3） 这 里的lS、lB和LM分别是 分布函数矩阵的局部应力场，应变场和位移场。lK的应力强度因子（ SIF）的阵列在裂纹前沿的任意一点，二次插值函数模拟元素中裂纹 前沿的应力强度因子是该点放置在裂纹前端的功能。然后下面的公式可给出： lKL（ 4） 12, Tn （ 5） , Ti il ill illlK K K （ 6） 其中 L 是应力强度因子的二次插值函数矩阵； 是应力强度因子在单元节点的阵列，i是 在裂纹前沿的应力强度因子的节点数组。局部解 式 （ 1） -（ 3）为等参单元，对任意形状的三维裂纹单元位移场如下： ( ) ；llN q M L N M L j mN q j m （ 7） 其中 q 是节点位移，并用奇异的元素数。 M 是 M 的节点单元边界， L 是元素边界节点。根据应变与位移之间的关系，为 ( )( )TEu ，可以得到以下方程 ；kB q B j mB q j m （ 8） 其中 和 为 应变向量应力矢量； ()TE 是微分算子； B 是应变矩阵；kB是奇异应变矩阵； D 是弹性系数矩阵。根据应力和应变之间的关系，可以得到以下方程： 第 38 卷增刊 上海交通大学学报 Vol.38.sup 2004 年 6 月 JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY Jun.2004 _；kD B q D B j mD B q j m （ 9） 3 三维裂纹应力强度因子 根据最小势能原理，可以得到以下方程 211()2 eNe T T Tej f d u f d （ 10） 当 是 整个研究区 时 ，2e是奇异的 区 。 f 和 f分别 是表面力的 奇异单元。将 式（ 7） -（ 9）代入式（ 10），可以得到以下方程： 1122T T T T Ts N sQ K Q Q F K Q K F (11) 其中 Q 是全球位移阵列。 F 和 F分别 是 整体 区域 研究和奇异区域边界力向量。 K 是整体刚度矩阵。NsK是耦合的整体刚度矩阵元素之间的定期和奇异的元素。SK是奇异单元的整体刚度矩阵。通过分别改变 Q 和 之间的独立性，可以得到以下方程 0NsK Q K F (12) 0s N sK Q K F (13) 当 TsN NsKK时，式（ 13） 可以改写为下列形式： 1 ()s N sK K Q F (14) 替代式（ 14）代入式（ 12），可以得到以下方程： 0KQ F (15) 当 1N s s s NK K K K K(16) 1N s sF F K K F (17) 4 数值例子 塔式起重机是桁架结构，其主 臂 是无缝钢管（ 219x20 ） 是 20 钢材料。裂纹的长度 为 5 500mm， 裂缝的深度 为 0.5 2.5mm。计算模型的建立应用二十节点三维有限元素。关节臂有限元网格与附近的封闭式纵向裂缝的主要危险如图 1第 38 卷增刊 上海交通大学学报 Vol.38.sup 2004 年 6 月 JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY Jun.2004 所示。裂缝计算长度为 500mm。分别计算裂缝深度， 为 1.5,2.5， 3 和 3.5 毫米。计算负荷（p）等于21700kgcm。 出于 焊接残余应力对主臂与封闭式纵向流的考虑，焊接残余应 力 （0）相当于等于 0.5 7p 。泊松分布，材料的杨氏模量的比分别等于 0.3 和 0.07 MN / 2m 。环拉伸应力在不同的径向截面无焊接残余应力，如图 2 所示。环拉伸应力不同轴向截面无焊接残余应力 如 图 3 和图 4 所示。通过上述计算，塔 式起重机 主臂的关节部分有环的拉伸应力在可能的流环附近传播。 从图 24 可知，主臂纵向闭合裂纹的环拉伸应力降低时有焊接残余应力。在不同的情况下，主臂关节部分应力强度因子的三维 裂纹附近无焊接残余应力和 对焊接残余应力进行了计算， 如 图 5,图 6 所示。从图 5 可以看出，最大应力强度因子等于3 21.9 MNm时裂纹的厚度为 3.5mm，关节 处 发生较少 ， 在 10cm 左右 。单边缘裂纹扩展门槛值 ， 它的长度是 0.025 到 5 毫米 ， 是从 4.0 至 3 26.4 MN m 的20 钢材料。主臂关节附近的最大应力强度因子小于单边缘裂纹扩展门槛值。 如图 6，最大应力强度因子等于 3 20.9 MN m 时，裂纹深度等于 3.5mm。最大应力强度因子与焊接残余应力小于无焊接残余应力。 第 38 卷增刊 上海交通大学学报 Vol.38.sup 2004 年 6 月 JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY Jun.2004 5 结论 本文 对 三维裂纹前端位移域的解决方案进行了介绍和三维裂纹问题的一种新的奇异单元的构造 ， 和提出了一种新方法 ， 计算三维裂纹问题的应力强度因子。第 38 卷增刊 上海交通大学学报 Vol.38.sup 2004 年 6 月 JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY Jun.2004 与 用 现有的方法， 对 结构强度与纵向闭合裂纹轴向张力进行了分析。这是可能的裂纹扩展，因为 是 有环的拉伸应力下的轴压塔式起重机主臂的关节附近的应力。焊接残余应力对焊接材料强度有害，但焊接残余应力使主臂环拉伸应力下 降 。在同一时间，应力强度因子的下降， 改进了 塔式起重机主臂的安全性。 参考文献 1 王 志超，张理苏 , 一种新的准谐调有限元列式方法及在断裂力学中的应用 1990.37（ 6） 1195-1201 第 38 卷增刊 上海交通大学学报 Vol.38.sup 2004 年 6 月 JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY Jun.2004 2 O.C.钦科维奇（ Zienkiewicz O.C.） ， 有限元法 ,麦格劳希尔图书有限公司第三版： 669， 1977 3Barsoum RS.Int J.Numer.Merh.Engng.1976,10:2537 4 Abdel Wahab MM, De Roeck G A Finite Element Solution for Elliptical Cracks Using the ICCI Method Engineering Fracture Mechanics, 1996,53:518-526 5 曹宗杰 三维奇异准谐调单元理论研究及其数值分析 ， 1994。