An analytical model for dynamic simulation of the composite rotor with internal damping

Ren, Yongsheng; Zhang, Xingqi; Liu, Yanghang; Chen, Xiulong

Journal of Vibroengineering

Browse Journal

Submit article

Published: 30 December 2014

Check for updates

An analytical model for dynamic simulation of the composite rotor with internal damping

Yongsheng Ren¹

Xingqi Zhang²

Yanghang Liu³

Xiulong Chen⁴

^{1, 2, 3, 4}College of Mechanical and Electronic Engineering, Shandong University of Science and Technology, Qindao, 266590, China

Corresponding Author:

Yongsheng Ren

Cite the article Download PDF

Downloads 1229

WoS Core Citations 1

Abstract

A theoretical model for the dynamics of composite rotor is presented. The composite shaft that includes rigid disks and is supported on rigid bearings is considered as a thin-walled Euler-Bernoulli beam. Internal damping of the composite shaft is taken into account. The equations of motion are derived using the thin-walled composite beam theory based on variational asymptotic method and Hamilton’s principle. The internal damping of shaft is introduced by adopting the multi-scale damping analysis method. Galerkin’s method is used to discretize and solve the governing equations. To demonstrate the validity of the present model, the convergence of the method is examined and the results are compared with those available in the literature. Numerical study shows the effect of design parameters on the natural frequencies, critical rotating speeds and instability thresholds of composite shaft. In addition, the free vibration responses due to the initial perturbations and the forced responses to unbalance for composite shaft are also presented.

1. Introduction

The rotating composite shaft is the primary component for helicopter power transmission and automotive drive system. This can be attributed to the improved properties of composites, such as high strength and stiffness as well as light-weight, compared to metals. On the other hand, composites have excellent damping characteristics, which make them play an important role on damping vibration of laminated plate and shell-like non-rotating structures [1]. Unfortunately, as the rotating speed is in the supercritical range, such damping may cause whirl instabilities of composite shaft. In the field of rotor dynamics, the damping in the composite shaft is called internal damping, which is also called rotating damping. Bucciarelli [2] demonstrated that increasing internal damping may reduce the instability threshold, whereas, the addition of external damping can improve rotor stability. Internal damping characteristics is very important for dynamical stability of composite shaft. Moreover, the destabilizing effect of internal damping can be easily influenced by the arrangement of the different composite layers: fiber orientation angles and number of layers. With the increasing demand for high speed composite rotors, the study of composite shaft based on accurate prediction of internal damping is becoming fundamental as it will provide more reasonable and safe rotor design.

The study of rotating shaft with internal damping was first presented by Kimball [3]. His’s results exhibited that internal damping causes shaft whirling above the first critical speed. Gunter [4] illustrated how above the first critical speed, internal viscous damping leads to instability in the rotating shaft. Vance and Lee [5] studied the stability of high speed rotors with internal viscous damping. The model they used included a single unbalance disk connected with a flexible shaft. The stability was predicted using mathematical methods and verified by numerical solutions of the equation of motion. Melanson and Zu [6] performed vibration analysis of an internally damped rotating shaft, modeled using Timoshenko shaft theory. Explicit analytical expressions for the complex frequencies and mode shapes were presented and the stability threshold was also determined. Montagnier and Hochard [7] studied dynamic instability of supercritical driveshaft with internal damping based on Euler-Bernoulli shaft model. The model taken into account the effect of translational and rotatory inertia. Simon and Flowers [8] presented an investigation of using adaptive control method to suppress instability and synchronous disturbance of a flexible shaft with internal viscous damping. However, all these foregoing literatures deal with only rotating metal shaft.

In recent years, dynamic analysis of the composite shaft with internal damping has also received some attention [9-11]. Singh and Gupta [9] presented the rotordynamic formulations which are applicable to the analysis of the composite shaft with internal damping, obtained by using a layerwise beam theory. Kim et al. [10] conducted the forced vibration analysis of a rotating tapered composite shaft. The stability of a rotating shaft driven through a universal joint has been investigated by Mazzei and Scott [11]. In the above studies, an equivalent viscous damping coefficient was introduced to account for the effect of composite internal damping, but no internal damping analysis has been proposed. To our best knowledge, very few models which incorporate damping estimations are available for predicting the dynamic behavior of composite shaft. Sino et al. [12] presented a dynamic analysis of composite rotating shaft. Internal damping was introduced by the complex constitutive relation of a viscoelastic composite.

The internal damping of composite shaft mainly comes from the dissipation of energy in the materials. Saravanos et al. [13] presented a hollow tubular beam finite element model to predict the damping properties of anisotropic thin-walled closed-section structures. However, this model is only applicable to non-rotating composite beams or blades.

An analytical model of rotating composite thin-walled shaft including internal damping is proposed. This model is based on the composite thin-walled beam theory, referred to as variational asymptotically method (VAM) by Berdichevsky et al. [14]. The internal damping of shaft is introduced via the multi-scale damping mechanics [15]. The flexible composite shaft is assumed supporting on on rigid bearings and containing of the rigid disks mounted on it. The equations of motion of the composite shaft are derived by Hamilton’s principle. Galerkin’s method is used then to discretize and solve the governing equations. The natural frequencies, critical rotating speeds, instability thresholds are obtained through numerical simulations. The effect of the ply angle is then assessed. The validity of the model is proved by comparing the results with those in literatures and convergence examination. Finally, the transient free vibration responses due to the initial perturbations and the steady state forced responses to unbalance are also obtained by using the time-integration method and the mode superposition method, respectively.

