康小方，柳景超，夏光輝，3，4，許慶虎，3，4

（1.安徽建筑大學(xué) 土木工程學(xué)院，安徽 合肥 230601；2.智能地下探測技術(shù)重點實驗室，安徽 合肥 230601；3.安徽省BIM 工程中心，安徽 合肥 230601；4.安徽省裝配式建筑研究院，安徽 合肥 230601）

泊松比是用作材料的泊松比效應(yīng)的度量系數(shù)，泊松比也被稱作材料的橫向變形系數(shù)，材料在力作用與力作用垂直方向上的橫向變形特性。泊松比是指材料在受力加載方向下，其橫向應(yīng)變與縱向應(yīng)變比值的相反數(shù)[1-3]。生活中所使用的材料泊松比取值一般不小于0，常見的金屬泊松比一般在0.3左右。液體是一種極為致密的材料且不可被壓縮，在壓力的作用下體積幾乎沒有變化，其泊松比接近0.5。

約160 年前，經(jīng)典彈性理論提出了材料出現(xiàn)負(fù)泊松比的可能性[4]。常見的內(nèi)凹結(jié)構(gòu)包括內(nèi)凹六邊形、雙箭頭模型、星型和手性等結(jié)構(gòu)。1982年，Gibson 等[5]首次提出了內(nèi)凹六邊形蜂窩模型；Yang S 等[6]提出了一種矩形穿孔的負(fù)泊松比結(jié)構(gòu)。采用拓?fù)鋬?yōu)化方法設(shè)計負(fù)泊松比超材料為超材料結(jié)構(gòu)優(yōu)化設(shè)計提供了一種新的途徑，該方法已廣泛應(yīng)用于負(fù)泊松比超材料結(jié)構(gòu)優(yōu)化[7-8]。

負(fù)泊松材料主要類型有[9]：①天然負(fù)泊松比材料；②胞狀負(fù)泊松比材料；③金屬負(fù)泊松比材料；④由多種普通基材所組成的多重負(fù)泊松比材料；⑤由負(fù)泊松比材料和非負(fù)泊松比材料復(fù)合形成的負(fù)泊松比復(fù)合材料。

與傳統(tǒng)泊松比材料相比，負(fù)泊松比超材料幾乎具有兩倍的抗斷裂性[10-12]，具有負(fù)泊松比效應(yīng)的復(fù)合材料抗斷裂性也大約是傳統(tǒng)復(fù)合材料的兩倍。橋梁的伸縮縫材料要求其性能滿足垂直和平行于橋梁軸線的兩個方向都可以自由變形。然而，負(fù)泊松比超材料具有受到壓力時橫向收縮、受到拉力時橫向膨脹的特性。負(fù)泊松比超材料安裝在梁體間隙之中擁有足夠的變形能力，同時負(fù)泊松比效應(yīng)還能夠?qū)Σ牧系臋M向和豎向承載能力有一定的增強。隨著負(fù)泊松比超材料結(jié)構(gòu)的快速發(fā)展，已經(jīng)涉及到人工假體[13]、智能傳感器[14-16]、紡織類材料、分子過濾器[17]、減震器[18]、隔聲器[19]和防護(hù)墊[20]等領(lǐng)域。1989 年，Bendsoe[21]采用微結(jié)構(gòu)[22]的拓?fù)鋬?yōu)化法對結(jié)構(gòu)進(jìn)行拓?fù)鋬?yōu)化。

本文將微結(jié)構(gòu)的設(shè)計轉(zhuǎn)換成了胞元拓?fù)鋬?yōu)化設(shè)計，采用胞元拓?fù)鋬?yōu)化法對內(nèi)凹六邊形微結(jié)構(gòu)進(jìn)行拓?fù)鋬?yōu)化，將優(yōu)化得到的胞元微結(jié)構(gòu)按一定周期進(jìn)行排列從而獲得負(fù)泊松比超材料結(jié)構(gòu)。

1 內(nèi)凹六邊形負(fù)泊松材料計算理論

1.1 負(fù)泊松比材料作用機理

泊松比是衡量材料在承受縱向力時的橫向尺寸變化力學(xué)參數(shù)。負(fù)泊松比是在考慮正負(fù)應(yīng)變的前提下，橫向應(yīng)變和縱向應(yīng)變的比值相反數(shù)，其數(shù)學(xué)表達(dá)式為，

