基于熱模態(tài)分析的熱誤差溫度測點(diǎn)優(yōu)化選擇

2012年4月機床與液壓Apr.2012第40卷第7期MACHINE TOOL HYDRAULICSVol. 40 No. 7DoI:10.3969/jisn.1001-3881.2012.07.001基于熱模態(tài)分析的熱誤差溫度測點(diǎn)優(yōu)化選擇張琨,張毅,侯廣鋒,楊建國(上海交通大學(xué)機械工程學(xué)院,上海200240)摘要:針對影響機床熱誤差建模的溫度場(chǎng)分布問(wèn)題,提出一種熱模態(tài)分析方法,對機床熱誤差建模溫度測點(diǎn)進(jìn)行優(yōu)化選擇。以數控機床主軸溫度場(chǎng)分析為例,利用熱模態(tài)方法得到主軸各模態(tài)的時(shí)間常數、溫度場(chǎng)及熱變形模態(tài)形狀,從而確定溫度測點(diǎn)的最優(yōu)位置。并通過(guò)實(shí)驗驗證了所建立模型的準確性與魯棒性。關(guān)鍵詞:數控機床;熱誤差;溫度測點(diǎn)選擇;熱模態(tài)中圖分類(lèi)號:TG502.15文獻標識碼:A文章編號:1001-3881(2012)7-001-3Selection of Sensor Placement for Thermal Error CompensationBased on Thermal Mode analysisZHANG Kun, ZHANG Yi, HOU Guangfeng, YANG JianguoSchool of Mechanical Engineering, Shanghai Jiaotong University, Shanghai 200240, China)Abstract: Aimed at the influence of temperature field distribution of machine tool on thermal error modeling, a method basedthermal mode analysis was proposed to select the optimum points of temperature sensors. Time constants, temperature field and thermalexpansion modal shapes of the spindle were calculated, so the optimum locations of temperature sensors were obtained. Satisfactory re-sults are achieved in terms of model accuracy and robustness through experimentKeywords: Machine tool; Thermal error; Temperature sensor placement; Thermal mode誤差補償技術(shù)是提髙機床加工精度的一項關(guān)鍵技1熱模態(tài)分析概述術(shù),在現代工業(yè)中越來(lái)越受到重視。有效實(shí)施誤差補模態(tài)分析是一種廣泛應用于結構動(dòng)力學(xué)問(wèn)題的方償的前提是建立合適的誤差補償模型。由于在實(shí)際加法,它能夠確定設計結構或機械部件的振動(dòng)特性,得工過(guò)程中很難直接測得刀尖點(diǎn)和工件之間的熱誤差,到結構的固有頻率和振型,由于它們是進(jìn)行動(dòng)力學(xué)研所以必須將熱誤差與一些容易測得的變量聯(lián)系起來(lái)建究的重要參數,所以模態(tài)分析也是更為詳細的動(dòng)力學(xué)立數學(xué)模型,最常見(jiàn)的就是熱誤差一溫度模型。由于分析如譜響應分析、瞬態(tài)動(dòng)力學(xué)分析、頻譜分析的起機床熱誤差的影響因素眾多,為了使建立的模型盡可點(diǎn)3。類(lèi)似地,模態(tài)分析方法也可以被引入到了熱傳能地具有魯棒性和精確性,如何對機床溫度測點(diǎn)進(jìn)行遞問(wèn)題的分析之中。模態(tài)分析的本質(zhì)是將復雜的動(dòng)優(yōu)化選擇成為熱誤差建模中的關(guān)鍵問(wèn)題。國內外有很態(tài)系統分解為簡(jiǎn)單的子系統,這也是有限元方法的本多學(xué)者在此方面進(jìn)行了廣泛深人的研究。最普遍的選擇溫度傳感器位置的方法是根據經(jīng)驗預先在機床多個(gè)質(zhì)屬性,因此也可以使用有限元方法對熱模態(tài)進(jìn)行分位置上布置溫度傳感器,然后使用統計學(xué)方法選出對析。機床結構熱變形的過(guò)程是首先由輸入的熱負載導方法大大加重了誤差測量和計算的工作量,且在實(shí)際致溫度場(chǎng)的變化,然后溫度場(chǎng)的不均勻變化最終導致加工中由于布線(xiàn)過(guò)多會(huì )影響機床正常工作。由于僅僅了熱變形。要從本質(zhì)上掌握機床的熱特性,就必須得從經(jīng)驗出發(fā),沒(méi)有充分考慮到機床溫度場(chǎng)的動(dòng)態(tài)特性到機床溫度場(chǎng)的分布,即對熱傳遞微分方程進(jìn)行求對熱誤差的影響,所建立的熱誤差模型很難保證很高解。瞬態(tài)熱傳遞問(wèn)題的基本微分方程為:的精度和魯棒性。作者從機床傳熱特性出發(fā),提出了CT(r)+K,T(t)=Q(r)一種基于熱模態(tài)分析的機床熱誤差測點(diǎn)優(yōu)化選擇方式中:Cr為熱容矩陣,K為熱傳遞系數矩陣,T(t)法,經(jīng)實(shí)驗驗證,提高了熱誤差模型的魯棒性為溫度向量,Q(t)為熱負載向量。v凵中國煤化工收稿日期:2011-03-11CNMHG基金項目:“高檔數控機床與基礎制造裝備”國家科技重大專(zhuān)項項目(2000X014-2作者簡(jiǎn)介:張琨(1988-),男,碩土研究生,主要從事數控機床熱誤差補償方面的研究。E- mail: wildbird2005@2機床與液壓第40卷式(1)的特征問(wèn)題是:上式中的點(diǎn),可以通過(guò)反向熱傳遞分析方法來(lái)Kr中=Cr中1A求解,使用奇異值分解方法可以得到r,的值。接式中;A為特征值對角矩陣,特征值為λ1、λ2、…、下來(lái)就可以由模態(tài)溫度來(lái)求出物理溫度、熱變形和熱λn;φr為相應的特征向量矩陣。誤差式(2)中的特征值A的物理意義是自然時(shí)間常T(t)=φ(t)數(代表第i個(gè)模態(tài)對于熱負載的響應時(shí)間)的倒U(t)=φ6(t)(10)數,即有:λ2=1/r;。類(lèi)似于結構動(dòng)態(tài)分析,每個(gè)特E(t)=φθ(t)征值或者自然時(shí)間常數都對應于一個(gè)模態(tài),最小的特為了表征對機床熱變形影響最大的熱模態(tài),引入征值或者說(shuō)最大的自然時(shí)間常數對應最低的模態(tài),從熱模態(tài)比重系數c=5r,c的大小代表了每個(gè)模態(tài)而可以得到以特征向量表示的溫度場(chǎng)模態(tài)形狀。每個(gè)所占的比重。通常熱變形過(guò)程中只有幾個(gè)模態(tài)起主要溫度場(chǎng)模態(tài)形狀反映的是在單獨的模態(tài)下溫度場(chǎng)的分作用,其余影響較小的模態(tài)可以被忽略掉。布情況,實(shí)際的瞬時(shí)溫度場(chǎng)則是這些模態(tài)的疊加。這3熱模態(tài)分析的應用些模態(tài)是機床結構的本質(zhì)屬性,相互之間不存在耦在如圖1所示的簡(jiǎn)化的車(chē)削中心主軸上進(jìn)行實(shí)驗合,且與外界輸入的熱負載無(wú)關(guān)。通過(guò)對自然時(shí)間常驗證,此車(chē)削中心主軸直接由主軸電機驅動(dòng),因此主數和模態(tài)形狀的分析,可以找到對機床結構熱變形影軸前后軸承是主要熱源,主軸熱誤差主要表現為軸右響最大的模態(tài)和關(guān)鍵點(diǎn),在相應的位置布置溫度傳感端的熱伸長(cháng)。器,從而達到對溫度測點(diǎn)進(jìn)行優(yōu)化選擇的目的。2熱模態(tài)分析過(guò)程不考慮溫度變化對材料熱特性的影響,瞬態(tài)熱傳遞的基本微分方程是一個(gè)一階耦合常系數微分方程,對于實(shí)際機床結構,其邊界條件往往很復雜,其解析解很難求得。比較簡(jiǎn)便和實(shí)際的解法是有限元模態(tài)解圖1CNC車(chē)削中心主軸主軸的一些特性參數如表1所示。在有限元模態(tài)解法中,使用φ1作為變換矩陣對表1主軸特性參數耦合微分方程進(jìn)行解耦,將實(shí)際物理溫度T(t)轉換對流傳熱系數h/(W·m2.(℃)-)2.7為模態(tài)溫度θ(t),即有:傳熱系數K(W·m-2.(℃)-)T(t)=中1O((3)密度p/(kg·m3)7.8×10由于在進(jìn)行解耦變換時(shí)與Cr相互正交,即:比熱容C/(Jkgt·(℃)1)434ΦnC,Φn=I4熱膨脹系數a((℃)-)1.06ΦK1Φ將上述三式代入式(1)有:長(cháng)度L6(t)+A0(t)=φQ(t)=6(t)為了進(jìn)行熱模態(tài)分析,對主軸進(jìn)行有限元建模。式中:5(t)為模態(tài)熱負載向量。由軸結構、熱源及傳熱條件的對稱(chēng)性,可以建立二維至此,熱負載矩陣Q(1)以有限元節點(diǎn)向量形式結構模型進(jìn)行模擬。機床先運行3h,接著(zhù)模擬中午休息停機0.5h,然后再運行3h后停機05h,直到表達出來(lái)。而且式(6)為一階常系數微分方程組,機床到達熱平衡狀態(tài)。當主軸轉速為20rymn空它可以寫(xiě)成一系列單獨變量的形式,即0(t)載切削時(shí),計算得到在前后軸承處產(chǎn)生的熱量分別為0(1)+“=(t)(7)102和85.4W?？紤]到模態(tài)初始條件θ(0)也可以由自然初始通過(guò)熱模態(tài)分析,可以得到各個(gè)模態(tài)的特征值以條件轉換而來(lái),即及熱模態(tài)比重系數。經(jīng)過(guò)計算可以發(fā)現:前4個(gè)模態(tài)(0)=T(0)的比重占到總比重的90%以上,說(shuō)明此4個(gè)模態(tài)在機床熱負載的變化可以用一系列的階躍輸入響應主軸熱變形中國煤化工模態(tài)的時(shí)間常來(lái)近似,解微分方程(8)可得數及比重如CNMHG形狀如圖2所6(t)=E,7,(1-e)示第7期張琨等:基于熱模態(tài)分析的熱誤差溫度測點(diǎn)優(yōu)化選擇表2主要4個(gè)模態(tài)時(shí)間常數及比置值為測量值與預測值之間的誤差值?？芍?熱誤差模熱模態(tài)時(shí)間常數/s比重/%型在機床預熱以及不同轉速加工情況下均能達到較好9145.4的預測效果。3242.62341318.8107642.7幅半世I測量值2一預測值0.505a-~~如“““時(shí)間/min020.40608置位置/m圖3熱誤差建模精度分析(n)模態(tài)I(b)模表Ⅱ4結論通過(guò)有限元建模,利用模態(tài)分析方法得到機床溫度場(chǎng)及熱變形模態(tài)形狀,找出與主軸熱變形相關(guān)性最大的位置,從而實(shí)現對機床熱誤差測點(diǎn)的優(yōu)化選擇。與傳統經(jīng)驗型熱誤差建模方法相比,熱模態(tài)分析方法能夠從本質(zhì)上揭示主軸的熱塑性關(guān)系。使用最小二乘法建立的線(xiàn)性熱誤差預測模型驗證了熱模態(tài)分析方法020.4060.810.20.40.60.8的準確性與魯棒性。位置/m(c)模態(tài)Ⅲ(d)模態(tài)Ⅳ參考文獻【1】楊建國潘志宏孫振勇,等.回歸正交設計在機床熱誤圖2主要4個(gè)溫度場(chǎng)模態(tài)形狀差建模中的應用[門(mén)]航空精密制造技術(shù),1999,35(5)通過(guò)分析時(shí)間常數的比重以及熱變形模態(tài)形狀圖【2】張宏韜基于神經(jīng)網(wǎng)絡(luò )的熱誤差組合建模技術(shù)及其實(shí)時(shí)形可以發(fā)現:模態(tài)Ⅱ、Ⅲ、Ⅳ占的比重更大,對應于補償應用[D].上海:上海交通大學(xué),2006【3】傅志方振動(dòng)模態(tài)分析與參數識別[M].北京:機械工業(yè)相應的模態(tài),在x1=0.15m、x2=0.5m、x=0.85m出版社,1990:17-50.處的位置溫度變化最大。因此,可以選擇這3個(gè)位置(4】 COUTINHOA LG A, LANDAU L, WROBEL L C,etal為熱誤差建模溫度測點(diǎn),將溫度傳感器布置在這3個(gè)Modal Solution of Transient Heat Conduction Utilizing位置,并使用渦流位移傳感器測量主軸熱誤差。從經(jīng)Lanczos Algorithm[ J]. Intemational Journal for Numerical驗分析出發(fā):x1、x靠近前后軸承的熱源位置,x2是Methods in Engineering, 1989, 28(1): 13-25[5] Ma Y Sensor Placement Optimization for Thermal Error Com-主軸上與熱誤差呈最大相關(guān)性的點(diǎn),由此也可以看出n Machine Tools[ D]. Ann Arbor: University of熱誤差模態(tài)分析的正確性。Michigan, 200為了驗證熱模態(tài)分析的魯棒性,在實(shí)驗中先后在【6】 FRASER S, ATTIA MH, OSMAN MOM. Control-Oriented2500、1500、3500r/min的轉速下運行120min,記Modeling of Thermal Deformation of Machine Tools Based錄三個(gè)溫度測點(diǎn)的溫度(°C)和主軸的熱誤差數值on Inverse Solution of Time- variant Thermal Loads with(μm)。根據實(shí)驗數據,采用最小二乘方法建立線(xiàn)Delayed Response[ J]. Journal of Manufacturing Scienceand Engineering,2004,126(2):286-296性模型如下:【7】傅建中陳子辰奇異值分解識別精密機械熱動(dòng)態(tài)特性Δx=5.564+0.37641+3.9454t2-0.4024t3參數的研學(xué)既,2004,38(4由該模型得到的機床熱誤差預測曲線(xiàn)如圖3所474-477中國煤化工示。其中測量值為渦流位移傳感器實(shí)際測得的主軸伸【8】楊建國CNMHG用[D].上海長(cháng)量,預測值為根據熱誤差模型計算所得的值,殘差上海交通大學(xué),1998