李 真，程立平，李衛(wèi)平

(1.中國商飛上海飛機設(shè)計研究院，上海 201210；2.中航工業(yè)中國飛機強度研究所，西安 710065)

碳纖維復(fù)合材料結(jié)構(gòu)具有比強度高、抗疲勞性能好[1]等優(yōu)點，近年來逐漸應(yīng)用于大型民機機身壁板結(jié)構(gòu)設(shè)計。壁板承受剪切載荷時，首先蒙皮發(fā)生剪切失穩(wěn)，之后壁板進入張力場狀態(tài)，仍然具有較強的承載能力，挖掘壁板失穩(wěn)后的承載能力對飛機減重設(shè)計有重要的意義。

Kuhn等[2-3]在NACA(National Advisory Committee for Aeronautics，美國國家航空咨詢委員會)支持的項目中研究了鋁合金機身單曲率壁板的張力場分析方法，并通過試驗進行了驗證。彭藝琳等[4]研究了鋁鋰合金材料加筋壁板的剪切屈曲性能。孫為民等[5]通過試驗和工程算法研究了金屬機身單曲率壁板在剪切載荷下的承載能力。汪厚冰等[6]通過分析和試驗研究了復(fù)材帽型加筋平直壁板的剪切屈曲性能。Zhang等[7]采用分析和試驗結(jié)合的方法研究了潮濕環(huán)境下復(fù)材壁板的屈曲和后屈曲特性。李真等[8]采用分析和試驗結(jié)合的方法研究了復(fù)材機身壁板在增壓載荷、拉伸、壓縮、剪切載荷作用下的疲勞和損傷容限性能。袁菲等[9]采用有限元軟件ABAQUS定義了壁板及連接界面的失效準則，研究了復(fù)材加筋壁板在壓縮載荷下的屈曲和破壞特性，分析結(jié)果和試驗數(shù)據(jù)對比良好。王喆等[10]采用內(nèi)聚力模型模擬復(fù)材機翼長桁與蒙皮界面失效，研究了機翼結(jié)構(gòu)變厚度長桁的承載能力，結(jié)果表明長桁變厚度設(shè)計可以有效提高結(jié)構(gòu)承載能力。陳昊等[11]設(shè)計了一種新型的試驗裝置，采用有限元分析和試驗方法研究了復(fù)材機身壁板對接結(jié)構(gòu)的靜力、疲勞和損傷容限性能。

上述研究工作主要采用有限元分析和試驗的方法對金屬機身壁板、復(fù)合材料平直壁板的失穩(wěn)載荷進行研究，未開展復(fù)材機身帶曲率壁板的剪切失穩(wěn)載荷研究和失穩(wěn)后承載能力的研究。確定復(fù)材壁板剪切失穩(wěn)載荷和挖掘失穩(wěn)之后的承載能力，對民機結(jié)構(gòu)減重有重要意義。通過不同邊界的工程方法和有限元特征值方法研究復(fù)材壁板的剪切屈曲載荷，將剪切張力場方法應(yīng)用到復(fù)材結(jié)構(gòu)中，研究復(fù)材壁板失穩(wěn)之后的承載能力和失效特性，最后通過3 m×2 m大型機身單曲率壁板試驗進行驗證。

1 壁板試驗

1.1 試驗件設(shè)計介紹

壁板試驗件為包括5個框、8個長桁的單曲率結(jié)構(gòu)，試驗件在飛機航向長度為3 100 mm，環(huán)向?qū)挾葹? 150 mm，半徑為2 960 mm，相鄰框之間距離為620 mm，相鄰長桁之間距離為210 mm，試驗件數(shù)模如圖1所示。試驗件長桁為帽型剖面，長桁與蒙皮通過膠膜共膠接工藝成型。試驗件由蒙皮、長桁、框和剪切剪片組成，各部分的鋪層信息見表1，材料為中模高強碳纖維預(yù)浸料，其材料力學(xué)性能見表2。

表1 蒙皮、長桁等結(jié)構(gòu)的鋪層信息

