01_断裂参数的数值计算方法_02.pptVIP

Fracture Mechanics 华中科技大学船海学院 袁锐 Fracture Mechanics 华中科技大学船海学院 袁锐 * HUST-XD-B-0111-2014 断裂参数的数值计算方法 * 断裂参数的数值计算方法 1.6.1 全局虚拟裂纹扩展法 研究裂纹扩展的规律，一种方法是分析裂纹尖端的应力应变场，从而得到表征裂纹尖端应力应变场强度的参量—应力强度因子K；另一种方法是用能量平衡的观点，考察裂纹扩展过程中物体能量的转化，从而得到表征裂纹扩展时能量变化的参数—能量释放率G。 考察如图所示的I型裂纹问题，板厚为B，裂纹长度为a。假设外载荷P慢慢增加，裂纹也随之沿延长线扩展。裂纹在失稳开裂前，扩展了面积dA=Bda。令2γ为裂纹扩展单位面积所需要的表面能，则产生dA面积需要外界提供的能量为dГ=2γdA。这时，外载荷P对裂纹板所作的工dW，一部分变成弹性应变能dV，另一部分由于形成裂纹新表面所消耗，转化成表面能dГ。根据能量守恒定律，有： 或在单位时间内，有 能量释放率是根据Griffith的能量释放观点而定义的量，其定义如下:能量释放率是指裂纹由某一端点向前扩展一个单位长度时(当裂纹不扩展时，外力功等于增加的应变能)，平板每单位厚度所释放出来的能量。为了纪念Griffith的功绩，用G来表示，其单位为MN/m（百万牛顿每米）。 * 断裂参数的数值计算方法 1.6.1 全局虚拟裂纹扩展法 * W，V是外载荷P和裂纹面积A=Ba的函数，表面能Г只是A的函数（对于均质材料而言），而P，A又是时间t的函数。如果考虑静态加载，即dP/dt=0，则式（1-48）可写成： 系统的位能II=V-W，因此上式变为： 上式的左边是裂纹扩展单位面积时，整个受力系统释放的能量（弹性位能），称为能量释放率，用G表示（I型裂纹用GI）： 它是与外载荷及结构形式（裂纹长度、形状、位置等）有关的一个力学参数，量纲为N.m-1.从单位来看，它可以被看做裂纹扩展单位长度时所需要的力，可看做企图驱动裂纹扩展的原动力，故又称作裂纹扩展力。式1.50的右边是裂纹扩展单位面积所需要的表面能（又称表面张力），是与材料有关的参数，可看做是材料的常数，故称之为临界能量释放率GIc。 进一步思考可能扩展时的情况！ 1.6.1 全局虚拟裂纹扩展法 断裂参数的数值计算方法 * 断裂参数的数值计算方法 在线弹性断裂力学中，从能量的观点考察裂纹的扩展过程是研究破坏机理的一条有效途径。如果裂纹扩展过程中系统所释放的能量，足以提供裂纹扩展所需要的能量，则裂纹发生扩展。 图 1 图 2(a) 图 2(b) * 方法一 对于图1所示的问题，我们可以把它看成是2(a)和图2(b)两种情况的叠加。即图1可以看成在图2(a)的基础上，逐渐松弛裂纹处的应力（即叠加图2(b)）而形成的。这是一个无裂纹到裂纹扩展为长度2a的过程。 图2(a) 无裂纹的体在无穷远处受均匀拉力，此状态下物体的应变能为U0。 图2(b) 大小形状完全相同的物体在裂纹表面受有均匀的后应力。 W分为两部分，一部分是松弛过程中，反向应力在裂纹表面所做的功，恒为负，另一部分就是在此过程中，无穷远处作用的应力由于位移的变化而做的功，恒为零，因为在弹性体内作用一对平衡力对无穷远处的位移没有影响。 图1 状态下物体的应变能为U。 * 方法二 应力强度因子准则和能量准则是从两种不同观点建立起来的准则，它们之间有一定的关系。假设裂纹沿其延长线扩展时，计算裂纹尖端应力松弛的功，可求得能量释放率和应力强度因子的关系。 设I 型裂纹尖端处于平面应变状态，在固定施力点的情况下裂纹沿其延长线扩展，裂纹尺寸从a增加到a+da，如图所示。此时，裂纹扩展段da上的应力消失，新表面在法线方向产生位移。比较裂纹长度为a+da和裂纹长度为a时板的能量差异就可以求出能量释放率。先把裂纹长度看做a+da，计算张开位移（x的原点在裂纹尖端）： * 正负号分别对应于上下表面。下角标a+da表示裂纹长度为a+da时的应力强度因子。在da段裂纹表面上加载，载荷的分布和大小按照表面的应力大小取值： 加载值由0增加到上式的最终值，da段裂纹逐渐靠拢知道完全闭合，恢复到裂纹为a的状态。因此从a+da状态恢复到a状态的加载过程中，外力对弹性体做了功，提高了应变能。由于任意x处的 和 之间存在线性关系，做功为： 式中，B为物体的厚度。反之，从a状态变为a+da状态，应力从 松弛到零，位移从 零变化到 ，松弛过程中释放的应变能dV=-dW。所以： 方法二 * 理解两种方法的共同性 * 1.6.1 全局虚拟裂纹扩展法 当GI=GIc时，裂纹可