發(fā)動(dòng)機(jī)導(dǎo)流裝置強(qiáng)度計(jì)算分析

湯禮軍

中航西安飛機(jī)工業(yè)集團(tuán)股份有限公司 陜西 西安 710089

引言

在飛機(jī)結(jié)構(gòu)設(shè)計(jì)中，一般情況下按照靜強(qiáng)度準(zhǔn)則已足以滿足使用要求。對(duì)于處于復(fù)雜環(huán)境中的結(jié)構(gòu)，內(nèi)外表面因存在巨大的溫度差，而產(chǎn)生很大的熱應(yīng)力，雖然結(jié)構(gòu)的常溫靜強(qiáng)度裕度系數(shù)較大，但常常會(huì)因?yàn)闊釕?yīng)力而導(dǎo)致結(jié)構(gòu)破壞。

發(fā)動(dòng)機(jī)導(dǎo)流裝置是用來(lái)保護(hù)機(jī)身蒙皮免受發(fā)動(dòng)機(jī)高溫氣流的熱蝕，承受飛機(jī)飛行時(shí)的氣動(dòng)載荷以及發(fā)動(dòng)機(jī)高溫氣流熱載荷，工作環(huán)境惡劣[1]。目前，對(duì)發(fā)動(dòng)機(jī)導(dǎo)流裝置的結(jié)構(gòu)強(qiáng)度分析，僅僅是考慮氣動(dòng)載荷的靜強(qiáng)度分析，尚未建立發(fā)動(dòng)機(jī)導(dǎo)流裝置的熱應(yīng)力有限元計(jì)算模型[2]。由于發(fā)動(dòng)機(jī)導(dǎo)流裝置存在大量不同材料的連接，受溫度影響將導(dǎo)致材料的熱變形和熱膨脹失配，進(jìn)而產(chǎn)生熱-應(yīng)力耦合可靠性問(wèn)題。因此發(fā)動(dòng)機(jī)導(dǎo)流裝置的熱應(yīng)力分析對(duì)結(jié)構(gòu)的強(qiáng)度校核至關(guān)重要。

為了獲取發(fā)動(dòng)機(jī)導(dǎo)流裝置在最大工作狀態(tài)下的溫度場(chǎng)分布和熱應(yīng)力狀況，應(yīng)用有限元的方法建立了有限元模型，根據(jù)發(fā)動(dòng)機(jī)導(dǎo)流裝置整體安裝結(jié)構(gòu)的溫度情況，確定了邊界條件，對(duì)發(fā)動(dòng)機(jī)導(dǎo)流裝置在最大工作狀態(tài)下的溫度場(chǎng)分布和熱應(yīng)力進(jìn)行了計(jì)算，同時(shí)與發(fā)動(dòng)機(jī)導(dǎo)流裝置在氣動(dòng)載荷下的應(yīng)力場(chǎng)分布進(jìn)行比較。對(duì)比結(jié)果為發(fā)動(dòng)機(jī)導(dǎo)流裝置結(jié)構(gòu)改進(jìn)設(shè)計(jì)提供了必要的依據(jù)。

1 有限元熱應(yīng)力計(jì)算描述

現(xiàn)有的各種結(jié)構(gòu)有限元素法分析軟件均有熱應(yīng)力計(jì)算功能，按指定的格式輸入必要的數(shù)據(jù)，即可執(zhí)行計(jì)算得到結(jié)果[3]。一般軟件提供兩種熱應(yīng)力計(jì)算途徑：一種方法為用同一個(gè)有限元模型作溫度計(jì)算和應(yīng)力分析，應(yīng)力分析時(shí)需要執(zhí)行兩個(gè)過(guò)程，先進(jìn)行熱分析，即按指定的熱條件和邊界條件計(jì)算結(jié)構(gòu)的節(jié)點(diǎn)溫度，然后將模型中各節(jié)點(diǎn)的溫度作為輸入文件，進(jìn)行有限元模型節(jié)點(diǎn)的位移和應(yīng)力的計(jì)算；另一種方法為結(jié)構(gòu)的熱分析和熱應(yīng)力計(jì)算用不同的模型或者不同的手段分別完成，應(yīng)力分析時(shí)由用戶提供有限元模型所有節(jié)點(diǎn)的溫度或者提供全部元素的平均溫度，由程序根據(jù)節(jié)點(diǎn)溫度或元素平均溫度計(jì)算等效熱載荷，最后完成應(yīng)力分析。

