基于k-nearest算法的增強型正交各向異性多孔結構設計

曹 偉，蔡和倫，謝卓尊，劉 斌，江開勇

(華僑大學 福建省特種能場制造重點實驗室/廈門市數字化視覺測量重點實驗室，福建 廈門 361021)

0 引言

多孔結構具有比強度高、比表面積大、輕量節材、饋能減振等優點，被廣泛應用于航空航天、汽車工業、生物醫療等眾多領域。近年來，隨著增材制造技術的飛速發展，使得具有復雜內部多孔結構的零件制造成為可能，因此多孔結構的輕量化、定制化設計和建模逐漸成為當前增材制造技術最具挑戰性的研究方向之一，主要表現為建模過程復雜、數據量龐大、計算和仿真成本較高等。此外，對于更加復雜的功能梯度多孔結構(Functionally Gradient Porous Materials, FGPMs)，要求根據零件局部應力大小和方向設計其內部結構，在保證幾何連續的情況下滿足梯度力學性能，極大地增加了設計復雜度。

根據建模方法的不同，多孔結構可以分為規則多孔結構和不規則的隨機多孔結構。由于設計方便、可控性高，早期的多孔結構大多基于規則胞元的晶格結構，即采用規則胞元陣列填充設計空間，主要包括蜂窩結構、菱形晶格結構、立方體結構、三周期極小曲面等[1-3]。蜂窩結構是一種由眾多柱狀孔穴胞元密鋪組成的多孔結構，這種結構在節省材料的同時，保證了較高的比強度與比剛度。邢昊等[4]針對二維蜂窩結構在共面方向集中載荷作用下承載性能不足的問題,提出一種基于拓撲優化密度映射的非均勻蜂窩結構建模方法。針對多孔結構內部支撐難以去除的問題，WU等[5]提出一種自支撐的菱形晶格結構，并用該結構填充零件內部，通過調節不同晶格的壁厚實現力學性能的梯度分布;WANG等[6]提出一種改進的面心立方多孔結構設計方法，通過選擇性激光融化(Selective Laser Melting, SLM)工藝打印后研究其抗沖擊行為。三周期極小曲面(Triply-Periodic Minimal Surface, TPMS)是一種數學上定義的空間結構，整個結構沒有自交點，具有無限延伸表面和復雜的開孔結構[7]。REN等[8]提出一種多尺度混合TPMS優化設計模型，該模型能較好地解決不同尺度下TPMS間邊界的連續性問題，并建立了應力與TPMS之間的映射關系。此外，在上述結構的基礎之上，還衍生出一些變異的晶格結構，例如對稱雙金字塔型的十二面體結構、星型結構等，通過改變連接桿的直徑或面的厚度實現力學性能的梯度分布[9-10]。雖然規則晶格多孔結構能夠降低設計復雜度、減少設計周期，但是由于其固有特性，其設計自由度較低，對復雜模型的適應性較差，難以進行各向異性設計，可能造成幾何結構、力學性能的不連續，限制了其使用范圍。

隨機多孔結構是一種不規則多孔結構，它不僅適用復雜模型，設計自由度高，還能完美解決規則多孔結構幾何、力學性能不連續等問題。典型的隨機多孔結構包括基于拓撲優化的多孔結構、基于維諾圖的多孔結構等。拓撲優化是一種根據給定約束，在設計域內對材料分布進行優化的方法[11]。DAYNES等[12]提出一種基于拓撲優化的功能梯度多孔結構建模方法，并通過實驗證明其強度和剛度遠高于同等質量的規則立方體晶格結構。然而，拓撲優化方法大多用于優化二維結構，對三維模型的適用性較差，無法滿足對復雜三維結構建模的需求。維諾多孔結構是一種模擬自然界物體的理想結構，如骨骼、海綿等。MARTNEZ等[13]提出一種基于直桿的三維維諾多孔結構的高效設計方法。在此基礎上，LIU等[14-15]將三維維諾圖與隱式曲面相結合，設計了一種全連通的光滑多孔結構，有效避免了應力集中。雖然三維維諾多孔結構可控性強、結構連續，但建模復雜度高，且難以實現各向異性多孔結構的建模。

