趙鵬鐸，盧芳云，陳 榮，李俊玲，林玉亮，譚多望

(1.國防科學(xué)技術(shù)大學(xué)理學(xué)院，湖南 長沙410073;2.中國工程物理研究院流體物理研究所沖擊波物理與爆轟物理國防科技重點實驗室，四川 綿陽621900)

在各類涉及爆炸和沖擊加載的領(lǐng)域中，材料的動態(tài)本構(gòu)關(guān)系始終是重要的研究課題。動態(tài)本構(gòu)關(guān)系涉及應(yīng)力張量、應(yīng)變張量，應(yīng)變率以及溫度等變量之間的關(guān)系。傳統(tǒng)的一維波加載實驗只能確定部分本構(gòu)參數(shù)，動態(tài)壓剪復(fù)合加載實驗則可以為本構(gòu)關(guān)系提供更加完整的實驗參數(shù)。除此之外，材料的屈服、損傷演化、失效、相變等均與剪切密切相關(guān)。材料動態(tài)壓-剪加載的研究不僅豐富了材料的動態(tài)力學(xué)現(xiàn)象，拓展了研究范圍，而且對于認(rèn)識材料動態(tài)響應(yīng)的演化規(guī)律和機(jī)理有重要的意義［1］。

盧芳云等［2］提出了分離式霍普金森壓剪桿實驗的思想，主要借鑒了傳統(tǒng)的SHPB 與動態(tài)壓剪復(fù)合平板沖擊加載實驗技術(shù)［3-6］，對傳統(tǒng)SHPB 進(jìn)行了拓展:入射桿端頭為楔形狀，同時使用了2 個透射桿。并討論了相對于已有的動態(tài)壓剪復(fù)合加載實驗［7-10］的優(yōu)點。

L.J.Griffiths 等［11］采用光通量法對碳纖維復(fù)合材料的動態(tài)性能進(jìn)行了研究。唐志平等［12］對光通量法用于動態(tài)力學(xué)量測試的可行性進(jìn)行了系統(tǒng)分析和研究，并與傳統(tǒng)的電阻應(yīng)變片法作了比較。結(jié)果表明:光通量法法能可靠地用于準(zhǔn)靜態(tài)直至應(yīng)變率達(dá)105s－1范圍內(nèi)質(zhì)點位移、速度，以及試件應(yīng)變、應(yīng)變率、應(yīng)力和泊松比等力學(xué)量的測量。K.T.Ramesh 等［13］利用光通量法測量了SHPB 實驗中試樣的橫向變形。李玉龍等［14］改進(jìn)了Ramesh 測試裝置，在激光器前加設(shè)了光學(xué)準(zhǔn)直部件，使激光光束成為平行光，并測試了動態(tài)拉伸試件的徑向應(yīng)變。

本文中利用半導(dǎo)體準(zhǔn)直激光器，采用光通量法，對SHPSB(split Hopkinson pressure bar)實驗技術(shù)中試樣的剪切應(yīng)變測量進(jìn)行實驗研究，并對測試系統(tǒng)進(jìn)行不確定度評定。最后對某含鋁PBX 試樣的剪切應(yīng)變測試值和文獻(xiàn)［2］給出的理論計算結(jié)果進(jìn)行對比分析。

1 實 驗

1.1 裝置及基本原理

實驗光學(xué)測試系統(tǒng)主要包括:半導(dǎo)體準(zhǔn)直激光器、光學(xué)狹縫、光電接收器以及光學(xué)擋板，如圖1 所示。其中準(zhǔn)直激光器輸出功率為20 mW，用以輸出波長為660 nm 的平行激光。激光光束與透射桿平行，激光光斑為圓形，經(jīng)過光學(xué)狹縫，光斑變?yōu)榫鶆驇?。光電接受器主要包括接受透鏡和光電傳感器。接收透鏡將激光匯聚到光電傳感器的感光面，這樣光電傳感器的輸出電壓就與激光的通光量成正比。光電傳感器的響應(yīng)頻率為2 MHz，系統(tǒng)噪音小于0.4 mV。光學(xué)擋板1 和光學(xué)擋板2 通過鋁箍，分別與入射桿和透射桿固定連接。由于光學(xué)擋板1 和光學(xué)擋板2 的遮擋，投射在光電接收器上的光斑呈矩形。實驗中，試樣分別與入射桿、透射桿固連，試樣界面的切向運動與實驗桿界面的切向運動是一致的，而光學(xué)擋板1 與入射桿固定連接，因此入射桿界面的切向運動和擋板1 的切向運動是一致的。同理，光學(xué)擋板2 與透射桿界面的切向運動一致。當(dāng)試樣受載剪切變形時，光學(xué)擋板1 與光學(xué)擋板2 出現(xiàn)切向的相對運動，于是造成光通量的變化，而與光電接收器連接的示波器記錄了由光通量變化對應(yīng)的電壓變化值。如果試樣僅沿著透射桿軸向運動，光通量將不會改變，即無變化電壓輸出。

