一種提高FDM三維打印系統成型效率的算法?

巢海遠，劉 景，童 晶，張 霞

（河海大學物聯網工程學院，常州213022）

一種提高FDM三維打印系統成型效率的算法?

巢海遠，劉 景，童 晶，張 霞

（河海大學物聯網工程學院，常州213022）

目前，FDM三維打印系統制件加工時間比較長，成型效率較低。在權衡加工時間和加工精度的前提下，從切片層厚和層掃描速度控制的角度，提出了一種基于面積變化率的自適應切片層厚控制和基于周長變化率的自適應層掃描速度控制相結合的算法。首先通過二維輪廓面積變化情況計算每層切片數據的層厚值；然后用獲得的層厚數據對模型進行切片，得到每層可打印的真實輪廓；最后計算每層二維多邊形輪廓的周長，并根據周長的變化情況確定該層合適的掃描速度。實際打印測試結果表明，該算法可以有效降低打印時間，提高成型效率，驗證了算法的有效性。

三維打印；切片層厚；掃描速度；成型效率

1 引 言

FDM（Fused Deposition Modeling）三維打印技術源于快速成形技術（Rapid Prototyping），采用離散／堆積原理［1］，在計算機控制下將三維模型直接制造成三維實體。其基本原理是加熱噴頭在計算機控制下，將熱塑性聚合物材料加熱熔化，使其在熔融狀態下從噴嘴擠出，靠高溫擠出絲材的自粘結性逐層堆積成型。整個制造過程中噴頭會根據零件的界面輪廓信息，作X—Y平面運動和Z方向的運動，快速冷卻后形成界面輪廓。一層成型完成后，噴頭沿Z方向上升一個層厚高度，用同樣的方法制造下一層面并與前一個層面熔結在一起，如此循環，直至完成整個模型實體的制造。

切片模塊是軟件控制系統的核心，切片生成的文件直接影響加工零件的規模、精度和復雜程度，最終影響加工成型的效率［2－3］。對加工成型時間有重要影響的因素主要有切片層厚控制和打印掃描速度控制。目前，商業化廣泛應用的FDM成型系統大都采用定層厚切片控制，在掃描速度控制方面，有的打印系統可以在打印過程中實時的人工調節，有的則采用切片軟件設定的掃描速度。因此，加工成型比較耗時，例如十幾厘米高的模型，利用Makerbot公司的Replicator 2三維打印機，在高精度模式下，經測試需要近十小時的打印時間，這已經成為三維打印大規模應用的一個瓶頸問題。

針對以上問題，提出了一種基于面積變化率的自適應層厚控制和基于周長變化率的自適應層掃描速度控制相結合的算法，降低了加工時間，提高了成型效率。

2 基于面積變化率的自適應層厚

2.1 定層厚切片和自適應層厚切片對成型效率的影響

切片層厚的控制技術主要分為：定層厚切片和自適應層厚切片［4－6］。其中，定層厚切片存在加工時間和表面質量難以協調的問題且臺階效應明顯，如圖1（a）所示。圖1中某加工零件由A和B兩部分組成，在定層厚切片條件下加工時兩部分的層厚都相等，但因為B部分外表面曲率大，容易出現較明顯的階梯效應。為了降低階梯效應，采用較小的切片層厚，如圖1（b）所示，但這樣使得A部分的層厚也減小，顯然降低了加工效率。為提高加工效率，采用較大的切片層厚，則B部分會臺階效應明顯，如圖1（c）所示。自適應層厚切片權衡了加工時間和加工精度，如圖1（d）所示，切片時A部分的層厚選擇較大層厚，而B部分的層厚根據形狀變化來調整大小。

自適應層厚是一種變層厚切片方法，根據制件的曲面特征信息調整合適的切片層厚進行切片。但是傳統的自適應切片算法大多是基于模型表面曲率特征，計算量大，實施比較復雜，不易集成至開源的切片軟件中［7－9］。目前，主流的商業化FDM成型系統所支持的切片軟件，如Makeware、Slice3r、Cura等都沒有自適應層厚切片的功能，而是在可打印范圍內讓用戶輸入指定的層厚值，即按定層厚切片。所以用戶在成型精度和成型效率之間很難權衡，為了保證成型精度，往往選用在允許范圍內的最小層厚進行打印，打印過程十分耗時。

圖1 定層厚與自適應層厚切片效果示意圖

2.2 基于面積變化率的層厚計算

