非对称矩阵特征值问题密集模态重分析方法

本文提出了一种非对称矩阵特征值问题的密集模态重分析方法。它将原密集特征值问题表达为与其临近的某一重特征值问题的小摄动,从而密集模态的重分析问题就转化为重频模态的重分析问题。     

任意非亏损矩阵特征灵敏度分析的模态展开法

把特征向量的各阶导数表示成所有模态的线性组合，并利用左模态与右模态间的双正交性，首先导出了任意非亏损矩阵的重特征值的一阶导数所满足的特征值问题，然后根据此特征值问题无、看重根的情况，再导出了异导重特征值和等导重特征值对应的可微特征向量、特征值和特征向量各阶导数的一般计算公式。算例显示了方法的正确性。     

一类结构的（n－3）重根及其非退化性

基于陀螺模态综合法，从总装阵的块式结构出发，构造性地证明了具有ｎ片桨叶的旋翼型结构陀螺特征值问题存在一系列的（ｎ－３）重特征根，得到了对应的（ｎ－３）个完备的振型．结论进而推广到有阻尼的旋翼型结构．继续研究证明过程表明：结论适用于更广泛的一类具有重复子结构的结构系统，结果表明这类结构几何上的重复性或对称性导致的重根不会引入退化性．不同类型的算例验证了所得到的解析结果．本文还试图说明动力子结构法的定性性质保持特性是值得继续探讨的课题     

Sensitivity Analysis of Multiple Eigenvalues*

ABSTRACT

This paper is devoted to sensitivity analysis of eigenvalues of nonsym-metric operators that depend on parameters. Special attention is given to the case of multiple eigenvalues. Due to the nondifferentiability (in the common sense) of multiple roots, directional derivatives of eigenvalues and eigenvectors in parametric space are obtained. Sensitivity analysis is based on the perturbation method of eigenvalues and eigenvectors. The generalized eigenvalue problem and vibrational systems are also investigated. Strong and weak interaction of eigenvalues are distinguished and interactions in two- and three-dimensional space are treated geometrically. It is shown that the strong interaction of eigenvalues is a typical catastrophe. Simple examples that illustrate the main ideas are presented. The results obtained are important for qualitative and quantitative study of mechanical systems subjected to static and dynamic instability phenomena.     

智能结构密频系统振动控制及其摄动分析

提出一种智能结构密频控制规律的设计方法。其主要工作为：将密频子空间转化为重频子空间,把密频系统作为重频系统上的摄动;为了解决重频密频子空间相对应的特征向量选取的敏感性问题,通过摄动分析得到与重频密频子空间相对应特征向量的线性组合并利用闭环系统特征值极点配置的方法得到重频密频系统的振动规律;把所设计的振动控制规律作用于原系统和摄动系统上,讨论了结构参数改变后对系统的动态特性的影响;最后通过算例证明该方法的有效性。     

亏损系统模态优化控制Riccati方程的求解方法

讨论亏损重复特征值的Ｒｉｃｃａｔｉ方程的模态控制算法．由于亏损特征值的Ｊｏｒｄａｎ块矩阵的阶数ｍ小于状态方程的阶数ｎ,即ｍ＜＜ｎ;所提出的模态控制算法可极大地减少计算工作量．其中的数值例子说明了该方法的有效性．     

亏损系统振动模态的可控可观性

利用奇异值分解方法来讨论系统广义模态的可控可观性的量度问题,得到了亏损系统广义模态可控可观性的量度指标,同时用实例说明了本文方法是有效的。     

任意非亏损系统特征灵敏度分析的直接摄动法

本文发展了一种任意非亏损系统特征灵敏度分析的二阶段摄动法，将未摄动的问题的解作为零阶近似，把摄动影响作为摄动后问题的高阶修正，经过严格的数学推导，得到了支配高阶修正量的完全方程组。本方法无损知道摄问题的全部特征向量，仅需知被摄模态的特征对。本方法可处理被摄问题具有重特征值，甚至具有等导重特征值这一高度退化极难处理的情况。算例显示了本方法的正确性。     

接近亏损系统的矩阵摄动法

随着系统参数的变化，带有重频的系统可转变成带有密集频率的系统，反之亦然．本文讨论了亏损系统与接近亏损系统之间的关系．并提出了接近亏损系统的平均移位的摄动方法．算例表明了此方法的有效性．     

一般实矩阵特征值问题的多维灵敏度分析

应用Kronecker代数和矩阵微分理论系统发展了一般实矩阵特征值问题的矩阵灵敏度分析方法，结构系统特征值和特征向量的导数很方便地排列成二维矩阵，所得的结果很容易形成计算机程序。     

Second-order sensitivity of eigenpairs in multiple parameter structures

This paper presents methods for computing a second-order sensitivity matrix and the Hessian matrix of eigenvalues and eigenvectors of multiple parameter structures. Second-order perturbations of eigenvalues and eigenvectors are transformed into multiple parameter forms,and the second-order perturbation sensitivity matrices of eigenvalues and eigenvectors are developed.With these formulations,the efficient methods based on the second-order Taylor expansion and second-order perturbation are obtained to estimate changes of eigenvalues and eigenvectors when the design parameters are changed. The presented method avoids direct differential operation,and thus reduces difficulty for computing the second-order sensitivity matrices of eigenpairs.A numerical example is given to demonstrate application and accuracy of the proposed method.     

ON THE STABILITY BOUNDARY OF HAMILTONIAN SYSTEMS