實驗采用精密測微絲桿對實驗測試系統(tǒng)進(jìn)行標(biāo)定。即利用精密測微絲桿在垂直于光路的方向，遮擋激光光斑來記錄光斑寬度變化量以及輸出電壓值。測微絲桿的變化范圍是0 ～10 mm，變化間隔為0.1 mm。絲桿的運動造成矩形光斑寬度S 的變化，即光斑寬度的變化(ΔS)就對應(yīng)著光電傳感器的輸出電壓變化(ΔU)。圖2 給出了兩次標(biāo)定實驗結(jié)果，在兩次標(biāo)定實驗中測微絲桿的位置相距4 cm。由于激光光束為光強(qiáng)均勻的平行光，兩次標(biāo)定結(jié)果基本一致，且ΔS 與ΔU 高度線性相關(guān)

圖1 實驗裝置示意圖Fig.1 Schematics of the experiment

圖2 實驗標(biāo)定結(jié)果Fig.2 Results of the calibration

由此標(biāo)定出相對于單位電壓變化所對應(yīng)的距離變化值，即電壓-距離轉(zhuǎn)換系數(shù)k(mm/V)。假設(shè)試樣的原始厚度為l0，而光束變化值與試樣兩個界面的相對切向位移一致，則試樣剪切應(yīng)變?yōu)?/p>

剪切應(yīng)變率可以由相對位移的導(dǎo)數(shù)，即相對速度Δv 表示

剪切應(yīng)變測量系統(tǒng)的不確定度主要來源是:電壓-距離轉(zhuǎn)換系數(shù)k 標(biāo)定過程引入的不確定度，光電接收器的系統(tǒng)噪聲引入的不確定度。由圖2 計算出k 引入的標(biāo)準(zhǔn)不確定度［15］u1=0.003 14 mm/V。光電接收器的最大噪聲為0.4 mV，以均勻分布考慮［15］，則由均勻分布的標(biāo)準(zhǔn)不確定度計算方法可知，由光電接收器的系統(tǒng)噪聲引入的標(biāo)準(zhǔn)不確定度，而合成標(biāo)準(zhǔn)不確定度可以表示為

式中:ai是第i 項不確定度的靈敏度系數(shù)，其數(shù)值等于對公式(2)關(guān)于k 與ΔU 分別求偏導(dǎo)數(shù)。例如某次實驗中，試樣的厚度是5 mm，實驗采集到的峰值電壓ΔU=500 mV，則測試的剪切應(yīng)變的標(biāo)準(zhǔn)不確定度uc=0.03%，這對于剪切破壞應(yīng)變較小的脆性材料(γc＜5%)來說，測量結(jié)果的不確定度依舊是小量，這表明本文提出的剪切應(yīng)變測試方法具有較高的可靠性。

1.2 性能驗證實驗

剪切應(yīng)變測試方法是否有效，關(guān)鍵是:(1)所選擇的準(zhǔn)直激光器的光斑在不同位置是否均勻;(2)激光接收器響應(yīng)是否滿足動態(tài)實驗要求;(3)擋板與實驗桿連接是否完全固連。

問題(1)已經(jīng)由靜態(tài)標(biāo)定實驗結(jié)果回答，如圖2 所示。為了解決其余問題，設(shè)計如圖3 所示的子彈撞擊長桿驗證實驗。在實驗桿端使用鋁箍將光學(xué)擋板與桿固連，利用準(zhǔn)直激光器與光電接收器對實驗桿端面的運動時程曲線進(jìn)行記錄。實驗桿上應(yīng)變片記錄了入射與反射應(yīng)變，利用一維應(yīng)力波理論計算實驗桿端面的位移時程曲線。圖4 給出了2 種方法所測得的桿端的位移時程曲線，其中實線是通過應(yīng)變片計算的桿端位移時程曲線，虛線為激光系統(tǒng)測試的結(jié)果，圖中曲線第2 次上升是由于應(yīng)力波反射產(chǎn)生二次加載造成的位移。

