黃小林 劉思奇 肖薇薇 吳 迪

(桂林電子科技大學(xué)建筑與交通工程學(xué)院，廣西桂林541004)

功能梯度材料(functionally graded materials,FGM)的板、殼構(gòu)件已經(jīng)廣泛應(yīng)用于航空航天、交通、機械等各種工程領(lǐng)域，研究此類構(gòu)件的動力特性對其設(shè)計有重要的理論意義和工程價值。

迄今為止，F(xiàn)GM圓錐殼的振動問題已引起了一些研究者的關(guān)注。Zhang等[1]基于Love的第一近似理論，用廣義微分求積法分析了不同邊界條件下FGM圓錐殼的自由振動，結(jié)果發(fā)現(xiàn)，除了材料組分指數(shù)外，邊界條件、幾何形狀等因素對FGM圓錐殼的基頻也有顯著影響。Nezhadi等[2]基于一階剪切變形理論，研究了FGM圓錐殼在四種脈沖負載類型的沖擊載荷作用下的自由振動和強迫振動，結(jié)果表明，對于不同錐角，隨著周向波數(shù)的增加，F(xiàn)GM圓錐殼的固有頻率均先減小后增大，且金屬材料在階躍脈沖作用下產(chǎn)生的動力響應(yīng)最大。

對圓錐殼與彈性地基的相互作用，Sofiyev等[3]用Pasternak雙參數(shù)地基模型對彈性地基上FGM圓錐殼的自由振動進行了分析，討論了彈性地基、圓錐殼的大小頭半徑比、材料體積分數(shù)等因素對量綱歸一化基頻的影響。結(jié)果顯示，Winkler-Pasternak地基對量綱歸一化頻率參數(shù)的影響十分顯著。Chan等[4]基于Pasternak雙參數(shù)地基模型研究了熱環(huán)境中彈性介質(zhì)包圍的偏心加筋功能梯度截錐殼的自由振動和非線性動力響應(yīng)，結(jié)果發(fā)現(xiàn)，隨著彈性地基系數(shù)的增大，ES-FGM圓錐殼的振幅值減小，而且Pasternak參數(shù)的影響比Winkler參數(shù)的大。通過對黏彈性地基上FGM板的振動分析，黃小林等[5]同樣發(fā)現(xiàn)Pasternak雙參數(shù)中的剪切和壓縮剛度對其振動頻率和動力響應(yīng)有顯著影響。

在實際工程中，由于材料制備技術(shù)的不足，F(xiàn)GM在生產(chǎn)過程中常出現(xiàn)一些孔隙，而這些孔隙對構(gòu)件力學(xué)特性的影響不可避免[6]。Cong等[7]通過研究彈性地基上多孔FGM板在機械、熱和熱機械載荷作用下的屈曲和后屈曲行為，發(fā)現(xiàn)孔隙體積分數(shù)對多孔FGM板的屈曲和后屈曲的影響顯著。Huang等[8]通過對非線性彈性地基上多孔sigmoid功能梯度板的非線性自由振動和強迫振動的研究，發(fā)現(xiàn)孔隙體積分數(shù)、線性Winkler參數(shù)和Pasternak地基參數(shù)對板的最大中心瞬態(tài)撓度有顯著影響。Farshid等[9]在彈性地基上雙向功能梯度多孔截錐殼在不同邊界條件下的動力屈曲問題中指出，均勻孔隙系數(shù)和不均勻孔隙系數(shù)對動力失穩(wěn)的影響取決于材料梯度指數(shù)的大小。他們的研究結(jié)果還表明，彈性地基使雙向FGM圓錐殼具有最佳動力性能。文獻[10-12]分別用不同的方法研究含孔隙的FGM板、殼結(jié)構(gòu)的動力響應(yīng)和動力穩(wěn)定特性，同樣發(fā)現(xiàn)了孔隙的影響不可忽略。

