非线性转子-轴承系统的周期解及近似解析表达式

通过对普通打靶方法进行改造提出一种确定非线性系统周期轨道及周期的新型打靶算法。首先通过改变系统的时间尺度，将非线性系统周期轨道的周期显式地出现在非线性系统的系统方程中，然后对传统打靶法进行改造，将周期也作为一个参数一起参与打靶法的迭代过程，迭代过程包含对周期轨道和周期的求解，迭代过程中的增量通过优化方法选择，从而能迅速确定出系统的周期轨道及其周期。应用所求的结果结合谐波平衡方法求得了非线性系统的周期轨道的近似解析表达式，理论上通过增加谐波的阶数任何精度的周期解都可以得到。最后将该方法应用于非线性转子轴承系统，求出了在某些参数下转子的周期解及其近似解析表达式，通过与四阶Runge-Kutta数值积分结果比较，验证了方法的有效性，计算结果对于转子系统运动的定量控制有重要理论指导意义。     

求非线性转子－轴承系统周期响应的一种计算方法

本文把求非线性转子－轴承系统瞬态响应的分块Ｎｅｗｍａｒｋ方法与打靶法结合求系统的周期解，该方法利用了分块Ｎｅｗｍａｒｋ方法求解速度快的优点和Ｊａｃｏｂｉ矩阵求解时每步不需迭代的特点。本文首先给出周期解的求解公式，然后用一个算例，讨论了周期解的稳定性及失稳后的分岔行为。     

求解非线性动力系统周期解推广的打靶法

提出一种确定非线性系统周期轨道及周期的改进打靶算法。首先通过改变系统的时间尺度，将非线性系统周期轨道的周期显式地出现在非线性系统的系统方程中，然后对传统打靶法进行改造，将周期也作为一个参数一起参入打靶法的迭代过程，从而能迅速确定出系统的周期轨道及其周期。该方法对初始迭代参数没有苛刻要求，可以用于分析强非线性系统，而且对参数激励系统同样有效，对高维系统也能迅速、准确地求得周期解。文中应用该方法对三维Rǒssler系统和八维非线性柔性转子-轴承系统的周期轨道和周期进行了求解，通过与四阶Runge-Kutta数值积分结果比较，验证了方法的有效性。     

旋转振动圆柱绕流周期解和Floquet稳定性

对低雷诺数旋转振动圆柱绕流问题运用低维Galerkin方法将N－S方程约化为一组非线性常微分方程组。运用打靶法数值求解了这组方程的周期解,并用Tloquet理论对周期解的稳定性进行了分析,确定了流动失稳的机制。     

悬索在外激励作用下的1:3内共振分析(Ⅱ):数值结果

本研究的第一部分已经推导了悬索在第一阶面内对称模态主共振和第三阶面内对称模态主共振下的平均方程,其中考虑了这两阶模态之间1∶3内共振.本文对平均方程的稳态解,周期解以及混沌解进行了研究.利用 Newton-Naphson 方法和拟弧长的延拓算法确定了主共振情况下的幅频响应曲线,通过利用 Jacobian 矩阵的特征值判断幅频响应曲线中解的稳定性.在这些幅频响应曲线中.都存在超临界 Hopf 分叉,导致平均方程的周期解.以这些超临界 Hopf 分叉为起点.利用打靶法和拟弧长的延拓算法确定了两种主共振情况下的周期解分支,同时通过利用 Floquet 理论判断这些周期解的稳定性.然后利用数值结果研究了两种主共振情况下的厨期解经过倍周期分叉通向混沌的过程.最后利用 Runge-Kutta 法研究了悬索两自由度离散模型的非线性响应.     

一种非线性系统周期解的延拓算法

提出了一种非线性系统周期解的延拓算法。指出了非线性系统周期解在分岔点处由于雅可比矩阵奇异而导致一般延拓方法延拓失败问题;然后基于推广的打靶法的思想,将普通延拓算法推广,提出了一种用于周期解延拓的算法。对于非线性动力系统,该算法可以在已知某一参数下的周期解的基础上,求解出在一定参数范围内非线性动力系统的解随参数的连续变化情况。应用该方法对非线性柔性转子-轴承系统的周期解与参数的依赖关系进行了求解,验证了方法的有效性。     

多自由度非线性系统的同伦分析方法

