李 杰，程怡豪,2，徐天涵，王明洋,3

（1. 陸軍工程大學(xué)爆炸沖擊防災(zāi)減災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室，江蘇 南京 210007；2. 陸軍工程大學(xué)野戰(zhàn)工程學(xué)院，江蘇 南京 210007；3. 南京理工大學(xué)機(jī)械工程學(xué)院，江蘇 南京 210094）

侵徹效應(yīng)是指彈體高速撞擊靶體產(chǎn)生的破壞效應(yīng)，若彈體穿透金屬靶體，人們又習(xí)慣稱為穿甲。天然巖石具有強(qiáng)度高、抗侵徹性能好的特點(diǎn)，因此重要的地下防護(hù)工程總是盡可能選址于強(qiáng)度高、完整性好的天然巖體中，許多防護(hù)工程也會(huì)利用塊石作為遮彈層的主要材料[1-3]。因此，彈體對(duì)巖石侵徹效應(yīng)歷來是防護(hù)工程領(lǐng)域關(guān)注的關(guān)鍵科學(xué)問題之一，其往往體現(xiàn)為相互“矛盾”的兩方面：一方面，武器裝備研究人員通過提升彈體速度、優(yōu)化彈體結(jié)構(gòu)和發(fā)展串聯(lián)戰(zhàn)斗部等手段增加侵徹深度[4]；另一方面，防護(hù)工程研究人員則需要合理設(shè)計(jì)以天然巖層為基礎(chǔ)的安全防護(hù)層厚度，通過增加掩體埋設(shè)深度、改進(jìn)防護(hù)層材料和配置方案等手段提高防護(hù)工程的抗侵徹能力[5-7]。

侵徹效應(yīng)的研究最早可以追溯到19 世紀(jì)以前[1,8]，在20 世紀(jì)取得長(zhǎng)足進(jìn)展。尤其最近三十年，因GBU-28 等精確制導(dǎo)鉆地武器在幾場(chǎng)局部戰(zhàn)爭(zhēng)中的杰出表現(xiàn)[1,4]，引起國(guó)內(nèi)外學(xué)者的高度關(guān)注，將巖石中侵徹效應(yīng)的研究向前推進(jìn)了一大步。從研究手段上看，侵徹效應(yīng)的研究主要有試驗(yàn)、數(shù)值模擬計(jì)算及理論分析等手段。早期的研究主要是通過對(duì)大量實(shí)驗(yàn)數(shù)據(jù)的總結(jié)，歸納出經(jīng)驗(yàn)或半經(jīng)驗(yàn)的計(jì)算公式，具有簡(jiǎn)單直接的特點(diǎn)，在世界各國(guó)有關(guān)防護(hù)結(jié)構(gòu)的設(shè)計(jì)規(guī)范或計(jì)算手冊(cè)中得到廣泛應(yīng)用[8-11]，其缺點(diǎn)在于適用范圍狹窄，這些公式在實(shí)驗(yàn)特定的范圍之內(nèi)具有較高的準(zhǔn)確性，但是介質(zhì)的特征參量和彈體的侵徹速度、幾何尺度等一旦超出實(shí)驗(yàn)范圍，就會(huì)造成預(yù)測(cè)預(yù)報(bào)的偏差[1,11-12]。近三四十年來，利用大型仿真程序進(jìn)行侵徹?cái)?shù)值計(jì)算的研究獲得快速發(fā)展[5,12-15]，數(shù)值模擬的誘人之處在于能提供彈體和靶體內(nèi)部響應(yīng)的詳細(xì)信息，這些信息通常是在實(shí)驗(yàn)中也難以觀察到的，然而其計(jì)算結(jié)果往往受本構(gòu)模型與相關(guān)參數(shù)選取的限制，具有較強(qiáng)的人為性。侵徹問題本質(zhì)上是一種復(fù)雜的沖擊局部問題，主要取決于撞擊物（彈體）與被撞擊物（靶體）的撞擊速度、撞擊姿態(tài)、材料性能和彈靶雙方的幾何與結(jié)構(gòu)特性[16-18]。為揭示問題實(shí)質(zhì)，根據(jù)彈靶的幾何和物理特性，以問題的主要影響因素為切入點(diǎn)，建立具有明晰的物理力學(xué)意義的簡(jiǎn)化理論模型仍然具有十分突出的理論意義和工程價(jià)值。

傳統(tǒng)鉆地武器普遍采用高硬度合金且對(duì)地打擊速度普遍不超過1 km/s，因此經(jīng)典巖石侵徹理論中常常將彈體假設(shè)為剛性彈，理論分析的核心在于彈體阻力函數(shù)，剛性侵徹是侵徹理論研究比較成熟的部分。目前，隨著科技的發(fā)展，彈體的發(fā)射速度更高，試驗(yàn)證實(shí)[19-30]，當(dāng)彈體速度由一般侵徹速度（＜900 m/s）向高速侵徹（1 200～1 700 m/s）轉(zhuǎn)移時(shí)，在侵徹過程中會(huì)發(fā)生彈體的質(zhì)量侵蝕，嚴(yán)重影響了彈體的侵徹性能，甚至出現(xiàn)侵徹深度的“逆減”效應(yīng)。近年來高超聲速技術(shù)已經(jīng)從概念和原理探索階段進(jìn)入了實(shí)質(zhì)性的技術(shù)開發(fā)階段，外軍正在研發(fā)的超高速動(dòng)能武器對(duì)地打擊速度3 000～5 000 m/s，此時(shí)的侵徹問題將具有顯著的超高速碰撞問題的特點(diǎn)，即侵徹過程是一個(gè)由動(dòng)能急劇釋放引起的極端高溫高壓過程，不僅會(huì)發(fā)生侵徹，還會(huì)引起顯著的成坑效應(yīng)并觸發(fā)強(qiáng)烈的地沖擊效應(yīng)，其作用范圍可遠(yuǎn)超侵徹深度，其破壞效應(yīng)接近于淺埋爆炸[13,15,31-34]，其侵徹機(jī)理與破壞模式與常規(guī)侵徹問題存在顯著差別。當(dāng)前對(duì)高速和超高速?gòu)楏w對(duì)巖石的侵徹效應(yīng)及其內(nèi)在機(jī)理研究尚不十分充分，這也是當(dāng)前理論研究的熱點(diǎn)問題。

本文中圍繞大速度范圍內(nèi)巖石類介質(zhì)的侵徹效應(yīng)經(jīng)驗(yàn)公式、侵徹理論模型、成坑效應(yīng)和地沖擊效應(yīng)的研究成果，綜述國(guó)內(nèi)外學(xué)者的相關(guān)理論研究成果，并對(duì)目前研究中尚有待解決的問題和下一步的研究方向進(jìn)行展望。需要說明的是，為達(dá)到最佳的侵徹效果，鉆地武器一般具有細(xì)長(zhǎng)的殼體外殼，除非特別指出，本文中主要論述理想垂直入射條件下長(zhǎng)桿彈對(duì)巖石厚靶的一維侵徹問題，不考慮貫穿問題和斜侵徹問題，同時(shí)，為了便于介紹理論模型的發(fā)展歷程，也會(huì)涉及到部分金屬、混凝土介質(zhì)的侵徹理論。

1 侵徹效應(yīng)經(jīng)驗(yàn)公式

鉆地武器侵徹巖石深度的研究是防護(hù)工程領(lǐng)域關(guān)注的問題之一，由于沖擊侵徹現(xiàn)象的復(fù)雜性，在很長(zhǎng)時(shí)期內(nèi)，該問題的實(shí)用計(jì)算都是采用以彈道試驗(yàn)為基礎(chǔ)的經(jīng)驗(yàn)公式。建立經(jīng)驗(yàn)公式的方法大體有兩種：一種是直接對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行回歸分析建立的經(jīng)驗(yàn)公式，又稱純經(jīng)驗(yàn)公式；另一種是預(yù)先假定作用在彈體上的阻力變化規(guī)律，運(yùn)用運(yùn)動(dòng)方程推導(dǎo)出計(jì)算公式，再利用試驗(yàn)數(shù)據(jù)資料修正公式中的常數(shù)，這種方法建立的公式稱為半經(jīng)驗(yàn)半理論公式。在防護(hù)工程領(lǐng)域，當(dāng)前應(yīng)用較多的巖石侵徹經(jīng)驗(yàn)公式主要包括美國(guó)圣迪亞國(guó)家實(shí)驗(yàn)室（SNL）的Young 公式[8]、美國(guó)陸軍水道實(shí)驗(yàn)站（WES）的Bernard 公式[9]（包含三個(gè)版本）和蘇聯(lián)的別列贊公式[10]等。

（1）Young 公式形式為：

式中：所有參數(shù)均采用英制單位，h 為侵徹深度，S 為可侵徹性指標(biāo)， S = f(fcQ)與巖石無側(cè)限抗壓強(qiáng)度fc和 巖石質(zhì)量 Q （巖石完整性程度）有關(guān)，Ns1為彈頭形狀系數(shù)，mp為彈體質(zhì)量，A 為彈體截面積， v0為彈體初始撞擊速度。

（2）Bernard 公式的三個(gè)版本的主體分別為：

（3）別列贊公式的形式為：

式中：λ1為彈形系數(shù)，λ2為彈徑系數(shù)， d為彈體直徑，α 為命中角（即彈體軸線與目標(biāo)表面法線方向夾角），Kα為彈體偏轉(zhuǎn)系數(shù)，Kq為靶體材料的侵徹系數(shù)（主要與材料無側(cè)向抗壓強(qiáng)度有關(guān)）。

