A Novel Methodology Using Simplified Approaches for Identification of Cracks in Beams 1

In this paper, natural frequency based forward and inverse methods are proposed for identifying multiple cracks in beams. Forward methods include simplified definition of the natural frequency drops caused by the cracks. The ratios between natural frequencies obtained from multi-cracked and un-cracked beams are determined by an approach that uses the local flexibility model of cracks. This approach does not consider nonlinear crack effects that can be easily neglected when the number of cracks is not excessive. In addition, an expression, which removes the necessity of repeating natural frequency analyses, is given for identifying the connection between the crack depths and natural frequency drops. These simplified approaches play crucial role in solving inverse problem using constituted crack detection methodology. Solution needs a number of measured modal frequency knowledge two times more than the number of cracks to be detected. Efficiencies of the methods are verified using the natural frequency ratios obtained by the finite element package. The crack detection methodology is also validated using some experimental natural frequency ratios given in current literature. Results show that the locations and depths ratios of cracks are successfully predicted by using the methods presented.


Latin American Journal of Solids and
crack identification methods take great attention throughout the researchers studying nondestructive testing and evaluation methods.
Changes in vibration characteristics can be observed by using modal parameters such as natural frequency and mode shape.In contrast to the mode shapes, natural frequencies can be measured easily, and they are not seriously affected by the experimental errors.In connection with this, crack detection systems, which use frequencies, can be inexpensive, non-invasive and automated for avoiding subjective operator differences.Therefore, frequency based inverse methods supported by the theoretical vibration models have been frequently proposed up to date.Dimarogonas (1996) and Doebling et al. (1998) present overviews to the methods examining the theoretical changes in dynamic behaviours and their agreement with measured vibration responses.Frequency based crack identification methods are summarised by Salawu (1997).Detailed review on vibration based identification of multiple cracks is given by Sekhar (2008).Recently, Jassim et al. (2013) review the studies on vibration analysis of damaged cantilever beams.
The identification of a single transverse crack in a beam is popularly studied using the lowest three natural frequencies which can be easily obtained (Chen et al. 2005;Chinchalkar 2001;Kim and Stubbs 2003;Liang et al. 1991;Nandwana and Maiti 1997;Owolabi et al. 2003).However, simultaneous detection of crack parameters is much more involved and complex than the identification of single crack.In many cases, supports of theoretical vibration analyses are inevitable for the detection.In addition, multiple crack detection can require the knowledge of additional parameters and thus many measurements in several test conditions.Number of unknown parameters is generally decreased by removing some of them which can have negligibly small effects.Theories of most studies are based upon the analytical method including local flexibility model which uses the conditions of compatibility and continuity at crack locations (Bakhtiari-Nejad et al. 2014;Caddemi and Calio 2009;Dado 1997;Khiem and Lien 2001;Nandwana and Maiti 1997;Ostachowicz and Krawczuk 1991;Shifrin and Ruotolo 1999).Local flexibility model neglects the distributed energy effects of crack and accumulates the effects of dynamical changes into rotational springs considered at the crack locations.In classical solution method that use local flexibility model, the equation set satisfying the boundary conditions at two ends of the beam is expanded with four new equations of continuity and compatibility conditions for each crack.With this method, to construct linear system is not a simple task for the beams with n cracks.Shifrin and Ruotolo (1999) extend this base approach by defining 2 + n equations for simplifying the vibration analysis of the beams with n cracks.However, their method can have the lack of exact- ness due to the some simplifications and linearization done for decreasing the matrix size.Some other methods, including exponentially decaying crack disturbance functions, have been proposed to develop vibration equations for continuous models (Christides and Barr 1984;Mazanoglu et al. 2009;Yang et al. 2001).When these methods are used in the inverse analyses, interaction effects of the cracks can be required to be neglected for being capable of detecting cracks.Accuracy of the vibration models is significantly important for the success of inverse methods presented for crack detection.
In the last two decades, several papers have been published to solve inverse problem of detecting multiple cracks with the knowledge of natural frequency changes.However, the number of these studies is still less than the analyses addressing the forward problem.Mazanoglu and Latin American Journal of Solids and Structures 12 (2015) 2460-2479 Sabuncu (2012) present an algorithm that uses natural frequencies for detection of cracks in beams and a process which minimises the errors in experimental results.Gillich and Praisach (2014) study on a crack detection method based upon natural frequency changes and a procedure for determining accurate natural frequencies in measurement.Khiem and Toan (2014) present Rayleigh quotient based explicit expression providing a tool for calculating natural frequencies of the beam with arbitrary number of cracks.A method that combines the frequency measurements and vibration modelling using transfer matrix method is presented by Patil and Maiti (2003).They observe notable amounts of errors like 10% in the prediction of locations and sizes of cracks.Khiem and Lien (2004) use the natural frequencies obtained by the dynamic stiffness matrix method and formulate the multiple crack detection as a non-linear optimization problem.Morassi and Rollo (2001) present a technique, which uses the changes in the first three natural frequencies, for a simply supported beam with two cracks having equal severity.Although the method is verified numerically, notable errors are seen in their predictions.Douka et al. (2004) use the antiresonance changes, complementary with natural frequency changes, in a prediction scheme for identification of two cracks in beams.In their work, crack detection is based upon the availability of many experimental results.Lee (2009) presents a simple method for detecting n cracks using 2n natural frequencies by means of the finite element and the Newton-Raphson methods.Their method reveals quite satisfactory results provided proper initial guesses are made for convergence of the solution.Ruotolo and Surace (1997) propose a solution procedure employing a genetic algorithm and the results of the finite element model for the detection of multiple cracks in beams.Krawczuk (2002) uses the wave propagation approach combined with an iterative searching strategy including two methods for damage detection in beam-like structures.
This paper presents simplified methods for determining natural frequency drops due to the multiple cracks in beams and a methodology for identifying location and depth ratio parameters of the cracks.First of all, theoretical background about the vibration of cracked beam is outlined.Theoretical simplifications are demonstrated to find natural frequencies of cracked beams.After that, methodology of crack identification, which employs the natural frequencies obtained by the forward analyses, is detailed.Then, several cracked beam scenarios are considered for discussing the results that are verified using the natural frequency ratios obtained by both the commercial finite element package (ANSYS©) and some experiments given by Kim and Stubbs (2003), Mazanoglu and Sabuncu (2012), Ruotolo and Surace (1997).Finally, achievements and shortcomings of the methods are clearly stated.

THEORETICAL BACKGROUND
Free bending vibration of a uniform Euler-Bernoulli beam is identified by following well known differential equation.where, E, I, A, and ρ represent modulus of elasticity, area moment of inertia, cross-section area, and density respectively.Flexural displacement is symbolised by w, and variables z, t are the position along the beam length and the time respectively.When this equation is solved by separating the variables of z and t, mode shape of the beam, ) (z W , is observed in the following solution form (Rao 1995).

0
where 1 C , 2 C , 3 C , and 4 C are the coefficients of harmonic and hyperbolic terms.Frequency parameter, β , which depends upon the natural frequency, ω , is written as: An equation set is formed by using mode shape function and its derivatives satisfying the end conditions of the beam.The equations appropriate for some general end conditions are given in Table 1.
In example, for a cantilever beam, the matrix obtained by the harmonic and hyperbolic terms of functions is formed as follows: where L is the length of beam.Zero determinant of the matrix, 0 M , gives singular values as uncracked beam's natural frequencies, 0 ω .

Vibration of the Beam with Cracks
When the cracks exist on a beam, local flexibility changes should be identified at the crack locations.Local flexibility changes are simulated by the rotational springs, which are joints of the sections separated by the cracks, as shown in Figure 1.Cracks are modelled as slots whose depths are assumed to be unchanged along the width of the beam.Crack widths are considered to be negligibly small and thus the cracked beam's vibration identification problem is handled by the Latin American Journal of Solids and Structures 12 (2015) 2460-2479 assumption of no mass loss from the beams.In classical approaches, the existence of n cracks requires the expression of n local flexibility changes for connecting 1 + n sections.Vibration form of each section can be expressed by appropriate harmonic and hyperbolic terms of the function written in Eq. ( 2).Connection at the crack location is provided by the continuity conditions with negligible effects of crack width.Deflection, bending moment and shear force are assumed to be equal at right hand and left hand sides of the crack as follow: In addition, compatibility condition relates bending moment with the difference of slopes between both sides of the crack as represented in following equation.
[ ] where Here, i k represents the local rotational stiffness caused by th i crack, and it is described by the theoretical and experimental expressions of the fracture mechanics (Ostachowicz and Krawczuk 1991;Tada et al. 1973) as follows: ( ) Solids and Structures 12 (2015) 2460-2479 where, b, h and a symbolise width and height of the beam and depth of the crack respectively.
is known as flexibility compliance function of th i crack that is formulated for the rectangular beam as follows: In conventional procedure, the equation set having size, 4 4 + n , is formed by 4n equations of continuity and compatibility conditions and 4 equations of the end conditions (Nandwana and Maiti 1997;Ostachowicz and Krawczuk 1991).Large sized matrix shaped by harmonic and hyperbolic terms of the equation set must be singular for determining natural frequencies.Instead, a simplified approach is presented in linear manner to find natural frequencies of the beams with cracks.

An Approach to Obtain Multi-cracked Beam's Natural Frequencies
Suppose that a cantilever beam considered as an example has a crack at the location 1 z .

8×
matrix, 1 M , is formed by the terms of functions providing continuity and compatibility conditions given in Eqs.(5, 6) and boundary conditions tabulated in Table 1.Here, subscript " 1 " is used for defining the parameters in cases of beams with one crack.Reduced form of the matrix is constructed as follows: Since the functions are continuous, infinite number of β can be available for producing zero determinant of 1 M .As is known, each β is correspond of a different natural frequency whose number is now represented by the superscript (m) .For each mode, single cracked beam's natural frequency,  and Structures 12 (2015) 2460-2479 Change in frequency can be examined by natural frequency ratio, . When the second crack initiates from the location, 2 z , m th mode natural frequency starts to decrease from ω .This proves that individual effect of the second crack should be equal to Thus, the natural frequency ratio in case of the existence of two cracks can be computed for m th vibration mode as follows (Mazanoglu and Sabuncu 2012): Minor errors arising from the linear system assumption are indicated by the approximate equivalence in the equation.If the crack at the location, 2 z , initiates first, individual effect of this crack should remain same.It is for this reason that the following expression is accurate.
In general aspect, Eqs.(11, 12) can be arranged for the cases of n cracked beams as follow: ω symbolise the natural frequencies for the beam with n cracks and solely th n crack respectively.As a result of the substitution of Eq. ( 14) into Eq.( 13), following expression is obtained for the natural frequency ratio of a beam with n cracks.
Each natural frequency ratio multiplied in Eq. ( 15) is obtained by making singularity analysis of 8 8× matrix shaped for every single cracked beam case.It is assumed here that each crack has a specific capacity of energy storing varying with the location and depth ratio.Since the energies consumed due to the cracks are locally stored on the springs in local flexibility model, each crack has independent influences on beams vibration.In the present work, this assumption is verified up to three cracks that need natural frequency ratios for six vibration modes.Existence of more cracks ( 3 > n ) requiring more natural frequency ratios (2n) is not taken into consideration since the errors of the frequency ratio multiplication increase as the order of natural frequency increases (Gillich et al. 2012).This is because the following crack detection methodology is effective for detecting three cracks at most.Structures 12 (2015) 2460-2479

A METHODOLOGY FOR DETECTING CRACKS
Novel crack detection procedure is now introduced step by step.
Step 1: Experimentally found modal frequency drop ratios, ) (m inp d , are employed as input data of proposed crack detection method.Instead of direct use of the natural frequency ratios, ) (m inp r , natu- ral frequency drop ratios are employed, since the methodology works better with this parameter. Step 2: Theoretical natural frequency drop ratios ) ( 1 m d should be calculated to have a prior knowledge about the frequency of the beam with a crack.Conventionally, this may require repeating singularity analysis for all crack depths to prepare prediction scheme of corresponding natural frequency drop ratios.This is quite inconvenient and time-consuming procedure.To avoid this, a simplified formula is proposed here to calculate the changes in natural frequency drop ratios reasoned by various crack depths at one position.For small cracks, fractional changes in modal strain energy are equal to the natural frequency drop ratios (Gillich and Praisach 2014;Gudmunson 1982;Kim and Stubbs 2003).By neglecting minor effects of crack based changes in second derivative of mode shape, approximate strain energy changes can directly be expressed by the flexibility increases which are the function of as given in Eq. ( 8).Once the drop ratio, arising from the crack with different depth, a, can be calculated by the following expression.Step 3: For the reference crack depth, theoretical natural frequency drop ratios are determined by shifting the crack location along the z coordinate axis of beam.By fitting a curve on these natural frequency drop ratios, reference polynomial functions, ) , ( , are constituted in following form: Eq. ( 17) is the set of reference polynomials constructed using natural frequency drop ratios for six modes of vibration ( 6 ,...., 2 , 1 = m ) that are necessary to inspect three cracks at most.k κ symbolises the coefficients of polynomial terms indexed by k.These coefficients can vary with the depth of crack and the geometric properties of beam.Values of coefficients and number of polynomial terms also depend upon the mode of vibration considered.Higher vibration modes require larger number of terms in the function.
Step 4: Reference polynomial functions should be generalised for determining natural frequency drop ratios caused by the crack with any depth and location.For this purpose, Eq. ( 16) is rearranged in the form of polynomial function as follows: symbolises the function set of natural frequency drop ratios varying with the location and depth of a crack.This set of polynomials satisfy the natural frequency drop ratios ( ) , ( ) when the unknown crack parameters, locations ( i z ) and depths ( i a ), are accurately substituted.
Step 5: It is definite that measured modal frequency ratios employed as input of the method should be equal to the theoretical natural frequency ratios.Accordingly, following equation set should be satisfied.
Eq. ( 19) is written in the form of the set of polynomial functions using the relations of Eq. ( 20) presenting the nonlinear equation set for six vibration modes allows to find unknown crack parameters, Step 6: In the final step, unknown crack parameters are initially predicted as , and then an optimisation problem of finding accurate crack parameters is solved by minimising following residual function set In this work, trust-region dogleg algorithm introduced in the package (Matlab © ) is carried out as a solution method of the equation set (Nocedal and Wright 1999).The beam is also modelled by the finite element package (ANSYS © ).In the program, cracks with the thicknesses, 0.4 mm, are modelled as slots formed by subtracting thin transverse blocks from the beam.Solid element called 'Solid95', which includes 20 dynamic nodes, is used as meshing element.Default edge length of the element is set to 5 mm using the "esize" command.Good convergence to exact results requires frequent meshing in the vicinity of any discontinuities.This necessity is provided by using the "smrtsize1" command that is the most refined mesh option of the free meshing procedure.Figure 3 shows meshing view around a slot considered as a crack.At result, natural frequencies are obtained by using "modal analysis" as the analysis type.It should be noted that variations of crack location and crack size, which lead to change of total number of elements, have negligible effects on sensitivity of the computations.

Verification of Forward Methods
The approaches formulated by Eq. ( 15) and Eq. ( 16) are verified by considering the example cases given in Table 2.It is seen that the crack modelled in Case 1 is fixed in Case 4 as one of the two cracks.Triple cracked beam scenario including these cracks is also considered in Case 10.Similar relationships are established throughout Cases 2,5,11 and Cases 3,6,12.Additionally, ex- traordinary cases like the existence of small or large distance between two cracks with incipient or advanced severity are handled in Cases 7,8,9.Eq. ( 15) is validated throughout Cases 4-12.Table 3 demonstrates the agreement between the natural frequency ratios computed by large sized matrix solution and the finite element package.In addition, for the beam cases with two and three cracks, comparison of the results obtained by large sized matrix solution and frequency ratio multiplication method are given in Table 4. Fractal numbers equal or very close to unity show the success of frequency ratio multiplication method.The largest deviations from unity are obtained in Case 7 due to the rise of nonlinearity with two cracks having advanced depth ratio as 0.4.The maximum deviation observed as 1.048 is still admissible for approximate detection of the cracks.It should be noted as an advantage of frequency ratio multiplication that the determination of singular values for the 8 8× matrix is easier than the solution of larger sized matrix con- structed for multi-cracked beams.Eq. ( 16), which states the change of natural frequency drop ratios caused by various depths of the cracks, also needs verification.Figure 4 shows natural frequency drop ratios caused by the crack whose location is shifted along the beam.It is seen that the results obtained by proposed formula- tion employing only the knowledge of reference crack depth ratio, 0.25, present good agreement with the results of matrix solution repeatedly done for the depth ratios, 0.1 and 0.4.When the crack depth ratio is 0.4, some minor discrepancies appear between the results of the matrix solution and presented formulation.These deviations increase with rising crack depth differences.This is because the selection of the reference at the average level of crack depth ratio.

Application of Crack Detection Procedure
In this subsection, proposed crack detection procedure is implemented on the cantilever beam whose properties are given in Section 4.1.First of all, a database including the natural frequency drop ratios due to the reference crack should be generated for all considered modes of vibration.
For this aim, a crack having reference depth ratio, 0.25, is modelled and it is shifted along the beam by taking proper spacing between two neighbour sample locations.Minimum normalised distances, which accurately reveal the shapes of natural frequency drop functions, are determined as 0.1 for the first and second vibration modes, 0.05 for the third and fourth vibration modes, and 0.025 for the fifth and sixth vibration modes of the cantilever beam model.Curve fitting procedure is applied through reference frequency drop ratios found at the sample crack locations.As a result, reference polynomial functions are generated for all vibration modes.The first mode coefficients ( ) 1 ( k κ ) of the fractal natural frequency drop polynomial are obtained as follows: Here, th 4 order polynomial function is adequate to define the change of natural frequency.As the vibration mode considered increases, shift of the crack along the beam result in more frequent modulation of the fractal natural frequency drops.Therefore, more sample locations and polynomial terms are required to fit a curve on drop ratios of higher natural frequencies.Following coefficients of the reference functions with the order of 8, 12, 16, 20, and 24 are employed to identify the natural frequency drop ratios obtained for the second, third, fourth, fifth, and sixth vibration modes respectively.Structures 12 (2015) 2460-2479 Figure 4: Natural frequency drop ratios obtained by the matrix solution for the crack depth ratios 0.1 (◊), 0.25 (o), 0.4 (□), and the drop ratios obtained by Eq. ( 16) for the crack depth ratios 0.1 (--), 0.25 (-), 0.4 (---).The ratios are illustrated on the (a) first, (b) second, (c) third, (d) fourth, (e) fifth, and (f) sixth mode of vibrations.

Latin American Journal of Solids and
Polynomial functions, which have coefficients proportional to the reference coefficients, are defined for unknown crack depths by using the expression given in Eq. ( 18).Thus, the equation set, given in Eq. ( 21), is generated to find unknown crack parameters.
The cracked beam cases given in Table 2 are also considered as the scenarios of the inverse analysis.The results obtained by the theoretical and numerical natural frequency inputs are represented in Table 5. Theoretical natural frequency inputs tabulated by Table 3 are employed for checking the amount of errors caused by the simplifications given in Eq. ( 15) and Eq. ( 16).Results show that the cracks are revealed with some minor deviations even if the beam model includes three cracks which are inspected by the natural frequency knowledge of six vibration modes.It is seen that the deviations generally rise as the number of cracks increases.Similar results are also obtained when the same cracked beam cases are handled by using the finite element inputs.As one would expect, the deviations generally increase due to the numerical errors originated from the finite element calculations.It should be noted that increasing number and depth of cracks brings some difficulties arise from the growth of the errors in linear modelling.The errors in natural frequency drop ratios obtained by simplified model may still be negligible in respect of natural frequency determination.However, large errors, which stem from the excessive number and depth ratio of cracks, lead deviations in polynomial coefficients which make difficult to converge towards accurate crack parameters.Besides, excessive number of cracks needs the use of higher modes of vibration frequencies which frequently modulate due to the crack shifted along the beam.Rise of the number of modulations in fractal natural frequency drops makes the method convergence difficult and necessitates suitable initial predictions for the crack parameters.Based upon these effects, the inverse method can be suggested to be efficiently applied up to three cracks with depth ratios less than 0.5.
It is also worth noting that, when three cracks are detected using six polynomial functions, it is not judged about the possibility of existence of more cracks.In other words, six polynomial functions are able to detect three cracks, while giving definite knowledge up to two cracks.If there are two cracks that are inspected using six polynomial functions, methodology detects the third crack at random position with negligibly small depth ratio.From a different point of view, the use of larger number of polynomial functions is helpful in assessment of cracks near the free end.Because, especially in measurement, lower mode natural frequencies may have insufficient sensitivity to the effects of a crack near the free end unless the crack has significant depth ratio.As recognised from Figure 4, changes of the first two natural frequencies may be inadequate to distinguish the effects of the crack having normalised distance less than 0.1 from the free end.This inadequacy can be diminished to 0.05 and 0.025 by considering the first four and six natural frequencies respectively.Note that all of these limitations and thus success of the method can vary with the sensitivity and accuracy of measurements.Structures 12 (2015) 2460-2479 4.4 Testing the method on data from literature Simplified theories and inverse method used are also validated by using several tests given in current literature.At first, selected experiments of single cracked free-free beam scenarios reported by Kim and Stubbs (2003) Mazanoglu and Sabuncu (2012), which is also specified in Section 4.1, is employed for testing the methods.Finally, cantilever beam experimented by Ruotolo and Surace (1997)  For all scenarios, actual and detected parameters of the cracks are tabulated in Table 7. Results show that the cracks are approximately identified for all cases.Applications using free-free and cantilever beam conditions demonstrate that boundary changes do not clearly affect the success of the methodology if the polynomial function is prepared for each different boundary condition.It is also shown that the methodology works well even if the crack is quite close to the free end as in the first case and its depth is advanced as in the third case.The largest error ratios are obtained as 1.6% for the location parameter of second and sixth scenarios and 4.6% for the depth ratio parameter of the eighth scenario.It can be observed that most of the errors are smaller than the errors given in the references.

CONCLUSIONS
This paper presents forward and inverse approaches to determine natural frequencies of cracked beams and to identify cracks.In forward analyses, two theoretical simplifications are introduced to find natural frequency drop ratios caused by cracks.One of them is helpful to compute natural frequency ratios of multiple cracked / un-cracked beams.The other one sets relation between the natural frequency drop ratios for various crack depths.Results are theoretically validated using well-known analytical approach.Proposed theoretical simplifications are also engaged in novel methodology of crack identification.It is demonstrated here that researchers do not have to use complex theoretical formulations to express crack effects and to detect cracks with acceptable proximity.In presented methodology, a set of polynomials, which have coefficients varying with unknown crack parameters, is produced using reference drop ratios obtained along the beam.Each crack to be detected needs two polynomial functions corresponding to two different natural frequency drop ratios.The function set is solved using the numerical optimisation to get singular values exhibiting the crack parameters.Proposed crack detection methodology is verified by numerically obtained frequency ratios and the experimental frequency changes given in literature.It is seen that locations and depth ratios of the cracks are identified with admissible deviations.The methodology is powerful in the detection of one and two cracks.Furthermore, it can efficiently be employed up to three cracks if the initial estimations done for the crack parameters are not far from the accurate parameters.Limitation on crack number mainly stems from the growth of nonlinearity and rising number of modulations in fractal natural frequency drops.Here, it must be borne in mind that the identification of single and double cracked beams satisfies the people in most applications of non-destructive testing.Structures 12 (2015) 2460-2479 Several achievements can be remarked when the methodology is evaluated in respect of nondestructive evaluation.First of all, it only needs a dynamic data acquisition system for determining natural vibration frequencies.A point to consider is requirement of the un-cracked beam's natural frequencies which are used as base data for zero setting procedure.This provides the robustness of method that is indispensable for the automated crack detection system.Simplified approaches for determining natural frequencies reduce the computational work load by eliminating integrations and differentiations.Simple relations also make the inverse method available and increase the adaptation speed of the method applied to the beams having different properties.In addition, if there is no measurement fail, the method does not depend on to the user experiences.As a result, it can be concluded that present methodology can be proposed due to its inexpensiveness, practicability, robustness, user independence, and convenience for automation.

4. 1
Properties of the Beam Tested NumericallyUniform aluminium alloy cantilever beam is simulated for checking the efficiency of the methods presented.Dimensions of the beam are given as The modulus of elasticity, density, and Poisson ratio of the beam are properties of a cantilever beam model with two cracks are shown in Figure2.As a result of the theoretical computations, un-cracked beam natural frequencies are found as .Several examples of cracked beam cases are considered for the theoretical verifications given in following subsection.

Figure 2 :
Figure 2: Geometric properties of the cantilever beam with cracks.

Figure 3 :
Figure 3: Meshing view of a cracked region of the beam modelled in the finite element package.
are utilised.The test beam has the material properties as 206 Secondly, double crack scenarios of the aluminium alloy cantilever beam experimented by

Table 4 :
Ratios between the natural frequencies obtained by analytical large size matrix solution and present approximate approach.

Table 5 :
The cracks detected by present methodology using frequency inputs obtained by theoretical and finite element analyses.

Table 6 :
is considered with double crack scenarios.Their tests were conducted on C30 steel beam with the geometric properties as Measured natural frequency ratios inputted to the methodology are given in Table6.Measured natural frequency ratios used as inputs of the methodology.

Table 7 :
Cracked beam scenarios tested by several researchers and the cracks detected by present methodology.