基于臨界傾角約束的樹形支撐結(jié)構(gòu)生成算法

沈振宏　戴　寧　李大偉　吳長友

南京航空航天大學(xué)，南京，210016

?

基于臨界傾角約束的樹形支撐結(jié)構(gòu)生成算法

沈振宏戴寧李大偉吳長友

南京航空航天大學(xué)，南京，210016

為節(jié)省三維打印的材料和時(shí)間，針對支撐結(jié)構(gòu)體積優(yōu)化的問題，提出了一種基于臨界傾角約束的樹形支撐結(jié)構(gòu)生成算法。首先識別待支撐區(qū)域，并基于面積自適應(yīng)采樣獲得待支撐點(diǎn)；再利用臨界傾角約束進(jìn)行樹形結(jié)構(gòu)節(jié)點(diǎn)的計(jì)算，并存儲同一棵樹上各節(jié)點(diǎn)的連接關(guān)系；采用掃掠和泊松曲面重建方法，生成樹形支撐實(shí)體，最后通過實(shí)例進(jìn)行了驗(yàn)證。實(shí)驗(yàn)結(jié)果表明，所提出的方法在保證打印穩(wěn)定性的前提下，有效減少了支撐材料，縮短了打印時(shí)間。

三維打??；臨界傾角約束；樹形支撐；泊松曲面重建

0　引言

三維打印工藝過程中普遍存在著零件變形、翹曲、坍塌等問題。例如在SLM(selective laser melting)工藝過程中，熱應(yīng)力超過材料強(qiáng)度時(shí)發(fā)生的塑性變形以及逐層鋪粉時(shí)的橫向擾動等因素都會導(dǎo)致成形質(zhì)量下降，因此需要設(shè)計(jì)支撐結(jié)構(gòu)進(jìn)行部分熱量傳導(dǎo)和零件固定。FDM(fused deposition modeling)工藝中，每個(gè)層面都是在上一層上堆積而成，上一層對當(dāng)前層起到定位和支撐的作用。當(dāng)前層的形狀發(fā)生較大變化，如出現(xiàn)懸臂結(jié)構(gòu)或較大傾斜面時(shí)，上層輪廓就不能給當(dāng)前層提供充分的定位和支撐，這時(shí)就需要設(shè)計(jì)一些輔助結(jié)構(gòu)用于支撐。因此，支撐結(jié)構(gòu)的生成，是三維打印成形工藝中的一項(xiàng)關(guān)鍵技術(shù)。

支撐生成算法首先需要尋找模型中的待支撐部位。一類方法是先獲得模型的分層輪廓數(shù)據(jù)，再利用分層數(shù)據(jù)進(jìn)行識別。Chalasani等[1]對各層多邊形進(jìn)行合并，得到模型截面的最小包絡(luò)矩形，進(jìn)而得到模型的包圍盒，該方法以SLC文件為輸入，可以很快得到截面包絡(luò)矩形的信息，但采用包圍盒的方法，需要人工刪除一些不需要添加支撐的區(qū)域。Huang等[2]自頂而下地對相鄰兩層實(shí)體切片進(jìn)行求差、并的多邊形布爾運(yùn)算，算法有效利用了材料的自支撐能力，但是由于涉及到任意多邊形的布爾運(yùn)算和任意多邊形的雙向偏置，具有一定的復(fù)雜度，隨著多邊形和節(jié)點(diǎn)的增多，算法的耗時(shí)會顯著增加。另一類方法是直接在三維模型上進(jìn)行識別。Lan等[3]先對模型中待支撐的邊進(jìn)行處理，然后查找待支撐面片，該方法針對的模型相對簡單，其適用范圍具有一定局限性。Alexander等[4]將模型表面質(zhì)量、加工時(shí)間、支撐體數(shù)量作為目標(biāo)函數(shù)，通過求多目標(biāo)問題的最優(yōu)解集確定模型的分層方向。該算法采用參數(shù)肩高進(jìn)行表面精度計(jì)算，并對肩高求值進(jìn)行了優(yōu)化，但是算法的時(shí)間復(fù)雜度較高。

