Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2019 / Article

Research Article | Open Access

Volume 2019 |Article ID 2541707 | 10 pages | https://doi.org/10.1155/2019/2541707

Thermal Analysis of 2D FGM Beam Subjected to Thermal Loading Using Meshless Weighted Least-Square Method

Academic Editor: Mohamed Shaat
Received06 Jan 2019
Revised18 Mar 2019
Accepted10 Apr 2019
Published30 Apr 2019

Abstract

The paper analyzed the thermal problem of the 2D FGM beam using meshless weighted least-square (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 clamped-clamped 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 Mori-Tanaka 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) [47], boundary element method (BEM) [8], or more recently developed meshless methods [923], 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, first-order shear deformation theory, high-order 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 element-free Galerkin (EFG) method [12, 16], the hp-clouds method, the meshless local Petrov-Galerkin (MLPG) method [9, 10], meshless Galerkin method using radial basis functions [11, 14, 21], the least-square method [20], and meshless point interpolation method [19]. The main advantage of the MLPG method compared with regular Galerkin-based methods is that no background mesh is used to evaluate various integrals appearing in the local weak formulation of problem, but it requires a high-order 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 temperature-dependent thermomechanical properties by the MLPG method. Andrew and Senthil [12] proposed a methodology for the two-dimensional 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 steady-state 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 two-dimensional (2D) higher-order 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 steady-state [29] and transient-state [30] heat conduction analysis of heterogeneous material using the meshless weighted least-square 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 Mori-Tanaka 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 least-squares approximation scheme which is originally developed for the smooth interpolation of irregularly distributed data.

2.1. The Moving Least-Square (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 steady-state 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 least-squares 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 steady-state heat conduction equations are written in a matrix formwhere and refer to node indices, N1, N2, and N3 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 least-squares 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 2-D 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 Mori-Tanaka 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 clamped-clamped 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 x-direction 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.


Cartesian gridMWLSFDM

10×100.15%0.35%
20×200.05%0.077%
30×300.024%0.033%

3.2. Case 2

A clamped-clamped 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 y-direction as .


Material propertyAlSiC

k(W/mK)23365
E(Pa)7e104.27e11
0.30.17
α (/K)2.34e-54.3e-6

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 Mori-Tanaka 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 x-displacements and y-displacement 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.


methodTemperature/KX displacement /mmY displacement/mm
maximumminimummaximumminimummaximumminimum

MWLS (31×15)361.8300.596-.8732.95-2.8
MWLS (61×30)361.3300.71-.723.88-3.84
COMSOL (948 triangle mesh)361.1300.639-.6393.81-3.81

For a=2, we analyzed the heat conduction and thermoelastic problem using MWLS method. Temperature field distribution, x-displacement, and y-displacement 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 y-displacement 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].


MWLS (31×15)Temperature/KTemperature/K
Q=5e5W/m3Q=5e6W/m3
maximumminimummaximumminimum

a=0361.8300910.8300
a=2425.63001558.2300
a=3436.13001660.9300

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 clamped-clamped 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 Mori-Tanaka 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 least-square 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. T2017-3) of Henan Polytechnic University. The correlative members of the projects are hereby acknowledged.

References

  1. 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
  2. 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
  3. J. R. Cho and D. Y. Ha, “Volume fraction optimization for minimizing thermal stress in Ni-Al2O3 functionally graded materials,” Materials Science and Engineering: A Structural Materials: Properties, Microstructure and Processing, vol. 334, no. 1-2, pp. 147–155, 2002. View at: Publisher Site | Google Scholar
  4. K.-S. Na and J.-H. Kim, “Thermal postbuckling investigations of functionally graded plates using 3-D finite element method,” Finite Elements in Analysis and Design, vol. 42, no. 8-9, pp. 749–756, 2006. View at: Publisher Site | Google Scholar
  5. X. Y. Kou and S. T. Tan, “A systematic approach for integrated computer-aided design and finite element analysis of functionally-graded-material objects,” Materials and Corrosion, vol. 28, no. 10, pp. 2549–2565, 2007. View at: Publisher Site | Google Scholar
  6. M. Lezgy-Nazargah, “Fully coupled thermo-mechanical analysis of bi-directional FGM beams using NURBS isogeometric finite element approach,” Aerospace Science and Technology, vol. 45, pp. 154–164, 2015. View at: Publisher Site | Google Scholar
  7. A. Najibi and R. Talebitooti, “Nonlinear transient thermo-elastic analysis of a 2D-FGM thick hollow finite length cylinder,” Composites Part B: Engineering, vol. 111, pp. 211–227, 2017. View at: Publisher Site | Google Scholar
  8. 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. 42-44, pp. 4511–4539, 2004. View at: Publisher Site | Google Scholar
  9. H. K. Ching and J. K. Chen, “Thermal stress analysis of functionally graded composites with temperature-dependent material properties,” Journal of Mechanics of Materials and Structures, vol. 2, no. 4, pp. 633–654, 2007. View at: Publisher Site | Google Scholar
  10. H. K. Ching and S. C. Yen, “Meshless local Petrov-Galerkin 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
  11. H. Wang and Q.-H. Qin, “Meshless approach for thermo-mechanical analysis of functionally graded materials,” Engineering Analysis with Boundary Elements, vol. 32, no. 9, pp. 704–712, 2008. View at: Publisher Site | Google Scholar
  12. A. J. Goupee and S. S. Vel, “Two-dimensional optimization of material composition of functionally graded materials using meshless analyses and a genetic algorithm,” Computer Methods Applied Mechanics and Engineering, vol. 195, no. 44-47, pp. 5926–5948, 2006. View at: Publisher Site | Google Scholar | MathSciNet
  13. 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
  14. D. F. Gilhooley, J. R. Xiao, R. C. Batra, M. A. McCarthy, and J. W. Gillespie Jr., “Two-dimensional 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
  15. G. Liu and Y. T. RGu, An Introduction to Meshfree Methods And Their Programming, Springer, Berlin, Germany, 2005.
  16. 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
  17. M. Mierzwiczak, W. Chen, and Z.-J. Fu, “The singular boundary method for steady-state nonlinear heat conduction problem with temperature-dependent thermal conductivity,” International Journal of Heat and Mass Transfer, vol. 91, pp. 205–217, 2015. View at: Publisher Site | Google Scholar
  18. 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
  19. V. N. Van Do and C.-H. Lee, “Nonlinear analyses of FGM plates in bending by using a modified radial point interpolation mesh-free 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
  20. 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
  21. J. Lin, S. Y. Reutskiy, and J. Lu, “A novel meshless method for fully nonlinear advection-diffusion-reaction problems to model transfer in anisotropic media,” Applied Mathematics and Computation, vol. 339, pp. 459–476, 2018. View at: Publisher Site | Google Scholar | MathSciNet
  22. J. Lin, C. Zhang, L. Sun, and J. Lu, “Simulation of seismic wave scattering by embedded cavities in an elastic half-plane 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
  23. J. Lin, C. S. Chen, C.-S. Liu, and J. Lu, “Fast simulation of multi-dimensional 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
  24. 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
  25. A. Alibeigloo, “Exact solution for thermo-elastic response of functionally graded rectangular plates,” Composite Structures, vol. 92, no. 1, pp. 113–121, 2010. View at: Publisher Site | Google Scholar
  26. 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
  27. K.-J. Sohn and J.-H. Kim, “Structural stability of functionally graded panels subjected to aero-thermal loads,” Composite Structures, vol. 82, no. 3, pp. 317–325, 2008. View at: Publisher Site | Google Scholar
  28. 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
  29. H. Zhou, Z. Liu, and B. Lu, “Heat conduction analysis of heterogeneous objects based on multi-color distance field,” Materials & Design, vol. 31, no. 7, pp. 3331–3338, 2010. View at: Publisher Site | Google Scholar
  30. H. M. Zhou, W. H. Zhou, G. Qin, Z. Y. Wang, and P. M. Ming, “Transient heat conduction analysis for distance-field-based irregular geometries using the meshless weighted least-square method,” Numerical Heat Transfer, Part B: Fundamentals, vol. 71, no. 5, pp. 456–466, 2017. View at: Publisher Site | Google Scholar
  31. X. Peng, N. Hu, H. Zheng, and H. Fukunaga, “Evaluation of mechanical properties of particulate composites with a combined self-consistent and Mori-Tanaka approach,” Mechanics of Materials, vol. 41, no. 12, pp. 1288–1297, 2009. View at: Publisher Site | Google Scholar
  32. H. Wang, Q. H. Qin, and Y. L. Kang, “A new meshless method for steady-state 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 © 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.


More related articles

461 Views | 240 Downloads | 0 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.