崔　立

(上海第二工業(yè)大學(xué)工學(xué)部，上海 201209)

?

基于預(yù)緊軸承動(dòng)剛度的高速電主軸動(dòng)特性分析*

上海市自然科學(xué)基金項(xiàng)目(14ZR1416800)；上海市聯(lián)盟計(jì)劃項(xiàng)目(LM201408)；國(guó)家自然科學(xué)基金資助(51305251、50905061)；上海第二工業(yè)大學(xué)重點(diǎn)學(xué)科“材料與工程”(XXKZD1601)

崔立

(上海第二工業(yè)大學(xué)工學(xué)部，上海 201209)

針對(duì)定位預(yù)緊的電主軸球軸承，基于球和套圈相互作用模型，推導(dǎo)了球軸承剛度求解模型，使用Newton-Raphson方法求解可到剛度矩陣?；赗iccati傳遞矩陣法，給出了考慮軸承動(dòng)剛度的電主軸臨界轉(zhuǎn)速計(jì)算方法。建立了銑削激振下考慮軸承動(dòng)剛度的電主軸動(dòng)力學(xué)方程，使用Runge-Kutta法求解得到電主軸動(dòng)力學(xué)響應(yīng)。對(duì)某銑削高速電主軸計(jì)算，結(jié)果發(fā)現(xiàn)軸承徑向剛度隨轉(zhuǎn)速增大而減小，采用定位預(yù)緊后剛度增大，且隨轉(zhuǎn)速增大軸心預(yù)緊剛度增大明顯?？紤]軸承動(dòng)剛度可得到更為準(zhǔn)確的電主軸動(dòng)特性結(jié)果，最后通過高速電主軸動(dòng)特性實(shí)驗(yàn)證明了該模型對(duì)動(dòng)特性預(yù)測(cè)的準(zhǔn)確性。

高速電主軸；動(dòng)剛度；球軸承；動(dòng)態(tài)特性；臨界轉(zhuǎn)速

電主軸的動(dòng)態(tài)特性指它抵抗受迫及自激振動(dòng)的能力。動(dòng)態(tài)特性包括固有振動(dòng)特性和響應(yīng)特性，這些特性對(duì)機(jī)床加工性能的影響很大，不僅影響機(jī)床的加工精度、使用壽命，也影響著其高效率性能的發(fā)揮。

為提高電主軸的精度，主軸上使用的滾動(dòng)軸承需預(yù)緊，軸承預(yù)緊后各滾動(dòng)體受力將更加均勻，可提高抗振性。目前關(guān)于考慮軸承預(yù)緊，研究電主軸性能的研究已有不少。Alfares等[1]建立了電主軸的五自由度動(dòng)力學(xué)方程，分析了軸承預(yù)緊力對(duì)電主軸動(dòng)力學(xué)性能影響，結(jié)果表明隨著初始預(yù)緊力的增大主軸的振動(dòng)量減小，但預(yù)緊力增大到一定程度時(shí)，主軸的振動(dòng)將不再明顯減小。Lin等[2]研究了在超高速情況下電主軸的動(dòng)力學(xué)特性，軸承預(yù)緊力對(duì)軸承剛度與主軸動(dòng)力學(xué)的影響，高轉(zhuǎn)速對(duì)軸承球離心力、陀螺力矩，進(jìn)而分析了對(duì)主軸整機(jī)動(dòng)力學(xué)的影響。Kang等[3]用有限元法，使用Ansys軟件分別完成機(jī)床主軸軸承系統(tǒng)靜態(tài)、動(dòng)力學(xué)特性與機(jī)床的顫振分析。趙志剛等[4]基于Ansys軟件對(duì)電主軸轉(zhuǎn)子軸承系統(tǒng)進(jìn)行了動(dòng)靜態(tài)性能分析，給出了軸承位置及跨距對(duì)電主軸動(dòng)特性影響。李松生等[5-6]用有限元法，建立了高速電主軸軸系轉(zhuǎn)子動(dòng)力學(xué)模型，分析了軸承支承特性隨速度變化的非線性特性和電主軸固有頻率。邱海飛等[7]設(shè)計(jì)了電主軸動(dòng)力學(xué)試驗(yàn)臺(tái)，確定了不同轉(zhuǎn)速段軸承最小預(yù)緊力。以上這些研究中，研究了預(yù)緊力對(duì)電主軸動(dòng)特性的影響，但未考慮預(yù)緊作用下軸承動(dòng)剛度的影響。

