基于曲率擬合及變刀觸點(diǎn)的防干涉四軸磨削算法

毛 磊 何 衛(wèi) 周 潔 駱廣進(jìn) 胡旭曉 周哲平

1.浙江大學(xué)，杭州，310027 2.杭州機(jī)床集團(tuán)有限公司，杭州，310022

0 引言

復(fù)雜曲面被廣泛應(yīng)用于航天航空、造船、汽車、玻璃儀器、陶瓷以及模具等行業(yè)[1]。目前，復(fù)雜曲面的加工以五軸數(shù)控銑床加工為主。在五軸數(shù)控銑床加工中，刀具干涉的檢查與避免算法是復(fù)雜曲面加工的關(guān)鍵問題。國內(nèi)外學(xué)者對此進(jìn)行了很多研究。Lin等[2]通過比較刀觸點(diǎn)處切平面上所有方向的加工表面的法向曲率和刀具包絡(luò)面的法向曲率來找到無干涉的加工路徑。Chen等[3]研究了R銑刀加工復(fù)雜曲面時避免兩種局部干涉的方法。

在復(fù)雜曲面磨削領(lǐng)域，相關(guān)的研究相對較少，目前較常見的復(fù)雜曲面磨削還是以手工拋磨為主，比較先進(jìn)的方法有六軸砂帶磨削[4]，此外機(jī)器人和機(jī)械手在曲面拋磨加工中也開始得到應(yīng)用[5－6]。然而這些先進(jìn)的磨削加工方法都需要很高的投資成本，且其加工算法也十分復(fù)雜。本文提出一種可以在四軸聯(lián)動成形數(shù)控磨床上加工無曲率突變的自由曲面的無干涉加工算法。

1 磨削加工設(shè)備及干涉

1.1 加工設(shè)備及圓弧面砂輪結(jié)構(gòu)

1.1.1 加工設(shè)備及其聯(lián)動關(guān)系

本文涉及的機(jī)床是合作公司的七軸五聯(lián)動數(shù)控磨床，其結(jié)構(gòu)如圖1所示。該成形磨床的五軸聯(lián)動采用4＋1方式，即4個直線軸X、Y、Z、V（其中V軸是砂輪修整軸，本算法中暫不涉及砂輪修整）、1個回轉(zhuǎn)軸（B軸），另外回轉(zhuǎn)臺傾角為定值，由A軸驅(qū)動，坐標(biāo)系滿足右手定則。

圖1 七軸五聯(lián)動數(shù)控磨床的坐標(biāo)系

加工時，將葉片固定在回轉(zhuǎn)臺上，保證葉片的一條長邊與Z軸平行。加工的運(yùn)動分配關(guān)系為：磨頭完成直線運(yùn)動ΔY，工件安裝在回轉(zhuǎn)臺上完成直線運(yùn)動ΔX′、ΔZ′以及旋轉(zhuǎn)運(yùn)動ΔB′。本算法不涉及砂輪修整軸，故實(shí)際上該機(jī)床結(jié)構(gòu)實(shí)現(xiàn)四軸聯(lián)動。

1.1.2 圓弧面砂輪結(jié)構(gòu)

采用的砂輪結(jié)構(gòu)簡圖見圖2。該砂輪采用圓弧母線圓弧面結(jié)構(gòu)。

圖2 圓弧母線砂輪結(jié)構(gòu)簡圖

在防干涉加工過程中，利用砂輪的圓弧表面來避免軸向干涉，同時以砂輪最大端面半徑作為過切干涉的基準(zhǔn)半徑。由于砂輪半徑越小越不容易發(fā)生干涉，所以本算法中不考慮砂輪修整的影響。

1.2 磨削干涉

