General Article

International Journal of Sustainable Building Technology and Urban Development. 31 December 2020. 209-221
https://doi.org/10.22712/susb.20200016

ABSTRACT


MAIN

  • Introduction

  • Problem Formulation

  • Determination of the equivalent axial spring stiffness Kx

  • Frequency equation of a cracked cantilever beam

  • Results interpretation

  •   Numerical example

  • Conclusion

Introduction

The presence of cracks or damage in structural elements alters its dynamic response. For the structural integrity, it is very important to study changes in the modal parameters due to the presence of cracks inside the structure. Existing damage in a beam or a bar can cause a stiffness reduction, and change the fundamental frequency of the beam. Because this damage can constitute an important toughening mechanism, problems dealing with dynamics behavior of cracked beams have received considerable research attention since they were introduced to fracture mechanics. As a result, a wide body of literature on this topic exists (Chabaat et al.2005, 2011) [1, 2]. In the last decade, Morassi (1993) [3]; Dimarogonas (1996) [4]; Cacciola et al. (2003) [5] and Yang et al. (2001) [6] used the energy method to study the effect of an existing crack on the dynamic behavior of a beam during the vibration. In their analysis, a numerical model is developed and the equivalent bending stiffness is computed. On the other hand, Metveev and Bovsunovsky(2002) [7] investigated the variation of the natural frequency of an Euler-Bernoulli cracked beam. Once they got expressions for bending vibrations, they have studied the effects of the ratio of the crack location to the length of the beam as well as the depth ratio of the crack to the height of the beam. Thus, effects of the crack location on the fundamental natural frequency have been discussed on a number of papers by researchers such as (Chondros et al. 1997) [8] and Collins (1992) [9]. The most common crack models in the literature to simulate the vibratory behavior of cracked beams are essentially divided into three main categories; a local stiffness reduction, a continuous model and a spring stiffness model. The first one is based on the local reduction of stiffness; it simply requires subdividing the beam with a sufficient number of elements and reducing the rigidity at the position where the damage has occurred (Vestroni and Danilo, 2000) [10]. This model uses the finite element method in order to introduce the reduction of stiffness in the damaged element to associate the real effects due to the presence of the crack. This model is able to solve the inverse problem by means of the determination of the location and the size of cracks in a damaged structure. Nevertheless, limitations of this model reside in the fact that the stiffness of the cracked section depends on the refinement of the model. Besides, a continuous model is considered by many researchers such as Chondros et al. (2001) [11] that developed consistent continuous theory on a single cracked bar under vibrations. The fixed-free bar with an edge transverse crack was chosen and the perturbation solution were able to determine the changes in longitudinal vibrations of natural frequencies. In another research work, (Chondros et al. 1998) [12] employed the same model to determine the natural frequency of a cracked cantilever bar and the obtained results agreed well with those from the experimental tests.

According to Castigliano’s theorem, the spring stiffness model is largely employed to simulate the dynamic behavior of cracked beams. In a previous study, Dimarogonas and Paipetis (1983) [13] developed massless rotational spring models based on the fracture mechanics, meanwhile; Adams et al. (1978) [14] considered a free-free cracked bar where the crack was modeled by a massless linear spring. The natural frequencies of longitudinal vibration were obtained experimentally and analytically. On the others hand, Springer et al. (1988) [15] have analyzed the problem of a cracked bar under a free longitudinal vibration. The modeling of the crack has been done in two different ways: One used linear springs and the other considered reductions in the cross-sectional area. The changes in natural frequencies are closer to those obtained from experiments. Loya and Rubio (2006) [16] employed the rotational massless spring model to analyze the lateral vibration of a cracked Timoshenko’s beam. The natural frequencies of the cracked beam were obtained by solving the differential equations of motion. Sayyad et al. (2012, 2013) [17, 18] were able to use an equivalent axial spring, connecting two segments of the beam to simulate the crack. In their study, the crack size and location can easily be detected from the measured longitudinal vibration data. Recently, Ayas et al (2019) [19] introduced a new numerical approach for studying a cantilever bar having a transverse crack. In their research work, the crack is modeled by an elastic longitudinal spring with a stiffness K. The extent of the damage in cracked bars can easily be estimated by the measurement of the difference between the dynamic response of an uncracked bar and a cracked bar.

