徐 磊，趙麗梅，余 炯

(貴州大學機械工程學院，貴陽 550025)

0 引言

隨著工業(yè)的發(fā)展，原有的零件加工方式因加工精度不高和效率低下，已無法滿足工業(yè)的需求，無心磨削方式的出現(xiàn)為工業(yè)零件的高精度和高效率生產(chǎn)提供了解決方案，目前，無心磨床已經(jīng)廣泛應(yīng)用于航空航天、汽車、冶金、電力等領(lǐng)域圓柱形、無中心孔的短軸和圓錐體等零件的生產(chǎn)[1]。我國對無心磨床的研究起步較晚，近年來，為突破國外的技術(shù)封鎖，縮小與國外在無心磨削領(lǐng)域的差距，國內(nèi)研究人員在無心磨床的自主研發(fā)方面做了大量的研究工作，但在磨削精度和效率方面與國外的無心磨床相比仍然存在較大的差距。其中造成磨床加工精度低下的主要原因在于砂輪磨削工件時發(fā)生顫振[2]。因此，從無心磨床磨削過程再生顫振機理出發(fā)，分析無心磨削的穩(wěn)定性，來提高無心磨床磨削穩(wěn)定性，從而提高加工精度是很有必要的[3]。

縱觀國內(nèi)外對機床切削顫振的研究，有很多研究人員做了很多卓有成效的工作。文獻[4]研究了單自由度再生顫振模型，通過磨削力表達式的推導，結(jié)合振動方程分析磨削穩(wěn)定性，并引入了同時考慮工件及砂輪再生顫振的雙再生效應(yīng)理論。文獻[5]通過建立的切入磨再生顫振動力學模型，利用工件表面波紋度幅值的增長率為穩(wěn)定性判據(jù)，研究再生顫振對磨削工件系統(tǒng)穩(wěn)定性的影響。文獻[6]在分析外圓磨削過程的顫振增長機理的基礎(chǔ)上，提出了再生顫振原理的外圓磨削動力學建模方法，來研究外圓磨削的穩(wěn)定性。文獻[7]推導了機床再生型切削顫振系統(tǒng)極限切削寬度隨主軸轉(zhuǎn)速變化的計算公式，提出了機床切削系統(tǒng)穩(wěn)定性極限預測方法。文獻[8]以外圓磨床磨削過程為研究對象，提出了一種用于雙再生外圓磨削過程穩(wěn)定性分析的系統(tǒng)方法，能夠直接從頻譜中獲取參數(shù)空間的穩(wěn)定特性，而不需要使用時域響應(yīng)。以上研究均取得了一定的研究成果，但在對機床顫振的研究中，主要是對車床、銑床和普通外圓磨床的研究，對無心外圓磨床的研究則涉及較少，而且，在建立顫振模型時，沒有考慮砂輪在磨削工件時產(chǎn)生的摩擦力，同時忽視了再生的磨削力和摩擦力對工件表面成形的影響，因此得出的結(jié)論具有很大的局限性，還無法指導實際生產(chǎn)。

本文以高液靜壓無心磨床無心磨削系統(tǒng)為研究對象，在分析高液靜壓無心磨床結(jié)構(gòu)工作原理和實際加工的基礎(chǔ)上，來研究無心磨削的顫振機理；抓住高液靜壓無心磨床磨削過程的特殊性，設(shè)定了特殊的坐標系，建立了同時考慮砂輪和工件的再生和摩擦效應(yīng)的無心磨削再生顫振數(shù)學模型，并計算穩(wěn)定性分析的特征值，結(jié)合數(shù)值延拓算法得到穩(wěn)定性葉瓣圖，為選擇合適的加工工藝參數(shù)來避免顫振提供理論依據(jù)。

1 高液靜壓無心磨床結(jié)構(gòu)原理分析及磨削過程再生顫振機理研究

1.1 無心磨床的結(jié)構(gòu)原理分析

