某飛行器復合材料薄壁加筋結(jié)構(gòu)艙段穩(wěn)定性分析*

劉文一　焦冀光

(91550部隊　大連　116023)

?

某飛行器復合材料薄壁加筋結(jié)構(gòu)艙段穩(wěn)定性分析*

劉文一焦冀光

(91550部隊大連116023)

摘要為了研究復合材料薄壁加筋結(jié)構(gòu)艙段的在特定載荷下的穩(wěn)定性,采用有限元方法仿真計算了在艙段薄壁厚度不變時,筋條截面積在200mm2、300mm2和400mm2下,艙段穩(wěn)定性隨著筋條高度變化情況,得到了臨界載荷,計算結(jié)果可作為工程設(shè)計參考依據(jù)。

關(guān)鍵詞艙段; 屈曲; 有限元法

Stability Analysis of A Aircraft’s Thin-walled Reinforced Composite Material Cabin

LIU WenyiJIAO Jiguang

(No. 91550 Troops of PLA, Dalian116023)

AbstractThe finite element method has been used to study the stability of thin-walled reinforced composite material cabin under a particular cabin load. The change behavior of cabin stability has been find when the rib’s height has been changed as it’s was 200mm2, 300mm2and 400mm2, though the thickness of the thin wall was constant. the critical load has been figured out, the results could be used as reference for engineering design.

Key Wordscabin, bulking, finite element method

Class NumberV212

1引言

結(jié)構(gòu)穩(wěn)定性分析是飛行器結(jié)構(gòu)設(shè)計中最重要的問題之一,因為復合材料薄壁加筋結(jié)構(gòu)的靜強度破壞中有很大一部分是由于喪失穩(wěn)定性引起的。由于結(jié)構(gòu)穩(wěn)定性的限制,飛行器結(jié)構(gòu)的設(shè)計應力往往遠小于結(jié)構(gòu)材料的極限強度。所以,保持穩(wěn)定性是不同復合材料薄壁加筋結(jié)構(gòu)形式的選擇和設(shè)計的主要依據(jù)[1]。

根據(jù)大量的試驗表明,結(jié)構(gòu)的失效除少數(shù)是由局部強度破壞外,多數(shù)是由結(jié)構(gòu)屈曲而引起的。在設(shè)計最初給出臨界屈曲載荷是飛行器結(jié)構(gòu)設(shè)計的輸入條件。對結(jié)構(gòu)進行屈曲分析,涉及到較復雜的彈(塑)性理論和數(shù)學計算,要通過求解高階偏微分方程組,才能求得精確的解析解。因此必須作大量簡化才能用數(shù)值方法求解。但是這種簡化勢必會帶來精度低、結(jié)果不可信等后果。隨著計算機和有限元方法的迅猛發(fā)展,形成了許多的實用分析程序,提高了對復雜結(jié)構(gòu)進行屈曲分析的能力和設(shè)計水平。分析計算結(jié)果時,應運用工程經(jīng)驗,避開實際中可能發(fā)生的屈曲模態(tài)所對應的屈曲載荷。

2計算模型

某飛行器艙段采用C/SiC復合材料薄壁加筋整體加工,基體三維編織復合材料克服了傳統(tǒng)層合板復合材料所固有的厚度方向的剛度和強度低、面內(nèi)和層間剪切強度低、易分層的弱點,具有結(jié)構(gòu)質(zhì)量輕,整體強度高等特點,是新型飛行器艙段的理想材料。復合材料力學性能參數(shù)如表1所示。

對于復合材料薄壁加筋結(jié)構(gòu)艙段,可能出現(xiàn)的屈曲現(xiàn)象有薄壁失穩(wěn)、筋條失穩(wěn)和整體失穩(wěn)。對于前兩種失穩(wěn)模式,可認為是沒有發(fā)揮薄壁加筋條結(jié)構(gòu)特有的穩(wěn)定性而需要避免的[2]。對于整體失穩(wěn)模式,則可認為是發(fā)揮了薄壁加筋條結(jié)構(gòu)穩(wěn)定性,但這種屈曲模態(tài)需避開實際飛行載荷,從而避免飛行器在工作過程中結(jié)構(gòu)失穩(wěn)[3]。

表1　C/SiC復合材料的纖維和基體主要力學性能

