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

    基于TET瞬態(tài)特征提取的滾動(dòng)軸承早期故障診斷研究*

    2021-06-22 08:16:50陳志剛車(chē)昊陽(yáng)
    機(jī)電工程 2021年6期
    關(guān)鍵詞:內(nèi)圈外圈瞬態(tài)

    陳志剛,趙 杰,張 楠,車(chē)昊陽(yáng)

    (1.北京建筑大學(xué) 機(jī)電與車(chē)輛工程學(xué)院,北京 100044;2.北京市建筑安全監(jiān)測(cè)工程技術(shù)研究中心, 北京 100044;3.中國(guó)石油集團(tuán)川慶鉆探工程有限公司 長(zhǎng)慶井下技術(shù)作業(yè)公司,陜西 西安 710021)

    0 引 言

    旋轉(zhuǎn)機(jī)械設(shè)備運(yùn)行環(huán)境復(fù)雜多變,滾動(dòng)軸承作為其重要部件,在其中發(fā)揮著不可或缺的作用。旋轉(zhuǎn)機(jī)械設(shè)備一旦發(fā)生事故,將會(huì)造成巨大經(jīng)濟(jì)損失和員工傷亡,所以對(duì)軸承進(jìn)行必要的故障診斷意義重大[1]。

    旋轉(zhuǎn)機(jī)械設(shè)備運(yùn)轉(zhuǎn)過(guò)程中通常夾雜著來(lái)自不同振源的噪聲,使得瞬態(tài)故障振動(dòng)信號(hào)難以被提取、識(shí)別,不能被有效地診斷和分析。

    在實(shí)驗(yàn)采集到的振動(dòng)信號(hào)中,通常夾雜著短時(shí)間的瞬變特征,而早期微弱故障往往發(fā)生在信號(hào)瞬變的時(shí)刻,在強(qiáng)噪聲背景下的特征提取困難,不易識(shí)別,所以研究高效的時(shí)頻分析方法尤為重要[2]。

    時(shí)頻分析(TFA)方法對(duì)于時(shí)變信號(hào)十分有效,在過(guò)去的幾十年里受到了業(yè)界極大的關(guān)注,其自身也得到了長(zhǎng)足的發(fā)展。但是,傳統(tǒng)的TFA技術(shù)依然存在較大的問(wèn)題,如信號(hào)兩邊的交叉項(xiàng)干擾、分解模態(tài)混疊、海森伯格不確定性原理等問(wèn)題。

    短時(shí)傅立葉變換[3](short-time Fourier transform, STFT)的時(shí)-頻分辨率較好,但對(duì)于不同的信號(hào)適應(yīng)性較差,需重新選擇窗函數(shù);小波變換[4](wavelet transform,WT)的窗口大小可根據(jù)信號(hào)頻率的不同做出改變,克服了STFT的缺點(diǎn),但是小波基的選擇受人為影響較大;Wigner-Ville[5]分布具有較高的時(shí)頻分辨率,對(duì)于多分量信號(hào)存在交叉干擾項(xiàng),分析效果較差。

    經(jīng)驗(yàn)?zāi)B(tài)分解(empirical mode decomposition, EMD)是HUANG N E[6]在1998年提出的一種自適應(yīng)時(shí)頻信號(hào)處理方法。該方法對(duì)于非線性、非平穩(wěn)的信號(hào)的處理效果較好,解決了WT人為選擇小波基函數(shù)的問(wèn)題;但是該方法也存在端點(diǎn)延拓[7]、模態(tài)混疊的問(wèn)題[8]。

    為了改善EMD的缺點(diǎn),聚合經(jīng)驗(yàn)?zāi)B(tài)分解[9](ensemble empirical mode decomposition, EEMD)被提了出來(lái)。該方法在分解過(guò)程中加入隨機(jī)白噪聲,這使得模態(tài)的混疊問(wèn)題在一定程度上得到了一些改善,但端點(diǎn)效應(yīng)問(wèn)題仍然存在。

    在總結(jié)了以上問(wèn)題的基礎(chǔ)上,TORRES M E[10]在EMD、EEMD的基礎(chǔ)上,提出了一種自適應(yīng)噪聲完備集合模態(tài)分解(CEEMDAN)。該方法具有更好的收斂性,且重構(gòu)誤差基本為0,重構(gòu)后具有良好的降噪效果。

    瞬態(tài)提取變換[11]是一種比較新的TFA方法。該方法的分辨率較高,能夠較好地提取出故障的瞬態(tài)特征。

    筆者進(jìn)行了大量文獻(xiàn)調(diào)研,發(fā)現(xiàn)目前基于TET的軸承瞬態(tài)故障診斷研究較少。因此,本文提出了一種瞬態(tài)特征提取的軸承故障診斷處理方法。首先,通過(guò)CEEMDAN對(duì)原始信號(hào)做自適應(yīng)降噪;然后,將重構(gòu)的信號(hào)利用瞬態(tài)特征提取算子(transient extraction operator,TEO)提取瞬態(tài)特征;最后,利用仿真信號(hào)和實(shí)驗(yàn)室軸承故障試驗(yàn)臺(tái)信號(hào)比較,對(duì)比STFT、同步擠壓變換(synchro squeezing transform,SST)、重分配法(reassignment method,RS)、解調(diào)TFA(demodulated TFA,DTFA)等方法。

    1 自適應(yīng)噪聲集合模態(tài)分解算法

    自適應(yīng)噪聲集合模態(tài)分解(CEEMDAN)是在EMD、EEMD的基礎(chǔ)上改進(jìn)得到的[12]。首先,該方法是將信號(hào)分解為一個(gè)模態(tài)分量IMF;然后,進(jìn)行總體平均計(jì)算,得到第一階IMF;最后,對(duì)剩余分量進(jìn)行上述操作。該方法有效地減小了模態(tài)混疊效應(yīng)和噪聲的遺留問(wèn)題。

    該算法具體步驟如下:

    (1)

    (2)

    將r1(t)作為輸入信號(hào),重復(fù)以上步驟,即可得到第K個(gè)模態(tài)分量IMFk(設(shè)總共有K個(gè)模態(tài)分量),最終剩余信號(hào)為R(t),則原始信號(hào)x(t)可得到K個(gè)IMF和一個(gè)R(t),即:

    (3)

    采用仿真信號(hào)可以驗(yàn)證CEEMDAN的降噪效果。由于軸承運(yùn)行過(guò)程中,原始信號(hào)中通常會(huì)具有調(diào)頻調(diào)幅成分,此處設(shè)仿真信號(hào)f(t)為:

    (4)

    式中:f1(t)—余弦信號(hào);f2(t),f3(t)—調(diào)頻信號(hào);f4(t)—隨機(jī)白噪聲,采樣頻率設(shè)置1 kHz。

    選用其中1 s的數(shù)據(jù),通過(guò)CEEMDAN進(jìn)行模態(tài)分解,可以得到9個(gè)模態(tài)分量IMF,如圖1所示。

    圖1 CEEMDAN分解結(jié)果

    然后,分別計(jì)算9個(gè)IMF的峭度值[13,14],如圖2所示。

    圖2 IMF峭度

    由圖2可以看出:第2、3、4個(gè)模態(tài)的峭度值較大,相關(guān)度較高;因此,此處選其進(jìn)行重構(gòu)。

    信號(hào)的預(yù)處理如圖3所示。

    圖3 信號(hào)預(yù)處理

    圖3中,重構(gòu)降噪之后的信號(hào)如圖3(b)所示。跟圖3(a)中的原始加噪信號(hào)相比,重構(gòu)降噪之后的信號(hào)濾除了部分噪聲,效果較好。

    2 瞬態(tài)提取算法

    由于該方法是基于短時(shí)傅里葉變化(STFT),不需擴(kuò)展參數(shù)和先驗(yàn)信息[11]。STFT表達(dá)式為:

    (5)

    式中:g(u-t)—可移動(dòng)窗口;s(u)—輸入信號(hào)。

    t0=0.5 s時(shí),δ(t)的時(shí)頻譜如圖4所示。

    圖4 t0=0.5 s時(shí)的時(shí)頻譜

    在數(shù)學(xué)上,δ(t)函數(shù)是一個(gè)沖擊函數(shù),在零點(diǎn)處的值無(wú)限大,非零處的值為零,且積分為1。此處令A(yù)=1,t0=0.5,則時(shí)域和頻域圖像如圖4(a,b)所示。由此,將其看作具有瞬態(tài)特征信號(hào)的理想模型。

    δ(t)通??杀硎緸?

    sδ(t)=A·δ(t-t0)

    (6)

    由于STFT的局限性,又根據(jù)海森堡(Heisenberg)測(cè)不準(zhǔn)原理,即使是δ(t)函數(shù)也不可能實(shí)現(xiàn)理想描述,為了探究STFT對(duì)δ(t)的時(shí)頻能量分布,筆者將式(6)代入式(5),可得:

    A·g(t0-t)·e-iωt0

    (7)

    因?yàn)閨e-iωt0|=1,δ(t)的STFT能量分布可以表示為:

    |G(t,ω)|=A·g(t0-t)

    (8)

    由式(8)可以看出:窗口函數(shù)g(t)在時(shí)域是緊湊的,頻域能量主要分布在t0時(shí)刻,并達(dá)到最大值A(chǔ)·g(0)。

    一個(gè)δ(t)函數(shù)的STFT由一系列具有相同群延遲(group delay, GD)的δ(t)函數(shù)構(gòu)成,通過(guò)計(jì)算G(t,ω)相對(duì)于頻率變量的導(dǎo)數(shù),可以得到精確GDs,即:

    ?ωG(t,ω)=?ω(A·g(t0-t)·e-iωt0)=
    -it0·A·g(t0-t)·e-iωt0=

    -it0·G(t,ω)

    (9)

    對(duì)于任意的(t,ω),若G(t,ω)≠0,那么二維GDt0(t,ω)可以定義為:

    (10)

    群延遲GD和提取算子TEO如圖5所示。

    圖5 群延遲GD和提取算子TEO

    為了更清楚地說(shuō)明GD,圖5(b)為ω0=50 Hz時(shí)GD的頻率片段。

    在t∈[t-Δ,t+Δ](Δ為時(shí)間偏移)時(shí)刻,所有二維GD的值都與t0=0.5 s時(shí)刻的值相同。對(duì)于理想的TFA方法,信號(hào)的能量應(yīng)該只出現(xiàn)在t0時(shí)刻,而不是在一定的范圍內(nèi)發(fā)生擴(kuò)散。

    為了消除這一能量擴(kuò)散帶來(lái)的影響,筆者只保留t0時(shí)刻的能量,提出了一種瞬態(tài)提取算子(TEO)的TFA方法:

    TEO(t,ω)=δ(t-t0(t,ω))

    (11)

    由此可得:

    (12)

    則:

    δ(t-t0(t,ω))=δ(t-t0)

    (13)

    由式(13)可以看出:TEO的值只有在t=t0時(shí)才可以提取時(shí)頻系數(shù)G(t,ω),如圖5(c)所示。

    因?yàn)槭?11)具有瞬態(tài)提取行為,筆者采用瞬態(tài)提取算子(TEO)的變換稱(chēng)為瞬態(tài)提取變換(TET),即:

    Te(t,ω)=G(t,ω)·TEO(t,ω)

    (14)

    實(shí)際診斷中,僅對(duì)信號(hào)做頻域處理往往是不夠的,有時(shí)還需要得到信號(hào)的時(shí)域信息。

    由STFT信號(hào)重構(gòu)表達(dá)式:

    (15)

    很容易可以得出TET重構(gòu)表達(dá)式:

    (16)

    3 仿真分析

    為驗(yàn)證筆者所提方法的瞬態(tài)特征提取效果,此處采用瞬變仿真信號(hào)進(jìn)行分析驗(yàn)證。

    仿真信號(hào)如圖6所示。

    圖6 仿真信號(hào)

    為了衡量不同方法的處理效果,筆者采用Renyi熵來(lái)估計(jì)信號(hào)的分散程度[15]。Renyi熵也是信息熵的一種,可以表示為:

    (17)

    在圖6的仿真信號(hào)中,加入信噪比1 dB~30 dB的高斯白噪聲,采用短時(shí)傅里葉變換(STFT)、瞬態(tài)提取變換(TET)、小波變換(WT)、同步擠壓變換(SST)、重分配變換(RS)、離散函數(shù)分析(DTFA)來(lái)對(duì)其加噪信號(hào)進(jìn)行分析處理。其中,TF結(jié)果的Renyi熵越低,表示該方法的效果越好,能產(chǎn)生更加集中的TF表示。

    在信噪比1 dB~30 dB時(shí),6種方法TF結(jié)果的Renyi熵水平如圖7所示。

    圖7 6種方法在信噪比1 dB~30 dB時(shí)TF結(jié)果的Renyi熵水平

    由圖7可以看出:在不同噪聲水平下,TET的效果最好。

    筆者選取Renyi熵較低的4種方法,即STFT、SST、RS、TET,查看其TF表示。4種方法的TF表示及其局部放大結(jié)果如圖8所示。

    圖8 4種方法TF表示及局部放大

    由圖8可知:STFT、SST、RS的TF表示結(jié)果比較分散;而TET結(jié)果能量集中,沒(méi)有發(fā)散現(xiàn)象,效果較好。

    4 實(shí)驗(yàn)驗(yàn)證

    4.1 信號(hào)采集

    筆者采用實(shí)驗(yàn)室軸承試驗(yàn)臺(tái)進(jìn)行數(shù)據(jù)采集。試驗(yàn)臺(tái)由交流電動(dòng)機(jī)、電機(jī)速度控制中心、支撐軸、測(cè)試軸承、傳感器等組成。

    傳感器選用美國(guó)PCB公司的352C33型ICP加速的傳感器,測(cè)試軸承為6105-SKF深溝球軸承,通過(guò)電火花加工的方式在軸承外圈和內(nèi)圈刻蝕直徑0.144 mm、深度0.232 mm的微小裂痕,模擬軸承運(yùn)行過(guò)程中外圈和內(nèi)圈產(chǎn)生的早期微弱故障。

    軸承實(shí)驗(yàn)臺(tái)及其部件如圖9所示。

    圖9 實(shí)驗(yàn)臺(tái)及其部件

    圖9中,采樣頻率Fs=12 kHz,轉(zhuǎn)速1 750 r/min。根據(jù)式(18,19),可計(jì)算出外圈和內(nèi)圈故障特征頻率為104.5 Hz和157.9 Hz。

    軸承外圈特征頻率為:

    (18)

    內(nèi)圈特征頻率為:

    (19)

    式中:r—轉(zhuǎn)速;n—滾珠個(gè)數(shù);d—滾動(dòng)體直徑;D—軸承節(jié)徑;α—滾動(dòng)體接觸角。

    4.2 軸承外圈信號(hào)分析

    軸承外圈信號(hào)的預(yù)處理如圖10所示。

    圖10 外圈信號(hào)預(yù)處理

    圖10中,筆者首先對(duì)外圈信號(hào)進(jìn)行處理,以展現(xiàn)軸承外圈的原始振動(dòng)信號(hào);傳感器在3點(diǎn)鐘方向采集??梢钥吹?瞬變信號(hào)完全被噪聲淹沒(méi),無(wú)法清晰地觀察;

    然后對(duì)信號(hào)進(jìn)行CEEMDAN分解,并計(jì)算各個(gè)模態(tài)分量IMF的峭度值,共有11個(gè)IMF分量,其中,第1、2、4個(gè)IMF峭度值較大,與原信號(hào)相關(guān)度較高,選其進(jìn)行重構(gòu)??梢钥闯?經(jīng)過(guò)重構(gòu)之后,濾除了大部分噪聲,能夠清晰地看到故障的瞬變特征。

    筆者選用第3節(jié)中效果較好的4種方法(STFT、SST、RS、TET)對(duì)降噪信號(hào)進(jìn)行處理,并生成TF表示,如圖11所示。

    圖11 4種方法TF表示及局部放大

    由圖11可以看出:STFT結(jié)果比較發(fā)散,無(wú)法準(zhǔn)確定位;RS和SST只提供了較為粗略的TF表示,依舊不能定位TF信息;TET則表現(xiàn)出了比較清晰的TF結(jié)果,效果較好。

    瞬態(tài)特征分量如圖12所示。

    圖12 瞬態(tài)特征分量

    圖12中,相比于降噪后的信號(hào)(圖10(c)),提取TET的結(jié)果(圖12(a))的瞬變特征有了明顯改善;

    圖12(a)的相鄰瞬變信號(hào)的間隔為9.5 ms,對(duì)其做Hilbert包絡(luò)譜分析后,可以清晰地觀察到故障頻率f0及其2倍、3倍、4倍頻率,可見(jiàn)其效果較好。

    4.3 軸承內(nèi)圈信號(hào)分析

    由于軸承內(nèi)圈、外圈的故障機(jī)理具有一定的差異性,來(lái)自外圈和滾動(dòng)體的振動(dòng)噪聲與內(nèi)圈本身的噪聲相疊加,會(huì)導(dǎo)致其故障特征更加難以識(shí)別[16]。為了驗(yàn)證本文所提方法的適用性,筆者對(duì)內(nèi)圈故障信號(hào)進(jìn)行分析。

    軸承內(nèi)圈的信號(hào)預(yù)處理如圖13所示。

    圖13 內(nèi)圈信號(hào)預(yù)處理

    圖13中,筆者首先采用與處理外圈信號(hào)相同的方法,對(duì)內(nèi)圈故障信號(hào)進(jìn)行了自適應(yīng)降噪處理,用CEEMDAN方法對(duì)其進(jìn)行了分解,通過(guò)調(diào)整CEEMDAN信噪比及迭代次數(shù),得到了11個(gè)IMF分量,并計(jì)算了每個(gè)分量的峭度值。

    由圖13(b)可知:第1、2、11分量峭度值較大,說(shuō)明與原信號(hào)相關(guān)度較高;剩余分量則為噪聲分量,通過(guò)重構(gòu)1、2、11分量達(dá)到降噪的目的,如圖13(c)所示。與原信號(hào)相比,其效果較好,濾除了大部分的噪聲。

    然后,筆者對(duì)降噪后的信號(hào)采用4種方法(STFT、SST、RS、TET)進(jìn)行了對(duì)比分析,并計(jì)算了每種方法的Renyi熵,如表1所示。

    表1 STFT、SST、RS、TET的Renyi熵

    由表1可以看出:RS方法的Renyi熵?cái)?shù)值最大,TET的數(shù)值最小,這說(shuō)明了TET的效果最好。

    4種方法的TF表示如圖14所示。

    圖14 4種方法的TF表示

    圖14中的結(jié)果表明:STFT的TF結(jié)果能量發(fā)散最為嚴(yán)重,無(wú)法準(zhǔn)確定位;SST結(jié)果能量受Heisenberg不確定性原理的影響也比較發(fā)散;RS結(jié)果相對(duì)較好,但還是沒(méi)有TET結(jié)果更清晰、準(zhǔn)確。

    最后,筆者提取了TET結(jié)果的TF分量,其瞬態(tài)特征分量分析結(jié)果如圖15所示。

    圖15 瞬態(tài)特征分量

    由圖15可以看出:圖15(a)中的瞬變信號(hào)非常明顯,分量的間隔基本為6.3 ms,與故障頻率157.9 Hz相對(duì)應(yīng);對(duì)瞬態(tài)分量在Hilbert進(jìn)行包絡(luò)譜分析,可以在圖中清楚地觀察到內(nèi)圈信號(hào)的故障特征頻率f0及其2、3、4、5倍頻,可見(jiàn)其效果很好。

    5 結(jié)束語(yǔ)

    本文提出了一種瞬態(tài)特征提取的軸承早期故障診斷方法,能夠有效地對(duì)故障信號(hào)自適應(yīng)進(jìn)行濾波降噪,并提取出了軸承故障瞬態(tài)特征,結(jié)果表明其抗噪性較強(qiáng),效果也較為明顯。

    主要結(jié)論如下:

    (1)利用CEEMDAN對(duì)仿真信號(hào)及實(shí)驗(yàn)信號(hào)進(jìn)行了自適應(yīng)分解,并通過(guò)峭度選擇相關(guān)度較大的分量重構(gòu)降噪,實(shí)測(cè)數(shù)據(jù)表明其效果較好,噪聲的魯棒性較強(qiáng);

    (2)對(duì)降噪之后的信號(hào)進(jìn)行了瞬變信號(hào)特征提取,然后做Hilbert包絡(luò)譜分析,能夠容易地找到軸承故障特征頻率及其多倍頻,并且選取軸承的外圈和內(nèi)圈進(jìn)行了分析對(duì)比,結(jié)果表明其皆可以有效地分析出其故障頻率,這說(shuō)明了該方法的適應(yīng)性較強(qiáng);

    (3)將該方法與STFT、SST、RS、TET等方法作了比較分析,結(jié)果表明TET優(yōu)于其他先進(jìn)TF方法,且TET基于STFT,運(yùn)算量相當(dāng),可以將其應(yīng)用于實(shí)際中的軸承早期微弱故障診斷之中。

    另外,對(duì)于CEEMDAN分量的選擇問(wèn)題,還有待于在后續(xù)的研究中進(jìn)一步進(jìn)行優(yōu)化,使其與原信號(hào)的相關(guān)度更高。

    猜你喜歡
    內(nèi)圈外圈瞬態(tài)
    深溝球軸承外圈表面凹坑缺陷分析
    哈爾濱軸承(2022年1期)2022-05-23 13:13:16
    特種復(fù)合軸承內(nèi)圈推力滾道磨削用工裝設(shè)計(jì)
    哈爾濱軸承(2021年4期)2021-03-08 01:00:48
    角接觸球軸承外圈鎖口高度自動(dòng)檢測(cè)規(guī)改進(jìn)
    哈爾濱軸承(2020年2期)2020-11-06 09:22:34
    高壓感應(yīng)電動(dòng)機(jī)斷電重啟時(shí)的瞬態(tài)仿真
    主軸軸承內(nèi)圈鎖緊用臺(tái)階套的裝配
    十億像素瞬態(tài)成像系統(tǒng)實(shí)時(shí)圖像拼接
    內(nèi)圈帶缺陷中介軸承的動(dòng)力學(xué)建模與振動(dòng)響應(yīng)分析
    基于瞬態(tài)流場(chǎng)計(jì)算的滑動(dòng)軸承靜平衡位置求解
    DC/DC變換器中的瞬態(tài)特性分析
    軸承內(nèi)圈與軸的配合過(guò)盈量分析
    軸承(2011年10期)2011-07-25 01:36:22
    亚洲内射少妇av| 老司机福利观看| 日韩 亚洲 欧美在线| 欧美成人a在线观看| 麻豆国产97在线/欧美| 久久人人爽人人爽人人片va| 亚洲av.av天堂| 在线免费观看的www视频| 高清日韩中文字幕在线| www.色视频.com| 午夜a级毛片| 久久婷婷人人爽人人干人人爱| 国产私拍福利视频在线观看| 99久久中文字幕三级久久日本| 亚洲精品色激情综合| 亚洲国产日韩欧美精品在线观看| 给我免费播放毛片高清在线观看| 国产精品电影一区二区三区| 别揉我奶头~嗯~啊~动态视频| 一进一出抽搐动态| 久久草成人影院| 国产成人av教育| 欧美日韩乱码在线| 国产成人影院久久av| 精品一区二区三区视频在线| 狂野欧美白嫩少妇大欣赏| 成人性生交大片免费视频hd| 国产真实乱freesex| 亚洲人成伊人成综合网2020| 国产精品伦人一区二区| 人妻久久中文字幕网| 真实男女啪啪啪动态图| 日韩一区二区视频免费看| 亚洲午夜理论影院| 熟妇人妻久久中文字幕3abv| 老熟妇乱子伦视频在线观看| 中文亚洲av片在线观看爽| 午夜福利在线观看免费完整高清在 | 亚洲最大成人手机在线| 欧美一区二区亚洲| 色播亚洲综合网| 毛片女人毛片| 夜夜夜夜夜久久久久| 老师上课跳d突然被开到最大视频| 久久99热这里只有精品18| 欧美精品国产亚洲| 日韩强制内射视频| 很黄的视频免费| 精品久久久久久久末码| 欧美色视频一区免费| 在线免费观看不下载黄p国产 | 日韩在线高清观看一区二区三区 | 身体一侧抽搐| 国产黄片美女视频| 久9热在线精品视频| 精品久久久久久久久亚洲 | 午夜福利在线在线| 极品教师在线免费播放| 久久人人爽人人爽人人片va| 精品久久国产蜜桃| 日日撸夜夜添| 观看免费一级毛片| 蜜桃久久精品国产亚洲av| 国产aⅴ精品一区二区三区波| 精品国内亚洲2022精品成人| 亚洲性久久影院| 国产视频一区二区在线看| 大型黄色视频在线免费观看| 久久精品国产自在天天线| 九九久久精品国产亚洲av麻豆| 精品人妻视频免费看| 少妇人妻一区二区三区视频| 日韩一本色道免费dvd| 国产精品,欧美在线| 夜夜看夜夜爽夜夜摸| 国产三级中文精品| 欧美日韩精品成人综合77777| 五月玫瑰六月丁香| 免费看美女性在线毛片视频| 亚洲aⅴ乱码一区二区在线播放| 国产成人影院久久av| 欧美xxxx性猛交bbbb| 干丝袜人妻中文字幕| 国内久久婷婷六月综合欲色啪| 国产欧美日韩精品一区二区| 少妇高潮的动态图| 热99在线观看视频| 高清毛片免费观看视频网站| 免费在线观看日本一区| 中文亚洲av片在线观看爽| 成人毛片a级毛片在线播放| 亚洲美女黄片视频| 久久精品久久久久久噜噜老黄 | 国内揄拍国产精品人妻在线| 久久久久久国产a免费观看| 国产伦人伦偷精品视频| 99riav亚洲国产免费| 国产人妻一区二区三区在| 久久精品国产99精品国产亚洲性色| 欧洲精品卡2卡3卡4卡5卡区| 国产午夜精品论理片| 国产精品98久久久久久宅男小说| 欧美日韩精品成人综合77777| 日日撸夜夜添| 99国产精品一区二区蜜桃av| 两人在一起打扑克的视频| 日本在线视频免费播放| 亚洲va日本ⅴa欧美va伊人久久| 国产精品综合久久久久久久免费| 少妇裸体淫交视频免费看高清| 国产精品久久久久久久久免| 中文亚洲av片在线观看爽| 久久人人爽人人爽人人片va| 国产精品伦人一区二区| 亚洲电影在线观看av| 少妇的逼好多水| 欧美xxxx性猛交bbbb| 91狼人影院| 欧美日韩精品成人综合77777| 精品久久久久久久末码| 少妇的逼水好多| 2021天堂中文幕一二区在线观| 一区二区三区四区激情视频 | 日本黄色片子视频| 国产视频一区二区在线看| 亚洲熟妇中文字幕五十中出| 久久99热6这里只有精品| 少妇熟女aⅴ在线视频| 久久亚洲精品不卡| 国产麻豆成人av免费视频| 亚洲精华国产精华液的使用体验 | 久久欧美精品欧美久久欧美| 国模一区二区三区四区视频| 在线国产一区二区在线| 热99在线观看视频| 51国产日韩欧美| av在线亚洲专区| 成人国产麻豆网| 日本欧美国产在线视频| 亚洲无线观看免费| 免费电影在线观看免费观看| 成年免费大片在线观看| 免费观看在线日韩| 一区二区三区四区激情视频 | 亚洲男人的天堂狠狠| 在线天堂最新版资源| 国产黄a三级三级三级人| 97超级碰碰碰精品色视频在线观看| 欧美日韩瑟瑟在线播放| 久久久久久久久久黄片| 成人鲁丝片一二三区免费| 婷婷精品国产亚洲av| 亚洲av美国av| 欧美丝袜亚洲另类 | www.www免费av| 亚洲av成人精品一区久久| 成人国产一区最新在线观看| 亚洲五月天丁香| 午夜精品一区二区三区免费看| 波多野结衣高清无吗| 欧美精品国产亚洲| 1024手机看黄色片| 黄色欧美视频在线观看| 白带黄色成豆腐渣| 一区福利在线观看| 日本爱情动作片www.在线观看 | 国产私拍福利视频在线观看| av女优亚洲男人天堂| 亚洲av一区综合| 亚洲熟妇熟女久久| 国产男人的电影天堂91| 免费av毛片视频| 他把我摸到了高潮在线观看| 成人国产麻豆网| 国产精品av视频在线免费观看| 一a级毛片在线观看| 一区福利在线观看| 国产精品av视频在线免费观看| 久久精品夜夜夜夜夜久久蜜豆| 精华霜和精华液先用哪个| 欧美日韩精品成人综合77777| 国产不卡一卡二| 色噜噜av男人的天堂激情| 成人国产麻豆网| 国产高清不卡午夜福利| 国内精品久久久久久久电影| 国产欧美日韩精品一区二区| 悠悠久久av| 久久中文看片网| 在线免费观看的www视频| 国产黄a三级三级三级人| www.www免费av| 欧美高清性xxxxhd video| 免费一级毛片在线播放高清视频| 婷婷亚洲欧美| 色在线成人网| 色哟哟·www| 欧美色视频一区免费| 国产色婷婷99| 免费观看在线日韩| 免费不卡的大黄色大毛片视频在线观看 | 免费黄网站久久成人精品| 麻豆一二三区av精品| 中文字幕av成人在线电影| a级毛片免费高清观看在线播放| 国产亚洲av嫩草精品影院| 国产大屁股一区二区在线视频| 女同久久另类99精品国产91| 一个人免费在线观看电影| 国产高清三级在线| 日韩中字成人| x7x7x7水蜜桃| 身体一侧抽搐| 欧美成人一区二区免费高清观看| 欧美+日韩+精品| 在线免费观看的www视频| 韩国av在线不卡| 国产老妇女一区| 噜噜噜噜噜久久久久久91| 欧美国产日韩亚洲一区| 亚洲成人精品中文字幕电影| 久久久国产成人免费| 亚洲欧美日韩东京热| 亚洲自偷自拍三级| 女的被弄到高潮叫床怎么办 | 狂野欧美白嫩少妇大欣赏| 婷婷精品国产亚洲av| 亚洲无线在线观看| 成人美女网站在线观看视频| 草草在线视频免费看| 可以在线观看的亚洲视频| 免费无遮挡裸体视频| 成人特级av手机在线观看| 一卡2卡三卡四卡精品乱码亚洲| 亚洲一区二区三区色噜噜| 搡老妇女老女人老熟妇| 中文字幕高清在线视频| 88av欧美| 亚洲第一电影网av| 欧美日韩精品成人综合77777| 麻豆国产av国片精品| 欧美不卡视频在线免费观看| 九九热线精品视视频播放| 久久精品国产鲁丝片午夜精品 | 直男gayav资源| 22中文网久久字幕| 午夜精品在线福利| 小说图片视频综合网站| 亚洲七黄色美女视频| 99精品在免费线老司机午夜| 色吧在线观看| 老司机福利观看| 老熟妇仑乱视频hdxx| 一区二区三区免费毛片| 国产伦在线观看视频一区| 国产亚洲91精品色在线| 日韩欧美在线二视频| 神马国产精品三级电影在线观看| 在线免费十八禁| 欧美最黄视频在线播放免费| 国产精品一区二区三区四区免费观看 | 亚洲美女黄片视频| 日韩欧美 国产精品| 999久久久精品免费观看国产| av福利片在线观看| 国产精品综合久久久久久久免费| 国产午夜精品久久久久久一区二区三区 | 丝袜美腿在线中文| 国产欧美日韩一区二区精品| 中出人妻视频一区二区| 又爽又黄a免费视频| 在线国产一区二区在线| 在线a可以看的网站| 在线观看av片永久免费下载| 男女边吃奶边做爰视频| 成年女人毛片免费观看观看9| 欧美精品啪啪一区二区三区| 久久精品国产清高在天天线| 亚洲四区av| 看片在线看免费视频| 色综合亚洲欧美另类图片| 久久久久性生活片| 日韩精品青青久久久久久| 国产亚洲精品久久久com| 白带黄色成豆腐渣| 日韩av在线大香蕉| 婷婷亚洲欧美| 丰满乱子伦码专区| 特大巨黑吊av在线直播| 真实男女啪啪啪动态图| 国内精品宾馆在线| 老女人水多毛片| 精品一区二区三区视频在线观看免费| 婷婷六月久久综合丁香| 高清毛片免费观看视频网站| 成人欧美大片| 91av网一区二区| 欧美成人一区二区免费高清观看| 大又大粗又爽又黄少妇毛片口| 欧美日本亚洲视频在线播放| 99久久无色码亚洲精品果冻| 亚洲在线自拍视频| av中文乱码字幕在线| 国产精品,欧美在线| 51国产日韩欧美| 一个人看视频在线观看www免费| 丰满的人妻完整版| 国产淫片久久久久久久久| 直男gayav资源| 午夜免费男女啪啪视频观看 | 国产av在哪里看| 69人妻影院| 精品人妻一区二区三区麻豆 | 亚洲18禁久久av| 久久久国产成人免费| 啦啦啦韩国在线观看视频| 久久久午夜欧美精品| av在线天堂中文字幕| 中文字幕免费在线视频6| 欧美最黄视频在线播放免费| 999久久久精品免费观看国产| 可以在线观看毛片的网站| 网址你懂的国产日韩在线| 国产男靠女视频免费网站| 日本爱情动作片www.在线观看 | 18禁裸乳无遮挡免费网站照片| 国产色婷婷99| 波多野结衣高清无吗| 亚洲av第一区精品v没综合| 中文在线观看免费www的网站| 波野结衣二区三区在线| 亚洲精品粉嫩美女一区| 极品教师在线视频| 亚洲成人免费电影在线观看| 琪琪午夜伦伦电影理论片6080| 91麻豆av在线| 国产精品久久久久久久电影| 久久久精品欧美日韩精品| 又爽又黄a免费视频| 亚洲精品在线观看二区| 中国美女看黄片| 在线国产一区二区在线| 国产91精品成人一区二区三区| 色视频www国产| 动漫黄色视频在线观看| 欧美性感艳星| 乱码一卡2卡4卡精品| 99久久久亚洲精品蜜臀av| 久久九九热精品免费| 在线观看免费视频日本深夜| 亚洲五月天丁香| 国语自产精品视频在线第100页| 日韩av在线大香蕉| 三级男女做爰猛烈吃奶摸视频| 色播亚洲综合网| 老熟妇乱子伦视频在线观看| 亚洲国产欧洲综合997久久,| 中文在线观看免费www的网站| www.色视频.com| 精品久久久久久久久久久久久| 精品一区二区免费观看| 中文字幕人妻熟人妻熟丝袜美| 久久久久久伊人网av| 亚洲第一电影网av| 免费不卡的大黄色大毛片视频在线观看 | 成熟少妇高潮喷水视频| 搡女人真爽免费视频火全软件 | 精品国内亚洲2022精品成人| 十八禁网站免费在线| 又黄又爽又刺激的免费视频.| 国产精品一及| 人妻丰满熟妇av一区二区三区| 十八禁网站免费在线| 日韩大尺度精品在线看网址| 免费高清视频大片| 亚洲中文日韩欧美视频| 免费观看人在逋| 成人永久免费在线观看视频| 听说在线观看完整版免费高清| 日日摸夜夜添夜夜添av毛片 | 麻豆成人av在线观看| 欧美黑人巨大hd| av天堂中文字幕网| 俺也久久电影网| 国产午夜福利久久久久久| 精品欧美国产一区二区三| 悠悠久久av| 精品欧美国产一区二区三| 亚洲经典国产精华液单| 露出奶头的视频| 国内精品久久久久精免费| av国产免费在线观看| 精品人妻熟女av久视频| 成人特级黄色片久久久久久久| 欧美激情在线99| 国产免费男女视频| 国产精品女同一区二区软件 | 深爱激情五月婷婷| 国产精品电影一区二区三区| 亚洲精品色激情综合| 国内精品久久久久久久电影| 久久久久国产精品人妻aⅴ院| 精品久久久久久久久av| 日韩欧美国产在线观看| 国产主播在线观看一区二区| aaaaa片日本免费| 黄色一级大片看看| 精品一区二区免费观看| 美女高潮的动态| 成人国产综合亚洲| 床上黄色一级片| 国产色爽女视频免费观看| 亚洲美女视频黄频| 亚洲欧美精品综合久久99| 成人毛片a级毛片在线播放| 国产亚洲精品久久久com| 波野结衣二区三区在线| 99国产精品一区二区蜜桃av| 亚洲,欧美,日韩| 成年女人永久免费观看视频| 亚洲欧美日韩高清在线视频| 免费av毛片视频| 欧美中文日本在线观看视频| 久久精品人妻少妇| 简卡轻食公司| 国产免费一级a男人的天堂| 性插视频无遮挡在线免费观看| 国产免费男女视频| 国产亚洲91精品色在线| 欧美另类亚洲清纯唯美| 嫩草影视91久久| 久久午夜亚洲精品久久| 日日摸夜夜添夜夜添av毛片 | 男人舔女人下体高潮全视频| 在线看三级毛片| 在线a可以看的网站| 国产精品精品国产色婷婷| 亚洲国产精品久久男人天堂| 12—13女人毛片做爰片一| 久久午夜福利片| 搡老熟女国产l中国老女人| 久久人人爽人人爽人人片va| 观看免费一级毛片| 国产精品av视频在线免费观看| 午夜福利视频1000在线观看| 免费看美女性在线毛片视频| 国产成人a区在线观看| 丰满的人妻完整版| 在线看三级毛片| 国产精品久久久久久久电影| 男女之事视频高清在线观看| 亚洲午夜理论影院| 亚洲不卡免费看| 国产精品久久久久久精品电影| 国产淫片久久久久久久久| 国产视频一区二区在线看| 夜夜爽天天搞| 国产色婷婷99| 国产乱人视频| 很黄的视频免费| 午夜影院日韩av| 免费看av在线观看网站| 国产一区二区在线观看日韩| 亚洲人成网站在线播| 久久这里只有精品中国| 精品无人区乱码1区二区| 国产男靠女视频免费网站| 国产美女午夜福利| 亚洲 国产 在线| 1024手机看黄色片| 两个人的视频大全免费| 久久久久九九精品影院| 国产亚洲精品av在线| 我要看日韩黄色一级片| 亚洲午夜理论影院| 淫秽高清视频在线观看| 欧美最黄视频在线播放免费| 99久久中文字幕三级久久日本| 成人二区视频| 国产三级在线视频| 亚洲欧美清纯卡通| 免费av不卡在线播放| 99九九线精品视频在线观看视频| 九色国产91popny在线| 他把我摸到了高潮在线观看| 亚洲av电影不卡..在线观看| 国产真实伦视频高清在线观看 | 亚洲经典国产精华液单| 一区二区三区四区激情视频 | 色视频www国产| 自拍偷自拍亚洲精品老妇| 麻豆av噜噜一区二区三区| 三级国产精品欧美在线观看| 精品久久久久久久末码| 欧美最新免费一区二区三区| 俄罗斯特黄特色一大片| 欧美高清成人免费视频www| 亚洲熟妇中文字幕五十中出| 亚洲一级一片aⅴ在线观看| 极品教师在线视频| 91av网一区二区| 又爽又黄无遮挡网站| 亚洲 国产 在线| 精品无人区乱码1区二区| 欧美一区二区精品小视频在线| 精品久久久久久久末码| 日韩精品中文字幕看吧| 精品免费久久久久久久清纯| 亚洲午夜理论影院| а√天堂www在线а√下载| 伦精品一区二区三区| 在线观看美女被高潮喷水网站| 美女被艹到高潮喷水动态| 久久久久久久久大av| 中文在线观看免费www的网站| 国产精品av视频在线免费观看| 亚洲男人的天堂狠狠| 日韩精品有码人妻一区| 夜夜夜夜夜久久久久| 国产真实伦视频高清在线观看 | 精品久久久久久久久亚洲 | 久久精品夜夜夜夜夜久久蜜豆| 亚洲不卡免费看| 国产私拍福利视频在线观看| 欧美高清性xxxxhd video| 国产蜜桃级精品一区二区三区| 十八禁国产超污无遮挡网站| 五月玫瑰六月丁香| 乱码一卡2卡4卡精品| 91麻豆精品激情在线观看国产| 又黄又爽又刺激的免费视频.| 波多野结衣高清作品| 国产av一区在线观看免费| 亚洲av免费在线观看| 久久久久久久久久久丰满 | 在线免费十八禁| 12—13女人毛片做爰片一| 久久久色成人| 亚洲精品亚洲一区二区| 很黄的视频免费| 精品人妻1区二区| 亚洲人成伊人成综合网2020| 久久精品人妻少妇| 给我免费播放毛片高清在线观看| 国产精品一区二区性色av| 一边摸一边抽搐一进一小说| 在线观看66精品国产| 亚洲av免费高清在线观看| 欧美三级亚洲精品| aaaaa片日本免费| 3wmmmm亚洲av在线观看| 久久精品国产鲁丝片午夜精品 | 成人一区二区视频在线观看| 亚洲无线在线观看| 亚洲,欧美,日韩| 国产精品久久久久久av不卡| 国产精品亚洲一级av第二区| 麻豆成人午夜福利视频| 亚洲精品一区av在线观看| 校园春色视频在线观看| 免费av毛片视频| 99热精品在线国产| 听说在线观看完整版免费高清| 欧美激情久久久久久爽电影| 俺也久久电影网| 日韩人妻高清精品专区| 欧美中文日本在线观看视频| 久久婷婷人人爽人人干人人爱| 成人综合一区亚洲| 超碰av人人做人人爽久久| 国语自产精品视频在线第100页| 精品欧美国产一区二区三| 久久人人精品亚洲av| 黄色丝袜av网址大全| 亚洲国产欧美人成| 久久久久国产精品人妻aⅴ院| 午夜精品在线福利| 亚洲欧美日韩无卡精品| 国产在视频线在精品| 一a级毛片在线观看| 日韩大尺度精品在线看网址| 少妇裸体淫交视频免费看高清| 不卡视频在线观看欧美| 男人的好看免费观看在线视频| 亚洲人成网站在线播放欧美日韩| 中亚洲国语对白在线视频| 精品不卡国产一区二区三区| 俺也久久电影网| 亚洲 国产 在线| 精品午夜福利视频在线观看一区| 香蕉av资源在线| 欧美最黄视频在线播放免费| 精品人妻视频免费看| а√天堂www在线а√下载| 亚洲色图av天堂| 一边摸一边抽搐一进一小说| 一级黄色大片毛片| 男女视频在线观看网站免费| 国产精品一区二区三区四区免费观看 | 欧美高清性xxxxhd video| 国产白丝娇喘喷水9色精品| 一本久久中文字幕| 久久久久久九九精品二区国产| 日本撒尿小便嘘嘘汇集6| 午夜福利欧美成人| 日本在线视频免费播放| 亚洲av成人av| 在线播放无遮挡| 欧美黑人巨大hd| 我的老师免费观看完整版| 久久久色成人| 国产真实伦视频高清在线观看 | 午夜福利在线在线| av黄色大香蕉| 亚洲国产欧美人成|