大型非经典阻尼系统的动态解耦方法

提出了大型非经典阻尼系统动力分析的动态解耦方法，同时在特征值问题的求解和动力响应的计算中完全避免复数运算．方法适用于无阻尼陀螺系统、非经典阻尼系统（包括阻尼陀螺系统），且不论系统亏损与否．得到了揭示阻尼作用的有趣的理论结果和实用的数值方法     

多自由度系统复模态理论的摄动方法——(一)一阶摄动

除了阻尼矩阵满足一定条件外,有阻尼多自由度线性系统运动方程,在一般情况下不能通过实模态变换而解耦。因此,许多情况下工程结构动力分析需要寻求系统的复模态和复特征值,为此如Foss.Frazy and Bishop等提出的惯用方法又太复杂和不经济。本文采用基于实模态理论的摄动方法,耒求解系统的复模态和复特征值,考虑到阻尼力比惯性力和弹性恢复力要小是符合工程实际的,把系统的模态和特征值按不同的量级展成级数,从而建立起各阶渐近方程,其零阶方程对应于无阻尼系统可按实模态理论求解,如果需要,可按高阶方程逐次求解得到复模态和复特征值各阶渐近修正。本方法不仅计算方便而且经济,其结果易于从零阶和一阶近似中得到复模态和复特征值,对于自由振动运动方程同样可以解耦。利用已得到的一阶复模态的结果,讨论了自由振动和强迫振动问题。文末给出了算例以说明本方法的计算精度。     

粘弹性阻尼修改后结构系统动特性的预测

无阻尼或比例阻尼结构系统经粘弹阻尼修改后,可变为粘弹性阻尼系统。要获得其动特性,需求解复特征值问题。但是,随之带来了计算量大、费用高等问题。尤其是粘弹性材料特性随频率而变化,需求解高阶非线性复特征值问题,这对于一个自由度较大的结构,计算量太大,通常较难实现。本文在特征值修改方法的基础上,提出粘弹性阻尼局部动力修改方法,即仅需已知原结构系统的实模态参数,就可求出粘弹阻尼修改后系统的复模态参数。还发展了特征值和修改量同时迭代方法,有效地解决了粘弹材料复模量随频率变化引起的非线性复特征值问题。     

求一般结构系统部分复模态的Ritz向量法及应用

构造了一种获取一般结构系统（刚度，阻尼矩阵可以是非对称的）Ｒｉｔｚ向量基的迭代方法，从而可以高维阻尼系统部分复模态解降价成为一个低维二次特征问题，将此降阶后的二次特征问题化成标准特征值后即可用ＱＲ方法求得系统的部分低阶持征解，作为应用，对一个实际转子－轴承系统的动力特性进行了计算分析。     

多自由度参激系统稳定性的数值分析

提出多自由度周期参激系统稳定性的数值直接法。通过将扰动方程表示成状态方程形式,再根据Flo-quet理论将扰动解表示成指数特征分量与周期分量之积,并将其周期分量与系统周期系数展成Fourier级数,导出一系列代数方程,建立矩阵特征值问题,从而由数值求解特征值可直接确定参激系统的稳定性。该方法可用于一般周期参激阻尼系统,特征值矩阵不含逆子阵。应用于斜拉索在支座周期运动激励下的参激振动不稳定性分析,数值结果表明该方法的有效性。     

关于非保守力学系统的稳定性

关于非保守力学系统的稳定性李俊峰王照林李铁成（清华大学工程力学系，北京１０００８４）１引言力学系统分为保守系统和非保守系统．对保守系统稳定性的研究开展得比较早，研究成果也比较多．非保守系统存在阻尼，线性阻尼有瑞利阻尼和约束阻尼．对于仅有瑞利阻尼作用的...     

随机结构系统的一般实矩阵特征值问题的概率分析

由于工程实际结构的复杂性和所用材料在统计上的离散性以及测量、加工、制造误差的存在，必然导致具有随机参数的随机结构振动系统，按结构参数的性质来划分，随机振动问题包括两方面内容：(1)确定结构问题；(2)随机结构问题。本文以现代数学理论为依托，研究了随机结构系统的一般实矩阵的特征值问题。根据Kronecker代数、向量值和矩阵值函数的灵敏度分析、一般二阶矩法和概率摄动技术给出了计算随机结构系统的一般实矩阵的特征值和特征向量的数值方法，可以有效地得出随机结构系统的一般实矩阵的特征向量的统计量，发展了2D矩阵值函数的随机结构系统的特征值问题概率分析理论。     

具有重频特性系统的动力响应分析