從參數(shù)選擇看，以上公式均以一定的方式考慮了彈體質(zhì)量、撞擊速度、彈體截面積（或彈體直徑）和巖石無側(cè)限抗壓強(qiáng)度，但不同之處在于：

（1） Young 公式和Bernard 公式均考慮了巖石的完整程度的影響，但別列贊公式未考慮，這導(dǎo)致后者無法合理預(yù)測(cè)抗壓強(qiáng)度相同但風(fēng)化程度不同的巖石的侵徹深度。

（2） Bernard 公式（III）相對(duì)于其早先版本最重要的修正之一是增加了對(duì)彈頭形狀的考慮。

（3） 別列贊公式的主要特色在于考慮了非垂直入射和彈體偏轉(zhuǎn)的影響，這是其他公式所未考慮的；但是別列贊公式對(duì)靶體性質(zhì)考慮較為簡(jiǎn)略，甚至沒有考慮巖石密度的影響。

從預(yù)測(cè)效果看[35]，Young 公式和Bernard 公式比較接近，別列贊公式的預(yù)測(cè)效果較以上公式均明顯偏大，其主要原因是別列贊公式的實(shí)驗(yàn)數(shù)據(jù)來源于對(duì)地表風(fēng)化巖石的原位侵徹實(shí)驗(yàn)，因此在相同條件下其侵徹深度較深，從防護(hù)工程設(shè)計(jì)角度而言偏于安全。

經(jīng)驗(yàn)公式能反映問題的主要影響因素，使用簡(jiǎn)便、計(jì)算可靠度較高，因而目前各國(guó)有關(guān)防護(hù)結(jié)構(gòu)的設(shè)計(jì)規(guī)范中，基本還是采用經(jīng)驗(yàn)公式計(jì)算彈體的侵徹深度，其缺點(diǎn)在于試驗(yàn)費(fèi)用高、適用范圍窄，公式只能反映彈體初始情況與最終結(jié)果間的關(guān)系，不能說明彈靶相互作用機(jī)理，不能準(zhǔn)確分析侵徹過程中彈體所受到的阻力。

此外，需要指出的是，以上公式均默認(rèn)相同條件下侵徹深度隨著彈體撞擊速度的增加而增加，但實(shí)際情況是，當(dāng)彈體打擊速度足夠高時(shí)，存在一個(gè)侵徹深度隨打擊速度增加而急劇下降的階段，與之伴隨的是彈道偏轉(zhuǎn)和彈體的侵蝕破壞，當(dāng)打擊速度進(jìn)一步增加時(shí)，由于彈靶接觸面的劇烈壓縮作用使介質(zhì)材料破壞、融化，侵徹深度隨著打擊速度的增加而緩慢增加并有趨于流體動(dòng)力學(xué)極限的趨勢(shì)[19-20,31, 36-37]，在長(zhǎng)桿彈侵徹鋁靶、混凝土靶的中也發(fā)現(xiàn)了類似的轉(zhuǎn)變現(xiàn)象[21-29,38-41]?？梢?，式(1)～(3)的適用條件主要局限于彈速不高、彈道偏轉(zhuǎn)不大、彈體變形不顯著的情形，不適用于超高速侵徹問題。

2 不同打擊速度下侵徹效應(yīng)的理論模型

2.1 對(duì)速度分區(qū)的討論

從桿形彈體對(duì)多種密實(shí)介質(zhì)（巖石、混凝土和鋁等）的侵徹實(shí)驗(yàn)結(jié)果看[19-26,36-41]，在近似理想垂直入射條件下，隨著撞擊速度增加，彈體均將經(jīng)歷從變形可忽略的階段（僅有少量質(zhì)量損失）向侵蝕階段（即彈體長(zhǎng)度嚴(yán)重縮短）的轉(zhuǎn)變，侵徹深度也將經(jīng)歷“迅速增加-逆減-緩慢增加-趨于流體動(dòng)力學(xué)極限”的過程，不同過程的物理機(jī)理不同。由于撞擊速度是決定侵徹效應(yīng)的首要因素，在速度分區(qū)問題上許多學(xué)者提出了許多建議。

Forrestal 等[41]通過鋼桿彈在0.5～3.0 km/s 范圍內(nèi)侵徹鋁靶的試驗(yàn)發(fā)現(xiàn)，隨撞擊速度提高，彈體經(jīng)歷了從剛體彈侵徹到侵蝕彈侵徹的轉(zhuǎn)變，伴隨侵深大幅下降的現(xiàn)象，Chen 等[39]根據(jù)對(duì)前述侵徹速度下彈體失效機(jī)理進(jìn)行分析，提出了“剛性彈侵徹-半流體侵徹-流體侵徹”三階段劃分模型，并利用撞擊函數(shù)I 給出了剛性彈侵徹上限和半流體侵徹下限 Ic2的臨界閾值（見圖1）。由侵蝕彈向半流體侵徹轉(zhuǎn)換的臨界值 Ic為：

圖 1 鋼桿彈撞擊鋁靶機(jī)理分區(qū)[39]Fig. 1 Schematic of penetration of steel projectiles into aluminum targets[39]

以上主要通過彈體的響應(yīng)對(duì)超高速侵徹問題進(jìn)行分區(qū)，更加完善的速度分區(qū)需要綜合考慮靶體的響應(yīng)。巖石介質(zhì)是由尺寸、形狀和礦物成分各不相同的許多顆粒連接在一起而組成的集合體，顆粒的排列方式不規(guī)則并且存在有缺陷。巖石構(gòu)造缺陷水平包括微觀上的原子層次（包括空缺和位錯(cuò)），細(xì)觀上的構(gòu)造缺陷（包括顆粒內(nèi)裂縫、沿過多個(gè)顆粒的粒間裂縫和沿粒界的裂縫、構(gòu)造上邊界碎片的分層、夾雜物等），甚至宏觀上的裂隙、節(jié)理以及層理等，其力學(xué)性質(zhì)受材料性質(zhì)、材料微結(jié)構(gòu)、荷載性質(zhì)、加載速率及材料應(yīng)力歷史和加載路徑等各種因素的制約，傳統(tǒng)的彈塑性力學(xué)本構(gòu)難以準(zhǔn)確描述。Shemyakin[47-49]研究指出在強(qiáng)爆炸沖擊作用近區(qū)巖石行為由彈性狀態(tài)向塑性狀態(tài)轉(zhuǎn)變時(shí)，不是轉(zhuǎn)向理想的塑性狀態(tài)，而是轉(zhuǎn)向顯著增強(qiáng)的塑性狀態(tài)，這種增強(qiáng)的本質(zhì)在于受限內(nèi)摩擦。對(duì)巖石介質(zhì)侵徹速度的準(zhǔn)確劃分需要準(zhǔn)確了解沖擊壓縮狀態(tài)下巖石介質(zhì)的動(dòng)力學(xué)行為及可能的力學(xué)狀態(tài)。Ahrens 等[50]基于沖擊波響應(yīng)過程的研究，認(rèn)為壓縮作用下巖石介質(zhì)的動(dòng)力響應(yīng)過程可劃分為五個(gè)階段（圖2）：0，彈性階段；1，塑性狀態(tài)；2，低壓力狀態(tài)；3，混合狀態(tài)，4，高壓力狀態(tài)。在彈性、塑性、低壓力階段，材料的強(qiáng)度起主要作用，表現(xiàn)為固體屬性；在高壓力狀態(tài)，材料的體積壓縮起主要作用，表現(xiàn)為流體動(dòng)力學(xué)屬性，可以認(rèn)為達(dá)到超高速侵徹機(jī)制的臨界條件是體積壓縮性在阻力項(xiàng)中占主導(dǎo)地位。

圖 2 巖石的動(dòng)力壓縮曲線[50]Fig. 2 Dynamic compression curves of rock[50]

王明洋等[32]在系統(tǒng)研究巖石介質(zhì)侵徹近區(qū)動(dòng)態(tài)可壓縮行為的基礎(chǔ)上，提出了流體彈塑性內(nèi)摩擦侵徹模型，認(rèn)為隨著打擊速度的提高，彈靶近區(qū)相互作用壓力增大，靶體的力學(xué)狀態(tài)發(fā)生從低應(yīng)力固體狀態(tài)、內(nèi)摩擦擬流體狀態(tài)至高應(yīng)力流體狀態(tài)的改變，在不同狀態(tài)區(qū)巖石的固體、流體屬性分配份額不同，導(dǎo)致巖體的應(yīng)力波狀態(tài)、應(yīng)力波衰減軌跡、應(yīng)力應(yīng)變狀態(tài)等力學(xué)行為發(fā)生顯著變化，具體如表1[51]所示。

表 1 沖擊壓縮作用下典型硬巖分區(qū)行為特征[51]Table 1 Dynamic behaviors of hard rock in different ranges under shock compression[51]

表1 數(shù)據(jù)與圖2 基本一致，表1 中低應(yīng)固體狀態(tài)對(duì)應(yīng)圖2 階段2，內(nèi)摩擦擬流體狀態(tài)對(duì)應(yīng)圖2 階段3，高應(yīng)力流體狀態(tài)對(duì)應(yīng)圖2 階段4。隨后，王明洋等[32]根據(jù)不同彈體速度打擊下靶體體積應(yīng)變曲線p/Ht偏離線性、剪切強(qiáng)度 τ/Ht趨近于 τp/Ht（ 見圖3）以及側(cè)壓力系數(shù) α趨近于1（見圖4）的程度，將撞擊速度大致分為三個(gè)區(qū)間：

