數(shù)控落地銑鏜床滑枕變形有限元分析及補償

余劍武 ++桂林 ++李民選 ++羅嗣春 ++張漢衛(wèi)

摘要：TK6920B型大型數(shù)控落地鏜銑床由于滑枕自質(zhì)量和主軸箱的重心偏移而導致滑枕變形，從而影響機床的加工精度。針對這一問題，提出使用撓曲線預(yù)應(yīng)力加工法和液壓拉桿法對滑枕變形進行補償，為了減少實驗次數(shù)和降低成本，本文使用Ansys WB有限元軟件對滑枕變形分析和變形補償效果進行預(yù)測。對滑枕變形量的仿真分析結(jié)果和實驗測量結(jié)果進行了對比分析，結(jié)果表明，滑枕變形量的誤差在5 μm以下，有限元分析結(jié)果是可靠的；補償后的滑枕最大變形量為20 μm左右，說明了補償方法的補償效果可以滿足實際生產(chǎn)中加工精度的要求。

關(guān)鍵詞：數(shù)控落地銑鏜床；滑枕變形；預(yù)應(yīng)力撓曲加工法；液壓補償法；有限元分析

中圖分類號：TG58 文獻標識碼：A

數(shù)控落地銑鏜床在加工過程中隨著滑枕伸出主軸箱，滑枕前端會產(chǎn)生偏離理想直線的誤差。產(chǎn)生這種誤差的主要原因有＼[1＼]：1）滑枕在主軸箱內(nèi)移動使滑枕和主軸箱整體的重心偏移，在垂直面內(nèi)產(chǎn)生向下的微小偏移，在水平方向產(chǎn)生前移，從而使主軸箱產(chǎn)生微小傾斜；2）數(shù)控落地銑鏜床在工作過程中滑枕伸出主軸箱的距離最多可以達到1 200 mm，由于滑枕、銑軸、鏜軸及相關(guān)附件自重的影響產(chǎn)生撓曲變形，而且隨著滑枕伸出主軸箱的距離增大其變形量也會相應(yīng)增大。數(shù)控落地銑鏜床滑枕前端的這種變形誤差不僅對零件的加工精度產(chǎn)生影響，還會影響機床的使用壽命，所以落地銑鏜床的滑枕變形成為亟待解決的關(guān)鍵技術(shù)難題。

湖南大學學報（自然科學版）2015年

第10期余劍武等：數(shù)控落地銑鏜床滑枕變形有限元分析及補償

數(shù)控落地銑鏜床滑枕變形誤差按照國家標準規(guī)定的精度要求是0。03 mm/500 mm＼[2＼]，即在滑枕外伸500 mm時，滑枕的變形量不能超過0。03 mm，但是用戶會提出更高的要求以達到高精度的加工目的。如何通過對數(shù)控落地銑鏜床滑枕變形和補償?shù)难芯縼頊p少滑枕的變形，許多學者提出了不同的解決思路和方法。萬東東[3]論文中介紹了意大利帕馬公司在滑枕變形方面的研究成果，設(shè)計出了在數(shù)控落地銑鏜床滑枕上部分別安裝兩個拉桿和油壓缸裝置，通過數(shù)控系統(tǒng)控制油缸壓力的變化來補償滑枕變形。王鴻博等人[4]采用伺服系統(tǒng)裝置和配重塊對MCMG 180 LG數(shù)控落地銑鏜床的主軸滑枕進行平衡補償?？但I民＼[5＼]設(shè)計了一種采用支撐輪、下壓輪、平衡條、支撐套組成的機械平衡補償裝置專門針對主軸滑枕“低頭”現(xiàn)象。以上方法對減少滑枕變形都具有一定的效果，但還不能完全滿足高精度的加工要求。隨著有限元分析技術(shù)的成熟和計算機技術(shù)的發(fā)展，使用有限元分析軟件來解決工程中的復雜問題變得越來越方便。有限元分析軟件Ansys WB由于為CAD軟件提供了良好的數(shù)據(jù)接口和兼容性，方便產(chǎn)品研發(fā)過程中的多平臺協(xié)作，在工業(yè)產(chǎn)品設(shè)計研發(fā)中被廣泛使用。