表2 壁板材料的力學(xué)性能

圖1 壁板試驗件數(shù)模示意圖

試驗件四周為加載區(qū)和過渡區(qū)，中間為考核區(qū)，大尺寸的試驗件保證了考核區(qū)應(yīng)變均勻和支持邊界真實。試驗件加載區(qū)通過粘貼復(fù)材板增強連接強度，再通過螺栓與試驗夾具連接。

1.2 試驗加載裝置

采用機身壁板多軸載荷試驗裝置[12-13]對試驗件進行加載，該裝置可以實現(xiàn)剪切載荷的均勻施加[8]。試驗裝置的長寬高約為6 m×4 m×5 m，如圖2所示。通過在試驗件的四周施加剪流實現(xiàn)剪切載荷的加載，在試驗件的兩個長邊通過連續(xù)的多個作動筒施加主動剪流，另外兩個短邊產(chǎn)生被動的剪流，最終實現(xiàn)在試驗件考核區(qū)施加均勻的剪切載荷，原理如圖3所示。

圖2 試驗裝置照片

圖3 剪切載荷加載原理圖

1.3 試驗過程及結(jié)果

通過應(yīng)變片測量試驗件考核區(qū)的應(yīng)變分布，試驗件考核區(qū)應(yīng)變片布置如圖4所示。以蒙皮預(yù)估失穩(wěn)載荷作為100%載荷進行逐級加載，載荷大小為575.57 kN。典型位置典型蒙皮和長桁應(yīng)變片布置如圖5所示，該圖為試驗件考核區(qū)的剖面圖，蒙皮上下表面背靠背粘貼花片，測量蒙皮航向、環(huán)向應(yīng)變及剪應(yīng)變。長桁上下表面背靠背粘貼單片，測量航向的應(yīng)變分布情況。試驗過程中實時監(jiān)控載荷及應(yīng)變片數(shù)據(jù)，通過應(yīng)變曲線的變化確定蒙皮及長桁的失穩(wěn)載荷。

圖4 試驗件考核區(qū)應(yīng)變片布置

圖5 典型蒙皮和長桁應(yīng)變片布置

首先施加較小載荷進行預(yù)試，確認試驗件安裝及考核區(qū)應(yīng)變的均勻性。根據(jù)預(yù)估失穩(wěn)載荷逐級加載剪切載荷，加載過程中對考核區(qū)應(yīng)變進行實時監(jiān)控，蒙皮應(yīng)變曲線出現(xiàn)分叉[14]，顯示蒙皮失穩(wěn)后進行卸載。載荷卸載后對試驗件進行全面目視檢查，局部采用A掃設(shè)備進行無損檢測，未發(fā)現(xiàn)新增損傷。

最后進行剪切載荷破壞試驗，逐級施加剪切載荷，直至試驗件破壞。試驗件破壞示意圖如圖6所示，壁板出現(xiàn)對角線破壞，失效模式包括蒙皮與長桁發(fā)生脫粘，蒙皮、長桁、框均發(fā)生分層及斷裂。

圖6 試驗件破壞示意圖

試驗件典型蒙皮位置應(yīng)變曲線如圖7所示，圖中給出蒙皮內(nèi)外表面背靠背花片0°、45°、90°的數(shù)據(jù)。典型長桁位置背靠背應(yīng)變曲線如圖8所示。由于加載純剪切載荷，航向、環(huán)向應(yīng)變較小，蒙皮失穩(wěn)之前0°、90°，以及長桁0°片應(yīng)變較小，蒙皮45°的應(yīng)變數(shù)據(jù)隨載荷線性增大。

圖7 典型蒙皮考核區(qū)應(yīng)變曲線

圖8 典型長桁考核區(qū)應(yīng)變曲線

