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

    基于SSWPT的轉(zhuǎn)子油膜失穩(wěn)故障分析

    2021-05-22 08:44:48唐貴基徐振麗
    動(dòng)力工程學(xué)報(bào) 2021年5期
    關(guān)鍵詞:渦動(dòng)波包時(shí)頻

    唐貴基,徐振麗

    (華北電力大學(xué) 河北省電力機(jī)械裝備健康維護(hù)與失效預(yù)防重點(diǎn)實(shí)驗(yàn)室,河北保定 071003)

    滑動(dòng)軸承作為旋轉(zhuǎn)機(jī)械的關(guān)鍵部件,因其能夠有效增強(qiáng)轉(zhuǎn)子系統(tǒng)的支撐能力而被廣泛使用[1]。當(dāng)轉(zhuǎn)子受到外界因素的干擾時(shí),軸承的油膜力與負(fù)載不平衡,使軸頸圍繞平衡位置旋轉(zhuǎn),從而產(chǎn)生油膜渦動(dòng)和油膜振蕩[2]。一旦轉(zhuǎn)子系統(tǒng)發(fā)生油膜振蕩,會(huì)嚴(yán)重?fù)p壞轉(zhuǎn)子,影響整個(gè)機(jī)械系統(tǒng)的正常運(yùn)行,因此研究轉(zhuǎn)子系統(tǒng)油膜失穩(wěn)故障對(duì)確保轉(zhuǎn)子正常安全運(yùn)行具有重要意義。實(shí)際工況下的轉(zhuǎn)子故障信號(hào)通常是非線性、非平穩(wěn)的多分量信號(hào),它們的頻率成分隨著時(shí)間變化。時(shí)頻分析方法是處理非平穩(wěn)信號(hào)的有效工具,被廣泛用于轉(zhuǎn)子故障診斷中[3-4]。目前,常用的時(shí)頻分析方法有短時(shí)傅里葉變換[5](Short-time Fourier Transform,STFT)、小波變換[6](Wavelet Transform,WT)、Wigner-Ville分布以及希爾伯特黃變換[7](Hilbert-Huang Transform,HHT)等。但是上述時(shí)頻分析方法存在一定的不足:STFT由于受時(shí)頻窗的影響,其高頻和低頻信號(hào)的分辨率不能同時(shí)兼顧;小波變換的時(shí)頻圖存在能量擴(kuò)散現(xiàn)象,使得頻率分辨率較低;Wigner-Ville分布由于受交叉項(xiàng)的干擾,限制了該方法的應(yīng)用。HHT是美國(guó)華裔科學(xué)家Huang N E在1998年提出的一種時(shí)頻分析方法,主要包括經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition,EMD)和希爾伯特變換兩部分。HHT方法具有良好的自適應(yīng)性,已被廣泛應(yīng)用于旋轉(zhuǎn)機(jī)械的故障診斷中[8-9]。然而,EMD方法存在模態(tài)混疊、端點(diǎn)效應(yīng)等問(wèn)題,嚴(yán)重影響了時(shí)頻分析的精度,對(duì)故障特征診斷造成了干擾。針對(duì)EMD方法的不足,近年來(lái)一些學(xué)者提出了變分模態(tài)分解(Variational Modal Decomposition,VMD)及經(jīng)驗(yàn)小波變換(Empirical Wavelet Transform,EWT)等分解方法,可用于旋轉(zhuǎn)機(jī)械故障診斷中,取得了豐碩的成果[10-12]。但是,VMD方法需要人為設(shè)定分解層數(shù)及懲罰因子,EWT方法對(duì)參數(shù)的選取較復(fù)雜,若參數(shù)設(shè)置不當(dāng)將導(dǎo)致故障特征提取不準(zhǔn)確。因此,上述時(shí)頻分析方法難以提取轉(zhuǎn)子在實(shí)際運(yùn)行工況下的故障特征。

    同步壓縮變換(Synchrosqueezing Transform,SST)是Daubechies等[13]提出的一種先進(jìn)的時(shí)頻后處理技術(shù),本質(zhì)是時(shí)頻重排與連續(xù)小波變換的結(jié)合,將小波變換得到的時(shí)間-尺度系數(shù)重新分配壓縮至?xí)r間-頻率軸上。SST對(duì)于非平穩(wěn)信號(hào)具有較高的頻率分辨率,是故障診斷領(lǐng)域的一個(gè)研究熱點(diǎn),目前已成功應(yīng)用到軸承故障診斷[14]。同步壓縮小波包變換[15](Synchrosqueezed Wave Packet Transform,SSWPT)是最近發(fā)展起來(lái)的一種新的時(shí)頻分析方法,該方法是SST與小波包變換的結(jié)合,能夠準(zhǔn)確提取信號(hào)的瞬時(shí)信息,比SST具有更高的時(shí)頻分辨率及更好的抗噪性能。目前,該方法在醫(yī)學(xué)信號(hào)分析方面取得了一定的成果。筆者將同步壓縮小波包變換引進(jìn)到轉(zhuǎn)子油膜失穩(wěn)故障分析中,提出了一種基于SSWPT的轉(zhuǎn)子油膜失穩(wěn)故障分析方法,在仿真信號(hào)、實(shí)驗(yàn)信號(hào)以及現(xiàn)場(chǎng)數(shù)據(jù)分析中取得了良好的效果。

    1 同步壓縮小波包變換原理

    給定一個(gè)母波包w(ξ),使其滿足:

    (1)

    其中,|ξ|>1,ε>0,m為非負(fù)整數(shù)。

    (2)

    式中:a為尺度因子;b為平移因子。

    小波包族的傅里葉變換表示為:

    (3)

    (4)

    式中:Wf(a,b)為小波包變換系數(shù);|a|>1,b∈R。

    信號(hào)通過(guò)小波包變換后可以得到小波包變換系數(shù),此時(shí)可求取瞬時(shí)頻率vf(a,b):

    (5)

    最后根據(jù)計(jì)算得到的瞬時(shí)頻率,建立(b,a)→(b,vf(a,b))的映射關(guān)系,將小波包變換系數(shù)從時(shí)間-尺度平面轉(zhuǎn)化到時(shí)間-頻率平面,使得信號(hào)中各頻率成分在頻域方向進(jìn)行壓縮,提高了小波包在頻域上的分辨率,使信號(hào)中各頻率分量清晰地顯示在時(shí)頻圖上,則同步壓縮小波包變換的表達(dá)式Tf(v,b)可以表示為:

    (6)

    其中,δ表示狄拉克函數(shù),Rvf(a,b)表示瞬時(shí)頻率的實(shí)部;v為信號(hào)的中心頻率。

    綜上所述,基于同步壓縮小波包變換的步驟如下:(1)利用小波包變換求取小波包變換系數(shù)Wf(a,b)和主偏導(dǎo)數(shù)?bWf(a,b);(2)計(jì)算瞬時(shí)頻率vf(a,b),并進(jìn)行同步壓縮得到Tf(v,b)。

    2 仿真信號(hào)分析

    為驗(yàn)證SSWPT的時(shí)頻分析能力,現(xiàn)構(gòu)造如下2組仿真信號(hào)x1(t)和x2(t):

    (7)

    (8)

    式中:x11(t)為調(diào)頻信號(hào);x12(t)為調(diào)幅調(diào)頻信號(hào);x13(t)為隨機(jī)噪聲;x21(t)、x22(t)為余弦信號(hào);x23(t)為隨機(jī)噪聲;fr1、fr2為轉(zhuǎn)頻;f1、f2和f3為中心頻率;N1、N2為信號(hào)點(diǎn)數(shù);t為時(shí)間;randn為隨機(jī)函數(shù),其中仿真信號(hào)的采樣頻率為1 024 Hz,采樣點(diǎn)數(shù)為2 048。

    仿真信號(hào)x1(t)由調(diào)頻信號(hào)x11(t)、調(diào)幅調(diào)頻信號(hào)x12(t)和隨機(jī)噪聲x13(t)3個(gè)分量疊加而成。其中fr1=2 Hz,f1=50 Hz,f2=100 Hz,f3=5 Hz。仿真信號(hào)x2(t)是1個(gè)由3個(gè)分量構(gòu)成的隨時(shí)間變化的非平穩(wěn)信號(hào),式(8)中包含轉(zhuǎn)頻fr2(t)=30+30t和3倍轉(zhuǎn)頻的諧波分量3fr2。

    x1(t)、x2(t)的時(shí)域波形圖及頻譜圖如圖1和圖2所示,其中f為頻率。從圖2(a)可以看出,仿真信號(hào)x1(t)存在2個(gè)調(diào)制頻帶,而2個(gè)中心頻率不易被識(shí)別。由圖2(b)可知,仿真信號(hào)x2(t)的頻譜存在較寬的頻帶,出現(xiàn)頻率涂抹現(xiàn)象。因此,僅依靠頻譜分析不能有效提取信號(hào)的特征頻率。

    (a) x1(t)的時(shí)域波形圖

    (b) x2(t)的時(shí)域波形圖

    (a) x1(t)的頻譜圖

    (b) x2(t)的頻譜圖

    分別對(duì)仿真信號(hào)x1(t)和x2(t)進(jìn)行同步壓縮小波包變換,結(jié)果如圖3所示。由圖3(a)可知,SSWPT方法能夠精確地把仿真信號(hào)x1(t)的頻率成分有效分離。由圖3(b)可知,SSWPT方法能夠檢測(cè)到各頻率成分隨時(shí)間的變化趨勢(shì)。

    (a) x1(t)的SSWPT時(shí)頻圖

    (b) x2(t)的SSWPT時(shí)頻圖

    圖4(a)和圖4(b)分別給出了HHT方法對(duì)仿真信號(hào)x1(t)和x2(t)的分析結(jié)果。由圖4(a)可知,HHT方法提取的頻率分量之間出現(xiàn)了不同程度的模態(tài)混疊現(xiàn)象,不能精確提取信號(hào)的頻率成分。從圖4(b)可以觀察到頻率分量在高頻部分聚集性較低。

    (a) x1(t)的HHT時(shí)頻圖

    (b) x2(t)的HHT時(shí)頻圖

    3 實(shí)驗(yàn)分析

    實(shí)驗(yàn)裝置采用圖5所示的Bently-RK4轉(zhuǎn)子實(shí)驗(yàn)臺(tái),分別進(jìn)行了恒定轉(zhuǎn)速和變轉(zhuǎn)速工況下的油膜失穩(wěn)故障實(shí)驗(yàn)。振動(dòng)信號(hào)通過(guò)電渦流傳感器采集,信號(hào)采樣頻率為1 280 Hz。

    圖5 Bently-RK4轉(zhuǎn)子實(shí)驗(yàn)臺(tái)

    3.1 恒定轉(zhuǎn)速工況下油膜渦動(dòng)故障分析

    首先將轉(zhuǎn)子升速至實(shí)驗(yàn)轉(zhuǎn)速2 800 r/min模擬油膜渦動(dòng)故障。油膜渦動(dòng)是轉(zhuǎn)子受到軸承油膜力的作用而產(chǎn)生的渦動(dòng),是一種自激振動(dòng)。當(dāng)轉(zhuǎn)子系統(tǒng)處于油膜渦動(dòng)狀態(tài)時(shí),其振動(dòng)信號(hào)中含有基頻及油膜渦動(dòng)頻率,由于渦動(dòng)頻率與基頻的二分頻接近,油膜渦動(dòng)又稱(chēng)為半頻渦動(dòng)。轉(zhuǎn)子油膜渦動(dòng)信號(hào)的時(shí)域波形及頻譜如圖6所示。由圖6可知,轉(zhuǎn)子油膜渦動(dòng)的時(shí)域波形發(fā)生變形。在頻譜圖中,二分頻振幅fo較為突出,且圖中出現(xiàn)了振幅較低的高次諧波成分。因此,僅根據(jù)頻譜分析不能有效診斷轉(zhuǎn)子是否發(fā)生了油膜渦動(dòng)故障。

    (a) 油膜渦動(dòng)信號(hào)的時(shí)域波形圖

    (b) 油膜渦動(dòng)信號(hào)的頻譜圖

    為準(zhǔn)確提取轉(zhuǎn)子故障特征,對(duì)振動(dòng)信號(hào)進(jìn)行SSWPT分析,結(jié)果見(jiàn)圖7。由圖7可知,SSWPT方法能夠準(zhǔn)確提取轉(zhuǎn)子油膜渦動(dòng)故障的特征頻率,時(shí)頻圖中能量高度集中在油膜渦動(dòng)頻率附近,且基頻存在明顯的調(diào)頻現(xiàn)象,而這些現(xiàn)象在頻譜圖中是觀察不到的。

    圖7 油膜渦動(dòng)信號(hào)的SSWPT時(shí)頻圖

    為突出SSWPT方法的優(yōu)越性,采用HHT方法對(duì)上述油膜渦動(dòng)信號(hào)進(jìn)行分析,結(jié)果如圖8所示。

    圖8 油膜渦動(dòng)信號(hào)的HHT時(shí)頻圖

    由圖8可知,通過(guò)HHT方法得到的時(shí)頻圖精度較低,能量擴(kuò)散,且模態(tài)混疊現(xiàn)象嚴(yán)重,與SSWPT的時(shí)頻圖相比,HHT時(shí)頻圖不能充分體現(xiàn)油膜渦動(dòng)故障的特征頻率。

    3.2 變轉(zhuǎn)速工況下油膜渦動(dòng)故障分析

    通過(guò)控制轉(zhuǎn)速調(diào)節(jié)裝置,使轉(zhuǎn)子轉(zhuǎn)速均勻上升,模擬油膜渦動(dòng)的演變過(guò)程。轉(zhuǎn)子的轉(zhuǎn)速由1 500 r/min上升至4 700 r/min,歷時(shí)100 s,采樣頻率為1 280 Hz。由于故障數(shù)據(jù)點(diǎn)數(shù)較多,為提高計(jì)算效率,對(duì)故障數(shù)據(jù)進(jìn)行降采樣處理,設(shè)置降采樣率為4,通過(guò)計(jì)算可以得到降采樣頻率為320 Hz。圖9為轉(zhuǎn)子升速過(guò)程中的振動(dòng)信號(hào)時(shí)域波形及頻譜圖。由圖9可知,轉(zhuǎn)子發(fā)生油膜渦動(dòng)時(shí),其過(guò)程比較復(fù)雜,很難從時(shí)域圖中發(fā)現(xiàn)油膜渦動(dòng)的規(guī)律,且頻譜圖中存在頻率模糊現(xiàn)象,不易識(shí)別油膜渦動(dòng)故障的頻率特征。

    (a) 油膜渦動(dòng)過(guò)程的時(shí)域波形圖

    (b) 油膜渦動(dòng)過(guò)程的頻譜圖

    為準(zhǔn)確分析轉(zhuǎn)子升速過(guò)程中的振動(dòng)信號(hào),對(duì)該振動(dòng)信號(hào)進(jìn)行SSWPT分析,結(jié)果產(chǎn)生2個(gè)頻率成分,分別為基頻和渦動(dòng)頻率,如圖10(a)所示。由圖10(a)可知,當(dāng)轉(zhuǎn)速較低時(shí),時(shí)頻圖的低頻處只存在基頻分量,說(shuō)明轉(zhuǎn)子系統(tǒng)沒(méi)有發(fā)生渦動(dòng)現(xiàn)象。10 s左右開(kāi)始出現(xiàn)半頻分量,油膜渦動(dòng)從此刻開(kāi)始;15 s左右油膜渦動(dòng)頻率出現(xiàn)振蕩現(xiàn)象,且隨著時(shí)間的變化基頻出現(xiàn)調(diào)頻現(xiàn)象;75 s后當(dāng)轉(zhuǎn)子轉(zhuǎn)速接近2倍臨界轉(zhuǎn)速時(shí),由于油膜渦動(dòng)頻率與轉(zhuǎn)子1倍臨界轉(zhuǎn)速的基頻重合,產(chǎn)生油膜振蕩,此時(shí)隨著轉(zhuǎn)子轉(zhuǎn)速的增加,油膜渦動(dòng)頻率不再增加,轉(zhuǎn)子系統(tǒng)處于油膜振蕩狀態(tài)。圖10(b)為HHT方法的分析結(jié)果,HHT時(shí)頻圖不易觀察到油膜渦動(dòng)頻率隨時(shí)間的變化規(guī)律。綜上所述,SSWPT方法可以成功提取油膜失穩(wěn)故障在變轉(zhuǎn)速工況下的演化特征。

    (a) 油膜渦動(dòng)過(guò)程的SSWPT時(shí)頻圖

    (b) 油膜渦動(dòng)過(guò)程的HHT時(shí)頻圖

    4 工程應(yīng)用

    采用某廠汽輪機(jī)轉(zhuǎn)子油膜渦動(dòng)故障數(shù)據(jù)進(jìn)行分析驗(yàn)證。該汽輪機(jī)為超臨界、一次中間再熱、單軸、雙缸雙排汽、凝汽式汽輪機(jī),型號(hào)C350/267-24.2/0.4/566/566。汽輪發(fā)電機(jī)組軸系共有7個(gè)支撐軸瓦,其中1號(hào)、2號(hào)軸瓦支撐高中壓轉(zhuǎn)子,3號(hào)、4號(hào)軸瓦支撐低壓轉(zhuǎn)子,5號(hào)、6號(hào)軸瓦支撐發(fā)電機(jī)轉(zhuǎn)子,7號(hào)軸瓦為滑環(huán)軸轉(zhuǎn)子外伸端的支撐軸瓦,滑環(huán)軸轉(zhuǎn)子為發(fā)電機(jī)轉(zhuǎn)子的外伸,各轉(zhuǎn)子之間均采用剛性聯(lián)軸器連接,軸系結(jié)構(gòu)圖見(jiàn)圖11。

    圖11 機(jī)組軸系結(jié)構(gòu)圖

    汽輪機(jī)的1號(hào)軸承在機(jī)組運(yùn)行期間經(jīng)常出現(xiàn)較大的不穩(wěn)定振動(dòng)。為分析機(jī)組故障,采用靈敏度為7.87 mV/μm的電渦流傳感器測(cè)量1號(hào)軸承的徑向振動(dòng),采集信號(hào)時(shí)機(jī)組處于額定工況附近。轉(zhuǎn)子轉(zhuǎn)速為3 000 r/min,采樣頻率為6 400 Hz,選取2 048個(gè)振動(dòng)數(shù)據(jù)點(diǎn)進(jìn)行分析,其波形及頻譜如圖12所示。從時(shí)域圖可知,信號(hào)波形發(fā)生周期性變化,頻譜圖中存在轉(zhuǎn)頻及二分頻,轉(zhuǎn)頻處幅值比較突出,但高次諧波幅值較低,僅根據(jù)頻譜圖診斷該機(jī)組轉(zhuǎn)子故障類(lèi)型的效果欠佳。

    (a) 振動(dòng)信號(hào)波形圖

    (b) 振動(dòng)信號(hào)頻譜圖

    為提取該機(jī)組轉(zhuǎn)子的故障特征,并進(jìn)行準(zhǔn)確判斷。采用SSWPT方法進(jìn)行分析,結(jié)果如圖13所示。圖13能夠清晰反映出機(jī)組振動(dòng)信號(hào)的轉(zhuǎn)頻及二分頻,且在轉(zhuǎn)頻、二分頻及高頻處存在調(diào)頻現(xiàn)象。因此,可以判斷出轉(zhuǎn)子發(fā)生了油膜渦動(dòng)故障。

    圖13 振動(dòng)信號(hào)的SSWPT時(shí)頻圖

    對(duì)采集的振動(dòng)信號(hào)進(jìn)行HHT方法分析,結(jié)果如圖14所示。圖14中轉(zhuǎn)子故障特征頻率不明顯,高頻處存在大量的噪聲干擾。因此,HHT方法不能有效提取轉(zhuǎn)子故障特征頻率,判別轉(zhuǎn)子故障類(lèi)型。

    圖14 振動(dòng)信號(hào)的HHT時(shí)頻圖

    5 結(jié) 論

    (1) 模擬仿真及實(shí)驗(yàn)信號(hào)分析結(jié)果表明,與HHT方法相比,基于SSWPT的時(shí)頻分析方法具有更高的時(shí)頻精度,有效避免了模態(tài)混疊現(xiàn)象,能夠準(zhǔn)確提取油膜失穩(wěn)故障在恒定轉(zhuǎn)速工況下的故障特征頻率及在變轉(zhuǎn)速工況下的故障演化特征。

    (2) 現(xiàn)場(chǎng)故障數(shù)據(jù)分析結(jié)果表明,SSWPT方法能夠去除噪聲干擾,可以精確地提取出轉(zhuǎn)子的故障特征頻率,更準(zhǔn)確地診斷轉(zhuǎn)子油膜失穩(wěn)故障,為轉(zhuǎn)子油膜失穩(wěn)故障診斷提供了一種新的時(shí)頻分析方法。

    猜你喜歡
    渦動(dòng)波包時(shí)頻
    基于小波包Tsallis熵和RVM的模擬電路故障診斷
    BTA鉆桿渦動(dòng)數(shù)學(xué)建模及實(shí)驗(yàn)研究
    基于小波包變換的電力系統(tǒng)諧波分析
    理想條件下BTA鉆鉆桿的渦動(dòng)分析
    小波包理論與圖像小波包分解
    基于時(shí)頻分析的逆合成孔徑雷達(dá)成像技術(shù)
    對(duì)采樣數(shù)據(jù)序列進(jìn)行時(shí)頻分解法的改進(jìn)
    雙線性時(shí)頻分布交叉項(xiàng)提取及損傷識(shí)別應(yīng)用
    基于小波包的全信息解調(diào)方法及其應(yīng)用
    淺析《守望燈塔》中的時(shí)頻
    亚洲精品美女久久久久99蜜臀| 精品国产亚洲在线| 久久影院123| 国产97色在线日韩免费| 99riav亚洲国产免费| 国产在线免费精品| 免费看a级黄色片| 亚洲色图av天堂| 我的亚洲天堂| 日本撒尿小便嘘嘘汇集6| 亚洲七黄色美女视频| 深夜精品福利| 在线观看一区二区三区激情| 亚洲三区欧美一区| 精品卡一卡二卡四卡免费| 黄色a级毛片大全视频| 国产97色在线日韩免费| 一边摸一边抽搐一进一小说 | 精品国产国语对白av| 国产单亲对白刺激| 黄片大片在线免费观看| 久久中文字幕一级| 少妇被粗大的猛进出69影院| 国产在线一区二区三区精| 欧美日韩精品网址| 男男h啪啪无遮挡| 国产一区二区 视频在线| 两人在一起打扑克的视频| 两性夫妻黄色片| 天堂中文最新版在线下载| 十八禁网站网址无遮挡| 免费看a级黄色片| 亚洲三区欧美一区| 五月天丁香电影| 免费一级毛片在线播放高清视频 | 69av精品久久久久久 | 亚洲免费av在线视频| 亚洲av国产av综合av卡| 国产精品av久久久久免费| 国产91精品成人一区二区三区 | 后天国语完整版免费观看| 欧美日韩中文字幕国产精品一区二区三区 | 在线天堂中文资源库| 精品视频人人做人人爽| 成人国产一区最新在线观看| 精品欧美一区二区三区在线| 成人手机av| 99国产精品一区二区蜜桃av | 在线观看免费视频网站a站| 亚洲精品在线美女| 国产老妇伦熟女老妇高清| 久久热在线av| 午夜激情久久久久久久| 香蕉国产在线看| av超薄肉色丝袜交足视频| 极品教师在线免费播放| 十八禁网站网址无遮挡| 丰满饥渴人妻一区二区三| 精品少妇内射三级| 久久影院123| 久久这里只有精品19| 亚洲人成电影观看| 国产精品久久久av美女十八| 999久久久精品免费观看国产| 成年版毛片免费区| 天天躁日日躁夜夜躁夜夜| 免费在线观看日本一区| 国产精品免费大片| av片东京热男人的天堂| 亚洲五月婷婷丁香| 一夜夜www| 成人亚洲精品一区在线观看| 美女高潮喷水抽搐中文字幕| 精品一品国产午夜福利视频| 两个人免费观看高清视频| 久久久国产欧美日韩av| 飞空精品影院首页| 国产日韩欧美视频二区| 欧美性长视频在线观看| 黄色怎么调成土黄色| 欧美成人免费av一区二区三区 | 国产欧美日韩精品亚洲av| 80岁老熟妇乱子伦牲交| 欧美在线黄色| 三上悠亚av全集在线观看| 亚洲欧美日韩另类电影网站| 欧美性长视频在线观看| 国产成人免费无遮挡视频| 国产日韩一区二区三区精品不卡| 狂野欧美激情性xxxx| 人人妻人人添人人爽欧美一区卜| 欧美黄色淫秽网站| 99精品久久久久人妻精品| 午夜精品久久久久久毛片777| 一区二区三区精品91| 99久久国产精品久久久| 日韩三级视频一区二区三区| 国产伦理片在线播放av一区| 亚洲视频免费观看视频| 在线观看一区二区三区激情| 乱人伦中国视频| 日韩欧美国产一区二区入口| 老汉色∧v一级毛片| 亚洲国产成人一精品久久久| 国产精品久久久久久精品古装| 久热爱精品视频在线9| 97人妻天天添夜夜摸| 在线观看66精品国产| 美女国产高潮福利片在线看| www日本在线高清视频| 免费日韩欧美在线观看| 大片免费播放器 马上看| 黄色毛片三级朝国网站| 日日夜夜操网爽| 免费在线观看日本一区| 黄片播放在线免费| 午夜免费鲁丝| 黄色a级毛片大全视频| 91大片在线观看| 色尼玛亚洲综合影院| svipshipincom国产片| 欧美国产精品一级二级三级| 国产97色在线日韩免费| 可以免费在线观看a视频的电影网站| 亚洲国产欧美在线一区| 人妻久久中文字幕网| 国产男女内射视频| 欧美另类亚洲清纯唯美| 最新在线观看一区二区三区| 国产精品亚洲av一区麻豆| 少妇的丰满在线观看| 亚洲国产中文字幕在线视频| 夫妻午夜视频| av视频免费观看在线观看| 精品视频人人做人人爽| 少妇被粗大的猛进出69影院| 国产伦理片在线播放av一区| 99国产精品免费福利视频| 久久毛片免费看一区二区三区| 国产又色又爽无遮挡免费看| 国产精品亚洲一级av第二区| 亚洲中文字幕日韩| 又紧又爽又黄一区二区| 一级a爱视频在线免费观看| 国产精品98久久久久久宅男小说| 成年版毛片免费区| 亚洲精品成人av观看孕妇| 热re99久久精品国产66热6| 天天躁夜夜躁狠狠躁躁| 亚洲精品在线观看二区| 搡老岳熟女国产| 日韩熟女老妇一区二区性免费视频| 啪啪无遮挡十八禁网站| 丝袜人妻中文字幕| 交换朋友夫妻互换小说| 久久免费观看电影| 国产精品成人在线| 日本av免费视频播放| 欧美日韩视频精品一区| 欧美日韩福利视频一区二区| 日本vs欧美在线观看视频| 首页视频小说图片口味搜索| 亚洲精品国产一区二区精华液| 在线播放国产精品三级| 最新美女视频免费是黄的| 久久人人爽av亚洲精品天堂| 国产成人免费无遮挡视频| 免费在线观看视频国产中文字幕亚洲| 一级a爱视频在线免费观看| 国产一区二区三区视频了| 男人舔女人的私密视频| 亚洲美女黄片视频| 久久精品熟女亚洲av麻豆精品| 国产精品免费一区二区三区在线 | 在线 av 中文字幕| 国产精品免费视频内射| 变态另类成人亚洲欧美熟女 | 免费在线观看日本一区| 日韩欧美一区二区三区在线观看 | 午夜成年电影在线免费观看| 色婷婷久久久亚洲欧美| 啦啦啦中文免费视频观看日本| 色尼玛亚洲综合影院| 男女高潮啪啪啪动态图| 欧美大码av| 久久午夜综合久久蜜桃| 国产aⅴ精品一区二区三区波| 久久久久久久久免费视频了| 啦啦啦在线免费观看视频4| 一本综合久久免费| 黄色成人免费大全| 麻豆成人av在线观看| 在线十欧美十亚洲十日本专区| 在线 av 中文字幕| 丰满饥渴人妻一区二区三| 久久久久久人人人人人| 日韩欧美国产一区二区入口| 日韩成人在线观看一区二区三区| 国产日韩欧美在线精品| 亚洲精品久久成人aⅴ小说| 国产不卡一卡二| kizo精华| 99在线人妻在线中文字幕 | 另类亚洲欧美激情| 黄片小视频在线播放| 大香蕉久久网| 久久人妻熟女aⅴ| 最近最新免费中文字幕在线| 亚洲熟女精品中文字幕| 性少妇av在线| 欧美一级毛片孕妇| 两个人看的免费小视频| 午夜免费鲁丝| 国产成+人综合+亚洲专区| 黄色视频在线播放观看不卡| 悠悠久久av| 国产在线观看jvid| 久久亚洲精品不卡| 成人免费观看视频高清| cao死你这个sao货| 亚洲国产欧美在线一区| 亚洲欧洲日产国产| 免费一级毛片在线播放高清视频 | 国产av国产精品国产| 侵犯人妻中文字幕一二三四区| 妹子高潮喷水视频| 国产精品久久久av美女十八| 精品乱码久久久久久99久播| 国产精品一区二区在线观看99| 久久免费观看电影| 国产一区二区在线观看av| 免费久久久久久久精品成人欧美视频| 50天的宝宝边吃奶边哭怎么回事| 99九九在线精品视频| 精品少妇一区二区三区视频日本电影| 天天躁夜夜躁狠狠躁躁| 久久精品熟女亚洲av麻豆精品| 成人免费观看视频高清| 一级,二级,三级黄色视频| 最黄视频免费看| 精品午夜福利视频在线观看一区 | 欧美性长视频在线观看| 丰满饥渴人妻一区二区三| 欧美av亚洲av综合av国产av| 丝袜美腿诱惑在线| 操出白浆在线播放| 久久久久久久久免费视频了| 成人永久免费在线观看视频 | 亚洲国产精品一区二区三区在线| 99国产精品一区二区蜜桃av | 69精品国产乱码久久久| 亚洲成国产人片在线观看| 精品少妇久久久久久888优播| 黄色视频在线播放观看不卡| 日韩 欧美 亚洲 中文字幕| 亚洲熟妇熟女久久| 老鸭窝网址在线观看| 国产区一区二久久| 黄色a级毛片大全视频| 欧美日韩精品网址| 精品亚洲成国产av| 日本黄色日本黄色录像| 免费看十八禁软件| 精品久久蜜臀av无| 91精品国产国语对白视频| 久久久久久免费高清国产稀缺| 国产真人三级小视频在线观看| 国产三级黄色录像| 香蕉国产在线看| 国产精品 欧美亚洲| 亚洲国产欧美网| 精品一区二区三区视频在线观看免费 | 999久久久精品免费观看国产| 日韩 欧美 亚洲 中文字幕| 国产成人精品在线电影| 国产av又大| 久久午夜综合久久蜜桃| www日本在线高清视频| 午夜免费成人在线视频| 捣出白浆h1v1| 亚洲av国产av综合av卡| 久久久国产一区二区| av电影中文网址| 天天躁夜夜躁狠狠躁躁| 亚洲av美国av| a级毛片在线看网站| 日韩中文字幕视频在线看片| 热re99久久国产66热| 久久久久久人人人人人| 午夜福利,免费看| 日韩人妻精品一区2区三区| 亚洲第一欧美日韩一区二区三区 | 纯流量卡能插随身wifi吗| 男女床上黄色一级片免费看| 又大又爽又粗| 亚洲欧洲精品一区二区精品久久久| 欧美日韩av久久| 国产xxxxx性猛交| 极品教师在线免费播放| 中文字幕人妻熟女乱码| 午夜免费鲁丝| 婷婷成人精品国产| 在线永久观看黄色视频| 交换朋友夫妻互换小说| 国产精品亚洲av一区麻豆| 怎么达到女性高潮| 亚洲五月色婷婷综合| 国产淫语在线视频| 性色av乱码一区二区三区2| 69精品国产乱码久久久| 国产精品一区二区免费欧美| 女人爽到高潮嗷嗷叫在线视频| 国产精品久久久久久精品古装| 久久精品熟女亚洲av麻豆精品| 中文字幕另类日韩欧美亚洲嫩草| 久久性视频一级片| 免费av中文字幕在线| 在线十欧美十亚洲十日本专区| 免费看a级黄色片| 欧美激情极品国产一区二区三区| 免费人妻精品一区二区三区视频| 国产亚洲午夜精品一区二区久久| 在线观看舔阴道视频| 好男人电影高清在线观看| 久久久水蜜桃国产精品网| 中文欧美无线码| 啦啦啦中文免费视频观看日本| 亚洲色图av天堂| 美女午夜性视频免费| 69av精品久久久久久 | 黄网站色视频无遮挡免费观看| 80岁老熟妇乱子伦牲交| 国产福利在线免费观看视频| 女人精品久久久久毛片| 天堂8中文在线网| 午夜激情av网站| 国产xxxxx性猛交| 国产免费福利视频在线观看| 中文字幕色久视频| 免费观看人在逋| 国产亚洲欧美精品永久| 麻豆成人av在线观看| 国产伦理片在线播放av一区| 免费人妻精品一区二区三区视频| 2018国产大陆天天弄谢| www日本在线高清视频| 麻豆成人av在线观看| av天堂在线播放| 精品久久久久久久毛片微露脸| 国产精品亚洲av一区麻豆| 精品欧美一区二区三区在线| 国产日韩一区二区三区精品不卡| 热re99久久国产66热| a级毛片黄视频| 波多野结衣av一区二区av| 亚洲人成77777在线视频| 免费在线观看黄色视频的| 50天的宝宝边吃奶边哭怎么回事| 欧美日韩黄片免| 久久久久精品人妻al黑| 成年人黄色毛片网站| 脱女人内裤的视频| 欧美日韩亚洲高清精品| 叶爱在线成人免费视频播放| 91麻豆av在线| 97人妻天天添夜夜摸| 国产高清激情床上av| 电影成人av| www.熟女人妻精品国产| 一边摸一边做爽爽视频免费| 91麻豆精品激情在线观看国产 | xxxhd国产人妻xxx| 亚洲精品中文字幕在线视频| 美女高潮喷水抽搐中文字幕| 久久香蕉激情| 亚洲国产欧美网| 一个人免费看片子| 亚洲国产精品一区二区三区在线| 桃花免费在线播放| 夜夜骑夜夜射夜夜干| 国产野战对白在线观看| 精品免费久久久久久久清纯 | 久久香蕉激情| 精品少妇内射三级| 亚洲成人免费av在线播放| 狠狠婷婷综合久久久久久88av| 可以免费在线观看a视频的电影网站| 日韩免费高清中文字幕av| 丁香欧美五月| cao死你这个sao货| 丰满人妻熟妇乱又伦精品不卡| 在线观看免费视频日本深夜| 考比视频在线观看| 啪啪无遮挡十八禁网站| 国产精品一区二区在线不卡| 黑人巨大精品欧美一区二区mp4| 亚洲久久久国产精品| 女警被强在线播放| 亚洲午夜精品一区,二区,三区| 精品视频人人做人人爽| 久久免费观看电影| 国产成人啪精品午夜网站| 中文亚洲av片在线观看爽 | 午夜福利乱码中文字幕| 十八禁人妻一区二区| 啦啦啦中文免费视频观看日本| av线在线观看网站| 久久人妻熟女aⅴ| 大码成人一级视频| 亚洲专区国产一区二区| 美女高潮喷水抽搐中文字幕| 无人区码免费观看不卡 | 美女高潮到喷水免费观看| 精品一区二区三区四区五区乱码| 一本久久精品| 国产淫语在线视频| 亚洲熟女精品中文字幕| 久久 成人 亚洲| 国产深夜福利视频在线观看| 亚洲欧美一区二区三区黑人| √禁漫天堂资源中文www| 激情视频va一区二区三区| 精品人妻熟女毛片av久久网站| 日本wwww免费看| 欧美国产精品va在线观看不卡| 久热这里只有精品99| 如日韩欧美国产精品一区二区三区| 国精品久久久久久国模美| 极品少妇高潮喷水抽搐| 夜夜爽天天搞| 久久青草综合色| 久久精品国产a三级三级三级| 亚洲av第一区精品v没综合| 法律面前人人平等表现在哪些方面| 中文亚洲av片在线观看爽 | 欧美午夜高清在线| 男人操女人黄网站| 免费一级毛片在线播放高清视频 | 涩涩av久久男人的天堂| 最近最新中文字幕大全免费视频| 91成人精品电影| 国产xxxxx性猛交| 桃红色精品国产亚洲av| 久久国产亚洲av麻豆专区| 日韩制服丝袜自拍偷拍| 亚洲久久久国产精品| 精品一区二区三区四区五区乱码| 国产成+人综合+亚洲专区| 国产主播在线观看一区二区| 亚洲伊人久久精品综合| 久久精品亚洲精品国产色婷小说| 在线永久观看黄色视频| 99热国产这里只有精品6| 色婷婷久久久亚洲欧美| 99国产精品免费福利视频| 一本色道久久久久久精品综合| 国产日韩欧美视频二区| 大香蕉久久网| av网站在线播放免费| 一区二区三区激情视频| 老司机影院毛片| 老司机午夜十八禁免费视频| 国产精品电影一区二区三区 | 女人被躁到高潮嗷嗷叫费观| 亚洲七黄色美女视频| 人人澡人人妻人| 免费黄频网站在线观看国产| 1024视频免费在线观看| 91精品国产国语对白视频| 一区福利在线观看| 99久久国产精品久久久| 午夜激情av网站| 国产成人精品无人区| 免费一级毛片在线播放高清视频 | 亚洲精品一二三| 国产成人啪精品午夜网站| 好男人电影高清在线观看| 亚洲精品自拍成人| 久久 成人 亚洲| 欧美变态另类bdsm刘玥| e午夜精品久久久久久久| 亚洲中文字幕日韩| 欧美激情高清一区二区三区| 日本av手机在线免费观看| 一个人免费看片子| 18在线观看网站| 一本大道久久a久久精品| 亚洲精品美女久久av网站| 欧美日韩亚洲国产一区二区在线观看 | 免费观看av网站的网址| 天堂8中文在线网| videos熟女内射| 嫁个100分男人电影在线观看| 成人永久免费在线观看视频 | 中文字幕最新亚洲高清| 丝袜在线中文字幕| 日韩欧美免费精品| 日日夜夜操网爽| 成年动漫av网址| 欧美黄色淫秽网站| 老汉色∧v一级毛片| a级毛片在线看网站| 成人精品一区二区免费| 久久中文看片网| 亚洲avbb在线观看| 国产日韩欧美亚洲二区| 香蕉国产在线看| 国产精品1区2区在线观看. | 亚洲国产欧美日韩在线播放| cao死你这个sao货| 一区二区三区乱码不卡18| 美女高潮喷水抽搐中文字幕| 午夜激情久久久久久久| 男人舔女人的私密视频| 汤姆久久久久久久影院中文字幕| 91九色精品人成在线观看| 国产午夜精品久久久久久| 国产片内射在线| 午夜福利欧美成人| 一本—道久久a久久精品蜜桃钙片| 黄色 视频免费看| 亚洲国产欧美一区二区综合| 亚洲自偷自拍图片 自拍| 90打野战视频偷拍视频| 一级毛片女人18水好多| 日韩欧美一区视频在线观看| 欧美成人免费av一区二区三区 | 久久久久精品国产欧美久久久| 色综合婷婷激情| av天堂在线播放| 亚洲av电影在线进入| 男女边摸边吃奶| 日本vs欧美在线观看视频| 777米奇影视久久| 高清欧美精品videossex| 欧美精品亚洲一区二区| 18禁美女被吸乳视频| 在线观看舔阴道视频| 国产一区二区激情短视频| 精品国产乱子伦一区二区三区| 日韩欧美三级三区| 国产aⅴ精品一区二区三区波| 精品一品国产午夜福利视频| 日韩中文字幕视频在线看片| 97在线人人人人妻| 久久久精品国产亚洲av高清涩受| 老汉色∧v一级毛片| 欧美精品一区二区免费开放| 亚洲精品成人av观看孕妇| 电影成人av| 侵犯人妻中文字幕一二三四区| 亚洲av欧美aⅴ国产| 欧美日韩av久久| 欧美在线一区亚洲| 成人永久免费在线观看视频 | 精品少妇久久久久久888优播| 精品乱码久久久久久99久播| 国产亚洲午夜精品一区二区久久| 国产一区有黄有色的免费视频| 国产区一区二久久| 丁香六月天网| 日韩大片免费观看网站| 国产野战对白在线观看| 桃花免费在线播放| 女人高潮潮喷娇喘18禁视频| 色婷婷av一区二区三区视频| 亚洲av成人不卡在线观看播放网| 精品熟女少妇八av免费久了| 日韩免费av在线播放| 欧美日韩福利视频一区二区| 亚洲精品国产区一区二| 香蕉国产在线看| 动漫黄色视频在线观看| 久久精品亚洲av国产电影网| 日本黄色视频三级网站网址 | 久久精品成人免费网站| 97人妻天天添夜夜摸| 国产精品 国内视频| av免费在线观看网站| 亚洲精华国产精华精| 在线观看66精品国产| 一边摸一边抽搐一进一小说 | 一边摸一边抽搐一进一出视频| 国产福利在线免费观看视频| 曰老女人黄片| 久久天躁狠狠躁夜夜2o2o| 在线观看免费视频日本深夜| 超色免费av| 久久久精品区二区三区| 成人黄色视频免费在线看| 三上悠亚av全集在线观看| 精品一区二区三区四区五区乱码| 国产男靠女视频免费网站| 最新在线观看一区二区三区| 亚洲国产毛片av蜜桃av| 人妻久久中文字幕网| 少妇精品久久久久久久| √禁漫天堂资源中文www| 国产91精品成人一区二区三区 | h视频一区二区三区| 精品午夜福利视频在线观看一区 | 久久精品国产亚洲av高清一级| 五月天丁香电影| 色老头精品视频在线观看| 免费观看人在逋| 下体分泌物呈黄色| 丝袜人妻中文字幕| 国产一区二区三区视频了| 日本一区二区免费在线视频| av天堂在线播放| 成年女人毛片免费观看观看9 | 国产精品熟女久久久久浪| 美女高潮到喷水免费观看| 最近最新中文字幕大全电影3 | 黄片小视频在线播放|