沖擊作用下CL-20/HMX 共晶力-熱-化學(xué)耦合響應(yīng)的動力學(xué)模擬

王 寧，蘇 晶，關(guān)紅波，成 杰，程立國，李 軍

（1. 湖北航天化學(xué)技術(shù)研究所，湖北 襄陽 441003；2. 航天化學(xué)動力技術(shù)重點實驗室，湖北 襄陽 441003；3. 湖北航鵬化學(xué)動力科技有限責(zé)任公司，湖北 襄陽 441003）

1 引言

含能材料的能量和安全性矛盾一直以來受到廣泛關(guān)注，含能共晶炸藥是解決這一矛盾的主要途徑之一，其中較為典型代表是CL-20/HMX 共晶炸藥［1-4］，該炸藥的爆轟性能優(yōu)于β-HMX，感度與β-HMX 相當(dāng)［2］，對其感度機理的深入研究有利于新型含能共晶炸藥的設(shè)計和開發(fā)。含能材料分子、晶體結(jié)構(gòu)與其感度的關(guān)系，國內(nèi)外開展了大量研究［4-8］。對于CL-20/HMX 共晶炸藥的撞擊感度，熱力學(xué)引發(fā)鍵鍵級判據(jù)和動力學(xué)活化能判據(jù)［9-10］在單質(zhì)含能材料中應(yīng)用較為廣泛。由于分子層級混合的含能共晶仍屬于不同分子結(jié)構(gòu)的混合，因此這些判據(jù)還在進(jìn)一步發(fā)展。對CL-20/HMX 共晶的穩(wěn)定性，Bolton 等［2］認(rèn)為與純組分相比，共晶中的C—H…O 氫鍵數(shù)目明顯增多，鍵長減小。陶俊等［11］根據(jù)最大鍵長和結(jié)合能，以范德華為主的作用力“敏化”了CL-20/HMX 共混體系的結(jié)構(gòu)，而根據(jù)內(nèi)聚能密度和徑向分布函數(shù)，得出CL-20/HMX 共晶體系由于存在長度相對較短的CH…O 氫鍵而導(dǎo)致其感度較低。若考慮晶體形貌或缺陷影響，用于預(yù)示感度的最大引發(fā)鍵鍵級判據(jù)和分子間弱相互作用強弱等將受到限制。隨著計算能力提高，從動力學(xué)角度模擬晶體結(jié)構(gòu)和感度的關(guān)系成為可能。

目前，立足于含能材料分子結(jié)構(gòu)，在晶體結(jié)構(gòu)、顆粒等空間上和介觀-微觀等時間尺度上，對強沖擊下含能材料發(fā)生的力-熱-化學(xué)變化耦合過程的研究已取得較大進(jìn)展。Jaramillo 等［12］、Cawkwell 等［13-14］和Bedrov 等［15］研究了RDX 不同晶面受到?jīng)_擊加載時發(fā)生的塑性行為，詳細(xì)分析了位錯、剪切帶和分子構(gòu)型的變化規(guī)律。Eason 等［16］分析了沖擊誘導(dǎo)PETN 時發(fā)生的彈塑性行為，而Shan 等［17］則分析了沿PETN 不同晶面沖擊加載時化學(xué)反應(yīng)過程。此外，孔洞塌縮形成熱點的機制也在不斷深入，Wood 等［18］分析了含孔洞缺陷的RDX 發(fā)生沖擊轉(zhuǎn)爆燃的過程，利用流體動力學(xué)方法和反應(yīng)分子動力學(xué)方法（ReaxFF-MD），詳細(xì)分析了HNS 晶體中孔洞塌縮形成熱點的過程［19］，而Shan 等［20］則在更長時間內(nèi)分析了孔洞塌縮后形成的非對稱熱點區(qū)域。然而含能共晶炸藥沖擊過程模擬未見報道，因此，本研究著重分析了沖擊作用下，有無缺陷CL-20/HMX 共晶發(fā)生的力-熱結(jié)構(gòu)變化和隨后的化學(xué)反應(yīng)過程，為進(jìn)一步分析沖擊波感度機理提供必要基礎(chǔ)。

2 計算模型與計算過程

2.1 CL-20/HMX 共晶初始模型及結(jié)構(gòu)優(yōu)化

