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

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

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

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

一种用于非线性振动系统的模态分析方法

本文提出了一种用于非线性振动系统的模态分析方法,将求解非线性系统模态的问题化为求解非线性特征值、特征向量的问题,利用模态研究系统的响应,文中分析了非线性保守系统、非线性自治系统和非线性非自治系统的线性模态,导出了三个模态包含原理。     

工程结构的随机特征问题研究及其在梁结构中的应用

采用子结构模态综合和摄动随机有限元相结合求解工程结构的随机特征问题。为求出随机特征对的方差，借助于模态截断概念推出诸特征值与特征向量对随机变量的偏导数。以染结构为典型算法，定量研究了子结构动模态的选取个数与随机特征对的计算精度间关系，以梁的长细比首次确定使用Timoshenko梁和Euler-Bernoulli梁两模型求解梁类工程结构随机特征问题的适用范围。     

复振型导数计算的移位多次模态加速法

本文研究了含粘阻尼结构的复振型导数计算问题，将导数计算问题成是一个简谐激谐的振的响应计算问题，采用多次模态加速法和移位法，导出了复振型导数计算的移位多次模态加速法，该方法具有确的数学和物理意义，可导出了已有的各种计算方法。算例表明本方法计算复振型导数只需用很少几个模态即可保证精度，计算量大大减少。     

复振型导数计算的移位多次模态加速法

本文研究了含粘性阻尼结构的复振型导数计算问题，将导数计算问题看成是一个简谐激振的响应计算问题，采用多次模态加速法和移位法，导出了复振型导数计算的移位多次模态加速法。该方法具有明确的数学和物理意义，可导出已有的各种计算方法。算例表明本方法计算复振型导数只需用很少几个模态即可保证精度，计算量大大减少。     

考虑粘弹性材料频变特性的复合结构频率响应分析

论文采用基于Layerwise离散层理论的有限元板单元建模粘弹性阻尼复合结构,发展了一种考虑粘弹性材料的频变特性的复合结构频率响应计算方法.该方法首先计算结构的质量矩阵和各个频率点下的结构复刚度矩阵,然后求解运动方程的若干阶复特征对并由此计算各阶模态的传递函数,最后将各模态传函线性叠加得到近似的总传函.为了提高计算效率,论文采用了一种高效的数值方法,即只计算若干频率点下的特征向量与特征值,并计算这些点处特征向量关于频率的高阶导数,通过泰勒系数展开逼近和Rayleigh商式,可求得附近若干频率点处的特征向量和特征值,从而避免了在各个频率点下求解大自由度结构特征方程的问题,可以极大地提高计算效率.对一端简支的三层约束阻尼梁算例进行了分析,并与文献中的结果作对比,结果验证了方法的有效性和计算效率.     

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

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

动柔度与特征向量导数

现有特征向量导数分析了“直接法”用于求解很多特征向量导数的问题时,显得计算效率过低,为解决这个问题,本文导出了一种过渡动柔度式,使得计算很多个特征向量的导数时,计算效率高于现有各种“直接法”。     

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

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

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.     

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

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

关于含参变量特征值轨迹偏转现象的探讨

含参变量的特征值轨迹的偏转现象出现在许多力学和理论物理问题中。本文通过两个实例,对特征值轨迹偏转现象进行了深入的探讨。先讨论了一个两维平板振动问题,通过求特征值的上、下界的方法,确定了该问题的精确特征值呈现了曲线(含一个参变量)偏转和曲面(含两个参变量)偏转现象;然后讨论了一个双震荡器的振动问题,发现当某个参数发生变化时的特征值曲线是由相交变为偏转的。且当一对特征值曲线偏转前后,两个特征向量几乎互相改变了方向;而当一对特征值曲面偏转前后,对应的特征向量的分量也发生偏转。     

Generalized-order perturbation with explicit coefficient for damage detection of modular beam

A general method is formulated to estimate damage location and extent from the explicit perturbation terms in specific set of eigenvectors and eigenvalues. At first, perturbed orthonormal equation is generated from the perturbation of eigenvectors and eigenvalues to obtain the k-th explicit perturbation coefficients. At second, perturbed eigenvalue equation is generated from the perturbation of eigenvector and eigenvalue, and first-order expansion of the stiffness matrix to obtain other explicit perturbation coefficients. Stiffness parameters are computed from these equations using an optimization method. The algorithm is iterative and terminates under certain criteria. A fixed–fixed modular beam with various numbers of elements is used as test structure to investigate the applicability of the developed approach. By comparison with the Euler–Bernoulli beam, discretization errors are analyzed. In six elements beam, first-order algorithm converges faster for small percentage damage. Second-order algorithm is more efficient for medium percentage damage. For large percentage damage, the second-order algorithm converges more effectively. Meanwhile, for eight elements large percentage damage and ten elements small percentage damage, second-order algorithm converges faster to the termination criterion.     

A spectral approach for the stability analysis of turbulent open-channel flows over granular beds

A novel Orr–Sommerfeld-like equation for gravity-driven turbulent open-channel flows over a granular erodible bed is here derived, and the linear stability analysis is developed. The whole spectrum of eigenvalues and eigenvectors of the complete generalized eigenvalue problem is computed and analyzed. The fourth-order eigenvalue problem presents singular non-polynomial coefficients with non-homogenous Robin-type boundary conditions that involve first and second derivatives. Furthermore, the Exner condition is imposed at an internal point. We propose a numerical discretization of spectral type based on a single-domain Galerkin scheme. In order to manage the presence of singular coefficients, some properties of Jacobi polynomials have been carefully blended with numerical integration of Gauss–Legendre type. The results show a positive agreement with the classical experimental data and allow one to relate the different types of instability to such parameters as the Froude number, wavenumber, and the roughness scale. The eigenfunctions allow two types of boundary layers to be distinguished, scaling, respectively, with the roughness height and the saltation layer for the bedload sediment transport.     

A cantilever conveying fluid: coherent modes versus beam modes

A semi-analytical approach to obtain the proper orthogonal modes is described for the non-linear oscillation of a cantilevered pipe conveying fluid. Theoretically, while the spatial coherent structures are the eigenfunctions of the time-averaged spatial autocorrelation functions, it emerges that once the Galerkin projection of the proper orthogonal modes is carried out using the uniform cantilever-beam modes, the spatial dependency of the integral eigenvalue problem can be eliminated by analytical manipulation which avoids any spatial discretization error. As the solution of the integral equation is obtained semi-analytically by linearly projecting the proper orthogonal modes on the cantilever-beam modes, any linear or non-linear operation can be carried out analytically on the proper orthogonal modes. Furthermore, the reduced-order eigenvalue problem minimizes the numerical pollution which leads to spurious eigenvectors, as may arise in the case of a large-scale eigenvalue problem (without the Galerkin projection of the eigenvectors on the cantilever-beam modes). This methodology can conveniently be used to study the convergence of the numerically calculated proper orthogonal modes obtained from the full-scale eigenvalue problem.     

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

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