閆萬方，楊 輝，尼文斌，賀麗慧

(中國(guó)航天空氣動(dòng)力技術(shù)研究院，北京 100074)

現(xiàn)代高空高速飛行器具有跨空域、寬速域、大機(jī)動(dòng)飛行的特點(diǎn)。傳統(tǒng)探出式大氣數(shù)據(jù)系統(tǒng)已難以滿足其飛行控制需要[1-3]，主要原因在于：① 高速飛行時(shí)，突出裝置難以適應(yīng)頭部的高溫，同時(shí)，與周圍大氣形成的激波干擾將嚴(yán)重影響飛行器氣動(dòng)性能；② 大機(jī)動(dòng)飛行時(shí)，前端突出裝置已成為引起頭部渦及橫側(cè)向失穩(wěn)的主要因素之一；③ 突出飛行器表面部件極不利于飛行器隱身。嵌入式大氣數(shù)據(jù)傳感(Flush Air Data Sensing,FADS)系統(tǒng)具有氣動(dòng)外形影響小、寬速域/大機(jī)動(dòng)及惡劣飛行環(huán)境適應(yīng)能力強(qiáng)的特點(diǎn)，其在可靠性、穩(wěn)定性方面也具有優(yōu)勢(shì)，是未來大氣數(shù)據(jù)傳感技術(shù)發(fā)展方向，也是此類飛行器研制的必然選擇。FADS技術(shù)主要依靠飛行器頭部不同位置的測(cè)壓點(diǎn)進(jìn)行壓力測(cè)量，通過壓力場(chǎng)模型和算法解算大氣參數(shù)。對(duì)于特定的飛行器，F(xiàn)ADS 的布局、模型和算法與飛行器復(fù)雜外形、飛行空域、速域等直接相關(guān)[4-5]。一般采用理論、數(shù)值方法結(jié)合風(fēng)洞試驗(yàn)手段獲得FADS 系統(tǒng)的壓力-氣參映射關(guān)系。秦永明[6]、王鵬[7]等研究確定了FADS系統(tǒng)的風(fēng)洞試驗(yàn)校準(zhǔn)參數(shù)，并針對(duì)不同布局形式提出了考慮攻角誤差、側(cè)滑角誤差和形壓誤差的風(fēng)洞試驗(yàn)校準(zhǔn)方法。由于理論外形與真實(shí)安裝的位置偏差、試驗(yàn)環(huán)境與真實(shí)環(huán)境差異以及試驗(yàn)?zāi)Ｐ团c真實(shí)外形尺寸差異等因素，還需借助飛行試驗(yàn)對(duì)解算算法進(jìn)行進(jìn)一步校準(zhǔn)和驗(yàn)證。目前，飛行校準(zhǔn)方法常采用拖錐法、GPS 法等[8-10]。其中，拖錐法只可校準(zhǔn)靜壓，GPS 法則需要在同一航線段上往返飛行以消除空中風(fēng)的影響，但對(duì)于高空高速飛行器飛行剖面，航線一般無法重復(fù)。劉朝君等[11]針對(duì)FADS系統(tǒng)的校準(zhǔn)問題，提出了一種基于探空氣球獲取試驗(yàn)空域大氣靜壓、靜溫、風(fēng)速風(fēng)向的空速校準(zhǔn)方法，但也只可對(duì)總壓、靜壓和馬赫數(shù)進(jìn)行校準(zhǔn)，無法校準(zhǔn)攻角、側(cè)滑角偏差。

本文針對(duì)某飛行器FADS的飛行試驗(yàn)校準(zhǔn)問題，研制了一套多氣參測(cè)量的壓力型基準(zhǔn)大氣數(shù)據(jù)系統(tǒng)。設(shè)計(jì)了可實(shí)現(xiàn)高線性度、解耦的壓力型姿態(tài)角測(cè)量空速管，配套研制了相應(yīng)的大氣數(shù)據(jù)計(jì)算機(jī)系統(tǒng)，并完成了基準(zhǔn)空速管的風(fēng)洞試驗(yàn)標(biāo)定和聯(lián)調(diào)試驗(yàn)校核。該系統(tǒng)具有小型化(對(duì)飛行器氣動(dòng)干擾小)、穩(wěn)定性好和測(cè)量精度高等優(yōu)點(diǎn)。通過將基準(zhǔn)大氣數(shù)據(jù)系統(tǒng)布置在飛行器翼尖，完成總壓、靜壓、馬赫數(shù)、攻角、側(cè)滑角等全氣參測(cè)量，實(shí)現(xiàn)了對(duì)FADS系統(tǒng)的飛行試驗(yàn)校準(zhǔn)和驗(yàn)證。

