基于ANSYS的应力强度因子计算.pdf

《基于ANSYS的应力强度因子计算.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、建精世器2010年第31卷第2期基于ANSYS的应力强度因子计算王伟(湖北超高压输变电公司，武汉430050)摘要：以平板表面裂纹为例，介绍和分析了运用有限元软件ANSYS计算应力强度因子的方法。通过对求得的应力强度因子值与解析解的比较，表明用有限元方法计算应力强度因子具有相当高的精度，并且操作简便。关键词：ANSYS；表面裂纹；应力强度因子CalculatetheStressIntensityFactorbytheFiniteElementSoftwareANSYSWANG耽i(HubeiEHVTra

2、nsmission&SubstationCompany，Wuhan430050，China)Abstract：Takingtheplanewithsurfacecrack勰anexample。itconductstheintroductionandanalysisoftheuseoffiniteelementsoftwareANSYScalculationofstressintensityfavorapproach．Itobtainedthroughthestr瞄intensityfactorvalue

3、swiththeanalyticsolutionofthecomparison，showingthatthefiniteelementmethodusingthestre鸽intensityfactorhasaveryhighprecisionandeasytou8e．Keywords：ANSYS；surfacecrack；stre鲳intensityfactor断裂力学是研究带裂纹材料或结构的强度以及裂纹扩展规律的一门学科。在结构有裂纹时，常规强度准则就不再使用。断裂力学提出了应力强度因子的概念，即裂

4、纹端部应力场强度由应力强度因子度量，因此应用断裂力学来解决实际工程问题时，必须首先确定裂纹尖端的应力强度因子。关于各种情况的应力强度因子，国内外有很多学者进行了大量研究工作，我国航空研究院在广泛收集国内外最新研究成果的基础上编制了一本新的应力强度因子手册[1l，然而对于许多不规则的结构形式，手册中并没有涉及。随着计算机技术和有限元理论的发展，使裂纹仿真成为可能，可以通过大型通用有限元软件计算裂纹尖端应力强度因子。该文以带三维表面裂纹平板为例，运用ANSYS有限元分析软件，创建了含三维表面裂纹的平板有限元

5、实体模型，并计算了裂纹尖端的应力强度因子。761裂纹体结构板宽W=100mm。t=12mm的板中有，口=1mm，c=5mm的半椭圆表面裂纹，受拉伸载荷口=600MPa作用。弹性模量E=207GPa，泊松比=0．3121oyt％ffff

6、ffIJIIIlll磊图1含中心表面裂纹平板建祷世界2010年第3l卷第2期2裂纹体建模放大图。AYSYS有限元裂纹体建模方法的主要思路是：将含裂纹构件中包含裂纹部分的小块体积视为局部裂纹体，分别创建整体模型和局部裂纹体模型；对整体模型和局部裂纹体模型进行布尔运算，将含有

7、裂纹部分的局部体积从整体结构中分离出来；再将含图5整个裂纹体裂纹面以外的实体粘接在一起，完成带裂纹体模型的创建；对含裂纹部分进行特殊的单元网格划分，余下的非裂纹体采用常规网格划分，从而生成整个裂纹体有限元模型【4J。具体实现方法如下：1)建立不含裂纹的1／2平板实体模型(见图2)图6扣除裂纹体模型的模型局部放大图6>选取图5所示的裂纹体中包含裂尖点的4块小体积，见图7，采用对裂尖所在表面上的4块面积进行二维网格划分[3，4】(裂纹尖端附近利用AN．SYS前处理命令生成二维的奇异单元)如图8所示。图2不含

8、裂纹的平板模型2)在裂纹尖端附近沿裂纹长度方向创建3块面积(见图3)图7裂尖位置示意图图3椭圆面．3)3块面积组合在一起，同时沿y轴正方向拉伸，生成含半边裂纹面的半裂纹体(见图4)图4半裂纹体4)同步骤2)。创建与之前面积相重合的面积，组合在一起沿y轴负方向拉伸，生成裂纹体的另半个裂纹面和裂纹体，整个裂纹体模型见图5。5)对不含裂纹的实体平板模型和裂纹体模型进行布尔运算，图6为扣除裂纹体模型的模型局部7飞霞震图8裂纹面7)将所生成的面网格沿裂纹深度方向进行俸扫掠网格划分生成三维单元，见图9。图9裂纹体划

9、分网格8)对裂纹体(图4)远离裂纹尖端的两块体积及图6所示扣除裂纹体部分的体积进行自由网格划分，生成带裂纹体有限元模型(见图10)建材世界2010年第31卷第2期图10有限元整体模型3应力强度因子计算对于表面裂纹模型，主要研究裂纹最深处和裂纹表面两点的应力强度因子。采用上节建立的有限元模型，最终计算所得两裂尖点的应力强度因子如表1所示。’。表1应力强度因子由表1可知，ansys计算结果与应力强度因子手册公式计算结果相差均在5％以内，分析结果