一種大口徑蝶閥流量特性的坐標(biāo)定位回歸算法

朱美印，張 松，但志宏，裴希同，王 信，王 曦

(1.北京航空航天大學(xué)能源與動(dòng)力工程學(xué)院，北京100191；2.先進(jìn)航空發(fā)動(dòng)機(jī)協(xié)同創(chuàng)新中心，北京100191；3.中國航發(fā)四川燃?xì)鉁u輪研究院，四川綿陽621703)

一種大口徑蝶閥流量特性的坐標(biāo)定位回歸算法

朱美印1，2，張 松3，但志宏3，裴希同3，王 信3，王 曦1，2

(1.北京航空航天大學(xué)能源與動(dòng)力工程學(xué)院，北京100191；2.先進(jìn)航空發(fā)動(dòng)機(jī)協(xié)同創(chuàng)新中心，北京100191；3.中國航發(fā)四川燃?xì)鉁u輪研究院，四川綿陽621703)

大口徑蝶閥是高空模擬試車臺(tái)的一個(gè)重要設(shè)備，其數(shù)學(xué)模型建立依賴于大口徑蝶閥的流量特性。為獲得大口徑蝶閥的流量特性，提出了一種基于特征坐標(biāo)，從大稀疏分散試驗(yàn)數(shù)據(jù)中提取連續(xù)特征信息的大口徑蝶閥流量特性的坐標(biāo)定位回歸算法。即先對無規(guī)律的試驗(yàn)數(shù)據(jù)過濾預(yù)處理，再基于特征坐標(biāo)對預(yù)處理數(shù)據(jù)在特征坐標(biāo)軸上進(jìn)行壓縮、合并。在此基礎(chǔ)上提出了迭代擴(kuò)充回歸流量特性系數(shù)的算法，從而獲得所需的大口徑蝶閥的流量特性。與高空模擬試車臺(tái)的試驗(yàn)數(shù)據(jù)對比表明，利用建立的大口徑蝶閥流量特性所計(jì)算的蝶閥流量與測量流量之間的穩(wěn)態(tài)誤差均在5%以內(nèi)。

高空模擬試車臺(tái)；大口徑蝶閥；流量特性；特征坐標(biāo)定位法；迭代擴(kuò)充算法

1 引言

航空發(fā)動(dòng)機(jī)高空模擬試車臺(tái)(簡稱高空臺(tái))[1-4]的各個(gè)子系統(tǒng)中，排氣管網(wǎng)是高空臺(tái)實(shí)現(xiàn)對試驗(yàn)發(fā)動(dòng)機(jī)飛行高度條件模擬的子系統(tǒng)，而該飛行高度條件的模擬依賴于排氣管網(wǎng)中出口調(diào)節(jié)閥的控制精度[5]。德國斯圖加特高空試車臺(tái)在排氣管網(wǎng)中采用直徑2 m的大口徑蝶閥實(shí)現(xiàn)了對空氣流量為140 kg/s的BR700民用發(fā)動(dòng)機(jī)的飛行環(huán)境模擬[6-7]。美國阿諾德工程發(fā)展中心(AEDC)建有國際一流的高空試車臺(tái)[8]，其用于航空發(fā)動(dòng)機(jī)試驗(yàn)的高空臺(tái)也采用大口徑蝶閥作為其排氣管網(wǎng)系統(tǒng)的出口調(diào)節(jié)閥來實(shí)現(xiàn)對飛行高度條件的模擬[9-10]。

