Next Article in Journal
A Robust Distributed Algorithm for Solving the Economic Dispatch Problem with the Penetration of Renewables and Battery Systems
Previous Article in Journal
A Digital 4D Information System on the World Scale: Research Challenges, Approaches, and Preliminary Results
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Applied Loads on Free Vibrations of Functionally Graded Material Plate–Shell Panels

by
José Simões Moita
1,
Victor Franco Correia
2,* and
Cristóvão Mota Soares
1
1
IDMEC, Instituto Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais, 1049-001 Lisboa, Portugal
2
ENIDH—Escola Superior Náutica Infante D. Henrique, Av. Eng. Bonneville Franco, 2770-058 Paço de Arcos, Portugal
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(5), 1993; https://doi.org/10.3390/app14051993
Submission received: 29 December 2023 / Revised: 24 February 2024 / Accepted: 27 February 2024 / Published: 28 February 2024
(This article belongs to the Section Materials Science and Engineering)

Abstract

:
An analysis of the influence of applied loads on free vibrations of plate and shell panels made of functionally graded materials is analyzed in the present work. Formulations for the static analysis considering geometrically nonlinear behavior, as well as linear buckling and free vibrations analyses are considered. A calculation of the through-thickness stress distribution is also performed. The finite element model is based on a higher order shear deformation theory using a non-conforming flat triangular plate/shell element with three nodes, and eight degrees of freedom per node is used in the numerical implementation. The results obtained with this numerical model are presented, discussed, and compared with alternative solutions published by other authors in some benchmark applications.

1. Introduction

Plate/shell panels have applications in a wide range of engineering situations, such as aeronautics and aerospace structural systems, among many others. Structures made of composite materials have been widely used to satisfy high performance demands. In such structures, stress singularities may occur at the interface between two different materials. In functionally graded structures, the smooth and continuous variation of the properties from one surface to the other eliminates abrupt changes in the stress and displacement distributions. In addition, a ceramic phase with low thermal conductivity can resist a high thermal environment, while the metal phase is strong with a mechanical load.
Thousands of studies of plate/shell panels composite structures have been published in the last five decades. More recently, a great number of works have also been published in functionally graded material (FGM) structures. Here, some representative works involving plate/shell panels and continuous variation of the material properties through the thickness direction, used to validate the static nonlinear and linear buckling and free vibrations analyses, are mentioned [1,2,3,4,5,6,7,8,9,10,11,12,13]. The aim and the novelty of this paper is to study the influence of external applied loads, which give the panels a state of stress and a slightly different geometry in comparison with the initial panel on the free vibration analysis of the FGM plates–shells, considering the geometrically nonlinear behavior.
Woo and Meguid [1] addressed the large deflections of FGM plates/shells, Reddy and Arciniega [2,3] studied the free vibrations and geometrically nonlinear analysis of FGM shells, Kim et al. [4] presented a formulation for small strains and large displacements, and Zhao and Liew [5] performed a geometrically nonlinear analysis of FGM shells. Tran et al. [6] presented the stability equations of FGM plates, while Thai and Kim [7] and Natarajan et al. [8] addressed the response of FGM plates in bending and free vibrations. Yin et al. [9], Moita et al. [10,11], and Long et al. [12] presented models for a geometrically nonlinear analysis of FGM plates and shells. Phan et al. [13] studied the buckling of FGM plates and shells.
In this paper, the influence of the applied loads on the natural frequencies of functionally graded material plate/shell panels subjected to general loadings is analyzed. The formulation considers the following: geometrically nonlinear deformation with the through-thickness stress distribution calculations, linear buckling, and free vibration analyses. A finite element model based on a higher order shear deformation theory using a non-conforming flat triangular plate/shell element with three nodes with eight degrees of freedom per node was adopted to study the plate/shell structures. This finite element model was revealed to be very efficient in the analysis of the proposed structures—it is easy to model the geometry, there is a low number of total degrees of freedom, there is good accuracy of the results, and the model is fast.
The solutions of some illustrative examples were performed, and the results are presented and discussed.

2. Formulation of the P-FGM Model

An FGM is made by mixing two distinct isotropic material phases, for example, a ceramic and a metal. In this work, the material properties of an FGM plate/shell structure are assumed to change continuously throughout the thickness (Figure 1), according to the volume fraction of the constituent materials, given by the power law function (Bao and Wang [14]). In addition, the continuous variation of the material mixture is approximated using a certain number of virtual layers, k, throughout the thickness direction-layer approach. The volume fraction of the ceramic and metal phases for each virtual layer is defined according to the power law:
V c k = 0.5 + z ¯ h   p ;   V m k = 1.0 V c k
where z ¯ is the thickness coordinate of the mid-surface of each layer and p is the gradient index of the power law.
Once the volume fractions of the ceramic and the metal, V c k and V m k , are defined, Young’s modulus E, the shear modulus G , the Poisson ratio ν , or the mass density ρ of each virtual layer k of the FGM plate/shell can be determined:
E k = V c k E c + V m k E m ;   G k = V c k G c + V m k G m ν k = V c k   ν c + V m k   ν m ;   ρ k = V c k   ρ c + V m k   ρ m

2.1. Displacement Fields and Strain-Displacement Relations

The displacement field of the plate/shell domain is based on Reddy’s third-order shear deformation theory [10,11]:
u x , y , z , ґ = u 0 x , y , ґ z   θ y x , y , ґ + z 3 c 1 θ y x , y , ґ w 0 x v x , y , z , ґ = v 0 x , y , ґ + z   θ x x , y , ґ + z 3 c 1 θ x x , y , w 0 y w x , y , z , ґ = w 0 x , y , ґ        
where u 0 ,   v 0 ,   w 0 are displacements of a generic point in the middle plane referred to the local axes—x, y, z directions; θ x ,   θ y are the rotations of the normal to the middle plane about the x axis (clockwise) and y axis (anticlockwise); w 0 x ,   w 0 y are the slopes of the tangents of the deformed mid-surface in the x , y directions; c 1 = 4 / 3   h 2 , with h denoting the total thickness of the plate/shell structure; and ґ is the time.
For the geometrically nonlinear behavior analysis, the strain tensor is conveniently represented in terms of the linear and nonlinear parts of the strain tensor as follows:
ε =   ε L + ε N L
The linear strain components associated with the displacement fields defined above can be written in the following form:
ε L = ε m + z ε b + z 3   ε * b   ε s + z 2   ε * s
Therefore, we can write in the detailed form:
ε x x =   u 0     x z     θ y     x + z 3 c 1 θ y x 2 w 0 x 2 + 1 2 u 0 x 2 + v 0 x 2 + w 0 x 2 ε y y =   v 0   y + z   θ x   y + z 3 c 1 θ x y 2 w 0 y 2 + 1 2 u 0 y 2 + v 0 y 2 + w 0 y 2 γ x y =   u 0   y +   v 0   x + z     θ y   y +   θ x   x + z 3 c 1   θ y y θ x x 2 2 w 0 x   y + u 0 x u 0 y + v 0 x v 0 y + w 0 x w 0 y γ x z = θ y + w 0 x + z 2 c 2 θ y w 0 x γ y z = θ x + w 0 y + z 2 c 2 θ x w 0 y
where c2 = 4/h2.

