黃象珊

(浙江經(jīng)濟(jì)職業(yè)技術(shù)學(xué)院 汽車技術(shù)學(xué)院，杭州 310018)

?

基于適應(yīng)性空間填充曲線生成刀具路徑的技術(shù)研究

黃象珊

(浙江經(jīng)濟(jì)職業(yè)技術(shù)學(xué)院 汽車技術(shù)學(xué)院，杭州 310018)

空間填充曲線還沒有被廣泛應(yīng)用于五軸數(shù)控系統(tǒng)中，因?yàn)樵诿枋鰳?biāo)準(zhǔn)空間填充曲線樣式時(shí)，由于尖角轉(zhuǎn)向會(huì)產(chǎn)生很大的誤差。為了消除較大運(yùn)動(dòng)誤差和尖角轉(zhuǎn)向的過(guò)度切削，將適應(yīng)性空間填充曲線生成思想融入路徑生成過(guò)程中，提出基于適應(yīng)性的空間填充曲線刀具路徑生成方法，在矩形網(wǎng)格上以最短Hamiltonian軌跡算法為指導(dǎo)來(lái)生成填充曲線，達(dá)到優(yōu)化適應(yīng)性填充曲線刀具路徑的目的，通過(guò)實(shí)例驗(yàn)證了方法的可行性。

刀具路徑；適應(yīng)性空間填充曲線；Hamiltonian

0 引言

應(yīng)用最為廣泛的空間填充曲線(SFC)需要考慮很多因素，包括刀具軌跡設(shè)計(jì)的遞歸Hilbert曲線。Hilbert曲線特別具有吸引力的是刀具軌跡設(shè)計(jì)，因?yàn)樗木?xì)性可以被用于適當(dāng)?shù)卦黾勇窂降拿芏?。然而，每條基于Hilbert曲線的刀具軌跡的精細(xì)性都在精確區(qū)域內(nèi)增加刀具軌跡密度，由于增加了刀具軌跡的長(zhǎng)度，從而導(dǎo)致了較低的加工效率。另外，Hilbert曲線還有一個(gè)容易使刀具不斷改變方向的特點(diǎn)，這就使得加工過(guò)程緩慢而且會(huì)引起很大的運(yùn)動(dòng)誤差。

為了彌補(bǔ)這些缺陷，本文提出利用適應(yīng)性空間填充曲線生成刀具路徑的方法。適應(yīng)性空間填充曲線被用于生成刀具路徑。對(duì)其特點(diǎn)描述如下：首先，適應(yīng)性空間填充曲線總是遵循局部?jī)?yōu)化方向；第二，相對(duì)于常規(guī)的SFC，適應(yīng)性空間填充曲線只在需要時(shí)轉(zhuǎn)向，也就是當(dāng)優(yōu)化方向改變的時(shí)候；第三，適應(yīng)性空間填充曲線消除了較大運(yùn)動(dòng)誤差的出現(xiàn)和由于尖角轉(zhuǎn)向而引起的過(guò)度切削；最后，適應(yīng)性空間填充曲線的局部精細(xì)化處理以與常規(guī)SFC的實(shí)現(xiàn)方式一樣。

適應(yīng)性空間填充曲線生成的一個(gè)問(wèn)題是以Hamiltonian路徑關(guān)于在(u，v)平面中包括CC點(diǎn)在內(nèi)的網(wǎng)格圖表為例進(jìn)行精確描述的。運(yùn)算法則包括三步：構(gòu)造網(wǎng)格、空間填充曲線的生成和更正刀具路徑。

1 空間填充曲線刀具路徑

1.1 填充曲線的描述

空間填充曲線[1-2]起源于1878年，George Cantor(1845—1918)證明了空間任意有限維的任何兩個(gè)光滑流形具有相同的基數(shù)。Cantor的發(fā)現(xiàn)表明了單位線段[0，1]可以被雙射到單位正方形[0，1]2上。然而1879年，Eugen Netto(1848—1919)證明了這種映射必須是不連續(xù)的，不能稱之為曲線。如果不滿足這個(gè)雙射條件，Giuseppe Peano(1858—1932)在1890年發(fā)現(xiàn)了一個(gè)區(qū)間到平面的映射。這是空間填充曲線的第一例子，如圖1所示；后來(lái)更多的例子被Hilbert(1891年，如圖2所示)。

圖1 Peano空間填充曲線的3次迭代

圖2 Hilbert空間曲線的6次迭代

1.2 空間填充曲線刀具路徑生成

當(dāng)確定好填充曲線的類型后，根據(jù)加工質(zhì)量的要求，選擇合適的填充曲線的階次，經(jīng)過(guò)干涉處理后 ，便可作為加工刀具的走刀路徑。