目前關(guān)于含孔隙構(gòu)件的研究都集中于梁和板，對彈性地基中含孔隙的FGM圓錐殼的振動鮮有研究。因此，本文假設(shè)FGM圓錐殼由金屬和陶瓷兩種材料構(gòu)成，考慮孔隙體積不可忽略，根據(jù)金屬和陶瓷材料的質(zhì)量組份計算相應(yīng)的體積組份，用改進的混合律模型分別計算均勻和非均勻孔隙分布的FGM物性參數(shù)，探討孔隙、彈性地基參數(shù)、半錐角等因素對FGM圓錐殼自由振動和動力響應(yīng)的影響，以期為進一步了解FGM圓錐殼動力特性，推廣FGM殼體在相關(guān)工程的應(yīng)用提供理論基礎(chǔ)。

1 基本方程及求解方法

1.1 物性參數(shù)

如圖1所示，假設(shè)受橫向動載荷q(s,θ,z)作用下的FGM截頂圓錐殼置于彈性地基中，其厚度為h，母線長度為L，半錐角為γ，R1和R2分別為小端和大端的平均半徑，Kw和Kp分別為彈性地基的壓縮剛度和剪切剛度參數(shù)。采用右手正交曲線坐標系(s,θ,z)，其中s為母線方向以錐頂為原點的坐標，θ為圓周方向坐標，z為厚度方向坐標(?h/2≤z≤h/2)。圖2給出了兩種孔隙分布類型的FGM圓錐殼截面圖，其中，圖2(a)為均勻分布孔隙，圖2(b)為非均勻分布孔隙。

圖1 彈性地基中含孔隙的FGM圓錐殼

圖2 孔隙分布類型

設(shè)Wc和Wm分別表示殼中陶瓷和金屬材料的質(zhì)量含量，有

兩種組分材料和孔隙的體積含量應(yīng)滿足

式中，，和α?分別表示整個圓錐殼的陶瓷、金屬和孔隙的體積含量。根據(jù)式(1)和式(2)，可計算陶瓷和金屬材料體積含量和分別為

其中，ρc和ρm分別為陶瓷和金屬材料的密度。

考慮孔隙沿圓錐殼厚度方向有均勻和非均勻分布兩種模式，本文非均勻分布類型為常見的沿殼體厚度方向中間孔隙多，頂部和底部孔隙少的情況，因此，其沿殼厚度方向的孔隙率α均勻分布時為

非均勻分布時為

設(shè)陶瓷材料體積組分Vc沿厚度方向冪律分布，則

式中，N為材料的組分指數(shù)(0≤N≤∞)。根據(jù)改進的混合律模型，可計算FGM圓錐殼的彈性模量、密度和泊松比等有效物性參數(shù)為

式中，下標c和m分別表示陶瓷和金屬相應(yīng)的物性參數(shù)。由式(5)和式(6)，可得孔隙均勻分布的FGM圓錐殼的物性參數(shù)為

孔隙非均勻分布的FGM圓錐殼的物性參數(shù)為

1.2 控制方程及求解

假定彈性地基在變形過程中始終與圓錐殼保持接觸，忽略基礎(chǔ)中阻尼和慣性力的影響，地基對圓錐殼的反力為

式中

其中，W為圓錐殼中面的撓度。引入應(yīng)力函數(shù)F，F(xiàn)與殼面內(nèi)力關(guān)系為

為了簡化計算過程，引入下列變換

基于經(jīng)典薄殼理論，依據(jù)哈密爾頓變分原理可導(dǎo)出彈性地基中此圓錐殼的運動方程為

假設(shè)殼的兩端為簡支邊界條件，滿足邊界條件的運動方程(12)和方程(13)的解可表示為[13]

式中w(t)為待定函數(shù)，m1=mπ/x0，x0=ln(S2/S1)，n1=n/sinγ，m，n分別為沿母線的半波數(shù)、沿平行圓的全波數(shù)。將式(14)代入式(13)，由諧波平衡法可求得應(yīng)力函數(shù)為