為了解決大型數(shù)控落地銑鏜床滑枕變形問題，本文以武漢重型機床集團有限公司生產(chǎn)的TK6920B大型數(shù)控落地銑鏜床為研究對象，采用預(yù)應(yīng)力撓曲加工法和液壓拉桿法對滑枕變形進行補償，使用Ansys WB有限元分析軟件研究滑枕的變形和補償以解決實驗驗證工作量大、成本高的問題。應(yīng)用Pro/E進行建模，導入Ansys WB中分析滑枕在不同外伸量時的變形量，并使用Ansys WB檢驗預(yù)應(yīng)力撓曲加工和液壓拉桿補償方法的補償效果。

1滑枕變形有限元分析

1。1滑枕三維建模

TK6920B大型數(shù)控落地銑鏜床是武漢重型機床集團有限公司研發(fā)的一種大型機床（圖1），其滑枕的截面尺寸為480 mm×520 mm，最大工作行程為1 200 mm，滑枕自身質(zhì)量約6 000 kg。

圖1TK6920B大型數(shù)控落地銑鏜床

Fig。1TK6920B CNC floor type boring

and milling machine tool

使用Pro/E軟件建立滑枕三維實體模型，如圖2所示，建模時，嚴格按照TK6920B數(shù)控落地銑鏜床滑枕的實際尺寸進行建模。在網(wǎng)格劃分時，為了提高Ansys WB分析速度，忽略倒角等細節(jié)。

圖2 滑枕三維模型

Fig。2 The 3D model of the ram

1。2滑枕材料屬性和網(wǎng)格劃分

數(shù)控落地銑鏜床滑枕材料采用QT6003，密度為7 300 kg/m3，彈性模量173 GPa，泊松比0。3＼[6＼]。

由于滑枕的三維模型比較復雜，采用Ansys WB自動劃分方法進行網(wǎng)格劃分，最終的網(wǎng)格劃分效果如圖3所示，當滑枕的外伸量為1 200 mm時，節(jié)點數(shù)為109 704，單元數(shù)為60 242。本文采用Ansys WB中自帶的單元畸變度（Skewness）檢測工具進行網(wǎng)格劃分的質(zhì)量檢測，經(jīng)檢測，單元畸變度平均值為0。42，屬于非常好級別的網(wǎng)格。

圖3滑枕網(wǎng)格劃分

Fig。3The meshed ram

1。3確定鏜削力的大小

鏜削力的大小在不同工況下不同，與吃刀量，進給量和切削對象的材料屬性都有很大的關(guān)系，所以很難確定在動態(tài)過程中的鏜削力大小。本文使用經(jīng)驗公式計算鏜削力F0的大?。躘7＼]：

F0=2。73×109×ap×f×η。（1）

根據(jù)數(shù)控落地銑鏜床的機械性能，在粗加工中鏜削力更大，其切削參數(shù)為：切削深度ap=10 mm，進給量f=0。9 mm/r，功率常數(shù)取η=0。96。經(jīng)計算可得鏜削力大小F0=23 545 N。

1。4確定功能附件的質(zhì)量

功能附件包括主軸或者安裝在滑枕前端的銑頭、平旋盤等。功能附件的自質(zhì)量在不同加工條件下是不同的，在進行滑枕變形有限元分析時通常采用經(jīng)驗原則，大體估算在大部分工況下附件的質(zhì)量大小。在不同加工工況下，數(shù)控落地銑鏜床功能附件的質(zhì)量約在500 kg左右，因此采用功能附件的重力F1=5 000 N。

1。5滑枕變形的有限元分析

在定義好相關(guān)參數(shù)，劃分網(wǎng)格后，加載邊界條件，啟用Ansys WB對滑枕變形進行有限元分析，得到圖4所示的滑枕總變形云圖。

滑枕的工作行程最大可達到1 200 mm，對滑枕每伸出100 mm做一次變形的有限元分析，可以得到不同伸出量時的滑枕變形量，見表1。有限元分析結(jié)果表明，在沒有補償?shù)那闆r下，滑枕在外伸量比較大時會產(chǎn)生較大的變形量，這種變形誤差很難滿足用戶的加工精度要求。

圖4滑枕總變形云圖

