Effects of a crack on the stability of a non-linear rotor system

The stability of a rotor system presenting a transverse breathing crack is studied by considering the effects of crack depth, crack location and the shaft's rotational speed. The harmonic balance method, in combination with a path-following continuation procedure, is used to calculate the periodic response of a non-linear model of a cracked rotor system. The stability of the rotor's periodic movements is studied in the frequency domain by introducing the effects of a perturbation on the periodic solution for the cracked rotor system.It is shown that the areas of instability increase considerably when the crack deepens, and that the crack's position and depth are the main factors affecting not only the non-linear behaviour of the rotor system but also the different zones of dynamic instability in the periodic solution for the cracked rotor. The effects of some other system parameters (including the disk position and the stiffness of the supports) on the dynamic stability of the non-linear periodic response of the cracked rotor system are also investigated.     

Six degrees of freedom coupled dynamic response of rotor with a transverse breathing crack

More coupling mechanisms between the six degrees of freedom (DOFs) are introduced by considering the contribution of the general transverse forces to stress intensity factor of mode I crack in predicting the crack additional flexibility matrix of the cracked rotor. And the obtained flexibility elements show a good agreement with the experiment result for a wide range of the crack depth ratio. Six DOFs coupled dynamic equations for cracked rotor are formulated by introducing three rotational DOFs. A response-dependent non-linear breathing crack model is applied to simulate the breathing behavior during operation in this paper. Numerical investigations are carried out to simulate the various parametric conditions on the dynamic characteristics of cracked rotor, including the crack depth, shaft slenderness ratio, and rotating speed ratio. A perturbation frequency component and its combinations with harmonic components are observed in the dynamic response obtained by six DOFs coupled models (SDCM). Some differences in evolution of whirling orbit obtained by SDCM and three DOFs coupled models are found, when the cracked rotor passes through the sub-critical speeds.     

A non-linear study of a cracked rotor

The influence of the presence of transverse cracks in a rotating shaft is analyzed. The paper addresses the influence of crack opening and closing on dynamic response during operation. The evolution of the orbit of the cracked rotor near half and one-third of the first critical speed is investigated. The dynamic response of the rotor with a breathing crack is evaluated by expanding the changing stiffness of the crack as a truncated Fourier series and then using the Harmonic Balance Method. This method is applied to compute various parametric studies including the effects of the crack depth and location on the dynamic of a crack rotor. The evolution of the first critical speed, associated amplitudes at the critical speed and half of the critical speed, and the resulting orbits during transient operation are presented and some distinguishing features of a cracked rotor are examined.     

非稳态非线性油膜力作用下柔性轴裂纹转子的动力特性分析

分析了在动载轴承非稳态非线性油膜力作用下，具有横向裂纹柔性轴Jeffcott转子在非线性涡动影响下的动力特性。通过数值计算表明，在油膜失稳转速前，随着裂纹轴刚度变化比的增大，系统在低转速区域内具有丰富的非线性动力行为，出现倍周期分叉及混沌现象，涡动振幅随转速升高而减小，直到非稳态非线性油膜失稳，在无裂纹转子油膜临界失稳点处发现了类Hopf分叉现象，系统运动由平衡变为拟周期运动；裂纹转子在油膜临界失稳时的系统运动亦为拟周期运动，裂纹转子轴刚度变化对油膜失稳点及油膜失稳之后转子的运动影响不大，转子系统作拟周期运动。     

裂纹转子的振动响应研究

本文研究了裂纹转子的振动响应。文章首先通过应力强度因子积分得到含裂纹轴单元的刚度矩阵；建立了裂纹转子的运动微分方程。进而研究了裂纹转子的振动响应，得出了裂纹转子的振动响应随裂纹位置和深度的变化关系。为工程上早期诊断微小裂纹提供了理论根据。     

Parametric instability of a Jeffcott rotor with rotationally asymmetric inertia and transverse crack

Both the rotationally asymmetric inertia and transverse crack frequently appear in the rotor system. The parametric excitations induced by this two features cause instability and severe vibration under certain operating conditions. Thus, the parametric instability of a Jeffcott rotor with asymmetric disk and open transverse crack is studied analytically. The vibration equations of four degrees-of-freedom of the system are established, and the stiffness coefficients of cracked rotor shaft are derived based upon the compliance method and strain energy release rate method. Then, utilizing the harmonic balance method and Taylor expansion technique, the unstable widths of simple and combination instability regions (SIR and CIR) are solved approximately. For a practical rotor system, the approximate unstable widths are verified by the Floquet numerical analysis. The effects of crack depth and position upon the unstable widths are discussed, and the conditions for zero unstable points (ZUPs) are given: Besides the asymmetric angle should be π/2 (for SIR) or 0 (for CIR), the relationships between the inertia asymmetry and crack parameters (depth and position) are also presented analytically. These results would be useful for crack detection and instability control of the asymmetric rotor-bearing system.     

