王 昱, 呂恩利, 王亞娟, 王飛仁, 林小娟

(華南農(nóng)業(yè)大學(xué) 工程學(xué)院,廣州 510642)

負(fù)泊松比負(fù)熱膨脹超材料微結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)

王 昱, 呂恩利, 王亞娟, 王飛仁, 林小娟

(華南農(nóng)業(yè)大學(xué) 工程學(xué)院,廣州 510642)

提出一種拓?fù)鋬?yōu)化設(shè)計(jì)方法,用于設(shè)計(jì)同時(shí)具有負(fù)泊松比和負(fù)熱膨脹系數(shù)的超材料.超材料是人工設(shè)計(jì)的復(fù)合材料,由微米或納米結(jié)構(gòu)的周期性陣列組成,并具有天然材料難以具備的超常物理性質(zhì).基于拓?fù)鋬?yōu)化理論,建立了一個(gè)表示多材料微觀結(jié)構(gòu)的多相邊界的水平集模型,并采用數(shù)值均勻化方法來(lái)計(jì)算微觀結(jié)構(gòu)的等效性能,然后通過(guò)參數(shù)化水平集優(yōu)化方法實(shí)現(xiàn)邊界形狀的演變和拓?fù)浣Y(jié)構(gòu)的設(shè)計(jì).通過(guò)兩個(gè)數(shù)值算例,給出超材料微結(jié)構(gòu)材料分布的優(yōu)化結(jié)果,證明了設(shè)計(jì)方法的可行性,為實(shí)現(xiàn)具備特殊性能的人工材料設(shè)計(jì)提供了一套系統(tǒng)的設(shè)計(jì)方法.

超材料; 拓?fù)鋬?yōu)化; 水平集方法; 負(fù)泊松比; 負(fù)熱膨脹

在工程應(yīng)用中,溫度變化會(huì)引起尺寸變化,所以必須考慮材料和結(jié)構(gòu)的膨脹與收縮.由于大多數(shù)材料的應(yīng)用環(huán)境將經(jīng)歷一系列的工作溫度,無(wú)論是技術(shù)還是基礎(chǔ)研究領(lǐng)域都會(huì)關(guān)注由溫度變化引起的結(jié)構(gòu)的熱膨脹[10].熱膨脹系數(shù)是一個(gè)重要的物理量,用來(lái)測(cè)量溫度對(duì)材料尺寸的改變程度.負(fù)熱膨脹材料(NTE)在一定的溫度升高范圍內(nèi)呈現(xiàn)收縮的熱彈性能.相關(guān)研究證明,由三相材料組成的復(fù)合材料的性能可以超越其中單相材料的性能[11].因此,極端的熱彈性能的材料通常是通過(guò)組合三相傳統(tǒng)的材料,如金屬或塑料等,得到一個(gè)新的具有周期結(jié)構(gòu)的復(fù)合材料.

泊松比[1]被定義為在單軸拉伸條件下,橫向收縮應(yīng)變與縱向拉伸應(yīng)變的比值,這是衡量彈性材料性能的基本指標(biāo).彈性理論允許泊松比為負(fù),但大多數(shù)傳統(tǒng)材料泊松比為正.與正泊松比的常規(guī)材料相比,具有負(fù)泊松比的材料(NPR)表現(xiàn)出非常規(guī)的屬性,即橫向拉伸時(shí)縱向膨脹,橫向收縮時(shí)縱向壓縮.由于廣泛的應(yīng)用潛力,自從早期研究工作以來(lái)[12],負(fù)泊松比材料已經(jīng)吸引了越來(lái)越多的關(guān)注.到目前為止,已有一些方法開(kāi)發(fā)設(shè)計(jì)拉脹材料[13-15],但鮮有系統(tǒng)的設(shè)計(jì)方法用于設(shè)計(jì)同時(shí)表現(xiàn)NPR和NTE的超材料.

