數(shù)控機床滾動導(dǎo)軌結(jié)合面剛度的有限元分析

范晉偉，王鴻亮，張?zhí)m清，袁 帥，紀 實

（北京工業(yè)大學(xué) 機械工程與應(yīng)用電子技術(shù)學(xué)院，北京 100124）

數(shù)控機床滾動導(dǎo)軌結(jié)合面剛度的有限元分析

范晉偉，王鴻亮，張?zhí)m清，袁 帥，紀 實

（北京工業(yè)大學(xué) 機械工程與應(yīng)用電子技術(shù)學(xué)院，北京 100124）

直線滾動導(dǎo)軌是數(shù)控機床的重要零部件，其結(jié)合面的接觸剛度直接影響對數(shù)控機床整機動態(tài)特性分析。對兩種不同直線滾動導(dǎo)軌進行結(jié)構(gòu)簡化，運用Pro/E對簡化后的滾動導(dǎo)軌進行幾何建模，將幾何模型轉(zhuǎn)化為IGES中間格式導(dǎo)入ANSYS Workbench有限元分析軟件中進行靜力學(xué)分析。通過對滑塊的法向和切向分別施加不同大小的載荷得到變形量的平均值，運用剛度計算公式得到了簡化后滑塊表面的法向剛度和切向剛度。提出了應(yīng)用滑塊表面剛度來估算結(jié)合面接觸剛度，并通過反求法進行了驗證，最后證明了該方法的正確性。

滾動導(dǎo)軌；結(jié)合面；有限元；接觸剛度；彈簧剛度；等效

0　引言

機床及各類機械一般都是由零部件組成，而零部件之間相互接觸的表面稱為結(jié)合面（或接觸面）。結(jié)合面包括固定結(jié)合面、半固定結(jié)合面和運動結(jié)合面。無論是哪種結(jié)合面，均是“柔性結(jié)合”，即結(jié)合面既表現(xiàn)出彈性又有阻尼，既儲存能量又消耗能量[1]。許多研究表明，機床中結(jié)合面的剛度約占機床總剛度的60%～80%，結(jié)合面引起的變形量約占機床總體靜變形量的85%～90%[2]。結(jié)合面剛度是機床總剛度的重要組成部分，甚至是機床的薄弱環(huán)節(jié)[3]，結(jié)合面的變形不僅會影響組成結(jié)構(gòu)的幾何關(guān)系，甚至?xí)绊憴C床的加工精度及壽命。在使用有限元軟件對機床部件或整機進行動態(tài)特性分析時，結(jié)合面模型建立的正確與否直接影響分析的結(jié)果。

導(dǎo)軌是機床的重要組成部分，其中滑動導(dǎo)軌和滾動導(dǎo)軌是最常見的形式，而導(dǎo)軌與滑塊之間的結(jié)合面是常見的運動結(jié)合面。滾動導(dǎo)軌是在滑動導(dǎo)軌的基礎(chǔ)上發(fā)展起來的，滾動導(dǎo)軌相對于滑動導(dǎo)軌具有以下優(yōu)點[4]：1）定位精度高；2）精度保持性好；3）運動速度高；4）工作效率高。隨著機床加工精度和工作效率的提高，滾動導(dǎo)軌正在逐漸取代滑動導(dǎo)軌。考慮到滑動導(dǎo)軌結(jié)合面剛度的研究理論十分復(fù)雜且不夠成熟，所得到的公式通用性也不強，目前為止尚未建立完整或通用的結(jié)合面動態(tài)特性參數(shù)數(shù)據(jù)庫[5]。滑動導(dǎo)軌結(jié)合面剛度與接觸表面粗糙度、所受載荷及移動速度等有關(guān)[6]，具有強烈的非線性，求解十分復(fù)雜。故這里只對滾動導(dǎo)軌結(jié)合面的接觸剛度進行了研究。隨著科學(xué)技術(shù)的發(fā)展，有限元分析分析軟件的應(yīng)用越來越廣泛，而且精度越來越高。常見的有限元分析軟件有ANSYS、ABAQUS、MSC/NASTRAN等，一般都能進行線性和非線性的結(jié)構(gòu)動力、靜力學(xué)分析。論文主要采用ANSYS Workbench對滾動導(dǎo)軌法向和切向進行靜剛度分析，將測量的變形量采用Excel曲線擬合得到結(jié)合面法向剛度和切向剛度。

1　滾動導(dǎo)軌的類型及結(jié)合面等效動力學(xué)模型