圖 3 侵徹速度界定及介質(zhì)壓縮狀態(tài)Fig. 3 Definition the scope of penetration velocity and medium compression state

圖 4 v /c與α 關(guān)系曲線Fig. 4 Curve between v /c and α

高速或超高速打擊條件下，彈體速度 vp與 彈靶接觸面處靶體巖石質(zhì)點(diǎn)運(yùn)動(dòng)速度 v 有以下?lián)Q算關(guān)系[32]：

式中：ρp和ρt分別為彈體和靶體密度，這樣，結(jié)合公式(5)～(6)就可以判斷不同的彈體打擊速度下靶體介質(zhì)所處的力學(xué)狀態(tài)。

由彈體變形和侵蝕所導(dǎo)致的侵深逆減現(xiàn)象可能出現(xiàn)在靶體的低應(yīng)力彈塑性區(qū)和內(nèi)摩擦擬流體區(qū)，具體視彈體形狀、彈靶材料形狀等因素制約。根據(jù)高強(qiáng)合金鋼對(duì)花崗巖等典型硬巖的侵徹試驗(yàn)結(jié)果[19]（見圖5），一般彈體在靶體的低應(yīng)力固體區(qū)即發(fā)生侵蝕，這樣，典型硬巖可分為：剛性彈固體侵徹（圖5 中區(qū)域I）、侵蝕彈固體侵徹(區(qū)域II)、擬流體侵徹(區(qū)域III)、流體動(dòng)力學(xué)侵徹四個(gè)階段，與Chen 等[39]分區(qū)不同，此處綜合考慮了彈體和靶體的力學(xué)狀態(tài)，剛性彈侵徹、侵蝕彈侵徹主要針對(duì)彈體狀態(tài)，固體侵徹、擬流體侵徹、流體動(dòng)力學(xué)侵徹則主要針對(duì)靶體力學(xué)狀態(tài)。

圖 5 花崗巖侵徹實(shí)驗(yàn)結(jié)果（撞擊速度1～2.5 km/s）[19]Fig. 5 Experimental results of penetration into granite targets with impact velocity from 1?2.5 km/s[19]

2.2 剛性彈體侵徹理論

剛性彈體侵徹理論的成果最豐富、發(fā)展最成熟，主要包括：空腔膨脹理論、微分面力法、速度勢(shì)理論、滑移線法等[1,52-54]。其中，空腔膨脹理論的固有特性決定了其不同彈頭形狀、不同靶體材料均有較好的適應(yīng)性，近十多年來，在Forrestal 等[55-56]的研究成果的基礎(chǔ)上，研究人員通過對(duì)材料本構(gòu)關(guān)系的修正，即將材料體積變形方程從不可壓縮向躍變壓縮、線性壓縮、非線性壓縮發(fā)展，強(qiáng)度模型則從理想彈塑性向線性硬化、冪次硬化、壓力相關(guān)強(qiáng)度準(zhǔn)則、分階段強(qiáng)度準(zhǔn)則發(fā)展并可考慮脆性斷裂效應(yīng)，使得空腔膨脹理論得到了新的推廣[57-65]。

空腔膨脹理論應(yīng)用于彈體侵徹分析中的關(guān)鍵是根據(jù)質(zhì)量和動(dòng)量守恒方程建立空腔膨脹速度與徑向應(yīng)力的關(guān)系，進(jìn)而計(jì)算彈體受到的侵徹阻力，為了便于計(jì)算，通常將侵徹阻力寫成如下顯式形式：

式中：v是瞬時(shí)侵徹速度，F(xiàn)z是侵徹阻力，α、β、γ 分別為速度二次項(xiàng)、一次項(xiàng)和零次項(xiàng)系數(shù)。特別地，當(dāng)β=0 時(shí)，式（7）被稱為Poncelet 型阻力函數(shù)并得到了最廣泛的應(yīng)用，其經(jīng)典形式為Forrestal 等[65]提出的：

式中：Rt是材料的固有阻力，N*是彈頭形狀系數(shù)。對(duì)于混凝土，Rt近似滿足如下經(jīng)驗(yàn)關(guān)系：

式中：fc為混凝土單軸抗壓強(qiáng)度。

在滿足Poncelet 型阻力函數(shù)時(shí)，由空腔膨脹理論得到的典型侵徹深度表達(dá)式為：

式中：K1和K2為靶體材料密度、強(qiáng)度和彈頭形狀決定的參數(shù)。值得注意的是，F(xiàn)rew 等[20]、張德志等[36]分別采用式（8a）對(duì)石灰?guī)r和花崗巖侵徹實(shí)驗(yàn)結(jié)果進(jìn)行擬合，發(fā)現(xiàn)Rt與彈體撞擊速度和彈體尺寸均密切相關(guān)，這說明式（8a）并不完全適合巖石介質(zhì)。錢七虎等[51]、王明洋等[66]采用應(yīng)力短波條件下的內(nèi)摩擦侵徹理論得到了彈體阻力的表達(dá)式，并根據(jù)參數(shù)分析發(fā)現(xiàn)恰恰是Poncelet 型阻力函數(shù)所忽略的速度線性阻力函數(shù)在硬巖介質(zhì)的中高速侵徹過程中起決定性作用，因此侵徹深度與撞擊速度之間近似滿足線性關(guān)系，這與經(jīng)驗(yàn)公式較為一致。

2.3 變形彈侵徹計(jì)算理論

實(shí)驗(yàn)研究表明，高速侵徹條件下，彈體可以在未出現(xiàn)侵蝕和顯著質(zhì)量損失的情況下發(fā)生磨蝕、鈍化、彎折等情況，嚴(yán)重影響彈道的穩(wěn)定性，引起侵徹深度陡然降低。上述情形被歸結(jié)為變形彈侵徹問題，此時(shí)仍然可以采用牛頓第二運(yùn)動(dòng)定律描述彈體運(yùn)動(dòng)，但必須計(jì)及彈頭的磨蝕和質(zhì)量的損失。

在彈體變形機(jī)理方面，何麗靈等[27,67]、Guo 等[68]研究了侵徹混凝土?xí)r彈體磨蝕和質(zhì)量損失的物理過程。何麗靈等[69-70]引入修正系數(shù)表征熔化彈體表面材料脫落和靶中骨料硬度對(duì)彈體質(zhì)量損失的影響，給出了質(zhì)量增量表達(dá)式，計(jì)算時(shí)，每個(gè)時(shí)步內(nèi)認(rèn)為彈頭形狀保持不變，計(jì)算完成后，再根據(jù)彈頭質(zhì)量增量確定彈頭形狀，計(jì)算認(rèn)為彈頭質(zhì)量損失對(duì)侵徹影響可以忽略，主要是彈頭形狀影響侵徹深度。鄭浩[71]采用彈靶分離的方法和修正后Archard 磨損彈體侵蝕模型，建立了彈體頭部鈍化演變差分計(jì)算方法，得到了不同撞擊速度下彈體頭部形狀。楊陽[72]利用量綱分析推導(dǎo)了考慮骨料硬度的彈體侵徹混凝土靶體的質(zhì)量損失無量綱公式。

在侵徹模型方面，Zhao 等[28]指出隨著初始撞擊速度的增加，彈頭將被磨損，并且明顯地從截卵形彈頭變?yōu)榘肭蛐危踔粮g的形狀，據(jù)此根據(jù)試驗(yàn)結(jié)果建立了彈頭形狀系數(shù)與瞬時(shí)侵徹速度的關(guān)系和彈體侵徹過程中的質(zhì)量，最后采用剛性彈侵徹計(jì)算方法得到變形彈侵徹計(jì)算模型，并發(fā)現(xiàn)在考慮磨損效應(yīng)后侵徹深度將在某一撞擊速度下出現(xiàn)極值，大致對(duì)應(yīng)于實(shí)驗(yàn)中發(fā)現(xiàn)的侵徹深度突降臨界速度。

Lan 等[73]假設(shè)撞擊的瞬間，彈體頭部立即發(fā)生變化（見圖6），并且橫截面積保持Ad不再變化，且變化后的頭部仍然保持球形，整個(gè)已變形區(qū)域的質(zhì)點(diǎn)速度均為 vm，在界面處存在質(zhì)點(diǎn)速度和橫截面積間斷，且壓力在界面處均勻分布，開坑截面積保持不變，得到了變形彈侵徹時(shí)彈尾速度與侵徹速度的關(guān)系：

式中：au由試驗(yàn)結(jié)果確定， v0為 初始撞擊速度， vg為從剛性彈侵徹進(jìn)入變形彈侵徹的臨界速度。最終得到了變形彈階段的侵徹深度計(jì)算公式：

式中：L 為彈體長(zhǎng)度。從以上研究成果總體來看，目前彈體變形過程中質(zhì)量磨蝕的影響還未定論，彈頭鈍粗與質(zhì)量磨蝕對(duì)侵徹的影響程度還需明晰，彈頭鈍粗的程度還很大程度上依賴于試驗(yàn)擬合結(jié)果。彈道穩(wěn)定性研究方面，還未見考慮彈體磨蝕引起的彈道偏移的定量理論研究成果。作者及其同事在花崗巖侵徹實(shí)驗(yàn)中未觀察到彈體墩粗的現(xiàn)象(見圖7)[19]，據(jù)此認(rèn)為侵徹深度的減小主要是彈體質(zhì)量損失引起的（圖5 中區(qū)域II 彈體質(zhì)量急劇降低伴隨著侵徹速度的急速逆轉(zhuǎn)），通過引入彈體質(zhì)量損失函數(shù)建立了變形彈侵徹理論[31]，這部分介紹詳見2.5.2 節(jié)。

