Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Development of a Shape Replicating Draping Unit for Continuous Layup of Unidirectional Non-Crimp Fabrics on Complex Surface Geometries
Next Article in Special Issue
Assessment of Replacement of Metal Parts by BFRP Composites into a Highly Efficient Electrical Prototype
Previous Article in Journal
Effect of Graphene Oxide as a Reinforcement in a Bio-Epoxy Composite
Previous Article in Special Issue
S-N Curve Characterisation for Composite Materials and Prediction of Remaining Fatigue Life Using Damage Function
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Objective Optimization of Functionally Graded Beams Using a Genetic Algorithm with Non-Dominated Sorting

Department of Civil Engineering, National Cheng Kung University, Tainan 70101, Taiwan
*
Author to whom correspondence should be addressed.
J. Compos. Sci. 2021, 5(4), 92; https://doi.org/10.3390/jcs5040092
Submission received: 13 March 2021 / Revised: 25 March 2021 / Accepted: 27 March 2021 / Published: 30 March 2021
(This article belongs to the Special Issue Feature Papers in Journal of Composites Science in 2021)

Abstract

:
A mixed layer-wise (LW) higher-order shear deformation theory (HSDT) is developed for the thermal buckling analysis of simply-supported, functionally graded (FG) beams subjected to a uniform temperature change. The material properties of the FG beam are assumed to be dependent on the thickness and temperature variables, and the effective material properties are estimated using either the rule of mixtures or the Mori–Tanaka scheme. The results shown in the numerical examples indicate the mixed LW HSDT solutions for critical temperature change parameters are in excellent agreement with the accurate solutions available in the literature. A multi-objective optimization of FG beams is presented to maximize the critical temperature change parameters and to minimize their total mass using a non-dominated sorting-based genetic algorithm. Some specific forms for the volume fractions of the constituents of the FG beam are assumed in advance, such as the one- and three-parameter power-law functions. The former is used in the thermal buckling analysis of the FG beams for comparison purposes, and the latter is used in their optimal design.

1. Introduction

Functionally graded (FG) beams, plates, and shells are emerging heterogeneous material structures, which are formed by mixing two- or multiple-phase materials with a pre-designed spatial distribution of volume fractions of the constituents [1,2,3,4]. Because the material properties of FG structures gradually and smoothly vary through their domain, some drawbacks can be prevented, including delamination and stress concentration, which usually occur at the interfaces between adjacent layers in the case of the laminated composite structures due to the material properties suddenly changing at these locations. On the other hand, FG structures provide a large optimal design space in engineering practice. Engineers can design the spatial distributions of the volume fractions of the constituents according to practice demands to obtain the best structural physical properties. Based on the above-mentioned advantages of FG structures, they are also becoming increasingly more popular in various high-end industries, such as aerospace, submarine, automobile, and nuclear industries, as well as for use in biomedical implants [5,6]. A variety of mechanical analyses of FG structures and single/multi-objective optimization of these structures, thus, are attracting considerable attention.
Some comprehensive literature surveys with regard to articles examining the mechanical analyses of FG structures can be found in the public literature [7,8,9,10,11,12,13,14]. Among these, the literature survey in this work focuses on articles investigating the thermal buckling behavior of FG beams using various beam theories and on articles related to the single- and multi-objective optimizations of FG beams.
Based on the Euler–Bernoulli theory (EBT), Kiani and Eslami [15] examined the thermal buckling behavior of FG beams under different boundary conditions, in which only the bending effects were considered in their analysis, rather than the shear deformation effects. The first-order shear deformation theory (FSDT) was used by Tauchert [16] to study the thermal buckling of simply supported, antisymmetric angle-ply laminates, in which the shear deformations of the laminates were simply assumed as a constant through their thickness direction. Reddy [17] developed a refined shear deformation theory (RSDT) to investigate the various mechanical behaviors of laminated composite (LC) structures. In Reddy’s formulation, the shear deformations induced in the LC structures were assumed to be a parabolic function distribution through the thickness direction; the shear correction factor commonly used in the FSDT was not needed, and the surface conditions regarding the traction stresses free were satisfied. In order to more precisely capture the shear deformation effects on the various mechanical behaviors of LC and FG structures, some advanced shear deformation theories have been proposed, including the third-order shear deformation theory (TSDT) [18,19], the sinusoidal shear deformation theory (SSDT) [20,21], the hyperbolic refined shear deformation theory (HRSDT) [22,23], and the exponential shear deformation theory (ESDT) [24]. Liu et al. [25] investigated the nonlinear response of FG structures using an iso-geometric continuum shell element method. Liu et al. [26] presented an analysis of thermo-mechanical snap-through instability of structures using a hybrid load- and displacement-controlled algorithm. Based on a shear deformation theory, Refrafi et al. [27] examined the hygro-thermo-mechanical buckling behavior of simply-supported FG sandwich plates resting on a Winkler–Pasternak foundation. Based on Carrera’s unified formulation (CUF) [28], the above-mentioned advanced and refined shear deformation theories can be hierarchically derived and included as special cases [29]. Some comparative studies for different mechanical behavior analyses of FG beams using assorted advanced and refined shear deformation beam theories have been carried out, including thermal buckling and post-buckling analyses [30], bending and free vibration analyses [31], thermal buckling analysis [32], and free vibration analyses [33,34]. The above-mentioned theories can be summarized as equivalent-single-layered (ESL) theories, in which a global coordinate system was introduced to describe the displacement field induced in LC and FG structures.
Because ESL theories are often inadequate for determining the three-dimensional (3D) stress field induced in loaded LC and FG beams, some layer-wise (LW) theories for these beams have been presented. Based on Hamilton’s theory, Tahani [35] developed a displacement-based LW theory for the static bending and free vibration analyses of LC beams. Lee and Saravanos [36] developed an LW FSDT for coupled thermo-piezo-electric composite beams under thermal loads. Shimpi and Ainapure [37] and Shimpi and Ghugal [38] developed an LW trigonometric shear deformation theory for the static analysis of LC beams. Based on an LW higher-order shear deformation theory (HSDT), Pandey and Pradyumna [39] developed a finite element method (FEM) for the static and dynamic analyses of FG sandwich plates. Bayat and Toussi [40] carried out thermal buckling and thermal postbuckling analyses of LC beams reinforced with shape memory alloy wires using a displacement-based LW theory. Based on the stationary principle of minimum potential energy [41,42] combined with the Lagrange multiplier method, Wu and Kuo [43] and Wu and Chen [44] developed a mixed LW HSDT for the bending, free vibration, and buckling behavior of LC plates. Wu and Xu [45] extended this approach to develop the strong and weak formulations of a mixed higher-order shear deformation theory for the static analysis of FG beams under thermo-mechanical loads. A comprehensive overview with regard to the development, numerical implementation, and application of LW theories was undertaken by Liew et al. [46].
Some articles related to the optimal design of LC structures and FG beams have been presented using different single-and multi-objective optimization algorithms combined with advanced and refined shear deformation beam theories. Walker et al. [47,48] developed an optimal procedure to select the best material combinations for minimum weight and cost of hybrid composite plates and for minimum mass of sandwiched composite cylindrical shells, respectively. Houmat [49] and Guenanou and Houmat [50] proposed a layer-wise optimization method to obtain an optimal lay-up design for maximum fundamental frequency of variable stiffness LC plates and symmetrically LC circular plates, respectively. Cho and Ha [51] presented volume fraction optimization for the purpose of minimizing thermal stresses induced in Ni-Al2O3 FG beams. Based on the Vlasov thin-walled theory [52], Nguyen and Lee [53] committed to the optimal design of thin-walled FG beams for buckling problems, in which an LW cubic interpolation for the through-thickness distribution of the volume fraction of the constituents was assumed, and a genetic algorithm (GA) was used as the optimal tool [54,55]. Kamarian et al. [56] studied the optimization of material compositions of FG beams to maximize their fundamental natural frequencies using the firefly algorithm and an adaptive neuro-fuzzy inference system. Yas et al. [57,58] engaged in the optimization of material compositions of a three-parameter FG beam in order to maximize its fundamental frequencies using the imperialist competitive algorithm, an artificial neural network [59], and a generalized differential quadrature (DQ) method [60,61,62]. Goupee and Vel [63] proposed a methodology to maximize the fundamental frequencies of FG structures by tailoring their material distributions. Na and Kim [64,65,66] investigated the volume fraction optimization of FG plates and panels under thermal loads to minimize the induced thermal stresses and to maximize the critical temperature parameters, in which a bi-objective optimization was considered and the optimal tool was a quasi-Newton method. Walker and Smith [67] presented a multi-objective optimization technique for minimizing the weighted sum of the mass and deflection of LC structures, in which a GA and an FEM were used. Tornabene and Ceruti [68] dealt with mixed static and dynamic optimization of four-parameter FG structures. The FSDT combined with the generalized DQ method was used for the static and dynamic analyses of the FG structures, and three different optimization schemes were used as the optimal tools, including the particle swarm optimization, the Monte Carlo, and the GA approaches.
In this work, the authors aim at investigating the material composition optimization of an FG beam under a uniform temperature change in order to maximize the critical temperature change parameter and minimize the total mass of the FG beam. A mixed LW HSDT [45] is developed for the thermal buckling analysis of the FG beam. The current solutions for the critical temperature change parameters of the FG beam obtained through considering the temperature-dependent (TD) and the temperature-independent (TI) material properties are examined, and solutions obtained using the von Kármán geometric nonlinear strains (GNS) and the full GNS are also estimated. A non-dominated sorting-based GA [69,70] is used for the current multi-objective optimization analysis in the current issue. The through-thickness distribution of the material properties of the FG beam is assumed to be a three-parameter power-law function of the volume fractions of the constituents [71,72], the material-property gradient indices of which are, thus, to be determined for the optimal material profile.

2. Effective Material Properties

In this work, two micromechanics models, the rule of mixtures [73,74] and the Mori–Tanaka scheme [75], are used to estimate the effective material properties of the FG beam, as described in the following sections.

2.1. The Rule of Mixtures

According to the rule of mixtures, the through-thickness distributions of the effective material properties of the FG beam are written in the following form,
F e f f z = Γ c z F c + Γ m z F m = F m + F c F m Γ c z ,
where Γ c and Γ m represent the volume fractions of the ceramic and metal materials of the constituents of the FG beam, respectively, such that Γ c + Γ m = 1 . F can be one of the engineering constants, including Young’s modulus E , Poisson’s ratio υ , thermal expansion coefficient α , and mass density ρ . The subscripts m and c represent the metal material and the ceramic material, respectively.

2.2. The Mori-Tanaka Scheme

Within the frame work of Eshelby’s elastic inclusion theory [73,74] restricted to a single inclusion in a semi-infinite elastic, homogeneous and isotropic medium, Mori and Tanaka [75] extended the original theory to more general cases of multiple inclusions embedded into a finite elastic medium, in which an equivalent inclusion-average stress method was used.
For a two-phase isotropic material, the effective material properties are written as the explicit formula [76,77] as follows:
K z = Γ c K c K m / [ 1 + 1 Γ c K c K m / K m + 4 / 3 G m ] + K m ,
G ( z ) = Γ c G c G m / [ 1 + 1 Γ c G c G m / G m + f m ] + G m ,
α ( z ) = Γ c α c α m / { K z [ 1 K z K c / K z K m ] } + α m ,
in which f m = G m 9 K m + 8 G m / 6 K m + 2 G m , and K and G represent the volume and the shear moduli, respectively.
In later sections of this paper, if no specific mention is given the effective material properties of the FG beam are estimated using the rule of mixtures.

3. The Mixed LW HSDT for FG Beams

In this work, the authors develop a unified formulation of the mixed LW HSDT for the thermal buckling analysis of simply supported FG beams subjected to a uniform temperature change, in which the material properties of the FG beam are considered to be temperature-dependent. The configuration and coordinates of the FG beam are shown in Figure 1a, in which h and L represent the thickness and the length of the FG beam, respectively. In the analysis, the FG beam is artificially divided into Nl layers as shown in Figure 1b, and the thickness of each individual layer constituting the beam is hm (m = 1 − Nl), such that m = 1 N l h m = h .

3.1. Strong-Form Formulation