1 增強型正交各向異性多孔結構建模

1.1 建模流程

本文提出的基于改進k-nearest算法的增強型正交各向異性多孔結構的設計流程如圖1所示，其初始輸入條件包括站點密度場ρ和方向場φ。首先，采用站點隨機采樣算法將密度場映射為設計域內的站點分布，控制正交各向異性多孔結構力學性能的梯度分布；在設計域內疊加方向場，將其映射為黎曼空間中的流形場，實現多孔結構局部正交各向異性的任意旋轉。需要說明的是，初始的密度場ρ和方向場φ可以從模型有限元分析的應力場提取，也可以交互式輸入。其次，計算鄰域內站點間的黎曼距離，采用改進的k-nearest算法連接相鄰的k(k≥2)個站點，生成2D多孔結構模型；最后，對生成的2D多孔結構在厚度方向進行拉伸，可以得到2.5D多孔結構模型，并最終通過3D打印機打印成型。

1.2 站點隨機采樣

站點的隨機采樣是建立多孔結構的基礎，本文提出一種密度場驅動的站點隨機采樣算法，先將設計域離散化為一系列正方形柵格，然后通過給定的密度場函數ρ在每一個網格中進行隨機采樣。以二維設計域為例，站點隨機采樣過程如圖2所示，具體算法如下：

步驟1將二維設計域離散為邊長為a的正方形柵格。

步驟2根據站點密度函數場ρ，計算柵格i內的站點數量ti=a2ρi。

步驟3將柵格i等分為4個正方形子柵格。

步驟4若柵格i內的站點數量ti≤4，則隨機選取ti個子柵格，并分別在每個子柵格內部任意位置隨機生成一個站點。

步驟5若ti>4，則轉步驟2，重新計算子柵格中的站點數，劃分柵格，并隨機生成站點(如圖2)。重復步驟2～步驟5，直至所有區域的站點分配完成。

對于3D模型，其站點隨機采樣算法和上述2D模型的算法類似，不同之處在于：需要將3D模型對應的三維設計域離散化為邊長為a的立方體柵格，然后將每個柵格平均劃分為8個邊張相等的子柵格，一般采用八叉樹結構表示。

從圖2可以看出，設計域內站點的分布是由密度場ρ驅動的，是多孔結構實現彈性性能梯度分布的根本因素，其原因在于站點分布的特征決定了材料局部分布，從而直接影響其局部力學性能。即站點分布密集的區域，材料分配相對較多，力學性能相對較高；相反地，站點分布稀疏的區域，材料分配較少，力學性能相對較低。因此，該算法可以通過站點密度分布函數從宏觀上控制多孔結構不同部位的力學性能，在應力較大的部位設置較大的站點密度，增加力學性能；反之，在應力較小的部位設置較小的站點密度，節約材料。此外，在微觀上又能保持局部隨機性，從而保證梯度多孔結構的幾何連續性。

隨機采樣算法中輸入的密度場函數ρ可以由用戶指定，也可以從模型的有限元分析得到，如圖3所示。

1.3 方向場

方向場主要用于控制正交各向異性多孔結構正交軸的旋轉方向。對于設計域中任意點鄰域內的多孔結構單元，其強度最高的正交軸方向始終與方向場保持一致，從而增強多孔結構的力學性能，降低材料消耗。方向場可以有多種形式，如場函數、流形場、應力方向場等。

1.4 改進的k-nearest算法

將方向場映射為黎曼空間中的流形場，采用黎曼距離計算站點間的“距離”，實現任意局部區域各向異性方向的旋轉，從而建立起具有梯度彈性性能的多孔結構。為方便理解，本節以二維黎曼空間為例闡述各向異性多孔結構的建模過程，且可以較為容易地拓展到三維黎曼空間。

1.4.1 黎曼度量

設點p是二維黎曼空間Ω中的一點，p點處的度規張量矩陣M(p)可視為一個正定矩陣，即

