基于廣義波阻抗梯度飛片的準(zhǔn)等熵壓縮技術(shù)*

陳子博，謝普初，劉東升，陳 偉，王永剛

（寧波大學(xué)沖擊與安全工程教育部重點(diǎn)實(shí)驗(yàn)室，浙江 寧波 315211）

梯度功能材料是兩種或多種材料復(fù)合且成分和結(jié)構(gòu)呈連續(xù)梯度變化的一種新型復(fù)合材料。波阻抗梯度功能材料（wave impedance functionally graded material，WIFGM）是一類較特殊的梯度功能材料，該材料引起人們的關(guān)注始于材料高壓本構(gòu)和狀態(tài)方程研究[1-2]。在輕氣炮實(shí)驗(yàn)中，波阻抗梯度功能材料制作的飛片在高速碰撞靶材后，可以在靶材中產(chǎn)生具有緩慢上升前沿的壓縮加載波，實(shí)現(xiàn)對(duì)靶材的準(zhǔn)等熵壓縮加載，熵增和溫升遠(yuǎn)小于沖擊絕熱壓縮情況。準(zhǔn)等熵壓縮加載填補(bǔ)了準(zhǔn)靜態(tài)等熵加載和沖擊絕熱加載之間空白，提供了一種考察物質(zhì)壓縮形態(tài)（尤其是偏離Hugoniot狀態(tài)）的新實(shí)驗(yàn)途徑，可以得到范圍更廣、壓縮度更高、精度更高的材料物態(tài)方程參數(shù)。另外，把波阻抗梯度飛片的準(zhǔn)等熵壓縮性和二級(jí)輕氣炮相結(jié)合可以發(fā)展三級(jí)炮超高速發(fā)射技術(shù)，彈丸速度可以達(dá)到10 km/s以上[3]，這在航天器空間碎片防護(hù)領(lǐng)域具有重要的應(yīng)用背景。

20世紀(jì)80年代，美國圣地亞核武器國家實(shí)驗(yàn)室首先開展了波阻抗梯度材料的制備技術(shù)研究。Barker等[1]采用顆粒共沉降成型結(jié)合高壓粘接方法制備了Ta/Cu/Al/TPX金屬-聚合物體系的波阻抗梯度飛片，可以在靶材中產(chǎn)生準(zhǔn)等熵壓縮波，但該波阻抗梯度飛片的力學(xué)強(qiáng)度較差，在發(fā)射過程中飛片易發(fā)生破壞。Chhabildas等[2]采用焊接方法制備了Ta/Cu/Ti/Al/Mg/TPX體系疊層飛片，但該飛片的波阻抗沿厚度方向是以階躍方式變化的，導(dǎo)致其產(chǎn)生的壓力剖面平順性較差。后來，美國勞倫斯利弗莫爾國家實(shí)驗(yàn)室Martin等[4]采用流延成型法制備了7種組分和19種組分的Mg/Cu體系波阻抗梯度材料，獲得了平滑的壓力剖面。近來，Yep等[5]進(jìn)一步發(fā)展了流延成型技術(shù)，制備了近100種組分的Mg/Cu/W體系梯度飛片，其具有更寬廣波阻抗變化范圍，可以實(shí)現(xiàn)500 GPa壓力值的準(zhǔn)等熵壓縮。從20世紀(jì)90年代開始，武漢理工大學(xué)和中國工程物理研究院流體物理研究所合作開展波阻抗梯度材料制備技術(shù)研究[6-8]，發(fā)展了放電等離子快速粉末燒結(jié)法，研制出的梯度飛片已成功應(yīng)用于準(zhǔn)等熵壓縮加載、超高速驅(qū)動(dòng)、太帕量級(jí)壓力的高精度平面加載以及復(fù)雜加、卸載等動(dòng)高壓實(shí)驗(yàn)技術(shù)領(lǐng)域。傳統(tǒng)的波阻抗梯度功能材料制備多是選擇多種不同密度不同波速的金屬材料進(jìn)行復(fù)合。為了得到材料波阻抗準(zhǔn)連續(xù)變化，需要進(jìn)行非常多組分配比，鋪層和燒結(jié)制作工序異常繁雜，質(zhì)量也難以精確控制，成品率低，價(jià)格昂貴，這大大制約了波阻抗梯度功能材料實(shí)際應(yīng)用，迫切需要發(fā)展新的制備工藝。最近，Taylar等[9]采用激光立體光刻技術(shù)制備了波阻抗梯度陶瓷飛片，實(shí)現(xiàn)對(duì)平板撞擊實(shí)驗(yàn)中沖擊加載波形的有效調(diào)控。

