Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Mechanisms of Component Degradation and Multi-Scale Strategies for Predicting Composite Durability: Present and Future Perspectives
Next Article in Special Issue
Physical, Mechanical and Microstructural Characteristics of Perlite-Based Geopolymers Modified with Mineral Additives
Previous Article in Journal
Enhancing Hydrogels with Quantum Dots
Previous Article in Special Issue
Impact Performance of 3D Orthogonal Woven Composites: A Finite Element Study on Structural Parameters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metaheuristic Optimization of Functionally Graded 2D and 3D Discrete Structures Using the Red Fox Algorithm

1
CIMOSM—Centro de Investigação em Modelação e Otimização de Sistemas Multifuncionais, ISEL, IPL—Instituto Politécnico de Lisboa, Av. Conselheiro Emídio Navarro 1, 1959-007 Lisboa, Portugal
2
IDMEC, Instituto Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais, 1, 1049-001 Lisboa, Portugal
*
Authors to whom correspondence should be addressed.
J. Compos. Sci. 2024, 8(6), 205; https://doi.org/10.3390/jcs8060205
Submission received: 5 April 2024 / Revised: 18 May 2024 / Accepted: 22 May 2024 / Published: 30 May 2024
(This article belongs to the Special Issue Characterization and Modelling of Composites, Volume III)

Abstract

:
The growing applicability of functionally graded materials is justified by their ability to contribute to the development of advanced solutions characterized by the material customization, through the selection of the best parameters that will confer the best mechanical behaviour for a given structure under specific operating conditions. The present work aims to attain the optimal design solutions for a set of illustrative 2D and 3D discrete structures built from functionally graded materials using the Red Fox Optimization Algorithm, where the design variables are material parameters. From the results achieved one concludes that the optimal selection and distribution of the different materials’ mixture and the different exponents associated with the volume fraction law significantly influence the optimal responses found. To note additionally the good performance of the coupling between this optimization technique and the finite element method used for the linear static and free vibration analyses.

1. Introduction

Functionally graded materials (FGMs) are usually particle composite materials, in which their constituent material volume fractions vary continuously according to given spatial coordinates [1,2]. This material concept emerged in Japan in 1984 for application as a thermal barrier material, capable of withstanding high-temperature gradients [2,3,4]. In a more general perspective, these composites’ properties and characteristics allow them to be used with increasing frequency in the transportation industry, namely in aircraft and spacecraft structures, as well as in the automotive and marine industries. Other applications include electronic and medical pieces of equipment and in a quite transversal manner their use as thermal coatings for engines and turbines, among others. FGM are in general made from different constituent materials, namely ceramic and metallic material phases which have shown to be appropriate to operate as a thermal barrier [5,6,7]. The composition of these materials often varies from one surface rich in material A to another surface rich in material B, with a certain distribution between the two surfaces, resulting in a continuous variation in the material microstructure. This continuous variation allows the effective properties of FGMs to have a gradual and smooth evolution, which makes it possible to minimize undesirable situations that occur in fibre-reinforced laminated structures due to their anisotropy, namely stress concentrations in the vicinity of material and geometric discontinuities and residual stresses, which can lead to delamination, cracks in the matrix and separation of adhesive bonds [8,9,10,11,12].
However, depending on the applications’ characteristics, other material mixtures have been utilized as illustrated in the studies developed by [13,14,15,16,17], where mixtures of different metallic, polymeric, and ceramic materials were used. The FGM concept extended to materials with specific biocompatible characterisics can also be found with an increasing frequency, namely in [18,19,20], among others.
Considering the wide range of potentially variable parameters involved in a generic structure design, obtaining better solutions according to certain objectives for a given structure subjected to specific loads and boundary conditions may be a complex task. The increasing complexity and non-linearity of real problems require the development of efficient optimization methods with high accuracy and response speed combined with low computational cost [21,22]. Metaheuristics are optimization techniques that are often based on natural phenomena [22] and have shown to be adequate to address such complex problems. They are generally based on rules translated from the observation of what occurs in nature considering randomness, to mimic some of the characteristics of biological behaviour, and natural or physical phenomena [23].
When one talks about FGM optimization [24], the parameters normally considered as design variables characterize, on one hand, the evolution of the volume fractions of the material (the power law exponent or other laws) and, on the other hand, the geometric characteristics of the cross-section. The objective functions usually aim at maximizing or minimizing mass, maximum stresses, resistance to fracture, behaviour in situations with temperature gradients, and dynamic response in free and forced vibration, or buckling loads [2,25]. To illustrate this, one can refer to the work conducted in [5] where a model was developed for the longitudinal variation of the composition of a metal/ceramic FGM, to optimize the heat flow through this material as a function of the material composition profile. Another work carried out in [2], the influence of the dimension for which the volume fraction is defined, in a thin cantilever beam with a hollow profile was studied, being found that the maximization of the natural frequencies, was benefited by a material mixture distribution in the longitudinal direction rather than considering it in the thickness direction, which has been more commonly studied. A very recent work on topology optimization of micro and macroscale topology optimization of FGM lattice structures was presented by [26]. Other studies focusing on optimizing the mechanical behaviour of composite structures have been studied by several researchers [27,28] considering a wide variety of optimization techniques.
The optimization algorithm considered in the present work consists of a very recent technique which has been used so far, mainly in the context of medical science applications, for example for lung x-ray image segmentation [29], electroencephalogram signals classification [30], in computer science for example, to optimize hyperparameters of deep-learning models [31], and very scarcely in engineering for example for the optimization of phase equilibrium and stability of chemical systems [32] or to estimate optimal model parameters of solid oxide fuel cells [33].
However, no published work was found considering the use of this technique in the optimization of the mechanical behaviour of structures in general, nor to optimize the innovative height functionally graded structures considered in the present work, in particular.
This work is a natural sequence of the study conducted in [34], which studied the influence of the power law exponent and the effective properties estimation, on the static and dynamic behaviour of two-dimensional planar structures in the free vibration regime. By generalizing and combining the previously developed structural analysis application with the present metaheuristic optimization technique, it becomes possible to optimizing the desired structure behaviour without extensive parametric studies and with less probability of achieving the effective best configuration.
The selection of the present optimization metaheuristic is mainly related to the conclusion drawn in [35] which stated that this technique shows better exploration and exploitation characteristics when compared with other techniques also analysed, namely the Genetic Algorithm, the Particle Swarm Optimization, the Grey Wolf Optimization, the Chimp Optimization Algorithm, the Butterfly Optimization Algorithm, the Whale Optimization Algorithm, the Fitness Dependent Optimizer and the Dragonfly Optimization Algorithm.
Hence, with this work, one aims to obtain the best performance in terms of maximum linear static displacements and fundamental frequencies for a set of two and three-dimensional structures in which the material gradient is a function of the vertical coordinate of the structure. From the authors’ knowledge, no published work focused on optimizing these types of FGM structures, or even on using the Red Fox algorithm for structural optimization, was found.
The remainder part of this work is organized as follows: Section 2 describes the concept of the FGMs that will be used in the present study, Section 3 presents the displacement field, constitutive relations and equilibrium equations associated with the linear static and free vibrations analyses to be developed for each individual considered in the context of the populations-based optimization technique that will be presented in Section 4. The methodology is presented in Section 5 and Section 6 contains the case studies considered. Finally, in Section 7 some conclusions are drawn based on the results achieved.

2. Functionally Graded Materials

