摘要：損傷容限評定工作是飛機(jī)結(jié)構(gòu)設(shè)計(jì)和結(jié)構(gòu)壽命估算工作的重要環(huán)節(jié)，本文立足于工程實(shí)際，針對機(jī)翼下壁板，開展了斷裂力學(xué)特性分析和損傷容限方法研究，為結(jié)構(gòu)選型提供基礎(chǔ)數(shù)據(jù)。本文針對飛機(jī)機(jī)翼下壁板對接結(jié)構(gòu)，提出了針對此類結(jié)構(gòu)的損傷容限評定的基本方法和流程，運(yùn)用ANSYS數(shù)值仿真及計(jì)算獲得了裂紋動(dòng)態(tài)擴(kuò)展過程中的應(yīng)力強(qiáng)度因子[2]。本文研究飛機(jī)機(jī)翼下壁板結(jié)構(gòu)損傷容限評定具有工程使用價(jià)值，為飛機(jī)結(jié)構(gòu)采用此類結(jié)構(gòu)的損傷容限評定提供了方法和技術(shù)[3]。

關(guān)鍵詞：應(yīng)力強(qiáng)度因子，損傷容限，機(jī)翼下壁板

1機(jī)翼下壁板結(jié)構(gòu)參數(shù)化設(shè)計(jì)思路

模擬典型T桁條和蒙皮（含凸臺）加筋結(jié)構(gòu)使用ANSYS參數(shù)化設(shè)計(jì)語言（ANSYSParametricDesignLanguage）——APDL。此為一種解釋性語言，可以用來自動(dòng)完成一些適用性強(qiáng)的任務(wù)，也可以用于根據(jù)參數(shù)來建立模型。APDL還包括其它很多特性，諸如重復(fù)某條命令，宏，if-then-else分支，do循環(huán)，標(biāo)量、向量及矩陣操作等[5]。

思路如下：

在ANSYS9.0中建立了機(jī)翼下壁板結(jié)構(gòu)的有限元模型，其中蒙皮和桁條均選取彈性八結(jié)點(diǎn)殼單元SHELL93模擬;鉚釘選用銷釘單元COMBIN7模擬，由于結(jié)構(gòu)的鉚釘非常多且間距與規(guī)格相同，建模過程中采用循環(huán)建立銷釘單元，簡化了建模過程。裂紋尖端選取奇異單元QPE（QuarterPointElement）模擬。

長桁用SHELL93單元（此常用單元參考ANSYS幫助文檔）進(jìn)行模擬，型材各個(gè)部分的厚度取為單元的厚度。蒙皮單元也用SHELL93單元，單元的厚度與蒙皮的實(shí)際厚度相同。有凸臺的位置蒙皮單元厚度為1.8mm，其它單元的厚度為1.5mm。

根據(jù)結(jié)構(gòu)的對稱性，選擇了一半結(jié)構(gòu)用于計(jì)算，參考應(yīng)力： ，P=100KN。蒙皮和長桁都加載，約束取對稱約束。在建模計(jì)算過程中，裂紋建成半裂紋模型，開裂的部分不約束，動(dòng)態(tài)計(jì)算了各種組合情況的應(yīng)力強(qiáng)度因子。

下面機(jī)翼下壁板結(jié)構(gòu)（不含裂紋面的）的靜力模型[7]。

從圖3可以看到，整個(gè)結(jié)構(gòu)的應(yīng)力還是比較均勻，左右應(yīng)力水平也比較對稱，由于模型長度較短，邊界效應(yīng)比較明顯，分析時(shí)候忽略邊界效應(yīng)的影響。蒙皮與桁條都施加載荷，整體結(jié)構(gòu)的傳力比較均勻，釘孔位置的應(yīng)力相對蒙皮與其他位置而言應(yīng)力水平較低，結(jié)構(gòu)設(shè)計(jì)較優(yōu)[6]。

圖4典型T桁條和蒙皮（含凸臺）加筋結(jié)構(gòu)（無裂紋）X方向上的位移云圖

圖4可以看到結(jié)構(gòu)的X向位移也比較對稱。這就為后邊的裂紋按照對稱擴(kuò)展分析奠定了基礎(chǔ)。

2機(jī)翼下壁板結(jié)構(gòu)斷裂力學(xué)參數(shù)計(jì)算

圖5典型T桁條和蒙皮（含凸臺）加筋結(jié)構(gòu)的開裂模式圖給出了結(jié)構(gòu)的六種開裂模式，針對初始開裂模式（含六個(gè)裂紋尖端）定義為開裂模式<1>，對應(yīng)于下圖的有限元模型開裂模式<1>，中間開裂模式<3>對應(yīng)與有限元模型開裂模式<2>。本文只給出兩個(gè)開裂模式有限元模型圖，計(jì)算數(shù)據(jù)還是給出每個(gè)裂紋的應(yīng)力強(qiáng)度因子，由于結(jié)構(gòu)和開裂模式比較對稱，經(jīng)過比較兩個(gè)裂紋尖端的應(yīng)力強(qiáng)度因子還是非常相近的[4]。

