馮吉路，孫志禮，孫洪哲，柴小冬

(東北大學 機械工程與自動化學院，沈陽　110819)

?

銑削顫振及其可靠度分析*

馮吉路，孫志禮，孫洪哲，柴小冬

(東北大學 機械工程與自動化學院，沈陽110819)

摘要：通過銑削顫振動力學建模，考慮了影響參數(shù)的隨機性，在解析法分析銑削系統(tǒng)阻尼比、剛度、固有頻率、銑削徑向切削寬度以及銑削力系數(shù)對銑削顫振穩(wěn)定性的影響的基礎上，應用蒙特卡洛算法，分別對不同阻尼比和剛度的銑削系統(tǒng)進行了銑削顫振可靠性分析。結果表明，主軸轉(zhuǎn)速和銑削顫振可靠度呈現(xiàn)明顯的非線性關系。在穩(wěn)定域按照極限切削深度進行銑削加工時，僅有在特定速度段能夠保障不發(fā)生顫振現(xiàn)象，剛度的增加幾乎不會影響銑削顫振可靠度，增加銑削系統(tǒng)阻尼可容易地得到銑削穩(wěn)定域所對應的主軸速度。

關鍵詞：銑削顫振；極限切削深度；可靠性；蒙特卡洛

0引言

銑削加工是一種最常見和最通用的加工方式，其具有良好的材料去除性能，因此在工業(yè)生產(chǎn)中扮演了重要的角色[1]。然而，在某些條件下銑削加工過程中會產(chǎn)生顫振現(xiàn)象，這不僅會引起銑削系統(tǒng)的自激振動和受迫振動，從而導致加工過程中的噪聲加劇，并且會影響被加工工件的表面質(zhì)量和刀具的使用壽命，嚴重的影響加工效率[2]。因此，研究人員對銑削加工過程中的顫振預測和顫振抑制進行了廣泛的研究。

Tobias和Tlusty等[3-4]首先提出了切削過程中的顫振穩(wěn)定性預測模型，建立了單自由度特征方程。Altintas和Budak[5]提出了銑削過程中二自由度的顫振解析模型，該模型建立過程中假設系統(tǒng)的動態(tài)參數(shù)不變且切削系數(shù)是時不變的。遲玉倫等[6]對銑削顫振穩(wěn)定域葉瓣圖確定方法進行了研究，得到了與實驗測量結果相一致的銑削顫振穩(wěn)定域葉瓣圖。劉戰(zhàn)強等[7]通過分析斜角切削的幾何關系，分析了銑刀參數(shù)對銑削穩(wěn)定性的影響。祝長生等[8]應用內(nèi)置力執(zhí)行器對銑削顫振進行了主動控制，并在抑制銑削顫振取得了明顯效果。李茂月等[9]通過變主軸轉(zhuǎn)速在線抑制顫振算法實現(xiàn)了抑制銑削顫振。Jui-Pin Hung等[10]分析了主軸系統(tǒng)對銑削刀具動態(tài)性能的影響，分析得到結構剛度和軸承剛度是影響銑削動態(tài)性能的主要因素之一。銑削系統(tǒng)的剛度、阻尼比、固有頻率以及材料切削力系數(shù)和徑向切削深度等值在實際實驗測試時，均為隨機變量。然而，上述研究并沒有考慮銑削系統(tǒng)動力學參數(shù)的隨機性對銑削顫振的影響，這樣會導致即便按照極限切削深度進行加工時，仍然發(fā)生顫振現(xiàn)象的情況，因此有必要對銑削系統(tǒng)進行顫振可靠性分析。

本文通過銑削動力學建模仿真，分別分析了各個參數(shù)對銑削顫振的影響，得到了銑削顫振穩(wěn)定性葉瓣圖。考慮了銑削系統(tǒng)動力學參數(shù)的隨機性，通過蒙特卡洛算法對銑削可靠性進行了分析，為實現(xiàn)銑削顫振的主動抑制提供了有效依據(jù)。

1銑削顫振動力學建模

在銑削過程中，前一個刀齒切削留下了波紋，由下一個刀齒進行去除，并在新的表面上形成波紋，兩個波紋之間的相位移會使得切屑的厚度產(chǎn)生變化，從而使得切削力發(fā)生變化。當銑削力變化的頻率接近銑削系統(tǒng)的某階模態(tài)時，系統(tǒng)會發(fā)生顫振現(xiàn)象。

