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

    基于粒子濾波狀態(tài)估計(jì)的滾動(dòng)軸承故障識(shí)別方法

    2017-12-15 02:30:30,
    計(jì)算機(jī)測(cè)量與控制 2017年11期
    關(guān)鍵詞:殘差故障診斷濾波

    ,

    (西北工業(yè)大學(xué) 動(dòng)力與能源學(xué)院,西安 710129)

    基于粒子濾波狀態(tài)估計(jì)的滾動(dòng)軸承故障識(shí)別方法

    史曉雪,吳亞鋒

    (西北工業(yè)大學(xué)動(dòng)力與能源學(xué)院,西安710129)

    提出了一種基于粒子濾波狀態(tài)估計(jì)的滾動(dòng)軸承故障識(shí)別方法,該方法主要包括故障模型建立和故障識(shí)別兩個(gè)步驟;在故障模型建立部分,首先依據(jù)滾動(dòng)軸承不同故障狀態(tài)下的振動(dòng)信號(hào),建立對(duì)應(yīng)的自回歸模型,作為故障模型;在故障識(shí)別部分,將正常狀態(tài)下對(duì)應(yīng)的模型,轉(zhuǎn)化為狀態(tài)空間模型,設(shè)計(jì)粒子濾波器,然后對(duì)不同的故障狀態(tài)進(jìn)行估計(jì),提取其殘差的相關(guān)特征,并結(jié)合模型參數(shù)特征應(yīng)用BP神經(jīng)網(wǎng)絡(luò)識(shí)別算法進(jìn)行故障識(shí)別;最后以美國(guó)凱斯西儲(chǔ)大學(xué)的滾動(dòng)軸承振動(dòng)數(shù)據(jù)為例,驗(yàn)證了該方法的有效性。

    滾動(dòng)軸承;粒子濾波;自回歸模型;狀態(tài)估計(jì)

    0 引言

    滾動(dòng)軸承廣泛應(yīng)用于各種旋轉(zhuǎn)機(jī)械設(shè)備中,而旋轉(zhuǎn)機(jī)械設(shè)備的故障大多數(shù)也是由滾動(dòng)軸承引起的。對(duì)滾動(dòng)軸承進(jìn)行故障識(shí)別,可提高機(jī)械設(shè)備的可靠性與維修性。如何從滾動(dòng)軸承故障振動(dòng)信號(hào)中提取有效的故障特征,實(shí)現(xiàn)早期故障的特征提取,一直是滾動(dòng)軸承故障識(shí)別所面臨的難題。

    機(jī)械設(shè)備故障診斷方法,一般可分為基于解析模型的方法和基于數(shù)據(jù)驅(qū)動(dòng)的方法兩大類(lèi)[1]。基于解析模型的故障診斷方法包括狀態(tài)估計(jì)法,參數(shù)估計(jì)法和等價(jià)空間法,其基本思想是利用系統(tǒng)的數(shù)學(xué)模型和可觀測(cè)的輸入輸出量構(gòu)造殘差信號(hào)來(lái)反映系統(tǒng)期望行為與實(shí)際運(yùn)行模式之間的不一致,然后給予對(duì)殘差信號(hào)的分析進(jìn)行故障診斷。基于數(shù)據(jù)驅(qū)動(dòng)的故障診斷方法可分為機(jī)器學(xué)習(xí)類(lèi)方法、多元統(tǒng)計(jì)分析類(lèi)方法、信號(hào)處理類(lèi)方法等,其基本思想是對(duì)過(guò)程運(yùn)行數(shù)據(jù)進(jìn)行分析處理,從而在不需知道解析模型的情況下完成系統(tǒng)的故障診斷?;诮馕瞿P偷墓收显\斷方法研究較早,對(duì)于線(xiàn)性系統(tǒng)的診斷也較為完善、成果較多。然而,非線(xiàn)型系統(tǒng)的故障診斷依然是當(dāng)前研究的難點(diǎn)和熱點(diǎn)問(wèn)題[3]。

    本文將基于數(shù)據(jù)驅(qū)動(dòng)的識(shí)別方法與基于解析模型的方法相結(jié)合,提出了一種基于粒子濾波的滾動(dòng)軸承故障識(shí)別方法。該方法包含兩個(gè)步驟,即故障模型建立和故障識(shí)別。故障模型建立是依據(jù)原始振動(dòng)信號(hào)建立故障模型,為后續(xù)狀態(tài)估計(jì)奠定基礎(chǔ)。在故障識(shí)別部分,設(shè)計(jì)正常狀態(tài)所對(duì)應(yīng)的粒子濾波器,對(duì)不同故障狀態(tài)進(jìn)行估計(jì),提取其殘差的相關(guān)特征,并結(jié)合模型參數(shù)特征應(yīng)用BP神經(jīng)網(wǎng)絡(luò)識(shí)別算法進(jìn)行故障識(shí)別。

    1 粒子濾波算法介紹

    粒子濾波算法(particle filtering,PF)源于蒙特卡洛(Monte carlo)的思想,通過(guò)非參數(shù)化的蒙特卡洛模擬方法來(lái)實(shí)現(xiàn)遞推貝葉斯濾波。適用于任何能用狀態(tài)空間模型描述的非線(xiàn)性系統(tǒng),精度可以逼近最優(yōu)估計(jì)。

    任一動(dòng)態(tài)系統(tǒng)的狀態(tài)空間模型可以表示為:

    xk=f(xk-1,vk)

    (1)

    yk=h(xk,nk)

    (2)

    上面兩式中,式(1)所描述的是狀態(tài)轉(zhuǎn)移方程,式(2)所描述的是觀測(cè)方程。xk為系統(tǒng)狀態(tài)在k時(shí)刻對(duì)應(yīng)的狀態(tài)變量,yk為觀測(cè)值,表示的是xk對(duì)應(yīng)的觀測(cè)值,vk為系統(tǒng)噪聲,nk為觀測(cè)噪聲,vk和nk相互獨(dú)立。

    (3)

    (4)

    假設(shè)重要密度可以分解為:

    q(x0:k|z1:k)=q(xk|x0:k-1,z1:k)q(x0:k-1|z1:k-1)

    (5)

    由于后驗(yàn)概率密度函數(shù)可表示為:

    ∝p(zk|xk)p(xk|xk-1)p(x0:k-1|z1:k-1)

    (6)

    將式(6)和式(5)帶入(4),可得重要性權(quán)值更新公式為:

    (7)

    (8)

    標(biāo)準(zhǔn)粒子濾波算法選擇易于實(shí)現(xiàn)的先驗(yàn)概率密度作為重要密度函數(shù),即:

    (9)

    將式(9)帶入式(8),重要性權(quán)值可簡(jiǎn)化為:

    (10)

    (11)

    而后驗(yàn)概率密度p(xk|z1:k)可表示為:

    (12)

    2 粒子濾波狀態(tài)估計(jì)滾動(dòng)軸承故障識(shí)別

    本文提出的基于粒子濾波狀態(tài)估計(jì)的故障識(shí)別方法流程圖如圖1所示,包括故障模型建立和故障識(shí)別兩個(gè)步驟。其基本思想是:滾動(dòng)軸承在正常狀態(tài)下,通過(guò)粒子濾波算法得到的狀態(tài)估計(jì)值與真實(shí)值的差值應(yīng)該比較??;當(dāng)處于故障狀態(tài)下時(shí),狀態(tài)估計(jì)值與真實(shí)值之間的差值較大。不同故障所獲得的差值不同,作為故障識(shí)別的依據(jù)。

    圖1 粒子濾波狀態(tài)估計(jì)故障識(shí)別方法流程圖

    2.1 故障模型建立

    在故障模型建立過(guò)程中,采用自回歸(AR)模型建立故障模型。其數(shù)學(xué)表達(dá)式為:

    x(n)=a1x(n-1)+a2x(n-2)+…+apx(n-p)+e(n)

    (13)

    其中:a為模型參數(shù);p為模型階數(shù);e(n)為模型殘差且服從標(biāo)準(zhǔn)正態(tài)分布。

    AR建模主要包括階數(shù)確定和參數(shù)求解兩個(gè)部分,其中,定階的主要方法有AIC,F(xiàn)PE,BIC等定階準(zhǔn)則;參數(shù)求解方法主要有最小二乘法,最大似然估計(jì)、矩估計(jì)、Levinson-Durbin算法等。本文選取AIC準(zhǔn)則確定模型階數(shù),LD算法計(jì)算模型參數(shù)。

    (1)AIC準(zhǔn)則模型定階。

    AIC準(zhǔn)則的一般表示形式為:

    (14)

    (2)模型參數(shù)求解。

    (15)

    記e(n)為預(yù)測(cè)值和真值之間的殘差,則:

    (16)

    因此總的預(yù)測(cè)誤差功率為:

    (17)

    采用Levinson-Durbin遞推算法計(jì)算模型參數(shù):

    (18)

    am(k)=am-1(k)+kmam-1(m-k)

    (19)

    (20)

    其中:rx(m)為x(n)的自相關(guān)函數(shù),ρm為m階時(shí)的前向預(yù)測(cè)的最小誤差功率。

    當(dāng)m=1時(shí),有:

    a1(1)=-rx(1)/rx(0)

    (21)

    (22)

    定義初始條件為:

    ρ0=rx(0)

    (23)

    則有:

    (24)

    上述條件為基于L-D算法求AR模型參數(shù)提供了初始值。L-D算法從低階開(kāi)始遞推,直到階次p,給出了在每一個(gè)階次時(shí)的所有參數(shù),即am(1),am(2),…,am(m),其中m=1,2,…,p。

    2.2 故障識(shí)別

    采用正常狀態(tài)所對(duì)應(yīng)的模型作為粒子濾波狀態(tài)空間方程,假設(shè)觀測(cè)值yk和xk存在線(xiàn)性關(guān)系,則根據(jù)狀態(tài)空間模型表達(dá)式(1)和(2),可構(gòu)造粒子濾波器如下:

    (25)

    yk=xk+vk

    (26)

    (1)狀態(tài)預(yù)測(cè)

    (4)更新粒子權(quán)值并歸一化。

    (5)粒子重采樣。

    在正常狀態(tài)下,殘差rk比較小,接近于0,在某種故障狀態(tài)下,殘差rk較大。根據(jù)不同故障狀態(tài)下殘差值的大小,計(jì)算其故障特征,其特征指標(biāo)如表1所示:

    表1 殘差相關(guān)特征

    其中,M為狀態(tài)估計(jì)步長(zhǎng)。

    3 仿真實(shí)驗(yàn)

    3.1 實(shí)驗(yàn)介紹

    為了驗(yàn)證本文提出方法的有效性,采用美國(guó)凱斯西儲(chǔ)大學(xué)電氣工程實(shí)驗(yàn)室的滾動(dòng)軸承實(shí)驗(yàn)平臺(tái)采集到的故障數(shù)據(jù)對(duì)該方法進(jìn)行驗(yàn)證。滾動(dòng)軸承故障模擬試驗(yàn)臺(tái)如圖所示。該裝置有一臺(tái)2英制馬力(hp)的電動(dòng)機(jī)在圖左部,一個(gè)扭矩編碼器/傳感器在圖中部,一個(gè)功率測(cè)試計(jì)在圖右部,以及電子控制器組成。實(shí)驗(yàn)軸承的轉(zhuǎn)速分別為1 797 rpm,1 772 rpm,1 750 rpm,1 730 rpm,振動(dòng)信號(hào)采樣頻率為12 kHz。軸承采用電火花加工的方式引入單點(diǎn)缺陷來(lái)模擬軸承故障,故障類(lèi)型包括內(nèi)圈故障,外圈故障和滾動(dòng)體故障這三種,損傷的尺寸分別為0.007、0.014、0.021、0.028英寸。

    圖2 凱斯西儲(chǔ)大學(xué)電氣工程實(shí)驗(yàn)室的滾動(dòng)軸承實(shí)驗(yàn)平臺(tái)

    3.2 對(duì)比分析

    采用1797 rpm下,損傷尺寸為0.007英寸對(duì)應(yīng)下的三種軸承故障數(shù)據(jù)進(jìn)行故障識(shí)別。對(duì)于每種故障狀態(tài),選取1~256個(gè)點(diǎn)進(jìn)行AR建模,依據(jù)AIC準(zhǔn)則確定模型階數(shù)為6階,計(jì)算得到模型參數(shù)如表2所示,將該參數(shù)帶入式(13)可得到四種故障模型。

    表2 不同狀態(tài)下模型參數(shù)

    將正常狀態(tài)下對(duì)應(yīng)的AR模型作為狀態(tài)空間模型,設(shè)計(jì)粒子濾波器,對(duì)不同故障狀態(tài)進(jìn)行估計(jì)。其中,粒子數(shù)N=100,估計(jì)步長(zhǎng)T=256?,F(xiàn)截取50步長(zhǎng)的殘差結(jié)果,如圖3所示。

    圖3 不同故障狀態(tài)殘差圖

    由表1所示公式,計(jì)算該殘差特征的絕對(duì)均值、標(biāo)準(zhǔn)差、方差及峰值,得到4種特征指標(biāo)。每種故障選取1組特征指標(biāo),如表3所示。

    表3 殘差相關(guān)特征

    將殘差相關(guān)特征結(jié)合參數(shù)特征作為BP神經(jīng)網(wǎng)絡(luò)的輸入進(jìn)行故障識(shí)別。采用3層的BP神經(jīng)網(wǎng)絡(luò),網(wǎng)絡(luò)結(jié)構(gòu)為N(10,6,3)。取顯示間隔為50,學(xué)習(xí)步長(zhǎng)0.3,最小均方誤差為0.1,最大訓(xùn)練次數(shù)為1000。

    為了驗(yàn)證該時(shí)域特征參數(shù)在故障識(shí)別中的優(yōu)勢(shì),現(xiàn)選取原始時(shí)域特征與該特征做比較。原始時(shí)域特征信息如表4所示。

    表4 原始時(shí)域特征指標(biāo)

    對(duì)表4所示原始時(shí)域特征采用相同的識(shí)別方法進(jìn)行故障識(shí)別,并將本文方法與其進(jìn)行對(duì)比,結(jié)果如表5所示。

    表5 不同特征的故障識(shí)別率

    從表5可以看出,三種故障的識(shí)別率由高到低依次為:外圈故障、內(nèi)圈故障和滾珠故障。結(jié)合圖3可知,殘差越大故障的識(shí)別率越高。對(duì)于內(nèi)圈故障,兩種方法所得的故障識(shí)別率基本一致,對(duì)于滾珠故障和外圈故障,本文方法所得的故障識(shí)別率較高??梢?jiàn),參數(shù)及殘差的相關(guān)特征相對(duì)于原始時(shí)域特征具有敏感性高的特點(diǎn)。因此,在滾動(dòng)軸承故障識(shí)別中,采用本文方法可以提高故障識(shí)別的準(zhǔn)確性。

    4 結(jié)論

    本文提出了一種基于粒子濾波狀態(tài)估計(jì)的滾動(dòng)軸承故障識(shí)別方法。對(duì)滾動(dòng)軸承的不同故障狀態(tài)建立故障模型,通過(guò)粒子濾波對(duì)每種故障狀態(tài)進(jìn)行估計(jì),進(jìn)而獲得不同故障狀態(tài)下的殘差的相關(guān)特征,將其與模型參數(shù)特征及BP神經(jīng)網(wǎng)絡(luò)相結(jié)合,應(yīng)用于滾動(dòng)軸的故障識(shí)別。為了驗(yàn)證本文所提方法的有效性,采用美國(guó)西儲(chǔ)大學(xué)滾動(dòng)軸承故障數(shù)據(jù)進(jìn)行分析,結(jié)果表明該方法具有識(shí)別率高的優(yōu)點(diǎn),可以在滾動(dòng)軸承故障識(shí)別中較好地應(yīng)用。

    [1] 張三同,郝晶晶. 多模型粒子濾波的故障診斷方法研究[J]. 控制工程, 2012, 19(5): 864-869.

    [2] 夏 飛, 郝碩濤, 張 濤, 等. 改進(jìn)粒子濾波在汽輪機(jī)故障診斷中的應(yīng)用[J]. 計(jì)算機(jī)測(cè)量與控制, 2016, 24(1): 35-38.

    [3] 張玲霞, 劉志倉(cāng), 王 輝,等. 非線(xiàn)性系統(tǒng)故障診斷的粒子濾波方法[J]. 電子學(xué)報(bào), 2015, 43(3): 615-619.

    [4] Alrowaie F, Kwok KE, Gopaluni RB. An algorithm for fault detection in stochastic non-linear state-space models using particle filters[J]. Control Engineering Practice, 2012, 20(10):1016-1032.

    [5] Zhang LX, Liu ZC, Wang H.Particle Filter Method for Fault Diagnosis in Nonlinear System[J]. Acta Electronica Sinica , 2015, 43(4): 615-619.

    [6] 杜京義,殷夢(mèng)鑫. 一種改進(jìn)粒子濾波算法應(yīng)用于故障診斷[J]. 系統(tǒng)仿真學(xué)報(bào), 2014, 26(1): 62-66.

    [7] 李宏坤, 楊 蕊, 任遠(yuǎn)杰,等. 利用粒子濾波與譜峭度的滾動(dòng)軸承故障診斷[J]. 機(jī)械工程學(xué)報(bào), 2017, 53(3): 63-72.

    [8] 姜洪開(kāi), 王 涵,郭 棟. 基于改進(jìn)粒子濾波的滾動(dòng)軸承故障識(shí)別[J]. 振動(dòng)與沖擊, 2014, 33(8): 132-134.

    FaultDiagnosisofRollingBearingBasedonParticleFilterStateEstimation

    Shi Xiaoxue,Wu Yafeng

    (College of Power and Energy, Northwestern Polytechnical University, Xi’an 710129, China)

    A fault diagnosis method of rolling bearing based on particle filter state estimation is proposed. The method mainly includes two steps: fault model establishment and fault identification. In the fault model establishment part, the corresponding autoregressive model is established according to the vibration signals of the rolling bearings in different fault states, which is used as the fault model; In the fault recognition part, the corresponding model in the normal state is transformed into the state space model, and the particle filter is designed. Then, the related features of the residuals are extracted by estimating the different fault states. Combined with the model parameters, the BP neural network identification algorithm is used to identify the faults. Finally, the validity of the method is verified by the vibration data of rolling bearing of Case Western Reserve University.

    rolling bearings; autoregressive model; particle filter; state estimation

    2017-08-15;

    2017-09-14。

    史曉雪(1992-),女,甘肅白銀人,碩士研究生,主要從事信號(hào)與信息處理方向的研究。

    吳亞鋒(1961-),男,陜西西安人,教授,博士生導(dǎo)師,主要從事現(xiàn)代信號(hào)處理理論與方法、計(jì)算機(jī)測(cè)控技術(shù)、振動(dòng)噪聲分析與控制方向的研究。

    1671-4598(2017)11-0030-04

    10.16526/j.cnki.11-4762/tp.2017.11.008

    TP391

    A

    猜你喜歡
    殘差故障診斷濾波
    基于雙向GRU與殘差擬合的車(chē)輛跟馳建模
    基于殘差學(xué)習(xí)的自適應(yīng)無(wú)人機(jī)目標(biāo)跟蹤算法
    基于遞歸殘差網(wǎng)絡(luò)的圖像超分辨率重建
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    RTS平滑濾波在事后姿態(tài)確定中的應(yīng)用
    基于線(xiàn)性正則變換的 LMS 自適應(yīng)濾波
    平穩(wěn)自相關(guān)過(guò)程的殘差累積和控制圖
    河南科技(2015年8期)2015-03-11 16:23:52
    基于LCD和排列熵的滾動(dòng)軸承故障診斷
    基于WPD-HHT的滾動(dòng)軸承故障診斷
    高速泵的故障診斷
    河南科技(2014年3期)2014-02-27 14:05:48
    最新中文字幕久久久久| 美女视频免费永久观看网站| 亚洲av免费高清在线观看| 春色校园在线视频观看| 狂野欧美白嫩少妇大欣赏| 韩国高清视频一区二区三区| 另类精品久久| 国产精品一区二区在线观看99| 亚洲成色77777| 精品人妻偷拍中文字幕| freevideosex欧美| 看非洲黑人一级黄片| 欧美日韩精品成人综合77777| 少妇人妻一区二区三区视频| 寂寞人妻少妇视频99o| 18禁在线播放成人免费| 亚洲激情五月婷婷啪啪| 欧美精品一区二区大全| 在线观看免费日韩欧美大片 | 精品人妻熟女毛片av久久网站| 性高湖久久久久久久久免费观看| 亚洲欧美精品自产自拍| 一级,二级,三级黄色视频| 观看av在线不卡| 亚洲伊人久久精品综合| 高清av免费在线| 天美传媒精品一区二区| 国内精品宾馆在线| 丝袜在线中文字幕| 国产成人精品一,二区| 成人美女网站在线观看视频| 中文字幕人妻丝袜制服| 国产男人的电影天堂91| 水蜜桃什么品种好| 日韩一区二区视频免费看| 久久精品国产亚洲av涩爱| 又黄又爽又刺激的免费视频.| 亚洲精品乱久久久久久| av福利片在线| 成人黄色视频免费在线看| 婷婷色麻豆天堂久久| 国产伦精品一区二区三区视频9| 男人舔奶头视频| 街头女战士在线观看网站| 少妇人妻一区二区三区视频| 一级毛片 在线播放| 国产精品欧美亚洲77777| 啦啦啦中文免费视频观看日本| 免费观看的影片在线观看| 国产日韩欧美视频二区| 亚洲图色成人| 国产日韩欧美亚洲二区| 久久久欧美国产精品| 国产白丝娇喘喷水9色精品| 亚洲精品aⅴ在线观看| 深夜a级毛片| www.色视频.com| 少妇人妻 视频| 色哟哟·www| 久久精品久久精品一区二区三区| 国产在视频线精品| 秋霞在线观看毛片| 新久久久久国产一级毛片| 成人午夜精彩视频在线观看| 女人久久www免费人成看片| 大话2 男鬼变身卡| 日本色播在线视频| 又粗又硬又长又爽又黄的视频| 校园人妻丝袜中文字幕| 我要看黄色一级片免费的| 国产无遮挡羞羞视频在线观看| 久久国产精品男人的天堂亚洲 | 一个人免费看片子| 精品国产一区二区三区久久久樱花| 亚洲伊人久久精品综合| 国产精品嫩草影院av在线观看| 在线观看一区二区三区激情| 国产免费一级a男人的天堂| 精品午夜福利在线看| 一本一本综合久久| 国产在线男女| 亚洲精品日韩av片在线观看| 最后的刺客免费高清国语| av在线app专区| 美女脱内裤让男人舔精品视频| 国产精品国产三级国产av玫瑰| 国产亚洲精品久久久com| 亚洲三级黄色毛片| av在线观看视频网站免费| 一级二级三级毛片免费看| 国产无遮挡羞羞视频在线观看| 亚洲欧美日韩东京热| 99热国产这里只有精品6| 美女福利国产在线| 亚洲国产精品国产精品| 国产女主播在线喷水免费视频网站| 国产国拍精品亚洲av在线观看| 成年人午夜在线观看视频| 日本91视频免费播放| 亚洲欧美日韩东京热| 国产黄色免费在线视频| 日韩精品有码人妻一区| 亚洲综合色惰| 一二三四中文在线观看免费高清| 一本久久精品| 国产深夜福利视频在线观看| 免费看日本二区| 国产免费一级a男人的天堂| 亚洲第一av免费看| 欧美老熟妇乱子伦牲交| 大陆偷拍与自拍| 黄片无遮挡物在线观看| 亚洲精品色激情综合| 最近中文字幕高清免费大全6| 久久久久网色| 成年av动漫网址| 久久狼人影院| 欧美激情极品国产一区二区三区 | 一本一本综合久久| 免费黄网站久久成人精品| a级片在线免费高清观看视频| 老司机影院毛片| 日韩精品有码人妻一区| 18禁在线播放成人免费| videos熟女内射| 制服丝袜香蕉在线| 日韩成人av中文字幕在线观看| 亚洲精品国产av蜜桃| 高清午夜精品一区二区三区| 亚洲国产精品一区三区| 亚洲av在线观看美女高潮| 国产av码专区亚洲av| 777米奇影视久久| 国产美女午夜福利| 成人亚洲精品一区在线观看| 爱豆传媒免费全集在线观看| 久久99热6这里只有精品| 国产在线视频一区二区| 中文精品一卡2卡3卡4更新| 国产黄色视频一区二区在线观看| 欧美日本中文国产一区发布| 自拍欧美九色日韩亚洲蝌蚪91 | 婷婷色综合大香蕉| 久久国内精品自在自线图片| 久久av网站| 最新中文字幕久久久久| 最近2019中文字幕mv第一页| 成人免费观看视频高清| 免费久久久久久久精品成人欧美视频 | 国产精品不卡视频一区二区| 美女国产视频在线观看| 少妇猛男粗大的猛烈进出视频| 国产黄色视频一区二区在线观看| 夫妻性生交免费视频一级片| 久久久久人妻精品一区果冻| 国产欧美日韩一区二区三区在线 | 黑人高潮一二区| 国内精品宾馆在线| 亚洲欧美日韩东京热| 女人精品久久久久毛片| 99热全是精品| 熟女电影av网| 大片电影免费在线观看免费| 美女cb高潮喷水在线观看| 又爽又黄a免费视频| 18+在线观看网站| 欧美 亚洲 国产 日韩一| 国产乱人偷精品视频| 色94色欧美一区二区| 亚洲天堂av无毛| 狂野欧美激情性xxxx在线观看| 国产黄色视频一区二区在线观看| 国产亚洲av片在线观看秒播厂| 久久久久久久久久久免费av| 看非洲黑人一级黄片| 国产综合精华液| 国产无遮挡羞羞视频在线观看| 国产成人一区二区在线| 亚洲电影在线观看av| 一本—道久久a久久精品蜜桃钙片| 欧美 亚洲 国产 日韩一| 久久狼人影院| av播播在线观看一区| 久久午夜福利片| 99热这里只有是精品50| 91精品一卡2卡3卡4卡| 欧美三级亚洲精品| 欧美精品国产亚洲| 亚洲av免费高清在线观看| 人人妻人人看人人澡| 国产精品偷伦视频观看了| 校园人妻丝袜中文字幕| 国产在线一区二区三区精| 成人免费观看视频高清| 亚洲av成人精品一二三区| 国产成人免费无遮挡视频| 欧美xxxx性猛交bbbb| 丁香六月天网| 国产精品一区二区三区四区免费观看| 国产精品偷伦视频观看了| 搡女人真爽免费视频火全软件| 97在线人人人人妻| 日韩欧美一区视频在线观看 | 国产女主播在线喷水免费视频网站| 街头女战士在线观看网站| 最新的欧美精品一区二区| 国产精品秋霞免费鲁丝片| 一边亲一边摸免费视频| 午夜福利影视在线免费观看| 熟妇人妻不卡中文字幕| 欧美日本中文国产一区发布| 久久99蜜桃精品久久| 欧美性感艳星| 亚洲欧美日韩另类电影网站| 久久婷婷青草| 久久热精品热| 又粗又硬又长又爽又黄的视频| 2021少妇久久久久久久久久久| av在线app专区| 欧美日韩综合久久久久久| 26uuu在线亚洲综合色| 一区二区三区四区激情视频| 18禁动态无遮挡网站| 欧美日本中文国产一区发布| 婷婷色综合大香蕉| 国产成人免费无遮挡视频| 丝瓜视频免费看黄片| 亚洲中文av在线| 女人精品久久久久毛片| 老司机影院成人| 在线观看免费视频网站a站| 亚洲精品久久久久久婷婷小说| 麻豆成人av视频| 色94色欧美一区二区| 自拍偷自拍亚洲精品老妇| 中国国产av一级| 一级黄片播放器| 97超碰精品成人国产| 男女边摸边吃奶| 亚洲精品乱码久久久久久按摩| 777米奇影视久久| 久久久久网色| 赤兔流量卡办理| 十八禁网站网址无遮挡 | 卡戴珊不雅视频在线播放| 亚洲色图综合在线观看| 在线 av 中文字幕| 亚洲三级黄色毛片| 国产女主播在线喷水免费视频网站| 有码 亚洲区| 国产成人aa在线观看| 成人午夜精彩视频在线观看| 天天操日日干夜夜撸| 18禁裸乳无遮挡动漫免费视频| 春色校园在线视频观看| 国产日韩欧美在线精品| 午夜福利在线观看免费完整高清在| 一本色道久久久久久精品综合| 国产精品熟女久久久久浪| 丝袜喷水一区| 高清黄色对白视频在线免费看 | 欧美97在线视频| 熟女电影av网| 色婷婷av一区二区三区视频| 亚洲va在线va天堂va国产| 久久久a久久爽久久v久久| 午夜福利在线观看免费完整高清在| 欧美精品国产亚洲| 少妇 在线观看| 国产精品不卡视频一区二区| 国产精品无大码| 中文字幕av电影在线播放| 观看免费一级毛片| 日韩人妻高清精品专区| 偷拍熟女少妇极品色| 欧美丝袜亚洲另类| 亚洲精品中文字幕在线视频 | 日韩中字成人| 亚洲精品日本国产第一区| 国产成人freesex在线| 亚洲av福利一区| 晚上一个人看的免费电影| 日本-黄色视频高清免费观看| 亚洲国产精品专区欧美| 观看美女的网站| 欧美日韩亚洲高清精品| 久久av网站| freevideosex欧美| 国精品久久久久久国模美| 久久久国产欧美日韩av| av网站免费在线观看视频| 中国国产av一级| 狠狠精品人妻久久久久久综合| 欧美三级亚洲精品| 在线观看国产h片| 亚洲精品乱码久久久v下载方式| 日本黄大片高清| 成人漫画全彩无遮挡| 高清在线视频一区二区三区| 亚洲精品自拍成人| 在线观看国产h片| 老熟女久久久| 欧美精品一区二区大全| 天堂8中文在线网| 精品亚洲成国产av| 国产精品国产三级专区第一集| 美女内射精品一级片tv| 国产av一区二区精品久久| 日本午夜av视频| 一级爰片在线观看| 国产精品一区二区性色av| 在线观看三级黄色| 国产成人aa在线观看| 亚洲av综合色区一区| 欧美老熟妇乱子伦牲交| 美女国产视频在线观看| 国模一区二区三区四区视频| 午夜影院在线不卡| 街头女战士在线观看网站| a级毛片在线看网站| 国产女主播在线喷水免费视频网站| 人人澡人人妻人| 乱人伦中国视频| .国产精品久久| av一本久久久久| 亚洲精品乱码久久久v下载方式| 欧美三级亚洲精品| 激情五月婷婷亚洲| 日韩一区二区三区影片| 国产精品不卡视频一区二区| 国产欧美日韩综合在线一区二区 | 亚洲国产精品成人久久小说| 日韩欧美一区视频在线观看 | 成人无遮挡网站| 日韩视频在线欧美| 久久久午夜欧美精品| 国产成人a∨麻豆精品| 国产色婷婷99| av播播在线观看一区| 99九九在线精品视频 | 免费av中文字幕在线| 亚洲国产欧美日韩在线播放 | 男人和女人高潮做爰伦理| 久热这里只有精品99| 桃花免费在线播放| 免费黄网站久久成人精品| av福利片在线| 亚洲成人一二三区av| 国产精品偷伦视频观看了| 麻豆乱淫一区二区| 午夜日本视频在线| 日韩 亚洲 欧美在线| 国产免费一级a男人的天堂| 精品卡一卡二卡四卡免费| 中国美白少妇内射xxxbb| 亚洲欧洲精品一区二区精品久久久 | a 毛片基地| 永久网站在线| 99久久精品国产国产毛片| 中文字幕人妻丝袜制服| 国语对白做爰xxxⅹ性视频网站| 成人亚洲欧美一区二区av| 午夜福利视频精品| 伦理电影大哥的女人| 久久6这里有精品| 十八禁网站网址无遮挡 | 夜夜看夜夜爽夜夜摸| 国产在视频线精品| 欧美日韩综合久久久久久| 成人无遮挡网站| 最近中文字幕2019免费版| 一级毛片我不卡| 亚洲成人av在线免费| 午夜91福利影院| 精品国产一区二区三区久久久樱花| av国产久精品久网站免费入址| 国产男女超爽视频在线观看| 国模一区二区三区四区视频| 欧美日韩综合久久久久久| 亚洲va在线va天堂va国产| 一级毛片黄色毛片免费观看视频| 国产在视频线精品| 精品熟女少妇av免费看| 热99国产精品久久久久久7| 亚洲第一区二区三区不卡| 热99国产精品久久久久久7| 国产免费福利视频在线观看| 男人和女人高潮做爰伦理| 亚洲av日韩在线播放| 男人爽女人下面视频在线观看| 大话2 男鬼变身卡| 久久久久精品性色| 日韩一本色道免费dvd| 亚洲精品成人av观看孕妇| 国产黄色免费在线视频| 建设人人有责人人尽责人人享有的| 久久99蜜桃精品久久| 免费看日本二区| 久久精品国产自在天天线| 精品99又大又爽又粗少妇毛片| 久久精品久久久久久噜噜老黄| 三上悠亚av全集在线观看 | 成人综合一区亚洲| 久久人人爽人人爽人人片va| 国产欧美日韩精品一区二区| 色哟哟·www| 午夜免费观看性视频| 日韩强制内射视频| 在线播放无遮挡| 观看美女的网站| 大香蕉久久网| 国产免费一区二区三区四区乱码| 美女脱内裤让男人舔精品视频| 欧美 亚洲 国产 日韩一| 永久免费av网站大全| 水蜜桃什么品种好| 蜜臀久久99精品久久宅男| 日本黄色片子视频| 免费观看a级毛片全部| 亚洲美女搞黄在线观看| 极品人妻少妇av视频| 乱码一卡2卡4卡精品| 在线观看三级黄色| 人妻 亚洲 视频| 国产亚洲最大av| 尾随美女入室| 久久6这里有精品| 成人亚洲欧美一区二区av| av不卡在线播放| 国产女主播在线喷水免费视频网站| 狂野欧美白嫩少妇大欣赏| 亚洲精品,欧美精品| 香蕉精品网在线| 免费播放大片免费观看视频在线观看| av在线app专区| 中文精品一卡2卡3卡4更新| 美女大奶头黄色视频| 人妻系列 视频| 亚洲天堂av无毛| 亚洲图色成人| 国产欧美另类精品又又久久亚洲欧美| 在线观看免费视频网站a站| 久久午夜综合久久蜜桃| 国产一区二区三区综合在线观看 | av天堂久久9| 免费看光身美女| 王馨瑶露胸无遮挡在线观看| 制服丝袜香蕉在线| 少妇 在线观看| 国产 一区精品| 亚洲精品国产av蜜桃| 女人精品久久久久毛片| 久久婷婷青草| 一级毛片久久久久久久久女| 亚洲中文av在线| 一本大道久久a久久精品| 国产永久视频网站| 麻豆乱淫一区二区| 乱系列少妇在线播放| 国产女主播在线喷水免费视频网站| 久久人人爽人人片av| 欧美日韩精品成人综合77777| 亚洲av在线观看美女高潮| 国产精品人妻久久久影院| 国产欧美日韩精品一区二区| 黄色怎么调成土黄色| 黄色日韩在线| 黄色视频在线播放观看不卡| 国产一区二区三区av在线| 九九爱精品视频在线观看| 日韩中文字幕视频在线看片| 午夜免费男女啪啪视频观看| 久久6这里有精品| 欧美日韩一区二区视频在线观看视频在线| 51国产日韩欧美| 亚洲真实伦在线观看| 久久久久久人妻| 久久久久视频综合| 国产成人午夜福利电影在线观看| 久久国产亚洲av麻豆专区| 永久网站在线| 欧美3d第一页| av免费在线看不卡| 国产午夜精品一二区理论片| 91精品国产国语对白视频| 肉色欧美久久久久久久蜜桃| 在线免费观看不下载黄p国产| 成人黄色视频免费在线看| 国产日韩欧美亚洲二区| 国产精品欧美亚洲77777| 亚洲欧美日韩东京热| 日韩免费高清中文字幕av| 亚洲一区二区三区欧美精品| 欧美成人午夜免费资源| 午夜久久久在线观看| 国产一区二区三区综合在线观看 | 高清不卡的av网站| 精品熟女少妇av免费看| 国产69精品久久久久777片| 日韩制服骚丝袜av| 亚洲人成网站在线观看播放| 午夜福利视频精品| 国产永久视频网站| 国产爽快片一区二区三区| 国产成人精品婷婷| 日本黄色日本黄色录像| 大香蕉97超碰在线| 只有这里有精品99| 精品亚洲乱码少妇综合久久| √禁漫天堂资源中文www| 日韩,欧美,国产一区二区三区| 久久国内精品自在自线图片| 黄片无遮挡物在线观看| tube8黄色片| 欧美性感艳星| 一级毛片黄色毛片免费观看视频| 青春草国产在线视频| 亚洲综合色惰| 97超碰精品成人国产| 久久午夜福利片| 日韩av在线免费看完整版不卡| 熟妇人妻不卡中文字幕| 亚洲欧洲日产国产| 久久久国产精品麻豆| 欧美精品亚洲一区二区| 少妇人妻 视频| 国产伦精品一区二区三区四那| 国产黄片视频在线免费观看| 亚洲欧美一区二区三区黑人 | 久久久久久伊人网av| 另类亚洲欧美激情| 卡戴珊不雅视频在线播放| 欧美精品一区二区大全| 成年av动漫网址| 午夜久久久在线观看| 伦理电影大哥的女人| 欧美精品亚洲一区二区| 精品国产一区二区久久| 人人妻人人看人人澡| 久久午夜福利片| 国产黄色免费在线视频| 国产日韩欧美亚洲二区| 美女xxoo啪啪120秒动态图| 成人亚洲欧美一区二区av| 免费在线观看成人毛片| 亚洲伊人久久精品综合| 国产精品不卡视频一区二区| 人妻人人澡人人爽人人| 在线观看免费日韩欧美大片 | 男女免费视频国产| 中文字幕亚洲精品专区| 国产精品免费大片| 亚洲经典国产精华液单| 欧美xxxx性猛交bbbb| 少妇人妻 视频| 赤兔流量卡办理| 国模一区二区三区四区视频| 自拍偷自拍亚洲精品老妇| 国产色爽女视频免费观看| 久久综合国产亚洲精品| 中文字幕人妻丝袜制服| 久热久热在线精品观看| 一级片'在线观看视频| 国产精品一区二区三区四区免费观看| 国产精品久久久久久精品电影小说| 精品国产乱码久久久久久小说| 亚洲丝袜综合中文字幕| 少妇猛男粗大的猛烈进出视频| 男女边摸边吃奶| 中文字幕人妻熟人妻熟丝袜美| 久久综合国产亚洲精品| 高清黄色对白视频在线免费看 | 国产亚洲91精品色在线| 成年人免费黄色播放视频 | 国产有黄有色有爽视频| 免费看av在线观看网站| 久久久国产一区二区| 国产黄片美女视频| 国产乱来视频区| 亚洲一级一片aⅴ在线观看| 免费看光身美女| 欧美精品亚洲一区二区| 亚洲av男天堂| 人人妻人人爽人人添夜夜欢视频 | 久久国产精品大桥未久av | 99久久中文字幕三级久久日本| 午夜免费男女啪啪视频观看| 六月丁香七月| 日韩av不卡免费在线播放| 亚洲国产精品国产精品| 伦理电影免费视频| 久久久久久久久久久免费av| 国产精品99久久99久久久不卡 | 在线看a的网站| 男人狂女人下面高潮的视频| 最近最新中文字幕免费大全7| 国产成人午夜福利电影在线观看| 亚洲婷婷狠狠爱综合网| 男人舔奶头视频| 丁香六月天网| 青春草亚洲视频在线观看| 欧美最新免费一区二区三区| 老司机亚洲免费影院| 哪个播放器可以免费观看大片| av线在线观看网站| 国产视频首页在线观看| 亚洲欧美日韩东京热| 国产淫片久久久久久久久| 亚洲真实伦在线观看| 久久久国产欧美日韩av| 最近中文字幕高清免费大全6| 久久99一区二区三区| 亚洲国产成人一精品久久久| 九色成人免费人妻av| 精品酒店卫生间|