蒙皮失穩(wěn)之后，背靠背的應(yīng)變片數(shù)據(jù)分叉，壁板進入張力場狀態(tài)，載荷發(fā)生重新分配，新增載荷由長桁和框承受。應(yīng)變曲線分叉后，長桁0°的應(yīng)變逐漸增大。蒙皮出現(xiàn)皺褶，發(fā)生局部彎曲變形，0°、45°、90°的應(yīng)變曲線均出現(xiàn)相反的趨勢。壁板破壞前，蒙皮90°片應(yīng)變急劇變大，隨后整個壁板發(fā)生破壞。

根據(jù)試驗結(jié)果，蒙皮失穩(wěn)時對應(yīng)的剪流為Fxyb=228.4 N/mm，剪切應(yīng)變?yōu)棣舃=3 722 με。對屈曲前應(yīng)變進行線性延伸，壁板破壞時剪流為Fxyf=305.3 N/mm，剪應(yīng)變?yōu)棣舊=4 978 με。εf/εb=4 978/3 722=1.34，蒙皮剪切失穩(wěn)后，壁板還具有34%的額外承載能力。

2 機身壁板在剪切載荷下的分析方法

2.1 蒙皮剪切失穩(wěn)工程分析

帽型長桁加筋壁板的蒙皮失穩(wěn)分析時，蒙皮的長度選取框間距，蒙皮的寬度選取帽型長桁R區(qū)之間的距離，如圖9所示。蒙皮剪切失穩(wěn)計算時，將單曲率的蒙皮簡化為平板。

圖9 蒙皮失穩(wěn)計算時寬度定義

采用工程方法計算蒙皮剪切失穩(wěn)載荷[15]，分為四邊簡支和四邊固支邊界條件，計算公式為

(1)

式中:Nxycr為失穩(wěn)載荷；D11、D22為板的彎曲剛度；b為蒙皮寬度；Ks為剪切失穩(wěn)系數(shù)，從文獻[15]中查取。

2.2 蒙皮剪切失穩(wěn)有限元分析

建立壁板試驗件及加載接頭的細化有限元模型，根據(jù)試驗裝置與試驗件連接形式施加載荷和約束支持，采用大型有限元分析軟件MSC.Nastran[16]對試驗件進行線性靜力分析，得到剪切載荷下壁板的位移、應(yīng)變分布，通過有限元特征值分析，得到壁板的剪切失穩(wěn)模態(tài)和失穩(wěn)載荷。

2.3 壁板張力場分析

NACA對鋁合金機身壁板的張力場特性進行了研究[2-3]，蒙皮剪切失穩(wěn)后，壁板進入半張力場狀態(tài)，載荷發(fā)生重新分配，新增載荷由長桁和框承擔(dān)。將其改進并推廣至復(fù)材機身壁板結(jié)構(gòu)，改進后的計算方法如下。

(1)張力場系數(shù)k計算。

(2)

式(2)中：t為蒙皮厚度；d為框距；R為機身半徑；h為桁距；τ為工作應(yīng)力；[τ]為蒙皮失穩(wěn)應(yīng)力。

(2)長桁平均應(yīng)變計算。

(3)

式(3)中:α為張力場角度；Astr為長桁橫截面積；Estr為長桁彈性模量。

(3)框應(yīng)變計算。

(4)

式(4)中:Afra為框橫截面積；Efra為框彈性模量。

(4)蒙皮應(yīng)變計算。

(5)

式(5)中:Eskin為蒙皮彈性模量；ν為蒙皮泊松比。

(5)張力場角度α計算。

(6)

采用迭代法計算張力場角度α:

①假定張力場角度α；

②通過式(3)～式(5)求出εfra、εstr、εskin；

③通過式(6)計算角度α；

④重復(fù)(2)～(5)直到角度α收斂。

蒙皮許用剪應(yīng)變根據(jù)經(jīng)驗曲線得到，該曲線與復(fù)材層板的拉伸強度和張力場系數(shù)k有關(guān)。圖10給出了本試驗件蒙皮鋪層對應(yīng)的許用應(yīng)變曲線。該鋪層在室溫環(huán)境下開孔拉伸失效應(yīng)變?yōu)?0 642 με。其他鋪層的層板對應(yīng)曲線可以根據(jù)該曲線插值得到。

