D0I:10.13374/j.issn1001053x.1994.03.021 第16卷第3期 北京科技大学学报 Vol.16 No.3 19946 Journal of University of Science and Technology Beijing Jume 1994 复合材料粘结层中具有双裂纹时的热应力 钱仁根) Shouetsu Itou2) 1)北京科技大学数力系,北京100083。 2)日本神奈川大学工学部 摘要采月拉普拉斯-傅利叶变换和施密特(Schmidt)方法,研究了复合材料因热流在裂纹周 围引起的热应力,并对用环氧树脂粘结的陶瓷-钢板进行了数值计算· 关键词陶瓷复合材料,热应力,应力强度因子/双裂纹 中图分类号TB331.0346.11 Thermal Stresses Around Two Cracks in An Adhesive Layer of Composite Materials Qian Rengen Shouetsu Itou?) I)Department of Mathematics and Mechanics.USTB.100083 Beijing,PRC 2)Kanagawa University,Japan ABSTRACT Application of the Laplace-Fourier transforms and the Schmidt method,thermal stresses around two cracks of composite materials due to heat flow are determined.As an exam- ple,the thermal stress intensity factors around two cracks in a ceramic and steel bonded com- posite are calculated numerically. KEY WORDS ceramic composites,thermal stress,stress intensity factor/two cracks 陶瓷一钢复合材料是一种新型材料,在高温环境下,必须确定它因热流干扰在裂纹周围 的应力强度因子,以便判断该材料能否正常工作.对于裂纹体的热应力问题Sh等学者都曾 研究过~,但对于带板的研究较少【56),而且从已发表的文献看,这些研究几乎都是对单裂 纹或无限个裂纹而言的, 本文研究了复合材料带板的粘结层中具有双裂纹时的热应力强度因子,并对用环氧树脂 粘结的陶瓷-钢复合材料进行了数值计算· 1 基本方程及边界条件 下面分析时忽略热耦合,并设定问题为平面应变,此时的波动方程式为: 1993-03-11收稿 第一作者男4岁散授
第 16 卷 第 3 期 1 9 9 4 年 6 月 北 京 科 技 大 学 学 报 oJ u r n a l o f U n i v e sr iyt o f S d en ce a n d T eC h n o l o g y eB ij in g V o l . 16 N o . 3 J山笼 1 99 4 复 合材 料粘结层 中具有 双裂 纹时 的热应力 钱仁根 ’ ) s ho ue st u Iot u ’ ) l) 北京 科技大学 数力 系 , 北 京 l仪幻83 2) 日本神 奈川 大学 工 学部 摘要 采 !目拉普拉斯 一 傅 利叶变换和 施密特 ( S e hi 面d t ) 方法 , 研究 了复合材 料 因 热流在 裂纹 周 围 引起的热应力 , 并对用 环氧树脂粘结的 陶瓷一 钢 板 进行了 数值计算 . 关 键词 陶瓷 复合材 料 , 热应力 , 应力强 度因 子 / 双裂 纹 中 图 分 类号 T B 3引 , 0 乡友1 . 1 hT e l , I刀a l S t res s es rA o t l n d T wo C r a c kS i n A n A d h es ive aL ye r o f C o m P o s it e M a te r i a l s Qi a n R州 g 曰 1 ’ ) hS ou e ts u sI ou ’ ) ! 、 D e P a r t m e n t o f M a t h e rn a t z c s a n d M e e h o n . e s , U S T B , 1 00 0 8 3 B e i j i n g , P R C Z ) K a n a g a wa U n ive sr l t y , J a P a n A BS T R A C T A P P lica t io n o f t h e aL p l a ce 一 F o iur er ha ns fo ns a n d t he S c l l而d t n 记 t h o d , t h e an l s t湍es a or u n d two car c k s o f co m P o s iet anr et ir a sI d ue ot h ae t fl o w a er d e t e mr i n de . sA a n ex a -m Ple , t h e t h e mar l s t n 泛洛 in t e n s iyt fa ct o sr a or im d two cat cks i n a cear imc a n d s曰 b o n d ed co -m P o s i et a er ca 1 C u l a回 n u n r ir ca l y . K E Y W O R I万 ce r a而e co m Po s i此 , t h e n 刀a l s t er S , s ter S in t e ns iyt af ct o r / t w o c r a e k s 陶瓷 一 钢 复合 材料是 一 种新型 材料 . 在高温 环境下 , 必须 确定 它 因热流 干扰在 裂纹周 围 的应力 强度 因 子 , 以便 判 断该材 料能否 正 常工 作 . 对于裂 纹体 的热应力 问题 5 1 等 学 者都 曾 研究 过 ! ’ 一 4 〕 , 但 对于带 板 的研究 较少 「’ , 6 】 , 而且从 已发表 的文献 看 , 这 些研究 几乎都是对单裂 纹 或无 限个 裂纹 而言 的 . 本 文研究 了复 合材 料带 板 的粘结 层 中具有 双裂 纹时 的热应力 强度 因子 , 并 对用环 氧树 脂 粘结 的陶瓷 一 钢复合材料 进行 了数值 计算 . 1 基本方程及边界条件 下 面 分 析时 忽略热 藕合 , 并 设定 问题 为平面应 变 . 此 时 的波动方 程式 为 : 19 3 一 03 一 1 收 稿 第一 作者 男 弘 岁 教授 DOI: 10. 13374 /j . issn1001 -053x. 1994. 03. 021
Vol.16No.3 钱仁根等:复合材料粘结层中具有双裂纹时的热应力 299. (2/ax2+a2/ay2-1/k0/at)(a2/0x2+2/ay2-1/ci.a2/0t2)Φ:=0 (2/0x2+0210y2-1/c,a2102)业,=0 (I) 勇、业:与相应的温度、位移和应力的关系为: T=1/m,(02/0x2+a2/0y2-1/c,02/0t2)Φ, (2) 4,=Φ,/0x-dψ,/0y,u,=Φ,/y+0业:/x (3) 0y=-2μ02D,/0x2+p,020,/t2+2μ,a2ψ:/0xdy) tm=2μ:020,/(0x0y)+4,(d2Φ,/0x2-02ψ:/0y) (4) 其中, k,=k:/(cP),m=x,(1+y)/(1-) (5) c=ec2,c=2(1-v)/(1-2v) (6) 上述式中,P:密度,℃:比热,k:热传动率,V:泊松比,a:线膨胀系数,4:剪切 弹性模量,c:压缩波速度,c24:剪切波速度.下标i=①,②,③,④,分别表示在裂纹 面上、下两侧的粘结层和带板中数值(如图1所示)· y 如在y=h,+h面上急剧地加热,在 y=-(h2+h,)面上急剧地冷却,并将温度 h ③ 分别保持为T,和一T,而裂纹面处于绝 ① 热状态,则在上、下表面上应力为零,温 -b-a ab② 度分别为THc)及-TH(t)其中,H(④) h ⊙ 为Heaviside阶跃函数.在y=h,y=-h 的界面上,位移、应力、温度及热流强 度均连续.在y=0的裂纹所在面上应力 图1复合材料及坐标系 和热流强度连续,而且: Fig.1 Composite plane and coordinate system y=0,ab时, 48=u8,v8=u8,T8=T2 (8) 式中的上标“0”表示是y=0面上的数值. 2理论分析 对式(1)进行拉普拉斯一付利叶变换并求解,则得变换后的:和亚:为: =Ausinh(y)+A2;cos h(iy)+A3i sin h(3y)+A4 cosh(y3y) y-Bsinh(Y2xy)+Bcosh(72y)(=①,②,③,④). (9) 其中A,B0=1,4,k=1,2)为待定常数,Y0=1,2,3)为特征根. 经推导粘结层①中y=0面上的变换后的应力、温度和位移可表示为 6,9/(2μo)=A1ok+A2Dk,”+A3Dk+mmT,·2πδ(传)/s·L T0/2μ@)=i(Aak,°+Ank,9+Aokg+m。T。·2π份/5·Lg) (10) T1y=Amk8+A20k0+A,ok8+moT,·2δ(5/s·L (11)
V o l · 1 6 o N · 钱 仁根等 3 复合材料粘结层 中具有双裂纹:时的热应力 · 卯 2 · 于 日 x 2 / ( 日 + ’ 日夕 / ’ 一 l / k ` 。 刁 / r 刁 ) ( , /刁x ’ + 刁, /即 , 一 l / c } : 口 ’ /刁t ’ ) 。 ` = o 丈日, /日x , + 日 , /口夕 , 一 l / e ; ` 日 , /刁t , )沙 : = O ( l ) 央 、 妙 ` 与相 应 的温 度 、 位移 和应力 的关 系为: 不= 1 / m `夕 ’ /日x Z + a Z /刁夕2一 l / e ; ` 日, / a t , ) 巾 ` ( 2 ) u ` = a 。 , / 刁x 一 日沙 ` / 刁夕 , v `二 日。 ` / 日夕+ 日妙 ` / 日x ( 3 ) 。 , ` = 一 2。 `。 ’ 。 ` / 刁x ’ + 。 `日 ’ 。 ` /刁t ’ + 2。 `刁 ’ 沙 ` /夕 x 口, ) T 二 , ` = 2 拼`日, 。 ` / (日x 日, ) + 拼`夕 , 。 ` / 刁x , 一 刁, 沙 ` / 刁, , ) ( 4 ) 其 中 , k ` = k ` / ( c ; p , ) , m , = : ` ( l + v ` ) / ( l 一 v ` ) ( 5 ) c卜 。子 c 圣 ` , 。卜 2 ( l 一 v ` ) / ( l 一 Z v ` ) ( 6 ) 上 述式 中 , p ` : 密 度 , :;c 比热 , k ` : 热 传 动 率 , v ` : 泊 松 比 , : ` : 线膨 胀 系 数 , 拜` : 剪 切 弹性模量 , c l` : 压缩波速度 , 几 ` : 剪切 波速度 . 下标 i 二 ① , ② , ③ , ④ , 分别 表 示 在 裂 纹 面上 、 下两侧 的粘结 层和 带板 中数值 (如 图 1 所示 ) . 如 在 y = h l + 气面上 急剧地 加热 , 在 ③一①④② y = 一 (气十 瓜)面上 急剧地 冷却 , 并将温 度 分别保 持 为 0T 和 一 0T , 而 裂纹 面处于 绝 热状态 , 则在 上 、 下表 面上应力 为零 , 温 度分别 为 OT H ()t 及一 OT H ( )t , 其 中 , (H )t 为 H ae v is kl e 阶跃函 数 . 在 y 二 气 , y 二 一 气 的界 面上 , 位 移 、 应 力 、 温度及 热流 强 度均 连续 . 在 夕= O的裂 纹所在 面上应 力 和热 流强度 连续 , 而且 : 夕= O , a I b 时 , 。 急= 。 急 , 。 签= 。 昌 , T 急二 T 急 一 b 一 a 七卜于仁 下| 二l ; 上 图 1 复合材料及坐标系 瑰 . 1 C o m po 滋et p l a ne a n d e o r 山n a et s y s et m 口T 脚即 二 口T 急/。y = O (7 ) ( 8 ) 式 中的上标 “ 0 ” 表示 是 y 二 O 面 上 的数值 2 理论分 析 对式 (l ) 进行拉 普拉 斯一付 利 叶变换并 求解 , 则得 变换后 的 示) 和 证: 为 : 币: = A l ` s ihn (下, : 夕) + A Z` co s h (下, , 夕) + A 3 、 s in h (下3 .夕) + A 4 : co s h (下3 ,夕) 叱 = B I : s in h (下2`夕) + B : co s h ( 下2 `夕) ( i = ① , ② , ③ , ④ ) . 其 中 A , , B ; ` 仃= l , … 4 , k = l , 2 ) 为待定 常数 , , , ` 仃= l , 2 , 3 ) 为特征根 . 经推 导粘结层 ① 中 夕= O 面 上 的变换后 的应力 、 温度 和位移 可表 示 为 司急/ (2 拜 ① ) 一 A l ① k尸+ A Z① 叮 + A 3 ① 衅 + 。 ① 0T · 2 二 占({) s/ · 衅 万x’@,0 /(2拜① ) = i l(A ① 畔 + A Z① 砰 + A 3① 畔 + 。 ① 0T · 2 兀 j ({) / s · 片 ) 行 ;0/ 勺一 A , ① k霭十 A Z① k浮+ A 3 ① 蜡 十 m ① 0T · 2动 ({) / 、 · 片 ( 9 ) ( 10 ) ( 1 1)
…300 北京科技大学学报 1994年No.3 T8=A1ok9+Aok0+A0kg+m0T。·2r6(⑤/5·Lg -iu8=A1ok8+A42ok0+A,ok8+moT%·2πδ(⑤)/s·Lg9 元8=A1ok8+Aok9+Aok8+mmT,·2π8()/s·Lg (12) 式中i为虚数符号,k⑧、…、k8、L⑧、…、L8均为已知函数,S,5为积分变换的变量. 同理,在粘结层②中y=0面上,也可以得到变换后的温度和位移为: T8=A1ok9+Aok+Aokg+2πmoT,·δ(⑤)/s·L号 -i=Awk+A2ok4+Aokg+2πmoT。·δ(G)/s·Lg 元”=A1ok6+A2mk+A0k8+2πmoT。·δ(⑤)/s·L8 (13) 在拉普拉斯空间中,将裂纹面上、下的温度差、位移差展开成如下的级数: 立c./2n)sin(sin't(a+b-x/b-a》-m/2 π(x8-t0)= (ax,|x>b) e,/(2n)sin(nsin-((a+b-Ixl)/(b-a)-n π(o9-v0)= (alx,Ixi>b) ∑d,/2m)sin(nsin'《a+b-x/b-a》-na/2)sgn π(u°-ua8)= (ax,Ixl>b) 式中ce,和dn为待定系数,sgn(x)为符号函数. 由式(12)、(13)消去待定系数A1m,A2,A,心,并与经付利叶变换后的(14)至(16)式 联立,可得拉普拉斯-付利叶空间中裂纹面上的应力及温度梯度.利用边界条件式(7), 经整理后可得: d+e.()R()=-U.(x.s) 喜d0hx.+gsG,-tc,) (17) 2c.swx,=-mDT15…R,0 n= 式中的2,R,H,G,W,U'。和R均为由推导所得的已知函数,式(17)利用施密特方 法7.可求得待定系数cn,dn和e.· 在拉普拉斯空间中,裂纹尖端a,b处的σ,8,t”呈现奇异性,因此,经数学推导, 可得应力强度因子为 )=2aa---豆e,2:051W2mb-0
· 0 3 · 北 京 科 技 大 学 学 报 望砰 年 1 N O . 3 雳二 A T ,① k 甲十 A Z① k 岁十 A J ① k 罗十 。 ① 0T · 2 二 占(匀/ 、 · L曾 一 俪 罄二 A I① k 乳十 A Z① k 霭+ A , ① k 界+ m ① 0T · 2 二 占(匀/ 、 · L 罗 石爹一 禹 ① k 霍十 凡 ① 灸含+ 凡 ① k 霍 + m ① 几 · 2 二 创灼s/ · 片 l( 2) 式 中 i 为虚 数符 号 , k 甲 、 … 、 k霭 、 L 甲 、 … 、 L 胃均为 已 知 函数 . 5 , 亡为积 分变 换 的变量 . 同理 , 在 粘 结层 ② 中 y 二 O 面上 , 也 可 以得到 变换后 的温度 和位移 为 : T 爹一 A . ① k 号+ 儿 ① k 罗+ A , ① k 罗十 2 7T m ① 0T · 占() / 、 · 研 一 俪 彗一 A I① k 散+ A Z① 嘴 十 A ; ① 蝶 +2 7T m ① 0T · 占({) / 、 · 蜡 。 罗= 街 ① k 篡十 儿 ① 畔 + 儿 ① 蝶 +2 二 ① 0T · 创{)/ 、 · 碍 (l3 ) 在拉 普拉 斯空 间 中 , 将 裂纹 面上 、 下 的温度 差 、 位 移差展 开成如下 的级数 : 艺 。 。 / (2 n ) s m ( 。 s血 一 ’ (( 。 + 吞一 l x l) / (b 一 a )) 一 。二 / 2 ) 叮 : 梦一 : 望 ( a } x } , ! x } > b ) ( 14 ) 艺 e 。 / (Z n ) s in ( n s i n 一 ’ (( a + 石一 l x l) / (乙一 a )) 一 n 二 / 2 ) 武 。 梦一 。 罗 a( l {xl , }xl > b) ( 1 5 ) 2了| 、、| J ! / … 、 一一一 艺J 。 / (Z n ) s 谊 ( n s in 一 ’ ( ( a + b 一 l x }) / (b 一 a )) 一 n 二 / 2 ) s , ( x ) a( I }xl , }xl > b) ( 16 ) 厂 | /、 .1 1 一 、,. ②.0 .0 一 U ① U 了.、 兀 式 中 味 、 气 和 氏为待 定 系数 , s gn (x ) 为符 号 函数 . 由式 ( 12 ) 、 (l 3) 消 去待 定系 数 A l ① , A Z ① , A 3 ① , 并 与经付 利叶 变换后 的 (l 4) 至 ( 16 ) 式 联立 , 可 得拉 普 拉斯 一 付利 叶空 间 中裂 纹 面上 的应 力 及 温 度梯 度 . 利 用 边 界 条 件 式 (7) , 经整 理后 可得 : 艺d 。 (S ) Q ( x , s ) + 艺 e 。 ( s )尺 ( x , s )一 u o ( 二 , 、 ) 艺峨( s ) 万 ( x , s ) + 艺 e 。 ( s ) G ( x , s )一 F 。 ( x , 、 ) 艺 c 。 ( s ) 体 ( x , s ) = 一 。 ① 0T / s · 尺 3 (0 ) ( 1 7 ) 式 中 的 Q , R , H , G , w , 0U , V0 和 R 3 均 为 由推 导 所 得 的 已 知 函 数 . 式 ( 17 ) 利 用 施 密 特 方 法 7[ , 8 ] 可 求得待 定系 数 吼 , 峨和 即 在拉普 拉斯 空 间中 , 裂 纹 尖 端 a , b 处 的 。 , .签 , : ;x0 ① 呈 现 奇 异 性 , 因此 , 经 数 学 推 导 , 可得应 力强 度 因子为 嵘 (s) 一 “ 忍丫“ 二 a( 一 二 ) 卜一互 “ , 2 ; ① 以 /寸2 二 (b 一 。 )
Vol.16 No.3 钱仁根等:复合材料粘结层中具有双裂纹时的热应力 301· K(o)=r品V2xa-习=424-Q)/W2xtb-a K6(s)=,8√2π(x-b) --宫e.(-ir24e-Q52✉6-。 (18) K(s)=t8√2π(x-b) --2d(-妒2e012ab-回 为求得物理空间中的应力强度因子,应将式(18)进行拉普拉斯逆变换,但目前尚不能 做到.因为 im,∫)=imf(e) 所以,利用式(19)可求得物理空间中,t→∞时的热应力强度因子, 3数值计算 首先,为了确认本分析的可靠性,将上、下带板及粘结层的材料常数设为一样,则当带 板的高度无限增大时,其计算结果应逐渐接近于同材料的无限板的解【乳.计算时设定相对 几何尺寸为a/(b-a)/2)=0.5,b/(b-a)/2)=2.5,(h,+h)=(h+h).计算结果表示在图2 中,由图可看出,当高度逐渐增大时,计算结果接近于用虚线表示的”解, 其次,将带板③,④分别设为陶瓷和钢,而①,②为环氧树脂.其材料带数表示在表1 中.设h,=h,h=h4,(h1+h)/(b-a)/2)=1.0,(h2+h,)/(b-a)/2)=1.0.计算所得的应力强 度因子Ka,K6随h,和a的变化而变化,其变化规律如图3、4所示.而K的数值极小,在 此不予表示, 表1材料常数 Table.1 Material properties 0.6 常 数 钢 陶瓷环氧树脂 4,/(GNm2) 79.2 120.0 0.893 : 0.3 0.27 0.35 g,×10-5/℃ 1.14 0.29 5.5 0.4 p,×10/kgm) 7.7 3.15 125 k/Wm'℃-) 48.6 15.5 0.188 K c×103/Nmkg1℃-)0.448 0.7- 0.712 1.05 0.0 4.0 8.0 12.0 由图3.图4可看出,对于a/(b-a)/2)=0.5, h1{6-a)/2} 1.0,1.5时的Kza,Kb的绝对值变化极小,几乎 可用同一曲线表示.说明当两个裂纹接近时,应 图2确认计算时的K。和飞 力强度因子的数值变化不大.但当两个裂纹接Fg.2 The stress intensity factors K.。andK,b 近到a/(b-a)/2)=0.1左右时,K的数值几乎 in verification of the usefulness
Vb l . 16 N b . 3 钱仁根等 : 复合材料粘结层中具有双裂纹时 的热应力 嵘 ( s ) = T 二 :急丫2 二 ( a 一 x ) = 艺d , 2 。 。 ( 一 Q会) /了2 二 ( 七一 a ) 戍 ( s ) = 。 得丫2 二 ( x 一 b ) X ~目 , 口 x ~ b + = 艺 。 , (一 l ) · + ’ 2 。 。 一 Q全/寸2 二 (卜 。 ) ( 18 ) 、 ( £ )一; ~ …一客 “ · `一 ` , · 2“ · “ “ /“ 2 · `卜 · , 为求得物理 空 间中 的应力强 度 因子 , 应 将式 ( 18) 进行 拉普拉 斯逆 变换 , 但 目前 尚不 能 做到 . 因为 夙厂(s ) 一 殃 f( ` ) 所 以 , 利用式 ( 19) 可求得 物理 空 间中 , t 一 ① 时 的热应力 强度 因子 3 数值计算 首 先 , 为 了确认本分析 的可靠性 , 将 上 、 下 带板及 粘结层 的材 料常数 设为 一样 , 则 当带 板 的 高度无 限增大 时 , 其计算结果 应逐渐 接 近于 同材 料 的 无 限板 的解 [ ’ 〕 . 计 算时设 定 相 对 几何尺寸为 a/ ( ( b 一 a) / 2 ) = .0 5 , b /《 b 一 a) / 2 ) = 2 . 5 , h( 1 + 气) = ( h Z + 凡) . 计算结果 表示在 图 2 中 . 由图可看 出 , 当高 度逐 渐增大 时 , 计 算结果 接近于 用虚 线表示 的 ” ,解 . 其 次 , 将 带板 ③ , ④ 分 别设 为陶瓷 和钢 , 而 ① , ② 为环氧树 脂 . 其 材料 常数 表 示 在 表 1 中 . 设 h , = 气 , h 3 = h 4 , ( h , + h 3 ) /((b 一 a ) / 2 ) = 1 . 0 , (气+ h 4 ) / ( b 一 a ) / 2 ) = 1 . 0 . 计算 所得 的应 力 强 度 因 子 凡 。 , 凡 。 随 h , 和 “ 的变化 而变 化 , 其 变化规律如 图 3 、 4 所 示 . 而 凡 的 数值 极 小 , 在 此不 予表示 . J ǔ .0 O 叫盆。味尼 衰 1 材料常数 1’a b卜 . 1 M a et ir a l 少。 伴币es 常 数 钢 陶 瓷 环氧树脂 .5姗.125.l05 拜` / ( G N · m 一 , ) 79 . 2 0 . 3 “ : x 10 一 , / ℃ 1 . 14 p ` x 10 , /倪 · m 一 , ) 7 . 7 k ` /W( · -m , ℃ 一 ’ ) 铭 . 6 c ` xl o ’ /伽m · kg 一 ’ ℃ 一` ) .0 48 1 20 . 0 0 . 27 .0 29 3 . 15 1 5 . 5 0 . 712 .0 89 3 0 3 5 … _ 井 必厂二冷、 、 、 ` 一 以 l 户 ( 】 由图 3 . 图 4可看 出 , 对于 a/ (( b 一 a) / 2 ) = .0 5, 1 . 0 , 1 . 5 时的 凡 。 , 凡 。 的绝 对值变 化极 小 , 几乎 可用 同一曲线表示 . 说明 当两 个裂纹 接近 时 , 应 力 强 度 因 子 的数值 变化 不大 . 但 当两个 裂纹接 近到 a/ 《 b 一 a) / 2) = 0 . 1 左 右 时 , K bZ 的数值 几乎 4 . 0 8 . 0 人/ {(b 一 a ) / 2 } ǎ泛兮! àǎ返q)ùN、ǎ . 9 · 匕六 圈 2 确认计算时的 凡 。 和 凡 。 瑰 . 2 hT e s tr es i n et 留i yt af e t o r s 凡 。 a n d K Z i n v e ir 五ca d o n o f ht e u s e fu in 砚粥
,302 北京科技大学学报 1994年No.3 不变,但K,的绝对值却急剧减小· -0.10r 0.30 ar(-a/2=010.51015 a(b-a/2j=0.1051.01 wn(/(-(41)/ -0.15 0.25 (2/e -0.20 0.20 -0.25 0,15 0.1 0.2 0.4 0.6 h,(b-a/2) -0.30 0.2 0.4 0.6 h,/b-a)/2} 图3K的数值 图4K的数值 Fig.3 The stress intensity factor K:. Fig.4 The stress intensity factor K2 4结论 从计算结果看,本文所介绍的方法是可行的,所得的数据和一些结论,对于处于热环境 中的陶瓷一钢复合材料的设计和使用是有意义的, 参考文献 1 Sih G C.On the Singular Character of Thermal Stresses Near a Crack tip.Trans ASME J Appl Mech,1962,29:587-589 2 Murakami Y.Stress Intensity Factors Handbook.Oxford:Pergamon Press,1987.1015~1105 3 Konishi Y.Atsumi A.The linear Thermoelastic Problem of Uniform Heat Flow Disturbed by a Two-dimensional Crack in a Strip.Int J Engng Sci,1973,11:1~7 4 Shindo Y,Atsumi A.Thermal Stresses in a Laminate Composite with Infinite Row of Parallel Cracks Normal to the Interfaces.Int J Engng Sci,1975,13 25~42 5 Itou S.Thermal Stresses Around An Insulated Crack in An Infinite Elastic Layer (in Japanese).Trans JSME.1991.57:1752~1758 6 Itou S,Imasato J.Thermal Stresses Around A Crack in An Adhesive Layer Between Two Dissimilar Elastic Layers (in Japanese).In:Proceedings of Conference on Materials and Mechanics(JSME),1991, A:525~527
3 0 2 北 京 科 技 大 学 学 报 1望辫 年 N o . 3 不变 , 但 K , 的绝 对值却 急剧 减小 . a 厂( 1 (b 一 a ) / 2…) }留 0 . 王 0 . 5】 { }{习 j { 飞月 · 9 9 声 ót ǎ闪ǎ/已。l à魂凶à沪(// 叼 入」 J z ( (b 一 a l/ : ) , 0 1…) ’ { “ !} , } { / } 洲 / 指屠 · 嵘昙ǎ闪健\色 。l 魂叫一(/) h (b( 一 a ) / 2 ) h . /【b( 一 a ) / 2 } 图 3 人 2 。 的数值 瑰 . 3 T h e s tr e s i n t e isn yt af e t o r K Z 。 图 4 K乃 的数值 瑰 . 4 hT e s tr e s i n t e sn i yt af e t o r K Z。 4 结 论 从计算结果 看 , 本 文所介 绍 的方法 是可行 的 , 所 得 的数 据和一 些结论 , 对于 处于热环 境 中的陶瓷 一 钢 复合材 料 的设计和使 用是有 意义 的 . 参 考 文 献 1 5山 G C . On ht e S in gu 址 C 山拍日沈r o f T卜e m 词 S位溺巴 N 比r a C ar e k it P . T ar ns AS M E J AP P I M ce h , 19 6 2 , 2 9 : 5 8 7 一 5 89 2 M l l用 ka n l j Y . S t比溺 I n et 比iyt F a Ct o sr H a nd 饮力 k . O xfo ld : eP 耳笋r 幻n P n 泛返 , 1987 . 10 15 一 1105 3 K o 比ih y , A ts 山 1刀 A . 卫l e Ul l ae r Th e n T比( 姚 ict P or b』。 n of U n jlb rm H ae t Fl o w D is 奴ir 饮劝 by a T叨 一山n r l 拐io anl C n ` k in a Sitr P . 玩t J E n gl 堪 劝 , 197 3 , 1 1 : l 一 7 4 5恤耐 O Y , A tS ~ A . 们比 m 司 S t l欲治图 in a U 口ina et 肋m p 粥 iet iw ht nI if n jet R o w of R川扭e 1 C 找比 ks No n 几al ot het ln 面血叱 . ntI J E列歹唱 Sd , 19 75 , 13 : 25 一 42 5 I ot u 5 . Th e m 司 S tr 已各留 A八〕四 d nA 加血edt C ar c k in nA 】nf 面 et E】as ict aL 界r ( in aPJ a n 已祀 ) . T alr 招 SJ M E , l卯1 , 57 : 1752 ~ 1758 6 I ot u S , 1 1n 巴祖ot J . 了比 m 司 S t l 巴妙路 A r 以川d A C ar c k in nA A id 鹉i记 加界r eB t节尼君n T认心 D i哭 im i】a r E姚宝 助碑” ( in J a P ~ ) . nI : P or eC 司in 邵 of Cb n fe 代泊ce on M aet ir a ls a nd M eC 抽而“ ( SJ M E ) , 1卯 1, A : 525 ~ 527