Next Article in Journal
Fuzzy Clustering Methods with Rényi Relative Entropy and Cluster Size
Next Article in Special Issue
Large Deformation Problem of Bimodular Functionally-Graded Thin Circular Plates Subjected to Transversely Uniformly-Distributed Load: Perturbation Solution without Small-Rotation-Angle Assumption
Previous Article in Journal
Figures of Graph Partitioning by Counting, Sequence and Layer Matrices
Previous Article in Special Issue
Closed-Form Solution for Circular Membranes under In-Plane Radial Stretching or Compressing and Out-of-Plane Gas Pressure Loading
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Refined Theory for Bending Vibratory Analysis of Thick Functionally Graded Beams

1
Laboratoire Interdisciplinaire Carnot de Bourgogne, Université Bourgogne Franche-Comté (UTBM), CEDEX, 90010 Belfort, France
2
Laboratoire LMC, Université Bourgogne Franche-Comté (UTBM), UMR-CNRS 5060, CEDEX, 90010 Belfort, France
3
Centre Borelli, ENS-University of Paris-Saclay, 91190 Gif-sur-Yvette, France
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(12), 1422; https://doi.org/10.3390/math9121422
Submission received: 16 April 2021 / Revised: 4 June 2021 / Accepted: 8 June 2021 / Published: 18 June 2021
(This article belongs to the Special Issue Mathematical Problems in Materials Science)

Abstract

:
A refined beam theory that takes the thickness-stretching into account is presented in this study for the bending vibratory behavior analysis of thick functionally graded (FG) beams. In this theory, the number of unknowns is reduced to four instead of five in the other approaches. Transverse displacement is expressed through a hyperbolic function and subdivided into bending, shear, and thickness-stretching components. The number of unknowns is reduced, which involves a decrease in the number of the governing equation. The boundary conditions at the top and bottom FG beam faces are satisfied without any shear correction factor. According to a distribution law, effective characteristics of FG beam material change continuously in the thickness direction depending on the constituent’s volume proportion. Equations of motion are obtained from Hamilton’s principle and are solved by assuming the Navier’s solution type, for the case of a supported FG beam that is transversely loaded. The numerical results obtained are exposed and analyzed in detail to verify the validity of the current theory and prove the influence of the material composition, geometry, and shear deformation on the vibratory responses of FG beams, showing the impact of normal deformation on these responses which is neglected in most of the beam theories. The obtained results are compared with those predicted by other beam theories. It can be concluded that the present theory is not only accurate but also simple in predicting the bending and free vibration responses of FG beams.

1. Introduction

Functionally graded materials (FGMs) are new types of composites obtained by mixing ceramic and metallic constituents [1,2,3,4]. Material properties vary continuously through the beam-thickness in function of the mixing proportion. This avoids the stress concentration observed in laminate composites. FGMs are reserved for specific uses, for example, coatings of thermal barriers for turbine blades, shielding for military applications, automotive, space and aerospace industries, biomedical materials.
In the past decades, studied on functionally graded material (FGM) beams [5,6,7,8,9,10], plates [11,12,13,14,15,16,17,18,19] and shells [20,21,22] have received substantial attention, and an extensive spectrum of beam and plate theories has been introduced, based on the classical and the shear deformation theories of beams and plates.
FGMs are currently in great demand by industries, requiring very specific models to analyze their behavior and predict their responses. Many researchers have been interested in different FGM structures analyses because of their wide application areas. Both main beam models, the Euler–Bernoulli model (CBT) for thin beams and the Timoshenko model (SDT) for thick beams, were introduced. The CBT model ignores the transverse shear deformation effect. It was modified to take into account the shear deformation into consideration, resulting in the SDT model. Nevertheless, this second model requires a shear correction to satisfy the top and bottom beams faces’ boundary conditions, which influences the results. The higher-order SDT aims to eliminate the failure of CBT and the first-order SDT by assuming a higher-order variation through FG beam thickness for transverse displacement without providing any shear correction.
Multiple models with various shear stress shapes have been proposed, e.g., the Reddy model [23]. Thai and Vo [24] have presented several refined theories of HSDT beams. They have shown that these models are very effective in the static and dynamic studies of FG and laminates beams. Recently, Ebrahimi et al. [25] analyzed thermo-mechanical vibration of temperature-dependent FGM Beams with porosities by using an HSDT. Aydogdu et al. [26] used the Euler–Bernoulli model and parabolic and exponential shear functions to examine the bending vibration-responses of a simply supported FG Beams. Ben Oumrane et al. [7] used different beam theories to investigate an FG thick Timoshenko-beam’s static behavior. A numerical solution for (TBT) and (HSDT) is presented by Simsek [27] using the Ritz method. The finite element method and HSDT are used by Rakesh et al. [28] to analyze the bending vibration responses of the thick FG beams. An analytical solution for the cantilevered thick FG beams is provided by Zhong and Yu [29] for various types of mechanical loads. Based on the neutral surface concept, Ould Larbi et al. [30,31] presented an efficient theory to study the bending free vibration of thick FG beams. Similarly, a new first order of SDT theory is developed by Bouremana et al. [32], based on the position of the neutral surface for thick FG beams.
The thickness-stretching impact was introduced first in analyzing the vibrational behavior of thick FG plates [33]. Osofero et al. [34] provided an analysis method of buckling in bending of FG sandwich beams, considering thickness-stretching and shear effects. Meradjah et al. [35] also integrated the thickness-stretching effects in a new shear strain theory to analyze the bending vibration of FG beams.
From the literature mentioned above, it is evident that there is no published work considering the impact of thickness-stretching on the mechanical vibration of thick FG beams. This problem is not well-investigated and there is a need for further studies. In this work, a refined theory is presented to analyze the bending vibration of the thick FG beam, with supported ends and under transverse loading. This theory provides a constant transverse displacement and higher-order variation of axial displacement through the depth of the beam so that there is no need for any shear correction factors. By superimposing the deflection on the bending, shear and thickness-stretching parts, the governing equations are derived from Hamilton’s principle. The equations system obtained is solved by using Navier’s solutions. The material characteristics are presumed to change through the beam thickness following the law of dosing. Detailed mathematical formulations are provided, and example results are proposed to show the relevance of the present theory and to verify this accuracy. Our approach presents the advantage of using less variables in comparation to other theories and it also proves the effects of thickness-stretching and the influences of many parameters such as material index and slenderness ratio on frequency response and stresses. The proposed higher-order normal deformation and shear theory is not only accurate, but also provides an elegant and easy-to-implement approach to simulating the bending and vibration behaviors of thick FG beams.