本文提出了考慮軸承動(dòng)剛度的電主軸動(dòng)特性計(jì)算方法，給出了電主軸臨界轉(zhuǎn)速和動(dòng)態(tài)響應(yīng)的計(jì)算方法，通過對(duì)高速電主軸動(dòng)特性測(cè)試驗(yàn)證模型的可靠性。

1　考慮動(dòng)剛度的電主軸臨界轉(zhuǎn)速分析

1.1考慮預(yù)緊的軸承剛度模型

軸承預(yù)緊方式有定位預(yù)緊和定壓預(yù)緊。由于定位預(yù)緊方式對(duì)增加軸承剛度更為有效，所以定位預(yù)緊廣泛應(yīng)用于大功率高速機(jī)床主軸上。

圖1a為球軸承初始狀態(tài)時(shí)球和內(nèi)外圈溝道曲率中心相對(duì)位置，內(nèi)外圈溝道的曲率中心分別為O2和O1，球的曲率中心為Ob，初始接觸角為α0。在外圈受到定位預(yù)緊作用下，外圈曲率中心從O1移動(dòng)到O1′，假設(shè)球軸承工作時(shí)受軸向載荷和徑向載荷作用，圖1b為受載狀態(tài)下，外圈曲率中心不變，內(nèi)圈曲率中心沿軸向和徑向移動(dòng)到O2′，球的曲率中心移動(dòng)到Ob′，α1i、α2i分別為受載后第i個(gè)滾動(dòng)體與外圈、內(nèi)圈的接觸角。

-(f1-0.5)Dw

(1)

-(f2-0.5)Dw

(2)

(3)

(4)

由赫茲接觸理論[8-9]可求得球和外圈、內(nèi)圈的接觸力Q1i、Q2i。

假設(shè)軸承所受載荷[Fi]={Fx,Fy,Fz}，建立球和內(nèi)圈的力平衡方程可得：

Q1isinα1i= Q2isinα2i

(5)

式中：ωoi是第i個(gè)球的公轉(zhuǎn)轉(zhuǎn)速，采用Newton-Raphson方法求解非線性方程組。軸承剛度可得

(6)

式中：Kb為剛度矩陣，求解方程組得到的雅可比矩陣即為軸承剛度矩陣。

1.2考慮動(dòng)剛度的電主軸臨界轉(zhuǎn)速計(jì)算

電主軸實(shí)際工作過程中應(yīng)避開臨界轉(zhuǎn)速，因此應(yīng)計(jì)算出電主軸的臨界轉(zhuǎn)速。而電主軸運(yùn)轉(zhuǎn)過程中由于軸承剛度會(huì)發(fā)生變化，為準(zhǔn)確計(jì)算電主軸臨界轉(zhuǎn)速，應(yīng)該考慮軸承動(dòng)剛度進(jìn)行計(jì)算。本文采用求解快速且精度較高的傳遞矩陣法。

將電主軸軸承系統(tǒng)劃分為若干桿單元和節(jié)點(diǎn)，建立動(dòng)力學(xué)方程如下：

(7)

式中：u為節(jié)點(diǎn)位移矩陣；M為質(zhì)量矩陣；C為阻尼矩陣；G為陀螺力矩陣；F為載荷矩陣；Ks為轉(zhuǎn)軸剛度矩陣[10]，Kb為軸承剛度矩陣。