為減少因形狀變化而造成的誤差，減小臺階效應，同時提高加工過程中的效率，提出基于截面面積變化率的自適應層厚控制。自適應層厚控制的原則是根據制件表面形狀的變化確定切片層厚，變化較緩的部分其層片厚度大，變化較快的部分層片厚度小。根據以上原則，基于面積變化率的自適應層厚是通過計算相鄰兩層面面積的變化率實現的，即計算出當前層面的面積和上一層面的面積，得到面積變化率，如果面積變化率不大于給定閾值，則滿足條件，得到當前層合適的切片層厚，否則為不滿足條件，繼續計算得到合適的層厚。其流程如圖2所示。

算法步驟為：

Step1：讀入STL文件，i＝0，給定初始切片層厚t0，切片一層計算該層的二維截面面積S0；

Step2：遞增一個初始層厚ti＝t0，同時計數器i＝i＋1，計算截面面積Si；

Step4：判斷是否ti≥tmin，若成立，則記錄該層切片層厚數據ti，轉Step5；否則ti＝tmin并記錄該層切片ti層厚數據，轉Step5；

Step5：判斷是否∑ti≤H，其中H為STL模型的高度，若成立，則轉Step2；否則，算法結束。

其中二維截面面積的計算方法為三角形分割法，是將二維多邊形切片輪廓分割成三角形，這些三角形的面積總和即是層面截面面積，三角形分割法得到的截面面積是比較精確的。

3 基于周長變化率的層掃描速度

3.1 掃描速度對成型效率影響

掃描速度指的是噴頭隨噴嘴運動的速度，每層的成型時間與該層的掃描速度及該層輪廓形狀的復雜度有關。若層面輪廓形狀簡單，填充速度快，則該層成型的時間就短；反之，時間就長［10］。在加工一些截面很小的層面時，若一層的成形時間太短，會導致前一層還來不及固化成形，下一層就接著堆積，將引起坍塌和拉絲現象。因此可以采用較慢的掃描速度，以增加層面成型時間。在加工一些截面很大的層面時，若一層的成形時間太長，會導致前一層截面已完全冷卻凝固后才開始堆積后一層，從而導致層間粘接不牢固現象，出現成型件的開裂傾向。因此可以采用較快的掃描速度，從而減少層面成型時間。在實際打印時，控制好每層的成型時間，才能獲得較高精度的制件。

圖2 基于面積變化率的層厚計算

在提高加工效率的前提下，為了得到較高精度的成型件，掃描速度還需要與擠出速度（絲材擠出噴嘴的速度）相協調。一般來說掃描速度增大，擠出速度也相應增大；掃描速度減小，擠出速度也應相應減小。掃描速度比擠出速度過快，則材料填充不足，會有“拉絲”、“斷絲”情況產生，甚至有時會出現表面空缺現象，難以成型。相反，掃描速度比擠出速度過慢，熔絲堆積在噴頭上，使成型面材料分布不均勻，表面會有疙瘩，同時引起噴嘴“碳化”，“粘附”現象，阻礙下一步的加工并影響最終成型件質量。因此，在實際打印過程中，掃描速度和擠出速度相互影響，且存在一個合理的匹配范圍，二者比值必須在此范圍內才會得到較高的制件精度。掃描速度與擠出速度匹配的合理范圍滿足［11］：

式中vj為擠出速度，vi為掃描速度，α為成型時出現斷絲現象的臨界值，β為出現粘附現象的臨界值。

3.2 基于周長變化率的自適應層掃描速度計算

為減少每層輪廓不必要的加工時間，提高加工過程中的效率，提出根據周長變化控制的自適應層掃描速度控制。自適應層掃描速度控制的原則是根據切片截面輪廓周長的大小確定層掃描速度，輪廓周長大的層掃描速度大，輪廓周長小的層掃描速度小。基于這一原則，周長變化率的自適應層掃描速度是通過計算一層輪廓周長與最小輪廓周長的比率來確定一層的掃描速度。首先根據基于截面面積變化率的自適應層厚算法計算得到所有的層厚值，以此層厚數據對輸入的STL文件進行切片；然后計算所有切片截面的二維輪廓周長，找出最小輪廓周長并確定與最小輪廓周長合適的掃描速度和擠出速度；最后每一層的層掃描速度根據該層輪廓周長與最小輪廓周長的比率來確定。其流程如圖3所示。

算法步驟為：

Step1：根據基于截面面積變化率的自適應層厚計算得到所有的層厚值ti（i＝0，1，...，n），以此層厚數據對輸入的STL文件進行切片；

Step2：計算所有切片截面的二維輪廓周長Pi（i＝0，1，...，n），找出最小輪廓周長Pmin；