國內(nèi)新建高空臺(tái)的排氣管網(wǎng)與國外的類似，也采用大口徑蝶閥實(shí)現(xiàn)對飛行高度的模擬，但目前尚未獲得該閥的流量特性，且有關(guān)大口徑蝶閥流量特性的參考文獻(xiàn)在公開資料中很少出現(xiàn)[11-15]。因此，如何獲取準(zhǔn)確的大口徑蝶閥流量特性成為一個(gè)技術(shù)瓶頸。針對這一技術(shù)難題，本文提出了一種基于特征坐標(biāo)從大稀疏分散試驗(yàn)數(shù)據(jù)中提取連續(xù)特征信息的大口徑蝶閥流量特性的坐標(biāo)定位回歸算法，獲得了所需的大口徑蝶閥的流量特性，并與高空臺(tái)調(diào)試試驗(yàn)數(shù)據(jù)進(jìn)行了對比驗(yàn)證。

2 大口徑蝶閥數(shù)據(jù)預(yù)處理

由大口徑蝶閥空氣質(zhì)量流量計(jì)算公式[11-13](式(1))可知，影響其質(zhì)量流量的因素有：蝶閥流量系數(shù)Φ，蝶閥流通面積S與管道流通截面積 A1的比值m，空氣密度 ρ，以及閥前壓力 p1。

根據(jù)氣體狀態(tài)方程和m=S/A1，可將式(1)轉(zhuǎn)換為：

式中：T1為閥前溫度，R為氣體常數(shù)?？梢姡狄酝?，參數(shù) p1、T1、S均可由測量計(jì)算得到。為了可采用式(2)計(jì)算流量，將不同壓比和不同開度下對應(yīng)的流量系數(shù)的插值表定義為閥門的流量特性。

流量系數(shù)定義[14-15]如下：

式中：pr為閥后壓力 p2與閥前壓力 p1的比值；k為氣體的絕熱指數(shù)；u為流束收縮系數(shù)，工程應(yīng)用中通常為常數(shù)。

蝶閥流量的調(diào)節(jié)參數(shù)為蝶閥開度VP，Φ用VP可表示為：

式中：r為蝶閥半徑。從上述可看出，流量系數(shù)僅與閥門開度和壓比有關(guān)。因此，準(zhǔn)確計(jì)算大口徑蝶閥的空氣質(zhì)量流量直接依賴于準(zhǔn)確的流量系數(shù)。

表1為新建高空臺(tái)4次調(diào)試試驗(yàn)中獲得的大口徑蝶閥167個(gè)穩(wěn)態(tài)點(diǎn)的數(shù)據(jù)。其中前5個(gè)參數(shù)Q、p1、p2、T1、VP均為試驗(yàn)測量所得，后3個(gè)參數(shù)m、pr、Φ則是通過測量參數(shù)計(jì)算所得。為了從蝶閥試驗(yàn)數(shù)據(jù)中提取不同工作狀態(tài)的信息，需要分段對表中試驗(yàn)數(shù)據(jù)進(jìn)行選擇處理，獲得較為完整的內(nèi)嵌蝶閥特征性能數(shù)據(jù)。

采用的過濾預(yù)處理篩選方法如下：首先根據(jù)蝶閥的流量計(jì)算公式可知，p1在流量計(jì)算中所占的權(quán)重最大，故以其作為排序基準(zhǔn)按照由小到大的順序重新排列；其次這8列數(shù)據(jù)中存在一些狀態(tài)比較接近的數(shù)據(jù)，對其取平均值篩選，分組存放。過濾預(yù)處理后可獲得反映不同穩(wěn)態(tài)工作點(diǎn)特性的114個(gè)去噪聲數(shù)據(jù)。

3 坐標(biāo)定位回歸算法

3.1 坐標(biāo)特性定位回歸的思想

