应力强度因子的计算

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

2、#39;=-选取复变解析函数：22 2(Z z b =-边界条件：a. , 0x y xy z =.b. , z a <出去z b =±处裂纹为自由表面上0, 0y xy =。c. 如切出xy 坐标系内的第一象限的薄平板，在x 轴所在截面上内力总和为p 。y '以新坐标表示：Z = ( K Z = 2. 在无限大平板中, 具有长度为2a 的穿透板厚的裂纹表面上，在距离1x a =±的范围内受均布载荷q 作用. 利用叠加原理:微段集中力qdxdK = 0K = 令cos cos x a a =, cos dx a d = 111sin ( 10cos 22(co

3、s a a a K d a -= 当整个表面受均布载荷时, 1a a .12( a K -= 3. 受二向均布拉力作用的无限大平板, 在x 轴上有一系列长度为2a , 间距为2b的裂纹.边界条件是周期的: a. , y x z =.b. 在所有裂纹内部应力为零. 0, , 22y a x a a b x a b =-<<-±<<±在区间内0, 0y x y =c. 所有裂纹前端y > 单个裂纹时 Z =又Z 应为2b 的周期函数 sinzZ =采用新坐标:z a =- sin(a Z +=当0时, sin,cos1222bbb=sin( sin

4、coscossin22222a a a bbbbb+=+ cossin222a a bbb=+2222sin(cos 2cos sin(sin2222222a a a a a bbbbbb b+=+22sin(sin 2cos sin22222a a a a bbbbb+-=sina= sinaK = = 取w M = 修正系数, 大于1, 表示其他裂纹存在对K 的影响. 若裂纹间距离比裂纹本身尺寸大很多(2125a b 可不考虑相互作用, 按单个裂纹计算.二、无限大平板、型裂纹问题应力强度因子的计算 1. 型裂纹应力强度因子的普遍表达形式(无限大板: (K Z =2. 无限大平板中的周期性的

5、裂纹, 且在无限远的边界上处于平板面内的纯剪切力作用. sin( zZ z = sin( a Z += ( K =3. 型裂纹应力强度因子的普遍表达形式(无限大板:( K = 4. 周期性裂纹: K =§2-2 深埋裂纹的应力强度因子的计算 1950年，格林和斯内登分析了弹性物体的深埋的椭圆形裂纹邻域内的应力和应变，得到椭圆表面上任意点, 沿y 方向的张开位移为:1222022(1 x z y y a c=-其中:202(1 ay E -=.为第二类椭圆积分. 有 = (于仁东书 22220sin( cos a d c=+ (王铎书 1962年,Irwin 利用上述结果计算在这种情况

6、下的应力强度因子 原裂纹面11cos , sin z x =又222222221111221x z c x a z a c a c+=+= 假设:椭圆形裂纹扩展时, 其失径的增值r 与成正比.r f = (f 远小于1rf = 边缘上任一点(, p x z ''', 有:1(sin (1 sin (1 x r f f x '=+=+=+1(cos (1 z r f z '=+=+11(, , (, p x z p x z '''均在0y =的平面内. 222242222(1 c x a z f a c a c ''

7、''''+=+=新的裂纹面仍为椭圆. 长轴(1 c f c '=+, 短轴(1 a f a '=+. y 向位移22002(1 2(1 (1 (1 a f a y f y E E '-+'=+原有裂纹面:222220( 1x z ya c y +=扩展后裂纹面:222220( 1x z y a c y '''+='''以1x x '=, 1z z '=, 代入原有裂纹面的边缘y 向位移y ', 有2222211112222222011(1 (1 x z x z

8、 y y a c f a f c'=-+=-'''+2222221111112222221(12 (12 12( x z x z x z f f f a c a c a c-=-+2f =2222200022(1 2y fy f f y fy ''=+又f = 2 y '=设各边缘的法向平面为平面应变, 有:31sin sin 22v k =+- 其中34k =- 当=时24(1 v K E -= 222216(1 2I r K E -= 2221E ( 41I K y ac =-又202(1 ay E -=122222 (sin cos

9、 I a K c a c=+ 在椭圆的短轴方向上，即2=，有I ImaxK K = =危险部位 椭圆片状深埋裂纹的应力强度因子当a c =时圆片状裂纹，2=2 I K =§2-3 半椭圆表面裂纹的应力强度因子计算一、表面浅裂纹的应力强度因子 当a B （板厚）线裂纹可以忽略后自由表面对A 点应力强度的影响 欧文假设：半椭圆片状表面线裂纹I K 与深埋椭圆裂纹的I K 之比等于边裂纹平板与中心裂纹平板的I K 值之比。I I I I K K K K =表边埋中又有：1220.1sin(1 tanI I AK A K W=+边中其中：A -裂纹长度；W-板宽度 当1A W 时22sin

10、A A W W ，tan A AW W1.1I I K K 边中 1.1I I K K =表埋1. 1. 1I I K K = 埋表 椭圆片状表面裂纹A 处的I K 值二、表面深裂纹的应力强度因子深裂纹：引入前后二个自由表面使裂纹尖端的弹性约束减少裂纹容易扩展I K 增大( I IK Me K =表面（埋藏） 其中：Me 弹性修正系数，应大于1，由实验确定一般情况下12Me M M =其中：1M 前自由表面的修正系数2M 后自由表面的修正系数 关于Me 表达式两种形式的论述 1. 巴里斯和薛a 0时接近于单边切口试样11.12M =b 1时接近于半圆形的表面裂纹11M =利用线性内插法110.

11、12(1 aM c=+-利用中心穿透裂纹弹性件的厚度校正系数 1222(tan 2B a M a B=B 板厚a 裂纹深度 c 裂纹长度当a B 时21M 浅裂纹不考后自由表面的影响 2. 柯巴亚希. 沙. 莫斯2110.12(1 2a M c=+-1222(tan 2B a M a B=表面裂纹的应力强度因子（应为最深点处） ：I K =§2-4 其他问题应力强度因子的计算一、 . 型复合问题应力强度因子的计算 复变数：iy x z +=，iy x z -=取复变解析函数：( x z p iq =+，11( z p iq =+取应力函数：2( ( ( ( z z zx z zx z

12、 =+或Re( (z zx z =+满足双调和方程分析第一应力不变量：22' 224Re(x y x z x y+=+= （推导过程略）对于. 型复合裂纹型：' Re Im x I I Z y Z =-， ' Re Im y I I Z y Z =+|0|0|0( 2Re x y IIZ += 型：'2Im Re x II II Z y Z =+ ' Re y II y Z =-000( |2Im |x y Z += 、型复合裂纹在裂纹前端处的不变量.000( |x y +=+ 0|K iK - 取复数形式的应力强度因子. K K iK =- 00( |

13、x y += 又( 4Re(x y x Z '+=( K Z '= 若采用z 坐标:( z aZ a K Z '=-= 选择( x z '满足具体问题的应力边界条件.这种方法利用普遍形式函数求解应力强度因子.1144( ( ( ( f F Z Z Z Z =+ (14(, ( F Z F Z 为解析函数-复变解析函数表达的双调和函数的普遍形式(或复变应力函数为普遍形式. 利用这个方法可以求解很多”无限大”平板中的穿透裂纹问题. 二、有限宽板穿透裂纹应力强度因子的计算实际情况:应看成有限宽计算. 必须考虑的自由边界对裂纹尖端应力场和位移场的影响. 在理论上得不到完

14、全解. 通过近似的简化或数值计算方法数值解.方法:边界配置法, 有限单元法等.针对有限宽板问题:寻找一个满足双调和方程和边界条件的应力函数或复变解析应力函数.边界配置法:将应力函数用无穷级数表达, 使其满足双调和方程和边界条件, 但不是满足所有的边界条件, 而是在有限宽板的边界上, 选足够多的点, 用以确定应力函数, 然后再由这样符合边界条件的应力函数确定K 值.边界配置法:计算平面问题的单边裂纹问题, 只限于讨论直边界问题. 以三点弯曲试样为例进行说明.(1威廉氏(Williams应力函数和应力公式Williams 应力函数:121(1 (, cos(1 cos(1 2212j j j j

15、j j j r C r j +=+-=-+满足双调和方程4(, 0r =. 边界条件:裂纹上、下表面(2=±, y 和xy 均为零. 上式满足. 在边界上的边界条件的满足如下确定:在有限宽板的边界上选取足够的点, 如图, 使这一点的边界条件满足j C (1(2为了计算方便引入无量纲量:2j j j D C BW p =其中：B -试件厚度, W -试件宽度.121(1 (, cos(1 cos(1 2212j jj j jpWr j j r D j BW +=+-=-+ +221(, y jjj pD A r x BW=122(1 cos(1 (1cos(3 22222j j j r

16、 j j j j jA W -=-+- 221(, x jjj pD B r y BW =21(, xy jjj p D E r x y BW=(2K 的计算针对型裂纹 :3(1sin sin 222x =- 3(1sin sin 222y =+当0=时 . y x =(0r 00|y r K =又因为当0=时, cos 1=, 当j =1时在 乘后与r 无关, 而当 2p2,3,4j =r 有关, 当0r 时都为零. 1210111lim( (21 1(1 1222r p r K D BW W -=- 1D = 应利用边界条件确定1D , 边界条件只个边界各点的应力, 可利用不同的边界条件,

17、a. 应力.b. ,n （n 为法向）.c. n , (为切向 （3）借用无裂纹体内的边界条件求系数j D取含裂纹三点弯曲试样的左半段的受力状态和不含裂纹的悬臂梁受力是一样的. 取m 个点分析, 以2m 有限级数代替无限级数精度足够.对于不同的点有:2111my jjy j p D ABW=12111mxyjjxy j p D EBW= 其中1j E 已知, 1xy 由材料力学计算.( p a K F BW W= 1357922222( 11.6( 18.4( 87.2( 150.4( 154.8( a a a a a a F W W W W W W=-+-+其中4s W =标准试件, 此式为

18、美国SEM-E399规范§2-5 确定应力强度因子的有限元法不同裂纹体在不同的开裂方式的应力强度因子是不同的. 一些实验方法、解析方法都有各自的局限性, 而有限元等数值解法十分有效地求解弹塑性体的应力和位移场, 而应力和位移场与K 密切相关, 所以, 可以通过有限元方法进行应力强度因子的计算.一、位移法求应力强度因子型:3(, 1cos cos 22u r k =- 3(, 1sin sin 22v r k =+- 有限元法裂纹尖端位移(, K r =，这种方法为外推法 二、应力法求应力强度因子型: (, ( iy iy r f = 有限元法(,0 y r K = K r 的关系曲线

19、外推K 的准确值.应力法与位移法比较:利用刚度法求应力时, 应力场比位移场的精度低(因应力是位移对坐标的偏导数.三、间接法求应力强度因子(应变能释放率法K G E ='利用有限元法确定G K . 四、J 积分法K K:围绕裂纹尖端的闭合曲线. T:积分边界上的力.u:边界上的位移.J 积分为:uJ Wdy T ds x =-其中12iy iy W =为应变能密度.线弹性问题:K J G E ='. 利用有限样方法计算回路积分K . §2-6 叠加原理及其应用一、K 的叠加原理及其应用 1. K 的叠加线弹性叠加原理:当n 个载荷同时作用于某一弹性体上时，载荷组在某一点

20、上引起的应力和位移等于单个载荷在该点引起的应力和位移分量之总和.叠加原理适用于K 证明:00|y r K = 设在1T 载荷作用下, 有:(1(1(1000, |y y r K = 设在2T 载荷作用下, 有: (2(2(2000, |y y r K = 由叠加原理有:(1(2000|y y y =+=(1(2K K K =+ 满足叠加原理计算复杂载荷下应力强度因子的方法:将复杂载荷分解成简单载荷, 简单载荷可查K 手册.2. 实例:铆钉孔边双耳裂纹的K 值 叠加原理:( ( ( ( ( ( ( 1( 2a b c d a b c K K K K K K K =+-=+其中 : ( ( 2b

21、aK = D 为圆孔直径, 可查应力强度因子手册.板有宽度:(a F W =板宽的修正.-( a( b( c ( d 这里:2f Da a =+ 即有效裂纹长度 . ( ( 2b a K =确定( c K :无限板宽中心贯穿裂纹受集中力p 作用.K =a 为有效裂纹长度 1(2 2f a D a =+ K = 有限板宽: (a F W = ( c K = ( (2a a K = + 二、应力场叠加原理及其应用1. 应力场叠加原理 0T :无裂纹时外边界约束在裂纹所处位置产生的内应力场.叠加原理: ( ( ( a b c c K K K K =+= 应力场叠加原理:在复杂的外界约束作用下, 裂纹

22、前端的应力强度因子等于没有外界约束,但在裂纹表面上反向作用着无裂纹时外界约束在裂纹出产生的内应=+( a( b( c力0T 所致的应力强度因子. 如图 2. 实例:旋转叶轮(或轴 内孔端裂纹的K 以等角速度运转的叶轮, 在内孔面有一长为2a 的贯穿裂纹, 求裂纹前段的应力强度因子.(1求解无裂纹时, 旋转体在无裂纹部位的内应力. 有弹性力学有：22222112222223(1 8r R R r f R R r R '+=+-2222211222222313(1 83R R r f R R r R ''+=+-'+其中：f 为叶轮密度, 为角速度, 1R 为叶轮内

23、径, 2R 为叶轮外径, r 为计算点的位置, 为泊松比.'= (平面应力1'=-(平面应变0uab一般情况下:12111050R R = 212( 1R R a 比较小:22(1r R . 22210223(1 8R T f R r'+=+(2根据类比原则: 比较( d 与( b :内孔半径一致, 裂纹大小及组态一样, 裂纹面上下受力一致, 外边界无约束, 唯一不同的是一个是有限体, 一个是无限体, 由于边界是自由的( d K K =(b(3.根据叠加原理 带中心孔的无限大板, 受双向拉应力220238f R '+=时, 孔边附近的应力(注意无裂纹时, 由弹性

24、力学知:( c( d21002(1 R T r=+( d K K =(c( 1( c a K K R =(a§2.7 实际裂纹的近似处理利用断裂力学进行安全评价时, 首先确定缺陷的大小, 部位和形状, 偏于安全考虑:夹杂、空洞、气孔、夹杂性裂纹裂纹应针对实际问题进行分析. 一、缺陷群的相互作用 1. 垂直外应力的并列裂纹并列裂纹的作用使K 下降工程上偏安全考虑(1并列裂纹作为单个裂纹考虑;(2对于密集的缺陷群, 假定它们在空间规则排列, 并可把空间裂纹简化成平面裂纹.2. 与外应力垂直的面内共线裂纹如裂纹中心间距大于缺陷尺寸五倍以上, 可做为单个裂纹处理, 否则必须考虑修正:W M

25、.二、裂纹形状的影响通过探伤手段缺陷的”当量尺寸”及其部位, 而缺陷的具体形状及实际尺寸难以确定裂纹形状的影响. 1. 探伤结果是面积当缺陷的面积相同时, 12a c =的椭圆裂纹K 最大以12a c =的椭圆裂纹分析是偏于安全的.2. 探伤的结果是最大线尺寸(1当最大直径相同时, 圆裂纹的K 比椭圆裂纹大以圆裂纹估算偏于安全.(2当缺陷长度一样时, 贯穿裂纹K 比其它裂纹的K 大以贯穿裂纹估算偏于安全.§2.8 塑性区及其修正小范围屈服:屈服区较小时(远远小于裂纹尺寸. 线弹性断裂力学仍可用. 一、塑性区的形状和大小 1. 屈服条件的一般形式屈服条件:材料超过弹性阶段而进入塑性阶段

26、的条件. a. 简单情况:单向拉压:12= 薄壁圆筒扭转:s =. b. 复杂情况:(, , , , , x y z xy xz yz f c = 用主应力表示123(, , f c =有：最大正应力条件，最大切应力条件，von.Mises 屈服条件(变形能条件,Tresca 屈服(切应力条件.2. 根据屈服条件确定塑性区形状大小a. 利用米塞斯(von.mises屈服条件.当复杂应力状态下的形状改变能密度等于单向拉伸屈服时的形状改变能密度, 材料屈服, 即:2222122331( ( ( 2s -+-+-=对于型裂纹的应力公式: 122x y+=±121sin 22=±

27、30=(平面应力, 薄板或厚板表面2222cos 13sin 222s K r =± -平面应力下, 型裂纹前端屈服区域的边界方程.当0=时, 201( 2sK r = 平面应变(厚板中心 312( z =+22222cos (12 3sin 222s K r *=-+ -平面应变下, 型裂纹前端屈服区的边界方程.当0=时, 210.16( (0.3 2sK r *= 221(12 ( 2sK =- b. 利用Tresca(屈雷斯加 屈服条件.在复杂受力下, 当最大切应力等于材料弹性拉伸时的屈服切应力, 材料即屈服.比较发现:平面应变塑性区尺寸小, 平面应变处于三向拉伸状态不易屈服.

28、 平面应变的有效屈服应力ys 比s 高, 塑性区中的最大应力1ys = 平面应变13ys s = 考虑实际情况3ys = 平面应力1ys s =3应力松弛的影响由于塑性变形引起应力松弛(应力松弛:应变量不变, 应力随时间降低应力松弛塑性区尺寸增大, 依据:单位厚含裂纹平板, 在外力作用下发生局部屈服后, 其净截面的内力应当与外界平衡. 虚线表示发生塑性变形前, 0=的平面内法向应力y 的分布规律.0|y =(图中虚线所示 此曲线下的面积为1( y F x dx =外力应力松弛后:2y F dx *=外力屈服区内的最大应力称为有效屈服应力ys , ( ( s ys s=平面应变平面应力ys r

29、为0|y ys =时的r 值, 21( 2ys ysK r = ( y y x dx dx *=又BD 与CE 下的面积应相等.FB 下的面积与ABC 下的面积相等. 即:( ys ysr r ys y x dx = 又201( 2ys ysK r r =(平面应力 ys s = 2201( 2( 8s sK K R r = 在平面应力条件下, 考虑应力松弛, x 轴的屈服区扩大1倍. 平面应变条件下:ys s = 可得2 ys sK r *= 2 sK R *= 注意:上述分析没有考虑材料强化。材料强化裂纹尖端塑性区的尺寸变小, 对于设计是偏于安全的.二、有效裂纹尺寸(讨论塑性区尺寸对应力强度因子的影响理论:线弹性理论. 修正:有效裂纹尺寸.y y A R r1 r B C ys o a ry o x 基本原理:设想裂纹的计算边界由 o 向右移到 o ( oo = ry 以便使弹性区域 内(即 x > R 的区域按线弹性理论所获得的应力 y | =0 和实际应力曲线 y 基本 符合. 有效裂纹尺寸 a有效 = a + ry 根据上述基本