式中，f(t)為待定函數(shù)，系數(shù)Ki定義見文獻[4]。將式(14)和式(15)代入式(12)，并在式(12)兩邊同乘exsin(m1x)sin(n1φ)，然后在殼內(nèi)進行Galerkin積分，將系數(shù)f(t)用w(t)表示再代入式(12)，得到關(guān)于w(t)的二階線性常微分方程組

式中，M和K分別為FGM圓錐殼的廣義質(zhì)量和剛度系數(shù)，且

當(dāng)q(s,φ,t)=0時，式(16)即為彈性地基上含孔隙的FGM圓錐殼的自由振動方程，令w(t)=ae?wt，代入式(16)，可求得自由振動頻率。當(dāng)q(s,φ,t)=0時，用Newmark數(shù)值積分法可求得橫向動載荷作用下的動力響應(yīng)。

2 算例比較與參數(shù)分析

以下計算實例中，F(xiàn)GM圓錐殼由陶瓷(Si3N4)和金屬(Ni)兩種材料組成，Ni物性參數(shù)為Em=205.098 GPa，νm=0.31，ρm=8900 kg/m3，Si3N4的物性參數(shù)為Ec=322.27 GPa，νc=0.24，ρc=2370 kg/m3。R1/h=100，L=2R1，量綱歸一化頻率參數(shù)

2.1 算例比較

表1將本文計算的FGM圓錐殼的最低量綱歸一化頻率參數(shù)與文獻[3]的結(jié)果進行了對比。取半錐角γ=π/6。由表1可知，本文的計算結(jié)果與文獻[3]接近，最大相對誤差約為3.6%。

表1 彈性地基上兩端簡支FGM截錐殼最低量綱歸一化頻率f的比較

2.2 參數(shù)分析

下面分析討論孔隙、材料組分指數(shù)、半錐角等因素對兩端簡支邊界條件下FGM截頂圓錐殼振動頻率和動力響應(yīng)的影響。取模態(tài)(m,n)為(1,6)，h=0.01 m，R1=1.0 m，L=2.0 m。

表2計算了不同孔隙率、材料組分指數(shù)、彈性地基參數(shù)對孔隙為均勻和非均勻分布兩種類型的FGM截錐殼自由振動的影響。取半錐角γ=π/6，陶瓷材料質(zhì)量組分Wc=0.5。從表中可以看出，當(dāng)無彈性地基時，含均勻分布和非均勻分布孔隙的頻率均隨孔隙率的增大而降低，且孔隙均勻分布的頻率均比孔隙非均勻分布的要高。這是由于孔隙率的增大一方面削弱了殼的整體剛度引起量綱歸一化頻率降低，另一方面減小了殼的質(zhì)量而引起量綱歸一化頻率提高，如果孔隙對整體剛度的影響比質(zhì)量的影響大，則量綱歸一化頻率降低，且孔隙為均勻分布對FGM圓錐殼整體剛度的影響比非均勻分布要大。在彈性地基作用下，含均勻分布孔隙的頻率隨孔隙率的增大而提高，含非均勻分布孔隙的頻率在地基參數(shù)(Kw,Kp)=(10 MN/m3，0)條件下，在N∈(0,1]區(qū)間內(nèi)隨著孔隙率的增大而降低，在(1,+∞)區(qū)間內(nèi)隨著孔隙率的增大而提高。在地基剪切參數(shù)Kp=10 MN/m3條件下，在N∈(0,0.2]區(qū)間內(nèi)隨孔隙率的增大而降低，在N∈(0.2,+∞)區(qū)間內(nèi)隨孔隙率的增大而提高。正如前所述，當(dāng)彈性地基壓縮參數(shù)和剪切參數(shù)為一定值時，兩種孔隙的頻率隨著孔隙率的增大而提高，并且孔隙均勻分布的頻率均比孔隙非均勻分布的頻率要高。這是由于彈性地基強化了殼的整體剛度，若孔隙對殼質(zhì)量的影響效果比對整體剛度的影響顯著，則引起量綱歸一化頻率提高。從表中還可以看出，量綱歸一化頻率隨著地基壓縮參數(shù)Kw、剪切參數(shù)Kp的增大而提高，隨著材料組份指數(shù)的增大而降低，地基剪切剛度Kp對頻率的影響比壓縮剛度Kw顯著，可見，彈性地基參數(shù)的增大強化了圓錐殼的整體剛度。