1 應(yīng)用背景及難點(diǎn)

基于某飛行器FADS系統(tǒng)的飛行試驗(yàn)校準(zhǔn)需求，提出了圖1所示的以壓力型姿態(tài)角測(cè)量空速管和大氣數(shù)據(jù)計(jì)算機(jī)為基準(zhǔn)系統(tǒng)的方法：兩根空速管分別安裝于飛行器翼尖位置，用于總壓、靜壓、攻角和側(cè)滑角的測(cè)量，總溫傳感器(圖1中未標(biāo)出)獲得總溫信息，最終通過大氣數(shù)據(jù)計(jì)算機(jī)解算獲得飛行器基準(zhǔn)大氣參數(shù)，實(shí)現(xiàn)FADS系統(tǒng)的飛行試驗(yàn)校準(zhǔn)和驗(yàn)證。

圖1 基準(zhǔn)大氣數(shù)據(jù)系統(tǒng)應(yīng)用方法示意圖

上述方法的難點(diǎn)和關(guān)鍵具體如下。

① 當(dāng)前，探出式大氣數(shù)據(jù)系統(tǒng)多采用“空速管+風(fēng)標(biāo)角度傳感器”的形式獲得總壓、靜壓、速度、高度、攻角和側(cè)滑角等氣參[12]，用于本文FADS飛行校準(zhǔn)試驗(yàn)時(shí)存在系統(tǒng)復(fù)雜、測(cè)量精度和可靠性差的問題，同時(shí)，會(huì)對(duì)飛行器局部流場(chǎng)產(chǎn)生影響。本文采用的壓力型多氣參測(cè)量的空速管具有體積小、穩(wěn)定性好和測(cè)量精度高等優(yōu)點(diǎn)，其關(guān)鍵難點(diǎn)在于有限空間內(nèi)多氣路壓力的同時(shí)精確測(cè)量，需突破多氣路之間的密封、同時(shí)精確取氣的關(guān)鍵技術(shù)。

② 由于FADS的存在，基準(zhǔn)空速管無法設(shè)置于飛行器頭部，而只能設(shè)置在流場(chǎng)相對(duì)干凈的翼尖兩側(cè)，但機(jī)翼依然會(huì)對(duì)空速管流場(chǎng)產(chǎn)生影響，常規(guī)的氣動(dòng)型面補(bǔ)償[13]和標(biāo)定方法已不適用。本文的另一關(guān)鍵難點(diǎn)在于風(fēng)洞試驗(yàn)標(biāo)定及氣參解算，需通過標(biāo)定和多參數(shù)擬合修正，獲得機(jī)體干擾下的空速管測(cè)量修正解算公式，解決壓力修正、氣參轉(zhuǎn)換問題，實(shí)現(xiàn)基準(zhǔn)氣參的高精度測(cè)量。

2 研制方案

2.1 壓力型姿態(tài)角測(cè)量空速管

2.1.1 測(cè)量原理