In the current LW HSDT, the displacement field for a typical individual layer is given as follows:
u ( m ) x , z m = u 0 ( m ) x + z m u 1 ( m ) x + z m 2 u 2 ( m ) x + z m 3 u 3 ( m ) x ,
w ( m ) x , z m = w 0 ( m ) x + z m w 1 ( m ) x + z m 2 w 2 ( m ) x + z m 3 w 3 ( m ) x ,
where m = 1, 2, …, Nl; u 0 ( m ) and w 0 ( m ) denote the mid-plane displacements of the layer in the x and z directions, and their ith-order expansion terms along the local thickness direction are u i ( m ) and w i ( m ) ( i = 1 , 2 , and 3 ) . The displacement components in the y direction are taken to be zero.
According to the perfect bonding assumptions at the interfaces between adjacent layers, the corresponding displacement continuity conditions at these places are given as
f x ( k ) = u 0 ( k + 1 ) h k + 1 / 2 u 1 ( k + 1 ) + h k + 1 2 / 4 u 2 ( k + 1 ) h k + 1 3 / 8 u 3 ( k + 1 ) u 0 ( k ) + h k / 2 u 1 ( k ) + h k 2 / 4 u 2 ( k ) + h k 3 / 8 u 3 ( k ) = 0 ,
f z ( k ) = w 0 ( k + 1 ) h k + 1 / 2 w 1 ( k + 1 ) + h k + 1 2 / 4 w 2 ( k + 1 ) h k + 1 3 / 8 w 3 ( k + 1 ) w 0 ( k ) + h k / 2 w 1 ( k ) + h k 2 / 4 w 2 ( k ) + h k 3 / 8 w 3 ( k ) = 0 ,
where k = 1, 2,…, (Nl − 1).
The strain-displacement relationship is given as
ε x ( m ) = ε x ( m ) l + ε x ( m ) n l = u 0 ( m ) , x + z m u 1 ( m ) , x + z m 2 u 2 ( m ) , x + z m 3 u 3 ( m ) , x + 1 / 2 u 0 ( m ) , x + z m u 1 ( m ) , x + z m 2 u 2 ( m ) , x + z m 3 u 3 ( m ) , x 2 + w 0 ( m ) , x + z m w 1 ( m ) , x + z m 2 w 2 ( m ) , x + z m 3 w 3 ( m ) , x 2 ,
ε z ( m ) = w ( m ) , z m = w 1 ( m ) + 2 z m w 2 ( m ) + 3 z m 2 w 3 ( m ) ,
γ x z ( m ) = u ( m ) , z m + w ( m ) , x = w 0 ( m ) , x + u 1 ( m ) + z m w 1 ( m ) , x + 2 u 2 ( m ) + z m 2 w 2 ( m ) , x + 3 u 3 ( m ) + z m 3 w 3 ( m ) , x ,
where the commas denote the derivative of the suffix variable, and remaining strains are zeroes, including ε y ( m ) , γ y z ( m ) , and γ x y ( m ) . ε x ( m ) n l denotes the full GNS, and when the von Kármán GNS is considered, it is reduced as ε x ( m ) n l = 1 / 2 w 0 ( m ) , x + z m w 1 ( m ) , x + z m 2 w 2 ( m ) , x + z m 3 w 3 ( m ) , x 2 . It is noted that the displacement continuity conditions at the interfaces between adjacent layers are imposed in this formulation, while the strain continuity conditions at these places are discontinuous.
The stress-strain relationship for an orthotropic material in thermal environment is given as
σ x ( m ) σ z ( m ) τ x z ( m ) = Q 11 ( m ) Q 13 ( m ) 0 Q 13 ( m ) Q 33 ( m ) 0 0 0 Q 55 ( m ) ε x ( m ) ε z ( m ) γ x z ( m ) Q 1 α ( m ) Q 3 α ( m ) 0 Δ T ,
where the symbol Δ T denotes the temperature change, measured from a room temperature of 300K. Q 11 ( m ) = 1 υ 23 υ 32 / E 2 E 3 Δ ( m ) , Q 13 ( m ) = υ 13 + υ 12 υ 32 / E 1 E 2 Δ ( m ) , Q 33 ( m ) = 1 υ 12 υ 21 / E 1 E 2 Δ ( m ) , Q 55 ( m ) = G 13 ( m ) , Δ = 1 υ 12 υ 21 υ 23 υ 32 υ 31 υ 13 2 υ 21 υ 32 υ 13 / E 1 E 2 E 3 ( m ) , Q 1 α ( m ) = Q 11 ( m ) α 1 ( m ) + Q 13 ( m ) α 3 ( m ) , Q 3 α ( m ) = Q 13 ( m ) α 1 ( m ) + Q 33 ( m ) α 3 ( m ) , in which the subscripts 1, 2, and 3 denote the principle axes of the material properties, and that E , υ , a n d α represent the Young’s modulus, Poisson’s ratio, and the thermal expansion coefficients. For isotropic materials, these stiffness coefficients will be reduced to Q 11 ( m ) = Q 33 ( m ) = E 1 υ / 1 + υ 1 2 υ ( m ) , Q 13 ( m ) = E υ / 1 + υ 1 2 υ ( m ) , Q 55 ( m ) = E / 2 1 + υ ( m ) , and α 1 ( m ) = α 3 ( m ) = α . Note that these engineering constants E , υ , and α in the analysis are considered to be dependent from the thickness and the temperature variables, and Δ T is a constant. In addition, other stress components are zero, including σ y ( m ) , τ y z ( m ) , and τ x y ( m ) , because a two-dimensional beam model is used.
The governing equations and associated boundary conditions are derived using the stationary principle of minimum potential energy combined with the Lagrange multiplier method, in which the displacement continuity conditions at the interfaces between adjacent layers given in Equations (7) and (8) are multiplied by the Lagrange multipliers and then substituted into the potential energy functional as the constraints, such that the extended potential energy functional of the Nl-layered FG beam is given, as follows:
Π p = m = 1 N l 0 L h m / 2 h m / 2 1 / 2 σ x ( m ) ε x ( m ) + 1 / 2 σ z ( m ) ε z ( m ) + 1 / 2 τ x z ( m ) γ x z ( m ) d z m d x Δ T m = 1 N l 0 L h m / 2 h m / 2 Q 1 α ( m ) ε x ( m ) n l d z m d x + k = 1 N l 1 0 L λ x ( k ) f x ( k ) + λ z ( k ) f z ( k ) d x ,
where λ x ( m ) and λ z ( m ) are the Lagrange multipliers, which are identical to the transverse shear and normal stress components (i.e., τ x z ( m ) and σ z ( m ) ) induced at the interfaces between adjacent layers, respectively.
Applying the stationary principle of minimum potential energy, following a standard variational process, and integrating the stress variables through the thickness direction, the authors obtain
δ Π p = m = 1 N l 0 L N x ( m ) N x t ( m ) δ u 0 ( m ) , x + M x ( m ) M x t ( m ) δ u 1 ( m ) , x + P x ( m ) P x t ( m ) δ u 2 ( m ) , x + R x ( m ) R x t ( m ) δ u 3 ( m ) , x + N z ( m ) N z t ( m ) δ w 1 ( m ) + 2 M z ( m ) M z t ( m ) δ w 2 ( m ) + 3 P z ( m ) P z t ( m ) δ w 3 ( m ) + N x z ( m ) δ w 0 ( m ) , x + δ u 1 ( m ) + M x z ( m ) δ w 1 ( m ) , x + 2 δ u 2 ( m ) + P x z ( m ) δ w 2 ( m ) , x + 3 δ u 3 ( m ) + R x z ( m ) δ w 3 ( m ) , x d z m d x + Δ T m = 1 N l 0 L S 1 ( m ) δ u 0 ( m ) + S 2 ( m ) δ u 1 ( m ) + S 3 ( m ) δ u 2 ( m ) + S 4 ( m ) δ u 3 ( m ) + S 5 ( m ) δ w 0 ( m ) + S 6 ( m ) δ w 1 ( m ) + S 7 ( m ) δ w 2 ( m ) + S 8 ( m ) δ w 3 ( m ) d x + k = 1 N l 1 0 L δ λ x ( k ) f x ( k ) + δ λ z ( k ) f z ( k ) + λ x ( k ) δ f x ( k ) + λ z ( k ) δ f z ( k ) d x + boundary terms = 0 ,
where
N x ( m ) N z ( m ) N x z ( m ) M x ( m ) M z ( m ) M x z ( m ) P x ( m ) P z ( m ) P x z ( m ) = h m / 2 h m / 2 1 z m z m 2 σ x ( m ) σ z ( m ) τ x z ( m ) d z m , R x ( m ) R x z ( m ) = h m / 2 h m / 2 z m 3 σ x ( m ) τ x z ( m ) d z m ,
N x t ( m ) N z t ( m ) M x t ( m ) M z t ( m ) P x t ( m ) P z t ( m ) = h m / 2 h m / 2 1 z m z m 2 Q 1 α ( m ) Δ T Q 3 α ( m ) Δ T d z m , R x t ( m ) = h m / 2 h m / 2 z m 3 Q 1 α ( m ) Δ T d z m ,
S 1 ( m ) = A 1 α ( m ) u 0 ( m ) , x x + B 1 α ( m ) u 1 ( m ) , x x + D 1 α ( m ) u 2 ( m ) , x x + F 1 α ( m ) u 3 ( m ) , x x ,
S 2 ( m ) = B 1 α ( m ) u 0 ( m ) , x x + D 1 α ( m ) u 1 ( m ) , x x + F 1 α ( m ) u 2 ( m ) , x x + H 1 α ( m ) u 3 ( m ) , x x ,
S 3 ( m ) = D 1 α ( m ) u 0 ( m ) , x x + F 1 α ( m ) u 1 ( m ) , x x + H 1 α ( m ) u 2 ( m ) , x x + J 1 α ( m ) u 3 ( m ) , x x ,
S 4 ( m ) = F 1 α ( m ) u 0 ( m ) , x x + H 1 α ( m ) u 1 ( m ) , x x + J 1 α ( m ) u 2 ( m ) , x x + L 1 α ( m ) u 3 ( m ) , x x ,
S 5 ( m ) = A 1 α ( m ) w 0 ( m ) , x x + B 1 α ( m ) w 1 ( m ) , x x + D 1 α ( m ) w 2 ( m ) , x x + F 1 α ( m ) w 3 ( m ) , x x ,
S 6 ( m ) = B 1 α ( m ) w 0 ( m ) , x x + D 1 α ( m ) w 1 ( m ) , x x + F 1 α ( m ) w 2 ( m ) , x x + H 1 α ( m ) w 3 ( m ) , x x ,
S 7 ( m ) = D 1 α ( m ) w 0 ( m ) , x x + F 1 α ( m ) w 1 ( m ) , x x + H 1 α ( m ) w 2 ( m ) , x x + J 1 α ( m ) w 3 ( m ) , x x ,
S 8 ( m ) = F 1 α ( m ) w 0 ( m ) , x x + H 1 α ( m ) w 1 ( m ) , x x + J 1 α ( m ) w 2 ( m ) , x x + L 1 α ( m ) w 3 ( m ) , x x ,
A 1 α ( m ) B 1 α ( m ) D 1 α ( m ) F 1 α ( m ) H 1 α ( m ) J 1 α ( m ) L 1 α ( m ) = h m / 2 h m / 2 Q 1 α ( m ) 1 z m z m 2 z m 3 z m 4 z m 5 z m 6 d z m .
Performing Equation (14) the integration by part yields
δ Π p = m = 1 N l 0 L N x ( m ) N x t ( m ) , x δ u 0 ( m ) M x ( m ) M x t ( m ) , x δ u 1 ( m ) P x ( m ) P x t ( m ) , x δ u 2 ( m ) R x ( m ) R x t ( m ) , x δ u 3 ( m ) + N z ( m ) N z t ( m ) δ w 1 ( m ) + 2 M z ( m ) M z t ( m ) δ w 2 ( m ) + 3 P z ( m ) P z t ( m ) δ w 3 ( m ) N x z ( m ) , x δ w 0 ( m ) + N x z ( m ) δ u 1 ( m ) M x z ( m ) , x δ w 1 ( m ) + 2 M x z ( m ) δ u 2 ( m ) P x z ( m ) , x δ w 2 ( m ) + 3 P x z ( m ) δ u 3 ( m ) R x z ( m ) , x δ w 3 ( m ) d x + Δ T m = 1 N l 0 L S 1 ( m ) δ u 0 ( m ) + S 2 ( m ) δ u 1 ( m ) + S 3 ( m ) δ u 2 ( m ) + S 4 ( m ) δ u 3 ( m ) + S 5 ( m ) δ w 0 ( m ) + S 6 ( m ) δ w 1 ( m ) + S 7 ( m ) δ w 2 ( m ) + S 8 ( m ) δ w 3 ( m ) d x + k = 1 N l 1 0 L δ λ x ( k ) f x ( k ) + δ λ z ( k ) f z ( k ) + λ x ( k ) δ f x ( k ) + λ z ( k ) δ f z ( k ) d x + ( boundary terms ) = 0 .
According to Equation (15), the Euler–Lagrange equations of the mixed LW HSDT can be obtained as follows:
δ u 0 ( m ) : N x ( m ) , x + λ x m 1 λ x ( m ) = Δ T S 1 ( m ) ,
δ u 1 ( m ) : M x ( m ) , x + N x z ( m ) + h m / 2 λ x m 1 + λ x ( m ) = Δ T S 2 ( m ) ,
δ u 2 ( m ) : P x ( m ) , x + 2 M x z ( m ) + h m 2 / 4 λ x m 1 λ x ( m ) = Δ T S 3 ( m ) ,
δ u 3 ( m ) : R x ( m ) , x + 3 P x z ( m ) + h m 3 / 8 λ x m 1 + λ x ( m ) = Δ T S 4 ( m ) ,
δ w 0 ( m ) : N x z ( m ) , x + λ z m 1 λ z ( m ) = Δ T S 5 ( m ) ,
δ w 1 ( m ) : M x z ( m ) , x + N z ( m ) + h m / 2 λ z m 1 + λ z ( m ) = Δ T S 6 ( m ) ,
δ w 2 ( m ) : P x z ( m ) , x + 2 M z ( m ) + h m 2 / 4 λ z m 1 λ z ( m ) = Δ T S 7 ( m ) ,
δ w 3 ( m ) : R x z ( m ) , x + 3 P z ( m ) + h m 3 / 8 λ z m 1 + λ z ( m ) = Δ T S 8 ( m ) ,
δ λ x ( k ) : f x ( k ) = 0 ,
δ λ z ( k ) : f z ( k ) = 0 ,
where m = 1, 2,…, Nl and k = 1, 2,…, (Nl − 1).
The possible boundary conditions are given as
Either   N x ( m ) = N x t ( m )   or   u 0 ( m ) = u ¯ 0 ( m )
Either   M x ( m ) = M x t ( m )   or   u 1 ( m ) = u ¯ 1 ( m ) ,
Either   P x ( m ) = P x t ( m )   or   u 2 ( m ) = u ¯ 2 ( m ) ,
Either   R x ( m ) = R x t ( m )   or   u 3 ( m ) = u ¯ 3 ( m ) ,
Either   N x z ( m ) = N x z t ( m )   or   w 0 ( m ) = w ¯ 0 ( m ) ,
Either   M x z ( m ) = M x z t ( m )   or   w 1 ( m ) = w ¯ 1 ( m ) ,
Either   P x z ( m ) = P x z t ( m )   or   w 2 ( m ) = w ¯ 2 ( m ) ,
Either   R x z ( m ) = R x z t ( m )   or   w 3 ( m ) = w ¯ 3 ( m ) ,
where u ¯ i ( m ) and w ¯ i ( m ) ( i = 0 , 1 , 2 , and 3 ) are the prescribed displacement components on the edges. The definition of each force resultant component and its relationship with the displacement components are given in Appendix A.
Substituting Equations (A1)–(A11) into Equations (16)–(25), the authors obtain the Euler–Lagrange equations of the current mixed LW HSDBT in terms of all displacement components, and they are given as follows:
δ u 0 ( m ) : A 11 ( m ) u 0 ( m ) , x x B 11 ( m ) u 1 ( m ) , x x D 11 ( m ) u 2 ( m ) , x x F 11 ( m ) u 3 ( m ) , x x A 13 ( m ) w 1 ( m ) , x 2 B 13 ( m ) w 2 ( m ) , x     3 D 13 ( m ) w 3 ( m ) , x + λ x m 1 λ x ( m ) = Δ T A 1 α ( m ) u 0 ( m ) , x x + B 1 α ( m ) u 1 ( m ) , x x + D 1 α ( m ) u 2 ( m ) , x x + F 1 α ( m ) u 3 ( m ) , x x ,
δ u 1 ( m ) : B 11 ( m ) u 0 ( m ) , x x D 11 ( m ) u 1 ( m ) , x x F 11 ( m ) u 2 ( m ) , x x H 11 ( m ) u 3 ( m ) , x x B 13 ( m ) w 1 ( m ) , x 2 D 13 ( m ) w 2 ( m ) , x 3 F 13 ( m ) w 3 ( m ) , x + A 55 ( m ) w 0 ( m ) , x + u 1 ( m ) + B 55 ( m ) w 1 ( m ) , x + 2 u 2 ( m ) + D 55 ( m ) w 2 ( m ) , x + 3 u 3 ( m ) + F 55 ( m ) w 3 ( m ) , x + h m / 2 λ x m 1 + λ x ( m ) = Δ T B 1 α ( m ) u 0 ( m ) , x x + D 1 α ( m ) u 1 ( m ) , x x + F 1 α ( m ) u 2 ( m ) , x x + H 1 α ( m ) u 3 ( m ) , x x ,
δ u 2 ( m ) : D 11 ( m ) u 0 ( m ) , x x F 11 ( m ) u 1 ( m ) , x x H 11 ( m ) u 2 ( m ) , x x J 11 ( m ) u 3 ( m ) , x x D 13 ( m ) w 1 ( m ) , x 2 F 13 ( m ) w 2 ( m ) , x 3 H 13 ( m ) w 3 ( m ) , x + 2 B 55 ( m ) w 0 ( m ) , x + u 1 ( m ) + 2 D 55 ( m ) w 1 ( m ) , x + 2 u 2 ( m ) + 2 F 55 ( m ) w 2 ( m ) , x + 3 u 3 ( m )           + 2 H 55 ( m ) w 3 ( m ) , x + h m 2 / 4 λ x m 1 + λ x ( m ) = Δ T D 1 α ( m ) u 0 ( m ) , x x + F 1 α ( m ) u 1 ( m ) , x x + H 1 α ( m ) u 2 ( m ) , x x + J 1 α ( m ) u 3 ( m ) , x x ,
δ u 3 ( m ) : F 11 ( m ) u 0 ( m ) , x x H 11 ( m ) u 1 ( m ) , x x J 11 ( m ) u 2 ( m ) , x x L 11 ( m ) u 3 ( m ) , x x F 13 ( m ) w 1 ( m ) , x 2 H 13 ( m ) w 2 ( m ) , x 3 J 13 ( m ) w 3 ( m ) , x + 3 D 55 ( m ) w 0 ( m ) , x + u 1 ( m ) + 3 F 55 ( m ) w 1 ( m ) , x + 2 u 2 ( m ) + 3 H 55 ( m ) w 2 ( m ) , x + 3 u 3 ( m ) + 3 J 55 ( m ) w 3 ( m ) , x + h m 3 / 8 λ x m 1 + λ x ( m ) = Δ T F 1 α ( m ) u 0 ( m ) , x x + H 1 α ( m ) u 1 ( m ) , x x + J 1 α ( m ) u 2 ( m ) , x x + L 1 α ( m ) u 3 ( m ) , x x ,
δ w 0 ( m ) : A 55 ( m ) w 0 ( m ) , x x + u 1 ( m ) , x B 55 ( m ) w 1 ( m ) , x x + 2 u 2 ( m ) , x D 55 ( m ) w 2 ( m ) , x x + 3 u 3 ( m ) , x F 55 ( m ) w 3 ( m ) , x x + λ z m 1 λ z ( m ) = Δ T A 1 α ( m ) w 0 ( m ) , x x + B 1 α ( m ) w 1 ( m ) , x x + D 1 α ( m ) w 2 ( m ) , x x + F 1 α ( m ) w 3 ( m ) , x x ,
δ w 1 ( m ) : B 55 ( m ) w 0 ( m ) , x x + u 1 ( m ) , x D 55 ( m ) w 1 ( m ) , x x + 2 u 2 ( m ) , x F 55 ( m ) w 2 ( m ) , x x + 3 u 3 ( m ) , x H 55 ( m ) w 3 ( m ) , x x + A 13 ( m ) u 0 ( m ) , x + B 13 ( m ) u 1 ( m ) , x + D 13 ( m ) u 2 ( m ) , x + F 13 ( m ) u 3 ( m ) , x + A 33 ( m ) w 1 ( m ) + 2 B 33 ( m ) w 2 ( m ) + 3 D 33 ( m ) w 3 ( m ) + h m / 2 λ z m 1 + λ z ( m ) = Δ T B 1 α ( m ) w 0 ( m ) , x x + D 1 α ( m ) w 1 ( m ) , x x + F 1 α ( m ) w 2 ( m ) , x x + H 1 α ( m ) w 3 ( m ) , x x ,
δ w 2 ( m ) : D 55 ( m ) w 0 ( m ) , x x + u 1 ( m ) , x F 55 ( m ) w 1 ( m ) , x x + 2 u 2 ( m ) , x H 55 ( m ) w 2 ( m ) , x x + 3 u 3 ( m ) , x J 55 ( m ) w 3 ( m ) , x x + 2 B 13 ( m ) u 0 ( m ) , x + 2 D 13 ( m ) u 1 ( m ) , x + 2 F 13 ( m ) u 2 ( m ) , x + 2 H 13 ( m ) u 3 ( m ) , x + 2 B 33 ( m ) w 1 ( m ) + 4 D 33 ( m ) w 2 ( m ) + 6 F 33 ( m ) w 3 ( m ) + h m 2 / 4 λ z m 1 λ z ( m ) = Δ T D 1 α ( m ) w 0 ( m ) , x x + F 1 α ( m ) w 1 ( m ) , x x + H 1 α ( m ) w 2 ( m ) , x x + J 1 α ( m ) w 3 ( m ) , x x ,
δ w 3 ( m ) : F 55 ( m ) w 0 ( m ) , x x + u 1 ( m ) , x H 55 ( m ) w 1 ( m ) , x x + 2 u 2 ( m ) , x J 55 ( m ) w 2 ( m ) , x x + 3 u 3 ( m ) , x L 55 ( m ) w 3 ( m ) , x x + 3 D 13 ( m ) u 0 ( m ) , x + 3 F 13 ( m ) u 1 ( m ) , x + 3 H 13 ( m ) u 2 ( m ) , x + 3 J 13 ( m ) u 3 ( m ) , x + 3 D 33 ( m ) w 1 ( m ) + 6 F 33 ( m ) w 2 ( m ) + 9 H 33 ( m ) w 3 ( m ) + h m 3 / 8 λ z m 1 + λ z ( m ) = Δ T F 1 α ( m ) w 0 ( m ) , x x + H 1 α ( m ) w 1 ( m ) , x x + J 1 α ( m ) w 2 ( m ) , x x + L 1 α ( m ) w 3 ( m ) , x x ,
δ λ x ( k ) : f x ( k ) = 0 ,
δ λ z ( k ) : f z ( k ) = 0 ,
where m = 1, 2,…, Nl and k = 1, 2,…, (Nl − 1).
The total number of Euler–Lagrange equations, i.e., Equations (34)–(43), is 10 N l 2 , 8 N l 2 , and 6 N l 2 for the mixed layer-wise third-order (LW3), second-order (LW2), and first-order (LW1) shear deformation theories, respectively, are in terms of the same number of unknowns as those of the corresponding Euler–Lagrange equations. The strong-form formulations of the mixed LW HSDT are, thus, obtained, including the Euler–Lagrange equations (Equations (34)–(43)) and the possible boundary conditions (Equations (26)–(33)).