高液靜壓無心磨床的工作原理分析對其工作過程的顫振機理研究和穩(wěn)定性分析具有至關(guān)重要的作用，因此，對高液靜壓無心磨床進行結(jié)構(gòu)分析和工作原理分析是十分必要的。本文以貴州省某廠生產(chǎn)的HFC-1808型高液靜壓無心磨床為研究對象，磨床的裝配結(jié)構(gòu)圖如圖1所示。高液靜壓無心磨床由磨削工件的砂輪，支承并帶動工件旋轉(zhuǎn)的導輪，支承工件的托板，對砂輪的尺寸、形狀、幾何角度進行修整的砂輪修整器，對導輪的形狀、尺寸和幾何角度等進行修整的導輪修整器，滑臺，導軌，機座等組成[9]。待加工的工件放在砂輪、導輪和托板之間，由托板和導輪進行支承，在導輪的帶動下旋轉(zhuǎn)，并通過砂輪對工件進行磨削；砂輪固定在機座上，導輪和托板安裝在滑臺上，滑臺由導軌支承，通過滑臺在導軌上的移動來實現(xiàn)工件的進給，進而達到加工精度，很明顯，高液靜壓無心磨床的磨削方式為切入式進給磨削。

1.砂輪修整器 2.砂輪 3.導輪修整器 4.滑臺5.導軌 6.導輪 7.工件 8.托板 9.機座圖1 高液靜壓無心磨床裝配結(jié)構(gòu)圖

1.2 無心磨削過程的再生顫振機理研究

在高液靜壓無心磨床切入式無心磨削過程當中，磨削輪用磨粒切割并再生工件表面，導輪帶動工件旋轉(zhuǎn)，并推動工件實現(xiàn)徑向進給，變相的使固定不動的砂輪完成徑向和圓周方向的進給。當工件被砂輪連續(xù)摩擦時，工件與砂輪之間產(chǎn)生的磨削力相互推離。根據(jù)再生顫振理論可知，磨削力和瞬時磨削深度之間存在著正比的關(guān)系[10]。磨削深度的任何變動都會引起磨削力的變化，因此時變的磨削力會激發(fā)砂輪與工件的振動，砂輪和工件的位移會影響磨削深度和磨削力，只要系統(tǒng)阻尼不足，這個狀態(tài)就會保持，這種自激振動稱為再生顫振。因此，為了更好了解再生效應(yīng)的磨削動力學，提出了一種基于再生理論的動力學數(shù)學模型。由于托板的表面是由耐磨損材料制作而成，故在高液靜壓無心磨床磨削過程當中，托板的磨損很小，對工件的再生效應(yīng)可以忽略不計。另外，導輪是調(diào)整輪，控制工件的運動，沒有參與工件的磨削，因此，對工件再生表面的影響忽略不計。

2 高液靜壓無心磨床磨削過程再生顫振動力學模型

2.1 無心磨削過程的數(shù)學模型

高液靜壓無心磨床磨削示意圖如圖2所示，Rg為砂輪的直徑，m；Mg為質(zhì)量，kg；Cg為阻尼系數(shù)，N/(m/s)；Kg為剛度，N/m；ωg為磨削工件的旋轉(zhuǎn)速度，rad/s；Rw為工件的直徑，m；Mw為質(zhì)量，kg；Cw為阻尼系數(shù)，N/(m/s)；Kw為剛度，N/m；ωw為旋轉(zhuǎn)速度，rad/s；工件在導輪和托板的推動下，以進給速度f(m/rev)向砂輪進給。以砂輪的回轉(zhuǎn)中心為原點，砂輪和工件的公法線方向為X軸，垂直于公法線方向為Y軸建立坐標系?；谒⒌淖鴺讼?，砂輪磨削工件的動力學控制方程可以離散為：

(1)

圖2 高液靜壓無心磨床磨削示意圖

2.2 切入式無心磨削磨削力

在高液靜壓無心磨床磨削工件的過程當中，給磨削系統(tǒng)帶來不穩(wěn)定的最主要的原因主要是砂輪和工件之間的磨削力，結(jié)合外圓磨床磨削力理論，提出了一種針對于切入式無心磨床砂輪磨削工件的磨削力模型；磨削力可以分解為垂直工件徑向的切向力Ft、沿工件徑向的徑向力Fn和沿工件軸向的軸向力Fa，其中軸向力最小，可以忽略不計，切向磨削力對砂輪的影響較小，為簡化問題，在砂輪的動力學模型中，將切向磨削力忽略，只考慮徑向磨削力對砂輪的影響，結(jié)合式(1)可知，F(xiàn)x和Fy分別表示工件在X方向和Y方向的激振力，并分別等于工件的徑向磨削力Fn和切向磨削力Ft。

Fx=Fn
Fy=Ft

(2)

