山東大學(xué)機(jī)械工程學(xué)院高效潔凈機(jī)械制造教育部重點(diǎn)實(shí)驗(yàn)室 路來驍 孫 杰 張 閣

中航工業(yè)成都飛機(jī)工業(yè)（集團(tuán)）有限責(zé)任公司數(shù)控加工廠 羅育果 熊青春

隨著現(xiàn)代飛機(jī)性能要求的不斷提高，飛機(jī)主承力構(gòu)件普遍采用整體化結(jié)構(gòu)設(shè)計(jì)，如整體的梁、框、壁板等。航空整體結(jié)構(gòu)件因尺寸大、薄壁特征多，在毛坯初始?xì)堄鄳?yīng)力、切削力、切削熱、裝夾力、環(huán)境溫度變化等因素的影響下，加工變形現(xiàn)象十分嚴(yán)重。環(huán)境溫度變化使工件在尺寸上伸長或縮短，嚴(yán)重制約了整體結(jié)構(gòu)件的尺寸精度。溫度補(bǔ)償成為解決工件線性變形、提高工件尺寸精度的有效途徑。

國內(nèi)外學(xué)者對變形補(bǔ)償做了大量研究工作，但針對航空整體結(jié)構(gòu)件變形補(bǔ)償?shù)难芯考杏?000年以后，研究思路可歸納為在線補(bǔ)償和離線補(bǔ)償2種方式。在線補(bǔ)償受在線監(jiān)測、實(shí)時(shí)補(bǔ)償技術(shù)的制約，應(yīng)用十分困難，目前，研究大都集中于離線補(bǔ)償方式。Law[1]研究了基于刀桿變形的腔槽加工過程中誤差補(bǔ)償方法；鄭聯(lián)語[2]提出了進(jìn)給量的局部優(yōu)化方法，定性地提出了優(yōu)化補(bǔ)償過程；王志剛[3]研究了切削力作用下側(cè)壁、腹板變形，提出了進(jìn)給量的局部優(yōu)化法與刀具偏擺數(shù)控補(bǔ)償技術(shù)；Lee[4]采用試驗(yàn)手段，對球頭立銑刀銑削側(cè)壁時(shí)的刀具路徑進(jìn)行了補(bǔ)償；S.Ratchev[5]將工件作為彈性體，建立了柔性切削力預(yù)測模型，預(yù)測了整體銑刀銑削側(cè)壁時(shí)的尺寸誤差并給出了誤差補(bǔ)償算法；孫杰[6]采用有限元軟件對薄壁框進(jìn)行了加工變形預(yù)測，基于加工變形值，在考慮一定回彈量基礎(chǔ)上對刀位點(diǎn)進(jìn)行修正以實(shí)現(xiàn)補(bǔ)償；Philippe Deépince[7]將刀具作為彈性體研究了整體立銑刀銑削加工過程零件表面尺寸誤差，在不降低生產(chǎn)率與加工精度的前提下，提出了一種考慮工件尺寸公差的補(bǔ)償算法；康永剛[8]給出了基于切削力信號的實(shí)際切深計(jì)算方法,并基于實(shí)際切深進(jìn)行了最大變形誤差的預(yù)測,為實(shí)現(xiàn)在線誤差補(bǔ)償打下了基礎(chǔ)；陳蔚芳[9]對多層補(bǔ)償算法進(jìn)行了優(yōu)化，解決了加工變形與補(bǔ)償量之間耦合問題；武輝[10]通過研究薄壁件加工過程中的彈性變形問題，在建立銑削力模型基礎(chǔ)上，提出了一種基于APDL變形預(yù)測的補(bǔ)償方法。

上述研究對變形的預(yù)測方法分為：刀具作為彈性體，工件作為剛體；刀具作為剛體，工件作為彈性體；刀具與工件都為彈性體。無論使用何種變形預(yù)測方法，對應(yīng)的補(bǔ)償方法大都針對局部讓刀變形，對整體結(jié)構(gòu)宏觀變形補(bǔ)償研究較少。針對零件局部特征的補(bǔ)償很難應(yīng)用于特征復(fù)雜、尺寸龐大的整體結(jié)構(gòu)件。

本文針對環(huán)境溫度變化對整體結(jié)構(gòu)件尺寸誤差影響問題，通過測定鋁合金毛坯材料熱膨脹系數(shù)，建立了溫度補(bǔ)償系統(tǒng)，實(shí)現(xiàn)了整體結(jié)構(gòu)件在特定溫度下的準(zhǔn)確尺寸重構(gòu)。