2 發(fā)動(dòng)機(jī)導(dǎo)流裝置有限元模型的建立

2.1 網(wǎng)格劃分

發(fā)動(dòng)機(jī)導(dǎo)流裝置的有限元建模及求解采用Patran程序完成，后處理采用Nastran程序完成[4]。網(wǎng)格劃分應(yīng)把有限元模型的節(jié)點(diǎn)、單元的分界線或分界面設(shè)置在相應(yīng)的不連續(xù)處，同時(shí)載荷的簡(jiǎn)化不能跨越主要受力構(gòu)件，網(wǎng)格過(guò)渡應(yīng)按單元尺度突然變化為最小的原則來(lái)處理。對(duì)于幾何形狀突變和過(guò)渡圓角應(yīng)采用較小的單元、較密集的網(wǎng)格，為了獲得較好的熱應(yīng)力結(jié)果，單元最大尺寸和最小尺寸之比不應(yīng)超過(guò)3。溫度梯度較小的方向，單元尺寸可以加大。

根據(jù)發(fā)動(dòng)機(jī)導(dǎo)流裝置真實(shí)結(jié)構(gòu)和單元類(lèi)型的力學(xué)特性，壁板蒙皮和隔框腹板簡(jiǎn)化為板殼單元，隔框的框緣簡(jiǎn)化為梁?jiǎn)卧诎彘L(zhǎng)桁簡(jiǎn)化為桿單元。為了減少計(jì)算過(guò)程中信息和數(shù)據(jù)的傳遞，方便將溫度場(chǎng)分析結(jié)果作為載荷輸入到熱應(yīng)力分析模型，溫度計(jì)算與應(yīng)力分析采用相同的有限元模型，有限元計(jì)算模型見(jiàn)圖1。

圖1 發(fā)動(dòng)機(jī)導(dǎo)流裝置有限元模型示意圖

2.2 邊界條件及載荷

發(fā)動(dòng)機(jī)導(dǎo)流裝置底部與機(jī)身蒙皮通過(guò)多排鉚釘連接，可以看作固支，約束3個(gè)平動(dòng)自由度和3個(gè)轉(zhuǎn)動(dòng)自由度；端部導(dǎo)流環(huán)與發(fā)動(dòng)機(jī)艙通過(guò)單排鉚釘連接，可以按簡(jiǎn)支處理，約束3個(gè)方向的平動(dòng)自由度。

發(fā)動(dòng)機(jī)導(dǎo)流裝置主要承受氣動(dòng)載荷和溫度載荷，應(yīng)力分析時(shí)分別對(duì)兩種載荷情況進(jìn)行計(jì)算。根據(jù)氣動(dòng)力分布計(jì)算結(jié)果，在發(fā)動(dòng)機(jī)導(dǎo)流裝置的蒙皮上施加0.01648MPa的均布?jí)毫Γ粶囟葓?chǎng)計(jì)算時(shí)，機(jī)身連接處的發(fā)動(dòng)機(jī)導(dǎo)流裝置蒙皮和端部導(dǎo)流環(huán)部分遠(yuǎn)離高溫氣流，溫度取40℃，中間蒙皮受高溫氣流輻射，溫度取120℃；參考溫度取25℃。

2.3 材料參數(shù)

表1 線膨脹系數(shù)

表2 導(dǎo)熱系數(shù)

3 計(jì)算結(jié)果與分析

熱應(yīng)力的計(jì)算時(shí)根據(jù)已有的溫度數(shù)據(jù),利用熱分析模塊對(duì)發(fā)動(dòng)機(jī)導(dǎo)流裝置的穩(wěn)態(tài)溫度場(chǎng)分布進(jìn)行計(jì)算[6]，將溫度場(chǎng)的計(jì)算結(jié)果建立場(chǎng)，換算成等效節(jié)點(diǎn)溫度載荷，導(dǎo)入到結(jié)構(gòu)分析模型中，作為熱應(yīng)力計(jì)算過(guò)程的載荷，即間接耦合的方式。

根據(jù)發(fā)動(dòng)機(jī)導(dǎo)流裝置溫度場(chǎng)計(jì)算結(jié)果,計(jì)算表明發(fā)動(dòng)機(jī)導(dǎo)流裝置結(jié)構(gòu)溫度分布梯度比較平緩，導(dǎo)流器中間蒙皮處溫度較高為120℃，導(dǎo)流器與機(jī)身連接處溫度較低為40℃。

