鄧景輝，陳平劍，付裕

中國(guó)直升機(jī)設(shè)計(jì)研究所 結(jié)構(gòu)強(qiáng)度研究室，景德鎮(zhèn) 333000

對(duì)于服役環(huán)境較為嚴(yán)酷的飛機(jī)(如艦載直升機(jī)和水上飛機(jī)等)常常存在腐蝕現(xiàn)象，比如，中國(guó)艦載直升機(jī)的主槳轂鋁制卡箍腐蝕情況比較嚴(yán)重，旋翼槳葉前緣存在銹蝕等，腐蝕損傷會(huì)嚴(yán)重影響機(jī)體材料的疲勞性能，加快疲勞裂紋的形成與擴(kuò)展，威脅飛機(jī)的飛行安全，因此，腐蝕對(duì)直升機(jī)結(jié)構(gòu)疲勞壽命的影響受到航空部門(mén)的重視和研究[1-4]。Bray等[5]在3.5%NaCl溶液中對(duì)鋁合金2024-T3和2524-T3進(jìn)行預(yù)腐蝕，然后進(jìn)行疲勞試驗(yàn)，試驗(yàn)研究發(fā)現(xiàn)：材料預(yù)腐蝕后的疲勞極限大幅度降低。May等[6]研究了腐蝕環(huán)境對(duì)高強(qiáng)度合金鋼疲勞極限的影響，研究發(fā)現(xiàn)：材料經(jīng)過(guò)預(yù)腐蝕處理后的疲勞極限比未腐蝕的疲勞極限降低了33%；預(yù)腐蝕試樣的疲勞源區(qū)存在腐蝕凹坑，且疲勞裂紋均從腐蝕凹坑處萌生。Kermanidis等[7]在NaCl溶液中對(duì)鋁合金2024-T351進(jìn)行預(yù)腐蝕，然后進(jìn)行疲勞試驗(yàn)，試驗(yàn)研究發(fā)現(xiàn)：預(yù)腐蝕后鋁合金2024-T351的脆性增加，而疲勞裂紋在腐蝕凹坑處的提前萌生是材料疲勞壽命大幅縮短的主要原因。Aydin和Savaskan[8]研究了不同腐蝕環(huán)境對(duì)鋅-鋁合金的影響，研究結(jié)果發(fā)現(xiàn)：酸性溶液和鹽水都能使鋅-鋁合金的疲勞強(qiáng)度降低，而酸性溶液對(duì)疲勞強(qiáng)度的影響更大。Wang等[9]試驗(yàn)研究了酸性、堿性和中性3種不同pH值下的溶液對(duì)鋼絲疲勞壽命的影響，試驗(yàn)結(jié)果表明：鋼絲對(duì)酸性溶液較為敏感，即材料在酸性溶液中的疲勞壽命下降最多，而在堿性和中性溶液中的疲勞壽命相差很小。

