Three-Dimensional Thermohydrodynamic Morton Effect Simulation—Part I: Theoretical Model

2014 ◽  
Vol 136 (3) ◽  
Author(s):  
Junho Suh ◽  
Alan Palazzolo

The present study is focused on accurate prediction of the Morton effect problem including journal asymmetric heating and the corresponding long period amplitude oscillations using a nonlinear time transient rotor-dynamic simulation. This paper presents a theoretical model of thermal induced synchronous instability problems in a nonlinear rotor–bearing system, and suggests a new computational algorithm for the nonlinear transient analysis of the Morton effect where the dynamic and thermal problems are combined. For the analysis of the Morton effect problem, a variable viscosity Reynolds equation and a 3D energy equation are coupled via temperature and viscosity, and solved simultaneously. Three-dimensional heat transfer equations of bearing and shaft are modeled by a finite element method, and thermally coupled with the fluid film via a heat flux boundary condition. Asymmetric heat flux into the synchronously whirling rotor is solved by the orbit time averaged heat flux from fluid film to the spinning shaft surface. The journal orbit is calculated by the nonlinear transient dynamic analysis of a rotor–bearing system with a variable time step numerical integration scheme. For the computation time reduction, modal coordinate transformation is adopted in dynamic and thermal transient analysis. Thermal bow effect makes a significant change to the dynamic behavior of a rotor–bearing system, and a thermal hysteresis bode plot, that is one of the characteristics of the Morton effect problem, is presented with time varying spin speed.

2016 ◽  
Vol 139 (1) ◽  
Author(s):  
Xiaomeng Tong ◽  
Alan Palazzolo

This paper performs the parametric studies corresponding with the theoretical Morton effect (ME) model explained in Part I of this paper, where the fully nonlinear transient analysis based on the finite element method is introduced. Operating parameters, such as oil supply temperature, bearing clearance, oil viscosity, etc., are perturbed from the testing conditions to investigate the shifting of critical speeds and ME instability onset speed (IOS). The ME is significantly affected by the rotor bending mode with large overhung deflections, and operating parameters should be adjusted to increase the separation margin between the operating speed and the corresponding critical speed for ME mitigation. Reducing the carryover flow ratio and using the asymmetric bearing pivot offset are capable to suppress the ME by reducing both the average and differential journal temperature. The heat barrier sleeve with air or ceramic isolation is designed to prevent the heat flux into the journal and can successfully mitigate the ME based on the simulations.


Author(s):  
Jianming Cao ◽  
Paul Allaire ◽  
Timothy Dimond

This paper provides a time transient method for solving coupled lateral and torsional analysis of a flexible rotor–bearing system including gyroscopic effects, nonlinear short journal bearings, nonlinear short squeeze film dampers (SFDs), and external nonlinear forces/torques. The rotor is modeled as linear, and the supporting components, including bearings and dampers, are modeled as nonlinear. An implicit Runge–Kutta method is developed to solve the nonlinear equations of motion with nonconstant operating speed since the unbalance force and the gyroscopic effect are related to both the rotational speed and the acceleration. The developed method is compared with a previous torsional analysis first to verify the nonlinear transient solver. Then the coupled lateral and torsional analysis of an example flexible three-disk rotor, perhaps representing a compressor, with nonlinear bearings and nonlinear dampers driven by a synchronous motor is approached. The acceleration effects on lateral and torsional amplitudes of vibration are presented in the analysis. The developed method can be used to study the rotor motion with nonconstant rotational speed such as during startup, shutdown, going through critical speeds, blade loss force, or other sudden loading.


2020 ◽  
Vol 143 (3) ◽  
Author(s):  
Joseph Oh ◽  
Baik Jin Kim ◽  
Alan Palazzolo

