裴少帥，胡迎春，高四良

?

基于MSC.Marc的固體發(fā)動機噴管非線性瞬態(tài)熱結(jié)構(gòu)耦合分析

裴少帥1，胡迎春2，高四良3

（1.西安航天動力技術研究所，陜西 西安 710025；2. 廣西師范大學職業(yè)技術師范學院，廣西 桂林 541004；3.西安航天動力機械廠，陜西 西安 710025）

利用MSC.Marc軟件對噴管進行了非線性瞬態(tài)熱結(jié)構(gòu)耦合分析。獲得了噴管結(jié)構(gòu)在工作過程中的溫度和應力分布、零部件粘接界面工作情況、預留熱膨脹間隙變化規(guī)律等。計算結(jié)果表明，依據(jù)常規(guī)設計方法設計的噴管方案結(jié)構(gòu)質(zhì)量有較大設計余量；喉襯與相鄰材料的膠粘面在噴管工作中很快失效；喉襯與擴張段絕熱層軸向端面對接間隙較小，在工作中會產(chǎn)生較大接觸應力；收斂段與擴張段螺紋聯(lián)接應力較小，有較高的設計可靠性。

噴管；熱結(jié)構(gòu)耦合；非線性瞬態(tài)分析

熱結(jié)構(gòu)設計是固體發(fā)動機結(jié)構(gòu)設計的重要部分，設計中需要考慮發(fā)動機結(jié)構(gòu)在內(nèi)部高溫高壓燃氣作用下的傳熱、燒蝕、碳化、強度、剛度等問題。設計中主要采用類比設計和試驗考核的傳統(tǒng)設計方法。基于成功經(jīng)驗的類比設計方法，對相近工況和相似結(jié)構(gòu)的設計，有一定的參考價值，但制約了新結(jié)構(gòu)的發(fā)展和應用；為了保證設計可靠性，設計裕度通常較大，不利于提高發(fā)動機的整體性能。試驗考核是對設計方案在特定載荷下的工作能力和可靠性的實物考核，是對設計方案評估的有效方法，但研制成本較高，研究周期較長。國外在固體發(fā)動機熱結(jié)構(gòu)設計過程中充分利于數(shù)值仿真技術，對發(fā)動機工作中的傳熱、燒蝕、碳化和結(jié)構(gòu)的強度、剛度等進行數(shù)值仿真分析計算，充分考慮熱與結(jié)構(gòu)的瞬態(tài)耦合作用。國內(nèi)也開展了相應的研究工作，對防熱材料在高溫下的力學及熱響應進行了討論，并應用于發(fā)動機結(jié)構(gòu)設計及性能分析[1-8]。相關研究結(jié)果表明，基于熱和結(jié)構(gòu)的完全耦合仿真分析技術的現(xiàn)代設計方法能夠設計計算出工作過程中每個時刻發(fā)動機結(jié)構(gòu)溫度、應力和剛度狀態(tài)，可以有效驗證新方案的可行性，改進設計方案，提高設計質(zhì)量。

本文在某發(fā)動機螺紋聯(lián)接式噴管設計中采用MSC.Marc通用有限元分析軟件熱結(jié)構(gòu)耦合分析模塊，對噴管方案進行瞬態(tài)熱結(jié)構(gòu)耦合分析計算，并依據(jù)計算結(jié)果對設計進行改進設計，優(yōu)化噴管結(jié)構(gòu)，顯著提高了噴管設計質(zhì)量。

1 螺紋聯(lián)接式噴管結(jié)構(gòu)

為了滿足設計空間、安裝工序和特殊維護等總體要求，發(fā)動機固定噴管與發(fā)動機殼體采用法蘭聯(lián)接。噴管擴張段與收斂段采用螺紋聯(lián)接，便于拆裝維護。噴管擴張段內(nèi)形面采用三次拋物線特殊形面，在保證擴張比的前提下最大限度減小擴張段長度。噴管結(jié)構(gòu)簡圖如圖1所示。噴管主要選材見表1。

