Investigation on the Steady-State Responses of Asymmetric Rotors

1992 ◽  
Vol 114 (2) ◽  
pp. 194-208 ◽  
Author(s):  
Yuan Kang ◽  
Yuan-Pin Shih ◽  
An-Chen Lee

This paper is to generalize the previous work by utilizing finite element formulation to accommodate the effects of both deviatoric inertia and stiffness due to asymmetry of flexible shaft and disk. A Timoshenko beam element is employed to simulate rotor-bearing systems by taking the gyroscopic moment, rotary inertia, shear deformation of shaft and, asymmetry of disk and shaft into account. Internal damping is not included but the extension is straightforward. Eulerian angles are used to describe the orientations of shaft element and disk, by which, in opposite to the vectorial approach, the mathematical formulation will be symmetric for angular displacements in two directions. The effects of the angle between the major axes of shaft and disk, deviatoric inertia of the asymmetric shaft, and characteristics of bearing on synchronous critical and subcritical speeds are estimated in conjunction with the harmonic balance method. Numerical examples show that the resonant speeds, at which peak responses occur, change due to various angles between major axes, asymmetry of shaft, stiffness, and damping of bearing.

1999 ◽  
Vol 121 (1) ◽  
pp. 78-83 ◽  
Author(s):  
H. Xie ◽  
G. T. Flowers ◽  
L. Feng ◽  
C. Lawrence

This paper investigates the steady-state responses of a rotor system supported by auxiliary bearings in which there is a clearance between the rotor and the inner race of the bearing. A simulation model based upon the rotor of a production jet engine is developed and its steady-state behavior is explored over a wide range of operating conditions for various parametric configurations. Specifically, the influence of rotor imbalance, clearance, support stiffness and damping is studied. Bifurcation diagrams are used as a tool to examine the dynamic behavior of this system as a function of the aforementioned parameters. The harmonic balance method is also employed for synchronous response cases. The observed dynamical responses is discussed and some insights into the behavior of such systems are presented.


Author(s):  
X. F. Wang ◽  
W. D. Zhu

A new spatial and temporal incremental harmonic balanced (STIHB) method is developed for obtaining steady-state responses of a one-dimensional continuous system. In the STIHB method, Galerkin procedure for a governing partial differential equation (PDE) in the spatial coordinate to obtain a set of ordinary differential equations (ODEs) and the harmonic balance procedure for the set of ODEs in the temporal coordinate to obtain the harmonic balanced residual are combined to be Galerkin procedures for the PDE in the spatial and temporal coordinates to simultaneously obtain the spatial and temporal harmonic balanced residual, and integrations in Galerkin procedures are replaced by the fast discrete sine transform (DST) or fast discrete cosine transform (DCT) in the spatial coordinate and the fast Fouriour transform (FFT) in the temporal coordinate, which is referred to as a DST-FFT or DCT-FFT procedure. The harmonic balanced residual for an arbitrary second-order PDE can be automatically and efficiently obtained by a computer program when the expression of the PDE is given, where numbers of basis functions in the spatial and temporal coordinates can be arbitrarily selected and no more extra derivations are needed. There are two versions of the STIHB method. In the simple version, the DST-FFT or DCT-FFT procedure to calculate the harmonic balanced residual and Broyden’s method that is a quasi-Newton method are combined to find solutions that make the residual vanish, which can be used to construct steady-state solutions of the PDE. In the complex version, the exact Jacobian matrix is derived and used in Newton-Raphson method to achieve faster convergence. While its derivation is complex, the exact Jacobian matrix for the arbitrary PDE can be automatically and efficiently obtained by following a calculation routine when the linearized expression of the PDE is given, and it can be easily implemented by a computer program. The exact Jacobian matrix can also be used to study stability of steady-state responses, where no more extra derivations are needed. The STIHB method is demonstrated by studying the transverse vibration of a string with geometric nonlinearity; its frequency-response curves with weak and strong nonlinearities and different numbers of trial functions are calculated, and stability of solutions on the curves is studied.


2017 ◽  
Vol 139 (5) ◽  
Author(s):  
Jocelyn Rebufa ◽  
Fabrice Thouverez ◽  
Erick Le Guyadec ◽  
Denis Mazuyer

A dynamic model of a rotating shaft on two textured hydrodynamic journal bearings is presented. The hydrodynamic mean pressure is computed using multiscale periodic homogenization and is projected on a flexible shaft with internal damping. Harmonic balance method (HBM) is used to study the limit cycles of unbalance response of the coupled system discretized by finite element method (FEM). Stability is analyzed with Floquet multipliers computation. An example of an isotropic texturing pattern representing laser dimples on a lightweight rotor is analyzed. Vibration amplitude and stability zone are compared with plain bearing lubrication. It is shown in an example that full surface texturing leads to relatively higher vibration amplitude compared to plain bearings.


2013 ◽  
Vol 23 (11) ◽  
pp. 1350177 ◽  
Author(s):  
A. Y. T. LEUNG ◽  
H. X. YANG ◽  
P. ZHU