3.2. Applications

Equations (34)–(43) associated with a set of boundary conditions (i.e., Equations (26)–(33)) can be composed as a well-post boundary value problem, and its Navier-type analytical solutions for the thermal buckling behavior of simply supported, multi-layered FG beams can be obtained using the Fourier series expansion method.
By satisfying the simply supported boundary conditions, various field variables of the mth-layer are, thus, expanded as
u i ( m ) = m ^ = 1 u i m ^ ( m ) cos m ^ π x / L ,
w i ( m ) = m ^ = 1 w i m ^ ( m ) sin m ^ π x / L ,
λ x ( k ) = m ^ = 1 λ x m ^ ( k ) cos m ^ π x / L ,
λ z ( k ) = m ^ = 1 λ z m ^ ( k ) sin m ^ π x / L ,
where the subscript i = 0, 1, 2, and 3; m = 1 , 2 , , N l ; k = 1 , 2 , , N l 1 .
Substituting Equations (44)–(47) in the Euler–Lagrange Equations (i.e., Equations (34)–(43)), the authors obtain
δ u 0 ( m ) : m ¯ 2 A 11 ( m ) u 0 m ^ ( m ) + m ¯ 2 B 11 ( m ) u 1 m ^ ( m ) + m ¯ 2 D 11 ( m ) u 2 m ^ ( m ) + m ¯ 2 F 11 ( m ) u 3 m ^ ( m ) m ¯ A 13 ( m ) w 1 m ^ ( m ) 2 m ¯ B 13 ( m ) w 2 m ^ ( m )   3 m ¯ D 13 ( m ) w 3 m ^ ( m ) + λ x m ^ m 1 λ x m ^ ( m ) = Δ T m ¯ 2 A 1 α ( m ) u 0 m ^ ( m ) + m ¯ 2 B 1 α ( m ) u 1 m ^ ( m ) + m ¯ 2 D 1 α ( m ) u 2 m ^ ( m ) + m ¯ 2 F 1 α ( m ) u 3 m ^ ( m ) ,
δ u 1 ( m ) : m ¯ 2 B 11 ( m ) u 0 m ^ ( m ) + m ¯ 2 D 11 ( m ) + A 55 ( m ) u 1 m ^ ( m ) + m ¯ 2 F 11 ( m ) + 2 B 55 ( m ) u 2 m ^ ( m ) + m ¯ 2 H 11 ( m ) + 3 D 55 ( m ) u 3 m ^ ( m )       + m ¯ A 55 ( m ) w 0 m ^ ( m ) + m ¯ B 13 ( m ) + m ¯ B 55 ( m ) w 1 m ^ ( m ) + 2 m ¯ D 13 ( m ) + m ¯ D 55 ( m ) w 2 m ^ ( m ) + 3 m ¯ F 13 ( m ) + m ¯ F 55 ( m ) w 3 m ^ ( m ) + h m / 2 λ x m ^ m 1 + λ x m ^ ( m ) = Δ T m ¯ 2 B 1 α ( m ) u 0 m ^ ( m ) + m ¯ 2 D 1 α ( m ) u 1 m ^ ( m ) + m ¯ 2 F 1 α ( m ) u 2 m ^ ( m ) + m ¯ 2 H 1 α ( m ) u 3 m ^ ( m ) ,
δ u 2 ( m ) : m ¯ 2 D 11 ( m ) u 0 m ^ ( m ) + m ¯ 2 F 11 ( m ) + 2 B 55 ( m ) u 1 m ^ ( m ) + m ¯ 2 H 11 ( m ) + 4 D 55 ( m ) u 2 m ^ ( m ) + m ¯ 2 J 11 ( m ) + 6 F 55 ( m ) u 3 m ^ ( m )         + 2 m ¯ B 55 ( m ) w 0 m ^ ( m ) + m ¯ D 13 ( m ) + 2 m ¯ D 55 ( m ) w 1 m ^ ( m ) + 2 m ¯ F 13 ( m ) + 2 m ¯ F 55 ( m ) w 2 m ^ ( m ) + 3 m ¯ H 13 ( m ) + 2 m ¯ H 55 ( m ) w 3 m ^ ( m )     + h m 2 / 4 λ x m ^ m 1 λ x m ^ ( m ) = Δ T m ¯ 2 D 1 α ( m ) u 0 m ^ ( m ) + m ¯ 2 F 1 α ( m ) u 1 m ^ ( m ) + m ¯ 2 H 1 α ( m ) u 2 m ^ ( m ) + m ¯ 2 J 1 α ( m ) u 3 m ^ ( m ) ,
δ u 3 ( m ) : m ¯ 2 F 11 ( m ) u 0 m ^ ( m ) + m ¯ 2 H 11 ( m ) + 3 D 55 ( m ) u 1 m ^ ( m ) + m ¯ 2 J 11 ( m ) + 6 F 55 ( m ) u 2 m ^ ( m ) + m ¯ 2 L 11 ( m ) + 9 H 55 ( m ) u 3 m ^ ( m )   + 3 m ¯ D 55 ( m ) w 0 m ^ ( m ) + m ¯ F 13 ( m ) + 3 m ¯ F 55 ( m ) w 1 m ^ ( m ) + 2 m ¯ H 13 ( m ) + 3 m ¯ H 55 ( m ) w 2 m ^ ( m ) + 3 m ¯ J 13 ( m ) + 3 m ¯ J 55 ( m ) w 3 m ^ ( m )   + h m 3 / 8 λ x m ^ m 1 + λ x m ^ ( m ) = Δ T m ¯ 2 F 1 α ( m ) u 0 m ^ ( m ) + H 2 F 1 α ( m ) u 1 m ^ ( m ) + m ¯ 2 J 1 α ( m ) u 2 m ^ ( m ) + m ¯ 2 L 1 α ( m ) u 3 m ^ ( m ) ,
δ w 0 ( m ) : m ¯ A 55 ( m ) u 1 m ^ ( m ) + 2 m ¯ B 55 ( m ) u 2 m ^ ( m ) + 3 m ¯ D 55 ( m ) u 3 m ^ ( m ) + m ¯ 2 A 55 ( m ) w 0 m ^ ( m ) + m ¯ 2 B 55 ( m ) w 1 m ^ ( m ) + m ¯ 2 D 55 ( m ) w 2 m ^ ( m )                 + m ¯ 2 F 55 ( m ) w 3 m ^ ( m ) + λ z m ^ m 1 λ z m ^ ( m ) = Δ T m ¯ 2 A 1 α ( m ) w 0 m ^ ( m ) + m ¯ 2 B 1 α ( m ) w 1 m ^ ( m ) + m ¯ 2 D 1 α ( m ) w 2 m ^ ( m ) + m ¯ 2 F 1 α ( m ) w 3 m ^ ( m ) ,
δ w 1 ( m ) : m ¯ A 13 ( m ) u 0 m ^ ( m ) + m ¯ B 13 ( m ) + m ¯ B 55 ( m ) u 1 m ^ ( m ) + m ¯ D 13 ( m ) + 2 m ¯ D 55 ( m ) u 2 m ^ ( m ) + m ¯ F 13 ( m ) + 3 m ¯ F 55 ( m ) u 3 m ^ ( m ) + m ¯ 2 B 55 ( m ) w 0 m ^ ( m ) + A 13 ( m ) + m ¯ 2 D 55 ( m ) w 1 m ^ ( m ) + 2 B 33 ( m ) + m ¯ 2 F 55 ( m ) w 2 m ^ ( m ) + 3 D 33 ( m ) + m ¯ 2 H 55 ( m ) w 3 m ^ ( m ) + h m / 2 λ x m ^ m 1 + λ x m ^ ( m ) = Δ T m ¯ 2 B 1 α ( m ) w 0 m ^ ( m ) + m ¯ 2 D 1 α ( m ) w 1 m ^ ( m ) + m ¯ 2 F 1 α ( m ) w 2 m ^ ( m ) + m ¯ 2 H 1 α ( m ) w 3 m ^ ( m ) ,
δ w 2 ( m ) : 2 m ¯ B 13 ( m ) u 0 m ^ ( m ) + 2 m ¯ D 13 ( m ) + m ¯ D 55 ( m ) u 1 m ^ ( m ) + 2 m ¯ F 13 ( m ) + 2 m ¯ F 55 ( m ) u 2 m ^ ( m ) + 2 m ¯ H 13 ( m ) + 3 m ¯ H 55 ( m ) u 3 m ^ ( m ) + m ¯ 2 D 55 ( m ) w 0 m ^ ( m ) + 2 B 33 ( m ) + m ¯ 2 F 55 ( m ) w 1 m ^ ( m ) + 4 D 33 ( m ) + m ¯ 2 H 55 ( m ) w 2 m ^ ( m ) + 6 F 33 ( m ) + m ¯ 2 J 55 ( m ) w 3 m ^ ( m ) + h m 2 / 4 λ z m ^ m 1 λ z m ^ ( m ) = Δ T m ¯ 2 D 1 α ( m ) w 0 m ^ ( m ) + m ¯ 2 F 1 α ( m ) w 1 m ^ ( m ) + m ¯ 2 H 1 α ( m ) w 2 m ^ ( m ) + m ¯ 2 J 1 α ( m ) w 3 m ^ ( m ) ,
δ w 3 ( m ) : 3 m ¯ D 13 ( m ) u 0 m ^ ( m ) + 3 m ¯ F 13 ( m ) + m ¯ F 55 ( m ) u 1 m ^ ( m ) + 3 m ¯ H 13 ( m ) + 2 m ¯ H 55 ( m ) u 2 m ^ ( m ) + 3 m ¯ J 13 ( m ) + 3 m ¯ J 55 ( m ) u 3 m ^ ( m ) + m ¯ 2 F 55 ( m ) w 0 m ^ ( m ) + 3 D 33 ( m ) + m ¯ 2 H 55 ( m ) w 1 m ^ ( m ) + 6 F 33 ( m ) + m ¯ 2 J 55 ( m ) w 2 m ^ ( m ) + 9 H 33 ( m ) + m ¯ 2 L 55 ( m ) w 3 m ^ ( m ) + h m 3 / 8 λ z m ^ m 1 + λ z m ^ ( m ) = Δ T m ¯ 2 F 1 α ( m ) w 0 m ^ ( m ) + m ¯ 2 H 1 α ( m ) w 1 m ^ ( m ) + m ¯ 2 J 1 α ( m ) w 2 m ^ ( m ) + m ¯ 2 L 1 α ( m ) w 3 m ^ ( m ) ,
δ λ x k : u 0 m ^ ( k + 1 ) h k + 1 / 2 u 1 m ^ ( k + 1 ) + h k + 1 2 / 4 u 2 m ^ ( k + 1 ) h k + 1 3 / 8 u 3 m ^ ( k + 1 ) u 0 m ^ ( k ) + h k / 2 u 1 m ^ ( k ) + h k 2 / 4 u 2 m ^ ( k ) + h k 3 / 8 u 3 m ^ ( k ) = 0 ,
δ λ z ( k ) : w 0 m ^ ( k + 1 ) h k + 1 / 2 w 1 m ^ ( k + 1 ) + h k + 1 2 / 4 w 2 m ^ ( k + 1 ) h k + 1 3 / 8 w 3 m ^ ( k + 1 ) w 0 m ^ ( k ) + h k / 2 w 1 m ^ ( k ) + h k 2 / 4 w 2 m ^ ( k ) + h k 3 / 8 w 3 m ^ ( k ) = 0 ,
where m ¯ = m ^ π / L , m = 1, 2,…, Nl, and k = 1, 2,…, (Nl − 1).
Equations (48)–(57) are a system of simultaneously linear algebraic equations representing a standard eigen-valued problem, in which the eigenvalues are the critical temperature changes, which can be obtained by letting the determinant of the coefficient matrix of these simultaneously linear algebraic equations be zero. Once the eigenvalues are obtained, their corresponding eigen-vectors (i.e., buckling mode shapes) can also be determined.

