劉繼凱 張乘虎 袁志玲 黃嘉奇

摘要：為了提高點陣材料結構件的力學性能，提出了一種基于Ordered SIMP方法的點陣-實體復合結構拓撲優(yōu)化方法.采用一種三維X型微結構單元作為點陣材料，通過數(shù)值擬合建立點陣材料相對密度與其等效物理屬性之間的函數(shù)關系.在宏觀結構拓撲優(yōu)化問題中，以點陣材料等效密度為設計變量，基于Ordered SIMP插值方法建立點陣材料和實體材料相結合的多材料插值模型，進而，以材料體積用量為約束，以結構柔度最小化為目標實現(xiàn)點陣-實體復合結構的多尺度拓撲優(yōu)化設計.通過數(shù)值算例和實驗測試表明，相比于僅使用點陣材料填充的設計，本方法能夠獲得更好的結構剛度.

關鍵詞：結構優(yōu)化;點陣結構;拓撲優(yōu)化;多尺度設計;Ordered SIMP

中圖分類號：TH122

文獻標志碼：A

收稿日期：2021-07-06

基金項目：江蘇省自然科學基金資助項目（BK20190198），Natural Science Foundation of Jiangsu Province（BK20190198）;山東省自然科學基金資助項目（ZR2020QE165），Natural Science Foundation of Shandong Province（ZR2020QE165）;發(fā)動機可靠性國家重點實驗室（skler-202001），State Key Laboratory of Engine Reliability（skler-202001）

作者簡介：劉繼凱（1987—），男，山東濰坊人，山東大學教授，博士研究生導師

?通信聯(lián)系人，E-mail：jikai_liu@sdu.edu.cn

拓撲優(yōu)化（Topology optimization）是以結構設計域內(nèi)的材料分布為研究對象，在滿足給定的邊界以及約束條件下，尋找設計域內(nèi)最優(yōu)材料分布以實現(xiàn)結構性能達到最優(yōu)的結構設計方法[1].而與實體材料相比，點陣材料作為一種先進新型輕質(zhì)材料，具有超輕、高能量吸收性、高孔隙率等功能屬性[2-6].因此，在拓撲優(yōu)化結構設計中采用點陣材料可以獲得更豐富的性能.

近年來，基于變密度點陣材料的多尺度結構拓撲優(yōu)化方法得到了廣泛的研究[7-9].Cheng等[10]提出了一種基于應力約束的變密度點陣材料結構拓撲優(yōu)化方法，通過實驗驗證了方法的有效性.Wang等[11]提出了一種控制點陣材料微觀結構中多個幾何參數(shù)的點陣材料拓撲優(yōu)化方法，相比于單一幾何參數(shù)變量控制的優(yōu)化方法，該設計方法優(yōu)化得到的結構性能具有明顯提升.Liu等[12]提出了一種多類點陣材料填充的拓撲優(yōu)化方法，每個單元中點陣材料的分布考慮兩個設計變量，即點陣材料類型和它的幾何參數(shù)，有效地擴大了結構的設計自由度.這些通過參數(shù)化表征點陣材料等效力學性質(zhì)的多尺度拓撲優(yōu)化方法具有很高的計算效率.然而，點陣材料的剛度和強度與實體材料相比仍有很大差距.因此，很多情況下僅使用點陣材料填充的設計無法提供足夠的力學性能.

為解決上述問題，本文基于Ordered SIMP方法[13]，建立了點陣-實體多材料插值模型，提出一種全新的點陣-實體復合結構拓撲優(yōu)化方法.將其應用于簡支梁的設計，并通過數(shù)值算例和實驗測試驗證了該方法的有效性.

1點陣-實體復合結構拓撲優(yōu)化模型

1.1點陣材料等效力學性能計算基于多尺度設計理念，將宏觀結構離散為若干