This study considers the optimization of FGM structures in which the material gradient is a function of the vertical coordinate of the structure (Y) [34]. The objective functions considered are minimizing the linear static displacements of selected points within the structures or maximizing their fundamental frequencies. Although it is possible to consider other design variables, for example, those of a geometrical nature, this work is focused on obtaining the material parameters’ configurations that will make it possible to achieve better structural responses without changing the structure’s geometry.
Figure 1 illustrates a two-dimensional frame-type structure where the material gradient through the Y direction, associated with the structure height, is represented.
If Y 0 is the ordinate of the lower end of the structure, H corresponds to the maximum height of the structure, Y is the vertical coordinate that varies between 0 and H, and ey is the exponent of the power law that adjusts how the property variation along the vertical coordinate occur, the generic expression of the volume fraction can be written as:
V f = Y Y 0 H e y
When the operating conditions of a given structure are known, it is possible to combine not only the most appropriate geometric characteristics, but also, as this work aims to demonstrate, the best constituent materials from within a given set, and the best parameters that control the distribution of the mixture of constituent materials mixture throughout the structure, to optimizing specific structural requirements [8,9].
The effective properties of FGMs are estimated according to the Voigt Rule of Mixtures [8], which makes it possible to determine the average value of a generic material property P F G M as a weighted average of the corresponding constituent materials’ properties ( P A ,   P B ), where the weights are their volume fractions, as in the expression [36]:
P F G M = P A V f + P B 1 V f
The sum of the volume fractions must be unitary in the case of dual-phase FGMs, considering the assumption that the material has no porosities. The expression that describes a material volume fraction ( V f ) evolution will reflect how that material will be distributed in a specific direction or multiple directions, which for the present work is given by Equation (1).

3. Displacement Field, Constitutive Relations, and Equilibrium Equations

Considering the structures that will be studied, one has to consider the first-order shear deformation to model their kinematical modelling [17,18,19,20,21]. The displacement field corresponding to this theory can be written as follows:
u x , y , z , t = u 0 x , y , t + z θ y 0 x , y , t y θ z 0 x , y , t v x , y , z , t = v 0 x , y , t z θ x 0 x , y , t w x , y , z , t = w 0 x , y , t + y θ x 0 x , y , t
The terms, u 0 , v 0   a n d   w 0 correspond, respectively, to the longitudinal (direction x ) and transverse (directions y and z ) displacements of the beam’s centreline in the x z and x y planes, θ x 0 corresponds to the torsion angle, and θ y 0 and θ z 0 correspond to the rotation of the beam’s centre plane in a direction perpendicular to the x z plane and x y planes, respectively. Considering the Elasticity Theory for small deformations, the corresponding deformation field is written as:
ε x x = u 0 x + z θ y 0 x y θ z 0 x γ x z = w 0 x + y θ x 0 x + θ z 0 γ x y = v 0 x z θ x 0 x θ y 0
The remaining strains ε y y , ε z z and γ y z are null, and the corresponding constitutive relations for a generic FGM are described as [36,37]:
σ x x σ x z σ x y = Q 11 x , y , z 0 0 0 Q 55 x , y , z k s 0 0 0 Q 66 x , y , z k s ε x x γ x z γ x y
where Q 11 = E F G M , and Q 55 = Q 66 = E F G M 2 1 + ν F G M . E F G M represents the FGM Young’s modulus and ν F G M the FGM Poisson’s ratio. The shear correction factor is set to k s = 5 / 6 .
Considering the aim of the present work, one has to use the principle of Hamilton to derive the equilibrium equations, which can be written as:
δ t 1 t 2 T Ω + U d t = 0
with T standing for the kinetic energy, U for the elastic strain energy and Ω for the work done by the forces applied to the system. These quantities can be written as:
T = 1 2 V ρ F G M u ˙ 2 + v ˙ 2 + w ˙ 2 d V U = 1 2 V σ x x ε x x + σ x z γ x z + σ x y γ x y d V Ω = k = 1 N f k q k
where ρ F G M represents the FGM density, f k stands for the punctual generalized forces and q k the generalized degrees of freedom. From the minimization of the Lagrangean functional, one achieves the equilibrium equations at the element level:
M e q ¨ e + K e q e = F e
which, for linear static analysis or for free vibration analysis will be written as:
K e q e = F e K e ω 2 M e q e = 0
These equilibrium equations were implemented using quadratic beam elements with six degrees of freedom per node—the axial ( u k ) and the transverse ( v k , w k ) displacements and the rotations ( θ x ,   θ y ,   θ z ), all to the midplane surface of the beam element. Figure 2 schematically shows the beam element in the three-dimensional space and the degrees of freedom associated with each node [38,39,40,41,42].
Further details on the implementation of the beam-bar element used to model these structures can be found in [34].

4. Structural Optimization

A generic optimization problem can be stated as:
min f ( x ^ ) s u b j e c t   t o g i x ^ 0 , i = 1 , , N h j x ^ = 0 , j = 1 , , M
where f ( x ^ ) is the objective function and x ^ is the vector of the design variables. Besides the side constraints associated with the design variables, the inequality and equality behavioral constraints are, respectively, designated as g i x ^ and h j x ^ .
In the present paper, the optimizations that are intended to be carried out are related to the minimization of the static displacements of specific points for a set of structures or the maximization of their fundamental frequencies. To this purpose one considers a metaheuristic optimization technique, the red fox optimization (RFO) algorithm.

Red Fox Metaheuristic Algorithm