2. Composite shaft

The slender thin-walled composite shaft subjected to a rotation along its longitudinal $x$ -axis at a constant rotating speed $Ω$ is shown in Fig. 1. The length, wall thickness and radius of the shaft are denoted by $L$ , $h$ and $r$ , respectively. To describe the motion of the shaft the following coordinate systems are used: (I) $(X, Y, Z)$ is inertial reference system whose origin $O$ is located in the geometric center, (II) $(x, y, z)$ is rotating reference system with the common origin $Ο$ . $(I, J, K)$ and $(i, j, k)$ denote the unit vectors of the reference system $(X, Y, Z)$ and $(x, y, z)$ , respectively. In addition, a local coordinate system $(ς, s, x)$ is used, where $ς$ and $s$ are measured along the direction normal and tangent to the middle surface, respectively.

The variation of the strain energy of the cross-section can be expressed as:

1

δ U_{s} = \int_{A}^{} δ ε_{x, y}^{T} {\bar{Q}}_{i j} ε_{x, y} d s d ς,

where $ε_{x, y}^{T} = {\begin{matrix} γ_{11} & γ_{12} & {2 γ}_{12} \end{matrix}}$ , $γ_{α β}$ $(α, β = 1, 2)$ are the in-plane tensorial strain components, $A$ is the cross-sectional area of the shaft and ${\bar{Q}}_{i j}$ equivalent off-axis stiffness matrix of a composite ply with respect to the system $ο s x ς$ .

The variation of the dissipated energy of the cross-section can be expressed as:

2

δ W_{s} = \int_{A}^{} δ ε_{x, y}^{T} {\bar{η}}_{i j} {\bar{Q}}_{i j} ε_{x, y} d s d ς,

where ${\bar{η}}_{i j}$ is equivalent off-axis damping matrix of a composite ply with respect to the system $ο s x ς$ .

Fig. 1Geometry and coordinate systems of composite thin-walled shaft

Further, the variation of the kinetic energy of the cross-section with rotating motion is:

3

δ T_{s} = \int_{A}^{} δ R^{T} d i a g (ρ) \ddot{R} d s d ς,

where $d i a g (ρ)$ is a diagonal matrix with components equal to the mass density $ρ$ of a ply.

The position, velocity and acceleration vectors for the deformed shaft are described as:

4

R = (y + u_{2}) i + (z + u_{3}) j + (x + u_{1}) k,

\dot{R} = [{\dot{u}}_{2} - Ω (z + u_{3})] i + [{\dot{u}}_{3} + Ω (y + u_{2})] j + {\dot{u}}_{1} k,

\ddot{R} = [{\ddot{u}}_{2} - 2 Ω {\dot{u}}_{3} - Ω^{2} (y + u_{2})] i + [{\ddot{u}}_{3} + 2 Ω {\dot{u}}_{2} - Ω^{2} (z + u_{3})] j + {\ddot{u}}_{1} k

where $u_{1}$ , $u_{2}$ and $u_{3}$ are the displacements of any point on the the cross-section of the shaft in the $x$ , $y$ and $z$ directions, dot indicates the time derivative of the variable.

2.1. Equivalent cross-section stiffness matrix

To derive the strain energy and dissipated energy of the cross-section, the composite thin-walled beam theory proposed by Berdichevsky et al. [14] is employed in the following.

For the case of no internal pressure acting on the shaft, Eq. (1) can be simplified by using free hoop stress resultant assumption as:

5

δ U_{s} = \frac{1}{2} \oint {\begin{matrix} δ γ_{11} & δ γ_{12} \end{matrix}} [\begin{matrix} A & B \\ B & C \end{matrix}] \{\begin{matrix} γ_{11} \\ γ_{12} \end{matrix}\} d s,

where $\oint (\cdot)$ denotes the integral around the loop of the mid-line cross-section, and the reduced axial, coupling and shear stiffness $A$ , $B$ and $C$ can be written as:

6

A (s) = A_{11} - \frac{A_{12}^{2}}{A_{22}}, B (s) = 2 [A_{16} - \frac{A_{12} A_{26}}{A_{22}}],

C (s) = 4 [A_{66} - \frac{A_{26}^{2}}{A_{22}}], A_{i j} = \sum_{k = 1}^{N} {\bar{Q}}_{i j}^{(k)} (z_{k} - z_{k - 1}), i, j = 1, 2, 6 .

In Eq. (5), $δ U_{s}$ can also be expressed with respect to the $u$ , $v$ , $w$ and $φ$ as:

7

δ U_{S} = δ Δ^{T} K Δ,

where $Δ$ is 4×1 column matrix of kinematic variables defined as $Δ^{T} = (\begin{matrix} u^{'} & φ^{'} & w^{''} & v^{''}) \end{matrix}$ , and $K$ is 4×4 symmetric stiffness matrix. Its components $k_{i j}$ are given in Ref. [14].

2.2. Equivalent cross-section damping matrix

Similar to the derivation of the previous cross-section stiffness formulations, the variation of the dissipated energy of the cross-section [15]:

8

