張文福，沙彥文

（1.南京工程學(xué)院建筑工程學(xué)院，江蘇南京 211167；2.南京工程學(xué)院機(jī)械工程學(xué)院，江蘇南京 211167）

目前，在機(jī)械工程領(lǐng)域，薄壁結(jié)構(gòu)的應(yīng)用十分廣泛。但是，薄壁結(jié)構(gòu)也有一定的缺點(diǎn)，由于截面較薄，在承受外力時(shí)容易發(fā)生變形。功能梯度材料（Functionally Gradient Materials，F(xiàn)GM）是由兩種或多種材料復(fù)合且成分和結(jié)構(gòu)呈連續(xù)梯度變化的一種新型復(fù)合材料，其概念最早在1984 年由日本科學(xué)家平井敏雄等提出[1]。近年來(lái)薄壁構(gòu)件逐漸受到重視，在材料制備中可采用功能梯度材料，通過(guò)將陶瓷、塑料、金屬等巧妙組合，充分發(fā)揮各類(lèi)材料的優(yōu)勢(shì)。

在薄壁結(jié)構(gòu)中，相比于開(kāi)口截面，閉口截面具有更好的抗扭和抗剪性能，但目前針對(duì)閉口截面薄壁構(gòu)件的研究多集中在箱型截面，對(duì)其他類(lèi)型的閉口截面如三角形截面的研究較少，且多集中在三角形桁架梁或桁架拱結(jié)構(gòu)。在靜力研究方面：文穎等人[2]提出了混合截面的約束扭轉(zhuǎn)有限元模型；何海玉等人[3]提出了倒三角形截面的鋼圓弧拱在全跨均布載荷作用下的極限承載力公式；楊則英等人[4]提出了計(jì)算平面曲線(xiàn)梁溫度變形的簡(jiǎn)單算法。在動(dòng)力研究方面：Dikaros等人[5]提出了任意單連接或多連接閉口截面梁的非均勻翹曲動(dòng)力學(xué)分析的一般公式；Wang 等人[6]確定了彎曲箱梁的7個(gè)基本位移，在曲線(xiàn)坐標(biāo)系下建立了一個(gè)三節(jié)點(diǎn)曲梁有限元，用于曲箱梁的振動(dòng)特性和動(dòng)力響應(yīng)分析；Carrera 等人[7]采用改進(jìn)的梁理論，對(duì)疊層復(fù)合材料箱形結(jié)構(gòu)的自由振動(dòng)進(jìn)行了討論，采用Carrera統(tǒng)一公式建立了高階模型；Jiang等人[8]建立了薄壁箱梁的改進(jìn)單元?jiǎng)偠群唾|(zhì)量矩陣。

近年來(lái)，功能梯度材料的國(guó)內(nèi)外研究成果逐漸增多。在靜力研究方面：Avhad等人[9]采用Navier解法得到了功能梯度組合梁在高度上彎曲的應(yīng)力和撓度；Sayyad等人[10]對(duì)功能梯度多孔圓形梁的靜態(tài)變形和振動(dòng)頻率進(jìn)行了分析，獲得了靜態(tài)和振動(dòng)問(wèn)題的Navier型封閉解；Sayyad等人[11]則提出了一種多項(xiàng)式型五階彎曲梁理論。在動(dòng)力研究方面：江希等人[12]利用Pei 研究的針對(duì)功能梯度材料梁的新型高階梁理論，將其由靜力學(xué)領(lǐng)域拓展到梁的模態(tài)分析,并通過(guò)構(gòu)造新型高階單元，采用有限元方法研究功能梯度梁的前三階自由振動(dòng)；Sayyad等人[13]提出了一種新的針對(duì)功能梯度曲梁自由振動(dòng)分析的高階剪切法向變形理論；Sharma等人[14]研究了軸向錐形功能梯度（AFG）梁的固有頻率。

然而，作為一種新型的復(fù)合材料，功能梯度材料現(xiàn)有的研究多集中在材料性質(zhì)、壽命及制造方面，由其制造的薄壁構(gòu)件的力學(xué)性能還需進(jìn)一步探討。本文基于板-梁理論，按材料性質(zhì)沿板件厚度梯度變化情形，對(duì)三角形截面懸臂梁的組合扭轉(zhuǎn)進(jìn)行研究；利用有限元軟件ANSYS 建立了相應(yīng)模型，并對(duì)研究結(jié)果進(jìn)行了驗(yàn)證。

1 FGM 三角形截面梁組合扭轉(zhuǎn)理論

1.1 基本假設(shè)

1.1.1 剛周邊假設(shè)

板件的橫截面在其面內(nèi)的剛度是無(wú)限大的，因此在構(gòu)件發(fā)生扭轉(zhuǎn)時(shí)，橫截面在其軸向投影的形狀保持不變。

1.1.2 變形假設(shè)

將板件的變形分解為平面內(nèi)變形和平面外變形兩部分，這樣可以簡(jiǎn)化分析。

1.1.3 板-梁假設(shè)

板件的平面內(nèi)、平面外縱向位移和應(yīng)變能分別由Timoshenko梁理論和Kirchhoff板理論確定。

1.2 計(jì)算簡(jiǎn)圖

如圖1所示等腰三角形中：bf為上翼緣寬度，tf為上翼緣厚度，bw為腹板寬度，tw為腹板厚度，C 為截面的形心，S 為截面的剪心，梁長(zhǎng)為L(zhǎng)。o-xyz 為截面整體坐標(biāo)系，原點(diǎn)位于形心，x 軸與翼緣平行，y 軸為截面的對(duì)稱(chēng)軸。o-nsz 為各塊板的局部坐標(biāo)系，均符合右手螺旋法則。hs1為上翼緣自身形心到剪心的距離，hs2為腹板底部到剪心的距離，ef1為上翼緣自身形心到整個(gè)截面形心的距離，ew為腹板形心到整個(gè)截面形心的距離，β為等腰三角形對(duì)稱(chēng)軸與腹板的夾角。

圖1 三角形截面尺寸圖

如圖2 所示，為三角形截面扭轉(zhuǎn)變形前后對(duì)比圖。

圖2 三角形截面扭轉(zhuǎn)變形圖

圖3 材料性質(zhì)沿板件厚度方向變化示意圖

1.3 截面的形心與剪心

根據(jù)文獻(xiàn)[15]，得到公式（1）和公式（2）。（1）截面的形心

其中，Aw和Af為三角形截面腹板和上翼緣的面積。

（2）截面的剪心

其中，Iy,w和Iy,f為三角形截面腹板和翼緣的慣性矩。

1.4 材料性質(zhì)沿厚度變化的組合扭轉(zhuǎn)理論

1.4.1 左腹板截面任意點(diǎn)位移

根據(jù)變形分解原理，可將左腹板任一點(diǎn)位移寫(xiě)為：

其中：等號(hào)后第1 項(xiàng)為左腹板平面內(nèi)的橫向位移；第2項(xiàng)為左腹板平面外的橫向位移。

1.4.2 左腹板平面內(nèi)彎曲應(yīng)變能

應(yīng)變能公式為：

其中，E(n)和G(n)分別為功能梯度材料的彈性模量和剪切模量，兩者關(guān)系為：

式中μ為泊松比，且彈性模量的變化模式為：

式中：pw=1，2，3…，為腹板功能梯度材料的梯度指數(shù)；Emw為腹板厚度方向中部的彈性模量；Eow為腹板厚度方向外邊緣的彈性模量；mw=EowEmw=GowGmw為腹板外邊緣與中部的模量比。

由平面內(nèi)位移的表達(dá)式（3）可得左腹板內(nèi)某一點(diǎn)的位移為：

沿著左腹板n 軸的位移

沿著左腹板s軸的位移

縱向位移根據(jù)Timoshenko梁理論可以表示為

其中：ψw1是閉口截面構(gòu)件中因剪切變形導(dǎo)致的截面轉(zhuǎn)角，后面需要用截面扭轉(zhuǎn)角θ 加以確定；b1是因左腹板的形心不在截面主軸上而引起的縱向位移。

線(xiàn)性應(yīng)變?yōu)椋?/p>

根據(jù)應(yīng)變能公式（4）及線(xiàn)性應(yīng)變公式（5）（6）（7），則左腹板平面內(nèi)應(yīng)變能為：

式中χ3是考慮功能梯度材料模量變化影響的修正系數(shù)，其表達(dá)式為：

1.4.3 左腹板平面外彎曲應(yīng)變能

同上方法，則左腹板平面外應(yīng)變能為：

式中χ4是考慮功能梯度材料模量變化影響的修正系數(shù)，其表達(dá)式為：

綜上，左腹板扭轉(zhuǎn)應(yīng)變能為：

由對(duì)稱(chēng)性，左右腹板的扭轉(zhuǎn)應(yīng)變能相同，即Uw1=Uw2。同上方法也可得到上翼緣的應(yīng)變能為：

式中χ1、χ2是考慮功能梯度材料模量變化影響的修正系數(shù)，其表達(dá)式為：

1.4.4 截面扭轉(zhuǎn)角與橫截面轉(zhuǎn)角的關(guān)系

在上述三角形截面的應(yīng)變能計(jì)算中，根據(jù)文獻(xiàn)[15]，截面扭轉(zhuǎn)角θ 與截面轉(zhuǎn)角ψf、ψw1、ψw2的關(guān)系如下：

1.4.5 截面總應(yīng)變能及總勢(shì)能

三角形截面總應(yīng)變能為：

當(dāng)翼緣和腹板的材料性質(zhì)變化一致，即Emf=Emw=Em時(shí)，應(yīng)變能表達(dá)式簡(jiǎn)化為：

式中：

為截面的翹曲慣性矩；

為截面的自由扭轉(zhuǎn)常數(shù)；

為考慮功能梯度材料模量變化影響的修正系數(shù)。

若在自由端施加集中扭矩Mt，則外力勢(shì)能為：

截面的總勢(shì)能為：

1.4.6 能量微分方程模型

依據(jù)能量變分原理，由δΠ=0可得：

采用分部積分法最終得到：

可得如下微分方程：

1.4.7 懸臂梁扭轉(zhuǎn)角的精確解析解

對(duì)式（13）微分方程模型求解，并結(jié)合邊界條件可得扭轉(zhuǎn)角表達(dá)式為：

所以，最大扭轉(zhuǎn)角即自由端的扭轉(zhuǎn)角為：

其中，K 為薄壁構(gòu)件的扭轉(zhuǎn)剛度參數(shù)（無(wú)量綱），是反映自由扭轉(zhuǎn)扭矩和翹曲扭轉(zhuǎn)在總扭矩中所占比例的重要系數(shù)。

2 FGM 三角形截面梁組合扭轉(zhuǎn)有限元驗(yàn)證

以圖4 所示的材料性質(zhì)沿板件厚度變化三角形截面懸臂梁扭轉(zhuǎn)角為例，進(jìn)行有限元驗(yàn)證。

圖4 懸臂梁在自由端作用集中扭矩

2.1 有限元模型的建立

2.1.1 基本假定

構(gòu)件的變形是微小的，即小變形；材料始終處于彈性階段；構(gòu)件發(fā)生扭轉(zhuǎn)時(shí)可以翹曲，但橫截面軸向投影的形狀不發(fā)生變化。

2.1.2 邊界條件

因?yàn)楸疚难芯康氖菓冶哿旱慕M合扭轉(zhuǎn)，所以約束條件為：在固定端約束所有自由度；在自由端施加集中扭矩；在每個(gè)截面處使用CERIG 命令建立剛性域，保證截面上所有節(jié)點(diǎn)的扭轉(zhuǎn)位移一致。

2.1.3 單元選擇

在A(yíng)NSYS中模擬功能梯度材料可以通過(guò)分層來(lái)解決，賦予每層不同的材料性質(zhì)，包括彈性模量E 和剪切模量G。ANSYS 中的SHELL99 可以模擬平板和曲殼一類(lèi)結(jié)構(gòu)，也適用于模擬薄壁結(jié)構(gòu)，一個(gè)單元最多可以分成250層，從而能夠較為精確地模擬功能梯度材料。

2.1.4 建模及求解過(guò)程

在SHELL99 中使用分層技術(shù)且每層等厚來(lái)模擬功能梯度材料。將第2 個(gè)關(guān)鍵選項(xiàng)（KEYOPT）設(shè)定為0，然后通過(guò)改變實(shí)常數(shù)（RMODIF）和循環(huán)命令（DO-ENDDO），賦予每層單元不同的性質(zhì)。

在建模時(shí)，翼緣和腹板選用同樣類(lèi)型的功能梯度材料：鋁合金-陶瓷-鋁合金型。鋁合金彈性模量為Em=70 GPa，陶瓷彈性模量為Eo=380 GPa，泊松比μ=0.3，梯度指數(shù)為2。彈性模量變化形式為二次函數(shù)形式：即板件的厚度方向中心為鋁合金，而最外部的兩側(cè)為陶瓷。

在對(duì)每層賦予材料性質(zhì)后，進(jìn)行幾何建模，施加集中扭矩1 N·m，最后求解。為保證薄壁結(jié)構(gòu)剛周邊特性，使用剛域（CERIG）命令，可以保證截面在扭轉(zhuǎn)時(shí)形狀不變。此外，當(dāng)懸臂梁在自由端施加扭矩進(jìn)行組合扭轉(zhuǎn)驗(yàn)證時(shí)，為避免自由端發(fā)生異常的過(guò)度變形，可以使用接觸單元（TARGE170和CONTA175）來(lái)保證模擬的正確性。

2.2 有限元驗(yàn)證

如圖5 所示，根據(jù)公式（17）計(jì)算了6 組不同尺寸的懸臂梁自由端扭轉(zhuǎn)角，并建立了相應(yīng)的有限元模型。如表1 所示，為理論值和有限元值的對(duì)比數(shù)據(jù)。

表1 材料性質(zhì)沿厚度變化的三角形截面懸臂梁扭轉(zhuǎn)角對(duì)比數(shù)據(jù)

圖5 三角形懸臂梁組合扭轉(zhuǎn)有限元模擬圖

3 結(jié)論

（1）基于板-梁理論，將材料性質(zhì)按照沿板件厚度變化的情況，推導(dǎo)了三角形截面懸臂梁組合扭轉(zhuǎn)的總勢(shì)能。建立能量變分模型和能量微分模型，求解得到三角形截面懸臂梁在端部扭矩作用下最大扭轉(zhuǎn)角的表達(dá)式。

（2）利用有限元軟件ANSYS 對(duì)公式進(jìn)行了驗(yàn)證，結(jié)果表明：轉(zhuǎn)角誤差在2.77%~4.29%。故公式推導(dǎo)具有一定的合理性，可為工程設(shè)計(jì)提供參考。