陀螺系统辛子空间迭代法.pdf

第19卷第1期 振 动 工 程 学 报 V01．19No．1 2006年3月 ofVibration Mar．2006 Journal Engineering 陀螺系统辛子空间迭代法 隋永枫，钟万勰 (大连理工大学工程力学系，辽宁大连116023) 摘要：转子系统的有限元分析可以导出陀螺系统的本征值问题，而陀螺本征值问题可在哈密顿体系下求解。基于辛 子空间迭代法的思想，提出了一种求解陀螺系统本征值问题的算法。首先引入对偶变量，将陀螺动力系统导入哈密 顿体系，将问题化为了哈密顿矩阵的本征值问题，由于稳定的陀螺系统其本征值必为纯虚数，利用这个特点，提出 了对应陀螺系统的辛子空间迭代法，从而可以求出系统任意阶的本征值及其振型。算例证明了这种算法的有效性。 关键词：陀螺系统；哈密顿矩阵；辛子空间；本征值；转子 中图分类号：0318文献标识码：A 文章编号：1004—4523(2006)01—0128—05 引 言 守陀螺系统的理论做了进一步的发展，在文献[5] 中，通过一组变换将陀螺系统本征值问题化为一对 陀螺系统特征值问题一直是转子动力学提出的 称矩阵的本征值问题，并给出了离散陀螺系统的 典型数学问题。经过多年的努力，发展出一系列的快 Rayleigh商定理，以及陀螺系统的最大最小定理和 速有效的求解方法[1]。但这些算法通常都是基于系 分离定理，完善了陀螺本征值的理论基础并给出了 统的刚度矩阵正定时的情况，但在许多高转速的情 一种有效的离散陀螺系统本征值的数值计算方法。 况下，经常会出现系统刚度矩阵不正定的情况[引，特 别是在自由度较多的情况下，其特征值问题是很难 础和算法，但是它始终没有克服由于计算过程中对 进行求解的。本文基于此种情况提出了陀螺系统的 矩阵求逆所带来的影响计算效率和计算误差的缺 辛子空间迭代法，这种方法较好地解决了当系统刚 点。而且需要指出，这种方法无法求解不正定陀螺系 度矩阵不正定(即对应哈密顿函数不正定)时陀螺系 统的本征值问题。 统的本征值问题。以下介绍陀螺系统特征值问题的 张文，包华丽在1989年提出了一种求解正定陀 发展概况。 螺系统本征值问题的高效算法[6]，这种方法先通过 早在1971年，Paardekooper首先提出了反对 平面镜象变换，把反对称阵化为缩减后的三对角形 称阵的Jacobi方法[3]，通过二维旋转，把反对称矩阵 式的反对称阵，再利用这个特殊形式的反对称阵， 化为典则型。这种方法虽然可以证明是收敛的，但其 把问题化为三对角正定对称阵的普通特征值问题， 精度较低，特别是只能无选择地算出全部的特征值， 然后对其进行求解，解决了正定陀螺系统本征值问 所以这种方法主要适用于低阶反对称矩阵的求解。 题。 1974年，Meirovitch提出了把陀螺特征值问题 另外在求解陀螺系统本征值问题方面比较常用 转换成对称阵的特征值问题的方法n]，这种方