圖3 實驗裝置示意圖Fig.3 Schematics of the valid experiment

圖4 桿端的位移時程曲線Fig.4 Displacement history of the end of the bar

由圖4 可知2 種方法所得的位移時程曲線一致，上升斜率相同。這表明接收器的頻響滿足動態(tài)實驗需要;擋板與實驗桿固連，未出現(xiàn)松動，表明擋板與實驗桿連接方法有效。問題(2)、(3)得到解決，因此本文中提出的剪切應(yīng)變測試技術(shù)是有效的。

2 實驗結(jié)果與討論

對密度為1.7 g/cm3的含鋁PBX 炸藥進(jìn)行SHPSB 實驗。實驗采用立方體試樣，其尺寸為13 mm×13 mm×5 mm。圖5 給出了典型的剪切應(yīng)力、應(yīng)變的原始波形，其中剪切應(yīng)力測試方法與文獻(xiàn)［2］相同。由圖5 可知剪切應(yīng)力與剪切應(yīng)變測試的原始信號的起跳時間相當(dāng)，前者的峰值時刻早于后者。圖6 是典型的剪切應(yīng)力、應(yīng)變時程曲線，由圖6 可知，2 種方法獲得的剪切應(yīng)變時程曲線不一致。

圖5 典型的實驗原始波形Fig.5 Oscilloscope signals of test

圖6 典型的剪切應(yīng)力、應(yīng)變時程曲線Fig.6 Shear stress and strain histories of the specimen

文獻(xiàn)［2］給出的剪切應(yīng)變率為

式中:l0為試樣的厚度，v1、v2分別為試樣兩端的切向質(zhì)點速度，εi、εr、εt分別為入射波、反射波和透射波產(chǎn)生的剪切應(yīng)變?yōu)榧羟胁ㄋ伲珿、ρ0分別為透射桿的剪切模量與密度。文獻(xiàn)［2］中是將一維彈性剪切應(yīng)力波中應(yīng)變與橫向速度的關(guān)系近似表示為透射剪切應(yīng)變與切向質(zhì)點速度的關(guān)系，即

在壓剪復(fù)合平板沖擊實驗中滿足公式(6)，即沒有考慮側(cè)向邊界對剪切波的影響。在長桿中，加載脈寬大于桿的直徑，側(cè)向邊界對剪切波傳播始終是有影響的。這類擾動被定義為彎曲波［16］。彎曲波的動力學(xué)方程不但要考慮動量守恒，同時要考慮動量矩守恒，因此公式(6)已不完全適用。圖7 中曲線1是利用LS-DYNA 模擬的彎曲波在細(xì)長桿中傳播的橫向速度曲線，代表了桿上某一橫截面積內(nèi)單元的平均橫向速度，曲線2 代表了該截面內(nèi)單元平均剪切應(yīng)變與剪切波速乘積結(jié)果，即公式(6)的右側(cè)項。

由圖7 可知，利用公式(6)計算的桿的橫向速度相對實際橫向速度明顯偏小。在SHPSB 透射桿上利用公式(6)建立橫向速度與剪切應(yīng)變的關(guān)系是不準(zhǔn)確的，透射桿與試樣端面的實際橫向速度大于文獻(xiàn)［2］計算結(jié)果。而試樣的剪切應(yīng)變可表示為

圖7 橫截面的橫向速度時程曲線Fig.7 Transverse velocity histories of the bar

式中:T 是加載波寬。透射桿的實際橫向速度v1＞c εt，則

公式(8)基于光通量法測量的剪切應(yīng)變值小于文獻(xiàn)［2］處理結(jié)果。由于在初始上升階段(圖7 中兩曲線的上升沿部分)，透射桿-試樣截面的橫向速度明顯小于入射桿-試樣截面的橫向速度，即v1＜v2，并且在此階段，v1與c εt相差較小，所以曲線上升階段γlaser與γcomp基本吻合，如圖6 所示。隨著v1與c εt的差值以及v1增大，γcomp曲線逐漸高于γlaser時程曲線。剪切應(yīng)變測試信號在達(dá)到峰值平臺之后，信號出現(xiàn)了上下波動，如圖5 所示。這可能是由于試樣/實驗桿界面連接出現(xiàn)了破壞，使得試樣界面的運動與實驗桿端面運動不再一致，此時波動的信號已經(jīng)不能反映試樣的剪切變形。信號出現(xiàn)波動的時刻對應(yīng)于實驗加載末段，剪切應(yīng)力已降低到1 MPa 左右，如圖6 所示，此階段已不是實驗關(guān)注的主要階段。

