• <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和排列熵的滾動軸承故障診斷
    少妇 在线观看| 岛国在线观看网站| 嫁个100分男人电影在线观看| 99在线人妻在线中文字幕| 亚洲熟妇中文字幕五十中出 | 亚洲第一青青草原| 欧美日韩亚洲高清精品| 亚洲欧美激情在线| 久热爱精品视频在线9| 亚洲中文日韩欧美视频| 久久久久久久久中文| 国产成人免费无遮挡视频| 可以在线观看毛片的网站| √禁漫天堂资源中文www| 国产成人啪精品午夜网站| 一区二区三区精品91| av超薄肉色丝袜交足视频| 欧美在线一区亚洲| 日本黄色日本黄色录像| 午夜a级毛片| 99国产极品粉嫩在线观看| 自线自在国产av| 免费在线观看影片大全网站| 精品久久久久久,| 黄色女人牲交| 91麻豆精品激情在线观看国产 | 日韩高清综合在线| 亚洲性夜色夜夜综合| 长腿黑丝高跟| 欧美老熟妇乱子伦牲交| 亚洲美女黄片视频| 亚洲欧洲精品一区二区精品久久久| 久久这里只有精品19| 亚洲成国产人片在线观看| 精品欧美一区二区三区在线| 成人国语在线视频| 在线观看免费午夜福利视频| 在线观看免费日韩欧美大片| 香蕉丝袜av| 淫秽高清视频在线观看| 宅男免费午夜| 国产三级在线视频| 黑人操中国人逼视频| 精品国产乱子伦一区二区三区| 亚洲一区二区三区欧美精品| 亚洲欧美一区二区三区黑人| 国产一卡二卡三卡精品| 亚洲第一av免费看| 侵犯人妻中文字幕一二三四区| 黄色成人免费大全| 欧美日韩视频精品一区| 中文字幕人妻丝袜制服| 丝袜人妻中文字幕| 欧美色视频一区免费| 日韩av在线大香蕉| 女生性感内裤真人,穿戴方法视频| 亚洲视频免费观看视频| 午夜免费观看网址| www.自偷自拍.com| 精品国产亚洲在线| 欧美激情久久久久久爽电影 | 999久久久精品免费观看国产| 国产精品乱码一区二三区的特点 | 韩国精品一区二区三区| 少妇粗大呻吟视频| 变态另类成人亚洲欧美熟女 | 交换朋友夫妻互换小说| 亚洲熟妇熟女久久| 亚洲国产毛片av蜜桃av| 精品久久久精品久久久| 精品久久久久久成人av| 这个男人来自地球电影免费观看| 国产免费av片在线观看野外av| 乱人伦中国视频| 少妇 在线观看| 亚洲一区中文字幕在线| 亚洲情色 制服丝袜| 成人18禁在线播放| 国产精品永久免费网站| 久久伊人香网站| svipshipincom国产片| 中出人妻视频一区二区| 日本撒尿小便嘘嘘汇集6| 亚洲国产欧美一区二区综合| 国产三级黄色录像| 亚洲第一青青草原| 99re在线观看精品视频| av欧美777| 亚洲精品成人av观看孕妇| 岛国在线观看网站| 咕卡用的链子| 一本大道久久a久久精品| 一级毛片女人18水好多| 国产激情久久老熟女| 欧美人与性动交α欧美精品济南到| 高清毛片免费观看视频网站 | 别揉我奶头~嗯~啊~动态视频| 久久天躁狠狠躁夜夜2o2o| 一夜夜www| 亚洲人成77777在线视频| 国产成人精品久久二区二区91| 久久精品成人免费网站| 一夜夜www| 久久久久久久午夜电影 | 精品福利永久在线观看| www.999成人在线观看| 国产精品九九99| 美女高潮喷水抽搐中文字幕| 免费av毛片视频| 成人国产一区最新在线观看| 好男人电影高清在线观看| 欧美+亚洲+日韩+国产| 亚洲精品国产区一区二| 麻豆成人av在线观看| 色婷婷久久久亚洲欧美| 在线观看舔阴道视频| av天堂久久9| 国产片内射在线| 国产高清videossex| 国产高清videossex| 老司机福利观看| 最近最新中文字幕大全免费视频| 欧美成人性av电影在线观看| 国产精品久久久久久人妻精品电影| 日韩欧美免费精品| 亚洲精品一二三| 好男人电影高清在线观看| 神马国产精品三级电影在线观看 | 欧美最黄视频在线播放免费 | 国产精品免费视频内射| 亚洲avbb在线观看| 99精品在免费线老司机午夜| 精品久久久精品久久久| 黑人操中国人逼视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产一卡二卡三卡精品| 亚洲一区二区三区不卡视频| 欧美乱码精品一区二区三区| 国产精品1区2区在线观看.| 99国产精品99久久久久| www日本在线高清视频| 一区二区三区国产精品乱码| 久久人人爽av亚洲精品天堂| 老司机午夜福利在线观看视频| 日韩欧美在线二视频| www国产在线视频色| 中出人妻视频一区二区| 久久久国产一区二区| 丝袜人妻中文字幕| 日本免费a在线| 亚洲一区二区三区欧美精品| 中国美女看黄片| 国产亚洲精品久久久久5区| 久久久久九九精品影院| 女警被强在线播放| 久久影院123| 日日摸夜夜添夜夜添小说| 国产精品国产高清国产av| 国产aⅴ精品一区二区三区波| 欧美成人免费av一区二区三区| 精品久久久久久成人av| 亚洲男人的天堂狠狠| 黄色怎么调成土黄色| 国产日韩一区二区三区精品不卡| 久久青草综合色| 国产99久久九九免费精品| 中文字幕av电影在线播放| 人妻丰满熟妇av一区二区三区| www.熟女人妻精品国产| 欧美中文日本在线观看视频| 国产欧美日韩一区二区精品| 69精品国产乱码久久久| 久久香蕉激情| 在线观看免费午夜福利视频| 久久久久九九精品影院| 欧美成狂野欧美在线观看| 欧美成狂野欧美在线观看| 精品国产一区二区久久| 在线观看免费午夜福利视频| 免费在线观看亚洲国产| 性欧美人与动物交配| 国产高清videossex| 久久人妻av系列| 看免费av毛片| 国产精品乱码一区二三区的特点 | 免费在线观看黄色视频的| 亚洲国产欧美一区二区综合| 又大又爽又粗| 欧美日韩乱码在线| 亚洲aⅴ乱码一区二区在线播放 | 国产精品98久久久久久宅男小说| 亚洲专区国产一区二区| 黑人巨大精品欧美一区二区mp4| 亚洲一区二区三区色噜噜 | 久久久久久免费高清国产稀缺| 激情视频va一区二区三区| 亚洲欧美激情在线| 999久久久国产精品视频| 动漫黄色视频在线观看| 国产精品99久久99久久久不卡| 亚洲国产欧美日韩在线播放| 亚洲成人免费av在线播放| 91字幕亚洲| 老司机靠b影院| 亚洲国产欧美日韩在线播放| 精品国产美女av久久久久小说| 叶爱在线成人免费视频播放| 国产精品成人在线| 成熟少妇高潮喷水视频| 淫妇啪啪啪对白视频| 老司机亚洲免费影院| 精品电影一区二区在线| 亚洲精品在线观看二区| 一区二区日韩欧美中文字幕| 99在线人妻在线中文字幕| 久久人妻熟女aⅴ| 十八禁网站免费在线| 99re在线观看精品视频| 亚洲成人精品中文字幕电影 | 国产精品日韩av在线免费观看 | 女人精品久久久久毛片| www.www免费av| 欧美日韩福利视频一区二区| 亚洲免费av在线视频| 免费av中文字幕在线| 可以在线观看毛片的网站| 757午夜福利合集在线观看| 黄色 视频免费看| 一进一出好大好爽视频| 国产精品av久久久久免费| 国产欧美日韩综合在线一区二区| 中文字幕色久视频| 欧美久久黑人一区二区| 亚洲午夜精品一区,二区,三区| 自线自在国产av| 午夜福利欧美成人| 人人妻人人爽人人添夜夜欢视频| 最近最新中文字幕大全免费视频| x7x7x7水蜜桃| 成人影院久久| 国产熟女xx| 久久香蕉精品热| 黄色女人牲交| 久久久国产精品麻豆| 这个男人来自地球电影免费观看| 中文字幕色久视频| 国产99久久九九免费精品| 欧美最黄视频在线播放免费 | 一级黄色大片毛片| 国产成人系列免费观看| 91成年电影在线观看| 国产精品美女特级片免费视频播放器 | 18美女黄网站色大片免费观看| av超薄肉色丝袜交足视频| 中文字幕色久视频| 脱女人内裤的视频| 精品福利观看| 日韩欧美国产一区二区入口| 国产欧美日韩精品亚洲av| 高清毛片免费观看视频网站 | 欧美人与性动交α欧美软件| 神马国产精品三级电影在线观看 | 国产一区二区三区在线臀色熟女 | 亚洲精品美女久久久久99蜜臀| 欧美日韩黄片免| 乱人伦中国视频| 亚洲,欧美精品.| 亚洲国产欧美一区二区综合| 97超级碰碰碰精品色视频在线观看| 黑人猛操日本美女一级片| 一级毛片精品| 不卡一级毛片| 精品国产一区二区三区四区第35| 久久精品aⅴ一区二区三区四区| 国产亚洲欧美在线一区二区| 性色av乱码一区二区三区2| 电影成人av| 老司机午夜十八禁免费视频| 天堂影院成人在线观看| 男男h啪啪无遮挡| 夜夜看夜夜爽夜夜摸 | 精品久久久久久,| 国产亚洲欧美在线一区二区| 国产精品偷伦视频观看了| av片东京热男人的天堂| www国产在线视频色| 国产人伦9x9x在线观看| netflix在线观看网站| 999久久久精品免费观看国产| 精品久久久久久,| 一a级毛片在线观看| 亚洲中文av在线| 丝袜人妻中文字幕| 国产极品粉嫩免费观看在线| 欧美大码av| 91精品国产国语对白视频| 9热在线视频观看99| 天天躁夜夜躁狠狠躁躁| 亚洲视频免费观看视频| 满18在线观看网站| 色尼玛亚洲综合影院| 亚洲国产精品合色在线| 色在线成人网| www.999成人在线观看| 一个人观看的视频www高清免费观看 | 国产一卡二卡三卡精品| av超薄肉色丝袜交足视频| 午夜成年电影在线免费观看| 脱女人内裤的视频| 中文亚洲av片在线观看爽| 亚洲,欧美精品.| 午夜免费激情av| 淫妇啪啪啪对白视频| 精品一区二区三区视频在线观看免费 | 50天的宝宝边吃奶边哭怎么回事| 好看av亚洲va欧美ⅴa在| 女警被强在线播放| 亚洲一区高清亚洲精品| 午夜日韩欧美国产| 又紧又爽又黄一区二区| 免费看a级黄色片| 成年人免费黄色播放视频| 久久精品成人免费网站| 老司机午夜福利在线观看视频| 一区二区三区国产精品乱码| 国产视频一区二区在线看| 一a级毛片在线观看| 一区二区三区精品91| 两个人看的免费小视频| 亚洲专区国产一区二区| 999久久久国产精品视频| 婷婷精品国产亚洲av在线| 91精品三级在线观看| 久久久久国内视频| 久久久久精品国产欧美久久久| 国产精品一区二区精品视频观看| 久久精品国产清高在天天线| tocl精华| 亚洲一区二区三区不卡视频| 亚洲精品久久成人aⅴ小说| 男人的好看免费观看在线视频 | 两人在一起打扑克的视频| 精品久久久精品久久久| 亚洲欧美日韩无卡精品| 99国产精品一区二区蜜桃av| 中文字幕最新亚洲高清| 12—13女人毛片做爰片一| 久久久水蜜桃国产精品网| x7x7x7水蜜桃| 欧美国产精品va在线观看不卡| 80岁老熟妇乱子伦牲交| 久久精品亚洲av国产电影网| 人妻丰满熟妇av一区二区三区| 亚洲成人久久性| 国产精品自产拍在线观看55亚洲| 性少妇av在线| 精品国产国语对白av| 日韩精品免费视频一区二区三区| 久久九九热精品免费| 韩国av一区二区三区四区| 亚洲视频免费观看视频| 免费在线观看黄色视频的| 搡老熟女国产l中国老女人| 美女大奶头视频| 国产成人影院久久av| 看片在线看免费视频| 老司机靠b影院| 欧美激情高清一区二区三区| 性欧美人与动物交配| 一级a爱视频在线免费观看| 欧美日本中文国产一区发布| 热99re8久久精品国产| 日韩免费高清中文字幕av| 女人高潮潮喷娇喘18禁视频| 亚洲av电影在线进入| 精品午夜福利视频在线观看一区| 成年版毛片免费区| 99国产精品免费福利视频| 亚洲一卡2卡3卡4卡5卡精品中文| 最近最新中文字幕大全电影3 | 欧美午夜高清在线| bbb黄色大片| 97人妻天天添夜夜摸| www.999成人在线观看| 最新在线观看一区二区三区| av免费在线观看网站| 人人妻人人爽人人添夜夜欢视频| 久久久精品国产亚洲av高清涩受| www日本在线高清视频| 亚洲视频免费观看视频| 精品国产一区二区久久| 超碰成人久久| 无遮挡黄片免费观看| 日本免费a在线| 色综合欧美亚洲国产小说| 久久久久国产精品人妻aⅴ院| 国产av一区二区精品久久| 一级片免费观看大全| 最新美女视频免费是黄的| 热re99久久精品国产66热6| 欧美日韩中文字幕国产精品一区二区三区 | 欧美激情久久久久久爽电影 | 日韩成人在线观看一区二区三区| av在线天堂中文字幕 | 午夜日韩欧美国产| 国产三级在线视频| 久久精品91无色码中文字幕| 亚洲午夜精品一区,二区,三区| 国产精品影院久久| 狠狠狠狠99中文字幕| av天堂久久9| 91字幕亚洲| 国内毛片毛片毛片毛片毛片| 校园春色视频在线观看| 亚洲色图av天堂| 日本黄色日本黄色录像| aaaaa片日本免费| 亚洲一区中文字幕在线| 国产一区二区在线av高清观看| 国产成人欧美在线观看| 日本黄色视频三级网站网址| 天堂中文最新版在线下载| 国产高清videossex| 亚洲专区中文字幕在线| 少妇裸体淫交视频免费看高清 | 国产精品国产av在线观看| 国产av又大| 亚洲自偷自拍图片 自拍| 国产av一区二区精品久久| 国产极品粉嫩免费观看在线| 成年版毛片免费区| 亚洲精品粉嫩美女一区| 成人三级做爰电影| 一级黄色大片毛片| 日韩三级视频一区二区三区| 老熟妇乱子伦视频在线观看| 天堂影院成人在线观看| 午夜91福利影院| 欧美在线黄色| 一区二区三区精品91| 咕卡用的链子| 亚洲,欧美精品.| 色婷婷av一区二区三区视频| a级毛片黄视频| 欧美激情 高清一区二区三区| 免费在线观看亚洲国产| 色老头精品视频在线观看| www.999成人在线观看| 亚洲欧美激情综合另类| 亚洲av五月六月丁香网| a级片在线免费高清观看视频| 在线观看日韩欧美| 欧美国产精品va在线观看不卡| 国产精品免费一区二区三区在线| 三级毛片av免费| 变态另类成人亚洲欧美熟女 | 日本精品一区二区三区蜜桃| 色婷婷av一区二区三区视频| 国产精品免费视频内射| 亚洲精品国产精品久久久不卡| 成人国语在线视频| 中文亚洲av片在线观看爽| 国产精品自产拍在线观看55亚洲| 在线播放国产精品三级| 国产精品影院久久| 母亲3免费完整高清在线观看| 一级毛片女人18水好多| 最好的美女福利视频网| 久久天躁狠狠躁夜夜2o2o| 国产亚洲精品一区二区www| 免费少妇av软件| 麻豆国产av国片精品| 手机成人av网站| 亚洲av电影在线进入| 精品久久久久久久久久免费视频 | 99在线人妻在线中文字幕| 91麻豆av在线| 少妇粗大呻吟视频| 国产精品 国内视频| 真人做人爱边吃奶动态| av视频免费观看在线观看| 麻豆国产av国片精品| 国产精品国产高清国产av| 国产午夜精品久久久久久| 欧美一级毛片孕妇| 18禁观看日本| 电影成人av| 黑人操中国人逼视频| 国产精品久久久人人做人人爽| 人人妻人人澡人人看| 欧美精品亚洲一区二区| 国产精品自产拍在线观看55亚洲| 丰满人妻熟妇乱又伦精品不卡| www.自偷自拍.com| 在线观看午夜福利视频| 久久久久久久久中文| 免费av毛片视频| 中文欧美无线码| 亚洲自偷自拍图片 自拍| www.999成人在线观看| 亚洲精品国产一区二区精华液| 好男人电影高清在线观看| 丝袜人妻中文字幕| 好看av亚洲va欧美ⅴa在| 天堂√8在线中文| 欧美中文综合在线视频| 纯流量卡能插随身wifi吗| 狠狠狠狠99中文字幕| 日日爽夜夜爽网站| 欧美午夜高清在线| 老司机午夜十八禁免费视频| 久久香蕉精品热| 嫩草影院新地址| 丁香六月欧美| 久久中文看片网| 国产麻豆成人av免费视频| 国产精品98久久久久久宅男小说| 最近视频中文字幕2019在线8| 久久人妻av系列| 深夜a级毛片| 午夜激情福利司机影院| 久久精品91蜜桃| 嫩草影院新地址| 韩国av一区二区三区四区| 天堂动漫精品| 亚洲国产精品久久男人天堂| 亚洲乱码一区二区免费版| 99精品久久久久人妻精品| 麻豆成人av在线观看| 日本免费a在线| av专区在线播放| 五月玫瑰六月丁香| 99久久成人亚洲精品观看| 国产真实伦视频高清在线观看 | 欧美激情在线99| 99久久99久久久精品蜜桃| 男女床上黄色一级片免费看| 特大巨黑吊av在线直播| 最好的美女福利视频网| 禁无遮挡网站| 在线a可以看的网站| 欧美xxxx黑人xx丫x性爽| 午夜免费男女啪啪视频观看 | 丰满人妻熟妇乱又伦精品不卡| 亚洲自拍偷在线| 精品无人区乱码1区二区| 日韩中文字幕欧美一区二区| 毛片一级片免费看久久久久 | 狂野欧美白嫩少妇大欣赏| 亚洲中文字幕一区二区三区有码在线看| 美女xxoo啪啪120秒动态图 | 亚洲激情在线av| 男人狂女人下面高潮的视频| 又粗又爽又猛毛片免费看| 99精品久久久久人妻精品| 精品不卡国产一区二区三区| 亚洲黑人精品在线| 日本熟妇午夜| 国产伦人伦偷精品视频| 国产精品电影一区二区三区| 一二三四社区在线视频社区8| 久久精品国产亚洲av天美| 国产精品影院久久| 国产日本99.免费观看| 老熟妇仑乱视频hdxx| 日本 欧美在线| 激情在线观看视频在线高清| 国产高清有码在线观看视频| 国产伦一二天堂av在线观看| 制服丝袜大香蕉在线| 国产熟女xx| 亚洲无线在线观看| 亚洲人成网站在线播| 亚洲国产高清在线一区二区三| 亚洲五月婷婷丁香| 欧美不卡视频在线免费观看| 毛片一级片免费看久久久久 | 性欧美人与动物交配| 国产野战对白在线观看| 久99久视频精品免费| 国产午夜福利久久久久久| 午夜精品在线福利| 动漫黄色视频在线观看| 国产精品嫩草影院av在线观看 | 亚洲第一欧美日韩一区二区三区| av中文乱码字幕在线| 亚洲成人久久爱视频| 成人毛片a级毛片在线播放| 最近在线观看免费完整版| 亚洲成人免费电影在线观看| 国产伦人伦偷精品视频| 特级一级黄色大片| 亚洲性夜色夜夜综合| 99久久99久久久精品蜜桃| 女生性感内裤真人,穿戴方法视频| 亚洲欧美日韩高清在线视频| 国产精品一区二区三区四区久久| 国产高清三级在线| 国产免费av片在线观看野外av| 亚洲不卡免费看| 少妇丰满av| 亚洲不卡免费看| 国产伦精品一区二区三区视频9| 国产精品亚洲av一区麻豆| 免费看a级黄色片| 又爽又黄无遮挡网站| 亚洲国产精品sss在线观看| 91麻豆av在线| 亚洲国产精品sss在线观看| 男人狂女人下面高潮的视频| 成人三级黄色视频| 亚洲av不卡在线观看| 搡女人真爽免费视频火全软件 | 欧美激情久久久久久爽电影| 69av精品久久久久久|