Method of lines for temperature field of functionally graded materials

Abstract

The finite element method (FEM) and the boundary element method (BEM) are often adopted. However, they are not convenient to spatially vary thermal properties of functionally graded material (FGM). Therefore, the method of lines (MOL) is introduced to solve the temperature field of FGM. The basic idea of the method is to semi-discretize the governing equation into a system of ordinary differential equations (ODEs) defined on discrete lines by means of the finite difference method. The temperature field of FGM can be obtained by solving the ODEs. The functions of thermal properties are directly embodied in these equations and these properties are not discretized in the domain. Thus, difficulty of FEM and BEM is overcome by the method. As a numerical example, the temperature field of a plane problem is analyzed for FGMs through varying thermal conductivity coefficient by the MOL.

Keywords

functionally graded materials / method of lines / temperature field

References

[1] SureshS, MortensenAFundamentals of Functionally Graded Materials [M], 1998, London, IOM Communications Ltd

[2] TanigawaY, AkaiT, KawamuraR, et al.. Transient heat conduction and thermal stress problems of a non-homogeneous plate with temperature-dependent material properties [J]. Journal of Thermal Stresses, 1996, 19(1): 77-102

[3] JabbariM, SohrabpourS, EslamiM R. Mechanical and thermal stresses in functionally graded hollow cylinder due to radially symmetric loads [J]. International Journal of Pressure Vessels and Piping, 2002, 79(4): 493-497

[4] ObataY, NodaN, TsujiT. Steady thermal stresses in a functionally gradient material plate[J]. Trans JSME, 1992, 58(12): 1689-1695

[5] KimK S, NodaN. Green’s function approach to unsteady thermal stresses in an infinite hollow cylinder of functionally graded material [J]. Acta Mech, 2002, 156(2): 145-161

[6] KimK S, NodaN. A green’s function approach to the deflection of a FGM plate under transient thermal loading[J]. Arch Appl Mech, 2002, 72(1): 127-137

[7] JinZ H, PaulinoG H. Transient thermal stress analysis of an edge crack in a functionally graded material [J]. Int J Fracture, 2001, 107(1): 73-98

[8] ShaoZ S. Mechanical and thermal stresses of a functionally graded circular hollow cylinder with finite length [J]. International Journal of Pressure Vessels and Piping, 2005, 82(2): 155-163

[9] FuchiyamaT, NodaN. Analysis of thermal stress in a plate of functionally gradient material [J]. JSAE Rev, 1995, 16(2): 263-268

[10] WangBao-lin, ZhenHui-tian. Application of finite element-finite difference method to the determination of transient temperature field in functionally graded materials [J]. Finite Elements in Analysis and Design, 2005, 41(3): 335-349

[11] PraveenG N, ReddyJ N. Nonlinear transient thermoelastic analysis of functionally graded ceramic-metal plates [J]. International Journal of Solids and Structures, 1998, 35(33): 4457-4476

[12] ShawR P, ManolisG D. Two dimensional heat conduction in graded materials using conformal mapping [J]. Commun Numer Methods Engrg, 2003, 19(2): 215-221

[13] SutradharA, PaulinoG H, GrayL J. Transient heat conduction in homogeneous and nonhomogeneous materials by the Laplace transform Galerkin boundary element method [J]. Engrg Anal Boundary Elem, 2002, 26(1): 119-132

[14] MatthewC W, GlaucioH P, RobertH, et al.. Stress-intensity factors for surface cracks in functionally graded material under mode-I thermomechanical loading [J]. International Journal of Solids and Structures, 2003, 41(9): 1081-1118

[15] AscherU, ChristiansenJ, RussellR D. Collocation software for boundary-value ODEs [J]. ACM Trans Mash Software, 1981, 7(2): 223-229

[16] KongXiang-qianApplication of the Finite Element Method to Heat Conduction Problems [M], 1986, Beijing, Science Press(in Chinese)