转动薄壁圆柱壳行波振动响应分析

本资料包含pdf文件1个，下载需要1积分

转动薄壁圆柱壳在诸多工程机械中应用广泛，如高速离心机、航空发动机高速转动鼓筒等。转动所带来的科氏力、离心惯性力及环向初应力，使得转动薄壁圆柱壳具有不 同于静态 圆柱壳 的动力学 特性。具体说来，-方面，离心惯性力使壳体结构产生环向初应力，导致壳体刚度增加，从而使转动薄壁圆柱壳频率随转速的升高而增加；另-方面，受转动速度矢量与变形速度矢量不-致而导致的科氏力的影响，转动薄壁圆柱壳频率随转速的变化发生分岔，产生不同频率的前、后行波，区别于静态圆柱壳频率所对应的驻波(振型)。转动薄壁圆柱壳这种行波振动的特点，使得传统的求解静态圆柱壳振动响应的方法不再适用。因此，对转动薄壁圆柱壳行波振动响应的研究，具有重要的理论意义和应用价值。

国内外众多学者通过解析、数值、实验等手段对转动薄壁圆柱壳的自由振动问题进行了大量的研究工作口叫 ，比较之下，对激振力作用下转动薄壁圆柱壳行波振动响应的研究则相对较少l⌒ 。在为数不多的行波振动响应研究中，以 Huang的工作最具代表性口。 。Huang将两端简支圆柱壳行波模态(traveling mode)代人基于 Love-Timoshenko壳体理论的转动薄壁圆柱壳的振动微分方程中，并利用三角函数的正交性将方程离散为常微分方程组，收稿 日期 ：2011-10-19；修订 日期 ：2013-Ol-10基金项目：国家自然科学基金资助项 目(91216106)通过数值求解常微分方程组，进而采用类似模态叠加的方法求得响应。基于此思路，Huang着重分析和探讨了科氏力对两端简支转动圆柱壳强迫振动的影响以及简谐移动载荷作用下的转动薄壁圆柱壳体的共振现象。需要指出的是，Huang的结果，均针对两端简支边界条件的转动柱壳，而对其他类型的边界条件并没有提及。两端简支约束的圆柱壳的振型可以解析地表达为三角 函数组合 的形式 ，而对于其他类型边界条件下圆柱壳的振型，解析表达式的获取则是困难的，因此有必要发展-种适用于求解各种类型边界约束的转动薄壁圆柱壳行波振动响应的方法。国内也有部分学者在转动薄壁圆柱壳振动响应 的求解方面做了-些工作。如李健等应用Donnel壳体理论ll ，采用复分析的方法研究了转动薄壁悬臂圆柱壳在法向激励作用下的行波振动。

因作者仅通过在静态薄壁圆柱壳上施加反向转动的激励来处理转动柱壳受迫振动问题，因此转动所带来的离心惯性力、科氏力等的影响未能在模型中反映出来，并且作者采用的 Donnel简化壳体理论尽管形式简单，但因其主要考虑法向弯曲变形，忽略了平面内两个方向的惯性力 ，不能准确反映 3个方 向运动的耦合，且对环向波数较小的模态，频率计算并不准确。

本文考虑由转动引起的科氏力、离心惯性力及环向初应力影响，利用 Hamilton原理，建立了基于Sanders壳体理论的转动薄壁圆柱壳振动微分方振 动 工 程 学 报 第 26卷程 。选取满足相应边界条件的轴向梁函数来逼近各类边界条件下的圆柱壳轴向振型分布 ▲而采用 Galerkin方法，对振动微分方程进行离散，得到 6自由度陀螺转子系统的动力学方程。然后，通过叠加各行波模态的响应，得到了任意典型边界条件下转动薄壁圆柱壳行波振动响应-般解的形式。最后，分别针对静坐标系下横向简谐力和恒力作用下的两端固支转动薄壁圆柱壳的行波振动响应，给出了相应算例，并对结果进行了简要分析。

1 基于 Sanders壳体理论的转动薄壁圆柱壳振动微分方程考虑长 L、厚 H、半径为 R，并以角速度 n绕中轴转动的各向同性薄壁圆柱壳，该圆柱壳的密度、弹性模量、泊松比分别为f0，E， 。如图 1所示，正交曲线坐标系(z，0，z)为建立在壳体中曲面上的动坐标系。相对于该坐标系，壳体中曲面上点的轴向、切向、法向的变形位移分别为 ， ，叫。在壳体中面微元上，-z，0，z三个方向的单位面积外载荷分量分别q ，q ，q 。

图 1 载荷作用下的转动薄壁圆柱壳模型图Fig.1 Sketch of the geometrical relations of a thin rota-ting cylindrical shel转动薄壁圆柱壳的动能可表示为，- r[碗z。z曲 22(wo～劬)厶n ( W )]RdxdO (1)由环向初应力所引起的初应变能为1l铷N (去 oqU)。 1(Or -b叫)fl(ow ) 胁∞式中 N；-0Hn R。为离心惯性力导致的环向初应力 。

由鼻振动理论可知，壳体上任-点的应变与中曲面的薄膜应变分量 e ，e ，o甜以及弯曲应变(曲率变化)分量 ， ，r。有关，具体表达形式如下 ㈣f ] 5。j J- 三 0 0 aIzo 1丽3 3 oRa az孙。 -筹a a0R a R a3a 2a2Raz R8x80(4) l j兰- 兰 csQ 。 ， Ql2- ， G Qzz z -万上 - 王王 1 - l ，u2 Jf-H j2， [ e r瑚e -q ，- ，q ) ，系数 -R (1- )/EH，算子 L-2 o o -第 3期 孙述鹏 ，等：转动薄壁圆柱壳行波振动响应分析其中，L L1-R。象 2( 鲁)导 筹，L2R12/1- ] - R ，L -R 2( )象c 著 著，L -皤LL- 尼 。 杀])， 嘉，这里忌-H。/(12R。)，y -P尺 (1- )/E，下同。算子L 形式如下r R ：著。ROO7-y。 。R 3dZ - 2y。 昙dt实际工程系统不可避免存在阻尼因素，受阻尼影响，系统瞬态响应将衰减，进而剩余工程中关心的稳态响应。通常情况下，阻尼是由试验测定的，这里假设阻尼为等效粘性阻尼，比例系数为 c ，并在建立的振动方程 中引入阻尼项。振动控制方程 (10)可改写为如下形式L L l- q c 疵 c 0 - c 由 (13)2 转动薄壁圆柱壳行波振动响应-般解对转动薄壁圆柱壳体 自由振动特性的研究表明[193，转动引起 的离心惯性力使壳体产生周向应力，导致壳体刚度增加，从而使得转动薄壁圆柱壳较静态圆柱壳频率有所增加；同时，受转动薄壁圆柱壳变形速度矢量与转动速度矢量方向不同所引起的科氏力影响，转动薄壁圆柱壳频率随转速变化发生分岔，产生频率不同的前、后行波。此时的模态已不是通常意义下的振型(驻波)，而是与时间有关的行波 。

基于此，将转动薄壁圆柱壳行波振动解的-般形式设为如下形式 ： - ∑∑[7 (t)cosnO- ( )sinnO]9(x)m 1 1- ∑∑[ (t)sinn8 。cosn0]q)；(z)仇 ； 1 1W - ∑∑[7 。(t)COsnO- 。(t)sinnO]9：(x)m - 1 n 1(14)式中 m， 分别为轴向振型的阶数和环向的波数。

u(z)， (z)和 (z)分别为 3个方向上的轴向振型分布，其形式与边界条件有关。

考虑到除两端简支约束薄壁圆柱壳外，具有其他边界条件的圆柱壳轴向振型分布很难用解析式表示出来，这里利用圆柱壳振型的轴向分布接近于相应边界条件梁振型函数 的特性 ，采用梁函数来逼近圆柱壳轴向振型函数。-般可设 3个方向的振型分布有如下形式口。](z)-d(x)/dx， (z)- ( )，w(z)- (z) (15)式中 (z)为相应边界条件下连续梁的振型函数，- 般形式如下(z)C1 sin(x) f2COS(z) f3 sinh(z)C4 cosh(x) (16)系数 C ( 1，，4)由边界条件决定。为后面算例中分析方便，这里给出两端固支连续梁的振型函数，其他边界条件的梁的振型函数可以参考振动理论的相关内容 。