結(jié)合高液靜壓無心磨床的工作原理可知，該無心磨床是利用砂輪磨削工件，而砂輪作為磨料磨削工具，由其上的結(jié)合劑固定的小顆粒磨料來研磨工件，從而產(chǎn)生時變的磨削力[11]。砂輪磨削工件表現(xiàn)出的外在的磨削力，實際上是由無數(shù)個單個磨粒磨削工件產(chǎn)生的切削力和滑擦力的合力，據(jù)此提出了圖3所示的單個磨粒磨削工件產(chǎn)生切削力和摩擦力的示意圖，由圖可知，在高液靜壓無心磨床磨削工件的過程中，砂輪的每個磨粒對工件表面進行切削和摩擦，并產(chǎn)生切向和法向的切削力和摩擦力[10]。

fct=ψfcn
fft=μffn

(3)

式(3)中，fct和fcn分別表示切削力在切向和徑向的分量，fft和ffn分別表示滑擦力在切向和徑向的分量，μ是摩擦系數(shù)，ψ是一種取決于磨粒平均尖端角的無量綱比。

(a) 砂輪磨削工件示意圖 (b) 單個磨粒磨削示意圖 (c) 磨削中各力分量圖3 單個磨粒產(chǎn)生切削力和摩擦力示意圖

在高液靜壓無心磨床磨削工件的過程中，將所有參與磨削的磨粒算在內(nèi)計算磨削力，得到總的徑向和切向磨削力如式(2)所示[13]：

(4)

其中，Kc是切削過程系數(shù)，N/m2；Kf是滑擦過程系數(shù)，N/m3/2；W為工件的寬度；Ds(t)為取決于磨削進給和砂輪-工件相對位移的瞬時磨削深度。

根據(jù)再生理論來確定磨削深度Ds(t)，同時考慮工件進給和工件與砂輪的相對位移，則磨削深度為：

Ds(t)=fcos(βg)+Xw(t)-Xg(t)-εYw(t)
(0<ε<1)

(5)

磨削是一個工件材料連續(xù)被砂輪移除生成新表面的過程，因此，在先前(t-Tw)時工件和砂輪相對位移的任何變動，都會影響瞬時磨削深度Ds(t)，則磨削深度Ds(t)變?yōu)椋?/p>

Ds(t)=fcos(βg)+(Xw(t)-Xg(t)-εYw(t))-

Xw(t-Tw)-Xg(t-Tw)-εYw(t-Tw)

(6)

而且，在磨削過程中，砂輪是不斷磨損生成新表面的過程，因此，在雙再生效應(yīng)的影響下，磨削深度Ds(t)變成：

Ds(t)=fcos(βg)+(Xw(t)-Xg(t)-εYw(t))-

Xw(t-Tw)-Xg(t-Tw)-εYw(t-Tw)-

g(Xw(t-Tg)-Xg(t-Tg)-εYw(t-Tg))

(7)

除了再生效應(yīng)，利用一個符號函數(shù)來考慮滑擦，則式(4)中砂輪和工件之間的摩擦系數(shù)μ是切向相對速度Vf的函數(shù)。

μ=sign(Vf)μd

(8)

如圖4所示，砂輪和工件之間的切向相對速度Vf取決于砂輪和工件的旋轉(zhuǎn)，二者的平面運動，即：

Vf=Vfg-Vfw

(9)

圖4 砂輪和工件在接觸區(qū)切向相對速度示意圖

式中，Vfg和Vfw分別是砂輪和工件在接觸區(qū)的切向速度：

(10)

3 高液靜壓無心磨床磨削過程穩(wěn)定性分析

在1965年，自Tobias首次提出穩(wěn)定性葉瓣圖的概念開始[14]，穩(wěn)定性葉瓣圖逐步發(fā)展并作為一種在切削加工規(guī)劃階段用來預測切削穩(wěn)定性和用來解釋機床動態(tài)切削過程行為的一種方法，并被大家廣泛接受。但是，在穩(wěn)定性葉瓣圖的構(gòu)建當中，人們只熱衷于車床或銑床切削過程穩(wěn)定性葉瓣圖的構(gòu)建，對切入式無心磨床磨削過程穩(wěn)定性葉瓣圖的構(gòu)建，幾乎沒有涉獵。因此本節(jié)嘗試從構(gòu)建切入式無心磨床磨削過程穩(wěn)定性葉瓣圖的角度出發(fā)，來研究高液靜壓無心磨床磨削過程穩(wěn)定性。

在高液靜壓無心磨床穩(wěn)定磨削的過程中，令X(t)=[xg(t)xw(t)yw(t)]T，將式(4)、式(8)～式(10)代入式(1)，可以得到如下磨削模型:

(11)

