热环境中超临界黏弹性输流管道自由振动分析.pdf

第期宿恒等:大变形柔性管道两相流流致振动研究文章编号:()D O I:/收到第稿,收到修改稿国家自然科学基金青年基金项目(),国家自然科学基金杰出青年基金(),上海市教育委员会科研创新计划项目(E ),上海 市 科 学 技 术 委 员 会 浦 江 人 才 计 划(A类 P J ),N a t i o n a lN a t u r a lS c i e n c eF o u n d a t i o no fC h i n a(),t h eN a t i o n a l S c i e n c eF u n df o rD i s t i n g u i s h e dY o u n gS c h o l a r s(),t h eP r o g r a mo fS h a n g h a iM u n i c i p a lE d u c a t i o nC o m m i s s i o n(E ),a n dt h eP u j i a n gP r o j e c to fS h a n g h a i S c i e n c ea n dT e c h n o l o g yC o mm i s s i o n(P J )通信作者E m a i l:x m a o s h u e d u c n热环境中超临界黏弹性输流管道自由振动分析高思禹毛晓晔丁虎陈立群(上海大学 力学与工程科学学院,上海 )摘要基于E u l e r B e r n o u l l i梁模型,本文研究了热环境中输流管道在超临界范围内流固耦合自由振动特性考虑温度增量以及初始轴向拉力作用,在两端简支边界条件下,利用广义H a m i l t o n原理建立输流管道横向振动偏微分积分控制方程通过解析方法得到输流管道非平凡静平衡位形及临界流速精确表达式,与微分求积单元法(D Q EM)数值结果吻合较好基于复模态法,结合伽辽金(G a l e r k i n)法离散系统偏微分积分控制方程,得到热环境下超临界输流管道的模态函数和固有频率结果表明,温度增量越大,临界流速越小,此时的管道越容易屈曲,但相同流速下超临界管道固有频率越大;初始拉力越大,临界流速越大,相同流速下超临界固有频率越小该研究可以为热环境中超临界状态下的管道系统振动设计提供理论指导关键词热环境,输流管道,超临界流速,非平凡静平衡位形中图分类号:O 文献标志码:AF r e eV i b r a t i o nA n a l y s i so nS u p e r c r i t i c a lV i s c o e l a s t i cP i p e sC o n v e y i n gF l u i d i nT h e r m a lE n v i r o n m e n tG a oS i y u M a oX i a o y eD i n gH u C h e nL i q u n(S c h o o l o fM e c h a n i c sa n dE n g i n e e r i n gS c i e n c e,S h a n g h a iU n i v e r s i t y,S h a n g h a i ,C h i n a)A b s t r a c t B a s e do n t h eE u l e r B e r n o u l l i b e a mt h e o r y,t h e c u r r e n tw o r k i n v e s t i g a t e s t h e f r e ev i b r a t i o no f as u p e r c r i t i c a lp i p e c o n v e y i n g f l u i d i n t h e t h e r m a l e n v i r o n m e n t C o n s i d e r i n g t h e t e m p e r a t u r e i n c r e m e n t a n dt h e i n i t i a l a x i a l f o r c e,t h ep a r t i a l d i f f e r e n t i a l i n t e g r a lg o v e r n i n ge q u a t i o nf o rt h et r a n s v e r s ev i b r a t i o no ft h es i m p l ys u p p o r t e dp i p e i s e s t a b l i s h e db y t h eg e n e r a l i z e dH a m i l t o nsp r i n c i p l e T h en o n t r i v i a l e q u i l i b r i u mc o n f i g u r a t i o n sa n dc o r r e s p o n d i n gc r i t i c a l f l o wv e l o c i t i e sa r eo b t a i n e da n a l y t i c a l l y,a n dv e r i f i e db yt h ed i f f e r e n t i a l q u a d r a t u r e e l e m e n tm e t h o d(D Q EM)B a s e do n t h e c o m p l e xm o d em e t h o da n d t h eG a l e r k i nm e t h o d,t h en a t u r a l f r e q u e n c i e so f t h es u p e r c r i t i c a l p i p ec o n v e y i n gf l u i d i nt h e t h e r m a l e n v i r o n m e n ta r eo b t a i n e d T h e r e s u l t s s h o wt h a t t h e l a r g e r t h e t e m p e r a t u r e i n c r e m e n t i s,t h e s m a l l e r t h e c r i t i c a l f l o wv e l o c i t yw i l lb e,w h i c ha l s om a k e s t h ep i p eb u c k l e dm o r ee a s i l y B u t t h e s u p e r c r i t i c a l n a t u r a l f r e q u e n c i e sb e c o m eb i g g e ra t t h es a m e f l u i dv e l o c i t y T h ed i s c u s s i o na l s of i n d st h a t t h e l a r g e r t h e i n i t i a l t e n s i o ni s,t h e l a r g e r t h e c r i t i c a l f l o wv e l o c i t yw i l l b e U n d e r t h i s c o n d i t i o n,t h en a t u r a l f r e q u e n c i e s i n t h e s u p e r c r i t i c a l r e g i o nb e c o m es m a l l e ra t t h es a m e f l o wv e l o c i t y T h i ss t u d yp r o v i d e sg u i d a n c e f o r t h ev i b r a t i o nd e s i g no fp i p e sc o n v e y i n gf l u i d i nt h es u p e r c r i t i c a l r e g i o nu n d e r t h e t h e r m a l e n v i r o n m e n t K e yw o r d s t h e r m a l e n v i r o n m e n t,p i p ec o n v e y i n gf l u i d,s u p e r c r i t i c a l f l u i ds p e e d,n o n t r i v i a le q u i l i b r i u mc o n f i g u r a t i o n动力学与控制学报 年第 卷引言输流管道流致振动的研究有广泛的应用前景,其研究成果可直接运用于海洋和土木工程、石油和化学加工工业、核电和电力工业、船舶和飞机,以及日常生活中管道内部流固耦合作用和外部激励引起的振动,会不断地传递到整个系统,引发系统大幅振动,从而导致管道出现疲劳损伤,甚至在极端振动条件下会因管道结构破裂,导致液体泄漏而产生不可估量的后果当管道中流体流速超过一定值时,管道将发生屈曲失稳,并在新的平衡位形附近振 动,这 被 称 为 输 流 管 道 系 统 的 超 临 界 状态,其特征为管道基频归零环境温度会显著改变系统共振特性,因此,将环境温度考虑进输流管道的研究是具有现实意义的;环境温度条件下的输流管道振动特性分析也是对已有理论的补充和发展,具有一定学术价值目前已有的研究主要聚焦于输流管道的振动特征,这些研究包含输流管道振动特性对管道寿命影响、流固耦合作用、建立精确的输流管道动力学模型及管道大幅度振动时的非线性动力学等杨晓东和金基铎发现因为假设的模态函数会使G a l e k i n截断法得到的固有频率产生误差,于是通过复模态法与其比较,结果表明截断阶数较低时,两种方法得到的固有频率拟合较好D i n g等针对非线性连续系统,建立了与两个非线性隔振器耦合的输流管道非线性受迫振动动力学模型周永兆等人研究了两端铰支输流管道的横向非线性自由振动,并且对偏微分方程做准静态假设得到偏微分积分方程,研 究 表 明 偏 微 分 控 制 方 程 有 较 强 的 非 线性颜雄和魏莎等研究了两端非对称弹性支承输流管道横向振动的固有特性,发现流体的流速越大,其对两端非对称弹性支承的管道固有频率的影响越明显当流速超过临界值时,管道会发生静态屈曲,即管道发生屈曲失稳,系统会重新稳定在曲线平衡位形D i n g等分别使用伽辽金截断法计算了轴向运动梁的固有频率且通过数值方法计算了轴向运动梁在超临界运动状态下的平衡位移 随岁寒和李成利用虚功原理建立了输流管道弯曲和振动的有限元方程,给出了两种边界条件下管道自由振动的前三阶固有频率与流体流速的关系,结果表明,前三阶固有频率随流速增大而减小 Z h o u等人研究了四种不同初始形状的悬臂式输送流体管道的非线性静力和动力特性结果证明了管道由流动引起的管道静态变形也可能非常大,管道颤振失稳的临界流速主要取决于刚刚达到临界流速时的静力平衡构型 缪旭和金基铎针对超临界状态下的输流管道的扰动方程,通过伽辽金法得到了超临界状态下系统非线性振动第一阶固有频率的解析表达式 M a o等研究了输流管道内共振为时的受迫振动响应,研究发现当流体的速度超过临界值时,直线平衡构型会发生弯曲,原来的单稳态系统转化为双稳态系统 T a n等人首次用T i m o s h e n k o梁理论研究了超临界范围内输流管道的振动特性,并且与E u l e r B e r n o u l l i管道相比,T i m o s h e n k o管道的平衡构型曲率更大,两种管道模型之间的相对差异受流速的影响 L u等人研究了输流管道在超临界绕流曲线处的内部共振和应力分布,以改善管道的机械疲劳性能 Y e等研究了超临界范围内微弯曲输送流体管道的振动特性,研究发现弯曲管道的固有频率可能会随着管道长度的增加而增加,但不一定是单调的 尽管近年来超临界高速输流管道的振动问题一直受到广泛关注,但目前的研究工作主要集中在理想环境的基础上在实际工程应用中,环境温度的改变也会诱发管道发生振动,降低管道寿命根据文献,Q i a n等人通过数值方法分析了温度和流速对输流管道的固有特性的影响结果表明,动力学特性由温度和流速共同决定,温度升高也会导致管道刚度的降低,即使内部没有流体也会导致管道发散不稳定 E b r a h i m i M a m a g h a n i等人研究了考虑热效应的轴向瑞利管的稳定性,发现温度越高时管道越容易以发散和颤振的形式失去稳定性 B a b a e i分析了热环境中功能梯度微曲管道在均匀侧压力作用下的非线性响应,结果发现温度升高提高了微曲管道的强度,同时也提高了微曲管道的上限载荷和下限载荷 但是以上工作没有涉及对热环境中输流管道超临界振动的研究基于上述文献调研,本文建立了具有线性热弹性关系的E u l e r B e r n o u l l i输流管道在超临界状态下的运动控制方程;对管道热环境下随流速分岔的位形进行了数值验证;利用复模态法和伽辽金法得到管道的模态函数并分析温度增量和初始轴力对第期高思禹等:热环境中超临界黏弹性输流管道自由振动分析管道超临界固有频率的影响运动方程及非平凡静平衡位形图为两端简支输流管道的物理模型简图管道材料密度为p,长度为L,外径为D,内径为d,P为初始轴力,弹性模量E,截面绕中性轴转动惯量Ip考虑K e l v i n黏弹性本构关系,黏弹性系数为管道内流体假设为不可压的黏性流体,密度为f,流速恒定管道均布Fc o s(t)的简谐横向激励力,只考虑激励平面内的横向位移v(x,t),其中t代表时间坐标,x表示沿管道轴向的空间坐标此外,考虑系统处于均匀温度场中,以室温为基础,温度增量为T图输流管道物理模型F i g P h y s i c a lm o d e l o f t h ep i p ec o n v e y i n gf l u i d根据能量法建立运动方程,E u l e r B e r n o u l l i输流管道横向运动的动能Te为TeLpApv,tdxLfAfv,ts i n(v,x)dxLfAfc o s(v,x)dx()式中第一项表示输流管道横向振动的动能,第二项及第三项表示流体的动能系统内力及外力做功W为WV(xT)xdV(UpUf)()式中V为输液管道的体积,x(x,t)和x(x,t)分别为正应力和正应变,T表示热应力,Up表示轴力P做功,Uf表示激励力做功,相应的表达式UpLPv,xvdx,UfLFc o s(t)vdx()式中管道材料应力和应变本构关系可表示为xE x x,t()当温度场发生变化时,管道会产生应变,沿x轴的热应变可表示为 TT()其中,是热膨胀系数因此,热应力可以表示为TE T()假设由振动变形引起的附加应力沿管道的轴向变化很小,可以将应力取作沿输流管道在长度上的平均值x(t)LLv,xdx()根据广义H a m i l t o n原理ttTedttt Wdt()从而得到热环境中黏弹性横向振动输流管道的偏微分积分模型pApv,t tfAf(v,t t v,x tv,x x)(PE ApT)v,x xE Ipv,x x x x Ipv,x x x x tE ApLv,x xLv,xdxApLv,x xLv,xv,x tdxApLv,x x tLv,xdxFc o s(t)()以及两端简支条件v(,t)v,x x(,t),v(L,t)v,x x(L,t)()表给出了本文研究的输流管道的物理参数值及几何参数值表输流管道物理参数T a b l eP h y s i c a l p a r a m e t e r so f t h ep i p ec o n v e y i n gf l u i dI t e mN o t a t i o nV a l u e管道长度(m)L管道外径(m)D 管道内径(m)d 弹性模量(G P a)E 初始拉力(N)P 管道材料密度(k g/m)p 流体密度(k g/m)f 热环境中输流管道非平凡静平衡位形非平凡位形解析解为求解输流管道静态屈曲构型,忽略式()中显含时间的项,得到只与空间坐标有关的位移函数(fAfPE ApT)v,x xE Ipv,x x x x动力学与控制学报 年第 卷E ApLv,x xLv,xdx()式()中,v即表示非平凡静平衡位形为求得解析解,先将定积分项假设为常数,该方程可以化简为(fAfPE ApT)v,x xE Ipv,x x x xE ApLv,x x()按照求导阶数,方程()可以改写为v,x x x xv,x x()其中fAfPE ApTE ApL/E Ip,方程()是一个四阶常微分方程,通解为vCCxCc o s(x)Cs i n(x)()其中Ci(i,)是常数将边界条件()代入常微分方程()得到特征方程s i n或k,k,即可得到输流管的非平凡静平衡位形 vCs i n(kx)()将表达式()代入式()可得到屈曲位形系数C kLfAfPE ApT()kE IpE Ap()从式()可以导出,只有当流速c r(k)(kE IpPE ApT)/fAf时,产生一对对称平衡曲线位形,c r(k)表示第k阶临界流速由于超临界流速具有复杂多变性,输流管道分析难度加大,因此,重点关注超临界流速范围在c r()c r()下,研究热环境中简支输流管道围绕正平衡位形v(x)的横向振动的非线性动力学特性非平凡位形仿真解微分 单 元 求 积 法(D Q EM)在 微 分 求 积 法(D QM)的基础上,将管道两端的旋转角度作为广义变量引入因此,它有N个节点基础上的(N)个变量基于(N)个内点方程和个边界条件,这(N)个变量可以被求解出来因此,D Q EM可以准确地满足四阶边界值问题由于篇幅所限,D Q EM对偏微分方程插值离散公式见 已有文献表明当节点采用非均匀分布时,计算结果精度及收敛性更高,因此采用C h e b y s h e v G a u s s L o b a t t o非均匀分布,即xic o s(i)NL,i,N()本文选择 个节点(N)进行模拟,以获得更准确的结果(a)基于D Q EM法的衰减振动曲线(a)D a m p e df r e ev i b r a t i o nv i aD Q EM(b)非平凡位形解析解与数值解对比验证(b)C o m p a r i s o no f t h ea n a l y t i c a ln o n t r i v i a l e q u i l i b r i u mc o n f i g u r a t i o nt oD Q EM图热环境中输流管道非平凡位形F i g N o n t r i v i a l e q u i l i b r i u mc o n f i g u r a t i o no f t h ep i p ei nt h e r m a l e n v i r o n m e n t图所示为热环境中输流管道的非平凡静平衡位形数值解和解析解的对比其中图(a)表示D Q EM中(N)个内部节点的自由衰减情况当它们趋于稳定时,输流管道静止于势能最小位置,该位置即屈曲平衡位形表格表示各节点对应的坐标值图(b)点表示微分求积单元法计算得到的数值解,实线表示理论推导得到的解析解从图中可以看出解析结果和数值结果吻合得很好,从而证明了非平凡静平衡位形解析解的正确性图为不同温度增量下非平凡静平衡位形图,结果表明温度增第期高思禹等:热环境中超临界黏弹性输流管道自由振动分析量越大,非平凡静平衡位形越大,管道越容易屈曲图不同温度增量下的非平凡静平衡位形F i g N o n t r i v i a l e q u i l i b r i u mc o n f i g u r a t i o n sw i t hd i f f e r e n t t e m p e r a t u r e i n c r e m e n t s热环境中超临界管道固有频率将管 道 的 非 平 凡 静 平 衡 位 形 代 入 式(),v(x,t)v(x,t)v(x),得到超临界输流管道非平凡平衡位形基础上的振动控制方程(pApfAf)v,t tfAf v,x t(fAfv,x xPE ApT)v,x xE Ipv,x x x xE ApLv,x xLv,xdxE ApLv,x xLv,xv,xdx()方程()的解可以写作v(x,t)n(x)eint()式中n为系统第n阶固有频率,n(x)为第n阶模态函数将式()代入局部线性化方程()可得(pApfAf)n(x)n ifAf nn(x)(fAfPE ApT)n(x)E Ip()n(x)E ApLn(x)Lv,xdxE ApLv,x xLn(x)v,xdx()选用静态梁模态函数为基函数,假设n(x)MjDn js i n(jx)()其中,Dn j为待定系数,形式解()满足简支边界条件()将模态函数()代入方程(),在方程两边同乘s i n(ix),i,M,并在,L 上对x进行积分,可得到关于n和dn j(j,M)的超静定齐次线性代数方程组 nMinGK dn dn MT()其中,dn jDn j/Dn(Dn),M为质量矩阵,G为陀螺矩阵,K为刚度矩阵式()存在非零解的条件为方程组的系数行列式为零|nMinGK|()基于此,可以得到不同离散阶数系统的固有频率图描述了用不同G a l e r k i n截断阶数计算得到的前两阶固有频率随超临界流速变化的趋势,从图中可以看出,随着流速的增大,超临界固有频率也增大从误差分析可以发现在四阶和六阶截断情况下,系统第一阶固有频率曲线拟合良好,第二阶固有频率的误差不大于因此认为四阶G a l e r k i n截断 达 到 了 收 敛,后 续 研 究 都 将 采 用 四 阶G a l e r k i n截断法对控制方程进行截断降维处理其中,截断阶数误差表达式为RGn(G a l e r k i n)n(G a l e r k i n)n(G a l e r k i n)()图超临界输流管道前两阶固有频率随流速的变化及截断误差分析F i g A c c u r a c yo f t h e f i r s t t w on a t u r a l f r e q u e n c i e sc h a n g i n gw i t hf l u i dv e l o c i t y i nt h es u p e r c r i t i c a l r e g i o n参数对热环境中固有频率的影响温度增量对超临界固有频率的影响基于收敛性结论,图对比了在室温基础下,不同温度增量对超临界输流管道固有频率的影响根据非平凡静平衡位形(),当环境温度升高时,管道的屈曲位形增大由于引入温度增量,刚度矩阵增大,此时的弯曲刚度增强文献,研究表动力学与控制学报 年第 卷(a)前两阶超临界固有频率随温度增量变化趋势(a)T h e f i r s t t w on a t u r e f r e q u e n c i e sc h a n g i n gw i t ht h e t e m p e r a t u r ei n c r e m e n t i nt h es u p e r c r i t i c a l r e g i o n(b)不同温度增量下第一阶超临界固有频率(b)T h e f i r s tn a t u r a l f r e q u e n c yw i t hd i f f e r e n t t e m p e r a t u r ei n c r e m e n t i nt h es u p e r c r i t i c a l r e g i o n(c)不同温度增量下第二阶超临界固有频率(c)T h es e c o n dn a t u r a l f r e q u e n c yw i t hd i f f e r e n t t e m p e r a t u r ei n c r e m e n t i nt h es u p e r c r i t i c a l r e g i o n图温度增量对超临界固有频率的影响F i g E f f e c to f t h e t e m p e r a t u r e i n c r e m e n t t on a t u r a lf r e q u e n c i e s i nt h es u p e r c r i t i c a l r e g i o n明,亚临界状态下简支输流管道在温度升高时,由于热膨胀导致管道产生轴向压力,从而温度升高会使管道固有频率降低,临界流速随温度升高而降低图(a)显 示 了 在 给 定 的 超 临 界 流 速(c r)下,管道固有频率随温度变化的增大而增大,与亚临界情况相反图(b)表明,超临界状态下流速越大,固有频率越大,与亚临界情况相反;温度增量越大,输流管道的临界流速越小,与亚临界结果一致图(c)表明当流速较大时,温度对管道固有频率的影响较小且温度升高对管道第二阶固有频率影响不明显轴力对热环境中超临界固有频率的影响图为初始轴力分别为 N、N、N时,热环境中超临界输流管道前两阶固有频率随流速的变化图(a)显示了在给定的超临界流速(c r)下,管道固有频率随轴力的增大而减小图(b)、(c)的结果表明,初始轴力越大,输流管道的临界流速越大也可得,在超临界情况下,温度增量与初始轴力对输流管道固有频率的影响相反(a)前两阶超临界固有频率随初始轴力变化趋势(a)T h e f i r s t t w on a t u r e f r e q u e n c i e sc h a n g i n gw i t ht h e i n i t i a la x i a l f o r c e i nt h es u p e r c r i t i c a l r e g i o n(b)不同初始轴力下第一阶超临界固有频率(b)T h e f i r s tn a t u r a l f r e q u e n c yw i t hd i f f e r e n t i n i t i a l a x i a lf o r c e i nt h es u p e r c r i t i c a l r e g i o n第期高思禹等:热环境中超临界黏弹性输流管道自由振动分析(c)不同初始轴力下第二阶超临界固有频率(c)T h es e c o n dn a t u r a l f r e q u e n c yw i t hd i f f e r e n t i n i t i a la x i a l f o r c e i nt h es u p e r c r i t i c a l r e g i o n图初始轴力对超临界固有频率的影响F i g E f f e c to f t h e i n i t i a l a x i a l f o r c e t on a t u r a l f r e q u e n c i e si nt h es u p e r c r i t i c a l r e g i o n结论本文以热环境中超临界自由振动的输流管道为研究对象,通过哈密顿原理建立了两端简支受温度影响的输流管道的控制方程得到非平凡静平衡位形及临界流速,基于D Q EM方法,进行数值验证利用复模态法进行模态分析,将模态函数作为权函数,使用四阶伽辽金截断法得到流固耦合作用下管道的固有频率及不同超临界流速下的模态函数,并分析了温度增量及初始轴力对管道固有频率的影响经过研究发现,非平凡静平衡位形近似解析解与数值解拟合良好,验证了其正确性当管道微幅振动时,超临界固有频率和模态函数不受非线性系数影响,随着超临界流速增加,高阶固有频率增大有趣的是,在亚临界状态下,温度增量越大或轴力越小时输流管道的固有频率越低,临界流速越低超临界状态下,温度增量越大或初始轴力越小,固有频率反而越大,此时临界流速越小,且温度增量和初始轴力的改变对管道第二阶固有频率影响较小此外,亚临界状态下,温度增量和流速的增大都会导致挠度增大,而轴力的增加会抑制最大挠度与本文温度增量、流速及轴力对静平衡位形的影响一致,验证了本文结果的正确性参考文献L IS,KA R N E YBW,L I UG F S I r e s e a r c h i np i p e l i n es y s t e m s Ar e v i e wo ft h el i t e r a t u r eJ J o u r n a lo fF l u i d sa n dS t r u c t u r e s,:GUO Y,Z HU B,YANG B,e ta l F l o w i n d u c e db u c k l i n ga n dp o s t b u c k l i n gv i b r a t i o nc h a r a c t e r i s t i c so fc o m p o s i t ep i p e si nt h e r m a le n v i r o n m e n tJ O c e a nE n g i n e e r i n g,:张凯凯,谭霞,丁虎,等超临界输流管道内共振下参激振动响应 J应用数学和力学,():Z HAN GKK,TANX,D I N G H,e t a l P a r a m e t r i cv i b r a t i o nr e s p o n s e so fs u p e r c r i t i c a lf l u i dc o n v e y i n gp i p e s i n:i n t e r n a l r e s o n a n c eJ A p p l i e dM e m a t i c sa n dM e c h a n i c s,():(i nC h i n e s e)丁虎输流管道振动的被动控制研究 J自然杂志,():D I NG H P a s s i v e c o n t r o l o f v i b r a t i o n i np i p e s c o n v e y i n gf l u i dJ C h i n e s eJ o u r n a lo fN a t u r e,():(i nC h i n e s e)杨晓东,金基铎输流管道流固耦合振动的固有频率分析 J振动与冲击,():YAN GXD,J I NJD A n a l y s i so f t h ei n h e r e n t f r e q u e n c yo f f l o w s o l i dc o u p l i n gv i b r a t i o n i n f l o wt r a n s f e rp i p e sJ J o u r n a l o fV i b r a t i o na n dS h o c k,():(i nC h i n e s e)D I NG H,J I JC,CHE NLQ N o n l i n e a r v i b r a t i o n i s o l a t i o nf o rf l u i d c o n v e y i n gp i p e su s i n gq u a s i z e r os t i f f n e s sc h a r a c t e r i s t i c sJ M e c h a n i c a l S y s t e m s a n dS i g n a lP r o c e s s i n g,:周永兆,杨晓东,金基铎输流管道非线性横向振动固有频率 分 析 J振 动、测 试 与 诊 断,(S):,Z HOU YZ,YANGXD,J I NJD N a t u r a l f r e q u e n c ya n a l y s i so fn o n l i n e a rt r a n s v e r s ev i b r a t i o n sf o rp i p e sc o n v e y i n gf l u i dJ J o u r n a lo f V i b r a t i o n,M e a s u r e m e n t&D i a g n o s i s,(S):,(i nC h i n e s e)颜雄,魏莎,毛晓晔,等两端弹性支承输流管道固有特性研究 J力学学报,():YANX,WE I S,MAOXY,e t a l S t u d yo nn a t u r a lc h a r a c t e r s t i c so ff l u i d c o n v e y i n gp i p e s w i t he l a s t i cs u p p p o r t sa t b o t he n d sJC h i n e s e J o u r n a l o fT h e o r e t i c a l a n dA p p l i e d M e c h a n i c s,():(i nC h i n e s e)黄慧春,张艳雷,陈立群受迫振动的超临界输液管G a l e r k i n数值模拟 J应用数学和力学,():动力学与控制学报 年第 卷HUANG H C,Z HAN G Y L,CHE N L QAG a l e r k i nn u m e r i c a lm e t h o df o rt h ep i p ec o n v e y i n gs u p e r c r i t i c a l f l u i du n d e r f o r c e dv i b r a t i o nJ A p p l i e dM e m a t i c s a n dM e c h a n i c s,():(i nC h i n e s e)D I N G H,CHE NLQ G a l e r k i nm e t h o d s f o rn a t u r a lf r e q u e n c i e so fh i g h s p e e da x i a l l ym o v i n gb e a m sJJ o u r n a lo fS o u n da n d V i b r a t i o n,():D I N G H,Z HAN GGC,CHE NLQ S u p e r c r i t i c a le q u i l i b r i u ms o l u t i o n so fa x i a l l ym o v i n gb e a m sw i t hh y b r i db o u n d a r yc o n d i t i o n sJ M e c h a n i c sR e s e a r c hC o mm u n i c a t i o n s,():张国策,丁虎,陈立群复模态分析超临界轴向运动梁横向非线性振动 J动力学与控制学报,():Z HAN GGC,D I N GH,CHE NLQ C o m p l e xm o d a la n a l y s i so ft r a n s v e r s a l l yn o n l i n e a rv i b r a t i o nf os u p e r c r i t