摘要
基于高阶剪切变形梁理论,推导了Winkler-Pasternak弹性基径向功能梯度空心圆柱管自由振动行为的控制方程。该方法无须引入剪切修正系数,自动满足空心圆柱管内外表面剪应力自由边界条件。通过引入辅助函数,将关于挠度和转角的耦合控制方程化为单一高阶微分控制方程。给出了典型边界条件下功能梯度空心圆柱管的频率和振型。将计算结果与已有文献结果对比,验证了所提理论的精度。可以为工程中常见的Winkler-Pasternak弹性基梁结构提供更高精度的一维弹性理论解,研究结果表明,功能梯度材料的梯度参数和弹性地基刚度系数对固有频率值影响显著。与高阶固有频率相比,刚度系数对低阶固有频率的影响更加明显。
关键词
Abstract
Based on the theory of high-order shear deformation beam, the governing equations for the free vibration behavior of Winkler-Pasternak elastic foundation radial function gradient hollow cylindrical pipe were derived. This method does not require the introduction of shear correction coefficient and automatically satisfies the free boundary condition of shear stress on the inner and outer surfaces of the hollow cylindrical pipe. By introducing auxiliary functions, the coupled equations for deflection and angle were transformed into a single high-order differential equation. The frequency and mode shapes of the function gradient hollow cylindrical pipe under typical boundary conditions were given. The calculation results were compared with the results in the existing literature to verify the accuracy of the proposed theory. It can provide higher precision one-dimensional elastic theoretical solutions for the common Winkler-Pasternak elastic beam structures in engineering. Research results show that the gradient parameter and elastic foundation stiffness coefficient of the function gradient material have a significant impact on the natural frequency value. Compared with the high-order natural frequency, the stiffness coefficient has a more significant impact on the low-order natural frequency.
圆柱形梁作为应用广泛的结构件,在工程中发挥着重要作用。与实心圆柱形梁相比,等质量的空心圆柱管可以承受更高荷载。圆柱管力学行为相关理论包括两个方面:一是基于三维弹性理论对空间变量进行求解,获得结构响应,这种方法精度高,但求解过程复杂[1-2]。二是基于一维弹性理论展开求解,这种方法变量少,计算效率高。
高阶剪切变形梁理论(简称高阶梁理论)最早由Levinson[3]提出。该理论同时考虑剪切变形和转动惯量,且无须剪切修正系数。现阶段,高阶梁理论已成为相关领域研究热点。Reddy等[4]提出的适用于弹性壳结构的高阶梁理论,成为许多研究的基础。Huang等[5-6]将高阶梁理论的适用对象从矩形截面梁扩展到圆形截面梁,Ma等[7]进一步将其扩展至环形截面空心圆柱管。
功能梯度材料是一种非均匀复合材料,其组分可以随着结构位置连续变化,因此具备优越的力学性能,得到科研人员的广泛关注[8]。对于功能梯度梁结构,材料属性可能沿轴向或径向变化[9]。在Euler-Bernoulli和Timoshenko梁理论的基础上,Li[10]提出了一种分析功能梯度梁弯曲和振动的方法。邓先琪等[11]提出了一种适用于功能梯度梁动静态力学行为的格林函数方法。王壮壮等[12]基于三阶剪切变形板理论和正弦剪切变形板理论,分析了石墨烯增强功能梯度板的屈曲和弯曲行为。已有关于功能梯度材料梁、板结构的研究大部分局限于自身变形、自由振动和屈曲等问题,关于弹性地基上功能梯度梁结构的研究鲜见报道。
梁在弹性地基中的力学行为是工程中的普遍问题。常见的地基模型有:Winkler模型、弹性连续介质模型、黏弹性模型、双参数模型等。双参数模型也被称为Winkler-Pasternak模型,克服了Winkler地基模型不能反映压力扩散的缺陷,数学处理难度也显著低于弹性连续介质模型和黏弹性模型。艾智勇等[13]采用分数阶Merchant模型,研究了饱和地基与梁共同作用的时效问题。周凤玺等[14]分析了初始轴向机械力作用下,三参数Winkler-Pasternak黏弹性地基上多孔功能梯度黏弹性材料梁在热环境中的自由振动特性,需要指出的是,其所研究的梁仍属于矩形截面梁,未涉及圆形或圆环形截面梁。蒲育等[15-16]对Winkler-Pasternak弹性基上梁结构等振动和屈曲问题开展了系统研究。现阶段,关于弹性地基上的功能梯度空心圆柱管振动行为的研究尚未开展。
本文基于高阶梁理论,研究双参数弹性基功能梯度空心圆柱管的自由振动问题。通过引入翘曲形状函数,建立关于挠度和转角的耦合控制方程。借助辅助函数,得到单一控制微分方程。给出了典型边界条件下的频率方程,讨论了梯度、刚度系数和长径比等参数对固有频率的影响。
1 理论公式推导
考虑在双参数弹性基之上放置长度为L,内外半径分别为Ri和Ro,材料属性沿径向梯度分布的空心圆柱管,如图1所示,M表示弯矩,Q表示剪力,q表示横向均布载荷。

图1Winkler-Pasternak弹性基功能梯度圆柱管示意图
Fig.1Schematic of a FG circular piperesting on a Winkler-Pasternak elastic foundation
功能梯度圆柱管的材料属性取决于半径。同时采用笛卡尔坐标系(x,y,z)和柱坐标(x,r,θ)。取中性轴为x轴,z轴向上。(u,v,w)和(u,wr,wθ)分别是笛卡尔坐标系和柱坐标下的位移分量,存在如下关系:
(1)
高阶梁理论中轴向位移分量u可以表示为挠度w和转角ψ的函数。
(2)
其中,u0代表轴力引起的轴向拉伸或压缩位移。本文不考虑轴力作用的面内横向振动问题,因此u0=0,挠度w是变量x和时间t的函数。f(y,z)为截面翘曲形状函数,用于描述横截面的翘曲形状。根据几何方程,可将应变分量表示为挠度w和转角ψ的函数。
(3)
(4)
(5)
本文选择多项式形式的翘曲形状函数:
(6)
将式(6)代入式(5),结合本构方程可知,当r=Ri,Ro时σxr=0,即空心圆柱管内外表面剪应力自由。根据应力与内力的关系,进一步获得M和Q的表达式。
(7)
(8)
其中,(EI)*、(EI)0、(GA)*分别是与弹性模量和剪切模量有关的积分常数,其积分表达式分别为
(9)
(10)
(11)
其中,E(r)和G(r)分别表示功能梯度圆柱管的弹性模量和剪切模量。对于双参数弹性基上的空心圆柱管,横向和纵向平衡方程分别为
(12)
(13)
其中:K1和K2分别为双参数弹性基的平移和旋转刚度。
(14)
(15)
(16)
ρ(r)表示功能梯度圆柱管的密度。若K2=0,则为Winkler弹性基。若K1=K2=0,则表示结构处于自由空间。将内力表达式(7)~(8)代入平衡方程(12)~(13),可得关于w和转角ψ的耦合控制方程。
(17)
(18)
为了将方程(19)和(20)解耦,引入辅助函数F(x),将挠度w和转角ψ分别写成如下形式:
(19)
(20)
将式(19)和式(20)分别代入式(17)和式(18)中,可将耦合控制方程化为下述单一控制方程。
(21)
其中
(22)
(23)
(24)
(25)
(26)
2 振动问题求解
对于自由振动,F(x)的解可以表示为
(27)
式中,ω为圆频率。为求解方便,引入下述无量纲量。
(28)
(29)
将式(27)~(29)中的无量纲量代入式(21),令q=0,可得
(30)
其中
(31)
(32)
(33)
本文重点关注工程中常见的低频振动问题,因此,可将式(30)解的形式表示为
(34)
其中,Cj(j=1,···,4)为未知系数。
(35)
低频振动的频率范围如式(36)所示:
(36)
也可写为
(37)
式中,截断圆频率Ωcut-off可表示为:
(38)
当圆频率Ω的值大于截断圆频率Ωcut-off时,圆柱管发生高频振动。λ2的值从实数变为虚数,式(34)中解的双曲函数变为三角函数。取特殊情况k1=k2=0,功能梯度空心圆柱管的截断圆频率可以表示为
(39)
结合式(29),也可将式(39)表示为
(40)
对于各向同性材料,式(39)和式(40)将分别退化为
(41)
(42)
其中,积分常数(GA)*和(EI)0分别退化为GA和EI1,I1、I和A的表达式如下:
(43)
与已有文献[8]中的结果一致。
至此,如式(34)所示,本文已求得辅助函数的表达式。将式(34)分别代入式(19)、式(20)、式(7)和式(8)中,可以获得挠度w、转角ψ、弯矩M和剪力Q的表达式。
(44)
(45)
(46)
(47)
其中
(48)
(49)
(50)
(51)
根据边界条件,可得关于位移和内力的方程组。为使未知系数Cj(j=1,···,4)存在非零解,需令系数矩阵行列式等于零,得到频率方程。频率方程通常是具有无限个正根的超越方程,对应无限个固有频率。作为例证,给出悬臂空心圆柱管的求解过程,边界条件如下:
(52)
利用式(52)的前三项,可得系数间关系如下:
(53)
(54)
(55)
将式(53)~(55)代入式(52)的第四项,可得:
(56)
其中
(57)
式(56)即为双参数弹性基上悬臂空心圆柱管的频率方程。求解该方程可得固有频率的精确值。一旦求得固有频率,就可通过挠度和转角表达式给出振型。
对于其他边界条件,可类似获得频率方程和振型结果。表1和表2分别给出了两端简支、固支-自由和两端固支三种边界条件下的频率方程和振型(包括挠度和转角)。需要指出的是,表1中的频率方程与文献[8]中相关特征方程形式类似。不同的是,本文的特征值λ1、λ2考虑了弹性基刚度系数K1和K2的影响。如果令K1=K2=0,本文所有计算结果可以退化到文献[8]中的相关结果。
表1Winkler-Pasternak弹性基空心圆柱管频率方程
Tab.1 Frequency equation for hollow cylindrical pipes on Winkler-Pasternak elastic foundation

表2Winkler-Pasternak弹性基空心圆柱管振型方程
Tab.2 Vibration mode equation for hollow cylindrical pipes on Winkler-Pasternak elastic foundation

3 结果与讨论
3.1 正确性验证
本文研究径向功能梯度圆柱管,材料属性服从幂函数规律,即
(58)
其中,Yi和Yo分别代表内、外表面的材料特性。N为梯度参数,代表幂函数形式的体积分数。下文数值计算中,内外表面材料分别设为铝和氧化锆,材料属性如下:
1)Al:Ei=70 GPa, Gi=27 GPa, ρi=2 702 kg/m3。
2)ZrO2: Eo=200 GPa, Go=77 GPa, ρo=5 700 kg/m3。
为了验证本文理论的准确性,将所提理论的结果与已有文献的结果对比。表3和表4分别给出不同k1、k2时,两端简支梯度圆柱管前六阶无量纲频率的值。表3和表4中计算结果对应模型参数分别为:k2=0,L/Ro=10,Ri/Ro=0.5,N=1;k1=1,L/Ro=10,Ri/Ro=0.5,N=1。需要指出的是,表3中k2=0工况所对应的结果为文献[8]中相关结果。可以发现,当k1、k2趋近于0时,本文理论计算结果逐渐趋近于文献[8]的计算结果,这与预期相符。
仔细观察表3和表4的数据可以发现,无量纲频率的值随着弹性地基刚度系数的增大而增大。但弹性地基刚度系数对低阶频率的影响更为显著,对基频的影响尤为显著。随着阶数增大,地基刚度系数对频率的影响逐渐减小。另外,Pasternak参数k2对固有频率的影响比Winkler参数k1更明显,这意味着旋转刚度对固有频率的影响比平移刚度更加显著。Najafi等[17]发现Pasternak参数对挠度的影响比Winkler参数更加明显。本文关于k2对振动行为影响的结论,与Najafi的研究结论并不冲突,二者均体现了地基刚度对于系统刚度的影响,与Winkler参数相比,Pasternak参数对弹性地基抗剪刚度的影响更大。所以k2对最大挠度和固有频率的影响更加显著。
表3不同k1时简支-简支梯度圆柱管前六阶无量纲频率
Tab.3 The first six dimensionless frequencies of simply supported-simply supported gradient cylindrical tubes for different k1

表4不同k2时简支-简支梯度圆柱管前六阶无量纲频率
Tab.4 The first six dimensionless frequencies of simply supported-simply supported gradient cylindrical tubes for different k2