Problem Formulation

In this study, a longitudinal vibration of a cantilever bar with a single crack at a clamped end as shown in Figures. 1 and 2 is considered. The length of the bar is L and the constant cross-section area is A. Its material proprieties, Young modulus E and the mass density ρ are given in Table 1. The effect of material damping on the natural longitudinal frequencies is negligible.

https://static.apub.kr/journalsite/sites/durabi/2020-011-04/N0300110402/images/Figure_susb_11_04_02_F1.jpg
Figure 1.

The real cracked cantilever bar.

https://static.apub.kr/journalsite/sites/durabi/2020-011-04/N0300110402/images/Figure_susb_11_04_02_F2.jpg
Figure 2.

Geometry of cracked section of cantilever bar.

Table 1.

Material properties

Material Aluminum
Modulus of elasticity 70 Gpa
Density 2700 Kg/m
Poisson's ratio 0.35

Considering a one-dimensional continuum media, Doyle (1997) [20] gives the equation of motion governing the longitudinal vibration of a uniform bar;

(1)
2u(x,t)x2=1c22u(x,t)t2

The transverse edge crack is modeled by a translational axial spring with stiffness Kx as shown in Figure 3.

https://static.apub.kr/journalsite/sites/durabi/2020-011-04/N0300110402/images/Figure_susb_11_04_02_F3.jpg
Figure 3.

Model of the cantilever bar with a crack at the clamped end.

There are several expressions of spring stiffness that can be found in literature for the transverse vibration (Rizos et al. 1990) [21] and (Ostachowicz and Krawczuk, 1991) [22]. Ruotolo and Surace (2004) gave concerning the longitudinal vibration, the equivalent axial spring stiffness expression for double edge cracked bar [23] and by analogy to the bending vibration of an edge cracked beam, Sayyad et al (2013) [18] were the ones who developed expressions of axial spring for longitudinal vibration. In the next section based on the Castigliano’s theorem and the strain energy release rate for an edge-cracked bar, we determine the expression of axial spring.

Determination of the equivalent axial spring stiffness Kx

The longitudinal spring stiffness equivalent to the transverse edge crack is based on stress intensity factor and Castigliano’s theorem. Under a remote load, the additional displacement ui along the direction of a force Fi due to the presence of the crack is determined using Castiglione’s theorem along with Paris equation (Tada et al., 2000) [24] and (Sih, 1973) [25].

If UT is the strain energy due to the crack; Castiglione’s theorem requires that the additional displacement is;

(2)
ui=UTFi

Along the force Fi, the strain energy takes the following form (1) ;

(3)
UT=0aUTada=b0aJ(a)da

where b is the bar thickness and Js the strain energy density function.

(4)
J(a)=1-v2KI2E

and KI is the stress intensity factor for fracture mode I.

Due to the stress σ, SIF in mode I is given by;

(5)
KI=σπaIah
(6)
with,Iah=1.12-0.231ah+10.55ah2-21.71ah3+30.39ah4

σ denotes the remote stress, a is the crack depth, h is the height of the element bar.

(7)
b0aJ(a)da=(1-v2)Eσ2πa.Iah2da
(8)
=bh2(1-v2)Eσ2πa.fahda

The local flexibility can be computed as follows;

(9)
C=UTFi=(1-v2)2πhf(a/h)EA

The equivalent axial spring stiffness Kx due to the crack can also be determined by the following formula;

(10)
Kx=1C=EA(1-v2)2πhf(a/h)

where;

(11)
f(a/h)=0.627(a/h)2-0.172(a/h)3+5.921(a/h)4-10.700(a/h)+31.564(a/h)6-67.445(a/h)7+139.048(a/h)8-146.577(a/h)9+92.306(a/h)10

Frequency equation of a cracked cantilever beam

Rao (2007) [26] gives the solution of Eq. 1 for a free longitudinal vibration of a bar;