擴張段與收斂段采用螺紋聯(lián)接，為了減小聯(lián)接螺紋結(jié)構(gòu)負荷，提高螺紋聯(lián)接結(jié)構(gòu)可靠性，必須減小擴張段重量；為了減小工作過程中擴張段外壁面溫升，必需增加噴管絕熱層厚度。擴張段絕熱層厚度的增加必將增加擴張段結(jié)構(gòu)重量，增大聯(lián)接螺紋結(jié)構(gòu)負荷。為了解決這一設計矛盾，在此對噴管結(jié)構(gòu)進行精確熱結(jié)構(gòu)耦合分析計算，合理設計擴張段絕熱層厚度，在保證噴管熱結(jié)構(gòu)可靠的前提下，減小螺紋聯(lián)接結(jié)構(gòu)載荷，提高聯(lián)接強度剛度。

圖1固體發(fā)動機噴管結(jié)構(gòu)簡圖

表1噴管主要材料表

名稱材料 喉襯高強石墨T705 喉襯背壁絕熱層高硅氧酚醛模壓 收斂段絕熱層碳纖維酚醛模壓 擴張段絕熱層碳布酚醛－高硅氧布酚醛復合纏繞 收斂段殼體高強鋼30CrMnSiA 擴張段殼體高強鋼30CrMnSiA

2 模型簡化及載荷分析

2.1 模型簡化

噴管計算模型中對噴管結(jié)構(gòu)進行了合理簡化。主要簡化如下：

（1）發(fā)動機噴管為軸對稱結(jié)構(gòu)，工作過程中高溫高壓燃氣在噴管內(nèi)軸對稱分布，因此計算采用軸對稱計算模型；

（2）收斂段和擴張段螺紋聯(lián)接結(jié)構(gòu)在圓周方向為非軸對稱結(jié)構(gòu)，但螺紋螺距為3 mm，螺紋直徑為232 mm，螺紋升角較小，簡化為軸對稱結(jié)構(gòu)對螺紋受力影響較小。噴管其他結(jié)構(gòu)和載荷軸對稱分布，因此計算模型采用軸對稱二維有限元模型。陰陽螺紋為幾何非線性接觸對；

（3）計算中不考慮絕熱層及喉襯的燒蝕；

（4）絕熱層在高溫燃氣作用下碳化分解等化學過程簡化為材料密度、比熱和熱導率隨溫度的非線性變化，不進行化學動力學分析計算；

（5）喉襯與喉襯背壁絕熱層、收斂段絕熱層、擴張段絕熱層等周邊非金屬材料采用膠粘工藝粘接組裝，工作過程中粘接界面上膠粘劑會分解碳化，膠粘面之間的連接關系從固結(jié)轉(zhuǎn)變?yōu)槟Σ两佑|。因此計算模型中喉襯周邊膠粘面處理為相互間可以傳遞熱量和力的多物理場非線性接觸。

2.2 載荷分析

噴管是固體發(fā)動機的能量轉(zhuǎn)換裝置，將固體推進劑藥柱燃燒后產(chǎn)生的高溫高壓燃氣，通過內(nèi)型面擴張、加速，并最終從擴張段大端口噴出，產(chǎn)生發(fā)動機前進的推力。高溫高壓燃氣在噴管中流動時，與噴管內(nèi)表面發(fā)生強烈的熱交換，向噴管內(nèi)表面?zhèn)鬟f熱量，沖刷、燒蝕噴管內(nèi)型面。作用噴管內(nèi)表面的燃氣壓力和燃氣與噴管內(nèi)表面的對流換熱是噴管工作過程中的主要熱結(jié)構(gòu)載荷。在此采用一維等熵定常流計算噴管內(nèi)燃氣壓力和溫度載荷，采用巴茲公式計算燃氣與噴管內(nèi)表面之間的對流換熱系數(shù)。計算公式如下：

式中：為溫度，K；為壓強，Pa，；為馬赫數(shù)；0為發(fā)動機燃燒室溫度，K；0為發(fā)動機平均壓強，Pa；為燃氣比熱比：為噴管內(nèi)計算截面處的通道面積，m2；A為噴管喉部截面處的通道面積，m2；h為對流換熱系數(shù)，w/（m2·K）；d為噴管喉徑，m；為燃氣動力粘性系數(shù)Pa·s；C為燃氣的定壓比熱容，J/（kg·K）；P為燃氣的普朗特數(shù)；為燃氣的特征速度，m/s；r為噴管喉部的曲率半徑，m；為對流換熱系數(shù)的修正因子，它隨馬赫數(shù)、燃氣的溫度以及壁面的溫度變化，取值在0.931和1.01之間，文中取1。