含能材料熱分解和沖擊壓縮過程的研究常采用ReaxFF-MD 方法。由于不同時期開發(fā)的力場有差異，因此首先需要計算研究對象超胞結(jié)構(gòu)來驗證力場的適用 性。本 研 究 依 據(jù)ε-CL-20 晶 體［21］、β-HMX［22］和CL-20/HMX 共 晶［2］的X 射 線 單 晶 衍 射 結(jié) 果，建 立 了4×2×2 的CL-20 超胞、5×3×4 的HMX 超胞和2×3×2的CL-20/HMX 共晶超胞，再采用ReaxFF-MD 方法［23］進(jìn)行不同晶體的分子動力學(xué)優(yōu)化，將計算結(jié)果和試驗結(jié)果對比以篩選合適力場。

可能滿足CL-20/HMX 共晶計算的力場有ReaxFF-2005［24］、ReaxFF-lg［25］、ReaxFF-2014［18］和ReaxFF-2018［19］等，其中ReaxFF-lg 因添加了低梯度修正項（lg），能夠更準(zhǔn)確描述晶格參數(shù)和密度，在ReaxFF-2014 的基礎(chǔ)上，ReaxFF-2018 增加了lg 項，并兼顧了多種含能材料反應(yīng)途徑和燃燒過程。因此，本研究利用ReaxFF-lg 和ReaxFF-2018 分別計算CL-20、HMX 和CL-20/HMX 晶體結(jié)構(gòu)參數(shù)，以驗證力場適用性。

分子動力學(xué)優(yōu)化過程分為三步，即先后運行20 ps的T=5 K、正則系綜NVT 和T=300 K、NVT 進(jìn)行弛豫，然后選擇等溫等壓系綜NPT（T=300 K、p=1.01×105Pa）進(jìn)行50 ps 分子動力學(xué)模擬。采用Nose-Hoover 控溫和控壓［26］，時間步長為0.1 fs。

2.2 不同晶向CL-20/HMX 共晶模型的建立及沖擊加載過程細(xì)節(jié)

利用Materials Studio 中morphology 模塊［27］預(yù)測真空狀態(tài)下CL-20/HMX 共晶形貌，發(fā)現(xiàn)（100）晶面面積最大。Ghosh 等［4］根據(jù)粉末XRD 和晶體形貌認(rèn)為（111）晶面是重要晶面。因此，選擇（100）和（111）晶面作為沖擊加載的主要方向。

考慮模擬體系大小、晶面的尺寸匹配和計算資源等因素，沿CL-20/HMX共晶的［100］晶向構(gòu)建了10×1×24的模型CL-20/HMX（100），沿［111］晶向構(gòu)建了5×1×48的模型CL-20/HMX（111）；沿CL-20 的［011］晶向構(gòu)建了1×10×24 的模型CL-20（011），沿HMX 的［011］晶向構(gòu)建了1×10×24 的模型HMX（011）。具體尺寸如表1 所示。為減少因周期性邊界條件而可能增加的非鍵弱相互作用，在z 向均增加了2 nm 的真空層。

表1 不同晶向CL-20/HMX 共晶模型詳細(xì)信息Table 1 Detailed information of various CL-20/HMX models

為分析CL-20/HMX 共晶在沖擊加載下發(fā)生的力-熱結(jié)構(gòu)變化和化學(xué)反應(yīng)，本研究先對構(gòu)建的CL-20/HMX 共晶模型結(jié)構(gòu)優(yōu)化，再按照動量鏡原理［15］開展了沖擊加載模擬。沖擊加載過程如下：?。?00］晶向的兩層和［111］晶向的四層重復(fù)單元作為固定壁，壁厚約2～3 nm，其余部分以質(zhì)點速度2 km·s-1［20，28］撞擊固定壁，形成反射沖擊波。當(dāng)反射沖擊波傳播到尾端時，采用吸收波陣面方法［14］，繼續(xù)在NVE 系綜下開展分子動力學(xué)模擬。整個模擬過程的時間步在0.1～0.01 fs之間調(diào)整，滿足每步電荷收斂平衡需要和捕捉到更全面的化學(xué)反應(yīng)產(chǎn)物。通過設(shè)定不同種類原子對的鍵級閾 值［29］區(qū) 分 分子碎片，后采用Sergey Zybin 提 供的mol_fra.c 源代碼［30］處理。所有反應(yīng)分子動力學(xué)模擬過程均通過LAMMPS 軟件［31］完成。

