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

    基于主元分析的邏輯回歸模型混合域協(xié)變量選取研究

    2018-07-26 08:28:56王萌王奉濤
    軸承 2018年12期
    關鍵詞:主元頻域時域

    王萌,王奉濤

    (大連理工大學 機械工程學院,遼寧 大連 116024)

    目前,滾動軸承壽命理論歸納為基于統(tǒng)計分析的壽命模型,基于斷裂力學分析的壽命模型以及基于狀態(tài)監(jiān)測的壽命模型,其中基于狀態(tài)監(jiān)測的壽命預測方法是當前研究的重點方向[1-3]。

    文獻[4]提出了基于邏輯回歸的性能退化評估模型,并用試驗進行了驗證。文獻[5]以均方根值和峭度值為協(xié)變量,比較了邏輯回歸模型和Weibull比例故障模型對軸承退化狀態(tài)的表征。文獻[6]將邏輯回歸模型應用在刀具的壽命預測中,其建立的刀具壽命預測模型的協(xié)變量選取的參數(shù)是小波包能量、小波包能量熵和時域特征。邏輯回歸模型與Weibull比例故障模型的主要區(qū)別是,邏輯回歸模型的假設條件和估計的參數(shù)較少,但能夠更精確的預測剩余壽命。

    在以往的研究中,邏輯回歸模型使用的協(xié)變量選取不一,通常選取幾個時域特征值作為協(xié)變量。然而,單一的評價指標并不能完全的表征軸承的退化過程。因此,提出一種基于PCA的邏輯回歸模型協(xié)變量選取方法,其可以在保留大多數(shù)特征值的情況下更好的表征軸承的退化狀態(tài),從而更加有效地對軸承的可靠性進行評估。

    1 理論基礎

    1.1 主元分析

    PCA是一種常用的數(shù)據(jù)分析方法,可以提取數(shù)據(jù)的主要特征分量,常用于高維數(shù)據(jù)的降維[7-8]。

    設有m個變量x1,x2,…,xm,每個變量有n個樣本,第i個變量可以表示為xi=(x1i,x2i,…,xni)T,則由此構成的n×m階矩陣為

    (1)

    將X記為X=[x1,x2,…,xm],其具體算法如下:

    1)對原始的特征矩陣X進行特征中心化,計算均值,即

    (2)

    2)計算特征矢量的協(xié)方差矩陣

    (3)

    3)計算C的特征值λi和特征向量vi

    Cvi=λivi;i=1,2,…,n。

    (4)

    4)計算貢獻率

    (5)

    5)將λi從大到小排列,取前k個特征值與其對應的特征向量組成Δ=(λ1,λ2,…,λk)和V=(v1,v2,…,vk)。

    6)計算k維特征矢量(k

    P=VTX。

    (6)

    1.2 邏輯回歸模型

    邏輯回歸模型中,因變量yt的取值為0和1,分別表示事件僅有的2種獨立情況;t為事件發(fā)生的時間坐標,t=1,2,…,n;協(xié)變量X(t)={x1(t),x2(t),x3(t),…,xm(t)},m為協(xié)變量個數(shù),即特征指標數(shù)。則事件不發(fā)生(yt=1)的條件概率為

    P(yt=1|X(t))=

    (7)

    式中:β0為截距(或稱常數(shù)項);βj為協(xié)變量xj(t)對應的回歸系數(shù)。

    滾動軸承當前的特征量為X(t),則滾動軸承的可靠度函數(shù)R(t|X(t))與累積失效分布函數(shù)F(t|X(t))=1-R(t|X(t))之比滿足

    exp[β0+β1x1(t)+β2x2(t)+…+βmxm(t)]。

    (8)

    由極大似然估計方法求出β0,β1,β2,…,βm后,滾動軸承的可靠度函數(shù)可以表示為

    P(t|X(t))=

    (9)

    2 方法步驟

    依據(jù)上述理論基礎,所設計的算法流程如圖1所示,具體實現(xiàn)步驟為:

    圖1 本文算法流程圖Fig.1 Flow chart of algorithm in this paper

    1)特征參數(shù)選擇。從滾動軸承數(shù)據(jù)集中提取全壽命周期的時域特征、頻域特征和時頻域特征參數(shù)。

    2)時域特征處理。對時域特征參數(shù)進行相對特征處理,求取相對時域特征值并與相對均方根值和相對峭度值進行相關分析,求取相關系數(shù)大于0.8的時域特征值組成時域特征集。

    3)時頻域特征處理。對軸承數(shù)據(jù)進行3層小波包分解,求取頻帶能量最大子帶對應的樣本熵作為時頻域特征值。

    4)PCA降維。對時域、頻域和時頻域組成的混合域特征集進行PCA降維,將高維特征集降為低維特征集,作為邏輯回歸模型的協(xié)變量。

    5)軸承可靠性評估。將基于混合域的邏輯回歸模型的可靠性與基于時域相對均方根值和相對峭度值的邏輯回歸模型的軸承可靠性進行對比,驗證基于混合域特征集的邏輯回歸模型的有效性。

    3 試驗分析

    3.1 試驗平臺

    采用由美國辛辛那提大學智能維護中心(IMS)提供的滾動軸承全壽命周期加速軸承性能退化試驗數(shù)據(jù)[9]。試驗臺中的軸承位置和傳感器布置如圖2所示,試驗軸承為雙列圓錐滾子軸承,每列16個滾子,滾子組節(jié)圓直徑為71.501 mm,接觸角為15.17°,軸承轉(zhuǎn)速為2 000 r/min,由彈性加載器施加26.689 kN的徑向載荷。油反饋管道安裝有磁性螺塞來收集潤滑油中的碎屑,用以驗證軸承的性能退化。系統(tǒng)電器開關的關閉由磁性螺塞所吸附的金屬碎屑量決定,隨著軸承性能不斷退化,當吸附的碎屑量達到預先設定的閾值時數(shù)據(jù)采集工作便會停止。傳感器為8個高靈敏度石英加速度傳感器PCB353B33,分別采集每個軸承x和y方向的加速度信號,采集時間間隔為20 min,采樣頻率為20 kHz,采樣長度為20 480點。

    圖2 試驗臺示意圖Fig.2 Diagram of test rig

    共進行了3次試驗,每次試驗中4套軸承均同時運行,選取第1次試驗中3#軸承的全壽命數(shù)據(jù)進行分析。第1次試驗結(jié)束時3#軸承內(nèi)圈的嚴重故障如圖3所示。

    圖3 軸承內(nèi)圈的嚴重故障Fig.3 Serious fault of bearing inner ring

    3.2 協(xié)變量選取方法研究

    根據(jù)滾動軸承振動信號的特性,選擇能準確反映滾動軸承運行狀態(tài)的特征值。

    3.2.1 時域特征

    而由于制造、安裝和實際工況的差異,同一工作環(huán)境下同型號軸承間的特征參數(shù)也會存在一定的差異,8套軸承正常工作期內(nèi)一段趨勢平穩(wěn)振動信號的平均時域特征如圖4所示。從圖中可以看出,不同軸承的時域特征值存在較大的差異,為避免軸承個體差異的影響,兼顧良好的上升趨勢并對初始損傷保持敏感,采用相對特征值進行后續(xù)處理。具體方法為:選取每個時域特征值正常運行期間中一段趨勢平穩(wěn)的數(shù)據(jù),將該段數(shù)據(jù)的平均值設定為標準值,計算原始特征值與其對應的標準值之比,即相對特征值[11]。

    圖4 軸承正常運行期間的平均時域特征參數(shù)Fig.4 Mean time domain characteristic parameters of bearing during normal operation

    時域特征值中常用的指標是均方根值和峭度值,取相對均方根值和相對峭度值作為進行相關分析的依據(jù),將19個初選時域特征指標分別與相對均方根值和相對峭度值進行比較,取相關系數(shù)大于0.8的時域指標組成時域特征集,結(jié)果見表1。由表可知,相關系數(shù)大于0.8的時域特征值有13個,包括均方根值、標準差、方差、絕對均值、平均功率、方根幅值、峰值、峭度、峰值指標、脈沖指標、裕度指標、峭度指標和THIKAT。將這些時域指標作為邏輯回歸模型的協(xié)變量參數(shù)。

    表1 相關系數(shù)Tab.1 Correlation coefficients

    3.2.2 頻域特征

    頻域特征值選取重心頻率和均方頻率,其可以描述功率譜重心位置的變化情況;以及頻率方差,用于描述功率譜能量分布的分散程度[12-13]。

    3.2.3 時頻域特征

    小波包分解既對信號中的高頻成分擁有較好的時頻分辨率,又對低頻信號的頻率分辨率好。因此,選取3層小波包分解頻帶能量最大子帶對應的樣本熵作為時頻域特征值。

    3.2.4 主元分析

    將選擇好的時域、頻域和時頻域特征值組成混合域特征集,然后采用PCA對高維相對特征參數(shù)進行降維,前3個主元的累計貢獻率見表2。由表可知:高維相對特征集經(jīng)過PCA降維后前2個主元的累計貢獻率已經(jīng)達到了90%以上,第3主元的貢獻率僅為3.09%,因此選擇前2個主元作為邏輯回歸模型的協(xié)變量。

    表2 主元分析結(jié)果Tab.2 Results of principal component analysis

    為觀察主元分析的效果,將標準高維訓練集投影到二維空間上,結(jié)果如圖5所示。從圖中可以看出:正常、早期故障、中期故障、嚴重故障從小到大依次排列,明顯反應出軸承的退化趨勢,其中恢復期[15]在第1主元上處于早期故障之間,這是由于軸承早期故障裂紋被磨平,振動沖擊減小,故障特征頻率不明顯所導致;從第1主元上可以看出各個故障區(qū)分布很明顯;從第2主元上可以看出特別明顯的早期故障趨勢,這是由于峭度值、峰值等特征指標對早期故障特別敏感的原因。

    圖5 第1,2主元滾動軸承狀態(tài)圖Fig.5 State diagram of rolling bearing of first and second principal component

    3.3 軸承可靠性評估

    為驗證所述算法的有效性,將基于混合域特征集邏輯回歸模型的預測結(jié)果與基于時域相對均方根值和相對峭度值的邏輯回歸模型結(jié)果進行對比分析。2種模型獲取的可靠度曲線如圖6所示。

    圖6 邏輯回歸模型的可靠度曲線Fig.6 Reliability curve of logistic regression model

    對比分析可知:

    1)軸承正常運行時期,混合域邏輯回歸模型的可靠度在軸承從接近100%的位置開始,并隨著軸承使用時間的增加平緩下滑;而時域邏輯回歸模型的可靠度起點不僅在88%的位置,而且在軸承正常運行期間出現(xiàn)較大波動;說明經(jīng)過PCA降維得到的特征值對軸承正常時期的表征要比僅僅依靠時域特征得到的可靠度要更加穩(wěn)定、準確。

    2)當軸承出現(xiàn)早期故障時,基于混合域特征的可靠度曲線出現(xiàn)小幅波動,并開始緩慢下降;而基于時域特征的可靠度曲線出現(xiàn)了劇烈的波動,這是由于峭度值對軸承早期故障特別敏感。

    3)到達中期故障時,基于混合域特征的可靠度曲線波動仍較小,而基于時域特征的可靠度曲線仍存在劇烈的波動。

    4)后期故障時,基于混合域的可靠度曲線能夠平緩地描述軸承的退化過程,而基于時域的可靠度曲線出現(xiàn)斷崖式的曲線。

    綜上所述,基于混合域的可靠度曲線在表征軸承的退化過程中要比基于單一時域特征值的曲線更加穩(wěn)定、準確,且不受單一特征值對故障狀態(tài)敏感性的影響。

    4 結(jié)論

    運用主元分析對邏輯回歸模型的協(xié)變量選取進行了優(yōu)化,通過試驗研究得到如下結(jié)論:

    1)相對時域特征值降低了軸承制造、安裝和實際工況差異的影響,使軸承間的各體差異對時域特征集構建的影響大大降低。

    2)將時域、頻域和時頻域組成的混合域特征集經(jīng)過PCA降維后得到的特征值比單純的時域特征值蘊含了更多能夠表征軸承退化狀態(tài)的信息。

    3)基于混合域特征集的邏輯回歸模型能更好的表征軸承的退化性能,比基于時域的模型更加穩(wěn)定,可靠。

    猜你喜歡
    主元頻域時域
    多元并行 誰主沉浮
    應用主元變換法分解因式
    基于時域信號的三電平逆變器復合故障診斷
    測控技術(2018年11期)2018-12-07 05:49:02
    頻域稀疏毫米波人體安檢成像處理和快速成像稀疏陣列設計
    雷達學報(2018年3期)2018-07-18 02:41:34
    運用結(jié)構的齊次化,選換主元解題
    文理導航(2018年2期)2018-01-22 19:23:54
    基于極大似然準則與滾動時域估計的自適應UKF算法
    基于改進Radon-Wigner變換的目標和拖曳式誘餌頻域分離
    基于時域逆濾波的寬帶脈沖聲生成技術
    一種基于頻域的QPSK窄帶干擾抑制算法
    基于頻域伸縮的改進DFT算法
    電測與儀表(2015年3期)2015-04-09 11:37:24
    国产黄频视频在线观看| 亚洲精品美女久久av网站| 国产精品一区二区在线不卡| 久久性视频一级片| 女人被躁到高潮嗷嗷叫费观| 王馨瑶露胸无遮挡在线观看| 成人三级做爰电影| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲激情五月婷婷啪啪| 满18在线观看网站| 十八禁网站网址无遮挡| www.999成人在线观看| 人妻一区二区av| 久久久久精品人妻al黑| 国产一卡二卡三卡精品| 亚洲视频免费观看视频| 亚洲精品久久成人aⅴ小说| 亚洲成人免费av在线播放| 少妇精品久久久久久久| 国产野战对白在线观看| 午夜影院在线不卡| 亚洲成人免费av在线播放| 美女高潮到喷水免费观看| 亚洲五月婷婷丁香| 久久人人爽av亚洲精品天堂| 国产亚洲一区二区精品| 宅男免费午夜| 国产色视频综合| 免费av中文字幕在线| 国产成人免费无遮挡视频| 国产日韩欧美视频二区| 国产淫语在线视频| 精品国产一区二区三区久久久樱花| 男人操女人黄网站| xxxhd国产人妻xxx| 欧美 日韩 精品 国产| 热99久久久久精品小说推荐| 又大又爽又粗| 久久影院123| 丝袜脚勾引网站| 日本91视频免费播放| 欧美人与善性xxx| 1024视频免费在线观看| 制服人妻中文乱码| 建设人人有责人人尽责人人享有的| 亚洲国产看品久久| 亚洲,欧美精品.| 久久av网站| 免费看不卡的av| 一区二区日韩欧美中文字幕| 亚洲国产精品国产精品| 国产成人系列免费观看| 亚洲国产精品999| 性色av乱码一区二区三区2| 欧美日韩成人在线一区二区| 国产成人一区二区三区免费视频网站 | 十分钟在线观看高清视频www| videosex国产| 欧美激情极品国产一区二区三区| 欧美成人精品欧美一级黄| 午夜精品国产一区二区电影| 女人被躁到高潮嗷嗷叫费观| 国产视频一区二区在线看| 我要看黄色一级片免费的| 亚洲色图综合在线观看| 欧美日韩视频高清一区二区三区二| 日本猛色少妇xxxxx猛交久久| 香蕉国产在线看| 久久九九热精品免费| 大香蕉久久成人网| 男的添女的下面高潮视频| 久久国产精品大桥未久av| 母亲3免费完整高清在线观看| 成在线人永久免费视频| 母亲3免费完整高清在线观看| 精品国产一区二区三区久久久樱花| 欧美激情高清一区二区三区| 女性生殖器流出的白浆| 极品少妇高潮喷水抽搐| 91麻豆精品激情在线观看国产 | 亚洲一区中文字幕在线| 亚洲国产精品成人久久小说| 国产精品久久久久成人av| 性色av乱码一区二区三区2| 另类精品久久| bbb黄色大片| av网站免费在线观看视频| 看十八女毛片水多多多| 在线观看免费视频网站a站| 少妇被粗大的猛进出69影院| 亚洲欧美色中文字幕在线| 黄色怎么调成土黄色| 久久精品久久精品一区二区三区| 亚洲色图综合在线观看| 建设人人有责人人尽责人人享有的| 日韩电影二区| 波野结衣二区三区在线| 一区在线观看完整版| 中文精品一卡2卡3卡4更新| 欧美中文综合在线视频| 精品久久久久久电影网| 女性生殖器流出的白浆| 亚洲一码二码三码区别大吗| 一区在线观看完整版| 中文字幕色久视频| 国产日韩一区二区三区精品不卡| 中文乱码字字幕精品一区二区三区| 免费久久久久久久精品成人欧美视频| 咕卡用的链子| 欧美激情极品国产一区二区三区| 一区二区av电影网| 9色porny在线观看| 亚洲五月色婷婷综合| 日韩视频在线欧美| 中文字幕人妻丝袜一区二区| 国产一区二区激情短视频 | 成人午夜精彩视频在线观看| 国产高清不卡午夜福利| 91国产中文字幕| 亚洲免费av在线视频| 亚洲精品国产一区二区精华液| 又粗又硬又长又爽又黄的视频| 黑人猛操日本美女一级片| 啦啦啦在线免费观看视频4| 亚洲成人免费电影在线观看 | 午夜激情久久久久久久| 一区二区三区乱码不卡18| www.自偷自拍.com| 国产亚洲欧美在线一区二区| 只有这里有精品99| 中国国产av一级| 黄色视频在线播放观看不卡| 亚洲国产看品久久| 成年美女黄网站色视频大全免费| 亚洲欧洲国产日韩| 丁香六月欧美| 热99国产精品久久久久久7| 亚洲图色成人| 1024香蕉在线观看| 精品国产乱码久久久久久男人| 亚洲黑人精品在线| 观看av在线不卡| 各种免费的搞黄视频| 桃花免费在线播放| 久久人人爽av亚洲精品天堂| a 毛片基地| 成年动漫av网址| 大型av网站在线播放| 中文字幕色久视频| 婷婷丁香在线五月| 亚洲成国产人片在线观看| 久久热在线av| 欧美精品av麻豆av| a级毛片在线看网站| 成人手机av| 亚洲五月婷婷丁香| www.999成人在线观看| 少妇裸体淫交视频免费看高清 | 女性生殖器流出的白浆| 1024视频免费在线观看| 可以免费在线观看a视频的电影网站| 男男h啪啪无遮挡| 日韩中文字幕欧美一区二区 | 国产亚洲欧美精品永久| 国产1区2区3区精品| 国产成人欧美| 80岁老熟妇乱子伦牲交| 叶爱在线成人免费视频播放| 麻豆av在线久日| 国产成人免费观看mmmm| 亚洲欧洲精品一区二区精品久久久| 男男h啪啪无遮挡| 一边摸一边做爽爽视频免费| 精品人妻1区二区| 99国产精品免费福利视频| 欧美日韩亚洲高清精品| 狠狠精品人妻久久久久久综合| 黑丝袜美女国产一区| 精品欧美一区二区三区在线| 欧美日韩黄片免| 午夜福利在线免费观看网站| 国产男人的电影天堂91| 国产高清视频在线播放一区 | 欧美精品一区二区免费开放| 亚洲自偷自拍图片 自拍| 欧美黑人精品巨大| 少妇人妻 视频| 男人爽女人下面视频在线观看| kizo精华| 成人国产av品久久久| 日韩av免费高清视频| 99热全是精品| 黑丝袜美女国产一区| 欧美中文综合在线视频| 国产高清视频在线播放一区 | 国产精品亚洲av一区麻豆| 国产精品99久久99久久久不卡| 人妻一区二区av| 国产成人91sexporn| 久久久久网色| 国产成人精品在线电影| 中国国产av一级| 久久精品国产a三级三级三级| 天天操日日干夜夜撸| 久久热在线av| 日韩视频在线欧美| 久久久精品国产亚洲av高清涩受| 99热全是精品| 色网站视频免费| 国产成人啪精品午夜网站| 啦啦啦啦在线视频资源| 满18在线观看网站| 99久久综合免费| 精品视频人人做人人爽| 日本wwww免费看| 丝袜脚勾引网站| 十八禁网站网址无遮挡| √禁漫天堂资源中文www| 悠悠久久av| 狠狠婷婷综合久久久久久88av| 国产高清videossex| 日韩制服丝袜自拍偷拍| 一级,二级,三级黄色视频| 国精品久久久久久国模美| 亚洲国产毛片av蜜桃av| 最黄视频免费看| 一本综合久久免费| 成年人免费黄色播放视频| 国产成人欧美在线观看 | 日韩伦理黄色片| 成年人午夜在线观看视频| 亚洲久久久国产精品| 最新在线观看一区二区三区 | 日本a在线网址| 欧美久久黑人一区二区| 美女午夜性视频免费| 国产一区二区三区av在线| 久久久久久久久免费视频了| 亚洲国产欧美网| 色综合欧美亚洲国产小说| 好男人视频免费观看在线| 国产日韩欧美视频二区| 男女无遮挡免费网站观看| 51午夜福利影视在线观看| 日本vs欧美在线观看视频| av国产精品久久久久影院| 999精品在线视频| 熟女av电影| 91老司机精品| 日日夜夜操网爽| 少妇 在线观看| 丝瓜视频免费看黄片| 18禁国产床啪视频网站| 国产一区二区三区av在线| 亚洲精品av麻豆狂野| 乱人伦中国视频| 欧美在线一区亚洲| 亚洲欧洲精品一区二区精品久久久| 亚洲av综合色区一区| 黄色视频不卡| 日韩一卡2卡3卡4卡2021年| 在线 av 中文字幕| 赤兔流量卡办理| 国产熟女午夜一区二区三区| 丁香六月欧美| 丝袜喷水一区| 最近中文字幕2019免费版| av又黄又爽大尺度在线免费看| 一边摸一边做爽爽视频免费| www.999成人在线观看| 性少妇av在线| 亚洲精品日本国产第一区| 国产精品国产av在线观看| 久久久久久久久久久久大奶| 亚洲精品国产一区二区精华液| 久久久久网色| 国产淫语在线视频| 天天操日日干夜夜撸| 夫妻午夜视频| 国产视频一区二区在线看| 亚洲国产精品999| 亚洲色图 男人天堂 中文字幕| 色视频在线一区二区三区| 天堂俺去俺来也www色官网| 热99久久久久精品小说推荐| 久久精品成人免费网站| 成人18禁高潮啪啪吃奶动态图| 亚洲精品在线美女| 国产欧美日韩精品亚洲av| 操出白浆在线播放| 亚洲欧美一区二区三区久久| 国产日韩欧美亚洲二区| 大型av网站在线播放| 亚洲国产精品一区三区| 午夜福利乱码中文字幕| 91麻豆av在线| 国产精品麻豆人妻色哟哟久久| 如日韩欧美国产精品一区二区三区| 热99国产精品久久久久久7| 国产免费现黄频在线看| 天堂8中文在线网| 日韩精品免费视频一区二区三区| 十分钟在线观看高清视频www| 亚洲国产成人一精品久久久| 国产色视频综合| 黄色视频在线播放观看不卡| 最近最新中文字幕大全免费视频 | a级片在线免费高清观看视频| 国产精品亚洲av一区麻豆| 久久久久久人人人人人| 久久九九热精品免费| 黄色视频在线播放观看不卡| 最近最新中文字幕大全免费视频 | 满18在线观看网站| av不卡在线播放| 看十八女毛片水多多多| 一级毛片女人18水好多 | 亚洲,一卡二卡三卡| 国产91精品成人一区二区三区 | 亚洲国产av影院在线观看| 热re99久久精品国产66热6| 久久99一区二区三区| 王馨瑶露胸无遮挡在线观看| av有码第一页| 久久 成人 亚洲| 男女床上黄色一级片免费看| 中国国产av一级| av欧美777| 久久久久久亚洲精品国产蜜桃av| 日韩制服丝袜自拍偷拍| 精品熟女少妇八av免费久了| 亚洲五月色婷婷综合| 一二三四社区在线视频社区8| 日韩中文字幕视频在线看片| 三上悠亚av全集在线观看| 国产亚洲精品久久久久5区| 99久久精品国产亚洲精品| 精品少妇久久久久久888优播| 啦啦啦啦在线视频资源| av视频免费观看在线观看| 亚洲欧洲国产日韩| 男人舔女人的私密视频| 亚洲一卡2卡3卡4卡5卡精品中文| 悠悠久久av| 国产黄频视频在线观看| 啦啦啦啦在线视频资源| 如日韩欧美国产精品一区二区三区| 少妇 在线观看| 亚洲av欧美aⅴ国产| 2021少妇久久久久久久久久久| 老司机影院毛片| 久久精品久久精品一区二区三区| 亚洲国产日韩一区二区| 黄色怎么调成土黄色| 国产伦理片在线播放av一区| 国产一区二区 视频在线| 少妇被粗大的猛进出69影院| 久久久久精品国产欧美久久久 | 热99国产精品久久久久久7| 成人免费观看视频高清| 国产精品 欧美亚洲| 亚洲专区中文字幕在线| 在线看a的网站| 夫妻午夜视频| 国产在线一区二区三区精| 欧美在线一区亚洲| 成人影院久久| 另类精品久久| 国产精品一区二区在线观看99| 一级毛片 在线播放| 久久久久国产精品人妻一区二区| 日韩,欧美,国产一区二区三区| 国产麻豆69| 秋霞在线观看毛片| 午夜福利在线免费观看网站| 大片电影免费在线观看免费| 欧美成人午夜精品| 亚洲人成77777在线视频| 一级片'在线观看视频| 另类精品久久| 久久精品国产亚洲av高清一级| 国产精品三级大全| 一区二区三区激情视频| 大陆偷拍与自拍| 欧美国产精品一级二级三级| 久久人人97超碰香蕉20202| 亚洲欧美精品自产自拍| 国产一区二区 视频在线| 精品亚洲成a人片在线观看| 日韩一本色道免费dvd| 免费一级毛片在线播放高清视频 | videos熟女内射| 欧美日韩成人在线一区二区| 人妻人人澡人人爽人人| 尾随美女入室| 女人精品久久久久毛片| 夫妻性生交免费视频一级片| 国产熟女午夜一区二区三区| 最黄视频免费看| 婷婷色综合www| 欧美日本中文国产一区发布| 各种免费的搞黄视频| 悠悠久久av| 久久天躁狠狠躁夜夜2o2o | 精品亚洲成a人片在线观看| 大香蕉久久成人网| 女性被躁到高潮视频| 岛国毛片在线播放| svipshipincom国产片| 亚洲国产欧美日韩在线播放| 亚洲黑人精品在线| 亚洲色图综合在线观看| 极品少妇高潮喷水抽搐| 在现免费观看毛片| 91麻豆精品激情在线观看国产 | 婷婷色av中文字幕| 性色av乱码一区二区三区2| 日本猛色少妇xxxxx猛交久久| 日韩制服丝袜自拍偷拍| 这个男人来自地球电影免费观看| 亚洲国产最新在线播放| 国产精品国产av在线观看| 欧美黄色淫秽网站| 婷婷色av中文字幕| 99国产精品一区二区蜜桃av | 下体分泌物呈黄色| 亚洲成色77777| 欧美久久黑人一区二区| 日本欧美视频一区| xxx大片免费视频| 国产淫语在线视频| 成人18禁高潮啪啪吃奶动态图| 亚洲国产精品999| 高潮久久久久久久久久久不卡| 在线天堂中文资源库| 亚洲色图综合在线观看| 亚洲av电影在线进入| av国产久精品久网站免费入址| 欧美日韩亚洲综合一区二区三区_| 男男h啪啪无遮挡| 国产高清不卡午夜福利| 欧美国产精品va在线观看不卡| 黑丝袜美女国产一区| 一本大道久久a久久精品| 大话2 男鬼变身卡| 9热在线视频观看99| 男女之事视频高清在线观看 | 韩国精品一区二区三区| 一区二区三区精品91| 欧美 亚洲 国产 日韩一| 日本黄色日本黄色录像| 伦理电影免费视频| 久久精品熟女亚洲av麻豆精品| 一区在线观看完整版| 精品国产一区二区三区四区第35| 久久久久久久久久久久大奶| 中国国产av一级| 午夜福利乱码中文字幕| 国产精品人妻久久久影院| 亚洲少妇的诱惑av| 国产成人91sexporn| 天天躁夜夜躁狠狠躁躁| 美女主播在线视频| 中文字幕制服av| 国产精品偷伦视频观看了| 久久国产精品影院| 日韩av不卡免费在线播放| 欧美精品啪啪一区二区三区 | 午夜激情av网站| 欧美在线一区亚洲| 夜夜骑夜夜射夜夜干| 免费观看av网站的网址| 9191精品国产免费久久| 久久精品人人爽人人爽视色| 亚洲色图综合在线观看| 人妻一区二区av| 男女国产视频网站| 少妇精品久久久久久久| 国产成人欧美在线观看 | 少妇裸体淫交视频免费看高清 | 国产精品九九99| 丰满迷人的少妇在线观看| 精品视频人人做人人爽| 国产麻豆69| 欧美成人精品欧美一级黄| 日韩大码丰满熟妇| 亚洲精品日本国产第一区| 少妇裸体淫交视频免费看高清 | 国产又爽黄色视频| 91麻豆av在线| 久久亚洲国产成人精品v| 最近中文字幕2019免费版| 国产日韩一区二区三区精品不卡| 久久99精品国语久久久| 国产成人一区二区在线| 亚洲第一青青草原| 成人国产一区最新在线观看 | 夫妻性生交免费视频一级片| 久久久亚洲精品成人影院| 女警被强在线播放| 亚洲精品国产一区二区精华液| 欧美 亚洲 国产 日韩一| 美国免费a级毛片| 在线亚洲精品国产二区图片欧美| 午夜福利在线免费观看网站| 99热网站在线观看| 国产免费视频播放在线视频| 免费在线观看影片大全网站 | a级毛片在线看网站| 欧美日韩福利视频一区二区| 男男h啪啪无遮挡| 亚洲精品在线美女| h视频一区二区三区| 亚洲一码二码三码区别大吗| 久久这里只有精品19| 99国产精品99久久久久| 搡老岳熟女国产| 中文字幕人妻丝袜制服| 亚洲综合色网址| 久久人人爽人人片av| 国产精品国产av在线观看| 女警被强在线播放| 丁香六月欧美| 妹子高潮喷水视频| 夫妻性生交免费视频一级片| 日本wwww免费看| 午夜免费鲁丝| 欧美日韩黄片免| 人体艺术视频欧美日本| 婷婷丁香在线五月| 亚洲精品日本国产第一区| 亚洲,欧美,日韩| 精品熟女少妇八av免费久了| 亚洲成人手机| 亚洲中文字幕日韩| 狠狠婷婷综合久久久久久88av| 在线av久久热| 午夜日韩欧美国产| 男女无遮挡免费网站观看| 精品人妻熟女毛片av久久网站| 大陆偷拍与自拍| 精品免费久久久久久久清纯 | cao死你这个sao货| 免费人妻精品一区二区三区视频| 久久精品亚洲av国产电影网| 夫妻午夜视频| 午夜激情av网站| 1024香蕉在线观看| 日本wwww免费看| 精品久久久精品久久久| 黑人巨大精品欧美一区二区蜜桃| 久久精品国产亚洲av涩爱| 色综合欧美亚洲国产小说| 九色亚洲精品在线播放| 一本一本久久a久久精品综合妖精| 欧美成狂野欧美在线观看| 美女国产高潮福利片在线看| 午夜福利一区二区在线看| 亚洲美女黄色视频免费看| 亚洲精品在线美女| 交换朋友夫妻互换小说| 不卡av一区二区三区| 午夜福利视频在线观看免费| 成年美女黄网站色视频大全免费| 69精品国产乱码久久久| 久久青草综合色| 精品一区二区三卡| 亚洲激情五月婷婷啪啪| 久久99热这里只频精品6学生| 精品一品国产午夜福利视频| 老汉色av国产亚洲站长工具| 亚洲精品久久久久久婷婷小说| 黄片播放在线免费| 建设人人有责人人尽责人人享有的| 欧美日韩综合久久久久久| 久久久国产一区二区| 国产亚洲av高清不卡| 亚洲七黄色美女视频| 大话2 男鬼变身卡| 国产日韩一区二区三区精品不卡| 亚洲,一卡二卡三卡| 欧美成人午夜精品| 免费观看人在逋| 国产精品一区二区免费欧美 | 婷婷成人精品国产| 亚洲五月婷婷丁香| 一二三四在线观看免费中文在| 日韩av不卡免费在线播放| 日韩中文字幕欧美一区二区 | 亚洲欧美清纯卡通| 久久久久久人人人人人| 男女之事视频高清在线观看 | 亚洲国产精品成人久久小说| 成年美女黄网站色视频大全免费| 欧美日韩一级在线毛片| 国产深夜福利视频在线观看| 国产精品成人在线| 亚洲精品国产区一区二| 18在线观看网站| 热re99久久精品国产66热6| 18在线观看网站| 欧美成人精品欧美一级黄| 免费观看a级毛片全部| 人人妻,人人澡人人爽秒播 | 免费在线观看视频国产中文字幕亚洲 | 久久久久网色| 又大又爽又粗| 建设人人有责人人尽责人人享有的| 免费在线观看黄色视频的| 人人妻,人人澡人人爽秒播 | 亚洲精品一卡2卡三卡4卡5卡 | 免费观看人在逋| 亚洲成国产人片在线观看| 国产一级毛片在线| 女人精品久久久久毛片|