Mode shapes analysis of a cracked beam and its application for crack detection

In this paper, mode shapes of a cracked beam with a rectangular cross section beam are analysed using finite element method. The 3D beam element is applied for this finite element analysis. The influence of the coupling mechanism between horizontal bending and vertical bending vibrations due to the crack on the mode shapes is investigated. Due to the coupling mechanism the mode shapes of a beam change from plane curves to space curves. Thus, the existence of the crack can be detected based on the mode shapes: when the mode shapes are space curves there is a crack in the beam. Also, when there is a crack, the mode shapes have distortions or sharp changes at the crack position. Thus, the position of the crack can be determined as a position at which the mode shapes exhibit such distortions or sharp changes. While in previous studies using 2D beam element, distortions in the mode shapes caused by a small crack could not be detected, these distortions in the case using the 3D beam element can be amplified and inspected clearly by using the projections of the mode shapes on appropriate planes. The quantitative analysis is also implemented to relate the size and position of the crack with the observed coupled modes. These results can be applied for crack detection of a beam. In this paper, the stiffness matrix of a cracked element obtained from fracture mechanics is presented and numerical simulations of three case studies are provided.     

Experimental investigations of the response of a cracked rotor to periodic axial excitation

In this paper, the coupling of lateral and longitudinal vibrations due to the presence of transverse surface crack in a rotor is explored. A crack in a rotor is known to introduce coupling between lateral and longitudinal vibrations. Steady state unbalance response of a cracked rotor with a single centrally situated crack subjected to periodic axial impulses is investigated experimentally. The cracked rotor is excited axially using an electrodynamic exciter at a frequency equal to its bending natural frequency in both non-rotating and rotating conditions. The resulting time domain and frequency domain signals of the cracked rotor are studied. Spectral response of the cracked rotor with and without axial excitation is found to be distinctively different. When excited axially, it shows prominent presence of rotor bending natural frequency. However for an uncracked rotor, the response is similar with or without axial excitation. It is thus proposed that the response of the rotor to axial impulse excitation could be used for more reliable diagnosis of rotor cracks.     

ANALYSIS OF THE RESPONSE OF A CRACKED JEFFCOTT ROTOR TO AXIAL EXCITATION

The coupling of lateral and longitudinal vibrations due to the presence of transverse surface crack in a rotor is explored. Steady state unbalance response of a Jeffcott rotor with a single centrally situated crack subjected to periodic axial impulses is studied. Partial opening of crack is considered and the stress intensity factor at the crack tip is used to decide the extent of crack opening. A crack in a rotor is known to introduce coupling between lateral and longitudinal vibrations. Therefore, lateral vibration response of a cracked rotor to axial impulses is studied in detail. Spectral analysis of response to periodic multiple axial impulses shows the presence of rotor bending natural frequency as well as side bands around impulse excitation frequency and its harmonics due to modulations caused by rotor running frequency. It is concluded that the above approach can prove to be a useful tool in detecting cracks in rotors.     

FUNDAMENTAL FREQUENCY OF CRACKED BEAMS IN BENDING VIBRATIONS: AN ANALYTICAL APPROACH

This paper presents an analytical approach to the fundamental frequency of cracked Euler-Bernoulli beams in bending vibrations. The flexibility influence function method used to solve the problem leads to an eigenvalue problem formulated in integral form. The influence of the crack was represented by an elastic rotational spring connecting the two segments of the beam at the cracked section. In solving the problem, closed-form expressions for the approximated values of the fundamental frequency of cracked Euler-Bernoulli beams in bending vibrations are reached. The results obtained agree with those numerically obtained by the finite element method.     

The use of roving discs and orthogonal natural frequencies for crack identification and location in rotors

A variety of approaches that have been developed for the identification and localisation of cracks in a rotor system, which exploit natural frequencies, require a finite element model to obtain the natural frequencies of the intact rotor as baseline data. In fact, such approaches can give erroneous results about the location and depth of a crack if an inaccurate finite element model is used to represent an uncracked model. A new approach for the identification and localisation of cracks in rotor systems, which does not require the use of the natural frequencies of an intact rotor as a baseline data, is presented in this paper. The approach, named orthogonal natural frequencies (ONFs), is based only on the natural frequencies of the non-rotating cracked rotor in the two lateral bending vibration x–z and y–z planes. The approach uses the cracked natural frequencies in the horizontal x–z plane as the reference data instead of the intact natural frequencies. Also, a roving disc is traversed along the rotor in order to enhance the dynamics of the rotor at the cracked locations. At each spatial location of the roving disc, the two ONFs of the rotor–disc system are determined from which the corresponding ONF ratio is computed. The ONF ratios are normalised by the maximum ONF ratio to obtain normalised orthogonal natural frequency curves (NONFCs). The non-rotating cracked rotor is simulated by the finite element method using the Bernoulli–Euler beam theory. The unique characteristics of the proposed approach are the sharp, notched peaks at the crack locations but rounded peaks at non-cracked locations. These features facilitate the unambiguous identification and locations of cracks in rotors. The effects of crack depth, crack location, and mass of a roving disc are investigated. The results show that the proposed method has a great potential in the identification and localisation of cracks in a non-rotating cracked rotor.     

