求解应力强度因子的样条虚边界元交替法.docVIP

求解应力强度因子的样条虚边界元交替法 　　摘要： 　　为求解平面裂纹问题的应力强度因子，提出基于Muskhelishvili基本解和样条虚边界元法的样条虚边界元交替法.该方法将平面内带裂纹有限域问题分解成带裂纹无限域问题与不带裂纹有限域问题的叠加.带裂纹无限域问题利用Muskhelishvili基本解法直接得出，不带裂纹有限域问题采用样条虚边界元法求解.利用该方法对复合型中心裂纹方板和I型偏心裂纹矩形板进行分析.数值结果表明该方法精度高且适用性强. 　　关键词： 　　平面裂纹； 样条虚边界元法； 交替法； 断裂力学； Muskhelishvili基本解； 应力强度因子 　　中图分类号： O302 　　文献标志码： A 　　Abstract： 　　To solve the stress intensity factor of plane crack problem， a spline fictitious boundary element alternating method is proposed based on Muskhelishvili fundamental solution and spline fictitious boundary element method. A finite field crack problem is transformed into the superposition of a simple finite field problem without crack and an infinite problem with crack. The Muskhelishvili fundamental solution method is used to solve the infinite problem with crack and the spline fictitious boundary element method is implemented to solve the simple finite field problem without crack. A square plate with a slant center crack and a Itype rectangular plate with an eccentric crack are solved using the method. The numerical results show that the method is of high accuracy and strong applicability. 　　Key words： 　　plane crack； spline fictitious boundary element method； alternating method； fracture mechanics； Muskhelishvili fundamental solution； stress intensity factor 　　0引言 　　任何材料和工程结构都会不同程度地存在裂纹缺陷，其产生和扩展对构件的承载能力会造成很大程度的破坏，因此断裂力学在现代强度理论中的地位越来越重要.应力强度因子是表征裂纹特性的重要参量，所以对其计算是断裂力学研究的重要环节. 　　在现阶段，求解应力强度因子的方法主要为改进的有限元法，包括奇异有限元法和扩展有限元法.[12]奇异有限元法通过移动节点使奇异点出现在1/4处而不是中点处，使得边界节点处出现奇异的应力场，但是该方法存在单元直接协调性和计算收敛性的问题.扩展有限元法改进单元的形函数，使之包含不连续性的基本成分，从而放松对网格密度的划分要求，但是其刚度矩阵存在病态问题并增加许多额外的未知量. 　　为更加高效精确地分析裂纹问题，有学者提出求解裂纹问题应力强度因子的SchwartzNeumann交替法.该方法将带裂纹的复杂结构分解成为一个不含裂纹的复杂结构与一个含裂纹的无限大域，运用迭代方法或者线性方程对分解后的结构进行求解.含裂纹无限域采用Muskhelishvili基本解[3]求解，能够直接利用表达式求解出平面内任意一点的响应和裂纹尖端的应力强度因子，因此具有精度高和计算量小的优点.不含裂纹的复杂结构可以采用数值方法求解，有限元法是最为常见的数值方法[410]，将交替法与有限元法相结合的方法称为有限元交替法.然而，有限元法的应力结果相较于位移结果来说精度较低，在循环迭代的过程中会造成误差的进一步增大. 　　样条虚边界元法是一种高效的间接边界元法，其只需要对边界进行划分，可降低问题求解的维度，使得计算效率大大提高，目前已经在工程