Ali Bakhshizade^{1} , Mohammad Reza Ashory^{2}
^{1}Department of Mechanical Engineering, Roshdedanesh Institute of Higher Education, Semnan, Iran
^{2}Modal Analysis Lab, School of Mechanical Engineering, Semnan University, Semnan, Iran
^{1}Corresponding author
Mathematical Models in Engineering, Vol. 1, Issue 2, 2015, p. 96101.
Received 17 September 2015; accepted 18 October 2015; published 30 December 2015
Copyright © 2015 JVE International Ltd.
Detection of structural damages using the vibration test data has been investigated by many researchers in the last decades. Identification of the defect in the early stages can prevent the unpredicted failure of structure. For beamlike structures, curvature techniques, e.g., mode shape curvature and flexibility curvature have been applied to localize the damage. In this paper, root mean square error using operational deflection shape curvature is used for damage detection. Differential Quadrature Method (DQM) is implemented to obtain the curvatures of the operational deflection shape. To demonstrate the effectiveness of method, the finite element model of a cantilever beam is used. The model is excited and the responses are measured in the simulated test. Enhanced frequency domain decomposition (EFDD) is used to extract the modal parameters. Also Mass change method is used to scale the mode shapes through which the scaled ODS parameters are computed. Using the RMSE the location of the damage is identified. The results show that this approach is more effective than the previous methods
Keywords: operational deflection shape, differential quadrature, RMSE.
Damage in the structural and mechanical systems is defined as changes in the material and/or geometric properties of these systems, including changes in the boundary conditions and system connectivity, which adversely affect the current or future performance of these systems” [1]. Vibration based techniques have been widely expanded in few decades. These methods localize the flaws and defects in structures based on changes in the modal parameters (natural frequencies, mode shapes and damping ratio).
These methods can be classified in many ways but the damage classification system presented by Rytter [2] is adopted here, which defines four levels of damage assessment:
1) (detection): determination that whether the damage is present in the structure;
2) (localization): determination of the geometric location of the damage;
3) (quantification): determination of the severity of damage;
4) (prediction): prediction of the remaining service life of the structure.
The aim of investigations in the last 20 years has been the global and online nondestructive evaluation. Amongst the techniques under development one can outline those using mechanical vibrations, smart materials, image processing, optical fibers, novelty based and statistical information. Vibrationbased methods have been applied successfully on large and complicated structures such as civil engineering structures (mostly bridges), aerospace structures and composite structures. During the last decades, research on vibrationbased structural damage identification methods has been expanding rapidly.
Some of the existing methods in vibration are based on the comparison between the parameters extracted from experimental observations and numerical models of the structure; there are also methods based on the comparison between parameters identified from experimentalonly observations [3].
The method presented in this paper is in the category of methods that use the change in the operational deflection shapes (ODSs) to detect, locate and relatively quantify the damage.
The Frequency Domain Decomposition (FDD) method is proposed by Brincker et al. [4, 5]. In this method first the Power Spectral Density (PSD) of the response is computed. Then the singular value decomposition is applied to obtain the modal parameters of system. The relation between the inputs and outputs is given by [5]:
where ${G}_{xx}$ is the PSD matrix of input, ${G}_{yy}$ is the PSD matrix of output and $H\left(j\omega \right)$ is the FRF matrix. If the input force is assumed to be a white signal, the PSD of the output can be given by:
where ${d}_{k}$ is a scalar, ${\varphi}_{k}$ is the $k$th mode shape vector, ${\gamma}_{k}$ is the modal participation vector, ${\lambda}_{k}$ is the $k$th complex resonance frequency, $\left(\stackrel{}{\bullet}\right)$ sign stands for the complex conjugate and $t$ indicates the transpose of a matrix.
The PSD of response in each frequency can be decomposed to the singular values and singular vectors using the following equation:
where ${U}_{i}$ is the $i$th singular vector, ${S}_{i}$ is the $i$th singular value, ${\omega}_{i}$ is the $i$th frequency and $H$ indicates the complex conjugate of a matrix.
The peaks of first singular values of system correspond to the natural frequencies of system. The singular vectors corresponding to the peaks of first singular values approximate the mode shapes.
In output only modal analysis, the mode shapes are not scaled. Mass change method is a method for scaling of mode shapes.
For a natural frequency ${\omega}_{n1}^{}$, and the corresponding mode shape $\left\{{\varphi}_{1}\right\}$ of a MDOF system, if the masses of vibrating system at some or all of the degrees of freedom are changed, the $r$th natural frequency and the corresponding mode shapes are changed to ${\omega}_{n2}^{}$ and $\left\{{\varphi}_{2}\right\}$ respectively. Considering small changes of mass matrix, it can be assumed that the mode shapes change negligibly [6]. Therefore:
The relation between the unscaled mode shape $\left\{\psi \right\}$ and the scaled mode shape $\left\{\phi \right\}$ is:
where $\alpha $ is the scaling factor.
It is shown in [6], that the best scaling factor is obtained when the mode shapes are scaled to the length and both the modified and unmodified mode shapes are used:
where $\left[\mathrm{\Delta}M\right]$ is the mass change matrix, $\left\{{\psi}_{1}\right\}$ and $\left\{{\psi}_{2}\right\}$ are the unscaled mode shapes before and after mass change.
Frequency response function can be obtained using the following equation [7]:
where $\eta $ is the damping ratio and $\alpha $ the receptance matrix. Damping ratio can be obtained from EFDD method. The column vector, k, of the receptance matrix, $\left\{{\alpha}_{k\left(\omega \right)}\right\}$, is the ODS at frequency $\omega $. Therefore, the ODS describes the shape exhibited by the structure at each excitation frequency $\omega $, given by the responses and normalized by the applied forces.
When a structure is damaged its stiffness and damping alter and, as a result does the receptance matrix:
where the superscript $d$ stands for “damaged”. It is reasonable to assume that the smaller the degree of correlation between the column vectors (ODSs), $\left\{\alpha j\left(\omega \right)\right\}$ and $\left\{d\alpha j\left(\omega \right)\right\}$, the larger the damage. Root mean square error using the spectral strain energy was first studied by Bayissa and Haritos [8]. In this paper is used the curvature of ODS instead of spectral strain energy in root mean square error criterion:
where $N$ is the number of measurement points. But central difference can not present the curvature accurately, two nodes are missing in the approximation and it needs to many points of measurement. In this work, differential quadrature method was used to obtain the curvature.
The idea of using differential quadrature method is to quickly compute the derivative of a function at any grid point within its bounded domain by estimating a weighted linear sum of values of function at a small set of points related to domain. According to this method nth derivative of $f\left(x\right)$ can be approximately obtained by:
In order to determine the weighting coefficients a specific test function is required. One of the best approaches has been proposed by Quan and Chang [9, 10]. They used the Lagrange interpolation function as test function and determined the weighting coefficients on this basis. Weighting coefficients for the firstorder derivatives are:
To evaluate the weighting coefficients of higher order derivatives recurrence formulae are derived as:
A study based on finite element model of a cantilever beam was conducted to evaluate the feasibility of applying curvature of ODS in RMSE criterion to identify the crack of a beamstructure. The geometry of the beam and its material properties are as follows: Young’s modulus $E=$206 GPa, density $\rho =$7800 kg/m^{3}, length $l=$800 mm, Cross sectional area $A=w\times h=$ 32 mm × 16 mm. The total number of elements is 8 and two damage cases are studied. Damage characteristics have been described in Table 1.
Table 1. Damage scenarios
${l}_{1}$ (mm)

Element number

$a$ (mm)


Case I

650

6

0.8

Case II

450

3

0.6

Fig. 1. Geometry of cantilever beam with a crack
Damage is simulated with a method that is proposed by Nahvi and Jabbari [11]. In this method the length of damaged element is smaller than the other elements.
The application of ODS curvature using DQM compared to the method presented in [8] for the above cases have been shown in the Figs. 25.
Fig. 2. Method in [8] (Case I)
Fig. 3. Proposed method (Case I)
Fig. 4. Method in [8] (Case II)
Fig. 5. Proposed method (Case II)
For case (I), the proposed method shows that the damage is between nodes 6 and 7 (Fig. 3). The result is concluded from Fig. 2. However, the difference in Fig. 3 is more indicative compared to the result in Fig. 2. The same conclusion can be found for case II by comparing Figs. 4 and 5.
In this study, a new method is developed to improve the damage localization of a beam structure which has been numerically investigated. Results show that the proposed approach can be used to localize the damage in the structure without any analytical model. When the number of measurement points is limited, DQM method has the good accuracy and the less time of computation. In comparison with method in [8] that need, a large number of measurement points to obtain an accurate result, the proposed method has inherent advantages of finding cracks where the method suggested in [8] is relatively weak.