• <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.

    freevideosex欧美| 18+在线观看网站| 欧美成人a在线观看| 爱豆传媒免费全集在线观看| 国产黄片美女视频| 女的被弄到高潮叫床怎么办| 三级毛片av免费| 午夜福利视频精品| 自拍偷自拍亚洲精品老妇| 日本黄色片子视频| 亚洲怡红院男人天堂| av.在线天堂| 免费黄色在线免费观看| 亚洲内射少妇av| 国产免费福利视频在线观看| 午夜久久久久精精品| 在线观看免费高清a一片| 三级经典国产精品| 国产在视频线在精品| 一级毛片久久久久久久久女| 两个人的视频大全免费| 中文在线观看免费www的网站| 欧美日韩一区二区视频在线观看视频在线 | 亚洲精品自拍成人| 男女那种视频在线观看| 最后的刺客免费高清国语| 美女黄网站色视频| 亚洲第一区二区三区不卡| 国产精品99久久久久久久久| 九九爱精品视频在线观看| 一级片'在线观看视频| 国产有黄有色有爽视频| 亚洲精品456在线播放app| 久久精品综合一区二区三区| 国产爱豆传媒在线观看| 日韩成人伦理影院| 99久国产av精品| 国内精品宾馆在线| 国产成年人精品一区二区| 美女脱内裤让男人舔精品视频| 国产91av在线免费观看| 欧美成人一区二区免费高清观看| 国产精品久久久久久精品电影| av.在线天堂| 欧美变态另类bdsm刘玥| 精品亚洲乱码少妇综合久久| 精品欧美国产一区二区三| 一级片'在线观看视频| 亚洲熟女精品中文字幕| 亚洲欧美日韩东京热| 如何舔出高潮| 一区二区三区乱码不卡18| 色哟哟·www| 日本wwww免费看| 高清视频免费观看一区二区 | 久久久精品免费免费高清| 婷婷色av中文字幕| 激情 狠狠 欧美| 草草在线视频免费看| 永久免费av网站大全| 欧美激情久久久久久爽电影| 精品人妻一区二区三区麻豆| 乱系列少妇在线播放| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产av新网站| 亚洲精品自拍成人| 日日摸夜夜添夜夜添av毛片| 国产精品麻豆人妻色哟哟久久 | 男女那种视频在线观看| 亚洲人与动物交配视频| 久久韩国三级中文字幕| 精品国产三级普通话版| 黄色日韩在线| 国产成人精品一,二区| 精品欧美国产一区二区三| 国产高清国产精品国产三级 | 国产精品三级大全| 久久久亚洲精品成人影院| 黄色一级大片看看| 免费看日本二区| 国产黄色视频一区二区在线观看| 日本免费a在线| 色吧在线观看| 肉色欧美久久久久久久蜜桃 | 午夜激情欧美在线| 免费观看性生交大片5| 91午夜精品亚洲一区二区三区| 国内精品一区二区在线观看| 午夜福利成人在线免费观看| 丰满少妇做爰视频| 一级a做视频免费观看| av免费观看日本| 国产爱豆传媒在线观看| 一个人看的www免费观看视频| 人人妻人人看人人澡| 久久久色成人| 啦啦啦中文免费视频观看日本| 亚洲自拍偷在线| 欧美性感艳星| 午夜福利高清视频| 精品国内亚洲2022精品成人| 精品一区二区三区视频在线| 久久6这里有精品| 欧美精品国产亚洲| 天天躁夜夜躁狠狠久久av| 国产毛片a区久久久久| 观看免费一级毛片| 久久久久久久大尺度免费视频| 噜噜噜噜噜久久久久久91| 久久国产乱子免费精品| 夜夜看夜夜爽夜夜摸| 精品少妇黑人巨大在线播放| 久99久视频精品免费| 天美传媒精品一区二区| 亚洲成人精品中文字幕电影| 国产一级毛片七仙女欲春2| 亚洲性久久影院| 国产大屁股一区二区在线视频| 高清欧美精品videossex| 亚洲av电影在线观看一区二区三区 | 啦啦啦中文免费视频观看日本| 婷婷六月久久综合丁香| 免费观看a级毛片全部| 亚洲精品亚洲一区二区| 男人舔奶头视频| 男女啪啪激烈高潮av片| 亚洲综合色惰| av卡一久久| 亚洲国产最新在线播放| 男女视频在线观看网站免费| 精品久久久久久久末码| 美女cb高潮喷水在线观看| 午夜老司机福利剧场| 亚洲av免费高清在线观看| 精品人妻熟女av久视频| 韩国av在线不卡| 亚洲美女搞黄在线观看| 超碰97精品在线观看| 女的被弄到高潮叫床怎么办| 国产久久久一区二区三区| 日本一本二区三区精品| 欧美一级a爱片免费观看看| 亚洲在线观看片| 国产在视频线在精品| 成人综合一区亚洲| 欧美变态另类bdsm刘玥| 久久久久久久久中文| 不卡视频在线观看欧美| 成人午夜精彩视频在线观看| 国产探花在线观看一区二区| 99热这里只有是精品50| 午夜视频国产福利| 91狼人影院| 亚洲经典国产精华液单| 国产老妇女一区| 女人久久www免费人成看片| 久久鲁丝午夜福利片| 亚洲色图av天堂| 肉色欧美久久久久久久蜜桃 | 99久久人妻综合| 狂野欧美白嫩少妇大欣赏| 亚洲精品中文字幕在线视频 | 国产成人免费观看mmmm| 亚洲熟妇中文字幕五十中出| xxx大片免费视频| 综合色av麻豆| 亚洲内射少妇av| 亚洲av电影在线观看一区二区三区 | 亚洲欧美精品自产自拍| 亚洲乱码一区二区免费版| 精品久久国产蜜桃| 18禁在线无遮挡免费观看视频| 有码 亚洲区| 午夜亚洲福利在线播放| 国产成年人精品一区二区| 免费不卡的大黄色大毛片视频在线观看 | 日韩成人av中文字幕在线观看| 亚洲精品国产成人久久av| 一个人免费在线观看电影| 女人被狂操c到高潮| 久久久欧美国产精品| 国产爱豆传媒在线观看| 亚洲国产av新网站| 激情 狠狠 欧美| 国产成人a∨麻豆精品| 麻豆成人午夜福利视频| 国产麻豆成人av免费视频| 亚洲欧美一区二区三区国产| 九九爱精品视频在线观看| 亚洲性久久影院| 国产精品99久久久久久久久| 欧美xxⅹ黑人| 国产午夜精品久久久久久一区二区三区| 成人国产麻豆网| 亚洲欧洲国产日韩| 有码 亚洲区| 少妇裸体淫交视频免费看高清| 久久99热这里只频精品6学生| 午夜福利视频精品| 日日摸夜夜添夜夜爱| 亚洲高清免费不卡视频| 久久久久精品性色| 国产成人a∨麻豆精品| 在现免费观看毛片| 亚洲精品第二区| 欧美日韩亚洲高清精品| 日韩欧美国产在线观看| 日韩欧美国产在线观看| 看黄色毛片网站| 精品酒店卫生间| 欧美区成人在线视频| 中文字幕av成人在线电影| 自拍偷自拍亚洲精品老妇| 亚洲国产精品sss在线观看| www.色视频.com| 欧美另类一区| 91久久精品国产一区二区成人| videossex国产| 亚洲精品日本国产第一区| 老司机影院成人| a级毛色黄片| 丝瓜视频免费看黄片| 国产在视频线在精品| 十八禁国产超污无遮挡网站| 日日摸夜夜添夜夜添av毛片| 欧美激情国产日韩精品一区| 大香蕉久久网| 一区二区三区高清视频在线| 久久99精品国语久久久| 男人舔女人下体高潮全视频| 高清av免费在线| 观看美女的网站| 好男人在线观看高清免费视频| 人妻少妇偷人精品九色| 亚洲一区高清亚洲精品| 久久精品国产亚洲av天美| 夜夜爽夜夜爽视频| 欧美+日韩+精品| 三级国产精品欧美在线观看| 亚洲国产精品sss在线观看| 亚洲av福利一区| 观看美女的网站| 美女cb高潮喷水在线观看| 韩国av在线不卡| 观看免费一级毛片| 啦啦啦韩国在线观看视频| 欧美+日韩+精品| 亚洲欧美一区二区三区黑人 | 久久99热这里只频精品6学生| av在线播放精品| 国产成人a∨麻豆精品| 美女cb高潮喷水在线观看| 亚洲av不卡在线观看| 伊人久久国产一区二区| 可以在线观看毛片的网站| 日韩国内少妇激情av| 熟女电影av网| 日韩欧美精品v在线| 亚洲av免费在线观看| 超碰av人人做人人爽久久| 一级毛片我不卡| 亚洲人成网站在线播| 亚洲国产日韩欧美精品在线观看| 亚洲av日韩在线播放| av免费观看日本| 我的女老师完整版在线观看| 少妇的逼好多水| 黑人高潮一二区| 狂野欧美激情性xxxx在线观看| 七月丁香在线播放| 搡老妇女老女人老熟妇| 国产淫片久久久久久久久| 久久久久久久大尺度免费视频| 国产成人a区在线观看| 亚洲精品一二三| 你懂的网址亚洲精品在线观看| 久久国产乱子免费精品| 97超视频在线观看视频| 成人特级av手机在线观看| 91精品伊人久久大香线蕉| 精品欧美国产一区二区三| 成人二区视频| 精品一区二区三卡| 国产高清不卡午夜福利| 成人毛片a级毛片在线播放| 精品一区二区免费观看| 亚洲熟妇中文字幕五十中出| 两个人视频免费观看高清| 熟女人妻精品中文字幕| 亚洲美女搞黄在线观看| 波多野结衣巨乳人妻| 美女脱内裤让男人舔精品视频| 亚洲精品自拍成人| 国产精品一及| 一级毛片久久久久久久久女| 一级黄片播放器| 久久久久久久亚洲中文字幕| 高清视频免费观看一区二区 | 久久久久久久久久黄片| 日韩av在线大香蕉| 水蜜桃什么品种好| 最近最新中文字幕大全电影3| 少妇高潮的动态图| 国产精品国产三级专区第一集| 日韩av免费高清视频| av卡一久久| 18禁动态无遮挡网站| 久久久久久久久久成人| 人体艺术视频欧美日本| www.色视频.com| 嫩草影院新地址| 午夜福利视频1000在线观看| 午夜精品国产一区二区电影 | 最近视频中文字幕2019在线8| 久久午夜福利片| 特大巨黑吊av在线直播| 午夜视频国产福利| 国产免费一级a男人的天堂| 精华霜和精华液先用哪个| 亚洲精品日韩av片在线观看| 国产欧美另类精品又又久久亚洲欧美| 精品久久久久久久久av| 男插女下体视频免费在线播放| 成年av动漫网址| 18+在线观看网站| av网站免费在线观看视频 | 亚洲内射少妇av| 搡老妇女老女人老熟妇| 超碰97精品在线观看| 中国美白少妇内射xxxbb| 午夜日本视频在线| 国产一区亚洲一区在线观看| 国产精品久久视频播放| 三级国产精品片| 国产伦精品一区二区三区四那| 一级毛片黄色毛片免费观看视频| 免费观看a级毛片全部| 青春草亚洲视频在线观看| 春色校园在线视频观看| 午夜精品国产一区二区电影 | 亚洲国产精品成人综合色| 国产亚洲一区二区精品| 亚洲国产日韩欧美精品在线观看| 日本-黄色视频高清免费观看| 九九久久精品国产亚洲av麻豆| 最近中文字幕高清免费大全6| 三级国产精品片| 91精品一卡2卡3卡4卡| 午夜福利在线在线| 99九九线精品视频在线观看视频| 国产中年淑女户外野战色| 国产亚洲精品av在线| 国内揄拍国产精品人妻在线| 亚洲国产精品成人综合色| 精品人妻一区二区三区麻豆| 大陆偷拍与自拍| 亚洲av男天堂| 国产老妇伦熟女老妇高清| 欧美变态另类bdsm刘玥| 日韩精品青青久久久久久| 91久久精品国产一区二区三区| 亚洲人与动物交配视频| 国产精品人妻久久久久久| 免费观看a级毛片全部| 国产高清有码在线观看视频| 国产91av在线免费观看| 国产熟女欧美一区二区| 国产精品一区二区三区四区久久| 天堂俺去俺来也www色官网 | 免费av观看视频| 国产精品无大码| 成人漫画全彩无遮挡| 国产视频内射| 亚洲,欧美,日韩| 人人妻人人澡欧美一区二区| 美女xxoo啪啪120秒动态图| 秋霞在线观看毛片| 国产探花极品一区二区| 美女大奶头视频| 毛片一级片免费看久久久久| 成人美女网站在线观看视频| ponron亚洲| 中国国产av一级| 欧美 日韩 精品 国产| 国产午夜福利久久久久久| 亚洲自偷自拍三级| 午夜激情欧美在线| 国产黄片美女视频| 亚洲一级一片aⅴ在线观看| 能在线免费观看的黄片| 可以在线观看毛片的网站| 日韩欧美 国产精品| 亚洲精品国产成人久久av| 国产v大片淫在线免费观看| 精品国产一区二区三区久久久樱花 | 国产黄频视频在线观看| 亚洲av一区综合| 国产精品麻豆人妻色哟哟久久 | 午夜福利在线观看吧| 国产精品麻豆人妻色哟哟久久 | 久久午夜福利片| 午夜精品一区二区三区免费看| 精品人妻偷拍中文字幕| 欧美zozozo另类| 精品一区二区三卡| 精品人妻视频免费看| 亚洲婷婷狠狠爱综合网| 又黄又爽又刺激的免费视频.| 国产白丝娇喘喷水9色精品| 高清日韩中文字幕在线| 成人av在线播放网站| 午夜免费男女啪啪视频观看| 国产精品无大码| 性色avwww在线观看| 天堂av国产一区二区熟女人妻| 国产女主播在线喷水免费视频网站 | 在现免费观看毛片| 亚洲图色成人| 国产精品综合久久久久久久免费| 国产高清有码在线观看视频| 亚洲在线自拍视频| 亚洲欧美日韩无卡精品| 99久久精品国产国产毛片| 精品一区二区三卡| 久久久亚洲精品成人影院| 久久97久久精品| 久久久久久久久久久丰满| 国产高潮美女av| 99九九线精品视频在线观看视频| 老师上课跳d突然被开到最大视频| 欧美不卡视频在线免费观看| 在线免费观看的www视频| 国产精品国产三级国产专区5o| 国产白丝娇喘喷水9色精品| 国产亚洲91精品色在线| 亚洲成人一二三区av| 夜夜看夜夜爽夜夜摸| 舔av片在线| 国产亚洲av片在线观看秒播厂 | 欧美日韩精品成人综合77777| 日韩欧美精品免费久久| 秋霞伦理黄片| 国产极品天堂在线| 色播亚洲综合网| 亚洲天堂国产精品一区在线| 久久久久久九九精品二区国产| 两个人视频免费观看高清| 免费黄频网站在线观看国产| 日韩av在线免费看完整版不卡| 色播亚洲综合网| 中文字幕久久专区| 久久精品夜夜夜夜夜久久蜜豆| 午夜激情欧美在线| 看十八女毛片水多多多| 亚洲成人中文字幕在线播放| 久久亚洲国产成人精品v| 久久久午夜欧美精品| 69人妻影院| 成人午夜精彩视频在线观看| 国内揄拍国产精品人妻在线| 国产黄色视频一区二区在线观看| 老司机影院毛片| 校园人妻丝袜中文字幕| 亚洲av成人精品一区久久| 亚洲欧美精品专区久久| 国产片特级美女逼逼视频| 亚洲成色77777| 亚洲最大成人av| 中文在线观看免费www的网站| 深夜a级毛片| 在现免费观看毛片| 成人高潮视频无遮挡免费网站| 久久久国产一区二区| 国产在视频线在精品| 亚洲真实伦在线观看| 日韩欧美精品免费久久| 久久久久国产网址| 国产精品一区二区三区四区久久| 欧美极品一区二区三区四区| freevideosex欧美| 日韩亚洲欧美综合| 日本熟妇午夜| 国产亚洲最大av| 国产日韩欧美在线精品| 哪个播放器可以免费观看大片| 午夜爱爱视频在线播放| 色网站视频免费| 久久精品久久精品一区二区三区| 国产伦精品一区二区三区四那| 99久久中文字幕三级久久日本| 婷婷色av中文字幕| 国产精品综合久久久久久久免费| av卡一久久| 高清视频免费观看一区二区 | 亚洲国产精品成人久久小说| 九草在线视频观看| 舔av片在线| 日韩电影二区| 日韩制服骚丝袜av| 久热久热在线精品观看| 亚洲av国产av综合av卡| 国语对白做爰xxxⅹ性视频网站| 伊人久久精品亚洲午夜| 床上黄色一级片| 内地一区二区视频在线| 亚洲乱码一区二区免费版| 亚洲国产色片| 大话2 男鬼变身卡| videossex国产| 国产毛片a区久久久久| 91久久精品国产一区二区三区| 最近最新中文字幕大全电影3| 亚洲精品成人久久久久久| 少妇熟女欧美另类| 纵有疾风起免费观看全集完整版 | 非洲黑人性xxxx精品又粗又长| 成人毛片60女人毛片免费| 成人无遮挡网站| 国产精品.久久久| 超碰av人人做人人爽久久| 国产探花在线观看一区二区| 亚洲怡红院男人天堂| 99久国产av精品| 国产真实伦视频高清在线观看| 免费黄色在线免费观看| 久久精品人妻少妇| 中文乱码字字幕精品一区二区三区 | 亚洲第一区二区三区不卡| 日日摸夜夜添夜夜爱| 97超碰精品成人国产| 色哟哟·www| 精品国内亚洲2022精品成人| 色吧在线观看| 久久久久久久久久久丰满| 亚洲欧美成人综合另类久久久| 一级av片app| 久久精品人妻少妇| 国产麻豆成人av免费视频| 一区二区三区高清视频在线| 亚洲va在线va天堂va国产| av福利片在线观看| 亚洲欧美日韩卡通动漫| 成人高潮视频无遮挡免费网站| 偷拍熟女少妇极品色| 国产成人一区二区在线| 舔av片在线| 在线播放无遮挡| 久久精品综合一区二区三区| 成人综合一区亚洲| 爱豆传媒免费全集在线观看| av黄色大香蕉| 91午夜精品亚洲一区二区三区| 一级二级三级毛片免费看| 赤兔流量卡办理| 国产成人福利小说| 亚洲丝袜综合中文字幕| 床上黄色一级片| 天堂网av新在线| 久久人人爽人人爽人人片va| 十八禁国产超污无遮挡网站| 国产男人的电影天堂91| 青春草视频在线免费观看| 精品久久久久久久久av| 国产色婷婷99| 午夜激情久久久久久久| 99久久精品一区二区三区| 综合色av麻豆| av在线亚洲专区| 亚洲av一区综合| 亚洲,欧美,日韩| 搡老妇女老女人老熟妇| www.av在线官网国产| 色播亚洲综合网| 高清午夜精品一区二区三区| 日本欧美国产在线视频| 久久久亚洲精品成人影院| 麻豆乱淫一区二区| 亚洲欧美日韩无卡精品| 18禁动态无遮挡网站| 婷婷六月久久综合丁香| 亚洲av在线观看美女高潮| www.色视频.com| 国产综合懂色| 国产 亚洲一区二区三区 | 久久人人爽人人片av| 国产亚洲5aaaaa淫片| 国产亚洲一区二区精品| 欧美日韩在线观看h| 国产精品蜜桃在线观看| 久久久成人免费电影| 极品少妇高潮喷水抽搐| 国产亚洲91精品色在线| 真实男女啪啪啪动态图| 婷婷色综合大香蕉| 中文字幕亚洲精品专区| 亚洲国产精品专区欧美| 亚洲久久久久久中文字幕| 亚洲乱码一区二区免费版| 国产v大片淫在线免费观看| 亚洲欧洲国产日韩| 亚洲图色成人| 丰满人妻一区二区三区视频av| 日韩电影二区| 欧美性猛交╳xxx乱大交人| 久久久久久久久久成人| 69av精品久久久久久| 国产精品1区2区在线观看.| 久久久久九九精品影院| 午夜福利网站1000一区二区三区| av女优亚洲男人天堂| 婷婷色麻豆天堂久久| 超碰97精品在线观看| 26uuu在线亚洲综合色| 日韩欧美国产在线观看| 777米奇影视久久|