基于密閉空間內(nèi)爆的核爆沖擊波載荷模擬試驗(yàn)和數(shù)值計(jì)算方法

李 懌，李 典，侯海量，李永清

(海軍工程大學(xué) 艦船與海洋學(xué)院，武漢 430033)

當(dāng)前，開展核爆炸試驗(yàn)已極為困難，但提高艦船在核爆炸沖擊波載荷下的防護(hù)能力勢(shì)在必行，找到一種可行的方式來模擬核爆炸沖擊波載荷對(duì)提高艦船在核爆炸沖擊波下的防護(hù)能力有著重要意義。

20世紀(jì)七十年代初，國外已開始核爆沖擊波載荷的模擬試驗(yàn)研究。英國Foulness多噸位空氣沖擊波模擬裝置總度達(dá)173 m，可產(chǎn)生超壓低，正壓作用時(shí)間長(zhǎng)的沖擊波，用于模擬遠(yuǎn)距核爆沖擊波；法國于1981年在Gramat研究中心建成了當(dāng)時(shí)世界上最大的爆炸沖擊波模擬裝置(LBS)，該裝置可利用高壓空氣來模擬遠(yuǎn)距離核爆沖擊波[1]；加拿大國防研究與發(fā)展中心(DRDC)的爆炸波模擬裝置總長(zhǎng)45 m，能夠產(chǎn)生峰值超壓為0.3 MPa，正壓作用時(shí)間為100 ms的沖擊波與核爆沖擊波相似[2]；美國的大型爆炸波/熱輻射模擬裝置長(zhǎng)約273 m，能夠生成高溫平面波，可進(jìn)行核爆熱效應(yīng)和沖擊波模擬[3]。但上述裝置均為核爆典型爆炸沖擊波的大型模擬裝置，均存在建造困難、試驗(yàn)成本高等問題。

國內(nèi)關(guān)于核爆沖擊波載荷模擬試驗(yàn)研究起步晚。早期，胡昌[4-5]通過使用快速加載試驗(yàn)機(jī)模擬核爆沖擊波動(dòng)壓，對(duì)短柱進(jìn)行了快速加載，但其快速加載曲線與實(shí)際核爆沖擊波動(dòng)壓曲線相差較大，試驗(yàn)有較大偏差。北京環(huán)境強(qiáng)度研究所[6]研制出了可用于小型試驗(yàn)的光敏炸藥和薄片炸藥，可產(chǎn)生平面載荷，使核爆載荷模擬方式從物理加載跨越到了化爆模擬，在一定程度上提高了模擬的精度，但正壓作用時(shí)間過短。在接下來的十年中，并未產(chǎn)生新的核爆模擬方式。

由于核爆與化爆沖擊波具有相似的載荷曲線，但前者在相同比例距離下正壓作用時(shí)間與爆距均大得多，可視為平面波。因此，如何得到長(zhǎng)脈寬平面爆炸沖擊波成了模擬核爆的關(guān)鍵技術(shù)。

劉平等[7]通過數(shù)值仿真模擬了一維坑道中多層裝藥、延時(shí)爆炸得到的長(zhǎng)脈寬爆炸沖擊波，但是由于要設(shè)置合理的起爆藥量和坑道長(zhǎng)度，以及足夠的短得起爆時(shí)間間隔，導(dǎo)致試驗(yàn)實(shí)施十分困難。張秀華等[8]利用可燃?xì)怏w爆炸產(chǎn)生沖擊波對(duì)結(jié)構(gòu)/構(gòu)件加載，實(shí)現(xiàn)了長(zhǎng)脈寬沖擊波模擬，但未產(chǎn)生與核爆載荷相似的平面載荷。對(duì)于核爆載荷數(shù)值仿真，國內(nèi)也研究較少。袁丙方等[9-10]將核爆載荷壓力曲線直接作用在數(shù)值模型上，并未考慮載荷與試件的流固耦合作用。

