Magazine of Civil Engineering. 2020. 100(8). Article No. 10002
Magazine of Civil Engineering
journal homepage: http://engstroy.spbstu.ru/
ISSN 2712-8172
DOI: 10.18720/MCE. 100.2
Free vibration of axially loaded multi-cracked Timoshenko beams
Y.S. Al Rjoub*a, A.G. Hamad*
a Civil Engineering Department, Jordan University of Science and Technology, Irbid, Jordan b Jordan University of Science and Technology, Irbid, Jordan * E-mail: ysalrjoub@just.edu.jo
Keywords: cracked Timoshenko beam, shear deformation, rotary inertia, free vibration, buckling load, transfer matrix.
Abstract. In this paper, the free vibration of axially-loaded, multi-cracked Timoshenko beams with differing boundary conditions, namely, hinged-hinged, fixed-fixed, fixed-hinged, and fixed-free is studied. The cracked beam system is represented as several beam segments connected by massless rotational springs with sectional flexibility. Each segment is assumed to obey the Timoshenko beam theory. A simple transfer matrix method is used to derive the characteristic equation of the axially-loaded, multi-cracked beam with differing boundary conditions. The characteristic equation and corresponding mode shapes are a function of natural frequency, crack size and location, and physical parameters of the beam. In this paper, the effects of crack depth, number of cracks, position of cracks, axial load, shear deformation and rotary inertia on the dynamic behavior of multi-cracked beams are studied in detail. It is found that there is good agreement between the results obtained in this study and results available in the literature. Additionally, interesting observations overlooked by other researchers are obtained.
Many engineering structures suffer from mechanical vibrations, environmental effects, corrosion, extended service and cyclic loading, etc. that create structural defects, including cracks. These cracks change the dynamic properties of the structure, such as natural frequency and mode shape. A crack in a beam element introduces local flexibility due to strain energy concentrations in the vicinity of the crack tip under load. This flexibility changes the dynamic behavior of the structural member.
There has been much experimental and analytical research performed on the effects of cracking and damage on the safety of structures [1-5]. In this paper, the local flexibility method is used, which considers the cracked beam as several beam segments with the same properties connected by massless rotational springs. Local flexibility is used by [6] to estimate the crack stress intensity factor. The authors of [7-9] used Castigliano's theorem and fracture mechanics to compute the flexibility of the cracked region of a uniform beam with a transverse surface crack, which is a function of the ratio of crack depth to the height of the cross-section. The function was obtained experimentally, according to the measurements provided by [10]. This factor can also be obtained from [11]. Many researchers studied the free vibration of cracked Euler-Bernoulli and Timoshenko beams, both analytically and experimentally, [12-19]. The free vibration of cracked Timoshenko beams with various boundary conditions is studied by [20] using ultraspherical polynomials. The authors of [21] analyzed the free vibration of Euler-Bernoulli and Timoshenko beams resting on Pasternak and generalized elastic foundations. The free and forced vibration of cracked beams is investigated by [22-24]. Other researchers have studied the free vibration of cracked, functionally-graded Euler-Bernoulli and Timoshenko beams [25-29]. A dynamic stiffness method is used by the authors of [30] to study the free and forced vibration of multi-cracked, continuous, functionally-graded beams. An analytical solution of the forced vibration of cracked Timoshenko beams with damping is obtained by [31]. The free vibration of functionally-graded cracked Timoshenko beams is studied by [32-33]. A differential quadrature method is used by [34] to study the free and forced vibration of multi-cracked Timoshenko beams carrying a moving mass. A method is proposed by [35] which is based on the measurements of natural frequencies to detect a single crack in functionally-graded Timoshenko beam. In this paper, the transfer matrix method (TMM) is used
Al Rjoub, Y.S., Hamad, A.G. Free vibration of axially loaded multi-cracked Timoshenko beams. Magazine of Civil Engineering. 2020. 100(8). Article No. 10002. DOI: 10.18720/MCE.100.2
© © I
—This work is licensed under a CC BY-NC 4.0
1. Introduction
Magazine of Civil Engineering, 100(8), 2020
to study the free vibration of axially-loaded, cracked beams. This method was first introduced by [36], and used by other researchers [37-41] to study the free vibration of cracked beams. Other researchers solved the problem using the finite element method, [42-44].
According to the authors' knowledge, the effects of four cracks on the dynamic behavior of axially-loaded Timoshenko beams have not been studied in detail in the existing literature. In this paper, the dynamic behavior of a multi-cracked beam under axial loading including the shear deformation and rotary inertia and their coupling effects with differing boundary conditions with an arbitrary number of transverse cracks is investigated. The simple transfer matrix method is used to ascertain the natural frequencies and corresponding mode shapes of the beam. Some results obtained in this paper are compared to those results available in literature.
2. Methods
2.1. Theoretical model
Consider a multi-cracked beam with length l, cross-sectional area of width b and height h, and a number of single-sided, open cracks located at points xc1, xc2, ..., xci with sizes ac,1, ac,2, ..., ac,i, respectively, as shown in Fig. 1.
ac2 U a
->xcl
"» Xc3
1
K2p TTT
K
K
R1
R2
b
x
c2
P
P
x
x
z
У
Figure 1. A multi-cracked, axially loaded beam with i edge open cracks, restrained by linear and rotational springs at both ends.
The boundary conditions are represented with linear and rotational springs at both ends (Ki, Kri, K2, and Kr2). Timoshenko (1921) proposed a beam theory which adds the effects of shear and rotation to the
Oy( x, t)
governing equation of motion of an Euler- Bernoulli beam. The total slope - is a combination of the
dx
bending slope x, t) and the shear slope y(x,t) [46], which can be written as:
= y( x, t ) + y (x, t) (1)
Ox V ' xyK '
where x, t), yxy (x, t) ,and t are the bending slope, the shear slope, and the time parameter, respectively.
In [47], it was assumed that the component of axial load P on the cross section is acting on the plane of the total slope. Accordingly, the shear force can be written as:
V ( x,t)+P dy ( ' >
dx
'СУ (x,' )-r( x,t)
dx
(2)
where V is the shear force (along the y-axis), A is the area of beam cross section, G is the shear modulus and k' is the shear correction factor (assumed equal to 5/6 in this paper). The dynamic equilibrium equations for the free-body diagram shown in Fig. 2.b becomes:
, . dM (x, t) dy (x, t) d2w -V (x,t) +-P7y ' + pI ^ = 0 (3)
y ; dx dx dt2
And,
dV ( x, t) d2y (x, t )
dx
dt2
(4)
where M is the bending moment (about the z-axis), p is the material density, and I is the moment of inertia. Using Eqns. (3) and (4), one can obtain the governing equation for the transverse vibration of a uniform Timoshenko beam subjected to constant axial force as:
ii,',
EI
k 'GA — P k 'GA
d4 y ( x, t )
dxq
+ m + P
d2y ( x, t ) _
cT
mr
mEI d4y ( x, t) m2r2 d4y ( x, t )
kGA dx2st2 + kGA dt4
k 'GA -P k 'GA
= 0
4
d4y ( x, t ) _
dx2dt2
(5a)
In a similar manner, Eqn. (5a) can be written in terms of the deflection slope alone as:
i u,
EI
k 'GA -P k 'GA
Ad V( x, t ) dV( x, t )
dx4
+ m + P-
dx2
■ — mr
k 'GA -P k 'GA
A d V( x, t )
dx2dt2
mEI d4v(x,t) m2r2 d4v(x,t) kGA dx2dt2 + kGA dt4 ~
(5b)
where m is the mass per unit length, r is the radius of gyration; K'A is the effective shear area of the cross section.
It is clear that Eqn. (5a) reduces to Euler beam theory (EBT) when the shear deformation and rotary inertia are equal to zero. Then Eqn. (5a) become:
d2 ( y, t y pd2 ( y, t ) EI—v „ 7 + m—^ t + P—^ ~ = 0
dx
4
dx2
dxz
A simple transfer matrix method is used to get the general form of the characteristic equation for the axially-loaded cracked beam, which is a function of frequency, axial load, the locations and sizes of the cracks, boundary conditions, and the geometrical and physical parameters of the beam. It is assumed that a beam with n segments is joined by massless, rotational springs, as shown in Fig. 2.
Figure 2. Beam elements joined by massless, rotational spring.
The free vibration solution for a Timoshenko beam (TBT) subjected to axial load can be found using the method of separation of variables as:
y( x, t ) = Y ( x)e
ia>t
(6a)
iot
y/( x, t) = p( x)eict (6b)
where, Y(x), \y(x) , cc, i, and x are the normal function of y(x, t) , the normal function of x, t) , the circular frequency, the imaginary number, and the spatial position on the x-axis, respectively. Substituting Eqn. (6a) into Eqn. (5a) leads to:
EI
k 'GA - P ^d 4 Y (x)
d 2 Y (x)
V
k 'GA
2 2 mr a
V k'GA J
+ maY (x) + P-
dx dx
22
r k'GA - P^ d2Y(x) mEIa2 d2Y(x) m2r2a4
dx2
k'GA dx2
k 'GA
Y (x ) = 0
And, substituting Eqn. (6b) into Eqn. (5b) yields:
r Ur^ *
EI
k 'GA - P V( x) 2 d V( x) v ; + ma + P-^
V
k 'GA
J
dx4
dx2
2 2 mr a
k'GA - P VV(x) mEIa2 dV(x) m2r2a4 ,, n
- -^---^ +-V x ) = 0
k 'GA - 2 - 2 v' V '
J
Ox2 k'GA Ox2 k'GA The general-form solutions of Eqn. (7a) and (7b) are:
Y(x) = A1 cos(ax) + A2 sin(ax) + A3 cosh(^x) + A4 sinh(^x)
y/(x) = m1A cos(ax) + m1 A2 sin(ax) + cosh(^x) + m2A4 sinh(^x)
where,
(7a)
(7b)
(8a) (8b)
a =
^22 2
mr a ma + ■
P
2
EI k' AG£ EI 4
2 2 2 mr2a2 ma2 P
-+-+-
EI k' AG£ EI 4
- 4
22 mr2a 2
ma
2
EIk' GA4 EI 4
ß =
mr2a2 ma2 P +-+ -
EI k' AGC EI 4
1
+ — 21
2 2 2
mr a ma P +-+ ■
v
EI k' AGC EI 4
4
^ 2 2 mr2 a2
ma
2 A
EIk' GA4 EI 4
e k' GA - P pAa where 4 = —- . -. .—. m =
pAaz
. =--a. m2 =-
k' GA 1 k' GAa k' GAß
+ ß and Ai to A4 are unknown constants.
The slope ^(x) , bending moment M(x), and shear force V(x) can be obtained as:
0( x) =
dY ( x) dx
= mi [-Ai sin(ax) + A2 cos(ax)] + m2 [-A3 sinh(ßx) - A4 cosh(ßx)]
M(x) = n1 [—A1 cos(ax) - A2 sin(ax)] + n2 [A3 cosh(^x) + A4 sinh(^x)] V(x) = b [A1 sin(ax)—A2 cos(ax)] + b2 [A3 sinh(^x) + A4 cosh(^x)]
9 9 9 9
where n1 = EIm1 a, n2 =— EIm2¡, b = EIm1 a — m^_pIa , b2 =—EIm2a — m2pIa . The [ S(x) ] matrix can be obtained from Eqn. (8a) and (9a-c) as:
(9a)
(9b) (9c)
2
2
2
[ s ( x)] =
cos(ax) -my sin(ax) -ny cos(ax) b sin(ax)
sin(ax) my cos(ax) - ny sin(ax) -by cos(ax)
cosh(^x) -m2 sinh(^x) n2 cosh(^x) b2 sinh(^x)
sinh(^x)
-m2cosh(^x) n2 sinh(^x) b2 cosh(^x)
Y+y( xc ) "y 0 0 0" Y ( xc )
°i+y( xc ) 0 y Ci 0 O ( xc )
Mi+y( xc ) 0 0 y 0 M ( xc )
JM( xc ) . _ 0 0 0 y Vi ( xc ) ,
The compatibility requirements at the crack position, xc, which are the continuity of the deflection, bending moment, shear force, and the discontinuity of slopes for segments i and i+1, can be written in matrix form as:
(11)
where Q is the flexibility coefficient at the ith crack, and depend on the load case of the cracked element.
2.2. Local flexibility due to a crack
An open crack on an elastic structure can be considered a source of local flexibility due to the strain energy concentration at the surrounding area of the crack tip. The idea of replacing a crack with a massless spring is presented to establish the relation between the strain energy concentration and applied loads. The flexibility coefficients are expressed in terms of stress intensity factors, utilizing Castigliano's theorem. Generalized loading conditions for abeam element of circular or rectangular cross section with a transverse surface crack are shown in Fig. 3.
Figure 3. Model of the cracked element.
The crack has a tip line parallel to the z-axis and the bar is loaded by shear forces P2, P3, Ps, and P9; bending moments P4, P5, P11, and P12; axial forces Pi and P7; and torsional moments P4 and P10. According to Castigliano's theorem, the additional displacement caused by the crack is given as:
ui =
dp "dp
J JdA
Ac
(12)
where U, Ac, and J are the strain energy due to the crack, the crack section, and strain energy density function, respectively.
The strain energy density function can be expressed according to [14] for plane stress as: 1
J =
V E J
f 6
Y
I Kn
V i=y J
6
Y
6
I Km + (y + V) I Km
V i =y J V i =y J
Y
(13)
where, v is Poisson's ratio and Kn, Km, Knn are the crack stress intensity factors for opening mode, shearing mode, and tearing mode, respectively; I = 1,2,...,6.
Stress intensity factors can be calculated following [14] as:
.__P
KI1 = CT1VnaF^a / A),ö"1 = —1
bh
I— 12 P zz
KI5 = <5\ naF1 (a / h), a5 — —r-2—
b h
KI6 = a647aF2 (a / h), < =
6P 6 bh2
KI2 = KI3 = KI4 = 0
KII2 = (a / hX <2 = kPl
bh
KII4 = a447aFII(a / h),<4 = —
bh
(14)
KII1 = kii 3 = kii 5 = kii 6 = 0
K^ 3 = <j3y[TaFIII (a / h), <3
kP3 bh
where
.__p
KIII4 = <4III > naFIII(a / h), <4III = TT
bh
KIII1 = KIII2 = KIII5 = KIII6 = 0
2h
F (a / h) =— tan I — 1 \7Ü V 2h
0.752 + 2.021 a | + 0.37 na I v h
f
1 - sin
na 2h
cos
na 2h
i
„ , ,12h (na
F2(a/h) = J—tan T7 V na V 2h
0.923 + 0.199
1 - sin
V
f
' na ^
v 2h jj
i
cos
na 2h
(15)
Fu (a / h) =
1.122-0.561| ^1 + 0.085(0\ + 0.18(a^
1 -
FIII(a / h) =
2h ' — tan na
na ~2h
From Eqns. (12). (14) and (15). the flexibility coefficients can be calculated as:
n b/2
Cn = (1 -v2)-^J dz Eb 0
-b/2
3
y— d b/2
Ç5 - C51 - (1 -v2)—-[JM2(^] J zdz
Eb 0 -b/2
d
2
d b/2 Ci6 - C61 - (1 -v2)—J^^d^] J dz
Eb h 0 -b/2
ju2 d b/2
C22 = (1 -v2)"^T[J^/2(0dfl J dz
0 -b/2
r>! d b/2
C24 - C42 - (1 -v2) —[JF(fldfl J dz
0 -b/2
9fr2 d b/2 C33 -(1 -v2)--^2-[J^F/2/(fldfl J dz
Eb2
0 -b/2 d b/2
C34 - C43 - (1 -v2)2^[J^F/2/(fldfl J dz (16)
Eb 0 -b/2
^ (p b/2 c44 -(1 -v2)-—(^))d^] j dz
0 -b/2
2 ■> ^
C55 - (1 - v2 ) [ j ^Fj2 (^)d^] bJ2 z2dz
Eb
0 -b/2
144.77- d b/2
C56 -C65 -(1 -v2)—— [JtF2(0)F2(tW] J zdz Eb h 0 -b/2
17TT d b/2
C66 -(1 -v2)—^J^22(^] J dz
Eb h 0 -b/2
1 . a ,
where ^ = 1 + v, <p = — and k is a shape coefficient for a rectangular cross section.
h
In this study, the flexural vibration of the beam is the only significant load in the equations above. Therefore, the contribution of the axial load component to the strain energy of the system is very small. Hence, the only component in the flexibility matrix is C66.
2.3. Transfer matrix method
Equation (11) can be written as:
{q+i( x )}=[q ]{q (x )} (17)
At the location of open edge crack, the two beam segments i and i+1 are connected by a massless, rotational springs, as shown in Fig. 2. Equation (6a) and Eqn. (8a) and Eqns. (9a) to (9c) can be written in matrix form for both segments i, and i+1, respectively, as:
' Y (x) A1,i'
0 (x) Mi (x) ' = [ S ( x) ]• A2,i A3,;
. V (x) . A4,i,
(18a)
' Y+1( x)" A1,i+1"
O+1( x) Mi+1( x ) ' = [ S (x) ]• A2,i+1 A3,i+1
Vi+1( x) _ A4,i+1 ^
(18b)
Equation (18a) can be rewritten for a segment, Xi-i < x < xt, as:
{Qii x)} = [S (x)]{4}
(19)
where,
{Qt (x)} =
and {A } =
1, i Al,i Ai
A4,i
Y+i( xc ) A
+l( X ) Mi+l( xc ) V+i( xc ) ,
The boundary conditions of the beam using equilibrium equations can be expressed as:
M (0) = Kr0(0)
V (0) = - KlY (0)
M (l) = - KR20(l)
V (l) = K2Y (l)
(20a) (20b) (20c) (20d)
In order to calculate the natural frequencies and corresponding mode shapes, the characteristic equation can be derived as in [48].
The continuity conditions of the deflection, slope, moment, and shear at x = 0 and x = l leads to:
(21)
YM(l) Y1(0)
Ol+1(l) • = [T ]• 0(0)
-KR 20i+1(l) Kr0(0)
K 2^(0 . - KY1(0)
where [T ] =
T1 To T1 T
11 T21 T31
'12 T22 T32
13 T23 T33
'14 T24
T34
T41 T42 TAi T
43
44
Rewriting equation (21) as:
if) ] \0i+1(l) J
T11 " K1T14
T12 + KR1T13
T 21 K1T24 T22 + KR1T23
(22a)
>
f-KR2ei+s ) ] I K2Yi+l(l) I
T31 - K1T34 T 41- K1T44
T32 + KR1T33 T42 + 1T43
Y1(0) 1 0(0) J
(22b)
The characteristic equation can be obtained by substitution of Equation (22a) into Equation (22b).
F1F4 - F2F3 = 0 (23)
where,
F1 — -KR2T21 + K1KR2T24 - T31 + K1T34
F2 - -KR2T22 - KR1KR2T23 - T32 - KR1T33
F3 - K2T11 - K1K2T14 - T41 + K1T44
F4 — K2T12 + K2KR1T13 - T42 - KR1T43 Characteristic Eqn. (23) can be derived for different boundary conditions as:
a. For fixed-free boundaries (^ — KRl — œ) ,
T34T43 - T33T44 = 0
b. For pinned-pinned boundaries (— — œ) ,
T34T12 _ T32T14 = 0
c. For fixed-pinned boundaries, (K — K2 — K^ — œ) ,
T34T13 ~ T33T14 = 0
d. For fixed-fixed boundaries
ies, | Kj —
— K — KR ! — Kp —
— œ) ,
T24T13 _ — 0
23T14
e. For free-free boundaries,
T31T42 _ — 0
32T41
2.4. The corresponding mode shapes
In order to obtain the mode shape of the axially-loaded, cracked Euler beam, the unknown constants {A } for each sub-segment of the cracked beam can be determined. At first, the unknown coefficients for the
first segment, {A} , are determined. Then, the unknown coefficients for the other segments may be obtained, using the following procedure.
1. At x = 0, Eqn. (19) can be written as:
(24)
' Y1(0) 1 1 0 1 0
0(0) 0 a 0 P A2,1
Kr10(0) - EI a2 0 EI 01 0 A3,1
- K1Y1(0)J 0 -(EIa3 - Pa) 0 ( EIP3 + PP . A4,1 j
Eqn. (24) can be rewritten as:
I
- EIa2
10 10" 0 a 0 ß
0
A
1,1
A
-2,1
A
3,1
0 -(EIa - Pa) Substituting Eqn. (25a) into Eqn. (25b) yields:
011 012 013 014 021 022 023 024
A4,1
EIß1 0
(25a)
0
(EIß3 + Pß)
A
1,1
A
2,1
A
3,1
A
4,1
A
A
1,1 2,1
A
A
3,1 4,1
= 0
where.
Q11 = n1 . 012 = KR1m1 . 013 = -»2 . Q14 = -KR1m2 .
Ö21 = Q22 = h . Q23 = -K1 . 024 = -b2 2. At X = l. substituting Eqns. (10). (20c). (20d). and (19) into Eqn. (21) yields:
Y+1(l) 3+1(l)
[-Kr 2^+1(l)
I
>11 - EIa2Tu (Tua-Tu(EIa3- Pa)) Tn + EITxß2 {TX2ß + ^(Elß3 + Pß)) T21 - EIT23a2 (T22a-T24(EIa3- Pa)) T2X + EIT2ß {T22ß + T24(EIß3 + Pß))
T31 - EIa2T33 (T32a-T34(EIa3- Pa)) T31 + E^ß2 (T32ß + T34(EIß3 + Pß)) T41 - EIT43a2 (T42a - T44(EIa3- Pa)) T41 + EIT43ß2 (T42ß + T44(EIß3 + Pß))
43a
Substituting Eqn. (38a) into Eqn. (38b) leads to:
Ö31 032 033 034 041 042 043 044
A
1,1
A
2,1
A
3,1
A
4,1
= 0
where.
Q31 = KR2 (T21 - n1T23 ) - T31 + n1T3
33
Ö32 = -KR2 (m1T22 - b1T24 ) - (m1T32 - b1T34 )
Q33 = KR2 (T21 + n2T23 ) - T31 - n2T3
33
Ö34 = KR2 (-m2T22 + b2T24 ) - (-m2T32 + b2T34 )
041 = K2 (T11 - n1T13) T41 + n1T4
43
042 = K2(m{r12 -bT14) -(m^ -bT4)
Ö43 = K2 (T11 + n2T13 ) - T41 - n2T
A1,1 A2,1 A3,1 A4,1
A1,1 A2,1 A3,1 A4,1
(25b)
(26)
(27)
(28a)
(28b)
(29)
43
044 = K2(-m{T12 + b2T14) - (-m2^2 + ¿¿Ttf)
Combining Eqn. (26) and Eqn. (29) yields:
Q11 Q12 Q13 Q14 41,1
Q21 Q22 Q23 Q24 42,1
Q31 Q32 Q33 Q34 43,1
Q41 Q42 Q43 Q44 _ 44,1
The unknown coefficients of Eqn. (31) correspond to the mode shapes of the first crack position, i.e. at x = Xc1. For other crack locations, the unknown coefficients can be obtained using the coefficients of the previous crack as:
{ A }— [S ( Xi-1)]-1 [Ci-1 ] [ S (Xi-1)] {4.-1} (32)
3. Results and Discussion
3.1. Validation of results
The present results obtained are compared with those found by [45]. Three different boundary conditions are used, namely, hinged-hinged, fixed-hinged, and fixed-fixed. Also, two different values of P/Pcr
are used, which are 0.2 and -0.2. Pcr is the critical buckling load of the intact beam and it is equal to
1 + - PE
k GA
, where Pe is the Euler buckling load. Values of l = 1 m, b = 0.1 m, h = 0.2 m, E = 200 GPa, p = 7850kg/m3, and v = 0.3 are used as the geometrical and material properties of the intact beam. Fig. 4 a-c show the first non-dimensional natural frequency ratio / &in obtained in this study versus the
non-dimensional crack location x/l for a beam with a/h = 0.1, a/h = 0.3, and a/h = 0.5, for hinged-hinged, fixed-hinged, and fixed-fixed beams, respectively, compared with those obtained by [45]. As you can see, the obtained results agree well with those reported by [45].
Figure 4. The first non-dimensional natural frequencies of a) hinged-hinged beam with single crack with a/h = 0.1, b) fixed-fixed beam with single crack with a/h = 0.3, and c) fixed-hinged beam with single crack with a/h = 0.5. The curves are obtained for P / Pcr values of -0.2, and 0.2
(from top to bottom).
3.2. Effects of One Crack and Axial Load on Natural Frequencies
In all subsequent sections, the geometrical and material properties of the intact beam are the same as those used in validation section, unless otherwise noted. The first three non-dimensional natural frequencies of a beam with a single crack versus crack location, for hinged-hinged, fixed-fixed, fixed-hinged, and fixed-free
Magazine of Civil Engineering, 100(8), 2020
boundary conditions, respectively, are shown in Fig. 5, 6, 7, and 8. In all figures, figure a) represents P/Pcr = 0, figure b) represents P/Pcr = -0.3, and figure c) represents P/Pcr = 0.3. Also, in all figures, eight values of aspect ratio a/h are used, 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, and 0.7. It is clear that the natural frequencies do not change with differing values of aspect ratio (a/h) at the location of zero-crossing of the mode shape and at critical points for a given frequency but do vary with differing boundary conditions. As an example, the zero-crossing points for hinged-hinged boundary conditions are located at 0.331, and 0.66/ for the third, non-dimensional natural frequency. The critical points are located at 0.11/ and 0.89/ for fixed-fixed boundary conditions and at 0.12/ for fixed-hinged boundary conditions for the second, non-dimensional natural frequency. An interesting observation overlooked by researchers is that the locations of critical points for the first non-dimensional frequencies are shifted toward the center of the beam when the beam is axially loaded in compression and toward the fixed supports of the beam when the beam is axially loaded in tension for fixed-fixed and fixed-hinged boundaries. Those locations are at 0.23/, 0.77/, and at 0.21/, 0.79/, for compression-and tension-loaded beams, respectively, for fixed-fixed boundaries, and 0.27/ and 0.25/ for tension- and compression-loaded beams, respectively, for fixed-hinged boundaries. As most researchers have found, the beam's natural frequencies are independent of crack location. Also, the natural frequency decreases with increasing aspect ratio. The first non-dimensional frequencies are most affected by the severity of cracks and the magnitude of the axial force. For axially loaded beams, since the compressive force softens the beam and tensile force stiffens the beam, the non-dimensional natural frequencies for all modes of the beam under compressive force are smaller than those of the same beam under tensile load for all boundary conditions. The first non-dimensional natural frequency shows a maximum 14 % and 20 % increase when the cracks are located near the supports and 16 % and 29 % decrease when the cracks are located at mid-span of the beam for P/Pcr = 0.3 and P/Pcr = -0.3, respectively, for an aspect ratio of 0.5, for an axially-loaded, hinged-hinged beam. for axially-loaded fixed-fixed beams, the first non-dimensional natural frequency shows a maximum 12 % increase when the cracks are located near the supports and 15 % decrease when the cracks are located at mid-span of the beam for P/Pcr = -0.3 and P/Pcr = 0.3 respectively, for aspect ratio of 0.1, and a maximum 15 % increase and 19 % decrease for P/Pcr = -0.3 and P/Pcr = 0.3 respectively, for an aspect ratio of 0.5. The values of percent increase or decrease are increased as the value of aspect ratio increases. For axially-loaded, fixed-hinged beam, the first non-dimensional natural frequency shows a maximum 13 % increase and 16 % decrease for P/Pcr =-0.3 and P/Pcr = 0.3 respectively, for an aspect ratio of 0.1, and a maximum 16 % increase and 20 % decrease for P/Pcr = -0.3 and P/Pcr = 0.3 respectively, for an aspect ratio of 0.6. Again, for axially-loaded fixed-free beam, the first non-dimensional natural frequency shows a maximum 15 % decrease and 13% increase for P/Pcr =-0.3 and P/Pcr = 0.3 respectively, for an aspect ratio of 0.1. A maximum 24 % decrease and 16% increase are obtained for P/Pcr = -0.3 and P/Pcr = 0.3 respectively, for an aspect ratio of 0.5. Another interesting point overlooked by researchers that for fixed-free beam, the first non-dimensional natural frequencies is not affected by crack location, aspect ratio of the crack, or the magnitude of axial load beyond 0.86/, 0.8/, and 0.92/ of the beam for zero axial load, tensile load and compressive load, respectively compressive load, respectively.
Figure 5. The first three non-dimensional natural frequencies of a hinged-hinged beam with a single crack for a) P/P«=0, b) P/P«=-0.3, and c) P/Pc=0.3. The curves are obtained for a/h values of 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, and 0.7 (from top to bottom), and c) P/Pc=0.3 with values of a/h= 0, 0.1, 0.2, 0.3,
0.4, and 0.5.
Figure 6. Same as Fig. 5 but with fixed-fixed boundary and values of a/h = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6 and 0.7 for c).
Figure 7. Same as Fig. 5 but with fixed-hinged boundary
and values of ah = 0, 0.1, 0.2, 0.3, 0.4, 0.5, and 0.6 for c).
a) —T—T—T—T— —T—T—T—T— —T—T—T—n
0.3-'-'-'-'- 0.3-'-'-'-'- 0.3-'-'-'-'-
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
Magazine of Civil Engineering, 100(8), 2020
b) 1.21-1-1-1-1- 1.21-1-1-1-1- 1.21-1-1-1-1-
Figure 8. Same as Fig. 5 but with fixed-free boundary and values of a/h = 0, 0.1, 0.2, 0.3, 0.4, and 0.5 for c).
3.3. Effects of mu/tip/e cracks and axia/ /oad on natura/ frequencies
Fig. 9 and 10 show the first three non-dimensional natural frequencies of a beam with four cracks versus crack location, for hinged-hinged and fixed-free boundaries, respectively. For brevity, the fixed-fixed and fixed-hinged results are not shown here. Five values of aspect ratio a/h are used in all figures (0, 0.1, 0.2, 0.3, and 0.4). Three cracks are fixed at locations of maximum variation of non-dimensional frequencies: 0.15/, 0.5/, and 0.85/, and at 0.01/, 0.6/, and 0.8/ for hinged-hinged and fixed-free boundaries, respectively, while the fourth crack moves along the beam. As seen, the presence of four cracks in the beam lowers the non-dimensional natural frequencies, with or without axial load relative to the case of a single crack. A maximum 17 % and 15 % decrease of the first non-dimensional frequencies is observed for an aspect ratio of 0.4 when the fourth crack is located closer to the supports and mid-span of the beam, for hinged-hinged boundary conditions-compare Fig. 5a) with Fig. 9a). For fixed-free beam, a maximum 14 % and 34 % decrease of the first non-dimensional frequency is observed for an aspect ratio of 0.4 when the fourth crack is closer to the fixed and free supports, respectively-compare Fig. 8 a) with Fig. 10 a). It can be concluded here that the strongest and weakest beam in terms of fundamental natural frequencies is fixed-fixed and fixed-free beam, respectively.
0 0.2 0.4 0.6 0.8 1 '0 0.2 0.4 0.6 0.8 1 ' 0 0.2 0.4 0.6 0.8 1
0 0.2 0.4 0.6 0.8 1 '0 0.2 0.4 0.6 0.8 1 ' 0 0.2 0.4 0.6 0.8 1
X 4 X 4 X 4
C C C
Figure 9. The first three non-dimensional natural frequencies of a hinged-hinged beam with four cracks for a) P/Pc=0, b) P/Pcr=-0.3, and c) P/Pc=0.3. The curves are obtained for a/h values of 0, 0.1, 0.2, 0.3, and 0.4 (from top to bottom).
Figure 10. Same as Fig. 9 but with fixed-free boundary.
3.4. Effects of number of cracks on natural frequencies and corresponding mode shapes.
Fig. 11 and 12 show the effect of number of cracks on the first three non-dimensional natural frequencies of a beam for fixed-fixed, and fixed-hinged boundaries, respectively. Again, for brevity, the results of hinged-hinged and fixed-free are not shown here. In all figures, a/h = 0.3 is used, figure a) represents P/Pcr = 0, figure b) represents P/Pcr = -0.3, and figure c) represents P/Pcr = 0.3, the solid line represents a beam with one moving crack, the dashed line represents a beam with two cracks: one crack at location 0.2l and one moving
crack along the beam, and the dotted line represents a beam with three cracks, cracks located at 0.2/ and 0.4/ and one moving crack along the beam. Increasing the number of cracks leads to a decrease in the non-dimensional natural frequency for all loading cases. Fig. 13 and 14 show the effects of the number of cracks on the first three mode shapes of a beam with hinged-hinged and fixed-free boundaries, respectively. The results of fixed-fixed and fixed-hinged are not shown. Again, in all figures, a/h = 0.3 is used, figure a) represents P/Pcr = 0, figure b) represents P/Pcr = -0.3, and figure c) represents P/Pcr = 0.3. It is clear that the number of cracks slightly affects the mode shapes.
Figure 11. Effects of number of cracks on the first three non-dimensional natural frequencies of fixed-fixed beam for a) P/Pc=0, b) P/Pcr=-0.3, and c) P/Pcr=0.3.
Figure 12. Same as Fig. 11 but with fixed-hinged boundary.
Figure 13. Effects of number of cracks on the first three mode shapes of hinged-hinged beam for a) P/Pcr=0, b) P/Pcr=-0.3, and c) P/Pcr=0.3. (Blue line represents one crack, green line represents
two cracks, and red line represents three cracks)
Figure 14. Same as Fig. 13 but with fixed-free boundary.
3.5. Effects of shear deformation and rotational inertia on natural frequencies
Fig. 15 and 16 show the effects of shear deformation and rotational inertia on the first three non-dimensional natural frequencies for fixed-fixed and fixed-hinged boundary conditions, respectively. The results
of hinged-hinged and fixed-free are not shown. The y-axis is Qy = coic^p/ E of the cracked beam and the x-axis is the position of the crack. In all figures, figure a) represents P/Pcr = 0, figure b) represents P/Pcr = -0.3, and figure c) represents P/Pcr = 0.3. Three values of slenderness ratio are used 5, 10, and 20. The value of a/h used is 0.3. It can be seen that the non-dimensional natural frequencies decrease as the value of slenderness increases. Also, it can be seen that the fundamental non-dimensional natural frequencies obtained by EBT and TBT gradually coincide at L/h greater than 10. A remarkable difference is observed on the fundamental non-dimensional natural frequencies obtained by EBT and TBT at L/h less than 10. At higher modes, the effects of L/h become more noticeable for all boundary conditions. In all figures, the results of L/h = 20 are not shown because the EBT and TBT results coincide with each other.
0.41-'-'-'-'--11-'-'-'-'--21-'-'-'-'-
0 0.2 0.4 0.6 0.8 1 0 0.2 D.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
Figure 15. The first three non-dimensional frequencies of fixed-fixed beam with a single crack obtained by Euler- Bernoulli model with L/h = 5, Timoshenko model with L/h = 5, Euler-Bernoulli model with L/h = 10, and Timoshenko model with L/h = 10 (represented by solid red, dashed red, solid green, and dashed green lines, respectively); a) P/Pcr=0, b) P/Pcr=-0.3, and c) P/Pcr=0.3.
4. Conclusions
A simple transfer matrix method is used to study the vibrational behavior of axially-loaded, multi-cracked Timoshenko beams with differing boundary conditions. The effects of slenderness ratio, axial load, aspect ratio, a/h, boundary conditions, and number of cracks on the non-dimensional natural frequencies and the corresponding mode shapes are investigated. The following conclusions can be drawn:
1. The tensile force stiffens the beam while the compressive force softens the beam.
2. For the first mode, the severity of the cracks and the magnitude of the axial load affect the non-dimensional frequencies.
3. A decrease in the natural frequency is observed as the number of cracks increases, due to increased flexibility of the beam.
4. For specific crack location, the natural frequency decreases as the aspect ratio, a/h, increases.
5. At the location of the zero-crossing of the mode shape, and at critical points along the beam, the natural frequencies do not change with differing aspect ratios.
6. For the first non-dimensional frequencies, the compressive force shifts the locations of critical points toward the center of the beam. The tensile force shifts them toward the supports for fixed-fixed and fixed-hinged boundaries.
7. The first non-dimensional natural frequency is not affected by the crack location, aspect ratio of the crack, or the value of axial load beyond a certain location of the fixed-free beam.
8. A remarkable difference of the non-dimensional natural frequencies obtained by the Euler-Beam theory and Timoshenko beam theory is observed at values of L/h less than 20.
References
1. Liebowitz, H., Claus, W.D. Failure of notched columns, J. Eng. Fract. Mech. 1986. 1. Pp. 379-383. DOI: 10.1016/0013-7944(68)90010-6
2. Dimarogonas, A.D. Vibration of cracked structures: a state of the art review. Engineering Fracture Mechanics, 1996. 55 (5). Pp. 831-857. DOI: 10.1016/0013-7944(94)00175-8
3. Zheng, T., Ji, T. Equivalent representations of beams with periodically variable cross-sections. Engineering Structures. 2011. 33 (3). Pp. 706-719. DOI: 10.1016/j.engstruct.2010.11.007
4. Christides, S., Barr, A.D.S. Torsional vibration of cracked beams of non-circular cross-section. International Journal of Mechanical Sciences. 1986. 28 (7). Pp. 473-490. DOI: 10.1016/0020-7403(86)90067-6
5. Chondros, T.G., Dimarogonas, A. D., Yao, J. A continuous cracked beam vibration theory. Journal of Sound and Vibration. 1998. 215(1). Pp. 17-34. DOI: 10.1006/jsvi.1998.1640
6. Dimarogonas, A.D. Vibration Engineering. West Publishing Company, St. Paul, 1976. 568 p.
7. Liebowitz, H. Vanderveldt, H., Harris, D.W. Carrying capacity of notched columns. International Journal of Solids and Structures. 1967. 3(4). Pp. 489-500. DOI: 10.1016/0020-7683(67)90003-0
8. Liebowitz, H., Claus Jr., W.D. Failure of notched columns. Engineering Fracture Mechanics. 1968. 1(2). Pp. 379-383. DOI: 10.1016/0013-7944(68)90010-6
9. Okamura, H., Liu, H.W., Chu, C.S., Liebowitz, H. A cracked column under compression. Engineering Fracture Mechanics. 1969. 1(3). Pp. 547-564. DOI: 10.1016/0013-7944(69)90011-3
10. Brown, W.F., Srawley, J.E. Plane Strain Crack Toughness Testing of High Strength Metallic Materials. ASTM International. 1966. Pp. 1-129. DOI: 10.1520/STP44663S
11. Tada, H., Paris, P.C., Irwin, G.R. The stress analysis of cracks handbook, 3rd edn. ASME Press, New York .2000. 698 p. DOI: 10.1115/1.801535
12. Fernandez-saez, J., Navarro, C. Fundamental frequency of cracked beams in bending vibrations: an analytical approach. Journal of Sound and Vibration. 2002. 256(1). Pp. 17-31. DOI: 10.1006/jsvi.2001.4197
13. Yoon, H., Son, I., Ahn, S. Free vibration analysis of Euler-Bernoulli beam with double cracks. Journal of Mechanical Science and Technology. 2007. 21(3). Pp. 476-485. DOI: 10.1007/BF02916309
14. Loya, J.A., Rubio, L., Fernandez-saez, J. Natural frequencies for bending vibrations of Timoshenko cracked beams. Journal of Sound and Vibration. 2006. 290(3-5). Pp. 640-653. DOI: 10.1016/j.jsv.2005.04.005
15. Zhong, S.C., Oyadiji, S.O. Analytical predictions of natural frequencies of cracked simply supported beams with a stationary moving mass. Journal of Sound and Vibration. 2008. 311(1-2). Pp. 328-352. DOI: 10.1016/j.jsv.2007.09.009
16. Naguleswaran, S. Vibration and stability of an Euler-Bernoulli beam with up to three-step changes in cross-section and in axial force. International Journal Mechanical Sciences. 2003. 45(9). Pp. 1563-1579. DOI: 10.1016/j.ijmecsci.2003.09.001
17. Gomes, H.M., de Almeida, F.J.F. An analytical dynamic model for single-cracked beams including bending, axial stiffness, rotational inertia, shear deformation and coupling effects. Applied Mathematical Modelling. 2014. 38(3). Pp. 938 -948. DOI: 10.1016/j.apm.2013.07.019
18. Aydin, K. Vibratory characteristics of Euler-Bernoulli beams with an arbitrary number of cracks subjected to axial load. Journal of Vibration and Control. 2008. 14(4). Pp. 485-510. DOI: 10.1177/1077546307080028
19. Neves, A.C., Simoes, F.M.F., da Costa, A.P. Vibrations of cracked beams: Discrete mass and stiffness models. Computers and Structures. 2016. 168. Pp. 68-77. DOI: 10.1016/j.compstruc.2016.02.007
20. Kwanghun, K., Sok, K.,Kyongjin, S., Pak, C., Kwangbok, H. A modeling method for vibration analysis of cracked beam with arbitrary boundary condition. Journal of Ocean Engineering and Science. 2018. 3(4). Pp. 367-381. DOI: 10.1016/j.joes.2018.11.003
21. Batihan, A.Q., Kadioglu, F.S. Vibration Analysis of a Cracked Beam on an Elastic Foundation. International Journal of Structu ral Stability and Dynamics. 2016. 16(05). 1550006 p. DOI: 10.1142/S0219455415500066
22. Khorram, A., Rezaeian, M., Bakhtiari-Nejad, F. Multiple cracks detection in a beam subjected to a moving load using wavelet analysis combined with factorial design. European Journal of Mechanics A/solids. 2013. 40. Pp. 97-113. DOI: 10.1016/j.euromechsol.2012.12.012
23. Heydari, M., Ebrahimi, A., Behzad, M. Forced vibration analysis of a Timoshenko cracked beam using a continuous model for the crack. Engineering Science and Technology, an International Journal. 2014. 17(4). Pp. 194-204. DOI: 10.1016/j.jestch.2014.05.003
24. Yang, J., Chen. Y., Xiang, Y., Jia, X.L. Free and forced vibration of cracked inhomogeneous beams under an axial force and a moving load. Journal of Sound and Vibration. 2008. 312(1-2). Pp. 166-181. DOI: 10.1016/j.jsv.2007.10.034
25. Khiem, N.T., Lien, T.V., Ninh, V.T.A. Natural Frequencies of Multistep Functionally Graded Beam with Cracks. Iranian Journal of Science and Technology, Transactions of Mechanical Engineering. 2019. 43, Pp. 881-916. DOI: 10.1007/s40997-018-0201-x
26. Lien, T.V., Duc, N.T., Khiem, N.T. Free Vibration Analysis of Multiple Cracked Functionally Graded Timoshenko Beams. Latin American Journal of Solids and Structures. 2017. 14(9). Pp. 1752-1766. DOI: 10.1590/1679-78253693
27. Kou, K.P., Yang, Y. A meshfree boundary-domain integral equation method for free vibration analysis of the functionally graded beams with open edged cracks. Composites Part B: Engineering. 2019. 156. Pp. 303-309. DOI: 10.1016/j.compositesb.2018.08.089
28. Khiem, N.T., Tran, H.T., Nam. D. Modal analysis of cracked continuous Timoshenko beam made of functionally graded material. Mechanics Based Design of Structures and Machines. 2019. DOI: 10.1080/15397734.2019.1639518
29. Chen, W., Chang, H. Vibration Analysis of Functionally Graded Timoshenko Beams. International Journal of Structural Stability and Dynamics. 2018. 18(01). 1850007 p. DOI: 10.1142/S0219455418500074
30. Lien, T.V., Duc, N.T., Khiem, N.T. Free and forced vibration analysis of multiple cracked FGM multi span continuous beams using dynamic stiffness method. Latin American Journal of Solids and Structures. 2019. 16(2). Pp. 1-26. DOI: 10.1590/1679-78255242
31. Chen, B., Zhao, X., Li, Y.H., Guo, Y. Forced vibration analysis of multi-cracked Timoshenko beam with the inclusion of damping by virtue of Green's functions. Applied Acoustics. 2019. 155. Pp. 477-491. DOI: 10.1016/j.apacoust.2019.06.016
32. Khiem, N.T., Huyen, N.N., Long, N.T. Vibration of Cracked Timoshenko Beam Made of Functionally Graded Material. In: Harvie J., Baqersad J. (eds) Shock & Vibration, Aircraft/Aerospace, Energy Harvesting, Acoustics & Optics. Conference Proceedings of the Society for Experimental Mechanics Series. Springer, Cham. 2017. 9. Pp. 133-143. DOI: 10.1007/978-3-319-54735-0_15
33. Lien, T.V., Duc, N.T., Khiem, N.T. Mode Shape Analysis of Multiple Cracked Functionally Graded Timoshenko Beams. Latin American Journal of Solids and Structures. 2017. 14(7). Pp. 1327-1344. DOI: 10.1590/1679-78253496
34. Chouiyakh, H. Azrar, L. Alnefaie, K., Akourri, O. Vibration and multi-crack identification of Timoshenko beams under moving mass using the differential quadrature method. International Journal of Mechanical Sciences. 2017. 120. Pp. 1-11. DOI: 10.1016/j.ijmecsci.2016.11.014
35. Khiem, N.T., Huyen, N.N. A method for crack identification in functionally graded Timoshenko beam. Nondestructive Testing and Evaluation . 2017. 32. Pp. 319-341. DOI: 10.1080/10589759.2016.1226304
36. Pestel, E.C., Leckie, F.A. Matrix methods in elasto mechanics, McGraw-Hill, London, 1963. 435 p.
37. Attar, M. A transfer matrix method for free vibration analysis and crack identification of stepped beams with multiple edge cracks and different boundary conditions. International Journal of Mechanical Sciences. 2012. 57(1). Pp. 19-33. DOI: 10.1016/j.ijmecsci.2012.01.010
38. Khiem, N.T., Lien, T.V. The dynamic stiffness matrix method in forced vibration analysis of multiple-cracked beam. Journal of Sound and Vibration. 2002. 254(3). Pp. 541-555. DOI: 10.1006/jsvi.2001.4109
39. Khiem, N.T., Lien, T.V. A simplified method for natural frequency analysis of a multiple cracked beam. Journal of Sound and Vibration. 2001 245(4). Pp. 737-751. DOI: 10.1006/jsvi.2001.3585
40. Bozyigit, B., Bozyigit, I., Yesilce, Y., Abdel Wahab, M. Crack Identification in Multi-Span Beams on Elastic Foundation by Using Transfer Matrix Method, 2020. Proceedings of the 13th International Conference on Damage Assessment of Structures. Lecture Notes in Mechanical Engineering. Springer, Singapore.
41. Tan, G., Liu, Y., Gong, Y., Shen, Y., Liu, Z. Free Vibration of the Cracked Non-uniform Beam with Cross Section Varying as Polynomial Functions. KSCE Journal of Civil Engineering. 2018. 22(1). Pp. 4530-4546. DOI: 10.1007/s12205-018-1833-5
42. Kisa, M. Vibration and stability of multi-cracked beams under compressive axial loading. International Journal of Physical Sciences. 2011. 6(11). Pp. 2681-2696. DOI: 10.5897/IJPS11.493
43. Nguyen, K.V. Mode shapes analysis of a cracked beam and its application for crack detection. Journal of Sound and Vibration. 2014. 333(3). Pp. 848-872. DOI: 10.1016 /j.jsv.2013.10.006
44. Zheng, D. Y., Kessissoglou, N.J. Free vibration analysis of a cracked beam by finite element method. Journal of Sound and Vibration. 2004. 273(3). Pp. 457-475. DOI: 10.1016 /S0022-460X(03)00504-2
45. Aydin, K. Vibratory characteristics of axially loaded Timoshenko beams with arbitrary number of cracks. Journal of Vibration and Acoustics. 2007. 129(3). Pp. 341-354. DOI: 10.1115/1.2731411
46. Weaver Jr.W., Timoshenko, S.P., Young, D.H. Timoshenko, S., Jr, W.W. Vibration Problems in Engineering. 5th edn. Wiley. New York, 1990. 624 p.
47. Cheng, F.Y. Matrix Analysis of Structural Dynamics: Applications and Earthquake Engineering. 1st edn. Marcel Dekker Inc., New York, 2000. 997 p.
48. Al Rjoub, Y.S., Hamad, A.G. Free Vibration of Axially Loaded Multi-Cracked Beams Using the Transfer 49. Matrix Method. International Journal of Acoustics and Vibration. 2019. 24(1). Pp. 119-138. DOI: 10.20855/ijav.2019.24.11274
Contacts:
Yousef Al Rjoub, ysalrjoub@just.edu.jo Azhar Hamad, co.azhar@yahoo.com
© Al. Rjoub, Y.S.,Hamad, A.G., 2020