系统与其子系统的稳定性及相互关系 - 中国力学学会.docVIP

系统与其子系统 张卫 暨南大学 理工学院 应用力学研究所 Email：tzhangw@jnu.edu.cn 摘要： 本文讨论复杂动力学系统及其子系统运动稳定性的关系问题。整体与部分是自然界的基本属性。本文通过一种系统的分解和集结方法，将系统处理成整体系统及其子系统两级结构进行稳定性分析，从中我们不仅能得到系统及其子系统的稳定性信息，还能获得子系统与系统间的稳定性关系的信息。显然，本文介绍的系统与其子系统稳定性关系的分析方法具有较普遍的应用领域。 关键词：稳定性，大系统，分解与集结，Liapunov 0．引言 复杂动力学系统的运动稳定性，特别是系统与其子系统的稳定性关系问题一直受理论研究者及工程师们的关注[1-3]，如对大型汽轮发电机组轴系的运动稳定性及多体动力学问题等。上世纪中叶发展起来的“大规模动力系统稳定性理论”为我们提供了理论分析基础。特别在国际知名学者A.N.Michel，R.K.Miller， D.D.Siljak等人的努力下, 这个理论得到了极大的发展[4][5]。这一理论要解决的核心问题，是如何从子系统的稳定性信息获得大系统的信息。我国学者高为炳、刘永清等人也已为大系统稳定性理论做出了突出贡献[6,7]。上世纪80年代，M.Ikeda和D.D.Siljak提出称为“重叠分解－集结”的分析方法，解决了某些大系统经典分析方法过于保守的问题，使一些“强耦合”问题得到了解决［8-10］，但其思路仍放在由子系统信息获取大系统信息的思路上。大系统理论及分析方法仍在发展中，但一般认为，至上世纪90年代，这种经典大系统稳定性问题基本趋于稳定发展阶段。 另一方面，虽然经典大系统分析的思路是通过获取子系统的信息，达到获得大系统信息的目的，但许多科学及工程中的问题，更需要我们了解子系统与大系统的关系及某子系统对大系统影响程度的信息。这样的复杂动力学系统我们可以列举出如控制系统、生化系统、大型转子－油膜轴承系统、多体动力学系统、生物系统、甚至经济及管理大系统等。针对稳定性问题，本文发展了一种解决上述问题的方法，即将复杂动力学系统视为两级：系统及其子系统，并通过集结得到一扩张系统。本文将显示，基于这一理论而发展出的分析方法，使我们不仅能得到子系统和其整体系统的稳定性信息，而且能够获得子系统与其大系统稳定性的关系，子系统对大系统稳定性影响的程度及子系统的变化对大系统稳定性敏感度的信息，解决上面提出的工程及科学实践中急需解决的问题。 [作者简介] 张卫（1959 -），男，研究员（博士），博导，研究方向：非线性动力学 非线性系统的原点稳定性 一个非线性动力学系统可以表示成如下方程： (1) 其中：为系统广义坐标，t指时间变量，为广义惯量矩阵，为陀螺力及循环力列向量， 为广义激励力列向量。 系统（1）的原点稳定性可由其一次近似系统 (2) 的原点去近似。这里，为对称矩阵，为反对称矩阵，不失一般性，上式可简写为： 其中，称为广义位移向量，且 ，分别称为广义质量、阻尼和测度矩阵。 则由Liapunov稳定性理论，系统（1）的原点稳定性与系统 (3) 的原点稳定性相同［11］，将上式改写为状态方程的形式： (4) 其中，，N=2n， 而系统矩阵 ， (5) 其中，I为单位矩阵。 则线性化系统（4）的稳定性可由构造二次型的Liapunov函数进行判断： (6) 其中，对称正定矩阵H满足Liapunov矩阵方程： (7) 复杂系统的二级结构 对一个实际的复杂动力学系统，子系统的识别一般可由问题本身给出，但一般要求子系统应具有式（4）、（5）所表达的完备动力学行为。譬如，多体系统的每一体的动力学关系，汽轮发电机组轴系中每一个油膜轴承及相应转轴等，可视为子系统。理论上，子系统的鉴列可根据问题的需要确定，对上述两类具备“树”或“链”式结构的系统，其系统矩阵A具有形式： (8) 其中 表示第i个子系统 ， (9) 将（8）写成 (10) 其中： ，，， (11) 　　　　现在，要构造一个系统，在这个系统中，原系统及各可识别的子系统皆包括其中，从而使原系统被处理成二级结构．为此，设一新状态变量