適用于自洽強度方法的沖擊加載-再加載實驗技術*

俞宇穎,譚 葉,譚 華,戴誠達,彭建祥,李雪梅,吳 強,王 翔

(1.中國工程物理研究院流體物理研究所沖擊波物理與爆轟物理重點實驗室，四川 綿陽 621999；2.中國工程物理研究院壓縮科學研究中心,四川 綿陽 621999)

適用于自洽強度方法的沖擊加載-再加載實驗技術*

俞宇穎1,2,譚 葉1,2,譚 華1,2,戴誠達1,2,彭建祥1,2,李雪梅1,2,吳 強1,2,王 翔1,2

(1.中國工程物理研究院流體物理研究所沖擊波物理與爆轟物理重點實驗室，四川 綿陽 621999；2.中國工程物理研究院壓縮科學研究中心,四川 綿陽 621999)

針對自洽強度方法存在的沖擊加載-再加載的難題，提出了一種采用較高硬度材料為支撐制作組合飛片的簡便方法。利用該方法獲得了鋁、錫和鋯基金屬玻璃較理想的沖擊加載-再加載粒子速度剖面，驗證了該方法的有效性。由本文獲得的沖擊加載-再加載粒子速度剖面，并根據自洽方法，計算得到了鋁、錫和鋯基金屬玻璃再加載過程剪應力變化數據。進一步分析表明，在本文涉及的壓力范圍內，僅由沖擊加載-卸載實驗得到的鋁、錫和鋯基金屬玻璃屈服強度將比實際結果降低20%～50%。因此，在采用自洽方法計算高壓強度時，沖擊加載-再加載數據不可或缺。

固體力學；高壓強度；沖擊加載-再加載；自洽方法；組合飛片

高壓下材料屈服強度特性是當前爆炸力學、沖擊動力學、高壓物理和材料科學等學科領域重點關注的基礎問題，研究成果在裝甲防護、航空航天器防護等方面有重要的應用價值。自20世紀60年代以來，先后發展了靜水壓線比較法[1]、壓剪方法[2]、橫向應力計方法[3]和自洽方法[4](也稱為雙屈服面方法或AC方法)等高壓強度研究方法。靜水壓線比較法受狀態方程精度影響，在高壓下精度并不高。壓剪方法和橫向應力計方法分別受到加載技術和測試技術的限制，往往局限于20 GPa以內。自洽方法原則上沒有壓力限制，但需要進行一組沖擊加載-卸載實驗和沖擊加載-再加載實驗才能獲取一個強度數據。其中沖擊加載-再加載實驗相對復雜，需要采用由高、低阻抗材料構成的組合飛片進行沖擊才能實現。由于組合飛片在氣炮驅動的過程中容易發生分離，導致在待測材料中形成沖擊加載-卸載-再沖擊加載，而非預期的沖擊加載-再加載，也就無法獲得強度數據。在現有的高壓強度研究方法中，自洽方法是目前應用最廣泛的一種方法，提供了包括金屬材料(如鋁及合金[4-10]、銅[6]、釩[11]、鎢[12-14]、鉭[15]、鈹[16]、鋼[17-18]等)、陶瓷材料(B4C[19]、SiC[20])、金屬玻璃(Zr51Ti5Ni10Cu25Al9)[21]和復合材料[22]等大量材料的高壓強度數據。但正如前面所指出的，由于沖擊加載-再加載實驗難度大，目前自洽方法提供的高壓強度數據大部分僅來自沖擊加載-卸載實驗，這類強度數據實際上并不完備。

