# A Hybrid SSA and SMA with Mutation Opposition-Based Learning for Constrained Engineering Problems

Authors: Shuang Wang; Qingxin Liu; Yuxiang Liu; Heming Jia; Laith Abualigah; Rong Zheng; Di Wu
Journal: Computational Intelligence and Neuroscience
 (2021)

---

## Abstract

Based on Salp Swarm Algorithm (SSA) and Slime Mould Algorithm (SMA), a novel hybrid optimization algorithm, named Hybrid Slime Mould Salp Swarm Algorithm (HSMSSA), is proposed to solve constrained engineering problems. SSA can obtain good results in solving some optimization problems. However, it is easy to suffer from local minima and lower density of population. SMA specializes in global exploration and good robustness, but its convergence rate is too slow to find satisfactory solutions efficiently. Thus, in this paper, considering the characteristics and advantages of both the above optimization algorithms, SMA is integrated into the leader position updating equations of SSA, which can share helpful information so that the proposed algorithm can utilize these two algorithms’ advantages to enhance global optimization performance. Furthermore, Levy flight is utilized to enhance the exploration ability. It is worth noting that a novel strategy called mutation opposition-based learning is proposed to enhance the performance of the hybrid optimization algorithm on premature convergence avoidance, balance between exploration and exploitation phases, and finding satisfactory global optimum. To evaluate the efficiency of the proposed algorithm, HSMSSA is applied to 23 different benchmark functions of the unimodal and multimodal types. Additionally, five classical constrained engineering problems are utilized to evaluate the proposed technique’s practicable abilities. The simulation results show that the HSMSSA method is more competitive and presents more engineering effectiveness for real-world constrained problems than SMA, SSA, and other comparative algorithms. In the end, we also provide some potential areas for future studies such as feature selection and multilevel threshold image segmentation.


---

## Body

## 1. Introduction

In recent years, metaheuristic algorithms have been widely concerned by a large number of scholars. Compared with other traditional optimization algorithms, the concept of metaheuristic algorithms is simple. Besides, they are flexible and can bypass local optima. Thus, metaheuristics have been successfully applied in different fields to solve various complex optimization problems in the real world [1–3].Metaheuristic algorithms include three main categories: evolution-based, physics-based, and swarm-based techniques. The inspirations of evolutionary-based methods are the laws of evolution in nature. The most popular evolution-based algorithms include Genetic Algorithm (GA) [4], Differential Evolution Algorithm (DE) [5], and Biogeography-Based Optimizer (BBO) [6]. Physics-based algorithms mimic the physical rules in the universe. There are some representative algorithms such as Simulated Annealing (SA) [7], Gravity Search Algorithm (GSA) [8], Black Hole (BH) algorithm [9], Multiverse Optimizer (MVO) [10], Artificial Chemical Reaction Optimization Algorithm (ACROA) [11], Ray Optimization (RO) [12], Curved Space Optimization (CSO) [13], Sine Cosine Algorithm (SCA) [14], Arithmetic Optimization Algorithm (AOA) [15], and Heat Transfer Relation-based Optimization Algorithm (HTOA) [16]. The third category algorithm is swarm-based techniques, which simulate the social behavior of creatures in nature. Some optimization techniques of this class include Particle Swarm Optimization (PSO) [17] Ant Colony Optimization Algorithm (ACO) [18], Firefly Algorithm (FA) [19], Grey Wolf Optimizer (GWO) [20], Cuckoo Search (CS) Algorithm [21], Whale Optimization Algorithm (WOA) [22], Bald Eagle Search (BES) algorithm [23], and Aquila Optimizer (AO) [24].It is worth noting that the most widely used swarm-based optimization algorithm is PSO [25]. PSO simulates the behavior of birds flying together in flocks. During the search, they all follow the best solutions in their paths. Cacciola et al. [26] discussed the problem of corrosion profile reconstruction starting from electrical data, in which PSO was utilized to obtain the image of the reconstructed corrosion profile. The result shows that PSO can obtain the optimal solution compared with LSM, and it takes the least time. This allows us to recognize the huge potential of the optimization algorithm.Salp Swarm Algorithm (SSA) [27] is a swarm-based algorithm proposed in 2017. SSA is inspired by swarm behavior, navigation, and foraging of salps in the ocean. Since SSA has fewer parameters and is easier to be realized in a program than other algorithms, SSA has been applied to many optimization problems, such as feature selection, image segmentation, and constrained engineering problems. However, like other metaheuristic algorithms, SSA may be easy to trap into local minima and lower population density. Therefore, many improved researches have been proposed to enhance the performance of SSA in many fields. Tubishat et al. [28] presented a Dynamic SSA (DSSA), which shows better accuracy than SSA in feature selection. Salgotra et al. [29] proposed a self-adaptive SSA to enhance exploitation ability and convergence speed. Neggaz et al. [30] proposed an improved leader in SSA using Sine Cosine Algorithm and disrupt operator for feature selection. Jia and Lang [31] presented an enhanced SSA with a crossover scheme and Levy flight to improve the movement patterns of salp leaders and followers. There are also other attempts on the hybrid algorithm of SSA. Saafan and El-gendy [32] proposed a hybrid improved Whale Optimization Salp Swarm Algorithm (IWOSSA). The IWOSSA achieves a better balance between exploration and exploitation phases and avoids premature convergence effectively. Singh et al. [33] developed a hybrid SSA with PSO, which integrated the advantages of SSA and PSO to eliminate trapping in local optima and unbalanced exploitation. Abadi et al. [34] proposed a hybrid approach by combining SSA with GA, which could obtain good results in solving some optimization problems.Slime Mould Algorithm (SMA) [35] is the latest swarm intelligence algorithm proposed in 2020. This algorithm simulates the oscillation mode and the foraging of Slime Mould in nature. SMA has a unique search mode, which keeps the algorithm from falling into local optima, and has superior global exploration capability. The approach has been applied in real-world optimization problems like feature selection [36], parameters optimization of the fuzzy system [37], multilevel threshold image segmentation [38], control scheme [39], and parallel connected multistacks fuel cells [40].Therefore, based on the capabilities of both above algorithms, we try to do a hybrid operation to improve the performance of SMA or SSA and then propose a new hybrid optimization algorithm (HSMSSA) to speed up the convergence rate and enhance the overall optimization performance. The specific method is that we integrate SMA as the leader role into SSA and retain the exploitation phase of SSA. At the same time, inspired by the significant performance of opposition-based learning and quasiopposition-based learning, we propose a new strategy named mutation opposition-based learning (MOBL), which switches the algorithm between opposition-based learning and quasiopposition-based learning through a mutation rate to increase the diversity of the population and speed up the convergence rate. In addition, Levy flight is utilized to improve SMA’s exploration capability and balance the exploration and exploitation phases of the algorithm. The proposed HSMSSA algorithm can improve both the exploration and exploitation abilities. The proposed HSMSSA is tested on 23 different benchmark functions and compared with other optimization algorithms. Furthermore, five constrained engineering problems are also utilized to evaluate HSMSSA’s capability on real-world optimization problems. The experimental results illustrate that the HSMSSA possesses the superior capability to search the global minimum and achieve less cost engineering design results than other state-of-the-art metaheuristic algorithms.The remainder of this paper is organized as follows. Section2 provides a brief overview of SSA, SMA, Levy flight, and mutation opposition-based learning strategy. Section 3 describes the proposed hybrid algorithm in detail. In Section 4, the details of benchmark functions, parameter settings of the selected algorithms, simulation experiments, and results analysis are introduced. Conclusions and prospects are given in Section 5.


## 2. Preliminaries

### 2.1. Salp Swarm Algorithm

In the deep sea, salps live in groups and form a salp chain to move and forage. In the salp chain, there are leaders and followers. The leader moves towards the food and guides the followers. In the process of moving, leaders explore globally, while followers thoroughly search locally [27]. The shapes of a salp and salp chain are shown in Figure 1.Figure 1
Individual salp and the swarm of salp [25].


#### 2.1.1. Leader Salps

The front salp of the chain is called the leader, so the following equation is used to perform this action to the salp leader:(1)Xj1=Fj+c1UB−LBr1+LB,r2<0.5,Fj−c1UB−LBr1+LB,r2≥0.5,(2)c1=2e−4t/T2,where Xj1 and Fj represent the new position of the leader and food source in the jth dimension and r1 and r2 are randomly generated numbers in the interval [0, 1]. It is worth noting that c1 is essential for SSA because it balances exploration and exploitation during the search process. t is the current iteration and T is the max iteration.


#### 2.1.2. Follower Salps

To update the position of the followers, the new concept is introduced, which is based on Newton’s law of motion as in the following equation:(3)Xji=12gt2+ω0t,i≥2,where Xji represents the position of ith follower salp in the jth dimension and g and ω0 indicate the acceleration and the velocity, respectively. The updating process of followers can be expressed as in the following equation:(4)Xji=12Xji+Xji−1.The pseudocode of SSA is presented in Algorithm1.Algorithm 1: Pseudocode of Salp Swarm Algorithm.
(1)
Initialize the population sizeN and max iteration T;(2)
Initialize the positions of salpXi (i = 1, 2, …, N)(3)
While (t ≤ T)(4)
Calculate fitness of each salp;(5)
Denote the best solution asF(6)
updatec1 by equation (2);(7)
Fori = 1 to Ndo(8)
if (i = = 1) then(9)
update position of leader salp by equation (1)(10)
Else(11)
update position of follower salp by equation (4)(12)
End if(13)
End for(14)
t = t + 1;(15)
End While(16)
Return the best solutionF;


### 2.2. Slime Mould Algorithm

The main idea of SMA is inspired by the behavior and morphological changes of Slime Mould in foraging. Slime Mould can dynamically change search mode based on the quality of food sources. If the food source has a high quality, the Slime Mould will use the region-limited search method. If the food concentration is low, the Slime Mould will explore other food sources in the search space. Furthermore, even if Slime Mould has found a high-quality food source, it still divides some individuals to explore another area in the region [35]. The behavior of Slime Mould can be mathematically described as follows:(5)Xt+1⟶=r4×UB−LB+LB,r3<z,Xbt⟶+vb⟶×W⟶⋅XAt⟶−XBt⟶,r5<p,vc⟶×Xt⟶,r5≥p,where parameters r3, r4, and r5 are random values in the range of 0 to 1. UB and LB indicate the upper and lower bound of search space. z is a constant. Xbt⟶ represents the best position obtained in all iterations, XAt⟶ and XBt⟶ represent two individuals selected randomly from the population, and Xt⟶ represents the location of Slime Mould. vc⟶ decreases linearly from one to zero, and vb⟶ is an oscillation parameter in the range [−a, a], in which a is calculated as follows:(6)a=arctanh−tT+1.The coefficientW⟶ is a very important parameter, which simulates the oscillation frequency of different food concentrations so that Slime Mould can approach food more quickly when they find high-quality food. The formula of W⟶ is listed as follows:(7)WSmellIndexi⟶=1+r6×logbF−SibF−wF+1,condition,1−r6×logbF−SibF−wF+1,others,(8)SmellIndex=sortS,where i ∈ 1, 2, …, N and S(i) represents the fitness of X. condition indicates that S(i) ranks first half of the Slime Mould, and r6 are random numbers uniformly generated in the interval of [0, 1]. bF represents the optimal fitness obtained in the current iterative process, wF represents the worst fitness value obtained in the iterative process currently, and SmellIndex denotes the sequence of fitness values sorted (ascends in the minimum value problem).Thep parameter can be described as follows:(9)p=tanhSi−DF,where DF represents the best fitness over all iterations. Figure 2 visualizes the general logic of SMA.Figure 2
The steps of SMA.The pseudocode of SMA is presented in Algorithm2.Algorithm 2: Pseudocode of Slime Mould Algorithm.
(1)
Initialize the population sizeN and max iteration T;(2)
Initialize the positions of Slime MouldXi (i = 1, 2, …, N)(3)
While (t ≤ T)(4)
Calculate fitness of each Slime Mould;(5)
updatebF, wF, and Xb;(6)
CalculateW by equation (3);(7)
Fori = 1 to Ndo(8)
updatep, vb, and vc;(9)
update positions by equation (1)(10)
End For(11)
t = t + 1;(12)
End While(13)
ReturnXb;


### 2.3. Levy Flight

Levy flight is an effective strategy for metaheuristic algorithms, successfully designed in many algorithms [41–44]. Levy flight is a class of non-Gaussian random processes that follows Levy distribution. It alternates between short-distance and occasionally long-distance walking, which can be inferred from Figure 3. The formula of Levy flight is as follows:(10)Levy=0.01×r7×σr81/β,σ=Γ1+β×sinπβ/2Γ1+β/2×β×2β−1/21/β,where r7 and r8 are random values in the range of [0, 1] and β is a constant equal to 1.5.Figure 3
Levy distribution and 2D Levy trajectory.
(a)(b)


### 2.4. Mutation Opposition-Based Learning

