论文部分内容阅读
有限元法是结构断裂分析的常用数值方法之一,随着计算机技术的飞速发展,利用有限元软件求解裂尖应力强度因子日趋成熟,但在求解过程中仍存在对网格离散技术的依赖性。广义参数Williams单元(简记:W单元)能以较少的离散单元数直接确定裂尖应力强度因子,且对奇异区尺寸不敏感。但现有研究成果仅限于裂尖正方形奇异区规则W单元,即裂尖位于奇异区几何中心,沿裂纹面与奇异区角点将奇异区离散为8个全等的扇形条元。基于规则W单元,本论文提出一种裂纹相对独立于整体网格的分析模型,以裂尖所在单元为奇异区并在该区域嵌入W单元。裂尖位置的任意性,使得沿裂尖周向构造的W单元不规则,进而研究了非规则W单元对奇异区形状和尺寸的敏感性以及周向离散数对计算精度的影响。根据裂尖坐标和奇异区四个角点坐标可自动形成非规则W单元,进一步分析了直线裂纹定步长扩展过程的裂尖应力强度因子,为裂纹扩展研究提供了新思路。主要研究内容及结论如下:1、推导了极坐标和直角坐标形式的Williams应力场和位移场,将W单元归纳为一种有限元内嵌特殊元,从理论证明了 W单元不依赖于奇异区形状,同时完善了 W单元求解应力强度因子的计算格式,进而编写了FORTRAN程序。算例分析表明:W单元解与奇异单元解及扩展有限元解吻合较好,且具有较高的计算精度。证明了 W单元理论推导更简单、网格离散更便捷、计算精度更高的优点。2、以平面裂纹问题为基本模型,对裂尖位于菱形、矩形、圆形和椭圆形奇异区中心的规则W单元开展研究,并对奇异区尺寸、a、m、n、N等影响因素进行了分析,建议:m=9,n=300,a=0.9,N=8,圆形奇异区半径S>1/16a,菱形、矩形和椭圆形奇异区建议ojB:ojA>0.4。算例分析表明:W单元的精度依赖于奇异区形状和尺寸。3、建立周向离散单元数为6和5的非规则W单元分析模型,并对其精度敏感性开展研究。给出了网格离散时裂尖相对于奇异区几何中心位置的取值范围:以5%为误差限,当N=6,建议|ζ|<0.72,|η|<0.48;当N=5,ζ∈[-0.56,0.4),η∈[-0.24,0.4],从而保证非规则W单元应用过程中的高精性。算例分析验证了裂尖位置参考范围的正确性。4、基于周向离散单元数为6和5的非规则W单元,建立了裂纹在单元内扩展与裂纹在单元间扩展的两种定步长扩展模式,通过算例模拟了平面裂纹的准静态扩展,验证了 W单元在裂纹扩展应用中的可行性与高精性,根据裂尖位置参考范围可很好地排查计算误差较大的点,并给出了非规则W单元研究裂纹扩展时的相关建议:裂纹在单元内扩展的模式可保持较高求解效率,适应于裂纹扩展趋势研究;裂纹在单元间扩展的模式可保持较高求解精度,适应于有高精度要求的研究。