Vibration of a rotating shaft with randomly varying internal damping

A simple Jeffcott rotor is considered with both external and internal damping. Coefficient of internal damping is subject to temporal random variations which may occasionally bring the rotor into the domain of dynamic instability. The corresponding sporadic outbreaks in the rotor's vibrational response (whirl) are studied by applying the Krylov-Bogoliubov averaging method to the complex equation of motion and using parabolic approximation for the random coefficient of the internal damping. This results in an explicit analytical solution for the radius of whirl which may be used for predicting reliability of the rotor. Furthermore, a convenient procedure is described for interpreting measured on-line test data for the rotor. Namely, the mean value of the coefficient of internal damping as well as its standard deviation and mean frequency of temporal variations may be estimated directly from the trace of whirl radius which exhibits spontaneous random outbreaks in response.     

Influence of load torque on stability of rotor driven by flexible shaft

The influence of the load torque on the stability of a symmetric rotor, driven by a flexible shaft, is studied. Both linear and angular displacements of the rotor are considered. The analysis—which is approximate, and with the deflection, the damping and the load torque assumed to be small—shows that the main destabilizing effect of the load torque is due to the transverse moments acting while the rotor is inclined. The reaction pattern at the ends of the shaft, determined by means of the Kirchhoff equations, indicates that the semitangential mode of loading (conservative) is operative.     

Detailed ball bearing model for magnetic suspension auxiliary service

Catcher bearings (CBs) provide backup protection for rotating machines with active magnetic bearings (AMBs). The CBs are required in the event of an AMB failure or high transient loads. Numerical simulations of a rotor drop on CBs in flywheel energy storage system are conducted with a detailed CB model which includes a Hertzian load–deflection relationship between mechanical contacts, speed-and-preload-dependent bearing stiffness due to centrifugal force, and a Palmgren's drag friction torque. The transient simulation results show the rotor shaft response variations with the design parameters: shaft/bearing friction coefficients, axial preload, support damping of damper liner, and side loads from magnetic bearings. The results reveal that friction coefficients, support damping, and side loads are critical parameters to satisfy CB design objectives and prevent backward (super) whirl.     

柔性基础对ICF装置稳定性的影响

 为探讨土-结构相互作用（SSI）对惯性约束聚变（ICF）装置在地脉动作用下响应的影响,建立了某ICF装置刚性地基和柔性地基两种有限元模型。通过设置人工虚拟边界的弹簧单元和阻尼装置来考虑SSI效应,并且通过设置具有概率意义的弹簧刚度和阻尼系数,来模拟岩土特性参数的不确定性。两个模型地脉动随机响应结果的对比分析,揭示了柔性地基条件下SSI效应及岩土地基参数不确定性对装置的影响。研究表明,考虑岩性地基的柔性时,ICF装置主要部位的位移响应和转角响应分别放大了1.24和1.49倍。     

宽弦高速跨音风扇颤振特性研究

颤振是航空发动机、燃气轮机等运行安全的重要威胁,但颤振稳定性与流动结构之间的关系尚不清晰。本文使用行波法和影响系数法,对某宽弦复合掠型高速跨音风扇转子的一阶模态进行了颤振特性研究,计算了在100%转速下从堵塞点到近失速点的颤振表现。使用影响系数法时,分析了不同通道数的计算域对气动阻尼计算的影响,并与行波法得到的结果进行了对比。研究了流动结构与叶片表面气动阻尼之间的关系,旨在提高对流动致颤机理的认识。结果表明影响系数法和行波法均能对叶片的气动阻尼进行较好的预测;流动结构方面,激波、激波附面层分离、叶尖泄漏流以及吸力面前缘叶顶附近的非定常压力波动,对叶片的气动阻尼分布有较大的影响。     

Nonlinear parametric instability of wind turbine wings

