于洋，鄧勇生，潘琳

(1.西安石油大學(xué)機(jī)械工程學(xué)院，陜西西安710065; 2.中國航天科技集團(tuán)公司第四研究院第四十四研究所，陜西西安710025)

自由曲面廣泛地應(yīng)用于現(xiàn)代工業(yè)設(shè)計(jì)中，例如航空、航天、汽車、模具等。5軸數(shù)控加工給復(fù)雜自由曲面的加工帶來了更多的靈活性、高效率以及良好的精度表現(xiàn)，可以處理三軸機(jī)床無法加工的曲面［1-2］。然而，5軸數(shù)控加工仍有許多問題有待解決和改進(jìn)，比如刀軸控制就被認(rèn)為是數(shù)控編程中最困難的問題之一，主要是由于5軸加工中刀具運(yùn)動的復(fù)雜性和復(fù)雜自由曲面曲率分布的不規(guī)則所致。在5軸加工中，由于刀具方向的不斷變化，當(dāng)?shù)毒咻S向過度依賴被加工曲面幾何時(shí)，容易產(chǎn)生刀具軸向變化劇烈，從而導(dǎo)致過切、碰撞以及切削誤差過大的情況。因此，要充分發(fā)揮5軸數(shù)控加工的優(yōu)勢，必須解決刀軸過度依賴曲面幾何導(dǎo)致變化劇烈這一瓶頸問題。

近年來，國內(nèi)外學(xué)者對此做了許多工作，然而大多數(shù)都集中在處理避免過切干涉的發(fā)生及其修正上［3-6］。對于加工精度的要求，通常都是通過增加刀具接觸點(diǎn)的方式達(dá)成，雖然的確能夠提高加工精度，但是會造成NC數(shù)據(jù)過于龐大的問題，很少有通過調(diào)整刀具軸向來達(dá)到避免過切、提高精度的文獻(xiàn)。MORISHIGE等［7］應(yīng)用2-Dimensional C-Space的方法，將傾斜角 (Lead Angle)及偏擺角 (Tilt Angle)映射至C-Space上，求取合理刀具軸向，但對于自由曲面，須將曲面離散成細(xì)小平面再做計(jì)算，效率很低。文獻(xiàn)［8-9］通過尋找可行性的刀具方向來控制曲面加工誤差，效率較低。文獻(xiàn) ［10］提出了通過設(shè)置權(quán)重的方式來進(jìn)行刀具軸向的均勻化分配，但對其中的干涉檢查和處理分析不透，并沒有從理論和方法上進(jìn)行創(chuàng)新和突破。基于此，作者提出了一種新的刀具軸向設(shè)計(jì)算法，該方法根據(jù)曲面的微分幾何特性，通過制定可能發(fā)生過切處的刀具軸向，藉以計(jì)算其他位置的軸向，同時(shí)進(jìn)行過切、干涉檢查及處理，提高切削效率，保證曲面的加工精度和加工質(zhì)量。

1 刀具軸向設(shè)計(jì)考量及初始位置

1.1 刀具軸向設(shè)計(jì)考量

考慮到5軸加工中，設(shè)計(jì)不合理的刀軸方向會導(dǎo)致機(jī)床震動、切削誤差過大、過切、碰撞以及加工不完整的情況發(fā)生，因此，設(shè)計(jì)刀具軸向時(shí)必須要滿足以下條件:(1)刀具軸向的設(shè)計(jì)必須考慮過切、碰撞的情況;(2)同一條刀具路徑上，兩個(gè)刀具接觸點(diǎn)之間的刀具軸向變化量應(yīng)該盡可能小。因此，作者提出了通過指定可能發(fā)生過切、干涉區(qū)域的刀具軸向，再依據(jù)四元數(shù)插補(bǔ)方法，計(jì)算出其余位置刀具軸向的方法。

1.2 刀具初始位置

設(shè)定當(dāng)前刀具接觸點(diǎn) (CC點(diǎn))坐標(biāo)為C(xC，yC，zC)，過點(diǎn)C的曲面切平面可表示如下:

其中［A，B，C］為切平面的單位法矢，也是刀具軸向的單位方向矢量，可表示如下:

其中:P(u，v)為曲面上任一點(diǎn)的矢量，u、v為變量，則曲面上任意一點(diǎn)可由直角坐標(biāo)系的分量x(u，v)、y(u，v)、z(u，v)來表示。