大部分工程实际问题可以用多自由度非线性系统来描述,这些系统的数学模型是许多个耦合的两阶常微分方程.一般地,要精确求解这些方程非常困难,因此可以考虑它们的解析近似解.同伦分析方法是解非线性系统响应的有用工具,本文将它应用于多自由度非线性系统的求解中.利用求两自由度耦合van del Pol振子周期解的实例,展示了同伦分析方法的有效性和巨大潜力.同时,把得到的解析近似解与系统的Runge-Kutta数值解作了比较,结果表明同伦分析方法是求解多自由度非线性系统的有效方法.     

悬索在外激励作用下的1∶3内共振分析(II)∶数值结果

本研究的第一部分已经推导了悬索在第一阶面内对称模态主共振和第三阶面内对称模态主共振下的平均方程,其中考虑了这两阶模态之间1∶3内共振。本文对平均方程的稳态解、周期解以及混沌解进行了研究。利用Newton-Naphson方法和拟弧长的延拓算法确定了主共振情况下的幅频响应曲线,通过利用Jacobian矩阵的特征值判断幅频响应曲线中解的稳定性。在这些幅频响应曲线中,都存在超临界Hopf分叉,导致平均方程的周期解。以这些超临界Hopf分叉为起点,利用打靶法和拟弧长的延拓算法确定了两种主共振情况下的周期解分支,同时通过利用Floquet理论判断这些周期解的稳定性。然后利用数值结果研究了两种主共振情况下的周期解经过倍周期分叉通向混沌的过程。最后利用Runge-Kutta法研究了悬索两自由度离散模型的非线性响应。     

含高阶余项的非线性动力系统数值计算法

建立了一种求解非线性动力系统高精度数值计算的新方法,重构了等价的非线性动力系统方程,该方程考虑了非线性函数的任意高阶项,并给出了该方程的Duhamel积分表达式,在时间步长内用Newton-Raphson法进行数值迭代求解,该方法能连续满足微分方程而不只是在离散的步长端点满足方程,从而打破了传统的Euler型有限差分法。计算实例表明,该方法计算精度高于传统的Runge-Kutta,Newmark-β和Wilson-θ等方法。     

齿轮－转子－滑动轴承系统时变非线性动力特性研究

本文应用求周期解的数值计算方法─—打靶法和判定周期解稳定性的Ｆｌｏｑｕｅｔ乘子研究了齿轮－转子－滑动轴承系统中齿轮啮合时变刚度，滑动轴承非线性特性对转子系统不平衡响应和失稳的影响，并比较了平衡位置失稳和不平衡响应周期解失稳，以及按双轴计算与单轴计算结果的差别，为工程设计理论计算提供基础。     

Numerical solutions of incompressible Euler and Navier-Stokes equations by efficient discrete singular convolution method

An efficient discrete singular convolution (DSC) method is introduced to the numerical solutions of incompressible Euler and Navier-Stokes equations with periodic boundary conditions. Two numerical tests of two-dimensional Navier-Stokes equations with periodic boundary conditions and Euler equations for doubly periodic shear layer flows are carried out by using the DSC method for spatial derivatives and fourth-order Runge-Kutta method for time advancement, respectively. The computational results show that the DSC method is efficient and robust for solving the problems of incompressible flows, and has the potential of being extended to numerically solve much broader problems in fluid dynamics. The project supported by the National Natural Science Foundation of China (No.19902010).     

A high-accuracy temporal–spatial pseudospectral method for time-periodic unsteady fluid flow and heat transfer problems

A temporal–spatial pseudospectral (TSP) method is proposed for the high-accuracy solutions of time-periodic unsteady fluid flow and heat transfer problems. In this method, both the spatial and temporal derivative terms in the governing equations are computed by pseudospectral method. The spatial derivatives are computed through Chebyshev and Lagrange polynomials while the time derivatives are computed by Fourier series. The TSP method is capable of directly finding out the periodic state solutions without the necessity to resolve the initial transient state solutions, hence holds high computational efficiency and high numerical accuracy properties for the time-periodic problems. This method is validated by three 2D benchmark problems: the time-periodic incompressible flow with exact solutions; the natural convection in enclosure with time-periodic temperature on one sidewall, and on both sidewalls. The TSP results fit well the exact solutions or the benchmark solutions and the TSP accuracy is much higher than the time marching spatial pseudospectral accuracy. Some time-dependent fluid flow and heat transfer characteristic parameters are analysed. The proposed TSP method could be further extended to more complex time-periodic unsteady fluid flow and heat transfer problems where high-accuracy results are required.     

