潛入式可拋噴管數(shù)值仿真范文

本站小編為你精心準備了潛入式可拋噴管數(shù)值仿真參考范文，愿這些范文能點燃您思維的火花，激發(fā)您的寫作靈感。歡迎深入閱讀并收藏。

《彈箭與制導(dǎo)學(xué)報》2016年第一期

摘要：

對一種潛入式可拋噴管進行強度分析和流場仿真，強度仿真結(jié)果揭示了燃氣壓強作用下可拋噴管的應(yīng)力分布，流場仿真結(jié)果揭示了燃氣經(jīng)由可拋噴管排出時的流場特征，結(jié)果表明可拋噴管整體結(jié)構(gòu)設(shè)計能夠滿足要求。研究中模擬了環(huán)形切割器工作對沖壓噴管被剪切部位的作用和影響，對模擬仿真結(jié)果與環(huán)形切割器功能驗證剪切試驗結(jié)果進行了對比分析，認為殼體剪切部位須設(shè)計一定厚度的臺階。數(shù)值仿真結(jié)果有助于可拋噴管的結(jié)構(gòu)優(yōu)化設(shè)計。

關(guān)鍵字：

可拋噴管；環(huán)形切割器；數(shù)值仿真；優(yōu)化設(shè)計

本文提出一種潛入式可拋噴管方案。該方案中可拋噴管通過連接件緊固連接于沖壓噴管之上，在助推器工作時穩(wěn)定工作，發(fā)動機轉(zhuǎn)級時依靠級間分離機構(gòu)可靠拋除。沖壓噴管為助推補燃室的一部分，可拋噴管拋除后，沖壓噴管滿足發(fā)動機沖壓級的工作要求。對潛入式可拋噴管進行了強度分析和流場仿真，強度分析和流場仿真的結(jié)果能夠為工程設(shè)計提供參考和依據(jù)。模擬分析了環(huán)形切割器的工作對剪切部位的作用和影響，并對比試驗結(jié)果進行數(shù)據(jù)分析，為可拋噴管結(jié)構(gòu)優(yōu)化提供依據(jù)。

1可拋噴管強度仿真

1.1可拋噴管模型潛入式可拋噴管的物理模型如圖1所示。模型包含金屬殼體、絕熱層、喉襯等部件。

1.2計算網(wǎng)格噴管模型為軸對稱結(jié)構(gòu)，對其進行網(wǎng)格劃分，網(wǎng)格數(shù)量為40萬，具體的網(wǎng)格劃分如圖2所示：

1.3材料屬性定義噴管模型由三種材料構(gòu)成，其屬性參數(shù)如表1所示。金屬殼體的材料為30CrMnSiA，絕熱層的材料為高硅氧/酚醛，沖壓噴管和可拋噴管喉襯的材料為石墨。

1.4載荷和邊界條件強度計算中，載荷為內(nèi)壓載荷，載荷的設(shè)置依據(jù)實際工況中內(nèi)壓的極值，在此基礎(chǔ)上給定足夠的安全系數(shù)f=1.5，噴管擴張段給定相對較小的壓強，外部壓強設(shè)為零，這樣的設(shè)置可以保證計算結(jié)果的相對保守。邊界條件的設(shè)置包括兩個部分：1）過軸線截面的對稱；2）模型前端面的固支載荷。

