断裂力学——3裂纹尖端应力场和位移场计算

1、1Shanghai University断裂力学Fracture Mechanics郭战胜郭战胜办公地点：延长校区力学所办公地点：延长校区力学所317室室平时答疑：每周一：平时答疑：每周一：5-6节节晚修答疑晚修答疑:每周一：每周一：18：00-20：30地点：地点：HE108或或HE104b2裂纹尖端附近的应力场和位移计算裂纹尖端附近的应力场和位移计算343cos(1 sinsin)2222xKr3cos(1 sinsin)2222yKr3cossincos2222xyKr0 xzyzReImxZyZReImyZyZ0z(平面应力) ()2RezxyZ (平面应变) RexyyZ 2Iiji

2、jKfr用张量标记可缩写成型裂纹求解53(21)coscos4222KrukG3(21)sinsin4222KrvkG0w平面应变 ()xywdzE 平面应力 3431k平面应变平面应力1(1)Re(1) ImuZyZE12Im(1) RevZyZE平面应力1(12 )ReImuZyZE12(1)ImRevZyZE平面应变型裂纹求解60ra IZ需要注意的是，推导过程中，使用了这个条件，所以。对于稍远处，应该用 所示的来确定应力分量和位移分量。前面得到的应力场和位移场公式只适用于裂纹尖端附近区域，即要求( )( )(2 )afZa （ + ）型裂纹求解7型裂纹求解 设无限大板含长2a的中心裂纹

3、，无穷远受剪应力作用8第一步：解第一步：解IIII型型WestergaardWestergaard应力函数应力函数求解方法与求解方法与I I型基本相同，主要差别是无穷远处边界上受力条件不型基本相同，主要差别是无穷远处边界上受力条件不同。选取应力函数同。选取应力函数 所以所以 ReIIIIyZzx ReImIIIIIIZzyZzy ReReZ zZzx ReImZ zZzy ImReZ zZzy因为因为 22ReIIIIyZzx 222ImReIIIIIIZzyZzy 2ReImIIIIIIZzyZzxy ReZIIy型裂纹求解 9得到得到IIII型裂纹问题各应力分量表达式为型裂纹问题各应力分量