根據(jù)滾動導(dǎo)軌的工作狀況的不同，常見的類型有三種，分別為[4]：徑向載荷型、四方向等載荷型及高載荷型，如圖1所示。徑向載荷型滾動導(dǎo)軌用于表面承受徑向載荷，適用于安裝在水平面上，接觸角為90°。四方向等載荷型滾動導(dǎo)軌四個方向具有相同承受載荷大小的能力，即四個方向結(jié)合面等效剛度相同（下述將對其進行驗證），能用于安裝在水平和垂直方向，適應(yīng)能力強，接觸角為45°。高載荷型滾動導(dǎo)軌用于對加工要求苛刻的機床，如施加載荷時伴隨很大的振動等。

圖1　常見滾動導(dǎo)軌類型

考慮到結(jié)合面“柔性結(jié)合”的特點，在對滾動導(dǎo)軌結(jié)合面的研究中，通常對滾動導(dǎo)軌結(jié)合面采用彈簧和阻尼器來等效。

等效后系統(tǒng)的動力學(xué)方程可統(tǒng)一表示為[7]：

式中：Ma、Mb是滑塊和導(dǎo)軌的質(zhì)量矩陣，Ca、Cb是滑塊和導(dǎo)軌阻尼矩陣，Ka、Kb是滑塊和導(dǎo)軌的剛度矩陣，Cj、Kj分別是結(jié)合間的阻尼矩陣和剛度矩陣，xa、xb是滑塊和導(dǎo)軌的坐標向量。其中Cj、Kj又可以分別用結(jié)合面間的法向阻尼、切向阻尼以及法向剛度、切向剛度表示為：

2　滾動導(dǎo)軌的幾何建模和有限元建模

由圖1和圖2可知徑向載荷型滾動導(dǎo)軌動力學(xué)模型可由四方向等載荷型滾動導(dǎo)軌與高載荷型滾動導(dǎo)軌模型聯(lián)合等效，故本文只對四方向等載荷型滾動導(dǎo)軌及高載荷型滾動導(dǎo)軌進行幾何建模和運用ANSYS Workbench進行有限元分析。分析結(jié)合面的剛度問題就是求解在載荷的作用下發(fā)生的彈性變形量，考慮到直線滾動導(dǎo)軌滾珠的直徑大小相同，為減少計算的重復(fù)性，可對導(dǎo)軌、滑塊和滾珠進行簡化[8]。下面分別以某高載荷型滾動導(dǎo)軌和四方向等載荷型滾動導(dǎo)軌為例進行剛度分析。運用Pro/E三維軟件進行幾何建模，為方便有限元分析，應(yīng)簡化幾何模型，去掉小圓孔及倒角等對分析結(jié)果影響不大的結(jié)構(gòu)。簡化后的幾何模型如圖2所示。

圖2　簡化后滾動導(dǎo)軌幾何模型

將幾何模型轉(zhuǎn)化為IGES中間格式導(dǎo)入ANSYS Workbench有限元分析軟件中，使用Static Structural模塊進行剛度分析。首先定義材料類型，由資料可知，滑塊和導(dǎo)軌材料都是HT300，楊氏模量為1.43×1011N/m2，泊松比0.27，滾珠材料是0Gr18Ni12Mo2Ti，楊氏模量為1.93×1011N/m2，泊松比0.27。滑塊與滾珠及導(dǎo)軌與滾珠之間的接觸面類型系統(tǒng)默認設(shè)置為Bonded，而為接近真實情況這里將接觸類型設(shè)置為No Separation。在有限元分析中，網(wǎng)格劃分的質(zhì)量好壞直接影響分析的結(jié)果，由于ANSYS Workbench具有強大的網(wǎng)格劃分功能，在結(jié)合有限元網(wǎng)格劃分的基本原則[9]的前提下，采用自動網(wǎng)格劃分與接觸面網(wǎng)格劃分同時使用。最后得到的有限元模型如圖3所示。導(dǎo)軌一般固定安裝在床身或立柱上，故邊界條件應(yīng)定義為導(dǎo)軌底面固定約束。

圖3　滾動導(dǎo)軌有限元模型

3　兩種不同類型滾動導(dǎo)軌的剛度分析

3.1高載荷型滾動導(dǎo)軌的剛度分析

