滕兆春 ，王偉斌，薛剛

（蘭州理工大學(xué) 理學(xué)院，甘肅 蘭州 730050）

功能梯度材料（Functionally Graded Materials，F(xiàn)GM）是一種新型的非均勻復(fù)合材料，其作為智能化材料的典型代表，在工程領(lǐng)域內(nèi)應(yīng)用廣泛［1］.隨著新型材料的興起，F(xiàn)GM 圓板/環(huán)板作為核工業(yè)、航天、機(jī)械和土木等工程中常見的重要功能結(jié)構(gòu)元件，在工程實(shí)際中具有較多的應(yīng)用背景，尤其是應(yīng)用于高溫工況中.為了準(zhǔn)確評(píng)估FGM 結(jié)構(gòu)的機(jī)械行為和為工程設(shè)計(jì)提供理論依據(jù)，許多學(xué)者先后展開了大量的研究工作.

目前，關(guān)于FGM 圓板結(jié)構(gòu)的力學(xué)行為研究已有較多的研究成果，具體可參考Boutahar 等［2］、Allahverdizadeh 等［3］、Ma 等［4］和李世榮等［5］諸多學(xué)者的系列研究工作.近年來，許多研究者對(duì)FGM 圓板/環(huán)板結(jié)構(gòu)的靜動(dòng)態(tài)響應(yīng)分析，如屈曲和振動(dòng)方面做了大量的工作，可見文獻(xiàn)［6-7］.由于FGM 圓板/環(huán)板是大多數(shù)結(jié)構(gòu)的關(guān)鍵元件，以下文獻(xiàn)集中利用不同理論和不同方法對(duì)FGM圓板/環(huán)板力學(xué)行為進(jìn)行研究.如Hosseini-Hashemi 等［8］基于一階剪切變形Mindlin 理論，給出了階梯厚度圓形和環(huán)形功能梯度板自由振動(dòng)的精確解.利用分離變量法求解了其自由振動(dòng)，研究了FGM 板的幾何參數(shù)和材料參數(shù)，如階梯厚度比、階梯位置和梯度指數(shù)對(duì)固有頻率的影響.Asemi等［9］首次采用三維彈性理論，研究了FGM 環(huán)形、扇形板的剪切屈曲問題，詳細(xì)研究了載荷、邊界條件、梯度指數(shù)和扇形角等對(duì)FGM 環(huán)形、扇形板剪切屈曲載荷和振型的影響.Malekzadeh 等［10］基于三維彈性理論，研究了熱環(huán)境下FGM 厚環(huán)板的自由振動(dòng)分析，應(yīng)用微分求積法（DQM）求解了熱彈性平衡方程和自由振動(dòng)方程，研究了溫升、材料和幾何參數(shù)對(duì)固有頻率的影響.Tajeddini 等［11］基于線性、小應(yīng)變和精確彈性理論，分析了Pasternak 地基上徑向變厚度的厚各向同性材料和FGM 圓形及環(huán)形板的三維自由振動(dòng)問題，并利用多項(xiàng)式Ritz法求解了特征值問題.

在此研究基礎(chǔ)上，考慮到FGM 構(gòu)件的實(shí)際應(yīng)用背景和材料自身情況，如FGM 微納米器件和多孔FGM 等，許多研究者相繼將微觀理論和多孔材料特性等應(yīng)用于FGM 圓板開展研究，如Goodarzi 等［12］利用修正應(yīng)變梯度理論（MSGT）和修正偶應(yīng)力理論（MCST），研究了黏彈性地基上圓形和環(huán)形納米板的熱-機(jī)振動(dòng)問題，使用微分求積法（DQM）和伽遼金法求解控制方程，研究了尺寸相關(guān)、平面預(yù)緊力、溫度變化和彈性介質(zhì)等因素對(duì)固有頻率的影響.Eshraghi等［13］給出了統(tǒng)一表示的位移場(chǎng)，基于修正的偶應(yīng)力理論介紹了熱載FGM 環(huán)形和圓形微板的靜力彎曲和自由振動(dòng)問題的求解方法，分析了熱載荷、材料和幾何參數(shù)對(duì)靜態(tài)變形、應(yīng)力和自振頻率的影響.Jabbari 等［14］研究了具有壓電驅(qū)動(dòng)器層的多孔材料徑向?qū)嵭膱A板的熱屈曲，利用變分得到了壓電多孔板的控制方程，給出了圓板在溫度載荷作用下的封閉解，分析了多孔板厚度、孔分布、壓電厚度、外加驅(qū)動(dòng)電壓和孔隙率變化對(duì)臨界溫度載荷的影響.Behravan等［15］基于三維彈性理論，分析了受非均勻面力和非對(duì)稱Kerr 彈性地基上非軸對(duì)稱變厚度FGM 多孔圓板的磁-彈問題，基于徑向和厚度方向上的微分求積和狀態(tài)空間矢量技術(shù)，討論了面力、邊界條件和彈性地基對(duì)位移、應(yīng)力、電磁應(yīng)力和磁擾動(dòng)量的影響.Rahmani 等［16］通過考慮兩種類型的孔隙率分布，首次應(yīng)用修正的高階夾層板理論，使用伽遼金法研究了兩種多孔FGM圓形夾層板的振動(dòng)特性.