將磨削加工中的干涉分為全局干涉和局部干涉。全局干涉是指砂輪端面或軸桿與加工表面以及加工環(huán)境中的夾具等發(fā)生碰撞。局部干涉又可分為過切干涉和軸向干涉。其中，過切干涉是由于局部加工表面曲率半徑小于砂輪半徑而引起的切除多余材料的情形，如圖3所示；軸向干涉是由砂輪厚度引起的，或者在本文涉及的圓弧面砂輪中，是由刀觸點(diǎn)法向量不一致引起的，如圖4所示。實(shí)際上將砂輪母線改成圓弧形式已在一定程度上減小了砂輪母線為直線母線形式時的軸向干涉。

圖3 過切干涉

圖4 軸向干涉

2 過切干涉的檢查及消除

根據(jù)引起過切干涉的原因，可以從兩個角度來避免干涉。最容易想到的就是減小砂輪半徑，但是自由曲面的曲率變化有時是很顯著的，要完成整個曲面的磨削加工，有時可用的砂輪半徑就要取得很小。砂輪半徑過小會帶來很多問題，如加工效率將大大降低、表面粗糙度由于殘余高度的加大將受到影響，此外如砂輪半徑過小，則要求主軸轉(zhuǎn)速很高才能滿足切削速度的要求，這就對砂輪主軸的剛度提出了很高的要求，因而提高了加工成本。為此本文提出一種新的方法來消除過切干涉。

2.1 過切干涉檢查

計(jì)算加工曲面各點(diǎn)的曲率半徑，理論計(jì)算方法可參考微分幾何中的曲面第一基本量和第二基本量相關(guān)內(nèi)容[7]，這里不贅述。

本文利用UG軟件通過二次開發(fā)獲取加工曲面各點(diǎn)參數(shù)，通過UF＿M(jìn)ODL＿ask＿curve＿points（）獲取等參數(shù)加工路徑的參數(shù)，包括該等參數(shù)線上各點(diǎn)的曲率半徑，記為rad＿of＿cur；利用UF＿M(jìn)ODL＿ask＿face＿props（）獲得磨削點(diǎn)主曲率，記為k1、k2，法向量記為n，加工點(diǎn)坐標(biāo)值記為Cpoint。

過切干涉僅發(fā)生在凹區(qū)域，因此根據(jù)曲面凹區(qū)域滿足的條件確定需要檢查干涉的磨削點(diǎn)。對于凸區(qū)域則不需要進(jìn)行過切干涉檢查，直接進(jìn)入刀位點(diǎn)進(jìn)行計(jì)算即可。判斷方式如下：

砂輪端面的初始位置與加工軌跡線平行，因此當(dāng)砂輪半徑大于加工軌跡線上磨削點(diǎn)的曲率半徑時會發(fā)生干涉，另外加工軌跡線上其他點(diǎn)的曲率半徑發(fā)生突變時也會引起干涉，如圖5所示。

圖5 過切干涉檢查示意圖

因此，干涉檢查過程如下：首先，對于每一條參數(shù)曲線，取 min｛rad＿of＿cur｝，記為ρmin，并獲取該最小曲率點(diǎn)坐標(biāo)，計(jì)算該坐標(biāo)與磨削點(diǎn)坐標(biāo)沿砂輪端面的距離dt；對于每一個磨削點(diǎn)，如果rad＿of＿cur＜Rgr（Rgr為砂輪的實(shí)際加工曲率半徑）則發(fā)生干涉，如果rad＿of＿cur＞Rgr則繼續(xù)進(jìn)行判斷，如果dt＜R（R為砂輪的最大加工曲率半徑）且min（rad＿of＿cur）＜Rgr則發(fā)生干涉，其他情況則通過干涉檢查。

2.2 過切干涉消除

過切干涉的消除主要基于曲率擬合的條件，通過調(diào)整砂輪端面與加工軌跡面的夾角來消除。