對高載荷型滾動導(dǎo)軌進行剛度分析，由于要求出結(jié)合面法向剛度和切向剛度，故需要在滑塊法向和切向分別施加載荷，為了更接近實際受載情況這里所施加的載荷均為面力。為減小誤差的影響，采用不同大小的載荷，在法向和切向分別施加100N、200N、300N、400N、500N的壓力。為更好地反映出法向接觸面和切向接觸面的變形量，在滑塊的上平面和左平面上均勻的取15的測點，測點的分布情況如圖4所示。通過力與滑塊表面變形量的關(guān)系來估算接觸面剛度，各測點的變形量通過坐標探針測量，測量結(jié)果如表1和表2所示。

圖4　高載荷型滾動導(dǎo)軌測點分布

表1　滑塊法向測點的變形量

表2　滑塊切向測點的變形量

滑塊在100N、200N、300N、400N、500N的壓力作用下法向和切向的平均變向量為將壓力與法向及切向平均變形量之間的關(guān)系通過Excel曲線擬合如圖5所示。

圖5　載荷與變形量的關(guān)系曲線

由剛度計算公式：

可知，擬合曲線的斜率即為剛度，故簡化后高載荷型滾動導(dǎo)軌的法向剛度Kn=184191N/mm、切向剛度Kt=72218N/mm。

在結(jié)合面的基礎(chǔ)參數(shù)中，剛度主要影響結(jié)構(gòu)系統(tǒng)的固有頻率和振動幅值，阻尼一般不影響系統(tǒng)固有頻率而影響振動幅值[10]。考慮到在靜力學(xué)分析中阻尼的存在并不影響分析的結(jié)果，為驗證所求剛度的正確性，可在導(dǎo)軌與滑塊結(jié)合面之間使用彈簧來等效滾珠，并重新進行有限元分析查看法向及切向的變形量。等效彈簧的布置如圖6所示，其中法向兩根彈簧剛度為92095.5N/mm，切向兩根彈簧剛度為36109N/mm。同樣，在法向和切向分別施加100N、200N、300N、400N、500N的壓力，并在相同位置布置測點探測變形量。有限元分析的過程與上述相同，分析結(jié)果如表3和表4所示。

圖6　高載荷型滾動導(dǎo)軌等效彈簧分布

表3　等效后滑塊測點的法向變形量

表4　等效后滑塊測點的切向變形量

同樣，采用Excel曲線擬合可得法向剛度Kn=171829N/mm和切向剛度Kt=68006N/mm，與先前所求滑塊法向剛度和切向剛度相比誤差分別為6.71%和5.83%，均在誤差范圍內(nèi)，故簡化后該直線滾動導(dǎo)軌結(jié)合面法向接觸剛度約為184191N/mm，切向接觸剛度約為72218N/mm。

3.2四方向等載荷型滾動導(dǎo)軌的剛度分析

四方向等載荷型滾動導(dǎo)軌結(jié)合面接觸剛度的分析與高載荷型滾動導(dǎo)軌剛度分析的過程相似，同樣運用力和滑塊表面變形量的關(guān)系來估算接觸面剛度。分別施加100N、200N、300N、400N、500N的壓力于滑塊的法向和切向，在滑塊的上平面和左平面分別均勻布置15測量點。在不同大小力的作用下可以求出法向平均變形量分別為7.2631×10-4mm、1.4526×10-3mm、2.1789×10-3mm、2.9052×10-3mm、3.6315×10-3mm，切向平均變形量分別為1.5209×10-4mm、3.0418×10-3mm、4.5627×10-3mm、6.0836×10-3mm、7.6171×10-3mm。此時考慮到力和幾何關(guān)系，計算接觸面剛度時力的大小為故法向剛度通過剛度計算公式并結(jié)合Excel曲線擬合可得導(dǎo)軌法向結(jié)合面接觸剛度為194718N/mm。由于該類型導(dǎo)軌的四方向剛度相等，故接觸面切向剛度也為194718N/mm。

運用彈簧剛度來等效滾珠剛度并反向驗證，等效彈簧的布置如圖7所示，6根彈簧剛度均為97359N/mm。在滑塊的法向和切向分別施加100N、200N、300N、400N、500N的壓力，并在同樣的位置布置測點。通過有限元分析的結(jié)果可得導(dǎo)軌法向平均變形量分別為7.5235×10-4mm、1.5047×10-3mm、2.2571×10-3mm、3.0094×10-3mm、3.7618×10-3mm，切向平均變形量分別為1.5903×10-3mm、3.1807×10-3mm、4.7711×10-3mm、6.3614×10-3mm、7.9518×10-3mm。通過反求法得到的滾動導(dǎo)軌法向剛度為187973N/mm，與先前所求結(jié)合面法向剛度相比誤差為3.46%，應(yīng)用彈簧剛度等效滾珠后滑塊切向變形量與等效前相比誤差約為4.35%，均在誤差范圍內(nèi)，故簡化后四方向等載荷型滾動導(dǎo)軌結(jié)合面的法向接觸剛度和切向接觸剛度約為194718N/mm。