δ W_{s} = \frac{1}{2} \oint {\begin{matrix} δ γ_{11} & δ γ_{12} \end{matrix}} [\begin{matrix} A_{d} & B_{d} \\ B_{d} & C_{d} \end{matrix}] \{\begin{matrix} γ_{11} \\ γ_{12} \end{matrix}\} d s,

where:

9

A_{d} = A_{d 11} + (\frac{A_{12}^{2}}{A_{22}^{2}}) A_{d 22} - (\frac{A_{12}}{A_{22}}) (A_{d 12} + A_{d 21}),

B_{d} = A_{d 16} + A_{d 61} + 2 (\frac{A_{12} A_{26}}{A_{22}^{2}}) A_{d 22} - (\frac{A_{26}}{A_{22}}) (A_{d 12} + A_{d 21}) - (\frac{A_{12}}{A_{22}}) (A_{d 26} + A_{d 62}),

C_{d} = 4 [A_{d 66} + (\frac{A_{26}^{2}}{A_{22}^{2}}) A_{d 22} - (\frac{A_{26}}{A_{22}}) (A_{d 26} + A_{d 62})],

A_{d i j} = \int_{- h / 2}^{h / 2} {\bar{ψ}}_{i l} {\bar{Q}}_{l j} d ς = 2 \sum_{k = 1}^{N / 2} {\bar{ψ}}_{i l}^{k} {\bar{Q}}_{l j}^{k} (h_{k} - h_{k - 1}), i, j, l = 1, 2, 6 .

The variation of the dissipated energy also be expressed in terms of the kinematic variables as:

10

δ W_{s} = δ Δ^{T} C Δ,

where $C$ is 4×4 symmetric damping matrix.

The formulation of its components $c_{i j}$ is analogous to stiffness components $k_{i j}$ , but the terms $A$ , $B$ and $C$ in Eq. (6) should be replaced by the terms $A_{d}$ , $B_{d}$ and $C_{d}$ in Eq. (9), respectively.

2.3. Equivalent cross-section damping matrix

Based on the displacement expressions presented in Ref. [14] and in view of Eqs. (3) and (4), the variation of the kinetic energy of the cross-section can be obtained as:

11

δ T_{s} = - (I_{1} δ u + I_{2} δ v + I_{3} δ w {+ I}_{4} δ φ),

where:

12

I_{1} = b_{1} \ddot{u}, I_{2} = b_{1} (\ddot{v} - 2 Ω \dot{w} - Ω^{2} v) - b_{2} (2 Ω \dot{φ} + Ω^{2}) - b_{3} (\ddot{φ} - Ω^{2} φ),

I_{3} = b_{1} (\ddot{w} + 2 Ω \dot{v} - Ω^{2} w) + b_{2} (\ddot{φ} - Ω^{2} φ) - b_{3} (2 Ω \dot{φ} + Ω^{2}),

I_{4} = b_{2} (\ddot{w} + 2 Ω \dot{v} - Ω^{2} w) - b_{3} (\ddot{v} - 2 Ω \dot{w} - Ω^{2} v) + (b_{4} + b_{5}) (\ddot{φ} - Ω^{2} φ),

b_{1} = \int_{A}^{} ρ d s d ς, b_{2} = \int_{A}^{} ρ y d s d ς, b_{3} = \int_{A}^{} ρ z d s d ς,

b_{4} = \int_{A}^{} ρ y^{2} d s d ς, b_{5} = \int_{A}^{} ρ z^{2} d s d ς .

3. Rigid disks

3.1. Kinetic energy of rigid disks

According to Eqs. (11) and (12) the expression of variation of the kinetic energy of the rigid disks fixed to the shaft is written as:

13

δ T_{D} = \sum_{l = 1}^{N_{D}} (I_{1 l}^{D} δ u + I_{2 l}^{D} δ v + I_{3 l}^{D} δ w + I_{4 l}^{D} δ φ) Δ (x - x_{D l}),

in which:

14

I_{1 l}^{D} = b_{1 l}^{D} \ddot{u,} I_{2 l}^{D} = b_{1 l}^{D} (\ddot{v} - 2 Ω \dot{w} - Ω^{2} v) - b_{2 l}^{D} (2 Ω \dot{φ} + Ω^{2}) - b_{3 l}^{D} (\ddot{φ} - Ω^{2} φ),

I_{3 l}^{D} = b_{1 l}^{D} (\ddot{w} + 2 Ω \dot{v} - Ω^{2} w) + b_{2 l}^{D} (\ddot{φ} - Ω^{2} φ) - b_{3 l}^{D} (2 Ω \dot{φ} + Ω^{2}),

I_{4 l}^{D} = b_{2 l}^{D} (\ddot{w} + 2 Ω \dot{v} - Ω^{2} w) - b_{3 l}^{D} (\ddot{v} - 2 Ω \dot{w} - Ω^{2} v) + (b_{4 l}^{D} + b_{5 l}^{D}) (\ddot{φ} - Ω^{2} φ),

b_{1 l}^{D} = \int_{A}^{} ρ_{l}^{D} d s d ς, b_{2 l}^{D} = \int_{A}^{} ρ_{l}^{D} y d s d ς, b_{3 l}^{D} = \int_{A}^{} ρ_{l}^{D} z d s d ς,

