非定常短轴承油膜力公式的变分修正

本文采用了变分方法对非定常短轴承的油膜压力分布公式进行了修正，既保留了短轴承公式的简洁形式，又使其适用于大长径比轴承。得出了具有足够精度、适合轴颈大扰动情况下的有限长圆柱轴承非定常油膜力的解析公式。与差分充零算法相比，短轴承公式的结果在轴承长径比为0．6时，误差已经超过百分之二百，而本方法计算结果的误差小于百分之五。因此采用本方法既提高了短轴承油膜力公式的计算精度，又保持了油膜力公式的简洁形式，不失为进行转子-轴承系统非线性动力分析的一种有效方法。     

考虑进油压力的滑动轴承非线性油膜力数据库

通过对ReyTlolds方程的非线性变换，提出了考虑进油压力边界条件时径向滑动轴承非线性油膜力数据库的建库方法，扩展了油膜力数据库计算方法的应用，通过引入2个有限数据域的新变量对转子轴心速度项和进油压力边界条件进行有限化处理，得到了更符合实际工况的连续性油膜力数据库及计算模型，同有限元法对比分析了非线性油膜力数据库的适用性．结果表明，非线性油膜力数据库模型的精度较高，所建立的非线性油膜力计算模型可用于对转子系统瞬态运动进行简便和快捷的分析．     

非线性振动系统周期运动及其稳定性的数值研究

§1引言确定型非线性振动系统的运动可分类如下: 1.非定常运动;2.定常运动:(1)周期运动,(2)各态历经运动,(3)浑沌运动。其中非定常运动是一暂态过程,会随着时间的增长逐步衰减乃至实际上消失。定常运动中的各态历经运动,指系统至少有两个互不通约(即其比值为无理数)的振动频率,因此运动虽然局      

大型转子-基础-地基系统的非线性动力分析

针对实际工程中的大型机组，在线性理论分析基础上，引入转子系统的非线性油膜力项，采用子结构模态综合法，形成一个比较接近实际大型汽轮发电机组的包括陀螺转子-非稳态非线性油膜转承-弹性基础-地基系统的非线性系统计算模型。通过对系统方程进行分块直接积分求解，得到了不同位置的轴承在不同转速和不同转子偏心量下引起的系统非线性动力学现象，为大机组的非线性分析和改进提供较完善的理论分析和计算的基础。     

用瑞利-李兹法求解瞬时非稳态滑动轴承油膜力的新算法

提出了一个利用瑞利一李兹方法求解Reynolds边界条件下非稳态滑动轴承油膜力的近似算法,充分利用油膜力分布的特性,用双曲余弦函数来表示油膜力的轴向压力分布,而用多项式函数插值法来求解油膜力的周向压力分布,并同时计算出油膜力的破裂边界。算例表明本算法达到了相当高的精度,可用于转子系统的非线性数值分析,能大大降低数值求解瞬态油膜力的计算时间。     

Parameter variation and data mining of oil-film bearings: a stochastic study on the Reynolds’s equation of lubrication

The modeling and simulation process of oil-film bearing dynamics constitutes a rather essential task integrated in the workflow of various mechanical products. Specifically, in the turbo charger industry, the correct capture and understanding of the associated nonlinear rotating dynamics is of utmost importance, since the system’s efficiency and lifetime span depends on it. The root cause of the nonlinear rotordynamic effects is the oil-film concentrated in the rotor’s journal bearings. Its behavior is highly coupled with both the system’s geometric and dynamic configuration. The dynamics of the oil-film are described by the well-known Navier–Stokes equation, which under a series of assumptions and simplifications results to the, so-called, Reynolds equation. In this paper, the Reynolds equation is numerically solved based on a finite difference scheme and several parameter variation studies are conducted in an effort to pinpoint the most influential parameters—journal bearing geometric dimensions, oil-film properties and rotor-velocity-driven inputs—with respect to designated responses—friction, oil-film pressure force, minimum oil-film thickness and boundary oil-flow—all of which are regarded as important in terms of the aforementioned system’s efficiency and lifetime span. Based on multivariate analysis algorithms, correlation outcomes and global sensitivity results are presented. In an effort to capture possible nonlinear phenomena, which might not be possible via linear data mining tools, the Spearman rank-order coefficient and self-organizing maps methodology are applied.     

