1. Introduction {#sec1}
===============

In structural engineering, most design optimization problems are highly nonlinear consisting of different design variables and complex constraints such as displacements, geometrical configuration, stresses, and load carrying capability. The design variables are normally grouped into two categories, namely, continuous variables and discrete variables. Optimization problems involving continuous and discrete variables generally require problem-specific search techniques \[[@B11]\]. Evolutionary multiobjective optimization techniques are examples of problem-specific search techniques. Several literatures have applied evolutionary multiobjective optimization techniques to solving multiobjective optimization problems to find a set of trade-off optimal solutions. Since most engineering problems involve multiobjective optimization, it is appropriate to apply an evolutionary optimization algorithm to solve them.

In the last two decades, different types of techniques aimed at effectively and efficiently exploring a search space by combining several basic heuristic methods have emerged \[[@B5]--[@B1]\]. These techniques currently referred to as "*Metaheuristics*" are used to describe heuristic methods applied to solving different practical problems.*Metaheuristics* can be considered as a global algorithmic framework used in solving several optimization problems with little changes, thereby making the algorithm adaptive to the specific problem \[[@B12]\].

Metaheuristic search techniques, such as simulated annealing (SA) \[[@B3]\], genetic algorithm (GA) \[[@B7]\], evolution strategies (ESs) \[[@B31]\], and particle swarm optimization (PSO) \[[@B27]\], which are generally developed based on natural phenomena have become the popular optimization techniques of recent years due to their capability of finding promising solutions for complicated optimization problems as well as their independence to the derivatives of objective functions.

Furthermore, metaheuristics can handle both discrete and real-valued variables and can be applied to a wide range of optimization problems effectively. Basically, both trajectory and population based metaheuristic approaches aim to locate the global optimum in the solution space through random moves. The key difference between the metaheuristics is in the way they propose the next move in the solution space. This motivates developers of optimization algorithms to find more efficient methodologies for originating robust optimization algorithms. However, sometimes this results in complicated approaches which are difficult to understand and implement. Hence, this study is an attempt to test the simplicity and efficiency methodology of GDE3 metaheuristics in solving engineering optimization purposes. [Section 2](#sec2){ref-type="sec"} describes the GDE3 metaheuristics briefly. Test cases are described and optimization results are discussed in [Section 3](#sec3){ref-type="sec"}. [Section 4](#sec4){ref-type="sec"} provides a clear conclusion of the study.

2. Generalized Differential Evolution Metaheuristic {#sec2}
===================================================

Several extensions of differential evolution \[[@B30]\] exist for solving constrained and nonconstrained multiobjective optimization problems \[[@B17], [@B19]\]. In comparison to the extension of differential evolution (DE), GDE3 makes differential evolution a suitable algorithm for multiobjective optimization as well as constrained optimization with little changes to the basic differential evolution algorithm. GDE3 extends DE/rand/1/bin strategy which exhibit slow convergence rates and strong exploration properties. GDE3 is a third version of generalized differential evolution modifying the selection process of the basic differential evolution algorithm \[[@B22]\]. The selection process in GDE3 is guided by these three rules:In a scenario where both the old vector and trial vector are infeasible, the old vector is selected if it dominates the trial vector, but if the trial vector weakly dominates the old vector, then the trial vector is selected.Feasible vector is selected in a situation where both feasible and infeasible vectors are generated.In a scenario where both the old vector and trial vector are feasible, the old vector is selected if it dominates the trial vector, but if the trial vector weakly dominates the old vector, then the trial vector is selected.

The whole GDE3 is presented in [Algorithm 1](#alg1){ref-type="fig"}. Parts that are new compared to previous GDE versions are framed in [Algorithm 1](#alg1){ref-type="fig"}. Without these parts, the algorithm is identical to GDE1. GDE3 can be seen as a combination of GDE2 and Pareto Differential Evolution Approach (PDEA). GDE3 is similar to differential evolution for multiobjective optimization (DEMO) except that DEMO does not contain constraint handling nor recede to basic DE in the case of a single objective because DEMO modifies the basic DE and does not consider weak dominance in the selection. Moreover, GDE3 has an improved diversity maintenance compared to DEMO. There are no constraints to be evaluated when *K* = 0 and *M* = 1, and the selection is simply$$\begin{matrix}
{x_{i,G + 1} = \begin{cases}
{u_{i,G},} & {\text{if}\,\, f\left\{ u_{i,G} \right. \leq f\left\{ x_{i,G} \right.,} \\
{x_{i,G},} & {\text{otherwise}.} \\
\end{cases}} \\
\end{matrix}$$This is the same as for the basic DE algorithm. The size of the population does not increase since this requires that *x* ~*i*,*G*~ and *u* ~*i*,*G*~ do not dominate each other even weakly, but in the case of a single objective, the reverse is the case. GDE3 performs the sorting of the vector by calculating the crowding distance of the vector. The selection process based on crowding distance gives GDE3 an advantage over NSGAII. In the case of comparing feasible, incomparable, and nondominating solutions, both offspring and parent vectors are saved for the population of the next generation \[[@B1]\]. There is no need to remove elements, since the population size does not increase. Hence, GDE3 is identical to basic DE in this case. GDE3 improves the ability to handle multiobjective optimization problems by giving a better distributed set of solutions and are less sensitive to the selection of control parameter values compared to the earlier GDE versions. As a result, this procedure reduces the computational costs of the metaheuristic and improves its efficiency. Readers interested in GDE3 should refer to the texts by \[[@B29], [@B21]\].

3. Implementation of Engineering Optimization Problems {#sec3}
======================================================

The metaheuristic optimization was implemented in NETBEAN v7.3; optimization runs were executed on an HP PC with a 2.30 GHz Intel Dual Core processor and 4 GB of RAM memory. Different examples taken from several optimization literatures were used to show the performance of GDE3 metaheuristic. These examples have been previously solved using a variety of other techniques, which is useful to show the validity and effectiveness of the GDE3 metaheuristic. The optimal results were compared with data recently published in literatures. An experiment has been performed to determine the best values of *F* and CR for better performance in GDE3 metaheuristic. For this purpose, both CR and *F* are varied from 0.1 to 1 with an increment of 0.1. The simulations were conducted for each value of *F* with respect to all values of CR. Hence, 100 such simulations were conducted. From the results, it was found that better Pareto optimal front is obtained by GDE3 with *F* = 0.5, CR = 0.9 and the termination condition is set to the 10,000 objective function evaluations.

Example 1 (welded beam design optimization problem).The welded beam problem is designed to minimize the fabrication cost by subjecting it to some constraints such as bending stress (*σ*), shear stress (*τ*), end deflection (*δ*), and buckling load (*P* ~*c*~). The design variables of the optimization problem are the thickness of the beam (*b*), the thickness of the weld (*h*), the welded joint length (*l*), and the beam width (*t*). [Figure 1](#fig1){ref-type="fig"} shows the welded beam design structure.The values of *l* and *h* must be integer multiples of 0.0065 in. Assuming *x* ~1~ =*h*, *x* ~2~ =*l*, *x* ~3~ =*t,*and *x* ~4~ =*b*as design variables, the optimization problem can be mathematically expressed as follows: $$\begin{matrix}
{\text{Minimize}\quad f\left( \overset{\rightarrow}{x} \right) = \left( 1 + C_{1} \right){x_{1}}^{2}x_{2}} \\
{\quad\quad\quad\quad\quad\quad\quad\quad\, + C_{2}x_{3}x_{4}\left( {14.0 + x_{2}} \right),} \\
{\text{Subject\ to}\quad g_{1}\left( \overset{\rightarrow}{x} \right) = \tau\left( \overset{\rightarrow}{x} \right) - \tau_{\max} \leq 0,} \\
{\,\quad\quad\quad\quad\quad g_{2}\left( \overset{\rightarrow}{x} \right) = \sigma\left( \overset{\rightarrow}{x} \right) - \sigma_{\max} \leq 0,} \\
{\,\quad\quad\quad\quad\quad g_{3}\left( \overset{\rightarrow}{x} \right) = x_{1} - x_{4} \leq 0,} \\
{\quad\quad\quad\quad\quad g_{4}\left( \overset{\rightarrow}{x} \right) = C_{1}\left( {x_{1}}^{2} \right)} \\
{\quad\quad\quad\quad\quad\quad\quad\quad\, + C_{2}x_{3}x_{4}\left( {14.0 + x_{2}} \right) - 5.0 \leq 0,} \\
{\,\quad\quad\quad\quad\quad g_{5}\left( \overset{\rightarrow}{x} \right) = 0.125 - x_{1} \leq 0,} \\
{\,\quad\quad\quad\quad\quad g_{6}\left( \overset{\rightarrow}{x} \right) = \delta\left( \overset{\rightarrow}{x} \right) - \delta_{\max} \leq 0,} \\
{\,\quad\quad\quad\quad\quad g_{7}\left( \overset{\rightarrow}{x} \right) = P - P_{c} \leq 0,} \\
\end{matrix}$$where$$\begin{matrix}
{\tau\left( \overset{\rightarrow}{x} \right) = \sqrt{\left( \tau^{\prime} \right)^{2} + \left( {2\tau^{\prime}\tau^{\prime\prime}} \right)\frac{x_{2}}{2R} + \left( \tau^{\prime\prime} \right)^{2}},} \\
{\tau^{\prime\prime} = \frac{MR}{J},} \\
{M = P\left( {L + \frac{x_{2}}{2}} \right),} \\
{R = \sqrt{\frac{{x_{2}}^{2}}{4} + \left( \frac{x_{1} + x_{3}}{2} \right)^{2}},} \\
{J = 2\left\{ {\sqrt{2}x_{1}x_{2}\left\lbrack {\frac{{x_{2}}^{2}}{12} + \left\lbrack \frac{x_{1} + x_{3}}{2} \right\rbrack^{2}} \right\rbrack} \right\},} \\
{\sigma\left( \overset{\rightarrow}{x} \right) = \frac{6PL}{x_{4}{x_{3}}^{2}},} \\
{\tau^{\prime} = \frac{P}{\sqrt{2}x_{1}x_{2}},} \\
{\delta\left( \overset{\rightarrow}{x} \right) = \frac{4PL^{3}}{E{x_{3}}^{3}x_{4}},} \\
{P_{c}\left( \overset{\rightarrow}{x} \right) = \frac{4.013E\sqrt{{x_{3}}^{2}{x_{4}}^{6}/36}}{L^{2}}\left( {1 - \frac{x_{3}}{2L}\sqrt{\frac{E}{4G}}} \right).} \\
\end{matrix}$$The simple bounds of the problem are *x* ~1~, *x* ~4~ ∈ \[0.1, 2.0\] and *x* ~2~, *x* ~3~ ∈ \[0.1, 10.0\]. The values of parameters involved in the formulation of the welded beam problem are also shown in [Table 1](#tab1){ref-type="table"}.The optimum design of the welded beam is executed using GDE3 metaheuristic, and the best solution is found as *x* ^\*^ = {*x* ~1~,  *x* ~2~,  *x* ~3~,  *x* ~4~} = {0.20572840999876, 3.47072911158159, 9.03661683005891, 0.20572540074781} which yields an objective function value of $f(\overset{\rightarrow}{x})$ = 1.7248496 as seen in [Table 2](#tab2){ref-type="table"}.The results obtained by GDE3 are presented in [Table 2](#tab2){ref-type="table"}. GDE3 found the global optimum requiring 400 iterations (i.e., 10,000 evaluations) per optimization run. [Table 3](#tab3){ref-type="table"} provides a comparison of this solution with the results of other optimization algorithms. It is apparent from the table that GDE3 metaheuristic finds a competitive solution using only 10,000 evaluations which is considerably lesser than those of other approaches. Further, a statistical evaluation of 100 independent runs of the GDE3 metaheuristic is tabulated in [Table 4](#tab4){ref-type="table"} considering the best, worst, average, and the standard deviation (std. dev.) of the obtained solutions. The ratio between the optimized costs corresponding to best and worst designs is 1.00042. Remarkably, GDE3 produced the overall best design result with a value of 1.724849. For continuous optimization problem, \[[@B26], [@B10]\] found a better design result with a value of 1.7248 at a higher function evaluation.

Example 2 (pressure vessel optimization problem).The pressure vessel problem is designed to minimize total cost which is comprised of the welding cost and forming material cost. The compressed air tank with a working pressure of 3000 psi and a minimum volume of 750 ft^3^ must be designed according to the ASME code on boilers and pressure vessels. The design variables of the optimization problem are the length of the cylindrical segment of the vessel (*L*), the thickness of the cylindrical skin (*T* ~*s*~), the inner radius (*R*), and the thickness of the spherical head (*T* ~*h*~).The variables *T* ~*s*~ and *T* ~*h*~ are discrete values which are integer multiples of 0.0625 inches. [Figure 2](#fig2){ref-type="fig"} shows the cylindrical pressure vessel capped at both ends by hemispherical heads.Assuming *x* ~1~ = *T* ~*s*~, *x* ~2~ = *T* ~*h*~, *x* ~3~ = *R*, and *x* ~4~ = *L* as the design variables, the optimization problem can be mathematically expressed as follows:$$\begin{matrix}
{\text{Minimize}\quad f\left( \overset{\rightarrow}{x} \right) = 0.6224x_{1}x_{3}x_{4} + 1.7781x_{2}{x_{3}}^{2}} \\
{\quad\quad\quad\quad\quad\quad\quad\quad + 3.1611{x_{1}}^{2}x_{4} + 19.8621{x_{1}}^{2}x_{3},} \\
{\text{Subject\ to}\quad g_{1}\left( \overset{\rightarrow}{x} \right) = 0.0193x_{3} - x_{1} \leq 0,} \\
{\,\quad\quad\quad\quad\quad g_{2}\left( \overset{\rightarrow}{x} \right) = 0.00954x_{3} - x_{2} \leq 0,} \\
{\,\quad\quad\quad\quad\quad g_{3}\left( \overset{\rightarrow}{x} \right) = x_{4} - 240 \leq 0,} \\
{\quad\quad\quad\quad\quad\, g_{4}\left( \overset{\rightarrow}{x} \right) = 750 \times 1728 - \pi{x_{3}}^{2}x_{4}} \\
{\quad\quad\quad\quad\quad\quad\,  \quad   - \frac{4}{3}\pi{x_{3}}^{3} \leq 0.} \\
\end{matrix}$$The simple bounds of the problem are *x* ~1~, *x* ~2~ ∈ \[1 × 0.0625, 99 × 0.0625\] and *x* ~3~, *x* ~4~ ∈ \[10.0, 240.0\]. Unlike the usual limit of 200 in considered in literatures, the upper bound of design variable*L* was increased to 240 in to expand the search space.Optimization results are presented in [Table 5](#tab5){ref-type="table"}. GDE3 produced a design result with a value of 6083.773 within 400 iterations (i.e., 10,000 evaluations). [Table 6](#tab6){ref-type="table"} compares the optimal design results produced by GDE3 with those reported in \[[@B11], [@B16], [@B26], [@B6], [@B14], [@B9]\]. Further, a statistical evaluation of 100 independent runs of the GDE3 metaheuristic is tabulated in [Table 7](#tab7){ref-type="table"} considering the best, worst, average, and the standard deviation (std. dev.) of the obtained solutions. The ratio between the optimized costs corresponding to worst and best designs is 1.00229. The best design result was produced by the Firefly algorithm. GDE3 metaheuristic produced the least performance compared to the other algorithms.

Example 3 (speed reducer design optimization problem).The speed reducer design problem \[[@B13]\] is designed to minimize the weight of the speed reducer subjecting it to some constraints such as shaft stresses, surface stress, gear teeth bending stress, and shafts crosswise deflections. The width of the gear face *x* ~1~, teeth module *x* ~2~, number of pinion teeth *x* ~3~, first shaft length between bearings *x* ~4~, second shaft length between bearings *x* ~5~, the diameter of the first shaft *x* ~6~, and diameter of the second shaft are the design variables of the optimization problem. [Figure 3](#fig3){ref-type="fig"} shows the schematic of the speed reducer.The mathematical expression for the speed reducer problem is as follows:$$\begin{matrix}
{\text{Minimize}\quad f\left( \overset{\rightarrow}{x} \right) = 0.7854x_{1}{x_{2}}^{2}} \\
{\quad\quad\quad\quad\quad\quad\quad  \,\, \cdot \left( {3.3333{x_{3}}^{2} + 14.9334x_{3}} \right.} \\
\left. {\quad\quad\quad\quad\quad\quad\quad\quad\quad - 43.0934\phantom{{x_{3}}^{2}}} \right) \\
{\quad\quad\quad\quad\quad\quad\quad    - 1.508x_{1}\left( {x_{6}}^{2} + {x_{7}}^{2} \right)} \\
{\quad\quad\quad\quad\quad\quad\quad\quad\, + 7.4777\left( {x_{6}}^{3} + {x_{7}}^{3} \right)} \\
{\quad\quad\quad\quad\quad\quad\quad\quad + 0.7854\left( {x_{4}{x_{6}}^{2} + x_{5}{x_{7}}^{2}} \right),} \\
{\text{Subject\ to}\quad g_{1}\left( \overset{\rightarrow}{x} \right) = \frac{27}{x_{1}{x_{2}}^{2}x_{3}} - 1 \leq 0,} \\
{\,\quad\quad\quad\quad\quad\, g_{2}\left( \overset{\rightarrow}{x} \right) = \frac{397.5}{x_{1}{x_{2}}^{2}{x_{3}}^{2}} - 1 \leq 0,} \\
{\,\quad\quad\quad\quad\,\quad g_{3}\left( \overset{\rightarrow}{x} \right) = \frac{1.93{x_{4}}^{3}}{x_{2}x_{3}{x_{6}}^{4}} - 1 \leq 0,} \\
{\,\quad\quad\quad\,\quad  \,\, g_{4}\left( \overset{\rightarrow}{x} \right) = \frac{1.93{x_{5}}^{3}}{x_{2}x_{3}{x_{7}}^{4}} - 1 \leq 0,} \\
{\quad\quad  \quad\quad\,\,\,\, g_{5}\left( \overset{\rightarrow}{x} \right) = \frac{1.0}{110{x_{6}}^{3}}} \\
{\quad\quad\quad  \quad\quad\quad\quad\,\,\,\,\, \cdot \sqrt{\left( \frac{745.0x_{4}}{x_{2}x_{3}} \right)^{2} + 16.9 \times 10^{6}}} \\
{\quad\quad\,\,\quad\quad\quad\quad\quad\quad - 1 \leq 0,} \\
{\quad\quad\quad\quad  \,\,\,\, g_{6}\left( \overset{\rightarrow}{x} \right) = \frac{1.0}{85{x_{7}}^{3}}} \\
{\quad\quad\quad\quad\quad\quad\quad\quad\,\,\, \cdot \sqrt{\left( \frac{745.0x_{5}}{x_{2}x_{3}} \right)^{2} + 157.5 \times 10^{6}}} \\
{\quad\quad\quad\quad\quad\quad\quad\quad\,\, - 1 \leq 0,} \\
{\,\quad\quad\quad  \quad\,\,\,\, g_{7}\left( \overset{\rightarrow}{x} \right) = \frac{x_{2}x_{3}}{40} - 1 \leq 0,} \\
{\,\quad\quad\quad  \quad\,\,\, g_{8}\left( \overset{\rightarrow}{x} \right) = \frac{5x_{2}}{x_{1}} - 1 \leq 0,} \\
{\,\quad\quad\quad\quad  \,\,\, g_{9}\left( \overset{\rightarrow}{x} \right) = \frac{x_{1}}{12x_{2}} - 1 \leq 0,} \\
{\,\quad\quad\quad\quad\,\quad\, g_{10}\left( \overset{\rightarrow}{x} \right) = \frac{1.5x_{6} + 1.9}{x_{4}} - 1 \leq 0,} \\
{\,\quad\quad\quad\quad  \,\,\,\, g_{11}\left( \overset{\rightarrow}{x} \right) = \frac{1.1x_{7} + 1.9}{x_{5}} - 1 \leq 0.} \\
\end{matrix}$$The simple bounds of the problem are *x* ~1~ ∈ \[2.6, 3.6\], *x* ~2~ ∈ \[0.7, 0.8\],  *x* ~3~ ∈ \[17, 28\], *x* ~4~ ∈ \[7.3, 8.3\],  *x* ~5~ ∈ \[7.8, 8.3\],  *x* ~6~ ∈ \[2.9, 3.0\], and *x* ~7~ ∈ \[5.0, 5.5\].The optimum design of the speed reducer is executed using GDE3 metaheuristic, and the best solution is found as *x* ^\*^ = {*x* ~1~,  *x* ~2~,  *x* ~3~,  *x* ~4~,  *x* ~5~,  *x* ~6~,  *x* ~7~} = {3.50000000047883, 0.7, 17.0, 7.3, 7.8, 3.35021466645262, 5.2866832298256} which yields an objective function value of $f(\overset{\rightarrow}{x})$ = 2996.34816529042 as seen in [Table 8](#tab8){ref-type="table"}.The results obtained by GDE3 are presented in [Table 8](#tab8){ref-type="table"}. GDE3 found the global optimum requiring 400 iterations per optimization run. [Table 9](#tab9){ref-type="table"} provides a comparison of this solution with the results of simple constrained particle swarm optimization. It is apparent from the table that GDE3 metaheuristic finds a competitive solution using only 10,000 objective function evaluations, which is considerably lesser than those of other approaches. Further, a statistical evaluation of 100 independent runs of the GDE3 metaheuristic is tabulated in [Table 10](#tab10){ref-type="table"} considering the best, worst, average, and the standard deviation (std. dev.) of the obtained solutions. The ratio between the optimized costs corresponding to best and worst designs is 1.0000003. Remarkably, GDE3 produced the overall best design result with a value of 2996.3481653.

Example 4 (tension/compression spring design optimization problem).The tension/compression spring problem is designed to minimize the weight of the spring subjecting it to some constraints such as shear stress, minimum deflection, outside diameter limits, and surge frequency. The design variables are the number of active coils *P*, the diameter of the mean coil *D*, and the diameter of the wire *d*. [Figure 4](#fig4){ref-type="fig"} shows the tension/compression spring design.Assuming *x* ~1~ = *d*, *x* ~2~ = *D*, and *x* ~3~ = *P*, as the design variables, the tension/compression spring design problem can be expressed as follows:$$\begin{matrix}
{\text{Minimize}\quad f\left( \overset{\rightarrow}{x} \right) = \left( {x_{3} + 2} \right)x_{2}{x_{1}}^{2},} \\
{\text{Subject\ to}\quad g_{1}\left( \overset{\rightarrow}{x} \right) = 1 - \frac{{x_{2}}^{3}x_{3}}{71,785{x_{1}}^{4}} \leq 0,} \\
{\quad\quad\quad\quad\quad\, g_{2}\left( \overset{\rightarrow}{x} \right) = \frac{4{x_{2}}^{2} - x_{1}x_{2}}{12,566\left( x_{2}{x_{1}}^{3} - {x_{1}}^{4} \right)}} \\
{\quad\quad\quad\quad\quad\quad\,  \quad   + \frac{1}{5,108{x_{1}}^{2}} - 1 \leq 0,} \\
{\,\quad\quad\quad\quad\quad\,\, g_{3}\left( \overset{\rightarrow}{x} \right) = 1 - \frac{140.45x_{1}}{{x_{2}}^{2}x_{3}} \leq 0,} \\
{\,\quad\quad\quad\quad  \,\, g_{4}\left( \overset{\rightarrow}{x} \right) = \frac{x_{2} + x_{1}}{1.5} - 1 \leq 0.} \\
\end{matrix}$$The simple bounds of the problem are *x* ~1~ ∈ \[0.05, 2.0\], *x* ~2~ ∈ \[0.25, 1.3\], and *x* ~3~ ∈ \[2.0, 15.0\].The optimum design of the tension/compression spring is carried out using GDE3 metaheuristic, and the best solution is found as *x* ^\*^ = {*x* ~1~, *x* ~2~, *x* ~3~} = {0.0517955276224998, 0.359283196922392, 11.1405163630287} which yields an objective function value of $f(\overset{\rightarrow}{x})$ = 0.0126658360085857 as seen in [Table 11](#tab11){ref-type="table"}.The results obtained by GDE3 are presented in [Table 11](#tab11){ref-type="table"}. GDE3 found the global optimum requiring 400 iterations per optimization run. [Table 12](#tab12){ref-type="table"} provides a comparison of this solution with the results of simple constrained particle swarm optimization. It is apparent from the table that GDE3 metaheuristic finds a competitive solution using only 10,000 objective function evaluations, which is considerably lesser than those of other approaches. Further, a statistical evaluation of 100 independent runs of the GDE3 metaheuristic is tabulated in [Table 13](#tab13){ref-type="table"} considering the best, worst, average, and the standard deviation (std. dev.) of the obtained solutions. The ratio between the optimized costs corresponding to worst and best designs is 1.000107.

4. Conclusion {#sec4}
=============

In the present study, the GDE3 algorithm is used as a simple and efficient optimization technique for handling engineering optimization problems. The GDE3 algorithm also uses a very simple mechanism to deal with constrained functions and results generated by the algorithm indicate that such mechanism, despite its simplicity, is effective in practice. From this study, performance evaluation of the GDE3 algorithm through benchmark design optimization examples reveals the efficiency of this technique in solving practical optimization problems. Although in the present study the algorithm is utilized only for solving engineering design optimization problems, GDE3 algorithm can easily be employed for solving other types of optimization problems as well.

Conflict of Interests
=====================

The authors declare that there is no conflict of interests regarding the publication of this paper.

![Schematic of the welded beam design problem \[[@B11]\].](TSWJ2015-936106.001){#fig1}

![Schematic of the pressure vessel design problem \[[@B11]\].](TSWJ2015-936106.002){#fig2}

![Schematic of the speed reducer design problem \[[@B13]\].](TSWJ2015-936106.003){#fig3}

![Schematic of the tension/compression spring design problem.](TSWJ2015-936106.004){#fig4}

![The GDE3 algorithm \[[@B22]\].](TSWJ2015-936106.alg.001){#alg1}

###### 

Values of parameters involved in the formulation of the welded beam problem \[[@B11]\].

  Constant item   Description                           Values
  --------------- ------------------------------------- ------------------
  *C* ~1~         The welded material                   0.10471 (\$/in3)
  *C* ~2~         The bar stock                         0.04811 (\$/in3)
  *τ* ~max⁡~      Shear stress of the welded material   13600 (psi)
  *σ* ~max⁡~      Normal stress of the bar material     30000 (psi)
  *δ* ~max⁡~      Bar end deflection                    0.25 (inch)
  *E*             Young\'s modulus of bar stock         30 × 10^6^ (psi)
  *G*             Shear modulus of bar stock            12 × 10^6^ (psi)
  *P*             Loading condition                     6000 (lb)
  *L*             Beam\'s projection length             14 (inch)

###### 

GDE3 solution vector for welded beam.

  ----------------------------------- ----------------------------------- ----------------------------------- ----------------------------------- ---------
  Best solution                       *x* ~1~                             *x* ~2~                             *x* ~3~                             *x* ~4~
  0.20572840999876                    3.47072911158159                    9.03661683005891                    0.20572540074781                    
  $g_{1}(\overset{\rightarrow}{x})$   $g_{2}(\overset{\rightarrow}{x})$   $g_{3}(\overset{\rightarrow}{x})$   $g_{4}(\overset{\rightarrow}{x})$   
  −0.66062798472194                   0.665171394633944                   3.00925094998*E* − 06               −3.43299575416113                   
  $g_{5}(\overset{\rightarrow}{x})$   $g_{6}(\overset{\rightarrow}{x})$   $g_{7}(\overset{\rightarrow}{x})$   $f(\overset{\rightarrow}{x})$       
  −0.08072840999876                   −0.235539990649711                  0.373971078704926                   1.72484969509211                    
  ----------------------------------- ----------------------------------- ----------------------------------- ----------------------------------- ---------

###### 

Welded beam problem: comparison of GDE3 results with other optimization methods.

  Researcher      Metaheuristic                                        *x* ~1~      *x* ~2~      *x* ~3~      *x* ~4~      *f*(*x*)       NE
  --------------- ---------------------------------------------------- ------------ ------------ ------------ ------------ -------------- ------------
  \[[@B8]\]       Genetic algorithm                                    0.2489       6.1730       8.1789       0.2533       2.4331         320,080
  \[[@B24]\]      Genetic algorithm                                    0.2489       6.1097       8.2484       0.2485       2.4000         6,273
  \[[@B2]\]       Social behavioral model                              0.2407       6.4851       8.2399       0.2497       2.4426         19,259
  \[[@B28]\]      Society and civilization algorithm                   0.2444       6.2380       8.2886       0.2446       2.3854         33,095
  \[[@B25]\]      Genetic algorithm                                    0.2443       6.2117       8.3015       0.2443       2.3816         320,000
  \[[@B15]\]      Particle swarm optimization                          0.2444       6.2175       8.2915       0.2444       2.3810         30,000
  \[[@B23]\]      Harmonic search                                      0.2442       6.2231       8.2915       0.2443       2.3810         110,000
  \[[@B16]\]      Simulated annealing---direct search                  0.2444       6.2158       8.2939       0.2444       2.3811         56,243
  \[[@B18]\]      Simulated annealing---genetic algorithm              0.2231       1.5815       12.8468      0.2245       2.2500         26,466
  \[[@B4]\]       Artificial Immune System---genetic algorithm         0.2444       6.2183       8.2912       0.2444       2.3812         320,000
  \[[@B26]\]      Harmonic search                                      0.2057       3.4705       9.0366       0.2057       1.7248         200,000
  \[[@B6]\]       Simple constrained particle swarm optimizer          0.2057       3.4705       9.0366       0.2057       1.7249         24,000
  \[[@B10]\]      Harmonic search---sequential quadratic programming   0.2057       3.4706       9.0368       0.2057       1.7248         90,000
  \[[@B32]\]      Differential evolution                               0.2444       6.2175       8.2915       0.2444       2.3810         24,000
  \[[@B31]\]      Evolutionary algorithm                               0.2443       6.2201       8.2940       0.2444       2.3816         28,897
  \[[@B11]\]      Firefly algorithm                                    0.2015       3.5620       9.0414       0.2057       1.7312         50,000
  \[[@B14]\]      Simple optimization                                  0.2057       3.4705       9.0366       0.2057       1.7246         10,000
  Present study   Generalized differential evolution 3                 **0.2057**   **3.4707**   **9.0366**   **0.2057**   **1.724849**   **10,000**

###### 

Statistical results of the GDE3 optimization.

  Best       Average    Worst      Std. dev.   Number of iterations
  ---------- ---------- ---------- ----------- ----------------------
  1.724849   1.725023   1.725569   0.0001018   400

###### 

GDE3 Solution vector for pressure vessel.

  ----------------------------------- ----------------------------------- ----------------------------------- ----------------------------------- ---------
  Best Solution                       *x* ~1~                             *x* ~2~                             *x* ~3~                             *x* ~4~
  0.74395291436715                    0.36774755668330                    38.5288195380221                    239.37719314082                     
  $g_{1}(\overset{\rightarrow}{x})$   $g_{2}(\overset{\rightarrow}{x})$   $g_{3}(\overset{\rightarrow}{x})$   $g_{4}(\overset{\rightarrow}{x})$   
  −0.00034669728332                   −0.00018261829057                   −0.62280685917400                   −42.436889517499                    
  $f(\overset{\rightarrow}{x})$                                                                                                                   
  6083.77328355025                                                                                                                                
  ----------------------------------- ----------------------------------- ----------------------------------- ----------------------------------- ---------

###### 

Pressure vessel problem: comparison of GDE3 results with optimization methods.

  Researcher      Metaheuristic                                     *x* ~1~   *x* ~2~   *x* ~3~   *x* ~4~    *f*(*x*)
  --------------- ------------------------------------------------- --------- --------- --------- ---------- -----------
  \[[@B16]\]      Simulated annealing---direct search               0.7683    0.3797    39.8096   207.2250   5868.76
  \[[@B9]\]       Particle swarm optimization---genetic algorithm   0.7500    0.3750    38.8601   221.3654   5850.383
  \[[@B26]\]      Harmonic search                                   0.7500    0.3750    38.8600   221.3600   5849.7
  \[[@B6]\]       Simple constrained particle swarm optimizer       0.8125    0.4375    42.0980   176.6360   6.059.714
  \[[@B11]\]      Firefly algorithm                                 0.7500    0.3750    38.8600   221.3600   5850.3
  \[[@B14]\]      Simple optimization                               1.1250    0.6250    58.2901   43.6927    7199.35
  Present study   Generalized differential evolution 3              0.74391   0.36774   38.5288   239.377    6083.773

###### 

Statistical results of the GDE3 optimization.

  Best       Average    Worst      Std. dev.   Number of iterations
  ---------- ---------- ---------- ----------- ----------------------
  6083.773   6092.318   6097.725   40.32205    400

###### 

GDE3 Solution vector for speed reducer.

  ------------------------------------ ------------------------------------ ----------------------------------- ----------------------------------- ---------
  Best Solution                        *x* ~1~                              *x* ~2~                             *x* ~3~                             *x* ~4~
  3.5000000004788                      0.7000000000000                      17.0000000000000                    7.3000000000000                     
  *x* ~5~                              *x* ~6~                              *x* ~7~                             $g_{1}(\overset{\rightarrow}{x})$   
  7.8000000000000                      3.3502146664526                      5.2866832298256                     −0.07391528052456                   
  $g_{2}(\overset{\rightarrow}{x})$    $g_{3}(\overset{\rightarrow}{x})$    $g_{4}(\overset{\rightarrow}{x})$   $g_{5}(\overset{\rightarrow}{x})$   
  −0.197998527251663                   −0.499172248315386                   −0.9014716976203                    −3.1892233298*E* − 10               
  $g_{6}(\overset{\rightarrow}{x})$    $g_{7}(\overset{\rightarrow}{x})$    $g_{8}(\overset{\rightarrow}{x})$   $g_{9}(\overset{\rightarrow}{x})$   
  −3.8408276559*E* − 11                −0.7025                              −1.3680001575*E* − 10               −0.5833333332763                    
  $g_{10}(\overset{\rightarrow}{x})$   $g_{11}(\overset{\rightarrow}{x})$   $f(\overset{\rightarrow}{x})$                                           
  −0.0513257534686439                  −0.0108523650245949                  2996.34816529042                                                        
  ------------------------------------ ------------------------------------ ----------------------------------- ----------------------------------- ---------

###### 

Speed reducer problem: comparison of generalized differential evolution 3 results with simple constrained particle swarm optimization.

  Solution                        Simple constrained particle swarm optimization \[[@B6]\]   Generalized differential evolution (Present study)
  ------------------------------- ---------------------------------------------------------- ----------------------------------------------------
  *x* ~1~                         3.5000                                                     3.5000
  *x* ~2~                         0.7000                                                     0.7000
  *x* ~3~                         17.0000                                                    17.0000
  *x* ~4~                         7.3000                                                     7.3000
  *x* ~5~                         7.8000                                                     7.8000
  *x* ~6~                         3.350214                                                   3.3502146
  *x* ~7~                         5.286683                                                   5.2866832
  $f(\overset{\rightarrow}{x})$   2996.348165                                                2996.3481653

###### 

Statistical results of the GDE3 optimization.

  Best           Average        Worst          Std. dev.   Number of iterations
  -------------- -------------- -------------- ----------- ----------------------
  2996.3481653   2996.3483815   2996.3491534   0.0000021   400

###### 

GDE3 Solution vector for tension/compression spring.

  ----------------------------------- ----------------------------------- ----------------------------------- ------------------------------- -----------------------------------
  Best Solution                       *x* ~1~                             *x* ~2~                             *x* ~3~                         $g_{1}(\overset{\rightarrow}{x})$
  0.0517955276224998                  0.359283196922392                   11.1405163630287                    −3.0601282864*E* − 05           
  $g_{2}(\overset{\rightarrow}{x})$   $g_{3}(\overset{\rightarrow}{x})$   $g_{4}(\overset{\rightarrow}{x})$   $f(\overset{\rightarrow}{x})$   
  −0.133636716257444                  −4.05865278285946                   −0.725947516970072                  0.01266583600858                
  ----------------------------------- ----------------------------------- ----------------------------------- ------------------------------- -----------------------------------

###### 

Tension/compression spring problem: comparison of GDE3 results with simple constrained particle swarm optimization.

  Solution                        Simple constrained particle swarm optimization \[[@B6]\]   Generalized differential evolution 3
  ------------------------------- ---------------------------------------------------------- --------------------------------------
  *x* ~1~                         0.051583                                                   0.0517955
  *x* ~2~                         0.354190                                                   0.3592831
  *x* ~3~                         11.438675                                                  11.140516
  $f(\overset{\rightarrow}{x})$   0.012665                                                   0.012665836

###### 

Statistical results of the GDE3 optimization.

  Best          Average       Worst         Std. dev.         Number of iterations
  ------------- ------------- ------------- ----------------- ----------------------
  0.012665836   0.012666648   0.012667194   3.97815*E* − 07   400

[^1]: Academic Editor: Fei Kang
