The influence of cracks in rotating shafts

Jean-Jacques Sinou, A.W. Lees

To cite this version:

Jean-Jacques Sinou, A.W. Lees. The influence of cracks in rotating shafts. Journal of Sound and
Vibration, Elsevier, 2005, 285 (4-5), pp.1015-1037. <10.1016/j.jsv.2004.09.008>. <hal-00207691>

HAL Id: hal-00207691

https://hal.archives-ouvertes.fr/hal-00207691
Submitted on 18 Jan 2008

HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est

archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents
entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non,
lished or not. The documents may come from émanant des établissements d’enseignement et de
teaching and research institutions in France or recherche français ou étrangers, des laboratoires
abroad, or from public or private research centers. publics ou privés.
Journal home page: http://www.sciencedirect.com/science/journal/0022460X

The influence of cracks in rotating shafts

Journal of Sound and Vibration, Volume 285, Issues 4-5, 6 August 2005, Pages 1015-1037
J-J. Sinou and A.W. Lees

THE INFLUENCE OF CRACKS IN ROTATING SHAFTS

J-J. Sinou1* and A.W. Lees2

Department of Mechanical Engineering, University of Wales Swansea,

Singleton Park, Swansea SA2 8PP, Wales, UK.

2
E-mail: a.w.lees@swansea.ac.uk
1
E-mail: jean-jacques.sinou@ec-lyon.fr
Permanent address: Ecole Centrale de Lyon, Laboratoire de Tribologie et Dynamique des Systèmes
UMR CNRS 5513, 36 avenue Guy de Collongue, 69134 Ecully, France.

ABSTRACT
In this paper, the influence of transverse cracks in a rotating shaft is analyzed. The paper addresses
the two distinct issues of the changes in modal properties and the influence of crack breathing on
dynamic response during operation. Moreover, the evolution of the orbit of a cracked rotor near half
of the first resonance frequency is investigated. The results provide a possible basis for an on-line
monitoring system.
In order to conduct this study, the dynamic response of a rotor with a breathing crack is evaluated by
using the Alternate Frequency/Time Domain approach. It is shown that this method evaluates the
non-linear behaviour of the rotor system rapidly and efficiently by modeling the breathing crack with
a truncated Fourier series. The dynamic response obtained by applying this method is compared with
that evaluated through numerical integration. The resulting orbit during transient operation is
presented and some distinguishing features of a cracked rotor are examined.

1. INTRODUCTION

The influence of a transverse crack in the shaft of rotating machines on the associated dynamic
behaviour, has been a focus of attention for many researchers [1-35]. The presence of a crack may
lead to a dangerous and catastrophic effect on the dynamic behaviour of rotating structures and cause
serious damage to rotating machinery. Therefore, the timely detection of a rotor crack would
potentially avoid severe damage and expensive repairs due to the failure of rotating machinery as
well as assuring the safety of personnel.
Generally, two different approaches are attempted to identify the presence of a crack in rotating
structures. The first approach is based on the fact that the presence of a crack in rotating shaft reduces
the stiffness of the structure, hence reducing the natural frequencies of the original uncracked shaft.
1
Various theoretical and experimental works [3-10] performed over the last three decades have
indicated that the change in modal properties (natural frequencies and modes shapes) may be useful
for the detection of a crack, as well as for the identification of both crack depth and location.
Moreover, the influence of opening and closing of crack due to the shaft self-weight for various
orientations of the shaft has been investigated [3,11] and showed the effectiveness of the change in
natural frequencies versus the orientation of the shaft to detect the orientation of the crack front.
Another approach of crack identification is based on the modification of the dynamic responses of
the crack rotor during its rotation. Indeed dynamic analysis of the cracked rotor based on theoretical
and experimental studies has been a subject of great interest for the last decades [12-22]. Wauer [1]
reviewed a literature survey on the state-of-art of the dynamics of cracked rotors. Mayes and Davies
[11] analyzed the response of a multi-rotor-bearing system containing a transverse crack in a rotor
both experimentally and theoretically. Gasch [12,13] studied the dynamic behaviour of a simple rotor
with a cross-sectional crack and the associated stability behaviour due to the crack and imbalance.
Henry and Okah-Avae [23] investigated the effects of the gravity and unbalance on the dynamic
behavior of a crack shaft. Moreover, many studies [3, 20, 22, 24-28] indicated the change in dynamic
responses and more particularly the tendency of the rotor to exhibit a harmonic component at twice
shaft speed close to half any resonance frequency. Indeed, due to the shaft self-weight and the
rotation of the rotor, the crack opens and closes during a complete revolution of the rotor; hence the
stiffness of the shaft varies. This opening and closing mechanism, which is called the breathing
effect, induces vibration of the second and higher harmonics of rotating speed in frequency domain.
Although the presence of this twice per revolution component can indicate the possibility of the
presence of a crack, it is also well known that the presence of such a component can be generated by
shaft misalignment, asymmetric shaft, looseness of bolts and nuts, or a range of other non-linearities
[26-27]. Therefore, this feature, in itself, is insufficient to indicate the presence of a transverse crack
in rotors. In this way, the additional observation of the orbits can be useful for revealing the presence
of a crack [24-29]. Effectively, near the speed range of half any resonant frequency, the orbit changes
from a single loop to double loops with speed, and then an internal loop may appear. This
observation is the signature of the presence of a crack that indicates the change in amplitude and
phase at half any resonance speed, and is also a characteristic for signals containing two vibration
components with the same direction of precession [24,27].
The aim of this paper is to investigate these effects by taking into account the non-linear dynamical
behaviour due to the breathing transverse crack in order to obtain some indications that might be
useful in detecting the presence of a crack in rotating system. Some parametric studies regarding the
location and depth of cracks are carried out in order to show their influence on the change in
frequencies. Next the influence of a transverse breathing crack on the response of a rotor model, is
investigated using the Alternate Frequency/Time Domain method [36-38] with a path following
procedure in order to predict the non-linear dynamics and also the twice per revolution component
close to half resonance response associated with the presence of a breathing crack. It is shown that
the non-linear behavior of the rotor system with a breathing crack may be obtained by modeling the
crack with a truncated Fourier series. In order to validate this approach, the results are compared
with those obtained through numerical integration. Moreover, it is demonstrated that a rotor with a
breathing crack, which opens and closes during rotations, shows a non-linear dynamic behavior due
to the variation of the rotor's stiffness during its rotation. Finally, the non-linear behavior of the crack
rotor at half the first shaft resonance speed is calculated, and the influence of the crack depth is
presented to provide information for an on-line identification of rotor cracks.