4. Optimal Design

4.1. Statement of the Optimization Problem

In this work, the authors consider the multi-objective optimization of the volume fractions of the constituents of FG beams subjected to a uniform temperature change in order to maximize the critical temperature change parameter of the FG beam and to minimize its total mass. The FG beam is considered to be formed by mixing a metal material and a ceramic material according to a three-parameter power-law function of the volume fractions of the constituents through the thickness coordinate of the FG beam, which is given as follows:
The three-parameter power-law function,
Γ c = [ 1 2 + z h + κ b 1 2 z h κ c ] κ p   and   Γ c + Γ m = 1 ,
where the symbols κ b , κ c , and κ p denote the material-property gradient indices. The material at the top surface of the FG beam (i.e., z = h/2) is ceramic rich, and when κ b = 0 , the material at the bottom surface of the FG beam (i.e., z = −h/2) is metal rich. Variations in the through-thickness distributions of the volume fraction given in Equation (58) with varying one of κ b , κ c , and κ p and holding the other two indices constants are shown in Figure 2.
In addition, when κ b = 0 , the three-parameter power-law function is reduced to the one-parameter power-law function, which is commonly used in the literature and is given as follows:
One-parameter power-law function,
Γ c = 1 2 + z h κ p   and   Γ c + Γ m = 1 ,
which is used in the thermal buckling cases in this work for a comparative study. Variations in the through-thickness distributions of the volume fraction given in Equation (59) with varying the κ p value are shown in Figure 3.
In the optimal design, a mass ratio R m is defined as follows:
R m = ρ ^ f ρ ^ c / ρ ^ m ρ ^ c ,
where ρ ^ f , ρ ^ c , and ρ ^ m are the total mass per unit area in x-y plane of the FG beam, the homogeneous ceramic material, and the homogeneous metal material, respectively. ρ ^ k = h / 2 h / 2 ρ k z d z , in which k = f, c, or m.
The critical temperature parameter Δ T ^ c r is defined as follows:
Δ T ^ c r = Δ T c r α c L / h 2 ,
where α c denotes the thermal expansion coefficient of a reference ceramic material.
The objective functions are defined as follows:
Objective   function   1 :   F 1 = R m ,
Objective   function   2 :   F 2 = 1 { [ Δ T c r Δ T c r c ] / Δ T c r m Δ T c r c } ,
where the ranges of F1 and F2 are 0 F 1 or F 2 1 .
Because non-dominated sorting is used in the current GA, the values of these objective functions given in Equations (62) and (63) are, thus, used to classify each design into its corresponding non-dominated front, which is also the assigned fitness value used for the sorting process. Minimization of the fitness function can be accomplished when the critical temperature change parameters of the FG beam are obtained. The detailed process of the current non-dominated sorting-based GA is described in the following sections.

4.2. The Non-Dominated Sorting-Based GA

