王英杰，閆鳳國(guó)，陳建營(yíng)，薛海恩，梁基冠

(1.陸軍工程大學(xué),江蘇 南京 210007; 2.中山大學(xué),廣東 珠海 519082)

1 概述

當(dāng)核武器在地面發(fā)生爆炸時(shí),核爆炸產(chǎn)生的能量分為兩種形式傳播,一是爆炸產(chǎn)生的能量直接耦合入土地面介質(zhì)中并形成直接的沖擊,是地下工程所要考慮的主要荷載;另外一種是以空氣沖擊波的形式在空氣中傳播,并會(huì)在地表面區(qū)域的介質(zhì)中耦合產(chǎn)生感生沖擊波[1-2]。同時(shí)觸地爆炸時(shí)還需要考慮地面的自由表面作用,而自由表面效應(yīng)和介質(zhì)強(qiáng)度的綜合效應(yīng),會(huì)導(dǎo)致接近自由表面處的峰值應(yīng)力和峰值速度低于下方軸線上等徑向距離處的峰值應(yīng)力與峰值速度。因此,不同于在空氣中爆炸的球?qū)ΨQ傳播模式,核武器觸地爆炸時(shí)沖擊波傳播較為復(fù)雜。本文首先通過(guò)數(shù)值模擬結(jié)果與經(jīng)驗(yàn)公式對(duì)比的方式,用TrueGrid軟件和LS-DYNA軟件建立更加復(fù)雜的工況并進(jìn)行模擬分析,分析地面材料為多種介質(zhì)時(shí),在觸地爆炸荷載作用下介質(zhì)中應(yīng)力波的傳播過(guò)程。

2 觸地爆效應(yīng)計(jì)算的經(jīng)驗(yàn)公式

目前對(duì)觸地爆的研究相對(duì)較少,但在少量的文獻(xiàn)中,依然可以找到相關(guān)的經(jīng)驗(yàn)公式,并且觸地爆的經(jīng)驗(yàn)公式總是與地下封閉爆炸存在一定的聯(lián)系,一般情況下地下爆炸的經(jīng)驗(yàn)公式是炸藥的質(zhì)量乘上某一折減系數(shù)。文獻(xiàn)[3]中就對(duì)觸地爆的經(jīng)驗(yàn)公式做了較為詳細(xì)的闡述,如沃克夫、徹爾內(nèi)、斯米爾諾夫、廖哈夫、克拉維茨、萊普遜和其他的學(xué)者給定了關(guān)于土在爆炸荷載作用下的相關(guān)參數(shù)。在本文中主要采用廖哈夫和克拉維茨提出的經(jīng)驗(yàn)公式進(jìn)行對(duì)比,分別表示為式(1)，式(2):

(1)

(2)

其中,Pi為爆炸沖擊波在土中形成的超壓,kg/cm2,約等于0.1 MPa;R為炸藥中心到測(cè)試點(diǎn)的距離,m;W為炸藥的質(zhì)量,kg。

另外,克拉維茨還提出了觸地爆時(shí),土中最大質(zhì)點(diǎn)速度的經(jīng)驗(yàn)公式,如式(3)所示：

(3)

3 材料本構(gòu)模型

3.1 空氣材料模型

空氣的成分復(fù)雜,在計(jì)算時(shí)我們需要根據(jù)空氣的特性將其簡(jiǎn)化。在大當(dāng)量的爆炸情況下,我們可以將空氣理想化為非黏性理想氣體,選用9號(hào)材料*MAT-NULL[4]來(lái)作為空氣的本構(gòu)模型,為了描述沖擊波壓力與空氣初始內(nèi)能密度的關(guān)系,我們引入狀態(tài)方程*EOS-LINEAR-POLYNOMIAL,如式(4)，式(5)所示:

Pa=C0+C1u+C2u2+C3u3+(C4+C5u+C6u2)Ea0

(4)

(5)

其中，Pa為空氣沖擊波壓力;Ea0為空氣初始單位體積內(nèi)能;ρ為擴(kuò)散密度;ρ0為標(biāo)準(zhǔn)密度;C0,C1,C2,C3,C4,C5,C6均為計(jì)算常數(shù),當(dāng)該狀態(tài)方程描述空氣時(shí),C4和C5可由式(6)計(jì)算:

C4=C5=λ-1

(6)