2
 2. MATHEMATICAL MODEL

This section describes the derivation of the system model which is used in the subsequent analysis.

2.1. ROTOR SYSTEM

The shaft is discretized into a number of Timoshenko beam finite elements having with four degrees
of freedoms at each node [39]:
(MTe + M eR )x ie + (CeR − ΩG eR )x ie + K eB x ie = Qie (1)
where M Te and M eR are the rotary and translational mass matrix of the shaft element. C eR , G eR , and
K eB are the external damping, gyroscopic, and stiffness matrices respectively. Q ie defines the
gravity force vector for the shaft. Ω is the rotational speed and the factor of damping for the shaft.
The damping is taken as classical for the sake of simplicity and CeR = βK eB where the β is a constant
factor of proportionality and internal rotor damping has been neglected.
The modeling of the rigid discs is given by
( )
M Td + M dR x id − ΩG d x id = Q id + Wid (2)
where M Td , M dR , and G dR are the translational mass, rotary mass and gyroscopic matrices
respectively. Q id consists of the weight of the disc.
Wid defines the unbalance forces due to disc
having mass m with an eccentricity e .
Finally, the discrete bearing stiffness coefficients are placed at the corresponding degrees of freedom
and the equation of motion for the complete rotor system is defined as follows
Mx + (C + ΩG )x + Kx = Q + W (3)
where M and G are the mass and gyroscopic matrices including mass and gyroscopic matrices of
the shaft and rigid discs. C and K are the external damping and stiffness matrices of the shaft. Q
and W define the vector of gravity force and imbalance force for the complete rotor system.

2.2. THE CRACKED ROTOR

In this section, the modeling of the crack and the breathing mechanism are discussed briefly. Mayes
and Davies [3] demonstrated that a transverse crack in a rotor shaft can be represented by the
reduction of the second moment of area ∆I of the element at the location of the crack. By using
Rayleigh's method, they obtained that the change in ∆I verified [3,5]
∆I I 0
= (1 − ν 2 ) F ( µ )
R
(4)
1 − ∆I I 0 l
where I 0 , R , l , ν , µ , and F ( µ ) are the second moment of area, the shaft radius, the length of the
section, the Poisson's ratio, the non-dimensional crack depth, and the compliance functions varied
with the non-dimensional crack depth µ , respectively. The non-dimensional crack depth µ is given
by µ = a R where a defines the crack depth of the shaft as illustrated in the Figure 1. Although
F ( µ ) can be derived from the appropriate stress factor, Mayes and Davies indicated that a good
approximation of this function can be achieved by considering the fractional change in the second
moment of area measured at the crack face. Their observations were validated by experimental tests
consisting of measuring the two functions FX ( µ ) and FY ( µ ) [3] for both the X and Y directions
indicated the direction of the crack front and the orthogonal direction associated as illustrated in
Figure 1. The values of the second moment of area with the new centroid associated are given in

3
Appendix A. Then the stiffness matrix due to the transversal crack K crack can be obtained at the crack
location, by using standard finite elements [5].
When the rotor is cracked, the opening and closing behavior due to the rotor rotation and shaft self-
weight results in a time dependent stiffness. To accurately predict the dynamic response of the rotor
system with a breathing crack, an appropriate crack model is essential. Many researchers have
studied this problem [11-13, 17] and various crack opening and closing models have been developed.
Penny and Friswell [16] compared different crack models: the hinge model (13) where the breathing
switches from it open to closed state abruptly during the rotation of the shaft; the crack model of
Mayes and Davies [11] where the opening and closing of the crack was described by a cosine
function; and finally the crack model of Jun et all [22] based on fracture mechanics where the
coupling stiffness as well as the direct stiffness are calculated as the shaft rotates (and so the crack
opens and closes). Penny and Friswell [16] demonstrated that for monitoring using low frequency
vibration, simple models of crack are adequate and sufficient for the prediction of the dynamic
behavior of a rotor with breathing crack, as well as the predicted whirl orbit at the steady-state
2 × harmonic of rotor speed. The extent of crack opening will be determined by the proportion of the
crack face which is subject to tensile axial stresses. By assuming that the gravity force is much
greater that the imbalance force, the function describing the breathing crack [1,3-4,11-13,16-17] may
be chosen as
1
f (t ) = (1 − cos(Ωt )) (5)
2
where Ω is the rotational speed of the rotor. As illustrated in the Figure 1, for f (t ) = 0 , the crack is
totally closed and the cracked rotor stiffness is equal to the uncracked rotor stiffness. For f (t ) = 1 ,
the crack is full open.
Finally, the dynamics equation of the rotor with a breathing crack, in this linearised approximation,
can be defined as
Mx + ( C + ΩG ) x + ( K − f ( t ).K crack ) x = Q + W (6)
where x , x , and x are the acceleration, velocity and displacement of the degree-of freedom of the
cracked rotor system. M , C , G and K define the mass, damping, gyroscopic and stiffness
matrices, respectively. Q and W are the vector of gravity force and imbalance force,
respectively. K crack is the stiffness matrix due to the crack and f (t ) the function representing the
breathing effect.
In this study, a rotor shaft of 1m length and 10mm diameter; two discs of 40mm diameter and 15mm
thickness are situated at 0.3m of each end of the shaft, as illustrated in Figure 2. A crack is added at
one third of the left end and a mass m at the eccentricity e is placed on the first disc. All the values of
the parameters are given in Table 1.