增材制造(additive manufacturing, AM)技術(shù)是基于計(jì)算機(jī)輔助設(shè)計(jì)數(shù)據(jù)，采用材料逐層累積的方法來制造零件[10-11]。材料結(jié)構(gòu)的計(jì)算機(jī)三維自由設(shè)計(jì)和增材制造技術(shù)的完美結(jié)合為新型功能材料或結(jié)構(gòu)材料的研發(fā)與制備提供新的途徑。本文中擬采用激光選區(qū)熔化（selective laser melting，SLM）金屬增材制造技術(shù)，通過幾何結(jié)構(gòu)設(shè)計(jì)來研發(fā)新型的廣義波阻抗梯度飛片材料。首先，基于應(yīng)力波傳播理論來探討波阻抗梯度飛片的幾何構(gòu)型設(shè)計(jì)；隨后對(duì)“結(jié)構(gòu)設(shè)計(jì)型”波阻抗梯度飛片的準(zhǔn)等熵壓縮加載特性進(jìn)行數(shù)值模擬，重點(diǎn)討論結(jié)構(gòu)設(shè)計(jì)的關(guān)鍵參數(shù)和碰撞速度對(duì)加載波形的影響規(guī)律；最后，采用增材制造技術(shù)制備波阻抗梯度飛片，開展氣炮實(shí)驗(yàn)，對(duì)比分析了實(shí)驗(yàn)與數(shù)值計(jì)算結(jié)果，驗(yàn)證新型波阻抗梯度飛片的設(shè)計(jì)可行性和適用性。

1 波阻抗梯度飛片的結(jié)構(gòu)設(shè)計(jì)

材料密度ρ0與波速c的乘積定義為材料波阻抗。在研究變截面桿中波傳播問題時(shí)人們引入廣義波阻抗的概念，定義為材料密度ρ0、波速c以及截面面積A的乘積。由此來看，除了改變材料密度和波速之外，還可以通過改變截面面積來實(shí)現(xiàn)廣義波阻抗的梯度變化。王禮立等[12]通過引入廣義波阻抗把錐形桿中的縱波傳播問題近似化為一維應(yīng)力平面縱波問題來處理，通過特征線方法探討了應(yīng)力波在錐形桿中的傳播規(guī)律，結(jié)果表明，錐形桿具有調(diào)節(jié)應(yīng)力波波形和強(qiáng)弱的特性。羅鑫等[13]對(duì)錐形桿共軸撞擊進(jìn)行數(shù)值分析，結(jié)果表明，撞擊所產(chǎn)生的加載波陣面具有較長的前沿升時(shí)，同時(shí)還消除了波形振蕩。陶俊林等[14]在SHPB實(shí)驗(yàn)中將圓柱形子彈改變?yōu)殄F形子彈，實(shí)現(xiàn)對(duì)一維應(yīng)力條件下加載波形的調(diào)整，上升沿時(shí)間顯著增大。受此啟發(fā)，在輕氣炮實(shí)驗(yàn)中也采用類似的設(shè)計(jì)思想。輕氣炮實(shí)驗(yàn)系統(tǒng)中采用薄板飛片作為加載手段，以便在樣品中產(chǎn)生一維沖擊波，因此不可以直接把薄板飛片改變?yōu)殄F形飛片。這里，我們把飛片分為兩部分，一部分作為飛片基座，然后在基座上密排疊加許多小正四棱錐臺(tái)，猶如基座上長出一根根變截面針刺，簡稱為“針床型”波阻抗梯度飛片，其幾何結(jié)構(gòu)示意圖如圖1所示，圖中e為基座高，h為四棱錐臺(tái)高，m為正四棱錐臺(tái)下底邊寬，n為正四棱錐臺(tái)上頂邊寬。當(dāng)“針床型”波阻抗梯度飛片高速撞擊樣品時(shí)，每一個(gè)小正四棱錐臺(tái)都可以看作“點(diǎn)”加載脈沖源，會(huì)產(chǎn)生一系列具有緩慢上升前沿的類球面波，隨著傳播距離增大球面波將相互疊加，最終可以得到近似平面的加載波形，從而滿足一維應(yīng)變加載條件。

圖1 “針床形”廣義波阻抗梯度飛片幾何結(jié)構(gòu)示意圖Fig.1 Schematic diagram of the “needle-bed” generalized wave impedance gradient flyer

2 計(jì)算模型

為了驗(yàn)證上述設(shè)計(jì)思想的可行性，首先開展“針床型”波阻抗梯度飛片高速撞擊靶板的數(shù)值模擬。在數(shù)值計(jì)算之前，還需要選擇合適的算法。如果采用拉格朗日算法，在“針床型”波阻抗梯度飛片高速碰撞靶板時(shí)，飛片和靶板都會(huì)產(chǎn)生由于局部接觸引起局部大變形，從而導(dǎo)致局部網(wǎng)格單元發(fā)生畸變問題，使得計(jì)算提前中止。通常的解決方案是采用侵蝕算法，把嚴(yán)重畸變的網(wǎng)格單元從模型中刪除，但網(wǎng)格單元?jiǎng)h除降低了計(jì)算精度。為了解決大變形網(wǎng)格畸變問題，人們常采用歐拉算法或者光滑粒子流體動(dòng)力學(xué)（SPH）算法。相比拉格朗日算法，雖然歐拉算法解決了網(wǎng)格畸變問題，但更加耗費(fèi)機(jī)時(shí)。本文中擬采用SPH算法來計(jì)算“針床型”波阻抗梯度飛片高速碰撞靶板問題。SPH算法是將連續(xù)的固體離散為一系列具有質(zhì)量、速度和能量的粒子，通過求解質(zhì)點(diǎn)組的動(dòng)力學(xué)方程和跟蹤每個(gè)質(zhì)點(diǎn)的運(yùn)動(dòng)軌道求得整個(gè)系統(tǒng)的力學(xué)行為。該方法是一種無網(wǎng)格拉格朗日數(shù)值計(jì)算方法，廣泛應(yīng)用于處理大變形問題。

“針床型”波阻抗梯度飛片上有非常多的小正四棱錐臺(tái)，周期性排列。小正四棱錐臺(tái)頂面尺寸非常小，為了保證足夠的計(jì)算精度，建模時(shí)需要?jiǎng)澐衷S多尺寸非常小粒子，整個(gè)模型的粒子數(shù)龐大，計(jì)算耗費(fèi)時(shí)間非常長。為此，這里采用了簡化模型，即選擇一個(gè)典型單胞來計(jì)算，單胞SPH計(jì)算模型如圖2所示，通過在x軸和y軸方向設(shè)置固定邊界條件來保證僅僅z軸方向發(fā)生應(yīng)變，以滿足一維應(yīng)變加載條件。典型計(jì)算單胞的邊長為2 mm，基座高2 mm，正四棱錐臺(tái)的底邊寬1 mm、頂邊寬0.2 mm、高6 mm，靶板邊長2 mm、高7 mm。模型中SPH粒子的直徑為100 μm。由于高速碰撞的物理過程，目前計(jì)算中常采用Johnson-Cook本構(gòu)模型，具體包含線性熱彈性、Mises屈服準(zhǔn)則、各向同性應(yīng)變硬化、應(yīng)變率硬化、各向同性損傷引起的軟化，其等效Mises屈服應(yīng)力表示如下：