Opposition-based learning (OBL) was proposed by Tizhoosh in 2005 [45]. The essence of OBL is selecting the best solution to the next iteration by comparing the current solution and its opposition-based learning solution. The OBL strategy has been successfully used in varieties of metaheuristic algorithms [46–51] to improve the ability of local optima stagnation avoidance, and the mathematical expression is as follows:(11)XOBLt+1=LB+UB−Xt.Quasiopposition-based learning (QOBL) [52] is an improved version from OBL, which applies quasiopposite points instead of opposite points. These points produced through QOBL have more likelihood of being unknown solutions than the points created by OBL. The mathematical formula of QOBL is as follows:(12)XQOBLt+1=CS+r9×MP−CS,if MP>CS,MP+r9×CS−MP,otherwise,CS=LB+UB2,MP=LB+UB−Xt.Considering the superior performance of the two kinds of opposition-based learning, we propose mutation opposition-based learning (MOBL) by combining the mutation rate with these two opposition-based learning. By selecting the mutation rate, we can give full play to the characteristics of the OBL and QOBL and effectively enhance the ability of the algorithm to jump out of the local optima. Figure4 is an MOBL example, in which Figure 4(a) shows an objective function and Figure 4(b) displays three candidate solutions and their OBL solutions or QOBL solutions. The mathematical formula is as follows:(13)Xt+1=XOBLt,ifr10<rate,XQOBLt,else,where rate is mutation rate, and we set it to 0.1.Figure 4
MOBL example: (a) objective function and (b) candidate solutions (red point), its OBL solution (blue point), and its QOBL solution (green point).
(a)(b)


## 2.1. Salp Swarm Algorithm

In the deep sea, salps live in groups and form a salp chain to move and forage. In the salp chain, there are leaders and followers. The leader moves towards the food and guides the followers. In the process of moving, leaders explore globally, while followers thoroughly search locally [27]. The shapes of a salp and salp chain are shown in Figure 1.Figure 1
Individual salp and the swarm of salp [25].


### 2.1.1. Leader Salps

The front salp of the chain is called the leader, so the following equation is used to perform this action to the salp leader:(1)Xj1=Fj+c1UB−LBr1+LB,r2<0.5,Fj−c1UB−LBr1+LB,r2≥0.5,(2)c1=2e−4t/T2,where Xj1 and Fj represent the new position of the leader and food source in the jth dimension and r1 and r2 are randomly generated numbers in the interval [0, 1]. It is worth noting that c1 is essential for SSA because it balances exploration and exploitation during the search process. t is the current iteration and T is the max iteration.


### 2.1.2. Follower Salps

To update the position of the followers, the new concept is introduced, which is based on Newton’s law of motion as in the following equation:(3)Xji=12gt2+ω0t,i≥2,where Xji represents the position of ith follower salp in the jth dimension and g and ω0 indicate the acceleration and the velocity, respectively. The updating process of followers can be expressed as in the following equation:(4)Xji=12Xji+Xji−1.The pseudocode of SSA is presented in Algorithm1.Algorithm 1: Pseudocode of Salp Swarm Algorithm.
(1)
Initialize the population sizeN and max iteration T;(2)
Initialize the positions of salpXi (i = 1, 2, …, N)(3)
While (t ≤ T)(4)
Calculate fitness of each salp;(5)
Denote the best solution asF(6)
updatec1 by equation (2);(7)
Fori = 1 to Ndo(8)
if (i = = 1) then(9)
update position of leader salp by equation (1)(10)
Else(11)
update position of follower salp by equation (4)(12)
End if(13)
End for(14)
t = t + 1;(15)
End While(16)
Return the best solutionF;


## 2.1.1. Leader Salps

The front salp of the chain is called the leader, so the following equation is used to perform this action to the salp leader:(1)Xj1=Fj+c1UB−LBr1+LB,r2<0.5,Fj−c1UB−LBr1+LB,r2≥0.5,(2)c1=2e−4t/T2,where Xj1 and Fj represent the new position of the leader and food source in the jth dimension and r1 and r2 are randomly generated numbers in the interval [0, 1]. It is worth noting that c1 is essential for SSA because it balances exploration and exploitation during the search process. t is the current iteration and T is the max iteration.


## 2.1.2. Follower Salps

To update the position of the followers, the new concept is introduced, which is based on Newton’s law of motion as in the following equation:(3)Xji=12gt2+ω0t,i≥2,where Xji represents the position of ith follower salp in the jth dimension and g and ω0 indicate the acceleration and the velocity, respectively. The updating process of followers can be expressed as in the following equation:(4)Xji=12Xji+Xji−1.The pseudocode of SSA is presented in Algorithm1.Algorithm 1: Pseudocode of Salp Swarm Algorithm.
(1)
Initialize the population sizeN and max iteration T;(2)
Initialize the positions of salpXi (i = 1, 2, …, N)(3)
While (t ≤ T)(4)
Calculate fitness of each salp;(5)
Denote the best solution asF(6)
updatec1 by equation (2);(7)
Fori = 1 to Ndo(8)
if (i = = 1) then(9)
update position of leader salp by equation (1)(10)
Else(11)
update position of follower salp by equation (4)(12)
End if(13)
End for(14)
t = t + 1;(15)
End While(16)
Return the best solutionF;


## 2.2. Slime Mould Algorithm

The main idea of SMA is inspired by the behavior and morphological changes of Slime Mould in foraging. Slime Mould can dynamically change search mode based on the quality of food sources. If the food source has a high quality, the Slime Mould will use the region-limited search method. If the food concentration is low, the Slime Mould will explore other food sources in the search space. Furthermore, even if Slime Mould has found a high-quality food source, it still divides some individuals to explore another area in the region [35]. The behavior of Slime Mould can be mathematically described as follows:(5)Xt+1⟶=r4×UB−LB+LB,r3<z,Xbt⟶+vb⟶×W⟶⋅XAt⟶−XBt⟶,r5<p,vc⟶×Xt⟶,r5≥p,where parameters r3, r4, and r5 are random values in the range of 0 to 1. UB and LB indicate the upper and lower bound of search space. z is a constant. Xbt⟶ represents the best position obtained in all iterations, XAt⟶ and XBt⟶ represent two individuals selected randomly from the population, and Xt⟶ represents the location of Slime Mould. vc⟶ decreases linearly from one to zero, and vb⟶ is an oscillation parameter in the range [−a, a], in which a is calculated as follows:(6)a=arctanh−tT+1.The coefficientW⟶ is a very important parameter, which simulates the oscillation frequency of different food concentrations so that Slime Mould can approach food more quickly when they find high-quality food. The formula of W⟶ is listed as follows:(7)WSmellIndexi⟶=1+r6×logbF−SibF−wF+1,condition,1−r6×logbF−SibF−wF+1,others,(8)SmellIndex=sortS,where i ∈ 1, 2, …, N and S(i) represents the fitness of X. condition indicates that S(i) ranks first half of the Slime Mould, and r6 are random numbers uniformly generated in the interval of [0, 1]. bF represents the optimal fitness obtained in the current iterative process, wF represents the worst fitness value obtained in the iterative process currently, and SmellIndex denotes the sequence of fitness values sorted (ascends in the minimum value problem).Thep parameter can be described as follows:(9)p=tanhSi−DF,where DF represents the best fitness over all iterations. Figure 2 visualizes the general logic of SMA.Figure 2
The steps of SMA.The pseudocode of SMA is presented in Algorithm2.Algorithm 2: Pseudocode of Slime Mould Algorithm.
(1)
Initialize the population sizeN and max iteration T;(2)
Initialize the positions of Slime MouldXi (i = 1, 2, …, N)(3)
While (t ≤ T)(4)
Calculate fitness of each Slime Mould;(5)
updatebF, wF, and Xb;(6)
CalculateW by equation (3);(7)
Fori = 1 to Ndo(8)
updatep, vb, and vc;(9)
update positions by equation (1)(10)
End For(11)
t = t + 1;(12)
End While(13)
ReturnXb;


## 2.3. Levy Flight

Levy flight is an effective strategy for metaheuristic algorithms, successfully designed in many algorithms [41–44]. Levy flight is a class of non-Gaussian random processes that follows Levy distribution. It alternates between short-distance and occasionally long-distance walking, which can be inferred from Figure 3. The formula of Levy flight is as follows:(10)Levy=0.01×r7×σr81/β,σ=Γ1+β×sinπβ/2Γ1+β/2×β×2β−1/21/β,where r7 and r8 are random values in the range of [0, 1] and β is a constant equal to 1.5.Figure 3
Levy distribution and 2D Levy trajectory.
(a)(b)


## 2.4. Mutation Opposition-Based Learning

Opposition-based learning (OBL) was proposed by Tizhoosh in 2005 [45]. The essence of OBL is selecting the best solution to the next iteration by comparing the current solution and its opposition-based learning solution. The OBL strategy has been successfully used in varieties of metaheuristic algorithms [46–51] to improve the ability of local optima stagnation avoidance, and the mathematical expression is as follows:(11)XOBLt+1=LB+UB−Xt.Quasiopposition-based learning (QOBL) [52] is an improved version from OBL, which applies quasiopposite points instead of opposite points. These points produced through QOBL have more likelihood of being unknown solutions than the points created by OBL. The mathematical formula of QOBL is as follows:(12)XQOBLt+1=CS+r9×MP−CS,if MP>CS,MP+r9×CS−MP,otherwise,CS=LB+UB2,MP=LB+UB−Xt.Considering the superior performance of the two kinds of opposition-based learning, we propose mutation opposition-based learning (MOBL) by combining the mutation rate with these two opposition-based learning. By selecting the mutation rate, we can give full play to the characteristics of the OBL and QOBL and effectively enhance the ability of the algorithm to jump out of the local optima. Figure4 is an MOBL example, in which Figure 4(a) shows an objective function and Figure 4(b) displays three candidate solutions and their OBL solutions or QOBL solutions. The mathematical formula is as follows:(13)Xt+1=XOBLt,ifr10<rate,XQOBLt,else,where rate is mutation rate, and we set it to 0.1.Figure 4
MOBL example: (a) objective function and (b) candidate solutions (red point), its OBL solution (blue point), and its QOBL solution (green point).
(a)(b)


## 3. The Proposed Algorithm

### 3.1. Details of HSMSSA

In SSA, the population is divided into leader salps and follower salps: leader salps are the first half of salps in the chain, and follower salps follow the leader. However, the leader salp has poor randomness and is easy to fall into local optima. For the SMA algorithm, Slime Mould selects different search modes according to the positive and negative feedback of the current food concentration and has a certain probability of isolating some individuals to explore other regions in search space. These mechanisms increase the randomness of Slime Mould and enhance the ability to explore. Thevb parameter is utilized to realize the oscillation mode of Slime Mould, which is in the range of [−a, a]. However, vb has the drawback of low randomness, which cannot effectively simulate the process of Slime Mould looking for food sources. Therefore, we introduce Levy flight into the exploration phase to further enhance the exploration ability. Next, we integrate SMA into SSA, change the position update method of leader salps, and further improve the randomness of the algorithm through Levy flight. For followers, we propose a mutation opposition-based learning to enhance its population diversity and increase the ability of the algorithm to jump out of the local optima. The mathematical formula of leader salps is as follows:(14)Xt+1⟶=r4×UB−LB+LB,r3<z,Xb+vb×W×XA−XB×Levy,r5<p,vc×Xi,r5≥p.The pseudocode of HSMSSA is given in Algorithm3, and the summarized flowchart is displayed in Figure 5. As shown in Algorithm 3, the position of the population is initially generated randomly. Then, each individual’s fitness will be calculated. For the entire population in each iteration, parameter W is calculated using equation (7). The search agents of population size N are assigned to the two algorithms, which can utilize the advantages of SSA and SMA, and realize the sharing of helpful information to achieve global optimization. If the search agent belongs to the first half of the population, the position will be updated using equation (14) in SMA with Levy flight. Otherwise, the position is determined using equation (4) and MOBL. Finally, if the termination criteria are satisfied, the algorithm returns the best solution found so far; else the previous steps are repeated.Figure 5
Flowchart of HSMSSA.Algorithm 3: Pseudocode of HSMSSA.
(1)
Set the initial values of the population sizeN and the maximum number of iterations T(2)
Initialize positions of the populationX(3)
Whilet ≤ T(4)
Check if the position goes out of the search space boundary and bring it back.(5)
Calculate the fitness of each search agent(6)
UpdateW, bF, and wF(7)
Fori = 1 to N(8)
Ifi = = 1(9)
Ifr3 < z(10)
Update position by equation (14)(11)
Else(12)
IFr5 < p(13)
Update position by equation (14)(14)
Else(15)
Update position by equation (14)(16)
End if(17)
End if(18)
Else(19)
Update position by equation (4)(20)
Ifr10 < 0.1(21)
Update the position of MOBL using equation (13)(22)
Else(23)
Update the position of MOBL using equation (13)(24)
End if(25)
End if(26)
Check if the position goes out of the search space boundary and bring it back.(27)
select the best position into the next iteration.(28)
t = t + 1(29)
End for(30)
End while(31)
ReturnXbest


### 3.2. Computational Complexity Analysis