非线性呼吸行为对带横向裂纹Jeffcott转子稳定性的影响

在分析现有力学模型局限性的基础上,考虑转子在失稳转速附近复杂的涡动状态,构建了不依赖于重力占优假设的含横向裂纹的Jeffcott转子的通用运动方程和一种拟合的裂纹呼吸模型。根据Floquet稳定性理论,用打靶法分析了该裂纹转子模型运动的稳定性。结果表明,刚度相对变化量、系统阻尼、裂纹相位角等系统参数对模型的运动稳定性都有较大影响,特别是裂纹相位角对稳定性的影响很复杂;在临界转速附近,裂纹的非线性呼吸行为可以提高转子的稳定性,且改变了失稳发散的行为;在亚临界转速附近,裂纹的非线性呼吸行为会降低转子的稳定性,使周期运动发生结构失稳。     

Analytical study of nonlinear synchronous full annular rub motion of flexible rotor–stator system and its dynamic stability

The nonlinear synchronous full annular rub motion of a flexible rotor induced by the mass unbalance and the contact-rub force with rigid and flexible stator is studied analytically. The nonlinear property is due to the dry friction force between stator and rotor. The exact solutions of the synchronous full annular rub motion and its run speed regions are obtained. The stability of the synchronous full annular rub motion is discussed analytically. The stability criterion and the stability regions of the synchronous full annular rub motion are obtained. A simplified approximate criterion formula for dynamic stability is also derived under the conditions of large impact stiffness, small damping and small friction. The simplified criterion formula can be used conveniently in engineering and matches the real situations of industry.     

Dynamic instability and steady-state response of an elliptical cracked shaft

An elliptical front crack has been found to be more accurate and realistic for modeling the transverse surface crack in rotating machinery compared with the widely used straight front crack. When the shaft rotates, the elliptical crack opens and closes alternatively, due to gravity, and thus, a “breathing effect” occurs. This variance in shaft stiffness is time-periodic, and hence, a parametrically excited system is expected. Therefore, the dynamic instability and steady-state response of a rotating shaft containing an elliptical front crack are studied in the paper. The local flexibility due to the crack is derived, and the governing equations of the crack shaft system are established using the assumed modes method. Utilizing the Bolotin’s method and harmonic balance method, the boundaries of two typical instability regions and maximum response amplitude of the cracked shaft could be computed numerically. The elliptical crack parameters (depth, shape factor and position) and damping are, respectively, considered and discussed for their effects on the dynamic behavior of the elliptical cracked shaft. Some research results might be helpful for the crack detection in rotating machinery.     

任意斜裂纹转子的耦合振动研究

对包含不同类型裂纹(横裂纹、横-斜裂纹以及任意斜裂纹)的转子的耦合振动进行研究,以揭示裂纹转子在不同方向上刚度参数的变化规律及其交叉耦合机理,特别是由此引发的振动特征. 对于包含不同类型裂纹的转子轴段,采用六自由度Timoshenko梁单元模型对其进行单元建模,并基于应变能理论推导计算柔度参数和刚度矩阵. 在此基础上, 采用纽马克-$\beta$数值算法求解裂纹转子的运动方程,获得裂纹转子在单故障或多故障激励(不平衡激励、扭转激励或不平衡激励加扭转激励)作用下的耦合振动响应,进而分析耦合振动谱特征. 与横裂纹和横-斜裂纹相比,任意斜裂纹使转子刚度矩阵的交叉耦合效应更显著,导致转子发生更强烈的弯-扭耦合甚至是纵-弯-扭耦合振动.无论是在不平衡激励还是扭转激励作用下, 弯曲振动与扭转振动幅度都更大. 而且,包含不同类型裂纹的转子的耦合振动特征频率,例如旋转基频与二倍频、扭转激励频率及其边带成分的幅值,对裂纹面方向角具有不同的敏感性. 所得的这些研究结果,可以为转子裂纹的特征参数辨识与诊断提供理论依据.      

A new efficient procedure to solve the nonlinear dynamics of a cracked rotor

The dynamic analysis of a cracked rotor with a breathing crack leads to the formulation of a nonlinear time-dependent problem. For the simple Jeffcott rotor model, this problem has been addressed using numerical integration methods that are very time consuming. A first simplification can be done assuming that the stiffness is a time-dependent function obtained with the quasi-static displacements of the shaft. In this study, we propose a new procedure to analyze the nonlinear dynamic of such a kind of cracked rotors using an iterative technique that transforms the full nonlinear problem in a succession of time-dependent linear ones. We show with different examples that this technique virtually gives the same results as the classical integration methods, but being much more efficient and achieving a significant saving of computation time. The calculations using the proposed method are over a 100 times faster than the corresponding to integrate the full nonlinear problem, being very helpful in on-line crack identification procedures. Also, this analysis shows that, in cases for which the vertical whirl amplitude is greater than the shaft weight static deflection, the use of simplified methods based on the quasi-static stiffness matrix could not be adequate.     

