应力强度因子.docx

1、应力强度因子第二章 应力强度因子的计算K-应力、位移场的度量=K的计算很重要,计算K值的几种方法：1.数学分析法：复变函数法、积分变换；2.近似计算法：边界配置法、有限元法；3.实验标定法：柔度标定法；4.实验应力分析法：光弹性法. 2-1 三种基本裂纹应力强度因子的计算一、无限大板I型裂纹应力强度因子的计算Ki = lim , Zi 计算K的基本公式，适用于型裂纹.1.在“无限大”平板中具有长度为2a的穿透板厚的裂纹表面上，距离X = _b处 各作用一对集中力p.；x 二 ReZi - y Im Zi；y 二 ReZi y Im Zixy = _y ReZi选取复变解析函数:2 pz a2

2、b2二(z2 _b2)边界条件:b.z va,出去z = b处裂纹为自由表面上=0,ixy =0。c.如切出xy坐标系内的第一象限的薄平板，在 x轴所在截面上内力总和为p以新坐标表示:2 p（匕 +a） Ja2 +孑二（a）2 -b2 ； （ 2a）Ki2p、a二（a2 -b2）2.在无限大平板中,具有长度为2a的穿透板厚的裂纹表面上，在距离x= _印的 范围内受均布载荷q作用.yb.110yqn rKq 1旺x J2 a1利用叠加原理：微段 集中力qdx dKi二2qa=dx金（a2 _x2）a 2qai . dxo - 2 2、二（a -x ）令x=acos：= a2-x2 = acosv

3、 , dx = acos当整个表面受均布载荷时，c a.=Ki = 2q：s in，(%)=qVa3.受二向均布拉力作用的无限大平板，在x轴上有一系列长度为2a，间距为2b 的裂纹.4.边界条件是周期的：a.门=；丁 .b.在所有 裂纹 内部 应力 为零.y =0, -a : x : a, -a _ 2b ： x ： a _2b在区间内-y =0 xy = 0采用新坐标：=z-a当 t 0时,sin匕=Jcos =12b 2b 2bJL 迟 JL 乜 JL JL 乜 JL= sin( a) =sin cos一a cos一 sin a 2b 2b 2b 2b 2b2a 1若裂纹间距离比裂纹本身尺

4、寸大很多（兰乞丄）可不考虑相互作用，按单个裂纹2b 5计算二、无限大平板nm型裂纹问题应力强度因子的计算1.u型裂纹应力强度因子的普遍表达形式（无限大板）：K lim Z（ 2.无限大平板中的周期性的裂纹，且在无限远的边界上处于平板面内的纯剪切 力作用.心=帆 J2 兀Z () = i VaJtan 舒3.川型裂纹应力强度因子的普遍表达形式（无限大板）:K皿叫、厂Z（）4.周期性裂纹：2b - aK = . a tan-na 2b 2-2 深埋裂纹的应力强度因子的计算y 卓x2C1950年，格林和斯内登分析了弹性物体的深埋的椭圆形裂纹邻域内的应力 和应变，得到椭圆表面上任意点，沿y方向的张开位

5、移为：2 2 1x z .2y =y0（1 2 2）a c2（1 -2）；a其中：yo =（丘丿.-为第二类椭圆积分有Ji | 2 2=2、1-C ；a sin2 d（于仁东书丿匹 a 2二 2sin2（-）2cos2 d （王铎书丿0 c1962年,lrwin利用上述结果计算在这种情况下的应力强度因子原裂纹面又N 二 Qcos : ,x,-sin :2 2 Xi 乙-2 2a cc2片2 a2zj 二 a2c2accsin2a2 cos2假设：椭圆形裂纹扩展时，其失径的增值r与成正比.边缘上任一点p （x,z）,有:xj(r)sin 炉=(1 f?sin=(1 f)xiz = r)cos 即