The red fox is a predator that hunts not only wild animals, but also farm animals, and in turn can be hunted by humans who intercept them in these attacks. The relationships between prey and predators and even between individuals in the same pack can be mimicked mathematically with a view to optimization, since only the fittest individuals survive to hunt in the different territories and to form a new pack elsewhere. The proposed model by [22] assumes that the optimization domain is comparable to hunting regions. Foxes have a hunting mechanism that can be considered in a first global search and then a local search. In addition to these two phases, a population control model has also been proposed, simulating the group leaving for another location or being killed by hunters. For this, the two best individuals are selected to reproduce and replace the weakest elements in each iteration. This mechanism allows the model to cover the entire domain to reach the global extreme, thus minimizing the risk of being confined to local extremes [22]. The algorithm is structured in a logical set of five steps as follows:
Step 1—Initial Population: The population is created with a selected number of foxes, which will be constant in each iteration. Each element is represented by x ¯ i j t , where i corresponds to the identification of the individual in the population, j represents the coordinate according to the dimension of the results space, and t is the current iteration. Based on the evaluation of the population generated through the objective function, the best individuals that provide the best solution are recorded, originating the population of the first iteration.
Step 2—Global Search: The members of the pack move to remote locations searching for food, and the information they collect is later shared with the other members to ensure their survival and development. The exploration of the surrounding territory is modelled on the basis of the objective function, and it is assumed that the best element explores the most interesting territories and shares this information with the group. Thus, the best individual in the population is recorded and the Euclidean distance of each of the other foxes is calculated using:
d x i ¯ t , x ¯ b e s t t = | | ( x ¯ i t x ¯ b e s t t | |
Individuals will move towards the fittest red fox through:
x ¯ i t = x ¯ i t + α s i g n x ¯ b e s t t x ¯ i t
where α is a randomly generated parameter for all the individuals in each iteration, belonging to the interval [0, d x ¯ i   t , x ¯ b e s t   t ]. At this point, the new population is evaluated using the objective function, and if the optimal values are better, the individuals take new positions, otherwise they return to the starting position. This simulates the path that the remaining members of the group will have to follow, after the information has been shared by the elements that have explored the surrounding territories, in which they must continue searching in promising areas, corresponding to step 3. The model also predicts the death of the worst individuals in the population, or the reproduction of the best individuals according to the equations in step 4.
Step 3—Local Search: When foxes spot prey, they try to approach it slowly in a circular motion so as not to be seen, until they are close enough to catch it. This behaviour is modelled using the parameter μ. This parameter is randomly generated in the interval [0, 1] and simulates the possibility of the fox being seen by the prey during the approach movement. For μ > 0.75 , the fox approaches the prey, otherwise it moves away. The movement of observation and approach of each individual to the prey was approximated using a modified Cochleoid equation, through the parameters a and ϕ 0 , where a is randomly generated for each iteration in the interval [0, 0.2], and characterizes the change in distance to the prey, and ϕ 0 models the angle of vision of each individual and is randomly generated in the interval [0, 2π] for all individuals. The radius of vision ( r v ) of the foxes is obtained by:
r v = a   sin ϕ 0 ϕ 0 ,     ϕ 0 0 θ ,     ϕ 0 = 0
where θ is a random value between [0, 1] generated at the start of the algorithm and characterises the influence of atmospheric conditions. For the individuals in the population, their movement is modelled using the system of Equation (14) for spatial coordinates, where each angular value ϕ n 1 is randomly generated in the interval [0, 2π].
x n e w 0 = a   r cos ϕ 1 + x a c t u a l 0 x n e w 1 = a   r sin ϕ 1 + a   r cos ϕ 2 + x a c t u a l 1 x n e w 2 = a   r sin ϕ 1 + a r sin ϕ 2 + a   r cos ϕ 3 + x a c t u a l 2 x n e w n 2 = a   r   k = 1 n 2 sin ϕ k + a   r cos ϕ n 1 + x a c t u a l n 2 x n e w n 1 = a   r sin ϕ 1 + a r sin ϕ 2 + + a   r sin ϕ n 1 + x a c t u a l n 1
Step 4—Reproduction and Pack Abandonment: For the dynamic control of the population, five per cent of the worst individuals in the population are selected according to the objective function, to introduce small changes to the pack. As these are the worst individuals, it is assumed that they have either moved territory or been killed by hunters; however, to ensure that the population is constant, the worst individuals are replaced by new individuals. In each iteration, the two best individuals (the alpha couple) are selected, x ¯ ( 1 ) t and x ¯ ( 2 ) t , and the centre of their habitat and diameter are calculated:
h a b i t a t c e n t r e t = x ¯ 1 t + x ¯ 2 t 2 h a b i t a t d i a m e t e r t = ( x ¯ ( 1 ) ) t ( x ¯ ( 2 ) ) t
For each iteration, a random parameter κ is generated within the interval [0, 1], which defines the replacements of individuals. For κ 0.45 , new individuals are introduced into the population, otherwise the alpha couple reproduces. In the first case, it is assumed that new family members will move to other regions, being able to reproduce. In the second case, new individuals are the result of the reproduction of the alpha couple, which results in:
x ¯ r e p r o d u c e d t = κ x ¯ 1 t + x ¯ 2 t 2
Step 5—Evaluating each individual according to the objective function—From the objective function, the optimal values x b e s t are recorded. For further details on the algorithm, please consult the corresponding literature [22].

5. Methodology

The generic goal of the optimization process is to determine the design parameter configuration that provide the best structure response, concerning a specific objective function to study.
In the present work, the studies developed have been structured into two types: the first considering only one design variable and the other considering four design variables. Both are, related to the material constitution of the FGM structures. In the first stage, the optimization process is related to determining the exponent of the power law that corresponds to the best behaviour of the structure for a given pair of materials. This design variable was assumed to vary in a continuous manner in the interval [0.1, 10].
The procedure is schematically represented in Figure 3 and highlights the construction of the finite element model for the specific data associated with each of the individuals of the population, within each optimization iteration.
As observed in this case, the FGM structure is made of aluminum and zirconia, and the design variable, the power law exponent, varies within the mentioned interval. The rule of mixtures considers that the power law’s volume fraction is directly associated with the distribution of metallic material along the height of the structure.
In the second stage, one considers a wider range of possibilities concerning the nature and number of design variables. Hence, in this case, it is possible to have different pairs of constituent materials to the resulting FGM, and to write the rule of mixtures associating the volume fraction directly to one or another of those constituent materials. A schematic diagram representation for the four design variable model is presented in Figure 4.
So, in addition to the power law exponent, another three variables were added, which can only take integer values of 1 or 2. The first one corresponds to the writing of the rule of mixtures, i.e., if the random number assumes 1, the volume fraction characterizes the distribution of metallic material; if it assumes 2, then the volume fraction describes the ceramic material distribution. The second additional design variable defines the metallic material to be used, in this case, aluminium and titanium, also randomly selected. The last additional design variable defines the possible ceramic material to use; zirconia and tungsten carbide were considered as options in the present study. For these three additional design variables introduced, that can only take the values {1, 2}, the expression that characterizes the update of the individual position i in each iteration of each algorithm is adjusted by rounding it to the nearest integer value.
The optimization procedure is illustrated in a summarized manner, via the flowchart depicted in Figure 5.
The present work thus considers the coupling between the red fox’s optimization technique and the finite element method that enables carrying out the necessary linear static and free vibration analyses for each individual, within each population.
In this work, which constitutes a first study in this context, a single-objective approach is considered. Also, besides the side constraints associated with each of the design variables space, possible behavioural constraints are not considered.

6. Results and Discussion

This section presents two verification cases, considering benchmark functions to illustrate the performance of the optimization technique. These cases are then followed by two additional performance assessment studies on FGM structures, where the results achieved via optimization processes, are compared with the trends obtained via parametric studies.
After, one conducts a set of case studies, to minimize specific points’ displacements in a certain structure or maximise their fundamental frequencies. The material properties used in the optimization cases next presented are the ones in Figure 3 and Figure 4.

6.1. Verification Cases

To verify the red fox optimization technique implemented, two studies are first presented using multimodal benchmark functions usually considered for this kind of verification purposes: the “Sum of Squares” function and the “Rastringin” function.
  • Case 1—“Sum of Squares” Benchmark function
The optimal values obtained through the implementation of this algorithm in the present work were checked against those obtained by [22], where each function was evaluated in 100 independent runs, 100 iterations and for 100 individuals.
The first 20 independent runs out of 100 (presented in Appendix A) were carried out for the “Sum of Squares” function, using 20 design variables, considering 100 iterations and 100 individuals. The results obtained and their characteristics are shown in Table 1.
The results obtained are shown in Table 2.
According to the results here presented and, in the Appendix A, one concludes that the behaviour of the optimization algorithm is very satisfactory, with average values for the global minima of the function in the different runs, close to the found in the literature and very close to the exact solution. The standard deviation shows a similar magnitude.
  • Case 2—“Rastringin ”Benchmark function
For the second verification case study, related to the “Rastringin” function, 10 independent runs were carried out, using two design variables, considering 100 iterations and 100 individuals. The characteristics of this function and the analysed domain are shown in Table 3.
As can be concluded from the results shown in Table 4, the implemented optimization technique presents a very good performance, with an average value very close to the exact solution. Again, the standard deviation shows a similar magnitude.

6.2. Performance Assessment Studies

In the previous sub-section, one considered two different benchmark functions, and compared the results obtained with the ones by other authors. In the present sub-section one intends to complement the previous one, applying now the optimization technique to two FGM structures. To be able to conclude on the performance of the technique for this type of problem, one has first conducted parametric studies for each structure to show the influence of the exponent ey, followed by optimization processes where that exponent is the design variable. The objective of this procedure is to demonstrate that the optimization results respect the known trends shown in the parametric studies.
In both cases, the transverse cross-section of the elements is square, with a 20 mm edge dimension, and an aspect ratio of 30. Considering the good performance of the optimization technique, in previous verification studies, a series of 10 runs was conducted for each optimization process, using as optimization parameters, a population of 5 individuals and 10 iterations.
  • Case 3—2D Truss—Analysis and Optimization
The first structure studied is a truss with 21 elements, subjected to a load P = 1 kN, as shown in Figure 6.
Figure 7 presents the minimum discretization of the analysed 2D truss, with quadratic beam elements.
For this truss, the parametric studies conducted allow to conclude that the node where the largest resultant displacement occurs is at node 19, which corresponds to the node where the load is applied. Table 5 shows the resultant displacements of the mentioned node, and the fundamental frequency of the structure, for a set of exponents of the volume fraction power law.
Now, considering the optimization metaheuristic to maximize the fundamental frequency and to minimize the maximum resultant displacement, taking into account the exponent of the volume fraction expression as a design variable, the results presented in Table 6 are obtained.
To note that in the second and the fourth column one has the exponents of the volume fraction law obtained for the two optimization processes. The first process corresponding to the minimization of the maximum displacement (fmin), and the second process related to the fundamental frequency maximization (fmax).
Figure 8 and Figure 9 show the evolution of Young’s modulus for the best solution found for minimizing the maximum resultant displacement and maximizing the fundamental frequency of the 2D truss, respectively, including a schematic representation of the structure’s configuration in relation to the constituent materials of the FGM.
Regarding the results obtained, it should be noted that in a situation where the maximum resulting displacement is to be minimized, the optimum solution in terms of the power law exponent is the value 10, corresponding to the truss composed mostly of zirconia in the lower zone up to a height of approximately 0.4 m. On the other hand, if the fundamental frequency of the structure is to be maximized, the best dynamic behavior in the free vibration regime is achieved with an exponent of approximately 0.19.
The optimal design configurations for this structure are consistent with the results obtained in the parametric studies.
  • Case 4—2D Frame—Analysis and Optimization
The second structure studied is a 2D frame constituted by 12 beams, the discretization of which, as well as other aspects such as geometrical configuration, boundary conditions and applied loading (Px = 50 N and q = 1000 N/m), are shown in Figure 10a,b.
Table 7 shows the resultant displacement of node 10, as this is the node where the greatest resultant displacement occurs, and the first natural frequency of the structure.
These results are again obtained by carrying out a parametric study for a set of pre-established exponent values. Now performing the optimization for each single-objective function, one obtains the results presented in Table 8.
Note that also in this table the second column contains the exponents of the volume fraction law obtained in the minimization of the maximum displacement, and the fourth column presents the exponents obtained in the fundamental frequency maximization.
Figure 11 and Figure 12 show the evolution of the modulus of elasticity for the best solutions found to minimize the maximum resultant displacement and to maximize the fundamental frequency of the frame structure, respectively, including a schematic representation of the structure’s configuration concerning the constituent materials of the FGM.
Similarly to the cases previously studied, the optimum solution that characterizes a maximum resulting displacement in frame 2 is a power law exponent of 10. From the perspective of maximizing the fundamental frequency of the structure, it should be noted that the best configuration is the one that corresponds to an approximately linear distribution between the two materials that make up the FGM. Also, for this frame structure, the optimal design configurations are consistent with the results obtained in the parametric studies.

6.3. Optimization Case Studies—One Design Variable

The current sub-section presents the results of the optimization studies conducted for two other structures in the three-dimensional space, also considering similar objective functions and the same design variable as in the two previous cases. Similar characteristics and optimization parameters continue to be used in the following cases.
  • Case 5—3D Structure I—Optimization
This case study considers a structure whose configuration can be observed in Figure 13. The structure is made of three bars joined at one point under a downward force of 1000 N.
Table 9 shows the spatial coordinates of the 3D truss supports.
The results obtained for each of the objective functions, are presented in Table 10. The two columns with the exponents ey are respectively related to the minimization of the maximum displacement fmin and the maximization of the fundamental frequency fmax.
Figure 14 and Figure 15 show the evolution of the modulus of elasticity characterized by the best solution found for minimizing the maximum resultant displacement and for maximizing the fundamental frequency of the 3D truss, respectively. These figures include a schematic representation of the structure’s configuration showing the constituent material mixture of the FGM throughout the structure height.
For this structure, it should be noted that for a situation in which the aim is to both minimize the resulting displacements and to maximize the fundamental frequency, the most advantageous distribution of properties is characterized by an exponent 10, in which zirconia makes up the entire structure up to a height of around 0.9 m.
  • Case 6—3D Structure II—Optimization
In this study, one considers a 3D frame-type structure as shown in Figure 16. This structure is subjected to optimization processes considering the same objective functions as previous cases. The design variable is also the exponent of the power law volume fraction function. The structure is discretized into 15 beam elements 1 m long (aspect ratio 20) and it is subjected to a concentrated load P = 1000 N in the static case.
The results obtained for each objective function are presented in Table 11, where the two columns with the exponents ey are respectively related to the minimization of the maximum displacement fmin and to the maximization of the fundamental frequency fmax.
Figure 17 and Figure 18 show the evolution of the modulus of elasticity associated with the best solution found for minimizing the maximum resultant displacement and for maximizing the fundamental frequency of the 3D frame, respectively. These figures include a schematic representation of the structure’s configuration with the evolution of the constituent material mixture throughout the structure height.
Once again, when the objective is to minimize the resulting displacement, the exponent that enables the best static behaviour of the structure is 10. If the objective function is the maximization of the fundamental frequency, the power law exponent that characterizes the best free vibration behaviour is, for the present structure, approximately 3.2.

6.4. Optimization Case Studies—Four Design Variables

In this subsection the optimization processes the structures will undergo have similar objective functions, but now consider four design variables as described in the Methodology section. Besides the exponent of the power law volume fraction, the next tables consider three additional columns whose acronyms are RM, M1 and M2. RM corresponds to how the rule of mixtures is considered, i.e., if a specific random variable assumes the value 1, the volume fraction characterizes the distribution of metallic material; if it assumes the value 2, then the volume fraction describes the ceramic material distribution. Regarding the acronym M1, this is related to the third design variable which defines the metallic constituent material to be used. The fourth design variable, M2, defines the ceramic material to be used.
  • Case 7—2D Frame-type Structure—Optimization
The optimization results for the 2D frame structure previously analyzed with only one design variable, but now considering four design variables, are presented in Table 12 and Table 13. These tables denote respectively the results achieved in the minimization of the maximum resultant displacement and the maximization of the fundamental frequency.
Figure 19 and Figure 20 show the evolution of Young’s modulus through the structure for the best solutions found in both optimization processes.
The optimum configuration found with four design variables that yield the lowest maximum resultant displacement corresponds to a frame made of titanium in the lower zone and tungsten carbide in the upper zone, reducing the maximum resultant displacement by around six times. The optimum solution found by the RFO for the fundamental frequency maximization corresponds to the frame made of tungsten carbide in the lower zone and aluminium in the upper zone of the structure, with an exponent of 0.1, and leads to a better solution when compared with the optimal solution using one design variable.
  • Case 8—3D Frame-type Structure—Optimization
This case study is the last of this work and considers the 3D frame already analysed in Case 4, but now with a broader perspective considering four design variables. Table 14 and Table 15 show the results obtained in the minimization process of the maximum resultant displacement and the maximization of the fundamental frequency of the 3D frame, respectively.
Figure 21 and Figure 22 depict the evolution of Young’s modulus for the best solution found for minimizing the maximum resultant displacement and maximizing the fundamental frequency of the 3D frame, respectively, including also schematic representations of the FGM evolution throughout the structure’s height.
The best solution found for the maximum displacement minimization corresponds to the combination of titanium next to the encasements and tungsten carbide in the upper area of the structure, obtaining a maximum displacement approximately four times lower than the result achieved with one design variable. To maximize the fundamental frequency, the optimum solution obtained with the four design variables corresponds to a combination of tungsten carbide next to the encasements and aluminum in the upper area of the structure, with an exponent of approximately 1.5. This leads to a fundamental frequency around 1.35 times higher than that found with the one-design variable optimization procedure.

7. Conclusions

In the present work, the optimization of a set of 2D and 3D discrete structures was conducted concerning the minimization of the maximum resulting displacement and the maximization of the fundamental frequency, in separate single-objective optimization processes. Only side constraints were considered.
Considering the first stage of the optimization studies where only one design variable was considered, the results achieved make it possible to conclude that, when the objective is the minimization of the maximum resulting displacement, in a situation where the structure is entirely metallic at the highest coordinate and entirely ceramic at the lowest coordinate, the progression between the two materials should correspond to a power law exponent of 10, which provides a stiffer structure. When the objective is to maximize the fundamental frequency, it is not possible to generalize, and it was found that for the 2D truss, the exponent should be approximately 0.18; while for the 2D frame, the best solution corresponds to an exponent of approximately 1.4. For the 3D truss it should be 10 and, for the 3D frame, the best solution found corresponds to an exponent of approximately 1.5.
When four design variables are considered, the benefits achieved in terms of reducing the resultant displacements are clear, in some cases very significantly, and in the frequencies’ increase for most cases as seen during the cases’ presentation. These improvements are possible through the optimal combination and distribution of the different materials, along with the exponent values assumed; however, it should not be generalized optimal solutions for all structures, as each structure should be considered on a case-by-case basis.
It is worth noting, the good performance of the implemented metaheuristic optimization technique and its coupling with the finite element method for FGM structures’ optimization.
Also pertinent to say that, even if a set of other non-functionally graded materials, such as isotropic and homogeneous materials are to be used in building the analysed structures, the present results can provide insightful indications about the disposition of such materials’ constant propertiesalong the structure’s height, for nearly optimized responses.
Following the promising optimal solutions found, will be relevant to pursue further optimization studies, considering a multiobjective approach for behavioural-constrained problems.

Author Contributions

Conceptualisation, M.A.R.L. and J.I.B. software, J.S.D.G. and M.A.R.L. formal analysis, J.S.D.G., M.A.R.L. and J.I.B. data curation, J.S.D.G. writing—original draft preparation, J.S.D.G. writing—review and editing, J.S.D.G., M.A.R.L. and J.I.B. visualisation, J.S.D.G. supervision, M.A.R.L. and J.I.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data generated or analysed are in the manuscript itself.

Acknowledgments

The authors acknowledge the support by FCT, through IDMEC, under LAETA, project UIDB/50022/2020.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Table A1. The complete set of 100 runs conducted for the second verification case.
Table A1. The complete set of 100 runs conducted for the second verification case.
RunfminRunfminRunfminRunfmin
11.4522 × 10−8262.4419 × 10−8512.2693 × 10−7763.0780 × 10−8
21.5391 × 10−10272.1011 × 10−6524.3243 × 10−10772.0209 × 10−10
33.8610 × 10−7285.2479 × 10−9533.4839 × 10−10782.0575 × 10−8
42.5061 × 10−10293.4899 × 10−8541.7969 × 10−6795.4393 × 10−8
59.6072 × 10−9302.1826 × 10−7551.5886 × 10−9801.2861 × 10−9
61.0757 × 10−8311.0316 × 10−6564.3608 × 10−7811.2005 × 10−7
71.1787 × 10−8321.0324 × 10−8575.3363 × 10−7825.4759 × 10−7
82.3982 × 10−9331.5102 × 10−9583.9452 × 10−9831.9531 × 10−7
96.9128 × 10−9344.4206 × 10−7591.4520 × 10−7841.1211 × 10−8
103.0236 × 10−7351.8854 × 10−8601.7240 × 10−8852.7878 × 10−10
118.3510 × 10−9361.6673 × 10−8616.4870 × 10−7866.8652 × 10−7
124.2426 × 10−11377.7937 × 10−10621.5427 × 10−8874.2135 × 10−9
131.1228 × 10−7381.2790 × 10−10638.1588 × 10−9881.1548 × 10−7
145.5395 × 10−7399.6050 × 10−9642.8892 × 10−11894.2666 × 10−8
159.7326 × 10−9401.9212 × 10−6652.6553 × 10−8903.5685 × 10−8
168.1542 × 10−9411.1275 × 10−8663.0041 × 10−8919.1252 × 10−8
172.2786 × 10−7427.4339 × 10−7676.5576 × 10−7923.1998 × 10−10
181.4118 × 10−11432.2533 × 10−7683.5742 × 10−9931.2050 × 10−9
199.3447 × 10−8441.2538 × 10−7696.6478 × 10−8941.6008 × 10−10
204.3996 × 10−7453.6096 × 10−9707.8763 × 10−11953.9891 × 10−8
211.2128 × 10−8467.1873 × 10−10713.1665 × 10−8964.4292 × 10−8
225.5121 × 10−9474.2522 × 10−10727.5925 × 10−9975.8937 × 10−8
234.6802 × 10−10481.2298 × 10−7737.1936 × 10−8989.8581 × 10−8
248.2569 × 10−10492.2679 × 10−10744.9475 × 10−10992.0966 × 10−10
251.2069 × 10−8504.2401 × 10−10751.0763 × 10−111001.0812 × 10−10
Present average value1.6230 × 10−7
Average value [22]7.5400 × 10−8
Present standard deviation3.7069 × 10−7
Standard deviation [22]2.7800 × 10−8

References

  1. Reddy, J.N. Analysis of functionally graded plates. Int. J. Numer. Methods Eng. 2000, 47, 663–684. [Google Scholar] [CrossRef]
  2. Maalawi, K. Optimization of Functionally Graded Material Structures: Some Case Studies. In Optimum Composite Structures; IntechOpen: London, UK, 2019. [Google Scholar] [CrossRef]
  3. Loja, M.A.R.; Barbosa, J.I.; Soares, C.M.M. A study on the modeling of sandwich functionally graded particulate composites. Compos. Struct. 2012, 94, 2209–2217. [Google Scholar] [CrossRef]
  4. Pasha, A.; Rajaprakash, B.M. Functionally graded materials (FGM) fabrication and its potential challenges & applications. Mater. Today Proc. 2021, 52, 413–418. [Google Scholar] [CrossRef]
  5. Markworth, A.J.; Saunders, J.H. A model of structure optimization for a functionally graded material. Mater. Lett. 1995, 22, 103–107. [Google Scholar] [CrossRef]
  6. Gasik, M.M. Micromechanical modelling of functionally graded materials. Comput. Mater. Sci. 1998, 13, 42–55. [Google Scholar] [CrossRef]
  7. Chalal, H.; Abed-Meraim, F. Quadratic Solid–Shell Finite Elements for Geometrically Nonlinear Analysis of Functionally Graded Material Plates. Materials 2018, 11, 1046. [Google Scholar] [CrossRef]
  8. Elishakoff, I.; Pentaras, D.; Gentilini, C. Mechanics of Functionally Graded Material Structures; World Scientific Pub Co. Pte Ltd.: Singapore, 2015. [Google Scholar]
  9. Wang, W.; Yuan, H.; Li, X.; Shi, P. Stress Concentration and Damage Factor Due to Central Elliptical Hole in Functionally Graded Panels Subjected to Uniform Tensile Traction. Materials 2019, 12, 422. [Google Scholar] [CrossRef]
  10. Sayyad, A.S.; Ghugal, Y.M. Bending, buckling and free vibration of laminated composite and sandwich beams: A critical review of literature. Compos. Struct. 2017, 171, 486–504. [Google Scholar] [CrossRef]
  11. Loja, M.R.; Soares, C.M. Modelling and design of adaptive structures using B-spline strip models. Compos. Struct. 2002, 57, 245–251. [Google Scholar] [CrossRef]
  12. Loja, M.A.R.; Soares, C.M.M.; Soares, C.A.M. Higher-order B-spline finite strip model for laminated adaptive structures. Compos. Struct. 2001, 52, 419–427. [Google Scholar] [CrossRef]
  13. Zhang, N.; Khan, T.; Guo, H.; Shi, S.; Zhong, W.; Zhang, W. Functionally Graded Materials: An Overview of Stability, Buckling, and Free Vibration Analysis. Adv. Mater. Sci. Eng. 2018, 2019, 1–18. [Google Scholar] [CrossRef]
  14. Mota, A.F.; Loja, M.A.R.; Barbosa, J.I.; Rodrigues, J.A. Porous Functionally Graded Plates: An Assessment of the Influence of Shear Correction Factor on Static Behavior. Math. Comput. Appl. 2020, 25, 25. [Google Scholar] [CrossRef]
  15. Neves, A.M.A.; Ferreira, A.; Carrera, E.; Cinefra, M.; Roque, C.; Jorge, R.M.N.; Soares, C.M.M. Static, free vibration and buckling analysis of isotropic and sandwich functionally graded plates using a quasi-3D higher-order shear deformation theory and a meshless technique. Compos. Part B Eng. 2013, 44, 657–674. [Google Scholar] [CrossRef]
  16. Loja, M.A.R.; Barbosa, J.I.; Soares, C.M.M. Dynamic behaviour of soft core sandwich beam structures using kriging-based layerwise models. Compos. Struct. 2015, 134, 883–894. [Google Scholar] [CrossRef]
  17. Merdaci, S.; Adda, H.M.; Hakima, B.; Dimitri, R.; Tornabene, F. Higher-Order Free Vibration Analysis of Porous Functionally Graded Plates. J. Compos. Sci. 2021, 5, 305. [Google Scholar] [CrossRef]
  18. Ergun, A.; Yu, X.; Valdevit, A.; Ritter, A.; Kalyon, D.M. Radially and Axially Graded Multizonal Bone Graft Substitutes Targeting Critical-Sized Bone Defects from Polycaprolactone/Hydroxyapatite/Tricalcium Phosphate. Tissue Eng. Part A 2012, 18, 2426–2436. [Google Scholar] [CrossRef]
  19. Muhammed, S.A.; Al-Khafaji, A.M.; Al-Deen, H.H.J.J. The Influence of Strontium Oxide on the Physio-Mechanical Properties of Biomedical-Grade Titanium in Ti-SrO Composites. J. Compos. Sci. 2023, 7, 449. [Google Scholar] [CrossRef]
  20. Yeslam, H.E. Flexural Behavior of Biocompatible High-Performance Polymer Composites for CAD/CAM Dentistry. J. Compos. Sci. 2023, 7, 270. [Google Scholar] [CrossRef]
  21. Hashim, F.A.; Hussain, K.; Houssein, E.H.; Mabrouk, M.S.; Al-Atabany, W. Archimedes optimization algorithm: A new metaheuristic algorithm for solving optimization problems. Appl. Intell. 2020, 51, 1531–1551. [Google Scholar] [CrossRef]
  22. Połap, D.; Woźniak, M. Red fox optimization algorithm. Expert Syst. Appl. 2020, 166, 114107. [Google Scholar] [CrossRef]
  23. Braik, M.S. Chameleon Swarm Algorithm: A bio-inspired optimizer for solving engineering design problems. Expert Syst. Appl. 2021, 174, 114685. [Google Scholar] [CrossRef]
  24. Nikbakht, S.; Kamarian, S.; Shakeri, M. A review on optimization of composite structures Part II: Functionally graded materials. Compos. Struct. 2019, 214, 83–102. [Google Scholar] [CrossRef]
  25. Birman, V.; Byrd, L.W. Modeling and Analysis of Functionally Graded Materials and Structures. Appl. Mech. Rev. 2007, 60, 195–216. [Google Scholar] [CrossRef]
  26. Santos, J.; Sohouli, A.; Suleman, A. Micro- and Macro-Scale Topology Optimization of Multi-Material Functionally Graded Lattice Structures. J. Compos. Sci. 2024, 8, 124. [Google Scholar] [CrossRef]
  27. Loja, M.A.R. On the use of particle swarm optimization to maximize bending stiffness of functionally graded structures. J. Symb. Comput. 2014, 61–62, 12–30. [Google Scholar] [CrossRef]
  28. 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. [Google Scholar] [CrossRef]
  29. Jaszcz, A.; Połap, D.; Damaševičius, R. Lung X-ray Image Segmentation Using Heuristic Red Fox Optimization Algorithm. Sci. Program. 2022, 2022, 1–8. [Google Scholar] [CrossRef]
  30. Bhuvaneshwari, M.; Kanaga, E.G.M.; Anitha, J. Bio-inspired Red Fox-Sine cosine optimization for the feature selection of SSVEP-based EEG signals for BCI applications. Biomed. Signal Process. Control. 2023, 80. [Google Scholar] [CrossRef]
  31. Gopi, P.S.S.; Karthikeyan, M. Red fox optimization with ensemble recurrent neural network for crop recommendation and yield prediction model. Multimedia Tools Appl. 2023, 83, 13159–13179. [Google Scholar] [CrossRef]
  32. Bamikole, J.O.; Narasigadu, C. Application of Pathfinder, Honey Badger, Red Fox and Horse Herd algorithms to phase equilibria and stability problems. Fluid Phase Equilibria 2023, 566. [Google Scholar] [CrossRef]
  33. Zhang, M.; Xu, Z.; Lu, X.; Liu, Y.; Xiao, Q.; Taheri, B. An optimal model identification for solid oxide fuel cell based on extreme learning machines optimized by improved Red Fox Optimization algorithm. Int. J. Hydrogen Energy 2021, 46, 28270–28281. [Google Scholar] [CrossRef]
  34. Gaspar, J.S.D.; Loja, M.A.R.L.; Barbosa, J.I. Static and Free Vibration Analyses of Functionally Graded Plane Structures. J. Compos. Sci. 2023, 7, 377. [Google Scholar] [CrossRef]
  35. Mohammed, H.; Rashid, T. FOX: A FOX-inspired optimization algorithm. Appl. Intell. 2023, 53, 1030–1050. [Google Scholar] [CrossRef]
  36. Reddy, J.N. Mechanics of Laminated Composite Plates and Shells: Theory and Analysis; CRC Press: Boca Raton, FL, USA, 2004. [Google Scholar]
  37. Loja, M.A.R.; Rzeszut, K.; Barbosa, J.I. Nonlocal Free Vibrations of Metallic FGM Beams. J. Compos. Sci. 2022, 6, 125. [Google Scholar] [CrossRef]
  38. Luo, Y. An efficient 3D Timoshenko beam element with consistant shape functions An Efficient 3D Timoshenko Beam Element with Consistent Shape Functions. Adv. Theor. Appl. Mech. 2008, 1, 95–106. Available online: https://www.researchgate.net/publication/259345543 (accessed on 7 March 2023).
  39. Bazoune, A.; Khulief, Y.; Stephen, N. Shape functions of three-dimensional timoshenko beam element. J. Sound Vib. 2003, 259, 473–480. [Google Scholar] [CrossRef]
  40. Sasa, G. Development of a New 3D Beam Finite Element with Deformable Section. Ph.D. Thesis, Université de Lyon, Lyon, France, 2017. Available online: http://theses.insa-lyon.fr/publication/2017LYSEI026/these.pdf (accessed on 5 March 2023).
  41. Gherlone, M.; Cerracchio, P.; Mattone, M.; Di Sciuva, M.; Tessler, A. Shape sensing of 3D frame structures using an inverse Finite Element Method. Int. J. Solids Struct. 2012, 49, 3100–3112. [Google Scholar] [CrossRef]
  42. Khebizi, M.; Guenfoud, H.; Guenfoud, M.; El Fatmi, R. Three-dimensional modelling of functionally graded beams using Saint-Venant’s beam theory. Struct. Eng. Mech. 2019, 72, 257–273. [Google Scholar] [CrossRef]
Figure 1. Schematic Representation of a Structure Built in FGM with a Variation in Properties as a Function of the Structure’s Height.
Figure 1. Schematic Representation of a Structure Built in FGM with a Variation in Properties as a Function of the Structure’s Height.
Jcs 08 00205 g001
Figure 2. Quadratic Beam Element and its Degrees of Freedom.
Figure 2. Quadratic Beam Element and its Degrees of Freedom.
Jcs 08 00205 g002
Figure 3. Schematic diagram representation of the interaction between the finite element calculation and the optimization algorithm in the one-design variable model.
Figure 3. Schematic diagram representation of the interaction between the finite element calculation and the optimization algorithm in the one-design variable model.
Jcs 08 00205 g003
Figure 4. Schematic diagram representation of the interaction between the finite element calculation and the optimization algorithm in the four-design variable model.
Figure 4. Schematic diagram representation of the interaction between the finite element calculation and the optimization algorithm in the four-design variable model.
Jcs 08 00205 g004
Figure 5. Flowchart of the Red Fox Algorithm.
Figure 5. Flowchart of the Red Fox Algorithm.
Jcs 08 00205 g005
Figure 6. Schematic Representation of 2D Truss.
Figure 6. Schematic Representation of 2D Truss.
Jcs 08 00205 g006
Figure 7. Minimum Discretization of Two-dimensional Truss.
Figure 7. Minimum Discretization of Two-dimensional Truss.
Jcs 08 00205 g007
Figure 8. Evolution of the Young’s Modulus of the 2D Truss for the Optimal Solution Found for Minimizing the Maximum Resulting Displacement using One Design Variable.
Figure 8. Evolution of the Young’s Modulus of the 2D Truss for the Optimal Solution Found for Minimizing the Maximum Resulting Displacement using One Design Variable.
Jcs 08 00205 g008
Figure 9. Evolution of the Young’s Modulus of the 2D Truss for the Optimal Solution Found for Maximizing the Fundamental Frequency using the One Design Variable.
Figure 9. Evolution of the Young’s Modulus of the 2D Truss for the Optimal Solution Found for Maximizing the Fundamental Frequency using the One Design Variable.
Jcs 08 00205 g009
Figure 10. (a) Schematic Representation of 2D Frame; (b) Minimum 2D Frame Discretization.
Figure 10. (a) Schematic Representation of 2D Frame; (b) Minimum 2D Frame Discretization.
Jcs 08 00205 g010
Figure 11. Evolution of the Modulus of Elasticity of the 2D Frame by the Best Solution Found for Minimizing the Maximum Resulting Displacement using One Design Variable.
Figure 11. Evolution of the Modulus of Elasticity of the 2D Frame by the Best Solution Found for Minimizing the Maximum Resulting Displacement using One Design Variable.
Jcs 08 00205 g011
Figure 12. Evolution of the Modulus of Elasticity of the 2D Frame by the Best Solution Found for Maximizing the Fundamental Frequency with One Design Variable.
Figure 12. Evolution of the Modulus of Elasticity of the 2D Frame by the Best Solution Found for Maximizing the Fundamental Frequency with One Design Variable.
Jcs 08 00205 g012
Figure 13. Schematic representation of three-dimensional structure under downward load.
Figure 13. Schematic representation of three-dimensional structure under downward load.
Jcs 08 00205 g013
Figure 14. Evolution of the Young’s Modulus of the 3D Truss for the Best Solution Found for Minimizing the Maximum Resulting Displacement using One Design Variable.
Figure 14. Evolution of the Young’s Modulus of the 3D Truss for the Best Solution Found for Minimizing the Maximum Resulting Displacement using One Design Variable.
Jcs 08 00205 g014
Figure 15. Evolution of the Young’s Modulus of the 3D Truss for the Best Solution Found for Maximizing the Fundamental Frequency using One Design Variable.
Figure 15. Evolution of the Young’s Modulus of the 3D Truss for the Best Solution Found for Maximizing the Fundamental Frequency using One Design Variable.
Jcs 08 00205 g015
Figure 16. Three-dimensional Frame.
Figure 16. Three-dimensional Frame.
Jcs 08 00205 g016
Figure 17. Evolution of the Modulus of Elasticity of the 3D Frame for the Best Solution Found for Minimizing the Maximum Resulting Displacement using One Design Variable.
Figure 17. Evolution of the Modulus of Elasticity of the 3D Frame for the Best Solution Found for Minimizing the Maximum Resulting Displacement using One Design Variable.
Jcs 08 00205 g017
Figure 18. Evolution of the Young’s Modulus of the 3D Frame for the Best Solution Found for Maximizing Natural Frequency using One Design Variable.
Figure 18. Evolution of the Young’s Modulus of the 3D Frame for the Best Solution Found for Maximizing Natural Frequency using One Design Variable.
Jcs 08 00205 g018
Figure 19. Evolution of Young’s Modulus of the 2D Frame for the Best Solution found for Minimizing the Maximum Resulting Displacement using four Design Variables.
Figure 19. Evolution of Young’s Modulus of the 2D Frame for the Best Solution found for Minimizing the Maximum Resulting Displacement using four Design Variables.
Jcs 08 00205 g019
Figure 20. Evolution of Young’s Modulus of 2D Frame for the Best Solution found for Maximizing Natural Frequency using four Design Variables.
Figure 20. Evolution of Young’s Modulus of 2D Frame for the Best Solution found for Maximizing Natural Frequency using four Design Variables.
Jcs 08 00205 g020
Figure 21. Evolution of Young’s Modulus of the 3D Frame for the Best Solution for Minimizing the Maximum Resulting Displacement using Four Design Variables.
Figure 21. Evolution of Young’s Modulus of the 3D Frame for the Best Solution for Minimizing the Maximum Resulting Displacement using Four Design Variables.
Jcs 08 00205 g021
Figure 22. Evolution of Young’s Modulus of the 3D Frame for the Best Solution Found for Maximizing the Fundamental Frequency using Four Design Variables.
Figure 22. Evolution of Young’s Modulus of the 3D Frame for the Best Solution Found for Maximizing the Fundamental Frequency using Four Design Variables.
Jcs 08 00205 g022
Table 1. Characteristics of the Benchmark Function “Sum of Squares”.
Table 1. Characteristics of the Benchmark Function “Sum of Squares”.
Function’s Name f Dim x m i n x m a x f m i n x ¯ *
Sum of Squares i = 1 N i · x i 2 20−10100 ( 0 , , 0 )
Table 2. Results of the Verification Study of the Sum of Squares Function with Dimension 20, using 20 Runs, 100 Iterations and Populations of 100 Individuals.
Table 2. Results of the Verification Study of the Sum of Squares Function with Dimension 20, using 20 Runs, 100 Iterations and Populations of 100 Individuals.
RunfminRunfmin
11.4522 × 10−8118.3510 × 10−9
21.5391 × 10−10124.2426 × 10−11
33.8610 × 10−7131.1228 × 10−7
42.5061 × 10−10145.5395 × 10−7
59.6072 × 10−9159.7326 × 10−9
61.0757 × 10−8168.1542 × 10−9
71.1787 × 10−8172.2786 × 10−7
82.3982 × 10−9181.4118 × 10−11
96.9128 × 10−9199.3447 × 10−8
103.0236 × 10−7204.3996 × 10−7
Present average value1.6230 × 10−7
Average value [22]7.5400 × 10−8
Present standard deviation3.7069 × 10−7
Standard deviation [22]2.7800 × 10−8
Table 3. Characteristics of the Benchmark Function “Rastrigin”.
Table 3. Characteristics of the Benchmark Function “Rastrigin”.
Function’s Name f Dim x m i n x m a x f m i n x ¯ *
Rastrigin 10 + i = 1 N [ x i 2 10 cos ( 2 π x i   ) ] 2−10100 ( 0 , , 0 )
Table 4. Results of the Verification Study of the Rastrigin Function with Dimension 2, using 10 Runs, 100 Iterations and Populations of 100 Individuals.
Table 4. Results of the Verification Study of the Rastrigin Function with Dimension 2, using 10 Runs, 100 Iterations and Populations of 100 Individuals.
Runx1x2fmin
12.1372 × 10−91.9420 × 10−109.1367 × 10−16
2−6.7385 × 10−8−8.2852 × 10−82.2627 × 10−12
37.1884 × 10−111.3883 × 10−104.8489 × 10−18
4−7.0627 × 10−106.3963 × 10−101.8013 × 10−16
55.6276 × 10−10−3.9101 × 10−109.3162 × 10−17
6−7.7343 × 10−114.0890 × 10−111.5185 × 10−18
7−2.4210 × 10−11−7.3662 × 10−111.1928 × 10−18
8−5.8023 × 10−9−8.7128 × 10−92.1740 × 10−14
9−6.1365 × 10−111.5509 × 10−105.5191 × 10−18
10−1.4728 × 10−14−2.3298 × 10−131.0812 × 10−23
Present average value2.2856 × 10−13
Average value [22]1.3740 × 10−7
Present standard deviation6.7808 × 10−13
Standard deviation [22]3.7200 × 10−8
Table 5. Resultant Displacement of Node 19 in 2D Truss and Fundamental Frequency of the Truss.
Table 5. Resultant Displacement of Node 19 in 2D Truss and Fundamental Frequency of the Truss.
eyDisplacement [μm]Frequency [rad/s]
0321.7973591.9599
0.1204.5508661.0307
0.2195.8514661.8260
0.5181.5643658.2207
1170.3850651.6240
2160.8003643.8601
5152.6575634.9667
10149.6053630.3244
Table 6. Results of the Optimization Study for 2D Truss using RFO with one design variable.
Table 6. Results of the Optimization Study for 2D Truss using RFO with one design variable.
Runeyfmin [μm]eyfmax [rad/s]
110149.60530.1609661.8009
210149.60530.2126661.7827
310149.60530.2430661.5994
410149.60530.2385661.6323
510149.60530.1764661.8437
610149.60530.1960661.8353
710149.60530.2310661.6837
810149.60530.2600661.4571
910149.60530.1641661.8135
1010149.60530.1895661.8450
Mean149.6053661.7294
SD3.0294 × 10−60.1246
Table 7. Resultant displacement of Node 10 of Frame 2 and The fundamental frequency of the structure.
Table 7. Resultant displacement of Node 10 of Frame 2 and The fundamental frequency of the structure.
eyDisplacement [μm]Frequency [rad/s]
043.54151240.7333
0.141.42791305.7998
0.240.09271346.3775
0.537.70151410.0339
135.55251444.5965
233.25061444.3502
530.21331393.7095
1028.13601356.8710
Table 8. Results of the optimization study for Frame 2 using RFO with one design variable.
Table 8. Results of the optimization study for Frame 2 using RFO with one design variable.
Runeyfmin [μm]eyfmax [rad/s]
11028.13600.91311441.6775
21028.13601.43881449.6560
37.311529.03572.24931440.2599
41028.13601.40551449.6710
51028.13601.45541449.6308
61028.13601.27421449.2112
71028.13600.96391443.4907
81028.13601.49031449.5424
91028.13601.56581449.1981
101028.13602.34011438.6473
Mean28.22601446.0985
SD0.26994.3013
Table 9. Spatial Coordinates of Three-dimensional Truss’s Simple Supports.
Table 9. Spatial Coordinates of Three-dimensional Truss’s Simple Supports.
SupportX [m]Y [m]Z [m]
10 1.5 −1
2−1 1.5 1
31 1.5 1
Table 10. Results of the Optimization Study: Minimum Displacement and Maximum Fundamental Frequency for One Design Variable.
Table 10. Results of the Optimization Study: Minimum Displacement and Maximum Fundamental Frequency for One Design Variable.
Runeyfmin [μm]eyfmax [rad/s]
1101.72668.8471255.1698
2101.72667.4426254.0038
3101.72669.8663255.8696
4101.72669.1000255.3535
5101.72668.7564255.1022
6101.72669.3028255.4959
7101.72669.1827255.4121
8101.72669.3112255.5017
9101.72667.9742254.4769
10101.72669.3355255.5184
Mean1.7266255.1904
SD00.5254
Table 11. Results of the Optimization Study of the 3D Frame using RFO with One Design Variable.
Table 11. Results of the Optimization Study of the 3D Frame using RFO with One Design Variable.
Runeyfmin [mm]eyfmax [rad/s]
1101.59303.1690121.0679
2101.59303.2185121.0680
3101.59303.2263121.0678
4101.59303.3891121.0598
5101.59303.0905121.0653
6101.59302.7019120.9957
7101.59304.4759120.8182
8101.59303.2125121.0681
9101.59303.0179121.0599
10101.59303.3001121.0657
Mean1.5930121.0336
SD00.0748
Table 12. Results of the Optimization Study for 2D Frame using the RFO Algorithm for Four design variables to Minimize the Maximum Resulting Displacement.
Table 12. Results of the Optimization Study for 2D Frame using the RFO Algorithm for Four design variables to Minimize the Maximum Resulting Displacement.
RuneyRMM1M2Optimum [μm]
11.09692226.3853
22.27322228.1455
30.12224.5157
40.12224.5157
50.12224.5157
61.82282227.5849
70.12224.5157
83.14612228.9572
90.12224.5157
100.12224.5157
Mean5.8167
SD1.6998
Table 13. Results of the Optimization Study for 2D Frame using RFO for four design variables to Maximize the Fundamental Frequency.
Table 13. Results of the Optimization Study for 2D Frame using RFO for four design variables to Maximize the Fundamental Frequency.
RuneyRMM1M2Optimum [rad/s]
12.40001221504.5482
20.43781121565.4877
31.21561221553.1407
41.04101221558.1536
50.12121574.5160
62.32041221507.6686
70.99871221559.0219
80.54961121564.4595
91.10621111447.1009
100.81121221560.2844
Mean1539.4382
SD38.2675
Table 14. Results of the Optimization Study of the 3D Frame using RFO for Four Design Variables to Minimize the Maximum Resulting Displacement.
Table 14. Results of the Optimization Study of the 3D Frame using RFO for Four Design Variables to Minimize the Maximum Resulting Displacement.
RuneyRMM1M2Optimum [μm]
110122546.1122
27.5461122551.0629
30.1222390.3244
410122546.1124
50.1222390.3247
60.1222390.3244
70.1222390.3244
80.1222390.3244
90.1222390.3244
105.5267122560.6505
Mean454.5884
SD78.79651
Table 15. Results of the Optimization Study of the 3D Frame using RFO for Four Design Variables to Maximize the Fundamental Frequency.
Table 15. Results of the Optimization Study of the 3D Frame using RFO for Four Design Variables to Maximize the Fundamental Frequency.
RuneyRMM1M2Optimum [rad/s]
11.2145112163.3825
22.4340122149.0320
31.9551122148.9432
42.9279112160.7510
52.9981112160.5673
61.5486112163.7520
75.5640122145.9512
82.2144112162.6061
92.1287111120.6251
102.8792112160.8786
Mean153.6489
SD12.7315
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gaspar, J.S.D.; Loja, M.A.R.; Barbosa, J.I. Metaheuristic Optimization of Functionally Graded 2D and 3D Discrete Structures Using the Red Fox Algorithm. J. Compos. Sci. 2024, 8, 205. https://doi.org/10.3390/jcs8060205

AMA Style

Gaspar JSD, Loja MAR, Barbosa JI. Metaheuristic Optimization of Functionally Graded 2D and 3D Discrete Structures Using the Red Fox Algorithm. Journal of Composites Science. 2024; 8(6):205. https://doi.org/10.3390/jcs8060205

Chicago/Turabian Style

Gaspar, J. S. D., M. A. R. Loja, and J. I. Barbosa. 2024. "Metaheuristic Optimization of Functionally Graded 2D and 3D Discrete Structures Using the Red Fox Algorithm" Journal of Composites Science 8, no. 6: 205. https://doi.org/10.3390/jcs8060205

Article Metrics

Back to TopTop