ＩｎｔｒｏｄｕｃｔｉｏｎＭａｎｙｍｅｃｈａｎｉｃａｌｓｙｓｔｅｍｓｃａｎｂｅｖｉｅｗｅｄａｓｌｉｎｅａｒＨａｍｉｌｔｏｎｉａｎｓｙｓｔｅｍｓ.Ｗｈｅｎｗｅｓｔｕｄｙｔｈｅｅｆｆｅｃｔｓｏｆｓｙｓｔｅｍｐａｒａｍｅｔｅｒｓｏｎｔｈｅｂｅｈａｖｉｏｒｏｆｔｈｅｓｙｓｔｅｍｓ,ｔｈｅｓｙｓｔｅｍｃａｎｂｅｒｅｇａｒｄｅｄａｓｔｈｅｓｙｓｔｅｍｄｅｐｅｎｄｉｎｇｏｎｐａｒａｍｅｔｅｒｓ.Ｖｅｒｙｉｍｐｏｒｔａｎｔｓｙｓｔｅｍｐａｒａｍｅｔｅｒｓ ,ｓｕｃｈａｓｃｒｉｔｉｃａｌｌｏａｄ ,ｃｒｉｔｉｃａｌａｎｇ…     

Eigenvalue-counting methods for non-proportionally damped systems

This paper presents numerical methods of counting the number of eigenvalues for non-proportionally damped system in some interested regions on the complex plane. Most of the eigenvalue analysis methods for proportionally damped systems use the well-known Sturm sequence property to check the missed eigenvalues when only a set of the lowest modes is used. However, in the case of the non-proportionally damped systems such as the soil–structure interaction system, the structural control system and composite structures, no counterpart of the Sturm sequence property for undamped systems has been established yet. In this study, a numerical method based on argument principle is explained with emphasis on the discretization of the contour and a new method based on Gleyse’s theorem is proposed. To verify the applicability of the methods, two numerical examples are considered.     

A new matrix perturbation method for analytical solution of the complex modal eigenvalue problem of viscously damped linear vibration systems

A new matrix perturbation analysis method is presented for efficient approximatesolution of the complex modal quadratic generalized eigenvalue problem of viscouslydamped linear vibration systems.First,the damping matrix is decomposed into the sum of aproportional-and a nonproportional-damping parts,and the solutions of the real modaleigenproblem with the proportional dampings are determined,which are a set of initialapproximate solutions of the complex modal eigenproblem.Second,by taking thenonproportional-damping part as a small modification to the proportional one and using thematrix perturbation analysis method,a set of approximate solutions of the complex modaleigenvalue problem can be obtained analytically.The result is quite simple.The new methodis applicable to the systems with viscous dampings-which do not deviate far away from theproportional-damping case.It is particularly important that the solution technique be alsoeffective to the systems with heavy,but not over,dampings.The solution formul     

Stochastic model reduction for robust dynamical characterization of structures with random parameters

In this paper, we characterize random eigenspaces with a non-intrusive method based on the decoupling of random eigenvalues from their corresponding random eigenvectors. This method allows us to estimate the first statistical moments of the random eigenvalues of the system with a reduced number of deterministic finite element computations. The originality of this work is to adapt the method used to estimate each random eigenvalue depending on a global accuracy requirement. This allows us to ensure a minimal computational cost. The stochastic model of the structure is thus reduced by exploiting specific properties of random eigenvectors associated with the random eigenfrequencies being sought. An indicator with no additional computation cost is proposed to identify when the method needs to be enhanced. Finally, a simple three-beam frame and an industrial structure illustrate the proposed approach.     

An acceleration for the eigensystem realization algorithm with partial singular values decomposition

The real-time identification of dynamic parameters is important for the control system of spacecraft. The eigensystem realization algorithm (ERA) is currently the typical method for such application. In order to identify the dynamic parameter of spacecraft rapidly and accurately, an accelerated ERA with a partial singular values decomposition (PSVD) algorithm is presented. In the PSVD, the Hankel matrix is reduced to dual diagonal form first, and then transformed into a tridiagonal matrix. The eigenvalues are computed by the bisection method in terms of the Sturm property, and the corresponding eigenvectors are obtained by the inverse iteration method. Finally, the eigenvalues and the eigenvectors are transformed into the singular values and the singular value vectors of the original matrix. An example for space station is presented to demonstrate the efficacy and accuracy of the proposed algorithm.     

Eigenvalue bounds for linear mechanical systems with nonmodal damping

Eigenvalue bounds for damped linear mechanical systems are presented based on a ‘modified equation of state’. The system exhibits nonmodal damping, by which we mean that classical normal modes do not obtain. The special case of an underdamped system is discussed, leading to bounds on a generalized damping ratio. The bounds are expressed in terms of the extreme eigenvalues of the constituent matrices. Geometrically, they define rectangular regions in the complex plane which contain the system eigenvalues.     

广义特征值问题中重特征值的特征向量导数

本文把重特征值的特征向量导数的计算方法，推广到非亏损矩阵的广义特征值问题，并给出了特征值导数也有重根时特征向量导数的计算式。本方法的优点是只需已知所考虑的重特征值的特征向量，因而计算量小，对于大型复杂结构更为适用。     

广义特征值问题中重特征值的特征向量导数

本文把重特征值的特征向量导数的计算方法，推广到非亏损矩阵的广义特征值问题，并给出了特征值导数也有重根时特征向量导数的计算式。本方法的优点是只需已知所考虑的重特征值的特征向量，因而计算量小，对于大型复杂结构更为适用。     

二类复特征值表达式的适用性研究

对于线性系统，根据其质量、刚度或阻尼等矩阵的性质了解特征值，是十分有意义的问题．对于有阻尼陀螺系统，本文首先通过模态正交性推导了有关矩阵的若干关系式．然后分析了复特征值的两类表达式的适用性，着重讨论了通过一元二次方程研究特征值性质的现行处理方法存在的困难．提出一种确定特征值的补充方法，指出可能存在虚数特征值的系统的性质．算例表明结论正确