本文設(shè)計(jì)了密閉爆坑內(nèi)化學(xué)爆炸模擬核爆炸沖擊波的試驗(yàn)，在典型測(cè)點(diǎn)處測(cè)得了反射壓力載荷特征曲線，并設(shè)計(jì)了一種新的考慮后燃燒效應(yīng)的數(shù)值計(jì)算方法，將數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果進(jìn)行對(duì)比，驗(yàn)證了數(shù)值計(jì)算方法的正確性和可行性。

1 理論分析

近年來，密閉空間內(nèi)爆炸的研究方興未艾[11-13]，其載荷作用過程包括爆炸沖擊波和準(zhǔn)靜態(tài)氣壓如圖1所示，其中準(zhǔn)靜態(tài)氣壓峰值相對(duì)較小，其衰減規(guī)律受密閉空間泄壓面積控制[14]，可以實(shí)現(xiàn)較長(zhǎng)的正壓作用時(shí)間。通過合理的密閉空間容積、裝藥、泄爆面積設(shè)計(jì)，控制初始沖擊波峰值和準(zhǔn)靜態(tài)氣壓衰減規(guī)律，可實(shí)現(xiàn)長(zhǎng)正壓作用時(shí)間核爆沖擊波的模擬。

1.1 試驗(yàn)準(zhǔn)靜態(tài)壓力分析

為明確試驗(yàn)藥量的大小，根據(jù)Calson公式[15]、勞氏軍規(guī)[16]、Mior[17]公式、JWL狀態(tài)方程、理想氣體狀態(tài)方程[18]以及文獻(xiàn)[19]中理論公式分別計(jì)算出達(dá)到150 kPa所需藥量如表1所示，考慮到黑索金與TNT之間能量的差異，需將爆炸黑索金的藥量換算成相同能量的TNT藥量，從表中可以看出不同方法計(jì)算結(jié)果差異較大，其中Calson公式、Mior公式、JWL狀態(tài)方程、理想氣體狀態(tài)方程均是基于理想狀態(tài)在密閉條件下推導(dǎo)而出，計(jì)算結(jié)果與試驗(yàn)可能有較大差異，而勞氏軍規(guī)與文獻(xiàn)[19]中公式為基于大量毀傷效應(yīng)實(shí)驗(yàn)數(shù)據(jù)(包括m/V<1 kg/m3)分析后給出的，已考慮實(shí)際各種因素的影響，出于試驗(yàn)要求，取兩者中的較大值，并考慮到細(xì)砂層的空隙對(duì)體積的影響，適量增大藥量，最終試驗(yàn)藥量定為750 g。

表1 不同公式150 kPa準(zhǔn)靜態(tài)壓力對(duì)應(yīng)藥量Tab.1 The corresponding dose of 150 kPa quasi-static pressure under different formulas

1.2 正壓作用時(shí)間分析

為確保其正壓時(shí)間達(dá)到200 ms以上，需對(duì)其泄壓面積進(jìn)行研究。根據(jù)文獻(xiàn)[14]中泄壓時(shí)間與泄壓面積的關(guān)系可知，當(dāng)正壓時(shí)間須達(dá)到200 ms以上時(shí)，泄壓面積應(yīng)小于0.061 16 m2，密閉爆坑上頂蓋的周長(zhǎng)為17.2 m，為滿足泄爆面積要求，頂蓋和側(cè)壁的縫隙需小于3.55 mm。

2 試驗(yàn)設(shè)計(jì)

2.1 爆坑與裝藥設(shè)計(jì)

采用內(nèi)爆的準(zhǔn)靜態(tài)壓載荷模擬長(zhǎng)時(shí)間作用的平面核爆沖擊波，爆坑長(zhǎng)寬均為4.3 m，深約3 m，底部鋪設(shè)細(xì)砂，試驗(yàn)試件埋設(shè)于爆坑底部，表面與細(xì)砂層平齊如圖2所示。細(xì)砂層與爆坑頂蓋間距1.4 m，利用軟木密封條對(duì)頂蓋與爆坑間空隙進(jìn)行密封，使頂蓋與側(cè)壁的縫隙不大于3 mm，以滿足泄爆面積要求，控制內(nèi)爆準(zhǔn)靜態(tài)氣壓衰減速率。

