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

    基于增量學習的數(shù)控機床故障診斷系統(tǒng)

    2022-12-20 07:57:26張煜瑩
    計量學報 2022年11期
    關(guān)鍵詞:增量故障診斷準確率

    張煜瑩, 陸 藝, 趙 靜

    (1.中國計量大學計量測試工程學院,浙江杭州310018;2. 杭州沃鐳智能科技股份有限公司,浙江杭州310018)

    1 引 言

    作為制造業(yè)的主要設(shè)備,數(shù)控機床的加工水平直接影響機械加工、航空等眾多產(chǎn)業(yè)[1]。主軸軸承作為數(shù)控機床中的關(guān)鍵零件,若長期于磨損、膠合等失效形式下工作會導(dǎo)致軸承溫度升高,故障加劇,最終影響零件加工精度或?qū)е略O(shè)備停機[2]。刀具作為與加工件直接接觸的部件,其磨損狀態(tài)直接影響加工件的精度[3]。軸承故障診斷常采用共振解調(diào)技術(shù),此方法通過對比包絡(luò)頻譜中的頻率和軸承故障特征頻率進行故障診斷[4~6]。針對刀具故障診斷,戴穩(wěn)將刀具振動數(shù)據(jù)轉(zhuǎn)化為小波尺度圖,將其輸入依照Alexnet網(wǎng)絡(luò)架構(gòu)搭建的卷積神經(jīng)網(wǎng)絡(luò)(convolutional neunral networks,CNN)模型中訓練及診斷,此方法對刀具磨損狀態(tài)判斷具有良好效果[7];陳炳旭借助DenseNet網(wǎng)絡(luò)和BILSTM記憶單元搭建卷積循環(huán)神經(jīng)網(wǎng)絡(luò)用于刀具狀態(tài)判斷,彌補卷積網(wǎng)絡(luò)對時序特征提取的不足[8]。但這些診斷方法存在2個問題:1)未考慮2種故障同時出現(xiàn)的狀況;2)這些方法均基于同一工況,實際工業(yè)現(xiàn)場中,由于工作任務(wù)的變化,機床可能處于不同轉(zhuǎn)速下,預(yù)訓練模型無法對所有轉(zhuǎn)速下的數(shù)據(jù)進行故障診斷。且預(yù)訓練模型對故障類別的涵蓋程度有限,因此需要考慮對預(yù)訓練模型進行更新優(yōu)化,使其滿足不同時期檢測需求。

    針對數(shù)控機床中主軸軸承和刀具同時出現(xiàn)故障或機床主軸轉(zhuǎn)速改變時的故障診斷問題,設(shè)計了1種基于增量學習的數(shù)控機床故障診斷系統(tǒng),該系統(tǒng)以主軸軸承和刀具作為整體研究對象。系統(tǒng)訓練過程涉及2個模型——教師模型和學生模型。教師模型為某一轉(zhuǎn)速下的故障診斷模型,學生模型為跨轉(zhuǎn)速故障診斷模型。通過增量學習使學生模型繼承教師模型的知識并學習跨轉(zhuǎn)速數(shù)據(jù)的知識,最終利用學生模型對測試數(shù)據(jù)進行測試,輸出結(jié)果。

    2 系統(tǒng)設(shè)計

    2.1 硬件設(shè)計

    系統(tǒng)硬件部分主要包括數(shù)控機床、振動加速度傳感器、恒流源、采集卡和邊緣端計算機。數(shù)控機床型號為vmc-850,傳感器型號為三軸加速度傳感器1A314E,用于采集x、y、z三個方向的振動信號;恒流源型號為YMC8204,用于輸出10 mA電流給傳感器供電;數(shù)據(jù)采集卡型號為PCI-1716,采樣頻率fs為20 kHz,最終采集所得的振動信號儲存在邊緣端計算機,用于后續(xù)模型訓練。硬件結(jié)構(gòu)示意圖如圖1所示。

    圖1 硬件結(jié)構(gòu)示意圖Fig.1 The hardware structure diagram

    2.2 軟件設(shè)計

    系統(tǒng)軟件部分主要包括數(shù)據(jù)采集模塊、故障仿真模塊以及故障診斷模塊。軟件流程圖如圖2所示。

    圖2 軟件結(jié)構(gòu)示意圖Fig.2 The software structure diagram

    數(shù)據(jù)采集模塊位于邊緣端計算機,通過傳感器等硬件利用此模塊采集主軸軸承和刀具的振動信號,并初步去除工頻干擾等影響。實際工業(yè)使用中,主軸轉(zhuǎn)速通常在1 000~6 000 r/min范圍內(nèi),而最常用轉(zhuǎn)速為3 000 r/min,因此將主軸轉(zhuǎn)速3 000 r/min下采集的數(shù)據(jù)記為數(shù)據(jù)集1,主軸轉(zhuǎn)速1 000,2 000,4 000,5 000,6 000 r/min時采集的數(shù)據(jù)記為數(shù)據(jù)集2。此模塊需要采集數(shù)據(jù)集1和數(shù)據(jù)集2,其中包括了訓練數(shù)據(jù)以及測試數(shù)據(jù)。

    試驗機床為工業(yè)現(xiàn)場機床,因此刀具磨損信號容易得到,但軸承無法人為破壞。為得到含有不同軸承故障和不同刀具磨損狀態(tài)的完整數(shù)據(jù)集,使用故障仿真模塊對主軸軸承故障進行仿真。故障仿真信號如式(1),式(2)所示[9]:

    Y(t)=S(t)+r(t)

    (1)

    (2)

    當y0=20,ζ=0.1,fn=1 100 Hz,fm=107.3 Hz,r(t)為信噪比為7 dB的高斯白噪聲時,生成故障特征頻率為107.3 Hz的軸承外圈故障仿真信號。機床實際振動信號、軸承外圈故障仿真信號Y(t)、添加仿真故障后機床振動信號及其包絡(luò)譜如圖3所示。

    圖3 機床實際振動及仿真信號Fig.3 Actual vibration and simulation signal of machine tool

    由圖3(d)可知,添加仿真故障后機床振動信號中存在頻率為107.3 Hz的信號分量及其倍頻,與預(yù)期相符,證明故障仿真模塊的有效性。后續(xù)仿真其余軸承故障只需修改對應(yīng)軸承故障特征頻率即可。

    故障診斷模塊主要實現(xiàn)模型的訓練以及對測試數(shù)據(jù)的模式識別。利用Anaconda在PyCharm中配置好試驗所需環(huán)境后,使用keras以tensorflow為后端搭建一維卷積神經(jīng)網(wǎng)絡(luò)對數(shù)據(jù)集1進行學習,得到教師模型;用教師模型對數(shù)據(jù)集2進行識別,通過模型泛化能力識別正確的部分記為數(shù)據(jù)集2-1,識別為未知類別的部分記為數(shù)據(jù)集2-2;將數(shù)據(jù)集 2-2 和數(shù)據(jù)集1輸入網(wǎng)絡(luò),通過知識蒸餾和混合損失函數(shù)的方式對教師模型的知識進行遷移和更新得到學生模型,達到增量學習的目的;最終用學生模型對測試數(shù)據(jù)進行診斷分類,判斷機床所處狀態(tài)。模型中所使用軟件及版本為Anaconda3-5.3.1,PyCharm 2019.3.3、Keras2.3.1、Tensorflow 2.1,處理器為Intel Core i5-4210U。

    3 基于增量學習的故障診斷模型

    3.1 算法描述

    目前深度學習領(lǐng)域的卷積神經(jīng)網(wǎng)絡(luò)常采用堆疊式3×3的卷積核,既利于加深網(wǎng)絡(luò)深度,又可用較少參數(shù)獲取較大感受野,從而抑制過擬合[10]。以此為鑒并結(jié)合一維信號的特點[11],提出應(yīng)用于診斷機床故障的模型—基于增量學習的深度卷積診斷模型(deep convolutional neural network based on incremental learning,ILDCNN),模型第1層采用大尺度卷積核,作用是提取短時特征。該卷積核大小由多次試驗優(yōu)化所得,作用是學習與故障診斷相關(guān)的特征,防止因噪聲干擾產(chǎn)生的過度學習。其余卷積層均采用3×1的卷積核。ILDCNN模型結(jié)合2種不同尺度的卷積核,既保證運算速度,同時在加深網(wǎng)絡(luò)深度前提下抑制過擬合。

    此外,模型針對訓練時中間層數(shù)據(jù)分布改變的問題引入了Batch Normalization算法,在每1層卷積之后都采用批量歸一化操作,提高網(wǎng)絡(luò)的泛化能力的同時防止梯度消失或爆炸、加快訓練速度。模型還針對出現(xiàn)未知故障的問題加入未知故障診斷模塊,當識別到未知類別時,由人工初步判斷故障類型,若查明原因,則將未知樣本打標簽后重新加入模型進行增量學習,以此完善模型。增量學習[12]流程如下:首先將教師模型作為獨立網(wǎng)絡(luò)單獨訓練,保存訓練結(jié)果最優(yōu)時的網(wǎng)絡(luò)結(jié)構(gòu)和模型參數(shù),通過設(shè)置蒸餾溫度T輸出教師模型提供的標簽,記為軟標簽。其次,將數(shù)據(jù)輸入與教師模型結(jié)構(gòu)相同的ILDCNN模型(即學生模型)中,而后分成2步計算:1)使用與教師模型相同的蒸餾溫度T輸出ILDCNN模型的預(yù)測值,記為軟預(yù)測;2)直接輸出ILDCNN模型的預(yù)測值,記為硬預(yù)測。軟預(yù)測與軟標簽通過計算交叉熵得到蒸餾損失Lold,硬預(yù)測與數(shù)據(jù)真實標簽通過計算交叉熵得到蒸餾損失Lnew,Lold與Lnew相加得到總損失函數(shù)L,通過梯度下降更新ILDCNN模型的參數(shù),最終實現(xiàn)對教師模型知識的遷移和對新數(shù)據(jù)的學習。ILDCNN模型的訓練方式如圖4所示。

    圖4 ILDCNN模型訓練方式Fig.4 The training method of ILDCNN model

    3.2 BN層

    (3)

    (4)

    (5)

    yi=γxi+β

    (6)

    式中:xi為mini-batch中第i個輸入;m為mini-batch長度;ε為保證數(shù)值穩(wěn)定的常數(shù)項;γ為縮放參數(shù);β為偏移參數(shù)。

    3.3 未知故障診斷模塊

    分類問題中常用Softmax函數(shù)作為模型最后一層的激活函數(shù)用于解決多分類、單標簽問題。該函數(shù)如式(7)所示,其作用是將輸入均映射到(0,1)區(qū)間內(nèi),且所有映射后的輸出和為1,選取概率最大值為預(yù)測類別。但未知故障可能不屬于任一類別或?qū)儆诙鄠€類別,此為多分類、多標簽問題,因此采用適合該問題的Sigmoid函數(shù)作為模型最后一層的激活函數(shù)[15]。該函數(shù)如式(8)所示,以(0,0.5)中心對稱,滿足事件發(fā)生與否的概率,且任一輸入經(jīng)映射后輸出范圍均為(0,1),因此映射后的輸出經(jīng)閾值劃分后滿足可能不屬于任一類別或?qū)儆诙鄠€類別的情況,即未知故障的識別。

    (7)

    (8)

    式中:zi為輸入數(shù)據(jù)中第i個元素;n為總輸入個數(shù);x為輸入數(shù)據(jù)。

    閾值選取以發(fā)揮模型最佳性能為目的,經(jīng)多次試驗所得。加入此模塊前,任何未知故障均被誤認為已知故障中的一種,此模塊為未知類別的分類提供了依據(jù)。識別到未知故障后,人工初步判斷故障類型并為其打標簽,而后將其作為新數(shù)據(jù)重新送入模型學習其特征以完善模型,此階段為增量學習,訓練方式如圖4所示。

    3.4 增量學習

    實際工業(yè)現(xiàn)場中會出現(xiàn)機床轉(zhuǎn)速改變以及未知故障等問題。如何使ILDCNN模型對未知領(lǐng)域數(shù)據(jù)依舊有良好診斷能力,是衡量智能診斷算法好壞的重要指標,為此將引入增量學習算法。

    此算法運用知識蒸餾和混合損失函數(shù)的方式對教師模型進行更新。知識蒸餾通過定義蒸餾溫度θ使ILDCNN模型的軟預(yù)測匹配教師模型的軟輸出,從而繼承更多有用信息,實現(xiàn)知識遷移[12,16]。

    L=Lold+Lnew

    (9)

    (10)

    (11)

    (12)

    (13)

    (14)

    4 故障診斷試驗

    4.1 試驗數(shù)據(jù)準備

    試驗數(shù)據(jù)采集于立式加工中心vmc-850,為充分利用數(shù)據(jù)并避免過擬合,利用重疊采樣進行數(shù)據(jù)增強。以訓練樣本長度2 048為采樣窗口,偏移量為1進行采樣,采樣方式如圖5所示。

    試驗涉及8類機床故障,即報廢刀具和新刀具各自對應(yīng)的軸承無故障、外圈故障、內(nèi)圈故障和滾動體故障。具體分類標簽如表1所示。每一類均包含數(shù)據(jù)集1的900個訓練樣本和100個測試樣本、數(shù)據(jù)集2的4 500個訓練樣本和500個測試樣本。

    圖5 重疊采樣Fig.5 Overlap sampling

    4.2 試驗?zāi)P图皡?shù)

    教師模型和ILDCNN模型參數(shù)均如表2所示,蒸餾溫度θ設(shè)為2。

    表1數(shù)據(jù)集分類標簽Tab.1 The classification label of data set

    表2模型結(jié)構(gòu)參數(shù)Tab.2 The structure parameters of model

    4.3 教師模型訓練

    教師模型用數(shù)據(jù)集1訓練后曲線如圖6所示。教師模型準確率穩(wěn)定在90%左右,損失值穩(wěn)定在0.25左右,且準確率上升迅速,波動較小,證明教師模型在此轉(zhuǎn)速下具有良好的訓練效果。

    圖6 教師模型訓練曲線Fig.6 Teacher model training curve

    4.4 閾值選取

    教師模型若要用于測試,則需先設(shè)置合適閾值以識別未知類別。為確定最佳閾值,發(fā)揮模型最佳性能,選擇多個閾值對模型進行試驗。使用未知故障樣本對教師模型進行測試,模型將未知故障樣本正確識別為未知類別的概率為模型準確率;使用數(shù)據(jù)集1的測試數(shù)據(jù)對教師模型進行測試,模型將其錯誤識別為未知類別的概率為模型誤檢率。準確率和誤檢率隨閾值變化如圖7所示。設(shè)置0.7、0.75、0.8、0.85、0.9五個閾值對模型準確率和誤檢率進行檢測,準確率或誤檢率數(shù)值均為20次試驗測得的均值。

    圖7 準確率和誤檢率隨閾值變化圖Fig.7 The accuracy and false detection rate change with threshold

    由圖7可知,準確率處于67.98%~87.08%范圍內(nèi),隨閾值增大而增大;誤檢率處于2.99%~7.57%范圍內(nèi),隨閾值增大而增大,且閾值大于0.8后誤檢率迅速上升。為使模型發(fā)揮最佳性能,需使準確率盡可能高,誤檢率盡可能低,因此取最佳閾值為0.8。

    4.5 教師模型測試

    為了檢測教師模型對不同轉(zhuǎn)速下數(shù)據(jù)的故障診斷能力,使用教師模型對數(shù)據(jù)集2進行故障識別,閾值為0.8,診斷準確率如圖8所示。

    圖8 教師模型對數(shù)據(jù)集2的診斷準確率Fig.8 The diagnostic accuracy of the teacher model on data set 2

    由圖8可知,診斷準確率均處于52.31%~65.76%,模型可檢測部分故障,但未達到工業(yè)現(xiàn)場實際需求,因此引入增量學習,將未識別出的樣本(即數(shù)據(jù)集2-2)重新打標簽后送入教師模型訓練優(yōu)化,得到ILDCNN模型。

    4.6 ILDCNN模型訓練

    ILDCNN模型以教師模型為基礎(chǔ),對數(shù)據(jù)集2-2進行增量學習,ILDCNN模型的訓練曲線如圖9所示。

    圖9 ILDCNN模型訓練曲線Fig.9 ILDCNN model training curve

    圖9(a)中曲線1~5,準確率穩(wěn)步上升,經(jīng)過約20次訓練后,準確率最終穩(wěn)定在77.40%~86.08%之間。圖9(b)中曲線1~5損失值走向幾乎一致,經(jīng)過約20次訓練后,損失值最終均穩(wěn)定在0.6左右。由此說明模型收斂情況較好,且針對不同轉(zhuǎn)速的數(shù)據(jù),增量學習均有成效。

    4.7 ILDCNN模型測試

    為了驗證增量學習的效果及ILDCNN模型是否出現(xiàn)災(zāi)難性遺忘的問題,使用ILDCNN模型對數(shù)據(jù)集1和數(shù)據(jù)集2中的測試數(shù)據(jù)分別進行20次診斷試驗并記錄準確率,結(jié)合圖8中教師模型對數(shù)據(jù)集2的診斷準確率,得出增量學習前后診斷準確率對比圖如圖10所示。

    圖10 增量學習前后準確率對比圖Fig.10 Comparison chart of accuracy before and after incremental learning

    由圖10可知,增量學習后,ILDCNN模型對數(shù)據(jù)集2的診斷準確率均處于76.49%~86.09%之間,相比增量學習前,準確率提高約22%,模型診斷效果顯著提高,證明增量學習可提高模型跨轉(zhuǎn)速故障診斷的能力。ILDCNN模型對數(shù)據(jù)集1測試數(shù)據(jù)的診斷準確率均處于83.53%~90.54%之間,證明ILDCNN模型對教師模型知識的遷移效果良好,未出現(xiàn)災(zāi)難性遺忘。

    為了進一步驗證ILDCNN模型在跨轉(zhuǎn)速故障診斷領(lǐng)域的優(yōu)化效果,將其與對比模型進行對比試驗,對比模型包括Fine Tuning(針對跨領(lǐng)域數(shù)據(jù)進行微調(diào))和Joint Training(各領(lǐng)域數(shù)據(jù)聯(lián)合訓練)這2種經(jīng)典的跨領(lǐng)域診斷算法。3種模型均以數(shù)據(jù)集1訓練原始模型,數(shù)據(jù)集2用于增量學習。試驗結(jié)果如圖11所示。

    圖11 模型對比試驗圖Fig.11 Model comparison test chart

    由圖11可知,Fine Tuning總體準確率均較低且在數(shù)據(jù)集1上的表現(xiàn)極差,這是由于Fine Tuning只針對新領(lǐng)域數(shù)據(jù)(即數(shù)據(jù)集2)進行優(yōu)化;Joint Training在各個轉(zhuǎn)速下的準確率均處于80%左右,且對數(shù)據(jù)集1有較好診斷能力,這是由于Joint Training將數(shù)據(jù)集1、2的所有數(shù)據(jù)均用于優(yōu)化模型,但數(shù)據(jù)量巨大,因此訓練耗時長,且該方法優(yōu)化模型時需要所有原始數(shù)據(jù),既不方便也不符合大數(shù)據(jù)時代的發(fā)展;ILDCNN模型達到與Joint Training相似的準確率,但該模型僅需數(shù)據(jù)集2和部分數(shù)據(jù)集1,將訓練用時縮短近1/6,隨著模型優(yōu)化次數(shù)增多,其優(yōu)勢將遠遠大于Joint Training。

    對比ILDCNN模型與傳統(tǒng)時頻域診斷方法的優(yōu)劣,得出如下分析:時域分析法在機床振動信號有明顯特征時可識別故障,但由于隨機信號混雜,準確率往往較低。而軸承故障更是難以從時域中分析,通常使用頻域分析[17],因此時域分析法不適合機床故障診斷領(lǐng)域。頻域分析法在機床故障領(lǐng)域應(yīng)用廣泛,其診斷準確率與ILDCNN模型接近,但該方法僅適合已知故障類別的判斷,無法對復(fù)合故障等未知故障進行診斷[18]。綜上,機床故障診斷模型需滿足明確分類故障類別、可診斷未知故障等條件,而ILDCNN模型滿足上述要求,體現(xiàn)ILDCNN模型在故障診斷領(lǐng)域的優(yōu)勢。

    5 結(jié) 論

    ILDCNN模型針對跨轉(zhuǎn)速領(lǐng)域中無法識別的未知故障采用增量學習的方式,對比傳統(tǒng)方法做了以下優(yōu)化:1)利用Batch Normalization算法加速了訓練過程;2)利用增量學習實現(xiàn)對教師模型的知識遷移和對跨轉(zhuǎn)速數(shù)據(jù)的特征學習,使得ILDCNN模型在不同轉(zhuǎn)速下的準確率均處于76.49%~86.09%,相比增量學習前提高約22%。

    ILDCNN模型不僅各方面優(yōu)于傳統(tǒng)時頻域診斷方法,且與Fine Tuning和Joint Training兩種跨領(lǐng)域算法相比,提高了跨轉(zhuǎn)速故障診斷的準確率,縮短了訓練用時,為大數(shù)據(jù)下的故障診斷模型提供新思路。

    猜你喜歡
    增量故障診斷準確率
    提質(zhì)和增量之間的“辯證”
    當代陜西(2022年6期)2022-04-19 12:12:22
    乳腺超聲檢查診斷乳腺腫瘤的特異度及準確率分析
    健康之家(2021年19期)2021-05-23 11:17:39
    不同序列磁共振成像診斷脊柱損傷的臨床準確率比較探討
    2015—2017 年寧夏各天氣預(yù)報參考產(chǎn)品質(zhì)量檢驗分析
    “價增量減”型應(yīng)用題點撥
    高速公路車牌識別標識站準確率驗證法
    基于均衡增量近鄰查詢的位置隱私保護方法
    電信科學(2016年9期)2016-06-15 20:27:25
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    德州儀器(TI)發(fā)布了一對32位增量-累加模數(shù)轉(zhuǎn)換器(ADC):ADS1262和ADS126
    基于LCD和排列熵的滾動軸承故障診斷
    国产 精品1| 日本色播在线视频| 精品一品国产午夜福利视频| 日日摸夜夜添夜夜爱| 1000部很黄的大片| 久久精品熟女亚洲av麻豆精品| 校园人妻丝袜中文字幕| 2022亚洲国产成人精品| 国产伦精品一区二区三区四那| 狂野欧美激情性bbbbbb| 97在线视频观看| 亚洲av欧美aⅴ国产| 日韩视频在线欧美| 简卡轻食公司| 麻豆国产97在线/欧美| 久久婷婷青草| av卡一久久| av国产免费在线观看| 性高湖久久久久久久久免费观看| 中文字幕亚洲精品专区| 国精品久久久久久国模美| 免费播放大片免费观看视频在线观看| 黄色一级大片看看| 国产爽快片一区二区三区| 日本av手机在线免费观看| 最新中文字幕久久久久| 亚洲经典国产精华液单| 久久精品国产鲁丝片午夜精品| 18禁在线无遮挡免费观看视频| 久久99热这里只有精品18| 美女脱内裤让男人舔精品视频| 丝袜喷水一区| 下体分泌物呈黄色| 亚洲国产精品一区三区| 欧美日韩在线观看h| 日韩,欧美,国产一区二区三区| 亚洲人成网站在线播| 最后的刺客免费高清国语| 一本色道久久久久久精品综合| 欧美三级亚洲精品| 人妻 亚洲 视频| 欧美精品亚洲一区二区| 成人漫画全彩无遮挡| 久久国产乱子免费精品| 伦理电影大哥的女人| 久热久热在线精品观看| 亚洲国产最新在线播放| 亚洲国产欧美在线一区| 亚洲av国产av综合av卡| 国产伦理片在线播放av一区| 如何舔出高潮| 一区二区三区免费毛片| 777米奇影视久久| 2022亚洲国产成人精品| 少妇高潮的动态图| 丝袜喷水一区| 一级黄片播放器| 国产精品久久久久久av不卡| 日日撸夜夜添| 中文字幕av成人在线电影| 国产精品99久久久久久久久| 蜜桃久久精品国产亚洲av| 日本欧美国产在线视频| 亚洲精品国产色婷婷电影| av免费观看日本| 一级毛片黄色毛片免费观看视频| 男女无遮挡免费网站观看| 三级国产精品片| 一区二区三区乱码不卡18| 有码 亚洲区| av免费观看日本| 精品熟女少妇av免费看| kizo精华| 老熟女久久久| 久久精品夜色国产| 男男h啪啪无遮挡| 高清av免费在线| 精品久久久久久久久亚洲| 精品人妻偷拍中文字幕| 在线观看av片永久免费下载| 免费观看在线日韩| 国产在线免费精品| 80岁老熟妇乱子伦牲交| 亚洲精品国产av蜜桃| 国产爱豆传媒在线观看| 91精品一卡2卡3卡4卡| 熟妇人妻不卡中文字幕| 大香蕉97超碰在线| 国产视频首页在线观看| 搡女人真爽免费视频火全软件| 纵有疾风起免费观看全集完整版| 内射极品少妇av片p| 国产亚洲av片在线观看秒播厂| 99热这里只有是精品在线观看| 少妇裸体淫交视频免费看高清| 国产精品一区二区性色av| 国产精品不卡视频一区二区| 涩涩av久久男人的天堂| 99热这里只有是精品在线观看| 日本黄大片高清| 亚洲性久久影院| 国产成人a区在线观看| 女的被弄到高潮叫床怎么办| 国产人妻一区二区三区在| 91精品国产国语对白视频| 久久久国产一区二区| 男女下面进入的视频免费午夜| 男女啪啪激烈高潮av片| 成人亚洲欧美一区二区av| 91狼人影院| 日韩视频在线欧美| 直男gayav资源| 亚洲欧美日韩另类电影网站 | 中国国产av一级| 18禁裸乳无遮挡动漫免费视频| 亚洲成人中文字幕在线播放| 成人毛片a级毛片在线播放| .国产精品久久| 人人妻人人爽人人添夜夜欢视频 | 国产一区二区在线观看日韩| 日日啪夜夜撸| 亚洲av日韩在线播放| .国产精品久久| 一个人看视频在线观看www免费| 久久婷婷青草| 国产一区有黄有色的免费视频| 亚洲av不卡在线观看| 少妇熟女欧美另类| 国产老妇伦熟女老妇高清| 91精品国产国语对白视频| 久久综合国产亚洲精品| 多毛熟女@视频| 亚洲国产欧美在线一区| 亚洲av.av天堂| 国产精品一及| 色哟哟·www| av又黄又爽大尺度在线免费看| 免费少妇av软件| 最近的中文字幕免费完整| 18禁裸乳无遮挡免费网站照片| 日本av手机在线免费观看| 欧美国产精品一级二级三级 | freevideosex欧美| av在线蜜桃| 亚洲成人av在线免费| 免费看光身美女| 亚洲精华国产精华液的使用体验| 国产成人精品一,二区| 亚洲av在线观看美女高潮| 国产视频内射| 国产成人aa在线观看| 99久久精品国产国产毛片| 毛片一级片免费看久久久久| 一级二级三级毛片免费看| 亚洲内射少妇av| 黑人高潮一二区| 精品一区二区免费观看| 国产真实伦视频高清在线观看| 日韩一本色道免费dvd| 亚洲电影在线观看av| 人人妻人人澡人人爽人人夜夜| 国产精品99久久久久久久久| 国产精品国产av在线观看| 久久久久久久亚洲中文字幕| 国产深夜福利视频在线观看| 午夜免费男女啪啪视频观看| 性高湖久久久久久久久免费观看| 亚洲电影在线观看av| 国产精品成人在线| 婷婷色麻豆天堂久久| 国产精品一区二区性色av| 国产在线免费精品| 国产探花极品一区二区| 中国美白少妇内射xxxbb| 久热这里只有精品99| 黄色欧美视频在线观看| 欧美高清性xxxxhd video| 51国产日韩欧美| 亚洲欧美清纯卡通| 国产成人a∨麻豆精品| 一级毛片久久久久久久久女| 一本—道久久a久久精品蜜桃钙片| 我的老师免费观看完整版| 亚洲成人av在线免费| 内射极品少妇av片p| 精品国产乱码久久久久久小说| 纵有疾风起免费观看全集完整版| 乱系列少妇在线播放| 在线观看一区二区三区激情| 夜夜看夜夜爽夜夜摸| 国产成人精品婷婷| 亚洲人与动物交配视频| 联通29元200g的流量卡| 纵有疾风起免费观看全集完整版| 亚洲av福利一区| 2021少妇久久久久久久久久久| 99热这里只有是精品50| 制服丝袜香蕉在线| 精品一区在线观看国产| 免费在线观看成人毛片| 蜜桃在线观看..| 久久久亚洲精品成人影院| 亚洲自偷自拍三级| 亚洲婷婷狠狠爱综合网| 妹子高潮喷水视频| av视频免费观看在线观看| 晚上一个人看的免费电影| av不卡在线播放| 亚洲av免费高清在线观看| 国产精品秋霞免费鲁丝片| 日韩中字成人| 免费观看a级毛片全部| 亚洲成人中文字幕在线播放| 免费高清在线观看视频在线观看| 亚洲av中文av极速乱| 亚洲av日韩在线播放| 成人毛片60女人毛片免费| 91久久精品电影网| 日产精品乱码卡一卡2卡三| 国产黄色免费在线视频| 成人毛片a级毛片在线播放| 国产免费一级a男人的天堂| 在线看a的网站| 永久免费av网站大全| 亚洲精品乱码久久久久久按摩| 我的老师免费观看完整版| 精品熟女少妇av免费看| 亚洲国产高清在线一区二区三| 1000部很黄的大片| 久久久亚洲精品成人影院| 亚洲欧美日韩无卡精品| 成人国产av品久久久| 最近手机中文字幕大全| 夜夜骑夜夜射夜夜干| 国产片特级美女逼逼视频| 久久精品夜色国产| 成人二区视频| 三级国产精品片| 精品99又大又爽又粗少妇毛片| 亚洲国产精品国产精品| 久久精品熟女亚洲av麻豆精品| 日韩欧美一区视频在线观看 | 在线亚洲精品国产二区图片欧美 | freevideosex欧美| 国产 一区精品| 亚洲av成人精品一二三区| 中国三级夫妇交换| xxx大片免费视频| 91久久精品国产一区二区三区| 国产免费一区二区三区四区乱码| 免费不卡的大黄色大毛片视频在线观看| 国产亚洲91精品色在线| a级毛片免费高清观看在线播放| 超碰av人人做人人爽久久| 天美传媒精品一区二区| 亚洲,一卡二卡三卡| 高清欧美精品videossex| 久久精品国产亚洲网站| 久久久久久久大尺度免费视频| 国产精品久久久久久久电影| 欧美三级亚洲精品| 交换朋友夫妻互换小说| 免费黄网站久久成人精品| 欧美成人午夜免费资源| 黄色欧美视频在线观看| 日日摸夜夜添夜夜添av毛片| 不卡视频在线观看欧美| 啦啦啦视频在线资源免费观看| 亚洲欧美清纯卡通| 亚洲国产精品成人久久小说| 免费看日本二区| videos熟女内射| 夜夜爽夜夜爽视频| 一级爰片在线观看| 日韩在线高清观看一区二区三区| 国产欧美亚洲国产| 纵有疾风起免费观看全集完整版| 亚洲一区二区三区欧美精品| 国产黄频视频在线观看| 亚洲伊人久久精品综合| 免费久久久久久久精品成人欧美视频 | 日本av免费视频播放| av免费在线看不卡| 国产精品免费大片| 美女xxoo啪啪120秒动态图| 美女内射精品一级片tv| 国产 精品1| 高清毛片免费看| 国产色婷婷99| 91精品一卡2卡3卡4卡| 黑丝袜美女国产一区| 亚洲精品乱码久久久久久按摩| 日本wwww免费看| 一边亲一边摸免费视频| 亚洲精品国产色婷婷电影| 亚洲国产最新在线播放| 免费在线观看成人毛片| 极品少妇高潮喷水抽搐| 尾随美女入室| av福利片在线观看| 我的老师免费观看完整版| av国产免费在线观看| 一级片'在线观看视频| 日本爱情动作片www.在线观看| 久久午夜福利片| 久久久久性生活片| 乱码一卡2卡4卡精品| 99国产精品免费福利视频| 亚洲av二区三区四区| 成人毛片60女人毛片免费| 99热国产这里只有精品6| 少妇的逼好多水| 女人十人毛片免费观看3o分钟| 日韩亚洲欧美综合| 国产亚洲一区二区精品| 丝袜喷水一区| 久久鲁丝午夜福利片| 久久精品人妻少妇| 视频区图区小说| 国产男女内射视频| 最后的刺客免费高清国语| 免费看av在线观看网站| 大又大粗又爽又黄少妇毛片口| 人人妻人人爽人人添夜夜欢视频 | 亚洲欧美一区二区三区国产| 美女视频免费永久观看网站| 国产精品麻豆人妻色哟哟久久| 色综合色国产| 只有这里有精品99| 国产欧美另类精品又又久久亚洲欧美| 成人影院久久| 久久国内精品自在自线图片| 久久久久久久大尺度免费视频| 国产午夜精品一二区理论片| 美女脱内裤让男人舔精品视频| 午夜日本视频在线| 各种免费的搞黄视频| 一区二区三区四区激情视频| 亚洲欧美日韩另类电影网站 | av专区在线播放| 国产黄片美女视频| 一区二区三区精品91| 国产男人的电影天堂91| 午夜福利视频精品| 久久久久国产精品人妻一区二区| 美女视频免费永久观看网站| 18+在线观看网站| 边亲边吃奶的免费视频| 日韩制服骚丝袜av| 久久久久人妻精品一区果冻| 午夜福利影视在线免费观看| 精品99又大又爽又粗少妇毛片| 日韩成人伦理影院| 欧美精品人与动牲交sv欧美| 久久久久精品久久久久真实原创| 我要看日韩黄色一级片| 老司机影院成人| 乱码一卡2卡4卡精品| 国产成人freesex在线| 日韩欧美 国产精品| 在线观看三级黄色| 久久韩国三级中文字幕| 99热这里只有是精品在线观看| 丝袜喷水一区| 九草在线视频观看| 久久精品国产自在天天线| 欧美性感艳星| 简卡轻食公司| 涩涩av久久男人的天堂| 三级经典国产精品| 最黄视频免费看| 中国国产av一级| 一级a做视频免费观看| 精品久久久精品久久久| 国产精品伦人一区二区| 国产欧美日韩精品一区二区| 午夜福利在线观看免费完整高清在| 国产成人91sexporn| 91久久精品国产一区二区成人| 少妇人妻精品综合一区二区| 国产永久视频网站| 国产白丝娇喘喷水9色精品| 久久久午夜欧美精品| 一边亲一边摸免费视频| 国国产精品蜜臀av免费| 六月丁香七月| av不卡在线播放| 蜜桃久久精品国产亚洲av| 伊人久久国产一区二区| 欧美激情国产日韩精品一区| 亚洲欧美精品专区久久| 日韩强制内射视频| 黄色视频在线播放观看不卡| 久久久精品免费免费高清| 午夜福利高清视频| 五月玫瑰六月丁香| 日韩欧美精品免费久久| 97在线人人人人妻| 国产精品久久久久久久电影| 日本wwww免费看| 久久 成人 亚洲| 久久影院123| 亚洲国产最新在线播放| 久久精品国产a三级三级三级| 国产成人精品福利久久| 久久精品国产自在天天线| 欧美日韩视频高清一区二区三区二| 国产精品99久久99久久久不卡 | 久久国内精品自在自线图片| 直男gayav资源| av天堂中文字幕网| 日韩制服骚丝袜av| 一级毛片我不卡| 国产一区有黄有色的免费视频| 亚洲国产高清在线一区二区三| 久久久久视频综合| 国内精品宾馆在线| 国产大屁股一区二区在线视频| 久久人妻熟女aⅴ| 五月玫瑰六月丁香| 日本欧美国产在线视频| av黄色大香蕉| 亚洲第一区二区三区不卡| 大片电影免费在线观看免费| 亚洲国产色片| 日本欧美视频一区| 国产成人精品久久久久久| 各种免费的搞黄视频| 男人和女人高潮做爰伦理| 亚洲av在线观看美女高潮| 伦理电影大哥的女人| 亚洲一级一片aⅴ在线观看| 日本欧美国产在线视频| 精品一区二区三卡| 寂寞人妻少妇视频99o| 男女边吃奶边做爰视频| 免费在线观看成人毛片| 亚洲性久久影院| 日本与韩国留学比较| 久久韩国三级中文字幕| 国产日韩欧美亚洲二区| 寂寞人妻少妇视频99o| 亚洲自偷自拍三级| 免费黄色在线免费观看| 久久亚洲国产成人精品v| 国产一区二区三区av在线| 天堂8中文在线网| 99热这里只有精品一区| 一级毛片黄色毛片免费观看视频| 亚洲自偷自拍三级| 热99国产精品久久久久久7| 国产午夜精品一二区理论片| 中文欧美无线码| 免费观看无遮挡的男女| 久久久久久久国产电影| 久久人妻熟女aⅴ| 欧美国产精品一级二级三级 | 日韩av在线免费看完整版不卡| 亚洲精品久久久久久婷婷小说| 国产高清国产精品国产三级 | 日韩av不卡免费在线播放| 噜噜噜噜噜久久久久久91| 国产乱来视频区| 高清视频免费观看一区二区| 久久精品久久久久久噜噜老黄| 亚洲综合精品二区| 国产色婷婷99| 男人舔奶头视频| 国产免费视频播放在线视频| 国产精品麻豆人妻色哟哟久久| 两个人的视频大全免费| 久久久久精品性色| 久久毛片免费看一区二区三区| 亚洲美女黄色视频免费看| 91久久精品国产一区二区成人| 精品人妻一区二区三区麻豆| freevideosex欧美| 女人十人毛片免费观看3o分钟| 99久久人妻综合| 你懂的网址亚洲精品在线观看| 激情 狠狠 欧美| 国产深夜福利视频在线观看| 欧美日韩一区二区视频在线观看视频在线| 自拍欧美九色日韩亚洲蝌蚪91 | 91在线精品国自产拍蜜月| 久久久久久久久久久丰满| av网站免费在线观看视频| 一级二级三级毛片免费看| 久热久热在线精品观看| 色网站视频免费| 黑人猛操日本美女一级片| 香蕉精品网在线| 欧美区成人在线视频| 男女国产视频网站| 99热这里只有是精品在线观看| 久久久久久久亚洲中文字幕| 成人亚洲欧美一区二区av| 欧美日韩视频精品一区| 午夜激情久久久久久久| 欧美xxxx黑人xx丫x性爽| 亚洲精品中文字幕在线视频 | 天堂中文最新版在线下载| 国产欧美另类精品又又久久亚洲欧美| 婷婷色综合www| 国产在线男女| 最新中文字幕久久久久| 51国产日韩欧美| 亚洲av二区三区四区| 精品久久国产蜜桃| 99久久人妻综合| 亚洲av中文av极速乱| 高清不卡的av网站| 久久精品国产a三级三级三级| 久久久久国产网址| av国产免费在线观看| 国产免费一区二区三区四区乱码| 日韩,欧美,国产一区二区三区| 青青草视频在线视频观看| 国产在线一区二区三区精| 国产在线一区二区三区精| 国产精品免费大片| 日韩av在线免费看完整版不卡| 午夜老司机福利剧场| 免费大片黄手机在线观看| 丰满迷人的少妇在线观看| 中文在线观看免费www的网站| 26uuu在线亚洲综合色| 精品人妻熟女av久视频| 国产乱人偷精品视频| 99热国产这里只有精品6| 国产乱来视频区| 天天躁夜夜躁狠狠久久av| 国产成人a∨麻豆精品| 美女脱内裤让男人舔精品视频| 日本欧美国产在线视频| 国产有黄有色有爽视频| 精品久久久久久久末码| 中文字幕制服av| 一级毛片aaaaaa免费看小| 一本色道久久久久久精品综合| 国产色爽女视频免费观看| 亚州av有码| 男人添女人高潮全过程视频| 日日撸夜夜添| 大又大粗又爽又黄少妇毛片口| 一级a做视频免费观看| 成人一区二区视频在线观看| 日韩av免费高清视频| 国产伦精品一区二区三区四那| 多毛熟女@视频| 久久毛片免费看一区二区三区| 久久ye,这里只有精品| 少妇的逼水好多| 国产国拍精品亚洲av在线观看| 99久久人妻综合| 2018国产大陆天天弄谢| 免费看不卡的av| 久久久精品94久久精品| 中文资源天堂在线| 99九九线精品视频在线观看视频| 国产精品99久久久久久久久| 国精品久久久久久国模美| 在线天堂最新版资源| 2021少妇久久久久久久久久久| 九草在线视频观看| 美女cb高潮喷水在线观看| av.在线天堂| 精品一区二区三区视频在线| 这个男人来自地球电影免费观看 | 国产成人a区在线观看| 亚洲怡红院男人天堂| 午夜精品国产一区二区电影| 黄色欧美视频在线观看| 一区二区av电影网| 日日啪夜夜爽| 亚洲一级一片aⅴ在线观看| 丝瓜视频免费看黄片| 亚洲一级一片aⅴ在线观看| 欧美极品一区二区三区四区| 人人妻人人看人人澡| a级毛片免费高清观看在线播放| 亚洲精品一二三| 观看美女的网站| 国产精品一区二区在线观看99| 五月玫瑰六月丁香| 少妇人妻精品综合一区二区| 美女主播在线视频| 国产视频内射| 精品国产三级普通话版| 亚洲av成人精品一区久久| 国产高清国产精品国产三级 | 午夜日本视频在线| 国产 精品1| 成人黄色视频免费在线看| a 毛片基地| 蜜臀久久99精品久久宅男| 精品99又大又爽又粗少妇毛片| 国产成人一区二区在线| 女人十人毛片免费观看3o分钟| 欧美日韩一区二区视频在线观看视频在线| 在线 av 中文字幕| 中文字幕av成人在线电影| 熟妇人妻不卡中文字幕| 18禁裸乳无遮挡动漫免费视频| 亚洲第一av免费看| 国产亚洲最大av| 在线播放无遮挡| 少妇裸体淫交视频免费看高清| 伦理电影大哥的女人| 国产成人a∨麻豆精品| 高清午夜精品一区二区三区| 又粗又硬又长又爽又黄的视频| 99热这里只有是精品在线观看| 亚洲激情五月婷婷啪啪| 女的被弄到高潮叫床怎么办|