徐　坤，陳美霞，謝　坤

（華中科技大學(xué) 船舶與海洋工程學(xué)院，武漢 430074）

?

正交各向異性功能梯度材料平板振動(dòng)分析

徐坤，陳美霞，謝坤

（華中科技大學(xué) 船舶與海洋工程學(xué)院，武漢 430074）

基于一階剪切理論，研究四邊簡(jiǎn)支正交各向異性功能梯度材料（FGM）板的自由振動(dòng)和受迫振動(dòng)。假設(shè)剪應(yīng)力沿厚度方向呈拋物線分布，利用剪切應(yīng)變能與剪切余能相等原理，得到正交各向異性功能梯度平板的剪切修正系數(shù)。利用虛位移原理得到功能梯度平板運(yùn)動(dòng)方程，并采用Navier方法對(duì)運(yùn)動(dòng)方程進(jìn)行求解。通過(guò)與有關(guān)文獻(xiàn)及有限元計(jì)算結(jié)果對(duì)比，驗(yàn)證該方法的正確性。在此基礎(chǔ)上，分析厚度方向上由纖維和基體按照不同體積梯度分布的三種（P-，S-，C-FGM）平板的固有振動(dòng)和受激振動(dòng)特性，結(jié)果表明纖維體積分?jǐn)?shù)變化區(qū)間越大，梯度型式及梯度指數(shù)對(duì)其振動(dòng)特性影響越顯著；纖維體積分?jǐn)?shù)關(guān)于平板中面反對(duì)稱(chēng)分布（S-FGM）時(shí)，平板振動(dòng)特性受梯度指數(shù)影響較小。

振動(dòng)與波；功能梯度材料；一階剪切理論；剪切修正系數(shù)

功能梯度材料（FGM）是一種新型復(fù)合材料，其具有空間連續(xù)梯度變化的微觀結(jié)構(gòu)，可引起材料屬性的梯度變化。功能梯度材料板相比層合板消除了材料連接處的宏觀界面，實(shí)現(xiàn)了材料內(nèi)部功能的漸變，具有緩和熱應(yīng)力、避免或降低層間應(yīng)力集中等優(yōu)點(diǎn)。此外，功能梯度材料具有良好的可設(shè)計(jì)性，可通過(guò)合理設(shè)計(jì)來(lái)改變材料的組成及梯度變化形式，實(shí)現(xiàn)不同功能的組合和優(yōu)良的力學(xué)特性。

目前，功能梯度材料板殼的彎曲和振動(dòng)問(wèn)題研究與層合板［1-2］類(lèi)似，主要有基于板理論的解析方法以及有限元方法等。Chi等［3-4］根據(jù)經(jīng)典板理論和Fourier級(jí)數(shù)展開(kāi)，分析了功能梯度平板彎曲問(wèn)題，認(rèn)為功能梯度平板彎曲剛度不同于均質(zhì)板，是關(guān)于材料參數(shù)分布的復(fù)雜組合。曹志遠(yuǎn)［5］和Loy［6］采用經(jīng)典板理論分別分析了矩形平板和簡(jiǎn)支功能梯度材料圓柱殼的固有頻率。Nguyen［7］基于1階剪切變形理論，分析了功能梯度材料平板的彎曲問(wèn)題。Reddy［8］和Akbarzadeh［9］基于3階剪切變形理論，分別得到了功能梯度平板的理論方程，以及通過(guò)Fourier級(jí)數(shù)展開(kāi)和Laplace變換，分析了梯度平板在動(dòng)載荷作用下的響應(yīng)。Ferreira［10］基于1階、3階剪切變形理論及MQ徑向基函數(shù)分析了功能梯度板的自由振動(dòng)。Qian［11］基于高階剪切變形理論及局部彼得洛夫-伽遼金法分析了功能梯度平板的靜力及動(dòng)力變形問(wèn)題。對(duì)于功能梯度材料夾芯板，李華東［12］將位移和橫向剪力用三角級(jí)數(shù)展開(kāi)，研究了四邊簡(jiǎn)支功能梯度夾層板的自由振動(dòng)問(wèn)題。