(z)-COS( z/L)-cosh(a z/L)rim[sin(A z/L)-sinh(a x/L)] (17)其 中f -- cos -cosh ]/[sin A -sinh ]lCOS ·cosh --10(18)将解的形式(14)代入转动薄壁圆柱壳振动微分方程(13)，并采用 Galerkin方法进行离散，对每个m， 组合的模态可得到-个6自由度系统动力学方程 ，写成矩阵形式有 /#H 0 0 ]J 0 c /pH 0 l，l 0 0 c /p/j(19)，K-[K ]，( ，J-1，2，3)，叫 -n-a-巩2 -旦- l-C， . . . 。 . . . . 、 、，，....。

、/ ) ) l 1 3 -G) ( 2 1， j f/ O 易 O f n0m驵 f/O O C n2 -462 振 动 工 程 学 报 第 26卷-1[R2案( 。卅 )]，112- [- ] Cm，～K13- 1 -tR ]景，121- [ny2gYR - ]鲁，- 1[( )豢(kn2 ]，- [c 鲁～- 1[/R- ]鲁，- [( )C鲁m 卅K3 [(1尼 nzyz -y n )-zkR2n z鲁 cjAi]， 等愚R4 fjQ1 -lQ q1g0q0qqq6- fL f 2pHn J。J。

q (z， ，t)COSn u (z)/Clmlq (z，0， )sinnO ( ) 21q (z， ，t)COSn (z)/c-q (z，0，t)sinnO u(z)/cq (z，0，t)COSnO (z) 21-q (z，0，t)sinnO w(z)/cdOdx，f -j。L u(z) u( )d山r，Cm 2-J ( ( )) u(z)dx，c -J ( ( )) (z)dx，c -j ( ( )) u (z) ，瞄 -J L0 v(z) v(x)dx，Cm 2-J Lo w(z) v(x)dx，fm 。-j ( ( ))Cm( )出，嵋-J ( (z)) (z)妇，cm25-j ，v (z)) v(z)(ix，c -j L0 w(z) W(x)dx，cm32-j Lo v(z) (x)dx，Cm 3-J ( (z)) (z)dx，c器- ( ( )) (x)dx，Cm -f ( (z)) (z) ，c3%-J ( (z)) w( )dx。

对方程(19)进行数值积分，可得每个 m，n组合的模态所对应的广义坐标的时间历程。理论上讲，将各个 m，n组合的模态所对应的广义坐标的时间历程，运用式(14)进行叠加，即可得出转动薄壁圆柱壳上各点的行波振动响应。

值得注意的是，环向波数 n-0时，转动圆柱壳的各阶模态均表现为驻波的形式，而考虑到本文主要探讨的是转动圆柱壳的行波振动响应的求解方法，因此在式(14)中没有考虑 -0这-特殊模态。

3 数值结果与讨论根据前面给出的方法，针对具有两端固支约束的转动薄壁圆柱壳模型，进行了计算和分析，算例所用几何参数及材料常数列于表 1〖虑到壳体的面外弯曲振动是对壳体动力学特性起主导作用的振动形式，故这里只分析对应于弯曲振动的行波频率特性，和法向的位移响应 W。

表 1 几何参数及材料常数Tab.1 Geometric parameter and material property长 千任 鼙厚 悍任俣亘 曼L/m R/m Him E/(G·Pa)p/(kg·m )2 1 0.002 70 2 6003.1 行波振动频率特性与模型验证忽略阻尼，则对每个m， 组合的模态所对应的6自由度系统动力学方程 (19)可退化为陀螺转子系统动力学方程的形式罡(20)因此，求解转动薄壁圆柱壳行波频率的问题，就转化为求解陀螺系统特征值问题。Meirovitch给出了求解此类特征值问题的方法 。需要指出，对于每个 m， 组合的模态，有 6个频率值，其中频率最低的两个对应于 W方向弯曲振动为主的行波频率。

为验证离散模型的正确性，本文将所求得的行波频率与由文献[6]中的解析表达式所得结果进行了比较 ，两者结果吻合 ，如表 2所示 。

图2所示为转速 6 000 r/min，轴向振型阶数- 1情况下，动坐标下中行波频率随环向波数的变化曲线。动坐标系下，后行波频率大于前行波频率，且前、后行波频率随环向波数的增加均呈现先减小后增加的趋势，并在环向波数大于 9时逐渐重合。

图 3所示为轴向振型阶数 m-1，环向波数 -3时，动坐标系中行波频率随转速的变化曲线。如图 3所示，转速为零时，前后行波频率值-致，产生的是驻波，即传统意义的振型。随着转速的增加，频率发生第 3期 孙述鹏，等 ：转动薄壁圆柱壳行波振动响应分析分岔，产生频率不同的前、后行波。前行波频率随转速升高先有微弱的减小趋势，继而随转速的升高而增加。而后行波频率随转速的变化，则-直呈现单调增加的趋势 。

表 2 本文所得行波频率与文献值的比较(转速6 000 r/rain)Tab.2 Comparison of travelling wave frequency between presentand literature study(rotating speed 6 000 r/min由文献[6]公式(12)求得 图 2 行波频率与环向波数的关系曲线Fig.2 Variation of the frequency parameter with respecttO the circumferential wave number图 3 行波频率随转速变化曲线Fig.3 Variation of the frequency parameter with respecttO rotating speed(Z1)式中 z和0分别为初始时刻激励点轴向和环向坐标 。

为计算和分析方便起见，这里仅取 m-1， -3所对应的行波模态，并根据式(14)进行叠加。

图4给出了初始时刻zL/2， -0位置，受静坐标系下幅值为 F-2 000 N，频率 fo-100 Hz简谐激励作用时，转速 -6 000 r/rain(100 Hz)的两端固支薄壁圆柱壳 32 -L/2，0 -兀/2位置行波振动响应的时间历程、稳态时间历程、稳态相图及频谱图。如图 4所示 ，受阻尼的影响 ，瞬态 响应逐 渐衰减，系统响应达到稳态。受转动影响，静坐标系下的单频激励 fo-100 Hz作用下的两端固支圆柱壳的稳态响应有两个频率成分，即 200，400 Hz，分别对应 于 392-. ， 392。

《0. 《O 2 4 6 8 10tS(a)响应时间历程(a)Response time historycurves9.95 9.96 9.979.98 9.99 10.00t/S(b)稳态响应时间历程Co)Steady-state time historyCBrves位移/mm(c)稳态相图(c)Steady-state phase diagram200 400 600 800f/Hz(d)稳态频谱图(d)Steady-state frequencyspectrogram图 4 静坐标中横向简谐力作用下(m-1， -3)组合行波模态的响应Fig.4 Response of(m- 1， - 3)traveling mode for arotating cylindrical shell under transverse har-monic loads in the stationary coordinate3.2 的两端 固2小q'标1 誊 l' -& 熹 动薄壁圆柱壳的行波振动响应H 用 时 ，嘉 ” 6 0 0川r/mi (o10v464 振 动 工 程 学 报 第 26卷图 5 静坐标系 中横 向简谐力作 用下 (m- 1， -3)组合行波模态响应的幅频曲线Fig.5 Amplitude-frequency curve of(m- 1， - 3)trav-elling mode for a rotating cylindrical shell undertransverse harmonic loads in the stationary coot-dinate置行波振动 响应 的幅频 曲线。由图 5中的标示可知，该系统前、后行波频率分别为 f，-294.4 Hz和db-410.0 Hz〖虑动、静坐标下行波频率关系，有 ，-f，3n-594.4 Hz，f -f -30-110.0Hz。 o 和-厂 分别为前、后行波共振频率，这与图 5幅频曲线所示峰值频率-致 。

综合上面的分析可知，-般地 ，当静坐标系下横向简谐力 的频率 满足如下条件时，会发生行波共振 ，-厂 -f，， ha2， fo-f ， -ha2 (22)具有转动薄壁圆柱壳结构的实际工程系统，在工作时应尽量避开这些共振点。

3.3 静坐标系中法向恒 力作用下两端固支转动薄壁圆柱壳的行波振动响应令 fo-0，由式(21)可得，受静坐标下幅值为 F的法向恒力作用的转动薄壁圆柱壳，其所受单位面积外载荷分量的表达式fq (z，0， )-0q z， ， -。 - -lq ( ，0， )--F/(R·L) (z-z) [ -( - )](23)若( ，0 )表示初始时刻动坐标系中转动薄壁圆柱壳上-点 P，此时该点位于静坐标下的空间点P。，两点重合。 时刻后 ，与静坐 标下的空 间点 P。

重合的转动薄壁圆柱壳上点的坐标则可表示为(z ，0 -n )。由式(14)，相应点的位移响应为叫 - ∑∑[ 。(f)cosn(O - )-m - 1 n- 11723(t)sinn(0 - )] (z ) (24)图 6给出了与静止坐标下 4个空间点位置重合的两端固支薄壁圆柱壳上点的行波振动响应时间历程图。圆柱壳 转动 速度 为 -6 000 r/min(100Hz)，且受静坐标系下幅值为F-2 000 N的法向恒力作用，初始时刻该恒力作用于z-L/2，0-0〖虑到该圆柱壳长径比较小，低阶行波模态的轴向振型阶数 聊通常取 1，并参照图 2给出的行波频率与环向波数的关系曲线，计算中相应点的位移响应采用了如下近似表达式1 10- ∑∑[ 。(t)cosn(O -g)-- l - 1叩23( )sinn(0 -g]t)] ( ) (25)O。

t/Sfa)x'L/2，口 OO.O- 0.4- O.8- 12- 1.6- 2.00 2 4 6 8 10t/sCo)x'L/2.07r焉 --- 2.0 . . . . . -2.0 . . . . 。

t/S(c)x'L/2，0z/2t/s(d) /2，0 3z/2图6 与静坐标下空间点对应的转动薄壁圆柱壳上相应点的响应Fig.6 Response of one point on a rotating cylindricalshell in stationary coordinate如图 6所示，4个空间点所对应的转动薄壁圆柱壳上点的响应随时间变化逐渐收敛。即，静坐标系下法向恒力作用下的两端固支转动薄壁圆柱壳，其形变形状在达到稳态后不随时间而改变 ，如 图 7所示。从波动的观点讲，两端固支转动薄壁圆柱壳，受静坐标系下法 向恒力作用 ，产生 了静坐标 系下 的驻波。受此影响，当形变导致 的应力水平超过材料破坏极限 ，可能产生破坏，此外也可能会给转动薄壁图7 静坐标系下法向恒力作用下的两端固支转动薄壁圆柱壳的形变Fig.7 Deformation of a rotating cylindrical shell underpoint forced嗣 呈第 3期 孙述鹏，等：转动薄壁圆柱壳行波振动响应分析 465圆柱壳带来疲劳加速等问题，这些都是值得进-步研究的课题。

4 结 论本文考虑由转动引起的科氏力、离心惯性力及环向初应力影响，利用 Hamilton原理，建立了基于Sanders壳体理论的转动薄壁圆柱壳振动微分方程，提出了-种适用于求解各种边界约束的转动薄壁圆柱壳行波振动响应的方法。基于此方法，分别针对静坐标系下横向简谐力和恒力作用下的两端固支转动薄壁圆柱壳的行波振动响应进行了求解和分析。得到以下结论：(1)通过选取满足相应边界条件的轴向梁函数来逼近各类边界条件下的圆柱壳轴向振型分布，可以克服除简支以外的各类边界条件下圆柱壳振型的解析表达式不易获取的困难。在此基础上发展的方法，适用于求解各种类型边界约束的转动薄壁圆柱壳行波振动响应。

(2)两端固支转动薄壁圆柱壳，随着转速的增加，其频率发生分岔，产生频率不同的前、后行波。

前行波频率随转速升高先有微弱的减小趋势，继而随转速的升高而增加。而后行波频率随转速的变化 ，则-直呈现单调增加的趋势。

(3)当静坐标系下横向简谐力的频率 满足fo-fs 72 或fo-fb-ng2时，两端固支转动薄壁圆柱壳会发生行波共振，具有转动薄壁圆柱壳结构的实际工程系统，在工作时应尽量避开这些共振点。

(4)受静坐标系下法向恒力作用，两端固支转动薄壁圆柱壳会产生不随时间改变的形变，即静坐标系下的驻波。受此影响，当形变导致的应力水平超过材料破坏极限，可能产生破坏，此外也可能会给转动薄壁圆柱壳带来疲劳加速等问题，这些都是值得进-步研究的课题 。