拓?fù)鋬?yōu)化已被確定為是在結(jié)構(gòu)優(yōu)化領(lǐng)域最有前途的技術(shù)之一[16],其中,水平集方法因其具有光滑的設(shè)計(jì)邊界,并集成了形狀優(yōu)化和拓?fù)鋬?yōu)化,已經(jīng)成為一種備受關(guān)注的形狀結(jié)構(gòu)拓?fù)鋬?yōu)化方法.特別是參數(shù)化水平集方法(PLSM)[17],已被證明是一個(gè)最強(qiáng)大的主流拓?fù)湫螤顑?yōu)化的結(jié)構(gòu)和機(jī)制.它能很好地繼承傳統(tǒng)水平集方法的優(yōu)點(diǎn),并使許多有效的基于梯度的優(yōu)化算法[18]直接應(yīng)用到水平集方法中[19].

因此,本文基于多相水平集模型來(lái)設(shè)計(jì)同時(shí)具有NPR和NTE的超材料拓?fù)湫螤?多相水平集模型是一個(gè)顯示多材料的邊界的水平集函數(shù),該函數(shù)能避免設(shè)計(jì)域出現(xiàn)材料重疊與空白.文中還將采用基于PLSM的超材料微結(jié)構(gòu)設(shè)計(jì)方法.在這種方法中,應(yīng)用了數(shù)值均勻化方法來(lái)計(jì)算復(fù)合材料的宏觀等效性能,并使用PLSM來(lái)實(shí)現(xiàn)形狀和多材料微結(jié)構(gòu)拓?fù)溲葑?最后,通過(guò)兩個(gè)典型的數(shù)值算例證明所提出方法的有效性.

1 拓?fù)鋬?yōu)化方法

1.1 參數(shù)化水平集方法

在水平集方法(LSM)方法中,設(shè)計(jì)邊界通過(guò)高維水平集函數(shù)的零水平集形成的隱式邊界表達(dá)[20-21].一個(gè)水平集函數(shù)Φ(x)要求在參考域D內(nèi)具有Lipschitz連續(xù)性.對(duì)于本文涉及的多相材料設(shè)計(jì)問(wèn)題,采用m個(gè)水平集函數(shù)相結(jié)合來(lái)描述m+1相材料,包括無(wú)材料相(即孔洞)和m相固體材料.固體材料、邊界及孔洞的水平集表達(dá)式分別為

(1)

式中:I=1,…,m為某一水平集函數(shù);m是水平集函數(shù)的總數(shù);D為設(shè)計(jì)參考域,包括所有允許的微觀結(jié)構(gòu)的形狀;ΩI是函數(shù)值為正的水平集函數(shù);?ΩI表示第I個(gè)水平集函數(shù)的設(shè)計(jì)邊界.

我最喜歡的蔬菜是西紅柿。西紅柿圓圓的，像個(gè)小燈籠似的。它的顏色各異：有的是粉色的，上面還隱隱約約的白點(diǎn)，媽媽說(shuō)這樣的西紅柿是沙瓤的，吃起來(lái)更美味；有的是紅紅的，像是小姑娘羞紅了的臉蛋，越看越可愛(ài)，我都有點(diǎn)不舍得吃了；還有黃綠相間的西紅柿，據(jù)說(shuō)是科學(xué)家叔叔們發(fā)明的新品種，聞起來(lái)有股淡淡的清香，吃起來(lái)酸甜可口，味道純正，科學(xué)家叔叔們真了不起！

在PLSM方法中,采用緊支徑向基函數(shù)(CS-RBF)插值得到高維水平集函數(shù),有利于平衡插值精度和效率.原來(lái)的水平集函數(shù)可以用CS-RBF在預(yù)先指定節(jié)點(diǎn)的插值表達(dá),即

形狀函數(shù)向量為

(3)

擴(kuò)展系數(shù)向量為

(4)