2. Theoretical Formulation

2.1. Model Definition

Figure 1 shows the proposed model for this study. It is a thick FG beam with length ( L ) , rectangular cross-section, width ( b ) and height ( h ) .
The FG beam is composed of a combination of metal and ceramics whose combination changes from the top surface purely ceramic to the bottom surface completely metallic. The beam material effective characteristic P is assumed to changes through the FG beam thickness about the volume ratio and the characteristics of the constituent materials. It is formulated by the law of mixing as follow:
P P m ϑ m + P c ϑ c ,   P = ( E , ρ , ν   )
E ,   ρ ,   ν are Young’s modulus, mass density and Poisson’s coefficient, respectively. Variation of ν is generally small, so it remains constant ϑ c and ϑ m are ceramic and metal volume proportions respectively, defined by [36]:
ϑ m + ϑ c = 1 ,   ϑ c = ( 0.5 + z h ) p ,   p 0
The gradient index ( p ) , with p 0 determines the profile of the material in the FG beam thickness direction. It can be modified to obtain the optimum component materials distribution. The plot in Figure 2 shows the distribution of ceramic volume proportion across the FG beam thickness for various material indexes.
Each effective characteristic of the FG beam can be expressed as follows:
P ( z ) = ( P c P m ) ( 0.5 + z h ) p + P m ,   P = E , ρ

2.2. Displacement and Strain Fields

