Anisotropic Behavior of Different Three-Dimensional Structures Materials under Thermal Stress Effect

Abstract

This work concerns the numerical modeling of stationary conduction heat transfer in a 3D three-dimensional anisotropic material subjected to an internal heat source, based on the finite element method MEF and using the Galerkin method. The field of study is a cube representing the seven crystalline systems subjected to an internal heat source and convective boundaries. The obtained equation system is solved by the LU method. The automatic mesh is managed for all the domain nodes via the program which we have written in FORTRAN language. This program allowed temperature field calculation and was applied for different crystalline systems: monoclinic, triclinic, orthorhombic, trigonal, cubic that are identified by their thermal conductivity tensors [kij]. The obtained temperature profiles obtained are in accordance with heat transfer theory and clearly illustrate the crystalline structure symmetry; this calculation permits to predict the possible thermal deformations in an anisotropic solid.