采用三階WENO有限差分格式,耦合求解爆轟產(chǎn)物質(zhì)量分?jǐn)?shù)輸運(yùn)方程與可壓縮歐拉方程,并通過建立內(nèi)能與熱容之間的關(guān)系求得瞬態(tài)溫度,自主開發(fā)密閉空間內(nèi)炸藥爆炸溫度場三維數(shù)值計(jì)算程序�；谒_發(fā)的程序開展密閉空間內(nèi)炸藥爆炸溫度場數(shù)值計(jì)算,初步探討內(nèi)爆炸溫度場演化過程及其分布規(guī)律。研究表明:密閉空間內(nèi)爆炸溫度載荷在爆炸初期呈現(xiàn)出多峰值特征,在爆炸后期,逐漸趨近于準(zhǔn)靜態(tài)平穩(wěn)值,且內(nèi)爆炸溫度場分布具有空間上的不均勻性;數(shù)值計(jì)算準(zhǔn)靜態(tài)溫度與理論計(jì)算吻合較好,相對(duì)誤差控制在1%以內(nèi),初步驗(yàn)證了所開發(fā)程序的可靠性和正確性。本文的研究可為進(jìn)一步探討考慮后燃燒效應(yīng)爆炸溫度場數(shù)值計(jì)算及毀傷評(píng)估提供一定的參考和指導(dǎo)。

【文章頁數(shù)】：8 頁

【部分圖文】：

圖1Newton迭代法計(jì)算溫度流程圖Fig.1ProcedureoftemperaturecalculationbyNewtoninteraction

、C、D為爆轟產(chǎn)物和空氣按照質(zhì)量加權(quán)的混合參數(shù)，具體求解過程參考文獻(xiàn)[15]。式(5)是關(guān)于瞬態(tài)溫度T的一元非線性方程，在程序計(jì)算每一個(gè)時(shí)間步內(nèi)，進(jìn)一步令0022330001123fTeeATTDDBTTCTTTT(6)22DfTABTCTT(7)采用Newton迭代法[6]進(jìn)行....

圖3爆炸初始條件(a)壓力場(單位：GPa)(b)溫度場(單位：K)(pressurefield)(unit：GPa)(temperaturefield)(unit：K)

夢(mèng)?招?邊界(不考慮艙室結(jié)構(gòu)的變形)[16]，正交均勻網(wǎng)格尺寸取為40mm，見圖2(b)。(a)密閉艙室尺寸(單位：mm)(b)網(wǎng)格劃分(sizeofclosedcabin)(unit：mm)(meshdistribution)圖2密閉艙室尺寸及網(wǎng)格劃分Fig.2Sizeofcl....

圖5不同藥量內(nèi)爆炸測(cè)點(diǎn)溫度-時(shí)間歷程曲線Fig.5Temperature-timehistorycurveofdifferentgaugingpointsfordifferentmassexplosives(c)7218g(d)12195g

刂圃?1%以內(nèi)，滿足實(shí)際的工程計(jì)算要求，也初步驗(yàn)證了本文所開發(fā)的內(nèi)爆炸溫度場數(shù)值計(jì)算程序的可靠性和正確性。表2溫度體積積分平均與理論解相對(duì)誤差Tab.2Relativeerrorbetweentemperaturevolumeintegralaverageandtheoretic....

圖6不同藥量內(nèi)爆炸溫度體積積分平均時(shí)間歷程曲線與理論解對(duì)比Fig.6Comparisonoftemperaturevolumeintegralaveragetimehistorycurvewiththeoreticalsolutionsfordifferentmassexplosives(c)7218g(d)12195g

tweentemperaturevolumeintegralaverageandtheoreticalsolution藥量(mass)/g體積積分平均(volumeintegralaverage)/K理論解(theoreticalsolution)/K相對(duì)誤差(relativee....

本文編號(hào)：3982228