Abstract Conventional rotordynamic analyses generally treat the rotor as a continuous body without considering effect of clamped joints. However, in modern rotating machines, rotors are often assembled with multiple complex-shaped parts and joints, which may significantly affect rotordynamic behavior. Several authors have proposed methods for implementing contact effects in rotordynamic analysis, but a more general modeling method for handling arbitrary contact geometries with various levels of surface roughness is needed. The present paper suggests a new contact model for rotordynamic analysis of an assembled rotor-bearing system with multiple parts connected by multiple joints. A contact element formulation is presented using solid finite elements and statistics-based contact theories. A test arrangement was developed to validate the proposed contact model for varying interface surface roughness and preloads. An iterative computation algorithm is introduced to solve the implicit relation between contact stiffness and stress distribution. Prediction results, using the contact model, are compared with measured natural frequencies for multiple configurations of a test rotor assembly. A case study is performed for an overhung type rotor-bearing system to investigate the effect of contact interfaces, between an overhung impeller and a rotor shaft, on critical speeds.


Author(s):  
Yi Liu ◽  
Heng Liu ◽  
Xin Wang ◽  
Minqing Jing

The nonlinear dynamic characteristics of three-dimensional rod-fastening rotor bearing system are investigated in this paper. The rod-fastening rotor includes discontinuous shaft, rotating disks, circumferentially distributed rods, and macrointerfaces between disks. The first three parts are discretized by three dimensional elements, and the macrointerfaces are connected by some springs whose stiffness is determined by a proposed linear partition method. For comparison, the three-dimensional dynamic model of a corresponding complete rotor bearing system is also built. After the rod-fastening and complete rotor bearing system are reduced by a component mode synthesis, periodic motions and stability margins are calculated by using the shooting method and path-following technique, and the local stability of system is obtained by using the Floquet theory. Comparative results show the both systems have a resemblance in the bifurcation features when mass eccentricity and rotating speed are changed. The vibration response has the identical frequency components when typical bifurcations occur. The dynamic stress is obtained by regarding the displacements of all nodes as load. Moreover, the unbalanced and insufficient of the pre-tightening forces lead to obvious disadvantageous influence on the stability and vibration of the both systems. Generally, this paper considers the interfacial effect of the rod-fastening rotor bearing system and the relative nonlinear dynamic features are obtained.


2005 ◽  
Vol 128 (2) ◽  
pp. 252-261 ◽  
Author(s):  
A. Choudhury ◽  
N. Tandon

In the present investigation, a theoretical model has been developed to obtain the vibration response due to a localized defect in various bearing elements in a rotor-bearing system under radial load conditions. The rotor-bearing system has been modeled as a three degrees-of-freedom system. The model predicts significant components at the harmonics of characteristic defect frequency for a defect on the particular bearing element. In the case of a defect on the inner race or a rolling element, the model predicts sidebands about the peaks at defect frequencies, at multiples of shaft and cage frequencies, respectively. The model has also predicted some additional components at harmonics of shaft and cage frequencies due to a local defect on the inner race and a rolling element, respectively. The expressions for all these spectral components have also been derived. Typical numerical results for an NJ 204 bearing have been obtained and plotted. The amplitude of the component at defect frequency, for an outer race defect, is found to be much higher as compared to those due to inner race defect or a rolling element defect of the same size and under similar conditions of load and speed. The results of vibration measurements on roller bearings with simulated local defects have also been presented to experimentally validate the theoretical model proposed. It can be observed from the results that the spectral components predicted by the theoretical model find significant presence in the experimental spectra. Comparison of the normalized analytical values of the spectral components with their experimental values shows fair agreement for most of the cases considered. Probable area of the generated excitation pulses has been calculated and the effects of pulse area variation on the experimental results have been studied.


Author(s):  
An Sung Lee ◽  
Byung Ok Kim

Turbomachinery such as turbines, pumps and compressors, which are installed in transportation systems such as warships, submarines and space vehicles, etc., often perform crucial missions and are exposed to potential dangerous impact environments such as base-transferred shock forces. To protect turbomachinery from excessive shock forces, it may be needed to accurately analyze transient responses of rotors, considering the dynamics of mount designs to be applied with. In this study a generalized FE transient response analysis model, introducing relative displacements, is firstly proposed to accurately predict transient responses of a flexible rotor-bearing system with mount systems to base-transferred shock forces. In the transient analyses the state-space Newmark method of a direct time integration scheme is utilized, which is based on the average velocity concept. Results show that for the identical mount systems considered, the proposed FE-based detailed flexible rotor model yields more reduced transient vibration responses to the same shocks than a conventional simple model or a Jeffcott rotor. Hence, in order to design a rotor-bearing system with a more compact light-weighted mount system, preparing against any potential excessive shock, the proposed FE transient response analysis model herein is recommended.