Nonlinear rotor dynamic is characterized by parametric excitation of both linear and nonlinear terms caused by centrifugal and Coriolis forces when formulated in a moving frame of reference. Assuming harmonically varying support point motions from the tower, the nonlinear parametric instability of a wind turbine wing has been analysed based on a two-degrees-of-freedom model with one modal coordinate representing the vibrations in the blade direction and the other vibrations in edgewise direction. The functional basis for the eigenmode expansion has been taken as the linear undamped fixed-base eigenmodes. It turns out that the system becomes unstable at certain excitation amplitudes and frequencies. If the ratio between the support point motion and the rotational frequency of the rotor is rational, the response becomes periodic, and Floquet theory may be used to determine instability. In reality the indicated frequency ratio may be irrational in which case the response is shown to be quasi-periodic, rendering the Floquet theory useless. Moreover, as the excitation frequency exceeds the eigenfrequency in the edgewise direction, the response may become chaotic. For this reason stability of the system has in all cases been evaluated based on a Lyapunov exponent approach. Stability boundaries are determined as a function of the amplitude and frequency of the support point motion, the rotational speed, damping ratios and eigenfrequencies in the blade and edgewise directions.     

Dynamic stability of two rigid rotors connected by a flexible coupling with angular misalignment

The effect of misalignment on the stability of two rotors connected by a flexible mechanical coupling subjected to angular misalignment is examined. The study performed is to understand the effect of angular misalignment on the stability of rotating machinery. The dimensionless stability criteria of the non-linear system of differential equations of two misaligned rigid rotors are derived using Liapunov's direct method. A rigid disk is attached at the middle of each rotor, where the rotor-disk assembly is mounted on two hydrodynamic bearings with four stiffness and four damping coefficients. Sets of dimensionless conditions for sufficient whirl stability of the two misaligned rotors are derived. The stability conditions are presented in graphical form for deeper understanding of the effect of the flexible mechanical coupling stiffness and angular misalignment on rotating machinery stability. The results show that an increase in angular misalignment or mechanical coupling stiffness terms leads to an increase of the model stability region.     

NON-LINEAR MODELLING OF ROTOR DYNAMIC SYSTEMS WITH SQUEEZE FILM DAMPERS—AN EFFICIENT INTEGRATED APPROACH

Squeeze film dampers used in rotor assemblies such as aero-engines introduce non-linear damping forces into an otherwise linear rotor dynamic system. The steady state periodic response of such rotor dynamic systems to rotating out-of-balance excitation can be efficiently determined by using periodic solution techniques. Such techniques are essentially faster than time marching techniques. However, the computed periodic solutions need to be tested for stability and recourse to time marching is necessary if no periodic attractor exists. Hence, an efficient integrated approach, as presented in this paper, is necessary. Various techniques have been put forward in order to determine the periodic solutions, each with its own advantages and disadvantages. In this paper, a receptance harmonic balance method is proposed for such a purpose. In this method, the receptance functions of the rotating linear part of the system are used in the non-linear analysis of the complete system. The advantages of this method over current periodic solution techniques are two-fold: it results in a compact model, and the receptance formulation gives the designer the widest possible choice of modelling techniques for the linear part. Stability of these periodic solutions is efficiently tested by applying Floquet Theory to the modal equations of the system and time marching carried out on these equations, when necessary. The application of this integrated approach is illustrated with simulations and an experiment on a test rig. Excellent correlation was achieved between the periodic solution approach and time marching. Good correlation was also achieved with the experiment.     

Swept sine testing of rotor-bearing system for damping estimation