考虑界面滑移效应的组合裂纹梁弯曲解析解

将上部子梁的裂纹等效为线性扭转弹簧,考虑组合梁连接面的滑移位移,建立了以组合裂纹梁挠度和滑移位移为基本未知量的组合裂纹梁弯曲变形一维数学模型．利用Laplace变换及其逆变换,给出了组合裂纹梁弯曲变形一维数学模型的解析通解．在此基础上,研究了均布载荷作用下简支组合裂纹梁的弯曲变形问题,数值分析了连接面剪切刚度、裂纹深度、数目和位置等参数对组合裂纹梁弯曲变形的影响,结果表明：在裂纹处,组合裂纹梁挠度曲线存在尖点,而横截面转角曲线存在跳跃,且随着裂纹数目和深度的增加,挠度和横截面转角跳跃值增大;随着连接面剪切刚度的增加,挠度和横截面转角减小,并最终趋于定值．并且,随着组合梁跨高比的增加,连接面剪切刚度对梁挠度影响逐渐减弱．     

FREE VIBRATION ANALYSIS OF SIMPLY SUPPORTED BEAM WITH BREATHING CRACK USING PERTURBATION METHOD

In this paper, a new analytical method for vibration analysis of a cracked simply supported beam is investigated. By considering a nonlinear model for the fatigue crack, the governing equation of motion of the cracked beam is solved using perturbation method. The solution of the governing equation reveals the superhaxmonics of the fundamental frequency due to the nonlinear effects in the dynamic response of the cracked beam. Furthermore, considering such a solution, an explicit expression is also derived for the system damping changes due to the changes in the crack parameters, geometric dimensions and mechanical properties of the cracked beam. The results show that an increase in the crack severity and approaching the crack location to the middle of the beam increase the system damping. In order to validate the results, changes in the fundamental frequency ratios against the fatigue crack severities are compared with those of experimental results available in the literature. Also, a comparison is made between the free response of the cracked beam with a given crack depth and location obtained by the proposed analytical solution and that of the numerical method. The results of the proposed method agree with the experimental and numerical results.     

考虑密封的悬臂转子系统的裂纹故障分析

以双盘悬臂立式转子-轴承系统为研究对象，建立了系统运动微分方程，并用数值方法分析了在非线性密封力和非线性油膜力作用下的裂纹转子的动力学特性。分析表明，在一定深度裂纹下，转子系统响应随不同角频率比表现出复杂的非线性现象，出现了周期k运动、拟周期运动和混沌运动等多种运动形式。在一定角速度时，工作在远离临界角速度区的转子系统对裂纹非常敏感，而工作在近临界角速度区的转子系统对裂纹不是特别敏感，但是裂纹对它的运动状态影响较大。该研究结果为该类转子-轴承系统的安全运行与故障诊断提供了一定的理论参考。     

On the beneficial effect of rotor suspension anisotropy on viscous-dry hysteretic instability

The destabilizing influence of the internal friction on the supercritical rotor whirl can be efficiently counterbalanced by other external dissipative sources and/or anisotropic suspension systems. The theoretical approach may take the internal dissipation into consideration either by dry or viscous models. Nevertheless, several numerical results and a new perturbation technique of the averaging type prove that similar rotor motions and stability limits are achievable by both models, whence the linear viscous assumption appears preferable. Thus, the internal hysteretic force may be expressed by the product of an equivalent viscous coefficient and the rotor centre velocity relative to a reference frame rotating with the shaft end sections. After calculating the natural frequencies and the response to dynamic imbalances, the stability of the steady motion is checked by the Routh-Hurwitz criterion, focusing the analysis on the individual influence of several characteristic properties, like the gyro structure, the stiffness anisotropy of the supports and their asymmetry, and searching for the external damping level needed for stability. A?fairly interesting result is that the benefit of the suspension anisotropy is most effective for a symmetric rotor mounted at the shaft mid-span and decreases significantly on increasing the configuration asymmetry, even moderately. It is also observed how the stability may somehow be associated with the coupling between progressive and retrograde precession motions.     

基于等效黏性弹簧的黏弹性Timoshenko裂纹梁弯曲解析解

