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

    基于TVAR-HMM 的滾動(dòng)軸承故障診斷

    2010-12-06 12:11:38王國(guó)鋒李玉波秦旭達(dá)李啟銘
    關(guān)鍵詞:特征向量頻譜故障診斷

    王國(guó)鋒,李玉波,秦旭達(dá),喻 秀,李啟銘

    (天津大學(xué)機(jī)械工程學(xué)院,天津 300072)

    機(jī)械故障診斷在現(xiàn)代工業(yè)生產(chǎn)中起著越來(lái)越重要的作用.對(duì)反映機(jī)械設(shè)備狀況的振動(dòng)信號(hào)進(jìn)行信號(hào)分析及特征提取,然后利用數(shù)學(xué)模型進(jìn)行故障識(shí)別是目前故障診斷的主要方法.

    滾動(dòng)軸承是機(jī)械設(shè)備中故障常發(fā)部件,多變的工況條件以及設(shè)備自身的非線性,使反映軸承故障的動(dòng)態(tài)信號(hào)表現(xiàn)出較強(qiáng)的非平穩(wěn)特性.目前信號(hào)主要的時(shí)頻分析方法有短時(shí)傅里葉變換、Wigner-Ville 分布、小波尺度分析等.這些方法的時(shí)頻分辨率都要受到不確定性原理的限制,因此無(wú)法同時(shí)在時(shí)域和頻域均得到很好的分辨率[1].

    故障識(shí)別模型以神經(jīng)網(wǎng)絡(luò)模型應(yīng)用最廣,但利用神經(jīng)網(wǎng)絡(luò)模型進(jìn)行故障識(shí)別需要大量的訓(xùn)練樣本,在有限的訓(xùn)練樣本情況下,泛化能力差,且優(yōu)化過(guò)程存在陷入局部極值的可能[2];另外,大部分神經(jīng)網(wǎng)絡(luò)模型處理的是靜態(tài)模式分類問(wèn)題,而旋轉(zhuǎn)機(jī)械的故障行為是一個(gè)動(dòng)態(tài)變化的過(guò)程.

    為此,筆者提出了應(yīng)用時(shí)變自回歸(time-varying autoregressive ,TVAR) 和 隱 馬 爾 科 夫 模 型 (hidden Markov model,HMM)的滾動(dòng)軸承故障診斷方法.時(shí)變自回歸模型是一種基于參數(shù)模型的時(shí)頻分析方法.它是在 AR(autoregressive)建模的基礎(chǔ)上發(fā)展起來(lái)的,TVAR 方法得出的時(shí)頻譜圖具有分辨率高、無(wú)交叉干擾項(xiàng)以及計(jì)算速度快等優(yōu)點(diǎn),非常適合處理非平穩(wěn)信號(hào),已在故障診斷領(lǐng)域有所應(yīng)用[3-5].隱馬爾科夫模型是信號(hào)動(dòng)態(tài)時(shí)間序列統(tǒng)計(jì)模型,能夠?qū)σ粋€(gè)時(shí)間跨度上的信息進(jìn)行統(tǒng)計(jì)建模和分類,具有嚴(yán)謹(jǐn)?shù)臄?shù)據(jù)結(jié)構(gòu)和可靠的計(jì)算性能.它能通過(guò)較少的樣本訓(xùn)練出可靠的模型,并按模式匹配原理,尋找與未知信號(hào)最相似的模式作為識(shí)別結(jié)果,現(xiàn)已成為語(yǔ)音識(shí)別的主流技術(shù)[6].傳統(tǒng)故障診斷方法一直停留在靜態(tài)分析基礎(chǔ)上,忽略了故障變化發(fā)展的動(dòng)態(tài)信息,而 HMM 用于故障診斷,可對(duì)故障變化的動(dòng)態(tài)信息進(jìn)行統(tǒng)計(jì)分析,能夠盡早地發(fā)現(xiàn)故障發(fā)展的跡象.與神經(jīng)網(wǎng)絡(luò)相比,保留了更多訓(xùn)練數(shù)據(jù)的統(tǒng)計(jì)信息,具有更高的識(shí)別率和魯棒性[7].

    筆者結(jié)合了兩種模型的優(yōu)點(diǎn),研究了非平穩(wěn)信號(hào)的時(shí)變參數(shù)建模以及隱馬爾科夫建模故障識(shí)別過(guò)程,并通過(guò)實(shí)驗(yàn)驗(yàn)證該方法在滾動(dòng)軸承故障診斷中的可靠性.

    1 時(shí)變自回歸模型

    {x(t),t=1,2,…,N }是一非平穩(wěn)時(shí)間序列,TVAR是利用先前值 x (t ? i)來(lái)估計(jì)現(xiàn)在值 x(t),可寫(xiě)成

    式中:x (t)為信號(hào)的采集值;ai(t)(i=1,2,…,p)為時(shí)變參數(shù);p 為模型階數(shù);e(t)為模型殘差,服從零均值正態(tài)分布.

    將時(shí)變參數(shù){ai(t),i=1,2,…,p}作為一組基時(shí)間函數(shù) gj(t )的線性組合,即

    式中:m 為擴(kuò)展維數(shù);aij為組合的時(shí)不變權(quán)值.

    將式(2)代入式(1),得

    式中

    這樣模型時(shí)變參數(shù)的估計(jì)問(wèn)題就變成常參數(shù) aij的估計(jì),從而把一個(gè)線性非平穩(wěn)問(wèn)題轉(zhuǎn)化為一個(gè)線性時(shí)不變問(wèn)題.

    文獻(xiàn)[8]闡述了利用最小二乘法估計(jì)模型參數(shù) A的過(guò)程,在估計(jì)出參數(shù) A 之后就可以根據(jù)式(2)計(jì)算出時(shí)變參數(shù) ai(t)(i =1,2,…,p),繼而求解模型方差,即

    瞬時(shí)頻率能夠通過(guò) TVAR 參數(shù) ai(t )提?。酶鱾€(gè)時(shí)刻的參數(shù)依次計(jì)算出每個(gè)時(shí)刻的參數(shù)譜圖,疊加在一起就可以構(gòu)造一個(gè)三維時(shí)頻譜圖,因此TVAR 譜函數(shù)定義為

    2 隱馬爾科夫模型

    HMM 是一個(gè)雙內(nèi)嵌式隨機(jī)過(guò)程,即 HMM 由兩個(gè)隨機(jī)過(guò)程組成:一是隱狀態(tài)轉(zhuǎn)移序列,它對(duì)應(yīng)一個(gè)單純的 Markov 過(guò)程;二是與隱狀態(tài)有關(guān)的觀測(cè)序列.其中隱狀態(tài)轉(zhuǎn)移序列是不可觀測(cè)的,只能通過(guò)另一個(gè)隨機(jī)過(guò)程的輸出觀測(cè)序列進(jìn)行推斷,所以稱之為隱馬爾科夫模型[9].

    HMM 可以記為 λ={Nλ,Mλ,πλ,Aλ,Bλ}或簡(jiǎn)寫(xiě)為λ={πλ,Aλ,Bλ}.其中 Nλ為模型中 Markov 鏈的狀態(tài)數(shù)目;Mλ為每個(gè)狀態(tài)對(duì)應(yīng)的可能的觀測(cè)值數(shù)目;πλ為初始概率分布矢量;Aλ為狀態(tài)轉(zhuǎn)移概率矩陣;Bλ為觀測(cè)值概率矩陣,對(duì)于連續(xù)的 HMM,Bλ是一組觀測(cè)值概率函數(shù).

    圖1為1 個(gè)三狀態(tài)的離散HMM 的例子,由圖可以看出,1 個(gè)HMM 可以由1 組狀態(tài)轉(zhuǎn)移概率以及與每一個(gè)狀態(tài)有關(guān)的觀測(cè)分布概率來(lái)定義.

    圖1 三狀態(tài)HMM示意Fig.1 Diagram of three-state HMM

    一個(gè)標(biāo)準(zhǔn)的HMM 模型需要解決概率計(jì)算、模型訓(xùn)練和隱狀態(tài)估計(jì)3 個(gè)基本問(wèn)題.

    (1)概率計(jì)算.給定模型參數(shù) λ={πλ,Aλ,Bλ}、觀測(cè)序列 O={o1,o2,…,ot},計(jì)算觀測(cè)序列 O 在給定模型下的概率,即計(jì)算似然 P(O︱λ)或者對(duì)數(shù)似然ln,P(O︱λ),它表示參數(shù) λ 擬合數(shù)據(jù) O 的準(zhǔn)確程度.可通過(guò)前向-后向算法對(duì)HMM 進(jìn)行遍歷實(shí)現(xiàn).

    (2)模型訓(xùn)練,即模型參數(shù)估計(jì)問(wèn)題.也就是說(shuō),對(duì)于初始模型和給定用于訓(xùn)練的觀測(cè)序列 O={o1,o2,…,ot},如何調(diào)整模型 λ={πλ,Aλ,Bλ}的參數(shù),使其能夠最好地?cái)M合觀測(cè)數(shù)據(jù) ? argmaxλ=f (O?)λ,這可以通過(guò)EM 算法或Baum-Welch 算法實(shí)現(xiàn).

    (3)隱狀態(tài)估計(jì).給定模型參數(shù) λ={πλ,Aλ,Bλ}和觀測(cè)序列 O={o1,o2,…,ot},基于某種最優(yōu)準(zhǔn)則估計(jì)最可能的隱狀態(tài)序列 Q={q1,q2,…,qt},即估計(jì)產(chǎn)生觀測(cè)序列的最可能經(jīng)過(guò)路徑,這可以通過(guò) Viterbi 算法估計(jì)實(shí)現(xiàn).

    3 TVAR-HMM故障診斷步驟

    圖2所示為利用 TVAR-HMM 進(jìn)行故障診斷的主要步驟.

    圖2 TVAR-HMM故障診斷流程Fig.2 Flow chart of fault diagnosis based on TVAR-HMM

    (1)對(duì)不同故障類型的振動(dòng)信號(hào)進(jìn)行 TVAR 建模,得到時(shí)頻譜圖.建模過(guò)程中,需要選定基函數(shù)gj(t)、模型階數(shù) p 和擴(kuò)展維數(shù)m ,文獻(xiàn)[10]詳細(xì)介紹了它們的選擇方法.

    (2)基于能量法提取故障特征.對(duì)信號(hào) TVAR 建模后,故障信息完全體現(xiàn)在時(shí)頻譜圖上.用能量法提取故障特征可以更加有效地提取時(shí)頻圖上的故障信息.文獻(xiàn)[11]已研究了這種方法的有效性.其主要過(guò)程是:將時(shí)頻空間劃分為若干區(qū)域,分別計(jì)算其時(shí)頻局部能量,將時(shí)頻局部能量作為元素,組成特征向量.

    (3)用特征向量訓(xùn)練 HMM 模型.特征向量為HMM 模型的觀測(cè)序列,故障類型為隱狀態(tài),1 種故障類型訓(xùn)練1 個(gè)HMM 模型.

    (4)利用訓(xùn)練好的模型對(duì)測(cè)試樣本進(jìn)行識(shí)別,即對(duì)于未知的觀測(cè)序列,分別計(jì)算在各個(gè)模型下的輸出概率,然后進(jìn)行比較,選取輸出概率最大的模型,即代表識(shí)別出的故障類型.

    4 實(shí)驗(yàn)研究

    為了驗(yàn)證本文方法在實(shí)際故障診斷中的應(yīng)用情況,筆者進(jìn)行了實(shí)驗(yàn)分析.實(shí)驗(yàn)在一個(gè)故障模擬平臺(tái)上進(jìn)行,包括驅(qū)動(dòng)電機(jī)、軸承座、扭矩加載器和采集設(shè)備,見(jiàn)圖 3.通過(guò)更換正常軸承和故障軸承進(jìn)行故障診斷實(shí)驗(yàn).實(shí)驗(yàn)測(cè)試的軸承是深溝球軸承,型號(hào)為6205-2RS,品牌為SKF,其內(nèi)圈直徑為 25,mm,外圈直徑為52,mm,滾子直徑為7.94,mm.軸承故障用電火花加工而成,故障直徑分別為 0.177,8,mm、0.355,6,mm、0.533,4,mm,3 個(gè)等級(jí),深度都為 0.279,4,mm.電機(jī)轉(zhuǎn)速為 1,797,r/min.振動(dòng)信號(hào)由加速度傳感器拾取,分別采集安裝無(wú)故障軸承、外圈故障軸承、內(nèi)圈故障軸承、滾子故障軸承時(shí)的振動(dòng)信號(hào),各測(cè)得數(shù)據(jù) 50 組,每組2,048 點(diǎn),采樣頻率為12,000,Hz.

    圖3 實(shí)驗(yàn)裝置Fig.3 Experimental device

    4.1 TVAR建模與特征提取

    對(duì)無(wú)故障、內(nèi)圈故障、外圈故障、滾子故障4 種狀態(tài)信號(hào) TVAR 建模時(shí),參照文獻(xiàn)[10]中的論述,基函數(shù)選擇傅里葉函數(shù),采用AIC準(zhǔn)則確定模型階數(shù)p.在確定模型階數(shù)前必須先確定擴(kuò)展維數(shù) m,為了保證信號(hào)時(shí)頻譜的高分辨率和計(jì)算速度,一般模型的擴(kuò)展維數(shù)取為 8~12.本次實(shí)驗(yàn)選取擴(kuò)展維數(shù) m=10,依據(jù)AIC 準(zhǔn)則確定軸承無(wú)故障、內(nèi)圈故障、外圈故障、滾子故障狀態(tài)下模型階數(shù) p 分別為 11、15、15、12.

    4 種類型的振動(dòng)信號(hào)TVAR 建模后,可以得到相應(yīng)的時(shí)頻譜,如圖 4(a)~(d)所示.由時(shí)頻譜中可以看出,不同故障類型軸承的振動(dòng)信號(hào)的頻率成分是不同的,并且軸承存在故障時(shí),時(shí)頻譜中存在時(shí)變性與沖擊性.圖 4(e)為軸承外圈故障信號(hào)時(shí)域波形,圖4(f)為軸承外圈故障信號(hào)經(jīng)傅里葉變換后得到的頻譜.通過(guò)圖4(d)~(f)的比較可以明顯看出圖4(d)與4(f)中信號(hào)的頻率成分保持一致,圖 4(d)與 4(e)中頻率成分發(fā)生沖擊的時(shí)間及時(shí)間間隔基本一致.這說(shuō)明經(jīng) TVAR 建模得到的時(shí)頻譜具有較高的時(shí)頻分辨率.

    基于能量法對(duì)時(shí)頻圖進(jìn)行特征提取時(shí),考慮到速度與精度雙重因素,經(jīng)過(guò)多次實(shí)驗(yàn)分析,選擇將時(shí)頻譜分為16 部分(過(guò)多會(huì)影響HMM 訓(xùn)練速度;過(guò)少會(huì)影響測(cè)試的準(zhǔn)確度),求每一部分能量的均值,將均值序列歸一化作為特征向量.通過(guò)特征提取,每種狀態(tài)得到 50 組 16 維的特征向量,4 種狀態(tài)共得到 4×50組特征向量.

    圖4 4種狀態(tài)的時(shí)頻譜Fig.4 Time-frequency spectra of four states

    4.2 HMM建模與故障識(shí)別

    建模時(shí),采用 4 個(gè)隱狀態(tài)代表無(wú)故障軸承、內(nèi)圈故障軸承、外圈故障軸承和滾子故障軸承,分別記作λ1、λ2、λ3、λ4.模型的觀測(cè)序列是 TVAR 建模后得到的 16 維特征向量.初始概率分布矢量、初始狀態(tài)轉(zhuǎn)移概率矩陣和初始觀測(cè)值概率矩陣均由 Matlab 中rand 函數(shù)隨機(jī)取得,然后歸一化.

    選用120 組特征向量(每種狀態(tài)30 組)作為訓(xùn)練樣本,訓(xùn)練前需使用Lloyd 算法對(duì)特征向量進(jìn)行標(biāo)量量化.訓(xùn)練算法為Baum-Welch 算法. HMM 的訓(xùn)練過(guò)程中,隨著迭代次數(shù)的增加,最大對(duì)數(shù)似然估計(jì)值也在不斷增加,直到達(dá)到收斂誤差為止.訓(xùn)練結(jié)束后,得到4 個(gè)隱狀態(tài)對(duì)應(yīng)的4 個(gè)HMM 識(shí)別模型.圖5 給出了4 種狀態(tài)的HMM 訓(xùn)練曲線,所有狀態(tài)在迭代次數(shù)20 都達(dá)到了收斂,收斂速度快.

    圖5 HMM訓(xùn)練曲線Fig.5 HMM training curves

    在訓(xùn)練好各狀態(tài) HMM 以后,將剩余的 80 組特征向量(每種狀態(tài)20 組)作為測(cè)試樣本輸入到各狀態(tài)HMM 模型中(測(cè)試前需使用 Lloyd 算法對(duì)特征向量進(jìn)行標(biāo)量量化編碼),在各個(gè)模型下能夠各自輸出 1個(gè)對(duì)數(shù)似然概率值.對(duì)數(shù)似然概率值反映了特征向量與各 HMM 的相似程度,對(duì)數(shù)似然概率值越大,特征向量越接近該狀態(tài) HMM,特征向量屬于使輸出對(duì)數(shù)似然概率值最大的模型所對(duì)應(yīng)的故障類型.識(shí)別算法選用 Viterbi 算法.表 1 給出了部分測(cè)試樣本的測(cè)試結(jié)果,通過(guò)得到的對(duì)數(shù)似然概率值,可以準(zhǔn)確識(shí)別出軸承的故障類型;對(duì)于同一種故障類型,故障程度越高,得到的對(duì)數(shù)似然概率值越大.

    為了比較基于HMM 和基于BP 網(wǎng)絡(luò)、RBF 網(wǎng)絡(luò)故障診斷模型的識(shí)別效果,將本實(shí)驗(yàn)的數(shù)據(jù)分別用BP 網(wǎng)絡(luò)、RBF 網(wǎng)絡(luò)故障診斷模型來(lái)識(shí)別.網(wǎng)絡(luò)的層數(shù)設(shè)置為3 層,輸入層神經(jīng)元個(gè)數(shù)為16 個(gè),輸出層神經(jīng)元個(gè)數(shù)為 4 個(gè).采用與基于 HMM 故障診斷模型同樣的訓(xùn)練數(shù)據(jù)和測(cè)試數(shù)據(jù)分別訓(xùn)練和測(cè)試 BP 網(wǎng)絡(luò)和RBF 網(wǎng)絡(luò),得出診斷結(jié)果.

    圖6為本次實(shí)驗(yàn)中 3 種模型的故障識(shí)別正確率比較.由圖6 可見(jiàn),基于 HMM 故障診斷模型的有效識(shí)別率明顯好于基于 BP 網(wǎng)絡(luò)或 RBF 網(wǎng)絡(luò)故障診斷模型.實(shí)驗(yàn)結(jié)果表明,基于 HMM 故障診斷模型可以成功地診斷旋轉(zhuǎn)機(jī)械的軸承故障,識(shí)別率較高,穩(wěn)定性好.

    表1 部分測(cè)試結(jié)果Tab.1 Some testing results

    圖6 HMM、BP、RBF識(shí)別有效性比較Fig.6 Diagnosis accuracies of HMM、BP、RBF

    5 結(jié) 論

    (1)TVAR 模型非常適合非平穩(wěn)信號(hào)的處理,時(shí)變參數(shù)體現(xiàn)了信號(hào)頻率的時(shí)變性和沖擊性.TVAR 建模分析得到的信號(hào)時(shí)頻譜的分辨率高,且無(wú)交叉項(xiàng)的干擾,計(jì)算速度快.

    (2)HMM 模型可以成功地對(duì)旋轉(zhuǎn)機(jī)械的軸承故障進(jìn)行診斷,識(shí)別率高.由于其為動(dòng)態(tài)時(shí)間序列統(tǒng)計(jì)模型,非常適合識(shí)別信息量大、非平穩(wěn)性強(qiáng)、特征重復(fù)再現(xiàn)性不佳的故障信號(hào),在旋轉(zhuǎn)機(jī)械故障診斷中有很好的應(yīng)用前景.

    [1]李 強(qiáng). 非平穩(wěn)信號(hào)特征提取新方法和實(shí)用診斷技術(shù)的研究[D]. 天津:天津大學(xué)機(jī)械工程學(xué)院,2003.Li Qiang. Research on New Method of Non-Stationary Signal Feature Extraction and Practical Diagnostic Techniques[D]. Tianjin:School of Mechanical Engineering,Tianjin University,2003(in Chinese).

    [2]王奉濤. 非平穩(wěn)信號(hào)故障特征提取與智能診斷方法的研究與應(yīng)用[D]. 大連:大連理工大學(xué)機(jī)械工程學(xué)院,2003.Wang Fengtao. Research on Non-Stationary Signal Feature Extraction and Intelligent Diagnostic Method[D].Dalian:School of Mechanical Engineering,Dalian University of Technology,2003(in Chinese).

    [3]陳 慧,張 龍,李 嶷. 基于 TVAR 的非平穩(wěn)工況轉(zhuǎn)子故障診斷技術(shù)研究[J]. 華東交通大學(xué)學(xué)報(bào),2006,23(5):116-118.Chen Hui,Zhang Long,Li Yi. TVAR based fault diagnosis for rotation machine under non-stationary condition[J]. Journal of East China Jiaotong University,2006,23(5):116-118(in Chinese).

    [4]Wang Guofeng,Luo Zhigao. Fault identification and classification of rolling element bearing based on timevarying autoregressive spectrum [J]. Mechanical Systems and Signal Progressing,2008,22(4):934-947.

    [5]張 龍,熊?chē)?guó)良,柳和生,等. 基于時(shí)變自回歸模型與支持向量機(jī)的旋轉(zhuǎn)機(jī)械故障診斷方法[J]. 中國(guó)電機(jī)工程學(xué)報(bào),2007,27(9):99-103.Zhang Long,Xiong Guoliang,Liu Hesheng,et al. TVAR and SVM based fault diagnosis method for rotation machine [J]. Proceedings of the CSEE,2007,27(9):99-103(in Chinese).

    [6]Abramovich Y I,Spencer N K,Turley M D E. Timevarying autoregressive(TVAR)models for multiple radar observations [J]. IEEE Transactions on Signal Processing,2007,55(4):1298-1311.

    [7]Tai A H,Ching Wai-Ki,Chan L Y. Detection of machine failure:Hidden Markov model approach[J]. Computers& Industrial Engineering,2009,57(2):608-619.

    [8]Hall M G,Oppenheim A V,Willsky A S. Time varying parametric modeling of speech [J]. Signal Processing,1983,5(3):267-285.

    [9]Ocak H,Loparo K A. A new bearing fault detection and diagnosis scheme based on hidden Markov modeling of vibration signals [C]//IEEE International Conference on Acoustics,Speech,and Signal Processing. Salt Lake City,USA,2001,5:3141-3144.

    [10]Zhan Y M,Jardine A K S. Adaptive autoregressive modeling of non-stationary vibration signals under distinct gear states(Part 1):Modeling [J]. Journal of Sound and Vibration,2005,286:429-450.

    [11]劉 偉. 基于能量的時(shí)頻特征提取研究[J]. 紅外與激光工程,2004,33(3):296-299.Liu Wei. Research on the time-frequency feature extraction based on energy[J]. Infrared and Laser Project,2004,33(3):296-299(in Chinese).

    猜你喜歡
    特征向量頻譜故障診斷
    二年制職教本科線性代數(shù)課程的幾何化教學(xué)設(shè)計(jì)——以特征值和特征向量為例
    克羅內(nèi)克積的特征向量
    一種用于深空探測(cè)的Chirp變換頻譜分析儀設(shè)計(jì)與實(shí)現(xiàn)
    一種基于稀疏度估計(jì)的自適應(yīng)壓縮頻譜感知算法
    一類特殊矩陣特征向量的求法
    EXCEL表格計(jì)算判斷矩陣近似特征向量在AHP法檢驗(yàn)上的應(yīng)用
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    認(rèn)知無(wú)線電頻譜感知技術(shù)綜述
    基于LCD和排列熵的滾動(dòng)軸承故障診斷
    基于WPD-HHT的滾動(dòng)軸承故障診斷
    国产精品久久久av美女十八| av视频免费观看在线观看| 免费观看无遮挡的男女| 日日摸夜夜添夜夜爱| 日本爱情动作片www.在线观看| 免费不卡的大黄色大毛片视频在线观看| 亚洲色图 男人天堂 中文字幕| 电影成人av| 国产熟女午夜一区二区三区| 亚洲伊人久久精品综合| 人体艺术视频欧美日本| 少妇人妻久久综合中文| 婷婷色麻豆天堂久久| 欧美最新免费一区二区三区| 老司机影院毛片| 黄频高清免费视频| 欧美日本中文国产一区发布| 高清视频免费观看一区二区| av视频免费观看在线观看| 亚洲欧美色中文字幕在线| 蜜桃在线观看..| 人人妻人人爽人人添夜夜欢视频| 飞空精品影院首页| 中文字幕精品免费在线观看视频| 亚洲综合精品二区| 亚洲欧美日韩另类电影网站| 老司机影院成人| 男人爽女人下面视频在线观看| 咕卡用的链子| 国产精品亚洲av一区麻豆 | 天天躁夜夜躁狠狠躁躁| 美女中出高潮动态图| 老司机影院毛片| 久久精品久久久久久噜噜老黄| 国产精品一区二区在线不卡| 国产精品免费视频内射| 国产av一区二区精品久久| 成年美女黄网站色视频大全免费| 啦啦啦在线观看免费高清www| 美国免费a级毛片| 黑人巨大精品欧美一区二区蜜桃| 亚洲经典国产精华液单| 亚洲美女视频黄频| 国产视频首页在线观看| 精品国产国语对白av| 精品99又大又爽又粗少妇毛片| 久久精品国产亚洲av涩爱| 欧美亚洲 丝袜 人妻 在线| 久久久精品94久久精品| 侵犯人妻中文字幕一二三四区| 一区福利在线观看| 美女大奶头黄色视频| 啦啦啦在线观看免费高清www| 在现免费观看毛片| 亚洲国产看品久久| 美女大奶头黄色视频| 国产亚洲午夜精品一区二区久久| 国产日韩欧美视频二区| 久久毛片免费看一区二区三区| 男的添女的下面高潮视频| 午夜激情久久久久久久| 欧美日韩精品网址| 国语对白做爰xxxⅹ性视频网站| 亚洲欧洲精品一区二区精品久久久 | 精品一区在线观看国产| 久久久久精品久久久久真实原创| 国产成人精品久久久久久| 国产精品一区二区在线观看99| 免费看不卡的av| 亚洲精品一二三| 久久人人97超碰香蕉20202| 黄色 视频免费看| 精品99又大又爽又粗少妇毛片| 久久青草综合色| 久久热在线av| 丰满迷人的少妇在线观看| 欧美日韩视频高清一区二区三区二| 久久精品久久精品一区二区三区| 中文字幕av电影在线播放| 久久精品国产综合久久久| av不卡在线播放| 亚洲欧美日韩另类电影网站| 亚洲精品在线美女| 九九爱精品视频在线观看| 老熟女久久久| h视频一区二区三区| 你懂的网址亚洲精品在线观看| 五月伊人婷婷丁香| 99久久综合免费| 久久精品国产a三级三级三级| 人人妻人人澡人人爽人人夜夜| 国产精品秋霞免费鲁丝片| 亚洲成人手机| 国产色婷婷99| av在线观看视频网站免费| 伦理电影免费视频| 国产精品久久久久久精品古装| 在线精品无人区一区二区三| 丝袜人妻中文字幕| 国产男人的电影天堂91| 在线 av 中文字幕| 另类亚洲欧美激情| 亚洲精品久久久久久婷婷小说| 美女主播在线视频| 看免费av毛片| 曰老女人黄片| 少妇的逼水好多| 最新的欧美精品一区二区| 青青草视频在线视频观看| 蜜桃国产av成人99| 色94色欧美一区二区| 桃花免费在线播放| 下体分泌物呈黄色| 国产成人一区二区在线| 成人黄色视频免费在线看| av在线播放精品| 免费黄网站久久成人精品| 亚洲人成77777在线视频| 久久久久久伊人网av| 亚洲综合精品二区| 亚洲人成网站在线观看播放| 少妇被粗大猛烈的视频| 国产在视频线精品| 最黄视频免费看| 2018国产大陆天天弄谢| 午夜激情av网站| 午夜免费观看性视频| 久久久久久久久久久久大奶| xxxhd国产人妻xxx| 丁香六月天网| freevideosex欧美| 国产成人免费无遮挡视频| 国产深夜福利视频在线观看| 国产精品一区二区在线观看99| 日韩三级伦理在线观看| 国产亚洲一区二区精品| 久久久久久久久久久免费av| 欧美少妇被猛烈插入视频| 人妻系列 视频| 亚洲在久久综合| 精品国产乱码久久久久久小说| 国语对白做爰xxxⅹ性视频网站| 美女中出高潮动态图| 欧美亚洲 丝袜 人妻 在线| 久久久久久久精品精品| 国产日韩欧美亚洲二区| a级毛片在线看网站| 午夜免费鲁丝| 欧美av亚洲av综合av国产av | 亚洲国产av新网站| 亚洲综合精品二区| 亚洲第一青青草原| 精品久久久久久电影网| 侵犯人妻中文字幕一二三四区| 91精品三级在线观看| 一区二区av电影网| 中文乱码字字幕精品一区二区三区| 欧美变态另类bdsm刘玥| 搡女人真爽免费视频火全软件| 春色校园在线视频观看| 亚洲精品一二三| 久久精品久久久久久噜噜老黄| 久久午夜福利片| 日韩 亚洲 欧美在线| 亚洲欧美一区二区三区久久| 两个人看的免费小视频| 日日爽夜夜爽网站| 国产精品国产三级国产专区5o| 久久精品国产鲁丝片午夜精品| a级片在线免费高清观看视频| 考比视频在线观看| 国产乱人偷精品视频| 极品少妇高潮喷水抽搐| a级片在线免费高清观看视频| 国产av精品麻豆| 国产精品偷伦视频观看了| 麻豆av在线久日| 午夜福利网站1000一区二区三区| 免费人妻精品一区二区三区视频| 久久精品国产自在天天线| 熟女少妇亚洲综合色aaa.| 美女福利国产在线| 一级黄片播放器| 巨乳人妻的诱惑在线观看| 在线天堂中文资源库| 婷婷色麻豆天堂久久| 2018国产大陆天天弄谢| 中文字幕精品免费在线观看视频| 如何舔出高潮| 国产伦理片在线播放av一区| 视频在线观看一区二区三区| 精品一区在线观看国产| 亚洲精品第二区| 免费看不卡的av| 夫妻性生交免费视频一级片| 亚洲精品久久久久久婷婷小说| 日本av手机在线免费观看| 女人高潮潮喷娇喘18禁视频| 成人黄色视频免费在线看| 亚洲欧美清纯卡通| 中文字幕人妻丝袜一区二区 | 国产片内射在线| 中文字幕亚洲精品专区| 日韩av在线免费看完整版不卡| 日本av手机在线免费观看| 欧美av亚洲av综合av国产av | 久久综合国产亚洲精品| 欧美av亚洲av综合av国产av | 青春草亚洲视频在线观看| 如何舔出高潮| 丝袜人妻中文字幕| 免费在线观看视频国产中文字幕亚洲 | 国产探花极品一区二区| 欧美精品国产亚洲| 韩国精品一区二区三区| 汤姆久久久久久久影院中文字幕| 久久国内精品自在自线图片| 国产在线免费精品| 国产片内射在线| 伊人久久国产一区二区| 久久99蜜桃精品久久| 欧美 亚洲 国产 日韩一| 99re6热这里在线精品视频| 日韩大片免费观看网站| 香蕉精品网在线| 午夜福利网站1000一区二区三区| 9191精品国产免费久久| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品99久久99久久久不卡 | 菩萨蛮人人尽说江南好唐韦庄| 午夜激情av网站| 国产免费现黄频在线看| 国产日韩一区二区三区精品不卡| 国产亚洲精品第一综合不卡| 国产国语露脸激情在线看| 亚洲欧美中文字幕日韩二区| 国产av精品麻豆| 青青草视频在线视频观看| 免费少妇av软件| 波多野结衣一区麻豆| 亚洲国产色片| 精品国产露脸久久av麻豆| 亚洲国产精品一区三区| 午夜福利乱码中文字幕| 女人精品久久久久毛片| 欧美 日韩 精品 国产| 亚洲一码二码三码区别大吗| av国产久精品久网站免费入址| 免费看不卡的av| 亚洲精品乱久久久久久| 久久精品国产a三级三级三级| 免费高清在线观看视频在线观看| 97人妻天天添夜夜摸| 伦理电影大哥的女人| 97在线视频观看| 精品亚洲乱码少妇综合久久| 国产熟女午夜一区二区三区| 在线观看www视频免费| 亚洲欧美精品自产自拍| 亚洲精品视频女| 国产成人免费无遮挡视频| 午夜免费男女啪啪视频观看| 有码 亚洲区| 97在线人人人人妻| 最黄视频免费看| 秋霞在线观看毛片| 成人毛片60女人毛片免费| 国产1区2区3区精品| 国产av一区二区精品久久| 一级毛片 在线播放| 亚洲av综合色区一区| 男女下面插进去视频免费观看| 精品卡一卡二卡四卡免费| 蜜桃在线观看..| 大码成人一级视频| 日日啪夜夜爽| 精品一区二区三区四区五区乱码 | 久久久久视频综合| 丝瓜视频免费看黄片| 黄片无遮挡物在线观看| 最近的中文字幕免费完整| 精品国产一区二区三区久久久樱花| 午夜免费鲁丝| 欧美+日韩+精品| 亚洲精品国产一区二区精华液| 十八禁高潮呻吟视频| 哪个播放器可以免费观看大片| 欧美日韩精品网址| 午夜福利在线观看免费完整高清在| 一区二区三区四区激情视频| 深夜精品福利| 中文字幕最新亚洲高清| 赤兔流量卡办理| 777米奇影视久久| 日韩制服丝袜自拍偷拍| 国产一区二区三区av在线| 欧美日韩精品网址| 交换朋友夫妻互换小说| 亚洲av免费高清在线观看| 女人久久www免费人成看片| 国产精品一区二区在线不卡| 国产不卡av网站在线观看| 中文精品一卡2卡3卡4更新| 熟妇人妻不卡中文字幕| 观看av在线不卡| 亚洲三区欧美一区| 欧美日韩成人在线一区二区| 久久久久国产精品人妻一区二区| 日韩一本色道免费dvd| 免费日韩欧美在线观看| 中文字幕另类日韩欧美亚洲嫩草| 9191精品国产免费久久| 黑丝袜美女国产一区| 国产1区2区3区精品| 欧美+日韩+精品| 多毛熟女@视频| 国产免费又黄又爽又色| 少妇的丰满在线观看| 免费黄网站久久成人精品| 男人操女人黄网站| 妹子高潮喷水视频| 在线观看免费日韩欧美大片| 国产日韩欧美视频二区| 91午夜精品亚洲一区二区三区| www日本在线高清视频| 亚洲人成77777在线视频| 亚洲精品av麻豆狂野| 国产精品二区激情视频| 青青草视频在线视频观看| 免费看不卡的av| 国产精品 欧美亚洲| 在线观看一区二区三区激情| 七月丁香在线播放| 女人高潮潮喷娇喘18禁视频| 日韩大片免费观看网站| 国产精品久久久av美女十八| xxx大片免费视频| 老女人水多毛片| 看非洲黑人一级黄片| 性色av一级| 亚洲精品久久久久久婷婷小说| 九色亚洲精品在线播放| 伊人亚洲综合成人网| av女优亚洲男人天堂| 黄频高清免费视频| 久久久久网色| 午夜福利网站1000一区二区三区| 欧美日韩精品网址| 男女啪啪激烈高潮av片| 欧美日韩精品网址| 国产黄色免费在线视频| 韩国av在线不卡| 美女大奶头黄色视频| 亚洲欧美中文字幕日韩二区| 麻豆av在线久日| 亚洲欧美一区二区三区国产| 久久久久久久久久久久大奶| 侵犯人妻中文字幕一二三四区| 久久久亚洲精品成人影院| 中文乱码字字幕精品一区二区三区| 日本-黄色视频高清免费观看| 最近手机中文字幕大全| 久久久国产精品麻豆| 精品一区在线观看国产| 蜜桃国产av成人99| 满18在线观看网站| 国精品久久久久久国模美| 国产日韩欧美亚洲二区| 午夜福利网站1000一区二区三区| 久久ye,这里只有精品| 18禁动态无遮挡网站| 黄频高清免费视频| 国产日韩一区二区三区精品不卡| 久久人人97超碰香蕉20202| 麻豆av在线久日| 亚洲国产欧美网| 日韩精品免费视频一区二区三区| 久久久亚洲精品成人影院| 在线观看国产h片| 国产深夜福利视频在线观看| 午夜福利在线观看免费完整高清在| 国产一区二区 视频在线| 久久精品国产亚洲av天美| 午夜福利网站1000一区二区三区| 亚洲国产精品一区三区| 十八禁高潮呻吟视频| 国产黄色免费在线视频| 高清不卡的av网站| 精品国产一区二区三区四区第35| 91aial.com中文字幕在线观看| 国产在视频线精品| 深夜精品福利| 日韩三级伦理在线观看| 成人二区视频| 视频在线观看一区二区三区| 久久久久久免费高清国产稀缺| 制服诱惑二区| 九九爱精品视频在线观看| 亚洲国产av影院在线观看| 黑人巨大精品欧美一区二区蜜桃| 色婷婷av一区二区三区视频| 日本91视频免费播放| 十分钟在线观看高清视频www| 国产一区有黄有色的免费视频| 777米奇影视久久| 一区福利在线观看| 亚洲国产精品一区三区| 精品少妇黑人巨大在线播放| 精品少妇内射三级| 国产熟女欧美一区二区| 免费在线观看视频国产中文字幕亚洲 | 青春草亚洲视频在线观看| 深夜精品福利| 91精品伊人久久大香线蕉| 久久久久久久久久人人人人人人| 又黄又粗又硬又大视频| 天天躁夜夜躁狠狠久久av| 99热全是精品| 少妇被粗大猛烈的视频| 国产野战对白在线观看| 啦啦啦在线观看免费高清www| 午夜免费观看性视频| 国产精品一区二区在线不卡| 久久久久国产网址| 精品卡一卡二卡四卡免费| 中文字幕另类日韩欧美亚洲嫩草| 黄色视频在线播放观看不卡| 最新的欧美精品一区二区| 蜜桃国产av成人99| 少妇人妻久久综合中文| 亚洲国产看品久久| 激情五月婷婷亚洲| 欧美+日韩+精品| 国产1区2区3区精品| 欧美日韩综合久久久久久| 中文字幕精品免费在线观看视频| 久久免费观看电影| 亚洲国产欧美日韩在线播放| av片东京热男人的天堂| 日本欧美视频一区| 免费久久久久久久精品成人欧美视频| 人人妻人人添人人爽欧美一区卜| 亚洲视频免费观看视频| 最近中文字幕高清免费大全6| 国产在线免费精品| 我要看黄色一级片免费的| 久久久久久久久久久久大奶| 超碰97精品在线观看| 99香蕉大伊视频| 97在线人人人人妻| 欧美中文综合在线视频| 国产成人一区二区在线| 精品卡一卡二卡四卡免费| 永久网站在线| 一边摸一边做爽爽视频免费| 波多野结衣一区麻豆| 亚洲美女搞黄在线观看| 色吧在线观看| 国产精品香港三级国产av潘金莲 | 成人国语在线视频| 少妇的逼水好多| 最近中文字幕高清免费大全6| 免费高清在线观看日韩| 人妻少妇偷人精品九色| 成人黄色视频免费在线看| 亚洲精品一区蜜桃| 午夜激情久久久久久久| 欧美+日韩+精品| 好男人视频免费观看在线| 寂寞人妻少妇视频99o| 国产精品免费视频内射| 精品国产乱码久久久久久小说| 亚洲av在线观看美女高潮| av福利片在线| 亚洲精品在线美女| 亚洲精品日韩在线中文字幕| 亚洲,一卡二卡三卡| 欧美黄色片欧美黄色片| 欧美少妇被猛烈插入视频| 日韩电影二区| 女人精品久久久久毛片| 国产乱人偷精品视频| 亚洲av男天堂| 18禁观看日本| 青草久久国产| 麻豆精品久久久久久蜜桃| 永久免费av网站大全| 国产精品久久久久成人av| 美女高潮到喷水免费观看| 精品一区在线观看国产| 91成人精品电影| 中文字幕色久视频| 亚洲国产毛片av蜜桃av| 国产成人91sexporn| 国产一级毛片在线| 成人毛片a级毛片在线播放| 2021少妇久久久久久久久久久| 丝袜美腿诱惑在线| 又黄又粗又硬又大视频| 最近中文字幕2019免费版| 亚洲成人手机| 精品酒店卫生间| 黄色毛片三级朝国网站| 亚洲国产日韩一区二区| 精品一区在线观看国产| 日韩视频在线欧美| 99久久中文字幕三级久久日本| 久久久久久久大尺度免费视频| 中文欧美无线码| 性色av一级| 91国产中文字幕| 精品国产超薄肉色丝袜足j| 免费高清在线观看视频在线观看| 亚洲精品成人av观看孕妇| 午夜日本视频在线| 免费黄色在线免费观看| 麻豆av在线久日| 一二三四在线观看免费中文在| 老鸭窝网址在线观看| 丝袜脚勾引网站| 成人毛片a级毛片在线播放| 久久久久精品性色| 国产精品 国内视频| 97精品久久久久久久久久精品| 国产精品 国内视频| 下体分泌物呈黄色| 母亲3免费完整高清在线观看 | 熟女av电影| 亚洲婷婷狠狠爱综合网| 久久免费观看电影| 99久久人妻综合| www.熟女人妻精品国产| 夫妻性生交免费视频一级片| 精品少妇一区二区三区视频日本电影 | 亚洲第一区二区三区不卡| 免费日韩欧美在线观看| 亚洲av欧美aⅴ国产| 欧美日韩视频高清一区二区三区二| 午夜日韩欧美国产| 成人午夜精彩视频在线观看| 精品久久蜜臀av无| 精品亚洲乱码少妇综合久久| 国产成人精品婷婷| 考比视频在线观看| av在线播放精品| 精品国产乱码久久久久久男人| 2021少妇久久久久久久久久久| 在线看a的网站| av女优亚洲男人天堂| 亚洲人成电影观看| 久久久久久久大尺度免费视频| 七月丁香在线播放| 男女午夜视频在线观看| av片东京热男人的天堂| 亚洲国产毛片av蜜桃av| 欧美精品国产亚洲| 午夜福利一区二区在线看| 亚洲国产精品一区三区| 老汉色av国产亚洲站长工具| av国产精品久久久久影院| 久久女婷五月综合色啪小说| 9191精品国产免费久久| 免费女性裸体啪啪无遮挡网站| 成年人午夜在线观看视频| 亚洲精品日韩在线中文字幕| 日韩一区二区三区影片| 午夜av观看不卡| 999久久久国产精品视频| 亚洲av男天堂| 久久国产精品男人的天堂亚洲| 大片电影免费在线观看免费| 亚洲人成网站在线观看播放| 天堂8中文在线网| 国产白丝娇喘喷水9色精品| 精品卡一卡二卡四卡免费| 制服丝袜香蕉在线| 欧美日韩精品网址| 免费在线观看黄色视频的| 亚洲人成网站在线观看播放| 午夜福利影视在线免费观看| 一区在线观看完整版| 久久这里有精品视频免费| 久久精品国产综合久久久| 两性夫妻黄色片| 天天影视国产精品| 久久精品久久精品一区二区三区| 一级片免费观看大全| av在线app专区| 黑人欧美特级aaaaaa片| 国产福利在线免费观看视频| 亚洲精品国产av蜜桃| 久久ye,这里只有精品| 免费观看在线日韩| 18禁观看日本| 亚洲国产精品一区三区| 男女免费视频国产| 丝瓜视频免费看黄片| 黑人欧美特级aaaaaa片| 边亲边吃奶的免费视频| 国产乱人偷精品视频| 人人妻人人添人人爽欧美一区卜| 亚洲av欧美aⅴ国产| 新久久久久国产一级毛片| 久久99热这里只频精品6学生| 亚洲精品aⅴ在线观看| 国产在线一区二区三区精| 久久热在线av| 天天操日日干夜夜撸| 男女午夜视频在线观看| 啦啦啦中文免费视频观看日本| 久久久久久久久久久久大奶| 精品亚洲成国产av| 久久 成人 亚洲| 国产精品偷伦视频观看了| h视频一区二区三区|