Next Article in Journal
Fast Switch and Spline Function Inversion Algorithm with Multistep Optimization and k-Vector Search for Solving Kepler’s Equation in Celestial Mechanics
Previous Article in Journal
Statistical Parameters Based on Fuzzy Measures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Constrained Mixed-Variable Design Optimization Based on Particle Swarm Optimizer with a Diversity Classifier for Cyclically Neighboring Subpopulations

1
Department of Mechanical Engineering, Chung-Ang University, 84 Heukseok-ro, Dongjak-gu, Seoul 06974, Korea
2
Department of Mechanical Engineering, University of Hawaii at Manoa, Honolulu, HI 96822, USA
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(11), 2016; https://doi.org/10.3390/math8112016
Submission received: 21 October 2020 / Revised: 6 November 2020 / Accepted: 7 November 2020 / Published: 12 November 2020
(This article belongs to the Section Engineering Mathematics)

Abstract

:
In this research, an easy-to-use particle swarm optimizer (PSO) for solving constrained engineering design problems involving mixed-integer-discrete-continuous (MIDC) variables that adopt two kinds of diversity-enhancing mechanisms to achieve superior reliability and validity was developed. As an initial diversity-boosting tool, the local neighborhood topology of each particle is set up such that information exchange is restricted to a limited number of consecutively numbered particles. This topological mechanism forces each particle to move in the search space while interacting only with its neighboring subpopulation. The second diversity-enhancing task is to ensure that the exploration behavior of each particle in the search space is governed such that it follows the diversity classifier decision applied to its subpopulation. This diversity classification iteratively adjusts the three-phase velocity-related mechanism of each particle such that it approaches or retreats from its previous best position/the current best position among the subpopulation. In summary, this PSO tool not only introduces the social interaction of the particle within its cyclically neighboring subpopulation but also exploits the three-phase velocity behavior law governed by the distributed diversity measures categorized for each neighboring subpopulation. This scheme has superior reliability, as well as high practicality for engineering optimization problems involving MIDC variables, which are handled by the widely adopted straightforward rounding-off technique used in most swarm-inspired metaheuristic search technologies.

1. Introduction

Several real-life engineering design tasks require the use of numerical optimization techniques that can handle highly nonlinear multimodal problems, with several complicated constraints on factors such as stress, deflection, load-carrying capability, and geometric configuration [1]. Although several gradient-based deterministic algorithms have been introduced during the past decades, their solution is a function of the initial search points and thus might not be the global optimum. This drawback of classical methods has induced the widespread adoption of biologically inspired metaheuristic algorithms; i.e., algorithms that mimic the behavior of different species or the social-cooperative behavior of swarms [1,2,3]. Some popular metaheuristic optimizers include the differential evolution (DE) algorithm, genetic algorithm (GA), particle swarm optimization (PSO), simulated annealing (SA) algorithm, and ant colony optimization. In this area of research, the PSO algorithm has attracted much attention as a reliable treatment for a wide range of large-scale, multimodal, non-convex, and nonlinear unconstrained optimization problems [4,5].
When handling engineering optimization problems via PSO, one frequently faces the following two obstacles.
(I)
Most engineering design tasks include multiple real-life physical constraints and thus necessitate PSO techniques that can account for these. Furthermore, such physical constraints are, in general, treated as hard constraints that should be satisfied by any feasible solution found via the optimization procedure.
(II)
A large fraction of design problems in engineering fields belong in the category of mixed-integer-discrete-continuous (MIDC) optimization problems and thus particular care has to be taken to find feasible solutions, which makes design tasks more challenging.
However, since PSO was initially devised not only to handle unconstrained optimization problems but also to work only with real-valued variables, the presence of non-continuous design variables results increases the difficulty of finding feasible or optimal solutions, as in most other evolutionary algorithms [6,7,8,9,10]. Trends and ways of tackling constrained mixed-variable problems in evolutionary algorithms can be found in several articles [9,11,12,13,14,15,16,17,18,19,20,21]. Nevertheless, a rigorous mechanism for handling integer/discrete variables is still challenging and the focus of intensive research, even in PSO applications for unconstrained optimization problems, as mentioned in [6].
Recently, [10] preliminarily investigated a straightforward procedure for solving constrained MIDC optimal design problems. In their research, they recommended a constrained PSO algorithm combined with the conventional technique of rounding the variable values to the nearest discrete or integer values for objective function evaluation, while all particles still freely evolve with floating-point values within a continuous search domain, regardless of the type of variable. Although this primitive technique seems to be rather unrefined, it has been widely introduced with promising application results in several studies on evolutionary algorithms [11,12,13,14,15,16,17,18,19,20,21]. Nevertheless, their constrained PSO-based method combined with the rounding-off mechanism may not guarantee acceptable reliability; i.e., the probability of finding the required design variables satisfying the optimality criteria may be insufficient. This performance degradation is certainly due to the insufficient diversification characteristics of swarms, which, in practical implementation, may accelerate the premature convergence phenomenon of particles to a local optimum. It should be noted that diversification helps the optimizer to explore the search space on a global scale, possibly resulting in the exploration of various diverse solutions for attaining global optimality [1]. However, most population-based evolutionary computation techniques, including PSO, depend on the basic underlying principle of reducing the search space in the progression toward a global optimum. Updating the position vector for each particle that outperforms the current global best position of the population could no longer be achieved after a certain iteration number owing to the swarm-movement pattern and, in this case, the whole swarm may remain trapped in a search region containing a local optimal solution [22,23]. In such a case, if a method to enhance population diversity is not introduced, particles evolving based on the conventional PSO mechanism may have difficulty reaching the true global optimal position. Such a diversity improvement is expected to motivate research studies of our diversity-guided PSO with a cyclic network mechanism.
The purpose of the present study is the realization of an easy-to-use diversity-guided PSO scheme for constrained mixed-variable optimization problems, which uses diversity classifiers for cyclically neighboring subpopulations and is characterized by superior reliability and validity. To achieve this aim, this paper develops a novel constrained PSO strategy that combines the following two kinds of diversity-enhancing mechanisms. The first diversity-boosting tool is to construct cyclically neighboring multi-subpopulations, one for each particle (host particle), where the subpopulation is composed of some successively numbered particles, with the host particles as the center. In this structure, the host particle only interacts with a limited number of neighboring particles (i.e., particles belonging to the subpopulation) during the search process within the problem-solution space. Thus, the “nearby neighbor” in this mechanism implies the particle within the subpopulation that transfers individual positional information and fitness values to the host particle [24,25,26]. Note that the host particle of a certain subpopulation can be one of several nearby neighbors of other subpopulations; thus, any information on the good fitness value of one host particle is cyclically propagated to other neighboring host particles during the evolution process. Therefore, this structure forces each particle to achieve social learning from its respective neighborhoods and not directly from the best position of the entire swarm, which results in enhancing the population diversity. Second, to achieve diversity improvement for each subpopulation, the aforementioned PSO mechanism is extended and combined with the three-phase velocity update law, which is introduced to governs the behavior of the host particles at each iteration. In this extended PSO algorithm, the velocity behavior of each host particle per iteration is governed by attraction, repulsion, and in-between phases [23,27], chosen according to the diversity classifier applied to the subpopulation of the host particle. In summary, our PSO scheme not only exploits constricted social learning of each particle through cyclically neighboring subpopulations but also incorporates the three-phase velocity behavior law, implemented by following the locally distributed diversity measures categorized for each subpopulation. This novel PSO scheme features high reliability, as well as superior practicality for engineering optimization problems involving MIDC variables, which are handled via the straightforward rounding-off technique that has been widely adopted in swarm-inspired metaheuristic search technologies. Note that the conventional PSO scheme proposed by [23,27] evaluates the sole diversity measure obtained from the position data of the entire swarm and targets optimization problems having continuous variables and no constraint functions. Therefore, their method may not reliably provide a sufficient level of population diversity to produce offspring that outperform their parents, in the case of MIDC nonlinear programming with multimodal characteristics; even when no constraint is involved. Several benchmark MIDC design problems are presented to clearly verify that the proposed PSO scheme is a highly useful tool, providing remarkable reliability.
The remainder of this paper is organized as follows. In the following section, typical MIDC engineering design problems and a brief sketch of the conventional evolutionary PSO algorithm are presented. In Section 3, the development of a novel PSO scheme with a diversity classifier for cyclically neighboring subpopulations is described, and its distinctive features are discussed in depth. The details of eight typical benchmark design examples are given in the successive sections, and then, simulation results obtained using the proposed PSO tool are compared in depth with those obtained using other PSO methods and various other evolutionary algorithms. Further discussions on the influence of variations in neighboring particle sizes on the performance of the proposed PSO scheme are also presented. Finally, we present our conclusions in Section 5.

2. Engineering Design Problem and Particle Swarm Optimizer

The classical PSO is a population-based stochastic computational technique with a powerful global search capability for resolving the following form of the optimization problem, without constraints:
min f ( x ) , over x R n and f ( x ) : R n R
where the linear or nonlinear objective function f : R n R is minimized with respect to the vector of design variables x R n . Therefore, a particle in the PSO algorithm represents a potential solution x in (1). Let D denote the limited sub-region of an entire n-dimensional Euclidean space and be assumed to contain the optimal solutions. The conventional evolutionary PSO searching mechanism is initiated with a swarm randomly generated over the space D . Thereafter, each particle moves in a coordinated way through the n-dimensional search space D . The behavior of each particle is mainly influenced by both its own best previous experience and a social compulsion to move toward a single best particle among the entire swarm, as follows:
x i k + 1 = x i k + v i k + 1 ,
v i k + 1 = c 0 v i k + c 1 r 1 , i k ( x pbest , i x i k ) + c 2 r 2 , i k ( x gbest x i k )
where i = 1 , 2 , , n p denotes the index of the particle where n p is the number of particles in the swarm; k = 1 , 2 , , k max represents the iteration number where k max is the maximum number of allowable iterations; x i k R n and v i k R n denotes the position and velocity vectors, respectively, for the ith particle at the k-th iteration; x pbest , i denotes the vector of the best previous position yielding the minimum fitness value f ( · ) for the ith particle; x gbest denotes the vector of the global best position found by the entire swarm; c 0 denotes the inertia weight and c 1 and c 2 represent the cognitive and social scaling parameters, respectively; and r 1 , i k and r 2 , i k are random parameters generated uniformly in the range of [ 0 , 1 ] .
However, the above classical PSO method frequently faces the following two obstacles. First, for engineering optimal design problems that include several physical constraints, it necessarily introduces a technique that is capable of handling the problem (1), subject to
F : = { x R n | h 1 ( x ) 0 , h 2 ( x ) 0 , h 3 ( x ) 0 , , h m ( x ) 0 } ,
where h ( x ) : R n R denotes the constraint function and a feasible region denoted by F represents the set of design variables that satisfy the given constraint conditions. It should be noted that h ( x ) 0 , = 1 , 2 , , m in (4) is generally treated as a hard constraint in most engineering design problems. It is therefore critical to handle such constraint functions within the framework of PSO for evaluating the fitness of each particle over the PSO iteration and eventually obtaining the resultant optimal vector x * . Second, in practical engineering design work, some or all optimization problem design variables very commonly have integer/discrete values, as
x : = ( x I , x D , x C ) T = ( x 1 I , x 2 I , , x n I I , x 1 D , x 2 D , , x n D D , x 1 C , x 2 C , , x n C C ) T
where the number of design variables becomes n : = n I + n D + n C . Let the feasible subsets of integers and discrete and continuous design variables, be denoted by R n I , R n D and R n C . Using this information, researchers have tried to modify the original PSO scheme to handle constraint conditions, as well as integers and discrete variables [6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21]. In the following section, we propose an easy-to-use PSO scheme with a diversity classifier for cyclically neighboring subpopulations, which allows superior reliability and validity for constrained MIDC optimization problems.