2.2. Constitutive Relations for the FGM Structures

The stress–strain relations are defined based on the stress–strain relations in each layer k, which can be written as follows:
σ k = Q k     ε k
where σ k = { σ x σ y σ x y τ x z τ y z } T is the stress vector, ε k = ε x ε y γ x y γ x z γ y z T is the strain vector, and Q k is the elasticity matrix of the k layer:
Q k = Q 11 Q 12 0 0 0 Q 21 Q 22 0 0 0 0 0 Q 66 0 0 0 0 0 Q 44 0 0 0 0 0 Q 55 k    
Q 1 1 k = Q 22 k = E k 1 ν k 2 ; Q 1 2 k = Q 21 k = ν k   E k 1 ν k 2 ; Q 4 4 k = Q 5 5 k = Q 6 6 k = E k 2   ( 1 + ν ) = G k
By integrating through the thickness, the linear elastic constitutive equation for each virtual layer, k, can be written as follows:
σ ^ k = D ^ k   ε k = N M M * Q Q * k = A B D 0 0 B C E 0 0 D E G 0 0 0 0 0 A s C s 0 0 0 C s E s k ε m ε b ε * b ε s ε * s k
where σ ^ k   is the resultant forces and moments, D ^ k   is   the constitutive matrix, and the sub-matrices of D ^ are given by the following:
A , B , C , D , E , G k = Q ¯ k   z k 1 z k ( 1 , z , z 2 , z 3 , z 4 , z 6 )   d z A s , C s , E s k = Q ¯ k   z k 1 z k ( 1 , z 2 , z 4 )   d z

3. Finite Element Approach

A non-conforming triangular plate/shell finite element model having three nodes and eight degrees of freedom per node is used. The degrees of freedom at each node are the displacements, the slopes, and the rotations, d i = u 0   v 0   w 0     w 0   y     w 0   x   θ x   θ y   θ z i T .
The rotation θ z i is introduced to consider a fictitious stiffness coefficient K θ z to eliminate the problem of a singular stiffness matrix for general shape structures [15]. The finite element formulation is obtained by the superposition of a membrane element with a bending element, and the side 1–2 of the triangle is coincident with the axis x of the local system attached to the finite element, with the origin at node 1 (Figure 2).
The element displacement vectors at any point of the mid-surface in the local coordinates (x,y) are obtained in terms of the corresponding nodal variables through linear shape functions L i , and the transverse displacement and slopes are expressed in terms of the corresponding nodal variables through cubic shape functions N i j , [15], which are also given in terms of area coordinates [15].
The displacement field can be represented in matrix form as follows:
u = Z   ( i = 1 3 N i   d i ) = Z   N   a e ;   a e = d 1 d 2 d 3 T
where the matrix of the shape functions and the appropriate matrix Z having powers of z k are as follows:
N i = L i 0 0 0 0 0 0 0 0 L i 0 0 0 0 0 0 0 0 N i 1 N i 2 N i 3 0 0 0 0 0 N i 1   y 2 N i   y N i 3   y 0 0 0 0 0 N i 1   x N i 2   x N i 3   x 0 0 0 0 0 0 0 0 L i 0 0 0 0 0 0 0 0 L i 0
Z = 1 0 0 0 z 3 c 1 0 z + z 3 c 1 0 0 1 0 z 3 c 1 0 z z 3 c 1 0 0 0 0 1 0 0 0 0 0

4. Virtual Work Principle

4.1. Static Analysis

The governing equations of the nonlinear problem are obtained from the virtual work principle in conjugation with an updated Lagrangian formulation [16,17]. A reference configuration is associated with time t, and the updated configuration is associated with the current time t = t + Δ t .
k = 1 N t A e h k 1 h k δ ε k L t Q k ε k L t d z d A e t + t A e h k 1 h k δ ε k N L t σ k L t d z d A e t = Π e t t + t k = 1 N t A e h k 1 h k δ ε k L t σ k L t d z d A e t
where Π e t t + t is the element external virtual work.
The membrane, bending, and shear strains, as well the higher order bending and shear strains can be represented by the following:
ε m = B m a e ;   ε b = B b z a e ;   ε s = B s a e ;   ε b * = B * b z 3 a e ;   ε s * = B * s z 2 a e
where the   B m , B b , B * b , B s , and B * s are components of the strain-displacement matrix B and are given explicitly in Moita et al. [18].
The virtual work principle applied to a finite element can be written in the following matrix form:
t A e B T D ^ B a e d A e t + t A e G T τ ^ G   a e d A e t = A e N T p d A e t A e   B T σ ^ t   d A e t
From this equation, the element linear stiffness matrix K L e , the element geometric stiffness matrix K σ e , the element external load vector F e x t e , and the element internal force vector F i n t e can be obtained [18].
The element geometric stiffness matrix   K σ e is given by the following:
K σ e = t A e G T τ ^ G dA e t
G i = L i x 0 0 0 0 0 0 0 L i y 0 0 0 0 0 0 0 0 L i x 0 0 0 0 0 0 0 L i y 0 0 0 0 0 0 0 0 N i 1 x N i 2 x N i 3 x 0 0 0 0 0 N i 1 y N i 2 y N i 3 y 0 0 0
τ ^ = N x x N x y 0 0 0 0 N x y N y y 0 0 0 0 0 0 N x x N x y 0 0 0 0 N x y N y y 0 0 0 0 0 0 N x x N x y 0 0 0 0 N x y N y y
To solve for general structures, these matrices and vectors, initially computed in the local coordinate system attached to the element, are transformed to a global coordinate system in the usual way [15]. Afterwards, by adding the contributions of all the elements in the domain, the system equilibrium equations can be obtained as follows:
K L + K σ t + Δ t t + Δ t   i 1   Δ q   i = F e x t t + Δ t F int t + Δ t t + Δ t   i 1
Using the Newton–Raphson incremental-iterative method, the incremental equilibrium path is obtained, and in the case of snap-through occurrence, the automatic arc-length method is used [19,20] via Equation (21) written in the following form:
K L + K σ t + Δ t t + Δ t   i 1   Δ q   i = λ i 1 + Δ λ i F e x t 0 F int t + Δ t t + Δ t   i 1
and an additional equation is employed to constrain the length of a load step:
Δ q i T Δ q i = Δ l 2
where F e x t 0 is a fixed (reference) load vector, λ is a load factor, Δ λ is the incremental load factor within the load step, and Δ l is the arc length.
In Crisfield [20], it is stated that the current stiffness parameter is a very useful index to give some scalar measure of the degree of non-linearity. In its unscaled form, it effectively measures the ‘stiffness’ of the system. In the present work, when using the arc length for any case studied here, this method calculates a stiffness parameter (herein called CSTIF) at any load level, starting with the value of CSTIF = 1.0 for the panel unloaded (stress-free state). In the incremental-iterative process, when the CSTIF became less than 1.0, it means the stiffness was decreasing, and if the CSTIF was greater than 1.0, the stiffness was increasing.

4.2. Linear Buckling Analysis

For the linear mechanical buckling analysis, only one load increment was used, and Equation (21) becomes:
K L 0 + K σ 0 q = F e x t 0
The linear elastic buckling standard eigenvalue problem is then given by the following:
K L + λ c r K σ 0 q = 0 ,         F e x t c r = λ c r   F e x t 0

