YTUP
Journals
About
Services
Guides
Sign InSubmit Article
HomeJournalsSigma Journal of Engineering and Natural Sciences10.62051/ytu.sigma-journal-of-engineering-and-natural-sciences-free-and-forced-whirling-analyses-of-a-single-disk-rotor-subject-to-axial-force
SJSigma Journal of Engineering and Natural Sciences
Get Alerted Download PDF
AbstractKeywords1. Introduction2. Governing Equations And Solution Procedure3. Numerical Results And Discussion4. ConclusionsShare and CiteRelated Articles
Article Open Access1 January 2018

Free and Forced Whirling Analyses of a Single-Disk Rotor Subject to Axial Force

Order Reprints Cite Share

Mohammad Reza ZEINOLABEDINI*, and Mansour RAFEEYAN

* Author to whom correspondence should be addressed.

Sigma Journal of Engineering and Natural Sciences 2018, Vol. 36, Issue 1, pp. 21-31; doi.org/10.62051/ytu.sigma-journal-of-engineering-and-natural-sciences-free-and-forced-whirling-analyses-of-a-single-disk-rotor-subject-to-axial-force

Download PDF

Abstract

Rotor vibrations and its control is an important subject in many industries such as power plants and gas stations. When lateral vibrations of rotors during operation exceed allowable level, a huge damage will be occurred. Surge and stall may be some common reasons of these vibrations. This paper aims to present a simple model for surge and stall and assumes that its effect as a distributed force acting on a rotor-disk system. This is a basic and conceptual model for future investigations in this area. Therefore, the effect of a distributed axial force exerted on an assembled disk on a rotating shaft is investigated theoretically. The set of equations for free and forced whirling analyses of a rotor- disc is derived. Also, transverse load composed of unbalanced mass and total weight of the system is considered. For forced whirling analysis, static deflection of the rotor is considered as the initial conditions and rotational speed of the rotor is considered as a time variable parameter which increases from zero to its nominal value in a limited period of time. For a simply supported rotor, the free whirling analysis is investigated using Galerkin method and using Galerkin and Newmark-beta methods, the forced whirling analysis is studied numerically. Forward and backward frequencies and Campbell diagrams are presented in free whirling analysis and variation of deflection, and shear force in any point of the rotor are depicted versus time in forced whirling analysis. The most advantages of the presented paper are consideration of time-dependency of rotating speed in forced whirling analysis. Results indicate that the axial load has a considerable effect on the forward and backward frequencies and lateral vibration amplitude of the rotor.

Keywords: Axial load; lateral vibration; whirling; campbell diagram.

1. Introduction

Rotor systems are widely used assemblies for power transmission in various kinds of machinery. The lateral vibrations of the rotors must be limited because excessive vibrations will affect the running status of it and may lead to catastrophic accidents. Due to several factors, which contribute to the energy transfer, the rotor rotations generate different modes of vibrations: lateral, torsional and axial modes. Among these modes, the lateral is of the greatest concern [1]. The simplest model that can be adopted to study the flexible behavior of rotors consists of a point mass attached to a massless shaft, often referred to as Jeffcott rotor. One of the most important vibration analysis parameter is a natural frequency. Equality of the system rotation speed and *

Corresponding Author: e-mail: zeinolabedini_mr@yahoo.com, tel: 00983531232498 21