A generalized Duffing–van der Pol oscillator with nonlinear fractional order damping is introduced and investigated by the residue harmonic homotopy. The cubic displacement involved in fractional operator is used to describe the higher-order viscoelastic behavior of materials and of aerodynamic damping. The residue harmonic balance method is employed to analytically generate higher-order approximations for the steady state responses of an autonomous system. Nonlinear dynamic behaviors of the harmonically forced oscillator are further explored by the harmonic balance method along with the polynomial homotopy continuation technique. A parametric investigation is carried out to analyze the effects of fractional order of damping and the effect of the magnitude of imposed excitation on the system using amplitude-frequency curves. Jump avoidance conditions are addressed. Neimark bifurcations are captured to delineate regions of instability. The existence of even harmonics in the Fourier expansions implies symmetry-breaking bifurcation in certain combinations of system parameters. Numerical simulations are given by comparing with analytical solutions for validation purpose. We find that all Neimark bifurcation points in the response diagram always exist along a straight line.


2016 ◽  
Vol 84 (1) ◽  
Author(s):  
X. F. Wang ◽  
W. D. Zhu

A new spatial and temporal harmonic balance (STHB) method is developed for obtaining periodic steady-state responses of a one-dimensional second-order continuous system. The spatial harmonic balance procedure with a series of sine and cosine basis functions can be efficiently conducted by the fast discrete sine and cosine transforms, respectively. The temporal harmonic balance procedure with basis functions of Fourier series can be efficiently conducted by the fast Fourier transform (FFT). In the STHB method, an associated set of ordinary differential equations (ODEs) of a governing partial differential equation (PDE), which is obtained by Galerkin method, does not need to be explicitly derived, and complicated calculation of a nonlinear term in the PDE can be avoided. The residual and the exact Jacobian matrix of an associated set of algebraic equations that are temporal harmonic balanced equations of the ODEs, which are used in Newton–Raphson method to iteratively search a final solution of the PDE, can be directly obtained by STHB procedures for the PDE even if the nonlinear term is included. The relationship of Jacobian matrix and Toeplitz form of the system matrix of the ODEs provides an efficient and convenient way to stability analysis for the STHB method; bifurcations can also be indicated. A complex boundary condition of a string with a spring at the boundary can be handled by the STHB method in combination with the spectral Tau method.


1979 ◽  
Vol 101 (3) ◽  
pp. 263-271 ◽  
Author(s):  
E. H. Chang ◽  
V. K. Garg ◽  
C. H. Goodspeed ◽  
S. P. Singh

A mathematical model for a six-axle locomotive is developed to investigate its dynamic response on tangent track due to vertical and/or lateral track irregularities. The model represents the locomotive as a system of thirty-nine degrees of freedom. The nonlinearities considered in the model are primarily associated with stiffness and damping characteristics of the primary suspension system. The transient and steady-state responses of the locomotive are obtained for the linear and nonlinear primary suspension systems. The response time-histories of the locomotive obtained by integrating the generalized equations of motion are presented. The potential uses of the model are indicated for studying the influence of different design parameters and predicting subsequent dynamic response.


2016 ◽  
Vol 139 (1) ◽  
Author(s):  
X. F. Wang ◽  
W. D. Zhu

A dynamic model of an automotive belt-drive system with a noncircular sprocket instead of a round sprocket is developed in this work to study the effect of reducing the angular variation of camshafts. There are two submodels in the belt-drive system, which are an engine model and a belt-drive model, and they are decoupled to simplify the analysis. When the belt-drive system operates at a steady-state, it is described by a nonlinear model with forced excitation, which can be approximated by a linear model with combined parametric and forced excitations. Steady-state responses of the engine and belt-drive models are calculated by a modified incremental harmonic balance method that incorporates fast Fourier transform and Broyden's method, which is efficient and accurate to obtain a periodic response of a multi-degree-of-freedom system. Steady-state responses of the angular variation of camshafts with different values of sprocket parameters are compared to investigate their optimal values to reduce the angular variation of camshafts. The optimal eccentricity and installation angle are larger and smaller than those from the kinematic model, respectively, which is consistent with published experimental results. This study first shows from a dynamic point of view why use of a noncircular sprocket can reduce the angular variation of camshafts in the operating speed of an engine. Simulation of a speed-up procedure for different sprocket parameters shows results that are consistent with steady-state responses. The belt-drive model developed in this work can be used to select sprocket parameters to minimize the angular variation of camshafts and numerically evaluate the dynamic performance of a belt-drive system with a given design of sprocket parameters.


Author(s):  
X. F. Wang ◽  
W. D. Zhu

A dynamic model of an automotive belt-drive system with a noncircular sprocket instead of a round sprocket is developed in this work to study the effect of reducing the angular variation of camshafts. There are two sub-models in the belt-drive system, which are an engine model and a belt-drive model, and they are decoupled to simplify the analysis. When the belt-drive system operates at a steady state, it is described by a nonlinear model with forced excitation, which can be approximated by a linear model with combined parametric and forced excitations. Steady-state responses of the engine and belt-drive models are calculated by a modified incremental harmonic balanced method that incorporates fast Fourier transform and Broyden’s method, which is efficient and accurate to obtain a periodic response of a multi-degree-of-freedom system. Steady-state responses of the angular variation of camshafts with different values of sprocket parameters are compared to investigate their optimal values to reduce the angular variation of camshafts. Optimal values of sprocket parameters obtained from the dynamic model of the belt drive differ from those from its kinematic model in some previous studies, which is consistent with published experimental results. This study first shows from a dynamic point of view why use of a noncircular sprocket can reduce the angular variation of camshafts in the operating speed of an engine. Simulation of a speed-up procedure for different sprocket parameters shows results that are consistent with steady-state responses.


Sign in / Sign up

Export Citation Format

Share Document