郎文杰, 陳國光，田曉麗，辛長范

1. 中北大學信息與通信工程學院，山西 太原 030051 2. 中北大學機電工程學院，山西 太原 030051

炮體表面瞬態(tài)溫度分布的散斑干涉測試研究

郎文杰1, 陳國光2，田曉麗2，辛長范2

1. 中北大學信息與通信工程學院，山西 太原 030051 2. 中北大學機電工程學院，山西 太原 030051

為了可以實現(xiàn)對炮體表面瞬態(tài)溫度分布變化的實時監(jiān)測，同時克服傳統(tǒng)瞬態(tài)高溫測試儀器單點探測以及熱慣性大等局限性，設計了基于散斑干涉與光譜頻域分析相結(jié)合的瞬態(tài)溫度分布測試系統(tǒng)。系統(tǒng)采用散斑干涉技術(shù)將炮體瞬態(tài)溫度變化引發(fā)的微小形變轉(zhuǎn)換成散斑干涉條紋，再由傅氏變換完成干涉條紋形變到光譜分布的函數(shù)轉(zhuǎn)換，從而通過光譜分布函數(shù)反演任意采樣時刻上的溫度分布。實驗采用ZX-FB1型光纖測溫儀測試單點位置上瞬態(tài)溫度作為標準值，再由555 nm激光器與面陣CCD采集散斑干涉條紋，分別使用圖像識別法與傅氏變換法完成干涉條紋與瞬態(tài)溫度的算法匹配，從而反演瞬態(tài)溫度。實驗結(jié)果顯示，兩種方法均能實現(xiàn)瞬態(tài)溫度檢測，但基于傅氏變換頻譜分析的散斑干涉法精度更高，并且可以有效地克服由表面瑕疵、漆面磨損等問題造成的粗大誤差。

瞬態(tài)溫度分布； 散斑干涉技術(shù)； 頻譜分析； 傅氏變換

引 言

火箭彈是陸軍機械化武器中的重要組成部分，由于其具有殺傷力強、成本低、機動性高的特點而被廣泛使用，雖然精度低于制導武器，但隨著彈體結(jié)構(gòu)的不斷完善，使其打擊精度也在不斷地提升[1]。在彈體可靠性測試中主要包括彈體的材料力學特性、彈體飛行過程中空氣動力學特性以及抗干擾能力等。而對于彈體瞬態(tài)溫度場分布的研究對測試控制飛行姿勢、調(diào)整打擊精度具有重要意義[2]。

瞬態(tài)溫度場分布特性的測試被應用于眾多領(lǐng)域，例如彈藥爆炸過程分析、炮管炮體溫度應變特性、材料合成溫度變化特性分析等[3-4]。在傳統(tǒng)的瞬態(tài)溫度場測試過程中，一般采用測溫網(wǎng)絡的方法，將大量的溫度傳感器至于被測表面或被測區(qū)域中，然后實時記錄每個位置上的溫度變化，通過高速數(shù)據(jù)采集系統(tǒng)將各點位的溫度數(shù)據(jù)傳給處理器進行分析。其中有熱電阻探測器、熱電偶探測器、數(shù)字溫度計以及光纖光柵探頭等[5-9]。不同的測溫探測器應用略有不同，大部分主要適用于低溫測試，或者溫度緩變的條件，并且以上檢測方法均屬接觸式測溫方法，限制了其應用范圍。

散斑干涉檢測技術(shù)是通過對被測物表面微小形變量的測定從而反演其瞬時溫度變化的一種非接觸式溫度檢測方法[10]，其被廣泛應用于器件加工、材料合成測試、機械結(jié)構(gòu)穩(wěn)定性等方面[11-12]。該技術(shù)是利用圖像采集與處理的手段將散斑干涉圖中反映被測表面形變的信息提取出來，再反演被測物溫度場瞬態(tài)變化的，具有非接觸、高實時性、高精度等優(yōu)點。本系統(tǒng)采用散斑干涉技術(shù)結(jié)合傅氏變換獲取其光譜分布函數(shù)，對炮體瞬態(tài)溫度場分布進行分析，相比單純采用散斑干涉與圖像識別的溫度測試方法具有更好的系統(tǒng)魯棒性。

1 溫度場分布測試系統(tǒng)設計

