劉　衛(wèi)　何周理 /

(1. 中國商飛北京民用飛機(jī)技術(shù)研究中心，北京 102209；2. 上海飛機(jī)設(shè)計(jì)研究院，上海 201210)

0　引言

由蒙皮、長桁、框組成的半硬殼式壁板是大型商用飛機(jī)常見的結(jié)構(gòu)。帽型長桁的截面尺寸較大，其兩邊突緣與蒙皮相連形成一個(gè)閉合截面，具有很高的受壓穩(wěn)定性[1],是壁板結(jié)構(gòu)中效率較高的縱向加強(qiáng)件。但帽型長桁為封閉結(jié)構(gòu)，內(nèi)部缺陷不易檢查，內(nèi)腔容易積液而產(chǎn)生腐蝕，因此在金屬半硬殼式壁板中，帽型長桁的使用受到了極大的限制。隨著比強(qiáng)度高、耐腐蝕性好的碳纖維復(fù)合材料大量應(yīng)用于現(xiàn)代商用飛機(jī)上，帽型長桁成為了設(shè)計(jì)師首選的縱向加強(qiáng)件。

國內(nèi)外針對(duì)金屬和復(fù)合材料壁板優(yōu)化已經(jīng)做了大量研究工作。Kidane等采用了一種典型單元體模型分析格柵加筋筒結(jié)構(gòu)總體屈曲[2]，Jaunky等采用改進(jìn)的等效法分析了不同類型壁板結(jié)構(gòu)的穩(wěn)定性，且結(jié)合職能算法對(duì)復(fù)合材料格柵加筋結(jié)構(gòu)的優(yōu)化問題進(jìn)行了研究[3]，常楠、張鐵亮等對(duì)加筋采用兩級(jí)方法進(jìn)行優(yōu)化分析[4-5]，張國柱等通過代理模型進(jìn)行金屬加筋板的布局優(yōu)化[6]。但在初步尺寸設(shè)計(jì)[7]時(shí)，帽型長桁的截面尺寸對(duì)設(shè)計(jì)師來說顯得尤為重要。模型的適當(dāng)簡化是縮短研制時(shí)間的重要途徑之一，本文通過局部屈曲和壓損工程算法來優(yōu)化帽型截面尺寸，并得出：(1)帽型長桁局部屈曲載荷與帽底寬度，帽腰和帽底夾角的曲線；(2)壓損載荷與帽底寬度，帽腰和帽底夾角的曲線，為設(shè)計(jì)師在初步尺寸設(shè)計(jì)時(shí)提供參考，縮短研發(fā)時(shí)間。

1　幾何模型

選取一段長200 mm的共固化蒙皮、帽型長桁作為研究對(duì)象，其截面幾何如圖1所示。

圖1　蒙皮、長桁截面

把帽底寬度W，帽腰和帽底夾角α作為變量，α取值范圍為50°～80°，W的取值范圍為60 mm～80 mm。

蒙皮鋪層為:

[45/0/0/- 45/0/90/90/0/- 45/0/0/45]。

長桁鋪層為:

[45/0/- 45/0/90/0/- 45/0/45]。

材料：選取典型的T300碳纖維材料作為研究對(duì)象，材料參數(shù)見表2。

2　局部屈曲

2.1 局部屈曲計(jì)算方法

長桁局部屈曲載荷計(jì)算方法采用參考文獻(xiàn)[8]5.7.1.3節(jié)和5.7.1.5節(jié)中介紹的分析方法，分析模型如圖2所示。

圖2　加筋壁板長桁中的腹板

開口薄壁剖面長桁的突緣可作一長邊自由，另一長邊簡支的長板處理，按式(1)計(jì)算突緣的軸壓局部屈曲載荷。

(1)

式中，Nxcr為單位寬度上的軸壓屈曲載荷；bt為突緣的寬度；L為長桁的長度；D11、D66為層壓板的縱向彎曲剛度和扭轉(zhuǎn)剛度。

對(duì)于薄壁長桁的腹板，可當(dāng)作兩長邊簡支的長板處理，計(jì)算局部屈曲載荷

(2)

式中，bf為腹板的寬度；D12、D22為—層壓板的泊松彎曲剛度和橫向彎曲剛度。

取突緣和腹板局部屈曲載荷的最小值作為長桁的局部屈曲載荷。

2.2 局部屈曲分析方法

把帽緣條看作突緣，其兩邊可作一邊自由，一邊簡支處理；把帽腰看作腹板，兩邊可作簡支處理；把中間蒙皮看作腹板，兩邊可作簡支處理，如圖3所示。

圖3　分析所對(duì)應(yīng)的6個(gè)單元

