应力强度因子

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、解析函数: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 的穿透板厚的裂纹表面上, 在距离1 x a =的 范围内受均布载荷 q 作用 .利用叠加原理 :微段 集中力 qdx dK = K = 令 cos cosx a a=, cosdx a d= 111sin (1cos22(cosaaaK da-=当整个表面受均布载荷时 ,1

3、a a. 12( aK -=3. 受二向均布拉力作用的无限大平板 , 在 x 轴上有一系列长度为 2a , 间距为 2b 的裂纹 . 边界条件是周期的 : a. , y x z =.b. 在 所有 裂纹 内部 应力 为零 . 0, , 22y a x a a b x a b =- 单个裂纹时 Z =又 Z 应为 2b 的周期函数 sinzZ =采用新坐标 :z a =- sin(a Z +=当 0时 , sin,cos1222bbb=sin( sincoscossin22222a a a bbbbb+=+ cossin222a a bbb=+2222sin(cos 2cos sin(sin22

4、22222a 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. 无限大平板中的周期性的裂纹 , 且在无限远的边界上处于平板面内的纯剪切 力作用 . sin( zZ z = sin( a Z += ( K = 3

5、. 型裂纹应力强度因子的普遍表达形式 (无限大板 :( 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 利用上述结果计算在这种情况下的应力强度因子 原裂纹面11cos , sin z x =又222222221111221x z c x a

6、 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 +=+=新的裂纹面仍为椭圆 . 长轴 (1 c f c =+, 短轴 (1 a f a =+. y 向位移22002(1 2(1 (1 (1 a f a y f y E E -+=+原有裂纹

7、面 :222220( 1x z ya c y +=扩展后裂纹面 :222220( 1x z y a c y +=以 1x x =, 1z z =, 代入 原有裂纹面的边缘 y 向位移 y , 有22222011(1 (1 x z x z y y a c f a f c=-+=-+222222-=-+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 a

8、c =-又202(1 ay E -=122222 (sin cos 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 A

9、K A K W=+边中其中:A -裂纹长度; W-板宽度 当 1A W 时 22sin 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. 巴

10、里斯和薛a . 0时 接近于单边切口试样 11.12M =b . 1时 接近于半圆形的表面裂纹 11M =利用线性内插法 110.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 -=

11、取复变解析函数:( x z p iq =+, 11( z p iq =+取应力函数:2( ( ( ( z z zx z zx z =+或 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 +=+

12、 0|K iK - 取复数形式的应力强度因子 . K K iK =- 00( |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 为解析函数 -复变解析函数表达的双调和函数的普遍形式 (或复变应力函数为普遍形式 . 利用这个方法可以求解很多”无限大”平板中的穿透裂纹问题 . 二、有限宽板穿透裂纹应力强度因子的计算实际情况 :应看成有限宽计算 . 必须考

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

14、ms应力函数和应力公式Williams 应力函数 :121(1 (, cos(1 cos(1 2212j j j j 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

15、 +=+-=-+ +221(, y jjj pD A r x BW=122(1 cos(1 (1cos(3 22222j j j r 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 时都为零 . 12101

16、11lim( (21 1(1 1222r p r K D BW W -=- 1D = 应利用边界条件确定 1D , 边界条件只个边界各点的应力 , 可利用不同的边界 条件 ,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

17、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 =标准试件 , 此式为美国 SEM-E399规范2-5 确定应力强度因子的有限元法不同裂纹体在不同的开裂方式的应力强度因子是不同的 . 一些实验方法、解 析方法都有各自的局限性 , 而有限元等数值解法十分有效地求解弹塑性体的应力 和位移场 , 而应力和位移场与 K 密切相关 , 所以 , 可以通过有限元方法进行应力强 度因子的计算 .一、位移法求应力强度因子型:3(, 1cos cos 22u r k =- 3(, 1s

18、in sin 22v r k =+- 有限元法 裂纹尖端位移 (, K r = ,这种方法为外推法 二、应力法求应力强度因子型 : (, ( iy iy r f = 有限元法 (,0 y r K = K r 的关系曲线外推 K 的准确值 .应力法与位移法比较 :利用刚度法求应力时 , 应力场比位移场的精度低 (因应力是位移对坐标的偏 导数 .三、间接法求应力强度因子 (应变能释放率法 K G E =利用有限元法确定 G K . 四、 J 积分法K K:围绕裂纹尖端的闭合曲线 . T:积分边界上的力 .u:边界上的位移 .J 积分为 :uJ Wdy T ds x =-其中 12iy iy W =

19、为应变能密度 .线弹性问题 :K J G E = . 利用有限样方法计算回路积分 K . 2-6 叠加原理及其应用一、 K 的叠加原理及其应用 1. K 的叠加线弹性叠加原理 :当 n 个载荷同时作用于某一弹性体上时,载荷组在某一点 上引起的应力和位移等于单个载荷在该点引起的应力和位移分量之总和 .叠加原理适用于 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 =+ 满足叠加原理计算复杂载

20、荷下应力强度因子的方法 :将复杂载荷分解成简单载荷 , 简单载 荷可查 K 手册 .2. 实例 :铆钉孔边双耳裂纹的 K 值 叠加原理 :( ( ( ( ( ( ( 1( 2a b c d a b c K K K K K K K =+-=+ 其中 : ( ( 2b 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 = 有限板宽

21、: (a F W = ( c K = ( (2a a K = + 二、应力场叠加原理及其应用1. 应力场叠加原理 0T :无裂纹时外边界约束在裂纹所处位置产生的内应力场 .叠加原理 : ( ( ( a b c c K K K K =+= 应力场叠加原理 :在复杂的外界约束作用下 , 裂纹前端的应力强度因子等于没 有外界约束 ,但在裂纹表面上反向作用着无裂纹时外界约束在裂纹出产生的内应=+( a( b( c力 0T 所致的应力强度因子 . 如图 2. 实例 :旋转叶轮 (或轴 内孔端裂纹的 K 以等角速度 运转的叶轮 , 在内孔面有一长为 2a 的贯穿裂纹 , 求裂纹前段的 应力强度因子 .(1

22、求解无裂纹时 , 旋转体在无裂纹部位的内应力 . 有弹性力学有:22222112222223(1 8r R R r f R R r R +=+-2222211222222313(1 83R R r f R R r R +=+-+其中:f 为叶轮密度 , 为角速度 , 1R 为叶轮内径 , 2R 为叶轮外径 , r 为计算点 的位置 , 为泊松比 .= (平面应力 1=-(平面应变 0uab一般情况下 :12111050R R = 212( 1R R a 比较小 :22(1r R . 22210223(1 8R T f R r+=+(2根据类比原则 : 比较 ( d 与 ( b :内孔半径一致

23、, 裂纹大小及组态一样 , 裂纹面上下受力一致 , 外 边界无约束 , 唯一不同的是一个是有限体 , 一个是无限体 , 由于边界是自由的( d K K = (b(3.根据叠加原理 带中心孔的无限大板 , 受双向拉应力 220238f R +=时 , 孔边附近的应力 (注意无裂纹时 , 由弹性力学知 :( c( d21002(1 R T r=+( d K K = (c( 1( c a K K R = (a2.7 实际裂纹的近似处理利用断裂力学进行安全评价时 , 首先确定缺陷的大小 , 部位和形状 , 偏于安全 考虑 :夹杂、空洞、气孔、夹杂性裂纹 裂纹应针对实际问题进行分析 . 一、缺陷群的相互

24、作用 1. 垂直外应力的并列裂纹并列裂纹的作用使 K 下降 工程上偏安全考虑(1并列裂纹作为单个裂纹考虑 ;(2对于密集的缺陷群 , 假定它们在空间规则排列 , 并可把空间裂纹简化成平 面裂纹 .2. 与外应力垂直的面内共线裂纹如裂纹中心间距大于缺陷尺寸五倍以上 , 可做为单个裂纹处理 , 否则必须考 虑修正 :W M .二、裂纹形状的影响通过探伤手段 缺陷的”当量尺寸”及其部位 , 而缺陷的具体形状及实际 尺寸难以确定 裂纹形状的影响 . 1. 探伤结果是面积当缺陷的面积相同时 , 12a c =的椭圆裂纹 K 最大 以 12a c =的椭圆裂纹分析是偏于安全的 .2. 探伤的结果是最大线尺

25、寸(1当最大直径相同时 , 圆裂纹的 K 比椭圆裂纹大 以圆裂纹估算偏于安 全 .(2当缺陷长度一样时 , 贯穿裂纹 K 比其它裂纹的 K 大 以贯穿裂纹估算 偏于安全 .2.8 塑性区及其修正小范围屈服 :屈服区较小时 (远远小于裂纹尺寸 . 线弹性断裂力学仍可用 . 一、塑性区的形状和大小 1. 屈服条件的一般形式屈服条件 :材料超过弹性阶段而进入塑性阶段的条件 . a. 简单情况 :单向拉压 :12= 薄壁圆筒扭转 :s =. b. 复杂情况 :(, , , , , x y z xy xz yz f c = 用主应力表示 123(, , f c =有:最大正应力条件,最大切应力条件, v

26、on.Mises 屈服条件 (变形能条 件 ,Tresca 屈服 (切应力条件 .2. 根据屈服条件确定塑性区形状大小a. 利用米塞斯 (von.mises屈服条件 .当复杂应力状态下的形状改变能密度等于单向拉伸屈服时的形状改变能密 度 , 材料屈服 , 即 :2222122331( ( ( 2s -+-+-=对于型裂纹的应力公式: 122x y+=121sin 22= 30=(平面应力 , 薄板或厚板表面 2222cos 13sin 222s K r = -平面应力下 , 型裂纹前端屈服区域的边界方程 .当 0=时 , 201( 2sK r = 平面应变 (厚板中心 312( z =+222

27、22cos (12 3sin 222s K r *=-+ -平面应变下 , 型裂纹前端屈服区的边界方程 .当 0=时 , 210.16( (0.3 2sK r *= 221(12 ( 2sK =- b. 利用 Tresca(屈雷斯加 屈服条件 .在复杂受力下 , 当最大切应力等于材料弹性拉伸时的屈服切应力 , 材料即屈 服 .比较发现 :平面应变塑性区尺寸小 , 平面应变处于三向拉伸状态不易屈服 . 平面应变的有效屈服应力 ys 比 s 高 , 塑性区中的最大应力 1ys = 平面应变 13ys s = 考虑实际情况 3ys = 平面应力 1ys s =3.应力松弛的影响由于塑性变形引起应力松

28、弛 (应力松弛 :应变量不变 , 应力随时间降低 应力松弛 塑性区尺寸增大 , 依据 :单位厚含裂纹平板 , 在外力作用下发生 局部屈服后 , 其净截面的内力应当与外界平衡 . 虚线表示发生塑性变形前 , 0=的平面内法向应力 y 的分布规律.0|y =(图中虚线所示 此曲线下的面积为1( y F x dx =外力应力松弛后 :2y F dx *=外力屈服区内的最大应力称为有效屈服应力 ys , ( ( s ys s=平面应变 平面应力ys r 为 0|y ys =时的 r 值 , 21( 2ys ysK r = ( y y x dx dx *=又 BD 与 CE 下的面积应相等 .FB 下的

29、面积与 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 根据上述基本原理有：