第05讲-应力强度因子的计算ppt课件.ppt
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_05.gif)
《第05讲-应力强度因子的计算ppt课件.ppt》由会员分享,可在线阅读,更多相关《第05讲-应力强度因子的计算ppt课件.ppt(95页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、1Shanghai University断裂力学Fracture Mechanics郭战胜郭战胜办公地点:延长校区力学所办公地点:延长校区力学所317室室平时答疑:每周一:平时答疑:每周一:5-6节节晚修答疑晚修答疑:每周一:每周一:18:00-20:30地点:地点:HE108或或HE104b2应力强度因子计算应力强度因子计算 预备知识预备知识:映射与广泛柯西积分公式映射与广泛柯西积分公式)(Sz.由已知解析函数经实轴或圆弧由已知解析函数经实轴或圆弧映射(反射)映射(反射)而得新的解析函数而得新的解析函数1.实轴映射实轴映射)(zf解析 Sz, 求)(zf也解析)(Sz定义)()(zfzf S
2、z设),(),()(yxiqyxpzf定义),(),( ),(),()()(11yxiqyxpyxiqyxpzfzf Sz用用1p,1q的柯西柯西黎曼条件,易证)(zf 也解析解析柯西黎曼条件柯西黎曼条件ReImImZZZyx ImReReZZZyx 32.单位圆上的映射单位圆上的映射/1ireier1/1若),(),()(iqpf,可导出:)/1 ()1()(fff,解析解析41.内内 内映射内映射 1)(kkkcZ2.外外 内映射内映射 0)(kkckcZ例例53.外外 外映射外映射 4.内内 外映射外映射 kkkcZ1)(0/)(kkcRRZ6)(/ 在 内不为零, 上, 本身可以是奇异
3、的,本身可以是奇异的,它对应它对应 平面上的角点平面上的角点)(Z k待定kknkieAZ1)()(1950,Darwin)5.76.7.)()(22mmiaZ2)11(1ln)(HZ8二二.柯西积分公式与广泛柯西积分公式柯西积分公式与广泛柯西积分公式F(t)F(z)L闭曲线,方向逆时针S内有限域,S无限域1.内域柯西公式内域柯西公式 )(zF在S内解析,在LS上连续S 0S )()(21zzzFdtzttFiC92.外域柯西公式外域柯西公式)(zF在 内解析,(包括 )S z z0)(z FS 0 S )()(21zzzFdtzttFiC3.含极点的广泛内域柯西公式含极点的广泛内域柯西公式)
4、(zF 在 内 处为 ,有n阶极点,除此以外,在 内解析Saz S)()()()()(001zFzgzFazAzFnsss则 S S )()()()(21zzzgzFdtzttgtFiC时,则104.外域广泛柯西积分公式外域广泛柯西积分公式)(zF 在 内解析, 处, ,则在 处展成级数有Sz0)(zF解析在主部 SzhnnzazazazazaazF221),(2210/)(则 S )(S )()()(21zzhzzhzFdtzttFiC1112Muskhelisvili穆什海里什维利数学弹性力学的几个基本问题Nikoloz (Niko) MuskhelishviliNikoloz (Niko
5、) Muskhelishvili (Georgian格鲁吉亚: February 16 1891 - July 16, 1976) was a notable Georgian and Soviet mathematician, one of the founders and first President (1941-1972) of the Georgian SSR Academy of Sciences (now Georgian Academy of Sciences) (then ), Doctor of Physical and Mathematical Sciences (193
6、4), Professor (1922). He is often referred by the Russian version of his name, Nikolai Ivanovich Nikolai Ivanovich Muskhelisvili.Muskhelisvili.是搞数学弹性理论的人必读的书。中文版是依据1953年出版的俄文第四版翻译的。1977年,Springer出版社根据当时最新的俄文修订版,推出了英文本:Muskhelishvili: Some Problems of the Mathematical Theory of Elasticity。中文本自五十年代出版后
7、,再没有修订过 13In 1914 he graduated from the St. Petersburg University (Russia).In 1917-1920 Muskhelishvili was Assistant Professor of this University, in 1920-1922 - Associate Professor of the Tbilisi State University (TSU), in 1922-1976 - a Professor of TSU, in 1941-1972 - first President of the Georgi
8、an SSR Academy of Sciences,in 1972-1976 - Honorary President of GAS.In 1939 Muskhelishvili was elected as Academician (Full Member) of the Academy of Sciences of the USSR (now the Russian Academy of Science.Muskhelishvili was author of outstanding scientific works in the fields of singular integral
9、equations, mathematical physics, theory of elasticity, etc. Muskhelisvili14由应力强度因子表达的脆性断裂准则为CKK11进行断裂安全分析时 1)需要计算构件的 值由构件的尺寸、形状和所受的载荷形式决定; 1K2)测定材料的 。 CK1用实验测定材料的 时,必须首先确定试件的标定式。 CK1因此,计算各种构件的应力强度因子应力强度因子,是线弹性断裂力学的一项重要任务。15计算 值的几种方法 K1.解析法:复变函数法、积分变换;2.数值解法:边界配置法、有限元法;3.实验标定法:柔度标定法;4.实验应力分析法:光弹性法.解析
10、法只能计算简单问题,大多数问题需要采用数值解法。工程中广泛采用有限元法,而且随着计算机技术的发展,能够计算越来越复杂的问题。其它求应力强度因子的方法,及工程估算和实验方法可查阅有关文献。o 对于一般一般的二维裂纹问题,可以用KolosovMuakhelishvili的方法程序性地程序性地求解应力和位移场以及应力强度因子,但这种方法求解过程需要数学的技巧。o 对于某些特殊情况特殊情况,可以采用Westergaard函数,即由需要求解两个复变解析函数 和 简化为确定一个复变函数 ,从而使问题简化。当然,Westergaard函数方法也是在少数情况下才能得出解析解。 Z解析法164,0U x y ,
11、 zxiyzxiyxzziyzz224z z 2240Uz z 记:,则KolosovMuakhelishvili应力函数法17应力函数 是实函数。 22zz 242zzz z 12z zzzz 12z zz zzz 积分之:待定函数两两共轭。KolosovMuakhelishvili应力函数法18 Re zzz 这就是著名的古萨应力函数,其中, , 为解析函数。 z z 所以求解双调和函数 的问题,归结为求解解析函数 , 的问题,称之为复应力函数复应力函数。 z zxKolosovMuakhelishvili应力函数法19应力的复变函数表示应力的复变函数表示取应力组合: 2222242xyU
12、xyUzzz z 4Rexyz zz22222222224yxxyiiUxyx yUiUxyz 22yxxyiz zz注意到 ,作第二个应力组合:KolosovMuakhelishvili应力函数法20位移的复变函数表示位移的复变函数表示 2G uivzzzzH34 , Plane Strain3 , Plane Stress1HKolosovMuakhelishvili应力函数法其他详见教材其他详见教材58-60页页21I-II复合型裂纹复合型裂纹 002 2limIIIzzKKiKzzzKolosovMuakhelishvili应力函数法要确定应力强度因子,就需要确定一个解析函数。对于复杂
13、结要确定应力强度因子,就需要确定一个解析函数。对于复杂结构或载荷条件,通常使用复变函数的保角映射原理。将构或载荷条件,通常使用复变函数的保角映射原理。将 平面内的几何图形,通过平面内的几何图形,通过 映射到平面映射到平面 中,简单中,简单的几何图形,从而使求解过程大为简化。的几何图形,从而使求解过程大为简化。IIIKKiK则根据应力场计算公式,可以求得则根据应力场计算公式,可以求得K的表达式的表达式22iyxz)(wz i范例1教材教材58-6023例:无限大板内长2a的穿透裂纹,集中力 作用在右上表面 ,求应力强度因子iQPFbz 解:取映射函数)1(2)(awz24解析法求解解析法求解I-
14、III-II复合型裂纹的应力强度因子复合型裂纹的应力强度因子复变数: iyxziyxz取复变解析函数: ( )x zpiq11( ) zpiq取应力函数 2( )( )( )( )zzzx zzx zRe ( )( )zzx z或 满足双调和方程 24范例125分析第一应力不变量 22224Re ( )xyx zxy对于.型复合裂纹型: ReImxIIZyZReImyIIZyZ| |0| |0| |0()2Re2Re2IxyIIKZ型: 2ImRexIIIIZyZReyIIyZ 000() |2Im|2Im|2xyKZ2526、型复合裂纹在裂纹前端处的不变量 000()|2Re|2Im|22x
15、yKK012Re()|2KiK取复数形式的应力强度因子 KKiK00()|2Re()|2xyK()4Re ( )xyx Z又 0lim2 2( )Kx Z2627若采用 2 2lim( )zaZaKzax z选择 满足具体问题的应力边界条件 ( )x z1144( )( )( )( )fF ZF ZZF ZZF Z-复变解析函数表达的双调和函数的普遍形式复变解析函数表达的双调和函数的普遍形式 或复变应力函数为普遍形式或复变应力函数为普遍形式 利用这个方法可以求解很多”无限大”平板中的穿透裂纹问题.2728三种基本裂纹应力强度因子的计算三种基本裂纹应力强度因子的计算一一. .无限大板无限大板型裂
16、纹应力强度因子的计算型裂纹应力强度因子的计算 0lim2KZ计算 的基本公式 K1.在“无限大”平板中具有长度为 的穿透板厚的裂纹表面上,距离 处各作用一对集中力P P 2axb ReImxZyZReImyZyZRexyyZ 选取复变解析函数: 2222222()Pz abZzbza29边界条件:边界条件:,0 xyxyz,za 除去 处裂纹自由 表面上 zb 0,0yxy如切出 坐标系内的第一象限的薄平板,在 轴所在截面上内力总和为P P xyx以新坐标表示 22222 ()()(2 )PaabZaba 2202lim2( )()P aKZab302.在无限大平板中,具有长度为 的穿透板厚的
17、裂纹表面上,在距离 的范围内受均布载荷q作用 2a1xa 利用叠加原理 集中力 qdx222()q adKdxax2202()aq aKdxax令 22sincosxaaxacosdxad 31111sin()10cos22sin ()cosaaaaaaaKqdqa当整个表面受均布载荷时 12sin ( )aaaKqqa3.受二向均布拉力作用的无限大平板,在 轴上有一系列长度为 ,间距为 的裂纹 x2a2b单个裂纹时 22zZza32边界条件边界条件是周期的:,yxz0,22yaxaabxab 0,0yxy22sin2(sin)(sin)22zbZzabb33采用新坐标: za22sin()2
18、()(sin)(sin)22abZaabb 当 时,0sin,cos1222bbbsin()sincoscossin22222aaabbbbbcossin222aabbb2222sin()() cos2cossin(sin)2222222aaaaabbbbbbb22sin()(sin)2cossin22222aaaabbbbb340sin22cossin222abZaabbb0sin2lim22 tan21cossin222aabKZbbaabbb2tan2baaab 2tan2wbaMab取 -修正系数修正系数,大于1,表示其他裂纹存在对 的影响 K 若裂纹间距离比裂纹本身尺寸大很多( )可
19、不考虑相互作用,按单个裂纹计算. 2125ab35二二. .无限大平板无限大平板、型裂纹问题应力强度因子的计算型裂纹问题应力强度因子的计算1.型裂纹应力强度因子的普遍表达形式(无限大板): 0lim( ) 2IIKZ2.无限大平板中的周期性的裂纹,且在无限远的边界上处于平板面内的纯剪切力作用. 22sin2( )(sin)(sin)22IIzbZzzabb22sin()2( )sin()(sin)22IIabZaabb3602lim2( )tan2IIbaKZaab 3.型裂纹应力强度因子的普遍表达形式(无限大板):0lim2( )KZ4.型周期性裂纹: 2tan2IIIbaKaab 积分变换
20、法积分变换法37取应力函数取应力函数满足双调和方程满足双调和方程:富里埃变换的(富里埃变换的(n n)阶导数:)阶导数:22222,7 1xyxyyxx y 22( , )072x y 二维双调和方程的Fourier Transforms( )( )1( )2( )()( )73nni xnnnd fFedxdxFiF 38将双调和方程(将双调和方程(7-2)作傅立叶变换)作傅立叶变换其中其中方程(方程(7-4)的一般解)的一般解2222220074dFdy ( , )( , )( , )75i xyFx yx y edx( , )()()76yyyABy eCDy e二维双调和方程的Four
21、ier 变换39应用反演公式:应用反演公式:及应力变换:及应力变换:1( , )( , )772i xx yy ed22222222( , )i xi xxxi xi xi xyyi xi xi xxyxydyedxedxydyedxedxedxxedxedxiedxx yy 二维双调和方程的Fourier 变换40得得:由反演公式,得:由反演公式,得:22278xyxyddydidy 22211221179221122i xi xxxi xi xyyi xi xxyxydededdyededdedieddy 二维双调和方程的Fourier Transforms41现讨论平面应变情形下位移的解
22、现讨论平面应变情形下位移的解作反演得:作反演得:若求得若求得 ,可得,可得 , 。222(1)(1)iduEdy2221(1)7 102i xdiuedEdy323211(1)(2)7 112i xddvedEdydx ( , )y( , )v x y( , )u x y二维双调和方程的Fourier 变换42半无限弹性平面的位移解现讨论受分布压力现讨论受分布压力 的半无限弹性平面问的半无限弹性平面问题题边界条件为:边界条件为:(1)0( ,0)( )( ,0)0(2),0yxyxyxyyxq xxy ( )q x43双调和方程的应力函数的傅立叶变换的一般解为:双调和方程的应力函数的傅立叶变换
23、的一般解为:由边界条件(由边界条件(2)可知:)可知: ,所以,所以由边界条件(由边界条件(1),确定),确定A、B:( , )()()76yyyABy eCDy e0CD( , )()yyABy e22000( ) ,( )( )1( )10( )i xyyyxyyyqqq x edxAqAdiBAqdy 半无限弹性平面的位移解44半无限弹性平面的位移解代入(代入(7-6)式应力函数的傅立叶变换)式应力函数的傅立叶变换得到应力解:得到应力解:21( , )( )(1)yyqy e()()()1( )(1)21( )(1)2( )2y ixxy ixyy ixxyqdy edqdy ediyq
24、ded 45半无限弹性平面的位移解对于平面应变问题对于平面应变问题将应力函数代入(将应力函数代入(7-10)、()、(7-11)得到位移表)得到位移表达式达式()( )y ixzqded 1( )(1 2 )21( )3(1)2y i xy i xuqey dEvqey dE 46裂纹问题的对偶积分方程现讨论裂纹边界受分布压力现讨论裂纹边界受分布压力 问题问题边界条件为:边界条件为:22(1)( ) ,0 ,00 ,0(2),0,0 ,yxyxyxyq xaxa yvxa yu vxy ( )q x47裂纹问题的对偶积分方程如果压力如果压力 分布对分布对 轴是对称的,则轴是对称的,则由边界条件
25、由边界条件 得:得:由边界条件由边界条件 得:得:o 引入代换:引入代换:o 式中式中 是贝塞尔函数是贝塞尔函数12,()( ), cos()()2qFJaa( )q xy( )yq x 02( )cos()( )(0)722qx dq xxa00()yvxa0cos()( )0()723xqdxaJ48裂纹问题的对偶积分方程利用上述代换,边界条件利用上述代换,边界条件(7-22)、(7-23)写为:写为:上式为对偶积分方程,由这一对方程决定函数上式为对偶积分方程,由这一对方程决定函数 ,于是便可求得,于是便可求得 。在求出。在求出 后,便可以后,便可以得到应力场和位移场的全部解。得到应力场和
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 05 应力 强度 因子 计算 ppt 课件
![提示](https://www.taowenge.com/images/bang_tan.gif)
限制150内