其中，λ為空氣的絕熱指數(shù);空氣的各參數(shù)取值見(jiàn)表1。

表1 空氣的主要計(jì)算參數(shù)

3.2 炸藥材料模型

由于重核裂變和輕核聚變等核武器的鏈?zhǔn)椒磻?yīng)的復(fù)雜性,目前國(guó)內(nèi)外關(guān)于沖擊波效應(yīng),特別是核武器產(chǎn)生的沖擊波效應(yīng)大多采用等效TNT當(dāng)量來(lái)近似替代[5]。因此我們用8號(hào)材料*MAT-HIGH-EX·PLOSIVE-BURN作為等效TNT的本構(gòu)模型,并通過(guò)引入狀態(tài)方程*EOS-JWL用于表示TNT爆轟波陣面壓力與初始相對(duì)體積V0、單位體積內(nèi)能E0之間的關(guān)系,如式(7)所示:

(7)

其中，P為爆轟壓力;V0為初始相對(duì)體積;E0為初始單位體積內(nèi)能;A，ω，R1，B，R2均為材料常數(shù)。炸藥的各參數(shù)取值如表2所示。

表2 炸藥的主要計(jì)算參數(shù)

3.3 巖石材料模型

國(guó)內(nèi)外對(duì)于巖石在大當(dāng)量爆炸理論的研究較少,為了方便計(jì)算,經(jīng)過(guò)綜合考慮，我們將巖石視為彈塑性材料,因?yàn)?號(hào)模型*MAT-PLASTIC-KINEMATIC能充分考慮巖石的應(yīng)變率影響,適合用于描述巖石材料模型的各向同性硬化和隨動(dòng)硬化塑性[6]。因此我們?cè)谀M中選用3號(hào)模型*MAT-PLASTIC-KINEMATIC作為巖石的材料模型。在爆炸荷載的作用下材料的應(yīng)變速率的影響不可忽略,材料模型的屈服強(qiáng)度可分為初始屈服強(qiáng)度和硬化部分強(qiáng)度,屈服強(qiáng)度可按式(8)計(jì)算:

(8)

(9)

表3 巖石的主要計(jì)算參數(shù)

3.4 土壤材料模型

土壤的成分復(fù)雜,在受到?jīng)_擊后會(huì)被迅速壓縮,為了能較好的模擬出土壤在受到?jīng)_擊波后的壓縮性能,我們一般都用5號(hào)材料模型*MAT-SOIL-AND-FOAM來(lái)作為土壤的本構(gòu)模型。并可通過(guò)定義一些參數(shù)以考慮壓力與體積應(yīng)變關(guān)系中熱效應(yīng)對(duì)其所產(chǎn)生的影響。該材料模型的屈服函數(shù)如式(10)所示:

φ=J2-(a0+a1p+a2p2)

(10)

其中，a0，a1，a2均為材料常數(shù),a0考慮了土壤摩擦角的影響,a1考慮了土壤黏聚力的影響,a2考慮土壤受動(dòng)荷載效應(yīng)的影響;p為靜態(tài)壓力。該屈服函數(shù)中引入第二個(gè)不變量J2來(lái)描述:

(11)

Sij=σij-δijp

(12)

(13)

其中，Sij為偏應(yīng)力張量;σij為應(yīng)力張量;δij為克羅內(nèi)克函數(shù)。土壤的各參數(shù)取值如表4所示[7]。

表4 土壤的主要計(jì)算參數(shù)

4 觸地爆的有限元模型

為了方便建模及提高計(jì)算效率,觸地爆采用1/4對(duì)稱模型,炸藥采用等效TNT方形炸藥,尺寸為85 m×85 m×85 m,所等效的TNT當(dāng)量約為100萬(wàn)t。建模過(guò)程結(jié)合TrueGrid軟件和LS-DYNA軟件,分為兩個(gè)部分,分別是空氣和炸藥、土壤與巖石。兩個(gè)部分均采用TrueGrid軟件建模,再將兩個(gè)部分分別導(dǎo)入LS-DYNA軟件中形成完整模型,如圖1所示。圖1(a)是炸藥與空氣模型,整體尺寸為2 600 m×1 300 m,圖1(b)是土、巖混合介質(zhì)模型,整體尺寸為1 200 m×1 200 m。該模型均采用Solid164單元,土與巖石接觸面采用自動(dòng)接觸。在建模過(guò)程中,炸藥與空氣的單元網(wǎng)格均采用共節(jié)點(diǎn)的方式建模[8-9]。