單元，每個單元代表一種待設計的微結構胞元，且微結構胞元的等效力學屬性可基于能量均勻化方法求得[14-18].在變密度法拓撲優(yōu)化設計中，微結構胞元的相對密度ρL表示為微觀尺度上微結構胞元實體部分體積相對于設計域體積之比，如式（1）所示：

式中：vL為微結構胞元實體部分的體積;v為微結構胞元的設計域體積.

本文中，采用一種三維X型微結構作為點陣材料胞元（如圖1），它具備抗拉壓和抗剪切能力，以及增材制造加工過程中的自支撐優(yōu)勢.通過改變結構中桿的直徑，點陣材料的相對密度隨之改變.點陣材料的等效彈性矩陣DH可表示為：

在不失一般性的前提下，為了便于計算式（2）中等效彈性矩陣DH，假設點陣材料的基材料彈性模量為E0=100MPa，泊松比為μ=0.3.基于能量均勻化方法[14-18]，可計算得到點陣材料的彈性矩陣，如圖2所示.顯然，為了更優(yōu)的結構性能，在優(yōu)化設計過程中需要設計一系列相對密度的點陣材料.然而，通過均勻化方法求解每個微觀結構彈性矩陣的方法相對耗時.為獲取一系列相對密度的點陣材料的彈性矩陣，通過采集一定數(shù)目樣本點數(shù)據(jù)，使用擬合的方法，建立點陣材料相對密度及其等效彈性矩陣的數(shù)學關系式.考慮擬合精度與計算效率，建立5次多項式函數(shù)作為擬合關系式：

式中：dij為彈性矩陣DH中的元素.

考慮點陣材料的可制造性，本文設計過程中點陣材料的相對密度變化范圍是15%～70%，在相對密度變化范圍內(nèi)取樣，并計算對應點陣材料的彈性矩陣.其彈性矩陣中有3個相互獨立的的參數(shù)，分別對應3個擬合函數(shù).基于24個樣本點，使用最小二乘法，得到擬合函數(shù)的相關參數(shù)如表1所示，擬合效果

如圖3所示.

1.2點陣-實體多材料插值模型

針對一定體積約束下的柔度最小化問題，在點陣-實體復合結構拓撲優(yōu)化設計方法中主要有以下兩種方案.其一，首先基于SIMP（solid isotropic mate?rial with penalization）方法[1]對結構進行小于體積分數(shù)約束的宏觀拓撲優(yōu)化設計，并將優(yōu)化得到的結果定義為實體.之后，在宏觀設計域進行剩余體積分數(shù)約束下的點陣材料填充的拓撲優(yōu)化設計.該方法的局限性在于將點陣-實體復合結構的拓撲優(yōu)化設計過程分離，造成設計結果偏離最優(yōu)解.另一種設計方案是將點陣材料的密度上限提高到實體材料密度，從而將實體材料加入到等效材料本構參數(shù)的擬合樣本中.然而，這種處理方法得到的優(yōu)化結果中會引入大量的高密度點陣材料，即密度接近于1而未形成完整實體的單胞，這部分高密度點陣材料單胞會形成明顯的封閉空腔，可制造性差.為了解決以上設計方案的不足，本節(jié)基于Ordered SIMP插值模型，建立了一套點陣-實體多材料插值模型.D4902AB4-8349-4BC7-B388-7D7244E450B4

多材料插值模型中的自變量是點陣材料和實體材料的相對密度，它們具有統(tǒng)一性，無需歸一化處理.因此，只需歸一化處理點陣材料和實體材料的物理屬性，作為多材料插值模型中的因變量.多材料插值模型可表示為：

式中：dAij為歸一化處理的材料彈性矩陣DA中的元素;d0ij為實體材料的彈性矩陣D0中的元素;ρe為單元密度;P為懲罰因子;ρLmin、ρLmax分別為點陣材料的密度下限和上限;ρSmin、ρSmax分別為實體材料的密度下限和上限;x1、x2分別為實體材料彈性矩陣歸一化系數(shù)，具體見式（5）和（6）.