式中，εx，εz分別為橫向應(yīng)變和縱向應(yīng)變；Δx，Δz分別為結(jié)構(gòu)受荷后的橫向變形和縱向變形。一般情況下，泊松比大小和材料內(nèi)部結(jié)構(gòu)相關(guān)，經(jīng)典彈性力學(xué)理論證明了具有各向同性的材料其泊松比值范圍為-1.0 ～0.5[23]。從泊松比取值范圍上來看，材料具有負(fù)泊松比效應(yīng)理論上是可行的。

1.2 內(nèi)凹六邊形拓?fù)鋬?yōu)化

本文從材料（微觀尺度）到結(jié)構(gòu)（宏觀尺度）進(jìn)行研究。所采用的方法是在保持宏觀材料分布不變的前提下，對微觀尺度材料下單胞結(jié)構(gòu)進(jìn)行拓?fù)鋬?yōu)化設(shè)計。將優(yōu)化得到的胞元微結(jié)構(gòu)按一定周期進(jìn)行排列從而獲得負(fù)泊松比超材料結(jié)構(gòu)。如圖1 所示：

圖1 微觀-宏觀-一體化

內(nèi)凹六邊形源自常規(guī)六邊形結(jié)構(gòu)，但是其兩側(cè)呈內(nèi)凹結(jié)構(gòu)，其中：a 為底部邊長，b 為斜邊長，c 為肋部邊長，a，c 為常量；θ0為初始角（＜ 90°）。如圖2 所示：

圖2 內(nèi)凹六邊形結(jié)構(gòu)胞元示意圖

根據(jù)負(fù)泊松比定義和內(nèi)凹六邊形結(jié)構(gòu)圖，負(fù)泊松比可表示為，

根據(jù)幾何分析，橫向應(yīng)變εa可表示為，

式中，θ0，θ和分別為胞元結(jié)構(gòu)的初始角度和變形后角度。

由圖2 可知胞元結(jié)構(gòu)尺寸b，a，c 和初始角度θ0之間的關(guān)系可表述為，

于是，胞元結(jié)構(gòu)的泊松比v為[24]，

由式（4）代入式（5）可得下式，

采用Matlab 軟件對式（3）和式（6）進(jìn)行計算，可得胞元結(jié)構(gòu)的橫向應(yīng)變εa（圖3）和泊松比v（圖4）。

圖3 結(jié)構(gòu)胞元橫向應(yīng)變εa 等高線

圖4 結(jié)構(gòu)胞元泊松比v 等高線

具有負(fù)泊松比效應(yīng)的內(nèi)凹六邊形拓?fù)鋬?yōu)化步驟可表述為：

（1）通過圖1 所示的微觀-宏觀一體化優(yōu)化思路，將宏觀結(jié)構(gòu)看作由一系列內(nèi)凹六邊形結(jié)構(gòu)胞元（圖2）所組成。

（2）對內(nèi)凹六邊形結(jié)構(gòu)胞元進(jìn)行結(jié)構(gòu)優(yōu)化，通過改變圖2 所示的初始角θ0，來實現(xiàn)結(jié)構(gòu)胞元的負(fù)泊松比效應(yīng)。

（3）根據(jù)步驟（2）獲得的負(fù)泊松比效應(yīng)結(jié)構(gòu)胞元，按一定周期進(jìn)行排列組成負(fù)泊松比超材料結(jié)構(gòu)。

（4）根據(jù)步驟（3）獲得的負(fù)泊松比超材料結(jié)構(gòu)，采用ANSYS 有限元軟件建立相應(yīng)的結(jié)構(gòu)模型。

2 內(nèi)凹六邊形胞元結(jié)構(gòu)力學(xué)行為

為驗證基于胞元拓?fù)鋬?yōu)化法的內(nèi)凹六邊形負(fù)泊松比超材料結(jié)構(gòu)受力性能，本節(jié)采用大型有限元軟件ANSYS 進(jìn)行建模，對優(yōu)化后的結(jié)構(gòu)胞元進(jìn)行受力分析和泊松比驗證。本文建立代表性內(nèi)凹六邊形結(jié)構(gòu)胞元，假定結(jié)構(gòu)胞元尺寸a 和c 為常量，初始角度θ0為變量，結(jié)構(gòu)胞元尺寸a=3 mm，c=1 mm，初始角度θ0=30°、45°和60°，材料彈性模量為1.2×105MPa，材料泊松比為0.29。

