爆炸膨脹環實驗數據處理方法討論＊

湯鐵鋼，桂毓林，李慶忠，陳永濤，童慧峰，劉倉理

（中國工程物理研究院流體物理研究所沖擊波物理與爆轟波物理國防科技重點實驗室，四川 綿陽621900）

1 引 言

爆炸膨脹環實驗是在1維應力狀態假定下研究材料高應變率拉伸加載時的動態性能，包括本構關系和斷裂特性，長期以來，爆炸膨脹環實驗數據的處理一直困擾著研究者們。P.C.Johnson 等［1］、C.R.Hoggatt等［2］運用爆炸膨脹環實驗技術測試了許多工程材料的本構關系數據，當時只能測量膨脹環的位移歷史，數據處理時需要對膨脹環的位移歷史進行2次微分，數據精度一直困擾著實驗數據的可用性。直到R.H.Warnes等［3］利用激光速度干涉儀（VISAR）直接測量了爆炸膨脹環的徑向膨脹速度，克服了求流動應力時位移關于時間2次微分的困難。在數據處理時，R.H.Warnes等對速度歷史曲線采用2次關系式進行了擬合，然后進行微分計算加速度。R.H.Warnes等同時指出，采用更高階的關系式進行擬合，可以對速度進行更好地近似擬合，但是求解出的應力應變關系卻令人費解；他們還指出，其他的一些光滑數據的方法也沒有得到實驗檢驗。F.Llorca等［4］通過數值模擬對膨脹環實驗中速度歷史的處理進行了分析，嘗試直接對測試獲得的離散速度數據進行處理，結果獲得的流動應力應變關系曲線振蕩非常劇烈，而采用線性和2項關系式擬合后，得到的流動應力應變關系與SCG 模型結果比較接近。

本文中，嘗試采用不同的數據處理方法對爆炸膨脹環實驗獲得的速度曲線進行處理，討論數據窗口選取對處理結果的影響，同時對數據處理方法中存在的一些問題進行分析。

2 膨脹環運動分析

對于自由飛行階段的膨脹環，膨脹環中的徑向應力很小，近似看作零，因此膨脹環僅在環向應力作用下作減速運動，如圖1所示，取其中一環向單元建立運動方程

金屬環在自由膨脹期間徑向應力為零，得到周向應力的運動方程

式中：r0為膨脹環初始半徑。將式（4）對時間t求導數，得

圖1 膨脹環運動分析示意圖Fig.1 Movement analysis of an expanding ring

3 膨脹環速度曲線

3.1 速度歷史特征及數據窗口的選取

在爆炸膨脹環實驗中，采用激光干涉技術可以獲得膨脹環的徑向速度歷史，以1 發無氧銅膨脹環實驗結果為例進行說明，速度歷史如圖2所示。從速度歷史可以看出，膨脹環具有較長一段自由膨脹過程，即速度歷史上的B 段，通常稱之為有用數據窗口，A 段為沖擊加載及應力波在膨脹環壁內的反射階段，C 段為膨脹環發生了失穩或者斷裂階段，對于處理膨脹環材料的應力－應變－應變率關系，A 段和C 段均為無用數據。為了揭示數據處理方法對結果的影響，本文中選取的這條速度曲線的B 段的斜率不單調變化。

圖2 膨脹環速度歷史及數據窗口選取Fig.2 Radial velocoty－time history of the expanding ring and the selection of data windows

如何合理地選取數據窗口B 段，是數據處理的前提。對于此發實驗，通過回收樣品測量確定膨脹環的斷裂應變約0.14，對應的斷裂時間約50μs，顯然，數據窗口的選取只能在50μs之前。從速度歷史上看（見圖2），在42～50μs之間速度有些振蕩，數據是否可用需要討論。為了更好地檢驗數據處理方法，擬在速度曲線上選取大部分區域重合的3個數據窗口，分別記為窗口1、窗口2、窗口3，如圖2所示。窗口1對應時間12～42μs，認為是最理想的數據窗口，窗口2對應時間10～42μs，窗口3對應時間10～50μs。從速度曲線看，窗口2相對窗口1沒有明顯的趨勢變化，而窗口3中的后8μs速度趨勢有明顯變化。

