Abstract
In the present paper, the free and forced vibration of multiple cracked multi span continuous beams made of functionally graded material (FGM) is investigated by the dynamic stiffness method. First, there are constructed dynamic stiffness matrix and nodal load vector of multiple cracked FGM beam element using Timoshenko beam theory and massless double spring model of crack. Then, effect of cracks parameters on free and forced vibration of the FGM continuous beams is examined. The theoretical developments are validated by numerical examples. The obtained results provide an efficient method to analyze free and forced vibration of multiple cracked FGM framed structures and assessment of the behavior of damaged structures.
Keywords
FGM; Multiple cracks; Timoshenko beam; Dynamic stiffness method; Vibration analysis
1. Introduction
Functionally graded materials (FGMs) have proved to be more benefit in comparison with the laminate composites and found widespread application in the high-tech industries such as aerospace, automobile, electronics, optics, chemistry, biomedical engineering etc,... For vibration analysis of FGM structures, a number of methods have been proposed such as analytical (Nguyen, Vo, Nguyen, & Lee, 2016Nguyen, T.-K., Vo, T. P., Nguyen, B.-D., & Lee, J. (2016). An analytical solution for buckling and vibration analysis of functionally graded sandwich beams using a quasi-3D shear deformation theory. Composite Structures, 156, 238-252.; Sina, Navazi, & Haddadpour, 2009Sina, S., Navazi, H., & Haddadpour, H. (2009). An analytical method for free vibration analysis of functionally graded beams. Materials & Design, 30(3), 741-747.), Rayleigh-Ritz method (Pradhan & Chakraverty, 2013Pradhan, K., & Chakraverty, S. (2013). Free vibration of Euler and Timoshenko functionally graded beams by Rayleigh-Ritz method. Composites Part B: Engineering, 51, 175-184.), Galerkin (Sheng & Wang, 2018Sheng, G., & Wang, X. (2018). Nonlinear vibration of FG beams subjected to parametric and external excitations. European Journal of Mechanics-A/Solids, 71, 224-234.), combined Fourier series - Galerkin method (Zhu & Sankar, 2004Zhu, H., & Sankar, B. (2004). A combined Fourier series-Galerkin method for the analysis of functionally graded beams. Journal of applied Mechanics, 71(3), 421-424.), differential quadrature method (DQM) (Sundaramoorthy Rajasekaran, 2013aRajasekaran, S. (2013a). Free vibration of centrifugally stiffened axially functionally graded tapered Timoshenko beams using differential transformation and quadrature methods. Applied Mathematical Modelling, 37(6), 4440-4463.; Sınır, Çevik, & Sınır, 2018Sınır, S., Çevik, M., & Sınır, B. G. (2018). Nonlinear free and forced vibration analyses of axially functionally graded Euler-Bernoulli beams with non-uniform cross-section. Composites Part B: Engineering, 148, 123-131.; Tang, Lv, & Yang, 2019Tang, Y., Lv, X., & Yang, T. (2019). Bi-directional functionally graded beams: asymmetric modes and nonlinear free vibration. Composites Part B: Engineering, 156, 319-331.), perturbation method (Sınır et al., 2018), asymptotic development method (Cao, Gao, Yao, & Zhang, 2018Cao, D., Gao, Y., Yao, M., & Zhang, W. (2018). Free vibration of axially functionally graded beams using the asymptotic development method. Engineering Structures, 173, 442-448.), discrete singular convolution and Taylor series expansion method (Wang & Yuan, 2017). While the aforementioned methods are limited to apply for analysis of simple FGM beam structures, the Finite Element Method (FEM) developed in (Alshorbagy, Eltaher, & Mahmoud, 2011Alshorbagy, A. E., Eltaher, M., & Mahmoud, F. (2011). Free vibration characteristics of a functionally graded beam by finite element method. Applied Mathematical Modelling, 35(1), 412-425.; Eltaher, Abdelrahman, Al-Nabawy, Khater, & Mansour, 2014; Mashat, Carrera, Zenkour, Al Khateeb, & Filippi, 2014Mashat, D. S., Carrera, E., Zenkour, A. M., Al Khateeb, S. A., & Filippi, M. (2014). Free vibration of FGM layered beams by various theories and finite elements. Composites Part B: Engineering, 59, 269-278.; Naccache & El Fatmi, 2018Naccache, F., & El Fatmi, R. (2018). Numerical free vibration analysis of homogeneous or composite beam using a refined beam theory built on Saint Venant’s solution. Computers & Structures, 210, 102-121.) could be used for more complex structures such as that are composed from different beam elements. However, because the FEM is developed on the base of the static shape functions, it requires a large amount of elements to study the high frequency components of vibration. This limitation of FEM could be overcome by using the Dynamic Stiffness Method (DSM) that was developed in (J. Banerjee, 2003Banerjee, J. (2003). Free vibration of sandwich beams using the dynamic stiffness method. Computers & Structures, 81(18-19), 1915-1922.; Deng & Cheng, 2016Deng, H., & Cheng, W. (2016). Dynamic characteristics analysis of bi-directional functionally graded Timoshenko beams. Composite Structures, 141, 253-263.; S Rajasekaran, 2013b; Su & Banerjee, 2015; Su, Banerjee, & Cheung, 2013; Van Lien, Khiem, & Duc, 2016Van Lien, T., Khiem, N. T., & Duc, N. T. (2016). Free vibration analysis of functionally graded Timoshenko beam using dynamic stiffness method. Journal of Science and Technology in Civil Engineering (STCE)-NUCE, 10(5), 19-28.) with the frequency-dependent shape functions. The most important advantage of the DSM is that allows investigating vibration of arbitrary high frequency with the minimum amount of elements.
Recently, vibrations of FGM beams with cracks have attracted attention of researchers because of the serious hazard that could be induced by the crack presence. The basic concepts of fracture mechanics for FGM structures have been given in (Erdogan & Wu, 1997Erdogan, F., & Wu, B. (1997). The surface crack problem for a plate with functionally graded properties. Journal of applied Mechanics, 64(3), 449-456.; Jin & Batra, 1996Jin, Z.-H., & Batra, R. (1996). Some basic fracture mechanics concepts in functionally graded materials. Journal of the Mechanics and Physics of Solids, 44(8), 1221-1235.) that are essentials to develop the models of crack in FGM structures. Namely, the continuous stiffness model (Christides & Barr, 1984Christides, S., & Barr, A. (1984). One-dimensional theory of cracked Bernoulli-Euler beams. International Journal of Mechanical Sciences, 26(11-12), 639-648.; E. C. Yang, Zhao, & Li, 2015Yang, E. C., Zhao, X., & Li, Y. H. (2015). Free vibration analysis for cracked FGM beams by means of a continuous beam model. Shock and Vibration, 2015.) proposes to treat a cracked FGM beam as the beam with continuously varying bending stiffness along its span. Therefore, the governing equation of vibration of the cracked FGM beam are similar to those of nonuniform beam that increases the difficulty to study beams made of FGM. An alternative model of crack that shows to be more appropriate for vibration analysis of cracked FGM beams was proposed by the use of rotational spring of stifness calculated in accordance with the concepts of fracture mechanics of FGM. Usefulness of the spring model of crack was validated in (A. Banerjee, Panigrahi, & Pohit, 2016Banerjee, A., Panigrahi, B., & Pohit, G. (2016). Crack modelling and detection in Timoshenko FGM beam under transverse vibration using frequency contour and response surface model with GA. Nondestructive Testing and Evaluation, 31(2), 142-164.) by using the 3D finite element formulation for an FGM beam with a crack. Using this model Yang et al. (J. Yang & Chen, 2008; J. Yang, Chen, Xiang, & Jia, 2008) have calculated natural frequencies, mode shapes and also response of cracked FGM Euler-Bernoulli beams to axial and moving loads. The influence of open edge cracks, material properties, beam slenderness ratio and end supports on free vibration of FGM Timoshenko beam was studied by Ke et al. (Ke, Yang, Kitipornchai, & Xiang, 2009Ke, L.-L., Yang, J., Kitipornchai, S., & Xiang, Y. (2009). Flexural vibration and elastic buckling of a cracked Timoshenko beam made of functionally graded materials. Mechanics of Advanced Materials and Structures, 16(6), 488-502.). The dynamic deflections of cracked FGM beam on an elastic foundation under a transverse moving load were studied by Yan et al. (Yan, Kitipornchai, Yang, & He, 2011). The authors showed that the elastic foundation makes dynamic deflections of FGM beam more sensitive to cracks. The effects of the geometric and elastic characteristics on the natural frequencies and mode shapes of a cracked FGM beam rested on Winkler-Pasternak foundation were investigated by Matbuly et al. (Matbuly, Ragb, & Nassar, 2009Matbuly, M., Ragb, O., & Nassar, M. (2009). Natural frequencies of a functionally graded cracked beam using the differential quadrature method. Applied mathematics and computation, 215(6), 2307-2316.) using differential quadrature method. The wave propagation in a cracked FGM cantilever beam under an impact force was considered by Akbas (Akbaş, 2016Akbaş, Ş. D. (2016). Wave propagation in edge cracked functionally graded beams under impact force. Journal of Vibration and Control, 22(10), 2443-2457.). The frequency equation in the form of third-order determinant that significantly simplifies calculating natural frequencies for FGM Euler-Bernoulli beam containing an arbitrary number of open edge cracks was established by Aydin (Aydin, 2013Aydin, K. (2013). Free vibration of functionally graded beams with arbitrary number of surface cracks. European Journal of Mechanics-A/Solids, 42, 112-124.). The natural frequencies and mode shapes of multiple cracked FGM Timoshenko beam were investigated in (Lien, Duc, & Khiem, 2017Khiem, N. T., Huyen, N. N., & Long, N. T. (2017). Vibration of cracked Timoshenko beam made of functionally graded material. In Shock & Vibration, Aircraft/Aerospace, Energy Harvesting, Acoustics & Optics, Volume 9 (pp. 133-143): Springer.a), (Lien, Đuc, & Khiem, 2017b) using the rotational spring model of cracks and actual position of neutral plane (Eltaher, Alshorbagy, & Mahmoud, 2013Eltaher, M., Alshorbagy, A., & Mahmoud, F. (2013). Determination of neutral axis position and its effect on natural frequencies of functionally graded macro/nanobeams. Composite Structures, 99, 193-201.). The obtained results show that, natural frequencies calculated with respect to neutral axis are over estimated than those calculated at the centroidal axis. Wei et al (Wei, Liu, & Xiang, 2012Wei, D., Liu, Y., & Xiang, Z. (2012). An analytical method for free vibration analysis of functionally graded beams with edge cracks. Journal of Sound and Vibration, 331(7), 1686-1700.) used spring model of crack to establish equations of motion of cracked FGM beam with rotary inertia and shear deformation included. Because of ignoring axial inertia, the bending vibration is uncoupled with axial one, so the model of cracks based only on rotational spring is inadequate for study of vibration in FGM beam where bending and longitudinal vibrations are usually coupled. To solve this problem, the two springs model of cracks were first instroduced by Sherafatnia et al. (Sherafatnia, Farrahi, & Faghidian, 2013Sherafatnia, K., Farrahi, G., & Faghidian, S. A. (2013). Analytic approach to free vibration and buckling analysis of functionally graded beams with edge cracks using four engineering beam theories. International Journal of Engineering-Transactions C: Aspects, 27(6), 979-990.), the cracks were modelled by two equivalent springs, one translational spring and one rotational spring with different stiffnesses evaluated based on the relations of the fracture mechanics. Using this model of cracks, the authors analyzed the natural frequencies and mode shapes of cracked FGM beam by different beam theories. Furthermore, Khiem et al. (N. Khiem & Huyen, 2017; N. T. Khiem, Huyen, & Long, 2017) also proposed two springs model of cracks which its stiffness calculated from crack depth. Unlikely to the previous studies, new formulas for determining crack magnitude of both the translational and rotational springs are proposed for arbitrary composition of FGM. These formulas allow one to use those developed in (Chondros, Dimarogonas, & Yao, 1998aChondros, T., Dimarogonas, A., & Yao, J. (1998a). A continuous cracked beam vibration theory. Journal of Sound and Vibration, 215(1), 17-34., 1998b) for cracked homogeneous beam. Using that crack model, the authors investigated the sensitivity of natural frequencies to crack of the Timoshenko FGM beam and proposed an exact procedure for identifying location and depth of single crack from given natural frequencies. The FEM has been also developed for modal analysis and identification of FGM Euler-Bernoulli beam with an open crack (Akbaş, 2013). The p-version FEM model was engaged by Yu and Chu (Yu & Chu, 2009Yu, Z., & Chu, F. (2009). Identification of crack in functionally graded material beams using the p-version of finite element method. Journal of Sound and Vibration, 325(1-2), 69-84.) to study change in frequencies and mode shapes of FGM beams due to crack location and size.
Although there is a widespread literature devoted to study vibration of FGM intact and cracked beam, most of them investigated simple beams except the Ref. (N. Khiem, Lien, & Ninh, 2018Khiem, N., Lien, T., & Ninh, V. (2018). Natural Frequencies of Multistep Functionally Graded Beam with Cracks. Iranian Journal of Science and Technology, Transactions of Mechanical Engineering, 1-36.) where multistep FGM beam with cracks was studied by the transfer matrix method. In this paper, the DSM is developed to study continuous FGM beam with arbitrary number of cracks based on the Timoshenko beam theory; spring model of cracks and power law distribution of FGM. The most important novelty of the present study is to use general solution for vibration shape of an FGM beam element with arbitrary number of cracks for constructing the dynamic stiffness matrix and load vector. This allows one to investigate crack continuous FGM beam with minimal number of elements equal to the number of spans in the beam. Moreover, the formulas proposed by Khiem et al. (N. Khiem & Huyen, 2017; N. T. Khiem et al., 2017) for calculating stiffness of springs representing cracks in dependence of varying material properties have been used with taking into account the actual position of neutral plane of FGM Timoshenko beam. The obtained results provide an efficient method for free and forced vibration analysis, as well as health monitoring of framed FGM structures with cracks.
2. Dynamic stiffness matrix and the nodal load vector of a multiple cracked FGM Timoshenko beam element
2.1. Governing equations
Consider a FGM beam of length L, cross sectional area A=b×h (Fig. 1) assuming that the material properties of FGM vary along the thickness direction by the power law distribution
where E, G, ρ and μ stand for Young’s, shear modulus, material density and Poisson’s ratio; subscripts t and b denote the top and bottom material; n is power law exponent; z is co-ordinate of point from the mid plane. In this study, we assumed that the value of Poisson’s ratios of the top and bottom material are constants. The beam is subject to axial and transverse loads. According to the Timoshenko beam theory, displacement field of the beam is expressed as
where are axial displacements, deflection of a point on neutral axis; is the distance between neutral axis and central axis x; θ is the rotation of the cross-section. Then the strain and stress field is obtained from displacement field using the following relations:
where κ is shear correction factor, κ =5/6 with rectangular section beam. Using the governing equations, kinetic and potential energy of the beam and work done by the external loads are calculated as
with the coefficients
Hence, using the Hamilton’s principle equations of motion in time domain are established as
with the physical boundary conditions
Using the Fourier transformation
the equations of motion in the frequency domain can be obtained in the form
where the following notations for matrices, vectors have been introduced
In case of free vibration equation (8) is reduced to
2.2. General solution of cracked FGM beam element
Seeking solutions of equation (9) in the form , one obtains characteristic equation
This is in fact a cubic algebraic equation with respect to with three roots η1,η2,η3 that allow for expressing solutions of above equation as (see Appendix 1 APPENDIX 1 A1. General solution of cub algebraic equation η 3 + a η 2 + b η + c = 0 where a = ω 2 [ I 11 A 33 + I 11 A 22 + I 22 A 11 A 11 A 22 ] ; b = ω 4 [ I 11 I 22 − I 12 2 A 11 A 22 + I 11 A 33 I 11 A 22 + I 22 A 11 A 11 A 22 ] − ω 2 I 11 A 22 c = ω 4 [ ω 2 I 11 A 33 I 11 I 22 − I 12 2 A 11 A 22 − I 11 2 A 11 A 22 ] Roots of cub algebraic equation are η1(ω),η2(ω),η3(ω) η 1 = − a / 3 + u − b 1 / u ; η 2,3 = − a / 3 − ( u − b 1 / u ) / 2 ± i 3 ( u + b 1 / u ) / 2 where u = ( a 1 + b 1 3 + c 1 2 − a 3 / 27 ) 1 / 3 ; a 1 = a b / 6 − c / 2 ; b 1 = b / 3 − a 2 / 9 ; c 1 = a 3 / 27 − a 1 A2. Nomenclature E, G, ρ Young’s, shear modulus and materila density n Power law exponent b, h, L width, height and lenght of the beam u ( x , t ) , θ ( x , t ) , w ( x , t ) Axial displacement, rotation and deflection U, Θ,W Amplitude of axial displacement, rotation and deflection z Vector z={U,Θ,W}T η 1 ,η 2, η 3 Three roots of characteristic equation N, Q, M Internal axial, shear forces and bending moment e 1 ,...,e n Crack positions γ j jth crack magnitudes G(x) Crack functions μj Recurrent coefficients z j Solution of Eq. (8) in subsequent interval (ej,ej+1) z0, zc General solution of Eq. (8) for intact and cracked beam ω, Ψ Natural frequency, mode shape )
Hence, general solution of the homogeneous equation (9) can be represented as
with being constants and
Based on (10) and (12) it is easy to verify that:
therefore, solution (11) can be rewritten in the form
where vector and matrix is
Assume, furthermore, that the beam has been cracked at position e and crack is modeled by two springs (N. Khiem & Huyen, 2017Khiem, N. T., Huyen, N. N., & Long, N. T. (2017). Vibration of cracked Timoshenko beam made of functionally graded material. In Shock & Vibration, Aircraft/Aerospace, Energy Harvesting, Acoustics & Optics, Volume 9 (pp. 133-143): Springer.), an axial spring of stiffness T and a rotational spring of stiffness R (Fig. 2). The conditions that must be satisfied at the crack are (Ke et al., 2009Ke, L.-L., Yang, J., Kitipornchai, S., & Xiang, Y. (2009). Flexural vibration and elastic buckling of a cracked Timoshenko beam made of functionally graded materials. Mechanics of Advanced Materials and Structures, 16(6), 488-502.)
So called crack magnitudes introduced in (15) can be calculated as functions of material parameters such as Young’s modulus, power exponent as
where (N. Khiem & Huyen, 2017Khiem, N. T., Huyen, N. N., & Long, N. T. (2017). Vibration of cracked Timoshenko beam made of functionally graded material. In Shock & Vibration, Aircraft/Aerospace, Energy Harvesting, Acoustics & Optics, Volume 9 (pp. 133-143): Springer.)
It is obviously that the crack magnitudes (16) include the case of homogeneous beam where
Now, seeking a particular solution of equation (9) that satisfied conditions
we obtain
where and [Gc(x)] is 3×3- matrix
with
and
Introducing the matrices
General solution of equation (9) satisfying conditions (15) at crack position can be rewritten as
By the manner that is employed to obtain solution (21) one is able to conduct solution of equation (9) for beam with single and double cracks as
So that solution of equation (9) for beam with n cracks can be expressed in the form
where so-called crack index vectors are calculated by the recurrent formulas
Substituting governing solution expression (13) into (22) - (23), we obtain
where
and
On the other hand, a particular solution of equation (8) can be found in the form
where [H(x,ω)] is matrix of transfer functions that satisfied equations
with boundary conditions
Hence, general solution of inhomogeneous equation (8) is composed of the general solution (24) and a particular solution (27)
2.3. Dynamic stiffness matrix and nodal load vector
Let’s consider an 2D-dimensional FGM beam element as shown in Fig. 3. Nodal displacements and forces of the element are introduced as
where
Substituting expression (30) into (32), we get
with [B F] is the matrix operator
Eliminating constant vector in equation (33) leads to
where and are respectively dynamic stiffness matrix and nodal load vector of the multiple cracked FGM beam element
For a given structure that consists of a number of FGM beam elements like above, by means of balancing all the internal forces at every nodes of the structure, there will be obtained total dynamic stiffness matrix and nodal load vector . Letting be the total DOF vector of the structure, equation of motion of which conducted by the dynamic stiffness method is
Therefore, natural frequencies can be found from of the equation
and mode shape related to the natural frequency ω j is
where
is an arbitrary constant and is the normalized solution of (38) corresponding to ω j.
In case of forced vibration, so-called frequency response in an element, say e-th, of the structure would be calculated by
3. Numerical results
3.1. Validation of the theoretical development
For validation of the proposed above theoretical development in this subsection, natural frequency ratios (cracked to intact), mode shapes of one-span FGM beam with different boundary conditions are computed and compared to those obtained in published studies.
Comparison of fundamental mode shape of a cracked FGM beam with different boundary conditions: (a) simple support, (b) clamped ends and (c) cantilever with (J. Yang & Chen, 2008Yang, J., Chen, Y., Xiang, Y., & Jia, X. (2008). Free and forced vibration of cracked inhomogeneous beams under an axial force and a moving load. Journal of Sound and Vibration, 312(1-2), 166-181.).
Figs. 4-5 show the variation of the first two frequency ratios (cracked/intact) of the FGM beam with an open edge crack of depth 20% along the crack position running from the left end to the right end of the beam. The results obtained in the present study are compared to those given by Yu & Chu (Yu & Chu, 2009Yu, Z., & Chu, F. (2009). Identification of crack in functionally graded material beams using the p-version of finite element method. Journal of Sound and Vibration, 325(1-2), 69-84.) (continuous lines) using p-version of the FEM with different boundary conditions as: a) Simple support, b) Clamped ends and c) Cantilever. Graphs given in the Figures demonstrate good agreement of developed above DSM with the improved FEM. The fundamental mode shape of cracked FGM beam given in Fig. 6 in comparison with that obtained by Yang & Chen in (J. Yang & Chen, 2008Yang, J., Chen, Y., Xiang, Y., & Jia, X. (2008). Free and forced vibration of cracked inhomogeneous beams under an axial force and a moving load. Journal of Sound and Vibration, 312(1-2), 166-181.) in different boundary conditions shows also excellent agreement of the DSM with the analytical method. This is a validation of not only the theory proposed above but also an advantage of the DSM developed for cracked FGM beam.
3.2. Free vibration of a multiple cracked FGM continuous beam
Consider a FGM multi-span continuous beam which has a rectangular section of b×h=0.1m×0.1m with material parameters: E t =70GPa; ρ t =2780kg/m 3 ; μ t =0.33; E b =350GPa; ρ b =7800kg/m 3 ; μ b =0.33 and power law index n=0.5 (Fig. 7).
Fig. 8 shows the variation of the first three frequency ratios (cracked to intact) of three-span FGM beam with a crack of depth 10%, 20%, 30%. Observing graphs given in the Figures allows one to make the following remarks:
-
a) Crack on different beam spans effects differently on natural frequencies: The first frequency is most sensitive to the crack appeared at the second (middle) span, while the second frequency is most sensitive to crack at the right span (near pinned end). The third frequency gets largest change for crack located near both the clamped and pinned ends of the beam.
-
b) There exist positions on continuous FGM beam, crack appeared at which makes no effect on certain natural frequency. Such positions are called here critical points for a given frequency. The critical point positions from the left end of the beam are 0.22m, 0.895m, 1.72m for the first frequency; 0.185m, 0.74m, 1.31m, 1.915m for the second frequency and 0.17m, 0.62m, 1.134m, 1.61m, 2.011m for the third one. The critical points do not depend upon the number of cracks that occurred in the beam and thay are a useful indication for localizing cracks if unchange of a certain frequency is recognized.
The variation of the fundamental frequency ratios of a FGM beam with various number of cracks from 1 to 10 on different beam spans
Fig. 9 shows the fundamental frequency ratios of a FGM beam varying with number of equidistant cracks from 1 to 10 in different cases of crack depths (10%, 20%, 30%). The graphs given in Figs. 9a-c show the frequency ratios for cracks on the first span, Figs. 9d-f and Figs. 9g-i for cracks on the second and third spans respectively. It is clear that increase of the number of cracks leads to more reduction of all the frequencies. It can be observed that variation of the frequencies is strongly dependent on which span the cracks are located.
The first two mode shapes of three-span continuous FGM beam in different number of cracks increasing from 1 to 4 on the second span.
Fig. 10a-b presents the first two mode shapes of a multi-span continuous FGM beam that has 1 to 4 equidistant cracks on the second span. The distances between cracks are 0.1m and the first crack is 0.1m far from the left node, crack depths are 30%. Fig 10c-d shows the deviation of the first two mode shapes between cracked beam to a corresponding intact one.
Deviation of the first two mode shapes of three-span continuous FGM beam with crack on the second span
Fig 11 shows the deviation of the first mode shape (Fig 11a-c) and the second mode shape (Fig 11d-f) of the beam with an open edge crack at positions 0.2m, 0.3m, 0.4m on the second span with depths 10%, 20%, 30% respectively. Observing graphs in the Figures allows one to make some remarks as follows:
-
The deviation of the mode shape has got sharp peak change at the crack position, but it is not maximum value. The increasing number of cracks leads to larger changes in the mode shapes of the beam, but largest change gets to be for cracked span. The deviation of mode shape is smooth along uncracked span. The mode shape deviation is should be be easily discriminated by using the wavelet analysis of mode shapes.
-
Likely to the natural frequencies, mode shapes of cracked continuous FGM beam has critical points where crack makes no change in certain mode shape like a simple span beam (Lien, Đuc, et al., 2017bLien, T. V., Đuc, N. T., & Khiem, N. T. (2017b). Mode Shape Analysis of Multiple Cracked Functionally Graded Timoshenko Beams. Latin American Journal of Solids and Structures, 14(7), 1327-1344.). For instance, crack at position of 0.2m from the left end of the second span does not change the first mode shape or crack at position 1.05m and 1.12m make no change in second mode shape (Fig. 11d). These mode shape nodes are useful in the problem of crack detection by measurement of mode shape.
3.3. Forced vibration of a multiple cracked FGM continuous beam
Consider an FGM continuous beam subjected to a harmonic distribution load q 0 =1000N/m and forced frequency ω as shown in Fig. 12.
Vibration amplitude of deflection of the beam with 1 to 4 equidistant cracks of depth 20% has been computed and presented in Fig. 13. Graphs given in Figs. 13a-c show the deflection diagrams and Figs. 13d-f show deviation of the beam deflection corresponding to the cases of cracked first, second and third spans respectively. Frequency of the load is 100 rad/s.
Vibration amplitude of deflection (a-c) and its deviation (d-f) of beam with 1 to 4 equidistant cracks
Amplitude diagrams of moment and shear force for three-span continuous FGM beam with 1 to 4 cracks on each span
Diagrams of moment and shear force are given in Fig. 14 for the beam with various number of cracks of equal depth 20%. Fig. 15 shows deflection diagrams of the beam in (a) various crack depth (10% - 40%) and (b) various load frequency (100rad/s - 400rad/s).
Vibration amplitude of deflection of FGM with various crack depth (a) and various load frequency (b)
Observing the Figs. 13 - 15 allows one to make the following remarks:
-
The vibration amplitude of deflection of cracked continuous FGM beam subjected to a harmonic distributed load gets a rough change at crack possitions, but the change is small and difficult to discriminate.
-
Increasing depth and number of cracks lead to larger change in vibration amplitude of deflection, so that can be easier to observe.
-
On the adjacent span, the variation of amplitude of deflection is continuous, uninterrupted, the variation amplitude can be higher at crack positions.
-
The vibration amplitude of deflection of the beam increases with growing load frequency so that can be properly chosen for crack localization.
-
All the chanes in deflection amplitude caused by cracks are useful indication to select a proper technique for crack detection by measurements of forced vibration amplitude.
4. Conclusions
In this paper, the dynamic stiffness method has been developed for vibration analysis of multiple cracked FGM continuous Timoshenko beams. Free and forced vibration of a three-span FGM continuous beam with multiple cracks have been examined to study effect of cracks and load on the beam natural frequencies, mode shapes and vibration amplitude. There are observed critical points for natural frequencies and mode shapes that are useful indicators in crack localization by measurements of the modal parameters. The crack-induced change in mode shapes and vibration amplitude of deflection demonstrates a possibility to propose an efficient procedure for crack detection by using wavelet analysis of the spatial vibration characteristics.
Acknowledgement
This research is funded by Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 107.02-2017.301.
References
- Akbaş, Ş. D. (2013). Free vibration characteristics of edge cracked functionally graded beams by using finite element method. International Journal of Engineering Trends and Technology, 4(10), 4590-4597.
- Akbaş, Ş. D. (2016). Wave propagation in edge cracked functionally graded beams under impact force. Journal of Vibration and Control, 22(10), 2443-2457.
- Alshorbagy, A. E., Eltaher, M., & Mahmoud, F. (2011). Free vibration characteristics of a functionally graded beam by finite element method. Applied Mathematical Modelling, 35(1), 412-425.
- Aydin, K. (2013). Free vibration of functionally graded beams with arbitrary number of surface cracks. European Journal of Mechanics-A/Solids, 42, 112-124.
- Banerjee, A., Panigrahi, B., & Pohit, G. (2016). Crack modelling and detection in Timoshenko FGM beam under transverse vibration using frequency contour and response surface model with GA. Nondestructive Testing and Evaluation, 31(2), 142-164.
- Banerjee, J. (2003). Free vibration of sandwich beams using the dynamic stiffness method. Computers & Structures, 81(18-19), 1915-1922.
- Cao, D., Gao, Y., Yao, M., & Zhang, W. (2018). Free vibration of axially functionally graded beams using the asymptotic development method. Engineering Structures, 173, 442-448.
- Chondros, T., Dimarogonas, A., & Yao, J. (1998a). A continuous cracked beam vibration theory. Journal of Sound and Vibration, 215(1), 17-34.
- Chondros, T., Dimarogonas, A., & Yao, J. (1998b). Longitudinal vibration of a continuous cracked bar. Engineering fracture mechanics, 61(5-6), 593-606.
- Christides, S., & Barr, A. (1984). One-dimensional theory of cracked Bernoulli-Euler beams. International Journal of Mechanical Sciences, 26(11-12), 639-648.
- Deng, H., & Cheng, W. (2016). Dynamic characteristics analysis of bi-directional functionally graded Timoshenko beams. Composite Structures, 141, 253-263.
- Eltaher, M., Abdelrahman, A., Al-Nabawy, A., Khater, M., & Mansour, A. (2014). Vibration of nonlinear graduation of nano-Timoshenko beam considering the neutral axis position. Applied mathematics and computation, 235, 512-529.
- Eltaher, M., Alshorbagy, A., & Mahmoud, F. (2013). Determination of neutral axis position and its effect on natural frequencies of functionally graded macro/nanobeams. Composite Structures, 99, 193-201.
- Erdogan, F., & Wu, B. (1997). The surface crack problem for a plate with functionally graded properties. Journal of applied Mechanics, 64(3), 449-456.
- Jin, Z.-H., & Batra, R. (1996). Some basic fracture mechanics concepts in functionally graded materials. Journal of the Mechanics and Physics of Solids, 44(8), 1221-1235.
- Ke, L.-L., Yang, J., Kitipornchai, S., & Xiang, Y. (2009). Flexural vibration and elastic buckling of a cracked Timoshenko beam made of functionally graded materials. Mechanics of Advanced Materials and Structures, 16(6), 488-502.
- Khiem, N., & Huyen, N. (2017). A method for crack identification in functionally graded Timoshenko beam. Nondestructive Testing and Evaluation, 32(3), 319-341.
- Khiem, N., Lien, T., & Ninh, V. (2018). Natural Frequencies of Multistep Functionally Graded Beam with Cracks. Iranian Journal of Science and Technology, Transactions of Mechanical Engineering, 1-36.
- Khiem, N. T., Huyen, N. N., & Long, N. T. (2017). Vibration of cracked Timoshenko beam made of functionally graded material. In Shock & Vibration, Aircraft/Aerospace, Energy Harvesting, Acoustics & Optics, Volume 9 (pp. 133-143): Springer.
- Lien, T. V., Duc, N. T., & Khiem, N. T. (2017a). Free Vibration Analysis of Multiple Cracked Functionally Graded Timoshenko Beams. Latin American Journal of Solids and Structures, 14(9), 1752-1766.
- Lien, T. V., Đuc, N. T., & Khiem, N. T. (2017b). Mode Shape Analysis of Multiple Cracked Functionally Graded Timoshenko Beams. Latin American Journal of Solids and Structures, 14(7), 1327-1344.
- Mashat, D. S., Carrera, E., Zenkour, A. M., Al Khateeb, S. A., & Filippi, M. (2014). Free vibration of FGM layered beams by various theories and finite elements. Composites Part B: Engineering, 59, 269-278.
- Matbuly, M., Ragb, O., & Nassar, M. (2009). Natural frequencies of a functionally graded cracked beam using the differential quadrature method. Applied mathematics and computation, 215(6), 2307-2316.
- Naccache, F., & El Fatmi, R. (2018). Numerical free vibration analysis of homogeneous or composite beam using a refined beam theory built on Saint Venant’s solution. Computers & Structures, 210, 102-121.
- Nguyen, T.-K., Vo, T. P., Nguyen, B.-D., & Lee, J. (2016). An analytical solution for buckling and vibration analysis of functionally graded sandwich beams using a quasi-3D shear deformation theory. Composite Structures, 156, 238-252.
- Pradhan, K., & Chakraverty, S. (2013). Free vibration of Euler and Timoshenko functionally graded beams by Rayleigh-Ritz method. Composites Part B: Engineering, 51, 175-184.
- Rajasekaran, S. (2013a). Free vibration of centrifugally stiffened axially functionally graded tapered Timoshenko beams using differential transformation and quadrature methods. Applied Mathematical Modelling, 37(6), 4440-4463.
- Rajasekaran, S. (2013b). Buckling and vibration of axially functionally graded nonuniform beams using differential transformation based dynamic stiffness approach. Meccanica, 48(5), 1053-1070.
- Sheng, G., & Wang, X. (2018). Nonlinear vibration of FG beams subjected to parametric and external excitations. European Journal of Mechanics-A/Solids, 71, 224-234.
- Sherafatnia, K., Farrahi, G., & Faghidian, S. A. (2013). Analytic approach to free vibration and buckling analysis of functionally graded beams with edge cracks using four engineering beam theories. International Journal of Engineering-Transactions C: Aspects, 27(6), 979-990.
- Sina, S., Navazi, H., & Haddadpour, H. (2009). An analytical method for free vibration analysis of functionally graded beams. Materials & Design, 30(3), 741-747.
- Sınır, S., Çevik, M., & Sınır, B. G. (2018). Nonlinear free and forced vibration analyses of axially functionally graded Euler-Bernoulli beams with non-uniform cross-section. Composites Part B: Engineering, 148, 123-131.
- Su, H., & Banerjee, J. (2015). Development of dynamic stiffness method for free vibration of functionally graded Timoshenko beams. Computers & Structures, 147, 107-116.
- Su, H., Banerjee, J., & Cheung, C. (2013). Dynamic stiffness formulation and free vibration analysis of functionally graded beams. Composite Structures, 106, 854-862.
- Tang, Y., Lv, X., & Yang, T. (2019). Bi-directional functionally graded beams: asymmetric modes and nonlinear free vibration. Composites Part B: Engineering, 156, 319-331.
- Van Lien, T., Khiem, N. T., & Duc, N. T. (2016). Free vibration analysis of functionally graded Timoshenko beam using dynamic stiffness method. Journal of Science and Technology in Civil Engineering (STCE)-NUCE, 10(5), 19-28.
- Wang, X., & Yuan, Z. (2017). Discrete singular convolution and Taylor series expansion method for free vibration analysis of beams and rectangular plates with free boundaries. International Journal of Mechanical Sciences, 122, 184-191.
- Wei, D., Liu, Y., & Xiang, Z. (2012). An analytical method for free vibration analysis of functionally graded beams with edge cracks. Journal of Sound and Vibration, 331(7), 1686-1700.
- Yan, T., Kitipornchai, S., Yang, J., & He, X. Q. (2011). Dynamic behaviour of edge-cracked shear deformable functionally graded beams on an elastic foundation under a moving load. Composite Structures, 93(11), 2992-3001.
- Yang, E. C., Zhao, X., & Li, Y. H. (2015). Free vibration analysis for cracked FGM beams by means of a continuous beam model. Shock and Vibration, 2015.
- Yang, J., & Chen, Y. (2008). Free vibration and buckling analyses of functionally graded beams with edge cracks. Composite Structures, 83(1), 48-60.
- Yang, J., Chen, Y., Xiang, Y., & Jia, X. (2008). Free and forced vibration of cracked inhomogeneous beams under an axial force and a moving load. Journal of Sound and Vibration, 312(1-2), 166-181.
- Yu, Z., & Chu, F. (2009). Identification of crack in functionally graded material beams using the p-version of finite element method. Journal of Sound and Vibration, 325(1-2), 69-84.
- Zhu, H., & Sankar, B. (2004). A combined Fourier series-Galerkin method for the analysis of functionally graded beams. Journal of applied Mechanics, 71(3), 421-424.
-
Available online December 10, 2018.
APPENDIX 1
A1. General solution of cub algebraic equation
where
Roots of cub algebraic equation are
where
A2. Nomenclature
- E, G, ρ Young’s, shear modulus and materila density
- n Power law exponent
- b, h, L width, height and lenght of the beam
- Axial displacement, rotation and deflection
- U, Θ,W Amplitude of axial displacement, rotation and deflection
- z Vector
- η 1 ,η 2, η 3 Three roots of characteristic equation
- N, Q, M Internal axial, shear forces and bending moment
- e 1 ,...,e n Crack positions
- γ j jth crack magnitudes
- G(x) Crack functions
- μj Recurrent coefficients
- z j Solution of Eq. (8) in subsequent interval
- z0, zc General solution of Eq. (8) for intact and cracked beam
- ω, Ψ Natural frequency, mode shape
Publication Dates
-
Publication in this collection
2019
History
-
Received
18 Aug 2018 -
Reviewed
08 Dec 2018 -
Accepted
08 Dec 2018