牟園偉，唐俊星，趙勇銘

(1.中國航空發(fā)動(dòng)機(jī)研究院，北京101304；2.西北工業(yè)大學(xué)，西安710072；3.中國航發(fā)湖南動(dòng)力機(jī)械研究所，湖南株洲412002)

外物損傷壓氣機(jī)葉片損傷容限分析

牟園偉1，唐俊星2，趙勇銘3

(1.中國航空發(fā)動(dòng)機(jī)研究院，北京101304；2.西北工業(yè)大學(xué)，西安710072；3.中國航發(fā)湖南動(dòng)力機(jī)械研究所，湖南株洲412002)

為研究外物損傷造成的初始裂紋對(duì)壓氣機(jī)葉片疲勞壽命的影響，開發(fā)了基于ANSYS平臺(tái)的三維平片裂紋擴(kuò)展整體參數(shù)化自動(dòng)模擬通用技術(shù)。通過應(yīng)力強(qiáng)度因子計(jì)算結(jié)果與文獻(xiàn)結(jié)果對(duì)比，證明所開發(fā)裂紋擴(kuò)展模型有較好的精度。利用該模型，研究了外物損傷初始裂紋位置、形態(tài)和方向?qū)簹鈾C(jī)葉片低周疲勞壽命的影響。研究表明，葉片后緣疲勞裂紋擴(kuò)展壽命最長；裂紋短/長軸比越大，疲勞裂紋擴(kuò)展壽命越長；當(dāng)初始裂紋面垂直于最大主應(yīng)力方向時(shí)，疲勞裂紋擴(kuò)展壽命最短。

航空發(fā)動(dòng)機(jī)；壓氣機(jī)葉片；外物損傷；三維裂紋；損傷容限；應(yīng)力強(qiáng)度因子；裂紋擴(kuò)展

1 引言

航空發(fā)動(dòng)機(jī)工作時(shí)，風(fēng)扇和壓氣機(jī)葉片常常會(huì)遭受到被高速氣流吸入的外來物體的撞擊而引起損傷，這種外物損傷常簡稱為FOD。據(jù)統(tǒng)計(jì)[1]，在提前翻修的發(fā)動(dòng)機(jī)零部件中，F(xiàn)OD葉片占有很大比重。外物損傷主要表現(xiàn)為在葉片局部區(qū)域出現(xiàn)缺口、撕裂、折斷和凹坑等形式，并產(chǎn)生應(yīng)力集中、殘余應(yīng)力及剪切帶等現(xiàn)象，甚至形成初始裂紋，顯著減小葉片疲勞強(qiáng)度。由于FOD可能對(duì)發(fā)動(dòng)機(jī)的安全性與經(jīng)濟(jì)性造成嚴(yán)重危害，有必要對(duì)此加強(qiáng)研究。

2001年P(guān)eters等[2]通過模擬高速FOD的沖擊作用，研究了FOD造成的微裂紋、應(yīng)力集中、殘余應(yīng)力、塑性損傷對(duì)早期疲勞裂紋萌生和疲勞裂紋擴(kuò)展的影響。2002年Hutchinson等[3]分析了FOD對(duì)疲勞壽命的影響，結(jié)果顯示對(duì)于比較深的凹坑，彈性應(yīng)力集中是影響疲勞壽命的最主要因素，其次是殘余應(yīng)力。Hudak等[4]發(fā)展了基于門檻逼近的小裂紋斷裂力學(xué)，用以處理在外物損傷缺口處和磨損疲勞邊緣接觸區(qū)發(fā)現(xiàn)微裂紋和微裂紋增長。2005年Chen[5]用有限元方法研究了硬物撞擊薄板邊緣，采用尺度分析對(duì)撞擊和材料變量進(jìn)行分類來表現(xiàn)最小的外物損傷，并著重關(guān)注殘余應(yīng)力或應(yīng)力集中的影響，給出了不同位置的疲勞裂紋。2007年Ding等[6]研究了FOD葉片在低周疲勞、高周疲勞和高、低周疲勞載荷共同作用下早期裂紋的萌生，用有限元計(jì)算了FOD葉片的殘余應(yīng)力，討論了殘余應(yīng)力對(duì)裂紋生長的影響。