3.2 數據擬合

在速度曲線上對選取的窗口數據截取后進行數據擬合，分別進行線性關系式、2次關系式、3次關系式擬合，擬合后進行微分求解加速度。線性關系式、2次關系式、3次關系式的表達式分別為

對于窗口1～3中的數據采用以上3種關系式進行擬合，求解應力時，首先要對擬合的速度公式進行求導得到加速度，所有關系式中的常數項a 求導后消失。線性擬合時，加速度（即速度斜率）為常數b，2次關系式擬合時，加速度為時間的線性函數，3次關系式擬合時，加速度為時間的2次函數。

4 數據處理結果分析

4.1 擬合方法對數據處理結果的影響

采用3種擬合關系式對實驗數據進行處理，首先對圖2 中窗口1（12～42 μs）數據進行處理，獲得的應力應變關系見圖3。線性擬合關系式獲得的應力應變關系也基本上呈線性關系，2 次關系式擬合時獲得的應力應變關系與線性擬合在所選的數據窗口中較接近，而采用3次關系式擬合時，應力應變關系呈現為上凹的曲線。

4.2 數據窗口選取對結果的影響

對于測試獲得的含有一定噪聲且離散的速度數據，很難準確界定有效數據窗口的邊界，因此，數據窗口的選取存在一定的人為性。通過改變數據窗口的寬度，來討論數據窗口對應力應變關系的影響。

首先采用線性關系式和2次關系式擬合對3個窗口數據進行處理，獲得的應力應變關系如圖4所示。采用線性擬合時，窗口寬度的改變直接影響速度的斜率（常數），即直接影響應力的幅值，因此獲得應力應變關系為3條相互平行的直線（實際上是與位移線性相關的曲線），其中窗口1和窗口2的結果比較接近，而窗口3的結果偏離稍大。采用2次關系式擬合時，窗口1、窗口2數據處理獲得的應力應變關系表現為略有下凹的曲線，而窗口3數據處理獲得的應力應變關系卻為明顯的上凹曲線，尤其是后半段與前2個數據窗口結果偏離較大，這與原始速度數據的后期振蕩是相關的。

采用3次關系式擬合對3個窗口數據進行處理，獲得的應力應變關系如圖5所示。對于3個不同的數據窗口，采用3次關系式擬合處理結果均為明顯的上凹曲線，數據窗口的選取對于3次關系式擬合處理結果的影響不及低階關系式明顯，但是處理出的應力應變關系卻不同于一般的拉伸應力應變關系曲線特征。

圖3 窗口1數據3種擬合方法處理結果Fig.3 Processing results of data window 1by three fitting expressions

圖4 3個數據窗口線性擬合及2次關系式擬合處理結果Fig.4 Processing results of three different data windows

4.3 分析討論

從應力計算表達式看，采用線性擬合時，加速度為常數，應力與位移（或應變）成線性增長關系，反映了材料的應變硬化效應。實際上，由公式（5）可知，隨著位移的增加，膨脹環的拉伸應變率下降。如果考慮材料的應變率效應，材料中的應力也應該隨著應變率的降低而降低，線性擬合不能反映材料的應變率效應。因此，對于應變率敏感的材料，線性擬合是不合適的。

以上研究發現，對于膨脹環速度曲線上人為選取的任何窗口數據，3次關系式擬合都是上凹曲線，說明高階關系式雖然能較好擬合膨脹環的速度歷史。但由此處理得到的應力應變關系卻與物理規律似乎背離，本文中推測R.H.Warnes等［3］在數據處理時可能也出現類似問題，即他所說的令人費解的原因。這到底是原始速度數據的問題，還是數據處理方法的問題，需要更進一步的分析。

圖5 3個數據窗口3次關系式擬合處理結果Fig.5 Processing results of three different data windows by trinomial fitting

圖6 5種不同擬合關系式對窗口1數據的處理結果Fig.6 Processing results of data window 1by five different fitting expressions

在不能明顯判斷數據是否可用的條件下，隨著窗口數據的增加（即增加窗口的時間長度），2次關系式擬合處理結果也出現了明顯的偏離，如圖4中的窗口3數據結果。到底是數據處理方法的問題還是數據窗口選取的問題呢？本文中嘗試了另幾種多次關系式，如

