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

对塔式起重机臂三维裂纹的有限元断裂分析 曹宗杰,匡振邦和李昌平 中国上海交通大学工程力学系， 200240长春航空空军大学机械工程系130022 zjcao96@,zbkiang@,Lichangpingl@ 摘要：本文介绍了在三维裂纹前端位移域的局部解和构造三维裂纹问题的一种，和提出了一种新方法用于计算三维裂纹问题的应力强度因子。现有的方法，结构强度与纵向闭合裂纹轴向张力进行了分析。裂纹扩展是可能的，因为轴压塔式起重机主臂的关节附近有环拉伸应力。应力和应力强度因子塔式起重机主臂关节计算。探讨了焊接残余应力和应力强度因子的影响。 关键词：应力强度因子；残余应力奇异有限元；裂纹1引言 电力系统塔式起重机塔式起重机类型最重的。一些灾难性事故发生时塔式起重机在使用工程中的。磁测量，封闭式纵向裂缝在塔式起重机主臂的表面常发现。裂缝影响结构的强度，因此纵向闭合裂纹塔机主臂的强度必须进行分析。因为裂纹尖端一般方法很难准确地确定裂纹如何影响塔式起重机结构强度。许多研究人员在裂纹尖端的奇异领域提出了各种模拟和处理。每一种方法都有其优点和缺点存在。他们中的大多数需要密集的网格计算的有限元方法。一些更为精确的元素，如杂交/混合有限元，有更精确的结果，但很难控制结果的稳定性。因此需要在工程上应用一种特殊的算法。 奇异有限元求解应力强度因子裂纹问题特殊优势。在本文中的三维裂纹的前端位移域解决方案，介绍了一种构建二维裂纹弹性问题和三维裂纹问题的新的奇异单元，提出了一种新方法用于计算应力强度因子三维裂纹问题。此方法不仅可以节省内存和计算时间，也增加了对涉及到应力强度因子有限域问题的计算解的精度。与现有的方法，纵向闭合裂纹结构强度轴向张力进行了分析。塔起重机的主臂关节应力和应力强度因子计算探讨了焊接残余应力和应力强度因子的影响。 2奇异有限元方法 由于应力裂纹奇异性的存在，常规有限元应力域精确的表达强度是非常困难的。因此，任意裂纹形状的三维物体，局部位移场，应力场和应变场可以表示为如下形式： （2） （3） 这里的、和分布函数矩阵的局部应力场，应变场和位移场。的应力强度因子（SIF）的阵列在裂纹前沿的任意一点，二次插值函数模拟元素中裂纹前沿的应力强度因子是该点放置在裂纹前端的功能。然后下面的公式可给出： （5） （6） 其中是应力强度因子的二次插值函数矩阵；是应力强度因子在单元节点的阵列，在裂纹前沿的应力强度因子的节点数组。局部解（）-（3）为等参单元，对任意形状的三维裂纹单元位移场如下： 其中q是节点位移，并用奇异的元素数。M是的节点单元边界，L是元素边界节点。根据应变与位移之间的关系，为，可以得到以下方程 其中和应变向量应力矢量；是微分算子；B是应变矩阵；是奇异应变矩阵；D是弹性系数矩阵。根据应力和应变之间的关系，可以得到以下方程： 三维裂纹应力强度因子根据最小势能原理，可以得到以下方程 是整个研究区，是奇异的。和分别奇异单元。式（7）-（9）代入式（），可以得到以下方程： 其中Q是全球位移阵列。F和分别整体研究和奇异区域边界力向量。K是整体刚度矩阵。是耦合的整体刚度矩阵元素之间的定期和奇异的元素。是奇异单元的整体刚度矩阵。通过分别改变Q和之间的独立性，可以得到以下方程 (12) (13) 当时，式（13）可以改写为下列形式： 替代式（14）代入式（12），可以得到以下方程： (15) 当 (16) (17) 4数值例子 塔式起重机是桁架结构，其主是无缝钢管（）20钢材料。裂纹的长度5～500mm裂缝的深度0.5～2.5。计算模型的建立应用二十节点三维有限元素。关节臂有限元网格与附近的封闭式纵向裂缝的主要危险如图1所示。裂缝计算长度为500mm。分别计算裂缝深度，1.5,2.5，3和3.5毫米。计算负荷（）等于。焊接残余应力对主臂与封闭式纵向流的考虑，焊接残余应（）相当于等于0.5。泊松分布，材料的杨氏模量的比分别等于0.3和0.07 /。环拉伸应力在不同的径向截面无焊接残余应力，如图2所示。环拉伸应力不同轴向截面无焊接残余应力图3和图4所示。通过上述计算，塔主臂的关节部分有环的拉伸应力在可能的流环附近传播。从图可知，主臂纵向闭合裂纹的环拉伸应力降低时有焊接残余应力。在不同的情况下，主臂关节部分应力强度因子的三维裂纹附近无焊接残余应力和焊接残余应力进行了计算，图所示。从图可以看出，最大应力强度因子等