這些研究加深了人們對(duì)FOD葉片疲勞強(qiáng)度弱化機(jī)理的認(rèn)識(shí)，但工程中尚難以確切把握FOD造成的初始裂紋對(duì)葉片剩余壽命的影響?；诖?，本文用損傷容限分析方法，研究了不同初始裂紋對(duì)壓氣機(jī)葉片低周疲勞壽命的影響規(guī)律。

2 結(jié)構(gòu)裂紋擴(kuò)展有限元模擬

2.1 含裂結(jié)構(gòu)有限元模型建模

為模擬各種復(fù)雜結(jié)構(gòu)的裂紋擴(kuò)展過程，本文構(gòu)造了獨(dú)立于被模擬對(duì)象的三維裂紋模型。該模型在實(shí)體層次上表達(dá)結(jié)構(gòu)的非連續(xù)性，并鑲嵌在被模擬對(duì)象中，而裂紋擴(kuò)展則通過不斷變換實(shí)體結(jié)構(gòu)的形狀來實(shí)現(xiàn)。如圖1所示，其中陰影部分為裂紋面。整個(gè)建模思想可簡單表述為：裂紋面包含于裂紋體中，而裂紋體又鑲嵌于非裂紋體中。

裂紋參量的計(jì)算精度與裂紋尖端網(wǎng)格質(zhì)量有較大關(guān)系，為提高計(jì)算效率，常在裂尖周向安排1/4邊中節(jié)點(diǎn)奇異元，如圖2、圖3所示。在圖2、圖3的裂紋體劃分中，1、2、3、4小子塊為裂紋尖端附近子塊，其拓?fù)浣Y(jié)構(gòu)簡單，幾何體規(guī)則，網(wǎng)格劃分較容易，且可獲得較好的網(wǎng)格質(zhì)量。以上整個(gè)操作過程步驟較多，但均可在ANSYS平臺(tái)中由程序自動(dòng)完成。含裂結(jié)構(gòu)有限元建模具體方法參見文獻(xiàn)[7]。

2.2 應(yīng)力強(qiáng)度因子計(jì)算及擴(kuò)展方向判斷

求解應(yīng)力強(qiáng)度因子的方法有位移法、J積分法、應(yīng)力法、節(jié)點(diǎn)力法等，本文采用公認(rèn)簡單又具有較高精度的位移法求解。

圖4示出了橢圓率a/c=0.5、裂紋深度與板厚比a/t=0.2時(shí)應(yīng)力強(qiáng)度因子隨角度的變化。可見，除裂紋表面附近區(qū)域外，本文結(jié)果基本介于Newman-Ra?ju與Lin X B&Smith R A[8]結(jié)果之間，與Newman-Raju解析解差別均在3%以內(nèi)，說明本文的裂紋模型及網(wǎng)格計(jì)算應(yīng)力強(qiáng)度因子具有良好的精度。

裂紋體形狀主要由裂紋前緣形狀決定，而裂紋前緣形狀常采用橢圓曲線或樣條曲線描述。為簡化建模過程，本文選用橢圓曲線描述。對(duì)于橢圓曲線，在裂紋平面內(nèi)僅需橢圓圓心(Xc，Yc)、長軸長度c、短軸長度a等少量參數(shù)即可表達(dá)，如公式(1)所示。

在裂紋擴(kuò)展計(jì)算時(shí)，僅計(jì)算長、短軸兩個(gè)方向的擴(kuò)展量就可以描述出新的裂紋前緣，如圖5所示。多軸應(yīng)力狀態(tài)下，裂紋前緣呈現(xiàn)三維擴(kuò)展特征。

