功能梯度材料是一种组分随空间连续变化且具有特定功能的新型复合材料,具有耐高温、抗疲劳、强度高、韧性好等优点。高性能功能梯度圆柱壳结构在航空发动机、航天飞机、火箭及化工设备等具有广泛应用前景,正日益取代单一材料的普通圆柱壳结构,成为该结构形式新的发展趋势[1]。
Loy等[2]利用Love壳体理论的应变-位移几何描述分析了功能梯度圆柱壳的自由振动特性,计算了两端简支功能梯度圆柱壳的频率,讨论了材料体积分数和组分材料配置对频率的影响。杜长城等[3]利用Donnell壳体理论,推导了功能梯度材料薄壁圆柱壳线性振动的简化控制方程,分析了自由振动特性。Haddadpour[4]和Sheng[5]等研究了热环境中功能梯度圆柱壳的自由振动。项松等[6]计算了简支各向同性和功能梯度圆柱壳的频率,研究了幂指数、x和θ方向的波数、厚径比对简支旋转功能梯度圆柱壳频率的影响。项爽[7]分析了各种参数对旋转功能梯度圆柱壳的振动特性的影响。
由文献[1, 2, 3, 4, 5, 6, 7, 8, 9]可知,前人对旋转功能梯度圆柱壳结构的研究集中考虑了经典边界条件的组合。事实上,边界条件作为影响结构振动特性的重要因素,有必要探究弹性边界约束对旋转功能梯度圆柱壳结构自由振动特性的影响。
本文采用一种改进傅里叶级数方法建立了弹性约束边界条件下旋转功能梯度圆柱壳结构振动分析模型。基于Love壳体理论,采用一种改进傅里叶级数方法[10, 11, 12, 13]和Rayleigh-Ritz法对旋转功能梯度圆柱壳自由振动进行建模,推导出频率计算方程。通过算例,分析了弹性边界约束对振动固有特性的影响。
1 理论建模图 1所示为弹性边界约束下旋转功能梯度圆柱壳结构,R、H、L和Ω为圆柱壳中曲面半径、厚度、长度和转速,在其中曲面上建立参考柱坐标系(x,θ,z),u、v、w分别为圆柱壳中曲面沿x、θ、z方向的位移。k1、k2、k3和k4分别为施加于圆柱壳左端x、θ、z方向的线性约束弹簧和z方向的旋转约束弹簧;k5、k6、k7和k8分别为施加于圆柱壳右端x、θ、z方向的线性约束弹簧和z方向的旋转约束弹簧。通过将弹簧刚度系数设置为零或无穷大,可以得到任意经典边界条件,简支、自由或固支,如表 1所示;将弹簧刚度设置为不同的数值,则可模拟任意弹性约束。
图 1截面显示圆柱壳材料为功能梯度材料,其物理性能参数表示为沿壳体厚度呈梯度变化形式:
$\left\{ \begin{array}{l} E=\left({{E_o} - {E_i}} \right){\left({\frac{{2z+H}}{{2H}}} \right)^N}+{E_i}\\ \mu=\left({{\mu _o} - {\mu _i}} \right){\left({\frac{{2z+H}}{{2H}}} \right)^N}+{\mu _i}\\ \rho=\left({{\rho _o} - {\rho _i}} \right){\left({\frac{{2z+H}}{{2H}}} \right)^N}+{\rho _i} \end{array} \right.$ | (1) |
式中:Ei和Eo分别为内、外表面材料杨氏模量,μi和μo分别为内、外表面泊松比,ρi和ρo分别为内、外表面质量密度,N为非负实数幂指数。可以发现:当z=-H/2时,E=Ei,μ=μi,ρ=ρi,圆柱壳材料为单一外表面材料;当z=H/2时,E=Eo,μ=μo,ρ=ρo,圆柱壳材料为单一内表面材料。
1.1 能量原理描述旋转功能梯度圆柱壳动能表达式为
$\begin{array}{l} T=\frac{1}{2}\int_{ - H/2}^{H/2} {\int_0^L {\int_0^{2\pi } {\rho \left[{{{\dot u}^2}+{{\dot v}^2}+{{\dot w}^2}+} \right.} } } \\ \left.{2\Omega \left({w\dot v - v\dot w} \right)+{\Omega ^2}\left({{v^2}+{w^2}} \right)} \right]R{\rm{d}}x{\rm{d}}\theta {\rm{d}}z \end{array}$ | (2) |
功能梯度圆柱壳结构自由振动的应变能描述为
$\begin{array}{l} {V_\varepsilon }=\frac{1}{2}\int_{ - H/2}^{H/2} {\int_0^L {\int_0^{2\pi } {\rho \left[{{\delta _x}\varepsilon _x^0+} \right.} } } {\delta _\theta }\varepsilon _\theta ^0+{\delta _{x\theta }}\varepsilon _{x\theta }^0+\\ \;\;\;\;\;\;\;\;\;\left.{z\left({{\delta _x}\varepsilon _x^0+{\delta _\theta }k_\theta ^0+{\delta _{x\theta }}{\tau ^0}} \right)} \right]R{\rm{d}}x{\rm{d}}\theta {\rm{d}}z \end{array}$ | (3) |
式中:σx、σθ、τxθ和εx、εθ、εxθ分别表示点(x,θ,z)处的应力分量和应变分量。
根据Love薄壳理论,位移-应变关系为
${\varepsilon _x}=\varepsilon _x^0+zk_x^0,{\varepsilon _\theta }=\varepsilon _\theta ^0+zk_\theta ^0,{\varepsilon _{x\theta }}=\varepsilon _{x\theta }^0+2z{\tau ^0}$ | (4) |
圆柱壳中曲面上的应变分量εx0、εθ0、εxθ0分别可表示为
$\varepsilon _x^0=\frac{{\partial u}}{{\partial x}},\varepsilon _\theta ^0=\frac{{\partial v}}{{R\partial \theta }}+\frac{w}{R},\varepsilon _{x\theta }^0=\frac{{\partial u}}{{R\partial \theta }}+\frac{{\partial u}}{{\partial x}}$ | (5) |
其中,曲面曲率分量κx0、κθ0、τ0为
$\left\{ {\begin{array}{*{20}{c}} {k_\theta ^0=- \frac{{{\partial ^2}w}}{{\partial {x^2}}}}\\ {k_\theta ^0=\frac{1}{{{R^2}}}\left({\frac{{\partial v}}{{\partial \theta }} - \frac{{{\partial ^2}w}}{{\partial {\theta ^2}}}} \right)}\\ {{\tau ^0}=\frac{1}{R}\left({\frac{{\partial v}}{{\partial x}} - \frac{{{\partial ^2}w}}{{\partial x\partial \theta }}} \right)} \end{array}} \right.$ | (6) |
平面应力-应变关系表示为
$\left[{\begin{array}{*{20}{c}} {{\sigma _x}}\\ {{\sigma _\theta }}\\ {{\sigma _{x\theta }}} \end{array}} \right]=\left[{\begin{array}{*{20}{c}} {{Q_{11}}}&{{Q_{12}}}&0\\ {{Q_{12}}}&{{Q_{22}}}&0\\ 0&0&{{Q_{66}}} \end{array}} \right]\left[{\begin{array}{*{20}{c}} {{\varepsilon _x}}\\ {{\varepsilon _\theta }}\\ {{\varepsilon _{x\theta }}} \end{array}} \right]$ | (7) |
其中
${Q_{11}}={Q_{22}}=\frac{E}{{1 - {\mu ^2}}},{Q_{12}}=\frac{{E\mu }}{{1 - {\mu ^2}}},{Q_{66}}=\frac{E}{{2\left({1+\mu } \right)}}$ | (8) |
由旋转引起的旋转功能梯度圆柱壳的初应变能表达式为
$\begin{array}{l} {V_h}=\frac{{{\Omega _2}}}{2}\int_{ - H/2}^{H/2} {\int_0^L {\int_0^{2\pi } {\rho \left[{{{\left({\frac{{\partial u}}{{\partial \theta }}} \right)}^2}+{{\left({\frac{{\partial v}}{{\partial \theta }}} \right)}^2}+{{\left({\frac{{\partial w}}{{\partial \theta }}} \right)}^2}} \right.} } }+\\ \;\;\;\;\;\;\;\;\;\left.{{v^2}+{w^2}+2w\frac{{\partial v}}{{\partial \theta }} - 2v\frac{{\partial w}}{{\partial \theta }}} \right]R{\rm{d}}x{\rm{d}}\theta {\rm{d}}z \end{array}$ | (9) |
储存于边界约束弹簧中的弹性势能表达式为
$\begin{array}{l} Vs=\frac{1}{2}{\int_0^{2\pi } {\left\{ {\left[{{k_1}{u^2}+{k_2}{v^2}+{k_3}{w^2}+{k_4}{{\left({\frac{{\partial w}}{{\partial x}}} \right)}^2}} \right]} \right.} _{x=0}}+\\ \;\;\;\;\;\;\;\;\;\left.{{{\left[{{k_5}{u^2}+{k_6}{v^2}+{k_7}{w^2}+{k_8}{{\left({\frac{{\partial w}}{{\partial x}}} \right)}^2}} \right]}_{x=L}}} \right\}R{\rm{d}}\theta \end{array}$ | (10) |
将下式:
${A_{ij}}=\int_{ - H/2}^{H/2} {{Q_{ij}}{\rm{d}}z,} {B_{ij}}=\int_{ - H/2}^{H/2} {{Q_{ij}}{\rm{d}}z} $ | (11) |
${D_{ij}}=\int_{ - H/2}^{H/2} {{Z^2}{Q_{ij}}{\rm{d}}z,} \rho t=\int_{ - H/2}^{H/2} {\rho {\rm{d}}z} $ | (12) |
代入动能表达式(2)、应变能表达式(3)和初应变能表达式(9)中,即可消除应变能、初应变能和动能表达式中关于径向坐标z的积分。式中:i,j=1,2,6。
1.2 位移场函数表达将旋转功能梯度圆柱壳结构的振动位移表示成环向模态和轴向模态相乘的形式:
$\left\{ {\begin{array}{*{20}{c}} {u{\rm{=cos}}\left({n\theta {\rm{+}}\omega t} \right){U_n}\left(x \right)}\\ {v{\rm{=sin}}\left({n\theta {\rm{+}}\omega t} \right){V_n}\left(x \right)}\\ {w{\rm{=cos}}\left({n\theta {\rm{+}}\omega t} \right){W_n}\left(x \right)} \end{array}} \right.$ | (13) |
其中:Un、Vn、Wn以如下改进傅里叶级数形式来构造:
$\left\{ {\begin{array}{*{20}{c}} {{U_n}=\sum\limits_{m{\rm{=0}}}^\infty {{A_m}{\rm{COS}}{\lambda _m}x+a{\xi _1}\left(x \right)+b{\xi _2}\left(x \right)} }\\ {{V_n}=\sum\limits_{m{\rm{=0}}}^\infty {{B_m}{\rm{COS}}{\lambda _m}x+c{\xi _1}\left(x \right)+d{\xi _2}\left(x \right)} }\\ {{W_n}=\sum\limits_{m{\rm{=0}}}^\infty {{C_m}{\rm{COS}}{\lambda _m}x+e{\xi _1}\left(x \right)+f{\xi _2}\left(x \right)+} }\\ {g{\xi _3}\left(x \right)+h{\xi _4}\left(x \right)} \end{array}} \right.$ | (14) |
式中:λm=mπ/L,ξ1(x)、ξ2(x)、ζ1(x)、ζ2(x)、ζ3(x)、ζ4(x)为辅助函数。
1.3 系统特征方程推导将位移表达式(13)代入1.1节得到的应变能、初应变能和动能表达式以及弹性势能表达式(2)、(3)、(9)、(10),并对含有θ的函数在[0,2π]进行积分,最终得到只含关于轴向坐标x的积分的能量表达式。根据能量守恒:
${T_{\max }}={\left\{ {{V_\varepsilon }+{V_h}+{V_s}} \right\}_{\max }}$ | (15) |
将能量表达式代入式(15),应用Rayleigh-Ritz法进行推导、变换,可得到关于傅里叶系数的系统特征方程:
$\left({K - {\omega ^2}M - \omega \Omega {M_1} - {\Omega ^2}{M_2}} \right)\left[{\begin{array}{*{20}{c}} {\bar U}\\ {\bar V}\\ {\bar W} \end{array}} \right]=0$ | (16) |
式中:K、M、M1和M2分别为刚度矩阵和质量矩阵,U、V和W分别为3个方向的傅里叶系数向量。
求解旋转功能梯度圆柱壳结构系统特征方程(16),可得到自由振动的固有频率。设置不同的边界弹簧刚度,可模拟任意弹性边界,进而可对弹性边界约束对固有特性的影响进行分析。需要指出的是,圆柱壳旋转产生的科氏力使得频率随转速的变化而发生分岔,产生前行波和后行波:前行波传播方向与旋转方向一致,求解得到的固有频率值ff小于0;后行波传播方向与旋转方向相反,求解得到的固有频率值fb大于0。
2 数值结果与分析除特别说明外,以下计算分析中研究外表面为不锈钢、内表面为镍的功能梯度圆柱壳结构。功能梯度材料参数沿壳体厚度呈指数变化,幂指数用N表示,圆柱壳轴向波数用m表示,所有计算中N=1, m=1。温度为300 K时,不同材料物理参数如表 2。
表 3给出了两端简支(S-S)条件下某种功能梯度圆柱壳的无量纲频率参数,并与文献[10, 11]中的结果进行了对比。圆柱壳长径比L/R=20,厚径比H/R=0.01,内外表面材料一致:杨氏模量E=1.878 9×1011 Pa,泊松比μ=0.308 9,质量密度ρ=6 93 3 kg/m3。结果用无量纲频率参数ω'表示为
$\omega '=\omega R\sqrt {\rho \left({1 - \mu '} \right)/E}$ | (17) |
n | 文献[9] | 文献[8] | 本文方法 |
1 | 0.016 052 | 0.015 949 | 0.016 052 |
2 | 0.009 369 | 0.009 081 | 0.009 355 |
3 | 0.022 102 | 0.021 952 | 0.022 081 |
4 | 0.042 093 | 0.042 001 | 0.042 071 |
5 | 0.068 007 | 0.067 939 | 0.067 984 |
6 | 0.099 729 | 0.099 675 | 0.099 705 |
7 | 0.137 238 | 0.137 192 | 0.137 210 |
8 | 0.180 527 | 0.180 486 | 0.180 500 |
9 | 0.229 593 | 0.229 556 | 0.229 570 |
10 | 0.284 435 | 0.284 400 | 0.284 410 |
观察表格可知,本文所用方法得到的结果与文献[8, 9]中的结果均吻合良好,故用本文方法求解旋转功能梯度圆柱壳结构自由振动固有频率的正确性得到了验证,后续分析结果是可信的。由于本文方法在对振动位移的表达中采用了一种改进的傅里叶级数进行展开,因此需要在实际计算中将展开项数有限截断,以将方程(18)的维数化为有限维。计算表明,当截断数为10时,结果已经足够精确。
2.2 单个弹性约束和转速对旋转功能梯度圆柱壳自由振动特性影响现有文献对功能梯度圆柱壳结构的研究大多只针对固支-固支(C-C)、固支-简支(C-S)、简支-简支(S-S)等经典边界条件,而在实际应用中圆柱壳所处的边界条件更为复杂,更多情况下是任意弹性边界。因而有必要进行弹性边界约束刚度对旋转功能梯度圆柱壳自由振动特性的影响分析。本节中约束弹簧刚度均以无量纲刚度的形式表示:
$\left\{ {\begin{array}{*{20}{c}} {{K_1}'={K_1}/{A_{11}},{K_2}'={K_2}/{A_{11}}}\\ {{K_3}'={K_3}/{A_{11}},{K_4}'={K_4}/{A_{11}}}\\ {{K_5}'={K_5}/{A_{11}},{K_6}'={K_6}/{A_{11}}}\\ {{K_7}'={K_7}/{A_{11}},{K_8}'={K_8}/{A_{11}}} \end{array}} \right.$ | (18) |
为了探讨弹性边界约束刚度和转速同时变化对功能梯度圆柱壳的影响,将长径比L/R=10,厚径比H/R=0.002的圆柱壳左右两端设置相同的弹性边界约束,依次设置为仅有轴向弹簧k'1、环向弹簧k'2、径向弹簧k'3或径向旋转弹簧k'4的情况,无量纲弹簧刚度从10-4变化到104,转速由0变化到120 rad/s。观察n=2时,功能梯度圆柱壳前、后行波频率ff和fb的变化情况。
观察图 2可知,4种不同弹性边界约束情况下,随着转速Ω的增大,均出现旋转功能梯度圆柱壳结构自由振动的前行波频率ff逐渐减小、后行波频率fb逐渐增大且后行波频率变化更大的现象;随着弹性约束刚度的增大,前后行波频率均呈增大趋势。而边界设置为仅有轴向弹性约束k'1、环向弹性约束k'2或径向弹性约束k'3时,弹性约束刚度的增大对前后行波频率的影响明显,当边界设置为径向旋转约束k'4时,刚度的变化对频率影响很小,几乎可以忽略。边界为轴向约束弹簧k'1时,刚度由10-2变化到101过程中,前后行波频率快速增大;边界为环向弹性约束k'2时,刚度在10-3~100区间变化引起的频率变化较大;边界为径向弹性约束k'3时,频率在刚度从10-4变化到10-1时有大幅增长。因此,旋转功能梯度圆柱壳结构自由振动行波频率对不同的弹性约束弹簧的敏感范围不同。而对比四种边界设置的频率变化量可知,环向弹性约束k2和径向弹性约束k3对频率的影响最大,轴向弹性约束k1的影响较小,径向旋转约束k4对频率几乎无影响。
2.3 双弹性约束对旋转功能梯度圆柱壳自由振动特性影响上述分析只针对单个弹性约束的刚度变化,对于双弹性约束同时变化的情况,分析了轴向约束弹簧k1和径向旋转约束弹簧k4同时变化或者环向约束弹簧k2 和径向约束弹簧k3同时变化对旋转功能梯度圆柱壳结构自由振动频率的影响,如图 3。
图 3(a)和图 3(c)为环向约束刚度k'2和径向约束刚度k'3均为1012时,轴向弹性约束刚度k'1和径向旋转弹性约束刚度k'4分别从10-4变化到103、从10-4变化到102对旋转功能梯度圆柱壳自由振动前后行波频率的影响。当k'1= k'4=10-4时,近似为简支边界(S-S);当k'1= k'4=102时,近似为固支边界(C-C)。图中,当径向旋转约束刚度k'4逐渐增大时,前后行波频率只在开始时略有增大,增幅很小,而随着轴向弹性约束刚度k'1逐渐增大,前后行波频率均有大幅增加。故轴向弹性约束k1对旋转功能梯度圆柱壳结构自由振动频率的影响远大于径向旋转弹性约束k4,且固支边界(C-C)行波频率大于简支边界(S-S)。
图 3(b)和图 3(d)中,圆柱壳左右两端边界设置相同,轴向弹性约束刚度k'1和径向旋转弹性约束刚度k'4均设置为0,环向弹性约束刚度k'2和径向弹性约束刚度k'3分别从10-4变化到102。当k'2= k3'=10-4时,近似为自由边界(F-F);当k'2= k3'=102时,近似为简支边界(S-S)。由图可知,当弹性约束刚度k'2、 k3'从10-4变化到10-1时,前后行波频率均快速增大,且随着k'3的增大而增大的幅度大于随着k'2的增大而增大的幅度。即环向弹性约束k2和径向弹性约束k3同时作用时,径向弹性约束k3对旋转功能梯度圆柱壳结构自由振动影响更大,且简支边界(S-S)行波频率大于自由边界(F-F)。
3 结论本文应用一种改进傅里叶级数建立了旋转功能梯度圆柱壳结构自由振动特性预报模型,通过在圆柱壳结构两端引入4种约束弹簧统一模拟边界条件。旋转功能梯度圆柱壳结构3个方向位移场采用改进傅里叶级数形式构造,辅助边界特性函数用于改善级数解的收敛性与精确性。基于Love壳体理论,建立了旋转功能梯度圆柱壳结构自由振动行波频率系统特征方程,运用MATLAB语言编程仿真,算例中本文结果与现有文献结果吻合良好,充分验证了本文模型的正确性与可靠性。
计算分析了单一弹性约束和转速同时变化对功能梯度圆柱壳结构自由振动频率的影响规律,发现:
1)随着转速的增大前行波频率逐渐减小而后行波频率逐渐增大,随着各种弹性约束的增大前后行波频率均有所增大,行波频率对不同的弹性约束弹簧的敏感区域不同。
2)环向弹性约束和径向弹性约束对频率影响较明显,轴向弹性约束对频率影响较小,而径向旋转弹性约束对频率几乎无影响。
3)讨论了双弹性约束同时变化对旋转功能梯度圆柱壳结构自由振动频率的影响,发现固支边界下的频率最大,简支边界其次,自由边界最小,径向弹性约束对频率的影响程度大于环向弹性约束。
[1] | 梁斌, 李戎, 张伟, 等. 功能梯度材料圆柱壳的振动特性研究[J]. 船舶力学, 2011, 15(1): 109-117. LIANG Bin, LI Rong, ZHANG Wei, et al. Vibration characteristics of functionally graded materials cylindrical shells[J]. Journal of ship mechanics, 2011, 15(1): 109-117. |
[2] | 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. |
[3] | 杜长城, 李映辉. 功能梯度薄壁圆柱壳的自由振动[J]. 动力学与控制学报, 2010, 8(3): 219-223. DU Changcheng, LI Yinghui. Free vibration of functionally graded cylindrical thin shells[J]. Journal of dynamics and control, 2010, 8(3): 219-223. |
[4] | HADDADPOURT H, MAHMOUDKHANI S, NAVAZI H M. Free vibration analysis of functionally graded cylindrical shells including thermal effects[J]. Thin-walled structures, 2007, 45(6): 591-599. |
[5] | SHENG G G, WANG X. Thermoelastic vibration and buckling analysis of functionally graded piezoelectric cylindrical shells[J]. Applied mathematical modelling, 2010, 34(9): 2630-2643. |
[6] | 项松, 李广超, 张魏, 等. 旋转功能梯度圆柱壳的固有频率计算[J]. 应用数学和与力学, 2012, 33(3): 332-341. XIANG Song, LI Guangchao, ZHANG Wei, et al. Natural frequency of rotating functionally graded cylindrical shells[J]. Applied mathematics and mechanics, 2012, 33(3): 332-341. |
[7] | 项爽. 旋转功能梯度材料圆柱壳的振动特性研究[D]. 郑州: 河南科技大学, 2013. XIANG Shuang. Free vibration of rotating functionally graded cylindrical shells[D]. Zhengzhou: University of Science and Technology of He'nan, 2013. |
[8] | ARSHAD S H, NAEEM M N, SULTANA N, et al. Effects of exponential volume fraction law on the natural frequencies of FGM cylindrical shells under various boundary conditions[J]. Archive of applied mechanics, 2011, 81(8): 999-1016. |
[9] | HOSSEINI-HASHEMI S H, ILKHANI M R, FADAEE M. Accurate natural frequencies and critical speeds of a rotating functionally graded moderately thick cylindrical shell[J]. International journal of mechanical science, 2013, 76: 9-20. |
[10] | DU Jingtao, LI W L, JIN Guoyong, et al. An analytical method for the in-plane vibration analysis of rectangular plates with elastically restrained edges[J]. Journal of sound and vibration, 2007, 306(3/4/5): 908-927. |
[11] | LI W L, ZHANG Xuefeng, DU Jingtao, et al. An exact series solution for the transverse vibration of rectangular plates with general elastic boundary supports[J]. Journal of sound and vibration, 2009, 321(1/2): 254-269. |
[12] | DAI Lu, YANG Tiejun, LI W L, et al. Dynamic analysis of circular cylindrical shells with general boundary conditions using modified Fourier Series method[J]. Journal of vibrationSound and acoustics vibration, 2012, 134(4): 041004. |
[13] | JIANG Shiliang, YANG Tiejun, LI W L, et al. Vibration analysis of doubly curved shallow shells with elastic edge restraints[J]. Journal of vibration and acoustics, 2013, 1356(3135): 034502. |