式中:N是設(shè)計(jì)域內(nèi)CS-RBF節(jié)點(diǎn)總數(shù);t為虛擬時(shí)間.

通過(guò)CS-RBFs插值可實(shí)現(xiàn)時(shí)間和空間的分離:形狀函數(shù)為空間函數(shù),膨脹系數(shù)僅與時(shí)間相關(guān).值得注意的是,用于插值每個(gè)隱式水平集函數(shù)的CS-RBFs是相同的.因此,解耦漢密爾頓-雅各比偏微分方程(H-J PDE)的時(shí)間相和空間相可得

(5)

(6)

由式(6)可知,正速度矢量所涉及的項(xiàng)都是通過(guò)設(shè)計(jì)域內(nèi)的節(jié)點(diǎn)求得,因此正速度擴(kuò)展到了整個(gè)設(shè)計(jì)領(lǐng)域.這樣,基于偏微分方程的水平集方程已轉(zhuǎn)化為代數(shù)方程.優(yōu)化唯一的未知數(shù)是插值的膨脹系數(shù),在優(yōu)化設(shè)計(jì)中定義為設(shè)計(jì)變量,并使用更高效的優(yōu)化算法進(jìn)行迭代.

1.2 三相復(fù)合材料的材料模型

通過(guò)材料模型[5],在設(shè)計(jì)域內(nèi)的任何計(jì)算點(diǎn)x的彈性常數(shù)C和熱應(yīng)變系數(shù)α可以由兩個(gè)水平集函數(shù)Φ1和Φ2計(jì)算得出.

(7)

(8)

1.3 超材料微結(jié)構(gòu)的形狀拓?fù)鋬?yōu)化

將在給定的設(shè)計(jì)域,單胞Y,進(jìn)行形狀拓?fù)鋬?yōu)化.復(fù)合材料的等效彈性張量和熱彈系數(shù)采用均勻化方法計(jì)算,得

(9)

(10)

(11)

(12)

(13)

由于水平集函數(shù)的參數(shù)化,許多更高效的基于梯度的優(yōu)化算法可用于求解優(yōu)化問(wèn)題,例如:移動(dòng)漸近線法(MMA)[18].采用MMA需要提供目標(biāo)函數(shù)和約束條件的關(guān)于設(shè)計(jì)變量的一階導(dǎo)數(shù).目標(biāo)函數(shù)關(guān)于設(shè)計(jì)變量的一階導(dǎo)數(shù)表示為

其中,

式中:χij為位移場(chǎng);αpq為局部熱應(yīng)變場(chǎng).

體積約束相對(duì)于設(shè)計(jì)變量的導(dǎo)數(shù)為

(18)

式中,δ為狄拉克函數(shù),是Heaviside函數(shù)的一階導(dǎo)數(shù).

2 數(shù)值算例

用來(lái)計(jì)算位移場(chǎng),設(shè)計(jì)域采用4個(gè)節(jié)點(diǎn)的四邊形單元離散,且水平集的節(jié)點(diǎn)與四邊形單元的節(jié)點(diǎn)一致.優(yōu)化過(guò)程中,水平集函數(shù)采用MMA方法更新設(shè)計(jì)變量.收斂準(zhǔn)則為連續(xù)兩次迭代的目標(biāo)函數(shù)值差小于10-4,或最大迭代次數(shù)達(dá)到200次(基于數(shù)值試驗(yàn)).設(shè)定兩相虛擬固體材料的彈性模量:E1=E2=1;泊松比:v1=v2=0.3;熱膨脹系數(shù):α1=1,α2=10.

2.1 算例1

算例1設(shè)計(jì)了一種正交各向異性材料,以實(shí)現(xiàn)負(fù)的泊松比(-0.5),以及其在垂直方向上的負(fù)熱應(yīng)變系數(shù)(-0.3).兩相固體材料體積分?jǐn)?shù)分別為0.3和0.1.拓?fù)湫螤畹难莼褪諗窟^(guò)程如圖1所示.紅色區(qū)域表示固體材料相1,藍(lán)色區(qū)域表示固體材料相2,而白色區(qū)域是孔洞.兩個(gè)水平集的表面和對(duì)應(yīng)于所得到的拓?fù)浣Y(jié)構(gòu)的兩個(gè)固相的零水平集輪廓.式(19)為優(yōu)化所得彈性張量與熱膨脹系數(shù).

(19)

由式(19)可知,優(yōu)化后材料的泊松比為-0.446 3,垂直方向熱應(yīng)變系數(shù)為-0.261 6.

從算例1的迭代過(guò)程可以注意到,形狀拓?fù)鋬?yōu)化設(shè)計(jì)方法的特點(diǎn)是擁有光滑的材料邊界和明顯的界線[20-21].此外,在優(yōu)化過(guò)程中,所提出的方法不僅可以合并初始設(shè)計(jì)中的孔洞,也可以創(chuàng)建新的孔洞,以實(shí)現(xiàn)最佳的拓?fù)浣Y(jié)構(gòu)的基本單元.基于所得到的最優(yōu)設(shè)計(jì)(如圖2所示),注意到所得的微結(jié)構(gòu)的拓?fù)浣Y(jié)構(gòu)與文獻(xiàn)[9]所獲得的結(jié)果相似.可知等效的彈性和熱彈性能是依賴于施加的負(fù)載后單胞的內(nèi)部結(jié)構(gòu)和變形的方式.因此,設(shè)計(jì)一種具有極限性能的復(fù)合材料,微結(jié)構(gòu)的設(shè)計(jì)將會(huì)出現(xiàn)類似旋轉(zhuǎn)剛性機(jī)構(gòu)[1-2]的結(jié)構(gòu).

圖1 目標(biāo)函數(shù)的收斂過(guò)程(算例1)

圖2 固體材料相輪廓圖(算例1)

2.2 算例2

算例2設(shè)計(jì)了一個(gè)正交各向異性材料,以實(shí)現(xiàn)負(fù)的泊松比(-0.65),以及在水平方向上的負(fù)熱應(yīng)變系數(shù),且熱應(yīng)變系數(shù)達(dá)到極小值.兩相固體材料體積分?jǐn)?shù)分別為0.25和0.10.拓?fù)湫螤畹难莼褪諗窟^(guò)程如圖3所示,式(20)為優(yōu)化所得彈性張量.

(20)

由式(20)可知,優(yōu)化后材料的泊松比為-0.644 2,水平方向熱應(yīng)變系數(shù)為-3.356.

式(20)結(jié)果表明,極端的熱膨脹的形成會(huì)導(dǎo)致一個(gè)非常低的結(jié)構(gòu)剛度,即實(shí)現(xiàn)一個(gè)方向的極端熱應(yīng)變將導(dǎo)致在另一個(gè)方向上的一個(gè)極端的剛度.因此,在一個(gè)特定的方向設(shè)計(jì)極端熱膨脹材料,應(yīng)在優(yōu)化約束中定義另一方向上的最小剛度值,使設(shè)計(jì)具有物理意義.在算例2中,設(shè)定垂直方向的體積模量的下限為0.09.

材料的等效彈性和熱彈性能取決于基本單元的內(nèi)部結(jié)構(gòu)和變形方式,為使微結(jié)構(gòu)在加載時(shí)發(fā)生旋轉(zhuǎn),凹角結(jié)構(gòu)的生成是使結(jié)構(gòu)發(fā)生較大旋轉(zhuǎn)效果的關(guān)鍵.在微結(jié)構(gòu)的特定布局下,可以通過(guò)類似旋轉(zhuǎn)剛性機(jī)構(gòu)的微結(jié)構(gòu)基本單元之間的角度變化,使多孔復(fù)合材料整體變形,達(dá)到負(fù)的極端力學(xué)性能指數(shù),即