圖10 試驗件蒙皮鋪層的張力場曲線

3 強度分析與試驗對比

3.1 蒙皮剪切失穩(wěn)分析

3.1.1 工程方法

蒙皮剪切失穩(wěn)分析時將曲板簡化為平板，蒙皮長度為620.0 mm，長桁R角之間蒙皮寬度為 145.9 mm。分別采用簡支、固支邊界條件計算蒙皮剪切失穩(wěn)載荷，計算結(jié)果如表3所示。

表3 蒙皮剪切失穩(wěn)工程和有限元分析與試驗結(jié)果對比

簡支邊界條件的分析結(jié)果比壁板試驗結(jié)果偏小29.2%。固支邊界條件的分析結(jié)果比試驗值大4.2%。結(jié)果表明長桁對蒙皮的邊界支持大幅強于簡支，略弱于固支。采用簡支、固支平均的剪切失穩(wěn)結(jié)果比試驗值偏保守12.5%。

3.1.2 有限元分析

壁板試驗件有限元建模采用shell單元，單元尺寸約10 mm，根據(jù)材料強度性能和鋪層順序確定單元屬性，約束和加載根據(jù)試驗設(shè)備加載情況施加，有限元模型示意圖如圖11所示。

圖11 試驗件有限元模型

采用MSC.Nastran進行剪切失穩(wěn)模態(tài)分析，失穩(wěn)載荷如表3所示，失穩(wěn)模態(tài)如圖12所示。分析值比試驗值大9.8%，分析結(jié)果不保守。由于有限元模型中試驗件為理想結(jié)構(gòu)，但試驗件制造時存在制造偏差，比如表面不平、材料屬性分散性、試驗件制造厚度偏差等，因此有限元分析結(jié)果需要考慮一定的折減系數(shù)。

3.2 張力場分析

采用張力場分析方法進行迭代分析，失穩(wěn)應(yīng)變?nèi)≡囼炇Х€(wěn)應(yīng)變，計算得到張力場系數(shù)k=0.184，設(shè)初始張力場角度α=30°，經(jīng)過4次迭代計算結(jié)果逐漸收斂，最終張力場角度為41.1°，迭代后的蒙皮分析應(yīng)變?yōu)棣臿=4 545 με，計算過程及結(jié)果如表4所示。計算結(jié)果表明，長桁的平均應(yīng)變和框的應(yīng)變均較小，蒙皮的應(yīng)變較大，結(jié)合試驗件最終的失效模式判斷蒙皮先發(fā)生失效，導(dǎo)致壁板試驗件發(fā)生失效，重點對蒙皮失效進行分析。

表4 張力場迭代計算結(jié)果

迭代計算后蒙皮的分析應(yīng)變εa比試驗失效應(yīng)變εf略小:εa/εf-1=4 545/4 978-1=-8.7%。根據(jù)圖10的曲線，張力場系數(shù)k=0.184時，對應(yīng)的蒙皮許用應(yīng)變?yōu)閇ε]=4 256 με。許用應(yīng)變小于試驗件實際失效應(yīng)變，[ε]/εf-1=4 256/4 978-1=-0.145，分析結(jié)果有14.5%的保守量。

4 結(jié)論

(1)采用先進大型機身壁板試驗設(shè)備完成了復(fù)材機身單曲率壁板的剪切失穩(wěn)及破壞試驗，研究了試驗件的剪切失穩(wěn)特性和張力場失效特性。

(2)采用工程分析方法計算蒙皮剪切失穩(wěn)載荷，簡支固支平均的分析結(jié)果比試驗結(jié)果略保守；采用有限元特征值方法計算的蒙皮剪切失穩(wěn)載荷大于試驗載荷，需要考慮一定的折減系數(shù)。

(3)壁板剪切失穩(wěn)后，還具有較強的承載能力，采用張力場分析方法分析復(fù)材機身壁板失效載荷可行，結(jié)果偏保守。