圖7　四方相等載荷型滾動導(dǎo)軌等效彈簧分布

4　結(jié)論

機床整機的動態(tài)特性直接影響機床的加工精度、可靠性及壽命等，在進行機床整機動態(tài)特性分析時結(jié)合面之間的接觸剛度模型建立的正確與否十分關(guān)鍵。結(jié)合面剛度是機床整機剛度的重要組成部分，尤其是導(dǎo)軌結(jié)合面剛度對工件的加工精度影響很大。對滾動導(dǎo)軌建立了等效動力學(xué)模型，在動力學(xué)模型中結(jié)合面之間的接觸剛度可以通過接觸面的法向剛度和切向剛度來等效。對高載荷型滾動導(dǎo)軌和四方向等載荷型滾動導(dǎo)軌進行了結(jié)構(gòu)簡化，應(yīng)用Pro/E軟件進行了幾何建模。將幾何模型轉(zhuǎn)化為中間格式導(dǎo)入ANSYS Workbench有限元分析軟件中進行靜力學(xué)分析，通過有限元分析求出了導(dǎo)軌滑塊的法向剛度和切向剛度。提出了應(yīng)用導(dǎo)軌滑塊的法向剛度和切向剛度來估算導(dǎo)軌結(jié)合面之間接觸面剛度的理論，并通過反向求解進行了驗證，最后證明了所提出的理論是正確的。該理論不僅能用于求解滾動導(dǎo)軌結(jié)合面接觸剛度，也為求解滾動導(dǎo)軌及機床整機動態(tài)特性提供了必要條件。

[1] 廖伯瑜,周新民,尹志宏.現(xiàn)代機械動力學(xué)及其工程應(yīng)用[M].北京:機械工業(yè)出版社,2003.

[2] 溫淑華.結(jié)合面接觸特性理論建模及仿真[M].北京:國防工業(yè)出版社,2012.

[3] 趙宏林,丁慶新,曾鳴,等.機床結(jié)合部特性的理論解析及應(yīng)用[J].機械工程學(xué)報,2008,(12):208-214.

[4] 苗青.數(shù)控機床導(dǎo)軌結(jié)合面的接觸剛度分析[D].洛陽:河南科技大學(xué),2014.

[5] 王維友.機床滑動導(dǎo)軌結(jié)合面動態(tài)特性參數(shù)識別實驗研究[D].南京:南京理工大學(xué),2011.

[6] 溫淑華,張學(xué)良,武美先,等.結(jié)合面法向接觸剛度分型模型建立與仿真[J].農(nóng)業(yè)機械學(xué)報,2009,40(11):197-202.

[7] 朱爽.超高速磨床結(jié)構(gòu)建模的研究[J].沈陽工程學(xué)院學(xué)報（自然科學(xué)版）,2014,10(3):275-279.

[8] 劉稱意.考慮導(dǎo)軌結(jié)合面影響的數(shù)控機床動態(tài)特性分析及結(jié)構(gòu)優(yōu)化[D].沈陽:東北大學(xué),2010.

[9] 杜平安.有限元網(wǎng)格劃分的基本原則[J].機械設(shè)計與制造,2000,(01):34-36.

[10] 崔中,文桂林,陳桂平,等.高速磨床整機動態(tài)特性研究[J].中國機械工程,2010,21(7):782-787.

Finite element analysis of NC machine rolling guide joint stiffness

FAN Jin-wei, WANG Hong-liang, ZHANG Lan-qing, YUAN Shuai, JI Shi

TH114；TG659

A

1009-0134(2016)02-0134-05

2015-11-29

國家科技重大專項：中高檔國產(chǎn)數(shù)控磨床可靠性規(guī)模化提升工程（2013ZX04011013）；國家自然科學(xué)基金項目：數(shù)控機床加工缺陷成因的逆向追蹤理論與實驗研究（51275014）

范晉偉（1965 -），男，河南西平人，教授，博士，研究方向為數(shù)控磨床可靠性與超精密加工。