式中：A、B、n、C和m是模型常數(shù)； εeq是等效塑性應(yīng)變是無量綱化等效塑性應(yīng)變率，為參考應(yīng)變率；T?=(T-T0)/(Tm-T0)為無量綱溫度，其中T0是室溫，Tm是材料的熔點(diǎn)。計(jì)算時(shí)波阻抗梯度飛片和靶板材料均為鋼，表1給出了數(shù)值模擬中采用的材料本構(gòu)模型參數(shù)。

圖2 波阻抗梯度飛片氣炮實(shí)驗(yàn)的SPH單胞計(jì)算模型Fig.2 SPH cell simulation model of impedance gradient flyer in gas gun experiments

表1 飛片和靶板材料本構(gòu)和狀態(tài)方程參數(shù)Table 1 Flyer and target material parameters

3 計(jì)算結(jié)果討論與分析

3.1 平面波產(chǎn)生和傳播特性分析

圖3 靶板內(nèi)不同時(shí)刻的應(yīng)力波傳播云圖Fig.3 Stress wave propagation contours in the specimen at different times

首先通過一個(gè)算例來討論一下“針床型”梯度飛片高速撞擊靶板后靶板中應(yīng)力波傳播特性。圖3給出了撞擊速度350 m/s下靶板中沿z軸方向傳播的不同時(shí)刻的壓力云圖。初始階段可以清楚地看到每一個(gè)小四棱錐臺(tái)接觸處產(chǎn)生了球面加載波，彼此獨(dú)立往前傳播；隨著傳播距離增大，球面加載波開始相互作用；隨著傳播距離進(jìn)一步增大，最后可以看到開始的球面波逐漸轉(zhuǎn)變?yōu)槠矫娌?，從而?shí)現(xiàn)對(duì)靶板平面加載。在實(shí)際應(yīng)用中，為了保證靶板承受的是平面加載，需要在靶板前增加背板。

3.2 撞擊速度和幾何結(jié)構(gòu)參數(shù)的影響

3.2.1 撞擊速度

首先討論飛片撞擊速度對(duì)準(zhǔn)等熵壓縮加載特性的影響。計(jì)算中，保持典型單胞幾何結(jié)構(gòu)參數(shù)不變，僅僅改變飛片撞擊速度。圖4給出了350 m/s速度下不同時(shí)刻飛片和靶板的變形過程和塑性應(yīng)變?cè)茍D，從圖中可以看到：(1) 飛片上的小四棱錐臺(tái)與靶板的接觸處首先出現(xiàn)了比較嚴(yán)重的塑性變形，塑性變形區(qū)隨著時(shí)間逐漸增大；(2) 小四棱錐臺(tái)之間間隙不斷被由于嚴(yán)重塑性變形而破損的物質(zhì)填充，小四棱錐臺(tái)的高度不斷減小，其整體壓縮量超過50%。圖5給出了200、350、500、800 m/s撞擊條件下靶板自由面速度時(shí)程曲線，從圖中看到：(1) 靶板中產(chǎn)生了緩慢上升的加載波形，在初始階段自由面速度波形出現(xiàn)了一些振蕩，但隨著速度增大，振蕩現(xiàn)象逐漸轉(zhuǎn)弱，速度時(shí)程曲線趨于光滑，同時(shí)還觀察到初始的振蕩現(xiàn)象還依賴于飛片撞擊速度，高速碰撞時(shí)振蕩現(xiàn)象比較微弱，且持續(xù)時(shí)間比較短；(2) 加載波形的上升沿時(shí)間隨著飛片撞擊速度的增大而不斷減小。自由面速度初始振蕩現(xiàn)象和上升沿時(shí)間對(duì)撞擊速度的依賴性主要?dú)w因于飛片的塑性變形過程及塑性波傳播。圖6給出了不同撞擊速度下飛片最終變形圖，隨著撞擊速度增大，小四棱錐臺(tái)壓縮量不斷增加。

圖7給出了自由面速度波形的上升沿時(shí)間隨著撞擊速度的變化曲線，結(jié)果顯示兩者之間呈指數(shù)衰減關(guān)系，擬合方程為：

圖4 不同時(shí)刻飛片和靶板的塑性應(yīng)變?cè)茍DFig.4 Plastic strain contours in the flyer and specimen at different times