為了保證結(jié)構(gòu)胞元可以兩端受拉，利用作用力與反作用力的關(guān)系將結(jié)構(gòu)胞元邊界條件定為下端固定約束，胞元上端施加5 MPa 大小的拉應(yīng)力。

2.1 初始角度θ0=30°內(nèi)凹六邊形

結(jié)構(gòu)胞元的縱向尺寸z=1.155 mm，x=3 mm。結(jié)構(gòu)胞元經(jīng)加荷載后，其縱向變形為Δz=-0.445 51 mm，結(jié)構(gòu)胞元的橫向變形Δx=-0.195 148 mm。將有限元仿真結(jié)果代入公式（1）可以得到材料的泊松比，泊松比=-0.167 9。

2.2 初始角度θ0=45°內(nèi)凹六邊形

結(jié)構(gòu)胞元初始角度θ0=45°，其縱向尺寸z=2 mm，x=3 mm。經(jīng)加荷載后，結(jié)構(gòu)胞元的縱向變形為Δz=-0.634 55 mm，結(jié)構(gòu)胞元的橫向變形Δx=-0.511 254 mm。將有限元仿真結(jié)果代入公式（1）可以得到材料的泊松比，泊松比v=-0.537 1。

2.3 初始角度θ0=60°內(nèi)凹六邊形

結(jié)構(gòu)胞元初始角度θ0=60°，結(jié)構(gòu)胞元的縱向尺寸z=3.464 mm，x=3 mm。在加載荷載后，其縱向變形為Δz=-0.717 286 mm，結(jié)構(gòu)胞元的橫向變形Δx=-0.554 977 mm。將有限元仿真結(jié)果代入公式（1）可以得到材料的泊松比，泊松比v=-1.492 36。

由圖5 ～圖10 可知，隨著結(jié)構(gòu)胞元初始角度的增大，內(nèi)凹六邊形結(jié)構(gòu)胞元負(fù)泊松比值也越大。

圖5 θ0=30°Z 方向變形圖

圖6 θ0=30°X 方向變形圖

圖7 θ0=45°Z 方向變形圖

圖8 θ0=45°X 方向變形圖

圖9 θ0=60°Z 方向變形圖

圖10 θ0=60°X 方向變形圖

為驗證有限元模型的準(zhǔn)確性，將受荷情況下結(jié)構(gòu)胞元的負(fù)泊松比效應(yīng)和理論公式進(jìn)行對比。根據(jù)結(jié)構(gòu)胞元尺寸a=3，c=1，初始角度θ0=30°、45°和60°，分別代入式（6）可得v30°=-0.168 6、v45°=-0.505 0、v60°=-1.517 5，將有限元仿真結(jié)果和理論計算結(jié)果進(jìn)行匯總，可得泊松比值如表1 所示：

表1 泊松比值

將有限元計算的結(jié)果與理論計算的結(jié)果進(jìn)行比對，發(fā)現(xiàn)誤差均在10%以內(nèi)，從側(cè)面反映出有限元仿真的準(zhǔn)確性及負(fù)泊松比理論計算的合理性。

3 內(nèi)凹六邊形超材料結(jié)構(gòu)性能仿真

3.1 內(nèi)凹六邊形超材料宏觀結(jié)構(gòu)陣列

為分析材料拓?fù)鋬?yōu)化后的超材料結(jié)構(gòu)宏觀力學(xué)性能，本文將經(jīng)過拓?fù)鋬?yōu)化后的第3 節(jié)內(nèi)凹六邊形按一定周期進(jìn)行排列得到相應(yīng)具有宏觀負(fù)泊松比效應(yīng)的超材料結(jié)構(gòu)。如圖11 所示：

圖11 內(nèi)凹六邊形宏觀結(jié)構(gòu)陣列示意圖

3.2 內(nèi)凹六邊形超材料結(jié)構(gòu)有限元仿真

根據(jù)圖9 所示，結(jié)構(gòu)胞元按照一定的規(guī)則進(jìn)行排列，可獲得相應(yīng)的超材料結(jié)構(gòu)。采用有限元軟件ANSYS 進(jìn)行建模仿真，超材料結(jié)構(gòu)邊界條件為下端固定約束，胞元上端施加5 MPa 大小的拉應(yīng)力。

3.2.1 初始角度θ0=30°超材料結(jié)構(gòu)

