張志鵬，劉康，曾祥兵，劉光磊，郭豐

(1.四川理工學(xué)院機械工程學(xué)院，四川自貢643000;2.四川長征機床集團(tuán)有限責(zé)任公司，四川自貢643000)

0 前言

隨著對數(shù)控機床加工精度的要求日益提高，誤差補償技術(shù)已經(jīng)成為提高機床精度的重要方法。幾何誤差、熱誤差、力誤差是數(shù)控機床的3大主要誤差，占總誤差的70% ～90%［1］，其中在精密加工過程中，熱誤差占總誤差的40% ～70%［2］，越是精密的機床，熱誤差所占比例越大。因此，進(jìn)行機床主軸熱誤差補償，對提高數(shù)控機床加工精度有重要意義。幾何誤差屬于靜態(tài)誤差［3］，經(jīng)過多年的研究，已經(jīng)取得了很多成果，熱誤差屬于動態(tài)誤差，簡單的誤差補償器無法隨著時間或其他因素的改變而實時在線補償。通用計算機可以實現(xiàn)一臺甚至多臺［4］機床的誤差補償，但實現(xiàn)成本較高。嵌入式裝置是以應(yīng)用為中心、軟硬件可裁剪的專用計算機，具有成本低、體積小、實時性高和界面友好的特點。隨著嵌入式設(shè)備應(yīng)用越來越廣泛［5］，嵌入式技術(shù)結(jié)合機床熱誤差補償研究取得了一些的進(jìn)展［6］，但只是針對某一特定型號機床和在特定環(huán)境下進(jìn)行的熱誤差補償研究，適用范圍有限。本文作者以ARM11開發(fā)板為主體，設(shè)計了3個模塊，在Fanuc機床上進(jìn)行了實驗，使用機床提供的外部機械原點偏移［7］功能，對Fanuc機床主軸熱誤差實時在線補償研究。

1 溫度采集分析模塊

該模塊主要功能是優(yōu)化溫度關(guān)鍵點，確定關(guān)鍵點后進(jìn)行實時補償，利用串口一將溫度值送給OK6410A開發(fā)板。通過對KVC650E立式加工中心的結(jié)構(gòu)特點和熱源分布進(jìn)行分析和初步試驗，在床身上對影響主軸熱誤差較大的位置分布了溫度傳感器，將數(shù)據(jù)采集到PC機上，進(jìn)行優(yōu)化分析，溫度變量的優(yōu)化選擇過程，其優(yōu)化過程如下:

(1)根據(jù)機床分析和基本原則初步布置測溫點;對溫度變量進(jìn)行模糊聚類，通過計算F統(tǒng)計量優(yōu)選閾值來確定最佳分組;

(2)對各溫度變量與熱誤差數(shù)據(jù)進(jìn)行灰色綜合關(guān)聯(lián)分析，選取各組中與熱誤差灰色綜合關(guān)聯(lián)度最大的溫度變量，并考慮所選溫度變量是否滿足實際情況和條件，若滿足，最終確定用于熱誤差建模的溫度變量，若不滿足，返回重選閾值。

(3)確定溫度關(guān)鍵點，溫度采集模塊一如圖1所示。

圖1 溫度采集分析模塊

通過溫度關(guān)鍵點優(yōu)化，可以消除溫度測點之間的相互耦合作用，增強熱誤差模型的魯棒性，而且還能減少溫度數(shù)據(jù)的檢測數(shù)量和簡化模型的復(fù)雜度。溫度采集模塊的硬件包括STC90C51RD芯片和帶有1602液晶顯示器的單片機，OK6410A實時顯示溫度傳感器數(shù)值，8路DS18B20溫度傳感器。經(jīng)過模塊化的設(shè)計，若要改變溫度傳感器的數(shù)目，單片機程序只需根據(jù)不同的溫度傳感器序列號，增減相應(yīng)的程序，并在計算機上進(jìn)行相關(guān)的優(yōu)化分析。

2 熱誤差建模及實時補償模塊

2.1 建模理論

熱誤差建模是熱誤差補償研究的關(guān)鍵技術(shù)，模型的預(yù)測精度直接影響到熱誤差補償系統(tǒng)的補償效果，國內(nèi)外學(xué)者進(jìn)行了大量研究，多元線性回歸 (MLR，Multiple Linear Regression)是行之有效的快速建模方法之一［8-9］。在不同的機床類型和不同的時間內(nèi)，每種建模都有各自的適應(yīng)范圍，模型大多不同，開發(fā)板的存儲空間比較大，并且可以擴(kuò)展存儲空間，這就為保存多個模型數(shù)據(jù)提供了條件。為了驗證補償器的可行性，首先將多元線性回歸算法移植到開發(fā)板，然后選擇具有更好補償效果的模型。多元線性回歸基于最小二乘原理，以因變量與自變量之間成線性關(guān)系的假設(shè)為基礎(chǔ)，通過實驗，得到多個測溫點溫度與機床單個方向上的熱誤差變化輸出關(guān)系。假設(shè)模型中熱誤差(因變量)Z與 m個溫度測點 (自變量)T1，T2，…，Tm存在線性的內(nèi)在聯(lián)系，通過實驗每個溫度測點得到n組數(shù)據(jù):T1k，T2k，…，Tmk，k＝1，2，…，n，則可得到如下關(guān)系式:

多元線性回歸模型關(guān)系式為:Z＝TC＋μ

其中:式中:Zi為機床某一方向第i次測得的熱誤差值;Ci為回歸系數(shù);Tki為機床上第k個傳感器第i次測得的溫度;μi為相互獨立且服從同一正態(tài)分布的隨機變量;Z與T為已知參數(shù)，利用最小二乘法可以求出計算參數(shù)C。

設(shè)c0，c1，…，cm分別為參數(shù) C0，C1，…，Cm的最小二乘估計，則多元線性回歸方程為:

2.2 硬件資源

(1)SamsungS3C6410處理器，ARM1176JZF-S內(nèi)核，主頻533 MHz/667 MHz;

(2)128M字節(jié)Mobile DDR內(nèi)存;(3)1G字節(jié)NAND Flash(MLC);(4)12 MHz、48 MHz、27 MHz、32.768 kHz 時鐘源;

(5)內(nèi)部實時鐘，帶有后備鋰電池座，斷電后系統(tǒng)時間不丟失;

(6)共4個串口，包括1個五線RS232電平串口和3個三線TTL電平串口及擴(kuò)展RS232串口;

(7)1個USB HOST插口，支持USB1.1，可插鼠標(biāo)、U盤等;

(8)1個100M網(wǎng)口，采用DM9000AE，帶連接和傳輸指示燈;

(9)3.5寸 LCD，支持 TFT LCD，支持10寸LVDS液晶屏;

(10)核心板為6層PCB設(shè)計，性能穩(wěn)定，經(jīng)過強電磁環(huán)境考驗。

以上的硬件資源表明:OK6410開發(fā)板為機床主軸熱誤差實時補償研究提供了良好的平臺。

2.3 系統(tǒng)平臺的搭建

系統(tǒng)平臺的搭建是實現(xiàn)一個信息處理終端，為OK6410A開發(fā)板搭建Linux系統(tǒng)，包括以下4個部分［10-11］:

(1)引導(dǎo)裝載程序 (BootLoader)，這里選用uboot1.1.6版，系統(tǒng)上電后首先被執(zhí)行，對CPU、內(nèi)存等進(jìn)行初始化，完成內(nèi)核映像的裝載和引導(dǎo);

(2)Linux內(nèi)核，采用官方的Linux3.0.1版本，配置編譯后生成zImage內(nèi)核映像文件，并制作yaffs2文件系統(tǒng)映像，利用SD卡一鍵安裝Linux系統(tǒng);

(3)圖形用戶界面，Qtopia是Trolltech公司為采用嵌入式Linux操作系統(tǒng)的消費電子設(shè)備而開發(fā)的綜合應(yīng)用平臺，選用功能強大且實用的Qtopia2.2.0版本。

(4)在Qtcreator軟件中編寫并交叉編譯應(yīng)用程序，QT體系開發(fā)模型如圖2所示。

圖2 QT的體系開發(fā)模型

2.4 軟件設(shè)計