In this work, the non-dominated sorting-based GA [69,70] is used as an optimal technique, the usual form of which is described as follows: The GA starts with an initial set of random solutions, which are called a population. Each individual in the population is called a chromosome, which represents a solution to the problem at hand. The chromosomes are made of discrete units called genes (or design variables). Each gene controls one or more features of the chromosome. The chromosomes evolve through successive iterations, called generations. Based on some measures of fitness, the chromosomes during each generation are evaluated. To create the next generation, new chromosomes in the next generation, called offspring, are formed using two operators, which are the crossover and mutation operators. In the former, some portions of two chromosomes in the current generation are merged together, and in the latter, some portions of chromosomes in the current generation are modified. Thus, the crossover operator leads the population to converge by means of making the chromosome in the population alike, and the mutation operator assists the search escaping from local optima by reintroducing genetic diversity back into the population. A new generation is formed by selecting and rejecting some of the parents and offspring according to their fitness values so as to keep the population size of each generation constant. After several generations, the GA converges to the best chromosome, representing the optimum solution to the problem considered. A flow chart of the non-dominated sorting-based GA is shown in Figure 4, and its related process is described as follows:
(a)
Initial population
The first population consists of chromosomes (i.e., solutions or designs), which is taken to be 200 in this work. Each chromosome is formed by several genes (i.e., design variables), each of which is represented by a randomly generated real number with four decimal places in a specific range.
(b)
Fitness value
As mentioned before, there are two objective functions considered in this work, the critical temperature change parameter of the FG beam and its total mass, the specific forms of which (i.e., F1 and F2) are, thus, defined as the fitness values and are used for the sorting process in the current GA.
(c)
Selection
In this work, the selection process consists of two approaches, the Pareto-ranking approach [55] and the crowding distance approach [69], which are described as follows:
(c-1) Pareto-ranking approach
The Pareto-ranking approach uses the Pareto dominance concept to evaluate fitness or assign selection probability to solutions. Each solution in the population of the parents is ranked according to a dominance rule, and then each solution is assigned a fitness value based on its rank in the population, rather than its actual objective function value. A dominance rule proposed by Goldberg [55] is adopted and is described as follows:
Step 1: Set i = 1 and n p 0 = the number of original solutions in the population.
Step 2: Identify non-dominated solutions xj (j = 1~ n n d f i ) in n p i , and assign them set to N D F i , in which N D F i denotes the ith non-dominated front, and the number of solutions in N D F i is n n d f i .
Step 3: Set i = i + 1, and then n p i + 1 = n p i 1 n p i . If n p i + 1 = 0 , go to Step 4; otherwise, go to Step 2.
Step 4: For every solution x k , k = 1 ~ n p 0 at generation t, assign rank r x k , t = i if x k N D F i .
In the above process, a solution A dominates the other solution B, which means both F 1 x A > F 1 x B and F 2 x A > F 2 x B when two objective functions are considered. Because all objective functions are to be minimized, a lower rank corresponds to a better solution, such that the population will be classified as several non-dominated fronts, where the first non-dominated front is also called the Pareto fronts of this population.
(c-2) Crowding distance approach
A crowding distance approach as proposed by Deb [69] is adopted in this work, which is aimed toward obtaining a uniform spread of solutions along the best-known Pareto front. The approach is described as follows:
Step 1: Rank the population and identify Pareto fronts as P F 1 , P F 2 , , P F R , using Equations (62) and (63). For each non-dominated front, repeat Steps 2 and 3.
Step 2: For each objective function F k , sort the solutions in P F j in the ascending order. If the number of solutions in P F j is l, then the crowding distance of the ith solution is measured as
c d k x i = F k x i + 1 F k x i 1 / F k max F k min
where i = 2, 3, …, (l−1), and c d k x 1 = c d k x l = .
Step 3: To obtain the total crowding distance c d x i = k = 1 2 c d k x i .
The selection rule is based on the rank assigned to the solutions, where a lower rank corresponds to a better solution. When the solutions are in the same rank, a greater crowding distance corresponds to a better solution.
(d)
Crossover
Any two solutions in the population are randomly paired for mating, in which a simulated binary crossover operator suggested by Deb and Goyal [78] is adopted in this work and is given as follows:
x j 1 t + 1 = 0.5 1 + γ j x j 1 t + 1 γ j x j 2 t ,
x j 2 t + 1 = 0.5 1 γ j x j 1 t + 1 + γ j x j 2 t ,
where γ j = 2 u j 1 / η j + 1 when u j 0.5 , and γ j = 1 / 2 1 u j 1 / η j + 1 when u j 0.5 , in which ui denotes a randomly generated real number between 0 and 1 and η > 0 is a user-defined index parameter, the value of which is suggested to be η j = 1 .
(e)
Mutation
Mutation occurs in each solution in the population when the randomly generated number between 0 and 1 for each solution is less than 0.1. A polynomial mutation operator suggested by Deb and Goyal [78] is adopted in this work and is given as follows:
x i t + 1 = x i t + β i ,
where β i = 2 u i 1 / η i + 1 1 when u i 0.5 and β i = 1 2 1 u i 1 / η i + 1 when u i 0.5 , in which ui denotes a randomly generated real number between 0 and 1 and η i is a user-defined index parameter, for which the value is suggested to be η i = 20 .
(f)
Elitist
Performing the above operators, including the non-dominated sorting, crowding distance, crossover, and mutation approaches, we can obtain the new generation (i.e., offspring) with the same number of solutions (i.e., np) as those in the original generation (i.e., parents). The elitist strategy in this work is as follows:
Step 1: Combine the solutions of parents and offspring to form a population with 2 np solutions.
Step 2: Perform the non-dominated sorting and crowding distance approaches to each solution of the population generated in Step 1. Then, select the first np solutions with smaller fitness values in the population to form the offspring.
(g)
End criteria
The end criteria of the current non-dominated sorting-based GA is given as
Criteria: when the number of generations reaches an assigned number, which is 100 in this work.

5. Illustrative Examples

5.1. Thermal Buckling Analysis of Laminated Composite Beams and FG Beams

Because no benchmark solutions for the critical temperature changes for a simply-supported FG beam can be found in the literature, for comparison purposes, the current mixed LW HSDT is applied to the thermal buckling analysis of simply-supported, cross-ply laminated beams under a uniform temperature, in which the material properties are regarded as a layer-wise constant distribution varying through the thickness coordinate. The current solutions for critical temperature change parameters are compared with those obtained using the EBT [79], FSDT [79], TSDT [79], HSDT [80], and HRSDT [22]. The material properties and geometric parameters are given as L / h = 10   and   50 , E 1 / E 2 = 10 , 20 , and 30 , G 12 / E 2 = G 13 / E 2 = 0.6 , υ 12 = υ 13 = 0.25 , α 2 / α 1 = 3 , in which the subscript 1 represents the material properties parallel to the reinforced fiber direction, and the subscripts 2 and 3 represent material properties perpendicular to the reinforced fiber direction. A dimensionless temperature change parameter is defined as Δ T ¯ c r = Δ T c r α 1 L / h 2 .
Table 1 and Table 2 show the convergent study for critical temperature change parameter solutions of [0°/90°/0°] laminated beams subjected to a uniform temperature change. It can be seen in Table 1 and Table 2 that the current LW1, LW2, and LW3 solutions converge rapidly. For a moderately thick beam (L/h = 10), the convergent solutions are yielded when the total number of layers are taken to be 9, 6, and 3 for the LW1, LW2, and LW3 theories, respectively, and these convergent solutions closely agree with one another. The convergence rate increases when the beam becomes thinner. The current solutions in Table 1 and Table 2 with a superscript a indicate that the von Kármán GNS, rather than the full GNS, is considered. The results show that the deviations between the solutions considering the von Kármán GNS and the full GNS are less than 1%, which is very minor, such that the von Kármán GNS is suitable for the current analysis of moderately thick laminated composite beams. The current solutions in Table 1 and Table 2 with a superscript b indicate that a narrow beam is considered, in which Poisson’s ratios are taken to be zeroes in the stiffness coefficients Qij (i, j = 1–3). The results show that the deviations between the solutions considering and those not-considering the narrow beam effects are significant, where the deviations as high as approximately to 5.5% for both a moderately thick (L/h = 10) and a thin LC beam (L/h = 50). Table 1 and Table 2 also show that the current convergent solutions considering the narrow beam effects closely agree with the HSRDT [22] and TSDT [79] solutions, and those not considering the narrow beam effects closely agree with the HSDT solutions [80]. In addition, the EBT solutions for the cases of L/h = 50 are the same as those for the cases of L/h = 10. In the former, the relative error between the EBT solutions and the current LW HSDT solutions is 1.2%, while it is 31.2% for the latter, such that the EBT is suitable for the current analysis of very thin laminated composite beams only.
The thermal buckling problems of simply-supported FG beams with either the TD or the TI material properties under a uniform temperature change are considered in Table 3. A flow chart for the solution process to determine the critical temperature change parameters of the FG beam is given in Figure 5 where the TD material properties are considered, in which the golden section method is used. The through-thickness distributions of the material properties of the FG beam are assumed to obey a single-parameter power-law function according to the volume fractions of the constituents, as shown in Equation (57). The effective material properties are estimated using either the rule of mixtures or the Mori–Tanaka scheme, which are given in Equations (1) and (2), respectively. Again, it can be seen in Table 3 that the current LW HSDT solutions converge rapidly. The solutions for the critical temperature change parameters obtained using the Mori–Tanaka scheme are approximately no more than 2% greater than those obtained using the rule of mixtures. The effects of different micromechanics models on the critical temperature change parameters of the FG beam are minor.
The results also show that the critical temperature change parameters obtained using the TI material properties always are greater than those obtained using the TD material properties. For a moderately thick beam (L/h = 10 and κ p = 1 ), deviations between the critical temperature change parameters obtained using the TI and the TD material properties are as high as approximately 36.5%, such that the current analysis using TI material properties is an unsafe analysis and, thus, is not recommended. The effects of TD material properties on the critical temperature change parameters are significant, and these must be considered for the current issue.

5.2. Optimization of Material Composition of FG beams

In this section, a multiple-objective optimization for material composition of a simply-supported FG beam subjected to a uniform temperature change is considered in order to maximize its critical temperature change and minimize its self-weight. In the optimal design, the FG beam is considered to be a two-phase composite material, where one phase is the ceramic material (ZrO2) and the other is the metal material (SUS304). The TD material properties of ZrO2 and SUS304 are given in Table 4 [4]. Variations in the material properties of ZrO2 and SUS304 with the temperature variable ranging from 300 K to 1100 K are shown in Figure 6. The material properties of the FG beam are assumed to obey a three-parameter power-law distribution of volume fractions of the constituents along the thickness of the FG beam, and the effective material properties are estimated using the rule of mixtures. The non-dominated sorting-based GA mentioned in Section 4.2 is used to determine some sets of Pareto-optimal solutions of the undetermined coefficients κ p , κ b , and κ c , in which 200 initial populations are randomly generated, in which the ranges of κ p , κ b , and κ c are taken as 0 κ p 50 , 0 κ b 1 , and 1 κ c 3 , respectively.
Figure 7 and Figure 8 show the populations at the initial, first, second, fifth, tenth, and 20th generations, where the TI and TD material properties are considered, respectively. It can be seen in Figure 7 and Figure 8 that the non-dominated, sorting-based GA converges rapidly and that the Poreto-optimal solutions can be yielded after 20 generations. Curve fitting for these Preto-optimal solutions is performed for the cases of TI and TD material properties, which is shown in Figure 9 with the solid and dashed lines, respectively, in which
F 2 = 1.3303 F 1 4 + 1.1843 F 1 3 + 0.4537 F 1 2 1.3 F 1 + 1.0005   for   the   TI   material   property   cases ,
F 2 = 1.8376 F 1 4 + 2.4524 F 1 3 1.2222 F 1 2 0.3734 F 1 + 0.9936 for   the   TD   material   property   cases .
Equations (68) and (69) can then be used to calculate the associated weight numbers (i.e., w 1 and w 2 ) for the objective functions F1 and F2, in which if we let d F 2 / d F 1 = λ , then w 1 = λ / 1 λ and w 2 = 1 / 1 λ .
Table 5 and Table 6 show every ten Pareto-optimal solutions sorted by the F2 function values from the largest to the smallest for the TI and TD material property cases, respectively. These Pareto-optimal solutions may provide design engineers with valuable information regarding what set of the material-property gradient indices ( κ p , κ b , and κ c ) they need according to the weight number ratio of (w2/w1). In addition, the results also show that the critical temperature changes of the FG beam for the TD material property cases are much less than those for the TI material property cases, such that the TD material property effects must be considered in TD problems.

6. Conclusions

In this work, the authors developed a mixed LW HSDT for the thermal buckling analysis of FG beams subjected to a uniform temperature change, and then they further developed a non-dominated sorting-based GA for multi-objectives optimization of the material composition of a three-parameter FG beam, in which the TI and TD material properties are considered. All of the results in the tables and figures are obtained using the self-developed MATLAB programs without any commercial toolbox.
In the thermal buckling analysis, the results show that the TD material properties must be considered in TD physical problems, in which for a moderately thick FG beam, the deviations in the critical temperature change parameters obtained using the TD material properties and the TI material properties are as high as approximately 36.5%, and the analysis based on the TI material materials is an unsafe analysis, due to the fact that it always overestimates the actual critical temperature change parameters.
In the optimal design cases, the results show the self-developed non-dominated sorting-based GA converges rapidly, where the Poreto-optimal solutions can be yielded after the 20th generation. The non-dominated sorting-based GA can also be extended to other optimal designs of FG beams with multiple objective functions.

Author Contributions

Conceptualization, C.-P.W.; methodology, C.-P.W.; software, K.-W.L.; validation, C.-P.W. and K.-W.L.; investigation: C.-P.W. and K.-W.L.; resources, C.-P.W.; data curation, K.-W.L.; writing—original draft preparation, C.-P.W.; writing—review and editing, C.-P.W.; visualization, C.-P.W.; supervision, C.-P.W.; project administration, C.-P.W. funding acquisition, C.-P.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Technology, Taiwan, grant number MOST 109-2221-E-006-015-MY3.

Institutional Review Board Statement

This study did not involve humans and animals.

Informed Consent Statement

This study did not involve humans and animals.

Data Availability Statement

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Relations between the Generalized Force/Moment Resultants and the Generalized Displacements

