0 前言
結(jié)合部的動(dòng)態(tài)特性參數(shù)是影響機(jī)床整機(jī)動(dòng)力學(xué)性能的關(guān)鍵因素。根據(jù)統(tǒng)計(jì)，機(jī)床中出現(xiàn)的振動(dòng)問(wèn)題有60%以上是源自結(jié)合部[1]，因此，結(jié)合部動(dòng)態(tài)特性參數(shù)分析一直是國(guó)內(nèi)外學(xué)術(shù)界研究熱點(diǎn)。直線(xiàn)導(dǎo)軌已廣泛地應(yīng)用于現(xiàn)代數(shù)控機(jī)床，圍繞其動(dòng)態(tài)參數(shù)建模，國(guó)內(nèi)外學(xué)者已做了大量的有益工作，其研究方法概括起來(lái)有三種:
① 試驗(yàn)測(cè)試法。該方法主要通過(guò)錘擊試驗(yàn)分別測(cè)定導(dǎo)軌豎直方向和水平方向的傳遞函數(shù)，再根據(jù)模態(tài)理論識(shí)別出兩個(gè)方向的剛度、阻尼系數(shù)，從而獲得直線(xiàn)導(dǎo)軌結(jié)合部的動(dòng)力學(xué)參數(shù)[2-3]。也有學(xué)者將直線(xiàn)導(dǎo)軌整體結(jié)構(gòu)離散為若干個(gè)子結(jié)構(gòu)或子系統(tǒng)，通過(guò)試驗(yàn)測(cè)試的方法得到各個(gè)子系統(tǒng)的剛度矩陣和阻尼矩陣，再將各個(gè)子結(jié)構(gòu)的運(yùn)動(dòng)方程綜合起來(lái)，從而得到整個(gè)系統(tǒng)的運(yùn)動(dòng)方程和動(dòng)力學(xué)特性[3-5]。② 理論計(jì)算法。該方法主要以結(jié)合部的基本特性參數(shù)為基礎(chǔ)，簡(jiǎn)化結(jié)合部模型，求得結(jié)合部各處的剛度和阻尼，或者通過(guò)有限元分析方法，得到導(dǎo)軌結(jié)合部的特性參數(shù)[1，6-9]。③ 試驗(yàn)測(cè)試與理論計(jì)算相結(jié)合方法。該方法首先通過(guò)試驗(yàn)方法測(cè)試導(dǎo)軌結(jié)合部的動(dòng)力學(xué)參數(shù)，然后再用有限元方法求解，為機(jī)床的動(dòng)態(tài)設(shè)計(jì)提供一種有效的方法[10-13]。然而，上述研究工作存在兩方面問(wèn)題：① 傳統(tǒng)直線(xiàn)導(dǎo)軌結(jié)合部模型，全部忽略了作為直線(xiàn)導(dǎo)軌必要部件——滾珠絲杠副對(duì)結(jié)合部動(dòng)態(tài)特性參數(shù)的影響，因此，所建立的模型不能很好地描述直線(xiàn)導(dǎo)軌副的動(dòng)態(tài)特性；② 試驗(yàn)測(cè)試法需要搭建較為復(fù)雜試驗(yàn)裝備，而有限元法建模較為繁瑣，且依賴(lài)于大型商用軟件。
本文以帶滾珠絲杠副直線(xiàn)導(dǎo)軌結(jié)合部為研究對(duì)象，借助于彈性力學(xué)中的赫茲接觸理論，分析計(jì)算直線(xiàn)滾動(dòng)導(dǎo)軌的線(xiàn)剛度、滾珠絲杠副和角接觸滾動(dòng)軸承的軸向剛度，建立帶滾珠絲杠副的機(jī)床直線(xiàn)導(dǎo)軌結(jié)合部的動(dòng)力學(xué)模型。通過(guò)一款帶滾珠絲杠副的直線(xiàn)導(dǎo)軌結(jié)合部的立式加工中心的整機(jī)動(dòng)態(tài)特性分析，驗(yàn)證了提出的結(jié)合部動(dòng)態(tài)特性模型的有效性和可行性。
1 結(jié)合部動(dòng)力學(xué)建模
圖1 為數(shù)控機(jī)床典型的直線(xiàn)滾動(dòng)導(dǎo)軌進(jìn)給系統(tǒng)的結(jié)構(gòu)示意圖，通過(guò)滾珠絲杠副和直線(xiàn)導(dǎo)軌副的組合來(lái)實(shí)現(xiàn)工作臺(tái)x 方向的進(jìn)給運(yùn)動(dòng)。
圖2 為帶滾珠絲杠副的直線(xiàn)滾動(dòng)導(dǎo)軌結(jié)合部的動(dòng)力學(xué)模型，結(jié)合部的x 向剛度為滾珠絲杠副的軸向剛度，y 向剛度為直線(xiàn)滾動(dòng)導(dǎo)軌的橫向剛度， z 向剛度為直線(xiàn)滾動(dòng)導(dǎo)軌的垂向剛度，并且每個(gè)方向上的剛度值均由四組彈簧等價(jià)模擬。因此，只須計(jì)算出滾珠絲杠副的軸向剛度、直線(xiàn)滾動(dòng)導(dǎo)軌的橫向剛度和垂向剛度，便可確定結(jié)合部的動(dòng)態(tài)特性參數(shù)。
現(xiàn)假設(shè)如下所述。
(1) 忽略表面粗糙度對(duì)結(jié)合部剛度特性的影響。因?yàn)閹L珠絲杠直線(xiàn)導(dǎo)軌結(jié)合部的接觸副均為高副，在預(yù)緊載荷或工作載荷的作用下，其接觸區(qū)彈性變形量可達(dá)到數(shù)個(gè)微米(一般在2～7 μm 之間)，遠(yuǎn)遠(yuǎn)大于滾道或滾珠的表面粗糙度(一般小于0.08μm)[14-15]。因此，在計(jì)算接觸剛度時(shí)，可不考慮表面粗糙度的影響，進(jìn)而可忽略材料塑性變形的影響。
(2) 滾珠與滾道接觸時(shí)只產(chǎn)生彈性變形，并服從Hooke 定理。在正常情況下，帶滾珠絲杠直線(xiàn)導(dǎo)軌結(jié)合部各接觸副變形均在材料彈性范圍內(nèi)，接觸點(diǎn)處產(chǎn)生的塑性變形量不超過(guò)滾動(dòng)體直徑的萬(wàn)分之一(遠(yuǎn)小于接觸副彈性變形量)[16]。因此，滾珠與滾道接觸只考慮彈性變形的假設(shè)，在分析帶滾珠絲杠副直線(xiàn)導(dǎo)軌結(jié)合部剛度時(shí)是適宜的。
1.1 滾珠絲杠副軸向剛度的計(jì)算
滾珠絲杠副的軸向剛度kx 為與滾珠絲杠副相關(guān)聯(lián)的零部件剛度的串聯(lián)總和[17]，其中以絲杠、螺母組件和支承軸承的軸向剛度影響最大，其他零部件的影響較小，可忽略不計(jì)，其動(dòng)力學(xué)模型如圖3所示。
顯然，進(jìn)給絲杠傳動(dòng)系統(tǒng)的軸向剛度kx 可表示為
式中 kx ——滾珠絲杠副軸向剛度
kS ——絲杠軸向剛度
kN ——螺母組件軸向剛度
kB ——支撐軸承軸向剛度
1.1.1 絲杠的軸向剛度
滾珠絲杠軸向剛度的計(jì)算，隨滾珠絲杠的支承方式不同而存在著差異。限于篇幅，本文只討論兩端固定支承結(jié)構(gòu)方式(圖4)，這也是滾珠絲杠副最為常見(jiàn)的一種支承形式，并不失一般性。
可由材料力學(xué)得到其軸向剛度為[18]
式中 d ——絲杠的螺紋底徑
l ——載荷作用點(diǎn)至右端軸承的距離
E ——絲杠材料的縱向彈性模量
當(dāng)螺母處于絲杠中間位置時(shí)，絲杠軸向剛度kS的值為最小，即
本文以絲杠的最小剛度kSmin 作為其軸向剛度值。
1.1.2 螺母組件軸向剛度的計(jì)算
螺母組件軸向剛度可由作用在螺母上的軸向載荷除以其軸向變形量求得，即
假設(shè)如下：① 負(fù)荷垂直于接觸表面，也就是說(shuō)，接觸表面完全光滑，不計(jì)及滾珠與滾道面間的摩擦力；② 滾珠與滾道接觸面的尺寸與其曲率半徑相比是很小的；③ 考慮到滾珠絲杠副的工作轉(zhuǎn)速一般都較低，在分析過(guò)程中不考慮滾珠離心力和陀螺力矩的影響；④ 不考慮由于制造誤差所產(chǎn)生的影響，即軸向工作載荷均勻地分配給每個(gè)滾珠。在以上假設(shè)下，螺母組件的軸向變形量完全由滾珠與滾道面的彈性接觸變形所引起，可利用赫茲彈性接觸理論來(lái)進(jìn)行理論計(jì)算。其計(jì)算過(guò)程主要可分為以下幾個(gè)步驟。
滾珠與螺母或絲杠滾道面的法向接觸變形量的計(jì)算。根據(jù)赫茲接觸理論，兩彈性體由于彈性變形引起的相對(duì)位移量(彈性接近量)為
滾珠與螺母滾道接觸點(diǎn)處的4 個(gè)主曲率分別為
當(dāng)已知τ值便可通過(guò)查表得到J 和ma 的值[20]。至此，只要將J、ma 以及Σρ 和p 的值代入式(5)，即可得到δ 的值。
(2) 單個(gè)滾珠軸向變形量的計(jì)算。單個(gè)滾珠在法向接觸壓力p 作用下而產(chǎn)生彈性變形如圖5 所示，由滾珠法向彈性接觸變形所產(chǎn)生的法向彈性位移量
δp 為
δ p =δ1 +δ 2 (15)
(3) 單個(gè)滾珠法向接觸力的計(jì)算。雙螺母預(yù)緊結(jié)構(gòu)(圖6)是滾珠絲杠最為常見(jiàn)的預(yù)緊方式，其螺母組件軸向剛度的計(jì)算方法如下。
假定滾珠絲杠副中螺母A上每個(gè)滾珠給予絲杠滾道面的法向作用力為pA，螺母B 上每個(gè)滾珠給予絲杠滾道面的法向作用力為pB，且預(yù)緊墊片通過(guò)螺母A、B 給予絲杠的預(yù)緊法向力為pP。則由絲杠靜力平衡可得
Fx − pA z sinα cosφ + pBz sinα cosφ = 0 (17)
在軸向載荷Fx 作用下，螺母A 相對(duì)絲杠所產(chǎn)生的軸向彈性接觸變形量δA 應(yīng)恰好等于螺母B 相對(duì)于絲杠所產(chǎn)生的軸向彈性恢復(fù)量。
由力學(xué)疊加原理可知，作用于同一物體的合力所產(chǎn)生的變形，應(yīng)等于各分力在同方向上所產(chǎn)生的變形量之和；又由赫茲接觸理論可知，兩彈性體的彈性趨近量與其法向壓力的2/3 次方成正比的關(guān)系，于是有
當(dāng)已知Fx 和pP 時(shí)，可由式(17)和式(19)聯(lián)合求解得到pA 和pB 的值。
得到了單個(gè)滾珠的法向接觸力后，便可由式 (5)～(16)求得螺母A 相對(duì)于絲杠所產(chǎn)生的軸向變形量δA，也即軸向載荷Fx 沿軸向的位移量δN，并最終由式(4)獲得螺母組件的軸向剛度kN。
1.1.3 支承軸承的軸向剛度計(jì)算
試驗(yàn)表明，軸承的變形量占整個(gè)傳動(dòng)系統(tǒng)總變形量的50%以上[20]。本文以數(shù)控機(jī)床常用的一對(duì)角接觸球軸承，背對(duì)背安裝的結(jié)構(gòu)為例(圖7)，進(jìn)行軸承剛度計(jì)算。之前螺母組件軸向剛度計(jì)算的假設(shè)條件，仍然適用于軸承軸向剛度的計(jì)算，在此不再贅述。與螺母組件軸向剛度的計(jì)算過(guò)程類(lèi)似，軸承軸向剛度的計(jì)算也可以分為三步。
(1) 利用式(5)～(14)計(jì)算滾珠與內(nèi)圈或外圈滾道面的法向接觸變形量。
(2) 利用式(15)和式(16)計(jì)算單個(gè)滾珠的軸向變形量。
(3) 單個(gè)滾珠法向接觸力的計(jì)算。
假定各軸承中單個(gè)滾珠所受的法向接觸力分別為p1、p2、p3 和p4，顯然p1= p3，p2= p4。由于軸承預(yù)緊而使其單個(gè)滾珠承受的法向力為p0。軸承1、2 中滾珠的受力情況與雙螺母預(yù)緊結(jié)構(gòu)中滾珠的受力情況類(lèi)似(圖6)，而軸承3、4 的受力情況與軸承1、2 的受力情況完全一樣，所以由靜力平衡條件可得
2 p1nsinβ − p2nsinβ = Fx (20)
式中 n ——單個(gè)軸承中的滾珠數(shù)
β ——軸承中滾珠與滾道面的接觸角
同樣由力學(xué)疊加原理和赫茲接觸理論可得到
當(dāng)已知n、β、Fx 和p0 時(shí)，可由式(20)和式(21)求得p1 和p2 的值。
得到了單個(gè)滾珠的法向接觸力后，由式(5)～ (16)可求得軸承外圈相對(duì)于軸承內(nèi)圈的軸向位移量δB，最后由式(4)求得軸承的軸向剛度kB。
在運(yùn)用式(5)～(16)時(shí)，應(yīng)將絲杠的結(jié)構(gòu)參數(shù)用軸承內(nèi)圈的結(jié)構(gòu)參數(shù)來(lái)替代，而將螺母的結(jié)構(gòu)參數(shù)用軸承外圈的結(jié)構(gòu)參數(shù)來(lái)替代，由于不涉及到滾珠絲杠副中螺紋升角φ 這一概念，須將公式中的φ 值設(shè)為零值。
通過(guò)以上的分析，分別得到了滾珠絲杠副中絲杠的軸向剛度kS，螺母組件的軸向剛度kN 和支承軸承的軸向剛度kB，可方便地由式(1)求得整個(gè)滾珠絲杠副系統(tǒng)的軸向剛度kx。
1.2 直線(xiàn)滾動(dòng)導(dǎo)軌線(xiàn)剛度的計(jì)算
一般數(shù)控機(jī)床上的直線(xiàn)滾動(dòng)導(dǎo)軌副都是由兩滑軌和四滑塊組成，因此只須計(jì)算單個(gè)滑塊的線(xiàn)剛度便可得到整個(gè)直線(xiàn)滾動(dòng)導(dǎo)軌副的剛度。圖8為單個(gè)直線(xiàn)滾動(dòng)導(dǎo)軌的結(jié)構(gòu)示意圖及動(dòng)力學(xué)模型圖。
以下分別計(jì)算單個(gè)直線(xiàn)滾動(dòng)導(dǎo)軌的橫向(y 向)剛度和垂向(z 向)剛度。
1.2.1 垂向剛度的計(jì)算
如圖9 所示，當(dāng)垂向力Fz 作用在滑塊上時(shí)，各列中單個(gè)滾珠的彈性力分別為F1、F2、F3 和F4，其中F1= F2 ，F3 =F4。γ 為滾珠與滾道面之間的接觸角。
由靜力平衡條件可以得到
2 (F1 − F3) msinγ = Fz (22)
由力學(xué)疊加原理和赫茲接觸理論可得到
式中 m ——單列滾道的接觸滾珠數(shù)
F0 ——由預(yù)壓載荷引起的單個(gè)滾珠的法向力
當(dāng)已知m、γ 、Fz 和F0 時(shí)，可由式(22)、(23)求得F1 和F3 的值。已知單個(gè)滾珠所受法向力，參照?qǐng)D5 所示的計(jì)算分析方法，便可求得滾珠的變形量，最后得到單個(gè)直線(xiàn)滾動(dòng)軸承的垂向剛度kg，整個(gè)直線(xiàn)滾動(dòng)導(dǎo)軌副由四個(gè)同樣的導(dǎo)軌單元并聯(lián)組成，因此總的垂向剛度kz＝4kg。
值得關(guān)注的是，由于直線(xiàn)滾動(dòng)導(dǎo)軌的結(jié)構(gòu)與滾珠絲杠、滾動(dòng)軸承不同，滾珠與滾道面接觸處的四個(gè)主曲率中，ρ22 的值應(yīng)取零。
1.2.2 橫向剛度的計(jì)算
直線(xiàn)滾動(dòng)導(dǎo)軌橫向(y 向)剛度的計(jì)算方法與垂向剛度的計(jì)算方法相同，不再贅述。
通過(guò)以上分析計(jì)算，得到帶滾珠絲杠副的直線(xiàn)滾動(dòng)導(dǎo)軌結(jié)合部空間三個(gè)相互垂直方向上的靜剛度，它們分別是滾珠絲杠副的軸向剛度kx，直線(xiàn)滾動(dòng)導(dǎo)軌副的橫向剛度ky 和垂向剛度kz。這樣，便建立帶滾珠絲杠副導(dǎo)軌結(jié)合部的動(dòng)態(tài)特性參數(shù)模型。
2 算例
2.1 基本參數(shù)
本文以一款立式加工中心為例，將本文建立的帶滾珠絲杠副直線(xiàn)滾動(dòng)導(dǎo)軌模型，應(yīng)用于加工中心動(dòng)力學(xué)特性分析，驗(yàn)證理論模型的有效性。表1～4列出了帶滾珠絲杠的直線(xiàn)滾動(dòng)導(dǎo)軌副的幾何、物理參數(shù)。
2.2 結(jié)果與討論
2.2.1 工作載荷對(duì)剛度的影響
為揭示工作載荷對(duì)結(jié)合部剛度特性參數(shù)的影響規(guī)律，作者運(yùn)用本文提出的計(jì)算方法，分析了結(jié)合部三個(gè)方向剛度與工作載荷的變化規(guī)律。圖10所示為結(jié)合部x 向剛度與其承受載荷的關(guān)系曲線(xiàn)。
查閱手冊(cè)，螺母組件的額定動(dòng)載荷pc=52.7 kN，按照技術(shù)要求，工作載荷F 在0.2pc～0.3pc 范圍內(nèi)變化?？梢钥闯?，結(jié)合部x 向剛度值隨著工作載荷的增大而略有上升，上升幅度在7%左右。
由于直線(xiàn)滾動(dòng)導(dǎo)軌滾珠與滾道面的接觸角為45°，所以結(jié)合部z 向接觸剛度與工作載荷的關(guān)系與y 向相等。
以上分析揭示了工作載荷對(duì)結(jié)合部剛度特性參數(shù)的影響規(guī)律，在工作載荷允許的變化范圍內(nèi)，剛度值的變化量一般為10%左右。
2.2.2 結(jié)合部動(dòng)態(tài)參數(shù)應(yīng)用
本文研究的立式加工中心為x、y、z 三軸聯(lián)動(dòng)機(jī)床，帶滾珠絲杠副的直線(xiàn)滾動(dòng)導(dǎo)軌結(jié)合部，分別位于立柱和銑頭、床身和床鞍、床鞍和工作臺(tái)之間 (圖12)。用彈簧單元模擬結(jié)合部接觸剛度，以靜止?fàn)顟B(tài)(假定零工作載荷)為對(duì)象，經(jīng)計(jì)算，可以得到此結(jié)合部x、y、z 三個(gè)方向上的剛度值分別是：
kx=0.419 7 MN/mm，ky=3.512 MN/mm，kz=3.512MN/mm。基于ANSYS 軟件，分析機(jī)床整機(jī)固有頻率特性，得到前四階固有頻率分別為85.75 Hz， 102.15 Hz，137.33 Hz，145.44 Hz，其振型如圖13所示。
表5 列出了結(jié)合部按照剛性連接簡(jiǎn)化處理與運(yùn)用本文理論模型的整機(jī)固有頻率特性比較。
可以看出，結(jié)合部動(dòng)態(tài)參數(shù)的描述是影響機(jī)床整機(jī)動(dòng)力學(xué)性能的重要因素，特別是對(duì)于高階固有頻率特性計(jì)算，影響更為明顯，因此，研究帶滾珠絲杠副直線(xiàn)滾動(dòng)導(dǎo)軌接觸剛度動(dòng)態(tài)模型，對(duì)機(jī)床整機(jī)動(dòng)態(tài)特性分析與設(shè)計(jì)有重要意義。
3 結(jié)論
(1) 運(yùn)用赫茲接觸理論，分析計(jì)算直線(xiàn)滾動(dòng)導(dǎo)軌的線(xiàn)剛度、滾珠絲杠副和角接觸球軸承的軸向剛度，在此基礎(chǔ)上，建立了帶滾珠絲杠副機(jī)床直線(xiàn)導(dǎo)軌結(jié)合部的動(dòng)態(tài)剛度特性模型，開(kāi)發(fā)了通用計(jì)算軟件。
(2) 以一臺(tái)立式加工中心為對(duì)象，分析計(jì)算了結(jié)合部三個(gè)方向動(dòng)態(tài)剛度隨工作載荷的變化規(guī)律，揭示了在工作載荷變化范圍內(nèi)，結(jié)合部動(dòng)態(tài)剛度的變化值一般小于10%。
(3) 考慮帶滾珠絲杠副的直線(xiàn)導(dǎo)軌結(jié)合部動(dòng)態(tài)影響，分析計(jì)算了立式加工中心的整機(jī)動(dòng)態(tài)特性，說(shuō)明結(jié)合部動(dòng)態(tài)參數(shù)是影響機(jī)床整機(jī)動(dòng)力學(xué)性能的重要因素，特別是對(duì)于高階固有頻率影響更為明顯。研究表明：帶滾珠絲杠副的直線(xiàn)導(dǎo)軌結(jié)合部動(dòng)態(tài)參數(shù)分析模型，具有物理概念清晰、計(jì)算簡(jiǎn)單等特點(diǎn)。
如果您有機(jī)床行業(yè)、企業(yè)相關(guān)新聞稿件發(fā)表，或進(jìn)行資訊合作，歡迎聯(lián)系本網(wǎng)編輯部， 郵箱：skjcsc@vip.sina.com