上述研究表明，腐蝕環(huán)境下的航空金屬材料疲勞性能降低，腐蝕疲勞的斷裂破壞機(jī)理表現(xiàn)為：金屬材料與腐蝕環(huán)境發(fā)生電化學(xué)腐蝕作用在材料表面產(chǎn)生腐蝕凹坑，疲勞初始裂紋從金屬表面腐蝕凹坑處萌生，在疲勞載荷的作用下裂紋進(jìn)一步發(fā)生擴(kuò)展，直至斷裂。對(duì)于金屬材料的預(yù)腐蝕疲勞，腐蝕凹坑形成階段僅限于材料與外界環(huán)境的電化學(xué)腐蝕作用，而后在疲勞載荷作用下腐蝕凹坑轉(zhuǎn)捩成疲勞裂紋，并擴(kuò)展至最終的斷裂破壞，因而，預(yù)腐蝕疲勞壽命可視為在疲勞載荷作用下，腐蝕凹坑轉(zhuǎn)捩成裂紋擴(kuò)展直至最終斷裂破壞的壽命。若將形狀不規(guī)則的腐蝕凹坑等效為初始的理想當(dāng)量裂紋，結(jié)合斷裂力學(xué)方法，并根據(jù)材料的裂紋擴(kuò)展性能，計(jì)算疲勞載荷作用下裂紋擴(kuò)展至臨界裂紋長(zhǎng)度時(shí)的疲勞壽命，就可以得到材料的預(yù)腐蝕疲勞壽命，因此，基于理想當(dāng)量裂紋的斷裂力學(xué)方法常常用于模擬航空金屬材料的預(yù)腐蝕疲勞性能[10-12]。Ghidini和Donne[13]將腐蝕疲勞過(guò)程分為7個(gè)階段，根據(jù)彈性斷裂力學(xué)方法和材料裂紋擴(kuò)展性能，預(yù)測(cè)了鋁合金材料在不同腐蝕環(huán)境中的疲勞壽命，預(yù)測(cè)結(jié)果與試驗(yàn)結(jié)果吻合良好。Codaro[14]和Ghali[15]等運(yùn)用斷裂力學(xué)方法，研究了腐蝕疲勞過(guò)程中腐蝕凹坑的大小、密度和形狀隨腐蝕時(shí)間的變化規(guī)律。Li和Akid[16]基于腐蝕凹坑擴(kuò)展模型和短裂紋理論，預(yù)測(cè)了碳素鋼材料在海水腐蝕環(huán)境中的疲勞壽命，預(yù)測(cè)結(jié)果與腐蝕疲勞試驗(yàn)結(jié)果吻合良好。Walde[17-18]和Gruenberg[19]等分析了預(yù)腐蝕2024-T3鋁合金的疲勞斷口表面腐蝕形貌特征，將腐蝕凹坑等效成表面橢圓裂紋，利用斷裂力學(xué)方法，預(yù)測(cè)了材料的疲勞壽命，預(yù)測(cè)結(jié)果相比于試驗(yàn)結(jié)果有足夠高的精度。Duquesnay等[20]研究了飛行載荷譜下7075-T6511鋁合金的預(yù)腐蝕疲勞性能，發(fā)現(xiàn)最大腐蝕凹坑深度能夠有效表征腐蝕損傷，將腐蝕區(qū)域的最大凹坑深度等效成表面橢圓裂紋，利用斷裂力學(xué)方法，能夠準(zhǔn)確預(yù)測(cè)材料在飛行載荷譜下的疲勞壽命。Newman[21]和Rokhlin[22]等將小裂紋理論和斷裂力學(xué)理論相結(jié)合，即裂紋從腐蝕凹坑處萌生階段使用小裂紋理論，初始裂紋擴(kuò)展至最后斷裂階段使用斷裂力學(xué)理論，提出了一種新的腐蝕疲勞壽命預(yù)測(cè)方法，對(duì)航空鋁合金材料的預(yù)腐蝕疲勞壽命進(jìn)行預(yù)測(cè)，發(fā)現(xiàn)這種新方法的預(yù)測(cè)結(jié)果與試驗(yàn)結(jié)果吻合得很好。

盡管?chē)?guó)內(nèi)外針對(duì)金屬材料腐蝕疲勞問(wèn)題開(kāi)展了大量研究，但是，這些工作主要集中于通過(guò)試驗(yàn)方法研究腐蝕對(duì)金屬材料疲勞和裂紋擴(kuò)展性能以及壽命的影響，事實(shí)上，由于金屬材料和腐蝕環(huán)境等的多樣性，僅僅通過(guò)試驗(yàn)，測(cè)定全部航空金屬材料在各種腐蝕環(huán)境下的疲勞和裂紋擴(kuò)展性能顯然是一項(xiàng)難以完成的艱巨任務(wù)，這將影響航空器的研制周期和飛行安全，為此，本文以直升機(jī)鋁合金材料LD2CS和LD10CS為研究對(duì)象，探討航空金屬材料預(yù)腐蝕疲勞性能的理論模擬方法，試圖以理論模擬代替實(shí)物試驗(yàn)，節(jié)省時(shí)間和費(fèi)用，為中國(guó)的直升機(jī)研制提供技術(shù)支持。

1 預(yù)腐蝕疲勞試驗(yàn)

鋁合金材料LD2CS和LD10CS主要用于直X型機(jī)主減速器前后撐桿等結(jié)構(gòu)，為了研究腐蝕損傷對(duì)其疲勞性能的影響，分別制備了鋁合金LD2CS和LD10CS的疲勞試樣，試樣的幾何形狀與尺寸如圖1所示。

根據(jù)ASTM G34-01腐蝕試驗(yàn)標(biāo)準(zhǔn)[23]，預(yù)腐蝕試驗(yàn)在ZJF-09G人工氣候腐蝕試驗(yàn)箱中進(jìn)行，全部試樣在5% NaCl溶液中按照當(dāng)量加速環(huán)境譜(如圖2所示)進(jìn)行預(yù)腐蝕周期浸潤(rùn)試驗(yàn)，當(dāng)量加速環(huán)境譜用于模擬近海海面上的實(shí)際腐蝕環(huán)境。試驗(yàn)前需要對(duì)試件進(jìn)行清洗烘干，試驗(yàn)時(shí)確保試件之間不發(fā)生相互接觸，每當(dāng)試件腐蝕到指定的腐蝕時(shí)間時(shí)應(yīng)先暫停試驗(yàn)，并從周期浸潤(rùn)試驗(yàn)箱中隨機(jī)取出試件，用蒸餾水對(duì)試件進(jìn)行清洗烘干后，再利用光學(xué)顯微鏡對(duì)試件表面進(jìn)行觀察。

