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 a, length L, height h, and width t. The singular stress distribution at the crack tip takes the form [3, 4, 10]σxs=KI2πr,
with the conditions that σxs acts at a distance r=b from the tip and KI is the stress intensity factor. The normal stress acting on the reduced cross-section passing through the crack tip is given as
σx=MIy,
where M is the bending moment and I=t(h-a)3/12 is the moment of inertia of the remaining part of the cracked beam. The distance y in (2.2) is found from Figure 1 asy=y̅-b=h2-a2-b,
where y̅ is the distance from the neutral axis of the reduced cross-section to the tip as shown in Figure 1.
Geometry of a beam with a single-edge crack.
The stress condition is considered as σxs=σx at the crack tip r=b. Substituting (2.3) into (2.2), σx and KI can be expressed asσx=12Mt(h-a)3(h2-a2-b),KI=2πbσx=12M2πbt(h-a)3(h2-a2-b).
The distance b can be determined from the equilibrium condition of forces along the x-axis:∫0bKI2πrdr=∫y̅-by̅σxdy.
The left-hand side of (2.6) is evaluated, using (2.5), to be∫0bKI2πrdr=24Mbt(h-a)3(h2-a2-b).
The right-hand side of (2.6) is evaluated, using (2.4), to be∫y̅-by̅σxdy=6Mbt(h-a)3(h-a-b).
Substitution of (2.7) and (2.8) into (2.6) leads to24Mbt(h-a)3(h2-a2-b)=6Mbt(h-a)3(h-a-b).
Sob=13(h-a).
Substituting (2.10) into (2.5), the stress intensity factor can be expressed asKI=6Mth3/2F(ah),
whereF(ah)=2π27⋅1(1-a/h)3=0.482(1-a/h)3.
3. Calculation of the Equivalent Flexibility in the Crack Beam
Let UT be the strain energy due to the crack. According to Castigliano’s theorem, the additional displacement is ui=∂UT/∂Pi under general loading Pi. In this work, the displacement will reduce toθ=∂UT∂M,
where the displacement ui is taken as the rotation θ since the bending moment M is the only load of the structure.
The strain energy has the form [4–6, 10]UT=∫0a∂UT∂ada=t∫0aJda,
where J is the strain energy density function. Therefore,ui=∂UT∂Pi=∂∂Pi[t∫0aJ(a)da](Paris'sequation).
The moment M substitutes for the generalized load Piθ=∂∂M[t∫0aJ(a)da].
The flexibility influence coefficient will be written asc=∂θ∂M=∂2∂M2[t∫0aJ(a)da].
The strain energy density function J has the formJ=KI2E′,
where E′=E/(1-υ2) for plane strain and E and υ are Young’s modulus and Poisson’s ratio, respectively. The flexibility scalar isc=∂2∂M2[t∫0aKI2E′da]=2π(1-υ2)h[1-(1-a/h)2]9EI(1-a/h)2,
where Poisson’s ratio is taken as υ=0.3 and the area moment of inertia is taken for the whole cross-section as I=th3/12. The nondimensional cracked section flexibility can be found from (3.7) asc*=EIcL=2π(1-υ2)[1-(1-a/h)2]9(1-a/h)2⋅hL.
4. Free Vibration of a Cracked Beam
A simple beam with length L and an open-edge crack at position X1 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 isEI∂4Yi(X,T)∂X4+ρA∂2Yi(X,T)∂T2=0Xi-1<X<Xi,i=1,2,
where ρ is the density of the material and A is the cross-section area of the rectangular beam. The boundary conditions of the simply supported beam areY(0,T)=Y(L,T)=0,Y′′(0,T)=Y′′(L,T)=0.
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 areY1(X1-,T)=Y2(X1+,T),Y1′′(X1-,T)=Y2′′(X1+,T),Y1′′′(X1-,T)=Y2′′′(X1+,T),
and the compatibility condition due to the rational flexibility isY2′(X1+,T)-Y1′(X1-,T)=EIcY2′′(X1+,T).
From the above equations, the following quantities are introduced for nondimensional analysis:y=YL,x=XL,t=T,xi=XiL,l1=L1L,l2=L2L.
Equation (4.1) can then be expressed in a nondimensional form asEIL4∂4yi(x,t)∂x4+ρA∂2yi(x,t)∂t2=0xi-1<x<xi,i=1,2.
The continuous conditions at the crack position as non-dimensional form arey1(x1-,t)=y2(x1+,t),y1′′(x1-,t)=y2′′(x1+,t),y1′′′(x1-,t)=y2′′′(x1+,t),y2′(x1+,t)-y1′(x1-,t)=c*y2′′(x1+,t),
where c* is the non-dimensional cracked section flexibility as (3.8).
Schematic diagram of a simply supported cracked beam.
Model for the cracked beam with sectional flexibility.
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, yi(x,t)=wi(x)ejωt, in (4.8), the differential equation for free vibration can be written aswi′′′′(x)-λ4wi(x)=0,xi-1<x<xi,i=1,2,
whereλ4=ρAω2L4EI.
From (4.9), the continuous conditions at the crack position arew1(x1-)=w2(x1+),w1′′(x1-)=w2′′(x1+),w1′′′(x1-)=w2′′′(x1+),w2′(x1+)-w1′(x1-)=c*w2′′(x1+).
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, iswi(x)=Aisinλ(x-xi-1)+Bicosλ(x-xi-1)+Cisinhλ(x-xi-1)+Dicoshλ(x-xi-1),xi-1<x<xi,i=1,2,
where Ai, Bi, Ci, and Di are constants associated with the ith segment (i=1,2). These constants of the second segment (A2, B2, C2, and D2) are related to those of the first segment (A1, B1, C1, and D1) through the continuous conditions in (5.3) and can be expressed as{A2B2C2D2}=[t11t12t13t14t21t22t23t24t31t32t33t34t41t42t43t44]{A1B1C1D1}=T̲4×4{A1B1C1D1},
where T̲4×4 is a 4×4 transfer matrix which depends on eigenvalue λ and the elements are derived from [8].
Using (5.5), the four constants of the first segment (A1, B1, C1, and D1) can be mapped into those of the second segment (A2, B2, C2, and D2); 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 asY(0,T)=0⟶w(0)=0,Y(L,T)=0⟶w(1)=0,Y′′(0,T)=0⟶w′′(0)=0,Y′′(L,T)=0⟶w′′(1)=0.
Due to (5.6) and (5.8), (5.4) yieldsB1=0,D1=0.
Satisfying the boundary conditions (5.7) and (5.9), (5.4) leads to the following equations:A2sinλl2+B2cosλl2+C2sinhλl2+D2coshλl2=0,-A2sinλl2-B2cosλl2+C2sinhλl2+D2coshλl2=0,
which can be expressed in matrix form as{00}=[sinλl2cosλl2sinhλl2coshλl2-sinλl2-cosλl2sinhλl2coshλl2]{A2B2C2D2}=B̲2×4{A2B2C2D2},
whereB̲2×4=[sinλl2cosλl2sinhλl2coshλl2-sinλl2-cosλl2sinhλl2coshλl2].
Substituting (5.5) into (5.12) and applying (5.10), one obtains{00}=B̲2×4{A2B2C2D2}=B̲2×4⋅T̲4×4{A1B1C1D1}=R̲2×4{A1B1C1D1},
whereR̲2×4=B̲2×4⋅T̲4×4=[r11r12r13r14r21r22r23r24].
A nontrivial solution of the simply supported beam requiresdet[r11(λ)r13(λ)r21(λ)r23(λ)]=0.
The characteristic equation of the cracked simply supported beam can be obtained asc*λnsinh(λnl1)sinh(λnl2)sin(λn(l1+l2))-c*λnsin(λnl1)sin(λnl2)sinh(λn(l1+l2))+2sinh(λn(l1+l2))sin(λn(l1+l2))=0,
where λn 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 20×20mm2, with modulus of elasticity E=2.06×1011N/m2, the density ρ=7800kg/m3, the crack is located at the position 240mm and crack depth a=8 mm.
A simplified stress intensity factor KI in (2.11) is expressed in this paper. The function F(a/h) in (2.12) can be compared with the expression given by [2]KI=6Mth3/2F2(ah),F2(ah)=πah2hπatan(πa2h)0.923+0.199(1-sin(πa/2h))4cos(πa/2h).
Another stress intensity factor is [12]KI=6Mth3/2F3(ah),F3(ah)=πah[1.122-1.44(ah)+7.33(ah)2-13.08(ah)3+14(ah)4].
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 a/h.
Stress intensity factor variations.
The nondimensional cracked section flexibility c* in (3.8) can be compared with the expression given by [5]c2=6π(1-υ2)⋅hL⋅Φ(ah),Φ(ah)=0.6272(ah)2-1.04533(ah)3+4.5948(ah)4-9.9736(ah)5+20.2948(ah)6-33.0351(ah)7+47.1063(ah)8-40.7556(ah)9+19.6(ah)10.
Another non-dimensional cracked section flexibility is [7]c3=6π(ah)2⋅hL⋅fJ(ah),fJ(ah)=0.6384-1.035(ah)+3.7201(ah)2-5.1773(ah)3+7.553(ah)4-7.332(ah)5+2.4909(ah)6.
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/h. A cracked beam with a/h 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 ω01=517.85, ω02=2071.4, ω03=4660.64, ω04=8285.58, and ω05=12946.22Hz. 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 (L1/L=0.8) for the first five modes is plotted as in Figure 6.
Natural frequency ratio with the crack depth at location 0.8.
Crack location
Crack depth
Natural frequency ratio
L1/L
a/h
ω1/ω01
ω2/ω02
ω3/ω03
ω4/ω04
ω5/ω05
0.8
0.0
1.0000
1.0000
0.9999
1.0000
0.9999
0.8
0.1
0.9966
0.9912
0.9913
0.9967
0.9999
0.8
0.2
0.9918
0.9792
0.9803
0.9927
0.9999
0.8
0.3
0.9850
0.9628
0.9660
0.9876
0.9999
0.8
0.4
0.9747
0.9394
0.9474
0.9814
0.9999
0.8
0.5
0.9581
0.9052
0.9237
0.9737
0.9999
0.8
0.6
0.9292
0.8544
0.8943
0.9648
0.9999
0.8
0.7
0.8735
0.7800
0.8605
0.9551
0.9999
0.8
0.8
0.7525
0.6809
0.8271
0.9459
0.9999
0.8
0.9
0.4791
0.5825
0.8021
0.9392
0.9999
Variation of natural frequency ratio with the crack depth of a simply supported beam (L1/L=0.8): (a) first frequency, (b) second frequency, (c) third frequency, (d) fourth frequency, and (e) fifth frequency.
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 (L1/L=0.5) is shown in Figure 7 for the first five modes.
Natural frequency ratio with the crack depth at location 0.5.
Crack location
Crack depth
Natural frequency ratio
L1/L
a/h
ω1/ω01
ω2/ω02
ω3/ω03
ω4/ω04
ω5/ω05
0.5
0.0
1.0000
1.0000
0.9999
1.0000
0.9999
0.5
0.1
0.9902
1.0000
0.9904
1.0000
0.9905
0.5
0.2
0.9770
1.0000
0.9778
1.0000
0.9786
0.5
0.3
0.9586
1.0000
0.9613
1.0000
0.9636
0.5
0.4
0.9321
1.0000
0.9391
1.0000
0.9448
0.5
0.5
0.8927
1.0000
0.9094
1.0000
0.9218
0.5
0.6
0.8314
1.0000
0.8699
1.0000
0.8948
0.5
0.7
0.7328
1.0000
0.8199
1.0000
0.8654
0.5
0.8
0.5723
1.0000
0.7634
1.0000
0.8377
0.5
0.9
0.3245
1.0000
0.7144
1.0000
0.8175
Variation of natural frequency ratio with the crack depth of a simply supported beam (L1/L=0.5): (a) first frequency, (b) second frequency, (c) third frequency, (d) fourth frequency, and (e) fifth frequency.
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 L1/L=0.8, crack ratio a/h=0.4, and normalized amplitude Y/Ymax. It is obvious that the mode shapes all show turning point at crack location L1/L=0.8.
(a) Normalized mode shape of the first mode with a crack location L1/L=0.8 and crack depth ratio a/h=0.4. (b) Normalized mode shape of the second mode with a crack location L1/L=0.8 and crack depth ratio a/h=0.4. (c) Normalized mode shape of the third mode with a crack location L1/L=0.8and crack depth ratio a/h = 0.4. (d) Normalized mode shape of the fourth mode with a crack location L1/L=0.8 and crack depth ratio a/h=0.4. (e) Normalized mode shape of the fifth mode with a crack location L1/L=0.8 and crack depth ratio a/h=0.4.
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 a/h is reduced, because, with a crack depth ratio a/h<0.5, 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 a/h 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 ωi/ω0i=1. For example, the crack at 0.8 L is also a nodal point of the fifth mode which gives the numerical result as ω5/ω05=1. 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.
SihG. C.1973Leiden, The NetherlandsNoordhoff International0449095ZBL0319.73053TadaH.ParisP. C.IrwinG. R.19852ndSt. Louis, Mo, USAParis ProductionsZBL0598.30060NobileL.lucio.nobile@mail.ing.unibo.itMixed mode crack initiation and direction in beams with edge crack200033210711610.1016/S0167-8442(00)00006-9DimarogonasA. D.PaipetisS. A.1983New York, NY, USAApplied ScienceChondrosT. G.DimarogonasA. D.YaoJ.A continuous cracked beam vibration theory19982151173410.1006/jsvi.1998.1640AnifantisN.DimarogonasA.Stability of columns with a single crack subjected to follower and vertical loads198319428129110.1016/0020-7683(83)90027-6ZBL0523.73042OstachowiczW. M.KrawczukM.Analysis of the effect of cracks on the natural frequencies of a cantilever beam1991150219120110.1016/0022-460X(91)90615-QLinH. P.ChangS. C.WuJ. D.Beam vibrations with an arbitrary number of cracks2002258598799910.1006/jsvi.2002.5184LinH. P.linhp@mail.dyu.edu.twChangS. C.Forced responses of cracked cantilever beams subjected to a concentrated moving load200648121456146310.1016/j.ijmecsci.2006.06.014ZBL1192.74182AlsabbaghA. S. Y.AbuzeidO. M.oabuzeid@ju.edu.joDadoM. H.Simplified stress correction factor to study the dynamic behavior of a cracked beam200933112713910.1016/j.apm.2007.10.018ZBL1167.74572LinH. P.linhp@mail.dyu.edu.twDirect and inverse methods on free vibration analysis of simply supported beams with a crack200426442743610.1016/j.engstruct.2003.10.014BrownW. F.SrawleyJ. E.1966American Society Testing and Materials