(1)

M(p)=ETUE,

(2)

(3)

(4)

定義1設坐標點p與q是黎曼空間Ω中任意兩點，M為p點對應的黎曼度量，則從p點到q點的有向距離dp(p,q)定義為：

(5)

由于M(p)≠M(q)，在黎曼空間中點p與點q的距離是非對稱的，會導致p到q的有向距離與q到p的有向距離在數值上不一致，即dp(p,q)≠dq(q,p)，因此本文將兩點間的“距離”定義為：

(6)

1.4.2k-nearest算法

如圖5所示為采用文獻[16]中的k-nearest算法設計的正交各向異性多孔結構：依次遍歷所有站點，對每個站點采用直線連接與其“距離”最近的k(k=2,3,4,…)個站點，并賦予每條直線寬度t(即2.5D多孔結構的壁厚，t≥2τ,τ為打印機的最小打印尺寸，如打印機噴頭直徑等)和厚度h，即可生成最終的2.5D正交各向異性多孔結構。

1.4.3 算法改進

從圖5可以看出，采用k-nearest算法生成的各向異性多孔結構存在兩個明顯缺陷：

(1)在站點分布相對稀疏的區域會沿方向場生成較大的孔洞，這種孔洞比期望的孔洞大一個量級，如圖5b中橙色圈所示。

(2)在站點分布相對密集的區域，多條線段會產生交叉現象。交叉點相當于增加了站點數量，導致實際站點密度和期望的站點密度ρi存在較大誤差，如圖5b中藍色圈中所示。

這種不均勻的孔洞分布以及新增的站點，不僅使多孔結構不符合設計要求，還會導致材料分布不均勻，產生應力集中，嚴重影響多孔結構的整體力學性能。此外，在后續模型3D打印成型時受限于打印機的精度，過于細小的孔洞(r≤2τ)幾乎無法打印。經過分析發現，產生上述兩種缺陷的根本原因是由于算法沒有考慮站點連線相交的情況。

基于此，本文提出一種改進的k-nearest算法，可以判斷并去除會產生相交的線段，確保每個站點均連接k個相對“最近”的站點。具體算法如下：

步驟1在歐式空間中，計算當前站點ni的鄰域，將鄰域內的站點存入存儲器Neighbour中。對2D/3D設計域，站點ni的鄰域半徑rni分別為：

(7)

(8)

其中：a為設計域內劃分柵格的邊長;hui、hvi和hwi分別為ni對應的黎曼橢圓/橢球各軸的長度。該方法保證鄰域內存在滿足條件的k個站點的同時，又避免了在黎曼空間中進行全局搜索，提高了算法效率。

步驟2在黎曼空間中，按式(6)計算ni與鄰域(Neighbour)內站點的“距離”，并從小到大排序。

步驟3依次從排序后的Neighbour取出站點，判斷ni與該點的連線與已生成的結構M是否相交。若不相交，則用線段連接兩站點，并更新M；若相交，則放棄連接。重復步驟3，直至與ni連接的站點數等于k為止。

步驟4遍歷設計域內所有站點，重復步驟1～步驟3。

圖6中上下兩欄分別為采用文獻[16]中的k-nearest算法和本文算法在站點分布相同、不同的各向異性度量hx∶hy(x和y為正交軸)下生成的2.5D多孔結構。不難發現，本文算法能夠明顯地減少上述兩種缺陷，使材料分布更均勻，從而增強了多孔結構的力學性能。

本文算法也能比較容易地拓展至三維空間，區別在于需要用直徑為d的圓柱體支桿代替二維空間中的線段，判斷當前生成的支桿是否會與已生成的結構M產生干涉，并采用布爾運算對新支桿與M的三角網格進行求交。算法偽代碼如表1所示，生成的三維多孔結構如圖7所示?？梢钥闯?，通過調整hu、hv和hw的比例，可以非常方便地控制多孔結構在不同正交軸向的各向異性程度。

表1 三維正交各向異性多孔結構生成算法

1.4.4k值的選取