電主軸運(yùn)轉(zhuǎn)過程中，當(dāng)軸承轉(zhuǎn)速或所受載荷變化時(shí)軸承剛度也將發(fā)生變化，因此應(yīng)考慮動(dòng)剛度，聯(lián)立式(6)與式(7)，使用Riccati傳遞矩陣法[10]求解可得到系統(tǒng)的頻率方程式：

(8)

式中：S為電主軸的Riccati傳遞矩陣。使用Newton-Raphson方法求解式(8)可得各階特征值：

s=λj+iωj

(9)

式中：λj為衰減因子；ωj為固有頻率。

在電主軸轉(zhuǎn)速范圍內(nèi)計(jì)算系統(tǒng)的特征值，得到固有頻率，當(dāng)固有頻率與電主軸轉(zhuǎn)速相等時(shí)，即為電主軸臨界轉(zhuǎn)速。

2　考慮動(dòng)剛度的電主軸響應(yīng)分析

銑床電主軸的激振力主要來(lái)自于銑刀的銑削力，圓柱銑刀的平均圓周切削力公式為：

(10)

式中：CF為銑削材料系數(shù)；ae為銑削深度；af為每齒進(jìn)給量；dq為銑刀直徑；ap為銑削寬度；Z為銑刀齒數(shù)。

銑削激振力可寫為：

q(t)=Fqcos(ωqt)

(11)

式中：ωq為激振力頻率，ωq=Zω；ω為軸頻。

銑削激振力是周期性的激振力，為分析銑削激振力作用下的動(dòng)態(tài)響應(yīng)，將式(11)中的銑削激振力g(t)替代式(7)，聯(lián)立式(5)求解軸承動(dòng)剛度矩陣Kb。使用Runge-Kutta法求解式(7)，得到電主軸各節(jié)點(diǎn)的動(dòng)力學(xué)響應(yīng)。

3　電主軸動(dòng)態(tài)性能分析與實(shí)驗(yàn)

3.1電主軸臨界轉(zhuǎn)速分析

以圖2所示的銑削電主軸為例進(jìn)行分析，電主軸總長(zhǎng)度為550 mm，前支承軸承采用3個(gè)型號(hào)相同的角接觸球軸承7014，后支承軸承采用兩個(gè)型號(hào)相同的角接觸球軸承7010，前后支承軸承都采用定位預(yù)緊，電主軸最高工作轉(zhuǎn)速為18 000 r/min。

對(duì)前支承軸承7014進(jìn)行剛度計(jì)算。圖3給出了在電主軸最高轉(zhuǎn)速18 000 r/min時(shí)，前支承軸承剛度隨預(yù)緊量變化曲線?？梢钥闯鲚S向剛度和徑向剛度均隨預(yù)緊量增大而增大，但增大幅度隨預(yù)緊量增大逐漸減小?？紤]到高速電主軸發(fā)熱量隨預(yù)緊量增大也會(huì)逐漸增大，選取較輕預(yù)緊，預(yù)緊量為10 μm。

當(dāng)定位預(yù)緊為10 μm時(shí)，圖4給出了在聯(lián)合載荷作用下球軸承徑向剛度隨隨轉(zhuǎn)速增大呈非線性減小，軸向剛度隨轉(zhuǎn)速增大而增大。軸承采用定位預(yù)緊后，軸向和徑向剛度與沒有預(yù)緊相比均增大，且隨轉(zhuǎn)速增大預(yù)緊剛度的增大更加明顯。

使用本文給出的動(dòng)剛度方法計(jì)算電主軸固有頻率，同時(shí)與不考慮動(dòng)剛度方法計(jì)算對(duì)比，圖5給出了電主軸的一階、二階臨界轉(zhuǎn)速計(jì)算結(jié)果對(duì)比，圖中虛線為同步激勵(lì)線，同步激勵(lì)線與各固有頻率曲線的交點(diǎn)即為電主軸的各階臨界轉(zhuǎn)速。