3. NATURAL FEQUENCIES

In this section, the changes in the natural frequencies of the rotor, which is the common first step in
the diagnosis of a crack, have been examined. The crack position, depth and the orientation of the
crack have been varied.
Tables 2 gives the values of the natural frequency for the uncracked and cracked shaft with the
variation of the non-dimensional crack depth µ . Table 3 gives the values of the natural frequency for
the uncracked and cracked shaft with the variation of the crack location (with µ = 1 ).
The natural frequencies associated with the vertical (first, third and fifth modes) and horizontal
(second, fourth, and sixth modes) frequencies are equal in the case of an uncracked shaft but are

4
different for the cracked rotor due to the presence of the crack. The highest changes in natural
frequencies occur in the vertical mode due to the orientation of the crack and the shaft self-height.
In order to compare the frequencies of the cracked and uncracked shaft, the percentage change in
natural frequencies is defined as follow
fi uncracked − fi cracked
%Ci = ×100 (7)
fi uncracked
where i defines the i th -frequency of the system.
Figure 3 shows the percentages changes in the first and second natural frequencies with the rotation
of the shaft for various non-dimensional crack depths µ . In this case the area of the open crack due
to the shaft's self-weight is calculated at each orientation of the shaft [5]. As observed previously by
Lees and Friswell [5], this simple procedure is sufficient to reflect the change in natural frequencies
due to the rotation of the shaft and the physics associated. If the non-dimensional crack depth µ is
equal 1(corresponding to the lost of half the shaft's area), it appears that the crack is fully open at a
orientation of 0 degree, and fully closed at 180 degrees, as illustrated in Figure 3(a). Moreover, one
observes that for a non-dimensional crack depth µ is less than 1, the maximum change in frequencies
is defined for ranges of shaft orientations that increases when the non-dimensional crack depth µ
decreases. This reflects the fact that the crack and the equivalent area of the cracked shaft due to the
shaft self-weight are completely open and closed for ranges of shaft orientations. Of course the
change in the natural frequencies increase when the crack depth increases.
The crack depth, the location of the crack, and the rotation of the shaft clearly effects the natural
frequencies of the rotor and these changes in modal properties can be used as an identification of
transverse crack in rotating machines. However this approach is time consuming and requires the
knowledge of the original natural frequencies of the rotor in its uncracked condition. This last point
may not be easy to obtain due to the fact that since the modal tests have been carried out, the natural
frequencies may have been modified during time rotating machines operation. Therefore the use of
the dynamic response in order to identify the presence of a transverse crack appears to be
unavoidable and the most convenient.

4. NON-LINEAR DYNAMICAL BEHAVIOUR

Due to the presence of a breathing crack and the mathematical model associated, the
determination of the dynamic behaviour of the rotor requires considerable computational resources
by using a classical numerical integration. In order to avoid these computational problems, one of the
most efficient and systematic approaches is the use of the harmonic balance method [38] that permit
the discretisation of the unknowns functions in time by using their Fourier components, which are
assumed to be constant with respect to time. Such methods include the incremental harmonic balance
method [41], the Fast Galerkin method [40,42], and the alternate frequency/time domain method [36-
38].
In this study, the Alternate Frequency/Time domain method (AFT method) is applied in order to
obtain the non-linear response of the rotor with the breathing crack.

4.1. ALTERNATE FREQUENCY/TIME DOMAIN METHOD

The system defined in equation (6) can be written more generally as
Mx + ( C + ΩG ) x + ( K − f ( t , x ).K crack ) x = Q + W (8)
where there is no restriction to the unbalance magnitude.

5
In the linearised crack model, it is usual to represent the portion of the crack below the neutral axis as
being open. There are a number of approaches differing only in detail but one such approach is that
of Mayes and Davies [3]. The area of the crack face below the neutral axis is calculated and this is
then represented with a chordal equivalent crack having the same area. The important approximation
here is that the neutral axis is usually taken as the centroid of the shaft, but this not exact.
With any non-zero deflection, the centre line of the shaft will be stretched and hence there will be a
tensile stress and this being so, for a shaft deflected under gravity, the true neutral axis will be above
the centre line. For small vibration amplitudes, the error introduced by this approximation is
relatively minor, simply reducing the effective crack by a constant factor. However, for appreciable
vibration amplitudes the error becomes significant and time dependent. It is beyond the scope of the
present paper to analyze the sensitivity of this effect but the correction (non-linear) terms are
important when the vibration amplitude is appreciable relative to the catenary.
In the present paper the position of the neutral axis at the crack location has been re-calculated at
each time step in order to give an accurate model of the crack dynamics.
This non-linear equation may be re-written
Mx + D(Ω )x + Kx + f (x, Ω, t ) − g(x, Ω, t ) = 0 (9)
where M , D and K are the mass, the damping and the stiffness matrices. f and g are the vectors
containing the non-linear expressions due to the breathing crack and the vector for the imbalance and
gravity force, respectively. In the absence of response to unbalance, the function f becomes a linear
function of x ; for non-negligible imbalance, the addition stresses interact with those arising from the
(gravity induced) catenary, so changing the neutral axis, and hence the active cross section of the
crack.
Setting x = x k + ∆x , x = x k + ∆x and x = x k + ∆x (where k defined the k th iteration process), one
considers the truncated Fourier series expansion
m
x = X0 + ∑ ⎡⎣ X 2i −1cos ( iΩt ) + X 2i sin ( iΩt ) ⎤⎦ (10)
i =1

where X 0 , X 2i −1 and X 2i are the Fourier coefficients of x , and the associated truncated Fourier
series expansion for ∆x
m
∆x = ∆X0 + ∑ ⎡⎣ ∆X 2i −1cos ( iΩt ) + ∆X 2i sin ( iΩt ) ⎤⎦ (11)
i =1

where ∆X 0 , ∆X 2i −1 and ∆X 2i are the Fourier coefficients of ∆x . The matrix of Fourier coefficients
of x and ∆ x are indicated and arranged as follows:
⎡ X 1,0 X i ,0 X n ,0 ⎤
⎢ ⎥
⎢ ⎥
X = ⎢ X 1,2 j X i,2 j X n,2 j ⎥ (12)
⎢ ⎥
⎢ ⎥
⎢ X 1,2 m X i ,2m X n , 2 m ⎥⎦
⎣
and

6
 ⎡ ∆X 1,0 ∆X i ,0 ∆X n , 0 ⎤
⎢ ⎥
⎢ ⎥
∆X = ⎢ ∆X 1,2 j ∆X i , 2 j ∆X n ,2 j ⎥ (13)
⎢ ⎥
⎢ ⎥
⎢∆X 1, 2 m ∆X i , 2 m ∆X n , 2 m ⎥⎦
⎣
where X i , j and ∆X i , j define the j th component of the Fourier coefficients for the i th degree-of-
freedom for the non-linear system. The number of harmonic coefficients m is selected in order to
consider only the significant harmonics expected in the solution. By replacing x and ∆x by their
Fourier series, one obtains (2m + 1) × n linear algebraic equations
AX + F − G + A.∆X = 0 (14)
with
⎡K ⎤
⎢ ⎥
⎢ ⎡−Ω2M + K ΩD ⎤ ⎥
⎢ ⎢ ⎥ ⎥
⎣ −Ω D −Ω 2
M + K ⎦
⎢ ⎥
⎢ ⎥
⎢ ⎥
A=⎢ ⎡− j Ω M + K
2 2
jΩD ⎤ ⎥
⎢ ⎢ ⎥ ⎥
⎣ − jΩD − j Ω
2 2
M + K ⎦
⎢ ⎥
⎢ ⎥
⎢ ⎥
⎢ ⎡−m Ω M + K
2 2
mΩD ⎤⎥
⎢⎢ ⎢ ⎥
⎣ ⎣ −mΩD −m 2 Ω2M + K⎦ ⎥⎥⎦
(15)
F and G represent the Fourier coefficients of the function f and g , respectively. F is difficult to
determine from the Fourier coefficients directly due to the dependence in time for this function.
Hence F can be calculated by using the following path (Alternating Frequency/Time domain AFT)
as follows
Γ −1 Γ
X ⎯⎯→
⎯ x(t ) ⎯
⎯→ f (t ) ⎯⎯→ F (16)
where Γ and Γ −1 define the Discrete Fourier Transform from time to frequency domain, and from
frequency to time domain, respectively. The DFT from time to frequency domain is given by
⎧ 1
⎪ for i = 1
⎪ 2m + 1
⎪ ⎛ ( j − 1) iπ ⎞
⎪ 2
Γij = ⎨ cos ⎜ ⎟ for i = 2, 4,..., 2m for j = 1, 2,..., 2m + 1 (17)
⎪ 2 m + 1 ⎝ 2 m + 1 ⎠
⎪
⎪ 2 sin ⎜⎛ ( j − 1)( i − 1) π ⎟⎞ for i = 3,..., 2m + 1
⎪⎩ 2m + 1 ⎝ 2m + 1 ⎠
and from frequency time domain

7
 ⎧
⎪
⎪1 for j =1
⎪ ⎛ i − 1 jπ ⎞
⎪
Γij −1 = ⎨cos ⎜
( ) for j = 2, 4,..., 2m for i = 1, 2,..., 2m + 1 . (18)
⎟
⎪ ⎝ 2m + 1 ⎠
⎪ ⎛ i −1 j −1 π ⎞
⎪sin ⎜ ( )( )
⎟ for j = 3,..., 2m + 1
⎪⎩ ⎝ 2m + 1 ⎠
Finally, the error vector R is given by
R = AX + F − G (19)
and the associated convergence values are defined by

δ 1 = R 02 + ∑ ( R 22 j −1 + R 22 j )
m
(20)
j =1

and

δ 2 = ∆X02 + ∑ ( ∆X 22 j −1 + ∆X 22 j )
m
(21)
j =1

The global procedure of the Alternate/Frequency Time domain method is illustrated in Figure 4.

4.2. PATH CONTINUATION

Usually, the dynamics of system and the solution associated have to be calculated at different
parameter values consecutively (in this study, the considered parameter is the speed of shaft
rotation). In order to reduce the time required for the calculation the path following technique can be
used (Narayanan and Sekar [38]). One considers the estimation of the neighboring point on the
solution branch by using the Lagrangian polynomial extrapolation method with four points. So, one
assumes that four points on the solution branch are obtained a priori in order to begin the
extrapolation scheme. The varied parameter is the speed rotation of the shaft Ω . Any point on the
solution branch is represented at X i , where X i is the Fourier coefficients of x . The arc length
between two consecutive points X i +1 and X i is given by δli +1 = (X i+1 − X i )T (X i+1 − X i ) for
i = 0,…,2 . Next, the arc length parameters are calculated as follows
⎧ L0 = 0
⎪
⎪ L1 = δl1
⎪
⎨ L2 = L1 + δl 2 (22)
⎪ L = L + δl
⎪ 3 2 3
⎪⎩ L4 = L3 + ∆l
Finally, by using the Lagrangian extrapolation scheme, the following estimated point at the distance
∆l might be defined by
⎛ 3 ⎞
3
⎜ ⎛ L3 − L j ⎞ ⎟
X4 = ∑ ∏ ⎜ ⎟⎟ Xi ⎟ for i = 0,… ,3 . (23)
⎜ ⎜
i =1 ⎜ j = 0 ⎝ Li − L j ⎠ ⎟
⎝ i≠ j ⎠

8
 4.3. RESULTS AND DISCUSSIONS
The dynamic response of the rotor with the breathing crack is obtained by using the AFT method and
the path following procedure. The Fourier series is truncated in order to consider only the significant
harmonics. All the values of the parameters (damping, non-dimensional crack depth µ , mass
unbalance, eccentricity of the mass unbalance, etc…) are given in Table 1. It may be noted that the
crack depth, the level of damping, and the level of unbalance are important parameters because these
three all influence the magnitude of the non-linear term
In order to demonstrate the capability of the Alternate/Frequency Time domain method for
determining a good approximation of the dynamic response of the rotor with a breathing crack, a
comparison is performed with results obtained by using direct numerical integration such as the
Runge-Kutta 4. Figure 5 shows the comparison between times histories of the vertical and horizontal
displacements of the cracked rotor around the crack position ( L = 0.375m ), evaluated through
numerical integration and the AFT method for various numbers of harmonics m . The associated
orbits are plotted in Figure 6. These results highlight that the periodic solution ( m = 1 ) does not
provide a good approximation of the solution but that the second harmonic coefficients ( m = 2 ) are
sufficient to obtain an accurate evaluation of the dynamic response of the rotor with a breathing
crack. Therefore, these results clearly demonstrate that the presence of a breathing crack results in a
non-linear dynamical behavior.
Figure 7 shows the horizontal and vertical steady-state responses of the cracked and uncracked
shafts, at the position of the shaft L = 0.375m . Moreover Figure 8 shows the horizontal and vertical
steady-state responses of the cracked shaft for each node of the shaft. Figure 7 indicates a decrease in
the first and second vertical natural frequencies as indicated at the marks 3 and 5. This change is due
to the reduction in system stiffness resulting from the presence of the transverse crack. Although the
first and second horizontal natural frequencies decrease due to the presence of the crack, it may be
observed that the differences between the cracked and uncracked frequencies are very small and very
difficult to detect in practice, as indicated at the marks 7 and 8 in the Figure 7.
In the case of the cracked shaft, 2X harmonics are observed in the horizontal and vertical directions
at one-half of the first vertical and horizontal frequencies, as illustrated at the marks 2 and 6 in
Figures 7 and 8. 3X harmonics in the vertical direction for the first frequency are also present at one-
third of the first vertical natural frequency as shown in Figures 7 and 8 at the mark 1. Moreover, 2X
harmonics in the vertical direction for the second vertical frequency at one-half of the second vertical
frequency. It should be noted that there are no visible 3X harmonics of the first horizontal natural
frequency and 2X harmonic of the second horizontal natural frequency. Considering the rotor without
a crack, no harmonics are predicted for the dynamic behavior of the rotor. In this case, the dynamic
response is purely synchronous.
In order to compare the dynamic behavior of the cracked and uncracked shaft near one-half and one-
third of the natural frequencies, the differences between the cracked and uncracked horizontal and
vertical amplitudes are calculated as follows
∆X = Max (X cracked − X uncracked ) (24)
and
∆Y = Max (Ycracked − Yuncracked ) (25)
at each element of the shaft. X cracked , Ycracked , X uncracked , and Yuncracked define the horizontal and
vertical displacement for the cracked and uncracked shaft respectively. Therefore the Figure 9 shows
clearly the differences ∆X and ∆Y between the cracked and uncracked rotor. As illustrated in
Figure 9(A) and 9(B), the differences ∆X and ∆Y increase at the one-half and one-third of the first
horizontal and vertical frequencies. Moreover, the differences ∆X and ∆Y at each element of the

9
shaft indicate clearly the first mode shape with the maximum difference obtained at the middle of the
shaft. Figure 9(C) shows the differences ∆X and ∆Y around one-half of the second horizontal and
vertical frequencies. In this case the differences ∆X and ∆Y also increase one-half of the second
horizontal and vertical frequencies and indicate the deformations of the second mode shape with a
node in the middle of the shaft and the maximums of ∆X and ∆Y are obtained at one fourth on each
end of the shaft.
As explained previously, the diagnosis of the presence of a crack in rotating machinery based only on
the appearance of a harmonics response at half the natural frequency in the spectrum may be
misleading. So one of the features to detect the presence of a transversal crack in a rotating shaft is
the use of the evolution of the orbits during time around one-half of the resonance frequencies.
Now, the attention is focussed on the evolution of the orbits around one-half of the natural
frequencies. Figure 10 shows the evolution of the rotor's orbits. It is shown that a distortion in the
orbit appears and increases when the speed of the rotor increases. Next, the shape of the orbit
changes from a form with only one loop to a double loop. And finally, when the speed of the cracked
rotor is passing through nearly half of the critical speed, the orbit changes from a double loop to a
inside loop (a loop containing another small loop inside) that indicates the change in both the phase
and the amplitude of the harmonic components. This dynamic behavior of the rotor through the
passage of half of the rotor speeds is the signature of a crack in the rotor. Moreover, Figure 11 shows
the orbit obtained close to a third of the first critical speed, half of the first critical speed, and half of
the second critical speed. It may be observed that the change in the dynamic behavior is also detected
near half of the second critical speed with the formation of double loops at one end of the shaft and a
distortion form on the other end (Figure 11(c)). Showing the dynamic behavior at third of the first
critical speed (Figure 11(a)), the change of the orbit appears clearly with a triple loops that indicates
the presence of the harmonics. All these changes of the dynamic behavior of the rotor near half and
third of critical speeds may be considered as a feature for detection of crack in the rotor shafts.
Finally a parametric study is presented with various crack depth are performed in order to show the
influence of the crack depth by considering the variation and change of orbit around one-half of the
first frequency. As illustrated in Figure 12, the orbits decrease when the non-dimensional crack depth
µ decreases. Moreover, the frequency at which 2X harmonic occurs increases as the crack depth
decreases, due to the reduced stiffness resulting from the crack, which changes the natural
frequencies. This is illustrated by plotting (in Figure 13) the differences ∆X and ∆Y at the middle
position of the shaft for the various crack depth: as it can be seen the position of the maximum
differences versus the frequency increase with the decreasing of the crack depth, as well as the
differences ∆X and ∆Y decrease. Therefore, the behavior of the cracked rotor around one-half of
the natural frequencies and the observation of the associated orbits could provide crack diagnostic
information and could permit the estimation of the crack depth.

5. CONCLUSION

In this study, the influence of transversal cracks has been investigated: the change of the shaft
frequencies, as well as the harmonic component of the dynamical system response and the evolution
of the orbits are the principal effects due to the presence of a crack in a rotating shaft. More
particularly, the changes in the non-linear dynamical behavior of the rotor system through half
resonance speeds appear to be the classical signature for detecting the presence of a breathing crack.
Indeed, the distortion of the orbit, and formation of a double loop and inside loop in the orbit could
be considered as one of the most practical indicators of the presence of a transversal crack for health
monitoring purposes. Moreover the observation of the orbit amplitudes at half resonance speeds
could also provide crack diagnostic information about the crack depth: when the crack depth
10
increases, an increase of the orbit amplitudes is observed, as well as a decrease of the shaft speed at
which the 2X harmonic component of the dynamic response is maximum
In this paper, the use of the Alternate Frequency/Time Domain method and a path following
procedure allows obtaining rapidly and efficiently the non-linear dynamical behavior of a rotating
shaft with a transverse breathing crack. This method has advantages in terms of computating time. It
is easily implemented.

APPENDIX A: MOMENTS ABOUT THE CENTROID FOR THE CRACKED ROTOR

As illustrated in Figure 1, the cross section of the shaft at the location of the crack has asymmetric
~ ~
area moments inertia about the neutral axis of bending. The area moments of inertia I X and I Y
about the X and Y-axes are defined as
~
I X = ∫ Y 2 dA = ∫ Y 2 dXdY (26)
A A
~
I Y = ∫ X 2 dA = ∫ X 2 dXdY (27)
A A
where A defines the uncracked area of the cross section. After integrating over the uncracked area,
one obtains
~
IX =
R4 ⎡
4 ⎢
⎣
( )
(1 − µ ) 1 − 4µ + 2µ 2 2µ - µ 2 + cos −1 (1 − µ )⎤⎥
⎦
(28)

~
IY =
πR 4
4
⎡2
( )
+ R 4 ⎢ (1 − µ ) 2µ - µ 2
⎣3
32 1
( )( ) ( ) ⎤
+ (1 − µ ) 1 − 4µ + 2µ 2 2µ - µ 2 + sin −1 ⎛⎜ 2µ - µ 2 ⎞⎟⎥
4 ⎝ ⎠⎦
(29)

where R and µ are the shaft radius and the nondimensional crack depth ( µ = a R ). Then, the
moment of inertia about the centroidal axes I X and I Y are obtained
~
IX = IX (30)
~ 2
I Y = I Y − AX (31)
where X defines the centroid of the cross section. The uncracked area of the cross section A, and the
distance from the axis X to the centroid of the cross section X are given by
A = R 2 ⎡(1 − µ ) 2µ - µ 2 + cos −1 (1 − µ )⎤ (32)
⎢⎣ ⎥⎦

X =
2 3
3A
(
R 2µ - µ 2
32
. ) (33)

APPENDIX B: NOMENCLATURE

µ non-dimensional crack depth

Ω rotational speed
e eccentricity of the unbalance mass
Lcrack position of the crack
f i cracked ith natural frequency for the cracked system
f i uncracked ith natural frequency for the uncracked system
m number of Fourier coefficients retained
%Ci percentage change between the ith cracked and uncracked natural frequencies
x displacement of the degree-of-freedom

11
x velocity of the degree-of-freedom
x acceleration of the degree-of-freedom
M mass matrix of the rotor system
K stiffness matrix of the rotor system
K crack stiffness matrix for the cracked element
C damping matrix of the rotor system
G gyroscopic matrix of the rotor system
D global damping matrix of the rotor system
Q vector of gravity force for the rotor system
W vector of imbalance force for the rotor system
X matrix of the Fourier coefficients associated with the displacements of the rotor
system

REFERENCES

1. J. Wauer 1990 Applied Mechanics Review, 43, 13-17. Dynamics of Cracked Rotors:
Literature Survey.
2. A.D. Dimarogonas 1996 Engineering Fracture Mechanics, 55, 831-857. Vibration of
Cracked Structures: a State of the Art Review.
3. I.W. Mayes and W.G.R. Davies 1984 Transactions of the ASME Journal of Vibration,
Acoustics, Stress, and Reliability in Design, 106, 139-145. Analysis of the Response of a
Multi-Rotor-Bearing System Containing a Transverse Crack in a Rotor.
4. W. G. R. Davies and I. W. Mayes 1984 Transactions of the ASME Journal of Vibration,
Acoustics, Stress, and Reliability in Design, 106, 146-153. The Vibrational Behaviour of a
Multi-Shaft, Multi-Bearing System in the Presence of a Propagating Transverse Crack.
5. A.W. Lees and M.I. Friswell 2001 Proc. of 19th International Modal Analysis Conference.
The Vibration Signature of Chordal Cracks in Asymmetric Rotors.
6. S.W. Doebling, C.R. Farrar, and M.B. Prime 1998 The Shock and Vibration Digest, 30(2),
91-105. A Summary Review of Vibration-Based Damage Identification Methods.
7. L. Hamidi, J.B. Piaud, and M.Massoud 1992 IMechE Conference on Vibrations in Rotating
Machinery, 283-288. A Study of Cracks Influence on the Modal Characteristics of Rotors.
8. Y.S.Lee, and M.J. Chung 2000 Computers and Structures, 77, 327-342. A Study on Crack
Detection Using Eigenfrequency Test Data.
9. T.C. Tsai, and Y.Z. Wang 1996 Journal of Sound and Vibration, 192(3), 607-620. Vibration
Analysis and Diagnosis of a Cracked Shaft.
10. A.S. Sekhar 1999 Journal of Sound and Vibration, 223(4), 497-512. Vibration
Characteristics of a Cracked Rotor with two Open Cracks.
11. I.W. Mayes and W.G.R. Davies 1976 IMechE Conference on Vibrations in Rotating
Machinery, 53-64. The Vibrational Behaviour of a Rotating System Containing a Transverse
Crack.
12. R. Gasch 1976 IMechE Conference on Vibrations in Rotating Machinery, 123-128.
Dynamic Behaviour of a simple rotor with a cross-sectional Crack.
13. R. Gasch 1993 Journal of Sound and Vibration, 160(2), 313-332. A Survey of the Dynamic
Behaviour of a Simple Rotating Shaft with a Transverse Crack.
14. A.S. Sekhar, and B.S. Prabhu 1998 Mech. Mach. Theory, 33(8), 1167-1175. Condition
Monitoring of Cracked Rotors Through Transient Response.

12
15. N. Pugno, C. Surace, and R. Ruotolo 2000 Journal of Sound and Vibration, 235(5), 749-
762. Evaluation of the Non-linear Dynamic Response to Harmonic Excitation of a Beam
with Several Breathing Cracks.
16. J.E.T.P. Penny, and M.I. Friswell 2003 Key Engineering Materials, 245, 223-232.
Simplified Modelling of Rotor Cracks.
17. M.I. Friswell , and J.E.T.P. Penny 2002 International Journal of Structural Health
Monitoring, 1(2), 139-148. Crack Modelling for Structural Health Monitoring.
18. A.K. Darpe, K. Gupta, and A. Chawla 2003 Journal of Sound and Vibration, 259(3), 649-
675. Dynamics of a Two-crack Rotor.
19. S. Prabhakar, A.S. Sekhar, and A.R. Mohanty 2002 Journal of Sound and Vibration, 256(4),
773-786. Crack Versus Coupling Misalignment in a Transient Rotor System.
20. M-C. Wu and S-C. Huang 1998 International Journal of Mechanical Science, 40(6), 545-
555. Vibration and Crack Detection of a Rotor with Speed-Dependent Bearings.
21. S.C. Huang, Y.M. Huang, S.M. Shieh 1993 Journal of Sound and Vibration, 162(3), 387-
401. Vibration and Stability of a Rotating Shaft Containing a Transverse Crack.
22. O.S. Jun, H.J. Eun, Y.Y. Earmme, and C.W. Lee 1992 Journal of Sound and Vibration,
155(2), 273-290. Modelling and Vibration Analysis of a Simple Rotor with a Breathing
Crack.
23. T.A. Henry, and B.E. Okah-Avae 1976 IMechE Conference on Vibrations in Rotating
Machinery, 15-19. Vibrations in Cracked Shafts.
24. A. Muszynska 1989 Orbit, 10(2), Bently Nevada Corporation, USA, September.
Misalignment and Shaft Crack-Related Phase Relationships for 1 × and 2 × vibration
Components of Rotor Responses.
25. A. Muszynska, P. Goldman, and D.E. Bently 1992 IMechE Conference on Vibrations in
Rotating Machinery, 257-262. Torsional/Lateral Vibration Cross-Coupled Responses due to
Shaft Anisotropy: a New Tool in Shaft Crack Detection.
26. S.A. Adewusi, and B.O. Al-Bedoor 2002 Shock and Vibration, 9, 91-104. Experimental
Study on the Vibration of an Overhung Rotor with a Propagating Transverse Crack.
27. D.E. Bently, 1986 Orbit, 18-23, Bently Nevada Corporation, USA, January. Vibration
Analysis Techniques for Detecting and Diagnosing Shaft Cracks.
28. D.E. Bently, 1982 Orbit, 2, 10-11, Bently Nevada Corporation, USA, July. Detecting
Cracked Shafts at Earlier Levels.
29. D. Southwick., 1994 Orbit, 11-15, Bently Nevada Corporation, USA, June. Using Full
Spectrum Plots, Part-2.
30. M. Liao, and R. Gasch 1992 IMechE Conference on Vibrations in Rotating Machinery, 289-
295. Crack Detection in Rotating Shafts – An Experimental Study.
31. S. Prabhakar, A.S. Sekhar, and A.R. Mohanty 2001 Mechanical Systems and Signal
Processing, 15(2), 447-450. Detection and Monitoring of Cracks in a Rotor-Bearing System
Using Wavelet Transforms.
32. L.W. Chen, and H.-K. Chen 1995 Journal of Sound and Vibration, 188(4), 497-573.
Stability Analyses of a Cracked Shaft Subjected to the End Load.
33. A.S. Sekhar, and J.K. Dey 2000 Mechanism and Machine Theory, 35, 1657-1674. Effects of
Cracks on Rotor System Instability.
34. G.D. Gounaris, C.A. Papadoloulos 2002 Engineering Fracture Mechanics, 69, 339-352.
Crack Identification in Rotating Shafts by Coupled Response Measurements.
35. P.C. Paris, and G.C. Sih 1965 ASTM , STP No. 381, p.30. Stress Analysis of Cracks.
36. J.J. Sinou, F. Thouverez, and L. Jezequel 2003 Proc. of 8th International Conference on
Recent Advances in Structural Dynamics. Non-Linear Stability Analysis Based on the
13
 Center Manifold Approach, the Padé Approximants and the Alternate/Frequency Time
domain Method.
37. T.M.Cameron, and J.H. Griffin 1989 Journal of Applied Mechanics, 56, 149-154. An
Alternating Frequency/Time Domain Method for Calculating the Steady State Response of
Nonlinear Dynamic.
38. S. Narayanan, and P. Sekar 1998 Journal of Sound and Vibration, 211(3), 409-424. A
Frequency Domain Based Numeric-Analytical Method for Non-linear Dynamical Systems.
39. M. Lalanne, and G. Ferraris 1990 Rotordynamics - Prediction in Engineering, John Wiley
and Sons.
40. A.H. Nayfeh, and D.T. Mook 1979 Nonlinear Oscillations, New-York : John Wiley & Sons.
41. Y.K. Cheung, S.H. Chen, S.L. Lau 1990 Journal of Sound and Vibration, 140(2), 273-286.
Application of the Incremental Harmonic Balance Method to Cubic Non-linearity Systems.
42. F.H. Ling, and X.X. Wu 1987 International Journal of Non-Linear Mechanics, 22, 89-98.
Fast Galerkin Method and its Application to Determine Periodic Solutions of Non-linear
Oscillators.

14
 Parameters Physical dimension
Length of the shaft 1m
Diameter of the shaft 0.01 m
Young's modulus of elasticity E 2.1 × 1011 N.m-2
Shear modulus G 7.7 × 1010 N.m-2
Poisson ratio ν 0.3
Density ρ 7800 kg.m-3
Position of disc 1 0.3 m
5.1. POSITION OF DISC 2 0.7 m
Outer diameter of discs 1 and 2 0.04 m
Inner diameter of discs 1 and 2 0.01 m
Thickness of discs 1 and 2 0.015 m
Location of the crack Lcrack 0.375 m
Non-dimensional crack depth µ 1
Mass unbalance 0.005 g
Phase unbalance 0 degree
Eccentricity of the mass unbalance 0.02 m
Coefficient of damping β 10-5

Table 1: Details of the rotor model

15
 uncracked µ = 0.25 µ = 0.5 µ = 0.75 µ =1

1st frequency (Hz) 16.597 16.581 16.549 16.479 16.257

2nd frequency (Hz) 16.597 16.582 16.559 16.543 16.541
3rd frequency (Hz) 65.367 65.290 65.131 64.787 63.755
4th frequency (Hz) 65.367 65.292 65.178 65.102 65.088
5th frequency (Hz) 176.038 176.032 176.015 175.980 175.866
6th frequency (Hz) 176.038 176.032 176.020 176.012 176.011

Table 2: Evolution of the natural frequencies with the non-dimensional crack depth ( Lcrack = 0.375m )

uncracked Position of the crack Lcrack

0.025m 0.225m 0.475m 0.725m 0.975m
st
1 frequency (Hz) 16.597 16.593 16.399 16.169 16.332 16.593
2nd frequency (Hz) 16.597 16.597 16.564 16.525 16.553 16.597
3rd frequency (Hz) 65.367 65.316 63.673 65.321 63.779 65.323
4th frequency (Hz) 65.367 65.359 65.078 65.360 65.094 65.360
5th frequency (Hz) 176.038 175.756 173.312 171.737 174.523 175.612
6th frequency (Hz) 176.038 175.993 175.559 175.295 175.775 175.970

Table 3: Evolution of the natural frequencies with the location of the crack Lcrack ( non-dimensional
crack depth µ = 1 )

16
Figure 1 : Model of stiffness variation for the breathing crack

Figure 2 : Rotor model with two discs and a cracked element

17
Figure 3 : Changes in the first and second natural frequencies with the rotation of the shaft and the
non-dimensional crack depth µ
( µ = 0.25 , µ = 0.5 , µ = 0.75 , µ = 1)
(a) first frequency (b) second frequency

Figure 4 : Alternate Frequency/Time domain method

18
Figure 5 : Evolution of the horizontal and vertical displacements at the position of the crack
( Lcrack = 0.35m ) for f = 8.285Hz
( Runge-Kutta 4, m=1, m=2)

Figure 6 : Evolution of the orbit at the position of the crack ( Lcrack = 0.35m ) for f = 8.285Hz
( Runge-Kutta 4, m=1, m=2)

19
Figure 7 : Comparison between the cracked and uncracked shafts - horizontal and vertical
displacements at the position L = 0.35m ( uncracked, cracked)

Figure 8 : Vertical and horizontal displacements for the cracked shaft

20
Figure 9 : Differences ∆X and ∆Y between the cracked and uncracked shaft for (A) the 3X
resonance of the first frequency (B) the 2X resonance of the first frequency and (C) the 2X
resonance of the second frequency

21
Figure 10 : Evolution of the orbit at the position of the crack ( Lcrack = 0.35m ) for the 2X resonance of
the first frequency

22
Figure 11 : Evolution of the orbits at various position on the shaft around the static deflection
(a) 3X resonance of the first frequency (b) 2X resonance of the first frequency
(c) 2X resonance of the third frequency

23
Figure 12 : Evolution of the orbit at the position of the crack ( Lcrack = 0.35m ) around the 2X
resonance of the first frequency, with the variation of the non-dimensional crack depth µ
( µ = 0.25 , µ = 0.5 , µ = 0.75 , µ = 1)

24
Figure 13 : Evolution of the differences ∆X and ∆Y between the cracked and uncracked shaft for
the 2X resonance of the first frequency at the middle of the shaft with the variation of the non-
dimensional crack depth µ
( µ = 0.25 , µ = 0.5 , µ = 0.75 , µ = 1)

25