基于散斑干涉的瞬態(tài)溫度測試系統(tǒng)包括激光器、準直系統(tǒng)、分束鏡、CCD以及處理電路等組成，如圖1所示。激光經(jīng)準直系統(tǒng)擴束準直到合理的光斑尺寸后進入分光棱鏡，分為兩束后一束照射被測物表面，另一束通過補償鏡照射在反射鏡上，兩束光反射后再次進入分束棱鏡后聚焦于CCD探測器上。在CCD上形成的散斑干涉條紋被采集并傳遞給處理系統(tǒng)，對數(shù)據(jù)進行頻域分析后最終計算得到被測物的瞬態(tài)溫度信息。

1.1 理論分析

對于特定的被測物體而言，其溫度的改變量會引起物體發(fā)生形變，雖然溫度引起的形變量通常很小，但是由于采用散斑干涉技術(shù)，非常微弱的形變也會在其干涉圖樣中有所體現(xiàn)，故依據(jù)熱力學定律可知，設壓強不變時，單位溫度所產(chǎn)生的物體變化可表示為

(1)

Fig.1 Transient temperature measurement system based on Speckle Pattern Interferometry

同理，當溫度不變時，單位壓強產(chǎn)生的物體變化可表示為

(2)

由此可得，被測物形變滿足以下方程

V(T，p)=[1+α(T-T0)-kp]V0(T0，0)

(3)

由于被測物炮體為鋼材，可近似看作彈性體，故有

(4)

將式(4)代入式(3)化簡可得

(5)

式中Δl為被測物表面形變量，ΔT為溫度改變量，l為初始長度，α為被測材料的溫度系數(shù)，kT為被測材料的壓強系數(shù)。

1.2 形變檢測的改進方法

傳統(tǒng)的散斑干涉測溫技術(shù)常常將采集得到的干涉圖像直接應用圖像處理軟件計算，通過圖像灰度變化量與溫度變化量的線性關(guān)系計算，該種方法運算簡單，速度快。但對于炮體瞬態(tài)溫度變化過程而言，其溫度變化范圍較大，該種方法引入的非線性明顯，對溫度測試精度影響大。

本系統(tǒng)對此進行改進，采用傅氏變換技術(shù)，將干涉條紋形變引起的頻譜變化與瞬態(tài)溫度建立函數(shù)關(guān)系，溫度與形變量是否線性不再影響系統(tǒng)對瞬態(tài)溫度的求解，而取決于時頻變換后中心波長的偏移量及振幅信息。故光譜分布函數(shù)為

(6)

式中，x表示形變量，I(x)表示x位置上的光強值，λ表示波長，ν表示波數(shù)。

2 瞬態(tài)溫度分布的計算

由于散斑干涉條紋與被測物表面的溫度相關(guān)，則在發(fā)生溫度變化前的反射波與參考波有以下函數(shù)關(guān)系

E(x，y)=a0eiφ0(x，y)+areiφr(x，y)

(7)

將式(7)及其共軛表達式相乘可得初始狀態(tài)下的光強表達式有

(8)

當炮體溫度發(fā)生瞬態(tài)改變時，其回波光的相位也發(fā)生相應的改變，其光強可表示為

(9)

通過將溫度瞬態(tài)變化過程中兩個不同時刻所產(chǎn)生的中心波長對應的幅值比求解出來就可以推導出其瞬態(tài)溫度了，其光強絕對值有

I(x，y)=|It-I0|=4a0ar|sin[(φ0-φr)+nπ]sin(nπ)|

(10)

其中，n=0, ±1, ±2, …。

而回波光相位可以表示為初始光相位為0時，兩光絕對值中的相位項，故有

(11)

綜上所述，將式(5)和式(10)代入式(11)，化簡可得

(12)

式中k表示壓強系數(shù)，α表示溫度系數(shù)，T0表示初始溫度，l表示材料長度，λ表示激光波長和E表示激光光波能量函數(shù)。由式(12)可知，當激光器選定時(實驗采用660 nm激光器)，其波長及功率為常數(shù)，故瞬態(tài)溫度主要由k和決定。

3 測溫實驗

3.1 設備及結(jié)果

系統(tǒng)采用法國Oxxius公司SLIM系列激光器(中心波長555nm)，探測器采用SI7000型面陣CCD，光譜分析算法及軟界面采用VS平臺獨立開發(fā)。實驗對某型炮管發(fā)射時炮體表面進行散斑干涉測試，干涉圖樣如圖2所示。

