圆孔裂纹板K值Ansys分析报告.docxVIP

含圆孔和裂纹板的应力强度因子分析一 问题描述含两圆孔和一倾斜裂纹的矩形板受垂直方向拉伸载荷的作用，如图1所示，孔的直径为D，裂纹长度为2a。试计算中心裂纹尖端的应力强度因子KＩ和KⅡ,并讨论其随几何参数L、h、a、D、θ等的变化规律，写一份分析报告。要求：报告中计算所用的分析方法和模型应阐述清楚，并写出必要的计算公式等。绘制应力强度因子随几何参数的变化曲线。列出必要的参考文献。图1 含两圆孔和倾斜裂纹矩形板受垂直拉伸载荷作用二 理论分析在线弹性断裂力学中，I型裂纹尖端的应力场为：I型裂纹尖端的位移场为：其中：同理，对II型裂纹尖端的应力场：显然，位移场和应力场均可以表示成应力强度因子的形式。通过对裂纹尖端的应力应变场分析来求解对应的应力强度因子，便是传统有限元求解应力强度因子的原理。而对于I、II复合型裂纹尖端的应力强度因子，可通过它们的叠加获得。计算应力强度因子常用方法有3大类：解析法、数值计算法和工程计算法。解析法分为复变函数法和积分变化法。解析法只能计算简单问题，大多数问题需要采用数值计算法。数值计算法可分为有限元法、加权残差法、边界元法和权函数法。当前工程中广泛采用的数值解法是有限单元法。随着有限元法的发展，有限元在断裂力学中的应用越来越普及。近年来，计算机技术得到了迅猛发展，许多功能强大的有限元软件也相继问世，大型通用有限元程序就是当前工程中应用最广泛的有限元软件之一。分析该问题采用有限元软件Ansys计算裂纹尖端的应力强度因子，通过阅读Ansys的帮助文件，得到Ansys基于有限元方法在线弹性范围内计算应力强度因子的原理。Ansys得出的应力强度因子使用的是位移外推法。三 模型建立图2 模型网格利用对称性，建立矩形板1/2的模型。定义模型参数，令W=10mm，H=8mm，a=2mm，θ=30°，L=5mm，h=4mm，D=2mm。为了模拟裂纹尖端奇异性，在裂纹尖端常使用1/4节点奇异单元。1 前处理新建plane183，选用6自由度三角形单元，板厚t=1mm。确定材料参数，弹性形模量E=2×105MPa，泊松比μ=0.3。建立模型，划分网格。裂尖和圆孔周围比较密集，其余部分较为稀疏，如图2.2 加载并求解由于模型建立的是矩形板的一半，因此模型左侧没有横向位移。在模型左侧施加位移约束。在矩形板上下两端施加拉应力σ=10MPa。载荷施加完成后进行求解（solution） 3 后处理使用通用后处理查看应力云图，如图3。由图3可以看出，裂纹尖端应力最大，圆孔两侧的应力值也较大。在裂纹尖端建立局部坐标，X轴与裂纹面平行。定义裂纹面路径，裂纹尖端以及上下裂纹面各取两点。General postproc Nodal calcs Stress Int Factr得出应力强度因子。图3 应力云图四 结果与讨论由于需要讨论应力强度因子与几何参数L、h、a、D、θ等的变化规律，则采用控制变量法建立对比模型。通过改变关键点的坐标，从而实现参数的改变。1 参数a对应力强度因子KⅠ、KⅡ的影响令θ=30°,L=5mm，h=4mm，D=2mm，a为变量，建立模型，记录应力强度因子KⅠ、KⅡ值,如表1。表1 a为变量应力强度因子值a / mm22.53KⅠ / MPamm1/225.16429.55134.492KⅡ / MPamm1/21.83652.13992.4120使用origin绘制应力强度因子随a的变化曲线，如图4。从曲线图中可以看出，应力强度因子KⅠ、KⅡ随a值增加而增大。说明应力强度因子与裂纹长度呈正相关，裂纹越长则裂纹尖端应力强度因子越大。图4 应力强度因子K随a变化曲线2参数θ对应力强度因子KⅠ、KⅡ的影响令a=2mm，L=5mm，h=4mm，D=2mm，θ为变量，建立模型，记录应力强度因子KⅠ、KⅡ值，如表2。表2 θ为变量应力强度因子值θ / °304560KⅠ / MPamm1/225.16419.94513.899KⅡ / MPamm1/21.83652.36822.2482使用origin绘制应力强度因子随θ的变化曲线，如图5。图5 应力强度因子K随θ变化曲线图中可以看出，KⅠ随θ值增大而减小，KⅡ随θ值变化而先增大后减小。当θ=0°时，则完全为Ⅰ型裂纹，KⅠ取得最大值，之后逐渐减小。KⅡ变化规律是先增后减，在45°时达到最大值。3参数L对应力强度因子KⅠ、KⅡ的影响令a=2mm，θ=30°,h=4mm，D=2mm，L为变量，建立模型，记录应力强度因子KⅠ、KⅡ值，如表3。表3 L为变量应力强度因子值L / mm55.56KⅠ / MPamm1/225.16424.63724.165KⅡ / MPamm1/21.83651.84131.7969使用origin绘制应力强度因子随L的变化曲线，如图6。图6 应力强度