Coupled transverse and axial vibratory behaviour of cracked beam with end mass and rotary inertia

In this paper the vibrational behaviour of a cracked cantilever beam carrying end mass and rotary inertia is investigated. The transverse and axial vibrations of the beam are coupled through the crack model. The values of the ratio between the cracked and uncracked beam natural frequencies, the frequency ratio, are examined and are shown to follow well-defined trends with respect to the crack parameters and end mass and rotary inertia. However, the coupling between the transverse and axial vibrations is shown to be weak for the first two modes for moderate values of crack depth ratio. High crack depth ratios appear to increase the coupling effects. Low aspect ratios are expected to show strong coupling effects and further investigation is recommended using Timoshenko beam theory.     

An insight into the breathing mechanism of a crack in a rotating shaft

The time history of local flexibilities associated with a breathing crack in a rotating shaft is the concern of this paper. Considering quasi-static approximation, the deflections of a circular cross-section beam presenting a crack of different depths, due to bending or torsion loads are analyzed with the aid of a refined nonlinear contact-finite element procedure in order to predict accurately the time-variant flexibility of the fractured shaft. This method predicts the partial contact of crack surfaces, and it is appropriate to evaluate the instantaneous crack flexibilities. The bending load is applied in several aperture angles, in order to simulate a rotating load on a fixed beam. Results obtained for the rotating beam can then be used for the analysis of cracked, horizontal axis rotors. The effect of friction is also considered in the cracked area. Portions of crack surfaces in contact are predicted, the direct and the cross-coupled flexibility coefficients are calculated by applying energy principles. The numerical results compared with relevant previously published results, show high consistency.     

Variation in the stability of a rotating blade disk with a local crack defect

The effect of a near root local blade crack on the stability of a grouped blade disk is investigated in this paper. A bladed disk comprised of periodically shrouded blades is used to simulate the coupled periodic structure. The blade crack is modeled using the local flexibility with coupling terms. The mode localization phenomenon introduced by the blade crack on the longitudinal and bending vibrations in the rotating blades has also been considered. Using the Galerkin's method, the imperturbation equations of a bladed disk in which one of the blades is cracked, subject to fluctuations in the rotation speed, can be derived. Employing the multiple scales method, the boundaries of the instability zones in the mistuned turbo blade system are approximated. Numerical results indicate that an additional unstable zone is introduced near the localization frequency and the regions of unstable zones are varied with the crack size and fluctuations in disk speed.     

The influence of crack breathing and imbalance orientation angle on the characteristics of the critical speed of a cracked rotor

By analyzing the limitations of weight dominance and by taking the complicated whirl of the rotor into account, general equations of motion have been developed in case of a Jeffcott rotor with a transverse crack. The angle between the crack direction and the shaft deformation direction is used to determine the closing and opening of the crack, allowing one to study the dynamic response without assuming weight dominance. Using the new equations, the dynamic response of a cracked rotor near its critical speed has been computed via a numerical method to investigate the influence of nonlinear breathing of the crack and that of the imbalance orientation angle β on the stability, critical speed and peak response of the rotor. The results show that nonlinear breathing can improve the stability of a rotor in contrast to a rotor with an open crack, and, with a reversed imbalance (70°<β<270°), that it can reduce the vibration response in contrast to an uncracked rotor. The basic characteristics of a cracked rotor near its critical speed are similar to those of an uncracked rotor. The critical speed can be determined by measuring the rotation of the center of gravity. The critical speed of a cracked rotor is located between the natural frequencies of the fully open crack and those of the fully closed crack and depends on the imbalance orientation angle. Its value is lowest at β≈90° and highest at β≈270°. The peak in the response at the critical speed is mainly determined by the imbalance orientation angle. At β≈0° and 180°, the peak corresponds to the maximum and minimum response, respectively.     

Three dimensional vibration generators with a single rotational input

This paper presents a novel device capable of generating three-dimensional vibrations with a single-axis of rotation. The device resembles a vibration motor with an eccentric weight, but the weight is allowed to move up and down in parallel to the axis of rotations. While spinning of the eccentric weight causes lateral vibrations, vertical vibrations are excited by superposing a periodic torque on the rotary shaft. The frequency and magnitude of vertical oscillations can be independently regulated. Since the vertical natural frequency is sensitive to rotational speeds, maximum vertical oscillations can be achieved by properly adjusting the rotational speed according to the excitation frequency. Relations between the excitation torque and the vertical shaking force are examined using the frequency response for the linearized system. Numerical simulations on the original nonlinear system are conducted to verify the performance of vertical oscillations.     

A method for determining the location and depth of cracks in double-cracked beams

The influence of two transverse open cracks on the antiresonances of a double cracked cantilever beam is investigated both analytically and experimentally. It is shown that there is a shift in the antiresonances of the cracked beam depending on the location and size of the cracks. These antiresonance changes, complementary with natural frequency changes, can be used as additional information carrier for crack identification in double cracked beams. Experimental results from tests on plexiglas beams damaged at different locations and different magnitudes are found to be in good agreement with theoretical predictions. Based on the results of the present work, an efficient prediction scheme for crack localization and characterization in double cracked beams is proposed.     

THE DYNAMIC STIFFNESS MATRIX METHOD IN FORCED VIBRATION ANALYSIS OF MULTIPLE-CRACKED BEAM

The dynamic behaviour of a beam with numerous transverse cracks is studied. Based on the equivalent rotational spring model of crack and the transfer matrix for beam, the dynamic stiffness matrix method has been developed for spectral analysis of forced vibration of a multiple cracked beam. As a particular case, when the excitation frequency is close to zero, the solution for static response of beam with an arbitrary number of cracks has been obtained exactly in an analytical form. In general case, the effect of crack number and depth on the dynamic response of beam was analyzed numerically.     

Analytical modeling and vibration analysis of internally cracked rectangular plates

This study proposes an analytical model for nonlinear vibrations in a cracked rectangular isotropic plate containing a single and two perpendicular internal cracks located at the center of the plate. The two cracks are in the form of continuous line with each parallel to one of the edges of the plate. The equation of motion for isotropic cracked plate, based on classical plate theory is modified to accommodate the effect of internal cracks using the Line Spring Model. Berger?s formulation for in-plane forces makes the model nonlinear. Galerkin?s method used with three different boundary conditions transforms the equation into time dependent modal functions. The natural frequencies of the cracked plate are calculated for various crack lengths in case of a single crack and for various crack length ratio for the two cracks. The effect of the location of the part through crack(s) along the thickness of the plate on natural frequencies is studied considering appropriate crack compliance coefficients. It is thus deduced that the natural frequencies are maximally affected when the crack(s) are internal crack(s) symmetric about the mid-plane of the plate and are minimally affected when the crack(s) are surface crack(s), for all the three boundary conditions considered. It is also shown that crack parallel to the longer side of the plate affect the vibration characteristics more as compared to crack parallel to the shorter side. Further the application of method of multiple scales gives the nonlinear amplitudes for different aspect ratios of the cracked plate. The analytical results obtained for surface crack(s) are also assessed with FEM results. The FEM formulation is carried out in ANSYS.     

CRACK IDENTIFICATION IN VIBRATING BEAMS USING THE ENERGY METHOD

An energy-based numerical model is developed to investigate the influence of cracks on structural dynamic characteristics during the vibration of a beam with open crack(s). Upon the determination of strain energy in the cracked beam, the equivalent bending stiffness over the beam length is computed. The cracked beam is then taken as a continuous system with varying moment of intertia, and equations of transverse vibration are obtained for a rectangular beam containing one or two cracks. Galerkin's method is applied to solve for the frequencies and vibration modes. To identify the crack, the frequency contours with respect to crack depth and location are defined and plotted. The intersection of contours from different modes could be used to identify the crack location and depth.     

New breathing functions for the transverse breathing crack of the cracked rotor system: Approach for critical and subcritical harmonic analysis