Fig。4Total deformation contours of ram

表1沒有補償情況下滑枕的變形量

Tab。1The deformation of ram without compensation

滑枕行程/mm

100

200

300

400

500

600

滑枕變形/μm

1。4

2。6

3。8

5。1

6。8

8。8

滑枕行程/mm

700

800

900

1 000

1 100

1 200

滑枕變形/μm

11。4

15。0

21。5

30。1

40。5

51。9

2滑枕變形補償方法

傳統(tǒng)中主要使用液壓機械補償?shù)姆绞綄碜冃芜M行補償，實踐證明這種方法的缺點是，當移動部件移動速度較快時，由于液壓響應(yīng)有滯后作用，這對加工精度產(chǎn)生了很大的影響。研究表明，綜合使用預(yù)應(yīng)力撓曲加工補償法和液壓拉桿補償法對滑枕變形進行補償，不僅補償精度高，而且穩(wěn)定可靠。

2。1預(yù)應(yīng)力撓曲加工補償法

預(yù)應(yīng)力撓曲加工法補償滑枕變形主要是針對滑枕由于自身的重力產(chǎn)生的變形。

2。1。1預(yù)應(yīng)力撓曲加工法補償原理

如圖5（a）所示，滑枕移動部件重心位置的下方有滾動支撐，滑枕由于自質(zhì)量產(chǎn)生撓曲變形，變形部分是圖中的陰影部分。使用數(shù)控加工的方法將陰影部分去除掉，滑枕裝配好之后如圖5（b）所示，由于滾動支撐始終隨滑枕的移動而移動，所以滑枕一直能保持平直的狀態(tài)＼[8＼]。

2。1。2預(yù)應(yīng)力撓曲加工示意圖

預(yù)應(yīng)力撓曲加工示意圖如圖6所示，根據(jù)滑枕的工作行程狀況，通過加工預(yù)處理達到減小滑枕變量的目的。由于滑枕變形的影響因素很多，撓曲線的變形近似看成直線，圖中實線表示滑枕的加工輪廓，為裝配前的形狀，虛線表示滑枕裝配后，在滑枕自身重力作用下的理想平直狀態(tài)。

（a） 滑枕撓曲變形示意圖

（b） 預(yù)應(yīng)力撓曲加工法補償效果示意圖

圖5預(yù)應(yīng)力撓曲加工法補償原理

Fig。5Principle of prestress flexural

deformation machining

圖6滑枕預(yù)應(yīng)力加工示意圖

Fig。6Schematic prestress machining of the ram

2。1。3撓曲性補償效果

使用Ansys WB導入模型，定義相關(guān)參數(shù)，劃分網(wǎng)格，加載邊界條件進行求解，得到的分析結(jié)果如表2所示。由表2中的分析結(jié)果可以看出，相對于表1中的分析結(jié)果，預(yù)應(yīng)力撓曲加工的方法對滑枕的變形起到了一定的補償作用，但是變形量依然較大。

表2預(yù)應(yīng)力撓曲加工補償法的補償效果

Tab。2The ram deformation after prestress

flexural deformation machining

滑枕行程/mm

100

200

300

400

500

600

滑枕變形/μm

1。4

2。5

3。7

4。9

6。4

8。2

滑枕行程/mm

700

800

900

1 000

1 100

1 200

滑枕變形/μm

10。5

13。6

19。3

26。5

35。2

44。6

2。2液壓拉桿補償法

在數(shù)控銑鏜床加工過程中會使用到不同的加工附件進行加工，在不同的工況條件下也會對滑枕產(chǎn)生不同變形，液壓拉桿補償法主要用于解決在不同加工附件和工況情況下的滑枕變形問題。

2。2。1液壓拉桿法補償原理

如圖7所示，滑枕變形液壓拉桿補償裝置在滑枕上部對稱設(shè)有兩根拉桿，拉桿置于拉桿孔內(nèi)，滑枕前部臺肩與拉桿前端相連，在滑枕的后端使用拉力油缸與拉桿相連，拉桿力作用于滑枕的工作端，在拉桿拉力的作用下，滑枕工作端的變形量可以得到不同程度的減小。拉力油缸通過NC實時控制液壓拉力大小，從而保證滑枕在不同外伸條件下，拉力油缸都能提供足夠的補償拉力＼[9＼]。

圖7滑枕液壓拉桿補償原理圖

Fig。7Compensation principle of hydraulic rod

2。2。2補償拉力的計算

根據(jù)材料力學的知識可知：鏜削力F0和功能附件重力F1產(chǎn)生的變形量ω0（向下）與拉桿拉力F產(chǎn)生的補償變形ω1（向上）的大小相等。

由材料力學公式可知：

ω0=F0l3/3EI+ql4/8EI ；（2）

ω1=Fhl2/EI。 （3）

式中：l為滑枕的外伸量；E為滑枕的彈性模量；I為滑枕的慣性矩；h為拉桿到滑枕中心線的垂直距離；q為滑枕自質(zhì)量等效的均布載荷。

由ω0=ω1可得到，補償拉力的計算公式為：

F=（8F0l+3ql2）/24h。（4）

由于在滑枕外伸量較小的情況下，滑枕的變形相對較小，能滿足加工精度，可以不用液壓拉桿補償，主要補償區(qū)間為600～1 200 mm。通過計算可得不同滑枕外伸量時補償拉力的大小，見表3。

表3補償拉力大小

Tab。3The magnitude of the compensation force

滑枕行程

/mm

600

700

800

900

1 000

1 100

1 200

補償拉力

/kN

31。7

37。2

42。8

48。4

54。1

59。9

65。8

2。2。3液壓拉桿補償法有限元分析

采用液壓拉桿法對滑枕的變形進行補償，對補償后的滑枕變形進行有限元分析，加載邊界條件進行求解，得到補償后的滑枕變形量，如表4所示。

表4液壓拉桿補償后的滑枕變形量

Tab。4The ram deformation with hydraulic

rod compensation

滑枕行程/mm

100

200

300

400

500

600

滑枕變形/μm

1。4

2。5

3。7

4。9

6。4

8。0

滑枕行程/mm

700

800

900

1 000

1 100

1 200

滑枕變形/μm

9。5

11。2

13。9

17。3

20。9

24。0

由表4中的分析結(jié)果可以看出，通過預(yù)應(yīng)力撓曲變形補償和液壓拉桿補償，滑枕的變形已經(jīng)變得很小了，在外伸1 200 mm時的變形也只有24。0 μm。

2。2。4實驗檢測結(jié)果

對補償后的滑枕變形量進行實驗檢測，實際工況條件和仿真分析的工況條件相同，使用角尺、平尺、千分表對不同滑枕外伸長量下的滑枕變形量進行了測量，最大變形量為20 μm。對補償后的滑枕變形量的仿真結(jié)果和實驗檢測結(jié)果進行了比較（圖8）。由圖8可知，有限元仿真分析得到的滑枕變形量與實驗測量的滑枕變形量數(shù)據(jù)相近?？傮w來說實驗結(jié)果和仿真結(jié)果變化趨勢一致，并且誤差不大，大概在5 μm 以下，說明了仿真分析結(jié)果是可靠的，補償方法是有效的。初步分析誤差主要來源于兩個方面：首先是由于實驗條件和檢測方法帶來的實驗測量誤差，其次是仿真時由于模型簡化等原因帶來的仿真誤差。

滑枕外伸量/mm

圖8補償后的仿真結(jié)果和實驗檢測結(jié)果比較

Fig。8Comparison of simulation results and experimental

measurement after deformation compensation

3結(jié)論

在大型數(shù)控落地銑鏜床的滑枕變形研究中，由于滑枕及附件形狀復雜，傳統(tǒng)上使用材料力學方法過度簡化模型去求滑枕的變形量往往不夠精準，通常需要進行大量的實驗，并測量采集相關(guān)數(shù)據(jù)，不但工作效率低，而且耗費大量的人力物力。本文采用有限元方法研究數(shù)控落地銑鏜床滑枕變形和補償，得到以下結(jié)論：

1） 采用Ansys WB有限元分析方法研究了數(shù)控落地銑鏜床滑枕的變形量和補償，有限元仿真分析得到的滑枕變形量與實驗測量的滑枕變形量變化趨勢一致，誤差在5 μm 以下。說明了仿真分析結(jié)果是可靠的。有限元分析方法可提高效率，為大型數(shù)控銑鏜床的研發(fā)設(shè)計提供了更好的設(shè)計分析方法，可大幅縮短產(chǎn)品研發(fā)周期，在市場競爭中占得先機。