壓力型姿態(tài)角測(cè)量空速管如圖2所示，即在常規(guī)空速管的基礎(chǔ)上設(shè)置用于攻角和側(cè)滑角測(cè)量的壓力孔，通過標(biāo)定和解算，實(shí)現(xiàn)飛行器姿態(tài)角測(cè)量。通過設(shè)計(jì)相對(duì)規(guī)則的外形，借助CFD手段對(duì)空速管表面壓力分布進(jìn)行分析，提取對(duì)總壓或靜壓最敏感、對(duì)氣動(dòng)姿態(tài)角最不敏感的區(qū)域，在此布置總、靜壓測(cè)壓點(diǎn)；提取出對(duì)氣動(dòng)姿態(tài)角變化最為敏感、隨位置變化壓力值平穩(wěn)的區(qū)域，布置姿態(tài)角測(cè)壓點(diǎn)。在此基礎(chǔ)上，結(jié)合CFD數(shù)值計(jì)算和真實(shí)模型風(fēng)洞試驗(yàn)標(biāo)定以及聯(lián)調(diào)驗(yàn)證，獲得精度較高的氣參解算模型。通過上述方法，可在較寬速域、姿態(tài)角范圍內(nèi)獲得線性度較高的測(cè)量結(jié)果，同時(shí)，所有參數(shù)均通過壓力測(cè)量和解算實(shí)現(xiàn)。該方法測(cè)量一致性好，系統(tǒng)簡(jiǎn)單，易實(shí)現(xiàn)小型化設(shè)計(jì)。

圖2 基準(zhǔn)空速管測(cè)量原理圖

2.1.2 氣動(dòng)外形設(shè)計(jì)

經(jīng)過優(yōu)化設(shè)計(jì)，基準(zhǔn)空速管的最終氣動(dòng)外形如圖3所示。采用“直桿型”設(shè)計(jì)，測(cè)量部分由橢圓形管鼻、等直段和錐段組成；空速管總長(zhǎng)約350 mm，前段直徑φ13 mm，后段直徑φ20 mm，測(cè)量錐長(zhǎng)20 mm，錐角10°。兩根空速管分別安裝于飛行器翼尖位置，空速管軸線與機(jī)體軸線平行。采用CFD方法(計(jì)算域包含了飛行器，以得到機(jī)翼對(duì)空速管壓力分布干擾下的真實(shí)模擬結(jié)果)獲得其表面的壓力分布信息，進(jìn)而依據(jù)各被測(cè)量的測(cè)量特點(diǎn)，選取壓力孔設(shè)置的敏感域，確定最終靜壓孔、攻角壓力孔和側(cè)滑角壓力孔的布置位置?？账俟芸倝嚎撞贾迷谄渥钋岸耍倝簻y(cè)量基本不受布局影響。

圖3 基準(zhǔn)空速管氣動(dòng)外形設(shè)計(jì)圖

2.1.3 結(jié)構(gòu)設(shè)計(jì)

空速管結(jié)構(gòu)設(shè)計(jì)、各測(cè)壓孔位置如圖4所示，分別包括總壓測(cè)量孔、靜壓測(cè)量孔、上攻角壓力測(cè)量孔、下攻角壓力測(cè)量孔、左側(cè)滑角壓力測(cè)量孔、右側(cè)滑角壓力測(cè)量孔和相對(duì)應(yīng)的引氣管路，同時(shí)，在總壓、靜壓測(cè)量孔附近還設(shè)有加熱系統(tǒng)。除總壓測(cè)量孔外，其他壓力測(cè)量孔均在其敏感域附近均勻地開有多個(gè)小孔。經(jīng)內(nèi)部氣路結(jié)構(gòu)設(shè)計(jì)，實(shí)現(xiàn)了6路壓力的獨(dú)立測(cè)量，且各氣路分別獨(dú)立密封，避免氣路間串氣。

圖4 基準(zhǔn)空速管結(jié)構(gòu)設(shè)計(jì)圖及各測(cè)壓孔位置

2.2 大氣數(shù)據(jù)計(jì)算機(jī)

基準(zhǔn)大氣數(shù)據(jù)計(jì)算機(jī)組成及工作原理如圖5所示，由IPT壓力傳感器、金屬殼體(結(jié)構(gòu)件)、電連接器、電路板組裝件和一些輔助材料組成。電路板上集成了A/D轉(zhuǎn)換器、DSP、通信電路等，可以同時(shí)完成12路壓力信號(hào)和1路溫度傳感器信號(hào)的采集和調(diào)理，并能在DSP內(nèi)完成動(dòng)壓、靜壓、總壓、真空速、表速、氣壓高度、下沉率、總溫、靜溫、密度、大氣機(jī)工作狀態(tài)字、空速管狀態(tài)字、攻角、側(cè)滑角、馬赫數(shù)等大氣參數(shù)的解算與發(fā)送。圖6為最終研制完成的基準(zhǔn)大氣數(shù)據(jù)系統(tǒng)實(shí)物圖。