圖 6 變形彈計(jì)算假設(shè)[73]Fig. 6 Assumption of deformable projectile[73]

圖 7 不同速度侵徹后回收彈體的形態(tài)[19]Fig. 7 The morphology of the recovered projectiles after penetration under different velocities[19]

2.4 半流體和流體侵徹理論模型

典型的超高速侵徹過程可以根據(jù)彈靶界面壓力隨時(shí)間的演化過程劃分為四個(gè)階段[74]（見圖8），即階段Ⅰ：瞬態(tài)激波階段，階段Ⅱ：穩(wěn)態(tài)侵徹階段（或主要侵徹階段），階段Ⅲ：二次侵徹階段（或慣性流動(dòng)階段），階段Ⅳ：彈性恢復(fù)階段。對(duì)于桿形彈而言，階段Ⅱ?qū)η謴剡^程起決定性作用，此時(shí)可以采用經(jīng)過修正的穩(wěn)態(tài)流體動(dòng)力學(xué)模型來分析機(jī)理；對(duì)于球形彈而言，階段Ⅱ的影響很小，一般需要綜合階段Ⅰ和階段Ⅲ進(jìn)行分析。本節(jié)重點(diǎn)綜述桿形彈的超高速侵徹理論模型，球形彈的超高速侵徹理論模型將在本文3.1 節(jié)“成坑效應(yīng)”中闡述。

圖 8 基于界面壓力的超高速侵徹階段劃分[74]Fig. 8 Phases during hypervelocity penetration based on interface pressure[74]

2.4.1 流體動(dòng)力學(xué)模型

桿形彈體超高速侵徹的理論模型最早來自高速射流的流體動(dòng)力學(xué)理論，當(dāng)彈體的速度極高，彈靶接觸面壓力極大導(dǎo)致可忽略彈體和靶體的強(qiáng)度時(shí)，可將桿形彈侵徹簡(jiǎn)化為聚能射流問題，彈靶接觸面的壓力平衡關(guān)系可由Birkhoff 等[75]建議的伯努利方程描述：

式中： vp為射流速度（彈體尾部速度）， u為彈靶界面速度，ρp為射流（彈體）密度，ρt為靶體密度。假設(shè)侵徹為定常過程，可以由公式(13)得到無量綱的侵徹深度：

式中：L 為彈體長(zhǎng)度。式（14）被稱為桿形彈的“流體動(dòng)力學(xué)侵徹極限”。需要指出的是，盡管上述計(jì)算結(jié)果經(jīng)常被視作連續(xù)射流和長(zhǎng)桿彈體在速度趨于無窮時(shí)的侵徹深度理論極限值，但在現(xiàn)有實(shí)驗(yàn)速度條件下的實(shí)際侵深往往與之存在顯著偏差，因此必須對(duì)上述模型進(jìn)行修正以使之更加符合實(shí)際。

2.4.2 修正流體動(dòng)力學(xué)模型

侵蝕彈侵徹時(shí)，只是彈靶接觸部分呈流體狀態(tài)，其余部分還處于剛體狀態(tài)[76]，因此實(shí)際上不能忽略彈體和靶體的強(qiáng)度特征。流體動(dòng)力學(xué)理論的最大不足在于未考慮材料強(qiáng)度，因而只適用于速度極高情況下的侵徹行為，對(duì)于半流體侵徹階段，彈靶相互作用的描述一般采用修正的流體動(dòng)力學(xué)模型。經(jīng)典的修正流體動(dòng)力學(xué)模型都是首先基于金屬靶體侵徹提出的，包括Allen-Rogers 模型[77]（A-R 模型）、Alekseevskii-Tate 模型[42-44]（A-T 模型）等。

（1）Allen-Rogers 模型

在聚能射流理論的流體動(dòng)力學(xué)模型基礎(chǔ)上，1961 年，Allen 和Rogers 在伯努利方程中加入強(qiáng)度項(xiàng)考慮靶體強(qiáng)度效應(yīng)的影響，形成Allen-Rogers 模型[77]：

式中： σ于靶體材料強(qiáng)度相關(guān)，將該式對(duì)時(shí)間積分可得無量綱侵徹深度表達(dá)式：

當(dāng) vp→∞（ 或忽略 σ）時(shí)，式(16)可退化成式(14)，Allen 和Rogers 成功用該模型解釋了鎂、鋁、錫等桿形彈高速撞擊鋁靶的試驗(yàn)數(shù)據(jù)，在高速侵徹作用下，侵徹深度趨近于流體動(dòng)力學(xué)極限。

（2）Alekseevskii-Tate 模型

A-T 模型是最經(jīng)典的桿形彈高速侵徹的理論模型，該模型由Alekseevskii[42]和Tate[43-44]各自提出，他們將彈體和靶體材料的強(qiáng)度（Yp和Rt）引入伯努利方程中，并聯(lián)合彈長(zhǎng)變化方程、侵徹方程和彈體減速運(yùn)動(dòng)方程，建立了侵徹計(jì)算的流體力學(xué)模型：

該模型假設(shè)彈體侵徹過程中僅彈體頭部較小區(qū)域和彈靶接觸面附近靶體處于流體狀態(tài)，其余彈體部分仍為剛體。對(duì)應(yīng)不同的彈靶組合，有兩種不同侵徹情形：當(dāng)Yp＜Rt時(shí)，彈體邊侵徹邊侵蝕，直到彈體速度 vp下 降到臨界速度時(shí)侵徹停止；當(dāng)Yp＞Rt時(shí)，彈體速度 vp下降到臨界速度，剩余彈體以剛性彈繼續(xù)侵徹。

對(duì)于A-T 模型求解，一般采用數(shù)值方法，由于Yp和Rt會(huì)對(duì)計(jì)算結(jié)果有較大影響，因此模型中Yp和Rt取值和理解一直是A-T 模型分析的重點(diǎn)和難點(diǎn)，目前仍難以達(dá)成一致性觀點(diǎn)[40]。Tate 最初曾建議將Yp取為Hugoniot 彈性極限（即HEL）[42]，而Rt取靶體材料HEL 的3.5 倍。后來在1986 年，通過實(shí)驗(yàn)數(shù)據(jù)擬合，重新評(píng)估了Yp和Rt[44]：

式中：σyp為彈體的動(dòng)態(tài)屈服應(yīng)力，σyt和Et分別為靶體的動(dòng)態(tài)屈服強(qiáng)度和彈性模量， δ一般取0.7。

（3）其他修正流體動(dòng)力學(xué)模型

在A-T 模型基礎(chǔ)上，國(guó)內(nèi)外學(xué)者提出了諸多改進(jìn)模型，如：孫庚辰-吳錦云-趙國(guó)志-史駿模型[78]（S-WZ-S 模型）、Rosenber-Marmor-Mayseless 模型[79]（R-M-M 模型）、Walker-Anderson 模型[80]（A-W 模型）、Zhang-Huang 模型[81]（Z-H 模型）、Lan-Wen 模型[73]（L-W 模型）等。這些模型最早都是針對(duì)金屬靶體提出的。樓建鋒[53]對(duì)各個(gè)模型的基本假設(shè)、控制方程和計(jì)算結(jié)果作了全面的比較，并認(rèn)為上述模型的關(guān)鍵控制方程均可以統(tǒng)一描述為如下形式，即：

式中： vp和u 分別表示彈體尾部速度和彈-靶界面速度，ρp為彈體密度，[Yp]和[Rt]分別為彈體名義強(qiáng)度和靶體名義阻力。盡管式（22）在形式上實(shí)現(xiàn)了模型的統(tǒng)一，但實(shí)際上各模型的基本假設(shè)、參數(shù)取值和預(yù)測(cè)效果差異很大，在文獻(xiàn)[1]中給出了以上各模型中[Yp]和[Rt]的比較（表2）。

表 2 不同修正流體動(dòng)力學(xué)模型中[Yp]和[Rt]值[1]Table 2 The value of [Yp] and [Rt] in different models[1]

從中可以看出，對(duì)于以A-T 模型為代表的修正流體動(dòng)力學(xué)模型，研究的矛盾和難點(diǎn)主要集中在模型中[Yp]和[Rt]的取值。

準(zhǔn)確獲知侵徹過程中的[Yp]和[Rt]值是一個(gè)十分艱難的任務(wù)。對(duì)于[Yp]，Rosenberg 等[82]對(duì)不同強(qiáng)度長(zhǎng)桿彈侵徹的模擬計(jì)算表明，[Yp]與彈靶強(qiáng)度、撞擊速度和長(zhǎng)徑比都相關(guān)，因此認(rèn)為[Yp]是A-T 模型中不能準(zhǔn)確定義的參數(shù)。由于[Yp]控制彈體侵蝕和減速，Anderson 等[83]曾建議通過實(shí)驗(yàn)測(cè)量彈尾的實(shí)時(shí)運(yùn)動(dòng)來推測(cè)[Yp]值，但目前為止，仍未見精確的侵徹過程中[Yp]實(shí)測(cè)數(shù)據(jù)。