3. PSO with a Diversity Classifier for Cyclically Neighboring Subpopulations

To improve the behavior performance of the swarm via a diversity-boosting mechanism, [23,27] studied a diversity-guided PSO (ATRE-PSO) strategy. Comparing to (3), the sole velocity update law for each particle is modified so that its velocity behavior follows one of the three-phase update laws, where the attraction, in-between, and repulsion phases are classified by the diversity measure DIV ( · ) per iteration, as follows:
v i k + 1 = c 0 v i k + c 1 r 1 , i k ( x pbest , i x i k ) + c 2 r 2 , i k ( x gbest x i k ) if DIV ( k ) > D high ; c 0 v i k c 1 r 1 , i k ( x pbest , i x i k ) c 2 r 2 , i k ( x gbest x i k ) if DIV ( k ) < D low ; c 0 v i k + c 1 r 1 , i k ( x pbest , i x i k ) c 2 r 2 , i k ( x gbest x i k ) otherwise ,
where D high and D low are chosen by the designer. The diversity classifier for each iteration is obtained by
DIV ( k ) = 1 n p i = 1 n p j = 1 n x i , j k x ^ j k 2 , x ^ j k : = 1 n p = 1 n p x , j k .
where x i , j denotes the jth entry of x i . They verified that their diversity-boosting mechanism allows particles to maintain an acceptable level of diversity for various benchmark problems with continuous variables and no constraint functions and thus helps to address the premature convergence problem. However, the above method may not reliably provide a level of population diversity that is sufficiently high to produce offspring that outperform their parents in the case of MIDC nonlinear programming with multimodal characteristics, even when no constraint is involved. Therefore, the PSO scheme that exploits the constricted social learning of each particle through cyclically neighboring subpopulations is presented in the following and such a scheme is extended for combination with the three-phase velocity update law to achieve diversity improvement for each subpopulation.
To consider the constraint functions of the optimization problem within the framework of our PSO algorithm, an efficient constraint handling mechanism introducing a so-called virtual objective function, f v ( x ) : R n R , is summarized below [28]. Note that any function can be such a virtual objective function if it guarantees that: (i) for any x R n satisfying all constraint conditions, f v ( x ) < 0 holds; and (ii) for any x i R n and x j R n , satisfying f ( x i ) < f ( x j ) , f v ( x i ) < f v ( x j ) holds. For example, the function arctan { f ( x ) } π / 2 can be a virtual function satisfying the above two properties (refer to [10,28] for details). Then, the given optimization problem (1) subject to (4) can be transformed into the following unconstrained optimization problem:
min x R n L ( x ) , L ( x ) : = h max ( x ) , if h max ( x ) 0 ; f v ( x ) , otherwise ,
where h max ( x ) is defined as h max ( x ) : = max h 1 ( x ) , h 2 ( x ) , , h m ( x ) . This modified objective function forces the swarm to keep moving into and remaining within the bounds of feasible space F during the initial iteration steps and then moving through the solution space to search for the optimum. Therefore, this constraint handling strategy enables one to obtain feasible design variables that, at least, satisfy the constraint conditions in a considerably more simple manner.
Next, to solve the unconstrained problem (8), the PSO scheme exploiting the constricted social learning of each particle through cyclically neighboring subpopulations is extended to incorporate the three-phase velocity behavior law (2) implemented by following the locally distributed diversity measures categorized for each subpopulation. In the canonical PSO strategy, the offspring of each particle is produced based on the information on its best previous position vector x pbest , i and the information on the current swarm best position vector x gbest , as shown in (3). This means that because x gbest : = arg min x { x pbest , i | i = 1 , , n p } L ( x ) , (i.e., x gbest is determined from all x pbest , i available at each evolutionary stage), x pbest , i may be mainly responsible for premature convergence [29]. Considering this viewpoint, the swarm flexibility is boosted by constructing cyclically neighboring multi-subpopulations, one for each particle (host particle), which results in the replacement of only one x gbest by different social-best positions for each subpopulation. Such a social-best position vector for the i-th particle x i at the k-th iteration, denoted by x sbest , i k , is determined only from its individual subpopulation, as follows:
x sbest , i k : = arg min x { x k | = i n s 2 , , i + n s 2 } L ( x )
where an even number of n s ( n p ) denotes the total number of neighboring particles within the subpopulation of x i , and x k : = x ( 1 mod n p ) + 1 k for < 1 or n p + 1 . This formulation verifies that the host particle x i interacts only with a limited number of particles belonging to the same subpopulation and those nearby neighbor particles transfer their positional information and fitness values to the host particle. Furthermore, the host particle of a certain subpopulation can be a nearby neighbor of some other subpopulation, which means that any information on the good fitness value of one host particle is cyclically propagated to other neighboring host particles during the evolution process. Therefore, this structure forces each particle to achieve social learning from its respective neighborhoods, not directly from the best position of the entire swarm, resulting in the enhancement of population diversity. It then follows that by replacing x gbest in (6) with x sbest , i k in (9) that
v i k + 1 = c 0 v i k ± c 1 r 1 , i k ( x pbest , i x i k ) ± c 2 r 2 , i k ( x sbest , i k x i k )
where the positive and negative signs are assigned according to the diversity measure DIV ( k ) at each iteration. This formulation implies that v i k enforces the accelerated movement of v i k toward x pbest , i and x sbest , i k . Conversely, the existing PSO method proposed by [23,27] governs the swarm behavior of all particles by following the sole diversity measure, DIV ( k ) in (7), derived from the position data of the entire population. Therefore, such a diversity classifier mechanism should be modified to be compatible with the swarm structure of cyclically neighboring multi-subpopulations. This can be achieved by developing multiple distributed diversity measures categorized for each local subpopulation as follows: for i = 1 , 2 , , n p ,
DIV i ( k ) = 1 n s + 1 = i n s 2 i + n s 2 j = 1 n x , j k x ^ j k 2 , x ^ j k : = 1 n s + 1 l = i n s 2 i + n s 2 x l , j k ,
where DIV i ( k ) refers to the diversity measure of x i k updated in each iteration. For n p = 9 and n s = 4 , the overall schematic diagram of our PSO with diversity classifiers applied respectively to cyclically neighboring multi-subpopulations is illustrated in Figure 1. In this figure, the first particle has DIV 1 ( k ) > D high and thus moves toward both pbest and sbest. In contrast, since DIV 2 ( k ) < D low , the second particle moves in the direction opposite to that of pbest and sbest. The ninth particle with D low DIV 9 ( k ) D high moves toward pbest but in the direction opposite to that of sbest.
Finally, a straightforward and reliable method for handling the MIDC design variables in (5) by applying the aforementioned PSO strategy in (2) and (10) with (11) is described. To achieve this objective, we introduce the widely used rounding-off technique, as follows.
(I)
Initialize the entire swarm. Then, the particles evolve in the D -dimensional search space, regardless of the variable types in x i ( : = ( x I , x D , x C ) T ) , of (5), via (2) and (10) with x sbest , i k . As a result, all elements (i.e., all design variables) in x i take floating-point values during the current stage.
(II)
For j = 1 , 2 , , n I , the variable x i , j , which is the jth entry of x i R , denotes the integer design variable x i , j I and thus must take an integer value. Let INT ( x i , j ) denote the nearest integer of x i , j R . i.e., x i , j is rounded to its nearest integer. It then follows that x i , j I INT ( x i , j ) is performed. Similarly, the discrete design variable x i , n I + D for = 1 , 2 , , n D takes the value of DIS ( x i , n I + ) as x i , n I + D DIS ( x i , n I + ) . Here, DIS ( x i , n I + ) indicates the nearest discrete value that x i , n I + takes in the given data set of discrete design values.
(III)
Fitness of individual particle with positional vector ( x i I , x i D , x i C ) , is evaluated based on the objective function L ( x ) , as in (8) subject to no constraint functions. Finally, x pbest , i and x sbest , i k in (9), are calculated.
For finding an optimal design variable vector, the aforementioned procedure for a set of n p particles keeps running until a certain termination criterion is met. Then, after a significant number of evolutions, optimal variables or approximate optimal variables are expected to be found. For handling constrained optimization problems with the MIDC variables, the flowchart of the overall PSO with a diversity classifier for cyclically neighboring subpopulations is illustrated in Figure 2, where the optimal design variable vector x * is determined as
x * : = arg min x { x i j | i = 1 , 2 , , n p ; j = 1 , 2 , , k } L ( x ) .
Conversely, the presented flowchart includes an optional mechanism such that, after a predetermined number of evolutions, k DIV , the velocity-related behavior of the particle is governed by the formula of (3), with x sbest , i k instead of x gbest . This mechanism only emphasizes the fine-searching ability via an attraction phase among particles belonging to the same subpopulation in the current search space that is sufficiently narrowed via the movement of the particles over the previous evolutions with the three-phase velocity update law.

4. Numerical Experimentation

This section introduces a wide variety of benchmark design problems that involve different types of MIDC variables and are subject to real-life design constraints for evaluating the performance of the proposed optimization algorithm. The 30 independent runs of the optimization task for each design problem were executed, where the user parameters were set as c 0 = 0.7298 and c 1 = c 2 = 1.4962 . The maximum PSO iteration number was set as k max = 1500 , for the car side impact design problem, and k max = 500 for the others. The swarm size was n p = 50 for the design problem of the reinforced concrete beam and n p = 100 for the others. To verify the practicality and reliability of the proposed PSO-based design method, the statistical results were primarily compared with several existing heuristic techniques [1,2,3,30], in which the authors extensively survey several evolutionary- and nature-inspired algorithms. Conversely, the relation between the neighborhood size, n s , and the robustness and preciseness of the PSO algorithm is studied in-depth, and one possible promising indicator for determining n s is then discussed.

