#### Abstract

This paper presents a new sectional flexibility factor to simulate the reduction of the stiffness of a single-edge open cracked beam. The structural model for crack of the beam is considered as a rotational spring which is related to the ratio of crack depth to the beam height, *a/h*. The mathematical model of this single-edge open crack beam is considered as an Euler-Bernoulli beam. The modified factor, *f(a/h)*, derived in this paper is in good agreement with previous researchers' results for crack depth ratio *a/h* less than 0.5. The natural frequencies and corresponding mode shapes for lateral vibration with different types of single-edge open crack beams can then be evaluated by applying this modified factor *f(a/h)*. Using the compatibility conditions on the crack and the analytical transfer matrix method, the numerical solutions for natural frequencies of the cracked beam are obtained. The natural frequencies and the mode shapes with crack at different locations are obtained and compared with the latest research literature. The numerical results of the proposed cracked beam model obtained by this method can be extended to construct frequency contour. The natural frequencies measured from field can be used in solving the inverse problem to identify cracks in structures.

#### 1. Introduction

A crack in the structure will reduce the structural strength and result in severe damage under critical loading conditions. The major issue of the structural health monitoring is to detect crack depth and location in the present study. The model with linear elastic fracture mechanics and Euler-Bernoulli beam theory are being widely used in the recent research literatures. The cracked beam is modeled as two-segment beam with the crack simulated as a rotational spring. The crack of the beam is considered as a local flexibility which is a function of the crack depth.

Sih [1] proposed strain energy density factor theory to discuss the all mixed-mode crack extension problems. The strain energy density factor is a linear elasticity function of the mixed-mode stress intensity factors. Tada et al. [2] presented the stress intensity factors of different modes due to general loading. The presented model introduced the local flexibility matrix from the stress intensity factors. Nobile [3] proposed a simple method for obtaining approximate stress intensity of straight cracked beams. The system takes account of the elastic crack tip stress singularity while using the elementary beam theory. Dimarogonas and Paipetis [4] introduced a general stiffness matrix for cracked structural members to model the respective dynamic system. The local flexibility can be derived further from the general stiffness matrix. Chondros et al. [5] developed a continuous cracked beam vibration theory for the lateral vibration of cracked Euler-Bernoulli beams with single-edge or double-edge open cracks. The crack was modeled as a continuous flexibility using the displacement field. Anifantis and Dimarogonas [6] studied the system stability of the cracked column with vertical load. The method developed a general flexibility matrix to express the local flexibility of a beam with a single-edge crack. Ostachowicz and Krawczuk [7] developed a new local flexibility which was derived from the stress intensity factor by Anifantis and Dimarogonas [6].

In order to obtain the natural frequencies of the crack beam, finite element method was used to compute the eigensolutions in the recent literatures. The order of the determinants increases as the degree of freedom increases in finite element method. In order to reduce the order of the determinants, Lin et al. [8] proposed using transfer matrix for beams with arbitrary number of cracks. The method uses only four unknown constants which can be solved through satisfying four boundary conditions. Lin and Chang [9] used the analytic transfer matrix method to solve eigensolutions of a cracked cantilever beam. The eigenfunctions obtained in this method are analytical solutions. The dynamic responses can be obtained by this method, and the solutions converge quite fast. Alsabbagh et al. [10] presented a new simplified formula for the stress correction factor by using strain energy density approach. A modified factor for local flexibility was used in solving the characteristic equation of the cracked beam. Lin [11] used the Timoshenko beam theory and transfer matrix method to solve the direct and inverse problems of simply supported beam with a single-edge open crack. The location and crack size of the beam can be determined by the method presented. The theoretical results are also validated by a comparison with experimental measurements.

#### 2. Derivation of Stress Intensity Factor

A prismatic beam is considered with an open and nonpropagating crack of depth , length , height , and width . The singular stress distribution at the crack tip takes the form [3, 4, 10] with the conditions that acts at a distance from the tip and is the stress intensity factor. The normal stress acting on the reduced cross-section passing through the crack tip is given as where is the bending moment and is the moment of inertia of the remaining part of the cracked beam. The distance in (2.2) is found from Figure 1 as where is the distance from the neutral axis of the reduced cross-section to the tip as shown in Figure 1.

The stress condition is considered as at the crack tip . Substituting (2.3) into (2.2), and can be expressed as The distance can be determined from the equilibrium condition of forces along the -axis: The left-hand side of (2.6) is evaluated, using (2.5), to be The right-hand side of (2.6) is evaluated, using (2.4), to be Substitution of (2.7) and (2.8) into (2.6) leads to So Substituting (2.10) into (2.5), the stress intensity factor can be expressed as where

#### 3. Calculation of the Equivalent Flexibility in the Crack Beam

Let be the strain energy due to the crack. According to Castigliano’s theorem, the additional displacement is under general loading . In this work, the displacement will reduce to where the displacement is taken as the rotation since the bending moment is the only load of the structure.

The strain energy has the form [4–6, 10] where is the strain energy density function. Therefore, The moment substitutes for the generalized load The flexibility influence coefficient will be written as The strain energy density function has the form where for plane strain and and are Young’s modulus and Poisson’s ratio, respectively. The flexibility scalar is where Poisson’s ratio is taken as and the area moment of inertia is taken for the whole cross-section as . The nondimensional cracked section flexibility can be found from (3.7) as

#### 4. Free Vibration of a Cracked Beam