根據(jù)上述多材料插值方法，單元密度ρe的取值范圍是[0.15～1].其中，點陣材料的密度設計上限即為實體材料的密度設計下限，點陣材料的相對密度取值范圍是[0.15～0.7]，實體材料的密度取值范圍是（0.7～1].在優(yōu)化迭代過程中，可通過施加密度懲罰的方式消除實體材料插值引起的中間密度問題.因此，建立點陣-實體多材料插值模型如圖4所示.

1.3結構柔度問題的拓撲優(yōu)化模型

本文以典型拓撲優(yōu)化中柔度最小化問題為例，展示所提出方法的求解流程，該問題的具體數(shù)學表達如下：

式中：x?為設計變量，包含每個單元的密度ρi;C為結構柔度;F為結構所受載荷矢量;U為結構位移;K為結構總剛度陣;V*為優(yōu)化過程中結構總體積;v0為材料單元體積;V為結構總體積上限;ρmax、ρmin分別為單元密度變化范圍上限和下限.

根據(jù)所建立點陣-實體多材料插值模型，得到單元密度ρe與其彈性矩陣之間的映射函數(shù)關系，其單元剛度矩陣ke如下：

式中：Ωe為單元設計域;B為單元的應變矩陣;DA·D0為DA與D0的哈達瑪積（Hadamard product）.

將單元剛度陣組裝為結構總剛度陣K，并通過有限元方法計算獲得結構位移U.根據(jù)鏈式法則，求解拓撲優(yōu)化問題的靈∫敏度信息.其中，目標靈敏度為：

式中：ue為單元位移;n為設計變量數(shù)目.

體積靈敏度為：

根據(jù)上述求解過程中所獲得的系統(tǒng)狀態(tài)和靈敏度信息，使用移動漸近線方法（MMA）[19]更新設計變量，直至滿足收斂條件，從而結束整個優(yōu)化過程.本文所提出的點陣-實體復合結構拓撲優(yōu)化方法具體流程如圖5所示.

2實施案例

本章對簡支梁結構進行點陣-實體復合結構的拓撲優(yōu)化設計.通過與其它點陣材料填充的設計方法對比，充分展示了該設計方法的有效性.進一步，通過三點彎曲實驗驗證了本文所提出的拓撲優(yōu)化方法的優(yōu)越性.

2.1數(shù)值算例

如圖6所示，長度48mm，寬度8mm，厚度為14mm的簡支梁由頂層、中間層和底層組成.梁的兩端被支撐，其中間位置承受壓力大小為60N的負載.優(yōu)化設計的要求如下：

1）簡支梁的中間層為設計域，頂層和底層為非設計域;

2）體積約束為設計域的40%.

使用本文提出的點陣-實體復合結構拓撲優(yōu)化設計方法設計簡支梁結構.考慮優(yōu)化設計問題的對稱性，只需對設計域的1/2進行優(yōu)化設計.首先進行網(wǎng)格劃分，使用八節(jié)點正六面體單元將設計域劃分為2688個單元，單元的大小為1mm×1mm×1mm，單元密度的變化范圍為[0.15～1]，在[0.15～0.7]范圍內(nèi)對應變密度點陣材料，其余范圍對應實體材料.拓撲優(yōu)化的初始設計為相對密度為0.4的單一點陣材料填充.拓撲優(yōu)化迭代過程如圖7所示，整個優(yōu)化過程中結構總體積基本保持不變，而結構柔度迭代曲線很快收斂，并最終收斂到189mJ，柔度減小到初始設計的63%.