主軸-銑刀-工件所構成的銑削動態(tài)模型如圖1所示。銑削系統(tǒng)的剛度和阻尼由主軸-刀具系統(tǒng)決定，模型中銑刀具有4個刀齒，x方向為刀具的進給方向。銑削系統(tǒng)簡化模型其動力學微分方程可表示：

(1)

式中，mx、my、ξx、ξy、kx、ky分別為銑削系統(tǒng)X、Y方向上的質(zhì)量、阻尼和剛度，F(xiàn)x(t)和Fy(t)分別為刀齒上銑削力在X、Y方向上的分量。

銑削過程中動態(tài)銑削力可以表示：

(2)

式中，a為軸向切削深度，Kt為切削系數(shù)，[A0]為切削方向系數(shù)矩陣，Δ(t)為切削厚度變化。

圖1　銑削系統(tǒng)動力學模型

考慮再生效應的動態(tài)切削力表達式：

(3)

式中，[G(iw)]刀具-工件接觸區(qū)傳遞函數(shù)矩陣。

通過求解方程的特征根可得到顫振頻率為wc時的極限切削深度：

(4)

式中，ΛRe為特征方程的實部，κ為特征方程虛部與實部的比值，N為銑刀的齒數(shù)。

根據(jù)特征方程虛部與實部的比值可以求得主軸對應的轉(zhuǎn)速：

(5)

式中，k為葉瓣數(shù)。

2銑削顫振結果分析

通過對上述方程分析結果可知，銑削顫振穩(wěn)定性分別由銑削系統(tǒng)的動力學參數(shù)、刀具參數(shù)和材料屬性決定。為分析銑削動力學參數(shù)對顫振穩(wěn)定性的影響，假設銑削系統(tǒng)的X和Y方向的動力學參數(shù)完全相同。仿真時的參數(shù)分別：Wn=600Hz、ξ=0.02、K=12500N/mm、N=4、Kt=800N/mm2、Kn=170N/mm2、ae=2.0mm。

保持其它影響參數(shù)不變的前提下，改變銑削系統(tǒng)阻尼，分析其對顫振穩(wěn)定性的影響，得到銑削顫振穩(wěn)定域葉瓣圖如圖2所示。銑削穩(wěn)定性葉瓣圖隨銑削系統(tǒng)阻尼的增加整體上移，極限切削深度隨阻尼比的增加而增大。系統(tǒng)阻尼比為0.0125、0.0150、0.0175和0.0200所對應的極限切削深度分別為2.85mm、3.42mm、3.99mm和4.57mm。

圖2　銑削系統(tǒng)阻尼比對顫振穩(wěn)定性的影響

保持其它影響參數(shù)不變的前提下，改變銑削系統(tǒng)剛度，分析其對顫振穩(wěn)定性的影響，得到銑削顫振穩(wěn)定域葉瓣圖如圖3所示。銑削穩(wěn)定性葉瓣圖隨銑削系統(tǒng)剛度的增加整體上移，但與圖2相比，極限切削深度的峰值變化更為明顯。極限銑削深度隨系統(tǒng)剛度的增加而增大。銑削系統(tǒng)剛度為12500N/mm、15000N/mm、17500N/mm和20000N/mm所對應的極限切削深度分別為4.57mm、5.48mm、6.40mm和7.31mm。

圖3　銑削系統(tǒng)剛度對顫振穩(wěn)定性的影響

改變銑削系統(tǒng)固有頻率進行銑削顫振分析，得到銑削顫振穩(wěn)定域葉瓣圖如圖4所示。銑削系統(tǒng)的固有頻率不會影響極限切削深度，但會使得穩(wěn)定域葉瓣圖整體右移，顫振穩(wěn)定域的面積逐漸減小。

圖4　銑削系統(tǒng)固有頻率對顫振穩(wěn)定性的影響

圖5為改變徑向銑削寬度得到的銑削顫振穩(wěn)定域葉瓣圖。通過分析可得，穩(wěn)定域葉瓣圖會隨著徑向銑削寬度的增加整體下移。徑向銑削寬度的增加會使極限銑削深度減小。即銑削寬度越大，越不容易實現(xiàn)銑削穩(wěn)定。徑向銑削寬度為2.0mm、2.4mm、2.8mm和3.2mm所對應的極限切削深度分別為4.57mm、3.73mm、3.16mm和2.73mm。

圖5　徑向銑削寬度對顫振穩(wěn)定性的影響