1.5計算結(jié)果分析載荷施加到模型上，模型就會產(chǎn)生變形，當載荷超過一定的限度時，模型就會產(chǎn)生過度變形或破壞，與變形和破壞相對應(yīng)的就是材料的剛度和強度。為了保證機構(gòu)正常工作，要對模型的強度和剛度提出如下的要求：1）強度要求，在規(guī)定載荷的作用下，結(jié)構(gòu)最大應(yīng)力應(yīng)當不大于許用應(yīng)力，即結(jié)構(gòu)應(yīng)當有足夠的抵抗破壞的能力；2）剛度要求，在載荷的作用下，結(jié)構(gòu)即使有足夠的強度，但是如果變形過大，仍然不能正常工作，即機構(gòu)應(yīng)當有足夠的抵抗變形的能力。這里的強度仿真主要考查模型的強度和剛度。殼體上的最大應(yīng)力值為710.7MPa（端面固支載荷處應(yīng)力值不計在內(nèi)），安全系數(shù)f=1.5已經(jīng)體現(xiàn)在內(nèi)壓載荷上，所以殼體的設(shè)計應(yīng)力σsj=710.7MPa。對于殼體材料30CrMnSiA，其抗拉強度極限為1080MPa，殼體計算富裕系數(shù)n=1.52。絕熱層上的最大壓應(yīng)力值為87.9MPa。同樣的，由于安全系數(shù)f=1.5已經(jīng)體現(xiàn)在了內(nèi)壓載荷上，即有絕熱層設(shè)計應(yīng)力σsj=87.9MPa。對于絕熱層材料高硅氧/酚醛，其抗壓強度極限為200MPa，絕熱層計算富裕系數(shù)n=2.28。石墨喉襯上最大應(yīng)力值為48.9MPa，同樣的，安全系數(shù)f=1.5已經(jīng)體現(xiàn)在內(nèi)壓載荷上，即有喉襯設(shè)計應(yīng)力σsj=48.9MPa。對于內(nèi)襯材料石墨，石墨的強度極限為70MPa，石墨喉襯的計算富裕系數(shù)n=1.43。由U1（軸向）、U2（徑向）兩個方向的總體位移圖，可知在X、Y兩個方向的總體位移的最大值均不大于1mm，整個機構(gòu)上沒有出現(xiàn)過大的變形，機構(gòu)能夠正常工作。最為關(guān)注的密封面上，相配合的兩個面幾乎不存在相對位移，密封圈能夠起到較好的密封效果。

2可拋噴管流場仿真

2.1計算網(wǎng)格對噴管模型進行流場計算時，我們主要關(guān)注燃氣的流道，而不過多關(guān)心噴管殼體的外部構(gòu)型。可拋噴管與沖壓噴管之間形成了一個背腔，背腔內(nèi)部的流動情況也是關(guān)注的一個方面。由于噴管模型的對稱性，取截面的一半作為二維模型，這樣結(jié)構(gòu)更加簡單，可以在同樣的網(wǎng)格密度的情況下減少網(wǎng)格的數(shù)量，既不影響計算結(jié)果的正確性又能降低計算成本。網(wǎng)格的劃分如圖9所示，采用結(jié)構(gòu)網(wǎng)格，近壁處以及速度梯度變化大的區(qū)域的網(wǎng)格進行局部加密處理。

2.2求解器和邊界條件本次計算中選擇密度基耦合求解器隱式方式，以矢量方程求解連續(xù)性方程、動量方程、能量方程和組分方程，通過狀態(tài)方程得到壓強，其他標量方程按照分離方式求解。由于網(wǎng)格為二維結(jié)構(gòu)化網(wǎng)格，對流項的差值方法采用一階精度既能保證計算的精度又易于收斂，梯度選項選擇Green-GaussCellBased，殘差收斂條件設(shè)置為10-5，湍流模型選擇Standardk-e模型。

2.3流場仿真及結(jié)果分析1）本次計算中，入口邊界壓強設(shè)置為9MPa，即有總壓P0=9MPa，噴管出口處的總壓均值約為8.05MPa，總壓恢復(fù)系數(shù)約為0.895。2）本次計算中，入口邊界溫度值設(shè)為3200K，由計算結(jié)果知，出口平均溫度約為1500K。3）由計算結(jié)果可知，噴管出口速度最大值約為2480m/s，平均速度約為2300m/s。4）由背腔內(nèi)部的速度矢量分布圖可知，所關(guān)心的背腔內(nèi)部渦流并沒有想象中的那么劇烈，背腔內(nèi)部的燃氣速度基本都低于5m/s。

3可拋噴管剪切試驗?zāi)M仿真

3.1可拋噴管剪切部位模型潛入式可拋噴管通過銷釘與助推補燃室沖壓噴管連接，剪切部位模型如圖11所示，通過仿真軟件仿真計算前對模型進行簡化，如圖12所示。通過調(diào)整保留部分和剪切部分的殼體厚度來改變可拋噴管剪切部位及其相鄰部位的狀態(tài)，針對不同狀態(tài)的殼體進行仿真計算，對比計算結(jié)果與試驗結(jié)果，為剪切部位結(jié)構(gòu)優(yōu)化提供依據(jù)。

