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

    銑削過程顫振穩(wěn)定性分析的研究進展

    2016-02-09 15:20:15盧曉紅王鳳晨王鑫鑫司立坤
    振動與沖擊 2016年1期
    關(guān)鍵詞:穩(wěn)定性分析

    盧曉紅, 王鳳晨, 王 華, 王鑫鑫, 司立坤

    (大連理工大學(xué) 精密與特種加工教育部重點實驗室, 大連 116024)

    ?

    銑削過程顫振穩(wěn)定性分析的研究進展

    盧曉紅, 王鳳晨, 王華, 王鑫鑫, 司立坤

    (大連理工大學(xué) 精密與特種加工教育部重點實驗室, 大連116024)

    摘要:綜述銑削過程顫振穩(wěn)定性分析的研究概況和進展。顫振建模和穩(wěn)定性分析是該方法兩個關(guān)鍵環(huán)節(jié)。依據(jù)顫振形成的物理條件,將其分為摩擦型顫振、振型耦合型顫振和再生型顫振。從切削過程的非線性和切削系統(tǒng)的非線性兩方面,重點介紹再生型顫振的非線性建模的研究成果。穩(wěn)定性分析方法根據(jù)對顫振模型的求解方法,分為頻域法、離散法及數(shù)值法,概括了各個方法的特點、效果及適用工況。最后介紹了近來興起的微細銑削研究領(lǐng)域中顫振穩(wěn)定性分析的研究成果。由于其尺度效應(yīng),微細銑削加工具有獨特的加工機理和特點,顫振建模中需考慮的因素與傳統(tǒng)銑削多有不同,但穩(wěn)定性分析方法仍大多沿用傳統(tǒng)銑削中的方法。

    關(guān)鍵詞:銑削;顫振;穩(wěn)定性分析;微細銑削

    銑削加工過程中刀具與工件間的相對振動,是降低產(chǎn)品的精度及表面質(zhì)量,影響生產(chǎn)效率的主要原因。顫振是在切削過程中,在無周期性外部激振力持續(xù)作用下由于加工系統(tǒng)本身特性所引起的一種切削振動,一般表現(xiàn)為金屬切削過程中刀具與工件之間強烈的相對振動。顫振的發(fā)生降低了切削用量和工件表面質(zhì)量,產(chǎn)生大量的噪聲,甚至?xí)?dǎo)致刀具的提前報廢。

    通過銑削過程顫振穩(wěn)定性分析,可獲取無顫振切削參數(shù)組合,以達到優(yōu)化工藝參數(shù),抑制振動的目的。國內(nèi)外學(xué)者在銑削過程顫振穩(wěn)定性分析領(lǐng)域的方法大多為從銑削過程動力學(xué)模型出發(fā)求得由切削參數(shù)構(gòu)成的加工穩(wěn)定域邊界。本文不僅介紹了銑削顫振穩(wěn)定性分析方法中幾個關(guān)鍵環(huán)節(jié)的國內(nèi)外研究成果,而且對近年來興起的微細銑削研究領(lǐng)域中顫振穩(wěn)定性分析的研究進展進行了詳細說明。

    1銑削過程顫振建模

    1.1銑削過程動力學(xué)模型

    顫振按照形成的物理條件可大體分為三種形式:摩擦型顫振、振型耦合型顫振和再生型顫振。

    (1) 摩擦型顫振

    摩擦型顫振是由切削速度方向上刀具與工件之間相互摩擦引起的,而產(chǎn)生這種相互摩擦的原因很多。Tlusty等[1]認為該原因是切削力相對于切削速度與刀具前角的動態(tài)變化的相位滯后;Sisson等[2]則通過研究刀具所受的切削力,得出切削速度增大時切削力下降的特性是主要原因的結(jié)論。國內(nèi)劉習(xí)軍等[3]建立了由非線性動態(tài)切削力相耦合的刀具彈性子系統(tǒng)與工件彈性子系統(tǒng)的多自由度非線性系統(tǒng)的顫振模型,以數(shù)值模擬計算方法,發(fā)現(xiàn)了系統(tǒng)存在內(nèi)共振現(xiàn)象;之后,又利用平均法求出了系統(tǒng)的1∶2內(nèi)共振的近似解析解,找出了切削速度影響切削顫振的動態(tài)切削條件[4]。

    (2) 振型耦合型顫振

    振型耦合型顫振是由于振動系統(tǒng)在兩個方向上剛度相差不大,導(dǎo)致兩個固有振型相耦合引起的。一般采用實驗?zāi)B(tài)分析方法即可求得系統(tǒng)的穩(wěn)定性條件。Tlusty等[5]提出振型耦合型顫振,研究對象一般為兩自由度線性振動系統(tǒng)。Gasparetto[6]根據(jù)系統(tǒng)特征方程,利用特征值分析方法,研究振型耦合型顫振的穩(wěn)定條件,作出了詳細的闡述。Iturrospe等[7]利用狀態(tài)空間分析將系統(tǒng)穩(wěn)定狀態(tài)與系統(tǒng)自然頻率定量地聯(lián)系到一起,并提出了一種以控制系統(tǒng)自然頻率變化來探測切削過程穩(wěn)定狀態(tài)的方法。

    振型耦合和再生效應(yīng)在多自由度的切削系統(tǒng)中往往同時存在。Huang等[8]建立了兩種效應(yīng)同時作用下的立銑顫振模型。研究表明:過程矩陣的特征值決定銑削穩(wěn)定性,而過程矩陣的特征向量影響結(jié)構(gòu)振動軌跡;通過分析振動軌跡,發(fā)現(xiàn)顫振會在振動能量周期性累積的方向上旋轉(zhuǎn)。Moradi等[9]利用多尺度法,研究了在再生效應(yīng)影響下XY方向的耦合模態(tài)之間能量傳遞。結(jié)果表明,當(dāng)能量傳遞給一個方向上的模態(tài)時,另一方向的振動會衰減,這樣可以通過將振動能量傳遞給表面質(zhì)量要求較低的方向上來抑制表面質(zhì)量要求較高的方向上的振動。林潔瓊等[10]綜合考慮再生效應(yīng)和振型耦合效應(yīng),建立了隨機擾動激勵下兩自由度切削加工系統(tǒng)的顫振時滯動力學(xué)模型,分析了剛度主軸方位對切削穩(wěn)定性的影響。

    (3) 再生型顫振

    再生型顫振是由于上次切削所形成的振紋與本次切削的振動位移之間的相位差導(dǎo)致刀具切削厚度的不同而引起的顫振。由于再生效應(yīng)的存在,切削過程中某一時刻顫振的狀態(tài)依賴于切削系統(tǒng)在過去某一時刻的狀態(tài),一般用時滯微分方程組來描述。Tlusty和星鐵太郎為代表所作的研究在線性模型范圍內(nèi)發(fā)展了一整套較為完善的理論與方法[11]。線性再生型顫振理論,一般將加工過程用反饋系統(tǒng)來表示,利用控制理論中的方法求解系統(tǒng)方程的穩(wěn)定性條件與穩(wěn)定性圖。

    切削系統(tǒng)本質(zhì)上是復(fù)雜的非線性系統(tǒng)。對于再生型顫振最初的研究目的是更容易分析其穩(wěn)定性,線性再生顫振理論對于切削系統(tǒng)做了大量的簡化假設(shè),不能描述顫振發(fā)生、發(fā)展和自穩(wěn)定的全過程,也無法解釋實際切削顫振中的很多現(xiàn)象。Hooke等考慮到切削過程和機床加工系統(tǒng)的非線性后,首先建立了非線性再生顫振模型。Stépán等[12], Balachandran等[13-15],師漢民[16-17]及龍新華等[18-19]也在非線性再生顫振理論方面做了很多有益的研究。

    1.2顫振模型的主要影響因素分析

    非線性動態(tài)銑削力被普遍認為是影響顫振模型的主要因素。切削過程中,切削力對于切削厚度存在非線性依賴關(guān)系。在再生顫振理論模型中,動態(tài)銑削力通常被表達為動態(tài)切削厚度的函數(shù),一般為動態(tài)切削厚度的三次冪函數(shù)和四分之三冪函數(shù)[12]。Hanna等[20]基于機床結(jié)構(gòu)的非線性和切削力的非線性,建立了常系數(shù)差分微分方程,其中以動態(tài)切削厚度的二次冪和三次冪函數(shù)描述切削力的非線性。大振幅的斷續(xù)切削也是導(dǎo)致顫振的重要因素。當(dāng)振幅較大時,刀具運動軌跡會有一部分越出工件范圍之外,與工件脫離接觸;此時,動態(tài)切削厚度為零,動態(tài)切削力也為零,這與銑削中刀齒周期性間斷切削工件的現(xiàn)象不同,是造成切削加工中非連續(xù)、非線性的主要因素之一。在考慮大振幅的斷續(xù)切削的同時,切削過程中的多重再生效應(yīng)也不可忽略。多重再生效應(yīng),簡而言之,就是在決定某次切削的動態(tài)切削厚度時,不僅要考慮上一次切削時刀具-工件系統(tǒng)間的動態(tài)位移,還需要考慮上一次乃至上上次切削時刀具留在工件上的波紋。由于大振幅的斷續(xù)切削和多重再生效應(yīng)的存在,動態(tài)切削厚度需表達為非光滑函數(shù)。Balachandran等[13]考慮到大振幅的斷續(xù)切削和多重再生效應(yīng),建立了適用于部分徑向切深銑削及全徑向切深銑削等多種工況的銑削顫振模型,利用龐加萊截面圖和數(shù)值方法預(yù)測了其穩(wěn)定及失穩(wěn)狀態(tài)下的運動,預(yù)測結(jié)果與實驗觀察結(jié)果很接近。Banihasan等[21]建立了詳盡的二自由度顫振模型,除兩個不同方向的狀態(tài)依賴時滯外,首次提出同時考慮切削區(qū)中多個不同滯后刀具位置的多重再生效應(yīng)的準確幾何模型,指出分段多重再生效應(yīng)是銑削過程中非線性振動的根本特征,最大Lyapunov指數(shù)和分岔的數(shù)值計算證明了高速切削時存在混沌振動。

    切削過程中,進給率對于時滯以至系統(tǒng)穩(wěn)定性的影響也為國內(nèi)外學(xué)者所廣泛研究。Zhao等[14-15]研究了進給率對時滯和大振幅斷續(xù)切削效應(yīng)的影響,建立了顫振模型。該模型為含周期系數(shù)矩陣的非線性非齊次非自治方程組,在進給方向和其他方向上有著不同的常數(shù)時滯項。研究結(jié)果表明,線性顫振模型對于大徑向銑削過程可以做出相當(dāng)準確的穩(wěn)定性預(yù)測,但對于小徑向銑削過程不能準確的預(yù)測其穩(wěn)定性。龍新華等[18-19]除了考慮進給率對于時滯的作用之外,還考慮了其對于切入角和切出角的影響,以含周期系數(shù)矩陣的非自治時滯微分方程組描述四自由度銑削系統(tǒng)顫振,得到了較精準的非線性顫振穩(wěn)定葉瓣圖。該顫振模型中,由于進給率的影響,時滯項不再是常數(shù)項而是時變項,并且在X和Y方向上不同。

    上述研究說明了銑削過程中的非線性,如切削厚度的時滯特性、大振幅的斷續(xù)切削、多重再生效應(yīng)等對于穩(wěn)定性預(yù)測及切削顫振的影響。此外,切削系統(tǒng)結(jié)構(gòu)的非線性,如刀具結(jié)構(gòu),機床結(jié)構(gòu)的動態(tài)特性,也都是影響切削顫振模型的重要因素。

    在刀具結(jié)構(gòu)對于切削系統(tǒng)穩(wěn)定性的影響方面,關(guān)于變齒距或不等齒距銑刀的研究較多,主要因為在切削過程中,使用變齒距的銑刀比使用均勻齒距銑刀穩(wěn)定性更好。銑刀的齒距是切削刃上的點到下一個切削刃上同一個點的距離。此類結(jié)構(gòu)的銑刀可以分為兩種,變螺旋角銑刀和變齒間角銑刀。Altintas等[22]針對變齒間角銑刀修改了銑削顫振解析模型,在頻域內(nèi)繪制了以機床主軸轉(zhuǎn)速和軸向切深為橫、縱坐標(biāo)的穩(wěn)定性葉瓣圖。研究結(jié)果表明,相比于等齒距銑刀,變齒距銑刀在切削過程中穩(wěn)定性更好,更不易發(fā)生顫振。因此,優(yōu)化設(shè)計這類變齒距銑刀也是國內(nèi)外研究者抑制切削過程中顫振的方法之一。Budak[23-24]從理論和試驗兩方面分析和證明了變齒間角銑刀在提高材料去除率和降低表面粗糙度等方面有明顯的作用,尤其在低主軸速度時,可以得到很高的無顫振軸向切深。Turner等[25]將Budak的分析方法應(yīng)用到了變螺旋角銑刀,并與等螺旋角銑刀及變齒間角銑刀的穩(wěn)定性預(yù)測進行了對比,結(jié)果表明變螺旋角銑刀同樣可以很大改善切削穩(wěn)定性。Sellmeier等[26]認為,因等齒間角銑刀的動態(tài)銑削過程需用單時滯項的非自治微分方程組來描述,而不等齒間角銑刀的動態(tài)銑削過程應(yīng)用多時滯量的微分方程組來描述;研究結(jié)果表明,在單自由度系統(tǒng)的穩(wěn)定性葉瓣圖中,有 “穩(wěn)態(tài)島”出現(xiàn)。Otto等[27]針對變齒距銑刀的分布式時滯,建立解析模型,在頻域內(nèi)拓展了零階求解法,通過求解多階頻率矩陣的特征值繪制了更為完整的穩(wěn)定性圖。Dombovari等[28]考慮銑刀螺旋角的連續(xù)變化,建立了分布式時滯微分方程組。他認為連續(xù)變化的螺旋角會引起切削過程中時滯的連續(xù)變化,因此軸向力的分布可以用關(guān)于滯后時間的加權(quán)分布函數(shù)來表示。該模型探究了變螺旋角銑刀的一些復(fù)雜穩(wěn)定特性,有助于優(yōu)化設(shè)計更高材料去除率的新型銑刀。Eksioglu等[29]提出了銑削動態(tài)系統(tǒng)的離散時間模型,任意形狀的銑刀沿刀具軸線被劃分為微分單元,適用于變齒間角和變螺旋角刀具。

    鋸齒銑刀由于其良好的抑振效果得國內(nèi)外學(xué)者的深入研究。Dombovari等[30]在時頻域內(nèi)建立了鋸齒銑刀切削穩(wěn)定性模型。由于沿切削刃每一點的時滯都與鋸齒結(jié)構(gòu)和進給率有關(guān),利用含周期系統(tǒng)矩陣的多時滯微分方程組來描述切削動態(tài)過程。研究表明:在切削相同體積的材料時,鋸齒銑刀需要的驅(qū)動扭矩更??;即使在切削鈦之類的難加工材料時,鋸齒銑刀可采用無顫振切削深度也更大,加工效率更高。

    機床結(jié)構(gòu),主要是刀具與主軸結(jié)合面,以及主軸本身的非線性動態(tài)特性,也是影響切削顫振模型的重要因素。Catania等[31]建立了銑床結(jié)構(gòu)的實驗?zāi)B(tài)模型,將機床結(jié)構(gòu)對顫振模型的影響考慮在內(nèi)。李勤良等[32]考慮結(jié)合面遲滯非線性的基礎(chǔ)上,建立含有非線性激振力的機床單自由顫振模型,利用多尺度法求解系統(tǒng)振幅和相位的分岔方程,對顫振系統(tǒng)的穩(wěn)定性進行研究,并在此基礎(chǔ)上分析轉(zhuǎn)速和滯回參數(shù)等對系統(tǒng)穩(wěn)定性產(chǎn)生的影響。Gao等[33]發(fā)現(xiàn),主軸-夾具-刀具系統(tǒng)的結(jié)構(gòu)特性對于切削過程穩(wěn)定性的影響在過往研究中往往以刀尖頻響函數(shù)的形式體現(xiàn),而主軸結(jié)構(gòu)本身往往受到忽視。在這方面,Tian等[34-36]考慮了旋轉(zhuǎn)主軸的陀螺效應(yīng)對于切削穩(wěn)定特性的影響,而Erturk等[37]探討了主軸-夾具-刀具裝配結(jié)構(gòu)特性的作用。Gao等[33]利用數(shù)值模擬方法研究了主軸中球軸承對于切削過程穩(wěn)定性的影響,他發(fā)現(xiàn)主軸中球軸承與滾道的非線性非光滑赫茲接觸剛度是復(fù)雜振動響應(yīng)的根源之一,尤其在大軸承間隙的情況下。

    主軸的動態(tài)特性對于切削顫振影響的研究主要集中在主軸速度變化上。Radulescu等[38-39]的研究指出變主軸速度切削可以降低切削振動振幅而且其切削系統(tǒng)的穩(wěn)定性相比于常主軸速度切削,對于切削系統(tǒng)模態(tài)參數(shù)變化的魯棒性更好。Takemura等[40]研究了主軸速度不同變化曲線下的切削穩(wěn)定性特點,變化曲線包括三角波,方波和正弦波。Sastry等[41]提出了一個針對變主軸速度下切削穩(wěn)定性的分析方法。該方法通過含時變周期系數(shù)及時變時滯的微分差分方程組描述動態(tài)切削系統(tǒng)的閉環(huán)模型,利用Floquet理論和傅里葉分析將穩(wěn)定問題簡化為確定無窮階特征方程根的位置?;赟astry等[41]的研究,Zatarain等[42]提出了一個頻域內(nèi)適用于任何主軸速度變化曲線的通用分析方法。龍新華等[43]利用半離散法,對變主軸速度下順銑和逆銑的穩(wěn)定性分別進行了研究,發(fā)現(xiàn)穩(wěn)定區(qū)內(nèi)的軸向切深深度與常主軸速度切削相比有明顯改善,另外變主軸速度下逆銑的穩(wěn)定性改善要比順銑更好。Seguy等[44]以三角波及正弦波的主軸速度變化曲線,研究主軸變化的幅度和頻率對于穩(wěn)定性的影響,結(jié)果表明低頻率、大幅度的主軸變速更利于改善穩(wěn)定性。

    2顫振穩(wěn)定性分析

    銑削穩(wěn)定性分析的研究內(nèi)容一般為預(yù)測由臨界加工參數(shù)構(gòu)成的銑削過程穩(wěn)定性邊界,從而劃分出穩(wěn)定切削區(qū)和不穩(wěn)定切削區(qū),繼而可以從穩(wěn)定性葉瓣圖中選擇適當(dāng)?shù)那邢鲄?shù),達到避免顫振,提高加工效率和加工質(zhì)量的目的。

    目前求解銑削加工穩(wěn)定性域的方法可以大體分為顫振模型求解和實驗法兩類。

    2.1顫振模型求解法

    顫振模型求解根據(jù)其對顫振模型的求解方法又可分為頻域法,離散法和數(shù)值法。

    (1) 頻域法

    頻域法一般為將顫振模型的時滯微分方程組利用傅里葉變換轉(zhuǎn)到頻域表示,基于控制理論,解析計算銑削穩(wěn)定邊界。Altintas等[45-46]首先提出了零階求解法。該方法求解基于再生效應(yīng)的兩自由度銑削顫振模型,將隨時間變化的定向動態(tài)銑削力系數(shù)矩陣,通過傅里葉變換到頻域,以傅里葉級數(shù)展開的平均量(即零次諧波分量)對其近似,并忽略了刀具-工件接觸區(qū)的交叉?zhèn)鬟f函數(shù),求解閉環(huán)系統(tǒng)特征方程的特征值,得到無顫振條件下的軸向切深。該方法所需要的計算條件較少,并且計算效率很高,但由于過于簡化,不能預(yù)報低徑向銑削工況時出現(xiàn)的倍周期分岔,精度較差。為解決這個問題,Merdol等[47]提出了多頻率法,該方法考慮了定向因子的高次諧波,在計算過程中需要迭代搜索顫振頻率,需求解多個特征值。Bachrathy等[48]將多頻率法擴展到適用于所有刀具結(jié)構(gòu)的穩(wěn)定性預(yù)測,包括可引入分布時滯的復(fù)雜刀具結(jié)構(gòu),把擴展多頻率法和多維二分法結(jié)合大幅提高了計算效率,并證明了在所測頻響函數(shù)質(zhì)量較差的情況下,擴展多頻率法依然可以得到可靠的穩(wěn)定預(yù)測結(jié)果。

    (2) 離散法

    再生顫振模型一般為時滯微分方程組,穩(wěn)定性由系統(tǒng)單值算子的特征值所決定,但單值算子由無窮維矩陣表示,造成求解閉合形式的穩(wěn)定預(yù)測結(jié)果十分困難。離散法通常是利用有限維轉(zhuǎn)換矩陣去近似無窮維的單值算子,從而降低求解難度,減少計算時間。離散法主要包括半離散法、全離散法和時間有限元法三種。

    Insperger等[49]為求解含周期時間系數(shù)矩陣的時滯微分方程組,從計算流體力學(xué)和有限元分析領(lǐng)域借鑒了半離散法;之后Insperger等[50]改進該算法,應(yīng)用于單自由度及兩自由度銑削系統(tǒng)的穩(wěn)定性預(yù)測。該方法用每個時間區(qū)段兩個相鄰滯后狀態(tài)值的加權(quán)和近似該時間區(qū)段的時滯項,同時對周期系數(shù)項做分片零階平均處理,從而構(gòu)造出單周期上逼近原微分動力系統(tǒng)的離散動力系統(tǒng),獲得單周期的Floquet狀態(tài)轉(zhuǎn)移矩陣,最后以狀態(tài)轉(zhuǎn)移矩陣的譜半徑大小來判斷是否穩(wěn)定。Insperger等[51]基于提高收斂性的目的,提出了一階半離散法:如周期時間系數(shù)項利用分段常值函數(shù)近似,對時滯項只需做不高于一階的逼近;為了取得更高階的收斂,周期時間系數(shù)需要用更高階的方法來逼近,如Magnus級數(shù)。Altintas等[52]比較了頻域法和半離散法,認為零階求解法對周期方向系數(shù)取平均,無需迭代,計算速度快,適用于大徑向切深比及每齒切削周期與結(jié)構(gòu)自然頻率接近的工況;多頻率法所需的頻響函數(shù)質(zhì)量要求不高,無需識別模態(tài)參數(shù),適用于低徑向切深比及刀齒數(shù)較少的間斷切削工況;半離散法的精度及計算速度取決于使用的離散時間間隔和節(jié)點數(shù)量,可考慮變齒距刀具和變主軸速度等復(fù)雜工況,不過半離散法和多頻率法相對于零階求解法來說,分析更準確,但需要迭代搜索穩(wěn)定性邊界,需要更多的計算時間。該方向上,后來的研究重點主要集中在利用不同的方法提高計算速度和收斂性上,如李中偉等[53]提出基于Magnus-Gaussian截斷的零階半離散穩(wěn)定分析法,Niu等[54]提出基于Runge-Kutta的分析方法等。

    Ding等[55]提出了一種基于數(shù)值積分的全離散法。該方法在對時間周期等距離散后,積分得到Floquet轉(zhuǎn)移矩陣,亦基于Floquet理論判穩(wěn)。全離散法與半離散法相似。二者的不同在于半離散基于微分方程而全離散法基于積分方程;半離散僅對時滯項及周期時間系數(shù)矩陣離散近似,而全離散法同步離散時滯項、當(dāng)前狀態(tài)和周期系數(shù)項。相同的離散時間間隔下,全離散法因矩陣指數(shù)函數(shù)只依賴轉(zhuǎn)速,因此計算時間更短,但收斂速度比一階半離散法低[56]?,F(xiàn)有的一些改進半離散法計算效率的方法也適用于全離散法[56]。全離散法同樣可用于多時滯項的穩(wěn)定性分級[57],并可進一步提高收斂性[58]。Li等[59]提出了另一類全離散法,該方法離散了所有時間相關(guān)項,利用數(shù)值迭代方法得到迭代公式,從而得到Floquet轉(zhuǎn)移矩陣。

    時間有限元法一般通過匹配近似切削中刀具運動和切入切出時刀具自由響應(yīng)的位移和速度,獲得切削時間區(qū)段上的刀具運動方程,再利用加權(quán)余量法得到單個刀齒切削周期的Floquet轉(zhuǎn)移矩陣,由Floquet理論判斷銑削穩(wěn)定性。Bayly等[60]在對單自由度銑削系統(tǒng)穩(wěn)定性預(yù)測中,針對間斷切削中刀具切削與非切削狀態(tài)描述方程不同等問題,使用了時間有限元法,隨后又將該方法擴展到兩自由度銑削系統(tǒng)。Garg等[61]將該方法擴展到含周期系數(shù)參數(shù)激勵的時滯系統(tǒng)穩(wěn)定性分析,并比較了使用單類型單元、提高多項式級數(shù)逼近精確解和使用多類型單元及三次多項式,增加單元數(shù)量逼近精確解這兩種方法的優(yōu)劣。Mann等[62]之后又將方法擴展到可應(yīng)用于所有可用空間狀態(tài)模型表示的系統(tǒng)。Ding等[63]基于時間有限元法,以積分方程技術(shù)求解銑削動態(tài)系統(tǒng)響應(yīng),提出數(shù)值積分法,并將其發(fā)展為具有指數(shù)收斂階的譜方法[64]和針對多時滯工況的變步長法[65]。

    綜上所述,可知離散法核心在于將有限維矩陣去近似無限維單值算子,方法之間的區(qū)別在于近似的具體手段,類似的還有Ulsoy等[66-67]提出的基于Lambert W函數(shù)的方法,Butcher等[68-69]提出的Chebyshev多項式法及Chebyshev配點法及Ding等[70]的微分求積法等。

    (3) 數(shù)值法

    數(shù)值法通常是直接求解時滯微分方程獲取銑削過程動態(tài)響應(yīng),然后由響應(yīng)構(gòu)造的穩(wěn)定判據(jù)或響應(yīng)振動幅值是否發(fā)散,來判斷銑削過程的穩(wěn)定性。Tlusty等[5, 71]提出了基于差分格式的時域系統(tǒng)動態(tài)響應(yīng)計算方法。Smith等[72]建立了切削力峰-峰值時域仿真模型,以切削力峰-峰值的變化作為顫振判據(jù)。Altintas等[73]針對低徑向切深的銑削建立了一種改進的時域模型,以仿真預(yù)測的動態(tài)切削厚度與靜態(tài)切削厚度比值,作為無量綱顫振判別系數(shù)。Li等[74]采用仿真預(yù)測的最大動態(tài)切削力與最大靜態(tài)切削力之比作為顫振判定標(biāo)準。Li等[75]提出了利用Runge-Kutta法的銑削穩(wěn)定性數(shù)值求解法。時域仿真方法能夠考慮非線性銑削力,大振幅斷續(xù)切削等非線性因素,其應(yīng)用相當(dāng)廣泛,但計算效率低。

    2.2實驗法

    實驗法是通過切削試驗,采集切削過程中的振動、切削力或噪音等信號,以適當(dāng)?shù)姆椒ǚ治鲂盘?、判別顫振,最后得到銑削穩(wěn)定性葉瓣圖。如Quintana等[76-77]通過采集加工過程中的噪聲信號,通過快速傅里葉變換(FFT)進行顫振識別,記錄不同主軸轉(zhuǎn)速發(fā)生顫振時的軸向切深值,繪制穩(wěn)定性曲線。遲玉倫等[78]進行了基于聲發(fā)射信號的顫振穩(wěn)定域葉瓣圖確定實驗方法研究。實驗法繪制穩(wěn)定性葉瓣圖在實際生產(chǎn)加工中具有一定的實用價值,但由于實驗誤差及多變的生產(chǎn)加工條件,其得到的無顫振臨界加工參數(shù)的可靠性及通用性都較差。實驗法多用來與顫振模型求解得到的穩(wěn)定性預(yù)測結(jié)果對比,以證明顫振模型及穩(wěn)定性分析方法的正確性。

    3微細銑削過程顫振穩(wěn)定性分析

    微銑削加工的零部件尺寸一般介于100 μm到10 mm之間,特征尺寸一般為10 μm~1 mm,加工過程中的振動尤其是刀具和工件之間的相對振動對零部件質(zhì)量和精度影響明顯。微銑削加工所用微銑刀直徑通常在1 mm以下,剛度低、易磨損、易折斷,振動會加速刀具磨損,嚴重時甚至造成刀具折斷,導(dǎo)致加工成本過高,制約微銑削加工技術(shù)應(yīng)用。因此抑制微銑削加工過程中的振動對于提高加工精度、延長微銑刀壽命、降低加工成本具有非常重要的意義。

    微細銑削顫振建模目前多是沿用傳統(tǒng)銑削的顫振建模方法和步驟,主要考慮再生效應(yīng)的影響,根據(jù)切削動力學(xué)系統(tǒng)的控制方程,分別求取控制方程中刀具的模態(tài)參數(shù)及切削力系數(shù),之后基于模型,以適當(dāng)方法分析其穩(wěn)定性。但微細銑削并不是宏觀銑削加工在尺度上的簡單縮小,其所具有的介觀尺度加工的一些特征(如尺度效應(yīng)、最小切削厚度及單齒切削現(xiàn)象等)以及工件材料微觀結(jié)構(gòu)的影響,使得微細銑削具有獨特的加工機理和特點。

    傳統(tǒng)銑削獲取刀具模態(tài)參數(shù)(一般以刀尖頻響函數(shù)的形式)的方法主要為實驗法,即基于力錘沖擊的模態(tài)實驗方法,而微細銑削的刀具結(jié)構(gòu)微小脆弱,力錘的錘擊很難直接施加在刀尖上,而且由于現(xiàn)有硬件限制,力錘模態(tài)實驗的激勵頻率通常僅能達到10 kHz左右,遠小于100 kHz以上的微銑刀自然頻率。Mascaydell等[79]根據(jù)動柔度耦合法,將微銑刀分為兩部分,以錘擊法求取刀柄部分的頻響函數(shù),以有限元分析獲取刀尖部分的頻響函數(shù),將兩部分耦合得到完整微銑刀的刀尖頻響函數(shù)。Filiz等[80]考慮了切削刃的幾何形狀,應(yīng)用Timoshenko梁理論求解刀具的頻響函數(shù)。Tajalli等[81-83]根據(jù)應(yīng)變梯度Timoshenko梁理論、擴展哈密頓原理建立了考慮陀螺效應(yīng)的旋轉(zhuǎn)刀具動態(tài)模型,利用精確動態(tài)剛度方法以期得到更準確的刀具模態(tài)參數(shù)。Uhlmann等[84]以Shi提出的壓電驅(qū)動器激振法得到刀具的前三階模態(tài)參數(shù),再基于建立了刀具的旋轉(zhuǎn)Euler-Bernouli梁模型。Mustapha等[85]建立可預(yù)估微銑刀橫向響應(yīng)的混合分析模型,該模型中不同類的結(jié)構(gòu)單元表示刀具的剛度、阻尼、幾何結(jié)構(gòu)等不同的特性。

    微細銑削中,與工件材料接觸并進行切削的是切削刃刃口圓弧半徑。由于切削厚度和切削刃刃口圓弧半徑一般在一個數(shù)量級,存在一個產(chǎn)生連續(xù)切削的臨界切削厚度值,即最小切削厚度。當(dāng)實際切削厚度大于最小切削厚度時,剪切效應(yīng)為主導(dǎo),工件材料發(fā)生剪切滑移變形;當(dāng)實際切削厚度小于最小切削厚度時,耕犁效應(yīng)為主導(dǎo),工件材料發(fā)生彈性變形。因此微細銑削中不同的切削厚度下切削力系數(shù)矩陣可能不同,而且在耕犁效應(yīng)為主導(dǎo)的情況下,工件材料的彈性變形會引起過程中很大的阻尼。過程阻尼力,尤其在低主軸轉(zhuǎn)速下,也是影響銑削過程穩(wěn)定性的重要因素。Rahnama等[86]考慮了切削刃圓弧半徑及銑削過程中的過程阻尼并利用刀具和工件間的等效接觸體積來定義過程阻尼參數(shù),但僅考慮了只有剪切效應(yīng)存在的工況,忽略了耕犁效應(yīng)及其帶來的非線性。Afazov等[87]考慮了剪切效應(yīng)、耕犁效應(yīng)、切削刃口圓弧半徑及刀具徑向跳動,有限元預(yù)測了正交切削下切削力,描述了切削力,切削速度和切削厚度的非線性關(guān)系,之后Afazov等[88]進一步研究了工件初始溫度的影響,比較了滑移-黏摩擦模型和黏滯摩擦模型在正交切削有限元建模中的準確性。Jin等[89]基于刀具結(jié)構(gòu)和材料本構(gòu)特性,建立了滑移線區(qū)域模型,得到了切削力系數(shù),過程阻尼力系數(shù)則是有限元分析獲得。以上的研究都沒有將剪切效應(yīng)和耕犁效應(yīng)下的切削力區(qū)分考慮,僅在力模型中考慮了耕犁效應(yīng)帶來的過程阻尼力,Song等[90]則根據(jù)以最小切削厚度為臨界點,剪切效應(yīng)及耕犁效應(yīng)下不同的切削力模型,利用時域仿真方法獲得了穩(wěn)定性預(yù)測圖。

    微細銑削的穩(wěn)定性分析方法仍采用傳統(tǒng)銑削使用的方法,如零階求解法,半離散法及時域法。其中時域法方面,Afazov等[87]利用Runge-Kutta方法求解,以x,y方向振動位移的統(tǒng)計方差作為顫振判別準則;Baschin等[91-92]在其研究中獲取銑刀運動軌跡,根據(jù)龐加萊截面和分岔理論來判別顫振;Park等[93]認為系統(tǒng)動態(tài)特性及切削參數(shù)可能在切削過程中發(fā)生變化,利用基于棱邊定理和零點排除法的魯棒性顫振穩(wěn)定性理論,繪制了魯棒性穩(wěn)定性葉瓣圖。

    4結(jié)論

    本文對銑削過程顫振穩(wěn)定性分析中動力學(xué)建模和穩(wěn)定性預(yù)測的研究進行了歸納和綜述,并對微細銑削領(lǐng)域的顫振穩(wěn)定性分析的研究進行了論述。

    顫振動力學(xué)模型方面,以再生型、振型耦合型、摩擦型為代表的顫振理論已相當(dāng)成熟。更精確的“機床-刀具-工件-夾具”系統(tǒng)動力學(xué)模型,針對難加工材料的低速切削過程動力學(xué)模型,多種顫振類型綜合的動力學(xué)模型,考慮工件材料微觀結(jié)構(gòu)的微細銑削動力學(xué)模型等是該領(lǐng)域有待進一步研究的課題。

    穩(wěn)定性分析方面,頻域法、離散法、數(shù)值法這三類主流方法發(fā)展已相當(dāng)成熟,但它們都有各自的局限性。頻域法計算效率高,由于其簡化,導(dǎo)致精度不高,忽略了很多復(fù)雜的非線性效應(yīng),不同工況的通用性較差;離散法通用性及精度較好,但計算效率不高;數(shù)值法可以考慮多種非線性效應(yīng),但計算效率最低,穩(wěn)定性判據(jù)通用性差。穩(wěn)定性分析方法仍待進一步的改進,以提高其通用性、精度及計算效率。微細銑削領(lǐng)域大多沿用銑削的穩(wěn)定性分析方法,針對微細銑削本身的加工機理及特點的方法仍有待研究。

    參 考 文 獻

    [1] Tlusty J, Spacek L. Self-excited vibration in machine tools[M]. Prague, 1954.

    [2] Sisson T R, Kegg R L. An explanation of low-speed chatter effects[J]. Journal of Engineering for Industry,1969,91(4): 951.

    [3] 劉習(xí)軍,陳予恕. 機床速度型切削顫振的非線性研究[J]. 振動與沖擊,1999,18(2): 7-11.

    LIU Xi-jun, CHEN Yu-shu.Nonlinear analysis of speed type cutting chatter of machine tools[J]. Journal of Vibration and Shock,1999,18(2): 7-11.

    [4] 劉習(xí)軍,王立剛,賈啟芬. 一種由干摩擦引起的車床切削顫振[J]. 工程力學(xué),2005(1): 107-112.

    LIU Xi-jun, WANG Li-gang, JIA Qi-fen.Cutting chatter of machine tools induced by dry friction[J]. Engineering Mechanics,2005(1): 107-112.

    [5] Tlusty J, Ismail F. Basic non-linearity in machining chatter[J]. CIRP Annals-Manufacturing Technology,1981, 30(1): 299-304.

    [6] Gasparetto A. Eigenvalue analysis of modecoupling chatter for machine-tool stabilization[J]. Journal of Vibration and Control, 2001(7): 181-197.

    [7] Iturrospe A, Atxa V, Abete J M. State-space analysis of mode-coupling in orthogonal metal cuttingunder wave regeneration[J]. International Journal of Adaptive Control and Signal Processing,2007(47): 1583-1590.

    [8] Huang C Y, Wang J. Mode coupling behavior in end milling[M]. New York:Amer Soc Mechanical Engineers, 2010, 875-884.

    [9] Moradi H, Movahhedy M R, Vossoughi G, et al. Investigation of the internal resonance and regenerative chatter dynamics in nonlinear milling process[J]. Proceedings of the Asme International Design Engineering Technical Conferences and Computers and Information in Engineering Conference 2011, Vol 1, Pts a and B: 23Rd Biennial Conference On Mechanical Vibration and Noise,2012: 141-150.

    [10] 林潔瓊,周曉勤,孔繁森,等. 剛度主軸方位對模態(tài)耦合再生切削系統(tǒng)動態(tài)響應(yīng)諧參數(shù)的影響[J]. 振動與沖擊,2009,28(5): 63-68.

    LIN Jie-qiong, ZHOU Xiao-qin, KONG Fan-sen, et al.Influence of orientations of principal stiffness axes on harmonic parameters of dynamic response in mode-coupled regenerative machining system[J]. Journal of Vibration and Shock, 2009,28(5):63-68.

    [11] 吳雅. 機床切削系統(tǒng)的顫振及其控制[M]. 北京: 科學(xué)出版社, 1993.

    [12] Stépán G, Kalmár-Nagy T. Nonlinear Regenerative Machine Tool Vibrations[C]//1997.

    [13] Balachandran B. Nonlinear dynamics of milling processes[J]. Philosophical Transactions of the Royal Society of London Series A-Mathematical Physical and Engineering Sciences,2001, 359(1781): 793-819.

    [14] Zhao M X, Balachandran B. Dynamics and stability of milling process[J]. International Journal of Solids and Structures,2001, 38(10): 2233-2248.

    [15] Balachandran B, Zhao M X. A mechanics based model for study of dynamics of milling operations[J]. Meccanica,2000, 35(2): 89-109.

    [16] 師漢民. 關(guān)于機床自激振動的一個非線性理論模型(第一部分)[J]. 應(yīng)用力學(xué)學(xué)報,1984(1): 1-14.

    SHI Han-min. Nonlinear chatter theory of machine tools (part 1: the stablity of machine tool chatter amplitude) [J]. Journal of Applied Mechanics,1984(1): 1-14.

    [17] 師漢民. 關(guān)于機床自激振動的一個非線性理論模型(第二部份)[J]. 應(yīng)用力學(xué)學(xué)報,1984(2): 75-88.

    SHI Han-min. Nonlinear chatter theory of machine tools (part 2: finite amplitude machine tool instability) [J]. Journal of Applied Mechanics, 1984(2): 75-88.

    [18] Long X H, Balachandran B. Stability analysis for milling process[J]. Nonlinear Dynamics,2007, 49(3): 349-359.

    [19] Long X H, Balachandran B, Mann B P. Dynamics of milling processes with variable time delays[J]. Nonlinear Dynamics,2007, 47(1/2/3): 49-63.

    [20] Hanna N H, Tobias S A. A theory of nonlinear regenerative chatter[J]. ASME Journal of Engineering for Industry, 1974, 96: 247-255.

    [21] Banihasan M, Bakhtiari-Nejad F. Chaotic vibrations in high-speed milling[J]. Nonlinear Dynamics, 2011, 66(4): 557-574.

    [23] Budak E. An analytical design method for milling cutters with nonconstant pitch to increase stability, part 1: theory[J]. Journal of Manufacturing Science and Engineering-Transactions of the Asme,2003, 125(1): 29-34.

    [24] Budak E. An analytical design method for milling cutters with nonconstant pitch to increase stability, part 2: application[J]. Journal of Manufacturing Science and Engineering-Transactions of the Asme,2003, 125(1): 35-38.

    [25] Turner S, Merdol D, Altintas Y, et al. Modelling of the stability of variable helix end mills[J]. International Journal of Machine Tools and Manufacture,2007, 47(9): 1410-1416.

    [26] Sellmeier V, Denkena B. Stable islands in the stability chart of milling processes due to unequal tooth pitch[J]. International Journal of Machine Tools & Manufacture, 2011, 51(2): 152-164.

    [27] Otto A, Radons G. Frequency domain stability analysis of milling processes with variable helix tools[C]//2012.

    [28] Dombovari Z, Stepan G. The effect of helix angle variation on milling stability[J]. Journal of Manufacturing Science and Engineering-Transactions of the Asme, 2012, 134(0510155).[29] Eksioglu C, Kilic Z M, Altintas Y. Discrete-time prediction of chatter stability, cutting forces, and surface location errors in flexible milling systems[J]. Journal of Manufacturing Science and Engineering, 2012, 134(6): 61006.

    [30] Dombovari Z, Altintas Y, Stepan G. The effect of serration on mechanics and stability of milling cutters[J]. International Journal of Machine Tools and Manufacture, 2010, 50(6): 511-520.

    [31] Catania G, Mancinelli N. A coupled theoretical-experimental dynamical model for chatter prediction in milling processes[J]. Proceedings of the Asme International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, 2010,1: 175-184.

    [32] 李勤良,汪博,趙斌,等. 考慮非線性遲滯力的機床顫振系統(tǒng)穩(wěn)定性研究[J]. 機械工程學(xué)報, 2013(11): 43-49.

    LI Qin-liang, WANG Bo, ZHAO Bin, et al.Research on the chatter stability of machine system taking the nonlinear hysteretic force into consideration[J]. Journal of Mechanical Engineering, 2013(11): 43-49.

    [33] Gao S H, Meng G, Long X H. Study of milling stability with hertz contact stiffness of ball bearings[J]. Archive of Applied Mechanics, 2011, 81(8): 1141-1151.

    [34] Movahhedy M R, Mosaddegh P. Prediction of chatter in high speed milling including gyroscopic effects[J]. International Journal of Machine Tools and Manufacture,2006, 46(9): 996-1001.

    [35] Xiong G L, Yi J M, Zeng C, et al. Study of the gyroscopic effect of the spindle on the stability characteristics of the milling system[J]. Journal of Materials Processing Technology,2003, 138(1): 379-384.

    [36] Tian J, Hutton S G. Chatter instability in milling systems with flexible rotating spindles—a new theoretical approach[J]. Journal of Manufacturing Science and Engineering,2001,123(1): 1-9.

    [37] Ertürk A, Budak E, ?zgüven H N. Selection of design and operational parameters in spindle-holder-tool assemblies for maximum chatter stability by using a new analytical model[J]. International Journal of Machine Tools and Manufacture,2007, 47(9): 1401-1409.

    [38] Radulescu R, Kapoor S G, Devor R E. An investigation of variable spindle speed face milling for tool-work structures with complex dynamics, part 1: simulation results[J]. Journal of Manufacturing Science and Engineering, 1997, 119(3): 266-272.

    [39] Radulescu R, Kapoor S G, Devor R E. An investigation of variable spindle speed face milling for tool-work structures with complex dynamics, part 2: physical explanation[J]. Journal of Manufacturing Science and Engineering, 1997, 119(3): 273-280.

    [40] Takemura T, Kitamura T, Hoshi T, et al. Active suppression of chatter by programmed variation of spindle speed[J]. Annals of the CIRP, 1974, 23(1): 121-122.

    [41] Sastry S, Kapoor S G, Devor R E. Floquet theory based approach for stability analysis of the variable speed face-milling process[J]. Journal of Manufacturing Science and Engineering-Transactions of the Asme,2002, 124(1): 10-17.

    [42] Zatarain M, Bediaga I, Munoa J, et al. Stability of milling processes with continuous spindle speed variation: analysis in the frequency and time domains, and experimental correlation[J]. Cirp Annals-Manufacturing Technology,2008, 57(1): 379-384.

    [43] Xinhua L, Balachandran B. Stability of up-milling and down-milling operations with variable spindle speed[J]. Journal of Vibration and Control, 2010, 16(7/8): 1151-1168.

    [44] Seguy S, Insperger T, Arnaud L, et al. Suppression of period doubling chatter in high-speed milling by spindle speed variation[J]. Machining Science and Technology, 2011,15(PII 9380210952): 153-171.

    [45] Budak E, Altintas Y. Analytical prediction of chatter stability in milling. I. general formulation[J]. Transactions of the ASME. Journal of Dynamic Systems, Measurement and Control, 1998, 120(1): 22-30.

    [47] Merdol S D, Altintas Y. Multi frequency solution of chatter stability for low immersion milling[J]. Transactions of the ASME. Journal of Manufacturing Science and Engineering,2004, 126(3): 459-466.

    [48] Bachrathy D, Stepan G. Improved prediction of stability lobes with extended multi frequency solution[J]. Cirp Annals-Manufacturing Technology, 2013, 62(1): 411-414.

    [49] Insperger T, Stepan G. Semi-discretization method for delayed systems[J]. International Journal for Numerical Methods in Engineering,2002, 55(5): 503-518.

    [50] Insperger T, Stepan G. Updated semi-discretization method for periodic delay-differential equations with discrete delay[J]. International Journal for Numerical Methods in Engineering,2004, 61(1): 117-141.

    [51] Insperger T, Stepan G, Turi J. On the higher-order semi-discretizations for periodic delayed systems[J]. Journal of Sound and Vibration,2008, 313(1/2): 334-341.

    [52] Altintas Y, Stepan G, Merdol D, et al. Chatter stability of milling in frequency and discrete time domain[J]. CIRP Journal of Manufacturing Science and Technology,2008,1(1): 35-44.

    [53] 李中偉,龍新華,孟光. 基于Magnus-Gaussian截斷的銑削系統(tǒng)穩(wěn)定性的半離散分析法[J]. 振動與沖擊,2009,28(5): 69-73.

    LI Zhong-wei, LONG Xin-hua, MENG Guang.Stability analysis of milling process by semi discretization method based on Magnus-Guassian truncation[J]. Journal of Vibration and Shock,2009,28(5): 69-73.

    [54] Niu J, Ding Y, Zhu L, et al. Runge-kutta methods for a semi-analytical prediction of milling stability[J]. Nonlinear Dynamics, 2014, 76(1): 289-304.

    [55] Ding Y, Zhu L, Zhang X, et al. A full-discretization method for prediction of milling stability[J]. International Journal of Machine Tools and Manufacture, 2010, 50(5): 502-509.

    [56] Insperger T. Full-discretization and semi-discretization for milling stability prediction: some comments[J]. International Journal of Machine Tools & Manufacture, 2010, 50(7): 658-662.

    [57] Zhang X, Xiong C, Ding Y. Improved full-discretization method for milling chatter stability prediction with multiple delays[M]. Intelligent Robotics and Applications, Springer, 2010:541-552.

    [58] Ding Y, Zhu L, Zhang X, et al. Second-order full-discretization method for milling stability prediction[J]. International Journal of Machine Tools and Manufacture, 2010, 50(10): 926-932.

    [59] Li M, Zhang G, Huang Y. Complete discretization scheme for milling stability prediction[J]. Nonlinear Dynamics, 2013, 71(1/2): 187-199.

    [60] Bayly P V, Halley J E, Mann B P, et al. Stability of interrupted cutting by temporal finite element analysis[J]. Journal of Manufacturing Science and EngineeringJournal of Manufacturing Science and Engineering,2003,125(2):220.

    [61] Garg N K, Mann B P, Kim N H, et al. Stability of a time-delayed system with parametric excitation[J]. Journal of Dynamic Systems Measurement and Control-Transactions of the Asme,2007, 129(2): 125-135.

    [62] Mann B P, Patel B R. Stability of delay equations written as state space models[J]. Journal of Vibration and Control, 2010, 16(7-8SI): 1067-1085.

    [63] Ding Y, Zhang X, Ding H, et al. Numerical integration method for prediction of milling stability[J].Journal of Manufacturing Science and Engineering,2011,133(3):31005.

    [64] Ding Y, Zhu L, Zhang X, et al. Milling stability analysis using the spectral method[J]. Science China Technological Sciences, 2011, 54(12): 3130-3136.

    [65] Zhang X, Xiong C, Ding Y, et al. Variable-step integration method for milling chatter stability prediction with multiple delays[J]. Science China Technological Sciences, 2011,54(12): 3137-3154.

    [66] Yi S, Nelson P, Ulsoy A. Delay differential equations via the matrix lambert W function and bifurcation analysis: application to machine tool chatter[J]. Mathematical Biosciences and Engineering,2007, 4(2): 355.

    [67] Maghami Asl F, Ulsoy A G. Analysis of a system of linear delay differential equations[J]. Journal of Dynamic Systems, Measurement and Control, Transactions of the ASME,2003, 125(2): 215-223.

    [68] Butcher E A, Bobrenkov O A, Bueler E, et al. Analysis of milling stability by the chebyshev collocation method: algorithm and optimal stable immersion levels[J]. Journal of Computational and Nonlinear Dynamics,2009, 4(3): 1-12.

    [69] Butcher E A, Ma H, Bueler E, et al. Stability of linear time-periodic delay-differential equations via chebyshev polynomials[J]. International Journal for Numerical Methods in Engineering,2004, 59(7): 895-922.

    [70] Ding Y, Zhu L, Zhang X, et al. Stability analysis of milling via the differential quadrature method[J]. Journal of Manufacturing Science and Engineering,2013,135(4):44502.

    [71] Tlusty J, Ismail F. Special aspects of chatter in milling[J]. Journal of Vibration, Acoustics Stress and Reliability in Design,1983, 105(1): 24-32.

    [72] Smith S, Tlusty J. Efficient simulation programs for chatter in milling[J]. CIRP Annals-Manufacturing Technology,1993, 42(1): 463-466.

    [73] Campomanes M L, Altintas Y. An improved time domain simulation for dynamic milling at small radial immersions[J]. Journal of Manufacturing Science and Engineering,2003,125(3): 416-422.

    [74] Li H Z, Li X P, Chen X Q. A novel chatter stability criterion for the modelling and simulation of the dynamic milling process in the time domain[J]. The International Journal of Advanced Manufacturing Technology,2003, 22(9/10): 619-625.

    [75] Zhongqun L, Qiang L. Solution and analysis of chatter stability for end milling in the time-domain[J]. Chinese Journal of Aeronautics,2008, 21(2): 169-178.

    [76] Quintana G, Ciurana J, Ferrer I, et al. Sound mapping for identification of stability lobe diagrams in milling processes[J]. International Journal of Machine Tools and Manufacture,2009, 49(3): 203-211.

    [77] Quintana G, Ciurana J, Teixidor D. A new experimental methodology for identification of stability lobes diagram in milling operations[J]. International Journal of Machine Tools and Manufacture,2008, 48(15): 1637-1645.

    [78] 遲玉倫, 李郝林. 銑削顫振穩(wěn)定域葉瓣圖確定方法研究[J]. 振動與沖擊, 2014,33(4):90-93.

    CHI Yu-lun,LI Hao-lin.Determination of chatter stability field lobe diagrams for a milling processing[J].Journal of Vibration and Shock,2014,33(4):90-93.

    [79] Mascardelli B A, Park S S, Freiheit T. Substructure coupling of microend mills to aid in the suppression of chatter[J]. Journal of Manufacturing Science and Engineering,2008,130(1):11010.

    [80] Filiz S, Ozdoganlar O B. Microendmill dynamics including the actual fluted geometry and setup errors—part I: model development and numerical solution[J]. Journal of Manufacturing Science and Engineering,2008, 130(3): 31119.

    [81] Tajalli S A, Movahhedy M R, Akbari J. Chatter instability analysis of spinning micro-end mill with process damping effect via semi-discretization approach[J]. Acta Mechanica, 2014, 225(3): 715-734.

    [82] Tajalli S A, Movahhedy M R, Akbari J. Size dependent vibrations of micro-end mill incorporating strain gradient elasticity theory[J]. Journal of Sound and Vibration, 2013, 332(15): 3922-3944.

    [83] Tajalli S A, Movahhedy M R, Akbari J. Investigation of the effects of process damping on chatter instability in micro end milling[J]. Procedia CIRP, 2012: 156-161.

    [84] Uhlmann E, Mahr F. A time domain simulation approach for micro milling processes[J]. Procedia CIRP, 2012, 4(0): 22-28.

    [85] Mustapha K B, Zhong Z W. A hybrid analytical model for the transverse vibration response of a micro-end mill[J]. Mechanical Systems and Signal Processing, 2013,34(1/2): 321-339.

    [86] Rahnama R, Sajjadi M, Park S S. Chatter suppression in micro end milling with process damping[J]. Journal of Materials Processing Technology,2009, 209(17): 5766-5776.

    [87] Afazov S M, Ratchev S M, Segal J, et al. Chatter modelling in micro-milling by considering process nonlinearities[J]. International Journal of Machine Tools and Manufacture, 2012: 28-38.

    [88] Afazov S M, Zdebski D, Ratchev S M, et al. Effects of micro-milling conditions on the cutting forces and process stability[J]. Journal of Materials Processing Technology, 2013, 213(5): 671-684.

    [89] Jin X, Altintas Y. Chatter stability model of micro-milling with process damping[J]. Journal of Manufacturing Science and Engineering, 2013, 135(3): 31011.

    [90] Song Q, Liu Z, Shi Z. Chatter stability for micromilling processes with flat end mill[J]. International Journal of Advanced Manufacturing Technology, 2014, 71(5/6/7/8): 1159-1174.

    [91] Biermann D, Baschin A. Influence of cutting edge geometry and cutting edge radius on the stability of micromilling processes[J]. Production Engineering,2009, 3(4/5): 375-380.

    [92] Baschin A, Kahnis P, Biermann D. Dynamic analysis of the micromilling process-influence of tool vibrations on the quality of microstructures[J]. Materialwissenschaft Und Werkstofftechnik,2008, 39(9): 616-621.

    [93] Park S S, Rahnama R. Robust chatter stability in micro-milling operations[J]. CIRP Annals-Manufacturing Technology, 2010, 59(1): 391-394.

    基金項目:國家自然科學(xué)基金(51305061 );中央高?;究蒲袠I(yè)務(wù)費專項資金資助(DUT13LAB13)

    收稿日期:2014-06-24修改稿收到日期:2014-09-30

    中圖分類號:TH113;TG54

    文獻標(biāo)志碼:A

    DOI:10.13465/j.cnki.jvs.2016.01.014

    Review about chatter stability analysis in milling process

    LU Xiao-hong, WANG Feng-chen, WANG Hua, WANG Xin-xin, SI Li-kun

    (Key Laboratory for Precision and Non-traditional Machining Technology of Ministy of Education, University of Technology, Dalian 116024, China)

    Abstract:Literatures about chatter stability analysis in milling process were reviewed. Chatter modeling and stability analysis were two key procedures of the studying method. The mechanisms leading to chatter were classified as dry fiction effect, mode coupling effect and regenerative effect. Nonlinear regenerative chatter modeling was mainly introduced according to nonlinearity in cutting processes and spindle-holder-tool systems. Stability analysis was divided into three categories, i.e., frequency domain method, discretization method and numerical simulation method based on the dynamic model solving methods. The features and applicable conditions of each method were described. The study achievements of chatter stability analysis in the micro-milling area were introduced. Due to its size effect, micro-milling had its unique cutting mechanism and characteristics. The factors to be considered in micro-milling chatter modeling were different from those in traditional milling, but the methods of stability analysis were mostly similar.

    Key words:milling; chatter; stability analysis; micro-milling

    第一作者 盧曉紅 女,博士,副教授,1978年生

    猜你喜歡
    穩(wěn)定性分析
    電廠灰渣庫穩(wěn)定性分析
    元壩某井場進場道路1號滑坡穩(wěn)定性分析及防治措施
    科技視界(2016年27期)2017-03-14 11:38:59
    高聳鋼結(jié)構(gòu)施工關(guān)鍵控制技術(shù)分析
    價值工程(2016年36期)2017-01-11 19:49:21
    框架預(yù)應(yīng)力錨桿邊坡支護結(jié)構(gòu)及其應(yīng)用分析
    低聚季銨鹽對聚驅(qū)采出水包油乳狀液破乳機理
    淺談邊坡穩(wěn)定性地質(zhì)問題的解決措施
    一種基于區(qū)間分割的時滯系統(tǒng)的鎮(zhèn)定控制
    民用飛機起落架擺振分析與仿真
    科技視界(2015年27期)2015-10-08 12:13:09
    有關(guān)軟弱結(jié)構(gòu)面的巖質(zhì)邊坡穩(wěn)定性分析
    居業(yè)(2015年6期)2015-09-05 19:58:14
    基于物理模擬的某邊坡危巖體穩(wěn)定性分析
    科技視界(2015年12期)2015-04-28 09:25:52
    av线在线观看网站| www.999成人在线观看| 久久国产精品影院| 又黄又粗又硬又大视频| av一本久久久久| 一本色道久久久久久精品综合| 香蕉国产在线看| av有码第一页| 久久狼人影院| 人人妻,人人澡人人爽秒播| 少妇粗大呻吟视频| 国产亚洲av高清不卡| 夜夜夜夜夜久久久久| 久久性视频一级片| 久久免费观看电影| 国产av精品麻豆| 一个人免费看片子| 日韩制服丝袜自拍偷拍| 777久久人妻少妇嫩草av网站| 在线 av 中文字幕| 少妇粗大呻吟视频| 亚洲欧美色中文字幕在线| 久久天堂一区二区三区四区| 亚洲av电影在线进入| 日韩视频一区二区在线观看| 国产在线免费精品| 国产真人三级小视频在线观看| 国产精品影院久久| 久久久久久人人人人人| 丰满饥渴人妻一区二区三| 曰老女人黄片| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲成人免费av在线播放| 国产一级毛片在线| 韩国高清视频一区二区三区| 国产精品免费视频内射| 热99国产精品久久久久久7| 精品久久久精品久久久| 一本—道久久a久久精品蜜桃钙片| 久久精品成人免费网站| 成年人免费黄色播放视频| 宅男免费午夜| 亚洲第一青青草原| 亚洲av成人不卡在线观看播放网 | 啪啪无遮挡十八禁网站| av超薄肉色丝袜交足视频| 日韩三级视频一区二区三区| 欧美黑人欧美精品刺激| 一二三四社区在线视频社区8| 久久亚洲精品不卡| 国产主播在线观看一区二区| 丝瓜视频免费看黄片| 黑人巨大精品欧美一区二区蜜桃| 老司机靠b影院| 丰满饥渴人妻一区二区三| 黄色视频在线播放观看不卡| 国产精品熟女久久久久浪| 色老头精品视频在线观看| 热99国产精品久久久久久7| 男人爽女人下面视频在线观看| 亚洲国产精品一区三区| 欧美国产精品一级二级三级| 欧美中文综合在线视频| 亚洲精品国产一区二区精华液| 一进一出抽搐动态| 18在线观看网站| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美精品自产自拍| 日本a在线网址| 啪啪无遮挡十八禁网站| 亚洲黑人精品在线| 99国产综合亚洲精品| 一个人免费在线观看的高清视频 | 亚洲av国产av综合av卡| 超色免费av| 性色av乱码一区二区三区2| 国内毛片毛片毛片毛片毛片| 天天添夜夜摸| 狠狠婷婷综合久久久久久88av| 91精品三级在线观看| 十八禁网站免费在线| 国产精品二区激情视频| 欧美在线黄色| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩 亚洲 欧美在线| 久久久久久免费高清国产稀缺| 丝袜在线中文字幕| av天堂在线播放| 宅男免费午夜| 国产精品国产av在线观看| 午夜福利一区二区在线看| 久久精品aⅴ一区二区三区四区| 国产精品国产av在线观看| 另类亚洲欧美激情| 久久av网站| 亚洲精品美女久久久久99蜜臀| 亚洲成av片中文字幕在线观看| 19禁男女啪啪无遮挡网站| 99精品欧美一区二区三区四区| www.999成人在线观看| av有码第一页| 免费高清在线观看日韩| 在线精品无人区一区二区三| www.999成人在线观看| 人人妻人人澡人人爽人人夜夜| 99精国产麻豆久久婷婷| 成年美女黄网站色视频大全免费| 精品国产一区二区三区久久久樱花| 90打野战视频偷拍视频| 精品国产乱码久久久久久男人| 欧美xxⅹ黑人| 中国美女看黄片| 久久免费观看电影| 青草久久国产| 亚洲人成电影免费在线| 视频区图区小说| 国产男人的电影天堂91| 欧美精品av麻豆av| 窝窝影院91人妻| 国产99久久九九免费精品| 大片免费播放器 马上看| 久久精品亚洲熟妇少妇任你| 久热爱精品视频在线9| 每晚都被弄得嗷嗷叫到高潮| 下体分泌物呈黄色| 一级a爱视频在线免费观看| 亚洲精品国产区一区二| 自线自在国产av| 免费黄频网站在线观看国产| 日韩 欧美 亚洲 中文字幕| 亚洲视频免费观看视频| 91精品三级在线观看| 女性生殖器流出的白浆| 亚洲一码二码三码区别大吗| 嫁个100分男人电影在线观看| 欧美日韩福利视频一区二区| 国产成人免费观看mmmm| 亚洲精品国产色婷婷电影| 欧美精品av麻豆av| 啦啦啦视频在线资源免费观看| netflix在线观看网站| 97人妻天天添夜夜摸| 一二三四在线观看免费中文在| 少妇粗大呻吟视频| 国产精品久久久av美女十八| 欧美变态另类bdsm刘玥| 狠狠狠狠99中文字幕| 亚洲人成电影观看| 免费在线观看黄色视频的| 91老司机精品| 国产精品一二三区在线看| av在线app专区| 青春草视频在线免费观看| 狠狠婷婷综合久久久久久88av| 国产成人精品久久二区二区91| 超色免费av| 亚洲精品乱久久久久久| 国产精品欧美亚洲77777| av国产精品久久久久影院| 日韩免费高清中文字幕av| 亚洲专区字幕在线| 两性午夜刺激爽爽歪歪视频在线观看 | 777久久人妻少妇嫩草av网站| 亚洲欧美日韩另类电影网站| 国产成人欧美| 精品一区二区三卡| 啦啦啦 在线观看视频| 亚洲人成77777在线视频| 国产成人精品无人区| 久9热在线精品视频| 亚洲国产av新网站| 亚洲国产中文字幕在线视频| av在线播放精品| 亚洲国产精品成人久久小说| 97精品久久久久久久久久精品| 另类精品久久| 91成人精品电影| 人人妻人人爽人人添夜夜欢视频| 国产av精品麻豆| 亚洲精品久久久久久婷婷小说| 久久性视频一级片| 日本vs欧美在线观看视频| 欧美在线黄色| h视频一区二区三区| 91成年电影在线观看| 免费观看av网站的网址| 最新的欧美精品一区二区| 亚洲av男天堂| 国产成人免费无遮挡视频| 亚洲五月色婷婷综合| 久久国产精品男人的天堂亚洲| 99国产精品一区二区蜜桃av | 在线永久观看黄色视频| 伊人亚洲综合成人网| 日韩有码中文字幕| 热99re8久久精品国产| 夫妻午夜视频| 亚洲国产欧美网| 麻豆乱淫一区二区| 精品人妻1区二区| 亚洲第一青青草原| 啪啪无遮挡十八禁网站| 欧美人与性动交α欧美精品济南到| 一本—道久久a久久精品蜜桃钙片| 人人妻人人添人人爽欧美一区卜| 在线 av 中文字幕| 中文字幕制服av| 久久精品国产亚洲av香蕉五月 | 高潮久久久久久久久久久不卡| 色婷婷久久久亚洲欧美| 热re99久久精品国产66热6| 亚洲av成人不卡在线观看播放网 | 久久久久国产一级毛片高清牌| 精品少妇久久久久久888优播| 色播在线永久视频| 麻豆乱淫一区二区| 亚洲第一av免费看| www.av在线官网国产| 免费看十八禁软件| 国产男人的电影天堂91| 午夜免费鲁丝| 97人妻天天添夜夜摸| 日本五十路高清| 欧美黑人欧美精品刺激| 亚洲av欧美aⅴ国产| av福利片在线| 久久这里只有精品19| 久久影院123| 精品高清国产在线一区| 久久国产精品影院| 日本欧美视频一区| 女性生殖器流出的白浆| 欧美成狂野欧美在线观看| 中文字幕色久视频| 精品少妇内射三级| 国产精品国产三级国产专区5o| 亚洲欧美色中文字幕在线| 大香蕉久久网| 两性夫妻黄色片| 日韩大片免费观看网站| 久久久久久久精品精品| av免费在线观看网站| 亚洲伊人色综图| 亚洲久久久国产精品| 交换朋友夫妻互换小说| 视频在线观看一区二区三区| 久久久久国产一级毛片高清牌| 在线观看免费高清a一片| 少妇 在线观看| 男女午夜视频在线观看| 19禁男女啪啪无遮挡网站| 亚洲五月婷婷丁香| 法律面前人人平等表现在哪些方面 | 黄色视频在线播放观看不卡| 老熟妇仑乱视频hdxx| 另类亚洲欧美激情| 狠狠狠狠99中文字幕| 日韩制服骚丝袜av| 国产99久久九九免费精品| 精品少妇一区二区三区视频日本电影| 亚洲国产毛片av蜜桃av| 性高湖久久久久久久久免费观看| 91国产中文字幕| 久久人人爽av亚洲精品天堂| 精品久久久精品久久久| 自线自在国产av| 成年人免费黄色播放视频| 中文欧美无线码| 亚洲国产中文字幕在线视频| 亚洲av片天天在线观看| 飞空精品影院首页| 亚洲成人国产一区在线观看| 欧美xxⅹ黑人| 国产麻豆69| 99九九在线精品视频| 欧美+亚洲+日韩+国产| 啦啦啦 在线观看视频| 免费人妻精品一区二区三区视频| 国产亚洲午夜精品一区二区久久| 超碰成人久久| 韩国精品一区二区三区| 在线天堂中文资源库| 制服诱惑二区| 久久人人爽av亚洲精品天堂| 亚洲成人免费电影在线观看| 久久精品aⅴ一区二区三区四区| 日韩欧美免费精品| 99精品久久久久人妻精品| 看免费av毛片| 麻豆乱淫一区二区| 1024香蕉在线观看| 国产精品 国内视频| 久久久国产一区二区| a级毛片黄视频| 在线观看人妻少妇| 国产精品自产拍在线观看55亚洲 | cao死你这个sao货| 美女扒开内裤让男人捅视频| 亚洲va日本ⅴa欧美va伊人久久 | a级片在线免费高清观看视频| 一级片免费观看大全| 美女扒开内裤让男人捅视频| 国产视频一区二区在线看| 人妻 亚洲 视频| 18在线观看网站| 欧美中文综合在线视频| 手机成人av网站| 欧美+亚洲+日韩+国产| 精品久久久久久久毛片微露脸 | 久久综合国产亚洲精品| 久久人人爽人人片av| 国产免费av片在线观看野外av| 99国产精品一区二区蜜桃av | avwww免费| a级毛片在线看网站| 亚洲国产av影院在线观看| 国产99久久九九免费精品| 每晚都被弄得嗷嗷叫到高潮| 精品人妻熟女毛片av久久网站| 在线观看免费日韩欧美大片| 久久午夜综合久久蜜桃| 日韩一区二区三区影片| 考比视频在线观看| 汤姆久久久久久久影院中文字幕| 久久影院123| 热re99久久精品国产66热6| 91精品三级在线观看| 超色免费av| 国产日韩一区二区三区精品不卡| 又紧又爽又黄一区二区| 精品乱码久久久久久99久播| 国产黄色免费在线视频| 亚洲欧洲精品一区二区精品久久久| 老熟妇乱子伦视频在线观看 | 天堂中文最新版在线下载| 十分钟在线观看高清视频www| 欧美 亚洲 国产 日韩一| 搡老岳熟女国产| 叶爱在线成人免费视频播放| 欧美黄色淫秽网站| 久久精品亚洲熟妇少妇任你| 国产xxxxx性猛交| 纵有疾风起免费观看全集完整版| 最新在线观看一区二区三区| 日韩免费高清中文字幕av| 少妇精品久久久久久久| 国产精品av久久久久免费| 欧美日韩亚洲国产一区二区在线观看 | 91精品国产国语对白视频| 亚洲精品国产av蜜桃| 成在线人永久免费视频| 中文字幕色久视频| 亚洲午夜精品一区,二区,三区| 精品一区二区三卡| 久久久久久久久久久久大奶| 国产精品 欧美亚洲| 亚洲中文字幕日韩| 国产欧美亚洲国产| 日本撒尿小便嘘嘘汇集6| 国产精品熟女久久久久浪| 久久精品熟女亚洲av麻豆精品| 高清黄色对白视频在线免费看| 免费久久久久久久精品成人欧美视频| 国产精品99久久99久久久不卡| 久久亚洲国产成人精品v| 欧美日韩福利视频一区二区| 亚洲中文日韩欧美视频| 99久久人妻综合| 中文字幕人妻丝袜制服| 成人三级做爰电影| 午夜两性在线视频| 中文字幕人妻丝袜一区二区| 操出白浆在线播放| 一本综合久久免费| 一区福利在线观看| 男人操女人黄网站| 91精品伊人久久大香线蕉| 亚洲成人免费av在线播放| 欧美在线一区亚洲| 自拍欧美九色日韩亚洲蝌蚪91| 欧美激情久久久久久爽电影 | 多毛熟女@视频| 动漫黄色视频在线观看| 久久久国产精品麻豆| 老司机午夜福利在线观看视频 | 国产精品 欧美亚洲| 婷婷色av中文字幕| 亚洲中文av在线| 两人在一起打扑克的视频| 国产精品久久久久久精品电影小说| 另类亚洲欧美激情| 黄色片一级片一级黄色片| 日韩免费高清中文字幕av| av欧美777| 后天国语完整版免费观看| 大片免费播放器 马上看| 日韩 欧美 亚洲 中文字幕| 黄片大片在线免费观看| 精品高清国产在线一区| 在线观看免费视频网站a站| 亚洲精品自拍成人| 国产成人av激情在线播放| 嫩草影视91久久| 高清黄色对白视频在线免费看| 欧美另类一区| 亚洲精品国产区一区二| 男女之事视频高清在线观看| 丰满人妻熟妇乱又伦精品不卡| 美女视频免费永久观看网站| 一个人免费在线观看的高清视频 | 一级毛片精品| 国产日韩欧美亚洲二区| 久久精品成人免费网站| 啦啦啦中文免费视频观看日本| 日本91视频免费播放| 亚洲人成电影观看| 中文字幕高清在线视频| 男女无遮挡免费网站观看| 亚洲色图 男人天堂 中文字幕| 91九色精品人成在线观看| 一区福利在线观看| 美女高潮到喷水免费观看| 99精品久久久久人妻精品| 91九色精品人成在线观看| 美女午夜性视频免费| 桃红色精品国产亚洲av| 国产免费一区二区三区四区乱码| 一区福利在线观看| 国产极品粉嫩免费观看在线| 亚洲欧美日韩另类电影网站| 五月开心婷婷网| 男女国产视频网站| 12—13女人毛片做爰片一| 午夜福利在线免费观看网站| 777米奇影视久久| 国产精品一二三区在线看| 欧美日韩亚洲综合一区二区三区_| 国产成人啪精品午夜网站| 天天躁夜夜躁狠狠躁躁| 99国产精品一区二区蜜桃av | 成年动漫av网址| 欧美人与性动交α欧美精品济南到| 天天躁日日躁夜夜躁夜夜| 91麻豆精品激情在线观看国产 | 国产成人影院久久av| 男男h啪啪无遮挡| 国产主播在线观看一区二区| 热99re8久久精品国产| 欧美日韩av久久| 成人免费观看视频高清| 免费av中文字幕在线| 精品亚洲乱码少妇综合久久| av又黄又爽大尺度在线免费看| 两个人看的免费小视频| 国产精品影院久久| 满18在线观看网站| 亚洲精品国产一区二区精华液| 午夜91福利影院| 国产淫语在线视频| 亚洲精品自拍成人| 亚洲欧洲日产国产| 97在线人人人人妻| 日韩大码丰满熟妇| 男男h啪啪无遮挡| 男女无遮挡免费网站观看| 午夜精品久久久久久毛片777| 老司机亚洲免费影院| 国产一区二区三区综合在线观看| 国产亚洲欧美精品永久| 好男人电影高清在线观看| 精品福利永久在线观看| 精品卡一卡二卡四卡免费| 日韩欧美一区二区三区在线观看 | 久久热在线av| 在线观看人妻少妇| 国产色视频综合| 美女国产高潮福利片在线看| 国产精品九九99| 亚洲综合色网址| 欧美+亚洲+日韩+国产| 男女边摸边吃奶| 一本综合久久免费| 捣出白浆h1v1| 日本猛色少妇xxxxx猛交久久| 韩国高清视频一区二区三区| 男女之事视频高清在线观看| 国精品久久久久久国模美| 在线永久观看黄色视频| 在线av久久热| 黄色 视频免费看| 丝袜人妻中文字幕| 五月天丁香电影| 亚洲精品国产精品久久久不卡| 大型av网站在线播放| 久久久久久久精品精品| 国产欧美日韩一区二区三 | 国产精品一区二区免费欧美 | 菩萨蛮人人尽说江南好唐韦庄| 激情视频va一区二区三区| 色94色欧美一区二区| 久久99热这里只频精品6学生| 国产av一区二区精品久久| 一本—道久久a久久精品蜜桃钙片| 欧美 日韩 精品 国产| 国产成人影院久久av| 亚洲午夜精品一区,二区,三区| 亚洲精品粉嫩美女一区| 成人影院久久| 国产一区有黄有色的免费视频| 99热国产这里只有精品6| 性色av一级| 亚洲黑人精品在线| 日韩免费高清中文字幕av| 国产精品九九99| 亚洲av电影在线观看一区二区三区| 国产一区有黄有色的免费视频| 一区二区av电影网| 国产成人欧美| 桃红色精品国产亚洲av| 黄色视频在线播放观看不卡| 男女之事视频高清在线观看| 黄色毛片三级朝国网站| 亚洲国产精品999| 欧美国产精品一级二级三级| av天堂在线播放| 亚洲国产欧美一区二区综合| av电影中文网址| 一二三四社区在线视频社区8| 婷婷丁香在线五月| 男女午夜视频在线观看| 视频区图区小说| 每晚都被弄得嗷嗷叫到高潮| 免费人妻精品一区二区三区视频| 曰老女人黄片| 肉色欧美久久久久久久蜜桃| 免费一级毛片在线播放高清视频 | 黄色怎么调成土黄色| 人妻 亚洲 视频| 欧美日韩亚洲国产一区二区在线观看 | 国产精品国产三级国产专区5o| 电影成人av| 热99国产精品久久久久久7| 久久精品人人爽人人爽视色| 精品少妇久久久久久888优播| 大陆偷拍与自拍| 国产一区二区 视频在线| av视频免费观看在线观看| 欧美另类亚洲清纯唯美| 精品国产一区二区三区四区第35| 少妇的丰满在线观看| 久久精品国产a三级三级三级| 在线观看免费午夜福利视频| 悠悠久久av| 后天国语完整版免费观看| 啦啦啦啦在线视频资源| 91国产中文字幕| 悠悠久久av| 欧美人与性动交α欧美精品济南到| 亚洲国产欧美日韩在线播放| 中文欧美无线码| 97精品久久久久久久久久精品| 国产成人精品无人区| 两性午夜刺激爽爽歪歪视频在线观看 | 日本vs欧美在线观看视频| 国产深夜福利视频在线观看| 欧美日韩一级在线毛片| 国产高清视频在线播放一区 | 777久久人妻少妇嫩草av网站| 天天操日日干夜夜撸| 久久久久网色| 大型av网站在线播放| 老司机亚洲免费影院| 亚洲欧美精品自产自拍| 欧美国产精品va在线观看不卡| 欧美激情高清一区二区三区| 搡老熟女国产l中国老女人| 黄色视频,在线免费观看| 9191精品国产免费久久| 天天添夜夜摸| 亚洲专区字幕在线| 男女之事视频高清在线观看| 老司机影院成人| 亚洲一码二码三码区别大吗| 岛国在线观看网站| cao死你这个sao货| 午夜激情av网站| 亚洲一区二区三区欧美精品| 欧美成人午夜精品| 一区二区三区精品91| 一边摸一边抽搐一进一出视频| 欧美xxⅹ黑人| 曰老女人黄片| 最近最新免费中文字幕在线| 一级黄色大片毛片| 91国产中文字幕| 每晚都被弄得嗷嗷叫到高潮| 精品国产乱码久久久久久小说| 国产深夜福利视频在线观看| 99热网站在线观看| 欧美av亚洲av综合av国产av| 久久久久国产一级毛片高清牌| 午夜影院在线不卡| 老司机福利观看| 考比视频在线观看| 天堂8中文在线网| 国产欧美亚洲国产| 亚洲av成人一区二区三| 亚洲九九香蕉| 色老头精品视频在线观看| 麻豆av在线久日| 在线观看www视频免费| 国产一区二区三区av在线| 欧美 日韩 精品 国产| 99久久国产精品久久久| 新久久久久国产一级毛片|