由于工藝等條件的限制，實(shí)際的功能梯度材料很少是各向同性的，而是大多表現(xiàn)出正交各向異性［13］。Ramirez［14］采用離散分層理論結(jié)合Ritz法分析了各向異性功能梯度板的彎曲問(wèn)題。Batra［15］基于1階剪切變性理論及有限元法，分析了纖維鋪角對(duì)各向異性功能梯度平板固有頻率的影響。Zhang［16］基于3階剪切變形理論，分析了正交各向異性功能梯度平板的非線性振動(dòng)問(wèn)題。陳偉球［17］通過(guò)引入位移函數(shù)和應(yīng)力函數(shù)，結(jié)合分層近似理論分析了橫觀各向同性功能梯度材料板的自由振動(dòng)。

文獻(xiàn)［3-11］的研究對(duì)象主要為各向同性材料功能梯度板，而且大都只假設(shè)材料的楊氏模量在厚度方向上梯度變化，未考慮梯度材料中密度等其他材料參數(shù)變化。此外，采用1階剪切變形理論求解功能梯度平板的振動(dòng)問(wèn)題同高階方法［8-16］相比，不僅具有足夠的精度［17］，且公式中變量的物理意義明確、步驟更加簡(jiǎn)便、求解效率更高。

文中基于1階剪切變形理論，假設(shè)了剪力沿厚度方向上的分布形式，導(dǎo)出了正交各項(xiàng)異性功能梯度板的剪切修正系數(shù)計(jì)算式及四邊簡(jiǎn)支功能梯度平板自由振動(dòng)及受迫振動(dòng)的理論公式，并分析了不同組分體積梯度分布的三種（P-，S-，C-FGM）正交各向異性功能梯度平板的自由振動(dòng)和受迫振動(dòng)。

1　基本公式

考察受到橫向載荷作用且材料參數(shù)沿厚度方向變化的正交各向異性梯度材料矩形板。圖1是一種由纖維和基體在厚度方向上按照不同比例混合得到的正交各項(xiàng)異性梯度平板。由于厚度方向上組分體積分?jǐn)?shù)梯度變化，平板沿厚度方向材料參數(shù)呈現(xiàn)梯度變化。此外，增強(qiáng)纖維的方向性使得平板呈正交各向異性。假設(shè)平板等厚度，平板材料彈性主軸方向與圖1中平板分析坐標(biāo)系一致，板長(zhǎng)為a，寬為b，厚度為h。

根據(jù)平板1階剪切變形的假設(shè)［18］，平板內(nèi)任意一點(diǎn)的位移為

圖1　纖維與基體梯度混合的功能梯度平板

沿厚度方向的剪應(yīng)變?yōu)?/p>

在1階變形假設(shè)下，板的應(yīng)力應(yīng)變關(guān)系滿(mǎn)足

式中Q11、Q12、Q22、Q44、Q55、Q66分別為板厚方向上的剛度系數(shù)，計(jì)算公式［18］為

G12、G13、G23為材料的剪切模量。截面合力N=（Nx，Ny，Nxy）T及截面彎矩M=（Mx，My，Mxy）T為

其中

矩陣元素Aij，Bij是關(guān)于Qij及厚度方向坐標(biāo)z的積分，表達(dá)式分別為

由于梯度材料平板的材料參數(shù)沿板厚方向呈梯度變化，沿厚度方向的剪切變形不能忽略。采用1階剪切變形假設(shè)分析功能梯度材料平板的動(dòng)力學(xué)問(wèn)題時(shí)，由于廣義胡克定律和橫向剪切剛度為有限值，該假設(shè)下橫向剪應(yīng)力沿厚度方向均勻分布，這與平板上下表面剪應(yīng)力為零的邊界條件相矛盾，為了解決這個(gè)矛盾，需要對(duì)剪切系數(shù)進(jìn)行修正［18］。這里對(duì)文獻(xiàn)［19］中的剪力分布假設(shè)進(jìn)行擴(kuò)展，假設(shè)梯度材料平板截面厚度方向上的應(yīng)力分布與均勻材料梁的剪應(yīng)力分布類(lèi)似，即假設(shè)沿厚度h方向的平板截面的剪應(yīng)力有式（8）中的拋物線型的分布形式