A numerical treatment of the periodic solutions of non-linear vibration systems

Direct numerical integration can be used to find the periodicsolutions for the equations of motion of nonlinear vibrationsystems.The initial conditions are iterated so that theycoincide With the terminal conditions.The time interval ofthe integration(i.e.,the period)and certain parameters ofthe equations of motion can be included in the iterations.Theintegration method has a variable stoplength.This Sbooting method can produce periodic solutions witha shorter computex time.The only error occurs in the numeri-cal integration and it can therefore be estimated and madesmall enough.Using this method one can treat a variety ofvibration problems.such as free conservative.forced.para-meter-excited and self-sustained vibrations with one or se-veral degrees-of-freedom.Unstable solutions and those Whichare sensitive to parameter Changes can also be calculated.Thestability of the solutions is investigated based on the thecryof differential equations with periodic coefficients.The ex-trapolation method and the proc     

一种新的求解非线性系统周期解方法

本文利用切比雪夫多项式的若干良好性质,对非自治强非线性动力系统进行分析。将状态矢量在主周期上先展开谐波级数的形式,再将各谐波展开为切比雪夫多项式的形式,从而将求周期解的问题转变为非线性代数方程组的求解问题,得出一种可以方便、迅速地获得近似周期解的解析方法。此方法不依赖于小参数假设,可以用于分析强非线性问题和高维问题,而且对参数激励系统同样有效。以Duffing系统周期解的计算为例,通过与标准谐波平衡方法和四阶Runge-Kutta数值积分结果比较,说明此方法的有效性。     

Subharmonic resonance of a clamped-clamped buckled beam with 1:1 internal resonance under base harmonic excitations

The subharmonic resonance and bifurcations of a clamped-clamped buckled beam under base harmonic excitations are investigated. The nonlinear partial integrodifferential equation of the motion of the buckled beam with both quadratic and cubic nonlinearities is given by using Hamilton's principle. A set of second-order nonlinear ordinary differential equations are obtained by spatial discretization with the Galerkin method. A high-dimensional model of the buckled beam is derived, concerning nonlinear coupling. The incremental harmonic balance (IHB) method is used to achieve the periodic solutions of the high-dimensional model of the buckled beam to observe the nonlinear frequency response curve and the nonlinear amplitude response curve, and the Floquet theory is used to analyze the stability of the periodic solutions. Attention is focused on the subharmonic resonance caused by the internal resonance as the excitation frequency near twice of the first natural frequency of the buckled beam with/without the antisymmetric modes being excited. Bifurcations including the saddle-node, Hopf, perioddoubling, and symmetry-breaking bifurcations are observed. Furthermore, quasi-periodic motion is observed by using the fourth-order Runge-Kutta method, which results from the Hopf bifurcation of the response of the buckled beam with the anti-symmetric modes being excited.     

强非线性动力系统的频率增量法

提出一类强非线性动力系统的暧时频率增量法,将描述动力系统的二阶常微分方程,化为以相位为自变量、瞬廛频率为未知函数的积分方程;用谐波平衡原理,将求解瞬时频率的积分问题,归结为求解以频率增量的Fourier系数为独立变量的线性代数方程组;给出了若干例子。     

Bifurcation and route-to-chaos analyses for Mathieu–Duffing oscillator by the incremental harmonic balance method

Bifurcations and route to chaos of the Mathieu–Duffing oscillator are investigated by the incremental harmonic balance (IHB) procedure. A new scheme for selecting the initial value conditions is presented for predicting the higher order periodic solutions. A series of period-doubling bifurcation points and the threshold value of the control parameter at the onset of chaos can be calculated by the present procedure. A sequence of period-doubling bifurcation points of the oscillator are identified and found to obey the universal scale law approximately. The bifurcation diagram and phase portraits obtained by the IHB method are presented to confirm the period-doubling route-to-chaos qualitatively. It can also be noted that the phase portraits and bifurcation points agree well with those obtained by numerical time-integration.     

一类线性Liénard方程周期边值问题解的存在唯一性的构造性证明

利用初值问题方法给出了一类线性Liénard方程周期边值问题解的存在唯一性的构造性证明,利用数值延拓方法,并给出了计算实例和一种大范围求解这类方程周期解的方法。