HSMSSA mainly consists of three components: initialization, fitness evaluation, and position updating. In the initialization phase, the computational complexity of positions generated is O(N × D), where D is the dimension size of the problem. Then, the computational complexity of fitness evaluation for the solution is O(N) during the iteration process. Finally, we utilize mutation opposition-based learning to keep the algorithm from falling into local optima; thus, the computational complexities of position updating of HSMSSA are O(2 × N × D). Therefore, the total computational complexity of the proposed HSMSSA algorithm is O(N × D + N + 2 × N × D).


## 3.1. Details of HSMSSA

In SSA, the population is divided into leader salps and follower salps: leader salps are the first half of salps in the chain, and follower salps follow the leader. However, the leader salp has poor randomness and is easy to fall into local optima. For the SMA algorithm, Slime Mould selects different search modes according to the positive and negative feedback of the current food concentration and has a certain probability of isolating some individuals to explore other regions in search space. These mechanisms increase the randomness of Slime Mould and enhance the ability to explore. Thevb parameter is utilized to realize the oscillation mode of Slime Mould, which is in the range of [−a, a]. However, vb has the drawback of low randomness, which cannot effectively simulate the process of Slime Mould looking for food sources. Therefore, we introduce Levy flight into the exploration phase to further enhance the exploration ability. Next, we integrate SMA into SSA, change the position update method of leader salps, and further improve the randomness of the algorithm through Levy flight. For followers, we propose a mutation opposition-based learning to enhance its population diversity and increase the ability of the algorithm to jump out of the local optima. The mathematical formula of leader salps is as follows:(14)Xt+1⟶=r4×UB−LB+LB,r3<z,Xb+vb×W×XA−XB×Levy,r5<p,vc×Xi,r5≥p.The pseudocode of HSMSSA is given in Algorithm3, and the summarized flowchart is displayed in Figure 5. As shown in Algorithm 3, the position of the population is initially generated randomly. Then, each individual’s fitness will be calculated. For the entire population in each iteration, parameter W is calculated using equation (7). The search agents of population size N are assigned to the two algorithms, which can utilize the advantages of SSA and SMA, and realize the sharing of helpful information to achieve global optimization. If the search agent belongs to the first half of the population, the position will be updated using equation (14) in SMA with Levy flight. Otherwise, the position is determined using equation (4) and MOBL. Finally, if the termination criteria are satisfied, the algorithm returns the best solution found so far; else the previous steps are repeated.Figure 5
Flowchart of HSMSSA.Algorithm 3: Pseudocode of HSMSSA.
(1)
Set the initial values of the population sizeN and the maximum number of iterations T(2)
Initialize positions of the populationX(3)
Whilet ≤ T(4)
Check if the position goes out of the search space boundary and bring it back.(5)
Calculate the fitness of each search agent(6)
UpdateW, bF, and wF(7)
Fori = 1 to N(8)
Ifi = = 1(9)
Ifr3 < z(10)
Update position by equation (14)(11)
Else(12)
IFr5 < p(13)
Update position by equation (14)(14)
Else(15)
Update position by equation (14)(16)
End if(17)
End if(18)
Else(19)
Update position by equation (4)(20)
Ifr10 < 0.1(21)
Update the position of MOBL using equation (13)(22)
Else(23)
Update the position of MOBL using equation (13)(24)
End if(25)
End if(26)
Check if the position goes out of the search space boundary and bring it back.(27)
select the best position into the next iteration.(28)
t = t + 1(29)
End for(30)
End while(31)
ReturnXbest


## 3.2. Computational Complexity Analysis

HSMSSA mainly consists of three components: initialization, fitness evaluation, and position updating. In the initialization phase, the computational complexity of positions generated is O(N × D), where D is the dimension size of the problem. Then, the computational complexity of fitness evaluation for the solution is O(N) during the iteration process. Finally, we utilize mutation opposition-based learning to keep the algorithm from falling into local optima; thus, the computational complexities of position updating of HSMSSA are O(2 × N × D). Therefore, the total computational complexity of the proposed HSMSSA algorithm is O(N × D + N + 2 × N × D).


## 4. Experimental Results and Discussion

This section compared the HSMSSA with some state-of-the-art metaheuristics algorithms on 23 benchmark functions to validate its performance. Moreover, five engineering design problems are employed as examples for real-world applications. The experimentations ran on Windows 10 with 24 GB RAM and Intel (R) i5-9500. All simulations were carried out using MATLAB R2020b.


### 4.1. Definition of 23 Benchmark Functions

To assess HSMSSA’s ability of exploration, exploitation, and escaping from local optima, 23 benchmark functions, including unimodal and multimodal functions, are tested [27]. The unimodal benchmark functions (F1–F7) are utilized to examine the exploitation ability of HSMSSA. The description of the unimodal benchmark function is shown in Table 1. The multimodal and fixed-dimension multimodal benchmark functions (F8–F23) shown in Tables 2 and 3 are used to test the exploration ability of HSMSSA.Table 1
Unimodal benchmark functions.
FunctionDimRangeFminF1x=∑i=1nxi230[−100, 100]0F2x=∑i=1nxi+∏i=1nxi30[−10, 10]0F3x=∑i=1n∑j−1ixj230[−100, 100]0F4x=maxixi,1≤i≤n30[−100, 100]0F5x=∑i=1n−1100xi+1−xi22+xi−1230[−30, 30]0F6x=∑i=1nxi+5230[−100, 100]0F7x=∑i=1nixi4+random0,130[−1.28, 1.28]0Table 2
Multimodal benchmark functions.
FunctionDimRangeFminF8x=∑i=1n−xisinxi30[−500, 500]−418.9829 × DimF9x=∑i=1nxi2−10cos2πxi+1030[−5.12, 5.12]F10x=−20exp−0.21/n∑i=1nxi2−exp1/n∑i=1ncos2πxi+20+e30[−32, 32]0F11x=1/4000∑i=1nxi2−∏i=1ncosxi/i+130[−600, 600]0F12x=π/n10sinπy1+∑i=1n−1yi−121+10sin2πyi+1+yn−12+∑i=1nuxi,10,100,4,whereyi=1+xi+1/4,uxi,a,k,m=kxi−am,xi>a,0,−a<xi<a,k−xi−am,xi<−a30[−50, 50]0F13x=0.1sin23πx1+∑i=1nxi−121+sin23πxi+1+xn−121+sin22πxn+∑i=1nuxi,5,100,430[−50, 50]0Table 3
Fixed-dimension multimodal benchmark functions.
FunctionDimRangeFminF14x=1/500+∑j=1251/j+∑i=12xi−aij6−12[−65, 65]0.998F15x=∑i=111ai−x1bi2+bix2/bi2+bix3+x424[−5, 5]0.00030F16x=4x12−2.1x14+1/3x16+x1x2−4x22+x242[−5, 5]−1.0316F17x=x2−5.1/4π2x12+5/πx1−62+101−1/8πcosx1+102[−5, 5]0.398F18x=1+x1+x2+1219−14x1+3x12−14x2+6x1x2+3x22×30+2x1−3x22×18−32x2+12x12+48x2−36x1x2+27x222[−2, 2]3F19x=−∑i=14ciexp−∑j=13aijxj−pij23[−1, 2]−3.86F20x=−∑i=14ciexp−∑j=16aijxj−pij26[0, 1]−3.32F21x=−∑i=15X−aiX−aiT+ci−14[0, 10]−10.1532F22x=−∑i=17X−aiX−aiT+ci−14[0, 10]−10.4028F23x=−∑i=110X−aiX−aiT+ci−14[0, 10]−10.5363In order to show the experimental results more representative, the HSMSSA is compared with the basic SMA [35] and SSA [27], AO [24], AOA [15], WOA [22], SCA [14], and MVO [10]. For all tests, we set the population size N = 30, dimension size D = 30, and maximum iteration T = 500, respectively, for all algorithms with 30 independent runs. The parameter settings of each algorithm are shown in Table 4. After all, average results and standard deviations are employed to evaluate the results. Note that the best results will be bolded.Table 4
Parameter settings for the comparative algorithms.
AlgorithmParametersSMA [35]z = 0.03SSA [27]c1 = [1, 0]; c2 ∈ [0, 1]; c3 ∈ [0, 1]AO [24]U = 0.00565; r1 = 10; ω = 0.005; α = 0.1; δ = 0.1; G1 ∈ [−1, 1]; G2 = [2, 0]AOA [15]α = 5; μ = 0.5;WOA [22]a1 = [2, 0]; a2 = [−1, −2]; b = 1SCA [14]a = [2, 0]MVO [10]WEP ∈ [0.2, 1]; TDR ∈ [0, 1]; r1, r2, r3 ∈ [0, 1]


#### 4.1.1. Evaluation of Exploitation Capability (F1–F7)

As we can see, unimodal benchmark functions have only one global optimum. These functions are allowed to evaluate the exploitation ability of the metaheuristic algorithms. It can be seen from Table5 that HSMSSA is very competitive with SMA, SSA, and other metaheuristic algorithms. In particular, HSMSSA can achieve much better results than other metaheuristic algorithms except F6. For F1–F4, HSMSSA can find the theoretical optimum. For all unimodal functions except F5, HSMSSA gets the smallest average values and standard deviations compared to other algorithms, which indicate the best accuracy and stability. Hence, the exploitation capability of the proposed HSMSSA algorithm is excellent.Table 5
Results of algorithms on 23 benchmark functions.
FunctionHSMSSASMASSAAOAOAWOASCAMVOF1Mean0.0000E + 003.3100E − 3032.1200E − 073.4900E − 1134.9600E − 063.0200E − 742.2039E + 011.2581E + 00Std0.0000E + 000.0000E + 002.3600E − 071.9100E − 1122.1600E − 061.6500E − 733.4334E + 014.4824E − 01F2Mean0.0000E + 008.7800E − 1611.9590E + 005.0900E − 591.9805E − 033.1700E − 502.6402E − 021.2927E + 00Std0.0000E + 004.8000E − 1601.4157E + 002.7900E − 581.9063E − 031.5800E − 492.9329E − 021.1310E + 00F3Mean0.0000E + 001.4159e − 3132.1110E + 033.7000E − 1019.6539E − 044.3306E + 048.0776E + 032.3785E + 02Std0.0000E + 000.0000E + 001.5514E + 032.0200E − 1008.1952E − 041.6945E + 044.7486E + 031.2902E + 02F4Mean0.0000E + 001.4800E − 1501.1761E + 011.7300E − 551.7772E − 025.2341E + 013.4937E + 011.9521E + 00Std0.0000E + 008.1300E − 1504.0647E + 009.5000E − 551.1446E − 022.8502E + 011.3547E + 017.5824E − 01F5Mean8.0753E − 025.0378E + 001.0827E + 025.2141E − 032.8034E + 012.8032E + 017.1843E + 044.7661E + 02Std2.7660E − 019.1628E + 001.4901E + 027.2117E − 032.0336E − 014.1099E − 011.2854E + 057.8524E + 02F6Mean8.2900E − 075.8093E − 031.3800E − 071.3859E − 043.1018E + 003.6069E − 011.1633E + 011.1004E + 00Std4.9600E − 072.4730E − 031.5300E − 072.1774E − 042.2067E − 012.2358E − 019.7843E + 003.6012E − 01F7Mean6.4400E − 051.6357E − 041.6394E − 017.8800E − 051.0695E − 043.0691E − 031.3682E − 013.6527E − 02Std6.2600E − 051.3341E − 046.1309E − 026.5100E − 051.0095E − 042.8218E − 032.0698E − 011.6353E − 02F8Mean−1.2569E + 04−1.2569E + 04−7.1073E + 03−9.3742E + 03−5.6142E + 03−1.0050E + 04−3.7684E + 03−7.6853E + 03Std2.0511E − 025.2108E − 019.0982E + 023.7928E + 033.9485E + 021.6846E + 032.8948E + 026.3161E + 02F9Mean0.0000E + 000.0000E + 005.3927E + 010.0000E + 001.5400E − 060.0000E + 003.6516E + 011.2767E + 02Std0.0000E + 000.0000E + 001.6886E + 010.0000E + 001.0800E − 060.0000E + 003.5214E + 013.5868E + 01F10Mean8.8800E − 168.8800E − 162.5258E + 002.1600E − 144.2571E − 044.3200E − 151.4566E + 012.5576E + 00Std0.0000E + 000.0000E + 006.8528E − 011.1400E − 131.8869E − 042.7200E − 158.3427E + 003.3443E + 00F11Mean0.0000E + 000.0000E + 002.0677E − 020.0000E + 005.1653E − 041.9620E − 029.3797E − 018.7497E − 01Std0.0000E + 000.0000E + 001.4808E − 020.0000E + 002.6954E − 035.4806E − 023.3652E − 017.4092E − 02F12Mean2.4700E − 055.9758E − 037.1903E + 003.3600E − 067.3782E − 011.7457E − 021.4353E + 032.1818E + 00Std5.1900E − 056.9167E − 032.7304E + 003.9400E − 062.7918E − 021.2428E − 027.0295E + 031.3736E + 00F13Mean4.6300E − 075.4042E − 031.2368E + 011.4900E − 052.9623E + 005.1556E − 011.8722E + 051.4761E − 01Std1.4800E − 077.9865E − 031.2274E + 012.0600E − 052.0518E − 022.3195E − 014.3886E + 056.2682E − 02F14Mean9.9800E − 019.9800E − 011.2294E + 002.4408E + 001.0093E + 012.9945E + 001.7920E + 009.9800E − 01Std5.0400E − 163.5900E − 136.7284E − 012.6435E + 003.9607E + 003.2280E + 009.8831E − 014.2300E − 11F15Mean3.9234E − 045.4854E − 042.8300E − 034.8202E − 047.5731E − 038.4669E − 041.0574E − 035.7644E − 03Std1.0289E − 042.5100E − 045.9501E − 031.0743E − 041.9781E − 026.5602E − 043.8453E − 041.2715E − 02F16Mean−1.0316E + 00−1.0316E + 00−1.0316E + 00−1.0311E + 00−1.0316E + 00−1.0316E + 00−1.0316E + 00−1.0316E + 00Std1.4300E − 136.6000E − 102.3300E − 124.6089E − 042.7600E − 111.0400E − 094.2000E − 053.8000E − 07F17Mean3.9789E − 013.9789E − 013.9789E − 013.9808E − 014.0051E − 013.9789E − 013.9921E − 013.9789E − 01Std2.6700E − 123.0000E − 081.3400E − 102.2131E − 041.4391E − 026.9400E − 061.2200E − 031.3300E − 07F18Mean3.9000E + 003.0000E + 003.0000E + 003.0365E + 001.2901E + 013.0003E + 003.0001E + 008.4000E + 00Std4.9295E + 001.8000E − 101.5000E − 133.6597E − 022.1835E + 011.5988E − 039.1200E − 052.0550E + 01F19Mean−3.8628E + 00−3.8628E + 00−3.8628E + 00−3.8560E + 00−3.7416E + 00−3.8575E + 00−3.8512E + 00−3.8628E + 00Std4.2800E − 118.4700E − 071.4600E − 107.4968E − 035.3652E − 015.7697E − 038.9408E − 031.9300E − 06F20Mean−3.2823E + 00−3.2426E + 00−3.2259E + 00−3.2075E + 00−3.2903E + 00−3.1996E + 00−2.9186E + 00−3.2655E + 00Std5.7048E − 025.7100E − 025.9965E − 027.2462E − 026.3492E − 021.2887E − 013.7689E − 016.1457E − 02F21Mean−1.0153E + 01−1.0153E + 01−6.3092E + 00−1.0142E + 01−7.4694E + 00−8.1869E + 00−2.4287E + 00−7.0425E + 00Std1.0300E − 071.9553E − 043.5233E + 002.6961E − 023.0153E + 002.6299E + 002.0334E + 003.0750E + 00F22Mean−1.0403E + 01−1.0403E + 01−7.8362E + 00−1.0392E + 01−7.3418E + 00−7.6682E + 00−3.4206E + 00−8.4065E + 00Std1.6600E − 073.1185E − 043.4923E + 001.8669E − 023.2433E + 003.4333E + 001.7953E + 002.9406E + 00F23Mean−1.0536E + 01−1.0536E + 01−9.4298E + 00−1.0525E + 01−7.9685E + 00−7.2480E + 00−3.6694E + 00−9.6618E + 00Std9.5800E − 082.5920E − 042.5607E + 001.6743E − 023.5174E + 003.2297E + 001.8159E + 002.3084E + 00


