牟秀軍，劉小忠，劉 濤，陳定海

（1.91049部隊，山東青島266102；2.海軍航空工程學(xué)院青島校區(qū)，山東青島266041）

武器裝備的主體結(jié)構(gòu)一般都是金屬材料，金屬構(gòu)件在機(jī)械加工、裝配、使用過程中，會不可避免產(chǎn)生各種缺陷，如劃傷、壓痕和腐蝕等，這些缺陷在交變載荷作用下形成疲勞裂紋源[1-3]，當(dāng)出現(xiàn)初始裂紋后，在疲勞載荷作用下，裂紋緩慢擴(kuò)展，達(dá)到臨界裂紋尺寸后，發(fā)生破壞，這將直接影響武器裝備可靠性、作戰(zhàn)性能。

虛擬裂紋閉合法（Virtual Crack Closure Technique，VCCT）最早被Rybicki 和Kanninen[4]在1977年用來解決線性裂紋問題，之后被Shivakumar Tan 和Newman[5]解決表面裂紋問題。Raju[6]從數(shù)學(xué)基礎(chǔ)理論對VCCT進(jìn)行研究，并將這種技術(shù)擴(kuò)展到任意類型的單元里。虛裂紋閉合法與有限元相結(jié)合[7-9]，其基礎(chǔ)理論是裂紋張開釋放出的能量與外力使裂紋閉合所做的功相等，其利用節(jié)點位移和節(jié)點力計算應(yīng)變能釋放率和應(yīng)力強(qiáng)度因子等裂紋擴(kuò)展參數(shù)[10-11]。

本文利用虛擬裂紋閉合法分析疲勞裂紋擴(kuò)展過程，為損傷容限評估方法提供一種技術(shù)支持。主要思想是：基于虛擬裂紋閉合法，利用有限元軟件Abaqus，編寫出用戶子單元程序，計算出不同裂紋擴(kuò)展階段的應(yīng)變能釋放率和裂紋尖端應(yīng)力強(qiáng)度因子的數(shù)值，仿真分析疲勞裂紋擴(kuò)展過程，并與預(yù)腐蝕鋁合金6A02 試驗件疲勞試驗得到的裂紋長度與疲勞壽命數(shù)據(jù)進(jìn)行對比分析。

1 疲勞裂紋擴(kuò)展試驗

疲勞裂紋擴(kuò)展試驗采用鋁合金6A02 材料，該鋁合金屬于Al-Mg-Si-Cu系，其主要成分，如表1所示，試驗件形狀尺寸如圖1所示。

表1 鋁合金6A02成分Tab.1 Chemical constitution of Aludur6A02

圖1 試驗件設(shè)計示意圖（h=7 mm）Fig.1 Schematic of the specimen（h=7 mm）

預(yù)腐蝕試驗：腐蝕溶液配制按95份蒸餾水與5份NaCl（分析純）配置濃度為5%的NaCl溶液，再加入適量的稀硫酸，使pH=4.0±0.2，將試驗件垂直放置在絕緣塑料托架的卡槽內(nèi)浸沒腐蝕溶液中，預(yù)腐蝕72 h后取出試驗件，清洗后試驗件表面如圖2所示。

圖2 預(yù)腐蝕后試驗件表面Fig.2 Configuration of corrosion specimen

疲勞裂紋擴(kuò)展試驗：利用MTS810 疲勞試驗機(jī)塊譜模式在3 種不同應(yīng)力水平進(jìn)行疲勞試驗，疲勞斷口如圖3所示。

依據(jù)斷口反推技術(shù)[12-13]得到不同應(yīng)力水平下裂紋長度與疲勞壽命數(shù)據(jù)，如表2所示。

圖3 試驗件疲勞斷口圖Fig.3 Fatigue fracture of specimen

表2 不同應(yīng)力水平下疲勞斷口裂紋(a,N)數(shù)據(jù)統(tǒng)計Tab.2 Data of fatigue fracture crack with different stress

2 VCCT疲勞單元建立

