鄭 星，黃?，?，毛勇建，張 軍，周 東

（中國工程物理研究院 總體工程研究所，四川 綿陽 621900）

1 引 言

爆炸波模擬裝置是一種對武器進行沖擊波壓力加載的試驗裝置，通常應用于抗爆試驗研究[1-2］。根據(jù)試驗需求，爆炸波模擬裝置可以產(chǎn)生平面波和曲面波，試驗時將被測目標置于爆炸波模擬裝置實驗段或管口，管內炸藥爆炸形成沖擊波作用在被測目標上，完成對被測目標的考核。在對被測目標進行抗爆考核時，沖擊波的幅值大小、正壓持續(xù)時間、形態(tài)、速度等參數(shù)是直接評估被測目標毀傷程度的重要指標。

目前，沖擊波參數(shù)的測試方法有等效靶法、電測法、薄膜測試法和光測法等[3，4］。爆炸波模擬加載試驗中，在實驗段大多采用平面波的加載方式，為了監(jiān)測爆炸波是否為平面波，通常在爆炸波模擬裝置不同段的同一圓周上布置多個壓力傳感器，通過壓力峰值來間接判斷沖擊波的狀態(tài)，但是該方法只能監(jiān)測管壁的壓力，無法對整個截面的壓力狀態(tài)進行判斷[5-7］。在對爆炸波速度監(jiān)測方面，通常是在爆炸波傳播的同一方向上布置多個壓力傳感器，通過壓力傳感器峰值時間來計算爆炸波的運動速度[8］。另外，為了確定加載在被測目標上的波是否為平面波，可在被測目標的迎爆面安裝多個傳感器，通過傳感器的峰值來判斷沖擊波狀態(tài)，但是由于壓力傳感器測試峰值受傳感器安裝高度，反射壓力等因素的影響從而引入較大的誤差，給判斷是否為平面波帶來困難。因此，雖然爆炸沖擊波的幅值、持續(xù)時間等參量可以利用沖擊波壓力傳感器進行獲取，其速度亦可通過布置多個測點進行計算，但是沖擊波的形態(tài)始終無法通過電測的方法進行直觀的反應，只能通過光學的方法進行觀測。高速紋影測試技術能夠實現(xiàn)空氣中流場圖像測量的目的，能夠直觀反映波紋的傳播和圖像顯示[9-12］。李斌等[13］利用光學紋影的方法對彈丸爆炸中的沖擊波威力進行了測量，實現(xiàn)了自由場中沖擊波流場的可視化。張雄星等[14］利用光學紋影方法實現(xiàn)了溫度場的定量測量，通過紋影技術獲取了不同溫度下加熱平臺上方的紋影圖像。胡洋等[15］利用“Z”字型光路的設計方案獲取了礦井瓦斯/空氣預混氣體爆燃過程的光學流程圖像。但是，針對爆炸波模擬裝置出口處的沖擊波紋影測量還鮮有文獻報道。由于爆炸波模擬裝置出口沖擊波形成時通常伴隨有火光和煙霧等，因此除了需要高靈敏度的高速相機外，還需要對出口處的其他非沖擊波產(chǎn)物進行處理。

本文通過高速CCD相機和一系列光學元件，利用反射式紋影測試原理，設計了高速、高靈敏度、遠距離激光紋影系統(tǒng)。通過該系統(tǒng)對簡易爆炸波發(fā)生裝置出口處的沖擊波演化過程及形態(tài)進行了測試和分析，為理解爆炸沖擊波的形成演化以及武器裝備在爆炸波沖擊作用下的考核評估提供圖像數(shù)據(jù)。

2 紋影測量原理

紋影測量法是利用光通過空氣擾動區(qū)域時光線發(fā)生折射，導致到達相機的光的照度發(fā)生變化，從而顯示擾動區(qū)狀態(tài)的一種方法。將紋影法與時間分辨率較高的高速相機相結合，便構成了高速紋影測量系統(tǒng)。

2.1 折射率梯度

光通過擾動區(qū)后，光線會發(fā)生一定程度的偏折。如圖1所示。

圖1 光波陣面經(jīng)擾動區(qū)的偏折Fig.1 Deflection of light wave front through disturbed region

平面光波陣面沿z軸傳播，光波陣面同時到達擾動初始界面O上的X和X′點，設擾動區(qū)折射率沿x方向變化，則經(jīng)過Δt后，X和X′點光線傳播的距離分別為Δz=cΔt n0，Δz′=cΔt n，其中n0和n分別為X和X′處的折射率。設：

2.2 基本原理

試驗中以激光光源作為點光源，通過準直擴束后將平行光束穿過含有可壓縮流場的試驗段，受擾動的光束發(fā)生偏折，通過紋影接收透鏡（或球面反射鏡）刀口系統(tǒng)成像，發(fā)生偏折的光束在成像面上發(fā)生照度的變化，從而得到?jīng)_擊波的紋影圖像。

如圖2所示，光源Ls通過透鏡L1和L2后成像在刀口面上，刀口位于第二球面鏡L2的焦平面上，圖中刀口垂直于圖面。其中，光源為點狀或與刀口平等的窄縫。照相透鏡Le將擾動區(qū)（試驗段）的像成像在成像裝置（CCD）上，從而記錄試驗區(qū)所需的紋影圖像。

由光學知識可知，圖2中刀口如果切去部分光源像，則在CCD上的光強就會減小。設剩余光源像的寬度為w，垂直圖面的高度為h，當試驗段沒有擾動時，到達CCD上任意點（x，y）的光強I是常數(shù)，其均勻照度可以精確表示為：

其中：I0為光源的初始強度，k為光源到刀口路徑上的光強損失的吸收系數(shù)，f為照相透鏡Le的焦距。（7）式忽略光學系統(tǒng)的幾何相差。

當試驗段空氣發(fā)生擾動，經(jīng)過試驗段的光線會由于試驗段折射率的變化發(fā)生偏折，設角度為ε，則這部分光線所成相應的光源像在垂直刀口方向上移動距離Δw。設偏折角沿x方向上的垂直分量為εx，則Δw=f2tanεx≈f2εx。

當平行光通過試驗段時，CCD上任一點的光強變化為ΔI=kI0(Δwh/f)，設吸收系數(shù)k不變，光強的相對變化為ΔI I=Δw w=εx(f2w)，將（5）式帶入，則可以通過CCD像中局部光強變化計算出試驗場中折射率梯度：

其中：w、f參數(shù)已知，光強的相對變化ΔI I可以通過CCD上的灰度變化讀取。因此，系統(tǒng)能夠檢測出的偏折角為：

由式（9）可知，系統(tǒng)能夠檢測出的偏折角與像寬成正比，與透鏡的焦距成反比。通常在一個實驗中透鏡的焦距固定，設系統(tǒng)能夠檢測到的最小光強相對變化為ΔEmin=(ΔI I)min，則系統(tǒng)能夠檢測出的最小偏折角為εmin=ΔEmin(w f)。因此，系統(tǒng)像寬越寬，能夠檢測的最小偏折角越大。設f w為紋影系統(tǒng)的靈敏度，靈敏度越大，系統(tǒng)能夠探測到的偏折角越小，系統(tǒng)的靈敏度越高。

圖2 紋影儀基本原理示意圖Fig.2 Principle diagram of schlieren instrument

3 測試系統(tǒng)設計方案

3.1 爆炸沖擊波加載系統(tǒng)

為了便于紋影系統(tǒng)觀察，利用鋼桶設計了一個能夠產(chǎn)生多個沖擊波陣面的簡易爆炸沖擊波發(fā)生裝置。如圖3所示，炸藥類型為TNT，藥量46 g。

圖3 簡易爆炸沖擊波發(fā)生裝置Fig.3 Simple explosion shock wave generator

圖3中，藥柱懸掛于鋼桶中部位置，藥柱距離鋼桶上部、下部和底部的距離均不相同。炸藥爆炸后，由于鋼桶內壁的反射作用，會產(chǎn)生多個沖擊波陣面，位于出口處的高速紋影系統(tǒng)即可觀察到多個爆炸沖擊波形態(tài)。

3.2 高速紋影測試系統(tǒng)

由于藥柱爆炸后形成的沖擊波有可能夾雜許多碎片損壞光學設備，因此本文采用的紋影測試原理為反射式，可以有效避免沖擊波或小飛濺物對光學器件的損傷，測量系統(tǒng)布局如圖4所示，兩個主反射鏡之間距離約50 m，主反射鏡焦距設置為5 m。系統(tǒng)采用波長520 nm綠色激光作為光源，光源通過3 mm狹縫和反射鏡后在兩個主反射鏡區(qū)間形成平行光，主反射鏡2通過反射鏡聚焦于刀口，最后通過CCD相機將檢測區(qū)域像成像于CCD陣面上，高速紋影測量系統(tǒng)布局如圖4中所示。由于爆炸沖擊波速度較快，因此CCD相機采用高速相機進行拍攝，拍攝幀頻4萬幅每秒，幅間間隔25μs。

圖4 高速紋影測量系統(tǒng)布局Fig.4 Layout of high speed schlieren measurement system

3.3 沖擊波壓力傳感器測量系統(tǒng)

為了對比紋影測量系統(tǒng)測量結果，在紋影成像范圍內布置自由場沖擊波壓力傳感器如圖4所示。其中，傳感器上升時間小于等于2μs，沖擊波壓力測量量程設置為2 MPa，數(shù)據(jù)采集系統(tǒng)采樣率為1 MSa/s，AD分辨率24 bit。

4 實驗結果與分析

4.1 沖擊波壓力傳感器測試結果

通過布置在簡易爆炸波發(fā)生裝置管口的壓力傳感器，測得爆炸沖擊波壓力曲線如圖5所示。由圖可見，由于藥柱爆炸以及壁面反射作用，產(chǎn)生了四個較為明顯的沖擊波峰值，分別用P1、P2、P3和P4表示，相鄰波峰值時間間隔如圖5所示，圖中以第一個沖擊波陣面（P1）峰值時刻為零時刻點進行計算。

圖5 沖擊波壓力曲線圖Fig.5 Curve of shock wave pressure

4.2 高速紋影系統(tǒng)測試結果

利用高速紋影系統(tǒng)對管口沖擊波形態(tài)進行測試，測得四個沖擊波陣面到達沖擊波壓力傳感器頂部時的狀態(tài)如圖6所示。圖中沖擊波陣面圖像為直接測量結果，未經(jīng)過任何處理。設第一個沖擊波陣面（P1）到達壓力傳感器頂端為零時 刻 點，經(jīng) 過0.55 ms、2.875 ms和4.375 ms后，后續(xù)三個沖擊波陣面依次到達壓力傳感器頂端。由圖可見四個沖擊波陣面形態(tài)在視場范圍內清晰可見，四個沖擊波陣面形態(tài)近似為平面波。

沖擊波壓力傳感器和高速紋影系統(tǒng)測得沖擊波陣面到達傳感器時刻點如表1所示。

圖6 四個沖擊波陣面狀態(tài)Fig.6 States of four shock front

表1 四個沖擊波陣面到達傳感器時間Tab.1 Time of four shock wave fronts arrival sensor （ms）

由表1可知，高速紋影系統(tǒng)測得四個沖擊波陣面到達壓力傳感器時間與沖擊波壓力傳感器測得結果基本吻合。測量誤差主要由圖像測量中視覺誤差和讀取誤差引起，沖擊波壓力傳感器測量系統(tǒng)采樣率為1 MSa/s，時間精度為1μs，而高速相機幅頻為40 000 FPS，相鄰兩張圖片之間時間精度為25μs，因此在圖片選取和像素讀取時引入誤差。由表1可見，兩者所測結果在時間量級上基本一致。因此，沖擊波壓力傳感器所測曲線中的四個壓力峰值所對應的圖像即為高速紋影系統(tǒng)所觀測到的四個沖擊波陣面圖像，壓力曲線四個峰值對應沖擊波陣面紋影圖像如圖7所示。

圖7 沖擊波壓力曲線四個峰值對應紋影圖像Fig.7 Schlieren image of four peaks in shock wave pressure curve

由式（8）以及圖6可知，受沖擊波波陣面影響，檢測區(qū)域沿x方向的折射率發(fā)生了變化。受折射率梯度影響，光線沿x方向發(fā)生了不同程度的偏折，因此從高速相機上觀測到不同程度的明暗圖像，即沖擊波陣面。

4.3 沖擊波速度測量分析

利用沖擊波的紋影圖像可以準確地獲取沖擊波在觀察區(qū)域內的沖擊波速度，實驗時，高速相機拍攝幀頻為40 000 FPS，則相鄰兩張圖片間隔0.025 ms，通過標定高速相機相鄰兩張圖片之間的距離，得到兩張圖片之間沖擊波移動距離約為9.1 mm。整個視場范圍內，第一個沖擊波陣面在相鄰兩張紋影圖像之間移動的距離基本相等，0.975 ms后第一個沖擊波陣面運動出紋影視場范圍內，移動距離約355 mm。第一個沖擊波陣面運動距離對應時間曲線如圖8所示。其中以第一個沖擊波陣面剛進入紋影視場為零時刻點，第一個沖擊波陣面運動軌跡如圖9所示，經(jīng)過0.5 ms和0.975 ms后第一個沖擊波陣面狀態(tài)分別如圖9中左數(shù)第三與第五幅圖片所示。由圖8可見，在紋影圖像范圍內，沖擊波陣面移動的距離與時間呈線性關系，在視場范圍內沖擊波各時刻點速度相同，經(jīng)計算，視場范圍內第一個沖擊波陣面速度約為364 m/s。

圖8 第一個波陣面移動距離-時間曲線Fig.8 Moving distance and time curve of the first wave front

圖9 第一個沖擊波陣面運動軌跡Fig.9 Trajectory of the first shock front

由沖擊波壓力曲線可知，在第二個沖擊波陣面（P2）和第三個沖擊波陣面（P3）之間形成了諸多較小的反射壓力波，由圖10中紋影圖像可見，該部分反射波清晰可見，主要是由沖擊波壓力傳感器支架與爆炸沖擊波發(fā)生裝置來流方向的反射引起。

圖10 反射沖擊波狀態(tài)Fig.10 State of reflected shock wave

5 結 論

針對簡易爆炸波發(fā)生裝置出口處的沖擊波圖像測量，采用反射式紋影測量技術，完整獲取了出口處多個沖擊波陣面的圖像狀態(tài)。通過與出口處沖擊波壓力傳感器壓力峰值進行對比，兩者所測得的波陣面數(shù)量吻合，并且通過對單個波陣面的運動歷程進行計算，獲取了測量區(qū)域內沖擊波的傳播速度。因此，本文設計的高速紋影測試系統(tǒng)能夠有效獲取爆炸沖擊波的運動圖像狀態(tài)，并且可以通過其運動歷程進行相應的速度計算。測量結果對理解沖擊波的傳播和演化規(guī)律以及對爆炸波模擬裝置的改進和武器裝備的評估具有重要意義。