圖1 鋁合金預(yù)腐蝕疲勞試樣幾何形狀與尺寸Fig.1 Geometric shape and size of pre-corroded fatigue specimen of aluminum alloys

圖2 加速腐蝕試驗(yàn)環(huán)境譜Fig.2 Environment spectra for accelerative corrosion tests

圖3給出了鋁合金材料LD10CS在不同腐蝕時(shí)間下的表觀腐蝕狀況照片。從圖3(a)可以看出，未腐蝕鋁合金試樣表面較為光滑且無(wú)明顯的缺陷和損傷，當(dāng)試樣的腐蝕時(shí)間達(dá)到980 h時(shí)，鋁合金表面出現(xiàn)小的腐蝕凹坑(如圖3(b)所示)，隨著腐蝕時(shí)間進(jìn)一步增加，腐蝕凹坑的面積、深度和密度會(huì)進(jìn)一步增加，同時(shí)，材料表面顏色變暗(如圖3(c)所示)，最后，當(dāng)腐蝕時(shí)間達(dá)到2 170 h時(shí)，材料表面出現(xiàn)暗斑且更加粗糙，表明材料開(kāi)始出現(xiàn)剝蝕，部分基體材料裸露在外(如圖3(d)所示)。

圖3 鋁合金LD10CS在不同腐蝕時(shí)間下的 腐蝕損傷照片 Fig.3 Pictures of corrosion damage with different corrosion times for LD10CS aluminium alloy

采用QBG-100kN高頻疲勞試驗(yàn)機(jī)，對(duì)鋁合金材料LD2CS和LD10CS的預(yù)腐蝕試樣進(jìn)行疲勞試驗(yàn)，測(cè)定其預(yù)腐蝕后的疲勞性能。疲勞試驗(yàn)采用拉-拉疲勞加載方式，加載波形為正弦波，加載應(yīng)力比為0.1，加載頻率為90～120 Hz。采用三級(jí)成組法，測(cè)定中、短疲勞壽命區(qū)的中值疲勞壽命，三級(jí)成組法的中值疲勞壽命分別對(duì)應(yīng)3×104～5×104、7×104～10×104、1×105～3×105次循環(huán)，每級(jí)成組法的有效試樣數(shù)不少于5件。疲勞試驗(yàn)數(shù)據(jù)如圖4所示，從圖4中可以看出，2種航空鋁合金材料LD2CS和LD10CS預(yù)腐蝕后的疲勞性能隨著腐蝕時(shí)間的增加而降低。

為了研究腐蝕環(huán)境對(duì)航空鋁合金材料疲勞性能的影響機(jī)理，利用JSM-6010LA型掃描電子顯微鏡(SEM)，對(duì)2種鋁合金材料LD2CS和LD10CS的典型疲勞斷口(如圖5所示)進(jìn)行觀測(cè)，并對(duì)試樣表面典型預(yù)腐蝕凹坑的深度a和寬度2c進(jìn)行測(cè)量(如表1所示)。從圖5可以看出，2種預(yù)腐蝕鋁合金材料的疲勞源區(qū)能觀察到明顯的腐蝕凹坑，且疲勞裂紋萌生于腐蝕凹坑，說(shuō)明在疲勞載荷作用下，因腐蝕凹坑存在應(yīng)力集中，疲勞裂紋更易在腐蝕凹坑處萌生，從而縮短了鋁合金材料的裂紋萌生壽命。

圖4 不同腐蝕時(shí)間下的疲勞試驗(yàn)結(jié)果Fig.4 Fatigue test results with different corrosion time

圖5 預(yù)腐蝕鋁合金材料的疲勞斷口形貌照片F(xiàn)ig.5 Fractographic pictures for pre-corroded aluminum alloy

表1 鋁合金材料關(guān)鍵腐蝕凹坑幾何尺寸

Table 1 Geometric size of key corrosion pit for aluminium alloy

材料試樣形式腐蝕時(shí)間t/h腐蝕凹坑深度a/μm腐蝕凹坑寬度2c/μm深寬比a/2c 980881590.553LD2CS無(wú)擦蝕15403633141.15621705434821.127 980681960.347LD10CS無(wú)擦蝕15402143680.58221703844250.904

2 腐蝕凹坑的當(dāng)量裂紋模型

