• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    ELM在航空發(fā)動機(jī)氣路部件故障診斷的應(yīng)用研究

    2016-08-10 08:02:09王小雷林學(xué)森楊欣毅
    系統(tǒng)仿真技術(shù) 2016年2期
    關(guān)鍵詞:極限學(xué)習(xí)機(jī)航空發(fā)動機(jī)故障診斷

    王小雷, 林學(xué)森, 楊欣毅

    (1.中國人民解放軍91715部隊82分隊,廣東 廣州 510750;2.海軍航空工程學(xué)院 飛行器工程系,山東 煙臺 264001)

    ?

    ELM在航空發(fā)動機(jī)氣路部件故障診斷的應(yīng)用研究

    王小雷1,林學(xué)森2,楊欣毅2

    (1.中國人民解放軍91715部隊82分隊,廣東 廣州510750;2.海軍航空工程學(xué)院 飛行器工程系,山東 煙臺264001)

    摘要:針對航空發(fā)動機(jī)故障呈現(xiàn)復(fù)雜性、多樣性、非線性等特點,運用傳統(tǒng)的BP、ELman神經(jīng)網(wǎng)絡(luò)進(jìn)行發(fā)動機(jī)氣路部件故障診斷存在網(wǎng)絡(luò)訓(xùn)練參數(shù)設(shè)置復(fù)雜,迭代次數(shù)多,訓(xùn)練速度慢,泛化能力欠缺等問題。為此提出利用極限學(xué)習(xí)機(jī)(Extreme Learning Machine,ELM)法識別渦扇發(fā)動機(jī)氣路部件故障。該方法通過為輸入權(quán)和隱藏層偏置隨機(jī)賦值,利用MP逆求解輸出權(quán)值。以某型渦扇發(fā)動機(jī)為對象,進(jìn)行ELM、 BP、ELman氣路部件故障診斷比較研究,實驗結(jié)果驗證了利用ELM識別渦扇發(fā)動機(jī)氣路部件故障的精確性、快速性、穩(wěn)定性。

    關(guān)鍵詞:航空發(fā)動機(jī); 氣路部件; 故障診斷; 極限學(xué)習(xí)機(jī)(ELM)

    1引言

    航空發(fā)動機(jī)氣路部件故障對飛機(jī)影響巨大,在飛行過程中,航空發(fā)動機(jī)氣路部件一旦發(fā)生故障,后果不堪設(shè)想,輕則對飛機(jī)的飛行性能造成一定的影響,重則有可能機(jī)毀人亡,造成重大損失,所以對航空發(fā)動機(jī)氣路故障診斷一直受到高度重視[1]。

    目前,應(yīng)用比較廣泛的航空發(fā)動機(jī)氣路部件數(shù)據(jù)驅(qū)動的方法是神經(jīng)網(wǎng)絡(luò)法,而在神經(jīng)網(wǎng)絡(luò)中應(yīng)用較多的是BP神經(jīng)網(wǎng)絡(luò)和ELman神經(jīng)網(wǎng)絡(luò)。傳統(tǒng)的BP神經(jīng)網(wǎng)絡(luò)和ELman神經(jīng)網(wǎng)絡(luò)等方法構(gòu)造信度函數(shù)矩陣,網(wǎng)絡(luò)訓(xùn)練參數(shù)設(shè)置復(fù)雜,迭代次數(shù)多,訓(xùn)練速度慢,泛化能力欠缺,容易產(chǎn)生局部最優(yōu)解,存在一定的隨機(jī)性和片面性,可信度較差。為了提高發(fā)動機(jī)氣路部件的安全性能和可靠性能,必須采用先進(jìn)的技術(shù)手段,來預(yù)防或者提前檢測航空發(fā)動機(jī)氣路部件的故障。針對發(fā)動機(jī)故障多樣性和復(fù)雜性等特點,本研究提出了一種基于極限學(xué)習(xí)機(jī)(ELM)[2]的航空發(fā)動機(jī)氣路部件故障診斷方法,極限學(xué)習(xí)機(jī)據(jù)有批處理,不迭代,訓(xùn)練和測試速度快的特點。

    為了驗證ELM的優(yōu)越性,在本文中將分別利用ELM、BP和ELman三種方法對渦扇發(fā)動機(jī)氣路部件故障進(jìn)行診斷,并對診斷結(jié)果進(jìn)行分析。

    2ELM的理論分析

    誤差反向傳播神經(jīng)網(wǎng)絡(luò)[3](back-propagationneuralnetworks),簡稱BP神經(jīng)網(wǎng)絡(luò),它是一種反向傳遞并修正誤差的多層映射網(wǎng),它由輸入層、中間層和輸出層組成,中間層也就是隱層,可以使一層或多層。BP神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)過程分為兩部分組成:正向傳播和反向傳播。當(dāng)正向傳播時,信息從輸入層經(jīng)隱單元層處理后傳向輸出層,每一層的神經(jīng)元的狀態(tài)只影響下一層神經(jīng)元的狀態(tài)。如果在輸出層得不到希望的輸出,則轉(zhuǎn)入反向傳播,將誤差信號沿原來的神經(jīng)元連接通路返回。返回過程中,逐一修改各層神經(jīng)元連接的權(quán)值。這種過程不斷迭代,最后使誤差信號達(dá)到允許的范圍之內(nèi)。

    ELman神經(jīng)網(wǎng)絡(luò)[4]是一種典型的局部回歸網(wǎng)絡(luò)(globalfeedforwardlocalrecurrent),ELman神經(jīng)網(wǎng)絡(luò)一般分為四層:輸入層、中間層(隱含層)、承接層和輸出層,,其輸入層、隱含層和輸出層的連接類似于前饋網(wǎng)絡(luò),輸入層的單元僅起信號傳輸作用,輸出層單元起線性加權(quán)作用,隱含層單元的傳遞函數(shù)可采用線性或非線性函數(shù),承接層又稱為上下文層或狀態(tài)層,用來記憶隱含層單元前一時刻的輸出值,可以認(rèn)為是一個反饋。

    在極限學(xué)習(xí)機(jī)(ELM)中,輸入權(quán)值和隱含層偏差被隨機(jī)初始化給定,輸出權(quán)值矩陣?yán)脧V義逆(MP)計算得到。相比傳統(tǒng)前饋神經(jīng)網(wǎng)絡(luò),ELM學(xué)習(xí)速度更快、精度更高、參數(shù)調(diào)整簡單,已經(jīng)得到不少學(xué)者的關(guān)注研究,但在實際應(yīng)用過程中,為了達(dá)到理想的誤差精度,ELM通常需要大量的隱含層節(jié)點。

    (1)

    式中:L為隱層節(jié)點個數(shù);wj=[ωj1,ωj2,…,ωjn]T為連接第j個隱層節(jié)點和輸入節(jié)點的權(quán)重向量;β=[βj1,βj1,…,βjm]T為連接第j個隱層節(jié)點和輸出節(jié)點的權(quán)重向量;bj是第j個隱層節(jié)點的偏置。式(1)可簡寫為

    (2)

    (3)

    式中:H是隱含層神經(jīng)輸出矩陣,H的第i列是關(guān)于輸入x1,x2,…,xn的第i個神經(jīng)元的輸出。由于隱層輸入權(quán)值和閾值不需要調(diào)整可以任意賦值,當(dāng)給定ωi和bi后,訓(xùn)練單隱藏層前饋神經(jīng)網(wǎng)絡(luò)就等于求線性系統(tǒng)Hβ=T的最小二乘解β=H+T=(HTH)-1HTH,這里H+表示隱含層輸出矩陣{(xi,ti)|xi∈Rm,i=1,2,…,N}的MoorPenrose廣義逆。

    ELM的學(xué)習(xí)算法為:給定一個訓(xùn)練樣{(xi,ti)|xi∈Rm,i=1,2,…,N},激活函數(shù)g(x),隱層神經(jīng)元L,具體步驟如下:

    步驟一:隨機(jī)設(shè)定隱層的權(quán)值ωi和閾值bi,i=1,2,…,L。

    步驟二:計算隱層輸出矩陣H。

    步驟三:計算隱層矩陣H的逆矩陣。

    步驟四:計算輸出層權(quán)值β:

    (4)

    其中:T=[t1,t2,…,tN]T。

    3航空發(fā)動機(jī)氣路故障診斷

    3.1渦扇發(fā)動機(jī)氣路部件故障診斷原理[5]

    航空發(fā)動機(jī)氣路部件主要分為:低壓壓器機(jī)、高壓壓氣機(jī)、燃燒室、高壓渦輪、低壓渦輪和尾噴管等。

    在這些部件之上,可以測量的參數(shù)有:低壓壓氣機(jī)轉(zhuǎn)速(低壓渦輪轉(zhuǎn)速)、高壓壓氣機(jī)轉(zhuǎn)速(高壓渦輪轉(zhuǎn)速)、高壓壓氣機(jī)后壓力、低壓渦輪后壓力、高壓壓氣機(jī)后溫度、低壓渦輪出口溫度、尾噴管出口溫度、耗油量等。一般將這些參數(shù)作為標(biāo)致航空發(fā)動機(jī)是否正常的重要參數(shù),在正常工作的情況下,這些參數(shù)都是穩(wěn)定在一定的值之內(nèi)的,當(dāng)這些參數(shù)超出正常工作的參數(shù)范圍時,則表明該部件出現(xiàn)故障,如圖1所示。

    圖1 故障判斷原理Fig.1 Diagram of fault decision

    3.2渦扇發(fā)動機(jī)氣路部件故障診斷

    為了研究采用極限學(xué)習(xí)機(jī)(ELM)判斷航空發(fā)動機(jī)氣路部件故障,本文將ELM與BP神經(jīng)網(wǎng)絡(luò)和ELman神經(jīng)網(wǎng)絡(luò)進(jìn)行比較,并驗證極限學(xué)習(xí)機(jī)在航空發(fā)動機(jī)氣路部件故障診斷中的優(yōu)勢,利用航空發(fā)動機(jī)氣路部件模擬軟件GSP[6],選擇渦扇發(fā)動機(jī)氣路部件模型,選取10個指標(biāo)作為數(shù)據(jù)測量點,如表1所示。

    表1 測量的參數(shù)Tab.1 Measurement Parameters

    在測量這10個指標(biāo)時,選擇衰減步長為-0.001,衰減范圍為[-0.001,-0.05],分別生成風(fēng)扇(低壓壓氣機(jī))的流量衰減數(shù)據(jù)50組,效率衰減數(shù)據(jù)50組,壓比衰減數(shù)據(jù)50組;高壓壓氣機(jī)的流量衰減數(shù)據(jù)50組,效率衰減數(shù)據(jù)50組,壓比衰減數(shù)據(jù)50組;高壓渦輪的流量衰減數(shù)據(jù)50組,效率衰減數(shù)據(jù)50組;低壓渦輪的流量衰減數(shù)據(jù)50組,效率衰減數(shù)據(jù)50組,總共生成500組衰減數(shù)據(jù),也就是得到500組有衰減的情況下,這10個指標(biāo)的測量數(shù)據(jù)。然后,生成一組在沒有衰減情況下,這10個指標(biāo)的測量數(shù)據(jù),也就是得到10個正常的的測量數(shù)據(jù),與之前的500組衰減數(shù)據(jù)共同構(gòu)成501組訓(xùn)練數(shù)據(jù)。

    將測量所得到的每一項數(shù)據(jù)進(jìn)行歸一化處理,如高壓壓氣機(jī)轉(zhuǎn)速NH,用高壓壓氣機(jī)測量轉(zhuǎn)速減去高壓壓氣機(jī)最小轉(zhuǎn)速所得的絕對值,除以高壓壓氣機(jī)最大轉(zhuǎn)速減去高壓壓氣機(jī)最小轉(zhuǎn)速,即

    (5)

    同理,分別對低壓壓氣機(jī)轉(zhuǎn)速NL,高壓壓氣機(jī)后壓力P3,低壓渦輪后壓力P6,高壓壓氣機(jī)后溫度T3,低壓渦輪出口溫度T6,尾噴管出口溫度T8,尾噴管出口壓力P8,發(fā)動機(jī)推力FN和耗油量Tsfc進(jìn)行歸一化處理。

    經(jīng)過歸一化處理后得到一個501組這10個指標(biāo)的數(shù)據(jù),即歸一化后500組衰減數(shù)據(jù),1組正常數(shù)據(jù),將前500組數(shù)據(jù)中,同一種測量指標(biāo)中采用的相同衰減方式的標(biāo)上同一個序號,依次標(biāo)上1,2,3,4,5,6,7,8,9,10,作為訓(xùn)練數(shù)據(jù)。

    通過實驗驗證得到ELM的測試精度與隱含層節(jié)點[7]的關(guān)系如圖2所示。

    圖2 ELM的測試精度與隱含層節(jié)點的關(guān)系Fig.2 Testing accuracy varies with the number of ELM hidden nodes

    分析可知當(dāng)隱含層節(jié)點為35時,測試精度最高為94.03%,所以在故障診斷中選取ELM的隱含層節(jié)點為35個。

    該數(shù)據(jù)均在MatlabR2013a環(huán)境下運行的,運行電腦處理器為AMDA6-3420MAPU,1.5GHz,內(nèi)存為4.00GB,64位操作系統(tǒng)。

    利用501組訓(xùn)練數(shù)據(jù)分別對ELM、BP和ELman進(jìn)行訓(xùn)練,然后對50組測試數(shù)據(jù)進(jìn)行測試,分別訓(xùn)練和測試20次,得到結(jié)果如圖3所示。

    圖3 訓(xùn)練時間比較圖Fig.3 Comparison of training time

    由圖3分析可知,ELM的訓(xùn)練時間最短,且訓(xùn)練時間波動最小,為[0.000 0,0.405 6],BP的訓(xùn)練時間最長,訓(xùn)練時間的波動最大為[17.752 9,92.227 8],ELman的訓(xùn)練時間和波動在兩者之間,波動范圍為[9.188 5,24.320 6]。

    圖4 測試精度比較圖Fig.4 Comparison of training testing accuracy

    由圖4分析可知,ELM,BP,ELman這三種方法的測試精度非常接近,其中ELM和BP的測試精度最高,且最穩(wěn)定,波動范圍為[0.980 0,1.000 0]和[0.900 0,1.000 0],而ELman的波動范圍最大,為[0.720 0,0.980 0]。

    圖5 測試時間比較圖Fig.5 Comparison of training testing time

    由圖5分析可以明顯發(fā)現(xiàn),ELM的測試時間最短,波動最小,為[0.000 0,0.124 8],其次是ELman的測試時間較短,波動為[0.124 8,0.327 6],BP的測試時間最長,波動最大,為[0.234 0,0.468 0]。

    將501組訓(xùn)練數(shù)據(jù)和50組測試數(shù)據(jù)分別用ELM,BP,ELman這三種方法進(jìn)行訓(xùn)練和測試,且采用三種方法各分別訓(xùn)練和測試20次,得到訓(xùn)練時間、測試精度和測試時間如表2所示。

    分析可知,ELM的訓(xùn)練時間最短,平均訓(xùn)練時間為0.099 8秒,且訓(xùn)練時間波動最小,訓(xùn)練時間方差為0.007 0;BP的訓(xùn)練時間最長,訓(xùn)練時間為ELM的411倍,訓(xùn)練時間的波動最大,方差為ELM的41390倍;ELman的訓(xùn)練時間和波動在兩者之間,訓(xùn)練時間為ELM的212倍,方差為ELM的1 333倍。由此可得,ELM的訓(xùn)練時間最短且最穩(wěn)定。

    表2 ELM、BP和ELman三種方法比較Tab.2 Performance comparison of ELM,BP and ELman

    ELM,BP,ELman這三種方法的測試精度非常接近,其中BP和ELM的測試精度最高,且最穩(wěn)定。

    分析可以發(fā)現(xiàn),ELM的測試時間最短,平均測試時間為0.042 9秒,測試時間方差為0.001 6;其次是ELman的測試時間較短0.179 4,為ELM的4倍,測試時間方差為0.001 6,與ELM相同;BP的測試時間最長,為0.269 1秒,為ELM的6倍,測試時間方差為0.002 5,為ELM的1.56倍。綜合可知ELM的測試時間最短且最穩(wěn)定。

    4結(jié)論

    通過應(yīng)用極限學(xué)習(xí)機(jī)對渦扇發(fā)動機(jī)氣路部件的故障診斷,并與應(yīng)用BP和ELman進(jìn)行渦扇發(fā)動機(jī)氣路部件故障診斷進(jìn)行比較,得出以下結(jié)論:

    (1) 運用ELM法進(jìn)行故障訓(xùn)練時,能夠在保持方差低于BP和ELman的前提下,大大縮短訓(xùn)練時間,平均訓(xùn)練時間為BP的平均訓(xùn)練時間的1/411,為ELman的平均訓(xùn)練時間的1/212,有利于進(jìn)行大數(shù)據(jù)訓(xùn)練時提高精度并減少時間。

    (2) 運用ELM法進(jìn)行故障診斷測試時,在保持比ELman精度高,與BP精度幾乎一樣的前提下,測試時間大大縮短了,為BP的1/6,為ELman的1/4,有利于在進(jìn)行大數(shù)據(jù)故障測試時節(jié)約大量時間,且能夠保持穩(wěn)定的精度。

    (3) 運用ELM法進(jìn)行故障測試時,與BP、ELman比較,多次測試所得的方差,即訓(xùn)練時間方差、測試精度方差和測試時間方差均是最小,模型的穩(wěn)定性最好。

    參考文獻(xiàn):

    [1]李應(yīng)紅,尉詢楷.航空發(fā)動機(jī)的智能診斷、建模與預(yù)測方法[M].北京:科學(xué)出版社,2013.

    LIYinghong,WEIXunkai.Intelligentdiagnosis,modelingandpredicationmethodforaeroengine[M].Beijing:SciencePress,2013.

    [2]HUANGGuangbin,ZHUQinyu,SiewCheekheong.Extremelearningmachine:theoryandapplications[J].Neurocomputing,2006,70(1/2/3):489- 501.

    [3]田雨波.混合神經(jīng)網(wǎng)絡(luò)技術(shù)[M].北京:科學(xué)出版社,2009.

    TIANYubo.Technologyonmixneuralnetwork[M].Beijing:SciencePress,2009.

    [4]任麗娜.基于Elman神經(jīng)網(wǎng)絡(luò)的中期電力負(fù)荷預(yù)測模型研究[D].蘭州理工大學(xué),2007.

    RENLina.MiddletermelectricpowerloadpredictionbasedonElmanneuralnetwork[D].Lanzhouuniversityofscienceandengineering,2007.

    [5]范作民,孫春林,白杰.航空發(fā)動機(jī)故障診斷導(dǎo)論[M].北京:科學(xué)出版社,2004.

    FANZuoming,SUNChunlin,Baijie.Introductiontoaeroenginefaultdiagnosis[M].Beijing:SciencePress,2004.[6]VisserWPJ,KogenhopO.Aqenerieapproachforgasturbineadaptivemodeling,[R],NLR-TP-2004-391.

    [7]張德豐.MATLAB神經(jīng)網(wǎng)絡(luò)應(yīng)用設(shè)計[M].北京:機(jī)械工業(yè)出版社,2011.

    ZHANGDefeng.ApplicationsanddesignofNeuralnetworkusingMatlab[M].Beijing:MechanicalIndustryPress,2011.

    王小雷男(1989-),湖南萊陽人,學(xué)士,主要研究方向為發(fā)動機(jī)維修與故障診斷。

    林學(xué)森男(1990-),山東濰坊人,博士研究生,主要研究方向為發(fā)動機(jī)狀態(tài)監(jiān)控與故障診斷。

    中圖分類號:TP 183

    文獻(xiàn)標(biāo)識碼:A

    Research on ELM Application in Aeroengine gas Path Components Fault Diagnosis

    WANG Xiaolei1,LIN Xuesen2,YANG Yinyi2

    (1.ChinaPLA91715Forces,82Unit,Guangzhou510750China;2.NavalAeronauticalEngineeringInstitute,Yantai264001,China)

    Abstract:The aeroengine gas path component faults have the characteristics of complexity,diversity and nonlinear.When applied traditional BP and ELman into aeroengine gas path components fault diagnosis,there exist some problems as complex setting of network parameters,multiple iterations,slow training speed and poor generation performance.In order to solve these problems,Extreme Learning Machine(Extreme Learning Machine,ELM)method was used to aeroengine gas path components fault diagnosis.This method randomly assigned input weights and biases and calculate output weights by Moore-penrose,pseudo- inverse.ELM,BP and ELman were studied and compared on a certain type of turbofan engine.The results verify the diagnostic accuracy,quickness and stability of ELM.

    Key words:aircraft engine; gas path components; fault diagnosis; extreme learning machine(ELM)

    猜你喜歡
    極限學(xué)習(xí)機(jī)航空發(fā)動機(jī)故障診斷
    基于判別信息極限學(xué)習(xí)機(jī)的高光譜遙感圖像分類
    基于極限學(xué)習(xí)機(jī)的玻璃瓶口缺陷檢測方法研究
    某型發(fā)動機(jī)喘振故障的研究
    “航空發(fā)動機(jī)強(qiáng)度與振動”教學(xué)改革的探索與思考
    基于PDM的航空發(fā)動機(jī)技術(shù)狀態(tài)管理
    淺談航空發(fā)動機(jī)起動與點火系統(tǒng)
    科技視界(2016年14期)2016-06-08 00:41:57
    極限學(xué)習(xí)機(jī)在圖像分割中的應(yīng)用
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    基于關(guān)節(jié)信息和極限學(xué)習(xí)機(jī)的人體動作識別
    基于LCD和排列熵的滾動軸承故障診斷
    国产伦精品一区二区三区四那| 久久久久久久亚洲中文字幕| 大香蕉97超碰在线| 丁香六月天网| 日韩欧美精品免费久久| 美女主播在线视频| 不卡视频在线观看欧美| 麻豆精品久久久久久蜜桃| 黄色怎么调成土黄色| 亚洲av欧美aⅴ国产| 久久狼人影院| 热re99久久国产66热| 只有这里有精品99| 久久久久久久久久成人| 亚洲精品第二区| 成人漫画全彩无遮挡| 日本-黄色视频高清免费观看| 少妇精品久久久久久久| 中文乱码字字幕精品一区二区三区| 亚洲av.av天堂| 高清午夜精品一区二区三区| 夜夜爽夜夜爽视频| 好男人视频免费观看在线| 一区在线观看完整版| 波野结衣二区三区在线| 亚洲第一av免费看| 成人二区视频| 菩萨蛮人人尽说江南好唐韦庄| 晚上一个人看的免费电影| 视频中文字幕在线观看| 九九久久精品国产亚洲av麻豆| 亚洲欧美精品自产自拍| 国产黄片美女视频| 久久婷婷青草| 久久热精品热| 成年女人在线观看亚洲视频| 黄色日韩在线| 国产av一区二区精品久久| av线在线观看网站| 黄色视频在线播放观看不卡| 天堂8中文在线网| www.色视频.com| 成人毛片a级毛片在线播放| 久久这里有精品视频免费| 男人爽女人下面视频在线观看| 少妇人妻久久综合中文| 精品熟女少妇av免费看| 18+在线观看网站| 最近2019中文字幕mv第一页| 久久6这里有精品| 日韩中文字幕视频在线看片| 日本与韩国留学比较| 在线天堂最新版资源| 日韩,欧美,国产一区二区三区| 久久久久久久亚洲中文字幕| 我要看黄色一级片免费的| 亚洲怡红院男人天堂| tube8黄色片| 亚洲精品日本国产第一区| 精品久久国产蜜桃| 亚洲欧美精品专区久久| 黄色一级大片看看| 日韩一区二区视频免费看| 国产欧美日韩一区二区三区在线 | 高清视频免费观看一区二区| 国产日韩欧美亚洲二区| 狂野欧美激情性bbbbbb| 国产极品粉嫩免费观看在线 | 伊人久久国产一区二区| 秋霞在线观看毛片| 国产精品嫩草影院av在线观看| 欧美国产精品一级二级三级 | 国产片特级美女逼逼视频| av卡一久久| 欧美日韩国产mv在线观看视频| 亚洲情色 制服丝袜| 我的女老师完整版在线观看| 男人舔奶头视频| 国产真实伦视频高清在线观看| 日本av免费视频播放| 久久久国产一区二区| 亚洲成人手机| 下体分泌物呈黄色| 亚洲精品一区蜜桃| 在线精品无人区一区二区三| 午夜福利,免费看| 黄色配什么色好看| 亚洲精品日韩在线中文字幕| www.色视频.com| 搡老乐熟女国产| 久久韩国三级中文字幕| av女优亚洲男人天堂| 午夜精品国产一区二区电影| 亚洲三级黄色毛片| 日韩在线高清观看一区二区三区| 高清午夜精品一区二区三区| 大话2 男鬼变身卡| 少妇被粗大猛烈的视频| 草草在线视频免费看| 亚洲av成人精品一区久久| 日本欧美国产在线视频| 欧美日韩视频高清一区二区三区二| 五月玫瑰六月丁香| 色婷婷av一区二区三区视频| 亚洲精品日韩av片在线观看| 深夜a级毛片| av免费观看日本| 一级,二级,三级黄色视频| 国产免费福利视频在线观看| 精品国产露脸久久av麻豆| 国产成人一区二区在线| 精品一品国产午夜福利视频| 又黄又爽又刺激的免费视频.| 中国美白少妇内射xxxbb| 中文字幕免费在线视频6| a 毛片基地| 狂野欧美白嫩少妇大欣赏| 少妇裸体淫交视频免费看高清| 亚洲av.av天堂| 亚洲国产精品专区欧美| 少妇丰满av| 3wmmmm亚洲av在线观看| 五月伊人婷婷丁香| 国产欧美日韩综合在线一区二区 | 精品久久久久久久久亚洲| 亚洲激情五月婷婷啪啪| av.在线天堂| av天堂中文字幕网| 熟女av电影| videos熟女内射| 色视频www国产| 亚洲va在线va天堂va国产| 国产亚洲91精品色在线| 色94色欧美一区二区| 国产成人精品福利久久| 黑丝袜美女国产一区| 久久av网站| 国产成人一区二区在线| 午夜福利,免费看| 国产欧美日韩精品一区二区| 黄色欧美视频在线观看| 丝袜喷水一区| 如何舔出高潮| 在线观看www视频免费| 一区二区三区乱码不卡18| 精品久久久精品久久久| 毛片一级片免费看久久久久| 久久久久久久精品精品| 夫妻午夜视频| 欧美少妇被猛烈插入视频| 桃花免费在线播放| 又大又黄又爽视频免费| 日本av免费视频播放| 男人舔奶头视频| 午夜福利网站1000一区二区三区| 精品熟女少妇av免费看| 成人免费观看视频高清| 亚洲人成网站在线观看播放| 欧美97在线视频| 精品人妻熟女av久视频| 大香蕉97超碰在线| 好男人视频免费观看在线| av在线播放精品| 色视频在线一区二区三区| 久久人人爽av亚洲精品天堂| 国产亚洲最大av| 18禁在线播放成人免费| 草草在线视频免费看| 精品国产乱码久久久久久小说| 久久久亚洲精品成人影院| 一级毛片黄色毛片免费观看视频| 成年人免费黄色播放视频 | 亚洲电影在线观看av| 天天操日日干夜夜撸| 国产永久视频网站| 天堂中文最新版在线下载| 高清午夜精品一区二区三区| 国产永久视频网站| 国产日韩欧美在线精品| av国产久精品久网站免费入址| 日本午夜av视频| 国产精品.久久久| 国产色婷婷99| 日本黄色片子视频| 日日爽夜夜爽网站| 99久久精品热视频| 精品国产一区二区三区久久久樱花| 日本黄色片子视频| 99九九线精品视频在线观看视频| 如日韩欧美国产精品一区二区三区 | 精品国产一区二区久久| 日韩在线高清观看一区二区三区| 亚洲精品第二区| 国产男女内射视频| 国产成人精品无人区| 日韩电影二区| 久久久久久久亚洲中文字幕| 一个人免费看片子| 一级爰片在线观看| 成年女人在线观看亚洲视频| 美女中出高潮动态图| 美女视频免费永久观看网站| 夜夜骑夜夜射夜夜干| 亚洲,一卡二卡三卡| 精品久久久久久久久av| 日韩中字成人| 国产 一区精品| 欧美bdsm另类| 最近最新中文字幕免费大全7| 极品人妻少妇av视频| 国产精品秋霞免费鲁丝片| 最新的欧美精品一区二区| 婷婷色麻豆天堂久久| 日产精品乱码卡一卡2卡三| 欧美亚洲 丝袜 人妻 在线| 亚洲精品亚洲一区二区| 国产黄频视频在线观看| 婷婷色综合大香蕉| 久久人人爽人人爽人人片va| 欧美少妇被猛烈插入视频| 精品少妇久久久久久888优播| 91精品伊人久久大香线蕉| 热99国产精品久久久久久7| 下体分泌物呈黄色| 国产精品久久久久久精品古装| 免费在线观看成人毛片| 免费黄色在线免费观看| av国产久精品久网站免费入址| av天堂久久9| √禁漫天堂资源中文www| 亚洲欧洲日产国产| 亚洲国产最新在线播放| 久久精品国产亚洲av天美| 在线播放无遮挡| 午夜福利网站1000一区二区三区| 高清av免费在线| 亚洲情色 制服丝袜| 插阴视频在线观看视频| 国产精品一二三区在线看| 亚洲丝袜综合中文字幕| 天天操日日干夜夜撸| 亚洲精品久久久久久婷婷小说| 亚洲经典国产精华液单| 国产精品.久久久| av天堂中文字幕网| 免费观看av网站的网址| 国产免费福利视频在线观看| 久久热精品热| 国产真实伦视频高清在线观看| 日韩精品有码人妻一区| www.色视频.com| 日日撸夜夜添| 日韩三级伦理在线观看| 久久久久人妻精品一区果冻| 卡戴珊不雅视频在线播放| 黄色怎么调成土黄色| 777米奇影视久久| 亚洲欧美精品专区久久| 亚洲av日韩在线播放| 国产伦理片在线播放av一区| 亚洲丝袜综合中文字幕| 99久久人妻综合| 99热这里只有是精品在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲国产精品一区三区| xxx大片免费视频| 日本wwww免费看| 国产av一区二区精品久久| 欧美最新免费一区二区三区| 汤姆久久久久久久影院中文字幕| 精品国产一区二区久久| 综合色丁香网| 在线免费观看不下载黄p国产| 日本91视频免费播放| 亚洲第一区二区三区不卡| 亚洲精品国产色婷婷电影| 国产老妇伦熟女老妇高清| 国产免费又黄又爽又色| 99视频精品全部免费 在线| 久久精品久久久久久噜噜老黄| 热re99久久国产66热| 老熟女久久久| 亚洲成人av在线免费| 日本wwww免费看| 精品少妇内射三级| 久久精品国产亚洲网站| 少妇被粗大的猛进出69影院 | 亚洲久久久国产精品| 看非洲黑人一级黄片| 久久午夜福利片| 国产色爽女视频免费观看| 成人毛片60女人毛片免费| 丰满乱子伦码专区| 亚洲国产色片| 国产高清三级在线| 我的老师免费观看完整版| 亚洲自偷自拍三级| 欧美 亚洲 国产 日韩一| 中文字幕免费在线视频6| 国产高清不卡午夜福利| 26uuu在线亚洲综合色| 久久韩国三级中文字幕| 国产精品久久久久久精品古装| 精品亚洲乱码少妇综合久久| 亚洲,欧美,日韩| 欧美日韩视频精品一区| √禁漫天堂资源中文www| 国产亚洲一区二区精品| 久久亚洲国产成人精品v| 亚洲国产精品国产精品| 午夜福利,免费看| .国产精品久久| 精品国产一区二区久久| 在线观看美女被高潮喷水网站| 久久人人爽av亚洲精品天堂| 天堂俺去俺来也www色官网| 秋霞伦理黄片| 国产白丝娇喘喷水9色精品| 黄片无遮挡物在线观看| 一边亲一边摸免费视频| 热99国产精品久久久久久7| 亚洲第一av免费看| 天堂中文最新版在线下载| 夫妻午夜视频| 熟女人妻精品中文字幕| 精品亚洲成国产av| kizo精华| 亚洲精品久久久久久婷婷小说| 一级片'在线观看视频| 下体分泌物呈黄色| av黄色大香蕉| 在线观看人妻少妇| 国产精品一区二区三区四区免费观看| 日本av免费视频播放| 国产精品女同一区二区软件| 另类精品久久| 在线精品无人区一区二区三| 国产熟女午夜一区二区三区 | 少妇人妻久久综合中文| 国产免费视频播放在线视频| 香蕉精品网在线| 老司机影院毛片| 精品午夜福利在线看| 中文字幕精品免费在线观看视频 | 国产黄频视频在线观看| 欧美少妇被猛烈插入视频| 成年av动漫网址| www.av在线官网国产| 久久国产乱子免费精品| 亚洲国产精品一区三区| 麻豆乱淫一区二区| 亚洲在久久综合| 少妇人妻 视频| 亚洲高清免费不卡视频| 亚洲av成人精品一区久久| 熟女电影av网| 一本—道久久a久久精品蜜桃钙片| 亚洲色图综合在线观看| 日韩大片免费观看网站| av女优亚洲男人天堂| 亚洲精华国产精华液的使用体验| 国产熟女午夜一区二区三区 | 成人特级av手机在线观看| 国产精品99久久99久久久不卡 | 寂寞人妻少妇视频99o| 人体艺术视频欧美日本| 九色成人免费人妻av| 麻豆成人午夜福利视频| 日本欧美视频一区| 日韩熟女老妇一区二区性免费视频| 亚洲av成人精品一二三区| 国产黄片视频在线免费观看| 国产乱来视频区| 人人妻人人爽人人添夜夜欢视频 | 日日撸夜夜添| 18禁在线无遮挡免费观看视频| 久久久久久久久久久丰满| 一级毛片黄色毛片免费观看视频| 亚洲美女视频黄频| 人妻人人澡人人爽人人| 在线观看一区二区三区激情| 插逼视频在线观看| 国产一区亚洲一区在线观看| 国产精品国产三级国产av玫瑰| 成人无遮挡网站| 亚洲精品456在线播放app| 国产在线免费精品| 草草在线视频免费看| 国产乱来视频区| 精品国产露脸久久av麻豆| 一区二区三区四区激情视频| 亚洲国产毛片av蜜桃av| 男人舔奶头视频| 午夜福利视频精品| 精品久久久噜噜| 国产一区二区三区av在线| 久久久精品免费免费高清| 日本免费在线观看一区| 自线自在国产av| 国产白丝娇喘喷水9色精品| 一级毛片 在线播放| 亚洲va在线va天堂va国产| 老女人水多毛片| 国产精品久久久久久久电影| 国产熟女欧美一区二区| 亚洲欧洲国产日韩| 热re99久久国产66热| 观看免费一级毛片| 日韩av在线免费看完整版不卡| 欧美国产精品一级二级三级 | 日韩,欧美,国产一区二区三区| 国产一区二区在线观看日韩| 久久韩国三级中文字幕| 又黄又爽又刺激的免费视频.| 2022亚洲国产成人精品| 最近中文字幕高清免费大全6| 色婷婷久久久亚洲欧美| 美女中出高潮动态图| 亚洲欧美精品专区久久| 最近中文字幕高清免费大全6| 内地一区二区视频在线| 边亲边吃奶的免费视频| 午夜久久久在线观看| 综合色丁香网| 91精品国产九色| 夫妻性生交免费视频一级片| 国产亚洲91精品色在线| 亚洲av综合色区一区| 黄色日韩在线| 国产成人91sexporn| 99久久综合免费| 黄色视频在线播放观看不卡| 日韩精品免费视频一区二区三区 | 99热国产这里只有精品6| √禁漫天堂资源中文www| 免费看日本二区| 亚洲电影在线观看av| 久久影院123| 十八禁高潮呻吟视频 | 国产片特级美女逼逼视频| 亚洲性久久影院| 日韩中文字幕视频在线看片| 另类精品久久| 制服丝袜香蕉在线| 精品少妇久久久久久888优播| 国产免费又黄又爽又色| 中文欧美无线码| 国产精品久久久久成人av| 人人妻人人爽人人添夜夜欢视频 | 一区二区三区免费毛片| 精品一品国产午夜福利视频| 哪个播放器可以免费观看大片| 国产精品久久久久成人av| 秋霞伦理黄片| 午夜福利网站1000一区二区三区| 国精品久久久久久国模美| av专区在线播放| 3wmmmm亚洲av在线观看| 亚洲综合精品二区| 亚洲欧美成人精品一区二区| 亚洲欧洲国产日韩| 免费黄色在线免费观看| 肉色欧美久久久久久久蜜桃| 在线观看美女被高潮喷水网站| 伊人久久国产一区二区| 日韩一本色道免费dvd| 免费观看无遮挡的男女| 国产伦理片在线播放av一区| 五月玫瑰六月丁香| 国产一级毛片在线| 自拍偷自拍亚洲精品老妇| 春色校园在线视频观看| 午夜视频国产福利| xxx大片免费视频| 成人影院久久| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品中文字幕在线视频 | 91久久精品电影网| 91久久精品国产一区二区成人| 人人妻人人爽人人添夜夜欢视频 | 欧美 亚洲 国产 日韩一| 亚洲精华国产精华液的使用体验| 成人国产麻豆网| 婷婷色av中文字幕| 日韩一本色道免费dvd| 少妇猛男粗大的猛烈进出视频| 免费看光身美女| 亚洲精品日韩av片在线观看| 2018国产大陆天天弄谢| 黄色视频在线播放观看不卡| 精品一品国产午夜福利视频| 我的老师免费观看完整版| 一级片'在线观看视频| 老司机亚洲免费影院| 国产一区有黄有色的免费视频| 精品熟女少妇av免费看| 国产欧美日韩综合在线一区二区 | 有码 亚洲区| 下体分泌物呈黄色| 在线观看一区二区三区激情| 亚洲成人av在线免费| 成人国产麻豆网| 亚洲欧美日韩卡通动漫| 在线观看www视频免费| 97超碰精品成人国产| 精品人妻偷拍中文字幕| 各种免费的搞黄视频| 性色avwww在线观看| 乱码一卡2卡4卡精品| .国产精品久久| 日本wwww免费看| 精品卡一卡二卡四卡免费| 欧美 日韩 精品 国产| 麻豆成人av视频| 26uuu在线亚洲综合色| 下体分泌物呈黄色| 成人免费观看视频高清| 男女国产视频网站| 久热这里只有精品99| 黑人高潮一二区| 又爽又黄a免费视频| 又大又黄又爽视频免费| 亚洲国产欧美日韩在线播放 | 黄色毛片三级朝国网站 | 国产精品免费大片| 十八禁高潮呻吟视频 | 人妻少妇偷人精品九色| .国产精品久久| 国产av国产精品国产| 午夜福利,免费看| 午夜激情久久久久久久| 欧美日韩亚洲高清精品| 精品国产露脸久久av麻豆| 国产免费一区二区三区四区乱码| 蜜桃在线观看..| 高清av免费在线| 中文在线观看免费www的网站| 亚洲欧美精品专区久久| 在线 av 中文字幕| 色94色欧美一区二区| 久久午夜综合久久蜜桃| 欧美一级a爱片免费观看看| 亚洲精品,欧美精品| 国产日韩欧美视频二区| 午夜福利网站1000一区二区三区| 高清欧美精品videossex| 人人妻人人澡人人爽人人夜夜| 观看美女的网站| 80岁老熟妇乱子伦牲交| 免费不卡的大黄色大毛片视频在线观看| 狠狠精品人妻久久久久久综合| 中文欧美无线码| 三级国产精品片| 亚洲精品视频女| 中文乱码字字幕精品一区二区三区| 三级经典国产精品| a级毛片免费高清观看在线播放| 在线观看人妻少妇| 久久久久久久久大av| 晚上一个人看的免费电影| 亚洲精品乱码久久久久久按摩| 国产精品久久久久成人av| 成人亚洲精品一区在线观看| 久久综合国产亚洲精品| av视频免费观看在线观看| 黑人高潮一二区| 又爽又黄a免费视频| 最新中文字幕久久久久| 啦啦啦啦在线视频资源| freevideosex欧美| 在现免费观看毛片| 美女脱内裤让男人舔精品视频| 日日摸夜夜添夜夜爱| www.av在线官网国产| 一二三四中文在线观看免费高清| 不卡视频在线观看欧美| 在线观看免费高清a一片| 三级经典国产精品| 王馨瑶露胸无遮挡在线观看| av.在线天堂| 最近手机中文字幕大全| 夫妻性生交免费视频一级片| 久久久久久久精品精品| 亚洲久久久国产精品| 日韩成人av中文字幕在线观看| 大又大粗又爽又黄少妇毛片口| 国产男人的电影天堂91| 欧美97在线视频| 美女国产视频在线观看| 中文字幕免费在线视频6| 国产黄频视频在线观看| 国产亚洲一区二区精品| 欧美亚洲 丝袜 人妻 在线| 日韩 亚洲 欧美在线| 美女大奶头黄色视频| 亚洲四区av| 熟女电影av网| 欧美日韩国产mv在线观看视频| 国产黄片美女视频| 午夜影院在线不卡| 久久国产精品大桥未久av | 少妇人妻一区二区三区视频| 99热这里只有是精品在线观看| 高清黄色对白视频在线免费看 | 久久久久久久久久人人人人人人| 午夜91福利影院| 色网站视频免费| 插阴视频在线观看视频| 人体艺术视频欧美日本| 午夜老司机福利剧场| av天堂中文字幕网| 国产成人a∨麻豆精品| 伦理电影大哥的女人| 免费观看无遮挡的男女| 岛国毛片在线播放|