本文以MK1808型高液靜壓無心磨床為研究對象，將要用于計算和分析的參數(shù)值列于表1。

表1 穩(wěn)定性仿真參數(shù)值

通過MATLAB編程對式(11)來進行特征值計算并結(jié)合數(shù)值延拓算法得到HFC-1808型高液靜壓切入磨無心磨床的穩(wěn)定性葉瓣圖(如圖5所示)，該圖以砂輪轉(zhuǎn)速為橫坐標，工件的徑向磨削深度為縱坐標。

圖5 高液靜壓無心磨床磨削穩(wěn)定性葉瓣圖

圖5所示的砂輪再生顫振穩(wěn)定性葉瓣圖，由不同的臨界點連接而成，并將整個平面分成兩個區(qū)域，每個區(qū)域具有不同的穩(wěn)定性性質(zhì)，依據(jù)穩(wěn)定性性質(zhì)的不同可分成4個穩(wěn)定性區(qū)域。

(1)切入式無心磨削顫振區(qū)：當砂輪轉(zhuǎn)速和徑向磨削深度匹配點完全位于曲線上方時，如圖中的G1點，砂輪要發(fā)生再生顫振，這時無論怎樣選擇砂輪轉(zhuǎn)速和磨削深度，砂輪磨削工件都會發(fā)生顫振，從而引發(fā)一系列的危害，因此，切入式無心磨削顫振區(qū)應(yīng)堅決避免。

(2)切入式無心磨削臨界穩(wěn)定區(qū)：當砂輪轉(zhuǎn)速和徑向磨削深度匹配點剛好落在曲線上時，如圖中G2點，為切入式無心磨削臨界穩(wěn)定區(qū)，此時，砂輪磨削工件可能會發(fā)生顫振，也可能穩(wěn)定的磨削，總之，只要顫振發(fā)生，就會嚴重影響工件的磨削精度，因此，切入式無心磨削臨界穩(wěn)定區(qū)也是應(yīng)該避免的區(qū)域。

(3)切入式無心磨削有條件穩(wěn)定區(qū)：當砂輪轉(zhuǎn)速和徑向磨削深度匹配點落在極限磨削點劃線hlim和葉瓣圖曲線之間時，如圖中G3點，為切入式無心磨削有條件穩(wěn)定區(qū)；當砂輪轉(zhuǎn)速固定不變，隨著徑向磨削深度的增加，匹配點先后經(jīng)過穩(wěn)定區(qū)和顫振區(qū)，砂輪先穩(wěn)定后發(fā)生顫振；當固定徑向磨削深度，隨著砂輪轉(zhuǎn)速的增加，穩(wěn)定區(qū)和顫振區(qū)相間出現(xiàn)。

(4)切入式無心磨削無條件穩(wěn)定區(qū)：當砂輪轉(zhuǎn)速和徑向磨削深度匹配點落在極限磨削點劃線以下時，如圖中G4點，為切入式無心磨削無條件穩(wěn)定區(qū)，在此區(qū)域，無論如何選擇砂輪轉(zhuǎn)速和徑向磨削深度，砂輪都不會發(fā)生顫振。

綜上所述，結(jié)合高液靜壓無心磨床穩(wěn)定性葉瓣圖，可以選取適當工藝參數(shù)，從而避免砂輪磨削工件時顫振的發(fā)生(顫振狀態(tài)即為不穩(wěn)定狀態(tài))，從而達到提高加工精度和生產(chǎn)效率的目的。并且，高液靜壓無心磨床磨削穩(wěn)定性葉瓣圖為企業(yè)在生產(chǎn)加工時的參數(shù)選擇提供了依據(jù)，同時，也為相關(guān)磨床的穩(wěn)定性分析提供了新的思路。

4 結(jié)論

(1)應(yīng)用再生顫振理論，得出高液靜壓無心磨床磨削工件時的不穩(wěn)定原因主要是因為砂輪和工件發(fā)生了再生顫振引起的。

(2)以再生顫振機理為基礎(chǔ)，建立了同時考慮砂輪和工件再生顫振的動力學模型。

(3)求解動力學模型矩陣方程，利用MATLAB編程得到高液靜壓無心磨床切入式無心磨削穩(wěn)定性葉瓣圖，并根據(jù)砂輪轉(zhuǎn)速和徑向磨削深度匹配點的不同，將穩(wěn)定性葉瓣圖劃分為4個區(qū)域，為磨削工藝參數(shù)的選擇提供了理論依據(jù)。