Analysis of Adaptive Fuzzy Technique for Multiple Crack Diagnosis of Faulty Beam Using Vibration Signatures

This paper discusses the multicrack detection of structure using fuzzy Gaussian technique. The vibration parameters derived from the numerical methods of the cracked cantilever beam are used to set several fuzzy rules for designing the fuzzy controller used to predict the crack location and depth. Relative crack locations and relative crack depths are the output parameters from the fuzzy inference system.Themethod proposed in the current analysis is used to evaluate the dynamic response of cracked cantilever beam. The results of the proposed method are in good agreement with the results obtained from the developed experimental setup.


Introduction
Beams are one of the most commonly used structural elements in numerous engineering applications and experience a wide variety of static and dynamic loads.Cracks may develop in beam-like structures due to such loads.Considering the crack as a significant form of such damage, its modeling is an important step in studying the behavior of damaged structures.As stated, beam type structures are being commonly used in steel construction and machinery industries.Studies based on structural health monitoring for crack detection deal with change in natural frequencies and mode shapes of the beam.
An analytical study has been performed by Yang et al. [1] on the free and forced vibration of inhomogeneous Euler-Bernoulli beams containing open edge cracks.Analytical solutions are obtained for cantilever, with different end conditions to evaluate the dynamic response of the beam due to the edge crack.Orhan [2] has performed a free and forced vibration analysis of a cracked beam in order to identify the cracks in a cantilever beam.Their study reveals that free vibration analysis provides more suitable information for the detection of cracks than the forced vibration analysis.Damage in a cracked structure has been analyzed using genetic algorithm technique by Vakil-Baghmisheh et al. [3].For modeling the cracked-beam structure, an analytical model of a cracked cantilever beam has been utilized, and natural frequencies are obtained through numerical methods.A genetic algorithm is utilized to monitor the possible changes in the natural frequencies of the structure.Theoretical and experimental dynamic behaviors of different multibeams systems containing a transverse crack have been performed by Saavedra and Cuitio [4].A new cracked stiffness matrix is deduced based on flexibility, and this can be used subsequently in the FEM analysis of crack systems.Bakhary et al. [5] used Artificial Neural Network (ANN) for damage detection.In his analysis, an ANN model is created by applying Rosenblueth's point estimate method verified by Monte Carlo simulation.The results have demonstrated that the statistical ANN approach gives more reliable identification of structural damage.Friswell et al. [6] have applied genetic algorithm to the problem of damage detection using vibration data.The objective is to identify the position of one or more damage sites in a structure and to estimate the extent of the damage.A comprehensive analysis of the stability of a cracked beam subjected to a follower compressive load is presented by Wang [7].The vibration analysis on such cracked beam has been conducted to identify the critical compression load for instability based on the variation of the first two resonant frequencies of the beam.Chondros et al. [8] have developed a continuous cracked beam vibration theory for the lateral vibration of cracked Euler-Bernoulli beams with single-edge or double-edge open cracks using the Hu-Washizu-Barr variational formulation.
A new method for natural frequency analysis of beam with an arbitrary number of cracks has been developed by Khiem and Lien [9] on the basis of the transfer matrix method and rotational spring model of crack.Cam et al. [10] have performed a study to obtain information about the location and depth of the cracks in cracked beam.Experimental and simulations results obtained are in good agreement.Zheng and Kessissoglou [11] have presented a method based on finite element method for detection of crack in faulty structural member.The results obtained from the proposed method are validated using experimental analysis.Tada et al. [12] have provided the basis for computation of compliance matrix for damage detection following fracture mechanics theory.Sekhar and Prabhu [13] have derived a method for crack detection in a cracked shaft using finite element analysis using correct expression for strain energy release rate function.Hossain et al. [14] have presented an investigation for comparative performance of intelligent system-like genetic algorithms (GAs) and adaptive neurofuzzy inference system (ANFIS) algorithms for identification of fault in an active vibration control (AVC) system.A comparative performance of the proposed method is presented and discussed through a set of experiments.Ranjbaran et al. [15] in their paper have formulated a method for vibration analysis of a beam assuming the beam to be nonuniform.The vibration characteristics of the beam are computed and the results are compared with other methods.Zhou and Biegalski [16] have proposed a method to analyze the vibration signatures of a deck truss bridge with cracks at the gusset plate connecting the lower lateral bracing.In-service monitoring has been performed to measure the vibration properties of the truss and the lateral bracing members to avoid resonance with the excitation frequency.Wada et al. [18] have proposed a fuzzy control method with triangular type membership functions using an image processing unit to control the level of granules inside a hopper.They have stated that the image processing technique can be used as a detecting element, and with the use of fuzzy reasoning methods, good process responses are obtained.Pawar et al. [17] have used a genetic fuzzy system to identify the crack depth location in a composite matrix cracking model.As described by them, the genetic fuzzy system combines the uncertainty characteristics of fuzzy logic with the learning ability of genetic algorithm.Parhi [19] has designed a mobile robot navigation control system using fuzzy logic.Fuzzy rules embedded in the controller of a mobile robot enable it to avoid obstacles in a cluttered environment that includes other mobile robots.
In the present study, a finite element model for a cracked beam element is developed, and the results from theoretical and finite element analyses have been used to set the fuzzy rules.The fuzzy rules are used for designing the fuzzy inference system based on Gaussian membership function which is subsequently applied for prediction of damage in a faulty structure.The theoretical, finite element and fuzzy analysis are done to study the response of a cantilever beam in the presence of cracks.Theoretical results are compared with the experimental, fuzzy, and FEM results.A close agreement between the results has been observed.