4.1. Optimal Design of Pressure Vessel

The design problem for a compressed air storage tank in Figure 3, where the cylindrical pressure vessels are capped at both ends by two hemispherical heads, is examined. Its minimum volume is 750 ( f t 3 ), and the working pressure is 3000 (psi). The design objective is to minimize the total manufacturing cost, which comprises a combination of material, forming, and welding costs. Four kinds of design variables are defined as (1) the inner radius of the shell, r; (2) the length of the cylindrical section, l; (3) the thickness of the head, t h ; (4) the thickness of the shell, t s . Let x : = ( x 1 , x 2 , x 3 , x 4 ) T be defined as x = ( l , r , t s , t h ) T . Here, l and r are continuous variables but the available thicknesses for t h and t s are integer multiples of 0.0625 . Then, the optimal design problem is formulated as
min x f ( x ) : = 0.6224 x 3 x 2 x 1 + 1.7781 x 4 x 2 2 + 3.1661 x 3 2 x 1 + 19.84 x 3 2 x 2
subject to the constraint conditions, which are in accordance with the ASME design codes, such as
h 1 ( x ) : = 0.0193 x 2 x 3 0
h 2 ( x ) : = 0.00954 x 2 x 4 0
h 3 ( x ) : = 750 × 1728 π x 2 2 x 1 4 3 π x 2 3 0 , ,
h 4 ( x ) : = x 1 240 0 ,
where h 3 ( x ) denotes the constraint function concerning the minimum volume of 750 ( f t 3 ). The lower bound of the design variable x 1 is fixed at 20 (in) whereas its upper bound is set at 240 (in), accounting for the constraint h 4 ( x ) . Note that this differs from the limit of 200 (in), as described in some articles [1,2]. When the limit of x 1 has been determined, the ranges of x 2 , x 3 , and x 4 are automatically derived through the constraint conditions as: x 2 [ 37.7 , 63 ] , x 3 [ 0.6875 , 1.25 ] , and x 4 [ 0.3125 , 0.625 ] .
The optimizing procedure presented in Figure 2 is executed to find the optimal decision variable vector x * that achieves the minimum cost function f ( x ) and simultaneously satisfies the four constraints on h i ( x ) . The modified unconstrained problem (8) is defined with f v ( x ) : = arctan { f ( x ) } π / 2 , and ( D low , D high , k DIV ) = ( 10 , 20 , 150 ) . The best solution and corresponding optimal objective function value are compared with several evolutionary and nature-inspired algorithms in Table 1. Among the compared approaches, a feasible best solution was obtained only by [1] and our PSO scheme. However, the worst, mean, and standard deviation values of the statistical results demonstrate that our scheme far surpassed the other reported optimization methods, including [1]. Despite this, the result of [31] seems to outperformours in terms of optimal cost value. However, their optimal design variables are not feasible as the constraint condition on h 3 ( x ) is violated; this is also the case in the method of [32]. Finally, the statistical result of our PSO method was obtained when the total number of fitness evaluations was 50,000 (= n p × k max ), which is twice the number of those conducted in [1]. However, even in the case of 25,000 function evaluations, with n p = 50 and k max = 500, our approach, which produced the best, worst, average, and standard deviation of 5850.3830603, 5850.4168780, 5850.3842066, and 0.0061715, respectively, outperformed the results given in [1].

4.2. Optimal Design of Reinforced Concrete Beam

The configuration of a reinforced concrete beam in Figure 4, which has a simply supported span of 30 (ft) subject to a live load of 2.0 (klbf) and a dead load of 1.0 (klbf). The yield stress of the reinforced steel σ y and the concrete compressive strength σ c are 50 (ksi) and 5 (ksi), respectively. The unit cost of steel and concrete are, respectively, 1.0 and 0.02 $ / in 2 per linear ft. The design objective is to achieve the minimum manufacturing cost, where the three design variables are: (1) the width b of the concrete beam, which is an integer variable; (2) the depth h of the concrete beam, which is a continuous variable; and (3) the cross-sectional area A s of the reinforcing bar. Here, A s is a discrete type variable that is one of the standardized dimensions listed in Table 2.
Let the design variable vector x : = ( x 1 , x 2 , x 3 ) T be defined as x = ( A s , b , h ) T . Then, the design problem is formulated as:
min x f ( x ) : = 29.4 x 1 + 0.6 x 2 x 3 .
The constraint conditions are stated as follows. The first constraint is for the width-to-depth ratio of the beam, which should be restricted to below 4, as
h 1 ( x ) : = x 2 x 3 4 0 .
The second constraint condition for guaranteeing the safety requirement defined in the ACC building code 318–77 is given as
M u = 0.9 x 1 σ y ( 0.8 x 2 ) 1.0 0.59 x 1 σ y 0.8 x 2 x 3 σ c 1.4 M d + 1.7 M l ,
where M u , M l , and M d denote the flexural strength, live load, and dead load moments of the beam, respectively. Let M d = 1350 (kip-in) and M l = 2700 (kip-in). The above condition is then reformulated as shown in [1]:
h 2 ( x ) : = 180 + 7.375 x 1 2 x 3 x 1 x 2 0 .
The bounds of the design variables are x 2 { 28 , 29 , , 40 } in, and 5 x 3 10 in. We set the user parameters as ( D low , D high , k DIV ) = ( 0.1 , 0.3 , 200 ) .
Table 3 compares the best solution and the statistical results for the proposed and other reported methods. This table shows that [1] also found the best solution but our PSO scheme offered the best objective function value, with a smaller number of function evaluations. Furthermore, the given statistical results show that the standard deviation on the optimized costs of 30 independent runs using our PSO method equals zero, demonstrating the remarkable optimization accuracy and reliability.

4.3. Optimal Design of Helical Compression Spring

The helical compression spring of Figure 5 is designed to support an axially-guided constant compressive load. This must be optimized to support a certain load without failure, with a minimum wire volume (minimum weight). The integer-valued number of spring coils (N), real-valued outside diameter of the spring (D), and spring wire diameter (d), which takes one of the allowable 42 discrete values given in Table 4, are included as optimization variables and are set as x = ( x 1 , x 2 , x 3 ) T : = ( D , N , d ) T . Accordingly, the objective function related to the minimization of the spring volume is formulated as
min x f ( x ) : = 1 4 π 2 x 1 x 3 2 ( x 2 + 2 ) .
The following constraint conditions specify the design restrictions:
h 1 ( x ) : = ( 8 c K F max ) / ( π x 3 2 ) S 0 , K : = ( 4 c 1 ) / ( 4 c 4 ) + 0.615 / c , c : = x 1 / x 3
h 2 ( x ) : = l l max 0
h 3 ( x ) : = d min x 3 0
h 4 ( x ) : = ( x 1 + x 3 ) D max 0
h 5 ( x ) : = 3.0 c 0 ,
h 6 ( x ) : = δ p δ p m 0 , δ p : = F p / k , k : = ( G x 3 4 ) / ( 8 x 2 c 3 )
h 7 ( x ) : = δ p + ( F max F p ) / k + 1.05 ( x 2 + 2 ) x 3 l 0 , l : = F max / k + 1.05 ( x 2 + 2 ) x 3
h 8 ( x ) : = δ w ( F max F p ) / k 0
The supplied numerical data included in the problem statement are listed in Table 5, and the constraint details can be found in [1]. Combining the constraints h 3 ( x ) - h 5 ( x ) , the boundaries of x 1 and x 2 become [ 0.6 , 3.0 ] and [ 1 , 70 ] , respectively, and the limits of x 3 are set automatically as [ 0.0090 , 0.5000 ] from Table 4. We set the user parameters as ( D low , D high , k DIV ) = ( 0.1 , 0.4 , 250 ) .
This design problem was previously examined using various optimization algorithms including FA [1], PSO [2,14], DE [21], GA [38] and HSIA (hybrid swarm intelligence approach) [39]. The reported optimization performances, including statistical results, were compared in-depth with our PSO scheme, as presented in Table 6. In this table, it can be seen that our PSO scheme detected the best solution, satisfying all constraints; this was also reported in [2,14,21]. Among these, [14] found the best solution with a smaller number of function evaluations. However, in terms of statistical results, our PSO algorithm surpassed other conventional optimization methods with a reasonable number of function evaluations.

4.4. Optimal Design of Belleville Spring

This benchmark example is to design a Belleville spring of minimum weight, as shown in Figure 6. The four design variables are as follows: (1) thickness t ( = : x 1 ), (2) height h ( = : x 2 ), (3) internal diameter D i ( = : x 3 ), and (4) external diameter D e ( = : x 4 ). The cost function is given as:
min x f ( x ) : = 0.07075 π ( x 4 2 x 3 2 ) x 1 .
The constraints related to compressive stress, deflection, height to deflection, height to maximum height, outer diameter, inner diameter, and slope are formulated as
h 1 ( x ) : = 1000 S 4 E δ max ( 1 μ 2 ) α x 4 2 β x 2 δ max 2 + γ x 1 0 ,
h 2 ( x ) : = 4 E δ max ( 1 μ 2 ) α x 4 2 x 2 δ 2 x 2 δ x 1 + x 1 3 δ = δ max P max 0 ,
h 3 ( x ) : = δ l δ max 0 ,
h 4 ( x ) : = H x 2 x 1 0 ,
h 5 ( x ) : = D max x 4 0 ,
h 6 ( x ) : = x 4 x 3 0 ,
h 7 ( x ) : = 0.3 x 2 x 4 x 3 0 ,
where
α = 6 π ln K K 1 K 2 , β = 6 π ln K K 1 ln K 1 , γ = 6 π ln K K 1 2 , P max = 5400 lb , δ max = 0.2 in , μ = 0.3 , S = 200 KPsi , H = 2 in , D max = 12.01 in , K = x 4 / x 3 , δ l = g ( a ) x 2 , a = x 2 / x 1 .
The variation in g ( a ) according to a ( = x 2 / x 1 ) can be determined as shown in Table 7. We set the user parameters as ( D low , D high , k DIV ) = ( 0.25 , 0.35 , 100 ) .
A comparison of the best solution and statistical results for this design problem is presented in Table 8. It can be verified from this table that although three conventional methods (MBA, ABC, and TLBO) also found the optimal solution, the remarkable superiority of our PSO scheme to these optimizers is demonstrated by the worst, mean, and standard deviation values. It should be noted that the results of [41] seem to outperform ours in terms of optimal cost value; however, the solution was not feasible because the second constraint ( h 2 ( x ) ) was violated.

4.5. Optimal Design of Speed Reducer