1 系統(tǒng)開發(fā)工具

1.1 CAA簡介

CAA(Component Application Architecture)組件應(yīng)用架構(gòu)是DS產(chǎn)品擴(kuò)展和客戶進(jìn)行二次開發(fā)的強(qiáng)有力工具。CAA建立在面向?qū)ο蟪绦蛟O(shè)計(jì)基礎(chǔ)上，使用COM（Component Object Model）、OLE（Object Linking and Embedding）技術(shù)，強(qiáng)化了開發(fā)功能。CAA使用RADE(Rapid Application Development Environment )快速開發(fā)平臺，以Visual Studio 2005為載體，通過加載RADE插件，在C++環(huán)境下對CATIA進(jìn)行二次開發(fā)，實(shí)現(xiàn)了與CATIA的無縫結(jié)合。

1.2 CAA補(bǔ)償函數(shù)

在室溫范圍內(nèi)，溫度變化對鋁合金整體結(jié)構(gòu)件尺寸誤差的影響是線性的，借助CATIA軟件變換特征中縮放命令可以實(shí)現(xiàn)對航空整體結(jié)構(gòu)件線性變形的補(bǔ)償?？s放命令作用為使用點(diǎn)、平面或曲面表面為參考對實(shí)體進(jìn)行放大或縮小操作。

CAA中 縮 放 函 數(shù) 為：CreateScaling(const CATISpecObject_var ihToScale, const CATISpecObject- _var ihReference, const CATICkeParm_var ihRatio)。對縮放函數(shù)中各參數(shù)具體說明，參數(shù)1 ihToScale為要進(jìn)行縮放的對象，該對象類型為CATISpecObject，系統(tǒng)中縮放操作的對象為工件的所有幾何特征；參數(shù)2 ihReference是縮放的參考元素，溫度補(bǔ)償系統(tǒng)中選取平面作為參考元素；參數(shù)3 ihRatio是縮放系數(shù)，縮放系數(shù)類型為CATICkeParm，在溫度補(bǔ)償系統(tǒng)中為了保證精度，縮放系數(shù)精確到小數(shù)點(diǎn)后8位。

2 溫度補(bǔ)償系統(tǒng)總體框架

圖1是溫度補(bǔ)償系統(tǒng)組成部分說明。溫度補(bǔ)償系統(tǒng)包括：熱膨脹系數(shù)查詢模塊、零件重構(gòu)模塊。2個模塊關(guān)系如圖1所示，查詢熱膨脹系數(shù)之后，首先在輸出界面將熱膨脹系數(shù)值返回給用戶，隨后該值被轉(zhuǎn)化為縮放系數(shù)傳遞到零件重構(gòu)模塊中。

在熱膨脹系數(shù)查詢模塊中，將試驗(yàn)測定的鋁合金7050-T7451和2124-T851兩種材料熱膨脹系數(shù)離散值擬合為6階多項(xiàng)式，建立溫度與熱膨脹系數(shù)之間函數(shù)關(guān)系，用戶只需輸入溫度值、選擇材料類型即可查詢到特定溫度下的熱膨脹系數(shù)值。

零件重構(gòu)模塊首先獲取熱膨脹系數(shù)查詢模塊的計(jì)算結(jié)果，計(jì)算縮放系數(shù)值，接收輸入條件如零件類型、縮放方向等，實(shí)現(xiàn)工件在特定溫度下的尺寸重構(gòu)，為數(shù)控編程提供準(zhǔn)確尺寸的CAD模型。

圖1 溫度補(bǔ)償系統(tǒng)流程Fig.1 Process of temperature compensation system

3 功能模塊實(shí)現(xiàn)

3.1 熱膨脹系數(shù)查詢模塊實(shí)現(xiàn)

3.1.1 熱膨脹系數(shù)表達(dá)

在溫度t下，與溫度變化1℃相應(yīng)的熱膨脹系數(shù)值定義為瞬時(shí)熱膨脹系數(shù)，即：

式（1）中，LT為溫度t時(shí)試件的長度。若一定溫度范圍內(nèi)熱膨脹系數(shù)的變化很小，公式（1）還可表示為：

式（2）中，LT0為工件初始長度，Δt為溫度變化量，ΔL為尺寸變化量，αm為平均熱膨脹系數(shù)。

為了保證整體結(jié)構(gòu)件尺寸重構(gòu)的精度，使用瞬時(shí)熱膨脹系數(shù)進(jìn)行計(jì)算。目前，使用瞬時(shí)熱膨脹系數(shù)計(jì)算時(shí)，無法直接進(jìn)行積分運(yùn)算，需將瞬時(shí)熱膨脹系數(shù)轉(zhuǎn)化為多項(xiàng)式進(jìn)行計(jì)算。采用公式（1）建立熱膨脹系數(shù)與溫度的數(shù)學(xué)理論模型，數(shù)據(jù)擬合后采用6階多項(xiàng)式的計(jì)算精度最高，2124-T851預(yù)拉伸鋁合金毛坯在0～400℃范圍內(nèi)瞬時(shí)熱膨脹系數(shù)擬合的多項(xiàng)式如表1所示：

表1中XCTE、YCTE、ZCTE分別為預(yù)拉伸方向、寬度方向、高度方向熱膨脹系數(shù)值，t為環(huán)境溫度值。

3.1.2 熱膨脹系數(shù)查詢模塊建立

在溫度補(bǔ)償系統(tǒng)中建立熱膨脹系數(shù)查詢模塊，以溫度值和材料類型為輸入條件，按照溫度與熱膨脹系數(shù)的數(shù)學(xué)理論模型進(jìn)行計(jì)算，輸出特定溫度下，某種鋁合金材料的熱膨脹系數(shù)值。

新建名為QueryDialog的Module，CAA中控件的添加方式與MFC類似。建立好的熱膨脹系數(shù)查詢模塊如圖2所示，模塊中包含Editor、Radio Button、Pushbutton控件。其中，Editor用于接收輸入的溫度值和輸出計(jì)算的熱膨脹系數(shù)值，Radio Button用于選擇材料類型，Pushbutton用于對熱膨脹系數(shù)查詢模塊的響應(yīng)。

表1 Al2124-T851預(yù)拉伸板鋁合金板材不同方向熱膨脹系數(shù)擬合公式

圖2 熱膨脹系數(shù)查詢模塊Fig.2 CTE querying module

3.2 零件重構(gòu)模塊實(shí)現(xiàn)

3.2.1 縮放系數(shù)表達(dá)

CATIA中縮放系數(shù)定義為工件在某方向上變形后尺寸與原尺寸的比值，即：

公式可變形為：

綜合公式（3）（4）得溫度與縮放系數(shù)的數(shù)學(xué)關(guān)系為：

式中，b為縮放系數(shù)，T為環(huán)境溫度，T0為工件設(shè)計(jì)溫度，f（t）為熱膨脹系數(shù)值，是溫度T的函數(shù)。

3.2.2 零件重構(gòu)模塊建立

零件重構(gòu)模塊流程如圖3所示，實(shí)現(xiàn)零件重構(gòu)的關(guān)鍵是：功能接口實(shí)現(xiàn)和零件特征遍歷。

表2列出了重構(gòu)模塊中所有使用過的接口及功能，這些接口的實(shí)現(xiàn)保證了零件重構(gòu)功能的實(shí)現(xiàn)。功能接口在定義完成后，必須馬上查詢是否成功得到該接口。如果成功獲取該接口，程序繼續(xù)向下執(zhí)行；若獲取接口失敗，則輸出一個自行定義的參數(shù)，方便供編程人員查錯。以管理參數(shù)與函數(shù)關(guān)系的接口CATICkeParm-Factory為例，說明實(shí)現(xiàn)功能接口的關(guān)鍵代碼：

圖3 零件重構(gòu)模塊流程圖Fig.3 Flowchart of dimensional reconstruction module

表2 功能接口表

零件重構(gòu)是對特征樹中所有的幾何特征進(jìn)行重構(gòu)。對于模型較為簡單的零件，手工編寫代碼也可完成重構(gòu)；但對于整體結(jié)構(gòu)件，其模型特征非常復(fù)雜，某些零件甚至達(dá)到上百個特征，手工編程幾乎不可能。為解決此問題，開發(fā)了如圖4所示的高效遍歷算法，提高了編程效率。使用CATIDescendants接口獲取特征集合，循環(huán)對特征賦值，執(zhí)行縮放操作。關(guān)鍵代碼為：

圖4 遍歷零件特征流程Fig.4 Process of traverse part feature

4 應(yīng)用實(shí)例

零件重構(gòu)模塊開發(fā)完成后如圖5所示，以典型航空整體結(jié)構(gòu)梁為例，演示環(huán)境溫度為30.0℃時(shí)，2124-T851鋁合金民機(jī)梁溫度補(bǔ)償效果。

鋁合金民機(jī)梁結(jié)構(gòu)如圖6所示，工件設(shè)計(jì)尺寸(20℃尺寸)為6216.146 mm×266.717mm×63.500mm,長寬比超過了20倍。主要結(jié)構(gòu)特點(diǎn)為雙面不對稱框體結(jié)構(gòu)，一面結(jié)構(gòu)簡單框少，一面結(jié)構(gòu)復(fù)雜框多，有緣條、筋條、腹板、下陷、腹板孔、外形側(cè)面孔等結(jié)構(gòu)特征。

工件測量室溫度恒定在20℃，但現(xiàn)場加工溫度受天氣條件影響在15～30℃范圍內(nèi)變化。在實(shí)際生產(chǎn)中，經(jīng)常會遇到因環(huán)境溫度變化導(dǎo)致現(xiàn)場加工時(shí)零件尺寸合格,但運(yùn)抵測量室測量時(shí)出現(xiàn)變形超差的情況。通過溫度變形補(bǔ)償系統(tǒng)的開發(fā)，能夠?qū)崿F(xiàn)對此變形量的預(yù)先補(bǔ)償。

圖5 零件重構(gòu)模塊Fig.5 Dimensional reconstruction module

圖6 典型航空整體結(jié)構(gòu)梁Fig.6 Typical aircraft integrated stucture beams

對于該民機(jī)梁，如果不進(jìn)行溫度變形補(bǔ)償，工件從20℃升高到30℃后，尺寸變?yōu)椋?217.398mm×266.768mm×63.512mm。寬度和高度方向的尺寸誤差小于0.1mm滿足公差范圍要求，可以忽略溫度變化對寬度和高度方向的尺寸誤差；預(yù)拉伸方向尺寸誤差為1.252mm嚴(yán)重超出了誤差范圍。采用溫度補(bǔ)償之后，對此變形量進(jìn)行預(yù)先補(bǔ)償，提高了工件的尺寸精度。

5 結(jié)論

（1）完成了兩種鋁合金材料熱膨脹系數(shù)值的精準(zhǔn)表達(dá)，并建立了熱膨脹系數(shù)、縮放系數(shù)與溫度的關(guān)系模型。

（2）采用CAA方式對CATIA進(jìn)行二次開發(fā)，建立了溫度補(bǔ)償系統(tǒng)，實(shí)現(xiàn)了特定溫度下工件尺寸重構(gòu)，為特定環(huán)境下數(shù)控編程提供精準(zhǔn)CAD模型。

[1]Law K M Y, Geddam A. Prediction of contour accuracy in the end milling of pockets. Journal of Materials Processing Technology,2001(113)：399-405.

[2]鄭聯(lián)語, 汪叔淳. 薄壁零件數(shù)控加工工藝質(zhì)量改進(jìn)方法. 航空學(xué)報(bào) , 2001, 22(5)：424-428.

[3]王志剛,何寧. 薄壁零件加工變形分析及控制方案. 中國機(jī)械工程 , 2002, 13(2)：114-117.

[4]Lee C M , Kim S W, Lee Y H. The optimal cutter orientation in ball end milling of cantilever shaped thin plate. Journal of Materials Processing Technology, 2004, 153-154.

[5]Ratchev S, Liu S, Becker A A. Error Compensation strategy in milling flexible thin-wall parts, Journal of Materials Processing Technology,2005, 162/163 ： 673-681.

[6]孫杰,柯映林. 殘余應(yīng)力對航空整體結(jié)構(gòu)件加工變形的影響分析. 機(jī)械工程學(xué)報(bào), 2005, 41(2)：117-122.

[7]Philippe.Dépincé, Jean-Yves Hascoёt. Active intergration of tool deflection effects in end milling. Part 2 Compensation of tool deflection.International Journal of Machine Tools & Manufactu-re, 2005(46)：945-956.

[8]康永剛,王仲奇,吳建軍. 基于實(shí)際切深的薄壁件加工變形誤差的預(yù)測. 西北工業(yè)大學(xué)學(xué)報(bào), 2007, 25(2)：251-256.

[9]陳蔚芳,樓佩煌,陳華. 薄壁件加工變形主動補(bǔ)償方法. 航空學(xué)報(bào) , 2009, 30(3)：570-576.

[10]武輝,劉維偉,李曉燕,等. 基于APDL的薄壁件加工變形補(bǔ)償方法. 航空精密制造技術(shù), 2012, 48(4)：35-38.