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

    調(diào)制信號雙譜在滾動軸承故障診斷中的應(yīng)用研究

    2021-06-02 08:59:52王建國谷豐收
    關(guān)鍵詞:雙譜邊帶譜分析

    張 超, 張 輝, 辛 闊, 王建國, 谷豐收

    (1. 內(nèi)蒙古科技大學(xué) 機(jī)械工程學(xué)院,內(nèi)蒙古 包頭 014010;2. 內(nèi)蒙古自治區(qū)機(jī)電系統(tǒng)智能診斷與控制重點實驗室,內(nèi)蒙古 包頭 014010;3. 哈德斯菲爾德大學(xué)效率與性能工程中心,英國哈德斯菲爾德HD1 3DH)

    引言

    在旋轉(zhuǎn)機(jī)械中,滾動軸承是支撐旋轉(zhuǎn)軸的重要零部件,并且在現(xiàn)代工業(yè)中,滾動軸承也得到廣泛的使用,其工作狀態(tài)的優(yōu)劣直接影響著整個機(jī)器設(shè)備運行的穩(wěn)定性和安全性[1-2]。因此,如果能及早的通過信號處理分析方法判別滾動軸承已產(chǎn)生的微小裂紋,將具有重要意義[3]。

    如今,分析軸承故障的常用手段仍然是基于振動信號分析。當(dāng)前,已經(jīng)開發(fā)了許多基于振動特征分析的軸承故障檢測的信號處理方法,例如Wigner-Viller分布(Wigner villedistribution, WVD)[4],經(jīng)驗小波變換(Empirical wavelet transform, EWT)[5],互補(bǔ)集合經(jīng)驗?zāi)B(tài)分解(Complementary ensemble empirical mode decomposition, EEMD)[6],總體局部均值分解(Ensemble local means decomposition,ELMD)[7]及變分模式分解(variational mode decomposition,VMD)[8]等。盡管這些方法可以有效地檢測滾動軸承的故障,但仍然存在一些局限性。此外,還有多種方法已用于軸承元件的故障診斷,但是這些方法大多數(shù)集中在降低噪聲上,而忽略信號在采集時就發(fā)生的調(diào)制現(xiàn)象,未對信號進(jìn)行解調(diào)處理,都可能降低診斷的精確性。

    對信號進(jìn)行頻譜分析是提取故障特征頻率的主要方法,例如:功率譜分析,雙譜分析。其中功率譜分析忽略了信號的相位信息,僅包含單個分量幅度信息,且不能解調(diào)調(diào)制信號;雙譜分析只考慮信號較高邊帶的信息,沒有考慮低邊帶信息。這些方法在提取故障特征時都有不足的地方,并且都不能充分的解調(diào)信號中的調(diào)制成分。

    基于上述原因,利用具有良好的信號解調(diào)特性的調(diào)制信號雙譜分析方法解調(diào)信號中的調(diào)制現(xiàn)象,并且調(diào)制信號雙譜有良好的降噪效果及對干擾頻率有很好的抑制作用。提出了調(diào)制信號雙譜分析方法并用于滾動軸承故障診斷。

    調(diào)制信號雙譜是一種基于傳統(tǒng)雙譜改進(jìn)的高級信號解調(diào)方法。作為一種考慮能上下兩個邊頻帶的雙譜分析方法,能對調(diào)制信號進(jìn)行有效分析,該方法可以較好地抑制隨機(jī)噪聲和非周期成分的干擾,反映出信號中的調(diào)制成分。自谷豐收等人在2011年提出了調(diào)制信號雙譜分析方法,并在下游機(jī)械設(shè)備的故障診斷中驗證這種方法的準(zhǔn)確性[9]。隨后,通過應(yīng)用調(diào)制信號雙譜分析方法來診斷轉(zhuǎn)子斷條的故障,得到好的診斷結(jié)果[10]。

    本文主要研究調(diào)制信號雙譜分析方法在滾動軸承故障診斷中的應(yīng)用。首先將采集到的振動信號,使用MSB方法進(jìn)行降噪和解調(diào),得到調(diào)制信號雙譜圖;但由于為準(zhǔn)確的量化邊帶幅度,使用調(diào)制信號雙譜?邊帶估計器;通過選擇幾個次優(yōu)的切片構(gòu)成最佳的切片;最終,通過復(fù)合切片譜實現(xiàn)故障特征頻率的提取,為進(jìn)一步驗證MSB的準(zhǔn)確性,并用調(diào)制信號雙譜相干性檢測。最終,通過對模擬信號及滾動軸承故障實驗數(shù)據(jù)進(jìn)行分析,來驗證該方法的有效性和優(yōu)越性。

    1 調(diào)制信號雙譜分析方法

    1.1 常規(guī)雙譜

    給定一離散的振動信號x(t),其離散傅里葉變換為[11]:

    其中X(f)是復(fù)數(shù),幅值|X(f)|和相位φf的表達(dá)為:

    通過傅里葉變換,雙譜B(fc,fx)在頻域內(nèi)可表示為:

    其中X*是X復(fù)共軛數(shù),E<>是數(shù)學(xué)期望表達(dá)式,fc,fx和fc+fx是3個不同的頻率。與二階累積量不同,常規(guī)雙譜的三級累積量是復(fù)數(shù),既包含原始信號的幅值信息又包含相位信息。如fc,fx和fc+fx處的3個頻率分別是獨立分量,每個頻率分量的相位符合均勻隨機(jī)分布規(guī)律,并分布在(-π,π)上。計算統(tǒng)計期望,由于隨機(jī)相位的分布,其雙譜將趨于零,因此應(yīng)用這種方法能夠抑制隨機(jī)噪聲。相反如果fc,fx,fc+fx3個頻率分量相互產(chǎn)生非線性耦合,其相位存在特定關(guān)系,即使3個頻率各自是隨機(jī)的,相位也符合式關(guān)系:

    因此,雙譜統(tǒng)計平均值將不會趨向于零。這種非線性耦合在雙譜分析中表現(xiàn)為雙頻B(fc,fx)的一個峰值。

    1.2 調(diào)制信號雙譜

    調(diào)制信號雙譜是一種考慮邊頻帶的雙譜分析方法,克服了常規(guī)雙譜在處理調(diào)制信號中的不足,能對調(diào)制信號進(jìn)行有效分析,可以很好地抑制隨機(jī)噪聲、非周期成分的干擾以及清晰反映出信號中的調(diào)制成分[12-13]。

    對于具有相應(yīng)離散傅里葉變換x(f)的離散時間信號x(t)的調(diào)制信號雙譜可被定義為:

    其中,BMS(fc,fx)為信號x(t)的雙譜,fc為載波頻率,fx為調(diào)制頻率,X*(fc)是X(fc)的復(fù)共軛,E<>表示統(tǒng)計期望,fc+fx為上邊帶頻率,fc?fx為下邊帶頻率。

    其中式5中調(diào)制信號雙譜的全部相位信息可表示為:

    為了更準(zhǔn)確的闡述MSB的有效性,使用調(diào)制信號雙譜相關(guān)性來評估成分之間的耦合程度[14],如下式所示:

    為了準(zhǔn)確的量化邊帶幅度,通過消除實質(zhì)影響來修改載波頻率fc的分量,為此對MSB歸一化處理,歸一化后的調(diào)制信號雙譜定義為調(diào)制信號雙譜邊帶估計器:

    BMS(fc,0)為fx= 0時的平方功率譜。

    在利用調(diào)制信號雙譜邊帶估計器進(jìn)行分析時,即必須選取出合適的fc切片位置。為此計算fx增量方向上MSB的平均值,得到fc切片。

    其中Δf是fx的分辨率。

    但是,故障特征較為明顯的切片位置有多個,選取多個切片可以綜合各切片處的故障特征信息,同時也可以減小存在于單個切片中的干擾成分。所以計算多個最優(yōu)切片的平均值,即復(fù)合切片譜,可以通過下式表示:

    N為所取切片的總數(shù)。

    2 故障特征提取流程

    為了減少噪聲的干擾及解決雙譜分析沒有考慮低邊帶的問題。研究了調(diào)制信號雙譜的滾動軸承故障特征提取方法,集體步驟如下,流程圖如圖1所示:

    圖1 診斷流程圖

    步驟1:利用加速度傳感器從內(nèi)蒙古科技大學(xué)機(jī)械工程學(xué)院的HZXT-DS-003雙跨雙轉(zhuǎn)子試驗臺采集軸承的外圈故障數(shù)據(jù);

    步驟2:對采集到的信號利用調(diào)制信號雙譜方法進(jìn)行分析;

    步驟3:為了更準(zhǔn)確的量化邊帶幅度,對MSB分析結(jié)果進(jìn)行歸一化處理;

    步驟4:計算MSB的fc切片,獲取幾個次優(yōu)的fc切片;

    步驟5:將次優(yōu)的fc切片構(gòu)成復(fù)合切片譜,提取故障特征頻率。

    3 仿真信號驗證

    本節(jié)利用仿真信號驗證所提方法的正確性和有效性。該信號包括3個典型過程:相位調(diào)制,幅值調(diào)制和相位耦合。

    本文的仿真信號為:

    其中,xm為相位調(diào)制成分,xm1l和xm1r分別為下邊帶幅值調(diào)制成分和上邊帶幅值調(diào)制成分,xma1為相位耦合成分。其中,載波頻率為fc0= 50 Hz,幅度為1,相位調(diào)制fpm= 24 Hz,其幅值為a,其中相位af= 0.5π,這意味著該相位調(diào)制過程是諧波相位耦合過程,因為頻率分量是諧波,并且相位通過求和或差相連接,還有非諧波相位耦合過程fa1= 10 Hz。幅值調(diào)制函數(shù)的幅值為0.2 a,其頻率fa2= 14 Hz,此外,還有一個附加的相位過程,其振幅為a,頻率仍為fa1且具有相同的相位af。另外,n(t)是添加到信號中的?1 dB的高斯噪聲,以顯示調(diào)制信號雙譜的噪聲抑制能力。

    當(dāng)不填加噪聲時,該仿真信號的時頻圖如圖2所示,圖2僅顯示了仿真信號的一部分。該信號具有409 600個采樣點,其采樣頻率為512 Hz,轉(zhuǎn)換到頻域中,在50 Hz處有明顯的峰值,此為載波頻率所對應(yīng)的峰值,也可以看出模擬信號發(fā)生調(diào)制;在10 Hz處出現(xiàn)諧波相位耦合過程;另外在26、74 Hz處,有相對高的峰值,表示相位調(diào)制過程產(chǎn)生的上下邊頻帶;在36、64 Hz處,有相對低的峰值,表示幅值調(diào)制過程產(chǎn)生的上下邊頻帶。

    圖2 仿真信號時頻圖(無噪聲)

    圖3為添加噪聲的仿真信號的時頻圖,圖3僅顯示了全部仿真信號的一部分,仿真信號同樣具有409 600個采樣點,其采樣率為512 Hz。首先觀察信號的時域圖,通過與未加噪聲的信號對比,可知時域信號中周期性的特征已完全被噪聲淹沒。然后轉(zhuǎn)到頻域圖中,只能看到在50 Hz處有明顯的峰值,其余的邊帶信息已被干擾頻率覆蓋。

    圖3 仿真信號時頻圖(含噪聲)

    為了驗證仿真信號的有效性,接下來使用功率譜分析方法對未加噪的仿真信號及帶噪信號進(jìn)行分析如圖4所示。

    圖4 功率譜圖

    由圖4可知,在無噪聲的情況下可以清晰的辨別信號中已發(fā)生調(diào)制信息產(chǎn)生的邊頻帶。當(dāng)加入噪聲后,在功率譜中仍可以觀察到調(diào)制分量產(chǎn)生邊帶信息。其中在50 ± 24 Hz處,表示相位調(diào)制過程;在50 ± 14 Hz處,表示幅值調(diào)制過程;另外在10 Hz處,為相位耦合過程。從圖中可知在功率譜中無法解調(diào)相位調(diào)制過程及相位耦合產(chǎn)生的高次諧波,并且對應(yīng)于不同頻率的幅度高于無噪聲功率譜的幅度。

    接下來,使用常規(guī)雙譜分析方法對加噪信號信息分析,如圖5所示,常規(guī)雙譜只在(50,50) Hz處出現(xiàn)唯一峰值,在(14,50)、(24,50) Hz處沒有出現(xiàn)明顯得調(diào)制成分,表明調(diào)制信號不能被解調(diào)。

    圖5 常規(guī)雙譜圖分析結(jié)果

    接下來,使用本文所提方法對填加噪聲的信號進(jìn)行分析,如圖6所示,其中在(10,50) Hz處沒有出現(xiàn)高的幅值,表明相位耦合過程被解調(diào),另外在(26,50)、(74,50)、(36,50)及(64,50) Hz處,沒有出現(xiàn)較高的峰值,說明由相位調(diào)制過程及幅值調(diào)制過程所產(chǎn)生的上下邊頻帶被解調(diào)。另外,在(14,50)、(24,50) Hz處出現(xiàn)明顯的峰值,分別為幅值和相位特征成份,可知信號被充分解調(diào)。

    圖6 調(diào)制信號雙譜分析結(jié)果

    為了進(jìn)一步驗證調(diào)制信號雙譜的準(zhǔn)確性,計算調(diào)制信號雙譜相干性去檢測原始信號中相位調(diào)制成分,相位耦合成分及幅值調(diào)制成分之間的耦合程度。從圖7中可以看出在(14,50)、(24,50) Hz處出現(xiàn)明顯的峰值,顯然調(diào)制信號雙相干性的結(jié)果可以有效的表明載波信號和調(diào)制信號之間的非線性作用。到此階段,成功的分離信號中的調(diào)制成分,找到相位調(diào)制成分和幅值調(diào)制成分。

    圖7 調(diào)制信號雙譜相關(guān)性分析結(jié)果

    4 實驗信號驗證

    為了近一步驗證此方法的有效性,對實際采集到的振動信號進(jìn)行調(diào)制信號雙譜分析,從而提取軸承故障的特征。該振動信號的采集源自于內(nèi)蒙古科技大學(xué)機(jī)械工程學(xué)院的HZXT-DS-003雙跨雙轉(zhuǎn)子試驗臺,如圖8所示。

    圖8 實驗裝置

    采集數(shù)據(jù)所用滾動軸承的型號為6205-2RS,其結(jié)構(gòu)參數(shù)如表1所示。

    表1 6205-2RS軸承的結(jié)構(gòu)參數(shù)

    滾動軸承內(nèi)圈和外圈溝道上分別用電火花加工直徑為0.178 mm的單點損傷,滾動軸承故障件如圖9所示。

    圖9 軸承故障件

    本次所采集的數(shù)據(jù)是用加速度傳感器從試驗臺的軸承座上采集,軸的轉(zhuǎn)速為1000 r/min,采樣頻率為12 000 Hz,選用時長為1 s的數(shù)據(jù)進(jìn)行分析。其故障特征頻率及轉(zhuǎn)頻如表2所示。其振動信號的時頻圖如圖10所示。

    表2 6205-2RS軸承的故障特征頻率及轉(zhuǎn)頻

    圖10 實驗信號時頻圖

    從圖10中可以看出,低頻信號的完全被淹沒,而高頻信號幅值明顯,經(jīng)計算所得,故障軸承的特征頻率應(yīng)該在低頻段,出現(xiàn)上述這一現(xiàn)象的原因是由于受到噪聲和沖擊引起的高頻共振信號??芍苯佑妙l譜分析的診斷效果非常不理想。

    首先使用對信號進(jìn)行功率譜分析,其頻譜圖如圖11所示,可以看出有很多干擾頻率,故障特征信息被干擾頻率嚴(yán)重影響,無法準(zhǔn)確的提取信號故障特征。因此基于功率譜的方法在提取故障頻率時會得到不準(zhǔn)確的結(jié)果。接下來對重構(gòu)的信號進(jìn)行MSB分析,為了獲取更可靠的結(jié)果,通過式(9),得到MSB的fc切片如圖12所示。選取幾個次優(yōu)MSB切片的平均值來改進(jìn)MSB的診斷結(jié)果,故選取455、617及677 Hz處的切片,構(gòu)造復(fù)合切片譜,從而提取故障特征頻率。如圖13所示,從中可以清楚的看到在60 Hz處出現(xiàn)明顯的峰值,與滾動軸承的外圈故障頻率非常接近,且存在二倍頻特征,特征頻率明顯,沒有出現(xiàn)調(diào)制現(xiàn)象,有效的提取故障特征頻率。

    圖11 功率譜圖

    圖12 MSB的fc切片

    圖13 MSB診斷結(jié)果圖

    為了驗證調(diào)制信號雙譜分析在故障特征提取的準(zhǔn)確性,使用常規(guī)雙譜分析方法對信號進(jìn)行處理如圖14所示,在43 Hz附近出現(xiàn)有外圈故障頻率和轉(zhuǎn)頻發(fā)生調(diào)制產(chǎn)生的邊帶信息,其幅值高于軸承外圈故障頻率。在104、135 Hz處出現(xiàn)二倍頻與轉(zhuǎn)頻發(fā)生調(diào)制現(xiàn)象產(chǎn)生的上下邊帶信息。說明信號發(fā)生了調(diào)制現(xiàn)象,使用常規(guī)雙譜不能充分解調(diào)信號中調(diào)制現(xiàn)象,影響故障特征信息的提取。

    圖14 常規(guī)雙譜診斷結(jié)果圖

    5 結(jié)論

    1)調(diào)制信號雙譜是一種考慮邊頻帶的雙譜分析方法,能對調(diào)制信號進(jìn)行有效分析,可以很好地抑制隨機(jī)噪聲、非周期成分的干擾以及清晰反映出信號中的調(diào)制成分。

    2)通過選取多個次優(yōu)的fc切片,可以結(jié)合多個切片處的故障特征信息,能更加準(zhǔn)確的提取故障特征信息。

    3)與常規(guī)雙譜進(jìn)行對比,通過分析得知,調(diào)制信號雙譜在抑制干擾頻率及解調(diào)調(diào)制信號由于常規(guī)雙譜。

    4)從模擬信號分析到實驗數(shù)據(jù)驗證,都證明此方法的有效性和準(zhǔn)確性。

    猜你喜歡
    雙譜邊帶譜分析
    納譜分析技術(shù)(蘇州)有限公司
    色譜(2022年5期)2022-04-28 02:49:10
    光晶格中鐿原子的拉曼邊帶冷卻
    雙譜圖在語音分析中的應(yīng)用
    Cr12MoV冷作模具鋼滲鉻層界面能譜分析
    基于雙譜特征融合的通信輻射源識別算法
    Rotenberg模型中一類遷移算子的譜分析
    沉香GC-MS指紋圖譜分析
    中成藥(2016年8期)2016-05-17 06:08:26
    一種HEVC樣點自適應(yīng)補(bǔ)償改進(jìn)方法
    基于小波包域雙譜的風(fēng)力機(jī)振動信號監(jiān)測方法研究
    基于邊帶相關(guān)置換的BDS抗窄帶干擾算法
    18禁国产床啪视频网站| 国产蜜桃级精品一区二区三区| 1024视频免费在线观看| 国产99白浆流出| 搡老乐熟女国产| 叶爱在线成人免费视频播放| 国产成人欧美| 亚洲自偷自拍图片 自拍| 久久久久久久久免费视频了| 黄色成人免费大全| 国产精品一区二区免费欧美| 国产97色在线日韩免费| 久99久视频精品免费| 国产av精品麻豆| 一个人免费在线观看的高清视频| 另类亚洲欧美激情| 在线观看一区二区三区激情| 岛国在线观看网站| 成人三级做爰电影| 男人的好看免费观看在线视频 | 中文字幕高清在线视频| 成人三级做爰电影| bbb黄色大片| 久久久久九九精品影院| 天天添夜夜摸| 在线免费观看的www视频| 日韩精品免费视频一区二区三区| 天天躁夜夜躁狠狠躁躁| 波多野结衣av一区二区av| 女生性感内裤真人,穿戴方法视频| 日韩人妻精品一区2区三区| 日本撒尿小便嘘嘘汇集6| 精品一区二区三区av网在线观看| 少妇 在线观看| 色在线成人网| 露出奶头的视频| 大型av网站在线播放| 少妇裸体淫交视频免费看高清 | 首页视频小说图片口味搜索| 亚洲国产欧美网| 国产无遮挡羞羞视频在线观看| 亚洲国产看品久久| 在线看a的网站| 精品一区二区三区四区五区乱码| 窝窝影院91人妻| 曰老女人黄片| 日韩 欧美 亚洲 中文字幕| 久久久久久大精品| 操出白浆在线播放| 日日爽夜夜爽网站| 老汉色∧v一级毛片| 国产主播在线观看一区二区| 人成视频在线观看免费观看| 91麻豆av在线| 欧美国产精品va在线观看不卡| 99精国产麻豆久久婷婷| 国产成人免费无遮挡视频| 国产伦一二天堂av在线观看| 91字幕亚洲| 亚洲成av片中文字幕在线观看| 一级片免费观看大全| 国产精品亚洲av一区麻豆| 日日爽夜夜爽网站| 伊人久久大香线蕉亚洲五| 99热只有精品国产| 久久久久国内视频| 中出人妻视频一区二区| 久久精品亚洲熟妇少妇任你| 国产真人三级小视频在线观看| 亚洲av电影在线进入| 男人操女人黄网站| 日日爽夜夜爽网站| 嫩草影视91久久| 欧美黑人欧美精品刺激| 久久精品亚洲av国产电影网| 人人妻人人添人人爽欧美一区卜| 国产99白浆流出| 国产一卡二卡三卡精品| 国产精品偷伦视频观看了| 亚洲国产中文字幕在线视频| 一级a爱视频在线免费观看| 成人永久免费在线观看视频| 男女之事视频高清在线观看| 香蕉丝袜av| 一进一出抽搐动态| 9色porny在线观看| 999久久久精品免费观看国产| 少妇被粗大的猛进出69影院| 一边摸一边抽搐一进一小说| 天堂影院成人在线观看| 国产精品亚洲一级av第二区| 午夜福利,免费看| 久久久久久久久久久久大奶| 正在播放国产对白刺激| 国产一区二区三区视频了| 日韩欧美一区二区三区在线观看| 免费看十八禁软件| xxx96com| 12—13女人毛片做爰片一| 亚洲视频免费观看视频| 亚洲色图 男人天堂 中文字幕| 伦理电影免费视频| 动漫黄色视频在线观看| 日本精品一区二区三区蜜桃| 亚洲专区中文字幕在线| 99re在线观看精品视频| а√天堂www在线а√下载| 日韩一卡2卡3卡4卡2021年| 亚洲精品av麻豆狂野| 亚洲av五月六月丁香网| 最好的美女福利视频网| 亚洲七黄色美女视频| 最新在线观看一区二区三区| 欧美 亚洲 国产 日韩一| 国产男靠女视频免费网站| 国产一区二区在线av高清观看| 欧美乱妇无乱码| 欧美久久黑人一区二区| 悠悠久久av| 久久精品影院6| 欧美日韩视频精品一区| 成人精品一区二区免费| 在线视频色国产色| 一个人观看的视频www高清免费观看 | 成人三级黄色视频| 交换朋友夫妻互换小说| 国产人伦9x9x在线观看| 国产精品综合久久久久久久免费 | 中文字幕人妻丝袜制服| 熟女少妇亚洲综合色aaa.| 欧美 亚洲 国产 日韩一| 国产精品久久久av美女十八| 母亲3免费完整高清在线观看| 老司机福利观看| 99国产精品一区二区三区| 亚洲欧美日韩无卡精品| 母亲3免费完整高清在线观看| 欧洲精品卡2卡3卡4卡5卡区| 不卡一级毛片| 日韩视频一区二区在线观看| 一级片'在线观看视频| 欧美日韩av久久| 欧美黑人精品巨大| tocl精华| 91在线观看av| avwww免费| 免费搜索国产男女视频| 国产成人av激情在线播放| 久久婷婷成人综合色麻豆| 国产精品二区激情视频| 18禁黄网站禁片午夜丰满| 久久热在线av| 欧美日韩精品网址| 亚洲精品美女久久av网站| 久久香蕉激情| 成人永久免费在线观看视频| 丰满的人妻完整版| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品亚洲精品国产色婷小说| 日本精品一区二区三区蜜桃| 99久久综合精品五月天人人| 丰满饥渴人妻一区二区三| 亚洲精品粉嫩美女一区| 老司机午夜福利在线观看视频| 少妇 在线观看| 久久久水蜜桃国产精品网| e午夜精品久久久久久久| 欧美日韩精品网址| 一个人免费在线观看的高清视频| 中文字幕最新亚洲高清| 欧美日本亚洲视频在线播放| 久久性视频一级片| 欧美性长视频在线观看| 亚洲精品美女久久av网站| 国产亚洲欧美98| 曰老女人黄片| 日韩三级视频一区二区三区| 久久精品国产99精品国产亚洲性色 | 国产成人免费无遮挡视频| 无遮挡黄片免费观看| 又大又爽又粗| 成人精品一区二区免费| 亚洲激情在线av| 中文字幕人妻熟女乱码| 黄色成人免费大全| 国产亚洲欧美在线一区二区| 校园春色视频在线观看| 亚洲全国av大片| 18禁美女被吸乳视频| 成人影院久久| 亚洲国产精品999在线| 在线观看午夜福利视频| 国产亚洲精品一区二区www| 精品人妻在线不人妻| 国产熟女午夜一区二区三区| 亚洲色图av天堂| 长腿黑丝高跟| 多毛熟女@视频| 国产成人欧美| 日本a在线网址| 精品久久久久久久久久免费视频 | 欧美日韩一级在线毛片| 欧美精品啪啪一区二区三区| 国产免费现黄频在线看| 国产高清videossex| 午夜精品国产一区二区电影| 色精品久久人妻99蜜桃| 久久精品亚洲精品国产色婷小说| 91成人精品电影| 叶爱在线成人免费视频播放| 国产国语露脸激情在线看| 精品卡一卡二卡四卡免费| 国产野战对白在线观看| 国产成人精品久久二区二区91| 久久这里只有精品19| 精品国内亚洲2022精品成人| 久久久久久免费高清国产稀缺| 啦啦啦 在线观看视频| 在线观看免费午夜福利视频| 一级,二级,三级黄色视频| 国产成人av教育| xxx96com| 侵犯人妻中文字幕一二三四区| 巨乳人妻的诱惑在线观看| 欧美激情高清一区二区三区| 午夜精品久久久久久毛片777| xxx96com| 麻豆成人av在线观看| 两个人免费观看高清视频| 久久99一区二区三区| 色婷婷av一区二区三区视频| 久久国产亚洲av麻豆专区| 亚洲 国产 在线| 久久久久久久午夜电影 | 99国产精品99久久久久| 电影成人av| 亚洲第一欧美日韩一区二区三区| 一级毛片精品| 国产高清视频在线播放一区| 久久天躁狠狠躁夜夜2o2o| 真人做人爱边吃奶动态| 午夜a级毛片| 欧美日韩福利视频一区二区| 热re99久久精品国产66热6| 精品国产乱码久久久久久男人| 久久人人爽av亚洲精品天堂| 80岁老熟妇乱子伦牲交| 最新在线观看一区二区三区| 不卡一级毛片| 99riav亚洲国产免费| 一本大道久久a久久精品| 91精品三级在线观看| 欧美另类亚洲清纯唯美| 日韩欧美国产一区二区入口| 久久香蕉激情| 动漫黄色视频在线观看| 免费看十八禁软件| 欧美人与性动交α欧美软件| 亚洲人成77777在线视频| 9色porny在线观看| 在线观看一区二区三区| 欧美丝袜亚洲另类 | 新久久久久国产一级毛片| 国产精品偷伦视频观看了| 757午夜福利合集在线观看| 女性生殖器流出的白浆| 夜夜躁狠狠躁天天躁| 天天影视国产精品| 操出白浆在线播放| 999精品在线视频| 一级片免费观看大全| 91在线观看av| 90打野战视频偷拍视频| 久久青草综合色| 国产主播在线观看一区二区| 91麻豆av在线| 丁香欧美五月| 国产精品亚洲av一区麻豆| 欧美日韩亚洲高清精品| 十八禁人妻一区二区| 高清欧美精品videossex| 十八禁网站免费在线| 欧美日韩亚洲国产一区二区在线观看| 神马国产精品三级电影在线观看 | 伊人久久大香线蕉亚洲五| 午夜免费激情av| 色老头精品视频在线观看| 视频区欧美日本亚洲| 制服人妻中文乱码| 另类亚洲欧美激情| 夫妻午夜视频| 国产一区在线观看成人免费| 最新在线观看一区二区三区| 波多野结衣高清无吗| 制服诱惑二区| 亚洲熟女毛片儿| 日日干狠狠操夜夜爽| 精品高清国产在线一区| 国产三级黄色录像| 淫妇啪啪啪对白视频| 国产精品久久久av美女十八| e午夜精品久久久久久久| 久久久久久人人人人人| 午夜精品久久久久久毛片777| 少妇 在线观看| 国产成人精品久久二区二区91| 日日爽夜夜爽网站| 国产xxxxx性猛交| 夜夜爽天天搞| 色尼玛亚洲综合影院| 亚洲欧美精品综合一区二区三区| 欧美日本中文国产一区发布| 无人区码免费观看不卡| 免费在线观看视频国产中文字幕亚洲| 亚洲激情在线av| 三级毛片av免费| 一边摸一边抽搐一进一出视频| 岛国视频午夜一区免费看| 亚洲色图 男人天堂 中文字幕| 我的亚洲天堂| 日本黄色视频三级网站网址| 亚洲熟妇中文字幕五十中出 | 国产一区二区三区视频了| 亚洲av电影在线进入| 欧美成人免费av一区二区三区| 亚洲三区欧美一区| 国产三级黄色录像| 狂野欧美激情性xxxx| 一边摸一边做爽爽视频免费| 91精品三级在线观看| 亚洲成a人片在线一区二区| 欧美不卡视频在线免费观看 | 大型av网站在线播放| 一本综合久久免费| 国产亚洲精品久久久久5区| 免费一级毛片在线播放高清视频 | 国产欧美日韩一区二区精品| 黄色 视频免费看| 高清毛片免费观看视频网站 | 国产有黄有色有爽视频| 亚洲一区二区三区不卡视频| 极品教师在线免费播放| 国产伦一二天堂av在线观看| 男女之事视频高清在线观看| www日本在线高清视频| 极品教师在线免费播放| 好看av亚洲va欧美ⅴa在| 久久久久久大精品| 亚洲欧美一区二区三区久久| 国产精品九九99| 久久青草综合色| 91精品三级在线观看| 久久国产精品人妻蜜桃| 天堂动漫精品| 国产精品电影一区二区三区| 国产精品九九99| 啦啦啦在线免费观看视频4| 免费高清视频大片| 脱女人内裤的视频| 久久香蕉国产精品| 丁香六月欧美| 亚洲人成77777在线视频| 亚洲中文日韩欧美视频| 久久久精品欧美日韩精品| 成人精品一区二区免费| 久久久精品欧美日韩精品| 色播在线永久视频| 国产精品av久久久久免费| 色播在线永久视频| av在线天堂中文字幕 | a级片在线免费高清观看视频| 黑人猛操日本美女一级片| avwww免费| 亚洲avbb在线观看| 他把我摸到了高潮在线观看| 亚洲情色 制服丝袜| 少妇被粗大的猛进出69影院| 久久影院123| 日韩欧美三级三区| 夜夜躁狠狠躁天天躁| 十八禁网站免费在线| 精品国产乱子伦一区二区三区| 天堂动漫精品| 中文字幕av电影在线播放| 成年人免费黄色播放视频| 亚洲色图综合在线观看| 五月开心婷婷网| 男男h啪啪无遮挡| 午夜福利在线免费观看网站| 黑人操中国人逼视频| 黑丝袜美女国产一区| 国产熟女xx| 国产有黄有色有爽视频| 老司机深夜福利视频在线观看| 免费女性裸体啪啪无遮挡网站| 纯流量卡能插随身wifi吗| 久久午夜亚洲精品久久| 精品电影一区二区在线| 夜夜夜夜夜久久久久| 欧美av亚洲av综合av国产av| 国产单亲对白刺激| 99re在线观看精品视频| 两人在一起打扑克的视频| 国产一区二区三区视频了| 99精品久久久久人妻精品| 自线自在国产av| 日韩一卡2卡3卡4卡2021年| 国产三级在线视频| 亚洲欧美精品综合久久99| av在线播放免费不卡| 精品久久久久久成人av| 精品国产美女av久久久久小说| 黄色丝袜av网址大全| 亚洲,欧美精品.| 怎么达到女性高潮| 天天添夜夜摸| 亚洲午夜精品一区,二区,三区| 最近最新免费中文字幕在线| 国产一区二区三区视频了| 国产熟女午夜一区二区三区| 国产深夜福利视频在线观看| 99久久国产精品久久久| 老汉色∧v一级毛片| 一级片免费观看大全| 欧美乱码精品一区二区三区| 啦啦啦 在线观看视频| 久久人妻av系列| 在线观看一区二区三区激情| 亚洲一区高清亚洲精品| 一区福利在线观看| 亚洲情色 制服丝袜| 国产午夜精品久久久久久| 看黄色毛片网站| 叶爱在线成人免费视频播放| 9热在线视频观看99| 99riav亚洲国产免费| 老司机靠b影院| 99热国产这里只有精品6| 久久欧美精品欧美久久欧美| 大码成人一级视频| a级片在线免费高清观看视频| 脱女人内裤的视频| 欧美日韩亚洲综合一区二区三区_| 超碰成人久久| 亚洲一区二区三区欧美精品| 悠悠久久av| 日本vs欧美在线观看视频| 国产在线观看jvid| 中文字幕最新亚洲高清| 久久久久久久午夜电影 | 69av精品久久久久久| 美国免费a级毛片| 90打野战视频偷拍视频| 淫妇啪啪啪对白视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品在线观看二区| 亚洲色图 男人天堂 中文字幕| 夜夜看夜夜爽夜夜摸 | 欧美一级毛片孕妇| 国产亚洲精品第一综合不卡| 欧美日韩福利视频一区二区| 精品国产亚洲在线| 久久天躁狠狠躁夜夜2o2o| 欧美国产精品va在线观看不卡| 黑丝袜美女国产一区| 国产精品野战在线观看 | 在线播放国产精品三级| 在线观看一区二区三区激情| 两人在一起打扑克的视频| 97碰自拍视频| 久久香蕉激情| 丝袜在线中文字幕| 久久精品亚洲av国产电影网| 国产亚洲av高清不卡| 亚洲精华国产精华精| 手机成人av网站| 性色av乱码一区二区三区2| 美国免费a级毛片| 黑丝袜美女国产一区| 国产伦一二天堂av在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 人人妻人人爽人人添夜夜欢视频| 两个人看的免费小视频| 18禁国产床啪视频网站| 天堂影院成人在线观看| 国产av一区二区精品久久| 最近最新中文字幕大全电影3 | 久久久久亚洲av毛片大全| 色综合欧美亚洲国产小说| 午夜福利,免费看| 美女福利国产在线| 色综合站精品国产| 亚洲中文日韩欧美视频| 国产一区二区激情短视频| 免费在线观看日本一区| 日本精品一区二区三区蜜桃| 国产成人av激情在线播放| 岛国视频午夜一区免费看| 精品国产乱子伦一区二区三区| 88av欧美| 久久国产精品人妻蜜桃| videosex国产| 又黄又爽又免费观看的视频| 亚洲色图 男人天堂 中文字幕| 欧美日本亚洲视频在线播放| 成人三级做爰电影| 日本黄色日本黄色录像| 又黄又粗又硬又大视频| 国产精品偷伦视频观看了| 精品人妻1区二区| 成人18禁在线播放| 亚洲av日韩精品久久久久久密| 欧美不卡视频在线免费观看 | 天天躁狠狠躁夜夜躁狠狠躁| 国产精品一区二区免费欧美| 欧美日韩亚洲综合一区二区三区_| 在线观看免费视频日本深夜| aaaaa片日本免费| 精品人妻在线不人妻| 久久国产亚洲av麻豆专区| 亚洲一码二码三码区别大吗| 欧美日韩亚洲高清精品| 老司机福利观看| 国产欧美日韩精品亚洲av| 不卡av一区二区三区| 超色免费av| 新久久久久国产一级毛片| av天堂久久9| 国产一区二区三区在线臀色熟女 | 在线观看日韩欧美| 国产又色又爽无遮挡免费看| 一级毛片精品| 国产精品 国内视频| 日韩免费av在线播放| 精品久久久久久成人av| 欧美激情极品国产一区二区三区| 手机成人av网站| 99热国产这里只有精品6| 午夜a级毛片| 免费一级毛片在线播放高清视频 | 国产成人啪精品午夜网站| 久久香蕉国产精品| 老熟妇乱子伦视频在线观看| 久久香蕉激情| 在线观看免费视频网站a站| 日日摸夜夜添夜夜添小说| 亚洲黑人精品在线| 欧美+亚洲+日韩+国产| 色尼玛亚洲综合影院| 国产免费现黄频在线看| 精品少妇一区二区三区视频日本电影| 99久久精品国产亚洲精品| 宅男免费午夜| 国产免费男女视频| 桃红色精品国产亚洲av| 亚洲国产精品sss在线观看 | 亚洲激情在线av| 好看av亚洲va欧美ⅴa在| 丝袜人妻中文字幕| 女人高潮潮喷娇喘18禁视频| av天堂久久9| 欧美激情高清一区二区三区| 亚洲精品中文字幕在线视频| 欧美日韩亚洲国产一区二区在线观看| 国产午夜精品久久久久久| 水蜜桃什么品种好| 免费不卡黄色视频| 欧美老熟妇乱子伦牲交| 国产激情欧美一区二区| 国产午夜精品久久久久久| 又大又爽又粗| 精品国产乱子伦一区二区三区| 国产精品久久久人人做人人爽| 脱女人内裤的视频| 久久久国产精品麻豆| 香蕉久久夜色| 免费观看人在逋| 日韩精品中文字幕看吧| a在线观看视频网站| 韩国精品一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 欧美激情 高清一区二区三区| 久久 成人 亚洲| 男女之事视频高清在线观看| 人人妻人人澡人人看| 国产精品影院久久| √禁漫天堂资源中文www| 成年人免费黄色播放视频| 国产精品久久久久成人av| a在线观看视频网站| 亚洲 欧美 日韩 在线 免费| 99香蕉大伊视频| 国产精品久久久av美女十八| 亚洲第一欧美日韩一区二区三区| 国产真人三级小视频在线观看| 正在播放国产对白刺激| 90打野战视频偷拍视频| 精品福利观看| 中文字幕人妻熟女乱码| 亚洲第一av免费看| 久久久久久久久免费视频了| 欧美日韩亚洲综合一区二区三区_| 欧美日韩精品网址| 国产精品成人在线| 国产三级在线视频| 高清在线国产一区| 啦啦啦 在线观看视频| 欧美国产精品va在线观看不卡| 国产精品偷伦视频观看了| 国产精品影院久久| 欧美另类亚洲清纯唯美| 中文字幕人妻丝袜一区二区| 久久精品国产综合久久久| 久久久久久久久久久久大奶| 色哟哟哟哟哟哟|