Abstract
The principal objective of this chapter is to introduce a new fractional-order theory for functionally graded anisotropic (FGA) structures. This theory called nonlinear uncoupled magneto-thermoelasticity theory involving three temperatures. Because of strong nonlinearity, it is very difficult to solve the problems related to this theory analytically. Therefore, it is necessary to develop new numerical methods for solving such problems. So, we propose a new boundary element model for the solution of general and complex problems associated with this theory. The numerical results are presented graphically in order to display the effect of the graded parameter on the temperatures and displacements. The numerical results also confirm the validity and accuracy of our proposed model.
Keywords
- boundary element method
- fractional-order heat transfer
- functionally graded anisotropic structures
- nonlinear uncoupled magneto-thermoelasticity
- three temperatures
1. Introduction
Functionally graded material (FGM) is a special type of advanced inhomogeneous materials. Functionally graded structure is a mixture of two or more distinct materials (usually heat-resisting ceramic on the outside surface and fracture-resisting metal on the inside surface) that have specified properties in specified direction of the structure to achieve a require function [1, 2]. This feature enables obtaining structures with the best of both material’s properties, and suitable for applications requiring high thermal resistance and high mechanical strength [3, 4, 5, 6, 7, 8, 9, 10, 11, 12].
Functionally Graded Materials have been wide range of thermoelastic applications in several fields, for example, the water-cooling model of a fusion reactor divertor is one of the most widely used models in industrial design, which is consisting of a tungsten (W) and a copper (Cu), that subjected to a structural integrity issue due to thermal stresses resulted from thermal expansion mismatch between the bond materials. Recently, functionally graded tungsten (W)–copper (Cu) has been developed by using a precipitation-hardened copper alloy as matrix instead of pure copper, to overcome the loss of strength due to the softening of the copper matrix.
The carbon nanotubes (CNT) in FGM have new applications such as reinforced functionally graded piezoelectric actuators, reinforced functionally graded polyestercalcium phosphate materials for bone replacement, reinforced functionally graded tools and dies for reduce scrap, better wear resistance, better thermal management, and improved process productivity, reinforced metal matrix functionally graded composites used in mining, geothermal drilling, cutting tools, drills and machining of wear resistant materials. Also, they used as furnace liners and thermal shielding elements in microelectronics.
There are many areas of application for elastic and thermoelastic functionally graded materials, for example, industrial applications such as MRI scanner cryogenic tubes, eyeglass frames, musical instruments, pressure vessels, fuel tanks, cutting tool inserts, laptop cases, wind turbine blades, firefighting air bottles, drilling motor shaft, X-ray tables, helmets and aircraft structures. Automobiles applications such as combustion chambers, engine cylinder liners, leaf springs, diesel engine pistons, shock absorbers, flywheels, drive shafts and racing car brakes. Aerospace applications rocket nozzle, heat exchange panels, spacecraft truss structure, reflectors, solar panels, camera housing, turbine wheels and Space shuttle. Submarine applications such as propulsion shaft, cylindrical pressure hull, sonar domes, diving cylinders and composite piping system. Biotechnology applications such as functional gradient nanohydroxyapatite reinforced polyvinyl alcohol gel biocomposites. Defense applications such as armor plates and bullet-proof vests. High-temperature environment applications such as aerospace and space vehicles. Biomedical applications such as orthopedic applications for teeth and bone replacement. Energy applications such as energy conversion devices and as thermoelectric converter for energy conservation. They also provide thermal barrier and are used as protective coating on turbine blades in gas turbine engine. Marine applications such as parallelogram slabs in buildings and bridges, swept wings of aircrafts and ship hulls. Optoelectronic applications such as automobile engine components, cutting tool insert coating, nuclear reactor components, turbine blade, tribology, sensors, heat exchanger, fire retardant doors, etc.
According to continuous and smooth variation of FGM properties throughout in depth, there are many laws to describe the behavior of FGM such as index [13], sigmoid law [14], exponential law [15] and power law [16, 17, 18, 19, 20, 21, 22, 23, 24].
There was widespread interest in functionally graded materials, which has developed a lot of analytical methods for analysis of elasticity [25, 26, 27, 28, 29, 30, 31, 32] and thermoelasticity [33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53] problems, some of which have become dominant in scientific literature. For the numerical methods, the isogeometric finite element method (FEM) has been used by Valizadeh et al. [54] for static characteristics of FGM and by Bhardwaj et al. [55] for solving crack problem of FGM. Nowadays, the boundary element method is a simple, efficient and powerful numerical tool which provides an excellent alternative to the finite element method for the solution of FGM problems, Sladek et al. [56, 57, 58] have been developed BEM formulation for transient thermal problems in FGMs. Gao et al. [59] developed fracture analysis of functionally graded materials by a BEM. Fahmy [60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72] developed BEM to solve elastic, thermoelastic and biomechanic problems in anisotropic functionally graded structures. Further details on the BEM are given in [73, 74] and the references therein.
In the present paper, we propose new FGA structures theory and new boundary element technique for modeling problems of nonlinear uncoupled magneto-thermoelasticity involving three temperatures. The boundary element method reduces the dimension of the problem, therefore, we obtain a reduction of numerical approximation, linear equations system and input data. Since there is strong nonlinearity in the proposed theory and its related problems. So, we develop new boundary element technique for modeling such problems. The numerical results are presented graphically through the thickness of the homogeneous and functionally graded structures to show the effect of graded parameter on the temperatures and displacements. The numerical results demonstrate the validity and accuracy of our proposed model.
A brief summary of the chapter is as follows: Section 1 outlines the background and provides the readers with the necessary information to books and articles for a better understanding of mechanical behaviour of magneto-thermoelastic FGA structures and their applications. Section 2 describes the formulation of the new theory and its related problems. Section 3 discusses the implementation of the new BEM for solving the nonlinear radiative heat conduction equation, to obtain the three temperature fields. Section 4 studies the development of new BEM and its implementation for solving the move equation based on the known three temperature fields, to obtain the displacement field. Section 5 presents the new numerical results that describe the through-thickness mechanical behaviour of homogeneous and functionally graded structures.
2. Formulation of the problem
We consider a Cartesian coordinate system for 2D structure (see Figure 1) which is functionally graded along the