3 MSC.Marc熱結(jié)構(gòu)耦合計算模塊

MSC.Marc是功能齊全的高級非線性有限元分析求解器，具有極強的分析能力，可以處理各種線性和非線性結(jié)構(gòu)分析[9-10]。熱結(jié)構(gòu)耦合分析模塊具有求解各類熱邊界載荷和結(jié)構(gòu)邊界載荷作用下結(jié)構(gòu)的溫度、變形、應力、應變等結(jié)構(gòu)響應。噴管內(nèi)表面燃氣壓力、溫度、對流換熱系數(shù)沿噴管軸線非均勻分布，MSC.Marc提供了用子程序film和forcem分別用于施加復雜對流換熱系數(shù)和表面壓力[11]。

在此將一維等熵定常流計算噴管內(nèi)燃氣壓力、溫度載荷的計算公式和計算燃氣與噴管內(nèi)表面間對流換熱的巴爾茲公式程序化，嵌入film和forcem用戶子程序中。在瞬態(tài)熱結(jié)構(gòu)耦合計算時，MSC.Marc在每個時間迭代步初先調(diào)用用戶子程序計算噴管內(nèi)形面上的壓力、溫度和對流換熱系數(shù)載荷，再進行基于上步結(jié)果的該步傳熱和結(jié)構(gòu)應力計算。計算流程結(jié)構(gòu)框圖如圖2所示。

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

通過分析計算，獲得了噴管結(jié)構(gòu)在固體發(fā)動機燃氣作用下的溫度、應力、界面狀態(tài)、預留熱脹間隙等結(jié)構(gòu)載荷瞬態(tài)響應。對計算結(jié)果分析如下。

圖2計算流程結(jié)構(gòu)框圖

4.1 結(jié)構(gòu)溫度響應分析

圖3～圖6分別給出了噴管及主要零部組件溫度響應。由計算結(jié)果可以看出，發(fā)動機工作前5 s，噴管內(nèi)表面各部分溫度迅速上升，5 s時已基本達到最終溫度。工作結(jié)束時各零部組件結(jié)構(gòu)溫度達到最高值。其中收段絕熱層內(nèi)表面溫度最高為3208℃。在發(fā)動機工作5s時喉襯與收斂段絕熱層粘接界面溫度已達到450℃以上，粘接界面開始碳化分解，喉襯與收斂段絕熱層的粘接關系變?yōu)榻佑|關系。金屬殼體在發(fā)動機工作過程中溫升較小。

4.2 結(jié)構(gòu)應力分析

噴管結(jié)構(gòu)應力響應如圖7所示，最大應力發(fā)生在螺紋聯(lián)接處，最大應力為464.3 MPa。殼體材料高強鋼30CrMnSiA極限強度為1080 MPa，設計安全系數(shù)大于2.3。

4.3 預留熱脹間隙分析

為了避免工作中喉襯發(fā)生熱膨脹與擴張段絕熱層發(fā)生擠壓，在喉襯與擴張段絕熱層配合界面處預留熱脹間隙，喉襯與擴張段絕熱層對接面間隙隨工作時間的變化曲線如圖8所示。噴管工作3.7 s時，喉襯與擴張段絕熱層對接面間隙已變?yōu)榱?，喉襯與擴張段絕熱層將產(chǎn)生相互擠壓。擠壓使喉襯向噴管出口變形受阻，導致喉襯內(nèi)部應力、喉襯與背壁的粘接界面應力復雜化。同時使擴張段絕熱層入口產(chǎn)生嚴重破壞。

圖3 噴管溫度分布云圖

圖3噴管內(nèi)表面溫度隨時間變化曲線

圖4喉襯膠粘面上的溫度變化曲線

圖5金屬殼體膠粘面上的溫度變化曲線

圖6噴管結(jié)構(gòu)應力分布圖

圖7喉襯與擴張段對接面間隙變化曲線

5 總結(jié)