Transverse and axial displacements of the FG beam are expressed, according to the quasi-3D theory [35,37] as follow:
{ U ( x , z ) = u 0 ( x ) z w b / x 𝒻 ( z ) w s / x W ( x , z ) = w b ( x ) + w s ( x ) + w s t ( x )
with,
w s t ( x ) = ( z ) ϕ ( x ) ,   ( z ) = 1 d 𝒻 ( z ) / d z
  • u 0 : Axial displacement;
  • w b : Bending transverse displacement;
  • w s : Shear transverse displacement;
  • w s t : Thickness-stretching displacement;
  • u 0 , w b , w s and ϕ are four unknowns to be determined;
  • 𝒻 ( z ) and ( z ) are the shape functions.
The strains are as follows:
{ ε x = u 0 / x z 2 w b / x 2 𝒻 ( z ) 2 w s / x 2 ε z = ( d ( z ) / d z ) ϕ γ x z = ( z ) [ w s / x + ϕ / x ]
The FG beam material follows Hooke’s law. So, linear elastic equation can be expressed as:
{ σ x σ z τ x z } = E ( z ) [ 1 ν 0 ν 1 0 0 0 1 2 ( 1 + ν ) ] { ε x ε z γ x z } = E ( z ) { u 0 x z 2 w b x 2 𝒻 ( z ) 2 w s x 2 + ν d ( z ) d z ϕ ν u 0 x ν z 2 w b x 2 ν 𝒻 ( z ) 2 w s x 2 + d ( z ) d z ϕ 1 2 ( 1 + ν ) ( w b x + w s x ) }

2.3. Calculation of Energies

  • Strain energy:
    { δ U = 0 L 0 b [ h 2 h 2 ( σ x δ ε x + τ x z δ γ x z + σ z δ ε z ) d z ] d y d x = 0 L ( N x δ u 0 x x b 2 δ w b x 2 x s 2 δ w s x 2 + Q x z ( δ w s x + δ ϕ x ) + N z δ ϕ ) d x
    where, N x ,   x b ,   x s ,   Q x z and N z are the stress resultants, specified as:
    { ( N x , x b , x s ) = h / 2 h / 2 [ 1 , z , 𝒻 ( z ) ] σ x b d z Q x z = h / 2 h / 2 τ x z ( z ) b d z N z = h / 2 h / 2 σ z d ( z ) d z b d z
    By using Equation (7), the stress resultants given in Equation (9) can be expressed as:
    { N x M x b M x s Q x z N z } = [ A B B s X 0 B D D s Y 0 B s D s H s Y s 0 X Y Y s Z 0 0 0 0 0 A s ] { u 0 x 2 w b x 2 2 w s x 2 ( w s x + ϕ x ) }
    ( A ,   , D , s , D s , s , X , Y , Y s , Z , A s ) are the FG beam stiffness expressed, and are given by the Appendix A.
  • Potential energy due the external transverse load applied:
    δ V = 0 L 𝓆 ( x ) δ W d x = 0 L 𝓆 ( x ) δ ( w b + w s + w s t ) d x
    𝓆 ( x ) : External transverse loading.
  • Kinetic energy:
    { δ K = 0 L 0 b [ h / 2 h / 2 ρ ( z ) ( U ˙ δ U ˙ + W ˙ δ W ˙ ) d z ] d y d x = 0 L ( I 0 [ u ˙ 0 δ u ˙ 0 + ( w ˙ b + w ˙ s ) ( δ w ˙ b + δ w ˙ s ) ] + J 0 [ ( w ˙ b + w ˙ s ) δ ϕ ˙ + ϕ ˙ δ ( w ˙ b + w ˙ s ) ] ) d x 0 L ( I 1 [ u ˙ 0 δ w ˙ b x + w ˙ b x δ u ˙ 0 ] + I 2 w ˙ b x δ w ˙ b x + J 1 [ u ˙ 0 δ w ˙ s x + w ˙ s x δ u ˙ 0 ] ) d x 0 L ( K 2 w ˙ s x δ w ˙ s x J 2 [ w ˙ b x δ w ˙ s x + w ˙ s x δ w ˙ b x ] + K 0 ϕ ˙ δ ϕ ˙ ) d x
    where ( I 0 ,   I 1 , I 2 , J 0 , J 1 , J 2 , K 0 , K 2 ) are mass inertias, and are given by the Appendix A.

2.4. Governing Equation

To obtain the beam governing-equation, the Hamilton principle is applied as follows:
t 1 t 2 ( δ U + δ V δ K ) dt = 0
Substituting the expressions for δ U , δ V and δ K from Equations (8), (11) and (12) into Equation (13), the following system is obtained by integrating the parts and bringing together the coefficients of δ u 0 , δ w b , δ w s and δ ϕ :
{ δ u 0 : N x x = I 0 u ¨ 0 I 1 w ¨ b x J 1 w ¨ s x δ w b : 2 M x b x 2 + 𝓆 = I 0 ( w ¨ b + w ¨ s ) + I 1 u ¨ 0 x I 2 2 w ¨ b x 2 J 2 2 w ¨ s x 2 + J 0 ϕ ¨ δ w s : 2 M x s x 2 + Q x z x + 𝓆 = I 0 ( w ¨ b + w ¨ s ) + J 1 u ¨ 0 x J 2 2 w ¨ b x 2 K 2 2 w ¨ s x 2 + J 0 ϕ ¨ δ ϕ : Q xz x N z = J 0 ( w ¨ b + w ¨ s ) + K 0 ϕ ¨
Equation (14) can be expressed in terms of displacement of u 0 , w b ,   w s and ϕ by using Equation (10) as follows:
{ δ u 0 : A 2 u 0 x 2 B 3 w b x 3 B s 3 w s x 4 + X ϕ x = I 0 u ¨ 0 I 1 w ¨ b x J 1 w ¨ s x δ w b : B 3 u 0 x 3 D 4 w b x 4 D s 4 w s x 4 + Y 2 ϕ x 2 + 𝓆 = I 1 u ¨ 0 x + I 0 ( w ¨ b + w ¨ s ) I 2 2 w ¨ b x 2 J 2 2 w ¨ s x 2 + J 0 ϕ ¨ δ w s : B s 3 u 0 x 3 D s 4 w b x 4 H s 4 w s x 4 + A s 2 w s x 2 + ( A s + Y s ) 2 ϕ x 2 + 𝓆 = J 1 u ¨ 0 x + I 0 ( w ¨ b + w ¨ s ) J 2 w ¨ b x 2 K 2 w ¨ s x 2 + J 0 ϕ ¨ δ ϕ : X u 0 x Y 2 w b x 2 + ( A s + Y s ) 2 w s x 2 A s 2 ϕ x 2 + Z ϕ = J 0 ( w ¨ b + w ¨ s ) + K 0 ϕ ¨

2.5. Analytical Solution for a Simple Supported Functionally Graded Beam (S-S FG Beam)

In this section, the analytical solution is given for the case of Simply Supported Functionally Graded beam presented in Figure 3 (S-S FG beam) under the uniform load distributed over its entire length.
Analytical solutions of the motion equations are provided, based on Navier type solutions. The following displacements u 0 ,   w b ,   w s and ϕ are assumed to be combinations of known trigonometric functions which satisfy the boundary conditions and unknown coefficients to be determined for each value of n .
[ u 0 ( x , t ) w b ( x , t ) w s ( x , t ) ϕ ( x , t ) ] = n = 1 [ U n cos ( λ x ) e i ω n t W b n sin ( λ x ) e i ω n t W s n sin ( λ x ) e i ω n t ϕ n sin ( λ x ) e i ω n t ]
ω n is the eigenfrequency associated with the nth eigenmode, λ = n π / L , and U n ,   W b n ,   W s n and ϕ n are the unknown coefficients.
The following boundary conditions are imposed for a beam with two ends simply supported.
U = W = x b = x s = 0  
The assumed mechanical transverse load 𝓆 ( x ) is developed in a sinusoidal Fourier series as:
𝓆 ( x ) = n = 1 Q n s i n ( λ x )
The coefficients Q n are provided below for some loads.
  • Sinusoidal distribution case:
    n = 1 Q 1 = 𝓆 0  
  • Uniform distribution case:
    Q n = 4 q 0   n π   ,   ( n = 1 ,   3 ,   5 )
So, analytical solutions may be reached from the eigenvalues system below for any fixed value of n :
( [ K ] ω n 2 [ M ] ) { Δ } = { F }
In the static problem case, Equation (24) becomes:
[ K ] { Δ } = { F }
[ K ] = [ 𝓀 11 𝓀 12 𝓀 13 𝓀 14 𝓀 12 𝓀 22 𝓀 11 𝓀 24 𝓀 13 𝓀 23 𝓀 33 𝓀 34 𝓀 14 𝓀 24 𝓀 34 𝓀 44 ] , [ M ] = [ 𝓂 11 𝓂 12 𝓂 13 0 𝓂 12 𝓂 22 𝓂 23 𝓂 24 𝓂 13 𝓂 23 𝓂 33 𝓂 34 0 𝓂 24 𝓂 34 𝓂 44 ] , { Δ } = { U n W b n W s n ϕ n } , { F } = { 0 Q n Q n 0 }
( 𝓀 11 , 𝓀 12 , 𝓀 23 , 𝓀 22 , 𝓀 13 , 𝓀 24 , 𝓀 33 , 𝓀 14 , 𝓀 34 , 𝓀 44 ) and ( 𝓂 11 ,   𝓂 12 , 𝓂 13 , 𝓂 22 , 𝓂 23 , 𝓂 24 , 𝓂 33 , 𝓂 34 , 𝓂 44   ) are given by the Appendix A.

3. Numerical Results and Discussion

In this part, a uniform transverse load is applied to the S-S FG beam. Numerical examples are proposed first to validate the model presented above and assess its accuracy. The FG beam material is composed by A l 2 O 3   (alumina) and A l (aluminum). The material characteristics of the corresponding components are listed in Table 1.
The dimensionless form is used as follows:
U ¯ = 100 ( E m h 3 q 0 L 4 ) U ( 0 , h 2 ) ; W ¯ = 100 ( E m h 3 q 0 L 4 ) W ( L 2 , h 2 ) , σ x ¯ = ( h q 0 L ) σ x ( L 2 , h 2 ) ; σ x z ¯ = ( h q 0 L ) σ x z ( 0 , 0 ) ; σ z ¯ = ( h q 0 L ) σ z ( L 2 , h 2 ) ; ω ¯ = ( ω L 2 / h ) ρ m E m
A numerical example set out in Table 2 is performed for various material indexes (p) and slenderness ration (L/h) to validate the present model. The results obtained by this theory concerning displacements and the stresses for (L/h = 5, 20) are compared with those of the analytical solution provided by Li et al. [39]. The following shape function based on Reddy beam theory is used:
f ( z ) = 4 z 3 3 h 2

3.1. Static Analysis

Table 2 reveals that the present theory is with good agreement compared to the results obtained by Li et al. [39] and thus confirm the validation of the proposed method. It can also be seen that the CBT model, which omits shear deformation effects under-estimates displacements and stresses of the thick FG beams.
The impact of the material index (p) and slenderness ratio (L/h) on transverse W ¯ and axial U ¯ displacements, axial σ ¯ x and shear σ ¯ x z stresses of S-S FG beam under uniform load are illustrated in Table 2 and Table 3. In general, all shear deformation beam models give almost identical results, except for the case of transverse shear stress σ ¯ x z . It can be explained by the different transverse shear strain shape functions g(z) used in each models. It can be seen that those both displacements and axial stress of S-S FGB increases with increasing the power-law exponent (p). This is due to the fact that higher values of power law index (p) correspond to high portion of metal in comparison with the ceramic part, thus makes such FG beams more flexible. In fact, when p = 0, beam is made from fully ceramic and has the smallest displacements. By increasing the exponent (p) of the power law, the constitution of the FGB changes from an all-metal beam to a beam mixing ceramic and metal, then to an all-metallic beam (p→∞). Therefore, the percentage of metal phase increases and the beam becomes more flexible, which leads to increased displacements. However, the effect of the material index on the shear stress is negligible since the value of this stress is negligible compared to the axial stress.
To investigate again the effects of the thickness stretching on displacements, a comparison between the non-dimensional displacements of the beam obtained from the present model with and without thickness stretching is made in Figure 4a,b for the transverse and axial displacements, respectively, for (L/h = 5) and (p = 5). The difference between the two curves is seen. It is large out in the middle of the beam and becomes zero at the ends for the transverse displacement and becomes zero out in the middle of the beam for the axial displacement, unlike for the axial displacement, which is large at the ends and vanishes out in the middle of the beam.
Effect of shear on the evolution of non-dimensional axial displacement along the beam on the upper and lower beam faces is evaluated from Figure 5a,b for (L/h = 5) and (p = 5). We notice that shear have a very slight effect on the axial displacement. It can also be seen that shear is maximum at the ends and zero at the mid-span of the beam.
In the second examination and analysis example, effects of shear and thickness stretching on the non-dimensional axial and transverse displacements of the S-S FG beam are evaluated in Figure 6, for (L/h = 5) and (p = 5). It is obvious from these figures that the shear effect is more important on the two displacements, and it is greatest than the effect of thickness stretching for the transverse displacement. It means that the inclusion of shear deformation effect leads to an increase in the deflections and more pronounced for short beams. So, the shear effect on the displacements cannot be neglected, especially for the thick beams.
In Figure 7a,b, the evolution of non-dimensional (a) in-plane displacement at the end U ¯ ( 0 ) and (b) transverse displacement W ¯ ( L / 2 ) at mid-span through the FG beam thickness under uniform load is presented for (L/h = 5) and (p = 5). A slight difference appears for this shortest beam. It is seen that the maximum displacement is at the bottom of the beam for axial displacement and the median plane for transverse displacement. This is due to the consideration of the thickness stretching ( ε z 0 ) .
Figure 8a illustrates the evolution of the non-dimensional transverse shear stress and Figure 8b shows the stress due to the thickness stretching versus the power-law index ( p ) for (L/h = 5). All curves display the material index dependence of the stresses. It is clear that stress due to the thickness stretching exhibits low but not negligible values compared with those of the transverse shear stress.
Variations of the non-dimensional axial stress σ x ¯ across the depth and versus non-dimensional length of the FGB are presented in Figure 9a,b for (L/h = 5) and (p = 5). The axial stress is zero at the end edges, but it is maximum at the middle of the beam, and the upper beam face is stretched. On the other hand, the lower face is compressed. Extension stress at the upper beam face is higher than the compressive stress at the lower face because at the upper face, the beam is ceramic-rich, whereas, at the lower face, it is metal-rich. It is observed that the neutral plan with zero axial stress is moved upwards relative to the middle position. This is due to the non-homogeneous material of the FG beam (p = 5).
Figure 10a,b plots the distributions of the non-dimensional shear stress σ x z ¯ through thethickness and versus the non-dimensional length of the FG beam, respectively, for (L/h = 5) and (p = 5). These figures reveal that the shear stress reaches its maximum value at the beam ends but with opposite signs. It is cancelled in the middle of the beam on the neutral plane and the lower and upper faces. As a result, the conditions of non-shearing on both the lower and upper faces of the beam are satisfied.
The thickness-stretching impact on the stresses is evaluated in Figure 11a,b for (L/h = 5) and (p = 5). It may be seen that σ z ¯ is equal to zero at both end edges. Its highest value is reached in the middle of the beam. The upper beam face is stretched. On the other hand, the lower face is compressed. Compressive stress at the FG beam’s lower face is higher than extension stress at the upper face.

3.2. Vibration Analysis

Table 4 summarizes dimensionless frequencies associated with the first, second and third mode shapes of the S-S FG beam for various (L/h) ratio and material parameter (p).
Again, the results obtained by the present model, when the thickness-stretching is neglected, correlate closely with these obtained by the HSDT solution [30]. HSDT model under-estimates the frequencies of thick FG beams. This is due to the stretching of de the beam thickness omitted by the HSDT formulation in the thick FG beams case. It is emphasized that in the HSDT [30] formulation, the unknowns number is greater than this provided by the present model.
The dimensionless frequency variation versus material parameter for an (L/h = 5) with and without taking a count of the thickness stretching is plotted in Figure 12a. The plot showed that increasing (p) leads to a decrease of the frequencies. The highest frequency is achieved for (p = 0, completely ceramic beam), and the lowest for (p→∞, completely metal beam). This can be explained by the fact that the increase in (p) leads to a decrease in the amount of ceramic in the mixture, which is replaced by the metal, resulting in the decrease of Young’s modulus that makes the beam more flexible. Additionally, in this kind of composition, the density of the metal is greater than the density of ceramic and the total mass of the FG beam increases by increasing the metal proportion in the beam, which makes it softer. Therefore, the frequency decreases by increasing the material index (p).
The offset between both curves shows difference between the frequencies of both models with and without thickness stretching is significant for higher modes and for small slenderness ratios (L/h) (sees also Table 4 and Table 5). This is due to the effects of normal deformation. This effect leads to an increase of the vibration frequencies, and this increase is amplified at higher vibration modes and for small slenderness ratios. It implies that the normal deformation beam model should be employed for a better prediction of the frequencies instead of the model which neglects the effects of transverse normal deformation. It is seen that the frequencies are under-estimated when thickness stretching is omitted.
The impact of the slenderness ratio (L/h) on the frequencies is shown in Figure 12b and Table 4 and Table 5. It is seen that an increase of slenderness ration (L/h) leads to an increase of frequencies, and this increase is amplified at higher vibration modes. The frequencies tend to increase when the beam becomes shorter (or thicker).

4. Conclusions

In this paper, a refined beam theory is performed for bending vibratory analysis of the thick FG beams, taking into consideration thickness-stretching. The transverse displacement is assumed to be the sum of three components, bending, shear and stretching of the thickness. This leads to reducing the unknown’s number of parameters, therefore the number of governing equations. According to a mixing law, the FG beam effective material characteristics are supposed to change continuously along the thickness direction depending on the volume proportion of the constituents. The governing equations are obtained from Hamilton’s principle and are solved by using Navier-solutions.
Both the analytical and numerical results obtained in this article agree with those obtained using other theories with more unknown parameters. Sensitivity analysis has been performed for the thickness-stretching, the material parameter and the beam-slenderness. For this investigation, the shortest beams exhibited the greatest thickness-stretching impact, and this needs to be taken into consideration in more physically realistic simulations. Both geometry and material parameters affect the vibrational responses of the FG beams. Finally, it is observed that the proposed higher order shear and normal deformation theory is not only accurate but also provides an elegant and easily implementable approach for simulating bending and vibration behaviors of FGM beams, of relevance for example in spacecraft thermo-structural design. The formulation lends itself particularly well to finite element simulations and also to other numerical methods employing symbolic computation for beam bending problems, which will be considered in the near future. Further applications of our refined FG beam are planned in another contribution to consolidate our approach.

Author Contributions

Conceptualization, Y.B. and N.L.; methodology, Y.B., N.L. and D.B.; software, Y.B. and N.L.; validation, N.L. and D.B.; formal analysis, Y.B.; investigation, Y.B., N.L. and D.B.; resources, Y.B. and N.L.; data curation, Y.B. and N.L.; writing—original draft preparation, Y.B. and D.B.; writing—review and editing, D.B.; visualization, Y.B. and N.L.; supervision, N.L. and D.B.; project administration, N.L. and D.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

( A : : D : s : D s : s ) = h / 2 h / 2 ( 1 , z , z 2 , 𝒻 ( z ) , z 𝒻 ( z ) , 𝒻 ( z ) 2 ) E ( z ) b d z
( X , Y , Y s , Z ) = h / 2 h / 2 ( ν , ν z , ν 𝒻 ( z ) , d ( z ) d z ) d ( z ) d z E ( z ) b d z
A s = h / 2 h / 2 ( z ) 2 E ( z ) 2 ( 1 + ν ) b d z
( I 0 , I 1 , I 2 ) = h / 2 h / 2 ( 1 , z , z 2 ) ρ ( z ) b d z
( J 0 , J 1 , J 2 ) = h / 2 h / 2 ( , 𝒻 ( z ) , z 𝒻 ( z ) ) ρ ( z ) b d z
( K 0 , K 2 ) = h / 2 h / 2 ( ( z ) 2 , 𝒻 ( z ) 2 ) ρ ( z ) b d z
𝓀 11 = A λ 2 ;   𝓀 12 = λ 3 ;   𝓀 23 = D 11 s λ 4
𝓀 22 = D 11 λ 4 ;   𝓀 13 = s λ 3 ;   𝓀 24 = Y λ 2
𝓀 33 = λ 4 + A s λ 2 ;   𝓀 14 = X λ ;   𝓀 34 = ( A s + Y s ) λ 2
𝓀 44 = Z + A s λ 2
𝓂 11 = I 0   ;   𝓂 12 = I 1 λ   ;   𝓂 13 = J 1 λ
𝓂 22 = I 0 + I 2 λ 2 ;   𝓂 23 = I 0 + J 2 λ 2   ;   𝓂 24 = J 0
𝓂 33 = I 0 + K 2 λ 2   ;   𝓂 34 = J 0   ;   𝓂 44 = K 0

References

  1. Koizumi, M. FGM activities in Japan. Compos. Part B Eng. 1997, 28, 1–4. [Google Scholar] [CrossRef]
  2. Moita, J.S.; Araújo, A.L.; Correia, V.F.; Soares, C.M.M. Mechanical and thermal buckling of functionally graded axisymmetric shells. Compos. Struct. 2021, 261. [Google Scholar] [CrossRef]
  3. Moleiro, F.; Madeira, J.F.A.; Carrera, E.; Reddy, J.N. Design optimization of functionally graded plates under thermo-mechanical loadings to minimize stress, deformation and mass. Compos. Struct. 2020, 245. [Google Scholar] [CrossRef]
  4. Gul, U.; Aydogdu, M.; Karacam, F. Dynamics of a functionally graded Timoshenko beam considering new spectrums. Compos. Struct. 2019, 207. [Google Scholar] [CrossRef]
  5. Li, X.F. A unified approach for analyzing static and dynamic behaviors of functionally graded Timoshenko and Euler–Bernoulli beams. J. Sound Vib. 2008, 318, 1210–1229. [Google Scholar] [CrossRef]
  6. Benatta, M.A.; Mechab, I.; Tounsi, A.; Bedia, E.A.A. Static analysis of functionally graded short beams including warping and shear deformation effects. Comput. Mater. Sci. 2008, 44, 765–773. [Google Scholar] [CrossRef]
  7. Ben-Oumrane, S.; Abedlouahed, T.; Ismail, M.; Mohamed Bachir, B.; Mustapha, M.; El Abbas, A.B. A theoretical analysis of flexional bending of Al/Al2O3 S-FGM thick beams. Comput. Mater. Sci. 2009, 44, 1344–1350. [Google Scholar] [CrossRef]
  8. Sina, S.A.; Navazi, H.M.; Haddadpour, H. An analytical method for free vibration analysis of functionally graded beams. Mater. Des. 2009, 30, 741–747. [Google Scholar] [CrossRef]
  9. Şimşek, M. Static analysis of a functionally graded beam under a uniformly distributed load by Ritz method. Int. J. Eng. Appl. Sci. 2009, 1, 1–11. [Google Scholar]
  10. Ghayesh, H.M. Vibration analysis of shear-deformable AFG imperfect beams. Compos. Struct. 2018, 200, 910–920. [Google Scholar] [CrossRef]
  11. Menaa, R.; Tounsi, A.; Mouaici, F.; Mechab, I.; Zidi, M.; Bedia, E.A.A. Analytical solutions for static shear correction factor of functionally graded rectangular beams. Mech. Adv. Mater. Struct. 2012, 19, 641–652. [Google Scholar] [CrossRef]
  12. Yaghoobi, H.; Yaghoobi, P. Buckling analysis of sandwich plates with FGM face sheets resting on elastic foundation with various boundary conditions: An analytical approach. Meccanica 2013, 48, 2019–2035. [Google Scholar] [CrossRef]
  13. Feldman, E.; Aboudi, J. Buckling analysis of functionally graded plates subjected to uniaxial loading. Compos. Struct. 1997, 38, 29–36. [Google Scholar] [CrossRef]
  14. Abrate, S. Functionally graded plates behave like homogeneous plates. Compos. Part B 2008, 39, 151–158. [Google Scholar] [CrossRef]
  15. Mahdavian, M. Buckling analysis of simply-supported functionally graded rectangular plates under non-uniform in-plane compressive loading. J. Solid Mech. 2009, 1, 213–225. [Google Scholar]
  16. Mohammadi, M.; Saidi, A.R.; Jomehzadeh, E. Levy solution for buckling analysis of functionally graded rectangular plates. Appl. Compos. Mater. 2010, 17, 81–93. [Google Scholar] [CrossRef]
  17. Della Croce, L.; Venini, P. Finite elements for functionally graded Reissner–Mindlin plates. Comput. Methods Appl. Mech. Eng. 2004, 193, 705–725. [Google Scholar] [CrossRef]
  18. Yang, J.; Liew, K.M.; Kitipornchai, S. Second-order statistics of the elastic buckling of functionally graded rectangular plates. Compos. Sci. Technol. 2005, 65, 1165–1175. [Google Scholar] [CrossRef]
  19. Zhao, X.; Lee, Y.Y.; Liew, K.M. Mechanical and thermal buckling analysis of functionally graded plates. Compos. Struct. 2009, 90, 161–171. [Google Scholar] [CrossRef]
  20. Loy, C.T.; Lam, K.Y.; Reddy, J.N. Vibration of functionally graded cylindrical shells. Int. J. Mech. Sci. 1999, 41, 309–324. [Google Scholar] [CrossRef]
  21. Pradhan, S.C.; Loy, C.T.; Lam, K.Y.; Reddy, J.N. Vibration characteristics of functionally graded cylindrical shells under various boundary conditions. Appl. Acoust. 2000, 61, 111–129. [Google Scholar] [CrossRef]
  22. Strozzi, M.; Pellicano, F. Nonlinear vibrations of functionally graded cylindrical shells. Thin Walled Struct. 2013, 67, 63–77. [Google Scholar] [CrossRef]
  23. Moleiro, F.; Correia, V.M.F.; Araújo, A.L.; Soares, C.M.M.; Ferreira, A.J.M.; Reddy, J.N. Deformations and stresses of multilayered plates with embedded functionally graded material layers using a layerwise mixed model. Compos. Part B Eng. 2019, 156. [Google Scholar] [CrossRef]
  24. Thai, H.T.; Vo, T.P. Bending and free vibration of functionally graded beams using various higher-order shear deformation beam theories. Int. J. Mech. Sci. 2012, 62. [Google Scholar] [CrossRef] [Green Version]
  25. Ebrahimi, F.; Jafari, A. A Higher-Order Thermomechanical Vibration Analysis of Temperature-Dependent FGM Beams with Porosities. J. Eng. 2016, 2016. [Google Scholar] [CrossRef] [Green Version]
  26. Aydogdu, M.; Taskin, V. Free vibration analysis of functionally graded beams with simply supported edges. Mater. Des. 2007, 28, 1651–1656. [Google Scholar] [CrossRef]
  27. Şimşek, M. Buckling of Timoshenko beams composed of two-dimensional functionally graded material (2D-FGM) having different boundary conditions. Compos. Struct. 2016, 149. [Google Scholar] [CrossRef]
  28. Patil, R.; Joladarashi, S.; Kadoli, R. Studies on free and forced vibration of functionally graded back plate with brake insulator of a disc brake system. Arch. Appl. Mech. 2020, 90. [Google Scholar] [CrossRef]
  29. Zhong, Z.; Yu, T. Analytical solution of a cantilever functionally graded beam. Compos. Sci. Technol. 2007, 67. [Google Scholar] [CrossRef]
  30. Larbi, L.O.; Kaci, A.; Houari, M.S.A.; Tounsi, A. An efficient shear deformation beam theory based on neutral surface position for bending and free vibration of functionally graded beams. Mech. Based Des. Struct. Mach. 2013, 41. [Google Scholar] [CrossRef]
  31. Larbi, L.O.; Hadji, L.; Meziane, M.A.A.; Bedia, E.A.A. An analytical solution for free vibration of functionally graded beam using a simple first-order shear deformation theory. Wind Struct. Int. J. 2018, 27. [Google Scholar] [CrossRef]
  32. Bouremana, M.; Houari, M.S.A.; Tounsi, A.; Kaci, A.; Bedia, E.A.A. A new first shear deformation beam theory based on neutral surface position for functionally graded beams. Steel Compos. Struct. 2013, 15. [Google Scholar] [CrossRef]
  33. Fekrar, A.; Houari, M.S.A.; Tounsi, A.; Mahmoud, S.R. A new five-unknown refined theory based on neutral surface position for bending analysis of exponential graded plates. Meccanica 2014, 49. [Google Scholar] [CrossRef]
  34. Osofero, A.I.; Vo, T.P.; Nguyen, T.K.; Lee, J. Analytical solution for vibration and buckling of functionally graded sandwich beams using various quasi-3D theories. J. Sandw. Struct. Mater. 2016, 18. [Google Scholar] [CrossRef] [Green Version]
  35. Meradjah, M.; Kaci, A.; Houari, M.S.A.; Tounsi, A.; Mahmoud, S.R. A new higher order shear and normal deformation theory for functionally graded beams. Steel Compos. Struct. 2015, 18. [Google Scholar] [CrossRef]
  36. Praveen, G.N.; Reddy, J.N. Nonlinear transient thermoelastic analysis of functionally graded ceramic-metal plates. Int. J. Solids Struct. 1998, 35. [Google Scholar] [CrossRef]
  37. Minh, P.P.; Manh, D.T.; Duc, N.D. Free vibration of cracked FGM plates with variable thickness resting on elastic foundations. Thin Walled Struct. 2021, 161. [Google Scholar] [CrossRef]
  38. Tornabene, F. Free vibration analysis of functionally graded conical, cylindrical shell and annular plate structures with a four-parameter power-law distribution. Comput. Methods Appl. Mech. Eng. 2009, 198. [Google Scholar] [CrossRef]
  39. Li, X.F.; Wang, B.L.; Han, J.C. A higher-order theory for static and dynamic analyses of functionally graded beams. Arch. Appl. Mech. 2010, 80. [Google Scholar] [CrossRef]
Figure 1. The geometry of the FG beam.
Figure 1. The geometry of the FG beam.
Mathematics 09 01422 g001
Figure 2. Ceramic volume proportion profile across the FG beam thickness, for different material indexes.
Figure 2. Ceramic volume proportion profile across the FG beam thickness, for different material indexes.
Mathematics 09 01422 g002
Figure 3. Simply Supported Functionally Graded beam (S-S FG beam) under uniform load.
Figure 3. Simply Supported Functionally Graded beam (S-S FG beam) under uniform load.
Mathematics 09 01422 g003
Figure 4. Thickness stretching effect on the non-dimensional transverse (a) and axial displacements (b).
Figure 4. Thickness stretching effect on the non-dimensional transverse (a) and axial displacements (b).
Mathematics 09 01422 g004
Figure 5. Shear effect on axial displacement a long of the FGB on: (a) lower beam face and (b) upper beam face.
Figure 5. Shear effect on axial displacement a long of the FGB on: (a) lower beam face and (b) upper beam face.
Mathematics 09 01422 g005
Figure 6. Shear and stretching effects on the transverse displacement.
Figure 6. Shear and stretching effects on the transverse displacement.
Mathematics 09 01422 g006
Figure 7. Non-dimensional: (a) axial displacement, U ¯ and (b) mid-span deflection W ¯ through the FG beam thickness.
Figure 7. Non-dimensional: (a) axial displacement, U ¯ and (b) mid-span deflection W ¯ through the FG beam thickness.
Mathematics 09 01422 g007
Figure 8. Non-dimensional: (a) Shear stress σ x z ¯ , (b) stress due to the thickness stretching σ z ¯ versus the material parameter.
Figure 8. Non-dimensional: (a) Shear stress σ x z ¯ , (b) stress due to the thickness stretching σ z ¯ versus the material parameter.
Mathematics 09 01422 g008
Figure 9. σ x ¯ distributions: (a) over the FG beam thickness, (b) versus the FG beam non-dimensional length.
Figure 9. σ x ¯ distributions: (a) over the FG beam thickness, (b) versus the FG beam non-dimensional length.
Mathematics 09 01422 g009
Figure 10. σ x z ¯ distributions: (a) over the FG beam thickness, (b) versus the FG beam non-dimensional length.
Figure 10. σ x z ¯ distributions: (a) over the FG beam thickness, (b) versus the FG beam non-dimensional length.
Mathematics 09 01422 g010
Figure 11. σ z ¯ distributions: (a) over the FG beam thickness; (b) versus the FG beam non-dimensional length.
Figure 11. σ z ¯ distributions: (a) over the FG beam thickness; (b) versus the FG beam non-dimensional length.
Mathematics 09 01422 g011
Figure 12. Effect of (a) the material parameter and (b) slenderness ratio (L/h) on the dimensionless fundamental frequency of the S-S FG beam.
Figure 12. Effect of (a) the material parameter and (b) slenderness ratio (L/h) on the dimensionless fundamental frequency of the S-S FG beam.
Mathematics 09 01422 g012
Table 1. Material characteristics of A l 2 O 3 and A l [38].
Table 1. Material characteristics of A l 2 O 3 and A l [38].
Components ν E   ( G P a ) ρ   ( k g / m 3 )
Ceramic (alumina A l 2 O 3 )0.33803960
Metal (aluminium A l )0.3702702
Table 2. Comparison of non-dimensional transverseand axial displacements, axial and shear stresses of S-S FG beam for various material indexes and slenderness ratio (L/h).
Table 2. Comparison of non-dimensional transverseand axial displacements, axial and shear stresses of S-S FG beam for various material indexes and slenderness ratio (L/h).
p Theory L / h = 5 L / h = 20
W ¯ U ¯ σ ¯ x σ ¯ x z W ¯ U ¯ σ ¯ x σ ¯ x z
0 Li et al. [39]3.16570.94023.80200.75002.89620.230615.01300.7500
CBT2.87830.92113.7500-2.87830.230315.0000-
Present3.16810.94063.79190.75032.89620.230615.01290.7429
0.5 Li et al. [39]4.82921.66034.99250.76764.46450.408719.70050.7676
CBT4.44011.63314.9206-4.44010.408319.6825-
Present4.82021.66534.98930.76744.46440.408719.70030.7599
1 Li et al. [39]6.25992.30455.88370.75005.80490.568623.20540.7500
CBT5.77462.27225.7959-5.77460.568023.1834-
Present6.24752.29035.87970.75035.80490.568523.20520.7429
5 Li et al. [39]9.78023.70898.10300.57908.81510.913331.81120.5790
CBT8.75083.64968.1329-8.75080.912431.7711-
Present9.77873.69558.10990.58678.81810.913431.81270.5998
10 Li et al. [39]10.89793.88609.70630.64369.68790.953638.13720.6436
CBT9.60723.80979.5228-9.60720.952438.0913-
Present10.88473.87809.70860.66459.69050.953638.13830.6572
Table 3. Effect of slenderness ratio (L/h) and material index (p) on non-dimensional transverse and axial displacements, axial and shear stresses of S-S FG beam.
Table 3. Effect of slenderness ratio (L/h) and material index (p) on non-dimensional transverse and axial displacements, axial and shear stresses of S-S FG beam.
p L / h = 2 L / h = 7 L / h = 12
W ¯ U ¯ σ ¯ x σ ¯ x z W ¯ U ¯ σ ¯ x σ ¯ x z W ¯ U ¯ σ ¯ x σ ¯ x z
0 4.72892.57281.65880.74523.04060.67035.29470.75042.94080.38889.02480.7507
0.5 6.94704.47402.18760.76254.66181.18656.95070.76744.52690.689011.84400.7677
1 8.90366.15512.58660.74526.05221.64928.18980.75045.88380.958313.95230.7507
5 15.66710.01373.68840.60099.35202.652111.26540.60708.97771.539819.14500.6073
10 18.14610.63284.38600.658110.3452.772713.49600.66489.88301.608222.94750.6652
Table 4. Comparison of the S-S FG beam frequencies ω i ¯ for some values of ( L / h ) and ( p ) .
Table 4. Comparison of the S-S FG beam frequencies ω i ¯ for some values of ( L / h ) and ( p ) .
L / h ModeTheory p
00.51510
51HSDT [30] ε z = 0 5.15304.41103.99003.40003.2810
Present ε z = 0 5.15274.41073.99043.40123.2816
ε z 0 5.15164.42304.01693.43103.2984
2HSDT [30] ε z = 0 17.884015.461014.012011.535011.0220
Present ε z = 0 17.881215.458814.010011.543111.0240
ε z 0 17.890015.505214.097811.634811.0785
3HSDT [30] ε z = 0 34.225029.849027.108021.699020.7530
Present ε z = 0 34.209729.838227.097921.715820.5561
ε z 0 34.297529.967027.281321.888420.6748
201HSDT [30] ε z = 0 5.46004.65104.20503.64803.5390
Present ε z = 0 5.46034.65114.20513.64853.5390
ε z 0 5.46024.66574.23513.68353.5595
2HSDT [30] ε z = 0 21.573018.396016.634014.373013.9260
Present ε z = 0 21.573218.396216.634414.374613.9263
ε z 0 21.571018.452016.751114.509414.0043
3HSDT [30] ε z = 0 47.594040.653036.769031.572030.5340
Present ε z = 0 47.593040.652636.767931.578030.5369
ε z 0 47.584140.770937.019231.864930.7005
Table 5. Effect of slenderness ratio (L/h) on non-dimensional t fundamental frequency of the S-S FG beam for various material indexes.
Table 5. Effect of slenderness ratio (L/h) on non-dimensional t fundamental frequency of the S-S FG beam for various material indexes.
L / h ModeTheory p
0 0.5 1 5 10
21Present ε z = 0 4.12293.58113.24882.63692.5066
ε z 0 4.12883.59393.26972.65782.5199
2Present ε z = 0 11.27559.95319.07767.03876.5825
ε z 0 11.360610.03119.15487.09866.6339
3Present ε z = 0 18.832916.755615.354011.704510.8540
ε z 0 19.059316.937015.504011.813310.9627
71Present ε z = 0 5.30514.53034.09723.52233.4072
ε z 0 5.30434.54364.12543.55453.4256
2Present ε z = 0 19.500916.762315.175212.746912.2452
ε z 0 19.498616.808215.272012.852712.3049
3Present ε z = 0 39.322834.042730.860925.312924.1409
ε z 0 39.351434.149531.055325.513124.2620
121Present ε z = 0 5.42024.61994.17723.61563.5045
ε z 0 5.41984.63424.20683.64983.5244
2Present ε z = 0 20.983617.935516.223213.899713.4318
ε z 0 20.979317.987116.333114.024613.5027
3Present ε z = 0 44.963138.579634.916629.511828.4027
ε z 0 44.954238.685235.142529.761828.5434
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Boutahar, Y.; Lebaal, N.; Bassir, D. A Refined Theory for Bending Vibratory Analysis of Thick Functionally Graded Beams. Mathematics 2021, 9, 1422. https://doi.org/10.3390/math9121422

AMA Style

Boutahar Y, Lebaal N, Bassir D. A Refined Theory for Bending Vibratory Analysis of Thick Functionally Graded Beams. Mathematics. 2021; 9(12):1422. https://doi.org/10.3390/math9121422

Chicago/Turabian Style

Boutahar, Youssef, Nadhir Lebaal, and David Bassir. 2021. "A Refined Theory for Bending Vibratory Analysis of Thick Functionally Graded Beams" Mathematics 9, no. 12: 1422. https://doi.org/10.3390/math9121422

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop