Next Article in Journal
The Effects of Thermoplastic Polyurethane on the Structure and Mechanical Properties of Modified Polypropylene Blends
Next Article in Special Issue
Effects of Order of Expansion for the Exponential Matrix and Number of Mathematical Layers in the Exact 3D Static Analysis of Functionally Graded Plates and Shells
Previous Article in Journal
Effect of Equal Channel Angular Pressing (ECAP) on Erosion-Corrosion of Pure Copper
Previous Article in Special Issue
Free Vibration Analysis of Functionally Graded Porous Doubly-Curved Shells Based on the First-Order Shear Deformation Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Closed Form Solutions for Thermal Buckling of Functionally Graded Rectangular Thin Plates

Institute of Solid Mechanics, Beihang University, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2017, 7(12), 1256; https://doi.org/10.3390/app7121256
Submission received: 31 October 2017 / Revised: 27 November 2017 / Accepted: 30 November 2017 / Published: 3 December 2017

Abstract

:
This work concerns the critical buckling temperature of functionally graded rectangular thin plates; the properties of functionally graded material vary continuously in accordance with the power law of thickness z. Closed form solutions for the critical thermal parameter are obtained for the plate with the following boundary condition combinations: simply supported, clamped and guided edges, under uniform, linear and nonlinear temperature fields via the separation-of-variable method. Furthermore, a new method is proposed to determine the critical buckling temperature from the critical thermal parameter. The present results coincide well with those in the literature, verifying the correctness of the present method. The influences of the length–thickness ratio, length–width ratio, power law index and initial temperature on critical buckling temperature are investigated.

1. Introduction

Designable Functionally Graded Materials (FGMs) are increasingly widely used in the aerospace, nuclear, and ship industries and other various industry departments, especially in thermal environments with ultra-high temperatures and a large temperature gradient. FGMs were proposed as thermal barrier materials by Japanese scientist Koizumi [1,2] to solve many problems in the new generation of aerospace airplane thermal protection systems. The volume content of each material in FGMs is continuously distributed in the spatial position on the macroscopic scale, and the physical properties change smoothly, so that the stress concentration of FGMs can be better avoided or reduced.
Since FGMs usually work in extremely high temperature environments, in order to ensure structural stability, the critical buckling temperature of FGM structures is very important in theory and application. Many researchers have studied the thermal buckling of FGM plates based on different theories in the past two decades.
Javaheri and Eslami [3] developed closed form solutions for the thermal buckling of simply supported FGM thin plates based on the Classical thin Plate Theory (CPT). Shariat and Eslami [4] obtained closed form solutions for the thermal buckling of simply supported imperfect FGM thin plates. Kiani et al. [5] investigated the critical buckling temperature of fully clamped FGM thin plates resting on elastic foundations; they used the Bubnov–Galerkin Solution (BGS), Power Series Solution (PSS) and Semi-Levy Solution (SLS) respectively to represent the lateral displacement which cannot be obtained by the Navier solution, and compared the solutions using these different approaches. Chu et al. [6] studied the critical buckling load and buckling model of FGM thin plate with in-plane material inhomogeneity by using radial basis function, but the thermal environment was not taken into account in their work.
Shear deformation plays a more prominent role for thick plates, so shear deformation theories are more accurate than CPT in this case. Wu [7] and Bouazza et al. [8] obtained a closed form solution for thermal buckling of simply supported FGM thick plates based on the first-order shear deformation theory (FSDT), and compared the results with those of CPT. Shariat and Eslami [9] developed closed form solutions for thermal buckling of imperfect simply supported FGM thick plates based on FSDT. Lee et al. [10] obtained a finite element solution for the critical buckling temperature based on neutral surface and FSDT. Park and Kim [11] presented thermal postbuckling and vibration analyses based on FSDT by using the nonlinear finite element method. Cong and Duc [12] studied the thermal stability of eccentrically stiffened sigmoid FGM thick plate with metal–ceramic–metal layers by the Galerkin method based on FSDT.
There are also some researches based on high-order shear deformation theories. Javaheri and Eslami [13] developed closed form solutions for thermal buckling of simply supported FGM thick plates based on the third-order shear deformation theory (TSDT). Najafizadeh and Heydari [14] investigated the critical buckling temperature of clamped FGM circular thick plates based on TSDT, and compared it with results based on CPT and FSDT. Matsunaga [15] investigated the critical buckling temperature of simply supported FGM thick plates based on the 2D higher-order deformation theory that he developed, and compared the results with CPT and other lower-order shear deformation theories such as FSDT and TSDT. Zenkour and Mashat [16] studied thermal buckling of simply supported FGM thick plates based on sinusoidal shear deformation plate theory (SPT) and obtained closed form solutions.
Some researchers have used three-dimensional linear elasticity theory to analyze the thermal buckling of FGM plates. Na and Kim [17] studied the critical buckling temperature of FGM plates with various boundary conditions under uniform, linear and sinusoidal temperature fields by establishing the finite element model via three-dimensional linear elasticity theory. Subsequently, Na and Kim [18] improved their previous work and studied the thermal buckling of FGM composite structures. Malekzadeh [19] investigated the thermal buckling of FGM arbitrary straight-sided quadrilateral plates via the differential quadrature method (DQM) based on three-dimensional linear elasticity theory. Asemi et al. [20] investigated buckling of functionally graded plates under biaxial compression, shear, tension-compression, and shear-compression load conditions by establishing a full compatible three-dimensional elasticity element, but they didn’t consider the thermal environment in their work. A review paper of various modeling techniques and solution methods in thermal analysis of FGM plates was also published [21].
From the above reviews, it can be seen that closed form solutions for thermal buckling of FGM plates are only available for a few boundary conditions such as fully simply supported parallel edges. Furthermore, there is currently no literature on the closed form solution for the critical buckling temperature of FGM thin plates for other boundary conditions. Furthermore, closed form solutions can play a reference role to study the convergence and precision of various numerical methods and serve the purpose of parametric designs of structures.
In the context of closed form solutions for the critical thermal parameter of simply supported FGM rectangular thin plates, guide and clamp boundary conditions are obtained via the separation-of-variable method. A novel method is presented to determine the critical buckling temperature from the critical thermal parameter. It is assumed that the FGM properties follow the power law distribution in the thickness direction. Two types of FGM thin plates and uniform, linear and nonlinear temperature fields are taken into account. The influences of the length–thickness ratio, length–width ratio, power law index and initial temperature on the critical buckling temperature are investigated in detail.

2. Governing Equation of Thermal Buckling