為計(jì)算裂紋前沿?cái)U(kuò)展角θ，如圖6所示，出現(xiàn)了一系列復(fù)合型斷裂準(zhǔn)則，如最大周向應(yīng)力準(zhǔn)則[9]，最大能量釋放率準(zhǔn)則[10]，應(yīng)變能密度準(zhǔn)則[11]及最大主應(yīng)力準(zhǔn)則[12]等。這些準(zhǔn)則最大的不同，在于是否考慮KIII與泊松比對(duì)裂紋擴(kuò)展方向的影響。如果KIII/KI較小，各準(zhǔn)則之間的差別不大，且算例數(shù)值結(jié)果說明，裂紋擴(kuò)展幾步后KIII/KI迅速減小[13]。本文應(yīng)用最大周向應(yīng)力準(zhǔn)則計(jì)算疲勞裂紋擴(kuò)展方向。裂紋角采用式(2)計(jì)算，等效應(yīng)力強(qiáng)度因子采用式(3)計(jì)算。

式中：KI、KII、KIII分別表示Ⅰ型、Ⅱ型、Ⅲ型應(yīng)力強(qiáng)度因子。

2.3 裂紋擴(kuò)展壽命預(yù)測(cè)方法

一般來說，疲勞裂紋擴(kuò)展速率與裂紋前沿應(yīng)力強(qiáng)度因子服從關(guān)系式：

對(duì)式(4)兩端積分，得到：

式中：裂紋擴(kuò)展量Δa設(shè)為定值，Δa越小模擬精度越高，但模擬循環(huán)次數(shù)也越多。

對(duì)于表面半橢圓裂紋，假定每模擬擴(kuò)展一次短軸端擴(kuò)展量為ΔaA，則長軸端擴(kuò)展量可由式(6)得到，進(jìn)而可確定新的裂紋前緣。

裂紋擴(kuò)展準(zhǔn)則采用工程中普遍認(rèn)為簡單實(shí)用的Paris公式(式(7))，該公式中與材料相關(guān)的參數(shù)C、m一般通過某一特定尺寸的CT試驗(yàn)獲得。嚴(yán)格說來，三維裂紋中裂紋前沿各點(diǎn)的應(yīng)力狀態(tài)均不一樣，也與特定尺度的CT試件不完全一樣。為考慮表面層對(duì)裂紋擴(kuò)展的延遲效應(yīng)，CA取接近平面應(yīng)變狀態(tài)下CT試件試驗(yàn)數(shù)據(jù)，[8]。

本模型已將含裂結(jié)構(gòu)幾何體生成、含裂有限元網(wǎng)格劃分、加載求解運(yùn)算、應(yīng)力強(qiáng)度因子計(jì)算及重新生成新的裂紋前緣等步驟參數(shù)化，應(yīng)用本文的程序模塊可自動(dòng)完成裂紋擴(kuò)展分析的整個(gè)過程。

3 初始裂紋位置對(duì)葉片疲勞壽命的影響

裂紋初始尺寸a0=1.5 mm，裂紋短長軸比a0/c0= 1.0。初始裂紋分別位于葉身中部、葉身前緣和葉身后緣，葉片實(shí)體模型如圖7所示。葉片材料為TC4，其密度4.44×103kg/m3，彈性模量109 GPa，泊松比0.334，斷裂韌度73.5 MN/m3/2。Paris公式參數(shù)C=5.25×10-8，m=2.85。葉片加載轉(zhuǎn)速38 000 r/min，卸載轉(zhuǎn)速0 r/min。葉身表面裂紋從1.5 mm深擴(kuò)展至3.1 mm深時(shí)，裂紋尖端應(yīng)力強(qiáng)度因子達(dá)到斷裂韌度，裂紋擴(kuò)展壽命10 065周。葉身前緣裂紋從初始1.5 mm深擴(kuò)展至11.5 mm深時(shí)，裂紋尖端應(yīng)力強(qiáng)度因子達(dá)到斷裂韌度，裂紋擴(kuò)展壽命9 310周。葉身后緣裂紋從初始1.5 mm深擴(kuò)展至11.0 mm深時(shí)，裂紋尖端應(yīng)力強(qiáng)度因子達(dá)到斷裂韌度，裂紋擴(kuò)展壽命12 850周。疲勞裂紋擴(kuò)展結(jié)果如圖8所示。