如圖6所示，銑削力系數(shù)變化對顫振穩(wěn)定性的影響，該圖顯示極限切削深度與銑削力系數(shù)成反比。銑削力系數(shù)反應了材料的難加工程度，切削力系數(shù)越大，材料越難加工。分析可知，銑削力系數(shù)越大，材料加工時越容易產(chǎn)生顫振現(xiàn)象。材料的銑削力系數(shù)為600N/mm2、800N/mm2、1000N/mm2和1200N/mm2所對應的極限切削深度分別為6.09mm、4.57mm、3.65mm和3.05mm。

圖6　銑削力系數(shù)對顫振穩(wěn)定性的影響

3銑削顫振可靠度分析

在銑削加工過程中，切削深度超過極限切削深度時會發(fā)生顫振現(xiàn)象，進而導致被加工工件的表面質(zhì)量和刀具壽命降低，并會使得切削噪聲加劇。在計算得到極限切削深度的情況下，由于機床的動力學參數(shù)、材料的銑削力系數(shù)和銑削徑向切削寬度等參數(shù)具有隨機性，在使用極限切削深度進行銑削加工時仍然有可能發(fā)生顫振現(xiàn)象。

設銑削系統(tǒng)X、Y方向的動力學參數(shù)相同，銑削的極限切削深度由下式表達：

flim ap=y(wn,ξ,K,Kt,ae)

(6)

式中，wn，ξ，K，Kt，ae為隨機參數(shù)向量。假設隨機參數(shù)服從標準正態(tài)分布，其均值為muX=(muwn，muξ，muK，muKt，muae)，根據(jù)6σ原則，其方差為sigmaX=0.005muX。

根據(jù)銑削顫振定義可得銑削顫振的極限狀態(tài)函數(shù)：

Z=Aplim-flim ap

(7)

當Z≤0時，視為發(fā)生銑削顫振現(xiàn)象，當Z>0時，視為正常銑削加工狀態(tài)。

應用MonteCarlo法對銑削極限狀態(tài)函數(shù)進行n直接抽樣計算，其中Z>0的次數(shù)為s，因此，銑削顫振可靠度：

(8)

應用MonteCarlo算法進行107抽樣模擬，其中，模擬參數(shù)muX1=(600Hz，0.02，12500N/mm，800N/mm2，2.0mm)， muX2=(600Hz，0.05，12500N/mm，800N/mm2，2.0mm)； muX3=(600Hz，0.02，12500N/mm，800N/mm2，2.0mm)， muX4=(600Hz，0.05， 25000N/mm，800N/mm2，2.0mm)，方差根據(jù)6σ原則進行取值。模擬結果分別如圖7～圖10所示。

如圖7所示，銑削系統(tǒng)阻尼均值muξ=0.02時，不同主軸轉(zhuǎn)速所對應的銑削顫振可靠度。該圖中主軸轉(zhuǎn)速和可靠度在低速狀態(tài)下并沒有一一對應的關系。在穩(wěn)定域按照極限切削深度進行銑削加工時，僅有在特定速度段能夠保障不發(fā)生顫振現(xiàn)象。在轉(zhuǎn)速為10500r/min～11500r/min時，銑削顫振可靠度等于1，在該條件下銑削工件會獲得較好的表面質(zhì)量。在其它轉(zhuǎn)速下，不能保障不發(fā)生顫振現(xiàn)象。

與如圖7相對應，圖8為銑削系統(tǒng)阻尼均值muξ=0.05時，不同主軸轉(zhuǎn)速所對應的銑削顫振可靠度。該圖中主軸轉(zhuǎn)速和可靠度在全轉(zhuǎn)速下存在著一一對應關系，且銑削顫振可靠度在多個轉(zhuǎn)速段等于1。對比圖7和圖8的結果可得，提高銑削系統(tǒng)的阻尼可以實現(xiàn)對銑削顫振的主動抑制。

通過對比圖9和圖10可知，銑削系統(tǒng)剛度的變化不會對銑削可靠性產(chǎn)生影響。在轉(zhuǎn)速為5000r/min～8500r/min時，對應的極限切削深度進行銑削加工能夠保障不發(fā)生顫振現(xiàn)象。

圖7　銑削系統(tǒng)阻尼均值muξ=0.02時的銑削可靠性

圖8　銑削系統(tǒng)阻尼均值muξ=0.05時的銑削可靠性

圖9　銑削系統(tǒng)剛度均值muK=12500N/mm時的銑削可靠性