考慮溫度場(chǎng)時(shí)，發(fā)動(dòng)機(jī)導(dǎo)流裝置計(jì)算結(jié)果見(jiàn)圖2。由圖2可知，發(fā)動(dòng)機(jī)導(dǎo)流裝置后部由于溫度場(chǎng)分布均勻和較好的邊界條件，應(yīng)力值在200MPa以下。發(fā)動(dòng)機(jī)導(dǎo)流裝置前端由于此區(qū)域溫度梯度大，導(dǎo)致熱變形不匹配，產(chǎn)生過(guò)大熱應(yīng)力，最大應(yīng)力為475MPa，出現(xiàn)在側(cè)壁板彎邊與端部導(dǎo)流環(huán)連接處，該材料在高溫下的強(qiáng)度極限為527MPa，其靜強(qiáng)度可以滿足設(shè)計(jì)要求，計(jì)算值與強(qiáng)度極限比較接近，若考慮高溫下的持久強(qiáng)度時(shí)，其強(qiáng)度極限可能會(huì)降低，該處結(jié)構(gòu)可能會(huì)發(fā)生破壞損傷。

考慮氣動(dòng)載荷時(shí)，發(fā)動(dòng)機(jī)導(dǎo)流裝置計(jì)算結(jié)果見(jiàn)圖3。由圖3可知，發(fā)動(dòng)機(jī)導(dǎo)流裝置應(yīng)力分布從整體上看是均勻合理的，結(jié)構(gòu)沒(méi)有出現(xiàn)局部應(yīng)力集中，最大應(yīng)力為192MPa，出現(xiàn)在側(cè)壁板底部靠近導(dǎo)流環(huán)區(qū)，該材料在常溫下的強(qiáng)度極限為569MPa，其靜強(qiáng)度可以滿足設(shè)計(jì)要求。

通過(guò)對(duì)圖2和圖3計(jì)算結(jié)果進(jìn)行對(duì)比，發(fā)動(dòng)機(jī)導(dǎo)流裝置熱應(yīng)力計(jì)算結(jié)果遠(yuǎn)遠(yuǎn)大于氣動(dòng)載荷作用下的計(jì)算結(jié)果，表明熱載荷是影響發(fā)動(dòng)機(jī)導(dǎo)流裝置應(yīng)力大小的主要因素。

圖2 發(fā)動(dòng)機(jī)導(dǎo)流裝置熱應(yīng)力計(jì)算云圖

圖3 發(fā)動(dòng)機(jī)導(dǎo)流裝置氣動(dòng)載荷計(jì)算云圖

4 結(jié)束語(yǔ)

本文采用有限元分析軟件，通過(guò)建立發(fā)動(dòng)機(jī)導(dǎo)流裝置有限元模型，調(diào)用熱分析及結(jié)構(gòu)分析模塊，對(duì)發(fā)動(dòng)機(jī)導(dǎo)流裝置結(jié)構(gòu)進(jìn)行了溫度場(chǎng)分布和熱應(yīng)力計(jì)算，得到以下主要結(jié)論：導(dǎo)流器溫度場(chǎng)呈現(xiàn)高溫邊界到低溫邊界的梯度分布；對(duì)兩種載荷工況下發(fā)動(dòng)機(jī)導(dǎo)流裝置應(yīng)力進(jìn)行計(jì)算，根據(jù)材料特性可知是滿足靜強(qiáng)度設(shè)計(jì)要求的，若考慮高溫下的持久強(qiáng)度時(shí)，側(cè)壁板彎邊與端部導(dǎo)流環(huán)連接處結(jié)構(gòu)可能會(huì)發(fā)生破壞損傷；比較兩種工況下應(yīng)力計(jì)算結(jié)果得出熱載荷時(shí)影響發(fā)動(dòng)機(jī)導(dǎo)流裝置應(yīng)力大小的主要因素。文章的計(jì)算結(jié)果可以為發(fā)動(dòng)機(jī)導(dǎo)流裝置的熱防護(hù)以及優(yōu)化設(shè)計(jì)提供可靠的依據(jù)，同時(shí)此計(jì)算分析方法為其他熱環(huán)境下工作的結(jié)構(gòu)強(qiáng)度校核提供了一定的工程參考價(jià)值。