尚金奎 ,馬曉光

(中國(guó)航空工業(yè)空氣動(dòng)力研究院,沈陽(yáng) 110034)

0 引 言

風(fēng)洞壓敏涂料技術(shù) PSP(Pressure Sensitive Paint)利用光學(xué)技術(shù)實(shí)現(xiàn)風(fēng)洞模型(試件)表面的壓力分布測(cè)量,它的基本原理是探針?lè)肿?壓敏涂料中的活性成份)的氧熒光猝滅效應(yīng)。它是一種“非介入式”的氣動(dòng)試驗(yàn)技術(shù),不需要對(duì)試件進(jìn)行特別加工,對(duì)風(fēng)洞流場(chǎng)不會(huì)產(chǎn)生影響,可以得到高品質(zhì)的試驗(yàn)數(shù)據(jù),并且由于PSP技術(shù)采用大面積噴涂的作業(yè)方式,所以能夠獲取整個(gè)試件表面的連續(xù)壓力分布數(shù)據(jù)。

在航空航天技術(shù)發(fā)達(dá)的美國(guó)和俄羅斯(前蘇聯(lián))等國(guó)家,早在上世紀(jì)80～90年代,PSP技術(shù)就已經(jīng)取得了突破性進(jìn)展,并應(yīng)用于型號(hào)風(fēng)洞測(cè)壓試驗(yàn)[2-3],在20世紀(jì)90年代末和21世紀(jì)初進(jìn)一步發(fā)展達(dá)到了比較成熟的狀態(tài),開(kāi)始了大規(guī)模的生產(chǎn)型應(yīng)用[4-6]。我國(guó)在這方面的研究雖然起步相對(duì)較晚,但是經(jīng)過(guò)不斷的努力,目前已經(jīng)初步具備了進(jìn)行生產(chǎn)型應(yīng)用的能力[7-8]。

PSP技術(shù)獲取數(shù)據(jù)的處理和修正工作直接關(guān)系到試驗(yàn)成敗和最終精度。國(guó)外相關(guān)方面的研究開(kāi)展較早且深入,技術(shù)相對(duì)成熟[5,6,9]。相對(duì)而言,國(guó)內(nèi)數(shù)據(jù)處理手段和方式較為單一,目前國(guó)內(nèi)主要是通過(guò)校準(zhǔn)箱獲取的校準(zhǔn)曲線得到試件表面壓力分布。筆者提出了一種試驗(yàn)圖像與中值濾波技術(shù)相結(jié)合的數(shù)據(jù)處理方法,在不使用校準(zhǔn)箱的情況下得到試件表面的壓力(Cp值)分布。

1 基本知識(shí)

1.1 壓敏涂料相關(guān)知識(shí)

壓敏涂料是PSP技術(shù)中的核心部分,由活性層和基層組成,其中的活性層是由發(fā)光體(熒光探針?lè)肿?和基質(zhì)構(gòu)成,基質(zhì)除了能固定發(fā)光體外,要求具有較高的氣體透過(guò)性,對(duì)熒光探針?lè)肿拥墓馕蘸桶l(fā)射無(wú)干擾,在較寬溫度范圍內(nèi)能夠彈性回復(fù)以及優(yōu)良的成膜性[7]。氧熒光猝滅效應(yīng)就發(fā)生在活性層中。

氧熒光猝滅效應(yīng)的強(qiáng)弱與環(huán)境中的氧分子濃度,探針?lè)肿蛹ぐl(fā)態(tài)的本征壽命以及激發(fā)態(tài)猝滅速率有關(guān),通?？梢杂肧tern-Volmer公式表述壓敏涂料氧熒光猝滅效應(yīng)的影響和氧分子濃度的關(guān)系[8]：

式中,I0和I分別表示無(wú)氧環(huán)境下和在氧濃度為[O2]環(huán)境下壓敏涂料的發(fā)光強(qiáng)度;kq表示雙分子碰撞猝滅速率常數(shù);τ0表示激發(fā)態(tài)的本征壽命;[O2]表示猝滅劑氧分子的濃度。

在風(fēng)洞試驗(yàn)中,要得到I0既不現(xiàn)實(shí),也不可能,因此引入?yún)⒄展鈴?qiáng) Iref,Iref不受氧猝滅的影響,同時(shí)氧氣濃度與壓力有正比關(guān)系,所以通常在風(fēng)洞試驗(yàn)中采用Stern-Volmer公式的另一種形式[10]