In this work, the FGM rectangular plate of length a, width b and thickness h is investigated as shown in Figure 1.
CPT is used with the following displacement fields
u = u 0 z w x ν = ν 0 z w y w = w 0
where u0 and v0 are the in-plane displacements on the midplane of the plate along the x and y direction, respectively; w is the midplane deflection. The strains associated with the displacement fields in Equation (1) are
[ ε x x ε y y ε x y ] = [ u 0 x z 2 w x 2 ν 0 y z 2 w y 2 u 0 y + ν 0 x 2 z 2 w x y ]
The constitutive relation for the FGM plate subjected to thermo-mechanical loadings is
[ σ x x σ y y σ x y ] = E 1 ν 2 [ 1 ν 0 ν 1 0 0 0 1 ν 2 ] ( [ ε x x ε y y ε x y ] Δ T [ α α 0 ] )
The stress resultants can be written as
( N x x , N y y , N x y ) = h / 2 h / 2 ( σ x x , σ y y , σ x y ) d z ( M x x , M y y , M x y ) = h / 2 h / 2 z ( σ x x , σ y y , σ x y ) d z
Substituting Equations (1)–(3) into Equation (4) yields
[ N x x M x x ] = [ h / 2 h / 2 { E 1 ν 2 [ u 0 x z 2 w x 2 + ν ( ν 0 y z 2 w y 2 ) ] } d z N T h / 2 h / 2 { E z 1 ν 2 [ u 0 x z 2 w x 2 + ν ( ν 0 y z 2 w y 2 ) ] } d z M T ] [ N y y M y y ] = [ h / 2 h / 2 { E 1 ν 2 [ ν 0 y z 2 w y 2 + ν ( u 0 x z 2 w x 2 ) ] } d z N T h / 2 h / 2 { E z 1 ν 2 [ ν 0 y z 2 w y 2 + ν ( u 0 x z 2 w x 2 ) ] } d z M T ] [ N x y M x y ] = [ h / 2 h / 2 [ E 2 ( 1 + ν ) ( u 0 y + ν 0 x 2 z 2 w x y ) ] d z h / 2 h / 2 [ E z 2 ( 1 + ν ) ( u 0 y + ν 0 x 2 z 2 w x y ) ] d z ]
where
N T = h / 2 h / 2 E α Δ T 1 ν d z M T = h / 2 h / 2 E α Δ T z 1 ν d z
From Equation (5), it can be seen that the in-plane loading and bending of the plate couple together because there is in-plane stretching in the geometrical midplane due to the materials’ asymmetry with respect to the midplane. However, the coupling can be eliminated by choosing the neutral surface instead of the geometrical midplane as a reference surface in which there are no in-plane displacements [22] (p. 1091).
The in-plane displacements of the midplane can be expressed as
u 0 = z 0 w x , ν 0 = z 0 w y
where z0 is the distance between the geometrical midplane and the neutral surface. Substituting Equation (7) into Equation (5) yields
[ N x x M x x ] = [ h / 2 h / 2 [ E ( z 0 z ) 1 ν 2 ( 2 w x 2 + ν 2 w y 2 ) ] d z N T h / 2 h / 2 [ E z ( z 0 z ) 1 ν 2 ( 2 w x 2 + ν 2 w y 2 ) ] d z M T ] [ N y y M y y ] = [ h / 2 h / 2 [ E ( z 0 z ) 1 ν 2 ( 2 w y 2 + ν 2 w x 2 ) ] d z N T h / 2 h / 2 [ E z ( z 0 z ) 1 ν 2 ( 2 w y 2 + ν 2 w x 2 ) ] d z M T ] [ N x y M x y ] = [ h / 2 h / 2 ( E ( z 0 z ) 1 + ν 2 w x y ) d z h / 2 h / 2 ( E z ( z 0 z ) 1 + ν 2 w x y ) d z ]
From Equation (8), it can be seen that the in-plane loading and bending of the plate have been decoupled.
The total potential energy of a plate in a thermal environment is [23,24]
U = 1 2 V { σ x x ( ε x x α Δ T ) + σ y y ( ε y y α Δ T ) + σ x y ε x y } d V
Using Equation (8) and employing the Euler equations to minimize the total potential energy function yields the equilibrium equations of an FGM plate
N x x x + N x y y = 0 N x y x + N y y y = 0 2 M x x x 2 + 2 M y y y 2 + 2 2 M x y x y + N x x 2 w x 2 + 2 N x y 2 w x y + N y y 2 w y 2 = 0
Substituting Equation (8) into Equation (10) determines the distance z0 and yields the thermal buckling that governs the equation of an FGM thin rectangular plate as
z 0 = B 11 + G 1 A 11 + G 0
( D 11 + G 2 ( B 11 + G 1 ) 2 A 11 + G 0 ) ( 4 w x 4 + 2 4 w x 2 y 2 + 4 w y 4 ) G 0 ( 2 w x 2 + 2 w y 2 ) = 0
or
( 4 w x 4 + 2 4 w x 2 y 2 + 4 w y 4 ) + G ^ 0 ( 2 w x 2 + 2 w y 2 ) = 0
where
G ^ 0 = G 0 D 11 + G 2 ( B 11 + G 1 ) 2 A 11 + G 0 G i = h 2 h 2 z i α E 1 ν ( T T 0 ) d z ( i = 0 , 1 , 2 ) ( A i j , B i j , D i j ) = h 2 h 2 Q i j ( 1 , z , z 2 ) d z ( i , j = 1 , 2 , 6 )
Q 11 = Q 22 = E 1 ν 2 Q 12 = Q 21 = ν E 1 ν 2 Q 66 = E 2 ( 1 + ν )
where the Young’s modulus E, mass density ρ, Poisson’s ratio ν, thermal conductivity k and thermal expansion coefficient α of FGM vary continuously with thickness z as
P = ( P t P b ) ( 1 2 + z h ) n + P b
where P stands for E, ν, α, k and ρ, respectively, and the subscripts t and b stand for the top and bottom surfaces of the plate, respectively. n is the power law index indicating the material variation profile through the thickness, and n = 0 represents a homogeneous plate. The material properties are also related to temperature by the following relation:
P P 0 = P 1 T ( z ) 1 + 1 + P 1 T ( z ) + P 2 T ( z ) 2 + P 3 T ( z ) 3
where P−1, P0, P1, P2 and P3 are the constants only related to the material, and T is the temperature function.
The temperature distribution of the uniform temperature field is
T ( z ) = T 0 + Δ T
which is independent of thickness z.
The temperature distribution of the linear temperature field is
T ( z ) = T 0 + Δ T ( 1 2 + z h )
The temperature changes from T0 to T0 + ΔT linearly as the thickness varies from −h/2 to h/2.
The nonlinear temperature field is from the one-dimensional thermal conduction equation:
d d z [ k ( z ) d T d z ] = 0
Note that the thermal conductivity of materials k(z) is independent of temperature T, so the analytical solution of Equation (20) can be obtained as
T ( z ) = T 0 + Δ T h 2 z 1 k ( z ) d z h 2 h 2 1 k ( z ) d z
Substituting Equation (16) into Equation (21) yields the nonlinear temperature field as
T ( z ) = T 0 + Δ T f ( z )
where
f ( z ) = j = 0 5 ( 1 ) j Δ k j ( j × n + 1 ) k b j ( z h + 1 2 ) j × n + 1 j = 0 5 ( 1 ) j Δ k j ( j × n + 1 ) k b j Δ k = k t k b
kt and kb are thermal conductivities at the top and bottom surfaces of the plate, respectively.

