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

    流形學(xué)習(xí)和M—KH—SVR的滾動(dòng)軸承衰退預(yù)測

    2018-01-04 12:08:10許迪葛江華王亞萍
    振動(dòng)工程學(xué)報(bào) 2018年5期
    關(guān)鍵詞:滾動(dòng)軸承特征提取故障診斷

    許迪 葛江華 王亞萍

    摘要: 針對滾動(dòng)軸承中存在數(shù)據(jù)樣本量大、非平穩(wěn)信號波動(dòng)復(fù)雜等問題,提出基于流形學(xué)習(xí)和M-KH-SVR(Multivariable-Krill Herd-Support Vector Regression )的滾動(dòng)軸承衰退預(yù)測方法。該方法首先提取了滾動(dòng)軸承的時(shí)域和頻域特征,組成初始特征向量;然后利用相關(guān)度量系數(shù)(Multiple Correlation Coefficient, MCC)對初始特征進(jìn)行篩選,得到相關(guān)程度較高的特征向量集,并通過局部線性嵌入(Locally Linear Embedding, LLE)方法進(jìn)行特征降維,進(jìn)而組成新的故障特征集;最后將磷蝦群算法引入到多變量支持向量回歸機(jī)中,并對其參數(shù)c和σ進(jìn)行優(yōu)化,利用磷蝦群局部尋優(yōu)和全局尋優(yōu)的能力,提高了參數(shù)選擇效率。通過對多變量特征進(jìn)行實(shí)驗(yàn)對比分析,結(jié)果表明該方法與傳統(tǒng)單一參數(shù)及多特征參數(shù)方法相比,具有良好的泛化性,大幅度提高了運(yùn)算效率和預(yù)測精度,對滾動(dòng)軸承的衰退階段劃分更加精確。

    關(guān)鍵詞: 故障診斷; 滾動(dòng)軸承; 特征提?。?磷蝦群算法; 支持向量回歸機(jī); 衰退預(yù)測

    中圖分類號:TH165+.3; TN911.7; TH133.33 文獻(xiàn)標(biāo)志碼: A 文章編號: 1004-4523(2018)05-0892-10

    DOI:10.16385/j.cnki.issn.1004-4523.2018.05.020

    引 言

    隨著大型旋轉(zhuǎn)機(jī)械設(shè)備自動(dòng)化、精密化、復(fù)雜化程度越來越高,對設(shè)備的工作環(huán)境要求越加嚴(yán)格。而滾動(dòng)軸承作為重要的機(jī)械零部件之一,其健康狀態(tài)直接影響著整臺(tái)機(jī)械設(shè)備乃至整條生產(chǎn)線能否正常運(yùn)行。由于滾動(dòng)軸承的實(shí)際壽命離散程度大,軸承輕微損傷對于精密機(jī)床而言,可能會(huì)導(dǎo)致非常嚴(yán)重的事故。因此對滾動(dòng)軸承進(jìn)行健康狀態(tài)監(jiān)測和退化趨勢預(yù)估極為必要,而關(guān)鍵在于建立合適的退化指標(biāo)及預(yù)測模型。

    振動(dòng)信號特征參數(shù)的選取對軸承的運(yùn)行狀態(tài)而言非常重要,衰退性能指標(biāo)可分為單一的時(shí)域特征參數(shù)、頻域特征參數(shù)及時(shí)頻域特征參數(shù)。其中,時(shí)域特征參數(shù)包括均方根值、峰峰值、裕度和峭度等[1-2],頻域特征參數(shù)則有頻譜方差、方根幅值和頻譜平均值等[3-4],但都存在信息不全面的問題,對某些故障情況敏感性差。單一的時(shí)頻域信號特征衰退評估通常采用經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition, EMD)、小波包分解等評估法[5-8],但這些方法都難以準(zhǔn)確地評價(jià)軸承非線性、非穩(wěn)定信號的退化過程。多特征參數(shù)性能退化評估方法通過提取振動(dòng)信號中多個(gè)特征參數(shù)來反映機(jī)械設(shè)備的運(yùn)行狀態(tài),定性評價(jià)機(jī)械設(shè)備的健康狀況。Lei Y[9]提出了一種新的自適應(yīng)多核關(guān)聯(lián)向量機(jī)的機(jī)器預(yù)測方法??凳貜?qiáng)、王玉靜[10]等提出了一種混沌優(yōu)化果蠅算法與多核超球體支持向量機(jī)的滾動(dòng)軸承健康狀態(tài)定量評估方法。對于多特征參數(shù)的評估,普遍存在信息量大,抗干擾能力強(qiáng)等特點(diǎn)、但存在信息冗余問題。

    流形學(xué)習(xí)是包含線性和非線性的一類數(shù)據(jù)降維方法的統(tǒng)稱。常用的線性維數(shù)約簡方法如主成分分析、獨(dú)立成分分析、Fisher判別分析、局部線性投影等,當(dāng)面對實(shí)際存在的海量高維數(shù)據(jù)樣本容量時(shí),其非線性特征則使現(xiàn)有的線性擬合模型不再適用。由此,非線性流形學(xué)習(xí)方法的典型代表如等距離特征映射、局部線性嵌入、拉普拉斯特征映射等算法應(yīng)運(yùn)而生,實(shí)現(xiàn)高維數(shù)據(jù)的維數(shù)約簡和可視化[11-12]。其中,LLE算法能夠?qū)W習(xí)任意維的低維特征,在數(shù)據(jù)進(jìn)行平移、展縮和旋轉(zhuǎn)變換時(shí),始終保持重構(gòu)誤差值最小。

    綜上所述,針對滾動(dòng)軸承性能衰退的預(yù)測從這兩方面都進(jìn)行了一定程度的研究。單一的性能特征參數(shù)包含信息量少、抗干擾性差等缺點(diǎn),而多特征參數(shù)存在信息量大和敏感程度差等缺點(diǎn),不能全面地對軸承退化過程進(jìn)行評價(jià)。流形學(xué)習(xí)方法可以將高維特征集映射到低維空間,把大量特征信息重新組合排序。因此本文提出基于流形學(xué)習(xí)和M-KH-SVR的滾動(dòng)軸承衰退預(yù)測方法,避免了因特征參數(shù)少而導(dǎo)致的單一化和特征參數(shù)過多而產(chǎn)生的信息冗余問題,證明了該方法具有很好的泛化性,對滾動(dòng)軸承的衰退過程區(qū)分更為精細(xì),提高了運(yùn)算效率和預(yù)測精度。

    通過表2中數(shù)據(jù)可得出,交叉驗(yàn)證法的預(yù)測精度為93.7%,平均相對百分比誤差為7.52%,遺傳算法的SVR預(yù)測精度為98.6%,平均相對百分比誤差為5.74%。運(yùn)用粒子群算法的SVR預(yù)測精度為98.8%,平均相對百分比誤差為4.88%,但是采用磷蝦群算法的SVR預(yù)測精度可以高達(dá)99.3%,平均相對百分比誤差為3.45%,擬合程度高,搜索能力也更強(qiáng),證明了新方法的可行性。

    3.2 實(shí)例驗(yàn)證Ⅱ

    采用辛辛那提大學(xué)[14]的滾動(dòng)軸承全壽命周期的試驗(yàn)數(shù)據(jù)為基礎(chǔ)進(jìn)行實(shí)驗(yàn)驗(yàn)證。如圖6所示,型號為Rexford ZA-2115的軸承4個(gè)安裝在試驗(yàn)臺(tái)上,滾動(dòng)體個(gè)數(shù)為16,直徑為0.331 cm,節(jié)圓直徑為2.815 cm,接觸角為15.17°,轉(zhuǎn)速為2000 r/min,PCB 353B33加速度傳感器分別安裝在軸承座的徑向方向和軸向方向。共運(yùn)行了33天收集到2000組試驗(yàn)數(shù)據(jù),采樣頻率是20 kHz。在實(shí)驗(yàn)后期階段,軸承3出現(xiàn)內(nèi)圈故障。

    集軸承3的部分原始實(shí)驗(yàn)數(shù)據(jù)。通過表3得到的原始特征指標(biāo)共計(jì)12類,但是由于原始特征參數(shù)的每種指標(biāo)衡量的方式不同,量綱不統(tǒng)一。因此把各類指標(biāo)歸一化后,通過MCC特征評價(jià)指標(biāo)進(jìn)行評價(jià),得到各個(gè)指標(biāo)與原始信號的相關(guān)程度值。通常相關(guān)度大于0.5以上,認(rèn)為關(guān)聯(lián)程度較高,因此篩選得到的特征指標(biāo)有均方根、峭度指標(biāo)、方根幅值、IMF能量熵P1~P4,共計(jì)7個(gè)指標(biāo)。采用實(shí)驗(yàn)采集到的滾動(dòng)軸承內(nèi)圈故障振動(dòng)信號,對比幾種不同預(yù)測方法。

    根據(jù)第1至32天的測試數(shù)據(jù),以此為預(yù)測模型的訓(xùn)練輸入值,而預(yù)測模型預(yù)測值是第2至33天的試驗(yàn)數(shù)據(jù)。通過比較預(yù)測數(shù)據(jù)與試驗(yàn)數(shù)據(jù)之間的差值得出預(yù)測結(jié)果的準(zhǔn)確與否。將預(yù)測輸入值各自分為單一特征指標(biāo)、多特征指標(biāo)、降維后的主分量指標(biāo),并對此三種不同的預(yù)測輸入值進(jìn)行評估實(shí)驗(yàn),具體實(shí)施步驟如下:

    (1)將滾動(dòng)軸承內(nèi)圈提取的均方根特征值,直接當(dāng)作預(yù)測的輸入。

    (2)將滾動(dòng)軸承內(nèi)圈提取18個(gè)故障特征值作為預(yù)測的輸入。

    (3)將滾動(dòng)軸承內(nèi)圈提取的故障特征值通過MCC特征評價(jià)后得到7個(gè)特征指標(biāo),經(jīng)LLE降維后得到的前3個(gè)特征分量作為預(yù)測的輸入。

    對比圖7(a)和7(b),可知由RMS預(yù)測的數(shù)據(jù)趨勢大體與原始數(shù)據(jù)相吻合,但由均方根自身特性所導(dǎo)致的,數(shù)據(jù)點(diǎn)在26.5天左右表現(xiàn)得較為平緩,波動(dòng)現(xiàn)象并不明顯,在26.5天到32.5天之間呈現(xiàn)增長緩慢的狀態(tài),而在32.5天以后開始表現(xiàn)為急劇增長的趨勢直至最后失效,可以看出單一的特征信號對滾動(dòng)軸承的退化趨勢無法起到準(zhǔn)確的預(yù)測效果。對比圖7(c)和7(d),數(shù)據(jù)點(diǎn)在26.5天與前者相似,可在26.5天到32.5天有逐漸增長的趨勢,波動(dòng)現(xiàn)象較為明顯,在32.5天以后呈現(xiàn)快速升高。因此,多特征值預(yù)測相較于單特征預(yù)測在整個(gè)衰退過程中表現(xiàn)更加明顯,存在的問題是在退化階段出現(xiàn)的波動(dòng)峰值預(yù)測并不準(zhǔn)確。對比圖7(e)和7(f), LLE降維預(yù)測特征值相較于前兩種方法,對滾動(dòng)軸承的全生命周期的退化趨勢表現(xiàn)更好,特別是在26.5天以后,滾動(dòng)軸承出現(xiàn)初期退化的趨勢,而在26.5到28天左右出現(xiàn)第一次衰退峰值,28到32天左右出現(xiàn)多次峰值,這是明顯的衰退加劇的表現(xiàn),32到33天退化趨勢急劇上升,一直到33天最后時(shí)刻滾動(dòng)軸承內(nèi)圈發(fā)生失效。

    根據(jù)表4數(shù)據(jù)可得,單特征參數(shù)相比多特征參數(shù)預(yù)測精度低,平均相對百分比誤差為9.4%,相比LLE特征參數(shù)作為輸入時(shí)高2倍有余。多特征參數(shù)雖說在精度上比單特征參數(shù)有略微上升但耗時(shí)長,平均相對百分比誤差為8.6%。LLE特征參數(shù)作為輸入時(shí)平均相對百分比誤差為4.2%,預(yù)測精度最高,計(jì)算時(shí)間相對較短。

    實(shí)驗(yàn)結(jié)果分析表明,LLE降維特征參數(shù)對滾動(dòng)軸承的各個(gè)退化階的段劃分更加細(xì)致,證明了LLE降維方法在衰退類別區(qū)分方面有很好的分類能力。磷蝦群算法對SVR中的參數(shù)進(jìn)行最優(yōu)選取,通過試驗(yàn)證明本文方法收斂速度快,預(yù)測精度高,驗(yàn)證了本文方法的實(shí)際應(yīng)用價(jià)值。

    4 結(jié) 論

    滾動(dòng)軸承的衰退趨勢預(yù)測對設(shè)備健康狀態(tài)起著重要的作用。本文提出了一種基于流形學(xué)習(xí)和M-KH-SVR的滾動(dòng)軸承衰退預(yù)測,避免了單一特征參數(shù)導(dǎo)致的信息缺失、故障敏感性差和多特征參數(shù)造成信息冗余的問題,通過評價(jià)篩選相關(guān)特征進(jìn)行局部線性嵌入降維,從而在保留原始特征的基礎(chǔ)上融合了更加全面的特征信息。并采用磷蝦群算法對多變量SVR預(yù)測模型中的c值和σ值進(jìn)行優(yōu)化,保證局部與全局最優(yōu),計(jì)算迭代速度更快。在提高模型預(yù)測精度和運(yùn)算效率的基礎(chǔ)上,對滾動(dòng)軸承的衰退過程劃分更為細(xì)致準(zhǔn)確,有效實(shí)現(xiàn)滾動(dòng)軸承的退化趨勢預(yù)測。

    參考文獻(xiàn):

    [1] Gu X, Yang S, Liu Y, et al. Rolling element bearing faults diagnosis based on kurtogram and frequency domain correlated kurtosis[J]. Measurement Science & Technology, 2016, 27(12):1—15.

    [2] Gangsar P, Tiwari R. Multiclass fault taxonomy in rolling bearings at interpolated and extrapolated speeds based on time domain bibration data by SVM algorithms[J]. Journal of Failure Analysis & Prevention, 2014, 14(6):826—837.

    [3] Xue X, Zhou J, Xu Y, et al. An adaptively fast ensemble empirical mode decomposition method and its applications to rolling element bearing fault diagnosis[J]. Mechanical Systems & Signal Processing, 2015, 62-63:444—459.

    [4] 鄭近德, 程軍圣. 改進(jìn)的希爾伯特-黃變換及其在滾動(dòng)軸承故障診斷中的應(yīng)用[J]. 機(jī)械工程學(xué)報(bào), 2015, 51(1):138—145.

    ZHENG Jinde, ChENG Jansheng. Improved Hilbert-Huang transform and its applications to rolling bearing fault diagnosis[J]. Journal of Mechanical Engineering, 2015, 51(1):138—145.

    [5] 王亞萍, 許 迪, 葛江華, 等. 基于SPWVD時(shí)頻圖紋理特征的滾動(dòng)軸承故障診斷[J]. 振動(dòng).測試與診斷, 2017, (1):115—119.

    WANG Yaping, XU Di, GE Jianhua, et al.Rolling bearing faults diagnostics based on SPWVD time-frequency distribution image texture feature[J]. Journal of Vibration, Measurement & Diagnosis, 2017, (1):115—119.

    [6] Alam S M S, Bhuiyan M I H. Detection of seizure and epilepsy using higher order statistics in the EMD domain[J]. IEEE Journal of Biomedical & Health Informatics, 2013, 17(2):312—318.

    [7] 孔德同, 劉慶超, 雷亞國,等. 一種改進(jìn)的EEMD方法及其應(yīng)用研究[J]. 振動(dòng)工程學(xué)報(bào), 2015, 28(6):1015—1021.

    Kong D T, Liu Q C, Lei Y G, et al. The improved EEMD method and its application[J]. Journal of Vibration Engineering, 2015, 28(6):1015—1021.

    [8] Yeh J R, Shieh J S, Huang N E. Complementary ensemble empirical mode decomposition: a novel noise enhanced data analysis method[J]. Advances in Adaptive Data Analysis, 2010 (2):135—156.

    [9] 雷亞國,陳 吳,李乃鵬,等.自適應(yīng)多核組合相關(guān)向量機(jī)預(yù)測方法及其在機(jī)械設(shè)備剩余壽命預(yù)測中的應(yīng)用[J].機(jī)械工程學(xué)報(bào),2016,52(1):87—93.

    Lei YYaoguo, Chen Wu, Li Naipeng, et al. A Relevance vector machine prediction method based on adaptive multi-kernel combination and its application to remaining useful life prediction of machinery[J]. Journal of Mechanical Engineering, 2016, 52(1):87—93.

    [10]康守強(qiáng), 王玉靜, 崔歷歷,等. 基于CFOA-MKHSVM的滾動(dòng)軸承健康狀態(tài)評估方法[J]. 儀器儀表學(xué)報(bào), 2016, 37(9):2029—2035.

    Kang Shouqiang, Wang Yujing, Cui Lili, et al.Health state assessment of a rolling bearing based on CFOA-MKHSVM method[J]. Chinese Journal of Scientific Instrument, 2016, 37(9):2029—2035.

    [11]肖 婷, 湯寶平, 秦 毅,等. 基于流形學(xué)習(xí)和最小二乘支持向量機(jī)的滾動(dòng)軸承退化趨勢預(yù)測[J]. 振動(dòng)與沖擊, 2015, 34(9):149—153.

    Xiao T, Tang B P, Qin Y, et al. Degradation trend prediction of rolling bearing based on manifold learning and least squares support vector machine[J]. Journal of Vibration & Shock, 2015, 34(9):149—153.

    [12]謝小欣, 胡建中, 許飛云,等. 一種基于多流形局部線性嵌入算法的故障診斷方法[J]. 機(jī)械工程學(xué)報(bào), 2013, 49(11):79—83.

    XIE Xiaoxin, HU Jianzhong, XU Feiyun, et al. A fault diagnosis method using multi-manifold learning based on locally linear embedding[J]. Journal of Mechanical Engineering, 2013, 49(11):79—83.

    [13]Loparo K A. Case Western Reserve University Bearing Data Center. [EB/OL]:[2017-09-10].http://www.eecs. cwru.edu/ laboratory/ bearing/download.

    [14]Qiu H, Lee J, Lin J, et al. Wavelet filter-based weak signature detection method and its application on rolling element bearing prognostics[J].Journal of Sound & Vibration, 2006, 289(4):1066-1090.

    Abstract: In order to deal with the big data sets and non-stationary signal of the rolling bearing, a degradation prediction method for rolling bearings is proposed combining manifold learning with multivariable support vector regression based on krill herd (M-KH-SVR). Firstly, the time-frequency domain features of the rolling bearing are extracted to construct the initial feature sets. Then the features are selected by multiple correlation coefficients to get the feature sets with high correlation, and the local linear embedding method is used to reduce the feature dimension and construct a new set of fault feature sets. Finally, krill herd algorithm is used to optimize parameter and in the multivariable support vector regression. By using its ability of local optimization and global optimization, the efficiency of parameter selection is improved. Compared with the single parameter and multi-features parameter method, the effectiveness of the proposed index is verified, and itscomputational efficiency and prediction accuracy is improved. It is shown that the division of the rolling bearing recession process is more accurate.

    Key words: fault diagnosis; rolling bearing; feature extraction; krill herd; support vector regression; performance degradation

    作者簡介:許 迪(1991—),女,博士研究生。電話:15114642186;E-mail:15114642186@163.com

    通訊作者:葛江華(1963—),女,教授。電話:13804550139;E-mail:gejianghua@sina.com

    猜你喜歡
    滾動(dòng)軸承特征提取故障診斷
    基于Gazebo仿真環(huán)境的ORB特征提取與比對的研究
    電子制作(2019年15期)2019-08-27 01:12:00
    一種基于LBP 特征提取和稀疏表示的肝病識(shí)別算法
    滾動(dòng)軸承裝配方法及注意事項(xiàng)
    滾動(dòng)軸承安裝使用中的幾個(gè)問題
    基于小波包絡(luò)分析的滾動(dòng)軸承典型故障診斷技術(shù)研究
    商情(2016年11期)2016-04-15 07:39:54
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    基于包絡(luò)分析法的軸承故障分析
    基于MED和循環(huán)域解調(diào)的多故障特征提取
    基于LCD和排列熵的滾動(dòng)軸承故障診斷
    基于WPD-HHT的滾動(dòng)軸承故障診斷
    我的女老师完整版在线观看| 丝袜喷水一区| 免费看日本二区| 三级国产精品片| 少妇 在线观看| 亚洲精品456在线播放app| 男人狂女人下面高潮的视频| 日本黄色片子视频| 中文字幕制服av| 成人美女网站在线观看视频| av在线app专区| 国产精品麻豆人妻色哟哟久久| 精品一区二区三卡| 午夜视频国产福利| 国产精品一区二区三区四区免费观看| 国产美女午夜福利| 久久精品国产亚洲网站| 国产v大片淫在线免费观看| 色哟哟·www| 国产黄色免费在线视频| 国产精品国产三级专区第一集| 91久久精品电影网| 欧美日韩综合久久久久久| 免费观看在线日韩| 免费高清在线观看视频在线观看| 最新中文字幕久久久久| 久久久国产一区二区| 久久综合国产亚洲精品| 久久 成人 亚洲| 在线 av 中文字幕| 超碰97精品在线观看| 内地一区二区视频在线| av专区在线播放| 伦理电影大哥的女人| 天堂中文最新版在线下载| 日产精品乱码卡一卡2卡三| 99热网站在线观看| 久久久精品94久久精品| 亚洲成人一二三区av| 在线播放无遮挡| 亚洲av中文av极速乱| 黑人猛操日本美女一级片| 亚洲高清免费不卡视频| 99热这里只有是精品在线观看| 久久久久久久亚洲中文字幕| 国产伦理片在线播放av一区| 亚洲国产高清在线一区二区三| 一级毛片黄色毛片免费观看视频| 亚洲精品日韩在线中文字幕| 最黄视频免费看| 大片电影免费在线观看免费| 黄色日韩在线| 亚洲av欧美aⅴ国产| 美女主播在线视频| 欧美精品国产亚洲| 国产精品一区二区在线不卡| 能在线免费看毛片的网站| 97精品久久久久久久久久精品| 高清毛片免费看| 国产一区有黄有色的免费视频| 少妇精品久久久久久久| 国产av码专区亚洲av| 久久久久精品久久久久真实原创| 天天躁日日操中文字幕| 亚洲图色成人| 大话2 男鬼变身卡| 舔av片在线| 在线免费观看不下载黄p国产| 夜夜爽夜夜爽视频| 肉色欧美久久久久久久蜜桃| 99久久精品热视频| 国产精品久久久久久av不卡| 国产亚洲av片在线观看秒播厂| av福利片在线观看| 亚洲精品乱久久久久久| 国产中年淑女户外野战色| 人妻 亚洲 视频| 午夜福利高清视频| 日韩强制内射视频| 免费看日本二区| 免费大片黄手机在线观看| 又大又黄又爽视频免费| 久久午夜福利片| 熟妇人妻不卡中文字幕| 熟女电影av网| 免费久久久久久久精品成人欧美视频 | 欧美日韩视频精品一区| 美女内射精品一级片tv| 大又大粗又爽又黄少妇毛片口| 边亲边吃奶的免费视频| 天堂8中文在线网| 国产精品久久久久久久久免| 亚洲经典国产精华液单| 久久人人爽人人爽人人片va| 亚洲av综合色区一区| 久久鲁丝午夜福利片| 日韩伦理黄色片| 国产成人a∨麻豆精品| 伊人久久国产一区二区| 久久久亚洲精品成人影院| 久久影院123| 秋霞在线观看毛片| 黄色日韩在线| 亚洲精品国产成人久久av| 啦啦啦视频在线资源免费观看| 丝瓜视频免费看黄片| 久久国产精品大桥未久av | 成人无遮挡网站| 国产 一区 欧美 日韩| 99精国产麻豆久久婷婷| 97在线视频观看| 亚洲精品成人av观看孕妇| 人人妻人人爽人人添夜夜欢视频 | 日韩大片免费观看网站| 秋霞在线观看毛片| av福利片在线观看| 国产精品三级大全| 干丝袜人妻中文字幕| 亚洲国产日韩一区二区| 精华霜和精华液先用哪个| 夫妻午夜视频| 97热精品久久久久久| 亚洲图色成人| 精品酒店卫生间| 日本与韩国留学比较| 免费久久久久久久精品成人欧美视频 | 精品酒店卫生间| 亚洲欧美一区二区三区黑人 | 天堂俺去俺来也www色官网| 亚洲成人一二三区av| 国产久久久一区二区三区| 亚洲国产精品国产精品| 久久精品国产亚洲av涩爱| 国产免费一级a男人的天堂| 精品国产三级普通话版| 亚洲国产精品一区三区| 简卡轻食公司| 国产日韩欧美亚洲二区| 国产一区二区三区综合在线观看 | 色哟哟·www| 久久国内精品自在自线图片| 国产伦精品一区二区三区四那| 亚洲欧美精品专区久久| 插阴视频在线观看视频| 久热这里只有精品99| 久久99精品国语久久久| 插阴视频在线观看视频| 国产成人一区二区在线| 亚洲自偷自拍三级| 中文字幕亚洲精品专区| 午夜福利高清视频| av在线播放精品| 又爽又黄a免费视频| 在线看a的网站| 成人漫画全彩无遮挡| 欧美极品一区二区三区四区| 亚洲电影在线观看av| 精品一区在线观看国产| 熟妇人妻不卡中文字幕| 亚洲av国产av综合av卡| 高清不卡的av网站| 青春草视频在线免费观看| 亚洲真实伦在线观看| 波野结衣二区三区在线| 夫妻午夜视频| av线在线观看网站| 午夜福利网站1000一区二区三区| 日日摸夜夜添夜夜爱| 亚洲国产毛片av蜜桃av| 日韩av在线免费看完整版不卡| 在线观看美女被高潮喷水网站| 99热这里只有精品一区| 2022亚洲国产成人精品| 我要看黄色一级片免费的| 赤兔流量卡办理| 在线观看一区二区三区| 欧美精品亚洲一区二区| 国产av精品麻豆| 男女国产视频网站| 蜜桃久久精品国产亚洲av| 亚洲国产高清在线一区二区三| a级毛片免费高清观看在线播放| 欧美 日韩 精品 国产| 夜夜看夜夜爽夜夜摸| 日韩强制内射视频| 国产日韩欧美在线精品| 精华霜和精华液先用哪个| 国产亚洲精品久久久com| 熟女av电影| 国内揄拍国产精品人妻在线| 久久久久性生活片| 国产精品国产av在线观看| 国产av精品麻豆| 国产亚洲欧美精品永久| 亚洲图色成人| 国产午夜精品久久久久久一区二区三区| 天美传媒精品一区二区| 国产v大片淫在线免费观看| 看十八女毛片水多多多| 国产亚洲欧美精品永久| 熟女人妻精品中文字幕| 免费大片黄手机在线观看| 黑人高潮一二区| 国产69精品久久久久777片| 国产黄片美女视频| 精品一区二区免费观看| 欧美zozozo另类| 99热这里只有是精品在线观看| h视频一区二区三区| 午夜福利网站1000一区二区三区| 亚洲自偷自拍三级| 婷婷色综合www| av黄色大香蕉| 国产成人精品福利久久| 午夜福利影视在线免费观看| 极品少妇高潮喷水抽搐| 黑丝袜美女国产一区| 一本色道久久久久久精品综合| 日本一二三区视频观看| 另类亚洲欧美激情| 亚洲欧美成人综合另类久久久| 五月天丁香电影| 91狼人影院| 国产精品国产三级国产av玫瑰| 99热全是精品| 全区人妻精品视频| 18禁在线播放成人免费| 夜夜爽夜夜爽视频| 国产成人91sexporn| 国产老妇伦熟女老妇高清| 午夜免费观看性视频| 美女主播在线视频| 简卡轻食公司| 99视频精品全部免费 在线| 亚洲av电影在线观看一区二区三区| av线在线观看网站| 哪个播放器可以免费观看大片| 美女cb高潮喷水在线观看| 亚洲性久久影院| 国产一区亚洲一区在线观看| 欧美成人精品欧美一级黄| 欧美xxⅹ黑人| 大片免费播放器 马上看| 亚洲三级黄色毛片| 最黄视频免费看| 人妻少妇偷人精品九色| 欧美三级亚洲精品| 高清不卡的av网站| 欧美激情国产日韩精品一区| 欧美国产精品一级二级三级 | 精品酒店卫生间| 一级a做视频免费观看| 如何舔出高潮| 亚洲av中文av极速乱| 国模一区二区三区四区视频| 九草在线视频观看| 亚洲国产精品一区三区| 80岁老熟妇乱子伦牲交| 97超视频在线观看视频| 久久精品国产亚洲av天美| 亚洲激情五月婷婷啪啪| 国产片特级美女逼逼视频| 亚洲中文av在线| 男女国产视频网站| 老女人水多毛片| 1000部很黄的大片| 各种免费的搞黄视频| 国产真实伦视频高清在线观看| 国产免费一级a男人的天堂| 高清av免费在线| 久久久色成人| 一级爰片在线观看| 国内少妇人妻偷人精品xxx网站| 在线观看一区二区三区| 美女cb高潮喷水在线观看| 欧美变态另类bdsm刘玥| 国产乱人视频| 国产av国产精品国产| 99热这里只有是精品50| 深夜a级毛片| videossex国产| 亚洲av中文字字幕乱码综合| 亚洲天堂av无毛| 欧美成人一区二区免费高清观看| 另类亚洲欧美激情| 草草在线视频免费看| 人妻少妇偷人精品九色| 日本黄色日本黄色录像| 男人舔奶头视频| 搡老乐熟女国产| 又爽又黄a免费视频| 中文字幕av成人在线电影| 一二三四中文在线观看免费高清| 久久久久久久国产电影| 欧美日韩一区二区视频在线观看视频在线| 91久久精品国产一区二区三区| 欧美日韩亚洲高清精品| 夫妻性生交免费视频一级片| 一级a做视频免费观看| 18禁动态无遮挡网站| 丝瓜视频免费看黄片| 婷婷色综合大香蕉| 亚洲av免费高清在线观看| 国产又色又爽无遮挡免| 少妇丰满av| 亚洲精品国产av成人精品| 另类亚洲欧美激情| 草草在线视频免费看| 18禁裸乳无遮挡免费网站照片| 高清在线视频一区二区三区| 亚洲欧美日韩东京热| av在线观看视频网站免费| 美女xxoo啪啪120秒动态图| 国产91av在线免费观看| 免费观看av网站的网址| 夫妻午夜视频| 一级黄片播放器| 国产黄片美女视频| 精品午夜福利在线看| 亚洲国产色片| 人妻制服诱惑在线中文字幕| 丰满人妻一区二区三区视频av| 99热6这里只有精品| 人人妻人人看人人澡| 又大又黄又爽视频免费| 97热精品久久久久久| 18禁在线播放成人免费| 大片免费播放器 马上看| 男女无遮挡免费网站观看| 中文欧美无线码| 欧美xxxx黑人xx丫x性爽| 日本欧美视频一区| 男女免费视频国产| 国产精品爽爽va在线观看网站| 91精品国产九色| 欧美日韩综合久久久久久| 成人国产av品久久久| av又黄又爽大尺度在线免费看| 麻豆精品久久久久久蜜桃| 国产av精品麻豆| 黑人猛操日本美女一级片| 直男gayav资源| 日日啪夜夜撸| videossex国产| 国产成人一区二区在线| 各种免费的搞黄视频| 国产精品久久久久久久久免| 欧美区成人在线视频| 一边亲一边摸免费视频| 97精品久久久久久久久久精品| 欧美人与善性xxx| 亚洲av.av天堂| 亚洲精品第二区| 久久国产精品大桥未久av | 成人漫画全彩无遮挡| 欧美成人午夜免费资源| 精品国产露脸久久av麻豆| 日本猛色少妇xxxxx猛交久久| 性色av一级| 国产有黄有色有爽视频| 黑人猛操日本美女一级片| 香蕉精品网在线| 夜夜看夜夜爽夜夜摸| 深爱激情五月婷婷| 亚洲精品一区蜜桃| 国产永久视频网站| 在线观看免费高清a一片| 日韩视频在线欧美| 在线观看免费日韩欧美大片 | 国产成人精品婷婷| 亚洲人成网站高清观看| 国产男人的电影天堂91| 国产日韩欧美在线精品| 色哟哟·www| 免费人成在线观看视频色| 国产男人的电影天堂91| 国产淫语在线视频| 精品一区在线观看国产| 国产精品不卡视频一区二区| 国产综合精华液| 中文资源天堂在线| 久久精品久久精品一区二区三区| www.av在线官网国产| 久久99热这里只频精品6学生| 尾随美女入室| 午夜免费鲁丝| 日韩,欧美,国产一区二区三区| av在线播放精品| 十八禁网站网址无遮挡 | 亚洲国产成人一精品久久久| 日韩精品有码人妻一区| 国产久久久一区二区三区| 亚洲伊人久久精品综合| 久久久a久久爽久久v久久| 免费不卡的大黄色大毛片视频在线观看| 内地一区二区视频在线| 亚洲欧美成人综合另类久久久| av线在线观看网站| 色吧在线观看| 精品视频人人做人人爽| 久久久精品94久久精品| 人人妻人人爽人人添夜夜欢视频 | 91精品伊人久久大香线蕉| 看免费成人av毛片| 少妇猛男粗大的猛烈进出视频| 欧美日韩在线观看h| 我要看黄色一级片免费的| 久久久久久久久久人人人人人人| 综合色丁香网| 小蜜桃在线观看免费完整版高清| 国产亚洲欧美精品永久| 中文在线观看免费www的网站| 亚洲欧洲国产日韩| 蜜桃亚洲精品一区二区三区| 免费人妻精品一区二区三区视频| 国产中年淑女户外野战色| 亚洲欧美精品自产自拍| 亚洲综合精品二区| 亚洲色图av天堂| 人人妻人人添人人爽欧美一区卜 | 少妇人妻 视频| 超碰97精品在线观看| 91久久精品电影网| 身体一侧抽搐| 国精品久久久久久国模美| 美女中出高潮动态图| 亚洲美女视频黄频| 99热6这里只有精品| 亚洲第一区二区三区不卡| 亚洲成人手机| 欧美老熟妇乱子伦牲交| 国产男人的电影天堂91| 建设人人有责人人尽责人人享有的 | 亚洲欧美日韩东京热| 欧美国产精品一级二级三级 | 亚洲av电影在线观看一区二区三区| 尾随美女入室| 啦啦啦视频在线资源免费观看| 亚洲成色77777| 97热精品久久久久久| 一级毛片黄色毛片免费观看视频| 下体分泌物呈黄色| 免费人妻精品一区二区三区视频| 欧美 日韩 精品 国产| 欧美bdsm另类| 最近最新中文字幕大全电影3| 2018国产大陆天天弄谢| 国产精品不卡视频一区二区| 一二三四中文在线观看免费高清| 大陆偷拍与自拍| 成年人午夜在线观看视频| 国产乱人视频| 嘟嘟电影网在线观看| 久久精品久久久久久噜噜老黄| 国产日韩欧美亚洲二区| 色哟哟·www| 欧美老熟妇乱子伦牲交| 美女福利国产在线 | 久久 成人 亚洲| 欧美激情极品国产一区二区三区 | a级毛色黄片| 欧美精品一区二区大全| 国内精品宾馆在线| 人体艺术视频欧美日本| 午夜福利视频精品| 日韩国内少妇激情av| 91久久精品国产一区二区三区| 国产成人a区在线观看| 国产精品av视频在线免费观看| 国产乱来视频区| 在线观看av片永久免费下载| 五月天丁香电影| 欧美区成人在线视频| 久久韩国三级中文字幕| 十八禁网站网址无遮挡 | 亚洲,一卡二卡三卡| 欧美高清成人免费视频www| 新久久久久国产一级毛片| 国产极品天堂在线| 在线亚洲精品国产二区图片欧美 | 午夜福利在线在线| 国产精品蜜桃在线观看| 亚洲国产av新网站| 精品视频人人做人人爽| 18禁在线播放成人免费| .国产精品久久| 国产精品熟女久久久久浪| 婷婷色麻豆天堂久久| 国产精品久久久久久精品古装| 亚洲国产精品专区欧美| 久久久久久久亚洲中文字幕| 亚洲国产精品一区三区| 国产精品.久久久| av不卡在线播放| 亚洲一区二区三区欧美精品| 精品99又大又爽又粗少妇毛片| 日日撸夜夜添| 18+在线观看网站| 国产黄片视频在线免费观看| 一级毛片黄色毛片免费观看视频| 欧美xxxx性猛交bbbb| 日本午夜av视频| 中文字幕免费在线视频6| 亚洲一级一片aⅴ在线观看| 一级毛片aaaaaa免费看小| 欧美日韩视频高清一区二区三区二| 网址你懂的国产日韩在线| 久久久久网色| 日本与韩国留学比较| 高清在线视频一区二区三区| 精品久久久久久久久av| 亚洲最大成人中文| 一区二区三区乱码不卡18| 中文天堂在线官网| 亚洲av不卡在线观看| 美女国产视频在线观看| 亚洲人成网站在线观看播放| 国产高清有码在线观看视频| 国产精品一区二区性色av| 人妻系列 视频| 成人国产麻豆网| 中文字幕久久专区| 人人妻人人添人人爽欧美一区卜 | 亚洲色图综合在线观看| 夜夜看夜夜爽夜夜摸| 亚洲欧美日韩另类电影网站 | 超碰97精品在线观看| 亚洲综合色惰| av视频免费观看在线观看| 卡戴珊不雅视频在线播放| av在线蜜桃| 免费看日本二区| 国产精品国产三级国产av玫瑰| 中文资源天堂在线| av国产久精品久网站免费入址| 男女边摸边吃奶| 嘟嘟电影网在线观看| 亚洲自偷自拍三级| 不卡视频在线观看欧美| 亚洲四区av| 久久国产乱子免费精品| 国产精品麻豆人妻色哟哟久久| 国产精品爽爽va在线观看网站| 成人免费观看视频高清| 亚洲欧美中文字幕日韩二区| 国产 一区精品| 最近最新中文字幕大全电影3| 日韩强制内射视频| 黑丝袜美女国产一区| 色婷婷久久久亚洲欧美| 精品人妻视频免费看| 欧美丝袜亚洲另类| 亚洲av男天堂| 91狼人影院| 午夜福利影视在线免费观看| 亚洲欧美精品专区久久| tube8黄色片| 国产精品精品国产色婷婷| 国产有黄有色有爽视频| 99热这里只有是精品在线观看| 高清日韩中文字幕在线| 中国美白少妇内射xxxbb| 人妻一区二区av| 麻豆国产97在线/欧美| 女性生殖器流出的白浆| 精品一区二区三区视频在线| 少妇丰满av| av专区在线播放| 亚洲高清免费不卡视频| 日韩av不卡免费在线播放| av线在线观看网站| 91久久精品国产一区二区三区| 大陆偷拍与自拍| 久久精品国产亚洲网站| 一本—道久久a久久精品蜜桃钙片| 人妻一区二区av| 婷婷色综合www| 哪个播放器可以免费观看大片| 国产色婷婷99| 久久国产亚洲av麻豆专区| 久久热精品热| av线在线观看网站| 日韩强制内射视频| 免费观看av网站的网址| 女的被弄到高潮叫床怎么办| 有码 亚洲区| 欧美最新免费一区二区三区| 乱系列少妇在线播放| www.色视频.com| 熟女av电影| 高清欧美精品videossex| 日本色播在线视频| 国产色婷婷99| 美女内射精品一级片tv| 这个男人来自地球电影免费观看 | 极品教师在线视频| 午夜激情久久久久久久| 国产白丝娇喘喷水9色精品| 少妇人妻精品综合一区二区| av免费在线看不卡| 在线精品无人区一区二区三 | 国产 一区精品| 中文字幕久久专区| 大片免费播放器 马上看| 日日啪夜夜撸| 亚洲第一av免费看| 女人十人毛片免费观看3o分钟| 男女边摸边吃奶| 一个人免费看片子| 久久99热这里只频精品6学生| 国产永久视频网站| 又爽又黄a免费视频| 亚洲精品久久午夜乱码| 亚洲av欧美aⅴ国产| 亚洲av在线观看美女高潮|