system frequency creates a critical state in the system. Also critical speed is an important parameter in the design of systems and their vibration analysis. The operation of rotating systems in the frequency range of the resonator can be very risky and cause serious damage to it. Parameters like mass and stiffness affect the natural frequency, and any change in the system that changes these two parameters will change the natural frequency of the system. In general, axial forces are indeed small therefore are not considered. In aircraft engines, the thrust generates axial force in the rotor system. In compressors stall and surge phenomena can generate horizontal force on disks. These forces may be constant, while in rotating machinery is specific, they may be harmonic or random. This is the main motivation of the current research. Grybos [2] studied the effect of shear deformations on the critical speeds of rotors. Jei and Lee [3] derived critical speeds of an un-symmetric rotor with rigid disks. Free and forced whirling analysis of viscoelastic rotors was investigated by Sturla and Argento [4]. Jun and Kim [5] studied the effect of torsional torque on whirling of rotors. Tiwari and Bhaduri [6] assessed the effect of boundary conditions on the natural frequencies of the rotor bearing system. Nelson et al. [7] studied the vibration analysis of the Timoshenko rotor with internal damping under axial load. Edney et al. [8] presented the dynamic analysis of the tapered Timoshenko rotor. Chen et al. [9] analyzed an exact and direct modeling for rotor bearing system subject to axial load. Choi et al. [10] presented the consistent derivation of a set of governing differential equations describing the vibration in two orthogonal planes and the torsional vibration of a straight rotor with dissimilar lateral principle moments of inertia, subjected to a constant compressive axial load. Huajian gouyang et al. [11] presented a dynamic model for the vibration of a rotating Timoshenko beam subjected to a three- directional load moving in the axial direction. Askarian et al. [12] assessed the effect of various parameters such as axial force, unbalance, and coupling misalignment on vibration of a rotor. Nawal H.Al et al. [13] derived the equation of motion that governs the transvers vibration of a beam loaded axially and compared the natural frequencies thereof. Motallebi et al. [14] adopted homotopy analysis based method to assess the vibration of a nonlinear beam subject to axial force. Torabi et al. [15] presented an analytical solution for whirling analysis of axial-loaded Timoshenko rotor and corresponding basic function. In this study the dynamics of a rotating shaft, modeled as a Timoshenko rotor, subject to distributed harmonic load acting on the surface of a disk on the rotor are investigated. For this, the dynamic equations of lateral vibrations of the rotor are derived using Timoshenko beam model and Newton’s second law. Then, these governing equations are solved by Galerkin and Newmark method and effect of the symmetric axial force acting on the disk on the Campbell diagram and the rotor vibration response are analyzed.The most advantages of the presented paper are consideration of time-dependency of rotating speed in forced whirling analysis.

2. Governing Equations And Solution Procedure

A simply supported Timoshenko rotor subject to distributed axial force acted on disk is shown in Figure 1. The rotational speed increases from zero to its nominal value ( ) in time as: Ω Ω

1 An element of rotor with forces and bending moment is shown in Figure 2.

(b) Figure 2. An element of rotor a) shear forces and bending moments, b) axial forces

By applying the equilibrium of forces and momentums, the set of equations of motion are extracted as follows: ,

Ω where, , , are displacement and rotation components in and directions, respectively; and ρ is the mass density. A, , are the cross- sectional area, moment of inertia around the x and y axes and polar moment of inertia, respectively, and and are the forces per unit length in x and y directions and axial force and the disk position. According to Timoshenko beam theory, components of bending moment (M) and shear force (F) in x and y directions are presented as follows [16]: ∂ EI ∂z ∂ (3) EI ∂z where E and G are the modulus of elasticity and shear modulus, respectively; k is the shear correction factor depending on the shape of the section and Poisson's material ratio [17]. and applying the following equation for a circular section: By neglecting the term of 2

where, and are the reaction of force in the support and axial force exerted on disk. By re-arranging the equations of motion Eq.5 is yield: ,

By introducing the following complex variables: (8) f and The external forces that consisting weight and unbalanced force written as: Ω

where, is the disk mass, is the unbalanced mass, is the unbalanced radius and θ is the angle of unbalanced mass with respect to horizontal axis, the following set of equations of motion is yield:

2 0 For a rotor with simple support the following equation must be met: 0∶

Inserting Eq.12 into Eq.10 and then running the required simplifications the following equations are obtained:

3. Numerical Results And Discussion

In order to present and analyze the results obtained in this study, some specifications of the subject rotor are selected from Table 1.