Many types of rotating components commonly operate above the first or second critical speed and they are subjected to run-ups and shutdowns frequently. The present study focuses on developing FRF of rotor bearing systems for damping estimation from swept-sine excitation. The principle of active vibration control states that with increase in angular acceleration, the amplitude of vibration due to unbalance will reduce and the FRF envelope will shift towards the right (or higher frequency). The frequency response function (FRF) estimated by tracking filters or Co-Quad analyzers was proved to induce an error into the FRF estimate. Using Fast Fourier Transform (FFT) algorithm and stationary wavelet transform (SWT) decomposition FRF distortion can be reduced. To obtain a theoretical clarity, the shifting of FRF envelope phenomenon is incorporated into conventional FRF expressions and validation is performed with the FRF estimated using the Fourier Transform approach. The half-power bandwidth method is employed to extract damping ratios from the FRF estimates. While deriving half-power points for both types of responses (acceleration and displacement), damping ratio (ζ) is estimated with different approximations like classical definition (neglecting damping ratio of order higher than 2), third order (neglecting damping ratios with order higher than 4) and exact (no assumptions on damping ratio). The use of stationary wavelet transform to denoise the noise corrupted FRF data is explained. Finally, experiments are performed on a test rotor excited with different sweep rates to estimate the damping ratio.     

The influence of crack breathing and imbalance orientation angle on the characteristics of the critical speed of a cracked rotor

By analyzing the limitations of weight dominance and by taking the complicated whirl of the rotor into account, general equations of motion have been developed in case of a Jeffcott rotor with a transverse crack. The angle between the crack direction and the shaft deformation direction is used to determine the closing and opening of the crack, allowing one to study the dynamic response without assuming weight dominance. Using the new equations, the dynamic response of a cracked rotor near its critical speed has been computed via a numerical method to investigate the influence of nonlinear breathing of the crack and that of the imbalance orientation angle β on the stability, critical speed and peak response of the rotor. The results show that nonlinear breathing can improve the stability of a rotor in contrast to a rotor with an open crack, and, with a reversed imbalance (70°<β<270°), that it can reduce the vibration response in contrast to an uncracked rotor. The basic characteristics of a cracked rotor near its critical speed are similar to those of an uncracked rotor. The critical speed can be determined by measuring the rotation of the center of gravity. The critical speed of a cracked rotor is located between the natural frequencies of the fully open crack and those of the fully closed crack and depends on the imbalance orientation angle. Its value is lowest at β≈90° and highest at β≈270°. The peak in the response at the critical speed is mainly determined by the imbalance orientation angle. At β≈0° and 180°, the peak corresponds to the maximum and minimum response, respectively.     

ADDING ONE BIFURCATION IN A ROTOR SYSTEM WITH CLEARANCES

A mathematical model of a rotor system with clearances is analysed by the application of modern nonlinear dynamic theory. From the bifurcation diagrams, it is discovered that the rotor system alternates between periodic and chaotic motions at a supercritical rotational speed, and after undergoing a chaotic region the periodic number of the motion will increase by one. At the same time, the periodic number is equal correspondingly to the integral multiple of the critical rotational speed. At the subcritical rotational speed, it is discovered that the chaotic bands among successive orders of superharmonic responses return to the period one through a reversed period-doubling bifurcation, as a result of a period-doubling bifurcation. It is shown that the increase of damping may reduce the width of the chaotic bands and the amplitude of the periodic response; the increase of nonlinear degree also leads to the reduction of chaotic bandwidth, but makes the amplitude of the subharmonic response increase. So it is suggested that proper damping and correct material selection by considering the dynamic characteristics of the rotor system may reduce the proportion of faults and enhance the dynamic characteristics when designing the rotor system. The working speed should not be selected at N times its natural frequency and should not be set in the chaotic bands among the successive orders of periodic motion for the same purpose.     

磁流变液阻尼器在转子振动控制中的应用

设计了一种转子振动控制用的剪切式磁流变液阻尼器,建立了磁流变液阻尼器－悬臂转子系统的分析模型,理论和实验研究了转子系统的不平衡响应特性。研究表明,随着施加磁场强度的增加,磁流变液阻尼和刚度增大,转子系统的临界振幅明显下降,系统的临界转速也明显提高。通过简单的开关控制,可抑止转子通过临界转速过程中的振动。     

Whirl and whip instabilities in rotor-bearing system considering a nonlinear force model