3. Closed Form Solutions of Thermal Buckling

After establishing the thermal buckling governing equation, the material model and temperature fields, the buckling solutions including thermal buckling mode function and critical buckling temperature are obtained below in separation-of-variable form.
Assuming that Equation (12) has the following separation-of-variable solution:
w ( x , y ) = ϕ ( x ) ψ ( y )
where w is the buckling mode function, and
ϕ ( x ) = A e ϑ x , ψ ( y ) = B e η y
where ϑ and η are unknown spatial eigenvalues. Substituting Equations (24) and (25) into Equation (12) yields the relation of ϑ, η and G ^ 0 as
( G ^ 0 + ϑ 2 + η 2 ) ( ϑ 2 + η 2 ) = 0
From Equation (26), we have the following relations
ϑ 1 , 2 = ± α 2 = ± η 2 ϑ 3 , 4 = ± i α 1 = ± i η 2 + G ^ 0 η 1 , 2 = ± β 2 = ± ϑ 2 η 3 , 4 = ± i β 1 = ± i ϑ 2 + G ^ 0
Therefore, the functions ϕ(x) and ψ(y) have the following closed forms:
ϕ ( x ) = A 2 cos α 1 x + B 2 sin α 1 x + C 2 cosh α 2 x + D 2 sinh α 2 x ψ ( y ) = A 1 cos β 1 y + B 1 sin β 1 y + C 1 cosh β 2 y + D 1 sinh β 2 y
where the five unknown variables α1, α2, β1, β2 and G ^ 0 are involved; then, four equations are needed for a unique solution. From Equation (27), one can obtain three relations of these five variables:
α 1 2 + α 2 2 = G ^ 0 β 1 2 + β 2 2 = G ^ 0
α 1 2 + β 1 2 = ϑ 2 + η 2 + 2 G ^ 0
Substituting ϑ = ±iα1, η = ±iβ1 into Equation (30) yields
α 1 2 + β 1 2 = G ^ 0
Another two transcendental eigenvalue equations can be derived by using the two pairs of boundary conditions of four edges. The boundary conditions used in this study are presented in Table 1 where n denotes the outer normal direction of the edge.
The buckling mode functions have the separation-of-variable form as shown in Equation (24); then, the two transcendental eigenvalue equations can be derived separately by considering the boundary conditions with respect to x and y, respectively. Substituting the two eigenfunctions in Equation (28) into the boundary conditions given in Table 1 yields the transcendental eigenvalue equations which are shown in Table 2 [25,26,27].
By using Equation (29), Equation (31) and two transcendental eigenvalue equations corresponding to specific boundary conditions of a rectangular plate, we can solve the five variables α1, α2, β1, β2 and G ^ 0 . For clarity of discussion, the solved G ^ 0 from eigenvalue equations is denoted by G ^ 0 cr .
Note that G ^ 0 cr is a critical thermal buckling parameter in which the critical temperature is involved, but in general it is hard or impossible to directly extract the explicit form of the critical buckling temperature Tcr from G ^ 0 cr . In the present study, a novel method of determining Tcr is proposed. To achieve this aim, an algebraic equation for the critical buckling temperature is given as
G ^ 0 ( T ) = G ^ 0 cr
where G ^ 0 (T) is given in Equation (14); it is a function of temperature, material properties and plate thickness, and Tcr can be solved numerically from Equation (31). For three simple cases, analytical solutions can be found; the details are given in the following.

3.1. Case 1 Temperature-Independent Homogeneous Thin Plates

For homogeneous isotropic plates, the material properties are independent of geometry and temperature. According to Equation (14), G ^ 0 can be written as
G ^ 0 = G 0 D 11 + G 2 ( B 11 + G 1 ) 2 A 11 + G 0 = ( 1 + ν ) α ( h / 2 h / 2 T d z T 0 h ) h 3 12 ( 1 + ν ) α ( h / 2 h / 2 T z 2 d z T 0 h 3 12 )
then, the analytical critical temperature can be solved from Equation (31). For a uniform temperature field, it is
Δ T cr , uni = G ^ 0 cr h 2 α ( 1 + ν ) ( G ^ 0 cr h 2 + 12 )
where ΔTcr = TcrT0; for a linear temperature field, it is
Δ T cr , line = 2 G ^ 0 cr h 2 α ( 1 + ν ) ( G ^ 0 cr h 2 + 12 )
It can be seen that, for the temperature-independent homogeneous isotropic thin plates, we have
Δ T cr , line = 2 Δ T cr , uni
The nonlinear temperature field does not exist for homogeneous isotropic plates since the thermal conductivities at the top and bottom surfaces are the same. Substituting kt = kb into Equation (23), we can find that the temperature distribution of the nonlinear temperature field reduces to the linear temperature field.

3.2. Case 2 Temperature-Independent FGM Thin Plates