Step3：賦值i＝0，設置掃描速度vmin＿i和擠出速度vmin＿j，并判斷是否在匹配的合理范圍（α，β）內，若不在，則重新設置vmin＿i和 vmin＿j；若在，則轉Step4；

Step4：計算第i層的掃描速度vi和擠出速度vj。

其中，掃描速度vi與擠出速度vj之間的比值始終滿足：，即始終在合理的匹配范圍之內，打印過程中可以保證較高的制件精度。

Step5：將掃描速度vi和擠出速度vj線性歸一化到打印系統能承受的打印速度范圍內：

其中，Vi＿min和Vi＿max分別為打印系統掃描速度能承受的最小值和最大值，Vj＿min和Vj＿max分別為打印系統擠出速度能承受的最小值和最大值。

Step6：判斷是否i≤n，若成立，則轉Step4，若不成立，則算法結束。

圖3 基于周長變化率的層掃描速度計算

4 實驗結果

CuraEngine是由Ultimaker團隊開發的以C＋＋語言寫成的跨平臺切片引擎，這是目前為止速度最快的切片引擎。CuraEngine完全開放源代碼，主要特點是切片速度快，可擴展性強［12］。在開源切片引擎CuraEngine的基礎上，以Visual Studio 2010為開發平臺，在Windows環境下采用C＋＋編程實現了所述算法。電腦硬件配置為Intel（R）Core（TM）2 Duo CPU T6670＠2.20GHz，2GM RAM。FDM成型設備為MakerBot Replicator 2，該設備支持物體最大打印尺寸285×153×155mm，層分辨率在100微米至400微米間，定位精度X／Y軸11微米，Z軸2.5微米，打印材料為PLA塑料。

測試模型如圖4（a）所示。模型的包圍盒大小為：長度75mm，寬度35mm，高100mm。實際打印測試結果如圖4（b）、（c）、（d）所示，拍攝效果對比不是很明顯，將切片仿真結果的局部細節在同一尺度下放大對比，可以看出算法有效降低了臺階效應。成型時間對比如表1所示，通過對比可以看出，運用基于面積變化率的自適應層厚算法切片層數明顯下降，同時運用基于周長變化率的自適應掃描速度算法在一定程度上降低了打印時間。相對于最小定層厚，自適應層厚打印成型時間節省了45.71%，采用自適應層厚和自適應掃描速度相結合的算法打印成型時間節省了49.52%。其中測試的模型形狀較簡單，對于工業設計中規則的CAD模型，算法是有一定優勢的，但是對于個性化制造領域及復雜的實際模型，算法有待進一步研究測試。

表1 成型時間對比

5 結束語

成型速度快和原型制作質量好是當今三維打印領域中兩個重要的發展方向［13］。提出了一種基于面積變化率的自適應切片層厚控制和基于周長變化率的自適應層掃描速度控制相結合的算法。為獲得切片層厚數據，通過計算輪廓面積變化率得出每層合適的層厚值；接著以該層厚數據文件對模型進行切片，得到每層的二維多邊形輪廓；最后利用周長變化率確定每層輪廓合適的層掃描速度。實際打印測試結果表明，該算法在一定程度上可以降低打印時間，提高加工成型效率。測試的實例形狀簡單且模型較小，所以自適應層掃描速度的作用體現的并不明顯。但是對于應用廣泛的工業設計、生物醫療以及機械模具加工等行業，實際打印的模型尺寸都非常大，加工成型往往需要十幾個小時，此時一定比例的時間節省會帶來很大的經濟效益，這對于推動三維打印的普及應用有重要意義。

圖4 測試結果

［1］朱林泉，白培康，朱江淼.快速成型與快速制造技術［M］.北京：國防工業出版社，2003.ZHU Lin－quan，BAI Pei－kang，ZHU Jiang－miao.Rapid prototyping and rapid manufacturing technology［M］.Beijing：National Defense Industry Press，2003.

［2］劉斌，黃樹槐.快速原型制造技術中實時切片算法的研究與實現［J］.計算機輔助設計與圖形學學報，1997，9（6）：488－493.LIU Bin，HUANGShu－huai.The dynamic controlling of spreading stream line in pump CAD system［J］.Journal of Computer－Aided Design＆Computer Graphics，1997，9（6）：488－493.

［3］馬巧梅，朱林泉.快速成型系統中自適應的直接切片算法的研究［J］.計算機工程與設計，2008，29（8）：2075－2077.MA Qiao－mei，ZHU Lin－quan.Research on adaptive direct slicing algorithm in rapid prototyping system［J］.Journal of Computer Engineering and Design，2008，29（8）：2075－2077.