考慮到蝶閥流量系數(shù)與特征參數(shù)的相關(guān)性，可通過建立特征坐標(biāo)的方法對反映流量特性的數(shù)據(jù)進(jìn)行定位。該特征坐標(biāo)由對流量系數(shù)貢獻(xiàn)最大的壓比和閥門開度這兩個(gè)特征參數(shù)確定，其大小反映了試驗(yàn)數(shù)據(jù)的特征分布，并將相近的數(shù)據(jù)點(diǎn)進(jìn)行基于特征坐標(biāo)的壓縮、合并。從而獲得了兩個(gè)特征參數(shù)一個(gè)不變而另一個(gè)變化時(shí)的函數(shù)關(guān)系，以補(bǔ)足后續(xù)采用擬合方法來回歸流量特性的必要條件。對二維坐標(biāo)在蝶閥的工作范圍進(jìn)行掃描遍歷，可構(gòu)造出完整的樣本數(shù)據(jù)空間。最后采用迭代擴(kuò)充回歸流量特性系數(shù)的算法獲得蝶閥的流量特性系數(shù)。

3.2 坐標(biāo)特性定位描述回歸處理算法

(1) 建立特征坐標(biāo)

影響蝶閥流量的特征參數(shù)主要有壓比、閥門開度及流量系數(shù)3個(gè)，而流量系數(shù)與壓比、閥門開度線性相關(guān)，故獨(dú)立變化的特征參數(shù)只有壓比和閥門開度2個(gè)。為此，建立以閥門開度為橫坐標(biāo)、壓比為縱坐標(biāo)的二元特征坐標(biāo)的流量系數(shù)函數(shù)關(guān)系：

上述建立流量系數(shù)函數(shù)關(guān)系的步驟如下：

Step1：以VP為基準(zhǔn)，在excel表格中將數(shù)據(jù)從小到大排列，構(gòu)建excel表格的橫軸數(shù)據(jù)。

Step2：以 pr為基準(zhǔn)，在excel表格中將數(shù)據(jù)從小到大排列，構(gòu)建excel表格的縱軸數(shù)據(jù)。

Step3：將上述橫軸數(shù)據(jù)和縱軸數(shù)據(jù)構(gòu)建流量系數(shù)函數(shù)關(guān)系表，如表2所示(由于數(shù)據(jù)量龐大，表中僅給出部分結(jié)果)。

(2) 壓縮、合并處理

pr的處理：考慮到壓力傳感器測量絕對精度為1 kPa，蝶閥前后壓力傳感器偏差在2 kPa以內(nèi)屬于正常范圍。按此邏輯推理，當(dāng) pr值的偏差在0.01以內(nèi)時(shí)，可壓縮、合并為一個(gè) pr值。

VP的處理：考慮到蝶閥開度的控制精度為0.5°及位移傳感器的測量精度，當(dāng)閥門的兩個(gè)不同開度測量值之差小于0.5°時(shí)，可壓縮、合并為一個(gè)VP值。

算例(針對表2)：首先處理 pr縱坐標(biāo)，表中 pr為0.319 4、0.321 8、0.322 9的 3個(gè)坐標(biāo)點(diǎn)偏差小于0.01，可將它們平均為 pr=0.321；再將原來這3個(gè)點(diǎn)對應(yīng)的流量系數(shù)值合并到 pr=0.321這一行，從而補(bǔ)充了這3個(gè)不同開度下的流量系數(shù)。其次處理VP橫坐標(biāo)，表中 VP為5.318°、5.468°、5.486°的這3個(gè)坐標(biāo)點(diǎn)偏差都小于0.5°，平均為VP=5.4°；再將原來這3個(gè)點(diǎn)對應(yīng)的流量系數(shù)值合并到VP=5.4°這一列，從而補(bǔ)充了這3個(gè)不同壓比下的流量系數(shù)。

通過上述縱、橫坐標(biāo)數(shù)據(jù)壓縮、合并后，可獲得足夠多的特征回歸數(shù)據(jù)。處理后，橫軸上的坐標(biāo)點(diǎn)由原來的114個(gè)點(diǎn)壓縮為41個(gè)，縱軸上的坐標(biāo)點(diǎn)由原來的114個(gè)點(diǎn)壓縮為29個(gè)，再加上 pr=1這一邊界點(diǎn)的特殊情況，同一壓比或同一開度下的數(shù)據(jù)就不止一個(gè)流量系數(shù)，滿足了曲線擬合條件。表3示出了經(jīng)上述處理后的部分結(jié)果。

