Next Article in Journal
Sonochemical Synthesis of Ce-doped TiO2 Nanostructure: A Visible-Light-Driven Photocatalyst for Degradation of Toluene and O-Xylene
Next Article in Special Issue
Bi-Axial Buckling of Laminated Composite Plates Including Cutout and Additional Mass
Previous Article in Journal
NDE Detection Techniques and Characterization of Aluminum Wires Embedded in Honeycomb Sandwich Composite Panels Using Terahertz Waves
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on the Buckling Behavior of Functionally Graded Plates with Stiffeners Based on the Third-Order Shear Deformation Theory

Modeling Evolutionary Algorithms Simulation and Artificial Intelligence, Faculty of Electrical & Electronics Engineering, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam
*
Authors to whom correspondence should be addressed.
Materials 2019, 12(8), 1262; https://doi.org/10.3390/ma12081262
Submission received: 19 March 2019 / Revised: 8 April 2019 / Accepted: 10 April 2019 / Published: 17 April 2019

Abstract

:
This paper presents a finite element formulation to study the mechanical buckling of stiffened functionally graded material (FGM) plates. The approach is based on a third-order shear deformation theory (TSDT) introduced by Guangyu Shi. The material properties of the plate were assumed to be varied in the thickness direction by a power law distribution, but the material of the stiffener was the same as that of the one of the bottom surface where the stiffener was placed. A parametric study was carried out to highlight the effect of material distribution, the thickness-to-width ratio, and stiffener parameters on the buckling characteristics of the stiffened FGM plates. Numerical results showed that the addition of stiffener to the FGM plate could significantly reduce the weight of the FGM plate but that both the FGM plates with and without stiffener had equally high strength in the same boundary condition and compression loading.

1. Introduction

The common way to achieve higher strength for functionally graded material (FGM) plates and shell without stiffeners is to either increase the thickness of this structure or to add stiffeners. The weight of the unstiffened structure will become higher with increasing thickness, but reinforcement with stiffener will reduce the weight as well as the cost of this structure. For this reason, using stiffeners is the best method in special cases such as ship building, bridge construction, aerospace, marine, and so on.
To give more useful information about the application in practice, the buckling behavior of composite structures has received much attention from scientists. Broekel and Gangadharaprusty [1] used experimental and theoretical solutions to study the mechanical responses of stiffened and unstiffened composite panels subjected to a uniform transverse loading. Liu and Wang [2] explored the elastic buckling of a plate reinforced by stiffener under in-plane loading. ANSYS modeling was applied to find the optimal height, number, and arrangement of stiffeners. Ueda et al. [3] proposed an analytical approach to research the buckling and deflection responses of a stiffened plate, which has a deflection under out-of-plane and in-plane loads. Danielson et al. [4] presented a combination of von Karman and nonlinear beam theories to predict the buckling behavior of a stiffened plate subjected to axial compression. Using the finite element method, Jiang and co-workers [5] found that the second-order 2D solid element gave accurate results for the buckling problem of unstiffened and stiffened plates. The ABAQUS solution was used by Hughes et al. [6] to examine the buckling behavior of T-stiffened panels subjected to uniaxial compression and lateral pressure. Pavlovcic et al. [7] investigated the buckling problem of imperfect stiffened panels using numerical simulation and tests.
Nowadays, with the development of science and technology, many fields require machine details and structures working in harsh environments, such as high temperature, large abrasion, complicated loads, and so on. Therefore, in order to satisfy this demand, the material industry must develop and find new materials specialized for practice areas. Thus far, many new materials have been created to meet this requirement, such as the typical functionally graded material invented by Japanese scientists in 1984. Due to its outstanding advantage of being able to perform in a high-temperature environment as well as its large load capacity, anti-radiation, and anti-corrosion properties, this material has been widely applied in many important areas, such as nuclear science, medicine, chemistry, mining, and so on.
Javaheri and Eslami [8] explored the thermal buckling behavior of rectangular FGM plates on the basis of classical plate theory and closed-form solutions. This approach was also used by Shariat and Eslami [9] to study the thermal buckling of imperfect FG plates. The thermal postbuckling response of FGM skew plates based on the finite formulation was adopted by Prakash et al. [10]. Recently, Moita et al. [11] presented a formulation for buckling and nonlinear analysis of FGM plates subjected to mechanical and thermal loadings. Based on a semianalytical approach, Dung and Nam [12] analyzed the nonlinear dynamics of eccentrically stiffened functionally graded circular cylindrical thin shells under external pressure and surrounded by an elastic medium. In [13], Yu and his co-workers investigated the thermal buckling of functionally graded plates with internal defects in which the extended isogeometric analysis was fully exploited. For this type of problem, Dung and Nga [14] also discovered the thermomechanical postbuckling of eccentrically stiffened sandwich plates on elastic foundations subjected to in-plane compressive loads, thermal loads, and thermomechanical loads at the same time. The buckling of parallel eccentrically stiffened functionally graded annular spherical segments were studied through the Donnell shell theory and smeared stiffeners technique by Nam et al. [15]. Bohlooly and Fard [16] introduced new results for buckling and postbuckling of concentrically stiffened piezo-composite plates on elastic foundations. Chi and Chung [17] found an analytical solution based on the classical plate theory for the mechanical behavior of fully simply supported FGM plates subjected to transverse loading. Nam et al. [18] presented finite modeling for the mechanical buckling of cracked stiffened FGM plates based on the first-order shear theory of Mindlin. Trabelsi et al. [19] explored the thermal postbuckling behavior of functionally graded plates and cylindrical shells using four-node element based on a modified first-order shear deformation theory. These authors continued using this approach to investigate the thermal buckling of functionally graded plates and cylindrical shells [20]. Chen et al. [21] studied the buckling and bending behavior of a functionally graded porous plate, with the formulations based on the first-order shear theory and the Chebyshev–Ritz method.
There are many plate theories (from classical to higher-order shear deformation theories) that we can apply to analyze the mechanical behavior of structures made of anisotropic and isotropic materials. The Shi shear deformation theory [22] is a higher-order shear deformation theory that has many advantages, as discussed in detail in the literature [22,23], and it gives an accurate solution for the analysis of shear flexible plates. This theory can be developed to solve the nonlinear constitutive behavior of materials [24] and problems of materials with misfitting inclusions [25].
In all of the above published works, many results and conclusions were achieved on the buckling behavior of unstiffened and stiffened plates. However, a detailed study on the percentage weight loss of a stiffened FGM plate compared with a unstiffened one (these two plates having the same buckling strength) has not been done despite this research being very important in structural design and manufacturing. This paper presents a finite element formulation for mechanical buckling responses of stiffened FGM plates based on the G. Shi shear deformation theory. New numerical results were computed to examine the effect of different parameters on the buckling problem of stiffened FGM plates. This study demonstrates clearly the decrease in the weight of a stiffened FGM plate compared with a unstiffened one. Furthermore, this work studied the effects of the distance between two stiffeners on the buckling loads of a stiffened plate in order to determine if the distance results in a buckling strength higher or lower than that of a plate with one central stiffener.
This paper is structured as follows. Section 2 shows the finite element method using the Shi shear deformation theory for the buckling problem of stiffened FGM plates. The numerical results of the buckling loads for stiffened FGM plates are computed and discussed in Section 3. Section 4 gives the main conclusions of this study.

2. Finite Element Formulation for Mechanical Buckling of Stiffened FGM Plates

We considered a functionally graded material stiffened plate composed of ceramic and metal phases. The material on the top surface of this plate was full of ceramic and was graded to metal at the bottom surface of the plate by the power law distribution. The stiffener was placed at the bottom surface, and it was full of metal. This meant that the material of the stiffener and the bottom surface was the same.
The thickness, length, and width of the plate have been noted as h, a, and b, respectively, while the depth and width of the stiffener are hs and bs, respectively, as sketched in Figure 1. The material was graded by the power law distribution, and it was used for describing the volume fraction of the ceramic (Vc) and the metal (Vm) as follows [23,26]:
V m + V c = 1
and
V c ( z ) = ( 1 2 + z h p ) n
where, h is the thickness of the plate; n is the gradient index (n ≥ 0); z is the thickness coordinate variable (−h/2 ≤ z ≤ h/2); and subscripts c and m represent the ceramic and metal constituents, respectively.
In this study, the Young’s modulus E and the Poisson’s ratio ν change through the z-direction as follows [23,26]:
E ( z ) = E m + ( E c E m ) V c ;   ν ( z ) = ν m + ( ν c ν m ) V c
Using the Shi shear deformation theory [22,23], the FGM plate model has the following displacement field (u, v, w):
u ( x , y , z ) = u 0 ( x , y ) + 5 4 ( z 4 3 h 2 z 3 ) ϕ x ( x , y ) + ( 1 4 z 5 3 h 2 z 3 ) w 0 , x
v ( x , y , z ) = v 0 ( x , y ) + 5 4 ( z 4 3 h 2 z 3 ) ϕ y ( x , y ) + ( 1 4 z 5 3 h 2 z 3 ) w 0 , y
w ( x , y , z ) = w 0 ( x , y )
where u0, v0, and w0 represent the displacements at z = 0 (the mid-plane of a plate); ϕ x and ϕ y are the transverse normal rotations of the y and x axes; the comma denotes the differentiation with respect to x and y coordinates.
Four nodes per element, seven degrees of freedom per node were used for this problem. The displacement vector of node i for plate element is as follows:
q e i = { u 0 i   v 0 i   w i   ϕ x i   ϕ y i   w 0 i x   w 0 i y } T ,   i = 1 ÷ 4
Because of the degree of freedom, w had the additional first derivative components w 0 i x ,   w 0 i y . Therefore, in order to guarantee the continuous condition of displacement w and its first derivative components at each node, we had to approximate the displacement w by Hermite interpolation functions. The other four degrees of freedom were approximated by Lagrange interpolation functions.
The displacements of the plate in this approach may be expressed as follows:
{ u 0 ,   v 0 ,   ϕ x ,   ϕ y } = { i = 1 4 N i u 0 i ,   i = 1 4 N i v 0 i ,   i = 1 4 N i ϕ x i ,   i = 1 4 N i ϕ y i }
w = H 1 w 01 + H 2 w 01 x + H 3 w 01 y + + H 10 w 04 + H 11 w 04 x + H 12 w 04 y
w x = x ( H 1 w 01 + H 2 w 01 x + H 3 w 01 y + + H 10 w 04 + H 11 w 04 x + H 12 w 04 y )
w y = y ( H 1 w 01 + H 2 w 01 x + H 3 w 01 y + + H 10 w 04 + H 11 w 04 x + H 12 w 04 y )
where Ni are Lagrange interpolating functions and Hi are Hermite interpolating functions.
The displacement vector is interpolated through the element’s nodal displacement vector as follows:
u 0 = B H · q e
where B H is the interpolation function matrix; u 0 and q e are expressed as follows:
u 0 = { u 0 , v 0 , w 0 , ϕ x , ϕ y , w 0 , x , w 0 , y } T
q e = { q e 1   q e 2   q e 3   q e 4 } T
The total strain of this plate in the case of the plate subjected to in-plane prebuckling stresses can be written as follows:
{ ε γ } = { ε ( 0 ) γ ( 0 ) } + z { ε ( 1 ) 0 } + z 2 { 0 γ ( 2 ) } + z 3 { ε ( 3 ) 0 } + { ε P G }
with
ε ( 0 ) = { u 0 x v 0 y u 0 y + v 0 x } , ε ( 1 ) = 1 4 { 5 ( ϕ x x + 2 w x 2 ) 5 ( ϕ y y + 2 w y 2 ) ( ϕ x y + 2 2 w x y + ϕ y x ) } , ε ( 3 ) = 5 3 h 2 { ϕ x x + 2 w x 2 ϕ y y + 2 w y 2 ϕ x y + 2 2 w x y + ϕ y x }
γ ( 0 ) = 1 4 { 5 ϕ x + w x 5 ϕ y + w y } ; γ ( 2 ) = 5 h 2 { ϕ x + w x ϕ y + w y } ; ε P G = { 1 2 ( w x ) 2 + z 2 2 ( ϕ x x ) 2 + z 2 2 ( ϕ y x ) 2 1 2 ( w y ) 2 + z 2 2 ( ϕ x y ) 2 + z 2 2 ( ϕ y y ) 2 w x w y + z 2 ϕ x x ϕ x y + z 2 ϕ y x ϕ y y 0 0 }
By substituting Equation (4) into Equations (13) and (14), the strain field can be obtained as follows:
ε = ( B 1 + B 2 + B 3 ) q e ; γ = ( B 4 + B 5 ) q e
with
B 1 = i 4 [ N i , x 0 0 0 0 0 N i , y 0 0 0 N i , y N i , x 0 0 0 ]
B 2 = 1 4 i = 1 4 [ 0 0 ( H ( 3 i 2 ) , x ) , x 5 N i , x 0 ( H ( 3 i 1 ) , x ) , x ( H ( 3 i ) , x ) , x 0 0 ( H ( 3 i 2 ) , y ) , y 0 5 N i , y ( H ( 3 i 1 ) , y ) , y ( H ( 3 i ) , y ) , y 0 0 2 ( H ( 3 i 2 ) , x ) , y 5 N i , y 5 N i , x 2 ( H ( 3 i 1 ) , x ) , y 2 ( H ( 3 i ) , x ) , y ]
B 4 = 5 4 i = 1 4 [ 0 0 H ( 3 i 2 ) , x N i 0 H ( 3 i 1 ) , x H 3 i , x 0 0 H ( 3 i 2 ) , y 0 N i H ( 3 i 1 ) , y H 3 i , y ]
B 5 = - 5 h 2 i = 1 4 [ 0 0 H ( 3 i 2 ) , x N i 0 H ( 3 i 1 ) , x H 3 i , x 0 0 H ( 3 i 2 ) , y 0 N i H ( 3 i 1 ) , y H 3 i , y ]
The constitutive relations are derived from Hooke’s law by the following equation:
σ = D m ( z ) ( ε ( 0 ) + z ε ( 1 ) + z 3 ε ( 3 ) ) ;   τ = D s ( z ) ( γ ( 0 ) + z 2 γ ( 2 ) )
with
σ = { σ x , σ y , τ x y } T ;   τ = { τ x z , τ y z } T
D m ( z ) = E ( z ) 1 v 2 ( z ) [ 1 v ( z ) 0 v ( z ) 1 0 0 0 ( 1 v ( z ) ) / 2 ] ,   D s ( z ) = E ( z ) 2 ( 1 + v ( z ) ) [ 1 0 0 1 ]
In this work, the stiffener was assumed to be parallel to the x-axis (See Figure 2). There was no delamilation phenomenon between the stiffener and the plate during the performance of the structure; the stiffener seemed to be a beam, and it just bent in the zx-plane. The displacement field of the x-stiffener can be expressed as follows:
u x s ( x , y , z , t ) = u 0 x s ( x , y , t ) + 5 4 ( z 4 3 h x s 2 z 3 ) ϕ x x s ( x , y , t ) + ( 1 4 z 5 3 h x s 2 z 3 ) w 0 , x x s
v x s ( x , y , z , t ) = 0
w x s ( x , y , z , t ) = w 0 x s ( x , y , t )
The displacement vector of node i for the x-stiffener element is written as follows:
q i x s = { u 0 i x s   0   w 0 i x s   ϕ x i x s   0   w 0 , x x s   0 } T ,   i = 1 2
The displacements of the stiffener can now be given in the following form:
{ u 0 x s ,   ϕ x x s } T = { i = 1 2 N i x s u 0 i x s , i = 1 2 N i x s ϕ x i x s } T
w o x s = H 1 x s w 01 x s + H 2 x s ( w x ) 1 + H 3 ( 0 ) 1 + H 4 x s w 02 x s + H 5 x s ( w x ) 2 + H 6 x s ( 0 ) 2
w 0 x s x = x ( H 1 x s w 01 x s + H 2 x s x ( w 0 x s x ) 1 + H 3 x s x ( 0 ) 1 + H 4 x s x w 02 x s + H 5 x s x ( w 0 x s x ) 2 + H 6 x s x ( 0 ) 2 )
where N i x s are Lagrange interpolating functions, and H i x s are Hermite interpolating functions. These functions can be obtained by substituting s = s0 into Ni and Hi of the plate.
The displacement vector of x-stiffener is interpolated as follows:
u 0 x s = B H x s · q e x s
with
u 0 x s = { u 0 x s , 0 , w 0 x s , ϕ x x s , 0 , w 0 , x x s , 0 } T
q e xs = { q 1 e xs   q 2 e xs } T
The strain components of the stiffener are as follows:
{ ε x s γ x s } = { ε x s ( 0 ) γ x s ( 0 ) } + z { ε x s ( 1 ) 0 } + z 2 { 0 γ x s ( 2 ) } + z 3 { ε x s ( 3 ) 0 } + { ε x s G }
with
ε x s ( 0 ) = { u 0 , x x s 0 0 } ;   ε x s ( 1 ) = 1 4 { 5 ( ϕ x , x x s + w , x x x s ) 0 0 } ;   ε x s ( 3 ) = 5 3 h x s 2 { ϕ x , x x s + w , x x x s 0 0 }
γ x s ( 0 ) = 1 4 { 5 ϕ x x s + w , x x s 0 } ; γ x s ( 2 ) = 5 h x s 2 { ϕ x x s + w , x x s 0 } ; ε x s G = { 1 2 ( w , x x s ) 2 + z 2 2 ( ϕ x , x x s ) 2 0 0 0 0 }
By substituting Equation (24) into Equations (27)–(29), the strain field is as follows:
ε x s = ( B 1 x s + B 2 x s + B 3 x s ) q e x s ; γ x s = ( B 4 x s + B 5 x s ) q e x s
with
B 1 x s = i 2 [ N i , x x s 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ]
B 2 x s = 1 4 i = 1 2 [ 0 0 H ( 3 i 2 ) , x x x s 5 N i , x x s 0 H ( 3 i 1 ) , x x x s H ( 3 i ) , x x x s 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ]
B 3 x s = 5 3 h x s 2 i = 1 2 [ 0 0 H ( 3 i 2 ) , x x x s N i , x x s 0 H ( 3 i 1 ) , x x x s H ( 3 i ) , x x x s 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ]
B 4 x s = 5 4 i = 1 2 [ 0 0 H ( 3 i 2 ) , x x s N i x s 0 H ( 3 i 1 ) , x x s H ( 3 i ) , x x s 0 0 0 0 0 0 0 ]  
B 5 x s = 5 h x s 2 [ 0 0 H ( 3 i 2 ) , x x s N i x s 0 H ( 3 i 1 ) , x x s H ( 3 i ) , x x s 0 0 0 0 0 0 0 ]
The relationship between stresses and strains obtained from Hooke’s law is as follows:
σ x s = D m x s ( ε x s ( 0 ) + z . ε x s ( 1 ) + z 3 ε x s ( 3 ) ) ; τ x s = D s x s ( γ x s ( 0 ) + z 2 γ x s ( 2 ) )
with
σ x s = [ σ x x s 0 0 ] T , τ x s = [ τ x z x s 0 ] T
D m x s = E x s 1 v x s 2 [ 1 v x s 0 v x s 1 0 0 0 ( 1 v x s ) / 2 ]
D s x s = E x s 2 ( 1 + v x s ) [ 1 0 0 1 ]
The x-stiffener is considered to place at the lower surface of the plate. The condition of displacement at the contact line is as follows:
u 0 | z = 0.5 h = u 0 | z = 0.5 h x s x s
Using Equations (3) and (19), Equation (34) becomes the following:
{ u 0 i x s = u 0 i + e 1 x s ϕ x i + e 2 x s w 0 i , x ϕ x i x s = ϕ x i ; w 0 i , x x s = w 0 i , x
where
e 1 x s = 5 12 ( h + h x s ) ; e 2 x s = 1 12 ( h + h x s )
Equation (35) can be rewritten as follows:
{ u 0 i x s 0 w 0 i x s ϕ x i x s 0 w 0 i , x x s 0 } = [ 1 0 0 e 1 x s 0 e 2 x s 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 ] { u 0 i v 0 i w 0 i ϕ x i ϕ y i w 0 i , x w 0 i , y }
or in shorter form:
u 0 i x s = T x · u 0 i
The nodal displacement vector of stiffener element is as follows:
q e x s = T x s · q e
in which
T x s = T x d i a g ( 4 , 4 )
The elastic strain energy of the stiffened plate is written as follows:
U = 1 2 N p V e ε T · σ · d V + 1 2 N x s V e x s ε x s T · σ x s · d V
or in matrix form:
U = 1 2 N p q e T K e p q e + 1 2 N x s ( q e ) T K e x s q e
where
K e p = S e [ B 1 T AB 1 + B 1 T BB 2 + B 1 T EB 3 + B 2 T BB 1 + B 2 T DB 2 + B 2 T FB 3 + B 3 T EB 1 + + B 3 T FB 2 + B 3 T HB 3 + B 4 T A B 4 + B 4 T B B 5 + B 5 T B B 4 + B 5 T D B 5 ] d S
K e x s = T x s T ( b x s l e x s [ ( B 1 x s ) T A x s B 1 x s + ( B 1 x s ) T B x s B 2 x s + ( B 1 x s ) T E x s B 3 x s + ( B 2 x s ) T B x s B 1 x s + + ( B 2 x s ) T   D ˜ x s B 2 x s + ( B 2 x s ) T F x s B 3 x s + ( B 3 x s ) T E x s B 1 x s + ( B 3 x s ) T F x s B 2 x s + + ( B 3 x s ) T H x s B 3 x s + ( B 4 x s ) T A x s B 4 x s + ( B 4 x s ) T B x s B 5 x s + + ( B 5 x s ) T B x s B 4 x s + ( B 5 x s ) T D x s B 5 x s ] d l ) T x s
in which
( A ,   B ,   D ,   E ,   F ,   H ) = h / 2 h / 2 ( 1 ,   z ,   z 2 ,   z 3 ,   z 4 ,   z 6 ) D m · d z
( A ,   B ,   D ) = h / 2 h / 2 ( 1 , z 2 ,   z 4 ) D s · d z
( A x s ,   B x s ,   D ˜ x s ,   E x s ,   F x s ,   H x s ) = h x s / 2 h x s / 2 ( 1 ,   z ,   z 2 ,   z 3 ,   z 4 ,   z 6 ) D m x s · d z
( A x s ,   B x s ,   D x s ) = h x s / 2 h x s / 2 ( 1 , z 2 ,   z 4 ) D s x s · d z
The geometric strain energy enforced by in-plane prebuckling stresses is then computed by the following:
U G = 1 2 N p V e   σ ^ 0 T · ε P G · d V + 1 2 N x s V e x s   σ ^ 0 x s T · ε x s G · d V
By substituting the geometric strain of the plate and the stiffeners into Equation (49), we get the following:
U G = 1 2 N p S e (   ε ¯ P G ) T σ 0 ·   ε ¯ P G · d S + 1 2 N x s S e x s (   ε ¯ x s G ) T σ 0 x s   ε ¯ x s G d S
where
  ε ¯ P G = { w 0 , x w 0 , y ϕ x , x ϕ x , y ϕ y , x ϕ y , y } = [ 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 x 0 0 0 0 0 0 y 0 0 0 0 0 0 0 x 0 0 0 0 0 0 y 0 0 ] { u v w ϕ x ϕ y w 0 , x w 0 , y } = L P G q e
σ 0 = d i a g ( h   σ ˜ ,   h 3 12   σ ˜ ,   h 3 12   σ ˜ ) ;   σ ˜ = [ σ x 0   τ x y 0 ; τ x y 0   σ y 0 ] T
ε ¯ x s G = { w 0 , x ϕ x , x } = [ 0 0 0 0 0 1 0 0 0 0 / x 0 0 0 ] { u x s 0 w 0 x s ϕ x x s 0 w 0 , x x s 0 } = L x s G q e x s
σ 0 x s = d i a g ( h x s σ ˜ x s ,   h x s 3 12 σ ˜ x s ,   h x s 3 12 σ ˜ x s ) ; σ ˜ x s = [ σ x 0   0 ; 0   0 ]
Equation (50) now becomes the following:
U G = 1 2 N p q e T K e G p q e + 1 2 N x s ( q e ) T K e G x s q e
with
K e G p = S e ( L p G ) T σ 0 L p G d S ;   K e G x s = T x s T ( b x s l e ( L x s G ) T σ 0 x s L x s G d l ) T x s
For the buckling problem, we get the following equation:
{ ( K p + K xs ) λ b ( K Gp + K Gxs ) } d =   0  
where K p ,   K x s and K G p ,   K G x s are the global stiffness matrix and global geometric stiffness matrix, respectively. d stands for the vector of unknowns. Equation (57) is solved to obtain the buckling load λ b and the buckling mode shape.

3. Numerical Results

3.1. Formulation Verification

First, we conducted a comparison of the critical buckling loads for a simply supported FGM plate with the analytical results of Meisam et al. [27] and Bodaghi et al. [28], as shown in Table 1. We considered a square plate with a = b = 1 m, the thickness h = a/100, the material properties Ec = 380 × 109 Pa, Em = 70 × 109 Pa, the Poisson ratio 0.3, and the plate as being under an axial compression at two opposite edges.
Next, we compared the buckling coefficient of the clamped plate with one central stiffener with a model, as shown in Figure 3. The plate had geometrical parameters a/b = 1 and thickness of h = a/200; the depth of the stiffener was hs = 10.483h, and the width was bs = hs/2.75. The material properties of the plate and the stiffener were E = 68.7 GPa, ν = 0.3. The buckling coefficient was compared with the results of Mukhopadhyay et al. [29] (semianalytical finite difference method), Peng et al. [30] (mesh-free method), and Rikards et al. [31] (finite element method). This comparison is listed in Table 2.
k b u k = λ b · a 2 / ( D · π 3 )
with D = E h 3 / 12 ( 1 ν 2 ) .
From Table 1 and Table 2, we can see clearly that the result of our work compared with other approaches has a very small error, so the calculation program used in this paper is verified.

3.2. Buckling of Rectangular FGM Plate with One Central Stiffener

3.2.1. Long FGM Plate with One Central Stiffener

In this work, two types of FGM plates (Si3N4/SUS304 plate, ZrO2/SUS304 plate) were employed. The material properties are presented in Table 3. We considered a rectangular plate with one central stiffener with b = 0.2 m, a/b = 2.5. Here, the stiffener was made of the same material as the material of the plate.
The buckling coefficient is calculated as follows:
k b u k = λ b u k · a 0 2 / ( E 0 · h 0 3 )
where E0 = 5 × 107 Pa, a0 = 10h0.
-Effect of the Depth of Stiffener (hs)
We evaluated how the depth of the stiffener affects the buckling behavior of the stiffened FGM plate. The volume fraction index n had a value in the range from 0.1 to 10. Considering three cases of plates with b/h = 100, 150, and 200, the geometrical parameters of the stiffener were bs = 2h, hs = h − 4h; the plate was fully simply supported (SSSS). We also determined the mass reduction of the stiffened structure compared with the plate without stiffeners with the same dimensions (the length and the width of this plate). For the stiffened plate with the thickness h = b/100, the buckling coefficients are represented by the horizontal dash-dot lines in the Figure 4. For each of the plates without stiffeners, we let the thickness of the plate change from b/95 to b/70. The buckling coefficients are represented by the solid lines in Figure 4 and are listed in Table 4. The thickness relationship between the Si3N4/SUS304 plate with and without a central stiffener is listed in Table 5 and Table 6. The red point, which is the intersection point between the dash-dot line and the solid line, represents the equivalent buckling coefficient of the plate without stiffeners and the stiffened plate. Following the same procedure for the stiffened plate with the thickness h = b/150 (the thickness of the unstiffened plate changed from b/140 to b/110) and h = b/200 (the thickness of the unstiffened plate changed from b/180 to h = b/155), we obtained the results given in Figure 5 and Table 7, Table 8, Table 9, Table 10, Table 11, Table 12, Table 13, Table 14, Table 15 and Table 16.
From Figure 4 and Figure 5 and Table 4, Table 5, Table 6, Table 7, Table 8, Table 9, Table 10, Table 11, Table 12, Table 13, Table 14, Table 15 and Table 16, we can come to the following conclusions:
-
The depth of stiffener (hs) has a strong influence on the buckling coefficient of the structure. For both Si3N4/SUS304 and ZrO2/SUS304 plates, when hs increased, the buckling coefficient increased, as we can see from the horizontal lines in the above figures.
-
For the Si3N4/SUS304 plate, with the three thicknesses h considered (h = b/100, h = b/150, h = b/200), when n was small, the addition of the stiffener enhanced the mass of the structure compared with the plate without stiffeners. In contrast, when n was large, the addition of stiffener significantly reduced the mass of the structure in order to obtain the same buckling coefficient between the stiffened plate (thinner plate) and the plate without a central stiffener (thicker plate). The maximum mass reduction was up to 13.17%, 18.16%, and 20% for cases h = b/200, h = b/150, h = b/100, respectively, herein hs = 4h.
-
For the ZrO2/SUS304 plate, for all values of n, the addition of stiffener decreased the mass of the structure compared with the unstiffened plate (thicker plate), but they still had the same buckling load. The larger the volume index fraction n got, the larger the mass reduction reached. The maximum mass reduction was up to 15.07%, 20.74%, and 23.26% for cases h = b/200, h = b/150, h = b/100, respectively, herein hs = 4h.
These numerical results are very interesting as it obviously demonstrates that the addition of stiffener still ensured the stability of the structure under a compressive load and that it could significantly reduce the mass of the structure due to the thickness of the plate decreasing. This point is really significant for many cases in practical engineering, where the mass of the structure must be reduced to make sure the process works under compressive loads, especially in aerospace engineering, shipbuilding and so on.
-Effect of the gradient index (n)
Next, the authors evaluated the influence of the gradient index n on the buckling response of the plate. We considered two cases of plates made of Si3N4/SUS304 or ZrO2/SUS304 with geometrical parameters hs = 2h, bs = b/100 and b/h ratio value of 10, 20, 50, and 100; the plate was fully simply supported (SSSS). By changing the volume fraction index n from 0.1 to 10, we had a variation of buckling coefficients of FGM plates with one central stiffener and without stiffeners, as shown in Figure 6.
Note that in Figure 6, the solid lines and dash-dot lines represent buckling coefficients of stiffened and unstiffened plates, respectively. The solid lines are always higher than the dash-dot lines, meaning that the buckling load of the stiffened plate was larger than that of the unstiffened plate. This proves that the stiffness of the plate will increase by the addition of the stiffener. We can see that when the thickness of the plate increased with a central stiffener, the plate became stiffer, leading to an enhancing of the buckling load of the plate. When the volume fraction index n was increased, the buckling load of the ZrO2/SUS304 plate got larger; in contrast, the buckling load of the Si3N4/SUS304 plate got smaller.
Figure 7 presents the first four buckling mode shapes of the FGM plate with and without a central stiffener. The figures show clearly that the stiffener had a strong influence on the buckling mode shape of the plate.
-Effect of the length-to-width ratio (a/b)
Next, we examined the effects of the length-to-width ratio (a/b) on the buckling load of the plate with one central stiffener. In this case, we constantly maintained the width of the plate as b = 0.2 m and the width of the stiffener as bs = b/100; the plate was under fully clamped (CCCC). By changing the length of the plate a, the a/b ratio reached the value from 2 to 4. We then obtained the diagram of the buckling load, as shown in Figure 8. The figure shows that when the length-to-width ratio (a/b) increased, the plate became “softer”, meaning the buckling load decreased.
Figure 9 presents the first four buckling mode shapes of the stiffened plate with a/b = 2 and a/b = 4. The results demonstrate that the length and width of the plate strongly affected the buckling mode shape of the structure.
-Effect of the boundary condition
Next, we investigated the influence of the boundary condition on the buckling load of the stiffened FGM plate. Figure 10 and Figure 11 present the buckling coefficient of the structure. The results in Figure 10 were calculated for plates with various b/h ratios, and the results in Figure 11 were calculated for plates with various a/b ratios. From these figures, we can see that the fully clamped supported plate had the largest buckling load, and the fully simply supported plate had the smallest buckling load.

3.2.2. Long FGM Plate with Two Stiffeners

We considered a rectangular plate with two parallel stiffeners, as shown in Figure 12, with b = 0.2 m, a/b = 2.5. The geometrical parameters of the stiffener were hs = 2h, bs = b/100. The buckling coefficient of the plate can be determined by the following formula:
k b u k = λ b u k a 0 2 / ( E 0 h 0 3 )
with E0 = 5 × 107 Pa, a0 = 10h0.
-Effect of the width of the stiffener and the distance between two stiffeners
The effects of the width of the stiffener and the distance between two stiffeners were studied. We compared the buckling loads between the plate with two stiffeners (bs = b/100, bs = b/200, Figure 13b) and the one with one central stiffener (bs = b/100, Figure 13a).
By varying the volume index fraction n and the distance between the two stiffeners d, we obtained the numerical results of the buckling coefficients of the plate, as plotted in Figure 14.
From Figure 14, we can reach the following conclusions. First, Figure 14 clearly shows that the distance between the two stiffeners strongly influences the buckling coefficient of the structure. In both cases, i.e., Si3N4/SUS304 and ZrO2/SUS304 plates, when d increased, the buckling coefficient of the plate decreased. This is explained by the fact that the strain energy focuses on the central area of the structure. Therefore, when the stiffeners are set there, the stiffness of this plate is larger than other places, thus leading to enhanced buckling coefficient of the plate.
Second, similar to the previous sections, when the volume index fraction n increased, the buckling coefficient of the Si3N4/SUS304 plate decreased, and the buckling coefficient of the ZrO2/SUS304 plate increased.
Third, we explored the different buckling coefficients between a plate with a central stiffener and a plate with two stiffeners. We considered a plate with two stiffeners and a plate with a central stiffener with the same geometrical parameters except the width of the stiffener bs. In case 1, we considered a plate with two stiffeners, which had the total bs of two stiffeners (each stiffener had bs/2) equal to the bs of the plate with a central stiffener (Figure 13a,b). For both computed results of Si3N4/SUS304 and ZrO2/SUS304 plates (Figure 14a,b), the buckling coefficient of the plate with one central stiffener was always larger than the buckling coefficient of the plate with two stiffeners. When the distance d was small, the plate with two stiffeners had the same buckling coefficient as a plate with one central stiffener. In case 2, we considered a plate with two stiffeners, with the width of each stiffener bs equal to bs of the plate with a central stiffener (Figure 13a,c). As can be seen in Figure 14b,d, for both Si3N4/SUS304 and ZrO2/SUS304 plates, when d = 0.5b, the buckling coefficient of the structure was the same as the plate with one stiffener. When the distance of stiffeners was d < 0.5b, the buckling coefficient of the plate with two stiffeners was larger than that of the plate with a central stiffener. When d > 0.5b, the buckling coefficient of the plate with two stiffeners was smaller than that of the plate with a central stiffener. These interesting results have a good significance in practical engineering. For case 1, all plates had the same volume of the stiffener, but using the plate with one central stiffener was much better due to its buckling load always being higher than that of the plate with two stiffeners. In case 2, the plate with two stiffeners was only better than that with one stiffener in the buckling problem when the distance between the two stiffeners was suitable.
-Effect of the plate thickness
In this study, the thickness of the plate (the plate model is given in Figure 12) had values from b/50 to b/5, and the numerical results are shown in Figure 15. We can again see the same phenomenon as in the case of the plate with a central stiffener. When the volume fraction index n increased, the buckling coefficient of Si3N4/SUS304 plate increased, and the buckling coefficient of ZrO2/SUS304 plate decreased. The reason for this is that the plate became softer.
Figure 16 presents the first four buckling mode shape of the plate with b/h = 5 and b/h = 50. We can see that the thickness of the structure strongly affected the buckling loads as well as buckling mode shapes.
-Effect of the length-to-width ratio (a/b)
We next explored the effects of the length-to-width ratio (a/b). The geometrical parameters were b = 0.2 m, h = b/10. We changed the length a and the volume fraction index n, and the results are plotted in Figure 17. As can be seen from the figure, the length-to-width ratio (a/b) had a robust influence on the buckling coefficient of the structure. This point can be explained by the fact that when the a/b ratio decreased, the length a correspondingly increased. This led to the plate becoming “softer”, and the buckling coefficient was therefore reduced.
Figure 18 presents the first four buckling mode shapes of the FGM plate for two cases, i.e., a/b = 2 and a/b = 5; the plate is fully simply supported (SSSS). We can see that both the boundary condition and the a/b ratio strongly affected not only the buckling loads but also the buckling mode shapes of the structure.

4. Conclusions

In this study, a linear finite element formulation based on the Shi shear deformation theory was employed to study the linear mechanical buckling of stiffened functionally graded material plates. Different parameters were examined, including the geometrical and mechanical properties of the plate. Based on the numerical results presented in the above sections, some major conclusions are listed as follows:
-
The volume index fraction n has a strong influence on the buckling load of the structure. Due to the different type of material, the buckling coefficient of the Si3N4/SUS304 plate decreased as the volume index fraction n increased. In contrast, the buckling coefficient of the ZrO2/SUS304 plate increased as the volume index fraction n increased.
-
The geometrical parameters of both the plate and the stiffener also strongly affect the buckling coefficient of the structure. For the Si3N4/SUS304 and ZrO2/SUS304 plates with one central stiffener or two stiffeners, the buckling load increased as the thickness of the stiffener increased.
-
The addition of stiffener to the plate can significantly reduce the total mass of the structure, but it still maintains the stability of the plate with the same buckling load as the plate without stiffeners. The maximum mass reduction could be reached at 20% for the Si3N4/SUS304 plate and 23.26% for the ZrO2/SUS304 plate. This point is very significant in practical engineering when we need to have lighter structures. Moreover, using a plate with two stiffeners was better than using a plate with one central stiffener when we had the right distance (d < 0.5b) between the two stiffeners.
-
Boundary conditions also deeply affect the buckling load of the structure. The plate with CCCC boundary condition had much more buckling load than the others.

Author Contributions

Investigation, H.-N.N.; Software, T.C.T., D.T.L.; Visualization, V.-D.P., P.V.M.; Writing—review & editing, supervision D.V.T.

Funding

This research was funded by Vietnam National Foundation for Science and Technology Development (NAFOSTED) grant number 107.02-2018.30.

Acknowledgments

D.V.T. gratefully acknowledges the supports of the Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 107.02-2018.30.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Broekel, J.; Prusty, B.G. Experimental and theoretical investigations on stiffened and un-stiffened composite panels under uniform transverse loading. Compos. Struct. 2004, 63, 293–304. [Google Scholar] [CrossRef]
  2. Liu, Y.; Wang, Q. Computational study of strengthening effects of stiffenerson regular and arbitrarily stiffened plates. Thin-Walled Struct. 2012, 59, 78–86. [Google Scholar] [CrossRef]
  3. Ueda, Y.; Rashed, S.M.H.; Paik, J.K. An incremental Galerkin method for plates and stiffened plates. Comput. Struct. 1987, 11, 147–156. [Google Scholar]
  4. Danielson, D.A.; Cricellx, A.S.; Frenzen, C.L.; Vasudevan, N. Buckling of stiffened plates under axial compression and lateral pressure. Int. J. Solids Struct. 1993, 30, 545–551. [Google Scholar] [CrossRef]
  5. Jiang, W.; Burt, G.; Roberts, J.C. Finite element modeling of stiffened and un-stiffened orthotropic plates. Comput. Struct. 1997, 63, 105–117. [Google Scholar] [CrossRef]
  6. Hughes, O.F.; Ghosh, B.; Chen, Y. Improved prediction of simultaneous local and overall buckling of stiffened panels. Thin-Walled Struct. 2004, 42, 827–856. [Google Scholar] [CrossRef]
  7. Pavlovcic, L.; Detzel, A.; Kuhlmann, U.; Beg, D. Shear resistance of longitudinallystiffened panels—Part 1: Tests and numerical analysis of imperfections. J. Constr. Steel Res. 2007, 63, 337–350. [Google Scholar] [CrossRef]
  8. Javaheri, R.; Eslami, M.R. Thermal Buckling of Functionally Graded Plates. AIAA J. 2002, 40, 161–169. [Google Scholar] [CrossRef]
  9. Shariat, B.A.; Eslami, M.R. Buckling of thick functionally graded plates under mechanical and thermal loads. Compos. Struct. 2007, 78, 433–439. [Google Scholar] [CrossRef]
  10. Prakash, T.; Singha, M.K.; Ganapathi, M. Thermal postbuckling analysis of FGM skew plates. Eng. Struct. 2008, 30, 22–32. [Google Scholar] [CrossRef]
  11. Moita, J.S.; Araújo, A.L.; Correia, V.F.; Mota Soares, C.M.; Herskovits, J. Buckling and nonlinear response of functionally graded plates under thermo-mechanical loading. Compos. Struct. 2018, 202, 719–730. [Google Scholar] [CrossRef]
  12. Van Dung, D.; Nam, V.H. Nonlinear dynamic analysis of eccentrically stiffened functionally graded circular cylindrical thin shells under external pressure and surrounded by an elastic medium. Eur. J. Mech. A Solids 2014, 46, 42–53. [Google Scholar] [CrossRef]
  13. Yu, T.; Bui, T.Q.; Yin, S.; Doan, D.H.; Wu, C.T.; Van Do, T.; Tanaka, S. On the thermal buckling analysis of functionally graded plates with internal defects using extended isogeometric analysis. Compos. Struct. 2016, 136, 684–695. [Google Scholar] [CrossRef]
  14. Van Dung, D.; Nga, N.T. Thermomechanical postbuckling analysis of eccentrically stiffened FGM sandwich plates with general Sigmoid and power laws based on TSDT. J. Sandw. Struct. Mater. 2018, 20, 907–945. [Google Scholar] [CrossRef]
  15. Nam, V.H.; Thi Phuong, N.; Huy Bich, D. Buckling analysis of parallel eccentrically stiffened functionally graded annular spherical segments subjected to mechanic loads. Mech. Adv. Mater. Struct. 2018. [Google Scholar] [CrossRef]
  16. Bohlooly, M.; Malekzadeh Fard, K. Buckling and postbuckling of concentrically stiffened piezo- composite plates on elastic foundations. J. Appl. Comput. Mech. 2018, 5, 128–140. [Google Scholar]
  17. Chi, S.; Chung, Y. Mechanical behavior of functionally graded material plates under transverse load—Part I: Analysis. Int. J. Solids Struct. 2006, 43, 3657–3674. [Google Scholar] [CrossRef] [Green Version]
  18. Hoai, N.V.; Doan, D.H.; Khoa, N.M.; Van Do, T.; Tran, H.T. Phase-field buckling analysis of cracked stiffened functionally graded plates. Compos. Struct. 2019, 217, 50–59. [Google Scholar] [CrossRef]
  19. Trabelsi, S.; Frikha, A.; Zghal, S.; Dammak, F. Thermal post-buckling analysis of functionally graded material structures using a modified FSDT. Int. J. Mech. Sci. 2018, 144, 74–89. [Google Scholar] [CrossRef]
  20. Trabelsi, S.; Frikha, A.; Zghal, S.; Dammak, F. A modified FSDT-based four nodes finite shell element for thermal buckling analysis of functionally graded plates and cylindrical shells. Eng. Struct. 2019, 178, 444–459. [Google Scholar] [CrossRef]
  21. Chen, D.; Yang, J.; Kitipornchai, S. Buckling and bending analyses of a novel functionally graded porous plate using Chebyshev-Ritz method. Arch. Civ. Mech. Eng. 2019, 19, 157–170. [Google Scholar] [CrossRef]
  22. Shi, G. A new simple third-order shear deformation theory of plates. Int. J. Solids Struct. 2007, 44, 4399–4417. [Google Scholar] [CrossRef] [Green Version]
  23. Bui, T.Q.; Thom, V.D.; That Ton, L.H.; Doan, D.H.; Tanaka, S.; Pham, D.T.; Nguyen, V.T.; Yu, T.; Hirose, S. On the high temperature mechanical behaviors analysis of heated functionally graded plates using FEM and a new third-order shear deformation plate theory. Compos. Part B 2016, 92, 218–241. [Google Scholar] [CrossRef]
  24. Lellep, J.; Majak, J. Nonlinear constitutive behavior of orthotropic materials. Mech. Compos. Mater. 2000, 36, 261–266. [Google Scholar] [CrossRef]
  25. Mori, T.; Tanaka, K. Average stress in matrix and average elastic energy of materials with misfitting inclusions. Acta Metall. 1973, 21, 571–574. [Google Scholar] [CrossRef]
  26. Reddy, J.N. Analysis of Functionally Graded Plates. Int. J. Numer. Method Eng. 2000, 47, 663–684. [Google Scholar] [CrossRef]
  27. Meisam, M.; Ali, R.S.; Emad, J. Levy solution for buckling analysis of functionally graded rectangular plates. Appl. Compos. Mater. 2010, 17, 81–93. [Google Scholar]
  28. Bodaghi, M.; Saidi, A.R. Levy-type solution for buckling analysis of thick functionally graded rectangular plates based on the higher-order shear deformation plate theory. Appl. Math. Model. 2010, 34, 3659–3673. [Google Scholar] [CrossRef]
  29. Mukhopadhyay, M. Vibration and stability analysis of stiffened plates by semi-analytical finite defference method-part I: Consideration of bending displacements only. J. Sound Vib. 1989, 130, 27–39. [Google Scholar] [CrossRef]
  30. Peng, L.X.; Liew, K.M.; Kitipornchai, S. Buckling and free vibration analyses of stiffened plates using mesh-free method. J. Sound Vib. 2006, 289, 421–449. [Google Scholar] [CrossRef]
  31. Rikards, R.; Chate, A.; Ozolinsh, O. Analysis for buckling and vibration of composite stiffened shells and plates. Compos. Struct. 2001, 51, 361–370. [Google Scholar] [CrossRef]
Figure 1. A functionally graded material (FGM) plate stiffened by an x-direction stiffener.
Figure 1. A functionally graded material (FGM) plate stiffened by an x-direction stiffener.
Materials 12 01262 g001
Figure 2. Plate element and stiffener element.
Figure 2. Plate element and stiffener element.
Materials 12 01262 g002
Figure 3. FGM plate with one central stiffener.
Figure 3. FGM plate with one central stiffener.
Materials 12 01262 g003
Figure 4. Buckling coefficient of FGM (Si3N4/SUS304) plate. Solid blue line: buckling coefficient of plate without stiffeners (b/h = 70–95); horizonal dash-dot line: buckling coefficient of plate with a central stiffener (b/h = 100) with various hs/h ratios. (a) n = 0.1; (b) n = 0.5; (c) n = 1; (d) n = 2; (e) n = 5; (f) n = 10.
Figure 4. Buckling coefficient of FGM (Si3N4/SUS304) plate. Solid blue line: buckling coefficient of plate without stiffeners (b/h = 70–95); horizonal dash-dot line: buckling coefficient of plate with a central stiffener (b/h = 100) with various hs/h ratios. (a) n = 0.1; (b) n = 0.5; (c) n = 1; (d) n = 2; (e) n = 5; (f) n = 10.
Materials 12 01262 g004aMaterials 12 01262 g004b
Figure 5. Buckling coefficient of FGM (ZrO2/SUS304) plate. Solid blue line: buckling coefficient of plate without stiffeners (b/h = 70–95); horizontal dash-dot line: buckling coefficient of plate that has central stiffener (b/h = 100) with various hs/h ratios. (a) n = 0.1; (b) n = 0.5; (c) n = 1; (d) n = 2; (e) n = 5; (f) n = 10.
Figure 5. Buckling coefficient of FGM (ZrO2/SUS304) plate. Solid blue line: buckling coefficient of plate without stiffeners (b/h = 70–95); horizontal dash-dot line: buckling coefficient of plate that has central stiffener (b/h = 100) with various hs/h ratios. (a) n = 0.1; (b) n = 0.5; (c) n = 1; (d) n = 2; (e) n = 5; (f) n = 10.
Materials 12 01262 g005
Figure 6. Variation of buckling coefficient of FGM plate with one central stiffener and without stiffeners (solid line: plate with one central stiffener, dash-dot line: plate without stiffeners). (a) Si3N4/SUS304; (b) ZrO2/SUS304.
Figure 6. Variation of buckling coefficient of FGM plate with one central stiffener and without stiffeners (solid line: plate with one central stiffener, dash-dot line: plate without stiffeners). (a) Si3N4/SUS304; (b) ZrO2/SUS304.
Materials 12 01262 g006
Figure 7. The first four buckling mode shapes of FGM plate with and without a central stiffener (b/h = 50, n = 0.5, Si3N4/SUS304).
Figure 7. The first four buckling mode shapes of FGM plate with and without a central stiffener (b/h = 50, n = 0.5, Si3N4/SUS304).
Materials 12 01262 g007
Figure 8. Variation of the buckling coefficient of the plate by the variation of the length-to-width ratio (a/b). (a) Si3N4/SUS304; (b) ZrO2/SUS304.
Figure 8. Variation of the buckling coefficient of the plate by the variation of the length-to-width ratio (a/b). (a) Si3N4/SUS304; (b) ZrO2/SUS304.
Materials 12 01262 g008
Figure 9. The first four buckling mode shape of stiffened FGM plate (b/h = 10, n = 0.5, Si3N4/SUS304).
Figure 9. The first four buckling mode shape of stiffened FGM plate (b/h = 10, n = 0.5, Si3N4/SUS304).
Materials 12 01262 g009
Figure 10. Variation of the buckling coefficient by the boundary condition and the thickness of the plate.
Figure 10. Variation of the buckling coefficient by the boundary condition and the thickness of the plate.
Materials 12 01262 g010
Figure 11. Variation of the buckling coefficient by the boundary condition and the length of the plate.
Figure 11. Variation of the buckling coefficient by the boundary condition and the length of the plate.
Materials 12 01262 g011
Figure 12. Model of the plate with two parallel stiffeners.
Figure 12. Model of the plate with two parallel stiffeners.
Materials 12 01262 g012
Figure 13. Plate with one stiffener and two stiffeners. (a) Plate with one stiffener; (b) plate with two stiffeners.
Figure 13. Plate with one stiffener and two stiffeners. (a) Plate with one stiffener; (b) plate with two stiffeners.
Materials 12 01262 g013
Figure 14. Buckling coefficient of the plate with one stiffener and two stiffeners. (a,c) Plate with two stiffeners bs = b/200 and plate with one central stiffener; (b,d) plate with two stiffeners bs = b/100 and plate with one central stiffener.
Figure 14. Buckling coefficient of the plate with one stiffener and two stiffeners. (a,c) Plate with two stiffeners bs = b/200 and plate with one central stiffener; (b,d) plate with two stiffeners bs = b/100 and plate with one central stiffener.
Materials 12 01262 g014aMaterials 12 01262 g014b
Figure 15. Variation of the buckling coefficient of the plate with two parallel stiffeners by the thickness of the structure. (a) Si3N4/SUS304; (b) ZrO2/SUS304.
Figure 15. Variation of the buckling coefficient of the plate with two parallel stiffeners by the thickness of the structure. (a) Si3N4/SUS304; (b) ZrO2/SUS304.
Materials 12 01262 g015
Figure 16. The first four buckling mode shapes of FGM plate with two stiffeners (n = 0.5, Si3N4/SUS304).
Figure 16. The first four buckling mode shapes of FGM plate with two stiffeners (n = 0.5, Si3N4/SUS304).
Materials 12 01262 g016
Figure 17. Variation of the buckling coefficient of the plate with two parallel stiffeners by changing of a/b ratio. (a) Si3N4/SUS304; (b) ZrO2/SUS304.
Figure 17. Variation of the buckling coefficient of the plate with two parallel stiffeners by changing of a/b ratio. (a) Si3N4/SUS304; (b) ZrO2/SUS304.
Materials 12 01262 g017
Figure 18. The first four buckling mode shapes of stiffened FGM plate (b/h = 10, n = 0.5, SSSS, Si3N4/SUS304).
Figure 18. The first four buckling mode shapes of stiffened FGM plate (b/h = 10, n = 0.5, SSSS, Si3N4/SUS304).
Materials 12 01262 g018aMaterials 12 01262 g018b
Table 1. Comparison of the critical buckling load (MN/m) for a simply supported FGM plate (a = b = 1 m, h = a/100).
Table 1. Comparison of the critical buckling load (MN/m) for a simply supported FGM plate (a = b = 1 m, h = a/100).
Buckling Loadn = 0n = 1n = 2
Meisam [27] (analytical method)1.37370.68470.5343
Bodaghi [28] (analytical method)1.37300.68440.5340
This work1.37630.68610.5353
Table 2. Buckling coefficient of the clamped plate with a single stiffener.
Table 2. Buckling coefficient of the clamped plate with a single stiffener.
Buckling CoefficientRikards [31]Rikards [31] AnsysMukhopadhyay [29]Peng [30]This Work
k b u k 24.8523.4425.4625.3326.26
Table 3. Properties of materials.
Table 3. Properties of materials.
MaterialSi3N4ZrO2SUS304
E (GPa)322.27168.06207.79
ν 0.240.2980.3218
Table 4. Buckling coefficient of FGM plate with and without a central stiffener (Si3N4/SUS304, b/h = 100).
Table 4. Buckling coefficient of FGM plate with and without a central stiffener (Si3N4/SUS304, b/h = 100).
nWith Stiffenerb/h = 100 (Si3N4/SUS304)
hs/h43.532.521.51
0.1 10.7699.00437.58576.48245.6695.11084.7654
0.59.97468.57297.17176.07855.27044.71454.3698
19.54568.33236.94365.85815.05434.50044.1564
29.22868.14516.76965.69234.89344.3423.9992
58.91377.95196.59185.52484.73214.18433.8432
108.65957.80216.45185.39124.60244.05693.717
Table 5. The thickness relationship between Si3N4/SUS304 plate with and without a central stiffener (b/h = 100).
Table 5. The thickness relationship between Si3N4/SUS304 plate with and without a central stiffener (b/h = 100).
nWith Stiffenerb/h = 100 (Si3N4/SUS304)
hs/h43.532.521.51
0.1without stiffenerb/h74.5 79.083.9 89.092.5--
0.574.1 78.083.188.092.1--
173.977.482.087.192.1--
273.6 76.981.687.091.6--
573.576.480.986.391.0--
1073.476.180.785.990.194.9 -
Table 6. The mass reduction between Si3N4/SUS304 plate with and without a central stiffener (b/h = 100).
Table 6. The mass reduction between Si3N4/SUS304 plate with and without a central stiffener (b/h = 100).
nWith Stiffenerb/h = 100 (Si3N4/SUS304)
hs/h43.532.521.51
0.1without stiffener 1.0302 1.01260.9917 0.96990.9636--
0.51.1199 1.09541.05811.02551.0029--
11.5861.13291.09521.05421.0172--
21.18921.16171.11761.06881.0325--
51.21101.18551.14011.08721.0466--
101.20001.19631.14781.09591.05951.0192-
Table 7. Buckling coefficient of FGM plate with and without a central stiffener (Si3N4/SUS304, b/h = 150).
Table 7. Buckling coefficient of FGM plate with and without a central stiffener (Si3N4/SUS304, b/h = 150).
nwith Stiffenerb/h = 150 (Si3N4/SUS304)
hs/h43.532.521.51
0.1 2.6252.25691.96231.73571.57031.45771.3884
0.52.49912.1351.84271.61751.45271.34041.2712
12.42942.0681.77751.55341.38921.27711.208
22.37582.01731.72871.50571.34221.23041.1614
52.32071.96561.67941.45791.29531.1841.1153
102.27751.92471.63991.41941.25741.14651.078
Table 8. The mass reduction between Si3N4/SUS304 plate with and without a central stiffener (b/h = 150).
Table 8. The mass reduction between Si3N4/SUS304 plate with and without a central stiffener (b/h = 150).
nwith Stiffenerb/h = 150 (Si3N4/SUS304)
hs/h43.532.521.51
0.1without stiffener 0.97820.96700.95010.9455---
0.51.06151.03691.01630.9945---
11.1077 1.06541.03901.0127---
21.13781.09761.06161.0336---
51.16531.12011.07891.0438---
101.18161.13421.08991.04881.0244--
Table 9. Buckling coefficient of FGM plate with and without a central stiffener (Si3N4/SUS304, b/h = 200).
Table 9. Buckling coefficient of FGM plate with and without a central stiffener (Si3N4/SUS304, b/h = 200).
nWith Stiffenerb/h = 200 (Si3N4/SUS304)
hs/h43.532.521.51
0.1 0.98070.86050.76490.69190.63880.60270.5807
0.50.92920.810.71490.64220.58920.55330.5312
10.9010.78240.68780.61530.56250.52660.5045
20.87960.76170.66760.59540.54270.50690.4849
50.85790.74090.64740.57560.52310.48740.4655
100.84060.72420.63110.55950.50720.47160.4497
Table 10. The mass reduction between Si3N4/SUS304 plate with and without a central stiffener (b/h = 200).
Table 10. The mass reduction between Si3N4/SUS304 plate with and without a central stiffener (b/h = 200).
nWith Stiffenerb/h = 200 (Si3N4/SUS304)
hs/h43.532.521.51
0.1without stiffener 0.93580.9347-----
0.51.02591.0044-----
11.06871.0306-----
21.09291.06171.0322----
51.11941.08161.0513----
101.13171.09601.0645----
Table 11. Buckling coefficient of FGM plate with and without a central stiffener (ZrO2/SUS304, b/h = 100).
Table 11. Buckling coefficient of FGM plate with and without a central stiffener (ZrO2/SUS304, b/h = 100).
nwith Stiffenerb/h = 100 (ZrO2/SUS304)
hs/h43.532.521.51
0.1 6.78756.29485.41744.40883.65193.12472.7944
0.57.10696.61275.59634.57793.81523.28482.9528
17.27526.78025.6914.66723.90133.36913.0362
27.41646.92085.77014.74183.97323.43963.1061
57.59187.09565.86474.83234.06163.52693.1932
107.72787.23145.93594.90134.12963.59463.2608
Table 12. The mass reduction between ZrO2/SUS304 plate with and without a central stiffener (b/h = 100).
Table 12. The mass reduction between ZrO2/SUS304 plate with and without a central stiffener (b/h = 100).
nwith Stiffenerb/h = 100 (ZrO2/SUS304)
hs/h43.532.521.51
0.1without stiffener 1.13651.14061.11751.07691.0410--
0.51.17941.17921.14011.09131.0567--
11.20111.19791.15191.10101.0645--
21.21781.21221.15911.10691.0671--
51.22961.22231.16661.11331.0647--
101.23261.22421.16501.11381.0614--
Table 13. Buckling coefficient of FGM plate with and without a central stiffener (ZrO2/SUS304, b/h = 150).
Table 13. Buckling coefficient of FGM plate with and without a central stiffener (ZrO2/SUS304, b/h = 150).
nwith Stiffenerb/h = 150 (ZrO2/SUS304)
hs/h43.532.521.51
0.1 1.85061.62011.34851.13660.97980.8720.8051
0.51.94471.67251.39851.18511.02740.91910.8520
11.99421.70011.42491.21061.05250.94390.8766
22.03581.72321.44691.23191.07350.96470.8973
52.08761.7511.47371.25821.09940.99050.9231
102.11521.7721.49411.27841.11951.01060.9432
Table 14. The mass reduction between ZrO2/SUS304 plate with and without a central stiffener (b/h = 150).
Table 14. The mass reduction between ZrO2/SUS304 plate with and without a central stiffener (b/h = 150).
nwith Stiffenerb/h = 150 (ZrO2/SUS304)
hs/h43.532.521.51
0.1without stiffener 1.11501.10231.06401.0339---
0.51.15471.12761.08711.0497---
11.17551.14471.10021.0589---
21.19231.15061.10671.0659---
51.20361.15741.11191.0695---
101.20741.15321.10721.0683---
Table 15. Buckling coefficient of FGM plate with and without a central stiffener (ZrO2/SUS304, b/h = 200).
Table 15. Buckling coefficient of FGM plate with and without a central stiffener (ZrO2/SUS304, b/h = 200).
nwith Stiffenerb/h = 200 (ZrO2/SUS304)
hs/h43.532.521.51
0.1 0.71140.60040.51090.44170.39090.35610.3347
0.50.73350.62160.53140.46180.410860.3760.3545
10.74520.63270.54230.47250.42130.38640.3649
20.7550.6420.55130.48130.43010.39510.3736
50.76670.65340.56240.49230.4410.4060.3845
100.77560.66210.5710.50080.44950.41450.3929
Table 16. The mass reduction between ZrO2/SUS304 plate with and without a central stiffener (b/h = 200).
Table 16. The mass reduction between ZrO2/SUS304 plate with and without a central stiffener (b/h = 200).
nwith Stiffenerb/h = 200 (ZrO2/SUS304)
hs/h43.532.521.51
0.1without stiffener 1.08551.05731.0300----
0.51.11481.0856-----
11.12991.0962-----
21.14271.1068-----
51.15281.1144-----
101.15071.1126-----

Share and Cite

MDPI and ACS Style

Nguyen, H.-N.; Tan, T.C.; Luat, D.T.; Phan, V.-D.; Thom, D.V.; Minh, P.V. Research on the Buckling Behavior of Functionally Graded Plates with Stiffeners Based on the Third-Order Shear Deformation Theory. Materials 2019, 12, 1262. https://doi.org/10.3390/ma12081262

AMA Style

Nguyen H-N, Tan TC, Luat DT, Phan V-D, Thom DV, Minh PV. Research on the Buckling Behavior of Functionally Graded Plates with Stiffeners Based on the Third-Order Shear Deformation Theory. Materials. 2019; 12(8):1262. https://doi.org/10.3390/ma12081262

Chicago/Turabian Style

Nguyen, Hoang-Nam, Tran Cong Tan, Doan Trac Luat, Van-Duc Phan, Do Van Thom, and Phung Van Minh. 2019. "Research on the Buckling Behavior of Functionally Graded Plates with Stiffeners Based on the Third-Order Shear Deformation Theory" Materials 12, no. 8: 1262. https://doi.org/10.3390/ma12081262

APA Style

Nguyen, H. -N., Tan, T. C., Luat, D. T., Phan, V. -D., Thom, D. V., & Minh, P. V. (2019). Research on the Buckling Behavior of Functionally Graded Plates with Stiffeners Based on the Third-Order Shear Deformation Theory. Materials, 12(8), 1262. https://doi.org/10.3390/ma12081262

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