In this case, the FGM thin plates have temperature-independent material properties which vary with the thickness following the power law in Equation (16). For a uniform temperature field, G ^ 0 follows from Equation (14) that
G ^ 0 = Δ T h / 2 h / 2 α E 1 ν d z h / 2 h / 2 E z 2 1 ν 2 d z Δ T h / 2 h / 2 α E z 2 1 ν d z ( h / 2 h / 2 E z 1 ν 2 d z Δ T h / 2 h / 2 α E z 1 ν d z ) 2 h / 2 h / 2 E 1 ν 2 d z Δ T h / 2 h / 2 α E 1 ν d z
Then, the critical buckling temperature is
Δ T cr = G ^ 0 cr J 0 J 1 + G ^ 0 cr J 2 J 5 2 G ^ 0 cr J 3 J 4 + J 0 J 5 Θ 2 ( G ^ 0 cr J 0 J 2 G ^ 0 cr J 4 2 + J 0 2 )
where
Θ = G ^ 0 cr 2 J 0 2 J 0 1 2 G ^ 0 cr 2 J 0 J 1 J 2 J 5 4 G ^ 0 cr 2 J 0 J 1 J 3 J 4 + 4 G ^ 0 cr 2 J 0 J 2 J 3 2 + 4 G ^ 0 cr 2 J 1 J 4 2 J 5 + G ^ 0 cr 2 J 2 2 J 5 2 4 G ^ 0 cr 2 J 2 J 3 J 4 J 5 2 G ^ 0 cr J 0 2 J 1 J 5 + 4 G ^ 0 cr J 0 2 J 2 3 + 2 G ^ 0 cr J 0 J 2 J 5 2 4 G ^ 0 cr J 0 J 3 J 4 J 5 + J 0 2 J 2 5
The critical buckling temperature for linear and nonlinear temperature fields can be obtained similarly. Table 3 presents the integrals J0 ~ J5 which can be calculated numerically, and f(z) is given in Equation (23). It is noteworthy that the critical buckling temperature is irrelevant to the initial temperature for temperature-independent materials.

3.3. Case 3 Temperature-Dependent Homogeneous Plates for a Uniform Temperature Field

Referring to Equations (14) and (17), in this case G ^ 0 can be written as
G ^ 0 = 12 α 0 ( 1 + ν 0 ϒ 1 ) ϒ 2 Δ T h 2 [ 1 α 0 ( 1 + ν 0 ϒ 1 ) ϒ 2 Δ T ]
where
ϒ 1 = ν 1 T 0 + Δ T + 1 + ν 1 ( T 0 + Δ T ) + ν 2 ( T 0 + Δ T ) 2 + ν 3 ( T 0 + Δ T ) 3 ϒ 2 = α 1 T 0 + Δ T + 1 + α 1 ( T 0 + Δ T ) + α 2 ( T 0 + Δ T ) 2 + α 3 ( T 0 + Δ T ) 3
It can be proved that ΔTcr is one of the roots of the following ninth-order algebraic equation
A 9 x 9 + A 8 x 8 + A 7 x 7 + A 6 x 6 + A 5 x 5 + A 4 x 4 + A 3 x 3 + A 2 x 2 + A 1 x + A 0 = 0
where the coefficients are given in Appendix A. Equation (42) can be solved by the numerical methods to obtain the critical buckling temperature for a uniform temperature field. In general, the critical buckling temperature should be the smallest positive solution of the nine-order algebraic equation. It can be seen that the initial temperature T0 affects the critical buckling temperature of temperature-dependent homogeneous plates under a uniform temperature field. The results provided later demonstrate the effect of T0 for temperature-dependent FGM plates under all kinds of thermal environments.

3.4. Case 4 Other Situations

Except for the above three cases, it is difficult to obtain the explicit expression of the critical buckling temperature, and Equation (32) must be solved numerically for Tcr.
In this study, the secant method is used to solve
G ^ 0 ( T ) G ^ 0 cr = 0
because it is difficult to determine the derivative d G ^ 0 /dT. Let x = T and
g ( x ) = G ^ 0 ( T ) G ^ 0 cr
Given two initial values x−1 and x0, the iterative scheme of the secant method is
x k + 1 = x k g ( x k ) ( x k x k 1 ) g ( x k ) g ( x k 1 ) , ( k = 0 , 1 , )
When xk and xk−1 satisfy
| x k x k 1 | | x k | η
the iteration stops, where the relative tolerance η > 0.

4. Numerical Results and Discussion

Two types of FGM thin plates are considered in this section. Type-I FGM plate is temperature-dependent; it is made from stainless steel (SUS304) and Silicon nitride (Si3N4). The material properties are shown in Table 4. Type-II FGM plate is temperature-independent; it is made from aluminum as its bottom surface and aluminum oxide as its top surface. Et = 70 GPa, αt = 23 × 10−6/K, Eb = 380 GPa, αb = 7.4 × 10−6/K, ν = 0.3. Plates with different boundary conditions investigated below are shown in Figure 2.
Firstly, the results of the present method are compared with those in the literature to validate the correctness of this method. Table 5 and Table 6 compare the critical buckling temperature ΔTcr of simply supported type-II FGM square plate based on different theories. Table 7 compares ΔTcr of fully clamped type-II FGM rectangular plate. It can be seen that the present results coincide well with available results.
Table 8 and Table 9 show ΔTcr of type-I FGM thin plate under a uniform temperature field with different boundary conditions. The results indicate that fully clamped (CCCC) plates have the largest ΔTcr, and ΔTcr decreases as the power law index n increases. Table 10 and Table 11 give ΔTcr of Type-I FGM thin plate under different temperature fields. Due to the small difference between the temperature distribution of the linear temperature field and nonlinear temperature field, ΔTcr under these two temperature fields is close. Furthermore, ΔTcr under a uniform temperature field is obviously lower than that under the other two temperature fields, so a uniform temperature field affects ΔTcr more significantly than the other two. Table 12 presents ΔTcr of Type-I FGM thin plate with different initial temperatures T0.
In addition, the critical temperatures are also plotted in Figure 3, Figure 4, Figure 5 and Figure 6 for different thickness–length ratios h/a, power law indexes n, length–width ratios a/b and initial temperatures T0 for type-I FGM plate. Figure 3 shows the relation between ΔTcr and h/a for type-I FGM plate. The critical buckling temperature ΔTcr increases as h/a increases since the thicker plates have larger bending rigidities than the thinner ones. One can also see that the critical buckling temperature of rectangular CCCC plate is the highest, whereas that of the CSCS plate is the lowest; this is because the CCCC plate has the highest bending rigidity compared with the other two cases. Furthermore, the uniform temperature field influences ΔTcr more remarkably than the other two temperature fields.
Figure 4 demonstrates that ΔTcr decreases as n decreases. It can be concluded from Equation (16) that the volume fraction of metal increases as n increases, and metallic materials have a lower critical buckling temperature than ceramic materials; this is because the thermal expansion coefficient α of ceramic materials is less than that of metallic materials.
Figure 5 shows the relation between ΔTcr and a/b for different n. It is observed that, with a/b increasing, ΔTcr increases gradually whatever n is; this indicates that the longer plates have a higher critical buckling temperature for different n. Homogeneous isotropic ceramic plate has the highest ΔTcr and metallic plate has the lowest, and the volume fraction of the metal component increases as n increases.
Figure 6 indicates that ΔTcr of type-I FGM plate decreases almost linearly as T0 increases. For temperature-independent material, the ΔTcr stays at a constant value, thus Tcr has the same increment as that of T0. However, for type-I FGM plate, which is made from temperature-dependent materials, ΔTcr decreases almost linearly as T0 increases, and Tcr increases almost linearly, but the increment of Tcr is less than that of T0.

5. Conclusions

In this study, the closed form analytical solutions were obtained via the separation-of-variable method to determine the thermal buckling of FGM thin plates with simply supported, clamped and guided edges. A novel method to find the critical buckling temperature was given; the explicit expressions of the critical buckling temperature were presented for three special cases.
Two types of FGM thin plates and uniform, linear and nonlinear temperature fields were taken into account in this work. The present results coincide well with those in the literature, which illustrates the correctness of this present method. It follows that the critical buckling temperature of FGM thin plates decreases as the length–thickness ratio and power law index increase, and increases as the length–width ratio increases, and the initial temperature affects the critical buckling temperature of temperature-dependent materials. The uniform temperature field has the greatest influence on the critical buckling temperature, and the critical buckling temperatures become very close for both linear and nonlinear temperature fields.

Acknowledgments

The work is supported by the National Natural Science Foundation of China (11372021, 11672019).

Author Contributions

Y.X. conceived the research ideas; Z.W. derived equations, calculated and analyzed numerical results; both authors participated in the writing of this article.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Coefficients of Equation (41).
Table A1. Coefficients of Equation (41).
A9 ν 0 α 3 ν 3
A8 ν 0 ( 8 T 0 γ 1 + γ 2 )
A7 ν 0 ( 28 T 0 2 γ 1 + 7 T 0 γ 2 + γ 3 )
A6 ν 0 ( 56 T 0 3 γ 1 + 21 T 0 2 γ 2 + 6 T 0 γ 3 + γ 4 ) + α 3
A5 ν 0 ( 70 T 0 4 γ 1 + 35 T 0 3 γ 2 + 15 T 0 2 γ 3 + 5 T 0 γ 4 + γ 5 ) + ( 5 T 0 α 3 + α 2 )
A4 ν 0 ( 56 T 0 5 γ 1 + 35 T 0 4 γ 2 + 20 T 0 3 γ 3 + 10 T 0 2 2 γ 4 + 4 T 0 γ 5 + γ 6 ) + ( 10 T 0 2 α 3 + 4 T 0 α 2 + α 1 )
A3 ν 0 ( 28 T 0 6 γ 1 + 21 T 0 5 γ 2 + 15 T 0 4 γ 3 + 10 T 0 3 γ 4 + 6 T 0 2 γ 5 + 3 T 0 γ 6 + γ 7 ) + ( 10 T 0 3 α 3 + 6 T 0 2 α 2 + 3 T 0 α 1 + 1 )
A2 ν 0 ( 8 T 0 7 γ 1 + 7 T 0 6 γ 2 + 6 T 0 5 γ 3 + 5 T 0 4 γ 4 + 4 T 0 3 γ 5 + 3 T 0 2 γ 6 + 2 T 0 γ 7 + γ 8 ) + ( 5 T 0 4 α 3 + 4 T 0 3 α 2 + 3 T 0 2 α 1 + 2 T 0 + α 1 ) G ^ 0 c r h 2 α 0 + ( G ^ 0 cr h 2 + 12 )
A1 ν 0 ( T 0 8 γ 1 + T 0 7 γ 2 + T 0 6 γ 3 + T 0 5 γ 4 + T 0 4 γ 5 + T 0 3 γ 6 + T 0 2 γ 7 + T 0 γ 8 + ν 1 α 1 ) + ( T 0 5 α 3 + T 0 4 α 2 + T 0 3 α 1 + T 0 2 + T 0 α 1 ) 2 G ^ 0 cr h 2 T 0 α 0 ( G ^ 0 cr h 2 + 12 )
A0 G ^ 0 c r h 2 T 0 α 0 ( G ^ 0 c r h 2 + 12 )
where
γ 1 = α 3 ν 3 γ 2 = α 2 ν 3 + α 3 ν 2 γ 3 = α 1 ν 3 + α 2 ν 2 + α 3 ν 1 γ 4 = α 1 ν 2 + α 2 ν 1 + α 3 + ν 3 γ 5 = α 1 ν 1 + ν 3 α 1 + ν 1 α 3 + α 2 + ν 2 γ 6 = ν 2 α 1 + ν 1 α 2 + ν 1 + α 1 γ 7 = ν 1 α 1 + ν 1 α 1 + 1 γ 8 = α 1 + ν 1

