Research Article  Open Access
ThreeDimensional Finite Element Modeling of Thermomechanical Problems in Functionally Graded Hydroxyapatite/Titanium Plate
Abstract
The composition of hydroxyapatite (HA) as the ceramic phase and titanium (Ti) as the metallic phase in HA/Ti functionally graded materials (FGMs) shows an excellent combination of high biocompatibility and high mechanical properties in a structure. Because the gradation of these properties is one of the factors that affects the response of the functionally graded (FG) plates, this paper is presented to show the domination of the grading parameter on the displacement and stress distribution of the plates. A threedimensional (3D) thermomechanical model of a 20node brick quadratic element is used in the simulation of the thermoelastic behaviors of HA/Ti FG plates subjected to constant and functional thermal, mechanical, and thermomechanical loadings. The convergence properties of the present results are examined thoroughly in order to assess the accuracy of the theory applied and to compare them with the established research results. Instead of the grading parameter, this study reveals that the loading field distribution can be another factor that reflects the thermoelastic properties of the HA/Ti FG plates. The FG structure is found to be able to withstand the thermal stresses while preserving the high toughness properties and thus shows its ability to operate at high temperature.
1. Introduction
The term functionally graded (FG) is used to refer to inhomogeneous advanced composite materials in which the composition percentage of at least two materials is arranged to vary from one surface to another surface in a structure [1]. The properties’ gradients produced from the volume fraction variation may be tailored to meet specific needs while providing the best utilization of the composite components. The notion of functionally graded materials (FGMs) was introduced in the early 1990s in Japan as a means of overcoming problems concerning heatresistant coatings in aerospace applications [2]. The development of the relatively advanced materials was then further extended to other fields of applications, including thermoelectric and piezoelectric devices [3], biomedical engineering [4], and nuclear fast breeder and fusion reactors [5], because of their phenomenological properties. This technology has somewhat mitigated performance issues arising from abrupt transitions between dissimilar materials, namely, sharp interfaces due to adhesion that induces substantial mechanical and thermal stresses [6]. Therefore, the optimal processing technique in which several methods are combined based on the properties of the constituent materials should be implemented [7]. The best combination will preserve the high performance of the fabricated FGMs and, at the same time, reduce production costs.
A key philosophy in selecting the composition for an FGM is that the combination should be composed of the attractive properties of each component. Hydroxyapatite (HA)/titanium (Ti) FGM is a good example of this, where the excellent biocompatibility property of HA is combined with the good mechanical properties of Ti to eliminate the weakness of both materials. The combination somehow provides a structure superior to either component alone. In this study, the thermoelastic properties of the HA/Ti FG plate are explored by looking at the atomic or microscopic scale of the structure. The response of the plate due to thermal and mechanical loadings is analyzed within the consideration of various grading parameter values. In a microscopic scale, the architecture of the FGMs is contrasted against the surrounding matrix of pure materials. This factor leads to various responses in displacement and stress fields when the materials are imposed on different loading conditions.
The effect of radial and normal loads under thermal excitation on the postbuckling response of a shear deformable FG cylindrical shell has been analyzed using a singular perturbation technique in which a higher shear deformation shell theory was applied. The study concluded that the properties’ variation along the thickness direction of the FG shell and the temperature field distribution are the two parameters that affect the postbuckling characteristics and the imperfection sensitivity of the shell [8]. Matsunaga [9] developed a twodimensional (2D) model using a higherorder deformation theory to analyze the response of FG plates subjected to thermal and mechanical loads. In the computation, equations that consider the effect of normal and transverse shear stress were derived using the power series expansion of the continuous displacement components method. The study revealed that the effect of the properties’ variation along the thickness of the FG plate is larger on the response of FG plates subjected to thermal loadings. The effects of thermal and mechanical loads on the response of simply supported rectangular FG plates have been analyzed using a threedimensional (3D) thermoelastic model based on the Fourier series and statespace methods [10]. The material properties of the FG plates were assumed to vary exponentially along the thickness direction except for the Poisson ratio, which was taken as constant. It is mentioned that the location of the neutral surface is different from that of the middle surface but reflected in the Young modulus variation. The research results in Alibeigloo [10] show that the contribution of the thermal loads on the thermoelastic response of the FG plates is bigger than that of mechanical loads. The effects of thermomechanical loading on the buckling response of FG beams were analyzed by Kiani and Eslami [11] using a derivation of the Timoshenko beam theory. In the analysis, various types of boundary conditions including clamps, simply supported, and roller combination were taken into account to observe the response field distribution. This study revealed that the critical buckling temperature field is greatly affected by the constituent temperature dependency. A local meshless method based on the PetrovGalerkin approach with moving Kriging interpolation has been used to solve large deformation problems in the geometrically nonlinear thermomechanical analysis of moderately thick FG plates. The study revealed that the thermal conductivity on temperature does affect the stress field analysis on the FG plates [12]. The same local Kriging meshless method has been used for mechanical and thermal buckling analyses of FG plates in an extended numerical work. The analyses show the dependency of the FG plates’ buckling performances on the volume fraction exponent under thermal and mechanical loadings with different combined essential conditions and loading conditions [13]. The elementfree kpRitz method combined with the firstorder shear deformation shell theory and von Kármán strains is another approach used for a postbuckling analysis of CNTRFG (carbon nanotube reinforced functionally graded) cylindrical panels. The influence of volume fraction, lengthtothickness ratio, radius, boundary conditions, and distribution types of CNTs on the postbuckling behavior of the CNTRFG panels are shown in this study [14]. An extended work on the static and dynamic analyses on the CNTRFG cylindrical panels revealed that the influencing parameters studied in the earlier paper have significant influences on the flexural strength and free vibration responses of the panels [15].
In this paper, the authors will take a closer look at HA/Ti FGMs and their behavior under different loading conditions. The mechanical, thermal, and coupled thermomechanical behavior of HA/Ti FGMs will be investigated using the finite element method (FEM). The thermomechanical behavior of HA/Ti FGM plates is examined for different loading conditions and different volume fraction indices. The present accomplished results are verified and compared with published works available in the literature. The obtained results are presented in the form of contours, tables, and figures to show the parametric effects such as grading parameter and load response, which can be treated as a benchmark for future and advanced research work. These results are significant from the point of view of developing and designing thermal barrier materials, which are capable of withstanding extreme surface temperature.
2. Theoretical Formulations
The basic concept of FGMs, which produces the compositional changes in the constituents, is illustrated in Figure 1 [16]. FGMs with high mechanical properties on one side of the structures and with perfect thermal resistance properties on the other side are advantageous for various thermalresistant applications. Several advantages involving the structural strength, energy conversion, and others have been identified using different FGM combinations [17]. The main advantage of FGMs in terms of thermomechanical problems is that the alternately changed properties within their structures act as a suspension of plastic yield point and failure due to thermomechanical loads [18, 19]. The gradation process is the key element to identify the utility of FGMs. The gradual change of properties in FGMs is caused by a positiondependent chemical composition, the geometric configuration, the pore size distribution, or the physical state arrangement.
A simply supported rectangular FG plate of length , width , and thickness with a Cartesian coordinate system () having the origin considered in this study is displayed in Figure 2. The material properties of the structure are assumed to vary in the thickness direction, which is represented by the coordinate while homogeneous in the other two remaining and directions. The ceramicrich and metalrich phases in the FGM are represented by the top location () and the bottom location (), respectively.
The typical material properties () at an arbitrary location within the structural domain are calculated based on the rule of mixture [20–23], which is expressed in the following equation:
In the equation, and are the top and the bottom surface material properties, respectively, while is the volume fraction of the top surface material of the FGM plate. The Poisson ratio of each layer is assumed to be constant [24] because the effect of the properties’ variation is negligible for a unit of structure. It is noticed that most FGMs are constructed with the ceramic phase at the top surface because metal with its high toughness could be more suitable to support the structure from the bottom side. The location of each composed material should be precisely identified in order to get accurate material properties according to the rule of mixture. The following equation is used for the calculation of the volume fraction of any phase of the composed materials: where is the coordinate at the direction of the plate and is the grading parameter that dictates the properties’ variation. The FGMs will have a homogeneous ceramic ingredient when and a homogeneous metal ingredient when . The properties’ variation along the thickness direction of the plate is only available when the grading parameter is between these two values. The volume fraction of the composed materials strongly depends on the grading parameter, , as summarized in Figure 3. This parameter somehow leads to either ceramic or metal phase domination in terms of the volume fraction in the FG structure. The grading parameter shows the balance or linear mixing ratio between the constituents. The mixing ratio of the composed materials can be interpreted in terms of the percentage according to the following expression: where and represent the top and the bottom surface material volume fraction, respectively.
In order to determine the stress distribution within the FG plate, it is necessary to analyze the deformation that takes place in any location of the structure. In this study, the HA/Ti FG plate geometry is meshed up using a 20node brick element to divide the structure into finite elements. The main reason for choosing this element type is that it can analyze curved boundaries using a piecewisequadratic approximation. The shape functions for the 20node quadratic element can be expressed as follows:
In these equations, , , and refer to the values of the local coordinates , , and at nodes. refers to the matrix of the shape functions, which can be written as follows:
The polynomials of these coordinates are the shape function used for the interpolation of the geometry and displacement field. The isoparametric element of the 20node brick type in the simulation of this study is shown in Figure 4. It is indicated that the nodes of this element type are located at the vertices and in the middle of the edges.
For 3D isoparametric elements including the 20node brick type, the displacement and the geometry field are related by interpolation using the same functions such as the following: where , , and are point coordinates, whereas , , and are coordinates of nodes: where , , and are the displacements at points with local coordinates , , and , and , , and are displacement values at nodes.
Because FGMs are commonly used in the application where high temperature gradients are encountered, basic knowledge on the temperature distribution within the structure must be known. The first part of the present simulation is devoted to the phenomenon of thermal conduction, which shows the temperature flow from the higher temperature surface to the lower temperature surface. The relationship between the heat flux vector and the temperature distribution is explained by the Fourier law in the following expression: where is the heat flux vector, is the temperature, and is the thermal conductivity matrix. As the temperature decreases from the higher temperature surface to the lower temperature surface, the temperature variation is only available on the plane of the FG plate. The and planes show homogeneous temperature because there is no temperature change along these directions. The onedimensional differential equation governed from the steadystate heat conduction of the FGM without heat source analysis performed in this study can be written as follows:
The temperature gradient along the thickness direction of the FGM plate constructed from the thermal conductivity values calculated based on the linear rule of mixtures is expressed as follows: where where and refer to the thermal conductivity of the upper and the lower surfaces of the plate, respectively. It should be noted that, in the case of homogeneous plates, the thermal conductivity does not depend on the thickness. This is why, in the pure metal and ceramic plates, there is only a similar temperature distribution that can be found.
It is obvious that the properties of each layer of the FGM must in some way represent an average of the properties of their individual components. The property averaging of the idealized model of the relatively isotropic nature of HA/Ti FG plates can be illustrated by considering the modulus of elasticity of the composed materials. When the matrix of one of the composed materials is ultimately bonded to another material, the strain of both matrixes will be the same. This isostrain condition is valid, although the elastic modulus of both materials will tend to be quite different. In other words, where and represent the strain and the stress, respectively, while is the Young modulus of the material. It is apparent that the load carried by the FG plate, , is the simple sum of loads carried by each component: where and refer to the load carried by ceramic and metal, respectively. Each load is equal, by definition, to a stress times crosssectional area, which is
is the area of the FGM plate, is the area of the ceramic, and is the area of the metal. Equally significant for the relative contribution of the ceramic phase to the FGM modulus is the fraction of the total FGM load, , in (13), carried by the axially loaded ceramic. Consider the following: where refers to the volume fraction of the ceramic phase. In most cases, higher modulus and strength of the ceramic are effectively transmitted to the FGM as a whole. As a result, the ductility of the matrix can produce a substantially less brittle material than the metal by itself. The concept of the thermomechanical analysis is based on the minimum potential energy principle, which states that where [] and [] are the element shear stiffness and bending matrices, respectively, while and are the element mechanical and thermal load vectors, respectively. Each of them can be defined by the following equations:
From the equations, [] and [] refer to the curvature and shear strain displacement matrices, respectively. and are the shear and bending material Young modulus, is the material matrices, and is the mechanical load given on the element. The thermal load vectors, which are represented by , are classified into two terms known as and , which are considered within the element shear stiffness and bending, respectively. The definition of the normal stress () and strain () of the structure is based on the following equations:
3. Finite Element Modeling
The numerical simulations of thermomechanical problems in FG plates are carried out by using the generalpurpose finite element software package ANSYS. The simulations are divided into three parts, where the first part shows the temperature field calculation using a thermal analysis preference. The second part of the simulation demonstrates the structural analysis in which the displacement and the stress fields are calculated, while the last part combines both analyses to a thermalstructural coupled analysis to consider the thermomechanical problems in the FG structures. The typical thermal and mechanical properties of the FGM constituents used to define the material properties for the calculation are listed in Table 1, where , , , and refer to thermal conductivity, Young modulus, Poisson ratio, and thermal expansion coefficient (CTE), respectively.
The model of the FG plate is designed to be in a 3D rectangular shape with eight overlapped layers. Because of the exponentially changed properties along the thickness of the FG plate, different properties calculated using the rule of mixture formulation are defined for each layer. However, because the effect of the variation of the Poisson ratio is negligible, this property is taken as a constant value upon which all the calculations are based on Apalak and Gunes [24]. The model is discretized into 400 elements and 3,003 nodes using a hexagonal mapped mesh function in mesh tool operation. The rectangular FG plate modeled throughout this FEM work is subjected to simply supported boundary condition where all nodes along direction edges are fixed in  and directions and all nodes along the direction edges are fixed in  and directions. Figures 5 and 6 display the geometric configuration and the property definition of the FG plate model, respectively.
From Figure 6, the bottom (material number 1) and the top layers (material number 8) have the pure ceramic and metal properties, respectively, whereas the intermediate layers (material numbers 2 to 7) possess the gradually changed material properties from the different stepwise compositions of both the ceramic and the metal materials. Considering the real service condition where the ceramic phase is typically taken as a thermal barrier on the metalbased structural components in various heatresistant coatings, the top and the bottom surfaces of the FG plate are assigned to be ceramic and metalrich layers, respectively. The concept of the thermomechanical analysis can be understood clearly using the illustration given in Figure 7.
It is indicated that, for the thermomechanical problem simulation, the FG plate is subjected to either a constant or a functional mechanical transverse load on the ceramicrich top surface in addition to a thermal load where the higher temperature is located at the top surface, whereas the lower temperature is located at the bottom.
For the simulation of the temperature field, the thermal conductivity of the plate can be input directly in the material model as thermal parameter. For constant loading, the high temperature of 300°C is loaded at the top surface of the FG plate as a thermal load while the bottom surface is maintained at 20°C of room temperature. The functional parameter of the sinusoidal thermal loading is defined as follows: where is completing as 300°C, is the coordinate at the plane, and is the coordinate at the plane. After finishing the simulation of the temperature field, the element is switched to its corresponding structural element in order to get the structural response of the FGMs. The elastic modulus, Poisson ratio, and the coefficient of thermal expansion are input directly in the material model as structural parameters. The steadystate simulation results of the temperature field are used in the structural response analyses as the thermal loads. The pressure load is applied to the top surface of the FG plate as a mechanical surface load. The function of the sinusoidal mechanical loading is defined as follows: where is the load intensity given to the FG plate.
4. Results and Discussion
Numerical results are presented in this section for simply supported HA/Ti FG rectangular plates with various grading indices. The FEM calculations have been performed by taking into consideration the loading parameters from thermal, mechanical, and coupled thermomechanical loads to observe the effects of different loading conditions on the response of the FG plates. In order to ensure that the requirements for the simulation model have been transferred to the computer model with sufficient accuracy, the first part of this section shows the verification of the FEM calculations through the comparison between the present results and the established research results. The results on the temperature, the displacement, and the stress fields of HA/Ti FG plates subjected to thermal, mechanical, and thermomechanical loadings are given and discussed deeply in the subsequent sections with the intention of using the FGM in thermal barrier applications.
4.1. Modeling Verification
Verification of the FEM model is addressed through three steps. In the early stage of the verification process, the reliability of the steadystate thermal analysis method is confirmed because the present values of temperature distributions in Al/ZrO_{2} shown in Figure 8 are highly consistent with those presented in Alshorbagy et al. [29]. It is indicated from the thermal analysis that the homogeneous plates’ thermal profiles exhibited higher temperature distribution than those that belong to the FG plates. The dependence of the temperature field on the thermal profile at the surfaces and the corresponding thermal conductivity properties of the constituents exhibit that the grading parameter has a dominant role on the thermal characteristic of the FG plates.
The result of the thermal analysis is defined as the thermal loads for thermalstructural coupled analysis in which the calculation of thermal stress is performed. The convergence studies of simply supported Monel/Zirconia FG plates subjected to thermal loads with respect to various mesh and layer numbers are summarized in Tables 2 and 3, respectively. The results in Table 2 indicate that as the mesh division increases, the FG plates’ response properties’ values converge at the established research results. Table 2 indicates that the increased layer number leads to better thermal stress relaxation because of the smaller mismatch in material properties of the constituents. The results from the present calculation are in excellent agreement with those reported in Srinivas et al. [30] and Reddy and Cheng [28].