有限长轴承非稳态油膜力自由-移动边界问题的有限条解法

讨论了不可压缩流体润滑的动载径向滑动轴承油膜压力分布的自由移动边界问题的有限条解法．将自由边界问题转化为全域(矩形域)的具有不等式约束的微分方程的边值问题，进一步化为具有不等式约束的泛函极值问题。借助有限条法在矩形域上离散这个泛函，得到了一个特殊的二次泛函的规化问题。通过变量平移变换，使其化为标准的二次规划问题。然后借助于牛顿非光滑算法，迭代求解非线性的互补方程。给出了有限长轴承真实的油膜应力分布。对于所求解方程的系数矩阵的高度稀疏性。给出了紧缩存储算法。节省了存储空间和减少了计算量。算例表明该方法是有效的。     

受电磁激励的连续转子-轴承系统非线性振动与分岔

研究了受横向不平衡电磁激励的转子．轴承系统的非线性振动响应。首先将转子．轴承系统简化为带有质量不平衡并受横向激励的连续梁,由于短轴承的油膜力和电磁力的共同激励,系统振动具有强非线性特性。用Galerkin方法把偏微分控制方程离散为常微分方程组,采用四阶Runge—Kutta法对该系统进行数值仿真研究。其次比较了转轴分别在电磁力、油膜力单独作用和两种力共同作用下的振动特性,研究表明电磁力和油膜力对转子系统的非线性振动和分岔有着不同的贡献：油膜力的存在抑制了拟周期运动的发生,延长了稳定运行区域;电磁力拉长了拟周期发生的区域,降低了转子系统发生突发性破坏的风险。最后给出了系统响应随转速、电磁参数、油膜粘度等控制参数变化的分岔图,表明：系统在两个方向的运动随控制参数的变化趋势基本相同,经历了周期、倍周期、拟周期等非线性运动交替出现的过程;且油膜粘度的增大有利于转子系统的安全运行。     

A highly-efficient method for stationary response of multi-degree-of-freedom nonlinear stochastic systems

Analytical and numerical studies of multi-degree-of-freedom(MDOF) nonlinear stochastic or deterministic dynamic systems have long been a technical challenge.This paper presents a highly-efficient method for determining the stationary probability density functions(PDFs) of MDOF nonlinear systems subjected to both additive and multiplicative Gaussian white noises. The proposed method takes advantages of the sufficient conditions of the reduced Fokker-Planck-Kolmogorov(FPK) equation when constructing the trial solution. The assumed solution consists of the analytically constructed trial solutions satisfying the sufficient conditions and an exponential polynomial of the state variables, and delivers a high accuracy of the solution because the analytically constructed trial solutions capture the main characteristics of the nonlinear system. We also make use of the concept from the data-science and propose a symbolic integration over a hypercube to replace the numerical integrations in a higher-dimensional space, which has been regarded as the insurmountable difficulty in the classical method of weighted residuals or stochastic averaging for high-dimensional dynamic systems. Three illustrative examples of MDOF nonlinear systems are analyzed in detail. The accuracy of the numerical results is validated by comparison with the Monte Carlo simulation(MCS) or the available exact solution. Furthermore, we also show the substantial gain in the computational efficiency of the proposed method compared with the MCS.     

Application of a modified Lindstedt–Poincaré method in coupled TDOF systems with quadratic nonlinearity and a constant external excitation