為了更好地展示最終優(yōu)化結構的細節(jié)，使用合并單元的方法，即將相鄰四個單元的密度取平均值等效為一個單元，重建了最終優(yōu)化結果對應的完整結構，如圖8所示.由圖8可知，機械性能更好的實體材料和高密度點陣材料分布在結構中施加邊界條件處，在受力較小的結構側(cè)上方，分布了低密度點陣材料.在相同工況和網(wǎng)格數(shù)目下，表2展示了本文所提出設計方法與另外兩種點陣材料填充設計方法的優(yōu)化結果對比.其中，設計方案A為本文所提出實體-點陣復合結構拓撲優(yōu)化設計方案，方案B為僅使用變密度點陣材料進行的填充優(yōu)化設計，方案C為均勻點陣材料填充.通過對比可以看出，本文所提出設計方法具有更優(yōu)異的力學性能.

2.2實驗驗證為了進一步驗證本文提出的點陣-實體復合結

構拓撲優(yōu)化方法的有效性，本節(jié)對3種設計方案對應的簡支梁模型進行了實驗測試.通常使用增材制造的加工方式打印形狀復雜的拓撲優(yōu)化模型[20].DLP光固化技術是采用紫外光在液態(tài)光敏樹脂表面進行掃描，從而逐層生成工件的增材制造加工方式[21].它具有加工效率高、成型精度高的優(yōu)點[22].使用DLP型3D打印機（RayshapeShape1）和標準灰色樹脂完成實驗模型的打印，并同時在紫外光照射下固化15min.對圖8中的設計結構進行等比例縮放后，模型的尺為240mm×70mm×20mm，其中每個點陣材料的尺寸為5mm×5mm×5mm.為了避免實驗偶然性，每組加工了兩個模型進行實驗.

在萬能拉伸試驗機上，以2mm/min的速度分別對實驗模型加載，并記錄載荷和位移曲線如圖9所示.由實驗曲線可知，設計方案A優(yōu)化后的模型具有更優(yōu)的剛度，這與數(shù)值算例的結果是相符的.

3結論

本文基于OrderedSIMP插值模型，提出了一種全新的點陣-實體復合結構拓撲優(yōu)化設計方法.總結全文可得到以下結論：

1）運用本文所提出方法可以實現(xiàn)點陣材料的相對密度和其等效力學屬性的高精度的擬合，從而降低多尺度結構拓撲優(yōu)化的計算成本.

2）本文將OrderedSIMP插值模型引入點陣-實體復合結構設計，解決了傳統(tǒng)設計方案中點陣材料和實體材料“分離”設計的缺陷.算法在尋優(yōu)過程自動分配實體材料和點陣材料的比例，避免了人為干預，更具嚴謹性.D4902AB4-8349-4BC7-B388-7D7244E450B4

3）本文所提出的點陣-實體復合結構設計方案可以得到機械性能更優(yōu)的結構，同時發(fā)揮點陣材料填充設計的優(yōu)勢，優(yōu)化結構滿足結構自支撐條件，且具有高能量吸收性和流體通透性.

4）本文所提出的點陣-實體復合結構拓撲優(yōu)化設計方法對不同點陣材料具有廣泛的適用性.

參考文獻

[1] BENDS?E M P，SIGMUND O.Extensions and applications[M]// Topology Optimization.Berlin，Heidelberg：Springer Berlin Hei?d e l b e r g ，2 0 0 4 ：7 1 - 1 5 8 .

[2]杜義賢，杜大翔，李涵釗，等.極限負泊松比的微結構拓撲優(yōu)化[J].機械設計，2018，35（4）：62-66.

[3]盧天健，何德坪，陳常青，等.超輕多孔金屬材料的多功能特性及應用[J].力學進展，2006，36（4）：517-535.

[4] ZHANGCH，LIUJK，YUANZL，etal.Anovellatticestructure topology optimization method with extreme anisotropic lattice prop?erties[J]. Journal of Computational Design and Engineering，2021， 8（5）： 1367-1390.

[5]范華林，楊衛(wèi).輕質(zhì)高強點陣材料及其力學性能研究進展[J].力學進展，2007，37（1）：99-112.

