Abstract

This work gives information about the development of refined plate theory to study the static bending behavior of functionally graded material (FGM) plates. The significant advantage of our proposed theory is that only one unknown variable exists in its displacement formula and governing equation. To illustrate the accuracy and effectiveness of this theory, an analytical approach based on the Navier solution is employed to obtain the solution for static bending of simply supported FGM plates. A good agreement for static bending of FGM plates with other literature results has been instituted. This work also investigates the deflection, in-plane normal, and shear stresses of sinusoidally loaded FGM rectangular plates with four simply supported edges. The influence of some parameters on the bending performance of FGM plates is also carefully considered.

1. Introduction

Because of many advantageous features of functionally graded material (FGM), it has had numerous applications in some fields of engineering, for example, transportation, mechanics, and other structural applications, and the use of this structure is growing very rapidly, so it has attracted a great amount of concern from many researchers. Numerous theories on the plate have been established to investigate the dynamic and static responses as well as buckling of plates made of advanced composite material.

Firstly, the classical plate theory (CPT) has been applied to investigate the thin FGM plate by many researchers. In this plate theory, the transverse shear stress is neglected, so it cannot be applied to analyze moderate and thick plates. Javaheri and Eslami [1] used CPT to research the buckling of FGM plates subjected to in-plane compressive load. Mohammadi and his coworkers [2] analyzed the buckling behavior of FGM plates using CPT and the Levy solution. Thermal buckling of FGM plates was investigated by Ghannadpour et al. [3], in which, the finite strip method based on CPT was used. Damanpack and his colleagues [4] used the boundary element method based on CPT to investigate static bending of thin FGM plates. Nevertheless, CPT theory accommodates to analyze thin plates only because it does not consider the effects of shear deformation.

Secondly, to overcome the disadvantage of CPT, the Reissner–Mindlin plate theory or the first shear deformation theory (FSDT), in which the transverse shear deformation is assumed to be constant through the thickness of the plate, was developed to be capable of analyzing moderate and thick plates. Croce and Venini [5] used the Reissner–Mindlin plate theory to develop a hierarchic family of finite elements to analyze the static bending of FGM plates in a thermal environment subjected to mechanical loadings. Nguyen and his coworkers [6] applied the FSDT and the Galerkin solution to investigate the postbuckling behavior of FGM plates under a combination of mechanical loadings and thermal loadings, in which, the material properties of the plates depend on the temperature. Hosseini-Hashemi and his coworkers employed FSDT [7] and the Reissner–Mindlin plate theory [8] to investigate FGM plates in case of free vibration. Nguyen and his colleagues [9] applied the FSDT to develop a plate model to study the FGM plates. Shimpi and his collaborators [10] utilized a new theory based on FSDT, and Thai and his coworkers [11] established a simple theory based on FSDT which can be used to analyze the free vibration and static bending of shear deformation FGM plates. In [12], an application of a simple theory based on FSDT for the investigation of FG sandwich plates was carried out. Nguyen et al. [13] developed a refined simple FSDT to analyze the FGM plate. In this work, Nguyen presented a new shape function to describe the distribution of shear strain as well as shear stress of the plate through its thickness. A combination of the FSDT and isogeometric analysis was developed by Yu and Yin et al. [14, 15] to analyze the FGM plates. Tan-Van et al. [16] employed a modern meshless method based on FSDT to research the FGM plates. Because of its constant shear stress, FSDT needs a correction coefficient. However, this coefficient depends on the material and boundary condition.