為得到平面壓力波載荷，降低初始爆炸沖擊波峰值超壓，保護(hù)爆坑上頂蓋，應(yīng)盡量采用平面裝藥。理想情況下，每平方米僅裝藥40.5 g，難以充分起爆，故采用平面“U”形布置的黑索金導(dǎo)爆索模擬平面裝藥加載如圖3所示，線裝藥量為25 g/m，爆速大于8 000 m/s，總長(zhǎng)30 m，總藥量為750 g，為防止爆坑頂蓋變形，導(dǎo)爆索布置平面與頂蓋的距離為0.95 m，與下方試件爆距為0.55 m。采用雷管在導(dǎo)爆索一端完成起爆。

1.空氣壓力傳感器；2.頂蓋；3.導(dǎo)爆索；4.試件；5.砂介質(zhì)；6.坑壁； 7.位移、應(yīng)變傳感器；8.鋼結(jié)構(gòu)支撐架；9.螺栓；10.側(cè)壁圖2 試驗(yàn)裝置Fig.2 Schematic diagram of the experimental device

(a)

2.2 測(cè)點(diǎn)布置

為考察受到載荷是否為平面載荷，在鋪設(shè)好的石英砂中對(duì)稱埋設(shè)4個(gè)壓力傳感器，且其承壓面與石英砂表面平齊，與爆坑頂面的距離1.4 m，距離爆坑邊壁≥5 cm，如圖4所示。

(a)

2.3 試驗(yàn)結(jié)果與分析

由于爆轟時(shí)間短(約3.75 ms)，爆轟產(chǎn)物的膨脹可視為是在爆轟完成后發(fā)生。爆轟產(chǎn)物初始膨脹是以導(dǎo)爆索為軸線的柱狀膨脹，極限膨脹半徑約為裝藥的30倍(約90 mm)，小于裝藥的布置間距，但以柱面形式傳播的沖擊波在相遇后會(huì)疊加形成平面波，向底部傳播如圖5所示。各壓力測(cè)點(diǎn)的載荷曲線如圖6所示，由圖可知，在該裝藥布置方式下，測(cè)點(diǎn)1,2,3,4的壓力曲線變化趨勢(shì)大致相同。在爆炸初期，壓力立刻上升至峰值超壓，隨后壓力降低，由于沖擊波在封閉腔體內(nèi)發(fā)生不斷的反射，壓力值仍較大。經(jīng)過不斷的反射，封閉空間內(nèi)形成較為均勻的準(zhǔn)靜態(tài)氣壓，隨著時(shí)間推移，準(zhǔn)靜態(tài)氣壓逐漸降低，最終歸為零，由圖可知準(zhǔn)靜態(tài)氣壓的持續(xù)時(shí)間約為230 ms。

圖5 爆轟產(chǎn)物與沖擊波傳播過程Fig.5 Schematic diagram of detonation products and shock wave changes

圖6 各測(cè)點(diǎn)壓力曲線Fig.6 Pressure curves of each test point

表2列出了各測(cè)點(diǎn)的反射壓力峰值，準(zhǔn)靜態(tài)氣壓值與正壓作用時(shí)間。由表2可知，測(cè)點(diǎn)1,2,3,4的反射超壓峰值差異較大，對(duì)四個(gè)測(cè)點(diǎn)值取平均值為728.13 kPa，四個(gè)測(cè)點(diǎn)可視為在均值上下震蕩，測(cè)點(diǎn)1、2、3、4分別與均值的相對(duì)偏差分別為8.26%、27.8%、37.75%、1.66%。由于在測(cè)點(diǎn)2，3傳感器上方存在導(dǎo)爆索爆炸，使其峰值偏大。四個(gè)測(cè)點(diǎn)正壓作用時(shí)間大致相同，其中測(cè)點(diǎn)1、2、3、4的正壓作用時(shí)間均為230 ms左右，測(cè)點(diǎn)2正壓作用時(shí)間最短為221 ms，測(cè)點(diǎn)4的正壓作用時(shí)間最長(zhǎng)為240 ms，相對(duì)偏差為9.1%。