5 純土介質(zhì)中的觸地爆效應(yīng)

為了驗(yàn)證數(shù)值模擬方式的準(zhǔn)確性,首先開(kāi)展純土介質(zhì)中核武器觸地爆的數(shù)值模擬,用于與經(jīng)驗(yàn)公式給出的觸地爆炸的結(jié)果進(jìn)行對(duì)比。由于與經(jīng)驗(yàn)公式進(jìn)行對(duì)比的數(shù)值模擬的純土介質(zhì)的有限元模型與章節(jié)4中所述模型基本相同,僅將土、巖石全部由土代替,因此不再贅述模型的建立過(guò)程,直接對(duì)數(shù)值模擬結(jié)果進(jìn)行對(duì)比分析。

5.1 土中應(yīng)力波的傳播過(guò)程分析

圖2呈現(xiàn)了土中應(yīng)力波傳遞的大致過(guò)程。在爆炸后的初始階段(0.06 s之前),應(yīng)力波近似以半球形對(duì)稱的方式向四周傳播,沿地表面方向和沿徑向深度方向的應(yīng)力波峰值差別不大,但從應(yīng)力波傳播區(qū)域來(lái)看,沿縱向深度方向的應(yīng)力波傳播明顯大于沿地表面?zhèn)鞑シ较?可以發(fā)現(xiàn)在觸地爆炸過(guò)程中更多的能量是直接耦合入土介質(zhì)中;隨著沖擊波的傳播(0.06 s～0.14 s),沖擊波波陣面以近似橢球型的方式繼續(xù)向四周傳播,沿著縱向深度方向的沖擊波傳播較為規(guī)律,沖擊波超壓峰值不斷衰減;但由于地表面的自由表面效應(yīng)和介質(zhì)衰減的綜合效應(yīng),地表面的沖擊波傳播過(guò)程較為復(fù)雜,會(huì)出現(xiàn)沖擊波的間斷和不規(guī)律傳播;在0.14 s以后,沖擊波傳播到所建立的模型之外的區(qū)域,在目標(biāo)土介質(zhì)區(qū)域,應(yīng)力波主要傳播至較深深度處,另外,地表面處由于空氣沖擊波影響而產(chǎn)生的感生地沖擊使得地表面處的應(yīng)力峰值相對(duì)較大,并且土層表面在爆炸荷載的作用下,會(huì)形成一個(gè)圓形的爆坑,爆坑的直徑隨著時(shí)間的增長(zhǎng)逐漸增大,在0.4 s后逐漸趨向于穩(wěn)定。

5.2 土中應(yīng)力波的衰減過(guò)程分析

為了更加清楚地觀測(cè)應(yīng)力波在土層中的衰減過(guò)程,并且對(duì)比數(shù)值模擬結(jié)果的準(zhǔn)確性,利用origin軟件對(duì)土層中不同比例距離處的入射超壓變化曲線進(jìn)行曲線擬合。圖3(a)中表示了土層中入射超壓在不同比例距離的情況下隨時(shí)間變化的過(guò)程,我們可以發(fā)現(xiàn),入射超壓在土層中傳遞時(shí),呈快速遞減到緩慢遞減的一個(gè)過(guò)程。結(jié)合圖3(b),當(dāng)比例距離為0.2 m/kg1/3～0.4 m/kg1/3時(shí),超壓峰值相較于其他段遞減趨勢(shì)更加明顯,這就更加反映了在近距離情況下爆炸荷載作用衰減較快,且其過(guò)程較為復(fù)雜。圖3(b)還對(duì)比了經(jīng)驗(yàn)公式計(jì)算得到的不同比例距離處的入射超壓峰值與數(shù)值模擬結(jié)果:從整體趨勢(shì)看,數(shù)值模擬計(jì)算得到的結(jié)果與廖哈夫和克勞維茨經(jīng)驗(yàn)公式的整體趨勢(shì)均較為吻合,但當(dāng)比例距離小于0.4 m/kg1/3時(shí),經(jīng)驗(yàn)公式計(jì)算得到的入射超壓峰值均稍小于數(shù)值模擬得到的結(jié)果,這可能是由于我們采用的土層材料模型的密度與經(jīng)驗(yàn)公式的材料模型的密度有所差別。