3 結(jié) 論

利用光通量法對SHPSB 實驗中的剪切應(yīng)變進(jìn)行了實驗測量，得到了試樣的剪切應(yīng)變與剪切應(yīng)變率。光學(xué)測試系統(tǒng)的標(biāo)定結(jié)果和性能測試實驗結(jié)果表明利用文中提出的基于光通量法的剪切應(yīng)變測量技術(shù)能夠準(zhǔn)確地測量SHPSB 試樣的剪切應(yīng)變。測量系統(tǒng)的不確定度分析結(jié)果表明該測量方法具有較高的可靠性。對比分析了基于光通量法測量的剪切應(yīng)變結(jié)果與理論近似處理結(jié)果，表明光通量法測量的剪切應(yīng)變相較應(yīng)用一維剪切波理論處理的剪切應(yīng)變值更加準(zhǔn)確。

［1］ 唐志平.壓剪復(fù)合平板沖擊加載技術(shù)進(jìn)展及其應(yīng)用［J］.力學(xué)進(jìn)展，2007，37(3):398-408.TANG Zhi-ping.Progress and application of combined compression and shear wave loading technique［J］.Advances in mechanics，2007，37(3):398-408.

［2］ 崔云霄，盧芳云，林玉亮，等.一種新的高應(yīng)變率復(fù)合壓剪實驗技術(shù)［J］.實驗力學(xué)，2006，21(5):584-590.CUI Yun-xiao，LU Fang-yun，LIN Yu-liang，et al.A new combined compression-shear loading technique at high strain rates［J］.Journal of Experimental Mechanics，2006，21(5):584-590.

［3］ Tong W.Pressure-shear stress wave analysis in plate impact experiments［J］.International Journal of Impact Engineering，1997，19(2):147-164.

［4］ Clifton R J，Klopp R W.Pressure-shear plate impact testing［Z］.American Society for Metals，1985:230-239.

［5］ Abou-Sayed A S，Clifton R J.Pressure-shear impact of 6061-T6 aluminum due to oblique-plate-impact［J］.Journal of Applied Mechanics，1977，44(1):85-88.

［6］ Gilat A，Clifton R J.Pressure-shear waves in 6061-T6 aluminum and alpha-tianium［J］.Journal of the Mechanics and Physics of Solids，1985，33(3):263-284.

［7］ Huang H，F(xiàn)eng R.A study of the dynamic tribological response of closed fracture surface pairs by Kolsky-bar compression-shear experiment［J］.International Journal of Solids and Structures，2004，41(11/12):2821-2835.

［8］ Meyer L W，Krüger L.Shear testing with hat-shaped specimen，metals handbook:Mechanical testing and evaluation［Z］.American Society for Metals，2000:994-1000.

［9］ Rittel D，Ravichandran G，Venkert A.The mechanical response of pure iron at high strain rates under dominant shear［J］.Materials Science and Engineering:A，2006，432(1/2):191-201.

［10］ Rittel D，Lee S，Ravichandran G.A shear-compression specimen for large strain testing［J］.Experimental Mechanics，2002，42(1):58-64.

［11］ Griffiths L J，Martin D J.Study of dynamic behavior of a carbon-fiber composite using split Hopkinson pressure bar［J］.Journal of Physics，D:Applied Physics，1974，7(17):2329-2344.

［12］ 唐志平，胡時勝，王禮立.動態(tài)力學(xué)量的光電測試法［J］.實驗力學(xué)，1991，6(4):401-406.TANG Zhi-ping，HU Shi-sheng，WANG Li-li.Photoelectric method for measuring dynamic mechanical variables［J］.Journal of experimental mechanics，1991，6(4):401-406.

［13］ Ramesh K T，Narasimhan S.Finite deformations and the dynamic measurement of radial strains in compression Kolsky bar experiments［J］.International Journal of Solids and Structures，1996，33(25):3723-3738.

［14］ Li Y，Ramesh K T.An optical technique for measurement of material properties in the tension Kolsky bar［J］.International Journal of Impact Engineering，2007，34(4):784-798.

［15］ 沙定國.誤差分析與測量不確定度評定［M］.北京:中國計量出版社，2003.

［16］ Meyers M A.Dynamic behavior of materials［M］.New York:John Wiley＆Sons，Inc，1994.