(12)
u(x,t)=A¯cosωxc+B¯sinωxc(Ccosωt+Dsinωt)

At the free end x = 0, we get;

(13)
EAu(o,t)x=0

Using both Eqs. 12 and 13 yields to;

(14)
B=0

Thus, Eq. 12 reduces to;

(15)
u(x,t)=Acosωxc(Ccosωt+Dsinωt)

The boundary condition at x = L can be expressed as;

(16)
EAu(L,t)x=-Kxu(L,t)

Substitution of Eqs.17 and 18 leads to the following expression;

(17)
-AωEAcsinωxc(Ccosωt+Dsinωt)=-KxAcosωLc(Ccosωt+Dsinωt)
(18)
EAωcKx=cotωLc

Introducing the mass of the bar as;

(19)
m=ρAL

Eq. 17 can be rewritten under the form;

(20)
λcotλ=β

where ; λ is the frequency parameter for both uncracked and cracked bar;

(21)
λ=ωLc
(22)
β=KKx=(1-v2)2πhLf(a/h)

By neglecting the higher order values in Eq. 11, Eq.22 becomes;

(23)
β=(1-v2)1.254πhLah2
(24)
λcotλ=(1-v2)1.254πhLah2

For a limited development of characteristics Eq. 24, λcotλ can be written;

(25)
λcotλ=1-λ23-λ445-2λ6945-o(8)Forλπ

In order to determine the influence of the higher order terms of Eq. 25 on the response frequency of a cracked cantilever bar, the equation is solved for considering 1, 2 and 3 higher order terms. Solutions for these last three cases are compared with the solution given by Eq. 24.

Case 1: One higher order term of Eq. 25;

Thus, Eq. 24 can be written as;

(26)
1-λ23=1-v21.254πhLah2
(27)
λ=31-1-v21.254πhLah2
(28)
ω=cL31-1-v21.254πhLah2

Case 2: Two higher order terms of Eq. 25:

Then, Eq. 24 can be written as;

(29)
1-λ23-λ445=β

After simplification Eq. 29 can be written as;

(30)
λ4+15λ2-45(1-β)=0
(31)
Let;y=λ2

Eq.30 becomes;

(32)
y2+15y-45(1-β)=0

The solution of this equation is given by;

(33)
y2+15y-45(1-β)=0
(34)
where;λ=3259-4β-5

Thus, the following equation is obtained;

(35)
ω=cL3259-41-v21.254πhLah2-5

Case 3: Three higher order terms of Eq. 25:

Ignoring higher order term of Eq. 25, the frequency Eq. 25 can be rewritten as follows;

(36)
1-λ23-λ445-2λ6945=β=(1-v2)1.254πhLah2

Let, y=λ2 then Eq. 36 becomes;

(37)
y3+10.5y2+157.5y-472.5(1-β)=0

Eq. 37 is a third-degree equation of the form:

(38)
y3+a1y2+a2y+a3=0

Solution to this equation is obtained as follows;

(39)
Q=3a2-a129=40.25
(40)
R=9a1a2-27a3-2a1354=469-236.25β

where;

(41)
S=R+Q3+R23
(42)
T=R-Q3+R23

Substitution of Eqs.39 and 40 into Eqs. 41 and 42, respectively, we obtained;

(43)
S=469-236.25β+(40.25)3+(469-236.25β)23
(44)
T=469-236.25β+(40.25)3+(469-236.25β)23

Solution of Eq.38 is given as follows;

(45)
y1=S+T-13a1

The frequency parameter λ can be expressed as;

(46)
λ=y1y10

The fundamental frequency can be computed by the following equation:

(47)
ω=cLy1

where;

(48)
y1=469-236.25β+(40.25)3+(469-236.25β)23+469-236.25β-(40.25)3+(469-236.25β)23-3.5

Values of frequency corresponding to the first mode for longitudinal vibration of an uncracked cantilever bar are summarized in Table 2.

Table 2.

Value of λ for uncracked bar (a/h=0)

Eq. 28Eq.35 Approximate Analytical solution Eq. 47 Numerical solution of Eq. 24
λ 1.732 1.601 1.578 1.571