在重要區(qū)域的刀具軸向的指定方法上，為了幾何上計(jì)算簡便，并且可以進(jìn)行嚴(yán)謹(jǐn)?shù)牡毒吲鲎驳臋z測，作者提出以曲面上CC點(diǎn)處的法線方向作為基準(zhǔn)，在刀具前進(jìn)方向，傾斜某一角度θ來得到刀具初始位置，如圖1所示。

圖1 初始刀具軸向

2 四元數(shù)內(nèi)插法

四元數(shù) (Quaternion)最初的理論是以一對三維空間向量之間的角度關(guān)系定義，近年來則常被應(yīng)用在電腦圖形學(xué)及機(jī)器人學(xué)上，用來描述運(yùn)動方式。

一個(gè)四元數(shù)q中包含4個(gè)元素，定義如下:

其中:qx，qy，qz，qw為實(shí)數(shù)，i，j，k則是符合以下規(guī)則的坐標(biāo)軸:

qw稱為純量部分，而(qx，qy，qz)稱為向量部分。

圖2 四元數(shù)插補(bǔ)法

四元數(shù)算法不僅支持相同的旋轉(zhuǎn)速度，而且在大圓弧上的插補(bǔ)向量連接兩個(gè)邊界向量，如圖2所示。給定兩個(gè)指定的刀軸方向V1(V1x，V1y，V1z)、V2(V2x，V2y，V2z)，對于從V1(V1x，V1y，V1z)、V2(V2x，V2y，V2z)的四元插補(bǔ)方程V(n)可以表示如下:

其中:參數(shù)n為插V1與V2之間的向量，θ= cos-1(V1·V2)

使用四元數(shù)插補(bǔ)法對圖3中所指定的刀具軸向進(jìn)行其余刀具軸向插補(bǔ)，插補(bǔ)后的刀具軸向分布如圖4所示，圖5顯示了刀具軸向變化角度的比較圖，可以看出:四元數(shù)插補(bǔ)法可以提供光順的刀軸變化，并且不會產(chǎn)生不需要的刀軸扭曲和跳動。

圖3 指定特定位置刀具軸向

圖4 四元插補(bǔ)法所生成的刀具軸向

圖5 四元數(shù)插補(bǔ)法刀具軸向變化

3 干涉刀位的修正

在完成刀具軸向設(shè)計(jì)之后，則需要對每一個(gè)刀具軸向來進(jìn)行刀具干涉?zhèn)蓽y。就空間幾何來看，直接計(jì)算檢查點(diǎn)到刀具的垂足距離d并判斷，若其垂足距離d小于刀具半徑D/2，則表示發(fā)生刀具碰撞;反之，若其垂足距離d大于刀具半徑D/2，則無碰撞。

圖6 碰撞示意及參數(shù)定義

假設(shè)刀具軸向的旋轉(zhuǎn)中心為P1(x1，y1，z1)，端點(diǎn)為P2(x2，y2，z2)，而檢查點(diǎn)為PZ(x3，y3，z3)，如圖6所示，可以用參數(shù)a表示成參數(shù)式:

令Δx1=x2-x1、Δy1=y2-y1、Δz1=z2-z1得:

接下來計(jì)算檢查點(diǎn)到刀具中心的垂足點(diǎn)P(x(a)，y(a)，z(a))的距離d:

將式(8)代入式(7)得:

令Δx2=x3-x1、Δy2=y3-y1、Δz2=z3-z1，式(9)可如下表達(dá):

又點(diǎn)PZ到刀具中心軸的垂足距離為最短距離，所以:

最后求得的a為:

將所得的a代入式 (7)求檢查點(diǎn)PZ到刀具中心軸的垂足點(diǎn)P，而PZ到刀具中心軸的垂足距離d亦可求得。最后判斷若垂足距離d小于刀具半徑D/2，則發(fā)生刀具碰撞亦檢查點(diǎn)PZ為碰撞點(diǎn);若垂足距離d大于刀具半徑D/2，則檢查點(diǎn)PZ不是碰撞點(diǎn)。

4 算例分析