2.3 含孔洞缺陷的CL-20/HMX 共晶模型的建立

孔洞缺陷是含能材料中最常見的缺陷類型之一。由于孔洞缺陷大小和沖擊強度對孔洞塌縮誘導(dǎo)的熱點形成有顯著影響，因此結(jié)合Shan 等［20］的研究結(jié)果，本研究建立了含孔洞直徑為20 nm、沿［100］晶向的模型CL-20/HMX（100）-void。為便于分析，還建立了含同樣孔洞尺寸的沿CL-20［011］晶向的模型CL-20（011）-void，沿CL-20/HMX 共晶［100］晶向的無缺陷模 型CL-20/HMX（100）-large，沿CL-20/HMX 共 晶［111］晶 向 的無缺陷模 型CL-20/HMX（111）-large。具體尺寸如表2 所示。

表2 ε-CL-20、β-HMX 和CL-20/HMX 晶胞參數(shù)的實驗值和模擬結(jié)果Table 2 Comparison of the experimental and calculated lattice parameters of CL-20，HMX，and CL-20/HMX

3 結(jié)果與討論

3.1 CL-20/HMX 力場適用性分析

為篩選適用于CL-20/HMX 共晶沖擊過程的力場，需正確描述CL-20/HMX 共晶結(jié)構(gòu)。為此，本研究通過ReaxFF-lg 和ReaxFF-2018 分 別 優(yōu) 化 了ε-CL-20、β-HMX 和CL-20/HMX 超胞結(jié)構(gòu)，獲取CL-20/HMX 共晶結(jié)構(gòu)，各物質(zhì)的晶胞參數(shù)和密度優(yōu)化結(jié)果見表2。由表2 可知，ReaxFF-lg 和ReaxFF-2018 這兩套 力 場 在描述含能材料晶體參數(shù)和密度時均具有較好的計算精度。

除優(yōu)化CL-20/HMX 的共晶結(jié)構(gòu)，為研究共晶的動態(tài)加載過程，還需計算相關(guān)高壓物理特性，包括靜水壓和體積壓縮比p-v/v0和Hugoniot 曲線相關(guān)的p-v/v0和T-v/v0及us-up關(guān)系等［20］。由 于CL-20/HMX 共晶 的 這些物理特性缺乏試驗數(shù)據(jù)，故本研究以高精度的色散修正的DFT-D［32］計算結(jié)果為參考，對比分析了ReaxFF-lg 和ReaxFF-2018 描述CL-20/HMX 共晶高壓物理特性的準(zhǔn)確性，結(jié)果如圖1～圖3所示。由圖1、圖2a和圖3 可知，與ReaxFF-lg 相比，采用ReaxFF-2018 計算值更接近于DFT 計算值。 由圖2b 可知，ReaxFF-2018 計算的高壓縮比下CL-20/HMX 共晶溫度偏離未反應(yīng)Hugoniot 曲線，這與部分物質(zhì)發(fā)生反應(yīng)有關(guān)，比如質(zhì)點速度超過1.787 km·s-1時，類似結(jié)果也反 映 在Islam 等［33］的 研 究 工 作 中。 因 此，選 取ReaxFF-2018 開展質(zhì)點速度為2 km·s-1時CL-20/HMX共晶力-熱-化學(xué)反應(yīng)過程研究。

圖1 CL-20/HMX 共晶的冷靜壓下體積壓縮比和靜水壓關(guān)系Fig.1 Relationship between compressive volume ratio and pressure for CL-20/HMX cocrystal

圖2 CL-20/HMX 共晶受沖擊后（a）p-v/v0和（b）T-v/v0關(guān)系Fig.2 Curves of p-v/v0（a）and T-v/v0（b）of CL-20/HMX cocrystal under shock loading

圖3 計算獲得的CL-20/HMX 共晶的us-up關(guān)系Fig.3 Calculated us-up relationship of CL-20/HMX cocrystal

3.2 CL-20/HMX 不同晶向響應(yīng)特性分析

為考察CL-20/HMX 共晶是否具有明顯的各向異性，分別沿［100］和［111］晶向沖擊加載了CL-20/HMX 共晶。初始質(zhì)點速度設(shè)定為2 km·s-1，其溫度和初始物質(zhì)（包括CL-20 和HMX）分解量隨時間的變化關(guān)系如圖4 所示。為對比分析，同時計算了沿［011］晶向沖擊加載的CL-20 和HMX。