Fig.2 Speckle interference fringe at the time t0 and tx

圖2(a)中是初始時刻(t0)CCD采集得到的干涉條紋數(shù)據(jù)，圖2(b)中是tx時刻后CCD采集得到的干涉條紋數(shù)據(jù)。為了盡量減小由于彈體發(fā)射所造成震動等對溫度測試的影響，在實際測試過程中，t0的選取并不是從0時刻起算的，而是從彈體射出后炮體穩(wěn)定的時刻計算的，每次的測試是隨反饋閾值給定的，故不具體標注單次的起始時間值。從干涉圖樣可以看出，散斑干涉條紋發(fā)生扭曲，在圖像中雖然能看出有形變，但形變量并不明顯，而在本方法中采用傅里葉變換對該干涉條紋進行頻譜分析時，微小的形變量會被光譜分布函數(shù)的中心波長振幅偏移量及其比值精確地反演出來，從而得到精確的瞬態(tài)溫度分布數(shù)據(jù)。

3.2 數(shù)據(jù)對比

實驗的測試溫度區(qū)間為800～1 800 ℃，標定采用ZX-FB1型光纖測溫儀，溫度測試精度為1 ℃，測溫范圍滿足實驗要求。分別采用基于圖像識別的傳統(tǒng)散斑干涉法(image recognition-speckle pattern interferometry, IR-SPI)和基于傅氏變換頻譜分析法(Fourier transform-speckle pattern interferometry, FT-SPI)對某型炮彈發(fā)射過程中炮體(發(fā)射筒前框)溫度進行測試，實驗數(shù)據(jù)如圖3所示。

Fig.3 Comparison of the test data for gun body transient temperature

由圖3中溫度測試數(shù)據(jù)經(jīng)計算分析可知，采用IR-SPI和FT-SPI兩種方法均能有效地將被測物表面的瞬態(tài)溫度數(shù)據(jù)檢測出。通過數(shù)據(jù)圖可以看出，隨著時間的推遲，檢測效果與標準值越接近。其中，IR-SPI的溫度測試誤差約為±2.0%，而FT-SPI的溫度測試誤差優(yōu)于±1.0%。FT-SPI方法優(yōu)于傳統(tǒng)方法，同時，采用基于圖像識別的IR-SPI方法得到的數(shù)據(jù)中存在個別偏差較大的測試點，分析其產(chǎn)生的原因是由于被測表面存在細紋及漆面磨損造成的，在采用光滑鋼板實驗時此現(xiàn)象消失，故表面的瑕疵在IR-SPI檢測方法下會導致溫度測試粗大誤差的出現(xiàn)。而采用FT-SPI方法由于其工作機理是整體條紋的頻譜分析，故不受個別點位的影響，所以被測表面的瑕疵、漆面磨損等對溫度檢測的影響可以忽略。在此基礎上，由于測試過程中的炮體材料基本均勻，不存在明顯溫度畸變點，故采用溫度標定的方法可以進一步提高系統(tǒng)溫度測試的準確度。標定方法是將標準溫度在相應測試區(qū)間中選取多個參考點，以此對FT-SPI方法得到的數(shù)據(jù)進行校正，得到校正系數(shù)，從而引入下一次測試中，由此系統(tǒng)的溫度測試誤差可達到±0.5%。

4 結(jié) 論

由于目前瞬態(tài)測試設備主要采用單點溫度檢測，不能對被測表面整體范圍監(jiān)測，并且熱慣性較大，瞬態(tài)變化過程難于捕獲，故設計開發(fā)了基于散斑干涉與光譜頻域分析相結(jié)合的瞬態(tài)溫度分布測試系統(tǒng)，該工作機理不受熱慣性影響，其測溫時間僅取決于系統(tǒng)散斑干涉條紋的采樣速度，大大提高了其瞬態(tài)測溫能力。與此同時，采用傅氏變換的頻譜分析技術(shù)使系統(tǒng)克服由表面瑕疵、漆面磨損等問題造成的粗大誤差。綜上所述，本系統(tǒng)可以實時、穩(wěn)定地獲取炮體瞬態(tài)溫度分布，符合設計要求。