Table 1. Considered rotor characteristics characteristics Rotor diameter

3.1. Free whirling

In order to determine the forward and backward rotor frequencies for free vibration, we have the following definitions: ,

In this regard, ω is the natural frequency of the whirling. By applied Galerkin method can be written: (17)

In order to obtain an unbiased zero response in (17), the determinants of the coefficients must be zero; in other words: |

By using equation (18), the natural frequencies will be obtained for the motion of the whirling. Given the coherent coefficients in the stiffness matrix, all the roots of the equation will be purely imaginary, and they will result in forward and backward frequencies. Initially, a rotor without axial force was investigated, the results of which are shown in the form of Campbell diagram for the first and second modes in Figures 3 and 4. As observed in Figures¸ when the rotor is stationary, the forward and backward frequencies are equal and, with increasing rotor speed, the forward frequencies increase and the backward frequencies decrease. The reason for this is the effect of the gyroscopic phenomenon.

Figure 3. The Campbell diagram of the first mode of rotor without axial force

Figure 4. The Campbell diagram of the second mode of rotor without axial force

Then the force of 10000N is applied to the disk. The Campbell diagram for the first and second modes is shown in Figures 5 and 6. The blue charts are related to the rotor without axial force, which is presented for comparison.

Figure 5. The Campbell diagram of the first mode of rotor with a force of 10,000 N

Figure 6. The Campbell diagram of the second mode of rotor with a force of 10,000 N

Applying this force will increase the frequency of the system in the first mode. Also, with increasing force, in the second mode the system frequencies begin to decrease. These trends depend on the position of the disk. There exists an appropriate agreement between results of this study and [12].

3.2. Forced Whirling

In order to examine the accuracy of the drawn up codes, the rotor is first modeled as a simple beam without rotation, and the results are compared with the accurate solution. These results are then obtained by assuming α=1/2, β=1/4 and ∆t=0.01 in Newmark-beta solution and by considering the first five sentences of Galerkin solution (Ni=5) by coding them in the MATLAB software. Given that there exist no reference to verify the results obtained here, the modeling method, drawn up codes, the axial force and the rotor’s rotation are removed, from this introduced system, leaving a beam with a simple support was taken into account. Therefore, the extent of motion, and the shear force along x axis in the middle and left support of the rotor are calculated and the results are shown in Figures 7 -8. 28

The results have reveals that a constant force of approximately 286 N is exerted on the left support. The present solution and analytical solution have been compared in Table 2. Analytical solution in Table 2 has been derived using strength of material formula. There exists an appropriate agreement between results and the drawn up codes of this study and [18]. Table 2. Comparison between present method and the analytical results Parameter (N)

In order to study the effects of axial load on the responses¸ a harmonic axial load on the disk is considered with an amplitude and frequency of 30000 N and 5 rad/s, respectively. The shear force at the left support for the unbalanced rotor with this axial harmonic force is shown in Figures 9-10.

As observed in this set up the value of the shear force in the support is oscillated due to both the unbalanced and horizontal force. As to the unbalanced rotor with axial harmonic force of p=300000sin (5t) the obtained results are shown in Figures 11-12.

As observed there exists a direct relation between horizontal force and direction shear force. Amplitude increasing of direction shear force can be easily seen in these figures. In general, these results indicate which forces are created horizontally the effects of stall and surge can increase the amplitude of rotor lateral vibrations. Consequently, by measuring the vibrations in the bearing of the system, it is possible to detect the occurrence of these phenomena because they can lead to system failure.

4. Conclusions

