损伤与断裂力学第4章(矿大)课件

第4章线弹性断裂力学的基本理论第4章线弹性断裂力学的基本理论14-1应力强度因子概念和能量释放观点的统一 假设不考虑塑性变形能、热能和动能等其它能量的损耗，则能量转换表现为所有能量在裂端释放以形成新的裂纹面积。下面以带有穿透板厚的I型裂纹的平板为例，来建立应力强度因子和能量释放率间的关系。4-1应力强度因子概念和能量释放观点的统一假设2裂纹尖端正前方的应力分布 裂纹长度(或裂纹半长度)为a的裂纹端点正前方r处有使裂纹面撑开的拉伸应力:裂纹尖端正前方的应力分布裂纹长度(或裂纹半长度)为3裂纹面上的位移 在初始应力如上式给出的情形下，设裂纹可以延长a长度，即把裂端前方撑开成长度为a+s的裂纹。此时在原坐标系的x=r处或离新裂纹端点s-r处，新裂纹上表面的位移v(s-r,)：裂纹面上的位移在初始应力如上式给出的情形下，设裂纹可4裂纹形成时外力做功当当裂裂纹纹表表面面张张开开至至上上式式给给出出的的位位移移值值时时，裂裂纹纹表表面面才才真真正正形形成成，此此时时裂裂纹纹表表面面已已无无应应力力作作用用。由由于于作作用用力力与与位位移移同同向向，当当裂裂纹纹长长度度延延长长s时时，作作用用力力对对裂裂纹纹上上表表面面所所做做的的功为：功为：B为平板的厚度裂纹形成时外力做功当裂纹表面张开至上5能量释放与应力强度因子 按照Griffith能量释放的观点，裂纹长度延长s时，此裂纹端所释放的能量将等于裂纹上下表面所做的功。因此，按照I型裂纹能量释放率GI的定义:当s0时，有KIa+sKI，经过积分得：能量释放与应力强度因子按照Griffith能量释6简化的统一形式对于平面问题，若取有效弹性模量对于平面问题，若取有效弹性模量E1和有效泊松比和有效泊松比1，而，而平面应力平面应力平面应变平面应变则平面应力和平面应变状态下的本构关系及其解均将完全相同则平面应力和平面应变状态下的本构关系及其解均将完全相同。于是：于是：这就是这就是著名的能量释著名的能量释放率与应力强度因子放率与应力强度因子的关系式。的关系式。简化的统一形式对于平面问题，若取有效弹性模量E1和有效泊松比7其余裂纹类型情况上式成立是基于裂纹沿原方向扩展的假设上。若假设II型裂纹和III型裂纹的扩展方向也是裂端正前方，则有：如果带裂纹的平板受到I，II，III三种载荷而成复合型裂纹时，若仍假设裂纹沿原方向扩展，则总能量释放率为：其余裂纹类型情况上式成立是基于裂纹沿原方向扩展的假设8提示实实验验结结果果指指出出，除除I型型裂裂纹纹可可以以沿沿原原方方向向扩扩展展外外，其其余余裂裂纹纹型型往往往往不不沿沿原原方方向向扩扩展展。因此总能量释放率只是近似估计式。如果要考虑裂纹真正的扩展方向来计算，这已不是解析的方法所能做到，必须要用数值解法，同时还要一套断裂理论指出裂纹开裂的方向。?提示实验结果指出，除I型9练习思考题思考题 为为何何不不同同型型裂裂纹纹的的能能量量释释放放率率可可以以线线性性叠叠加加?而而应应力力强强度因子却不能度因子却不能?习题习题 试试建建立立IIII型型和和IIIIII型型裂裂纹纹的的应应力力强强度度因因子子和和能能量量释释放放率率之间的关系之间的关系。练习思考题104-2柔度法 柔度法是通过柔度随裂纹长度而改变这个性质，用测量的方法来得到G，然后再利用G与K的关系来得到K值。由于I型裂纹的G与K的关系式是精确的，并且I型裂纹容易施加载荷，所以柔度法一般只用在I型裂纹。I型裂纹是最危险的。4-2柔度法柔度法是通过柔度随裂纹长度而改变这11恒载荷和恒位移时G的表达式当边界是给定位移时，外界对系统不做功，则这是恒位移时能量释放率的表达式。若系统边界某范围是给定载荷，这是恒载荷时的能量释放率表达式。柔度法一般限制在二维问题，尤其是I型裂纹，柔度法通常用来做应力强度因子的标定恒载荷和恒位移时G的表达式12恒载荷柔度法一块很长的矩形板，板厚为B，板下边固定，上边某点有拉力P，载荷点位移为。拉力P方向垂直裂纹面。在裂长为a时，拉力P可产生位移(a)，当裂纹增至时，位移也增至:因为位移和拉力有如下关系,C是柔度,于是：恒载荷柔度法一块很长的矩形板，板厚为B，板下边固定13恒位移柔度法长矩形板如图，一边固定，另一边强迫作位移 后也加以固定。假设裂纹长度 由增至 ，则应变能的改变是OBC减去OAC，即等于负的OAB的面积，。此时，载荷改变量P也是负值。恒位移柔度法长矩形板如图，一边固定，另一边强迫作位移后14应该指出，载荷与位移之间的正比关系是建立在材料服从虎克定律的基础上，因此要求裂纹尖端的塑性区是微小的，如果裂端塑性区大到不可忽略时，则载荷与位移之间的线性关系不再成立，这时属于弹塑性断裂力学的范围。应该指出，载荷与位移之间的正比关系是15柔度法一般应用于恒载荷时平板的I型裂纹问题，要求裂纹前沿整齐，有相同的能量释放率。整个应力强度因子标定的步骤如下(1)选定一标准试件长条板单边裂纹试件，用薄刀片加工，制成长为a1的I型裂纹。然后材料试验机上拉伸，画出拉力和加载点位移关系线。此时关系应是线性的。柔度法一般应用于恒载荷时平板的I型裂纹16拉力值不得大到使P-关系产生非线性。下一步，再度使裂纹稍稍延长至长度a2，有在同一张纪录纸上记下此时的P-关系。如此进行至少十多次，裂纹长度已相当长时才停止。拉力值不得大到使P-关系产生非线性17(2)求出不同裂纹长度下的柔度C，柔度C是直线P-斜率的倒数。把柔度与裂长的关系画在图中。若是数据点足够多，可用最小二乘法把数据点拟合成一条多项式表示的曲线。(2)求出不同裂纹长度下的柔度C，柔度C是直线18损伤与断裂力学第4章(矿大)课件19习题1.求恒拉力下双悬臂梁试件的应力强度因子。2.有限大小弹性圆柱体含有位于中心的圆裂纹时，试求在拉力作用下类似的能量释放率表达式。3.如图所示，在上下表面，有位移v=v0和xy=0，即被刚体固接。此时横向位移u不受约束。求在平面应力时的应力强度因子。习题204-3断裂判据断裂过程区 断裂总是始于裂端的极小区域，当其损伤达到临界程度时才发生的。在此小区域中材料的微结构起决定影响，也是宏观力学不适用的地方。这个小区域就叫做断裂过程区FPZ(fractureprocesszone)。4-3断裂判据断裂过程区21K场区 在第三章中，给出各型裂纹的裂端应力场时，已忽略掉高次项，因此也仅适合裂纹尖端的小区域内，此区域称为K场区。K场区内的应力应变强度可用应力强度因子来度量；场区外则须加上高次项。K场区在第三章中，给出各型裂纹的裂端应力场时，已22关于K场区和断裂过程区 如果K场区尺寸小于断裂进行区尺寸，则计算应力强度因子已失掉意义，此时宏观力学在裂端区是不实用的。反过来，若K场区尺寸比断裂进行区尺寸大几倍以上，则断裂进行区是否会发生断裂，受其外部的K场区强度所制约，因此，断裂判据可建立在K场区强度是否达到临界条件这个基础上。关于K场区和断裂过程区如果K场区尺寸小于断裂进23 由于无限大应力实际上不存在，裂端总有个塑性区，而塑性区内的应力是有界的。因此，应力强度因子断裂判据成立的条件是，塑性区尺寸比K场区小几倍，也要比裂纹长度小几倍以上。幸运的是，许多高强度合金和工程材料在发生脆性断裂时，多是K场区强度起支配作用的。因此，应力强度因子断裂判据适合于这些材料的脆性断裂。由于无限大应力实际上不存在，裂端总有个塑性区，而塑24断裂判据 对于一个单独型的裂纹，利用应力强度因子和能量释放率的对于一个单独型的裂纹，利用应力强度因子和能量释放率的关系，可有断裂判据：关系，可有断裂判据：I I型型裂裂纹纹是是最最常常见见的的裂裂纹纹型型，其其失失稳稳断断裂裂开开始始的的临临界界点点K Kcrcr，通通常常与与试试件件(或或构构件件)的的厚厚薄薄、大大小小有有关关。当当试试件件(或或构构件件)厚厚到到某某一一程程度度和和大大到到某某一一程程度度，脆脆性性材材料料的的K Kcrcr值值达达到到极极小小值值，以以后后尺尺寸寸厚厚度度再再增增加加，K Kcrcr仍仍维维持持此此极极小小值值，此此极极小小值值用用符符号号代代表表即即为为K KICIC ，其其相相应应的的G GICIC值值称称为为平平面面应应变变的的断断裂裂韧韧度度。因因此此，I I型型裂裂纹纹保保守守的的判据为判据为:断裂判据对于一个单独型的裂纹，利用应力强度因子和25断裂判据可以解决下列两个问题(1)当知道工作载荷时，可以计算出断裂时的临界裂纹尺寸；(2)当知道裂纹尺寸和位置时，可计算出可能引起断裂的载荷。断裂判据可以解决下列两个问题(1)当知道工作载荷时，可以计算26例题例例题题34GrNi3Mo钢所制成的粗轴，探伤检查发现主要的缺陷是内部有一半径为40mm的圆裂纹，裂纹面的法线方向与轴向平行。已知轴半径远大于裂纹尺寸，同时测得钢的KIC为99.2MN/m3/2,试问要是发生断裂，轴向拉伸应力至少有多大?例题例题34GrNi3Mo钢所制成的粗轴，探伤检查27解答因为粗轴半径远大于圆裂纹半径，可采用无限大弹性体有圆裂纹的应力强度因子的解。于是临界条件为：于是：解答因为粗轴半径远大于圆裂纹半径，可采用无限大28习题1 1 有有一一平平均均半半径径为为800mm800mm，壁壁厚厚10mm10mm的的圆圆柱柱壳壳形形高高压压容容器器，所所用用钢钢材材具具有有K KICIC值值130MN/m130MN/m3/23/2,现现发发现现有有一一长长为为12mm12mm、最最深深为为4mm4mm的的内内表表面面半半椭椭圆圆轴轴向向裂裂纹纹。问问(A)(A)内内压压多多达达时时才才会会发发生生断断裂裂?(B)?(B)若若内内压压只只是是临临界界内内压压的的1/41/4，并并假假设设裂裂纹纹形形状状不不变变，问问裂裂纹纹扩扩展展至至多多大大才才会会发发生生断裂断裂?2 2 在在习习题题1 1中中，若若发发现现的的裂裂纹纹是是长长为为200 200 mmmm，深深为为4 4 mmmm的的内内表表面面横横向向(圆圆周周向向)裂裂纹纹。试试求求出出习习题题1 1中中的的两两个个问题。问题。习题1有一平均半径为800mm，壁厚10mm的294-4阻力曲线能量释放率可做为裂纹是否扩展的倾向能力的度量，又称为裂纹扩展力。裂纹扩展力必须大于裂纹扩展阻力，裂纹才有可能扩展。对平面应变的脆性断裂来说，裂纹扩展阻力由KIC确定，是个常数值，不随裂纹增长而变。但对不同厚度的平板，尤其是厚度小于平面应变所要求的厚度时，裂纹扩展阻力不再是常数。为了说明裂纹扩展阻力的观念，现在以平面应变无限大平板I型中心裂纹为例，4-4阻力曲线能量释放率可做为裂纹30脆性断裂阻力曲线当拉伸应力保持定值时，裂纹扩展力G随a增加而线性上升。在1时，裂纹半长度为a1就达到裂纹扩展阻力值GIC。超过a1，就发生失稳断裂；低于a1，则裂纹不扩展。脆性断裂阻力曲线当拉伸应力保持定值时，裂纹扩31脆性断裂阻力曲线以小于1的拉伸应力2作用时，必须超过较长的a2才会发生断裂。图中带箭头的直线代表裂纹扩展力，只有当裂纹扩展力大于常数值的阻力R=KIC，才会发生失稳断裂。脆性断裂阻力曲线以小于1的拉伸应力2作用时32脆性断裂阻力曲线如果将x轴改为代表裂纹扩展增量a，则可以改画成下图。a0部分才是真正扩展。a0部分即表示不扩展，而以负方向离原点的距离表示裂纹半长度的大小。脆性断裂阻力曲线如果将x轴改为代表裂纹扩展增量33韧性断裂阻力曲线在板厚较薄而不合乎平面应变条件时，裂纹扩展阻力R随a增加而增加。图(49)的例子仍是Griffith裂纹，此时裂纹扩展阻力是一曲线。此曲线叫做阻力曲线或R曲线。图4-9非平面应变的R曲线韧性断裂阻力曲线在板厚较薄而不合乎平面应变条件时34韧性断裂阻力曲线例如在韧性断裂时，裂纹扩展阻力往往是呈曲线的。一旦达到并稍为超过裂纹开始扩展的条件时，若外力仍维持不变，则较长的裂纹(例如图中的裂长a2受到2作用时)有可能稍为扩展，然后很快地停止下来。只有当外力较大时，才有可能引起失稳扩展。图4-9非平面应变的R曲线韧性断裂阻力曲线例如在韧性断裂时，裂纹扩展阻力往35韧性断裂阻力曲线阻力曲线的测定一般是针对裂纹扩展阻力不为常数值时才实施。脆性材料平面应变的恒载荷试验时，试件一启裂就立即失稳扩展。但对阻力随裂纹扩展增量而变的情形，达到启裂点后不一定会发生扩展，即使扩展也不一定是失稳扩展。当扩展力稍稍超过启裂点时，往往有一段稳定扩展(也叫做亚临界裂纹扩展)。当达到失稳断裂时，这时的a量已达到不可忽略了，对于有稳定扩展阶段的断裂韧度测试中，若监测启裂点不容易时，可以用阻力曲线的测量，然后用外推法得出启裂点。韧性断裂阻力曲线阻力曲线的测定一般36思考题1.在带裂纹平板的表面铆接或焊接一块平板，把裂纹覆盖起来(此称为加筋板)。试根据阻力曲线的原理，讨论裂纹可能的扩展行为。2.如果已测量到一组(R,a)数据，这里R是裂纹扩展阻力。你如何得到启裂时的R值?3.图(49)中，若裂纹很短(例如为a1)，裂纹扩展力为G1的概念在什么情况下是正确的?什么情况下是不正确的?思考题374-5应变能密度因子 考虑二维的裂纹问题，受到I、II、III型三种载荷中的任一种或两种以上载荷的作用。裂纹前缘是平直的，即整个前缘各点的应力强度因子值都相同，如图所示，裂纹端点区附近的一点P处有体积元，其应力场为三种裂纹应力场的叠加：4-5应变能密度因子考虑二维的裂纹问题，受到I38一般情况下的裂纹尖端应力场一般情况下的裂纹尖端应力场39平面应力平面应变记记平面应力记40应变能密度公式于是，平面应变于是，平面应变时在时在P P点的应变能点的应变能密度为密度为 ：应变能密度公式于是，平面应变时在P点的应变能密度为：41式中：式中：42FPZ 设裂端有个以裂端为原点、半径为r0的圆形损伤核(或叫断裂进行区)。设r0值远小于K场区尺寸，如此，在脆性断裂时，断裂是否会发生要由K场区应力应变强度来决定。对于复合型裂纹，应变能密度可综合度量K场区应力应变的强度；更重要的是，应变能密度是个单参数，它代替两个以上的应力强度因子(多参数)。因此，根据应变能密度的概念而建立的断裂判据，显然要比多参数判据简单。FPZ设裂端有个以裂端为原点、半径为r0的圆形损43K场区 损伤核周界是损伤核与K场区的交界。在K场区所有位置的应变能密度中，周界上的应变能密度对断裂是否发生，起着决定性的作用。K场区损伤核周界是损伤核与K场区的交界。在K场区44应变能密度因子S 此应变能密度因子只是极坐标的函数，与另一变数r无关。应变能密度因子S45Sih(薛昌明)假说Sih(薛昌明)提出下列两个假说（1）裂纹扩展的方向为S的一个局部极小值的方向，即这里0为裂纹扩展角，或叫做开裂角。（2）当此S极小值，即Smin=S(0)，达到或超过一临界值Scr时，就发生失稳断裂。Sih(薛昌明)假说Sih(薛昌明)提出下46薛昌明介绍国际著名学者、美籍华裔科学家、美国Lehigh大学资深教授、现代工程断裂力学的奠基人之一，Lehigh大学断裂与固体力学研究所所长。现担任华东理工大学特聘兼职教授等。薛昌明介绍国际著名学者、美籍华裔科学家、美国Le47生于1953年获美国Porland大学机械工程学士学位，1957年获美国纽约大学机械工程硕士学位，1960年获美国里海大学机械工程博士学位。1965年，任美国里海大学教授。从1965年开始，薛博士先后兼任美国加州理工学院客座教授（1965）、美国Hahnemann医学院和费城医院兼职教授（1972）、西安交通大学名誉教授（1995）、清华大学客座教授（1995）、中国科学院荣誉研究员（1995）、西安交通大学高科技研究院院长（1996）、中国科技大学客座教授（1997）等。生于1953年获美国Porland48主要兼职有：国际理论与应用断裂力学杂志主编；国际断裂力学杂志地区编委；美国ASME应用力学杂志副主编；工程断裂力学杂志、断裂杂志、热传导杂志和先进复合材料杂志的编辑顾问委员会成员；美国机械工程师协会（ASME）理事；国际断裂协会名誉理事；国际断裂联合会和国际细观力学学会的创始人。薛昌明教授编写出版了六本学术专著，编撰出版了6个系列31卷的科学技术系列丛书，编撰出版了24次学术会议论文集，获得了3项美国专利，发表了350余篇学术论文。主要兼职有：国际理论与应用断裂力49疑问 上述的假说给出了裂纹的开裂角和失稳断裂发生的条件，但损伤核究竟有多大?真的是圆形吗?这些都是薛昌明避而不谈的地方。此外，有局部极小值的角度可能不只一处，裂纹开裂的方向究竟其中哪一方向?是不是最小的那个方向?从物理意义讲，应变能密度最小的地方，应该是最稳定的地方，因为可以释放的多余的能量较少。而应变能密度最大处则有更多的能量可释放，因此应该比较不稳定，比较容易从该处开裂。基于上述的争论，薛昌明对此理论做了一番补充。疑问上述的假说给出了裂纹的开裂角和失稳断50薛昌明的补充说明考虑到损伤核尺寸甚小，该区域的塑性变形相当大，应变能密度可以分为两部分，一部分是体积膨胀能密度，另一部分是歪形能密度。决定脆性开裂的是体积膨胀能密度，决定塑性失稳的则是歪形能密度。因此，在S的几个局部极小值的角度处，裂纹应沿体积膨胀能密度较大的方向开裂。在S最大值处，可以验证此处歪形能密度是占支配地位的。薛昌明的补充说明考虑到损伤核尺寸甚小51体积膨胀能和歪形能考察应变能密度，可以分静水应力引起的体积膨胀能密度和歪变形引起的歪形能密度两部分，如图所示。因此应变能密度可分为：下标V代表体积膨胀部分，D代表歪形部分。体积膨胀能和歪形能考察应变能密度，可以分静水应力引52由弹性力学公式可得：薛昌明认为I型裂纹开裂方向，虽说是Smin在所处方向，但此时(dU/dV)V占支配地位。由弹性力学公式可得：53S因子理论的验证现在考虑平面应变I型裂纹来验证S因子理论是否恰当。对于I型裂纹：即：于是：时，S有极小值。时，S有极大值。S因子理论的验证现在考虑平面应变I型裂纹来验证S因子理论是否54在损伤核周界，应变能密度因子分为膨胀和歪变形相关的两部分，即 对I型裂纹，在 时：这里若取=1/3，则有:故平面应变I型裂纹的在薛昌明应变能密度因子理论所确定的开裂角处膨胀变形能占支配地位。在损伤核周界，应变能密度因子分为膨胀和歪变形相关的两55习题1.试用应变能密度因子的假说求平面应变II型裂纹的可能扩展方向。问此时SV是否仍占支配地位?在无限大平板有相同的中心裂纹时，问纯I型裂纹和II型裂纹的临界应力比cr/cr为多少?2.试求平面应变I型裂纹的SV表达式，问其极大值是否在=0的方向?习题试用应变能密度因子的假说求平面应变II型裂纹的可564-6平面I-II复合型裂纹S因子判据：4-6平面I-II复合型裂纹S因子判据：57S因子理论由Sih的因子理论可确定裂纹扩展方向。现在以无限大平板有穿透板厚的斜裂纹为例，来确定可能的裂纹扩展方向和断裂时机。S因子理论58S因子理论在xy坐标系下，无穷远处的应力分量为:y=,x=xy=0。在新坐标系xy下，应力分量变为:这是I、II型复合裂纹问题，其应力强度因子为：S因子理论在xy坐标系下，无穷远处的应力分量为:y=,59S因子理论确定裂纹扩展方向的方程：此时：由S因子第二个假说，断裂判据为:S因子理论确定裂纹扩展方向的方程：60S因子理论 对于不同的角，表给出上述例题的开裂角0，表中数据是在=1/3的情形下得到的（单位:度）。S因子理论对于不同的角，表给出上述例题的开裂角0，表61最大周向应力理论 与S因子假说比较，最大周向应力理论显得比较简单，且容易被接受，此理论假设裂纹是沿着裂端区圆形损伤周界的最大周向应力所处位置的方向开裂。现在仍以I和II型复合型来说明：裂端区周向应力和剪切应力可由应力转换而得到：最大周向应力理论与S因子假说比较，最大周向应力理论显得比62对求偏导，得：确定开裂角的方程为：即：所得结果与S因子理论相差不多。失稳判据为损伤核周界的最大周向应力大于等于材料的临界应力（强度）：对求偏导，得：63最大能量释放理论 在Griffith能量释放理论中，已假设裂纹是沿着原来的方向扩展，这对I型裂纹是正确的，但对II型裂纹或复合型裂纹，实验已证明裂纹并不沿原来方向扩展。因此前面所建立的能量释放率与应力强度因子的关系式，只有在I型的情况下才正确，其余都只是近似的结果。最大能量释放理论在Griffith能量释放理论中，已假设64 按照物理观点，裂纹扩展的方向应是有最多能量可释放的方向：由此确定以及发生的方位。失稳断裂判据取为：按照物理观点，裂纹扩展的方向应是有最多能量可释放的方向65最大能量释放率理论的物理意义是很明显的，提法是比较恰当的。问题是对II型或复合型，要找出的方向只有通过复杂的数值解法来进行，而不象S因子理论和最大周向应力理论那样，只要知道应力强度因子，问题就已解决。因此在工程应用中，在I、II复合型时，使用保守的断裂准则：再保守一点：最大能量释放率理论的物理意义是很明显的，提法是比较恰66复合型裂纹工程简化断裂判据复合型裂纹工程简化断裂判据67习题1有一半径很大的薄壁圆柱壳形压力容器，发现有一穿透壁厚、长度为2a的裂纹。裂纹位向与轴向成角。设内压为p，壁厚为t，半径为R。(1)试用最大周向应力理论给出决定开裂角的方程。(2)设KIC=183MN/m3/2，=45o，a=100mm，R=1000mm，t=30mm，试用最大周向应力理论给出可能的断裂时内压。习题1有一半径很大的薄壁圆柱壳形压力容器，发现68