4.3. Free Vibration Analysis

The Hamilton principle is used to formulate the governing equations for vibration analysis:
ґ 1 ґ 2 k = 1 N A e   h k 1 h k   δ ε ¯ k L     Q ¯ k   ε ¯ k L   d z   d t A e A e   h k 1 h k   δ u ˙ T   ρ k   u ˙   d z   d t A e   d ґ = 0
ґ 1 ґ 2 k = 1 N   A e h k 1 h k δ a e T B T Q ¯   B   a e   d z   d A e   A δ a ˙ e T N T   ρ k h k 1 h k Z T Z   d z   N a ˙ e   d A e   d ґ = 0
where the element mass matrix Me is given by the following:
M e = A N T   k = 1 N ρ k h k 1 h k Z T Z   d z     N   d A e
After assembling all the elements, the natural frequencies and the respective mode shapes are obtained by solving the standard eigenvalue problem:
K L ω 2   M q = 0
The free vibrations of the loaded panels are obtained considering the corresponding stiffness matrix by the following equation:
K L + K σ t t ω 2   t t M q = 0

5. Applications

First, a few numerical applications are presented, which serve for model validation against alternative results published by other researchers in different analyses, and also allow us to calculate load-displacement paths, which will give the information about the normal and shear stresses, limit points, and inflection points. The later applications, which address the aim of this work, show the non-linear analysis by using Equations (22) and (23) to obtain the influence of the applied loading at several load levels on the free vibrations of the FGM plates and shells. The updated geometry and the tangential stiffness K L + K σ matrix at each load level is calculated, and the eigenvalue problem in Equation (30) is solved for that loading condition to obtain the natural frequencies.

5.1. Nonlinear Analysis of a Simply Supported and Clamped FGM Square Plate under Uniform Transverse Pressure

Simply supported (SS) and clamped (CL) FGM square plates subjected to a uniform pressure load are considered in this first numerical application. The geometry is defined by the side dimension a = 1.0 m and the thickness h = 0.01 m. Three gradient indexes are considered: p = 0.5, p = 1, and p = 2. The plate, discretized through the thickness using 20 virtual layers, is made of constituents zirconia and aluminum ( E c = 151 × 109 Pa, ν c = 0.3, ρc = 5700 kg/m3, E m =70 × 109 Pa, ν m = 0.3, ρm = 2702 kg/m3). An incremental-iterative loading is applied, and the load-displacement paths obtained for the center of the plate are presented in Figure 3. This analysis had been validated in Moita et al. [11], although for different geometry and material properties. It is observed that the displacements increase when the gradient index increases. This happens because E c > E m , and as p increases, the volume fraction of ceramic decreases.
For the last load increment, considering the simply supported plate, the through-the-thickness distributions of the total normal and transverse shear stresses at the Gauss point nearest to the center point of the plate are shown in Figure 4 and Figure 5 for two different gradient indexes. From these figures, the influence of the gradient index p in the distributions of the stresses can be observed because the volume fractions of the constituents are related with the gradient index, and as expected, there are no abrupt changes in the stresses, in contrast with the structures made of composite materials. Thus, the results obtained for displacements, stresses, linear buckling, and natural frequencies of the plate–shell panels are affected by the gradient index p.
In this work, the variation through the thickness of FGM is approximated considering a certain number of virtual layers (20 layers in the present work). All the quantities are obtained at the middle of each layer, including the stresses. The upper and lower surfaces correspond to z ¯ / h = 0.5 and z ¯ / h = −0.5, respectively, where the shear stresses are zero. However, as mentioned above, the first and last virtual layers are obtained at the middle of these layers, i.e., z ¯ / h = 0.475 and z ¯ / h = −0.475.

5.2. Nonlinear Analysis of Hinged and Clamped FGM Cylindrical Panels under a Center Point Load

Two FGM cylindrical shell panels (Figure 6), one with the straight sides simply supported (hinged) and the curved sides free and the other with all sides clamped, are analyzed in this subsection. The geometry is defined by the following: R = 2.54 m, L = 0.508 m, and subtended angle 2 θ = 0.2 rad, H = 0.0127 m, and it is modelled by 392 triangular finite elements. The constituents are ZrO2 and Ti–6Al–4V with the following mechanical properties, respectively: E c = 168.0 × 109 N/m2, ν c = 0.298, ρc = 5700 kg/m3, E m = 105.7 × 109 N/m2,   ν m = 0.298, ρm =4429 kg/m3. For this case of geometry and the boundary conditions of the cylindrical panel, the snap-through can occur. The results obtained using the present model for the transverse displacement paths are shown in Figure 7a,b and Figure 8a,b, respectively, for the hinged cylindrical panel and for the clamped cylindrical panel. These analyses have been validated in Moita et al. [11], although for different geometry and material properties. It is observed from the paths obtained for the hinged cylindrical panel that the snap-through does occur, while for the clamped FGM panels, the snap-through does not occur. Additionally, as expected, the displacements increase when the gradient index increases, considering that E c > E m . Considering the same cylindrical panel with both boundary conditions, but under radial pressure load, the corresponding load-displacement paths are shown in Figure 7b and Figure 8b.

5.3. Critical Loads of a Clamped FGM Square Plate under Uniaxial In-Plane Applied Load

The linear buckling loads are obtained first in order to know which values should be considered for the calculations of the free vibrations under uniaxial in-plane applied load.
For the purpose of comparison, a square (a × a) FGM plate with all the edges movable and clamped, made of a mixture of zirconia and aluminum (Al/ZrO2), is considered. The material properties for the zirconia and aluminum are, respectively, E c = 151.0 × 109 N/m2, ν c = 0.3 , E m = 70.0 × 109 N/m2, and ν m = 0.3 ; the side dimensions are a = 0.2 m and thickness H = 0.01 m. The FGM plates are under uniaxial compression, and the power law index-p ranges from full ceramic to metal. The critical buckling parameters, defined as λ c r = N x x   a 2 / D 0 , with D 0 = E m   h 3 / 12 ( 1 ν 2 ) , are shown in Table 1. From this table, an excellent agreement between the results (for the first mode) obtained is observed with the present model (PM) compared with those obtained by Nguyen-Xuan et al. [21] using the MITC4 model and by Pham et al. [13]. For benchmarking purposes, the eigenvalues of some higher modes are also presented. Again, as Ec > Em, there is a decrease in the buckling loads when the p-index increases.
For comparison purposes, a square FGM plate with all edges movable and clamped is considered, but now with H = 0.01 m and sides with lengths of a = 1.0 m. The material properties are the following: E c = 380.0 × 109 N/m2, ν c = 0.3 , E m = 70.0 × 109 N/m2, ν m = 0.3 . The FGM plates are under uniaxial compression, and the power law index-p ranges from full ceramic to metal. The critical buckling parameter is now defined as λ c r = N x x   a 2 / E c H 3 , and the results obtained with the present model (PM) are compared in Table 2 with the alternative solutions for the first mode, obtained analytically by Wu et al. [22] using the first-order shear deformation theory. Figure 9 shows the variation of the non-dimensional buckling loads with the gradient index p.

