溫慶榮,李 想,魏夢琦,劉宏旭,胡建飛
(華北光電技術研究所,北京 100015)
現有軍用光電設備中,紅外成像、紅外制導、紅外對抗、戰艦導航、車載夜視、紅外探測和紅外瞄準等以紅外應用為主的技術所占的比重越來越重要,而紅外光學元件是紅外光電系統中的關鍵部件。紅外材料不僅需要具有良好的紅外透過性能,要同時承受在高速飛行中雨水、沙塵和冰雹等顆粒物的沖擊損傷和長期海面鹽霧的化學腐蝕。
由于ZnS玻璃可透過可見光、激光到紅外多個波段,因此在紅外光電設備中得到廣泛的應用。而ZnS屬于特殊的脆性材料,力學性能參數并不唯一,而且不同批次的材料其力學性能差異也較大;并且在高速動態下,ZnS材料的力學性能會隨著碰撞速度的增加發生變化。因此,本文通過設計試驗方案對本批次的ZnS樣件開展力學性能的測試,并針對試驗結果進行數據分析,得出實際的力學參數,為后續的仿真分析提供真實的數據。
為了能夠與試驗的機器匹配,本文針對拉伸和壓縮試驗分別選用直徑10 mm和直徑5 mm的ZnS圓柱玻璃,高度均為4.2 mm,樣件如圖1所示。

圖1 玻璃樣件
根據試驗目標,設計的試驗矩陣如表1所示。針對玻璃材料分別開展準靜態、中高應變率下測試其壓縮工況下極限應變、極限應力的試驗。

表1 試驗矩陣
2.2.1 準靜態試驗方法
ZnS玻璃屬于脆性材料,極限應變非常小,并且試驗樣本尺寸又比較小,在進行準靜態壓縮試驗時,由于試驗的壓縮機的壓頭與被壓縮的樣本之間有一定的間隙,從試驗機的位移數據來獲得樣件的應變數據不準確。因此,為了獲得玻璃樣件較為準確的應變數據,采用數字圖像處理技術(DIC)的方法,首先要對玻璃樣件噴射散斑。圖2為玻璃樣件噴射散斑以后的照片。

圖2 玻璃樣件噴射散斑
準靜態試驗是在萬能試驗機上進行的,電子萬能材料試驗機簡稱電子萬能試驗機,是材料力學性能測試的專用設備,主要用于材料的拉伸、壓縮、彎曲、剪切等力學性能試驗,并采集力-位移曲線。由于試驗機在壓縮過程中存在間隙,因此壓縮位移的采集是通過在試驗機壓頭噴射散斑然后使用非接觸測量的方式來獲取,進而可獲取更為準確的力-位移曲線。圖3是準靜態試驗的試驗機以及樣本數據采集裝置。

圖3 試驗機及數據采集裝置
2.2.2 準靜態試驗數據處理方法
根據獲取的力-位移曲線,按照下式來計算工程應力σE和工程應變εE,其計算表達式如下所示。
σE=P/As
(1)
εE=ΔL/L
(2)
εT=ln(1+εE)
(3)
σT=σE(1+εE)
(4)
將(1)和(2)式代入(3)和(4)式中,可以計算得到真實的應力和應變。式中,As和L分別表示試件的初始截面積和長度,ΔL為試件的壓縮量。
2.2.3 中高速試驗方法
中高速試驗是在霍普金森壓桿試驗機上進行的,分離式霍普金森壓桿(split Hopkinson pressure bar,簡稱SHPB)[1-2],是現今最為廣泛使用并且被認為有效的測試材料剛應變率下力學特性的試驗裝置。其主要特點是可以實現試件的高應變率變形,同時保持試件內的動態應力平衡或者說試件內軸向應力梯度接近于零。將玻璃樣件貼上應變片,并通過潤滑劑夾緊在輸入桿和輸出桿中間,同時將應變片與數據采集系統連接。本試驗的電磁桿是單向加載,根據試驗目標應變率范圍來估算氣壓。圖4為霍普金森壓桿試驗機試驗過程圖片。

(a) (b)
2.2.4 中高速試驗數據處理方法
霍普金森壓桿工作原理示意圖如圖5所示,在分離式Hopkinson 壓桿系統中,加載力是入射應力波,通過子彈撞擊入射桿產生,并在試驗過程中測定輸入桿和輸出桿的應變。