針對上述自洽方法的問題，胡建波等[9-10]通過采用平面焊接技術加強組合飛片之間的結合力來避免其分離。為了簡化波系作用，利于數據處理分析，組合飛片的第一層通常為待測樣品材料。胡建波等[9-10]利用由此制作的TC4/LY12鋁組合飛片成功地加速到3.2 km/s，并對LY12鋁樣品進行了對稱碰撞，獲得了34 GPa壓力下(指第一次沖擊壓力，下同)LY12鋁的沖擊加載-再加載粒子速度剖面。在此之前，鋁的沖擊加載-再加載實驗的最高壓力僅為22 GPa[5]。類似地，M.D.Furnish等[18]則采用爆炸焊接技術制作了Ta/2169鋼組合飛片。該飛片被成功地加速到2.0 km/s，并由此獲得了48 GPa壓力下2169鋼的沖擊加載-再加載粒子速度剖面。上述焊接技術的應用極大提升了沖擊加載-再加載實驗能力，拓展了自洽方法的實驗壓力范圍。然而，上述焊接技術比較復雜，更重要的是焊接處理時會經歷高溫、高壓作用過程，可能改變組合飛片中的待測樣品材料初始狀態，從而對高壓強度結果造成影響。值得注意的是，大量的輕氣炮實驗數據表明：在輕氣炮驅動飛片加速過程中，飛片會發生如弓形或馬鞍形等復雜的彎曲變形，而且不同材料飛片的變形量及形狀也不同[23-24]。事實上，上述因素極易造成組合飛片在氣炮驅動的過程中發生分離。

基于上述分析，本文中提出一種實現沖擊加載-再加載實驗的簡便方法—采用較高硬度材料為支撐，通過環氧樹脂與待測樣品粘接制成組合飛片，由此減小氣炮加載下飛片的彎曲變形來避免分離，并通過鋁、錫和鋯基金屬玻璃的沖擊實驗驗證該方法的有效性。在此基礎上，獲得鋁、錫和鋯基金屬玻璃再加載過程剪應力的變化情況，并初步分析其對強度結果的影響。

1 實驗方法

實驗采用如圖1所示的反向碰撞方式，即由組合飛片直接撞擊透明的單晶LiF窗口。組合飛片的第一層為待測材料(即樣品)，其后為用于支撐樣品的具有較高硬度的材料，兩層飛片之間采用環氧樹脂進行粘接，環氧樹脂厚度約10 μm。上述組合飛片的制作過程均在常溫、常壓下進行，不會改變組合飛片中的待測樣品材料初始狀態。LiF窗口尺寸為?28 mm×12 mm，碰撞面鍍有1 μm厚的鋁膜作為光學測試的反射面，同時為保護長歷時測量過程中鋁膜不受破壞，鋁膜前粘接了8 μm銅箔。一種激光干涉測速技術—DISAR(displacement interferometer system for any reflector)技術[25]用于測量樣品/LiF窗口界面粒子速度剖面，通過剖面可以直觀判斷樣品中是否經歷預期的沖擊加載-再加載過程。組合飛片沖擊LiF窗口的速度由磁測速技術進行測量。根據實測的飛片速度以及樣品和窗口材料的Hugoniot參數，采用阻抗匹配法可以計算得到樣品的沖擊壓力和沖擊波速度，結合圖1所示的波系作用可以進一步計算得到再加載彈塑性波的聲速cL。

組合飛片中支撐材料除具有較高硬度外，其阻抗應高于樣品材料的阻抗，以實現對樣品的沖擊加載-再加載。當然，為確保再加載過程出現彈塑性特征信息，再加載壓力幅度不能過高，以避免形成沖擊波。根據上述條件，針對LY12鋁樣品采用了TC4鈦合(成分為Ti90-Al6-V4)作為支撐材料，錫和鋯基金屬玻璃(Zr51Ti5Ni10Cu25Al9，詳細參數參見文獻[26])則采用45鋼作為支撐材料。

圖1 實驗裝置及波系作用示意圖Fig.1 Schematic diagrams of experimental setup and wave interactions

2 實驗結果與初步分析