The design objective is to minimize the weight of the speed reducer of Figure 7, which is subject to constraints, including restrictions concerning the bending and surface stresses of the gear teeth, transverse deflection of the two shafts caused by the transmitted force, and stresses in the shafts. The decision variables, x 1 to x 7 , represent, respectively, b (face width), m (tooth module), n (number of teeth on the pinion), l 1 (length of shaft 1 between bearings), l 2 (length of shaft 2 between bearings), d 1 (shaft diameter 1), and d 2 (diameter of shaft 2). Here, x 3 should take an integer value, while all the other variables are continuous.
The cost function then becomes
min x f ( x ) : = 0.7854 x 1 x 2 2 ( 3.3333 x 3 2 + 14.9334 x 3 43.0934 ) 1.508 x 1 ( x 6 2 + x 7 2 ) + 7.4777 ( x 6 3 + x 7 3 ) + 0.7854 ( x 4 x 6 2 + x 5 x 7 2 ) .
The following 11 constraint conditions result in the problem being highly complex:
h 1 ( x ) : = 27 x 1 1 x 2 2 x 3 1 1 0 ,
h 2 ( x ) : = 397.5 x 1 1 x 2 2 x 3 2 1 0 ,
h 3 ( x ) : = 1.93 x 2 1 x 3 1 x 4 3 x 6 4 1 0 ,
h 4 ( x ) : = 1.93 x 2 1 x 3 1 x 5 3 x 7 4 1 0 ,
h 5 ( x ) : = 745 x 4 x 2 x 3 2 + 16.9 × 10 6 0.1 x 6 3 1100 0 ,
h 6 ( x ) : = 745 x 5 x 2 x 3 2 + 157.5 × 10 6 0.1 x 7 3 850 0 ,
h 7 ( x ) : = x 2 x 3 40 0 ,
h 8 ( x ) : = 5 x 1 / x 2 0 ,
h 9 ( x ) : = x 1 / x 2 12 0 ,
h 10 ( x ) : = ( 1.5 x 6 + 1.9 ) x 4 1 1 0 ,
h 11 ( x ) : = ( 1.1 x 7 + 1.9 ) x 5 1 1 0 ,
where the limits of the variables are set as 2.6 x 1 3.6 , 0.7 x 2 0.8 , 17 x 3 28 , 7.3 x 4 8.3 , 7.3 x 5 8.3 , 2.9 x 6 3.9 , and 5.0 x 7 5.5 . For this minimization problem, the performance of our PSO algorithm is compared with the extensive survey results presented in some articles [3,44,45] of the following six optimization methods: (i) DEDS: differential evolution with dynamic stochastic selection (ii) DELC: differential evolution with level comparison (iii) HEAA: hybrid evolutionary algorithm and adaptive constraint handling technique (iv) MDE: modified differential evolution (v) PSO-DE: PSO with differential evolution and (vi) MBA: mine blast algorithm (see [3] for details.). The user parameters were set as ( D low , D high , k DIV ) = ( 0.25 , 0.35 , 100 ) .
The optimization results are presented in Table 9, which demonstrates that our PSO scheme outperformed some other optimization techniques in terms of precision (i.e., the lowest cost function value). For the detected best solution, the constraint values become h 1 ( x * ) = 0.073915 , h 2 ( x * ) = 0.197999 , h 3 ( x * ) = 0.107955 , h 4 ( x * ) = 0.904644 , h 5 ( x * ) = 534.994106 , h 6 ( x * ) = 0.000000 , h 7 ( x * ) = 28.100000 , h 8 ( x * ) = 0.000000 , h 9 ( x * ) = 7.000000 h 10 ( x * ) = 0.143836 , and h 11 ( x * ) = 0.000000 , which demonstrates that all the constraint conditions are satisfied. Furthermore, it is observed in Table 9 that the proposed optimizer yields better results than those reported in the previous studies for the best, mean, and worst solutions, demonstrating the remarkable consistency of our optimization method.

4.6. Optimal Design of Stepped Cantilever Beam

This example considers a stepped cantilever beam subject to the end load shown in Figure 8 where minimizing the volume of the beam is the design objective. The height and width of the cross-sectional area in all five segments are taken as the design variables, which results in { x 1 , x 3 , x 5 , x 7 , x 9 } : = { b 1 , b 2 , b 3 , b 4 , b 5 } and { x 2 , x 4 , x 6 , x 8 , x 10 } : = { h 1 , h 2 , h 3 , h 4 , h 5 } . Then, the objective function is formulated as:
min x f ( x ) : = 100 ( x 1 x 2 + x 3 x 4 + x 5 x 6 + x 7 x 8 + x 9 x 10 ) ,
where x 1 { 1 , 2 , 3 , 4 , 5 } , x 3 , x 5 { 2.4 , 2.6 , 2.8 , 3.1 } , x 2 , x 4 { 45 , 50 , 55 , 60 } , x 6 { 30 , 31 , , 65 } , 1 x 7 , x 9 5 , and 30 x 8 , x 10 65 . The constraint conditions are described as follows. First, the bending stress in all segments must be less than the allowable limit of 14,000 ( N / cm 2 ) :
h 1 ( x ) : = 600 P x 9 x 10 2 14 , 000 0 ,
h 2 ( x ) : = 1200 P x 7 x 8 2 14 , 000 0 ,
h 3 ( x ) : = 1800 P x 5 x 6 2 14 , 000 0 ,
h 4 ( x ) : = 2400 P x 3 x 4 2 14 , 000 0 ,
h 5 ( x ) : = 3000 P x 1 x 2 2 14 , 000 0 , ,
where the concentrated load P = 50,000 (N). Second, the tip deflection of the cantilever beam must be smaller than the allowable limit of 2.7 ( cm ):
h 6 ( x ) : = 100 3 P 3 E 1 I 5 + 7 I 4 + 19 I 3 + 37 I 2 + 61 I 1 2.7 0 ,
where the elastic modulus of the material is E = 2 × 10 7 ( N / cm 2 ). Finally, the aspect ratio between the height and width of the cross section in each segment must be less than a value of 20:
h 7 ( x ) : = x 10 x 9 20 0 ,
h 8 ( x ) : = x 8 x 7 20 0 ,
h 9 ( x ) : = x 6 x 5 20 0 ,
h 10 ( x ) : = x 4 x 3 20 0 ,
h 11 ( x ) : = x 2 x 1 20 0 , .
The user parameters of our PSO method were set as ( D low , D high , k DIV ) = ( 5 , 10 , 100 ) .
The best design variables and statistical results obtained via eight optimization methods are compared in Table 10. This demonstrates that the proposed PSO scheme surpasses all other reported optimizers for finding the lowest cost function value. Furthermore, the mean and worst values achieved via our optimizer are considerably better than those reported in other studies. It is noteworthy that the zero standard deviation demonstrates the remarkable search precision and reliability of our optimization engine.

4.7. Optimal Design of Rolling Element Bearing

The design objective is to increase, as much as possible, the dynamic load-carrying capacity of the rolling element bearing in Figure 9. Five design variables are defined as follows; ball diameter ( D b = : x 1 ), pitch diameter ( D m = : x 2 ), number of balls ( Z = : x 3 ), and inner and outer raceway curvature coefficients ( f i = : x 4 and f o = : x 5 ). Note that other variables, such as K D min ( = : x 6 ), K D max ( = : x 7 ), ε ( = : x 8 ), e ( = : x 9 ), and ζ ( = : x 10 ), appear only in the constraint conditions. Here, the number of balls x 3 should take an integer value whereas the remainders are taken as continuous values. The optimization problem is stated as
max x f ( x ) : = C d = f c x 5 2 / 3 x 1 1.8 , if x 1 25.4 ( mm ) , C d = 3.647 f c x 5 2 / 3 x 1 1.4 , if x 1 > 25.4 ( mm ) ,
subject to the following constraints that are derived from kinematic and manufacturing considerations:
h 1 ( x ) : = ϕ 0 2 sin 1 ( x 1 / x 2 ) x 5 + 1 0 ,
h 2 ( x ) : = 2 x 1 x 6 ( D d ) 0 ,
h 3 ( x ) : = x 7 ( D d ) 2 x 1 0 ,
h 4 ( x ) : = x 10 B w x 1 0 ,
h 5 ( x ) : = x 2 0.5 ( D + d ) 0 ,
h 6 ( x ) : = ( 0.5 + x 9 ) ( D + d ) x 2 0 ,
h 7 ( x ) : = 0.5 ( D x 2 x 1 ) x 8 x 1 0 ,
h 8 ( x ) : = x 3 0.515 0 ,
h 9 ( x ) : = x 4 0.515 0 ,
where
f c = 37.91 1 + 1.04 1 γ 1 + γ 1.72 x 3 ( 2 x 4 1 ) x 4 ( 2 x 3 1 ) 0.41 10 / 3 0.3 × γ 0.3 ( 1 γ ) 1.39 ( 1 + γ ) 1 / 3 × 2 x 3 2 x 3 1 0.41 , ϕ 0 = 2 π 2 cos 1 [ { ( D d ) / 2 3 ( T / 4 ) } 2 + { D / 2 T / 4 x 1 } 2 { d / 2 + T / 4 } 2 ] 2 { ( D d ) / 2 3 ( T / 4 ) } { D / 2 T / 4 x 1 } , γ = x 1 x 2 , x 3 = r i x 1 , x 4 = r o x 1 , T = D d 2 x 1 , D = 160 , d = 90 , B w = 30 , r i = r o = 11.033 , 0.5 ( D + d ) x 2 0.6 ( D + d ) , 0.15 ( D d ) x 1 0.45 ( D d ) , 4 x 5 50 , 0.515 x 3 , x 4 0.6 , 0.4 x 6 0.5 , 0.6 x 7 0.7 , 0.3 x 8 0.4 , 0.02 x 9 0.1 , 0.6 x 10 0.85 .
The simulation results using our PSO scheme are presented and compared with the results of other conventional methods in Table 11. The user parameters of the PSO method were set as ( D low , D high , k DIV ) = ( 0.3 , 0.5 , 100 ) . Although the reported best objective function value of TLBO in [30] is identical to ours, their solution does not satisfy the given constraint conditions. It should be noted that the constraint function values given in [30] cannot be obtained from their optimal solution. However, [3] reported the best cost function value of 85535.9611 but their optimal solution provides f ( x * ) = 81843.68625 , which is a lower value than ours. Furthermore, the constraint function value presented by [3] demonstrates that their method produced an infeasible solution. Among the compared studies, a feasible solution was obtained by [52]; however, our optimization scheme surpassed their method, having remarkable precision and an acceptable number of function evaluations.

4.8. Car Side Impact Design Problem

This problem addresses the case of an automobile side impact. Its objective is to improve the side impact crash performance while reducing the total weight of the automobile. Figure 10 illustrates the finite element model of the automobile, which is adopted from [1]. This problem contains eleven variables, such as the thickness of the B-Pillar inner ( x 1 ), B-Pillar reinforcement ( x 2 ), floor side inner ( x 3 ), cross members ( x 4 ), door beam ( x 5 ), door belt line reinforcement ( x 6 ), roof rail ( x 7 ), materials of B-pillar inner and floor side inner ( x 8 and x 9 ), and barrier height and hitting position ( x 10 and x 11 ). Nine design variables are continuous ( x , = 1 , 2 , , 7 , 10 , 11 ), and x 8 and x 9 are discrete. This optimal design problem is formulated as
min x f ( x ) : = 1.98 + 4.90 x 1 + 6.67 x 2 + 6.98 x 3 + 4.01 x 4 + 1.78 x 5 + 2.73 x 7 ,
subject to
h 1 ( x ) = F a ( load in the abdomen ) 1 ( kN ) ,
h 2 ( x ) : = V C u ( dummy upper chest ) 0.32 ( m / s ) ,
h 3 ( x ) : = V C m ( dummy middle chest ) 0.32 ( m / s ) ,
h 4 ( x ) : = V C l ( dummy lower chest ) 0.32 ( m / s ) ,
h 5 ( x ) : = Δ u r ( upper rib deflection ) 32 ( mm ) ,
h 6 ( x ) : = Δ m r ( middle rib deflection ) 32 ( mm ) ,
h 7 ( x ) : = Δ l r ( lower rib deflection ) 32 ( mm ) ,
h 8 ( x ) : = F p ( public force ) 4 ( kN ) ,
h 9 ( x ) : = V M B P ( velocity of V-Pillar at the middle point ) 9.9 ( mm / ms ) ,
h 10 ( x ) : = V F D ( velocity of front door at V-Pillar ) 15.7 ( mm / ms ) ,
where
F a = 1.16 0.3717 x 2 x 4 0.00931 x 2 x 10 0.484 x 3 x 9 + 0.01343 x 6 x 10
V C u = 0.261 0.0159 x 1 x 2 0.188 x 1 x 8 0.019 x 2 x 7 + 0.0144 x 3 x 5 + 0.0008757 x 5 x 10 + 0.08045 x 6 x 9 + 0.00139 x 8 x 11 + 0.00001575 x 10 x 11
V C m = 0.214 + 0.00817 x 5 0.131 x 1 x 8 0.0704 x 1 x 9 + 0.03099 x 2 x 6 0.018 x 2 x 7 + 0.0208 x 3 x 8 + 0.121 x 3 x 9 0.00364 x 5 x 6 + 0.0007715 x 5 x 10 0.0005354 x 6 x 10 + 0.00121 x 8 x 11 + 0.00184 x 9 x 10 0.02 x 2 2
V C l = 0.74 0.61 x 2 0.163 x 3 x 8 + 0.001232 x 3 x 10 0.166 x 7 x 9 + 0.227 x 2 2
Δ u r = 28.98 + 3.818 x 3 4.2 x 1 x 2 + 0.0207 x 5 x 10 + 6.63 x 6 x 9 7.7 x 7 x 8 + 0.32 x 9 x 10
Δ m r = 33.86 + 2.95 x 3 + 0.1792 x 10 5.057 x 1 x 2 11.0 x 2 x 8 0.0215 x 5 x 10 9.98 x 7 x 8 + 22.0 x 8 x 9
Δ l r = 46.36 9.9 x 2 12.9 x 1 x 8 + 0.1107 x 3 x 10
F p = 4.72 0.5 x 4 0.19 x 2 x 3 0.0122 x 4 x 10 + 0.009325 x 6 x 10 + 0.000191 x 11 2
V M B P = 10.58 0.674 x 1 x 2 1.95 x 2 x 8 + 0.02054 x 3 x 10 0.0198 x 4 x 10 + 0.028 x 6 x 10
V F D = 16.45 0.489 x 3 x 7 0.843 x 5 x 6 + 0.0432 x 9 x 10 0.0556 x 9 x 11 0.000786 x 11 2 .
The simple bounds of the design variables are 0.5 x 1 , x 2 , x 3 , x 4 , x 5 , x 6 , x 7 1.5 , 30 x 10 , x 11 30 , and two discrete type variables are x 8 , x 9 { 0.192 , 0.345 } .
Table 12 compares the best solution and the statistical results of the proposed method and other reported methods. The user parameters of our PSO method were set as ( D low , D high , k DIV ) = ( 0.01 , 0.03 , 200 ) . This table shows that our optimization scheme produced the lowest cost function value together with a highly reliable performance, which is confirmed from the worst, mean, and standard deviation values. In particular, the superiority of our optimizer to the recently developed FA method in [1] is in the objective function values and constrained accuracy performance.

4.9. Discussions

The statistical results corresponding to the different percent ratio of n s to n p for the numerical examples in the previous sections are given in Figure 11. This figure verifies that a remarkable optimization accuracy and reliability can be achieved when the ratio n s / n p is approximately 10 % . However, the worst solution and standard deviation were most strongly dependent on the value of n s (see the cases of n s / n p 4 % and n s / n p 32 % ). The above findings can be explained as follows. As illustrated in Figure 1, each particle evolves based only on the local information of the neighbors. The amount of information becomes more limited as the ratio of n s to n p decreases. Therefore, the low n s / n p results in the tardy propagation of the successful social-best search result by a certain particle to particles belonging to other subpopulations. Such a case may lead to a slow convergence rate of the swarm toward an optimal solution while reducing the risk of the particles prematurely converging. This fact implies that more iterations could be required to guarantee the convergence of particles to a global optimum. In contrast, a high ratio of n s to n p implies that a particle obtains a large amount of information from several other neighboring particles. Therefore, a high convergence rate may be achieved despite the high risk of a fatal premature convergence phenomenon occurring, which usually causes the performance consistency of the optimizer to deteriorate. Figure 12 illustrates the convergence characteristics of f min k : = min f ( x 1 k ) , f ( x 2 k ) , , f ( x n p k ) for all 30 independent runs. Here, “Modified ATRE-PSO" refers to the method in which the original ATRE-PSO of [23,27] and our constraint handling scheme are combined, which is identical to the proposed PSO method with n s / n p = 100 % (i.e., star topology). The figure shows that the proposed cyclic network-based PSO method converged to the near-optimal solution in every run of the algorithm as compared with the star topology-based method. This improved reliability was achieved in all the considered design problems and thus can be considered a superiority of the optimization capability of the proposed method.

5. Conclusions

This study presented the realization of an easy-to-use diversity-guided PSO methodology for engineering optimal design problems with MIDC variables and subject to various real-life physical constraints. This scheme uses diversity classifiers for cyclically neighboring subpopulations and is characterized by superior reliability and validity. This PSO strategy combined two kinds of diversity-enhancing mechanisms. The first diversity-boosting tool was to construct cyclically neighboring multisubpopulations, one for each particle (host particle), where the host particle interacts only with a limited number of nearby neighbor particles during the search process. In this structure, information on the good fitness value of one host particle is cyclically propagated to other neighboring host particles during the evolution process. Therefore, this structure enforces each particle to achieve social learning from its respective neighborhoods and not directly from the best position of the entire swarm, which results in the enhancement of population diversity. The second diversity improvement for each subpopulation was to extend the above PSO mechanism to be combined with the idea of the three-phase velocity update law for governing the behavior of the host particle. In this extended PSO algorithm, the velocity behavior of each host particle per iteration is governed by either attraction, repulsion, or in-between phases, chosen according to the diversity classifier applied to the subpopulation of the host particle. In summary, our PSO scheme not only exploits the constricted social learning of each particle through cyclically neighboring subpopulations but also incorporates the three-phase velocity behavior law implemented by following the locally distributed diversity measures categorized for each subpopulation. This novel PSO scheme features high reliability, as well as superior practicality for engineering optimization problems involving MIDC variables, which are handled via the straightforward rounding-off technique that has been widely adopted in swarm-inspired metaheuristic search technologies. Several benchmark MIDC design problems were investigated in-depth, clearly verifying that the proposed PSO scheme is a highly useful tool, providing remarkable reliability.

Author Contributions