#### 4.1.2. Evaluation of Exploration Capability (F8–F23)

Unlike unimodal functions, multimodal functions have many local optima. Thus, this kind of test problem turns very useful to evaluate the exploration capability of an optimization algorithm. The results shown in Table5 for functions F8–F23 indicate that HSMSSA also has an excellent exploration capability. In fact, we can see that HSMSSA can find the theoretical optimum in F9, F11, F16–F17, and F19–F23. These results reveal that HSMSSA can also provide superior exploration capability.


#### 4.1.3. Analysis of Convergence Behavior

The convergence curves of some functions are selected and shown in Figure6, which show the convergence rate of algorithms. It can be seen that HSMSSA shows competitive performance compared to other state-of-the-art algorithms. The HSMSSA presents a faster convergence speed than all other algorithms in F7–F13, F15, and F19–F23. For other benchmark functions, HSMSSA shows a better capability of local optima avoidance than other comparison algorithms in F5 and F6.Figure 6
Convergence curves of 23 benchmark functions.


#### 4.1.4. Qualitative Results and Analysis

Furthermore, Figure7 shows the results of several representative test functions on search history, trajectory, average fitness, and convergence curve. From search history maps, we can see the search agent’s distribution of the HSMSSA while exploring and exploiting the search space. Because of the fast convergence, the vast majority of search agents are concentrated near the global optimum. Inspecting trajectory figure in Figure 5, the first search agent constantly oscillates in the first dimension of the search space, which suggests that the search agent investigates the most promising areas and better solutions widely. This powerful search capability is likely to come from the Levy flight and MOBL strategies. The average fitness presents if exploration and exploitation are conducive to improve the first random population, and an accurate approximation of the global optimum can be found in the end.Figure 7
Parameter space, search history, trajectory, average fitness, and convergence curves of HSMSSA.Similarly, it can be noticed that the average fitness oscillates in the early iterations and then decreases abruptly and begins to level off. The average fitness maps also show the significant improvement of the first random population and the final global optimal, accurate approximation acquisition. At last, the convergence curves reveal the best fitness value found by search agents after each iteration. By observing this, the HSMSSA shows breakneck convergence speed.


#### 4.1.5. Wilcoxon Signed-Rank Test

Because the algorithm results are random, we need to carry out statistical tests to prove that the results have statistical significance. We use Wilcoxon signed-ranks (WSR) test results to evaluate the statistical significance of the two algorithms at 5% significance level [53]. The WSR is a statistical test that is applied to two different results for searching the significantly different. As is well-known, a p-value less than 0.05 indicates that it is significantly superior to other algorithms. Otherwise, the obtained results are not statistically significant. The calculated results of the Wilcoxon signed-rank test between HSMSSA and other algorithms for each benchmark function are listed in Table 6. HSMSSA outperforms all other algorithms in varying degrees. This superiority is statistically significant on unimodal functions F2 and F4–F7, which indicates that HSMSSA possesses high exploitation. HSMSSA also shows better results on multimodal function F8–F23, suggesting that HSMSSA has a high capability of exploration. To sum up, HSMSSA can provide better results for almost all benchmark functions than other comparative algorithms.Table 6
The results of Wilcoxon’s sign rank test for all functions.
FunctionHSMSSA versusSMASSAAOAOAWOASCAMVOF1NaN6.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F26.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F3NaN6.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F46.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F56.1035E − 056.1035E − 054.2725E − 036.1035E − 056.1035E − 056.1035E − 056.1035E − 05F66.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F72.1545E − 026.1035E − 051.6882E − 017.1973E − 011.8311E − 046.1035E − 056.1035E − 05F86.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F9NaN6.1035E − 05NaN6.1035E − 05NaN6.1035E − 056.1035E − 05F10NaN6.1035E − 05NaN6.1035E − 051.9531E − 036.1035E − 056.1035E − 05F11NaN6.1035E − 05NaN6.1035E − 05NaN6.1035E − 056.1035E − 05F126.1035E − 056.1035E − 053.0518E − 046.1035E − 056.1035E − 056.1035E − 056.1035E − 05F136.1035E − 056.1035E − 053.0518E − 046.1035E − 056.1035E − 056.1035E − 056.1035E − 05F146.1035E − 055.0000E − 016.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F151.5259E − 036.1035E − 051.2207E − 044.3721E − 021.8311E − 046.1035E − 051.2207E − 04F166.1035E − 051.3184E − 026.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F176.1035E − 056.1035E − 056.1035E − 054.2725E − 046.1035E − 056.1035E − 056.1035E − 05F181.7334E − 016.4697E − 034.3252E − 023.5339E − 024.3252E − 024.3252E − 024.3252E − 02F196.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F202.5574E − 021.0254E − 021.1597E − 032.2931E − 014.5359E − 026.1035E − 053.5913E − 01F216.1035E − 052.5574E − 026.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F226.1035E − 054.2120E − 016.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F236.1035E − 054.2120E − 016.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05(W|L|T)(17|1|5)(20|3|0)(19|1|3)(21|2|0)(21|0|2)(23|0|0)(22|1|0)


### 4.2. Experiments on Engineering Design Problems

In this section, HSMSSA is evaluated to solve five classical engineering design problems: pressure vessel design problem, tension spring design problem, three-bar truss design problem, speed reducer problem, and cantilever beam design. To address these problems, we set the population sizeN = 30 and maximum iteration T = 500. The results of HSMSSA are compared to various state-of-the-art algorithms in the literature. The parameter settings are the same as previous numerical experiments.


#### 4.2.1. Pressure Vessel Design Problem

The pressure vessel design problem [53] is to minimize the total cost of cylindrical pressure vessel to match pressure requirements and form the pressure vessel shown in Figure 8. Four parameters in this problem need to be minimized, including the thickness of the shell (Ts), the thickness of head (Th), inner radius (R), and the length of the cylindrical section without the head (L), as shown in Figure 8. The constraints and equation are as follows.Figure 8
Pressure vessel design problem.Consider(15)x⟶=x1x2x3x4=TsThRL.Minimize(16)fx⟶=0.6224x1x3x4+1.7781x2x32+3.1661x12x4+19.84x12x3,subject to(17)g1x⟶=−x1+0.0193x3≤0,g2x⟶=−x3+0.00954x3≤0,g3x⟶=−πx32x4−43πx33+1296000≤0,g4x⟶=x4−240≤0.Variable range is(18)0≤x1≤99,0≤x2≤99,10≤x3≤200,10≤x4≤200.From the results in Table7, we can see that HSMSSA can obtain superior optimal values compared with SMA, SSA, AO, AOA, WOA, SCA, and MVO.Table 7
Comparison of HSMSSA results with other competitors for the pressure vessel design problem.
AlgorithmOptimum variablesOptimum costTsThRLHSMSSA0.85339920.418395645.8059135.53855961.5318SMA0.87448080.427383546.83887125.58486008.4161SSA0.89347740.438732747.53594119.21566047.6543AO0.88280920.431252447.26462121.64496028.1585AOA0.8071050.442651544.63354147.56596052.8917WOA0.83100910.364667144.00895154.31826047.0417SCA0.88200380.433508447.24144125.79226139.5293MVO0.83806770.499210145.82367135.36236253.5397


#### 4.2.2. Tension Spring Design Problem

This problem [27] tries to minimize the weight of the tension spring, and there are three parameters that need to be minimized, including the wire diameter (d), mean coil diameter (D), and the number of active coils (N). Figure 9 shows the structure of the tension spring. The mathematical of this problem can be written as follows.Figure 9
Tension spring design problem.Consider(19)x⟶=x1x2x3x4=dDN.Minimize(20)fx⟶=x3+2x2x12,subject to(21)g1x⟶=1−x23x371785x14≤0,g2x⟶=4x22−x1x212566x2x13−x14+15108x12≤0,g3x⟶=1−140.45x1x22x3≤0,g4x⟶=x1+x21.5−1≤0.Variable range is(22)0.05≤x1≤2.00,0.25≤x2≤1.30,2.00≤x3≤15.00.Results of HSMSSA for solving tension spring design problem are listed in Table8, which are compared with SMA, SSA, AO, AOA, WOA, SCA, and MVO. It is evident that HSMSSA obtained the best results compared to all other algorithms.Table 8
Comparison of HSMSSA results with other competitors for the tension spring design problem.
AlgorithmOptimum variablesOptimum weightdDNHSMSSA0.050.34809910.64860.011007SMA0.0573720.572324.14940.011584SSA0.059360.635423.47410.012256AO0.0569120.558334.32710.011442AOA0.0525790.410088.98430.012453WOA0.0552920.472776.9290.012906SCA0.05880.61723.650.012057MVO0.052510.3760210.335130.012790


#### 4.2.3. Three-Bar Truss Design Problem

Three-bar truss design is a complex problem in the field of civil engineering [49]. The goal of this problem is to achieve the minimum weight in truss design. Figure 10 shows the design of this problem. The formula of this problem can be described as follows.Figure 10
Three-bar truss design problem.Consider(23)x⟶=x1x2=A1A2.Minimize(24)fx⟶=22x1+x2∗l,subject to(25)g1x⟶=2x1+x22x12+2x1x2P−σ≤0,g2x⟶=x22x12+2x1x2P−σ≤0,g3x⟶=12x2+x1P−σ≤0.Variable range is(26)0≤x1,x2≤1,where l=100cm,P=2KN/cm2,andσ=2KN/cm2.Results of HSMSSA for solving the three-bar design problem are listed in Table9, which are compared with SMA, SSA, AO, AOA, WOA, SCA, and MVO. It can be observed that HSMSSA has an excellent ability to solve the problem in confined space.Table 9
Comparison of HSMSSA results with other competitors for the three-bar truss design problem.
AlgorithmOptimum variablesOptimum costx1x2HSMSSA0.788420.40811263.8523SMA0.791090.40022263.8668SSA0.778230.43736263.9363AO0.778990.43541263.9197AOA0.763420.48382264.3549WOA0.783570.42252263.886SCA0.780580.43766264.5463MVO0.754920.51216264.7851


#### 4.2.4. Speed Reducer Problem

In this problem [15], the total weight of the reducer is minimized by optimizing seven variables. Figure 11 shows the design of this problem, and the mathematical formula is as follows.Figure 11
Speed reducer problem.Minimize(27)fx⟶=0.7854x1x223.3333x32+14.9334x3−43.0934−1.508x1x62+x72+7.4777x63+x73,subject to(28)g1x⟶=27x1x22x3−1≤0,g2x⟶=397.5x1x22x32−1≤0,g3x⟶=1.93x43x2x3x64−1≤0,g4x⟶=1.93x53x2x3x74−1≤0,g5x⟶=745x4/x2x32+16.9×106110.0x63−1≤0,g6x⟶=745x4/x2x32+157.5×10685.0x63−1≤0,g7x⟶=x2x340−1≤0,g8x⟶=5x2x1−1≤0,g9x⟶=x112x2−1≤0,g10x⟶=1.5x6+1.9x4−1≤0,g11x⟶=1.1x7+1.9x5−1≤0.Variable range is(29)2.6≤x1≤3.6,0.7≤x2≤0.8,17≤x3≤28,7.3≤x4≤8.3,7.8≤x5≤8.3,2.9≤x6≤3.9,5.0≤x7≤5.5.The comparison results are listed in Table10, which shows the advantage of HSMSSA in realizing the minimum total weight of the problem.Table 10
Comparison of HSMSSA results with other competitors for the speed reducer design problem.
AlgorithmOptimum variablesOptimum weightx1x2x3x4x5x6x7HSMSSA3.49760.7177.37.83.350065.285532995.4374SMA3.514430.7177.324447.805273.352355.284943002.5941SSA3.497670.7177.877978.094013.39435.285743018.644AO3.51380.7177.414617.812913.377095.284573009.9097AOA3.559890.7177.499978.33.46235.285123061.8731WOA3.497390.7177.870398.07693.453095.28543034.0596SCA3.60.7177.38.33.379325.27993063.9102MVO3.60.7178.38.33.382765.360413111.6609


#### 4.2.5. Cantilever Beam Design

Cantilever beam design is a type of concrete engineering problem. This problem aims to determine the minimal total weight of the cantilever beam by optimizing the hollow square cross-section parameters [24]. Figure 12 illustrates the design of this problem, and the mathematical described is as follows.Figure 12
Cantilever beam design [24].Consider(30)x=x1x2x3x4x5.Minimize(31)fx⟶=0.6224x1+x2+x3+x4+x5,subject to(32)gx⟶=60x13+27x23+19x33+7x43+1x53−1≤0.Variable range is as follows:0.01≤x1,x2,x3,x4,x5≤100.The results are shown in Table11. From this table, we can see that the performance of HSMSSA is better than all other algorithms and the obtained total weight is minimized.Table 11
Comparison of HSMSSA results with other competitors for the cantilever beam design problem.
AlgorithmOptimum variablesOptimum weightx1x2x3x4x5HSMSSA5.99795.32894.49993.48172.1661.3400SMA6.08465.23094.56023.47252.13491.3405SSA6.10565.10014.3033.73652.31831.3456AO5.82365.38134.51783.47682.31591.3426AOA5.70145.37224.65773.66432.15511.3448WOA5.84925.05874.89173.6822.15641.3469SCA5.91655.91874.51333.31671.98231.3508MVO7.10885.01614.22033.40642.1051.3647As a summary, this section demonstrates the superiority of the proposed HSMSSA algorithm in different characteristics and real case studies. HSMSSA is able to outperform the basic SMA and SSA and other well-known algorithms with very competitive results, which are derived from the robust exploration and exploitation capabilities of HSMSSA. Excellent performance in solving industrial engineering design problems indicates that HSMSSA can be widely used in real-world optimization problems.


## 4.1. Definition of 23 Benchmark Functions

To assess HSMSSA’s ability of exploration, exploitation, and escaping from local optima, 23 benchmark functions, including unimodal and multimodal functions, are tested [27]. The unimodal benchmark functions (F1–F7) are utilized to examine the exploitation ability of HSMSSA. The description of the unimodal benchmark function is shown in Table 1. The multimodal and fixed-dimension multimodal benchmark functions (F8–F23) shown in Tables 2 and 3 are used to test the exploration ability of HSMSSA.Table 1
Unimodal benchmark functions.
FunctionDimRangeFminF1x=∑i=1nxi230[−100, 100]0F2x=∑i=1nxi+∏i=1nxi30[−10, 10]0F3x=∑i=1n∑j−1ixj230[−100, 100]0F4x=maxixi,1≤i≤n30[−100, 100]0F5x=∑i=1n−1100xi+1−xi22+xi−1230[−30, 30]0F6x=∑i=1nxi+5230[−100, 100]0F7x=∑i=1nixi4+random0,130[−1.28, 1.28]0Table 2
Multimodal benchmark functions.
FunctionDimRangeFminF8x=∑i=1n−xisinxi30[−500, 500]−418.9829 × DimF9x=∑i=1nxi2−10cos2πxi+1030[−5.12, 5.12]F10x=−20exp−0.21/n∑i=1nxi2−exp1/n∑i=1ncos2πxi+20+e30[−32, 32]0F11x=1/4000∑i=1nxi2−∏i=1ncosxi/i+130[−600, 600]0F12x=π/n10sinπy1+∑i=1n−1yi−121+10sin2πyi+1+yn−12+∑i=1nuxi,10,100,4,whereyi=1+xi+1/4,uxi,a,k,m=kxi−am,xi>a,0,−a<xi<a,k−xi−am,xi<−a30[−50, 50]0F13x=0.1sin23πx1+∑i=1nxi−121+sin23πxi+1+xn−121+sin22πxn+∑i=1nuxi,5,100,430[−50, 50]0Table 3
Fixed-dimension multimodal benchmark functions.
FunctionDimRangeFminF14x=1/500+∑j=1251/j+∑i=12xi−aij6−12[−65, 65]0.998F15x=∑i=111ai−x1bi2+bix2/bi2+bix3+x424[−5, 5]0.00030F16x=4x12−2.1x14+1/3x16+x1x2−4x22+x242[−5, 5]−1.0316F17x=x2−5.1/4π2x12+5/πx1−62+101−1/8πcosx1+102[−5, 5]0.398F18x=1+x1+x2+1219−14x1+3x12−14x2+6x1x2+3x22×30+2x1−3x22×18−32x2+12x12+48x2−36x1x2+27x222[−2, 2]3F19x=−∑i=14ciexp−∑j=13aijxj−pij23[−1, 2]−3.86F20x=−∑i=14ciexp−∑j=16aijxj−pij26[0, 1]−3.32F21x=−∑i=15X−aiX−aiT+ci−14[0, 10]−10.1532F22x=−∑i=17X−aiX−aiT+ci−14[0, 10]−10.4028F23x=−∑i=110X−aiX−aiT+ci−14[0, 10]−10.5363In order to show the experimental results more representative, the HSMSSA is compared with the basic SMA [35] and SSA [27], AO [24], AOA [15], WOA [22], SCA [14], and MVO [10]. For all tests, we set the population size N = 30, dimension size D = 30, and maximum iteration T = 500, respectively, for all algorithms with 30 independent runs. The parameter settings of each algorithm are shown in Table 4. After all, average results and standard deviations are employed to evaluate the results. Note that the best results will be bolded.Table 4
Parameter settings for the comparative algorithms.
AlgorithmParametersSMA [35]z = 0.03SSA [27]c1 = [1, 0]; c2 ∈ [0, 1]; c3 ∈ [0, 1]AO [24]U = 0.00565; r1 = 10; ω = 0.005; α = 0.1; δ = 0.1; G1 ∈ [−1, 1]; G2 = [2, 0]AOA [15]α = 5; μ = 0.5;WOA [22]a1 = [2, 0]; a2 = [−1, −2]; b = 1SCA [14]a = [2, 0]MVO [10]WEP ∈ [0.2, 1]; TDR ∈ [0, 1]; r1, r2, r3 ∈ [0, 1]


### 4.1.1. Evaluation of Exploitation Capability (F1–F7)

As we can see, unimodal benchmark functions have only one global optimum. These functions are allowed to evaluate the exploitation ability of the metaheuristic algorithms. It can be seen from Table5 that HSMSSA is very competitive with SMA, SSA, and other metaheuristic algorithms. In particular, HSMSSA can achieve much better results than other metaheuristic algorithms except F6. For F1–F4, HSMSSA can find the theoretical optimum. For all unimodal functions except F5, HSMSSA gets the smallest average values and standard deviations compared to other algorithms, which indicate the best accuracy and stability. Hence, the exploitation capability of the proposed HSMSSA algorithm is excellent.Table 5
Results of algorithms on 23 benchmark functions.
FunctionHSMSSASMASSAAOAOAWOASCAMVOF1Mean0.0000E + 003.3100E − 3032.1200E − 073.4900E − 1134.9600E − 063.0200E − 742.2039E + 011.2581E + 00Std0.0000E + 000.0000E + 002.3600E − 071.9100E − 1122.1600E − 061.6500E − 733.4334E + 014.4824E − 01F2Mean0.0000E + 008.7800E − 1611.9590E + 005.0900E − 591.9805E − 033.1700E − 502.6402E − 021.2927E + 00Std0.0000E + 004.8000E − 1601.4157E + 002.7900E − 581.9063E − 031.5800E − 492.9329E − 021.1310E + 00F3Mean0.0000E + 001.4159e − 3132.1110E + 033.7000E − 1019.6539E − 044.3306E + 048.0776E + 032.3785E + 02Std0.0000E + 000.0000E + 001.5514E + 032.0200E − 1008.1952E − 041.6945E + 044.7486E + 031.2902E + 02F4Mean0.0000E + 001.4800E − 1501.1761E + 011.7300E − 551.7772E − 025.2341E + 013.4937E + 011.9521E + 00Std0.0000E + 008.1300E − 1504.0647E + 009.5000E − 551.1446E − 022.8502E + 011.3547E + 017.5824E − 01F5Mean8.0753E − 025.0378E + 001.0827E + 025.2141E − 032.8034E + 012.8032E + 017.1843E + 044.7661E + 02Std2.7660E − 019.1628E + 001.4901E + 027.2117E − 032.0336E − 014.1099E − 011.2854E + 057.8524E + 02F6Mean8.2900E − 075.8093E − 031.3800E − 071.3859E − 043.1018E + 003.6069E − 011.1633E + 011.1004E + 00Std4.9600E − 072.4730E − 031.5300E − 072.1774E − 042.2067E − 012.2358E − 019.7843E + 003.6012E − 01F7Mean6.4400E − 051.6357E − 041.6394E − 017.8800E − 051.0695E − 043.0691E − 031.3682E − 013.6527E − 02Std6.2600E − 051.3341E − 046.1309E − 026.5100E − 051.0095E − 042.8218E − 032.0698E − 011.6353E − 02F8Mean−1.2569E + 04−1.2569E + 04−7.1073E + 03−9.3742E + 03−5.6142E + 03−1.0050E + 04−3.7684E + 03−7.6853E + 03Std2.0511E − 025.2108E − 019.0982E + 023.7928E + 033.9485E + 021.6846E + 032.8948E + 026.3161E + 02F9Mean0.0000E + 000.0000E + 005.3927E + 010.0000E + 001.5400E − 060.0000E + 003.6516E + 011.2767E + 02Std0.0000E + 000.0000E + 001.6886E + 010.0000E + 001.0800E − 060.0000E + 003.5214E + 013.5868E + 01F10Mean8.8800E − 168.8800E − 162.5258E + 002.1600E − 144.2571E − 044.3200E − 151.4566E + 012.5576E + 00Std0.0000E + 000.0000E + 006.8528E − 011.1400E − 131.8869E − 042.7200E − 158.3427E + 003.3443E + 00F11Mean0.0000E + 000.0000E + 002.0677E − 020.0000E + 005.1653E − 041.9620E − 029.3797E − 018.7497E − 01Std0.0000E + 000.0000E + 001.4808E − 020.0000E + 002.6954E − 035.4806E − 023.3652E − 017.4092E − 02F12Mean2.4700E − 055.9758E − 037.1903E + 003.3600E − 067.3782E − 011.7457E − 021.4353E + 032.1818E + 00Std5.1900E − 056.9167E − 032.7304E + 003.9400E − 062.7918E − 021.2428E − 027.0295E + 031.3736E + 00F13Mean4.6300E − 075.4042E − 031.2368E + 011.4900E − 052.9623E + 005.1556E − 011.8722E + 051.4761E − 01Std1.4800E − 077.9865E − 031.2274E + 012.0600E − 052.0518E − 022.3195E − 014.3886E + 056.2682E − 02F14Mean9.9800E − 019.9800E − 011.2294E + 002.4408E + 001.0093E + 012.9945E + 001.7920E + 009.9800E − 01Std5.0400E − 163.5900E − 136.7284E − 012.6435E + 003.9607E + 003.2280E + 009.8831E − 014.2300E − 11F15Mean3.9234E − 045.4854E − 042.8300E − 034.8202E − 047.5731E − 038.4669E − 041.0574E − 035.7644E − 03Std1.0289E − 042.5100E − 045.9501E − 031.0743E − 041.9781E − 026.5602E − 043.8453E − 041.2715E − 02F16Mean−1.0316E + 00−1.0316E + 00−1.0316E + 00−1.0311E + 00−1.0316E + 00−1.0316E + 00−1.0316E + 00−1.0316E + 00Std1.4300E − 136.6000E − 102.3300E − 124.6089E − 042.7600E − 111.0400E − 094.2000E − 053.8000E − 07F17Mean3.9789E − 013.9789E − 013.9789E − 013.9808E − 014.0051E − 013.9789E − 013.9921E − 013.9789E − 01Std2.6700E − 123.0000E − 081.3400E − 102.2131E − 041.4391E − 026.9400E − 061.2200E − 031.3300E − 07F18Mean3.9000E + 003.0000E + 003.0000E + 003.0365E + 001.2901E + 013.0003E + 003.0001E + 008.4000E + 00Std4.9295E + 001.8000E − 101.5000E − 133.6597E − 022.1835E + 011.5988E − 039.1200E − 052.0550E + 01F19Mean−3.8628E + 00−3.8628E + 00−3.8628E + 00−3.8560E + 00−3.7416E + 00−3.8575E + 00−3.8512E + 00−3.8628E + 00Std4.2800E − 118.4700E − 071.4600E − 107.4968E − 035.3652E − 015.7697E − 038.9408E − 031.9300E − 06F20Mean−3.2823E + 00−3.2426E + 00−3.2259E + 00−3.2075E + 00−3.2903E + 00−3.1996E + 00−2.9186E + 00−3.2655E + 00Std5.7048E − 025.7100E − 025.9965E − 027.2462E − 026.3492E − 021.2887E − 013.7689E − 016.1457E − 02F21Mean−1.0153E + 01−1.0153E + 01−6.3092E + 00−1.0142E + 01−7.4694E + 00−8.1869E + 00−2.4287E + 00−7.0425E + 00Std1.0300E − 071.9553E − 043.5233E + 002.6961E − 023.0153E + 002.6299E + 002.0334E + 003.0750E + 00F22Mean−1.0403E + 01−1.0403E + 01−7.8362E + 00−1.0392E + 01−7.3418E + 00−7.6682E + 00−3.4206E + 00−8.4065E + 00Std1.6600E − 073.1185E − 043.4923E + 001.8669E − 023.2433E + 003.4333E + 001.7953E + 002.9406E + 00F23Mean−1.0536E + 01−1.0536E + 01−9.4298E + 00−1.0525E + 01−7.9685E + 00−7.2480E + 00−3.6694E + 00−9.6618E + 00Std9.5800E − 082.5920E − 042.5607E + 001.6743E − 023.5174E + 003.2297E + 001.8159E + 002.3084E + 00