本文算法的核心是將站點p與其相對“最近”的k個站點相連，因此，k值的選取十分關鍵，對多孔結構的連通性和力學性能都有較大的影響。圖8所示為不同k值對應的二維結構，可以發現：當k=2,3時，幾何連通性很差，無法生成多孔結構；當k=4,5時，雖然能夠生成多孔結構，但是仍存在局部區域不連通，影響了多孔結構結構的力學性能；當k≥7時，雖然結構幾何連通性較好，但過于緊湊，降低了多孔結構的設計空間。而BALISTER等[17]已經證明當k=6時候，生成的結構最優，即具有相對較好的幾何連通性和較大的設計空間。因此，在默認情況下本文后續章節中所建立的多孔結構其k值均為6。

2 多孔結構的性能分析

2.1 各向異性多孔結構的數值均質化

數值均質化是面向周期性或近周期性材料宏觀性能的有限元分析方法，廣泛應用于復合材料研究。由于本文算法中的站點生成器是均勻分布的隨機序列，所生成的站點具有隨機性，從概率意義上，當站點密集時，即模型宏觀尺寸相對于平均孔徑足夠大時，模型整體結構呈現近周期性。因此，數值均質化方法仍然適用于本文多孔結構的分析[15-16]。這種達到近周期性體積的多孔結構模型被稱作代表性體積單元(Representative Volume Element, RVE)，是數值均質化的分析對象。

數值均質化能有效地確定周期性復合材料的彈性張量。設材料單元為V，體積為|V|，在材料彈性變形階段，應力σ與應變ε成線性關系[18]：

σ=Cε。

(9)

其中C是彈性矩陣，用于表征材料的彈性行為。

為方便編程實現，本文采用ANDREASSEN等[19]提出的方法，將多孔結構模型劃分為正方形網格，對網格單元節點編號，并施加周期性邊界條件，通過線彈性方程kχ=f求解出χij，可得到均質化彈性張量為：

(10)

其中:i、j為對應矩陣的第i、j行/列;k為剛度矩陣;f和χ分別是等效節點的載荷列陣與位移列陣；ke為單元剛度矩陣;χe表示對應網格單元節點的整體位移列陣，可以從χ的對應節點中提??；Ve表示對應網格單元的體積;N為RVE單元劃分的有限網格單元總數量。需要說明的是，上式僅適用于描述多孔結構材料的彈性變形行為，而無法描述材料達到屈服極限后的塑性變形或斷裂行為。

2.2 彈性性能分析

2.2.1 數值均質化分析

為分析本文正交各向異性多孔結構的彈性性能，對圖6中k-nearest算法和本文算法生成的4組2.5D多孔結構進行數值均質化分析。所有結構的尺寸均為50×50×7 mm3，壁厚t=0.4 mm。每組模型的站點分布相同，且站點是均勻隨機分布的。4組多孔結構的各向異性度量hx∶hy分別為1∶1、1∶3、1∶5和1:7，其中x、y為正交軸。從圖中可以看出，當hx∶hy=1∶1時，所生成的多孔結構呈各向同性；當hx∶hy≠1時，對應的多孔結構呈明顯的正交各向異性，且比例越大(或越小)，各向異性程度越高。

(11)

(12)

表2 不同多孔結構的彈性性能分析與對比

由表2可以看出，本文算法生成的多孔結構的相對密度比k-nearest算法生成的有所增加，其增量分別為16.2%、14.3%、17.5%，平均增加16.0%。主要原因是k-nearest算法連接的是絕對“最近”的k個站點，而本文算法為避免產生相交，連接的是相對“近”的站點，因此桿的平均長度要長于k-nearest算法，導致相對密度增大。雖然本文算法生成的多孔結構相對密度較大，但是其彈性性能相對于k-nearest算法大幅增強，以hx∶hy=1∶3為例，本文算法生成的多孔結構在x、y方向上的相對彈性模量分別增強了133.3%和26.9%。

2.2.2 壓縮實驗