The actual breathing mechanism of the transverse breathing crack in the cracked rotor system that appears due to the shaft weight is addressed here. As a result, the correct time-varying area moments of inertia for the cracked element cross-section during shaft rotation are also determined. Hence, two new breathing functions are identified to represent the actual breathing effect on the cracked element stiffness matrix. The new breathing functions are used in formulating the time-varying finite element stiffness matrix of the cracked element. The finite element equations of motion are then formulated for the cracked rotor system and solved via harmonic balance method for response, whirl orbits and the shift in the critical and subcritical speeds. The analytical results of this approach are compared with some previously published results obtained using approximate formulas for the breathing mechanism. The comparison shows that the previously used breathing function is a weak model for the breathing mechanism in the cracked rotor even for small crack depths. The new breathing functions give more accurate results for the dynamic behavior of the cracked rotor system for a wide range of the crack depths. The current approach is found to be efficient for crack detection since the critical and subcritical shaft speeds, the unique vibration signature in the neighborhood of the subcritical speeds and the sensitivity to the unbalance force direction all together can be utilized to detect the breathing crack before further damage occurs.     

Echolocation transmitting beam of the Atlantic bottlenose dolphin

The transmitting beam patterns of echolocation signals emitted by an Atlantic bottlenose dolphin Tursiops truncatus were measured in the vertical and horizontal planes with an array of seven hydrophones. Particular emphasis was placed on accurately verifying the animal's position on a bite-plate/tail-rest stationing device using underwater video monitoring equipment. The major axis of the vertical beam was directed at an angle of 5 degrees above the plane defined by the animal's lips. This angle was 15 degrees lower than previously measured. The vertical beam measurements indicate that the major axis of the transmitting beam is aligned with the major axis of the receiving beam. The horizontal beam was directed forward. The directivity index of 26.5 dB calculated from the beam pattern measured in both planes agreed well with previous calculation of 25.4 dB.     

Vibration based algorithm for crack detection in cantilever beam containing two different types of cracks

In this paper, a simple method for detection of multiple edge cracks in Euler–Bernoulli beams having two different types of cracks is presented based on energy equations. Each crack is modeled as a massless rotational spring using Linear Elastic Fracture Mechanics (LEFM) theory, and a relationship among natural frequencies, crack locations and stiffness of equivalent springs is demonstrated. In the procedure, for detection of m cracks in a beam, 3m equations and natural frequencies of healthy and cracked beam in two different directions are needed as input to the algorithm.     

Free vibration analysis of beams and frames with multiple cracks for damage detection

The problem of calculating the natural frequencies of beams with multiple cracks and frames with cracked beams is studied. The natural frequencies are obtained using a new method in which a rotational spring model is used to represent the cracks. For beams, dynamic stiffness matrices of order 4 are obtained in a recursive manner, according to the number of cracks, by applying partial Gaussian elimination. The Wittrick–Williams algorithm is used to compute the natural frequencies in the resulting transcendental eigenvalue problem. Published numerical examples for cracked beams are used for validation. The global dynamic stiffness matrix of a frame with multiply cracked members is then assembled. A published two bay frame example is used to evaluate the new method. The effect of changing the location of a crack in a two bay two storey frame is studied numerically, giving insight into the inverse problem of damage detection.     

Dynamic response of cracked rotor-bearing system under time-dependent base movements

Dynamic response of cracked rotor-bearing system under time-dependent base movements is studied in this paper. Three base angular motions, including the rolling, pitching and yawing motions, are assumed to be sinusoidal perturbations superimposed upon constant terms. Both the open and breathing transverse cracks are considered in the analysis. The finite element model is established for the base excited rotor-bearing system with open or breathing cracks. Considering the time-varying base movements and transverse cracks, the second-order differential equations of the system will not only have time-periodic gyroscopic and stiffness coefficients, but also the multi-frequency external excitations. An improved harmonic balance method is introduced to obtain the steady-state response of the system under both base and unbalance excitations. The response spectra, orbits of shaft center and frequency response characteristics, are analyzed accordingly. The effects of various base angular motions, frequency and amplitude of base excitations, and crack depths on the system dynamic behaviors are considered in the discussions.     

Chaotic torsional vibration of imbalanced shaft driven by a limited power supply

In this paper, torsional vibrations of imbalanced shaft driven by a limited power supply are studied. It is shown that mutual interaction of shaft and power supply may in particular result in chaotic self-oscillations that correspond to the strange attractors in the phase space of the coupled dynamical system “shaft–power supply”. In this particular model, strange attractors represent classical Lorenz and Feigenbaum attractors. Rotation characteristic of the power supply and resonance characteristic of the shaft rotational motion in one of the resonance zones are studied. It is shown that at certain intervals, these characteristics may be non-unique, which corresponds to the case of chaotic dynamics. Such non-trivial properties of the coupled system “shaft–power supply” could be used for a better understanding of complex vibrational phenomena in real applied systems such as problems related to the damping of the torsional vibrations.