5.4. Free Vibrations of a Simply Supported FGM Square Plate

The free vibrations of a simply supported FGM square plate are now investigated. The material properties for the zirconia and aluminum are, respectively, the following: E c = 151.0 × 109 N/m2, ν c = 0.3 , ρc = 5700 kg/m3, and   E m = 70.0 × 109 N/m2, ν m = 0.3 , ρm = 2702 kg/m3. The side dimension is a = 0.05 m, and the thickness is H = 0.01 m. The non-dimensional frequencies ω ¯ = ω h ρ m / E m   obtained with the present model (PM) are presented in Table 3 and compared with the alternative solutions obtained by Tran et al. [6] using a NURBS-based isogeometric approach and an HSDT model. For pure ceramic and pure metal, the results are very similar. The discrepancies between the results for the FGM plates are attributed to the different homogenization considered by Tran et al. [6] following the Mori–Tanaka scheme, while the present model uses the rule of mixtures. The corresponding graphical representation of the non-dimensional frequencies obtained by the present model is presented in Figure 10 for vibration modes 1 and 2.
The next investigation is to achieve how the free vibrations of an FGM plate with all the edges clamped and a power law index p = 1.0 are affected by two types of loading: the in-plane loading and transverse loading. The material properties are as follows: E c = 380.0 × 109 N/m2, ν c = 0.3 , ρm = 3000 kg/m3, E m = 70.0 × 109 N/m2, ν m = 0.3 , ρm = 2702 kg/m3. The geometry is defined by a = 1.0 m and H = 0.01 m. The results obtained are given in Table 4a–c. First, the free vibrations are calculated for the unloaded plate, and afterwards, for the three cases of loading:
(a)
Uniaxial in-plane uniform load from px = 0 to a final px = 2.0, 4.0, and 10.0 MPa, and the geometrically nonlinear behavior is considered;
(b)
Transverse uniform load from pz = 0 to a final pz = 20.0, 40.0, and 60.0 kPa;
(c)
Uniaxial in-plane uniform load combined with transverse uniform load, ranging from zero to the previously defined magnitudes.
Here, p z is uniformly distributed pressure applied at the lower surface of the plate with the positive direction of z. p x   is pressure applied at the edge x = 0 with the positive direction of x axis, in the usual way. In this context, P c will be a concentrated load applied to the center of the plate aligned with the positive direction of axis z.
From Table 4 a–c, it is seen that the uniform compressive loading decreases the natural frequencies, while the uniform transverse pressure increases the natural frequencies. In addition, it is observed when both types of loads are applied, the first frequency increases, but less so than when only the transverse pressure is applied, and the other frequencies increase or decrease depending on the values of two types of loadings. In addition, it is shown that the size of the decreasing/increasing of frequency is directly connected with the increase of the loading.
By analyzing the stiffness parameter for the plates, it can be verified that for the simply supported plate with axial pressure, p x , the stiffness parameter decreases, and for a uniform transverse pressure, p z , the stiffness parameter increases. For the clamped plate with axial pressure, p x , the stiffness parameter decreases, and for a uniform transverse pressure, p z , the stiffness parameter increases.
As for the cases of center point load and transverse pressure load, in Figure 11a,b, the load-displacement paths for the case of axial pressure loads applied to clamped and simply supported plates are presented. For the case of the clamped plates, an imperfect initial plate is considered, which consists of, in general, an initial center displacement w c 0 = h / 100 .
From Figure 11a, it is observed that the minimum stiffness occurs for the load level of 138 MPa for the clamped plate, corresponding to the nonlinear buckling, and from Figure 11b, an inflection point near the load level of 35 MPa is observed.
Considering the gradient index p = 1.0, the variation of the fundamental frequency versus applied type load is shown in Figure 11c,d for the case of the clamped plate. The variations match accordingly with Table 3 and Table 4. Proceeding in the same way, the influence of magnitude and the type of applied loads on the fundamental frequency of simply supported plates with all the edges restrained in u and v displacements is shown in Figure 11e,f. The free vibration frequencies of the first and second modes are influenced by the magnitude and type of loading, as shown in Figure 11f, for the case of the simply supported plate.
From Figure 11c–f, it can be verified that the variation of the fundamental frequency of the simply supported and clamped plates decreases as the stiffness of plate decreases under increasing axial pressure. For the case of transversal uniform pressure, the fundamental frequency of the simply supported and clamped plates increases due to the increasing plate stiffness with the increase of transverse pressure.

5.5. Critical Loads of Hinged and Clamped FGM Cylindrical Panels under Uniaxial Compression

To validate the present model, a buckling analysis of a square FGM panel made of a mixture of aluminum and zirconia with the following properties is considered: E m = 70.0 × 109 N/m2, ν m = 0.3, E c = 151.0 × 109 N/m2, ν c = 0.3. The geometry of the panel is as follows: length L = 0.1 m, angle 2θ = 0.2 rad, radius R = 0.5 m, and thickness h = 0.1 m. The FGM panel is assumed to have the curved edges fixed and the straight edges simply supported (FSFS). The panel is under uniaxial compression applied on the curved edges. A 20 × 20 finite element mesh is used (800 triangular elements). The results of the first three buckling load parameters N ¯ c r = N c r R / E m h 2 are given in Table 5a for the different power law indexes considered. From Table 5a, a good agreement with the results obtained by Zhao and Liew [23] is observed for the first buckling load, but with discrepancies varying from 4.5% to 9.5% for the second and third buckling loads, respectively. These discrepancies can be related to the different models—HSDT and a classical finite element used by present model, and FSDT and a mesh-free method used by Zhao and Liew.
Hinged and clamped FGM cylindrical shell panels are now considered to obtain the critical loads. The geometry and mechanical properties are the same as in Section 5.2. The results for the non-dimensional buckling loads, λ c r = N y c r R / E m H 2 , obtained using the present model, are shown in Table 5 b,c, and it is observed that λ c r   decreases when the gradient index p increases. This happens because Ec > EFGM > Em. The corresponding graphic representation is shown in Figure 12 and Figure 13, respectively.

5.6. Free Vibrations of Hinged and Clamped FGM Cylindrical Panels