[1] Kim Kyeong-Suk, Chang Ho-Seob, Naseem Akhter. International Journal of Precision Engineering and Manufacturing, 2010, 11(6): 979.

[2] DENG Xiang-yang, LIU Shou-xian, PENG Qi-xian(鄧向陽, 劉壽先, 彭其先). Explosion and Shock Waves(爆炸與沖擊波)，2015, 35(2): 191.

[3] Bingleman L W, Schajer G S. Experimental Mechanics，2011, 51: 1207.

[4] LI Xiang-yu, HUANG Zhan-hua, ZHU Meng(李翔宇, 黃戰(zhàn)華, 朱 猛). Acta Optica Sinica(光學學報)，2015, 39(6): 129.

[5] CAI Chang-qing, ZHANG Yong-shan, WANG Da-yang(蔡長青, 張永山, 汪大洋). Laser & Optoelectronics Progress(激光與光電子學進展)，2014, 51(11): 141.

[6] CHEN Man-tang, LIU Yun-peng, LIU Yang(陳滿堂, 劉云朋, 劉 楊). Fire Control & Command Control(火力與指揮控制)，2014, 39(10): 179.

[7] WU Jin-hui, YANG Rui-feng, WANG Gao(武錦輝, 楊瑞峰, 王 高). Opto-Electronic Engineering(光電工程)，2012, 39(9): 132.

[8] Gaetano Restivo, Gabriel A Isaicu, Gary L Cloud. J. Nondestruct Eval., 2008, 27: 135.

[9] YUAN Zong-heng, SONG Mei-jie, XIONG Xian-ming(袁縱橫, 宋美杰, 熊顯名). Journal of Optoelectronics·Laser(光電子·激光)，2011, 22(3): 28.

[10] ZHANG Jia-ming, WANG Wen-rui, NIE Shuai(張佳明, 王文瑞, 聶 帥). China Measurement & Test(中國測試)，2014, 40(5): 39.

[11] Ruiz P D, Huntley J M, Wildman R D. Appl. Opt., 2005, 44(19): 3945.

[12] TANG Cheng, REN Hong-wei, CHEN Xia(唐 晨, 任宏偉, 陳 霞). Laser Optoelectronics Progress(激光與光電子學進展)，2010, 47: 021201-1.

Testing Research of Transient Temperature Distribution for the Barrel Surface by Speckle Pattern Interferometry

LANG Wen-jie1, CHEN Guo-guang2, TIAN Xiao-li2, XIN Chang-fan2

1. School of Information and Communication Engineering, North University of China, Taiyuan 030051, China 2. School of Mechatronics Engineering, North University of China, Taiyuan 030051, China

There are some problems in the traditional transient temperature test equipment. The thermal inertia is great, and can only be a single point of detection. To be able to achieve real-time monitoring for transient temperature distribution change of the gun body surface, the test system for transient temperature distribution was designed based on Speckle Pattern Interferometry (SPI) and spectroscopy. In the system, transient temperature change of the barrel led to slight deformation, and it was converted into speckle interference fringes by SPI technology. Spectral distribution function was obtained by the interference fringes by the Fourier transform, so the information of interference fringe deformation was incorporated into the frequency domain. The data of temperature distribution can be inverted on any sampling time by spectral distribution function. In experiments, the ZX-FB1 fiber optic thermometer was used to test transient temperature on a single point as the standard value. The center wavelength of the laser was 555 nm, and the speckle pattern interference fringes were collected by area array CCD. Image Recognition-Speckle Pattern Interferometry (IR-SPI) and Fourier Transform-Speckle Pattern Interferometry (FT-SPI) were used in experiments, the calculation of transient temperature was completed through two methods. Experimental results are that both methods can achieve transient temperature detection. But the FT-SPI is higher in terms of accuracy, and it can effectively overcome the gross error caused by the surface defects, paint wear and other similar problems.

Transient temperature distribution; SPI; Spectral analysis; Fourier transform

Aug. 12, 2015; accepted Dec. 21, 2015)

2015-08-12，

2015-12-21

國家自然科學基金項目(51305409)和兵科院縱向課題(JZ2009100)資助

郎文杰，1977年生，中北大學信息與通信工程學院副教授 e-mail: lwjnuc@163.com

O432

A

10.3964/j.issn.1000-0593(2016)02-0355-04