Figure 1.
Geometry of the FGA structure.
The fractional-order governing equations of three temperatures nonlinear uncoupled magneto-thermoelasticity in FGA structures can be written as follows [6].
where
The nonlinear time-dependent two dimensions three temperature (2D-3 T) radiation diffusion equations coupled by electron, ion and phonon temperatures may be written as follows
where
and
The total energy per unit mass can be expressed as follows
where
Initial and boundary conditions can be written as
3. BEM numerical implementation for temperature field
This section outlines the solution of 2D nonlinear time-dependent three temperatures (electron, ion and phonon) radiation diffusion equations using a boundary element method.
Now, let us consider
By using a finite difference scheme of Caputo fractional time derivative of order
Where
According to Eq. (18), the fractional order heat Eq. (4) can be replaced by the following system
According to Fahmy [60], and using the fundamental solution which satisfies the system (21), the boundary integral equations corresponding to nonlinear three temperature heat conduction-radiation equations can be written as
which can be written in the absence of internal heat sources as follows
Time temperature derivative can be written as
where
We suppose that
Then, Eq. (23) yields the following boundary integral equation
where
and
In which the entries of
Using the standard boundary element discretization scheme for Eq. (26) and using Eq. (28), we have
The diffusion matrix can be defined as
with
In order to solve Eq. (30) numerically the functions
where
By differentiating Eq. (34) with respect to time we get
The substitution of Eqs. (34)–(36) into Eq. (30) leads to
By using initial and boundary conditions, we get
This system yields the temperature, that can be used to solve (1) for the displacement.
4. BEM numerical implementation for displacement field
Based on Eqs. (2) and (3), Eq. (1) can be rewritten as
where
when the temperatures are known, the displacement can be computed by solving (39) using BEM. By choosing
The first term in (41) can be integrated partially using Gau β theory yields
The last term in (42) can be integrated partially twice using Gau β theory yields
Based on Eq. (43), Eq. (42) can be rewritten as
which can be written as
The boundary tractions are
By using the symmetry relation of elasticity tensor, we obtain
Using Eqs. (46)–(48), the Eq. (45) can be reexpressed as
We define the fundamental solution
By modifying the weighting functions, Eq. (49) can be written as
From (39), (50) and (51), the representation formula may be written as
Let
The displacement particular solution may be defined as
Differentiation of (54) leads to.
Now, we obtain the traction particular solution
The domain integral may be approximated as follows
The use of (57) together with the dual reciprocity
Leads to
From (50), we can write
By using (52), (59) and (60), we obtain
According to Fahmy [9, 10, 11], the right-hand side integrals of (61) can be reexpressed as
and
According to Fahmy [12], Guiggiani and Gigante [77] and Mantič [78] Eqs. (62) and (63) can respectively be expressed as
By using (64) and (65), the dual reciprocity boundary integral equation becomes
On the basis of isoparametric concept, we can write
By implementing the point collocation procedure and using (67) and (68), Eq. (66) may be reexpressed as
Let us suppose that
We can write (69) as follows
By using the point collocation procedure,
Now, from (74), we may derive
where
By considering the following known
Hence (76) may be written as
From the first row of (79), we can calculate the unknown fluxes
From the second row of (79) and using (80) we get
where
Eq. (81) can be written at
where
In order to solve (83), The implicit backward finite difference scheme has been applied based on the Houbolt’s algorithm and the following approximations
By using (85) and (86), we have from (83)
In which.
We implement the successive over-relaxation (SOR) of Golub and Van Loan [79] for solving (87) to obtain
5. Numerical results and discussion
The BEM that has been used in the current paper can be applicable to a wide variety of FGA structures problems associated with the proposed theory of three temperatures nonlinear uncoupled magneto-thermoelasticity. In order to evaluate the influence of graded parameter on the three temperatures and displacements, the numerical results are carried out and depicted graphically for homogeneous (
Figures 2–4 show the distributions of the three temperatures

Figure 2.
Variation of the electron temperature Te through the thickness coordinate x.

Figure 3.
Variation of the ion temperature Ti through the thickness coordinate x.

Figure 4.
Variation of the photon temperature Tp through the thickness coordinate x.
Figures 5 and 6 show the distributions of the displacements

Figure 5.
Variation of the displacement u1 through the thickness coordinate x.

Figure 6.
Variation of the displacement u2 through the thickness coordinate x.
Figures 7 and 8 show the distributions of the displacements

Figure 7.
Variation of the displacement u1 with time τ.

Figure 8.
Variation of the displacement u2 with time τ.
6. Conclusion
The main aim of this article is to introduce a new fractional-order theory called nonlinear uncoupled magneto-thermoelasticity theory involving three temperatures for FGA structures and new boundary element technique for solving problems related to the proposed theory. Since the nonlinear three temperatures radiative heat conduction equation is independent of the displacement field, we first determine the temperature field using the BEM, then based on the known temperature field, the displacement field is obtained by solving the move equation using the BEM. It can be seen from the numerical results that the graded parameter had a significant effect on the temperatures and displacements through the thickness of the functionally graded structures. Since there are no available results for the considered problem. So, some literatures may be considered as special cases from the considered problem based on replacing the heat conduction by three temperatures radiative heat conduction. The numerical results demonstrate the validity and accuracy of our proposed model. From the proposed BEM technique that has been performed in the present paper, it is possible to conclude that the proposed BEM should be applicable to any FGM uncoupled magneto-thermoelastic problem with three-temperature. BEM is more efficient, accurate and easy to use than FDM or FEM, because it only needs to solve the unknowns on the boundaries and BEM users need only to deal with real geometry boundaries. Also, BEM is reducing the computational cost of its solver. The present numerical results for our complex problem may provide interesting information for computer scientists, designers of new FGM materials and researchers in FGM science as well as for those working on the development of new functionally graded structures.