利用圖1所示的實驗裝置，在?30 mm二級輕氣炮上進行了5發鋁、錫和鋯基金屬玻璃為樣品的沖擊加載-再加載實驗，沖擊速度范圍為2.49～4.39 km/s，在樣品中產生的壓力范圍為28～48 GPa，樣品和支撐材料的直徑均為28 mm,詳細實驗參數見表1,表中Hs為樣品厚度，hs為支撐材料厚度，D為沖擊波速度，p為沖擊壓力，τH和τc分別為Hugoniot狀態剪應力和臨界剪應力。為了盡量減小氣炮加載下組合飛片的彎曲變形，支撐材料總厚度不小于3 mm。尤其是在最高沖擊速度為4.39 km/s的實驗中(實驗2)，除2 mm厚的TC4外另增加了3 mm厚的鉭作為支撐材料。

表1 平板沖擊實驗條件及結果Table 1 Experimental conditions and results for planar plate-impact experiments

圖2 LY12鋁樣品/窗口界面粒子速度剖面Fig.2 Particle velocity profiles measured at LY12 aluminum sample/window interface

圖3 錫樣品/窗口界面粒子速度剖面Fig.3 Particle velocity profiles measured at Sn sample/window interface

圖4 鋯基金屬玻璃樣品/窗口界面粒子速度剖面Fig.4 Particle velocity profile measured at Zr-based bulk metallic glass sample/window interface

由DISAR測得的鋁、錫和鋯基金屬玻璃為樣品的沖擊加載-再加載粒子速度剖面如圖2～4所示。從圖中可以看到，5發實驗均獲得了預期的沖擊加載-再加載波剖面，這表明采用較高硬度材料為支撐以減小組合飛片彎曲變形來避免分離的方法是簡便、有效的。此外，鋁、錫和鋯基金屬玻璃為樣品的再加載波剖面均出現一定程度的彈塑性波剖面特征，表明相對應支撐材料的阻抗選取是合適的。圖中uw,H和uw,1分別為再加載彈性波和塑性波起始對應的粒子速度。

圖5 LY12鋁從沖擊壓縮態再加載和卸載過程的拉格朗日縱波和體波聲速隨粒子速度的變化Fig.5 Longitudinal and bulk Lagrangian wave speeds during reloading and unloading from shocked state of LY12 aluminum

根據圖1所示的波傳播特性，由圖2～4的粒子速度剖面可以計算得到沿著再加載過程的拉格朗日縱波聲速：

(1)

式中：Ds為碰撞時樣品中產生的沖擊波的速度(由實測的飛片速度以及樣品和窗口材料的Hugoniot參數，通過阻抗匹配法計算得到沖擊波速度)，t為如圖1所示的來自樣品/支撐材料界面的再加載波到達樣品/窗口界面時間與碰撞時間之差。

圖5給出了LY12鋁再加載過程拉格朗日縱波聲速隨粒子速度u變化的典型結果(實驗1)。其中，粒子速度u由樣品/窗口界面粒子速度uw結合增量型阻抗匹配法[12,27]計算得到，由此得到的粒子速度計及了加載波在樣品/窗口界面反射造成的影響。uH和u1分別為Hugoniot狀態對應的粒子速度和再加載進入塑性屈服時對應的粒子速度，由圖2對應的uw,H和uw,1根據增量型阻抗匹配法計算得到。圖5中同時給出了由LY12鋁沖擊加載-卸載粒子速度剖面計算得到的卸載過程拉格朗日縱波聲速的結果[28]。從圖中可以看到，再加載過程也呈現與卸載過程相同的準彈性行為特征，即彈、塑性波速為光滑過渡，而沒有發生突降。如圖中虛線所示，將卸載過程塑性段聲速線性外延可得到再加載過程的拉格朗日體波聲速cb。

在獲得縱波和體波的基礎上，根據J.R.Asay等[4]提出的自洽強度方法，對再加載過程的聲速進行如下計算可得到剪應力的變化：

(2)