軟件設(shè)計主要在Linux系統(tǒng)下完成，虛擬機選用VMware Workstation10，系統(tǒng)使用Ubuntu 12.04版本，安裝QtCreator4.8.1后，即可編寫QT程序，主要是基于QT下的熱誤差建模編程、通信編程和應(yīng)用程序的交叉編譯及移植［1-2］，同時還需在開發(fā)板上設(shè)置環(huán)境變量，可以利用觸摸筆打開程序或在開發(fā)板上設(shè)置開機啟動熱誤差補償程序，點擊“開始補償”按鈕后，便可以進(jìn)行熱誤差補償，開發(fā)流程如圖3所示。

圖3 QT開發(fā)流程圖

為了增強該模塊的通用性，機床類型包含F(xiàn)anuc和Siemens 2類機床，誤差補償類型設(shè)計了幾何誤差、熱誤差、切削力誤差和綜合誤差模型算法，模型包括:經(jīng)典的多元線性回歸算法、應(yīng)用廣泛且具有較強魯棒性和容錯性的BP神經(jīng)網(wǎng)絡(luò)［13］。經(jīng)過移植QT程序，開發(fā)板具有良好的人機交互界面，便于觀察傳感器溫度值和X、Y、Z三軸的補償量。設(shè)計的圖形用戶界面 (GUI，Graphic User Interface)如圖4所示。

圖4 GUI界面結(jié)構(gòu)圖

3 機床接收數(shù)據(jù)模塊

經(jīng)分析配置Fanuc 0i Mate-MC系統(tǒng)的 KVC650E立式加工中心和配置Siemens 840D系統(tǒng)的GMC4000H/2龍門加工中心的通信接口，分別制定了兩種通信補償方案:西門子機床需要解決的是內(nèi)嵌式S7-300PLC與ARM11以太網(wǎng)通信，是基于QT下的Socket通信［5］。文中主要介紹Fanuc 0i Mate-MC系統(tǒng)KVC650E立式加工中心的通信接口，經(jīng)查閱資料，利用擴(kuò)展板將X輸入開關(guān)量信號與串口控制器進(jìn)行點對點連接，ARM11經(jīng)串口二送出的誤差值送給X地址，根據(jù)機床功能手冊，在LADDER3軟件編寫對應(yīng)的讀取外部數(shù)據(jù)部分梯形圖，利用M-Card將梯形圖與機床進(jìn)行梯形圖轉(zhuǎn)換，模塊設(shè)計流程如圖5所示。

圖5 模塊三流程圖

4 熱誤差實時補償實驗

為了驗證該補償器的可行性及補償效果，以多元線性回歸算法建立數(shù)學(xué)模型，以主軸軸向為例，在KVC650E立式加工中心上進(jìn)行實驗，整個過程流程如圖6所示。

圖6 整體流程圖

將采集8路溫度數(shù)據(jù)及檢測Z軸軸向熱變形值采集到電腦上進(jìn)行分析，利用MATLAB數(shù)學(xué)工具，得到各溫度傳感器的變化曲線如圖7所示，主軸熱變形曲線如圖8所示。

圖7 各溫度傳感器的溫度變化曲線

圖8 主軸軸向熱變形曲線

首先計算出各溫度變量間的相關(guān)系數(shù)矩陣R，然后計算R的傳遞包t(R)，得到模糊等價矩陣，借助F統(tǒng)計量對溫度變量進(jìn)行分組，經(jīng)分析取閾值λ＝0.987，此時8個溫度變量共分為4組，分組結(jié)果如表1所示。

表1 溫度變量的分組

根據(jù)灰色綜合關(guān)聯(lián)度分析，由兩序列間的灰色綜合關(guān)聯(lián)度定義

分別計算出各溫度變量與主軸軸向熱誤差的灰色綜合關(guān)聯(lián)系數(shù)。按灰色綜合關(guān)聯(lián)系數(shù)大小進(jìn)行排序，如表2所示。

根據(jù)模糊聚類的分組結(jié)果，從每組中選擇溫度測點數(shù)據(jù)與熱誤差灰色綜合關(guān)聯(lián)最大的值，則應(yīng)選擇T1、T2、T6、T8作為溫度關(guān)鍵點。通過表2，最終確定T2、T6、T8作為溫度關(guān)鍵點，為了減小了模型復(fù)雜度，提升模型的精確性和魯棒性，采用實時相對溫度作為自變量，建立熱誤差補償模型，確定多元線性回歸公式為:

利用MATLAB工具分析熱誤差測量值、MLR模型預(yù)測值及殘差，結(jié)果如圖9所示。

圖9 MLR補償誤差

圖9 中模型的最大殘差絕對值為2.687 8μm。計算殘差均方根為S＝1.245 1μm。

圖10為熱誤差補償實驗現(xiàn)場，可以在機床的顯示屏上看到傳入的G地址數(shù)據(jù)，利用千分表測量機床軸向位移隨著數(shù)據(jù)的改變而移動，從而證明了補償器的可行性。

圖10 熱誤差補償實驗現(xiàn)場

利用千分表測量主軸軸向熱誤差，補償預(yù)測值和實際測得Z軸變形量如圖11所示，經(jīng)過熱誤差實時補償實驗，誤差值控制在-10～4μm范圍內(nèi)，軸向熱誤差從補償前的24μm降低到14μm，補償效果達(dá)到40%，從而證明了該補償器能夠取得良好的補償精度。

圖11 實驗對比熱誤差補償效果圖

5 結(jié)論

(1)介紹了一種基于ARM11機床熱誤差補償器，將多元線性回歸預(yù)測模型嵌入到補償器中，在KVC650E立式加工中心上實驗，證明了該補償器的可行性及良好的補償效果。

(2)實驗表明，采用多元線性回歸建模的方法，能夠達(dá)到一定的補償效果。為了進(jìn)一步提高誤差控制精度，可以根據(jù)實驗情況將幾個線性回歸模型適當(dāng)組合或補償器中嵌入預(yù)測效果更好的模型，比如神經(jīng)網(wǎng)絡(luò)模型，將誤差控制在更小的范圍內(nèi)。

(3)該補償器具有成本低、實時性高、通用性強的特點，隨著研究不斷的完善，能夠進(jìn)一步開發(fā)成商業(yè)產(chǎn)品。因此，該補償器的研究具有重要參考價值和實際意義。

［1］張宏韜.雙轉(zhuǎn)臺五軸數(shù)控機床誤差的動態(tài)實時補償研究［D］.上海:上海交通大學(xué)，2011.

［2］BRYAN J.International Status of Thermal Error Research［J］.Annals of the CIRP，1990，39(2):645-656.

［3］曹永潔.基于激光測試技術(shù)的數(shù)控機床誤差辨識與補償研究［D］.杭州:浙江大學(xué)，2006.

［4］張毅.數(shù)控機床誤差測量、建模及網(wǎng)絡(luò)群控實時補償系統(tǒng)研究［D］.上海:上海交通大學(xué)，2013.

［5］時巖.嵌入式技術(shù)發(fā)展?fàn)顩r綜述［J］.電子技術(shù)，2012(5):4-5.

［6］任兵，任小紅，徐衛(wèi)東，等.基于ARM_Linux的數(shù)控機床熱誤差補償控制器設(shè)計［J］.組合機床與自動化加工技術(shù)，2010(4):25-31.

［7］任永強，楊建國，羅磊等.基于外部機床坐標(biāo)系偏移的熱誤差實時補償［J］.中國機械工程，2003，14(14):1243-1246.

［8］RAMESH R，MANNAN M A.Error Compensation in Machine Tools a Review(Part II):Thermal Errors.International［J］.Journal of Machine Tool ＆ Manufacture，2000，40(9):1257-1284.

［9］YUAN J，NI J.The real-time Error Compensation Technique for CNC Machining Systems［J］.Mechatronics，1998，8(4):359-380.

［10］飛凌嵌入式.《OK6410 QT移植百科全書》［EB/OL］.［2014-1］.www.witech.com.cn.

［11］張霄.基于ARM11的嵌入式linux系統(tǒng)移植與實現(xiàn)［D］.濟(jì)南:山東大學(xué)，2013.

［12］蔡志明，盧傳富，李立夏.精通QT4編程［M］.2版.北京:電子工業(yè)出版社，2008.

［13］王健偉，宋執(zhí)環(huán).基于嵌入式系統(tǒng)的神經(jīng)網(wǎng)絡(luò)在線訓(xùn)練平臺實現(xiàn)［J］.傳感器與微系統(tǒng)，2010，29(8):100-103.