由一系列初始角度θ0=30°的結(jié)構(gòu)胞元組成的超材料結(jié)構(gòu)，其縱向尺寸z=3.462 mm，橫向尺寸x=11 mm。經(jīng)加荷載后，超材料結(jié)構(gòu)的縱向變形為Δz=-0.380 01 mm，超材料結(jié)構(gòu)的橫向變形Δx=-0.136 8 mm。經(jīng)公式（1）計算可得超材料結(jié)構(gòu)的泊松比v=-0.113。

3.2.2 初始角度θ0=45°超材料結(jié)構(gòu)

根據(jù)初始角度θ0=45°的結(jié)構(gòu)胞元，經(jīng)過一定的排列規(guī)律組成宏觀超材料結(jié)構(gòu)，其縱向尺寸z=6 mm，橫向尺寸x=11 mm。經(jīng)加荷載后，超材料結(jié)構(gòu)的縱向變形為Δz=-1.13555 mm，超材料結(jié)構(gòu)的橫向變形Δx=-0.598 61 mm。經(jīng)公式（1）計算可得超材料結(jié)構(gòu)的泊松比v=-0.288。

3.2.3 初始角度θ0=60°超材料結(jié)構(gòu)

采用初始角度θ0=60°的結(jié)構(gòu)胞元，經(jīng)組合成超材料結(jié)構(gòu)，其縱向尺寸z=6.928 mm，橫向尺寸x=11 mm。經(jīng)加荷載后，超材料結(jié)構(gòu)的縱向變形為Δz=-1.538 47 mm，超材料結(jié)構(gòu)的橫向變形Δx=-1.299 89 mm。經(jīng)公式（1）計算可得超材料結(jié)構(gòu)的泊松比v=-0.532。

由圖12～圖17 可得，內(nèi)凹六邊形超材料結(jié)構(gòu)的負(fù)泊松比效應(yīng)跟結(jié)構(gòu)胞元特性有關(guān)。隨著結(jié)構(gòu)胞元初始角度增大，超材料結(jié)構(gòu)的負(fù)泊松比值也越大，其變化規(guī)律與結(jié)構(gòu)胞元負(fù)泊松比效應(yīng)類似。

圖12 =30°Z 方向變形圖

圖13 =30°X 方向變形圖

圖14 =45°Z 方向變形圖

圖15 =45°X 方向變形圖

圖16 =60°Z 方向變形圖

圖17 =60°X 方向變形圖

4 結(jié)論

本文采用胞元拓?fù)鋬?yōu)化的方法，以超材料負(fù)泊松比效應(yīng)增強為目的，推導(dǎo)得出了內(nèi)凹角越大其負(fù)泊松比效應(yīng)越強的理論，并應(yīng)用有限元軟件ANSYS 對其進(jìn)行仿真，結(jié)果表明：

（1）基于胞元拓?fù)鋬?yōu)化的方法，得到內(nèi)凹六邊形負(fù)泊松比超材料在底邊長和縱向長度保持不變時，內(nèi)凹角越大其負(fù)泊松比效應(yīng)越為明顯。并且結(jié)構(gòu)胞元按一定周期進(jìn)行排列的宏觀結(jié)構(gòu)的負(fù)泊松比效應(yīng)和結(jié)構(gòu)胞元的負(fù)泊松效應(yīng)改變規(guī)律保持一致。

（2）相較于傳統(tǒng)的對結(jié)構(gòu)胞元的形狀進(jìn)行改變來達(dá)到拓?fù)鋬?yōu)化的結(jié)果，本文采取了對內(nèi)凹六邊形初始角度進(jìn)行改變，將負(fù)泊松比效應(yīng)與內(nèi)凹六邊形的角度用函數(shù)關(guān)系關(guān)聯(lián)了起來，更為有效地確定了負(fù)泊松比效應(yīng)與形狀的具體關(guān)系。

（3）內(nèi)凹六邊形胞元按一定周期排列后整體結(jié)構(gòu)的負(fù)泊松比效應(yīng)將略有減小，在其受到縱向拉（壓）應(yīng)力時，其橫向?qū)l(fā)生膨脹（收縮）變形。橋梁伸縮縫要求其在垂直和平行于橋梁的軸線上都能自由變形，這種負(fù)泊松比結(jié)構(gòu)較為符合其需求。拓?fù)鋬?yōu)化可以得到大量的負(fù)泊松比超材料結(jié)構(gòu)，對負(fù)泊松比超材料結(jié)構(gòu)的獲取和強化有著極為重要的作用。