式中:uH和u1分別為圖5所示的Hugoniot狀態對應的粒子速度和再加載進入塑性屈服時對應的粒子速度，ρ0為材料的初始密度，cL和cb分別為圖5所示的再加載過程對應的拉格朗日縱波和體波聲速。由此計算得到的LY12鋁的τc-τH結果見表1。結合胡建波等[29]和俞宇穎等[21]分別給出的錫和鋯基金屬玻璃卸載過程聲速數據，采用上述相同方法計算得到了錫和鋯基金屬玻璃再加載過程的τc-τH，計算結果見表1。

在獲得τc-τH的基礎上，如能獲得相同沖擊壓力下卸載過程的聲速則可計算得到τc+τH，進而可以確定在該沖擊壓力下材料的屈服強度Y=(τc-τH)+(τc+τH)。由于沖擊加載-再加載實驗的困難，先前大量文獻中將屈服強度Y≈τc+τH，即認為再加載過程的τc-τH可以忽略。就LY12鋁而言，本文再加載實驗得到的38.3和48.5 GPa沖擊壓力下的τc-τH分別為0.73和0.77 GPa；胡建波等[10]通過卸載實驗得到的32.2和54.7 GPa沖擊壓力下的τc+τH分別為0.85和1.05 GPa。顯然，在上述壓力范圍內僅由沖擊加載-卸載實驗得到的LY12鋁強度Y≈τc+τH將比實際結果Y=2τc降低約45%。對于錫，本文再加載實驗得到的28.7和38.1 GPa沖擊壓力下的τc-τH分別為0.08和0.19 GPa；依據胡建波等[29]的卸載實驗剖面數據計算得到的28.3和39.2 GPa沖擊壓力下的τc+τH分別為0.07和0.17 GPa。僅由沖擊加載-卸載實驗得到的錫強度Y≈τc+τH將比實際結果Y=2τc降低約50%。同樣地，就鋯基金屬玻璃而言，本文再加載實驗得到的39.1 GPa沖擊壓力下的τc-τH為0.53 GPa，俞宇穎等[21]通過卸載實驗得到的37.3 GPa沖擊壓力下的τc+τH為1.73 GPa，僅由沖擊加載-卸載實驗得到的強度Y≈τc+τH將比實際結果Y=2τc降低約20%。綜上所述，在采用自洽方法計算高壓強度時沖擊加載-再加載數據不可或缺。

3 結 論

針對自洽高壓強度方法存在的因組合飛片分離而難以實現沖擊加載-再加載的難題，提出了一種簡便方法—采用較高硬度材料為支撐，通過環氧樹脂與待測樣品粘接制成組合飛片，由此減小氣炮加載下飛片的彎曲變形來避免分離。采用該方法，在二級輕氣炮上進行了沖擊速度2.49～4.39 km/s、沖擊壓力28～48 GPa范圍內鋁、錫和鋯基金屬玻璃為樣品的驗證實驗，5發實驗均獲得了較理想的沖擊加載-再加載粒子速度剖面，表明了該方法的有效性。由本文獲得的沖擊加載-再加載粒子速度剖面，根據自洽方法計算得到了鋁、錫和鋯基金屬玻璃再加載過程剪應力變化數據。結合已有的沖擊加載-卸載過程剪應力變化數據的分析表明，在本文涉及的壓力范圍內，僅由沖擊加載-卸載實驗得到的鋁、錫和鋯基金屬玻璃屈服強度將比實際結果降低約20%～50%。因此，在采用自洽方法計算高壓強度時沖擊加載-再加載數據不可或缺。

[1] Fowles G R. Shock wave compression of hardened and annealed 2024 aluminum[J]. Journal of Applied Physics, 1961,32(8):1475-1487.

[2] Clifton R J, Klopp R W. Pressure-shear plate impact testing[M]∥Metals Handbook: Mechanical Testing. OH: American Society for Metals, 1985:230-239.

[3] Rosenberg Z, Partom Y, Yaziv D. The use of in-material stress gauges for estimating the dynamic yield strength of shock-loaded solids[J]. Journal of Applied Physics, 1984,56(1):143-146.