First, for comparison purposes, we utilized a clamped FGM cylindrical panel with the following geometry: R/10, a/h = 100, Ec = 380.0 × 109 N/m2, ν c = 0.3 , ρc = 3000 kg/m3, Em = 70.0 × 109 N/m2, ν m = 0.3 , ρc = 2707 kg/m3. The non-dimensional fundamental frequency is given by ω ¯ = ω a 2 ρ m / D m , with D m = E m H 3 / 12 ( 1 ν m 2 ) . The results obtained with the present model are presented and compared in Table 6, with good agreement, with the alternative solutions of Pradyumna and Bandyopadhyay [24] obtained using a higher-order formulation and a C0 finite element and the solutions published by Neves et al. [25] obtained using a higher-order shear deformation theory and radial basis functions. The graphical representation of the non-dimensional fundamental frequency versus the p-index is presented in Figure 14.
The natural vibrations of the hinged and clamped FGM cylindrical shells from Section 5.2, with the same geometry and material properties, are now analyzed. The first three natural frequencies are obtained for the initial geometry (unloaded panel) considering different gradient indexes p. Next, the cylindrical panel is loaded by three kinds of loads. The applied loads produced displacements—new geometry—and a stress distribution over the plate/shell structure, which depends on the gradient index p. The new geometry and the state of stresses, which influences the stiffness of the structure at each load level, both have an influence on the values of the natural frequencies. In Table 7 a,b, Table 8 a,b and Table 9 a,b, the first three natural frequencies are shown for different gradient indexes. The graphic representation of the influence of applied loads on the fundamental frequency is shown in Figure 14 for the gradient index p = 1. The loadings are defined in Figure 6.
  • Hinged FGM panel loaded:
    (a)
    By a uniform pressure on the curved side, from p y py = 0 to a final p y = 20.0; 40.0 MPa. For this case of loading, the displacement v at y = L must be constrained to v = 0;
    (b)
    By a center point load, from P c = 0 to a final P c = 20.0; 40.0 kN;
    (c)
    By a uniform pressure on the curved side, from p y = 0 to a final p y = 20.0; 40.0 MPa, combined with a center point load, from P c = 0 to a final P c = 20.0; 40.0 kN.
The geometrically non-linear behavior is considered, and an incremental-iterative process is used. The results obtained with the present model are shown in Table 7 a,b. From these tables, the following can be observed:
(a)
The uniform compressive load minimally increases the first two natural frequencies but minimally decreases the third natural frequency.
(b)
The transverse center point load decreases the first three natural frequencies.
(c)
When both types of loads are applied, the first two frequencies decrease, but less than when only the center point load is applied.
(d)
It is shown that the frequencies decrease/increase more with the increase of the loadings.
2.
Clamped FGM panel loaded:
(a)
By a uniform pressure on the curved side, from p y = 0 to a final p y = 25.0; 50.0 MPa;
(b)
By a center point load, from P c = 0 to a final P c = 20.0; 40.0 kN, or by an external radial pressure, form pr = 0 to a final pr = 2 MPa;
(c)
By a uniform pressure on the curved side, from p y = 0 to a final p y = 25.0; 50.0 kPa, combined with a center point load, from P c = 0 to a final P c = 20.0; 40.0 kN, or combined with an external radial pressure, from p r = 0 to a final p r = 1 MPa; 2 MPa.
The corresponding results are shown in Table 8 a,b and Table 9 a,b. From Table 8 a,b, the following is observed:
(a)
The uniform compressive load, the center point load, and the radial pressure load decrease the first three natural frequencies.
(b)
When any two types of loads are applied, the first three frequencies decrease much more.
(c)
It is shown that the frequencies decrease more with the increase of the loadings.
By comparing Table 9 a,b, it is observed that the external radial pressure decreases the natural frequencies more than the center point load and uniform compressive load, as well as their combinations.
Analyzing the stiffness parameter for the cylindrical panels, it can be verified that for the cylindrical hinged panel, subject to an axial pressure on the free curved edge, p y , the stiffness parameter decreases. For a central point load, P c , the stiffness parameter decreases, and for a uniform radial pressure, p r , the stiffness parameter decreases. For the clamped panel subject to an axial pressure on the curved edge, p y , the stiffness parameter decreases. For a central point load, P c , t h e s t i f f n e s s parameter decreases, and for a uniform radial pressure, p r ,   t h e s t i f f n e s s   parameter decreases.
Considering the gradient index p = 1.0, the variation of the fundamental frequency versus applied type load and boundary conditions is shown in Figure 15a–c for the case of hinged cylindrical panels, and Figure 16a–c for the case of clamped cylindrical panels.
For the case of the hinged panel subjected to a center point load or radial pressure, the stiffness decreases, thus the fundamental frequency decrease is very low when the center point load is very close to the limit load. Observing Figure 15c, the fundamental frequency of the hinged shell panel increases under the axial pressure load, although the stiffness decreases, which can be attributed to the effect of the deformed shape on the frequency. For the uniaxial pressure, linear buckling has been achieved at 279.6 MPa. However, using the nonlinear analysis, the maximum transverse displacement is obtained at 200 MPa. At this pressure, the panel has the highest natural frequency.
For the axial compression load applied to the clamped cylindrical panel, shown in Figure 16a, the stiffness decreases and the fundamental frequency decreases under axial pressure load until the inflection point. For the same clamped cylindrical panel subjected to a center point load or a radial pressure load, the stiffness decreases, thus the frequency decreases until the inflection point.

6. Conclusions

A finite element model for the free vibration analysis of functionally graded material plate/shell panels under mechanical loading is presented. The loading is achieved by an incremental-iterative process to consider the geometrically nonlinear deformation of the structures.
The finite element model is based on Reddy’s third-order shear deformation theory applied to a non-conforming flat triangular plate/shell element with three nodes and eight degrees of freedom per node.
From the present applications, it is observed that the FGM plate/shell structures have a smooth variation of stresses through the thickness, and the free vibrations are influenced by the applied loads.
For the clamped and simply supported (with immovable edges) FGM plates, the uniaxial compressive pressure load leads to a decrease in the natural frequencies, while the uniform transverse pressure load leads to an increase in the frequencies. For the clamped FGM cylindrical panels, both types of loads lead to a decrease in the natural frequencies of the free vibration. In contrast, for the hinged FGM cylindrical panels, a uniform uniaxial compressive load leads to an increase in the first two frequencies, while the center point load and radial pressure load decrease the first three frequencies.
From these observations, it can be concluded that loads applied to plate/shell FGM structures can have significant impact on the natural frequencies, depending on the magnitude of the loads, the type of the applied loads, the boundary conditions, and the gradient index p.
In this work, only the gradient index p = 1 has been considered to build the numerical cases that show the influence of the applied loads on the fundamental frequency of plates and cylindrical panels (Figure 11c–f, Figure 15a–c and Figure 16a–c). It is observed that the stiffness is related to the frequency at each load level. In the case of the hinged cylindrical panel, the deformed shape seems to have an important role in influencing the frequency.
In future work, in order to develop a more comprehensive study, the effect on the variation of natural frequencies of higher modes (as shown in Figure 11f) with an arbitrary p-index could be studied. The present work provides results that can be used as a benchmark for interested readers to compare with new or available alternative theories in the development of new finite element models for FGM plate and shell structures with general geometry.

Author Contributions

Conceptualization, J.S.M.; methodology, J.S.M.; software, J.S.M.; validation, J.S.M.; formal analysis, all authors; investigation, all authors; resources, all authors; data curation, J.S.M.; writing—original draft preparation, J.S.M. and V.F.C.; writing—review and editing, V.F.C.; supervision, C.M.S.; project administration, C.M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by FCT, Fundação para a Ciência e Tecnologia, through IDMEC, under LAETA, project UIDB/50022/2020.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors on request.

Acknowledgments