飛行器在飛行工作過程中主要受軸向壓力載荷,對艙段屈曲分析時將其底部端框面的六個自由度全部約束,將壓力載荷施加在前端框面上,建立的有限元模型如圖1所示。

圖1　復合材料艙段有限元模型

3計算理論

在有限元方法中,線形穩(wěn)定性(屈曲)理論通過提取使線性系統(tǒng)剛度矩陣奇異的特征值來獲得結(jié)構(gòu)的失穩(wěn)載荷即失穩(wěn)模態(tài)[4]。線性屈曲分析的特點是變形為小變形,單元應力必須是彈性的,載荷引起的內(nèi)部單元力分布是不變的。離散后有限元模型的屈曲分析控制方程為[5]

Ku=λKdu

(1)

式中:K為剛度矩陣,Kd為剛度微分矩陣,λ為特征值,u為特征向量。

這個方程可簡化為一種標準形式求解。采用Lanczos法對其實特征值求解,將式(1)改寫為[6～7]

(2)

選擇初始迭代矢量:

(3)

依次取j=1,2,3…,直到出現(xiàn)收斂的迭代根,迭代函數(shù)為

rj+1=(K-σKd)-1Kdqj-αjqj-βjqj-1

(4)

其中:

(5)

將αj、βj合并為一個變換矩陣Tj,如下所示:

(6)

設(shè)θ為特征方程Tjs=sθ的特征值根,則式(2)的特征值為

(7)

所以λ為使線性系統(tǒng)剛度矩陣奇異的特征值,即為所求艙段的復合材料薄壁加筋結(jié)構(gòu)的臨界壓力[8～9]。

4計算結(jié)果及分析

為了研究復合材料薄壁加筋結(jié)構(gòu)艙體的結(jié)構(gòu)參數(shù)對屈曲載荷的影響,把握薄壁加筋結(jié)構(gòu)的失效規(guī)律,在薄壁厚度和加筋肋的截面積不變的情況下,通過改變筋條的高度和寬度研究結(jié)構(gòu)的屈曲載荷和屈曲形式的規(guī)律[10～11]。

表2　筋條幾何參數(shù)變化對結(jié)構(gòu)失穩(wěn)形式的影響

根據(jù)設(shè)計要求,艙段的薄壁厚度不變,要求得到在筋條橫截面積為200mm2、300mm2和400mm2時的最大穩(wěn)定載荷。計算時設(shè)置了在三種橫截面積下,筋條高度變化分別為10mm、13mm、15mm、18mm和20mm。表2給出了三種不同橫截面面積筋條幾何參數(shù)變化時,臨界失穩(wěn)載荷及失穩(wěn)模式的計算結(jié)果。

從表2中可以看出:在筋條橫截面積不變的情況下,隨著筋條高度的增加,薄壁失穩(wěn)、整體失穩(wěn)和筋條失穩(wěn)這三種失穩(wěn)形式將依次出現(xiàn)[12]。

筋條橫截面積為200mm2時,筋條失穩(wěn)和薄壁失穩(wěn)均出現(xiàn)兩次,而整體失穩(wěn)僅在筋高為15mm時出現(xiàn)一次;筋條高度為18mm時,出現(xiàn)筋條局部失穩(wěn),而筋條高度為20mm時,則出現(xiàn)筋條全部失穩(wěn),而臨界載荷則由389.26MPa降低至323.45MPa,如圖2所示。

圖2　筋條的兩種失穩(wěn)形式

筋條橫截面積為300mm2時筋條失穩(wěn)僅出現(xiàn)一次,整體失穩(wěn)在筋高為18mm和15mm時均出現(xiàn),薄壁失穩(wěn)在筋高13mm和10mm時均出現(xiàn);

筋條橫截面積為400mm2時,筋條失穩(wěn)出現(xiàn)一次,整體失穩(wěn)在筋高為18mm時僅出現(xiàn)一次,而薄壁失穩(wěn)則出現(xiàn)三次,但失穩(wěn)形式均為局部失穩(wěn)

在所有的艙段整體失穩(wěn)模式中,失穩(wěn)形式差別不大;在所有的薄壁失穩(wěn)模式中,均出現(xiàn)薄壁局部失穩(wěn)形式,如圖3所示。

圖3　艙段失穩(wěn)形式

