• <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有码第一页| 69精品国产乱码久久久| 中文字幕最新亚洲高清| av不卡在线播放| 成人av一区二区三区在线看| 亚洲中文av在线| 一级片免费观看大全| 人成视频在线观看免费观看| 中文字幕另类日韩欧美亚洲嫩草| 午夜免费观看网址| 亚洲男人天堂网一区| 久久久久久久精品吃奶| 最新在线观看一区二区三区| 免费在线观看视频国产中文字幕亚洲| 天天躁狠狠躁夜夜躁狠狠躁| а√天堂www在线а√下载 | 亚洲成人国产一区在线观看| 国产精华一区二区三区| 国产高清国产精品国产三级| 精品国产亚洲在线| a级片在线免费高清观看视频| 最近最新免费中文字幕在线| 国产精品一区二区在线不卡| 妹子高潮喷水视频| 日韩欧美一区视频在线观看| 一边摸一边抽搐一进一小说 | 成年人免费黄色播放视频| 黄片大片在线免费观看| 成人影院久久| 国产一区二区三区视频了| 亚洲,欧美精品.| 色精品久久人妻99蜜桃| 高清av免费在线| 热99国产精品久久久久久7| 色在线成人网| 国产成人精品久久二区二区91| 成年女人毛片免费观看观看9 | 国产麻豆69| 亚洲精品国产一区二区精华液| 又黄又粗又硬又大视频| 国产亚洲精品第一综合不卡| 大陆偷拍与自拍| 无人区码免费观看不卡| 91精品国产国语对白视频| 欧美 日韩 精品 国产| 热99久久久久精品小说推荐| 成人手机av| av一本久久久久| 一二三四在线观看免费中文在| 在线观看免费视频网站a站| 亚洲av第一区精品v没综合| aaaaa片日本免费| 看黄色毛片网站| 韩国精品一区二区三区| 国产三级黄色录像| 久99久视频精品免费| 亚洲色图综合在线观看| 欧美 日韩 精品 国产| 天堂√8在线中文| 99香蕉大伊视频| 一级毛片女人18水好多| 日韩中文字幕欧美一区二区| 国产一区二区三区综合在线观看| 99精品久久久久人妻精品| 看免费av毛片| 精品亚洲成a人片在线观看| 久久中文字幕人妻熟女| 少妇的丰满在线观看| 欧美 亚洲 国产 日韩一| 精品久久久久久,| а√天堂www在线а√下载 | 久久久久精品人妻al黑| 午夜精品久久久久久毛片777| 欧美最黄视频在线播放免费 | 精品视频人人做人人爽| 欧美精品人与动牲交sv欧美| 成人影院久久| 自线自在国产av| 日韩熟女老妇一区二区性免费视频| 美女 人体艺术 gogo| 男男h啪啪无遮挡| 亚洲全国av大片| 露出奶头的视频| 国产成人免费观看mmmm| av不卡在线播放| 午夜福利乱码中文字幕| 又大又爽又粗| 脱女人内裤的视频| 12—13女人毛片做爰片一| 国产精品九九99| 无人区码免费观看不卡| 国产免费男女视频| 成年动漫av网址| 久久精品成人免费网站| 亚洲情色 制服丝袜| 精品国产亚洲在线| 国产午夜精品久久久久久| av福利片在线| 欧美性长视频在线观看| 国产精品香港三级国产av潘金莲| 成人手机av| 99精品在免费线老司机午夜| 国产精品久久久av美女十八| 大型av网站在线播放| 欧美中文综合在线视频| 国产亚洲av高清不卡| 一本一本久久a久久精品综合妖精| 69精品国产乱码久久久| 一级毛片高清免费大全| 精品亚洲成国产av| 国产成人精品久久二区二区91| 伦理电影免费视频| 国产一区二区三区视频了| 岛国毛片在线播放| 亚洲成a人片在线一区二区| 欧美人与性动交α欧美软件| 亚洲男人天堂网一区| 一进一出抽搐动态| 一本大道久久a久久精品| 国产淫语在线视频| 桃红色精品国产亚洲av| 看片在线看免费视频| 新久久久久国产一级毛片| 亚洲一码二码三码区别大吗| 50天的宝宝边吃奶边哭怎么回事| 精品亚洲成国产av| 一级片免费观看大全| √禁漫天堂资源中文www| 热re99久久国产66热| 高清欧美精品videossex| 99久久人妻综合| 香蕉丝袜av| 国产男女内射视频| 亚洲精品在线美女| 午夜福利免费观看在线| 成人三级做爰电影| 中文亚洲av片在线观看爽 | 国产免费男女视频| 国产精品久久久人人做人人爽| 久久精品亚洲av国产电影网| 天天躁日日躁夜夜躁夜夜| 啦啦啦免费观看视频1| 日韩大码丰满熟妇| 久久 成人 亚洲| 精品久久久久久,| 最近最新免费中文字幕在线| 亚洲精品av麻豆狂野| 亚洲少妇的诱惑av| 欧美精品人与动牲交sv欧美| 亚洲精品久久午夜乱码| 色在线成人网| 一区在线观看完整版| 欧美最黄视频在线播放免费 | 久久精品国产清高在天天线| 国产欧美日韩综合在线一区二区| 国产精品1区2区在线观看. | xxxhd国产人妻xxx| 色在线成人网| 亚洲一卡2卡3卡4卡5卡精品中文| 激情在线观看视频在线高清 | 一本大道久久a久久精品| 少妇被粗大的猛进出69影院| 精品久久久久久电影网| 欧美不卡视频在线免费观看 | 欧美乱妇无乱码| 国产成人精品久久二区二区免费| 久久国产精品男人的天堂亚洲| 最新的欧美精品一区二区| 啦啦啦在线免费观看视频4| 在线观看一区二区三区激情| 日本一区二区免费在线视频| 国产亚洲精品久久久久久毛片 | 丰满饥渴人妻一区二区三| 国产一区二区激情短视频| 中文字幕色久视频| 久久久精品区二区三区| 精品国产一区二区三区久久久樱花| 国产激情久久老熟女| 侵犯人妻中文字幕一二三四区| 美女高潮到喷水免费观看| 亚洲综合色网址| 999精品在线视频| 狠狠婷婷综合久久久久久88av| 久久影院123| 亚洲av成人不卡在线观看播放网| 亚洲人成电影观看| 91九色精品人成在线观看| 久久久精品免费免费高清| 日本一区二区免费在线视频| 亚洲三区欧美一区| 久久人妻福利社区极品人妻图片| 怎么达到女性高潮| 国产真人三级小视频在线观看| 国产av又大| 欧美 亚洲 国产 日韩一| 99精品欧美一区二区三区四区| 国产区一区二久久| 欧美日韩黄片免| 免费观看人在逋| 最近最新中文字幕大全电影3 | 欧美日韩福利视频一区二区| 激情在线观看视频在线高清 | 欧美国产精品va在线观看不卡| 精品乱码久久久久久99久播| 久久久国产成人精品二区 | 午夜福利视频在线观看免费| 日韩人妻精品一区2区三区| 人人澡人人妻人| 亚洲欧美激情在线| 欧美日韩中文字幕国产精品一区二区三区 | 丝瓜视频免费看黄片| 免费人成视频x8x8入口观看| 一本综合久久免费| 老汉色av国产亚洲站长工具| 搡老乐熟女国产| 精品亚洲成国产av| 久9热在线精品视频| 久久久久国产精品人妻aⅴ院 | 中文字幕精品免费在线观看视频| 久久香蕉国产精品| 免费av中文字幕在线| 欧美+亚洲+日韩+国产| 人成视频在线观看免费观看| 婷婷丁香在线五月| www日本在线高清视频| 日本五十路高清| 一二三四社区在线视频社区8| 在线观看免费日韩欧美大片| 欧美乱码精品一区二区三区| 五月开心婷婷网| 亚洲精品国产色婷婷电影| 精品人妻1区二区| 国产精品成人在线| 欧美精品啪啪一区二区三区| 亚洲av电影在线进入| 国产高清国产精品国产三级| 国产主播在线观看一区二区| 夜夜夜夜夜久久久久| 日本黄色视频三级网站网址 | 久久久精品区二区三区| 大型av网站在线播放| 美女福利国产在线| 高清视频免费观看一区二区| 人妻久久中文字幕网| 天天躁夜夜躁狠狠躁躁| a在线观看视频网站| 国产精品久久久久成人av| 亚洲,欧美精品.| 久久精品国产99精品国产亚洲性色 | 一进一出抽搐gif免费好疼 | 久久精品人人爽人人爽视色| 热99re8久久精品国产| 色播在线永久视频| tube8黄色片| 最近最新中文字幕大全免费视频| 国产成人啪精品午夜网站| 黄色成人免费大全| 天堂动漫精品| 国产精品偷伦视频观看了| 久久 成人 亚洲| 国产亚洲精品第一综合不卡| 一夜夜www| 欧洲精品卡2卡3卡4卡5卡区| 一边摸一边做爽爽视频免费| 欧美午夜高清在线| 国产欧美日韩一区二区精品| 亚洲午夜理论影院| 亚洲国产毛片av蜜桃av| 老鸭窝网址在线观看| 久久午夜综合久久蜜桃| 免费在线观看影片大全网站| 天天添夜夜摸| 中文字幕色久视频| 亚洲三区欧美一区| 国产一区二区三区在线臀色熟女 | 亚洲在线自拍视频| 两个人免费观看高清视频| 成人国语在线视频| 韩国精品一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 激情视频va一区二区三区| 欧美激情高清一区二区三区| 夫妻午夜视频| 成人av一区二区三区在线看| 日韩视频一区二区在线观看| 777米奇影视久久| 久久久久久免费高清国产稀缺| 国产一区二区激情短视频| 国产亚洲欧美精品永久| 一级a爱视频在线免费观看| 十八禁网站免费在线| 国产精品亚洲av一区麻豆| 午夜激情av网站| 国产蜜桃级精品一区二区三区 | 97人妻天天添夜夜摸| 国产不卡av网站在线观看| 在线观看免费日韩欧美大片| 99国产精品99久久久久| 美女午夜性视频免费| 亚洲一区二区三区欧美精品| 久久午夜综合久久蜜桃| 欧美日韩亚洲高清精品| 成人亚洲精品一区在线观看| 中文字幕色久视频| 国产精品98久久久久久宅男小说| 99热只有精品国产| 精品国产一区二区久久| 国产精品久久久久成人av| 黄色女人牲交| 大码成人一级视频| 国产亚洲精品第一综合不卡| 一边摸一边做爽爽视频免费| 高清黄色对白视频在线免费看| 韩国av一区二区三区四区| 亚洲精品粉嫩美女一区| 亚洲在线自拍视频| 亚洲av美国av| 精品少妇久久久久久888优播| 欧美午夜高清在线| svipshipincom国产片| 午夜福利欧美成人| 十分钟在线观看高清视频www| 成人av一区二区三区在线看| 日韩有码中文字幕| 中文字幕制服av| 成人黄色视频免费在线看| 亚洲三区欧美一区| 一区二区三区精品91| 侵犯人妻中文字幕一二三四区| cao死你这个sao货| 色老头精品视频在线观看| 大陆偷拍与自拍| 国产精品一区二区在线不卡| 亚洲少妇的诱惑av| 国产成人影院久久av| 久久国产亚洲av麻豆专区| 丝瓜视频免费看黄片| 欧美精品人与动牲交sv欧美| 免费一级毛片在线播放高清视频 | av网站在线播放免费| 国产无遮挡羞羞视频在线观看| 大香蕉久久网| 脱女人内裤的视频| 免费看a级黄色片| 亚洲五月天丁香| 日日爽夜夜爽网站| 国产日韩欧美亚洲二区| 麻豆av在线久日| 亚洲 欧美一区二区三区| av天堂在线播放| 1024视频免费在线观看| 亚洲色图综合在线观看| 国产又色又爽无遮挡免费看| 国产精品成人在线| 91九色精品人成在线观看| 热99久久久久精品小说推荐| bbb黄色大片| 午夜激情av网站| 亚洲av电影在线进入| 看免费av毛片| 日韩熟女老妇一区二区性免费视频| 亚洲精品av麻豆狂野| 国产三级黄色录像| 亚洲午夜理论影院| 少妇裸体淫交视频免费看高清 | 欧美 亚洲 国产 日韩一| 国产伦人伦偷精品视频| 久久久久国产精品人妻aⅴ院 | 美女视频免费永久观看网站| 18禁黄网站禁片午夜丰满| 男人舔女人的私密视频| 色婷婷av一区二区三区视频| 一级a爱片免费观看的视频| 亚洲精品乱久久久久久| av在线播放免费不卡| 亚洲精品国产一区二区精华液| 精品一品国产午夜福利视频| 伦理电影免费视频| 人人妻人人爽人人添夜夜欢视频| 亚洲在线自拍视频| 国产精品九九99| 看黄色毛片网站| 成在线人永久免费视频| 亚洲一区二区三区不卡视频| 亚洲av熟女| 欧美日韩精品网址| 亚洲av日韩在线播放| 91精品国产国语对白视频| 啦啦啦免费观看视频1| 中文字幕人妻熟女乱码| 9热在线视频观看99| 日韩熟女老妇一区二区性免费视频| 高清av免费在线| 国产高清视频在线播放一区| 黄频高清免费视频| 国产精品电影一区二区三区 | 久久久精品国产亚洲av高清涩受| 在线十欧美十亚洲十日本专区| 欧美午夜高清在线| 欧美黄色淫秽网站| 欧美 日韩 精品 国产| 丝袜美腿诱惑在线| 精品一区二区三区av网在线观看| www.自偷自拍.com| 久久久久国内视频| 欧美日韩国产mv在线观看视频| 两个人免费观看高清视频| 国产国语露脸激情在线看| 一级a爱视频在线免费观看| 俄罗斯特黄特色一大片| 亚洲精品国产色婷婷电影| 亚洲国产精品sss在线观看 | 亚洲av熟女| 两性午夜刺激爽爽歪歪视频在线观看 | 99riav亚洲国产免费| 欧美日韩黄片免| 在线观看免费午夜福利视频| 涩涩av久久男人的天堂| 国产激情久久老熟女| 91麻豆精品激情在线观看国产 | 午夜视频精品福利| 女人精品久久久久毛片| 久热这里只有精品99| 两个人看的免费小视频| videosex国产| 欧美国产精品va在线观看不卡| 电影成人av| 麻豆乱淫一区二区| 99香蕉大伊视频| av中文乱码字幕在线| 人妻久久中文字幕网| 他把我摸到了高潮在线观看| 窝窝影院91人妻| 日韩视频一区二区在线观看| 欧美日韩av久久| 亚洲五月婷婷丁香| 99久久人妻综合| 自拍欧美九色日韩亚洲蝌蚪91| 国产高清激情床上av| 在线观看舔阴道视频| 精品国产国语对白av| 国内毛片毛片毛片毛片毛片| 国产精品1区2区在线观看. | 中文亚洲av片在线观看爽 | 午夜福利欧美成人| 久久性视频一级片| 三上悠亚av全集在线观看| 女人精品久久久久毛片| 在线永久观看黄色视频| 啦啦啦视频在线资源免费观看| 99精国产麻豆久久婷婷| 水蜜桃什么品种好| av福利片在线| 搡老岳熟女国产| 久久ye,这里只有精品| 黄色视频不卡| 欧美 亚洲 国产 日韩一| 亚洲国产欧美一区二区综合| 成人av一区二区三区在线看| 国产单亲对白刺激| 午夜两性在线视频| 精品无人区乱码1区二区| 搡老岳熟女国产| 久久影院123| 久久久国产精品麻豆| 色婷婷久久久亚洲欧美| 欧美精品一区二区免费开放| 亚洲少妇的诱惑av| 黄色丝袜av网址大全| 久久久水蜜桃国产精品网| 国产在视频线精品| 这个男人来自地球电影免费观看| 自线自在国产av| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲成人国产一区在线观看| 亚洲av片天天在线观看| 最近最新免费中文字幕在线| 999久久久国产精品视频| 欧美一级毛片孕妇| 人人妻人人爽人人添夜夜欢视频| 女人爽到高潮嗷嗷叫在线视频| 亚洲精华国产精华精| 午夜两性在线视频| 99久久精品国产亚洲精品| 久久亚洲精品不卡| 国产精品电影一区二区三区 | 在线观看www视频免费| 美女 人体艺术 gogo| 成年动漫av网址| 午夜免费鲁丝| 18禁美女被吸乳视频| 久久久精品区二区三区| 久久久久久亚洲精品国产蜜桃av| 18禁国产床啪视频网站| 韩国av一区二区三区四区| 久久天躁狠狠躁夜夜2o2o| 欧美 日韩 精品 国产| 老司机午夜福利在线观看视频| 国产91精品成人一区二区三区| 91精品国产国语对白视频| 欧美在线一区亚洲| 国产欧美日韩一区二区三区在线| 欧美另类亚洲清纯唯美| 99在线人妻在线中文字幕 | 制服诱惑二区| 欧美av亚洲av综合av国产av| 精品一区二区三区av网在线观看| 国产伦人伦偷精品视频| 国产一区二区三区视频了| 久久精品国产99精品国产亚洲性色 | 国产精品电影一区二区三区 | 久久婷婷成人综合色麻豆| 欧美黄色淫秽网站| 男女床上黄色一级片免费看| 无限看片的www在线观看| 国产精品偷伦视频观看了| 亚洲美女黄片视频| 人人澡人人妻人| 欧美日韩亚洲综合一区二区三区_| 国产精品亚洲一级av第二区| 日本vs欧美在线观看视频| 亚洲视频免费观看视频| 91九色精品人成在线观看| 18禁裸乳无遮挡免费网站照片 | 日韩有码中文字幕| 久久亚洲精品不卡| 成人精品一区二区免费| 午夜福利在线观看吧| 99久久人妻综合| 国产高清videossex| 亚洲av电影在线进入| 麻豆成人av在线观看| 99热网站在线观看| 黄片小视频在线播放| 女人爽到高潮嗷嗷叫在线视频| 黄片大片在线免费观看| 亚洲国产毛片av蜜桃av| 久久青草综合色| 精品福利永久在线观看| 国产深夜福利视频在线观看| svipshipincom国产片| 午夜福利在线免费观看网站| 精品少妇一区二区三区视频日本电影| 免费不卡黄色视频| 亚洲男人天堂网一区| 99国产精品一区二区三区| 久9热在线精品视频| 久久精品熟女亚洲av麻豆精品| 啦啦啦视频在线资源免费观看| 母亲3免费完整高清在线观看| 丝瓜视频免费看黄片| 亚洲欧洲精品一区二区精品久久久| 国产精品自产拍在线观看55亚洲 | 久久久精品国产亚洲av高清涩受| 他把我摸到了高潮在线观看| 女人爽到高潮嗷嗷叫在线视频| 成人手机av| 欧美在线黄色| 黄色毛片三级朝国网站| av超薄肉色丝袜交足视频| 国产在线一区二区三区精| 精品第一国产精品| 满18在线观看网站| 国产欧美日韩一区二区三区在线| 丰满人妻熟妇乱又伦精品不卡| 久久国产精品大桥未久av| 日韩中文字幕欧美一区二区| 777米奇影视久久| 成人免费观看视频高清| 国产亚洲av高清不卡| 国产成人欧美在线观看 | 欧美人与性动交α欧美软件| 午夜精品久久久久久毛片777| 一区二区日韩欧美中文字幕| 精品国产国语对白av| 成年人黄色毛片网站| 日韩精品免费视频一区二区三区| 色综合婷婷激情| 超色免费av| 99re在线观看精品视频| 国产成人精品久久二区二区免费| 国产精品偷伦视频观看了| 日韩 欧美 亚洲 中文字幕| 热re99久久国产66热| 久久久久国内视频| 好男人电影高清在线观看| 国产精品免费大片| 99国产综合亚洲精品| 亚洲精品成人av观看孕妇| 欧美黄色淫秽网站| 国产三级黄色录像| 日韩中文字幕欧美一区二区| 一级片'在线观看视频| 一个人免费在线观看的高清视频| 国产亚洲一区二区精品| 亚洲avbb在线观看| 午夜福利在线免费观看网站| 999久久久精品免费观看国产| 法律面前人人平等表现在哪些方面| 国产av精品麻豆| 飞空精品影院首页| 欧美精品av麻豆av| 99精国产麻豆久久婷婷| 欧美 日韩 精品 国产| 搡老乐熟女国产|