Research Article  Open Access
H. M. Zhou, X. M. Zhang, Z. Y. Wang, "Thermal Analysis of 2D FGM Beam Subjected to Thermal Loading Using Meshless Weighted LeastSquare Method", Mathematical Problems in Engineering, vol. 2019, Article ID 2541707, 10 pages, 2019. https://doi.org/10.1155/2019/2541707
Thermal Analysis of 2D FGM Beam Subjected to Thermal Loading Using Meshless Weighted LeastSquare Method
Abstract
The paper analyzed the thermal problem of the 2D FGM beam using meshless weighted leastsquare (MWLS) method. The MWLS as a meshless method is fully independent of mesh, and an approximate function was used to construct a series of linear equations to solve the unknown field variable, which avoided the troublesome task of numerical integration. The effectiveness and accuracy of the approach were illustrated by a clampedclamped FGM beam which was subjected with interior heat source. The volume fraction of FGM beam was assumed to be given by a simple power law distribution. The effective material properties of the FGM beam were assumed to be temperature independent and calculated by MoriTanaka method. The results showed that a good agreement was achieved between the proposed meshless method and commercial COMSOL Multiphysics.
1. Introduction
FGMs can resist high temperatures and are proficient in reducing the thermal stress and have received more attention from the researchers [1]. Most of these researches on FGMs have been restricted to heat conduction analyses, thermal stress analyses, thermal buckling analyses, thermal vibration, and optimization problem. Various numerical techniques, such as the finite difference method (FDM) [2, 3], finite element method (FEM) [4–7], boundary element method (BEM) [8], or more recently developed meshless methods [9–23], have been developed for analyzing thermal related problems and other problems. Because of the complexity of the relevant governing equation, analytical solutions are usually difficult to obtain for those arbitrary geometry and complex boundary conditions, and the exact solutions are usually obtained based on classical plate theory, firstorder shear deformation theory, highorder shear deformation theory, and so on [24, 25]. Compared with FEM, FDM, and BEM, the meshless methods are associated with a class of numerical techniques that approximate a given differential equation or a set of differential equations using global interpolations on the discrete nodes or background mesh, exhibiting the advantages of avoiding mesh generation, simple data preparation, easy postprocessing, and so forth.
Liu and Gu [15] introduced meshless methods and their programming, such as the elementfree Galerkin (EFG) method [12, 16], the hpclouds method, the meshless local PetrovGalerkin (MLPG) method [9, 10], meshless Galerkin method using radial basis functions [11, 14, 21], the leastsquare method [20], and meshless point interpolation method [19]. The main advantage of the MLPG method compared with regular Galerkinbased methods is that no background mesh is used to evaluate various integrals appearing in the local weak formulation of problem, but it requires a highorder quadrature rule to obtain converged results and thus needs much more computational effort in terms of CPU time than that for the FEM. Ching and Chen [9] proposed the thermoelastic analysis of a functionally graded composite considering temperaturedependent thermomechanical properties by the MLPG method. Andrew and Senthil [12] proposed a methodology for the twodimensional simulation and optimization of material distribution of functionally graded materials for thermomechanical processes using a genetic algorithm. Katsikadelis [13] employed the meshless analog equation method to solve the 2D elastostatic problem for inhomogeneous anisotropic. A meshless algorithm of fundamental solution coupling with radial basis functions based on analog equation theory was proposed to simulate the static thermal stress distribution in 2D FGMs [14]. Mierzwiczak et al. [17] presented the singular boundary method for steadystate nonlinear heat conduction problems. Bhavani et al. [26] solved thermoelastic equilibrium equations for a functionally graded beam to obtain the axial stress distribution. Sohn and Kim [27] analyzed static and dynamic stabilities of FG panels which are subjected to thermal and aerodynamic loads. Hiroyuki [28] presented a twodimensional (2D) higherorder deformation theory presented for the evaluation of displacements and stresses in functionally graded (FG) plates subjected to thermal and mechanical loading.
Zhou et al. presented steadystate [29] and transientstate [30] heat conduction analysis of heterogeneous material using the meshless weighted leastsquare method. In this paper the pure meshless method (MWLS) was then extended to solve problems of thermoelastic analysis for the FGM beam with interior heat source. The volume fractions of constituent materials composing the FGM beam are assumed to be given by a simple power law distribution. Material properties of the FGM are obtained by MoriTanaka method. The paper is divided as follows: firstly, we give problem description and MWLS analysis about the thermal problem. Then, in order to demonstrate the efficiency and accuracy of the proposed method, numerical implementation is given in the next section. The last section includes some conclusions.
2. MWLS Analysis of the Thermal Problem
The solution of MWLS analysis to the thermal problem is described in this section. The shape functions in MWLS analysis is a moving leastsquares approximation scheme which is originally developed for the smooth interpolation of irregularly distributed data.
2.1. The Moving LeastSquare (MLS) Approximation Scheme
Construct the local approximate function of an unknown field variable function expressed aswhere is the basis function and the quadratic basis is used in this paper. In 2D space, the number of basis function m=6; are unknown coefficients, which are solved by minimizing a weighted discrete residual given asThe minimum value of JJ may be achieved through differentiating with respect to in which are the positions of the N nodes, refers to the nodal parameter of the field variable at node , . is the weighting function and usually a compactly supported function that is only nonzero in a small neighborhood called the “support domain” of node . The exponential function is used in this study. where for circular support domain and is a constant. denotes the radius of the circular support domain.
To obtain , (3) can be rewritten in the matrix formwhere , , .
is chosen to make be the nonsingular matrix everywhere in the whole domain; however, the circular support domain must have enough neighborhood nodes. Through finding out the kkth nearest points of the evaluation point , the smallest support domain radius including these points can be obtained. The value of kk is gained by comparing some numerical examples with their analytical solutions in Zhou et al. [29, 30].
Solving from (5), coefficients can be obtainedSubstituting (6) back into (1) and removing unknown , we haveSet , where is the shape function.
2.2. Heat Conduction Analysis of FGM Object
The steadystate heat conduction equation and the thermal boundary conditions of FGM objects are as follows, respectively,where is the temperature field on a fixed domain surrounded by a closed boundary . The variable denotes the physical dimensions expressed in Cartesian coordinates, : . is the outward surface normal. The parameters and are the thermal conductivity, ambient temperature, the heat transfer coefficient, and the heat resource, respectively.
Substituting the unknown temperature field variable of (8) with the approximate function of (7), the residuals are minimized in a leastsquares manner,
We use an alternative discrete equation to avoid integration and (10) can be rewritten as
The system equations of the MWLS method for solving steadystate heat conduction equations are written in a matrix formwhere and refer to node indices, N_{1}, N_{2,} and N_{3} is the number of interior nodes in the boundary , and , respectively.
2.3. Thermoelastic Analysis of FGM Object
Consider the 2D FGM anisotropic linear elastic body defined in the domain bounded by . The governing equation and boundary condition can be written in the following form disregarding the body forces in which is the components of the Cauchy stress tensor. A comma followed by index denotes the partial differentiation with respect to coordinate of a material point. is the unit outward normal to . are the displacement components, and are the prescribed displacements on . are the prescribed traction on . and are the complementary parts of the boundary . MWLS analysis requires a discretization of the domain ; (15)(16) becomes where , , and is the number of interior nodes in the domain , in the boundary and , respectively.
Substituting of (17)~(19) with the approximate function of (7), where , , and denote shape function similar to of (7)where , ; No is the normal vector of any point.
Substituting the unknown field variable of (15)~(16) into (1), the residuals are minimized in a leastsquares manner,Similar to (12) the system equations of the MWLS method for solving thermoelastic problem is written in the following matrix formwhere where , , and are obtained by (23)~(25) and denotes the displacement of x, y, .
The thermal stresses is written in the matrix form where is the stiffness matrix for a linearly elastic, isotropic 2D solid. ɛ is the infinitesimal strain vector. in which , for plane stress with and α denoting the Young’s modulus, Poisson’s ratio, and coefficient of thermal expansion, respectively, and for plane stain.
2.4. Material Properties
Two homogenization methods are often used to evaluate the effective material properties for FGMs. One is the rule of mixtures, and the other is the micromechanical model. The former is simply a linear rule of mixtures and the effective value can be determined bywhere the volume fractions satisfy , and may be elastic modulus , bulk modulus , Poisson’s ratio , coefficient of thermal expansion a, thermal conductivity , and shear modulus μ.
The most widely used micromechanical model is MoriTanaka model [9, 31]. It is the modified rule of mixtures and the effective material properties can be defined using the following relation,where .
3. Numerical Results and Discussions
In order to demonstrate the efficiency and accuracy of the presented method, firstly, we choose a isotropic square region (Case 1) with defined boundary conditions; through heat conduction analysis the results are compared with the analytical solutions and FDM. Then a clampedclamped FGM beam (Case 2) which was subjected with interior heat source is analyzed using MWLS method.
3.1. Case 1
A 100×100m isotropic square region is shown in Figure 1. The top and bottom boundaries are insulated. The left and right boundaries are assigned a temperature of 200°C and 100°C, respectively. The spatial variation of the thermal conductivity is taken to be cubic in the xdirection as .
An analytical solution is given as
The design region is discretized as 20×20, in (4), and the material properties are a linear rule of mixtures (Eq. (32)). The temperature field distribution is computed in Eq. (12)(14) and shown in Figure 2. According to Wang et al. [32] the estimation rule is maximum of relative error% = %, in which is the analytical solution and is numerical solution. Then, our MWLS method is compared with FDM in different resolutions. The results are shown in Table 1. Obviously our method has a high precision compared with the analytical solution in spite of any resolution.

