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

    基于集合經(jīng)驗(yàn)?zāi)B(tài)分解與希爾伯特-黃變換的有軌電車(chē)異常振動(dòng)分析

    2022-11-15 06:59:34閆轉(zhuǎn)芳張會(huì)杰季元進(jìn)
    科技創(chuàng)新與應(yīng)用 2022年32期
    關(guān)鍵詞:軸箱區(qū)段幅值

    閆轉(zhuǎn)芳,張會(huì)杰,季元進(jìn)

    (1.中車(chē)青島四方機(jī)車(chē)車(chē)輛股份有限公司,山東 青島 266000;2.同濟(jì)大學(xué),上海 200000)

    有軌電車(chē)在運(yùn)行時(shí),若車(chē)輛發(fā)生故障,一些故障信息會(huì)在車(chē)輛受到線路沖擊時(shí),通過(guò)異常的振動(dòng)信號(hào)體現(xiàn)出來(lái)。異常振動(dòng)使振動(dòng)信號(hào)在時(shí)頻特性上發(fā)生相應(yīng)的變化,異常振動(dòng)信號(hào)與正常信號(hào)的頻率、能量分布以及能量大小都有所不同。對(duì)有軌電車(chē)振動(dòng)信號(hào)的分析有諸多時(shí)域及頻域方法,然而以傅里葉變換為代表的傳統(tǒng)時(shí)頻分析方法在處理復(fù)雜的故障信號(hào)時(shí)有一定的局限性[1-2]。有軌電車(chē)異常振動(dòng)信號(hào)多為非平穩(wěn)非線性數(shù)據(jù),在對(duì)這類(lèi)信號(hào)的處理上,希爾伯特-黃變換(HHT)具有明顯的優(yōu)勢(shì)。HHT的分解依據(jù)為信號(hào)本身的時(shí)間尺度特性,與小波分解和傅里葉分解等方法的本質(zhì)區(qū)別在于HHT免去了預(yù)先設(shè)定基函數(shù)這一步[3]。

    由于信號(hào)的內(nèi)部特征可以由HHT準(zhǔn)確地展現(xiàn)出來(lái),且HHT的時(shí)頻域分辨率相對(duì)較高,因此HHT在振動(dòng)分析、故障檢測(cè)、參數(shù)識(shí)別等領(lǐng)域具有重大意義,被廣泛地應(yīng)用。陳雙喜等[4]提出利用一種改進(jìn)的經(jīng)驗(yàn)?zāi)B(tài)分解方法,有效地提取車(chē)輛-軌道耦合系統(tǒng)的動(dòng)力學(xué)特性。李再幃等[5]利用HHT方法對(duì)車(chē)輛-軌道系統(tǒng)垂向振動(dòng)進(jìn)行時(shí)頻分析,對(duì)軌道高低不平順與車(chē)輛垂向振動(dòng)加速度的關(guān)系進(jìn)行定量研究。蘇鵬等[6]將HHT應(yīng)用于橋梁的振動(dòng)分析中,對(duì)列車(chē)過(guò)橋時(shí)的振動(dòng)響應(yīng)進(jìn)行處理,為橋梁振動(dòng)分析與健康檢測(cè)提供了有力的工具。

    本文以某型有軌電車(chē)為研究對(duì)象,提出了基于EEMD和HHT的有軌電車(chē)異常振動(dòng)分析研究方法。模態(tài)混疊問(wèn)題難以避免地發(fā)生在傳統(tǒng)HHT采用的經(jīng)驗(yàn)?zāi)B(tài)分解當(dāng)中,針對(duì)此問(wèn)題,本文提出利用集合經(jīng)驗(yàn)?zāi)B(tài)分解(EEMD)方法對(duì)原始信號(hào)進(jìn)行分解,并利用相關(guān)系數(shù)法選取真實(shí)分量,剔除無(wú)意義的虛假分量。對(duì)真實(shí)分量進(jìn)行Hilbert變換得到異常振動(dòng)信號(hào)的Hilbert譜,與正常振動(dòng)狀態(tài)對(duì)比分析,實(shí)現(xiàn)有軌電車(chē)異常振動(dòng)與正常振動(dòng)狀態(tài)的時(shí)頻特征區(qū)分。對(duì)Hilbert積分得振動(dòng)信號(hào)的邊際譜,通過(guò)邊際譜得到異常振動(dòng)的主頻,分析異常振動(dòng)區(qū)段對(duì)應(yīng)的線路問(wèn)題及故障模式。

    1 基本原理

    1.1 鐵路接頭平順性評(píng)判方法

    鋼軌接頭引起的瞬時(shí)脈沖經(jīng)過(guò)傅里葉變換后,在頻域上是一個(gè)無(wú)限帶寬的信號(hào),在理論上應(yīng)當(dāng)存在所有頻率的成分,并且各頻率成分的能量沿整個(gè)頻率軸均布。因此該信號(hào)成分中必定含有輪對(duì)-軸箱系統(tǒng)的固有頻率,也就必定會(huì)引起系統(tǒng)發(fā)生廣義共振?;谏鲜鲈?,共振解調(diào)被用作鐵路接頭平順性的評(píng)判方法。

    鐵路接頭平順性評(píng)判方法可描述如下。

    (1)采集軸箱加速度。

    (2)對(duì)軸箱垂向振動(dòng)加速度進(jìn)行帶通濾波,帶寬為[20,450]Hz。

    (3)計(jì)算濾波后的軸箱垂向振動(dòng)加速度有效值。(4)將有效值進(jìn)行歸一化操作得到接縫不良指數(shù)如下

    式中:mRMS是有效值的平均值;n是RMS點(diǎn)數(shù)。

    (5)若接縫不良指數(shù)WJI大于4,也就是該處的沖擊能量是平均沖擊能量的4倍時(shí),則判定接頭平順性不良。

    1.2 集合經(jīng)驗(yàn)?zāi)B(tài)分解

    集合經(jīng)驗(yàn)?zāi)B(tài)分解(EEMD)數(shù)據(jù)分析方法運(yùn)用噪聲輔助,對(duì)信號(hào)數(shù)據(jù)疊加高斯白噪聲,進(jìn)行多次經(jīng)驗(yàn)?zāi)J椒纸?,有效解決了經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)的模態(tài)混疊問(wèn)題[7]。EMD原理如下:通過(guò)對(duì)原始信號(hào)x(t)的上、下包絡(luò)線計(jì)算出平均包絡(luò)線m1(t)。用原始信號(hào)x(t)減去m1(t)得到一個(gè)去掉低頻的新數(shù)據(jù)序列h1(t),通過(guò)h1(t)的極值點(diǎn)個(gè)數(shù)與過(guò)零點(diǎn)個(gè)數(shù)相差不超過(guò)1及h1(t)的上、下包絡(luò)線均值恒為0兩點(diǎn)依據(jù)判定h1(t)是否為一個(gè)本征模分量(IMF)。若不是則重復(fù)操作直至得到滿(mǎn)足條件的本征模分量,記作c1(t)。c1(t)表示原始信號(hào)中的最高頻成分。用原始信號(hào)x(t)減去c1(t),得新序列r1(t):將r1(t)作為新的“原始信號(hào)”重復(fù)以上步驟,直至殘余信號(hào)分量rn為周期函數(shù)、單調(diào)函數(shù)、常數(shù)或僅含一個(gè)極值的簡(jiǎn)單信號(hào)函數(shù)。原始信號(hào)經(jīng)EMD分解為

    式中:ci為所得到的各個(gè)IMF分量,包含了信號(hào)高頻段至低頻段的不同成分;rn(t)為殘余信號(hào)。

    EEMD是指在待分解的信號(hào)中加入高斯白噪聲,將融合白噪聲的信號(hào)通過(guò)EMD分解為多個(gè)IMF。由于高斯白噪聲擁有頻率均勻分布和幅值零均值的特點(diǎn),因此運(yùn)用不相同但幅值相等的高斯白噪聲,將其加入信號(hào)當(dāng)中,以此來(lái)改變信號(hào)的極值點(diǎn)特性[8]。由于噪聲經(jīng)過(guò)總體平均后被互相抵消,最終結(jié)果來(lái)自于多次EMD分解以后得到的IMF集成均值,以防止模態(tài)混疊的發(fā)生。

    1.3 相關(guān)性檢驗(yàn)

    EEMD分解得到的所有IMF并非都是有意義的,其中可能包含部分虛假、無(wú)意義的分量,如噪聲分量和趨勢(shì)項(xiàng)分量等。各IMF與原信號(hào)的相關(guān)程度評(píng)估是基于相關(guān)系數(shù)法,相關(guān)系數(shù)較低的被稱(chēng)為虛假分量,需將其除去。保留真實(shí)IMF進(jìn)行重構(gòu),達(dá)到降噪的目的[9]。信號(hào)中存在一些幅值較小的IMF,其真實(shí)存在,但容易被當(dāng)作虛假分量而被除去。為了防止誤刪的情況發(fā)生,在IMF與原信號(hào)計(jì)算相關(guān)系數(shù)前,需要對(duì)它們進(jìn)行歸一化處理[10]。歸一化相關(guān)系數(shù)為

    1.4 希爾伯特-黃變換

    希爾伯特-黃變換(HHT)將非平穩(wěn)信號(hào)進(jìn)行平穩(wěn)化處理,也適用于非線性信號(hào)[8]。包含有不同尺度的趨勢(shì)或波動(dòng)的信號(hào)被逐級(jí)分解成為多個(gè)IMF,然后將各IMF進(jìn)行HHT得到時(shí)間與頻率聯(lián)合分布的Hilbert譜。原信號(hào)的頻率含量由瞬時(shí)頻率與瞬時(shí)能量?jī)芍笜?biāo)來(lái)表征。

    其步驟如下:首先,對(duì)每個(gè)本征模分量ci(t)作Hilb ert變換得到di(t)

    構(gòu)造解析信號(hào)

    得相應(yīng)IMF的瞬時(shí)幅值為

    相位為

    進(jìn)一步可以得瞬時(shí)頻率

    因此,原始信號(hào)x(t)可展成

    式中:Re表示取實(shí)部。式(9)表征了x(t)的時(shí)間-頻率-幅值的三維關(guān)系,將信號(hào)幅度在時(shí)頻上的分布稱(chēng)為Hilbert時(shí)頻譜,簡(jiǎn)稱(chēng)Hilbert譜,以H(w,t)表示,其數(shù)學(xué)表達(dá)式為

    對(duì)H(w,t)積分可得Hilbert邊際譜h(w)

    式中:T為信號(hào)的總長(zhǎng)度。一方面,邊際譜反映了信號(hào)瞬時(shí)頻率總幅值的大小,在數(shù)理統(tǒng)計(jì)角度來(lái)說(shuō),其表示的是全部所測(cè)數(shù)據(jù)的累加幅值;另一方面,其反映了在信號(hào)整個(gè)時(shí)間段內(nèi),每個(gè)頻率值對(duì)累加幅值的貢獻(xiàn)。

    H(w,t)是一種變化規(guī)律,用來(lái)代表所測(cè)數(shù)據(jù)的幅值隨時(shí)間和頻率,在整個(gè)頻段上的變化規(guī)律,而h(w)反映的是其變化情況。

    2 有軌電車(chē)運(yùn)行試驗(yàn)

    本文所進(jìn)行特征提取與異常振動(dòng)分析的有軌電車(chē)振動(dòng)數(shù)據(jù)來(lái)自某型有軌電車(chē),整列車(chē)包括3節(jié)車(chē)體與4個(gè)轉(zhuǎn)向架,轉(zhuǎn)向架全部為獨(dú)立輪對(duì)轉(zhuǎn)向架,其中車(chē)體連接處采用鉸接式獨(dú)立車(chē)輪轉(zhuǎn)向架。被測(cè)車(chē)輛的一、二系彈簧、減振器狀態(tài)均正常。車(chē)輪踏面選用LO67588P踏面。正線直線為50 kg工字軌,軌底坡1/40;部分線路為Ri60R2槽型軌及護(hù)軌。

    在有軌電車(chē)的三位軸軸箱上布置垂向加速度計(jì),壓電式,型號(hào)為L(zhǎng)A0107T-100,量程100 g,測(cè)試精度為0.5%。有軌電車(chē)編組及測(cè)點(diǎn)布置如圖1所示。采樣頻率1000 Hz。

    圖1 有軌電車(chē)編組及軸箱加速度測(cè)點(diǎn)布置

    有軌電車(chē)運(yùn)行試驗(yàn)工況為AW0空載,運(yùn)行速度為40 km/h。運(yùn)行途中某段異常振動(dòng)區(qū)間的軸箱垂向加速度數(shù)據(jù)如圖2所示。由加速度圖像可以看出該區(qū)段行駛時(shí)有一些明顯的沖擊位置,相應(yīng)地引起有軌電車(chē)的異常振動(dòng)。

    圖2 試驗(yàn)全程軸箱垂向加速度

    采用鐵路接頭平順性評(píng)判方法對(duì)振動(dòng)數(shù)據(jù)進(jìn)行初步研究。該區(qū)段接縫不良指數(shù)如圖3所示。由圖3可以看出,從1.7 s開(kāi)始,大約每過(guò)2.3 s出現(xiàn)一次沖擊,沖擊位置的接縫不良指數(shù)遠(yuǎn)遠(yuǎn)超過(guò)合格值4。該區(qū)段鐵軌采用有縫軌道,根據(jù)沖擊點(diǎn)間隔初步判斷沖擊是由軌縫處不平順引起的。異常振動(dòng)的特征頻率分析及原因分析見(jiàn)下文。

    圖3 接縫不良指數(shù)圖像

    3 異常振動(dòng)信號(hào)EEMD與HHT分析

    對(duì)異常振動(dòng)信號(hào)的分析流程如圖4所示。首先,對(duì)有軌電車(chē)軸箱的異常振動(dòng)信號(hào)進(jìn)行EEMD分解,得到各個(gè)IMF分量及其頻譜圖如圖5所示。IMF1-IMF8是從高頻到低頻的IMF分量,RES為殘余分量。

    圖4 算法流程

    圖5 異常振動(dòng)信號(hào)EEMD各分量時(shí)域圖與頻域圖

    采用相關(guān)系數(shù)法計(jì)算各IMF與原始信號(hào)之間的相關(guān)系數(shù),見(jiàn)表1,其中IMF1—IMF4的相關(guān)系數(shù)較高,并且是原始振動(dòng)信號(hào)的主要成分,集中了原始信號(hào)中最顯著的特征,對(duì)整體響應(yīng)的影響最大;而IMF5—IMF8的相關(guān)系數(shù)均小于0.1,可認(rèn)定為噪聲分量或趨勢(shì)量分量,對(duì)原信號(hào)的特征提取無(wú)意義,故剔除。

    表1 異常振動(dòng)信號(hào)各IMF與原始信號(hào)的相關(guān)系數(shù)

    原始信號(hào)與由IMF1、IMF2、IMF3、IMF4組成的重構(gòu)信號(hào)對(duì)比如圖6所示。由圖可見(jiàn),原信號(hào)與去噪后的信號(hào)基本重合,原始信號(hào)的主要特征得以保留。

    圖6 原始信號(hào)與重構(gòu)信號(hào)對(duì)比

    對(duì)軸箱垂向振動(dòng)信號(hào)的IMF1—IMF4分量進(jìn)行瞬時(shí)處理,得到瞬時(shí)頻率、瞬時(shí)相位及瞬時(shí)能量,進(jìn)一步得到振動(dòng)信號(hào)的Hilbert譜圖,如圖7所示。由圖7可知,該信號(hào)能量幾乎全部分散在0~150 Hz的寬頻帶中,頻率高于450 Hz的部分幾乎沒(méi)有能量分布。其中50~100 Hz頻帶內(nèi)的能量較大,其他頻段能量比較弱。在1.1、3.3、5.5、8.5、10.7、13.25 s前后在50~100 Hz頻帶能量出現(xiàn)集中放大現(xiàn)象,瞬時(shí)能量遠(yuǎn)大于其他時(shí)段的瞬時(shí)能量。振動(dòng)全程中存在4處明顯的沖擊。

    圖7振動(dòng)信號(hào)的Hilbert譜

    圖8 直觀地展示了振動(dòng)信號(hào)的時(shí)間-瞬時(shí)頻率-瞬時(shí)能量的三維關(guān)系,由圖8可知,該區(qū)段的振動(dòng)信號(hào)存在6處沖擊,且沖擊間隔平均在2.2 s左右,且沖擊能量逐次增強(qiáng),在5.5 s,70 Hz達(dá)到能量最高點(diǎn),最高瞬時(shí)能量為21.87。

    圖8 時(shí)間-瞬時(shí)頻率-瞬時(shí)能量圖像

    選取上述振動(dòng)信號(hào)的5.4~5.6 s的沖擊區(qū)段和5.8~6.0 s的平穩(wěn)區(qū)段進(jìn)行希爾伯特-黃變換,并通過(guò)HHT后的任意瞬時(shí)頻率對(duì)時(shí)間的積分,得到異沖擊區(qū)段和平穩(wěn)區(qū)段的邊際譜,用以表示信號(hào)各瞬時(shí)頻率的總幅值的大小,如圖9所示。由圖9可知,沖擊區(qū)段信號(hào)的能量主要分布在60 Hz附近的低頻區(qū)域,信號(hào)的主振幅為0.0704,高頻區(qū)域總的振動(dòng)量較小。平穩(wěn)區(qū)段信號(hào)的能量同樣是主要分布在60 Hz附近的低頻區(qū)域,信號(hào)的主振幅為0.0278。

    圖9 振動(dòng)信號(hào)的邊際譜

    4 軌道線路勘察

    該段線路采用25 m標(biāo)準(zhǔn)鋼軌,各段鋼軌之間的軌縫較寬。在該異常振動(dòng)區(qū)段的軌縫處,軌道連接板較高,導(dǎo)致有軌電車(chē)經(jīng)過(guò)時(shí)鋼輪與軌道連接板的頂面產(chǎn)生沖擊。該速度下,約2.25 s經(jīng)過(guò)一次軌縫并產(chǎn)生一次沖擊,造成短暫時(shí)段內(nèi)的異常振動(dòng)。由上節(jié)HHT結(jié)果可知,車(chē)輪與軌道連接板頂面的沖擊激發(fā)了有軌電車(chē)與軌道的共振,共振頻率為60 Hz,沖擊區(qū)段產(chǎn)生很大的瞬時(shí)能量。圖10顯示了沖擊處軌道連接板的磨損光帶。

    圖10 軌道連接板的頂面沖擊

    5 結(jié)論

    本文采用HHT進(jìn)行有軌電車(chē)異常振動(dòng)分析。得到了以下結(jié)論。

    (1)當(dāng)有軌電車(chē)發(fā)生異常振動(dòng)時(shí),其振動(dòng)信號(hào)的Hilbert譜與正常振動(dòng)是不同的,可以利用振動(dòng)信號(hào)的Hilbert譜來(lái)進(jìn)行有軌電車(chē)異常振動(dòng)分析及故障診斷。

    (2)邊際譜可以反映有軌電車(chē)振動(dòng)信號(hào)中各個(gè)頻率與總幅值的關(guān)系,明顯區(qū)分出異常振動(dòng)與正常振動(dòng),找到異常振動(dòng)的真實(shí)頻率值,通過(guò)異常頻率成分對(duì)有軌電車(chē)進(jìn)行故障診斷。

    猜你喜歡
    軸箱區(qū)段幅值
    中老鐵路雙線區(qū)段送電成功
    某型車(chē)輛軸箱在線振動(dòng)測(cè)試分析
    基于SPS模式的轉(zhuǎn)向架軸箱裝配線仿真研究
    站內(nèi)特殊區(qū)段電碼化設(shè)計(jì)
    站內(nèi)軌道區(qū)段最小長(zhǎng)度的探討
    雙牙棒消隙C軸箱的結(jié)構(gòu)設(shè)計(jì)
    基于S變換的交流電網(wǎng)幅值檢測(cè)系統(tǒng)計(jì)算機(jī)仿真研究
    電子制作(2017年7期)2017-06-05 09:36:13
    淺析分路不良區(qū)段解鎖的特殊操作
    正序電壓幅值檢測(cè)及諧波抑制的改進(jìn)
    209 P型轉(zhuǎn)向架軸箱彈簧斷裂的原因及預(yù)防
    免费女性裸体啪啪无遮挡网站| 日韩一本色道免费dvd| 免费av不卡在线播放| 99国产综合亚洲精品| 中文天堂在线官网| 男女免费视频国产| 永久免费av网站大全| 在线看a的网站| a级毛片在线看网站| 午夜日本视频在线| 男人操女人黄网站| 草草在线视频免费看| 亚洲国产av新网站| 少妇的丰满在线观看| 一级a做视频免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久视频综合| 多毛熟女@视频| 97在线视频观看| 国产成人aa在线观看| 午夜91福利影院| 一级毛片我不卡| 久久午夜福利片| 久久精品国产综合久久久 | 十八禁网站网址无遮挡| 少妇高潮的动态图| 亚洲经典国产精华液单| 美女脱内裤让男人舔精品视频| 人成视频在线观看免费观看| 亚洲av国产av综合av卡| 日产精品乱码卡一卡2卡三| 国产一区二区三区av在线| 国产熟女午夜一区二区三区| 国产成人免费观看mmmm| 在线观看三级黄色| 欧美少妇被猛烈插入视频| 国产片内射在线| 亚洲,一卡二卡三卡| 午夜福利影视在线免费观看| www日本在线高清视频| 国产亚洲精品久久久com| a 毛片基地| 久久久精品区二区三区| 国产一级毛片在线| 韩国精品一区二区三区 | 精品一区二区免费观看| 国产一区亚洲一区在线观看| 男女高潮啪啪啪动态图| 久久ye,这里只有精品| 看免费av毛片| 国产在视频线精品| 久久国产亚洲av麻豆专区| 欧美另类一区| 久久久久人妻精品一区果冻| 美女福利国产在线| 欧美人与性动交α欧美软件 | 少妇高潮的动态图| 免费观看在线日韩| freevideosex欧美| 免费观看a级毛片全部| 99热网站在线观看| 国产精品免费大片| xxxhd国产人妻xxx| 美女内射精品一级片tv| 国产精品99久久99久久久不卡 | 五月玫瑰六月丁香| 狂野欧美激情性xxxx在线观看| 高清不卡的av网站| 亚洲精品日本国产第一区| 亚洲五月色婷婷综合| 久久久精品区二区三区| 18禁观看日本| 国产又爽黄色视频| 亚洲成av片中文字幕在线观看 | 美女中出高潮动态图| 国产精品人妻久久久久久| 国产白丝娇喘喷水9色精品| 欧美xxⅹ黑人| 欧美日韩视频精品一区| 美女国产视频在线观看| 91久久精品国产一区二区三区| 久久国产亚洲av麻豆专区| 人妻少妇偷人精品九色| 久久精品夜色国产| 999精品在线视频| 欧美精品高潮呻吟av久久| 国产免费又黄又爽又色| 中国美白少妇内射xxxbb| 久久久久精品人妻al黑| 人体艺术视频欧美日本| 亚洲av中文av极速乱| 哪个播放器可以免费观看大片| 一级毛片 在线播放| 中国国产av一级| 精品少妇内射三级| 美女国产高潮福利片在线看| 妹子高潮喷水视频| 国产亚洲精品第一综合不卡 | 免费久久久久久久精品成人欧美视频 | 99久久中文字幕三级久久日本| 91久久精品国产一区二区三区| 久久 成人 亚洲| 午夜激情av网站| 老熟女久久久| 日韩欧美一区视频在线观看| 免费大片黄手机在线观看| 久久久久久久久久久免费av| 国产熟女欧美一区二区| 最新的欧美精品一区二区| 久久精品久久久久久噜噜老黄| √禁漫天堂资源中文www| 亚洲成人手机| 五月天丁香电影| 最近的中文字幕免费完整| 婷婷成人精品国产| 在线观看人妻少妇| 亚洲精品一二三| 国产白丝娇喘喷水9色精品| 伊人亚洲综合成人网| 老司机影院毛片| 亚洲经典国产精华液单| 亚洲第一区二区三区不卡| 丝袜人妻中文字幕| 久久99精品国语久久久| 最新中文字幕久久久久| 高清av免费在线| 人人妻人人添人人爽欧美一区卜| 咕卡用的链子| 熟女人妻精品中文字幕| 香蕉精品网在线| 亚洲三级黄色毛片| 99国产精品免费福利视频| 全区人妻精品视频| 国产精品成人在线| 日韩制服骚丝袜av| 国产激情久久老熟女| 综合色丁香网| 性色avwww在线观看| 老司机影院成人| 黄色 视频免费看| 少妇精品久久久久久久| 久久亚洲国产成人精品v| 国产成人精品在线电影| 欧美xxⅹ黑人| 制服诱惑二区| 中国三级夫妇交换| 少妇的丰满在线观看| 伦理电影免费视频| 午夜91福利影院| 日韩一区二区三区影片| 人人澡人人妻人| 五月玫瑰六月丁香| 久久热在线av| 亚洲av福利一区| 九色成人免费人妻av| 成人国语在线视频| 日本爱情动作片www.在线观看| 欧美人与性动交α欧美软件 | 少妇高潮的动态图| 母亲3免费完整高清在线观看 | 日韩,欧美,国产一区二区三区| 插逼视频在线观看| 免费人妻精品一区二区三区视频| 亚洲综合色网址| 亚洲,一卡二卡三卡| 曰老女人黄片| 天天躁夜夜躁狠狠躁躁| 97超碰精品成人国产| 国产精品久久久久久久电影| 日本av免费视频播放| av又黄又爽大尺度在线免费看| 黄色毛片三级朝国网站| 久久精品夜色国产| 熟妇人妻不卡中文字幕| 人妻 亚洲 视频| 国产一区二区激情短视频 | 欧美 日韩 精品 国产| 亚洲国产看品久久| 久久久久久久久久成人| av网站免费在线观看视频| 一区二区av电影网| 久久久久人妻精品一区果冻| 色婷婷久久久亚洲欧美| 亚洲美女搞黄在线观看| 成年动漫av网址| videosex国产| 日韩中字成人| 亚洲欧洲国产日韩| www.熟女人妻精品国产 | 久久av网站| 久久久久久久久久久免费av| 欧美日韩成人在线一区二区| 日韩中文字幕视频在线看片| 91在线精品国自产拍蜜月| 韩国av在线不卡| 久久精品aⅴ一区二区三区四区 | 激情视频va一区二区三区| 精品亚洲成a人片在线观看| 亚洲久久久国产精品| 久久久亚洲精品成人影院| 免费不卡的大黄色大毛片视频在线观看| 91精品伊人久久大香线蕉| 妹子高潮喷水视频| 成年女人在线观看亚洲视频| 国产精品一区二区在线观看99| 中文字幕人妻熟女乱码| 国产精品欧美亚洲77777| 免费日韩欧美在线观看| 国产欧美日韩一区二区三区在线| 秋霞伦理黄片| 欧美日韩一区二区视频在线观看视频在线| 制服诱惑二区| 国产片内射在线| 97人妻天天添夜夜摸| 国产精品麻豆人妻色哟哟久久| 国产一区二区激情短视频 | 免费少妇av软件| 黑丝袜美女国产一区| 黄片无遮挡物在线观看| 在线观看三级黄色| 99九九在线精品视频| 如何舔出高潮| 国产又色又爽无遮挡免| 人成视频在线观看免费观看| 国产精品.久久久| 如日韩欧美国产精品一区二区三区| 色哟哟·www| 又黄又粗又硬又大视频| 精品人妻一区二区三区麻豆| 在线天堂中文资源库| 91国产中文字幕| 男女高潮啪啪啪动态图| 老司机亚洲免费影院| 十八禁高潮呻吟视频| 久久精品人人爽人人爽视色| 久久精品国产自在天天线| 久久久久人妻精品一区果冻| 热99国产精品久久久久久7| 97在线人人人人妻| 少妇高潮的动态图| 国产精品成人在线| 激情五月婷婷亚洲| 人妻少妇偷人精品九色| 在线观看人妻少妇| 国产男女超爽视频在线观看| 黑丝袜美女国产一区| 丁香六月天网| 男女高潮啪啪啪动态图| 日韩中字成人| 全区人妻精品视频| 国产又爽黄色视频| 热re99久久精品国产66热6| 在线观看www视频免费| 人妻 亚洲 视频| 欧美xxⅹ黑人| 色婷婷av一区二区三区视频| 亚洲欧美日韩卡通动漫| 香蕉精品网在线| 亚洲国产日韩一区二区| 欧美精品一区二区免费开放| 丰满迷人的少妇在线观看| 欧美丝袜亚洲另类| 国产毛片在线视频| 91国产中文字幕| 久久久久人妻精品一区果冻| av不卡在线播放| 乱码一卡2卡4卡精品| 熟女电影av网| 啦啦啦在线观看免费高清www| 一本大道久久a久久精品| 丝袜喷水一区| 国产有黄有色有爽视频| 日本午夜av视频| 日本av手机在线免费观看| 久久久精品区二区三区| 18禁国产床啪视频网站| 最新的欧美精品一区二区| www.熟女人妻精品国产 | 免费女性裸体啪啪无遮挡网站| 成人国语在线视频| 又黄又粗又硬又大视频| 免费少妇av软件| 亚洲av日韩在线播放| 少妇的丰满在线观看| 性色avwww在线观看| av在线app专区| 国产男女超爽视频在线观看| 亚洲色图 男人天堂 中文字幕 | 国产黄色视频一区二区在线观看| 国产xxxxx性猛交| 亚洲av日韩在线播放| 大码成人一级视频| 亚洲av免费高清在线观看| 日本91视频免费播放| 免费观看在线日韩| 性色av一级| h视频一区二区三区| 天天躁夜夜躁狠狠久久av| 国产精品人妻久久久影院| 亚洲欧美中文字幕日韩二区| 精品一区二区三卡| www.av在线官网国产| 香蕉国产在线看| 啦啦啦中文免费视频观看日本| 亚洲精品456在线播放app| 久久精品国产亚洲av涩爱| 国产极品天堂在线| 搡老乐熟女国产| 激情视频va一区二区三区| 免费在线观看完整版高清| 午夜av观看不卡| 在线观看美女被高潮喷水网站| 欧美亚洲日本最大视频资源| 99国产精品免费福利视频| 看非洲黑人一级黄片| 一区二区三区四区激情视频| 纯流量卡能插随身wifi吗| 黑丝袜美女国产一区| 十分钟在线观看高清视频www| 欧美国产精品va在线观看不卡| 成年动漫av网址| av在线观看视频网站免费| 日本爱情动作片www.在线观看| 高清毛片免费看| 寂寞人妻少妇视频99o| 亚洲内射少妇av| 亚洲成人一二三区av| 国产免费现黄频在线看| 亚洲美女视频黄频| 最近手机中文字幕大全| 午夜免费男女啪啪视频观看| 美女主播在线视频| 精品国产国语对白av| 日韩视频在线欧美| 亚洲精品久久久久久婷婷小说| 免费观看性生交大片5| 亚洲图色成人| 欧美xxⅹ黑人| 亚洲av福利一区| 韩国高清视频一区二区三区| 看十八女毛片水多多多| 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩av久久| 80岁老熟妇乱子伦牲交| av播播在线观看一区| 毛片一级片免费看久久久久| 熟妇人妻不卡中文字幕| 色网站视频免费| 中文精品一卡2卡3卡4更新| 搡老乐熟女国产| 女的被弄到高潮叫床怎么办| 日本免费在线观看一区| 午夜久久久在线观看| 国产精品久久久久久精品电影小说| 午夜福利,免费看| 亚洲av男天堂| 欧美xxxx性猛交bbbb| 成年人免费黄色播放视频| 十分钟在线观看高清视频www| 水蜜桃什么品种好| 深夜精品福利| 美女国产视频在线观看| 97人妻天天添夜夜摸| 岛国毛片在线播放| 一本—道久久a久久精品蜜桃钙片| 少妇 在线观看| 色94色欧美一区二区| 男女边吃奶边做爰视频| 女人被躁到高潮嗷嗷叫费观| 中文字幕人妻丝袜制服| 午夜福利影视在线免费观看| videos熟女内射| 日产精品乱码卡一卡2卡三| 18禁在线无遮挡免费观看视频| 精品国产一区二区三区四区第35| 亚洲美女黄色视频免费看| av卡一久久| 男人添女人高潮全过程视频| 波多野结衣一区麻豆| 亚洲精品国产av成人精品| 欧美激情 高清一区二区三区| 亚洲精品乱久久久久久| 欧美3d第一页| 午夜福利在线观看免费完整高清在| 黑人猛操日本美女一级片| 国产成人欧美| 美女脱内裤让男人舔精品视频| 中文字幕最新亚洲高清| 久久久久久人妻| 中文天堂在线官网| 亚洲精品av麻豆狂野| 精品一品国产午夜福利视频| 哪个播放器可以免费观看大片| av不卡在线播放| 亚洲,欧美,日韩| 一本—道久久a久久精品蜜桃钙片| a级毛色黄片| 插逼视频在线观看| 青春草亚洲视频在线观看| 亚洲图色成人| 日韩一本色道免费dvd| 99久久综合免费| 亚洲精华国产精华液的使用体验| 免费播放大片免费观看视频在线观看| 国产片特级美女逼逼视频| 黄色配什么色好看| 一边亲一边摸免费视频| 韩国高清视频一区二区三区| 国产男女超爽视频在线观看| 欧美少妇被猛烈插入视频| 中文字幕免费在线视频6| 男的添女的下面高潮视频| 日日爽夜夜爽网站| 亚洲国产看品久久| 欧美日韩精品成人综合77777| 人妻少妇偷人精品九色| 五月开心婷婷网| 黑人猛操日本美女一级片| 亚洲精品一区蜜桃| 国产极品天堂在线| 日韩免费高清中文字幕av| 不卡视频在线观看欧美| 2018国产大陆天天弄谢| 在线亚洲精品国产二区图片欧美| 久久久亚洲精品成人影院| 国产免费一区二区三区四区乱码| 国产成人午夜福利电影在线观看| 精品久久久久久电影网| 韩国av在线不卡| 中文字幕人妻丝袜制服| 免费大片18禁| 精品国产一区二区三区四区第35| 国产国语露脸激情在线看| 日韩 亚洲 欧美在线| 99精国产麻豆久久婷婷| 色视频在线一区二区三区| 日韩电影二区| 99香蕉大伊视频| 免费黄色在线免费观看| 欧美日韩一区二区视频在线观看视频在线| 国产成人av激情在线播放| 免费人成在线观看视频色| 亚洲伊人久久精品综合| 久久精品久久久久久久性| 男女下面插进去视频免费观看 | 免费看不卡的av| 亚洲av男天堂| 女人被躁到高潮嗷嗷叫费观| 国产极品天堂在线| 青青草视频在线视频观看| 亚洲第一av免费看| 视频区图区小说| 久久久精品区二区三区| 交换朋友夫妻互换小说| 国产免费一区二区三区四区乱码| 亚洲国产精品999| 大片电影免费在线观看免费| 成年av动漫网址| av在线观看视频网站免费| 春色校园在线视频观看| 国产成人精品福利久久| 久久久久人妻精品一区果冻| 亚洲丝袜综合中文字幕| 99热网站在线观看| 久久久精品区二区三区| 欧美精品一区二区大全| 好男人视频免费观看在线| 99国产精品免费福利视频| 国产亚洲午夜精品一区二区久久| 丝袜喷水一区| 日日撸夜夜添| 国产精品久久久久久久久免| 久久免费观看电影| 热re99久久国产66热| 人妻一区二区av| 午夜久久久在线观看| 少妇精品久久久久久久| 久久国产精品大桥未久av| 国产欧美另类精品又又久久亚洲欧美| 狠狠精品人妻久久久久久综合| 亚洲一区二区三区欧美精品| 国产精品熟女久久久久浪| 精品国产一区二区三区久久久樱花| 亚洲欧美日韩卡通动漫| 伦精品一区二区三区| 一级片'在线观看视频| 极品少妇高潮喷水抽搐| 日韩中字成人| 成人二区视频| 国产极品天堂在线| 最近2019中文字幕mv第一页| 制服丝袜香蕉在线| 性色avwww在线观看| 久久久久久久久久久免费av| 美女内射精品一级片tv| 99re6热这里在线精品视频| av一本久久久久| 久久精品aⅴ一区二区三区四区 | 国产免费福利视频在线观看| 久久精品国产综合久久久 | 青春草亚洲视频在线观看| 久久这里有精品视频免费| 国产国拍精品亚洲av在线观看| 日日啪夜夜爽| 国产麻豆69| 久久久久久人人人人人| 一级,二级,三级黄色视频| 亚洲精品国产av成人精品| 色94色欧美一区二区| 少妇猛男粗大的猛烈进出视频| 黄色一级大片看看| 久久精品久久精品一区二区三区| 免费看不卡的av| √禁漫天堂资源中文www| 巨乳人妻的诱惑在线观看| 日本vs欧美在线观看视频| 香蕉精品网在线| 日本色播在线视频| 婷婷成人精品国产| 成人黄色视频免费在线看| 国产精品久久久久久精品古装| 亚洲美女黄色视频免费看| 美女视频免费永久观看网站| 制服诱惑二区| 久久精品久久久久久久性| 亚洲精品,欧美精品| a级毛色黄片| 在线观看免费日韩欧美大片| 麻豆精品久久久久久蜜桃| 免费高清在线观看视频在线观看| 熟女人妻精品中文字幕| 纵有疾风起免费观看全集完整版| av女优亚洲男人天堂| 男女免费视频国产| 五月天丁香电影| 免费观看在线日韩| 2018国产大陆天天弄谢| 丝袜喷水一区| 又粗又硬又长又爽又黄的视频| 成年女人在线观看亚洲视频| 中国三级夫妇交换| 国产精品麻豆人妻色哟哟久久| 国产一区亚洲一区在线观看| 91成人精品电影| 日韩精品有码人妻一区| 熟女人妻精品中文字幕| 女性生殖器流出的白浆| 久久热在线av| 国产探花极品一区二区| 综合色丁香网| 午夜福利视频在线观看免费| 欧美精品一区二区免费开放| 免费黄色在线免费观看| av网站免费在线观看视频| 9191精品国产免费久久| 熟女av电影| 啦啦啦在线观看免费高清www| 日韩精品免费视频一区二区三区 | 国内精品宾馆在线| 亚洲欧美精品自产自拍| 在线 av 中文字幕| 国产免费一级a男人的天堂| 国产成人aa在线观看| 亚洲精品自拍成人| 国产精品国产三级国产专区5o| 只有这里有精品99| 热re99久久精品国产66热6| 国产精品一国产av| 中文字幕亚洲精品专区| 亚洲欧美成人综合另类久久久| 久久精品国产鲁丝片午夜精品| 99久久中文字幕三级久久日本| 国产一区有黄有色的免费视频| 久久久久视频综合| 亚洲丝袜综合中文字幕| 亚洲av电影在线观看一区二区三区| 国产又色又爽无遮挡免| 男女高潮啪啪啪动态图| 久久精品熟女亚洲av麻豆精品| 男女边吃奶边做爰视频| kizo精华| 最近中文字幕高清免费大全6| 男女午夜视频在线观看 | 丰满饥渴人妻一区二区三| 日产精品乱码卡一卡2卡三| 免费看av在线观看网站| 免费播放大片免费观看视频在线观看| 亚洲国产精品国产精品| 自拍欧美九色日韩亚洲蝌蚪91| 久久这里有精品视频免费| 成人影院久久| 黄色配什么色好看| 王馨瑶露胸无遮挡在线观看| 插逼视频在线观看| 老司机影院成人| 亚洲精品乱码久久久久久按摩| 免费看av在线观看网站| 日韩精品有码人妻一区| 日日啪夜夜爽| 777米奇影视久久| 午夜福利乱码中文字幕| 国产高清国产精品国产三级| 亚洲精品美女久久av网站| 少妇猛男粗大的猛烈进出视频| 人体艺术视频欧美日本| 久久久久精品性色| 精品人妻偷拍中文字幕| 免费看av在线观看网站| 日本wwww免费看| 亚洲第一av免费看| 中国国产av一级| av有码第一页| 久久人人爽av亚洲精品天堂| 国产在线视频一区二区| 精品久久久精品久久久|