Linear models and synchronous response are generally adequate to describe and analyze rotors supported by hydrodynamic bearings. Hence, stiffness and damping coefficients can provide a good model for a wide range of situations. However, in some cases, this approach does not suffice to describe the dynamic behavior of the rotor-bearing system. Moreover, unstable motion occurs due to precessional orbits in the rotor-bearing system. This instability is called “oil whirl” or “oil whip”. The oil whirl phenomenon occurs when the journal bearings are lightly loaded and the shaft is whirling at a frequency close to one-half of rotor angular speed. When the angular speed of the rotor reaches approximately twice the natural frequency (first critical speed), the oil whip phenomenon occurs and remains even if the rotor angular speed increases. Its frequency and vibration mode correspond to the first critical speed. The main purpose of this paper is to validate a complete nonlinear solution to simulate the fluid-induced instability during run-up and run-down. A flexible rotor with a central disk under unbalanced excitation is modeled. A nonlinear hydrodynamic model is considered for short bearing and laminar flow. The effects of unbalance, journal-bearing parameters and rotor arrangement (vertical or horizontal) on the instability threshold are verified. The model simulations are compared with measurements at a real vertical power plant and a horizontal test rig.     

Control of hysteretic instability in rotating machinery by elastic suspension systems subject to dry and viscous friction

Most of the undesired whirling motions of rotating machines can be efficiently reduced by supporting journal boxes elastically and controlling their movement by viscous dampers or by dry friction surfaces normal to the shaft axis, which rub against the frame. In the case of dry dampers, resonance ranges of the floating support configuration can be easily cut off by planning a motionless adhesive state of the friction surfaces. On the contrary, the dry friction contact must change automatically into sliding conditions when the fixed support resonances are to be feared. Moreover, the whirl amplitude can be restrained throughout the speed range by a proper choice of the suspension-to-shaft stiffness ratio and of the support-to-rotor mass ratio.This theoretical research deals firstly with the natural precession speeds and looks for Campbell plots in dependence on the shaft angular speed, for several rotor-suspension systems. Then, the steady response to unbalance is investigated, in terms of rotor and support orbits and of conical path of the rotor axis. In this search, the ranges of adhesive or sliding contact are identified in particular for system with dry friction damping. At last, the destabilizing influence of the shaft hysteresis in the supercritical regime is focalized and the counterbalancing effect of the other dissipative sources is verified. In the nonlinear case of dry friction dampers, the control of linear stability is fulfilled by a perturbation procedure, checking the magnitude of Floquet characteristic multipliers on the complex plane. Moreover, the nonlinear stability far from steady motion is tested by the direct numerical solution of the full motion equations. The comparison configuration of suspension systems with viscous dampers and no dry friction is examined through an analytical first approximation approach and closed-form results for stability thresholds are derived in particular for the symmetric case.     

On the notion of a rotating fluid force induced by swirling flow

The Bently/Muszynska (B/M) model shows that oil whirl and oil whip are both self-sustained vibrations associated with two unstable modes of a rotor–fluid system. The model includes a rotating fluid damping and inertia force. In certain configurations, the rotating damping force overcomes the frictional internal damping of the rotor and pushes the rotor into a stable limit cycle of circular orbiting. Such a notion of a rotating fluid force is based on bulk-flow models of fluid-filled clearances that could be approximated as narrow since the tangential velocity of the fluid then translates to one angular velocity at a certain radial distance defined by an average radius. This paper scrutinizes the assumption of a rotating fluid inertia force and pinpoints the additional inertial effects of the swirling flow as the gap width increases. These effects are clarified by deriving the equation of motion of a body with a mass subjected to motion-induced fluid forces of a confined swirling flow. We show that the inertial effects of the swirling flow counteract the destabilizing effect of the rotating damping force. However, if the body mass is larger than the displaced fluid mass, instability follows. The frequency of the unstable mode is unchanged by the additional inertial effects and is always equal to the frequency of the damping that induces the instability.     

Occurrence of stable periodic modes in a pendulum with cubic damping

