沿板平面材料梯度變化的二維梯度薄板撓度分析

許楊健，龐宇飛，劉 碩

沿板平面材料梯度變化的二維梯度薄板撓度分析

許楊健1，龐宇飛1，劉 碩2

(1.河北工程大學 土木工程學院，河北 邯鄲 056038；2.哈爾濱工業大學 航天學院， 黑龍江 哈爾濱 150001)

采用有限單元法，分析了由Ti-6Al-4V、Al 1100和ZrO2組成的二維功能梯度材料（2D-FGM）薄板的小撓度問題，計算了在不同溫度環境下，四邊簡支受均布荷載的2D-FGM薄板撓度值和四邊固支受均布荷載的2D-FGM薄板撓度值，并對本方法進行了正確性檢驗，給出了2D-FGM薄板撓度分布圖。進一步討論了溫度和孔隙率對2D-FGM薄板撓度的影響。結果表明，溫度和孔隙率對受載2D-FGM薄板撓度變形影響較為明顯，在環境溫度不變的情況下，四邊固支和四邊簡支的2D-FGM薄板撓度都隨孔隙率的增大而逐漸變大，在孔隙率不變情況下，四邊固支和四邊簡支的2D-FGM薄板撓度都隨環境溫度的不斷升高而逐漸變大。

功能梯度材料薄板；撓度；溫度；孔隙率；有限元法

FGM是材料科學領域發展起來的一種新一代非均勻復合材料，它是各國在航空航天高技術領域里的競爭產物。薄板構件在機械、土建、造船、航空等工程領域中應用都很廣泛，因此，研究薄板結構性能對各個工程領域都很重要，目前，研究學者們就薄板結構性能做出了大量的研究[1-5]，并獲得了豐富的研究成果。Ping Tan等[6]分析了在彈性約束邊界下的變剛度復合環形薄板的自由和強迫振動。董文堂[7]導出了四邊固支邊界條件下的矩形薄板納維葉解法。Rui Li等[8]得出了基于多點支撐下的矩形薄板自由振動新的解析解。Shilang Xu等[9]對由不同的膠凝材料組成的織物增強混凝土（TRC）薄板進行了耐高溫性能試驗研究。許楊健等[10]采用層合板有限元法對FGM薄板的瞬態熱彈性彎曲問題進行了分析研究。Amin Joodaky[11]對不同邊界條件下FGM薄斜板的撓度和應力進行了分析。孫建等[12]計算出了構造撓度函數的系數。目前，大多數研究只在材料屬性沿板厚度方向變化上，而沿著板平面材料梯度變化的研究相對很少，曹志遠等[13-14]得到了沿板平面材料特性梯度變化的功能梯度板件力學量三維分布狀態。本文基于經典薄板理論，利用有限元方法結合Matlab編程求解出2D-FGM薄板受均布荷載下的撓度值，并分別討論了四邊固支和四邊簡支兩種邊界條件下，溫度和孔隙率對受均布荷載2D-FGM薄板撓度值的影響?？晒┬滦筒牧蠐隙茸冃我罂量痰暮娇蘸教?、精密儀器制造等領域作參考，為FGM的設計和應用提供了理論計算依據。

1 模型建立

如圖1所示，建立由金屬相Ti-6Al-4V、Al 1100和陶瓷相ZrO2組成的2D-FGM矩形薄板模型，考慮材料組分物性值隨溫度變化，假設材料屬性沿x、y軸方向變化，沿薄板厚度無變化。

圖1 2D-FGM 薄板模型Fig.1 Model of 2D-FGM thin plate

2 FGM薄板彎曲變形有限元方程

為了研究薄板的彎曲變形，則采用經典薄板理論，構造矩形薄板單元。矩形薄板單元共有四個節點，每個節點有3個基本未知數：撓度w，法線繞x軸的轉動θx和繞y軸的轉動θy。

根據微小位移的假定，位移可表示為：

單元節點位移向量為：

取位移函數為：

由位移w在薄板單元各個節點處的取值，可以確定w中的12個參量，

式中｛N｝為形函數。

根據虛功原理得：：

式中：｛δ*｝代表虛位移，表示FGM的廣義彈性矩陣。而由于｛δ*｝e具有任意性，故得：

式中：［K］ε表示功能梯度材料薄板的單元剛度矩陣。

3 數值算例與分析

3.1 FGM薄板物性值

根據FGM材料研究現狀，本文FGM薄板采用由金屬相Ti-6Al-4V、Al 1100和陶瓷相ZrO2組成的二維FGM薄板作為研究對象，不考慮組分的影響，考慮物性系數隨溫度的變化，依據文獻[15]，金屬Al 1100、Ti-6Al-4V以及陶瓷ZrO2組分材料的變物性物性系數分別如下：

金屬Al 1100：

金屬Ti-6Al-4V：

陶瓷ZrO2：

表1 四端固支均質材料薄板受均勻荷載時中心點的撓度Tab.1 Def l ection at center point of clamped supported homogeneous sheet under uniformly distributed loads

表2 四端簡支均質材料薄板受均勻荷載時中心點的撓度Tab.2 Def l ection at center point of simply supported homogeneous sheet under uniformly distributed loads

3.2 正確性檢驗

為驗證本文推導的表達式的正確性，取薄板厚度t為4 mm，a、b都取為100 mm，橫向均布荷載取q=0.1 N/mm2，同時將FGM薄板退化為Ti-6Al-4V、Al 1100和ZrO2三種均勻材料薄板，在單元內對材料屬性取平均值。將本文有限元計算程序計算的四端固支均質薄板中心撓度與文獻[7]的3.7式計算結果對比，四端簡支均質薄板的與文獻[16]的3.43式計算結果對比。本文方法和文獻方法計算薄板中心撓度的結果如表1和表2所示：

通過對比兩種方法的計算結果，誤差都不超過±0.8%，則本文結果正確。

3.3 溫度、孔隙率影響

圖2 固支2D-FGM薄板的撓度曲線Fig.2 Deflection curve of clamped supported 2D-FGM thin plate

圖3 簡支2D-FGM薄板的撓度曲線Fig.3 Deflection curve of simply supported 2D-FGM thin plate

圖2 和圖3表明：不同孔隙率、和不同溫度下的四邊固支和四邊簡支2D-FGM薄板的撓度分布圖。分析可知：(1)薄板的撓度值從板的邊緣到中間逐漸變大，且薄板的最大撓度值靠近板的中心處。(2)在不考慮薄板孔隙率變化的情況下，薄板隨環境溫度的增加撓度逐漸增大，并且變化率越來越大。如：在圖2中Ax=0.0, Ay=1.0時，T=500 K薄板中心點處撓度比T=300 K時增加10.7%，而T=700 K薄板中心點處撓度比T=500 K時增加了14.2%。在圖3中Ax=0.0, Ay=1.0時，T=500 K薄板中心點處撓度比T=300 K時增加10.7%，而T=700 K薄板中心點處撓度比T=500 K時增加了14.2%。(3)在不考慮環境溫度變化的情況下，薄板撓度隨孔隙率的增加不斷變大。如：在圖2中T=700 K時，Ay=0.0，Ax=1.0薄板中心點撓度比Ay=0.0, Ax=0.0時增加了8.3%。取T=500 K時，Ax=0.0，Ay=1.0薄板中心點撓度比Ax=0.0，Ay=0.0為時增加了43.9 %。在圖3中T=300 K時，Ay=0.0，Ax=1.0薄板中心點撓度比Ay=0.0，Ax=0.0時增加了7.4 %，T=900 K時，Ax=0.0, Ay=1.0薄板中心點撓度比Ax=0.0，Ay=0.0時增加了45.0%。

4 結論

1）通過對由Ti-6Al-4V、Al 1100和ZrO2組成的2D-FGM薄板算例的撓度計算，計算結果和文獻給出的結果是一致的，說明了本文方法正確可行，給出了不同孔隙率下2D-FGM薄板撓度分布。

2）根據算例計算結果分析可知，無論FGM薄板四端簡支還是四端固支，環境溫度和孔隙率對受載FGM薄板撓度影響均較為明顯，薄板撓度隨著環境溫度的升高而增大且變化率也在增大，隨著孔隙率變大撓度也在逐漸增大。

3）在本文的相同條件下，四邊固支薄板比四邊簡支薄板撓度值小，更適合承受大的荷載。

[1]MIRSALEHI M，AZHARI M，AMOUSHAHI H.Stability of thin FGM microplate subjected to mechanical and thermal loading based on the modified couple stress theory and saline finite strip method[J].Aerospace Science and Technology,2015;47(8)：356-366.

[2]TAN P, NIE GJ.Free and forced vibration of variable stiffness composite annular thin plates with elastically restrained edges[J].Composite Structures,2016;149：398-407.

[3]LIN R, TIAN Y, WANG P C.New analytic free vibration solutions of rectangular thin plates resting on multiple point supports[J]. International Journal of Mechanical Sciences,2016;110：53-61.

[4]IIDA J, HASEBE N.Stress intensity factors of a rhombic hole with symmetric cracks under uniform transverse thin plate bending[J].Engineering Fracture Mechanics,2016,156：16-24.

[5]王春玲,高 典,劉俊卿.橫觀各向同性彈性半空間地基上四邊自由各向異性矩形薄板彎曲解析解[J].力學季刊,2015，36(1)：95-104.

[6]PING TAN, G.J. NIE.Free and forced vibration of variable stiffness composite annular thin plates with elastically restrained edges[J].Composite Structures,2016,149：398-407.

[7]董文堂.固支邊矩形薄板的納維葉解法[J].黃石高等?？茖W校學報, 1999,1(15)： 1-4.

[8]LI R, TIAN Y, WANG P C.New analytic free vibration solutions of rectangular thin plates resting on multiple point supports[J].International Journal of Mechanical Sciences,2016,110：53-61.

[9]XU S L, SHEN L H, WANG J Y.The high-temperature resistance performance of TRC thin-plates with different cementitious materials： Experimental study[J]. Construction and Building Materials,2016,115：506-519.

[10]許楊健, 趙志崗.梯度功能材料薄板瞬態熱彈性彎曲有限元分析[J].工程力學,2001,18(1)：71-81.

[11]JOODAKY A, JOODAKY I, HEDAYATI M.Deflection and stress analysis of thin FGM skew plates on Winkler foundation with various boundary conditions using extended Kantorovich method[J].Composites： Part B,2013;51(3)：191-196.

[12]孫 建,胡 洋.均布和靜水壓力作用下固支矩形薄板力學特性[J].應用力學學報，2015,32(6)： 908-914.

[13]曹志遠,程紅梅.沿板平面材料梯度變化功料板件分析[J].力學季刊,2007,28(2)：203-208.

[14]曹志遠,程紅梅.沿板平面材料組分變異功梯度板件分析[J].同濟大學學報,2007,35(11)：1460-1465.

[15]王 飛.變物性二維梯度板熱彈性行為[D]. 邯鄲：河北工程大學, 2014.

[16]曲慶璋.彈性板理論[M].北京：人民交通出版社出版,1999.

（責任編輯 王利君）

Deflection analysis of two-dimensional graded thin plates varying along plane direction

XU Yangjian1,PANG Yufei1,LIU Shuo2

(1 College of Civil Engineering, Hebei University of Engineering, Hebei Handan, 056038;2 School of Astronautics, Harbin Institute of Technology, Heilongjiang Harbin, 150001 )

The small deflection of the two dimensional functionally graded material (2D-FGM) thin plate (Ti-6Al-4V, Al 1100 and ZrO2) is analyzed by the finite element method(FEM) .The deflection value of 2D-FGM thin plate subjected to uniform distributed load with four edges simply-supported or clampedsupported is calculated under different temperatures. At the same time, the accuracy of this method is examined. The deflection distribution diagrams are obtained. Effects of temperature and porosity on the deflection of 2D-FGM thin plate are discussed. The results show that the influence of temperature and porosity on the load 2D-FGM thin plate is obvious. In the case of constant ambient temperature, the deflection of 2D-FGM thin plate with four edges clamped-supported and four edges simply-supported increases with the increasing of porosity. Under the condition of constant porosity, the deflection of 2D-FGM thin plate with four edges clamped-supported and four edges simply-supported increases with the increasing of temperature.

functionally gradient material thin plate; deflection; temperature; porosity; finite element method

TB333

A

1673-9469（2017）01-0034-04

10.3969/j.issn.1673-9469.2017.01.008

2016-09-08 特約專稿

河北省自然科學基金資助項目 (E2011402033)

許楊健(1956 -)，男，四川富順人，教授，從事固體力學計算、有限元應用和梯度功能材料方面的研究。