表2 彈性地基上含孔隙FGM截錐殼的量綱歸一化頻率f的比較

圖3～圖7討論了孔隙率、陶瓷材料質(zhì)量組分、半錐角等因素對孔隙均勻與非均勻分布的功能梯度材料圓錐殼在s=(S1+S2)/2中面圓周處動力響應(yīng)的影響。均布突加載荷q=3000 Pa。

圖3顯示了孔隙率α?=0.05，0.1，0.2時孔隙均勻和非均勻分布的動撓度W隨時間t變化的曲線，由圖3可看出，兩種孔隙分布類型的動撓度峰值均隨孔隙率的增大而增大，其中孔隙均勻分布的動撓度隨孔隙率變化更為敏感，并且在相同孔隙率的條件下，孔隙均勻與非均勻分布的振幅值最大相差約35%。

圖4顯示了不同陶瓷材料質(zhì)量組分對孔隙均勻和非均勻分布的FGM圓錐殼動撓度的影響，從圖4可看出，動撓度峰值隨陶瓷材料質(zhì)量組分的增大而減小。

圖4 陶瓷材料質(zhì)量組分對圓錐殼動力響應(yīng)的影響(γ=π/4,N=2,α?=0.1,Kw=0,KP=0)

圖5給出了三種不同的半錐角對孔隙均勻和非均勻分布的FGM圓錐殼動撓度的影響，從圖中可看出，F(xiàn)GM截錐殼動撓度峰值隨半錐角的增大而明顯增大。

圖5 半錐角對圓錐殼動力響應(yīng)的影響(N=2,Wc=0.5,α?=0.1,Kw=0,KP=0)

圖6計算了不同彈性地基參數(shù)對孔隙均勻和非均勻分布的FGM圓錐殼動撓度的影響，從圖6可分析出，殼體動撓度峰值隨著地基壓縮參數(shù)、剪切參數(shù)的增大而減小，其中地基剪切參數(shù)對殼動力響應(yīng)的影響比壓縮參數(shù)的影響更為顯著。

圖6 彈性地基參數(shù)對圓錐殼動力響應(yīng)的影響(γ=π/4,N=2,Wc=0.5,α?=0.1)

圖7比較了不同材料組分指數(shù)對孔隙均勻和非均勻分布的FGM圓錐殼動撓度的影響。從圖7可觀察到，殼動撓度峰值隨材料組分指數(shù)的增大而增大。

圖7 材料組分指數(shù)對圓錐殼動力響應(yīng)的影響(γ=π/4,Wc=0.5,α?=0.1,Kw=0,KP=0)

3 結(jié)論

本文基于改進的混合律模型及經(jīng)典薄殼理論，建立了彈性地基中含均勻和非均勻孔隙的FGM圓錐殼的振動方程，求出其在兩端簡支邊界條件下自由振動和動力響應(yīng)的解析解，計算結(jié)果表明：

(1)自振頻率隨地基壓縮和剪切參數(shù)的增大而增大，隨材料組分指數(shù)的增大而減小，動力響應(yīng)隨地基壓縮和剪切參數(shù)、陶瓷材料質(zhì)量組分的增大而減小，隨孔隙率、半錐角、材料組分指數(shù)的增大而增大。

(2)孔隙對殼自振頻率的影響比較復(fù)雜，不僅與孔隙率的大小和分布形式有關(guān)，還與材料的質(zhì)量組分和彈性地基參數(shù)有關(guān)，與非均勻分布孔隙殼體相比，均勻分布孔隙殼體的自振頻率和動力響應(yīng)隨孔隙率的變化更為敏感。

(3)增大彈性地基參數(shù)或半錐角可有效減小圓錐殼的動撓度，且地基的剪切參數(shù)比壓縮參數(shù)的影響更為顯著。