為了更好地驗(yàn)證文中所提出的刀軸設(shè)計(jì)算法，作者使用主軸旋轉(zhuǎn)式5軸機(jī)床進(jìn)行切削實(shí)驗(yàn)。實(shí)驗(yàn)將針對兩種不同的刀具軸向設(shè)計(jì)方式:傳統(tǒng)法向法與四元數(shù)插補(bǔ)法，對比兩種加工方法的切削誤差，并進(jìn)行探討。加工條件如表1所示，曲面模型如圖7所示，采用等參數(shù)加工法，機(jī)床型號為Mazak5軸立式加工中心500-5X。

表1 加工條件

圖8為使用四元數(shù)插補(bǔ)法所設(shè)計(jì)的刀具軸向分布圖，圖9為使用傳統(tǒng)的法向刀具軸向設(shè)計(jì)方法所生成的刀具軸向?？梢悦黠@看出:使用傳統(tǒng)方法所設(shè)計(jì)的刀軸變化較為劇烈，刀軸角度平均變化量大概在25°左右，最大的變化量達(dá)到了58°;使用四元數(shù)插補(bǔ)法設(shè)計(jì)的刀軸角度變化不超過4°，很適合用來進(jìn)行刀具軸向規(guī)劃。

圖7 實(shí)驗(yàn)所使用模型

圖8 四元數(shù)法設(shè)計(jì)刀具軸向

圖9 法向法所生成刀具軸向

圖10顯示了采用法向刀軸方法加工輪廓時(shí)切削誤差變化圖;圖11顯示了采用四元數(shù)刀軸方法加工輪廓時(shí)切削誤差變化圖?？梢钥闯?法線加工方法誤差較大，已經(jīng)超過設(shè)定誤差0.02 mm，加工效果不理想，這是由于刀具軸向變化較為劇烈所引起;而四元數(shù)法對誤差控制較好，均小于設(shè)定誤差emax=0.02 mm，較傳統(tǒng)方法有大幅度的改進(jìn)。

圖10 法向刀軸切削誤差

圖11 四元數(shù)插補(bǔ)刀軸切削誤差

5 結(jié)論

通過使用四元數(shù)插補(bǔ)算法，解決了5軸自由曲面加工中刀軸變化劇烈的問題，并且根據(jù)刀具半徑與檢查點(diǎn)到刀具的垂足距離判斷刀具是否與曲面發(fā)生干涉，成功地解決了刀軸光順與曲面干涉的檢測與修正這一5軸數(shù)控加工中的瓶頸問題。算法的運(yùn)用結(jié)果表明:該算法效率高，精度易于控制，可用于葉輪類對刀軸變化敏感且干涉較為嚴(yán)重的5軸加工中。

【1】LEE Y S，CHANG T C.2-Phase Approach to Global Tool Interference Avoidance in 5-axis Machining［J］.Computer-Aided Design，1995，27(10):715-729.

【2】YOU C F，CHU C H.Tool-path Veri Cation in Ve-axis Machining of Sculptured Surfaces［J］.The International Journal of Advanced Manufacturing Technology，1997，13(4): 248-255.

【3】HWANG Y R，HO M C，HU C H.Five-axis Tool Path Generation Based on Discrete Point Data［C］//Proc Precision Machinery Manufacturing Conference，Taipei，2001:149-153.

【4】HO M C，HWANG Y R.A New Decision Algorithm of Maximum Allowable Step Length for 5-axis Table-tilting Type Machining［J］.Journal of the Chinese Institute of Engineers，2002，25(2):233-242.

【5】BALASUBRAMANIAM M，HO S，SARMA S，et al.Generation of Collision-free 5-axis Tool Paths Using a Haptic Surface［J］.Computer-Aided Design，2002，34:267-279.

【6】張和明，張玉云，熊光楞.復(fù)雜曲面5坐標(biāo)數(shù)控加工干涉檢查及刀位修正［J］.清華大學(xué)學(xué)報(bào):自然科學(xué)版，1998 (2):65-68.

【7】MORISHIGE K，KASE K，TAKEUCHI Y.Tool Path Generation Using C-Space for 5-Axis Control Machining［J］.Trans of the ASME，J of Manufacturing Science and Engineering，1999，121(1):144-149.

【8】LEE Y S.Mathematical Modeling Using Different End Mills and Tool Placement Problems for 4-and 5-axis NC Complex Surface Machining［J］.International Journal of Production Research，1998，36(3):785-814.

【9】LEE Y S.Non-isoparametric Tool Path Planning by Machining Strip Evaluation for 5-axis Sculptured Surface Machining［J］.Computer-Aided Design，1998，30(7):559-570.