The definition of generalized force and moment resultants and their relations with the displacement components are given as follows:
N x ( m ) = h m / 2 h m / 2 σ x ( m ) d z m = A 11 ( m ) u 0 ( m ) , x + B 11 ( m ) u 1 ( m ) , x + D 11 ( m ) u 2 ( m ) , x + F 11 ( m ) u 3 ( m ) , x + A 13 ( m ) w 1 ( m ) + 2 B 13 ( m ) w 2 ( m ) + 3 D 13 ( m ) w 3 ( m ) ,
M x ( m ) = h m / 2 h m / 2 σ x ( m ) z m d z m = B 11 ( m ) u 0 ( m ) , x + D 11 ( m ) u 1 ( m ) , x + F 11 ( m ) u 2 ( m ) , x + H 11 ( m ) u 3 ( m ) , x + B 13 ( m ) w 1 ( m ) + 2 D 13 ( m ) w 2 ( m ) + 3 F 13 ( m ) w 3 ( m ) ,
P x ( m ) = h m / 2 h m / 2 σ x ( m ) z m 2 d z m = D 11 ( m ) u 0 ( m ) , x + F 11 ( m ) u 1 ( m ) , x + H 11 ( m ) u 2 ( m ) , x + J 11 ( m ) u 3 ( m ) , x + D 13 ( m ) w 1 ( m ) + 2 F 13 ( m ) w 2 ( m ) + 3 H 13 ( m ) w 3 ( m ) ,
R x ( m ) = h m / 2 h m / 2 σ x ( m ) z m 3 d z m = F 11 ( m ) u 0 ( m ) , x + H 11 ( m ) u 1 ( m ) , x + J 11 ( m ) u 2 ( m ) , x + L 11 ( m ) u 3 ( m ) , x + F 13 ( m ) w 1 ( m ) + 2 H 13 ( m ) w 2 ( m ) + 3 J 13 ( m ) w 3 ( m ) ,
N z ( m ) = h m / 2 h m / 2 σ z ( m ) d z m = A 13 ( m ) u 0 ( m ) , x + B 13 ( m ) u 1 ( m ) , x + D 13 ( m ) u 2 ( m ) , x + F 13 ( m ) u 3 ( m ) , x + A 33 ( m ) w 1 ( m ) + 2 B 33 ( m ) w 2 ( m ) + 3 D 33 ( m ) w 3 ( m ) ,
M z ( m ) = h m / 2 h m / 2 σ z ( m ) z m d z m = B 13 ( m ) u 0 ( m ) , x + D 13 ( m ) u 1 ( m ) , x + F 13 ( m ) u 2 ( m ) , x + H 13 ( m ) u 3 ( m ) , x + B 33 ( m ) w 1 ( m ) + 2 D 33 ( m ) w 2 ( m ) + 3 F 33 ( m ) w 3 ( m ) ,
P z ( m ) = h m / 2 h m / 2 σ z ( m ) z m 2 d z m = D 13 ( m ) u 0 ( m ) , x + F 13 ( m ) u 1 ( m ) , x + H 13 ( m ) u 2 ( m ) , x + J 13 ( m ) u 3 ( m ) , x + D 33 ( m ) w 1 ( m ) + 2 F 33 ( m ) w 2 ( m ) + 3 H 33 ( m ) w 3 ( m ) ,
N x z ( m ) = h m / 2 h m / 2 τ x z ( m ) d z m = A 55 ( m ) w 0 ( m ) , x + u 1 ( m ) + B 55 ( m ) w 1 ( m ) , x + 2 u 2 ( m ) + D 55 ( m ) w 2 ( m ) , x + 3 u 3 ( m ) + F 55 ( m ) w 3 ( m ) , x ,
M x z ( m ) = h m / 2 h m / 2 τ x z ( m ) z m d z m = B 55 ( m ) w 0 ( m ) , x + u 1 ( m ) + D 55 ( m ) w 1 ( m ) , x + 2 u 2 ( m ) + F 55 ( m ) w 2 ( m ) , x + 3 u 3 ( m ) + H 55 ( m ) w 3 ( m ) , x ,
P x z ( m ) = h m / 2 h m / 2 τ x z ( m ) z m 2 d z m = D 55 ( m ) w 0 ( m ) , x + u 1 ( m ) + F 55 ( m ) w 1 ( m ) , x + 2 u 2 ( m ) + H 55 ( m ) w 2 ( m ) , x + 3 u 3 ( m ) + J 55 ( m ) w 3 ( m ) , x ,
R x z ( m ) = h m / 2 h m / 2 τ x z ( m ) z m 3 d z m = F 55 ( m ) w 0 ( m ) , x + u 1 ( m ) + H 55 ( m ) w 1 ( m ) , x + 2 u 2 ( m ) + J 55 ( m ) w 2 ( m ) , x + 3 u 3 ( m ) + L 55 ( m ) w 3 ( m ) , x ,
where A i j ( m ) B i j ( m ) D i j ( m ) F i j ( m ) H i j ( m ) J i j ( m ) L i j ( m ) T = h m / 2 h m / 2 Q i j ( m ) 1 z m z m 2 z m 3 z m 4 z m 5 z m 6 T d z m , in which the superscript T refers to the transport of the vector.

References

  1. Miyamoto, Y.; Kaysser, W.A.; Rabin, B.H.; Kawasaki, A.; Ford, R.G. Functionally Graded Materials: Design, Processing, and Applications; Springer: Manhattan, NY, USA, 1999. [Google Scholar]
  2. Reddy, J.N. Analysis of functionally graded plates. Int. J. Numer. Methods Eng. 2000, 47, 663–684. [Google Scholar] [CrossRef]
  3. Reddy, J.N.; Chin, C.D. Thermomechanical analysis of functionally graded cylinders and plates. J. Therm. Stresses 1998, 21, 593–626. [Google Scholar] [CrossRef]
  4. Shen, H.S. Functionally Graded Materials: Nonlinear Analysis of Plates and Shells; CRC Press: Boca Raton, FL, USA, 2009. [Google Scholar]
  5. Koizumi, M. The concept of FGM. Ceram. Tran. 1993, 34, 3–10. [Google Scholar]
  6. Koizumi, M. FGM activities in Japan. Compos. Part B 1997, 28, 1–4. [Google Scholar] [CrossRef]
  7. Jha, D.K.; Kant, T.; Singh, R.K. A critical review of recent research on functionally graded plates. Compos. Struct. 2013, 96, 833–849. [Google Scholar] [CrossRef]
  8. Liew, K.M.; Ferreira, A.J.M. A review of meshless methods for laminated and functionally graded plates and shells. Compos. Struct. 2011, 93, 2013–2041. [Google Scholar] [CrossRef]
  9. Punera, D.; Kant, T. A critical review of stress and vibration analyses of functionally graded shell structures. Compos. Struct. 2019, 210, 787–809. [Google Scholar] [CrossRef]
  10. Swaminathan, K.; Naveenkumar, D.T.; Zenkour, A.M.; Carrera, E. Stress, vibration and buckling analyses of FGM plates-A state-of-the-art review. Compos. Struct. 2015, 120, 10–31. [Google Scholar] [CrossRef]
  11. Sayyad, A.S.; Ghugal, Y.M. Modeling and analysis of functionally graded sandwich beams: A review. Mech. Adv. Mater. Struct. 2019, 26, 1776–1795. [Google Scholar] [CrossRef]
  12. Thai, H.T.; Kim, S.E. A review of theories for the modeling and analysis of functionally graded plates and shells. Compos. Struct. 2015, 128, 70–86. [Google Scholar] [CrossRef]
  13. Wu, C.P.; Chiu, K.H.; Wang, Y.M. A review on the three-dimensional analytical approaches of multilayered and functionally graded piezoelectric plates and shells. CMC-Comput. Mater. Continua 2008, 8, 93–132. [Google Scholar]
  14. Wu, C.P.; Liu, Y.C. A review of semi-analytical numerical methods for laminated composite and multilayered functionally graded elastic/piezoelectric plates and shells. Compos. Struct. 2016, 147, 1–15. [Google Scholar] [CrossRef]
  15. Kiani, Y.; Eslami, M.R. Thermal buckling analysis of functionally graded material beams. Int. J. Mech. Mater. Des. 2010, 6, 229–238. [Google Scholar] [CrossRef]
  16. Tauchert, T.R. Thermal buckling of antisymmetric angle-ply laminates. J. Therm. Stresses 1987, 10, 113–124. [Google Scholar] [CrossRef]
  17. Reddy, J.N. A simple higher-order theory for laminated composite plates. J. Appl. Mech. 1984, 51, 745–752. [Google Scholar] [CrossRef]
  18. Tran, T.T.; Nguyen, N.H.; Do, T.V.; Minh, P.V.; Duc, N.D. Bending and thermal buckling of unsymmetric functionally graded sandwich beams in high-temperature environment based on a new third-order shear deformation theory. J. Sandw. Struct. Mater. 2019. [Google Scholar] [CrossRef]
  19. Wattanasakulpong, N.; Prusty, B.G.; Kelly, D.W. Thermal buckling and elastic vibration of third-order shear deformable functionally graded beams. Int. J. Mech. Sci. 2011, 53, 734–743. [Google Scholar] [CrossRef]
  20. Soldatos, K.P.; Elishakoff, I. A transverse shear and normal deformable orthotropic beam theory. J. Sound Vibr. 1992, 155, 528–533. [Google Scholar] [CrossRef]
  21. Touratier, M. An efficient standard plate theory. Int. J. Eng. Sci. 1991, 29, 901–916. [Google Scholar] [CrossRef]
  22. Bouazza, M.; Benseddiq, N.; Zenkour, A.M. Thermal buckling analysis of laminated composite beams using hyperbolic refined shear deformation theory. J. Therm. Stresses 2019, 42, 332–340. [Google Scholar] [CrossRef]
  23. Nguyen, T.K.; Nguyen, T.T.P.; Vo, T.P.; Thai, H.T. Vibration and buckling analysis of functionally graded sandwich beams by a new higher-order shear deformation theory. Compos. Part B 2015, 76, 273–285. [Google Scholar] [CrossRef] [Green Version]
  24. Mantari, J.L.; Oktem, A.S.; Soares, C.G. A new higher order shear deformation theory for sandwich and composite laminated plates. Compos. Part B 2012, 43, 1489–1499. [Google Scholar] [CrossRef]
  25. Liu, N. An isogeometric continuum shell element for modeling the nonlinear response of functionally graded material structures. Compos. Struct. 2020, 237, 111893. [Google Scholar] [CrossRef]
  26. Liu, N.; Plucinsky, P.; Jeffers, A.E. Combining load-controlled and displacement-controlled algorithms to model thermal-mechanical snap-through instabilities in structures. J. Eng. Mech. 2017, 143, 04017051. [Google Scholar] [CrossRef]
  27. Refrafi, S.; Bousahla, A.A.; Bousahla, A.; Menasria, A.; Bourada, F.; Tounsi, A.; Bedia, E.A.A.; Mahmoud, S.R.; Benrahou, K.H.; Tounsi, A. Effects of hygro-thermo-mechanical conditions on the buckling of FG sandwich plates resting on elastic foundations. Comput. Concrete 2020, 25, 311–325. [Google Scholar]
  28. Carrera, E. Theories and finite elements for multilayered plates and shells: A unified compact formulation with numerical assessment and benchmarking. Arch. Comput. Methods Eng. 2003, 10, 215–296. [Google Scholar] [CrossRef]
  29. Carrera, E.; Giunta, G.; Petrolo, M. Beam Structures: Classical and Advanced Theories; John Wiley and Sons: Hoboken, NJ, USA, 2011. [Google Scholar]
  30. She, G.L.; Yuan, F.G.; Ren, Y.R. Thermal buckling and post-buckling analysis of functionally graded beams based on a general higher-order shear deformation theory. Appl. Math. Modell. 2017, 47, 340–357. [Google Scholar] [CrossRef]
  31. Thai, H.T.; Vo, T.P. Bending and free vibration of functionally graded beams using various higher-order shear deformation beam theories. Int. J. Mech. Sci. 2012, 62, 57–66. [Google Scholar] [CrossRef] [Green Version]
  32. Aydogdu, M. Thermal buckling analysis of cross-ply laminated composite beams with general boundary conditions. Compos. Sci. Tech. 2007, 67, 1096–1104. [Google Scholar] [CrossRef]
  33. Simsek, M. Fundamental frequency analysis of functionally graded beams by using different higher-order beam theories. Nucl. Eng. Des. 2010, 240, 697–705. [Google Scholar] [CrossRef]
  34. Timarci, T.; Soldatos, K.P. Comparative dynamic studies for symmetric cross-ply circular cylindrical shells on the basis of a unified shear deformable shell theory. J. Sound Vib. 1995, 187, 609–624. [Google Scholar] [CrossRef]
  35. Tahani, M. Analysis of laminated composite beams using layerwise displacement theories. Compos. Struct. 2007, 79, 535–547. [Google Scholar] [CrossRef]
  36. Lee, H.J.; Saravanos, D.A. Coupled layerwise analysis of thermopiezoelectric composite beams. AIAA J. 1996, 34, 1231–1237. [Google Scholar] [CrossRef]
  37. Shimpi, R.P.; Ainapure, A.V. A beam finite element based on layerwise trigonometric shear deformation theory. Compos. Struct. 2001, 53, 153–162. [Google Scholar] [CrossRef]
  38. Shimpi, R.P.; Ghugal, Y.M. A new layerwise trigonometric shear deformation theory for two-layered cross-ply beams. Compos. Sci. Technol. 2001, 61, 1271–1283. [Google Scholar] [CrossRef]
  39. Pandey, S.; Pradyumna, S. Analysis of functionally graded sandwich plates using a higher-order layerwise theory. Compos. Part B 2018, 153, 325–336. [Google Scholar] [CrossRef]
  40. Bayat, Y.; Toussi, H.E. Exact solution of thermal buckling and post buckling of composite and SMA hybrid composite beam by layerwise theory. Aerosp. Sci. Technol. 2017, 67, 484–494. [Google Scholar] [CrossRef]
  41. Washizu, K. Variational Methods in Elasticity and Plasticity; Pergamon Press: New York, NY, USA, 1982. [Google Scholar]
  42. Reddy, J.N. Energy and Variational Methods in Applied Mechanics: With an Introduction to the Finite Element Method; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 1984. [Google Scholar]
  43. Wu, C.P.; Kuo, H.C. Interlaminar stresses analysis for laminated composite plates based on a local high order lamination theory. Compos. Struct. 1992, 20, 237–247. [Google Scholar] [CrossRef]
  44. Wu, C.P.; Chen, W.Y. Vibration and stability of laminated plates based on a local high order plate theory. J. Sound Vib. 1994, 177, 503–520. [Google Scholar] [CrossRef]
  45. Wu, C.P.; Xu, Z.R. Strong and weak formulations of a mixed higher-order shear deformation theory for the static analysis of functionally graded beams under thermo-mechanical loads. J. Compos. Sci. 2020, 4, 158. [Google Scholar] [CrossRef]
  46. Liew, K.M.; Pan, Z.Z.; Zhang, L.W. An overview of layerwise theories for composite laminates and structures: Development, numerical implementation and application. Compos. Struct. 2019, 216, 240–259. [Google Scholar] [CrossRef]
  47. Walker, M.; Reiss, T.; Adali, S. A procedure to select the best material combinations and optimally design hybrid composite plates for minimum weight and cost. Eng. Opt. 1997, 29, 65–83. [Google Scholar] [CrossRef]
  48. Walker, M.; Smith, R. A procedure to select the best material combinations and optimally design composite sandwich cylindrical shells for minimum mass. Mater. Des. 2006, 27, 160–165. [Google Scholar] [CrossRef]
  49. Houmat, A. Optimal lay-up design of variable stiffness laminated composite plates by a layer-wise optimization technique. Eng. Opt. 2018, 50, 205–217. [Google Scholar] [CrossRef]
  50. Guenanou, A.; Houmat, A. Optimum stacking sequence design of laminated composite circular plates with curvilinear fibres by a layer-wise optimization method. Eng. Opt. 2018, 50, 766–780. [Google Scholar] [CrossRef]
  51. Cho, J.R.; Ha, D.Y. Volume fraction optimization for minimizing thermal stress in Ni-Al2O3 functionally graded materials. Mater. Sci. Eng. A 2002, 334, 147–155. [Google Scholar] [CrossRef]
  52. Vlasov, V.Z. Thin-Walled Elastic Beams; Program for Scientific Translation: Jerusalem, Israel, 1961. [Google Scholar]
  53. Nguyen, T.T.; Lee, J. Optimal design of thin-walled functionally graded beams for buckling problems. Compos. Struct. 2017, 179, 459–467. [Google Scholar] [CrossRef]
  54. Gen, M.; Cheng, R. Genetic Algorithms & Engineering Design; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 1997. [Google Scholar]
  55. Goldberg, D.E. Genetic Algorithms in Search, Optimization, and Machine Learning; Addison-Wesley Publishing Company, Inc.: White Plains, NY, USA, 1989. [Google Scholar]
  56. Karama, M.; Afaq, K.S.; Mistou, S. Mechanical behavior of laminated composite beam by the new multi-layered laminated composite structures model with transverse shear stress continuity. Int. J. Solids Struct. 2003, 40, 1525–1546. [Google Scholar] [CrossRef]
  57. Yas, M.H.; Kamarian, S.; Jam, J.E.; Pourasghar, A. Optimization of functionally graded beams resting on elastic foundations. J Solid Mech. 2011, 3, 365–378. [Google Scholar]
  58. Yas, M.H.; Kamarian, S.; Pourasghar, A. Application of imperialist competitive algorithm and neural networks to optimise the volume fraction of three-parameter functionally graded beams. J. Exper. Theoret. Art. Intell. 2014, 26, 1–12. [Google Scholar] [CrossRef]
  59. Hagan, M.T.; Demuth, H.B.; Beale, M. Neural Network Design; PWS Publishing Company: Boston, MA, USA, 1995. [Google Scholar]
  60. Bert, C.W.; Malik, M. Differential quadrature: A powerful new technique for analysis of composite structures. Compos. Struct. 1997, 39, 179–189. [Google Scholar] [CrossRef]
  61. Du, H.; Lim, M.K.; Lin, R.M. Application of generalized differential quadrature method to structural problems. Int. J. Numer. Methods Eng. 1994, 37, 1881–1896. [Google Scholar] [CrossRef]
  62. Wu, C.P.; Lee, C.Y. Differential quadrature solution for the free vibration analysis of laminated conical shells with variable stiffness. Int. J. Mech. Sci. 2001, 43, 1853–1869. [Google Scholar] [CrossRef]
  63. Goupee, A.J.; Vel, S. Optimization of natural frequencies of bidirectional functionally graded beams. Struct. Multidiscipl. Opt. 2006, 32, 473–484. [Google Scholar] [CrossRef]
  64. Na, K.S.; Kim, J.H. Optimization of volume fractions for functionally graded panels considering stress and critical temperature. Compos. Struct. 2009, 89, 509–516. [Google Scholar] [CrossRef]
  65. Na, K.S.; Kim, J.H. Volume fraction optimization of functionally graded composite panels for stress reduction and critical temperature. Fin. Elem. Anal. Des. 2009, 45, 845–851. [Google Scholar] [CrossRef]
  66. Na, K.S.; Kim, J.H. Volume fraction optimization for step-formed functionally graded plates considering stress and critical temperature. Compos. Struct. 2010, 92, 1283–1290. [Google Scholar] [CrossRef]
  67. Walker, M.; Smith, R.E. A technique for the multiobjective optimization of laminated composite structures using genetic algorithms and finite element analysis. Compos. Struct. 2003, 62, 123–128. [Google Scholar] [CrossRef]
  68. Tornabene, F.; Ceruti, A. Mixed static and dynamic optimization of four-parameter functionally graded completely doubly curved and degenerate shells and panels using GDQ method. Math. Problems Eng. 2013, 2013, 867079. [Google Scholar] [CrossRef] [Green Version]
  69. Deb, K. Multi-Objective Optimization Using Evolutionary Algorithms; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2002. [Google Scholar]
  70. Konak, A.; Coit, D.W.; Smith, A.E. Multi-objective optimization using genetic algorithms: A tutorial. Reliab. Eng. Syst. Saf. 2006, 91, 992–1007. [Google Scholar] [CrossRef]
  71. Tornabene, F.; Viola, E. Vibration analysis of spherical structural elements using the GDQ method. Comput. Math. Appl. 2007, 53, 1538–1560. [Google Scholar] [CrossRef] [Green Version]
  72. Tornabene, F.; Viola, E. 2-D solution for free vibrations of parabolic shells using generalized differential quadrature method. Eur. J. Mech. A/Solids 2008, 27, 1001–1025. [Google Scholar] [CrossRef]
  73. Eshelby, J.D. The elastic field outside an ellipsoidal inclusion. Proc. Roy. Soc. A 1959, 252, 561–569. [Google Scholar]
  74. Eshelby, J.D. The determination of the elastic field of an ellipsoidal inclusion and related problems. Proc. Roy. Soc. A 1959, 241, 376–396. [Google Scholar]
  75. Mori, T.; Tanaka, K. Average stress in matrix and average elastic energy of materials with misfitting inclusions. Acta Metall. 1973, 21, 571–574. [Google Scholar] [CrossRef]
  76. Ramirez, F.; Heyliger, P.R. Discrete layer solution to free vibrations of functionally graded magneto-electro-elastic plates. Mech. Adv. Mater. Struct. 2006, 13, 249–266. [Google Scholar] [CrossRef]
  77. Ramirez, F.; Heyliger, P.R.; Pan, E. Static analysis of functionally graded elastic anisotropic plates using a discrete layer approach. Compos. Part B 2006, 37, 10–20. [Google Scholar] [CrossRef]
  78. Deb, K.; Goyal, M. A combined genetic adaptive search (GeneAS) for engineering design. J. Comput. Sci. Inform. 1996, 26, 30–45. [Google Scholar]
  79. Khdeir, A.A. Thermal buckling of cross-ply laminated composite beams. Acta Mech. 2001, 149, 201–213. [Google Scholar] [CrossRef]
  80. Li, Z.M.; Qiao, P. Thermal postbuckling analysis of anisotropic laminated beams with different boundary conditions resting on two-parameter elastic foundation. Eur. J. Mech. A/Solids 2015, 54, 30–43. [Google Scholar] [CrossRef]
