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

    基于EEMD和模糊BP神經(jīng)網(wǎng)絡(luò)的滾動(dòng)軸承故障診斷

    2017-11-02 04:07:07張國銀王海瑞郝家驥宋怡然
    化工自動(dòng)化及儀表 2017年1期
    關(guān)鍵詞:故障診斷模態(tài)振動(dòng)

    張國銀 王 雪 王海瑞 郝家驥 宋怡然

    (昆明理工大學(xué)信息工程與自動(dòng)化學(xué)院)

    基于EEMD和模糊BP神經(jīng)網(wǎng)絡(luò)的滾動(dòng)軸承故障診斷

    張國銀 王 雪 王海瑞 郝家驥 宋怡然

    (昆明理工大學(xué)信息工程與自動(dòng)化學(xué)院)

    針對滾動(dòng)軸承故障振動(dòng)信號的非平穩(wěn)特征,提出一種基于總體經(jīng)驗(yàn)?zāi)B(tài)分解(EEMD)和模糊BP神經(jīng)網(wǎng)絡(luò)的故障診斷方法。首先對滾動(dòng)軸承的振動(dòng)信號采用總體經(jīng)驗(yàn)?zāi)B(tài)分解方法進(jìn)行分解,得到若干個(gè)本征模態(tài)函數(shù)分量(IMF);然后提取各分量的均方差、峭度和能量,把這些特征參數(shù)作為學(xué)習(xí)集和訓(xùn)練集,將學(xué)習(xí)集輸入到模糊BP神經(jīng)網(wǎng)絡(luò)中進(jìn)行學(xué)習(xí);最后把訓(xùn)練集輸入到特征參數(shù)經(jīng)過學(xué)習(xí)訓(xùn)練后的模糊BP神經(jīng)網(wǎng)絡(luò)中進(jìn)行故障類型識別,并與BP神經(jīng)網(wǎng)絡(luò)進(jìn)行比較。實(shí)驗(yàn)結(jié)果表明:所提方法能有效地應(yīng)用于滾動(dòng)軸承故障診斷,而且比BP神經(jīng)網(wǎng)絡(luò)具有更高的精確度。

    滾動(dòng)軸承 故障診斷 EEMD 模糊BP神經(jīng)網(wǎng)絡(luò) 能量 峭度

    目前,在滾動(dòng)軸承的故障診斷中,主要以研究振動(dòng)信號中所產(chǎn)生的故障特征進(jìn)行故障診斷。當(dāng)軸承中的一個(gè)或幾個(gè)元件發(fā)生局部故障時(shí),相接觸的部件會(huì)因撞擊而產(chǎn)生能量集中的周期性脈沖信號,從而激起軸承系統(tǒng)各部件的高頻固有振動(dòng),產(chǎn)生幅值調(diào)制現(xiàn)象,不同的故障會(huì)有不同的故障信息特征[1]。一般來說,滾動(dòng)軸承故障診斷技術(shù)分為4個(gè)過程:首先對選取的能反映滾動(dòng)軸承狀態(tài)信息的信號進(jìn)行預(yù)處理;接著對信號進(jìn)行特征提取,選取能夠反映其故障征兆的特征值;然后根據(jù)特征值選取恰當(dāng)?shù)奶幚矸椒▽μ卣髦颠M(jìn)行分析,從而得出其故障狀態(tài)診斷與分析;最后按照結(jié)果做出相應(yīng)的決策,及時(shí)對設(shè)備進(jìn)行相應(yīng)的處理。而提取故障特征的有效性決定了能否準(zhǔn)確地判斷故障類型。因此,判斷滾動(dòng)軸承故障的關(guān)鍵就是準(zhǔn)確全面地提取故障特征信號。

    文獻(xiàn)[2]應(yīng)用模糊神經(jīng)網(wǎng)絡(luò)對軌道電路故障進(jìn)行診斷,取得了較好的效果,但沒有對原始信號進(jìn)行預(yù)處理,原始信號中摻雜了干擾信號,會(huì)影響對故障類型的判斷。文獻(xiàn)[3]應(yīng)用經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)和神經(jīng)網(wǎng)絡(luò)進(jìn)行輪軌故障噪聲診斷,但是EMD在分解過程中會(huì)產(chǎn)生模態(tài)混疊現(xiàn)象,同樣會(huì)對故障分類的識別產(chǎn)生影響。EMD可用于非線性與非平穩(wěn)信號處理,可是EMD分解本身存在模態(tài)混疊等不足,為了減少模態(tài)混疊對實(shí)驗(yàn)結(jié)果帶來的影響,Hang N E和Wu Z H提出利用高斯白噪聲的總體經(jīng)驗(yàn)?zāi)B(tài)分解(EEMD)方法[4,5]。模糊神經(jīng)網(wǎng)絡(luò)的關(guān)鍵是對于模糊規(guī)則及其隸屬度函數(shù)的改進(jìn),是基于神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)能力實(shí)現(xiàn),此方法將模糊控制與神經(jīng)網(wǎng)絡(luò)聯(lián)系在一起,既改進(jìn)了不足,又融合了它的優(yōu)點(diǎn)[6]。基于上述研究,筆者將EEMD方法和模糊BP神經(jīng)網(wǎng)絡(luò)結(jié)合起來,對滾動(dòng)軸承信號進(jìn)行分析與故障診斷,首先對采集到的軸承噪聲信號進(jìn)行EEMD分解,計(jì)算各IMF分量的均方差值、峭度和能量,并將其各均值作為模糊BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練和測試的特征參數(shù)輸入,對滾動(dòng)軸承3種不同狀況下的故障進(jìn)行分類識別。

    1 EEMD方法

    滾動(dòng)軸承振動(dòng)信號中的噪聲干擾嚴(yán)重影響故障診斷的準(zhǔn)確性,為此需要使用消噪方法對信號進(jìn)行處理以減少噪聲干擾。具體的消噪方法包括小波濾噪、EMD分解濾噪等。小波消噪具有多分辨率等優(yōu)點(diǎn),但小波參數(shù)的選擇將直接影響診斷結(jié)果,并且不能自適應(yīng)地進(jìn)行信號分解。EMD作為新的消噪方法,具有自適應(yīng)分解信號的能力,但EMD分解存在模態(tài)混疊、端點(diǎn)效應(yīng)等問題。因此Huang N E等提出了EEMD分解方法[4,5]。EEMD的本質(zhì)是在原有信號中加入高斯白噪聲的多次EMD,利用白噪聲的頻率平均分布的統(tǒng)計(jì)特性,自適應(yīng)地從滾動(dòng)軸承振動(dòng)信號中分離出高頻調(diào)制信號,從而有效抑制模態(tài)混疊問題。EEMD方法既克服了小波基函數(shù)和閾值的選擇,又解決了EMD的模態(tài)混疊問題。

    EEMD算法分解過程的具體步驟如下:

    a. 對分析信號x(t)加入隨機(jī)高斯白噪聲序列(一般白噪聲標(biāo)準(zhǔn)差是原始信號標(biāo)準(zhǔn)差的0.1~0.4倍);

    b. 對加入白噪聲的信號進(jìn)行EMD分解,得到K個(gè)IMF分量cj(t)(j=1,…,K)和余項(xiàng)r(t);

    c. 計(jì)算分解后IMF的均值,把各個(gè)IMF的均值作為最終的結(jié)果。

    將步驟a~b重復(fù)N次,每次加入不同的白噪聲序列,所得的第i次加入白噪聲后的信號xi(t)為:

    (1)

    式中cij(t)——第i次加入高斯白噪聲后分解所得的第j個(gè)IMF。

    步驟c對原信號經(jīng)過EEMD分解后所得的第j個(gè)IMF分量的最終結(jié)果cj(t)為:

    (2)

    (3)

    2 模糊BP神經(jīng)網(wǎng)絡(luò)

    模糊BP神經(jīng)網(wǎng)絡(luò)整合了模糊邏輯系統(tǒng)和神經(jīng)網(wǎng)絡(luò)的優(yōu)點(diǎn),模糊邏輯系統(tǒng)善于通過表達(dá)和抽取知識來處理結(jié)構(gòu)化的知識,但是不具備學(xué)習(xí)能力;神經(jīng)網(wǎng)絡(luò)則善于通過樣本學(xué)習(xí)來處理非結(jié)構(gòu)的知識,但是不能對模糊相似信息進(jìn)行有效區(qū)分。這樣的結(jié)合方式使模糊BP神經(jīng)網(wǎng)絡(luò)不僅具有知識易于理解的優(yōu)點(diǎn),還無需建立精確的數(shù)學(xué)模型,適合處理不確定性和非線性問題[7]。模糊BP神經(jīng)網(wǎng)絡(luò)的兩個(gè)組成部分為前件網(wǎng)絡(luò)和后件網(wǎng)絡(luò),前件網(wǎng)絡(luò)的作用為匹配模糊規(guī)則,后件網(wǎng)絡(luò)的作用則是產(chǎn)生模糊規(guī)則[8]。筆者采用BP算法對各節(jié)點(diǎn)的隸屬度函數(shù)進(jìn)行優(yōu)化,具體的網(wǎng)絡(luò)結(jié)構(gòu)有5層,如圖1所示。

    圖1 模糊BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)

    第1層為輸入層,有3個(gè)神經(jīng)元(均方差、峭度和能量)作為輸入信號(x1,x2,x3)傳送到下一層。

    3 基于EEMD與模糊BP神經(jīng)網(wǎng)絡(luò)的滾動(dòng)軸承故障診斷方法

    峭度值反映了信號中沖擊成分的大小[9],沖擊成分所占的比率越大,峭度值越大。由于高頻IMF信號中含有故障引起的固有振動(dòng)成分,因此故障信息最容易提取。當(dāng)軸承出現(xiàn)不同的故障時(shí),故障信號在各個(gè)IMF分量中的能量分布會(huì)發(fā)生相應(yīng)的改變,所以可以選擇不同頻帶的能量作為故障識別的特征參量。均方差能反映信號幅值概率密度分布的特性,即幅值分布的離散程度,當(dāng)滾動(dòng)軸承出現(xiàn)故障時(shí),由于沖擊的增強(qiáng)將導(dǎo)致信號的幅值分布特性發(fā)生變化。根據(jù)上文所述,可將各IMF分量的均方差值、峭度值和能量值作為模糊BP神經(jīng)網(wǎng)絡(luò)的特征參數(shù),能夠進(jìn)行較好的故障分類識別。

    通過上述分析,利用EEMD與模糊BP神經(jīng)網(wǎng)絡(luò)的滾動(dòng)軸承故障診斷方法的具體步驟為:

    a. 對內(nèi)圈故障、外圈故障和滾動(dòng)體故障狀態(tài)下的信號進(jìn)行采集;

    b. 對每一種狀態(tài)下的信號進(jìn)行EEMD分解,得到若干個(gè)IMF分量,按其頻率大小從高到低依次排列,由于故障信息主要存在于信號的高頻部分,故選取前6個(gè)IMF分量進(jìn)行分析;

    d. 計(jì)算每一個(gè)信號中分解的6個(gè)IMF分量的峭度均值、能量均值和均方差均值,把這些均值作為模糊BP神經(jīng)網(wǎng)絡(luò)的參數(shù)輸入;

    e. 將得到的均方差、峭度和能量的均值中的30組數(shù)據(jù)作為訓(xùn)練樣本對模糊BP神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練,訓(xùn)練結(jié)束后再用其他21組數(shù)據(jù)作為測試樣本進(jìn)行故障識別。

    4 實(shí)驗(yàn)結(jié)果

    為了證明筆者所提方法的可行性,用美國凱斯西儲大學(xué)振動(dòng)實(shí)驗(yàn)室提供的數(shù)據(jù)進(jìn)行實(shí)驗(yàn)。實(shí)驗(yàn)臺包括一個(gè)電動(dòng)機(jī)、一個(gè)扭矩傳感器、一個(gè)功率測試計(jì)和電子控制器。使用加速度傳感器采集軸承的振動(dòng)加速度信號,分別放置在風(fēng)扇端軸承座和驅(qū)動(dòng)端的軸承座上。采集此振動(dòng)信號的采樣頻率為12kHz。實(shí)驗(yàn)中,使用驅(qū)動(dòng)端的加速度數(shù)據(jù)進(jìn)行滾動(dòng)軸承的故障特征提取,進(jìn)而判斷出故障類型。

    由于信號中的特征信息主要集中在高頻部分,而低頻IMF分量所含的特征信息過少,所以保留EEMD分解后的IMF1~I(xiàn)MF6。采用的滾動(dòng)軸承故障信號包含內(nèi)圈故障、外圈故障和滾動(dòng)體故障3種狀態(tài)。圖2為不同狀態(tài)下的EEMD分解結(jié)果。

    圖2 不同狀態(tài)下的EEMD分解結(jié)果

    圖3中,模糊BP神經(jīng)網(wǎng)絡(luò)在經(jīng)過225次訓(xùn)練后,誤差精度為0.001 0,網(wǎng)絡(luò)訓(xùn)練結(jié)束。圖4中,BP神經(jīng)網(wǎng)絡(luò)在經(jīng)過358次訓(xùn)練后,誤差精度為0.022 3,無法達(dá)到期望的精度,網(wǎng)絡(luò)訓(xùn)練結(jié)束。兩者相比,模糊BP神經(jīng)網(wǎng)絡(luò)不但在訓(xùn)練次數(shù)上比BP神經(jīng)網(wǎng)絡(luò)少,而且訓(xùn)練誤差也達(dá)到了期望值,訓(xùn)練精度也有所提高。

    圖3 模糊BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練誤差變化曲線

    圖4 BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練誤差變化曲線

    用訓(xùn)練好的模糊BP神經(jīng)網(wǎng)絡(luò)和BP神經(jīng)網(wǎng)絡(luò)分別對測試樣本進(jìn)行分類識別,每種情況各7個(gè)測試樣本,一共21組數(shù)據(jù),見表1。

    為了更好地證明筆者所提方法的可行性,表2中列出了兩種方法的實(shí)際輸出結(jié)果,并與期望輸出結(jié)果和故障類型進(jìn)行對比,BP神經(jīng)網(wǎng)絡(luò)的輸出設(shè)定為3個(gè)值,分別代表3個(gè)故障,值最大的一個(gè)就對應(yīng)相應(yīng)的故障。表2中BP神經(jīng)網(wǎng)絡(luò)的輸出樣本“1”代表滾動(dòng)軸承的內(nèi)圈故障,“2”代表滾動(dòng)軸承的外圈故障,“3”代表滾動(dòng)軸承處于滾動(dòng)體故障狀態(tài)。同時(shí),為了更加直觀地判斷兩種方法的優(yōu)劣,在表3中列出了兩種方法的故障識別率對比。

    表1 測試數(shù)據(jù)

    表2 模糊BP神經(jīng)網(wǎng)絡(luò)與BP神經(jīng)網(wǎng)絡(luò)診斷結(jié)果

    (續(xù)表2)

    表3 模糊BP神經(jīng)網(wǎng)絡(luò)與BP神經(jīng)網(wǎng)絡(luò)故障識別率對比

    表2中,故障信號經(jīng)過EEMD分解后再使用上述兩種方法都可以根據(jù)故障信號中的信息進(jìn)行故障類型的識別,根據(jù)對比結(jié)果可知:模糊BP神經(jīng)網(wǎng)絡(luò)比BP神經(jīng)網(wǎng)絡(luò)具有更高的故障識別率(表3)。在此次實(shí)驗(yàn)過程中,模糊BP神經(jīng)網(wǎng)絡(luò)出現(xiàn)了一個(gè)識別錯(cuò)誤,準(zhǔn)確率達(dá)到95.2%;但BP神經(jīng)網(wǎng)絡(luò)中出現(xiàn)了3個(gè)識別錯(cuò)誤,準(zhǔn)確率只有85.7%。而且BP神經(jīng)網(wǎng)絡(luò)的錯(cuò)誤故障識別都是出現(xiàn)在內(nèi)圈故障中,這是因?yàn)锽P神經(jīng)網(wǎng)絡(luò)不能處理和描述模糊信息,當(dāng)數(shù)據(jù)過于接近時(shí)就無法正確區(qū)分?jǐn)?shù)據(jù)而導(dǎo)致產(chǎn)生錯(cuò)誤的結(jié)果;而模糊BP神經(jīng)網(wǎng)絡(luò)不但具有模糊處理能力,還具有學(xué)習(xí)能力,彌補(bǔ)了BP神經(jīng)網(wǎng)絡(luò)的不足。

    5 結(jié)束語

    針對滾動(dòng)軸承的故障分析,筆者提出一種將EEMD和模糊BP神經(jīng)網(wǎng)絡(luò)結(jié)合的方法,實(shí)驗(yàn)結(jié)果表明該方法不僅對噪聲產(chǎn)生了一定的抑制效果,還可以實(shí)現(xiàn)多種故障的分類,達(dá)到了較好的精度,表明該方法可以有效地區(qū)分滾動(dòng)軸承故障的類別。但是在此次實(shí)驗(yàn)中,模糊BP神經(jīng)網(wǎng)絡(luò)仍然出現(xiàn)了識別錯(cuò)誤,這是因?yàn)槟:鼴P神經(jīng)網(wǎng)絡(luò)的模糊規(guī)則具有冗余和不精簡性,當(dāng)輸入?yún)?shù)和對應(yīng)的模糊子集較多時(shí),神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)將會(huì)變得比較復(fù)雜,影響計(jì)算效率。因此,改進(jìn)模糊BP神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)和參數(shù)將是下一步研究的重點(diǎn)。

    [1] 胡愛軍,馬萬里,唐貴基.基于集成經(jīng)驗(yàn)?zāi)B(tài)分解和峭度準(zhǔn)則的滾動(dòng)軸承故障特征提取方法[J].中國電機(jī)工程學(xué)報(bào),2012,32(11):106~111.

    [2] 黃贊武,魏雪業(yè),劉澤.基于模糊神經(jīng)網(wǎng)絡(luò)的軌道電路故障診斷方法研究[J].鐵道學(xué)報(bào),2012,34(11):54~59.

    [3] 江航,尚春陽,高瑞鵬.基于EMD和神經(jīng)網(wǎng)絡(luò)的輪軌故障噪聲診斷識別方法研究[J].振動(dòng)與沖擊,2014,33(17):34~38.

    [4] Huang N E,Shen Z,Long S R,et al.The Empirical Mode of Decomposition and the Hilbert Spectrum for Nonlinear and Non-stationary Time Series Analysis[J].Proc Royal Society,1988,454:903~905.

    [5] Wu Z H,Huang N E.Ensemble Empirical Mode Decomposition:A Noise-Assisted Data Analysis Method[J].Advances in Adaptive Data Analysis,2009,1(1):1~4.

    [6] 彭廣平,李泳鮮.模糊神經(jīng)網(wǎng)絡(luò)在汽車發(fā)動(dòng)機(jī)電子點(diǎn)火系統(tǒng)故障診斷中的應(yīng)用[J].汽車科技,2011,(3):47~50.

    [7] 李恒嵬.模糊神經(jīng)網(wǎng)絡(luò)研究現(xiàn)狀綜述[J].遼寧科技學(xué)院學(xué)報(bào),2010,12(2):15~17.

    [8] 周潤景.基于MATLAB與fuzzyTECH的模糊與神經(jīng)網(wǎng)絡(luò)設(shè)計(jì)[M].北京:電子工業(yè)出版社,2010:197~199.

    [9] 吳小濤,楊錳,袁曉輝,等.基于峭度準(zhǔn)則EEMD及改進(jìn)形態(tài)濾波方法的軸承故障診斷[J].振動(dòng)與沖擊,2015,34(2):38~44.

    RollingBearingFaultDiagnosisBasedonEEMDandFuzzyBPNeuralNetwork

    ZHANG Guo-yin, WANG Xue, WANG Hai-rui, HAO Jia-ji, SONG Yi-ran

    (FacultyofInformationEngineeringandAutomation,KunmingUniversityofScienceandTechnology)

    Considering non-stationary characteristics of bearing’s fault vibration signals, an EEMD and fuzzy neural network-based fault diagnosis method was proposed, in which, having EEMD method adopted to decompose the bearing’s bearing vibration signals into several intrinsic mode function components(IMF); and then, extracting each component’s mean square error, kurtosis and energy and taking them as learning set and training set; and finally, having the learning set input into the fuzzy BP neural network for learning and the training set into the fuzzy BP neural network for fault type identification. Experimental results show that, the method proposed can be effectively applied to the rolling bearing’s fault diagnosis and it outperforms the BP neural network in the accuracy.

    rolling bearing, fault diagnosis, EEMD, fuzzy BP neural network, energy, kurtosis

    TH133.33

    A

    1000-3932(2017)01-0034-06

    國家自然科學(xué)基金項(xiàng)目(61263023)。

    張國銀(1979-),副教授,從事計(jì)算機(jī)應(yīng)用和智能算法方面的研究。

    聯(lián) 系 人:王雪(1990-),碩士研究生,從事模式識別和自動(dòng)化技術(shù)方面的研究,1067834094@qq.com。

    2016-03-23,

    2016-10-10)

    猜你喜歡
    故障診斷模態(tài)振動(dòng)
    振動(dòng)的思考
    振動(dòng)與頻率
    中立型Emden-Fowler微分方程的振動(dòng)性
    國內(nèi)多模態(tài)教學(xué)研究回顧與展望
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識別
    UF6振動(dòng)激發(fā)態(tài)分子的振動(dòng)-振動(dòng)馳豫
    由單個(gè)模態(tài)構(gòu)造對稱簡支梁的抗彎剛度
    基于LCD和排列熵的滾動(dòng)軸承故障診斷
    基于WPD-HHT的滾動(dòng)軸承故障診斷
    bbb黄色大片| 熟女少妇亚洲综合色aaa.| 69精品国产乱码久久久| 成人黄色视频免费在线看| 91国产中文字幕| 久久亚洲真实| 男女之事视频高清在线观看| 国产男女超爽视频在线观看| 久久精品国产亚洲av高清一级| 亚洲一码二码三码区别大吗| 国产成人影院久久av| 国产黄色免费在线视频| 男人操女人黄网站| 久久性视频一级片| 美女主播在线视频| 夜夜夜夜夜久久久久| 十分钟在线观看高清视频www| 欧美日韩亚洲国产一区二区在线观看 | 久久国产精品人妻蜜桃| 黄色成人免费大全| 人妻 亚洲 视频| 亚洲人成77777在线视频| 亚洲熟女精品中文字幕| xxxhd国产人妻xxx| 国产成人欧美| 欧美日本中文国产一区发布| 欧美人与性动交α欧美精品济南到| 伦理电影免费视频| 中文亚洲av片在线观看爽 | 人人妻人人澡人人爽人人夜夜| 成人手机av| 狂野欧美激情性xxxx| 美女福利国产在线| av天堂在线播放| 亚洲熟妇熟女久久| 三级毛片av免费| 久久影院123| 国产精品一区二区精品视频观看| 久久精品aⅴ一区二区三区四区| 亚洲人成77777在线视频| 国产色视频综合| 成人18禁高潮啪啪吃奶动态图| 亚洲一区中文字幕在线| 久热这里只有精品99| 亚洲av日韩在线播放| 十分钟在线观看高清视频www| 成人精品一区二区免费| 国产欧美日韩一区二区三区在线| 日韩熟女老妇一区二区性免费视频| 男人操女人黄网站| 日日爽夜夜爽网站| av线在线观看网站| 五月开心婷婷网| 香蕉国产在线看| 欧美大码av| 97人妻天天添夜夜摸| 天天影视国产精品| 啪啪无遮挡十八禁网站| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩av久久| 日韩成人在线观看一区二区三区| 狠狠婷婷综合久久久久久88av| 国产欧美日韩一区二区精品| 国产1区2区3区精品| 少妇粗大呻吟视频| 视频区欧美日本亚洲| 69精品国产乱码久久久| 精品高清国产在线一区| 国产99久久九九免费精品| 麻豆乱淫一区二区| 一本一本久久a久久精品综合妖精| 中亚洲国语对白在线视频| 在线av久久热| 岛国在线观看网站| 男女无遮挡免费网站观看| 亚洲第一青青草原| 日韩精品免费视频一区二区三区| 色尼玛亚洲综合影院| 中文字幕人妻熟女乱码| 免费人妻精品一区二区三区视频| 香蕉国产在线看| 国产成人免费观看mmmm| 亚洲综合色网址| 夫妻午夜视频| 99精品在免费线老司机午夜| 欧美性长视频在线观看| 国产精品99久久99久久久不卡| 少妇粗大呻吟视频| 午夜福利欧美成人| 亚洲av美国av| 制服诱惑二区| 日日爽夜夜爽网站| 极品少妇高潮喷水抽搐| 极品教师在线免费播放| 午夜激情久久久久久久| 欧美精品亚洲一区二区| 十八禁网站网址无遮挡| bbb黄色大片| 国产成人精品久久二区二区免费| 日韩视频在线欧美| 最近最新中文字幕大全免费视频| 亚洲国产欧美在线一区| 曰老女人黄片| 日韩欧美一区视频在线观看| 欧美日韩成人在线一区二区| 国产精品电影一区二区三区 | 91麻豆精品激情在线观看国产 | 久久精品亚洲av国产电影网| 亚洲第一青青草原| 成人18禁在线播放| 欧美黄色淫秽网站| 熟女少妇亚洲综合色aaa.| 一级片'在线观看视频| 国产人伦9x9x在线观看| 欧美日韩国产mv在线观看视频| 久热爱精品视频在线9| 嫩草影视91久久| 一区二区日韩欧美中文字幕| 欧美激情极品国产一区二区三区| 亚洲av电影在线进入| 免费看十八禁软件| 精品国产乱码久久久久久小说| 高清在线国产一区| 无限看片的www在线观看| 亚洲中文字幕日韩| videosex国产| 美女扒开内裤让男人捅视频| 免费在线观看影片大全网站| 免费高清在线观看日韩| 男女高潮啪啪啪动态图| 精品一区二区三区视频在线观看免费 | 深夜精品福利| 99久久精品国产亚洲精品| 又黄又粗又硬又大视频| 1024香蕉在线观看| 99久久精品国产亚洲精品| 中文字幕精品免费在线观看视频| 亚洲第一青青草原| 日韩一卡2卡3卡4卡2021年| 免费日韩欧美在线观看| 国产欧美日韩精品亚洲av| 99久久国产精品久久久| av不卡在线播放| 国产精品香港三级国产av潘金莲| 欧美日韩一级在线毛片| 91九色精品人成在线观看| 精品一区二区三卡| 亚洲天堂av无毛| 欧美乱码精品一区二区三区| 久久久精品免费免费高清| 制服诱惑二区| 色婷婷久久久亚洲欧美| 日韩欧美三级三区| 十分钟在线观看高清视频www| 曰老女人黄片| 人人妻人人澡人人爽人人夜夜| av又黄又爽大尺度在线免费看| 水蜜桃什么品种好| 国产一区二区三区综合在线观看| 欧美黄色淫秽网站| 热re99久久精品国产66热6| 精品久久久久久久毛片微露脸| 国产区一区二久久| 麻豆成人av在线观看| 啦啦啦 在线观看视频| 亚洲成a人片在线一区二区| 亚洲色图av天堂| 国产av一区二区精品久久| 国产欧美日韩一区二区精品| 亚洲午夜精品一区,二区,三区| avwww免费| 后天国语完整版免费观看| 最近最新中文字幕大全电影3 | 老汉色∧v一级毛片| 免费高清在线观看日韩| 两个人看的免费小视频| 狠狠精品人妻久久久久久综合| 欧美日韩亚洲综合一区二区三区_| 熟女少妇亚洲综合色aaa.| 在线观看66精品国产| 亚洲国产看品久久| bbb黄色大片| 高清av免费在线| 亚洲av国产av综合av卡| 亚洲欧洲精品一区二区精品久久久| 欧美成人午夜精品| 女人久久www免费人成看片| 成人三级做爰电影| 老司机靠b影院| 高清视频免费观看一区二区| 欧美激情 高清一区二区三区| 一区二区三区激情视频| 日韩大片免费观看网站| 99热网站在线观看| 国产亚洲午夜精品一区二区久久| 99久久99久久久精品蜜桃| 好男人电影高清在线观看| 在线观看免费视频网站a站| 国产在线免费精品| 亚洲av国产av综合av卡| 18禁黄网站禁片午夜丰满| 免费不卡黄色视频| 欧美大码av| 欧美 亚洲 国产 日韩一| 男人舔女人的私密视频| av有码第一页| 国产成人av教育| netflix在线观看网站| 亚洲三区欧美一区| 18禁国产床啪视频网站| 久久精品亚洲av国产电影网| 国产免费视频播放在线视频| 国产一区有黄有色的免费视频| 一夜夜www| av天堂久久9| 热re99久久精品国产66热6| 欧美日本中文国产一区发布| 久久久久视频综合| 精品福利观看| 国产在线视频一区二区| 亚洲五月婷婷丁香| 中文字幕av电影在线播放| 伊人久久大香线蕉亚洲五| 久久人妻福利社区极品人妻图片| 精品少妇久久久久久888优播| 亚洲 国产 在线| 叶爱在线成人免费视频播放| 亚洲成国产人片在线观看| 午夜激情av网站| 欧美日韩中文字幕国产精品一区二区三区 | 在线天堂中文资源库| 日韩中文字幕视频在线看片| 热re99久久国产66热| 91字幕亚洲| 欧美激情 高清一区二区三区| 一区二区三区国产精品乱码| 中文字幕精品免费在线观看视频| 亚洲国产成人一精品久久久| 日本黄色日本黄色录像| 久久久国产成人免费| 免费在线观看日本一区| 欧美成人免费av一区二区三区 | 成人国产av品久久久| 2018国产大陆天天弄谢| 色婷婷av一区二区三区视频| 丝袜人妻中文字幕| 国产精品影院久久| 久久精品国产a三级三级三级| 国产精品电影一区二区三区 | 人妻 亚洲 视频| 久久av网站| 狠狠狠狠99中文字幕| 一区二区三区精品91| 老鸭窝网址在线观看| 亚洲熟女精品中文字幕| 欧美一级毛片孕妇| 久久国产精品人妻蜜桃| 夫妻午夜视频| 老鸭窝网址在线观看| 国产av精品麻豆| 99精品欧美一区二区三区四区| 国产精品久久久人人做人人爽| 19禁男女啪啪无遮挡网站| 久久久久久久精品吃奶| 老司机午夜十八禁免费视频| 看免费av毛片| 久久毛片免费看一区二区三区| 久久久久国产一级毛片高清牌| 亚洲人成电影观看| 成人手机av| 国产色视频综合| 美女高潮喷水抽搐中文字幕| www.自偷自拍.com| 久久国产精品男人的天堂亚洲| 十八禁网站免费在线| 19禁男女啪啪无遮挡网站| 亚洲人成伊人成综合网2020| 看免费av毛片| 亚洲成人免费av在线播放| 男女午夜视频在线观看| 电影成人av| www.熟女人妻精品国产| 亚洲国产成人一精品久久久| 中文字幕人妻丝袜制服| 精品午夜福利视频在线观看一区 | 亚洲中文字幕日韩| 十八禁网站免费在线| 欧美黄色淫秽网站| 国产精品亚洲一级av第二区| 免费看a级黄色片| 精品一品国产午夜福利视频| 男女下面插进去视频免费观看| 欧美中文综合在线视频| 91字幕亚洲| 天天添夜夜摸| 丰满人妻熟妇乱又伦精品不卡| 久久久久久久国产电影| 久久久久久久精品吃奶| 色播在线永久视频| 91麻豆精品激情在线观看国产 | 97在线人人人人妻| 交换朋友夫妻互换小说| 一区二区三区国产精品乱码| 国产av一区二区精品久久| 麻豆国产av国片精品| 亚洲精品久久午夜乱码| 男女午夜视频在线观看| 日韩大码丰满熟妇| 美国免费a级毛片| 成年动漫av网址| tocl精华| 超碰成人久久| 亚洲精品一二三| 后天国语完整版免费观看| 美国免费a级毛片| 国产精品免费一区二区三区在线 | 啦啦啦 在线观看视频| 成人国产av品久久久| 91麻豆精品激情在线观看国产 | 日韩免费高清中文字幕av| 久久精品国产亚洲av香蕉五月 | 老司机午夜十八禁免费视频| 午夜激情av网站| 美女高潮喷水抽搐中文字幕| 大码成人一级视频| 丝袜在线中文字幕| 国产福利在线免费观看视频| 日日爽夜夜爽网站| 啦啦啦免费观看视频1| 亚洲五月婷婷丁香| 久久精品国产综合久久久| h视频一区二区三区| 一夜夜www| 成人三级做爰电影| www.自偷自拍.com| 十八禁人妻一区二区| 精品卡一卡二卡四卡免费| 男女无遮挡免费网站观看| 老汉色∧v一级毛片| 性高湖久久久久久久久免费观看| 国产亚洲午夜精品一区二区久久| av天堂在线播放| 一级片免费观看大全| 免费在线观看视频国产中文字幕亚洲| 老汉色av国产亚洲站长工具| 日韩人妻精品一区2区三区| 丝袜美足系列| 黄片小视频在线播放| 免费人妻精品一区二区三区视频| 国产av国产精品国产| 精品国产一区二区三区四区第35| 视频区图区小说| 超色免费av| 午夜日韩欧美国产| 国产欧美日韩一区二区三| 国产精品亚洲av一区麻豆| 亚洲专区字幕在线| 一区二区日韩欧美中文字幕| 精品一区二区三区四区五区乱码| 麻豆成人av在线观看| 亚洲九九香蕉| 日本av手机在线免费观看| 搡老岳熟女国产| 91av网站免费观看| 人人妻人人澡人人看| 午夜91福利影院| 午夜福利乱码中文字幕| 国产欧美日韩综合在线一区二区| 国产精品99久久99久久久不卡| 90打野战视频偷拍视频| 成人18禁高潮啪啪吃奶动态图| 97在线人人人人妻| 久久香蕉激情| 91九色精品人成在线观看| 日本精品一区二区三区蜜桃| 亚洲精品国产一区二区精华液| 激情在线观看视频在线高清 | 啦啦啦 在线观看视频| 午夜老司机福利片| 又黄又粗又硬又大视频| 啦啦啦 在线观看视频| 日韩三级视频一区二区三区| 国产在线免费精品| 国产色视频综合| 国产精品电影一区二区三区 | 免费不卡黄色视频| 一边摸一边抽搐一进一出视频| 十八禁网站网址无遮挡| 一级片'在线观看视频| 久久久精品区二区三区| 国产精品成人在线| 成年版毛片免费区| 中国美女看黄片| 搡老熟女国产l中国老女人| 一进一出好大好爽视频| av又黄又爽大尺度在线免费看| 国产精品久久久久久精品古装| 精品免费久久久久久久清纯 | 久久久精品94久久精品| 日本黄色视频三级网站网址 | 欧美 亚洲 国产 日韩一| 飞空精品影院首页| 十八禁网站网址无遮挡| 午夜老司机福利片| 色播在线永久视频| 建设人人有责人人尽责人人享有的| 99精品欧美一区二区三区四区| 国产精品秋霞免费鲁丝片| 黑人操中国人逼视频| 亚洲男人天堂网一区| 欧美日韩亚洲国产一区二区在线观看 | 老司机靠b影院| 香蕉丝袜av| 中文字幕av电影在线播放| 亚洲精品一二三| 法律面前人人平等表现在哪些方面| 自线自在国产av| 18禁美女被吸乳视频| 亚洲中文日韩欧美视频| 波多野结衣av一区二区av| 一级,二级,三级黄色视频| 国产99久久九九免费精品| 欧美黑人精品巨大| 在线av久久热| 国产深夜福利视频在线观看| 成年人黄色毛片网站| 757午夜福利合集在线观看| 国产精品一区二区在线观看99| 亚洲情色 制服丝袜| 嫩草影视91久久| 亚洲人成电影观看| 成人精品一区二区免费| 99国产精品免费福利视频| 亚洲avbb在线观看| 国产在线一区二区三区精| 女性生殖器流出的白浆| 制服人妻中文乱码| 久久久久久免费高清国产稀缺| 亚洲精品粉嫩美女一区| 久久热在线av| 宅男免费午夜| 法律面前人人平等表现在哪些方面| 免费少妇av软件| 国产区一区二久久| 精品亚洲成国产av| 黄色片一级片一级黄色片| 国产精品久久久久久精品电影小说| 国产xxxxx性猛交| 亚洲人成伊人成综合网2020| 不卡一级毛片| 香蕉国产在线看| 国产精品久久电影中文字幕 | 波多野结衣一区麻豆| 亚洲精品一卡2卡三卡4卡5卡| 国产精品九九99| 91国产中文字幕| 波多野结衣一区麻豆| 亚洲va日本ⅴa欧美va伊人久久| 日韩视频一区二区在线观看| 看免费av毛片| 日韩欧美一区二区三区在线观看 | 国产欧美日韩一区二区三区在线| 国产av精品麻豆| 久久久久精品国产欧美久久久| 精品一区二区三区视频在线观看免费 | 国产淫语在线视频| 首页视频小说图片口味搜索| av线在线观看网站| 一区二区av电影网| 国产欧美日韩一区二区三区在线| av福利片在线| 韩国精品一区二区三区| bbb黄色大片| 天堂俺去俺来也www色官网| 一级,二级,三级黄色视频| 日韩免费高清中文字幕av| 深夜精品福利| 人人澡人人妻人| 欧美激情 高清一区二区三区| 视频在线观看一区二区三区| 亚洲专区国产一区二区| 久久人妻熟女aⅴ| 日韩免费高清中文字幕av| 人人妻人人澡人人爽人人夜夜| 一本大道久久a久久精品| 99精国产麻豆久久婷婷| 欧美日韩亚洲国产一区二区在线观看 | 久久精品aⅴ一区二区三区四区| 视频区图区小说| 成人国产一区最新在线观看| 一级,二级,三级黄色视频| 一二三四在线观看免费中文在| 亚洲一区二区三区欧美精品| 国产日韩欧美视频二区| 超色免费av| 建设人人有责人人尽责人人享有的| 这个男人来自地球电影免费观看| 欧美日韩精品网址| 色精品久久人妻99蜜桃| 侵犯人妻中文字幕一二三四区| 国产有黄有色有爽视频| 一区二区三区国产精品乱码| 国产三级黄色录像| 免费一级毛片在线播放高清视频 | 国产精品久久久人人做人人爽| 国产精品 欧美亚洲| 亚洲精品久久成人aⅴ小说| 一区在线观看完整版| 亚洲欧美激情在线| 18在线观看网站| 麻豆成人av在线观看| 性高湖久久久久久久久免费观看| 精品国产乱子伦一区二区三区| 一二三四社区在线视频社区8| 伦理电影免费视频| 欧美黑人欧美精品刺激| 久久99热这里只频精品6学生| 日韩欧美一区二区三区在线观看 | aaaaa片日本免费| 亚洲 国产 在线| 别揉我奶头~嗯~啊~动态视频| 欧美国产精品va在线观看不卡| 国产有黄有色有爽视频| 色婷婷久久久亚洲欧美| 亚洲精品国产精品久久久不卡| av国产精品久久久久影院| 国产一区二区激情短视频| 王馨瑶露胸无遮挡在线观看| 精品久久久久久久毛片微露脸| 搡老熟女国产l中国老女人| 性色av乱码一区二区三区2| 不卡av一区二区三区| 日韩人妻精品一区2区三区| 日韩一区二区三区影片| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲专区字幕在线| 亚洲欧美一区二区三区黑人| 国产精品久久久久久精品电影小说| 精品免费久久久久久久清纯 | 精品福利观看| 久久人妻福利社区极品人妻图片| 久久精品国产综合久久久| videosex国产| 亚洲欧美一区二区三区久久| 香蕉丝袜av| 一本大道久久a久久精品| 交换朋友夫妻互换小说| 亚洲五月色婷婷综合| 桃红色精品国产亚洲av| 另类精品久久| 久久久国产成人免费| 自线自在国产av| aaaaa片日本免费| 免费看a级黄色片| 夜夜夜夜夜久久久久| 久久狼人影院| 久久久精品94久久精品| 国产精品免费一区二区三区在线 | 老司机影院毛片| 99九九在线精品视频| 狠狠婷婷综合久久久久久88av| 老司机亚洲免费影院| 91成年电影在线观看| 国产在视频线精品| 肉色欧美久久久久久久蜜桃| 精品少妇久久久久久888优播| 少妇粗大呻吟视频| 麻豆国产av国片精品| 国产欧美日韩一区二区三区在线| 在线av久久热| 国产淫语在线视频| 国产精品一区二区免费欧美| 美女高潮到喷水免费观看| 一区二区三区国产精品乱码| 最新的欧美精品一区二区| 飞空精品影院首页| 亚洲avbb在线观看| 深夜精品福利| 青青草视频在线视频观看| 国产欧美亚洲国产| 国产麻豆69| 97在线人人人人妻| 亚洲精品国产一区二区精华液| 欧美变态另类bdsm刘玥| 精品国产乱子伦一区二区三区| 久久久国产精品麻豆| av国产精品久久久久影院| 黄色视频,在线免费观看| 色视频在线一区二区三区| 美女高潮到喷水免费观看| 男男h啪啪无遮挡| 叶爱在线成人免费视频播放| 真人做人爱边吃奶动态| 久久中文字幕人妻熟女| av国产精品久久久久影院| 国产免费现黄频在线看| 女人爽到高潮嗷嗷叫在线视频| 国产精品麻豆人妻色哟哟久久| 国产亚洲欧美精品永久| 无遮挡黄片免费观看| 法律面前人人平等表现在哪些方面| 好男人电影高清在线观看| 亚洲成a人片在线一区二区| videos熟女内射| 久久精品亚洲av国产电影网| 热re99久久精品国产66热6| 国产一区二区在线观看av| 亚洲精品在线观看二区| 热99re8久久精品国产| 新久久久久国产一级毛片| 亚洲精品国产一区二区精华液| 在线观看免费视频日本深夜| 亚洲欧美一区二区三区黑人| 亚洲人成伊人成综合网2020| 蜜桃国产av成人99|