1第八章疲劳裂纹扩展8.1疲劳裂纹扩展速率8.2疲劳裂纹扩展寿命预测8.3影响疲劳裂纹扩展的若干因素8.4疲劳裂纹扩展速率试验返回主目录2第八章疲劳裂纹扩展(Fatiguecrackgrowth)研究问题:含裂纹体的疲劳裂纹扩展规律,疲劳裂纹扩展寿命预测方法。研究方法裂纹尖端的应力应变场LEFM:KEPFM:d断裂力学法初始条件:初始裂纹尺寸a0?破坏条件:临界裂纹尺寸ac3Thefatiguelifeofcomponentismadeupofinitiationandpropagationstages.Thesizeofthecrackattransitionfrominitiationtopropagationisusuallyunknownandoftendependsonthepointofviewoftheanalystandthesizeofthecomponentbeinganalyzed.构件的疲劳寿命由起始和扩展二部分组成。从起始到扩展转变时的裂纹尺寸通常未知且往往取决于分析的着眼点和被分析构件的尺寸。4Forexample,foraresearcherequippedwithmicroscopicequipmentitmaybeontheorderofacrystalimperfection,dislocation(位错),ora0.1mmcrack,whiletotheinspectorinthefielditmaybethesmallestcrackthatisreadilydetectablewithnondestructiveinspectionequipment.例如,对于有显微设备的研究者,上述尺寸可能是晶粒缺陷、位错或0.1mm的量级,而对于现场检验者,则是无损检测设备可检出最小的裂纹。5问题:有缺陷怎么办?发现裂纹,能否继续使用?剩余寿命?如何控制检修?均匀、无缺陷材料循环载荷作用S-N曲线e-N曲线裂纹萌生寿命理论基础:线弹性断裂力学(1957)计算手段:计算机迅速发展;实验手段:高倍电镜、电液伺服疲劳机,电火花切割机等研究可能疲劳裂纹扩展研究需求我们已经讨论过应力寿命方法和应变寿命方法,现在讨论用于疲劳裂纹扩展估计的断裂力学法。6给定a,,da/dN;给定,a,da/dN。8.1疲劳裂纹扩展速率讨论张开型(I型)裂纹。a>>rp,LEFM力学可用。一.aN曲线二、疲劳裂纹扩展控制参量aN曲线的斜率,就是裂纹扩展速率da/dN。a(mm)a0NCCTCT123R=0K,a故K,da/dN标准试样预制疲劳裂纹恒幅疲劳实验记录a,N7裂纹只有在张开的情况下才能扩展,故控制参量K定义为:K=Kmax-KminR>0K=KmaxR<0疲劳裂纹扩展速率da/dN的控制参量是应力强度因子幅度K=f(,a),即:da/dN=(K,R,…)应力比R=Kmin/Kmax=min/max=Pmin/Pmax;与K相比,R的影响是第二位的。8三.疲劳裂纹扩展速率FCGRFatigueCrackGrowthRateR=0时的da/dN-K曲线,是基本曲线。实验a=a0R=0=consta(mm)a0NR=0iadadNaN曲线ai,(da/dN)i,ai,Kida/dN-K曲线lgda/dN10-5-610-9lg(K)~9lgda/dN12310-5-610-9lg(K)~1.da/dN-K曲线低、中、高速率三个区域:cK=(1-R)K=(1-R)KmaxthK低速率区:有下限或门槛值KthK<Kth,裂纹不扩展。高速率区:有上限Kmax=Kc,扩展快,寿命可不计。中速率区:有对数线性关系。可表达为:da/dN=C(K)mC、m和Kth,是描述疲劳裂纹扩展性能的基本参数。微解理为主微孔聚合为主条纹为主10lgda/dN12310-5-610-9lg(K)~thK微解理为主微孔聚合为主条纹为主三种破坏形式:微解理型低速率条纹型稳定扩展微孔聚合型高速率11Paris公式:da/dN=C(K)m2.裂纹扩展速率公式K是疲劳裂纹扩展的主要控制参量;疲劳裂纹扩展性能参数C、m由实验确定。3.扩展速率参数C,m的确定实验a=a0R=0记录ai、Ni(K)i=f(,ai,)(da/dN)i=(ai+1-ai)/(Ni+1-Ni)ai=(ai+1-ai)/2lg(da/dN)=lgC+mlg(K)最小二乘法C,m?128.2疲劳裂纹扩展寿命预测1.基本公式应力强度因子:afK=p中心裂纹宽板f=1;单边裂纹宽板f=1.12临界裂纹尺寸aC:有线弹性断裂判据:或CCKafK=pmaxmax2max)(1pfKaCC=疲劳裂纹扩展公式:得到裂纹扩展方程:(f,,R,a0,ac)=Ncf一般是裂纹尺寸的函数,通常需要数值积分。),,,,(),(RafRKdNda==13=CCaaNmdNafCda00)(p得到:--=--)ln()(1]11[)15.0()(1015.015.00aafCaamfCNCmmCmmCppm=2m2da/dN用Paris公式表达时的裂纹扩展方程对于无限大板,f=const.,在=const.作用下,由Paris公式da/dN=C(K)m积分有:14已知a0,ac,给定寿命Nc,估算在使用工况(R)下所允许使用的最大应力Smax。2.Paris公式的应用抗疲劳断裂设计计算已知载荷条件S,R,初始裂纹尺寸a0,估算临界裂纹尺寸ac,剩余寿命Nc.已知载荷条件S,R,给定寿命Nc,确定ac及可允许的初始裂纹尺寸a0。断裂判据:CCKafK=pmaxmax裂纹扩展方程:Nc=(f,,R,a0,ac)基本方程15解:1.边裂纹宽板K的表达式:K=1.12(pa)1/2例1:边裂纹板a0=0.5mm,载荷为max=200Mpa。R=0,材料参数ys=630MPa,u=670MPa,Kth=5.5MPa,Kc=104MPa,裂纹扩展速率为da/dN=6.9×10-12(K)3,试估算其寿命。4.临界裂纹长度ac?由断裂判据有:Kc=1.12max(pac)1/2;ac=68mm3.长度为a0的初始裂纹是否扩展?K=1.12(pa)1/2=9MPa>Kth=5.52.K=Kmax-Kmin=1.12(max-min)=1.12165.估算裂纹扩展寿命Nc:由裂纹扩展速率方程得:Nc=189500次循环a0(mm)Kc(MPam)ac(mm)Nc(千周)%0.510468189.51001.510468101.953.82.51046874.939.50.5208272198.41050.55217171.790.6讨论1:a0和Kc对疲劳裂纹扩展寿命的影响控制a0,可大大提高疲劳裂纹扩展寿命。高强脆性材料Kc低,ac、Nc小,扩展寿命可不计。17“若疲劳寿命完全由裂纹扩展所贡献,则S-N曲线可由da/dN-K关系获得且指数与Paris公式相同”。对于含有缺陷或裂纹的焊、铸件,是特别真实的。讨论2:da/dN-K曲线与S-N曲线之关系上例中,若以aL(aLaC)定义寿命,=const.,由paris公式:maWafCdNda]),,([p=积分有:Const.afCdaNLaamm==0][p此即S-N曲线:1CNSma=18讨论3:Miner理论用于裂纹扩展阶段假设尺寸为a0的裂纹,在S1、S2、S3下经n1、n2、n3循环后,扩展到aL。S1下循环n1次从a0扩展到a1;S1mn1=10)(aaadaS2下循环n2次从a1扩展到a2;S2mn2=21)(aaadaS3下循环n3次从a2扩展到aL;S3mn3=L2)(aaadaS1mN1=L0)(aaadaS2mN2=L0)(aaadaS3mN3=L0)(aaada在Si下从a0到aL的裂纹扩展寿命为N1、N2、N3。19此即Miner理论。若不计加载次序影响,Miner理论也可用于裂纹扩展阶段。在Si下循环ni次的损伤为ni/Ni,所以总损伤为:n1/N1+n2/N2+n3/N3=(++)/=110)(aaada21)(aaadaL2)(aaadaL0)(aaada若a0=0.5,aL=30mm,每年载荷谱如表。先算各Si下的裂纹扩展寿命Ni,再算ni/Ni。Si(MPa)ni(103)Ni(103)ni/Ni15030426.60.070320020180.00.11112501092.10.1086300553.30.0938设寿命为年,则有:n/N=1,=1/n/N=2.6年20例2中心裂纹宽板,作用应力max=200MPa,min=20MPa。Kc=104MPa,工作频率0.1Hz。为保证安全,每1000小时进行一次无损检验。试确定检查时所能允许的最大裂纹尺寸ai。[da/dN=4×10-14(K)4m/c]解:1.计算临界裂纹尺寸ac:对于中心裂纹宽板f=1.0,有:ac==0.086m2max)(1pcK2.检查期间的循环次数:N=0.1×3600×1000=3.6×105次213.尺寸ai的裂纹,在下一检查期内不应扩展至ac。本题m=4,由裂纹扩展方程有:--=--]11[)15.0()(115.015.00aamfCNmCmmCp注意=max-min=180Mpa,有:=160.8得到:ai=1/160.8=0.0062m=6.2mmcmciaCNa1)(1+=p讨论:若检查发现ai>6.2mm,则不安全。要继续使用,降低应力水平或缩短检查期。22如:检查时发现裂纹ai=10mm,若不改变检查周期继续使用,则应满足:)11(])([112/12/---mcmimmaaaYCNp注意,改变,临界裂纹尺寸ac不再为0.086m,而应写为:ac=解得:<159MPa,max=/(1-R)<176Mpa22max))1((1)(1pp-=ccKRK如缩短检修周期,同样可求得由ai=10mm到ac=86mm的循环次数为:N<213238次,检查期周为:T<N/(0.1×3600)=592小时。234)计算任一时刻的裂纹长度ai及其对应Ki。:ai=ao+ai,(i=1,…n)iiiaWafKp=),,(3.恒幅载荷下,裂纹扩展的数值计算方法由Paris公式有:da/dN=C(K)m=Cm(a)已知a0,,参数C、m,则数值计算方法为:1)裂纹是否会扩展?thaaKaWafK==00),,(0p2)临界裂纹尺寸ac。即:2max)(1pfKaCC=3)选取增量ai。如ai=0.01ai-1;ai越小精度越高246)假定在ai-1-ai内,da/dN不变,且:裂纹增长ai的循环数:与ai对应的累计循环次数Ni为:(/)();()/dadNCKKKKiimiii==+-12NadadNiii=/(/)iiNN=5)如(Ki-Ki-1)/Ki<e(=0.01),满足精度,继续。否则,令ai=ai/2,返回4。重复3)-6),直到ai=a0+ai=ac时,停止。由算得的(ai,Ni)数据,可作a-N曲线,且从ai扩展到ac的寿命为:Nc=Ni25Fracturemechanicsapproachesrequirethataninitiationcracksizebeknoworassumed.Forcomponentswithimperfectionsordefects(suchasweldingporous(疏松),inclusionsandcastingdefects,et