圖4 沿［100］和［111］沖擊加載CL-20/HMX 共晶時（a）溫度和（b）未分解分子隨時間的變化規(guī)律Fig.4 Evolution of temperature（a）and percentage of decomposed molecules（b）vs time along the［100］and［111］shocks in CL-20/HMX cocrystal

由圖4a 可知，整個快速沖擊響應(yīng)過程持續(xù)了約5 ps，不同模型沖擊加載時間與模型加載方向長度有關(guān)?？焖贈_擊響應(yīng)后發(fā)生緩慢化學(xué)分解過程，表現(xiàn)在不同模型的溫度均呈現(xiàn)不同程度的升高。CL-20、HMX 和CL-20/HMX 共晶的溫度增長速率大小順序為：CL-20＞CL-20/HMX（100）＞HMX＞CL-20/HMX（111）。由 圖4b 可 知，CL-20、HMX 和CL-20/HMX 共 晶 中CL-20 分子的分解速率大小順序為CL-20＞CL-20/HMX（100）＞CL-20/HMX（111），HMX 分 子 則 為CL-20/HMX（100）＞HMX＞CL-20/HMX（111）。因此單純CL-20 的分解速度最快，而CL-20/HMX 和HMX的分解速度相當(dāng)，CL-20/HMX（111）模型的分解速度甚至比單純HMX 的還慢，可以初步解釋CL-20/HMX共晶的沖擊波感度和HMX 更接近的實驗現(xiàn)象。與［111］晶向相比，沿［100］晶向加載時CL-20/HMX 共晶分解更快，表明出明顯的沖擊各向異性。

另外，當(dāng)快速沖擊加載結(jié)束時模型會達(dá)到最大壓縮點，此時若只施加周期性邊界條件繼續(xù)開展分子動力學(xué)模擬，將因為末端產(chǎn)生的反向稀疏波而使得體系變得松散，因此在最大壓縮點時額外施加一靜止活塞［14］，可以有效阻止稀疏波進(jìn)入波陣面后區(qū)域，從而能夠進(jìn)一步分析結(jié)構(gòu)塑性變形和化學(xué)反應(yīng)過程。

為研究CL-20/HMX 共晶沖擊各向異性產(chǎn)生機理，統(tǒng)計分析了最大壓縮點時CL-20/HMX（100）和CL-20/HMX（111）模型中完整的CL-20 和HMX 分子的分布，計算結(jié)果如圖5 所示。由分子堆積過程可知，每兩層CL-20 分子和一層HMX 分子交替堆疊形成共晶，HMX 和CL-20 之間因為氫鍵作用而較穩(wěn)定的組裝在一起。由圖5a 中沿［100］晶向沖擊加載時，作用力將垂直于CL-20 和HMX 分子層而逐層傳遞，而由圖5b 可知，沿［111］晶向沖擊加載時，CL-20 分子層和HMX 分子層之間產(chǎn)生滑移。與［100］晶向的這種顯著差異，直觀表現(xiàn)為即使在［111］晶向沖擊加載端也有較多保持完整的CL-20 和HMX 分子。結(jié)合圖4a 可知，與沿［100］晶向沖擊加載相比，通過不同活性分子的交替排布和一定量的滑移等能夠有效分散沖擊加載能量，致使CL-20/HMX（111）模型的響應(yīng)速度變慢，表現(xiàn)出明顯的沖擊各向異性。

圖5 沿［100］（a）和［111］（b）沖擊加載且處于最大壓縮狀態(tài)的CL-20/HMX 共晶結(jié)構(gòu)（紅色表示完整的CL-20 分子，綠色表示完整的HMX 分子，其余屬于碎片分子）Fig.5 Structures of CL-20/HMX models at the maximum moment along［100］（a）and［111］（b）shock directions（Color ranks：red-intact CL-20 molecules，green- intact HMX molecules，others-fragments）