一般有阻尼线性系统出现重特征值时,基于振型正交性的复振型分解法将不再适用.本文综合运用高等数学、线性代数和复变函数理论,对具有重频特性的一般有阻尼线性多自由度系统给出了系统动力响应在时域中的计算方法.该方法充分利用复振型分解法和留数矩阵解耦法的优点,不仅概念清晰,而且易于理解和掌握,适合于大型复杂系统的动力响应分析.此外,本文给出了双自由度体系产生重特征值的条件,对典型实例进行了地震响应分析,并通过与Newmark-β法计算结果的对比,论证了文中所给计算公式的正确性.本文提出的分析方法具有普适性,对线性结构、机电和控制系统也都是适用的.     

动力学平衡方程的Euler中点辛差分求解格式

给出了动力学方程${\pmb M}\ddot {\pmb x} + {\pmb C}\dot {\pmb x} + {\pmb K \pmb x} = {\pmb R}$的二阶Euler中点隐式差分求解格式，分保守系统、无 阻尼受迫振动系统和阻尼系统3种情况, 讨论了算法中Jacobi矩阵${\pmb A}$的性质，譬 如${\pmb A}$是否为辛矩阵以及谱半径等. 对于无阻尼系统，证明了无论是否存在外 载荷，Jacobi 矩阵都是辛矩阵. 证明了辛矩阵的所有本征值的模为1，其谱半径永远 为1, 以及$\delta = 0.5$和$\alpha = 0.25$的Newmark算法就是Euler中点隐式差 分格式，对保守系统它们都是辛算法. 严格证 明了Euler中点辛格式是严格保持系统能量的. 通过算例详细讨论了保辛算法用于求解非保 守系统动态特性的优越性，如广义保结构特性等；分析了保辛算法的相位误差以及由其引起 的系统的附加能量特性；分析了保辛算法和$\delta \ne 0.5$的Newmark算法的精度随着激励频率与系统固有频率比的变化情况等     

不确定参数闭环控制系统特征值的区间分析

运用区间理论,讨论了区间参数的结构振动控制问题,给出了求解闭环系统区间特征值的一种方法．基于区间参数导出了区间刚度矩阵和质量矩阵,然后利用矩阵摄动理论和区间扩张理论,推导了复区间特征值上下界估计的算法．这些结果是从二阶系统的左右特征向量出发得到的．将该文方法应用到悬臂梁的控制问题,数值结果表明它是有效的．     

First-order perturbation solution to the complex eigenvalues

The matrix perturbation method is extended to discrete linear nonconservative system with unsymmetrical matrices in this article. By introducing the concept of the adjoint complex eigenvector and by making use of the orthogonality relationship in the complex mode theory, the first-order perturbation solution to the complex eigenvalues is derived. Numerical example shows that this method is efficient and practicable.     

Linear and non‐linear stability analysis for finite difference discretizations of high‐order Boussinesq equations

This paper considers a method of lines stability analysis for finite difference discretizations of a recently published Boussinesq method for the study of highly non‐linear and extremely dispersive water waves. The analysis demonstrates the near‐equivalence of classical linear Fourier (von Neumann) techniques with matrix‐based methods for formulations in both one and two horizontal dimensions. The matrix‐based method is also extended to show the local de‐stabilizing effects of the non‐linear terms, as well as the stabilizing effects of numerical dissipation. A comparison of the relative stability of rotational and irrotational formulations in two horizontal dimensions provides evidence that the irrotational formulation has significantly better stability properties when the deep‐water non‐linearity is high, particularly on refined grids. Computation of matrix pseudospectra shows that the system is only moderately non‐normal, suggesting that the eigenvalues are likely suitable for analysis purposes. Numerical experiments demonstrate excellent agreement with the linear analysis, and good qualitative agreement with the local non‐linear analysis. The various methods of analysis combine to provide significant insight into the numerical behaviour of this rather complicated system of non‐linear PDEs. Copyright © 2004 John Wiley & Sons, Ltd.     

On the computation of the main eigen-pairs of the continuous-time linear quadratic control problem

The degeneration of the eigenvalue equation of the discrete-time linear quadraticcontrol problem to the continuous-time one when△t→0?is given first.When thecontinuous-time n-dimensional eigenvalue equation,which has all the eigenvalues located inthe left half plane,has been reduced from the original2n-dimensional one,the present paperproposes that several of the eigenvalues nearest to the imaginary axis be obtained by thematrix transformation A_e=e~A.All the eigenvalues of A_e are in the unit circle,with theeigenvectors unchanged and the original eigenvaiues can be obtained by a logarithmoperation.And several of the eigenvalues of A_e nearest to the unit circle can be calculated bythe dual subspace iteration method.     

Dispersion relation in the limit of high frequency for a hyperbolic system with multiple eigenvalues