利用MSC.Marc通用有限元軟件對噴管結(jié)構(gòu)進行了非線性瞬態(tài)熱－結(jié)構(gòu)耦合分析計算，獲得了工作過程中噴管在高溫燃氣作用下熱結(jié)構(gòu)響應，對設計方案的評估提供了較準確的理論依據(jù)。計算結(jié)果表明依據(jù)常規(guī)設計方法設計的噴管方案有較大設計余量；喉襯與擴張段絕熱層端面對接間隙較小，會產(chǎn)生較大的相互擠壓。計算結(jié)果為噴管設計改進指明的方向，利于提高設計質(zhì)量。

[1]鄭曉亞，王衛(wèi)祥，邵明玉. 復合噴管傳力路徑及熱結(jié)構(gòu)分析[J]. 應用力學學報，2015，32（4）：618-623.

[2]白俊華，胡春波，李佳明. 固體發(fā)動機噴管喉襯流固耦合換熱規(guī)律數(shù)值研究[J]. 西北工業(yè)大學學報，2014，32（1）：23-28.

[3]韓珺禮，蒙上陽，楊軍輝，等. 固體火箭發(fā)動機噴管熱防護設計[J]. 計算機輔助工程，2013，22（a01）：160-163.

[4]鄭曉亞，王衛(wèi)祥，蔡飛超，等. 固體火箭發(fā)動機噴管熱結(jié)構(gòu)一體化計算[J]. 航空動力學報，2012，27（9）：2122-2127.

[5]孫林，閻海生，鮑福廷. 基于統(tǒng)計學的旋轉(zhuǎn)火箭發(fā)動機噴管熱結(jié)構(gòu)分析[J]. 固體火箭技術，2015，38（3）：356-371.

[6]龐越華，班大偉. 耦合場下固體火箭發(fā)動機噴管的動力學分析[J]. 機械工程師，2013（10）：77-78.

[7]孫林，鮑福廷，時偉. 旋轉(zhuǎn)固體火箭發(fā)動機噴管熱結(jié)構(gòu)分析[J]. 航空計算技術，2014，44（3）：62-65.

[8]張曉光，王長輝，劉宇，等. 固體火箭發(fā)動機喉襯流場及熱結(jié)構(gòu)耦合分析[J]. 固體火箭技術，2011，34（5）：579-583.

[9]王家昆，施慶華，李廣沈，等. 壓力鍋體有限元熱應力耦合分析[J]. 機械，2015，42（7）：36-39.

[10]陳真. 基于有限元的圓錐滾子軸承建模研究[J]. 2012（10）：21-24.

[11]陳火紅，尹偉奇，薜小香. MSC.Marc二次開發(fā)指南[M]. 北京：科學出版社，2004.

Nonlinear Transient Thermal-Structure Coupling Analysis of Solid Rocket Motor Nozzle Based on MSC.Marc

PEI Shaoshuai1，HU Yingchun2，GAO Siliang3

( 1.The Institute of Xi’an Aerospace Solid Propulsion Technology, Xi’an 710025, China; 2.Teachers College for Vocational and Technical Education, Guangxi Normal University, Guilin 541004, China; 3.Xi’an Aerospace Power Machinery Factory, Xi’an 710025, China )

The nonlinear transient thermal-structure coupling analysis of nozzle is carried out by using MSC.MARC software. The temperature and stress distribution of nozzle structure in the working process and the bonding interface of the parts and the variation law of the reserved thermal expansion gap are obtained. The calculation results show that the design of nozzle scheme has the large designing allowance according to the conventional design method. The adhesive surfaces of the throat and the adjacent material are quickly ineffective in the working process. The axial expansion gap of the throat and the expansion section is smaller and so that the greater contact stress. The thread joint between the convergent section and the expansion section has the small stress of and the high design reliability.

nozzle；coupling thermal structure；nonlinear transient analysis

TN751.1；V214.4+3

A

1006－0316 (2018) 04－0006－05

10.3969/j.issn.1006-0316.2018.04.002

2017－11－16

國家自然科學基金資助項目（51565007）

裴少帥（1978－），男，陜西西安人，碩士，高級工程師，主要研究方向為固體發(fā)動機設計；胡迎春（1971－），女，江蘇南京人，博士，教授，主要研究方向為算法與多學科優(yōu)化設計、計算機圖像處理等；高四良（1976－），男，陜西西安人，碩士，高級工程師，主要研究方向為機械系統(tǒng)設計。