如前所述，雖然腐蝕凹坑當(dāng)量裂紋的半橢圓[13,17,24]和半圓模型[25-27]得到廣泛使用，但由于材料和腐蝕環(huán)境的不確定性等因素，腐蝕凹坑的幾何形狀差異性很大，導(dǎo)致半橢圓和半圓模型與實(shí)際情況存在較大偏差。圖6(a)和圖6(b)分別給出了鋁合金材料LD2CS和LD10CS典型疲勞斷口的腐蝕凹坑幾何形狀，從圖中可以看出，腐蝕凹坑的幾何截面形狀更接近拋物線形狀(如圖6(c)所示)，其空間表面形狀更接近拋物面， 因此，本文將腐蝕凹坑表面形狀理想化為拋物面(如圖6(d)所示的實(shí)線表面形狀)，拋物面與坐標(biāo)軸x、y和z的交點(diǎn)分別為A(c，0，0)、B(0，a，0)和C(0，0，z0)。根據(jù)Green和Sneddon的假定[28]：矢徑方向裂紋擴(kuò)展量與矢徑之比為常數(shù)，由此推知，在疲勞載荷作用下拋物面形狀的腐蝕凹坑當(dāng)量裂紋擴(kuò)展之后的幾何形狀仍為拋物面(如圖6(d)所示的虛線表面形狀)，與坐標(biāo)軸x、y和z的交點(diǎn)分別為A1(c+Δc，0，0)、B1(0，a+Δa，0)和C1(0，0，z0+Δz0)，且

圖6 鋁合金材料腐蝕凹坑拋物線模型圖Fig.6 Parabolic model topography of corrosion pit for aluminum alloy

(1)

式中:f為常數(shù)，且f?1；c、a和z0分別為腐蝕凹坑當(dāng)量裂紋沿坐標(biāo)軸x、y和z方向的初始長(zhǎng)度；Δc、Δa和Δz0分別為腐蝕凹坑當(dāng)量裂紋沿坐標(biāo)軸x、y和z方向的擴(kuò)展長(zhǎng)度。

根據(jù)拋物面關(guān)于y軸的對(duì)稱(chēng)性原理，圖6(d)所示的拋物面方程可表示為

y=b1x2+d1z2+b0

(2)

式中：b0、b1和d1為拋物面的形狀參數(shù)。

由A、B、C這3點(diǎn)的坐標(biāo)，可得初始腐蝕凹坑當(dāng)量裂紋拋物面ABC的方程為

(3)

由A1、B1、C1這3點(diǎn)的坐標(biāo)，可得擴(kuò)展后的腐蝕凹坑當(dāng)量裂紋拋物面A1B1C1的方程為

(4)

圖7 腐蝕凹坑在平面z=0上投影的微面積Fig.7 Micro area of projection of corrosion pit for z=0

拋物面方程式(3)和式(4)在平面z=0的投影即拋物線AB和A1B1(如圖7所示)。又由于拋物線關(guān)于y軸對(duì)稱(chēng)，因此，只需選取拋物線AB和A1B1的1/2作為研究對(duì)象，其表達(dá)式分別為

(5)

(6)

為了計(jì)算裂紋擴(kuò)展區(qū)ABB1A1的面積(見(jiàn)圖7)，對(duì)裂紋擴(kuò)展區(qū)ABB1A1取微單元PQQ1P1，其面積dA′可表示為

dA′=hds

(7)

式中：h為微單元PQQ1P1的高；ds為弧長(zhǎng)微元。

根據(jù)式(5)，弧長(zhǎng)微元ds可表示為

(8)

令點(diǎn)P的坐標(biāo)為(x,y),P1的坐標(biāo)為(x1,y1),則x1和y1可表示為

x1=(1+f)x

(9)

y1=(1+f)y

(10)

由式(9)和式(10)，以及過(guò)點(diǎn)P的切線方程，可得點(diǎn)P1到切線L的距離為

(11)

將式(8)和式(11)代入式(7)，可得微單元PQQ1P1的面積dA′為

(12)

聯(lián)立式(4)和式(5)，并忽略f的高階項(xiàng)，可以得到

(13)

又由式(11)，可得

(14)

將式(14)代入式(13)，得到

(15)

腐蝕凹坑當(dāng)量裂紋擴(kuò)展后的體積改變量可表示為

ΔV=V2-V1=?Ω2dxdydz-?Ω1dxdydz

(16)

式中：V1和V2分別為初始和擴(kuò)展后的腐蝕凹坑當(dāng)量裂紋空間體積；Ω1和Ω2分別為V1和V2對(duì)應(yīng)的空間積分區(qū)域。