對于支撐體的結(jié)構(gòu)設(shè)計(jì)，最簡單的方法是包圍盒法，即求出模型的最小包絡(luò)長方體。為了方便加工后支撐體的去除，Strano等[5]采用較稀疏的內(nèi)部結(jié)構(gòu)進(jìn)行填充，Kritchman等[6]使用了可溶材料打印支撐體。但包圍盒法沒有考慮零件原型的形態(tài)特征，因此它也是一種最耗費(fèi)支撐材料和成形時(shí)間的方法。為了減小包圍盒法中支撐體的體積，陳巖等[7]采用直臂支撐結(jié)構(gòu)，該算法能自動添加支撐點(diǎn)和支撐桿，但算法中的支撐傳遞參數(shù)無法自動選取，僅能通過實(shí)驗(yàn)估算。Huang等[8]提出了一種斜壁結(jié)構(gòu)，該支撐結(jié)構(gòu)分為三種類型，通過計(jì)算待支撐區(qū)域多邊形的穩(wěn)定系數(shù)，進(jìn)而確定添加的支撐類型，但該算法減小的支撐體積有限，支撐與模型的接觸面積較大。Dumas等[9]采用了一種橋梁支撐結(jié)構(gòu)，該結(jié)構(gòu)由一系列水平橋梁和豎直支柱構(gòu)成，材料的耗費(fèi)較少，但是該方法總體路線較為復(fù)雜，水平橋梁中錨線段的選取不明確。

本文提出一種基于臨界傾角約束的樹形支撐結(jié)構(gòu)生成算法，該算法對于待支撐的點(diǎn)集，利用下層支撐點(diǎn)的臨界傾角約束獲得下層支撐點(diǎn)位置，構(gòu)建節(jié)點(diǎn)結(jié)構(gòu)體，生成樹形結(jié)構(gòu)，最后通過掃掠和泊松曲面重建生成樹實(shí)體。

1　技術(shù)流程

基于臨界傾角約束的樹形支撐結(jié)構(gòu)生成算法的技術(shù)流程如圖1所示。首先識別模型中的待支撐區(qū)域，根據(jù)面片方向與加工方向(z軸正方向)的夾角保存模型中的待支撐面片，對待支撐面片集合進(jìn)行相鄰面片的合并，獲得若干獨(dú)立的待支撐區(qū)域。然后基于面積自適應(yīng)采樣獲得待支撐點(diǎn)。之后，進(jìn)入算法的核心部分，即利用臨界傾角約束條件自上而下逐層生成樹形支撐結(jié)構(gòu)。以待支撐點(diǎn)為輸入，利用下層支撐點(diǎn)的位置約束條件進(jìn)行樹形結(jié)構(gòu)節(jié)點(diǎn)的計(jì)算，構(gòu)建節(jié)點(diǎn)結(jié)構(gòu)體，生成若干獨(dú)立的樹結(jié)構(gòu)。最后，通過截面單元掃掠和泊松曲面重建生成支撐實(shí)體。

圖1　基于臨界傾角約束的樹形支撐結(jié)構(gòu)算法流程圖

2　基于臨界傾角約束的樹形支撐結(jié)構(gòu)生成

2.1待支撐區(qū)域識別

零件待支撐區(qū)域的自動識別，首先要讀入零件的模型數(shù)據(jù)。三維打印系統(tǒng)多用STL格式文件表示零件的數(shù)據(jù)，即用大量的三角面片近似整個(gè)零件的表面。根據(jù)零件模型中面片方向與加工方向的夾角α與臨界傾角β的大小關(guān)系，識別需要添加支撐的面片。α≤β時(shí)，由于材料自身的粘力作用，無需添加輔助支撐結(jié)構(gòu)；α>β時(shí)，需要對面片添加支撐(圖2粗實(shí)線區(qū)域)。同時(shí)對向下傾斜且不需要支撐區(qū)域的局部最低的點(diǎn)和邊進(jìn)行識別。

圖2　待支撐區(qū)域識別

對于識別的待支撐面片集合Tsup，需要將集合中共邊的面片進(jìn)行合并，得到若干獨(dú)立的待支撐區(qū)域A0，A1，…，An，面片合并結(jié)果為

(1)

Ai∩Aj=?i，j=1，2，…，n且i≠j