The results of a previous paper (Muracchini et al., 1992) are generalized by considering a hyperbolic system in one space dimension with multiple eigenvalues. The dispersion relation for linear plane waves in the high-frequency limit is analyzed and the recurrence formulas for the phase velocity and the attenuation factor are derived in terms of the coefficients of a formal series expansion in powers of the reciprocal of frequency. In the case of multiple eigenvalues, it is also verified that linear stability implies λ$\lambda$-stability for the waves of weak discontinuity. Moreover, for the linearized system, the relationship between entropy and stability is studied. When the nonzero eigenvalue is simple, the results of the paper mentioned above are recovered. In order to illustrate the procedure, an example of the linear hyperbolic system is presented in which, depending on the values of parameters, the multiplicity of nonzero eigenvalues is either one or two. This example describes the dynamics of a mixture of two interacting phonon gases.     

Several solution methods for the generalized complex eigenvalue problem with bounded uncertainties

The aim of this paper is to evaluate the effects of uncertain-but-bounded parameters on the complex eigenvalues of the non-proportional damping structures. By combining the interval mathematics and the finite element analysis, the mass matrix, the damping matrix and the stiffness matrix were represented as the interval matrices. Firstly, with the help of the optimization theory, we presented an exact solution—the vertex solution theorem, for determining the exact upper bounds or maximum values and exact lower bounds or minimum values of complex eigenvalues of structures, where the extreme values are reached on the boundary of the interval mass, damping and stiffness matrices. Then, an interval perturbation method was proposed, which needs less computational efforts. A numerical example of a seven degree-of-freedom spring-damping-mass system was used to illustrate the computational aspects of the presented vertex solution theorem and the interval perturbation method in comparison with Deif’s method.     

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

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

Roe‐type Riemann solver for gas–liquid flows using drift‐flux model with an approximate form of the Jacobian matrix

This work presents an approximate Riemann solver to the transient isothermal drift ‐ flux model. The set of equations constitutes a non‐linear hyperbolic system of conservation laws in one space dimension. The elements of the Jacobian matrix A are expressed through exact analytical expressions. It is also proposed a simplified form of A considering the square of the gas to liquid sound velocity ratio much lower than one. This approximation aims to express the eigenvalues through simpler algebraic expressions. A numerical method based on the Gudunov's fluxes is proposed employing an upwind and a high order scheme. The Roe linearization is applied to the simplified form of A . The proposed solver is validated against three benchmark solutions and two experimental pipe flow data. Copyright © 2015 John Wiley & Sons, Ltd.     

基于奇异值分解的复模态矩阵摄动法

对非自伴随系统的振动重分析问题,提出了一种简单的通用方法。从子空间缩聚出发,基于复矩阵的奇异值分解定理,推导了同时适用于孤立 特征值,相重特征值和相近特征值三种复特征值情况的一阶和二阶摄动公式。算例表明,该方法通用性好,且具有足够的精度。     

非半单分叉的Normal Form计算

在文［１］基础上，提出一种仅知道派生线性系统零实部特征值时求解非线性系统非半单分叉ＮｏｒｍａｌＦｏｒｍ的方法．通过适当的分类，将要求解的线性代数方程组分为若干相互独立的方程组．将所求系数向量按字典序列排列后，各独立方程组的系数矩阵是上三角矩阵．在非共振情形，各系数向量可按反字典序列递推求出．在共振情形，根据文中的二个定理，巧妙地由一简单的常数矩阵的最大秩子矩阵，定位其系数矩阵的满秩子矩阵，解决了这类方程组的降维简化．通过消元法，把简化后的方程化成类似于半单分叉ＮｏｒｍａｌＦｏｒｍ求解过程中方程的形式，其解法也类似．该方法非常易于在计算机代数软件平台上程序化．     

Linear stability analysis of flow in a periodically grooved channel

We have conducted the linear stability analysis of flow in a channel with periodically grooved parts by using the spectral element method. The channel is composed of parallel plates with rectangular grooves on one side in a streamwise direction. The flow field is assumed to be two‐dimensional and fully developed. At a relatively small Reynolds number, the flow is in a steady‐state, whereas a self‐sustained oscillatory flow occurs at a critical Reynolds number as a result of Hopf bifurcation due to an oscillatory instability mode. In order to evaluate the critical Reynolds number, the linear stability theory is applied to the complex laminar flow in the periodically grooved channel by constituting the generalized eigenvalue problem of matrix form using a penalty‐function method. The critical Reynolds number can be determined by the sign of a linear growth rate of the eigenvalues. It is found that the bifurcation occurs due to the oscillatory instability mode which has a period two times as long as the channel period. Copyright © 2003 John Wiley & Sons, Ltd.