An analytical approach is developed for the nonlinear oscillation of a conservative, two-degree-of-freedom (TDOF) mass-spring system with serial combined linear–nonlinear stiffness excited by a constant external force. The main idea of the proposed approach lies in two categories, the first one is the transformation of two nonlinear differential equations of a two-mass system using suitable intermediate variables into a single nonlinear differential equation. Another is the treatment a quadratic nonlinear oscillator (QNO) by the modified Lindstedt–Poincaré (L-P) method presented recently by the authors. The first-order and second-order analytical approximations for the modified L-P method are established for the QNOs with satisfactory results. After solving the nonlinear differential equation, the displacements of two-mass system can be obtained directly from the governing linear second-order differential equation. Unlike the common perturbation method, the modified L-P method is valid for weak as well as strong nonlinear oscillation systems. On the other hand, the new approach yields simple approximate analytical expressions valid for small as well as large amplitudes of oscillation. In short, this new approach yields extended scope of applicability, simplicity, flexibility in application, and avoidance of complicated numerical integration as compared to the previous approaches such as the perturbation and classical harmonic balance methods. Two examples of nonlinear TDOF mass-spring systems excited by a constant external force are selected and the approximate solutions are verified with the exact solutions derived from the Jacobi elliptic function and also the numerical fourth-order Runge–Kutta solutions.     

Nonlinear time transformation method for strong nonlinear oscillation systems

In this paper, a nonlinear time transformation method is presented for the analysis of strong nonlinear oscillation systems. This method can be used to study the limit cycle behavior of the autonomous systems and to analyze the forced vibration of a strong nonlinear system. The project partly supported by the Foundation of Zhongshan University Advanced Research Center     

Study of differential control method for solving chaotic solutions of nonlinear dynamic system

In this paper, we focus on the need to solve chaotic solutions of high-dimensional nonlinear dynamic systems of which the analytic solution is difficult to obtain. For this purpose, a Differential Control Method (DCM) is proposed based on the Mechanized Mathematics-Wu Elimination Method (WEM). By sampling, the computer time of the differential operator of the nonlinear differential equation can be substituted by the differential quotient of solving the variable time of the sample. The main emphasis of DCM is placed on substituting the differential quotient of a small neighborhood of the sample time of the computer system for the differential operator of the equations studied. The approximate analytical chaotic solutions of the nonlinear differential equations governing the high-dimensional dynamic system can be obtained by the method proposed. In order to increase the computational efficiency of the method proposed, a thermodynamics modeling method is used to decouple the variable and reduce the dimension of the system studied. The validity of the method proposed for obtaining approximate analytical chaotic solutions of the nonlinear differential equations is illustrated by the example of a turbo-generator system. This work is applied to solving a type of nonlinear system of which the dynamic behaviors can be described by nonlinear differential equations.     

NONLINEAR DYNAMICAL BIFURCATION AND CHAOTIC MOTION OF SHALLOW CONICAL LATTICE SHELL

The nonlinear dynamical equations of axle symmetry are established by the method of quasi-shells for three-dimensional shallow conical single-layer lattice shells. The compatible equations are given in geometrical nonlinear range. A nonlinear differential equation containing the second and the third order nonlinear items is derived under the boundary conditions of fixed and clamped edges by the method of Galerkin. The problem of bifurcation is discussed by solving the Floquet exponent. In order to study chaotic motion, the equations of free oscillation of a kind of nonlinear dynamics system are solved. Then an exact solution to nonlinear free oscillation of the shallow conical single-layer lattice shell is found as well. The critical conditions of chaotic motion are obtained by solving Melnikov functions, some phase planes are drawn by using digital simulation proving the existence of chaotic motion.     

Reduced augmented Lagrangian bi-conjugate gradient method for impact-contact problems

To avoid the numerical oscillation of the penalty method and non-compatibility with explicit operators of conventional Lagrange multiplier methods used in transient contact problems to enforce surface contact conditions, a new approach to enforcing surface contact constraints for the transient nonlinear finite element problems, referred to as “the reduced augmented Lagrangian bi-conjugate gradient method (ALCG)”, is developed in this paper. Based on the nonlinear constrained optimization theory and is compatible with the explicit time integration scheme, this approach can also be used in implicit scheme naturally. The new surface contact constraint method presented has significant advantages over the widely adopted penalty function methods and the conventional Lagrangian multiplier methods. The surface contact constraints are satisfied more accurately for each step by the algorithm, so the oscillation of numerical solution for the explicit scheme is depressed. Through the development of new iteration strategy for solving nonlinear equations, ALCG method improves the computational efficiency greatly. Project supported by State Education Commission Doctoral Foundation and Natural Science Foundation of Liaoning Province.     