Figure 1. (a) Configuration and coordinates of a simply-supported functionally graded (FG) beam; (b) cross section of an Nl-layered FG beam.
Figure 1. (a) Configuration and coordinates of a simply-supported functionally graded (FG) beam; (b) cross section of an Nl-layered FG beam.
Jcs 05 00092 g001
Figure 2. Through-thickness distributions of the volume fractions of the ceramic material for a three-parameter power-law material model when (a) κ b = 0.5 , κ p = 5 , and κ c = 1 , 1.5 , 2 , and 3 ; (b) κ b = 0.5 ,   κ c = 1.5 , and κ p = 1 , 5 , 10 , and 20 ; and (c) κ p = 5 , κ c = 1.5 , and κ b = 0 , 0.025 , 0.5 , and 1 .
Figure 2. Through-thickness distributions of the volume fractions of the ceramic material for a three-parameter power-law material model when (a) κ b = 0.5 , κ p = 5 , and κ c = 1 , 1.5 , 2 , and 3 ; (b) κ b = 0.5 ,   κ c = 1.5 , and κ p = 1 , 5 , 10 , and 20 ; and (c) κ p = 5 , κ c = 1.5 , and κ b = 0 , 0.025 , 0.5 , and 1 .
Jcs 05 00092 g002aJcs 05 00092 g002b
Figure 3. Through-thickness distributions of volume fractions of the ceramic material for a one-parameter power-law material model when κ p = 0.2 , 0.5 , 1 , 2 , and 5 .
Figure 3. Through-thickness distributions of volume fractions of the ceramic material for a one-parameter power-law material model when κ p = 0.2 , 0.5 , 1 , 2 , and 5 .
Jcs 05 00092 g003
Figure 4. Flow chart of the non-dominated sorting-based genetic algorithm (GA).
Figure 4. Flow chart of the non-dominated sorting-based genetic algorithm (GA).
Jcs 05 00092 g004
Figure 5. Flow chart of the golden section method used for the cases of temperature-dependent (TD) material properties.
Figure 5. Flow chart of the golden section method used for the cases of temperature-dependent (TD) material properties.
Jcs 05 00092 g005
Figure 6. Variations of material properties of the ceramic and the metal materials with the temperature: (a) Young’s modulus; (b) Poisson’s ratio; (c) Thermal expansion coefficient.
Figure 6. Variations of material properties of the ceramic and the metal materials with the temperature: (a) Young’s modulus; (b) Poisson’s ratio; (c) Thermal expansion coefficient.
Jcs 05 00092 g006
Figure 7. Two hundred populations for different generations shown in the objective spaces, in which the temperature-independent (TI) material properties are considered; (a) initial generation; (b) 1st generation; (c) 2nd generation; (d) 5th generation; (e) 10th generation; (f) 20th generation.
Figure 7. Two hundred populations for different generations shown in the objective spaces, in which the temperature-independent (TI) material properties are considered; (a) initial generation; (b) 1st generation; (c) 2nd generation; (d) 5th generation; (e) 10th generation; (f) 20th generation.
Jcs 05 00092 g007
Figure 8. Two hundred populations for different generations shown in the objective spaces, in which the TD material properties are considered; (a) initial generation; (b) 1st generation; (c) 2nd generation; (d) 5th generation; (e) 10th generation; (f) 20th generation.
Figure 8. Two hundred populations for different generations shown in the objective spaces, in which the TD material properties are considered; (a) initial generation; (b) 1st generation; (c) 2nd generation; (d) 5th generation; (e) 10th generation; (f) 20th generation.
Jcs 05 00092 g008
Figure 9. Curve fitting for the Pareto-optimal fronts for the cases of TI and TD material properties.
Figure 9. Curve fitting for the Pareto-optimal fronts for the cases of TI and TD material properties.
Jcs 05 00092 g009
Table 1. Convergent study of critical temperature parameter solutions of simply supported, laminated (0°/90°/0°) beams under a uniform temperature change ( α 2 / α 1 = 3 , E 1 / E 2 = 20 ).
Table 1. Convergent study of critical temperature parameter solutions of simply supported, laminated (0°/90°/0°) beams under a uniform temperature change ( α 2 / α 1 = 3 , E 1 / E 2 = 20 ).
L/hTheories
Δ T ¯ c r
10Current LW1 (nl = 1)0.8315
Current LW1 (nl = 3)0.8123
Current LW1 (nl = 6)0.7995
Current LW1 (nl = 9)0.7970
Current LW2 (nl = 1)0.8306
Current LW2 (nl = 3)0.7953
Current LW2 (nl = 6)0.7950
Current LW3 (nl = 1)0.8010
Current LW3 (nl = 3)0.7950
Current LW3 (nl = 6)0.7950
a Current LW3 (nl = 6)0.8001
b Current LW3 (nl = 6)0.8383
a, b Current LW3 (nl = 6)0.8438
HSDT [80]0.78912
HRSDT [22]0.8230
TSDT [79]0.8229
FSDT [79]0.8281
EBT [79]1.1072
50Current LW1 (nl = 1)1.0411
Current LW1 (nl = 3)1.0370
Current LW1 (nl = 6)1.0359
Current LW1 (nl = 9)1.0357
Current LW2 (nl = 1)1.0392
Current LW2 (nl = 3)1.0355
Current LW2 (nl = 6)1.0355
Current LW3 (nl = 1)1.0373
Current LW3 (nl = 3)1.0355
Current LW3 (nl = 6)1.0355
a Current LW3 (nl = 6)1.0360
b Current LW3 (nl = 6)1.0931
a, b Current LW3 (nl = 6)1.0936
HSDT [80]1.04656
HRSDT [22]1.0921
TSDT [79]1.0921
FSDT [79]1.0925
EBT [79]1.1072
a von Karman geometric nonlinear terms are considered. b The Poisson ratio is considered to be zero in Qij (i = 1 and 3).
Table 2. Convergent study of critical temperature parameter solutions of simply supported, laminated (0°/90°/0°) beams under a uniform temperature change ( α 2 / α 1 = 3 , L / h = 10 ).
Table 2. Convergent study of critical temperature parameter solutions of simply supported, laminated (0°/90°/0°) beams under a uniform temperature change ( α 2 / α 1 = 3 , L / h = 10 ).
E 1 / E 2
Theories
Δ T ¯ c r
10Current LW1 (nl = 1)0.8335
Current LW1 (nl = 3)0.8162
Current LW1 (nl = 6)0.8088
Current LW1 (nl = 9)0.8073
Current LW2 (nl = 1)0.8302
Current LW2 (nl = 3)0.8063
Current LW2 (nl = 6)0.8062
Current LW3 (nl = 1)0.8134
Current LW3 (nl = 3)0.8061
Current LW3 (nl = 6)0.8061
a Current LW3 (nl = 6)0.8123
b Current LW3 (nl = 6)0.8896
a, b Current LW3 (nl = 6)0.8965
HSDT [80]0.81683
HRSDT [22]0.8833
TSDT [79]0.8832
FSDT [79]0.8868
EBT [79]1.0370
30Current LW1 (nl = 1)0.7828
Current LW1 (nl = 3)0.7622
Current LW1 (nl = 6)0.7463
Current LW1 (nl = 9)0.7432
Current LW2 (nl = 1)0.7824
Current LW2 (nl = 3)0.7411
Current LW2 (nl = 6)0.7407
Current LW3 (nl = 1)0.7454
Current LW3 (nl = 3)0.7406
Current LW3 (nl = 6)0.7406
a Current LW3 (nl = 6)0.7445
b Current LW3 (nl = 6)0.7680
a, b Current LW3 (nl = 6)0.7720
HSDT [80]0.72608
HRSDT [22]0.7472
TSDT [79]0.7471
FSDT [79]0.7528
EBT [79]1.1329
a von Karman geometric nonlinear terms are considered. b The Poisson ratio is considered to be zero in Qij (i = 1 and 3).
Table 3. Critical temperature parameter solutions for simply supported, FG beams under a uniform temperature change.
Table 3. Critical temperature parameter solutions for simply supported, FG beams under a uniform temperature change.
L/h κ p Micromechanical Models Theories Δ T ¯ c r ( TI ) Δ T ¯ c r ( TD )
101Mori–TanakaCurrent LW1 (nl = 2)0.51510.3757
Current LW1 (nl = 4)0.50430.3692
Current LW1 (nl = 8)0.50150.3674
Current LW2 (nl = 2)0.50070.3674
Current LW2 (nl = 4)0.50060.3667
Current LW3 (nl = 2)0.50060.3667
Current LW3 (nl = 4)0.50060.3667
101Rule of mixturesCurrent LW1 (nl = 2)0.51220.3678
Current LW1 (nl = 4)0.50140.3614
Current LW1 (nl = 8)0.49870.3600
Current LW2 (nl = 2)0.49780.3595
Current LW2 (nl = 4)0.49780.3595
Current LW3 (nl = 2)0.49770.3595
Current LW3 (nl = 4)0.49770.3595
105Mori–TanakaCurrent LW1 (nl = 2)0.53150.4525
Current LW1 (nl = 4)0.52020.4436
Current LW1 (nl = 8)0.51720.4412
Current LW2 (nl = 2)0.51640.4405
Current LW2 (nl = 4)0.51630.4405
Current LW3 (nl = 2)0.51630.4405
Current LW3 (nl = 4)0.51630.4405
105Rule of mixturesCurrent LW1 (nl = 2)0.53070.4467
Current LW1 (nl = 4)0.51930.4381
Current LW1 (nl = 8)0.51640.4358
Current LW2 (nl = 2)0.51550.4350
Current LW2 (nl = 4)0.51540.4350
Current LW3 (nl = 2)0.51540.4350
Current LW3 (nl = 4)0.51540.4350
205Mori–TanakaCurrent LW1 (nl = 2)0.54320.5194
Current LW1 (nl = 4)0.53220.5081
Current LW1 (nl = 8)0.52950.5066
Current LW2 (nl = 2)0.52860.5042
Current LW2 (nl = 4)0.52850.5042
Current LW3 (nl = 2)0.52850.5042
Current LW3 (nl = 4)0.52850.5042
205Rule of mixturesCurrent LW1 (nl = 2)0.54230.5170
Current LW1 (nl = 4)0.53140.5066
Current LW1 (nl = 8)0.52860.5042
Current LW2 (nl = 2)0.52770.5026
Current LW2 (nl = 4)0.52770.5026
Current LW3 (nl = 2)0.52770.5026
Current LW3 (nl = 4)0.52770.5026
Table 4. Temperature dependent material properties of the metal material (SUS304) and the ceramic material (ZrO2) [4].
Table 4. Temperature dependent material properties of the metal material (SUS304) and the ceramic material (ZrO2) [4].
MaterialsP0P−1P1P2P3P at 300K
ZrO2
E244.27 × 1090−1.371 × 10−31.214 × 10−6−3.681 × 10−10168.06 × 109
α 12.766 × 10−60−1.491 × 10−31.006 × 10−5−6.778 × 10−1118.591 × 10−6
υ 0.288201.133 × 10−4000.298
ρ 365700003657
SUS304
E201.04 × 10903.079 ×10−4−6.534 × 10−70207.79 × 109
α 12.330 × 10−608.086 × 10−40015.321 × 10−6
υ 0.32620−2.002 × 10−43.797 × 10−700.318
ρ 816600008166
Table 5. Pareto-optimal solutions for simply supported FG beams under a uniform temperature change sorted according to the F2 function values from the largest to the smallest, in which temperature-independent material properties are considered.
Table 5. Pareto-optimal solutions for simply supported FG beams under a uniform temperature change sorted according to the F2 function values from the largest to the smallest, in which temperature-independent material properties are considered.
Order Number of Pareto-
Optimal Solutions
(F1, F2) ρ ^ Δ T c r ( κ b , κ c , κ p ) ( w 1 , w 2 )
1(0.000, 1.000)3658.1300.7326(0.358, 1.829, 0.000)(0.565, 0.435)
10(0.037, 0.953)3820.7303.1011(0.000, 1.000, 0.037)(0.558, 0.442)
20(0.072, 0.910)3975.6305.2840(0.000, 2.392, 0.076)(0.549, 0.451)
30(0.113, 0.861)4157.9307.7473(0.000, 3.000, 0.125)(0.537, 0.463)
40(0.154, 0.815)4337.3310.0472(0.000, 2.831, 0.178)(0.523, 0.477)
50(0.206, 0.760)4565.6312.7782(0.000, 2.759, 0.253)(0.502, 0.498)
60(0.251, 0.716)4768.5315.0162(0.000, 2.714, 0.328)(0.482, 0.518)
70(0.309, 0.665)5023.2317.5743(0.000, 3.000, 0.435)(0.456, 0.544)
80(0.380, 0.610)5335.0320.3524(0.000, 2.967, 0.591)(0.424, 0.576)
90(0.433, 0.572)5570.0322.2288(0.002, 2.958, 0.736)(0.403, 0.597)
100(0.496, 0.531)5851.2324.2938(0.001, 2.958, 0.945)(0.385, 0.615)
110(0.576, 0.482)6204.6326.7625(0.004, 2.820, 1.300)(0.381, 0.619)
120(0.644, 0.439)6502.7328.9122(0.000, 3.000, 1.723)(0.399, 0.601)
130(0.695, 0.403)6731.2330.7254(0.000, 2.745, 2.190)(0.428, 0.572)
140(0.741, 0.367)6934.5332.5478(0.653, 1.122, 7.805)(0.459, 0.541)
150(0.801, 0.311)7199.1335.3457(0.864, 1.065, 22.95)(0.506, 0.494)
160(0.844, 0.264)7387.2337.7239(0.864, 1.082, 25.69)(0.545, 0.455)
170(0.883, 0.213)7559.9340.2622(0.881, 1.179, 24.36)(0.580, 0.420)
180(0.926, 0.148)7751.7343.5581(0.954, 1.146, 50.00)(0.621, 0.379)
190(0.967, 0.076)7930.0347.1474(0.537, 1.000, 45.04)(0.661, 0.339)
200(1.000, 0.008)8077.5350.5544(0.000, 2.005, 50.00)(0.684, 0.316)
Table 6. Pareto-optimal solutions of simply supported FG beams under a uniform temperature change sorted according to the F2 function values from the largest to the smallest, in which temperature-dependent material properties are considered.
Table 6. Pareto-optimal solutions of simply supported FG beams under a uniform temperature change sorted according to the F2 function values from the largest to the smallest, in which temperature-dependent material properties are considered.
Order Number of Pareto-
Optimal Solutions
(F1, F2) ρ ^ Δ T c r ( κ b , κ c , κ p ) ( w 1 , w 2 )
1(0.000, 0.994)3657.0172.5(0.465, 1.620, 0.000)(0.272, 0.728)
10(0.040, 0.977)3833.8174.4(0.000, 3.000, 0.041)(0.315, 0.685)
20(0.098, 0.947)4090.2177.7(0.000, 3.000, 0.105)(0.354, 0.646)
30(0.151, 0.917)4326.6181.1(0.000, 1.799, 0.173)(0.375, 0.625)
40(0.199, 0.888)4535.2184.4(0.000, 1.787, 0.241)(0.385, 0.615)
50(0.246, 0.858)4743.0187.7(0.000, 3.000, 0.318)(0.390, 0.610)
60(0.311, 0.816)5031.9192.4(0.000, 1.000, 0.440)(0.391, 0.609)
70(0.364, 0.781)5268.3196.2(0.000, 2.613, 0.560)(0.391, 0.609)
80(0.424, 0.743)5532.8200.6(0.000, 1.561, 0.716)(0.393, 0.607)
90(0.483, 0.705)5791.1204.8(0.000, 1.474, 0.901)(0.400, 0.600)
100(0.531, 0.672)6002.2208.5(0.000, 2.792, 1.086)(0.411, 0.589)
110(0.585, 0.633)6242.9212.9(0.013, 1.444, 1.364)(0.431, 0.569)
120(0.651, 0.579)6533.8218.8(0.000, 2.239, 1.746)(0.465, 0.535)
130(0.697, 0.537)6736.6223.6(0.000, 1.443, 2.132)(0.496, 0.504)
140(0.746, 0.484)6956.2229.5(0.046, 2.894, 2.754)(0.535, 0.465)
150(0.793, 0.425)7164.5236.1(0.343, 1.047, 5.655)(0.575, 0.425)
160(0.836, 0.362)7353.8243.1(0.373, 1.592, 5.894)(0.612, 0.388)
170(0.884, 0.281)7562.9252.1(0.428, 1.210, 10.03)(0.652, 0.348)
180(0.915, 0.219)7703.5259.1(0.798, 1.231, 21.23)(0.677, 0.323)
190(0.965, 0.105)7924.0271.9(0.879, 1.297, 34.51)(0.713, 0.287)
200(1.000, 0.013)8077.6282.1(0.117, 1.000, 49.74)(0.736, 0.264)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wu, C.-P.; Li, K.-W. Multi-Objective Optimization of Functionally Graded Beams Using a Genetic Algorithm with Non-Dominated Sorting. J. Compos. Sci. 2021, 5, 92. https://doi.org/10.3390/jcs5040092

AMA Style

Wu C-P, Li K-W. Multi-Objective Optimization of Functionally Graded Beams Using a Genetic Algorithm with Non-Dominated Sorting. Journal of Composites Science. 2021; 5(4):92. https://doi.org/10.3390/jcs5040092

Chicago/Turabian Style

Wu, Chih-Ping, and Kuan-Wei Li. 2021. "Multi-Objective Optimization of Functionally Graded Beams Using a Genetic Algorithm with Non-Dominated Sorting" Journal of Composites Science 5, no. 4: 92. https://doi.org/10.3390/jcs5040092

Article Metrics

Back to TopTop