The deflections of the Ni/ZrO_{2} FG plates subjected to mechanical loads calculated in the present simulation are shown in Table 4. The displacement and the stress field results are obtained using structural analysis that considers the typical properties such as the elastic modulus, the CTE, and the Poisson ratio. The simulation method is verified because the present values are highly consistent with that reported in Alshorbagy et al. [29].

The thermalstructural coupled model that was developed in order to calculate the thermoelastic behaviors of the FG plates subjected to thermomechanical loadings is verified in this section because the present results are in excellent agreement with those presented in Alshorbagy et al. [29], which are shown in Table 5. In response to the simulation method verification, the FEM formulation used in this study is considered valid for determining the solutions for thermomechanical problems in FGMs.

4.2. Thermoelastic Properties of HA/Ti FGM
The FEM simulation verified in the previous section is used for the analysis of thermomechanical problems in HA/Ti FG plates with various grading indices, . The results of the temperature profile, the displacement, and the stress distribution of HA/Ti FG plates imposed to thermal, mechanical, and coupled thermomechanical loadings simulated from the present FE analysis are given and discussed in this section. The structural responses of the FG plates are given in terms of the following nondimensional parameters: nondimensional load intensity, , nondimensional normal stress, , nondimensional transverse shear stress, , nondimensional central deflection, , nondimensional thickness, , nondimensional coordinate at direction, ,where , , , and are referring to the pressure load, plate length, plate thickness, and Young modulus of the bottom surface material, respectively.
4.2.1. Deformation by Thermal Loadings
In the steadystate heat conduction analysis, the HA/Ti FG plate is subjected to a temperature gradient along the thickness where the temperature of the HA top surface is 300°C and the temperature of the Ti bottom surface is 20°C. Figure 9 depicts the temperature distribution profile achieved from the FE simulation for the HA/Ti FG plates. It shows the temperature contour for heat conduction from the higher temperature 300°C surfaces to the lower temperature 20°C surfaces in FG plates with various grading indices (). From the analysis, it is observed that the temperature distribution within the homogeneous plates () is usually larger than those of FG plates (). The dependence of the temperature profile on the thermal conductivity of the constituents’ mismatch, which thus affects the thermal loading, shows the domination of the grading parameter on the thermal field of the HA/Ti FG plates. A comparison between the temperature profile of FG plates subjected to constant loading and that subjected to functional loading shows no significant difference as shown in Figure 10.
(a)
(b)
(c)
(d)
(a) Constant loadings
(b) Sinusoidal loadings
The grading number, , is taken as the optimum grading index for HA/Ti FG plates because it leads to the maximum thermal relaxation with the lowest temperature distribution within the FG plate. The temperature profile obtained from the thermal analysis is transferred to the structural analysis as the thermal loading. The deflection of the HA/Ti FG plates with various grading indices under thermal excitation is shown in Figure 11.
(a) Constant loadings
(b) Sinusoidal loadings
It is indicated that the deflection in FG plates is in intermediate layers between the full ceramic and the full metal plates. The elongation of the FG plates subjected to thermal loading depends on the CTE properties. In this study, it is known that the CTE of HA is higher than that of Ti. As the volume fraction of the ceramic phase increases within the smaller grading index number, the properties of the ceramic will dominate the deformation behavior of the FG plates. The results shown in Figure 11 are in excellent agreement with the theory and thus lead to the conclusion that the grading parameter has a dominant effect on the structural response of the HA/Ti FG plates subjected to thermal loading.
The normal stress field distribution in HA/Ti FG plates induced by the temperature profile shown in Figures 10 and 11 is shown in Figure 12. The smaller area of maximum normal stress contours for FG plates () illustrated in Figure 12 is very advantageous to show the thermal stress relaxation in FGMs. The similar normal stress contour distribution of full ceramic and full metal plates () shows that the variation of thermal conductivity is negligible for homogeneous plates with heat conduction calculation. The nondimensional normal stress distribution along the thickness of HA/Ti FG plates subjected to thermal loading is given in Figure 13. The maximum tensile stress is located at the Ti bottom surface, while the maximum compressive stress appears at the HA top surface of the FG plates. Both stress components are lower than that of homogeneous plates. The results reveal that grading parameter, , does affect the response of FG plates under thermal loadings.
(a)
(b)
(c)
(d)
(a) Constant loadings
(b) Sinusoidal loadings
Figure 14 depicts the transverse shear stress field distribution in HA/Ti FG plates with various grading parameters under thermal loadings. The different shear stress contours between homogeneous and FG plates lead to the conclusion that the properties’ gradation in FGMs has impressed on the transverse shear stress response of FGMs subjected to thermal loading. The nondimensional transverse shear stress distribution in the middle plane for various grading parameters of HA/Ti FG plates is indicated in Figure 15.
(a)
(b)
(c)
(d)
(a) Constant loadings
(b) Sinusoidal loadings
The shear stress field is significant for the area near the edges of the FG plates because there are shear stress fields observed at this location. The shear stress at the locations away from the surface of the FG plates is similar for various grading parameters, . The difference between the shear stress field distribution of homogeneous and FG plates displayed in Figure 15 shows that the grading parameter does affect the shear stress response of the FG plates under thermal loads.
4.2.2. Deformation by Mechanical Loadings
In order to calculate the elastostatic behavior of the FG plates with various grading indices, the pressure load that takes place in the interval of is applied to the HA top surface. The structural responses of FG plates subjected to mechanical loads will be shown and discussed in this section. To study the displacement response, the nondimensional central deflection of HA/Ti FG plates with various grading parameters, , subjected to mechanical loadings summarized in Table 6 can be considered.