It can be seen from Table 1 that the proposed approach can easily give the value of the frequency parameter λ for an uncracked bar (a/h = 0). The accuracy of the approximate solution depends mainly on the terms of higher order considered for calculation. It is noticed that the accuracy value of λ is obtained when one considers for the calculation three terms of higher order (Eq.47).

Results interpretation

Numerical example

A cracked cantilever bar under free longitudinal vibration is analyzed. The crack is located near the clamped end of the bar. The dimensions of the bar are 0.22 m x 0.0254 m x 0.01 m, the material properties are given in Table 2. The effect of the terms of higher order on the fundamental mode shape of vibration for both uncracked and cracked bars are studied (see Figures 4 and 5), while effects of crack depth ratio on the fundamental longitudinal frequency for a cracked cantilever bar are presented in Table 3 and shown in Figure 6.

https://static.apub.kr/journalsite/sites/durabi/2020-011-04/N0300110402/images/Figure_susb_11_04_02_F4.jpg
Figure 4.

First mode shape of uncracked cantilever bar for (β=0 and a/h=0).

https://static.apub.kr/journalsite/sites/durabi/2020-011-04/N0300110402/images/Figure_susb_11_04_02_F5.jpg
Figure 5.

First mode shape of a cracked cantilever bar for (β=0.0997 and a/h=0.5).

Table 3.

Values of First longitudinal frequency (Hz) for different crack depth ratio

a/h Β Approximate Analytical solution (Eq.47) Numerical solution of equation.24 Error %
0.0 0,0000 5923,517 5896,114 -0,46
0.1 0,0040 5912,255 5885,978 -0,45
0.2 0,0160 5882,225 5857,825 -0,42
0.3 0,0359 5833,425 5809,025 -0,42
0.4 0,0638 5762,103 5739,580 -0,39
0.5 0,0997 5668,257 5647,611 -0,37
0.6 0,1436 5548,135 5531,618 -0,30

https://static.apub.kr/journalsite/sites/durabi/2020-011-04/N0300110402/images/Figure_susb_11_04_02_F6.jpg
Figure 6.

Longitudinal natural frequency ratio vs. crackdepth ratio for cracked cantilever bar.

Effect of high terms of frequency on the mode shape for a cracked bar

In order to determine the influence the high order terms of the frequency parameter on the fundamental longitudinal mode shape; Figures 4 and 5 shows the modes shapes for an uncracked and a cracked cantilever bar obtained by the both approximate analytical solution (Eqs. 35 and 47) and numerical solution of analytical equation (Eq. 24). It is observed that the three curves are almost identical either in the case of the cracked bar (a/h = 0.5) or the case of the intact bar (a/h = 0).

Determination of the longitudinal frequency of a cantilever cracked bar

Longitudinal natural frequencies of the cracked bar for the first mode of vibration have been studied and are given in Table 3.

Figure 6 shows the fundamental normalized longitudinal frequencies (w/w0) determined by the approximate analytical solution (Eqs. 35 and Eq47) and also by the numerical solution (Eq. 24) for different crack depth ratios (a/h). It can be seen that the estimated values of the first longitudinal frequency of the cracked cantilever bar are very close to those calculated from the analytical solution. The crack is modeled by a longitudinal spring. According to Eq. (11), it is noticed that when the crack depth ratio increases, normalized longitudinal frequencies of a cracked cantilever bar (w/w0) decreases. In this case, the spring stiffness undergoes the same reduction and involves a reduction in the stiffness of the cracked bar. On the other hand, the proposed approach gives values of the first longitudinal frequency for different crack depth ratios which differ only by about -0.47% from those calculated from the numerical solution of the analytical solution.

Mode shapes analysis of a cracked cantilever bar

The modes shapes for longitudinal vibration of cracked and uncracked cantilever bar are presented and compared in Figures 7, 8, 9.

https://static.apub.kr/journalsite/sites/durabi/2020-011-04/N0300110402/images/Figure_susb_11_04_02_F7.jpg
Figure 7.

First mode of longitudinal vibration for cracked cantilever bar.

