应力强度因子.doc

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