由修正流體動(dòng)力學(xué)模型的控制方程（公式(22)）可知，影響彈體侵徹效應(yīng)的是模型中[Rt]與[Yp]的差值[Rt]- [Yp]，目前研究人員普遍認(rèn)為其對(duì)長(zhǎng)桿彈高速侵徹能力的影響較小，因此通常取[Yp]為定值或零強(qiáng)度進(jìn)而主要研究[Rt]的規(guī)律。目前常采用三種方式確定[Rt]值[84]：(1) 通過空腔膨脹等理論模型進(jìn)行推導(dǎo)；（2）通過侵徹實(shí)驗(yàn)數(shù)據(jù)反向擬合；（3）在數(shù)值模擬中獲得瞬時(shí)壓力，再對(duì)時(shí)間或位置積分獲得侵徹中的平均值。由于采用的模型不同，不同的修正理論獲得的[Rt]通常具有顯著的差異[1,84]，由于A-T 及其修正模型主要針對(duì)穩(wěn)態(tài)侵徹階段，同時(shí)[Rt]在侵徹過程中劇烈變化，因此利用最終侵徹深度反向擬合[Rt]的方法也不盡合理。Anderson 等[83]發(fā)現(xiàn)在A-T 模型中，無法同時(shí)匹配侵徹速度和侵徹深度，隨后Anderson 等[85]詳細(xì)比較了利用侵徹深度反向擬合的靶體阻力和對(duì)應(yīng)數(shù)值模擬中按時(shí)間平均、按侵深平均及僅考慮穩(wěn)態(tài)侵徹階段的靶體阻力，發(fā)現(xiàn)二者差異顯著，在超過4.5 km/s 的撞擊速度下，用侵徹深度反向擬合的[Rt]為負(fù)值，這顯然違背了客觀物理規(guī)律。在巖石類介質(zhì)侵徹效應(yīng)的理論研究，如何準(zhǔn)確確定[Rt]和[Yp]值，仍將是一個(gè)艱難的工作。

2.5 統(tǒng)一的內(nèi)摩擦流體彈塑性侵徹理論模型

目前關(guān)于侵徹理論計(jì)算模型主要分空腔膨脹理論及射流理論兩種，其中空腔膨脹理論主要適用于研究固體彈塑性侵徹問題；而射流理論則主要適用于研究流體動(dòng)力學(xué)侵徹問題，目前尚缺乏一種涵蓋從低速至高速、超高速侵徹的全過程理論模型。王明洋等[31-32]、李杰等[86-87]在系統(tǒng)總結(jié)爆炸和沖擊加載作用下巖石動(dòng)態(tài)壓縮試驗(yàn)數(shù)據(jù)的基礎(chǔ)上，指出：在固體彈塑性侵徹區(qū)域與流體動(dòng)力學(xué)侵徹區(qū)域之間，還存在一個(gè)擬流體過渡區(qū)，在這一區(qū)域材料的行為兼具固體和流體屬性，隨后創(chuàng)新提出流體彈塑性內(nèi)摩擦侵徹理論，表征了材料從低應(yīng)力固體彈塑性至高應(yīng)力流體之間的應(yīng)力狀態(tài)，推導(dǎo)出了從固體侵徹至流體侵徹全過程的阻抗演變公式：

式中： v?和 v??分別為巖石介質(zhì)進(jìn)入擬流體侵徹和流體動(dòng)力學(xué)侵徹對(duì)應(yīng)的彈體臨界速度，按公式(5)～(6)確定，dv為擬合系數(shù)，由試驗(yàn)確定，隨彈速增加，公式（24）中不同項(xiàng)應(yīng)力影響分配份額發(fā)生變化，根據(jù)侵徹壓力狀態(tài)遞進(jìn)過程中，不同參數(shù)演化趨向極限的程度，將侵徹過程分為固體侵徹、半流體侵徹和流體動(dòng)力學(xué)侵徹。下式給出了3種侵徹情況下，阻抗的計(jì)算公式和速度閾值：

在不同的侵徹狀態(tài)下，通過對(duì)彈頭進(jìn)行受力分析得到彈體的最終侵徹深度。

2.5.1 剛性彈侵徹階段

在剛性彈侵徹時(shí)，巖石介質(zhì)處于固體侵徹狀態(tài)，根據(jù)牛頓第二定律得到彈體運(yùn)動(dòng)微分方程[32,51,66]：

式中：F 為彈頭阻力，受彈頭形狀影響，應(yīng)通過分析心角為θ0（ θ0=arctan((s?r0)/s)=arctan(1?1/2ψ)，ψ=s/2r0為彈頭頭部曲率），彈體垂直侵入靶體，初始侵入速度為 v0，中途侵入速度為 vp，則作用在彈彈頭微面積上的阻力積分得到。例如，對(duì)于弧形彈頭（見圖9），設(shè)彈桿半徑為r0，彈頭表面某點(diǎn)的法線方向與彈軸的夾角為θ，彈頭圓弧半徑為s，圓體頭部微面積上的法向阻力和切向阻力分別為：

圖 9 弧形彈頭的幾何參數(shù)Fig. 9 Ogive-nose projectile geometry

式中： σr按公式(25)固體侵徹階段取值， μs為彈靶間摩擦因數(shù)，則彈頭軸向合力為：

王明洋等[51,66]分析了不同彈體形狀下剛性彈侵徹階段彈體侵徹深度，統(tǒng)一用下列公式表示：

對(duì)于弧形彈：

當(dāng)0.1≤v0/cp≤0.2時(shí)，式(29a)中對(duì)數(shù)項(xiàng)影響小于5%，可簡(jiǎn)化為：

2.5.2 侵蝕彈侵徹階段

隨著侵徹速度進(jìn)一步增加，當(dāng)巖體介質(zhì)阻抗超過彈體的動(dòng)力屈服強(qiáng)度時(shí)，彈頭將發(fā)生屈服、磨蝕，造成由于彈體變形和質(zhì)量損失帶來的侵徹深度的急劇下降，仍然采用式（26）所示的控制方程，但需要引入彈體質(zhì)量的損失與侵徹速度的關(guān)系[32]：

式中： mp0為彈體的初始質(zhì)量， vcr為發(fā)生質(zhì)量損失的臨界速度，由條件σr=Yp來確定，當(dāng) vp＞vcr時(shí)，彈體侵徹時(shí)出現(xiàn)質(zhì)量侵蝕，侵徹過程中，彈體速度不斷減小，被侵蝕的質(zhì)量不斷增加，當(dāng)速度降低至 vp≤vcr時(shí)，彈體重新恢復(fù)剛體侵徹，彈體質(zhì)量不再發(fā)生變化，αe為質(zhì)量損失參數(shù)，可通過實(shí)驗(yàn)確定。

此時(shí)將公式（31）代入彈體運(yùn)動(dòng)微分方程（26），可得到侵蝕彈侵徹階段彈體侵徹深度計(jì)算公式：

2.5.3 內(nèi)摩擦擬流體侵徹階段

對(duì)于堅(jiān)硬巖石，當(dāng)靶體進(jìn)入內(nèi)摩擦擬流體狀態(tài)時(shí)，彈體一般也已進(jìn)入流體狀態(tài)，此時(shí)須采用修正的流體動(dòng)力學(xué)方程描述彈體行為，即[32]：

式中： σr按公式(25)擬流體侵徹階段取值。若彈體強(qiáng)度Yp可忽略，在理想定常侵徹條件下，得到侵徹深度計(jì)算公式：

式中： κ按公式(24)函數(shù)給出。對(duì)于非定常侵徹，可取公式(33)為控制方程，并聯(lián)合A-T 模型的彈長(zhǎng)變化方程、侵徹方程和彈體減速運(yùn)動(dòng)方程進(jìn)行數(shù)值求解計(jì)算。

2.5.4 流體動(dòng)力學(xué)侵徹階段

在公式(34)中，若κ→1，進(jìn)入流體動(dòng)力學(xué)階段，侵徹深度計(jì)算公式演化為[32]：

隨著 Ma的 增大， θ →1，于是式(35)退化成為：

內(nèi)摩擦侵徹理論模型實(shí)現(xiàn)了由低速至高速、超高速侵徹的全過程計(jì)算，并得到了1.0～5.0 km/s鋼彈侵徹花崗巖的試驗(yàn)驗(yàn)證[32]（見圖10）。

圖 10 花崗巖侵徹深度的實(shí)驗(yàn)結(jié)果與理論計(jì)算結(jié)果的預(yù)測(cè)效果[32]Fig. 10 Comparison of calculation results with experimental results of penetration depth in granite[32]

3 超高速侵徹成坑和地沖擊效應(yīng)

研究表明，隨著彈體速度提高，超高速侵徹深度逐漸趨向極限，但彈體動(dòng)能急劇釋放引起的極端高溫高壓過程致使彈坑半徑呈現(xiàn)非線性擴(kuò)增現(xiàn)象，同時(shí)產(chǎn)生類似于爆炸的強(qiáng)地沖擊現(xiàn)象，其本質(zhì)是一種幅值高、作用時(shí)間長(zhǎng)的應(yīng)力波[13,33]。由成坑效應(yīng)和地沖擊效應(yīng)帶來的附加毀傷效應(yīng)值得引起工程防護(hù)設(shè)計(jì)和武器研發(fā)人員的重視，同時(shí)對(duì)研究隕石撞擊效應(yīng)和破巖技術(shù)具有重要的意義。

3.1 成坑效應(yīng)

