汪思國,湯更生 ,楊 輝 ,張 龍,張 俊
(中國空氣動力研究與發(fā)展中心,四川綿陽 621000)
風(fēng)洞試驗中旋轉(zhuǎn)模型的掃描角和轉(zhuǎn)速定量實時測量一直是比較棘手的問題。國內(nèi)風(fēng)洞常規(guī)的測量手段有兩種:一種是在模型內(nèi)部安裝轉(zhuǎn)速記錄儀,這種方式不僅會改變彈體的實際結(jié)構(gòu)和重量,破壞模型的真實性和完整性,而且難以測量模型的瞬態(tài)轉(zhuǎn)速和掃描角,測量精度和準(zhǔn)確度也不高;另一種是雙目圖像測量方法,這種測量方式不改變模型整體結(jié)構(gòu),只需要在模型表面粘貼標(biāo)記物,但在圖像處理過程中高速旋轉(zhuǎn)體的特征點難以提取,且要進(jìn)行繁瑣的坐標(biāo)變換,測量誤差也較大。
投影圖像測量方法采用非接觸測量方式,風(fēng)洞試驗中不需要在模型內(nèi)部安裝測量器件和在表面粘貼任何標(biāo)識物,有效保證了試驗?zāi)P偷恼鎸嵭院屯暾?減少了試驗過程的繁瑣準(zhǔn)備工作;圖像處理過程方便快捷,有效保證了測量結(jié)果的準(zhǔn)確度和精度。
根據(jù)西安某所風(fēng)洞試驗傘-彈系統(tǒng)結(jié)構(gòu)的特殊性和復(fù)雜性,以及試驗過程中彈體的快速旋轉(zhuǎn)等特點,使用投影圖像處理這種非接觸測量方法穩(wěn)妥可靠。為此,在實驗室開展了投影圖像法測量研究,對傘-彈模型進(jìn)行了模擬測量,并首次將這種方法應(yīng)用于國內(nèi)立式風(fēng)洞試驗測量上,取得了良好的測量結(jié)果,測量精度達(dá)到了用戶的要求。
彈體的掃描角β是彈體的軸線與風(fēng)速方向(或者設(shè)備裝置軸線)的夾角。在傘-彈系統(tǒng)中,由于彈體是一個回轉(zhuǎn)體,因此,彈體的掃描角實際上就是彈體的母線與風(fēng)速方向(或者設(shè)備裝置軸線)的夾角。圖1給出了傘-彈系統(tǒng)掃描角的定義。
圖1 傘-彈系統(tǒng)掃描角定義Fig.1 Defining for sweep angle of parachute-ball system
由于彈體是一個回轉(zhuǎn)體,因此,彈體的掃描角β實際上就是彈體的母線AC1與設(shè)備裝置軸線 AA1的夾角。在圖 2所示的長方體 ABCD-A1B1C1D1中,假設(shè)長方體的對角線 AC1是彈體的母線或母線的一部分,AA1是風(fēng)速方向或設(shè)備裝置軸線,則立體角∠A1AC1就是彈體的掃描角。令A(yù)A1=1,則有:
彈體的掃描角β可以分解到兩個正交平面內(nèi),分量分別為β1和β2。只要測量出掃描角在兩個正交平面內(nèi)的角分量β1和β2,很容易計算出彈體的掃描角β。
圖2 掃描角分解原理圖Fig.2 Principle of analyzing sweep angle
由于掃描角可以分解到兩個正交平面內(nèi),所以對彈體的掃描角測量,可以采用光學(xué)投影照相測量法,光學(xué)投影和照相的目標(biāo)就是彈體本身,圖3所示就是光學(xué)投影照相測量法的原理圖。具體方法是用兩束相互正交的平行光束,分別將β1和β2所在的兩個正交平面投影到立式風(fēng)洞試驗段后面一定位置兩個正交的投影屏PY和PX上,用兩個相互正交的CCD成像系統(tǒng)將投影圖像記錄下來,通過實時分析處理得到角分量 β1和 β2的投影角 βY和βX。
對于投影角βY和βX的測量,必須在試驗前用設(shè)備裝置的軸線對兩個CCD成像系統(tǒng)分別進(jìn)行定標(biāo),以這一軸線作為測量基準(zhǔn),在投影圖像的采集過程中必須保持兩個CCD在上下和前后位置不動且不產(chǎn)生偏擺。從投影屏PY上的投影圖像進(jìn)行投影角βY的計算時,只需要根據(jù)投影屏PY上投影圖像中彈體母線上任意兩點的坐標(biāo)計算出投影圖像的斜率,根據(jù)斜率取余角就可以得到投影角βY,同理也可以得到投影角 βX。
對彈體旋轉(zhuǎn)角θ的計算:測量出投影屏PX和PY上投影圖像中對應(yīng)的投影角βX和βY后,可以知道掃描角β在兩個正交平面內(nèi)的角分量β1和β2,根據(jù)幾何關(guān)系可以計算出彈體的旋轉(zhuǎn)角θ。在圖2所示的長方體ABCD-A1B1C1D1中:
根據(jù)兩個高速高分辨率CCD記錄到的投影圖像,可以判斷出旋轉(zhuǎn)角θ的象限范圍,從而實現(xiàn)對旋轉(zhuǎn)角θ在360°旋轉(zhuǎn)范圍內(nèi)的角度測量。將每一時刻的旋轉(zhuǎn)角θ處理并記錄下來,再根據(jù)兩個高速高分辨率CCD的同步攝像幀頻,可以得到相鄰兩個時刻的旋轉(zhuǎn)角θ1和 θ2的時間間隔 ΔT,并且在每一次試驗中的ΔT是固定的,那么傘彈旋轉(zhuǎn)的瞬時角速度為:
圖3 掃描角測試系統(tǒng)圖Fig.3 Principle of testing sweep angle
根據(jù)上面介紹的掃描角測試方法,設(shè)計了投影圖像法測量旋轉(zhuǎn)體掃描角的測量系統(tǒng),在實驗廠房進(jìn)行了模擬測量,該系統(tǒng)由兩路正交且相同的測量系統(tǒng)構(gòu)成,每一路測量系統(tǒng)包括投影系統(tǒng)、投影屏、圖像采集系統(tǒng)和運動控制系統(tǒng)等分系統(tǒng)組成。2007年10月采用該系統(tǒng)在西安某所立式風(fēng)洞中對某旋轉(zhuǎn)彈體實驗進(jìn)行了彈體的掃描角和轉(zhuǎn)速測量,并得到了良好的測試結(jié)果。圖4給出了系統(tǒng)設(shè)計的三維效果圖。
圖4 掃描角測試系統(tǒng)三維效果圖Fig.4 Three-dimensional effect of testing sweep angle
根據(jù)系統(tǒng)中高速高分辨率面陣CCD的分辨率(1280×1024)和掃描角測量分辨率(0.1°)的要求,結(jié)合現(xiàn)場試驗距離,計算出CCD的拍攝范圍為893.6mm×714.9mm,于是選定兩路光源截面尺寸為900mm×750mm的矩形投影平行光束作為投影光源。
采用2路高分辨率CCD(1280×1024)對傘-彈系統(tǒng)彈體在投影屏上的投影圖像進(jìn)行正交同步圖像拍攝,兩個水平方向上的拍攝區(qū)域大小為893.6mm左右,確保角度測量精度達(dá)到規(guī)定要求,CCD所拍攝的圖像經(jīng)支持64位總線、采集速度能達(dá)到74幀/秒的圖像采集卡采集,存放于當(dāng)?shù)赜嬎銠C內(nèi)存(8G)中,采集結(jié)束后兩路圖像數(shù)據(jù)經(jīng)光纖網(wǎng)絡(luò)傳輸?shù)礁咚賵D像處理服務(wù)器上進(jìn)行保存并進(jìn)行圖像處理,給出掃描角和角速度計算結(jié)果。掃描角測量必須由同一時刻兩個正交方向的圖像數(shù)據(jù)計算得出,這對兩路CCD采集同步要求較高,同步觸發(fā)采集采用美國Epix公司T TL同步調(diào)理模塊控制(同步時間達(dá)到50ns),達(dá)到兩路圖像采集系統(tǒng)采集精確同步要求。
由于CCD在兩個水平方向上所拍攝的有效范圍為893.6mm,當(dāng)傘-彈系統(tǒng)完全飄離CCD視場時,高分辨率CCD、投影光源必須移動相應(yīng)位置,確保圖像在CCD視場范圍內(nèi)。系統(tǒng)采用兩路低分辨率面陣CCD在兩個方向3.2m范圍進(jìn)行圖像拍攝,從而可確定傘-彈系統(tǒng)的當(dāng)前位置,位置工控機可以通過電機控制卡控制安裝在軌道上的圖像采集CCD、投影光源在水平方向左右同步平移適當(dāng)位置,使傘-彈系統(tǒng)進(jìn)入CCD拍攝視場范圍內(nèi)以及確保投影光源能在拍攝位置進(jìn)行投影。
采用該系統(tǒng)在西安某所立式風(fēng)洞對某傘彈系統(tǒng)進(jìn)行了彈體的掃描角和轉(zhuǎn)速測量。圖像采集方式為外同步觸發(fā)方式,采集容量為500幅圖像,曝光時間為0.8ms,試驗大廳兩臺圖像采集機數(shù)據(jù)通過光纖網(wǎng)絡(luò)存儲到測控間服務(wù)器的硬盤上,利用圖像處理和角度計算軟件,對兩路圖像分別進(jìn)行角度序列計算,然后進(jìn)行兩路角度合成,得到彈體的掃描角和轉(zhuǎn)速曲線。以100幅圖像(1.3513s)為計算時長,得到的轉(zhuǎn)速測試結(jié)果和掃描角測試結(jié)果(圖5,6)。測試結(jié)果表明:圖像處理系統(tǒng)具備轉(zhuǎn)速和掃描角測量功能,能夠可靠地給出穩(wěn)態(tài)掃描試驗過程中傘彈系統(tǒng)的掃描角曲線和彈體旋轉(zhuǎn)的角速度。
圖5 穩(wěn)態(tài)掃描試驗中傘彈轉(zhuǎn)速隨時間變化曲線Fig.5 Parachute-ball rotator rotating speed-time curve
圖6 穩(wěn)態(tài)掃描試驗中傘彈掃描角隨時間變化曲線Fig.6 Parachute-ball rotator sweep angle-time curve
該系統(tǒng)誤差主要來源于系統(tǒng)的安裝誤差和圖像邊緣判讀誤差。系統(tǒng)的安裝誤差包括兩路正交系統(tǒng)的夾角誤差、各光源的偏心差以及CCD的安裝誤差。對于夾角誤差,在對兩路掃描角分量進(jìn)行掃描角合成時,對系統(tǒng)的夾角誤差進(jìn)行了修正;對于各光源的偏心差和CCD的安裝誤差,采用零角度懸掛靜態(tài)模型試驗,對每一個小口徑光源引起的變形誤差和CCD安裝進(jìn)行了補償修正;對于圖像邊緣判讀誤差,判讀過程中圖像邊緣放寬到2個像素點(判讀誤差0.12°),有效克服了投影圖像邊緣模糊對特征點提取的影響。
風(fēng)洞試驗中測量傘彈的轉(zhuǎn)速時,必須要獲得相鄰兩個時刻的模型旋轉(zhuǎn)角,而模型的旋轉(zhuǎn)角直接由該時刻的兩個掃描角分量決定,所以掃描角測量誤差直接導(dǎo)致了傘彈的轉(zhuǎn)速計算誤差。計算轉(zhuǎn)速時取的兩個時刻的時間間隔越長,轉(zhuǎn)速計算誤差就越小,這時的轉(zhuǎn)速就是該時間段內(nèi)的平均轉(zhuǎn)速。一般而言,轉(zhuǎn)速計算誤差的直接來源是由于某一時刻掃描角分量的誤差導(dǎo)致該時刻模型旋轉(zhuǎn)角的計算誤差,所以采用較長時段計算轉(zhuǎn)速的誤差很小,計算時段太短引起的計算誤差相對稍大一些。
(1)半導(dǎo)體激光器安裝過程中采用準(zhǔn)值鏡校準(zhǔn)方法,減小了矩形光源中各半導(dǎo)體激光器的光軸偏差,從而使模型投影圖像變形得到了較大程度的改善;
(2)通過在試驗段位置懸掛0°模型,分別對各方格光源進(jìn)行了角度測量修正。處理過程中軟件自動調(diào)用相應(yīng)方格的修正值,消除了由于CCD的安裝誤差、光源安裝誤差及模型投影圖像變形帶來的角度測量誤差;
(3)圖像處理過程中,通過彈體邊界特征點 x,y的坐標(biāo),判斷出角度符號的正負(fù);利用每一時刻兩個掃描角分量的角度的符號變化,判斷出每一時刻掃描角的象限角度,從而準(zhǔn)確測量出了傘彈系統(tǒng)一定周期內(nèi)的平均轉(zhuǎn)速和瞬態(tài)轉(zhuǎn)速;
(4)圖像處理過程中,采用線性回歸算法,大大提高了投影圖像彈體邊界特征點提取的準(zhǔn)確性,排除彈體分段后以及因為各投影矩形光源不平行引起的模型變形對邊界特征點提取的干擾,大大提高了掃描角計算的精度;
(5)圖像處理軟件通過灰度閾值動態(tài)調(diào)節(jié),有效解決了圖像采集過程中外界光對投影圖像的干擾,確保彈體邊界特征點的有效提取;通過對模型尺寸的動態(tài)定義,配合圖像放大倍數(shù),可以計算不同尺寸不同型號彈體的掃描角度和轉(zhuǎn)速。
投影圖像法測量系統(tǒng)采用了非接觸測量方式,提高了測量精度和準(zhǔn)確度,并將其應(yīng)用在西安某所立式風(fēng)洞進(jìn)行傘彈風(fēng)洞試驗。經(jīng)對多種型號的旋轉(zhuǎn)彈體進(jìn)行了掃描角和轉(zhuǎn)速的動態(tài)測量的試驗驗證,結(jié)果表明:對于風(fēng)洞試驗中快速旋轉(zhuǎn)模型的轉(zhuǎn)速和角度測量來說,投影圖像法測量手段先進(jìn),穩(wěn)妥可靠,有效的保證了試驗?zāi)P偷耐暾院驼鎸嵭?圖像處理系統(tǒng)具備轉(zhuǎn)速和掃描角測量功能,能夠可靠地給出穩(wěn)態(tài)掃描試驗過程中傘彈系統(tǒng)的掃描角曲線和彈體旋轉(zhuǎn)的角速度。
[1]GONZALEZ R C,WOODS R E.Digital image processing[M].2nd end.,Prentice Hall,Upper Saddle River,NJ.2002.
[2]Image Processing Toolbox,Users Guide,Version 4[Z].The MathWorks,Inc.Natick,MA.2003.
[3]HARTLEY R,ZISSERMAN A.計算機視覺中的多視圖幾何[M].合肥:安徽大學(xué)出版社,2002.
[4]劉政崇.穩(wěn)態(tài)掃描模擬實驗裝置研制技術(shù)總結(jié)報告[R].四川綿陽:中國空氣動力研究與發(fā)展中心設(shè)備設(shè)計及測試技術(shù)研究所,2007.
[5]馬頌德,張正友.計算機視覺-計算機理論與算法基礎(chǔ)[M].北京:科學(xué)出版社,1998,56-59.