### 4.1.2. Evaluation of Exploration Capability (F8–F23)

Unlike unimodal functions, multimodal functions have many local optima. Thus, this kind of test problem turns very useful to evaluate the exploration capability of an optimization algorithm. The results shown in Table5 for functions F8–F23 indicate that HSMSSA also has an excellent exploration capability. In fact, we can see that HSMSSA can find the theoretical optimum in F9, F11, F16–F17, and F19–F23. These results reveal that HSMSSA can also provide superior exploration capability.


### 4.1.3. Analysis of Convergence Behavior

The convergence curves of some functions are selected and shown in Figure6, which show the convergence rate of algorithms. It can be seen that HSMSSA shows competitive performance compared to other state-of-the-art algorithms. The HSMSSA presents a faster convergence speed than all other algorithms in F7–F13, F15, and F19–F23. For other benchmark functions, HSMSSA shows a better capability of local optima avoidance than other comparison algorithms in F5 and F6.Figure 6
Convergence curves of 23 benchmark functions.


### 4.1.4. Qualitative Results and Analysis

Furthermore, Figure7 shows the results of several representative test functions on search history, trajectory, average fitness, and convergence curve. From search history maps, we can see the search agent’s distribution of the HSMSSA while exploring and exploiting the search space. Because of the fast convergence, the vast majority of search agents are concentrated near the global optimum. Inspecting trajectory figure in Figure 5, the first search agent constantly oscillates in the first dimension of the search space, which suggests that the search agent investigates the most promising areas and better solutions widely. This powerful search capability is likely to come from the Levy flight and MOBL strategies. The average fitness presents if exploration and exploitation are conducive to improve the first random population, and an accurate approximation of the global optimum can be found in the end.Figure 7
Parameter space, search history, trajectory, average fitness, and convergence curves of HSMSSA.Similarly, it can be noticed that the average fitness oscillates in the early iterations and then decreases abruptly and begins to level off. The average fitness maps also show the significant improvement of the first random population and the final global optimal, accurate approximation acquisition. At last, the convergence curves reveal the best fitness value found by search agents after each iteration. By observing this, the HSMSSA shows breakneck convergence speed.


### 4.1.5. Wilcoxon Signed-Rank Test

Because the algorithm results are random, we need to carry out statistical tests to prove that the results have statistical significance. We use Wilcoxon signed-ranks (WSR) test results to evaluate the statistical significance of the two algorithms at 5% significance level [53]. The WSR is a statistical test that is applied to two different results for searching the significantly different. As is well-known, a p-value less than 0.05 indicates that it is significantly superior to other algorithms. Otherwise, the obtained results are not statistically significant. The calculated results of the Wilcoxon signed-rank test between HSMSSA and other algorithms for each benchmark function are listed in Table 6. HSMSSA outperforms all other algorithms in varying degrees. This superiority is statistically significant on unimodal functions F2 and F4–F7, which indicates that HSMSSA possesses high exploitation. HSMSSA also shows better results on multimodal function F8–F23, suggesting that HSMSSA has a high capability of exploration. To sum up, HSMSSA can provide better results for almost all benchmark functions than other comparative algorithms.Table 6
The results of Wilcoxon’s sign rank test for all functions.
FunctionHSMSSA versusSMASSAAOAOAWOASCAMVOF1NaN6.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F26.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F3NaN6.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F46.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F56.1035E − 056.1035E − 054.2725E − 036.1035E − 056.1035E − 056.1035E − 056.1035E − 05F66.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F72.1545E − 026.1035E − 051.6882E − 017.1973E − 011.8311E − 046.1035E − 056.1035E − 05F86.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F9NaN6.1035E − 05NaN6.1035E − 05NaN6.1035E − 056.1035E − 05F10NaN6.1035E − 05NaN6.1035E − 051.9531E − 036.1035E − 056.1035E − 05F11NaN6.1035E − 05NaN6.1035E − 05NaN6.1035E − 056.1035E − 05F126.1035E − 056.1035E − 053.0518E − 046.1035E − 056.1035E − 056.1035E − 056.1035E − 05F136.1035E − 056.1035E − 053.0518E − 046.1035E − 056.1035E − 056.1035E − 056.1035E − 05F146.1035E − 055.0000E − 016.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F151.5259E − 036.1035E − 051.2207E − 044.3721E − 021.8311E − 046.1035E − 051.2207E − 04F166.1035E − 051.3184E − 026.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F176.1035E − 056.1035E − 056.1035E − 054.2725E − 046.1035E − 056.1035E − 056.1035E − 05F181.7334E − 016.4697E − 034.3252E − 023.5339E − 024.3252E − 024.3252E − 024.3252E − 02F196.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F202.5574E − 021.0254E − 021.1597E − 032.2931E − 014.5359E − 026.1035E − 053.5913E − 01F216.1035E − 052.5574E − 026.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F226.1035E − 054.2120E − 016.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F236.1035E − 054.2120E − 016.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05(W|L|T)(17|1|5)(20|3|0)(19|1|3)(21|2|0)(21|0|2)(23|0|0)(22|1|0)


## 4.1.1. Evaluation of Exploitation Capability (F1–F7)

As we can see, unimodal benchmark functions have only one global optimum. These functions are allowed to evaluate the exploitation ability of the metaheuristic algorithms. It can be seen from Table5 that HSMSSA is very competitive with SMA, SSA, and other metaheuristic algorithms. In particular, HSMSSA can achieve much better results than other metaheuristic algorithms except F6. For F1–F4, HSMSSA can find the theoretical optimum. For all unimodal functions except F5, HSMSSA gets the smallest average values and standard deviations compared to other algorithms, which indicate the best accuracy and stability. Hence, the exploitation capability of the proposed HSMSSA algorithm is excellent.Table 5
Results of algorithms on 23 benchmark functions.
FunctionHSMSSASMASSAAOAOAWOASCAMVOF1Mean0.0000E + 003.3100E − 3032.1200E − 073.4900E − 1134.9600E − 063.0200E − 742.2039E + 011.2581E + 00Std0.0000E + 000.0000E + 002.3600E − 071.9100E − 1122.1600E − 061.6500E − 733.4334E + 014.4824E − 01F2Mean0.0000E + 008.7800E − 1611.9590E + 005.0900E − 591.9805E − 033.1700E − 502.6402E − 021.2927E + 00Std0.0000E + 004.8000E − 1601.4157E + 002.7900E − 581.9063E − 031.5800E − 492.9329E − 021.1310E + 00F3Mean0.0000E + 001.4159e − 3132.1110E + 033.7000E − 1019.6539E − 044.3306E + 048.0776E + 032.3785E + 02Std0.0000E + 000.0000E + 001.5514E + 032.0200E − 1008.1952E − 041.6945E + 044.7486E + 031.2902E + 02F4Mean0.0000E + 001.4800E − 1501.1761E + 011.7300E − 551.7772E − 025.2341E + 013.4937E + 011.9521E + 00Std0.0000E + 008.1300E − 1504.0647E + 009.5000E − 551.1446E − 022.8502E + 011.3547E + 017.5824E − 01F5Mean8.0753E − 025.0378E + 001.0827E + 025.2141E − 032.8034E + 012.8032E + 017.1843E + 044.7661E + 02Std2.7660E − 019.1628E + 001.4901E + 027.2117E − 032.0336E − 014.1099E − 011.2854E + 057.8524E + 02F6Mean8.2900E − 075.8093E − 031.3800E − 071.3859E − 043.1018E + 003.6069E − 011.1633E + 011.1004E + 00Std4.9600E − 072.4730E − 031.5300E − 072.1774E − 042.2067E − 012.2358E − 019.7843E + 003.6012E − 01F7Mean6.4400E − 051.6357E − 041.6394E − 017.8800E − 051.0695E − 043.0691E − 031.3682E − 013.6527E − 02Std6.2600E − 051.3341E − 046.1309E − 026.5100E − 051.0095E − 042.8218E − 032.0698E − 011.6353E − 02F8Mean−1.2569E + 04−1.2569E + 04−7.1073E + 03−9.3742E + 03−5.6142E + 03−1.0050E + 04−3.7684E + 03−7.6853E + 03Std2.0511E − 025.2108E − 019.0982E + 023.7928E + 033.9485E + 021.6846E + 032.8948E + 026.3161E + 02F9Mean0.0000E + 000.0000E + 005.3927E + 010.0000E + 001.5400E − 060.0000E + 003.6516E + 011.2767E + 02Std0.0000E + 000.0000E + 001.6886E + 010.0000E + 001.0800E − 060.0000E + 003.5214E + 013.5868E + 01F10Mean8.8800E − 168.8800E − 162.5258E + 002.1600E − 144.2571E − 044.3200E − 151.4566E + 012.5576E + 00Std0.0000E + 000.0000E + 006.8528E − 011.1400E − 131.8869E − 042.7200E − 158.3427E + 003.3443E + 00F11Mean0.0000E + 000.0000E + 002.0677E − 020.0000E + 005.1653E − 041.9620E − 029.3797E − 018.7497E − 01Std0.0000E + 000.0000E + 001.4808E − 020.0000E + 002.6954E − 035.4806E − 023.3652E − 017.4092E − 02F12Mean2.4700E − 055.9758E − 037.1903E + 003.3600E − 067.3782E − 011.7457E − 021.4353E + 032.1818E + 00Std5.1900E − 056.9167E − 032.7304E + 003.9400E − 062.7918E − 021.2428E − 027.0295E + 031.3736E + 00F13Mean4.6300E − 075.4042E − 031.2368E + 011.4900E − 052.9623E + 005.1556E − 011.8722E + 051.4761E − 01Std1.4800E − 077.9865E − 031.2274E + 012.0600E − 052.0518E − 022.3195E − 014.3886E + 056.2682E − 02F14Mean9.9800E − 019.9800E − 011.2294E + 002.4408E + 001.0093E + 012.9945E + 001.7920E + 009.9800E − 01Std5.0400E − 163.5900E − 136.7284E − 012.6435E + 003.9607E + 003.2280E + 009.8831E − 014.2300E − 11F15Mean3.9234E − 045.4854E − 042.8300E − 034.8202E − 047.5731E − 038.4669E − 041.0574E − 035.7644E − 03Std1.0289E − 042.5100E − 045.9501E − 031.0743E − 041.9781E − 026.5602E − 043.8453E − 041.2715E − 02F16Mean−1.0316E + 00−1.0316E + 00−1.0316E + 00−1.0311E + 00−1.0316E + 00−1.0316E + 00−1.0316E + 00−1.0316E + 00Std1.4300E − 136.6000E − 102.3300E − 124.6089E − 042.7600E − 111.0400E − 094.2000E − 053.8000E − 07F17Mean3.9789E − 013.9789E − 013.9789E − 013.9808E − 014.0051E − 013.9789E − 013.9921E − 013.9789E − 01Std2.6700E − 123.0000E − 081.3400E − 102.2131E − 041.4391E − 026.9400E − 061.2200E − 031.3300E − 07F18Mean3.9000E + 003.0000E + 003.0000E + 003.0365E + 001.2901E + 013.0003E + 003.0001E + 008.4000E + 00Std4.9295E + 001.8000E − 101.5000E − 133.6597E − 022.1835E + 011.5988E − 039.1200E − 052.0550E + 01F19Mean−3.8628E + 00−3.8628E + 00−3.8628E + 00−3.8560E + 00−3.7416E + 00−3.8575E + 00−3.8512E + 00−3.8628E + 00Std4.2800E − 118.4700E − 071.4600E − 107.4968E − 035.3652E − 015.7697E − 038.9408E − 031.9300E − 06F20Mean−3.2823E + 00−3.2426E + 00−3.2259E + 00−3.2075E + 00−3.2903E + 00−3.1996E + 00−2.9186E + 00−3.2655E + 00Std5.7048E − 025.7100E − 025.9965E − 027.2462E − 026.3492E − 021.2887E − 013.7689E − 016.1457E − 02F21Mean−1.0153E + 01−1.0153E + 01−6.3092E + 00−1.0142E + 01−7.4694E + 00−8.1869E + 00−2.4287E + 00−7.0425E + 00Std1.0300E − 071.9553E − 043.5233E + 002.6961E − 023.0153E + 002.6299E + 002.0334E + 003.0750E + 00F22Mean−1.0403E + 01−1.0403E + 01−7.8362E + 00−1.0392E + 01−7.3418E + 00−7.6682E + 00−3.4206E + 00−8.4065E + 00Std1.6600E − 073.1185E − 043.4923E + 001.8669E − 023.2433E + 003.4333E + 001.7953E + 002.9406E + 00F23Mean−1.0536E + 01−1.0536E + 01−9.4298E + 00−1.0525E + 01−7.9685E + 00−7.2480E + 00−3.6694E + 00−9.6618E + 00Std9.5800E − 082.5920E − 042.5607E + 001.6743E − 023.5174E + 003.2297E + 001.8159E + 002.3084E + 00


