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
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
The constitutive relation for the FGM plate subjected to thermo-mechanical loadings is
The stress resultants can be written as
Substituting Equations (1)–(3) into Equation (4) yields
where
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
where
z0 is the distance between the geometrical midplane and the neutral surface. Substituting Equation (7) into Equation (5) yields
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]
Using Equation (8) and employing the Euler equations to minimize the total potential energy function yields the equilibrium equations of an FGM plate
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
or
where
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
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:
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
which is independent of thickness
z.
The temperature distribution of the linear temperature field is
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:
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
Substituting Equation (16) into Equation (21) yields the nonlinear temperature field as
where
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:
where
w is the buckling mode function, and
where
ϑ and
η are unknown spatial eigenvalues. Substituting Equations (24) and (25) into Equation (12) yields the relation of
ϑ,
η and
as
From Equation (26), we have the following relations
Therefore, the functions
ϕ(
x) and
ψ(
y) have the following closed forms:
where the five unknown variables
α1,
α2,
β1,
β2 and
are involved; then, four equations are needed for a unique solution. From Equation (27), one can obtain three relations of these five variables:
Substituting
ϑ = ±i
α1,
η = ±i
β1 into Equation (30) yields
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 . For clarity of discussion, the solved from eigenvalue equations is denoted by .
Note that
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
. 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
where
(
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),
can be written as
then, the analytical critical temperature can be solved from Equation (31). For a uniform temperature field, it is
where Δ
Tcr =
Tcr −
T0; for a linear temperature field, it is
It can be seen that, for the temperature-independent homogeneous isotropic thin plates, we have
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,
follows from Equation (14) that
Then, the critical buckling temperature is
where
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
can be written as
where
It can be proved that Δ
Tcr is one of the roots of the following ninth-order algebraic equation
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
because it is difficult to determine the derivative d
/d
T. Let
x =
T and
Given two initial values
x−1 and
x0, the iterative scheme of the secant method is
When
xk and
xk−1 satisfy
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 (Si
3N
4). 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.