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

    基于EMD的金剛石砂輪磨損狀態(tài)聲發(fā)射監(jiān)測

    2019-04-13 13:23郭力霍可可郭君濤

    郭力 霍可可 郭君濤

    摘 ???要:針對磨削金剛石砂輪磨損狀態(tài)聲發(fā)射信號小波分析中存在的問題,根據(jù)工程陶瓷部分穩(wěn)定氧化鋯磨削過程中聲發(fā)射信號非線性非平穩(wěn)性的特點(diǎn),采用經(jīng)驗?zāi)B(tài)分解方法將磨削聲發(fā)射信號分解為多個平穩(wěn)的固有模態(tài)函數(shù)之和,并提取其有效值、方差和能量系數(shù)等特征值.在磨削金剛石砂輪從輕度磨損狀態(tài)轉(zhuǎn)變?yōu)閲?yán)重磨損狀態(tài)時,固有模態(tài)函數(shù)的有效值(IMFrms)和方差(IMFvar)增大,而能量系數(shù)(IMFpe)發(fā)生明顯的變化;將其做為最小二乘支持向量機(jī)的輸入?yún)?shù),對金剛石砂輪的輕度磨損狀態(tài)和嚴(yán)重磨損狀態(tài)成功地進(jìn)行了智能監(jiān)測.

    關(guān)鍵詞:氧化鋯磨削;金剛石砂輪磨損狀態(tài)監(jiān)測;聲發(fā)射;經(jīng)驗?zāi)B(tài)分解;最小二乘支持向量機(jī)

    中圖分類號:TG58,TP183 ??????????????????????文獻(xiàn)標(biāo)志碼:A

    Acoustic Emission Intelligent Monitoring of Diamond Grinding

    Wheel Wear Based on Empirical Mode Decomposition

    GUO Li?覮 ,HUO Keke,GUO Juntao

    (College of Mechanical and Vehicle Engineering, Hunan University, Changsha 410082,China)

    Abstract: In view of the existing problem in the wavelet analysis of acoustic emission signals in wear state of diamond grinding wheel, because engineering ceramics partially stabilized zirconia grinding acoustic emission signals have nonlinear and nonstationary characteristics, using empirical mode decomposition method the acoustic emission signals were decomposed into several stationary intrinsic mode functions and then the root mean squares, variances and energy coefficients were extracted. When the wear state of diamond grinding wheel changes from mild wear to severe wear, the root mean squares(IMFrms) and variances(IMFvar) of the intrinsic mode function increase, and the energy coefficients(IMFpe) change significantly. As the input parameter of the least squares support vector machine, the wear state of diamond grinding wheel was successfully monitored.

    Key words: partially stabilized zirconia grinding;diamond grinding wheel wear state monitoring;acoustic emission;empirical mode decomposition;least squares support vector machine

    工程陶瓷部分穩(wěn)定氧化鋯PSZ以高耐磨性、耐高溫、高硬度、耐腐蝕、絕緣和不導(dǎo)熱等優(yōu)異性能,在航空航天和電子、化學(xué)、機(jī)械工業(yè)等重要領(lǐng)域的應(yīng)用越來越廣泛.一般工程陶瓷要精密加工后才能使用,要用金剛石砂輪來磨削工程陶瓷,由于工程陶瓷是硬脆難加工材料,使昂貴的金剛石砂輪磨削中容易產(chǎn)生磨損導(dǎo)致磨削力磨削溫度增大,磨削去除率和磨削質(zhì)量下降,導(dǎo)致磨削加工費(fèi)用上升.為了準(zhǔn)確地預(yù)測金剛石砂輪的磨損狀態(tài),以便及時修整砂輪,提高金剛石砂輪的利用率,降低金剛石砂輪的損耗,提高工程陶瓷的磨削質(zhì)量,有必要在陶瓷磨削加工中對金剛石砂輪的磨損狀態(tài)進(jìn)行監(jiān)測.金剛石砂輪磨損狀態(tài)測量難以采用直接方法,即便用直接法測量會導(dǎo)致磨削過程中斷.而聲發(fā)射(acoustic emission,AE)等非接觸間接監(jiān)測砂輪磨損減少了停機(jī)時間和成本,提高砂輪壽命和磨削質(zhì)量,又不依賴磨床操作者的技能,是智能制造監(jiān)測的主要的研究方向;它包括磨削AE信號采集與處理、特征提取與智能決策.

    文獻(xiàn)[1]應(yīng)用經(jīng)驗?zāi)B(tài)分解(Empirical Mode Decomposition, EMD)將車削加工聲發(fā)射信號分解為幾個固有模態(tài)函數(shù)(Intrinsic Mode Function, IMF),再提取IMF分量的歸一化能量值并將其作為特征向量,應(yīng)用最小二乘支持向量機(jī)LS-SVM對車刀磨損狀態(tài)成功地進(jìn)行識別.文獻(xiàn)[2]對陶瓷磨削砂輪磨損進(jìn)行了聲發(fā)射智能監(jiān)測,應(yīng)用了自適應(yīng)神經(jīng)模糊推理系統(tǒng)、多層感知器神經(jīng)網(wǎng)絡(luò)、徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)和廣義回歸神經(jīng)網(wǎng)絡(luò)等幾種先進(jìn)的方法,砂輪小切深條件下砂輪磨損分類準(zhǔn)確率達(dá)90%以上.文獻(xiàn)[3]應(yīng)用總體經(jīng)驗?zāi)B(tài)分解 EEMD得到無模態(tài)混疊的固有模態(tài)函數(shù)IMF分量,再應(yīng)用相應(yīng)的算法從所有IMF分量中提取敏感的IMF;提取敏感IMF的能量作為支持向量機(jī)(SVM)分類器的輸入,車刀磨損狀態(tài)分類預(yù)測準(zhǔn)確率在93%以上.文獻(xiàn)[4]研究了基于聲發(fā)射信號EMD分解與IMF奇異值熵的車刀磨損狀態(tài)監(jiān)測方法,由于鋒利車刀的IMF奇異值熵大于磨損車刀的,所以可以用來監(jiān)測車刀的磨損.文獻(xiàn)[5]把經(jīng)驗?zāi)B(tài)分解EMD應(yīng)用到磨削燒傷聲發(fā)射信號的分析中,提出了基于平均能量的固有模態(tài)函數(shù)IMF分量選擇方法.文獻(xiàn)[6]認(rèn)為希爾伯特黃變換(HHT)能夠自適應(yīng)地高效分解信號,反應(yīng)信號的局部特征,可以分析非線性非平穩(wěn)振動信號.所提取的銑刀磨損特征固有模態(tài)函數(shù)IMF3、IMF4、IMF5和IMF6的振幅均值和邊際譜中的最大幅值能準(zhǔn)確可靠地反應(yīng)銑刀磨損狀態(tài).文獻(xiàn)[7]切削聲發(fā)射信號的頻率高,而切削加工中的環(huán)境噪聲信號頻率都比較低,所以聲發(fā)射信號的抗干擾能力強(qiáng).切削刀具磨損的聲發(fā)射信號具有非線性特性,所以采用混沌理論對不同刀具磨損程度的聲發(fā)射信號進(jìn)行了分析,分析發(fā)現(xiàn)關(guān)聯(lián)維數(shù)可作為特征量用于刀具磨損檢測.文獻(xiàn)[8]發(fā)現(xiàn)可以通過使用激光照射試驗產(chǎn)生的熱聲發(fā)射信號來訓(xùn)練神經(jīng)網(wǎng)絡(luò),對磨削高溫合金的砂輪磨損狀態(tài)進(jìn)行監(jiān)測,有效地減少了磨削實(shí)驗而節(jié)省了費(fèi)用.文獻(xiàn)[9]認(rèn)為自適應(yīng)最稀疏時頻分析方法在抗噪聲性能、抑制端點(diǎn)效應(yīng)和模態(tài)混淆以及提高分量的準(zhǔn)確性等方面要優(yōu)于經(jīng)驗?zāi)B(tài)分解方法,是信號分析的好方法.文獻(xiàn)[10]提出了兩種蟻群優(yōu)化算法和著名序列前向浮動選擇法,作為聲發(fā)射信號特征提取的最優(yōu)化方法,作為磨削砂輪磨損狀態(tài)智能監(jiān)測.

    由于小波分解方法只能分析平穩(wěn)信號,所以使用小波分解方法來準(zhǔn)確分析金剛石砂輪磨損聲發(fā)射信號這種非平穩(wěn)信號是比較困難的[1].本文結(jié)合工程陶瓷部分穩(wěn)定氧化鋯PSZ磨削過程中AE信號非線性非平穩(wěn)性的特點(diǎn),采用經(jīng)驗?zāi)B(tài)分解EMD方法將其AE信號分解為多個平穩(wěn)的固有模態(tài)函數(shù)之和,再提取其有效值、方差和能量系數(shù)等特征值,將其做為最小二乘支持向量機(jī)LS-SVM的輸入?yún)?shù)對磨削金剛石砂輪的輕度磨損和嚴(yán)重磨損狀態(tài)進(jìn)行智能監(jiān)測.

    1 ??磨削聲發(fā)射信號的采集

    在一臺MGS-250AH精密平面磨床上進(jìn)行試驗,試驗中使用青銅結(jié)合劑金剛石砂輪粒度為60#、寬度為10 mm、直徑為150 mm. 試驗中工作臺速度不變,砂輪轉(zhuǎn)速固定為1 600 r/min即線速度12 m/s,磨削保持切深為5 μm,即實(shí)驗中磨削參數(shù)保持不變.工程陶瓷部分穩(wěn)定氧化鋯PSZ磨削試件安裝在虎鉗夾具中,夾具用電磁吸力安裝在磨床工作臺上,聲發(fā)射傳感器用磁力安裝在虎鉗工件夾具上盡可能靠近磨削工件.試驗不用磨削液,消除了磨削液對磨削聲發(fā)射信號分析的影響.本試驗采用的部分穩(wěn)定氧化鋯PSZ磨削試件尺寸為35 mm×15 mm×10 mm, 磨削在試件35 mm×15 mm的平面內(nèi)沿15 mm方向進(jìn)行.PSZ材料密度為5.97 g/cm3; 抗彎強(qiáng)度大于800 MPa; 微觀硬度為12.5 GPa;斷裂韌性為7.5 MPa·m1/2;彈性模量為205 GPa;PSZ材料掃描電鏡照片顯示其晶粒尺寸小于1 μm,其中存在一定量的氣孔.可見PSZ工程陶瓷是典型的硬脆難加工材料.

    PSZ磨削過程產(chǎn)生的AE信號,傳播到放在平面磨床工作臺夾具上的美國物理聲學(xué)公司PAC生產(chǎn)的壓電AE傳感器WSa,經(jīng)過AE壓電連接器的前置放大,再到PAC生產(chǎn)的PCI-2聲發(fā)射信號采集卡和AEwin采集系統(tǒng)進(jìn)行信號采集,并輸入到計算機(jī)進(jìn)行保存等待MATLAB軟件進(jìn)行后續(xù)的處理.為了不失真,AE信號采樣頻率為2 MHz,前置放大增益為40 dB,門檻為固定40 dB.WSa聲發(fā)射傳感器的工作頻率范圍在100 ~ 1 000 kHz,諧振頻率為125 kHz.

    采集金剛石砂輪在輕度磨損和嚴(yán)重磨損[10-11]程度下的多次磨削PSZ聲發(fā)射信號,每次實(shí)驗?zāi)ハ鲄?shù)相同,即實(shí)驗中工作臺速度不變,砂輪轉(zhuǎn)速固定為1 600 r·min-1即線速度12 m/s,磨削保持切深為5 μm.砂輪輕度磨損實(shí)驗是采用修整后的金剛石砂輪磨削5次PSZ,此時砂輪磨粒磨刃鋒銳,所以用測力儀測量到的法向磨削力增加很小,因而采集到的磨削AE信號可作為砂輪輕度磨損的AE信號.砂輪嚴(yán)重磨損實(shí)驗是在砂輪輕度磨損的5次PSZ磨削實(shí)驗結(jié)束后繼續(xù)進(jìn)行,一直磨削PSZ到20次,砂輪隨著磨削次數(shù)增加磨損加劇,所以用測力儀測量到的法向磨削力增加明顯,采集到的磨削AE信號作為砂輪嚴(yán)重磨損的AE信號.因此實(shí)驗中每磨削20次對金剛石砂輪進(jìn)行修整,再進(jìn)行下次砂輪輕度磨損和嚴(yán)重磨損磨削AE實(shí)驗,實(shí)驗中只測量砂輪磨損的磨削AE信號.

    2 ??磨削聲發(fā)射信號的處理

    圖1~圖2可以看出金剛石砂輪在輕度磨損變到嚴(yán)重磨損時磨削AE原始信號幅值增加,但是AE信號頻率主要分布范圍沒有發(fā)生明顯的變化,都是主要集中在0~600 kHz之內(nèi).AE信號在頻率峰值上有稍微的變化:砂輪輕度磨損時主要峰值頻率29 kHz、95 kHz和201 kHz;砂輪嚴(yán)重磨損時主要峰值頻率30 kHz、88 kHz和220 kHz.金剛石砂輪表面狀態(tài)對磨削過程有重要的影響,砂輪磨損后磨粒的切削刃鈍化,磨削力加大磨削溫度升高.

    經(jīng)驗?zāi)B(tài)分解EMD[1]方法是將非線性非平穩(wěn)的砂輪磨損聲發(fā)射信號自適應(yīng)分解成一系列平穩(wěn)的固有模態(tài)函數(shù)IMF,所分解出來的各IMF分量包含了原聲發(fā)射信號的不同時間尺度的局部特征信號,并且無須預(yù)先設(shè)定任何基函數(shù),優(yōu)于傳統(tǒng)的傅里葉變換和小波變換聲發(fā)射信號分析方法.在磨削金剛石砂輪與陶瓷工件的聲發(fā)射信號較低頻率段中,一般是磨削液和機(jī)床運(yùn)動部件摩擦產(chǎn)生的噪聲;這樣就可以把分布在較高頻率段的金剛石砂輪磨粒與陶瓷工件的切削、滑擦和犁耕以及切屑斷裂、磨粒破裂等產(chǎn)生的聲發(fā)射信號[12-13]清晰地提取出來進(jìn)行分析.

    采集到的PSZ磨削實(shí)驗聲發(fā)射波形流信號轉(zhuǎn)化為文本文件,并導(dǎo)入MATLAB軟件中的EMD工具箱將采集到的聲發(fā)射信號進(jìn)行EMD分解,每個聲發(fā)射信號可以分解為12~16個固有模態(tài)函數(shù)(IMF)和一個剩余信號.固有模態(tài)函數(shù)和剩余信號與原始信號的相關(guān)系數(shù)很大程度上可以反映它們之間的相關(guān)程度,相關(guān)系數(shù)越大說明IMF含原信號中的有效成分越高[1].PSZ磨削聲發(fā)射信號經(jīng)驗?zāi)B(tài)分解信號的相關(guān)系數(shù)如圖3所示,可以看出imf1~imf6與原始信號的相關(guān)系數(shù)較大,其中又以imf4、imf5、imf3的相關(guān)系數(shù)最大;imf7之后固有模態(tài)函數(shù)和剩余信號與原始信號的相關(guān)系數(shù)接近零即可以剔除.因此固有模態(tài)函數(shù)imf1~imf6包含了原始信號中主要的有效成分,可以作為原始信號的代表信號,從中可以提取信息用于砂輪磨損狀態(tài)的監(jiān)測.

    砂輪不同磨損狀態(tài)時的聲發(fā)射信號經(jīng)過經(jīng)驗?zāi)B(tài)分解得到的imf1~imf6信號的波形和頻譜,如圖4和圖5.從頻譜圖5可見,不同固有模態(tài)函數(shù)的頻率范圍有明顯的區(qū)別:imf1信號的主要頻率在550 kHz附近;imf2的主要頻率在210 kHz附近;imf3的主要頻率在100 kHz附近;imf4的主要頻率在50

    kHz附近;imf5 和imf6的主要頻率在30 kHz附近.砂輪輕度磨損與嚴(yán)重磨損時的聲發(fā)射信號經(jīng)過EMD處理的各階固有模態(tài)函數(shù)頻率主要分布范圍和主要的峰值頻率上沒有明顯的變化;但是從表1中看出,砂輪嚴(yán)重磨損時的波形幅值明顯大于砂輪輕度磨損時的波形幅值.隨著砂輪磨損加劇,固有模態(tài)函數(shù)imf1~imf6的幅值增大,也就是各個分量的能量出現(xiàn)了不同程度的增大,其中高頻分量imf1~imf3對應(yīng)幅值增大明顯在兩倍以上;而低頻分量imf4~imf6的增加較小.說明隨著金剛石砂輪磨損程度的加深,AE信號高頻成分在不斷增加.砂輪在輕度磨損到嚴(yán)重磨損imf1~imf6的幅值都有明顯增大,所以它們都是砂輪磨損的特征量.

    原始磨削聲發(fā)射信號三個峰值頻率分別在30 kHz、100 kHz和210 kHz附近,這與經(jīng)驗?zāi)B(tài)分解得到的固有模態(tài)函數(shù)imf2~imf6的主要頻率(如表1)相似,因而可以說imf2~imf6中分別包含了原始聲發(fā)射信號中的一個重要頻率段的主要信息,imf2包含原信號中220 kHz頻率附近的大量信息,imf3包含原信號中100 kHz頻率附近的大量信息,imf4~imf6包含原信號中30 kHz頻率附近的大量信息,而imf1則包含了400 kHz以上的高頻部分的主要信息,因而固有模態(tài)函數(shù)imf1~imf6與原始AE信號密切相關(guān),所以它們的相關(guān)系數(shù)較大.

    分析還可見固有模態(tài)函數(shù)imf1~imf6的能量占據(jù)了AE信號總能量的很大比重,imf1~imf6信號的總能量占據(jù)了原始AE信號總能量的全部(98%以上),其中imf3~imf5(頻率主要在100 kHz以下)起主要作用,占據(jù)了原始AE信號總能量的85%以上,可以看出磨削氧化鋯PSZ陶瓷的聲發(fā)射信號頻率主要集中在100 kHz以內(nèi).同時由圖3可見imf3~imf5的相關(guān)系數(shù)最大,所以相關(guān)系數(shù)與聲發(fā)射信號能量密切相關(guān).

    3 ??砂輪磨損聲發(fā)射信號特征提取

    對EMD分解得到的固有模態(tài)函數(shù)imf1~imf6進(jìn)行分析,利用式(1)~(3)計算每個IMF的有效值(IMFrms)、方差(IMFvar)、能量系數(shù)(IMFpe).金剛石砂輪的磨損加劇,磨削AE信號增強(qiáng)則其有效值增大.隨著金剛石砂輪的磨損,砂輪表面的金剛石磨粒鈍化,切除PSZ材料的能力下降,與之對應(yīng)的AE聲發(fā)射信號也會出現(xiàn)更大的波動,即AE信號方差增大.AE信號能量系數(shù),即每個固有模態(tài)函數(shù)的能量所占原始信號總能量的百分比,由于每個固有模態(tài)函數(shù)所對應(yīng)的AE信號主要頻率不同,能量系數(shù)也在一定程度代表了不同頻率段所占據(jù)的能量比重.砂輪的磨損聲發(fā)射信號頻率有從低頻向高頻移動的趨勢,這一特征可以用能量系數(shù)來表征.砂輪輕度磨損與砂輪嚴(yán)重磨損聲發(fā)射信號特征值如圖6~8所示.

    固有模態(tài)函數(shù)有效值 的計算公式:

    式(1)中: i表示數(shù)據(jù)點(diǎn); N表示信號長度即信號總的數(shù)據(jù)點(diǎn)數(shù);j表示第j個固有模態(tài)函數(shù);xj(i)表示第j個固有模態(tài)函數(shù)的第i個數(shù)據(jù)點(diǎn);IMFrms(j)表示第j個固有模態(tài)函數(shù)的有效值.

    固有模態(tài)函數(shù)方差I(lǐng)MFvar的計算公式:

    式(2)中:i表示數(shù)據(jù)點(diǎn); N表示信號長度即信號總的數(shù)據(jù)點(diǎn)數(shù);j表示第j個固有模態(tài)函數(shù);xj(i)表示第 j個固有模態(tài)函數(shù)的第i個數(shù)據(jù)點(diǎn);μj表示第 j個固有模態(tài)函數(shù)的平均值;IMFvar(j)表示第j個固有模態(tài)函數(shù)的方差.

    固有模態(tài)函數(shù)能量系數(shù) 的計算公式:

    式(3)中: i表示數(shù)據(jù)點(diǎn); N表示信號長度即信號總的數(shù)據(jù)點(diǎn)數(shù);j表示第j個固有模態(tài)函數(shù); m表示經(jīng)驗?zāi)B(tài)分解得到的信號數(shù)目;xj(i)表示第j個固有模態(tài)函數(shù)的第i個數(shù)據(jù)點(diǎn);μj表示第j個固有模態(tài)函數(shù)的平均值;IMFpe(j)表示第j個固有模態(tài)函數(shù)的能量系數(shù).

    從圖6可以明顯看出:砂輪磨損嚴(yán)重時,各固有模態(tài)函數(shù)imf1~imf6的有效值IMFrms要大于輕度砂輪磨損的有效值,所以有效值可以作為砂輪磨損狀態(tài)判別的特征量.砂輪的磨損加劇使磨粒切削刃變鈍,磨?;磷饔迷鰪?qiáng)磨削能量增加,聲發(fā)射信號的幅值增大,導(dǎo)致每個固有模態(tài)函數(shù)的IMFrms增加.但是由于磨削PSZ的AE信號主要頻率在100 kHz以內(nèi),因而每個固有模態(tài)函數(shù)的IMFrms增大程度不同.其中imf3~imf5的有效值的增大最為明顯, imf3~imf5的相關(guān)系數(shù)最大,其主要頻率在100 kHz以下并且占據(jù)了總能量的85%以上,是磨削氧化鋯陶瓷時聲發(fā)射信號的主要頻率范圍.因而在砂輪磨損狀態(tài)變化引起聲發(fā)射源變化時,在imf3~imf5上的反應(yīng)更為明顯,而imf1、imf2和imf6的變化較小.

    從圖7可以明顯看出:砂輪磨損嚴(yán)重時,各固有模態(tài)函數(shù)imf1~imf6的方差I(lǐng)MFvar要大于砂輪輕度磨損方差值,是砂輪磨損狀態(tài)的特征量.金剛石砂輪的磨損使金剛石磨粒切削刃變鈍,PSZ陶瓷材料切除變得困難,磨削變得不穩(wěn)定,聲發(fā)射信號出現(xiàn)更大的波動,因而固有模態(tài)函數(shù)方差I(lǐng)MFvar有變大趨勢.但是因為磨削氧化鋯的聲發(fā)射信號的主要頻率不同,在不同的固有模態(tài)函數(shù)方差I(lǐng)MFvar的增大程度不一樣.由于imf3~imf5的主要頻率在100 kHz以下并且占據(jù)了總能量的85%以上,集中了磨削PSZ陶瓷時聲發(fā)射信號的大部分能量,imf3~imf5的方差值的變化最為明顯,而imf1、imf2和imf6則增大相對較小.

    固有模態(tài)函數(shù)的能量系數(shù)IMFpe表示了每個固有模態(tài)函數(shù)的能量所占AE原始信號總能量的百分比,能量系數(shù)越大,則表明該固有模態(tài)函數(shù)對AE原始信號的貢獻(xiàn)越多.從圖8可以看出:在不同砂輪磨損狀態(tài)下固有模態(tài)函數(shù)imf1~imf6的能量系數(shù)各不相同,但是它們都主要集中在了imf3~imf5上,并且砂輪嚴(yán)重磨損與砂輪輕度磨損相比,imf3和imf4的能量系數(shù)增大,imf5的能量減小.這是因為imf3~imf5的相關(guān)系數(shù)大,其主要頻率在100 kHz以下,是磨削氧化鋯陶瓷時聲發(fā)射信號的主要頻率范圍.當(dāng)砂輪磨損嚴(yán)重時,主要頻率有輕微地向高頻部分移動的趨勢,因而imf3和imf4能量比重增加,imf5(主要頻率分布范圍比imf3和imf4低)的能量比重減少.至于imf1、imf2頻率較高,所占原信號的比重不大;而imf6的頻率低.

    4 ??基于LS-SVM的砂輪磨損狀態(tài)監(jiān)測

    支持向量機(jī)(SVM, Support Vector Machine)克服了過擬合、泛化能力差等問題,在解決小樣本、非線性、局部極小和高維模式識別問題中優(yōu)勢明顯.但是標(biāo)準(zhǔn)支持向量機(jī)SVM訓(xùn)練樣本時因為求解二次規(guī)劃問題而訓(xùn)練速度慢,在磨削實(shí)時監(jiān)測中對大量的聲發(fā)射信號數(shù)據(jù)可能來不及處理,為此應(yīng)用最小二乘支持向量機(jī)(LS-SVM, Least Squares Support Vector Machine)[1]來提高訓(xùn)練速度,為工程陶瓷磨削金剛石砂輪磨損狀態(tài)的智能在線監(jiān)測打下基礎(chǔ).

    實(shí)驗獲得不同砂輪磨損狀態(tài)的磨削聲發(fā)射信號共40組,其中砂輪輕度磨損20組,砂輪嚴(yán)重磨損20組.經(jīng)過EMD分解從固有模態(tài)函數(shù)imf1~imf6中分別提取有效值(IMFrms)、方差(IMFvar)、能量系數(shù)(IMFpe).有效值(IMFrms),方差(IMFvar),能量系數(shù)(IMFpe)各為一組,分別作為LS-SVM的輸入特征向量數(shù)據(jù)集,見表2~4.每組輸入特征向量數(shù)據(jù)集的原始數(shù)據(jù)集被分為兩部分:一部分共30組用于LS-SVM的訓(xùn)練,其中砂輪嚴(yán)重磨損和砂輪輕度磨損狀態(tài)各15組;另一部分10組用于砂輪磨損狀態(tài)的預(yù)測,其中砂輪嚴(yán)重磨損和砂輪輕度磨損狀態(tài)各5組.每個輸入特征向量由6個特征值組成,分別為固有模態(tài)函數(shù)imf1~imf6的有效值(IMFrms)或方差(IMFvar)或能量系數(shù)(IMFpe).

    為了減少AE信號數(shù)據(jù)大小差距過大對LS-SVM識別精度的影響,預(yù)先對磨削AE信號原始數(shù)據(jù)進(jìn)行歸一化處理,采用最小-最大標(biāo)準(zhǔn)化方法,將原始數(shù)據(jù)映射到[0 1]之間,即

    式(4)中: X為原始數(shù)據(jù);Xmax為原始數(shù)據(jù)中的最大值;Xmin為原始數(shù)據(jù)中的最小值; X*為歸一化的數(shù)據(jù).砂輪的磨損狀態(tài)為兩類,用“0”表示砂輪輕度磨損,用“1”表示砂輪嚴(yán)重磨損.由以上得到LS-SVM的樣本集{(xi,yi)},其中xi是由AE信號歸一化后的LS-SVM輸入特征組成的矩陣,代表第i個輸入特征向量;yi∈{0,1},為分類標(biāo)簽代表第i個輸入特征向量對應(yīng)的砂輪磨損狀態(tài).選用徑向基核函數(shù)作為LS-SVM的核函數(shù);正則化參數(shù)c和核函數(shù)參數(shù)σ是需要優(yōu)化的兩個重要參數(shù),采用交叉驗證的方法得到,解決了LS-SVM欠學(xué)習(xí)和過學(xué)習(xí)的問題.

    從金剛石砂輪磨損狀態(tài)預(yù)測結(jié)果(見表5)來看,輸入特征分別為有效值(IMFrms)或方差(IMFvar)或能量系數(shù)(IMFpe)時,磨削金剛石砂輪磨損狀態(tài)的LS-SVM預(yù)測正確率都達(dá)到了100%,很好地實(shí)現(xiàn)了金剛石砂輪磨損狀態(tài)監(jiān)測;同時也說明了所提取磨削聲發(fā)射信號特征對砂輪磨損狀態(tài)監(jiān)測的合理性.

    5 ??結(jié) ??論

    砂輪磨損加劇會導(dǎo)致磨粒摩擦作用增強(qiáng)而切削作用下降,反映在PSZ陶瓷磨削AE信號高頻率分量增強(qiáng)明顯,并且有從低頻向高頻轉(zhuǎn)移的趨勢.結(jié)合工程陶瓷部分穩(wěn)定氧化鋯PSZ磨削過程中AE信號非線性非平穩(wěn)性的特點(diǎn),采用經(jīng)驗?zāi)B(tài)分解EMD方法將其AE信號分解為多個平穩(wěn)的固有模態(tài)函數(shù)之和.隨著金剛石砂輪從輕度磨損轉(zhuǎn)變?yōu)閲?yán)重磨損,即砂輪磨損加劇時,其磨削聲發(fā)射信號各固有模態(tài)函數(shù)的有效值和方差不同程度的增大,能量系數(shù)在imf3~imf5上亦表現(xiàn)出由低頻向高頻移動的特點(diǎn).分別用磨削聲發(fā)射信號的固有模態(tài)函數(shù)imf1~imf6的有效值(IMFrms)、方差(IMFvar)和能量系數(shù)(IMFpe)作為特征向量,用最小二乘支持向量機(jī)智能監(jiān)測PSZ磨削中金剛石砂輪輕度磨損和嚴(yán)重磨損狀態(tài)取得了很好的效果.

    參考文獻(xiàn)

    [1] ???聶鵬,董慧,李正強(qiáng),等. 基于改進(jìn) EMD 和 LS-SVM 的刀具磨損狀態(tài)識別[J]. 北京工業(yè)大學(xué)學(xué)報,2013,39(12):1784-1790.

    NIE P,DONG H,LI Z Q,et al. State recognition of tool wear based on improved empirical mode decomposition and least squares support vector machine[J]. Journal of Beijing University of Technology,2013,39(12):1784-1790.(In Chinese)

    [2] ???NAKAI M E,AGUIAR P R,GUILLARDI H,et al. Evaluation of neural models applied to the estimation of tool wear in the grinding of advanced ceramics[J]. Expert Systems with Applications,2015,42:7026-7035.

    [3] ???聶鵬,徐洪圭,劉新宇,等. EEMD方法在刀具磨損狀態(tài)識別的應(yīng)用[J]. 傳感器與微系統(tǒng),2012,31(5):147-149.

    NIE P,XU H Y,LIU X Y,et al. Application of EEMD method in state recognition of tool wear[J]. Transducer and Microsystem Technologies,2012,31(5):147-149.(In Chinese)

    [4] ???李慶,宋萬清. IMF奇異值熵的刀具磨損狀態(tài)診斷[J]. 制造業(yè)自動化,2013,35(12):52-55.

    LI Q,SONG W Q. Tool wear state diagnosis with IMF singular value entropy[J]. Manufacturing Automation,2013,35(12):52-55.(In Chinese)

    [5] ???楊振生. 面向磨削燒傷問題的間接監(jiān)測技術(shù)研究[D]. 杭州:浙江大學(xué),2013.

    YANG Z S. ?Research on the indirect monitoring technique of grinding burn and its interrelated problems in precision grinding process[D]. Hangzhou: Zhejiang University,2013.(In Chinese)

    [6] ???孫惠斌,牛偉龍,王俊陽. 基于希爾伯特黃變換的刀具磨損特征提取[J]. 振動與沖擊,2015,34(4):158-164.

    SUN H B,NIU W L,WANG J Y. Tool wear feature extraction based on Hilbert-Huang transformation[J]. Journal of Vibration and Shock,2015,34(4):158-164.(In Chinese)

    [7] ???關(guān)山,彭昶. 刀具磨損聲發(fā)射信號的混沌特性分析[J]. 農(nóng)業(yè)工程學(xué)報,2015,31(11):60-65

    GUAN S,PENG C. Chaotic characteristic analysis of tool wear acoustic emission signal[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE),2015,31(11): 60-65. (In Chinese)

    [8] ???CHEN X,?魻P?魻Z T T. Effect of different parameters on grinding efficiency and its monitoring by acoustic emission[J]. Production & Manufacturing Research,2016,4(1):190-208.

    [9] ???彭延峰,劉貞濤,程軍圣,等. 基于初值優(yōu)化的自適應(yīng)最稀疏時頻分析方法[J]. 湖南大學(xué)學(xué)報(自然科學(xué)版),2017, 44(8):50-56.

    PENG Y F,LIU Z T,CHENG J S,et al. Adaptive and sparsest time-frequency analysis method based on initial value optimization [J]. Journal of Hunan University(Natural Sciences),2017,44(8):50-56.(In Chinese)

    [10] ?LIAO T W. Warren Liao. Feature extraction and selection from acoustic emission signals with an application in grinding wheel condition monitoring[J]. Engineering Applications of Artificial Intelligence, 2010,23:74-84.

    [11] ?石建,丁寧. 基于聲發(fā)射技術(shù)的砂輪磨損狀況在線檢測[J]. 長春大學(xué)學(xué)報,2013,23(8):931-936.

    SHI J,DING N. On-line detection of the state of grinding wheel wear based on acoustic emission technique[J]. Journal of Changchun University,2013,23(8):931-936.(In Chinese)

    [12] ?MOHAMMED A,F(xiàn)OLKES J,CHEN X. ?Detection of grinding temperatures using laser irradiation and acoustic emission sensing technique[J]. Materials and Manufacturing Processes,2012,27(4):395-400.

    [13] ?CHEN X,GRIFFIN J,LIU Q. ?Mechanical and thermal behaviors of grinding acoustic emission[J]. International Journal of Manufacturing Technology and Management,2007,12(1/2/3):184-199.

    老司机影院成人| 偷拍熟女少妇极品色| 久久99精品国语久久久| 国内精品美女久久久久久| 日本免费在线观看一区| 日本免费在线观看一区| 日韩av不卡免费在线播放| 黄色一级大片看看| 亚洲av二区三区四区| 中国国产av一级| av黄色大香蕉| 黄片wwwwww| 韩国av在线不卡| 国产成人精品福利久久| 精品人妻偷拍中文字幕| 高清欧美精品videossex| 噜噜噜噜噜久久久久久91| 亚州av有码| 国产欧美亚洲国产| 乱码一卡2卡4卡精品| 97热精品久久久久久| 中文欧美无线码| 欧美老熟妇乱子伦牲交| 久久这里有精品视频免费| 午夜视频国产福利| 久久影院123| 一级毛片aaaaaa免费看小| 丰满人妻一区二区三区视频av| 嘟嘟电影网在线观看| 久久久亚洲精品成人影院| 亚洲va在线va天堂va国产| 五月开心婷婷网| 国产综合精华液| 看黄色毛片网站| 亚洲av福利一区| 又黄又爽又刺激的免费视频.| 国产精品熟女久久久久浪| 亚洲在线观看片| 草草在线视频免费看| 啦啦啦啦在线视频资源| 欧美精品人与动牲交sv欧美| 久久久久久久精品精品| 欧美高清成人免费视频www| 自拍欧美九色日韩亚洲蝌蚪91 | 午夜激情久久久久久久| 欧美高清性xxxxhd video| 亚洲精品,欧美精品| 精品一区二区三区视频在线| 九九在线视频观看精品| 深爱激情五月婷婷| 午夜免费男女啪啪视频观看| 丰满乱子伦码专区| av在线老鸭窝| 少妇人妻久久综合中文| 99热这里只有精品一区| 日本色播在线视频| 成年免费大片在线观看| 国产乱人视频| 好男人在线观看高清免费视频| 三级国产精品欧美在线观看| 久久久久网色| 国产高潮美女av| 欧美xxxx性猛交bbbb| 国产乱来视频区| 久久久久久久国产电影| 亚洲人成网站在线播| 少妇 在线观看| 日日啪夜夜撸| 嘟嘟电影网在线观看| 亚洲美女视频黄频| 国产 一区 欧美 日韩| 18禁裸乳无遮挡动漫免费视频 | 日日啪夜夜撸| 大码成人一级视频| 99久久精品一区二区三区| 国产视频首页在线观看| 国产毛片在线视频| 男男h啪啪无遮挡| 夜夜看夜夜爽夜夜摸| 水蜜桃什么品种好| 国产精品福利在线免费观看| 亚洲不卡免费看| 自拍欧美九色日韩亚洲蝌蚪91 | 国产精品一区二区性色av| 亚洲欧美日韩另类电影网站 | 女人被狂操c到高潮| 九九久久精品国产亚洲av麻豆| 观看免费一级毛片| 免费av不卡在线播放| 久久国内精品自在自线图片| 三级国产精品欧美在线观看| 午夜激情福利司机影院| 日本色播在线视频| tube8黄色片| 久久ye,这里只有精品| 成年版毛片免费区| 国产精品嫩草影院av在线观看| 欧美潮喷喷水| 国产v大片淫在线免费观看| 赤兔流量卡办理| 欧美极品一区二区三区四区| 乱系列少妇在线播放| 午夜福利视频精品| 亚洲欧洲日产国产| 在线观看一区二区三区| 另类亚洲欧美激情| 美女内射精品一级片tv| 亚洲一区二区三区欧美精品 | 寂寞人妻少妇视频99o| 干丝袜人妻中文字幕| 丰满人妻一区二区三区视频av| 亚洲精品一二三| 美女内射精品一级片tv| 神马国产精品三级电影在线观看| 国产高潮美女av| 国产伦在线观看视频一区| 高清在线视频一区二区三区| 国产爽快片一区二区三区| 中文字幕久久专区| 一个人观看的视频www高清免费观看| 国内精品美女久久久久久| 成人亚洲精品av一区二区| 又粗又硬又长又爽又黄的视频| 搞女人的毛片| 在线 av 中文字幕| 少妇人妻 视频| kizo精华| 高清欧美精品videossex| 日本午夜av视频| 成人欧美大片| 天天躁日日操中文字幕| 日日啪夜夜撸| 在线观看一区二区三区激情| 一区二区三区四区激情视频| 交换朋友夫妻互换小说| 免费黄网站久久成人精品| 国产男女内射视频| 亚洲成人av在线免费| 天堂俺去俺来也www色官网| 我的女老师完整版在线观看| 伦精品一区二区三区| 久久久久网色| www.色视频.com| 国产成人freesex在线| 少妇 在线观看| 中文精品一卡2卡3卡4更新| 国产毛片a区久久久久| 日韩 亚洲 欧美在线| 丰满人妻一区二区三区视频av| 99re6热这里在线精品视频| 下体分泌物呈黄色| 精品一区二区免费观看| 国产精品人妻久久久久久| 色婷婷久久久亚洲欧美| 一区二区av电影网| 干丝袜人妻中文字幕| 一本色道久久久久久精品综合| 国产一区有黄有色的免费视频| 日日撸夜夜添| 精品亚洲乱码少妇综合久久| 欧美日韩视频高清一区二区三区二| 国产午夜福利久久久久久| 久久久久久久国产电影| 国产精品女同一区二区软件| 国产免费一级a男人的天堂| 亚洲美女搞黄在线观看| a级毛片免费高清观看在线播放| 国产av码专区亚洲av| 国产精品蜜桃在线观看| 国产高清国产精品国产三级 | 永久免费av网站大全| 精品人妻一区二区三区麻豆| 国精品久久久久久国模美| 婷婷色麻豆天堂久久| 精品视频人人做人人爽| 国产亚洲av嫩草精品影院| 视频区图区小说| 免费大片黄手机在线观看| 又爽又黄无遮挡网站| 91久久精品国产一区二区三区| 国产女主播在线喷水免费视频网站| 97在线人人人人妻| 嫩草影院入口| 少妇人妻 视频| 岛国毛片在线播放| 亚洲精品中文字幕在线视频 | xxx大片免费视频| 欧美日韩一区二区视频在线观看视频在线 | 最近最新中文字幕免费大全7| 少妇高潮的动态图| 在线亚洲精品国产二区图片欧美 | 最近中文字幕2019免费版| 中文字幕免费在线视频6| 国精品久久久久久国模美| 美女主播在线视频| 日本熟妇午夜| 亚洲国产成人一精品久久久| 最近中文字幕高清免费大全6| 国产精品秋霞免费鲁丝片| 天美传媒精品一区二区| 免费观看性生交大片5| 欧美日韩视频高清一区二区三区二| 欧美性感艳星| 黄色怎么调成土黄色| 嫩草影院新地址| 日韩不卡一区二区三区视频在线| 人人妻人人澡人人爽人人夜夜| 尤物成人国产欧美一区二区三区| 最近手机中文字幕大全| 777米奇影视久久| 国产乱人视频| 麻豆成人av视频| 欧美日韩亚洲高清精品| av卡一久久| www.色视频.com| 午夜福利在线观看免费完整高清在| 免费少妇av软件| 日韩,欧美,国产一区二区三区| 免费观看在线日韩| 国产精品一区二区性色av| av.在线天堂| 纵有疾风起免费观看全集完整版| 欧美精品一区二区大全| 午夜福利高清视频| 99久久精品热视频| 精品国产露脸久久av麻豆| 建设人人有责人人尽责人人享有的 | 好男人在线观看高清免费视频| 国产精品爽爽va在线观看网站| 日本午夜av视频| 亚洲电影在线观看av| 精品久久久久久久人妻蜜臀av| 国产精品国产三级国产专区5o| 久久97久久精品| 久久久久久久久大av| 纵有疾风起免费观看全集完整版| 日韩在线高清观看一区二区三区| 中文字幕制服av| 日韩欧美精品免费久久| 一级a做视频免费观看| 国产精品成人在线| 国产高清不卡午夜福利| 国产淫语在线视频| 久久精品久久久久久久性| 亚洲高清免费不卡视频| 中文精品一卡2卡3卡4更新| 日日啪夜夜撸| 一区二区三区免费毛片| 国产一区二区在线观看日韩| 免费观看在线日韩| 国产伦理片在线播放av一区| 别揉我奶头 嗯啊视频| 成年人午夜在线观看视频| 色视频在线一区二区三区| 久久久久久伊人网av| 亚洲激情五月婷婷啪啪| 大香蕉久久网| av专区在线播放| 中文字幕免费在线视频6| 国产老妇伦熟女老妇高清| 久久久久久久精品精品| 免费电影在线观看免费观看| 男女边摸边吃奶| 黄色日韩在线| 18+在线观看网站| 尤物成人国产欧美一区二区三区| 久久精品夜色国产| 王馨瑶露胸无遮挡在线观看| av一本久久久久| 少妇人妻久久综合中文| 久久久午夜欧美精品| 色网站视频免费| 国产毛片在线视频| 免费少妇av软件| 中文字幕亚洲精品专区| 久久国内精品自在自线图片| 亚洲av一区综合| 国产av不卡久久| 亚洲欧美成人综合另类久久久| 精品人妻偷拍中文字幕| 三级男女做爰猛烈吃奶摸视频| 18+在线观看网站| 男女那种视频在线观看| 欧美 日韩 精品 国产| 丝袜脚勾引网站| 国产免费一区二区三区四区乱码| 激情 狠狠 欧美| 爱豆传媒免费全集在线观看| 亚洲精品成人久久久久久| 交换朋友夫妻互换小说| av天堂中文字幕网| 国产精品99久久久久久久久| 亚洲在线观看片| 国产久久久一区二区三区| 国产乱来视频区| 中文乱码字字幕精品一区二区三区| 国产免费视频播放在线视频| 久久影院123| 国产男女内射视频| 我的女老师完整版在线观看| 久久国产乱子免费精品| 天天躁夜夜躁狠狠久久av| 国产精品久久久久久久电影| 黄色一级大片看看| 亚洲无线观看免费| a级一级毛片免费在线观看| 一级爰片在线观看| 亚洲伊人色综图| 飞空精品影院首页| 国产亚洲一区二区精品| 黄色视频在线播放观看不卡| 亚洲精品国产av成人精品| 国产1区2区3区精品| 国产成人av激情在线播放| www.熟女人妻精品国产| 狂野欧美激情性xxxx| 菩萨蛮人人尽说江南好唐韦庄| 国产 一区精品| 免费高清在线观看视频在线观看| 国产成人精品无人区| 国产男女超爽视频在线观看| 啦啦啦 在线观看视频| 亚洲在久久综合| 国产成人一区二区在线| 少妇被粗大的猛进出69影院| 国产淫语在线视频| a级毛片在线看网站| 亚洲av电影在线观看一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 欧美亚洲日本最大视频资源| 欧美黑人精品巨大| 国产av精品麻豆| 日本午夜av视频| 天堂8中文在线网| xxxhd国产人妻xxx| 国产免费视频播放在线视频| 久久韩国三级中文字幕| 热re99久久国产66热| 午夜日韩欧美国产| 亚洲婷婷狠狠爱综合网| 亚洲婷婷狠狠爱综合网| 久久精品国产综合久久久| 日韩欧美精品免费久久| 国产成人精品福利久久| 少妇被粗大的猛进出69影院| 欧美变态另类bdsm刘玥| 午夜福利乱码中文字幕| 考比视频在线观看| 这个男人来自地球电影免费观看 | 女的被弄到高潮叫床怎么办| 国产成人午夜福利电影在线观看| 免费观看人在逋| 欧美日韩福利视频一区二区| 搡老乐熟女国产| 最近2019中文字幕mv第一页| 国产精品国产av在线观看| 亚洲,欧美,日韩| 19禁男女啪啪无遮挡网站| 亚洲成色77777| 狠狠婷婷综合久久久久久88av| 国产福利在线免费观看视频| 久久天躁狠狠躁夜夜2o2o | 亚洲精品一二三| 捣出白浆h1v1| 国产高清不卡午夜福利| 精品免费久久久久久久清纯 | 亚洲一卡2卡3卡4卡5卡精品中文| 操出白浆在线播放| 精品一区二区免费观看| 国产一区亚洲一区在线观看| www.精华液| 青青草视频在线视频观看| 国产成人免费无遮挡视频| 午夜日韩欧美国产| 免费在线观看视频国产中文字幕亚洲 | 美女主播在线视频| 国产成人免费观看mmmm| 又粗又硬又长又爽又黄的视频| 精品一区在线观看国产| 看免费av毛片| 婷婷色综合www| 亚洲欧美日韩另类电影网站| 如何舔出高潮| 亚洲av日韩在线播放| 国产精品一区二区在线不卡| 亚洲欧美色中文字幕在线| 一区二区三区精品91| 啦啦啦在线免费观看视频4| 最近手机中文字幕大全| 国产精品免费视频内射| 丝袜在线中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 国产又色又爽无遮挡免| 女人久久www免费人成看片| 中文字幕另类日韩欧美亚洲嫩草| 麻豆av在线久日| 人人妻人人澡人人看| 日韩中文字幕视频在线看片| 中国国产av一级| 王馨瑶露胸无遮挡在线观看| 国产精品久久久久久精品电影小说| 啦啦啦啦在线视频资源| 乱人伦中国视频| 午夜免费观看性视频| 亚洲精品美女久久av网站| 国产有黄有色有爽视频| 这个男人来自地球电影免费观看 | 亚洲精品日本国产第一区| 国产精品久久久久久精品古装| 国产亚洲精品第一综合不卡| 国产又色又爽无遮挡免| 免费高清在线观看日韩| 亚洲av国产av综合av卡| 一区二区av电影网| 亚洲精品国产色婷婷电影| 欧美日韩国产mv在线观看视频| 777久久人妻少妇嫩草av网站| 女的被弄到高潮叫床怎么办| 精品一区二区三卡| 岛国毛片在线播放| 国产精品 国内视频| 国产一区二区在线观看av| 国产精品 欧美亚洲| 中文字幕高清在线视频| 亚洲欧美中文字幕日韩二区| 女的被弄到高潮叫床怎么办| 日韩电影二区| 1024香蕉在线观看| 国产精品香港三级国产av潘金莲 | 别揉我奶头~嗯~啊~动态视频 | kizo精华| 精品国产一区二区三区四区第35| 成年av动漫网址| 超碰97精品在线观看| 欧美精品亚洲一区二区| 亚洲自偷自拍图片 自拍| 黄色毛片三级朝国网站| 久久热在线av| 黑人猛操日本美女一级片| 精品福利永久在线观看| 久久女婷五月综合色啪小说| 国产探花极品一区二区| 99国产精品免费福利视频| 纵有疾风起免费观看全集完整版| 欧美日韩福利视频一区二区| 激情五月婷婷亚洲| 夜夜骑夜夜射夜夜干| 悠悠久久av| 国产不卡av网站在线观看| av.在线天堂| 日本wwww免费看| 国产成人av激情在线播放| a 毛片基地| 精品人妻一区二区三区麻豆| 精品人妻熟女毛片av久久网站| 在线看a的网站| 国产成人91sexporn| 精品少妇久久久久久888优播| 97人妻天天添夜夜摸| 无遮挡黄片免费观看| 黄色视频不卡| 久久性视频一级片| 叶爱在线成人免费视频播放| 亚洲精品乱久久久久久| 下体分泌物呈黄色| 免费在线观看完整版高清| 人人妻人人澡人人看| 日本欧美国产在线视频| 免费看不卡的av| 久久精品久久久久久噜噜老黄| 午夜福利影视在线免费观看| 亚洲国产精品国产精品| 90打野战视频偷拍视频| 日韩不卡一区二区三区视频在线| 综合色丁香网| 婷婷色av中文字幕| 亚洲欧美成人精品一区二区| 天天操日日干夜夜撸| 成人三级做爰电影| 亚洲四区av| 欧美精品高潮呻吟av久久| 天天添夜夜摸| 国产成人91sexporn| 丝袜美足系列| 99久久人妻综合| 成人亚洲精品一区在线观看| 日本一区二区免费在线视频| av在线播放精品| 欧美 日韩 精品 国产| 国产一区二区三区综合在线观看| 欧美精品av麻豆av| 国产淫语在线视频| 亚洲国产av新网站| 久久久久精品性色| 看免费成人av毛片| 2018国产大陆天天弄谢| 欧美 日韩 精品 国产| 国产av码专区亚洲av| 制服人妻中文乱码| 极品少妇高潮喷水抽搐| 只有这里有精品99| 午夜91福利影院| 天天躁狠狠躁夜夜躁狠狠躁| 欧美xxⅹ黑人| 午夜久久久在线观看| 侵犯人妻中文字幕一二三四区| 国产色婷婷99| 天天躁夜夜躁狠狠久久av| 久久久精品94久久精品| 少妇人妻久久综合中文| 男女边吃奶边做爰视频| 久久青草综合色| 2021少妇久久久久久久久久久| 国产精品一区二区精品视频观看| 国产成人欧美| 亚洲国产看品久久| 久久鲁丝午夜福利片| 久久99一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲一级一片aⅴ在线观看| 精品一区在线观看国产| 嫩草影院入口| 韩国av在线不卡| 波多野结衣av一区二区av| 99re6热这里在线精品视频| 国产国语露脸激情在线看| 欧美乱码精品一区二区三区| 国产免费又黄又爽又色| 1024视频免费在线观看| 久久青草综合色| 国产野战对白在线观看| 老鸭窝网址在线观看| 精品一区二区三区av网在线观看 | 天天躁夜夜躁狠狠久久av| 青青草视频在线视频观看| 久久人人97超碰香蕉20202| 黄频高清免费视频| av天堂久久9| 黑人巨大精品欧美一区二区蜜桃| 成年美女黄网站色视频大全免费| 国产片内射在线| 亚洲专区中文字幕在线 | 国产日韩一区二区三区精品不卡| 999久久久国产精品视频| 天堂中文最新版在线下载| 久久国产亚洲av麻豆专区| 日日啪夜夜爽| 亚洲激情五月婷婷啪啪| 丁香六月欧美| 欧美97在线视频| 香蕉国产在线看| 最近手机中文字幕大全| 国产精品国产av在线观看| 国产欧美日韩综合在线一区二区| 99久国产av精品国产电影| 国产不卡av网站在线观看| 国产精品秋霞免费鲁丝片| 国产 精品1| 女人久久www免费人成看片| 各种免费的搞黄视频| 丰满饥渴人妻一区二区三| 欧美日韩综合久久久久久| 国产精品三级大全| 伊人久久大香线蕉亚洲五| 久久久久久久久久久久大奶| 日日爽夜夜爽网站| 久久婷婷青草| 亚洲欧洲国产日韩| av女优亚洲男人天堂| 黄色一级大片看看| 美女扒开内裤让男人捅视频| 久久久久久人人人人人| 欧美精品人与动牲交sv欧美| 啦啦啦视频在线资源免费观看| 久久精品熟女亚洲av麻豆精品| 99香蕉大伊视频| 在现免费观看毛片| 国产在线免费精品| 男女高潮啪啪啪动态图| 极品少妇高潮喷水抽搐| 亚洲美女黄色视频免费看| 老司机亚洲免费影院| 大香蕉久久网| 国产视频首页在线观看| 丁香六月天网| 观看av在线不卡| 少妇 在线观看| 中文字幕另类日韩欧美亚洲嫩草| 黄片播放在线免费| 日韩精品免费视频一区二区三区| 又大又黄又爽视频免费| 欧美成人午夜精品| 欧美xxⅹ黑人| 精品少妇内射三级| 高清视频免费观看一区二区| 午夜激情久久久久久久| 国产深夜福利视频在线观看| 一区二区三区精品91| 男人爽女人下面视频在线观看| 日韩熟女老妇一区二区性免费视频| 搡老乐熟女国产| 美女午夜性视频免费| 视频区图区小说| 日韩一区二区视频免费看| 欧美日韩一级在线毛片| 韩国av在线不卡| 人成视频在线观看免费观看| 午夜91福利影院| 黄色一级大片看看| 欧美久久黑人一区二区| 亚洲精品久久成人aⅴ小说| 十分钟在线观看高清视频www| 在线亚洲精品国产二区图片欧美| 大香蕉久久成人网| 国产精品一区二区精品视频观看|