可以看出，考慮動(dòng)剛度計(jì)算得到的電主軸臨界轉(zhuǎn)速為22 500 r/min、31 500 r/min，而不考慮動(dòng)剛度計(jì)算得到的電主軸臨界轉(zhuǎn)速則為25 000 r/min、38 000 r/min，兩種計(jì)算方法得到的一階和二階臨界轉(zhuǎn)速差別分別為11.1%和20.6%，這是由于軸承剛度隨轉(zhuǎn)速增大而改變，引起臨界轉(zhuǎn)速值的改變。

3.2電主軸動(dòng)態(tài)響應(yīng)分析

CF取68．2，銑削深度ae取2 mm，每齒進(jìn)給量af取0.2 mm，銑刀直徑dq取32 mm，銑削寬度ap為25 mm，銑刀齒數(shù)Z為4。在最大轉(zhuǎn)速18 000 r/min時(shí)，銑削激振力為：q(t)=1935cos(1200t)，將電主軸劃分為10個(gè)節(jié)點(diǎn)，求電主軸各節(jié)點(diǎn)位移響應(yīng)。

圖6為考慮動(dòng)剛度和未考慮動(dòng)剛度時(shí)，電主軸各節(jié)點(diǎn)位移響應(yīng)最大值曲線，可看出最大位移在電主軸前端，考慮動(dòng)剛度和未考慮動(dòng)剛度時(shí)最大位移分別為19 μm和15 μm，未考慮動(dòng)剛度時(shí)電主軸位移響應(yīng)整體偏小。

3.3電主軸動(dòng)態(tài)特性實(shí)驗(yàn)

在自主開發(fā)的高速大功率電主軸實(shí)驗(yàn)臺(tái)上測(cè)試，實(shí)驗(yàn)臺(tái)包括電主軸、變頻器、水冷裝置、空氣壓縮機(jī)，以及溫度、振動(dòng)、位移、LMS數(shù)據(jù)采集等測(cè)試系統(tǒng)。圖7是電主軸性能測(cè)試實(shí)驗(yàn)臺(tái)，電主軸最高轉(zhuǎn)速18 000 r/min，功率45 kW，支承軸承為角接觸球軸承，預(yù)緊量為10 μm。

使用該實(shí)驗(yàn)臺(tái)進(jìn)行高速電主軸前端振動(dòng)位移測(cè)試，轉(zhuǎn)速測(cè)試范圍為2 000 r/min到18 000 r/min，加載轉(zhuǎn)矩20 N·m，圖8為電主軸前端位移實(shí)驗(yàn)結(jié)果和計(jì)算結(jié)果對(duì)比，隨轉(zhuǎn)速增大振動(dòng)位移逐漸增大，考慮動(dòng)剛度計(jì)算與測(cè)試結(jié)果較為接近，高速時(shí)有一定誤差，是由于電主軸內(nèi)部零件存在不平衡量導(dǎo)致的振動(dòng)，實(shí)驗(yàn)結(jié)果驗(yàn)證了文中方法的可靠性。

4　結(jié)語(yǔ)

建立了考慮預(yù)緊的球軸承剛度求解模型，基于Riccati傳遞矩陣法，給出了考慮軸承動(dòng)剛度的電主軸臨界轉(zhuǎn)速計(jì)算方法和動(dòng)特性計(jì)算方法，對(duì)某高速電主軸動(dòng)特性進(jìn)行了分析，得到如下結(jié)果：

(1)球軸承徑向剛度隨轉(zhuǎn)速增大而減小，軸向剛度隨轉(zhuǎn)速增大而增大，采用定位預(yù)緊后，軸向和徑向剛度與沒有預(yù)緊相比均增大，且隨轉(zhuǎn)速增大預(yù)緊軸向剛度增大更加明顯。

(2)考慮滾動(dòng)軸承的動(dòng)剛度計(jì)算轉(zhuǎn)子系統(tǒng)的臨界轉(zhuǎn)速，與定剛度值計(jì)算相比有較大差別，為準(zhǔn)確計(jì)算電主軸臨界轉(zhuǎn)速應(yīng)考慮動(dòng)剛度。