(3) 回歸擬合

表3數(shù)據(jù)分布不均勻，有的行/列中數(shù)據(jù)點(diǎn)較多，而有的行/列中數(shù)據(jù)點(diǎn)較少。為此，采用擴(kuò)充迭代回歸處理算法以獲得較高的回歸精度。具體步驟如下：

Step 1：對表3中數(shù)據(jù)點(diǎn)占據(jù)最多的行和列進(jìn)行標(biāo)記。

Step 2：對標(biāo)記行進(jìn)行一維數(shù)據(jù)樣條曲線擬合，再對網(wǎng)格點(diǎn)通過擬合出的函數(shù)進(jìn)行插值補(bǔ)充，獲得該行網(wǎng)格點(diǎn)的流量系數(shù)值，并以不同顏色對其進(jìn)行標(biāo)注；采用相同處理方法對標(biāo)記列進(jìn)行處理。將這些行和列上補(bǔ)充的新數(shù)據(jù)點(diǎn)加入到其余未進(jìn)行擬合運(yùn)算的行和列中，從而擴(kuò)充了這些行和列的數(shù)據(jù)點(diǎn)。

表2 流量系數(shù)函數(shù)關(guān)系表Table 2 The functional relationship of flow coefficient

表3 對表2的壓縮、合并處理結(jié)果Table 3 Thecompressed and merged processed results of table 2

Step 3：重新標(biāo)記剩余行和列中數(shù)據(jù)點(diǎn)占據(jù)最多的行和列，重復(fù)Step 2，直至所有的行和列擬合完畢。

對表3采用擴(kuò)充迭代回歸處理算法獲得非標(biāo)準(zhǔn)蝶閥流量特性，見表4(表中僅為部分結(jié)果，其中壓比低于0.32的不在大口徑蝶閥的工作范圍內(nèi)，不需要計(jì)算)。表中，pr=0的數(shù)據(jù)是根據(jù)回歸曲線插值所得的數(shù)據(jù)，無填充色的數(shù)據(jù)為計(jì)算出的流量系數(shù)值，用黃色填充的數(shù)據(jù)為同一壓比下的數(shù)據(jù)擬合值，用藍(lán)色填充的數(shù)據(jù)為同一開度下的數(shù)據(jù)擬合值，用灰色填充的數(shù)據(jù)為邊界情況下的流量系數(shù)值。

(4) 規(guī)格化處理

對表4采用VP在30°之前按步長為2°的增量劃分，在30°以上按5°劃分；pr按0.05的等間隔劃分，以此進(jìn)行規(guī)格化處理。具體做法是，用表4對橫、縱坐標(biāo)劃分后形成的網(wǎng)格點(diǎn)進(jìn)行插值計(jì)算，算出相應(yīng)點(diǎn)上的函數(shù)值，獲得的最終的蝶閥流量特性如表5所示，其三維曲面如圖1所示。

表4 大口徑蝶閥非標(biāo)準(zhǔn)的流量特性Table 4 The nonstandard flow characteristics of a large butterfly valve

表5 大口徑蝶閥流量特性表Table 5 The flow characteristics of a large butterfly valve

圖1 大口徑蝶閥流量特性圖Fig.1 The flow characteristics diagram of a large butterfly valve

4 算法驗(yàn)證

高空臺(tái)試驗(yàn)數(shù)據(jù)分為兩部分，分別在不同試驗(yàn)條件下獲得。其中第一部分?jǐn)?shù)據(jù)用于上述算法進(jìn)行蝶閥流量特性的回歸，第二部分?jǐn)?shù)據(jù)用于下述對比驗(yàn)證。