圖3 目標(biāo)函數(shù)的收斂過(guò)程(算例2)

負(fù)泊松比和負(fù)熱膨脹系數(shù).在上述兩個(gè)數(shù)值例子中,可發(fā)現(xiàn)最佳的泊松比和熱應(yīng)變系數(shù)不能準(zhǔn)確地接近規(guī)定的值.由于拓?fù)湫螤顑?yōu)化最初是連續(xù)體結(jié)構(gòu)的一種設(shè)計(jì)方法,它實(shí)際上是利用結(jié)構(gòu)的變形來(lái)近似模擬剛體的旋轉(zhuǎn)效應(yīng),因此難以使微結(jié)構(gòu)實(shí)現(xiàn)如剛性連桿機(jī)構(gòu)般的旋轉(zhuǎn)變形.在優(yōu)化過(guò)程中會(huì)發(fā)生在折返結(jié)構(gòu)部分的類似鉸鏈結(jié)構(gòu)斷開(kāi)的現(xiàn)象,為保證優(yōu)化所得的結(jié)果為連續(xù)體結(jié)構(gòu),所得結(jié)果無(wú)法準(zhǔn)確地達(dá)到設(shè)定的泊松比和熱應(yīng)變系數(shù).盡管如此,該方法仍然提供了一個(gè)系統(tǒng)的復(fù)合材料微結(jié)構(gòu)設(shè)計(jì)方法,使材料同時(shí)達(dá)到預(yù)期的負(fù)泊松比和負(fù)熱膨脹系數(shù).

圖4 固體材料相輪廓圖(算例2)

3 結(jié) 論

利用數(shù)值均勻化方法和參數(shù)化的水平集方法,提出了同時(shí)具有NPR和NTE性質(zhì)復(fù)合材料的系統(tǒng)設(shè)計(jì)方法.數(shù)值算例的結(jié)果驗(yàn)證了所提出的拓?fù)湫螤钤O(shè)計(jì)方法的有效性.材料微結(jié)構(gòu)設(shè)計(jì)可以從二維推廣到三維,但計(jì)算成本將顯著增加.同時(shí),該方法具有拓展至其他類型超材料設(shè)計(jì)的潛力.

[1] EVANS K E,ALDERSON A.Auxetic materials:functional materials and structures from lateral thinking![J].Advanced Materials,2000,12(9):617-628.

[2] LAKES R.Cellular solid structures with unbounded thermal expansion[J].Journal of Materials Science Letters,1996,15(6):475-477.

[3] GAO Y,ZHOU N,YANG F,et al.P-phase precipitation and its effect on martensitic transformation in (Ni,Pt)Ti shape memory alloys[J].Acta Materialia,2012,60(4):1514-1527.[4] GAO Y,ZHOU N,WANG D,et al.Pattern formation during cubic to orthorhombic martensitic transformations in shape memory alloys[J].Acta Materialia,2014,68:93-105.

[5] SIGMUND O.A new class of extremal composites[J].Journal of the Mechanics and Physics of Solids,2000,48(2):397-428.

[6] GUEST J K,PRéVOST J H.Design of maximum permeability material structures[J].Computer Methods in Applied Mechanics and Engineering,2007,196(4/5/6):1006-1017.

[7] 張衛(wèi)紅,汪雷,孫士平.基于導(dǎo)熱性能的復(fù)合材料微結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)[J].航空學(xué)報(bào),2006,27(6):1229-1233.

[8] 孫士平,張衛(wèi)紅.多相材料微結(jié)構(gòu)多目標(biāo)拓?fù)鋬?yōu)化設(shè)計(jì)[J].力學(xué)學(xué)報(bào),2006,38(5):633-638.

[9] 孫良奎,于哲峰,黃潔.基于超材料的平板二維定向傳熱結(jié)構(gòu)設(shè)計(jì)[J].物理學(xué)報(bào),2015,64(22):0224401.