考虑裂纹的缝隙和黏性效应,将梁中横向裂纹等效为黏弹性扭转弹簧,利用广义Delta函数,给出了Laplace变换域内裂纹梁的等效抗弯刚度,得到了具有任意开闭裂纹数目且满足标准线性固体黏弹性本构的Timoshenko梁在时间域内的弯曲变形显式解析通解．在此基础上,通过两个数值算例,分析了时间、梁跨高比和裂纹深度等参数对黏弹性Timoshenko开裂纹梁弯曲变形的影响．结果表明：裂纹黏性对Timoshenko裂纹梁的弯曲具有显著的影响．相比于裂纹的弹性扭转弹簧模型,考虑裂纹黏性效应的黏弹性Timoshenko裂纹梁在裂纹处挠度尖点和转角跳跃现象十分明显．另外,由于横向剪切引起的附加变形,Timoshenko裂纹梁的稳态挠度与Euler-Bernoulli梁挠度的差值为常数,其大小与裂纹模型、梁跨高比或裂纹深度无关,这些结果对梁裂纹无损检测具有指导意义．     

多种非线性力作用下不平衡弹性转子的分岔特性

研究了4自由度不平衡弹性转子在非线性油膜力、非线性内阻力和非线性弹性力联合作用下的动力学特性。结果表明，当只有非线性油膜力作用时，转子只存在由于油膜失稳而导致的倍周期分岔。而当非线性油膜力与非线性内阻力共同作用时，在油膜失稳后，转子产生低频振动。转速继续增加，还会诱发内阻失稳，产生概周期运动。在倍周期分岔中，存在分岔激变现象。本文发现的由于油膜涡动而导致的内阻失稳(概周期运动)是一种未见报道的转子失稳模式(组合失稳)，它与油膜失稳(倍周期运动)一起可作为转子故障诊断的典型失稳模式。     

Nonlinear dynamic response of a functionally graded plate with a through-width surface crack

A nonlinear vibration analysis of a simply supported functionally graded rectangular plate with a through-width surface crack is presented in this paper. The plate is subjected to a transverse excitation force. Material properties are graded in the thickness direction according to exponential distributions. The cracked plate is treated as an assembly of two sub-plates connected by a rotational spring at the cracked section whose stiffness is calculated through stress intensity factor. Based on Reddy’s third-order shear deformation plate theory, the nonlinear governing equations of motion for the FGM plate are derived by using the Hamilton’s principle. The deflection of each sub-plate is assumed to be a combination of the first two mode shape functions with unknown constants to be determined from boundary and compatibility conditions. The Galerkin’s method is then utilized to convert the governing equations to a two-degree-of-freedom nonlinear system including quadratic and cubic nonlinear terms under the external excitation, which is numerically solved to obtain the nonlinear responses of cracked FGM rectangular plates. The influences of material property gradient, crack depth, crack location and plate thickness ratio on the vibration frequencies and transient response of the surface-racked FGM plate are discussed in detail through a parametric study.     

Resonant dynamics in a rotordynamic system with nonlinear inertial coupling and shaft anisotropy

The response of a nonlinear, damped Jeffcott rotor with anisotropic stiffness is considered in the presence of an imbalance. For sufficiently small external torque or large imbalance, resonance capture or rotordynamic stall can occur, whereby the rotational velocity of the shaft is unable to increase beyond the fundamental resonance between the rotational and translational motion. This phenomena provides a mechanism for energy transfer from the rotational to the translational mode. Using the method of averaging a reduced-order model is developed, valid near the resonance, that describes this resonant behavior. The equilibrium points of these averaged equations, which correspond to stationary solutions of the original equations and rotordynamic stall, are described as the applied torque, damping, and anisotropy vary. As the anisotropy increases, assumed to arise from increasing shaft cracks, the torque required to eliminate the possibility of stall increases. However, when the system is started with zero initial conditions, the minimum torque required to pass through the resonance is approximately constant as the anisotropy increases. The predictions from the reduced-order model are verified against numerical simulations of the original equations of motion.     

基于均匀荷载面曲率的简支裂纹梁损伤识别

为了采用模态参数对结构裂纹进行定位与定量,基于集中柔度模型,采用无质量的扭转弹簧模拟裂纹,建立简支裂纹梁的振动微分方程。针对现有柔度曲率指标仅能判断裂纹的大致范围,基于线性插值理论,建立裂纹位置与相邻测点均匀荷载面曲率差的关系,提出裂纹进一步定位公式,实现裂纹位置的精确定位。针对现有大多数损伤识别方法无法实现裂纹的损伤定量,基于位移曲率与结构刚度和弯矩的关系,理论推导了均匀荷载面曲率的结构刚度损伤程度识别方法,基于弹簧串联原理和线刚度思想,首次提出串联等效线刚度模型,建立裂纹深度与均匀荷载面曲率的关系,实现裂纹深度的定量。通过简支裂纹梁数值算例,考虑多裂纹的损伤情况,验证了新方法对裂纹定位与定量的有效性。