由公式(7)可知，影響裂紋擴(kuò)展速率的關(guān)鍵參數(shù)是裂紋尖端應(yīng)力強(qiáng)度因子。應(yīng)力強(qiáng)度因子越大，裂紋擴(kuò)展速率越快，裂紋擴(kuò)展壽命越短。如果對(duì)無初始裂紋葉片施加相同載荷，葉身表面初始裂紋位置應(yīng)力相對(duì)前、后緣初始裂紋位置應(yīng)力最大，如圖9所示。但葉身表面裂紋尖端的應(yīng)力集中效應(yīng)相對(duì)前、后緣裂紋較小，所以葉身表面裂紋擴(kuò)展速率介于葉身前緣裂紋與葉身后緣裂紋之間。葉身前緣初始裂紋位置應(yīng)力大于葉身后緣初始裂紋位置應(yīng)力，而前后緣裂紋形態(tài)一致即裂紋尖端應(yīng)力集中效應(yīng)相近，導(dǎo)致葉身前緣裂紋應(yīng)力強(qiáng)度因子高于后緣裂紋應(yīng)力強(qiáng)度因子，進(jìn)而前緣裂紋擴(kuò)展壽命低于后緣裂紋擴(kuò)展壽命。

4 初始裂紋形態(tài)對(duì)葉片疲勞壽命的影響

以葉身表面半橢圓裂紋為研究對(duì)象，改變初始半橢圓裂紋的短、長軸比a0/c0=0.2、0.4、0.6、0.8、1.0，裂紋初始尺寸a0=1.0 mm，裂紋最終尺寸a0=3.0 mm(接近穿透葉背)。圖10展示了不同形態(tài)初始半橢圓裂紋的擴(kuò)展形貌，表1列出了裂紋擴(kuò)展壽命隨裂紋初始形態(tài)的變化結(jié)果。由表中可知，表面裂紋初始短、長軸比越大，即表面裂紋初始深度(短軸)一致，長度(長軸)越短，初始裂紋短軸端應(yīng)力強(qiáng)度因子越小，導(dǎo)致裂紋擴(kuò)展至同一深度的擴(kuò)展壽命越長。

5 初始裂紋方向?qū)θ~片疲勞壽命的影響

以葉身表面半橢圓裂紋為研究對(duì)象，改變初始裂紋方向與水平方向夾角θ(如圖11所示)，分別為20°、15°、10°、5°、0°、-5°、-10°、-15°、-20°。裂紋擴(kuò)展過程中裂紋平面與水平面夾角不斷變化。裂紋初始尺寸a0=1.0 mm，裂紋最終尺寸a0=3.0 mm(接近穿透葉背)。葉片疲勞壽命隨初始裂紋方向變化曲線見圖12，疲勞裂紋擴(kuò)展前后方向變化見表2。計(jì)算得出θ=-15°時(shí)，葉片疲勞裂紋擴(kuò)展壽命最短。由于葉片上下具有不同的扭角，在離心力作用下裂紋所在位置最大主應(yīng)力方向與垂直方向夾角為-15°左右，裂紋擴(kuò)展過程中裂紋向最有利于擴(kuò)展的方向偏轉(zhuǎn)。所以當(dāng)θ=-15°時(shí)，初始裂紋面幾乎與該位置最大主應(yīng)力方向垂直，裂紋擴(kuò)展速度最快。

表1 裂紋擴(kuò)展壽命隨裂紋初始形態(tài)的變化Table 1 Crack propagation life varies with morphological changes of initial crack

表2 疲勞裂紋擴(kuò)展方向變化Table 2 The fatigue crack propagation direction change

6 結(jié)論

本文用損傷容限分析方法，模擬研究了外物損傷造成的初始裂紋位置、形態(tài)和方向?qū)簹鈾C(jī)葉片低周疲勞壽命的影響。研究發(fā)現(xiàn)，相同尺寸、相同短/長軸比的初始裂紋，后緣裂紋疲勞擴(kuò)展壽命最長，葉身表面裂紋疲勞擴(kuò)展壽命次之，前緣裂紋疲勞擴(kuò)展壽命最短；初始深度一致的葉身表面裂紋，裂紋短/長軸比越大，疲勞裂紋擴(kuò)展壽命越長；葉身表面裂紋擴(kuò)展過程中向最有利于擴(kuò)展的方向偏轉(zhuǎn)。當(dāng)初始裂紋面垂直于最大主應(yīng)力方向時(shí)，疲勞裂紋擴(kuò)展壽命最短。下一步將探索用損傷容限方法模擬外物損傷葉片的高周疲勞壽命，進(jìn)而為受損壓氣機(jī)葉片的剩余安全壽命評(píng)估打下基礎(chǔ)。