[10] SIGMUND O,TORQUATO S.Composites with extremal thermal expansion coefficients[J].Applied Physics Letters,1996,69(21):3203-3205.

[11] LAKES R.Materials with structural hierarchy[J].Nature,1993,361(6412):511-515.

[12] LAKES R.Foam structures with a negative Poisson's ratio[J].Science,1987,235(4792):1038-1040.

[13] YANG W,LI Z M,SHI W,et al.Review on auxetic materials[J].Journal of Materials Science,2004,39(10):3269-3279.

[14] 曹先凡,劉書(shū)田.基于拓?fù)涿枋龊瘮?shù)的特定性能材料設(shè)計(jì)方法[J].固體力學(xué)學(xué)報(bào),2006,27(3):217-222.

[15] 宋迎東,孫杰,孫志剛,等.基于遺傳算法的復(fù)合材料細(xì)觀結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)[J].固體力學(xué)學(xué)報(bào),2009,30(4):416-423.

[16] BENDSOE M P,SIGMUND O.Topology optimization:theory,methods and applications[M].Berlin Heidelberg:Springer-Verlag,2004.

[17] LUO Z,TONG L Y,WANG M Y,et al.Shape and topology optimization of compliant mechanisms using a parameterization level set method[J].Journal of Computational Physics,2007,227(1):680-705.

[18] 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.

[19] BELYTSCHKO T,XIAO S P,PARIMI C.Topology optimization with implicit functions and regularization[J].International Journal for Numerical Methods in Engineering,2003,57(8):1177-1196.

[20] ALLAIRE G,JOUVE F,TOADER A M.Structural optimization using sensitivity analysis and a level-set method[J].Journal of Computational Physics,2004,194(1):363-393.

[21] WANG M Y,WANG X M,GUO D M.A level set method for structural topology optimization[J].Computer Methods in Applied Mechanics and Engineering,2003,192(1/2):227-246.

(編輯:丁紅藝)

Topology Optimization Design of Microstructured Metamaterials with Negative Poisson’s Ratio and Negative Thermal Expansion

WANG Yu, Lü Enli, WANG Yajuan, WANG Feiren, LIN Xiaojuan

(CollegeofEngineering,SouthChinaAgriculturalUniversity,Guangzhou510642,China)

A level set-based method was introduced to design the metamaterials which can exhibit negative Poisson’s ratio and negative thermal expansion simultaneously.Negative Poisson’s ratio (also called auxeticity) and negative thermal expansion are phenomena that may rarely be found in the nature materials.It is now known that these unusual properties are mainly dependent on the particular features in the material’s micro- or nano- structural deformation when subjected to a thermal or mechanical load.A numerical homogenization method was applied to evaluate the effective properties of the microstructure,and the multi-phase level set method was used to implement the shape and topology evolutions of the microstructure.Two numerical examples were given to demonstrate the feasibility of the proposed topology optimization method,which provides a systematic design method for the design of artificial materials with special properties.

metamaterial;topologicalshapeoptimization;letsetmethod;negativePoisson’sratio;negativethermalexpansion

1007-6735(2016)06-0551-06

10.13255/j.cnki.jusst.2016.06.008

2016-09-08

廣東省農(nóng)產(chǎn)品冷鏈物流裝備工程技術(shù)研究中心建設(shè)項(xiàng)目(2016B090920092);華南農(nóng)業(yè)大學(xué)高水平青年科技人才培育基金項(xiàng)目(216262)

王 昱(1987-),女,講師.研究方向:結(jié)構(gòu)拓?fù)鋬?yōu)化、汽車輕量化設(shè)計(jì)等.E-mail:yu-wang@scau.edu.cn

呂恩利(1979-),男,教授.研究方向:果蔬冷鏈物流技術(shù)與裝備、交通運(yùn)輸、汽車維修檢測(cè)等.E-mail:enlilv@scau.edu.cn

TH 122

A