To deal with the shortcoming of the FSDT, many scientists focus on developing higher-order shear deformation plate theories (HSDT). The first good example is the HSDT which was developed by Reddy [17] which was applied to study the plates made of FGM. The success of this theory is that it needs no correction coefficient and also the shear stresses are equal to zeros on the upper surface and the lower surface. Some other noteworthy HSDT has been carried out by many researchers such as Javaheri and Eslami [18], Bodaghi and Saidi [19], Ferreira et al. [20], and Talha and Singh [21] which have been employed to analyze FGM plates. Besides, HSDT has been combined with isogeometric for analysis of FGM plates by Tran et al. [22], composite sandwich plates by Nguyen-Xuan [23], and FGM plates by Thai et al. [24]. Zenkour [25, 26] developed generalized shear deformation theory and a HSDT with trigonometric function as well as 3D elasticity solutions to investigate thick FGM plates. Bui and his associates [27] applied HSDT and the finite element method (FEM) to investigate mechanical behaviors of the plate made of FGM in the high-temperature environment. Do et al. [28] used HSDT to analyze the bidirection FGM plate by FEM. Mantari and Guedes Soares [2932] presented some other plate theories such as generalized HSDT and generalized hybrid HSDT to exponentially study FGM plates and FGM shells. Arefi and his coworkers [33] applied a two-variable sinusoidal shear deformation theory and a physical neutral surface to investigate FGM plates with the piezoelectric layer resting on the Winkler–Pasternak foundation. Tornabene and his colleagues [34] studied double-curved, singly curved, revolution shells and plates using a generalized HSDT in combination with the local generalized differential quadrature method. Mechab and his coworkers [35] applied a new theory to inspect the static bending and dynamic response of FGM plates, in which, the number of variable unknowns of the theory is four. Benachour et al. [36] investigated the free vibration of plates made of FGM by applying a new plate theory which contains four variable unknowns. Using a novel method where the transverse displacement was separated into two parts, Shimpi [37] developed a new two-variable and single-variable refined plate theory (RPT) to study isotropic plates. His idea was developed and applied to investigate orthotropic plates by Shimpi and Patel [38], Thai and Kim [39], and Mechab et al. [40] with two unknown variables. Another extended work by Shimpi et al. was presented in [41], in which a single-variable refined theory with an inertia associated term in its displacement field is applied to study free vibrations of isotropic plates. Filippi et al. [42] developed a number of plate elements based on an improved theory based on a zig-zag power function to analyze metallic and composite layered structures with viscoelastic layers. In another extended work by Carrera and his colleagues [43], a comprehensive of classical, higher-order, zig-zag, and variable kinematic shell elements were established for analysis of composite multilayered structures. Alaimo et al. [44] presented a development of an advanced analytical formulation for damped free vibration and frequency response analysis of composite plate structures embedding viscoelastic layers model. In this work, the governing equations were derived by using the principle of virtual displacement and layer-wise models which associated to linear up to fourth-order variations of the unknown variables in the thickness direction.

Because the RPT is very simple in its formulas, it has had a large number of appliances in a lot of works to investigate many kinds of the plate including isotropic and orthotropic plates. However, according to the experiences of the authors, there are no works using this theory to investigate the plate with a continuous varying of properties such as FGM. In the current work, a development is operated to modify RPT, so this theory consists of only one unknown variable, and then it is employed to investigate the static bending problem of FGM plates. The Navier solution is occupied to solve the governing equation of fully simple supported FGM plates. The proposed theory is verified by validity studies. Besides, the investigation about the effects of some parameters on the static bending behavior of the plate is also considered and debated.

2. Material Properties of FGM Plates

The FGM plate was made by mixing two or more different ingredients with a continuous variation in the thickness of the plate. In this study, FGM plates with the power-law function (P-FGM) and exponential function (E-FGM) as shown in Figure 1 were considered.

For a P-FGM plate, the volume of ceramic is obtained using the following formula:in which is the power-law index and is the thickness of the plate. The material properties of a P-FGM can be determined aswhere are, respectively, Young’s modulus or density of the ceramic and metal.

The material properties of E-FGM can be determined aswhere is Young’s modulus or density of the material at the lower surface of the FGM plate and is the material parameter.

3. New Single-Variable Refined Plate Theory

3.1. Assumption of New Single-Variable Refined Plate Theory

The assumption of new single-variable refined plate theory based on the RPT of Shimpi [37] is given as follows:(1)The displacement is separated into two parts, the first part is the bending part and the second part is the shear part :(2)The normal stress is very small in contrast with other normal stresses. Thus, the normal stress is negligible. Therefore, by applying Hooke’s law for a linearly elastic material, the relation between normal stresses and strains is written as(3)The displacements in the -direction and in the -direction are divided into two parts, which are the bending part and the shear part:

The first parts as well as are similar to the displacements given by the CPT, which are

The second parts and give rise to shear strain and therefore to shear stresses which have distribution through the depth and are equal to zero at the top and bottom surfaces of the plate. Consequently, the expression for and can be obtained aswhere is the shear distributed profile function, and this function describes the spreading of the shear stresses throughout the thickness. In this study, the shear distributed shape function of Shimpi [37] as given in the following formula is used:

It is noticed that do not provide to moment , and .

3.2. Expressions for Displacement of Proposed Theory

By using the assumptions which are discussed above, for the case of the plate, the expressions of displacement are

The formulas for strain fields arewhere .

The formulas for normal stresses and are obtained by using equation (11) and equation (5). The formulas for , and are calculated by using equation (11) and the following constitutive equations:

Subsequently, the expressions for normal stresses and shear stresses are

Equation (13) and equation (14) can be rewritten as

The moments and shear forces are obtained as

Based on the assumption of the present theory, the moments and shear forces are

After integrating equation (19) and equation (20) over the thickness , the moments and shear forces are expressed as the following formulas:

The coefficients are calculated as

We can see that if does not depend on -direction, then . The expressions for moments do not include , and the expressions for shear forces do not include .

The equilibrium equations are

Multiply the first two equations with respect to and then integrating with respect to through the thickness, we get

Substituting equation (21) and equation (22) into equation (27), one getsand the first two equations of equation (28) becomeorin which and is the Laplace operator, .

It is clear that equation (30) reveals the relation between the bending part and the shear part of transverse displacement, and it is similar to Shimpi’s theory. However, in this theory, the coefficient not only depends on the thickness but also depends on the varying material properties, and it is an implicit integral expression. This is a significant particular different point of this theory in comparison with Shimpi’s theory. Consequently, this theory is compatible to investigate FGM plates. If Young’s modulus and Poisson ratio of the material are constant, the implicit integral expression becomes an explicit expression of Young’s modulus, Poisson ratio, and the thickness so it is completely identical with the plate theory of Shimpi.

By means of introducing equation (30) into equation (10), the new displacement formulas of the proposed theory are

It can see that the formulas of displacement of the proposed theory consist of only one unknown variable, the bending component .

3.3. Expressions for Strains, Stresses, Moments, and Shear Forces

By the way of introducing equation (31) into equation (11), new expressions for strains of proposed theory are obtained as

By introducing equation (31) into equation (15) and equation (16), the expressions for stresses of the proposed theory can be obtained as

By introducing equation (31) into equation (21) and equation (22), the expressions for moments as well as shear forces are taken as

3.4. Governing Equation

Substituting equation (35) and equation (36) into equation (27), it becomes

This is a simple governing differential equation of the static bending problem of a plate. It is clear that the governing equation of the plate consists of only one unknown variable which is the bending component . The coefficient is an implicit integral expression of Young’s modulus, Poisson ratio, and the thickness . In case of the homogeneous plate, coefficient becomes an explicit expression of Young’s modulus, Poisson ratio, and the thickness , and the governing differential equation of this theory is completely identical to that of Shimpi.

4. Analytical Solutions

In this work, the Navier procedure is applied to analyze a rectangular FGM plate, and the boundary condition of the plate is simply supported at all edges. The solution is implicit as the following formula:where , and are the quantities which must be determined. The distributed load can expand in the form of double trigonometric series asin which

Substituting equation (38) and equation (39) into equation (37), the analytical solution for the bending component of transverse displacement isand with sinusoidally distributed load, the coefficient isand with uniformly distributed load, the coefficient is

5. Illustrative Examples and Discussion

Continuously, some examples will be considered to demonstrate the accurateness and effectiveness of the proposed theory which will be applied to investigate some problems of FGM plates. The following nondimensional formulas are used in the rest of this work:

5.1. Validity Study

In this first example, a P-FGM plate is investigated, and the plate is made of two components: the first component is aluminum (Al) and the second component is alumina (Al2O3). The plate has the length of , the width of , and length-to-thickness ratio of The load which acts upon the plate are uniformly distributed load or sinusoidally load. The material properties of aluminum are  = 70 GPa and  = 0.3 and those of alumina are  = 380 GPa and  = 0.3. The nondimensional values of transverse displacement and nondimensional values as well as distribution of stresses are compared with those of the study by Zenkour [26]. The results of Zenkour are calculated using the sinusoidal shear deformation theory (SSDT). Table 1 illustrates the comparison of the results of the plate under uniform load, and Table 2 indicates the comparison of the results in the case of the sinusoidally loaded plate. According to Tables 1 and 3, it can see that the results of the current work are very similar to the solutions of Zenkour using SSDT.