圖5 大氣數(shù)據(jù)計(jì)算機(jī)設(shè)計(jì)原理圖

圖6 基準(zhǔn)大氣數(shù)據(jù)系統(tǒng)實(shí)物圖

3 風(fēng)洞試驗(yàn)標(biāo)定研究

3.1 風(fēng)洞試驗(yàn)標(biāo)定方法

基準(zhǔn)空速管研制完成后，對(duì)其進(jìn)行了風(fēng)洞試驗(yàn)標(biāo)定，如圖7所示，以獲得空速管的修正公式和固化修正迭代次數(shù)，同時(shí)對(duì)空速管測(cè)量、加熱、密封等性能進(jìn)行檢驗(yàn)。一般空速管風(fēng)洞試驗(yàn)標(biāo)定是將成品實(shí)物置于風(fēng)洞設(shè)備中，給定來流馬赫數(shù)、攻角、側(cè)滑角、總/靜壓等狀態(tài)參數(shù)，測(cè)量空速管對(duì)應(yīng)量的輸出值，與真值比較并進(jìn)行修正使其最大限度逼近真值。由于本文基準(zhǔn)大氣數(shù)據(jù)系統(tǒng)設(shè)置在飛行器兩側(cè)翼尖，壓力測(cè)量受機(jī)翼干擾較大，勢(shì)必與攻角、側(cè)滑角耦合。因此，綜合考慮風(fēng)洞設(shè)備尺寸限制，采用“飛行器半模+空速管成品實(shí)物”作為試驗(yàn)?zāi)Ｐ瓦M(jìn)行標(biāo)定。另外，考慮到基準(zhǔn)大氣數(shù)據(jù)系統(tǒng)的高精度測(cè)量需求，對(duì)風(fēng)洞進(jìn)行了降噪處理。

圖7 風(fēng)洞試驗(yàn)標(biāo)定

3.2 標(biāo)定結(jié)果分析

風(fēng)洞試驗(yàn)完成了20余車次、共計(jì)160個(gè)狀態(tài)點(diǎn)的標(biāo)定工作。典型試驗(yàn)結(jié)果如圖8所示，圖8中CPt、CP∞、ΔCP分別為總壓系數(shù)、靜壓系數(shù)和姿態(tài)角壓差系數(shù)(攻角壓差ΔPα=P下-P上，側(cè)滑角壓差ΔPβ=P右-P左)。圖8(a)、圖8(b)分別為基準(zhǔn)狀態(tài)下空速管總壓系數(shù)、靜壓系數(shù)隨攻角變化的重復(fù)性試驗(yàn)，由圖可知，3次試驗(yàn)重復(fù)性良好，系統(tǒng)測(cè)量性能可靠。圖8(c)、圖8(d)分別為Ma=0.8時(shí)空速管總壓系數(shù)、靜壓系數(shù)隨姿態(tài)角的變化，由圖可知，總壓基本不隨攻角和側(cè)滑角變化，靜壓則主要隨側(cè)滑角呈現(xiàn)出一定的變化規(guī)律。圖8(e)為Ma=0.8時(shí)，姿態(tài)角壓差系數(shù)隨攻角的變化，由圖可知，隨著攻角增加，攻角壓差呈現(xiàn)單調(diào)線性增加，不同側(cè)滑角下變化趨勢(shì)相同，但仍有微小差別；側(cè)滑角壓差呈現(xiàn)減小趨勢(shì)，且隨著β=-5°～5°變化，該變化趨勢(shì)逐漸增大。圖8(f)為Ma=0.8時(shí)，姿態(tài)角壓差系數(shù)隨側(cè)滑角變化，由圖可知，隨側(cè)滑角增加，側(cè)滑角壓差呈現(xiàn)單調(diào)線性增加，不同攻角下變化趨勢(shì)相同，但仍有微小差別；攻角壓差呈現(xiàn)增加趨勢(shì)，且隨著α=-2.7°～12.3°變化，該變化趨勢(shì)逐漸增大。圖8(g)為β=0°時(shí)不同Ma下姿態(tài)角壓差隨攻角變化，由圖可知，不同Ma下攻角壓差的變化規(guī)律呈現(xiàn)出了良好的一致性，但各曲線之間數(shù)值仍有微小差別。圖8(h)為α=0.3°時(shí)，不同Ma下姿態(tài)角壓差隨側(cè)滑角的變化，由圖可知，不同Ma下側(cè)滑角壓差的變化規(guī)律呈現(xiàn)出了良好的一致性，但各曲線之間數(shù)值也呈現(xiàn)出一定差別。