綜上所述，已有關(guān)于FGM 圓板的文獻(xiàn)表明，對(duì)多孔FGM 圓板的研究則相對(duì)較少，同時(shí)一些研究結(jié)果也很好地驗(yàn)證了微分變換法（Differential Transformation Method，DTM）在求解特征值問題時(shí)的實(shí)用性和有效性［17］.目前關(guān)于熱環(huán)境下彈性地基上多孔FGM 圓板自由振動(dòng)問題的研究在國內(nèi)外還鮮見有文獻(xiàn)報(bào)道，因此本文考慮均勻分布孔隙引入力學(xué)模型，采用DTM 對(duì)不同溫度場(chǎng)中多孔FGM 圓板自由振動(dòng)展開研究，分析梯度指數(shù)、孔隙率、升溫、邊界條件、厚度與半徑比和彈性地基系數(shù)等對(duì)多孔FGM 圓板無量綱固有頻率以及相關(guān)各參數(shù)對(duì)臨界溫升值的影響，為今后多孔FGM圓板的研究提供數(shù)據(jù)支撐.

1 多孔FGM 圓板的物性參數(shù)

1.1 材料物性描述

如圖1（a）所示，選取柱坐標(biāo)(r，θ，z)分別表示徑向、環(huán)向和橫向坐標(biāo)，建立彈性地基上多孔FGM 圓板的力學(xué)模型.多孔FGM 圓板半徑為R，厚度為h，kw為Winkler 彈性剛度系數(shù)，上表面為純陶瓷，下表面為純金屬，含均勻孔隙FGM 各物性參數(shù)P(z，T)（彈性模量E、密度ρ和熱膨脹系數(shù)α）可用改進(jìn)的混合冪律公式表示［18］：

圖1 彈性地基上多孔FGM圓板的幾何模型和橫截面上均勻分布的孔隙Fig.1 Geometric model of porous FGM circular plate on elastic foundation and uniformly distributed pores in cross section

式中：n為梯度指數(shù)；Pc、Pm分別表示無孔隙時(shí)陶瓷和金屬材料的物性參數(shù)；β為孔隙率，均勻分布的孔隙如圖1（b）所示.

上述金屬和陶瓷的某一物性參數(shù)P隨溫度T的變化可采用 Touloukian 非線性函數(shù)［19］統(tǒng)一表述為：

式中：P0、P-1、P1、P2和P3為溫度相關(guān)的系數(shù)，一般由實(shí)驗(yàn)給出，見表4.

1.2 溫度場(chǎng)描述

在不考慮材料導(dǎo)熱率情況下，本文給出一種統(tǒng)一的溫度場(chǎng)表達(dá)式，可以更好地描述溫度場(chǎng)

式中：ΔT=Tc-Tm表示多孔FGM 圓板上下表面溫度差，Tc和Tm分別為上下表面溫度值；η為溫度指數(shù).當(dāng)Tm=T0=300 K時(shí)，表示無應(yīng)力狀態(tài)的參考溫度，此時(shí)當(dāng)η=0 時(shí)為均勻升溫，當(dāng)η=1 時(shí)為線性升溫，當(dāng)η取其他值時(shí)均為非線性升溫.

將式（2）、式（3）代入式（1），即可得到不同孔隙率下多孔FGM圓板的物性參數(shù).

2 控制微分方程的推導(dǎo)與DTM變換

2.1 無量綱控制微分方程

基于小撓度假設(shè)和一階剪切變形理論，多孔FGM圓板的位移分量如下：

式中：ur、uθ和uz分別表示圓板徑向、環(huán)向和橫向位移；φ為法線的轉(zhuǎn)角；w為橫向位移；t為時(shí)間；z0為物理中面.物理中面內(nèi)應(yīng)力分量與應(yīng)變分量為零，其表達(dá)式如下［20］：

圓板考慮軸對(duì)稱，其幾何方程為：

式中：εr、εθ和γrz分別表示徑向應(yīng)變、環(huán)向應(yīng)變和切應(yīng)變.物理方程為

式中：σr，σθ和τrz分別表示徑向應(yīng)力、環(huán)向應(yīng)力和切應(yīng)力；ν為泊松比；T0表示無應(yīng)力狀態(tài)的參考溫度，本文取T0=300 K.

多孔FGM 圓板單位面積的薄膜力和彎矩分別為：

式中：B、C、D為剛度系數(shù)；κs為剪切修正系數(shù)，且取κs=π212；NT和MT分別為熱軸力和熱彎矩，可定義為

應(yīng)用Hamilton 變分原理可推導(dǎo)出熱環(huán)境下彈性地基上多孔FGM 圓板的自由振動(dòng)問題的控制微分方程.Hamilton原理如下：

式中：Π 為Hamilton 作用量；S、U和V分別為系統(tǒng)的動(dòng)能、應(yīng)變能和外力做的功；δ為變分符號(hào)；t1和t2分別表示系統(tǒng)運(yùn)動(dòng)的初始時(shí)刻和終止時(shí)刻.

可得多孔FGM圓板的動(dòng)能的變分為：

慣性系數(shù)為：

其中，

多孔FGM圓板的應(yīng)變能的變分為：

多孔FGM圓板的外力做功的變分為：

式中：kw為Winkler彈性剛度系數(shù).利用上述Hamilton原理，可得到多孔FGM 圓板的熱環(huán)境下彈性地基上多孔FGM圓板自由振動(dòng)的控制微分方程.

將式（8）～式（12）代入式（17）和式（18），并假設(shè)多孔FGM 圓板做簡(jiǎn)諧振動(dòng)，令w(r，t)=wˉ(r)cosωt，φ(r，t)=-φ(r)cosωt，其中ω是多孔FGM圓板的固有頻率，wˉ(r)和-φ(r)是形狀函數(shù).代入式（17）、式（18），可以得到多孔FGM圓板的運(yùn)動(dòng)控制微分方程為：

引入如下無量綱參數(shù)：

式中：λ為厚度與半徑比；Kw為無量綱Winkler 彈性剛度系數(shù)為無量綱熱載荷；Ω為無量綱固有頻率.

綜上可得到熱環(huán)境下彈性地基上多孔FGM 圓板自由振動(dòng)的無量綱控制微分方程

對(duì)于熱環(huán)境下彈性地基上多孔FGM 圓板自由振動(dòng)問題，本文考慮了兩種邊界條件，其形式如下：

ξ=0，圓板中心邊界為：

2.2 DTM變換

對(duì)于熱環(huán)境下彈性地基上多孔FGM 圓板自由振動(dòng)的無量綱控制微分方程組（20）、（21）的特征值問題，這里采用一種半解析方法——微分變換法（DTM）進(jìn)行求解.該方法在求解特征值問題時(shí)具有過程簡(jiǎn)單、編程方便的特點(diǎn)，且對(duì)于特征值無需經(jīng)逆變換過程即可得到［17，21］.按照DTM的求解過程及原理［22］，將其無量綱控制微分方程組轉(zhuǎn)換為如下的迭代代數(shù)方程組可表示為：

邊界條件的DTM變換如下：

ξ=0，圓板中心邊界為：

3 算例分析與數(shù)值討論

應(yīng)用MATLAB 軟件編寫相關(guān)程序，由此可得到通過DTM 求解多孔FGM 圓板自由振動(dòng)時(shí)的無量綱固有頻率.為了驗(yàn)證DTM 的正確性，不考慮材料受溫度依賴，經(jīng)算例1 退化為孔隙率為零時(shí)的FGM 圓板，得到表1、表2 和表3 結(jié)果，并與文獻(xiàn)［23］和［24］中的結(jié)果進(jìn)行比較.考慮材料受溫度依賴，經(jīng)算例2得到圖2～圖9新的結(jié)果.

表1 固支和簡(jiǎn)支邊界條件下一階剪切變形理論和經(jīng)典板理論的無量綱固有頻率對(duì)比（λ=0.01）Tab.1 Comparison of dimensionless natural frequencies between first-order shear deformation theory and classical plate theory under C and S boundary conditions（λ=0.01）

表2 固支和簡(jiǎn)支邊界條件下不同梯度指數(shù)對(duì)FGM 圓板無量綱固有頻率的影響（λ=0.05）Tab.2 Effects of different gradient indices on the dimensionless natural frequencies of FGM circular plates under C and S boundary conditions（λ=0.05）

表3 固支和簡(jiǎn)支邊界條件下不同Winkler系數(shù)對(duì)FGM 圓板無量綱固有頻率的影響Tab.3 Effect of different Winkler parameter on the dimensionless natural frequency of FGM circular plate under C and S boundary conditions

圖2 均勻升溫下梯度指數(shù)對(duì)無量綱前三階固有頻率的影響曲線Fig.2 The influence curve of gradient index on the first three dimensionless natural frequencies under uniform temperature rise

算例1：表1 給出了在固支（C）和簡(jiǎn)支（S）邊界條件下，F(xiàn)GM 圓板取不同梯度指數(shù)時(shí)，基于經(jīng)典板理論（CPT）的前兩階固有頻率［23］與本文基于一階剪切變形理論（FSDT）的DTM 解對(duì)比，發(fā)現(xiàn)DTM 求出的解與文獻(xiàn)解在準(zhǔn)確性和精度上完全吻合.其中取ρc=3 800 kg/m3，ρm=2 700 kg/m3，Ec=380 GPa，Em=70 GPa，ν=0.3，λ=0.01，β=0.

表2 和表3 為固支（C）和簡(jiǎn)支（S）邊界條件下，不同梯度指數(shù)和Winkler 系數(shù)下，基于一階剪切變形理論（FSDT）得到FGM 圓板和彈性地基上FGM圓板振動(dòng)的無量綱前兩階固有頻率，若依據(jù)文獻(xiàn)［24］對(duì)剛度無量綱過程處理方法，得到本文解?與之對(duì)照，發(fā)現(xiàn)其數(shù)值結(jié)果完全吻合.本文在剛度無量綱過程中與文獻(xiàn)［24］不同，引入了物理中面，雖然文中拉-彎耦合剛度為零，但對(duì)結(jié)構(gòu)的耦合影響被帶入彎曲剛度D和慣性剛度I2中，所以得到本文解??與表3 中無量綱前兩階固有頻率略有差別，但也進(jìn)一步證明了DTM 的精度滿足要求，其中取ρc=3 800 kg/m3，ρm=2 702 kg/m3，Ec=380 GPa，Em=70 GPa，ν=0.3，β=0.

算例2：圖2～圖9 中各物性參數(shù)均遵循表4 中的實(shí)驗(yàn)數(shù)據(jù)，由于泊松比的取值對(duì)振動(dòng)頻率幾乎沒有影響，故在此ν取常數(shù).

表4 陶瓷和金屬材料物性參數(shù)的溫度相關(guān)系數(shù)［19］Tab.4 Temperature correlation coefficient of physical parameters of ceramic and metal materials

圖2 和圖3 為多孔FGM 圓板在固支（C）和簡(jiǎn)支（S）邊界下，當(dāng)β=0.1，λ=1，Kw=10，ΔT=100 K時(shí)，不同梯度指數(shù)對(duì)無量綱前三階固有頻率Ω的影響曲線.可以看出在溫度指數(shù)一定時(shí)，隨著梯度指數(shù)n的增加，多孔FGM 圓板的無量綱固有頻率Ω逐漸下降，當(dāng)n在小值范圍取值時(shí)（1 附近），無量綱固有頻率Ω（減小的）波動(dòng)很劇烈；當(dāng)n在較大值范圍取值時(shí)，無量綱固有頻率變化趨于平緩.合理解釋了FGM中陶瓷材料向金屬材料過渡的特性.溫度指數(shù)的變化表示升溫類型的變化，均勻升溫（η=0）和線性升溫（η=1）類型下，梯度指數(shù)對(duì)固有頻率的影響趨勢(shì)基本一致，當(dāng)梯度指數(shù)n一定時(shí)，隨著溫度指數(shù)η增大，無量綱固有頻率在增大，邊界約束越強(qiáng)，無量綱固有頻率越大.

圖3 線性升溫下梯度指數(shù)對(duì)無量綱前三階固有頻率的影響曲線Fig.3 The influence curve of gradient index on the first three dimensionless natural frequencies under linear temperature rise

圖4為多孔FGM 圓板在固支（C）和簡(jiǎn)支（S）邊界下，當(dāng)β=0.1，n=1，Kw=10，ΔT=100 K時(shí)，不同溫度指數(shù)和厚度與半徑比λ對(duì)無量固有綱頻率Ω的影響曲線.可以看出隨著厚度與半徑比λ的增大，多孔FGM圓板的無量綱固有頻率Ω在非線性單調(diào)減小，λ的取值在0.15～0.35 之間，C 邊界下一階固有頻率減小的速度較快，S 邊界下減小得較緩慢，而兩種邊界下二階固有頻率減小的速度均較快.當(dāng)厚度與半徑比λ一定時(shí)，隨著溫度指數(shù)η增大，無量綱固有頻率增大，邊界約束越強(qiáng)，無量綱頻率越大.

圖4 不同厚度與半徑比，不同溫度指數(shù)對(duì)無量綱前兩階固有頻率的影響曲線Fig.4 The influence curve of different thickness-to-diameter ratio and different temperature index on the first two order dimensionless natural frequencies

圖5為多孔FGM 圓板在固支（C）和簡(jiǎn)支（S）邊界下，當(dāng)β=0.1，n=1，Kw=0，λ=0.1和不同溫度指數(shù)時(shí)，溫升值對(duì)無量綱基頻Ω1的影響曲線.由結(jié)構(gòu)穩(wěn)定性理論［25］可知，當(dāng)結(jié)構(gòu)自振基頻為零時(shí)，意味著結(jié)構(gòu)構(gòu)件發(fā)生失穩(wěn)，所加載荷為最危險(xiǎn)載荷即臨界屈曲載荷，此時(shí)對(duì)應(yīng)的溫度即為臨界屈曲溫度.由圖5可見，隨著溫升值的增加，多孔FGM 圓板的無量綱頻率逐漸下降，具體表現(xiàn)在無量綱基頻Ω1均隨溫升值的增大而減小，由于材料受溫度依賴性溫度300 K ≤T≤1 100 K［26］，當(dāng)溫升值取800 K 以內(nèi)時(shí)，η=0 時(shí)，表示均勻升溫下，無量綱基頻Ω1減小至臨界升溫ΔTcr（或臨界屈曲溫度Tcr=T0+ΔTcr），（C 邊界下，ΔTcr=679 K，S 邊界下，ΔTcr=278 K），結(jié)構(gòu)發(fā)生熱屈曲而失穩(wěn)，此時(shí)達(dá)到臨界平衡狀態(tài)，即Ω=0.當(dāng)η=1時(shí)，表示線性升溫，S邊界同樣減小至臨界升溫，此時(shí)ΔTcr=490 K，發(fā)生熱屈曲，C 邊界只有減小的趨勢(shì).同樣，邊界約束越弱，達(dá)到臨界熱屈曲時(shí)的溫升值越小.

圖5 不同升溫類型對(duì)無量綱基頻的影響曲線Fig.5 The influence curve of different heating types on dimensionless fundamental frequency

圖6 和圖7 反映了固支（C）和簡(jiǎn)支（S）邊界和不同升溫類型下，當(dāng)ΔT=100 K，n=1，β=0.1，λ=0.1時(shí)，孔隙率β對(duì)多孔FGM 圓板無量綱前兩階固有頻率Ω的影響曲線.圖中顯示：均勻升溫（η=0）和線性升溫（η=1）類型下，孔隙率β對(duì)無量綱固有頻率Ω的影響趨勢(shì)基本一致，當(dāng)升溫類型和溫升值一定時(shí)，隨著孔隙率β的增大，多孔FGM 圓板的無量綱固有頻率逐漸上升，是因?yàn)榭紫兜拇嬖诓粌H削弱了FGM圓板的整體剛度，同時(shí)也降低了等效質(zhì)量，此時(shí)質(zhì)量的弱化強(qiáng)于剛度弱化.當(dāng)孔隙率β一定時(shí)，隨著溫升值增大，無量綱固有頻率減小，邊界約束越弱，無量綱固有頻率越小.

圖6 均勻升溫下孔隙率對(duì)無量綱前兩階固有頻率的影響曲線Fig.6 The influence curve of porosity on the first two dimensionless natural frequencies under uniform temperature rise

圖7 線性升溫下孔隙率對(duì)無量綱前兩階固有頻率的影響曲線Fig.7 The influence curve of porosity on the first two dimensionless natural frequencies under linear temperature rise

圖8 和圖9 給出了在固支（C）和簡(jiǎn)支（S）邊界條件下，當(dāng)ΔT=100 K，n=1，β=0.1，λ=0.1 時(shí)，無量綱Winkler 彈性剛度系數(shù)Kw對(duì)無量綱固有頻率Ω的影響.結(jié)果表明：當(dāng)無量綱彈性剛度系數(shù)增大時(shí)，無量綱固有頻率增大，且無量綱二階固有頻率幾乎呈線性變化.當(dāng)剛度系數(shù)增大時(shí)，增強(qiáng)了結(jié)構(gòu)的整體剛度，且當(dāng)彈性剛度系數(shù)Kw一定時(shí)，邊界約束越強(qiáng)，無量綱固有頻率越大.

圖8 均勻升溫下無量綱Winkler彈性剛度系數(shù)對(duì)無量綱前兩階固有頻率的影響曲線Fig.8 The influence curve of dimensionless Winkler elastic stiffness coefficient on the first two dimensionless natural frequencies under uniform temperature rise

圖9 線性升溫下無量綱Winkler彈性剛度系數(shù)對(duì)無量綱前兩階固有頻率的影響曲線Fig.9 The influence curve of dimensionless Winkler elastic stiffness coefficient on the first two dimensionless natural frequencies under linear temperature rise

4 結(jié)論

本文基于一階剪切變形理論，考慮修正的混合冪律多孔FGM 圓板模型，在統(tǒng)一溫度場(chǎng)下應(yīng)用Hamilton 原理導(dǎo)出該熱環(huán)境下該結(jié)構(gòu)自由振動(dòng)的控制微分方程，使用微分變換法（DTM）對(duì)自由振動(dòng)和屈曲的運(yùn)動(dòng)控制微分方程及邊界條件進(jìn)行變換，通過編寫MATLAB 程序計(jì)算多孔FGM 圓板的無量綱固有頻率和臨界溫升值，并且與已有文獻(xiàn)對(duì)比以驗(yàn)證結(jié)果的正確性.分析了不同邊界、孔隙率、梯度指數(shù)、溫度指數(shù)、厚度與半徑比、溫升值和Winkler 彈性剛度系數(shù)對(duì)多孔FGM 圓板無量綱固有頻率的影響和相關(guān)參數(shù)各對(duì)臨界溫升值的影響.主要結(jié)論如下：

1）無量綱固有頻率隨著梯度指數(shù)n的增大而減小，隨著梯度指數(shù)n趨向于無窮大，無量綱固有頻率趨于不變，合理地解釋了FGM 中陶瓷材料向金屬材料過渡的特性.

2）隨著厚度與半徑比的增大，無量綱固有頻率減小，隨升溫類型η值的增大，無量綱固有頻率增大.

3）無量綱固有頻率隨著孔隙率β增大而增大，因?yàn)榇嬖诳紫叮档土硕嗫譌GM 圓板的整體剛度和有效質(zhì)量，而對(duì)質(zhì)量的弱化強(qiáng)于剛度弱化的原因?qū)е铝诉@一結(jié)果.

4）隨溫升值的增大，無量綱固有頻率減小，減小直至發(fā)生熱屈曲而失穩(wěn)，使得無量綱基頻變?yōu)榱悖藭r(shí)達(dá)到臨界平衡狀態(tài).

5）無量綱Winkler 彈性剛度系數(shù)Kw越大，無量綱固有頻率Ω越大，且二階固有頻率幾乎呈線性變化.