[4] Asay J R, Chhabildas L C. Determination of the shear strength of shock compressed 6061-T6 aluminum[C]∥Meyers M M, Murr L E. Shock Waves and High-strain-rate Phenomena in Metals: Concepts and Application. New York: Plenum Publishing Corp., 1981:417-431.

[5] Asay J R, Chhabildas L C, Kerley G I, et al. High pressure strength of shocked aluminum[M]∥Gupta Y M. Shock Waves in Condensed Matter. New York: Plenum Press,1986:145-149.

[6] Morris C E, Fritz J N, Holian B L. Quasi-elastic high pressure waves in 2024 Al and copper[C]∥Nellis W J. Shock Waves in Condensed Matter. AIP, 1982:382-386.

[7] Huang H, AsayJ R. Compressive strength measurements in aluminum for shock compression of the stress range of 4-22 GPa[J]. Journal of Applied Physics, 2005,98(3):033524.

[8] Huang H, Asay J R.Reshock and release response of aluminum single crystal[J]. Journal of Applied Physics, 2007,101(6):063550.

[9] 胡建波,戴誠達,俞宇穎,等.雙屈服面法測量金屬材料動高壓屈服強度的若干改進[J].爆炸與沖擊,2006,26(6):516-521. Hu Jianbo, Dai Chengdai, Yu Yuying, et al. Some improvements of the self-consistent method for measuring the dynamic yield strength of ductile metals[J]. Explosion and Shock Waves, 2006,26(6):516-521.

[10] 胡建波,譚華,俞宇穎,等.鋁的動態屈服強度測量及再加載彈性前驅波的形成機理分析[J].物理學報,2008,57(1):405- 410. Hu Jianbo, Tan Hua, Yu Yuying, et al. Measurements of dynamic yield strength of aluminum alloy and mechanism analysis of elastic precursor during reloading[J]. Acta Physica Sinica, 2008,57(1):405- 410.

[11] Chhabildas L C, Hills C R. Dynamic shock studies of vanadium[C]∥Murr L E. Metallurgical Applications of Shock-wave and High-strain-rate Phenomena. New York: Marcel Dekker, 1985:429-448.

[12] Asay J R, Chhabildas L C, Danderkar D P. Shear strength of shock-loaded polycrystalline tungsten[J]. Journal of Applied Physics, 1980,51(9):4774-4783.

[13] 華勁松.高溫高壓下鎢合金的本構方程研究[D].北京:中國工程物理研究院北京研究生部,1999:1-118.

[14] 張江躍,譚華,虞吉林.雙屈服法測定93W合金的屈服強度[J].高壓物理學報,1997,11(4):254-259. Zhang Jiangyue, Tan Hua, Yu Jilin. Determination of the yield strength of 93W alloys by using AC techniques[J]. Chinese Journal of High Pressure Physics, 1997,11(4):254-259.

[15] Chhabildas L C, Barker L M, Asay J R, et al. Relationship of fragment size to normalized spall strength for materials[J]. International Journal of Impact Engineering, 1990,10(1):107-124.

[16] Chhabildas L C, Wise J L, Asay J R. Reshock and release behavior of beryllium[C]∥Nellis W J. Shock Waves in Condensed Matter. AIP, 1982:422-426.

[17] Duffy T S, Ahrens T J. Dynamic compression of an Fe-Cr-Ni alloy to 80 GPa[J]. Journal of Applied Physics, 1997,82(9):4259-4269.

[18] Furnish M D, Alexander C S, Brown J L, et al. 2169 steel waveform measurements for equation of state and strength determination[J]. Journal of Applied Physics, 2014,115(3):033511.

[19] Vogler T J, Reinhart W D, Chhabildas L C. Dynamic behavior of boron carbide[J]. Journal of Applied Physics, 2004,95(8):4173-4183.

[20] Vogler T J, Reinhart W D, Chhabildas L C, et al. Hugoniot and strength behavior of silicon carbide[J]. Journal of Applied Physics, 2006,99(2):023512.