b_{4 l}^{D} = \int_{A}^{} ρ_{l}^{D} y^{2} d s d ς, b_{5 l}^{D} = \int_{A}^{} ρ_{l}^{D} z^{2} d s d ς,

where, $ρ_{l}^{D}$ denotes the mass density of the disk $l$ . The symbol $Δ (x - x_{D l})$ denotes Dirac delta function, $N_{D}$ the number of the disks mounted on the shaft and $x_{D l}$ the location of the $l$ th disk.

3.2. Work of external forces

The external forces include the centrifugal force resulting from the unbalanced rigid disks. The virtual work done by the centrifugal force of rigid disks can be given by:

15

δ W_{e} = p_{x}^{u} δ u + p_{v}^{u} δ v + p_{z}^{u} δ w + G_{x}^{u} δ φ,

where, $p_{x}^{u}$ , $p_{y}^{u}$ and $p_{z}^{u}$ denote centrifugal force per unit length, $G_{x}^{u}$ denotes centrifugal moment per unit length. They are given by:

16

p_{x}^{u} = 0, p_{y}^{u} = Ω^{2} \cos Ω t \sum_{l = 1}^{N_{D}} b_{1 l}^{D} e_{l} Δ (x - x_{D l}),

p_{z}^{u} = Ω^{2} \sin Ω t \sum_{l = 1}^{N_{D}} b_{1 l}^{D} e_{l} Δ (x - x_{D l}), G_{x}^{u} = 0,

where, the mass unbalance associated with disks is assumed as the concentrated masses $b_{1 l}^{D} = \iint_{A}^{} ρ_{l}^{D} d A$ at points with small distances of eccentricity $e_{l}$ .

4. Equations of motion and approximate solution method

The equations of motion of the composite rotor can be derived by employing the following Hamilton’s principle:

17

\int_{t_{0}}^{t_{1}} \int_{0}^{L} (δ U_{s} - δ T_{s} - δ T_{D} + δ W_{s} - δ W_{e}) d x d t = 0 .

In order to find the approximate solution, the quantities $u (x, t)$ , $v (x, t)$ , $w (x, t)$ and $φ (x, t)$ are assumed in the form:

18

u (x, t) = \sum_{j = 1}^{N} X_{j u} (t) α_{j} (x), φ (x, t) = \sum_{j = 1}^{N} X_{j φ} (t) θ_{j} (x),

v (x, t) = \sum_{j = 1}^{N} X_{j v} (t) ψ_{j} (x), w (x, t) = \sum_{j = 1}^{N} X_{j w} (t) ψ_{j} (x),

where $α_{j} (x)$ , $θ_{j} (x)$ and $ψ_{j} (x)$ are mode shape functions which fulfill all the boundary conditions of the composite shaft, $X_{j u} (t)$ , $X_{j φ} (t)$ , $X_{j v} (t)$ and $X_{j w} (t)$ are the generalized coordinates.

Substituting Eq. (18) into the governing equations of motion, and applying Galerkin’s procedure, the following governing equations in matrix form can be found:

19

M \ddot{X} + (C + G) \dot{X} + (K + H) X = f,

where $X^{T} = (X_{1 u} (t), \dots, X_{N u} (t), X_{1 φ} (t), \dots, X_{N φ} (t), X_{1 w} (t), \dots, X_{N w} (t), X_{1 v} (t), \dots, X_{N v} (t)),$ $M$ is the mass matrix, $C$ is the damping matrix, $G$ is the gyroscopic matrix, $K$ the stiffness matrix. $H$ is the rotating softening matrix, $f$ is the generalized force vector. Their detailed expressions are follows:

20

M = [\begin{matrix} - b_{1} H_{i j} & 0 & 0 & 0 \\ 0 & - (b_{4} + b_{5}) L_{i j} & - b_{2} M_{i j} & b_{3} M_{i j} \\ 0 & b_{2} Q_{i j} & b_{1} R_{i j} + b_{5} S_{i j} + \sum_{i = 1}^{N_{D}} b_{1 l}^{D} R_{i j}^{D l} + \sum_{i = 1}^{N_{D}} b_{5 l}^{D} S_{i j}^{D i} & 0 \\ 0 & - b_{3} Q_{i j} & 0 & b_{1} R_{i j} + b_{4} S_{i j} + \sum_{i = 1}^{N_{D}} b_{1 l}^{D} R_{i j}^{D l} + \sum_{i = 1}^{N_{D}} b_{4 l}^{D} S_{i j}^{D l} \end{matrix}],

21

C = [\begin{matrix} c_{11} E_{i j} & c_{12} F_{i j} & c_{13} G_{i j} & c_{14} G_{i j} \\ c_{12} I_{i j} & c_{22} J_{i j} & c_{23} K_{i j} & c_{24} K_{i j} \\ c_{13} N_{i j} & c_{23} O_{i j} & c_{33} P_{i j} & c_{34} P_{i j} \\ c_{14} N_{i j} & c_{24} O_{i j} & c_{34} P_{i j} & c_{44} P_{i j} \end{matrix}]

22

G = [\begin{matrix} 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 2 Ω b_{1} R_{i j} + 2 Ω \sum_{i = 1}^{N_{D}} b_{1 l}^{D} R_{i j}^{D l} \\ 0 & 0 & - 2 Ω b_{1} R_{i j} - 2 Ω \sum_{i = 1}^{N_{D}} b_{1 l}^{D} R_{i j}^{D l} & 0 \end{matrix}],