References

  1. Koizumi, M. The concept of FGM. Funct. Gradient Mater. 1993, 34, 3–10. [Google Scholar] [CrossRef]
  2. Koizumi, M. FGM activities in Japan. Compos. Part B Eng. 1997, 28, 1–4. [Google Scholar] [CrossRef]
  3. Javaheri, R.; Eslami, M.R. Thermal buckling of functionally graded plates. AIAA J. 2002, 40, 162–184. [Google Scholar] [CrossRef]
  4. Shariat, B.S.; Eslami, M.R. Thermal buckling of imperfect functionally graded plates. Int. J. Solids Struct. 2006, 43, 4082–4096. [Google Scholar] [CrossRef]
  5. Kiani, Y.; Bagherizadeh, E.; Eslami, M.R. Thermal buckling of clamped thin rectangular FGM plates resting on Pasternak elastic foundation (three approximate analytical solutions). J. Appl. Math. Mech. 2011, 91, 581–593. [Google Scholar] [CrossRef]
  6. Chu, F.; He, J.; Wang, L.; Zhong, Z. Buckling analysis of functionally graded thin plate with in-plane material inhomogeneity. Eng. Anal. Bound. Elem. 2016, 65, 112–125. [Google Scholar] [CrossRef]
  7. Wu, L. Thermal buckling of a simply supported moderately thick rectangular FGM plate. Compos. Struct. 2004, 64, 211–218. [Google Scholar] [CrossRef]
  8. Bouazza, M.; Tounsi, A.; Adda-Bedia, E.A.; Megueni, A. Thermoelastic stability analysis of functionally graded plates: An analytical approach. Comput. Mater. Sci. 2010, 49, 865–870. [Google Scholar] [CrossRef]
  9. Shariat, B.S.; Eslami, M.R. Effect of initial imperfections on thermal buckling of functionally graded plates. J. Therm. Stress. 2005, 28, 1183–1198. [Google Scholar] [CrossRef]
  10. Lee, Y.H.; Bae, S.I.; Kim, J.H. Thermal buckling behavior of functionally graded plates based on neutral surface. Compos. Struct. 2016, 137, 208–214. [Google Scholar] [CrossRef]
  11. Park, J.S.; Kim, J.H. Thermal postbuckling and vibration analyses of functionally graded plates. J. Sound Vib. 2006, 289, 77–93. [Google Scholar] [CrossRef]
  12. Cong, P.H.; Duc, N.D. Thermal stability analysis of eccentrically stiffened sigmoid-FGM plate with metal–ceramic–metal layers based on FSDT. Cogent Eng. 2016, 3. [Google Scholar] [CrossRef]
  13. Javaheri, R.; Eslami, M.R. Thermal buckling of functionally graded plates based on higher order theory. J. Therm. Stress. 2002, 25, 603–625. [Google Scholar] [CrossRef]
  14. Najafizadeh, M.M.; Heydari, H.R. Thermal buckling of functionally graded circular plates based on higher order shear deformation plate theory. Eur. J. Mech. A/Solids 2004, 23, 1085–1100. [Google Scholar] [CrossRef]
  15. Matsunaga, H. Thermal buckling of functionally graded plates according to a 2D higher-order deformation theory. Compos. Struct. 2009, 90, 76–86. [Google Scholar] [CrossRef]
  16. Zenkour, A.M.; Mashat, D.S. Thermal buckling analysis of ceramic-metal functionally graded plates. Nat. Sci. 2010, 2, 968–978. [Google Scholar] [CrossRef]
  17. Na, K.S.; Kim, J.H. Three-dimensional thermal buckling analysis of functionally graded materials. Compos. Part B Eng. 2004, 35, 429–437. [Google Scholar] [CrossRef]
  18. Na, K.S.; Kim, J.H. Three-dimensional thermomechanical buckling analysis for functionally graded composite plates. Compos. Struct. 2006, 73, 413–422. [Google Scholar] [CrossRef]
  19. Malekzadeh, P. Three-dimensional thermal buckling analysis of functionally graded arbitrary straight-sided quadrilateral plates using differential quadrature method. Compos. Struct. 2011, 93, 1246–1254. [Google Scholar] [CrossRef]
  20. Asemi, K.; Shariyat, M.; Salehi, M.; Ashrafi, H. A full compatible three-dimensional elasticity element for buckling analysis of FGM rectangular plates subjected to various combinations of biaxial normal and shear loads. Finite Elem. Anal. Des. 2013, 74, 9–21. [Google Scholar] [CrossRef]
  21. Swaminathan, K.; Sangeetha, D.M. Thermal analysis of FGM plates—A critical review of various modeling techniques and solution methods. Compos. Struct. 2017, 160, 43–60. [Google Scholar] [CrossRef]
  22. Xu, T.F.; Xing, Y.F. Closed-form solutions for free vibration of rectangular FGM thin plates resting on elastic foundation. Acta Mech. Sin. 2016, 32, 1088–1103. [Google Scholar] [CrossRef]
  23. Bouhadra, A.; Benyoucef, S.; Tounsi, A.; Bernard, F.; Bouiadjra, R.B.; Sid Ahmed Houari, M. Thermal buckling response of functionally graded plates with clamped boundary conditions. J. Therm. Stress. 2015, 38, 630–650. [Google Scholar] [CrossRef]
  24. Thornton, E.A. Thermal buckling of plates and shells. Appl. Mech. Rev. 1993, 46, 485–506. [Google Scholar] [CrossRef]
  25. Xing, Y.F.; Liu, B. New exact solutions for free vibrations of rectangular thin plates by symplectic dual method. Acta Mech. Sin. 2009, 25, 265–270. [Google Scholar] [CrossRef]
  26. Xing, Y.F.; Liu, B. New exact solutions for free vibrations of thin orthotropic rectangular plates. Compos. Struct. 2009, 104, 187–195. [Google Scholar] [CrossRef]
  27. Xing, Y.F.; Xu, T.F. Solution methods of exact solutions for free vibration of rectangular orthotropic thin plates with classical boundary conditions. Compos. Struct. 2013, 104, 187–195. [Google Scholar] [CrossRef]