2 空間填充曲線刀具路徑生成

2.1 網(wǎng)絡(luò)構(gòu)造

圖3所示為由兩個(gè)等參數(shù)刀具路徑疊加而得的網(wǎng)格[3-4]m×n。等參數(shù)路徑的生成是通過(guò)計(jì)算最小刀具軌跡間隔并將其作為下一個(gè)刀具路徑的連續(xù)偏移。刀具路徑間隔依賴于局部表面形狀、切削刀具的形狀和尺寸以及允許的皺褶高度。此外，兩個(gè)相鄰切削刀具路徑的加工帶寬度必須重疊，以保證加工表面的誤差(皺褶高度)在公差范圍內(nèi)。此外，網(wǎng)格被看作是一個(gè)每?jī)蓚€(gè)相鄰頂點(diǎn)都用邊緣線連接起來(lái)的無(wú)向圖，如圖3d所示。圖的頂點(diǎn)相當(dāng)于最初的一系列在必需平面上的CC點(diǎn)，而兩個(gè)相連接頂點(diǎn)間的距離就是在三維空間R3內(nèi)對(duì)應(yīng)的CC點(diǎn)之間的距離。注意切削刀具沿著在任意兩個(gè)相連接頂點(diǎn)之間的路徑滿足皺褶深度約束。這個(gè)特征允許依靠空間填充曲線來(lái)優(yōu)化刀具路徑。

(a)在V方向上的等參數(shù)刀具路徑(b)在u方向上的等參數(shù)刀具路徑(c)兩個(gè)等參數(shù)刀具路徑的迭代(d)最終的網(wǎng)格圖3 兩個(gè)等參數(shù)刀具路徑疊加而得的網(wǎng)絡(luò)

2.2 空間填充曲線刀具路徑的生成

2.2.1 Hamiltonian回路形成原因

網(wǎng)格構(gòu)造被看作是Hamiltonian路徑問(wèn)題來(lái)確切描述的。找一條距離很短的路徑只不過(guò)是個(gè)旅問(wèn)題。既然這是個(gè)NP-hard問(wèn)題，尋找最優(yōu)化方案的算法求解較慢且效率低。Hamiltonian軌跡運(yùn)算是一種簡(jiǎn)單且計(jì)算效率高的法則，該算法是基于由Dafner等為了二維圖像掃描而開發(fā)的覆蓋合并算法，Hamiltonian軌跡算法被應(yīng)用擴(kuò)展到任意矩形網(wǎng)格中[5]。

2.2.2 矩形網(wǎng)格上Hamiltonian路徑的求解算法

算法思路：首先，不相交的回路經(jīng)過(guò)所有頂點(diǎn)，然后回路被合并為單獨(dú)的Hamiltonian回路。最初的回路是由小矩形環(huán)狀軌跡通過(guò)相鄰4個(gè)頂點(diǎn)連接起來(lái)的，然后將偶數(shù)行和列的頂點(diǎn)和奇數(shù)行和列上的頂點(diǎn)連接起來(lái)，如圖4a所示。另外，如果m或n是奇數(shù)，則將頂點(diǎn)沿著邊界用短劃線連接成虛線回路，如圖4b所示。相鄰回路可連接一個(gè)較大的回路。合并的代價(jià)按照下式計(jì)算。

Cost(A，B)=∣s∣+∣t∣-∣e∣-∣f∣

其中，∣e∣表示在R3中由邊緣線e連接兩個(gè)頂點(diǎn)間的距離。

合并兩個(gè)虛線回路的代價(jià)設(shè)為-∞，即所有的虛線回路開始時(shí)都被連接起來(lái)。這樣就保證了刀具路徑從Hamiltonian軌跡移除虛線邊后依然是連續(xù)的。而且，實(shí)線回路A可與一個(gè)虛線回路D合并，僅當(dāng)A與位于相反方向的實(shí)線回路C合并時(shí)。為了實(shí)現(xiàn)這個(gè)合并屬性，將A與D的合并代價(jià)等于A與C的合并代價(jià)。這個(gè)合并屬性被用于減少沿著邊界有大量拐點(diǎn)的不合適的狹窄的之字形刀具路徑生成的可能性。

對(duì)偶圖G′首先被構(gòu)造出來(lái)以合并所有小回路，如圖4b所示。G中的每個(gè)小回路在G′中都定義了一個(gè)頂點(diǎn)，而且G中連接A與B兩個(gè)小回路的兩條邊s和t在G′中也定義了一個(gè)邊v′。