Local Flexibility of a Cracked Beam under Bending and
Axial Loading.The presence of a transverse surface crack of depth " 1 " and " 2 " on beam of width "" and height "" introduces a local flexibility, which can be defined in matrix form, the dimension of which depends on the degrees of freedom.Here, a 2 × 2 matrix is considered.A cantilever beam is subjected to axial force ( 1 ) and bending moment ( 2 ), shown in Figure 1(a), which gives coupling with the longitudinal and transverse motion.The cross-sectional view and the front view of the beam are shown in Figures 1(b) and 3, respectively.
The strain energy release rate at the fractured section can be written as follows [12]: where 1 and  2 are the stress intensity factors of mode I (opening of the crack) for loads  1 and  2 , respectively.The values of stress intensity factors from earlier studies [12] are where expressions for  1 and  2 are as follows: Let   be the strain energy due to the crack.Then, from Castigliano's theorem, the additional displacement along the force   is    The strain energy will have the form where  =   / is the strain energy density function.
From ( 1) and ( 4), thus we have The flexibility influence coefficient   will be, by definition, and can be written as From ( 9), calculating  11 ,  12 (= 21 ), and  22 , we get The local stiffness matrix can be obtained by taking the inversion of compliance matrix.That is, Figure 2 shows the variation of dimensionless compliances to that of relative crack depth.

Analysis of Vibration Characteristics of the Cracked Beam.
A cantilever beam of length "", width "", and depth "", with a crack of depth " 1 " and " 2 " at a distance " 1 " and " 2 ", respectively, from the fixed end, is considered (shown in Figure 1).Taking  1 (, ),  2 (, ), and  3 (, ) as the amplitudes of longitudinal vibration for the sections before, in-between, and after the crack,  1 (, ),  2 (, ), and  3 (, ) are the amplitudes of bending vibration for the same sections.
The normal function for the system can be defined as where  = /,  = /,  = /,  1 =  1 /,  2 =  2 /,   = /  ,   = (/) 1/2 ,   = ( 2 /  ) 1/2 ,   = (/) 1/2 , and  = ,   , ( = 1, . . ., 18).Constants are to be determined from boundary conditions.The boundary conditions of the cantilever beam in consideration are At the cracked section Also, at the cracked section  1 , we have Multiplying both sides of the previous equation by / 11  12 , we get Similarly, Multiplying both sides of the previous equation by / 2  22  21 , we get where    The results of the theoretical analysis for the first three mode shapes for uncracked and cracked beams are shown in Figure 4.