It is noticed that the nondimensional deflection of FG plates subjected to constant loadings is higher than those subjected to sinusoidal loadings. When subjected to constant loadings, the whole area at the same plane of the plate is imposed to a uniformly distributed load. However, when imposed on sinusoidal loadings, a certain coordinate is imposed on a certain loading magnitude. As a result, the lower deflection is obtained in FG plates subjected to functional loadings. Figure 16 depicts the nondimensional central deflections of the FG plates with various grading parameters induced by the different load intensities.
(a) Constant loadings
(b) Sinusoidal loadings
It can be seen that the central deflection, , decreases with the increase in the grading index number, . It is known that the displacement response of a material strongly depends on the elastic moduli of its ingredient. The FG plates with a higher grading index number contain higher ceramic volume fractions. Because the elastic moduli of the HA are lower than that of the Ti, the high toughness property of the metal dominates the elastic properties of the FG plates. The rigidity of FG plates with higher Ti volume fraction will be higher as well. The property averaging of the relatively isotropic nature of HA/Ti FG plates leads to the intermediate deflection between the full ceramic and full metal deflection given in Figure 16.
The normal stress field in HA/Ti FG plates with various grading parameters is illustrated in Figure 17. The normal stress induced by mechanical loadings is not significant because there is no variation in the stress contours. However, the differences in normal stress values of homogeneous and FG plates are advantageous to explain the domination of the grading parameter .
(a)
(b)
(c)
(d)
4.2.3. Deformation by Thermomechanical Loading
In this section, the thermomechanical problem with HA/Ti FG plates with various grading parameters, , is considered. For the thermoelastic behavior simulation, the FG plates are imposed on a nondimensional load intensity that is the interval of simultaneously with thermal loadings. The thermoelastic displacement field of HA/Ti FG plates under various load intensity values is illustrated in Figure 18.
(a)
(b)
(c)
(d)
It is observed that when subjected to load intensity of and 2, the maximum deflection contour appears at the location near the surface of the FG plates. However, when the higher load intensity () is applied, the maximum deflection contour moves to the center location of the FG plates. These results show that the displacement response of HA/Ti FG plates subjected to thermomechanical loadings differs from those subjected to thermal or mechanical individual loading. The thermal excitation on the pressureloaded plates leads to either the shrinkage or expansion of the FG plates. The nondimensional central deflections of HA/Ti FG plates with various grading parameters, , under thermomechanical loading are summarized in Table 7.