Conceptualization and methodology, M.C., S.S. and T.-H.K.; software and experiments, T.-H.K. and M.C.; validation and formal analysis, M.C. and S.S.; writing—original draft preparation, M.C.; writing—review and editing, T.-H.K.; funding acquisition, T.-H.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Basic Science Research Program through the National Research Foundation (NRF) of Korea funded by the Ministry of Education (NRF-2016R1D1A1B03935288).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gandomi, A.H.; Yang, X.-S.; Alavi, A.H. Mixed variable structural optimization using firefly algorithm. Comput. Struct. 2011, 89, 2325–2336. [Google Scholar] [CrossRef]
  2. Datta, D.; Figueira, J.R. A real-integer-discrete-coded particle swarm optimization for design problems. Appl. Soft Comput. 2011, 11, 3625–3633. [Google Scholar] [CrossRef]
  3. Sadollaha, A.; Bahreininejada, A.; Eskandarb, H.; Hamdia, M. Mine blast algorithm: A new population based algorithm for solving constrained engineering optimization problems. Appl. Soft Comput. 2013, 13, 2592–2612. [Google Scholar] [CrossRef]
  4. Banks, A.; Vincent, J.; Anyakoha, C. A review of particle swarm optimization. Part I: Background and development. Nat. Comput. 2007, 6, 467–484. [Google Scholar] [CrossRef]
  5. Banks, A.; Vincent, J.; Anyakoha, C. A review of particle swarm optimization. Part II: Hybridisation, combinatorial, multicriteria and constrained optimization, and indicative applications. Nat. Comput. 2007, 7, 109–124. [Google Scholar] [CrossRef]
  6. Tomassetti, G. A cost-effective algorithm for the solution of engineering problems with particle swarm optimization. Eng. Optim. 2010, 42, 471–495. [Google Scholar] [CrossRef] [Green Version]
  7. Pulido, G.T.; Coello, C.A.C. A constraint-handling mechanism for particle swarm optimization. In Proceedings of the 2004 Congress on Evolutionary Computation, Portland, OR, USA, 19–23 June 2004; pp. 1396–1403. [Google Scholar]
  8. Parsopoulos, K.E.; Vrahatis, M.N. Unified particle swarm optimization for solving constrained engineering optimization problems. In Advances in Natural Computation; Springer: Berlin/Heidelberg, Germany, 2005; pp. 582–591. [Google Scholar]
  9. Mezura-Montes, E.; Coello Coello, C.A. Constraint-handling in nature-inspired numerical optimization: Past, present and future. Swarm Evol. Comput. 2011, 1, 173–194. [Google Scholar] [CrossRef]
  10. Chun, S.; Kim, Y.-T.; Kim, T.-H. A diversity-enhanced constrained particle swarm optimizer for mixed integer-discrete-continuous engineering design problems. Adv. Mech. Eng. 2013, 2013, 130750-1–130750-11. [Google Scholar] [CrossRef] [Green Version]
  11. Das, S.; Suganthan, P.N. Differential evolution: A survey of the state-of-the-art. IEEE Trans. Evol. Comput. 2011, 15, 4–31. [Google Scholar] [CrossRef]
  12. Parsopoulos, K.E.; Vrahatisy, M.N. Recent approaches to global optimization problems through particle swarm optimization. Nat. Comput. 2002, 1, 235–306. [Google Scholar] [CrossRef]
  13. Coelho, L.S. An efficient particle swarm approach for mixed-integer programming in reliability-redundancy optimization applications. Reliab. Eng. Syst. Saf. 2009, 94, 830–837. [Google Scholar] [CrossRef]
  14. He, S.; Prempain, E.; Wu, Q.H. An improved particle swarm optimizer for mechanical design optimization problems. Eng. Optim. 2004, 36, 585–605. [Google Scholar] [CrossRef]
  15. Venter, G.; Sobieski, J.S. Particle swarm optimization. AIAA J. 2003, 41, 1583–1589. [Google Scholar] [CrossRef] [Green Version]
  16. Venter, G.; Sobieski, J.S. Multidisciplinary optimization of a transport aircraft wing using particle swarm optimization. Struct. Multidiscip. Optim. 2004, 26, 121–131. [Google Scholar] [CrossRef] [Green Version]
  17. Akhtar, S.; Tai, K.; Ray, T. A socio-behavioural simulation model for engineering design optimization. Eng. Optim. 2002, 34, 341–354. [Google Scholar] [CrossRef]
  18. Cagnina, L.C.; Esquivel, S.C.; Coello Coello, C.A. Solving engineering optimization problems with the simple constrained particle swarm optimizer. Informatica 2008, 32, 319–325. [Google Scholar]
  19. Coello, C.A.C. Use of a self-adaptive penalty approach for engineering optimization problems. Comput. Ind. 2000, 41, 113–127. [Google Scholar] [CrossRef]
  20. Hu, X.; Eberhart, R.C.; Shi, Y. Engineering optimization with particle swarm. In Proceedings of the 2003 the IEEE Swarm Intelligence Symposium, Indianapolis, IN, USA, 24–26 April 2003; pp. 53–57. [Google Scholar]
  21. Lampinen, J.; Zelinka, I. Mixed-integer-discrete-continuous optimization by differential evolution. Part 2. A practical example. In Fifth International Mendel Conference on Soft Computing (MENDEL’99); Brno University of Technology: Brno, Czech Republic, 1999; pp. 77–81. [Google Scholar]
  22. Liu, H.; Abraham, A.; Zhang, W. A fuzzy adaptive turbulent particle swarm optimization. Int. J. Innov. Comput. Appl. 2007, 1, 39–47. [Google Scholar] [CrossRef]
  23. Pant, M.; Thangaraj, R.; Abraham, A. Particle swarm optimization: Performance tuning and empirical analysis. In Foundations of Computational Intelligence Volume 3: Global Optimization; Springer: Berlin/Heidelberg, Germany, 2009; pp. 101–128. [Google Scholar]
  24. Lane, J.; Engelbrecht, A.; Gain, J. Particle swarm optimization with spatially meaningful neighbours. In Proceedings of the 2008 IEEE Swarm Intelligence Symposium, St. Louis, MO, USA, 21–23 September 2008; pp. 1–8. [Google Scholar]
  25. Richards, M.; Ventura, D. Dynamic sociometry in particle swarm optimization. In Proceedings of the Sixth International Conference on Computational Intelligence and Natural Computing, Cary, NC, USA, 26–30 September 2003; pp. 1557–1560. [Google Scholar]
  26. Mendes, R.; Neves, J. What makes a successful society? Experiments with population topologies in particle swarms. Lect. Notes Comput. Sci. 2004, 3171, 346–355. [Google Scholar]
  27. Pant, M.; Radha, T.; Singh, V.P. A simple diversity guided particle swarm optimization. In Proceedings of the 2007 IEEE Congress on Evolutionary Computation, Singapore, 25–28 September 2007; pp. 3294–3299. [Google Scholar]
  28. Maruta, I.; Kim, T.-H.; Song, D.H.; Sugie, T. Synthesis of fixed-structure robust controllers using a constrained particle swarm optimizer with cyclic neighborhood topology. Expert Syst. Appl. 2013, 40, 3595–3605. [Google Scholar] [CrossRef]
  29. Zavala, A.E.M.; Aguirre, A.H.; Diharce, E.V. Robust PSO-based constrained optimization by perturbing the particle’s memory. In Swarm Intelligence: Focus on Ant and Particle Swarm Optimization; I-Tech Education and Publishing: Chicago, IL, USA, 2007; pp. 57–76. [Google Scholar]
  30. Rao, R.V.; Savsani, V.J.; Vakharia, D.P. Teaching-learning-based optimization: A novel method for constrained mechanical design optimization problems. Comput.-Aided Des. 2011, 43, 303–315. [Google Scholar] [CrossRef]
  31. Mahdavi, M.; Fesanghary, M.; Damangir, E. An improved harmony search algorithm for solving optimization problems. Appl. Math. Comput. 2007, 188, 1567–1579. [Google Scholar] [CrossRef]
  32. Dimopoulos, G.G. Mixed-variable engineering optimization based on evolutionary and social metaphors. Comput. Methods Appl. Mech. Eng. 2007, 196, 803–817. [Google Scholar] [CrossRef]
  33. Hedar, A.R.; Fukushima, M. Derivative-free filter simulated annealing method for constrained continuous global optimization. J. Glob. Optim. 2005, 35, 521–649. [Google Scholar] [CrossRef]
  34. Amir, H.M.; Hasegawa, T. Nonlinear mixed-discrete structural optimization. J. Struct. Eng. 1898, 115, 626–646. [Google Scholar] [CrossRef]
  35. Shih, C.J.; Yang, Y.C. Generalized Hopfield network based structural optimization using sequential unconstrained minimization technique with additional penalty strategy. Adv. Eng. Softw. 2002, 33, 721–729. [Google Scholar] [CrossRef]
  36. Mezura-Montes, E.; Hernandez-Ocana, B. Modified bacterial foraging optimization for engineering design. In Intelligent Engineering Systems through Artificial Neural Networks; ASME Press (eBook): New York, NY, USA, 2009; Volume 19, pp. 357–364. [Google Scholar]
  37. Yun, Y.S. Study on Adaptive Hybrid Genetic Algorithm and Its Applications to Engineering Design Problems. Master’s Thesis, Waseda University, Tokyo, Japan, 2005. [Google Scholar]
  38. Wu, S.-J.; Chow, P.-T. Genetic algorithms for nonlinear mixed discrete-integer optimization problems via meta-genetic parameter optimization. Eng. Optim. 1995, 24, 137–159. [Google Scholar] [CrossRef]
  39. Guo, C.; Hu, J.; Ye, B.; Cao, Y. Swarm intelligence for mixed-variable design optimization. J. Zhejiang Univ. Sci. 2004, 5, 851–860. [Google Scholar] [CrossRef] [Green Version]
  40. Sandgren, E. Nonlinear integer and discrete programming in mechanical design optimization. ASME J. Mech. Des. 1990, 112, 223–229. [Google Scholar] [CrossRef]
  41. Siddall, J.N. Optimal Engineering Design: Principles and Applications; CRC Press: New York, NY, USA, 1982. [Google Scholar]
  42. Coello, C.A.C. Treating constraints as objectives for single-objective evolutionary optimization. Eng. Optim. 2000, 32, 275–308. [Google Scholar] [CrossRef]
  43. Deb, K. GeneAS: A robust optimal design technique for mechanical component design. In Evolutionary Algorithms in Engineering Applications; Springer: Berlin/Heidelberg, Germany, 1997; pp. 497–514. [Google Scholar]
  44. Liu, H.; Cai, Z.; Wang, Y. Hybridizing particle swarm optimization with differential evolution for constrained numerical and engineering optimization. Appl. Soft Comput. 2010, 10, 629–640. [Google Scholar] [CrossRef]
  45. Wang, L.; Li, L.P. An effective differential evolution with level comparison for constrained engineering design. Struct. Multidiscip. Optim. 2010, 41, 947–963. [Google Scholar] [CrossRef]
  46. Chen, T.Y.; Chen, H.C. Mixed-discrete structural optimization using a rank-niche evolution strategy. Eng. Optim. 2009, 41, 39–58. [Google Scholar] [CrossRef]
  47. Lamberti, L.; Pappalettere, C. Move limits definition in structural optimization with sequential linear programming. Part II: Numerical examples. Comput. Structres 2003, 81, 215–238. [Google Scholar] [CrossRef]
  48. Thanedar, P.B.; Vanderplaats, G.N. Survey of discrete variable optimization for structural design. J. Struct. Eng. 1995, 121, 301–306. [Google Scholar] [CrossRef]
  49. Erbatur, F.; Hasancebi, O.; Tutuncu, I.; Kilic, H. Optimal design of planar and space structures with genetic algorithms. Comput. Structres 2000, 75, 209–224. [Google Scholar] [CrossRef]
  50. Lemonge, A.C.C.; Barbosa, H.J.C. An adaptive penalty scheme for genetic algorithms in structural optimization. Int. J. Numer. Methods Eng. 2004, 59, 703–736. [Google Scholar] [CrossRef]
  51. Bernardino, H.S.; Barbosa, H.J.C.; Lemonge, A.C.C. A hybrid genetic algorithm for constrained optimization problems in mechanical engineering. In Proceedings of the 2007 IEEE Congress on Evolutionary Computation, Singapore, 25–28 September 2007; pp. 646–653. [Google Scholar]
  52. Gupta, S.; Tiwari, R.; Shivashankar, B.N. Multi-objective design optimization of rolling bearings using genetic algorithm. Mech. Mach. Theory 2007, 42, 1418–1443. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the proposed diversity-guided particle swarm optimizer (PSO) with a cyclic network mechanism.