其中Qx和Qy分別為平板垂直于x和y單位長(zhǎng)度截面上的剪力的合力。由于功能梯度平板材料厚度方向上的梯度變化，實(shí)際剪應(yīng)力一般不再呈拋物線型式分布，但工程上采用該拋物線分布假設(shè)分析梯度平板得到的結(jié)果有足夠精度。根據(jù)剪切余能的定義，單位面積功能梯度板的剪切余應(yīng)變能為

將式（8）代入式（9），得到應(yīng)變能的表達(dá)式

假設(shè)C44、C55為修正后的剪切剛度，用C44、C55表示功能梯度平板的剪切余應(yīng)變能密度為

按照剪切余應(yīng)變能密度相等得到

修正前單位長(zhǎng)度截面剪力表達(dá)式為

修正后的截面剪力的表達(dá)式為

由剪切修正系數(shù)的定義，得到剪切修正系數(shù)為

由式（15）可知，功能梯度材料平板的剪切修正系數(shù)是關(guān)于假設(shè)的剪應(yīng)力分布函數(shù)以及功能梯度平板厚度方向的材料剪切模量的分布函數(shù)共同作用的表達(dá)式，在已知材料剪切模量分布的情況下，通過(guò)式（15）得到功能梯度平板修正后的剪切剛度。

根據(jù)虛位移原理，類(lèi)比層合板的振動(dòng)平衡方程［16］，得到功能梯度平板的動(dòng)力學(xué)方程為

對(duì)于四邊簡(jiǎn)支功能梯度板，平衡方程可采用雙三角級(jí)數(shù)方法求解［18］，假設(shè)功能梯度平板的位移滿(mǎn)足式（18）中的雙三角級(jí)數(shù)表達(dá)式

其中umn、νmn、φxmn、φymn和wmn為廣義位移幅值，式中的m和n分別為平板振動(dòng)時(shí)在x方向及y方向的半波數(shù)。同樣，將橫向載荷展開(kāi)成

截取位移函數(shù)表達(dá)式（18）前N項(xiàng)，代入平衡方程（16）中，并在平衡方程等式兩邊分別乘上對(duì)應(yīng)的三角函數(shù)表達(dá)式，利用三角函數(shù)正交性解耦，得到關(guān)于廣義位移幅值的線性方程組

其中m，n，k，l=1，2…N，等式左邊是關(guān)于廣義位移為未知變量的線性組合，等式右邊為橫向載荷廣義力。求解得到功能梯度平板各個(gè)廣義位移幅值，得到平板在橫向力作用下的響應(yīng)。

對(duì)于功能梯度平板的固有頻率，令qmn=0，由任意階振型（m，n）下廣義位移幅值不全為零，由式（21）系數(shù)矩陣的行列式為零，可得對(duì)應(yīng)振型下的功能梯度板固有頻率值。

其中系數(shù)矩陣中的各元素即為平衡方程式（16）中對(duì)應(yīng)振型廣義位移幅值的系數(shù)。

2　數(shù)值分析

2.1方法驗(yàn)證

首先將本文方法的結(jié)果同文獻(xiàn)［20］各向同性功能梯度材料平板的固有頻率進(jìn)行對(duì)比，如表1所示（無(wú)量綱化及p見(jiàn)文獻(xiàn)［20］）。本文方法結(jié)果同文獻(xiàn)結(jié)果差別較小，表明本方法能夠較準(zhǔn)確退化分析各向同性功能梯度材料的固有振動(dòng)。

對(duì)于正交各向異性功能梯度板，考察圖1所示的纖維和基體沿平板厚度方向按照一定梯度變化的比例混合而成的簡(jiǎn)支平板，板長(zhǎng)a=1.2 m，板寬b= 0.8 m，厚0.02 m。在平板厚度方向材料參數(shù)與纖維和基體的比例存在如式（22）的關(guān)系［19］。

其中E1為梯度板縱向楊氏模量，E2為橫向楊氏模量，ν12為主泊松比，G12為剪切模量。假設(shè)G23=G12，Gf和Gm分別為纖維和基體的剪切模量，纖維和基體均為各向同性材料，取纖維楊氏模量Ef=72 000 Mpa，基體楊氏模量Em=6 900 Mpa，纖維泊松比νf=0.22，基體的泊松比νm=0.4，纖維的密度ρf=2 500 kg/m3，基體的密度ρm=1 200 kg/m3，且剪切模量滿(mǎn)足式（22）中的關(guān)系。

假設(shè)基體材料體積分?jǐn)?shù)Vm沿厚度方向呈梯度變化，且存在式（24）中的P、S、C三種型式的體積分?jǐn)?shù)分布函數(shù)

S分布：

三種分布型式的表達(dá)式（24）中，V1、V2分別是功能梯度平板中基體體積分?jǐn)?shù)Vm的最小值和最大值。特別地，對(duì)于P、C兩種分布，V1、V2分別為功能梯度平板下表面和上表面處的基體體積分?jǐn)?shù)。式（27）中的指數(shù)γ選取不同的值，基體體積分?jǐn)?shù)Vm的曲線也會(huì)呈現(xiàn)不同的變化趨勢(shì)，稱(chēng)γ為梯度指數(shù)。圖2是V1=0.25，V2=0.75時(shí)，P、S、C三種分布型式功能梯度板厚度方向Vm隨γ變化曲線。

如圖2（a），P型分布下功能梯度平板厚度方向上基體材料的組分呈現(xiàn)冪律變化，梯度指數(shù)較大或較小時(shí)，基體體積分?jǐn)?shù)會(huì)在平板一側(cè)表面出現(xiàn)較大梯度的變化。S型分布情況下，基體體積分?jǐn)?shù)曲線關(guān)于平板中面呈類(lèi)似“S”反對(duì)稱(chēng)分布，C型分布情況基體材料體積分?jǐn)?shù)關(guān)于平板中面對(duì)稱(chēng)。

分別計(jì)算P分布形式下不同γ值下梯度平板的固有頻率，并將固有頻率同有限元方法結(jié)果進(jìn)行對(duì)比。鑒于目前尚沒(méi)有功能梯度材料有限元分析的商用軟件，采用層合平板來(lái)近似等效功能梯度平板，即在有限元建模中，將梯度平板沿厚度方向上離散分成多層，每一層視為均勻材料，各相鄰層之間實(shí)現(xiàn)材料參數(shù)階梯變化，來(lái)近似模擬材料參數(shù)沿厚度方向上的連續(xù)的梯度變化。文中分層數(shù)為20時(shí)有限元結(jié)果已收斂，故有限元模型均采用20層等效。通過(guò)有限元方法同本文方法的對(duì)比驗(yàn)證本文方法的正確性。表2是本文方法與有限元方法計(jì)算功能梯度板不同振型對(duì)應(yīng)的固有頻率結(jié)果對(duì)比。

表1　不同厚寬比及p值下Al/ZrO2功能梯度方板基頻對(duì)比

為驗(yàn)證梯度平板動(dòng)響應(yīng)分析正確性，計(jì)算梯度平板在受到點(diǎn)力作用下的頻響。激勵(lì)力作用于距離平板長(zhǎng)邊和短邊各1/4處，如圖1，方向垂直于板平面，大小為1 N。圖3（a）和圖3（b）分別是本文方法與有限元方法計(jì)算P分布下激勵(lì)點(diǎn)速度頻響和平板均方振速的對(duì)比曲線。

模態(tài)分析和諧響應(yīng)分析結(jié)果表明1階剪切方法與有限元方法吻合得較好，本文方法適用于梯度平板的動(dòng)力學(xué)分析。

2.2功能梯度平板固有振動(dòng)分析

對(duì)2.1中的模型，當(dāng)V1=0.25，V2=0.75時(shí)P、S、C梯度型式平板的基頻隨γ變化如圖4。

由圖4可知，P、C分布梯度平板的基頻隨著梯度指數(shù)γ的增加而增大。這是由于γ越大，富含纖維的材料沿厚度方向的分布范圍越大。纖維的楊氏模量比基體的楊氏模量相比要大很多，兩者按照比例混合后，纖維的體積分?jǐn)?shù)主要決定了梯度平板的剛度。纖維豐度高的材料沿厚度方向的分布范圍越大，平板的剛度也就越大，相應(yīng)振型的固有頻率也就越高。對(duì)于S分布梯度板，隨著梯度指數(shù)變化，關(guān)于中面對(duì)稱(chēng)的位置纖維體積分?jǐn)?shù)總是呈現(xiàn)相反的變化，板的質(zhì)量及剛度變化較小，因此γ對(duì)其固有振動(dòng)特性影響不大。

表2　功能梯度板各階振型固有頻率同有限元結(jié)果對(duì)比

圖2　不同梯度指數(shù)γ下基體體積分?jǐn)?shù)的分布曲線

圖3　本文方法與有限元結(jié)果對(duì)比

圖4　功能梯度板基頻隨梯度指數(shù)變化曲線

圖5（a）—圖5（c）所示的是V1和V2取不同組合時(shí)P、S、C梯度型式平板基頻隨γ的變化曲線。

由圖5（a）可知，P分布功能梯度平板的基頻總是隨著梯度指數(shù)的增加而增大，且（V2-V1）的差值越大，這種增加的趨勢(shì)越明顯。由圖5（b）可知，對(duì)于S型分布功能梯度平板，基頻隨著梯度指數(shù)的增加先增大，后減小。由圖5（c）可知，對(duì)于C型分布梯度平板，（V2-V1）對(duì)基頻的影響同P分布具有類(lèi)似的規(guī)律。P、S分布中，（V2-V1）決定了平板上下兩個(gè)表面基體材料體積分?jǐn)?shù)的差異，這表明功能梯度材料平板兩側(cè)材料差異越小，梯度指數(shù)對(duì)基頻的影響越小。C分布中，（V2-V1）決定了功能梯度材料板厚度中心位置與上下表面的基體材料體積的差異，這種差異越小，γ對(duì)基頻的影響越小。

2.3功能梯度平板受迫振動(dòng)分析

對(duì)2.1中的模型，分析平板在點(diǎn)力作用下的均方振速，圖6（a）—圖6（c）是不同γ值下平板均方振速。

由圖6（a）可知，P型分布梯度平板的均方振速峰值頻率隨著γ值的增加而增大，這與平板固有頻率隨γ的變化規(guī)律一致。此外，大部分頻率下均方振速隨著γ值的增大而減小，這與平板中纖維體積分?jǐn)?shù)增加，剛度有所增加有關(guān)。由圖6（b）可知，S分布梯度板的均方振速變化隨γ變化并不明顯，這是由于平板的材料參數(shù)總是關(guān)于平板的中面呈現(xiàn)相反規(guī)律的變化，對(duì)平板的振動(dòng)特性影響起到了相互削減的作用。由圖6（c）可知，當(dāng)γ=1與γ=1/5時(shí)梯度平板的均方振速差異較小，若P繼續(xù)增加，當(dāng)γ=5時(shí)，平板的均方振速峰值頻率增大，且大部分頻率下均方振速水平有所下降。這是由于γ＜1時(shí)，材料參數(shù)隨γ的變化并不明顯，均方振速隨γ變化不大。γ從1到5，厚度方向上纖維的比例顯著增加，使得均方振速水平有所降低。

圖5　平板基頻隨體積分?jǐn)?shù)變化曲線

圖6　功能梯度板均方振速隨梯度指數(shù)變化曲線

3　結(jié)語(yǔ)

基于1階剪切理論得到簡(jiǎn)支正交各向異性梯度平板的自由振動(dòng)及受橫向載荷激勵(lì)下的諧響應(yīng)分析方法，通過(guò)對(duì)特定體積梯度混合得到的正交各向異性梯度復(fù)合材料平板的振動(dòng)特性的數(shù)值分析，可以得到以下結(jié)論：

（1）采用1階剪切方法分析正交各向異性功能梯度平板的振動(dòng)特性，通過(guò)與有限元方法的對(duì)比，表明本文方法具有較高的精度，利用該方法能夠較準(zhǔn)確地分析功能梯度材料平板的固有振動(dòng)和受迫振動(dòng)問(wèn)題；

（2）功能梯度材料平板的固有振動(dòng)特性受到組分材料體積分?jǐn)?shù)變化區(qū)間、梯度型式、梯度指數(shù)等因素的綜合影響。平板中纖維的體積分?jǐn)?shù)的變化區(qū)間越大，梯度型式及梯度指數(shù)對(duì)其振動(dòng)特性影響作用越顯著。當(dāng)纖維材料的體積分?jǐn)?shù)變化區(qū)間一定時(shí)，類(lèi)似S型反對(duì)稱(chēng)分布型式的平板受到梯度指數(shù)的影響較小，而冪律（P-）及對(duì)稱(chēng)分布型式（C-）功能梯度平板的振動(dòng)特性受梯度指數(shù)的影響較大；

（3）功能梯度材料平板的受迫振動(dòng)特性受到組分體積分?jǐn)?shù)、梯度型式以及梯度指數(shù)的影響，影響規(guī)律同固有振動(dòng)特性規(guī)律一致，表現(xiàn)為平板厚度上纖維分布范圍越大，平板剛度變大，均方振速有減小的趨勢(shì)。

［1］胡昊灝，商德江.復(fù)合層合矩形板水下聲輻射解析計(jì)算［J］.噪聲與振動(dòng)控制，2014，01：205-208+217.

［2］趙飛，吳錦武，趙龍勝.采用分層理論計(jì)算層合板的固有頻率和振型［J］.噪聲與振動(dòng)控制，2014，02：34-40.

［3］ CHI S H，CHUNG Y L.Mechanical behavior of functionally graded material plates under transverse load-Part I:Analysis［J］.International Journal of Solids and Structures，2006，43（13）:3657-3674.

［4］ CHI S H，CHUNG Y L.Mechanical behavior of functionally graded material plates under transverse loadpart II:numerical results［J］.International Journal of Solids and Structures，2006，43（13）:3675-3691.

［5］曹志遠(yuǎn).不同邊界條件功能梯度矩形板固有頻率解的一般表達(dá)式［J］.復(fù)合材料學(xué)報(bào)，2005，22（5）：172-177.

［6］ LOY C T，LAM K Y，REDDY J N.Vibration of functionally graded cylindrical shells［J］.International Journal of Mechanical Sciences，1999，41（3）:309-324.

［7］NGUYEN T K，SAB K，BONNET G.First-order shear deformation plate models for functionally graded materials［J］.Composite Structures，2008，83（1）:25-36.

［8］REDDY J N.Analysis of functionally graded plates［J］. InternationalJournalforNumericalMethodsin Engineering，2000，47（1-3）:663-684.

［9］AKBARZADEH A H，ZAD S K H，ESLAMI M R，et al. Mechanical behaviour of functionally graded plates under staticanddynamicloading［C］.Proceedingsofthe Institution of Mechanical Engineers，Part C:Journal of Mechanical Engineering Science，2011，225（2）:326-333.

［10］FERREIRA A J M，BATRA R C，ROQUE C M C，et al. Natural frequencies of functionally graded plates by a meshless method［J］.Composite Structures，2006，75（8）: 593-600.

［11］QIAN L F，BATRA R C，CHEN L M.Static and dynamic deformations of thick functionally graded elastic plates by using higher-order shear and normal deformable plate theory and meshless local Petrov-Galerkin method［J］. Composites Part B Engineering，2004，35:685-697.

［12］李華東，朱錫，梅志遠(yuǎn)，等.正交各向異性功能梯度夾層板的自由振動(dòng)分析［J］.哈爾濱工程大學(xué)學(xué)報(bào)，2014，（3）.

［13］馬力，吳林志.各向異性功能梯度材料中裂紋對(duì)彈性波的散射［C］.中國(guó)力學(xué)學(xué)會(huì)學(xué)術(shù)大會(huì) 2005論文摘要集（下），2005.

［14］RAMIREZ F，HEYLIGER P R，PAN E.Static analysis of functionally graded elastic anisotropic plates using a discretelayerapproach［J］.CompositesPartB Engineering，2006，37（1）:10-20.

［15］BATRA R C，JIN J.Natural frequencies of a functionally graded anisotropic rectangular plate［J］.Journal of Sound &Vibration，2005，282（1-2）:509-516.

［16］ZHANG W，YANG J，HAO Y.Chaotic vibrations of an orthotropic FGM rectangular plate based on third-order shear deformation theory［J］.Nonlinear Dynamics，2010，59（4）:619-660.

［17］陳偉球，葉貴如.橫觀各向同性功能梯度材料矩形板的自由振動(dòng)［J］.振動(dòng)工程學(xué)報(bào)，2001，14（3）：263-267.

［18］Reddy J N.Mechanics of laminated composite plates and shells:theory and analysis［M］.CRC press，2004.

［19］王震鳴，劉國(guó)璽，呂明身.各向異性多層扁殼的大撓度方程［J］.應(yīng)用數(shù)學(xué)和力學(xué)，1982，3（1）：49-65.

［20］FERREIRA A J M，BATRA R C，ROQUE C M C，et al. Natural frequencies of functionally graded plates by a meshless method［J］.Composite Structures，2006，75（8）: 593-600.

VibrationAnalysis of Orthotropic Functionally Graded Plates

XUKun，CHEN Mei-xia，XIEKun

（School of NavalArchitecture and Ocean Engineering，Huazhong University of Science and Technology，Wuhan 430074，China）

Based on the first order shear deformation theory（FSDT），free and forced vibrations of orthotropic functionally graded material（FGM）plates with simply supported boundary conditions are analyzed.Utilizing the principle of shear strain energy and assuming that the shear stress distribution has the parabola form across the plate thickness，the shear correction coefficient of the orthotropic FGM plate is obtained.The dynamic equations of the plate are deduced by the principle of virtual displacement and solved by the Navier’s double Fourier series method.The present method is verified by finite element solutions.On this basis，the free and forced vibration characteristics of the plates with three different functional gradient types indicated by P-，S-and C-FGM respectively，are obtained and analyzed.These different functional gradients across the thickness of the plate are realized by different fibre-to-matrix volume ratios.The results show that the effects of gradient type and gradient exponent on FGM plates’vibration characteristics are more obvious when the volume fraction of the fibre varies in a larger range.The gradient exponent has little effect when the volume fraction of the fibre is antisymmetrically distributed（indicated by S-FGM）across the plate thickness.

vibration and wave；functionally graded material；first order shear deformation theory（FSDT）；shear correction coefficient

O327

ADOI編碼：10.3969/j.issn.1006-1335.2016.04.003

1006-1355（2016）04-0014-07

2016-01-02

國(guó)家自然科學(xué)基金資助項(xiàng)目（51179071）

徐坤（1989-），男，湖北省隨州市人，碩士生，主要研究方向?yàn)檎駝?dòng)與噪聲控制。

陳美霞，女，碩士生導(dǎo)師。E-mail:chenmx26@hust.edu.cn