圖4給出了三種筋條截面積下,臨界失穩(wěn)載荷隨筋條高度的變化趨勢。

圖4　臨界失穩(wěn)載荷隨筋條高度的變化

從圖2可以看出,隨著筋條高度的增大,三種截面積艙段的臨界失穩(wěn)載荷均出現(xiàn)先升高后降低的變化趨勢。

截面積為200mm2時,最大臨界載荷出現(xiàn)在筋高為15mm時,其值為423.21MPa,此時艙段的失穩(wěn)形式為整體失穩(wěn)。

截面積為300mm2時,最大臨界載荷出現(xiàn)在筋高為18mm時,其值為501.98MPa,此時艙段的失穩(wěn)形式為整體失穩(wěn)。

截面積為400mm2時,最大臨界載荷出現(xiàn)在筋高為18mm時,其值為553.24MPa,此時艙段的失穩(wěn)形式為整體失穩(wěn)。

5結(jié)語

1) 在同一筋條橫截面積下,隨著筋條高度的增加,依次出現(xiàn)薄壁失穩(wěn)、整體失穩(wěn)和筋條失穩(wěn)這三種形式,截面積為300mm2時艙段整體屈曲形式出現(xiàn)兩次,其余兩種截面積下艙段整體屈曲形式各出現(xiàn)一次;

2) 隨著筋條截面積的增加,艙體發(fā)生整體失穩(wěn)時的臨界失穩(wěn)載荷將隨之增大,這表明增大筋條截面積,有利于提高艙段整體穩(wěn)定性;

3) 綜合計算結(jié)果分析,可將筋條截面積為200mm2,筋高為15mm;筋條截面積為300mm2,筋高為15mm和18mm;筋條截面積為400mm2,筋高為18mm時的計算結(jié)果作為工程設(shè)計參考依據(jù)。

參 考 文 獻

[1] 吳存利,師成梅,段世慧,等.應用MSC.Nastran計算飛機壁板穩(wěn)定性問題[C]//2005年MSC.Software中國用戶論文集,2005:253-254.

[2] 龍凱,賈長治,李寶峰.Patran與Nastran有限元分析從入門到精通[M].北京:機械工業(yè)出版社,2011:281-292.

[3] 羅忠,陳志堅.特殊耐壓殼體艙段穩(wěn)定性有限元分析[J].2005年MSC.Software中國用戶論文集,2005:12-13.

[4] 陳志堅,郭日修.加肋軸對稱旋轉(zhuǎn)殼非線性穩(wěn)定性分析[J].中國造船,2004(1):5-6.

[5] 陳鐵云,陳伯真.彈性薄殼力學[M].武漢:華中工學院出版社,1983:125-126.

[6] 黃英蘭,徐元銘.基于MSC.Nastran的復合材料格柵壁板屈曲數(shù)值仿真[J].2007年MSC.Software中國用戶論文集,2007:78-79.

[7] 李金龍,楊萬均,等.復合材料發(fā)射器透濕率試驗研究[J].彈箭與制導學報,2014(3):45-46.

[8] 朱晉生,王卓,歐峰.先進復合材料在航空航天領(lǐng)域的應用[J].新技術(shù)新工藝,2012(10):6-8.

[9] 何子奇,周緒紅,劉占科,等.冷彎薄壁卷邊槽鋼軸壓構(gòu)件畸變與局部相關(guān)屈曲試驗研究[J].建筑結(jié)構(gòu)學報,2013(11):5-7.

[10] 田春雷,皮愛國,黃風雷.碳纖維復合材料力學性能及其復合彈體侵徹試驗研究[J].彈箭與制導學報學報,2013(3):15-16.

[11] 王亮,商霖,牛智玲,等.基于NASTRAN導彈儀器艙噪聲環(huán)境預示研究[J].導彈與航天運載技術(shù),2013(2):16-17.

[12] 崔釗,李建波,王俊超,等.加裝格尼襟翼的旋翼機配平及穩(wěn)定性研究[J].飛行力學,2013(5):113-114.

中圖分類號V212

DOI:10.3969/j.issn.1672-9730.2016.02.030

作者簡介:劉文一,男,碩士,工程師,研究方向:航天器結(jié)構(gòu)設(shè)計和強度分析。

*收稿日期:2015年8月7日,修回日期:2015年9月23日