应力强度因子.docx

应力强度因子 断裂与损伤力学 应力强度因子 数值计算方法综述 2013年6月 第一章 应力强度因子求解方法概述 含有裂纹的工程结构的断裂力学分析一直是一个重要问题，在断裂力学理 论中应力强度因子是线弹性断裂力学中最重要的参量。它是由构件的尺寸、形 状和所受的载荷形式而确定。由于裂尖应力场强度取决于应力强度因子，因此 在计算各种构件或试件的应力强度因子是线弹性断裂力学的一项重要任务。 由于应力强度因子在裂纹体分析中的中心地位，它的求解自断裂力学问世 以来就受到了高度的重视。迄今为止，已经产生了众多的理论和致值解法。 70 年代中期以前的有关工作在文献中已有相当全面的总结，近 20 年来，求解的方 法又得刭了明显的发展与完善。 下文将穿透裂纹问题 (二维 )与部分穿透裂纹问题 （三维）分开讨论。 第二章 二维裂纹问题 2.1 复变函数法 由 Muskhelishvili 的复变函数法，应力函数为： Re[ z ( z) (z)] 1/ 2[ z ( z) z ( z) ( z) _] 平面应变情况下的应力与位移为： 2 2 x y x2 y2 4 Re[ ( z)] x y i 2 xy 2[ z (z) ( z)] u iv 3 4 (z) 1 [ z ( z) (x)] 2 2 可以证明，在裂纹尖端区域： K K I  iK  2 2 lim[  z z0  ( z)] 由上式可见。由于  k 仅与  (z) 有关，因此只需确定一个解析函数  ( z)  ，就 能求得 k I ，这一方法一般只能用来解无限体裂纹问题。对于含孔边裂纹的无限 大板，通常可利用复变函数的保角映射原理来简化解题过程。如采用复变（解 析）变分方法，则可求解具有复杂几何形状的含裂纹有限大板的应力强度因子。 2.2 积分方程法 弹性边值问题可以变为求解下列形式的积分方程： M (r ,t ).P(t) dt f ( r ) (t a)(b t ) 由积分方程解出沿裂纹的坐标的函数， 便能直接求出应力强度因子 k 。这个 积分方程在有些特殊情况下可用普通的 Gauss-Chebyshellr 积分或它的修正形式 来求解。 2.3 边界配置法 边界配置法是求解各类边值问题的一种半解析半数值方法。用应力函数法求解二维裂纹问题，关键是选择合适的满足全部边界条件的双调和应力函数，而对有限体或裂纹分布较复杂的情况，封闭形式的应力函数是很难选取的。边界配置法克服了这一困难，它的基本思路是选择以级数展开形式的函数作为满足双调和方程和裂纹面边界条件的应力函数，通过边界条件来确定含有限项的级数中的待定系数。这些待定系数可以通过求解满足边界上的应力，载葡或位移的一组线性代数方程而确定。求解中可以在指定点上精确地满足，也可以在最小均方差的意义上满足边界条件。这样得到的级致解一般能精确满足域内的给定条件，并且近似地满足其余边界上的条件。 在裂纹问题的边界配置法中有两种基本的应力函数可供选择，即 Williams 的应力函数和 Muskhelishyili 的复变应力函数， 从发展过程看， 前者一般用在边缘裂纹问题中，后者可用于内埋裂纹与边缘裂纹的情况。 边界配置法的求解精度较高。它的不足之处是：对于不同类型的裂纹问题，应力函数必须改变。而建立这些新的应力函数的工作量将是很大的，对于较复杂的几何与载荷情况，应力函数所应满足的边界条件很难确定，另外，解的收敛性还没有得到严格的证明。 2.4 边界力法 边界力法通过利用无限体中有限数量的集中力和集中力矩的叠加来求解边 值问题。这种解法以无限体中集中力和集中力矩的弹性解为基本解，对于不含 裂纹的板，基本解取 Muskhelishyili 的解，对于含裂纹的板， 则取 Erdogan 的解 作为基本解。由于 Erdogan 的解精确地满足了裂纹面应力为零的条件，所以裂 纹面就不再需要作为边界的一一部分加以考虑。 因为基本解满足了物体内部的所有弹性力学方程，余下所需满足的条件只 是边界条件。这些边界条件则是通过在相应于真实裂纹体的假想边界上施加一 系列的集中力和集中力矩来满足的，先把假想的边界离散化为一组线段，在每 一段的中心，在离开假想边界处加上一对集中力和力矩，这些力和力矩的值可 通过近似地满足边界条件得以确定。 与其他数值方法相比，边界力法有其明显的优点。由于这一方法已精确地 满足了裂纹面上的边界条件，所以它不需要像边界元法那样把裂纹面视为边界 的一部分。另外，它也克服了边界配位法中所需要的对每一类裂纹问题都要建 立新的应力函数的缺点。这种解法只要较小的自由度就能达到相当高的精度。 因此它在求解几何形体复杂的裂纹向题中有着明显的优点，但在处理复杂载荷 的能力方面，则远非如权函数法那样灵活。 2.5