采用切片法，對(duì)式(16)進(jìn)行積分變換，可得

(17)

對(duì)式(17)進(jìn)行積分，并略去f的高階項(xiàng)，可得腐蝕凹坑當(dāng)量裂紋體積改變量為

(18)

在沿z軸方向的拉伸載荷作用下，金屬板表面腐蝕凹坑當(dāng)量裂紋為I型裂紋問(wèn)題(如圖8所示)。根據(jù)斷裂力學(xué)知識(shí)，可以得到平面應(yīng)變狀態(tài)下無(wú)限大板I型裂紋尖端沿z方向的張開(kāi)位移為[29]

(19)

式中：E為彈性模量；KI為I型裂紋應(yīng)力強(qiáng)度因子；υ為泊松比。

圖8 腐蝕凹坑擴(kuò)展斷裂力學(xué)模型Fig.8 Fracture mechanics model for corrosion pit extension

聯(lián)立式(15)和式(19)，可以得到腐蝕凹坑I型當(dāng)量裂紋的應(yīng)力強(qiáng)度因子為

(20)

平面應(yīng)變狀態(tài)下，能量釋放率G與I型裂紋應(yīng)力強(qiáng)度因子KI的關(guān)系可寫(xiě)為[29]

(21)

根據(jù)斷裂力學(xué)中能量釋放率的定義，可以得到腐蝕凹坑當(dāng)量裂紋從拋物面ABC擴(kuò)展到拋物面A1B1C1時(shí)的能量改變量ΔU為

(22)

將式(12)和式(21)代入式(22)，可得

(23)

根據(jù)斷裂力學(xué)知識(shí)和式(18)，應(yīng)力σ做的功ΔW可寫(xiě)為

(24)

由能量平衡方程，可得

ΔW=ΔU

(25)

將式(20)、式(23)和式(24)代入式(25)，并積分，得到

(26)

式中：

(27)

(28)

(29)

再將式(26)代入式(20)，可得無(wú)限大平板腐蝕凹坑的I型當(dāng)量裂紋的應(yīng)力強(qiáng)度因子表達(dá)式為

(30)

為考慮幾何形狀和裝配效應(yīng)對(duì)應(yīng)力強(qiáng)度因子的影響，工程上，常對(duì)式(30)進(jìn)行修正[30]，則有限寬板腐蝕凹坑I型當(dāng)量裂紋的應(yīng)力強(qiáng)度因子可表示為

(31)

式中：

(32)

(33)

式中：t0為板的厚度；M1、M2、M3、g1和fw均為幾何修正系數(shù)。

3 腐蝕凹坑疲勞壽命計(jì)算

金屬材料裂紋擴(kuò)展性能不僅與材料本身的斷裂韌度KC、斷裂門(mén)檻值ΔKth和厚度有關(guān)，還受到應(yīng)力比R和環(huán)境等外部因素的影響。目前，金屬材料的裂紋擴(kuò)展性能常用裂紋擴(kuò)展速率da/dN來(lái)表征，常用的裂紋擴(kuò)展速率表達(dá)式有Paris公式、Walker公式、Forman公式和四參數(shù)Forman公式等，其中四參數(shù)Forman公式同時(shí)考慮了應(yīng)力比、斷裂韌度和斷裂門(mén)檻值對(duì)裂紋擴(kuò)展速率的影響，且適用于裂紋擴(kuò)展全范圍，因此，得到日益廣泛的應(yīng)用，其表達(dá)式為[31]

(34)

式中：C、n、p和q為材料常數(shù)；ΔK為應(yīng)力強(qiáng)度因子變程;fop為裂紋張開(kāi)函數(shù)。

在沿z軸方向拉伸載荷作用下，金屬板表面的腐蝕凹坑當(dāng)量裂紋分別沿x和y軸方向擴(kuò)展(見(jiàn)圖8)，由于預(yù)腐蝕凹坑當(dāng)量裂紋沿y軸方向擴(kuò)展導(dǎo)致斷裂(見(jiàn)圖6(a)和圖6(b))，因此，本文只考慮腐蝕凹坑當(dāng)量裂紋沿y軸方向的擴(kuò)展速率。根據(jù)式(34)，腐蝕凹坑當(dāng)量裂紋沿y軸方向的擴(kuò)展速率可表示為

(35)

式中:ΔKa為腐蝕凹坑當(dāng)量裂紋沿y軸方向擴(kuò)展的應(yīng)力強(qiáng)度因子變程。

由式(31)和ΔK的定義，可得

(36)

根據(jù)式(35)和式(36)，可得疲勞載荷譜作用下第i個(gè)應(yīng)力循環(huán)引起的腐蝕凹坑當(dāng)量裂紋沿y軸方向的擴(kuò)展增量為

