蒲賽虎 郭毅 雷廷萬
(成都飛機設(shè)計研究所 四川省成都市 610000)
靜壓信號是飛機大氣數(shù)據(jù)系統(tǒng)需要采集的重要信息?;陟o壓信號及其變化,可以解算出氣壓高度、升降速度等大氣參數(shù),再結(jié)合總壓信號、總溫信號等,還可以解算出馬赫數(shù)、真空速等大氣參數(shù)[1]。因此,靜壓測量的準(zhǔn)確性對大氣參數(shù)的準(zhǔn)確性有直接影響。國軍標(biāo)對靜壓誤差的允許范圍有明確規(guī)定[2]。
現(xiàn)代民用及軍用飛機常采用嵌入式靜壓傳感器來測量靜壓[3][4]:在機身適當(dāng)位置嵌入靜壓傳感器,并使傳感器外表面與機身蒙皮外表面共形、齊平,在傳感器外表面中心開測壓孔并通過管路與壓力傳感器連接,從而測量出機身當(dāng)?shù)仂o壓,再通過事先標(biāo)定的機身當(dāng)?shù)仂o壓與來流靜壓的換算關(guān)系,就可以得到來流靜壓。從該測量原理上看,嵌入式靜壓傳感器的外表面應(yīng)當(dāng)與機身蒙皮的外表面完全共形、齊平。但實際上,由于制造誤差、安裝誤差、蒙皮外表面涂覆涂料等因素,可能使嵌入式靜壓傳感器外表面低于蒙皮外表面,即兩者之間存在階差(考慮飛機設(shè)計通常遇到的情況,本文僅討論嵌入式靜壓傳感器外表面“低于”蒙皮外表面的情況)。從伯努利方程[5]分析,該階差將導(dǎo)致嵌入式靜壓傳感器測量到的靜壓大于當(dāng)?shù)仂o壓(即靜壓多指),這一靜壓誤差將使得后續(xù)解算出的氣壓高度等大氣參數(shù)有誤差。如果能得到階差對嵌入式靜壓傳感器測量的具體影響量,一方面可據(jù)此在大氣數(shù)據(jù)系統(tǒng)設(shè)計過程中對階差高度提出控制要求,另一方面也可以在階差存在時,對其導(dǎo)致的靜壓誤差進行修正。然而,基于伯努利方程的分析只能得到定性結(jié)果,且只適用于低速情況,而現(xiàn)代飛機的飛行包線通常包含了亞、跨、超聲速范圍,因此,得到亞、跨、超聲速條件下,階差對嵌入式靜壓傳感器測量的影響的定量結(jié)果,對保證靜壓測量的準(zhǔn)確性有重要意義。但尚未見國內(nèi)外公開文獻對此詳細討論。
為此,本文采用計算流體力學(xué)(Computational Fluid Dynamics, 簡稱CFD)方法對不同高度階差情況下,嵌入式靜壓傳感器在亞、跨、超聲速條件下的流動進行了仿真,并提取出了嵌入式靜壓傳感器在各種條件下的靜壓測量值,從而得到了各種高度的階差對嵌入式靜壓傳感器測量的影響的定量結(jié)果,在此基礎(chǔ)上分析了階差導(dǎo)致的靜壓誤差以及由此帶來的氣壓高度誤差。
在本文中,階差對靜壓測量的影響是用有、無階差時的“嵌入式靜壓傳感器感受到的靜壓”的差值來表征的。差值的絕對值越大,則表示階差對靜壓測量的影響越大,反之,則影響越小。所謂“嵌入式靜壓傳感器感受到的靜壓”,在本文中是CFD 計算得到的嵌入式靜壓傳感器測壓孔所在位置的當(dāng)?shù)仂o壓。
CFD 計算針對平板流動模型進行,整體計算區(qū)域如圖1(1)所示:計算區(qū)域是一個長方體,長方體底面是平板,平板正中間是一個嵌入式靜壓傳感器(圖1(2)示出了嵌入式靜壓傳感器附近的局部放大圖,注意到嵌入式靜壓傳感器外表面取為圓形——大多數(shù)飛機上采用此外形。方便起見,用Φ 表示其直徑),長方體計算區(qū)域的左面是來流邊界,右面是出流邊界,前后和上表面是滑移邊界。
圖1:CFD 計算模型
圖2:嵌入式靜壓傳感器附近的壓力云圖
本文對嵌入式靜壓傳感器不同Φ 值(10mm、20mm、40mm、80mm),不同高度階差(0mm、0.2mm、0.4mm、0.6mm、0.8mm),在不同來流馬赫數(shù)(0.2、0.5、0.8、1.1、1.5、2.0、3.0、4.0)條件下的流動進行了仿真。
需要說明的是,本文在CFD 計算所涉及的網(wǎng)格劃分過程中,在平板及嵌入式靜壓傳感器附近生成了較密的粘性網(wǎng)格以模擬邊界層流動[6]。同時,為了保證各算例的可比性,特別注意盡量保證各算例的網(wǎng)格的一致性。
本文中的壓力都是以無量綱化的壓力系數(shù)形式給出的,計算式為:
階差為0mm 時,氣流平順流過嵌入式靜壓傳感器附近區(qū)域,而階差大于0mm 時,相當(dāng)于在嵌入式靜壓傳感器周圍存在一圈臺階,這必然使得氣流不再平順。圖2 示出了Φ 值為20mm、階差為0.8mm時,典型馬赫數(shù)條件下,計算區(qū)域x向中間剖面的壓力云圖(流向是從圖左到圖右)??梢婋A差的確使得嵌入式靜壓傳感器附近流場的壓力出現(xiàn)了明顯的梯度:左側(cè)相當(dāng)于一個后臺階流動,在臺階后出現(xiàn)了一個低壓區(qū)——壓力低于主流壓力(也就低于無階差時壓力);右側(cè)相當(dāng)于一個前臺階流動,在臺階前出現(xiàn)了一個高壓區(qū)——壓力高于主流壓力(也就高于無階差時壓力)。隨著馬赫數(shù)的增加,低壓區(qū)范圍逐漸右移而擴大,而高壓區(qū)范圍則逐漸減小。在亞、跨聲速范圍(圖2(1)),測壓孔處在高壓區(qū),而在超聲速范圍(圖2(2)),測壓孔處在低壓區(qū)。
為了觀察階差對靜壓測量影響的具體數(shù)值(即靜壓誤差),圖3 示出了有、無階差時壓力系數(shù)的差值隨馬赫數(shù)、階差高度、嵌入式靜壓傳感器Φ 值的變化曲線。整體看來,階差將導(dǎo)致百分之幾的壓力系數(shù)誤差,這說明其影響還是可觀的,不能忽略。另外,曲線整體上還呈現(xiàn)如下現(xiàn)象:
(1)Cp 差值和馬赫數(shù)的關(guān)系:Cp 差值在亞聲速時是正值,且隨著馬赫數(shù)的增大而增大,在跨聲速范圍達到最大值,之后隨著馬赫數(shù)的增大而減小直至變成負值,在M>2.0 后,該差值隨馬赫數(shù)變化較小;
(2)Cp 差值和階差高度的關(guān)系:在亞、跨聲速范圍,Cp 差值隨階差高度的增大而增大,在超聲速范圍,該差值隨階差高度變化較?。?/p>
(3)Cp 差值和嵌入式靜壓傳感器Φ 值的關(guān)系:Cp 差值隨著Φ 值的增大而減小。
如引言中所述,靜壓誤差將導(dǎo)致與靜壓相關(guān)的大氣參數(shù)出現(xiàn)誤差。以氣壓高度(與靜壓的關(guān)系式可詳見文獻[7])為例,圖4 示出了在氣壓高度為5000m 時,圖3 所示靜壓誤差導(dǎo)致的氣壓高度誤差(即有階差時解算出的氣壓高度與5000m 之差)。整體看來,階差導(dǎo)致的氣壓高度誤差最大可達百米級,這說明其對氣壓高度測量精度影響是比較大的,不能忽略。同時,與靜壓誤差的現(xiàn)象相對應(yīng),氣壓高度誤差整體還呈現(xiàn)如下現(xiàn)象:
(1)氣壓高度誤差和馬赫數(shù)的關(guān)系:在亞、跨聲速范圍,階差將導(dǎo)致氣壓高度少指,少指量隨著馬赫數(shù)的增大而增大,在跨聲速范圍達到最大值,之后隨著馬赫數(shù)的增大而減小直至變成氣壓高度多指,多指量隨著馬赫數(shù)增大而增大;
(2)氣壓高度誤差和階差高度的關(guān)系:在整個馬赫數(shù)范圍,氣壓誤差隨階差高度的增大而增大;
(3)氣壓高度誤差和和嵌入式靜壓傳感器Φ 值的關(guān)系;隨著Φ 值的增大,氣壓高度誤差逐漸減小。
圖3:有無階差的壓力系數(shù)之差
圖4:階差造成的氣壓高度誤差
本文通過CFD 計算,觀察到階差造成的嵌入式靜壓傳感器附近流場的改變情況,并得到了其對靜壓測量影響的定量結(jié)果。從結(jié)果來看,階差將導(dǎo)致壓力系數(shù)最大有百分之幾的誤差,由此導(dǎo)致氣壓高度最大有百米級的誤差。該分析結(jié)果已在某型號試飛過程中得到驗證??梢婋A差導(dǎo)致的靜壓誤差是較大的,在設(shè)計過程中應(yīng)考慮盡量避免在嵌入式靜壓傳感器周圍出現(xiàn)階差,或者設(shè)計嵌入式靜壓傳感器時應(yīng)在飛機結(jié)構(gòu)允許的條件下適當(dāng)增大其外表面直徑(Φ值),以降低加工、安裝、涂層等因素的影響。若階差已經(jīng)產(chǎn)生,則需考慮相應(yīng)的靜壓誤差修正。本文提供的定量結(jié)果可為這種靜壓誤差修正提供參考。