Nonlinear modeling and dynamic analysis of the rotor-bearing system

To overcome the shortcomings of extreme time-consuming in solving the Reynolds equation, two efficient calculation methods, based on the free boundary theory and variational principles for the unsteady nonlinear Reynolds equation in the condition of Reynolds boundary, are presented in the paper. By employing the two mentioned methods, the nonlinear dynamic forces as well as their Jacobians of the journal bearing can be calculated saving time but with the same accuracy. Of these two methods, the one is called a Ritz model which manipulates the cavitation region by simply introducing a parameter to match the free boundary condition and, as a result, a very simple approximate formulae of oil-film pressure is being obtained. The other one is a one-dimensional FEM method which reduces the two-dimensional variational inequality to the one-dimensional algebraic complementary equations, and then a direct method is being used to solve these complementary equations, without the need of iterations, and the free boundary condition can be automatically satisfied. Meanwhile, a new order reduction method is contributed to reduce the degrees of freedom of a complex rotor-bearing system. Thus the nonlinear behavior analysis of the rotor-bearing system can be studied time-sparingly. The results in the paper show the high efficiency of the two methods as well as the abundant nonlinear phenomenon of the system, compared with the results obtained by the usual numerical solution of the Reynolds equation.     

Series solutions of annular axisymmetric stagnation flow and heat transfer on moving cylinder

The steady, laminar, incompressible flow and heat transfer of a viscous fluid between two circular cylinders for two different types of thermal boundary conditions are investigated. The governing Navier-Stokes and thermal equations of the flow are reduced to a nonlinear system of ordinary differential equations. The equations are solved analyt- ically using the homotopy analysis method （HAM）. Convergence of the HAM solutions is discussed in detail. These solutions are then compared with recently obtained numericM and perturbative solutions. Plots of the velocity and temperature profiles are provided for various values of the relevant parameters.     

PMUCR方法在高维非线性动力学系统中的应用

对现有的若干种胞映射方法应用于高维非线性动力学系统全局分析时存在的局限性进行了分析,在此基础上,提出了胞映射方法应用于高维系统时需遵循的几个原则:选择合适的分析平面;减少数据量;确定吸引子在相空间中的位置。根据这些原则,对胞参考点映射法(PMUCR)进行了改进,并对该改进方法应用于2维系统和4维系统时消耗的相对CPU时间进行了比较,结果表明,该改进方法能有效地应用于高维非线性动力学系统中。最后以船舶机械非线性隔振系统为例,分别分析了该改进方法在系统呈现周期运动和混沌运动时的应用。     

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.     

Robust control for MIMO hybrid dynamical system of underwater vehicles by composite adaptive fuzzy estimation of uncertainties

An organization structure of global oscillation with respect to a cracked rotor system with oil-film force is investigated in this paper. We profit from GPU cluster parallel computing to present a number of high-quality phase diagrams, and exhibit global dynamic characteristics of the system. An interesting scenario, “eye” of chaos, is discovered in this cracked rotor system, emerging as the accumulation limit of forward and reverse period-doubling bifurcation cascades. In this system, it is a common phenomenon that the vibration response of the rotor presents three typical characteristics in parameter space with the rotation speed increasing. Moreover, these phase diagrams assist us to identify multi-attractor coexisting that makes the dynamics behavior of this system become more enrich and complex. These results we represent get us better to understand the nonlinear response of the cracked rotor system and are beneficial to control and diagnose the crack.     

1/3 subharmonic solution of elliptical sandwich plates

IntroductionInrecentyears,withtheessentialadvantageoflightweightandhighrigidity ,sandwichplatesandshellshavebeenusedasanimportantpatternofstructuralelementsinaeronautical,astronauticalandnavalengineering .However,nonlinearproblemsforsandwichplatesandshellsareonlyinvestigatedbyafewbecauseofthedifficultiesofnonlinearmathematicalproblems.LiuRen_huaiandXuJia_chu[1,2 ]andothershavemadesomeinvestigationsinthisfield .Bifurcationofnonlinearvibrationforsandwichplateshasnotyetbeeninvestigated .Inthisp…