(37)

則經(jīng)過(guò)j個(gè)應(yīng)力循環(huán)后，腐蝕凹坑當(dāng)量裂紋沿y軸方向的總長(zhǎng)度可寫(xiě)為

(38)

根據(jù)平面應(yīng)變I型裂紋脆斷準(zhǔn)則，可知

(39)

式中：KIC為材料的平面應(yīng)變斷裂韌度；acr為腐蝕凹坑當(dāng)量裂紋沿y軸開(kāi)始失穩(wěn)擴(kuò)展時(shí)的臨界裂紋長(zhǎng)度。

根據(jù)式(35)～式(39)，可建立腐蝕凹坑處裂紋擴(kuò)展的累計(jì)求和算法流程：① 測(cè)量金屬材料表面關(guān)鍵腐蝕凹坑的幾何尺寸；② 根據(jù)構(gòu)件或試樣的幾何形狀尺寸、疲勞載荷、材料裂紋擴(kuò)展性能參數(shù)和腐蝕凹坑幾何尺寸，結(jié)合式(35)～式(38)，計(jì)算腐蝕凹坑當(dāng)量裂紋在疲勞載荷作用下應(yīng)力強(qiáng)度因子和沿y軸的擴(kuò)展長(zhǎng)度；③ 根據(jù)斷裂判據(jù)式(39)，確定腐蝕凹坑當(dāng)量裂紋沿y軸擴(kuò)展的臨界裂紋長(zhǎng)度，當(dāng)腐蝕凹坑當(dāng)量裂紋沿y軸的裂紋長(zhǎng)度等于臨界裂紋長(zhǎng)度時(shí)，即可得材料發(fā)生破壞時(shí)的疲勞壽命。

4 模型驗(yàn)證

鋁合金材料LD2CS和LD10CS的力學(xué)性能參數(shù)如表2[32]所示。根據(jù)試樣的幾何形狀尺寸(見(jiàn)圖1)、腐蝕凹坑幾何尺寸(見(jiàn)表1)、材料裂紋擴(kuò)展性能參數(shù)(見(jiàn)表2)和疲勞載荷，由式(34)～式(39)，利用裂紋擴(kuò)展的累計(jì)求和法，計(jì)算鋁合金材料LD2CS和LD10CS腐蝕凹坑當(dāng)量裂紋在疲勞載荷作用下的應(yīng)力強(qiáng)度因子和沿y軸的擴(kuò)展長(zhǎng)度，獲得腐蝕凹坑當(dāng)量裂紋沿y軸擴(kuò)展至臨界裂紋長(zhǎng)度時(shí)的疲勞壽命(如圖4所示)。采用類(lèi)似模擬方法，可分別獲得半橢圓模型和半圓模型的預(yù)腐蝕疲勞性能模擬結(jié)果(如圖4所示)。從圖4中可以看出，3種模型均能有效模擬鋁合金材料的預(yù)腐蝕疲勞性能，其中拋物線模型模擬精度最佳，半橢圓模型次之，半圓模型精度最差?；诟g疲勞斷口形貌特征而建立的拋物線模型能更精確地描述腐蝕凹坑形貌特征，計(jì)算的應(yīng)力強(qiáng)度因子更接近真實(shí)狀態(tài)，因此，拋物線模型預(yù)測(cè)疲勞性能的結(jié)果精度最高；而半圓模型僅僅使用一個(gè)參數(shù)(即腐蝕凹坑深度)描述腐蝕凹坑的形貌特征，顯然與實(shí)際情況有較大出入，預(yù)測(cè)的疲勞性能精度最差；半橢圓模型雖然也采用了2個(gè)參數(shù)(腐蝕凹坑深度和寬度)描述腐蝕凹坑的形貌特征，比半圓模型的精度高，但其數(shù)學(xué)性質(zhì)決定了所描述幾何形狀的局限性，難以精確描述腐蝕凹坑的拋物面形狀，因此，半橢圓模型預(yù)測(cè)疲勞性能的精度雖高于半圓模型，但低于拋物線模型。

表2 鋁合金材料力學(xué)性能參數(shù)[32]Table 2 Mechanical properties of aluminum alloys[32]

從圖4中還可以看出，半圓模型模擬的疲勞性能比試驗(yàn)結(jié)果保守，這是因?yàn)楦g凹坑的底部往往較為尖銳，而半圓模型描述的腐蝕凹坑底部過(guò)于平緩，因此，計(jì)算的應(yīng)力強(qiáng)度因子偏小，導(dǎo)致模擬的疲勞性能偏于保守。

5 結(jié) 論