漢代服飾是中國傳統(tǒng)服飾長(zhǎng)河中一顆璀璨的明珠，漢代服飾種類的多樣及其發(fā)展變化對(duì)中國傳統(tǒng)服飾的發(fā)展有著重要影響，漢代的開放和經(jīng)濟(jì)的發(fā)達(dá)以及漢代博大精深的文化使得漢代服飾也具有豐富的文化內(nèi)涵，不僅體現(xiàn)了漢代的時(shí)代特色，等級(jí)制度的森嚴(yán)而且反映出人們審美意識(shí)。它豐富的文化內(nèi)涵使得它在中國傳統(tǒng)服飾歷史中具有承上啟下的作用，我們應(yīng)該通過研究漢代服飾了解漢代的文化，從小的縮影中窺探整個(gè)中國的傳統(tǒng)服飾文化，同時(shí)應(yīng)該保護(hù)漢代服飾等一系列的傳統(tǒng)文化代表。

表2 各測(cè)點(diǎn)特征參數(shù)試驗(yàn)結(jié)果Tab.2 Each test point feature parameter test results

對(duì)測(cè)點(diǎn)1,2,3,4準(zhǔn)靜態(tài)壓力進(jìn)行分析，關(guān)于準(zhǔn)靜態(tài)壓力的取法，國內(nèi)外學(xué)者有著不同的看法，考慮到在準(zhǔn)靜態(tài)壓力階段，壓力值在不斷震蕩，本文采用在該階段壓力的平均值作為該階段的準(zhǔn)靜態(tài)壓力值[20]，如圖7所示。

(a) 測(cè)點(diǎn)2壓力曲線

根據(jù)以上方法測(cè)得四個(gè)測(cè)點(diǎn)準(zhǔn)靜態(tài)壓力值基本相同，準(zhǔn)靜態(tài)壓力值點(diǎn)的最大值為測(cè)點(diǎn)4處為157.58 kPa，最小值為測(cè)點(diǎn)處2為154.90 kPa，相對(duì)偏差僅為1.7%。取四個(gè)測(cè)點(diǎn)準(zhǔn)靜態(tài)壓力的平均值作為該密閉空間的準(zhǔn)靜態(tài)壓力，其大小為156.47 kPa。綜上所述，各測(cè)點(diǎn)的壓力載荷變化趨勢(shì)大致相同，且載荷作用時(shí)間與準(zhǔn)靜態(tài)氣壓值也大致相同，僅反射超壓載荷的峰值存在差別，故各測(cè)點(diǎn)的壓力載荷曲線具有較好的一致性。