以及他們的組合關系式，對速度窗口中的數據進行擬合處理。結果發現，在窗口1中，多種關系式擬合結果基本一致，見圖6，當最高次不大于2時，得到的結果偏差相對較小，而最高次項大于2時，獲得應力應變關系為上凹曲線，從形式上可能偏離了物理規律。而對窗口2和窗口3中數據進行非線性擬合處理得到的應力應變關系均出現不同程度的發散，由此看來，數據窗口的準確選取是采用擬合方法進行膨脹環實驗數據處理的基礎。

5 數據光滑與差分方法

在膨脹環實驗數據處理時為什么總是先對速度歷史進行擬合呢？主要原因是測試的速度歷史是離散的數據，而且測試的速度數據包含了噪聲，速度曲線本身不光滑，直接差分求解加速度會帶來極大的偏差。實際上，目前對于數據的光滑技術已經很成熟了，可以先對速度數據進行光滑，然后直接進行差分求解加速度。理論上講，直接進行求解，能更好地反映每個時刻速度與應力的對應關系。另一個好處是，直接差分計算對窗口選取的依賴性較弱，因為差分方法可以計算出瞬時加速度，只依賴于瞬時的速度數據，不再需要擬合關系式。但由于膨脹環實驗是在1維應力假定條件下進行的，數據處理時仍熱需要確定合理的數據窗口，以確保所選取的膨脹環拉伸過程滿足1維應力狀態。

為了更好地澄清3次擬合關系式處理結果形式上異常的原因，采用Mathcad中的樣條函數對速度曲線進行高階光滑，然后進行直接差分計算，處理得到的應力應變關系與擬合關系式處理結果同列于圖7中。令人驚奇的是，在窗口1中，3次關系式擬合處理結果與直接差分計算結果吻合很好。

由此看來，采用高階關系式擬合進行處理時，出現結果異常的原因并非數據處理方法的問題，主要原因是來自原始數據的測試精度。采用VISAR 進行速度測試時，是典型的點測試，對于環向尺寸較大的膨脹環，VISAR 探頭對準的任意一點，有可能是頸縮、斷裂等現象發生的點，獲得的速度歷史就不能反映膨脹環的均勻膨脹過程。因此，在膨脹環實驗中，速度測試最好是多點，采用平均速度進行數據處理會提高數據精度，單點速度測試會給數據處理結果帶來很大的不確定性。

圖7 直接差分方法與擬合方法處理結果的比較Fig.7Comparison of results by the fitting methods and the direct difference method

6 結 論

（1）對于爆炸膨脹環實驗，速度數據窗口的合理選取是進行數據處理的前提，數據窗口的選取既要避開前期由于加載引起的速度振蕩，也要避開后期由于膨脹環發生頸縮、斷裂引起的速度振蕩；

（2）速度曲線在數據窗口內斜率單調變化時，可以采用低階關系式擬合進行處理，擬合關系式不超過2階時，處理結果會比較接近；

（3）采用數據光滑技術和差分方法相結合的數據處理方法，能較好地反映速度歷史和瞬時應力的對應關系，但同時要求速度測試具有較高的精度；同時澄清了在數據處理時采用高階關系式擬合出現異常現象的原因是速度點測試方法帶來的。

［1］ Johnson P C，Stein B A，Davis R S.Measurement of dynamic plastic flow properties under uniform stress［C］∥Symposium on the Dynamic Behavior of Materials.1963：195－198.

［2］ Hoggatt C R，Recht R F.Stress－strain data obtained at high rates using an expanding ring［J］.Experimental Mechanics，1969，9（10）：441－448.

［3］ Warnes R H，Duffey T A，Karpp R R，et.al.An improved technique for determining dynamic material properties using the expanding ring［C］∥Meyers M A，Murr L E.Shock Waves and High－strain－rate Phenomena in Metals.New York：Plenum Press，1981：23－36.

［4］ Llorca F，Juanicotena A.Expanding ring test：Numerical simulation－application to the analysis of experimental data［J］.Journal De Physique IV，1997，7：325－340.