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

    水中爆炸作用機(jī)理及毀傷效應(yīng)研究綜述

    2016-01-04 12:16:35莊鐵栓
    火炸藥學(xué)報(bào) 2016年1期
    關(guān)鍵詞:作用機(jī)理綜述

    伍 俊,楊 益,莊鐵栓,2

    (1. 總參工程兵科研三所,河南 洛陽(yáng) 471023; 2. 解放軍理工大學(xué)爆炸沖擊防災(zāi)減災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210007)

    ?

    水中爆炸作用機(jī)理及毀傷效應(yīng)研究綜述

    伍俊1,楊益1,莊鐵栓1,2

    (1. 總參工程兵科研三所,河南 洛陽(yáng) 471023; 2. 解放軍理工大學(xué)爆炸沖擊防災(zāi)減災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210007)

    摘要:為深入揭示水中爆炸的作用機(jī)理及毀傷效應(yīng)特性,從不同炸藥水中爆炸機(jī)理、水中爆炸沖擊波傳播、水中爆炸氣泡脈動(dòng)、水中爆炸結(jié)構(gòu)破壞效應(yīng)與動(dòng)態(tài)響應(yīng)、水中爆炸測(cè)試試驗(yàn)技術(shù)、水中爆炸數(shù)值模擬與毀傷評(píng)估等6個(gè)方面,綜述了國(guó)內(nèi)外對(duì)水中爆炸作用機(jī)理及毀傷效應(yīng)的研究進(jìn)展。提出急需解決的關(guān)鍵技術(shù)問(wèn)題為:非理想炸藥水中爆炸能量釋放與傳播、沖擊波/氣泡耦合作用對(duì)結(jié)構(gòu)的毀傷、水中爆炸結(jié)構(gòu)毀傷評(píng)估方法、水中爆炸多尺度模擬技術(shù)等。附參考文獻(xiàn)86篇。

    關(guān)鍵詞:水中爆炸;作用機(jī)理;毀傷效應(yīng);氣泡脈動(dòng);非理想炸藥;毀傷評(píng)估;綜述

    引言

    水與空氣兩種介質(zhì)的物理性質(zhì)存在較大差異,導(dǎo)致炸藥水中爆炸與空氣中的爆炸特性有許多不同之處,不僅爆炸的物理現(xiàn)象和機(jī)理不同,爆炸荷載作用也各有特點(diǎn)。對(duì)水中爆炸的研究始于19世紀(jì)初,1860年,美國(guó)的Albort將軍指導(dǎo)進(jìn)行了早期艦船爆炸試驗(yàn),并于1881年發(fā)表了對(duì)試驗(yàn)情況的總結(jié)報(bào)告[1]。水中爆炸研究在兩次世界大戰(zhàn)中得到了快速發(fā)展,在此期間,美國(guó)Cole[2-5]根據(jù)大量的水中爆炸試驗(yàn),于1948年對(duì)水中爆炸的物理現(xiàn)象、爆炸機(jī)理和水中爆炸遵循的基本定律進(jìn)行了研究和總結(jié),使水中爆炸研究形成了比較系統(tǒng)的理論體系。在之后的幾十年間,人們?cè)谒斜ɑA(chǔ)理論研究方面并未取得突破性進(jìn)展,僅局限于對(duì)Cole理論的修正和補(bǔ)充。20世紀(jì)70年代以后,隨著計(jì)算機(jī)技術(shù)和計(jì)算方法的發(fā)展進(jìn)步,水中爆炸數(shù)值模擬推動(dòng)了水中爆炸研究的快速發(fā)展,當(dāng)前的水中爆炸數(shù)值模擬結(jié)果可以指導(dǎo)并應(yīng)用于工程實(shí)踐[6-7]。

    水中爆炸研究主要關(guān)注爆炸在水中產(chǎn)生、演變、傳播以及與水下目標(biāo)相互作用的理論、方法和手段,涉及兵器科學(xué)與技術(shù)、船舶與海洋工程、力學(xué)、材料科學(xué)等多個(gè)學(xué)科,現(xiàn)已成為支撐水中兵器及其對(duì)目標(biāo)毀傷能力、水中工程技術(shù)創(chuàng)新發(fā)展的關(guān)鍵基礎(chǔ)。國(guó)內(nèi)外也有多個(gè)專職研究機(jī)構(gòu)圍繞水中爆炸現(xiàn)象、傳播規(guī)律、結(jié)構(gòu)破壞效應(yīng)、試驗(yàn)技術(shù)及應(yīng)用等方面開(kāi)展了大量的研究工作,使水中爆炸研究得到了迅猛發(fā)展,本文對(duì)此進(jìn)行總結(jié)及評(píng)述。

    1不同炸藥水中爆炸機(jī)理研究

    水中爆炸及其沖擊波物理學(xué)現(xiàn)象十分復(fù)雜,只有掌握不同炸藥水中爆炸現(xiàn)象的機(jī)理才能正確理解動(dòng)態(tài)沖擊載荷在水中的發(fā)展和傳播。張志江等[8]對(duì)高能炸藥水中爆炸的能量輸出特性進(jìn)行了數(shù)值模擬研究,采用位標(biāo)函數(shù)捕捉物質(zhì)界面,用修正的虛擬流體方法處理界面問(wèn)題,其中含能材料爆轟產(chǎn)物JWL狀態(tài)方程系數(shù)通過(guò)Kihara-Hikita-Tanaka計(jì)算的等熵膨脹數(shù)據(jù)擬合得到,并實(shí)現(xiàn)了程序化,通過(guò)對(duì)TNT和PETN兩種高能炸藥在水中爆炸的能量輸出特性進(jìn)行計(jì)算,結(jié)果表明,計(jì)算結(jié)果與試驗(yàn)相似率計(jì)算結(jié)果比較吻合,誤差在10%以內(nèi)。

    牛國(guó)濤等[9]對(duì)含納米鋁和微米鋁炸藥進(jìn)行了水下爆炸試驗(yàn),對(duì)比分析了二者在水下爆炸的沖擊波峰值壓力、沖擊波能、氣泡能。結(jié)果表明,在以RDX為基的含鋁炸藥水下爆炸試驗(yàn)中,微米鋁質(zhì)量分?jǐn)?shù)在25%左右時(shí),達(dá)到最大沖擊波壓力值;含納米鋁的炸藥沖擊波壓力隨鋁粉含量呈線性變化,且最大值高于含微米鋁炸藥,在水下炸藥納米鋁配方設(shè)計(jì)中,納米鋁粉質(zhì)量分?jǐn)?shù)不大于20%,可得到較高的沖擊波能,納米鋁和微米鋁的級(jí)配可以明顯提高沖擊波壓力,但是含微米鋁炸藥沖擊波能始終占優(yōu)。史銳等[10]對(duì)TNT和3種新設(shè)計(jì)的含鋁炸藥進(jìn)行了水中爆炸實(shí)驗(yàn),比較了各炸藥的爆炸性能。結(jié)果表明,沖擊波峰值超壓、比例沖量和比例能流密度等參數(shù)能較好地符合爆炸相似率;盡管含鋁炸藥的沖擊波能隨距離變化存在跳動(dòng)性,但仍可以發(fā)現(xiàn),設(shè)計(jì)的3種PBX炸藥的沖擊波能和氣泡能較理想炸藥都有大幅提高。

    馮曉軍等[11]以RDX為主炸藥,通過(guò)添加不同含量的Al粉和AP制成6種RDX基復(fù)合炸藥,采用水中爆炸實(shí)驗(yàn),從沖擊波參數(shù)和氣泡參數(shù)等方面分析了Al粉和AP對(duì)復(fù)合炸藥水中爆炸性能的影響,結(jié)果表明,主炸藥含量不變時(shí)隨著Al與AP摩爾比的增大,沖擊波峰值壓力、時(shí)間常數(shù)、沖量和能流密度都逐漸減小,當(dāng)Al與AP的摩爾比約為3.8時(shí),氣泡半徑或氣泡脈動(dòng)周期達(dá)到最大值。宋浦等[12]通過(guò)小型爆炸水池實(shí)驗(yàn),研究了端面點(diǎn)起爆和中心點(diǎn)起爆狀態(tài)下TNT裝藥水中爆炸能量輸出結(jié)構(gòu)的變化規(guī)律,并將兩者進(jìn)行了對(duì)比,結(jié)果表明,裝藥端面點(diǎn)起爆后,能量輸出結(jié)構(gòu)-沖擊波壓力峰值和比沖量有較大的變化,靠近起爆端方向的沖擊波壓力峰值和比沖量變化均較其他方位提高,端面點(diǎn)起爆狀態(tài)比中心點(diǎn)起爆特定方位的沖擊波壓力峰值和比沖量增大約10%。得出在裝藥形狀基本不變的條件下,改變起爆方式即可實(shí)現(xiàn)在水下特定方位處的爆炸能量輸出結(jié)構(gòu)變化。

    周霖等[13]通過(guò)典型炸藥爆炸的性能參數(shù)分析,計(jì)算了不同類型水中兵器用混合炸藥水中爆炸的沖擊波能、氣泡能、沖擊波峰值壓力方程中的系數(shù)k和α以及水中爆炸沖擊波峰值壓力,研究表明,鋁氧比在0.38左右時(shí),炸藥水中爆炸的沖擊波能達(dá)到最大,氣泡能則隨著鋁氧比的增加而增大,此時(shí)提高炸藥爆熱,可獲得較大的沖擊波能,從而得到較大的作用沖量,實(shí)現(xiàn)對(duì)目標(biāo)的有效毀傷;鋁氧比接近1,即鋁粉質(zhì)量分?jǐn)?shù)約40%時(shí),氣泡能達(dá)到最大值,再增加鋁粉,氣泡會(huì)降低;炸藥水中爆炸的沖擊波峰值壓力依賴于超壓方程中系數(shù)k和α,k越大α越小,炸藥水中爆炸沖擊波峰值壓力就越高;爆熱和密度的乘積決定k值大小,爆熱和密度乘積越大,k值就越大,反之k值越小;沖擊波損失系數(shù)、爆熱和爆壓影響α值大小,爆熱越小、爆壓越大、沖擊波損失系數(shù)越大,α值就越大,反之α值就越小。

    閆偉杰等[14]通過(guò)測(cè)定某PBX型高能炸藥的水中爆炸能量輸出參數(shù),給出了水中兵器用該P(yáng)BX型炸藥的水下爆炸能量輸出特性;對(duì)大藥量帶殼PBX炸藥水中爆炸沖擊波的遠(yuǎn)場(chǎng)壓力進(jìn)行了測(cè)量,探討了鋼質(zhì)殼體對(duì)水中爆炸的沖擊波能及氣泡能產(chǎn)生的影響,并用數(shù)值計(jì)算結(jié)果進(jìn)行了驗(yàn)證。結(jié)果表明,同種炸藥質(zhì)量的裝藥,帶鋼殼炸藥所產(chǎn)生的沖擊波峰值壓力在近場(chǎng)大于裸炸藥,而在遠(yuǎn)場(chǎng)兩者的峰值壓力的差距逐漸變?。粠т撡|(zhì)殼體炸藥所產(chǎn)生的沖擊波能要大于裸炸藥,而氣泡能則小于裸炸藥。

    2水中爆炸沖擊波傳播研究

    水的密度比空氣的大很多,因此爆炸產(chǎn)物的膨脹在水中比在空氣中慢很多。因?yàn)樗拿芏却?,可壓縮性小,這就使得水下沖擊波的初始?jí)毫Ρ瓤諝庵写蟮枚?,由于水的特殊性質(zhì),開(kāi)展水中爆炸沖擊波傳播規(guī)律研究對(duì)深化毀傷效應(yīng)認(rèn)識(shí)具有重要作用。趙根等[15]對(duì)水中沖擊波傳播規(guī)律的公式形式進(jìn)行了探討,通過(guò)對(duì)4次水中沖擊爆破試驗(yàn)及對(duì)監(jiān)測(cè)資料的分析,得到了高能、普通乳化炸藥水中爆炸沖擊波的傳播規(guī)律。結(jié)果表明,相同藥量條件下,炸藥水中爆炸產(chǎn)生的沖擊波大小與炸藥類型有關(guān),大小關(guān)系是:TNT炸藥>高能乳化炸藥>普通乳化炸藥。李金河等[16]采用PCB138壓力傳感器測(cè)量了TNT、RS211、HLZY-1和HLZY-3等幾種炸藥水中爆炸沖擊波遠(yuǎn)場(chǎng)的壓力—時(shí)間歷程,計(jì)算得到這幾種炸藥水中爆炸沖擊波性能參數(shù)及其相似常數(shù)。結(jié)果表明,含鋁炸藥水中爆炸沖擊波遠(yuǎn)場(chǎng)的傳播服從指數(shù)變化的相似率,含鋁炸藥的水中爆炸沖擊波時(shí)間常數(shù)大,比沖量和比能流密度較高,其沖擊波性能比TNT優(yōu)越。劉志[17]通過(guò)試驗(yàn)研究了水中爆炸沖擊波的傳播特性,通過(guò)現(xiàn)場(chǎng)試爆及其現(xiàn)象分析,得出了藥包布置原則;利用高速攝影技術(shù)再現(xiàn)了水中爆炸沖擊波波陣面的動(dòng)態(tài)傳播過(guò)程,得出波陣面?zhèn)鞑ニ俣燃捌鋫鞑ヒ?guī)律;根據(jù)水沖擊波波陣面?zhèn)鞑ニ俣龋贸霾煌嚯x處的峰值壓力,并對(duì)水中沖擊波峰值壓力、傳播距離及藥量關(guān)系進(jìn)行分析。

    楊莉等[18]研究了沉底裝藥水中爆炸沖擊波的傳播規(guī)律,結(jié)果表明,自由場(chǎng)和沉底爆炸沖擊波壓力均呈現(xiàn)出指數(shù)衰減規(guī)律,但自由場(chǎng)水中爆炸由于氣泡脈動(dòng)作用存在明顯的二次壓力波,其峰值壓力約為沖擊波峰值壓力的1/4,而沉底爆炸則無(wú)明顯的二次壓力波情況,但在裝藥正上方位置處較自由場(chǎng)爆炸出現(xiàn)顯著波動(dòng)。項(xiàng)大林等[19]為了研究裝藥殼體厚度對(duì)水中爆炸沖擊波特性的影響,對(duì)1.0 kg柱形含鋁炸藥分別在厚度為6.0 mm的鋼殼或硬鋁殼裝藥下進(jìn)行了水中爆炸實(shí)驗(yàn)與數(shù)值模擬研究,結(jié)果表明,相對(duì)硬鋁殼,鋼殼裝藥的沖擊波沖量、衰減時(shí)間常數(shù)以及沖擊波能都偏大,而沖擊波峰值相差不大;針對(duì)不同厚度鋼殼裝藥的數(shù)值模擬表明,隨著殼體厚度的增加,沖擊波參數(shù)明顯增強(qiáng),當(dāng)殼體厚度超過(guò)最優(yōu)值時(shí)這種效應(yīng)減弱;金屬殼的存在導(dǎo)致沖擊波峰值爬升產(chǎn)生滯后效應(yīng);對(duì)一定質(zhì)量的炸藥,存在可有效提高沖擊波壓力峰值的最優(yōu)殼厚,填裝比可作為衡量效果的重要指標(biāo)。

    盛振新等[20]研究了兩發(fā)炸藥水中爆炸沖擊波之間的相互作用,利用AUTODYN模擬軟件,計(jì)算了同時(shí)起爆和延時(shí)起爆兩種情況下兩發(fā)炸藥水中爆炸沖擊波峰值壓力,并與單發(fā)炸藥水中爆炸沖擊波進(jìn)行對(duì)比,通過(guò)分析峰值壓力和沖擊波到達(dá)延時(shí),得出在兩發(fā)炸藥同時(shí)起爆的情況下,炸藥爆炸產(chǎn)生的氣泡會(huì)對(duì)臨近炸藥的爆炸沖擊波產(chǎn)生衰減作用;兩發(fā)炸藥延時(shí)起爆且不會(huì)引起殉爆的情況下,在距先起爆炸藥較近處,先起爆炸藥產(chǎn)生的氣泡會(huì)對(duì)后起爆炸藥的沖擊波產(chǎn)生衰減作用,并且使沖擊波的傳播速度變慢;在距后起爆炸藥較近處,后起爆炸藥的爆炸會(huì)對(duì)先起爆炸藥的沖擊波產(chǎn)生增強(qiáng)作用,并使沖擊波的傳播速度變快。

    師華強(qiáng)等[21]通過(guò)求解球形一維流體動(dòng)力學(xué)方程,模擬了水下爆炸沖擊波的產(chǎn)生和傳播過(guò)程。研究表明,水下爆炸沖擊波近場(chǎng)壓力變化劇烈,對(duì)水的狀態(tài)方程要求高,Two-phase狀態(tài)方程不僅能描述水或水蒸氣的單一狀態(tài),還能很好地描述水和水蒸氣的混合狀態(tài),適合模擬水下爆炸沖擊波的產(chǎn)生和傳播過(guò)程。當(dāng)網(wǎng)格大小為1.5mm時(shí),計(jì)算精度滿足要求,并驗(yàn)證了沖擊波在近場(chǎng)也滿足相似定律。

    3水中爆炸氣泡脈動(dòng)研究

    由于沖擊波形成的氣泡射流速度高達(dá)每秒幾千米,對(duì)其附近結(jié)構(gòu)會(huì)造成嚴(yán)重?fù)p傷,研究水中爆炸氣泡脈動(dòng)對(duì)工程設(shè)計(jì)及應(yīng)用具有重要意義。黃超等[22]采用高速攝影技術(shù)研究了水下爆炸氣泡的動(dòng)態(tài)特性,分別從定性和定量的角度分析了水箱邊界影響產(chǎn)生的誤差,研究了小當(dāng)量PETN 的水下爆炸氣泡現(xiàn)象,結(jié)果表明,在鋼制實(shí)驗(yàn)水箱中建立的實(shí)驗(yàn)室尺度水下爆炸氣泡的實(shí)驗(yàn)方法是科學(xué)可行的,考慮到邊界對(duì)氣泡運(yùn)動(dòng)的影響隨著藥量的增大而增加,因此,實(shí)驗(yàn)時(shí)需要控制使用的藥量。汪斌等[23]在實(shí)驗(yàn)水箱中采用高速攝影技術(shù)獲得了炸藥水中爆炸氣泡脈動(dòng)過(guò)程圖像,結(jié)果表明,短長(zhǎng)徑比為1.05的PENT炸藥柱水中爆炸產(chǎn)生的氣泡形狀近似球形,氣泡脈動(dòng)周期和半徑隨炸藥質(zhì)量的增加而增大;在膨脹階段,氣泡中心位置保持不變,隨著氣泡半徑的不斷縮小,氣泡上浮越來(lái)越明顯;受重力影響,氣泡下表面收縮速度高于上表面,在收縮至氣泡半徑最小時(shí)下表面向上沖頂,氣泡潰滅形成水射流。

    張阿漫等[24]研究了自由場(chǎng)水中爆炸氣泡的物理特性,將水中爆炸氣泡運(yùn)動(dòng)階段周圍流場(chǎng)假設(shè)為無(wú)粘、無(wú)旋、不可壓縮的理想流體,運(yùn)用邊界元法模擬自由場(chǎng)中氣泡運(yùn)動(dòng);在氣泡運(yùn)動(dòng)模擬過(guò)程中引入數(shù)值光順技術(shù)及彈性網(wǎng)格技術(shù),避免因網(wǎng)格扭曲而導(dǎo)致的數(shù)值發(fā)散,并開(kāi)發(fā)了計(jì)算程序。計(jì)算值與試驗(yàn)值吻合良好,誤差小于10%。張阿漫等[25]進(jìn)一步研究了三維水下爆炸氣泡動(dòng)力學(xué),結(jié)果表明,位于目標(biāo)結(jié)構(gòu)附近的水下爆炸氣泡,在多種情況下可向結(jié)構(gòu)形成噴射,噴射會(huì)對(duì)目標(biāo)結(jié)構(gòu)產(chǎn)生嚴(yán)重破壞;在氣泡附近正確的位置,在自由表面和重力的作用下,氣泡產(chǎn)生的噴射可能直接撞擊到結(jié)構(gòu)上,或者完全錯(cuò)開(kāi),因此研究氣泡動(dòng)力學(xué)對(duì)于提高水下結(jié)構(gòu)生存能力具有重要意義。黃超等[26]研究了柱形裝藥水中爆炸氣泡的形態(tài)演變,采用高速攝影方法在實(shí)驗(yàn)水箱中獲得了長(zhǎng)徑比為3.35~6.75的柱形PETN炸藥水中爆炸氣泡脈動(dòng)的圖像,進(jìn)而結(jié)合真實(shí)的爆轟過(guò)程和Rayleigh氣泡運(yùn)動(dòng)方程,研究了氣泡的形態(tài)演變規(guī)律。結(jié)果表明,柱形裝藥在水中爆炸時(shí),形成的初始?xì)馀菪螤顬榉乔驅(qū)ΨQ形,這種非對(duì)稱特征隨著長(zhǎng)徑比的增加而增大。氣泡表面的運(yùn)動(dòng)也表現(xiàn)出明顯的非球?qū)ΨQ特征,氣泡表面在裝藥徑向的膨脹運(yùn)動(dòng)呈指數(shù)衰減,在裝藥軸向兩端的膨脹運(yùn)動(dòng)更接近于分段線性衰減;氣泡表面的不對(duì)稱運(yùn)動(dòng)與柱形裝藥水下爆炸的能量輸出結(jié)構(gòu)不均有關(guān),與沖擊波的有效能量分布規(guī)律相似;當(dāng)柱形裝藥水平安裝時(shí),水中爆炸氣泡形成偏射流,射流的方向指向起爆端上方,而垂直安裝時(shí)氣泡射流與球形裝藥的類似,在氣泡底部中心形成,方向垂直向上,但射流速度可能更高。

    趙生偉等[27]研究了小當(dāng)量水中爆炸氣泡的脈動(dòng)現(xiàn)象,設(shè)計(jì)了用于小當(dāng)量(不大于10g TNT)水中爆炸研究的爆炸水箱,在爆炸水箱裝置中進(jìn)行了多次0.125、1.000、3.370、8.000g TNT當(dāng)量的PETN球形炸藥水中爆炸試驗(yàn),采用高速攝影系統(tǒng)獲得了氣泡脈動(dòng)圖片;試驗(yàn)結(jié)果給出了試驗(yàn)范圍內(nèi)氣泡脈動(dòng)過(guò)程中氣泡直徑、速度及加速度隨時(shí)間變化的擬合曲線,通過(guò)對(duì)氣泡脈動(dòng)周期進(jìn)行分析,對(duì)氣泡直徑進(jìn)行了歸一化處理,數(shù)據(jù)能較好地符合爆炸相似率。

    金輝等[28]采用高速攝影光學(xué)測(cè)量和自由場(chǎng)壓力電學(xué)測(cè)量?jī)煞N方法對(duì)球形裝藥的水下爆炸氣泡脈動(dòng)過(guò)程進(jìn)行實(shí)驗(yàn)研究,基于勢(shì)流理論,假設(shè)爆炸氣泡周圍是不可壓縮、無(wú)粘、無(wú)旋的流場(chǎng),采用邊界積分法對(duì)炸藥水中爆炸氣泡脈動(dòng)情況進(jìn)行數(shù)值模擬,實(shí)驗(yàn)結(jié)果和模擬結(jié)果吻合,炸藥在自由場(chǎng)工況下爆炸時(shí),氣泡脈動(dòng)規(guī)律與現(xiàn)有水下爆炸理論一致;炸藥在近水底爆炸時(shí),水底對(duì)氣泡脈動(dòng)影響較大,爆炸產(chǎn)生的氣泡能大部分被水底吸收,二次脈動(dòng)壓力較自由場(chǎng)時(shí)小,二次脈動(dòng)周期增加;近水底爆炸氣泡在邊界的作用下,在膨脹階段下表面呈扁平狀,在坍塌階段指向壁面的射流形成較早,回彈階段會(huì)形成背離壁面的反向射流。金輝等[29]針對(duì)不同邊界條件下炸藥水下爆炸的能量輸出結(jié)果,設(shè)置自由場(chǎng)、沉底、近水底和近水面等不同邊界條件的水下爆炸實(shí)驗(yàn),結(jié)果表明,炸藥在距水底、水面一定距離時(shí)爆炸,邊界條件對(duì)沖擊波影響較小,對(duì)氣泡脈動(dòng)影響較大;當(dāng)炸藥與水底距離小于1倍最大氣泡半徑處爆炸時(shí),由于水底邊界影響,不能形成完整的氣泡脈動(dòng);當(dāng)炸藥在距水面小于1倍最大氣泡半徑處爆炸時(shí),氣泡直接滋出水面,不能形成氣泡脈動(dòng)壓力;當(dāng)炸藥在距水面1~1.5倍最大氣泡半徑處爆炸時(shí),部分氣泡滋出水面,氣泡脈動(dòng)壓力較自由場(chǎng)爆炸時(shí)小;氣泡在膨脹和收縮過(guò)程中,由于氣泡不斷上升,水面對(duì)氣泡脈動(dòng)的影響范圍要大于水底。黃毅等[30]基于勢(shì)流理論,建立考慮氣泡遷移效應(yīng)、自由表面效應(yīng)的水下爆炸氣泡脈動(dòng)模型,通過(guò)不同藥量和氣泡深度的計(jì)算,對(duì)遷移效應(yīng)和自由表面效應(yīng)這兩種影響因素進(jìn)行比較,討論兩種因素對(duì)氣泡脈動(dòng)的影響,將氣泡脈動(dòng)載荷的數(shù)值計(jì)算結(jié)果與經(jīng)驗(yàn)公式進(jìn)行對(duì)比和驗(yàn)證,結(jié)果吻合良好,計(jì)算結(jié)果也表明,考慮遷移效應(yīng)、自由表面效應(yīng)的氣泡脈動(dòng)模型可以對(duì)水下爆炸氣泡脈動(dòng)載荷提供較準(zhǔn)確的預(yù)測(cè)。

    4水中爆炸結(jié)構(gòu)破壞效應(yīng)與動(dòng)態(tài)響應(yīng)研究

    近場(chǎng)水中爆炸結(jié)構(gòu)破壞效應(yīng)與動(dòng)態(tài)響應(yīng)研究,對(duì)于目標(biāo)結(jié)構(gòu)威脅分析、毀傷評(píng)估具有重要作用,也是開(kāi)展結(jié)構(gòu)防護(hù)設(shè)計(jì)的重要依據(jù),因而加強(qiáng)近場(chǎng)水中爆炸結(jié)構(gòu)破壞效應(yīng)與動(dòng)態(tài)響應(yīng)研究很有必要。顧文彬等[31]用水下爆炸沖擊分析(USA)方法,研究淺層水中爆炸沖擊波與混凝土墩的相互作用。結(jié)果表明:(1)在沖擊波作用下,迎爆面壓力響應(yīng)是入射沖擊波壓力的2倍以上,最大可達(dá)到2.26倍,在迎爆面上沿表面,由于水中沖擊波入射角增大,表面壓力響應(yīng)的增加倍數(shù)更高;(2)爆炸作用過(guò)程中,混凝土墩頂面、側(cè)面及背爆面也會(huì)受到?jīng)_擊波壓力作用,頂面及側(cè)面壓力響應(yīng)可達(dá)3MPa以上,該壓力對(duì)墩體的破壞也有貢獻(xiàn);(3)在迎爆面及墩體內(nèi)部,節(jié)點(diǎn)的加速度變化是隨機(jī)的,有時(shí)發(fā)生很大震蕩,表明在強(qiáng)沖擊波載荷作用下,結(jié)構(gòu)內(nèi)部及表面材料的受力狀態(tài)是復(fù)雜的,這是流-固耦合算法把爆炸沖擊波、流體介質(zhì)、混凝土結(jié)構(gòu)變形與運(yùn)動(dòng)作為一個(gè)系統(tǒng)來(lái)研究的結(jié)果;(4)在平行于水平面的橫截面上、沿沖擊波作用方向,混凝土材料的加速變化很大,表明沖擊波作用引起的墩體內(nèi)部壓應(yīng)力波衰減較快;(5)由于迎爆面及墩體內(nèi)部加速度的振蕩,其速度響應(yīng)也有振蕩,但振蕩頻率、幅度等變化規(guī)律不同,這是由結(jié)構(gòu)材料動(dòng)態(tài)阻尼特性決定的,因此,結(jié)構(gòu)物的動(dòng)態(tài)破壞不僅與材料的強(qiáng)度有關(guān),還與材料的微觀結(jié)構(gòu)、加載速率等因素有關(guān)。

    賈憲振等[32]研究了深水環(huán)境中水下爆炸沖擊波作用下圓柱殼動(dòng)態(tài)響應(yīng),基于通用有限元程序ABAQUS,運(yùn)用靜力分析和動(dòng)力分析相結(jié)合的方法,分析了深水中圓柱殼動(dòng)態(tài)響應(yīng)的影響因素,研究了圓柱殼的破壞隨深度而變化的規(guī)律。結(jié)果表明,在同樣的沖擊波載荷作用下,圓柱殼的毀傷程度隨著所處深度的增大而逐漸加重;爆距不變時(shí),不同爆心方位所對(duì)應(yīng)的圓柱殼的響應(yīng)不同;預(yù)應(yīng)力的存在使圓柱殼毀傷程度有所增加,并且預(yù)應(yīng)力對(duì)圓柱殼響應(yīng)的影響與爆心方位有關(guān);能夠?qū)е聢A柱殼破壞的炸藥質(zhì)量隨深度增加而迅速減??;圓柱殼在淺水中和深水中的破壞機(jī)制不同,前者是由于沖擊波的作用,而后者是沖擊波和靜水壓力聯(lián)合作用的結(jié)果。

    徐俊祥等[33]針對(duì)水下爆炸的物理過(guò)程以及由此引起的混凝土重力壩的動(dòng)力響應(yīng)的復(fù)雜性,建立了全耦合的模型,考慮了炸藥的爆炸、爆炸波的水中傳播、水介質(zhì)與壩體結(jié)構(gòu)的相互作用,對(duì)混凝土重力壩的動(dòng)力響應(yīng)進(jìn)行了模擬,結(jié)果表明,計(jì)算的爆炸荷載下水中的峰值壓力和Cole的經(jīng)驗(yàn)公式具有較好的一致性,表明全耦合分析模型是正確的;壩基加速度遠(yuǎn)小于壩體加速度值,在簡(jiǎn)化計(jì)算中忽略由于爆炸沖擊荷載引起的基礎(chǔ)震動(dòng)對(duì)結(jié)構(gòu)動(dòng)力響應(yīng)的影響是可行的;水中爆炸荷載下,壩體水平方向的加速度均大于壩體垂直方向的加速度,水平方向的作用對(duì)壩體的動(dòng)力響應(yīng)起控制作用。

    Zhu Feng等[34]利用LS-DYNA軟件對(duì)拱形壩水中爆炸響應(yīng)進(jìn)行了數(shù)值模擬研究。結(jié)果表明,水下爆炸沖擊波荷載作用下,壩體水平方向的響應(yīng)大于垂直方向的響應(yīng);響應(yīng)程度與裝藥點(diǎn)的距離有關(guān),距離裝藥點(diǎn)越遠(yuǎn),響應(yīng)程度越小。顧文彬等[35]研究了淺層水中爆炸沖擊波對(duì)混凝土墩斜碰撞作用,得出:(1)單個(gè)裝藥爆炸對(duì)混凝土墩作用時(shí),迎爆面以反射超壓破壞作用為主,背爆面以繞射作用沖量為主,JHL-3裝藥在比例距離小于5.3時(shí),迎爆面反射作用因數(shù)大于2.7,且隨比例距離的增加而減??;(2)兩個(gè)裝藥關(guān)于混凝土墩對(duì)稱設(shè)置同步起爆時(shí),混凝土迎爆面先后受到直達(dá)沖擊波、背爆面沖擊波的透射波和繞射波作用;裝藥間隔足夠大時(shí),迎爆面首先受前驅(qū)波作用;繞射波使迎爆面受到負(fù)壓拉伸作用;(3)回歸得到的迎爆面沖擊波反射超壓計(jì)算半理論半經(jīng)驗(yàn)公式,在試驗(yàn)比距離范圍內(nèi)具有足夠精度;(4)兩個(gè)裝藥對(duì)稱和非對(duì)稱設(shè)置同時(shí)爆炸作用時(shí),在兩個(gè)裝藥對(duì)混凝土墩側(cè)面和頂面的先后作用仍然有顯著的破壞疊加效應(yīng)。

    賴?guó)Q等[36]采用LS-DYNA軟件對(duì)不同加筋強(qiáng)度和不同加筋位置筋板結(jié)構(gòu)的水中接觸爆炸進(jìn)行了數(shù)值模擬,定義描述不同加筋強(qiáng)度、位置的強(qiáng)度因子和距離因子,建立不同強(qiáng)度因子的6種加筋模型和不同距離因子的3種加筋模型;通過(guò)對(duì)水中接觸爆炸下典型加筋結(jié)構(gòu)破損過(guò)程分析和不同模型縱、橫破口長(zhǎng)度和破口面積的對(duì)比,得到了加筋板在水中接觸爆炸下強(qiáng)度因子、距離因子對(duì)破口形狀和大小的影響規(guī)律。

    李裕春等[37]利用試驗(yàn)方法研究淺層水中沖擊波與混凝土結(jié)構(gòu)的相互作用,探討混凝土結(jié)構(gòu)受水中沖擊波作用的破壞機(jī)制。結(jié)果表明,淺層水中爆炸沖擊波對(duì)混凝土結(jié)構(gòu)的作用體現(xiàn)在沖擊壓縮和拉伸兩個(gè)方面。由于水介質(zhì)的流動(dòng)性和各向同性,混凝土結(jié)構(gòu)受沖擊波的圍壓作用,使混凝土材料處于多向應(yīng)力狀態(tài),多軸應(yīng)力狀態(tài)是導(dǎo)致混凝土材料失效和斷裂的最根本原因。李萬(wàn)等[38]基于某水下目標(biāo)的抗水下爆炸實(shí)驗(yàn)數(shù)據(jù),利用小波變換良好的時(shí)頻局部化性質(zhì)對(duì)監(jiān)測(cè)的水下目標(biāo)內(nèi)部裝置沖擊加速度信號(hào)進(jìn)行時(shí)頻特征分析,得到了不同頻帶上沖擊信號(hào)振動(dòng)分量的時(shí)間歷程曲線和能量分布。結(jié)果表明,基于沖擊信號(hào)小波時(shí)頻特征分析的頻帶能量可以同時(shí)反映沖擊振動(dòng)的強(qiáng)度、頻率和持續(xù)時(shí)間對(duì)目標(biāo)損傷的影響,同時(shí)結(jié)合不同頻帶的時(shí)間歷程曲線可以獲得沖擊波和二次壓力波在不同頻帶上分布和衰減的細(xì)節(jié)信息,由此可確定沖擊波和二次壓力波對(duì)目標(biāo)毀傷的影響。韋灼彬等[39]采用數(shù)值模擬方法研究了港口水下爆炸荷載的沖擊特性。研究表明,港口水下爆炸氣泡膨脹荷載不可忽略,由比沖量可知,大部分區(qū)域氣泡膨脹荷載大于沖擊波荷載;自由水面對(duì)沖擊波荷載與氣泡膨脹荷載均有較大影響,越接近水面二者比沖量越小,氣泡膨脹荷載衰減越快,在水面附近其沖量甚至?xí)∮跊_擊波荷載。水底對(duì)氣泡膨脹荷載影響不大,但對(duì)沖擊波荷載影響較大;受水底反射波影響,沖擊波作用時(shí)間縮短,導(dǎo)致比沖量迅速減??;水底有淤泥層時(shí),沖擊波在泥層與水體交界面反射不顯著,水底反射沖擊波主要來(lái)自淤泥底部與巖石層分界面。

    C.F.Hung等[40]對(duì)圓筒殼結(jié)構(gòu)水下爆炸的動(dòng)態(tài)響應(yīng)進(jìn)行了研究,3種圓筒殼結(jié)構(gòu)分別是未強(qiáng)化、內(nèi)部強(qiáng)化和外部強(qiáng)化,圓筒殼結(jié)構(gòu)水中爆炸響應(yīng)試驗(yàn)在爆距35~210cm下進(jìn)行,沖擊波荷載由小藥量爆炸產(chǎn)生,爆炸距離小于50cm的試驗(yàn)?zāi)軌蛴^察到圓筒殼結(jié)構(gòu)的塑性變形,其他條件下的測(cè)試用于觀察圓筒的線性響應(yīng)。Mark Riley等[41]采用一種Eulerian計(jì)算流體動(dòng)態(tài)程序?qū)傂越Y(jié)構(gòu)附近水下爆炸荷載進(jìn)行了數(shù)值計(jì)算,將模擬結(jié)果與試驗(yàn)測(cè)試結(jié)果進(jìn)行了對(duì)比。研究表明,二維分析可以獲得良好的氣泡歷程以及相對(duì)于試驗(yàn)數(shù)據(jù)更容易理解的脈沖荷載;氣泡塌陷和射流出現(xiàn)的時(shí)間與網(wǎng)格關(guān)系顯著,同時(shí)網(wǎng)格的細(xì)化并不一定獲得好的計(jì)算結(jié)果,二維計(jì)算方法可以更好地理解計(jì)算中存在的問(wèn)題,三維計(jì)算可以獲得改善的脈沖預(yù)測(cè)結(jié)果,同時(shí)對(duì)數(shù)值計(jì)算獲得的氣泡擴(kuò)展時(shí)程與實(shí)際時(shí)程進(jìn)行了對(duì)比。

    Kenji Murata等[42]研究了水下沖擊波和氣泡脈沖荷載對(duì)模型鋼圓筒的作用,結(jié)果表明,第一次氣泡脈沖荷載對(duì)目標(biāo)的影響要高于沖擊波荷載,雖然沖擊波荷載峰值要高于第一次氣泡脈沖壓力,目標(biāo)變形程度與沖擊波脈沖和第一次氣泡脈沖呈線性增加關(guān)系。Erin Gauch等[43]通過(guò)數(shù)值模擬對(duì)水下爆炸載荷作用下預(yù)應(yīng)力薄復(fù)合材料板的響應(yīng)進(jìn)行了研究,考慮了直接壓應(yīng)力、非直接壓應(yīng)力、直接拉應(yīng)力3種平面預(yù)應(yīng)力的影響。通過(guò)板中心點(diǎn)曲度、材料破壞和剝落來(lái)評(píng)估預(yù)應(yīng)力條件的影響,結(jié)果表明,水下沖擊波荷載作用下,薄復(fù)合材料板受預(yù)應(yīng)力的影響最小,板的響應(yīng)或損傷和剝落量基本保持不變。

    Jian Li等[44]對(duì)水下爆炸作用下結(jié)構(gòu)的動(dòng)態(tài)響應(yīng)進(jìn)行了試驗(yàn)及數(shù)值模擬研究,結(jié)果表明,人工容積黏性系數(shù)對(duì)沖擊波峰值壓力影響顯著;沖擊波峰值壓力、氣泡脈動(dòng)周期以及圓柱的變形位移試驗(yàn)結(jié)果和有限元分析結(jié)果相似;討論了圓柱爆炸的長(zhǎng)徑比和沖擊波峰值壓力角的影響;測(cè)量了不同爆距下的塑性變形量,獲得了總的變形曲線。

    5水中爆炸測(cè)試試驗(yàn)技術(shù)研究

    水中爆炸測(cè)試試驗(yàn)技術(shù)可為毀傷機(jī)理研究提供科學(xué)有效的數(shù)據(jù)支持和研究手段,其中水中爆炸參數(shù)測(cè)量的傳感器技術(shù)影響至關(guān)重要。美國(guó)NOL實(shí)驗(yàn)室、英國(guó)NCRE實(shí)驗(yàn)室[45]及J.B.Dempesy[46]等在20世紀(jì)50年代至70年代對(duì)晶體的涂覆材料進(jìn)行了研究,發(fā)現(xiàn)蠟、硅膠、橡膠等會(huì)不同程度地影響傳感器的線性度及頻響,研究發(fā)現(xiàn)油封效果最好,可使傳感器的性能得到很大改善,也消除了蠟封傳感器得出的尺寸效應(yīng)和峰壓平臺(tái)效應(yīng)。

    ROUND.B.Tussing[47]對(duì)電氣石傳感器的涂覆材料作了系統(tǒng)的研究和闡述,提出采用裝滿硅油的聚乙烯塑料管進(jìn)行密封被證明是最穩(wěn)定、連續(xù)的可靠方法。美國(guó)PCB公司作為全球最大的傳感器生產(chǎn)商之一,生產(chǎn)的138系列水下爆炸沖擊波傳感器延續(xù)了這項(xiàng)研究成果。138系列傳感器具有較好的密封性,且具有體積敏感性,對(duì)安裝方向沒(méi)有要求,現(xiàn)在已經(jīng)被眾多國(guó)家和機(jī)構(gòu)廣泛用于水下爆炸壓力測(cè)試中。目前除了PCB公司138系列外,英國(guó)DRA生產(chǎn)的AWETR2-5000、美國(guó)NSWC生產(chǎn)的NOL系列及德國(guó)MOD生產(chǎn)的WTD-4000也是自由場(chǎng)壓力傳感器。

    GerryRude等[48]用AWETR2-5000G傳感器在深度為2m的小型爆炸容器中進(jìn)行了爆炸測(cè)試。雖然得到了完整的壓力時(shí)程曲線,但水面截?cái)嗪髩毫ω?fù)值達(dá)到了近-1MPa,認(rèn)為是傳感器自身原因引起的非正常現(xiàn)象,但沒(méi)有具體分析。John E.Slater等[49]對(duì)上述各種水下爆炸自由場(chǎng)壓力傳感器進(jìn)行了測(cè)試比較,測(cè)試的內(nèi)容包括熱幅射靈敏度、熱傳導(dǎo)靈敏度、靈敏度壓力標(biāo)定、水下爆炸沖擊波參數(shù)等;發(fā)現(xiàn)John E.SlaterPCB138系列、US/NOL系列和UK/AWETR2-5000G系列測(cè)試結(jié)果具有較好的重復(fù)性和可比性,德國(guó)的WTD-4000傳感器較其他傳感器測(cè)量峰值低10%~15%,并且出現(xiàn)共振現(xiàn)象,說(shuō)明碳阻傳感器在自由場(chǎng)壓力測(cè)試中的局限性。

    Kenji Murata和Daiki Hasegawa等[50-51]對(duì)PVDF壓電薄膜應(yīng)用于自由場(chǎng)壓力測(cè)量進(jìn)行了一系列研究。測(cè)得水下爆炸自由場(chǎng)壓力波形與使用電氣石測(cè)得的結(jié)果有較好的一致性。國(guó)內(nèi)利用PVDF壓電薄膜進(jìn)行近場(chǎng)壓力測(cè)量,并取得了一定的成果。池家春等[52]利用在高壓范圍內(nèi)能有效工作的錳銅壓阻壓力計(jì)和PVDF壓電壓力計(jì)兩種測(cè)試系統(tǒng)進(jìn)行了TNT/RDX炸藥球水中爆炸的探索性研究,在10倍藥包半徑的范圍內(nèi)測(cè)量了沖擊波峰值壓力的衰減規(guī)律。趙繼波等[53]用PVDF計(jì)測(cè)量了水中爆炸條件下0.7~3.0倍藥包半徑范圍內(nèi)的近場(chǎng)壓力,在沖擊波峰壓小于4GPa的情況下,得到了沖擊波的壓力衰減歷程和近場(chǎng)壓力剖面曲線。

    在高速攝影測(cè)試技術(shù)方面,Akio Kira[54]利用高速攝影技術(shù),獲得了大劑量球形爆炸物水中爆炸時(shí)的現(xiàn)象及水下沖擊波的傳播軌跡,并利用非線性曲線擬合技術(shù)和Arbitary Lagrangian-Eulerian方法進(jìn)行模擬實(shí)驗(yàn),獲得了很好的模擬結(jié)果。H.G.Snay[55]通過(guò)判讀水下爆炸高速攝影所獲得的氣泡產(chǎn)生和脈沖過(guò)程圖片,對(duì)半無(wú)限水介質(zhì)水下爆炸的流體力學(xué)進(jìn)行了深入研究;由于沖擊波的傳播速度與波頭壓力存在著蘭金-茍果尼歐條件關(guān)系,因此可以通過(guò)光學(xué)方法測(cè)量沖擊波傳播速度來(lái)獲得沖擊波壓力值。Bruceton[56]利用照相機(jī)進(jìn)行過(guò)壓力測(cè)量,測(cè)量結(jié)果和理論計(jì)算結(jié)果較接近。Menon S[57]、Lee J J[58]、Klaseboer E[59]和Huang C F[60]等人的研究表明:采用小當(dāng)量裝藥,在有限水域(如水箱或水池)中進(jìn)行水下爆炸氣泡的實(shí)驗(yàn)是可行的,而且,裝藥方式直接采用凝聚態(tài)炸藥比采用玻璃球充可燃?xì)怏w的效果要好,與實(shí)際情況更相符。

    在水中爆炸實(shí)驗(yàn)裝置研制方面,伍俊等設(shè)計(jì)和研制了具有多種功能的水中爆炸大型實(shí)驗(yàn)裝置[61-62],運(yùn)用一維彈性平面波理論,對(duì)該實(shí)驗(yàn)裝置結(jié)構(gòu)設(shè)計(jì)的主要問(wèn)題和抗爆隔震性能進(jìn)行了分析,并采用氣泡帷幕消波技術(shù),對(duì)降低水中爆炸沖擊波防護(hù)措施的可行性進(jìn)行了初步試驗(yàn)和探討;驗(yàn)證了該水中爆炸實(shí)驗(yàn)裝置結(jié)構(gòu)設(shè)計(jì)的合理性和降低水中沖擊波防護(hù)措施的有效性。

    6水中爆炸數(shù)值模擬與毀傷評(píng)估研究

    與試驗(yàn)研究相比,數(shù)值計(jì)算不僅省時(shí)省力,節(jié)約高昂的試驗(yàn)費(fèi)用,而且很容易避免試驗(yàn)研究中由于外界干擾及人為因素的影響而出現(xiàn)意想不到的誤差,對(duì)于推動(dòng)水中爆炸作用機(jī)理及毀傷效應(yīng)的研究都至關(guān)重要。

    A.Abe等[63]采用流體程序ANSYS AUTODYN模擬了水下爆炸及氣泡脈沖,結(jié)果表明,利用AUTODYN程序可以精確模擬氣泡脈沖現(xiàn)象,JWL-Miller模型可以成功給出非理想爆轟的相關(guān)現(xiàn)象。陳娟[64]以艦船近場(chǎng)水下爆炸為研究背景,對(duì)近場(chǎng)和近邊界水下爆炸問(wèn)題進(jìn)行了研究。近場(chǎng)和近邊界水下爆炸問(wèn)題涉及流場(chǎng)或結(jié)構(gòu)的極大變形、運(yùn)動(dòng)物質(zhì)交界面、沖擊波與邊界的相互作用和自由表面等特性,應(yīng)用傳統(tǒng)的基于網(wǎng)格的數(shù)值方法對(duì)此類問(wèn)題進(jìn)行模擬研究可能會(huì)由于網(wǎng)格發(fā)生嚴(yán)重的畸變而導(dǎo)致計(jì)算崩潰。陳衛(wèi)東等[65]認(rèn)為水下爆炸數(shù)值模擬中大變形、運(yùn)動(dòng)物質(zhì)交界面和多介質(zhì)耦合等問(wèn)題,給基于網(wǎng)格的數(shù)值算法帶來(lái)了很大的困難和挑戰(zhàn),并通過(guò)對(duì)TNT球形炸藥一維水下爆炸和TNT炸藥二維水下爆炸過(guò)程進(jìn)行了數(shù)值模擬,計(jì)算了沖擊波強(qiáng)度,分析了水下爆炸流場(chǎng)參量的空間分布和變化規(guī)律,并分別與AUTODYN、COLE經(jīng)驗(yàn)公式和SPH方法計(jì)算結(jié)果進(jìn)行了比較,結(jié)果吻合較好,為工程應(yīng)用提供了參考。

    楊文山等[66]對(duì)SPH方法模擬艦船水下爆炸的主要技術(shù)問(wèn)題及其研究現(xiàn)狀進(jìn)行了分析,結(jié)果表明,SPH方法作為一種新興的數(shù)值方法,在模擬大變形問(wèn)題時(shí)較網(wǎng)格方法有重大優(yōu)勢(shì),已成功應(yīng)用于水下爆炸的模擬。S.K.Chan[67]通過(guò)改進(jìn)有限元網(wǎng)絡(luò)模型,在計(jì)算無(wú)限水下爆炸的沖擊效應(yīng)時(shí),提高了計(jì)算精度,縮短了計(jì)算時(shí)間;張振華等[68]應(yīng)用商業(yè)有限元程序MSC.DYTRAN數(shù)值模擬了球形藥包在無(wú)限水域中爆炸產(chǎn)生的沖擊波,采用合理的計(jì)算參數(shù)和有限元模型,較好地模擬了水下爆炸沖擊波的傳播過(guò)程,通過(guò)調(diào)整水的狀態(tài)方程參數(shù),達(dá)到了提升沖擊波峰值應(yīng)力的效果,有效降低了單元數(shù)和計(jì)算時(shí)間。J.R.Britt等[69]運(yùn)用線性球體波動(dòng)理論編寫(xiě)了計(jì)算水下爆炸沖擊波在彈性水底反射變化的BOTREF程序;R.Willianson[70]通過(guò)修正CONGRATS程序和TRANSV程序,編寫(xiě)了計(jì)算水下爆炸沖擊波折射效應(yīng)和最大移動(dòng)速度的計(jì)算程序;對(duì)于半無(wú)限水介質(zhì)水下爆炸情況,主要存在自由水面對(duì)沖擊波的反射影響,Whitham方法使計(jì)算自由水面對(duì)水下爆炸沖擊波反射影響變得比較簡(jiǎn)捷。M.Kamegai等[71]采用CALE程序?qū)IGMAM水下爆炸實(shí)驗(yàn)和近自由面水下爆炸現(xiàn)象進(jìn)行數(shù)值模擬,計(jì)算結(jié)果均與實(shí)驗(yàn)結(jié)果一致;符松等[72]采用位標(biāo)函數(shù)計(jì)算法及NND計(jì)算格式對(duì)近水面水下爆炸進(jìn)行數(shù)值模擬,比跟蹤自由面的方法簡(jiǎn)潔實(shí)用;鄭哲敏[73]采用變分法計(jì)算淺水中爆炸時(shí)水的運(yùn)動(dòng)過(guò)程,計(jì)算結(jié)果與實(shí)際相符。田躍華等[74]根據(jù)實(shí)驗(yàn)結(jié)果,應(yīng)用最小二乘法、沖量等效、能量等效方式,得到兩種淺水爆炸有效破壞半徑的計(jì)算公式,柏勁松等[75]采用三階精度的PPM方法和VOF方法相結(jié)合,運(yùn)用Lagrange-Remapping算法,計(jì)算爆轟產(chǎn)物、水、空氣多種介質(zhì)之間的相互作用過(guò)程,建立了近水面水下爆炸3種介質(zhì)相互作用的二維數(shù)值計(jì)算模型,該方法可以用來(lái)處理界面兩邊高密度比可壓縮流動(dòng)以及強(qiáng)剪切滑移運(yùn)動(dòng)等問(wèn)題。

    陳菲等[76]在Fedkiw研究的基礎(chǔ)上采用二階迎風(fēng)TVD求解歐拉方程得到流場(chǎng)解,采用五階WENO求解lever set方程追蹤多流體界面,采用GFM方法處理流體內(nèi)界面,在此基礎(chǔ)上,對(duì)運(yùn)動(dòng)激波和氣泡串相互作用過(guò)程進(jìn)行了數(shù)值模擬,得到了不同時(shí)刻運(yùn)動(dòng)激波和圓管內(nèi)的兩個(gè)氣泡作用后的演化圖像。計(jì)算結(jié)果表明,運(yùn)動(dòng)激波經(jīng)過(guò)氣泡串后,激波強(qiáng)度會(huì)不斷衰減,說(shuō)明利用多個(gè)低密度氣泡會(huì)有效降低作用在目標(biāo)表面的載荷。

    吳國(guó)民等[77]研究了水中爆炸的數(shù)值模擬技術(shù),基于MSC.Dytran軟件,通過(guò)一維球?qū)ΨQ數(shù)值模型,模擬了水下爆炸沖擊波傳遞以及第一次氣泡脈動(dòng)過(guò)程,其計(jì)算結(jié)果與經(jīng)驗(yàn)公式計(jì)算結(jié)果吻合較好;在此基礎(chǔ)上討論了網(wǎng)格劃分方式、網(wǎng)格密度以及計(jì)算區(qū)域大小對(duì)計(jì)算結(jié)果精度的影響。周章濤等[78]應(yīng)用AUTODYN 軟件對(duì)沖擊波與水背圓柱殼、氣背圓柱殼的相互作用進(jìn)行了分析,在此基礎(chǔ)上對(duì)雙層圓柱殼內(nèi)外沖擊波場(chǎng)的演化過(guò)程進(jìn)行了數(shù)值模擬,分別在內(nèi)外殼表面附近水域設(shè)置了多個(gè)測(cè)點(diǎn)以記錄各點(diǎn)的壓力時(shí)程曲線,對(duì)沖擊波場(chǎng)的演化過(guò)程及衰減、匯聚現(xiàn)象進(jìn)行了研究。在水中爆炸毀傷效應(yīng)評(píng)估方面,侯曉峰等[79]為合理確定水中爆炸條件下作用在水中結(jié)構(gòu)上的荷載及其毀傷特性,通過(guò)分析比較現(xiàn)有理論及方法,提出了不同工況下的荷載計(jì)算方法、水中結(jié)構(gòu)毀傷等級(jí)劃分標(biāo)準(zhǔn)以及兩種毀傷評(píng)估工程方法。

    李萬(wàn)等[80]針對(duì)水下目標(biāo)抗水下爆炸實(shí)驗(yàn)數(shù)據(jù),利用小波包良好的時(shí)頻局部化性質(zhì),對(duì)被監(jiān)測(cè)水下目標(biāo)內(nèi)部裝置的沖擊加速度信號(hào)進(jìn)行了能量分析,得到了沖擊信號(hào)的時(shí)頻分布和不同頻帶上的能量分布。沖擊信號(hào)的頻帶能量分布與目標(biāo)毀傷的關(guān)系密切,選用沖擊信號(hào)峰值、沖擊信號(hào)主振頻帶能量、水下目標(biāo)內(nèi)部裝置自振頻率所在頻帶能量作為判別因子,建立距離判別模型,對(duì)水下目標(biāo)毀傷情況進(jìn)行了預(yù)測(cè)。陸路等[81]用經(jīng)驗(yàn)公式計(jì)算了大當(dāng)量炸藥水下爆炸時(shí)作用在壩體上的沖擊波壓強(qiáng),用ANSYS 軟件模擬了大壩作為彈性體在沖擊波作用下的應(yīng)力和應(yīng)變場(chǎng),從計(jì)算結(jié)果判斷壩體在所有高度將可能貫穿性斷裂,同時(shí)對(duì)斷裂后壩體運(yùn)動(dòng)作了定量估計(jì)。張倫平等[82]對(duì)某金屬材料制成的圓筒形試件進(jìn)行水下近距爆炸試驗(yàn),得到水下近距爆炸條件下材料的動(dòng)態(tài)斷裂應(yīng)變,結(jié)合能量方法,給出水下近距爆炸時(shí)材料的動(dòng)態(tài)屈服強(qiáng)度,并與靜態(tài)拉伸試驗(yàn)及霍普金森拉桿試驗(yàn)(SHTB)結(jié)果對(duì)比,比較屈服強(qiáng)度、極限拉伸應(yīng)變和能量吸收率等參數(shù),進(jìn)而探討了SHTB試驗(yàn)所得失效判據(jù)在水下爆炸條件下的適用性。R.Rajendran等[83]對(duì)水下爆炸作用下高強(qiáng)低合金鋼的性能進(jìn)行了評(píng)估,通過(guò)增加沖擊波強(qiáng)度,開(kāi)展了直徑290mm、厚度4mm圓板和300mm×250mm矩形板的水下爆炸隆起測(cè)試,以研究材料的塑性變形和斷裂誘因,并利用經(jīng)驗(yàn)?zāi)P蛠?lái)預(yù)測(cè)板隆起的高度;研究了具有不同裝藥和爆距的新型測(cè)試方法,證明了該方法與模型具有很好的一致性。F.Vanden Abeele等[84]對(duì)海底管道水下爆炸的損傷情況進(jìn)行了有限元分析,通過(guò)建立有限元模型來(lái)模擬水下爆炸條件下海底管道的瞬時(shí)響應(yīng),評(píng)估海底管道在水下爆炸作用下的行為和安全性,并將數(shù)值結(jié)果與試驗(yàn)結(jié)果進(jìn)行了對(duì)比。

    7結(jié)語(yǔ)及展望

    由以上研究現(xiàn)狀可知,水中爆炸經(jīng)過(guò)多年的研究發(fā)展,在水中爆炸作用機(jī)理、傳播規(guī)律、氣泡脈動(dòng)、破壞效應(yīng)、試驗(yàn)技術(shù)、數(shù)值模擬及毀傷評(píng)估等方面都有了很大進(jìn)展。但由于水中爆炸作用機(jī)理和毀傷效應(yīng)是一個(gè)復(fù)雜的物理過(guò)程,其爆炸沖擊波具有壓力峰值大、持續(xù)時(shí)間短等特點(diǎn),同時(shí)又涉及氣泡脈動(dòng)荷載及氣泡射流耦合等問(wèn)題,因而水中爆炸產(chǎn)生的破壞作用比在空氣中強(qiáng)烈得多,更容易造成水中結(jié)構(gòu)的破壞,因此,對(duì)水中爆炸作用機(jī)理及毀傷效應(yīng)研究還存在如下問(wèn)題迫切需要解決:

    (1)新型非理想炸藥水中爆炸能量釋放與傳播規(guī)律研究。

    目前,水中爆炸能量釋放與傳播規(guī)律的相關(guān)研究主要存在以下問(wèn)題:炸藥水中爆炸試驗(yàn)技術(shù)還比較單一,試驗(yàn)的有效性和準(zhǔn)確性尚需提高;對(duì)水中炸藥非理想爆轟的研究還不充分;對(duì)炸藥能量輸出結(jié)構(gòu)與炸藥組成及裝藥結(jié)構(gòu)關(guān)系不明確;對(duì)炸藥和水中兵器戰(zhàn)斗部在不同爆炸界面的爆炸效應(yīng)、水中爆炸沖擊波傳播和氣泡脈動(dòng)規(guī)律不清楚;對(duì)傳統(tǒng)炸藥在水中自由場(chǎng)狀態(tài)下的能量輸出特性研究比較充分,而對(duì)不同邊界條件下水中爆炸特性研究相對(duì)較少[85]。

    (2)沖擊波/氣泡耦合作用對(duì)結(jié)構(gòu)毀傷及動(dòng)態(tài)響應(yīng)研究。

    研究沖擊波/氣泡的相互作用對(duì)工程應(yīng)用具有重要意義。如水中爆炸,尤其是中近場(chǎng)爆炸所產(chǎn)生的沖擊波會(huì)在結(jié)構(gòu)物、海底以及自由面處發(fā)生反射,反射波在其傳播歷程中會(huì)遭遇爆炸氣泡,研究反射沖擊波與氣泡的相互耦合作用有利于分析爆炸產(chǎn)生的兩種主要載荷間的耦合作用,進(jìn)而深入揭示水中爆炸載荷對(duì)結(jié)構(gòu)的破壞作用[86]。由于這類耦合作用影響因素復(fù)雜、邊界條件不確定,目前對(duì)該方面深入的研究報(bào)道較少,爆炸現(xiàn)象難以觀察和測(cè)量,并且由于與目標(biāo)強(qiáng)烈的耦合效應(yīng),近場(chǎng)水中爆炸混凝土等非均質(zhì)材料結(jié)構(gòu)動(dòng)態(tài)響應(yīng)的試驗(yàn),以及精細(xì)化數(shù)值模擬研究也較少。因而加強(qiáng)近場(chǎng)水中爆炸混凝土結(jié)構(gòu)動(dòng)態(tài)響應(yīng)研究很有必要,可為水中爆炸毀傷機(jī)理研究提供有效的分析工具。

    (3)水中爆炸結(jié)構(gòu)毀傷評(píng)估方法研究。

    水中爆炸毀傷研究不僅在軍事上十分重要,而且對(duì)民用工業(yè)也很有意義,因而必須充分了解水中爆炸毀傷現(xiàn)象及機(jī)理。目前還沒(méi)有建立科學(xué)、統(tǒng)一的水中爆炸毀傷評(píng)估驗(yàn)證試驗(yàn)方法和評(píng)估標(biāo)準(zhǔn),長(zhǎng)期以來(lái)缺少炸藥水中爆炸、武器戰(zhàn)斗部水中爆炸威力和水中目標(biāo)毀傷效果的試驗(yàn)方法及其評(píng)估標(biāo)準(zhǔn),嚴(yán)重影響了水中爆炸毀傷評(píng)估研究。需要加強(qiáng)研究的關(guān)鍵問(wèn)題包括:水中爆炸毀傷作用及水中爆炸荷載的能量輸出規(guī)律研究;目標(biāo)對(duì)于炸藥水中爆炸近場(chǎng)能量輸出結(jié)構(gòu)荷載分布和動(dòng)態(tài)響應(yīng)過(guò)程的影響研究;水中爆炸對(duì)目標(biāo)毀傷判據(jù)的確定等。

    (4)水中爆炸多尺度數(shù)值模擬技術(shù)研究。

    水中爆炸試驗(yàn)研究消耗巨大,而爆炸試驗(yàn)也會(huì)對(duì)人員及周圍環(huán)境產(chǎn)生很大危險(xiǎn),考慮到水中爆炸的復(fù)雜性,水中爆炸試驗(yàn)研究的準(zhǔn)確性和效率也有待提高,開(kāi)發(fā)精確高效的水中爆炸研究方法成為水中爆炸研究面臨的問(wèn)題之一。從當(dāng)前我國(guó)水中爆炸數(shù)值模擬的研究現(xiàn)狀來(lái)看,主要存在以下幾個(gè)問(wèn)題:缺乏自行開(kāi)發(fā)的、成熟有效的大型水中爆炸數(shù)值模擬軟件以及多尺度混合單元高效算法,材料本構(gòu)模型及狀態(tài)方程的構(gòu)建還不完善,主要涉及樁基混凝土材料損傷模型、新型炸藥燃燒模型及能量輸出特性;水中爆炸邊界效應(yīng)的處理方法等還需深入研究。

    總之,水中爆炸作用機(jī)理及毀傷效應(yīng)研究涉及學(xué)科領(lǐng)域面廣,需要解決的關(guān)鍵技術(shù)問(wèn)題較多。長(zhǎng)期以來(lái),由于研究重點(diǎn)的不同,水中爆炸研究重點(diǎn)多以艦船安全防護(hù)為主,而針對(duì)水中工程結(jié)構(gòu)毀傷機(jī)理與效應(yīng)評(píng)估的問(wèn)題研究還相對(duì)較少。因此,進(jìn)一步加強(qiáng)水中爆炸作用機(jī)理及結(jié)構(gòu)破壞效應(yīng)研究,對(duì)于拓展研究思路,創(chuàng)新水中結(jié)構(gòu)防護(hù)措施和方法,提高水中結(jié)構(gòu)的保障能力及工程應(yīng)用水平具有重要意義。

    參考文獻(xiàn):

    [1]Pusey H C. Technical information support for survivability[J]. The Shock and Vibration Bulletin, 1983, 53(3):21-31.

    [2]Zong Z, Lan K Y. Dynamic plastic response of a submarine oil pipeline to an underwater explosion bubble[J].Journal of Applied Mechanics, 2000,67(4):758-762.

    [3]Zong Z. Dynamic plastic response of a submerged free-free beam to underwater gas bubble[J]. Acta Mechanica, 2003,161(3):179-194.

    [4]Zong Z, Lam K Y. The flexural response of a submarine pipeline to an underwater explosion bubble[J]. Journal of Offshore Mechanics and Engineering, 2000,122(3):194-199.

    [5]庫(kù)爾 P.水中爆炸[M].羅耀杰,等譯.北京:國(guó)防工業(yè)出版社,1960.

    Cole P. Underwater explosion [M]. LUO Yaojie, et al. Translate. Beijing: Publish of National Defense, 1960.

    [6]金輝,蔣國(guó)巖.某爆炸試驗(yàn)場(chǎng)論證報(bào)告[R].大連:中國(guó)人民解放軍91439部隊(duì),2008.

    JIN Hui, JIANG Guo-yan. Demonstration reports of certain explosion testing field [R]. Dalian: Army 91439 of PLA, 2008.

    [7]張慶明,金輝. 某爆炸國(guó)防重點(diǎn)實(shí)驗(yàn)室論證報(bào)告[R].北京:北京理工大學(xué),2008.

    ZHANG Qing-ming, JIN hui. Demonstration reports of certain major defense laboratory[R]. Beijing: Beijing Institute of Technology, 2008

    [8]張志江,徐更光.高能炸藥水中爆炸能量輸出特性數(shù)值分析[J].含能材料,2008,16(2):171-174.

    ZHANG Zhi-jiang, XU Geng-guang. Numerical analysis on energy output of underwater explosion for high energetic explosive [J]. Chinese Journal of Energetic Materials, 2008,16(2):171-174.

    [9]牛國(guó)濤,封雪松,曹少庭,等.含納米鋁炸藥水下爆炸試驗(yàn)研究[J].火工品,2013 (4):49-52.

    NIU Guo-tao, FENG Xue-song, CAO Shao-ting, et al. Underwater explosion test on aluminized explosive nano-aluminum [J]. Initiators and Pyrotechnics, 2013 (4):49-52.

    [10]史銳,徐更光,劉德潤(rùn),等.炸藥爆炸能量的水中測(cè)試與分析[J].火炸藥學(xué)報(bào),2008,31(4):1-5.

    SHI Rui, XU Geng-guang, LIU De-run, et al. Underwater test and analysis for explosion energy of explosives[J]. Chinese Journal of Explosives and Propellants, 2008,31(4):1-5.

    [11]馮曉軍,韋國(guó)平,徐洪濤,等.Al/AP對(duì)RDX基復(fù)合炸藥水中爆炸參數(shù)的影響[J].火炸藥學(xué)報(bào),2012,35(4):10-15.

    FENG Xiao-jun, WEI Guo-ping, XU Hong-tao, et al. Influence of Al/AP on the underwater explosion parameters of RDX-based composite explosives [J]. Chinese Journal of Explosives and Propellants, 2012,35(4):10-15.

    [12]宋浦,肖川,梁安定,等.不同起爆方式對(duì)TNT水中爆炸作用的影響[J].火炸藥學(xué)報(bào),2008,31(2):75-78.

    SONG Pu, XIAO Chuan, LIANG An-ding, et al. Effect of different initiating modes on TNT underwater explosion [J]. Chinese Journal of Explosives and Propellants, 2008,31(2):75-78.

    [13]周霖,徐少輝,徐更光.炸藥水下爆炸能量輸出特性研究[J].兵工學(xué)報(bào),2006,27(2):235-238.

    ZHOU Lin, XU Shao-hui, XU Geng-guang. Research on energy output characteristics for underwater explosion of explosives [J]. Acta Armamentarii, 2006,27(2):235-238.

    [14]閆偉杰,郭小溪.炸藥水中爆炸能量特性的試驗(yàn)及仿真[J].水面兵器,2010,19(1):43-47.

    YAN Wei-jie, GUO Xiao-xi. Experiment and simulation on characteristic of underwater explosion energy of explosive[J]. Surface Weapons, 2010,19(1):43-47.

    [15]趙根,季榮,鄭曉寧,等.乳化炸藥水中爆炸沖擊波傳播規(guī)律試驗(yàn)研究[J].爆破,2011,28(2):1-4.

    ZHAO Gen, JI Rong, ZHENG Xiao-ning, et al. Experimental investigation on propagation rule of shock wave by emulsion explosives underwater blasting[J]. Blasting, 2011,28(2):1-4.

    [16]李金河,趙繼波,譚多望,等. 炸藥水中爆炸的沖擊波性能[J].爆炸與沖擊,2009,29(2):172-176.

    LI Jin-he, ZHAO Ji-bo, TAN Duo-wang, et al. Underwater shock wave performances of explosives [J]. Explosion and Shock Waves, 2009,29(2):172-176.

    [17]劉志.水下爆炸沖擊波的傳播特性試驗(yàn)研究[D]. 重慶:西南交通大學(xué),2007.

    LIU Zhi. Research on transmission property of shock wave test underwater explosion [D]. Chongqing:Southwest Jiaotong University, 2007.

    [18]楊莉,汪玉,汪斌,等.沉底裝藥水中爆炸現(xiàn)象的實(shí)驗(yàn)研究[J].爆炸與沖擊,2013,33(2):175-180.

    YANG Li, WANG Yu, WANG Bin, et al. Experimental investigation on loading characteristics of underwater explosion from a bottom charge [J]. Explosion and Shock Waves, 2013,33(2):175-180.

    [19]項(xiàng)大林,榮吉利,李健,等.金屬殼體裝藥水下爆炸的沖擊波特性[J].爆炸與沖擊,2012,32(1):67-72.

    XIANG Da-lin, RONG Ji-li, LI Jian, et al. Shock wave features of underwater explosion of explosives with metal shell [J]. Explosion and Shock Waves, 2012,32(1):67-72.

    [20]盛振新,劉榮忠,郭銳.水下爆炸沖擊波相互作用的仿真分析[J].火工品,2012(3):25-29.

    SHENG Zhen-xin, LIU Rong-zhong, GUO Rui. Study on the shock wave interaction of underwater explosions[J]. Initiators and Pyrotechnics, 2012(3):25-29.

    [21]師華強(qiáng),宗智,賈敬蓓.水下爆炸沖擊波的近場(chǎng)特性[J].爆炸與沖擊,2009,29(2):125-130.

    SHI Hua-qiang, ZONG Zhi, JIA Jing-bei. Short-range characters of underwater blast waves[J]. Explosion and Shock Waves, 2009,29(2):125-130.

    [22]黃超,汪斌,姚熊亮,等.實(shí)驗(yàn)室尺度水下爆炸氣泡實(shí)驗(yàn)方法[J].傳感器與微系統(tǒng),2011,30(12):75-78.

    HUANG Chao, WANG Bin, YAO Xiong-liang, et al. Laboratory-scale underwater explosion bubble experiment method[J]. Transducer and Microsystem Technologies, 2011,30(12):75-78.

    [23]汪斌,張遠(yuǎn)平,王彥平.水中爆炸氣泡脈動(dòng)現(xiàn)象的實(shí)驗(yàn)研究[J].爆炸與沖擊,2008,28(6):272-276.

    WANG Bin, ZHANG Yuan-ping, WANG Yan-ping. Experimental study on bubble oscillation formed during underwater explosions [J]. Explosion and Shock Waves, 2008,28(6):272-276.

    [24]張阿漫,姚熊亮,聞雪友.自由場(chǎng)水中爆炸氣泡的物理特性[J].爆炸與沖擊,2008,28(5):391-400.

    ZHANG A-man, YAO Xiong-liang, WEN Xue-you. Physical behaviors of an underwater explosion bubble in a free field[J]. Explosion and Shock Waves, 2008,28(5):391-400.

    [25]Zhang A M ,Yao A L,Yu X B.The dynamics of three-dimensional underwater explosion bubble[J].Journal of Sound and Vibration, 2008,311:1196-1212.

    [26]黃超,汪斌,劉倉(cāng)理,等.柱形裝藥水中爆炸氣泡的形態(tài)演變[J].高壓物理學(xué)報(bào),2011,25(3):235-241.

    HUANG Chao, WANG Bin, LIU Cang-li, et al. On the shape evolution of underwater explosion bubble by cylindrical charges[J]. Chinese Journal of High Press Physics, 2011,25(3):235-241.

    [27]趙生偉,周剛,王占江,等.小當(dāng)量水中爆炸氣泡的脈動(dòng)現(xiàn)象[J].爆炸與沖擊,2009,29(2):213-216.

    ZHAO Sheng-wei, ZHOU Gang, WANG Zhan-jiang, et al. Bubble pulses of small-scale underwater explosion[J]. Explosion and Shock Waves, 2009,29(2):213-216.

    [28]金輝,周華,沈曉樂(lè),等.炸藥水下爆炸氣泡運(yùn)動(dòng)特性研究[C]∥第十一屆全國(guó)沖擊動(dòng)力學(xué)學(xué)術(shù)會(huì)議論文集. 北京:中國(guó)力學(xué)學(xué)會(huì),2011.

    JIN Hui, ZHOU Hua, SHEN Xiao-le, et al. Research on the dynamics of bubble induced by underwater explosions [C]∥The Eleventh National Shock Dynamics Academic Symposium. Beijing: Chinese Society of Theoretical and Applied Mechanics, 2011.

    [29]金輝,李兵,權(quán)琳,等.不同邊界條件下炸藥水中爆炸的能量輸出結(jié)構(gòu)[J].爆炸與沖擊,2013,33(3):325-329.

    JIN Hui, LI Bing, QUAN Lin, et al. Configuration of explosive energy output in different underwater boundary conditions[J]. Explosion and Shock Waves, 2013,33(3):325-329.

    [30]黃毅,張弩.水下爆炸氣泡脈動(dòng)荷載影響因素分析[J].中國(guó)艦船研究,2013,8(6):33-39.

    HUANG Yi, ZHANG Nu. Analysis of the impact parameters of the bubble pulsation load in underwater explosion[J]. Chinese Journal of Ship Research, 2013,8(6):33-39.

    [31]顧文彬,葉序雙,王自力.水中爆炸沖擊波作用下混凝土墩動(dòng)態(tài)響應(yīng)初步分析[J].爆炸與沖擊,1999,19(1):72-77.

    GU Wen-bin, YE Xu-shuang, WANG Zi-li. Dynamic response analysis of concrete body for underwater shock waves[J]. Explosion and Shock Wave, 1999,19(1):72-77.

    [32]賈憲振,胡毅亭,董明榮,等. 深水環(huán)境中水下爆炸沖擊波作用下圓柱殼動(dòng)態(tài)響應(yīng)的數(shù)值模擬研究[J].振動(dòng)與沖擊,2008(5):165-170.

    JIA Xian-zhen, HU Yi-ting, DONG Ming-rong, et al. Numerical simulation of dynamic response of cylindrical shell under the effect of explosion wave in deep water condition[J]. Journal of Vibration and Shock,2008(5):165-170.

    [33]徐俊祥,劉西拉.水中爆炸沖擊下混凝土壩動(dòng)力響應(yīng)的全耦合分析[J].上海交通大學(xué)學(xué)報(bào),2008,42(6):1001-1004.

    XU Jun-xiang, LIU Xi-la. Full coupled simulation of concrete dams subjected to underwater explosion[J].Journal of Shanghai JiaoTong University, 2008,42(6):1001-1004.

    [34]ZHU Feng,ZHU Wei-hua, ZHU Xiao-xi, et al. Numerical simulation of arch dam withstand underwater explosion[C]∥ 2012 International Conference on Modeling, Identification and Control.Wuhan:[s.n.],2012:24-26.

    [35]顧文彬,鄭向平,劉建青,等.淺層水中爆炸沖擊波對(duì)混凝土墩斜碰撞作用試驗(yàn)研究[J].爆炸與沖擊,2006,26(4):361-366.

    GU Wen-bin, ZHENG Xiang-ping, LIU Jian-qing, et al. Experimental investigation of the oblique collision effects of explosion shock wave on concrete frustum in shallow water[J]. Explosion and Shock Wave, 2006,26(4):361-366.

    [36]賴?guó)Q,馮順山,黃廣炎,等.不同加筋結(jié)構(gòu)在水中接觸爆炸下的破壞規(guī)律[J].爆炸與沖擊,2012,32(6):599-604.

    LAI Ming, FENG Shun-shan, HUANG Guang-yan, et al. Damage of different reinforced structure subjected to underwater contact explosion[J]. Explosion and Shock Wave, 2012,32(6):599-604.

    [37]李裕春,程克明,沈蔚,等.水中沖擊波對(duì)混凝土結(jié)構(gòu)破壞的實(shí)驗(yàn)研究[J].材料工程,2008(12): 19-23.

    LI Yu-chun, CHENG Ke-ming, SHEN Wei, et al. Damage analysis of concrete structure by underwater shock[J]. Material Engineer, 2008(12): 19-23.

    [38]李萬(wàn),張志華,周峰,等.水下目標(biāo)在水下爆炸作用下沖擊響應(yīng)的時(shí)頻特征[J].爆炸與沖擊,2012,32(3):309-315.

    LI Wan, ZHANG Zhi-hua, ZHOU Feng, et al. Time-frequency characteristics of shock responses of underwater target to underwater explosion[J]. Explosion and Shock Wave, 2012,32(3):309-315.

    [39]韋灼彬,唐廷,王立軍.港口水下爆炸荷載沖擊特性研究[J].振動(dòng)與沖擊,2014,33(6):18-25.

    WEI Zhuo-bin, TANG Ting, WANG Li-jun. Shock characteristics of underwater explosion in port[J]. Journal of Vibration and Shock, 2014,33(6):18-25.

    [40]Hung C F, Lin B J, Hwang-Fuu J J, et al.Dynamic response of cylindrical shell structures subjected to underwater explosion[J].Ocean Engineering, 2009,36:564-577.

    [41]Ma R, Ma S, van Aanhold J E, et al. Loading on a rigid target form close proximity underwater explosions[J].Shock and Vibration, 2012,19:555-571.

    [42]Murata K,Takahashi K,Kato Y. Underwater shock and bubble pulse loading against model steel cylinder[J]. Materials Science Forum, 2004,465:283-288.

    [43]Gauch E,Leblanc J,Shukla A. Response of thin composite panels subjected to underwater explosion loading[J]. Computers and Structures,2012(112/113): 342-353.

    [44]Li Jian,Rong Ji-li. Experimental and numerical investigation of the dynamic response of structure subjected to underwater explosion[J].European Journal of Mechanics B/Fluids, 2012,32:59-69.

    [45]Swift E, Slifko J J P. Recent work on measurement explosion pressures[J]. Ncltr, 1968,7:1-24.

    [46]Dempesy J B,Price R S.Reduction of scatter in underwater shock wave measurements made with piezoelectric gages, AD 739319[R]. Springfield: NTIS,1972.

    [47]Tussing R B. Accuracy and response of tourmaline gages for measurement of underwater explosion phenomena, ADA128734[R]. Springfield: NTIS, 1982: 1-41.

    [48]Rude G, Slate J E. Small-scale tank facility for studying underwater explosion phenomena[C]∥69th Shock and Vibration Symposium.Arvonia:HI-TEST Laboratories, 1998,10:663-673.

    [49]Slate J E, Rude G. Underwater pressure gauge comparison test study[C]∥ 69th Shock and Vibration Symposium.Arvonia:HI-TEST Laboratories, 1998,10:1-16.

    [50]Murata K, Takahashi K, Kato Y. Precise measurement of underwater explosion phenomena by pressure sensor using fluoropolymer [J]. Journal of Materials Processing Technology, 1999,85:39-42.

    [51]Hasegawa D, Inoue H, Murata K. Impulsive ultrasonic measurement for underwater shockwave detection[J]. Japanese Journal of Applied Physics, 2005,44(6B): 4694-4695.

    [52]池家春,馬冰. TNT/RDX(40/60)炸藥球水中爆炸波研究[J].高壓物理學(xué)報(bào),1999,13(3):199-204.

    CHI Jia-chun, MA Bing. Study on blast wave of underwater explosion of TNT/RDX(40/60)explosives[J]. Chinese Journal of High Pressure Physics, 1999,13(3):199-204.

    [53]趙繼波,譚多望,張遠(yuǎn)平,等.PVDF計(jì)在水中爆炸近場(chǎng)壓力測(cè)試中的應(yīng)用[J].火炸藥學(xué)報(bào),2009,32(3):1-4.

    ZHAO Ji-bo, TAN Duo-wang, ZHANG Yuan-ping, et al. Application of PVDF quantifier in the test of near field of underwater explosion[J]. Chinese Journal of Explosives and Propellants, 2009,32(3):1-4.

    [54]Akio Kira.Underwater explosion of spherical explosive[J]. Journal of Materials Processing Technology, 1999,85:64-68.

    [55]Snay H G. Hydrodynamics of underwater explosion[C]∥ Symposium Naval Hydrodyn1st. Washington D.C: [s.n.], 1965.

    [56]GASPIN J B. Depth scaling of underwater explosion phenomena, AD-A020473[R]. Springfield: NTIS, 1975.

    [57]Menon S, Lal M. On the dynamics and instability of bubbles formed during underwater explosions[J]. Experimental Thermal and Fluid Science,1998,16: 305-321.

    [58]Lee J J,Gregson J,Rude G, et al. Underwater shock and bubble interactions from twin explosive charges[C]∥26th International Symposium on Shock Waves. Man chester: Europeam Office of Aerospace Researchand Development,2009: 1549-1554.

    [59]Klaseboer E,Hung K C,Wang C, et al. Experimental and numerical investigation of the dynamics of an underwater explosion bubble near a resilient /rigid structure[J].Journal of Fluid Mech, 2005,537:387-413.

    [60]Huang C F,Hwangfu J J. Experimental study of the behaviors of mini-charge underwater explosion bubbles near different boundaries[J].Fluid Mech, 2010,605: 1-26.

    [61]伍俊,莊鐵栓,閆鵬,等.多功能水中爆炸實(shí)驗(yàn)裝置抗爆性能分析與試驗(yàn)研究[J].防護(hù)工程,2013,35(4):11-16.

    WU Jun, ZHUANG Tie-shuan, YAN Peng, et al. Development of a large-sized experiment device for water explosion[J]. Protective Engineering, 2013,35(4):11-16.

    [62]伍俊,莊鐵栓,閆鵬,等.水中爆炸實(shí)驗(yàn)裝置結(jié)構(gòu)設(shè)計(jì)與防護(hù)研究[J].振動(dòng)與沖擊,2013,32(11):131-136.

    WU Jun, ZHUANG Tie-shuan, YAN Peng, et al. Structure design and protection investigation of experiment device for water explosion[J]. Journal of Vibration and Shock, 2013,32(11):131-136.

    [63]Abe,Katayama,M. Murata,K. et al. Numerical study of underwater explosions and following bubble pulses[C]∥15th Topical Conference on Shock Compression of Condensed Matter.Hawaii: [s.n.], 2007.

    [64]陳娟.水下爆炸沖擊載荷的SPH算法研究[D].哈爾濱:哈爾濱工程大學(xué),2010.

    CHEN Juan. Study on the SPH algorithm of shock load of underwater explosion[D]. Harbin:Harbin Engineering University,2010.

    [65]陳衛(wèi)東,楊文淼,張帆.基于物質(zhì)點(diǎn)法的水下爆炸沖擊波數(shù)值模擬[J].高壓物理學(xué)報(bào),2013,27(6):813-820.

    CHEN Wei-dong, YANG Wen-miao, ZHANG Fan. Material point method for numerical simulation of underwater explosion blast wave[J]. Journal of High Press Physics, 2013,27(6):813-820.

    [66]楊文山,孟曉宇,王祖華.SPH方法模擬水下爆炸研究進(jìn)展[J].艦船科學(xué)技術(shù),2012,34(12):3-7.

    YANG Wen-shan ,MENG Xiao-yu, WANG Zu-hua. Review of applied researches in warship underwater explosion[J]. Ship Science and Technology, 2012,34(12):3-7.

    [67]Chan S K. An improvement in the modified finite element procedure for underwater shock analysis[C]∥Proceeding of 62nd Shock and Vibration Symposium.Arvonia:HI-TEST Laboratories,1992.

    [68]張振華,朱錫,白雪飛.水下爆炸沖擊波的數(shù)值模擬研究[J].爆炸與沖擊,2004,24(2):182-188.

    ZHANG Zhen-hua, ZHU Xi, BAI Xue-fei. Numerical simulation study on underwater shock waves, Explosion and Shock Wave[J]. 2004,24(2):182-188.

    [69]Britt J R. Bottom reflection of underwater explosion shock waves: computer program, AD-A204130[R].Springfield:NTIS,1971.

    [70]Williamson R. Prediction techniques for refraction of underwater explosion shock waves: programs for computing refraction effects and peak transnational velocity, AD-A003435[R]. Springfield:NTIS, 1974.

    [71]Kamegai M. A study of near surface and underwater explosions by computer simulations, UCRL-ID-116360[R].Springfield:NTIS,1994.

    [72]符松,王智平,張兆順,等.近水面水下爆炸的數(shù)值研究[J].力學(xué)學(xué)報(bào),1995,27(3):267-276.

    FU Song, WANG Zhi-ping, ZHANG Zhao-shun, et al. Numerical simulation on underwater explosion near free surface[J]. Chinese Journal of Mechanics, 1995,27(3):267-276.

    [73]鄭哲敏,楊振聲,金镠.爆炸處理水下海淤軟基[C]∥工程爆破文集.北京:冶金工業(yè)出版社,1993.

    ZHENG Zhe-min, YANG Zhen-sheng, JIN Liu. Cured silt-bottomed soft base by explosion[C]∥ Symposium of Engineering Blast. Beijing: Metallurgical Industry press, 1993.

    [74]田躍華,張志才.淺水爆炸對(duì)破壞半徑的影響[J].火炸藥學(xué)報(bào),2002(2):28-60.

    TIAN Yue-hua, ZHANG Zhi-cai. Influence of underwater explosion in shallow water to damage radius[J]. Chinese Journal of Explosives and Propellants, 2002(2):28-60.

    [75]柏勁松,陳森華,李平,等.水下爆炸過(guò)程的高精度數(shù)值計(jì)算[J].應(yīng)用力學(xué)學(xué)報(bào),2003,20(1):103-106.

    BAI Jin-song, CHEN Sen-hua, LI Ping, et al. High-accuracy numerical simulation of underwater explosion process[J]. Chinese Journal of Applied Mechanics, 2003,20(1):103-106.

    [76]陳菲,張夢(mèng)萍,徐勝利.運(yùn)動(dòng)激波和氣泡串相互作用的初步數(shù)值模擬[J].計(jì)算物理,2004,21(5):443-448.

    CHEN Fei, ZHANG Meng-ping, XU Sheng-li. Numerical studies of moving shock interacting with bubble string[J].Chinese Journal of Computation Physics, 2004,21(5):443-448.

    [77]吳國(guó)民,周心桃,肖漢林,等.水中爆炸數(shù)值仿真[J].艦船科學(xué)技術(shù),2012,34(9):20-26.

    WU Guo-min, ZHOU Xin-tao, XIAO Han-lin, et al. Numerical simulation of underwater explosion[J]. Ship Science and Technology, 2012,34(9):20-26.

    [78]周章濤,張倫平,何斌,等.水下爆炸沖擊波在圓柱殼內(nèi)外水域傳播的數(shù)值模擬[C]∥全國(guó)沖擊動(dòng)力學(xué)學(xué)術(shù)會(huì)議.北京:中國(guó)力學(xué)學(xué)會(huì),2006.

    ZHOU Zhang-tao, XU Lun-ping, HE Bin, et al. Numerical simulation of the near water area of cylindrical shell ship under the effect of underwater shock wave[C]∥ Symposium of Underwater Explosion. Beijing: Chinese Society of Theoretical and Applied Mechanics, 2006.

    [79]侯曉峰,王全勝,錢(qián)展芃,等.水中爆炸條件下結(jié)構(gòu)毀傷評(píng)估方法研究[J].防護(hù)工程,2014,36(1):40-45.

    HOU Xiao-feng, WANG Quan-sheng, QIAN Zhan-peng, et al. Study on the damage assessment methods of structure subjected to underwater explosion[J]. Protective Engineering, 2014,36(1):40-45.

    [80]李萬(wàn),張志華,陳滄海,等.水下爆炸毀傷水下目標(biāo)的能量分布特征[J].高壓物理學(xué)報(bào),2012,26(5):537-544.

    LI Wan, ZHANG Zhi-hua ,CHEN Cang-hai, et al. Features of energy distribution of underwater target by underwater explosion[J]. Chinese Journal of High Pressure Physics, 2012,26(5):537-544.

    [81]陸路,李昕,周晶.水下爆炸引起的混凝土壩危險(xiǎn)性估計(jì)[J].水資源與水工程學(xué)報(bào),2010,21(5):6-9.

    LU Lu, LI Xin, ZHOU Jing. Risk analysis of high concrete dam result from underwater explosion[J]. Journal of Water Resources and Water Engineering, 2010,21(5):6-9.

    [82]張倫平,嚴(yán)波,潘建強(qiáng),等.水下爆炸作用下金屬材料動(dòng)態(tài)失效判據(jù)[J].噪聲與振動(dòng)控制,2012(6): 54-57.

    ZHANG Lun-ping, YAN Bo, PAN Jian-qiang, et al. Dynamic invalidation criterion of metallic material during underwater explosion[J]. Nose and Shock Control, 2012(6): 54-57.

    [83]Rajendran R, Narasimhan K. Performance evaluation of HSLA steel subjected to underwater explosion[J].Journal of Materials Engineering and Performance, 2001,10(1):66-74.

    [84]Van den A F, Verleysen P. Finite element analysis of subsea pipelines subjected to underwater explosion[C]∥ Proceedings of the ASME 2013,32nd International Conference on Ocean. France: Offshore and Arctic Engineering,2013.

    [85]金輝,張慶明,高春生,等.不同邊界條件水下爆炸氣泡脈動(dòng)對(duì)比的實(shí)驗(yàn)研究[J].兵工學(xué)報(bào),2009,30(2):213-218.

    JIN Hui, ZHANG Qing-ming, GAO Chun-sheng, et al. Contrast test of explosion bubble pulsation in different underwater boundary condition[J]. Acta Armamentarii, 2009,30(2):213-218.

    [86]姚熊亮,汪玉,張阿漫.水下爆炸氣泡動(dòng)力學(xué)[M].哈爾濱:哈爾濱工程大學(xué)出版社,2011.

    YAO Xiong-liang, WANG Yu, ZHANG A-man. Dynamics of underwater explosion bubbles[M]. Harbin: Publish of Harbin Engineering University,2011.

    A Review of Research on Action Mechanism and Damage Effect of Underwater Explosion

    WU Jun1, YANG Yi1, ZHUANG Tie-shuan1,2

    (1.The Third Engineer Scientific Research Institute of the Headquarters of the General Staff, Luoyang Henan 471023, China; 2. State Key Laboratory of Disaster Prevention and Mitigation of Explosion and Impact,PLA University of Science and Technology, Nanjing 210007, China)

    Abstract:To further reveal the action mechanism and damage effect characteristics of underwater explosion, the progresses of study on the action mechanism and damage effect of underwater explosion were systematically reviewed from six aspects as mechanism, shock wave propagation, bubble pulsation, damage effect and dynamic response of structure, test technology and numerical simulation and damage assessment for underwater explosion of different explosives. Some key technical problems which should be investigated or solved in future including energy release and diffusion of the underwater explosion of non-ideal explosive, damage of coupling action of shock wave and air bubble to structure and damage assessment methods for underwater explosion, multi-scale simulation technology for underwater explosion were proposed with 86 references.

    Keywords:underwater explosion; action mechanism; damage effect; bubble pulsation; non-ideal explosive; demage assessment; review

    中圖分類號(hào):TJ55; TD235.2+1

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1007-7812(2016)01-0001-13

    作者簡(jiǎn)介:伍俊(1963-),男,高級(jí)工程師,從事防災(zāi)減災(zāi)及防護(hù)工程技術(shù)研究。E-mail: wujun61489@sina.com

    基金項(xiàng)目:國(guó)家自然科學(xué)基金資助(51578543)

    收稿日期:2015-11-09;修回日期:2015-11-22

    DOI:10.14077/j.issn.1007-7812.2016.01.001

    猜你喜歡
    作用機(jī)理綜述
    SEBS改性瀝青綜述
    石油瀝青(2018年6期)2018-12-29 12:07:04
    NBA新賽季綜述
    NBA特刊(2018年21期)2018-11-24 02:47:52
    近代顯示技術(shù)綜述
    電子制作(2018年14期)2018-08-21 01:38:34
    媒體關(guān)注影響企業(yè)債務(wù)融資的研究綜述與展望
    甲狀腺內(nèi)注射奧曲肽治療Graves病的臨床效果研究
    益生菌在畜禽生產(chǎn)中的作用及其研究進(jìn)展
    科技視界(2016年23期)2016-11-04 19:56:51
    生物強(qiáng)化技術(shù)在污水處理中的作用機(jī)理及應(yīng)用現(xiàn)狀
    義利之辨對(duì)企業(yè)社會(huì)責(zé)任影響的作用機(jī)理
    商(2016年23期)2016-07-23 11:02:24
    JOURNAL OF FUNCTIONAL POLYMERS
    農(nóng)業(yè)信息化帶動(dòng)農(nóng)業(yè)經(jīng)濟(jì)發(fā)展的作用機(jī)理研究
    日本免费一区二区三区高清不卡| 亚洲av免费高清在线观看| 91av网一区二区| 久久久午夜欧美精品| 国产毛片a区久久久久| 免费在线观看成人毛片| 国产精品蜜桃在线观看 | 哪里可以看免费的av片| 亚洲av免费在线观看| 黄片无遮挡物在线观看| 午夜a级毛片| 大香蕉久久网| av.在线天堂| 看黄色毛片网站| 国产在线男女| 秋霞在线观看毛片| 欧美成人精品欧美一级黄| 成人性生交大片免费视频hd| 精品无人区乱码1区二区| 一区二区三区四区激情视频 | 高清在线视频一区二区三区 | 五月玫瑰六月丁香| 天堂√8在线中文| 久久精品影院6| 国产一区二区三区在线臀色熟女| 国产高清激情床上av| 1000部很黄的大片| 人人妻人人澡欧美一区二区| 五月伊人婷婷丁香| 在线观看一区二区三区| 国产精品免费一区二区三区在线| 亚洲成人中文字幕在线播放| avwww免费| 久久久午夜欧美精品| 美女国产视频在线观看| 亚洲一级一片aⅴ在线观看| 高清日韩中文字幕在线| 色噜噜av男人的天堂激情| 国产国拍精品亚洲av在线观看| 黄片无遮挡物在线观看| 熟女电影av网| 成人毛片a级毛片在线播放| 日韩欧美 国产精品| 我的老师免费观看完整版| 成人毛片60女人毛片免费| 日本与韩国留学比较| 成人欧美大片| 国产探花在线观看一区二区| 亚洲av电影不卡..在线观看| 欧美一级a爱片免费观看看| 国产探花在线观看一区二区| 尾随美女入室| 日本熟妇午夜| 国产精品1区2区在线观看.| 国产精品1区2区在线观看.| 91午夜精品亚洲一区二区三区| 91午夜精品亚洲一区二区三区| 免费人成在线观看视频色| 日本撒尿小便嘘嘘汇集6| 男女那种视频在线观看| 国产综合懂色| 成人一区二区视频在线观看| 欧美激情国产日韩精品一区| 亚洲国产高清在线一区二区三| 成人毛片a级毛片在线播放| 欧美激情国产日韩精品一区| 亚洲经典国产精华液单| 国产av麻豆久久久久久久| 又黄又爽又刺激的免费视频.| 亚洲精品456在线播放app| 国产精品麻豆人妻色哟哟久久 | 熟女电影av网| 日本熟妇午夜| 亚洲第一区二区三区不卡| 亚洲成人久久性| 婷婷精品国产亚洲av| 99国产极品粉嫩在线观看| 色哟哟哟哟哟哟| 精品人妻偷拍中文字幕| av在线观看视频网站免费| 亚洲欧美成人综合另类久久久 | 黄色欧美视频在线观看| 午夜久久久久精精品| 中国美白少妇内射xxxbb| 高清午夜精品一区二区三区 | 国产真实乱freesex| 秋霞在线观看毛片| 成人漫画全彩无遮挡| 美女cb高潮喷水在线观看| 亚洲成人久久爱视频| 性欧美人与动物交配| 日韩视频在线欧美| 网址你懂的国产日韩在线| 91aial.com中文字幕在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲精品乱码久久久久久按摩| 色噜噜av男人的天堂激情| 亚洲欧洲日产国产| 国产成人影院久久av| 国产单亲对白刺激| 99久久成人亚洲精品观看| 内射极品少妇av片p| 国产爱豆传媒在线观看| 亚洲美女视频黄频| 超碰av人人做人人爽久久| 国产精品嫩草影院av在线观看| 桃色一区二区三区在线观看| 真实男女啪啪啪动态图| 久久久久国产网址| 日本免费一区二区三区高清不卡| 国产不卡一卡二| 激情 狠狠 欧美| 最后的刺客免费高清国语| 婷婷色av中文字幕| 看片在线看免费视频| 网址你懂的国产日韩在线| 精华霜和精华液先用哪个| 99久久中文字幕三级久久日本| 久久国产乱子免费精品| 美女大奶头视频| 97在线视频观看| 日本一二三区视频观看| 欧美成人精品欧美一级黄| 亚洲欧美成人精品一区二区| 亚洲电影在线观看av| 26uuu在线亚洲综合色| 国产精品爽爽va在线观看网站| 国产成人影院久久av| 如何舔出高潮| 18+在线观看网站| 精品一区二区三区人妻视频| 久久婷婷人人爽人人干人人爱| 激情 狠狠 欧美| 熟女人妻精品中文字幕| 亚洲国产精品久久男人天堂| 99久久精品一区二区三区| 国产精品久久久久久久电影| 深夜a级毛片| 我要看日韩黄色一级片| 亚洲最大成人手机在线| 国产单亲对白刺激| 国产精品人妻久久久久久| 国产高清不卡午夜福利| 婷婷色综合大香蕉| 日韩强制内射视频| 精品一区二区三区视频在线| 午夜精品一区二区三区免费看| 中国美白少妇内射xxxbb| 在线观看av片永久免费下载| 在线免费观看的www视频| 国产精品久久久久久av不卡| 午夜激情欧美在线| 青春草视频在线免费观看| 国语自产精品视频在线第100页| 亚洲精品久久久久久婷婷小说 | 欧美+亚洲+日韩+国产| 亚洲国产欧美在线一区| 岛国在线免费视频观看| 少妇人妻精品综合一区二区 | 欧美区成人在线视频| 久久婷婷人人爽人人干人人爱| 97在线视频观看| 成人性生交大片免费视频hd| 国产乱人偷精品视频| 国产精品一区二区三区四区免费观看| 91午夜精品亚洲一区二区三区| 国产成人91sexporn| 欧美激情在线99| 老司机福利观看| 小蜜桃在线观看免费完整版高清| 国产精品乱码一区二三区的特点| 成年女人看的毛片在线观看| 国产精品av视频在线免费观看| av在线天堂中文字幕| 乱系列少妇在线播放| 精品久久久噜噜| 一级黄色大片毛片| 久久精品夜夜夜夜夜久久蜜豆| 亚洲无线在线观看| 成人av在线播放网站| 国产色婷婷99| 国产熟女欧美一区二区| 少妇被粗大猛烈的视频| 少妇裸体淫交视频免费看高清| 亚洲国产精品成人久久小说 | 久久婷婷人人爽人人干人人爱| 高清日韩中文字幕在线| 日日干狠狠操夜夜爽| 成人鲁丝片一二三区免费| 尾随美女入室| 女同久久另类99精品国产91| av免费在线看不卡| 能在线免费看毛片的网站| 99久久人妻综合| 波野结衣二区三区在线| 最新中文字幕久久久久| 丝袜美腿在线中文| 麻豆精品久久久久久蜜桃| 真实男女啪啪啪动态图| 国产精品久久久久久久电影| 亚洲av.av天堂| 国产成人影院久久av| 国产成人精品一,二区 | 一区二区三区四区激情视频 | 午夜福利在线在线| 久久精品影院6| 日本欧美国产在线视频| 亚洲三级黄色毛片| 国产午夜福利久久久久久| 午夜福利在线在线| 久久精品影院6| 亚洲欧美清纯卡通| 亚洲人成网站在线播| 网址你懂的国产日韩在线| 菩萨蛮人人尽说江南好唐韦庄 | 中文精品一卡2卡3卡4更新| 欧美色视频一区免费| 天堂av国产一区二区熟女人妻| 国产精品久久久久久久电影| 又粗又硬又长又爽又黄的视频 | 大又大粗又爽又黄少妇毛片口| 久久精品91蜜桃| 爱豆传媒免费全集在线观看| 亚洲av免费高清在线观看| av专区在线播放| 中文字幕av成人在线电影| 国产精品爽爽va在线观看网站| 99在线视频只有这里精品首页| 男女那种视频在线观看| 永久网站在线| 亚洲熟妇中文字幕五十中出| 99久久精品国产国产毛片| 亚洲在线观看片| 亚洲一区二区三区色噜噜| 精品久久国产蜜桃| 国产美女午夜福利| 国产在线男女| 寂寞人妻少妇视频99o| 天堂√8在线中文| 亚洲人与动物交配视频| 久久亚洲国产成人精品v| 日韩欧美精品v在线| 99在线人妻在线中文字幕| 国产老妇伦熟女老妇高清| 免费一级毛片在线播放高清视频| 黄色配什么色好看| 97人妻精品一区二区三区麻豆| 丰满乱子伦码专区| 一级黄片播放器| 国产精品不卡视频一区二区| 欧美日韩在线观看h| 亚洲aⅴ乱码一区二区在线播放| 国产视频内射| 内射极品少妇av片p| 亚洲美女视频黄频| 一本久久中文字幕| 精品日产1卡2卡| 久久鲁丝午夜福利片| 日韩中字成人| 嫩草影院入口| 老师上课跳d突然被开到最大视频| 婷婷精品国产亚洲av| 爱豆传媒免费全集在线观看| 一本一本综合久久| 国产激情偷乱视频一区二区| 自拍偷自拍亚洲精品老妇| 天美传媒精品一区二区| 亚洲美女搞黄在线观看| 日本-黄色视频高清免费观看| 丰满的人妻完整版| 国产精品嫩草影院av在线观看| 免费观看精品视频网站| 国产三级中文精品| 国产精品女同一区二区软件| 成年av动漫网址| 夜夜看夜夜爽夜夜摸| 亚洲精品粉嫩美女一区| 久久精品久久久久久久性| 看十八女毛片水多多多| 国产高清激情床上av| 99久国产av精品国产电影| 天堂av国产一区二区熟女人妻| 91久久精品电影网| 国产一区二区激情短视频| 日本黄色片子视频| 三级男女做爰猛烈吃奶摸视频| 久久中文看片网| 美女 人体艺术 gogo| 日本熟妇午夜| 少妇被粗大猛烈的视频| 久久九九热精品免费| 欧美高清成人免费视频www| 亚洲欧美精品专区久久| 99久久无色码亚洲精品果冻| 精品久久久噜噜| 日日啪夜夜撸| 99国产精品一区二区蜜桃av| 最近最新中文字幕大全电影3| 不卡一级毛片| 综合色av麻豆| 特大巨黑吊av在线直播| 久久精品91蜜桃| 欧美性猛交黑人性爽| 人妻久久中文字幕网| 午夜福利在线观看免费完整高清在 | 菩萨蛮人人尽说江南好唐韦庄 | 超碰av人人做人人爽久久| 国产久久久一区二区三区| 小说图片视频综合网站| 亚洲精品色激情综合| 人妻系列 视频| 久久亚洲国产成人精品v| 26uuu在线亚洲综合色| 在现免费观看毛片| 亚洲精华国产精华液的使用体验 | 99热6这里只有精品| 免费看光身美女| 亚洲一级一片aⅴ在线观看| 日韩 亚洲 欧美在线| 看免费成人av毛片| 国产高清视频在线观看网站| 亚洲成a人片在线一区二区| 网址你懂的国产日韩在线| 麻豆国产97在线/欧美| 免费看av在线观看网站| 在线观看免费视频日本深夜| 精品午夜福利在线看| a级毛色黄片| 在线观看66精品国产| 久久久久久久午夜电影| 99久国产av精品| 变态另类成人亚洲欧美熟女| 又粗又硬又长又爽又黄的视频 | 国产精品乱码一区二三区的特点| 菩萨蛮人人尽说江南好唐韦庄 | 天天躁日日操中文字幕| 国产av麻豆久久久久久久| 精品久久久久久久久亚洲| 麻豆国产av国片精品| 不卡视频在线观看欧美| 国产私拍福利视频在线观看| 秋霞在线观看毛片| 一级黄色大片毛片| 日韩强制内射视频| 久久久久久久久久久丰满| 亚洲精品乱码久久久久久按摩| 欧美激情在线99| 成人一区二区视频在线观看| 在线观看一区二区三区| 久99久视频精品免费| 国产伦在线观看视频一区| 国产精品人妻久久久影院| 麻豆乱淫一区二区| 精品不卡国产一区二区三区| 我要看日韩黄色一级片| 国产综合懂色| 中文欧美无线码| 一区二区三区免费毛片| 免费av观看视频| 免费av毛片视频| 亚洲成a人片在线一区二区| 亚洲电影在线观看av| 青青草视频在线视频观看| 男女那种视频在线观看| 哪个播放器可以免费观看大片| 久久久久九九精品影院| 尾随美女入室| 国产精品一区二区三区四区久久| 99热全是精品| 你懂的网址亚洲精品在线观看 | 午夜福利在线观看吧| 尾随美女入室| 我的女老师完整版在线观看| 美女高潮的动态| 国产精品久久久久久久电影| 午夜亚洲福利在线播放| 国产色婷婷99| 青春草国产在线视频 | 久久久成人免费电影| 亚洲成a人片在线一区二区| 国产成人精品久久久久久| а√天堂www在线а√下载| 啦啦啦韩国在线观看视频| 久久这里只有精品中国| 久久人人爽人人片av| 日本-黄色视频高清免费观看| 免费大片18禁| 99热精品在线国产| 久久精品国产亚洲av香蕉五月| 少妇人妻精品综合一区二区 | 国产伦精品一区二区三区视频9| av女优亚洲男人天堂| 黄色日韩在线| 欧美一区二区亚洲| 99精品在免费线老司机午夜| 97人妻精品一区二区三区麻豆| 一本久久中文字幕| 免费看光身美女| 亚洲无线在线观看| 99热这里只有是精品50| 久久精品国产自在天天线| 能在线免费观看的黄片| 一个人看的www免费观看视频| 色综合亚洲欧美另类图片| 国模一区二区三区四区视频| 内射极品少妇av片p| 久久久成人免费电影| 国产亚洲av嫩草精品影院| 亚洲国产精品合色在线| 亚洲最大成人av| 亚洲av成人精品一区久久| 男女做爰动态图高潮gif福利片| 久久国内精品自在自线图片| 伦精品一区二区三区| 亚洲精品粉嫩美女一区| 国产精品麻豆人妻色哟哟久久 | 久久久久久久久久成人| 国产成人福利小说| 晚上一个人看的免费电影| 亚洲精品亚洲一区二区| 日本一二三区视频观看| 99久久九九国产精品国产免费| 国产精品日韩av在线免费观看| 九九爱精品视频在线观看| 自拍偷自拍亚洲精品老妇| 三级毛片av免费| 国产高潮美女av| 春色校园在线视频观看| .国产精品久久| 亚洲欧美日韩卡通动漫| 久久久久免费精品人妻一区二区| 国产精品国产三级国产av玫瑰| 晚上一个人看的免费电影| 久久99精品国语久久久| 日韩亚洲欧美综合| 一级黄色大片毛片| 国产精品久久久久久av不卡| 日本熟妇午夜| 我要看日韩黄色一级片| 小蜜桃在线观看免费完整版高清| 热99在线观看视频| 精品久久久久久久久久久久久| 女同久久另类99精品国产91| 99在线人妻在线中文字幕| 亚洲图色成人| 国产不卡一卡二| 91精品国产九色| 91久久精品电影网| 一本一本综合久久| 日本-黄色视频高清免费观看| 搡女人真爽免费视频火全软件| 久久久久性生活片| 一区二区三区高清视频在线| 亚洲精品国产av成人精品| 波野结衣二区三区在线| 69av精品久久久久久| 色5月婷婷丁香| 狂野欧美激情性xxxx在线观看| 日韩欧美国产在线观看| avwww免费| 欧美极品一区二区三区四区| 2022亚洲国产成人精品| 可以在线观看的亚洲视频| 色哟哟哟哟哟哟| 久久九九热精品免费| 国产老妇女一区| 午夜福利在线观看吧| 成年版毛片免费区| 国产高清有码在线观看视频| 欧美+日韩+精品| av在线天堂中文字幕| 久久久欧美国产精品| 97在线视频观看| 精品久久久久久久人妻蜜臀av| 亚洲美女视频黄频| 亚洲不卡免费看| 村上凉子中文字幕在线| 在现免费观看毛片| 在线天堂最新版资源| 最后的刺客免费高清国语| 一级黄色大片毛片| 91久久精品国产一区二区成人| 国产午夜精品论理片| 偷拍熟女少妇极品色| 国产伦理片在线播放av一区 | 亚洲一区高清亚洲精品| 男人舔奶头视频| 中文字幕精品亚洲无线码一区| 天堂中文最新版在线下载 | av女优亚洲男人天堂| av在线蜜桃| 欧美色视频一区免费| 悠悠久久av| 欧美色视频一区免费| 欧美精品一区二区大全| 午夜老司机福利剧场| 黄色日韩在线| 一级毛片我不卡| 你懂的网址亚洲精品在线观看 | 此物有八面人人有两片| 日本欧美国产在线视频| 一个人看视频在线观看www免费| 97超视频在线观看视频| 中文在线观看免费www的网站| 国产精品av视频在线免费观看| 午夜福利成人在线免费观看| 26uuu在线亚洲综合色| 久久九九热精品免费| 国产免费一级a男人的天堂| av女优亚洲男人天堂| 人妻夜夜爽99麻豆av| 国产亚洲av嫩草精品影院| av.在线天堂| 久久九九热精品免费| 国产高潮美女av| 99热6这里只有精品| 日日啪夜夜撸| 蜜桃亚洲精品一区二区三区| 热99在线观看视频| 中文在线观看免费www的网站| 日韩精品有码人妻一区| 久久精品国产亚洲av天美| 国产精品久久久久久精品电影| 国产v大片淫在线免费观看| 亚洲成人中文字幕在线播放| 免费看av在线观看网站| 久久久国产成人免费| 亚洲最大成人中文| 色播亚洲综合网| 人体艺术视频欧美日本| 国内精品宾馆在线| 男的添女的下面高潮视频| 国产中年淑女户外野战色| 亚洲欧美精品专区久久| 男人狂女人下面高潮的视频| 成人国产麻豆网| 又爽又黄a免费视频| 亚洲国产欧美在线一区| 看片在线看免费视频| 18禁在线播放成人免费| 久久久久国产网址| 卡戴珊不雅视频在线播放| 变态另类成人亚洲欧美熟女| 亚洲四区av| 一进一出抽搐gif免费好疼| 日韩一区二区三区影片| 国产一区二区在线观看日韩| 别揉我奶头 嗯啊视频| 99久久中文字幕三级久久日本| 久久久a久久爽久久v久久| av免费在线看不卡| 可以在线观看毛片的网站| or卡值多少钱| 一个人看视频在线观看www免费| 日日干狠狠操夜夜爽| 婷婷精品国产亚洲av| 国内精品宾馆在线| 亚洲激情五月婷婷啪啪| 免费看日本二区| 久久精品国产亚洲av香蕉五月| 国产亚洲av片在线观看秒播厂 | 九九爱精品视频在线观看| 色综合亚洲欧美另类图片| 日韩视频在线欧美| 嫩草影院精品99| 日韩,欧美,国产一区二区三区 | 欧美激情在线99| 夜夜看夜夜爽夜夜摸| 国产人妻一区二区三区在| 只有这里有精品99| 99热这里只有精品一区| 日韩三级伦理在线观看| 看片在线看免费视频| 69人妻影院| 女同久久另类99精品国产91| 精品人妻熟女av久视频| 免费看a级黄色片| 精品久久久久久成人av| 日韩三级伦理在线观看| 国产精品人妻久久久久久| 99国产精品一区二区蜜桃av| 亚洲第一区二区三区不卡| 少妇熟女欧美另类| 一区福利在线观看| 成人综合一区亚洲| 欧美一区二区亚洲| 成人高潮视频无遮挡免费网站| 91久久精品国产一区二区成人| 婷婷精品国产亚洲av| 欧美色视频一区免费| 亚洲自拍偷在线| 中文字幕人妻熟人妻熟丝袜美| 99热只有精品国产| 久久久欧美国产精品| 毛片一级片免费看久久久久| 免费观看在线日韩| 精品无人区乱码1区二区| 色噜噜av男人的天堂激情| 中文字幕久久专区| 一级毛片久久久久久久久女| 成人亚洲欧美一区二区av| 能在线免费观看的黄片| 婷婷色综合大香蕉| 精品久久久久久久久久免费视频| 久久久久久大精品| 老熟妇乱子伦视频在线观看| 国产片特级美女逼逼视频| 久久国产乱子免费精品| 久久久久久久久久久丰满| 成年女人永久免费观看视频| 欧美激情久久久久久爽电影| 少妇裸体淫交视频免费看高清| 欧美激情久久久久久爽电影| 午夜福利视频1000在线观看| 久久精品人妻少妇|