https://static.apub.kr/journalsite/sites/durabi/2020-011-04/N0300110402/images/Figure_susb_11_04_02_F8.jpg
Figure 8.

cond mode of longitudinal vibration for cracked cantilever bar.

https://static.apub.kr/journalsite/sites/durabi/2020-011-04/N0300110402/images/Figure_susb_11_04_02_F9.jpg
Figure 9.

Third mode of longitudinal vibration for cracked cantilever bar.

Figure 7 presents the influence of the crack depth ratio on the fundamental longitudinal mode shape for cracked cantilever bar which the frequency parameter λ is calculated by the approximate analytical approach. The results indicate that there is a deviation for the first mode shapes of the cantilever cracked bar compared with those of the uncracked bar. It is observed that this deviation increases with an increase of the crack depth ratio. For this case, different crack depth ratio (a/h = 0 to a/h = 0.6) are considered. For a crack depth ratio a/h < 0.3, small changes in mode shapes of a cracked cantilever bar are noticed and for a crack depth ratio (0.3 > a/h < 0.6), the longitudinal mode shape of a cracked cantilever bar undergoes changes.

Figures 8 and 9 show the second and the third mode shape of a cracked and an uncracked beam, respectively. It can be noticed from these figures the deviation between a cracked and an uncracked beam at the position of the crack increases with an increase in the mode of vibration. It is found that for all mode shapes, the difference between the mode shapes of a cracked and an uncracked beam is always sharply changed at the crack position with a depth ratio a/h greater than 0.3.

Conclusion

In this research work, the fundamental longitudinal natural frequencies and mode shapes of a cantilever bar with a crack near at clamped end are analyzed. A frequency equation corresponding to this particular problem is obtained and developed in series of terms, neglecting higher order terms. This equation is converted into a third degree equation. The solution for this last gives an expression of the fundamental longitudinal frequency of a cracked cantilever bar. The influence of different crack’s depths ratio on the normalized longitudinal natural frequencies and mode shapes of a cracked cantilever bar has been shown. It is found that the increase in the crack depth ratio involves a decrease in the normalized fundamental longitudinal frequency. It is also shown in this study that this approach for a cantilever cracked bar by the approximate analytical approach is suitable and useful in the determination of the fundamental longitudinal frequency. Obtained results are very close to those determined by the numerical solution from the analytical equation.

References