Incessantly, an E-FGM plate which has a length-to-thickness ratio of and is considered. Young’s modulus of it is obtained by a function of the exponent, while Poisson’s ratios are assumed constant and equal to 0.3. The plate is subjected to bisinusoidal load. Table 2 indicates the numerical results which are calculated by the proposed theory and those of Zenkour [25] using a solution of 3D elasticity theory as well as Mantari and Guedes Soares [29] using HSDT. According to this illustration, the numerical results of the plates using the proposed theory are similar with published data for both moderate and thick plates.

According to above two examples, it can be concluded that a new theory which is developed by authors is accurate and effective, and also it can be applied to investigate FGM plates in both cases of thin and thick plates.

5.2. Static Bending Behavior of FGM Plates

sthis section, a P-FGM plate which is made of aluminum (Al) as the metal and zirconia (ZrO2) as the ceramic subjected to sinusoidally distributed load is investigated. The material properties of constituents areAl:  = 70 GPa,  = 0.3,ZrO2:  = 200 GPa,  = 0.3.

The numerical results which are computed by using the new theory are publicized in Table 4. In this example, both thick and thin plates are considered for several values of other parameters.

The results in Table 4 show that the deflection of the P-FGM plate rises as there is an increase the power-law index, especially when it varies in a range of 0 to 2. This also shows that the risk ceramic plate is harder than the risk metal plate. It is true that the highest value of the deflection is achieved in the case of the metallic plate, but the smallest value of the deflection is achieved in the case of the ceramic plate. The reason is that the elastic modulus of the metal is less than that of ceramic. In addition, when the aspect ratio increases, the deflection decreases. Moreover, the deflection of the plate increases when a/h is increased. Figures 2 and 3 show more clearly the effects of a/b and a/h on the deflection of the P-FGM plate. Figure 2 shows the nondimensional results of deflection of the P-FGM plate gathering to the ratio of a/h. According to Figure 2, the side-to-thickness ratio has principal responsibility to the bending behavior of the plate. As a/h ratio rises in a range of 2 to 4, the deflection decreases rapidly. When it is greater than 4, the deflection approximate unchanged with its increase. The exhibition in Figure 3 shows that the deflection decreases when the aspect ratio increases especially when this ratio varies from 0 to 1.

Figures 47 performed the distributions of the stresses in -direction of the plate. As exhibited in Figures 4 and 5, the in-plane normal stresses and the longitudinal tangential stress are nonlinearly distributed along with the thickness of the P-FGM plate; the neutral surface is not placed at the midsurface of the P-FGM plate, and it differs from a homogeneous plate. On the upper surface, the normal stress is tensile, while it is compressive stress on the lower surface. In addition, the maximum value of the tensile longitudinal tangential stress occurs at a point on the lower surface, while the maximum compressive longitudinal tangential stress occurs at a point on the upper surface. Figures 6 and 7 represent the dissemination of the shear stresses throughout the depth of the plate. It can easily see that the plate is free of shear strains and stresses at both upper and lower surfaces. Nevertheless, the shear stresses of an FGM plate are not parabolic, and they are not highest on the midplane as in the case of the homogeneous plate. When the aspect ratio increases, the shear stress increases as given in Figure 7.

The nondimensional distribution of normal stress in-plane longitudinal tangential stress shear stress , and shear stress in the plane of the plate are shown in Figure 8. It shows that the normal stress and shear stresses are symmetric, but the longitudinal tangential stress is unsymmetrical. The in-plane normal stress is maximum at the central point and is minimum at four corners of the plate, and this also appears with shear stress . It is clear that is compressive stress over the upper surface of the plate. In opposite sides, the shear stress is maximum at four corners of the plate and is minimum at the central point. Not to mention, the maximum longitudinal tangential stress appears at two opposite corners of the plate while the minimum value of longitudinal tangential stress appears at two other corners.

6. Conclusions

In this study, a new single-variable refined plate theory which consists of only one unknown variable in its displacement formula and its governing equation was developed. The proposed theory was successfully verified against the available literature in many cases of P-FGM and E-FGM. The proposed theory was applied to analyze P-FGM plates subjected to mechanical load in several cases of thick and thin plates. Furthermore, a large parametric investigation was aimed at checking the sensitivity of the static bending of P-FGM plates to different mechanical and geometrical properties. This comprehensive parametric investigation was presented as benchmark results for future works.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

Hoang Nam Nguyen gratefully acknowledges the support of the Ton Duc Thang University Foundation for Science and Technology Development, under grant number 2019/TDT-PLHDLV.