3.2 梯度参数和刚度系数对频率的影响
本节研究梯度参数和弹性地基刚度系数对固有频率的影响。根据表1中的频率方程,表5~7分别给出两端简支、两端固支和固支-自由功能梯度空心圆柱管的无量纲固有频率。需表5~7中计算结果对应的模型参数为:k1=1,k2=1,L/Ro=10,Ri/Ro=0.5。根据表中结果可知,随着梯度参数N的增大,无量纲频率逐渐降低。这是因为N增大时,功能梯度空心圆柱管的主要组分材料从氧化锆逐渐过渡为铝。为了详细比较梯度参数对不同阶频率的影响,表5~7提供了不同梯度参数与梯度参数N=0对应无量纲固有频率之间的相对误差εp。可以发现,梯度参数N对基频的影响是最小的。随着模态数的增大,影响效果略有增大,这一现象对于文中每种边界条件均一致。
基频对于振动问题分析至关重要。图2和图3(对应模型参数:L/Ro=10,Ri/Ro=0.5,N=1)分别绘制了无量纲基频随刚度系数k1和k2的变化情况,ω1代表基频。可以发现基频随着刚度系数增大而增大,与文献[18]中的结论一致。比较图2和图3,可以发现Pasternak系数k2对基频的影响大于Winkler系数k1。
表5不同梯度参数时简支-简支梯度圆柱管前三阶无量纲频率与相对误差εp
Tab.5 The first three dimensionless frequencies and relative error εp of double simply supported FG cylindrical tubes for different gradient index N

表6不同梯度参数时固支-固支梯度圆柱管无量纲频率与相对误差εp
Tab.6 The dimensionless frequency and the relative error εp of clamped-clamped gradient cylindrical tubes for different gradient parameters

表7不同梯度参数时固支-自由梯度圆柱管无量纲频率与相对误差εp
Tab.7 The dimensionless frequency and the relative error εp of clamped-free gradient cylindrical tubes for different gradient parameters


图2简支-简支梯度圆柱管无量纲基频随k1变化情况
Fig.2The variation of the dimensionless fundamental frequency with k1 for simply supported-simply supported gradient cylindrical tubes

图3简支-简支梯度圆柱管无量纲基频随k2变化情况
Fig.3The variation of the dimensionless fundamental frequency with k2 for simply supported-simply supported gradient cylindrical tubes
图4和图5分别绘制了不同刚度系数k1和k2时,两端简支功能梯度空心圆柱管的基频ω1随长度L的变化情况。对应模型参数分别为:k2=0,Ro=1,Ri/Ro=0.5,N=1和k1=1,Ro=1,Ri/Ro=0.5,N=1。显然,随着长度L的增大,基频迅速降低,而刚度系数k1和k2则会不同程度地提高基频的值。
获得固有频率后,可以根据表2中的公式确定挠度和转角。图6和图7给出了不同系数k1和k2时的模态振型,对应模型参数分别为:k2=0,Ri/Ro=0.5,L/Ro=10,N=1和k1=1,Ri/Ro=0.5,L/Ro=10,N=1。将每个振型归一化处理,可以发现,Pasternak刚度系数k2对振型的影响比Winkler刚度系数k1更加明显。原因是k2对剪切刚度的影响高于k1,剪切刚度越高,对应的模态位移越小。随着阶数的提高,刚度系数对振型的影响更加明显。

图4不同k1时简支-简支梯度圆柱管基频随长度变化情况
Fig.4Variation of the fundamental frequency with length for simply supported-simply supported gradient cylindrical tubes at different values of k1

图5不同k2时简支-简支梯度圆柱管基频随长度变化情况
Fig.5Variation of the fundamental frequency with length for simply supported-simply supported gradient cylindrical tubes at different values of k2

图6不同k1对应的前三阶振型
Fig.6The first three modes of vibration corresponding to different k1

图7不同k2对应的前三阶振型
Fig.7The first three modes of vibration corresponding to different k2
4 结论
本文基于高阶梁理论,研究了双参数弹性基功能梯度空心圆柱管自由振动行为。通过引入自动满足表面剪应力自由边界条件的翘曲形状函数,推导了双参数弹性基上功能梯度空心圆柱管振动行为的控制方程,求得了几种典型边界条件下自由振动频率方程。通过数值算例验证了模型的正确性。得出下述一些结论:
1)弹性地基刚度系数的增大会使频率增大。
2)无量纲频率值随着地基刚度系数的增大而增大。但弹性地基刚度系数对低阶频率的影响更为明显,是对基频的影响最为显著。
3)Pasternak系数k2对固有频率和振型的影响比 Winkler系数k1更显著,即旋转刚度对振动行为的影响比平移刚度更显著。
4)梯度参数N显著影响无量纲频率,可以通过调整梯度参数获得期望的固有频率。
5)长径比较小时基频受刚度系数的影响更大。随着长径比增大,弹性刚度系数的影响逐渐减小。