传递矩阵—多项式法的改进及其应用

本文提出了一种计算转子－滑动轴承－挤压油膜阻尼器转子系统的阻尼有频率及稳定性的方法，即将传递矩阵－多项式法改进后，应用于这类具有分叉结点的系统，直接得到系统的特征多项式，然后用Ｂａｉｒｓｔｏｗ－Ｎｅｗｔｏｎ人僻因子法求了出全部特征根。本方法具有占内存少、计算速度快、数值稳定等特点，文中对一转子模型进行了计算，其结果与实测基本相符。     

挤压油膜阻尼器－滑动轴承－刚性转子系统的稳定性及分岔行为

对挤压油膜阻尼器－滑动轴承－转子系统的稳定性及分岔行为进行了研究，由于该动力系统为一强非线性系统，具有复杂的非线性现象。本文采用Ｆｌｏｑｕｅｔ理论对其周期解的稳定性进行了计算分析：随着系统参数的变化，该系统将出现稳态周期解、准周期分岔、倍周期分岔。文中也对系统平衡点的稳定性进行了分析，讨论了其Ｈｏｐｆ分岔行为     

传递矩阵法在材料力学中的应用

将传递矩阵法应用于材力中的梁、轴问题,优点是能同时求出内力和变形,无需区分静定与静不定问题,易于计算机计算     

非线性系统全局动态特性分析的PCM法及其在转子轴承系统中的应用

本文将Ｐｏｉｎｃａｒｅ映射的思想与胞映射法相结合，提出了可用于高维非线性动力系统全局稳定性分析的新型数值方法：ＰＣＭ（Ｐｏｉｎｃａｒｅ－Ｃｅｌｌ－Ｍａｐｐｉｎｇ）法，和胞映射法相比，新方法在实用上具有明显的优点。为说明ＰＣＭ法的有效性，本文应用此方法对平衡转子轴承非线性动力系统进行了全局稳定性分析，同时给出了一确定状态空间中存在的所有周期解及其吸引域。     

Riccati传递矩阵法在流体动力学中的应用

传统的传递矩阵法在求解管内流体的固有频率时,因奇点干扰出现漏根现象.本文结合气体压缩机管道,运用转子动力学中的Riccati方法计算管内气体的固有频率及振型,从根本上消除了漏根现象.算例表明了本方法的可靠性.     

转子系统阻尼固有频率的传递矩阵—多项式计算法

本文提出了转子——支承系统阻尼固有频率的一种计算方法,即采用Riccati分块矩阵递推,直接得到系统的特征多项式L(S)=0,然后用Bairstow-Newton法迭代解出全部特征根。和传统的传递矩阵法比,该法计算速度快,占用内存少,避免了数值不稳定和丢根,且可依次求得系统的全部解。文中对一模型转子和国产20万千瓦汽轮机转子系统做了计算,计算结果与实测值较好地接近。这个方法将有助于诸如汽轮发电机组这样的高速大型旋转机械的合理设计和可靠运行。     

求解转子—轴承系统动力特性的实模态与复模态结合算法

本文基于固定界面法推导了一种适于转子——轴承系统动力分析的实模态与复模态结合算法,并进行了理论分析与计算、结果表明,本文所推导的模态综合法在计算过程中保留了陀螺力矩的信息,同时计及了基础及滑动轴承油膜刚度、阻尼的影响,精度高,耗时少,具有工程实用价值。     

等步长高次多项式插值加速度法的局限性

动力时程分析中,在几个相邻的等长时间步之间对加速度的变化规律用多项式插值来描叙,经过推导可求解得到整个时间域上的动力方程的解答.根据泰勒展开原理分析表明,随着所取多项式次数的增加,收敛精度增高,计算步长适当放大,截断误差仍能在容许的范围之内.但是随着所取多项式次数增大,其算法的稳定域减少, 计算步长受到了此小稳定域的限制,收敛精度不再是所取计算步长宽度的决定因素,稳定域大小成了所取计算步长宽度的决定因素.因为一旦步长超出了此小稳定域范围,虽然在每个时间步内的截断误差不大,其传递的误差却会被放大到很多倍,最后导致计算结果严重失真.分析结果显示,多项式插值次数采用到步长的三次时,与一次多项式插值(对应线性加速度法)和二次多项式插值(对应二次加速度法)的分析方法相比,算法的稳定域急剧变窄,为h/T≤0.0099(h为计算步长,T为结构的固有周期),此小稳定域限制了计算步长的选择范围,其收敛精度很高因此可放大计算步长的优势无法施展.本文推导了三次加速度法的求解过程,进行了一个理想单自由度系统的动力时程分析计算,验证了结论的正确性.表明同时考虑收敛精度和稳定域来确定计算步长的宽度时,二次加速度法为优.     

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

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

一种中高频激励下动力系统响应的级数展开改进算法

针对频率响应函数的级数展开法在中高频激励时计算发散的问题,提出一种新的级数展开改进算法.将系统的结构模态划分为低阶和截断的高阶模态,在模态叠加分析的基础上,将频率响应函数进行泰勒级数展开.根据高低阶模态对质量矩阵和刚度矩阵的耦合特性,用低阶模态及系统矩阵表达高阶模态对响应的影响.研究结果表明,该算法将频率响应函数的级数展开法扩展到高频激励和中频激励范围阶段,在非完备模态条件下提高了频率响应函数的计算精度,数值计算检验了该方法准确可靠并有很好的收敛性.     

转子轴承系统电磁阻尼器的位移反馈振动控制

