沈建邦 肖俊華

燕山大學(xué)河北省重型裝備與大型結(jié)構(gòu)力學(xué)可靠性重點實驗室, 秦皇島, 066004

0 引言

負泊松比材料和結(jié)構(gòu)受到軸向壓縮時會發(fā)生橫向收縮，或受到軸向拉伸時會發(fā)生橫向膨脹，EVANS等[1]將這一類具有負泊松比效應(yīng)的材料和結(jié)構(gòu)命名為auxetics(意為拉伸膨脹)，因此負泊松比材料又被稱為拉脹材料。國內(nèi)外學(xué)者已對負泊松比材料和結(jié)構(gòu)進行了廣泛和深入的研究，研究發(fā)現(xiàn)，負泊松比材料和結(jié)構(gòu)壓痕阻力大[2]、抗沖擊性能好[3]、吸聲性能良好[4],且具有同向彎曲性能[5]和抗裂紋擴展性[6]等優(yōu)良性能，被廣泛應(yīng)用于汽車、航空、船舶等工業(yè)領(lǐng)域。典型的負泊松比材料和結(jié)構(gòu)包括泡沫聚合物、蜂窩結(jié)構(gòu)、聚丙烯纖維及某些晶體結(jié)構(gòu)等[7]。本文所研究的負泊松比蜂窩結(jié)構(gòu)不但具有質(zhì)量小、比強度高、結(jié)構(gòu)簡單和可設(shè)計性強等傳統(tǒng)蜂窩結(jié)構(gòu)的優(yōu)良性能[8]，還具有其他獨特的力學(xué)性能，近年來已成為一個新興的研究熱點。

LAKES[9]在實驗室制造出了第一個人工負泊松比材料。CHOI等[10]對材料負泊松比效應(yīng)的產(chǎn)生機理、微觀結(jié)構(gòu)形成方法和結(jié)構(gòu)受外力作用時的變形情況進行了系統(tǒng)的研究。NKANSAH等[11]發(fā)現(xiàn)了聚四氯乙烯具有負泊松比特性。LAKES等[12]采用對傳統(tǒng)三維結(jié)構(gòu)進行轉(zhuǎn)變的方法，指出了三維內(nèi)凹單元的變形行為。YANG等[13]在微極彈性理論的基礎(chǔ)上，對一種二維三角蜂窩結(jié)構(gòu)的幾何參數(shù)與泊松比之間的關(guān)系進行了研究。WAN等[14]基于大撓度理論對負泊松比內(nèi)凹蜂窩結(jié)構(gòu)的力學(xué)性能進行了研究。盧子興等[15]基于旋轉(zhuǎn)機制建立了一種具有負泊松比效應(yīng)的、由部分內(nèi)凹及部分規(guī)則六邊形組成的二維多胞材料力學(xué)模型。GASPAR[7]基于平均線應(yīng)變假設(shè)和一階非均勻性對新型非常規(guī)負泊松比顆粒狀材料進行了研究。顏芳芳等[16]將柔性蜂窩結(jié)構(gòu)應(yīng)用于無人機外翼，探討了蜂窩結(jié)構(gòu)各項參數(shù)的改變對蜂窩結(jié)構(gòu)力學(xué)性能的影響，進而對無人機外翼的性能進行了優(yōu)化。SUN等[17]提出了一種具有負泊松比子結(jié)構(gòu)的多功能分層蜂窩結(jié)構(gòu)，基于歐拉梁理論推導(dǎo)出了結(jié)構(gòu)的彈性模量。郭磊[18]基于旋轉(zhuǎn)剛性單元模型建立了新型多尺寸剛性矩形單元組合模型。魯超等[19]利用柔性懸臂梁模型，給出了蜂窩壁板大變形條件下面內(nèi)等效彈性模量的理論計算公式。金愛兵等[20]從微觀結(jié)構(gòu)角度分析了巖體特殊的負泊松比效應(yīng)，研究了具有負泊松比效應(yīng)的晶體材料。LI等[21]設(shè)計了一種新型負泊松比空心骨架胞元結(jié)構(gòu)，采用有限元法與實驗方法對結(jié)構(gòu)進行了研究。蔣偉等[22]改進了傳統(tǒng)蜂窩材料的結(jié)構(gòu)設(shè)計, 提出了一種環(huán)形負泊松比結(jié)構(gòu), 運用能量法推導(dǎo)出了環(huán)形蜂窩芯等效彈性模量公式，并采用有限元法對環(huán)形蜂窩芯結(jié)構(gòu)進行了力學(xué)性能仿真。LU等[23]設(shè)計了一種在內(nèi)凹蜂窩結(jié)構(gòu)中增加一個窄肋的新型蜂窩結(jié)構(gòu)。FU等[24]依據(jù)梁理論推導(dǎo)出了一種新型負泊松比手性三維材料的等效彈性模量和泊松比的解析解。

本文基于環(huán)形蜂窩結(jié)構(gòu)[22]提出了一種具有負泊松比效應(yīng)的可變弧角曲邊內(nèi)凹蜂窩結(jié)構(gòu)，該結(jié)構(gòu)的曲邊內(nèi)凹程度可以改變，當(dāng)曲邊弧角為180°時即為環(huán)形蜂窩結(jié)構(gòu)。利用能量法推導(dǎo)出了結(jié)構(gòu)的等效彈性模量和等效泊松比的解析表達式，并與已有理論解和有限元結(jié)果進行了比較；討論了等效彈性模量和等效泊松比與結(jié)構(gòu)幾何參數(shù)之間的關(guān)系。

1 結(jié)構(gòu)建模

圖1為本文所提出的曲邊內(nèi)凹蜂窩結(jié)構(gòu)示意圖，其構(gòu)思源于環(huán)形蜂窩結(jié)構(gòu)[22]，所提結(jié)構(gòu)的曲邊弧角θ可以改變。由圖1可知，胞元結(jié)構(gòu)在面內(nèi)呈中心對稱，豎直壁板長度為2h，兩豎直壁板間寬度為2l，曲邊圓弧半徑為r，彎曲壁板為圓的部分圓弧，其對應(yīng)的曲邊弧角為θ，連接部分長度為m，壁板厚度為t，結(jié)構(gòu)垂直于所在平面的厚度為d，當(dāng)弧角θ取180°時，l=r。胞元受到橫向或縱向載荷作用時，結(jié)構(gòu)會產(chǎn)生變形。當(dāng)作用載荷較小時，壁板發(fā)生彈性變形，結(jié)構(gòu)的彈性模量和泊松比等力學(xué)性能由胞元的參數(shù)確定；當(dāng)作用載荷較大時，壁板發(fā)生塑性變形，此時結(jié)構(gòu)的力學(xué)性能隨載荷大小變化。本文采用解析方法研究胞元內(nèi)各幾何參數(shù)對結(jié)構(gòu)宏觀等效力學(xué)性能的影響。

圖1 曲邊內(nèi)凹蜂窩胞元結(jié)構(gòu)示意圖Fig.1 Schematic diagram of honeycomb cell structure with curved concave sides

2 變形分析

2.1 Y軸方向受載

根據(jù)文獻[22]的思路，本文對曲邊內(nèi)凹蜂窩模型進行受力分析(圖2)。令胞元結(jié)構(gòu)在上下胞元連接點處作用圖2a所示的豎直方向的對稱載荷FS，使得胞元內(nèi)部各壁板產(chǎn)生拉伸、剪切和彎曲變形。本文研究的結(jié)構(gòu)為梁架結(jié)構(gòu)，結(jié)構(gòu)中的拉伸應(yīng)變能和剪切應(yīng)變能相對于彎曲應(yīng)變能很小，因此采用能量法求解結(jié)構(gòu)各方向的位移變形時，只考慮結(jié)構(gòu)的彎曲應(yīng)變[22]。因不考慮拉伸應(yīng)變能，故將胞元間的連接桿視為剛性桿，桿兩端的位移相等，可直接研究胞元主體。對稱載荷直接作用在E、F處，胞元結(jié)構(gòu)為封閉梁結(jié)構(gòu)。該力學(xué)問題為超靜定問題，且封閉梁結(jié)構(gòu)關(guān)于中心對稱，因此可取結(jié)構(gòu)的1/4部分進行研究，如圖2b所示，自由端受豎直載荷FS/2和未知彎矩M0的作用。

(a)整體胞元結(jié)構(gòu)受載示意圖 (b)1/4結(jié)構(gòu)的內(nèi)力示意圖

(c)只作用F/2載荷時1/4結(jié)構(gòu) (d)只作用單位彎矩時1/4結(jié)構(gòu)

(e)只作用水平單位力時1/4結(jié)構(gòu)圖2 Y軸方向受載時的受力分析Fig.2 Force analysis under Y-axis load

因彎矩M0未知，故將其在下式中記作MR1，則自由端H處截面的轉(zhuǎn)角變形協(xié)調(diào)條件滿足：

δ11MR1+δ1F=0

(1)

式中,δ1F為僅作用載荷FN=FS/2時H所在截面的轉(zhuǎn)角，如圖2c所示；δ11為在H處作用單位彎矩時H所在截面的轉(zhuǎn)角,如圖2d所示。

(2)

(3)

Im=dt3/12

式中,Em為材料彈性模量;Im為截面對中性軸的慣性矩。

由式(1)可得

(4)

由式(4)可得

(5)

運用卡氏定理可求出在Y軸方向作用對稱載荷FS時，胞元結(jié)構(gòu)在Y軸方向的位移變形ΔYY，即為結(jié)構(gòu)應(yīng)變能對Y軸方向作用力FS求偏導(dǎo)，其表達式如下:

(6)

(7)

2.2 X軸方向受載

如圖3所示，令結(jié)構(gòu)只受X方向的水平對稱載荷FS作用，取其1/4結(jié)構(gòu)研究，可分別求得在X軸和Y軸方向的位移變形，其表達式分別如下：

圖3 X軸方向受力時結(jié)構(gòu)示意圖Fig.3 Schematic diagram of structure under X-axis load

(8)

(9)

3 等效彈性模量與等效泊松比

依據(jù)Y軸受力時的位移變形公式(式(6)和式(7))，可求得Y軸方向的等效彈性模量和等效泊松比分別為

(10)

(11)

式中,εYY、εYX分別為Y方向受力時Y方向和X方向的應(yīng)變；σY為Y方向應(yīng)力。

同理可求得X軸方向的彈性模量和等效泊松比分別為

(12)

(13)

式中，εXX、εXY分別為X方向受力時X方向和Y方向的應(yīng)變；σX為X方向應(yīng)力。

4 結(jié)果分析與討論

4.1 本文解析解與已有理論解及有限元結(jié)果比較

本文所提出的可變弧角曲邊內(nèi)凹蜂窩結(jié)構(gòu)(圖1)，在弧角θ=180°時將退化為環(huán)形蜂窩結(jié)構(gòu)[22]。有限元計算模型采用5×5陣列的周期胞元結(jié)構(gòu)，如圖4所示，材料選取較為常見的鋁合金，彈性模量Em=71 GPa，密度ρ=2 770 kg/m3，泊松比ν=0.33。算例中厚度t取定值1 mm。設(shè)置位移約束條件如下：①約束與X軸平行的結(jié)構(gòu)對稱中心線Y方向的位移,②約束與Y軸平行的對稱中心線X方向的位移,③約束整個結(jié)構(gòu)Z方向的位移。分別作用X軸方向?qū)ΨQ載荷與Y軸方向?qū)ΨQ載荷，求得結(jié)構(gòu)位移變形后，計算等效彈性模量。

圖4 有限元計算模型Fig.4 Finite element calculation model

選取與文獻[22]相同的幾何參數(shù)，圖5給出了環(huán)形蜂窩結(jié)構(gòu)時的Y軸方向等效彈性模量EY的本文解析解、文獻[22]理論解以及有限元結(jié)果的比較曲線圖。從圖5中可以看出，本文解析解、文獻[22]理論解和有限元結(jié)果三者完全吻合，驗證了本文理論推導(dǎo)和解析結(jié)果的正確性,同時也表明了本文在采用能量法求解時忽略拉伸應(yīng)變能和剪切應(yīng)變能這一假設(shè)的合理性。

圖5 環(huán)形蜂窩結(jié)構(gòu)時Y軸方向等效彈性模量與曲邊圓弧半徑的關(guān)系曲線Fig.5 The relationship between Y-axis equivalent elastic modulus and curved arc radius in circular honeycomb structure

4.2 弧角θ對等效彈性模量和等效泊松比的影響

當(dāng)弧角θ取為任意角度時，尚無其他理論解答可以進行對比。為校驗本文理論推導(dǎo)和解析結(jié)果的有效性，在研究弧角變化對結(jié)構(gòu)等效彈性模量和等效泊松比的影響時，同時給出了本文解析解和有限元結(jié)果。

(a)Y軸方向等效彈性模量

(b)X軸方向等效彈性模量圖6 等效彈性模量隨弧角的變化Fig.6 Variations of the equivalent elastic modulus with arc angles

(a)Y向受載時等效泊松比

(b)X向受載時等效泊松比圖7 等效泊松比隨弧角的變化Fig.7 Variations of the equivalent Poisson’s ratio with arc angles

4.3 結(jié)構(gòu)各參數(shù)對等效泊松比的影響

利用所得理論公式研究結(jié)構(gòu)中幾何參數(shù)對等效泊松比的影響規(guī)律。值得注意的是，不同弧角下的h和l中某一參數(shù)固定時，另一參數(shù)取值范圍有所限制，因此在各圖中不同曲線的橫坐標起止值不同。

(a)X向受載時等效泊松比

(b)Y向受載時等效泊松比圖8 等效泊松比隨豎直壁板半長的變化Fig.8 Variations of the equivalent Poisson’s ratio with half-length of vertical wall

(a)X向受載時等效泊松比

(b)Y向受載時等效泊松比圖9 等效泊松比隨兩豎直壁板間半寬的變化Fig.9 Variations of the equivalent Poisson’s ratio with half-width between two vertical walls

(a)X向受載時等效泊松比

(b)Y向受載時等效泊松比圖10 等效泊松比隨連接部分長度的變化Fig.10 Variations of the equivalent Poisson’s ratio with connected past lengths

由上述討論可知，改變結(jié)構(gòu)參數(shù)時，結(jié)構(gòu)的等效泊松比既能取正值也能取負值，不同方向的等效泊松比隨結(jié)構(gòu)幾何參數(shù)的改變有較大差異。

5 結(jié)論