Analysis of Cracked Beam Using Finite Element Method (FEM)
In the following section, FEM is analyzed for vibration analysis of a cantilever cracked beam (Figure 5).The relationship between the displacement and the forces can be expressed as where overall flexibility matrix  ovl can be expressed as The displacement vector in ( 23) is due to the crack.The forces acting on the beam element for FEM analysis are shown in Figure 5.
Under this system, the flexibility matrix  intact of the intact beam element can be expressed as where  intact = [ / 0 0 / ].The displacement vector in (25) is for the intact beam.The total flexibility matrix  tot of the cracked beam element can now be obtained by Through the equilibrium conditions, the stiffness matrix   of a cracked beam element can be obtained as follows [13]: where  is the transformation matrix and is expressed as The results of the finite element analysis for the first three mode shapes of the cracked beam are compared with that of the numerical analysis of the cracked beam and are presented in Figure 6.

Analysis of the Fuzzy Controller
The fuzzy controller developed has got six input parameters and two output parameters.The linguistic terms used for the inputs are as follows: relative first natural frequency = "fnf "; relative second natural frequency = "snf "; relative third natural frequency = "tnf "; average relative first mode shape difference = "fmd"; average relative second mode shape difference = "smd"; average relative third mode shape difference = "tmd".
The linguistic terms used for the outputs are as follows: first relative crack location = "rcl1" and first relative crack depth = "rcd1"; second relative crack location = "rcl2" and second relative crack depth = "rcd2".
The fuzzy controller used in the present text is shown in Figure 7(a).The Gaussian membership functions are shown pictorially in Figure 7(b).The linguistic terms for the Gaussian membership functions, used in the fuzzy controller, are described in Table 1.

Rule no. 6 of
where  = 1 to 10,  = 1 to 10,  = 1 to 10,  = 1 to 10,  = 1 to 10, and  = 1 to 10 because "fnf, " "snf, " "tnf, " "fmd, " "smd, " and "tmd" have ten membership functions each.From expression (29), two sets of rules can be written: If fnf is fnf  and snf is snf  and tnf is tnf  and fmd is fmd  and smd is smd  and tmd is tmd  , then rcd1 is rcd1  and rcd2 is rcd2  ; If fnf is fnf  and snf is snf  and tnf is tnf  and fmd is fmd  and smd is smd  and tmd is tmd  , then rcl1 is rcl1  and rcl2 is rcl2  .
According to the usual fuzzy logic control method [19], a factor   is defined for the rules as follows: where freq  , freq j , and freq  are the first, second, and third relative natural frequencies of the cantilever beam with crack, respectively; moddif  , moddif  , and moddif  are the average first, second, and third relative mode shape differences of the cantilever beam with crack, respectively.By applying the composition rule of inference [19], the membership values of the relative crack location and relative crack depth, (location) rcl and (depth) rcd ( = 1, 2), can be computed as The overall conclusion by combining the outputs of all the fuzzy rules can be written as follows: The crisp values of relative crack location and relative crack depth are computed using the centre of gravity method [19] as relative crack location

Experimental Setup
Experiments are performed to determine the natural frequencies and mode shapes for different crack depths on Aluminum beam specimen (800 × 38 × 6 mm).The experimental setup is shown in Figure 9.The amplitude of transverse vibration at different locations along the length of the Aluminum beam is recorded by positioning the vibration pickup and tuning the vibration generator at the corresponding resonant frequencies.These results for first three modes are plotted in Figure 10.Corresponding numerical results are also presented in the same graph for comparison.

Discussion
The results from fuzzy controller and the information obtained from theoretical, finite element and experimental analysis of the cracked cantilever beam are depicted later.Fuzzy logic systems promise an efficient way for damage assessment in a cracked structure.They are able to treat uncertain and imprecise information; they make use of knowledge in the form of linguistic rules.At first, the theoretical expressions for the cracked and uncracked beams are developed for calculation of natural frequencies, mode shapes following correct expression of strain energy release rate.From Figure 2, it is observed that the compliances increase with the increase in relative crack depth.Finite element analysis is being performed on the cracked beam element (Figure 5) to find out the vibration characteristics.The comparison between results obtained from theoretical analysis for cracked and uncracked beams with the magnified view at the vicinity of crack location is presented in Figure 4. Results from finite element analysis (FEA) and numerical analysis are compared for both uncracked and cracked beams and are shown in Figure 6.For validation of the information obtained from various methodologies for analysis of the cracked cantilever beam, an experimental setup is developed as shown in Figure 9. Experiments are performed on the Aluminum beam specimen (800 × 38 × 6 mm) to estimate the first three mode shapes which are compared with the mode shapes obtained from the analysis as mentioned earlier.
The vibration signatures (natural frequencies, mode shapes) are used to establish the fuzzy rules for designing fuzzy controller based on Gaussian membership functions that are depicted in Figures 7(a) and 7(b).The linguistic terms of fuzzy rules in the present fuzzy controller are given in Table 1.Some of the actual rules made for the fuzzy controller of the present investigation are listed in Table 2.The outputs of the Gaussian fuzzy controller obtained by activating rule 6 and rule 19 from Table 2 are presented in Figure 8.The mode shapes obtained from experimental and numerical, finite element, fuzzy analysis for cracked and uncracked beams are compared graphically in Figure 10.Some of the predicted outputs of the developed fuzzy controller and corresponding numerical finite element, and experimental results are presented in Table 3.It is observed that the results of all analyses are in good agreement.

Conclusions
In the present section, results obtained from the different analyses have the following conclusions.
As clear cut deviation of the mode shapes, natural frequencies can be detected for the cracked and uncracked beams at the vicinity of crack location.The comparison of results derived from theoretical and finite element method (FEM) for the cracked structure shows a good agreement.The fuzzy controller developed with Gaussian membership functions is designed with the help of the vibration signatures obtained from numerical and finite element analyses.
The first three relative natural frequencies and the first three mode shapes in dimensionless forms are the input parameters to the fuzzy inference system.Relative crack location and relative crack depth are the outputs from the system.The comparison of results between experimental and fuzzy analyses shows a close agreement.From the comparison of results, it is observed that the developed fuzzy inference system can predict the relative crack location and relative crack depth in a faster and accurate way, thereby decreasing a considerable amount of computational time.The proposed method can be used as an online condition monitoring tool, and in the future, hybrid technique can be developed for a faster and more efficient way for fault detection in the domain of dynamic vibrating structure.

Figure 5 :
Figure 5: View of a crack beam element subjected to axial and bending forces.

Figure 7 :
Figure 7: (a) Fuzzy controller.(b1) Membership functions for relative natural frequency for first mode of vibration.(b2) Membership functions for relative natural frequency for second mode of vibration.(b3) Membership functions for relative natural frequency for third mode of vibration.(b4) Membership functions for relative mode shape difference for first mode of vibration.(b5) Membership functions for relative mode shape difference for second mode of vibration.(b6) Membership functions for relative mode shape difference for third mode of vibration.(b7a) Membership functions for relative crack depth 1. (b7b) Membership functions for relative crack depth 2. (b8a) Membership functions for relative crack location 1. (b8b) Membership functions for relative crack location 2.

Table 1 :
Description of fuzzy linguistic terms.

Table 2
Resultant values of relative crack depth and relative crack location when rules 6 and 19 of Table 2 are activated.If fnf is fnf  and snf is snf  and tnf is tnf  and fmd is fmd  and smd is smd  and tmd is tmd  , then rcl1 is rcl1  and rcd1 is rcd1  and rcl2 is rcl2  and rcd2 is rcd2  ,

Table 3 :
Comparison of results between fuzzy Gaussian controller, numerical FEM analysis, and experimental setup., average relative first mode shape difference, average relative second mode shape difference, and average relative third mode shape difference.The outputs from the fuzzy controller are relative crack depth and relative crack location.Twenty numbers of the fuzzy rules out of several hundreds of fuzzy rules are being listed in Table2.The fuzzy controller results when rule 6 and rule 19 are activated from Table2are shown in Figure8.