(a) 2D display
(b) Top boundary
3.2. Case 2
A clampedclamped FGM beam is shown in Figure 3, length L=1m, width D=0.5m, and thickness H=0.1m; material property is shown in Table 2, interior heat source Q=5e5W/m3; the spatial variation of the volume fraction of Al is taken to be a power law distribution in the ydirection as .

The beam is assumed to be in a state of plane strain normal to the xy plane, and the design region is discretized as 31×15 and 61×30. The effective material properties are determined by the MoriTanaka model (Eq. (33)~(36)). In order to verify the proposed computational method, we do some comparisons between the MWLS and the commercial COMSOL Multiphysics for a homogeneous material (a=0), relevant results are shown in Figure 4 and listed in Table 3. The results obtained with the two methods are in good agreement in temperature field aspect; however, the xdisplacements and ydisplacement have a little difference in 31×15 grid and in 61×30 the results are in good agreement in Table 3. From Figure 4 we also can know that our method and COMSOL Multiphysics have the same distribution trend. The maximum temperature 361.3K is in the center () of the beam.

(a) Temperature field of our method
(b) Temperature field of COMSOL
(c) xdisplacement of our method
(d) xdisplacement of COMSOL
(e) ydisplacement of our method
(f) ydisplacement of COMSOL
For a=2, we analyzed the heat conduction and thermoelastic problem using MWLS method. Temperature field distribution, xdisplacement, and ydisplacement are plotted in Figures 5, 6, and 7, respectively. Figure 5 indicates that the maximum temperature 425.6K is higher than the homogeneous material of Figure 2(a), in the Cartesian coordinates () of the beam. Figures 6 and 7 show that when subjected to temperature rise, the beam expands and the maximum ydisplacement is located at the top middle of the beam. Then we do heat conduction analysis in different material distribution in different heat source; the result is listed in Table 4. From Table 4, we can know that the maximum temperature of FGM model is higher than that of the fully metal model. Moreover, as the volume fraction index is increased, the maximum temperature increases. This is because for FGMs, when the volume fraction index is increased, the contained quantity of ceramic increases. Finally, to make a comparison, we do thermoelastic analysis and obtain thermal stresses in the neutral axis of the beam among a=0, a=2, and a=3, as shown in Figure 8. In Figures 8(a) and 8(c), the volume fraction of Al is gradually decreased from a=0 to a=3, the σx and σy stresses are in an upward trend. The maximum thermal stress always occurred in the vicinity of neutral axis of the beam from Figures 8(a), 8(c), and 8(d). The results also agreed well with the presented elasticity solutions of [26].