目前成坑效應(yīng)的研究主要集中于球形彈撞擊效應(yīng)，研究背景主要來源于破片或太空垃圾對(duì)航天航空飛行器的破壞以及隕石對(duì)類地星體的撞擊成坑等。對(duì)于地質(zhì)類材料，地質(zhì)學(xué)和行星學(xué)家圍繞小行星和隕石撞擊成坑效應(yīng)，開展了大量室內(nèi)模型實(shí)驗(yàn)，并綜合采用參數(shù)擬合方法以及半解析手段得到了金屬介質(zhì)中的成坑尺寸計(jì)算公式。如Burchell 等[88]根據(jù)花崗巖的撞擊試驗(yàn)擬合獲得了經(jīng)驗(yàn)公式：

式（37）要求采用mm、km/s 作為長(zhǎng)度和速度單位；Takagi 等[89]根據(jù)玄武巖的撞擊試驗(yàn)得到了成坑參數(shù)和彈體動(dòng)能之間的關(guān)系

更早時(shí)候，Donald[90-91]基于100 多次玄武巖和花崗巖的試驗(yàn)結(jié)果歸納出了成坑深度和成坑直徑的經(jīng)驗(yàn)公式

式（39）要求采用cm、g、s 作為長(zhǎng)度、質(zhì)量和時(shí)間單位進(jìn)行計(jì)算，其中i 是入射角度，以垂直入射為i = 90°且要求i＞15°。

以上關(guān)系大多是零星的經(jīng)驗(yàn)公式，缺乏針對(duì)不同因素、不同尺度下相似關(guān)系的考慮。經(jīng)典的超高速撞擊成坑的相似關(guān)系在Kinslow[92]和張慶明等[93]的著作中得到了較為充分的總結(jié)，而專門針對(duì)地質(zhì)類材料的相似關(guān)系的系統(tǒng)理論成果則首推Holsapple、Housen 和Schmidt 等學(xué)者等在過去二十多年積累的成果[94-101]，其形式總體滿足表2 的結(jié)果，其中ρp為彈體密度，a0為彈體半徑，Yt為靶體強(qiáng)度，ε 為靶體孔隙率參數(shù)，含有不同下標(biāo)的f 為ε 的函數(shù)且與對(duì)應(yīng)的幾何參數(shù)和控制區(qū)域有關(guān)，g 為重力加速度，ξ1，ξ2為擬合常數(shù)。其中ξ2一般取為1/3，1/3＜ξ1＜2/3 且ξ1與ε 有關(guān)：當(dāng)ε = 30%～35%時(shí)，ξ1=0.4；當(dāng)介質(zhì)孔隙率較小時(shí)，ξ1=0.6。當(dāng)彈體直徑不超過～10 m 級(jí)別時(shí)，成坑相似率在強(qiáng)度控制區(qū)域；而當(dāng)彈體直徑超過～10 m 級(jí)別時(shí)，成坑相似率在重力控制區(qū)域。表3 中的公式得到了海量實(shí)驗(yàn)數(shù)據(jù)的驗(yàn)證并形成了適用于寬廣條件的計(jì)算參數(shù)，因此得到了較為廣泛的應(yīng)用并有專門的網(wǎng)站供讀者進(jìn)行估算[102]。近年來MEMIN 課題組利用上述相似關(guān)系開展了大量的室內(nèi)模型實(shí)驗(yàn)[103-108]，并從宏觀和微觀角度加深了對(duì)成坑機(jī)理的理解。

表 3 地質(zhì)類材料成坑效應(yīng)的相似關(guān)系Table 3 Similarity laws of cratering effects in geological material

目前對(duì)于桿形彈成坑效應(yīng)的研究不多。但目前隨著高超聲速動(dòng)能武器的發(fā)展，已逐漸引起世界各國(guó)學(xué)者的關(guān)注。

完整的成坑效應(yīng)包括不同侵徹深度處彈坑橫向尺寸的全部信息，一般來說，地質(zhì)材料和金屬介質(zhì)的撞擊成坑外觀差異十分顯著（見圖11）。在垂直撞擊下，金屬介質(zhì)中將形成內(nèi)壁光滑、輪廓清晰的圓形截面彈坑（圖12(c)），在彈坑邊緣可以觀察到顯著外翻的“唇沿”，附近介質(zhì)顯著隆起[93,109]。相比之下，巖石介質(zhì)中的成坑一般難以用簡(jiǎn)單形狀函數(shù)描述，其典型特征是在中央坑周圍存在一個(gè)邊緣極不規(guī)則的剝裂區(qū)域，其內(nèi)壁凹凸不平。這一不規(guī)則區(qū)域主要是由于反射應(yīng)力波引起的表面剝離效應(yīng)，一般認(rèn)為是與巖石介質(zhì)抗拉強(qiáng)度低、脆性顯著、內(nèi)部缺陷分布有關(guān)。

圖 11 巖石與金屬靶體超高速撞擊成坑的典型外觀[106, 109]Fig. 11 Typical appearance of craters formed by hypervelocity impact[106, 109]

圖 12 不同條件下的成坑效應(yīng)Fig. 12 Cratering effects under different conditions

當(dāng)撞擊速度較低時(shí)，典型巖石介質(zhì)成坑效應(yīng)包括表面淺碟形的開坑區(qū)和隧道區(qū)（圖12（a））[110-111]；當(dāng)撞擊速度較高時(shí)，上述兩個(gè)區(qū)域的邊界趨于模糊，隧道區(qū)逐漸演化為中央彈坑（圖12（b））[19, 21-22, 103-108, 112]。

現(xiàn)有研究主要針對(duì)開坑區(qū)范圍計(jì)算。Forrestal 等[23、Frew 等[113]對(duì)一般彈道速度下開坑區(qū)深度hc和開坑直徑Dc進(jìn)行了分析，認(rèn)為混凝土hc是彈體直徑的2 倍，Dc為彈體直徑的2～8 倍。Li 等[111]采用滑移線理論給出了混凝土hc的深度的簡(jiǎn)化計(jì)算方法，并認(rèn)為其主要取決于彈頭形狀和彈體直徑，與撞擊速度無關(guān)。劉海鵬等[110]給出了混凝土hc和Dc的解析計(jì)算方法，發(fā)現(xiàn)開坑區(qū)傾斜角度與材料的力學(xué)性能和撞擊條件關(guān)系不大并接近于24.7°。

由于長(zhǎng)桿彈成坑效應(yīng)機(jī)理復(fù)雜，對(duì)于長(zhǎng)桿彈成坑效應(yīng)的計(jì)算主要采用半理論半經(jīng)驗(yàn)的計(jì)算公式，許多學(xué)者將類似于表2 中的相似關(guān)系引入長(zhǎng)桿彈侵徹，并通過試驗(yàn)建立了成坑效應(yīng)計(jì)算公式，如錢秉文等[21]利用Holsapple-Housen 模型得到了鎢合金短桿彈超高速侵徹混凝土成坑效應(yīng)中Dc和Vc的歸一化表達(dá)式：

式中：k 為擬合參數(shù)，Yt取混凝土單軸抗壓強(qiáng)度。程怡豪等[114]得到了涵蓋文獻(xiàn)[21-22,29]實(shí)驗(yàn)數(shù)據(jù)的混凝土成坑深度的分段計(jì)算公式（見圖13）：

式中：Yt建議取為式（8b）中的Rt。

以上半理論半經(jīng)驗(yàn)公式具有簡(jiǎn)單直接的特點(diǎn)，在實(shí)驗(yàn)特定的范圍之內(nèi)具有較高的準(zhǔn)確性，其缺點(diǎn)在于介質(zhì)的特征參量和彈體的侵徹速度、幾何尺度等一旦超出實(shí)驗(yàn)范圍，就會(huì)很容易造成預(yù)測(cè)預(yù)報(bào)的偏差，因此需要從巖石的成坑機(jī)理出發(fā)，建立物理意義明確的理論模型。

圖 13 式（41）與混凝土超高速侵徹深度實(shí)驗(yàn)結(jié)果的對(duì)比[114]Fig. 13 Comparison between Eq. (41) and experimental results of hyper-velocity penetration depth into concrete[114]

圖 14 高速?gòu)楏w侵徹巖石擴(kuò)孔范圍計(jì)算簡(jiǎn)圖[116]Fig. 14 Calculation diagram of cavitation induced by high-velocity projectile penetration into rocks[116]

目前對(duì)于長(zhǎng)桿彈成坑效應(yīng)的理論模型，較為經(jīng)典的是Slepyan 模型。桿形彈超高速撞擊巖石過程中，彈靶近區(qū)巖石的超高壓力致使巖石破碎，并經(jīng)由彈靶接觸邊界向外噴射[21,115]。Slepyan[116]采用如圖14所示的計(jì)算模型研究了擴(kuò)孔范圍，圖中Ω1（實(shí)質(zhì)上是一個(gè)外直徑2R0的變截面空心圓柱）為擬流體區(qū)，該區(qū)域的破碎巖石可以視為無黏性不可壓縮理想流體；Ω2為裂紋區(qū)，該區(qū)域內(nèi)的質(zhì)點(diǎn)位移很小并可視為剛性；Ω3為彈性區(qū)，仍舊保持著巖石的初始物理力學(xué)特征。擬流體區(qū)與裂紋區(qū)的邊界滿足P=Ht，并且隨著彈體的侵徹，該邊界不斷沿著x1軸移動(dòng)。因此，上述問題轉(zhuǎn)化為擬流體破碎介質(zhì)以強(qiáng)度為邊界的管道中遇到彈體阻礙時(shí)的流動(dòng)問題。

