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

    基于諧波小波包和BSA優(yōu)化LS-SVM的銑刀磨損狀態(tài)識別研究

    2017-09-14 00:57:02董彩云張超勇孟磊磊肖鵬飛林文文
    中國機械工程 2017年17期
    關(guān)鍵詞:波包正確率特征提取

    董彩云 張超勇 孟磊磊 肖鵬飛 羅 敏 林文文

    1.華中科技大學(xué)機械學(xué)院數(shù)字制造裝備與技術(shù)國家重點實驗室,武漢,4300742.湖北汽車工業(yè)學(xué)院電氣與信息工程學(xué)院,十堰,4420023.寧波大學(xué)機械工程與力學(xué)學(xué)院,寧波,315211

    基于諧波小波包和BSA優(yōu)化LS-SVM的銑刀磨損狀態(tài)識別研究

    董彩云1張超勇1孟磊磊1肖鵬飛1羅 敏2林文文3

    1.華中科技大學(xué)機械學(xué)院數(shù)字制造裝備與技術(shù)國家重點實驗室,武漢,4300742.湖北汽車工業(yè)學(xué)院電氣與信息工程學(xué)院,十堰,4420023.寧波大學(xué)機械工程與力學(xué)學(xué)院,寧波,315211

    針對銑削刀具磨損狀態(tài)識別問題,提出諧波小波包和最小二乘支持向量機(LS-SVM)的狀態(tài)識別方法。為克服傳統(tǒng)小波包分解的頻帶交疊問題,采用諧波小波包提取不同磨損狀態(tài)下銑削力信號的各頻段信號能量,歸一化處理后,輸入LS-SVM多類分類器,實現(xiàn)銑削刀具磨損狀態(tài)的識別。針對LS-SVM的懲罰因子和核參數(shù)對模型識別精度影響較大的問題,提出回溯搜索算法(BSA)進行自動參數(shù)尋優(yōu)。實驗結(jié)果表明,諧波小波包比小波包在刀具磨損狀態(tài)特征提取時具有更好的識別效果。與粒子群算法進行比較,證明BSA優(yōu)化LS-SVM具有更高的識別精度。

    刀具磨損;諧波小波包;回溯搜索算法;最小二乘支持向量機

    0 引言

    刀具智能狀態(tài)監(jiān)測技術(shù)在先進制造生產(chǎn)中具有重要意義。刀具磨損在線監(jiān)測作為其主要內(nèi)容,可以保證產(chǎn)品質(zhì)量,提高生產(chǎn)效率,降低生產(chǎn)成本,縮短生產(chǎn)時間,避免因刀具失效造成的零件加工表面精度低、零件報廢、機床故障等問題,對加工過程自動化具有重要作用。國內(nèi)外有許多學(xué)者研究了刀具磨損智能監(jiān)測和刀具壽命預(yù)算方法,如TETI等[1]與ROTH等[2]概括了刀具磨損狀態(tài)監(jiān)測技術(shù)在傳感信號選擇、信號處理與模式識別等方面的發(fā)展現(xiàn)狀。

    目前刀具磨損狀態(tài)監(jiān)測的方法主要是間接法,通過引入模式識別的間接測量方法實現(xiàn)高準確率快速識別[3]。間接法監(jiān)測過程大致分為信息采集、特征提取、狀態(tài)識別、決策控制四個階段,其中特征提取和模式識別階段尤為關(guān)鍵。

    信號特征提取方法有時域法、頻域法、時頻域法、時間序列分析、分形特征等。時域和頻域法都不能兼顧信號的時域和頻域中的總體和局部變化,不能合理地處理非平穩(wěn)信號。時頻域法包括STFT、小波分析、HHT等。文獻[4]研究了小波分析技術(shù)在刀具磨損狀態(tài)識別中的應(yīng)用。文獻[5]研究了小波包特征提取技術(shù)在鈦合金車削刀具磨損狀態(tài)監(jiān)測中的應(yīng)用。文獻[6]利用分形理論對刀具磨損信號進行處理,發(fā)現(xiàn)切削力動態(tài)分量的分形維數(shù)和刀具磨損狀態(tài)具有很強的相關(guān)性。文獻[7]采用具有嚴格盒形譜特性的諧波小波包,能夠?qū)⑿盘柌唤化B且無遺漏地分解到不同的頻段,獲得信號在不同頻段上的特征信息。諧波小波具有優(yōu)良的濾波特性,在微弱信號、突變信號等非平穩(wěn)信號特征提取中具有很強的優(yōu)勢[8-9]。文獻[10-13]研究了諧波小波在振動信號、弱信號、碰摩故障以及旋轉(zhuǎn)機械故障等方面的應(yīng)用,均取得了很好的效果。

    在模式識別方面,文獻[14-16]研究了最小二乘支持向量機(least squares support vector machine,LS-SVM)在刀具磨損狀態(tài)監(jiān)測中的應(yīng)用,證明了該方法與其他算法相比,具有更高的識別精度,由于LS-SVM的模型參數(shù)對其識別精度有較大影響,不同參數(shù)下的刀具磨損量識別精度差別大,故如何選擇合適的懲罰因子和核參數(shù)在模式識別中非常關(guān)鍵。文獻[17]通過改進的粒子群優(yōu)化算法對最小二乘支持向量機的懲罰因子和徑向基核函數(shù)參數(shù)進行了優(yōu)化選擇,對車刀磨損量進行模式識別,通過與標準粒子群算法優(yōu)化的模型對比,發(fā)現(xiàn)識別誤差更小,收斂更快。文獻[18] 采用回溯搜索算法(backtracking search optimization algorithm,BSA)對多對不同Benchmark函數(shù)進行測試,與PSO、ABC等算法對比,證明了該算法操作簡單,在全局搜索能力和收斂速度上占優(yōu)。本文嘗試利用BSA與LS-SVM相結(jié)合的模式識別方法進行銑刀刀具磨損狀態(tài)識別,BSA用于優(yōu)化LS-SVM的懲罰因子和核參數(shù)。特征提取結(jié)果分為兩組,一組用于輸入LS-SVM中進行訓(xùn)練,得到優(yōu)化后的模式識別模型,另一組用于刀具磨損狀態(tài)檢測。

    1 刀具磨損機理

    1.1刀具磨損形式

    切削時由于切屑和刀具前面、加工表面與刀具后刀面之間會產(chǎn)生持續(xù)的熱量并存在摩擦,造成刀具的磨損。刀具磨損分為正常磨損(前刀面磨損、后刀面磨損、前后刀面同時磨損)和非正常磨損(脆性破損、塑性破損)。

    1.2刀具磨損過程

    刀具磨損量會隨切削時間變化而逐漸變大,磨損的快慢與加工條件、加工材料有關(guān)。刀具磨損過程在不同條件下基本類似,一般典型的刀具磨損曲線圖見圖1,刀具磨損過程一般分為初期磨損(Ⅰ)、正常磨損(Ⅱ)和急劇磨損(Ⅲ)三個階段。

    圖1 刀具磨損曲線Fig.1 Tool wear curve

    球頭銑削力模型通常都是基于刀具坐標系建立的,原點設(shè)定為刀尖點位置,有三個方向的切削分力(Fx,F(xiàn)y,F(xiàn)z)。切削力作為與刀具在加工時的磨損、破損狀態(tài)關(guān)系非常密切的物理特征,磨損量增大會導(dǎo)致切削力增大,刀具破損則引起切削力突變,所以對切削力變化的監(jiān)測能夠間接反映刀具磨損的狀態(tài)。對切削力造成影響的原因比較多,如:工件的材料、切削用量和刀具幾何參數(shù),但是刀具的材質(zhì)、是否選用切削液對刀具的磨損也有影響。

    2 諧波小波變換

    2.1諧波小波的概念

    諧波小波具有嚴格的矩形頻譜,頻域表達式為

    (1)

    (2)

    (3)

    諧波小波的時域表達式為

    ψ(t)=ψe(t)+iψo(t)=
    (exp(i4πt)-exp(i2πt))/(i2πt)

    (4)

    ψ(t)的實部和虛部如圖2所示,可以看出,諧波小波的實部是偶函數(shù),虛部是奇函數(shù),由信號分析理論可知,諧波小波具有“鎖定”信號相位的功能[19]。

    圖2 諧波小波的實部和虛部Fig.2 Real and imaginary parts of harmonic wavelet

    對諧波小波基進行伸縮、平移則可生成L2(R)空間的規(guī)范正交基,也即諧波小波函數(shù)族:

    ψ(t)=[exp(i4πt(2j-k))-
    exp(i2πt(2j-k))]/[i2πt(2j-k)]

    (5)

    其中,j為收縮因子,反映了諧波小波的分解層數(shù),j∈Z;k為平移因子。j=0時,在頻域內(nèi)諧波小波只在2π~4π內(nèi)有值,在第j層時只在2j+1π~2j+2π有值。諧波小波頻譜分布如圖3所示。

    圖3 諧波小波頻譜分布Fig.3 Harmonic wavelet spectrum distribution

    2.2廣義諧波小波

    由圖3可知,諧波小波同傳統(tǒng)小波一樣,頻率分辨率在低頻時高,在高頻時低。但是對實際信號進行分析時,同樣需要在高頻時具有較高分辨率,因此NEWLAND[7]構(gòu)造出了廣義諧波小波。

    取諧波小波每層的頻率段ω=[2πm,2πn],則諧波小波的頻域表達式為

    (6)

    其中,m,n滿足:

    m值取2j,j=0,1…,N-1,N為分解層數(shù)。

    對式(6)進行逆傅里葉變換可以得到時域表達式:

    ψm,n(t)=[exp(i2πtn)-exp(i2πtm)]/[i2πt(n-m)]

    (7)

    (8)

    式(8)就是頻率帶寬為2π(n-m)、時間中心為k/(n-m)的廣義諧波小波的一般形式。可以看出:不同頻帶所對應(yīng)的小波是相互正交的,相同頻帶(k≠0時)對應(yīng)的小波也是相互正交的。以諧波小波系數(shù)ψ(t)作為L2(R)的一組正交基,對信號f(t)進行諧波小波分解,就可以將信號分解到相互獨立的頻段,能量較弱的細節(jié)信號也能準確地顯現(xiàn)出來,有利于信號特征的提取。

    對于離散信號f(t),t=0,1,…,N-1,廣義諧波小波變換的時域表達式為

    (9)

    對式(9)進行離散傅里葉變換,則可得到頻域表達式:

    (10)

    2.3諧波小波包

    由式(8)可知,參數(shù)m、n決定了小波函數(shù)的帶寬和頻帶中心的位置,可以通過調(diào)節(jié)m、n的值來調(diào)節(jié)諧波小波的帶寬和中心頻率,從而在任意分解層上將信號分解到不同的頻段,實現(xiàn)諧波小波包變換,文獻[7,19]給出了諧波小波包變換實現(xiàn)的具體過程。

    3 回溯搜索算法優(yōu)化LS-SVM基本理論

    3.1LS-SVM基本原理

    支持向量機(support vector machine,SVM)是VAPNIK等[20-21]根據(jù)統(tǒng)計學(xué)中的結(jié)構(gòu)風(fēng)險最小化原則提出的一種機器學(xué)習(xí)方法,其基本原理是尋找一個使分類間隔最大化的最優(yōu)分類面。針對非線性問題,SVM引入核函數(shù),將輸入數(shù)據(jù)變換到高維空間,從而實現(xiàn)對問題的線性分類。SVM對小樣本、非線性、高維數(shù)等問題同樣可以獲得較好的分類效果。

    LS-SVM是標準SVM的變形,是將SVM求解二次規(guī)劃問題轉(zhuǎn)換成求解線性方程組,同時省去了不敏感損失系數(shù),使算法復(fù)雜程度大大降低[22]。用于分類的LS-SVM算法的詳細推導(dǎo)過程如下[23]。

    給定N個訓(xùn)練樣本{(x1,y1),(x2,y2)…,(xN,yN)},其中xi為n維的訓(xùn)練樣本輸入,yi為訓(xùn)練樣本輸出,i=1,2,…,N,yik∈{-1,+1}。LS-SVM算法的目標優(yōu)化函數(shù)為

    (11)

    i=1,2,…,N

    式中,φ(·)為核函數(shù),用來將輸入樣本映射到高維特征空間;ω為權(quán)矢量;ei為誤差;b為偏置量;γ為懲罰因子。

    LS-SVM優(yōu)化問題對應(yīng)的 Lagrange函數(shù):

    (12)

    其中,αi為拉格朗日乘子,αi≥0。

    (13)

    其中,Y=[y1y2…yN]T,1v=[11…1]T,α=[α1α2…αN]T,I為N×N階的單位矩陣,Z=[y1φ(x1)y2φ(x2) …yNφ(xN)]T,核函數(shù)矩陣Ω=ZZT,Ωkl=ykylφ(xk)Tφ(xl)=ykylK(xk,xl),核函數(shù)K(xk,xl)是滿足Mercer條件的對稱函數(shù)。

    LS-SVM的分類決策函數(shù)為

    (14)

    徑向基(RBF)核函數(shù)是最常用的核函數(shù)K(x,xi)=exp(-‖x-xi‖2/σ2),σ為徑向基函數(shù)的寬度。基于徑向基核函數(shù)的LS-SVM,需確定懲罰因子γ和核參數(shù)σ兩個參數(shù)。

    LS-SVM方法本身是針對二分類問題而設(shè)計的,但是平時要解決的問題大多是多分類問題,因此文獻[24]提出了將多分類問題轉(zhuǎn)化為多個二分類問題的解決方案。目前常用的將多分類問題轉(zhuǎn)化為多個二分類問題的方法主要有[25]:一對一方法(one versus one,1vs1)、一對多方法(one versus all,1vsA)、有向無環(huán)圖方法(directed acyclic graph,DAG)、糾錯輸出編碼方法(error-correlating output codes,ECOC)、最小輸出編碼(minimum output codes,MOC)等。

    3.2回溯搜索算法

    BSA是CIVICIOGLU[18]于2013年提出的一種新種群進化算法,該算法的流程和差分進化算法類似,但在變異操作(擾動策略)和交叉操作(混合策略)上,兩者有本質(zhì)的區(qū)別。BSA采用新型的擾動策略和混合策略,優(yōu)化效率大大提高,該算法僅有一個控制參數(shù)(混合比例參數(shù))使得操作更加簡單。

    BSA的算法流程有兩次選擇操作,可稱為選擇Ⅰ和選擇Ⅱ[26-27]。選擇Ⅰ用來選擇歷史種群,選擇Ⅱ用于更新種群。它的算法流程可分為初始化種群、選擇Ⅰ、變異、交叉、選擇Ⅱ五部分。

    (1)種群初始化。BSA的初始化包括兩個種群的初始化,分別為種群Pop和歷史種群oldPop,初始化按式(15)進行:

    (15)

    其中,i=1,2,…,popsize;j=1,2,…,D,popsize是種群規(guī)模,D是個體長度(問題維數(shù)),lowj、upj分別是第j維分量的最小值和最大值,U(·)表示均勻分布,Popi,j、oldPopi,j是(lowj,upj)上服從均勻分布的隨機數(shù),兩次初始化操作相互獨立。為了避免第一次執(zhí)行選擇Ⅰ操作時,oldPop值為空,這里對oldPop也進行了初始化,以保證算法的可行性。

    (2)選擇Ⅰ。BSA的選擇Ⅰ用來為每一次迭代選擇一個新的歷史種群oldPop,其算法思想如下:

    ifa

    (16)

    式中,a、b為(0,1)上服從均勻分布的隨機數(shù)。

    (3)變異。在oldPop的值確定后,對oldPop中的個體進行隨機排序,生成種群oldPop'。利用下式進行變異:

    Mutant=Pop+F(oldPop'-Pop)

    (17)

    這里F是變異尺度系數(shù),用以控制變異的程度;F=3rand(n),rand(n)是服從標準正態(tài)分布的隨機數(shù)。

    (4)交叉。BSA使用了一種新的種群交叉策略,通過設(shè)置混合比例參數(shù)mr來控制種群間交叉的粒子個數(shù)[28],具體表達式如下:

    (18)

    其中,mapi,j為N×D的二元整數(shù)矩陣中的元素,元素初始賦值全為1,具體計算表達式如下:

    (19)

    (5)選擇Ⅱ。通過貪婪選擇機制,比較種群T與種群Pop對應(yīng)位置的適應(yīng)度,用T中所有適應(yīng)度比Pop較好的個體代替Pop的對應(yīng)位置個體,完成Pop的更新,一次迭代完成。選擇Ⅱ如下式所示:

    (20)

    記錄當前最優(yōu)解和對應(yīng)的解向量,重復(fù)步驟(2)~步驟(5),直到滿足循環(huán)終止條件,最后輸出最優(yōu)解。

    3.3回溯搜索算法優(yōu)化LS-SVM

    本文提出運用BSA優(yōu)化算法對LS-SVM的懲罰因子γ和核參數(shù)σ進行自動選擇,以達到最佳組合值,使識別模型準確率達到最高。

    假設(shè)已建立了系統(tǒng)的樣本集{(x1,y1),(x2,y2),…,(xN,yN)},BSA優(yōu)化LS-SVM過程中的LS-SVM模型的優(yōu)劣采用最小分類錯誤率(minimal classification error rate,MCER)來衡量:

    (20)

    式中,m1為被錯誤分類的樣本數(shù);N為訓(xùn)練樣本數(shù)量。

    如果在BSA優(yōu)化過程中直接采用MCER作為適應(yīng)度值,則可能會造成識別模型對訓(xùn)練樣本的過擬合,但對于測試樣本的泛化能力很差,從而失去了模型的實用性。為了提高識別模型的泛化能力,本文在參數(shù)優(yōu)化時采用了K-折交叉驗證(K-fold cross validation,K-CV),即在BSA的每次迭代過程中,對當前γ和σ所確定的LS-SVM進行交叉驗證,將交叉驗證的平均MCER作為BSA迭代的適應(yīng)度值。

    K-折交叉驗證,即將訓(xùn)練樣本分成K組(一般為均分),每個子樣本集分別作為一次驗證集,則可得到K個模型,用K個模型的平均分類錯誤率作為此K-CV下分類器的性能指標。K一般從3開始取。使用K-CV交叉驗證方法,可以有效避免識別模型的過學(xué)習(xí)以及欠學(xué)習(xí)。

    基于BSA優(yōu)化LS-SVM的刀具磨損狀態(tài)識別步驟如下:①提取特征向量,利用諧波小波包將銑削力信號分解到各個頻段,提取各頻段能量特征,并進行歸一化處理,作為LS-SVM模型輸入,刀具磨損狀態(tài)類別序號作為輸出。②初始化BSA種群,將LS-SVM的懲罰因子γ和核參數(shù)平方σ2作為每個個體的坐標值,用訓(xùn)練樣本訓(xùn)練LS-SVM,并進行K-折交叉驗證,計算種群的適應(yīng)度f。③按照BSA流程進行迭代尋優(yōu)。對每個個體,將適應(yīng)度f(xi)與自身最優(yōu)值進行比較,更新個體最優(yōu)值;將個體最優(yōu)值與種群(全局)最優(yōu)值進行比較,更新種群最優(yōu)值。④判斷是否滿足結(jié)束條件,若滿足則輸出最優(yōu)化參數(shù)γ和σ2,若不滿足則跳轉(zhuǎn)到步驟②。結(jié)束條件一般設(shè)置為算法達到最大迭代次數(shù)、適應(yīng)度小于給定精度等,本文以算法達到最大迭代次數(shù)為結(jié)束條件。⑤選出最優(yōu)化參數(shù)γ和σ2,建立LS-SVM刀具磨損狀態(tài)識別模型。⑥刀具磨損狀態(tài)識別驗證,利用測試樣本檢驗?zāi)P妥R別能力。

    4 刀具磨損監(jiān)測實驗及特征提取

    4.1刀具磨損監(jiān)測實驗

    試驗數(shù)據(jù)來源于美國紐約預(yù)測與健康管理學(xué)會(Prognostic and Health Management Society,PHM)2010年高速數(shù)控機床刀具健康預(yù)測競賽的開放數(shù)據(jù)。實驗條件如表1所示。

    表1 實驗條件

    試驗在上述切削條件下重復(fù)進行6次全壽命周期試驗。端面銑削材料為正方形,每次走刀端面銑的長度為108 mm且每次走刀時間相等,每次走刀后測量刀具的后刀面磨損量。試驗監(jiān)測數(shù)據(jù)有x、y、z三向銑削力信號,x、y、z三向銑削振動信號以及聲發(fā)射均方根值。

    本文選取第4組試驗的x方向銑削力信號進行分析。本組試驗總共315次走刀,每次走刀采樣點數(shù)在220 000以上,為避免切入、切出的影響,本文提取每次走刀的50 001~100 000共50 000個數(shù)據(jù)點進行分析。

    圖4所示為第4組試驗刀具磨損變化曲線,圖中共有flute1、flute2、flute3、average四條曲線,分別表示x、y、z方向切削刃和其平均磨損情況,圖中橫坐標為走刀次數(shù),縱坐標為后刀面磨損量。本文采用三個切削刃磨損量的均值來衡量球頭銑刀的磨損情況。可以看出,前30次走刀為初期磨損階段,30~225次走刀為正常磨損階段,225次走刀后為急劇磨損階段。

    圖4 球頭銑刀磨損變化曲線Fig.4 Ball milling cutter wear curve

    根據(jù)圖4的銑削磨損情況,本文采用K均值聚類分析方法將磨損狀態(tài)分為5類,狀態(tài)1:1~30次走刀;狀態(tài)2:31~130次走刀;狀態(tài)3:131~200次走刀;狀態(tài)4:201~260次走刀;狀態(tài)5:261~315次走刀。

    4.2諧波小波包頻帶能量敏感特征提取

    (a)第15次走刀 (b)第80次走刀

    (c)第165次走刀 (d)第240次走刀

    (e)第270次走刀圖5 走刀次數(shù)為15、80、165、240、270的x方向銑削力信號的原始信號Fig.5 Original signal of the x-direction milling force signal for 15, 80, 165, 240,270 feeding times

    (a)第15次走刀 (b)第80次走刀

    (c)第165次走刀 (d)第240次走刀

    (e)第270次走刀圖6 走刀次數(shù)為15、80、165、240、270的x方向銑削力信號的頻譜Fig.6 The frequency of the x-direction milling force signal for 15, 80, 165, 240,270 feeding times

    對前8個頻段分別求信號的能量:

    (22)

    j=1,2,…,8

    式中,n為采樣點數(shù)。

    以這8個頻段分別對應(yīng)的能量值為元素可以構(gòu)成一個特征向量:T=(E1,E2,…,E8),由于各頻段能量數(shù)值大小相差比較大,對特征向量進行歸一化處理:

    (23)

    其中,E為總能量,T′=(E1/E,E2/E,…,E8/E),歸一化特征向量的各元素為 0~1之間的數(shù)。

    圖7所示為對走刀次數(shù)為15、80、165、240、270的x方向銑削力信號分別進行 6層諧波小波包分解后得到的頻段能量分布圖。從圖中可以發(fā)現(xiàn),在不同磨損狀態(tài)下的力信號在同一頻段的能量存在差異,因此可將其作為LS-SVM的輸入向量,對刀具磨損狀態(tài)進行識別。

    (a)第15次走刀 (b)第80次走刀

    (c)第165次走刀 (d)第240次走刀

    (e)第270次走刀圖7 諧波小波包歸一化頻段能量分布圖Fig.7 Harmonic wavelet packet normalized band energy distribution map

    為了展示諧波小波包特征提取相對于小波包的優(yōu)越性,將其與常規(guī)的小波包dB9的特征提取結(jié)果輸入改進LS-SVM中進行狀態(tài)識別,比較識別結(jié)果。圖8所示為對走刀次數(shù)為15、80、165、240、270的x方向銑削力信號分別進行6層小波包分解后得到的前8個歸一化頻段能量分布圖。圖8中,頻段3、5、6、7、8數(shù)值很小,且變化不大,為減少計算量,可忽略不計。最終特征向量T′=(E1/E,E2/E,E4/E)。

    (a)第15次走刀 (b)第80次走刀

    (c)第165次走刀 (d)第240次走刀

    (e)第270次走刀圖8 小波包歸一化頻段能量分布圖Fig.8 Wavelet packet normalized band energy distribution map

    圖9為頻段1、頻段2、頻段4的歸一化能量分布隨走刀數(shù)的變化趨勢圖??梢钥闯觯C波小波包能量分布和小波包能量分布雖然變化趨勢相似,仍存在較大差異,且諧波小波包能量分布波動小。

    4.3模型的訓(xùn)練和測試

    本實驗總共315次走刀,每次走刀取一個樣本,總共315個樣本,劃分為5個狀態(tài)(同4.1節(jié)),每種狀態(tài)按80%為訓(xùn)練樣本,20%為測試樣本劃分。狀態(tài)1:訓(xùn)練樣本24個,測試樣本6個;狀態(tài)2:訓(xùn)練樣本80個,測試樣本20個;狀態(tài)3:訓(xùn)練樣本56個,測試樣本14個;狀態(tài)4:訓(xùn)練樣本48個,測試樣本12個;狀態(tài)5:訓(xùn)練樣本44個,測試樣本11個;則訓(xùn)練樣本總共為252個,測試樣本總共為63個。

    (a)諧波小波包能量分布

    (b)小波包能量分布圖9 頻段能量變化趨勢Fig.9 Band energy trends

    這一5分類問題,采用荷語魯汶天主教大學(xué)編寫的LSSVMlabv1.8版本工具箱,運行條件為:Dell臺式機,Windows 7操作系統(tǒng),雙核英特爾i5處理器,8G內(nèi)存,MATLAB 2014a。通過多次試驗對比1vs1、1vsA、ECOC、MOC四種方法,發(fā)現(xiàn)四種方法的分類準確率差異很小,而運行速度相差比較大。由于MOC采用最小數(shù)目的兩類分類器,只需要構(gòu)建3個兩類分類器,運行速度最快,因此,本文采用基于MOC編碼的多分類LS-SVM來對刀具磨損狀態(tài)進行識別。

    確定了訓(xùn)練數(shù)據(jù)后,進行回溯搜索算法的初始化工作。進化代數(shù)為30;種群規(guī)模為20;混合比例參數(shù)mr=1;γ,σ2的最大值和最小值都為1000和0.01。為了避免過學(xué)習(xí)以及欠學(xué)習(xí)狀態(tài)的發(fā)生,提高模型的泛化能力,采用3-CV交叉驗證方法。

    在參數(shù)一致的條件下,分別將諧波小波包特征提取結(jié)果與小波包特征提取結(jié)果輸入優(yōu)化后的LS-SVM中進行狀態(tài)識別,圖10為某次迭代過程,兩種特征提取方法得到的種群的最佳正確率與平均正確率隨進化代數(shù)的變化圖。

    (a)諧波小波包

    (b)小波包圖10 不同特征提取方法的最優(yōu)正確率與平均正確率變化圖Fig.10 The optimal accuracy and average correct rate with different feature extractians

    從圖10中可以看出,諧波小波包分析的平均正確率進化代數(shù)為16時達到最優(yōu),而小波包在20代之后達到最優(yōu),基于諧波小波包特征提取方法的狀態(tài)識別結(jié)果能更快收斂達到最大平均正確率。但是一次迭代過程可能存在隨機性,由于適應(yīng)度函數(shù)的選擇,達到最優(yōu)值的解可能不止一個,同時由于LS-SVM的最佳γ,σ2的組合可能也不止一個,因此運行10次,每次迭代結(jié)束時對比結(jié)果。由表2和表3可知,10次運行結(jié)果中,σ2=0.01,γ∈[20,30],測試樣本正確率存在區(qū)別。通過測試實驗結(jié)果比較發(fā)現(xiàn),表2和表3中訓(xùn)練樣本最佳正確率都可以達到100%,而表2中的訓(xùn)練樣本平均正確率和測試樣本正確率都比表3的高,且表2中正確率都在95%以上,說明采用諧波小波包提取特征向量的方法比小波包的正確率更高。

    表2 諧波小波包特征提取的識別結(jié)果

    表3 小波包特征提取的識別結(jié)果

    為了說明本文采用的回溯搜索算法在優(yōu)化LS-SVM參數(shù)上的優(yōu)越性,與粒子群優(yōu)化算法進行對比。粒子群算法的參數(shù)、進化代數(shù)、種群規(guī)模、γ及σ2的最大化值、最小值與回溯搜索算法相同,迭代速度范圍為[-0.2×1000,0.2×1000]。采用線性遞減權(quán)重,ω的最大值和最小值分別為0.8和0.4,學(xué)習(xí)因子c1=c2=2,同樣采用3-CV交叉驗證方法。圖11為某次迭代過程,兩種算法的種群的最優(yōu)正確率與平均正確率隨進化代數(shù)的變化圖。從圖中可以看出,在測試樣本分類的平均正確率上,BSA優(yōu)化LS-SVM的狀態(tài)識別模型可以達到100%,明顯比PSO優(yōu)化LS-SVM的模型效果更好,且收斂迅速。

    (a)BSA優(yōu)化

    (b)PSO優(yōu)化圖11 不同優(yōu)化算法的最優(yōu)正確率與平均正確率變化圖Fig.11 The optimal accuracy and average correct rate with different optimization algorithm

    由于適應(yīng)度函數(shù)的選擇,達到最優(yōu)值的解可能不止一個,同時由于LS-SVM的最佳γ,σ2的組合可能也不止一個,本文運行10次,每次迭代結(jié)束時對比結(jié)果。由表4和表5可知,10次運行結(jié)果中,σ2=0.01,γ∈[20,30],表4的測試樣本正確率都在90%以上。經(jīng)過10次試驗對比,兩種優(yōu)化方法的群體最佳正確率都在20代以內(nèi)達到最優(yōu),但是平均正確率BSA比PSO效果好得多。說明BSA在LS-SVM參數(shù)尋優(yōu)方面更有優(yōu)勢,驗證了采用BSA優(yōu)化LS-SVM的狀態(tài)模式識別方法具有較高的正確率。

    表4 BSA優(yōu)化LS-SVM結(jié)果

    表5 PSO優(yōu)化LS-SVM結(jié)果

    5 結(jié)論

    本文提出了通過諧波小波包對銑削力信號進行分頻段能量特征提取,并采用回溯搜索算法對LS-SVM參數(shù)進行尋優(yōu),對銑削刀具磨損狀態(tài)進行識別的方法。實驗結(jié)果表明,對于采用銑削力信號進行刀具磨損狀態(tài)識別,利用諧波小波包提取不同頻段能量,將歸一化處理后樣本作為LS-SVM的輸入樣本相比小波包提取特征具有較高的識別正確率。在LS-SVM參數(shù)尋優(yōu)方面,通過對比回溯搜索算法和經(jīng)典的粒子群優(yōu)化算法,實驗測試結(jié)果顯示了回溯優(yōu)化算法在LS-SVM參數(shù)尋優(yōu)方面的優(yōu)越性,以及對于銑削刀具磨損狀態(tài)識別,將諧波小波包分解和回溯搜索算法優(yōu)化的LS-SVM相結(jié)合的方法可以獲得良好的效果。

    [1] TETI R, JEMIELNIAK K, O’DONNELL G, et al. Advanced Monitoring of Machining Operations[J]. CIRP Annals—Manufacturing Technology,2010,59(2):717-739.

    [2] ROTH J T, DJURDJANOVIC D, YANG X, et al. Quality and Inspection of Machining Operations: Tool Condition Monitoring[J]. Journal of Manufacturing Science & Engineering, 2010, 132(4):575-590.

    [3] SNR D E D. Sensor Signals for Tool-wear Monitoring in Metal Cutting Operations—a Review of Methods[J]. International Journal of Machine Tools & Manufacture, 2000, 40(8):1073-1098.

    [4] 叢龍慧,韓玉杰.小波分析在刀具磨損狀態(tài)檢測中的應(yīng)用[J].林業(yè)機械與木工設(shè)備,2009,37(4):51-52. CONG Longhui, HAN Yujie. Application of Wavelet Analysis in Tool Wear Status Detection[J]. Forestry Machinery & Woodworking Equipment, 2009,37(4):51-52.

    [5] 陳侃,傅攀,李威霖,等.鈦合金車削加工過程中刀具磨損狀態(tài)監(jiān)測的小波包子帶能量變換特征提取新方法[J].組合機床與自動化加工技術(shù),2011(1):35-38. CHEN Kan, FU Pan, LI Weilin, et al. The Tool Wear Monitoring Based on SUB-WPA Feature Selection in Process of Titanium Alloy Machining[J]. Modular Machine Tool & Automatic Manufacturing Technique, 2011(1):35-38.

    [6] 林穎, 劉亞俊, 陳忠. 基于分形理論和神經(jīng)網(wǎng)絡(luò)的刀具磨損監(jiān)測[J]. 中國機械工程, 2004, 15(16):20-22. LIN Ying, LIU Yajun, CHEN Zhong. Tool Wear Detection Using Fractal Theory and Neural Network[J].China Mechanical Engineering, 2004, 15(16):20-22.

    [7] NEWLAND D E. Harmonic Wavelet Analysis[J]. Proceedings of the Royal Society A: Mathematical Physical & Engineering Sciences, 1993, 443:203-225.

    [8] NEWLAND D E. Wavelet Analysis of Vibration: Part 1—Theory[J].Journal of Vibration and Acoustics,1994,116(4):409-416.

    [9] NEWLAND D E. Wavelet Analysis of Vibration: Part 2—Wavelet Maps[J].Journal of Vibration and Acoustics,1994,116(4):417-425.

    [10] 高強,何正嘉.諧波小波及其時頻剖面圖在旋轉(zhuǎn)機械診斷中的應(yīng)用[J].西安交通大學(xué)學(xué)報,2000,34(9):62-66. GAO Qiang, HE Zhengjia. Application of Harmonic Wavelet and Its Time-frequency Profile Plot to Diagnosis of Rotating Machinery[J].Journal of Xi’an Jiaotong University,2000,34(9):62-66.

    [11] 陳潔,沈遠彤,李宏偉.基于諧波小波的弱信號提取[J].信號處理,2005(增刊1):40-43. CHEN Jie, SHEN Yuantong, LI Hongwei. Applying Harmonic Wavelet Extracting Weak Signal[J].Signal Processing,2005(S1):40-43.

    [12] 萬方義,許慶余.裂紋轉(zhuǎn)子—滑動軸承系統(tǒng)動靜件碰摩故障的諧波小波分析[C]//2002年全國振動工程及應(yīng)用學(xué)術(shù)會議論文集.北京,2002:111-116. WAN Fangyi, XU Qingyu. Harmonic Wavelet Analysis on Cracked Rotor Sliding Bearing System with Rotor and Stator Rubbing[C]//Proceedings of the 2002 National Conference on Vibration Engineering and Application.Beijing,2002:111-116.

    [13] 肖忠會.諧波小波分析及其在旋轉(zhuǎn)機械信號分析中的應(yīng)用[J].風(fēng)機技術(shù),2001(1):43-46. XIAO Zhonghui. Harmonic Wavelet Analysis and Its Application to Analyzing Signals of Rotating Machinery[J].Compressor Blower & Fan Technology,2001(1):43-46.

    [14] 關(guān)山,王龍山,聶鵬.基于EMD與LS-SVM的刀具磨損識別方法[J].北京航空航天大學(xué)學(xué)報,2011,37(2):144-148. GUAN Shan, WANG Longshan, Nie Peng. Identification Method of Tool Wear Based on Empirical Mode Decomposition and Least Squares Support Vector Machine[J].Journal of Beijing University of Aeronautics and Astronautics,2011,37(2):144-148.

    [15] SHI D F, GINDY N N. Tool Wear Predictive Model Based on Least Squares Support Vector Machines[J].Mechanical Systems and Signal Processing,2007,21(4):1799-1814.

    [16] WIDODO A, YANG B S. Support Vector Machine in Machine Condition Monitoring and Fault Diagnosis[J].Mechanical Systems and Signal Processing,2007,21(6):2560-2574.

    [17] 李威霖,傅攀,張爾卿.基于粒子群優(yōu)化LS-SVM的車刀磨損量識別技術(shù)研究[J].計算機應(yīng)用研究,2014,31(4):1094-1097. LI Weilin, FU Pan, ZHANG Erqing. Application of Particle Swarm Optimization-Least Square Support Vector Machine in Tool Wear Monitoring[J].Application Research of Computers,2014,31(4):1094-1097.

    [18] CIVICIOGLU P. Backtracking Search Optimization Algorithm for Numerical Optimization Problems[J].Applied Mathematics & Computation,2013,219(15):8121-8144.

    [19] 張文斌,周曉軍,林勇,等.基于諧波小波包方法的旋轉(zhuǎn)機械故障信號提取[J].振動與沖擊,2009,28(3):87-89. ZHANG Wenbin, ZHOU Xiaojun, LIN Yong, et al. Harmonic Wavelet Packet Method Used to Extract Fault Signal of a Rotation Machinery[J].Journal of Vibration and Shock,2009,28(3):87-89.

    [20] VLADIMIR N, VAPNIK V. Statistical Learning Theory[M].New York:Wiley,1998.

    [21] VAPNIK V.The Nature of Statistical Learning Theory[M].Dordrecht: Springer Science & Business Media,2013.

    [22] SUYKENS J A K, GESTEL T V, Brabanter J D, et al. Least Squares Support Vector Machines[J]. Euphytica, 2002, 2(2):1599-1604.

    [23] SUYKENS J A K. Least Squares Support Vector Machine Classifiers[J].Neural Processing Letters,1999,9(3):293-300.

    [24] HSU C W, LIN C J. A Comparison of Methods for Multiclass Support Vector Machines[J]. IEEE Transactions on Neural Networks,2002,13(2):415-425.

    [25] 董曉睿.基于支持向量機的多分類模型的研究和設(shè)計[D].南昌:南昌大學(xué),2013. DONG Xiaorui. Research and Design of Multiclassification Model Based on Support Vector Machine[D]. Nanchang:Nanchang University,2013.

    [26] 王曉娟,劉三陽,田文凱.帶高效變異尺度系數(shù)和貪婪交叉策略的回溯搜索優(yōu)化算法[J].計算機應(yīng)用,2014,34(9):2543-2546. WANG Xiaojuan, LIU Sanyang, TIAN Wenkai. Improved Backtracking Search Optimization Algorithm with New Effective Mutation Scale Factor and Greedy Crossover Strategy[J].Journal of Computer Applications,2014,34(9):2543-2546.

    [27] 田文凱.一種高效自學(xué)習(xí)性回溯搜索優(yōu)化算法[J].電子科技,2015,28(2):41-45. TIAN Wenkai. Effective Self-learning Backtracking Search Optimization Algorithm[J].Electronic Science and Technology,2015,28(2):41-45.

    [28] 李牧東,趙輝,翁興偉.具有廣泛學(xué)習(xí)策略的回溯搜索優(yōu)化算法[J].系統(tǒng)工程與電子技術(shù),2015,37(4):958-963. LI Mudong, ZHAO Hui, WENG Xingwei. Backtracking Search Optimization Algorithm with Comprehensive Learning Strategy[J].Journal of Systems Engineering and Electronics,2015,37(4):958-963.

    (編輯王旻玥)

    StateRecognitionofMillingToolWearsBasedonHarmonicWaveletPacketandBSAOptimizationLS-SVM

    DONG Caiyun1ZHANG Chaoyong1MENG Leilei1XIAO Pengfei1LUO Min2LIN Wenwen3

    1.State Key Lab of Digital Manufacturing Equipment & Technology,Huazhong University of Science and Technology,Wuhan,430074 2.School of Electrical and Information Engineering,Hubei Automotive Industries Institute,Shiyan, Hubei,442002 3.School of Mechanical Engineering and Mechanics,Ningbo University,Ningbo, Zhejiang,315211

    Aiming at the problems of milling tool wear state recognitions, a state recognition method was proposed based on harmonic wavelet packet and LS-SVM. To overcome the band overlapping problems in traditional wavelet packet decompositions, the milling force signal energies of each bands were extracted in different wear states by harmonic wavelet packet, which were brought in multi-class LS-SVM classifier after normalizing, then the classification recognition of different cutting tool states was achieved. BSA was proposed to search the optimal values of the kernel functional parameters and error penalty factors which affected the precision of the LS-SVM significantly. Experimental results show that harmonic wavelet packet is more effective and feasible than wavelet packet, and the proposed milling tool wear recognition method has higher accuracy.

    tool wear; harmonic wavelet packet; backtracking search algorithm(BSA); least squares support vector machine(LS-SVM)

    2016-09-13

    國家自然科學(xué)基金資助項目(51575211,51421062);國家自然科學(xué)基金國際(地區(qū))合作與交流項目(51561125002);湖北省自然科學(xué)基金資助項目(2014CFB348);高等學(xué)校學(xué)科創(chuàng)新引智計劃資助項目(B16019)

    TH164

    10.3969/j.issn.1004-132X.2017.17.011

    董彩云,女,1993年生。華中科技大學(xué)機械科學(xué)與工程學(xué)院碩士研究生。主要研究方向為刀具磨損狀態(tài)識別。張超勇(通信作者),男,1973年生。華中科技大學(xué)機械科學(xué)與工程學(xué)院副教授。E-mail:zcyhust@hust.edu.cn。孟磊磊,男,1991年生。華中科技大學(xué)機械科學(xué)與工程學(xué)院博士研究生。肖鵬飛,男,1993年生。華中科技大學(xué)機械科學(xué)與工程學(xué)院碩士研究生。羅敏,男,1967年生。湖北汽車工業(yè)學(xué)院電氣與信息工程學(xué)院教授。林文文,男,1989年生。寧波大學(xué)機械工程與力學(xué)學(xué)院講師。

    猜你喜歡
    波包正確率特征提取
    門診分診服務(wù)態(tài)度與正確率對護患關(guān)系的影響
    基于小波包Tsallis熵和RVM的模擬電路故障診斷
    基于Daubechies(dbN)的飛行器音頻特征提取
    電子制作(2018年19期)2018-11-14 02:37:08
    Bagging RCSP腦電特征提取算法
    生意
    品管圈活動在提高介入手術(shù)安全核查正確率中的應(yīng)用
    天津護理(2016年3期)2016-12-01 05:40:01
    生意
    故事會(2016年15期)2016-08-23 13:48:41
    基于小波包變換的電力系統(tǒng)諧波分析
    小波包理論與圖像小波包分解
    基于MED和循環(huán)域解調(diào)的多故障特征提取
    我的女老师完整版在线观看| 波野结衣二区三区在线| 成人毛片a级毛片在线播放| 视频区图区小说| 欧美一区二区亚洲| 丝袜脚勾引网站| 亚洲国产精品国产精品| 国产白丝娇喘喷水9色精品| 久久99蜜桃精品久久| 2018国产大陆天天弄谢| 18禁在线无遮挡免费观看视频| 国产无遮挡羞羞视频在线观看| 赤兔流量卡办理| 三级国产精品欧美在线观看| 丰满人妻一区二区三区视频av| 久久国产亚洲av麻豆专区| 三级国产精品片| videos熟女内射| 网址你懂的国产日韩在线| 免费观看a级毛片全部| 青春草亚洲视频在线观看| 成年人午夜在线观看视频| 极品少妇高潮喷水抽搐| 色网站视频免费| 久久精品国产自在天天线| 观看美女的网站| 久久久久性生活片| 亚洲无线观看免费| 亚洲天堂av无毛| 青春草亚洲视频在线观看| 熟女人妻精品中文字幕| 蜜桃在线观看..| 日韩人妻高清精品专区| 精品久久久久久久久av| 亚洲欧美一区二区三区国产| 人妻 亚洲 视频| av在线蜜桃| 免费不卡的大黄色大毛片视频在线观看| 国产精品久久久久久av不卡| 亚洲国产成人一精品久久久| 十八禁网站网址无遮挡 | 黑丝袜美女国产一区| 五月伊人婷婷丁香| 一二三四中文在线观看免费高清| 直男gayav资源| 伊人久久精品亚洲午夜| 日韩精品有码人妻一区| 18禁在线无遮挡免费观看视频| 久久久久久人妻| 人人妻人人爽人人添夜夜欢视频 | 身体一侧抽搐| 少妇的逼好多水| 免费观看的影片在线观看| 观看免费一级毛片| 久久精品国产亚洲av涩爱| 国产白丝娇喘喷水9色精品| 久久久久久久久久成人| 亚洲精品亚洲一区二区| av又黄又爽大尺度在线免费看| 久久ye,这里只有精品| 国产 精品1| 精品一区二区三卡| 成人毛片60女人毛片免费| 亚洲精品色激情综合| av国产久精品久网站免费入址| 亚洲精品乱码久久久v下载方式| 伊人久久精品亚洲午夜| 亚洲精品,欧美精品| 少妇的逼好多水| videossex国产| 精品亚洲乱码少妇综合久久| 亚洲av日韩在线播放| 大码成人一级视频| 老司机影院成人| 国产亚洲av片在线观看秒播厂| 成人亚洲欧美一区二区av| 搡女人真爽免费视频火全软件| 国产精品久久久久久av不卡| 久久国产精品男人的天堂亚洲 | 97热精品久久久久久| 国产在线一区二区三区精| 久久女婷五月综合色啪小说| 大码成人一级视频| 国产黄色视频一区二区在线观看| 国产欧美日韩一区二区三区在线 | 欧美丝袜亚洲另类| 99久久精品一区二区三区| 欧美精品人与动牲交sv欧美| 偷拍熟女少妇极品色| 男人爽女人下面视频在线观看| 亚洲三级黄色毛片| 亚洲精品久久久久久婷婷小说| 久久ye,这里只有精品| 色视频www国产| 人妻一区二区av| 国产黄频视频在线观看| 人人妻人人添人人爽欧美一区卜 | 国产有黄有色有爽视频| 国产男女内射视频| 国产无遮挡羞羞视频在线观看| 亚洲国产日韩一区二区| 少妇人妻久久综合中文| 久久国内精品自在自线图片| 国产午夜精品一二区理论片| 免费人妻精品一区二区三区视频| 成人影院久久| 欧美成人精品欧美一级黄| 亚洲国产精品999| 少妇人妻一区二区三区视频| 国产精品久久久久久av不卡| 久久精品熟女亚洲av麻豆精品| 国产高潮美女av| 一边亲一边摸免费视频| 性色av一级| 99国产精品免费福利视频| 精品国产乱码久久久久久小说| 嫩草影院入口| 老女人水多毛片| 国产有黄有色有爽视频| 国产精品精品国产色婷婷| 国产一级毛片在线| 精品酒店卫生间| 亚洲伊人久久精品综合| 乱系列少妇在线播放| 91精品一卡2卡3卡4卡| 国产一区二区三区av在线| 亚洲欧美日韩卡通动漫| 国产大屁股一区二区在线视频| 免费观看av网站的网址| 搡老乐熟女国产| 久久精品国产亚洲网站| 国产欧美亚洲国产| 啦啦啦在线观看免费高清www| 国产精品99久久99久久久不卡 | 亚洲美女搞黄在线观看| 亚洲av欧美aⅴ国产| 激情 狠狠 欧美| 香蕉精品网在线| 久久久久久久久久人人人人人人| 夫妻午夜视频| 欧美丝袜亚洲另类| 亚洲国产高清在线一区二区三| 少妇的逼好多水| 美女高潮的动态| 女性生殖器流出的白浆| av又黄又爽大尺度在线免费看| 狂野欧美激情性bbbbbb| h日本视频在线播放| 亚洲精品色激情综合| 久久久久国产精品人妻一区二区| 国产精品一区二区在线观看99| 亚州av有码| 午夜激情福利司机影院| 久久久久久久久久人人人人人人| 国产亚洲午夜精品一区二区久久| 嘟嘟电影网在线观看| 欧美最新免费一区二区三区| 精品一区在线观看国产| 久久精品国产a三级三级三级| 亚洲欧美精品自产自拍| 日本爱情动作片www.在线观看| 直男gayav资源| 亚洲av日韩在线播放| 丰满迷人的少妇在线观看| 中文字幕免费在线视频6| 精品国产露脸久久av麻豆| 大片电影免费在线观看免费| 欧美极品一区二区三区四区| 国产av精品麻豆| xxx大片免费视频| 韩国av在线不卡| 久久久久国产网址| 日本猛色少妇xxxxx猛交久久| 建设人人有责人人尽责人人享有的 | 国产爽快片一区二区三区| av.在线天堂| 国产亚洲最大av| 极品教师在线视频| 欧美日韩一区二区视频在线观看视频在线| 18禁在线播放成人免费| 亚洲国产精品国产精品| 欧美精品一区二区免费开放| 亚洲国产精品成人久久小说| 亚洲中文av在线| 一区二区av电影网| 国产av国产精品国产| 婷婷色av中文字幕| 交换朋友夫妻互换小说| 国产黄片视频在线免费观看| 秋霞伦理黄片| 国产成人91sexporn| 狂野欧美白嫩少妇大欣赏| 国产在线免费精品| 国产男人的电影天堂91| 欧美成人午夜免费资源| 亚洲真实伦在线观看| 久久午夜福利片| 日韩欧美一区视频在线观看 | 免费观看性生交大片5| 九九爱精品视频在线观看| 91久久精品电影网| 超碰av人人做人人爽久久| 97在线人人人人妻| 国产精品无大码| 亚洲综合精品二区| 十八禁网站网址无遮挡 | 婷婷色综合www| 久久亚洲国产成人精品v| 菩萨蛮人人尽说江南好唐韦庄| 日韩 亚洲 欧美在线| 黄色配什么色好看| 夫妻午夜视频| 国产精品麻豆人妻色哟哟久久| 久久久午夜欧美精品| 美女内射精品一级片tv| 国产毛片在线视频| 日韩av免费高清视频| 香蕉精品网在线| 男女免费视频国产| 欧美成人a在线观看| 中文字幕人妻熟人妻熟丝袜美| 91精品一卡2卡3卡4卡| 免费av中文字幕在线| 一级二级三级毛片免费看| 亚洲欧美中文字幕日韩二区| 又黄又爽又刺激的免费视频.| 久久av网站| 久久精品人妻少妇| av又黄又爽大尺度在线免费看| 麻豆成人午夜福利视频| 成人一区二区视频在线观看| 永久网站在线| 夫妻午夜视频| 干丝袜人妻中文字幕| 亚洲欧美日韩东京热| 边亲边吃奶的免费视频| 一级毛片 在线播放| 最近手机中文字幕大全| 丝袜喷水一区| 亚洲欧美日韩东京热| 免费高清在线观看视频在线观看| 亚洲欧洲日产国产| 久久久久精品性色| 免费黄频网站在线观看国产| 久久国产精品大桥未久av | 最近最新中文字幕大全电影3| 国产午夜精品一二区理论片| h日本视频在线播放| 蜜臀久久99精品久久宅男| 久久ye,这里只有精品| 建设人人有责人人尽责人人享有的 | 国产精品不卡视频一区二区| 久久久午夜欧美精品| 91在线精品国自产拍蜜月| 国产91av在线免费观看| 国产男女超爽视频在线观看| 色吧在线观看| 国产片特级美女逼逼视频| 国产在线男女| 一区二区三区精品91| 少妇猛男粗大的猛烈进出视频| 精品人妻一区二区三区麻豆| 午夜福利在线观看免费完整高清在| 国产精品免费大片| 久久精品夜色国产| 精品少妇黑人巨大在线播放| 最新中文字幕久久久久| av国产久精品久网站免费入址| 精品一区二区三卡| av国产免费在线观看| 大香蕉97超碰在线| 精品亚洲乱码少妇综合久久| 国产精品国产三级专区第一集| 亚洲综合精品二区| 夫妻性生交免费视频一级片| 一二三四中文在线观看免费高清| 精品少妇黑人巨大在线播放| 欧美性感艳星| 18+在线观看网站| 18禁裸乳无遮挡动漫免费视频| 特大巨黑吊av在线直播| freevideosex欧美| 久久久久国产精品人妻一区二区| 成人毛片a级毛片在线播放| 97热精品久久久久久| 免费看光身美女| 亚洲国产欧美人成| 亚洲av综合色区一区| 日韩,欧美,国产一区二区三区| 汤姆久久久久久久影院中文字幕| 黑人高潮一二区| 天天躁日日操中文字幕| 国产午夜精品一二区理论片| 欧美日韩精品成人综合77777| 大香蕉97超碰在线| 亚洲在久久综合| 亚洲av.av天堂| 麻豆精品久久久久久蜜桃| 免费看光身美女| 一本色道久久久久久精品综合| 中文字幕精品免费在线观看视频 | 亚洲精品成人av观看孕妇| 久久鲁丝午夜福利片| 一区在线观看完整版| 91精品伊人久久大香线蕉| 国产日韩欧美亚洲二区| 午夜福利高清视频| av又黄又爽大尺度在线免费看| 美女高潮的动态| 亚洲中文av在线| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲国产精品成人久久小说| 国产淫语在线视频| 国内揄拍国产精品人妻在线| 国内精品宾馆在线| 亚洲av中文字字幕乱码综合| 在线观看一区二区三区| 久久午夜福利片| a级一级毛片免费在线观看| 汤姆久久久久久久影院中文字幕| 亚洲高清免费不卡视频| 黄色欧美视频在线观看| 自拍偷自拍亚洲精品老妇| 啦啦啦中文免费视频观看日本| 国产大屁股一区二区在线视频| 国产亚洲5aaaaa淫片| 插阴视频在线观看视频| 国产精品不卡视频一区二区| 国产高潮美女av| 18禁在线播放成人免费| av在线app专区| 寂寞人妻少妇视频99o| 舔av片在线| 久久久久国产精品人妻一区二区| 九草在线视频观看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美清纯卡通| 久久毛片免费看一区二区三区| 日本色播在线视频| 热99国产精品久久久久久7| 亚洲美女搞黄在线观看| 日韩人妻高清精品专区| 啦啦啦视频在线资源免费观看| 丝瓜视频免费看黄片| 精品熟女少妇av免费看| 久久韩国三级中文字幕| 国产av一区二区精品久久 | 91久久精品国产一区二区三区| 建设人人有责人人尽责人人享有的 | 内射极品少妇av片p| 国产精品秋霞免费鲁丝片| 久久97久久精品| 在线亚洲精品国产二区图片欧美 | videossex国产| 建设人人有责人人尽责人人享有的 | 国产av一区二区精品久久 | 激情 狠狠 欧美| 亚洲av成人精品一区久久| 寂寞人妻少妇视频99o| 久久综合国产亚洲精品| 晚上一个人看的免费电影| 制服丝袜香蕉在线| 国产精品爽爽va在线观看网站| 国产亚洲一区二区精品| 久久女婷五月综合色啪小说| a级毛片免费高清观看在线播放| 国产精品伦人一区二区| av在线蜜桃| 日本av手机在线免费观看| 亚洲人成网站高清观看| 国产精品国产三级专区第一集| 一本久久精品| 欧美人与善性xxx| 亚洲av成人精品一二三区| 国内精品宾馆在线| 九草在线视频观看| 日本猛色少妇xxxxx猛交久久| 婷婷色综合www| 午夜免费男女啪啪视频观看| 亚洲欧美精品专区久久| 97热精品久久久久久| 久久av网站| 亚洲一级一片aⅴ在线观看| 国产乱来视频区| 亚洲av二区三区四区| 日韩av免费高清视频| 天美传媒精品一区二区| 2018国产大陆天天弄谢| 午夜激情福利司机影院| 精品酒店卫生间| 一级a做视频免费观看| 97在线人人人人妻| 自拍偷自拍亚洲精品老妇| 2021少妇久久久久久久久久久| 欧美激情极品国产一区二区三区 | 色网站视频免费| 成人漫画全彩无遮挡| 免费少妇av软件| 久久久久久久久大av| 久久精品国产亚洲网站| 亚洲欧美日韩东京热| 国产亚洲91精品色在线| 91精品国产九色| 国产爱豆传媒在线观看| 最后的刺客免费高清国语| 亚洲欧美日韩东京热| 国产精品一区二区在线不卡| 一级毛片aaaaaa免费看小| 最新中文字幕久久久久| 男人爽女人下面视频在线观看| 十八禁网站网址无遮挡 | 亚洲av福利一区| 少妇人妻精品综合一区二区| 成人高潮视频无遮挡免费网站| 一级a做视频免费观看| 亚洲国产精品国产精品| 精品国产乱码久久久久久小说| 久久精品久久久久久久性| 国产乱人视频| 久久久久精品久久久久真实原创| 久久久国产一区二区| 久久99热6这里只有精品| 交换朋友夫妻互换小说| 久久午夜福利片| 国产亚洲av片在线观看秒播厂| 日韩不卡一区二区三区视频在线| 亚洲av不卡在线观看| 欧美精品一区二区免费开放| a级毛色黄片| 国产乱人视频| 99久国产av精品国产电影| 久久精品国产亚洲av涩爱| 一级毛片电影观看| 久久久久久久久久久丰满| 日本av免费视频播放| 国产乱来视频区| 亚洲av中文av极速乱| 国产在视频线精品| 黄色配什么色好看| 校园人妻丝袜中文字幕| 久久av网站| 91久久精品国产一区二区三区| 乱码一卡2卡4卡精品| 这个男人来自地球电影免费观看 | av视频免费观看在线观看| 日本黄大片高清| 嫩草影院入口| 成人亚洲精品一区在线观看 | 在线亚洲精品国产二区图片欧美 | 直男gayav资源| 日韩不卡一区二区三区视频在线| 免费在线观看成人毛片| 一级毛片黄色毛片免费观看视频| 久久99热这里只频精品6学生| 尤物成人国产欧美一区二区三区| 国产精品偷伦视频观看了| 国产精品熟女久久久久浪| 国产成人精品福利久久| 成人特级av手机在线观看| 亚洲av成人精品一区久久| 婷婷色麻豆天堂久久| 国产精品久久久久成人av| 网址你懂的国产日韩在线| 综合色丁香网| 91狼人影院| 亚洲欧洲日产国产| 在线观看人妻少妇| 国产 一区精品| 久久久久性生活片| 国产乱人视频| 欧美xxxx性猛交bbbb| 在线看a的网站| 六月丁香七月| 少妇的逼水好多| 亚洲精品乱久久久久久| 国产久久久一区二区三区| 久久99精品国语久久久| 大香蕉97超碰在线| 欧美区成人在线视频| 国精品久久久久久国模美| 亚洲性久久影院| 在线观看国产h片| 久久婷婷青草| 新久久久久国产一级毛片| 亚洲欧美日韩无卡精品| 国产精品一二三区在线看| 久久精品人妻少妇| 下体分泌物呈黄色| 男女免费视频国产| 日本欧美国产在线视频| 你懂的网址亚洲精品在线观看| videossex国产| 夫妻性生交免费视频一级片| 六月丁香七月| 免费av不卡在线播放| 大片免费播放器 马上看| 亚洲色图av天堂| 久久99热这里只频精品6学生| 欧美人与善性xxx| 国产精品久久久久久av不卡| 成人毛片a级毛片在线播放| 一级av片app| 99热网站在线观看| av视频免费观看在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产有黄有色有爽视频| 国产精品国产三级国产专区5o| 国产精品秋霞免费鲁丝片| 精品一品国产午夜福利视频| 久久97久久精品| 男女免费视频国产| 美女视频免费永久观看网站| 黑人高潮一二区| 亚洲国产精品999| av国产久精品久网站免费入址| 中文字幕av成人在线电影| 欧美极品一区二区三区四区| 2018国产大陆天天弄谢| 久久久久国产网址| 国产淫片久久久久久久久| 深夜a级毛片| 国产精品麻豆人妻色哟哟久久| 国产精品精品国产色婷婷| 大陆偷拍与自拍| 国产精品99久久99久久久不卡 | 成人亚洲欧美一区二区av| 少妇人妻精品综合一区二区| 国产精品秋霞免费鲁丝片| 久久这里有精品视频免费| 亚洲aⅴ乱码一区二区在线播放| 纵有疾风起免费观看全集完整版| 亚洲国产成人一精品久久久| 日韩av免费高清视频| 免费大片黄手机在线观看| 2021少妇久久久久久久久久久| 亚洲欧美精品专区久久| 国产一区有黄有色的免费视频| 久久综合国产亚洲精品| 天堂中文最新版在线下载| 啦啦啦啦在线视频资源| 亚州av有码| 99国产精品免费福利视频| 亚洲四区av| 国产一区二区三区av在线| 久久久久久久亚洲中文字幕| 日本av手机在线免费观看| 寂寞人妻少妇视频99o| 嘟嘟电影网在线观看| 国产亚洲91精品色在线| 岛国毛片在线播放| 日韩 亚洲 欧美在线| 国产伦精品一区二区三区视频9| 午夜免费观看性视频| av卡一久久| freevideosex欧美| 最近最新中文字幕大全电影3| 一级毛片久久久久久久久女| 高清午夜精品一区二区三区| 超碰97精品在线观看| 日韩中文字幕视频在线看片 | 春色校园在线视频观看| 欧美日韩视频精品一区| 男女无遮挡免费网站观看| 亚洲欧美日韩卡通动漫| 日韩人妻高清精品专区| 国产在线一区二区三区精| 国产精品一区二区在线不卡| 欧美日韩一区二区视频在线观看视频在线| 亚洲av综合色区一区| 啦啦啦在线观看免费高清www| 国产成人freesex在线| 久久久久网色| 国内揄拍国产精品人妻在线| 国产精品国产av在线观看| 99久久综合免费| 亚洲av国产av综合av卡| 中文在线观看免费www的网站| 在线观看美女被高潮喷水网站| 亚洲第一区二区三区不卡| 欧美成人精品欧美一级黄| 国产亚洲午夜精品一区二区久久| 欧美最新免费一区二区三区| 建设人人有责人人尽责人人享有的 | 一区二区三区四区激情视频| 精品人妻偷拍中文字幕| 高清在线视频一区二区三区| 少妇裸体淫交视频免费看高清| 成人毛片60女人毛片免费| 高清毛片免费看| 亚洲国产精品一区三区| 欧美xxⅹ黑人| 国产爽快片一区二区三区| 日韩强制内射视频| 免费观看在线日韩| 女的被弄到高潮叫床怎么办| 精品久久国产蜜桃| 亚洲成人av在线免费| 成年免费大片在线观看| 国产男人的电影天堂91| 久久久久久久久久久免费av| a级毛色黄片| 午夜激情久久久久久久| h日本视频在线播放| 色视频www国产| 一边亲一边摸免费视频| 久久国产精品大桥未久av | 麻豆成人av视频| 精品一区二区三区视频在线| 亚洲国产av新网站| 欧美人与善性xxx| 乱系列少妇在线播放| 我的女老师完整版在线观看| 人妻制服诱惑在线中文字幕| 国产美女午夜福利| 日日撸夜夜添|