(a)被初始回路覆蓋的無(wú)向圖G (b)對(duì)偶圖G′

如圖5所示除了構(gòu)造對(duì)應(yīng)的最小生成樹之外，還闡明了合并算法[6-7]。在所有回路都被合并為Hamiltonian回路后，通過(guò)移除所有回路邊緣而生成了刀具軌跡，具體如圖5c所示。令T為一系列邊緣最小生成樹，合并算法如下。

1)設(shè)T=φ，將G′的邊按合并代價(jià)以升序排序。如果出現(xiàn)一個(gè)結(jié)，將連接兩個(gè)頂點(diǎn)的邊排在前面。

2)考慮第一個(gè)未通過(guò)的邊線。將其包括在T中，當(dāng)且僅當(dāng)它其他T的邊在對(duì)偶圖G′中沒形成回路，且也不違背合并屬性。

3)若邊放入T中，合并G中的兩個(gè)回路，其相當(dāng)于用新增加的邊連接G′中的兩個(gè)頂點(diǎn)，然后進(jìn)入第4步，否則返回第2步。

4)如果T包括n-1條邊(其中n代表G′中頂點(diǎn)的數(shù)量)，停止并輸出Hamiltonian回路的結(jié)果，否則返回到第2步。

2.2.3 刀具路徑更正

按照上面的步驟生成的刀具路徑還有兩方面需要修改[8-9]：第一，刀具路徑的軌跡應(yīng)該被修改以消除底切；第二，當(dāng)?shù)毒吒淖兎较驎r(shí)，刀具方向也要作仔細(xì)調(diào)整。

刀具路徑的調(diào)整是必要的，因?yàn)楫?dāng)?shù)毒吒淖兎较驎r(shí)，空間填充曲線刀具路徑上的轉(zhuǎn)向會(huì)致使刀具錯(cuò)過(guò)某些加工表面區(qū)域。在每一個(gè)轉(zhuǎn)向處，兩條相鄰刀具路徑沿不同方向切削形成的切削帶不會(huì)重疊而形成底切，或重疊部分不足而產(chǎn)生了很明顯的褶皺(如圖6所示)。

(a)通過(guò)小回路合并構(gòu)造的Hamiltonian回路(b)最小生成樹 (c)生成的刀具軌跡圖5 刀具軌跡的生成

圖6 在相鄰刀具路徑上通過(guò)使用空間填充曲線生成的加工帶(短畫線)

3 舉例應(yīng)用

一個(gè)曲面控制頂點(diǎn)給出，要求加工表面公差h=0.01mm ，用半徑R=3mm 的平頭刀進(jìn)行加工，應(yīng)用上述算法實(shí)現(xiàn)的刀具路徑結(jié)果。

例： 如圖7所示，一個(gè)三次Bezier表面，既有凸形區(qū)域又有凹形區(qū)域。這個(gè)表面被描述為：

(a)例子中表面在(u，v)區(qū)域兩個(gè)等參數(shù)刀具路徑的重疊

其中，Px，Py，Pz是16個(gè)控制點(diǎn)的x，y，z坐標(biāo)。利用重疊兩個(gè)等參數(shù)軌跡構(gòu)造的網(wǎng)格如圖7所示。刀具軌跡的結(jié)果如圖8和圖9所示。刀具軌跡更正用于切削加工前、后的加工結(jié)果，如圖10和圖11所示。很明顯，在保證表面質(zhì)量的同時(shí)，這種方法能將未被移除的多余材料去除。表1所示為就刀具軌跡長(zhǎng)度而言，空間填充曲線(SFC)刀具路徑對(duì)等參數(shù)刀具路徑的情況。

(a)例子中表面在(u，v)區(qū)域中的SFC刀具路徑(b)例子中的表面在工件坐標(biāo)系中的SFC刀具路徑圖8 對(duì)于例子中的表面在(u，v)區(qū)域和工件坐標(biāo)系中的SFC刀具路徑

圖9 對(duì)于例子中刀具路徑的合并和不合并的情況(a)例子中刀具路徑的合并情況(b)例子中刀具路徑的不合并情況

刀具路徑刀具軌跡長(zhǎng)度/mm例在v方向上的等參數(shù)3917．31在u方向上的等參數(shù)2648．12SFC刀具路徑2637．54

總之例如加工一個(gè)像以上例子中的簡(jiǎn)單表面，不需要很多帶有銳角變化的拐點(diǎn)。因此，相應(yīng)的SFC刀具軌跡無(wú)論是從刀具軌跡長(zhǎng)度方面還是加工時(shí)間方面都比常規(guī)的等參數(shù)刀具軌跡要好。但是對(duì)于復(fù)雜形狀的表面，當(dāng)角速度未被約束時(shí)，尖角拐點(diǎn)會(huì)減緩加工進(jìn)程。對(duì)于這樣的表面，SFC方法更適用于有高速旋轉(zhuǎn)軸的銑床上。