Slepyan 模型將破碎巖石可以視為無黏性不可壓縮理想流體，實(shí)際上由于巖石力學(xué)性質(zhì)特殊性，很難達(dá)到如此力學(xué)狀態(tài)。王明洋等[31]在Slepyan 模型基礎(chǔ)上，通過內(nèi)摩擦理論（2.5 節(jié)）對(duì)模型中流體項(xiàng)的修正，得到Ω1內(nèi)的伯努利方程和連續(xù)方程為：

Ma=vp/c?0， 的計(jì)算可借助Gurevich[117]的研究結(jié)果：對(duì)于彈體侵徹這類圓錐體的軸對(duì)稱問題求解極其復(fù)雜，但是試驗(yàn)與數(shù)值分析表明，平面問題與軸對(duì)稱問題的壓縮射流系數(shù)基本相等，可以將等效楔形體所致的壓縮射流系數(shù)視為相應(yīng)圓錐體的壓縮射流系數(shù)。平面問題壓縮射流的系數(shù)由下列公式確定：

式中：α 是楔塊頂端的半角（0 ＜α≤π/2） ， x=α/π， ζx=Ma·(1+Ma2)?0.5。

通過粉碎區(qū)邊界壓力做功與裂紋增長(zhǎng)能量之間的關(guān)系，得到徑向裂紋區(qū)半徑 Rc簡(jiǎn)單的相似關(guān)系：

圖15 給出了擬流體侵徹范圍徑向裂紋區(qū)半徑理論計(jì)算結(jié)果與花崗巖中實(shí)驗(yàn)結(jié)果對(duì)比。結(jié)果表明，相較于Slepyan 模型，利用內(nèi)摩擦理論修正的計(jì)算模型撞擊成坑計(jì)算公式與實(shí)驗(yàn)結(jié)果吻合更好。

以上公式及表2 的結(jié)果均是針對(duì)超高速撞擊成坑最終形態(tài)的計(jì)算公式，不能再現(xiàn)成坑的全過程。Maxwell 的Z 模型[118]是較早成功用于描述超高速撞擊成坑和核爆成坑過程的理論模型，該模型對(duì)彈坑邊緣和前驅(qū)沖擊波陣面之間的速度場(chǎng)進(jìn)行矢量化描述，并在一定的邊界條件和初始條件下求解最終成坑尺寸和介質(zhì)拋射過程，其基本控制方程為：

圖 15 徑向裂紋區(qū)半徑計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果對(duì)比[32]Fig. 15 Comparison of crater radius between calculation results and experimental results[32]

式中： R為 徑向方向向量，αz和Z 近似為常數(shù)， U為速度場(chǎng)向量。式（46）描述了沖擊波后不可壓縮條件下流場(chǎng)的運(yùn)動(dòng)，當(dāng)Z = 2，3，4 時(shí)，沖擊波后流場(chǎng)如圖16 所示。求解Z 模型尚需材料本構(gòu)模型和邊界條件，過程較為復(fù)雜，一般必須借助于數(shù)值方法，但對(duì)于特殊情形可以求得解析解。例如文獻(xiàn)[34]的模型可以視作Z=2 時(shí)的特例，此時(shí)流場(chǎng)滿足球?qū)ΨQ條件，因此可以采用球形空腔膨脹理論進(jìn)行分析。

圖 16 Z 模型在不同Z 值下的速度場(chǎng)[118]Fig. 16 Streamlines with different values of Z[118]

綜上，超高速打擊情況下，彈靶界面處的沖擊壓力遠(yuǎn)大于材料的動(dòng)態(tài)強(qiáng)度，材料呈現(xiàn)出準(zhǔn)流體或流體狀態(tài)，在接觸面周圍流動(dòng)，最終形成遠(yuǎn)大于彈體口徑的彈坑，并在彈坑周圍形成破碎區(qū)和徑向裂紋區(qū)，巖石的成坑特性反映了彈體打擊過程中能量分配的份額，從而決定了輻射出來的應(yīng)力波基本參數(shù)，對(duì)于評(píng)判超高速動(dòng)能武器打擊的毀傷效應(yīng)具有重要價(jià)值，但是相關(guān)研究尚不完善。目前在實(shí)驗(yàn)上僅能獲得打擊結(jié)束后靶體的成坑形貌，很難連續(xù)測(cè)得彈坑的形成規(guī)律，因此難以建立起能精確描述成坑效應(yīng)的理論模型，而且對(duì)影響成坑的主要因素尚存在分歧。目前迫切需要對(duì)超高速侵徹過程中彈靶接觸面的介質(zhì)壓碎過程、塑性流動(dòng)軌跡、應(yīng)力波反射剝離效應(yīng)等進(jìn)行更加詳細(xì)的理論分析、試驗(yàn)量測(cè)和模擬研究工作，從而建立起更加具有物理意義的力學(xué)模型。

3.2 地沖擊演化規(guī)律

爆炸與撞擊都是能量的高速釋放轉(zhuǎn)化過程，屬于高壓力、高應(yīng)變率、高溫度的作用過程，鄧國(guó)強(qiáng)等[13,119]對(duì)超高速武器對(duì)地打擊效應(yīng)進(jìn)行了數(shù)值仿真，結(jié)果表明超高速對(duì)地打擊將形成塑性沖擊波，在彈坑形狀、沖擊波波形及衰減規(guī)律上高能裝藥的淺埋爆炸效應(yīng)均具有一定相似性，所不同的是超高速撞擊過程中彈體動(dòng)能具有定向性。

Oberbeck[120]通過實(shí)驗(yàn)表明：在動(dòng)能一定的情況下，可以用爆炸能量等于彈體動(dòng)能、且爆炸沖擊波壓力峰值等于超高速撞擊沖擊波峰值的淺埋爆炸來模擬超高速撞擊，實(shí)驗(yàn)利用長(zhǎng)4.6 mm、直徑6.4 mm、質(zhì)量0.435 g 的圓柱形鋁質(zhì)彈丸，以2.00 km/s 的速度撞擊石英砂靶體，彈體的動(dòng)能為870 J，彈丸對(duì)石英砂的沖擊壓力峰值為8 300 MPa。使用和彈丸形狀一模一樣、裝填密度為1 g/cm3的0.15 g TNT 炸藥（爆炸沖擊波壓力峰值為8 300 MPa），進(jìn)行淺埋爆炸，裝藥中心埋深分別為0、3.2、6.3、9.5、14.3 mm。對(duì)比發(fā)現(xiàn)，炸藥埋深為6.3 mm 左右（比例埋深）時(shí)，爆炸彈坑和撞擊彈坑匹配最好。

Baldwin[121]認(rèn)為，對(duì)于特定的彈靶組合和炸藥，等效模擬存在速度限制。以鋁質(zhì)彈丸高速撞擊石英砂為例，如果彈丸撞擊速度高達(dá)9.6 km/s，則彈丸在石英砂中產(chǎn)生的沖擊波壓力高達(dá)50 GPa，而裝填密度為1.6 g/cm3的TNT 炸藥爆炸沖擊波壓力峰值為23 GPa，同時(shí)要求二者的單位能量和沖擊波壓力峰值都相等有很大困難，用爆炸模擬將不再合適。

Holsapple[122]以彈坑體積相同為標(biāo)準(zhǔn)，提出了化學(xué)爆炸模擬超高速撞擊問題中“等效埋深”的概念，多種撞擊條件和化學(xué)爆炸條件下，石英砂中的等效埋深實(shí)驗(yàn)結(jié)果，并通過試驗(yàn)得出了定性結(jié)論：撞擊速度固定的情況下，隨著撞擊能量的增加，等效埋深減??；撞擊能量固定的條件下，隨著撞擊速度的增加，等效埋深減小。

王明洋等[32-33]認(rèn)為，從力學(xué)本質(zhì)上講，巖石中超高速撞擊與淺埋爆炸都屬于強(qiáng)動(dòng)載作用下材料動(dòng)力學(xué)行為和過程的問題，二者近源區(qū)產(chǎn)生強(qiáng)沖擊波，隨著傳播距離增加迅速衰減為短波、彈塑性波，波的傳播衰減規(guī)律相同，可以用相同的物理力學(xué)方程描述。隨后王明洋等以彈坑體積和形態(tài)為等效指標(biāo)（基于沖擊成坑深度確定等效裝藥埋深，基于沖擊成坑半徑確定彈坑拋擲指數(shù)），建立超高速撞擊與標(biāo)準(zhǔn)裝藥爆炸的能量等效關(guān)系：

式中： kδ為爆破多方指數(shù)，與炸藥和巖石性質(zhì)有關(guān)； Nδ=R/h 為彈坑形狀系數(shù)； R為表面成坑半徑。

與化爆、核爆一樣，超高速撞擊引起的地沖擊也可在一定范圍內(nèi)采用以下公式描述[32,123]：

式中：σpk為應(yīng)力波峰值，r 為撞擊點(diǎn)到考察點(diǎn)的距離，n 為衰減指數(shù)。衰減系數(shù) n的取值與巖石壓力狀態(tài)密切相關(guān)。Melosh[123]較為系統(tǒng)地闡述了地質(zhì)材料的超高速撞擊成坑及沖擊波傳播的過程，如圖1 所示，在速度10 km/s 的彈體的撞擊下，彈體和巖石靶體之間形成峰值壓力50 GPa 以上的沖擊波并向地下傳播，波的傳播規(guī)律如圖17 所示，在50 GPa 以上的強(qiáng)沖擊區(qū)波的壓力衰減指數(shù)約為3.6，在5～20 GPa 的過渡壓力范圍內(nèi)衰減指數(shù)約為1.4～1.8，在小變形彈塑性變形區(qū)衰減指數(shù)約為1.1～1.2（見圖18）。