Figure 1. Schematic diagram of the proposed diversity-guided particle swarm optimizer (PSO) with a cyclic network mechanism.
Mathematics 08 02016 g001
Figure 2. Flowchart of the PSO scheme with a diversity classifier for cyclically neighboring subpopulations for solving the constrained MIDC optimization problems.
Figure 2. Flowchart of the PSO scheme with a diversity classifier for cyclically neighboring subpopulations for solving the constrained MIDC optimization problems.
Mathematics 08 02016 g002
Figure 3. Configuration of a pressure vessel.
Figure 3. Configuration of a pressure vessel.
Mathematics 08 02016 g003
Figure 4. Configurationof a reinforced concrete beam.
Figure 4. Configurationof a reinforced concrete beam.
Mathematics 08 02016 g004
Figure 5. Configuration of a helical compression spring.
Figure 5. Configuration of a helical compression spring.
Mathematics 08 02016 g005
Figure 6. Configuration of a Belleville spring.
Figure 6. Configuration of a Belleville spring.
Mathematics 08 02016 g006
Figure 7. Configuration of a speed reducer.
Figure 7. Configuration of a speed reducer.
Mathematics 08 02016 g007
Figure 8. Configuration of a stepped cantilever beam.
Figure 8. Configuration of a stepped cantilever beam.
Mathematics 08 02016 g008
Figure 9. Configuration of a rolling element bearing.
Figure 9. Configuration of a rolling element bearing.
Mathematics 08 02016 g009
Figure 10. Finite element model utilized in the car side impact problem [1].
Figure 10. Finite element model utilized in the car side impact problem [1].
Mathematics 08 02016 g010
Figure 11. Distributions of f ( x * ) corresponding to n s / n p for eight numerical examples.
Figure 11. Distributions of f ( x * ) corresponding to n s / n p for eight numerical examples.
Mathematics 08 02016 g011
Figure 12. Convergence behaviors of f ( x * ) for eight numerical examples.
Figure 12. Convergence behaviors of f ( x * ) for eight numerical examples.
Mathematics 08 02016 g012
Table 1. Optimization results for the pressure vessel design problem.
Table 1. Optimization results for the pressure vessel design problem.
Method[32] PSO-GA[31] HS[33] SA-DS[1] FA[2] PSOPresent Study PSO ( n s = 16 )
x 1 (l)221.365487221.36553207.22555221.36547221.365548221.3654714
x 2 (r)38.86010238.8601039.8096238.8601038.86009938.8601036
x 3 ( t s )0.75000.750.76830.750.750000.7500000
x 4 ( t h )0.37500.3750.37970.3750.375000.3750000
h 1 ( x * ) −0.0000−0.0000−0.0000−0.0000−0.000000−0.0000000
h 2 ( x * ) −0.0043−0.0043−0.0000−0.0043−0.004275−0.0042746
h 3 ( x * ) 0.04460.2713−10.7065−0.0134−0.000190−0.0255540
h 4 ( x * ) −18.6345−18.6345−32.7744−18.6345−18.634452−18.6345290
Best objective value5850.3830645849.761695868.764845850.383065850.383765850.38306
Objective deviation a 1.0 + 6.837 × 10 10 1.0 1.062 × 10 4 1.0 + 3.142 × 10 3 1.0 1.0 + 1.197 × 10 7 1.0
FeasibilityInfeasibleInfeasibleFeasibleFeasibleFeasibleFeasible
Worst objective valueN/A bN/A6804.3281006258.968255850.5917975850.38308
Average objective valueN/AN/A6164. 5858675937.33790N/A5850.38306
Standard deviationN/AN/A257.473670164.54747N/A 3.7 × 10 6
Function evaluations100,000200,000N/A25,00031,436–124,96850,000
a The ratio of the best objective value of each method to the lowest one among all methods. b N/A denotes non-available.
Table 2. Discrete values of the reinforcing bars.
Table 2. Discrete values of the reinforcing bars.
Bar Type A s ( in 2 )Bar Type A s ( in 2 )Bar A s ( in 2 )Bar Type A s ( in 2 )
1#40.26#51.869#63.9612#77.2
1#50.3110#4, 2#924#9413#77.8
2#40.47#52.1713#54.0310#87.9
1#60.4411#4, 5#62.27#74.28#98
3#4, 1#70.63#82.3714#54.3414#78.4
2#50.6212#4, 4#72.410#64.411#88.69
1#80.798#52.4815#54.6515#79
4#40.813#42.66#84.7412#89.48
2#60.886#62.648#74.813#810.27
3#50.939#52.7911#64.8411#911
5#4, 1#9114#42.85#9514#811.06
6#4, 2#71.215#4, 5#7, 3#9312#65.2815#811.85
4#51.247#63.089#75.412#912
3#61.3210#53.107#85.5313#913
7#41.44#83.1613#85.7214#914
5#51.5511#53.4110 7, 6#9615#915
2#81.588#63.5214#66.16
8#41.66#73.68#86.32
4#61.7612#53.7215#6, 11#7, 7#96.6
9#4, 3#71.85#83.959#87.11
Table 3. Optimization results of the reinforced concrete beam design problem.
Table 3. Optimization results of the reinforced concrete beam design problem.
Reference Method[34] SD-RC a[35][36] BFO[37][1] FAPresent Study PSO ( n s = 8 )
GHN-ALM bGHN-EP cGAGA-FL
x 1 ( A s )7.86.66.32N/A7.206.166.326.32
x 2 (b)313334N/A32353434
x 3 (h)7.798.4952278.637180N/A8.04518.75008.50008.5000
h 1 ( x * ) −0.0205−0.1155−0.0635N/A−0.0224000
h 2 ( x * ) −4.20120.0159−0.7745N/A−2.8779−3.6173−0.2241−0.22409
Best objective value374.2362.2455362.00648376.2977366.1459364.8541359.2080359.2080
Objective deviation 1.0 + 4.174 × 10 2 1.0 + 8.456 × 10 3 1.0 + 7.791 × 10 3 1.0 + 4.758 × 10 2 1.0 + 1.931 × 10 2 1.0 + 1.572 × 10 2 1.01.0
FeasibilityFeasibleInfeasibleFeasibleN/AFeasibleFeasibleFeasibleFeasible
Worst objective valueN/AN/AN/AN/AN/AN/A669.150359.2080
Average objective valueN/AN/AN/AN/A371.5417365.8046460.706359.2080
Standard deviationN/AN/AN/AN/AN/AN/A80.738700.0000
Function evaluations396N/AN/A100,000100,00030,00025,00020,000
a Hybrid discrete steepest descent and rotating coordinate directions methods. b Generalized Hopfield network-based augmented Lagrange multiplier approach. c GHN-based extended penalty approach.
Table 4. Discrete allowable values of the spring wire diameter d.
Table 4. Discrete allowable values of the spring wire diameter d.
Allowable Wire Diameter (in)
0.00900.00950.01040.01180.01280.01320.0140
0.01500.01620.01730.01800.02000.02300.0250
0.02800.03200.03500.04100.04700.05400.0630
0.07200.08000.09200.10500.12000.13500.1480
0.16200.17700.19200.20700.22500.24400.2630
0.28300.30700.33100.36200.39400.43750.5000
Table 5. Parameter values used in the formulation of the helical compression spring problem.
Table 5. Parameter values used in the formulation of the helical compression spring problem.
NotationDescriptionValue
F max Maximum working load1000.0 (lb)
SMaximum allowable shear stress 1.89 × 10 5 (psi)
l max Maximum free length14.0 (in)
d min Minimum wire diameter0.2 (in)
D max Maximum outside spring diameter3.0 (in)
F p Preload compression force300.0 (lb)
δ p m Maximum allowable deflection6.0 (in)
under preload
δ w Deflection from preload position to1.25 (in)
maximum load position
GShear modulus of the material 1.15 × 10 7 (psi)
Table 6. Optimization results of the helical compression spring design problem.
Table 6. Optimization results of the helical compression spring design problem.
Reference Method[40] NLPA a[38] GA[39] HSIA[21] DE[14] PSO[1] FA[2] PSOPresent Study PSO ( n s = 16 )
x 1 (D)1.1807011.2274111.2231.2230411.2230411.2230491.2230411.223041
x 2 (N)109999999
x 3 (d)0.2830.2830.2830.2830.2830.2830.2830.283
h 1 ( x * ) 5430.9550.9931008.811008.81141008.81141008.021008.80591008.811398
h 2 ( x * ) 8.81878.92648.9468.945648.94568.9468.9456358.945636
h 3 ( x * ) 0.082980.08300.0830.083000.0830.0830.0830000.083000
h 4 ( x * ) 1.81931.77261.776961.776961.7771.7771.4939591.493959
h 5 ( x * ) 1.17231.33711.321701.321701.32171.3221.3217001.321700
h 6 ( x * ) 5.46435.45855.46435.464295.46435.4645.4642865.464286
h 7 ( x * ) 0.00.00.00.00.000000.0000000.000000
h 8 ( x * ) 0.00.01340.00.00.00000.00000.0000100.000000
Best objective value2.79952.66812.6592.658562.658562.6585762.6585592.658559
Objective deviation 1.0 + 5.301 × 10 2 1.0 + 3.589 × 10 3 1.0 + 1.659 × 10 4 1.0 + 3.761 × 10 7 1.0 + 3.761 × 10 7 1.0 + 6.394 × 10 6 1.01.0
Worst objective valueN/AN/AN/AN/AN/A7.8162919N/A2.660784
Average objective valueN/AN/AN/AN/A2.7380244.3835958N/A2.658890
Standard deviationN/AN/AN/AN/AN/A4.6076313N/A0.000611
Function evaluationsN/AN/AN/A26,00015,00075,0004784–98,99250,000
a Nonlinear programming algorithm.
Table 7. Variation of g ( a ) with a.
Table 7. Variation of g ( a ) with a.
a≤1.41.51.61.71.81.922.12.22.32.42.52.62.7≥2.8
g ( a ) 0.5810.850.770.710.660.630.60.560.550.530.520.510.510.5
Table 8. Optimization results of the Belleville spring design problem.
Table 8. Optimization results of the Belleville spring design problem.
Reference Method[42][43][41] OPTIVAR [3] NBV[30]Present Study PSO ( n s = 16 )
GeneAS-IGeneAS-IIABCTLBO
x 1 (t)0.2080.2050.2100.2040.204143N/A0.2041430.204143
x 2 (h)0.20.2010.2040.2000.2N/A0.20.200000
x 3 ( D i )8.7519.5349.26810.03010.0304732N/A10.0304710.030473
x 4 ( D e )11.06711.62711.49912.01012.01N/A12.0112.010000
h 1 ( x * ) 2145.4109−10.33962127.2624134.0816 4.58 × 10 4 N/A 1.77 × 10 6 6.93 × 10 11
h 2 ( x * ) 39.750182.8062194.222554−12.5328 3.04 × 10 7 N/A 7.46 × 10 8 8.85 × 10 10
h 3 ( x * ) 0.000000.00100.00400.0000 9.24 × 10 10 N/A 5.8 × 10 11 8.81 × 10 13
h 4 ( x * ) 1.5921.59401.58601.59601.595856N/A1.5958571.595857
h 5 ( x * ) 0.9430.38300.51100.00000N/A 2.35 × 10 9 1.58 × 10 10
h 6 ( x * ) 2.3162.09302.23101.98001.979526N/A1.9795271.979527
h 7 ( x * ) 0.213640.203970.208560.198990.198965N/A0.1989660.198966
Best objective value2.1219642.018072.162561.9787151.9796751.9796751.9796751.979675
Objective deviation 1.0 + 7.187 × 10 2 1.0 + 1.939 × 10 2 1.0 + 9.238 × 10 2 1.0 4.849 × 10 4 1.01.01.01.0
FeasibilityFeasibleInfeasibleFeasibleInfeasibleFeasibleN/AFeasibleFeasible
Worst objective valueN/AN/AN/AN/A2.0054312.1042971.9797571.979675
Average objective valueN/AN/AN/AN/A1.9846981.9954751.9796881.979675
Standard deviationN/AN/AN/AN/A7.78 × 10 3 0.070.450.000000
Function evaluationsN/AN/AN/AN/A15,000150,000150,00050,000
Table 9. Optimization results of the speed reducer design problem.
Table 9. Optimization results of the speed reducer design problem.
Reference MethodDEDSDELCHEAAMDEPSO-DEMBAPresent Study PSO ( n s = 16 )
x 1 (b)3.53.53.5000223.5000103.50000003.5000003.500000
x 2 (m)0.70.70.700000390.700000.7000000.7000000.700000
x 3 (n)171717.0000121717.00000017.00000017.000000
x 4 ( l 1 )7.37.37.3004277.3001567.3000007.3000337.300000
x 5 ( l 2 )7.7153197.7153197.7153777.8000277.8000007.7157727.800000
x 6 ( d 1 )3.3502143.3502143.3502303.3502213.3502143.3502182.900000
x 7 ( d 2 )5.2866545.2866545.2866635.2866855.28668325.2866545.286683
Best objective value2994.4710662994.4710662994.4991072996.3566892996.3481672994.4824532896.259285
Objective deviation 1.0 + 3.3910 × 10 2 1.0 + 3.3910 × 10 2 1.0 + 3.392 × 10 2 1.0 + 3.4561 × 10 2 1.0 + 3.4558 × 10 2 1.0 + 3.3914 × 10 2 1.0
Worst objective value2994.4710662994.4710662994.7523112996.3901372996.3482042999.6524442896.259380
Average objective value2994.4710662994.4710662994.6133682996.3672202996.3481742996.7690192896.259292
Standard deviation 3.6 × 10 12 1.9 × 10 12 7.0 × 10 2 8.2 × 10 3 6.4 × 10 6 1.560.000017
Function evaluations30,00030,00040,00024,00054,35025,00050,000
Table 10. Optimization results of the stepped cantilever beam design problem.
Table 10. Optimization results of the stepped cantilever beam design problem.
Method x 1 x 2 x 3 x 4 x 5 x 6 x 7 x 8 x 9 x 10 Objective FunctionFunction
( b 1 )( h 1 )( b 2 )( h 2 )( b 3 )( h 3 )( b 4 )( h 4 )( b 5 )( h 5 )BestDeviationWorstMeanStd. Dev.Evaluations
[46]RNES 1 a3603.1552.6502.31143.1081.82234.30764269.594 1.0 + 5.887 × 10 3 N/A12,000
RNES 23603.1552.6502.26743.7971.84934.28264322.433 1.0 + 6.714 × 10 3 N/A12,000
RNES 33603.1552.6502.34842.8041.78334.75364299.108 1.0 + 6.349 × 10 3 N/A12,000
RNES 43603.1552.6502.49141.512.11333.23165416.896 1.0 + 2.384 × 10 2 N/A12,000
[47]DOTN/A65391.59 1.0 + 2.345 × 10 2 N/AN/A
SLP bN/A65451.50 1.0 + 2.439 × 10 2 N/AN/A
MLD c-SLPN/A65352.20 1.0 + 2.283 × 10 2 N/AN/A
[48]C/RU d4623.1602.6552.20544.091.75135.0373555 1.0 + 1.512 × 10 1 N/AN/A
PD e3603.1552.6502.27645.5281.7534.99564537 1.0 + 1.007 × 10 2 N/AN/A
LAD f3603.1552.6502.26245.2331.7534.99564403 1.0 + 7.975 × 10 3 N/AN/A
CAD g3603.1552.6502.27945.5531.7535.00464558 1.0 + 1.040 × 10 2 N/AN/A
[49]GAOS Level 13603.1552.6502.345.51.83564815 1.0 + 1.442 × 10 2 N/A10,000
GAOS Level 23603.1552.6502.2745.251.753564447 1.0 + 8.664 × 10 3 N/A10,000
[50]GA-APM h3603.1552.6502.289445.62561.793134.59364698.56 1.0 + 1.260 × 10 2 73931.35968107.046N/A35,000
[51]AIS-GA3603.1552.6502.23544.3952.00432.87965559.60 1.0 + 2.608 × 10 2 77272.7870857.12N/A35,000
AIS-GA-C i3603.1602.6502.31143.1862.22531.25066533.47 1.0 + 4.132 × 10 2 76852.8671821.69N/A35,000
[1]FA3603.1552.6502.20544.0911.75034.99563893.52 1.0 + 1.409 × 10 6 64262.9942064144.75312175.9187950,000
Present studyPSO ( n s = 16 )3603.1552.6502.2045644.091111.7497634.9951463893.431.063893.4308063893.430800.0000050,000
a Rank-niche evolution strategy. b Sequential linear programming. c Move limit definition. d Continuous/round up. e Precise discrete. f Linear approximate discrete. g Conservative approximate discrete. h Adaptive penalty method. i Clearing.
Table 11. Optimization results of the rolling element bearing design problem.
Table 11. Optimization results of the rolling element bearing design problem.
Reference Method[52] GA[3] MBA[30]Present Study PSO (ns = 16)
ABCTLBO
x 1 ( D b )125.7171125.7153N/A125.7191125.719056
x 2 ( D m )21.42321.423300N/A21.4255921.425590
x 3 (Z)1111.000N/A1111.000000
x 4 ( f i )0.5150.515000N/A0.5150.515000
x 5 ( f o )0.5150.515000N/A0.5150.515000
x 6 ( K D min )0.41590.488805N/A0.4242660.411776
x 7 ( K D max )0.6510.627829N/A0.6339480.613510
x 8 ( ε )0.3000430.300149N/A0.30.300000
x 9 (e)0.02230.097305N/A0.0688580.059359
x 10 ( ζ )0.7510.646095N/A0.7994980.667473
h 1 ( x * ) 0.0008210N/A0 a0.000000
h 2 ( x * ) 13.7329998.630183N/A13.1525714.026828
h 3 ( x * ) 2.7240001.101429N/A1.52520.094509
h 4 ( x * ) −1.1072.040448N/A0.719056−1.401405
h 5 ( x * ) 0.7170000.715366N/A16.495440.719056
h 6 ( x * ) 4.85789923.611002N/A014.120649
h 7 ( x * ) 0.00210.000480N/A00.000000
h 8 ( x * ) 0.0000070N/A2.5593630.000000
h 9 ( x * ) 0.0000070N/A00.000000
Best objective value81,843.385,535.9611 b81,859.741681,859.7481,859.7416
Objective deviation 1.0 2.009 × 10 4 1.0 + 4.491 × 10 2 1.0 1.0 1.955 × 10 8 1.0
FeasibilityFeasibleInfeasibleN/AInfeasibleFeasible
Worst objective valueN/A84,440.194878,897.8180,807.855181,859.7401
Average objective valueN/A85,321.403081,49681,438.98781,859.7415
Standard deviationN/A211.520.690.660.0003
Function evaluations225,00050,00010,00010,00050,000
a Some constraint function values do not coincide with those calculated using the provided optimal solution. Their solution provides h 1 ( x * ) = 0.000004 , h 2 ( x * ) = 13.152560 , h 3 ( x * ) = 1.525180 , h 4 ( x * ) = 2.559350 , h 5 ( x * ) = 0.719100 , h 6 ( x * ) = 16.495400 , h 7 ( x * ) = 0.000022 , h 8 ( x * ) = 0.0 , and h 9 ( x * ) = 0.0 . b This best value does not coincide with that calculated using the provided optimal solution. Their solution provides f ( x * ) = 81843.68625 .
Table 12. Optimization results of the car side impact design problem.
Table 12. Optimization results of the car side impact design problem.
Reference Method[1]Present Study PSO ( n s = 16 )
PSODEGAFA
x 1 0.500000.500000.500050.500000.500000
x 2 1.116701.116701.280171.360001.116366
x 3 0.500000.50000.500010.500000.500000
x 4 1.302081.302081.033021.202001.302197
x 5 0.500000.500000.500010.500000.500000
x 6 1.500001.500000.500001.120001.500000
x 7 0.500000.500000.500000.500000.500000
x 8 0.345000.345000.34994 a0.345000.345000
x 9 0.192000.192000.192000.192000.192000
x 10 −19.54935−19.5493510.31198.87307−19.561544
x 11 −0.00431−0.004310.00167−18.99808−0.000190
Best objective value22.8447422.84298 b22.8565322.84298 c22.842969
Objective deviation 1.0 + 7.753 × 10 5 1.0 + 4.815 × 10 7 1.0 + 5.937 × 10 4 1.0 + 4.815 × 10 7 1.0
FeasibilityFeasibleFeasibleInfeasibleInfeasibleFeasible
Worst objective value23.2135424.1260626.24057824.0662322.846465
Average objective value22.8942923.2282823.5158522.8937622.843136
Standard deviation0.150170.344510.665550.166670.000649
a x 8 { 0.192 , 0.345 } . b The optimal solution is identical to that of PSO but different objective function values are provided in Gandomi et al. [1]. c This best value does not coincide with that calculated using the provided optimal solution, f ( x * ) = 24.06622 . Furthermore, their optimal solution may not guarantee the satisfaction of two constraint conditions, h 8 ( x * ) = 4.02129 > 4 kN and h 10 ( x ) = 15.84839 > 15.7 mm/ms.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kim, T.-H.; Cho, M.; Shin, S. Constrained Mixed-Variable Design Optimization Based on Particle Swarm Optimizer with a Diversity Classifier for Cyclically Neighboring Subpopulations. Mathematics 2020, 8, 2016. https://doi.org/10.3390/math8112016

AMA Style

Kim T-H, Cho M, Shin S. Constrained Mixed-Variable Design Optimization Based on Particle Swarm Optimizer with a Diversity Classifier for Cyclically Neighboring Subpopulations. Mathematics. 2020; 8(11):2016. https://doi.org/10.3390/math8112016

Chicago/Turabian Style

Kim, Tae-Hyoung, Minhaeng Cho, and Sangwoo Shin. 2020. "Constrained Mixed-Variable Design Optimization Based on Particle Swarm Optimizer with a Diversity Classifier for Cyclically Neighboring Subpopulations" Mathematics 8, no. 11: 2016. https://doi.org/10.3390/math8112016

APA Style

Kim, T. -H., Cho, M., & Shin, S. (2020). Constrained Mixed-Variable Design Optimization Based on Particle Swarm Optimizer with a Diversity Classifier for Cyclically Neighboring Subpopulations. Mathematics, 8(11), 2016. https://doi.org/10.3390/math8112016

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop