專利名稱：一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法
一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法技術(shù)領(lǐng)域
本發(fā)明屬于固體推進劑煙霧特性評價技術(shù)領(lǐng)域，尤其是涉及一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法。
背景技術(shù)：
隨著機械化、信息化技術(shù)的發(fā)展，現(xiàn)代武器系統(tǒng)越來越多的采用雷達、光電等自動跟蹤、控制和制導(dǎo)手段，武器系統(tǒng)的自動化程度、打擊精度等方面得到顯著提高，極大的提升了彈箭武器系統(tǒng)的作戰(zhàn)效能。但與此同時，彈箭武器發(fā)射和運載過程中的煙焰問題危害性日益顯現(xiàn)，對武器系統(tǒng)信息化、智能化、隱身化性能影響很大，對武器系統(tǒng)快速、精確打擊、隱身等功能造成了嚴重的不利影響，甚至已成為制約我軍彈箭武器發(fā)揮正常效能、新裝備訂貨以及新型武器裝備發(fā)展的主要技術(shù)瓶頸。尤其是近年來，重型反坦克導(dǎo)彈研制過程中出現(xiàn)推進劑裝藥煙霧遮擋制導(dǎo)信號的問題，解決問題耗費巨大。
固體推進劑是彈箭武器系統(tǒng)的發(fā)射能源，其低特征信號性能是保證武器系統(tǒng)戰(zhàn)場使命完成的重要因素?；鸺艢庥鹆鞯奶卣餍盘?Exhaust Plume Signature)是一種包含有系統(tǒng)或火箭發(fā)動機排氣全部性能或特性的術(shù)語，此性能或特性可被用作探測、識別、或攔截執(zhí)行任務(wù)的發(fā)射平臺或?qū)?。羽流的特征信號特性主要包括煙、輻射能的散發(fā)、能見度(視程/能見距離)和雷達波吸收等4項參數(shù)。固體推進劑煙霧是羽流特征信號的重要表征參數(shù)，其從表觀特征可以細分為:一次煙霧和二次煙霧，一次煙霧主要由火箭發(fā)動機噴管噴出的凝相顆粒構(gòu)成，主要存在形式為金屬氫氧化物、氧化物和氯化物等，復(fù)合推進劑中的金屬粉(鋁、鎂等)燃燒產(chǎn)生的金屬氧化物(Al203、Mg0等)是一次煙霧的主要來源；二次煙霧由燃氣中可凝氣體組成，是燃氣與大氣相互作用的結(jié)果。復(fù)合推進劑中的AP和含氟組分燃燒分解后形成HCl、HF，在特定的溫度和濕度條件下，與空氣相結(jié)合形成水和HC1、HF的共沸液滴煙云。
二次煙作為特征信號性能的主要組成部分，是排氣羽流與大氣相互作用的結(jié)果，其表觀特征是彈箭武器在噴管明亮的后燃區(qū)后部拖有長長得煙云尾跡。為評價固體推進劑二次煙霧特性，世界上一些發(fā)達國家都建立了相關(guān)測試評估方法，建設(shè)了各種實驗測試設(shè)施來對固體推進劑特征信號性能進行檢測及表征。尤其是近年來CRAFT Tech公司獲得了美國國防部SBIR項目資助開展的二次煙形成模型的研究。國內(nèi)在對國外測試設(shè)施研究的基礎(chǔ)上，已建立起一套具有獨立知識產(chǎn)權(quán)的低成本、易控制的固體推進劑特征信號檢測系統(tǒng)，起草完成了相關(guān)行業(yè)標準，發(fā)表了系列的研究論文，申請了一些國防專利。但在低特征信號推進劑研制及評價過程中，解決問題需要耗費大量的人力物力，研制周期大幅度延長，帶來巨大的經(jīng)濟損失和政治影響，因而現(xiàn)如今迫切需要一種虛擬試驗方法對固體推進劑的二次煙特性進行快速評價。虛擬試驗是借助于計算機的高速解算，按實際試驗要求對基于推進劑燃燒和燃燒產(chǎn)物流動的數(shù)學(xué)模型進行模擬試驗，不僅可以作為真實試驗的前期準備工作或在一定程度上替代傳統(tǒng)的試驗(如一些極限工況)；并且能大幅減少真實試驗次數(shù)，降低試驗費用，縮短試驗周期；同時，具有較好的交互性，使得各種試驗信息及時反饋；另外，虛擬試驗不受氣象條件、場地、時間和次數(shù)的限制，試驗過程可方便實現(xiàn)回放、再現(xiàn)和重復(fù)。發(fā)明內(nèi)容
本發(fā)明所要解決的技術(shù)問題在于針對上述現(xiàn)有技術(shù)中的不足，提供一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，其方法步驟簡單、設(shè)計合理且實現(xiàn)方便、使用效果好，能有效解決現(xiàn)有固體推進劑二次煙信號性能評價時存在的人力物力投入大、周期長、試驗量大等缺陷。
為解決上述技術(shù)問題，本發(fā)明采用的技術(shù)方案是:一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，其特征在于該方法包括以下步驟:
步驟一、初始參數(shù)設(shè)定與存儲:
首先，通過與數(shù)據(jù)處理器相接的參數(shù)輸入單元，輸入配制所評價固體推進劑所用的組分數(shù)量N和各組分的組分信息，并將所輸入的信息同步存儲至與所述數(shù)據(jù)處理器相接的數(shù)據(jù)存儲單元內(nèi)；其中，各組分的組分信息均包括化學(xué)式和質(zhì)量配比IV i=l、2、…、N，N為配制所評價固體推進劑所用組分的數(shù)量；其中，0 < Hii < 100，N ^ 2 ；
之后，通過所述參數(shù)輸入單元自預(yù)先建立的燃燒產(chǎn)物數(shù)據(jù)庫中，選出所評價固體推進劑燃燒后所產(chǎn)生的所有燃燒產(chǎn)物；所述燃燒產(chǎn)物數(shù)據(jù)庫中存儲有多種燃燒產(chǎn)物的屬性信息；其中，各燃燒產(chǎn)物的屬性信息均包括化學(xué)式、相對分子質(zhì)量和相態(tài)，其中相態(tài)為氣相或凝聚相；同時，通過所述參數(shù)輸入單元對所評價固體推進劑燃燒后所產(chǎn)生燃燒產(chǎn)物的數(shù)量m和所產(chǎn)生燃燒產(chǎn)物中氣相產(chǎn)物的數(shù)量Q進行設(shè)定；
步驟二、能量特性參數(shù)計算，其計算過程如下:
步驟201、建模:根據(jù)最小自由能原理，建立最小自由能數(shù)學(xué)模型和燃燒室溫度計算模型；
步驟202、平衡組成計算:所述數(shù)據(jù)處理器調(diào)用步驟201中所建立的最小自由能數(shù)學(xué)模型，并結(jié)合步驟一中所設(shè)定的初始參數(shù)，計算得出所評價固體推進劑燃燒后的燃燒產(chǎn)物平衡組成；
步驟203、絕熱燃燒溫度計算:所述數(shù)據(jù)處理器調(diào)用步驟201中所建立的燃燒室溫度計算模型，計算得出所評價固體推進劑燃燒后處于化學(xué)平衡狀態(tài)時的絕熱燃燒溫度；
步驟三、羽流場計算，其計算過程如下:
步驟301、燃燒室工作參數(shù)以及噴管幾何參數(shù)與噴流計算域范圍設(shè)定:首先，通過所述參數(shù)輸入單元對發(fā)動機`元對發(fā)動機噴管的幾何參數(shù)和噴流計算域范圍進行設(shè)定；其中，發(fā)動機噴管的幾何參數(shù)包括入口半徑^、喉部半徑r2、出口半徑r3、收斂段長度dl、喉部圓柱段長度d2、擴張角a和擴張段長度d3，其中入口半徑ri為發(fā)動機噴管的入口半徑，喉部半徑r2為發(fā)動機噴管的喉部半徑，出口半徑r3為發(fā)動機噴管的出口半徑，收斂段長度dl為發(fā)動機噴管入口到噴管喉部前端的長度，喉部圓柱段長度d2為發(fā)動機噴管的喉部長度，擴張角a為發(fā)動機噴管擴張段的壁面與其軸線之間的夾角，且擴張段長度d3為發(fā)動機噴管喉部末端到噴管出口的長度；噴流計算域范圍包括噴流軸向長度Xm和噴流徑向長度ym，其中噴流軸向長度xm為發(fā)動機噴管出口到噴流下游出口的長度，噴流徑向長度ym為發(fā)動機噴管軸線到噴流徑向外邊界的長度；
步驟302、發(fā)動機羽流場計算:所述數(shù)據(jù)處理器調(diào)用FULENT軟件，并結(jié)合步驟二中計算得出的能量特性參數(shù)、步驟301中所設(shè)定的燃燒室工作參數(shù)以及噴管幾何參數(shù)與噴流計算域范圍，對所評價固體推進劑進行羽流場計算，并相應(yīng)計算得出所評價固體推進劑燃燒后所產(chǎn)生的各氣相產(chǎn)物在發(fā)動機噴管軸線上各位置處的摩爾百分含量fz和發(fā)動機噴管軸線上各位置處的壓強P氣相；其中，z為正整數(shù)，且z=l、-,Q；
步驟四、二次煙特征值計算，其計算過程如下:
步驟401、建模:采用所述數(shù)據(jù)處理器，建立固體推進劑的二次煙特征值計算模型，且所建立的二次煙特征值計算模型為二次煙相對濕度計算模型:RHamb=IOOX [K-ftotalX0.16589)] (I)，式⑴中/toM; =Zffj0 + fHC, +/冊，其中 ftQtal 為所評價固體推進劑燃燒后所產(chǎn)生燃燒產(chǎn)物中凝聚物在發(fā)動機噴管軸線上的摩爾百分含量，為所評價固體推進劑燃燒后所產(chǎn)生燃燒產(chǎn)物中的水在發(fā)動機噴管軸線上的摩爾百分含量，fHci為所評價固體推進劑燃燒后所產(chǎn)生燃燒產(chǎn)物中的HCl在發(fā)動機噴管軸線上的摩爾百分含量，且fHF為所評價固體推進劑燃燒后所產(chǎn)生燃燒產(chǎn)物中的HF在發(fā)動機噴管軸線上的摩爾百分含量；，式中Ptl為通過所述參數(shù)輸入單元預(yù)先輸入的0°C時水的飽和蒸氣壓，P=PHa+PHF，其中Prel為氣相產(chǎn)物HCl的氣體分壓且PHa=fHaXP，，Phf為氣相產(chǎn)物HF的氣體分壓且PHF=fHFXPn#，為步驟302中計算得出的發(fā)動機噴管軸線上的壓強；
步驟402、二次煙特征值計算:所述數(shù)據(jù)處理器調(diào)用步驟401中所建立的二次煙相對濕度計算模型，并結(jié)合步驟302中計算得出的所評價固體推進劑燃燒后所產(chǎn)生的各氣相產(chǎn)物在發(fā)動機噴管軸線上各位置處的摩爾百分含量fz，計算得出所評價固體推進劑燃燒后的發(fā)動機羽流場中噴管軸線上各位置處的二次煙相對濕度RH-。
上述一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，其特征是:步驟402計算得出所評價固體推進劑燃燒后的發(fā)動機羽流場中噴管軸線上各位置處的二次煙相對濕度RHamb后，所述數(shù)據(jù)處理器還需調(diào)用二次煙特征值評價模塊對發(fā)動機羽流場中噴管軸線上各位置處的二次煙等級進行評價，且RHamb值越大，二次煙等級越高。
上述一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，其特征是:步驟402計算得出所評價固體推進劑燃燒后的發(fā)動機羽流場中噴管軸線上各位置處的二次煙相對濕度RHamb后，所述數(shù)據(jù)處理器還需調(diào)用曲線繪制模塊繪制出所評價固體推進劑燃燒后所產(chǎn)生二次煙的相對濕度RHamb的軸向分布曲線，并通過所述顯示單元對所繪制的相對濕度RHamb的軸向分布曲線進行同步顯示；其中，相對濕度RHamb的軸向分布曲線為所評價固體推進劑燃燒后所產(chǎn)生二次煙的相對濕度RHamb隨X軸的取值變化而變化的曲線；其中X軸為發(fā)動機噴管軸線所在的坐標軸，且發(fā)動機噴管出口處x=0。
上述一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，其特征是:步驟302中計算得出的所評價固體推進劑燃燒后所產(chǎn)生的各氣相產(chǎn)物在發(fā)動機噴管軸線上各位置處的摩爾百分含量匕后，所述數(shù)據(jù)處理器還需曲線繪制模塊繪制出所評價固體推進劑燃燒后所產(chǎn)生各氣相產(chǎn)物在發(fā)動機噴管軸線上的摩爾百分含量分布曲線Iz ;其中，摩爾百分含量分布曲線Iz為所評價固體推進劑燃燒后所產(chǎn)生的第z種氣相產(chǎn)物的摩爾百分含量fz隨X軸的取值變化而變化的曲線；其中X軸為發(fā)動機噴管軸線所在的坐標軸，且發(fā)動機噴管出口處x=0。
上述一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，其特征是:步驟401 中 p0=610.78Pa。
上述一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，其特征是:所述數(shù)據(jù)處理器調(diào)用所述二次煙特征值評價模塊對發(fā)動機羽流場中噴管軸線上各位置處的二次煙等級進行評價時，當RHamb值在區(qū)間(0.90，I]時，二次煙等級為微煙；當1^_值在區(qū)間(0.55，0.90]時，二次煙等級為少煙；且當RHamb值在區(qū)間
時，二次煙等級為有煙。
上述一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，其特征是:步驟201中所建立的最小自由能數(shù)學(xué)模型為
權(quán)利要求
1.一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，其特征在于該方法包括以下步驟: 步驟一、初始參數(shù)設(shè)定與存儲: 首先，通過與數(shù)據(jù)處理器相接的參數(shù)輸入單元，輸入配制所評價固體推進劑所用的組分數(shù)量N和各組分的組分信息，并將所輸入的信息同步存儲至與所述數(shù)據(jù)處理器相接的數(shù)據(jù)存儲單元內(nèi)；其中，各組分的組分信息均包括化學(xué)式和質(zhì)量配比IV i=l、2、…、N，N為配制所評價固體推進劑所用組分的數(shù)量；其中，O < Hii < 100，N彡2 ; 之后，通過所述參數(shù)輸入單元自預(yù)先建立的燃燒產(chǎn)物數(shù)據(jù)庫中，選出所評價固體推進劑燃燒后所產(chǎn)生的所有燃燒產(chǎn)物；所述燃燒產(chǎn)物數(shù)據(jù)庫中存儲有多種燃燒產(chǎn)物的屬性信息；其中，各燃燒產(chǎn)物的屬性信息均包括化學(xué)式、相對分子質(zhì)量和相態(tài)，其中相態(tài)為氣相或凝聚相；同時，通過所述參數(shù)輸入單元對所評價固體推進劑燃燒后所產(chǎn)生燃燒產(chǎn)物的數(shù)量m和所產(chǎn)生燃燒產(chǎn)物中氣相產(chǎn)物的數(shù)量Q進行設(shè)定； 步驟二、能量特性參數(shù)計算，其計算過程如下: 步驟201、建模:根據(jù)最小自由能原理，建立最小自由能數(shù)學(xué)模型和燃燒室溫度計算模型； 步驟202、平衡組成計算:所述數(shù)據(jù)處理器調(diào)用步驟201中所建立的最小自由能數(shù)學(xué)模型，并結(jié)合步驟一中所設(shè)定的初始參數(shù)，計算得出所評價固體推進劑燃燒后的燃燒產(chǎn)物平衡組成； 步驟203、絕熱燃燒溫度計算:所述數(shù)據(jù)處理器調(diào)用步驟201中所建立的燃燒室溫度計算模型，計算得出所評價固體推進劑燃燒后處于化學(xué)平衡狀態(tài)時的絕熱燃燒溫度； 步驟三、羽流場計算，其計算過程如下: 步驟301、燃燒室工作參數(shù)以及噴管幾何參數(shù)與噴流計算域范圍設(shè)定:首先，通過所述參數(shù)輸入單元對發(fā)動機的燃燒室內(nèi)壓強P。、環(huán)境壓力和環(huán)境溫度Tjf分別進行設(shè)定；之后，再通過所述參數(shù)輸入單元對發(fā)動機噴管的幾何參數(shù)和噴流計算域范圍進行設(shè)定；其中，發(fā)動機噴管的幾何參數(shù)包括入口半徑^、喉部半徑r2、出口半徑r3、收斂段長度dl、喉部圓柱段長度d2、擴張角a和擴張段長度d3，其中入口半徑ri為發(fā)動機噴管的入口半徑，喉部半徑1~2為發(fā)動機噴管的喉部半徑，出口半徑r3為發(fā)動機噴管的出口半徑，收斂段長度dl為發(fā)動機噴管入口到噴管喉部前端的長度，喉部圓柱段長度d2為發(fā)動機噴管的喉部長度，擴張角a為發(fā)動機噴管擴張段的壁面與其軸線之間的夾角，且擴張段長度d3為發(fā)動機噴管喉部末端到噴管出口的長度；噴流計算域范圍包括噴流軸向長度Xm和噴流徑向長度ym，其中噴流軸向長度xm為發(fā)動機噴管出口到噴流下游出口的長度，噴流徑向長度ym為發(fā)動機噴管軸線到噴流徑向外邊界的長度； 步驟302、發(fā)動機羽流場計算:所述數(shù)據(jù)處理器調(diào)用FULENT軟件，并結(jié)合步驟二中計算得出的能量特性參數(shù)、步驟301中所設(shè)定的燃燒室工作參數(shù)以及噴管幾何參數(shù)與噴流計算域范圍，對所評價固體推進劑進行羽流場計算，并相應(yīng)計算得出所評價固體推進劑燃燒后所產(chǎn)生的各氣相產(chǎn)物在 發(fā)動機噴管軸線上各位置處的摩爾百分含量fz和發(fā)動機噴管軸線上各位置處的壓強P氣相；其中，z為正整數(shù)，且z=l、…、Q ; 步驟四、二次煙特征值計算，其計算過程如下: 步驟401、建模:采用所述數(shù)據(jù)處理器，建立固體推進劑的二次煙特征值計算模型，且所建立的二次煙特征值計算模型為二次煙相對濕度計算模型=RHamb =IOOX [K-ftotalX0.16589) ] (I)，式⑴中 /_/ = fHlo +/hc +Ihf，其中 ftQtal 為所評價固體推進劑燃燒后所產(chǎn)生燃燒產(chǎn)物中凝聚物在發(fā)動機噴管軸線上的摩爾百分含量，^匆所評價固體推進劑燃燒后所產(chǎn)生燃燒產(chǎn)物中的水在發(fā)動機噴管軸線上的摩爾百分含量，fHCi為所評價固體推進劑燃燒后所產(chǎn)生燃燒產(chǎn)物中的HCl在發(fā)動機噴管軸線上的摩爾百分含量，且fHF為所評價固體推進劑燃燒后所產(chǎn)生燃燒產(chǎn)物中的HF在發(fā)動機噴管軸線上的摩爾百分含量'K=S，式中Ptl為通過所述參數(shù)輸入單元預(yù)先輸入的0°C時水的飽和蒸氣壓， PoP=PHa+PHF，其中Pm為氣相產(chǎn)物HCl的氣體分壓且PHa=fHaXP，，PHF為氣相產(chǎn)物HF的氣體分壓且PHF=fHFXPn#，為步驟302中計算得出的發(fā)動機噴管軸線上的壓強； 步驟402、二次煙特征值計算:所述數(shù)據(jù)處理器調(diào)用步驟401中所建立的二次煙相對濕度計算模型，并結(jié)合步驟302中計算得出的所評價固體推進劑燃燒后所產(chǎn)生的各氣相產(chǎn)物在發(fā)動機噴管軸線上各位置處的摩爾百分含量fz，計算得出所評價固體推進劑燃燒后的發(fā)動機羽流場中噴管軸線上各位置處的二次煙相對濕度RH-。
2.按照權(quán)利要 求1所述的一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，其特征在于:步驟402計算得出所評價固體推進劑燃燒后的發(fā)動機羽流場中噴管軸線上各位置處的二次煙相對濕度RHamb后，所述數(shù)據(jù)處理器還需調(diào)用二次煙特征值評價模塊對發(fā)動機羽流場中噴管軸線上各位置處的二次煙等級進行評價，且RHamb值越大，二次煙等級越高。
3.按照權(quán)利要求1或2所述的一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，其特征在于:步驟402計算得出所評價固體推進劑燃燒后的發(fā)動機羽流場中噴管軸線上各位置處的二次煙相對濕度RHamb后，所述數(shù)據(jù)處理器還需調(diào)用曲線繪制模塊繪制出所評價固體推進劑燃燒后所產(chǎn)生二次煙的相對濕度RHamb的軸向分布曲線，并通過所述顯示單元對所繪制的相對濕度RHamb的軸向分布曲線進行同步顯示；其中，相對濕度RHamb的軸向分布曲線為所評價固體推進劑燃燒后所產(chǎn)生二次煙的相對濕度RHamb隨X軸的取值變化而變化的曲線；其中X軸為發(fā)動機噴管軸線所在的坐標軸,且發(fā)動機噴管出口處x=0。
4.按照權(quán)利要求1或2所述的一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，其特征在于:步驟302中計算得出的所評價固體推進劑燃燒后所產(chǎn)生的各氣相產(chǎn)物在發(fā)動機噴管軸線上各位置處的摩爾百分含量匕后，所述數(shù)據(jù)處理器還需曲線繪制模塊繪制出所評價固體推進劑燃燒后所產(chǎn)生各氣相產(chǎn)物在發(fā)動機噴管軸線上的摩爾百分含量分布曲線Iz ;其中，摩爾百分含量分布曲線Iz為所評價固體推進劑燃燒后所產(chǎn)生的第z種氣相產(chǎn)物的摩爾百分含量fz隨X軸的取值變化而變化的曲線；其中X軸為發(fā)動機噴管軸線所在的坐標軸，且發(fā)動機噴管出口處X=O。
5.按照權(quán)利要求1或2所述的一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，其特征在于:步驟401中pQ=610.78Pa。
6.按照權(quán)利要求2所述的一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，其特征在于:所述數(shù)據(jù)處理器調(diào)用所述二次煙特征值評價模塊對發(fā)動機羽流場中噴管軸線上各位置處的二次煙等級進行評價時，當RHamb值在區(qū)間(0.90，I]時，二次煙等級為微煙；當RHamb值在區(qū)間(0.55，0.90]時，二次煙等級為少煙；且當RHamb值在區(qū)間
時，二次煙等級為有煙。
7.按照權(quán)利要求1或2所述的一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，其特征在于:步驟201中所建立的最小自由能數(shù)學(xué)模型為
8.按照權(quán)利要求7所述的一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，其特征在于:步驟202中進行平衡組成計算時，所述數(shù)據(jù)處理器調(diào)用參數(shù)計算模塊，根據(jù)步驟一中所輸入的配制所評價固體推進劑所用各組分的化學(xué)式和質(zhì)量配比Hii,并結(jié)合各組分的相對分子質(zhì)量，對\.和1^.進行計算；之后，所述數(shù)據(jù)處理器結(jié)合預(yù)先輸入的^和」， 且根據(jù)公式
9.按照權(quán)利要求1或2所述的一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，其特征在于:步驟302中進行發(fā)動機羽流場計算時，所采用的發(fā)動機羽流場為二維軸對稱簡化 模型。
全文摘要
本發(fā)明公開了一種基于羽流場的固體推進劑二次煙霧特性虛擬評價方法，包括以下步驟一、初始參數(shù)設(shè)定與存儲；二、能量特性參數(shù)計算；三、發(fā)動機羽流場計算；四、二次煙特征值計算，數(shù)據(jù)處理器調(diào)用預(yù)先建立的二次煙相對濕度計算模型，并結(jié)合步驟三中的羽流場計算結(jié)果，計算得出所評價固體推進劑燃燒后的發(fā)動機羽流場中噴管軸線上各位置處的二次煙相對濕度RHamb。本發(fā)明方法步驟簡單、設(shè)計合理且實現(xiàn)方便、使用效果好，能有效解決現(xiàn)有固體推進劑煙霧信號性能評價時存在的人力物力投入大、周期長、試驗量大等缺陷。
文檔編號G06F19/00GK103186714SQ20131014155
公開日2013年7月3日 申請日期2013年4月22日 優(yōu)先權(quán)日2013年4月22日
發(fā)明者李猛, 趙鳳起, 張曉宏, 徐司雨, 羅陽, 郝海霞 申請人:西安近代化學(xué)研究所