圖5中給出了其中兩種開裂模式的有限元模型圖，圖<1>是有六個(gè)裂紋尖端的網(wǎng)格圖，按照對稱分析時(shí)候只需要分析三個(gè)裂紋的擴(kuò)展情況和相互影響。圖<2>因?yàn)锽和C裂紋已經(jīng)和中間對接縫貫穿，所以形成以A和D為裂紋尖端的長裂紋。裂紋長度計(jì)數(shù)法和試驗(yàn)中相對應(yīng)，由于蒙皮在有限元模型中簡化為一個(gè)平面，模型中的角裂紋按照1mm的穿透裂紋處理。

3機(jī)翼下壁板結(jié)構(gòu)的應(yīng)力強(qiáng)度因子計(jì)算

求解應(yīng)力強(qiáng)度因子的方法有多種，但對于比較復(fù)雜的結(jié)構(gòu)中的裂紋問題，有效的方法還是有限元法。本文采用有限元求解應(yīng)力強(qiáng)度因子的方法直接法，基于1/4節(jié)點(diǎn)單元的位移法，這種特殊單元能夠反映裂紋尖端的奇異性，不需要用過細(xì)的網(wǎng)格就可以得到比較精確的結(jié)果。通過ANSYS計(jì)算的應(yīng)力強(qiáng)度因子。相關(guān)曲線如圖6是蒙皮無裂紋的不同長桁裂紋的無量綱應(yīng)力強(qiáng)度因子曲線;圖7是長桁上無裂紋蒙皮上裂紋B和C已經(jīng)貫穿，只有A和D裂紋的無量綱的應(yīng)力強(qiáng)度因子曲線。

4下壁板結(jié)構(gòu)不同開裂模式下的應(yīng)力強(qiáng)度因子分析

圖7給出的曲線是蒙皮不含裂紋的不同長桁裂紋的無量綱的應(yīng)力強(qiáng)度因子曲線。當(dāng)桁條裂紋E和F還沒有貫穿的時(shí)候，E和F裂紋尖端的應(yīng)力強(qiáng)度因子基本相近，便于作曲線只選取較大值得一端。在曲線中前三點(diǎn)表示的就是還有兩個(gè)裂紋尖端的無量綱應(yīng)力強(qiáng)度因子，曲線凸峰最高點(diǎn)是桁條裂紋擴(kuò)展到型材拐角處并且已經(jīng)貫通成一個(gè)裂紋G的數(shù)據(jù)點(diǎn)。曲線的趨勢和《民機(jī)損傷容限設(shè)計(jì)手冊》[1]典型蒙皮和“工”型桁條加筋結(jié)構(gòu)的趨勢一致。

5結(jié)束語

對有限元方法中在裂尖使用奇異單元計(jì)算應(yīng)力強(qiáng)度因子的方法進(jìn)行了歸納和總結(jié)，應(yīng)用有限元分析軟件ANSYS計(jì)算了新型渦扇飛機(jī)機(jī)翼下壁板結(jié)構(gòu)不同開裂模式的應(yīng)力強(qiáng)度因子，本章給出了典型“T”字型桁條和蒙皮（含凸臺）加筋結(jié)構(gòu)的開裂模式，并且計(jì)算了兩種開裂模式下蒙皮裂紋和長桁裂紋組合情況下的應(yīng)力強(qiáng)度因子，并將其無量綱化和繪出曲線圖。為后邊的裂紋擴(kuò)展分析和剩余強(qiáng)度分析提供了基礎(chǔ)的數(shù)據(jù)。

參考文獻(xiàn)：

[1]《民機(jī)耐久性及損傷容限設(shè)計(jì)手冊》，航空工業(yè)出版社，2003

[2]《損傷容限設(shè)計(jì)手冊》，西北工業(yè)大學(xué)出版社，1989

[3]《飛機(jī)耐久性及損傷容限設(shè)計(jì)手冊》，航空航天部技術(shù)研究院

[4]傅祥炯，《結(jié)構(gòu)疲勞與斷裂》，西北工業(yè)大學(xué)出版社，1995

[5]《FRACTUREANALYSISOFTHEFAA/NASAWIDESTIFFENEDPANELS》

[6]《飛機(jī)結(jié)構(gòu)強(qiáng)度》，北京航空航天大學(xué)出版社，2005

[7]《TheMechanicsofFractureandFatigue》，

作者簡介：張妮娜，1982.4.20出生，籍貫：會寧，民族，漢，性別：女，學(xué)歷：碩士研究生。畢業(yè)院校：西北工業(yè)大學(xué)，研究方向：固體力學(xué)