6、=(1 f )z1= p (x；z), p(M,Zi)均在 y=0 的平面内. ,:2 2 -2 2 4 2 2 2 2=c x a z (i f)ac a c=新的裂纹面仍为椭圆长轴c=(i f)c,短轴a=(i f)a.=y向位移2 2原有裂纹面：二 二,上）2 =ia c yo2 2 2 2 2 2二 1 一（1 2门笃 一（1 2f）% =1 _笃一刍 2f （笃 刍）a c a c a c=2f二 y2 =2fyf=2f （1+f）2y。2 E2fy。2又 f =二 Jc2 sin2 + a2 cos2 申 ac 2 cos2 :Jc2s in2 +a2 ac设各边缘的法向平面为平面

7、应变，有:= KiJ厂昱）2 y2 ,c2sin2 : a2cos24 1 4 acyo2（1二 2）；aE二 K = （a）2（c2sin2a2 cos2$ c在椭圆的短轴方向上，即，J,有K二KImax二亠1危险部位2 0椭圆片状深埋裂纹的应力强度因子当a = c时r圆片状裂纹，二：K| 2”八旳.a2 兀 2-3半椭圆表面裂纹的应力强度因子计算、表面浅裂纹的应力强度因子半椭圆片状表面线裂纹 Ki与深埋椭圆裂纹的Ki之比等于边裂纹平板与中心裂纹平板的KI值之比K i表Ki边K i 埋 K i 中又有:O.lsi n丄兀AtanW仏=(1 WK冲其中：A-裂纹长度；W-板宽度,A 2兀A 2

8、兀A 兀A当 一 1 时 sin : ， tan :W WW W WKi 1.1Ki 中Ki表-1.1Ki埋1. 16 aKi表 一 1. Ki埋一椭圆片状表面裂纹A处的Ki值、表面深裂纹的应力强度因子深裂纹：引入前后二个自由表面=使裂纹尖端的弹性约束减少=裂纹容易扩展=Ki增大= KI（表面）=Me K（埋藏）其中：Me 弹性修正系数，应大于1,由实验确定一般情况下Me = M“ M 2其中：Mi 前自由表面的修正系数M2 后自由表面的修正系数关于Me表达式两种形式的论述1.巴里斯和薛a.0时二接近于单边切口试样M1 =1.12b. ac1时=接近于半圆形的表面裂纹 M1 =1利用线性内插法

9、M1 =1 0.12(1 -空)c利用中心穿透裂纹弹性件的厚度校正系数12B 兀a ?=M 2 = ( tan )2兀a 2BB 板厚a裂纹深度c裂纹长度当a B时M 2 : 1 =浅裂纹不考后自由表面的影响2.柯巴亚希.沙.莫斯a 2M1 =1 0.12(1 )22c2B 兀a 3M 2 = ( tan )2兀a 2B二表面裂纹的应力强度因子(应为最深点处)：Ki=Mea 2-4其他问题应力强度因子的计算一、 I . u型复合问题应力强度因子的计算复变数：z=xiy， z=x-iy取复变解析函数：x(z)二p iq ,=p1 iq1取应力函数：2丨! 呼 zx(z) zx(z)或 =Ref-

10、： (z) - zx(z)= 满足双调和 方程分析第一应力不变量:(推导过程略)J F 二匕 =4Rex(z) ex cy对于I . U型复合裂纹I 型：二x 二 ReZ|y Im Z； , ；y 二 ReZ, y Im Z,k(二x 二 y)u L 0 二 2Im Zu I = 2Im 7=I、U型复合裂纹在裂纹前端处的不变量k k(J r),十2Re2；丨”加一2 “。(Ki -iKu)| Q取复数形式的应力强度因子.K二Ki -iKu k二(JC | 广 2Re( )|飞又(J 二)=4Rex(Z)-K =吧2.、云-x(Z)若米用 z 坐标：=Z - a = K = 2-. 2二 li

11、m Z - ax (Z)选择x (z)满足具体问题的应力边界条件=这种方法利用普遍形式函数求解应力强度因子f =Fi(Z) Fi(Z) ZF4(Z) ZF4(Z) ( Fi(Z),F4(Z)为解析函数)-复变解析函数表达的双调和函数的普遍形式(或复变应力函数为普遍形式). 利用这个方法可以求解很多”无限大”平板中的穿透裂纹问题 二、有限宽板穿透裂纹应力强度因子的计算实际情况：应看成有限宽计算，必须考虑的自由边界对裂纹尖端应力场和 位移场的影响 在理论上得不到完全解通过近似的简化或数值计算方法 数值解.方法:边界配置法，有限单元法等.针对有限宽板问题：寻找一个满足双调和方程和边界条件的应力函数或

12、复变 解析应力函数.边界配置法：将应力函数用无穷级数表达，使其满足双调和方程和边界条件， 但不是满足所有的边界条件，而是在有限宽板的边界上，选足够多的点，用以确定 应力函数,然后再由这样符合边界条件的应力函数确定 K值.边界配置法：计算平面问题的单边裂纹问题，只限于讨论直边界问题. 以三点弯曲试样为例进行说明.(1)威廉氏(Williams)应力函数和应力公式Williams应力函数：j i . (1) -(r,巧八 Cj r2 cos(j-1户七 cos(. 1)可2 t满足双调和方程I4 (rj)=O.HT边界条件：裂纹上、下表面(V = _)，和y均为零.=上式满足.2在边界上的边界条件

13、的满足如下确定：在有限宽板的边界上选取足够的点, 如图,使这一点的边界条件满足=Cj1PFjL WI_P244 I_P2kJ - P2_s2其中：B -试件厚度,W-试件宽度.(2) K的计算 卫DjBj(r,RBW j j jqQ又因为当二=0时,cost =1 ,当j =1时在乘.2二r后与r无关,而当j =2,3,4：:时在乘.、2二r后与r有关,当r 0时都为零.，二（为切向）应利用边界条件确定Di,边界条件只个边界各点的应力，可利用不同的边界条件,a.应力.b. ,（ n为法向）.c. dn（3）借用无裂纹体内的边界条件求系数 Dj取含裂纹三点弯曲试样的左半段的受力状态和不含裂纹的悬

14、臂梁受力是 样的.取m个点分析，以2m有限级数代替无限级数精度足够.对于不同的点有：2mxyi DjE1j = xy1其中E1 j已知，xy1由材料力学计算.Bv j=Ki F(旦)BW W1 3 5 7 9a a； a； a： a：F( ) =11.6( )2 -18.4( )2 87.2( )2 -150.4( )2 154.8( )2W W W W W W其中s =4W标准试件，此式为美国SEM-E399规范 2-5 确定应力强度因子的有限元法不同裂纹体在不同的开裂方式的应力强度因子是不同的 .一些实验方法、解析方法都有各自的局限性，而有限元等数值解法十分有效地求解弹塑性体的应力 和位移

15、场,而应力和位移场与K密切相关,所以,可以通过有限元方法进行应力强 度因子的计算一、位移法求应力强度因子I型：u(rp) =0J丄(2k1)cos?cos岂4G Y 2兀 2 2Ki4G、应力法求应力强度因子有限元法=；y（r,0）= Ki -；y.2二 rKi r的关系曲线外推=Ki的准确值.应力法与位移法比较：利用刚度法求应力时，应力场比位移场的精度低（因应力是位移对坐标的偏 导数）三、间接法求应力强度因子（应变能释放率法）利用有限兀法确定G= Ki .四、J积分法-:围绕裂纹尖端的闭合曲线T :积分边界上的力u :边界上的位移.-hjJ 积分为：J 二Wdy -T dsp x1其中W =

16、 Gy ；iy为应变能密度.2线弹性问题：J二G = 0 .E利用有限样方法计算回路积分二K. 2-6 叠加原理及其应用Ki的叠加原理及其应用 1. &的叠加线弹性叠加原理：当n个载荷同时作用于某一弹性体上时，载荷组在某一点 上引起的应力和位移等于单个载荷在该点引起的应力和位移分量之总和 叠加原理适用于K证明： 二 lim、2二 r；y |些设在Ti载荷作用下,有：巧|世,Ki(1lim)/(1)y 0 设在T2载荷作用下，有：6|曲，Ki=四72y 由叠加原理有：；y l：_o = ；y二(2)y |甘=Ki二Ki(1) - Ki 满足叠加原理计算复杂载荷下应力强度因子的方法：将复杂载荷分解

17、成简单载荷，简单载 荷可查手册.2.实例:铆钉孔边双耳裂纹的 值叠加原理：Ki(a) =Ki(b) - Ki -Ki(d)= Ki二丄(心(5 Ki(c)2其中：Ki(b (昔) D为圆孔直径,可查应力强度因子手册2板有宽度：F(a)sec - 板宽的修正.W V W这里：af = D a 即有效裂纹长度.二(a )22=K 滿(詈)2确定Kz(c):无限板宽中心贯穿裂纹受集中力 p作用.a)、应力场叠加原理及其应用1.应力场叠加原理To:无裂纹时外边界约束在裂纹所处位置产生的内应力场叠加原理：Ki(a) = Ki(b) - Ki(c) = Ki(c)-应力场叠加原理：在复杂的外界约束作用下，

18、裂纹前端的应力强度因子等于没 有外界约束，但在裂纹表面上反向作用着无裂纹时外界约束在裂纹出产生的内应力To所致的应力强度因子如图2.实例：旋转叶轮（或轴）内孔端裂纹的以等角速度运转的叶轮，在内孔面有一长为2a的贯穿裂纹，求裂纹前段的 应力强度因子（1）求解无裂纹时，旋转体在无裂纹部位的内应力有弹性力学有：R22其中：f为叶轮密度， 为角速度，R1为叶轮内径，R2为叶轮外径，r为计算点- （平面应力）口（平面应变）一般情况下：旦 1 1 =(旦)2 R2 10 50 R,a比较小：(匸)2二1.R2二 T。f 直(1 R2)8 r根据类比原则：比较(d)与(b):内孔半径一致，裂纹大小及组态一样

19、，裂纹面上下受力一致，外边界无约束，唯一不同的是一个是有限体，一个是无限体，由于边界是自由的二宀 K(d)I I(3) 根据叠加原理(注意无裂纹时),由弹性力学知:二.K 十(d)二 K 二 K:(c) ：af (旦) 2.7 实际裂纹的近似处理利用断裂力学进行安全评价时，首先确定缺陷的大小，部位和形状，偏于安全 考虑：夹杂、空洞、气孔、夹杂性裂纹 =裂纹应针对实际问题进行分析.一、缺陷群的相互作用1.垂直外应力的并列裂纹并列裂纹的作用使下降=工程上偏安全考虑(1)并列裂纹作为单个裂纹考虑；对于密集的缺陷群，假定它们在空间规则排列，并可把空间裂纹简化成平 面裂纹.2.与外应力垂直的面内共线裂纹

20、如裂纹中心间距大于缺陷尺寸五倍以上 ，可做为单个裂纹处理，否则必须考虑修正：MW .二、裂纹形状的影响通过探伤手段=缺陷的”当量尺寸”及其部位，而缺陷的具体形状及实际 尺寸难以确定二裂纹形状的影响.1.探伤结果是面积当缺陷的面积相同时，？=丄的椭圆裂纹 心最大=以旦二丄的椭圆裂纹分析c 2 c 2是偏于安全的.2.探伤的结果是最大线尺寸(1)当最大直径相同时，圆裂纹的K比椭圆裂纹大=以圆裂纹估算偏于安 全.(2)当缺陷长度一样时,贯穿裂纹K比其它裂纹的K大=以贯穿裂纹估算 偏于安全. 2.8 塑性区及其修正小范围屈服：屈服区较小时（远远小于裂纹尺寸）.=线弹性断裂力学仍可用 一、塑性区的形状和

21、大小1.屈服条件的一般形式屈服条件:材料超过弹性阶段而进入塑性阶段的条件.a.简单情况：薄壁圆筒扭转：.一s b.复杂情况：f （二x,；yz, xy,xz, yz） =C 用主应力表示 f（G,6,6）=C有：最大正应力条件，最大切应力条件， VOn.Mises屈服条件（变形能条件）,Tresca屈服（切应力条件）.2.根据屈服条件确定塑性区形状大小a.利用米塞斯（von.mises）屈服条件.当复杂应力状态下的形状改变能密度等于单向拉伸屈服时的形状改变能密 度,材料屈服，即：2 2 2（二1-匚2）（匚2 -6） 匸3-匚1）对于I型裂纹的应力公式e ecos1_sin 2 26 =0（平

22、面应力,薄板或厚板表面）Ki2 2 8 2 0=r 2 cos 1 二 3sin s2-平面应力下，I型裂纹前端屈服区域的边界方程.当八0时,r0二丄Ki 2（1 ）22 二-s平面应变（厚板中心）二3 二J 塑性区尺寸增大，依据:单位厚含裂纹平板，在外力作用下发生 局部屈服后，其净截面的内力应当与外界平衡虚线表示发生塑性变形前，v -0的平面内法向应力的分布规律.打二K（图中虚线所示）J2兀r此曲线下的面积为R = c.-y（x）dx =外 力应力松弛后：F2 二 Cy dx=外力屈服区内的最大应力称为有效屈服应力 二ys,二ys= 、2 匚$ （平面应变）Bs （平面应力）rys为二yl7

23、=：；ys 时的 r 值,rys ）2- 2 -y.二 ：y（X）dX= ；y dX又BD与CE下的面积应相等.=FB下的面积与ABC下的面积相等.即：又“*（竺）2沙（平面应力）ys1 K 2 KI 2-R= （） =20二；（一） 兀S 8 J=在平面应力条件下，考虑应力松弛，X轴的屈服区扩大1倍.KiCs注意:上述分析没有考虑材料强化。材料强化裂纹尖端塑性区的尺寸变小 对于设计是偏于安全的、有效裂纹尺寸（讨论塑性区尺寸对应力强度因子的影响） 理论:线弹性理论.修正:有效裂纹尺寸.基本原理：设想裂纹的计算边界由o向右移到0（ 00、ry）以便使弹性区域内（即x R的区域）按线弹性理论所获得

24、的应力甘和实际应力曲线、二y基本 符合=有效裂纹尺寸a有效=a ry根据上述基本原理有：Qy也,z仝心22 = ys心2(R-ry)平面应力：1 Ki 2 1 Ki 2R ( )fyss= ry ()兀S 2兀%平面应变：屮 1 KI 2 尸 1 KI 2R(),=226= ry ()2 S 4/2兀 SR=ry2裂纹的计算边界正好在塑性区的中心、应力强度因子的计算用a有效代替a,进行K的计算1. K】表达式简单的可用解析式a.无限宽板中心穿透裂纹线弹性：Ki -匚、- a小范围屈服：K-匚 - （a ry）平面应力：厂盘巴2二 KJ 二 Mp K 其中 Mpb.深埋裂纹（椭圆片状）c.表面浅

25、裂纹二 KJ 二平面应变：ry = （-Ki）2Ki、2 . 2 心-KJ T：(a )-增大因子（塑性区修正因子）.二 ：二 21 -()a cr平面应变：ry线弹性：K小范围屈服：$二 K1 U )2CF二、二 a3 1(a 42二 JK、22-0.18( )22Ki 二1.1； a(a1 (K1)221.1 a令 Q 二 2 -0.2122 -0.212( )22-形状因子-s二 KI 二1.1 ad.表面深裂纹匚 Ki 二MiM (a 4；(?)22M1M.na二 Ki 2 T冲2 (M1M2)()24.2 -s诗2很小,令辔=0.212M1M2 r2 _0.212(二)222. K】表达式复杂一般用图解法实际有限尺寸的裂纹试样.其中：丫（a）,F（a） 般是3的复杂函数.W W W可用逐次逼近法，以a代入求K a有效 修正的Ki“当Kn-1）与Ki（n）之差满足一定的要求为止.可用图解法（参考清华大学.断裂损伤理论与应用）.