基于鋁合金材料預(yù)腐蝕疲勞斷口的形貌特征，建立了新的鋁合金材料腐蝕凹坑當(dāng)量裂紋的拋物線模型，推導(dǎo)了腐蝕凹坑當(dāng)量裂紋的應(yīng)力強(qiáng)度因子新公式，發(fā)展了疲勞載荷作用下航空鋁合金材料的預(yù)腐蝕疲勞壽命估算方法，并模擬了預(yù)腐蝕鋁合金材料LD2CS和LD10CS的疲勞性能，得到以下結(jié)果：

1) 建立的腐蝕凹坑當(dāng)量裂紋拋物線模型、應(yīng)力強(qiáng)度因子公式和剩余壽命估算方法，能有效模擬預(yù)腐蝕鋁合金材料LD2CS和LD10CS的疲勞性能，為鋁合金材料腐蝕疲勞性能模擬提供有益參考，可大幅度減少腐蝕疲勞性能試驗(yàn)數(shù)據(jù)和成本。

2) 與傳統(tǒng)的半橢圓模型和半圓模型相比，本文提出的腐蝕凹坑當(dāng)量裂紋拋物線模型具有更高的模擬精度，實(shí)用價(jià)值更大。

[1] BUDNICK J. NAVAIR air vehicle corrosion challenges[C]∥Tri-Service Corrosion Conference. Washington, D.C.: US Department of Defence, 2003.

[2] 張有宏, 呂國(guó)志, 陳躍良. LY12CZ 鋁合金預(yù)腐蝕及疲勞損傷研究[J]. 航空學(xué)報(bào), 2005, 26(6): 779-782.

ZHANG Y H, LV G Z, CHEN Y L. Predicting fatigue life from pre-corroded LY-12 CZ aluminum test[J]. Acta Aeronautica et Astronautica Sinica, 2005, 26(6): 779-782 (in Chinese).

[3] PETROYIANNIS P V, KERMANIDIS A T, AKID R, et al. Analysis of the effects of exfoliation corrosion on the fatigue behaviour of the 2024-T351 aluminum alloy using the fatigue damage map[J]. International Journal of Fatigue, 2005, 27(7): 817-827.

[4] 譚曉明, 張丹峰, 卞貴學(xué). 腐蝕對(duì)新型高強(qiáng)度鋁合金疲勞裂紋萌生機(jī)制及擴(kuò)展行為的作用[J]. 機(jī)械工程學(xué)報(bào), 2014, 50(22): 76-83.

TAN X M, ZHANG D F, BIAN G X. Effect of corrosion damage on fatigue crack initiation mechanism and growth behavior of high strength aluminum alloy[J]. Journal of Mechanical Engineering, 2014, 50(22): 76-83 (in Chinese).

[5] BRAY G H, BUCCI R J, COLVIN E L, et al. Effect of prior corrosion on theS/Nfatigue performance of aluminum sheet alloys 2024-T3 and 2524-T3[J]. ASTM Special Technical Publication, 1997, 1298(5): 89-103.

[6] MAY M E, LUC T P, SAINTIER N, et al. Effect of corrosion on the high cycle fatigue strength of martensitic stainless steel X12CrNiMoV12-3[J]. International Journal of Fatigue, 2013, 47(2): 330-339.

[7] KERMANIDIS A T, PETROYIANNIS P V, PANTELAKIS S G. Fatigue and damage tolerance behavior of corroded 2024 T351 aircraft aluminum alloy[J]. Theory and Applied Fracture Mechanics, 2005, 43(1): 121-132.

[8] AYDIN M, SAVASKAN T. Fatigue properties of zinc-aluminum alloys in 3.5% NaCl and 1% HCl solutions[J]. International Journal of Fatigue, 2004, 26(1): 103-110.

[9] WANG S Q, ZHANG D K, CHEN K, et al. Corrosion fatigue behaviors of steel wires used in coalmine[J]. Materials and Design, 2014, 53(1): 58-64.

[10] WANG Q Y, PIDAPARTI R M, PALAKAL M J. Comparative study of corrosion-fatigue in aircraft materials[J]. American Institute of Aeronautics and Astronautics, 2001, 39(2): 325-330.

[11] ISHIHARA S, SAKA S, NAN Z Y, et al. Prediction of corrosion fatigue lives of aluminum alloy on the basis of corrosion pit growth law[J]. Fatigue Fracture & Engineering of Materials& Structure, 2006, 29(6): 472-480.

[12] SRIAMAN M R, PIDAPATI R A. Life prediction of aircraft aluminum subjected to pitting corrosion under fatigue condition[J]. Journal of Aircraft, 2009, 46(4): 1253-1259.