（1）計(jì)算砂輪的實(shí)際加工曲率半徑Rgr。由于本文算法適用的機(jī)床為固定刀軸方向——沿Z軸負(fù)方向，而回轉(zhuǎn)臺帶傾角回轉(zhuǎn)的數(shù)控磨床結(jié)構(gòu)，因此，磨削點(diǎn)法向量與砂輪端面存在夾角，對砂輪半徑有放大作用，如圖6所示。根據(jù)該法向截面上投影橢圓的曲率半徑即可求得砂輪實(shí)際加工半徑：

圖6 砂輪實(shí)際加工半徑

（2）計(jì)算磨削點(diǎn)等參數(shù)方向與最小曲率方向夾角。根據(jù)微分幾何知識，對于自由曲面上任意一點(diǎn)都有兩個主曲率，記為k1、k2（k1＜k2），它們的切向即為兩個主方向。由于曲面是連續(xù)的，故將曲面凹區(qū)域的磨削點(diǎn)定義為凹橢圓點(diǎn)。根據(jù)歐拉定理，曲面上任意點(diǎn)在其切平面上任意方向具有不同的曲率半徑，其范圍在1／k1～1／k2之間連續(xù)變化，如圖7所示，則磨削點(diǎn)等參數(shù)方向的曲率可根據(jù)下式計(jì)算[8]：

圖7 磨削點(diǎn)曲率的連續(xù)變化

式中，θ為kt所在方向與k1所在方向的夾角。

根據(jù)2.1節(jié)獲得的參數(shù)可以求得磨削點(diǎn)等參數(shù)方向的kt與最小曲率方向的夾角：

根據(jù)2.1節(jié)干涉檢查結(jié)果，發(fā)生過切干涉有兩種情況，前者取后者取kt＝分別代入式（2）求得θ。

（3）計(jì)算滿足曲率擬合的最佳磨削截面。設(shè)最佳截面方向的曲率為kfit，則最佳截面方向與最小曲率方向的夾角為

代入式（1）得

因此，等參數(shù)加工軌跡截面與最佳磨削截面的夾角φ＝θ－λ。

這樣就找到了無干涉加工截面。通過回轉(zhuǎn)臺的旋轉(zhuǎn)運(yùn)動令該截面與砂輪端面平行就能消除過切干涉。接下來根據(jù)機(jī)床旋轉(zhuǎn)運(yùn)動的形式反解計(jì)算回轉(zhuǎn)臺旋轉(zhuǎn)量。

3 軸向干涉的檢查與消除

上述過切干涉消除確定了刀軸方向，消除了曲面k2方向上的干涉，但在k1方向上還可能存在干涉。對于圖2所示砂輪結(jié)構(gòu)，由于在k1方向上砂輪圓弧面的曲率半徑一般小于曲面的曲率半徑，所以不考慮過切干涉，只需考慮軸向干涉。軸向干涉是由于加工點(diǎn)法矢量與刀觸點(diǎn)法矢量不一致引起的，如圖8所示。

消除該類干涉的關(guān)鍵就是找到合適的刀觸點(diǎn)。每一個磨削點(diǎn)會對應(yīng)不同的刀觸點(diǎn)，根據(jù)磨削點(diǎn)的法向量求得調(diào)整后的刀觸點(diǎn)，然后通過平移將合適的刀觸點(diǎn)平移到與磨削點(diǎn)重合，即可實(shí)現(xiàn)無干涉的加工，如圖9所示。該方法只需要四軸聯(lián)動即可完成加工。

圖8 法向量不一致引起的軸向干涉

圖9 調(diào)整刀觸點(diǎn)避免軸向干涉

4 葉片曲面磨削仿真實(shí)例

4.1 葉片磨削仿真軟件

葉片加工仿真軟件界面見圖10，主要由刀位點(diǎn)計(jì)算模塊、加工參數(shù)設(shè)置及仿真模塊組成。

圖10 葉片磨削仿真軟件界面