[6] ZHANG L，SONG B，F(xiàn)U J J，et al.Topology-optimized lattice structures with simultaneously high stiffness and light weight fabri?cated by selective laser melting：Design，manufacturing and char?acterization[J].Journal of Manufacturing Processes，2020，56： 1166-1177.

[7] CHENG L，LIU J K，TO A C.Concurrent lattice infill with feature evolution optimization for additive manufactured heat conduction design[J].Structural and Multidisciplinary Optimization，2018， 5 8（ 2 ）：5 1 1 - 5 3 5 .

[8]廖中源，王英俊，王書亭.基于拓撲優(yōu)化的變密度點陣結構體優(yōu)化設計方法[J].機械工程學報，2019，55（8）：65-72.

[9] WU Z J，XIA L，WANG S T，et al.Topology optimization of hier?archical lattice structures with substructuring[J]. Computer Methods in Applied Mechanics and Engineering，2019，345： 602-617.

[10] CHENG L，BAI J X，TO A C.Functionally graded lattice struc?ture topology optimization for the design of additive manufactured components with stress constraints[J].Computer Methods in Ap?plied Mechanics and Engineering，2019，344：334-359.

[11] WANG C，GU X J，ZHU J H，et al.Concurrent design of hierar?chical structures with three-dimensional parameterized lattice mi?crostructures for additive manufacturing[J].Structural and Multi?disciplinary Optimization，2020，61（3）：869-894.

[12] LIU Z，XIA L，XIA Q，et al.Data-driven design approach to hier?archical hybrid structures with multiple lattice configurations[J]. Structural and Multidisciplinary Optimization，2020，61（6）：2227- 2235.

[13] ZUO W J，SAITOU K.Multi-material topology optimization using ordered SIMP interpolation[J]. Structural and Multidisciplinary O p t i m i z a t i o n ，2 0 1 7 ，5 5（ 2 ）：4 7 7 - 4 9 1 .

[14] CHENG G D，CAI Y W，XU L.Novel implementation of homog?enization method to predict effective properties of periodic materi?als[J].Acta Mechanica Sinica，2013，29（4）：550-556.D4902AB4-8349-4BC7-B388-7D7244E450B4

[15] GAO J，LI H，GAO L，et al.Topological shape optimization of 3D micro-structured materials using energy-based homogenization method[J]. Advances in Engineering Software，2018，116： 89-102.

[16] DONG G Y，TANG Y L，ZHAO Y F.A 149 line homogenization code for three-dimensional cellular materials written in Matlab [J].Journal of Engineering Materials and Technology，2019，141 （ 1 ）：0 1 1 0 0 5 .

[17] XIA L，BREITKOPF P.Concurrent topology optimization design of material and structure within FE2 nonlinear multiscale analysis framework[J].Computer Methods in Applied Mechanics and En?g i n e e r i n g ，2 0 1 4 ，2 7 8 ：5 2 4 - 5 4 2 .

[18] SIVAPURAM R，DUNNING P D，KIM H A.Simultaneous mate?rial and structural optimization by multiscale topology optimization [J].Structural and Multidisciplinary Optimization，2016，54（5）： 1267-1281.

[19] SVANBERG K. The method of moving asymptotes—a new method for structural optimization[J]. International Journal for Numerical Methods in Engineering，1987，24（2）：359-373.

[20] LIU J K，GAYNOR A T，CHEN S K，et al.Current and future trends in topology optimization for additive manufacturing[J]. Structural and Multidisciplinary Optimization，2018，57（6）：2457- 2483.

[21]劉雨，陳張偉.陶瓷光固化3D打印技術研究進展[J].材料工程，2020，48（9）：1-12.

[22]吳甲民，楊源祺，王操，等.陶瓷光固化技術及其應用[J].機械工程學報，2020，56（19）：221-238.D4902AB4-8349-4BC7-B388-7D7244E450B4