In this study, the effect of applied horizontal distributed force on a disk of a single disk rotor was studied using Campbell's diagram. The rotor was assumed to be uniform and Timoshenko theory was used. Partial differential equations of motion were derived by considering equilibrium equations for an element of the rotor. The rotor was considered simply supported and GalerkinNewmark method was applied directly to the partial differential equations of motion.The results show that the application of horizontal force in a rotor causes a change in the critical speed of the system. This is due to the effect of horizontal force on the stiffness of the system. The greatest effect is seen in the first mode of the system and increases with the increase of the applied horizontal force. However, with the increase of the force of the second mode, the second mode also changes the critical velocities, which is less than the first mode and decreases. The results of this study show that the effect of this force should be considered in designing rotary systems that are observing the presence of a horizontal force. Also¸ the effects of axial symmetrically force acted on disk on lateral vibration of a single disk rotor was analysed. Responses of lateral vibration for without and with axial symmetrically force were a analysed. From the results it is observed that by increasing the amplitude of force the amplitude of direction shear force increases.

Share and Cite

ZEINOLABEDINI, M.R.; RAFEEYAN, M. Free and Forced Whirling Analyses of a Single-Disk Rotor Subject to Axial Force. Sigma Journal of Engineering and Natural Sciences 2018, Vol. 36, pp. 21-31. https://doi.org/10.62051/ytu.sigma-journal-of-engineering-and-natural-sciences-free-and-forced-whirling-analyses-of-a-single-disk-rotor-subject-to-axial-force

Export:

Related Articles

Examination of strength reduction factor in CFRP-confined columns under axial compressionAysun TEKİN ÖZER, Sema ALACALI, 1 January 2024A Combined Decision Algorithm for Diagnosing Bearing Faults Using Artificial Intelligent TechniquesHüseyin Metin ERTUNÇ, 1 January 2018A Research on the Effect of the Reinforcement Geometry on the Soil Bearing Capacity Under the StripPayam MAJEDI, Semet ÇELİK et al., 1 January 2018An Empirical Approach for Propeller Tip Vortex Cavitation NoiseSavaş SEZEN, Ali DOĞRUL et al., 1 January 2018
Publication History
Published1 January 2018
Versionv1
AccessOpen Access
10.62051/ytu.sigma-journal-of-engineering-and-natural-sciences-free-and-forced-whirling-analyses-of-a-single-disk-rotor-subject-to-axial-force
Article Figures (9)
Figure 1Figure 2Figure 3Figure 4Figure 5Figure 6Figure 7Figure 8Figure 9
Related Articles
Examination of strength reduction factor in CFRP-confined columns under axial compressionAysun TEKİN ÖZER, Sema ALACALISigma Journal of Engineering and Natural Sciences, 1 January 2024A Combined Decision Algorithm for Diagnosing Bearing Faults Using Artificial Intelligent TechniquesHüseyin Metin ERTUNÇSigma Journal of Engineering and Natural Sciences, 1 January 2018A Research on the Effect of the Reinforcement Geometry on the Soil Bearing Capacity Under the StripPayam MAJEDI, Semet ÇELİK et al.Sigma Journal of Engineering and Natural Sciences, 1 January 2018
Sigma Journal of Engineering and Natural Sciences coverSigma Journal of Engineering and Natural Sciences Download PDF

Subscribe to YTUP

Stay connected and receive the latest research updates directly in your inbox.

YTUP — Yıldız Technical University Publishing

Advancing knowledge and fostering innovation through high-quality, peer-reviewed academic publications.

About YTU

Discover

  • ›Articles
  • ›Journals
  • ›Research Topics
  • ›Open Access Policy

Guidelines

  • ›Author guidelines
  • ›Services for authors
  • ›Policies and publication ethics
  • ›Editor guidelines
  • ›Fee policy

Explore

  • ›Articles
  • ›Research Topics
  • ›Journals
  • ›How we publish

Support

  • ›Help center
  • ›Emails and alerts
  • ›Contact us
  • ›Submit
  • ›Career opportunities
YTU Logo

© 2026 Yıldız Technical University (Istanbul, Turkey)

Terms and ConditionsTerms of UsePrivacy PolicyPrivacy SettingsDisclaimer
Like this platform? Join our teamHave feedback or questions?
Supervisor