5.3 質(zhì)點(diǎn)運(yùn)動(dòng)速度分析

在對(duì)地表面中爆炸引起的動(dòng)力響應(yīng)進(jìn)行研究時(shí),質(zhì)點(diǎn)的運(yùn)動(dòng)速度也是一個(gè)需要關(guān)注的重要變量,因此除了對(duì)土層中的入射超壓峰值進(jìn)行驗(yàn)證,我們還對(duì)比分析了不同比例距離深度處測(cè)點(diǎn)的運(yùn)動(dòng)速度。圖4(a)給出了選定的不同深度土層中質(zhì)點(diǎn)運(yùn)動(dòng)速度變化的情況,以某一深度土層質(zhì)點(diǎn)來(lái)看:其質(zhì)點(diǎn)運(yùn)動(dòng)速度表現(xiàn)為快速上升并緩慢下降的過(guò)程,這是因?yàn)楸_擊波到達(dá)土層中某一點(diǎn)時(shí),會(huì)引起該點(diǎn)處的速度瞬間發(fā)生變化,當(dāng)沖擊波離開(kāi)以后,該點(diǎn)處速度逐漸減小,而由于土層中距離炸藥越近的部分產(chǎn)生的塑性變形越大,因此土層中某一點(diǎn)處的速度減小時(shí),會(huì)受到前面土層的阻礙,消耗的時(shí)間也相對(duì)增長(zhǎng)。結(jié)合圖4(a)、圖4(b)也可以發(fā)現(xiàn),隨著比例距離深度的增加,其最大質(zhì)點(diǎn)速度也迅速衰減,當(dāng)比例距離達(dá)到0.6 m/kg1/3之后,其最大質(zhì)點(diǎn)速度逐漸保持穩(wěn)定。由于現(xiàn)有的觸地爆研究資料較少,關(guān)于觸地爆在土中引起的最大質(zhì)點(diǎn)速度成熟的經(jīng)驗(yàn)公式也較少,因此在圖4(b)中僅用克拉維茨提出的最大質(zhì)點(diǎn)速度經(jīng)驗(yàn)公式作為對(duì)比參考。從圖4中可以看出,數(shù)值模擬得到的最大質(zhì)點(diǎn)速度與經(jīng)驗(yàn)公式的計(jì)算結(jié)果吻合較好。

綜合土介質(zhì)中爆炸應(yīng)力波的傳播過(guò)程和對(duì)應(yīng)力波超壓峰值、質(zhì)點(diǎn)最大速度的分析來(lái)看,數(shù)值模擬得到的應(yīng)力波傳播過(guò)程和經(jīng)驗(yàn)研究較為符合,且計(jì)算得到的應(yīng)力比超壓峰值和質(zhì)點(diǎn)最大速度與經(jīng)驗(yàn)公式擬合較好,驗(yàn)證了所采用的數(shù)值模擬方法的合理性和數(shù)值模型參數(shù)的準(zhǔn)確性。

6 土、巖混合介質(zhì)中的觸地爆效應(yīng)

從上述純土介質(zhì)中的結(jié)果對(duì)比分析來(lái)看,此數(shù)值模擬方法對(duì)觸地爆的研究分析的結(jié)果是合理的。所以,我們進(jìn)一步將該數(shù)值模擬方法推廣至土、巖混合介質(zhì)的情形,并分析爆炸荷載作用下土、巖混合介質(zhì)的動(dòng)力響應(yīng)。土層與巖層之間的接觸均采用自動(dòng)面面接觸。

6.1 土、巖混合介質(zhì)中應(yīng)力波的傳播過(guò)程分析

