XFEM在高压管道裂纹扩展计算中应用.docVIP

XFEM在高压管道裂纹扩展计算中应用 　　摘要： 介绍XFEM的基本理论，针对常规有限元计算困难的三维裂纹，运用XFEM完成三维高压管道中的裂纹扩展计算.根据管道裂纹扩展的计算结果，总结得到管道裂纹扩展的规律，得到径向裂纹比轴向裂纹更严重的结论. 　　关键词： 高压管道； 疲劳； 裂纹扩展； 三维裂纹； 有限元 　　中图分类号： U4文献标志码： B 　　Abstract： The basic theory of XFEM is introduced. As to the problem that it is difficult to calculate 3D crack by the normal finite element， XFEM is used to perform the crack propagation calculation on a 3D highpressure pipe. According to the crack propagation calculation results of the pipe， the law of pipe crack propagation is summarized and it is concluded that the crack in radial direction is more critical than that in axial direction. 　　Key words： highpressure pipe； fatigue； crack propagation； 3D crack； finite element 　　收稿日期： 2014[KG*9〗05[KG*9〗10修回日期： 2014[KG*9〗10[KG*9〗24 　　作者简介： 陈星文（1983―），男，山东青岛人，工程师，硕士，研究方向为管道力学，（Email）chxw@0引言 　　在用有限元法处理工程问题时，由于常规有限元法的位移函数是连续函数，所以利用常规有限元法所建立的单元其内部位移必须是连续的.因此，常规有限元法在处理如裂纹这类不连续问题时，必须把裂纹面作为单元边界面，裂纹尖端作为为单元的节点. 　　在裂纹尖端附近区域，应力集中使其应力急剧变化.常规有限元法的位移函数是多项式，在裂纹尖端的网格划分中，若多项式的次数不变，要模拟裂纹尖端附近应力的急剧变化，使有限元解逼近于精确解，则必须将单元尺寸划分得非常小.另外，如果单元尺寸不变，可通过增加多项式的次数来反映裂纹尖端单元的应力变化.实践证明，在处理裂纹尖端应力集中这类问题时，高次多项式的逼近性能非常差，求解结果不能满足精度要求.[1] 　　在模拟裂纹扩展问题上，常规有限元显得更加无能为力，由于必须把裂纹面设置为单元界面，裂纹尖端设置成单元节点[12]，因此当裂纹扩展后，就必须对结构重新划分网格，而由于裂纹尖端存在应力集中，网格必须细化，所以会极大地增加计算量，对于一些复杂结构问题甚至无法实现. 　　由于常规有限元法在计算裂纹问题上存在上述缺陷，1996年德克萨斯大学教授MELENK等[3]提出单位分解概念，通过单位分解，扩大有限元形函数的形式，使得不连续函数和非多项式函数可以添加到有限元的位移函数中，使得许多在相关问题上逼近性很强的函数可以用于有限元单元位移函数，从而使有限元的计算变得简洁，结果更加精确.1999年，以BELYTSCHKO为首的美国西北大学力学研究组提出XFEM理论：在单位分解法的基础上改变常规有限元的位移函数.[45]其在常规有限元位移函数的基础上加上裂纹尖端位移函数因子，[68]从而无须将网格在裂纹尖端细化就可以准确地逼近结构的精确解.而且，XFEM不需要将结构内部存在的几何或物理界面作为网格划分的边界[9]，其所使用的网格可以与结构内部的形状特点无关.这样，在裂纹扩展的过程中就不需要对结构重新进行网格划分，为裂纹扩展计算提供可行性和便捷性. 　　1XFEM基本理论 　　1.1裂纹界面位移改进函数 　　在常规有限元中，单元位移函数通常表示为u=Niui 　　v=Nivi 　　w=Niwi（1）式中：ui，vi和wi为节点i的位移；Ni为结构的形函数. 　　由于常规有限元的单元位移函数是连续函数，因此对结构进行网格划分时裂纹面必须作为单元的边界，以保证有限元单元内部位移的连续性.对于带裂纹的平板（图1），按常规有限元方法，在不考虑裂纹尖端应力集中的情况下，网格划分见图2. 　　按图2的网格，结构内任意点的位移可表示为u=10i=1Niui（2）令a=u9+u102，b=u10-u92，N11=N9+N10，u11=a，代入式（2），得u=8i=1Niui+N11u11+bN1