虛擬裂紋閉合法的基本思想是裂紋從a擴(kuò)展到a+Δa所需的功與裂紋從a+Δa閉合到a所需的功相等。圖4 顯示長度為a的裂紋和其增長量Δa的虛擬裂紋擴(kuò)展。

圖4 虛擬裂紋擴(kuò)展圖Fig.4 Figure of virtual crack propagation

Irwin的裂紋閉合積分為：

式（1）、（2）中：σyy和τxy為長度為a裂紋尖端的應(yīng)力分量；Δu和Δv為虛擬裂紋面上點的相對滑離位移和張開位移；B為裂紋厚度常量；Δa為長度為a裂紋擴(kuò)展量；r為裂紋擴(kuò)展積分變量；GI和GII為裂紋模式I和模式II下的應(yīng)變能釋放率分量。

采用虛擬裂紋擴(kuò)展技術(shù)（見圖5），考慮裂紋尖端附近的有限元網(wǎng)格，虛擬裂紋線上分布的應(yīng)力做的功等于節(jié)點力在節(jié)點位移上做的功，即

式（3）中：Fy1表示節(jié)點1上y方向的節(jié)點力；v1,1′表示節(jié)點1和1′間y方向的變化。

利用VCCT 得到裂紋擴(kuò)展應(yīng)變能釋放率后，應(yīng)力強(qiáng)度因子可計算得到平面應(yīng)力和平面應(yīng)變：

式（4）、（5）中：E為彈性模量；μ為泊松比；KΙ和KΙΙ為裂紋模式I和模式II下的應(yīng)力強(qiáng)度因子。

基于虛擬裂紋閉合法進(jìn)行疲勞裂紋擴(kuò)展壽命預(yù)測，結(jié)合有限元軟件Abaqus，編寫出VCCT 的用戶子單元，從而實現(xiàn)疲勞裂紋擴(kuò)展壽命分析。思路是在每一個裂紋增量步中計算出裂尖應(yīng)力強(qiáng)度因子，再根據(jù)疲勞裂紋擴(kuò)展速率公式，計算確定裂紋擴(kuò)展一個單元長度的循環(huán)次數(shù)。計算分析流程如圖6所示。

圖6 用戶子單元計算流程圖Fig.6 Flow chart of UEL

對于裂紋長度與裂紋擴(kuò)展壽命的差分關(guān)系是

在有限元模型中裂紋增長量Δa為界面2 節(jié)點的距離，裂紋擴(kuò)展一個單元長度對應(yīng)的裂紋擴(kuò)展壽命是

由于應(yīng)力強(qiáng)度因子ΔK與裂紋長度a有關(guān)，疲勞裂紋擴(kuò)展壽命對有限元分析中網(wǎng)格大小比較敏感，越精密的網(wǎng)格，其Δa很小，這樣預(yù)測的疲勞裂紋擴(kuò)展壽命越精確，但其建立模型很繁瑣，計算的代價很高。利用較大的網(wǎng)格獲得較高的計算精度，即對應(yīng)力強(qiáng)度因子利用形狀修正函數(shù)進(jìn)行插值，經(jīng)過仿真計算分析發(fā)現(xiàn)，插值點取中間值時就能顯著降低網(wǎng)格敏感性影響。因此，插值點的應(yīng)力強(qiáng)度因子為

3 疲勞裂紋擴(kuò)展壽命仿真分析

通過對試驗件疲勞斷口觀察發(fā)現(xiàn)，疲勞裂紋萌生于孔的周圍腐蝕坑斑處（如圖7所示，其中A處發(fā)亮為鋁合金的氧化產(chǎn)物）。因此，利用有限元軟件Abaqus軟件模擬試驗件，無須對整個試驗件進(jìn)行網(wǎng)格細(xì)化，只須對孔的周圍進(jìn)行網(wǎng)格細(xì)化（如圖8 所示），該網(wǎng)格細(xì)化在滿足精度要求下，大大節(jié)省了分析運(yùn)行的時間，單元類型選取為CPS4。

圖7 疲勞斷口的腐蝕疲勞源圖Fig.7 Fatigue source of fatigue fracture