(a)
(b)
(c)
(d)
4. Conclusion
In this paper, a novel thermoelastic analysis of FGM beam based on MWLS method was presented. We do thermoelastic and heat conduction analysis aimed at a clampedclamped thick beam which is subjected with interior heat source. The FGM beam is assumed to be given by a simple power law distribution. Material properties of the FGM beam are obtained by MoriTanaka method. Through being compared with analytical solution and the commercial software of COMSOL Multiphysics, the effectiveness and accuracy are verified. We also listed the comparison of thermal stresses with the variation of power law index. The present method of analysis will be also useful in the design and optimization of FGM objects.
Nomenclature
Coefficient  
Matrices of computation  
:  Radius of the circular support domain 
:  The nodal parameter of the field variable at node 
Node indices  
Thermal conductivity  
Number of neighbor points  
Number of nodes  
Outward surface normal  
Shape functions  
Number of interior nodes in the boundary and  
Basis function  
Heat source  
:  Components of the Cauchy stress tensor 
:  Complementary parts of the boundary 
Stiffness matrix  
Young’s modulus  
Coefficient of thermal expansion  
Volume fraction  
Temperature  
:  Ambient temperature 
The heat transfer coefficient  
:  The moving leastsquare approximation function 
:  The positions of the nodes 
Fixed domain  
Closed boundary of  
:  Dirichlet boundary 
:  Neumann boundary 
:  Mixed boundary 
:  Weight functions 
Normal heat flux  
:  Given traction on 
:  Displacement components 
Displacement of x, y  
Infinitesimal strain vector  
Poisson’s ratio  
Shear modulus  
Bulk modulus. 
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
The work described in this paper was supported by a grant from the National Natural Science Foundation of China (Projects No. 51505131, U1504106) and Program for Innovative Research Team (No. T20173) of Henan Polytechnic University. The correlative members of the projects are hereby acknowledged.
References
 K. Swaminathan and D. Sangeetha, “Thermal analysis of FGM plates – A critical review of various modeling techniques and solution methods,” Composite Structures, vol. 160, pp. 43–60, 2017. View at: Publisher Site  Google Scholar
 S. Turteltaub, “Optimal control and optimization of functionally graded materials for thermomechanical processes,” International Journal of Solids and Structures, vol. 39, no. 12, pp. 3175–3197, 2002. View at: Publisher Site  Google Scholar
 J. R. Cho and D. Y. Ha, “Volume fraction optimization for minimizing thermal stress in NiAl2O3 functionally graded materials,” Materials Science and Engineering: A Structural Materials: Properties, Microstructure and Processing, vol. 334, no. 12, pp. 147–155, 2002. View at: Publisher Site  Google Scholar
 K.S. Na and J.H. Kim, “Thermal postbuckling investigations of functionally graded plates using 3D finite element method,” Finite Elements in Analysis and Design, vol. 42, no. 89, pp. 749–756, 2006. View at: Publisher Site  Google Scholar
 X. Y. Kou and S. T. Tan, “A systematic approach for integrated computeraided design and finite element analysis of functionallygradedmaterial objects,” Materials and Corrosion, vol. 28, no. 10, pp. 2549–2565, 2007. View at: Publisher Site  Google Scholar
 M. LezgyNazargah, “Fully coupled thermomechanical analysis of bidirectional FGM beams using NURBS isogeometric finite element approach,” Aerospace Science and Technology, vol. 45, pp. 154–164, 2015. View at: Publisher Site  Google Scholar
 A. Najibi and R. Talebitooti, “Nonlinear transient thermoelastic analysis of a 2DFGM thick hollow finite length cylinder,” Composites Part B: Engineering, vol. 111, pp. 211–227, 2017. View at: Publisher Site  Google Scholar
 A. Sutradhar and G. H. Paulino, “The simple boundary element method for transient heat conduction in functionally graded materials,” Computer Methods Applied Mechanics and Engineering, vol. 193, no. 4244, pp. 4511–4539, 2004. View at: Publisher Site  Google Scholar
 H. K. Ching and J. K. Chen, “Thermal stress analysis of functionally graded composites with temperaturedependent material properties,” Journal of Mechanics of Materials and Structures, vol. 2, no. 4, pp. 633–654, 2007. View at: Publisher Site  Google Scholar
 H. K. Ching and S. C. Yen, “Meshless local PetrovGalerkin analysis for 2D functionally graded elastic solids under mechanical and thermal loads,” Composites Part B: Engineering, vol. 36, no. 3, pp. 223–240, 2005. View at: Publisher Site  Google Scholar
 H. Wang and Q.H. Qin, “Meshless approach for thermomechanical analysis of functionally graded materials,” Engineering Analysis with Boundary Elements, vol. 32, no. 9, pp. 704–712, 2008. View at: Publisher Site  Google Scholar
 A. J. Goupee and S. S. Vel, “Twodimensional optimization of material composition of functionally graded materials using meshless analyses and a genetic algorithm,” Computer Methods Applied Mechanics and Engineering, vol. 195, no. 4447, pp. 5926–5948, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 J. T. Katsikadelis, “The 2D elastostatic problem in inhomogeneous anisotropic bodies by the meshless analog equation method (MAEM),” Engineering Analysis with Boundary Elements, vol. 32, no. 12, pp. 997–1005, 2008. View at: Publisher Site  Google Scholar
 D. F. Gilhooley, J. R. Xiao, R. C. Batra, M. A. McCarthy, and J. W. Gillespie Jr., “Twodimensional stress analysis of functionally graded solids using the MLPG method with radial basis functions,” Computational Materials Science, vol. 41, no. 4, pp. 467–481, 2008. View at: Publisher Site  Google Scholar
 G. Liu and Y. T. RGu, An Introduction to Meshfree Methods And Their Programming, Springer, Berlin, Germany, 2005.
 S. Senthil, V. J. Goupee, and A. J. Goupee, “A methodology for the optimization of material composition of functionally graded materials,” in Proceedings of the NSF DMII Grantees’ Conference, Scottsdale, Ariz, USA, 2005. View at: Google Scholar
 M. Mierzwiczak, W. Chen, and Z.J. Fu, “The singular boundary method for steadystate nonlinear heat conduction problem with temperaturedependent thermal conductivity,” International Journal of Heat and Mass Transfer, vol. 91, pp. 205–217, 2015. View at: Publisher Site  Google Scholar
 Z.J. Fu, W. Chen, and Q.H. Qin, “Boundary knot method for heat conduction in nonlinear functionally graded material,” Engineering Analysis with Boundary Elements, vol. 35, no. 5, pp. 729–734, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 V. N. Van Do and C.H. Lee, “Nonlinear analyses of FGM plates in bending by using a modified radial point interpolation meshfree method,” Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 57, pp. 1–20, 2018. View at: Publisher Site  Google Scholar  MathSciNet
 L. Sator, V. Sladek, and J. Sladek, “Bending of FGM plates under thermal load: classical thermoelasticity analysis by a meshless method,” Composites Part B: Engineering, vol. 146, pp. 176–188, 2018. View at: Publisher Site  Google Scholar
 J. Lin, S. Y. Reutskiy, and J. Lu, “A novel meshless method for fully nonlinear advectiondiffusionreaction problems to model transfer in anisotropic media,” Applied Mathematics and Computation, vol. 339, pp. 459–476, 2018. View at: Publisher Site  Google Scholar  MathSciNet
 J. Lin, C. Zhang, L. Sun, and J. Lu, “Simulation of seismic wave scattering by embedded cavities in an elastic halfplane using the novel singular boundary method,” Advances in Applied Mathematics and Mechanics, vol. 10, no. 2, pp. 322–342, 2018. View at: Publisher Site  Google Scholar  MathSciNet
 J. Lin, C. S. Chen, C.S. Liu, and J. Lu, “Fast simulation of multidimensional wave problems by the sparse scheme of the method of fundamental solutions,” Computers & Mathematics with Applications, vol. 72, no. 3, pp. 555–567, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 S. S. Vel and R. C. Batra, “Exact solution for thermoelastic deformations of functionally graded thick rectangular plates,” AIAA Journal, vol. 40, no. 7, pp. 1421–1433, 2002. 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
 B. V. Sankar and J. T. Tzeng, “Thermal stresses in functionally graded beams,” AIAA Journal, vol. 40, no. 6, pp. 1228–1232, 2002. View at: Publisher Site  Google Scholar
 K.J. Sohn and J.H. Kim, “Structural stability of functionally graded panels subjected to aerothermal loads,” Composite Structures, vol. 82, no. 3, pp. 317–325, 2008. View at: Publisher Site  Google Scholar
 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
 H. Zhou, Z. Liu, and B. Lu, “Heat conduction analysis of heterogeneous objects based on multicolor distance field,” Materials & Design, vol. 31, no. 7, pp. 3331–3338, 2010. View at: Publisher Site  Google Scholar
 H. M. Zhou, W. H. Zhou, G. Qin, Z. Y. Wang, and P. M. Ming, “Transient heat conduction analysis for distancefieldbased irregular geometries using the meshless weighted leastsquare method,” Numerical Heat Transfer, Part B: Fundamentals, vol. 71, no. 5, pp. 456–466, 2017. View at: Publisher Site  Google Scholar
 X. Peng, N. Hu, H. Zheng, and H. Fukunaga, “Evaluation of mechanical properties of particulate composites with a combined selfconsistent and MoriTanaka approach,” Mechanics of Materials, vol. 41, no. 12, pp. 1288–1297, 2009. View at: Publisher Site  Google Scholar
 H. Wang, Q. H. Qin, and Y. L. Kang, “A new meshless method for steadystate heat conduction problems in anisotropic and inhomogeneous media,” Archive of Applied Mechanics, vol. 74, no. 8, pp. 563–579, 2005. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 H. M. Zhou 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.