算法采用區(qū)域增長法進(jìn)行相鄰面片的合并。首先任取集合中的一個(gè)面片作為種子，判斷該種子的共邊三角面片是否在集合Tsup中。若共邊面片在集合中且未存為種子面片，則將其存到種子集中。繼續(xù)尋找與其他種子共邊的三角面片，直到種子集中的所有種子都查找過為止。重復(fù)上述步驟，完成各個(gè)區(qū)域的面片合并，流程如圖3所示。

圖3　相鄰面片合并算法流程圖

2.2基于面積自適應(yīng)采樣的待支撐點(diǎn)選擇

獲得待支撐區(qū)域后，需要在區(qū)域內(nèi)采樣一系列的待支撐點(diǎn)。算法首先獲得區(qū)域投影的凸包，再對凸包進(jìn)行柵格劃分，根據(jù)柵格中面片的大小選取待支撐的點(diǎn)。

獲得模型投影的凸包有兩種方法。一種方法是先提取區(qū)域邊界，再求邊界點(diǎn)的凸包。邊界提取，即先保存面域中的非公用邊，得到組成邊界的各條線段，再利用種子擴(kuò)散法將相鄰線段連接起來。另一種方法是直接將面域中的各頂點(diǎn)往某一平面投影，得到這些頂點(diǎn)投影的凸包。

得到投影的凸包后，采用一定的柵格間距d對凸包進(jìn)行網(wǎng)格劃分。柵格間距d的選取與材料的臨界打印長度(材料在沒有支撐的情況下保持不變形的極限打印長度)L有關(guān)。為了保證支撐的穩(wěn)定性，同時(shí)避免選擇的支撐點(diǎn)過于密集，選取經(jīng)驗(yàn)值d=L/2。然后從柵格的左下頂點(diǎn)開始，逐行往上判斷每個(gè)柵格中的頂點(diǎn)個(gè)數(shù)。當(dāng)柵格中的點(diǎn)個(gè)數(shù)為0時(shí)，插入相應(yīng)的點(diǎn)并保存；當(dāng)柵格中的點(diǎn)個(gè)數(shù)大于1時(shí)，計(jì)算這些點(diǎn)的中心點(diǎn)并保存。最后通過射線求交，求得這些點(diǎn)在模型上的交點(diǎn)，得到一系列的待支撐點(diǎn)。

2.3基于臨界傾角約束的支撐結(jié)構(gòu)生成

獲得待支撐點(diǎn)后，算法基于自上而下的生長策略，利用下層支撐點(diǎn)的臨界傾角約束條件，生成樹形支撐結(jié)構(gòu)。對于待支撐點(diǎn)集P中的點(diǎn)p，假設(shè)通過下方的點(diǎn)s進(jìn)行支撐，則向量sp與加工方向的夾角不大于臨界傾角β。因此s可落于以點(diǎn)p為頂點(diǎn)、錐角為2β的圓錐體上或體內(nèi)的任意位置。圖4中，點(diǎn)p1、p2為零件原型上待支撐的點(diǎn)，c1為點(diǎn)p1對應(yīng)的圓錐體，c2為點(diǎn)p2對應(yīng)的圓錐體，則同時(shí)支撐p1和p2的下層節(jié)點(diǎn)可取圓錐體c1和c2相交區(qū)域中的任意點(diǎn)。本文選取該區(qū)域的最高點(diǎn)s，因?yàn)榇藭r(shí)線段p1s、p2s的長度之和是最小的，即打印的支撐材料最少。

圖4　錐約束求下層支撐點(diǎn)示意圖

通過待支撐點(diǎn)的配對求交，并將交點(diǎn)插入到點(diǎn)集P中，直到點(diǎn)集P中所有點(diǎn)都計(jì)算過為止。該算法的輸入為待支撐的點(diǎn)集P(按z從大到小存儲)和以點(diǎn)集P中各點(diǎn)為頂點(diǎn)生成的對應(yīng)的圓錐集C。整體支撐結(jié)構(gòu)的生成通過遍歷點(diǎn)p∈P和其對應(yīng)的圓錐c∈C執(zhí)行以下操作：

(1)查找圓錐c和模型m的最近交點(diǎn)sm，圓錐c與其他圓錐{ci|ci∈C且ci≠c}的最近交點(diǎn)；

(2)在圓錐c與其他圓錐{ci|ci∈C且ci≠c}的交點(diǎn)集合中，選擇交點(diǎn)到相應(yīng)兩圓錐頂點(diǎn)p、pi距離之和L1最小的交點(diǎn)smin。如果線段psmin的長度大于線段psm的長度，刪除點(diǎn)p和相應(yīng)的圓錐c，保存線段psm，對下個(gè)點(diǎn)重復(fù)執(zhí)行步驟1；

(3)將交點(diǎn)smin插入到點(diǎn)集P中，保存線段psmin、pismin；

(4)生成點(diǎn)smin對應(yīng)的圓錐，刪除點(diǎn)p、圓錐c和ci；

(5)當(dāng)P非空時(shí)，對下一個(gè)點(diǎn)重復(fù)上述步驟。

通過配對求交獲得樹形結(jié)構(gòu)中上下節(jié)點(diǎn)的連接關(guān)系，建立node結(jié)構(gòu)體來存儲樹形結(jié)構(gòu)中的各個(gè)節(jié)點(diǎn)，其定義如下：

Struct node

{intelement；

intlevel；

boolbvisit；

node*parent， lchild， rchild；

}

其中，成員變量element、level、bvisit分別為節(jié)點(diǎn)的索引值、節(jié)點(diǎn)在樹結(jié)構(gòu)中的層數(shù)、節(jié)點(diǎn)的訪問標(biāo)志，同時(shí)，每個(gè)節(jié)點(diǎn)還包含指向父節(jié)點(diǎn)和左右子節(jié)點(diǎn)的指針。通過node結(jié)構(gòu)體將同一棵樹上的節(jié)點(diǎn)關(guān)聯(lián)起來，同時(shí)按照樹結(jié)構(gòu)層數(shù)從小到大的順序逐層地把節(jié)點(diǎn)存在列表中。

2.4基于泊松曲面重建的支撐實(shí)體生成

2.4.1掃掠獲得初始曲面頂點(diǎn)

得到樹形支撐結(jié)構(gòu)后，需要把它轉(zhuǎn)化為實(shí)體，然后將帶有支撐結(jié)構(gòu)的原型零件轉(zhuǎn)化為STL文件，經(jīng)后續(xù)分層處理得到零件和支撐的截面輪廓，再對其進(jìn)行路徑規(guī)劃，從而得到零件及支撐體的加工路徑[10-11]。

算法通過分層保存的節(jié)點(diǎn)信息，使用水平截面單元逐層掃掠各層各節(jié)點(diǎn)到其子節(jié)點(diǎn)的路徑，獲得初始曲面頂點(diǎn)。然后對這些頂點(diǎn)進(jìn)行編號和存儲，生成網(wǎng)格曲面。此時(shí)生成的曲面在相交的區(qū)域沒有進(jìn)行相應(yīng)的處理(圖5三角形區(qū)域)，如果直接打印的話，將產(chǎn)生振動，影響打印的穩(wěn)定性。

圖5　掃掠示意圖

2.4.2泊松曲面重建

由于掃掠生成的曲面是相互獨(dú)立的，并且存在相交區(qū)域，因此需要對掃掠曲面進(jìn)行曲面重建，使每棵樹成為一個(gè)獨(dú)立的整體。本文采用的泊松曲面重建方法[12]首先計(jì)算三維指標(biāo)函數(shù)χ(點(diǎn)在模型內(nèi)部時(shí)，χ=1；點(diǎn)在模型外部時(shí)，χ=0)，然后采用Marching Cubes算法提取一個(gè)適當(dāng)?shù)牡戎得鎭慝@得重建曲面。

重建算法首先計(jì)算模型的指標(biāo)函數(shù)。指標(biāo)函數(shù)的梯度是一個(gè)在表面附近非零，在其他任意位置都為零的向量場。因此，帶法矢的表面點(diǎn)樣本可以看作是指標(biāo)函數(shù)梯度的樣本。計(jì)算指標(biāo)函數(shù)的問題因此可簡化為梯度算子的反算，即找到標(biāo)量函數(shù)χ，使其梯度最佳逼近由采樣點(diǎn)定義的向量場V，即求minχ‖χ-V‖。如果采用散度算子，該問題就可以轉(zhuǎn)化為一個(gè)標(biāo)準(zhǔn)的泊松問題，即計(jì)算拉普拉斯算子等于向量場V的散度的標(biāo)量函數(shù)χ:

Δχ=·χ=·V

(2)

泊松曲面重建的重要步驟如下(重建結(jié)果如圖6所示)：

(1)根據(jù)采樣點(diǎn)的位置信息定義一個(gè)八叉樹?，然后給八叉樹的每一個(gè)節(jié)點(diǎn)o∈?，定義一個(gè)節(jié)點(diǎn)函數(shù)F0，之后定義函數(shù)空間f?,F≡Span{F0}，其結(jié)構(gòu)類似于傳統(tǒng)的小波表示；

(2)選擇基函數(shù)F；

(3)定義向量場V；

(4)求解minχ‖χ-V‖的最小值；

(5)抽取等值面。

圖6　泊松曲面重建示意圖

另外，在進(jìn)行泊松曲面重建之前，算法對樹形支撐結(jié)構(gòu)進(jìn)行了局部優(yōu)化：①在支撐結(jié)構(gòu)與模型的接觸區(qū)域，采用小圓接觸。小圓的半徑r=0.3 mm，圓平面法矢方向與待支撐點(diǎn)的法向量同向。采用小圓接觸可以減小支撐體與模型表面的接觸面積，方便后處理中支撐結(jié)構(gòu)的剝離。②在支撐結(jié)構(gòu)與底板的接觸處采用較大的底面圓，底面圓半徑為4.0 mm，從而更好地黏住底板，增加支撐結(jié)構(gòu)穩(wěn)定性。

3　實(shí)驗(yàn)結(jié)果與分析

3.1實(shí)驗(yàn)條件

提出的算法在MicrosoftVisualStudio2008軟件開發(fā)平臺上實(shí)現(xiàn)，并通過OpenGL進(jìn)行模型顯示。建模后的模型在美國MakerBot公司的MakerBotReplicator2打印機(jī)上完成打印，打印材料為ABS塑料絲，打印速度為40mm/s，移動速度為60mm/s，層高為0.2mm。

3.2臨界傾角β標(biāo)定實(shí)驗(yàn)

由于臨界傾角β并不是一個(gè)定值，不同的打印機(jī)和打印材料對應(yīng)的臨界傾斜角可能不同。為此設(shè)計(jì)了不同傾角的樹枝并進(jìn)行實(shí)驗(yàn)，如圖7所示。得到結(jié)果如表1所示，表1中的編號與圖7的編號相對應(yīng)。

圖7　臨界傾角標(biāo)定試驗(yàn)

編號(1,1)(1,2)1,3)(1,4)(1,5)傾角(°)010203040變形情況無無無無無編號(2,1)(2,2)(2,3)(2,4)(2,5)傾角(°)5060708090變形情況輕微輕微稍重稍重嚴(yán)重編號(3,1)(3,2)(3,3)(3,4)(3,5)傾角(°)3031323334變形情況無無無無無編號(4,1)(4,2)(4,3)(4,4)(4,5)傾角(°)3536373839變形情況無無無無無編號(5,1)(5,2)(5,3)(5,4)(5,5)傾角(°)4041424344變形情況無無無無無編號(6,1)(6,2)(6,3)(6,4)(6,5)傾角(°)4546474849變形情況輕微輕微輕微輕微輕微

由表1可以看出，當(dāng)傾角大于45°的時(shí)候，樹枝開始發(fā)生變形，為了保證打印過程中的穩(wěn)定性，保留一定的余量，取臨界傾角β為40°。

3.3支撐建模及實(shí)驗(yàn)

獲得了打印機(jī)對應(yīng)的臨界傾斜角β后，基于臨界傾角約束逐層生成樹形支撐結(jié)構(gòu)的過程如圖8所示，首先識別待支撐區(qū)域，然后基于面積自適應(yīng)采樣獲得待支撐點(diǎn)，再根據(jù)下層支撐點(diǎn)的臨界傾角約束計(jì)算樹形結(jié)構(gòu)中節(jié)點(diǎn)的位置，基于自上而下的生長策略生成樹形支撐結(jié)構(gòu)，最后通過泊松重建獲得樹形結(jié)構(gòu)實(shí)體。

(a)識別待支撐區(qū)域　　　(b)計(jì)算待支撐點(diǎn)

(c)生成支撐結(jié)構(gòu)(狀態(tài)1)　(d)生成支撐結(jié)構(gòu)(狀態(tài)2)

(e)生成支撐結(jié)構(gòu)(狀態(tài)3)　　(f)泊松重建生成實(shí)體

由于打印過程中存在擾動，影響樹形結(jié)構(gòu)整體的打印穩(wěn)定性，因此算法對于圖8e中的支撐結(jié)構(gòu)底層進(jìn)行了輔助支撐，再進(jìn)行泊松曲面重建，得到圖8f所示的結(jié)果。

為了客觀評價(jià)提出的算法，將算法生成的支撐模型與打印機(jī)內(nèi)置軟件MakerWare和Magics RP v16.03(比利時(shí)Materialise公司)生成的支撐模型作對比，圖9所示為懸吊球模型在不同方法支撐下的打印結(jié)果。圖10所示為添加了樹形支撐結(jié)構(gòu)的牙頜托盤模型和人體骨頭模型。3個(gè)模型的打印結(jié)果對比如表2所示。

(a)MakerWare　　　　　　　(b)Magics RP

(c)本文算法　　　　　(d)去除支撐后的模型圖9　模型打印結(jié)果對比

圖10　添加支撐后的牙頜托盤和骨頭模型

對比參量懸吊球牙頜托盤人體骨頭打印時(shí)間(min)MakerWare9493102MagicsRP788085本文算法747468材料消耗(g)MakerWare8.016.714.6MagicsRP2.914.112.4本文算法2.712.310.0

由表2可以得出，相比MakerWare，本文算法降低了20%左右的時(shí)間消耗和近40%的材料消耗；相比Magics RP，本文算法在時(shí)間和材料耗費(fèi)上也均占有優(yōu)勢。在模型表面成形質(zhì)量方面，由于算法生成的支撐與物體模型的接觸面積很小，所以生成的支撐對模型表面損害較小，同時(shí)降低了剝除支撐體的難度。

4　結(jié)論

(1)針對包圍盒法和直臂支撐法生成的支撐結(jié)構(gòu)體積大、打印時(shí)間較長的問題，提出一種基于臨界傾角約束的樹形支撐結(jié)構(gòu)自動生成算法。算法利用下層支撐點(diǎn)的臨界傾角約束條件，計(jì)算樹形支撐結(jié)構(gòu)節(jié)點(diǎn)的位置。

(2)針對支撐體生成、去除和支撐體穩(wěn)定性的問題，本文通過掃掠獲得初始支撐結(jié)構(gòu)曲面，采用泊松曲面重建生成樹形支撐實(shí)體。在支撐結(jié)構(gòu)與模型的接觸位置采用小圓接觸的方式進(jìn)行優(yōu)化，同時(shí)對結(jié)構(gòu)中承載較重的底部樹干進(jìn)行了輔助支撐。實(shí)驗(yàn)結(jié)果表明，該方法在保證模型原型正常打印的前提下，有效地減小了支撐體的體積和數(shù)量，同時(shí)縮減了打印時(shí)間。

(3)本文方法仍有不足之處，對于需要添加支撐的區(qū)域仍不能實(shí)現(xiàn)最小化，可能生成了一些額外的支撐。因此后續(xù)的研究重點(diǎn)是對區(qū)域的優(yōu)化分析。

[1]Chalasani K，Jones L，Roscoe L. Support Generation for Fused Deposition Modeling[C]//Proceedings of Solid Freeform Fabrication Symposium. Austin：Texas，USA， 1995：229-241.

[2]Huang Xiaomao，Ye Chunsheng，Mo Jianhua，et al. Slice Data Based Support Generation Algorithm for Fused Deposition Modeling[J]. Tsinghua Science and Technology，2009，14(S1)：223-228.

[3]Lan Po-Ting，Chou Shuo-Yan，Chen Lin-lin，et al. Determining Fabrication Orientation for Rapid Prototyping with Stereolithography Apparatus[J]. Computer-Aided Design，1997，29(1)：53-62.

[4]Alexander P，Allen S，Dutta D. Part Orientation and Build Cost Determination in Layered Manufacturing[J]. Computer-Aided Design，1998，30(5)：343-356.

[5]Strano G，Hao L，Everson R，et al. A New Approach to the Design and Optimization of Support Structures in Additive Manufacturing[J]. Advanced Manufacturing Technology，2013，66(9/12)：1247-1254.

[6]Kritchman E，Gothait H，Miller G. System and Method for Printing and Supporting Three Dimensional Objects：US，7364686[P]. 2008-04-29.

[7]陳巖，王士瑋，楊周旺，等. FDM三維打印的支撐結(jié)構(gòu)的設(shè)計(jì)算法[J]. 中國科學(xué)，2015，45(2)：259-269.

Chen Yan，Wang Shiwei，Yang Zhouwang，et al. Construction of Support Structure for FDM 3D Printers[J]. Scientia Sinica，2015，45(2)：259-269.

[8]Huang Xiaomao，Ye Chunsheng，Wu Siyu，et al. Sloping Wall Structure Support Generation for Fused Deposition Modeling[J]. Advanced Manufacturing Technology，2009，42(11/12)：1074-1081.

[9]Dumas J，Hergel J，Lefebvre S. Bridging the Gap：Automated Steady Scaffolding for 3d Printing[J]. ACM Transactions on Graphics，2014，33(4)：1-10.

[10]卞宏友，劉偉軍，王天然. 基于STL模型支撐算法的研究[J]. 機(jī)械設(shè)計(jì)與制造，2005(7)：49-52.

Bian Hongyou，Liu Weijun，Wang Tianran. Study on Algorithm for Generate Support Based on STL File[J]. Machinery Design & Manufacture，2005(7):49-52.

[11]洪軍，王崴，張宇紅. 光固化快速成形自動支撐技術(shù)研究[J]. 機(jī)械工程學(xué)報(bào)，2004，40(11)：134-138.

Hong Jun，Wang Wei，Zhang Yuhong，et al. Research on Auto-support in Stereolithography[J]. Journal of Mechanical Engineering，2004，40(11):134-138.

[12]Kazhdan M，Bolitho M，Hoppe H. Poisson Surface Reconstruction[C]//Proceedings of the fourth Eurographics Symposium on Geometry Processing. Aire-la-Ville，Switzerland：Eurographics Association，2006：61-70.

(編輯張洋)

Generation of Branching Support Structures Based on Critical Angle Constraint

Shen ZhenhongDai NingLi DaweiWu Changyou

Nanjing University of Aeronautics & Astronautics，Nanjing，210016

In order to optimize support structures’ volume, and save the material and printing time, an automatic generation algorithm was proposed based on critical angle constraint. Firstly, the algorithm detected areas that required support, and the support points were obtained by adaptive sampling based on the area. Secondly, through critical angle constraint, the nodes in branching support structures were calculated, and all nodes information in a tree were saved. Finally, the resulting structure was obtained by the sweep surface after Poisson reconstruction. The algorithm was evaluated by experiments, and the experimental results demonstrate that the algorithm provides a reliable method to ensure the printing stability as well as reducing the printing material and time.

3D printing; critical angle constraint；branching support；Poisson surface reconstruction

沈振宏，男，1990年生。南京航空航天大學(xué)機(jī)電學(xué)院碩士研究生。主要研究方向數(shù)字化設(shè)計(jì)與制造、三維打印技術(shù)。戴寧(通信作者)，男，1978年生。南京航空航天大學(xué)機(jī)電學(xué)院副教授。李大偉，男，1989年生。南京航空航天大學(xué)機(jī)電學(xué)院碩士研究生。吳長友 ，男，1989年生。南京航空航天大學(xué)機(jī)電學(xué)院碩士研究生。

2015-06-15

國家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)資助項(xiàng)目(SS2013AA040802)；江蘇省科技支撐項(xiàng)目(BE2014009-3)；航空科學(xué)基金資助項(xiàng)目(20151652024)；江蘇省三維打印裝備與制造重點(diǎn)實(shí)驗(yàn)室開放課題資助項(xiàng)目(BM2013006)

TH164

10.3969/j.issn.1004-132X.2016.08.020