[21] 俞宇穎,習鋒,戴誠達,等.動高壓加載下鋯基金屬玻璃強度測量[J].爆炸與沖擊,2014,34(1):1-5. Yu Yuying, Xi Feng, Dai Chengdai, et al. Measurement of strength in a Zr-based bulk metallic glass under dynamic high-pressure loading[J]. Explosion and Shock Waves, 2014,34(1):1-5.

[22] Yuan F P, Tsai L, Prakash V, et al. Dynamic shear strength of S2 glass fiber reinforced polymer composites under shock compression[J]. Journal of Applied Physics, 2008,103(10),103537.

[23] 王翔.金屬材料狀態方程精密實驗測量技術研究[D].北京:中國工程物理研究院北京研究生部,2004:1-94.

[24] Mitchell A C, Nellis W J. Shock compression of aluminum, copper and tantalum[J]. Journal of Applied Physics, 1981,52(5):3363-3374.

[25] Weng Jidong, Tan Hua, Wang Xiang, et al. Optical-fiber interferometer for velocity measurements with picosecond resolution[J]. Appled Physics Letters, 2006,89(11):111101.

[26] 俞宇穎,習鋒,戴誠達,等.沖擊加載下Zr51Ti5Ni10Cu25Al9金屬玻璃的塑性行為[J].物理學報,2012,61(19):382-387. Yu Yuying, Xi Feng, Dai Chengda, et al. Plastic behavior of Zr51Ti5Ni10Cu25Al9metallic glass under planar shock loading[J]. Acta Physica Sinica, 2012,61(19):382-387.

[27] 譚華.實驗沖擊波物理導引[M].北京:國防工業出版社,2007:163-167.

[28] Yu Yuying, Tan Hua, Hu Jianbo, et al. Determination of effective shear modulus of shock-compressed LY12 Al from particle velocity profile measurements[J]. Journal of Applied Physics, 2008,103(10):103529.

[29] Hu Jianbo, Zhou Xianming, Dai Chengda, et al. Shock-induced bct-bcc transition and melting of tin identified by sound velocity measurements[J]. Journal of Applied Physics, 2008,104(8):083520.

(責任編輯 張凌云)

A shock-reload wave technique for dynamic strength study of materials at high pressure by self-consistent method

Yu Yuying1,2, Tan Ye1,2, Tan Hua1,2, Dai Chengda1,2, Peng Jianxiang1,2,Li Xuemei1,2, Wu Qiang1,2, Wang Xiang1,2

(1.LaboratoryforShockWaveandDetonationPhysicsResearch,InstituteofFluidPhysics,ChinaAcademyofEngineeringPhysics,Mianyang621999,Sichuan,China；2.CenterforCompressionScience,ChinaAcademyofEngineeringPhysics,Mianyang621999,Sichuan,China)

A convenient method for fabricating a layered impactor with a sample backed up by high hardness materials was developed to overcome the obstacle of the shock-reload experiments in the self-consistent yield strength technique. This method was validated by a series of ideal shock-reload particle velocities of aluminum, tin, and Zr-based bulk metallic glass obtained from the reverse-impact experiments at the peak shock stresses from 28 GPa to 48 GPa. The sum of the shear stresses for these materials in the reload process from the shocked state was estimated, and compared with the previously reported data in the release process. It is shown that the yield strength under a high pressure for the materials investigated will underestimate 20%-50% if without the reload data. Thus, shock-reload experiments are essential for the self-consistent yield strength method.

solid mechanics; high-pressure yield strength; shock-reload wave technique; self-consistent method; layered impactor

10.11883/1001-1455(2016)04-0491-06

2014-11-28；

2015-05-27

國家自然科學基金項目(10972206,11172281)

俞宇穎(1976— )，男，博士，副研究員,yuyinyu@caep.cn。

O347;O521.2國標學科代碼：13015

A