式中：t為自由面速度波形的上升沿時(shí)間，v0為飛片的撞擊速度。另外，靶板自由面速度峰值也依賴于飛片撞擊速度。圖8給出了自由面速度的峰值與撞擊速度的比值隨撞擊速度的變化曲線，隨著撞擊速度增大，自由面速度的峰值更接近于撞擊速度，比值趨近于1，兩者之間呈指數(shù)增長關(guān)系，擬合方程為：

式中：vmax/v0為靶板自由面速度峰值與飛片撞擊速度比，v0為飛片的撞擊速度。

圖5 不同撞擊速度下靶板自由面速度時(shí)程曲線Fig.5 Free surface velocity profiles of specimen at different impact velocities

圖6 不同撞擊速度下飛片變形的最終形態(tài)Fig.6 Final deformation of the flyer at different impact velocities

圖7 壓縮波上升沿時(shí)間隨著撞擊速度的變化曲線Fig.7 Compression wave front time vs.impact velocities

圖8 自由面速度峰值與撞擊速度比值隨著撞擊速度的變化曲線Fig.8 Ratios between peak velocity and impact velocity vs.impact velocity

3.2.2 幾何結(jié)構(gòu)參數(shù)

除了上述的飛片撞擊速度，“針床型”飛片的幾何結(jié)構(gòu)參數(shù)也對(duì)準(zhǔn)等熵壓縮加載特性有顯著的影響。首先討論小四棱錐臺(tái)高度的影響。其他幾何參數(shù)不變條件下，四棱錐臺(tái)高(h)分別選擇6、8、10 mm，以相同撞擊速度(350 m/s)撞擊靶板。圖9給出了靶板自由表面速度時(shí)程曲線，從中看到加載波形的上升沿時(shí)間隨著錐臺(tái)高度增加而增大，但速度峰值隨著錐臺(tái)高度增加而減小。接著，討論小四棱錐臺(tái)錐角的影響，這里定義正四棱錐傾斜面與底面夾角稱為錐角，如圖10所示，保持小四棱錐臺(tái)上底邊寬度不變，通過改變下底邊邊長來實(shí)現(xiàn)。下底邊邊長分別取 1、1.5、2和 3 mm，對(duì)應(yīng)錐角分別為 86.2°、83.8°、81.4°、76.9°。圖11給出了相同撞擊速度(350 m/s)不同錐角條件下靶板自由面速度時(shí)程曲線的對(duì)比，圖中顯示：上升沿時(shí)間隨著錐角增大而減小，但峰值速度隨著錐角增大而增大。

圖9 不同四棱錐臺(tái)高條件下靶板自由面速度時(shí)程曲線對(duì)比Fig.9 Comparison of the free surface velocity profiles for the specimen with different heights of the pyramid

圖12(a)～(b)分別給出了上升沿時(shí)間和峰值速度隨錐角的變化關(guān)系曲線，上升沿時(shí)間隨錐角增大呈線性減小關(guān)系，擬合關(guān)系式為：

圖10 四棱錐臺(tái)的錐角定義示意圖Fig.10 Schematic diagram of the cone angle of the pyramid

圖11 不同錐角條件下靶板自由面速度時(shí)程曲線對(duì)比Fig.11 Comparison of free surface velocity profiles for the pyramid specimen with different cone angles

圖12 壓縮波上升沿時(shí)間和峰值速度隨著正四棱錐臺(tái)錐角的變化曲線Fig.12 Compression wave front time and peak velocity vs.the cone angle of the pyramid

而峰值速度與錐角之間呈線性增長關(guān)系，擬合關(guān)系式為：

式中：t表示波陣面前沿上升沿時(shí)間，μs；β表示梯度飛片正四棱錐傾斜面與底面夾角形成錐角；vmax為靶板自由面速度峰值，m/s。