2009 ◽  
Vol 1158 ◽  
Author(s):  
Bang-Hung Tsao ◽  
Katie Sondergelt ◽  
Jacob Lawson ◽  
James Scofield ◽  
Levi Elston

AbstractA three dimensional thermal-fluid and stress model of a single chip SiC power sub-module was generated using ANSYS in order to determine the maximum temperature and deformation under various conditions. The effects of heat flux, working fluid temperature and differential pressure on temperature and thermal stress contours were of particular concern. Steady state analysis with water as the working fluid, a simulated heat flux of 11.12×104 W/m2, an interface coupling film coefficient of either 30 or 200 W/m2-K between the cooling plate and fluid, and ambient film coefficients from 6 W/m2-K to 300 W/m2-K, predicts maximum device junction temperatures between 374 and 316 K, and corresponding deformations from .0351% to .0293%. Under the same boundary and loading conditions, but with air as the working fluid, the deformations reached .0405% to .0296%, with temperatures between 427 and 316 K. Transient analysis also showed junction temperatures in the predicted range and determined the time to reach steady state to be between 150 and 2500 seconds depending on the boundary conditions. Experiments were conducted in order to validate ANSYS results.


2021 ◽  
Author(s):  
Yuqi Li ◽  
Zhong Luo ◽  
Jinwen Wang ◽  
Hui Ma ◽  
Dongsheng Yang

Abstract Bolted joints are widely used in industrial rotating machines to fasten the adjacent disks together and affect system dynamic properties. Therefore, there is a strong need to study their influence on the response of such systems. This paper investigated the effect of the eccentric phase difference of the disk and bolted-disk joint on rotor dynamics, which takes into account the time-varying bending stiffness of the bolted joint. The bolted-disk joint is modeled as a two-node element based upon the energy theorem and Lagrange’s principle, where the relative lateral displacement stiffness, relative bending stiffness, and coupling stiffness between the adjacent disks are considered. And then the dynamic model of the rotor-bearing system is derived based on the proposed bolted joint element and lumped mass modeling method. Combining with the Newmark-β integration scheme, the established model allows the dynamic response characteristics of the rotor-bearing system with the bolted joint to be predicted, and the impact of eccentric phase difference in the rotor system on the response to be investigated. The validity of the simulation results was confirmed by experiment. Through the modeling method proposed in this paper and obtained results, the bifurcation characteristics of the bolted joint rotor system can be predicted.


2016 ◽  
Vol 139 (1) ◽  
Author(s):  
Yi Liu ◽  
Heng Liu ◽  
Nanshan Wang

The effects of typical machining errors on the dynamic features of rod-fastened rotor bearing system (RBS) are studied in this paper. Three micron-sized machining errors are considered in a three-dimensional (3D) rod-fastened model. The static effects of machining errors are investigated by applying finite element method. Results demonstrate that machining errors not only bring about mass eccentricity but also cause obvious rotor bending due to large pretightening force. Then, nonlinear dynamic features such as stability and bifurcation are analyzed by using target-shooting technique, track-following method, and Floquet theory. Analysis data indicate that rotor bending originated from machining errors reduces the system stability evidently. It is also observed that the vibration value continues to go up after critical speed as rotating speed increases. It is a particular property compared with integral rotor. It explains the reason why the machining precision of rod-fastened rotor is much higher than that of the corresponding integral rotor to some extent. Moreover, differences between machining errors are compared and the results show that the machining precision of axial assembly interfaces should be paid more attention in the rod-fastened rotor design.


Sign in / Sign up

Export Citation Format

Share Document