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

    基于融合評價指標和神經網絡的刀具磨損預測

    2021-11-01 09:32:50秦國華高杰葉海潮姜國杰黃帥賴曉春
    兵工學報 2021年9期
    關鍵詞:融合信號

    秦國華, 高杰, 葉海潮, 姜國杰, 黃帥, 賴曉春

    (1.南昌航空大學 航空制造工程學院, 江西 南昌 330063;2.北京航空材料研究院, 北京 100095;3.江西省教育廳 江西省教育國際合作與教師發(fā)展中心, 江西 南昌 330083)

    0 引言

    刀具磨損過快必將增加刀具的消耗,影響加工質量,降低生產效率和提高加工費用[1-2]。在研究車刀刀具磨損試驗中,需要對車削過程中的車刀狀態(tài)進行監(jiān)控,通過接收安裝在刀柄上的加速度傳感器和聲發(fā)射傳感器收集的信號來預測刀具的磨損狀態(tài)。但是這種信號是轉換信號換算得來的。這種間接獲得的信號無法直接用于控制端的輸出[3],因為它包含了外界環(huán)境、機床本身產生的噪聲等干擾。為此,許多專家學者傾注精力盡可能多地去除這些噪聲信號,保留有用信號。

    小波變換是一種信號的時間- 尺度(時間- 頻率)分析方法,具有多分辨率分析的特點,而且在時域和頻域都具有表征信號局部特征的能力[4]。近年來,利用小波變換進行分解重構去噪成為研究的一大熱點。王普等[5]提出一種分層自適應小波閾值降噪法:將振動信號進行小波分解,獲取各分解層的小波系數,保留低頻信號的小波系數;對高頻信號的小波系數進行分層自適應閾值處理,最后對處理后的小波系數重構得到降噪后的信號。楊恢先等[6]提出一種介于軟、硬閾值函數之間的閾值函數,但它對高頻信號會產生過度去噪的現(xiàn)象。黎鎖平等[7]根據噪聲香農熵與含噪信號香農熵之比,提出了一種最優(yōu)分解層數的判斷方法。Bonda等[8]采用改進的小波變換和希爾伯特- 黃變換相結合的方法對車削信號進行去噪處理。Ge等[9]在采用小波閾值去噪時,根據信號幅頻特性的不同對信號進行降噪處理,抑制白噪聲的能力較強。Beale等[10]提出一種小波包自適應算法,通過融合非傳統(tǒng)的噪聲估計方法、閾值選擇和閾值應用構造,增強去噪性能。綜上所述可見,在信號降噪領域中小波閾值去噪應用最廣泛。

    影響小波閾值去噪的主要因素在于小波閾值函數、小波閾值以及分解層數的選擇。目前研究大多集中在小波閾值以及閾值函數的改進或選擇上,沒有把閾值選取與分解層數和小波基函數相聯(lián)系。針對這一問題,本文提出一種改進方案,即在選取閾值的同時找尋與其最適分解層數及小波基函數,并在均方根誤差RMSE、信噪比SNR、平滑度r等傳統(tǒng)評價參數基礎上引入熵值法求權重的方法,構造出新的融合指標對降噪效果進行評價,以此獲取分解層數、小波基函數及閾值規(guī)則最優(yōu)組合。在提取出刀具磨損信號特征值后,建立基于切削工藝參數的刀具磨損神經網絡預測方法。

    1 小波去噪方法

    由于人為以及外界等各種因素干擾下,儀器采集到的信號中包含了真實信號和噪聲信號。從而采集的信號數據可以表示為

    f(i)=s(i)+y(i),i=1,2,…,M,

    (1)

    式中:f(i)為采集的信號;s(i)為真實信號;y(i)為噪聲信號;M為信號長度。

    在小波閾值去噪過程和方法[11]中,分解層數、小波基函數和閾值選擇函數是影響去噪效果的重要因素[12]。如果閾值太小,則信號降噪不完全;如果閾值過大,則會刪除有用的信號,影響準確度[13]。經典的閾值有以下4種:

    1)通用閾值(sqtwolog規(guī)則)。對小波分解采用相同的閾值,得到各個層次的小波系數。通用閾值法的閾值Ts描述為

    (2)

    式中:σ為噪聲信號的均方差;n為小波系數向量長度。

    2)Stein無偏風險閾值(rigrsure規(guī)則)。對每個閾值求出對應的風險值,選出其中最小的。設w為一向量,其元素wk∈w為小波系數的平方并按照從大到小的順序排列,再設一個風險向量R,其元素rt∈R為

    (3)

    以R元素中的最小值rb=min {r1,r2,…,rn}作為風險值,然后根據rb的下標變量b求出對應的wmin=wb,則Stein無偏風險閾值Tr為

    (4)

    3)啟發(fā)式閾值(heursure規(guī)則)。啟發(fā)式閾值法是前兩種閾值規(guī)則的綜合,以此選擇最優(yōu)的預測變量閾值。啟發(fā)式閾值Th的計算方法如下:

    (5)

    4)極大極小閾值(minimaxi規(guī)則)。采用極大極小的原則來選擇閾值,產生一個最小均方誤差的極值。極大極小閾值Tm的計算方法如下:

    (6)

    2 去噪評價指標

    判斷小波去噪效果的指標主要包括均方根誤差RMSE、信噪比SNR、平滑度r等[14]。它們的具體計算公式分別如下:

    (7)

    (8)

    (9)

    值得注意的是,在信號處理中,對于機械類故障所產生的非穩(wěn)態(tài)信號,分解尺度的選擇一般從j=2開始,這是因為j=1適用于穩(wěn)態(tài)信號的分解。

    傳統(tǒng)評價方式是單個或多個指標分別判斷降噪效果,在真實信號未知情況下,往往誤差很大。因此,應找尋一個更加全面準確描述信號去噪效果的評價指標,以此評價指標來選擇小波分解層數和閾值的最佳組合。

    從不同角度去判斷去噪效果,是融合指標選擇的關鍵。在實際采集的信號中,真實信號是未知的,導致互相關系數計算的準確度不高,且互相關系數與均方根誤差的定義有很大的重疊性。另一方面,由于無法獲取真實信號,也就無法直接描述采集信號的去噪效果,故真實信號一般近似為采集信號與背景噪聲信號之差。因此,選用均方根誤差、信噪比及平滑度3個指標來構造融合指標。

    由于每個指標都有不同的單位和意義,為了方便處理,采用極大極小標準化方法對各指標進行處理,使得各指標值均控制在[0.1,0.9]范圍內,具體如下:

    (10)

    (11)

    (12)

    然后根據熵值法確定各指標的權重,即

    (13)

    (14)

    dψ=1-Eψ,

    (15)

    (16)

    式中:Pψ(j)表示概率;ψ表示均方根誤差RMSE、信噪比SNR和平滑度r;Eψ表示指標熵值;m表示指標數量;dψ表示信息熵冗余值;Wψ表示權重比。

    由此,可構造第j尺度下的融合指標T為

    T(j)=WRMSECRMSE(j)+WSNRCSNR(j)+WrCr(j),

    (17)

    式中:T值代表所選參數下降噪效果的好壞,T值越大表示降噪效果越好,即T值最大所對應的為最優(yōu)參數。

    3 模擬信號的去噪實驗

    3.1 仿真信號的構建

    通過3個不同頻率的正弦信號和一個低頻的趨勢信號,構造出一個原始不含噪聲信號s(i),其表達式為

    s(i)=3sin (3 400πi)+2sin (6 000πi)+
    sin (10 000πi)+0.001i.

    (18)

    然后,在原始不含噪聲信號s(i)的基礎上添加白噪聲,生成仿真含噪聲信號。這里采用MATLAB軟件的awgn()函數進行加噪,得到SNR=5 dB的含噪聲信號f(i)。這樣,原始不含噪聲信號與含噪聲信號如圖1(a)、圖1(b)所示。

    圖1 不含噪聲信號與含噪聲信號Fig.1 Signals without and with noise

    3.2 小波基函數對降噪的影響

    采用N階小波基函數dbN對含噪聲信號f(i)進行小波閾值降噪,其他降噪參數設置如表1所示。

    表1 小波閾值降噪參數設置

    表2所示為選用表1各項指標參數計算得到的SNR、r、RMSE和構造的融合指標T數據。表2中,dbN表示N階小波基函數。從表2中可見各評價參數的變化趨勢:在傳統(tǒng)指標評價方式下,小波基函數在db7時降噪效果最佳,而構造的融合指標T也在小波基函數db7處取得最佳值。表明在此降噪參數下降噪效果最佳,與傳統(tǒng)評價指標結果一致,符合實際情況,從側面證明了融合指標T的準確性。

    表2 不同小波基函數dbN下去噪指標

    3.3 分解尺寸對降噪的影響

    在小波基函數為db7、其他參數與表1保持一致的基礎上,再一次對加入噪聲強度為SNR=5 dB的含噪聲信號Sj(i)進行去噪,結果如表3所示。

    表3 不同分解尺度下去噪指標

    從信噪比指標評價定義可知:分解尺度j=9時降噪效果最佳;在平滑度定義下,分解尺度j=9時降噪效果為最佳;在均方根誤差定義下,分解尺度j=8時的降噪效果最佳。因此,無法明確指出最佳的分解尺度。此時需要借助去噪后的圖像與不含噪聲信號圖進行辨別,但存在研究者主觀意識的判斷,導致誤差存在。

    3.4 最佳降噪參數的選擇

    從構造的融合指標T可知,分解尺度j=9時,其值最大。根據融合指標的定義可知,j=9為最佳分解尺度。分別畫出不含噪聲信號的時域圖和頻域圖,以及選取j=8和j=9時去噪后的信號圖與頻域圖,如圖2~圖4所示。

    圖2 不含噪聲信號Fig.2 Signal without noise

    圖3 分解尺度j=8的去噪聲信號Fig.3 Denoised signal for decomposition scale j=8

    圖4 分解尺度j=9的去噪聲信號Fig.4 Denoised signal for decomposition scale j=9

    通過比較時域圖和頻域圖中的細節(jié)部分可知:當分解尺度j=9時,降噪后的信號更接近真實信號;當分解尺度j=8時,降噪的效果不如分解尺度j=9.這一結果表明融合指標T能夠有效且準確地在一組數據中挑選出一個最佳的數據,即依據融合指標T對去噪效果進行評估,能夠獲得一組消噪效果最優(yōu)的降噪參數。

    3.5 基于融合指標的降噪參數優(yōu)選方法

    采用融合指標T能夠優(yōu)選出降噪參數,其優(yōu)選方法和流程構造如圖5所示。

    圖5 最優(yōu)降噪參數選擇流程圖Fig.5 Selection flowchart of optimal noise reduction parameters

    去噪參數的優(yōu)選步驟具體如下:

    步驟1初始化信號的分解尺度j=2.

    步驟2初始化小波基函數的階數N=2.

    步驟3對含噪聲信號f(i)采用小波基函數dbN、分解尺度j進行小波分解。

    步驟4對分解所得到的細節(jié)系數分別進行sqtwolog等4種閾值以及軟閾值函數處理。

    步驟5輸出4個去噪聲后的信號Sj(i),計算相應的信噪比、平滑度、均方根誤差3個傳統(tǒng)評價指標以及融合指標T.

    步驟6判斷當前小波基函數的階數是否超出最大階數,即N>9?若是則轉步驟7,否則N=N+1,并轉步驟3.

    步驟7在分解尺度j下,從db2~db9中依次選取小波基函數的最佳階數。

    步驟8判斷分解尺度是否超出最大值,即j>9?若否,則j=j+1,轉步驟2,否則轉步驟9.

    步驟9優(yōu)選過程結束,獲得最佳閾值、階數和分解尺度。

    針對圖1(b)的含噪聲信號,利用圖5的優(yōu)選過程,可得表4中的數據。由于T=0.734 942 99為最大值,故最佳參數分別為小波基函數db9,分解尺度j=9,閾值規(guī)則選用為rigrsure,尺度系數為sln.

    表4 每一分解尺度下最佳參數

    由于圖1(b)中仿真實驗添加的高斯白噪聲信號過于單一,而監(jiān)測刀具磨損采集的信號數據中含有的噪聲較為復雜,需要從刀具磨損試驗中進一步驗證本文方法的可行性和有效性。

    4 實際信號的去噪應用

    4.1 測量磨損量

    為進一步驗證本文方法對含噪聲信號進行去噪聲處理的可行性和有效性,在切削加工過程中,通過安裝在刀柄上的振動傳感器采集信號,對采集到的信號進行相應的去噪聲,即信號預處理,通過提取去噪聲后的信號特征向量,建立刀具磨損量與特征向量之間的函數關系。

    切削試驗前,首先要測量出刀具后刀面的磨損量VB. 本文試驗選用日本基恩士公司生產的VHX- 600型超景深三維顯微鏡來測量刀具磨損量,如圖6所示。為了準確地獲取刀具磨損量,避免測量過程中可能存在的誤差,每次均測量3次后取平均值,作為刀具后刀面的磨損量。

    圖6 刀具磨損的測量Fig.6 Measurement of tool wear

    4.2 采集信號

    工件為45號鋼棒料,切削過程使用切削液。車刀選用YT15可轉位硬質合金刀片,初始磨損量為0.276 mm,前角為5°,后角為6°,主偏角為45°,副偏角為45°,刃傾角為0°.

    試驗采用大連機床廠生產的CKA6150型數控車床,加工時選用的工藝參數分別為主軸轉速520 r/min、進給量0.03 mm/r、背吃刀量1.0 mm. 值得一提的是,所有切削試驗中車削的長度均為100 mm. 傳感器選用江蘇聯(lián)能電子公司生產的CA-YD-186型加速度傳感器,該傳感器帶有磁性底座,方便拆卸。傳感器位置一般盡量靠近切削部位,通常選擇刀柄的側面和上方,這里將傳感器安裝在靠近可拆卸刀片的壓板附近,具體如圖7所示。

    圖7 刀具磨損振動信號采集試驗裝置Fig.7 Acquisition experimental device of vibration signal for tool wear

    振動信號由江蘇聯(lián)能電子公司生產的YE6231型四通道動態(tài)數據采集儀采集,如圖7(a)所示,采樣頻率設置為12 kHz,采集到的原始信號如圖8所示。

    圖8 刀具磨損原始信號Fig.8 Time-frequency diagram of original signal

    4.3 降噪處理

    由基于融合指標的降噪參數優(yōu)選方法可知,當分解層次為j=2、小波基函數為db5、閾值為sqtwolog時,降噪效果最佳。降噪(簡稱融合指標降噪法)后的刀具磨損信號的時域與頻域圖如圖9所示。

    圖9 融合指標法降噪聲后的信號Fig.9 Denoised signal obtained by composite index method

    根據文獻[3]建立的基于信噪比評價的閾值優(yōu)選小波去噪法,得到的最佳降噪參數分別為分解層次j=7,db9小波基,以及sqtwolog通用閾值。同樣,采用該組參數對采集到的刀具磨損原始信號進行去噪聲(簡稱信噪比降噪法),得到的時域與頻域圖如圖10所示。

    圖10 信噪比法降噪聲后的信號Fig.10 Denoised signal obtained by SNR method

    對比圖9和圖10可知,融合指標法和信噪比法均能夠有效地過濾掉高頻噪聲信號部分;但相比于信噪比法,融合指標法更加完整地保留了低頻信號部分。顯然,融合指標法的降噪效果更為明顯。

    5 刀具磨損的預測方法

    5.1 網絡結構

    利用融合指標法對刀具磨損的采集信號進行去噪聲后,能夠更精確地選取出特征值,進而可以通過建立基于特征值的刀具磨損神經網絡預測模型,實現(xiàn)切削過程中工藝參數對磨損量的影響。

    利用選定的刀具車削工件時,影響刀具磨損的重要因素為切削速度、進給量、背吃刀量等工藝參數,從而信號特征值也不一樣。因此,刀具磨損的BP神經網絡結構如圖11所示,其中:輸入層神經元x=[x1,x2,…,xi,…,xu]T,包括提取的2個信號特征以及3個切削用量,即u=5;輸出層即為刀具磨損量,故輸出層僅有v=1個神經元;根據Kolmogorov定理,綜合考慮網絡的學習速度和泛化能力,根據通用的經驗公式[15-16]確定隱含層的k個神經元,即

    圖11 網絡結構圖Fig.11 Network structure

    k=2u+1.

    (19)

    5.2 訓練樣本

    工件初始直徑為50 mm,車出光面后直徑為48.80 mm. 隨后進行后續(xù)的16組車削試驗,切削過中使用切削液。由于主軸轉速為320 r/mm≤x1≤1 120 r/mm,進給量為0.05 r/mm≤x2≤0.15 r/mm、背吃刀量為0.3 mm≤x3≤1.0 mm,采用正交試驗設計方法確定輸入樣本,如表5所示。

    表5 刀具磨損車削試驗正交表

    根據表5的數據,可獲得16組刀具磨損的車削試驗,采集到各自的刀具磨損信號。從每組信號中選取2組樣本信號,對其使用融合指標降噪法進行去噪聲,采用小波包分解技術對降噪聲后的信號提取與刀具磨損有關特征向量。同時,測量出16組試驗后的刀具磨損量作為訓練樣本,從而可得輸入樣本和訓練樣本如表6所示。

    表6 輸入樣本及其相應的訓練樣本

    為了使各個輸入層的變量都處于[0,1],需采用(20)式對輸入樣本進行歸一化處理,以提高神經網絡的識別精度,即

    (20)

    式中:x′i為第i個切削參數的輸入樣本;xmax和xmin分別為樣本的最大值和最小值。

    5.3 網絡訓練

    輸入樣本歸一化后,即可開始進行網絡訓練。訓練中,隱藏層11個神經元的激活函數選擇雙曲正切函數,實現(xiàn)不同樣本的空間分割,而輸入層1個神經元的激活函數選擇線性函數,用于輸出網絡的識別結果。

    初始化網絡權系數在[0,1]之間,學習誤差為0.000 000 1,最大訓練周期10 000. 將32個樣本隨機輸入網絡進行訓練,網絡在107步收斂到要求精度,神經網絡訓練過程誤差曲線如圖12所示。

    圖12 網絡訓練過程Fig.12 Training process of network

    網絡訓練好后,隨機抽取11組不同磨損狀態(tài)的刀具車削試驗采集的待測樣本,提取其特征量,對其歸一化后,輸入網絡進行磨損量的預測,其目標輸出和實際輸出部分結果如表7所示。

    通過表7可以看出,刀具磨損量的模型預測結果非常接近于試驗測量值,最大的預測誤差也沒超出6.0%. 由此可見,經融合指標降噪法挑選出來的降噪參數對采集到的振動信號進行去噪聲處理后,將其輸入到刀具磨損量的預測模型中,其模型預測精度高,可靠性強。

    表7 網絡測試集

    6 結論

    1)在3個傳統(tǒng)指標SNR、r和RMSE的基礎上,引入熵值法獲得3個指標的權重比,構造了融合指標T. 相比單個傳統(tǒng)指標的去噪效果,經融合指標去噪聲后的信號更符合實際信號。

    2)提出了降噪參數的優(yōu)選方法。對于含有復雜噪聲的刀具磨損信號,基于多指標融合評價的降噪參數優(yōu)選方法能夠更加客觀準確地挑選出最佳降噪參數,不僅有效地過濾噪聲的高頻部分,也能比較完整地保留低頻部分。

    3)通過刀具磨損的信號特征,建立了工藝參數與刀具磨損之間的神經網絡預測模型。試驗結果表明,模型預測精度高,最大誤差不超過6.0%,基于多指標融合評價的最佳降噪參數更加有利于提取與刀具磨損密切相關的信號特征量。

    猜你喜歡
    融合信號
    一次函數“四融合”
    村企黨建聯(lián)建融合共贏
    融合菜
    信號
    鴨綠江(2021年35期)2021-04-19 12:24:18
    從創(chuàng)新出發(fā),與高考數列相遇、融合
    完形填空二則
    寬窄融合便攜箱IPFS500
    《融合》
    孩子停止長個的信號
    基于LabVIEW的力加載信號采集與PID控制
    av天堂在线播放| 精品熟女少妇八av免费久了| 丰满迷人的少妇在线观看| www日本在线高清视频| 久久人人97超碰香蕉20202| kizo精华| 一级,二级,三级黄色视频| 日韩中文字幕视频在线看片| 不卡av一区二区三区| 97在线人人人人妻| 国产亚洲精品久久久久5区| 久久午夜亚洲精品久久| 黄色怎么调成土黄色| 亚洲色图av天堂| 免费黄频网站在线观看国产| 久久人妻福利社区极品人妻图片| 久久久久视频综合| 不卡av一区二区三区| 精品亚洲成a人片在线观看| 美女高潮到喷水免费观看| 亚洲色图 男人天堂 中文字幕| 国产高清videossex| 9热在线视频观看99| 欧美在线黄色| 最近最新免费中文字幕在线| 久久影院123| 一边摸一边抽搐一进一出视频| 欧美精品高潮呻吟av久久| 天天添夜夜摸| 日韩中文字幕视频在线看片| 水蜜桃什么品种好| 日韩欧美三级三区| 国产精品偷伦视频观看了| 久久ye,这里只有精品| 美女福利国产在线| 中文字幕av电影在线播放| 俄罗斯特黄特色一大片| 国产精品自产拍在线观看55亚洲 | 亚洲av美国av| 国产精品自产拍在线观看55亚洲 | av福利片在线| 欧美日韩亚洲国产一区二区在线观看 | 国产淫语在线视频| 亚洲精品美女久久av网站| 纵有疾风起免费观看全集完整版| 免费观看人在逋| 久久中文字幕人妻熟女| 久久99一区二区三区| 免费看十八禁软件| 精品熟女少妇八av免费久了| 亚洲国产欧美在线一区| 黄色丝袜av网址大全| 国产成人av激情在线播放| 久久青草综合色| 一本久久精品| 九色亚洲精品在线播放| 狠狠精品人妻久久久久久综合| 丁香六月天网| 欧美老熟妇乱子伦牲交| 国产熟女午夜一区二区三区| 9191精品国产免费久久| 国产亚洲精品第一综合不卡| 男女之事视频高清在线观看| 法律面前人人平等表现在哪些方面| 成年人午夜在线观看视频| 欧美在线黄色| 久久久久久久久久久久大奶| 一级片'在线观看视频| 美女国产高潮福利片在线看| 亚洲,欧美精品.| 99riav亚洲国产免费| 国产高清videossex| 亚洲熟女毛片儿| 久久中文看片网| 久久性视频一级片| 亚洲成人手机| 少妇精品久久久久久久| 久久精品91无色码中文字幕| 9色porny在线观看| 欧美亚洲日本最大视频资源| 国产成人精品在线电影| 国产精品亚洲av一区麻豆| 高清欧美精品videossex| 无遮挡黄片免费观看| 亚洲欧美一区二区三区黑人| 亚洲国产中文字幕在线视频| 热re99久久精品国产66热6| 亚洲精品美女久久久久99蜜臀| 成人黄色视频免费在线看| 中文字幕av电影在线播放| 丝瓜视频免费看黄片| 亚洲av日韩精品久久久久久密| 国产国语露脸激情在线看| 精品卡一卡二卡四卡免费| 在线观看人妻少妇| 午夜福利在线观看吧| 国产精品一区二区在线不卡| 国产熟女午夜一区二区三区| 亚洲色图综合在线观看| 日韩成人在线观看一区二区三区| 久久精品成人免费网站| 国产成人免费无遮挡视频| 黑人巨大精品欧美一区二区蜜桃| 免费女性裸体啪啪无遮挡网站| 老司机福利观看| 亚洲一区中文字幕在线| 美女高潮喷水抽搐中文字幕| 一区二区av电影网| 精品福利观看| 亚洲熟妇熟女久久| 一本色道久久久久久精品综合| 激情在线观看视频在线高清 | av超薄肉色丝袜交足视频| 国产aⅴ精品一区二区三区波| 精品国产亚洲在线| 91成年电影在线观看| 日韩大片免费观看网站| 精品少妇一区二区三区视频日本电影| 19禁男女啪啪无遮挡网站| 亚洲国产av影院在线观看| 亚洲自偷自拍图片 自拍| 精品视频人人做人人爽| 麻豆国产av国片精品| 成人精品一区二区免费| avwww免费| 久久人妻熟女aⅴ| 国产成人免费观看mmmm| 欧美精品av麻豆av| 成人国产av品久久久| 另类精品久久| 午夜福利,免费看| 汤姆久久久久久久影院中文字幕| 欧美中文综合在线视频| 我要看黄色一级片免费的| av天堂在线播放| 色94色欧美一区二区| 两性夫妻黄色片| 亚洲欧洲日产国产| 久久青草综合色| 久久久久网色| 精品国产超薄肉色丝袜足j| 两个人看的免费小视频| 岛国在线观看网站| 我要看黄色一级片免费的| 狠狠精品人妻久久久久久综合| 美女高潮喷水抽搐中文字幕| 最近最新中文字幕大全免费视频| 99久久精品国产亚洲精品| 成人精品一区二区免费| 999久久久国产精品视频| 捣出白浆h1v1| 欧美精品一区二区免费开放| 国产精品 欧美亚洲| 又黄又粗又硬又大视频| 91成人精品电影| 天堂动漫精品| 高清毛片免费观看视频网站 | 亚洲专区国产一区二区| 亚洲美女黄片视频| 国产一区二区三区在线臀色熟女 | 日本欧美视频一区| 香蕉丝袜av| 午夜福利在线免费观看网站| 一进一出抽搐动态| 9色porny在线观看| 97人妻天天添夜夜摸| 国产一区有黄有色的免费视频| 亚洲男人天堂网一区| 亚洲中文av在线| 男人操女人黄网站| 亚洲精品国产色婷婷电影| 夜夜骑夜夜射夜夜干| 日本欧美视频一区| 少妇被粗大的猛进出69影院| 极品教师在线免费播放| bbb黄色大片| 国产精品久久久久成人av| 黑人巨大精品欧美一区二区蜜桃| 亚洲性夜色夜夜综合| 国产极品粉嫩免费观看在线| 午夜激情久久久久久久| 水蜜桃什么品种好| 午夜福利视频在线观看免费| 老司机影院毛片| 亚洲天堂av无毛| 一区二区三区国产精品乱码| 亚洲欧美激情在线| 日韩成人在线观看一区二区三区| 久久久久精品国产欧美久久久| 亚洲精品av麻豆狂野| 欧美 日韩 精品 国产| 国产精品免费一区二区三区在线 | 亚洲男人天堂网一区| 最黄视频免费看| 我的亚洲天堂| 亚洲熟女毛片儿| 天堂动漫精品| 99精品久久久久人妻精品| 最近最新中文字幕大全电影3 | 日韩大码丰满熟妇| 伦理电影免费视频| 亚洲午夜精品一区,二区,三区| 丝瓜视频免费看黄片| 最近最新中文字幕大全免费视频| 久久国产精品人妻蜜桃| 丝瓜视频免费看黄片| 亚洲一码二码三码区别大吗| 最近最新中文字幕大全电影3 | 国产老妇伦熟女老妇高清| 一区福利在线观看| 欧美日韩亚洲综合一区二区三区_| 欧美黄色淫秽网站| av不卡在线播放| 在线观看免费视频网站a站| 极品教师在线免费播放| 欧美日韩福利视频一区二区| 色视频在线一区二区三区| 日韩中文字幕欧美一区二区| 十八禁网站网址无遮挡| 亚洲一区二区三区欧美精品| 在线播放国产精品三级| 999久久久精品免费观看国产| 90打野战视频偷拍视频| 免费久久久久久久精品成人欧美视频| 变态另类成人亚洲欧美熟女 | 99在线人妻在线中文字幕 | 日韩一卡2卡3卡4卡2021年| 在线观看舔阴道视频| 婷婷成人精品国产| 两个人免费观看高清视频| 一夜夜www| av网站在线播放免费| 757午夜福利合集在线观看| 日韩欧美一区视频在线观看| 亚洲七黄色美女视频| 日本vs欧美在线观看视频| 在线观看免费高清a一片| 12—13女人毛片做爰片一| 成年人免费黄色播放视频| 十八禁人妻一区二区| av网站在线播放免费| 久久国产精品影院| 日韩欧美一区二区三区在线观看 | 欧美精品高潮呻吟av久久| 免费高清在线观看日韩| 一本色道久久久久久精品综合| 女性被躁到高潮视频| 亚洲 国产 在线| 国产男女超爽视频在线观看| 久久久精品国产亚洲av高清涩受| 久久久国产一区二区| 超碰97精品在线观看| 欧美日本中文国产一区发布| 在线观看免费日韩欧美大片| 国产97色在线日韩免费| 欧美人与性动交α欧美软件| 一区福利在线观看| 亚洲人成电影免费在线| 国产成+人综合+亚洲专区| tube8黄色片| 亚洲性夜色夜夜综合| 大型av网站在线播放| svipshipincom国产片| 在线观看免费视频日本深夜| 国产不卡av网站在线观看| 18在线观看网站| 啦啦啦视频在线资源免费观看| 三上悠亚av全集在线观看| videosex国产| 一本色道久久久久久精品综合| av免费在线观看网站| 91麻豆精品激情在线观看国产 | 日韩欧美国产一区二区入口| 99国产综合亚洲精品| 国产男女内射视频| 妹子高潮喷水视频| 午夜福利免费观看在线| 18禁黄网站禁片午夜丰满| 午夜免费鲁丝| 蜜桃国产av成人99| 十八禁网站网址无遮挡| 亚洲熟女精品中文字幕| 国精品久久久久久国模美| 老熟女久久久| 人成视频在线观看免费观看| av天堂久久9| 亚洲三区欧美一区| a级片在线免费高清观看视频| 亚洲五月婷婷丁香| 午夜91福利影院| 两个人免费观看高清视频| 日韩大码丰满熟妇| 色播在线永久视频| 无遮挡黄片免费观看| 一夜夜www| 在线播放国产精品三级| 91字幕亚洲| 亚洲av第一区精品v没综合| 国产日韩一区二区三区精品不卡| 18禁观看日本| 国产成人精品在线电影| 一区福利在线观看| 国产色视频综合| av电影中文网址| 日韩欧美三级三区| av片东京热男人的天堂| 国产精品1区2区在线观看. | 国产一区有黄有色的免费视频| 亚洲色图 男人天堂 中文字幕| 啦啦啦在线免费观看视频4| 精品亚洲乱码少妇综合久久| 狠狠婷婷综合久久久久久88av| 高清av免费在线| 国产成人欧美| 性色av乱码一区二区三区2| 久久ye,这里只有精品| 精品国产超薄肉色丝袜足j| 99精品欧美一区二区三区四区| 亚洲人成电影免费在线| 在线观看舔阴道视频| 在线永久观看黄色视频| 国产精品美女特级片免费视频播放器 | 久久久水蜜桃国产精品网| aaaaa片日本免费| 欧美精品啪啪一区二区三区| 久久精品国产亚洲av香蕉五月 | 久久久精品国产亚洲av高清涩受| 欧美精品啪啪一区二区三区| 亚洲男人天堂网一区| 日韩一卡2卡3卡4卡2021年| av免费在线观看网站| 99精品欧美一区二区三区四区| 菩萨蛮人人尽说江南好唐韦庄| 欧美日韩精品网址| 国产在线免费精品| 一级片'在线观看视频| 成人国产一区最新在线观看| 久久国产精品大桥未久av| 大码成人一级视频| 丰满饥渴人妻一区二区三| 欧美激情高清一区二区三区| 国产aⅴ精品一区二区三区波| 涩涩av久久男人的天堂| 亚洲一区二区三区欧美精品| 国产av国产精品国产| 丝瓜视频免费看黄片| 天堂中文最新版在线下载| 午夜成年电影在线免费观看| 精品国产亚洲在线| 国产精品av久久久久免费| 精品第一国产精品| 涩涩av久久男人的天堂| 在线观看免费日韩欧美大片| 欧美 亚洲 国产 日韩一| 国产精品免费一区二区三区在线 | 精品久久久久久久毛片微露脸| 国产精品欧美亚洲77777| 午夜视频精品福利| 亚洲七黄色美女视频| 色综合婷婷激情| 久久精品国产亚洲av香蕉五月 | 女人精品久久久久毛片| 美女视频免费永久观看网站| 交换朋友夫妻互换小说| 99精品在免费线老司机午夜| 日韩视频在线欧美| 久9热在线精品视频| 女性被躁到高潮视频| 五月开心婷婷网| 久久久久久久久免费视频了| 免费少妇av软件| 如日韩欧美国产精品一区二区三区| www日本在线高清视频| 19禁男女啪啪无遮挡网站| 国产激情久久老熟女| 在线观看人妻少妇| 国产精品偷伦视频观看了| 亚洲专区字幕在线| 捣出白浆h1v1| www.精华液| 久久人人97超碰香蕉20202| 最新的欧美精品一区二区| 脱女人内裤的视频| av有码第一页| 香蕉久久夜色| 久久久久久久精品吃奶| 精品少妇一区二区三区视频日本电影| 成人手机av| 久久婷婷成人综合色麻豆| 久久久国产成人免费| 老熟妇仑乱视频hdxx| 精品熟女少妇八av免费久了| 精品亚洲乱码少妇综合久久| av又黄又爽大尺度在线免费看| 亚洲成人免费av在线播放| 成人国产av品久久久| 亚洲人成伊人成综合网2020| 丝袜喷水一区| 自拍欧美九色日韩亚洲蝌蚪91| 日本欧美视频一区| 丰满人妻熟妇乱又伦精品不卡| 蜜桃在线观看..| 国产aⅴ精品一区二区三区波| 男女床上黄色一级片免费看| 成人精品一区二区免费| 中文字幕人妻丝袜制服| a级毛片黄视频| 欧美激情高清一区二区三区| 美女福利国产在线| a级片在线免费高清观看视频| 欧美人与性动交α欧美精品济南到| 午夜91福利影院| 久久精品熟女亚洲av麻豆精品| 黄片大片在线免费观看| 99久久99久久久精品蜜桃| 宅男免费午夜| 可以免费在线观看a视频的电影网站| 亚洲成人免费电影在线观看| 亚洲av电影在线进入| 两个人看的免费小视频| 女人高潮潮喷娇喘18禁视频| tocl精华| 热99国产精品久久久久久7| 亚洲精品美女久久av网站| 正在播放国产对白刺激| 国产日韩欧美亚洲二区| 亚洲天堂av无毛| 少妇 在线观看| 香蕉丝袜av| 啦啦啦视频在线资源免费观看| 精品人妻熟女毛片av久久网站| 国产精品 欧美亚洲| 国产精品麻豆人妻色哟哟久久| 中文欧美无线码| 在线十欧美十亚洲十日本专区| 午夜福利在线免费观看网站| 欧美久久黑人一区二区| 午夜福利欧美成人| 久久久精品国产亚洲av高清涩受| 免费在线观看影片大全网站| 色精品久久人妻99蜜桃| 成人永久免费在线观看视频 | 成人影院久久| 精品一品国产午夜福利视频| 中文字幕色久视频| 亚洲专区中文字幕在线| 少妇裸体淫交视频免费看高清 | 精品免费久久久久久久清纯 | 女人精品久久久久毛片| 亚洲欧美精品综合一区二区三区| 国内毛片毛片毛片毛片毛片| xxxhd国产人妻xxx| 国产不卡av网站在线观看| 一本—道久久a久久精品蜜桃钙片| 国产高清视频在线播放一区| 日韩人妻精品一区2区三区| 欧美人与性动交α欧美软件| 美女高潮到喷水免费观看| 人人妻人人添人人爽欧美一区卜| 久久久欧美国产精品| 欧美日韩亚洲国产一区二区在线观看 | 日韩大片免费观看网站| 91字幕亚洲| 黄片大片在线免费观看| 激情在线观看视频在线高清 | 交换朋友夫妻互换小说| 999久久久精品免费观看国产| 亚洲av电影在线进入| 国产成人免费无遮挡视频| 日韩三级视频一区二区三区| 国产人伦9x9x在线观看| 久久精品熟女亚洲av麻豆精品| 亚洲成国产人片在线观看| 国产在线一区二区三区精| 精品熟女少妇八av免费久了| 亚洲精品在线美女| 亚洲国产av新网站| 日韩欧美一区视频在线观看| 窝窝影院91人妻| 免费观看人在逋| 久久国产精品影院| 三级毛片av免费| 亚洲国产看品久久| 黑人欧美特级aaaaaa片| 一个人免费在线观看的高清视频| 黄色视频在线播放观看不卡| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美av亚洲av综合av国产av| 国产一区二区 视频在线| 90打野战视频偷拍视频| 97人妻天天添夜夜摸| 69精品国产乱码久久久| 日本欧美视频一区| 色视频在线一区二区三区| 俄罗斯特黄特色一大片| 久久精品亚洲熟妇少妇任你| 男女之事视频高清在线观看| 亚洲自偷自拍图片 自拍| 免费女性裸体啪啪无遮挡网站| 一区二区三区国产精品乱码| 欧美激情极品国产一区二区三区| 黄色怎么调成土黄色| 99国产精品免费福利视频| 国产99久久九九免费精品| 日韩欧美三级三区| 两个人看的免费小视频| 欧美黄色片欧美黄色片| 国产精品欧美亚洲77777| 国产精品香港三级国产av潘金莲| 法律面前人人平等表现在哪些方面| 精品视频人人做人人爽| 国产成人欧美在线观看 | 久久久久久久大尺度免费视频| 一区在线观看完整版| 日本黄色视频三级网站网址 | 欧美人与性动交α欧美精品济南到| 国产无遮挡羞羞视频在线观看| 久久精品亚洲av国产电影网| 动漫黄色视频在线观看| 国产黄频视频在线观看| 一区二区三区激情视频| 欧美成狂野欧美在线观看| 国产精品亚洲一级av第二区| 久久精品成人免费网站| 少妇 在线观看| 成人手机av| 亚洲熟女毛片儿| 在线观看免费日韩欧美大片| 中文亚洲av片在线观看爽 | 国产黄频视频在线观看| 一区二区三区激情视频| 久久性视频一级片| 国产无遮挡羞羞视频在线观看| 国产欧美日韩一区二区精品| 纯流量卡能插随身wifi吗| 操出白浆在线播放| 老熟女久久久| 一本色道久久久久久精品综合| 亚洲国产毛片av蜜桃av| 欧美成人免费av一区二区三区 | 日韩欧美国产一区二区入口| 夜夜夜夜夜久久久久| 男女无遮挡免费网站观看| 久久ye,这里只有精品| 午夜福利乱码中文字幕| 法律面前人人平等表现在哪些方面| 亚洲精品久久成人aⅴ小说| 久久影院123| 一二三四社区在线视频社区8| 国产成人精品久久二区二区免费| 日韩欧美免费精品| 9热在线视频观看99| 高清毛片免费观看视频网站 | 国产精品美女特级片免费视频播放器 | 亚洲avbb在线观看| 久久精品国产亚洲av香蕉五月 | 精品少妇久久久久久888优播| 亚洲人成77777在线视频| 亚洲va日本ⅴa欧美va伊人久久| 精品免费久久久久久久清纯 | 少妇 在线观看| 汤姆久久久久久久影院中文字幕| 一区二区三区激情视频| av福利片在线| 人人妻,人人澡人人爽秒播| 美女扒开内裤让男人捅视频| 蜜桃在线观看..| 精品国内亚洲2022精品成人 | 亚洲av成人不卡在线观看播放网| 日韩有码中文字幕| 亚洲全国av大片| 最新在线观看一区二区三区| 国产又色又爽无遮挡免费看| 久久影院123| 国产一区二区三区综合在线观看| 精品少妇内射三级| 男女之事视频高清在线观看| 亚洲精品一卡2卡三卡4卡5卡| 制服诱惑二区| 国产无遮挡羞羞视频在线观看| 日韩 欧美 亚洲 中文字幕| 亚洲精品国产一区二区精华液| 午夜福利一区二区在线看| 成人影院久久| 十八禁网站网址无遮挡| 搡老熟女国产l中国老女人| 肉色欧美久久久久久久蜜桃| 久久国产精品影院| 操出白浆在线播放| 在线观看人妻少妇| 久久久国产一区二区| 国产精品国产av在线观看| 国产精品电影一区二区三区 | 亚洲男人天堂网一区| 99国产精品99久久久久| 肉色欧美久久久久久久蜜桃| 欧美日韩精品网址| 最黄视频免费看| av天堂在线播放| 免费久久久久久久精品成人欧美视频| 18禁国产床啪视频网站| 黑丝袜美女国产一区| 亚洲av成人一区二区三| 亚洲七黄色美女视频| 欧美成人午夜精品| 亚洲人成伊人成综合网2020| 国产又爽黄色视频| 黄色a级毛片大全视频| 老司机影院毛片| 97人妻天天添夜夜摸| 亚洲色图 男人天堂 中文字幕| 日韩欧美国产一区二区入口| 成年女人毛片免费观看观看9 | 不卡一级毛片| 亚洲成国产人片在线观看|