本文研究在周期激振力作用下转子轴承系统电磁阻尼器的位移反馈振动控制。首先给出受电磁阻尼器控制的一般线性多自由度转子轴承系统的振动方程，控制力和电磁阻尼器结构参数，控制电流等之间的关系式，然后提出在保证系统稳定的前提下，通过位移反馈控制其电流位移增益矩阵应满足的条件；其次以单盘对称转子轴承系统为例，求出转子系统受电磁阻器控制的幅频响应和相频响应；最后给出算例，打印出电磁阻尼器的静态电流，电流位移增益     

An Improved IMPES Method for Two-Phase Flow in Porous Media

In this paper we develop and numerically study an improved IMPES method for solving a partial differential coupled system for two-phase flow in a three-dimensional porous medium. This improved method utilizes an adaptive control strategy on the choice of a time step for saturation and takes a much larger time step for pressure than for the saturation. Through a stability analysis and a comparison with a simultaneous solution method, we show that this improved IMPES method is effective and efficient for the numerical simulation of two-phase flow and it is capable of solving two-phase coning problems.     

Conditions for noncircular whirling of nonlinear isotropic rotors

Nonlinear rotors are often considered as potential sources of chaotic vibrations. The aim of the present paper is that of studying in detail the behaviour of a nonlinear isotropic Jeffcott rotor, representing the simplest nonlinear rotor. The restoring and damping forces have been expanded in Taylor series obtaining a Duffing-type equation. The isotropic nature of the system allows circular whirling to be a solution at all rotational speeds. However there are ranges of rotational speed in which this solution is unstable and other, more complicated, solutions exist.The conditions for stability of circular whirling are first studied from closed form solutions of the mathematical model and then the conditions for the existence of solutions of other type are studied by numerical experimentation. Although attractors of the limit cycle type are often found, chaotic attractors were identified only in few very particular cases. An attractor supposedly of the last type reported in the literature was found, after a detailed analysis, to be related to a nonchaotic polyharmonic solution.As the typical unbalance response of isotropic nonlinear rotors has been shown to be a synchronous circular whirling motion, the convergence characteristics of Newton-Raphson algorithm applied to the solution of the set of nonlinear algebraic equations obtained from the differential equations of motion are studied in some detail. c damping coefficient i imaginaty unit (i=% MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbqfgBHr% xAU9gimLMBVrxEWvgarmWu51MyVXgaruWqVvNCPvMCG4uz3bqefqvA% Tv2CG4uz3bIuV1wyUbqee0evGueE0jxyaibaieYlf9irVeeu0dXdh9% vqqj-hEeeu0xXdbba9frFf0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea% 0dXdar-Jb9hs0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabe% aadaabauaaaOqaamaakaaabaGaeyOeI0IaaGymaaWcbeaaaaa!3E66!\[\sqrt { - 1}\]) k stiffness m mass t time x _istate variables i=1, 4 z complex co-ordinate (z=x+iy)[J] Jacobian matrix Oxyz inertial co-ordinate frame Oz rotating co-ordinate frame perturbation term eccentricity complex co-ordinate (=+i) system eigenvalues nonlinearity parameter nondimensional time phase spin speed u nonrotating t rotating0 amplitude t nondimensional termsNomenclature     

非线性动力系统全局分析的变胞胞映射法与转子／轴承系统的全局稳定性

在Ｐｏｉｎｃａｒｅ映射及胞映理论的基础上，提出了一种非线性动力系统全局分析的新方法－－变胞胞映射法，这种新方法改变了原胞映射法中胞在胞空间分布的不合理性及运算逻辑的不合理性，更适用于高维、大求解域非线性动力系统的求解。应用此方法，对具有非线性油膜力的Ｊｅｆｆｃｏｔ转子轴承系统进行了全局分析，绘制了系统分岔后的全局吸引域图，解释了一些工程中常见的非线性现象。     

滑动轴承挤压油膜阻尼器的试验研究

本文对滑动轴承挤压油膜阻尼器的关键技术问题进行了大量的试验研究,结合理论分析,设计出结构简单合理、性能可靠的阻尼器。并提出了适合于工程设计的结构参数。     

外弹性支承滑动轴承—刚性转子系统非线性动力稳定性的研究

研究外弹性支承滑动轴承－刚性转子系统的动力稳定性问题。建立了弹性支承圆瓦轴承－Jeffcot转子系统的力学模型,针对滑动轴承水平转子,在油膜力非线性的情况下,采用Newmark方法计算其响应。对于一个具体的算例,通过改变其参数,确定其稳定区,并与滑膜力线性化时理论分析的结果进行了比较。     

求解多体系统动力学方程的子系统迭代法

提出一种求解多体系统动力学方程的迭代法，将多体系统分割成一个主子系统和一个或多个辅子系统，各子系统通过分割点的位移、速度、加速度和力的协调相联系，从而大大降低了方程的耦合程度、非线性和刚性。文中详细讨论了迭代法的收敛条件。     

变截面梁动力响应的频域传递矩阵计算方法

为提高变截面梁振动分析的计算效率,提出了基于频域传递矩阵法的动力计算算法．首先,选择线速度、角速度、弯矩和剪力作为求解变量,通过Laplace变换将变截面梁的动力响应时域偏微分方程转换为频域常微分方程;然后,通过求解频域方程并结合协调和边界条件建立变截面梁的频域传递矩阵;通过构造傅里叶级数展开形式的时域响应函数,对变截面梁传递矩阵方法求解的频响函数进行Laplace逆变换,建立了变截面梁的固有特性计算和时域瞬态响应计算方法,最后,借助数值仿真软件,开发了变截面梁动力响应分析的计算程序．完成对算例的仿真计算和分析,并与有限元计算结果进行对比,数值结果验证了该方法的正确性和有效性．