## 4.1.2. Evaluation of Exploration Capability (F8–F23)

Unlike unimodal functions, multimodal functions have many local optima. Thus, this kind of test problem turns very useful to evaluate the exploration capability of an optimization algorithm. The results shown in Table5 for functions F8–F23 indicate that HSMSSA also has an excellent exploration capability. In fact, we can see that HSMSSA can find the theoretical optimum in F9, F11, F16–F17, and F19–F23. These results reveal that HSMSSA can also provide superior exploration capability.


## 4.1.3. Analysis of Convergence Behavior

The convergence curves of some functions are selected and shown in Figure6, which show the convergence rate of algorithms. It can be seen that HSMSSA shows competitive performance compared to other state-of-the-art algorithms. The HSMSSA presents a faster convergence speed than all other algorithms in F7–F13, F15, and F19–F23. For other benchmark functions, HSMSSA shows a better capability of local optima avoidance than other comparison algorithms in F5 and F6.Figure 6
Convergence curves of 23 benchmark functions.


## 4.1.4. Qualitative Results and Analysis

Furthermore, Figure7 shows the results of several representative test functions on search history, trajectory, average fitness, and convergence curve. From search history maps, we can see the search agent’s distribution of the HSMSSA while exploring and exploiting the search space. Because of the fast convergence, the vast majority of search agents are concentrated near the global optimum. Inspecting trajectory figure in Figure 5, the first search agent constantly oscillates in the first dimension of the search space, which suggests that the search agent investigates the most promising areas and better solutions widely. This powerful search capability is likely to come from the Levy flight and MOBL strategies. The average fitness presents if exploration and exploitation are conducive to improve the first random population, and an accurate approximation of the global optimum can be found in the end.Figure 7
Parameter space, search history, trajectory, average fitness, and convergence curves of HSMSSA.Similarly, it can be noticed that the average fitness oscillates in the early iterations and then decreases abruptly and begins to level off. The average fitness maps also show the significant improvement of the first random population and the final global optimal, accurate approximation acquisition. At last, the convergence curves reveal the best fitness value found by search agents after each iteration. By observing this, the HSMSSA shows breakneck convergence speed.


## 4.1.5. Wilcoxon Signed-Rank Test

Because the algorithm results are random, we need to carry out statistical tests to prove that the results have statistical significance. We use Wilcoxon signed-ranks (WSR) test results to evaluate the statistical significance of the two algorithms at 5% significance level [53]. The WSR is a statistical test that is applied to two different results for searching the significantly different. As is well-known, a p-value less than 0.05 indicates that it is significantly superior to other algorithms. Otherwise, the obtained results are not statistically significant. The calculated results of the Wilcoxon signed-rank test between HSMSSA and other algorithms for each benchmark function are listed in Table 6. HSMSSA outperforms all other algorithms in varying degrees. This superiority is statistically significant on unimodal functions F2 and F4–F7, which indicates that HSMSSA possesses high exploitation. HSMSSA also shows better results on multimodal function F8–F23, suggesting that HSMSSA has a high capability of exploration. To sum up, HSMSSA can provide better results for almost all benchmark functions than other comparative algorithms.Table 6
The results of Wilcoxon’s sign rank test for all functions.
FunctionHSMSSA versusSMASSAAOAOAWOASCAMVOF1NaN6.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F26.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F3NaN6.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F46.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F56.1035E − 056.1035E − 054.2725E − 036.1035E − 056.1035E − 056.1035E − 056.1035E − 05F66.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F72.1545E − 026.1035E − 051.6882E − 017.1973E − 011.8311E − 046.1035E − 056.1035E − 05F86.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F9NaN6.1035E − 05NaN6.1035E − 05NaN6.1035E − 056.1035E − 05F10NaN6.1035E − 05NaN6.1035E − 051.9531E − 036.1035E − 056.1035E − 05F11NaN6.1035E − 05NaN6.1035E − 05NaN6.1035E − 056.1035E − 05F126.1035E − 056.1035E − 053.0518E − 046.1035E − 056.1035E − 056.1035E − 056.1035E − 05F136.1035E − 056.1035E − 053.0518E − 046.1035E − 056.1035E − 056.1035E − 056.1035E − 05F146.1035E − 055.0000E − 016.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F151.5259E − 036.1035E − 051.2207E − 044.3721E − 021.8311E − 046.1035E − 051.2207E − 04F166.1035E − 051.3184E − 026.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F176.1035E − 056.1035E − 056.1035E − 054.2725E − 046.1035E − 056.1035E − 056.1035E − 05F181.7334E − 016.4697E − 034.3252E − 023.5339E − 024.3252E − 024.3252E − 024.3252E − 02F196.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F202.5574E − 021.0254E − 021.1597E − 032.2931E − 014.5359E − 026.1035E − 053.5913E − 01F216.1035E − 052.5574E − 026.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F226.1035E − 054.2120E − 016.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05F236.1035E − 054.2120E − 016.1035E − 056.1035E − 056.1035E − 056.1035E − 056.1035E − 05(W|L|T)(17|1|5)(20|3|0)(19|1|3)(21|2|0)(21|0|2)(23|0|0)(22|1|0)


## 4.2. Experiments on Engineering Design Problems

In this section, HSMSSA is evaluated to solve five classical engineering design problems: pressure vessel design problem, tension spring design problem, three-bar truss design problem, speed reducer problem, and cantilever beam design. To address these problems, we set the population sizeN = 30 and maximum iteration T = 500. The results of HSMSSA are compared to various state-of-the-art algorithms in the literature. The parameter settings are the same as previous numerical experiments.


### 4.2.1. Pressure Vessel Design Problem

The pressure vessel design problem [53] is to minimize the total cost of cylindrical pressure vessel to match pressure requirements and form the pressure vessel shown in Figure 8. Four parameters in this problem need to be minimized, including the thickness of the shell (Ts), the thickness of head (Th), inner radius (R), and the length of the cylindrical section without the head (L), as shown in Figure 8. The constraints and equation are as follows.Figure 8
Pressure vessel design problem.Consider(15)x⟶=x1x2x3x4=TsThRL.Minimize(16)fx⟶=0.6224x1x3x4+1.7781x2x32+3.1661x12x4+19.84x12x3,subject to(17)g1x⟶=−x1+0.0193x3≤0,g2x⟶=−x3+0.00954x3≤0,g3x⟶=−πx32x4−43πx33+1296000≤0,g4x⟶=x4−240≤0.Variable range is(18)0≤x1≤99,0≤x2≤99,10≤x3≤200,10≤x4≤200.From the results in Table7, we can see that HSMSSA can obtain superior optimal values compared with SMA, SSA, AO, AOA, WOA, SCA, and MVO.Table 7
Comparison of HSMSSA results with other competitors for the pressure vessel design problem.
AlgorithmOptimum variablesOptimum costTsThRLHSMSSA0.85339920.418395645.8059135.53855961.5318SMA0.87448080.427383546.83887125.58486008.4161SSA0.89347740.438732747.53594119.21566047.6543AO0.88280920.431252447.26462121.64496028.1585AOA0.8071050.442651544.63354147.56596052.8917WOA0.83100910.364667144.00895154.31826047.0417SCA0.88200380.433508447.24144125.79226139.5293MVO0.83806770.499210145.82367135.36236253.5397


### 4.2.2. Tension Spring Design Problem

This problem [27] tries to minimize the weight of the tension spring, and there are three parameters that need to be minimized, including the wire diameter (d), mean coil diameter (D), and the number of active coils (N). Figure 9 shows the structure of the tension spring. The mathematical of this problem can be written as follows.Figure 9
Tension spring design problem.Consider(19)x⟶=x1x2x3x4=dDN.Minimize(20)fx⟶=x3+2x2x12,subject to(21)g1x⟶=1−x23x371785x14≤0,g2x⟶=4x22−x1x212566x2x13−x14+15108x12≤0,g3x⟶=1−140.45x1x22x3≤0,g4x⟶=x1+x21.5−1≤0.Variable range is(22)0.05≤x1≤2.00,0.25≤x2≤1.30,2.00≤x3≤15.00.Results of HSMSSA for solving tension spring design problem are listed in Table8, which are compared with SMA, SSA, AO, AOA, WOA, SCA, and MVO. It is evident that HSMSSA obtained the best results compared to all other algorithms.Table 8
Comparison of HSMSSA results with other competitors for the tension spring design problem.
AlgorithmOptimum variablesOptimum weightdDNHSMSSA0.050.34809910.64860.011007SMA0.0573720.572324.14940.011584SSA0.059360.635423.47410.012256AO0.0569120.558334.32710.011442AOA0.0525790.410088.98430.012453WOA0.0552920.472776.9290.012906SCA0.05880.61723.650.012057MVO0.052510.3760210.335130.012790


### 4.2.3. Three-Bar Truss Design Problem

Three-bar truss design is a complex problem in the field of civil engineering [49]. The goal of this problem is to achieve the minimum weight in truss design. Figure 10 shows the design of this problem. The formula of this problem can be described as follows.Figure 10
Three-bar truss design problem.Consider(23)x⟶=x1x2=A1A2.Minimize(24)fx⟶=22x1+x2∗l,subject to(25)g1x⟶=2x1+x22x12+2x1x2P−σ≤0,g2x⟶=x22x12+2x1x2P−σ≤0,g3x⟶=12x2+x1P−σ≤0.Variable range is(26)0≤x1,x2≤1,where l=100cm,P=2KN/cm2,andσ=2KN/cm2.Results of HSMSSA for solving the three-bar design problem are listed in Table9, which are compared with SMA, SSA, AO, AOA, WOA, SCA, and MVO. It can be observed that HSMSSA has an excellent ability to solve the problem in confined space.Table 9
Comparison of HSMSSA results with other competitors for the three-bar truss design problem.
AlgorithmOptimum variablesOptimum costx1x2HSMSSA0.788420.40811263.8523SMA0.791090.40022263.8668SSA0.778230.43736263.9363AO0.778990.43541263.9197AOA0.763420.48382264.3549WOA0.783570.42252263.886SCA0.780580.43766264.5463MVO0.754920.51216264.7851


