03_线状裂纹的虚拟裂纹闭合法.pptVIP

Fracture Mechanics 华中科技大学船海学院 袁锐 Fracture Mechanics 华中科技大学船海学院 袁锐 * HUST-XD-B-0111-2014 线状裂纹的虚拟裂纹闭合法 * * 3.1 虚拟裂纹闭合法的数学解释 目前有两种方式对虚拟裂纹闭合法进行数学解释： 1 Raju方法 基本假设：虚拟裂纹尖端后面的张开位移和实际裂纹尖端后面的张开位移近似相等。 2 Xie-Waas方法 基本假设：虚拟子裂纹尖端前沿的应力分布和实际主裂纹尖端前沿的应力分布规律近似相同。 * 3.1 虚拟裂纹闭合法的数学解释 3.1.1 Raju 方法 图3.1显示了一个长度为a的主裂纹，虚拟扩展了一个增量为?a的微小裂纹。Irwin的裂纹闭合积分可以写为： 图3.1 虚拟裂纹闭合法示意图 * 将裂纹从a扩展到a+?a所需要的功与将裂纹从a+?a闭合到a所需的功是相等的。如图3.2所示，考虑裂纹尖端附近的有限元网格，虚拟裂纹线上分布的应力做的功等于节点力在节点位移上做的功。即： 3.1.1 Raju 方法 3.1 虚拟裂纹闭合法的数学解释 分布于裂纹尖端前沿且趋近于裂纹尖端区域的应力可以表示为： * 对于线性四边形单元，位移可以插值为： 3.1.1 Raju 方法 将公式3.4和3.3代入公式3.2中，并消去 得： 依据基本假设：虚拟裂纹尖端后面的张开位移和初始裂纹尖端后面的张开位移近似相等，公式3.1可重写为： 公式3.6的优点是它可以在一步有限元分析中计算出来。如图3.1所示，裂纹尖端前面的应力分布为： 基于内插值，裂纹尖端后面的张开位移为： 将公式3.7和公式3.8代入公式3.6中，得到： 3.1 虚拟裂纹闭合法的数学解释 * 3.1.1 Raju 方法 3.1 虚拟裂纹闭合法的数学解释 公式3.9括号中部分简化为： 公式3.10中的极限表明裂纹尖端的单元需要很小，并且裂纹尖端处的单元越小，所得的GI就越精确，通过数值实验可以确定必要的网格细化程度。然而，大量的数值实验表明公式3.10的值对有限元网格尺寸并不敏感。因此，公式3.10中的极限可以很好的近似的表达为： 类似的，可以得到II型的计算公式如下： * 3.1.2 Xie-Waas 方法 3.1 虚拟裂纹闭合法的数学解释 如图3.3所示，一个主裂纹（parent crack），其应力强度因子为KI。主裂纹前端有一个虚拟的微小子裂纹（child crack），其应力强度因子为 。根据线弹性断裂力学： 上述公式对于主裂纹和子裂纹都是成立的。 图3.3 主裂纹和虚拟示意图子裂纹 * 3.1.2 Xie-Waas 方法 3.1 虚拟裂纹闭合法的数学解释 通过虚拟裂纹扩展法（参看公式1.26a），并参考图3.2，能量释放率可用下式计算： 第一步计算的节点力 第二步计算的节点张开位移 对于主裂纹而言，上平面（ ）和下平面（ ）的垂直位移分别为： 主裂纹后部的张开位移很容易得到，为： 同样地，对虚拟子裂纹而言，有 1 5 IV III I II 3 4 Δa Δc X Y 1 5 3 4 X Y 1’ IV III I II 注意上标不同！ * 3.1 虚拟裂纹闭合法的数学解释 3.1.2 Xie-Waas 方法 接下来建立 和 之间的关系： 对于主裂纹而言，沿x轴方向，裂纹前端应力分布为： 依据基本假设：虚拟子裂纹尖端前沿的应力分布和初始裂纹尖端前沿的应力分布规律近似相同。式3.16所描述的应力分布同样适用于裂纹。因此，根据应力强度因子的定义，有 注意上标不同！ 当裂纹尖端使用相同的网格尺寸时，即 ，公式3.19退化为： 由于F和v均只是在第一步中计算，因此可以省去上标，这也是公式3.20的优点！ * 3.1 虚拟裂纹闭合法的数学解释 3.1.3 节点力的计算 使用虚拟裂纹闭合法时涉及节点力和节点位移。由于节点位移是有限元分析中的基本变量，因而可以直接获得，相对简单，但在计算节点力时则应该注意。注意到计算中使用的节点力是内力，内力将裂纹尖端的单元共享节点约束在一起。因此，裂纹尖端周围的单元对这些节点力都有贡献。图3.4显示了在裂纹尖端周围有四个单元的情况下的自由体受力图。如果没有外力作用在裂纹面上，节点力在每个坐标方向上的合力为0，即： 上标代表图中的单元编号。因此，节点力既可以通过上面的两个单元（单元1和单元2）计算，也可以通过下面的两个单元计算： I IV II III