為進一步研究上述正交各向異性多孔結構的力學性能，驗證數值均質化結果的準確性，將圖6中的4組多孔結構采用熔融沉積成型(Fused Deposition Modeling, FDM)工藝進行3D打印成型，打印材料選擇為聚乳酸(Polylactic Acid, PLA)，基體彈性模量為1 200 MPa。每種多孔結構打印2組，分別沿x、y方向進行壓縮實驗。實驗設備為萬測電子萬能試驗機(型號：TSE504D)，速率設置為1 mm/min，實驗結果如圖9所示。

從圖9b中各組試樣的力—位移曲線可以看出，在多孔結構發生塑性變形之前，本文算法生成的多孔結構能夠承受更大的載荷；隨著多孔結構各向異性程度增加，y方向能承受的最大載荷值逐漸增大，同時x方向能承受的最大載荷值逐漸降低，該結果和上一節數值均質化的結果保持一致。

此外，通過壓縮實驗，還可以測得試樣的彈性模量。表3所示為壓縮實驗測得的試樣彈性模量與數值均質化計算的相對彈性模量之間的對比情況?？梢钥闯?，對不同拉伸比例下的2.5D各向異性多孔結構，數值均值化方法相對于壓縮實驗結果的誤差較小，最大相對誤差僅為9.38%，表明數值均質化方法能較為準確的預測本文提出的各向異性多孔結構的彈性性能。

表3 壓縮試驗與均質化結果對比

對于三維各向異性多孔結構，采用光固化(Stereolithography，SLA)工藝打印了3組試樣，材料為光敏樹脂，試樣在沿三個正交軸方向的拉伸比例分別為1∶1∶1、1∶3∶3和1∶1∶3，尺寸均為40×40×40 mm3，桿直徑為d=0.4 mm。對3組試樣分別沿著x、y、z方向進行壓縮實驗，實驗結果如圖10和表4所示。

表4 三維各向異性多孔結構壓縮實驗結果 MPa

實驗結果表明，隨著y、z方向各向異性程度的增加，對應方向的彈性性能也相應提升，而其他方向上的彈性性能逐漸降低。同時，從1∶3∶3和1∶1∶3樣品的實驗結果還可以看出，僅對一個方向進行拉伸的多孔結構的彈性性能遠高于在兩個方向同時拉伸的多孔結構。

2.2.3 失效分析

通過上述壓縮實驗，測得不同多孔結構在y方向的抗壓強度如表5所示?？梢钥闯?，隨著y方向拉伸比例的增加，兩種算法生成的多孔結構在y方向的抗壓強度逐漸增大，但是本文算法生成多孔結構的抗壓強度相比于k-nearest算法有明顯的增強，增強幅度均達到200%以上。

表5 各向異性多孔結構的抗壓強度

在壓縮實驗中，當試樣承受的壓應力超過抗壓強度后，試樣會發生局部屈服或破裂，進而發生失效，如表6所示(壓縮12%)?？梢钥闯?，在應變較大時，hx∶hy=1∶1的各向同性多孔結構未表現出明顯的失效，且本文結構要優于k-nearest算法生成的結構；而1∶3、1∶5和1∶7的各向異性多孔結構均會發生不同程度的失效，如桿狀結構發生較大彎曲變形甚至斷裂，且失效大多發生在孔洞較大的位置，與1.4.3小節的分析一致。其原因主要是k-nearest算法會導致材料分布不均勻，在孔洞較大的位置材料分布較少，從而產生應力集中。雖然本文設計的增強型多孔結構能在一定程度上使材料分布更均勻，提高了力學性能，但是并不能從根本上解決這一問題，當多孔結構各向異性程度增大時，材料分布的不均勻性也隨之增大。

表6 各向異性多孔結構的失效分析(壓縮12%)

2.2.4 隨機性分析

由于站點是在密度場驅動下隨機采樣的，這種隨機性可能會對多孔結構的彈性性能產生影響，造成性能的不穩定。為研究隨機性的影響程度，在保持站點密度ρ不變的情況下，對hx∶hy=1∶1、1∶3、1∶5和1∶7四種拉伸比例各隨機生成了3組多孔結構，尺寸為50×50×7 mm3，壁厚t=0.4 mm，如表7所示。從模型外觀上看，站點分布的隨機性對多孔結構的影響不大。