[1]康繼東，陳士煊，徐志懷，等.壓氣機(jī)葉片外物損傷及其維修性的研究進(jìn)展[J].燃?xì)鉁u輪試驗(yàn)與研究，1998，11 (1)：59—62.

[2]Peters J O，Ritchie R O.Foreign-object damage and high-cycle fatigue of Ti-6Al-4V[J].Materials Science and Engineering，2001，319：597—601.

[3]Hutchinson J W，Chen X.Particle impact on metal sub?strates with application to foreign object damage to aircraft engines[J].Journal of the Mechanics and Physics of Solids，2002，50：2669—2690.

[4]Hudak S J，Chan K S，Chell G G，et al.A damage toler?ance approach for predicting the threshold stresses for high cycle fatigue in the presence of supplemental damage [C]//.TMS Annual Meeting.2002：107—120.

[5]Chen X.Foreign object damage on the leading edge of a thin blade[J].Mechanics of Materials，2005，37(4)：447—457.

[6]Ding J，Hall R F，et al.Fatigue crack growth from foreign object damage under combined low and high cycle load?ing.PartⅠ:Experimental studies[J].International Journal of Fatigue，2007，29：1339—1349.

[7]唐俊星，陸山.三維裂紋整體參數(shù)化?；椒╗J].航空動(dòng)力學(xué)報(bào)，2008，23(4)：737—741.

[8]Lin X B，Smith R A.Finite element modeling of fatigue crack growth of surface cracked plates–Part III:Stress in?tensity factor and fatigue crack growth life[J].J Engng Fracture Mech，1999，63：541—556.

[9]崔振源.表面裂紋理論及其應(yīng)用[M].西安：西北工業(yè)大學(xué)出版社，1982.

[10]Nuismer R J.An energy release rate criterion for mixed mode fracture[J].Int J Frac，1975，11：245—50.

[11]Sih G C,Cha B C K.Fracture criterion for three-dimen?sional crack problems[J].J Engng Fract Mech，1974，6：699—732.

[12]Sch?llmann M，Richard H，et al.A new criterion for the prediction of crack development in multiaxially loaded structures[J].Int J Frac，2002，117：129—41.

[13]唐俊星.含夾雜粉末盤概率壽命三維仿真方法與應(yīng)用[D].西安：西北工業(yè)大學(xué)，2008.

Tolerance analysis of foreign object damage on compressor blade

MU Yuan-wei1，TANG Jun-xing2，ZHAO Yong-ming3
(1.Aero Engine Academy of China，Beijing 101304，China；2.Northwestern Polytechnical University，Xi’an 710072，China；3.AECC Hunan Prowerplant Research Institute，Zhuzhou 412002，China)

In order to study the influence of initial crack caused by foreign object damage on compressor blade fatigue life,a total parameter-based and geometry independent general method was developed for au?tomatically simulating of 3D crack planar propagation on ANSYS platform.The accuracy of this crack mod?el was demonstrated by comparison with other analytical or experimental stress intensity factor results avail?able in the literature.The influence law of initial crack location,crack size,and crack direction on compres?sor blade fatigue life was studied by the model.The results show that the fatigue crack propagation life is longer in trailing edge than in other location and is longer as the short/length axis ratio increases.When the initial crack surface was perpendicular to the direction of maximum principal stress,fatiguecrack propaga?tion life is the shortest.

aero-engine；compressor blade；foreign object damage；3D crack；FOD(foreign object damage)tolerance；stress intensity factor；crack propagation

V232.3

：A

：1672-2620(2017)01-0052-06

2016-04-18；

：2016-10-23

牟園偉(1984-)，男，河北保定人，高級(jí)工程師，博士，主要從事航空發(fā)動(dòng)機(jī)零構(gòu)件強(qiáng)度、壽命及可靠性研究。