23

K = [\begin{matrix} k_{11} E_{i j} & k_{12} F_{i j} & k_{13} G_{i j} & k_{14} G_{i j} \\ k_{12} I_{i j} & k_{22} J_{i j} & k_{23} K_{i j} & k_{24} K_{i j} \\ k_{13} N_{i j} & k_{23} O_{i j} & k_{33} P_{i j} & k_{34} P_{i j} \\ k_{14} N_{i j} & k_{24} O_{i j} & k_{34} P_{i j} & k_{44} P_{i j} \end{matrix}],

24

H = [\begin{matrix} 0 & 0 & 0 & 0 \\ 0 & (b_{4} + b_{5}) L_{i j} Ω^{2} & 0 & 0 \\ 0 & 0 & - b_{1} R_{i j} Ω^{2} - Ω^{2} \sum_{l = 1}^{N_{D}} b_{1 l}^{D} R_{i j}^{D l} & 0 \\ 0 & 0 & 0 & - b_{1} R_{i j} Ω^{2} - Ω^{2} \sum_{l = 1}^{N_{D}} b_{1 l}^{D} R_{i j}^{D l} \end{matrix}],

25

f^{T} = (0, \dots, 0,0, \dots, 0, Ω^{2} \sum_{l = 1}^{N_{D}} b_{1 l}^{D} e_{l} B_{1}^{D l} \sin Ω t, \dots, Ω^{2} \sum_{l = 1}^{N_{D}} b_{1 l}^{D} e_{l} B_{N}^{D l} \sin Ω t, Ω^{2} \sum_{l = 1}^{N_{D}} b_{1 l}^{D} e_{l} B_{1}^{D l} \cos Ω t, \dots, Ω^{2} \sum_{l = 1}^{N_{D}} b_{1 l}^{D} e_{l} B_{N}^{D l} \cos Ω t),

where:

26

E_{i j} = \int_{0}^{L} α_{i} α_{j}^{''} d x, F_{i j} = \int_{0}^{L} α_{i} θ_{j}^{''} d x, G_{i j} = \int_{0}^{L} α_{i} ψ_{j}^{'''} d x, H_{i j} = \int_{0}^{L} α_{i} α_{j} d x,

I_{i j} = \int_{0}^{L} θ_{i} α_{j}^{''} d x, J_{i j} = \int_{0}^{L} θ_{i} θ_{j}^{''} d x, K_{i j} = \int_{0}^{L} θ_{i} ψ_{j}^{'''} d x, L_{i j} = \int_{0}^{L} θ_{i} θ_{j} d x,