圖8 標(biāo)定數(shù)據(jù)分析

由以上分析可知，基準(zhǔn)空速管的總壓Pt、靜壓P∞、攻角α、側(cè)滑角β的測(cè)量均由各對(duì)應(yīng)測(cè)壓孔的測(cè)量主量所主導(dǎo)，機(jī)翼對(duì)各壓力測(cè)量產(chǎn)生了一定量的微小干擾，這些干擾均在預(yù)期的范圍內(nèi)，可依據(jù)3.3節(jié)中的數(shù)據(jù)處理方法進(jìn)行修正。

3.3 數(shù)據(jù)處理分析

基準(zhǔn)空速管直接測(cè)量獲得總壓、靜壓、上攻角壓力、下攻角壓力、左側(cè)滑角壓力和右側(cè)滑角壓力，共計(jì)6路。標(biāo)定的目的是通過風(fēng)洞試驗(yàn)建立目標(biāo)真值(總壓、靜壓、攻角、側(cè)滑角)與空速管直接測(cè)量的6路壓力之間的函數(shù)關(guān)系，即獲得修正解算公式。空速管直接測(cè)量到的壓力數(shù)據(jù)與真實(shí)值(總壓、靜壓、攻角、側(cè)滑角)比較，應(yīng)考慮這些測(cè)壓孔位受Ma、α、β這3個(gè)參數(shù)影響，需通過回歸分析方法處理實(shí)際測(cè)量值以得到目標(biāo)值的修正解算公式。

為防止回歸過程中出現(xiàn)過擬合現(xiàn)象，同時(shí)進(jìn)一步提高標(biāo)定公式的預(yù)測(cè)精度，采用多元逐步回歸的方法獲得解算公式。對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行分析后，依據(jù)經(jīng)驗(yàn)，找出對(duì)被解算量影響較大的參數(shù)，建立以此參數(shù)為中心的目標(biāo)方程，方程階次擴(kuò)展至包含純3次項(xiàng)和交叉3次項(xiàng)，后以“擬合誤差最小”為原則進(jìn)行逐步回歸，最終確定修正解算方程?；陲L(fēng)洞標(biāo)定數(shù)據(jù)和上述數(shù)據(jù)處理方法，得到的標(biāo)定公式為

ΔCP∞=0.04326-0.08025Ma+0.10995Ma2-0.00005α2-

0.00322β-0.00049β2

(1)

α=2.27216+46.295ΔCPα+5.76761Ma2-6.92605Ma-

12.84476ΔCPαMa

(2)

β=-2.7583+58.32468ΔCPβ+4.07995Ma-3.89947Ma2-

15.60627ΔCPβMa+1.28437ΔCPβα

(3)

其中：總壓由直接測(cè)量獲得，馬赫數(shù)、動(dòng)壓等參數(shù)則直接由計(jì)算公式給出，這里不再羅列。

從以上公式可以看出，雖然考慮了三次權(quán)重項(xiàng)，但最優(yōu)回歸結(jié)果只包含部分二次項(xiàng)，同時(shí)，從系數(shù)值看，各被測(cè)量均由各對(duì)應(yīng)測(cè)壓孔的主量所主導(dǎo)。這說明，一方面空速管的測(cè)量線性度高，達(dá)到了設(shè)計(jì)預(yù)期；另一方面，回歸方法在一定程度上避免了風(fēng)洞試驗(yàn)數(shù)據(jù)源誤差引入的過擬合問題。圖9為擬合數(shù)據(jù)集的誤差帶分布，由圖9可知，各測(cè)量的誤差帶較小，Ma誤差集中在0.005以內(nèi)，總壓、靜壓誤差集中在250 Pa以內(nèi)。

圖9 擬合數(shù)據(jù)集誤差分布