3.2環(huán)形切割器剪切行為模擬模型中剪切部分厚度均為2.5mm，保留部分與剪切部分間設(shè)有臺階，臺階厚度為0mm-10mm不等，通過仿真軟件模擬環(huán)形切割器對不同狀態(tài)殼體的作用和剪切部位的影響。對比分析圖13-圖18認為，當殼體剪切部分厚度為2.5mm，臺階的厚度越厚，則環(huán)形切割器工作時對殼體保留部分的影響越小，保留部分殼體沿徑向的變形越小。

3.3環(huán)形切割器剪切試驗環(huán)形切割器剪切試驗共完成了4發(fā)，均為常溫狀態(tài)。試驗中工裝殼體切割部分厚度為2.5mm、3mm和3.5mm，臺階厚度均為9mm。

3.4剪切試驗?zāi)M仿真結(jié)果分析對比分析試驗結(jié)果和模擬仿真結(jié)果，認為試驗殼體臺階厚度越大，環(huán)形切割器工作對殼體造成的變形越小。由于可拋噴管保留部分殼體上模壓粘接著高硅氧/酚醛復(fù)合絕熱材料，環(huán)形切割器工作對殼體造成的變形過大時，會導(dǎo)致復(fù)合絕熱材料與保留部分殼體脫粘，影響沖壓發(fā)動機的正常工作。四次剪切試驗中，殼體臺階厚度為9mm時，環(huán)形切割器工作時保留部分殼體變形量非常微小，可滿足設(shè)計要求。參考不同厚度的殼體臺階在模擬仿真中變形量（如表3所示），分析認為臺階厚度設(shè)計為6mm-10mm時，殼體變形量可滿足設(shè)計要求，這一結(jié)論為可拋噴管結(jié)構(gòu)優(yōu)化提供了依據(jù)。

4結(jié)論

本文通過對所設(shè)計的潛入式可拋噴管進行仿真分析，得出如下的結(jié)論：1）強度仿真結(jié)果揭示了可拋噴管的應(yīng)力分布，結(jié)果表明噴管本體及其連接部位均滿足結(jié)構(gòu)強度設(shè)計要求，在設(shè)置足夠的安全系數(shù)的情況下，仍具有一定的計算富裕系數(shù)。2）由流場仿真結(jié)果可知，噴管型面結(jié)構(gòu)能夠滿足性能要求。潛入式噴管的背腔內(nèi)部存在渦流，渦流速度較小，對噴管性能影響不大。3）噴管剪切位置附近臺階厚度設(shè)計為6mm-10mm時，殼體變形量可滿足設(shè)計要求，可拋噴管優(yōu)化設(shè)計時臺階厚度須大于6mm。

參考文獻：

[1]傅德彬.數(shù)值仿真及其在航天發(fā)射技術(shù)中的應(yīng)用[M].北京：國防工業(yè)出版社，2011.

[2]王定軍，宋會玲，劉偉.冷氣推進器流場仿真與分析[J].火箭推進，2010，36（2）：36~39.

[3]張碩，王寧飛.固體火箭發(fā)動機噴管及出口處流場特性的數(shù)值分析[J].彈箭與制導(dǎo)學(xué)報，2007，27（1）：177~180.

[4]宋大明，周長省.固體火箭發(fā)動機瞬態(tài)內(nèi)流暢數(shù)值仿真[J].彈箭與制導(dǎo)學(xué)報，2010，30（6）.

[5]劉敦啟，張澤遠.固體火箭發(fā)動機長尾噴管內(nèi)襯燒蝕流場分析[J].彈箭與制導(dǎo)學(xué)報，2011，31（2）：27~29.

[6]曹軍，房雷，吳川.某固體火箭沖壓發(fā)動機絕熱層傳熱燒蝕數(shù)值仿真研究[J].航空兵器，2011，（3）.

[7]孫得川，由旭.發(fā)動機引流推力矢量方案的內(nèi)流暢分析[J].航空兵器，2014，（5）.

[8]王定軍，宋會玲，白少卿，魏超.減壓閥節(jié)流口流場仿真和分析[J].火箭推進，2009，35（6）：37~4.

作者：商慧增 崔金平 單位：中國空空導(dǎo)彈研究院