M_{i j} = \int_{0}^{L} θ_{i} ψ_{j} d x, N_{i j} = \int_{0}^{L} ψ_{i} α_{j}^{'''} d x, O_{i j} = \int_{0}^{L} ψ_{i} θ_{j}^{'''} d x, P_{i j} = \int_{0}^{L} ψ_{i} ψ_{j}^{''''} d x,

Q_{i j} = \int_{0}^{L} ψ_{i} θ_{j} d x, R_{i j} = \int_{0}^{L} ψ_{i} ψ_{j} d x, S_{i j} = \int_{0}^{L} ψ_{i} ψ_{j}^{''} d x, R_{i j}^{D l} = \int_{0}^{L} ψ_{i} ψ_{j} Δ (x - x_{D l}) d x,

S_{i j}^{D l} = \int_{0}^{L} ψ_{i} ψ_{j}^{''} Δ (x - x_{D l}) d x, B_{i}^{D l} = \int_{0}^{L} ψ_{i} Δ (x - x_{D l}) d x, i, j = 1, \dots, N .

The solution for free vibration of Eq. (19) can be written in the form $\{X\} = {X_{0}} e^{λ t}$ , from Eq. (19), one can obtain the following characteristic equation:

27

\det (- λ^{2} M + i λ (G + C) + K + H) = 0

Complex eigenvalue $λ$ can be expressed in the form:

28

λ = σ + i ω .

The damping natural frequency or whirl frequency of the system is the imaginary part $ω$ , whereas its real part $σ$ gives the decay or growth of the amplitude of vibration. A negative value of $σ$ indicates a stable motion, whereas a positive value indicates an unstable motion, growing exponentially in time.

It should be pointed out that the govering Eq. (19) is derived based on the rotating frame $O x y z$ . In addition, for a CUS (Circumferentially Uniform Stiffness) configuration [16], Eq. (19) involving in terms of displacements can be split into two independent equation systems associated with both flapping bending-sweeping bending and extension-twist motions. The latter is important to helicopter blades. However, we will focus on the flapping bending-sweeping bending coupling only and give all numerical examples with reference to the inertial frame $O X Y Z$ in this paper, the resulting govering equation are not shown for the sake of simplicity.

5. Numerical results

5.1. Model verifications

In order to examine the influence of the number of mode shape functions used in the solution of the equation on the accuracy of the results, the numerical calculations are performed by considering the shaft made of graphite-epoxy whose elastic characteristics are listed in Table 1. The shaft has rectangular cross-section of fixed width $a =$ 0.32 m, wall thickness $h =$ 0.01016 m and length $L =$ 4.5952 m, whose layup is ${[θ]}_{16}$ with clamped-free boundary conditions.

Table 1Mechanical properties of composite material [13]

$ρ$ (kg/m³)	$E_{11}$ (GPa)	$E_{22}$ (GPa)	$G_{12}$ (GPa)	$υ_{12}$	$η_{l 1}$ (%)	$η_{l 2}$ (%)	$η_{l 6}$ (%)
1672	25.8	8.7	3.5	0.34	0.65	2.34	2.89

The numerical results of natural frequencies and modal dampings are shown in Tables 2 and 3 for an increasing number of mode shape functions. From these tables, it can be seen that to obtain the accurate results of the first two natural frequencies and dampings of flapping, sweeping and torsional mode no more than five mode shape functions are required. This indicates clearly that the convergence of the present model is quite good.

A comparison of predictions using the present model with those obtained in Ref. [13] are also shown in Table 2 and Table 3. A perfect agreement of numerical results with those in Ref. [13] can be seen.

Table 2Modal frequencies and damping of cantilever composite box beam: L/a= 14.36, a/b= 5, [90]16

Mode	Natural frequency (Hz)				Damping (%)
	Present			Ref. [13]	Present			Ref. [13]
	$N =$ 1	$N =$ 3	$N =$ 5	Ref. [13]	$N =$ 1	$N =$ 3	$N =$ 5	Ref. [13]
First flapping	1.80	1.81	1.81	1.80	2.39	2.38	2.38	2.35
Second flapping	–	11.36	11.36	11.5	–	2.37	2.37	2.35
First sweeping	5.67	6.20	6.20	6.50	3.03	2.53	2.53	2.35
Second sweeping	–	39.21	39.21	39.7	–	2.36	2.36	2.37
First torsional	37.83	37.84	37.84	37.7	2.89	2.89	2.89	2.89
Second torsional	112.90	113.01	113.01	113.3	2.82	2.92	2.92	2.89

Table 3Modal frequencies and damping of cantilever composite box beam: L/a= 14.36, a/b= 5, [0]16

Mode	Natural frequency (Hz)				Damping (%)
	Present			Ref. [13]	Present			Ref. [13]
	$N =$ 1	$N =$ 3	$N =$ 5	Ref. [13]	$N =$ 1	$N =$ 3	$N =$ 5	Ref. [13]
First flapping	3.12	3.13	3.13	3.10	0.64	0.65	0.65	0.65
Second flapping	–	19.02	19.02	19.8	–	0.69	0.69	0.67
First sweeping	10.80	10.81	10.81	11.0	0.68	0.68	0.68	0.68
Second sweeping	–	68.38	68.38	65.6	–	0.85	0.85	0.90
First torsional	37.83	37.84	37.84	37.7	2.89	2.89	2.89	2.89
Second torsional	112.91	113.01	113.01	113.3	2.33	2.92	2.92	2.89

Table 4Comparison of the non-dimensional natural frequencies of a rotating cantilever composite shaft

$Ω^{*}$	$ω_{1}^{*}$	$ω_{2}^{*}$	$ω_{3}^{*}$	$ω_{4}^{*}$	$ω_{5}^{*}$	$ω_{6}^{*}$
0	3.5160	3.5160	22.0345	22.0345	61.6973	61.6973
Ref. [17]	3.5160	3.5160	22.0340	22.0340	61.6970	61.6970
2	1.5160	5.5160	20.0345	24.0345	59.6973	69.6973
Ref. [17]	1.5160	5.5160	20.0340	24.0340	59.6970	63.6970
3.5	0.0160	7.0160	18.5345	25.5345	58.1973	65.1973
Ref. [17]	0.0000	7.0160	18.5340	25.5340	58.1970	65.1970
4	–	7.5160	18.0345	26.0345	57.6973	65.6973
Ref. [17]	–	7.5160	18.0340	26.0340	57.6970	65.6970
8	–	11.5160	14.0345	3.30345	53.6973	69.6973
Ref. [17]	–	11.5160	14.0340	30.0340	53.6970	69.6970

The non-dimensional natural frequencies of a rotating cantilever composite shaft obtained using the present model together with those obtained in Ref. [17] are shown in Table 4 for different non-dimensional rotating speeds. A perfect agreement of numerical results with those in Ref. [17] can be seen.

5.2. Results and discussion

In the following numerical example, the composite rotor system is a composite shaft with one rigid disk supported by two rigid bearings at the ends, the disk is at the mid-point of the shaft. The material properties of the composite shaft are those in Table 1. The stacking sequence is ${[\pm θ]}_{5}$ . The length, the radius and the wall thickness of the shaft are assumed to be $L =$ 2.47 m, $r =$ 0.0135 m and $h =$ 0.001321 m, respectively whereas those of a uniform rigid disk are $b_{1 l}^{D} =$ 2.4364 kg, $e =$ 5×10^-5 m, $b_{4 l}^{D} = b_{5 l}^{D} =$ 0.1901 kg/m². The rotor system is shown in Fig. 2.

Fig. 2A composite shaft-rigid disk with rigid bearings

Fig. 3 shows the variation of the first two flexural natural frequencies vs. rotating speed for various ply angles, where, 1F denotes the first forward mode whereas 2F denotes the second forward mode. The synchronous whirl line $ω = Ω$ is also presented in the same figure which is the well-known Campbell diagram. The intersection point of the line $ω = Ω$ and the flexural natural frequency curves is referred to as the critical speed. The motion corresponding to positive flexural natural frequencies are in forward mode. The backward mode branch of the Campbell diagram with negative value of frequencies is on the lower half $ω Ω$ -plane and is symmetrical with respect to the $Ω$ axes. The curves corresponding to negative frequencies are not shown in this figure for the sake of simplicity.

Fig. 3The Campbell diagram of a composite rotor for various ply angles

From Fig. 3, it is seen that in the inertial frame, natural frequencies of the rotating composite rotor become function of rotating speed as flapping bending and sweeping bending motions are coupled due to the presence of internal damping (there is no the Corilois effect in the inertial frame). In fact, natural frequencies of a rotating shaft are constant at all rotating speeds in the inertial frame if the influence of internal damping is not considered. The reason is that for the rotating shaft of circular section, the displacement and acceleration vectors and the elastic restoring force always have same direction [18]. As seen in Fig. 3, the effect of the ply angle and rotating speed appear to be more significant for the higher-mode ones.

Fig. 4 shows the variation of the first two flexural dampings vs. rotating speed for various ply angles. It can be seen clearly that as the rotating speed is increased, the dampings of the backward modes decrease and remain negative for all rotating speed, so the backward modes are stable. From the results of Fig. 4, it can be also observed that the dampings corresponding to the forward modes are negative at low rotating speed and increase with increasing rotating speed, and at certain value of rotating speed the dampings vanish and then become positive. Transformation of dampings from negative to positive values marks the onset of unstable motion. The rotating speed corresponding to zero damping is the threshold of instability of the composite rotor. From Fig. 4 it seems that the effect of the ply angle on the damping of the mode 2 is more significant.

Fig. 4The decay plot of a composite rotor for various ply angles

a) Mode 1

b) Mode 2

Fig. 5Root locus of a composite rotor with internal damping

Fig. 5 shows the root locus plot of the complex eigenvalue as a function of rotating speed for various ply angles. For brevity, only the first mode is represented. In the figure the rotating speed is varied from $Ω =$ 0.0 to $Ω =$ 200 rpm in increments of 8.0 rpm. It is readily seen that the curve crosses the imaginary axis into the right hand side of the locus (the unstable region) with the increase of the rotating speed. It is also important to note that in the stable region, the decrease of ply angle will exert a stabilizing effect on the composite rotor, as seen by the curves with lower ply angles being further into the negative region of the locus. In contrast, in the unstable region, the decrease of ply angle promotes further instability, shown by the curves with lower ply angles being further into the right hand side of the locus.

Fig. 6 shows the effect of ply angle on the first critical rotating speed. It can be seen that as the ply angle increases, the critical rotating speed decreases and the critical speed is maximum at $θ =$ 0°.

Fig. 7 shows the effect of ply angle on the threshold of instability for the flexural mode. It is evident that the general effect of the ply angle on the threshold of instability is similar to that associated with the critical rotating speeds. By comparing Fig. 6 with Fig. 7, it may be noted that for the same ply angle the threshold of instability is larger than the critical rotating speed. This implies that the onset of instability always occurs after the critical rotating speed.

Fig. 6The variation of the first critical speed with ply angle

Fig. 7The variation of thresholds of instability with ply angle

Fig. 8The time history of deflections and the trajectories of the geometric center of the rigid disk (Ω= 20 rpm)

a)

b)

c)

The transient responses for the geometric center of the rigid disk due to the initial perturbations are calculated by integrating Eq. (20) (with $\{f\} = 0$ ) using the fourth order Runge-Kutta method. Figs. 8-10 present the time history of deflections and the trajectories of the geometric center of the rigid disk for rotating speeds 20 rpm, 27.89 rpm and 32.5 rpm, respectively. In the simulation, the ply angle $θ =$ 45°, and the initial conditions $X_{j v} (0) =$ 10^-5, $X_{j w} (0) =$ 10^-5, ${\dot{X}}_{j v} (0) = 0$ and ${\dot{X}}_{j w} (0) = 0$ are specified. Fig. 8 shows that at rotating speed below threshold of instability the oscillation in the two direction can be damped out by internal damping. Fig. 9 shows that when rotating speed is equal to the threshold of instability a whirling response with constant amplitude occurs as internal damping equals to zero in this case. It can also be noted from Fig. 10 that the forward mode of a supercritical composite rotor is unstable since the internal damping becomes positive in this case.

Fig. 9The time history of deflections and the trajectories of the geometric center of the rigid disk (Ω= 27.89 rpm)

a)

b)

c)

Finally, the steady state response of the composite rotor under the action of unbalance rigid disk is studied. The effect of internal damping on the steady state response is illustrated. In the simulation, the mode superposition method is used to obtain frequency response plots. Fig. 11 shows the two frequency response plots, the plot in dashed line corresponding to the case in which the internal damping are considered and the solid line corresponding to the case in which they are not. As can be seen from Fig. 11, internal damping tends to reduce resonance peak amplitude at critical speeds. This shows that the effect of internal damping on the the forced response is the same as external damping. Fig. 12 shows the effect of ply angle on unbalance steady state response with internal damping. It can be observed that the resonance peak amplitude of the lower mode is much stronger affected by internal damping, on the other hand, the effect of internal damping on the resonance frequency of the higher mode is more significant.

Fig. 10The time history of deflections and the trajectories of the geometric center of the rigid disk (Ω= 32.5 rpm)

a)

b)

c)

Fig. 11The frequency response of composite rotor with and without internal damping

6. Conclusions

A model was presented for the study of the dynamical behavior of the composite rotor with internal damping. The presented model was used to predict the natural frequencies, critical rotating speeds, instability thresholds, the transient free vibration responses and the steady state unbalance responses. Theoretical solutions of the free vibration of the system were determined by applying Galerkin’s method. From the present analysis and the numerical results, the following main conclusions were drawn:

1) The developed model provides means of predicting of the dynamic behavior of rotating composite rotors with internal damping.

Fig. 12The frequency response of composite rotor with various internal damping

2) The ply angle affect the steady state unbalance responses and instability behavior of composite rotors significantly.

3) There is an obvious decrease in the first critical rotating speed and instability threshold as the ply angle is increased. And the critical rotating speed and threshold of instability have their maximum values at $θ =$ 0°.

4) The onset of instability always occurs after the critical rotating speed.

References

Chandra R., Singh S. P., Gupta K. Damping studies in fiber-reinforced composites-a review. Composite Structures, Vol. 46, Issue 1, 1999, p. 41-51.

Search CrossRef
Bucciarelli L. L. On the instability of rotating shafts due to internal damping. Journal of Applied Mechanics, Vol. 49, 1982, p. 425-434.

Search CrossRef
Kimball Jr A. L. Internal friction theory of shaft whirling. General Electric Review, Vol. 27, 1924, p. 244-251.

Search CrossRef
Gunter E. J. Rotor-bearing stability. Proceedings of the First Turbo-machinery Symposium, 1972, p. 119-141.

Search CrossRef
Vance J. M., Lee J. Stability of high speed rotors with internal friction. ASME Journal of Engineering for Industry, Vol. 96, 1974, p. 960-968.

Search CrossRef
Melanson J., Zu J. W. Free vbration and stability analysis of internally damped rotating shafts with general boundary conditions. ASME Journal of Vibration and Acoustics, Vol. 120, 1998, p. 776-783.

Search CrossRef
Montagnier O., Hochard Ch. Dynamic instability of supercritical driveshafts mounted on dissipative supports-effects of viscous and hysteretic internal damping. Journal of Sound and Vibration, Vol. 305, 2007, p. 378-400.

Search CrossRef
Simon A., Flowers G. T. Adaptive disturbance rejection and stabilization for rotor systems with internal damping. Proceedings of the ASME International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, IDETC/CIE, 2009, San Diego, California, USA.

Search CrossRef
Singh S. P., Gupta K. Composite shaft rotordynamic analysis using a layerwise theory. Journal of Sound and Vibration, Vol. 91, Issue 5, 1996, p. 739-756.

Search CrossRef
Kim W., Argento A., Scott R. A. Forced vibration and dynamic stability of a rotating tapered com posite Timoshenko shaft: bending motions in end-milling operations. Journal of Sound and Vibration, Vol. 246, Issue 4, 2001, p. 563-600.

Search CrossRef
Mazzei A. J., Scott R. A. Effects of internal viscous damping on the stability of a rotating shaft driven through a universal joint. Journal of Sound and Vibration, Vol. 265, Issue 4, 2003, p. 863-885.

Search CrossRef
Sino R., Baranger T. N., Chatelet E., Jacque T. G. Dynamic analysis of a rotating composite shaft. Composites Science and Technology, Vol. 68, 2008, p. 337-345.

Search CrossRef
Saravanos D. A., Varelis D., Plagianakos T. S., et al. A shear beam finite element for the damping analysis of tubular laminated composite beams. Journal of Sound and Vibration, Vol. 291, 2006, p. 802-823.

Search CrossRef
Berdichevsky V. L., Armanios E., Badir A. M. Theory of anisotropic thin-walled closed-cross-section beams. Composites Engineering, Vol. 2, Issue 5-7, 1992, p. 411-432.

Search CrossRef
Ren Y. S., et al. Structure damping of thin-walled composite one-cell beams. Journal of Vibration and Shock, Vol. 31, Issue 3, 2012, p. 141-146,152, (in Chinese).

Search CrossRef
Smith E. C., Chopra I. Formulation and evaluation of an analytical model for composite box-beams. Journal of American Helicopter Society, Vol. 36, Issue 3, 1991, p. 23-25.

Search CrossRef
Banerjee J. R., Su H. Development of a dynamic stiffness matrix for free vibration analysis of spinning beams. Computers and Structures, Vol. 82, 2004, p. 2189-2197.

Search CrossRef
Kim J. Rotation effects on vibration of structures seen from a rotating beam simply supported off the rotation axis. Journal of Vibration and Acoustics, Vol. 128, 2006, p. 328-337.

Search CrossRef

About this article

Received

03 September 2014

Accepted

19 November 2014

Published

30 December 2014

Keywords

composite shaft

internal damping

vibration and stability

Acknowledgements

The research is funded by the National Natural Science Foundation of China (Grant Nos. 11272190), Shandong Provincial Natural Science Foundation of China (Grant Nos. ZR2011EEM031) and Graduate Innovation Project of Shandong University of Science and Technology of China (Grant Nos. YC130210).

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.