圖10 應(yīng)用沒修正過(guò)的SFC刀具路徑的實(shí)際加工 圖11 應(yīng)用修正過(guò)的SFC刀具路徑的實(shí)際加工

4 結(jié)論

在參數(shù)線法刀具路徑生成的基礎(chǔ)上結(jié)合適應(yīng)性空間填充曲線生成思想，提出適應(yīng)性填充曲線刀具路徑成方法[10]；以加工要求為指導(dǎo)生成，兩方向上的參數(shù)線路徑，然后疊加形成網(wǎng)格，在Hamiltonian路徑生成算法的引導(dǎo)下，完成網(wǎng)格的依次連接，最后在網(wǎng)上自然的呈現(xiàn)出適應(yīng)性空間填充曲線路徑，該路徑即為適應(yīng)性空間填充曲線刀具路徑；通過(guò)實(shí)例應(yīng)用，證實(shí)所提出方法的可行性及設(shè)計(jì)算法的正確性。

[1] 趙玉剛，李偉，郭峰,等. 基于Cox-de Boor 遞推的任意次NURBS曲線插補(bǔ)算法的研究與仿真[J].組合機(jī)床與自動(dòng)化加工技術(shù)，2012(6):45-49.

[2] Sagan,H.(1994). Space-Filling Curves. Springer-Verlag, New York.

[3] 何晶晶，夏銘，黃志勇. 三次樣條曲線回轉(zhuǎn)的車削加工[J].組合機(jī)床與自動(dòng)化加工技術(shù).2011(4):106-109.

[4] 周永情，李德明，孫軍偉. 自由曲線輪廓數(shù)控加工路徑生成及加工[J]. 組合機(jī)床與自動(dòng)化加工技術(shù).2011(4):94-97.

[5] 馬哈諾夫，安諾泰派布.五軸數(shù)控銑床切削加工先進(jìn)數(shù)值優(yōu)化方法[M].北京：機(jī)械工業(yè)出版社，2013.

[6] 淡卜綢，李德信.基于Hilbert填充曲線的自由曲面刀具路徑規(guī)劃研究[J]. 機(jī)械設(shè)計(jì)與制造,2010(12)：233-236

[7] 周濟(jì)，周艷紅.數(shù)控加工技術(shù)[M]. 北京：國(guó)防工業(yè)出版社，2002.

[8] 張文博. 基于填充曲線的曲面數(shù)控刀具軌跡自動(dòng)生成算法研究[D].長(zhǎng)春：長(zhǎng)春理工大學(xué)，2006.

[9] 叢萌.數(shù)控刀具軌跡自動(dòng)生成技術(shù)的研究[D]長(zhǎng)春：長(zhǎng)春理工大學(xué)，2004.

[10] 淡卜綢，趙亞寧，魯開講. 基于參數(shù)線法的填充曲線刀具路徑生成方法研究[J].機(jī)械設(shè)計(jì)與制造,2010(12)：230-232.

(編輯 李秀敏)

Research of Technology Tool Path Based on the Adaptability Space Filling Curve Comes into being Cutter Route′s

HUANG Xiang-shan

(College of Automotive Technical ，Zhejiang Technology Institute of Economy，Hangzhou 310018，China)

Space- filling curve， not having a quilt broad apply to among the five scrolls of numerical control system, when because of the form filling up a curve in the space describing a standard since that's the point angle gets lost may produce very big error. For removing excessive cutting that motion error and the point angle change to more, the curve come into being thought adaptability space is filled up melts to enter route come into being process middle , suggests that method comes to generate the filling curve being guided by the shortest Hamiltonian trajectory algorithm owing to that the adaptability space fills up the curve cutter route come into being on rectangle net lattice , reaches the purpose optimizing adaptability filling up curve cutter route , the feasibility by the fact that the example has verified method .

the tool path；adaptability space- filling curve；hamiltonian

1001-2265(2014)01-0053-04

10.13462/j.cnki.mmtamt.2014.01.015

2013-04-21;

2013-05-30

黃象珊(1976—)，男，浙江溫州人，浙江經(jīng)濟(jì)職業(yè)技術(shù)學(xué)院講師，工學(xué)碩士，主要從事機(jī)電一體化方面教學(xué)與研究，(E-mail)hxshan98@163.com。

TH162;TG71

A