### 4.2.4. Speed Reducer Problem

In this problem [15], the total weight of the reducer is minimized by optimizing seven variables. Figure 11 shows the design of this problem, and the mathematical formula is as follows.Figure 11
Speed reducer problem.Minimize(27)fx⟶=0.7854x1x223.3333x32+14.9334x3−43.0934−1.508x1x62+x72+7.4777x63+x73,subject to(28)g1x⟶=27x1x22x3−1≤0,g2x⟶=397.5x1x22x32−1≤0,g3x⟶=1.93x43x2x3x64−1≤0,g4x⟶=1.93x53x2x3x74−1≤0,g5x⟶=745x4/x2x32+16.9×106110.0x63−1≤0,g6x⟶=745x4/x2x32+157.5×10685.0x63−1≤0,g7x⟶=x2x340−1≤0,g8x⟶=5x2x1−1≤0,g9x⟶=x112x2−1≤0,g10x⟶=1.5x6+1.9x4−1≤0,g11x⟶=1.1x7+1.9x5−1≤0.Variable range is(29)2.6≤x1≤3.6,0.7≤x2≤0.8,17≤x3≤28,7.3≤x4≤8.3,7.8≤x5≤8.3,2.9≤x6≤3.9,5.0≤x7≤5.5.The comparison results are listed in Table10, which shows the advantage of HSMSSA in realizing the minimum total weight of the problem.Table 10
Comparison of HSMSSA results with other competitors for the speed reducer design problem.
AlgorithmOptimum variablesOptimum weightx1x2x3x4x5x6x7HSMSSA3.49760.7177.37.83.350065.285532995.4374SMA3.514430.7177.324447.805273.352355.284943002.5941SSA3.497670.7177.877978.094013.39435.285743018.644AO3.51380.7177.414617.812913.377095.284573009.9097AOA3.559890.7177.499978.33.46235.285123061.8731WOA3.497390.7177.870398.07693.453095.28543034.0596SCA3.60.7177.38.33.379325.27993063.9102MVO3.60.7178.38.33.382765.360413111.6609


### 4.2.5. Cantilever Beam Design

Cantilever beam design is a type of concrete engineering problem. This problem aims to determine the minimal total weight of the cantilever beam by optimizing the hollow square cross-section parameters [24]. Figure 12 illustrates the design of this problem, and the mathematical described is as follows.Figure 12
Cantilever beam design [24].Consider(30)x=x1x2x3x4x5.Minimize(31)fx⟶=0.6224x1+x2+x3+x4+x5,subject to(32)gx⟶=60x13+27x23+19x33+7x43+1x53−1≤0.Variable range is as follows:0.01≤x1,x2,x3,x4,x5≤100.The results are shown in Table11. From this table, we can see that the performance of HSMSSA is better than all other algorithms and the obtained total weight is minimized.Table 11
Comparison of HSMSSA results with other competitors for the cantilever beam design problem.
AlgorithmOptimum variablesOptimum weightx1x2x3x4x5HSMSSA5.99795.32894.49993.48172.1661.3400SMA6.08465.23094.56023.47252.13491.3405SSA6.10565.10014.3033.73652.31831.3456AO5.82365.38134.51783.47682.31591.3426AOA5.70145.37224.65773.66432.15511.3448WOA5.84925.05874.89173.6822.15641.3469SCA5.91655.91874.51333.31671.98231.3508MVO7.10885.01614.22033.40642.1051.3647As a summary, this section demonstrates the superiority of the proposed HSMSSA algorithm in different characteristics and real case studies. HSMSSA is able to outperform the basic SMA and SSA and other well-known algorithms with very competitive results, which are derived from the robust exploration and exploitation capabilities of HSMSSA. Excellent performance in solving industrial engineering design problems indicates that HSMSSA can be widely used in real-world optimization problems.


## 4.2.1. Pressure Vessel Design Problem

The pressure vessel design problem [53] is to minimize the total cost of cylindrical pressure vessel to match pressure requirements and form the pressure vessel shown in Figure 8. Four parameters in this problem need to be minimized, including the thickness of the shell (Ts), the thickness of head (Th), inner radius (R), and the length of the cylindrical section without the head (L), as shown in Figure 8. The constraints and equation are as follows.Figure 8
Pressure vessel design problem.Consider(15)x⟶=x1x2x3x4=TsThRL.Minimize(16)fx⟶=0.6224x1x3x4+1.7781x2x32+3.1661x12x4+19.84x12x3,subject to(17)g1x⟶=−x1+0.0193x3≤0,g2x⟶=−x3+0.00954x3≤0,g3x⟶=−πx32x4−43πx33+1296000≤0,g4x⟶=x4−240≤0.Variable range is(18)0≤x1≤99,0≤x2≤99,10≤x3≤200,10≤x4≤200.From the results in Table7, we can see that HSMSSA can obtain superior optimal values compared with SMA, SSA, AO, AOA, WOA, SCA, and MVO.Table 7
Comparison of HSMSSA results with other competitors for the pressure vessel design problem.
AlgorithmOptimum variablesOptimum costTsThRLHSMSSA0.85339920.418395645.8059135.53855961.5318SMA0.87448080.427383546.83887125.58486008.4161SSA0.89347740.438732747.53594119.21566047.6543AO0.88280920.431252447.26462121.64496028.1585AOA0.8071050.442651544.63354147.56596052.8917WOA0.83100910.364667144.00895154.31826047.0417SCA0.88200380.433508447.24144125.79226139.5293MVO0.83806770.499210145.82367135.36236253.5397


## 4.2.2. Tension Spring Design Problem

This problem [27] tries to minimize the weight of the tension spring, and there are three parameters that need to be minimized, including the wire diameter (d), mean coil diameter (D), and the number of active coils (N). Figure 9 shows the structure of the tension spring. The mathematical of this problem can be written as follows.Figure 9
Tension spring design problem.Consider(19)x⟶=x1x2x3x4=dDN.Minimize(20)fx⟶=x3+2x2x12,subject to(21)g1x⟶=1−x23x371785x14≤0,g2x⟶=4x22−x1x212566x2x13−x14+15108x12≤0,g3x⟶=1−140.45x1x22x3≤0,g4x⟶=x1+x21.5−1≤0.Variable range is(22)0.05≤x1≤2.00,0.25≤x2≤1.30,2.00≤x3≤15.00.Results of HSMSSA for solving tension spring design problem are listed in Table8, which are compared with SMA, SSA, AO, AOA, WOA, SCA, and MVO. It is evident that HSMSSA obtained the best results compared to all other algorithms.Table 8
Comparison of HSMSSA results with other competitors for the tension spring design problem.
AlgorithmOptimum variablesOptimum weightdDNHSMSSA0.050.34809910.64860.011007SMA0.0573720.572324.14940.011584SSA0.059360.635423.47410.012256AO0.0569120.558334.32710.011442AOA0.0525790.410088.98430.012453WOA0.0552920.472776.9290.012906SCA0.05880.61723.650.012057MVO0.052510.3760210.335130.012790


## 4.2.3. Three-Bar Truss Design Problem

Three-bar truss design is a complex problem in the field of civil engineering [49]. The goal of this problem is to achieve the minimum weight in truss design. Figure 10 shows the design of this problem. The formula of this problem can be described as follows.Figure 10
Three-bar truss design problem.Consider(23)x⟶=x1x2=A1A2.Minimize(24)fx⟶=22x1+x2∗l,subject to(25)g1x⟶=2x1+x22x12+2x1x2P−σ≤0,g2x⟶=x22x12+2x1x2P−σ≤0,g3x⟶=12x2+x1P−σ≤0.Variable range is(26)0≤x1,x2≤1,where l=100cm,P=2KN/cm2,andσ=2KN/cm2.Results of HSMSSA for solving the three-bar design problem are listed in Table9, which are compared with SMA, SSA, AO, AOA, WOA, SCA, and MVO. It can be observed that HSMSSA has an excellent ability to solve the problem in confined space.Table 9
Comparison of HSMSSA results with other competitors for the three-bar truss design problem.
AlgorithmOptimum variablesOptimum costx1x2HSMSSA0.788420.40811263.8523SMA0.791090.40022263.8668SSA0.778230.43736263.9363AO0.778990.43541263.9197AOA0.763420.48382264.3549WOA0.783570.42252263.886SCA0.780580.43766264.5463MVO0.754920.51216264.7851


## 4.2.4. Speed Reducer Problem

In this problem [15], the total weight of the reducer is minimized by optimizing seven variables. Figure 11 shows the design of this problem, and the mathematical formula is as follows.Figure 11
Speed reducer problem.Minimize(27)fx⟶=0.7854x1x223.3333x32+14.9334x3−43.0934−1.508x1x62+x72+7.4777x63+x73,subject to(28)g1x⟶=27x1x22x3−1≤0,g2x⟶=397.5x1x22x32−1≤0,g3x⟶=1.93x43x2x3x64−1≤0,g4x⟶=1.93x53x2x3x74−1≤0,g5x⟶=745x4/x2x32+16.9×106110.0x63−1≤0,g6x⟶=745x4/x2x32+157.5×10685.0x63−1≤0,g7x⟶=x2x340−1≤0,g8x⟶=5x2x1−1≤0,g9x⟶=x112x2−1≤0,g10x⟶=1.5x6+1.9x4−1≤0,g11x⟶=1.1x7+1.9x5−1≤0.Variable range is(29)2.6≤x1≤3.6,0.7≤x2≤0.8,17≤x3≤28,7.3≤x4≤8.3,7.8≤x5≤8.3,2.9≤x6≤3.9,5.0≤x7≤5.5.The comparison results are listed in Table10, which shows the advantage of HSMSSA in realizing the minimum total weight of the problem.Table 10
Comparison of HSMSSA results with other competitors for the speed reducer design problem.
AlgorithmOptimum variablesOptimum weightx1x2x3x4x5x6x7HSMSSA3.49760.7177.37.83.350065.285532995.4374SMA3.514430.7177.324447.805273.352355.284943002.5941SSA3.497670.7177.877978.094013.39435.285743018.644AO3.51380.7177.414617.812913.377095.284573009.9097AOA3.559890.7177.499978.33.46235.285123061.8731WOA3.497390.7177.870398.07693.453095.28543034.0596SCA3.60.7177.38.33.379325.27993063.9102MVO3.60.7178.38.33.382765.360413111.6609


## 4.2.5. Cantilever Beam Design

Cantilever beam design is a type of concrete engineering problem. This problem aims to determine the minimal total weight of the cantilever beam by optimizing the hollow square cross-section parameters [24]. Figure 12 illustrates the design of this problem, and the mathematical described is as follows.Figure 12
Cantilever beam design [24].Consider(30)x=x1x2x3x4x5.Minimize(31)fx⟶=0.6224x1+x2+x3+x4+x5,subject to(32)gx⟶=60x13+27x23+19x33+7x43+1x53−1≤0.Variable range is as follows:0.01≤x1,x2,x3,x4,x5≤100.The results are shown in Table11. From this table, we can see that the performance of HSMSSA is better than all other algorithms and the obtained total weight is minimized.Table 11
Comparison of HSMSSA results with other competitors for the cantilever beam design problem.
AlgorithmOptimum variablesOptimum weightx1x2x3x4x5HSMSSA5.99795.32894.49993.48172.1661.3400SMA6.08465.23094.56023.47252.13491.3405SSA6.10565.10014.3033.73652.31831.3456AO5.82365.38134.51783.47682.31591.3426AOA5.70145.37224.65773.66432.15511.3448WOA5.84925.05874.89173.6822.15641.3469SCA5.91655.91874.51333.31671.98231.3508MVO7.10885.01614.22033.40642.1051.3647As a summary, this section demonstrates the superiority of the proposed HSMSSA algorithm in different characteristics and real case studies. HSMSSA is able to outperform the basic SMA and SSA and other well-known algorithms with very competitive results, which are derived from the robust exploration and exploitation capabilities of HSMSSA. Excellent performance in solving industrial engineering design problems indicates that HSMSSA can be widely used in real-world optimization problems.


## 5. Conclusion

In this paper, a Hybrid Slime Mould Salp Swarm Algorithm (HSMSSA) is proposed by combining the whole SMA as leaders and the exploitation phase of SSA as followers. At the same time, two strategies, including Levy flight and mutation opposition-based learning, are incorporated to enhance the capabilities of exploration and exploitation of HSMSSA. The 23 standard benchmark functions are utilized to evaluate this algorithm for analyzing its exploration, exploitation, and local optima avoidance capabilities. The experimental results show competitive advantages compared to other state-of-the-art metaheuristic algorithms, proving that HSMSSA has better performance than others. Five engineering design problems are solved as well to verify the superiority of the algorithm further, and the results are also very competitive with other metaheuristic algorithms.The proposed HSMSSA can produce very effective results for complex benchmark functions and constrained engineering problems. In the future, HSMSSA can be applied to real-world optimization problems such as multiobjective problems, feature selection, multithresholding image segmentation, convolution neural network, or any problem that belongs to NP-complete or NP-hard problems.



---

*Source: 6379469-2021-09-07.xml*