朱珮旗, 樊紅日, 錢波, 茅健
(上海工程技術(shù)大學(xué)機(jī)械與汽車工程學(xué)院, 上海 201620)
熔 融 沉 積 制 造 ( Fused Deposition Modeling,F(xiàn)DM) 是熱塑性聚合物最常用的增材制造技術(shù)之一。在FDM 工藝中, 將復(fù)合材料送入加熱的噴頭中熔化或液化, 然后擠出并沉積在構(gòu)建模型的基板上, 沉積層固化并與相鄰層粘合/焊接, 形成所需的3D 幾何形狀。 該方法最大優(yōu)點(diǎn)是成形效率高、 零件制造周期短, 可以獲得高性能的成形產(chǎn)品, 具有十分廣闊的應(yīng)用前景[1-3]。 FDM 經(jīng)過(guò)多年的發(fā)展, 已應(yīng)用于航空航天、 醫(yī)療、 建筑等眾多領(lǐng)域。
軌跡規(guī)劃在FDM 中非常關(guān)鍵, 軌跡規(guī)劃策略的優(yōu)劣可能會(huì)影響產(chǎn)品的表面粗糙度、 尺寸精度及使用性能等。 此外, 不同的路徑意味著相應(yīng)打印頭的移動(dòng)策略不同, 進(jìn)而影響零件打印的效率[4]。 因此, 更好的軌跡規(guī)劃策略可以獲得更優(yōu)的制造性能、 提高零件制造質(zhì)量和效率。
目前主要的軌跡成形算法有: (1) 光柵填充路徑。 該算法是一種結(jié)構(gòu)簡(jiǎn)單、 易于實(shí)現(xiàn)的基礎(chǔ)算法,但其制造精度不高, 具有大量空行程以及成形效率低, 尤其對(duì)于復(fù)雜曲面構(gòu)件, 該算法很難滿足制造要求[5]。 (2) Zigzag 填充路徑。 該算法首先對(duì)整體切片進(jìn)行掃描, 然后從起始點(diǎn)開始采用往復(fù)的平行直線等距分割二維水平輪廓。 這種填充方式的優(yōu)點(diǎn)是規(guī)劃過(guò)程簡(jiǎn)單且經(jīng)過(guò)的空行程較少, 但路徑的邊緣處容易形成階梯效應(yīng), 近距離的正反向往返易因材料密集而導(dǎo)致散熱減緩、 應(yīng)力集中, 從而導(dǎo)致產(chǎn)品出現(xiàn)變形、 裂紋等現(xiàn)象[6]。 Zigzag 填充路徑算法是目前應(yīng)用最多的填充方式之一。 (3) 輪廓偏置路徑。 該算法基于切片的外輪廓和內(nèi)輪廓, 以固定距離向模型的內(nèi)部偏置,從而產(chǎn)生填充路徑。 因?yàn)樘畛渎窂骄鶠榉忾]的多邊形, 不會(huì)出現(xiàn)穿越截面輪廓的情況, 縮短了空行程,提高了加工效率。 但當(dāng)切片輪廓比較復(fù)雜時(shí), 可能會(huì)造成內(nèi)輪廓或外輪廓的自相交或內(nèi)輪廓和外輪廓相互交錯(cuò)的情況, 極大地增加了偏置輪廓計(jì)算的復(fù)雜度[7]。 (4) 螺旋偏置路徑。 該算法一般先提取零件輪廓骨架線, 再得到偏置軌跡, 會(huì)出現(xiàn)填充路徑自相交或者互相交情況, 且現(xiàn)有螺旋加工路徑生成方法均存在路徑不均勻問(wèn)題, 尤其是對(duì)偏移路徑處理生成螺旋線過(guò)程復(fù)雜, 效率低下, 難以保證實(shí)際加工質(zhì)量[8]。 (5) 分形線填充路徑。 該算法利用分形曲線局部與整體相似的原則不斷進(jìn)行迭代、 重復(fù)嵌套,從而生成可以充滿整個(gè)截面輪廓的填充路徑。 由于分形曲線局部與整體相似, 所以截面輪廓內(nèi)各處的填充路徑規(guī)律相同、 密度均勻, 打印出的分層面平整, 材料分布均勻, 提高了整體的打印精度和性能, 目前較為常用的分形填充曲線為Hilbert 曲線。 分形填充曲線雖然在垂直方向上的力學(xué)性能較好, 但它對(duì)于模型側(cè)向的支撐較弱, 在實(shí)際打印過(guò)程中易產(chǎn)生變形的情況; 且填充曲線由短線段組成, 需要打印頭頻繁地變化移動(dòng)的方向, 導(dǎo)致打印機(jī)使用壽命降低[9]。
復(fù)合材料熔絲成形增材制造技術(shù)對(duì)成形零部件的成形質(zhì)量以及成形效率的要求越來(lái)越高, 研究人員提出了各種各樣的新型路徑規(guī)劃方法。 如YANG 等[10]提出了基于模型單層輪廓的等距偏置路徑規(guī)劃, 切片輪廓向內(nèi)均勻偏置形成填充路徑。 呂鵬輝、 張起貴[11]提出了一種基于Voronoi 圖的多邊形路徑填充算法, 能夠直接獲取三維模型切片輪廓的漸進(jìn)填充掃描線。 DING等[12]提出了一種減少內(nèi)部空隙的中軸轉(zhuǎn)換(Medial Axis Transform, MAT) 路徑規(guī)劃方法。 李浩等人[13]提出了一種固定可變寬度路徑規(guī)劃方法, 通過(guò)在不同區(qū)域設(shè)定不同的沉積寬度, 達(dá)到更高的填充率。
單一的軌跡填充算法各有其優(yōu)點(diǎn), 但不具有普適性, 且很難保證同時(shí)滿足成形質(zhì)量與制造效率的要求。 本文作者針對(duì)單一軌跡填充算法局限性問(wèn)題, 提出一種復(fù)合軌跡成形算法, 將輪廓偏置路徑與Zigzag掃描線填充路徑相結(jié)合, 并在拐角處使用圓弧過(guò)渡對(duì)軌跡進(jìn)行優(yōu)化, 以減小高速打印時(shí)加減速過(guò)程造成的慣性沖擊, 在保證零件成形質(zhì)量的同時(shí), 進(jìn)一步提高制造效率。 建立復(fù)合軌跡的算法理論, 在此基礎(chǔ)上提出一種基于圓弧過(guò)渡打印噴頭自適應(yīng)速度策略與時(shí)間模型, 對(duì)打印過(guò)程中打印頭的速度進(jìn)行前饋控制。 最后通過(guò)MATLAB 進(jìn)行模擬與仿真, 以驗(yàn)證提出的復(fù)合軌跡算法的有效性。
復(fù)合軌跡成形算法概念及其原理如下:
(1) 導(dǎo)入模型并生成切片后, 基于非均勻有理B樣條(Non-Uniform Rational B-Splines, NURBS) 曲線生成該切片的輪廓偏置曲線[14], 根據(jù)材料特性及工藝條件, 輪廓自邊界向內(nèi)偏置3-4 次, 以保持原始CAD 模型的精度, 提高模型表面質(zhì)量;
(2) 在切片的輪廓偏置曲線內(nèi), 采用Zigzag 路徑進(jìn)行填充。 Zigzag 路徑填充算法有效減少了打印過(guò)程中的空行程, 消除了打印過(guò)程中擠出頭出絲和斷絲引起的路徑缺陷;
(3) 將Zigzag 路徑拐角處的直線連接優(yōu)化為圓弧連接, 以減少打印頭在拐角處產(chǎn)生速度突變, 提高制造精度。
算法流程如圖1 所示。
圖1 復(fù)合軌跡成形算法流程Fig.1 Compound trajectory forming algorithm flow
此算法生成模型切片填充路徑的流程如圖2 所示。
圖2 復(fù)合軌跡成形算法切片路徑生成過(guò)程Fig.2 Slice path generation process of compound path forming algorithm
獲取加工零件CAD 模型的某一切片層后, 利用切片層與模型邊界上的相交點(diǎn)集建立一系列基于NURBS 的 輪 廓 曲 線[15]。 非 均 勻 有 理 B 樣 條(NURBS) 曲線是在一組(多個(gè)) 區(qū)間上分別定義的B 樣條曲線, 其特殊之處在于不需要維護(hù)分段之間的連續(xù)性,n階NURBS 擁有Cn-1級(jí)連續(xù)性, 因此它可以在較低的階次下表示出復(fù)雜的圖形[16]。 基于NURBS 的輪廓曲線的一般形式如下:
定義2, 如果基于NURBS 生成的輪廓曲線中所有控制點(diǎn)都是凸的, 則該曲線定義為I 型輪廓曲線,如圖3 (a) 所示; 否則, 該曲線為II 型輪廓曲線,如圖3 (b) 所示。
圖3 輪廓曲線示例Fig.3 Examples of profile curves: (a) type I; (b) type II
I 型輪廓曲線的偏移曲線生成。 如圖4 所示, 基于原本輪廓曲線控制點(diǎn)的中心, 生成NURBS 的偏移輪廓曲線的新控制點(diǎn), 控制點(diǎn)生成的公式如下所示:
圖4 Ⅰ型輪廓曲線的偏移曲線Fig.4 Offset curve of type I profile curve
Ⅱ型輪廓曲線的偏移曲線生成。 如圖5 所示, 將Ⅱ型輪廓曲線中的所有凹面點(diǎn)與切片輪廓中的另一個(gè)控制點(diǎn)相連, 將切片劃分成2 個(gè)或多個(gè)區(qū)域, 使得每個(gè)區(qū)域中的控制點(diǎn)都是凸面的, 圖5 中的D1、D2即為分割后的2 個(gè)子區(qū)域。 在每個(gè)區(qū)域中使用公式(3) 生成偏移曲線的新控制點(diǎn), 然后連接所有新控制點(diǎn), 即可生成Ⅱ型輪廓的偏移曲線。
圖5 Ⅱ型輪廓曲線Fig.5 Type II profile curve
Zigzag 路徑填充算法采用往復(fù)的平行直線等距分割二維水平輪廓, 這種填充方式的優(yōu)點(diǎn)是規(guī)劃過(guò)程較為簡(jiǎn)單、 易于實(shí)現(xiàn), 消除了打印過(guò)程中起弧和斷弧造成的路徑缺陷, 使打印過(guò)程更加連貫[19]。 但在路徑拐角處時(shí), 打印頭加速度變化較大, 可能出現(xiàn)過(guò)沖現(xiàn)象, 故文中將拐角處直線段變?yōu)閳A弧曲線, 以減小速度的突變, 提高制造精度。
基于圓弧過(guò)渡的Zigzag 填充路徑生成算法如下:
(1) 生成填充掃描線。 以輪廓偏置曲線為邊界,在其內(nèi)部以固定掃描線間距生成填充掃描線, 將掃描線與輪廓曲線的交點(diǎn)坐標(biāo)依次存入二維數(shù)組ScanLine[m][n]中, 其中m代表所生成的掃描線的數(shù)量,n代表每條掃描線與輪廓偏置曲線的交點(diǎn)數(shù)量。
(2) 掃描線分區(qū)。 從m=1 開始, 比較相鄰2 條掃描線的交點(diǎn)數(shù)量, 如果ScanLine[m]=ScanLine[m+1],則這2 條掃描線屬于同一區(qū)域; 如果不相等, 則ScanLine[m+1]屬于下一區(qū)域。
(3) 生成Zigzag 填充線。 分別在每個(gè)區(qū)域內(nèi)生成Zigzag 填充線, 不同區(qū)域的掃描線不連接。
(4) 將掃描線拐角處的直線連接優(yōu)化為圓弧連接, 以減小高速打印時(shí)加減速過(guò)程造成的柔性沖擊,提高零件成形質(zhì)量。
基于圓弧過(guò)渡的Zigzag 填充路徑生成流程如圖6所示。
圖6 基于圓弧過(guò)渡的Zigzag 填充路徑生成示意Fig.6 Schematic of Zigzag filling path generation based on arc transition: (a) generate fill scan path and partition;(b) generate Zigzag fill path; (c) arc optimization by connecting straight path at corners
在文中所提切片復(fù)合軌跡成形算法的基礎(chǔ)上, 為了進(jìn)一步提高制造精度及效率, 提出一種熔絲擠出成形過(guò)程中打印頭速度控制策略, 并建立相應(yīng)的時(shí)間模型以提前對(duì)整個(gè)制造過(guò)程所需的總體時(shí)間進(jìn)行預(yù)估及分析, 以達(dá)到對(duì)制造過(guò)程進(jìn)行前饋控制。
基于輪廓偏移曲線的打印示意圖如圖7 所示, 構(gòu)建基于輪廓的加工時(shí)間模型如下:
圖7 輪廓偏移曲線打印示意Fig.7 Profile offset curve print schematic
其中:si,j,l表示Ci,j中的第l段曲線;m為Ci,j中曲線段的總數(shù)。
為了優(yōu)化打印頭沿軌跡運(yùn)動(dòng)的控制, 提出一種自適應(yīng)速度策略, 對(duì)于每段si,j,l, 打印頭從零速度開始, 然后加速到額定的最大速度, 在Ci,j的末端減速到零。 最大速度、 最小速度、 加速度和減速度由所用打印機(jī)器及噴頭的規(guī)格、 性能決定[20]。 該設(shè)計(jì)旨在通過(guò)控制打印速度來(lái)減小打印頭的過(guò)沖, 以提高打印精度。
假設(shè)打印頭工作時(shí)是勻加速或勻減速運(yùn)動(dòng), 加速度為a, 在每段si,j,l中, 根據(jù)曲線長(zhǎng)度的不同, 速度的變化有3 種情況:
(1) 如圖8 (a) 所示,si,j,l的距離不足以使打印頭加速到最大速度, 為了保證打印頭在末端減速到零, 需預(yù)先計(jì)算出打印頭在si,j,l上可加速到的最大速度, 使得打印頭在當(dāng)前曲線段末端減速到零;
圖8 基于輪廓的打印頭自適應(yīng)速度變化情況Fig.8 Adaptive speed change of print head based on contour: (a) case 1; (b) case 2; (c) case 3
(2) 如圖8 (b) 所示,si,j,l的距離剛好可以使打印頭加速到額定的最大速度, 然后以大小相同、 方向相反的加速度進(jìn)行減速運(yùn)動(dòng);
(3) 如圖8 (c) 所示,si,j,l的距離可以使打印頭加速到最大速度后, 以最大速度勻速打印一段距離,再進(jìn)行對(duì)應(yīng)的勻減速運(yùn)動(dòng)。
基于以上分析及物理運(yùn)動(dòng)基本公式, 可求得打印頭從靜止以加速度a加速到最大速度vmax所需距離為
當(dāng)打印頭運(yùn)動(dòng)到曲線段末端時(shí), 速度需要減小為零, 且減速階段所需的距離和時(shí)間與加速階段相同。
當(dāng)si,j,l<2s1時(shí), 如圖8 (a) 所示, 打印頭無(wú)法加速到最大速度, 則打印頭在si,j,l上的加速距離為
減速階段所需的距離和時(shí)間與加速階段相等。
當(dāng)si,j,l>2s1時(shí), 如圖8 (c) 所示, 打印頭先加速到最大速度, 然后以最大速度勻速運(yùn)動(dòng)一段距離, 再進(jìn)行對(duì)應(yīng)的勻減速運(yùn)動(dòng), 其中, 加速和減速階段的距離為s1, 時(shí)間為t1, 則勻速運(yùn)動(dòng)的距離為
將生成的基于圓弧過(guò)渡的Zigzag 填充路徑進(jìn)一步劃分為Ⅰ型和Ⅱ型, 其中Ⅰ型曲線是Zigzag 填充路徑的主要組成部分, Ⅱ型曲線是相鄰的2 條I 型線之間的圓弧連接曲線。 圖9 所示為2 種類型曲線的示意。
圖9 Zigzag 填充路徑中2 種曲線類型Fig.9 Two types of Zigzag filling path
圖10 k=1 時(shí)打印頭自適應(yīng)速度變化情況Fig.10 Adaptive speed change of print head when k=1: (a)case 1; (b) case 2; (c) case 3
圖11 k>1 時(shí)打印頭自適應(yīng)速度變化情況Fig.11 Adaptive speed change of print head when k>1:(a) case 1; (b) case 2; (c) case 3
圖12 Ⅱ型軌跡打印頭自適應(yīng)速度變化情況Fig.12 Adaptive speed change of type II track print head:(a) case 1; (b) case 2; (c) case 3
構(gòu)建基于圓弧過(guò)渡的Zigzag 填充路徑的加工時(shí)間模型如下:
2.2.2 Ⅱ型軌跡自適應(yīng)速度策略
其中:ωmax為最大角速度;aω為額定角加速度;Rk為的旋轉(zhuǎn)半徑。
基于不同輪廓和曲線類型分別建立了相應(yīng)的時(shí)間模型, 以便提前對(duì)整個(gè)制造過(guò)程所需的總體時(shí)間進(jìn)行預(yù)估及分析, 對(duì)制造過(guò)程進(jìn)行前饋控制。
為驗(yàn)證所提出的復(fù)合路徑填充算法的可行性, 基于VC++6.0 平臺(tái), 在現(xiàn)有工藝軟件的基礎(chǔ)上進(jìn)行復(fù)合軌跡算法的開發(fā), 并增加了STL 路徑及G 代碼文件輸出功能; 使用MATLAB 編寫的程序處理STL 路徑文件, 得到路徑數(shù)據(jù)并顯示模型的三維路徑, 再根據(jù)所提出的打印頭速度策略模型進(jìn)行仿真制造, 并得到預(yù)估的制造時(shí)間; 然后基于不同算法對(duì)模型切片進(jìn)行軌跡填充, 并將G 代碼文件導(dǎo)入XXX3D 打印機(jī)進(jìn)行實(shí)際打印實(shí)驗(yàn), 以驗(yàn)證復(fù)合填充算法的優(yōu)化性。
使用一個(gè)滾針軸承座模型生成復(fù)合填充軌跡, 如圖13 所示, 該滾針軸承的尺寸為長(zhǎng)38.16 mm、 寬44.93 mm、 高14.50 mm。 分別采用光柵填充算法、Zigzag 填充算法以及文中提出的復(fù)合路徑填充算法生成該模型切片的填充路徑, 結(jié)果如14 圖所示; 輸出該模型的路徑文件并通過(guò)MATLAB 程序處理, 得到該模型的仿真結(jié)果如圖15 所示。
圖13 滾針軸承座模型Fig.13 Model of needle bearing block
圖14 不同算法生成的切片填充路徑Fig.14 Slice filling paths generated by different algorithms: (a) raster filling; (b) Zigzag filling; (c) composite algorithm filling
圖15 MATLAB 模型仿真結(jié)果Fig.15 MATLAB model simulation results: (a) raster filling; (b) Zigzag filling; (c) composite algorithm filling
根據(jù)提出的自適應(yīng)速度策略, 分別對(duì)滾針軸承座模型使用不同軌跡填充算法進(jìn)行單層切片和整體模型(文中以5 層切片為例) 的模擬制造, 其中打印頭運(yùn)動(dòng)最大速度(熔絲填充速度/熔絲輪廓速度) 為100 mm/s、 加速度為100 mm/s2、 路徑填充間距為0.2 mm, 模擬構(gòu)建的時(shí)間對(duì)比結(jié)果如圖16 所示。
圖16 不同算法模擬制造結(jié)果Fig.16 Manufacturing simulation results of different algorithms
路徑的起落點(diǎn)個(gè)數(shù)可以一定程度上反映該填充路徑的連貫程度, 起落點(diǎn)越少, 說(shuō)明制造過(guò)程中打印頭出絲、 斷絲的次數(shù)就越少, 則實(shí)際制造過(guò)程中的空行程也越少, 打印過(guò)程更加連貫, 成形效率更高。 起落點(diǎn)及空行程示意如圖17 所示。
圖17 起落點(diǎn)及空行程示意Fig.17 Schematic of landing point and empty travel
對(duì)不同算法模擬制造后生成的填充軌跡的起落點(diǎn)個(gè)數(shù)進(jìn)行統(tǒng)計(jì), 結(jié)果如圖18 所示。
圖18 不同填充路徑起落點(diǎn)個(gè)數(shù)Fig.18 Number of landing points of different filling paths
由圖16、 18 可知: 光柵填充算法的起落點(diǎn)遠(yuǎn)遠(yuǎn)多于Zigzag 填充和復(fù)合路徑填充算法, 過(guò)多的起落點(diǎn)會(huì)使得打印頭頻繁變換材料擠出的起止?fàn)顟B(tài), 進(jìn)而增加出絲與斷絲的次數(shù)和行進(jìn)過(guò)程中的空行程數(shù)量, 降低材制造效率; 而Zigzag 填充算法雖然起落點(diǎn)個(gè)數(shù)較少, 但其制造時(shí)間長(zhǎng), 且由于滾針軸承座模型輪廓為曲線, 在邊緣處軌跡搭接寬度越大, 鋸齒狀越明顯,影響零件輪廓精度和形狀; 而文中提出的復(fù)合路徑填充算法保留了Zigzag 填充算法起落點(diǎn)少的優(yōu)點(diǎn), 相比光柵填充平均減少了約85%的起落點(diǎn), 且在拐角處進(jìn)行了圓弧優(yōu)化, 有效地緩解了軌跡邊界處的階梯效應(yīng), 相較于Zigzag 填充算法減少了28.3%~36.2%的制造時(shí)間, 進(jìn)一步提高了制造效率。
實(shí)驗(yàn)所使用的打印機(jī)為Creality CR?100 3D 打印機(jī), 如圖19 所示。 其機(jī)器參數(shù)如表1 所示。
表1 打印機(jī)參數(shù)Tab.1 Printer parameters
圖19 Creality CR?100 3D 打印機(jī)實(shí)體Fig.19 3D printer entity of Creality CR?100
實(shí)驗(yàn)采用單擠出噴頭, 噴嘴直徑為0.4 mm, 設(shè)置層高為0.2 mm, 噴頭擠出溫度為200 ℃, 熱床溫度為40 ℃, 擠出材料為PLA 絲材, 絲材直徑為1.75 mm。 將光柵填充算法、 Zigzag 填充算法, 以及文中提出的復(fù)合路徑填充算法所生成的G 代碼文件輸入到3D 打印機(jī)中打印, 比較不同路徑軌跡填充算法下的打印精度。 模型單層切片打印效果如圖20 所示。
圖20 模型單層切片打印效果Fig.20 Model single-layer slice printing effect: (a)raster filling; (b) Zigzag filling; (c) composite algorithm filling
由圖20 可以得出: 使用光柵填充算法打印出的模型切片邊界處“毛刺” 較多, 材料分布不均勻,成形精度不高; 使用Zigzag 填充算法打印出的模型切片邊界處產(chǎn)生了階梯效應(yīng), 因打印頭近距離往返運(yùn)動(dòng)且材料密集導(dǎo)致產(chǎn)品的變形較大; 使用復(fù)合路徑填充算法打印出的模型切片精度較高, 邊界形狀完整, 材料分布較為均勻, 成形效果良好, 證明了所提復(fù)合路徑填充算法的可行性。
文中提出了一種復(fù)合軌跡成形算法和自適應(yīng)速度策略來(lái)優(yōu)化增材制造中打印頭軌跡的生成, 旨在提高復(fù)雜模型制造的精度和成形效率。 主要結(jié)論如下:
(1) 引入NURBS 輪廓曲線表示切片層的輪廓,以保持原始模型切片的幾何精度, 為進(jìn)一步生成內(nèi)部填充軌跡建立了良好的原始數(shù)據(jù)源。
(2) 在Zigzag 填充路徑的基礎(chǔ)上, 將拐角處的直線連接變?yōu)閳A弧優(yōu)化, 以減小高速打印時(shí)加減速過(guò)程引起的過(guò)沖現(xiàn)象, 提高零件成形質(zhì)量。
(3) 提出自適應(yīng)速度策略, 為整個(gè)打印過(guò)程建立時(shí)間分析模型, 對(duì)制造過(guò)程中打印頭的速度變化進(jìn)行前饋控制, 制造時(shí)間減少了28.3%~36.2%, 提高了制造效率。