2個(gè)帽緣條①和⑤，2個(gè)帽腰②和④，1個(gè)帽頂③，以及帽底中間蒙皮⑥，共6個(gè)單元屈曲載荷，將最小值作為長桁的局部屈曲載荷。

各單元的彎、扭剛度如表1所示。

表1　各單元彎、扭剛度

2.3 計(jì)算結(jié)果

當(dāng)α=80°時(shí)，局部屈曲載荷取得最小值，如圖4所示。在這種情況下，局部屈曲載荷隨變量W的變化曲線如圖4所示，顯然，當(dāng)α=80°，W=80 mm時(shí)，局部屈曲載荷取得最小值4 047 N，如圖4所示。

圖4　局部屈曲載荷隨W的變化曲線

給定α為50°～80°和W為60 mm～80 mm的條件下，計(jì)算出局部屈曲載荷隨W、α變化如圖5所示。

圖5　局部屈曲載荷隨W,α的變化曲線

3　壓損

3.1 壓損計(jì)算方法

長桁壓損載荷計(jì)算方法采用參考文獻(xiàn)[8]5.7.2.4節(jié)中介紹的分析方法，計(jì)算方法如下：

(3)

σcu=εcrExc

(4)

(6)

式中，b為腹板或者突緣的寬度；εcr為層壓板纖維方向的壓縮許用值；Exc為層壓板x方向的面內(nèi)剛度；Eyc為層壓板y方向的面內(nèi)剛度；μxy，μyx為層壓板面內(nèi)等效泊松比。

對(duì)長桁的組成單元逐一計(jì)算：

對(duì)于一邊自由、一邊簡支的組成單元：m=0.575，e=-0.797。

對(duì)于沒有自由端的組成單元：m=0.868，e=-0.869。

取各組成單元的σcc和σcu中的最低值，求其加權(quán)平均值得到壓損許用值：

(7)

壓損載荷Fcc計(jì)算公式為：

(8)

3.2 壓損分析

分析2個(gè)帽底、2個(gè)帽腰、1個(gè)帽頂和帽底中間蒙皮共6個(gè)組成單元的壓損許用值，再進(jìn)行加權(quán)平均得到整個(gè)長桁的壓損許用值，如圖2所示。該帽型長桁的材料力學(xué)參數(shù)見表2。

表2　材料力學(xué)參數(shù)

3.3 計(jì)算結(jié)果

壓損許用值σcc的計(jì)算結(jié)果如圖6所示。當(dāng)α=50°，w=61.62 mm時(shí)，壓損許用值取得最大值145.9 MPa；當(dāng)α=80°,w=80 mm時(shí)壓損許用值取得最小值120.77 MPa。

圖6　壓損許用值隨W，α的變化曲線

壓損載荷Fcc的計(jì)算結(jié)果如圖7所示。當(dāng)α=80°，W=80 mm時(shí)，壓損載荷取得最大值48.69 kN；當(dāng)α=50°，W=60 mm時(shí)，壓損載荷取得最小值46.68 kN。

圖7　壓損載荷隨W，α的變化曲線

4　結(jié)果對(duì)比

給定角度α，局部屈曲載荷和壓損許用值的計(jì)算結(jié)果如圖8、圖9所示。

圖8　壓損、局部屈曲計(jì)算結(jié)果對(duì)比(α=60°)

圖9　壓損、局部屈曲計(jì)算結(jié)果對(duì)比(α=70°)

從圖8中可以看出，當(dāng)α為60°時(shí)，影響長桁截面設(shè)計(jì)的載荷是局部屈曲載荷。局部屈曲載荷在W為60 mm～70 mm的區(qū)間內(nèi)變化微弱，曲線表現(xiàn)為平直；局部屈曲載荷在W為70 mm～80 mm的區(qū)間內(nèi)變化較為明顯，曲線有一定的斜率，且隨著W的增大，局部屈曲載荷變小。

從圖9中可以看出，當(dāng)α為70°時(shí)，影響長桁截面設(shè)計(jì)的載荷也是局部屈曲載荷。局部屈曲載荷在W為60 mm～80 mm的區(qū)間內(nèi)變化較為明顯，曲線有一定的斜率，且隨著W的增大，局部屈曲載荷變小。

4　結(jié)論

1) 影響帽型長桁截面設(shè)計(jì)的載荷是局部屈曲載荷；

2) 在局部屈曲狀況下，屈曲載荷隨帽型長桁底寬的增大而減??；

3) 在局部屈曲狀況下，屈曲載荷先隨帽腰和帽底的夾角α增大而增大，當(dāng)α為63°時(shí)屈曲載荷為最大，隨后α變大屈曲載荷變?。?/p>

4) 在壓損狀況下，壓損許用值隨帽型長桁底寬、帽腰和帽底的夾角α增大而減小。