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

    基于LS-DYNA 的液電效應(yīng)沖擊波數(shù)值模擬*

    2022-03-17 07:25:08楊睿智
    爆炸與沖擊 2022年2期
    關(guān)鍵詞:沖擊波炸藥等離子體

    余 慶,張 輝,楊睿智

    (中國石油大學(xué)(北京)石油工程學(xué)院,北京 102249)

    自蘇聯(lián)Yutkin 把液電效應(yīng)作為沖擊動(dòng)力源并廣泛應(yīng)用于工業(yè)加工以來,液電效應(yīng)研究進(jìn)入迅速發(fā)展階段,學(xué)者們開展了大量的研究工作,強(qiáng)大的沖擊波已經(jīng)被用于體外碎石、液電成形、油田解堵、壓裂等領(lǐng)域。截至目前,液電效應(yīng)的產(chǎn)生機(jī)理仍是人們研究的重點(diǎn),液相放電擊穿理論主要分為場(chǎng)致電離理論和氣泡理論,其中“氣泡擊穿理論”更為普遍接受?;谠摾碚?,可以將液相放電等離子體的形成過程概括為:電極加熱→氣層(泡)形成→氣層擊穿→形成等離子體通道,推動(dòng)周圍液體介質(zhì)形成沖擊波。實(shí)際應(yīng)用中,通常把氣層擊穿之前的過程稱為預(yù)擊穿過程,之后的過程稱為主放電過程。由于在液體擊穿過程中伴隨著強(qiáng)烈的熱、電、機(jī)械等物理化學(xué)變化,因而對(duì)于液電效應(yīng)擊穿機(jī)理的研究并不充分,至今仍缺乏一個(gè)被研究者普遍接受的擊穿機(jī)理解釋理論。而從實(shí)驗(yàn)角度出發(fā)研究液相等離子體的特性和機(jī)理,通常受到測(cè)量?jī)x器的限制,不能準(zhǔn)確地獲取等離子體通道內(nèi)部的特性,并且受環(huán)境因素影響嚴(yán)重,這使得測(cè)量數(shù)據(jù)的準(zhǔn)確性進(jìn)一步下降。

    綜上所述,由于實(shí)驗(yàn)條件的限制以及擊穿理論的不成熟,至今鮮有成熟的商業(yè)數(shù)值模擬軟件能夠有效模擬液電效應(yīng)的整個(gè)過程,一般模擬方式是通過假設(shè)等離子體通道阻抗模型,并結(jié)合等離子體通道能量守恒方程、LCR 等效電路方程、沖擊波近似方程進(jìn)行一維的數(shù)值模擬計(jì)算。然而,由于通道電阻隨通道內(nèi)溫度、數(shù)密度等因素的變化而變化,準(zhǔn)確估計(jì)等離子體通道的電阻變化十分困難。工程上通常利用數(shù)值模擬軟件對(duì)一些具體工況進(jìn)行模擬,分析液電效應(yīng)產(chǎn)生的沖擊波對(duì)結(jié)構(gòu)物的作用效果,而僅僅是一維的模擬計(jì)算數(shù)據(jù)尚不足以滿足工程需要,因此,開展將液電效應(yīng)產(chǎn)生的沖擊波運(yùn)用于已有數(shù)值模擬軟件中的研究具有重要意義。

    本文中,介紹兩種基于商用顯示動(dòng)力學(xué)軟件LS-DYNA 間接模擬液電效應(yīng)產(chǎn)生沖擊波的方法:水下爆炸等效方法和理想氣體等效方法,并對(duì)計(jì)算結(jié)果進(jìn)行比較,給出改進(jìn)建議,以期為液電效應(yīng)產(chǎn)生沖擊波的工程應(yīng)用提供參考依據(jù)。

    1 液電效應(yīng)沖擊波特性分析

    在介紹液電效應(yīng)沖擊波等效方法之前,有必要對(duì)液電效應(yīng)沖擊波的形成與能量特性進(jìn)行分析。圖1為液電效應(yīng)產(chǎn)生沖擊波的示意圖??梢钥闯?,液電效應(yīng)產(chǎn)生沖擊波的過程是一個(gè)瞬時(shí)、連續(xù)且復(fù)雜的能量轉(zhuǎn)換過程。首先向電容器充電,當(dāng)電容充電電壓達(dá)到要求后,觸發(fā)控制開關(guān),這樣存儲(chǔ)在電容器內(nèi)的電能迅速注入到位于液電介質(zhì)中的正負(fù)極兩端,位于正負(fù)電極之間的水介質(zhì)逐漸被加熱、氣化,形成等離子體通道并在通道與液體介質(zhì)界面處向外輻射沖擊波,并且過程中還伴隨著光、熱輻射和熱傳導(dǎo)。在間隙擊穿之后,放電回路可以視為RLC 等效放電回路,圖2 為液相放電等效電路,圖中為儲(chǔ)能電容,和分別為回路電阻和等離子體通道電阻,和分別為回路電感和等離子體通道電感,=+為電路總電感。這里將看作常數(shù),并忽略的影響,基于基爾霍夫定律就可列出微分方程:

    圖1 沖擊波形成示意圖Fig. 1 Schematic diagram of shock wave generation

    圖2 液相放電等效電路Fig. 2 Equivalent circuit of high-voltage discharge in a liquid

    在形成沖擊波的整個(gè)過程中主要涉及4 個(gè)能量概念:(1)儲(chǔ)存在電容器內(nèi)的初始能量;(2)液體介質(zhì)擊穿時(shí)儲(chǔ)存在電容器內(nèi)的能量;(3)注入到等離子體通道內(nèi)的沉積能量;(4)由沉積能量轉(zhuǎn)換為等離子體通道內(nèi)的機(jī)械能,即沖擊波能量。這些能量之間的關(guān)系14]可以表示為:

    式中:η為充電電容儲(chǔ)能轉(zhuǎn)換為擊穿時(shí)電容電能的效率,η為放電能量轉(zhuǎn)換為等離子體通道沉積能量的效率,η為等離子通道沉積能量轉(zhuǎn)換為沖擊波能量的效率,η 為總能量轉(zhuǎn)換效率,為初始電壓,為擊穿電壓,()為等離子體通道電功率,τ 為主放電過程持續(xù)時(shí)間,()為放電間隙兩端測(cè)量的電流,為壓力傳感器到放電通道的水平距離,ρ為液體密度,為液體介質(zhì)中沖擊波的傳播速度,為測(cè)得的沖擊壓力。

    為了能夠?qū)ι鲜龈拍钣星宄恼J(rèn)識(shí),這里使用文獻(xiàn)[15-16]中的典型實(shí)驗(yàn)波形,如圖3 所示,圖3(a)描述了高壓探頭距離放電間隙中心17 cm 處等離子體通道的電流、電壓、沖擊波變化,圖3(b)描述了以擊穿時(shí)刻作為初始時(shí)刻,由式(3) 計(jì)算得到的沉積功率與沉積能量變化。通過圖像數(shù)據(jù),利用式(1)~(6)分別計(jì)算,可得到η=73.38%,η=47.6%,η=28.2%,η=9.85%。由于不同的實(shí)驗(yàn)條件下得到的效率存在明顯差異,因此應(yīng)根據(jù)具體的實(shí)驗(yàn)數(shù)據(jù)來確定。

    圖3 典型的實(shí)驗(yàn)結(jié)果Fig. 3 Typical experimental results

    而關(guān)于液電效應(yīng)沖擊波峰值壓力的大小一般有近似擬合公式:

    式中:為液電效應(yīng)沖擊波的峰值壓力;主要由傳感器距離決定;α 為常數(shù),根據(jù)實(shí)際的實(shí)驗(yàn)結(jié)果擬合得到。Touya 等根據(jù)實(shí)驗(yàn)數(shù)據(jù)給出了參數(shù)與的關(guān)系,得到:

    式中:的單位為MPa;的單位為mm;的單位為kJ。

    2 水下爆炸等效

    水下爆炸等效主要是基于能量等效的原則,不同的是一種是爆炸能量等效,一種是沖擊波能量等效。

    2.1 爆炸能量等效

    爆炸能量等效主要是基于注入到等離子體通道內(nèi)的能量與炸藥爆炸能量相等的原則來建立相關(guān)關(guān)系,具體關(guān)系式為:

    式中:為炸藥爆炸的能量;為TNT 炸藥的質(zhì)量,kg;為TNT 炸藥的爆炸熱值,一般取值為4.52 MJ/kg。

    根據(jù)該方法,只要測(cè)得放電間隙兩端的電流與電壓值,就可以根據(jù)波形得到液電效應(yīng)注入到間隙內(nèi)的能量,然后求得炸藥質(zhì)量。質(zhì)量給定后,根據(jù)等離子體通道的具體形態(tài),將炸藥設(shè)置成圓柱形或球形,即可通過LS-DYNA 軟件模擬不同爆心距下水下爆炸沖擊波的峰值壓力。關(guān)于水下爆炸沖擊波的研究較多,其中Cole 等、Zamyshlyaev 等關(guān)于無限水中爆炸激波的峰值壓力經(jīng)驗(yàn)公式應(yīng)用較廣:

    式中:為炸藥中心距離測(cè)點(diǎn)的距離,m;為炸藥質(zhì)量,kg;為沖擊波峰值壓力,MPa;為炸藥包半徑。

    假設(shè)為0.472 kJ,通過式(9)計(jì)算得到TNT 炸藥質(zhì)量約為1.044 2×10kg。假定炸藥形態(tài)為球形,則可以建立如圖3 所示的有限元模型。為了減小計(jì)算量,采用1/8 模型,水域的四周設(shè)置無反射邊界,水域的對(duì)稱面設(shè)置法向約束,水域?yàn)殚L(zhǎng)方體, 尺寸為12 cm×12 cm×24 cm,單元使用多物質(zhì)ALE 算法;球形炸藥通過關(guān)鍵字*INITIAL-VOLUME-FRACTION-GEOMETRY 添加到水域中;水域和炸藥采用ALE 單元構(gòu)建。

    圖4 水下爆炸數(shù)值模型Fig. 4 Numerical models of underwater explosion

    TNT 炸藥采用Mat_High_Explosive_Burn 材料模型和JWL 狀態(tài)方程,該方程可以描述為:

    式中:為爆轟產(chǎn)物的壓力;為單位體積的內(nèi)能;為比容,即單位體積裝藥產(chǎn)生的爆轟產(chǎn)物的體積;、、、、ω 為JWL 狀態(tài)方程參數(shù)。具體參數(shù)為:密度ρ=1 630 kg/cm,C-J 爆轟壓力=27 GPa,=374 GPa,=7.33 GPa,=4.15,=0.95,ω=0.3,=7 GJ/m。

    水采用空物質(zhì)材料本構(gòu)模型和Grüneisen 狀態(tài)方程,表達(dá)式為:

    式中:、、為-曲線斜率的系數(shù);為對(duì)一階體積的修正;ρ為初始密度;γ為Grüneisen 狀態(tài)方程參數(shù);=ρ/ρ–1。具體參數(shù)為:ρ=1 020 kg/m,=1 484 m/s,、、、γ、分別為1.979、0、0、0.11、3。

    2.2 沖擊波能量等效

    不同于爆炸能量等效,沖擊波能量等效主要是基于爆炸產(chǎn)生的沖擊波能量與液電效應(yīng)產(chǎn)生的沖擊波能量相等來建立相關(guān)關(guān)系的,具體表達(dá)式如下:

    式中:為液電效應(yīng)產(chǎn)生的沖擊波能量,右端項(xiàng)為爆炸產(chǎn)生的沖擊波能量,單位為kJ。

    根據(jù)該方法,只要利用實(shí)驗(yàn)測(cè)得液中沖擊波的壓力時(shí)程曲線,或者通過式=η=ηη計(jì)算得到?jīng)_擊波的能量,即可得到炸藥質(zhì)量。

    同樣假設(shè)為0.472 kJ,η為28.2%,通過式(13)計(jì)算得到TNT 炸藥質(zhì)量為3.194 5×10kg。就可以利用有限元模型模擬水下爆炸產(chǎn)生沖擊波的過程。

    3 理想氣體等效

    不同于將等離子體通道等效于炸藥,等離子體通道還可以看作絕熱膨脹的理想氣體。則等離子體通道內(nèi)的壓力可以由下面的狀態(tài)方程描述:

    式中:為等離子體通道內(nèi)的壓力;ρ為當(dāng)前等離子體通道內(nèi)的密度;ρ為初始等離子體通道內(nèi)的密度;γ 為絕熱指數(shù),這里取1.25。

    為了在LS-DYNA 中描述等離子體通道的材料與狀態(tài)方程,這里分別使用了關(guān)鍵字*MAT_NULL 和關(guān)鍵字*EOS_LINEAR_POLYNOMIALWITH_ENERGY_LEAK,這樣就可以表征等離子體通道內(nèi)的壓力變化,具體形式為:

    式中:、、、、均為0;==γ-1,使其形式滿足式(14);為初始單位體積內(nèi)能,可以根據(jù)沉積功率時(shí)程曲線來自定義;=ρ/ρ-1;初始等離子體通道密度ρ=1 000 kg/m。

    水的材料與狀態(tài)方程的選取以及有限元模型均與水下爆炸等效部分中的模型相同。為了保證計(jì)算精度的同時(shí)節(jié)省計(jì)算成本,等離子體通道的半徑選擇為1 mm。同樣假設(shè)為0.472 kJ,為了便于輸入到LS-DYNA 中,沉積功率曲線常常簡(jiǎn)化為三角波,如圖5 所示,這樣就可以模擬不同爆心距下的沖擊波壓力。

    圖5 等離子體通道內(nèi)的沉積功率Fig. 5 Electrodeposition power-time curves in the plasma part

    4 分析與討論

    4.1 網(wǎng)格尺寸的選取

    一般而言,基于LS-DYNA 的水下爆炸沖擊波峰值壓力的準(zhǔn)確模擬與所建立的數(shù)值模型密切相關(guān),網(wǎng)格密度、人工黏性以及狀態(tài)方程參數(shù)等因素均會(huì)對(duì)峰值壓力產(chǎn)生影響,其中網(wǎng)格尺寸對(duì)峰壓的影響尤為顯著,網(wǎng)格太大會(huì)影響計(jì)算精度,網(wǎng)格太小則會(huì)產(chǎn)生不必要的計(jì)算成本,因此選擇合適的網(wǎng)格尺寸進(jìn)行數(shù)值模擬至關(guān)重要。本文中主要根據(jù)水下爆炸經(jīng)驗(yàn)公式(式(10))與數(shù)值模擬結(jié)果的擬合程度來確定網(wǎng)格尺寸,值得注意的是,這種確定網(wǎng)格尺寸的方式主要基于爆炸等效方法,而未考慮理想氣體等效方法,但這里仍然給出了網(wǎng)格尺寸對(duì)基于理想氣體等效方法模擬結(jié)果的影響規(guī)律,計(jì)算結(jié)果如圖6 所示。從圖6(a)~(b)中可以看到,隨著網(wǎng)格尺寸減小,不同爆心距下,水下爆炸產(chǎn)生的沖擊波峰壓逐漸增大,與理論公式越來越吻合,說明網(wǎng)格尺寸減小使得計(jì)算精度逐漸提高。當(dāng)網(wǎng)格尺寸從0.10 cm減小到0.08 cm 后,兩種情況下的計(jì)算結(jié)果差別不大,但在實(shí)際的計(jì)算過程中,計(jì)算網(wǎng)格數(shù)目和計(jì)算時(shí)間卻分別增長(zhǎng)了1.95 倍和2.50 倍,因此選擇0.10 cm 的網(wǎng)格尺寸進(jìn)行數(shù)值模擬。不同于水下爆炸等效,基于理想氣體等效方法計(jì)算得到的沖擊波峰壓隨著網(wǎng)格尺寸的減小呈先增大后減小的趨勢(shì),峰壓計(jì)算結(jié)果也比采用水下爆炸等效方法小2 個(gè)數(shù)量級(jí),該特性多半與沉積能量的輸入方式和傳遞效率有關(guān)。

    圖6 不同網(wǎng)格大小下峰值壓力的比較Fig. 6 Comparison of the peak pressures under different mesh sizes

    4.2 不同等效方法的比較

    從圖7 中可以看出,沖擊波以等離子體通道中心為原點(diǎn),以球面波的形式迅速向外傳播,初始?jí)毫^大(呈現(xiàn)紅色),隨著沖擊波向外擴(kuò)散,沖擊波壓力逐漸衰減(逐漸變?yōu)榫G色)。對(duì)比圖7(a) 和(c)、圖7(b)和(c)可以看出,采用爆炸等效方法產(chǎn)生的沖擊波壓力顯著高于理想氣體等效方法產(chǎn)生的沖擊波壓力,約高1~2 個(gè)量級(jí)。通過圖7(a)~(b)可以發(fā)現(xiàn),采用爆炸能量等效方法比采用沖擊波能量等效方法產(chǎn)生的沖擊波壓力值大。而從沖擊波的傳播速度來看,爆炸等效方法產(chǎn)生的沖擊波速度大于理想氣體等效方法,這點(diǎn)可通過=155 μs 時(shí)沖擊波的移動(dòng)距離觀察到。

    為了定量分析采用3 種等效方法模擬沖擊波的效果,給出了距放電間隙中心17 cm 處的沖擊波壓力時(shí)程曲線和距放電中心不同位置處的沖擊波峰值壓力變化曲線,分別如圖8 和圖9 所示。從圖8 中可以看出,采用爆炸等效方法模擬的沖擊波波速較高,峰值更大,衰減更快,而采用理想氣體等效方法模擬的沖擊波則相反,波速較低,峰值小,衰減慢;采用爆炸能量等效方法比采用沖擊波能量等效方法模擬得到的沖擊波峰值壓力高,接近于其2 倍,這主要是由于前者計(jì)算得到的炸藥量更大,根據(jù)式(10)可知,峰值壓力的大小與炸藥質(zhì)量成正相關(guān),炸藥質(zhì)量越大,峰值壓力更高,這些特點(diǎn)也與圖7 觀察到的現(xiàn)象一致;而且可以發(fā)現(xiàn),3 種等效方法模擬得到的波形在衰減階段都存在明顯的振蕩,說明數(shù)值模擬夸大了實(shí)際壓力波動(dòng),這可能是由于數(shù)值模擬中對(duì)水的黏滯性和導(dǎo)熱性造成的能量耗散考慮不足導(dǎo)致的,可引入人工黏性來改進(jìn),但是引入人工黏性又會(huì)引起峰壓下降,本文中主要以沖擊波峰壓作為主要分析因素,因此不再調(diào)整人工黏性系數(shù)。從圖9 中可以看出,隨著與放電中心距離的增大,不管采用哪種方法模擬得到的沖擊波峰值壓力都是逐漸減小的。同一位置處,爆炸能量等效的峰值壓力最高,沖擊波能量等效的峰值壓力次之,理想氣體等效的峰值壓力最低,而且對(duì)比Touya 經(jīng)驗(yàn)公式(式(8)),發(fā)現(xiàn)采用等效爆炸方法模擬出的沖擊波峰值壓力普遍要比Touya 經(jīng)驗(yàn)公式要高,這與文獻(xiàn)[20]中的模擬結(jié)果是一致的。只有當(dāng)爆心距較大時(shí),采用沖擊波能量等效方法模擬得到的峰值壓力才越來越趨近于Touya 經(jīng)驗(yàn)公式計(jì)算得到的峰值壓力。

    圖7 不同等效方法下的沖擊波傳播過程Fig. 7 Propagation process of the shock wave using different equivalent methods

    圖8 不同等效方法下沖擊波壓力時(shí)程曲線Fig. 8 Shock wave pressure-time history curves using different equivalent methods

    圖9 距離放電中心不同位置處的沖擊波峰值壓力Fig. 9 Shock wave peaks at different points from the discharge center

    4.3 方法改進(jìn)與修正

    通過對(duì)比,發(fā)現(xiàn)上述兩類等效方法各有優(yōu)缺點(diǎn),采用爆炸等效,關(guān)鍵在于準(zhǔn)確估計(jì)出炸藥質(zhì)量,而炸藥質(zhì)量的準(zhǔn)確估計(jì)取決于能量轉(zhuǎn)換效率以及等離子體通道內(nèi)沉積能量的準(zhǔn)確計(jì)算。在實(shí)際工程作業(yè)中,由于工況的復(fù)雜性、測(cè)量手段的局限性以及放電擊穿的隨機(jī)性,會(huì)導(dǎo)致沉積能量和效率的計(jì)算誤差較大;采用理想氣體等效方法模擬雖然能夠詳細(xì)地描述等離子體通道內(nèi)的注入能量,在液電成形方向模擬效果良好,但是模擬得到的峰值壓力卻不理想。為此,以在相同爆心距下的峰值壓力相等為基準(zhǔn),對(duì)上述等效方法進(jìn)行改進(jìn)。

    對(duì)于爆炸等效方法,可按照以下步驟進(jìn)行改進(jìn)。

    (1)通過調(diào)整數(shù)值模型的網(wǎng)格大小、人工黏性系數(shù)等參數(shù),使得數(shù)值模擬得到的水下爆炸沖擊波峰值壓力與式(10)較吻合;

    (2)聯(lián)立Touya 經(jīng)驗(yàn)公式(式(8))(如果條件允許也可通過實(shí)驗(yàn)得到水中不同位置處的峰壓與放電能量之間的關(guān)系式來代替)和爆炸經(jīng)驗(yàn)公式(式(10)),得到如下關(guān)系式:

    (3)根據(jù)沖擊波與結(jié)構(gòu)物的作用區(qū)域來確定式(16)中的或。例如:當(dāng)爆炸激波與結(jié)構(gòu)物的作用區(qū)域在5~15 cm 之間時(shí),和的取值范圍分別為0.37~0.57 和0.70~0.82,可分別取0.47 和0.76。

    (4)計(jì)算某一擊穿能量下的炸藥質(zhì)量。求得和,給定擊穿能量即可根據(jù)式(16)計(jì)算得到炸藥質(zhì)量和球形炸藥半徑,例如:=472 J,則=/η≈992 J,計(jì)算得到炸藥質(zhì)量為1.565 8×10kg,炸藥半徑為2.483 2 mm。

    (5)將炸藥半徑輸入到LS-DYNA 中,即可計(jì)算不同爆心距下的沖擊波峰壓。如圖9 所示,可以看到,修改后的等效爆炸模型在10~15 cm 處的峰值壓力略底于Touya 經(jīng)驗(yàn)公式,但是整體十分吻合。

    而對(duì)于理想氣體等效方法,可采取增加輸入沉積功率的方式來進(jìn)行改進(jìn),如圖9 所示。通過大量的試算,發(fā)現(xiàn)當(dāng)輸入到LS-DYNA 中的沉積功率為測(cè)量計(jì)算得到的沉積功率的30 倍,甚至100 倍時(shí),計(jì)算得到的沖擊波峰值壓力有明顯提高,但也依然小于Touya 經(jīng)驗(yàn)公式計(jì)算得到的沖擊波峰壓,而此時(shí)輸入的峰值沉積功率達(dá)到4 056 MW,與實(shí)際差距較大,說明采用理想氣體等效方法模擬沖擊波峰值壓力效果不太理想。

    4.4 沉積能量對(duì)峰壓的影響

    上述模型的建立過程中,僅僅給出了等離子體通道沉積能量為471 J 時(shí)采用不同等效方法模擬沖擊波峰壓的結(jié)果對(duì)比,而沉積能量是決定沖擊波峰壓大小的關(guān)鍵因素之一,因此有必要探究不同沉積能量下各等效方法的結(jié)果差異。依然采用文獻(xiàn)[15-16]中的實(shí)驗(yàn)數(shù)據(jù),利用LS-DYNA 計(jì)算了沉積能量分別為176 和301 J 時(shí)的沖擊波峰壓,計(jì)算結(jié)果如圖10 所示。結(jié)合圖9 和圖10 可以看到,不同沉積能量下,采用爆炸能量等效方法計(jì)算得到的峰值壓力相較于其他等效方法和Touya 經(jīng)驗(yàn)公式總是最高的;隨著沉積能量減小,采用沖擊波能量等效方法計(jì)算得到的峰值壓力經(jīng)歷了先高于、然后趨近、再低于Touya 經(jīng)驗(yàn)公式計(jì)算峰壓的變化過程,這主要是由液電效應(yīng)放電特性決定的,電極間距離減小導(dǎo)致的沉積能量的減小會(huì)導(dǎo)致η減小,但η基本保持不變,而Touya 經(jīng)驗(yàn)公式主要是由放電能量(和η)確定的,沖擊波能量等效方法計(jì)算得到的炸藥質(zhì)量主要由和η確定,在放電能量變化不大時(shí),Touya 經(jīng)驗(yàn)公式計(jì)算得到的峰壓變化也不大,但是由于電極間距離減小導(dǎo)致的沉積能量的減小會(huì)使得沖擊波能量減小,進(jìn)而計(jì)算得到的炸藥質(zhì)量也會(huì)減小,因而會(huì)出現(xiàn)沖擊波能量等效方法計(jì)算得到的峰壓逐漸小于Touya 經(jīng)驗(yàn)公式的現(xiàn)象;采用理想氣體等效方法相比其他等效方法模擬計(jì)算得到的峰壓總是最低的,小1~2 個(gè)數(shù)量級(jí);采用改進(jìn)的等效爆炸方法在不同沉積能量下均能與經(jīng)驗(yàn)公式擬合較好,說明改進(jìn)方法是可行的。

    圖10 不同沉積能量下沖擊波峰值壓力曲線Fig. 10 Shock wave peaks under different deposited energies

    5 結(jié) 論

    為了有效模擬液電效應(yīng)產(chǎn)生的沖擊波以滿足工程需要,介紹了兩種基于LS-DYNA 的數(shù)值模擬方法,水下爆炸等效方法(分為爆炸能量等效與沖擊波能量等效)和理想氣體等效方法,并在等離子體通道沉積能量為472 J 的情況下對(duì)模擬方法進(jìn)行比較與改進(jìn),分析了不同沉積能量下各等效方法峰壓計(jì)算結(jié)果的差異,得到以下主要結(jié)論。

    (1)模擬得到?jīng)_擊波峰值壓力排序?yàn)?,爆炸能量等效的沖擊波峰值壓力最高,沖擊波能量等效的沖擊波峰值壓力次之,理想氣體等效沖擊波峰值壓力最低,理想氣體等效模擬的峰壓相較于前2 種方法小1~2 個(gè)數(shù)量級(jí);沖擊波波速大小排序?yàn)椋芰康刃c沖擊波能量等效的沖擊波波速相等,且高于理想氣體等效的沖擊波波速。

    (2)放電能量變化不大時(shí),電極間距離減小引起的沉積能量減小會(huì)使得3 種等效方法(爆炸能量等效、沖擊波能量等效、理想氣體等效)模擬得到的峰壓均有不同程度的減?。? 種等效方法得到的峰壓的大小順序不發(fā)生變化;改進(jìn)后的等效爆炸方法能夠適應(yīng)沉積能量的變化,與Touya 經(jīng)驗(yàn)公式擬合較好。

    (3)基于LS-DYNA 對(duì)液電效應(yīng)沖擊波峰值壓力進(jìn)行準(zhǔn)確模擬,除了選取合適的等效方法,還應(yīng)結(jié)合具體的放電條件,建立合適的數(shù)值模型,在滿足計(jì)算要求的條件下實(shí)現(xiàn)沖擊波峰值壓力的快速計(jì)算。

    猜你喜歡
    沖擊波炸藥等離子體
    “炸藥”驚魂
    議論火炸藥數(shù)字化制造
    連續(xù)磁活動(dòng)對(duì)等離子體層演化的影響
    基于低溫等離子體修飾的PET/PVC浮選分離
    武漢沖擊波
    中國公路(2019年10期)2019-06-28 03:05:08
    能源物聯(lián)網(wǎng)沖擊波
    能源(2018年10期)2018-12-08 08:02:34
    等離子體種子處理技術(shù)介紹
    醫(yī)生集團(tuán)沖擊波
    超聲雙探頭聯(lián)合定位法在體外沖擊波碎石術(shù)中的應(yīng)用
    Al粉對(duì)炸藥爆炸加速能力的影響
    亚洲av免费高清在线观看| 在线观看午夜福利视频| 97人妻精品一区二区三区麻豆| 国产午夜福利久久久久久| 成年女人看的毛片在线观看| 国产探花在线观看一区二区| 99在线视频只有这里精品首页| 日韩av在线大香蕉| 波多野结衣高清无吗| 岛国在线观看网站| 国产精品亚洲av一区麻豆| 久久久国产成人免费| 午夜福利视频1000在线观看| 观看美女的网站| 国产亚洲欧美在线一区二区| 亚洲片人在线观看| 午夜免费成人在线视频| 日本黄色视频三级网站网址| 欧美最新免费一区二区三区 | 淫秽高清视频在线观看| 国产高清视频在线观看网站| 国产91精品成人一区二区三区| 国产伦精品一区二区三区视频9 | 麻豆国产97在线/欧美| 亚洲国产欧洲综合997久久,| 老汉色av国产亚洲站长工具| 国产伦精品一区二区三区视频9 | a在线观看视频网站| 天天躁日日操中文字幕| 亚洲七黄色美女视频| 日本免费a在线| 国产精品精品国产色婷婷| 成年女人永久免费观看视频| 91av网一区二区| 在线视频色国产色| 97人妻精品一区二区三区麻豆| 色在线成人网| 在线视频色国产色| 成人三级黄色视频| 日韩 欧美 亚洲 中文字幕| 亚洲成人中文字幕在线播放| 欧美+亚洲+日韩+国产| 国产免费一级a男人的天堂| 韩国av一区二区三区四区| 日日夜夜操网爽| 夜夜夜夜夜久久久久| 日本a在线网址| 国产成人系列免费观看| 757午夜福利合集在线观看| 人妻久久中文字幕网| 久久6这里有精品| 国产精品免费一区二区三区在线| 欧美最新免费一区二区三区 | 亚洲自拍偷在线| 久久精品亚洲精品国产色婷小说| 首页视频小说图片口味搜索| 久久久精品欧美日韩精品| 国产精品1区2区在线观看.| 精品人妻一区二区三区麻豆 | 高潮久久久久久久久久久不卡| 日韩高清综合在线| 亚洲最大成人中文| 亚洲av一区综合| 国内毛片毛片毛片毛片毛片| 国内毛片毛片毛片毛片毛片| 午夜福利免费观看在线| 国产国拍精品亚洲av在线观看 | 日日干狠狠操夜夜爽| 一夜夜www| 在线观看美女被高潮喷水网站 | 国产一区二区三区视频了| 日韩精品中文字幕看吧| 欧美日韩综合久久久久久 | 2021天堂中文幕一二区在线观| 好男人电影高清在线观看| 嫩草影视91久久| 精品人妻1区二区| 日韩免费av在线播放| 欧美一区二区精品小视频在线| 日韩精品青青久久久久久| 女人高潮潮喷娇喘18禁视频| 精华霜和精华液先用哪个| 免费av毛片视频| 极品教师在线免费播放| 男插女下体视频免费在线播放| 色视频www国产| 三级毛片av免费| 19禁男女啪啪无遮挡网站| 欧美一区二区国产精品久久精品| 免费人成视频x8x8入口观看| 国产美女午夜福利| 两个人看的免费小视频| 欧美乱码精品一区二区三区| 91久久精品电影网| 亚洲成人中文字幕在线播放| 国产国拍精品亚洲av在线观看 | 亚洲国产精品成人综合色| 露出奶头的视频| 欧美日韩国产亚洲二区| 脱女人内裤的视频| 夜夜躁狠狠躁天天躁| 亚洲在线自拍视频| 午夜免费男女啪啪视频观看 | 亚洲欧美日韩卡通动漫| 欧美乱码精品一区二区三区| 狠狠狠狠99中文字幕| 亚洲乱码一区二区免费版| 男女那种视频在线观看| 好男人在线观看高清免费视频| 精华霜和精华液先用哪个| 国产午夜精品论理片| 在线观看日韩欧美| 色吧在线观看| 亚洲天堂国产精品一区在线| 国产色爽女视频免费观看| 亚洲欧美激情综合另类| 好男人在线观看高清免费视频| 免费看十八禁软件| 亚洲精品久久国产高清桃花| 18禁裸乳无遮挡免费网站照片| 久久6这里有精品| 久久九九热精品免费| av片东京热男人的天堂| 国产成人啪精品午夜网站| 又爽又黄无遮挡网站| 欧美色欧美亚洲另类二区| 精品午夜福利视频在线观看一区| 在线播放国产精品三级| www.www免费av| 亚洲aⅴ乱码一区二区在线播放| 成人特级av手机在线观看| 嫩草影院精品99| 少妇裸体淫交视频免费看高清| 美女高潮喷水抽搐中文字幕| 国产不卡一卡二| h日本视频在线播放| 亚洲一区二区三区不卡视频| 天堂av国产一区二区熟女人妻| 日本在线视频免费播放| 人妻丰满熟妇av一区二区三区| 国产欧美日韩精品一区二区| 美女被艹到高潮喷水动态| 欧美中文综合在线视频| 特大巨黑吊av在线直播| 亚洲欧美激情综合另类| 欧美3d第一页| 岛国在线免费视频观看| 叶爱在线成人免费视频播放| 一区二区三区免费毛片| 亚洲五月婷婷丁香| 亚洲av免费高清在线观看| 一本久久中文字幕| 国产69精品久久久久777片| 欧美色欧美亚洲另类二区| 蜜桃亚洲精品一区二区三区| 九色国产91popny在线| 夜夜躁狠狠躁天天躁| 日日摸夜夜添夜夜添小说| 国产 一区 欧美 日韩| 国产在线精品亚洲第一网站| 国产在线精品亚洲第一网站| 男女那种视频在线观看| 国产精品99久久久久久久久| 欧美+日韩+精品| 国语自产精品视频在线第100页| 国产亚洲精品一区二区www| 日本熟妇午夜| 国产探花极品一区二区| 国产精品98久久久久久宅男小说| 亚洲中文字幕一区二区三区有码在线看| 久久精品91蜜桃| 999久久久精品免费观看国产| 免费av不卡在线播放| 欧美日韩黄片免| 日韩精品青青久久久久久| 亚洲乱码一区二区免费版| 无遮挡黄片免费观看| 少妇丰满av| 国产真实伦视频高清在线观看 | 波多野结衣高清作品| 91麻豆精品激情在线观看国产| 88av欧美| 久久久国产精品麻豆| 18禁在线播放成人免费| 国产高清有码在线观看视频| 嫩草影院入口| 免费看美女性在线毛片视频| 色综合站精品国产| 91久久精品国产一区二区成人 | 亚洲精品国产精品久久久不卡| 亚洲精品色激情综合| 天堂√8在线中文| 国产一区二区亚洲精品在线观看| 久久久久久久久久黄片| 五月玫瑰六月丁香| 嫩草影院精品99| 国产伦人伦偷精品视频| 国产淫片久久久久久久久 | 日本与韩国留学比较| 欧美av亚洲av综合av国产av| 啦啦啦免费观看视频1| 91av网一区二区| 国产私拍福利视频在线观看| 19禁男女啪啪无遮挡网站| 麻豆国产97在线/欧美| 亚洲av日韩精品久久久久久密| 亚洲黑人精品在线| 午夜精品久久久久久毛片777| bbb黄色大片| 亚洲国产精品合色在线| 国产精品国产高清国产av| 亚洲成人免费电影在线观看| 动漫黄色视频在线观看| 在线天堂最新版资源| 日韩欧美 国产精品| 美女免费视频网站| av欧美777| 天堂动漫精品| 精品久久久久久久末码| tocl精华| 亚洲欧美日韩高清专用| 99热6这里只有精品| 亚洲久久久久久中文字幕| 啦啦啦韩国在线观看视频| 少妇的逼水好多| 国产探花极品一区二区| 夜夜看夜夜爽夜夜摸| 日韩人妻高清精品专区| 亚洲18禁久久av| 欧美日韩中文字幕国产精品一区二区三区| 香蕉丝袜av| 欧美成人一区二区免费高清观看| 国产精品嫩草影院av在线观看 | 美女 人体艺术 gogo| 校园春色视频在线观看| 国产精品野战在线观看| 国内精品一区二区在线观看| 日本成人三级电影网站| 国内少妇人妻偷人精品xxx网站| 五月伊人婷婷丁香| 波多野结衣高清作品| 性色av乱码一区二区三区2| 高清毛片免费观看视频网站| 日韩有码中文字幕| 欧美成人免费av一区二区三区| 日韩欧美国产在线观看| 国产成人av激情在线播放| 免费av观看视频| 亚洲成av人片在线播放无| 人妻久久中文字幕网| av天堂在线播放| 午夜精品久久久久久毛片777| 老熟妇乱子伦视频在线观看| 男人的好看免费观看在线视频| 99久久无色码亚洲精品果冻| 午夜福利视频1000在线观看| 精品一区二区三区视频在线 | 国产免费一级a男人的天堂| 亚洲av免费在线观看| 69av精品久久久久久| 法律面前人人平等表现在哪些方面| 岛国在线观看网站| 神马国产精品三级电影在线观看| 国产高潮美女av| 国产欧美日韩精品亚洲av| 久久精品91无色码中文字幕| 久久伊人香网站| 国产伦一二天堂av在线观看| 日本五十路高清| 一二三四社区在线视频社区8| 亚洲五月天丁香| e午夜精品久久久久久久| 久久久久精品国产欧美久久久| 一本综合久久免费| 国产野战对白在线观看| 99精品欧美一区二区三区四区| 制服丝袜大香蕉在线| 精品一区二区三区视频在线 | 麻豆久久精品国产亚洲av| 一级a爱片免费观看的视频| 黄片小视频在线播放| 亚洲av第一区精品v没综合| 亚洲欧美日韩高清专用| 亚洲不卡免费看| 国产精品一区二区三区四区免费观看 | 午夜免费成人在线视频| 婷婷六月久久综合丁香| 观看美女的网站| x7x7x7水蜜桃| 51国产日韩欧美| 亚洲午夜理论影院| 色哟哟哟哟哟哟| 亚洲国产精品999在线| 亚洲无线在线观看| 国产精华一区二区三区| 两人在一起打扑克的视频| 久久精品国产综合久久久| 亚洲国产高清在线一区二区三| 国产高清激情床上av| 岛国视频午夜一区免费看| 久久6这里有精品| 国产精华一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久久久中文| 国产极品精品免费视频能看的| 免费电影在线观看免费观看| 日韩欧美免费精品| 一级作爱视频免费观看| 国产aⅴ精品一区二区三区波| h日本视频在线播放| 女人高潮潮喷娇喘18禁视频| 免费观看人在逋| 99久久精品国产亚洲精品| svipshipincom国产片| 国产亚洲精品av在线| 久久香蕉精品热| 久久99热这里只有精品18| 禁无遮挡网站| 成人18禁在线播放| 757午夜福利合集在线观看| 国产精品久久久久久亚洲av鲁大| 88av欧美| 在线观看免费午夜福利视频| 国产在线精品亚洲第一网站| 亚洲国产精品合色在线| 国产野战对白在线观看| 90打野战视频偷拍视频| 久久九九热精品免费| 高潮久久久久久久久久久不卡| 日韩中文字幕欧美一区二区| 男女下面进入的视频免费午夜| 精品免费久久久久久久清纯| 国产美女午夜福利| 亚洲在线自拍视频| 亚洲欧美日韩东京热| 午夜福利18| 深夜精品福利| 免费搜索国产男女视频| 国产久久久一区二区三区| 手机成人av网站| 成年女人看的毛片在线观看| 国产欧美日韩一区二区三| 亚洲国产精品成人综合色| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产私拍福利视频在线观看| 免费看十八禁软件| 亚洲人成网站在线播| 搡老熟女国产l中国老女人| 日韩欧美国产一区二区入口| 国产精品三级大全| 久99久视频精品免费| 免费观看精品视频网站| 色在线成人网| 床上黄色一级片| 在线观看午夜福利视频| 国产成人啪精品午夜网站| 婷婷亚洲欧美| 88av欧美| a在线观看视频网站| 精品久久久久久成人av| 免费看光身美女| 3wmmmm亚洲av在线观看| 岛国在线免费视频观看| 可以在线观看的亚洲视频| 午夜久久久久精精品| 在线观看66精品国产| 啦啦啦免费观看视频1| 国产伦人伦偷精品视频| 欧美国产日韩亚洲一区| 一个人看的www免费观看视频| 国产不卡一卡二| 免费人成视频x8x8入口观看| 日本一本二区三区精品| 中文资源天堂在线| 亚洲乱码一区二区免费版| 99国产综合亚洲精品| 亚洲欧美激情综合另类| 久久久久久久久中文| 国产一区二区在线观看日韩 | 偷拍熟女少妇极品色| 午夜激情欧美在线| 欧美高清成人免费视频www| 1000部很黄的大片| 精品免费久久久久久久清纯| 国内精品久久久久久久电影| 久久久精品大字幕| 老熟妇仑乱视频hdxx| 99精品欧美一区二区三区四区| 国产亚洲精品久久久com| 亚洲国产欧美网| 欧美+亚洲+日韩+国产| 欧美性猛交黑人性爽| 国产亚洲欧美98| 91九色精品人成在线观看| 91av网一区二区| 欧美在线一区亚洲| 亚洲久久久久久中文字幕| 国产精品一区二区三区四区免费观看 | 熟女人妻精品中文字幕| 夜夜躁狠狠躁天天躁| 日韩人妻高清精品专区| 丰满乱子伦码专区| 成熟少妇高潮喷水视频| 国产精品亚洲美女久久久| 欧美bdsm另类| 男女视频在线观看网站免费| 特大巨黑吊av在线直播| 我的老师免费观看完整版| 变态另类丝袜制服| 18禁黄网站禁片午夜丰满| 脱女人内裤的视频| 久久精品国产自在天天线| 日韩欧美精品免费久久 | 欧美日韩一级在线毛片| 啦啦啦观看免费观看视频高清| 久久精品人妻少妇| 90打野战视频偷拍视频| 高清在线国产一区| 国产一区二区三区视频了| 舔av片在线| 日本一本二区三区精品| 欧美日韩黄片免| 19禁男女啪啪无遮挡网站| 国内精品久久久久久久电影| 日本成人三级电影网站| 99热只有精品国产| 真人一进一出gif抽搐免费| 9191精品国产免费久久| 久久中文看片网| 午夜老司机福利剧场| 免费一级毛片在线播放高清视频| 午夜福利在线在线| www.999成人在线观看| 美女cb高潮喷水在线观看| 一本精品99久久精品77| 亚洲成人久久性| 亚洲内射少妇av| 99riav亚洲国产免费| 成年免费大片在线观看| 欧美日韩国产亚洲二区| 18禁国产床啪视频网站| 女警被强在线播放| 熟妇人妻久久中文字幕3abv| 精品久久久久久,| 日本 欧美在线| 好看av亚洲va欧美ⅴa在| 日韩欧美精品v在线| 亚洲黑人精品在线| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 午夜福利在线在线| 波多野结衣巨乳人妻| 亚洲乱码一区二区免费版| 亚洲性夜色夜夜综合| 国产v大片淫在线免费观看| 欧美日韩精品网址| 国产精品自产拍在线观看55亚洲| 中出人妻视频一区二区| 波野结衣二区三区在线 | 国产真人三级小视频在线观看| www.www免费av| 国产精品电影一区二区三区| 国产一区二区三区视频了| 嫩草影院精品99| 一卡2卡三卡四卡精品乱码亚洲| 在线播放无遮挡| 日韩人妻高清精品专区| 一个人免费在线观看的高清视频| 最近最新中文字幕大全电影3| 国产99白浆流出| 长腿黑丝高跟| 在线看三级毛片| 欧美日韩综合久久久久久 | 99热只有精品国产| 欧美成狂野欧美在线观看| 99热精品在线国产| 亚洲片人在线观看| 1024手机看黄色片| 亚洲国产日韩欧美精品在线观看 | 国产欧美日韩一区二区精品| 最新美女视频免费是黄的| 中文字幕av成人在线电影| 欧美激情在线99| 国产激情偷乱视频一区二区| 亚洲欧美日韩无卡精品| 色精品久久人妻99蜜桃| 夜夜躁狠狠躁天天躁| 人人妻,人人澡人人爽秒播| 亚洲av日韩精品久久久久久密| 欧美中文日本在线观看视频| 国产乱人视频| 国产一区二区亚洲精品在线观看| 久久精品国产99精品国产亚洲性色| 小说图片视频综合网站| 欧美在线一区亚洲| 日韩大尺度精品在线看网址| 久久精品人妻少妇| 亚洲五月天丁香| 久久久久精品国产欧美久久久| 日韩成人在线观看一区二区三区| 日本a在线网址| 国产伦一二天堂av在线观看| 国产精品久久久久久精品电影| 久久精品91无色码中文字幕| 午夜视频国产福利| 深夜精品福利| 国产单亲对白刺激| 在线a可以看的网站| 国产一区二区三区在线臀色熟女| 蜜桃久久精品国产亚洲av| 一个人看的www免费观看视频| 综合色av麻豆| 欧美一级a爱片免费观看看| 亚洲专区中文字幕在线| 变态另类丝袜制服| 人妻久久中文字幕网| 2021天堂中文幕一二区在线观| 久久久久久久精品吃奶| 中文字幕久久专区| 国产精品av视频在线免费观看| 一本一本综合久久| 好看av亚洲va欧美ⅴa在| 亚洲av免费在线观看| 国产精品永久免费网站| 在线a可以看的网站| 国产精品精品国产色婷婷| 99热这里只有精品一区| 成人三级黄色视频| 天天躁日日操中文字幕| 波多野结衣巨乳人妻| 亚洲激情在线av| 色尼玛亚洲综合影院| 1000部很黄的大片| 欧美zozozo另类| 久久精品国产综合久久久| 午夜免费观看网址| 成人av在线播放网站| 久久国产精品影院| 在线观看美女被高潮喷水网站 | 亚洲 国产 在线| 在线观看免费午夜福利视频| 三级男女做爰猛烈吃奶摸视频| 露出奶头的视频| 国产伦精品一区二区三区四那| 国产伦精品一区二区三区视频9 | 蜜桃亚洲精品一区二区三区| 欧美三级亚洲精品| 变态另类丝袜制服| 成人欧美大片| 香蕉av资源在线| 老司机福利观看| 精品人妻偷拍中文字幕| 偷拍熟女少妇极品色| 亚洲欧美日韩高清在线视频| 久久婷婷人人爽人人干人人爱| 欧美3d第一页| 91久久精品电影网| 51国产日韩欧美| 国产精品久久久久久久电影 | av片东京热男人的天堂| 大型黄色视频在线免费观看| 黄片大片在线免费观看| 伊人久久精品亚洲午夜| 亚洲成人精品中文字幕电影| 给我免费播放毛片高清在线观看| 91在线观看av| 国产成人福利小说| 亚洲avbb在线观看| 久久久成人免费电影| 亚洲第一欧美日韩一区二区三区| 日韩欧美国产一区二区入口| 99riav亚洲国产免费| 国产精品亚洲一级av第二区| 国产一区二区在线av高清观看| 色av中文字幕| 男女午夜视频在线观看| 97超级碰碰碰精品色视频在线观看| 国产一级毛片七仙女欲春2| 日本免费a在线| 亚洲国产欧美网| 国内毛片毛片毛片毛片毛片| 亚洲欧美日韩高清在线视频| 少妇人妻一区二区三区视频| av女优亚洲男人天堂| 成人国产一区最新在线观看| 一本一本综合久久| 在线观看66精品国产| 舔av片在线| 久久精品国产亚洲av涩爱 | 18禁美女被吸乳视频| 波多野结衣巨乳人妻| 深爱激情五月婷婷| 两个人看的免费小视频| 草草在线视频免费看| 亚洲av电影不卡..在线观看| 亚洲真实伦在线观看| 亚洲国产色片| 伊人久久精品亚洲午夜| 香蕉av资源在线| 欧美一级毛片孕妇| 一个人看的www免费观看视频| 亚洲欧美日韩卡通动漫| 午夜日韩欧美国产| 长腿黑丝高跟| 99热这里只有精品一区| 麻豆久久精品国产亚洲av| 午夜福利在线观看吧| 久久久久国产精品人妻aⅴ院| 在线观看美女被高潮喷水网站 | av天堂在线播放| 性色av乱码一区二区三区2| 少妇熟女aⅴ在线视频| 国产精品日韩av在线免费观看| 脱女人内裤的视频| 国产主播在线观看一区二区| 欧美3d第一页|