首先建立如圖2、圖3所示的對比驗(yàn)證平臺(tái)，其中圖3中的模塊與本文流量計(jì)算公式對應(yīng)一致。其次驗(yàn)證蝶閥流量特性的正確性，需對試驗(yàn)數(shù)據(jù)進(jìn)行選擇——該數(shù)據(jù)應(yīng)包含穩(wěn)態(tài)和動(dòng)態(tài)過程數(shù)據(jù)。

圖2 蝶閥流量特性對比驗(yàn)證平臺(tái)的頂層結(jié)構(gòu)Fig.2 The top structure of butterfly valve flow characteristics contrastive verification platform

圖3 對比驗(yàn)證平臺(tái)的內(nèi)嵌蝶閥流量特性模塊Fig.3 The embedded butterfly valve flow characteristics module of contrastive verification platform

對比過程如下：構(gòu)建圖2中的4個(gè)一維插值表，將第二部分的輸入數(shù)據(jù)按對應(yīng)的時(shí)間關(guān)系建立在該表中，當(dāng)起動(dòng)運(yùn)行平臺(tái)仿真模型后，將在大口徑蝶閥流量計(jì)算模型的輸出端口上生成流量仿真數(shù)據(jù)；同時(shí)，構(gòu)建圖2中試驗(yàn)數(shù)據(jù)流量的一維插值表，將第二部分的輸出數(shù)據(jù)也按對應(yīng)的時(shí)間關(guān)系建立在該表中，當(dāng)同步起動(dòng)運(yùn)行平臺(tái)仿真模型后，該模塊生成真實(shí)的流量試驗(yàn)數(shù)據(jù)；在示波器中同時(shí)輸入流量仿真數(shù)據(jù)和流量試驗(yàn)數(shù)據(jù)的曲線進(jìn)行對比，結(jié)果如圖4所示?？梢?，穩(wěn)態(tài)過程中的相對誤差均在4%以內(nèi)，動(dòng)態(tài)過程中的相對誤差均在5%以內(nèi)(不包括試驗(yàn)毛刺干擾情況)。圖5、圖6分別為大口徑蝶閥閥門開度曲線和壓比曲線。可見，在試驗(yàn)過程中，隨著蝶閥開度和壓比的變化，蝶閥的流量也隨之變化。

圖4 大口徑蝶閥流量對比結(jié)果Fig.4 The comparison results of a large butterfly valve flow

圖5 大口徑蝶閥閥門開度曲線Fig.5 The curve of a large butterfly valve opening

圖6 大口徑蝶閥壓比曲線Fig.6 The pressure ratio curve of a large butterfly valve

5 結(jié)論

針對大口徑蝶閥存在數(shù)據(jù)冗余、分散、稀疏和無規(guī)律而導(dǎo)致的無法采用直接擬合獲取流量特性的問題，基于特征坐標(biāo)定位方法從試驗(yàn)數(shù)據(jù)中提取了大口徑碟閥的流量特性，并用不同的試驗(yàn)數(shù)據(jù)對該特性進(jìn)行對比驗(yàn)證。結(jié)果表明，穩(wěn)態(tài)過程中的相對誤差不大于4%，動(dòng)態(tài)過程中的相對誤差不大于5%，滿足大口徑蝶閥工程應(yīng)用要求。

[1]侯敏杰.高空模擬試驗(yàn)技術(shù)[M].北京：航空工業(yè)出版社，2014.

[2]周漢忠，劉大響.即將建成投產(chǎn)的SB101高空試車臺(tái)[J].燃?xì)鉁u輪試驗(yàn)與研究，1990，3(1)：1—10.

[3]侯敏杰.SB101高空模擬試車臺(tái)非等流量試驗(yàn)方法研究[J].燃?xì)鉁u輪試驗(yàn)與研究，1995，8(2)：18—20.

[4]侯敏杰，劉冬根.航空發(fā)動(dòng)機(jī)高空臺(tái)的發(fā)展與展望[J].航空科學(xué)技術(shù)，2012，(3)：1—4.