(3)實(shí)驗(yàn)驗(yàn)證了考慮軸承動(dòng)剛度可得到更為準(zhǔn)確的電主軸動(dòng)特性結(jié)果，為電主軸動(dòng)態(tài)特性的準(zhǔn)確分析提供了方法。

[1]Mohammed A Alfares, Abdallah A Elsharkawy. Effects of axial preloading of angular contact ball bearings on the dynamics of a grinding machine spindle system[J]. Journal of Materials Processing Technology, 2003,136: 48-59.

[2]Lin Chi-wei, Tu Jay F, Kamman Joe. An integrated thermo mechanical dynamic model to characterize motorized machine tool spindles during very high speed rotation[J]. International Journal of Machine Tool & Manufacture, 2003, 43: 1035-1050.

[3]Kang Y, Chang Y P. Integrated CAE strategies for the design of machine tool spindle-bearing systems[J]. Finite Elements in Analysis and Design,2001, 37: 485-511.

[4]趙志剛, 巫少龍. 基于ANSYS的機(jī)床電主軸動(dòng)態(tài)性能分析[J]. 制造技術(shù)與機(jī)床, 2010(8): 149-152.

[5]李松生，楊柳欣，王兵華.高速電主軸軸系轉(zhuǎn)子動(dòng)力學(xué)特性分析[J].軸承, 2002(2):15-17.

[6]李蓓智，吳喜如，楊建國(guó)．精密數(shù)控車床電主軸一軸承-轉(zhuǎn)子系統(tǒng)動(dòng)態(tài)性能計(jì)算分析[J]．機(jī)械設(shè)計(jì)與制造，2010(3): 173-175．

[7]邱海飛. 高速陶瓷電主軸動(dòng)力學(xué)特性研究[J]. 制造技術(shù)與機(jī)床, 2014(12): 101-107.

[8]Harris T A. Rolling bearing analysis[M]. New York: Wiley, 2009.

[9]王黎欽，崔立，鄭得志. 航空發(fā)動(dòng)機(jī)高速球軸承動(dòng)態(tài)特性分析[J]. 航空學(xué)報(bào), 2007, 28(6): 1461-1467.

[10]聞邦椿. 高等轉(zhuǎn)子動(dòng)力學(xué)-理論、技術(shù)與應(yīng)用[M]. 北京：機(jī)械工業(yè)出版社, 2000.

(編輯李靜)

(收修改稿日期：2016-05-11)

如果您想發(fā)表對(duì)本文的看法，請(qǐng)將文章編號(hào)填入讀者意見調(diào)查表中的相應(yīng)位置。

Dynamic characteristics analysis of high speed motorized spindle based on dynamic stiffness of rolling bearing with preload

CUI Li

(Faculty of Engineering, Shanghai Polytechnic University, Shanghai 201209, CHN)

For spindle bearings with position preload, bearing stiffness model is presented based on the interactions model between the ball and rings. Newton-Raphson method is used to solve the stiffness matrix. Based on Riccati transfer matrix method, calculation method of critical speed of motorized spindle considering the dynamic stiffness of bearing is presented. The dynamic equation of the motorized spindle considering dynamic stiffness of bearing under milling excitation is given, and the dynamic responses are solved by using Runge-Kutta method. Taking a high speed milling spindle as an example, results show that radial stiffness of the bearing decrease with increase of rotating speed, the bearing stiffness increases with the positioning preload, and the increase becomes obvious with the increases of speed. More accuracy of dynamic responses of the spindle system can be obtained by considering dynamic stiffness. Finally the dynamic characteristic experiment of high-speed electric spindle demonstrates the accuracy of the proposed model.

high speed motorized spindle; dynamic stiffness; ball bearing; dynamic characteristics; critical speed

TH133.3

A

10.19287/j.cnki.1005-2402.2016.07.014

崔立，男，1981年生，博士，副教授，研究方向?yàn)楦咚佥S承技術(shù)、摩擦學(xué)、轉(zhuǎn)子動(dòng)力學(xué)，已發(fā)表論文30余篇。

160723