哈努卡耶夫[124]總結(jié)了值與側(cè)壓力系數(shù)α 間的簡(jiǎn)單關(guān)系：

式中：正號(hào)對(duì)應(yīng)沖擊波傳播區(qū)域（高應(yīng)力流體狀態(tài)），負(fù)號(hào)對(duì)應(yīng)固體彈塑性波傳播區(qū)域（低應(yīng)力固體區(qū)）。在高應(yīng)力流體區(qū)， α ≈1， 故 n ≈3 ，在低應(yīng)力固體區(qū)， α ≈0.2～0.6，上述關(guān)系大致描述了n 值從近區(qū)（n＞2）向遠(yuǎn)區(qū)（n＜2）突躍式減小的現(xiàn)象。

借助于室內(nèi)縮比撞擊試驗(yàn)與預(yù)埋壓力傳感器的方法可以直接測(cè)得不同位置的壓力變化，進(jìn)而得到壓力隨距離衰減的真實(shí)過程。程怡豪[125]總結(jié)了不同文獻(xiàn)中實(shí)測(cè)或數(shù)值計(jì)算得到的衰減指數(shù)（見表4）。

對(duì)于巖土中爆炸應(yīng)力波波形計(jì)算，應(yīng)用最多的是文獻(xiàn)[126]提供的計(jì)算公式：

式中： α=1/ta, ta=r/ct為達(dá)到時(shí)間，tr為升壓時(shí)間，f 耦合系數(shù)，ρt為靶體材料密度，cP縱波速度，Q 為等效TNT 當(dāng)量；r 為爆心至測(cè)點(diǎn)距離， n為衰減指數(shù)。

圖 17 球形彈超高速撞擊下介質(zhì)中壓力分布[123]Fig. 17 Pressure distribution in medium under hypervelocity spherical projectile[123]

圖 18 峰值壓力隨距離衰減曲線[123]Fig. 18 Peak pressure decay with impact of distance[123]

表 4 地質(zhì)類材料地沖擊衰減指數(shù)nTable 4 Power exponent n for attenuation of ground shock in geological material

花崗巖中超高速撞擊地沖擊試驗(yàn)證實(shí)，按公式(50)計(jì)算的應(yīng)力波形與實(shí)測(cè)波形基本一致（見圖19～20）。

綜上，在超高速?gòu)楏w對(duì)地撞擊時(shí)，產(chǎn)生類似于淺埋裝藥爆炸的強(qiáng)地沖擊效應(yīng)，在進(jìn)行防護(hù)工程抗超高速動(dòng)能彈打擊最小安全防護(hù)層厚度設(shè)計(jì)時(shí)，除侵徹近區(qū)的局部破壞效應(yīng)（侵徹、成坑）外，還應(yīng)考慮由于超高速撞擊所引起的地沖擊荷載。近年來，在超高速動(dòng)能武器研制及工程防護(hù)需求的推動(dòng)下，不同學(xué)者對(duì)下地沖擊演化及壓力峰值衰減規(guī)律進(jìn)行了研究，建立了超高速撞擊巖石地沖擊效應(yīng)與淺埋爆炸的等效計(jì)算方法，但相關(guān)研究還未達(dá)到機(jī)理明晰的地步：（1）對(duì)于衰減系數(shù)n，目前仍主要采用基于試驗(yàn)擬合和數(shù)值模擬所得到的經(jīng)驗(yàn)系數(shù)，還未建立系數(shù)n 與應(yīng)力波、巖石力學(xué)參數(shù)的關(guān)聯(lián)關(guān)系；（2）由圖18 和公式（49）可知，參數(shù)n 隨距撞擊源（或爆心）比例距離的增加而急劇減小，這也造成衰減系數(shù)n 擬合的主觀性，在試驗(yàn)過程中測(cè)點(diǎn)位置、測(cè)點(diǎn)數(shù)目的變化均可造成n 值的較大差異；（3）目前超高速?gòu)楏w侵徹巖石的地沖擊效應(yīng)試驗(yàn)主要采用室內(nèi)縮比撞擊試驗(yàn)與預(yù)埋壓力傳感器的方法，由于巖石尺度小、地沖擊衰減快造成了試驗(yàn)測(cè)量的困難，同時(shí)在巖石中預(yù)埋壓力傳感器的方法破壞了巖石的整體性，也會(huì)對(duì)試驗(yàn)結(jié)果造成誤差，如何實(shí)現(xiàn)超高速撞擊過程中地沖擊傳播的高分變率全過程觀測(cè)仍是未來發(fā)展的方向；（4）超高速侵徹地沖擊效應(yīng)的數(shù)值模擬研究發(fā)現(xiàn)[13,119]，超高速成坑地沖擊具有定向傳播特性，同時(shí)由于彈靶界面移動(dòng)速度接近應(yīng)力波傳播速度，其追趕前驅(qū)應(yīng)力波會(huì)產(chǎn)生壓力波形的疊加倍增效應(yīng)，相關(guān)現(xiàn)象采用等效計(jì)算的方法難以準(zhǔn)確描述，需要建立更加準(zhǔn)確的計(jì)算模型。

圖 19 彈速3 558 m/s 時(shí)實(shí)測(cè)地沖擊壓力時(shí)程曲線[32]Fig. 19 The experiment time history curve of ground shock with impact velocity 3 558 m/s[32]

圖 20 彈速3 558 m/s 時(shí)按公式（50）計(jì)算的地沖擊壓力時(shí)程曲線[32]Fig. 20 the calculated time history curve of ground shock with impact velocity 3 558 m/s[32]

4 結(jié)論與展望

超高速武器是當(dāng)今軍事強(qiáng)國(guó)致力發(fā)展的下一代武器系統(tǒng)，近年來，國(guó)內(nèi)外學(xué)者圍繞超高速動(dòng)能武器打擊巖石的毀傷效應(yīng)與工程防護(hù)問題進(jìn)行了廣泛而細(xì)致的研究，取得了豐富的研究成果。本文中圍繞大速度范圍內(nèi)巖石介質(zhì)的侵徹效應(yīng)經(jīng)驗(yàn)公式、侵徹理論模型、成坑效應(yīng)和地沖擊效應(yīng)的研究成果進(jìn)行綜述，通過對(duì)相關(guān)文獻(xiàn)的查閱分析，得到的主要結(jié)論和建議如下：

（1）大速度范圍內(nèi)桿形彈的侵徹效應(yīng)具有明顯的多階段特征。隨著撞擊速度增加，彈體將經(jīng)歷從變形可忽略的階段（僅有少量質(zhì)量損失）向侵蝕階段（即彈體長(zhǎng)度嚴(yán)重縮短）的轉(zhuǎn)變，靶體材料經(jīng)歷從固體侵徹、內(nèi)摩擦擬流體侵徹項(xiàng)動(dòng)力學(xué)侵徹行為的轉(zhuǎn)變，侵徹深度經(jīng)歷“迅速增加—逆減—緩慢增加—趨于流體動(dòng)力學(xué)極限”的過程，不同階段的物理機(jī)理不同，在進(jìn)行侵徹效應(yīng)計(jì)算時(shí)需要考慮所采用的物理力學(xué)模型對(duì)不同彈靶狀態(tài)的適用性問題。

（2）近年來圍繞巖石和混凝土介質(zhì)的侵徹深度計(jì)算問題開展了大量研究，建立了空腔膨脹、修正流體動(dòng)力學(xué)、內(nèi)摩擦流體彈塑性等理論模型，實(shí)現(xiàn)了對(duì)大速度范圍內(nèi)侵徹深度變化規(guī)律的描述。但相關(guān)研究仍或多或少存在表征不夠精細(xì)的問題，一些參數(shù)的獲取仍然具有明顯的經(jīng)驗(yàn)性，甚至存在一定爭(zhēng)議，需要更進(jìn)一步的深入寬廣應(yīng)變率范圍內(nèi)巖石的動(dòng)態(tài)力學(xué)行為，盡可能的了解不同加載范圍材料的實(shí)際性能。

（3）超高速侵徹過程下的成坑和地沖擊效應(yīng)不可忽視，但目前相關(guān)理論研究仍處于起步階段，主要采用半理論半經(jīng)驗(yàn)的計(jì)算公式，需要從超高速侵徹成坑、地沖擊演化的三維效應(yīng)和力學(xué)機(jī)制出發(fā)，構(gòu)建更加完善、更加細(xì)致的演化模型。

（4）由于巖石的工程地質(zhì)屬性，其變形破壞具有顯著的尺度效應(yīng)，由于發(fā)射手段限制，目前超高速侵徹效應(yīng)理論的研究缺乏大尺寸彈體的超高速侵徹實(shí)驗(yàn)驗(yàn)證，需要發(fā)展更為先進(jìn)的試驗(yàn)設(shè)備和技術(shù)，針對(duì)不同彈靶材料開展大尺度、速度范圍的室內(nèi)相似模擬與原型試驗(yàn)，建立寬廣尺度范圍內(nèi)彈體侵徹的相似規(guī)律。