亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于模糊指數(shù)融合和正交基分解的發(fā)動機性能監(jiān)控

        2013-07-01 23:43:00李冬李本威孫濤宋漢強趙勇
        燃氣渦輪試驗與研究 2013年5期
        關(guān)鍵詞:滑動規(guī)則發(fā)動機

        李冬,李本威,孫濤,宋漢強,趙勇

        基于模糊指數(shù)融合和正交基分解的發(fā)動機性能監(jiān)控

        李冬1,李本威2,孫濤2,宋漢強1,趙勇1

        (1.海軍航空工程學(xué)院研究生管理大隊,山東煙臺264001;2.海軍航空工程學(xué)院飛行器工程系,山東煙臺264001)

        針對發(fā)動機性能監(jiān)控過程中出現(xiàn)的單參數(shù)監(jiān)控信息量不足和多參數(shù)容易矛盾的情況,提出一種發(fā)動機性能監(jiān)控模糊指數(shù)融合方法。該方法基于滑動窗口參數(shù)均值和熵,采取有效的信息融合技術(shù),建立表征發(fā)動機性能的模糊融合指數(shù)。利用神經(jīng)網(wǎng)絡(luò)方法,依據(jù)已確定的模糊融合規(guī)則推導(dǎo)出剩余決策規(guī)則。實例表明,模糊融合指數(shù)能很好地跟蹤發(fā)動機性能緩慢衰退和突變兩種情況。用正交基分解的方法對模糊指數(shù)進行重構(gòu)預(yù)測,與線性擬合和二次擬合相結(jié)合的方法相比,其預(yù)測精度更高,能準確預(yù)測發(fā)動機的性能變化。

        航空發(fā)動機;性能監(jiān)控;模糊信息融合;正交基分解;決策規(guī)則;信息熵

        orthogonal basis decomposition;decision-making logic;information entropy

        1 引言

        航空發(fā)動機長時間工作后,各部件老化將導(dǎo)致其做功能力下降,因此有必要開展發(fā)動機性能監(jiān)控研究,對性能下降的發(fā)動機及時跟蹤,為發(fā)動機視情維修提供指導(dǎo)[1,2]。文獻[3]基于神經(jīng)網(wǎng)絡(luò)優(yōu)化算法對發(fā)動機燃氣溫度裕度(EGTM)進行預(yù)測,并以EGTM的變化表征發(fā)動機性能變化,但EGTM只作為發(fā)動機性能下降的一個重要方面。單參數(shù)評估由于忽略了其它參數(shù)的影響,或參數(shù)間存在矛盾時,容易出現(xiàn)誤判。相比之下,多參數(shù)綜合評估法則相對客觀、合理,可減少誤判機率。文獻[4]~[6]采用綜合加權(quán)法,融合發(fā)動機多個參數(shù)得到性能綜合指數(shù),以此監(jiān)控發(fā)動機性能狀況。文獻[7]和[8]提出基于支持向量機對發(fā)動機性能狀況進行預(yù)測,但其參數(shù)選取過于依賴經(jīng)驗,不易直接獲得。文獻[9]采用回歸樹和模糊神經(jīng)網(wǎng)絡(luò)結(jié)合的方法對性能指數(shù)多步預(yù)測,取得了不錯效果;文獻[10]以最終誤差預(yù)報準則優(yōu)選嵌入維數(shù),通過回歸支持向量機預(yù)測發(fā)動機性能綜合指數(shù),并以此監(jiān)控發(fā)動機性能。但兩種方法同樣存在回歸參數(shù)難以確定的問題。文獻[11]利用組合模型預(yù)測發(fā)動機性能指數(shù),以免疫粒子群算法確定組合模型權(quán)重,較好地預(yù)測了短期的發(fā)動機性能指數(shù)變化。

        本文在上述分析的基礎(chǔ)上,通過構(gòu)建表征發(fā)動機性能變化的參數(shù)數(shù)組和模糊隸屬度函數(shù),融合多參數(shù)信息對發(fā)動機狀況進行評估。利用神經(jīng)網(wǎng)絡(luò)方法,依據(jù)部分已確定的規(guī)則推導(dǎo)出剩余決策規(guī)則,并將所有決策規(guī)則并入到模糊指數(shù)合成中。采用多參數(shù)模糊指數(shù)融合的方法,對發(fā)動機性能緩慢衰退及突變兩種狀況進行監(jiān)控,判斷發(fā)動機所處狀況。采用正交基分解的方法對模糊指數(shù)進行重構(gòu)預(yù)測,并與線性擬合和二次擬合相結(jié)合的方法進行對比。

        2 構(gòu)造性能參數(shù)監(jiān)控的統(tǒng)計量

        在實際發(fā)動機性能監(jiān)控中,反映發(fā)動機運行狀態(tài)的參數(shù)很多,但實際條件的復(fù)雜性,使得這些熱力參數(shù)與發(fā)動機性能狀態(tài)的對應(yīng)關(guān)系具有一定的不確定性和模糊性。由于發(fā)動機性能衰退伴隨使用時間的增加而逐漸增加,所以選擇與飛行時間相關(guān)的性能參數(shù)作為表征發(fā)動機性能狀況的參數(shù)。發(fā)動機監(jiān)控參數(shù)的選擇受傳感器可安裝數(shù)量及位置限制,本文確定低壓轉(zhuǎn)子轉(zhuǎn)速n1、高壓轉(zhuǎn)子轉(zhuǎn)速n2、渦輪后燃氣溫度T4作為發(fā)動機性能監(jiān)控參數(shù)[11]。

        為監(jiān)控發(fā)動機性能狀態(tài),將各類參數(shù)按照滑動窗口進行處理,求取每一滑動窗口內(nèi)數(shù)據(jù)的均值和熵,以此構(gòu)造樣本參數(shù)的統(tǒng)計量,即:

        式中:k=n-r+1,n為參數(shù)時間序列長度,r為滑動窗口長度,mean()為取均值,entropy()為取熵運算,xik為i種參數(shù)(i=3,分別為n1、n2和T4)k個樣本點的數(shù)值,sik為i種參數(shù)k個滑動窗口內(nèi)參數(shù)的平均值,tik為i種參數(shù)k個滑動窗口內(nèi)參數(shù)的熵。平均值反映發(fā)動機性能狀況;熵表示滑動窗口內(nèi)數(shù)據(jù)波動情況,對應(yīng)性能變化程度。熵的計算式為:

        式中:nikj為i種參數(shù)在第k個滑動窗口內(nèi)j組內(nèi)對應(yīng)的個數(shù),hik為i種參數(shù)在k個滑動窗口的信息熵,m為k個滑動窗口內(nèi)參數(shù)的分組數(shù),參數(shù)按等分分組。

        當發(fā)動機性能緩慢衰退時,滑動窗口內(nèi)樣本均值和熵緩慢變化;當發(fā)動機性能突變時,樣本均值和熵相應(yīng)會出現(xiàn)突變。當參數(shù)非均勻變化時,對應(yīng)的熵減小。直觀看,這兩個統(tǒng)計量能反映發(fā)動機的性能狀況。

        隨著使用時間的增加,發(fā)動機會出現(xiàn)性能下降。一般來講,n1和n2隨飛行時間的增加而逐漸偏離正常狀態(tài),T4則逐漸增加,數(shù)據(jù)變化趨勢不同。如果對上述三個參數(shù)進行融合,會使融合結(jié)果出現(xiàn)矛盾,因此有必要對數(shù)據(jù)進行規(guī)范化處理。引入下列變換形式:

        (1)效益型數(shù)據(jù)

        式中:aik可為sik、tik,ai,mink為第i種參數(shù)的最小值,ai,maxk為第i種參數(shù)的最大值,a′ik為處理后的值。定性看,n1、n2平均值屬成本型數(shù)據(jù),T4平均值屬效益型數(shù)據(jù),三種參數(shù)的熵屬成本型數(shù)據(jù)。

        由上述分析可知,構(gòu)造某一樣本點表征發(fā)動機性能的參數(shù)數(shù)組形式為[s1k,t1k,s2k,t2k,s3k,t3k],s1k、t1k、s2k、t2k、s3k、t3k分別表示在第k個樣本點的平均值和熵。

        3 多參數(shù)模糊指數(shù)融合

        采用多參數(shù)模糊指數(shù)融合方法,監(jiān)控發(fā)動機性能狀況。通過構(gòu)建輸入和輸出隸屬度函數(shù),將輸入?yún)?shù)轉(zhuǎn)化為反映發(fā)動機性能的指標。每個輸入對應(yīng)四種狀態(tài),輸出對應(yīng)三種狀態(tài)。構(gòu)造的隸屬度輸入、輸出函數(shù)[12]如圖1所示。

        分析圖1可知,輸入模糊集定義為:ZV:零值,SV:小值,MV:適中值,BV:大值。輸出模糊集定義為:NR:正常范圍,MR:中間狀態(tài),AR:不正常狀態(tài)。其中,輸入?yún)?shù)0~0.33表示ZV~SV,0.33~0.66表示SV~MV,0.66~1.00表示MV~BV。輸出參數(shù)0~0.50表示NR~MR,0.50~1.00表示MR~AR。

        圖1 輸入、輸出參數(shù)的隸屬度函數(shù)Fig.1 Membership functions for input and output parameters

        定義輸入和輸出隸屬度函數(shù)后,結(jié)合模糊決策規(guī)則(依據(jù)學(xué)習(xí)和領(lǐng)域?qū)<医?jīng)驗得到),才能將輸入?yún)?shù)融合,得到模糊融合指數(shù)。研究中,所得診斷規(guī)則基本與專家經(jīng)驗相一致[13]。本文總結(jié)了相關(guān)決策規(guī)則,這些規(guī)則以IF(輸入?yún)?shù)的表示范圍)、THEN(發(fā)動機所處的性能狀況)表示。輸入模糊集中每個參數(shù)對應(yīng)四種不同的狀態(tài)。表1列出了部分決策規(guī)則。

        表1 部分決策規(guī)則Table 1 Part of decision-making logic

        以表1第1條決策規(guī)則為例進行分析可知,當s1、s2、s3分別處于ZV、ZV和SV狀態(tài)時,發(fā)動機性能處于正常狀態(tài)。這些決策規(guī)則構(gòu)成判斷發(fā)動機性能的基礎(chǔ)。模糊輸入集有6個參數(shù),這樣一共會產(chǎn)生46個決策規(guī)則。通過以往經(jīng)驗及理論分析,能得到大部分決策規(guī)則,只有一部分規(guī)則未知??紤]到神經(jīng)網(wǎng)絡(luò)方法處理非線性問題及泛化的能力[14],即使不是樣本集中的輸入也能給出合適的輸出。以確定的決策規(guī)則作為訓(xùn)練樣本,其中輸入標記集為6參數(shù)(形式同性能變化參數(shù)數(shù)組),輸出標記集為3參數(shù)。此時,神經(jīng)網(wǎng)絡(luò)的輸入和輸出標記集中的參數(shù)為標記符號。本文約定:輸入標記集中的ZV:0,SV:1, MV:2,BV:3;輸出標記集中的NR:0,MR:0.5,AR:1。注意到圖1(b),輸出標記集的參數(shù)對應(yīng)兩個隸屬度,即輸出參數(shù)具有模糊性。但結(jié)合確定的模糊規(guī)則,只存在唯一的隸屬度。利用下述模糊融合算法求解模糊指數(shù)。

        由圖2可知,模糊輸入集中每個參數(shù)對應(yīng)2個數(shù)值,6個參數(shù)對應(yīng)12個點。這12個點至少會產(chǎn)生26=64個規(guī)則(如果存在模糊規(guī)則)。模糊指數(shù)融合及預(yù)測流程如圖3所示。

        圖2 輸入?yún)?shù)的模糊集合Fig.2 Fuzzy set of input parameters

        圖3 模糊指數(shù)融合和預(yù)測流程Fig.3 Flow chart of fuzzy index fusion and prediction

        根據(jù)圖3中的計算流程,模糊指數(shù)融合及其預(yù)測算法流程總結(jié)為:

        Step1:原始數(shù)據(jù)經(jīng)滑動窗口和規(guī)范化處理,確定參數(shù)數(shù)組[s1、t1、s2、t2、s3、t3]。

        Step2:模糊指數(shù)集合有64種組合,可產(chǎn)生64個決策規(guī)則。對于某一確定的規(guī)則,依據(jù)模糊合成運算,得到模糊輸出參數(shù)位置。其中,某一規(guī)則對應(yīng)輸出參數(shù)的模糊合成計算式為:

        式中:a1∈{h1,h2},a2∈{h5,h6},a3∈{h9,h10},a4∈{h3,h4},a5∈{h7,h8},a6∈{h11,h12},t為模糊輸出參數(shù)。

        Step3:對于某一輸出模糊集,根據(jù)決策規(guī)則都會產(chǎn)生相應(yīng)的隸屬度值,再根據(jù)式(6)得到模糊融合指數(shù)。

        式中:index為模糊融合指數(shù),b為輸出參數(shù)對應(yīng)的隸屬度。

        Step4:在得到的模糊融合指數(shù)基礎(chǔ)上,利用正交基分解[15]的方法對指數(shù)進行重構(gòu)預(yù)測。取部分數(shù)據(jù)進行訓(xùn)練,建立預(yù)測模型,并利用剩余數(shù)據(jù)驗證所建模型。

        由式(4)、式(5)可看出,模糊融合指數(shù)在0~1之間變化。接近0代表發(fā)動機性能處于正常狀態(tài),接近1則代表性能接近故障狀態(tài)。

        4 實例分析

        a、b兩臺發(fā)動機的性能參數(shù)變化過程如圖4所示,圖中縱坐標經(jīng)濾波后除以最大值進行歸一化處理。實際中,發(fā)動機a性能逐漸衰退,無突變情況。而發(fā)動機b約在第45個樣本點性能出現(xiàn)較嚴重問題,后經(jīng)一段時間排除,性能得到恢復(fù)。

        分析可知,圖4(a)中參數(shù)隨樣本點緩慢增加,表征發(fā)動機性能緩慢衰退;圖4(b)中低壓轉(zhuǎn)子轉(zhuǎn)速約在第45個樣本點附近出現(xiàn)突變,但又很快恢復(fù)到平穩(wěn)狀態(tài),其它2個參數(shù)變化較平穩(wěn),未出現(xiàn)突變。這樣就出現(xiàn)了相互矛盾的情況,各參數(shù)之間看似相悖,不易對發(fā)動機性能做出準確判斷。

        圖5從熵的角度分析發(fā)動機性能狀況。其中,發(fā)動機a 3種參數(shù)的熵基本保持恒定,表征性能變化平穩(wěn);而發(fā)動機b的低壓轉(zhuǎn)子轉(zhuǎn)速熵在第45個樣本點附近驟降,其它2種熵基本不變。3種參數(shù)的熵同樣出現(xiàn)了矛盾,不能正確判斷發(fā)動機性能狀況。

        以模糊融合指數(shù)分析發(fā)動機性能狀況。首先求取滑動窗口內(nèi)三種參數(shù)的均值和熵。經(jīng)規(guī)范化處理后,作為模糊指數(shù)合成隸屬度函數(shù)的輸入。其中,n1、n2的平均值采用式(4)規(guī)范化處理,T4的平均值采用式(3)規(guī)范化處理;3種參數(shù)的熵均用式(4)規(guī)范化處理。按照模糊指數(shù)融合算法,得到模糊指數(shù)??紤]發(fā)動機性能緩慢衰退和運行過程中性能突變兩種情況,跟蹤對應(yīng)模糊融合指數(shù)的變化,結(jié)果見圖6。

        圖4 發(fā)動機原始參數(shù)變化過程Fig.4 Original parameter variation of engine

        由圖6(a)可知,根據(jù)發(fā)動機性能模糊融合指數(shù)判斷,發(fā)動機在1~40時間點性能衰退速率較快,近似成直線趨勢;在45~75時間點內(nèi)變化速率平穩(wěn),指數(shù)維持在0.6附近,表征發(fā)動機性能在這一段時間內(nèi)基本穩(wěn)定。實際中,發(fā)動機在初始階段呈加速衰退趨勢,但部件經(jīng)一段時間磨合后,性能變化較穩(wěn)定。圖中曲線走勢與實際性能參數(shù)緩慢衰退過程相符。圖6(b)的模糊融合指數(shù)表示,發(fā)動機性能約在第45個樣本點附近出現(xiàn)突變,且指數(shù)很大,接近0.8,預(yù)示出現(xiàn)較嚴重問題,與實際情況相符。經(jīng)維修后,相對以前性能得到一定改善,因此第45個樣本點后狀態(tài)較平穩(wěn)。

        圖5 發(fā)動機參數(shù)的熵Fig.5 Stable entropy of engine parameters

        圖6 發(fā)動機性能變化情況Fig.6 Changes of engine performance

        采用正交分解的方法對發(fā)動機性能狀態(tài)進行預(yù)測,同時引入文獻[16]的線性擬合與二次擬合相結(jié)合的預(yù)測方法作比較分析。文獻[16]預(yù)測問題的關(guān)鍵是選取線性預(yù)測和二次預(yù)測轉(zhuǎn)折點,借助偏度統(tǒng)計量,當數(shù)據(jù)偏度出現(xiàn)連續(xù)增大或減小時,應(yīng)改變擬合模型(擬合系數(shù))。以下面兩種判別標準判斷預(yù)測精度。

        (1)誤差的標準差(rsm)

        (2)相對誤差的最大值(rmax)

        式中:n為預(yù)測數(shù)據(jù)個數(shù),yi為模糊融合指數(shù)實際值,為預(yù)測值。

        本文建立100組統(tǒng)計樣本,前80組進行訓(xùn)練,建立預(yù)測模型,后20組數(shù)據(jù)驗證模型。由于模糊融合指數(shù){yi}包含近似線性趨勢項,通過一階差分法將趨勢項直接剔除,從而得到平穩(wěn)時序{?yi},其表達式為?xi=xi-xi-1。對{?xi}建立預(yù)測模型,再根據(jù)差分關(guān)系還原得到{xi}的預(yù)測模型[17]。兩種預(yù)測結(jié)果如圖7、表2所示。

        由圖表中可知,正交基分解方法在誤差的標準差和在相對誤差最大值兩項指標上都優(yōu)于文獻[16]中方法。據(jù)此可認為,正交基分解方法能對模糊融合指數(shù)進行準確預(yù)測,從而更好地監(jiān)控發(fā)動機性能。

        圖7 模糊融合指數(shù)預(yù)測結(jié)果Fig.7 Prediction of fuzzy fusion index

        表2 預(yù)測結(jié)果對比Table 2 Comparison of prediction results

        5 結(jié)論

        本文基于收集的發(fā)動機測量數(shù)據(jù),綜合各測量數(shù)據(jù)的信息,將其融合成反映發(fā)動機性能變化的指數(shù),并基于該指數(shù)監(jiān)控發(fā)動機性能。該方法能克服單參數(shù)監(jiān)控信息量不足及多參數(shù)監(jiān)控容易矛盾的缺點,對發(fā)動機實施有效監(jiān)控。主要研究結(jié)論為:

        (1)基于模糊融合算法,構(gòu)建了輸入和輸出參數(shù)的隸屬度函數(shù),并以此為基礎(chǔ)進一步進行模糊融合,得到當指數(shù)接近0表征發(fā)動機性能正常,接近1表征發(fā)動機接近故障狀態(tài)。模糊融合指數(shù)能很好地監(jiān)控發(fā)動機性能緩慢衰退和突變兩種情況。

        (2)采用正交基分解方法對模糊融合指數(shù)進行重構(gòu)預(yù)測,得到的融合指數(shù)能很好地判斷發(fā)動機的性能狀況,且該預(yù)測方法優(yōu)于文獻[16]中的線性擬合和二次擬合相結(jié)合的方法。

        [1]李冬,浦鵬,譚巍,等.高壓壓氣機性能老化預(yù)測和影響分析[J].燃氣渦輪試驗與研究,2011,24(4):1—5.

        [2]李冬,黃帥,曹明川,等.基于支持向量機的發(fā)動機性能衰退指標分類和預(yù)測[J].燃氣渦輪試驗與研究,2012,25(3):20—25.

        [3]張瑩,左洪福,任淑紅,等.一種航空發(fā)動機性能衰退預(yù)測的優(yōu)化算法[J].中國機械工程,2011,22(9):1009—1013.

        [4]謝壽生.某型飛機發(fā)動機故障診斷與性能趨勢監(jiān)測研究[D].西安:西北工業(yè)大學(xué),1998.

        [5]胡金海,謝壽生.基于遺傳算法的發(fā)動機性能監(jiān)控與故障診斷[J].推進技術(shù),2003,24(3):198—200.

        [6]胡金海,謝壽生,胡劍鋒,等.基于粗糙集理論的航空發(fā)動機性能綜合評判[J].系統(tǒng)工程與電子技術(shù),2006,28 (5):704—707.

        [7]尉詢楷,李應(yīng)紅,王碩,等.基于支持向量機的航空發(fā)動機滑油監(jiān)控分析[J].航空動力學(xué)報,2004,19(3):392—397.

        [8]徐啟華,師軍.基于支持向量機的航空發(fā)動機故障診斷[J].航空動力學(xué)報,2005,20(2):298—302.

        [9]Tran V T,Yang B S,Tan A C C.Multi-Step Ahead Direct Prediction for the Machine Condition Prognosis Using Re?gression Trees and Neuro-Fuzzy Systems[J].Expert Sys?tems with Applications,2009,36(5):9378—9387.

        [10]胡金海,謝壽生,駱廣琦,等.基于支持向量機方法的發(fā)動機性能趨勢預(yù)測[J].推進技術(shù),2005,26(3):260—264.

        [11]李冬,馬力,樊照遠,等.基于免疫粒子群優(yōu)化的不等時間間隔發(fā)動機性能綜合指數(shù)組合預(yù)測[J].燃氣渦輪試驗與研究,2013,26(2):42—45.

        [12]Boutros T,Liang M.Mechanical Fault Detection Using Fuzzy Index Fusion[J].Machine Tools and Manufacture,2007,47(11):1702—1714.

        [13]陳恬,孫健國.粗糙集與神經(jīng)網(wǎng)絡(luò)在航空發(fā)動機氣路故障診斷中的應(yīng)用[J].航空動力學(xué)報,2006,21(1):207—212.

        [14]徐士良.數(shù)值分析與算法[M].北京:機械工業(yè)出版社,2007.

        [15]侯鳳陽.基于人工神經(jīng)網(wǎng)絡(luò)的燃氣輪機氣路故障診斷研究[D].南京:南京航空航天大學(xué),2007.

        [16]Li Y G,Nilkitsaranont P.Gas Turbine Performance Prog?nostic for Condition-Based Maintenance[J].Applied Ener?gy,2009,86(11):2152—2161.

        [17]楊叔子,吳雅.時間序列分析的工程應(yīng)用[M].武漢:華中理工大學(xué)出版社,1996.

        Engine Performance Monitoring Based on Fuzzy Index Fusion and Orthogonal Basis Decomposition

        LI Dong1,LI Ben-wei2,SUN Tao2,SONG Han-qiang1,ZHAO Yong1
        (1.Graduate Students Brigade,Naval Aeronautical Engineering University,Yantai 264001;2.Department of Airborne Vehicle Engineering,Naval Aeronautical Engineering University,Yantai 264001)

        Aiming at information shortage of single parameter and easy contradiction of multi-parameter in the process of engine performance monitor,a method of fuzzy index fusion in the engine performance moni?tor was presented.With effective information fusion technology,based on parameter mean and entropy of sliding window,performance fuzzy fusion index was established.Residual decision-making logic was de?duced by neural network according to pre-determined fuzzy fusion rule.Example indicated fuzzy index could track the engine performance slow deterioration and sudden changes.Fuzzy fusion index was recon?structed and predicted by adopting orthogonal basis decomposition method.Compared with combination of linear and quadratic prediction method,orthogonal basis decomposition method gets higher accuracy,and predicts performance change accurately.

        aero-engine;performance monitoring;fuzzy information fusion;

        V235.13

        A

        1672-2620(2013)05-0035-06

        2012-10-17;

        2013-10-11

        李冬(1984-),男,遼寧葫蘆島人,博士研究生,主要從事航空發(fā)動機狀態(tài)監(jiān)控、性能評估、性能預(yù)測技術(shù)研究。

        猜你喜歡
        滑動規(guī)則發(fā)動機
        撐竿跳規(guī)則的制定
        數(shù)獨的規(guī)則和演變
        發(fā)動機空中起動包線擴展試飛組織與實施
        一種新型滑動叉拉花鍵夾具
        Big Little lies: No One Is Perfect
        讓規(guī)則不規(guī)則
        Coco薇(2017年11期)2018-01-03 20:59:57
        TPP反腐敗規(guī)則對我國的啟示
        新一代MTU2000發(fā)動機系列
        滑動供電系統(tǒng)在城市軌道交通中的應(yīng)用
        一種基于變換域的滑動聚束SAR調(diào)頻率估計方法
        熟女丝袜美腿亚洲一区二区三区 | 国产av91在线播放| 亚洲精品成人久久av| 久久青青草原一区网站| 久久一区二区三区少妇人妻| 亚洲理论电影在线观看| 日韩精品无码中文字幕电影| 热の国产AV| 国产精品九九热| 久久国产A∨一二三| 亚洲一区二区三区乱码在线| 色婷婷亚洲一区二区三区在线| 国产香港明星裸体xxxx视频| 一本久久a久久精品vr综合| 思思久久96热在精品国产| av天堂精品久久久久| 国产av午夜精品一区二区入口 | 麻豆人妻无码性色AV专区 | 久久综合伊人77777麻豆| 天天爽夜夜爽人人爽| 日日碰狠狠丁香久燥| 亚洲色大成网站www在线观看| 校花高潮一区日韩| 国产视频一区二区三区观看| 国产一区二区视频免费在| 成人网站在线进入爽爽爽| 男女18禁啪啪无遮挡| 99国产精品丝袜久久久久| 日本一区中文字幕在线播放| 国产一区二区三区在线视频观看| 人人妻人人做人人爽| 亚洲av中文无码字幕色三 | 色两性网欧美| 国产大片中文字幕| 美女福利视频网址导航| 在厨房拨开内裤进入毛片| 久久精品国产视频在热| 国产成人无码区免费网站| 久久久久久久尹人综合网亚洲| 日本高清在线播放一区二区| 包皮上有一点一点白色的|