Investigation of the influence of functionally graded parameter,
The main purpose of this chapter is to propose a new boundary element formulation for the modeling and optimization of three-temperature nonlinear generalized magneto-thermoelastic functionally graded anisotropic (FGA) composite microstructures’ problems, which is the gap of this study. Numerical results show that anisotropy and the functionally graded material have great influences on the nonlinear displacement sensitivities and nonlinear thermal stress sensitivities of composite microstructure optimization problem. Since, there are no available data for comparison, except for the problems with one-temperature heat conduction model, we considered the special case of our general study based on replacing three-temperature radiative heat conductions with one-temperature heat conduction. In the considered special case, numerical results demonstrate the validity and accuracy of the proposed technique. In order to solve the optimization problem, the method of moving asymptotes (MMA) based on the bi-evolutionary structural optimization method (BESO) has been implemented. A new class of composite microstructures problems with holes or inclusions was studied. The two-phase magneto-thermoelastic composite microstructure which is studied in this chapter consists of two different FGA materials. Through this chapter, we investigated that the optimal material distribution of the composite microstructures depends strongly on the heat conduction model, functionally graded parameter, and shapes of holes or inclusions.
- boundary element method
- modeling and optimization
- nonlinear generalized magneto-thermoelasticity
- functionally graded anisotropic
- composite microstructures
In the last few years, there is significant interest in using advanced composite structures, and among the oldest examples of them, reinforced concrete, mixing concrete and steel, and plastics laminated with wood. The main benefit of the composite structures which consist of two or more different materials is that the properties of each material can be combined to form a single unit that performs better than the separate component parts. The most common form of a composite structure in construction is a steel and concrete composite, where concrete works well in pressure but has less resistance to tension. However, steel is extremely strong in tension, and when tied together, it results in a highly efficient and lightweight unit usually used for structures such as buildings and multistory bridges. Although fiberglass and carbon/epoxy composites are not yet as important as the oldest advanced composite structures in terms of tonnage or total revenue, they are very important in engineering, aerospace, transportation, bioengineering, optics, electronics, commodities, chemical plant, and energy industries, especially for the new airplanes that will concentrate on achieving major improvements in the fuel use, emissions, noise, transportation energy consumption, and other important issues to conserve the environment [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21].
Microstructure has been known to play a major role in determining the behavior of material. Therefore, material engineers strive to control the microstructure by improving their properties with the aim of producing a uniform microstructure throughout the material. They also produced FGMs whose microstructures depend on the position by treating the microstructure as a position-dependent variable; the properties of different materials can be combined into one component to achieve an optimum performance in a specific application [22, 23].
In recent years, great attention has been directed toward the study of nonlinear generalized magneto-thermoelastic interactions in functionally graded anisotropic (FGA) structures due to its many applications in physics, geophysics, earthquake engineering, astronautics, aeronautics, mining engineering, military technologies, plasma, robotics, high-energy particle accelerators, nuclear reactors, automobile industries, nuclear plants, soil dynamics, and other engineering and industrial applications. Duhamel  and Neuman  proposed the classical thermoelasticity (CTE) theory which has the following two paradoxes: first, the infinite propagation speeds of thermal signals are predicted, and second, there is no any elastic term included in heat equation. Biot  invented the classical coupled thermoelasticity (CCTE) theory to beat the first paradox in CTE, but CTE and CCTE share the second paradox. Then, numerous generalized thermoelasticity theories have been introduced to overcome the two paradoxes inherent in CTE, such as the extended thermoelasticity (ETE) theory of Lord and Shulman ; temperature-rate-dependent thermoelasticity (TRDTE) theory of Green and Lindsay ; three linear generalized thermoelasticity theories of Green and Naghdi (GN) [29, 30]; namely I, II, and III, respectively [where, GN theory I is based on Fourier’s law of heat conduction and is identical to CTE theory, GN theory II characterizes the thermoelasticity without energy dissipation (TEWOED), and GN theory III characterizes the thermoelasticity with energy dissipation (TEWED)]; dual phase-lag thermoelasticity (DPLTE) [31, 32]; and three-phase-lag thermoelasticity (TPLTE) .
A large amount of research has been done on the generalized problems of thermoelasticity [34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44]. Our interest in studying the three-temperature thermoelasticity [45, 46, 47, 48, 49] has increased due to its important low-temperature and high-temperature applications. Due to the computational difficulties, inherent in solving three-temperature nonlinear generalized magneto-thermoelastic problems of FGA composite microstructures, the problems become too complicated with no general analytical solution. Therefore, we propose a new boundary element modeling technique which has recently been successfully developed and implemented to obtain the approximate solutions for such problems. Now, the boundary element method (BEM), which is also called boundary integral equation method, has been widely adopted in a large variety of engineering and industrial applications. In the BEM, only the boundary of the solution domain needs to be discretized, so, it has a major advantage over other methods which require the whole domain discretization, such as the finite difference method (FDM) [50, 51, 52], discontinuous Galerkin method (DGM) , and finite element method (FEM) [54, 55, 56, 57]. This advantage of BEM over domain methods has significant importance for modeling of nonlinear generalized thermoelastic problems which can be implemented using BEM with little cost and less input data [58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71]. Recently, scientists were convinced that only the FEM method could solve complex engineering problems. But now after the huge achievements of the BEM and its ability to solve complex engineering problems with high efficiency, it gets them to change their conviction. Also, they tried to combine FEM and BEM in the solution of their complex problems.
The main aim of this chapter is to propose a novel boundary element formulation for modeling and optimization of three-temperature nonlinear generalized thermoelastic problems of functionally graded anisotropic (FGA) composite microstructures. The proposed boundary element technique has been
A brief summary of the chapter is as follows: Section 1 introduces an overview of the historical background for a better understanding of the nonlinear generalized magneto-thermoelastic problems and composite materials applications. Section 2 describes the physical modeling of the three-temperature nonlinear generalized thermoelastic problems of FGA composite microstructures. Section 3 outlines the BEM implementation for solving the governing equations of the considered problem to obtain the three temperatures and displacement fields. Section 4 outlines the topology optimization technique used to obtain the optimal composite microstructure with and without holes or inclusions of various shapes. Section 5 presents the new numerical results that describe the effects of anisotropy and functionally graded parameters on the problem’s fields’ sensitivities during the optimization process. Section 6 outlines the significant findings of this chapter.
2. Formulation of the problem
Consider a Cartesian coordinates system as shown in Figure 1 . We shall consider a functionally graded anisotropic composite microstructure of a finite thickness placed in a primary magnetic field acting in the direction of the -axis. The considered composite microstructure occupies the region with functionally graded material properties in the thickness direction.
where , , , , and are the mechanical stress tensor, Maxwell’s electromagnetic stress tensor, displacement vector, temperature, and reference temperature, respectively; and are, respectively, the constant elastic moduli and stress-temperature coefficients of the anisotropic medium; is the magnetic permeability; is the perturbed magnetic field; are the thermal conductivity coefficients; are specific heat coefficients; , and denote electron, ion, and phonon, respectively; is the second order tensor associated with the TEWED and TEWOED theories; is the electron-ion energy coefficient; is the electron-phonon energy coefficient; are constants; , , and are the density, time, and unified parameter which introduced to consolidate all theories into a unified equations system, respectively; , , and are the relaxation times; and is a functionally graded parameter. Also, , , , and are suitably prescribed functions; are the tractions defined by ; and and are the Kronecker delta functions.
A superposed dot denotes the differentiation with respect to the time, and a comma followed by a subscript denotes partial differentiation with respect to the corresponding coordinates.
The unit mass total energy can be written as
By using the following initial and boundary conditions:
By using the fundamental solution that satisfies the following equation:
where are singular points.
3. BEM implementation
The field equations can be written in the following operator form:
where the operators and are defined above in Eq. (23), and the operators and are defined as follows:
By applying the weighted residual method (WRM) to the differential Eq. (24), we obtain
Now, we can choose the fundamental solution as weighting function as
The corresponding traction field can be expressed as
The traction vector can be expressed as
The fundamental solution can be expressed as
By implementing the WRM and integration by parts, we can write Eq. (25) in the following form:
Based on the sifting property, we can express Eq. (34) as follows:
The generalized thermoelastic vectors and tensors can be written in contracted notation as follows:
The thermoelastic representation formula (38) can be written in contracted notation as
The vector can be splitted as
The thermoelastic representation formula (38) can also be expressed as follows:
In order to transform the domain integral in Eq. (44) to the boundary, we approximate the source vector by a series of given known functions and unknown coefficients :
Thus, the thermoelastic representation formula (44) can be expressed as
By applying the WRM to the following elastic and thermal equations:
Now, the weighting functions were chosen as the elastic and thermal fundamental solutions and .
Then, the representation formulae of elastic and thermal fields are given as follows:
The elastic and thermal representation formulae can be combined in one single equation as
By differentiation of Eq. (59) with respect to , we obtain
According to the technique of Partridge et al. , the displacements and velocities can be approximated as
where are known functions, and and are unknown coefficients.
The gradients of the displacement and velocity can be approximated as
Solving the system (70) for , , and yields
Now, we can write the coefficients in terms of nodal values of the displacements, , velocities, , and accelerations, as follows:
where , , , , , and are represent the volume, mass, damping, stiffness, capacity, and conductivity matrices, respectively; , , , , and represent the acceleration, velocity, displacement, temperature, and external force vectors, respectively, is a Green and Lindsay material constants vector, and and are coupling matrices.
Hence, the governing equations lead to the following coupled system of differential-algebraic equations (DAEs) as in Farhat et al. :
where and .
From Eq. (78) we obtain
From Eq. (83), we have
First step. Predict the displacement field: .
Third step. Correct the displacement field (81) by using the computed temperature.
4. Design sensitivity and optimization
According to Fahmy [77, 78], the design sensitivities of the nonlinear temperature field and nonlinear displacement field can be performed by the implicit differentiation of Eqs. (76) and (77), respectively, which describe the structural response with respect to the design variables, then we can compute the nonlinear thermal stresses sensitivities.
In order to solve our topology optimization problem, the method of moving asymptotes (MMA)  has been implemented as an optimizer in our topology optimization program. The benefit of MMA algorithm is that it replaces the original nonlinear, non-convex optimization problem by a sequence of approximating convex subproblems which are much easier to solve. The implemented MMA is based on the bi-directional evolutionary structural optimization (BESO), which is the evolutionary topology optimization approach that allows modification of the structure by either adding efficient material or removing inefficient material to or from the structure design [88, 89, 90, 91, 92, 93, 94, 95, 96]. This addition or removal depends upon the sensitivity analysis. Sensitivity analysis is the estimation of the response of the structure to the modification of the input design variables and is dependent upon the calculation of derivatives.
The homogenized vector of thermal expansion coefficients can be written in terms of the homogenized elastic matrix and homogenized stress-temperature coefficients vector as follows:
For the material design, the derivative of the homogenized thermal expansion coefficients vector can be expressed as
where and for any th material phase, can be calculated using the adjoint variable method  as
where, is the volume of the base cell
5. Numerical examples, results, and discussion
The proposed technique used in the current chapter should be applicable to any three-temperature nonlinear generalized magneto-thermoelastic problem. The application is for the purpose of illustration.
The two anisotropic materials considered in the calculation are monoclinic graphite-epoxy and North Sea sandstone reservoir rock, where the physical data of monoclinic graphite-epoxy material is given as follows:
Mechanical temperature coefficient:
Tensor of thermal conductivity:
Mass density and heat capacity J/(kg⋅K), Oersted, Gauss/Oersted, , and .
The physical data of the North Sea sandstone reservoir rock is given as follows:
Mechanical temperature coefficient:
Tensor of thermal conductivity:
Mass density and heat capacity J/(kg⋅K), Oersted, Gauss/Oersted, , and .
The initial and boundary conditions considered in the calculations are
In order to study the effects of anisotropy and functionally graded materials on composite microstructure, we consider the following four cases, namely, isotropic homogeneous (IH), isotropic functionally graded (IF), anisotropic homogeneous (AH), and anisotropic functionally graded (AF). Also, we considered total temperature as the considered temperature field in all calculations of this study.
Figure 2 shows the variations of the nonlinear three-temperature , , and and total temperature , with the time through composite microstructure.
Figures 3 and 4 show the variation of the nonlinear displacement sensitivities and , with time for different cases IH, IF, AH, and AF. It was shown from these figures that the anisotropy and functionally graded material have great effects on the nonlinear displacement sensitivities through the FGA composite microstructure.
Figures 5 – 7 show the variation of the nonlinear thermal stress sensitivities , , and , respectively, with time for different cases IH, IF, AH, and AF. It was noted from these figures that the anisotropy and functionally graded material have great influences on the nonlinear thermal stress sensitivities through the FGA composite microstructure.
For comparison purposes with those of other studies, we only considered one-dimensional numerical results of the considered three-temperature problem. In the considered special case, the nonlinear displacement and nonlinear thermal stress results are plotted in Figures 8 and 9 , respectively. It can be noticed from these that the BEM results, which are based on replacing one-temperature heat conduction with three-temperature heat conduction, are in excellent agreement when compared to results obtained from the finite difference method of Pazera and Jędrysiak  and the finite element method (FEM) of Xiong and Tian . We thus demonstrate the validity and accuracy of our proposed BEM technique.
Three numerical examples of BESO topological optimization of composite microstructures are performed to illustrate the optimization results of this study . In order to obtain the functionally graded parameter effects during the optimization process of the considered composite microstructure, we consider the following values , and in the one-temperature heat conduction model and the three-temperature radiative heat conduction model.
Example 1. Composite microstructures without holes or inclusions.
The mean compliance has been minimized, to obtain the maximum stiffness for the composite microstructures made from two competitive materials and without holes or inclusions. Investigation of the effect of the functionally graded parameter on the optimal composite microstructure has been shown in Table 1 for the model and in Table 2 for the model. It is noticed from these tables that the heat conduction model and functionally graded parameter have a significant effect on the topology optimization process of the multi-material FGA composite microstructures.
Example 2. Composite microstructures with circular or square holes.
The mean compliance has been minimized to obtain the maximum stiffness for the composite microstructures made from two competitive materials and with circular or square holes. Investigation of the effect of the functionally graded parameter on the optimal composite microstructure with circular holes has been shown in Table 3 for the model and in Table 4 for the model. Also, the investigation of the effect of the functionally graded parameter on the optimal composite microstructure with square holes has been shown in Table 5 for the model and in Table 6 for the model. It is noticed from these tables that the heat conduction model, functionally graded parameter, and holes shape have a significant effect on the topology optimization process of the multi-material FGA composite microstructures.
Example 3. Composite microstructures with circular or square inclusions.
The mean compliance has been minimized to obtain the maximum stiffness for the composite microstructures made from two competitive materials and with circular or square inclusions. Investigation of the effect of the functionally graded parameter on optimal composite microstructure with circular inclusions has been shown in Table 7 for the model and in Table 8 for the model. Also, the investigation of the effect of the functionally graded parameter on the optimal composite microstructure with square inclusions has been shown in Table 9 for the model and in Table 10 for the model. It is noticed from these tables that the heat conduction model, functionally graded parameter, and inclusions shape have a significant effect on the topology optimization process of the multi-material FGA composite microstructures.
The BESO topology optimization problem implemented in the numerical examples to find the distribution of the two materials in the design domain that minimize the compliance of the structure subject to a volume constraint in both phases can be stated as
where is the design variable; is the volume of the material phase, where and denote the element which is made of material; is the mean compliance; is the total load on the structure, which is the sum of mechanical and thermal loads; is the displacement vector; is the volume of the solid material; is the total number of elements; is the global stiffness matrix; is a small value (e.g., 0.0001), which guarantees that none of the elements will be removed completely from design domain; is the mechanical load vector; and is the thermal load vector. Also, the BESO parameters considered in these examples can be seen in Table 11 . The validity of our implemented BESO topology optimization technique has been demonstrated in our recent reference .
|Variable name||Variable description||Variable value|
|Final volume fraction of the material 1 for both interpolations||0.10|
|Final volume fraction of the material 2 for both interpolations||0.20|
|Evolutionary ratio for interpolation 1||2%|
|Evolutionary ratio for interpolation 2||3%|
|Volume addition ratio for interpolation 1||3%|
|Volume addition ratio for interpolation 2||2%|
|Filter ratio for interpolation 1||4 mm|
|Filter ratio for interpolation 2||3 mm|
|Convergence tolerance for both interpolations||0.01%|
|Convergence parameter for both interpolations||5|
Example 4. Laminated composite microstructure with three different sets of boundary conditions are considered in this example to validate the BEM formulation of the current study. These boundary conditions are called: simply—simply supported (SS), clamped—clamped (CC), and clamped—simply supported (CS). One-temperature (1T) and three-temperature (3T) models of nonlinear thermal stresses sensitivities results have been compared with the finite element method (FEM) results of Rajanna et al.  as well as with the finite volume method (FVM) results of Fallah and Delzendeh , which are tabulated in Table 12 for different types of boundary conditions and different methods. It can be observed that the BEM results for all the three types of boundary conditions are in excellent agreement with FEM results of  and the FVM results of .
The main aim of this chapter is to describe a new boundary element formulation for the modeling and optimization of the three-temperature nonlinear generalized magneto-thermoelastic functionally graded anisotropic (FGA) composite microstructures. The governing equations of the considered model are very difficult to solve analytically because of the nonlinearity and anisotropy. To overcome this, we propose a new boundary element formulation for solving such equations, where we used the three-temperature nonlinear radiative heat conduction equations combined with electron, ion, and phonon temperatures. Numerical results show the three-temperature distributions through composite microstructure. The effects of anisotropy and functionally graded material on the three-temperature nonlinear displacement sensitivities and nonlinear thermal stress sensitivities through the composite microstructure are very significant and pronounced. Because there are no available results in the literature to confirm the validity and accuracy of our proposed technique except for one-temperature heat conduction, we replace the three-temperature radiative heat conduction with one-temperature heat conduction as a special case from our current general study. In the considered special case, the BEM results have been compared graphically with the FDM results and FEM results, and it can be noticed that the BEM results are in excellent agreement with the FDM and FEM results. These results thus demonstrate the validity and accuracy of our proposed technique.
Numerical examples are solved using the method of moving asymptotes (MMA) algorithm based on the bi-evolutionary structural optimization method (BESO), where we used the topological optimization to manufacture three-temperature magneto-thermoelastic composite microstructures to obtain the required specific engineering properties. A new class of FGA composite microstructures consisting of two competitive materials has been studied, taking into account existence of holes or inclusions. The effects of the heat conduction model, functionally graded parameter, and holes shape and inclusions shape on the optimal composite microstructure are investigated through the considered examples with great practical interest.
The ability to understand and manipulate composite microstructures has been fundamental to our technical development over time. Today, scientists and engineers recognize the importance of composite microstructures use for economic and environmental reasons. Based on the BEM implementation and its results, this study concluded that the boundary element technique is the most suitable technique for the manufacturing of FGA composite microstructures in the future works. This technique aimed to describe the behavior of FGA composite microstructures and achieves improvement in the composition optimization and mechanical properties of the resulting FGA composite microstructures.
Due to three-temperature and numerous low-temperature and high-temperature applications in laminated composites microstructures, as a future work and based on the findings obtained in the present study, we would suggest further research to develop numerical techniques for solving the three-temperature nonlinear thermoelastic wave propagation problems and for manufacturing of advanced laminated composites. The numerical results of our considered study can provide data references for mechanical engineers, computer engineers, geotechnical engineers, geothermal engineers, technologists, new materials designers, physicists, material science researchers, and those who are interested in novel technologies in the area of three-temperature magneto-thermoelastic FGA composite microstructures. Application of three-temperature theories in advanced manufacturing technologies, with the development of soft machines and robotics in biomedical engineering and advanced manufacturing, and nonlinear generalized magneto-thermoelastic problems will be encountered more often where three-temperature radiative heat conduction will turn out to be the best choice for thermomechanical analysis in the design and analysis of advanced composite microstructures.