其中,Iref、Pref,表示在選定參照條件的溫度下,壓敏涂料的參照光強(qiáng)和參考?jí)毫?A(T)和B(T),通常稱為Stern-Volmer系數(shù)。

在不考慮溫度影響,且采用壓力系數(shù)Cp=2(pp∞)/p∞的情況下,公式(2)可以寫成如下形式 ,

1.2 圖像中值濾波

一幅圖像的灰度值或亮度值可以看作一個(gè)二維函數(shù)f(x,y)[11],x,y分別表示圖像平面的橫坐標(biāo)和縱坐標(biāo)。一般情況下,信號(hào)在圖像中屬于低頻成份,噪聲在圖像中屬于高頻成份,所以可以利用它們之間特性的差異對(duì)圖像進(jìn)行適當(dāng)處理,消除或減少圖像中的噪聲。

減少圖像噪聲的方法在空間域中又稱為圖像的平滑空間濾波,圖像中值濾波就是其中的一種,它是典型的非線性低通濾波器。顧名思義,所謂圖像中值濾波就是將圖像中各點(diǎn)的值用其鄰域的中值代替,即,若用矩陣[xij]表示一幅數(shù)字化的圖像(i,j表示圖像像素坐標(biāo)),W[xij]表示對(duì)圖像[xij]中的點(diǎn) xij做窗口操作,med(W[xij])表示對(duì)窗口內(nèi)的所有點(diǎn)取中值,用[yij]表示經(jīng)中值濾波后的圖像,那么[12]

由于圖像中值濾波能夠有效的去除圖像中的噪聲,并且能夠克服線性濾波對(duì)圖像造成的模糊,所以得到了廣泛的應(yīng)用。

2 試驗(yàn)圖像的處理

試驗(yàn)采用的試件是半翼展機(jī)翼模型,試驗(yàn)馬赫數(shù)為0.8,迎角為12°,試驗(yàn)采用的CCD相機(jī)分辨率為512×512,12位灰度等級(jí)。

在半翼展試件上表面共9行測(cè)壓孔,可以通過(guò)常規(guī)測(cè)壓技術(shù)獲得測(cè)壓孔當(dāng)?shù)氐腃p值,其中靠近翼根的為第一行,共10個(gè)測(cè)壓孔,其余測(cè)壓孔個(gè)數(shù)沿翼展方向依次為 9、8 、8、7 、6、6、5、3,靠近前緣的測(cè)壓孔為每行測(cè)壓孔中的第一個(gè),靠近后緣的為最后一個(gè)。試件表面除測(cè)壓孔外,還有標(biāo)志點(diǎn),標(biāo)志點(diǎn)在進(jìn)行多幅圖像配準(zhǔn)時(shí)使用。試件原始光強(qiáng)比圖像和MATLAB中值濾波后的圖像如圖1、2所示。

圖1 半翼展原光強(qiáng)比灰度圖像Fig.1 Original light intensity ratio image

圖2 濾波后的光強(qiáng)比圖像Fig.2 Filtered light intensity ratio image

在圖1中,試驗(yàn)獲得的光強(qiáng)比圖像中存在測(cè)壓孔(如圖中以紅色菱形框框定的位置)和標(biāo)志點(diǎn)(如圖中以藍(lán)色圓形圈定的位置),這兩種位置點(diǎn)在參考圖像和試驗(yàn)圖像中均應(yīng)表現(xiàn)為黑色(像素灰度值為0,即沒(méi)有熒光產(chǎn)生),但由于干擾和噪聲的存在,在表現(xiàn)為光強(qiáng)比(試驗(yàn)圖像和參考圖像光強(qiáng)之比)的圖1中測(cè)壓孔和標(biāo)志點(diǎn)處本身像素灰度值波動(dòng)很大,與其鄰域灰度值差別也很大,這與其周圍信號(hào)像素灰度值變化比較平緩的特點(diǎn)形成了鮮明的對(duì)比。并且由于測(cè)壓孔占用像素較少,距離測(cè)壓孔1～2個(gè)像素的點(diǎn)其灰度與周圍像素有較大差異,但距離為2～3個(gè)像素的地方,其鄰域像素受到的干擾較小,所以能夠?qū)⑵淇醋鲌D像噪聲,從而可以應(yīng)用圖像濾波的方法將測(cè)壓孔所呈現(xiàn)的噪聲在圖像上濾除,恢復(fù)測(cè)壓孔處的光強(qiáng)比數(shù)據(jù)。

