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

    基于階次分析的風(fēng)機(jī)增速箱故障診斷

    2016-03-21 02:31:54陳長(zhǎng)征王海童孫自強(qiáng)
    重型機(jī)械 2016年2期
    關(guān)鍵詞:角域階次時(shí)域

    陳長(zhǎng)征,王海童,孫自強(qiáng)

    (沈陽(yáng)工業(yè)大學(xué) 機(jī)械工程學(xué)院,遼寧 沈陽(yáng) 110178)

    ·故障診斷·

    基于階次分析的風(fēng)機(jī)增速箱故障診斷

    陳長(zhǎng)征,王海童,孫自強(qiáng)

    (沈陽(yáng)工業(yè)大學(xué) 機(jī)械工程學(xué)院,遼寧 沈陽(yáng) 110178)

    風(fēng)力發(fā)電機(jī)組工作環(huán)境惡劣,振動(dòng)信號(hào)往往具有非平穩(wěn)的特性。為了提取故障信息的特征頻率,提出了基于重采樣的階次分析和小波閾值去噪相結(jié)合的診斷方法。該方法與經(jīng)典的計(jì)算階次跟蹤算法(COT)相比較,提高了重采樣精度,與傳統(tǒng)的FFT方法相比較,消除了因轉(zhuǎn)速不穩(wěn)定引起的“頻率模糊”現(xiàn)象。并通過(guò)風(fēng)場(chǎng)實(shí)驗(yàn)表明該方法有效,具有工程應(yīng)有價(jià)值。

    風(fēng)力發(fā)電機(jī);小波去噪;階次分析;重采樣

    0 前言

    齒輪箱作為風(fēng)機(jī)最主要的結(jié)構(gòu)之一,其發(fā)生故障的頻率最高,造成的停機(jī)時(shí)間也最長(zhǎng)。同時(shí)由于受風(fēng)機(jī)自身結(jié)構(gòu)和風(fēng)速時(shí)刻變化的雙重影響,所采集到風(fēng)機(jī)齒輪箱的故障信號(hào)往往具有非平穩(wěn)的特性,而傳統(tǒng)的基于FFT的普通頻譜分析方法只能提取機(jī)械在轉(zhuǎn)速穩(wěn)定狀況下的特征頻率。若將此方法用于提取風(fēng)機(jī)故障時(shí)的特征頻率,則將出現(xiàn)難以辨認(rèn)的“頻率模糊”現(xiàn)象[1-6]。階次分析是以角域采樣的理論為基礎(chǔ)建立起來(lái)的,其重要的思想就是將振動(dòng)信號(hào)從時(shí)域轉(zhuǎn)換到角域,因此階次分析的關(guān)鍵技術(shù)就是實(shí)現(xiàn)等角度采樣[7]。要實(shí)現(xiàn)對(duì)風(fēng)機(jī)旋轉(zhuǎn)不穩(wěn)定狀態(tài)下的角域采樣,可利用階次跟蹤的方法,常用的階次跟蹤方法有兩種:一種是硬件階次跟蹤方法,一種是計(jì)算階次跟蹤方法[8]。硬件階次跟蹤不僅成本高而且在風(fēng)機(jī)上安裝比較困難。計(jì)算階次跟蹤可利用數(shù)值插值的方法實(shí)現(xiàn)角域采樣,降低了成本。而基于等角度重采樣的階次分析是最常用的一種階次跟蹤方法但在實(shí)際工程中,由于受各種因素的影響,所采集到振動(dòng)信號(hào)信息不可避免的會(huì)存在各種噪聲,從而會(huì)影響重采樣的精度,本文采用小波硬閾值去噪的方法對(duì)其重采樣后的角域信號(hào)進(jìn)行去噪處理,再進(jìn)行FFT變換得到階次譜,通過(guò)和經(jīng)典的階次跟蹤方法(COT)相比較,提高了其重采樣精度。

    1 基于階次分析和小波閾值去噪處理非平穩(wěn)信號(hào)的原理

    1.1 階次分析和經(jīng)典的計(jì)算階次跟蹤(COT)

    (1)

    式中,O為階次;F為振動(dòng)信號(hào)的頻率;n為參考軸的轉(zhuǎn)速。

    為了能得到等間隔角度的時(shí)間采樣點(diǎn),假定在短時(shí)間內(nèi),轉(zhuǎn)軸的角加速度是保持不變的,轉(zhuǎn)過(guò)的角度和時(shí)間的關(guān)系可表示為

    φ(t)=b0+b1t+b2t2

    式中,t為時(shí)間點(diǎn);并假設(shè)每2個(gè)脈沖之間的相位差為Δφ;b0~b2表示未知系數(shù),則

    由式(2)、式(3)可得

    (4)

    式中,Δφ為采樣角度;k為插值系數(shù),其與采樣階次之間的關(guān)系為

    (5)

    根據(jù)式(4)求出的時(shí)間點(diǎn),可得到等角度間信號(hào)序列,從而實(shí)現(xiàn)時(shí)域信號(hào)到角域信號(hào)的轉(zhuǎn)變,對(duì)重采樣后的信號(hào)進(jìn)行FFT變換,得到信號(hào)的階次譜,并依據(jù)齒輪階次和頻率的關(guān)系,可以得出故障信號(hào)幅值與階次的函數(shù)關(guān)系[9]。

    在定量選礦廢水中加入次氯酸鈣粉末、次氯酸鈉溶液、雙氧水、高錳酸鉀、氯酸鉀等氧化劑各水體的COD濃度變化情況見(jiàn)表2。

    1.2 小波閾值去噪

    小波閾值去噪算法是一種對(duì)小波系數(shù)進(jìn)行非線性處理的方法,通過(guò)設(shè)置一個(gè)適當(dāng)?shù)拈撝?,達(dá)到濾除噪聲信號(hào)的目的。目前閾值的選取方式包括硬閾值法和軟閾值法兩種。其中硬閾值函數(shù)為

    (6)

    軟閾值函數(shù)

    (7)

    2 風(fēng)機(jī)增速箱故障信號(hào)重采樣的實(shí)現(xiàn)及采樣率的設(shè)置

    重采樣需要通過(guò)同步采集振動(dòng)信號(hào)與參考軸的轉(zhuǎn)速信號(hào)的方式實(shí)現(xiàn),并通過(guò)對(duì)鍵相信號(hào)的計(jì)算和插值,實(shí)現(xiàn)對(duì)振動(dòng)信號(hào)的等角度采樣。計(jì)算鍵相脈沖到達(dá)時(shí)刻通常有兩種方法:一種是在上升沿設(shè)定一個(gè)觸發(fā)閾值;一種是在下降沿設(shè)定一個(gè)觸發(fā)閾值。在定好觸發(fā)閾值之后,至少選取靠近閾值的一點(diǎn)作為鍵相時(shí)標(biāo)。其觸發(fā)閾值的取值范圍為

    (8)

    式中,fh為鍵相脈沖的高電平穩(wěn)態(tài)值;fi為鍵相脈沖的低電平穩(wěn)態(tài)值;fmax和fmin分別表示觸發(fā)閾值的最大值和最小值。在其范圍內(nèi)選取一個(gè)觸發(fā)閾值f,并設(shè)靠近觸發(fā)閾值前后的兩點(diǎn)分別記為fj和fj+1,與之對(duì)應(yīng)的時(shí)間標(biāo)記為tj和tj+1,采樣的周期為T,則可利用兩點(diǎn)線性插值法計(jì)算鍵相脈沖觸發(fā)閾值處的時(shí)標(biāo)ti為

    (9)

    由式(9)可得出觸發(fā)閾值處的時(shí)標(biāo),由此可求得每段脈沖間隔的平均轉(zhuǎn)速[11]

    (10)

    式中,k為每轉(zhuǎn)的脈沖數(shù)。

    由于重采樣過(guò)程中角域信號(hào)受時(shí)域采樣頻率性質(zhì)的影響以及采樣后仍然需要經(jīng)過(guò)FFT變換得到階次譜,為了避免階次混疊,在滿足奈奎斯特采樣定律的同時(shí),角域采樣的階次應(yīng)滿足:

    σs≥2 omax

    (11)

    式中,omax表示需要分析的最大階次。即采樣的角度間隔為

    (12)

    在對(duì)角域進(jìn)行采樣時(shí),角域的采樣率隨轉(zhuǎn)速的變化而變化,設(shè)某一時(shí)域信號(hào)X(ω,t)的采樣率為fs,則依據(jù)采樣定律,可得抗混濾波器的最大頻率:

    (13)

    在角域可將此信號(hào)看作是隨時(shí)域ω變化的變采樣率采樣信號(hào),即

    (14)

    由式(11)~(13)可得與角域?qū)?yīng)的時(shí)域抗混濾波最大階次

    (15)

    (16)

    綜上所述,在時(shí)域內(nèi)若想保留信息,依據(jù)式(16)及采樣定理可得fs滿足:

    (17)

    在角域采樣時(shí),若不進(jìn)行抗混濾波,由式(10)、(15)可得

    (18)

    由式(17)、(18)可知重采樣率是由時(shí)域采樣頻率和最小轉(zhuǎn)速共同決定的,所以為了避免混疊現(xiàn)象,最終采樣頻率要比預(yù)設(shè)值大。

    3 診斷實(shí)例分析

    為了驗(yàn)證上述方法的優(yōu)越性和可行性,對(duì)瓦房店某風(fēng)場(chǎng)存在潛在故障的1.5MW風(fēng)機(jī)(CCWE-1500/70.DF)增速箱進(jìn)行了長(zhǎng)時(shí)間的數(shù)據(jù)采集,在其增速箱的低速端和高速端安裝加速度傳感器,并采用Matlab軟件對(duì)所測(cè)到的故障信號(hào)進(jìn)行處理對(duì)比分析。其風(fēng)機(jī)增速箱的結(jié)構(gòu)簡(jiǎn)圖如圖1所示。

    圖1 增速箱的結(jié)構(gòu)簡(jiǎn)圖Fig.1 Structural diagram of gearbox

    由圖1可知,此型號(hào)的增速箱由一級(jí)行星輪系加二級(jí)平行軸組成。此次采集的振動(dòng)信號(hào)以輸入軸為參考軸,采樣頻率為2 048 Hz,采樣時(shí)間5 s,依據(jù)主軸轉(zhuǎn)速和各部件的嚙合情況[12],可計(jì)算出二級(jí)齒輪額定工作狀態(tài)下的1倍頻、2倍頻、3倍頻分別為97.39 Hz、194.78 Hz、389.56 Hz,轉(zhuǎn)動(dòng)頻率為9.8 Hz。此處齒輪點(diǎn)蝕的明顯故障特征就是齒輪齒面磨損嚴(yán)重時(shí)會(huì)出現(xiàn)高階振動(dòng),且隨轉(zhuǎn)軸頻率的增大,振動(dòng)能量也會(huì)明顯增大。在增速箱的輸出軸徑向安裝加速度傳感器,其采集的故障振動(dòng)信號(hào)的時(shí)域波形如圖2所示。

    圖2 故障齒輪的原始信號(hào)Fig.2 Original signal of failure gear

    由圖2的時(shí)域信號(hào)可以看出,振動(dòng)沖擊比較明顯,且信號(hào)中含有大量的噪聲信號(hào),由于轉(zhuǎn)速的不平穩(wěn)造成齒輪振動(dòng)信號(hào)具有非平穩(wěn)性,沖擊的時(shí)間間隔也不均勻。對(duì)其先進(jìn)行FFT分析,得到頻譜圖,如圖3所示。

    圖3 故障齒輪振動(dòng)信號(hào)的FFT頻譜Fig.3 FFT spectrum of fault gear vibration signal

    圖3中該信號(hào)在275 Hz處波動(dòng),F(xiàn)FT譜是非常模糊的,無(wú)法有效判斷其故障特征。本文分別利用經(jīng)典計(jì)算階次跟蹤(COT)、階次分析和小波閾值去噪相結(jié)合的方法對(duì)該故障信號(hào)進(jìn)行處理分析。對(duì)原始信號(hào)進(jìn)行sinc插值,得到等角度采樣信號(hào)如圖4所示,但由于得到的等角度采樣信號(hào)同樣會(huì)受噪聲和調(diào)制的雙重影響,使得階次分析產(chǎn)生的分析效果不是太明顯,所以對(duì)含有噪聲信號(hào)的重采樣信號(hào)選用db4小波對(duì)其進(jìn)行分解,分解層數(shù)為5層,采用硬閾值方法去噪,所得的等角度采樣信號(hào)如圖5所示,可見(jiàn)去噪后的信號(hào)逼近程度高,無(wú)相位失真和信號(hào)損失,頻譜更加清晰。

    圖4 等角度采樣信號(hào)Fig.4 Angle sampling signal

    圖5 硬閾值去噪后的等角度采樣信號(hào)Fig.5 The Angle sampling signal after the hard threshold de-noising

    得到等角度采樣信號(hào)后,對(duì)其分別進(jìn)行FFT變換,得到了故障信號(hào)的階次譜,如圖6、7所示。

    圖6 基于COT重采樣的階次譜Fig.6 Order spectrum of re-sampling based on the COT

    圖7 硬閾值去噪后基于COT重采樣的階次譜Fig.7 Order spectrum of re-sampling based on the COT after the hard threshold de-noising

    圖8 故障信號(hào)階次細(xì)化譜Fig.8 Zoom spectrum of fault signal order

    由圖4~7可以看出等角度重采樣的信號(hào)比較穩(wěn)定,得到的階次譜與FFT譜相比較,消除了頻譜的模糊現(xiàn)象。由圖7可知,信號(hào)在10.84階次處有明顯的波峰,邊頻帶少且稀疏,對(duì)應(yīng)嚙合頻率的1倍頻,故依據(jù)此處齒輪的故障特征初步定為磨損或點(diǎn)蝕。為進(jìn)一步明確故障類型,將5~20階次處的邊頻帶局部細(xì)化,所得細(xì)化譜如圖8所示,由細(xì)化譜可知,在嚙合階次10.84兩端的變頻較小,幅值的下降趨勢(shì)比較明顯,再結(jié)合時(shí)域波形的沖擊振動(dòng)比較明顯來(lái)看,此故障應(yīng)為點(diǎn)蝕故障。在本次的檢修過(guò)程中,對(duì)風(fēng)機(jī)齒輪箱進(jìn)行了拆卸,發(fā)現(xiàn)風(fēng)機(jī)齒輪箱內(nèi)部太陽(yáng)輪軸的齒輪齒面有金屬塊脫落的現(xiàn)象。此外為了驗(yàn)證階次分析和小波閾值去噪相結(jié)合的方法與經(jīng)典階次跟蹤方法(COT)重采樣精度的差別,利用均方根誤差來(lái)評(píng)價(jià)重采樣精度。

    (19)

    由式(19)可以求出基于經(jīng)典階次跟蹤方法重采樣的均方根誤差為0.9526,,基于小波閾值去噪和階次分析相結(jié)合的方法重采樣的均方根誤差為0.3728,由此可知對(duì)信號(hào)進(jìn)行經(jīng)典階次跟蹤的重采樣,再對(duì)其重采樣后的信號(hào)進(jìn)行小波閾值去噪,所得的重采樣精度更高。

    4 結(jié)論

    大型風(fēng)力機(jī)齒輪箱在實(shí)際運(yùn)行過(guò)程中,由于受環(huán)境和自身結(jié)構(gòu)的雙重影響,對(duì)其所采集到的時(shí)域信號(hào)往往具有非平穩(wěn)特性且含有大量的噪聲信號(hào),為了能更好的改善分析效果,提取故障信號(hào)的有效特征頻率,本文將階次分析和小波閾值去噪的診斷方法相結(jié)合,解決了因轉(zhuǎn)速波動(dòng)引起的頻率模糊和階次分析不適合在強(qiáng)噪聲環(huán)境條件下工作的問(wèn)題,并提高了重采樣的精度,通過(guò)實(shí)踐工程中的應(yīng)用和理論上的對(duì)比分析,證明了此方法在處理旋轉(zhuǎn)機(jī)械故障中的有效性,且能快速準(zhǔn)確的診斷出風(fēng)機(jī)齒輪箱的故障類型,在實(shí)際的故障診斷中具有廣闊的工程前景和意義。

    [1]JacekUrbanek,TomaszBarszcz,NaderSawalhi,etal.Comparisonofamplitude-basedandphase-basedmethodsforspeedtrackinginapplicationtowindturbines[J].MetrologyAndMeasurementSystems, 2011,XVIII(2): 295-304.

    [2]BartelmusW. ,ZimrozR.Anewfeatureformonitoringtheconditionofgearboxesinnon-stationaryoperatingconditions[J].MechanicalSystemsandSignalProcessing, 2009, 23(5):1528-1534.

    [3] 張磊,陳長(zhǎng)征,劉杰.自適應(yīng)EEMD法提取大型風(fēng)機(jī)主軸承早期故障特征[J].重型機(jī)械,2015(3):5-8.

    [4] 彭富強(qiáng),于德介,武春燕.基于自適應(yīng)時(shí)變?yōu)V波階比跟蹤的齒輪箱故障診斷[J].機(jī)械工程學(xué)報(bào),2012,48(7):77-85.

    [5] 陳向民,于德介,羅潔思.基于線調(diào)頻小波路徑追蹤階比循環(huán)平穩(wěn)解調(diào)的齒輪故障診斷[J].機(jī)械工程學(xué)報(bào),2012,48(3):95-101.

    [6] 程軍圣,李寶慶,楊宇.基于廣義解調(diào)時(shí)頻分析和瞬時(shí)頻率計(jì)算的階次譜方法在齒輪故障診斷中的應(yīng)用[J].振動(dòng)與沖擊,2011,30(9):30-34.

    [7] 叢華,吳廣平,繞國(guó)強(qiáng),等.計(jì)算階次分析中避免階次混疊的濾波定階方法及其應(yīng)用[J].振動(dòng)與沖擊,2012,31(12):42-45.

    [8] 宋寶玉,解志杰,張鋒,等.基于角度同步平均和階次分析的低速斜齒輪故障診斷[J].吉林大學(xué)學(xué)報(bào),2015,45(2):454-459.

    [9] 岳曉峰,朱成偉.階次分析在變速箱故障診斷中的應(yīng)用及研究[J].自動(dòng)化儀表,2015,36(6):8-11.

    [10]景新幸,冼燦嬌,楊海燕.基于改進(jìn)小波閾值去噪算法的確定[J].電聲技術(shù),2015,39(5):80-83.

    [11]李蓉,于德介,陳向民,等.基于階次分析與循環(huán)平穩(wěn)解調(diào)的齒輪箱復(fù)合故障診斷方法[J].中國(guó)機(jī)械工程,2013,24(10):1320-1326.

    [12]Amirat Y, Choqueuse V, Benbouzid M.EEMD-based wind turbine bearing failure detection using the generator stator current homopolar component[J] Mechanical Systems and Signal Processing,2013,41(1):667-678.

    Fault diagnosis of wind turbine gearbox based on order analysis

    CHEN Chang-zheng, WANG Hai-tong, SUN Zi-qiang

    (School of Mechanical Engineering, Shenyang University of Technology, Shenyang 110178, China)

    Wind turbine run under the poor working conditions, vibration signals often have the character of non-stationary. In order to extract the characteristic frequency of fault information, a new diagnosis algorithm is proposed based on the re-sampling order analysis and wavelet threshold de-noising. Compared with the algorithm of classical calculation order tracking (COT), the accuracy of re-sampling is improved, and compared with the methods of traditional FFT to eliminate the phenomenon of “frequency aliasing” cased by speed instability. Then through the wind field experiment shows that the method is effective and valuable in engineering.

    wind turbines; wavelet de-noising; order analysis; re-sampling

    2015-08-21;

    2015-09-25

    遼寧省自然科學(xué)基金(2014028017)

    陳長(zhǎng)征(1964),男,遼寧省沈陽(yáng)市,教授,博士生導(dǎo)師,主要從事振動(dòng)噪聲的研究。

    王海童(1988-),男,沈陽(yáng)工業(yè)大學(xué)碩士研究生。

    TH132.46

    A

    1001-196X(2016)02-0086-05

    猜你喜歡
    角域階次時(shí)域
    機(jī)載巡航導(dǎo)彈外形隱身改進(jìn)的電磁散射影響
    階次分析在驅(qū)動(dòng)橋異響中的應(yīng)用
    基于Vold-Kalman濾波的階次分析系統(tǒng)設(shè)計(jì)與實(shí)現(xiàn)*
    高空高速偵察機(jī)電磁散射特性研究
    基于時(shí)域信號(hào)的三電平逆變器復(fù)合故障診斷
    一類線性微分方程解的增長(zhǎng)性及Borel方向
    基于齒輪階次密度優(yōu)化的變速器降噪研究
    基于極大似然準(zhǔn)則與滾動(dòng)時(shí)域估計(jì)的自適應(yīng)UKF算法
    基于時(shí)域逆濾波的寬帶脈沖聲生成技術(shù)
    無(wú)限級(jí)整函數(shù)在角域內(nèi)的取值和增長(zhǎng)性
    国产成人精品福利久久| 亚洲精品国产av蜜桃| 成人18禁高潮啪啪吃奶动态图 | 精品国产一区二区三区久久久樱花| 亚洲,欧美,日韩| 欧美人与善性xxx| 国产在线视频一区二区| 夫妻性生交免费视频一级片| 你懂的网址亚洲精品在线观看| 制服丝袜香蕉在线| 亚洲国产精品一区二区三区在线| 亚洲电影在线观看av| 国产精品一区二区性色av| 久久热精品热| 日本av免费视频播放| 免费久久久久久久精品成人欧美视频 | 另类亚洲欧美激情| 国产精品国产三级专区第一集| 久久久国产一区二区| 免费观看的影片在线观看| 亚洲国产精品国产精品| 欧美日韩一区二区视频在线观看视频在线| 亚洲自偷自拍三级| 欧美日韩av久久| 在线观看人妻少妇| 日本黄大片高清| 色吧在线观看| 日韩免费高清中文字幕av| 高清毛片免费看| 十八禁网站网址无遮挡 | 久久精品国产亚洲av天美| 欧美激情极品国产一区二区三区 | 久久久久久久久久久久大奶| 欧美另类一区| 亚洲精品国产av成人精品| 国产成人a∨麻豆精品| 看非洲黑人一级黄片| 免费人妻精品一区二区三区视频| 大片电影免费在线观看免费| av在线播放精品| 日韩 亚洲 欧美在线| 国产亚洲av片在线观看秒播厂| 亚洲久久久国产精品| 我的老师免费观看完整版| 成人国产麻豆网| 亚洲av福利一区| 久久 成人 亚洲| 亚洲av综合色区一区| 自线自在国产av| 欧美97在线视频| 日日摸夜夜添夜夜添av毛片| 高清黄色对白视频在线免费看 | 色视频在线一区二区三区| 天堂俺去俺来也www色官网| 免费看不卡的av| 亚洲精品视频女| 国产亚洲午夜精品一区二区久久| 国产精品一二三区在线看| 国产精品不卡视频一区二区| av免费观看日本| 黄色怎么调成土黄色| 成人二区视频| 国产黄片美女视频| 日韩三级伦理在线观看| 曰老女人黄片| av专区在线播放| 在线观看av片永久免费下载| 国产综合精华液| 日韩一区二区三区影片| 伦理电影大哥的女人| 免费播放大片免费观看视频在线观看| 最后的刺客免费高清国语| av在线app专区| 在线天堂最新版资源| 国产精品国产三级国产av玫瑰| 人妻人人澡人人爽人人| 丝袜喷水一区| 亚洲精品乱码久久久久久按摩| 黑人猛操日本美女一级片| 最近的中文字幕免费完整| 在线播放无遮挡| 欧美精品一区二区免费开放| 国精品久久久久久国模美| 亚州av有码| 国产精品一区二区在线不卡| 亚洲精品国产av蜜桃| 黄色日韩在线| 2018国产大陆天天弄谢| 下体分泌物呈黄色| 午夜福利,免费看| 91aial.com中文字幕在线观看| 欧美国产精品一级二级三级 | 五月伊人婷婷丁香| 精品一区在线观看国产| 在线观看国产h片| 插逼视频在线观看| 寂寞人妻少妇视频99o| 蜜臀久久99精品久久宅男| 一本色道久久久久久精品综合| 春色校园在线视频观看| 久久亚洲国产成人精品v| 街头女战士在线观看网站| 久久人人爽av亚洲精品天堂| 亚洲精品一区蜜桃| 久久久久久久久久久久大奶| 校园人妻丝袜中文字幕| 自线自在国产av| 国产有黄有色有爽视频| 少妇的逼好多水| 最近的中文字幕免费完整| 午夜福利,免费看| av福利片在线| 99久久精品国产国产毛片| 美女中出高潮动态图| 亚洲性久久影院| 久久久a久久爽久久v久久| 成年人午夜在线观看视频| 韩国av在线不卡| 国产免费又黄又爽又色| 十八禁网站网址无遮挡 | 亚洲国产色片| 日韩,欧美,国产一区二区三区| 国内揄拍国产精品人妻在线| 欧美日韩av久久| 色5月婷婷丁香| 26uuu在线亚洲综合色| 国产成人免费无遮挡视频| 亚洲国产精品国产精品| 国产永久视频网站| 91精品国产国语对白视频| 国产日韩欧美亚洲二区| 男女啪啪激烈高潮av片| 久久韩国三级中文字幕| 两个人免费观看高清视频 | 亚洲精品一区蜜桃| 交换朋友夫妻互换小说| 久久97久久精品| 国产日韩欧美亚洲二区| 自拍偷自拍亚洲精品老妇| 日日摸夜夜添夜夜爱| 中文字幕人妻熟人妻熟丝袜美| 三级国产精品欧美在线观看| 日本爱情动作片www.在线观看| 十八禁高潮呻吟视频 | 超碰97精品在线观看| 伦理电影大哥的女人| 黑人高潮一二区| 99热6这里只有精品| 国产伦在线观看视频一区| 国产精品秋霞免费鲁丝片| av不卡在线播放| av视频免费观看在线观看| 欧美日韩在线观看h| 美女内射精品一级片tv| 少妇被粗大的猛进出69影院 | 99热这里只有是精品在线观看| 日韩一区二区视频免费看| 一级毛片黄色毛片免费观看视频| 偷拍熟女少妇极品色| 久久精品国产鲁丝片午夜精品| 日韩欧美 国产精品| 亚洲国产精品一区三区| 国产色爽女视频免费观看| 亚洲av国产av综合av卡| 亚洲成人手机| 人人妻人人看人人澡| 看非洲黑人一级黄片| 99久国产av精品国产电影| kizo精华| 草草在线视频免费看| 国产精品嫩草影院av在线观看| 97超碰精品成人国产| 看非洲黑人一级黄片| 国产精品熟女久久久久浪| 国产又色又爽无遮挡免| 十分钟在线观看高清视频www | 亚洲国产精品一区二区三区在线| 9色porny在线观看| 久久精品国产鲁丝片午夜精品| 91久久精品电影网| av国产精品久久久久影院| 成人美女网站在线观看视频| 只有这里有精品99| 日韩亚洲欧美综合| 伦理电影大哥的女人| 2018国产大陆天天弄谢| 高清毛片免费看| 欧美精品人与动牲交sv欧美| 好男人视频免费观看在线| 精品一品国产午夜福利视频| 国产成人91sexporn| 水蜜桃什么品种好| 国产av精品麻豆| 亚洲精品国产成人久久av| 熟妇人妻不卡中文字幕| 自拍偷自拍亚洲精品老妇| 亚洲国产最新在线播放| 精品人妻熟女毛片av久久网站| 亚洲av.av天堂| 狂野欧美白嫩少妇大欣赏| 亚洲中文av在线| 亚洲av福利一区| 国产视频内射| 人妻 亚洲 视频| 人妻夜夜爽99麻豆av| 久久久国产一区二区| 男人和女人高潮做爰伦理| av在线app专区| 国产精品久久久久久久电影| 婷婷色综合大香蕉| 成人漫画全彩无遮挡| 久久99热6这里只有精品| 国产毛片在线视频| 一区在线观看完整版| 夜夜骑夜夜射夜夜干| 少妇被粗大的猛进出69影院 | 黄色日韩在线| 久久99热这里只频精品6学生| 欧美一级a爱片免费观看看| 国产白丝娇喘喷水9色精品| 国产男人的电影天堂91| 蜜桃久久精品国产亚洲av| 久久久久久久久久久久大奶| 一本大道久久a久久精品| 欧美+日韩+精品| 国产欧美另类精品又又久久亚洲欧美| 熟妇人妻不卡中文字幕| 国产一区二区三区av在线| 国产在视频线精品| 高清视频免费观看一区二区| 国产精品免费大片| 又爽又黄a免费视频| 蜜臀久久99精品久久宅男| 亚洲精品国产av成人精品| 内射极品少妇av片p| 自拍欧美九色日韩亚洲蝌蚪91 | 精品午夜福利在线看| 嘟嘟电影网在线观看| 国产免费一级a男人的天堂| 国精品久久久久久国模美| 国产欧美日韩精品一区二区| 国产精品国产三级国产av玫瑰| 成人二区视频| 蜜臀久久99精品久久宅男| 久久国内精品自在自线图片| 看十八女毛片水多多多| 天堂8中文在线网| 成人综合一区亚洲| 亚洲欧美成人综合另类久久久| 男人爽女人下面视频在线观看| 国产伦精品一区二区三区四那| 亚洲一区二区三区欧美精品| 男女边吃奶边做爰视频| 51国产日韩欧美| 亚洲人成网站在线观看播放| 成人毛片60女人毛片免费| 又粗又硬又长又爽又黄的视频| 女性被躁到高潮视频| 在线观看人妻少妇| 国产成人精品无人区| 精品一区在线观看国产| 狂野欧美白嫩少妇大欣赏| 久久毛片免费看一区二区三区| 亚洲怡红院男人天堂| 秋霞在线观看毛片| 日韩一本色道免费dvd| 亚洲综合精品二区| 国产精品一区二区三区四区免费观看| 校园人妻丝袜中文字幕| 日韩一区二区视频免费看| 国产成人免费观看mmmm| 久久久久久久久久久丰满| 蜜臀久久99精品久久宅男| 老司机影院成人| 国产黄片美女视频| 啦啦啦视频在线资源免费观看| 午夜免费男女啪啪视频观看| 精品久久久久久久久av| av在线老鸭窝| 菩萨蛮人人尽说江南好唐韦庄| 另类精品久久| 亚洲精品久久久久久婷婷小说| 在现免费观看毛片| 欧美三级亚洲精品| 人妻人人澡人人爽人人| 久久国产精品大桥未久av | 日韩熟女老妇一区二区性免费视频| 国产精品一区二区在线不卡| 免费人成在线观看视频色| 亚洲国产av新网站| 丝瓜视频免费看黄片| 欧美日韩视频精品一区| 青春草国产在线视频| 国产免费一级a男人的天堂| videossex国产| 97超碰精品成人国产| 国产国拍精品亚洲av在线观看| 久久久久久久久久久丰满| av免费观看日本| 麻豆精品久久久久久蜜桃| 成人毛片60女人毛片免费| 一本久久精品| 在线天堂最新版资源| 国产成人a∨麻豆精品| 一个人看视频在线观看www免费| 有码 亚洲区| 涩涩av久久男人的天堂| 精品一区二区免费观看| 亚洲人成网站在线播| 91aial.com中文字幕在线观看| 欧美xxxx性猛交bbbb| 我的老师免费观看完整版| 一级毛片我不卡| 九草在线视频观看| 啦啦啦视频在线资源免费观看| 高清视频免费观看一区二区| 免费大片18禁| 人妻一区二区av| √禁漫天堂资源中文www| tube8黄色片| 日本wwww免费看| 亚洲国产欧美日韩在线播放 | 丁香六月天网| 国产一区二区在线观看av| 三上悠亚av全集在线观看 | 国产亚洲av片在线观看秒播厂| 亚洲欧美成人综合另类久久久| 欧美 日韩 精品 国产| 自拍偷自拍亚洲精品老妇| 久久久久国产精品人妻一区二区| 一级,二级,三级黄色视频| 亚洲欧美中文字幕日韩二区| 国产免费福利视频在线观看| 亚洲欧美清纯卡通| 亚洲欧美精品自产自拍| 成年人午夜在线观看视频| 久久人妻熟女aⅴ| 人妻人人澡人人爽人人| av.在线天堂| 精品卡一卡二卡四卡免费| 女性生殖器流出的白浆| 在线观看免费高清a一片| 日韩,欧美,国产一区二区三区| 天堂8中文在线网| 熟女av电影| 黑人猛操日本美女一级片| 国产淫语在线视频| 日韩av免费高清视频| 精品一区二区三卡| √禁漫天堂资源中文www| 日韩精品免费视频一区二区三区 | 成人毛片a级毛片在线播放| av不卡在线播放| 天堂中文最新版在线下载| 一级毛片我不卡| 国产av精品麻豆| 成人无遮挡网站| 色94色欧美一区二区| 91精品国产国语对白视频| 国产成人精品福利久久| 高清不卡的av网站| 大香蕉久久网| 国产精品久久久久久av不卡| 最近中文字幕2019免费版| av国产久精品久网站免费入址| .国产精品久久| 一区二区三区免费毛片| 人体艺术视频欧美日本| 2021少妇久久久久久久久久久| 肉色欧美久久久久久久蜜桃| 欧美 亚洲 国产 日韩一| 亚洲第一区二区三区不卡| av不卡在线播放| 成年美女黄网站色视频大全免费 | 成人综合一区亚洲| 成人特级av手机在线观看| 国产午夜精品一二区理论片| 国产精品一区二区在线不卡| 久久精品国产亚洲av涩爱| 在线观看美女被高潮喷水网站| 亚洲欧美成人精品一区二区| 亚洲av电影在线观看一区二区三区| 熟女av电影| 免费大片18禁| 女人久久www免费人成看片| 毛片一级片免费看久久久久| 黑丝袜美女国产一区| 国语对白做爰xxxⅹ性视频网站| 久久国产精品大桥未久av | 亚洲精品久久久久久婷婷小说| 三级国产精品片| av播播在线观看一区| 国产无遮挡羞羞视频在线观看| 嫩草影院入口| 在线天堂最新版资源| 国内精品宾馆在线| 久久热精品热| 啦啦啦啦在线视频资源| 久久久久精品久久久久真实原创| videos熟女内射| 亚洲精品乱久久久久久| 日韩欧美一区视频在线观看 | 国产精品一区二区三区四区免费观看| 久久青草综合色| 国产一区二区三区综合在线观看 | 国产综合精华液| 国产在线视频一区二区| 91久久精品国产一区二区成人| 国产精品不卡视频一区二区| 亚洲精品日韩av片在线观看| 亚洲美女黄色视频免费看| 国模一区二区三区四区视频| 在线精品无人区一区二区三| 亚洲欧美成人综合另类久久久| 黄色视频在线播放观看不卡| 亚洲成色77777| 精品少妇久久久久久888优播| 欧美日本中文国产一区发布| 能在线免费看毛片的网站| 高清午夜精品一区二区三区| av女优亚洲男人天堂| 一级,二级,三级黄色视频| 韩国高清视频一区二区三区| 国产精品无大码| 女的被弄到高潮叫床怎么办| 一级毛片 在线播放| 十八禁高潮呻吟视频 | 免费黄频网站在线观看国产| 亚洲欧洲日产国产| 99国产精品免费福利视频| 日日摸夜夜添夜夜爱| 国产av精品麻豆| 国产精品嫩草影院av在线观看| 日韩一本色道免费dvd| 欧美丝袜亚洲另类| 国内揄拍国产精品人妻在线| 成人美女网站在线观看视频| 日韩强制内射视频| av有码第一页| 2021少妇久久久久久久久久久| 在线观看av片永久免费下载| 丝袜在线中文字幕| 一级毛片 在线播放| 国产淫语在线视频| 少妇裸体淫交视频免费看高清| 三级经典国产精品| 国产片特级美女逼逼视频| 一级,二级,三级黄色视频| 国产无遮挡羞羞视频在线观看| 国产极品粉嫩免费观看在线 | 另类亚洲欧美激情| 久久6这里有精品| 国产av码专区亚洲av| 老司机影院成人| 2022亚洲国产成人精品| 国产无遮挡羞羞视频在线观看| 亚洲美女视频黄频| 亚洲三级黄色毛片| 综合色丁香网| av福利片在线观看| 国产精品99久久99久久久不卡 | av在线播放精品| 在线观看免费视频网站a站| 黑人猛操日本美女一级片| 欧美激情国产日韩精品一区| 能在线免费看毛片的网站| 中文天堂在线官网| 亚洲av综合色区一区| 精品久久久久久电影网| 夜夜爽夜夜爽视频| 亚洲情色 制服丝袜| 亚洲中文av在线| 国内精品宾馆在线| av播播在线观看一区| 丝袜在线中文字幕| 亚洲精品日韩av片在线观看| 欧美丝袜亚洲另类| 黄片无遮挡物在线观看| 亚洲精品中文字幕在线视频 | 丰满少妇做爰视频| 国产av精品麻豆| 亚洲va在线va天堂va国产| av.在线天堂| 热99国产精品久久久久久7| 成年美女黄网站色视频大全免费 | 男人狂女人下面高潮的视频| 亚洲真实伦在线观看| 建设人人有责人人尽责人人享有的| 美女大奶头黄色视频| 久久久久网色| 寂寞人妻少妇视频99o| 免费观看的影片在线观看| 成年美女黄网站色视频大全免费 | 国产精品蜜桃在线观看| av免费观看日本| 久久久久久久国产电影| 亚洲av成人精品一二三区| 久久精品国产亚洲av天美| 久久久久久久久大av| 亚洲国产精品国产精品| 国产一区二区在线观看av| av线在线观看网站| 亚洲国产精品专区欧美| 三上悠亚av全集在线观看 | 亚洲国产成人一精品久久久| 亚洲av.av天堂| 国产午夜精品一二区理论片| 精品少妇内射三级| 久久免费观看电影| 国产爽快片一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 最近手机中文字幕大全| 国产精品国产三级专区第一集| 黑人巨大精品欧美一区二区蜜桃 | 亚洲av男天堂| 免费看不卡的av| 毛片一级片免费看久久久久| a级毛片在线看网站| 99热这里只有精品一区| 日本vs欧美在线观看视频 | 精品一区二区免费观看| 中文在线观看免费www的网站| 十分钟在线观看高清视频www | 色94色欧美一区二区| 91精品国产国语对白视频| 色网站视频免费| 狂野欧美激情性bbbbbb| 国产伦精品一区二区三区视频9| 99久久精品一区二区三区| 中文字幕av电影在线播放| 边亲边吃奶的免费视频| 精品99又大又爽又粗少妇毛片| 国产成人精品婷婷| 国产在线免费精品| 麻豆成人午夜福利视频| 日韩精品有码人妻一区| 99久久精品一区二区三区| 少妇猛男粗大的猛烈进出视频| 男的添女的下面高潮视频| 狠狠精品人妻久久久久久综合| 午夜久久久在线观看| 成人免费观看视频高清| 97精品久久久久久久久久精品| 亚洲美女视频黄频| kizo精华| 中文字幕精品免费在线观看视频 | 日韩av不卡免费在线播放| 黑人巨大精品欧美一区二区蜜桃 | 一级片'在线观看视频| 免费看光身美女| 男人和女人高潮做爰伦理| 精华霜和精华液先用哪个| 99热这里只有是精品50| 精品少妇久久久久久888优播| 亚洲av日韩在线播放| 中文字幕制服av| 午夜久久久在线观看| 色94色欧美一区二区| 少妇熟女欧美另类| 一级毛片 在线播放| 久久久久久久久久成人| 成年人午夜在线观看视频| 免费播放大片免费观看视频在线观看| 欧美变态另类bdsm刘玥| 午夜91福利影院| av播播在线观看一区| 久久国产精品男人的天堂亚洲 | 国产精品一区二区在线不卡| 全区人妻精品视频| 亚洲精品中文字幕在线视频 | 老女人水多毛片| 99九九线精品视频在线观看视频| 欧美日韩一区二区视频在线观看视频在线| 亚洲av综合色区一区| 久久久国产一区二区| h日本视频在线播放| 少妇人妻精品综合一区二区| 99久久中文字幕三级久久日本| 麻豆成人av视频| 麻豆成人午夜福利视频| 亚洲怡红院男人天堂| 国语对白做爰xxxⅹ性视频网站| 国产午夜精品久久久久久一区二区三区| 国产免费又黄又爽又色| 久久韩国三级中文字幕| 欧美亚洲 丝袜 人妻 在线| 一级爰片在线观看| 卡戴珊不雅视频在线播放| 亚洲精品亚洲一区二区| 国产色爽女视频免费观看| 亚洲色图综合在线观看| 亚洲精品自拍成人| 久久久久久久亚洲中文字幕| 一区二区三区四区激情视频| 亚洲一区二区三区欧美精品| 又黄又爽又刺激的免费视频.| 精品少妇久久久久久888优播| 欧美精品人与动牲交sv欧美| av在线app专区| 亚洲精品久久午夜乱码| 亚洲精品aⅴ在线观看| 伊人久久国产一区二区| 天美传媒精品一区二区| 日韩一本色道免费dvd| 亚洲美女搞黄在线观看| av天堂久久9| 大话2 男鬼变身卡| 亚洲精品456在线播放app| 午夜久久久在线观看| 熟女电影av网| 香蕉精品网在线| 亚洲欧洲精品一区二区精品久久久 | 成人午夜精彩视频在线观看| 国模一区二区三区四区视频| 一区在线观看完整版|