為進(jìn)一步分析沿不同晶向沖擊CL-20/HMX 共晶模型時能否發(fā)生滑移，計算了沿x 向和z 向的局部位移分布（如圖6）和z 向的局部速度分布（如圖7）。圖6a和圖6b 中CL-20/HMX（100）模型處于3.2 ps，波陣面位于20 nm，而圖6c 和圖6d 中CL-20/HMX（100）處于3.8 ps，波陣面位于17.5 nm，兩模型的壓縮量均為最大壓縮量的60%。在CL-20/HMX（111）模型中，z向起始端屬于固定壁位置，沿z 向的2.5～17.5 nm 范圍屬于波后受壓區(qū)域，17.5～32.5 nm 屬于未壓縮的CL-20/HMX 共晶區(qū)域，32.5～40 nm 范圍屬于空隙區(qū)。由圖6a 可知，沿［100］晶向加載時，受壓縮區(qū)沿z 向位移方向和沖擊波傳播方向一致，而圖6b 中沿x 向位移方向并無明顯規(guī)律。由圖6c 和圖6d 可知，沿［111］晶向加載時，受壓縮區(qū)沿z 向和x 向位移方向均有正有負(fù)，而且位移數(shù)值也比沿［100］晶向的大。相似規(guī)律也反映在圖7 中。這些結(jié)果表明，沿［111］晶向沖擊加載后，受壓縮區(qū)的確發(fā)生了明顯的滑移，最終延緩了反應(yīng)進(jìn)程，而沿［100］晶向沖擊加載時滑移并不顯著。

圖7 沿不同晶向沖擊加載的CL-20/HMX（100）和CL-20/HMX（111）模型的局部速度變化規(guī)律Fig.7 Variations of local velocities along different shock directions in CL-20/HMX（100）and CL-20/HMX（111）models

3.3 孔洞缺陷對沖擊過程的影響

在分析無缺陷CL-20/HMX 模型的基礎(chǔ)上，進(jìn)一步考慮孔洞缺陷對沖擊過程中力-熱耦合過程影響。Φ20 nm 孔洞存在下溫度、壓力以及初始物質(zhì)分解量隨時間的變化關(guān)系分別如圖8 和圖9 所示。由圖4a、圖8a、圖4b、圖9 可知，沿［100］和［111］沖擊加載無缺陷CL-20/HMX 共晶時，溫度和不同分子分解量相對變化規(guī)律一致，即沿［111］沖擊加載時溫度（主要和粒子質(zhì)點速度相關(guān)）降低最多，壓力增加緩慢，CL-20 和HMX 分子分解速率最低，表明不同規(guī)模模型對計算結(jié)果的影響較小。然而，當(dāng)存在Φ20 nm 孔洞缺陷時，伴隨孔洞的壓縮，體系溫度出現(xiàn)先增后減的“鼓包”現(xiàn)象，體系壓力則出現(xiàn)先減后增，對CL-20 和HMX 分子的分解產(chǎn)生一定的影響。下面結(jié)合孔洞壓縮過程瞬態(tài)圖做進(jìn)一步分析。

沿［100］沖擊加載含Φ20 nm 孔洞的CL-20/HMX共晶時，孔洞受壓縮初期，模型中各物理量分布如圖10 所示。沿z 軸正向，含缺陷CL-20/HMX 共晶分別由固定壁、含孔洞的CL-20/HMX 共晶受壓區(qū)和運動區(qū)、空隙區(qū)等四部分組成。由圖10a 可知，波陣面在碰到孔洞時被隔斷，然而以沖擊波傳播方向為參考，孔洞上游壁面和下游壁面附近區(qū)域溫度均較高（圖10b），這是由于這兩部分局部區(qū)域?qū)?yīng)的粒子質(zhì)點速度和位移均因孔洞存在而增加了自由度，受周圍分子的約束減弱而動能增加的緣故。

圖8 沿［100］和［111］沖擊加載含Φ20 nm 缺陷CL-20/HMX 共晶模型溫度和壓力隨時間的變化規(guī)律Fig.8 Evolution of temperature and pressure along the［100］and［111］shocks in CL-20/HMX cocrystal with 20 nm diameter void

圖9 不同CL-20/HMX 模型中未分解CL-20 和HMX 分子數(shù)的變化規(guī)律Fig.9 Evolution of percentage of undecomposed CL-20 and HMX molecules in different CL-20/HMX models

為詳細(xì)分析沖擊加載CL-20/HMX［100］模型的結(jié)構(gòu)變化，根據(jù)孔洞變形階段，計算了幾組典型時刻的密度、溫度和壓力的瞬態(tài)分布，結(jié)果如圖11 所示。選取的典型時刻已標(biāo)注于圖8 中，即2.86 ps 對應(yīng)于圖10中時刻，5.44 ps 對應(yīng)于孔洞臨近閉合時刻，6.62 ps 對應(yīng)于近球形孔洞塌縮“熱點”產(chǎn)生的應(yīng)力波正追趕沖擊波的時刻，8.30 ps 對應(yīng)于孔洞塌縮熱點產(chǎn)生的應(yīng)力波和之前沖擊波復(fù)合形成的加強型沖擊波繼續(xù)向前傳播的某一時刻。由圖8～圖11a 所展示的孔洞塌縮過程和后效分析可知，孔洞塌縮過程導(dǎo)致孔洞上游壁面和下游壁面獲得較高的質(zhì)點速度，兩者是相對運動，孔洞塌縮的過程產(chǎn)生較高的溫度、降低局部壓力、減弱CL-20 和HMX 分子的分解。在圖11b 中，孔洞閉合后，以孔洞閉合位置為中心產(chǎn)生新熱點，此時CL-20 和HMX 分子的分解速率明顯提高，形成了更高的溫度和壓力區(qū)并向周圍擴散。在圖11c 中，當(dāng)孔洞塌縮熱點產(chǎn)生的應(yīng)力波和之前沖擊波復(fù)合時形成了加強型沖擊波，繼續(xù)向前傳播。從新熱點到兩波平齊，含Φ20 nm孔洞的CL-20 模型需要1.61 ps，而CL-20/HMX 模型則需要1.92 ps。

圖10 沿［100］方向沖擊加載含Φ20 nm 缺陷CL-20/HMX 共晶時孔洞受壓縮初期各物理量的瞬態(tài)變化Fig.10 Evolution of various parameters of CL-20/HMX with 20 nm diameter void along the［100］shock direction at initial stage，which includes（a）density；（b）temperature；（c）pressure；（d）displacement

由孔洞塌縮形成熱點過程可知，孔洞缺陷首先造成受力結(jié)構(gòu)變化，以2 km·s-1質(zhì)點速度沖擊下Φ20 nm 孔洞周圍結(jié)構(gòu)沒有產(chǎn)生射流現(xiàn)象，而是粘塑性收縮過程，然后孔洞閉合產(chǎn)生的高溫高壓條件和結(jié)構(gòu)上的粘塑性變形均有效促使CL-20 和HMX 分子發(fā)生快速的化學(xué)分解，孔洞塌縮新熱點的形成反過來進(jìn)一步增強了沖擊加載過程。

4 結(jié)論

采用反應(yīng)分子動力學(xué)方法，探討了有無缺陷CL-20/HMX 共晶在沖擊加載下發(fā)生的力-熱結(jié)構(gòu)變化和隨后的初始化學(xué)反應(yīng)過程。結(jié)果表明：

（1）ReaxFF-lg 和ReaxFF-2018 均能夠較好的描述CL-20、HMX 以及CL-20/HMX 共晶的晶格參數(shù)、密度。與ReaxFF-lg 相比，利用ReaxFF-2018 計算的CL-20/HMX 共晶高壓物理特性值更接近于DFT-D 計算值。

（2）利用動量鏡原理產(chǎn)生平面反射波后，在最大壓縮點增加一靜止活塞，可有效阻止稀疏波進(jìn)入波后區(qū)域，有利于分析結(jié)構(gòu)塑性變形情況或化學(xué)反應(yīng)過程，不僅可實現(xiàn)高速沖擊加載過程模擬，也能夠?qū)崿F(xiàn)低速撞擊過程模擬，具有較大應(yīng)用潛力。

（3）沖擊加載無缺陷CL-20、HMX 和CL-20/HMX共晶時，CL-20 比HMX 的分解速度快，CL-20/HMX 共晶的分解速度接近HMX。和［111］晶向相比，沿［100］晶向沖擊加載CL-20/HMX 共晶分解更快，這與CL-20 和HMX 分子層的交替排布以及能否產(chǎn)生滑移等因素有關(guān)。

（4）以2 km·s-1質(zhì)點速度沖擊下Φ20 nm 孔洞周圍結(jié)構(gòu)沒有產(chǎn)生射流現(xiàn)象，而是粘塑性收縮過程?？锥此s形成的高溫高壓條件和結(jié)構(gòu)上的粘塑性變形有效促使CL-20 和HMX 分子發(fā)生快速的分解，孔洞塌縮新熱點的形成進(jìn)一步增強了沖擊加載過程。