圖3 圖像橫坐標(biāo)x-Cp值關(guān)系曲線圖Fig.3 x-Cprelational graphs

為了得到PSP圖像中測(cè)壓孔處的光強(qiáng)比數(shù)據(jù),結(jié)合測(cè)壓孔處像素灰度值的變化特點(diǎn),采用的濾波器為5×5的空心中值濾波器。為了盡量獲得最接近真實(shí)的數(shù)據(jù),在計(jì)算式(3)中系數(shù)m和n時(shí),應(yīng)該小心的選取參考點(diǎn),可以根據(jù)試件的空氣動(dòng)力學(xué)特性知識(shí)以及風(fēng)洞試驗(yàn)的經(jīng)驗(yàn)積累在風(fēng)洞試驗(yàn)前事先確定,也可以在試驗(yàn)完成后通過(guò)分析PSP光強(qiáng)比圖像確定。本文選取的參考點(diǎn)是圖1中由菱形框框定的3個(gè)測(cè)壓孔。

濾波后的光強(qiáng)比圖像見(jiàn)圖2,從圖2中可以看到空心中值濾波很好地去除了測(cè)壓孔處的"噪聲",對(duì)標(biāo)志點(diǎn)處灰度值也有一定的"恢復(fù)"作用,并且試件邊緣及其構(gòu)件邊緣以及光強(qiáng)比數(shù)據(jù)發(fā)生較大變化形成的邊界被較好地保留了下來(lái)。

3 試驗(yàn)結(jié)果與分析

以圖1中所示的3個(gè)測(cè)壓孔為參考點(diǎn),結(jié)合濾波后得到的參考點(diǎn)光強(qiáng)比數(shù)據(jù)和常規(guī)測(cè)壓技術(shù)實(shí)測(cè)Cp值,采用直線擬合方式得到的直線為

其中,m=6.1449,n=-5.2604。公式(5)就是光強(qiáng)比與Cp值關(guān)系曲線。

為了更加清晰的表示PSP試驗(yàn)方法和常規(guī)測(cè)壓方法所得的Cp值之間的關(guān)系,下面列出了Cp值與圖像橫坐標(biāo)曲線關(guān)系圖,見(jiàn)圖3。圖中PSP技術(shù)數(shù)據(jù)曲線是一行測(cè)壓孔中相鄰測(cè)壓孔之間連線所在位置的PSP Cp值與其所在位置橫坐標(biāo)x形成的曲線。

圖3中,第1、2和3行測(cè)壓孔處PSP方法獲得的Cp值形成的曲線在測(cè)壓孔處與常規(guī)測(cè)壓方法所獲得的測(cè)壓孔處Cp值吻合較好。第4、5行測(cè)壓孔處,PSP方法獲得的Cp值在機(jī)翼前緣與常規(guī)測(cè)壓方法所獲得的數(shù)據(jù)相比明顯偏小,其中第4行第4個(gè)測(cè)壓孔附近PSP方法獲得的數(shù)據(jù)形成的曲線突然向上突起,與常規(guī)測(cè)壓方法獲得的Cp值形成較大的差值,在第5個(gè)測(cè)壓孔處兩種方法獲得的Cp值又很好的吻合,這個(gè)突起的存在不符合空氣動(dòng)力學(xué)常識(shí),從原始光強(qiáng)比圖像來(lái)看(圖1),此處,局部圖像灰度突然變小,形成了一個(gè)像素灰度相對(duì)不連續(xù)區(qū),由此可以推斷這一區(qū)域的原始光強(qiáng)比數(shù)據(jù)存在錯(cuò)誤,而造成數(shù)據(jù)錯(cuò)誤的原因可能是多方面的,比如壓敏涂料噴涂后的局部污染、局部區(qū)域壓敏涂料噴涂質(zhì)量較差以及不恰當(dāng)?shù)膲好敉苛蠂娡亢筇幚淼?。?、7、8和9行測(cè)壓孔處,雖然在第6、7行測(cè)壓孔機(jī)翼前緣位置兩種方法的偏差相對(duì)于其后的其他位置較大,但整體來(lái)說(shuō)這4行測(cè)壓孔處PSP數(shù)據(jù)與常規(guī)測(cè)壓方法測(cè)得的Cp值吻合較好。

從以上的分析中可以看出,除去個(gè)別錯(cuò)誤的數(shù)據(jù)和局部的機(jī)翼前緣外,該方法處理所得的PSP數(shù)據(jù)與常規(guī)測(cè)壓方法能夠較好的吻合。至于局部機(jī)翼前緣位置偏離較大的原因則仍需要進(jìn)行進(jìn)一步的研究,采用校準(zhǔn)箱校準(zhǔn)曲線的方法獲得的PSP數(shù)據(jù)同樣存在這樣的問(wèn)題。

4 結(jié)論與展望

目前,國(guó)內(nèi)通常采用校準(zhǔn)箱校準(zhǔn)曲線的方法得到試件壓力數(shù)據(jù)。筆者提出的方法結(jié)合傳統(tǒng)測(cè)壓孔得到的壓力數(shù)據(jù)和圖像中值濾波技術(shù)獲得試件表面的壓力分布,加入了風(fēng)洞試驗(yàn)中試件自身因素對(duì)試驗(yàn)結(jié)果的影響,可以在一定程度上彌補(bǔ)校準(zhǔn)箱方法的缺陷,提高試驗(yàn)精度,具有較大的發(fā)展?jié)摿蛻?yīng)用前景。

相對(duì)擁有先進(jìn)航空技術(shù)的歐美等國(guó)家而言,PSP技術(shù)在我國(guó)發(fā)展還不夠成熟,存在大量有待解決的問(wèn)題,需要進(jìn)行更加深入的研究,但是PSP技術(shù)由于不需要對(duì)試件進(jìn)行特別的加工,而且能夠獲得試件表面的壓力場(chǎng)數(shù)據(jù),具有天然的優(yōu)勢(shì),是風(fēng)洞試驗(yàn)技術(shù)的發(fā)展方向之一,只有真正的掌握它才能適應(yīng)未來(lái)飛行器設(shè)計(jì)的需要,為未來(lái)飛行器的設(shè)計(jì)提供有力的支持。

[1]ENGLER R H,K LEIN C.360°PSP measurements in transonic flow[C].Proceedings of the 19th International Congress on Instrumentation in Aerospace Simulation Facilities,149-158,2001.

[2]MCLACHLAN B G,BELL J H.Pressure-sensitive paint in aerodynamic testing[J].Experimental Thermal and Fluid Science,10(4)：470-485,1995.

[3]BOROVOY V,BYKOV A,MOSHAROV V,et al.Pressure sensitive paint application in shock wind tunnel[C].International Congress on Instrumentation in Aerospace Simulation Facilities,'95 Record,34.1-34.4,1995.

[4]BELL J H,BURNER A W.Data fusion in wind tunnel testing;Combined pressure paint and model deformation measurements[C].20th Advanced Measurement and Ground Testing Technology Conference,AIAA-98-2500,1998.

[5]ERICKSON G E,GONZALEZ H A.Pressure-sensitive paint investigation of double-delta wing vortex flow manipulation[C].43rd AIAA Aerospace Sciences Meeting&Exhibit.AIAA 2005-1059,,2005.

[6]LIU Tian-shu,GUILLE M,SULLIVAN J P.Accuracy of pressure sensitive paint[J].AIAA Journal,39(1),2001.

[7]程厚梅,張永存,陳柳生.光學(xué)涂層測(cè)壓技術(shù)[M].空氣動(dòng)力學(xué)前沿研究,北京：中國(guó)宇航出版社,108-115,2003.

[8]張永存.在高速風(fēng)洞中使用熒光壓力傳感器技術(shù)對(duì)飛機(jī)模型壓力場(chǎng)的試驗(yàn)研究[J].流體力學(xué)試驗(yàn)與測(cè)量,2000,14(2)：21-25.

[9]KAMMEYER M,DONOVAN J,KELBLE C,et al.Accuracy assessment of a pressure-sensitive paint measurement system[C].40th Aerospace Sciences Meeting&Exhibit,AIAA 2002-0530,2002.

[10]WATKINS A N,JORDAN J D,LEIGHT Y B D,et al.Development of next generation lifetime PSP imaging systems[C].Proceedings of the 20th International Congress on Instrumentation in Aerospace Simulation Facilities,2003.

[11]GONZALEZ R C,WOODS R E,EDDINS S L.數(shù)字圖像處理(MATLAB版)[M].北京：電子工業(yè)出版社,2005.

[12]邢藏菊,王守覺(jué),鄧浩江等.一種基于極值中值的新型濾波算法[J].中國(guó)圖象圖形學(xué)報(bào),6(A)：533-536,2001.