銑削過程顫振穩定性分析的研究進展

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

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

?

銑削過程顫振穩定性分析的研究進展

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

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

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

關鍵詞：銑削；顫振；穩定性分析；微細銑削

銑削加工過程中刀具與工件間的相對振動，是降低產品的精度及表面質量，影響生產效率的主要原因。顫振是在切削過程中，在無周期性外部激振力持續作用下由于加工系統本身特性所引起的一種切削振動，一般表現為金屬切削過程中刀具與工件之間強烈的相對振動。顫振的發生降低了切削用量和工件表面質量，產生大量的噪聲，甚至會導致刀具的提前報廢。

通過銑削過程顫振穩定性分析，可獲取無顫振切削參數組合，以達到優化工藝參數，抑制振動的目的。國內外學者在銑削過程顫振穩定性分析領域的方法大多為從銑削過程動力學模型出發求得由切削參數構成的加工穩定域邊界。本文不僅介紹了銑削顫振穩定性分析方法中幾個關鍵環節的國內外研究成果，而且對近年來興起的微細銑削研究領域中顫振穩定性分析的研究進展進行了詳細說明。

1銑削過程顫振建模

1.1銑削過程動力學模型

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

(1) 摩擦型顫振

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

(2) 振型耦合型顫振

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

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

(3) 再生型顫振

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

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

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

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

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

上述研究說明了銑削過程中的非線性，如切削厚度的時滯特性、大振幅的斷續切削、多重再生效應等對于穩定性預測及切削顫振的影響。此外，切削系統結構的非線性，如刀具結構，機床結構的動態特性，也都是影響切削顫振模型的重要因素。

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

鋸齒銑刀由于其良好的抑振效果得國內外學者的深入研究。Dombovari等[30]在時頻域內建立了鋸齒銑刀切削穩定性模型。由于沿切削刃每一點的時滯都與鋸齒結構和進給率有關，利用含周期系統矩陣的多時滯微分方程組來描述切削動態過程。研究表明：在切削相同體積的材料時，鋸齒銑刀需要的驅動扭矩更小；即使在切削鈦之類的難加工材料時，鋸齒銑刀可采用無顫振切削深度也更大，加工效率更高。

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

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

2顫振穩定性分析

銑削穩定性分析的研究內容一般為預測由臨界加工參數構成的銑削過程穩定性邊界，從而劃分出穩定切削區和不穩定切削區，繼而可以從穩定性葉瓣圖中選擇適當的切削參數，達到避免顫振，提高加工效率和加工質量的目的。

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

2.1顫振模型求解法

顫振模型求解根據其對顫振模型的求解方法又可分為頻域法，離散法和數值法。

(1) 頻域法

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

(2) 離散法

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

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

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

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

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

(3) 數值法

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

2.2實驗法

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

3微細銑削過程顫振穩定性分析

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

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

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

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

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

4結論

本文對銑削過程顫振穩定性分析中動力學建模和穩定性預測的研究進行了歸納和綜述，并對微細銑削領域的顫振穩定性分析的研究進行了論述。

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

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

參 考 文 獻

[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] 劉習軍，陳予恕. 機床速度型切削顫振的非線性研究[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] 劉習軍，王立剛，賈啟芬. 一種由干摩擦引起的車床切削顫振[J]. 工程力學,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] 林潔瓊，周曉勤，孔繁森，等. 剛度主軸方位對模態耦合再生切削系統動態響應諧參數的影響[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] 吳雅. 機床切削系統的顫振及其控制[M]. 北京: 科學出版社, 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] 師漢民. 關于機床自激振動的一個非線性理論模型(第一部分)[J]. 應用力學學報,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] 師漢民. 關于機床自激振動的一個非線性理論模型(第二部份)[J]. 應用力學學報,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] 李勤良，汪博，趙斌，等. 考慮非線性遲滯力的機床顫振系統穩定性研究[J]. 機械工程學報, 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截斷的銑削系統穩定性的半離散分析法[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] 遲玉倫, 李郝林. 銑削顫振穩定域葉瓣圖確定方法研究[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.

基金項目：國家自然科學基金(51305061 )；中央高校基本科研業務費專項資金資助(DUT13LAB13)

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

中圖分類號:TH113；TG54

文獻標志碼：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年生