The authors dedicate this paper to the memory of Professor Rolands Rikards (1942–2022), who sadly passed away on 2 February 2022. He was a world academic and distinguished researcher expert in the fields of characterization of material properties, multidisciplinary optimization, solid mechanics, finite element techniques, advanced composite materials, and computational methods. We extend our condolences to his wife and family, his colleagues at Riga Technical University, and his friends in Latvia.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Woo, J.; Meguid, S.A. Nonlinear analysis of functionally graded plates and shallow shells. Int. J. Solids Struct. 2001, 38, 7409–7421. [Google Scholar] [CrossRef]
  2. Reddy, J.N.; Arciniega, R.A. Free vibration analysis of functionally graded plates. In Analysis and Design of Plated Structures: Dynamics; Woodhead Publishing: Cambridge, UK, 2006. [Google Scholar]
  3. Arciniega, R.A.; Reddy, J.N. Large deformation analysis of functionally graded shells. Int. J. Solids Struct. 2007, 44, 2036–2052. [Google Scholar] [CrossRef]
  4. Kim, K.D.; Lomboy, G.R.; Han, S.C. Geometrically nonlinear analysis of functionally graded material (FGM) plates and shells using a four-node quasi-conforming shell element. J. Compos. Mater. 2008, 42, 485–511. [Google Scholar] [CrossRef]
  5. Zhao, Z.; Liew, K.M. Geometrically nonlinear analysis of functionally graded shells. Int. J. Mech. Sci. 2009, 51, 131–144. [Google Scholar] [CrossRef]
  6. Tran, L.V.; Ferreira, A.J.M.; Nguyen-Xuan, H. Isogeometric analysis of functionally graded plates using higher-order shear deformation theory. Compos. Part B Eng. 2013, 51, 368–383. [Google Scholar] [CrossRef]
  7. Thai, H.-T.; Kim, S.-E. A simple higher-order shear deformation theory for bending and free vibration analysis of functionally graded plates. Compos. Struct. 2013, 96, 165–173. [Google Scholar] [CrossRef]
  8. Natarajan, S.; Ferreira, A.J.M.; Bordas, S.; Carrera, E.; Cinefra, M.; Zenkour, A.M. Analysis of Functionally Graded Material Plates Using Triangular Elements with Cell-Based Smoothed Discrete Shear Gap Method. Math. Probl. Eng. 2014, 247932. [Google Scholar] [CrossRef]
  9. Yn, S.; Yu, T.; Bui, T.Q.; Nguyen, M.N. Geometrically nonlinear analysis of functionally graded plates using isogeometric analysis. Eng. Comput. 2015, 32, 519–558. [Google Scholar] [CrossRef]
  10. Moita, J.S.; Araújo, A.L.; Mota Soares, C.M.; Mota Soares, C.A.; Herskovits, J. Material and Geometric Nonlinear Analysis of Functionally Graded Plate-Shell Type Structures. Appl. Compos. Mater. 2016, 23, 537–554. [Google Scholar] [CrossRef]
  11. Moita, J.S.; Araújo, A.L.; Correia, V.F.; Mota Soares, C.M.; Herskovits, J. Higher-order finite element models for the static linear and nonlinear behaviour of functionally graded material plate-shell structures. Compos. Struct. 2019, 212, 465–475. [Google Scholar] [CrossRef]
  12. Long, N.V.; Tu, T.M.; Truong, H.K.; Hai, L.T.; Trang, V.T.T. Displacement-based and stress-based analytical approaches for nonlinear bending analysis of functionally graded porous plates resting on elastic substrate. Acta Mech. 2022, 233, 1689–1714. [Google Scholar] [CrossRef]
  13. Pham, V.V.; Nguyen, V.C.; Hadji, L.; Mohamed-Ouejdi, B.; Ömer, C. A comprehensive analysis of in-plane functionally graded plates using improved first-order mixed finite element model. Mech. Based Des. Struct. Mach. 2023. [Google Scholar] [CrossRef]
  14. Bao, G.; Wang, L. Multiple cracking in functionally graded ceramic/metal coatings. Int. J. Solids Struct. 1995, 32, 2853–2871. [Google Scholar] [CrossRef]
  15. Zienkiewicz, O.C. The Finite Element Method; McGraw-Hill: New York, NY, USA, 1977. [Google Scholar]
  16. Bathe, K.J.; Ho, L.W. A simple and effective element for analysis of general shell structures. Comput. Struct. 1981, 13, 673–681. [Google Scholar] [CrossRef]
  17. Bathe, K.J. Finite Element Procedures in Engineering Analysis; Prentice-Hall Inc.: Englewood Cliffs, NJ, USA, 1982. [Google Scholar]
  18. Moita, J.S.; Araújo, A.L.; Martins, P.G.; Mota Soares, C.M.; Mota Soares, C.A. Analysis of active-passive plate structures using a simple and efficient finite element model. Mech. Adv. Mater. Struct. 2011, 18, 159–169. [Google Scholar] [CrossRef]
  19. Crisfield, M.A. A fast incremental/iterative solution procedure that handles snap-through. Comput. Struct. 1980, 62, 13–55. [Google Scholar]
  20. Crisfield, M.A. Non-Linear Finite Element Analysis of Solid and Structures, Volume 1: Essentials; John Wiley & Sons: Chichester, UK, 1991. [Google Scholar]
  21. Nguyen-Xuan, H.; Tran, L.V.; Thai, C.H.; Nguyen-Thoi, T. Analysis of functionally graded plates by an efficient finite element method with node-based strain smoothing. Thin-Walled Struct. 2012, 54, 1–18. [Google Scholar] [CrossRef]
  22. Wu, T.-L.; Shukla, K.K.; Huang, J.H. Post-buckling analysis of functionally graded rectangular plates. Compos. Struct. 2007, 81, 1–10. [Google Scholar] [CrossRef]
  23. Zhao, X.; Liew, K.M. A mesh-free method for analysis of the thermal and mechanical buckling of functionally graded cylindrical shell panels. Comput. Mech. 2010, 45, 297–310. [Google Scholar] [CrossRef]
  24. Pradyumna, S.; Bandyopadhyay, J.N. Free vibration analysis of functionally graded curved panels using a higher-order finite element formulation. J. Sound Vib. 2008, 318, 176–192. [Google Scholar] [CrossRef]
  25. Neves, A.M.A.; Ferreira, J.M.; Carrera, E.; Cinefra, M.; Roque, C.M.C.; Jorge, R.M.N.; Soares, C.M.M. Free vibration analysis of functionally graded shells by a higher-order shear deformation theory and radial basis functions collocation, accounting for through-the-thickness deformations. Eur. J. Mech. A Solids 2013, 37, 24–34. [Google Scholar] [CrossRef]