1
M. Chabaat, S. Djouder, and M. Touati, Semi Empirical Stress Analysis of a Brittle Material in a Vicinity of a Stress Concentrator. Applied Mechanics and Materials, Trans Tech Publications, 2005, pp. 243-252. 10.4028/www.scientific.net/AMM.3-4.243
2
M. Chabaat and H. Ayas, Complex variable Green's functions for crack-microcracks interactions. Key Engineering Materials, Trans Tech Publications, 2011. pp. 123-128. 10.4028/www.scientific.net/KEM.465.123
3
A. Morassi, Crack-induced changes in eigenparameters of beam structures. Journal of Engineering Mechanics. 119(9) (1993), pp. 1798-1803. 10.1061/(ASCE)0733-9399(1993)119:9(1798)
4
A.D. Dimarogonas, Vibration of cracked structures: a state of the art review. Engineering fracture mechanics. 55(5) (1996), pp. 831-857. 10.1016/0013-7944(94)00175-8
5
P. Cacciola, N. Impollonia, and G. Muscolino, Crack detection and location in a damaged beam vibrating under white noise. Computers & structures. 81(18-19) (2003), pp. 1773-1782. 10.1016/S0045-7949(03)00201-3
6
X.F. Yang, A.S.J. Swamidas, and R. Seshadri, Crack identification in vibrating beams using the energy method. Journal of sound and vibration. 244(2) (2001), pp. 339-357. 10.1006/jsvi.2000.3498
7
V.V. Matveev and A.P. Bovsunovsky, Vibration-based diagnostics of fatigue damage of beam-like structures. Journal of Sound and Vibration. 249(1) (2002), pp. 23-40. 10.1006/jsvi.2001.3816
8
T.G. Chondros, A.D. Dimarogonas, and J. Yao, A consistent cracked bar vibration theory. Journal of Sound and Vibration. 200(3) (1997), pp. 303-313. 10.1006/jsvi.1996.0718
9
K.R. Collins, R.H. Plaut, and J. Wauer, Free and forced longitudinal vibrations of a cantilevered bar with a crack. Journal of vibration and acoustics, 114(2) (1992), pp. 171-177. 10.1115/1.2930246
10
F. Vestroni and D. Capecchi, Damage detection in beam structures based on frequency measurements. Journal of Engineering Mechanics. 126(7) (2000), pp. 761-768. 10.1061/(ASCE)0733-9399(2000)126:7(761)
11
T.G. Chondros, A.D. Dimarogonas, and J. Yao, Vibration of a beam with a breathing crack. Journal of Sound and vibration. 239(1) (2001), pp. 57-67. 10.1006/jsvi.2000.3156
12
T.G. Chondros, A.D. Dimarogonas, and J. Yao, A continuous cracked beam vibration theory. Journal of sound and vibration. 215(1) (1998), pp. 17-34. 10.1006/jsvi.1998.1640
13
A.D. Dimarogonas and C.A. Papadopoulos, Vibration of cracked shafts in bending. Journal of sound and vibration. 91(4) (1983), pp. 583-593. 10.1016/0022-460X(83)90834-9
14
R.D. Adams, P. Cawley, C.J. Pye, and B.J. Stone, A vibration technique for non-destructively assessing the integrity of structures. Journal of Mechanical Engineering Science. 20(2) (1978), pp. 93-100. 10.1243/JMES_JOUR_1978_020_016_02
15
W.T. Springer, K.L. Lawrence, and T.J. Lawley, Damage assessment based on the structural frequency-response function. Experimental Mechanics, 28(1) (1988), pp. 34-37. 10.1007/BF02328993
16
J.A. Loya, L. Rubio, and J. Fernández-Sáez, Natural frequencies for bending vibrations of Timoshenko cracked beams. J. of Sound and Vibration, 290(3-5) (2006), pp. 640-653. 10.1016/j.jsv.2005.04.005
17
F.B. Sayyad and B. Kumar, Approximate analytical method for damage detection in beam. International Journal of Damage Mechanics. 21(7) (2012), pp. 1064-1075. 10.1177/1056789511429142
18
F.B. Sayyad, B. Kumar, and S.A. Khan, Approximate analytical method for damage detection in free-free beam by measurement of axial vibrations. International Journal of Damage Mechanics, 22(1) (2013), pp. 133-142. 10.1177/1056789512440897
19
H. Ayas, M. Chabaat, and L. Amara, Dynamic analysis of a cracked bar by the method of characteristics. International Journal of Structural Integrity. 10(4) (2019), pp. 438-453. 10.1108/IJSI-01-2018-0001
20
J.F. Doyle, Wave Propagation in Structures. 1997, Springer, New York, NY.
21
P.F. Rizos, N. Aspragathos, and A.D. Dimarogonas, Identification of crack location and magnitude in a cantilever beam from the vibration modes. Journal of sound and vibration. 138(3) (1990), pp. 381-388. 10.1016/0022-460X(90)90593-O
22
W.M. Ostachowicz and M. Krawczuk, Analysis of the effect of cracks on the natural frequencies of a cantilever beam. Journal of sound and vibration. 150(2) (1991), pp. 191-201. 10.1016/0022-460X(91)90615-Q
23
R. Ruotolo and C. Surace, Natural frequencies of a bar with multiple cracks. Journal of Sound and Vibration. 272(1-2) (2004), pp. 301-316. 10.1016/S0022-460X(03)00761-2
24
H. Tada, P. Paris, and G. Irwin, The analysis of cracks handbook. 2000, New York: ASME Press, vol. 2, pp. 1.
25
G.C. Sih, Some basic problems in fracture mechanics and new concepts. Engineering fracture mechanics. 5(2) (1973), pp. 365-377. 10.1016/0013-7944(73)90027-1
26
S.S. Rao, Vibration of continuous systems. 2007, New York : Wiley.
페이지 상단으로 이동하기