[13] GHIDINI T， DONNE C D. Fatigue life predictions using fracture mechanics methods[J]. Engineering Fracture Mechanics, 2009, 76(1): 134-148.

[14] CODARO E N, NAKAZATO R Z, HOROVISTIZ A L, et al. An image processing method for morphological characterization and pitting corrosion evaluation[J]. Material Science and Engineering: A, 2002, 334(1-2): 298-306.

[15] GHALI E, DIETZEL W, KAINER K U. Testing of general and localized corrosion of magnesium alloys: A critical review[J]. Journal of Materials Engineering and Performance, 2004, 13(5): 7-23.

[16] LI S X, AKID R. Corrosion fatigue life prediction of a steel shaft material in seawater[J]. Engineering Failure Analysis, 2013, 34(8): 324-334

[17] WALDE K, HILLBERRY B M. Characterization of pitting damage and prediction of remaining fatigue life[J]. International Journal of Fatigue, 2008, 30(1): 106-118.

[18] WALDE K, BROCKENBROUGH J R, CRAIG B A, et al. Multiple fatigue crack growth in pre-corroded 2024-T3 aluminum[J]. International Journal of Fatigue, 2005, 27(10-12): 1509-1518.

[19] GRUENBERG K M, CRAIG B A, HILLBERRY B M, et al. Predicting fatigue life of pre-corroded 2024-T3 aluminum from breaking load tests[J]. International Journal of Fatigue, 2004, 26(6): 615-627.

[20] DUQUESNAY D L, UNDERHILL P R, BRITT H J. Fatigue crack growth from corrosion damage in 7075-T6511 aluminum alloy under aircraft loading[J]. International Journal of Fatigue, 2003, 25(5): 371-377.

[21] NEWMAN J C, ABBOTT W. Fatigue life calculations on pristine and corroded open-hole specimens using small-crack theory[J]. International Journal of Fatigue, 2009, 31(8-9): 1246-1253.

[22] ROKHLIN S I, KIM J Y, NAGY H, et al. Effect of pitting corrosion on fatigue crack initiation and fatigue life[J]. Engineering Fracture Mechanics, 1999, 62(4-5): 425-444.

[23] American Society for Testing Materials International. Standard test method for exfoliation corrosion susceptibility in 2XXX and 7XXX series aluminum alloys: ASTM G34-01[S]. West Conshohocken, PA: American Society for Testing and Materials International, 2007: 2-7．

[24] WALDE K, BROCKENBROUGHB J R, CRAIGC B A, et al. Multiple fatigue crack growth in pre-corroded 2024-T3 aluminum[J]. International Journal of Fatigue, 2005, 27(10-12): 1509-1518.

[25] NAN Z Y, ISHIHARA S, GOSHIMA T. Corrosion fatigue behavior of extruded magnesium alloy AZ31 in sodium chloride solution[J]. International Journal of Fatigue, 2008, 30(7): 1181-1188.

[26] MEDVED J J, BRETON M, IRVING P E. Corrosion pit size distributions and fatigue lives—A study of the EIFS technique for fatigue design in the presence of corrosion[J]. International Journal of Fatigue, 2004, 26(1): 71-80.

[27] CHAUSSUMIER M, MABRU C, SHAHZAD M, et al. A predictive fatigue life model for anodized 7050 aluminium alloy[J]. International Journal of Fatigue, 2013, 48(1): 205-213.

[28] GREEN A E, SNEDDON I N. The stress distribution in the neighborhood of a flat elliptic crack in an elastic solid[J]. Proceedings of the Royal Society A, 1946, 187(1009): 229-260.

[29] IRWIN G R. Crack extension force for a part-through crack in a plate[J]. Journal of Applied mechanics, 1962, 29(4): 651-654.

[30] 中國(guó)航空研究院. 應(yīng)力強(qiáng)度因子手冊(cè)[M]. 北京: 科學(xué)出版社， 1993: 162-168.

Chinese Aeronautical Establishment. Handbook of stress intensity factors[M]. Beijing: Science Press, 1993: 162-168 (in Chinese).

[31] XIONG J J, SHENOI R A. Fatigue and fracture reliability engineering[M]. London: Springer, 2011: 86-89.

[32] 北京航空材料研究所. 航空金屬材料疲勞裂紋擴(kuò)展速率手冊(cè)[M]. 北京： 北京航空材料研究所, 1984: 251-256.

Beijing Institue of Aeronautical Material. Handbook of fatigue crack growth rates of aeronautical metallic materials[M]. Beijing: Beijing Institute of Aeronautical Material Press, 1984: 251-256 (in Chinese).