圖10　銑削系統(tǒng)剛度均值muK=25000N/mm時的銑削可靠性

4結論

通過數(shù)值模擬，分析了各個因素對銑削顫振穩(wěn)定性的影響，應用MonteCarlo進行銑削顫振可靠度的求解，得到如下結論：

(1)銑削顫振穩(wěn)定性隨銑削系統(tǒng)的阻尼和剛度的增加而增大，隨徑向銑削寬度和銑削力系數(shù)的增加而減小。銑削系統(tǒng)的固有頻率對銑削的極限切削深度沒有影響，但會使銑削顫振葉瓣圖整體右移。

(2)銑削系統(tǒng)剛度的變化不會對銑削顫振可靠度產(chǎn)生影響。

(3)在穩(wěn)定域按照極限切削深度進行銑削加工時，僅有在特定速度段能夠保障不發(fā)生顫振現(xiàn)象，且增加銑削系統(tǒng)阻尼可容易地得到銑削穩(wěn)定域所對應的主軸速度。

[參考文獻]

[1] Grossi N, Sallese L, Scippa A, et al. Chatter Stability Prediction in Milling Using Speed-varying Cutting Force Coefficients[J]. Procedia CIRP. 2014, 14: 170-175.

[2] Peng Z K, Jackson M R, Guo L Z, et al. Effects of bearing clearance on the chatter stability of milling process[J]. Nonlinear Analysis: Real World Applications, 2010, 11(5): 3577-3589.

[3] J Tlusty, M Polacek.The stability of the machine tool against self-excited vibration in machining[J]. Proceedings of the ASME International Research in Production Engineering,1963: 465-474.

[4] S A Tobias, W Fishwick.The chatter of lathe tools under orthogonal cutting conditions[J].Transactions of the ASME ,1958 (80): 1079-1088.

[5] Altinta. Analytical Prediction of Stability Lobes in Milling[J]. CIRP Annals - Manufacturing Technology. 1995, 44(1): 357-362.

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

[7] 湯愛君，劉戰(zhàn)強. 銑刀參數(shù)對薄壁零件銑削穩(wěn)定性的影響[J].華南理工大學學報(自然科學版), 2009, 37(2): 29-34.

[8] 喬曉利，祝長生. 基于內(nèi)置力執(zhí)行器的銑削顫振的主動控制[J]. 機械工程學報,2012,48(1): 185-192.

[9] 李茂月，韓振宇，富宏亞，等. 基于開放式控制器的銑削顫振在線抑制[J]. 機械工程學報,2012, 48(17): 172-182.

[10] Hung J, Lai Y,et al. Analysis of the machining stability of a milling machine considering the effect of machine frame structure and spindle bearings: experimental and finite element approaches[J]. The International Journal of Advanced Manufacturing Technology. 2013, 68(9-12): 2393-2405.

(編輯趙蓉)

Investigation of the Milling Chatter and its Probability

FENG Ji-lu，SUN Zhi-li，SUN Hong-zhe，CHAI Xiao-dong

(School of Mechanical Engineering & Automation, Northeastern University, Shenyang 110819, China)

Abstract：The analytical method was applied to account for the influence of the damping ratio，stiffness，natural frequency，radial milling width and milling force coefficient on the milling stability by the dynamical model of milling chatter. Taking the randomness above parameters into consideration，the Monte Carlo methoad was used to explain the milling chatter reliability in the different damping ratios and stiffness.Current results show that the relationship between the spindle speed and milling chatter reliability is clearly nonlinear. The increasing stiffness of milling system has no effect on the milling chatter reliability. The chatter phenomenon just only can not be occurred in the specific speed zone when machining at the zone of milling stability，which associated with the spindle speed can be easily obtained by increasing the damping ratio of milling system.

Key words：milling chatter;limit cutting depth;reliability;Monte Carlo

中圖分類號：TH166;TG506

文獻標識碼：A

作者簡介：馮吉路(1988—)，男，天津人，東北大學博士研究生，研究方向機床主軸動力學及加工動力學，(E-mail)fengjilu123@sina.com; 孫志禮(1957—),男,山東巨野縣人,東北大學教授，博導，研究方向為高檔數(shù)控機床性能分析與可靠性設計。

*基金項目：國家科技重大項目(2013ZX04011-011)

收稿日期：2015-03-02

文章編號：1001-2265(2015)12-0008-04

DOI:10.13462/j.cnki.mmtamt.2015.12.003