A simple beam with length and an open-edge crack at position is considered as shown in Figure 2. Euler-Bernoulli beam bending theory was used in solving the free vibration problem. According to [8, 9], the differential equation of motion for each segment is where is the density of the material and is the cross-section area of the rectangular beam. The boundary conditions of the simply supported beam are The beam at the crack position is simulated as a rational spring with sectional flexibility as shown in Figure 3. The continuous conditions at the crack position are and the compatibility condition due to the rational flexibility is From the above equations, the following quantities are introduced for nondimensional analysis: Equation (4.1) can then be expressed in a nondimensional form as The continuous conditions at the crack position as non-dimensional form are where is the non-dimensional cracked section flexibility as (3.8).

#### 5. Calculation of Natural Frequencies

The equation of motion for the cracked simple beam system can be expressed as (4.8). Using the method of separation of variables, , in (4.8), the differential equation for free vibration can be written as where From (4.9), the continuous conditions at the crack position are A closed-form solution to this eigenvalue problem can be obtained by employing transfer matrix methods [8, 9]. The general solution of (5.1), for each segment, is where , , , and are constants associated with the th segment . These constants of the second segment (, , , and ) are related to those of the first segment (, , , and ) through the continuous conditions in (5.3) and can be expressed as where is a transfer matrix which depends on eigenvalue and the elements are derived from [8].

Using (5.5), the four constants of the first segment (, , , and ) can be mapped into those of the second segment (, , , and ); thereby, the number of independent constants can be reduced to four. For the case of a simply supported beam, the corresponding boundary conditions of (4.2) and (4.3) can be written as Due to (5.6) and (5.8), (5.4) yields Satisfying the boundary conditions (5.7) and (5.9), (5.4) leads to the following equations: which can be expressed in matrix form as where Substituting (5.5) into (5.12) and applying (5.10), one obtains where A nontrivial solution of the simply supported beam requires The characteristic equation of the cracked simply supported beam can be obtained as where is the eigenvalues of the system. This characteristic equation can be solved by using the Newton-Raphson method to obtain the eigenvalues and corresponding eigenfunctions.

#### 6. Numerical Results

In order to verify the procedure presented in this paper, results obtained by applying this method are compared with the available data for single-edge open cracked beam. A 300 mm simple supported beam of cross-section , with modulus of elasticity , the density , the crack is located at the position and crack depth mm.

A simplified stress intensity factor in (2.11) is expressed in this paper. The function in (2.12) can be compared with the expression given by [2] Another stress intensity factor is [12]

In this research, small crack depth ratio which implies early stage of structure damage is considered. The comparison of data obtained by the proposed method and previous research is shown in Figure 4. The stress intensity factor of the present model and those in [2, 12] are rather close to each other for small crack depth ratio .

The nondimensional cracked section flexibility in (3.8) can be compared with the expression given by [5] Another non-dimensional cracked section flexibility is [7]

Three generalized loading conditions, bending, tension, and torsion, on a cracked beam were considered to evaluate sectional flexibility in the past literature. Beams are mainly affected by bending moment in most loading cases; therefore only bending effects are considered in evaluating the simplified cracked section flexibility. The results of the proposed method of non-dimensional cracked section flexibility are compared with those of previous research in Figure 5. The results of this simplified method and those of [5, 7] are in good agreement for small crack depth ratio . A cracked beam with greater than 0.5 is already severely damaged and is not suitable for applying this sectional flexibility model.

The first five natural frequencies of the uncracked beam are calculated as , , , , and Hz. The ratios of natural frequencies between cracked beam and uncracked beam are listed in Table 1. The variation of natural frequency ratio with the crack depth of this simply supported beam with a crack located at the section () for the first five modes is plotted as in Figure 6.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

With a crack located at the center of the beam, the ratios of natural frequencies between cracked beam and uncracked beam are listed in Table 2. The variation of natural frequency ratio with the crack depth of a simply supported beam with a crack located at () is shown in Figure 7 for the first five modes.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

It is quite obvious that the natural frequencies decrease due to the existence of cracks. That is due to the cracked beam becoming more flexible due to the reduction of moment of inertia of the section property.

Figures 8(a)–8(e) show the first five mode shapes of a cracked simply supported beam with single open crack at , crack ratio , and normalized amplitude . It is obvious that the mode shapes all show turning point at crack location .

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

#### 7. Conclusion

The simplified stress intensity factor and flexibility were derived utilizing the crack beam theorem of Nobile [3] and Dimarogonas [4, 5]. The order of polynomial functions for crack depth ratio is reduced, because, with a crack depth ratio , the high-order terms will approach to 0. Therefore, these higher-order arithmetic terms can be neglected and, in order to predict the early stage of structural damage, crack depth ratio should be less than 0.5. The simplified stress intensity factor and flexibility are compared with the recent literature, and the numerical results are found to be in good agreement.

If the crack is right on the position of nodal point of certain modes, frequencies ratio shows no difference with . For example, the crack at 0.8 *L* is also a nodal point of the fifth mode which gives the numerical result as . The numerical results obtained by this method are in good agreement with the actual vibration response of a simply supported beam. The turning points of certain mode shape function reveal the information about crack location. For the case with a crack at 0.8 *L*, it can be seen obviously from shape function of mode two, mode three, and mode four with a turning point at 0.8 *L*.

The simplified stress intensity factor and flexibility of this method can be further extended to construct frequencies ratio contours for beams with cracks. The natural frequencies obtained by applying this model can be used to verify the experimental measurements in a similar way to that in [11]. The location and crack depth of a beam can then be identified as an inverse problem by matching up field measurement of frequencies of a cracked beam.