The negative values indicated in Table 7 show the difference between the displacement response induced by the mechanical loadings and that induced by the thermomechanical loadings. The property averaging in FG plates leads to the intermediate deflection values given in Table 7. Figure 19 shows the relationship between the load intensity given on the top surface of the FG plates and the displacement response of the plates. The positive and negative values of deflection given in Table 6 lead to the intersection point shown in Figure 19.
(a) Constant loadings
(b) Sinusoidal loadings
The variation of normal stress field distribution along the thickness of HA/Ti FG plates with various grading parameters, , is shown in Figure 20. The maximum compressive stress is located at the HA top surface and increases as the ratio increases. The maximum tensile stress that is found at the Ti bottom surface of FG plates is lower than that of homogenous plates. Figure 21 displays the normal stress field distribution in FG plates with various grading parameters. The most significant aspect of the contour is that the normal stress distributions of FG plates differ from those in isotropic plates.
(a) Constant loadings
(b) Sinusoidal loadings
(a)
(b)
(c)
(d)
Figure 22 depicts the transverse shear stress distribution in HA/Ti FG plates with various grading parameters under thermomechanical loading. The different shear stress contours between homogeneous and FG plates lead to the conclusion that the properties’ gradation in FGMs has impression on the transverse shear stress response of FGMs subjected to coupled thermomechanical loadings. The nondimensional transverse shear stress distribution along the thickness and at the middle plane of HA/Ti FG plates with various grading parameters subjected to coupled thermomechanical loading is indicated in Figures 23 and 24, respectively.
(a)
(b)
(c)
(d)
(a) Constant loadings
(b) Sinusoidal loadings
(a) Constant loadings
(b) Sinusoidal loadings
It is obvious that the transverse shear stress distribution in the middle plane is more consistent than that along the thickness of the FG plates. This is why the thickness parameter of the FGM layer is often considered during the design of an FGM structure. In this study, the stepwise built connection type of the FG plate layers leads to significant shear stress distribution, which consists of both negative and positive values. The larger layer number is required in order to minimize the shear stress field in the FG plates.
5. Conclusion
FEM modeling and analysis of thermomechanical problems of a rectangular HA/Ti FG plate are presented. The effects of different grading indices of constituents on displacement and stress responses of HA/Ti FG plates under thermal, mechanical, and thermomechanical loads have been investigated. From the numerical results obtained, it can be concluded that the present analysis is adequate to observe some relevant fundamental characteristics of temperature profile and stress distributions. Moreover, the following conclusions can be made from the study.(i)Different loading conditions (thermal, mechanical, or coupled thermomechanical) lead to the different temperature profile, deflection, and stress distribution in HA/Ti FG plates.(ii)Different values of grading indices show a considerable effect on temperature profile, displacement, and stress fields of HA/Ti FG plates.(iii)The temperature distribution in FG plates is smaller than that of homogeneous plates.(iv)Deflection in FG plates subjected to sinusoidal loadings is smaller than those subjected to constant loadings.(v)The properties’ gradation leads to the relaxation in thermal and thermomechanical stress in FG plates.
This study revealed that HA/Ti FGM has a potential application in high thermal resistance materials because it can withstand the extreme surface temperature and the thermal stress.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This paper was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, under Grant no. 135008D1434. The authors, therefore, acknowledge and thank DSR’s technical and financial support. The authors from King Abdulaziz University would like to thank Universiti Malaysia Pahang (UMP) and the team members of the FGMs Research Group, UMP, for their support and for conducting the above research work.
References
 B. Kieback, A. Neubrand, and H. Riedel, “Processing techniques for functionally graded materials,” Materials Science and Engineering A, vol. 362, no. 12, pp. 81–105, 2003. View at: Publisher Site  Google Scholar
 H. H. Ibrahim, H. H. Yoo, and K.S. Lee, “Supersonic flutter of functionally graded panels subject to acoustic and thermal loads,” Journal of Aircraft, vol. 46, no. 2, pp. 593–600, 2009. View at: Publisher Site  Google Scholar
 A. Hosseinzadeh and M. T. Ahmadian, “Application of piezoelectric and functionally graded materials in designing electrostatically actuated micro switches,” Journal of Solid Mechanics, vol. 2, no. 2, pp. 179–189, 2010. View at: Google Scholar
 G. Mattei, A. Tirella, and A. Ahluwalia, “Functionally graded materials (FGMs) with predictable and controlled gradient profiles: computational modeling and realization,” Computer Modeling in Engineering and Sciences, vol. 87, no. 6, pp. 483–504, 2012. View at: Google Scholar
 L.S. Liu, Q.J. Zhang, and P.C. Zhai, “The optimization design of metal/ceramic FGM armor with neural net and conjugate gradient method,” Materials Science Forum, vol. 423–425, pp. 791–796, 2003. View at: Google Scholar
 S. N. S. Jamaludin, F. Mustapha, D. M. Nuruzzaman, T. M. Y. T. Ya, and S. Basri, “Finite element calculation of residual thermal stresses for functionally graded hydroxyapatitetitanium plate design,” Academic Platform Journal of Engineering and Science, vol. 1, no. 2, pp. 1–10, 2013. View at: Google Scholar
 S. N. S. Jamaludin, F. Mustapha, D. M. Nuruzzaman, and S. Basri, “A review on the fabrication techniques of functionally graded ceramicmetallic materials in advanced composites,” Scientific Research and Essays, vol. 8, pp. 828–840, 2013. View at: Google Scholar
 H.S. Shen and N. Noda, “Postbuckling of FGM cylindrical shells under combined axial and radial mechanical loads in thermal environments,” International Journal of Solids and Structures, vol. 42, no. 1617, pp. 4641–4662, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 H. Matsunaga, “Stress analysis of functionally graded plates subjected to thermal and mechanical loadings,” Composite Structures, vol. 87, no. 4, pp. 344–357, 2009. View at: Publisher Site  Google Scholar
 A. Alibeigloo, “Exact solution for thermoelastic response of functionally graded rectangular plates,” Composite Structures, vol. 92, no. 1, pp. 113–121, 2010. View at: Publisher Site  Google Scholar
 Y. Kiani and M. R. Eslami, “Thermomechanical buckling of temperaturedependent FGM beams,” Latin American Journal of Solids and Structures, vol. 10, no. 2, pp. 223–246, 2013. View at: Google Scholar
 P. Zhu, L. W. Zhang, and K. M. Liew, “Geometrically nonlinear thermomechanical analysis of moderately thick functionally graded plates using a local PetrovGalerkin approach with moving Kriging interpolation,” Composite Structures, vol. 107, pp. 298–314, 2014. View at: Google Scholar
 L. W. Zhang, P. Zhu, and K. M. Liew, “Thermal buckling of functionally graded plates using a local Kringing meshless method,” Composite Structures, vol. 108, pp. 472–492, 2014. View at: Google Scholar
 K. M. Liew, Z. X. Lei, J. L. Yu, and L. W. Zhang, “Postbuckling of carbon nanotubereinforced functionally graded cylindrical panels under axial compression using meshless approach,” Computer Methods in Applied Mechanics and Engineering, vol. 268, pp. 1–17, 2014. View at: Google Scholar
 L. W. Zhang, Z. X. Lei, K. M. Liew, and J. L. Yu, “Static and dynamic of carbon nanotube reinforced functionally graded cylindrical panels,” Composite Structures, vol. 111, pp. 205–212, 2014. View at: Google Scholar
 T. Nomura, H. Moriguchi, K. Tsuda, K. Isobe, A. Ikegaya, and K. Moriyama, “Material design method for the functionally graded cemented carbide tool,” International Journal of Refractory Metals and Hard Materials, vol. 17, no. 6, pp. 397–404, 1999. View at: Publisher Site  Google Scholar
 E. Carrera, S. Brischetto, and A. Robaldo, “Variable kinematic model for the analysis of functionally graded material plates,” AIAA Journal, vol. 46, no. 1, pp. 194–203, 2008. View at: Publisher Site  Google Scholar
 Y. Miyamoto, W. A. Kaysser, B. H. Rabin, A. Kawasaki, and R. G. Ford, Functionally Graded Materials: Design, Processing and Applications, Kluwer Academic Publishers, London, UK, 1999.
 V. Birman and L. W. Byrd, “Modeling and analysis of functionally graded materials and structures,” Applied Mechanics Reviews, vol. 60, no. 1–6, pp. 195–216, 2007. View at: Publisher Site  Google Scholar
 J. N. Reddy, “Analysis of functionally graded plates,” International Journal for Numerical Methods in Engineering, vol. 47, no. 1–3, pp. 663–684, 2000. View at: Google Scholar  Zentralblatt MATH
 A. Ghosh and I. Reimanis, Functionally Graded Materials: Manufacture, Properties and Applications, vol. 76 of Ceramic Transactions, American Ceramic Society, 1997.
 C. S. Lee, S.H. Ahn, L. C. DeJonghe, and G. Thomas, “Effect of functionally graded material (FGM) layers on the residual stress of polytypoidally joined Si_{3}N_{4}Al_{2}O_{3},” Materials Science and Engineering A, vol. 434, no. 12, pp. 160–165, 2006. View at: Publisher Site  Google Scholar
 C. S. Lee, S.G. Kim, S.H. Ahn, L. C. DeJonghe, and G. Thomas, “Three dimensional analysis of thermal stress and prediction of failure of polytypoidally joined Si_{3}N_{4}Al_{2}O_{3} Functionally Graded Material (FGM),” Materials Transactions, vol. 48, no. 9, pp. 2489–2493, 2007. View at: Publisher Site  Google Scholar
 M. K. Apalak and R. Gunes, “Thermal residual stress analysis of NiAl_{2}O_{3}, NiTiO_{2}, and TiSiC functionally graded composite plates subjected to various thermal fields,” Journal of Thermoplastic Composite Materials, vol. 18, no. 2, pp. 119–152, 2005. View at: Publisher Site  Google Scholar
 C. Chu, J. Zhu, Z. Yin, and P. Lin, “Optimal design and fabrication of hydroxyapatiteTi asymmetrical functionally graded biomaterial,” Materials Science and Engineering A, vol. 348, no. 12, pp. 244–250, 2003. View at: Publisher Site  Google Scholar
 H. H. Moroi, K. Okimoto, R. Moroi, and Y. Terada, “Numeric approach to the biomechanical analysis of thermal effects in coated implants,” The International Journal of Prosthodontics, vol. 6, no. 6, pp. 564–572, 1993. View at: Google Scholar
 G. Ruseska, E. Fidancevska, and J. Bossert, “Mechanical and thermalexpansion characteristics of Ca_{10}(PO_{4})_{6}(OH)_{2}Ca_{3}(PO_{4})_{2} composites,” Science of Sintering, vol. 38, no. 3, pp. 245–253, 2006. View at: Publisher Site  Google Scholar
 J. N. Reddy and Z.Q. Cheng, “Threedimensional thermomechanical deformations of functionally graded rectangular plates,” European Journal of Mechanics—A/Solids, vol. 20, no. 5, pp. 841–855, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. E. Alshorbagy, S. S. Alieldin, M. Shaat, and F. F. Mahmoud, “Finite element analysis of the deformation of functionally graded plates under thermomechanical loads,” Mathematical Problems in Engineering, vol. 2013, Article ID 569781, 13 pages, 2013. View at: Publisher Site  Google Scholar
 G. Srinivas, U. S. Prasad, M. Manikandan, and A. P. Kumar, “Simulation of traditional composites under thermal loads,” Research Journal of Recent Sciences, vol. 2, pp. 273–278, 2013. View at: Google Scholar
Copyright
Copyright © 2014 S. N. S. Jamaludin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.