表7 隨機生成的站點與對應的多孔結構

表8 站點隨機性對多孔結構彈性性能的影響

續表8

2.2.5 彈性設計空間

為了研究本文提出的各向異性多孔結構的彈性設計空間，本文隨機生成了大量相同尺寸的2.5D多孔結構(約2 400個)，其拉伸比例hx∶hy范圍選取為10∶1到1∶10，旋轉角度θ=0。對這些多孔結構逐一進行數值均質化分析，將其結果繪制成散點圖如圖11所示?？梢钥闯?，各向異性多孔結構對應的點均勻分布在對角線的兩側，且呈對稱分布，即hx∶hy>1對應的多孔結構分布在對角線的下方，而hx∶hy<1對應的多孔結構分布在對角線上方。同時，隨著拉伸比例的增大，各向異性程度增大，對應點距離對角線越遠。當hx=hy時，多孔結構呈各向同性，對應點分布在對角線上。

圖11中散點覆蓋的區域即為本文所提出的多孔結構(2.5D)的彈性設計空間，在多孔結構壁厚為常量時，其彈性性能主要受兩個設計參數的影響，即拉伸比例hx∶hy和站點密度ρ。其中，hx∶hy與多孔結構的各向異性呈正相關，ρ與多孔結構的宏觀彈性性能呈正相關。需要說明的是，圖11所示的多孔結構設計空間僅適用于線彈性階段，多孔結構發生塑性變形時的力學性能并不在本文的討論范圍之內。

2.3 案例

本節主要討論正交各向異性多孔結構的應用情況，建模過程均采用Visual Studio 2013的MFC框架，搭配CUDA程序并行計算，在硬件環境為Inter Core2 i5-10400 CPU，16 GB內存與NVIDIA GeForce GTX 1650顯卡的計算機上實現，最后通過3D打印機打印成型。

2.3.1 應力驅動的梯度多孔結構

如圖12所示為采用本文方法設計應力驅動的梯度各向異性多孔梁結構的建模過程。首先對梁結構進行有限元分析，得到如圖12a所示的應力場。將應力大小映射為站點密度場，通過1.2節的隨機采樣算法可以生成如圖12b所示的站點分布；同時，將應力方向映射為如圖12c所示的方向場。最后采用1.4節提出的改進k-nearest算法生成如圖12d所示的2.5D正交各向異性多孔結構。梁模型的尺寸為250×50×30 mm3，壁厚t=1 mm，體積分數為43%。為對比分析，本文還設計了如圖12e所示的蜂窩結構，其尺寸、壁厚均與圖12d中的多孔結構相同。此外，通過調整蜂窩結構的胞元大小，使其體積分數也為43%。

采用SLA3D打印機將上述兩種結構打印成型，打印材料為光敏樹脂。分別對其做三點彎曲破壞性實驗，如圖13a所示，實驗設備為萬測電子萬能試驗機(型號：TSE504D)，選用彎曲夾具，速率設置為1 mm/min。實驗測得的力—位移曲線如圖13b所示，其他測試結果如表9所示。

表9 兩種結構的三點彎曲實驗結果比較

項目重量/g最大載荷/kN彎曲模量/MPa彎曲強度/MPa蜂窩結構1801.46828 258.1880.1各向異性多孔結構1782.35366 454.81 413.3提升率/%-1.1160.29135.1760.58

不難看出，在兩種結構的重量在幾乎相等的條件下，相比于規則的蜂窩結構，本文方法設計的應力驅動的梯度多孔結構的力學性能大幅提升，最大載荷、彎曲模量和彎曲強度分別提升了60.29%、135.17%和60.58%。這是由于本文方法能夠通過站點密度分布調整局部材料的分布，在應力較大的部位分配較多材料，應力較少的部位減少材料分配，實現在相同材料用量的情況下，提升多孔結構力學性能；或在承受相同載荷的情況下，降低材料的用量，達到綠色節材的目的。