圖8 有限元模型Fig.8 Finite element model

考慮到有限元網(wǎng)格劃分和計算效率，在結(jié)合虛擬裂紋閉合法對裂紋擴(kuò)展能量釋放率計算模擬時，設(shè)初始裂紋預(yù)置位0.25 mm，每步擴(kuò)展0.25 mm，共擴(kuò)展15步，通過仿真模擬得到鋁合金6A02 裂紋擴(kuò)展不同階段的應(yīng)力云圖，如圖9所示。

圖9 不同裂紋擴(kuò)展階段應(yīng)力云圖Fig.9 Stress cloud of different crack propagatio

疲勞裂紋擴(kuò)展壽命進(jìn)行仿真預(yù)測，其疲勞裂紋擴(kuò)展速率公式運(yùn)用Paris公式，并利用VCCT方法建立的用戶子單元，分別計算不同應(yīng)力水平下裂紋每步擴(kuò)展的裂紋擴(kuò)展能量釋放率（如圖10所示），通過能量釋放率計算出裂紋擴(kuò)展每步的應(yīng)力強(qiáng)度因子幅值（如圖11所示），結(jié)合式（6）~（8），得到裂紋每步擴(kuò)展對應(yīng)的循環(huán)次數(shù)，也就獲得不同應(yīng)力水平下的裂紋擴(kuò)展曲線，并與裂紋擴(kuò)展試驗過程進(jìn)行對比，如圖12所示。

圖10 在不同應(yīng)力水平下裂紋擴(kuò)展能量釋放率Fig.10 Release rate of crack propagation energy with different stress

圖11 不同應(yīng)力水平下裂紋尖端應(yīng)力強(qiáng)度因子幅值Fig.11 Amplitude of SIF on crack tip with different stress

圖12 試驗與仿真的裂紋擴(kuò)展曲線對比圖Fig.12 Comparison chart of crack growth curve with experiment and simulation

在圖12中，疲勞試驗得到的裂紋擴(kuò)展數(shù)據(jù)是以循環(huán)次數(shù)為計數(shù)單位測量裂紋長度，圖中在相同循環(huán)次數(shù)下應(yīng)力越大裂紋擴(kuò)展越快；而VCCT疲勞單元仿真得到的裂紋擴(kuò)展數(shù)據(jù)是以裂紋擴(kuò)展長度為計數(shù)單位計算裂紋擴(kuò)展壽命，圖中在相同裂紋擴(kuò)展長度下應(yīng)力越小裂紋擴(kuò)展壽命越長。因此，從圖中可以對比試驗值與仿真值分析得到，利用VCCT疲勞單元仿真得到裂紋擴(kuò)展趨勢與試驗結(jié)果基本吻合，建立的基于虛擬裂紋閉合法的有限元裂紋擴(kuò)展模型，可仿真再現(xiàn)疲勞裂紋擴(kuò)展過程，預(yù)測裂紋擴(kuò)展壽命。

4 小結(jié)

結(jié)合虛擬裂紋閉合法，編寫用戶單元子程序，計算出不同裂紋擴(kuò)展階段的應(yīng)變能釋放率和裂紋尖端應(yīng)力強(qiáng)度因子，并利用Paris公式建立有限元裂紋擴(kuò)展模型，并與通過對預(yù)腐蝕鋁合金6A02 材料進(jìn)行裂紋擴(kuò)展試驗得到疲勞裂紋長度與疲勞壽命數(shù)值進(jìn)行對比分析，得到基于虛擬裂紋閉合法建立的有限元裂紋擴(kuò)展模型可仿真再現(xiàn)疲勞裂紋擴(kuò)展過程，預(yù)測裂紋擴(kuò)展壽命，為損傷容限評估方法提供一種技術(shù)支持。

[1] 李旭東，劉治國，穆志韜.基于短裂紋的LD10CZ鋁合金腐蝕預(yù)疲勞裂紋擴(kuò)展[J]. 海軍航空工程學(xué)院學(xué)報，2013，28（1）：51-56.

