断裂与损伤力学大作业二.docVIP

《断裂与损伤力学》大作业二 含孔板和裂纹板的应力强度因子分析 姓 名 学 号 学 院 航空宇航学院 专 业 工程力学 邮 箱 任课教师 郭树祥 2013年6月 含圆孔和裂纹板的应力强度因子分析 模型与要求 1.1 模型 含两圆孔和一倾斜裂纹的矩形板受到垂直方向拉伸载荷的作用，如图1所示，矩形板的水平、垂直距离分别为2W、2H，孔的直径为D，裂纹长度为2a，裂纹与水平方向的夹角为θ，两圆孔的水平、垂直距离分别为2L、2h。 图1 含两圆孔和倾斜裂纹矩形板受垂直拉伸载荷作用 1.2 要求 计算中心裂纹尖端的应力强度因子KⅠ和KⅡL、h、a、D、θ等的变化规律； 绘制应力强度因子随几何参数的变化曲线； 列出必要的参考文献。 方案 理论分析 2.1.1 线弹性断裂力学 Ⅰ型裂纹尖端的应力场为： （1） Ⅰ型裂纹尖端的位移场为： （2） 其中， （3） 从上面可以看出，对于Ⅰ型裂纹而言，裂纹尖端的应力场和位移场均可以表示成应力强度因子的形式，所以可以通过裂纹尖端的应力应变场求其裂尖应力强度因子。这也正是传统有限元求解应力强度因子的原理。但是从上面的表达式同样可以看出，在裂纹尖端应力具有1/2的奇异性，其趋向于无穷大；而位移则趋于零，所以在裂纹尖端应力场和位移场均具有很大的梯度，所以就需要划分很精细的网格来求解应力位移场。 有限元在计算裂纹尖端的应力应变场时，通常在裂纹尖端通过引入奇异单元来模拟应力的奇异性，这样即使在单元数目有限的情况仍然能很好地求解出裂纹尖端的应力应变场。 2.1.2 J积分求解裂纹尖端的应力强度因子 传统的有限元在计算裂纹尖端的应力强度因子的时候，无可避免地遇到裂尖复杂应力场和位移场的计算，J积分则可以完全避免这种复杂的处理过程。 1968年由Rice和Cherepanov提出了一个围绕裂纹尖端的围线积分，该积分与路径无关，保持为一个常数，并且可以反映出裂纹尖端的应力应变场。 （4） 图2 J积分的定义 坐标变换： 得到 则 在线弹性情况下： （5） 将（1）式带入到式（5）中，并令r=R，得到： （6） （7） （8） 将式（2），（6），（7），（8）带入到（4）式中可以得到： （9） 从（9）式中可以得出，Ⅰ型裂纹裂尖应力强度因子KⅠ与J积分具有一一对应的关系，所以可以通过求解J积分得到裂尖的应力强度因子。 为了简便，上面的推导过程均是针对纯Ⅰ型裂纹而言。对于Ⅰ型和Ⅱ型复合型裂纹，推导过程类似。 在Ⅰ型和Ⅱ型复合裂纹情况下，求解出来的J积分是KⅠ和KⅡ的函数，可以通过J积分求解出Ⅰ型和Ⅱ型裂纹的应力强度因子KⅠ和KⅡ。但是只求解出一个J积分是无法得到两个应力强度因子的，Ishikawa将Ⅰ型和Ⅱ型复合型裂纹的J积分进行了如下分解（图3）： （i＝1，2） 在计算出和后，就可以求解出Ⅰ型和Ⅱ型复合裂纹的应力强度因子和。 图3 复合型裂纹J积分的分解 1、假定材料尺寸为120×100，边界条件为在左边下一点限定x、y两方向位移，在底边限定y方向位移，这样可以体现对称变形，符合实际情况。 裂纹区域设置如图4。 图4 裂纹区域设置图 载荷边界图5 图5 载荷边界图 4、板的有限元网格划分在裂纹的尖端网格需细分，而在其他的区域网格的划分可以稀疏一些，且以四节点等参单元为主，裂纹附近网格的划分见下图6。 图6 网格划分图 5、图7是变形后板的应力场分布。 图7 板应力分布图 6、裂尖处应力分布图8 图8 裂尖应力分布图 3 结果与分析 实际计算时，选用的材料板，宽2W=120mm，高2H=100mm，厚度为b=5mm，E