有限元法在应力强度因子计算中研究与应用.docVIP

有限元法在应力强度因子计算中的应用 马源 myeric90@126.com（原创） 摘要：本文构建了含裂纹平板的二位1/4模型及三维1/2模型，分别用于计算张开型、滑移型和撕开型裂纹尖端应力强度因子。有限元分析的结果在误差范围内可以较好地与解析解吻合。计算结果表明，裂纹周向单元的划分会严重影响有限元计算的结果，三维模型沿厚度方向提高划分密度可以有效提高计算精度。 关键词：有限元 应力强度因子 精度分析 1前言 工程分析中，材料中的裂纹会对结构可靠性带来很大的影响。历史上有很多航空航天事故、建筑事故都是由于裂纹引起的断裂导致结构失效的。为了检验结构是否能够一般用于判断裂纹是否延伸的重要判据就是应力强度因子K（Stress Intensity Factor，SIF）。对于任何材料，其应力强度因子极限KC只与材料本身的属性有关，而和裂纹尺寸、裂纹周围应力强度无关。在具体的工程分析中，评估含裂纹结构稳定性，只需要计算含裂纹结构在要求的工况下的裂纹尖端应力强度因子K值，若KKc，则裂纹会发生扩展，导致结构失效。 具体工况下，应力强度因子K的计算可以通过多种方法完成。弹性力学给出了三种基本断裂模式（分别为张开型、滑移型和撕开型，见图1）的应力强度因子解析解。但是对于一般几何结构而言，求解析解的复杂程度会随结构的复杂程度成倍增加，很多情况甚至无法求出解析解。有限单元法及边界元法可以用于具有复杂几