Dynamical systems with nonlinear damping show interesting behavior in the periodic and chaotic phases. The Froude pendulum with cubical and linear damping is a paradigm for such a system. In this work the driven Froude pendulum is studied by the harmonic balancing method; the resulting nonlinear response curves are studied further for resonance and stability of symmetric oscillations with relatively low damping. The stability analysis is carried out by transforming the system of equations to the linear Mathieu equation.     

OPTIMUM DESIGN OF A HELICOPTER ROTOR FOR LOW VIBRATION USING AEROELASTIC ANALYSIS AND RESPONSE SURFACE METHODS

An aeroelastic analysis based on finite elements in space and time is used to model the helicopter rotor in forward flight. The rotor blade is represented as an elastic cantilever beam undergoing flap and lag bending, elastic torsion and axial deformations. The objective of the improved design is to reduce vibratory loads at the rotor hub that are the main source of helicopter vibration. Constraints are imposed on aeroelastic stability, and move limits are imposed on the blade elastic stiffness design variables. Using the aeroelastic analysis, response surface approximations are constructed for the objective function (vibratory hub loads). It is found that second order polynomial response surfaces constructed using the central composite design of the theory of design of experiments adequately represents the aeroelastic model in the vicinity of the baseline design. Optimization results show a reduction in the objective function of about 30 per cent. A key accomplishment of this paper is the decoupling of the analysis problem and the optimization problems using response surface methods, which should encourage the use of optimization methods by the helicopter industry.     

基于碰撞能量模型的非线性碰摩评估研究

提出了一种用于评估非线性碰摩的碰撞能量模型(IEM),根据模型特性可得出IEM指标η用于评估碰摩故障的严重程度或发生概率. 通过对在不同阻尼比ζ状态下IEM指标η随转速比Ω的响应变化趋势研究,给出了η随ζ和Ω的响应联合分布. 结合IEM模型特性提出了"敏感区域"的概念. 此区域概念的提出,不仅对于转子设备的操作运行具有较大的参考价值,对于转子系统的设计更有重要的指导意义. 关键词： 碰撞能量模型 碰摩 转子系统 故障诊断     

Error field amplification and rotation damping in tokamak plasmas

Toroidal rotation is normally very weakly damped in plasmas that are magnetically confined in the nominally toroidally symmetric tokamak. However, a strong damping of toroidal rotation is observed as such plasmas approach marginal stability for perturbations that produce a kinklike distortion of the plasma. It is shown that the damping of toroidal rotation by very small departures of the magnetic field from toroidal symmetry is greatly enhanced as marginal stability is approached. The response of a plasma to perturbations is studied using a set of electrical circuit elements, which provide an equation for the rotational damping that requires minimal information about the plasma.     

IDENTIFICATION OF SPEED-DEPENDENT BEARING PARAMETERS

Bearing dynamic characteristics have been a major unknown in the modelling and analysis of large turbo-generators. An identification algorithm for bearing dynamic characterization by using unbalance response measurements is developed for multi-degree-of-freedom (m.d.o.f.) flexible rotor-bearing systems. The algorithm identifies the bearing dynamic parameters, consisting of four effective stiffness and four damping coefficients for each bearing, utilizing frequency domain synchronous unbalance response measurements from the accelerometers attached to the bearing housings in the horizontal and vertical directions, for a minimum two different unbalance configurations. The procedure of identifying bearing dynamic coefficients by using the proposed algorithm is presented and demonstrated through a numerical example. Adding noise to the simulated signal checks the robustness of the algorithm against measurement noise. Combinations of regularization and the generalized singular value decomposition (SVD) are used to tackle an ill-posed problem due to the nearly circular orbit of the rotor at the bearings, as a special case for nearly isotropic bearings. It is demonstrated that by measuring noisy bearing responses with the direction of rotation of the rotor both in the clockwise and anticlockwise directions, the bearing estimation problem for circular orbit becomes well-conditioned. The regularization algorithm is tested for an experimental rotor-bearing rig. The response reproduction capabilities are excellent even in the presence of measurement noise.