Figure 1. Cross-section representation of an FGM plate.
Figure 1. Cross-section representation of an FGM plate.
Applsci 14 01993 g001
Figure 2. Finite element geometry.
Figure 2. Finite element geometry.
Applsci 14 01993 g002
Figure 3. Load-displacement paths for different gradient indexes.
Figure 3. Load-displacement paths for different gradient indexes.
Applsci 14 01993 g003
Figure 4. Distribution of total normal stress σx across the thickness.
Figure 4. Distribution of total normal stress σx across the thickness.
Applsci 14 01993 g004
Figure 5. Distribution of total transverse shear stress τxz across the thickness.
Figure 5. Distribution of total transverse shear stress τxz across the thickness.
Applsci 14 01993 g005
Figure 6. Cylindrical shell panel geometry and loadings.
Figure 6. Cylindrical shell panel geometry and loadings.
Applsci 14 01993 g006
Figure 7. (a) Load-displacement paths for the hinged cylindrical panel under point load. (b) Load-displacement paths for the hinged cylindrical panel under radial pressure load (p = 1).
Figure 7. (a) Load-displacement paths for the hinged cylindrical panel under point load. (b) Load-displacement paths for the hinged cylindrical panel under radial pressure load (p = 1).
Applsci 14 01993 g007
Figure 8. (a) Load-displacement paths for the clamped cylinder panel under center point load. Inflection point at 140 kN for p = 1. (b) Load-displacement paths for the clamped cylinder panel under radial pressure load. Inflection point at 2.25 MPa for p = 1.
Figure 8. (a) Load-displacement paths for the clamped cylinder panel under center point load. Inflection point at 140 kN for p = 1. (b) Load-displacement paths for the clamped cylinder panel under radial pressure load. Inflection point at 2.25 MPa for p = 1.
Applsci 14 01993 g008
Figure 9. Variation of the non-dimensional buckling loads versus the gradient index p for a clamped plate with immovable edges.
Figure 9. Variation of the non-dimensional buckling loads versus the gradient index p for a clamped plate with immovable edges.
Applsci 14 01993 g009
Figure 10. Variation of the non-dimensional frequency ω ¯ of the simply supported plate versus the gradient index p.
Figure 10. Variation of the non-dimensional frequency ω ¯ of the simply supported plate versus the gradient index p.
Applsci 14 01993 g010
Figure 11. (a) Load-displacement path for the clamped plate. (b) Load-displacement path for the SSSS plate. (c) Clamped plate under px. (d) Clamped plate under pz. (e) Simply supported plate under px. (f) Simply supported plate under pz.
Figure 11. (a) Load-displacement path for the clamped plate. (b) Load-displacement path for the SSSS plate. (c) Clamped plate under px. (d) Clamped plate under pz. (e) Simply supported plate under px. (f) Simply supported plate under pz.
Applsci 14 01993 g011
Figure 12. Variation of the non-dimensional buckling loads versus the gradient index p for a hinged cylindrical panel.
Figure 12. Variation of the non-dimensional buckling loads versus the gradient index p for a hinged cylindrical panel.
Applsci 14 01993 g012
Figure 13. Variation of the non-dimensional buckling loads versus the gradient index p for a clamped cylindrical panel.
Figure 13. Variation of the non-dimensional buckling loads versus the gradient index p for a clamped cylindrical panel.
Applsci 14 01993 g013
Figure 14. Variation of the non-dimensional frequency ω ¯ versus the gradient index p for a clamped cylindrical panel.
Figure 14. Variation of the non-dimensional frequency ω ¯ versus the gradient index p for a clamped cylindrical panel.
Applsci 14 01993 g014
Figure 15. (a) Hinged cylindrical panel. (b) Hinged cylindrical panel. Limit load: 89.92 KN. Limit load: 0.55 MPa. (c) Hinged cylindrical panel.
Figure 15. (a) Hinged cylindrical panel. (b) Hinged cylindrical panel. Limit load: 89.92 KN. Limit load: 0.55 MPa. (c) Hinged cylindrical panel.
Applsci 14 01993 g015
Figure 16. (a) Clamped cylindrical panel. (b) Clamped cylindrical panel. (c) Clamped cylindrical panel. Inflection point at 140 kN. Inflection point at 2.25 MPa.
Figure 16. (a) Clamped cylindrical panel. (b) Clamped cylindrical panel. (c) Clamped cylindrical panel. Inflection point at 140 kN. Inflection point at 2.25 MPa.
Applsci 14 01993 g016
Table 1. Non-dimensional buckling loads for plates with movable and clamped edges and different gradient indexes, with a/H = 20.
Table 1. Non-dimensional buckling loads for plates with movable and clamped edges and different gradient indexes, with a/H = 20.
Power Law Index
Mode SourceCeramicp = 0.2p = 0.5p = 1.0p = 2.0p = 5.0Metal
1 PM203.833181.190160.204142.854130.327120.06694.492
1 Ref. [21]203.497180.331159.529142.533130.342120.471-
1 Ref. [13]202.738-158.889141.962-119.773-
2 PM230.012204.622180.998161.502146.936135.029106.628
3 PM372.296331.653293.574261.502237.458217.276172.587
4 PM473.132421.549373.224332.490301.889276.107219.332
Table 2. Non-dimensional buckling loads for square plates with clamped edges and different gradient indexes, with a/H = 100.
Table 2. Non-dimensional buckling loads for square plates with clamped edges and different gradient indexes, with a/H = 100.
Power Law Index
Mode SourceCeramicp = 1.0p = 2.0p = 5.0Metal
1 PM 9.1874.6063.5713.0201.690
1 Ref. [22]9.1584.6183.5793.034-
2 PM10.6885.3424.1553.4931.969
3 PM18.0999.0477.0365.9083.334
Table 3. Non-dimensional frequencies ω ¯ for simply supported plate with different gradient indexes.
Table 3. Non-dimensional frequencies ω ¯ for simply supported plate with different gradient indexes.
ModeSourceCeramicp = 1.0p = 2.0Metal
1PM0.24620.22750.22540.2112
Ref. [6]0.24610.21850.21900.2113
2PM0.45390.43420.42370.3899
Ref. [6]0.45390.41180.40390.3897
3PM0.45390.43420.42370.3899
Ref. [6]0.45390.41180.40390.3897
4PM0.53790.50130.49300.4620
Ref. [6]0.53850.47940.47680.4623
5PM0.53790.50100.49300.4620
Ref. [6]0.53850.47940.47680.4623
Table 4. (a) Natural frequencies [Hz] for a fully clamped plate under different types of loading (p = 1). (b) Natural frequencies [Hz] for a fully clamped plate under different types of loading (p = 1). (c) Natural frequencies [Hz] for a fully clamped plate under different types of loading (p = 1).
Table 4. (a) Natural frequencies [Hz] for a fully clamped plate under different types of loading (p = 1). (b) Natural frequencies [Hz] for a fully clamped plate under different types of loading (p = 1). (c) Natural frequencies [Hz] for a fully clamped plate under different types of loading (p = 1).
(a)
ModePlate Unloadedpx = 2 MPapz = 20 kPapx = 2 MPa + pz = 20 kPa
1 142.325141.442144.252143.185
2291.271290.132292.259290.786
3292.066291.098293.134292.019
4432.917431.828433.395432.224
(b)
ModePlate Unloadedpx = 4 MPapz = 40 kPapx = 4 MPa + pz = 40 kPa
1 142.325140.552149.526147.391
2291.271288.854295.034291.357
3292.066290.042295.806293.672
4432.917430.741435.552432.328
(c)
ModePlate Unloadedpx = 10 MPapz = 60 kPapx = 10 MPa + pz = 60 kPa
1 142.325137.843157.010154.310
2291.271284.797299.130292.997
3292.066287.177299.768296.824
4432.917427.461438.402433.154
Table 5. (a) The critical buckling load N ¯ c r of an FGM cylindrical panel under uniaxial compression and FSFS boundary conditions. (b) Non-dimensional buckling loads for the hinged cylindrical panel with different gradient indexes. (c) Non-dimensional buckling loads for the clamped cylindrical panel with different gradient indexes.
Table 5. (a) The critical buckling load N ¯ c r of an FGM cylindrical panel under uniaxial compression and FSFS boundary conditions. (b) Non-dimensional buckling loads for the hinged cylindrical panel with different gradient indexes. (c) Non-dimensional buckling loads for the clamped cylindrical panel with different gradient indexes.
(a)
SourceModep = 0p = 0.5p = 1.0p = 2.0p = 5.0
Ref. [23]11.71951.37001.22291.10250.9949
PM1.73801.38841.23161.10650.9973
Ref. [23]21.84161.45751.30011.17961.0750
PM1.94801.53851.36691.23651.1247
Ref. [23]32.39131.88631.68371.40171.4018
PM2.23001.73101.53311.38361.2675
(b)
Power Law Index
ModeCeramicp = 0.5p = 1.0p = 5.0Metal
10.13300.11370.10600.09430.0836
20.15950.13690.12770.11320.1004
30.32942.96640.26800.24050.2073
(c)
Power Law Index
ModeCeramicp = 0.5p = 1.0p = 5.0Metal
1 2.06681.78091.66381.47791.3004
22.16791.86041.74041.56261.3640
33.45752.96642.77472.49192.1754
Table 6. Non-dimensional fundamental frequency ω ¯ for the clamped cylindrical panel with different gradient indexes.
Table 6. Non-dimensional fundamental frequency ω ¯ for the clamped cylindrical panel with different gradient indexes.
 Power Law Index