需要說明的是，由于圖12b中的站點分布是由設計者根據應力場和設計經驗綜合設置的密度場函數隨機采樣得到的，由站點分布驅動的梯度正交各向異性多孔結構只是一個較優的結構，還存在一定的優化空間。若要獲得最優的多孔結構，使材料分布更符合應力分布，則需要通過實驗、仿真或機器學習等方法建立應力場和站點密度場之間定量的映射模型，從應力場得到最優站點分布，疊加方向場后生成最優的正交各向異性多孔結構，這也是后續工作需要進一步研究的內容。此外，由于站點是隨機采樣得到的，對外觀對稱的模型無法實現內部多孔結構的對稱(如圖12d內部多孔結構左右不對稱)，從而會在一定程度影響其美學特性。

表10列出了k-nearest算法和本文算法在2.5D機械零件輕量化設計應用情況的對比分析。通過對零件受力情況進行有限元分析，可以得到零件的應力分布。將應力大小映射為設計域內的站點分布，應力方向映射為方向場，分別采用k-nearest算法和本文算法生成了兩組多孔結構?？梢钥闯?，兩種算法生成的多孔結構均能在滿足材料屈服極限的條件下實現零件的輕量化，但本文算法生成的多孔結構的最大Mises應力明顯小于k-nearest算法，說明本文算法改善了k-nearest算法材料分布不均勻、易產生應力集中的現象，增強了多孔結構的力學性能。此外，若后續工作中能建立應力場和站點密度之間的映射模型，則可以進一步對多孔結構進行優化，實現零件最佳輕量化的目標。

表10 本文算法在機械零件輕量化設計的應用情況

2.3.2 其他應用

本文正交各向異性多孔結構的設計方法還可以應用于其他場景，如圖14所示。圖中模型均采用FDM工藝打印，材料為PLA。其中，圖14a為2.5D椅子模型，通過設置方向場，可以實現座位、靠背和椅腿等部位的正交各向異性。圖14b為3D牙齒模型，其內部被3D正交各向異性多孔結構填充。根據牙齒的受力情況，可以發現多孔結構沿豎直方向和徑向呈正交各向異性，在實現牙齒模型輕量化的同時，又能保證在上述兩個方向上的強度。圖14c為采用正交各向異性多孔結構填充的3D車輪模型，多孔結構在輪轂部分沿徑向拉伸，使輪轂能夠承受較強的徑向力；在輪胎部分沿周向拉伸，使輪胎相對柔軟，從而在運動時實現緩震的功能。

3 結束語

本文提出一種基于k-nearest算法的增強型正交各向異性多孔結構設計方法，通過站點密度場驅動站點的隨機采樣，通過改進k-nearest算法，疊加全域方向場，生成增強型正交各向異性多孔結構，并通過3D打印機將多孔結構打印成型。最后，分別通過數值均質化和壓縮實驗對打印的多孔結構進行力學性能、失效狀態和隨機性進行了分析和討論。實驗結果表明，采用本文方法設計的正交各向異性多孔結構具有更優良的彈性性能和更大的彈性設計空間。此外，本文提出的正交各向異性多孔結構設計方法可控性高，不僅可以通過黎曼橢圓/球的拉伸比例控制局部正交各向異性程度，通過方向場控制各向異性的朝向，還可以通過站點密度場調控多孔結構力學性能的梯度分布。

雖然多孔結構性能明顯增強，但本文方法在實際工程應用中仍存在一定的局限性，即并不能從根本上解決局部材料分布不均勻的情況，尤其是各向異性較大的區域。其次，由于本文未建立起應力分布與站點密度函數之間的數學映射關系，導致在設計時需要根據設計師經驗生成站點的隨機分布，無法實現設計過程的智能化和結果的最優化。后續工作包括：1)解決局部材料分布不均勻問題；2)建立應力場與多孔結構設計參數之間的數學模型，將應力大小映射為站點密度，將應力方向映射為多孔結構的方向場，從而建立應力驅動的梯度多孔結構模型。