最后討論正四棱錐臺(tái)的上底邊寬的影響。保持下底邊和錐臺(tái)高度不變，上底邊寬分別為0.2、0.4和0.6 mm，以相同初始速度（350 m/s）沖擊靶板。圖13給出了不同上底邊寬帶條件下的靶板自由面速度的時(shí)程曲線，從圖中可知：自由面速度曲線初始階段出現(xiàn)了明顯的臺(tái)階，隨著上底邊寬度增大，臺(tái)階高度也增大，為了獲得比較平順的準(zhǔn)等熵壓縮加載條件，四棱錐臺(tái)上底邊寬度不易過大；同時(shí)上底邊寬度也影響著上升沿時(shí)間，寬度增大，上升沿時(shí)間變小。

圖13 不同上底邊寬度條件下靶板自由面速度時(shí)程曲線對(duì)比Fig.13 Comparison of free surface velocity profiles for the specimen with different widths of the upper edge

4 實(shí) 驗(yàn)

4.1 廣義波阻抗梯度飛片制備

利用Pro-E軟件建立了廣義波阻抗梯度飛片三維模型，如圖14(a)所示，飛片直徑53 mm，厚8 mm，基座高2 mm，四棱錐臺(tái)高6 mm，正四棱錐臺(tái)下底邊寬2 mm。采用德國EOS—M280激光選區(qū)熔化金屬3D打印機(jī)進(jìn)行制備，粉末材料采用GP1不銹鋼粉末，制備完成的飛片如圖14(b)所示。

圖14 廣義波阻抗梯度飛片的3D幾何模型和樣品Fig.14 3D model and product of generalized wave impedance gradient flyer

4.2 氣炮實(shí)驗(yàn)

實(shí)驗(yàn)在一級(jí)輕氣炮上進(jìn)行，實(shí)驗(yàn)裝置示意圖如圖15所示，為了保證擊靶平面性，飛片安裝于鋁合金彈托上，通過高壓氣體驅(qū)動(dòng)彈托和飛片高速撞擊靶板，飛片擊靶速度由測(cè)速探針來測(cè)量，采用全光纖激光位移干涉系統(tǒng)（DISAR）來監(jiān)測(cè)靶板自由面速度時(shí)程曲線。實(shí)驗(yàn)中測(cè)得飛片撞擊速度為348 m/s。圖16給出了實(shí)測(cè)的靶板自由面速度時(shí)程曲線。按照實(shí)驗(yàn)條件開展了相關(guān)的數(shù)值模擬工作，計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果的對(duì)比也如圖16所示，兩者基本一致，從而驗(yàn)證了設(shè)計(jì)思想和數(shù)值計(jì)算的有效性。

圖15 輕氣炮實(shí)驗(yàn)裝置示意圖Fig.15 Schematic of the gas gun device

圖16 實(shí)驗(yàn)與數(shù)值計(jì)算得到靶板自由面速度時(shí)程曲線對(duì)比Fig.16 Experimental and numerical free surface velocities of the target

5 結(jié) 論

(1) “針床型”小四棱錐臺(tái)高速碰撞靶板后，在靶板中產(chǎn)生密集分布的球面加載波，隨著傳播距離的增大，球面加載波相互疊加形成緩慢上升前沿的平面加載波，從而實(shí)現(xiàn)對(duì)靶板的準(zhǔn)等熵壓縮加載；

(2) 準(zhǔn)等熵壓縮加載的特性強(qiáng)烈依賴于飛片擊靶的速度，隨著速度的增加，加載波上升沿時(shí)間減小，而峰值速度增大；

(3) 廣義波阻抗飛片的關(guān)鍵幾何參數(shù)也對(duì)準(zhǔn)等熵壓縮加載特性影響顯著，隨著四棱錐臺(tái)高增大，加載波上升沿時(shí)間增長，而峰值速度減小；相反，保持四棱錐臺(tái)高不變，只改變錐角的大小，隨著錐角的增大加載波上升沿時(shí)間縮短，而峰值速度增大；小四棱錐臺(tái)變鈍（即上底邊寬度增大）會(huì)導(dǎo)致加載波初始階段出現(xiàn)明顯的臺(tái)階，為了獲得比較平順的準(zhǔn)等熵壓縮加載條件，四棱錐臺(tái)需要保持尖銳度。