ModeSourceCeramicp = 0.5p = 1.0
1PM102.80786.28577.186
Ref. [24]102.92387.54577.077
Ref. [25]102.78785.47877.638
Table 7. (a) Natural frequencies [Hz] for hinged FGM cylindrical panel, considering different gradient indexes and loading types. (b) Natural frequencies [Hz] for hinged FGM cylindrical panel, considering different gradient indexes and loading types.
Table 7. (a) Natural frequencies [Hz] for hinged FGM cylindrical panel, considering different gradient indexes and loading types. (b) Natural frequencies [Hz] for hinged FGM cylindrical panel, considering different gradient indexes and loading types.
(a)
Index pFrequencyHinged
Unloaded
Hinged
py = 20.0 MPa
Hinged
Pc = 20 kN
Hinged
py = 20 MPa
+Pc = 20 kN
1313.274316.566295.679300.890
0.52351.139355.694340.132345.626
3479.608479.055465.256465.001
1305.210310.625286.167291.501
1.02343.077348.121331.984337.146
3473.809473.269458.445458.245
1297.077302.268275.787281.273
5.02335.446339.251321.954326.915
3471.027470.347454.304454.022
(b)
Index pFrequencyHinged
Unloaded
Hinged
py = 40.0 MPa
Hinged
Pc = 40 kN
Hinged
py = 40 MPa
+Pc = 40 kN
1313.274320.973276.674282.668
0.52351.139353.198335.139332.018
3479.608478.479449.294448.897
1305.210313.219261.988271.357
1.02343.077345.654316.605322.333
3473.809472.696440.668440.884
1297.077305.536249.127258.280
5.02335.446336.296305.175309.792
3471.027469.658434.684434.823
Table 8. (a) Natural frequencies [Hz] for clamped FGM cylindrical panel, considering different gradient indexes and loading types. (b)Natural frequencies [Hz] for clamped FGM cylindrical panel, considering different gradient indexes and loading types.
Table 8. (a) Natural frequencies [Hz] for clamped FGM cylindrical panel, considering different gradient indexes and loading types. (b)Natural frequencies [Hz] for clamped FGM cylindrical panel, considering different gradient indexes and loading types.
(a)
Index pFrequencyClamped
Unloaded
Clamped
py = 25 MPa
Clamped
Pc = 20 kN
Clamped
py = 25 MPa
+Pc = 20 kN
1530.354525.088517.085512.525
0.52910.949906.395902.049897.850
3950.677940.616944.500935.275
1523.455518.643509.337504.472
1.02899.626894.811890.297885.837
3938.582927.958932.117922.348
1515.247509.832500.414494.868
5.02892.631887.328882.877877.858
3929.081917.333922.329911.439
(b)
Index pFrequencyClamped
Unloaded
Clamped
py = 50 MPa
Clamped
Pc = 40 kN
Clamped
py = 50 MPa
+Pc = 40 Kn
1530.354521.489504.152493.762
0.52910.949901.861893.429884.175
3950.677930.478938.769919.380
1523.455513.874495.884484.824
1.02899.626890.029881.296871.642
3938.582927.243926.165905.783
1515.247509.832486.434473.780
5.02892.631881.968873.564862.640
3929.081905.471916.220893.425
Table 9. (a) Natural frequencies [Hz] for clamped FGM cylindrical panel, considering different gradient indexes and loading types. (b) Natural frequencies [Hz] for clamped FGM cylindrical panel, considering different gradient indexes and loading types.
Table 9. (a) Natural frequencies [Hz] for clamped FGM cylindrical panel, considering different gradient indexes and loading types. (b) Natural frequencies [Hz] for clamped FGM cylindrical panel, considering different gradient indexes and loading types.
(a)
Index pFrequencyClamped
Unloaded
Clamped
py = 25 MPa
Clamped
pr = 1 MPa
Clamped
py = 25 MPa
+pr = 1 MPa
1530.354525.746487.428483.450
0.52910.949906.135878.406875.414
3950.677940.051928.638921.602
1523.455518.375478.441474.085
1.02899.626894.541865.624862.355
3938.582927.362915.678908.076
1515.247509.531468.468463.246
5.02892.631886.980857.417853.489
3929.081916.673905.549896.673
(b)
Index pFrequencyClamped
Unloaded
Clamped
py = 50 MPa
Clamped
pr = 2 MPa
Clamped
py = 50 MPa
+pr = 2 MPa
1530.354521.489453.350442.269
0.52910.949901.862849.881842.178
3950.677930.474914.233895.260
1523.455513.874445.354432.959
1.02899.626890.029837.117828.118
3938.582917.243902.145880.988
1515.247504.466439.559423.800
5.02892.631881.968830.819819.155
3929.081905.471894.629869.139
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Moita, J.S.; Correia, V.F.; Soares, C.M. Influence of Applied Loads on Free Vibrations of Functionally Graded Material Plate–Shell Panels. Appl. Sci. 2024, 14, 1993. https://doi.org/10.3390/app14051993

AMA Style

Moita JS, Correia VF, Soares CM. Influence of Applied Loads on Free Vibrations of Functionally Graded Material Plate–Shell Panels. Applied Sciences. 2024; 14(5):1993. https://doi.org/10.3390/app14051993

Chicago/Turabian Style

Moita, José Simões, Victor Franco Correia, and Cristóvão Mota Soares. 2024. "Influence of Applied Loads on Free Vibrations of Functionally Graded Material Plate–Shell Panels" Applied Sciences 14, no. 5: 1993. https://doi.org/10.3390/app14051993

APA Style

Moita, J. S., Correia, V. F., & Soares, C. M. (2024). Influence of Applied Loads on Free Vibrations of Functionally Graded Material Plate–Shell Panels. Applied Sciences, 14(5), 1993. https://doi.org/10.3390/app14051993

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