• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    爆炸二極管的數(shù)值模擬設(shè)計(jì)

    2017-05-07 06:06:48馬宏昊沈兆武任麗杰
    含能材料 2017年4期

    崔 宇, 馬宏昊, 沈兆武, 任麗杰, 楊 明, 余 勇

    (1. 中國科學(xué)院材料力學(xué)行為和設(shè)計(jì)重點(diǎn)實(shí)驗(yàn)室, 安徽 合肥 230026; 2. 中國科學(xué)技術(shù)大學(xué), 安徽 合肥 230026)

    工程爆破隨著經(jīng)濟(jì)發(fā)展而起步、發(fā)展、壯大并且在采礦、交通等領(lǐng)域取得了巨大的經(jīng)濟(jì)效益[1],工程中發(fā)生的事故,特別是炮孔、起爆網(wǎng)路的早爆、誤爆問題嚴(yán)重危害了工人生命安全并延誤工程進(jìn)度,這促使研究人員密切關(guān)注爆破安全技術(shù)[2-3]。張英梅[4]等對煤礦井下的雜散電流進(jìn)行數(shù)學(xué)建模,分析了雜散電流引發(fā)雷管早爆的機(jī)理。龍昌軍[5]等在拆除爆破中采用隔熱結(jié)構(gòu)來防止高溫導(dǎo)致的炮孔早爆。曹春杰[6]等為防止高溫火區(qū)環(huán)境中炮孔早爆,將雷管和炸藥進(jìn)行隔熱處理。束學(xué)來[7]等提出耐熱炸藥以應(yīng)對高溫環(huán)境下的早爆問題。

    上述研究多是采用隔熱保護(hù)手段來防止早爆問題,并沒有提出一種針對起爆網(wǎng)路安全性的手段。如果能設(shè)計(jì)出一種裝置: 當(dāng)激發(fā)信號從主網(wǎng)絡(luò)正向傳遞過來能順利起爆炮孔裝藥,而單個(gè)炮孔裝藥早爆,誤爆的激發(fā)信號不能傳遞給主網(wǎng)絡(luò),不會引起整個(gè)起爆網(wǎng)路的起爆,這將極大降低事故損失。王飛[3]等設(shè)計(jì)一種利用不同感度多級裝藥結(jié)構(gòu)的單向起爆元件,結(jié)構(gòu)較為復(fù)雜且利用鈍感炸藥作為隔爆材料保護(hù)起爆網(wǎng)路也會引發(fā)新的安全問題。為此,本研究設(shè)計(jì)了一種簡單易行的爆炸二極管,該結(jié)構(gòu)利用飛片沖擊起爆實(shí)現(xiàn)正向傳爆,利用炸藥在鉛質(zhì)線延期體[8]元件處爆炸時(shí),透射沖擊波反向隔爆。同時(shí)采用數(shù)值計(jì)算以及實(shí)驗(yàn)針對該爆炸二極管的激發(fā)藥藥量以及延期體長度等關(guān)鍵因素進(jìn)行了分析。

    2 爆炸二極管的仿真設(shè)計(jì)

    2.1 爆炸二極管的結(jié)構(gòu)

    設(shè)計(jì)的爆炸二極管結(jié)構(gòu)如圖1所示,爆轟信號從起爆網(wǎng)絡(luò)主網(wǎng)路導(dǎo)爆索輸入,經(jīng)由延期體、激發(fā)藥,剪切內(nèi)帽底部形成高速飛片,撞擊傳爆藥,最后傳遞給輸出端導(dǎo)爆索,最終起爆炮孔裝藥。當(dāng)炮孔裝藥意外起爆時(shí),爆轟信號從輸出端導(dǎo)爆索逆向輸入,經(jīng)由傳爆藥起爆激發(fā)藥,最后沖擊壓縮鉛質(zhì)延期體,雖然會點(diǎn)燃延期體中藥柱,但是能量已經(jīng)不足以起爆輸出端導(dǎo)爆索,必然可以在主網(wǎng)絡(luò)之外阻斷爆轟信號,達(dá)到隔爆的目的。其中延期體、內(nèi)帽、管殼外力擠壓固連,輸入端、輸出端導(dǎo)爆索也與管殼固連。

    圖1 爆炸二極管的結(jié)構(gòu)圖

    1—固連處, 2—導(dǎo)爆索, 3—內(nèi)帽, 4—延期體, 5—激發(fā)藥, 6—定位臺階, 7—傳爆藥, 8—管殼

    Fig.1 The structure of the explosive diode

    1—fixed points, 2—detonating fuse, 3—internal cylinder, 4—delay element, 5—firing charge, 6—step positioned, 7—booster charge, 8—tube

    2.2 正向傳爆中激發(fā)藥極限藥量仿真設(shè)計(jì)

    針對圖1設(shè)計(jì)的爆炸二極管,現(xiàn)在采用點(diǎn)火增長模型[9]結(jié)合LS-DYNA軟件,對激發(fā)藥藥量進(jìn)行正向傳爆穩(wěn)定性分析。

    激發(fā)藥采用高能燃燒計(jì)算模型和JWL狀態(tài)方程[10],如式1所示。

    (1)

    式中,p為爆轟壓力,GPa,具體參數(shù)如表1所示,激發(fā)藥為太安(PETN)。

    表1 激發(fā)藥爆轟參數(shù)[10]

    Table 1 The detonation parameters of firing charge[10]

    ρ/g·cm-3D/km·s-1pCJ/GPaABR1R2ω0.885.176.23.4860.1137.02.00.24

    內(nèi)帽與殼體為鋼材料,密度為7.83 g·cm-3,在正向傳爆模擬中可以將延期體處理為封閉端與殼體固連。傳爆藥采用JWL狀態(tài)方程以及Lee&Traver三項(xiàng)式點(diǎn)火與增長反應(yīng)速率模型[11]。該模型可以反映受沖擊的炸藥內(nèi)部的爆轟發(fā)展過程,具體形式如式2所示。

    G21-λeλgpz

    (2)

    式中,λ為反應(yīng)度,即為炸藥的反應(yīng)百分?jǐn)?shù),I、a、b、c、d、e、g、x、y、z、G1、G2為可調(diào)參數(shù),采用1.64 g·cm-3的TNT作為傳爆藥,具體參數(shù)如表2所示。

    表2 TNT點(diǎn)火與增長反應(yīng)速率模型參數(shù)[12]

    Table 2 Ignition and Growth reactive flow model parameters for TNT[12]

    IG1G2abcdegxyz5036010000.6671.00.6671.00.11141.21

    采用1/4模型計(jì)算,建模單位g-cm-μs,簡化模型如圖2所示。由于涉及炸藥大變形,采用Lagrange方法出現(xiàn)網(wǎng)格畸變,計(jì)算失敗,所以采用ALE方法。激發(fā)藥(50 mg時(shí)109512個(gè)單元)、傳爆藥(298480個(gè)單元)和空氣域(435136個(gè)單元)采用ALE網(wǎng)格,內(nèi)帽(20176個(gè)單元)與殼體(113464個(gè)單元)采用Lagrange網(wǎng)格,兩種網(wǎng)格間采用流固耦合算法,內(nèi)帽與殼體間共節(jié)點(diǎn)作用。均為映射六面體網(wǎng)格,共976768個(gè)單元。

    圖2 正向傳爆計(jì)算結(jié)構(gòu)模型

    Fig.2 The calculation model for transmitting positive detonation

    調(diào)整激發(fā)藥藥量10,20,30,40,50 mg,控制傳爆藥長度均為5 mm,飛片厚0.2 mm,殼體壁厚0.4 mm,自上而下起爆。其中50 mg藥量模擬過程如圖3所示(1)激發(fā)藥的起爆; (2)高速飛片生成; (3)飛片撞擊傳爆藥; (4)傳爆藥中爆轟發(fā)展。圖4、圖5分別為50,10 mg激發(fā)藥時(shí)傳爆藥中軸線上爆轟發(fā)展過程。對比兩圖,可以看出隨著距撞擊面長度的增長,壓力均是先增大。3.3 μs之后,50 mg藥量下傳爆藥的峰值壓力維持18 GPa左右,對比1.64 g·cm-3的TNT的C-J壓力約21 GPa[12],判定傳爆成功; 而10 mg藥量下,峰值壓力在4.2 μs后反而開始降低,爆轟信號衰減,且軸線上峰值中最大值為5.2 GPa,遠(yuǎn)小于C-J壓力,判定傳爆失敗。表明激發(fā)藥藥量決定了正向傳爆的穩(wěn)定性。

    不同激發(fā)藥藥量結(jié)果見表3。由表3可見,隨著激發(fā)藥藥量的增大,飛片最大速度逐步遞增,飛片撞擊傳爆藥的初始壓力也在遞增,傳爆藥藥柱中的最大壓力也在逐步遞增。當(dāng)撞擊壓力達(dá)到6.3 GPa時(shí),可以看到傳爆藥藥柱完全爆轟,峰值壓力達(dá)到18.6 GPa,此時(shí)峰值對應(yīng)的位置為距離撞擊端面2.6 mm。而當(dāng)藥量小于10 mg時(shí),壓力峰值將會在4.7 mm后出現(xiàn)衰減,即爆轟信號衰減傳爆失敗??芍?當(dāng)激發(fā)藥藥量大于50 mg,傳爆藥完全爆轟,小于10 mg傳爆失敗。

    a. 0.39 μs b. 1.13 μs

    c. 2.26 μs d. 2.71 μs

    圖3 50 mg激發(fā)藥下的計(jì)算過程

    Fig.3 The calculation process under the 50 mg firing charge

    圖4 50 mg激發(fā)藥下的爆轟成長過程

    Fig.4 The detonation development process under the 50 mg firing charge

    2.3 反向隔爆中延期體的臨界尺寸仿真設(shè)計(jì)

    反向隔爆設(shè)計(jì)中,鉛質(zhì)延期體起主要隔爆作用。與正向傳爆的順序相反,當(dāng)炮孔意外引爆時(shí),會將傳爆藥起爆,傳爆藥將會反向沖擊壓縮內(nèi)帽起爆激發(fā)藥,最后一起沖擊壓縮鉛質(zhì)延期體。采用如圖6的1/4簡化模型,建模單位g-cm-μs,仍采用ALE方法。輸入端導(dǎo)爆索(25350個(gè)單元)與激發(fā)藥(8112個(gè)單元)均被簡化為帶有點(diǎn)火與增長反應(yīng)速率模型(表2所示)的TNT[12],5 mm長傳爆藥(25350個(gè)單元)則直接采用高能炸藥燃燒模型以及JWL狀態(tài)方程[12]。內(nèi)帽、延期體、殼體采用共節(jié)點(diǎn)方式固連在一起,采用Lagrange網(wǎng)格,炸藥與空氣則采用ALE網(wǎng)格,兩者之間采用流固耦合算法,延期體鉛(5 mm下為25350個(gè)單元)使用SCG模型[13],具體參數(shù)如文獻(xiàn)[13]所示,所有網(wǎng)格均為映射六面體網(wǎng)格,共181428個(gè)單元。

    圖5 10 mg激發(fā)藥下的爆轟成長過程

    Fig.5 The detonation development process under the 10 mg firing charge

    表3 不同激發(fā)藥藥量下的計(jì)算物理量

    Table 3 The calculated physical quantity with the various mass of firing charge

    massoffiringcharge/mgmaximumvelocityofflyingplanes/m·s-1collisionpressure/GPapressurepeak/GPapeakposition/mm109001.45.34.72011402.010.25.03012302.913.35.04015004.214.45.05015506.318.62.6

    圖6 反向隔爆計(jì)算結(jié)構(gòu)模型

    Fig.6 The calculation structure model for stopping back-ward detonation

    控制激發(fā)藥30 mg,調(diào)節(jié)延期體長度為1.0,2.0,3.0,4.0,5.0 mm,自下而上起爆。其中延期體1.0 mm時(shí)的計(jì)算過程如圖7所示: (1)傳爆藥爆轟,作為反向爆轟輸入信號; (2)內(nèi)帽中裝藥起爆; (3)延期體變形; (4)、(5)輸入端裝藥的爆轟發(fā)展過程。圖8、圖9分別為1、5 mm延期體時(shí),軸線上激發(fā)藥(虛線)和輸入端裝藥(實(shí)線)的爆轟發(fā)展過程。對比發(fā)現(xiàn),激發(fā)藥壓力均近似15 GPa,說明激發(fā)藥直接被起爆。而1 mm下輸入端裝藥壓力也達(dá)到15 GPa,說明該處裝藥也被起爆,隔爆失敗; 5 mm時(shí)輸入端裝藥卻由于初始壓力僅有2.4 GPa,而逐步衰減至熄爆說明隔爆成功。這是因?yàn)閮?nèi)帽底僅0.2 mm厚,遠(yuǎn)小于延期體長度。表明: 延期體長度對隔爆起決定性作用。

    c. 0.34 μs b. 1.08 μs c. 1.38 μs d. 1.79 μs e. 2.56 μs

    圖7 1.0 mm延期體下的計(jì)算過程

    Fig.7 The calculation process under the 1.0 mm delay element

    圖8 1 mm延期體下的爆轟成長過程

    Fig.8 The detonation development process under the 1.0 mm delay element

    圖9 5 mm延期體下的爆轟成長過程

    Fig.9 The detonation development process under the 5.0 mm delay element

    結(jié)果如表4所示: 隨著延期體長度增加,激發(fā)藥初始壓力基本不變,且均直接爆轟而輸入端裝藥初始壓力則逐步降低,輸入端裝藥中的壓力峰值也由15.8 GPa降到2.5 GPa,峰值出現(xiàn)的位置也由距撞擊端面5 mm至撞擊端面上??梢钥闯鲅悠隗w鉛柱長度低于4 mm不能隔爆,5 mm以上可以反向阻爆。

    表4 不同延期體長度下的計(jì)算物理量

    Table 4 The calculated physical quantity with the various length of the delay elements

    lengthofdelayelements/mminitialpressureincharge/GPamaximalvelocityofdelayelements/m·s-1collisionpressure/GPapressurepeak/GPapeakposition/mm1.014.88705.8615.85.02.014.67504.9314.35.03.014.56483.7111.65.04.014.75333.025.72.65.014.74272.482.50

    3 實(shí) 驗(yàn)

    3.1 正反對稱驗(yàn)證實(shí)驗(yàn)

    考慮到正向傳爆以及反向隔爆本身具有對稱性,而且同時(shí)進(jìn)行正反實(shí)驗(yàn)簡單高效、節(jié)省材料以及對比鮮明,所以本文設(shè)計(jì)一種對稱裝藥結(jié)構(gòu)驗(yàn)證模擬結(jié)果以及設(shè)計(jì)的可行性。如圖10所示整體結(jié)構(gòu)左右對稱,可以任選一端如A作為輸入端,則另外一端為輸出端。從A輸入爆轟信號,經(jīng)由延期體、激發(fā)藥、傳爆藥,信號傳遞到B管殼中,對應(yīng)正向傳爆試驗(yàn); 在B管殼中經(jīng)由傳爆藥、激發(fā)藥、延期體,對應(yīng)反向隔爆實(shí)驗(yàn)。如圖11a所示,第一段管殼為爆轟信號發(fā)生器,產(chǎn)生爆轟信號輸入第二段管殼,其中第二段管殼與第三段管殼為對稱裝藥結(jié)構(gòu),第三段管殼后的導(dǎo)爆索綁在驗(yàn)證板上,為薄鐵板。PETN作為激發(fā)藥不如RDX[14]故采用50 mg的低密度RDX,考慮工程上可靠因子的作用,傳爆藥采用PETN(1.0 g·cm-3)長度為10.0 mm,延期體長度約5.0 mm。引爆爆轟信號發(fā)生器,結(jié)果如圖11b所示,第二段管殼完全爆轟,第三段管殼碎裂至延期體位置,驗(yàn)證板基本完好,最后的輸出端導(dǎo)爆索也完好。重復(fù)實(shí)驗(yàn)10次,結(jié)果均如此。對比模擬結(jié)果50 mg激發(fā)藥正向穩(wěn)定傳爆,5 mm延期體反向可靠隔爆,實(shí)驗(yàn)驗(yàn)證了模擬結(jié)果,而且實(shí)驗(yàn)的重復(fù)性好,結(jié)果可信,即設(shè)計(jì)的爆炸二極管能達(dá)到正向穩(wěn)定傳爆反向可靠隔爆的效果。

    圖10 對稱裝藥結(jié)構(gòu)

    Fig.10 The symmetrical charge structure

    a. before experiment b. after experiment

    圖11 實(shí)驗(yàn)前后對照圖

    Fig.11 Comparison diagrams of before and after the experiment

    3.2 飛片速度測量

    由表3數(shù)據(jù)可以看出,飛片速度與激發(fā)藥藥量、撞擊壓力等呈正相關(guān)。而且由炸藥對飛片的一維拋射理論[15]有:

    (4)

    式中,Vm為飛片的極限速度,m·s-1;η=16m/(27M),D為炸藥爆速,m·s-1; 利用公式(4)對η求導(dǎo)得:

    (5)

    對(5)式恒大于零,也表明當(dāng)飛片質(zhì)量不變時(shí)隨著裝藥質(zhì)量的增大,飛片的速度越大。由此可以通過測量飛片速度來反映不同激發(fā)藥藥量下的起爆能量。采用斷通法[16],裝置如圖12所示,去除傳爆藥,激發(fā)藥爆轟剪切內(nèi)帽形成飛片,飛片打穿雙面鋁箔時(shí)產(chǎn)生一個(gè)電信號,由兩個(gè)電信號間的時(shí)間差與兩層雙面鋁箔間的距離即可得到飛片速度。結(jié)果如表5所示。隨著激發(fā)藥藥量增大,飛片速度先升高再降低,與數(shù)值計(jì)算結(jié)果以及理論有一定差距。分析如下: 一方面由于激發(fā)藥并不是直接爆轟而是燃燒轉(zhuǎn)爆轟,以及空氣阻力的影響,必然會使實(shí)際測量結(jié)果偏小; 另一方面由于飛片被打碎以及殼體破裂,側(cè)向稀疏波的影響,也會導(dǎo)致速度偏小,正如圖13所示,分別代表藥量由10 mg增加到50 mg,1#~5#飛片破損越嚴(yán)重,殼體側(cè)向膨脹和破碎越明顯。

    圖12 測速實(shí)驗(yàn)裝置圖

    Fig.12 The measuring equipment of detecting velocity

    表5 不同藥量下的飛片速度

    Table 5 The velocity of flying planes with the various mass of firing charge

    massoffiringcharge/mglengthofthedelayelement/mmaveragevelocityofflyingplanes/km·s-1105.510.79205.511.03305.521.2405.490.97505.500.98

    圖13 不同藥量下收集到的管殼和飛片

    Fig.13 The tubes and flying planes collected with the various mass of firing charge

    4 結(jié) 論

    (1) 通過數(shù)值模擬與正反對稱裝藥實(shí)驗(yàn)相結(jié)合的方法設(shè)計(jì)了一種爆炸二極管,實(shí)現(xiàn)了正向穩(wěn)定傳爆,反向可靠隔爆的目標(biāo)。

    (2) 模擬結(jié)果顯示,當(dāng)藥量大于50 mg,激發(fā)藥正向穩(wěn)定傳爆,小于10 mg則會傳爆失敗; 30 mg激發(fā)藥藥量下,延期體長度小于4 mm時(shí),反向隔爆失敗,大于5 mm時(shí)可以可靠隔爆。實(shí)驗(yàn)結(jié)果顯示50 mg激發(fā)藥,5 mm延期體可以達(dá)到100%整體成功率。

    (3) 飛片測量結(jié)果顯示,由于激發(fā)藥燃燒轉(zhuǎn)爆轟、空氣阻力殼體飛片碎裂等影響,飛片速度較模擬值偏低,但是仍可以看出10 mg激發(fā)藥速度只有790 m·s-1,必然會導(dǎo)致正向傳爆失敗。

    參考文獻(xiàn):

    [1] 馮叔瑜, 鄭哲敏. 讓工程爆破技術(shù)更好地服務(wù)社會、造福人類——我國工程爆破60年回顧與展望[J]. 中國工程科學(xué), 2014(11): 5-13.

    FENG Shu-yu, ZHENG Zhe-ming. Let the engineering blasting technology serving for society and benefit of mankind better:Sixty years of engineering blasting in China, review and prospect[J].EngineeringScience, 2014(11): 5-13.

    [2] 王付景.論城市控制爆破中的安全管理[J]. 低碳世界, 2015 (23): 264-265.

    WANG Fu-jing. The safety management in the urban controlled blasting[J].LowCarbonWorld, 2015 (23): 264-265.

    [3] 王飛, 馬宏昊, 沈兆武.一種應(yīng)用于高溫火區(qū)爆破中的不可逆起爆網(wǎng)路的耐溫與傳爆性能[J]. 含能材料, 2016, 24(1): 106-110

    WANG Fei, MA Hong-hao, SHEN Zhao-wu. Heat resistance and explosive transfer performance in the irreversible detonating network applied in high temperature mine blasting[J].ChineseJournalofEnergeticMaterials(HannengCailiao), 2016, 24(1): 106-110.

    [4] 張英梅, 袁曉光. 煤礦井下直流牽引網(wǎng)絡(luò)中雜散電流分布[J]. 太原理工大學(xué)學(xué)報(bào), 1998(3): 237-239.

    ZHANG Ying-mei, YUAN Xiao-guang. The stray current distribution in DC traction power supply system in the coal mine[J].Journaloftaiyuanuniversityoftechnology, 1998(3): 237-239.

    [5] 龍昌軍, 農(nóng)青鳳, 梁建銘, 等. 硅錳冶煉爐爐底凝結(jié)物拆除控制爆破[J]. 現(xiàn)代礦業(yè), 2015(10): 49-51.

    LONG Chang-jun, NONG Qing-feng, LIANG Jian-ming. Demolition blasting of the smelting furnaces with the concretion for silicon and manganese[J].ModernMining, 2015(10): 49-51.

    [6] 曹春杰. 露天煤礦采空區(qū)上方作業(yè)安全厚度的確定及爆破設(shè)計(jì)[J]. 煤炭工程, 2013, 45(9): 25-27.

    CAO Chun-jie. Determination of security coping thickness and blasting design over goafs in surface mine[J].CoalEngineering, 2013, 45(9): 25-27.

    [7] 束學(xué)來, 鄭炳旭, 郭子如, 等. 耐熱炸藥機(jī)理分析與優(yōu)化淺析[J]. 工程爆破, 2014, 20(5): 59-63.

    SU Xue-lai, ZHEN Bing-xue, GUO Zi-ru. Mechanism and optimization analysis of heat-resistant explosives[J].EngineeringBlasting, 2014, 20(5): 59-63.

    [8] 沈兆武, 馬宏昊. 小直徑低線密度裝藥的導(dǎo)爆索或雷管延期線的制作方法,中國: CN,CN101556130 A[P]. 2009.

    [9] Lee E L, Tarver C M. Phenomenological model of shock initiation in heterogeneous explosives[J].PhysicsofFluids(1958-1988), 1980, 23(12): 2362-2372.

    [10] 奧爾連科.ЛЛ.爆炸物理學(xué)[M].孫承緯,譯.北京:科學(xué)出版社, 2011: 1038-1039

    [11] Tarver C M, Hallquist J O, Erickson L M. Modeling short pulse duration shock initiation of solid explosives[C]∥Presented at the 8th Intern. Symp. on Detonation, Albuquerque, N. Mex., 15 Jul. 1985.

    [12] White B W, Sullivan K T, Gash A E, et al. Modeling the Detonation Wave Dynamics in Reactive Materials[R]. Lawrence Livermore National Laboratory (LLNL), Livermore, CA, 2016.

    [13] 李茂生, 陳棟泉. 高溫高壓下材料的本構(gòu)模型[J]. 高壓物理學(xué)報(bào), 2001, 15(1): 23-31.

    LI Mao-sheng, CHEN Dong-quan. A constitutive model for materials under high-temperature and pressure[J].ChineseJournalofHighPressurePhysics, 2001, 15(1): 23-31.

    [14] 馬宏昊, 沈兆物, 陳文川, 等. PETN 和 RDX 在無起爆藥雷管中作為激發(fā)藥的對比性研究[J]. 含能材料, 2008, 16(3): 285-289.

    MA Hong-hao, SHEN Zhao-wu, CHEN Wen-chuan, et al. Research on the performance of PETN and RDX as the excitation powder[J].ChineseJournalofEnergeticMaterials(HannengCailiao), 2008, 16(3): 285-289.

    [15] 張寶坪, 等, 爆轟物理學(xué)[M]. 北京: 兵器工業(yè)出版社, 2001: 291-296.

    [16] 馬宏昊. 高安全雷管機(jī)理與應(yīng)用的研究[D]. 合肥: 中國科學(xué)技術(shù)大學(xué), 2008.

    MA Hong-hao. Research on the principle and application of the high-safed etonator[D]. Hefei: University of science and technology of China, 2008.

    欧美激情极品国产一区二区三区| 99久久99久久久精品蜜桃| 国产无遮挡羞羞视频在线观看| 亚洲成人免费av在线播放| 91精品三级在线观看| 永久免费av网站大全| 免费不卡黄色视频| 午夜91福利影院| 欧美日韩黄片免| 国产欧美日韩一区二区三 | 国产免费福利视频在线观看| 中文字幕色久视频| 9191精品国产免费久久| 黑人巨大精品欧美一区二区蜜桃| 精品久久久久久久毛片微露脸 | 午夜av观看不卡| 国产亚洲欧美精品永久| 我要看黄色一级片免费的| 黑人巨大精品欧美一区二区蜜桃| 国产熟女欧美一区二区| 精品久久蜜臀av无| 亚洲欧美中文字幕日韩二区| 午夜免费观看性视频| 亚洲欧美一区二区三区黑人| 久久久国产一区二区| 国产精品久久久久成人av| 亚洲伊人色综图| 性高湖久久久久久久久免费观看| 国产成人91sexporn| 中文字幕av电影在线播放| 亚洲av综合色区一区| 久久av网站| 中国美女看黄片| 一级片'在线观看视频| 高清欧美精品videossex| 精品人妻1区二区| 夫妻性生交免费视频一级片| 国产精品亚洲av一区麻豆| 99九九在线精品视频| 国产爽快片一区二区三区| 99re6热这里在线精品视频| 18在线观看网站| www.精华液| 曰老女人黄片| 久久久欧美国产精品| 老鸭窝网址在线观看| 亚洲精品乱久久久久久| 国产精品99久久99久久久不卡| 青春草视频在线免费观看| 亚洲av日韩在线播放| 中文字幕色久视频| 久久久久国产精品人妻一区二区| 尾随美女入室| 免费人妻精品一区二区三区视频| 亚洲精品国产色婷婷电影| 亚洲自偷自拍图片 自拍| 久久久国产欧美日韩av| 天天操日日干夜夜撸| 又粗又硬又长又爽又黄的视频| 99国产精品一区二区三区| 国产精品久久久av美女十八| 中文字幕人妻丝袜制服| 久久久精品94久久精品| 免费在线观看影片大全网站 | 国产精品 欧美亚洲| 9热在线视频观看99| 国产淫语在线视频| 国产精品偷伦视频观看了| 一级毛片电影观看| 晚上一个人看的免费电影| 亚洲精品美女久久久久99蜜臀 | 人人妻,人人澡人人爽秒播 | 国语对白做爰xxxⅹ性视频网站| 精品亚洲成国产av| 亚洲av欧美aⅴ国产| √禁漫天堂资源中文www| 黄色视频在线播放观看不卡| 中文字幕人妻丝袜一区二区| 精品国产一区二区久久| 亚洲成人手机| 精品一区二区三卡| 欧美精品一区二区免费开放| 亚洲国产最新在线播放| 91成人精品电影| 欧美激情高清一区二区三区| 国产精品欧美亚洲77777| 97人妻天天添夜夜摸| 午夜免费观看性视频| 亚洲精品日韩在线中文字幕| 国精品久久久久久国模美| 久久久久久久国产电影| 久久人妻福利社区极品人妻图片 | 国产不卡av网站在线观看| 十八禁人妻一区二区| 两个人看的免费小视频| 成人国产一区最新在线观看 | 国产精品免费视频内射| 一二三四社区在线视频社区8| 成人黄色视频免费在线看| 啦啦啦中文免费视频观看日本| 久久人人97超碰香蕉20202| 国产精品一区二区在线观看99| 精品久久久久久电影网| 9热在线视频观看99| 国产又色又爽无遮挡免| 香蕉国产在线看| 亚洲精品av麻豆狂野| 日韩制服丝袜自拍偷拍| 久久99一区二区三区| 午夜久久久在线观看| 久久人人爽人人片av| 人妻一区二区av| 免费观看a级毛片全部| 久久国产亚洲av麻豆专区| 乱人伦中国视频| 性色av乱码一区二区三区2| 国产男人的电影天堂91| 韩国高清视频一区二区三区| 国产男女超爽视频在线观看| 国产激情久久老熟女| 亚洲精品第二区| 人人妻人人爽人人添夜夜欢视频| 亚洲av男天堂| 亚洲欧美激情在线| 人人妻人人澡人人爽人人夜夜| 久久久久久亚洲精品国产蜜桃av| 亚洲欧美一区二区三区国产| 国产一区二区三区综合在线观看| 精品欧美一区二区三区在线| 美女视频免费永久观看网站| 国产精品一区二区在线观看99| 国产精品人妻久久久影院| 久久国产精品影院| 久久鲁丝午夜福利片| 日本一区二区免费在线视频| 下体分泌物呈黄色| www.av在线官网国产| 国产成人91sexporn| 国产片特级美女逼逼视频| 1024香蕉在线观看| 色婷婷av一区二区三区视频| 狂野欧美激情性bbbbbb| 视频区欧美日本亚洲| 色婷婷av一区二区三区视频| 黄片播放在线免费| 一级a爱视频在线免费观看| 免费黄频网站在线观看国产| 超色免费av| 免费一级毛片在线播放高清视频 | 成年av动漫网址| 欧美日韩亚洲综合一区二区三区_| 成人三级做爰电影| 少妇粗大呻吟视频| 伊人亚洲综合成人网| 精品一区二区三区av网在线观看 | 婷婷色综合www| 国产av一区二区精品久久| 成人黄色视频免费在线看| 9191精品国产免费久久| 亚洲国产精品999| 亚洲av成人精品一二三区| 在线 av 中文字幕| 成年人免费黄色播放视频| 桃花免费在线播放| 热re99久久精品国产66热6| 国产精品一区二区在线观看99| 一级片'在线观看视频| 丰满饥渴人妻一区二区三| 亚洲av日韩精品久久久久久密 | 精品亚洲乱码少妇综合久久| 十八禁人妻一区二区| 中文字幕高清在线视频| 免费看不卡的av| 99香蕉大伊视频| 日韩免费高清中文字幕av| 80岁老熟妇乱子伦牲交| 国产成人av教育| 老司机影院成人| 男女床上黄色一级片免费看| 高清av免费在线| 天堂俺去俺来也www色官网| 国产免费现黄频在线看| 男女高潮啪啪啪动态图| 涩涩av久久男人的天堂| 成人黄色视频免费在线看| 人妻人人澡人人爽人人| 精品少妇黑人巨大在线播放| 久久久国产一区二区| 999精品在线视频| 国产有黄有色有爽视频| 高潮久久久久久久久久久不卡| 人人妻,人人澡人人爽秒播 | 亚洲欧美日韩另类电影网站| 国产成人免费无遮挡视频| 午夜免费成人在线视频| 亚洲激情五月婷婷啪啪| 午夜福利视频精品| 国产伦理片在线播放av一区| av网站在线播放免费| 老汉色∧v一级毛片| 欧美激情高清一区二区三区| 日日摸夜夜添夜夜爱| 日韩电影二区| 日本vs欧美在线观看视频| 亚洲精品成人av观看孕妇| 久久国产精品大桥未久av| 免费av中文字幕在线| 成年av动漫网址| 日韩人妻精品一区2区三区| av在线播放精品| 国产淫语在线视频| 国产成人免费观看mmmm| 纯流量卡能插随身wifi吗| 国产成人精品无人区| 免费在线观看日本一区| 久久天堂一区二区三区四区| 狂野欧美激情性xxxx| 激情视频va一区二区三区| 最近中文字幕2019免费版| 国产精品一区二区在线不卡| 欧美在线黄色| bbb黄色大片| 少妇 在线观看| 午夜福利在线免费观看网站| 欧美久久黑人一区二区| 欧美av亚洲av综合av国产av| 欧美日韩视频精品一区| 老司机在亚洲福利影院| 精品国产乱码久久久久久男人| av天堂久久9| 丝袜喷水一区| 又大又黄又爽视频免费| 亚洲av综合色区一区| 在线看a的网站| 啦啦啦在线免费观看视频4| 日韩电影二区| 爱豆传媒免费全集在线观看| 欧美日韩黄片免| 99国产精品一区二区蜜桃av | 男女之事视频高清在线观看 | 18禁国产床啪视频网站| 可以免费在线观看a视频的电影网站| 国产成人欧美| 日韩大码丰满熟妇| 亚洲精品国产av成人精品| 人人妻人人爽人人添夜夜欢视频| 男女国产视频网站| 久久精品熟女亚洲av麻豆精品| 亚洲欧洲国产日韩| 欧美 日韩 精品 国产| 精品福利永久在线观看| 国产免费视频播放在线视频| 国产成人一区二区在线| 伦理电影免费视频| 免费黄频网站在线观看国产| 精品久久久久久电影网| 91字幕亚洲| 国产黄色免费在线视频| 无限看片的www在线观看| 在线天堂中文资源库| www.精华液| 男女午夜视频在线观看| 两性夫妻黄色片| a 毛片基地| 天天添夜夜摸| 精品少妇内射三级| 在线天堂中文资源库| 国产视频首页在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 一区二区av电影网| 男人舔女人的私密视频| 大香蕉久久成人网| 又紧又爽又黄一区二区| 欧美精品一区二区大全| 不卡av一区二区三区| 性高湖久久久久久久久免费观看| 国产av国产精品国产| 赤兔流量卡办理| 亚洲精品乱久久久久久| 国产爽快片一区二区三区| 大香蕉久久成人网| 777久久人妻少妇嫩草av网站| 叶爱在线成人免费视频播放| 99国产精品99久久久久| 黄色视频不卡| 亚洲欧美日韩高清在线视频 | 国产一区有黄有色的免费视频| 别揉我奶头~嗯~啊~动态视频 | 亚洲精品一二三| 国产一区亚洲一区在线观看| 久久99一区二区三区| 啦啦啦视频在线资源免费观看| 宅男免费午夜| 男女边吃奶边做爰视频| 午夜福利乱码中文字幕| 午夜久久久在线观看| 一个人免费看片子| 在线av久久热| 日本91视频免费播放| 欧美日韩福利视频一区二区| 久久久久久人人人人人| 国产色视频综合| 午夜激情av网站| 男女之事视频高清在线观看 | 日韩,欧美,国产一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 日韩视频在线欧美| 两人在一起打扑克的视频| 午夜视频精品福利| 亚洲人成电影免费在线| 午夜精品国产一区二区电影| 日本午夜av视频| 麻豆乱淫一区二区| 视频在线观看一区二区三区| 97人妻天天添夜夜摸| 男女之事视频高清在线观看 | www.av在线官网国产| 黑人欧美特级aaaaaa片| 亚洲精品一二三| 国产精品欧美亚洲77777| 免费少妇av软件| 欧美黄色淫秽网站| 免费黄频网站在线观看国产| 各种免费的搞黄视频| 女人高潮潮喷娇喘18禁视频| 日韩制服骚丝袜av| 亚洲专区中文字幕在线| 99热国产这里只有精品6| 婷婷色av中文字幕| 亚洲色图 男人天堂 中文字幕| 亚洲成人免费电影在线观看 | 中文字幕亚洲精品专区| 亚洲人成电影免费在线| 亚洲专区中文字幕在线| 久久这里只有精品19| 欧美日韩亚洲高清精品| 99香蕉大伊视频| 蜜桃国产av成人99| 免费观看a级毛片全部| 欧美日韩黄片免| 只有这里有精品99| 一区二区三区激情视频| 久久性视频一级片| 激情视频va一区二区三区| 婷婷色麻豆天堂久久| 国产精品99久久99久久久不卡| 交换朋友夫妻互换小说| 久热爱精品视频在线9| 欧美在线一区亚洲| 97精品久久久久久久久久精品| 欧美精品人与动牲交sv欧美| 国产精品99久久99久久久不卡| 纵有疾风起免费观看全集完整版| 手机成人av网站| 国产欧美日韩精品亚洲av| 一级毛片电影观看| 亚洲国产看品久久| 午夜激情久久久久久久| 久久久精品免费免费高清| 美女扒开内裤让男人捅视频| 国产精品久久久人人做人人爽| 熟女少妇亚洲综合色aaa.| 国产成人系列免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 丝袜喷水一区| 大片电影免费在线观看免费| 久久精品aⅴ一区二区三区四区| 色综合欧美亚洲国产小说| 久久国产精品男人的天堂亚洲| 亚洲成国产人片在线观看| 无遮挡黄片免费观看| 久久久久久亚洲精品国产蜜桃av| 国产又爽黄色视频| 高清av免费在线| 成年动漫av网址| 国产亚洲午夜精品一区二区久久| 欧美变态另类bdsm刘玥| 亚洲情色 制服丝袜| 亚洲精品一二三| 国产一区二区 视频在线| 亚洲国产精品999| 青青草视频在线视频观看| 亚洲视频免费观看视频| 日本av手机在线免费观看| 捣出白浆h1v1| 99久久综合免费| 99精国产麻豆久久婷婷| 在现免费观看毛片| 亚洲激情五月婷婷啪啪| 国产无遮挡羞羞视频在线观看| 如日韩欧美国产精品一区二区三区| 久久久精品区二区三区| 免费在线观看黄色视频的| 久久人人97超碰香蕉20202| 精品一区二区三区四区五区乱码 | 别揉我奶头~嗯~啊~动态视频 | 亚洲人成77777在线视频| netflix在线观看网站| 天天添夜夜摸| 五月天丁香电影| 国产福利在线免费观看视频| 亚洲五月色婷婷综合| 最新的欧美精品一区二区| 在线观看www视频免费| 一区在线观看完整版| 久久久久久久国产电影| 国产97色在线日韩免费| av欧美777| 国产高清videossex| 亚洲精品一二三| 欧美日韩视频高清一区二区三区二| 99国产精品免费福利视频| 久久青草综合色| av有码第一页| 欧美人与性动交α欧美软件| 亚洲国产欧美一区二区综合| 热re99久久精品国产66热6| 最黄视频免费看| av线在线观看网站| 中文精品一卡2卡3卡4更新| avwww免费| 国产高清videossex| 欧美黑人欧美精品刺激| 黄色毛片三级朝国网站| 中文字幕另类日韩欧美亚洲嫩草| 九色亚洲精品在线播放| 欧美成狂野欧美在线观看| 狠狠精品人妻久久久久久综合| 亚洲av成人精品一二三区| 日本91视频免费播放| 熟女少妇亚洲综合色aaa.| 国产欧美日韩一区二区三 | 国产在视频线精品| 国产精品 国内视频| 久久精品久久久久久久性| 波多野结衣一区麻豆| 久久精品国产a三级三级三级| 首页视频小说图片口味搜索 | xxxhd国产人妻xxx| 成在线人永久免费视频| 五月天丁香电影| 丝袜人妻中文字幕| 亚洲专区中文字幕在线| 国产一区亚洲一区在线观看| 国产欧美日韩一区二区三 | 老司机亚洲免费影院| 久久久亚洲精品成人影院| 妹子高潮喷水视频| 菩萨蛮人人尽说江南好唐韦庄| 少妇精品久久久久久久| 欧美少妇被猛烈插入视频| 青青草视频在线视频观看| 国产伦理片在线播放av一区| 精品少妇内射三级| 天天操日日干夜夜撸| 日韩中文字幕欧美一区二区 | 亚洲成人免费av在线播放| 午夜久久久在线观看| 深夜精品福利| 午夜日韩欧美国产| av网站免费在线观看视频| 激情视频va一区二区三区| 人人妻人人澡人人看| 午夜福利视频在线观看免费| 激情五月婷婷亚洲| 黄色视频不卡| 在现免费观看毛片| 国产色视频综合| 久久精品国产a三级三级三级| 婷婷丁香在线五月| 深夜精品福利| 午夜福利影视在线免费观看| 日韩av不卡免费在线播放| 妹子高潮喷水视频| 十分钟在线观看高清视频www| 久热这里只有精品99| 一级毛片女人18水好多 | 久久久久国产精品人妻一区二区| 精品熟女少妇八av免费久了| 搡老岳熟女国产| 制服人妻中文乱码| 免费高清在线观看日韩| 亚洲欧美一区二区三区国产| 欧美日韩福利视频一区二区| 欧美在线黄色| 精品少妇久久久久久888优播| 人体艺术视频欧美日本| 一级黄色大片毛片| 一区二区av电影网| 人妻 亚洲 视频| 亚洲,一卡二卡三卡| 天堂中文最新版在线下载| 侵犯人妻中文字幕一二三四区| 亚洲人成电影免费在线| 日韩 欧美 亚洲 中文字幕| 黄片播放在线免费| 一级毛片黄色毛片免费观看视频| 久久国产精品影院| 国产精品秋霞免费鲁丝片| av欧美777| av国产精品久久久久影院| 免费观看人在逋| 国产日韩欧美视频二区| 免费不卡黄色视频| 国产一卡二卡三卡精品| 久久精品国产综合久久久| 色精品久久人妻99蜜桃| 91字幕亚洲| 2018国产大陆天天弄谢| 黄色视频不卡| 手机成人av网站| 99香蕉大伊视频| 欧美在线黄色| 亚洲 欧美一区二区三区| 性色av乱码一区二区三区2| 一区二区三区激情视频| 日韩大片免费观看网站| 亚洲av美国av| 又黄又粗又硬又大视频| 欧美成狂野欧美在线观看| 无遮挡黄片免费观看| 丰满迷人的少妇在线观看| 叶爱在线成人免费视频播放| 建设人人有责人人尽责人人享有的| 精品久久久久久电影网| 久久天堂一区二区三区四区| 亚洲国产最新在线播放| 免费高清在线观看视频在线观看| 午夜免费成人在线视频| 亚洲国产精品一区三区| 精品国产乱码久久久久久男人| 亚洲熟女毛片儿| 菩萨蛮人人尽说江南好唐韦庄| 国产精品久久久人人做人人爽| 热re99久久国产66热| 美女大奶头黄色视频| 亚洲精品国产av蜜桃| 国产精品三级大全| 五月开心婷婷网| 国产精品99久久99久久久不卡| 视频区图区小说| 欧美变态另类bdsm刘玥| 午夜91福利影院| 亚洲五月婷婷丁香| 国产精品久久久av美女十八| 久久青草综合色| 人成视频在线观看免费观看| 亚洲国产成人一精品久久久| www.精华液| 人人妻人人澡人人爽人人夜夜| 一级黄色大片毛片| 亚洲免费av在线视频| 深夜精品福利| 五月天丁香电影| av福利片在线| 考比视频在线观看| 天天躁日日躁夜夜躁夜夜| 国产精品成人在线| 免费少妇av软件| 啦啦啦在线免费观看视频4| 母亲3免费完整高清在线观看| 成人国产一区最新在线观看 | 90打野战视频偷拍视频| 又紧又爽又黄一区二区| 首页视频小说图片口味搜索 | 免费在线观看视频国产中文字幕亚洲 | 国产高清不卡午夜福利| 国产亚洲av高清不卡| 首页视频小说图片口味搜索 | 97精品久久久久久久久久精品| 极品人妻少妇av视频| 亚洲av美国av| 亚洲一区中文字幕在线| 视频区图区小说| 18禁国产床啪视频网站| 不卡av一区二区三区| 日韩伦理黄色片| 国产午夜精品一二区理论片| 亚洲国产最新在线播放| 亚洲色图 男人天堂 中文字幕| 日韩电影二区| 好男人视频免费观看在线| 欧美日韩成人在线一区二区| 免费少妇av软件| 99国产精品99久久久久| 嫁个100分男人电影在线观看 | 人人妻人人添人人爽欧美一区卜| 国产福利在线免费观看视频| 黑人巨大精品欧美一区二区蜜桃| 亚洲国产精品一区三区| 久久99一区二区三区| 韩国精品一区二区三区| 两个人看的免费小视频| 99国产精品一区二区蜜桃av | 成人三级做爰电影| 久久精品人人爽人人爽视色| 建设人人有责人人尽责人人享有的| xxx大片免费视频| 国产三级黄色录像| 国产精品国产三级国产专区5o| 亚洲成国产人片在线观看| 欧美激情极品国产一区二区三区| 侵犯人妻中文字幕一二三四区| 一本色道久久久久久精品综合| 岛国毛片在线播放| 欧美成人午夜精品| 色综合欧美亚洲国产小说| 久久女婷五月综合色啪小说| 亚洲中文字幕日韩| 亚洲精品国产一区二区精华液| 免费高清在线观看日韩| 国产成人欧美| 亚洲欧洲精品一区二区精品久久久| 免费观看人在逋| 亚洲 国产 在线|