2） 綜合使用預(yù)應(yīng)力加工補償法和液壓拉桿補償法對滑枕變形進行補償，能夠大幅減少滑枕的變形量，最大變形量為20 μm左右，高于國家標準規(guī)定的精度要求0。03 mm/500 mm，有效提高了數(shù)控落地銑鏜床的加工精度，可以滿足實際生產(chǎn)中加工精度的要求。

參考文獻

[1]齊齊哈爾二機床（集團）有限公司。 滑枕撓曲變形雙向補償裝置：中國，CN101913106A[P]。2010-12-15。

Qiqihar No。2 Machine Tool （Group） Company Limited。 Bidirectional compensation device of ram deflection： China，CN101913106A[P]。2010-12-15。 （In Chinese）

[2]JB/T8490。1-96數(shù)控落地銑鏜床、落地銑鏜加工中心精度檢驗[S]。北京：中華人民共和國機械工業(yè)部， 1996。

JB/T8490。1-96CNC floor type boring and milling machine & machining center precision test[S]。Beijing： People's Republic of China Ministry of Machinery Industry， 1996。 （In Chinese）

[3]萬東東。 大型數(shù)控落地銑鏜床主軸滑枕變形分析與誤差補償技術(shù)的研究[D]。蘇州： 蘇州大學機電工程學院， 2012。

WAN Dongdong。 The research on the spindle ram deformation analysis and error compensation technology of the large NC floor type boring and milling machine[D]。 Suzhou： College of Mechanical and Electric Engineering，Soochow University， 2012。 （In Chinese）

[4]王鴻博， 阮衛(wèi)平， 王寶平，等。MCMG 180 LG落地銑鏜床主軸滑枕平衡補償控制系統(tǒng)的設(shè)計[J]。制造技術(shù)與機床， 2009（5）：92-95。

WANG Hongbo， RUAN Weiping， WANG Baoping， et al。Design of balance compensation control system of chief axis ram for the floortype boring and milling machine MCMG 180 LG[J]。 Manufacturing Technology & Machine Tool， 2009（5）：92-95。 （In Chinese）

[5]康獻民。 滑枕移動重心位置變動及補償?shù)难芯縖J]。 機械設(shè)計與制造， 2000（3）：72-73。

KANG Xianming。 Study on the gravitational center and compensation during ram movement[J]。 Machinery Design & Manufacture， 2000（3）： 72-73。 （In Chinese）

[6]李民選。 TK6920B數(shù)控落地銑鏜床滑枕變形補償精度實驗研究[D]。 長沙： 湖南大學機械與運載工程學院，2013。

LI Minxuan。 The experimental research for deformation precision compensation of ram for TK6920B CNC floor type boring and milling machine[D]。 Changsha： College of Mechanical and Vehicle Engineering，Hunan University， 2013。 （In Chinese）

[7]浦林祥。 金屬切削機床夾具設(shè)計手冊[M]。2版。北京：機械工業(yè)出版社，1995： 545-547。

PU Linxiang。 The fixture design manual of metal cutting machine tool[M]。 2nd ed。 Beijing： China Machine Press， 1995： 545-547。 （In Chinese）

[8]戴晨，劉小鵬，張文橋。 TK6916B數(shù)控落地銑鏜床補償系統(tǒng)分析[J]。 湖北工業(yè)大學學報， 2007，22（4）：5-13。

DAI Chen， LIU Xiaopeng， ZHANG Wenqiao。 Analysis of compensation systems of floor type NC boring and milling machine based on TK6916B[J]。 Journal of Hubei University of Technology， 2007， 22（4）：5-13。 （In Chinese）

[9]蕪湖恒升重型機床股份有限公司。一種滑枕撓度變形補償裝置：中國，CN 102380802A[P]。2012-03-21。

Wuhu Hengsheng Heavy Machine Tool Company Limited。 A compensation device of ram deflection：China， CN102380802A[P]。2012-03-21。 （In Chinese）