對(duì)于爆炸中受到長(zhǎng)時(shí)間平面載荷的艦船結(jié)構(gòu)來說，超壓峰值對(duì)于整個(gè)爆炸過程的影響作用較小，對(duì)艦船結(jié)構(gòu)毀傷作用起其決定性影響的是爆炸中整個(gè)超壓的過程。為體現(xiàn)整個(gè)超壓過程對(duì)艦船結(jié)構(gòu)毀傷作用的大小，通過沖量的大小來反映整個(gè)超壓過程的毀傷作用的大小。表3是各個(gè)測(cè)點(diǎn)測(cè)得的壓力通過計(jì)算所得的沖量值，其中各個(gè)測(cè)點(diǎn)壓力值的相對(duì)零點(diǎn)是通過對(duì)在爆炸前各測(cè)點(diǎn)所測(cè)的震蕩波值求平均值所得，計(jì)算沖量時(shí)壓力作用時(shí)間的起點(diǎn)和終點(diǎn)分別為各個(gè)測(cè)點(diǎn)的壓力發(fā)生突變的時(shí)間點(diǎn)。通過表可知，各測(cè)點(diǎn)在爆炸過程中所受的沖量大致相同，其中測(cè)點(diǎn)3所受沖量最大，為21.618 kpa·s，測(cè)點(diǎn)1所受沖量最小，為16.914 kpa·s，其四個(gè)測(cè)點(diǎn)的平均值為19.966 kpa·s，測(cè)點(diǎn)1、2、3、4與平均值的相對(duì)偏差分別為15.3%、3.2%、10.2%、8.3%，根據(jù)圖6分析，造成各測(cè)點(diǎn)沖量偏差的主要階段為準(zhǔn)靜態(tài)壓力下降階段，根據(jù)試驗(yàn)環(huán)境分析，由于爆坑頂蓋密封性在四條邊長(zhǎng)處存在差異，即爆坑四邊的泄爆面積不同，且爆坑內(nèi)容積較大，靠近不同邊的壓力傳感器處的壓力下降受各邊泄爆面積影響，造成泄爆的速率不同，故在泄爆階段，各壓力傳感器所測(cè)得的壓力曲線有所差異，最后爆坑內(nèi)壓力達(dá)到一致。

表3 各測(cè)點(diǎn)沖量Tab.3 Impulse of each test point

3 數(shù)值計(jì)算

在仿真過程中發(fā)現(xiàn)用高能炸藥燃燒模型和JWL狀態(tài)方程來表征RDX炸藥的材料模型，當(dāng)計(jì)算時(shí)長(zhǎng)達(dá)到230 ms時(shí)，計(jì)算時(shí)間過長(zhǎng)無法計(jì)算，且準(zhǔn)靜態(tài)壓力值均值為98.68 kPa如圖8所示，與試驗(yàn)偏差為36.9%。根據(jù)文獻(xiàn)[21]研究結(jié)果，藥量體積比較小情況下，后燃燒效應(yīng)對(duì)密閉空間爆炸準(zhǔn)靜態(tài)氣壓影響較大。為驗(yàn)證試驗(yàn)中后燃燒放熱對(duì)準(zhǔn)靜態(tài)壓力的影響，利用ANSYS-LSDYNA軟件有限元軟件對(duì)該封閉空間導(dǎo)爆索爆炸進(jìn)行數(shù)值模擬。

圖8 未考慮后燃燒效應(yīng)數(shù)值計(jì)算曲線Fig.8 The numerical calculation curve of the post-combustion stage is not considered

3.1 模型建立

利用ANSYS-LSDYNA有限元軟件對(duì)試驗(yàn)工況進(jìn)行建模。由于試驗(yàn)爆坑過大且爆坑關(guān)于x、y軸對(duì)稱，故建立1/4對(duì)稱模型，并對(duì)導(dǎo)爆索所在區(qū)域進(jìn)行局部網(wǎng)格加密處理，網(wǎng)格大小為2 mm，將計(jì)算時(shí)長(zhǎng)設(shè)為230 ms，時(shí)間步長(zhǎng)間隔為0.01 ms；模型中的空氣域采用歐拉網(wǎng)格建模，邊界采用固定邊界，為模擬實(shí)際試驗(yàn)情況，空氣域上表面四周邊緣采用無反射的OUTFLOW邊界，以模擬頂蓋和側(cè)壁的泄爆縫隙。

忽略導(dǎo)爆索的爆轟過程與爆轟產(chǎn)物的膨脹過程，可將裝藥等效為片狀高壓氣團(tuán)。根據(jù)爆炸過程中質(zhì)量守恒和能量守恒，由排布方式可將導(dǎo)爆索定義為3 000 mm×2 600 mm×2 mm的高壓氣團(tuán)如圖9所示。

圖9 1/4數(shù)值計(jì)算模型Fig.9 Numerical calculation model

采用MAT_NULL模型和EOS_LINEAR_POLYNOMIAL狀態(tài)方程描述高壓氣團(tuán)和空氣域模型，線性多項(xiàng)式狀態(tài)方程如式(1)所示

P=C0+C1μ+C2μ2+(C4+C5μ+C6μ2)E0

(1)

式中：C0～C6為該線性多項(xiàng)式的系數(shù)，對(duì)于理想氣體而言，C0=C1=C2=C3=C6=0,C4=C5=γ-1,其中γ取1.4；E為單位體積下的初始內(nèi)能，對(duì)于空氣，E0=2.525×105J·m-3，對(duì)于高壓氣團(tuán)，由能量守恒定可知，其總內(nèi)能與750 g導(dǎo)爆索在爆炸階段和后燃燒效應(yīng)能量相同。

為計(jì)算后燃燒效應(yīng)放熱，假設(shè)密閉空間體積為V，炸藥的質(zhì)量為WRDX，為計(jì)算爆炸后產(chǎn)物反應(yīng)是否完全，引入一參數(shù)δ=WRDX/V,當(dāng)δ≤x時(shí)，RDX炸藥在密閉空間內(nèi)能夠完全反應(yīng)，生成反應(yīng)產(chǎn)物CO2、H2O、N2。通過計(jì)算可得：x=0.636 1 kg/m3，當(dāng)δ≤0.636 1時(shí)，RDX炸藥在密閉空間內(nèi)能夠完全反應(yīng)，在本試驗(yàn)中δ=0.028 97，故本試驗(yàn)需考慮后燃燒效應(yīng)產(chǎn)生熱量，對(duì)于RDX炸藥，通過查詢文獻(xiàn)[22]可知，其爆熱值為5 820 kJ·kg-1,其由文獻(xiàn)[23]可知，依據(jù)最大放熱規(guī)則，對(duì)于RDX炸藥爆轟產(chǎn)生C的物質(zhì)的量nC與RDX物質(zhì)的量1.5nRDX的關(guān)系為nC=1.5nRDX,其后燃燒反應(yīng)為C→CO2，反應(yīng)放出熱量為393.5×103J，則有

EAV=393.5nC×103

(2)

ETV=ERDX+EAV=5.820×106WRDX+

393.5nC×103

(3)

(4)

式中:EAV為后燃燒能量；ERDX為RDX炸藥爆熱；ETV為密閉爆坑內(nèi)爆炸釋放的總能量。故高壓氣團(tuán)單位體積的初始內(nèi)能為E0=4.076×108J·m-3；空氣的密度取1.280 kg·m-3；高壓氣團(tuán)的質(zhì)量為750 g，經(jīng)計(jì)算可得其密度為48.07 kg·m-3。

3.2 結(jié)果分析

對(duì)試驗(yàn)工況進(jìn)行數(shù)值計(jì)算，取典型測(cè)點(diǎn)的超壓隨時(shí)間變化曲線如圖10所示，其中數(shù)值計(jì)算測(cè)點(diǎn)的選取是根據(jù)實(shí)際測(cè)點(diǎn)的相同位置選取且由于是對(duì)稱模型故只取兩個(gè)測(cè)點(diǎn)即可。從圖10和表4可知，兩個(gè)測(cè)點(diǎn)沖擊波超壓曲線變化趨勢(shì)一致，峰值壓力均為659.531 kPa。為求得測(cè)點(diǎn)測(cè)得的準(zhǔn)靜態(tài)壓力情況將測(cè)點(diǎn)30 ms測(cè)得曲線取平均值，可得出測(cè)點(diǎn)1，2處的準(zhǔn)靜態(tài)壓力峰值分別為148.43 kPa和146.66 kPa，兩者的相對(duì)偏差僅為3.79%。

(a) 數(shù)值計(jì)算測(cè)點(diǎn)1、2壓力曲線