Figure 1. Functionally Graded Material (FGM) thin plate and Cartesian coordinates.
Figure 1. Functionally Graded Material (FGM) thin plate and Cartesian coordinates.
Applsci 07 01256 g001
Figure 2. Plates with different boundary conditions.
Figure 2. Plates with different boundary conditions.
Applsci 07 01256 g002
Figure 3. The relation between the critical buckling temperature ΔTcr and the length–thickness ratio h/a of type-I FGM rectangular thin plate (a/b = 2, n = 2, T0 = 300 K).
Figure 3. The relation between the critical buckling temperature ΔTcr and the length–thickness ratio h/a of type-I FGM rectangular thin plate (a/b = 2, n = 2, T0 = 300 K).
Applsci 07 01256 g003
Figure 4. The relation between the critical buckling temperature ΔTcr and power law index n of type-I FGM rectangular thin plate (a/h = 100, a/b = 2, T0 = 300 K).
Figure 4. The relation between the critical buckling temperature ΔTcr and power law index n of type-I FGM rectangular thin plate (a/h = 100, a/b = 2, T0 = 300 K).
Applsci 07 01256 g004
Figure 5. The relation between the critical buckling temperature ΔTcr and the length–width ratio a/b of type-I FGM, ceramic and metal rectangular thin plate (a/h = 100, n = 2, T0 = 300 K).
Figure 5. The relation between the critical buckling temperature ΔTcr and the length–width ratio a/b of type-I FGM, ceramic and metal rectangular thin plate (a/h = 100, n = 2, T0 = 300 K).
Applsci 07 01256 g005
Figure 6. (a) The relation between ΔTcr and the initial temperature T0 of type-I FGM rectangular thin plate (a/h = 100, a/b = 2, n = 2); (b) The relation between Tcr and the initial temperature T0 of type-I FGM rectangular thin plate (a/h = 100, a/b = 2, n = 2).
Figure 6. (a) The relation between ΔTcr and the initial temperature T0 of type-I FGM rectangular thin plate (a/h = 100, a/b = 2, n = 2); (b) The relation between Tcr and the initial temperature T0 of type-I FGM rectangular thin plate (a/h = 100, a/b = 2, n = 2).
Applsci 07 01256 g006
Table 1. Boundary conditions of a rectangular plate.
Table 1. Boundary conditions of a rectangular plate.
Simple Support (S) w = 0 M n = 0
Guide (G) w n = 0 Q n = 0
Clamp (C) w = 0 w n = 0
Table 2. Eigenvalue equations for various boundary conditions.
Table 2. Eigenvalue equations for various boundary conditions.
B.C.Eigenvalue Equations for y = 0 and bEigenvalue Equations for x = 0 and a
S-S sin β 1 b = 0 sin α 1 a = 0
G-G sin β 1 b = 0 sin α 1 a = 0
C-C 1 cos β 1 b cosh β 2 b sin β 1 b sinh β 2 b = β 1 2 β 2 2 2 β 1 β 2 1 cos α 1 a cosh α 2 a sin α 1 a sinh α 2 a = α 1 2 α 2 2 2 α 1 α 2
S-G cos β 1 b = 0 cos α 1 a = 0
S-C β 2 tan β 1 b = β 1 tanh β 2 b α 2 tan α 1 a = α 1 tanh α 2 a
G-C β 1 tan β 1 b = β 2 tanh β 2 b α 1 tan α 1 a = α 2 tanh α 2 a
Table 3. The integrals in Equation (38).
Table 3. The integrals in Equation (38).
IntegralUniform
Temperature Field
Linear
Temperature Field
Nonlinear
Temperature Field
J0 h / 2 h / 2 α E 1 ν d z h / 2 h / 2 α E 1 ν ( 1 2 + z h ) d z h / 2 h / 2 α E 1 ν f ( z ) d z
J1 h / 2 h / 2 E z 2 1 ν 2 d z h / 2 h / 2 E z 2 1 ν 2 ( 1 2 + z h ) d z h / 2 h / 2 E z 2 1 ν 2 f ( z ) d z
J2 h / 2 h / 2 α E z 2 1 ν d z h / 2 h / 2 α E z 2 1 ν ( 1 2 + z h ) d z h / 2 h / 2 α E z 2 1 ν f ( z ) d z
J3 h / 2 h / 2 E z 1 ν 2 d z h / 2 h / 2 E z 1 ν 2 ( 1 2 + z h ) d z h / 2 h / 2 E z 1 ν 2 f ( z ) d z
J4 h / 2 h / 2 α E z 1 ν d z h / 2 h / 2 α E z 1 ν ( 1 2 + z h ) d z h / 2 h / 2 α E z 1 ν f ( z ) d z
J5 h / 2 h / 2 E 1 ν 2 d z h / 2 h / 2 E 1 ν 2 ( 1 2 + z h ) d z h / 2 h / 2 E 1 ν 2 f ( z ) d z
Table 4. The material properties of type-I FGM plate.
Table 4. The material properties of type-I FGM plate.
PropertyMaterial TypeP−1P0P1P2P3
ESi3N40348.43 × 109−3.070 × 10-42.160 × 10−7−8.946 × 10−11
SUS3040201.04 × 1093.079 × 10-4−6.534 × 10−70
νSi3N400.2400000
SUS30400.3262−2.002 × 10-43.797 × 10−70
αSi3N405.8723 × 10-69.095 × 10-400
SUS304012.330 × 10-68.086 × 10-400
ρSi3N402370000
SUS30408166000
kSi3N409.19000
SUS304012.04000
Table 5. Critical buckling temperature ΔTcr of SSSS type-II FGM square plate with different a/h under a uniform temperature field (a/b = 1).
Table 5. Critical buckling temperature ΔTcr of SSSS type-II FGM square plate with different a/h under a uniform temperature field (a/b = 1).
nMethoda/h = 10a/h = 20a/h = 40a/h = 60a/h = 80a/h = 100
1Present782.016197.81349.60022.05612.4097.943
CPT [15]794.377198.59449.64822.06612.4127.943
TSDT [15]757.891196.25749.50022.03712.4027.939
HDT5 [15]749.260195.62349.40422.02912.4007.939
5Present713.894180.84145.36120.17311.3507.265
CPT [15]726.571181.64245.41020.18211.3527.265
TSDT [15]678.926178.52845.21320.14411.3407.260
HDT5 [15]669.402177.80545.16620.13411.3377.259
10Present733.217185.86346.62920.73711.6677.468
CPT [15]746.927186.73146.68220.74711.6707.469
TSDT [15]692.519183.14146.45520.70311.6577.462
HDT5 [15]683.211182.43246.40920.69411.6547.462
Table 6. Critical buckling temperature ΔTcr of SSSS type-II FGM plate with different a/b under a uniform temperature field (a/h = 100).
Table 6. Critical buckling temperature ΔTcr of SSSS type-II FGM plate with different a/b under a uniform temperature field (a/h = 100).
nMethoda/b = 1a/b = 2a/b = 3a/b = 4a/b = 5
1Present7.942619.851639.687667.4319103.0568
CPT [16]7.943719.859439.718867.5220103.2690
FPT [16]7.940019.835839.624867.2506102.6355
TSDT [16]7.940019.835839.624867.2506102.6356
SSDT [16]7.940019.835939.624967.2510102.6365
2Present7.041517.599535.185159.781691.3644
CPT [16]7.042617.606535.213059.862191.5538
FPT [16]7.039217.585335.128559.618491.9850
TSDT [16]7.039017.584035.123459.603790.9508
SSDT [16]7.039017.584035.123359.603490.9501
5Present7.264518.156336.296461.665994.2363
CPT [16]7.265718.164236.328561.758594.4542
FPT [16]7.261518.138036.223661.455993.7481
TSDT [16]7.260618.132736.202561.395193.6069
SSDT [16]7.260618.132436.201461.392193.5999
10Present7.467918.664537.311563.388496.8646
CPT [16]7.469218.673137.346363.488897.1005
FPT [16]7.464418.642737.224663.137896.2820
TSDT [16]7.463418.636637.200663.068796.1213
SSDT [16]7.463418.636537.200163.067396.1183
Table 7. Critical buckling temperature ΔTcr of CCCC type-II FGM rectangular plate under a uniform temperature field (a/h = 100, a/b = 2).
Table 7. Critical buckling temperature ΔTcr of CCCC type-II FGM rectangular plate under a uniform temperature field (a/h = 100, a/b = 2).
Method n = 0n = 0.5n = 1n = 2n = 5n = 10
Present32.35518.33215.03213.32613.74814.133
Kiani [5]33.62819.05315.62313.85014.28914.690
Table 8. Critical buckling temperature ΔTcr of type-I FGM rectangular thin plate with different a/h and various B.C. under a uniform temperature field (a/b = 2, T0 = 300 K).
Table 8. Critical buckling temperature ΔTcr of type-I FGM rectangular thin plate with different a/h and various B.C. under a uniform temperature field (a/b = 2, T0 = 300 K).
B.C.na/h = 20a/h = 30a/h = 40a/h = 50a/h = 60a/h = 70a/h = 80
CSCS1583.7802306.6174186.0370123.948488.133965.727850.8301
2547.3748283.8704171.2913113.782680.759360.157646.4854
5517.0574264.5071158.8189105.229074.576555.499342.8584
CSCC1783.2272431.1364267.9431181.1536130.034197.595975.8125
2742.2713401.3746247.6355166.7565119.405189.472969.4245
5713.3132376.0498230.3741154.5835110.456782.655964.0761
CCCC11057.5407640.8390413.1423285.7974208.3541158.0883123.7707
21013.9656602.6029384.3119264.3487192.0395145.3645113.6185
51030.1122571.4283359.7632246.1036178.2185134.6260105.0764
Note: ‘B.C.’ stands for Boundary Condition.
Table 9. Critical buckling temperature ΔTcr of type-I FGM rectangular thin plate with different a/b and various B.C. under a uniform temperature field (a/h = 100, T0 = 300 K).
Table 9. Critical buckling temperature ΔTcr of type-I FGM rectangular thin plate with different a/b and various B.C. under a uniform temperature field (a/h = 100, T0 = 300 K).
B.C.na/b = 2a/b = 3a/b = 4a/b = 5a/b = 6a/b = 7
CSCS132.9286 57.753792.8929136.4766187.0319243.2765
230.0846 52.836985.1409125.3608172.2154224.5874
527.7157 48.728878.6384115.9968159.6822208.7213
CSCC149.3892 100.7232 168.6840 248.9947 338.3012 433.8633
245.1642 92.3546 155.1861 229.9260 313.6459 403.9628
541.6377 85.3292143.7844213.7327292.6429378.5233
CCCC181.4595 176.7382295.2208428.0121568.5417711.2512
274.6176 162.6580273.1802398.4099532.6899671.3248
568.8856 150.7570254.4250373.2176502.7141640.2733
Table 10. Critical buckling temperature ΔTcr of type-I FGM rectangular thin plate with different a/h and various B.C. (a/b = 2, n = 2, T0 = 300 K).
Table 10. Critical buckling temperature ΔTcr of type-I FGM rectangular thin plate with different a/h and various B.C. (a/b = 2, n = 2, T0 = 300 K).
B.C.T.F.a/h = 20a/h = 30a/h = 40a/h = 50a/h = 60a/h = 70
CSCSuniform547.3748283.8704171.2913113.782680.759360.1576
linear1133.6559581.2390352.4874235.3624167.7132125.2884
nonlinear1182.5836607.9852369.1144246.6213175.8071131.3707
CSCCuniform742.2713401.3746247.6355166.7565119.405189.4729
linear1594.4982822.6864507.5211343.2718246.8452185.6015
nonlinear1655.7434859.5887531.0572359.4812258.6364194.5373
CCCCuniform1013.9656602.6029384.3119264.3487192.0395145.3645
linear2230.30481256.9658787.2933541.4994394.6294299.7612
nonlinear2275.46431310.1371822.7439566.5211413.1558313.9882
Note: ‘T.F.’ represents Temperature Field.
Table 11. Critical buckling temperature ΔTcr of type-I FGM rectangular thin plate with different a/b and various B.C. (a/h = 100, n = 2, T0 = 300 K).
Table 11. Critical buckling temperature ΔTcr of type-I FGM rectangular thin plate with different a/b and various B.C. (a/h = 100, n = 2, T0 = 300 K).
B.C.T.F.a/b = 2a/b = 3a/b = 4a/b = 5a/b = 6a/b = 7
CSCSuniform30.084652.836985.1409125.3608172.2154224.5874
linear62.9589110.1631176.7119258.9990354.3651460.7078
nonlinear66.0438115.5227185.2298271.3522371.0771482.1803
CSCCuniform45.164292.3546155.1861229.9260313.6459403.9628
linear94.2791191.5108319.7465471.5480642.0001828.0680
nonlinear98.8767200.7238334.8863493.5001671.3553865.1898
CCCCuniform74.6176162.6580273.1802398.4099532.6899671.3248
linear155.0859334.9406559.4693816.52631101.51781417.3660
nonlinear162.5833350.7718585.2726853.17691149.27821475.2611
Table 12. Critical buckling temperature ΔTcr of type-I FGM rectangular thin plate with different initial temperature T0 (a/h = 100, a/b = 2, n = 2).
Table 12. Critical buckling temperature ΔTcr of type-I FGM rectangular thin plate with different initial temperature T0 (a/h = 100, a/b = 2, n = 2).
B.C.T.F.T0 = 300T0 = 350T0 = 400T0 = 450T0 = 500T0 = 550
CSCSuniform30.084629.129328.227127.373626.564625.7966
linear62.958961.037259.190757.413155.698654.0411
nonlinear66.043864.030762.095060.230458.430656.6896
CSCCuniform45.164243.754542.421041.157539.958438.8184
linear94.279191.476988.777986.174183.657481.2201
nonlinear98.876795.942793.114890.384687.744085.1846
CCCCuniform74.617672.363170.224468.192566.25964.4162
linear155.0859150.7014146.4609142.3542138.3707134.4997
nonlinear162.5833157.9968153.5573149.2545145.0775141.0154

Share and Cite

MDPI and ACS Style

Xing, Y.; Wang, Z. Closed Form Solutions for Thermal Buckling of Functionally Graded Rectangular Thin Plates. Appl. Sci. 2017, 7, 1256. https://doi.org/10.3390/app7121256

AMA Style

Xing Y, Wang Z. Closed Form Solutions for Thermal Buckling of Functionally Graded Rectangular Thin Plates. Applied Sciences. 2017; 7(12):1256. https://doi.org/10.3390/app7121256

Chicago/Turabian Style

Xing, Yufeng, and Zekun Wang. 2017. "Closed Form Solutions for Thermal Buckling of Functionally Graded Rectangular Thin Plates" Applied Sciences 7, no. 12: 1256. https://doi.org/10.3390/app7121256

APA Style

Xing, Y., & Wang, Z. (2017). Closed Form Solutions for Thermal Buckling of Functionally Graded Rectangular Thin Plates. Applied Sciences, 7(12), 1256. https://doi.org/10.3390/app7121256

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