4 性能分析與討論

4.1 聯(lián)調(diào)試驗(yàn)方法

基準(zhǔn)空速管風(fēng)洞試驗(yàn)標(biāo)定獲得修正公式后，用于大氣數(shù)據(jù)計(jì)算機(jī)系統(tǒng)解算軟件的編制。在大氣數(shù)據(jù)計(jì)算機(jī)研制和軟件調(diào)試完成后，與空速管合并組成完整的基準(zhǔn)大氣數(shù)據(jù)系統(tǒng)在風(fēng)洞中進(jìn)行了聯(lián)調(diào)試驗(yàn)，如圖10所示。其目的是在風(fēng)洞中進(jìn)行全系統(tǒng)正常工況下的飛行參數(shù)測(cè)量精度及系統(tǒng)穩(wěn)定性考核。風(fēng)洞只提供流場(chǎng)、給定空速管姿態(tài)，并記錄相應(yīng)風(fēng)洞氣參，空速管壓力測(cè)量及相應(yīng)的氣參解算由大氣數(shù)據(jù)計(jì)算機(jī)完成。

圖10 風(fēng)洞試驗(yàn)聯(lián)調(diào)

4.2 聯(lián)試性能分析

共完成了10車次的聯(lián)調(diào)風(fēng)洞試驗(yàn)。圖11為典型狀態(tài)下，大氣數(shù)據(jù)計(jì)算機(jī)解算獲得的氣參對(duì)比，其中，紅線為本文采用逐步回歸法獲得的氣參值，黑線為采用考慮了所有二次項(xiàng)的全元回歸法獲得的氣參值。風(fēng)洞氣參(真值)為Ma=0.8；β=5°；α=-2.7°，-0.7°，0.3°，1.3°，3.3°，5.3°，7.3°，9.3°，11.3°。由圖11可知，采用本文方法可獲得比全元回歸精度相對(duì)較高的氣參值，尤其是對(duì)于攻角和側(cè)滑角的解算。表1給出了聯(lián)調(diào)試驗(yàn)的標(biāo)準(zhǔn)偏差，為方便比較，將擬合樣本集誤差也列在其中。由表1可知，基準(zhǔn)大氣數(shù)據(jù)系統(tǒng)各飛行參數(shù)測(cè)量及解算精度滿足均較高，達(dá)到了研制目標(biāo)。

圖11 典型狀態(tài)(Ma=0.8，β=5°)下解算后氣參對(duì)比

表1 風(fēng)洞試驗(yàn)標(biāo)準(zhǔn)偏差

FADS飛行校準(zhǔn)試驗(yàn)于近期圓滿完成，獲得了滿意的試驗(yàn)效果。

5 結(jié)束語(yǔ)

針對(duì)某飛行器FADS系統(tǒng)的飛行試驗(yàn)校準(zhǔn)問題，研制了一套壓力型姿態(tài)角測(cè)量的多氣參基準(zhǔn)大氣數(shù)據(jù)系統(tǒng)：

① 設(shè)計(jì)了可實(shí)現(xiàn)高線性度、解耦的壓力型姿態(tài)角測(cè)量空速管，解決了有限空間內(nèi)多氣路之間密封、同時(shí)精確取氣和壓力精確測(cè)量的難題，具有集成度高、體積小、穩(wěn)定性和測(cè)量精度高的優(yōu)點(diǎn)。

② 配套研制了大氣數(shù)據(jù)計(jì)算機(jī)，解決了飛行器機(jī)體干擾下的風(fēng)洞試驗(yàn)標(biāo)定和多參數(shù)擬合修正、氣參高精度解算的難題，獲得了機(jī)體干擾下的空速管測(cè)量修正解算公式，實(shí)現(xiàn)了基準(zhǔn)氣參的高精度測(cè)量。

③ 通過基準(zhǔn)空速管的風(fēng)洞試驗(yàn)標(biāo)定和聯(lián)調(diào)試驗(yàn)，實(shí)現(xiàn)了對(duì)全系統(tǒng)正常工況下的飛行參數(shù)測(cè)量精度及系統(tǒng)穩(wěn)定性考核，最終成功應(yīng)用于飛行試驗(yàn)中。