表4 數(shù)值計(jì)算測(cè)點(diǎn)1、2數(shù)據(jù)Tab.4 Numerical calculation data of test points 1 and 2

3.3 偏差分析

圖11是試驗(yàn)各測(cè)點(diǎn)壓力曲線和數(shù)值計(jì)算值壓力曲線。由圖可知，爆炸沖擊波階段，由于試驗(yàn)過程中各種因素影響，導(dǎo)致測(cè)點(diǎn)1、2、3、4沖擊波峰值不同，其平均值為728.13 kPa，數(shù)值計(jì)算所測(cè)得沖擊波峰值為659.53 kPa,相對(duì)偏差為9.42%；準(zhǔn)靜態(tài)壓力階段，數(shù)值計(jì)算前30 ms曲線與試驗(yàn)曲線擬合較好如圖11所示，試驗(yàn)中各測(cè)點(diǎn)所測(cè)得準(zhǔn)靜態(tài)壓力平均值為156.46 kPa，數(shù)值計(jì)算所得平均值為147.55 kPa，相對(duì)偏差為5.70%；在泄爆階段，數(shù)值計(jì)算壓力曲線與試驗(yàn)各測(cè)點(diǎn)壓力曲線相差較大，是由于數(shù)值計(jì)算與實(shí)際工況的泄爆面積不符，數(shù)值計(jì)算的泄爆面積較大，導(dǎo)致準(zhǔn)靜態(tài)壓力下降較快。結(jié)果說明本文采用用高壓氣團(tuán)代替裝藥的方法進(jìn)行數(shù)值計(jì)算能夠很好地模擬出裝藥的壓力曲線，具有一定的可靠性；對(duì)于當(dāng)質(zhì)量體積比δ≤0.636 1的RDX炸藥，考慮其后燃燒效應(yīng)放出熱量的影響，計(jì)算偏差減小了31.7%。

(a) 試驗(yàn)與數(shù)值計(jì)算壓力對(duì)比曲線

4 結(jié) 論

本文設(shè)計(jì)了密閉爆坑內(nèi)化學(xué)爆炸模擬核爆炸沖擊波的試驗(yàn)方法，在典型測(cè)點(diǎn)處測(cè)得了反射超壓的壓力載荷特征，并設(shè)計(jì)了一種考慮后燃燒效應(yīng)的新式數(shù)值計(jì)算方法，將數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果進(jìn)行對(duì)比，驗(yàn)證了數(shù)值計(jì)算方法的有效性和正確性，得到有關(guān)結(jié)論如下：

(1) 密閉爆坑內(nèi)小當(dāng)量U型導(dǎo)爆索爆炸下各測(cè)點(diǎn)壓力載荷有較好的一致性，能有效形成平面長(zhǎng)脈寬沖擊載荷，與核爆載荷特性相近。

(2) 在密閉空間中，對(duì)于質(zhì)量體積比δ≤0.636 1的RDX炸藥，考慮其后燃燒效應(yīng)能量考慮其后燃燒效應(yīng)放出熱量較不考慮后燃燒偏差減小了31.7%,能夠有效地減小數(shù)值計(jì)算的偏差。

(3) 在數(shù)值計(jì)算中，通過用高壓氣團(tuán)代替?zhèn)鹘y(tǒng)的炸藥進(jìn)行數(shù)值計(jì)算，試驗(yàn)和數(shù)值計(jì)算值對(duì)比發(fā)現(xiàn)，試驗(yàn)所測(cè)得超壓峰值與數(shù)值計(jì)算峰值相對(duì)偏差在10%以內(nèi)，試驗(yàn)測(cè)得準(zhǔn)靜態(tài)壓力值與數(shù)值計(jì)算準(zhǔn)靜態(tài)壓力值在5.7%，偏差在合理范圍內(nèi)，說明該數(shù)值計(jì)算方法具有一定的可靠性和正確性。