刀位點(diǎn)計(jì)算模塊包含了上述防干涉磨削算法的相關(guān)代碼，由該模塊生成加工軌跡；然后用“生成刀具”按鈕建立砂輪的參數(shù)化模型；加工參數(shù)設(shè)置及磨削仿真模塊完成砂輪模型與葉片表面模型的快速差集運(yùn)算。加工仿真過程如圖11所示，加工過程數(shù)據(jù)（刀軌）如圖12所示。本仿真軟件結(jié)合機(jī)床的運(yùn)動形式，通過砂輪的平移運(yùn)動和葉片繞回轉(zhuǎn)臺固定回轉(zhuǎn)軸的轉(zhuǎn)動，四軸聯(lián)動完成加工過程的仿真。圖12中前三列數(shù)據(jù)顯示的是實(shí)時平移量，最后一列顯示的是實(shí)時回轉(zhuǎn)臺旋轉(zhuǎn)量。

圖11 仿真時砂輪正磨削葉片表面

圖12 葉片磨削加工仿真實(shí)時數(shù)據(jù)截圖

4.2 葉片表面形狀精度分析

利用UG偏差測量模塊分析磨削仿真加工后的葉片表面。磨削仿真采用直徑50mm、圓弧半徑5mm的砂輪，并設(shè)置U向進(jìn)給量為0.5mm，V向進(jìn)給量為1mm。仿真磨削加工后葉片表面的形狀偏差最大為0.0508mm，最小為0.000 25mm，由此可知磨削仿真精度達(dá)到0.05mm，且消除了加工中的過切和干涉。

5 結(jié)論

（1）提出了曲面磨削中過切干涉和軸向干涉的概念及其形式，并給出其各自的干涉消除方法。

（2）提出尋找曲面合適的曲率方向，通過旋轉(zhuǎn)運(yùn)動使其與砂輪端面平行來消除過切干涉的方法，這樣就可以采用較大直徑的砂輪而不需規(guī)定砂輪半徑小于加工曲面的最小半徑，有效地減小了殘留高度，提高了加工表面質(zhì)量，并提高了加工效率。

（3）設(shè)計(jì)了一種圓弧母線圓弧面結(jié)構(gòu)的砂輪，使得整個圓弧面作為工作表面，可以方便地找到能與磨削點(diǎn)法向量一致的刀觸點(diǎn)從而完成無干涉擬合，使得在四軸聯(lián)動數(shù)控磨床上進(jìn)行扭轉(zhuǎn)葉片的加工成為可能。

[1]Wu Fuzhong.Five－axis Tool Path Planning by Free Form Curve Offset in Surface[C]／／Proceedings of the IEEE International Conference on Automation and Logistics.Qingdao，2008：559-563.

[2]Lin Than，Lee Jae－Woo，Bohez E L J.A New Accurate Curvature Matching and Optimal Tool Based Five－axis Machining Algorithm[J].Journal of Mechanical Science and Technology，2009，23：2624-2634.

[3]Chen Tao，Ye Peiqing，Wang Jinsong.Local Interference Detection and Avoidance in Five－axis NC Machining of Sculptured Surfaces[J].International Journeral Advantage Manufacture Technology，2005，25：343-349.

[4]石璟.面向葉片加工的數(shù)控砂帶拋磨系統(tǒng)的研究[D].無錫：江南大學(xué)，2009.

[5]Zhao Yang，Zhao Ji，Zhang Lei.Path Planning for Automatic Robotic Blade Grinding[C]／／International Conference on Mechatronics and Automation.Changchun，2009：1556-1560.

[6]Choi B O，Lee M K.Development of a Hybrid Robot Arm for Propeller Grinding[C]／／Proceedings of SPIE－The International Society for Optical Engineering.Intelligent Robots and Computer Vision XV：Algorithms，Techniques，Active Vision，and Materials Handling.Boston，1996，2904：517-527.

[7]陳維恒.微分幾何[M].北京：北京大學(xué)出版社，2006.

[8]蘇高峰，張秋菊.基于UG的葉片數(shù)控砂帶磨削刀位數(shù)據(jù)的計(jì)算[J].機(jī)床與液壓，2006（1）：61-63.