管道振动问题中一种辛变换的计算方法.pdfVIP

2009年4月 沈阳航空工业学院学报 Apr．2009 第26卷第2期 JournalofShenyangInstituteofAeronauticalEngineering Vo1．26 No．2 文章编号：1007—1385(2009)02—0005—04 管道振动问题 中一种辛变换的计算方法 张艳雷 金基铎 马 斌 (沈阳航空工业学院航空宇航工程学院，辽宁 沈阳 110136) 摘 要：以符合哈密顿 (Hamiltonian)特性的管道振动方程的简化问题给出了一种具体的计算方 法。为了保持原方程哈密顿结构，这种简化过程中使用的矩阵必须是辛矩阵(sympletiematrix)。 针对管道振动方程给出一种辛矩阵的计算方法，该方法利用了哈密顿方程线性部分的特征空间， 用特征空间的一些性质来完成辛矩阵的计算。 关键词：哈密顿特性；辛矩阵；特征空间 中图分类号：TB533 文献标识码 ：A 一 个哈密顿系统的方程 ，通过适当的变换往 面；：(u一)c控+A：+ 6 往可以得到实质性的简化。已有的管道振动方程 模型_l 大多是在拉各朗日体系下建立的，并不 一 0 0 ．，=( )c=(为靴 3 符合哈密顿特性。用哈密顿力学方法建立的管道 式中： 为质量 比，g为广义位移，P 为广义动 振动方程 J，其线性部分保持有哈密顿结构 (JHJ 量，7为管道预警力，‰为流速。 = H )。众所周知，在保持哈密顿结构前提下的 对于两端铰支管道 ： 简化必须是通过辛变换才能完成。辛变换是指符 b=8／3，Cll= 一1T，C22= 一4ax，Al ，A2=2,tro 合辛条件 ( J。T= )的矩阵变换，这种变换能够 对于两端固支管道 ： 使简化后的方程仍保持其哈密顿结构。由辛条件 b=3．3420，C1l= 一12．3026，22= 一46．0501，Al= 我们也可以看到，这种辛矩阵并不是唯一的。本 文针对输流管道振动方程 ，在临界流速条件下给 4．73o0，A2=7．8532。 令 ： 出了具体的辛矩阵计算方法。该方法具有一定的 M，“0b l 0 通用性 。 O 0 1 1 管道振动方程 的哈密顿结构及临 H = O 0 M，“0b 界情况 一 “，1 一 M，0b 0 首先分析和计算临界流速u，并在临界流速