LI XUDONG，LIU ZHIGUO，MU ZHITAO. Fatigue crack growth of LD10CZ aluminum alloy with pre-corrosion based on micro-crack[J].Journal of Naval Aeronautical and Astronautical University，2013，28（1）：51-56.（in Chinese）

[2] 陳定海，穆志韜，朱做濤，等.腐蝕坑對疲勞裂紋擴(kuò)展的影響分析[J].裝備環(huán)境工程，2012，26（4）：4-8.

CHEN DINGHAI，MU ZHITAO，ZHU ZUOTAO，et al.Analysis of corrosion pits on fatigue crack growth[J].Equipment Environmental Engineering，2012，26（4）：4-8.（in Chinese）

[3] 金平，楊凱，薛慶増.腐蝕環(huán)境下材料的疲勞缺口系數(shù)[J].海軍航空工程學(xué)院學(xué)報，2010，25（5）：485-487.

JIN PING，YANG KAI，XUE QINGZENG. Fatigue notch factor of corrosion material[J].Journal of Naval Aeronautical and Astronautical University，2010，25（5）：485-487.（in Chinese）

[4] RYBICKI E F，KANNINEN M F.A finite element calculation of stress intensity factors by a modified crack closure integral[J]. Engineering Fracture Mechanics，1977，9：931-938.

[5] SHIVAKUMAR K N，TAN P W，NEWMAN J J.A virtual crack-closure technique for calculating stress intensity factors for cracked three dimensional bodies[J]. International Journal of Fracture，1988，36：43-50.

[6] RAJU I S. Calculation of strain energy release rates with higher order and singular finite elements[J]. Engineering Fracture Mechanics，1987，28：251-274.

[7] KRUEGER R. The virtual crack closure technique history，approach and applications，2002-10[R].USA：ICASE，2002.

[8] XIE D，SHERRILL B，BIGGERS J. Strain energy release rate calculation for a moving delamination front of arbitrary shape on the virtual crack closure technique. Part I：Formulation and validation[J]. Engineering Fracture Mechanics，2006，73（6）：771-775.

[9] XIE D，SHERRILL B，BIGGERS J. Strain energy release rate calculation for a moving delamination front of arbitrary shape on the virtual crack closure technique.Part II：Sensitivity study on modeling details[J]. Engineering Fracture Mechanics，2006，73（6）：771-775.

[10]肖濤，左正興，劉棟，等.基于虛擬裂紋閉合法計算裂紋擴(kuò)展的能量釋放率[J]. 北京理工大學(xué)學(xué)報，2010，30（1）：37-41.

XIAO TAO，ZUO ZHENGXING，LIU DONG，et al.Computation of the crack propagation engergy release rate based on the virtual crack closure technique[J].Journal of Beijing Institute of Technology，2010，30（1）：37-41.（in Chinese）

[11]魯國富，劉勇，張呈林.基于虛擬裂紋閉合技術(shù)的應(yīng)變能釋放率分析[J].復(fù)合材料學(xué)報，2009，26（2）：210-216.

LU GUOFU，LIU YONG，ZHANG CHENGLIN. Analysis of strain energy release rate based on virtual crack closure technique[J]. Journal of Acta Materiae Compositae Sinica，2009，26（2）：210-216.（in Chinese）

[12]劉新靈，張崢，陶春虎.疲勞斷口定量分析[M].北京：國防工業(yè)出版社，2010：54-83.

LIU XINLING，ZHANG ZHENG，TAO CHUNHU. Fatigue fractography quantitative analysis[M]. Beijing：National Defense Industry Press，2010：54-83.（in Chinese）

[13]王勝霞，竇松柏.斷口定量分析在直升機(jī)關(guān)鍵動部件疲勞試驗分析中的應(yīng)用[J].直升機(jī)技術(shù)，2012，2：19-25.

WANG SHENGXIA，DOU SONGBAI. Application of quantitative fractography analysis in the helicopter parts fatigue test[J]. Helicopter Technique，2012，2：19-25.（in Chinese）