[5]王惠儒.大型航空發(fā)動(dòng)機(jī)試驗(yàn)及試驗(yàn)設(shè)備研究[J].燃?xì)鉁u輪試驗(yàn)與研究，2008，21(1)：13—17.

[6]盧傳義.國外高空臺(tái)的發(fā)展現(xiàn)狀與趨勢(上)[J].國際航空，1997，(9)：37—38.

[7]盧傳義.國外高空臺(tái)發(fā)展的現(xiàn)狀與趨勢(下)[J].國際航空，1997，(10)：38—39.

[8]黃春峰，劉 麟，劉志友，等.美國阿諾德工程發(fā)展中心(AEDC)發(fā)動(dòng)機(jī)試驗(yàn)設(shè)備的發(fā)展[C]//.中國航空學(xué)會(huì)第七屆動(dòng)力年會(huì)論文摘要集.2010.

[9]Davis M，Montgomery P.A flight simulation vision for aero-propulsion altitude ground test facilities[J].Journal of Engineering for Gas Turbines&Power，2005，127(1)：21—31.

[10]Montgomery P，Burdette R，Klepper J，et al.Evolution of a turbine engine test facility to meet the test needs of fu?ture aircraft systems[R].ASME GT-2002-30605，2002.

[11]孫淮清，王建中.流量測量節(jié)流裝置設(shè)計(jì)手冊[M].北京：化學(xué)工業(yè)出版社，2005.

[12]馬卡洛夫 А Н，舍爾曼 М Я.節(jié)流裝置計(jì)算[M].姚文化，譯.北京：冶金工業(yè)出版社，1957.

[13]楊源泉.閥門設(shè)計(jì)手冊[M].北京：機(jī)械工業(yè)出版社，2000.

[14]楊世忠，邢麗娟.調(diào)節(jié)閥流量特性分析及應(yīng)用選擇[J].閥門，2006，(5)：33—36.

[15]朱美印，裴希同，張 松，等.一種輪盤式特種調(diào)節(jié)閥流量特性的修正算法[J].燃?xì)鉁u輪試驗(yàn)與研究，2016，29(5)：40—45.

A coordinate positioning and regression algorithm for the flow characteristics of a large butterfly valve

ZHU Mei-yin1，2，ZHANG Song3，DAN Zhi-hong3，PEI Xi-tong3，WANG Xin3，WANG Xi1，2
(1.School of Energy and Power Engineering，Beijing University of Aeronautics and Astronautics，Beijing 100191，China；2.Collaborative Innovation Center for Advanced Aero-Engine，Beijing 100191，China；3.AECC Sichuan Gas Turbine Establishment，Mianyang 621703，China)

Large butterfly valve is an important equipment in the new Altitude Test Facilities(ATF).The establishment of the mathematical model of ATF depends on the flow characteristics of the large butterfly valve.In order to acquire the flow characteristics of a large butterfly valve,a coordinate positioning and re?gression algorithm which extracts continuous feature information from the large sparse experiment data based on the feature coordinate was proposed.First of all,the irregular data was preprocessed by filtering.Then based on the feature coordinates,the preprocessed data was compressed and merged on the feature ax?is.Finally,an iterative algorithm was proposed to extend the flow coefficient and then the expected flow characteristics of large butterfly valve were acquired.Compared with the experimental data of ATF,the re?sults show that the steady state error between the flow calculated by the method mentioned above and the measurement flow is less than 5%.

Altitude Test Facilities；large butterfly valve；flow characteristics；feature coordinate positioning method；iterative expansion algorithm

V217+.2

A

1672-2620(2017)04-0039-06

2016-08-26；

2017-07-26

航空發(fā)動(dòng)機(jī)高空模擬技術(shù)重點(diǎn)實(shí)驗(yàn)室基金(201504026)

朱美印(1991-)，男，四川瀘州人，博士研究生，主要從事航空發(fā)動(dòng)機(jī)控制、高空臺(tái)數(shù)字仿真平臺(tái)研究及魯棒控制等。