［4］畢曉亮，朱昌明，侯麗雅，等.快速成型中的自適應切片方法研究［J］.計算機應用研究，2002，19（4）：21－22.BIXiao－liang，ZHU Chang－ming，HOU Li－ya，et al.Application of an adaptive slicing procedure in RPM software［J］.Journal of Application Research of Computers，2002，19（4）：21－22.

［5］趙吉賓，劉偉軍.快速成型技術中分層算法的研究與進展［J］.計算機集成制造系統，2009，15（2）：209－221.ZHAO Ji－bin，LIU Wei－jun.Recent progress in slicing algorithm of rapid prototyping technology［J］.Journal of Computer Integrated Manufacturing Systems，2009，15（2）：209－221.

［6］王天明，金燁，習俊通，等.FDM工藝過程中絲材的粘結機理與熱學分析［J］.上海交通大學學報，2006，40（7）：1230－1233，1238. WANG Tian－ming，JIN Ye，XI Jun－tong，et al.The adhesivemechanism and thermal analysis of fibers in the FDM process［J］.Journal of Shanghai Jiao Tong University，2006，40（7）：1230－1233，1238.

［7］彭安華，張劍峰，張江林，等.FDM工藝參數對制件精度影響的實驗研究［J］.淮海工學院學報（自然科學版），2008，17（1）：21－24.PENG An－hua，ZHANG Jian－feng，ZHANG Jianglin，et al.Experimental research into the influence of technical parameters on pats precision［J］.Journal of Huaihai Institute of Technology（Natural Science Edition），2008，17（1）：21－24.

［8］何新英，陶明元，葉春生，等.FDM工藝成形過程中影響成形件精度的因素分析［J］.機械與電子，2004（9）：77－78.HE Xin－ying，TAO Ming－yuan，YE Chun－sheng，et al.Analysis for factors that affect the precision of the shaped pieces on fused deposition modeling process［J］.Journal of Machinery＆Electronics，2004，（9）：77－78.

［9］胡汝霞.基于RE／RP直接集成的模型切片處理方法研究［D］.江蘇：江蘇大學機械制造及其自動化專業，2005.HU Ru－xia.Research of slicing method directly integrated in RE／RP［D］.Jiangsu：Automation and machinery manufacturing professional of Jiangsu university，2005.

［10］A Dolenc，I M?kel?.Slicing procedures for layered manufacturing techniques［J］.Computer－Aided Design，1994，26（2）：119－126.

［11］Prashant Kulkarni，Debasish Dutta.An accurate slicing procedure for layered manufacturing［J］.Computer－Aided Design，1996，28（9）：683－697.

［12］K Mani，PKulkarni，D Dutta.Region－based adaptive slicing［J］.Computer－Aided Design，1999，31（5）：317－333.

［13］CURA.Website of the source code of cura［EB／OL］.https：／／github.com／Ultimaker／CuraEngine，2013.

An Algorithm for Im proving Efficiency of Fused Deposition Molding 3D Printing System

CHAO Hai－yuan，LIU Jing，TONG Jing，ZHANG Xia
（College of Internet of Things Engineering，HoHai University，Changzhou 213022，China）

Currently，the FDM（Fused Deposition Modeling）3D printing system is conducted with long time and low efficiency.Considering process time and process accuracy，in control perspective of slice thickness and scan speed，an algorithm is proposed to reduce printing time and improve molding efficiency，which combining adaptive slice thickness based on the rate of change of the area and adaptive scan speed based on the rate of change of the perimeter.Firstly，the rate of change of the area is used to calculate the value of slicing thickness of each layer.According to the slicing data，then the mesh is sliced to obtain actual printing contours of each layer.Finally，the perimeter of two－dimensional polygon of each layer is computed and the rate of change of the perimeter is used to determine the appropriate scanning speed of each layer.The actual printing test results validate the feasibility of the proposed algorithm and show that it can effectively reduce the printing time and improve the efficiency of forming.

3D printing；Slice thickness；Scan speed；Molding efficiency

10.3969／j.issn.1002－2279.2014.06.012

TP391.7；TH166

：A

：1002－2279（2014）06－0038－05

國家自然科學基金資助項目（61202284）；常州市應用基礎研究項目（CJ20130045）

巢海遠（1989－），男，江蘇常州人，碩士研究生，主研方向：計算機圖形學、三維打印。

2014－09－10