圖5 霍普金森壓桿加載示意圖
在試驗中,粘貼在入射桿上的應變片用于采集入射波和反射波,粘貼在透射桿上的應變片用于采集透射波。根據一維應力波理論,樣件內部的應力、應變和應變率可以利用所采集到的入射波、反射波和透射波計算得到[3-6]。
假定E是加載桿的彈性模量,A和As分別為加載桿和試驗件的截面積,L為試驗件的厚度,C0為波速。εI(t),εR(t)和εT(t)分別代表了應變片測得的入射應變、反射應變和透射應變,按照下述推導可以得到應力、應變及應變率的表達式[7]。
設P1,P2代表樣件兩端面的壓力,則:P1=EA[εI(t)+εR(t)],P2=EAεR(t),那么式樣的平均應力可表示為:
(5)
樣件兩端面的質點速度V1,V2可以分別表示成V1=C0[εI(t)-εR(t)],V2=C0εR(t),樣件的應變率可表示為:
(6)
對(5)和(6)進一步簡化可得式(7)和(8),應變率對時間積分即可得工程應變(9)。
(7)
(8)
(9)
若由樣件均勻性假定則有εI(t)+εR(t)=εT(t),則上式可以進一步簡化:
(10)
(11)
(12)
本文中的試驗裝置的參數如下:
壓桿的橫截面積:
A=πR2=π·12.552=494.8(mm2)
試件的橫截面積:
As=πr2=π·2.52=19.63(mm2)
壓桿的彈性模量:E=123000(MPa)
試件長度:L=4.2(mm)

由于準靜態試驗尚未采集到較為準確的彈性模量,因此采用中高速試驗測量所得的平均彈性模量(64.08 GPa)來對準靜態試驗的應變數據進行計算。計算后的準靜態試驗的應力-應變曲線如圖6所示。

圖6 準靜態壓縮試驗應力-應變曲線
在中高應變率下玻璃的應變太小并且速度太快,應變片有可能提前脫落從而產生無效的數據。本文從試驗中取6次應變數據,通過配合壓桿信號處理得到的應力數據來計算樣件的彈性模量。
3.2.1 樣件的彈性模量
圖7是在壓縮試驗中,通過6次有效的試驗數據得出的樣件的應力-應變曲線,通過這些數據可以計算出彈性模量。

圖7 測量彈性模量的應力-應變曲線
根據彈性模量的計算公式E=σ/ε,得出的圖7中6次試驗的彈性模量在表2中所示。彈性模量的范圍在48.7 GPa到88.51 GPa之間。

表2 壓縮試驗測得的彈性模量
3.2.2 樣件的極限應力

表3 壓縮試驗應變率和極限應力
3.2.3 樣件的極限應變
在計算極限應變時,本文假設樣件是理想彈性體,取3.2.1中計算的彈性模量值作為輸入。其中003、006和013編號的試件未測到彈性模量數值,則根據其應變率的值來選擇相近應變率下試件的彈性模量值作為其值。最終不同試驗編號樣件的極限應變如表4所示。

表4 壓縮試驗的極限應變
3.2.4 應力-應變曲線圖
中高應變率下修正過后的應力-應變曲線如圖8~圖10所示。圖11是應力與應變率的關系圖,從圖中可以看出,在中高應變率的范圍內樣件的極限應力沒有單純隨著應變率的增高而增高,而是呈現浮動變化的,說明樣件在100~500 s-1應變率范圍內的應變率效應不明顯。

圖8 應變率100~200 s-1范圍的樣件應力-應變曲線

圖9 應變率200~300 s-1范圍的樣件應力-應變曲線

圖10 應變率400~500 s-1范圍的樣件應力-應變曲線

圖11 樣件應力與應變率的關系圖
根據處理后的試驗數據,計算的力學參數主要包括:極限應力,壓縮應力-應變曲線中第一個最大應力值;彈性模量,壓縮應力與應變片測得的應變曲線斜率獲得;極限應變,極限應力比上彈性模量獲得。表5為樣件的力學參數計算結果。

表5 樣件的力學參數計算結果
本文針對ZnS樣件開展了壓縮試驗,并獲取了有效的數據,通過數據分析得出如下結論:
1)該批次的ZnS材料準靜態下的極限應力為248.59 MPa、極限應變為0.00385;
2)該批次的ZnS材料中高速狀態下的極限應力為518.89 MPa、彈性模量為65.24 GPa、極限應變為0.0083;
3)在中高應變率的范圍內樣件的極限應力沒有單純隨著應變率的增高而增高,而是呈現浮動變化的,說明樣件在100~500/s應變率范圍內的應變率效應不明顯。
因此,根據本文所得出的試驗據,后續針對該機載設備的力學仿真分析將能夠更加接近真實狀況,可以有效的縮短產品的研制周期,提升其機載環境適應性能。