圖5描述了爆炸沖擊波在多種介質(zhì)中的傳遞情況。從傳播過(guò)程可以看出,觸地爆炸產(chǎn)生的應(yīng)力波整體傳播規(guī)律和在純土介質(zhì)中的傳播規(guī)律差別不大,主要呈現(xiàn)出以下幾個(gè)特點(diǎn):爆炸產(chǎn)生的沖擊波同時(shí)沿著深度方向和地表面方向傳播,且主要沿深度方向傳播;在沿地表面深度方向,依然表現(xiàn)出入射超壓遞減的趨勢(shì)和沿地表面?zhèn)鞑シ较?會(huì)受到自由表面效應(yīng)的影響;但和純土介質(zhì)不同的是,由于不同介質(zhì)性質(zhì)的差異,在應(yīng)力波傳播過(guò)程中在土層與巖層的接觸面處會(huì)存在應(yīng)力積累,這是因?yàn)橥翆优c巖層在爆炸荷載作用下產(chǎn)生的塑性變形不連續(xù),造成了該現(xiàn)象的產(chǎn)生,使得混合介質(zhì)中應(yīng)力波的分布結(jié)果相對(duì)于純土介質(zhì)較為復(fù)雜。

圖6顯示了混合地質(zhì)的質(zhì)點(diǎn)速度隨時(shí)間變化的云圖?？梢钥闯?起爆后,和炸藥接觸的土壤層的質(zhì)點(diǎn)速度迅速上升,在比例爆距0.1 m/kg1/3可達(dá)到接近100 m/s,隨著離地越遠(yuǎn),質(zhì)點(diǎn)速度逐漸下降,當(dāng)碰到地下巖層的時(shí)候,質(zhì)點(diǎn)速度開(kāi)始出現(xiàn)急劇下降,這是因?yàn)閹r層的彈性模量較大,其相對(duì)于土層堅(jiān)硬了許多,因此也導(dǎo)致了沖擊波經(jīng)過(guò)巖層時(shí)質(zhì)點(diǎn)速度發(fā)生銳減。隨著距離地面深度的增加,在沖擊波傳至第二層深埋土?xí)r,其最大質(zhì)點(diǎn)速度已經(jīng)衰減至1 m/s以下。

6.2 土、巖混合介質(zhì)中應(yīng)力波的衰減過(guò)程分析

圖7中給出了應(yīng)力波在土層與巖層混合介質(zhì)地層中傳播的整體過(guò)程,并將數(shù)值模擬得到的入射超壓峰值與純土介質(zhì)的計(jì)算結(jié)果進(jìn)行對(duì)比。從整體上看,在多種介質(zhì)混合的情況下,入射超壓在傳播過(guò)程中會(huì)發(fā)生較大的變化,但總體的應(yīng)力衰減較慢,特別是在比例距離為0.202 m/kg1/3時(shí),該深度處的應(yīng)力水平會(huì)長(zhǎng)時(shí)間保持較高的水平,這是因?yàn)樵撋疃惹『迷谕翆优c巖層的接觸位置附近。從圖7(b)中我們可以發(fā)現(xiàn),在土層、巖層混合的地質(zhì)中,不同比例距離處的質(zhì)點(diǎn)的入射超壓峰值均大于純土地質(zhì)中的相應(yīng)結(jié)果,并且這一現(xiàn)象在比例距離較小時(shí)更為顯著,結(jié)合壓力傳播云圖來(lái)看,這也是由于不同巖層接近位置處的應(yīng)力積累造成的。綜合以上對(duì)入射超壓的分析,可以得出,混合地質(zhì)中的入射超壓比單一地質(zhì)衰減要慢,其超壓峰值相對(duì)較高,且在不同巖層的接觸位置會(huì)出現(xiàn)應(yīng)力積累的現(xiàn)象。

6.3 土、巖混合介質(zhì)中質(zhì)點(diǎn)速度和介質(zhì)變形分析

圖8給出了土巖混合介質(zhì)工況中不同深度的質(zhì)點(diǎn)速度及與單一介質(zhì)計(jì)算結(jié)果的對(duì)比,可以明顯看出,混合地質(zhì)中質(zhì)點(diǎn)速度的衰減速率要遠(yuǎn)大于在純土地質(zhì)中質(zhì)點(diǎn)速度的衰減。由于選點(diǎn)多集中于巖層中,并結(jié)合對(duì)入射超壓的描述,可以得出以下結(jié)論:應(yīng)力波在巖石中的傳播速度要大于在土介質(zhì)中的傳播速率,并且在巖石中質(zhì)點(diǎn)速度衰減較快,這一現(xiàn)象在比例距離較小時(shí)差別較為顯著,而隨著比例距離的增加,二者差別逐漸減小。