4、表达式为 ReIm2ZyZxReZyyZyZxyImRe进而可得到位移分量进而可得到位移分量ZyZEvZyZEuImRe)21 ()1 (ReIm)1 (2)1 (平面应变平面应变 型裂纹求解 10第二步：选第二步：选IIII型裂纹的型裂纹的 ( )Zz边界条件：边界条件： 0 xyy0yax 0 zyxyz， 在 处在处选取选取 22( )zZzza能够满足全部边界条件。能够满足全部边界条件。型裂纹求解 112223/222lim( )limlim( )lim0zzzzzZzzaaZzza在裂纹表面在裂纹表面 处处 0yax 2222( )zxZzzaxa虚数虚数 Re( )0Zz0 xyy

5、ReIm2ZyZxReZyyZyZxyImRez 只有实部且为一常数 0IIZz0 xyxy满足平板周围的边界条件 满足裂纹表面处的边界条件 型裂纹求解 12将坐标原点移到右裂尖，采用新坐标将坐标原点移到右裂尖，采用新坐标 az ( )( )2afZa 当当 0)(f趋于常数趋于常数, ,设设: : ， 00lim( )lim( )2KfZ右裂尖附近右裂尖附近, , 在很小范围内时在很小范围内时 0lim2( )KZ用解析函数求解II型裂纹尖端应力强度因子的定义式 型裂纹求解 13第三步：用第三步：用 求求IIII型裂尖附近的应力场和位移场型裂尖附近的应力场和位移场 ( )Zz 应力强度因子是

6、在裂尖时应力强度因子是在裂尖时 存在极限，若考虑裂尖附近存在极限，若考虑裂尖附近的一个微小区域，则有：的一个微小区域，则有：02( )KZ( )2KZ若以极坐标表示复变量若以极坐标表示复变量 )sin(cosirrei( )(cossin)222KZir则可得到则可得到 332233cossin22222 2IIIIIIIIKKKZrir sin2 sincos22yrr型裂纹求解 143sin(2coscos)2222xKr 3cossincos2222yKr 3cos(1 sinsin)2222xyKr0 xzyz()zxy 平面应变 0z平面应力 3(23)sinsin4222KrukG

7、3(22)coscos4222KrvkG 把上面两式代入前面应力表达式中，应力和位移场得表达式把上面两式代入前面应力表达式中，应力和位移场得表达式3134k平面应力平面应变型裂纹求解 15对于I型和II型裂纹来说，是属于平面问题。但对于III型裂纹，由于裂纹面是沿z方向错开，因此平行于xy平面的位移为零，只有z方向的位移不等于零 型裂纹求解对于此类反平面问题，前面给出的平面问题的基本方程已不适用，因此不能沿用Airy应力函数求解，需要从弹性力学的一般（空间）问题出发，推导公式。弹性力学一般问题的基本方程，可以仿照平面问题的方法导出 16反平面(纵向剪切)问题, 其位移 ( , ),0ww x

8、y uv根据几何方程和物理方程:1xzxzwrxG1yzyzwryG0 xyxyz型裂纹求解问题描述:无限大板,中心裂纹(穿透) ,无限远处受与 方向平行的 作用.2az17单元体的平衡方程:0yzxzxy位移函数满足Laplace方程,所以为调和函数. 解析函数性质:任意解析函数的实部和虚部都是解析的.1( , )Im( )w x yZzGImImxzZwGZxxImReyzZwGZyy边界条件边界条件: :0,0yzyxa,0,xzyzz型裂纹求解222220wwwxy非零应力分量18选取函数 22( )lzZzza满足边界条件 型裂纹求解在裂纹表面在裂纹表面 处，处， 0yax IIIZ

9、z只有实部而无虚部，有 0yz满足裂纹表面处满足裂纹表面处的边界条件的边界条件 y x IITlZz ReIIIlZz Im0IIIZz ,0yzlxz 当或，都有，即由非零应力分量公式知，满足平板周围的边界条件。满足平板周围的边界条件。 19取新坐标 za ()1( )(2 )IIIaZfa 型裂纹求解同样，为计算方便，将坐标原点从裂纹的中心移到裂纹的右尖端 当当 0)(f趋于常数趋于常数, ,设设: : ， II00lim( )lim( )2KfZ右裂尖附近右裂尖附近, , 在很小范围内时在很小范围内时 II0lim2( )KZ 用解析函数求解III型裂纹尖端应力强度因子的定义式 20 应

10、力强度因子是在裂尖时应力强度因子是在裂尖时 存在极限，若考虑裂尖附近存在极限，若考虑裂尖附近的一个微小区域，则有：的一个微小区域，则有：0II2( )KZ II( )2KZ若以极坐标表示复变量若以极坐标表示复变量 )sin(cosirrei则可得到则可得到 IIIIIIcossin222KZir Recos22Imsin22IIIIIIIIIIIIKZrKZr sin22cos22IIIxzIIIyzKrKr 这就是III型裂纹问题在裂纹尖端附近的应力场表达式 型裂纹求解21则可得到则可得到 IIIIIIcossin222KZir这就是III型裂纹问题在裂纹尖端附近的位移场表达式 1222co

11、ssin2222IIIIIIIIIIIIKKrZdKi 2Recos22Imsin2IIIIIIIIIIIIrZzKrZzK2sin2IIIKrwG型裂纹求解22( )(2 )aZa （ + ） 0 lim2IIKZa应力强度因子( )2aZa 0lim2( )KZa ()( )(2 )laZa II0lim2( )lKZa 注意：以上三种类型求解方法，仅适用于含贯穿裂纹的无限大板在注意：以上三种类型求解方法，仅适用于含贯穿裂纹的无限大板在载荷或位移对裂纹中点的坐标轴对称或反对称的情况。载荷或位移对裂纹中点的坐标轴对称或反对称的情况。23值得指出的是，上述三种裂纹问题的应力场表达式，虽然是根据

12、无限大半具有中心穿透裂纹且在均匀外加应力作用下获得的。进一步的分析表明，这些解具有普遍的意义，也就是说，对于其他有限尺寸板的穿透裂纹（包括中心裂纹和边裂纹），在非均匀受力条件下，裂纹尖端附近的应力场（更确切地说是应力场的奇异项）表达式也是相同的，其不同之处仅仅是应力强度因子的不同，因此，对于特定的含裂纹结构只需要确定相应的应力强度因子就可以了。24通过前面的推导，各种类型裂尖应力和位移场可表示为通过前面的推导，各种类型裂尖应力和位移场可表示为)(2) I () I (ijijfrK3 , 2 , 1,ji)() I () I (iigrKu3 , 2 , 1i若上标写成若上标写成IIII、II

13、IIII，代表，代表IIII型或型或IIIIII型裂纹。型裂纹。裂纹尖端应力场是渐进解，仅仅适合于裂纹尖端附近裂纹尖端应力场是渐进解，仅仅适合于裂纹尖端附近25线弹性裂尖场特点线弹性裂尖场特点v三种变形情况下裂纹尖端应力场和应变场都具有奇异奇异性性，即在裂纹尖端处，应力和应变为无穷大，这种不真实的性质是由于所采用的本构关系所决定的，即认为材料能承受无限大的应力，且应变与应力呈线性关系。另外，在上述的分析中，裂纹假设成理想的尖裂纹，即裂纹尖端曲率为无穷大。实际上，裂纹尖端不可避免地会出现塑性区，并且裂纹尖端地曲率是有限的，但是在塑性区很小的情况下，在围绕裂尖的一个环状区域环状区域内K场场是适用的

14、。vK场内的位移与 成线性比例关系。12r26线弹性裂尖场特点线弹性裂尖场特点v三种情况下的K场有相似的形式，分别由应力强度因子决定着其场的强度。SIF取决于外加载荷，而且与构件几何、裂纹尺寸有关，但是与( )坐标无关。在K场范围内，应力和应变均正比于SIF，所以SIF是裂纹尖端附近应力、应变场强度的表征，是描述裂尖场强度的参数。v裂尖场与角分布函数成比例。角分布函数仅与角 有关，而与r无关，对于同一种变形模式，角分布函数是相同的。所以，无论构件的形状、尺寸以及裂纹的尺寸如何，裂尖场都是相同的。 r27o 应力不适宜作为判断含裂纹材料承载能力的力学参量应力不适宜作为判断含裂纹材料承载能力的力学

15、参量裂裂尖场应力具有奇异性，只要存在载荷，应力就趋于无穷大。尖场应力具有奇异性，只要存在载荷，应力就趋于无穷大。依照传统强度理论，含裂纹结构必定破坏。即依照传统强度理论，含裂纹结构必定破坏。即传统的强度传统的强度条件判断准则失去意义。条件判断准则失去意义。o 应力强度因子作为判定裂纹尖端应力场强度的物理参量引入。应力强度因子作为判定裂纹尖端应力场强度的物理参量引入。 线弹性断裂力学的主要任务之一就是确定含裂纹构件的线弹性断裂力学的主要任务之一就是确定含裂纹构件的应力强度因子。应力强度因子。应力强度因子是有限量，它是代表应应力强度因子是有限量，它是代表应力场强度的物理量，用其作为参量建立破坏条件

16、是力场强度的物理量，用其作为参量建立破坏条件是合适的。合适的。 应力强度因子应力强度因子28aYK名义应力，即裂纹位置上按无裂纹计算的应力名义应力，即裂纹位置上按无裂纹计算的应力aY裂纹尺寸，即裂纹长或深裂纹尺寸，即裂纹长或深 形状系数，与裂纹大小、位置有关形状系数，与裂纹大小、位置有关应力强度因子一般写为：应力强度因子一般写为：应力强度因子单位：应力强度因子单位：N.mN.m-3/2-3/2应力强度因子应力强度因子29应力强度因子 鉴于应力强度因子应力强度因子的重要性，在断裂力学这门科学近半个世纪的快速发展中，应力强度因子的分析计算一直是一个经久不衰的研究课题，这可从这方面的专著（如二十世纪

17、七十年代Sih的专著和近期的专著）和专门的应力强度因子手册可见一斑。从研究方法上，从解析的Westergaard stress function、 Muskhelishvili stress function 到解析的或半解析的Green Function、Singular Integral Equation、Conforming Mapping（保形映射）, 及数值方法如Boundary Collocation Method， Finite Element Method (有限元法)和Boundary Element Method (边界元法)。 30脆性断裂的脆性断裂的K K准则准则应力强

18、度因子与应变能释放率的关系应力强度因子与应变能释放率的关系 根据前面所述的应变能释放率公式 与应力强度因子 可以发现它们之间应有一定关系。这关系将进一步揭示应力强度因子的物理意义。EaG21aK 以张开型裂纹为例，由于应变能释放率代表裂纹扩展以张开型裂纹为例，由于应变能释放率代表裂纹扩展单位面积所释放的应变能。那么逆向思维一下单位面积所释放的应变能。那么逆向思维一下31左图a所示裂纹原长为a，扩展微小长度 （图b）后，释放出的能量可用从图b状态闭合到图c状态所作的功来计算。闭合时作用在裂纹上表面上x位置的应力由图b中的0值，逐渐增加到图a中的 a)(xy利用上节的裂尖附近应力和位移场，可以计算

19、使裂纹闭合单位面积所利用上节的裂尖附近应力和位移场，可以计算使裂纹闭合单位面积所作的功，显然这部分功应该等于裂纹扩展单位面积所释放的能量。作的功，显然这部分功应该等于裂纹扩展单位面积所释放的能量。32由I型裂纹的应力表达式， 当 ， 时 xr 0 xKxy2)(由图b看出，闭合时的位移最初为 其中 ， ),(rvxar注意：图注意：图b b与图与图a a的坐标原点不同。的坐标原点不同。由I型裂纹的位移表达式： )22(24),(kxaKxav闭合后，位移为0。 闭合过程中，应力在 段所作的功为 avBdxya03(21)sinsin4222Krvk3cos1 sinsin2222yKr33闭合

20、单位面积所作的功裂纹扩展单位面积所释放的能量=由于：2200114aayvBdxKkaxGKdxB aaxE 其中， （平面应力）， （平面应变） EE 21EE可见，应力强度因子与应变能释放率有对应关系： 不仅表示裂尖附近弹性应力场的强度，也可确定裂纹扩展释放的能量率，故：对于线弹性断裂问题， 与 等价KG0/ 2aaxdxax 34 同理，对于II型和III型裂纹同样可得到类似关系2EKG2(1)KGE 需要注意：对于需要注意：对于I I型和型和IIII型裂纹问题可分为平面应力和平面应变问题，型裂纹问题可分为平面应力和平面应变问题，而对于三型裂纹问题只是一种反平面问题。而对于三型裂纹问题只

21、是一种反平面问题。脆性断裂的脆性断裂的K K准则准则我们已经讲了脆性材料裂纹失稳扩展的临界条件为：CGG1135可以得到以应力强度因子表示的裂纹失稳扩展的临界条件为：可以得到以应力强度因子表示的裂纹失稳扩展的临界条件为：IICKK表示裂尖的应力强度因子表示裂尖的应力强度因子 达到达到 时，裂纹失稳扩展。时，裂纹失稳扩展。 与与 都是材料常数，称为材料的都是材料常数，称为材料的平面应变断裂韧度平面应变断裂韧度。在线弹性条件下在线弹性条件下1KCK1CG12EKGCC强调：强调： 与与 概念不同，概念不同， 是表示裂尖应力场强度的一个参量，可用弹性理论是表示裂尖应力场强度的一个参量，可用弹性理论方

22、法进行计算，由载荷及裂纹体形状和尺寸决定，方法进行计算，由载荷及裂纹体形状和尺寸决定， 断裂韧度，材料具有的一种机械性能，表示材料抵断裂韧度，材料具有的一种机械性能，表示材料抵抗脆性断裂的能力，由试验测定。抗脆性断裂的能力，由试验测定。IKICKICK脆性断裂的脆性断裂的K K准则准则CK1IK36注意：对于线弹性断裂问题，采用G准则和K准则所得的结果是一样的。但是由于利用弹性理论可直接计算应力强度因子，而且试验测定 比 测定方便，故工程一般常用K准则。CK1CG1 根据K准则，可以计算剩余强度（临界应力）和临界裂纹长度，进行断裂安全分析。 例如：对具中心裂纹无限大板，受双轴拉应力 CKaK1

23、11/CCKa)/(221CKa 对于其它结构， 表达式不同。1K可得 37 根据实验和理论分析，断裂韧度随试件厚度增加而下降，根据实验和理论分析，断裂韧度随试件厚度增加而下降，如下图。这是由于：如下图。这是由于：1 1）薄板的裂尖处于平面应力状态，断裂韧度较高，裂纹不易）薄板的裂尖处于平面应力状态，断裂韧度较高，裂纹不易扩展，用扩展，用 表示；表示；2 2）随板厚增加，裂尖处于平面应变状态的部分增加，裂纹较）随板厚增加，裂尖处于平面应变状态的部分增加，裂纹较易扩展，断裂韧度降低，当厚度降至一定值后，断裂韧度降易扩展，断裂韧度降低，当厚度降至一定值后，断裂韧度降至最小，称为平面应变断裂韧度，用

24、至最小，称为平面应变断裂韧度，用 表示。表示。CKICK断裂韧度与板厚的关系断裂韧度与板厚的关系需要注意：金属在平面应力条件下裂尖产生较大塑性变形，需要注意：金属在平面应力条件下裂尖产生较大塑性变形，K K准则（建立在线弹准则（建立在线弹性断裂力学基础上）不适用，而要采用第三章的弹塑性断裂力学的断裂准则。性断裂力学基础上）不适用，而要采用第三章的弹塑性断裂力学的断裂准则。但是当裂尖塑性变形区较小时，通过下一节的修正后，仍可用但是当裂尖塑性变形区较小时，通过下一节的修正后，仍可用K K准则。准则。38线弹性断裂力学在小范围屈服时的推广线弹性断裂力学在小范围屈服时的推广39屈服条件s1scf),(

25、321cf 单向拉压：薄壁圆筒扭转：在应力空间 在主应力主应力空间谓之屈服条件或屈服面方程单向应力复杂应力* 0, ffc或谓之屈服函数c塑性约束系数yssc有效屈服应力，材料屈服点(,)xyzxyxzyzfc 40特雷斯卡（Tresca）假设最大剪应力最大剪应力是屈服的控制因素 122331 , ccc即,0)()()(221322322221*cccf材料屈服，屈服函数为：在主应力空间是六棱柱，在12平面是六边形 时，41在 平面是六角形 C12c 12c1c1c 2c2c 122331 ,ccc即,12CC-C-C1242 米泽斯 （Mises）假设)()()(121213232221G

26、U形控制因素是形状改变比能形状改变比能（歪形能、畸变能） Mises屈服条件为：UC形即：2221223311()()() 12CG222122331()()()c或：即Mises屈服条件或屈服方程。4312 0, 3,在平面有2221212()c在主应力空间，屈服面是圆柱222122331()()()c221122:c 化简 是椭圆方程（屈服曲线） 3 2 1n 平 面Mises屈服条件44在主应力空间屈服面是圆柱221122:c 椭圆方程 C可由简单实验求出 与六棱柱外接CC-C-C1245 C可由简单实验求出 如：由Mises屈服条件：222122331()()()C123,0s2222

27、sssC单向拉伸屈服时，22132322212)()()(ss)()()(21213232221即：或：46纯剪切屈服如由Mises屈服条件：222122331()()()C132 0 ss2226sssC纯剪切屈服时，2221223312()()()6s即：22, 62 3ssssCConstC47 线弹性裂纹尖端场，其应力场具有 r -1/2 的奇异性， 该奇异性的幅值大小可用应力强度因子应力强度因子来表征。 但从物理学物理学的角度来看，真正奇异的应力是不存在的，也就是说在裂纹尖端附近很小的范围说在裂纹尖端附近很小的范围内，内，K K场是不适用的场是不适用的。在裂纹尖端附近的材料必定发生屈

28、服。在外加载荷作用下裂纹尖端的应力有限有两种原因。其一是裂尖附近由于应力集中,裂尖的材料会发生不同程度的塑性变形，其二是裂纹尖端并不是理想的曲率无穷大的形状，而总是有钝化的。 Irwin 小范围屈服理论48那么线弹性断裂力学能否继续使用呢？如果裂纹尖端的塑性区（或者说偏离K场的区域）很小（这种情况我们称为小范围屈服）小范围屈服），从而对裂纹尖端场的总体影响不大。Irwin 通过研究认为在该情况下应力强度因子K K 仍有意义，仍然可以认为是K K场主导场主导着裂纹的行为。如塑性区尺寸比裂纹长度小一个数量级，工程中一般仍用线弹性理论计算应力强度因子，不过要对应力强度因子进行修正。4950小范围屈服

29、条件小范围屈服条件 o 在线弹性情况下，裂纹尖端场完全由应力强度因子K来主导，称为K主导区主导区，其尺寸 取取决于裂纹和构件的几何形状决于裂纹和构件的几何形状。o 但如果考虑裂纹尖端的弹塑性性质，则在裂纹尖端存在一个塑性区塑性区，其尺寸记为 ，显然，随着载荷的增大随着载荷的增大，越来越多的材料发生屈服，即 越来越越来越大大。 prprKr51o 塑性区的存在会改变其相邻区域的场，使之偏离K场，这一明显偏离K场，但仍属于线弹性的区域将裂尖的塑性区和K场连接起来，称为过渡区过渡区。o 如果这一过渡区的尺寸与过渡区的尺寸与 相当相当，同样就不能再将K作为主导参数，K场即失去了其场即失去了其主导地位主

30、导地位。因此，要认为K仍然是裂纹断裂形为的主导参数，必须满足：0.3 0.5ptKrrra建议的小范围屈服条件建议的小范围屈服条件 Kr52Irwin 小范围屈服理论KrK场KrK场过渡区pr塑性区53o K主导区大小即 是与载荷没有明显关系与载荷没有明显关系的，而塑性区尺寸 是载荷的单调函数是载荷的单调函数。随着外载外载的增增大大，塑性区不断长大，并使K场失去其主导地位场失去其主导地位。o 工程处理上,一般认为，当外加载荷P小于0.5P0时可以认为是小范围屈服,其中是P0裂纹体达到全面屈服时的载荷。对于理想塑性材料，P0即是塑性极限载荷。o 因此，当 时认为裂纹尖端场仍由裂纹尖端场仍由 K场

31、所场所主导主导，所有外载及几何信息仍可通过K来反映，它决定着裂尖附近的塑性区尺寸和塑性变形的大小。Krpr00.5PP54前面已指出，裂纹尖端应力场是渐进解渐进解，仅仅适合于裂纹尖端附近： 小范围屈服是指： o塑性区形状的估算： 作下列假定： （1）忽略裂纹尖端材料屈服后对塑性区外K场的影响； （2）材料为理想塑性，且遵循 Von-Mises 屈服条件。 arrK02. 0arrKp02. 0) 1 ()2(Irwin 小范围屈服理论00时，按线弹性解2IyKr55CRACK TIP PLASTICITYFirst approximation:Better approaches-selecte

32、d shape: better size estimation-Irwin-Dugdale-Better shape but first order approximation for the sizeIrwin approach: stress redistribution; elastic plastic; plane stress221YSyKryprr*256确定塑性区尺寸的确定塑性区尺寸的Irwin理论理论 o在不考虑裂纹尖端塑性影响在不考虑裂纹尖端塑性影响的情况下，线弹性裂纹尖端的K场分布为：o K场的表达式 2Kr先假设先假设裂纹尖端塑性区的存在不致改变其周围的应力场，不引起应力

33、松驰，即没有过渡区没有过渡区的存在。则只要将上式代入屈服条件，即可以得到塑性区的尺寸和形状。1122123311223cos1 sinsin2223cos1 sinsin2223cossincos2220KrKrKr 222平面应力平面应变57Von Mises屈服条件屈服条件 o几种最常用形式为一般形式：平面应力：平面应变：2222212233132ijijssS S1或213ijijkkijS 为偏斜应力12, 3为主应力屈服应力22221122223333111223316662s2211221122123s 22222112211221211 232s 58根据材料力学主应力求解公式得

34、到型裂纹的应力公式 122()22xyxyxy12cos1 sin222Kr30平面应力 59Irwin 小范围屈服理论o可以得到塑性区尺寸可以得到塑性区尺寸rp 为为： 对平面应力情况 对平面应变情况 o在裂纹延长线上塑性区尺寸在裂纹延长线上塑性区尺寸ro为为： 对平面应力情况对平面应变情况)2sin31 ()2cos(2122sIpKr)2sin3)21()2cos(21222sIpKr2)(21sIoKr22)(2)21 (sIoKr)3()4()5()6(据上式画出据上式画出 曲线，如下图中实线所示。这条闭合曲线表示裂曲线，如下图中实线所示。这条闭合曲线表示裂尖附近塑性区的周边形状，曲

35、线上各点的相当应力等于屈服极限，尖附近塑性区的周边形状，曲线上各点的相当应力等于屈服极限，内部各点超出屈服极限，未考虑应力松弛效应。内部各点超出屈服极限，未考虑应力松弛效应。r3()2cos22xyKr 60Irwin 小范围屈服理论由图可见，平面应变的塑性区远比平面应力的小，原因是：平面应由图可见，平面应变的塑性区远比平面应力的小，原因是：平面应变状态下，沿厚度方向约束所产生的变状态下，沿厚度方向约束所产生的 是拉应力，在三向拉伸应是拉应力，在三向拉伸应力状态下，材料不易屈服而变脆。力状态下，材料不易屈服而变脆。 z61确定塑性区尺寸的确定塑性区尺寸的Irwin理论理论 62说明说明o 平面

36、应力与平面应变的塑性区形状不同塑性区形状不同。这样的形状容易从其应力状态的差异想象出来。o 平面应变的塑性区尺寸（在同样的 下）小于小于平面应力的塑性区尺寸。o 例如， 平面应变情况下的 仅是平面应力的16%。这是因为在平面应变情况下，裂尖材料承受的是三轴拉伸应力状态三轴拉伸应力状态，而Von Mises屈服条件（以及Tresca条件）认为静水应力不影响屈服。K0.30r0r63确定塑性区尺寸的确定塑性区尺寸的Irwin理论理论 64First Order Aproximations of Plastic Zone Shapes Plastic zone shape from Von Mise

37、s yield criterionThrough-thickness plastic zone in a plate of intermediate thicknessEmpirical Rules to estimating Plane Stress vs. Plane Strain conditions: -Plane Stress: 2.ry B -Plane Strain: 2.ry 1/10 B65利用Tresca屈服条件 在复杂受力下,当最大切应力等于材料弹性拉伸时的屈服切应力,材料即屈服.比较发现:平面应变塑性区尺寸小,平面应变处于三向拉伸状态不易屈服.平面应变的有效屈服应力 比

38、 高 yss塑性区中的最大应力 1ys平面应变 13yss32 2ys平面应力 1yssTrescaTresca屈服条件下的塑性区尺寸屈服条件下的塑性区尺寸66Irwin 小范围屈服理论o可以得到塑性区尺寸可以得到塑性区尺寸rp 为为： 67平面应力平面应变68Plastic zone shapes for sliding mode and tearing modes69o It was mentioned that it is extremely difficult to properly describe size and shape of the plastic zone at the

39、same time.Plastic zone appearance on the front surface, back surface and a normal section of a notched silicon iron specimen in plane stress70应力松弛的修正应力松弛的修正在上面的分析中，我们假设塑性区不影响其周围的应力分布。即未考虑塑性区内塑性变形引起的应力松弛，即应力再分布影响。应力再分布影响。这样，就相当于将奇异的K场在裂纹前的塑性区的简单地用Von Mises屈服应力代替。因此，上面给出的塑性区尺寸塑性区尺寸的解显然无法无法满足总体静力平衡方程满足

40、总体静力平衡方程。Irwin认为，我们可以将塑性区尺寸的增大到某一值将塑性区尺寸的增大到某一值，使总体的静力平衡方程得到满足。711x2xaAABBCC0r0Rys22eff或BC BCACA B C BC BC22A B C AB000002200022dd22rrysKKrrKRrrKr场如图所示，虚线 表示线弹性裂纹尖端场即K场，曲线 表示考虑塑性区引起应力松驰后的应力分布，其中近似认为 是 段的简单平衡。在小范围屈服条件下在小范围屈服条件下，认为认为 下方下方的面积等于的面积等于 下方的面积下方的面积。因此，要使裂纹前方延长线上的应力与外载相平衡，就要求应力松驰后的曲线 与线弹性的K场

41、下面 的面积相等，即： 应力松弛的修正应力松弛的修正72对于平面应力情况，当 0ysrK22103跟据Mises条件 sys，即单向拉伸时的屈服极限。 把 rKy2/12202sKr， 代入上面的积分，得到在考虑应力可见，应力松弛使塑性区尺寸增加一倍。20012sssKKKRr松驰条件下，平面应力I型裂纹尖端的塑性区尺寸应力松弛的修正应力松弛的修正73同理，对于平面应变情况，当同理，对于平面应变情况，当 ， 0据据MisesMises条件条件 ysrK221ys2)(213sys211，把，把 rKy2/12220212sKr代入上面的积分，代入上面的积分，可以得到在考虑应力松驰时，平面应变

42、I型可见，平面应变状态下，若考虑塑性区应力松弛影响，塑可见，平面应变状态下，若考虑塑性区应力松弛影响，塑性区尺寸同样增加一倍。上述结果，是偏安全的近似解。性区尺寸同样增加一倍。上述结果，是偏安全的近似解。裂纹的尖端塑性区尺寸为22000220121222122sssKKrRrKr应力松弛的修正应力松弛的修正74 关于塑性区的尺寸和形状，两点补充说明：o 分析没有考虑材料强化,材料强化使裂纹尖端塑性区的尺寸变小,对于设计是偏于安全的.o 一种非常简化的分析，实际上裂纹尖端的塑性区尺寸和形状与上面的结果都有所偏差。在平面应力情况下，还有其他的塑性区形状，如窄条屈服区。o 在三维情况下，例如核电站的

43、压力容器和管道中的一个穿透裂纹，塑性区是一个三维的复杂形状。75考虑应力松弛修正后的考虑应力松弛修正后的Ko在考虑塑性区修正后，裂纹前方的应力场变成了的 分布，也就是说在塑性区以外的应力场相等于向前移动了（ ）的距离。因为 ，所以Irwin建议将裂纹尺寸进行如下的修正:o等效裂纹的尖端在屈服区的中心，它由修正裂尖的K场所包围。如果在线弹性情况下，K表示为: A B C 00Rr002Rreff0aar等效裂纹尺寸或等效裂纹尺寸或当量裂纹尺寸当量裂纹尺寸 处称为物理裂纹尖端，处称为物理裂纹尖端， 处称为虚设裂纹尖端处称为虚设裂纹尖端 0r 0rr,KKa G0,KKar G修正后修正后76yra

44、a 称为等效裂纹长度等效裂纹模型法等效裂纹模型法指以 代替原裂纹长，对应力强度因子进行修正。这说明，塑性区的存在相当于裂纹长度增加，即裂纹体的柔度增加，因而裂纹的应变能释放率也增加。a在引入小范围屈服情况下等效裂纹长度等效裂纹长度的概念后，线弹性断裂力学中的应力强度因子理论仍然有效。只要将应力强度因子K中的裂纹长度用等效裂纹长度代替即可等效裂纹长度与应力强度因子等效裂纹长度与应力强度因子77等效裂纹长度与应力强度因子等效裂纹长度与应力强度因子应力强度因子应力强度因子 裂纹尖端应力场强弱的标志。裂纹尖端应力场强弱的标志。 取等效裂纹长度取等效裂纹长度 ，令等效裂尖附近应力场的，令等效裂尖附近应力

45、场的线弹性理论分布曲线在原裂纹塑性区边界线弹性理论分布曲线在原裂纹塑性区边界C1 C1 yraa即在即在 处的应力等于处的应力等于 yrRrysysy又因为又因为 ，所以有，所以有 rKy2ysyrRK)(2222ysyKRrK应力松弛后的应力强度因子。应力松弛后的应力强度因子。 78平面应力下平面应力下 ，有，有sys221sKR代入上式，并近似设代入上式，并近似设 KK得得 2221syKr平面应变下，按平面应变下，按 222)21 (1sKRsys)21 (1则则 222)21 (21syKr 等效裂纹长度与应力强度因子等效裂纹长度与应力强度因子2yRr裂纹的计算边界正好在塑性区的中心

46、79另外，若按一般采用的公式另外，若按一般采用的公式 2221sKRsys22则：则： 2241syKr继而按等效裂纹长度计算等效应力强度因子继而按等效裂纹长度计算等效应力强度因子 ，一，一般工程应用中，取般工程应用中，取 ，又因，又因 ，用等，用等效裂长效裂长 代替代替 ，则有：，则有：K KKaYKyra ayraYK对于平面应力情况，代入相应的对于平面应力情况，代入相应的 ，得，得yr2221sYaYK等效裂纹长度与应力强度因子等效裂纹长度与应力强度因子80同理可得平面应变状态下应力强度因子同理可得平面应变状态下应力强度因子22241sYaYK 可见两种状态下应力强度因子都扩大。上述结论

47、都可见两种状态下应力强度因子都扩大。上述结论都是近似的，我们假设了是近似的，我们假设了 ，且未考虑等效裂长对形，且未考虑等效裂长对形状因子状因子Y Y的影响。对于复杂问题要用逐次逼近法求的影响。对于复杂问题要用逐次逼近法求 ，具，具体步骤见书。体步骤见书。 KKK等效裂纹长度与应力强度因子等效裂纹长度与应力强度因子81o前面我们已经有 I 型裂纹应力强度因子 KI 的表达式 其中Y(g，a) 称为几何影响因子。o引入小范围屈服情况下的等效裂纹长度， o由于o而塑性区长度 ro又依赖于应力强度因子 KI，所以，考虑小范围屈服修正后的应力强度因子需要迭代计算。 aagYagYKI),()(effe

48、ffIaagYK),(oeffraa00,0.5 0.5saRPP或Irwin的上述理论是在小范围屈服的条件下建立的，即要求:等效裂纹长度与应力强度因子等效裂纹长度与应力强度因子82如何而来？如何而来？因为平面应变下，沿板厚在第三向拉应力因为平面应变下，沿板厚在第三向拉应力 ,三向拉伸三向拉伸应力作用下应力作用下,材料不易屈服材料不易屈服,即材料的有效屈服应力即材料的有效屈服应力 比比单向拉伸屈服应力要高单向拉伸屈服应力要高,而平面应力条件下而平面应力条件下,有效屈服应有效屈服应力力 . 3ys33ys下面进行证明：下面进行证明：设设 是最大主应力是最大主应力, , ,代入代入mises准则准

49、则1max21n31m22222max(1)()(1) 2snnmm设塑性约束系数设塑性约束系数 ,代入上式有代入上式有maxscss=ys122max2(1)scnmmnnm83对对型裂纹平面应变型裂纹平面应变: 211 sin21 sin2n31211()21 sin2m 在x轴上, 01,2nmmax121sc若取 ,则 ,即133c max3s对平面应力对平面应力30,0,1mnmax1scmaxs84把塑性区中最大应力把塑性区中最大应力 叫做有效屈服应力用叫做有效屈服应力用 表示表示,1maxys13()1 2()ssyssplane strainplane stress表面平面应变

50、在表面平面应变在 的平面上的平面上,屈服区内最大应力屈服区内最大应力 是是 的三倍的三倍.0yss实际一般试件表面是处于平面应力实际一般试件表面是处于平面应力,只有中心部分才是平面应只有中心部分才是平面应变变,故平均约束系数故平均约束系数 ,实验测定实验测定 ,用环形切口圆用环形切口圆棒试件所做的拉伸试验棒试件所做的拉伸试验,在三向拉伸状态下：在三向拉伸状态下：3c 1.52.0c 一般取一般取1.72 2ysss2 2()()syssplane strainplane stress , ,85复习复习- -线弹性断裂力学线弹性断裂力学 线弹性材料的断裂准则线弹性材料的断裂准则-应力强度因子断

51、裂准则：应力强度因子断裂准则： 条件：塑性区比条件：塑性区比K场区小得多，而场区小得多，而K场区又比裂纹长场区又比裂纹长度小得多度小得多 , crcrKKGG或失稳断裂86o 用柔度法确定临界应变能释放率用柔度法确定临界应变能释放率 柔度：变形与载荷的比值柔度：变形与载荷的比值 总应变能总应变能柔度：柔度： 应变能释放率：应变能释放率： 临界应变能释放率：临界应变能释放率：CGcF21122VFcF212cVcGFaa212crcrFcGba复习复习-线弹性断裂力学线弹性断裂力学87Linear Elastic Crack-tip Fields Mode I:(general case)88Mode III:Mode II:89Angular distributions of crack-tip stresses for the three modes (rectangular: lef