圖9是混合地質(zhì)的質(zhì)點(diǎn)位移隨時(shí)間變化的云圖。從中可以看出,其變化趨勢(shì)和質(zhì)點(diǎn)速度變化趨勢(shì)大致相同,起爆后,炸藥產(chǎn)生的壓縮波迅速向地層傳播,與炸藥表面接觸的土層位移較大,當(dāng)沖擊波傳遞至巖層時(shí),質(zhì)點(diǎn)位移出現(xiàn)明顯的下降。從質(zhì)點(diǎn)速度和質(zhì)點(diǎn)位移的云圖來(lái)看,到最終應(yīng)力波破壞效果達(dá)到穩(wěn)定后,土介質(zhì)整體位移較大,地面形成爆坑,而巖石介質(zhì)中開(kāi)始爆炸導(dǎo)致的質(zhì)點(diǎn)位移有部分逐漸恢復(fù),綜合來(lái)看,土巖混合巖層的不同介質(zhì)差異對(duì)爆炸沖擊波會(huì)產(chǎn)生的明顯削弱作用[10]。

圖10是混合地質(zhì)的有效塑性應(yīng)變隨時(shí)間變化的云圖?？梢钥闯?土壤和巖層的有效塑性應(yīng)變變化情況截然不同,與炸藥接觸的土層,僅在與炸藥接觸的較小范圍內(nèi)有明顯的塑性應(yīng)變,而與爆炸中心水平距離較遠(yuǎn)處的塑性應(yīng)變不大,這是由于觸地爆的沖擊波主要是沿著地層垂直向下傳播的原因;在巖層中,也僅與第一層土層接觸的200 m范圍內(nèi)有較明顯的塑性應(yīng)變,隨著深度的增加,沖擊波在巖層中所產(chǎn)生的有效塑性應(yīng)變發(fā)生驟減。整體來(lái)看,僅在觸地爆較近的位置處有一定的等效塑性應(yīng)變,而當(dāng)距離較遠(yuǎn)時(shí),在地介質(zhì)中幾乎無(wú)塑性變形產(chǎn)生。

7 結(jié)語(yǔ)

1)通過(guò)LS-DYNA軟件,對(duì)核武器觸地爆的過(guò)程進(jìn)行了數(shù)值模擬,并將數(shù)值模擬得到的結(jié)果與前蘇聯(lián)廖哈夫和克拉維茨等人根據(jù)實(shí)驗(yàn)數(shù)據(jù)擬合提出的經(jīng)驗(yàn)公式進(jìn)行對(duì)比,數(shù)值模擬得到的結(jié)果和經(jīng)驗(yàn)公式得到的結(jié)果整體規(guī)律是一致的,數(shù)據(jù)結(jié)果差別不大,驗(yàn)證了數(shù)值模擬方法的可靠性。

2)開(kāi)展了純土介質(zhì)中核武器觸地爆的數(shù)值模擬,分析了土介質(zhì)中爆炸應(yīng)力波的傳播過(guò)程和對(duì)應(yīng)力波超壓峰值、質(zhì)點(diǎn)最大速度等,數(shù)值模擬得到的應(yīng)力波傳播過(guò)程和經(jīng)驗(yàn)研究較為符合,且計(jì)算得到的應(yīng)力比超壓峰值和質(zhì)點(diǎn)最大速度與經(jīng)驗(yàn)公式擬合較好,驗(yàn)證了所采用的數(shù)值模擬方法的合理性和數(shù)值模型參數(shù)的準(zhǔn)確性。

3)通過(guò)對(duì)觸地爆工況下所得實(shí)驗(yàn)結(jié)果的驗(yàn)證,建立了土巖介質(zhì)混合的模型并對(duì)不同的工況結(jié)果進(jìn)行了對(duì)比分析,通過(guò)對(duì)應(yīng)力波傳播過(guò)程的分析以及對(duì)爆炸沖擊波所產(chǎn)生的峰值超壓、質(zhì)點(diǎn)最大速度、最大位移、有效塑性應(yīng)變分布等特征參數(shù)的分析比較來(lái)看,在地下設(shè)置不同特性的巖層,會(huì)在巖層交界處存在應(yīng)力積累線性,但其最大質(zhì)點(diǎn)速度和質(zhì)點(diǎn)位移均會(huì)衰減較快,因此發(fā)現(xiàn)不同巖層混合能夠大幅降低爆炸產(chǎn)生的沖擊波的動(dòng)態(tài)破壞效應(yīng)。