Recently developed chaotic charged system search was combined to feasible-based method to solve constraint engineering optimization problems. Using chaotic maps into the CSS increases the global search mobility for a better global optimization. In the present method, an improved feasible-based method is utilized to handle the constraints. Some constraint design examples are tested using the new chaotic-based methods, and the results are compared to recognize the most efficient and powerful algorithm.
1. Introduction
The charged system search (CSS) is known as one of the efficient optimization algorithms among metaheuristic algorithms in recent years. The large number of researches on applying and enhancing this method proves this [1–11]. This algorithm uses the governing laws of electrostatics in physics and the governing laws of motion from the Newtonian mechanics [1].
In physics, the space surrounding an electric charge has a property known as the electric field. This field exerts a force on other electrically charged objects. The electric field surrounding a point charge is specified by Coulomb’s law. Coulomb confirmed that the electric force between any two small charged spheres is inversely proportional to the square of the separation distance between the particles directed along the line joining them and proportional to the product of the charges of the two particles. Also, the magnitude of the electric field at a point inside a charged sphere can be obtained using Gauss’s law that is proportional to the separation distance between the particles. In the CSS, charged particles (CPs) or solution candidates are treated as a charged sphere that can exert electrical forces on each other according to the Coulomb and Gauss laws of electrostatics. The resultant force acts on each CP creating an acceleration according to Newton's second law by which, in combination with Newtonian mechanics, the position of each CP can be determined [1].
Recently, Talatahari et al. [7] developed chaotic charged system search algorithms (CCSS) based on the idea of using chaotic systems instead of random processes. In the CCSS algorithms, the role of randomness can be played by a chaotic dynamics. Experimental studies show the benefits of using chaotic signals instead of random signals. Nine chaos-based CSS methods were developed, and then, for each variant, the performance of ten different chaotic maps was investigated to identify the most powerful variant. As suggested in [7], here more elaborated experiments are performed to discover the adaptive algorithm for solving optimization problems. Also, this paper combines the most powerful CCSS algorithms with feasible-based method in order to solve constraint engineering design problems. In order to evaluate these algorithms, some well-studied numerical engineering examples are analyzed using some best CCSS methods and compared to the other methods.
2. Standard and Chaotic-Based Charged System Search Algorithm2.1. Standard Charged System Search Algorithm
The charged system search contains a number of charged particle (CPs) where each one is treated as a charged sphere and can insert an electric force to the others. The pseudocode for the CSS algorithm is summarized as follows [1].
Step 1 (initialization).
The magnitude of charge for each CP is defined as
(1)qi=fit(i)-fitworstfitbest-fitworsti=1,2,…,N,
where fitbest and fitworst are the best and the worst fitness of all the CPs, fit(i) represents the fitness of the agent i, and N is the total number of CPs. The separation distance rij between two charged particles is defined as follows:
(2)rij=∥Xi-Xj∥∥(Xi+Xj)/2-Xbest∥+ε,
where Xi and Xj are the positions of the ith and jth CPs, respectively, Xbest is the position of the best current CP, and ε is a small positive number. The initial positions of CPs are determined randomly as
(3)xi,j(o)=xi,min+randij·(xi,max-xi,min),i=1,2,…,N,
where xi,j(o) determines the initial value of the ith variable for the jth CP; xi,min and xi,max are the minimum and the maximum allowable values for the ith variable.
Step 2 (CM creation).
A number of the best CPs and the values of their corresponding fitness functions are saved in the charged memory (CM).
Step 3 (the probability function and kind of the forces determination).
The probability of moving each CP toward the others is determined using the following function:
(4)pij={1,fit(i)-fitbestfit(j)-fit(i)>rand∨fit(j)>fit(i),0,otherwise.
The kind of the forces can be attractive or repelling determined by using the force parameter arij defined as
(5)arij={+1,kt<randij,-1,kt>randij,
where arij determines the type of the force, in which +1 represents the attractive force and −1 denotes the repelling force, and kt is a parameter to control the effect of the kind of the force.
Step 4 (forces determination).
The resultant force vector for each CP is calculated as
(6)Fj=qj∑i,i≠j(qia3rij·i1+qirij2·i2)arijpij(Xi-Xj),∑i,i≠j(qia3rij·i1+i2i2)j=1,2,…,Ni1=1,i2=0⟺rij<ai1=0,i2=1⟺rij≥a.
Step 5 (solution construction).
Each CP moves to the new position as
(7)Xj,new=randj1·ka·Fjmj·Δt2+randj2·kv·Vj,old·Δt+Xj,old,Vj,new=Xj,new-Xj,oldΔt,
where ka and kv are the acceleration and the velocity coefficients, respectively, and randj1 and randj2 are two random numbers uniformly distributed in the range (0, 1).
Step 6 (CP position correction).
If each CP swerves off the predefined bounds, its position is corrected using the harmony search-based handling approach.
Step 7 (CM updating).
The better new vectors are included in the CM, and the worst ones are excluded from the CM.
Step 8 (terminating criterion control).
Steps 3–7 are repeated until a terminating criterion is satisfied.
2.2. Chaotic Charged System Search Algorithm
To the best of our knowledge, random initialization of the CSS and the adjusted limit parameters may affect the performance of the algorithm and reduce or increase its convergence speed. In other words, the parameters of the algorithm such as kt, ka, and kv are the key factors to control the balance of the exploration and exploitation of the algorithm. This problem was solved by introducing chaotic charged system search [7] where the suitable values or the initial population can be generated chaotically by using chaotic maps. Chaos is a deterministic, random-like process found in nonlinear, dynamical system, which is nonperiod, nonconverging, and bounded [12]. The nature of chaos looks to be random and unpredictable, possessing an element of regularity and randomness of a simple deterministic dynamical system, and chaotic system may be considered as the source of randomness [13, 14].
The chaotic CSS algorithms, denoted by CCSS, can be obtained by using the values generated by a chaotic map instead of one or more random parameters needed in the CSS algorithm. Therefore, considering which parameter is defined chaotically, we can specify different algorithms. Table 1 summarizes different chaotic-based CSS algorithms. There are different chaotic maps to be used as listed in Table 2.
Different chaotic-based CSS algorithms.
Algorithm
Equation
Condition
CCSS-1
xi,j=xi,min+cmij·(xi,max-xi,min)
The initial positions of CPs are determined chaotically.
In [7], by hybridizing the charged system search and chaos, different methods were developed to solve numerical global optimization problems. The proposed approaches utilized different chaotic maps to adapt the parameters of the CSS algorithm. Nine chaotic CSS algorithms by using ten different chaotic maps are defined resulting in 90 different methods. These methods were then analyzed for the benchmark mathematical functions. The simulation resulting demonstrated that some tested CCSS approaches are efficient methods to explore the search space and discover the global solution. From the statistical investigation, it was shown that the CCSS-8, utilizing the chaotic maps, instead of aij, pij, kv, and ka, has better performance than other approaches. The second best methods were CCSS-6 and CCSS-9. In order to simplify determining the most adaptive maps, for each algorithm three maps are chosen as the best ones as summarized in Table 3. In this paper, the selected maps are utilized for solving constraint engineering design examples.
Best maps for each CCSS.
Best maps
CCSS
1
2
3
4
5
6
7
8
9
1
Tent
Tent
Sinusoidal
Tent
Liebovtech
Liebovtech
Liebovtech
Sinus
Sinus
2
Circle
Gauss
Ikeda
Ikeda
Sinus
Sinusoidal
Zaslavski
Liebovtech
Tent
3
Sinus
Ikeda
Sinus/Gauss
Circle
Tent
Zaslavski
Sinus
Tent
Gauss
On the other hand, it is necessary to handle the constraints of the problem by using a suitable method [20–22]. The feasible-based approach [23] deals with constrained search spaces by using the separation of constraints and objectives. In this method, the idea is to avoid the combination of the value of the objective function and the constraints of a problem to assign fitness, like when using a penalty function [23–26]. Kaveh and Talatahari [27] have presented a modified feasible-based method which employs the following four rules as.
Rule 1.
Any feasible solution is preferred to any infeasible solution.
Rule 2.
Infeasible solutions with slight violations of the constraints are treated as feasible ones.
Rule 3.
Between two feasible solutions, the one with better objective function value is preferred.
Rule 4.
Between two infeasible solutions, the one having smaller sum of constraint violations is preferred.
By using the first and fourth rules, the search tends to the feasible region rather than infeasible region, and by employing the third rule, the search tends to the feasible region with good solutions [28]. For most of the engineering optimization problems, the global minimum locates on or close to the boundary of a feasible design space. By applying Rule 2, the CPs approach to the boundaries and can go near to the global minimum in a great probability.
If the location of CPs become out of the variable boundaries, the solutions cannot be used. In this paper, using the harmony search-based handling approach, this problem is dealt with. According to this mechanism, any component of the CP’s vector violating the variable boundaries can be generated randomly from CM as
(8)xi,j={w.p.HMCR==>selectanewvalueforvariablefromCM==>w.p.(1-PAR)donothing==>w.p.PARchooseaneighboringvaluew.p.(1-HMCR)==>selectanewvaluerandomlyfromtheallowablelist.
4. Numerical Engineering Problems
The population size N=20 is sufficient for most of the problems while the maximum number of function evolutions for the examples is set to 4,000 (200 iterations). The explanation of the examples is presented in Section 4.1, and the performance of the CCSS algorithms to optimize these functions is investigated in the next subsection.
4.1. Description of the Examples
Three engineering design problems which have been previously solved using a variety of other techniques are considered to perform investigation on efficiency of the proposed algorithms. The description of these examples is as the following.
4.1.1. A Tension/Compression String Design Problem
This problem consists of minimizing the weight of a tension/compression spring subject to constraints on shear stress, surge frequency, and minimum deflection as shown in Figure 1. The design variables are the mean coil diameter D(=x1), the wire diameter d(=x2), and the number of active coils N(=x3). The problem is defined with details in [1].
Schematic of tension/compression string.
4.1.2. A Welded Beam Design Problem
The welded beam structure, shown in Figure 2, is a practical design problem that often has been used as a benchmark problem for testing different optimization methods. The objective is to find the minimum fabricating cost of the welded beam subject to constraints on shear stress (τ), bending stress (σ), buckling load (Pc), end deflection (δ), and side constraint. There are four design variables, namely, h(=x1), l(=x2), t(=x3), and b(=x4). The detailed information about the constraints and objective function is presented in [1].
Schematic of welded beam structure.
4.1.3. A Pressure Vessel Design Problem
A cylindrical vessel is capped at both ends by hemispherical heads as shown in Figure 3. The objective is to minimize the total cost, including the cost of material, forming, and welding.
Schematic of pressure vessel.
The variables contain x1 as the thickness of the shell (Ts), x2 as the thickness of the head (Th), x3 as the inner radius (R), and x4 as the length of cylindrical section of the vessel, not including the head (L). Details related to the objective and constraint functions are available in [1].
4.2. Experimental Results4.2.1. Results for the Tension/Compression String Problem
For this example, performances are assessed on the basis of the best fitness values and the statistics results of the new approaches as reported in Figure 4. Simulation results show that all the proposed methods perform satisfactorily. However, the best feasible solutions obtained by the methods CCSS-8 (with Sinus map), CCSS-6 (with Zaslavski map), and CCSS-5 (with Sinus map) are the best results among the others. Meanwhile, the CSS-8 and CCSS-6 methods (almost with all the three maps) have better performance in relation to the mean and worst values. Moreover, the standard deviation of the results by the CSS-6 in 50 independent runs for this problem is the smallest among the others.
The results obtained by the CCSS methods for the string problem: (a) best results (b) mean of results, (c) worst result, and (d) Standard deviation.
4.2.2. Results for the Welded Beam Problem
Figure 5 presents the statistical information of this problem obtained by the proposed algorithms. Clearly, the worst solutions belong to the CCSS-7, while the best ones belong to the CCSS-6 and CCSS-8. The CCSS-4 as well as CCSS-9 stand in the second place. Although some new methods have small differences in relation to best, mean, or worst results, obviously many of them perform very well in relation to the standard deviation. Similar to the previous example, the CCSS-6 has the best standard deviation compared to the other methods.
The results obtained by the CCSS methods for the welded beam problem: (a) best results (b) mean of results, (c) worst result, (d) Standard deviation.
4.2.3. Results for the Pressure Vessel Problem
The obtained results using 9 variants of the presented algorithms are shown in Figure 6. Similar to the previous examples, many of the CCSS methods can find a suitable result successfully; however, the best solution found by CCSS-6 and CCSS-4 is better than the best solutions found by the other techniques. Also, it can be seen that the average searching quality of the CCSS-8 (with Sinus and Tent maps) and CCSS-6 (with Sinusoidal map) algorithms is better than those of other methods. Many of the proposed methods improve the reliability of the algorithm by reducing the standard deviation values.
The results obtained by the CCSS methods for the pressure vessel problem: (a) best results (b) mean of results, (c) worst result, and (d) Standard deviation.
5. Conclusion Remarks
As suggested in [7] due to the superiorities of the CCSS methods, here more elaborated experiments are performed to discover the better methods which can be utilized in solving engineering problems. From 90 different chaotic charged system search algorithms (obtained by using 10 chaotic maps and 9 different methods), 27 most efferent methods (3 different chaotic maps for each algorithm) are chosen and investigated to solve constraint engineering problems.
Performances are assessed on the basis of the best fitness values and the statistics results of the new approaches from 50 runs with different seeds. Simulation results show that for all examples, the proposed methods perform satisfactorily. Almost all of the proposed methods improve the reliability of the algorithm by reducing the standard deviation values. From numerical results, it is clear that CCSS-6, in which the coefficient of the velocity is determined chaotically, is the most reliable algorithm having the smallest standard deviation values, while, the algorithm with chaotic coefficients for the force and velocity (CCSS-7) is the worst one. Meanwhile, the CSS-8 and CCSS-6 methods have better performance in relation to the best, mean, and worst values. In CCSS-8, all parameters of the algorithm are determined chaotically, but the initial positions of agents are defined randomly. To sum up, the coefficient of the velocity plays a key role in reliability of the algorithm and the results show that a chaotic velocity coefficient can improve the performance of the algorithm; in addition, using all chaotic parameters of the algorithm can improve the performance, as well; however, chaotic initialization has very small or even no influence on the final results. This chaotic charged model, which is a coupled system, for constraint optimization problems in future could also benefit from other coupled systems such as the kinetic models of competition and the corresponding hybrid competition models [29, 30].
KavehA.TalatahariS.A novel heuristic optimization method: charged system search20102133-4267289KavehA.TalatahariS.A charged system search with a fly to boundary method for discrete optimum design of truss structures20101132772932-s2.0-77957969896KavehA.TalatahariS.Optimal design of skeletal structures via the charged system search algorithm20104168939112-s2.0-7795295200510.1007/s00158-009-0462-5KavehA.TalatahariS.Charged system search for optimum grillage system design using the LRFD-AISC code20106667677712-s2.0-7754908323110.1016/j.jcsr.2010.01.007KavehA.TalatahariS.Geometry and topology optimization of geodesic domes using charged system search20114322152292-s2.0-7995181421210.1007/s00158-010-0566-yKavehA.TalatahariS.An enhanced charged system search for configuration optimization using the concept of fields of forces20114333393512-s2.0-7955169120010.1007/s00158-010-0571-1TalatahariS.KavehA.SheikholeslamiR.An efficient charged system search using chaos for optimization problems201112305325TalatahariS.SheikholeslamiR.ShadfaranM.PorbabaM.Charged system search algorithm for optimum design of gravity retaining walls subject to seismic loading201220121030162810. 1155/2012/301628TalatahariS.KavehA.SheikholeslamiR.Engineering design optimization using chaotic enhanced charged system search algorithms2012223102269228510.1007/s00707-012-0704-2TalatahariS.KavehA.Mohajer RahbariN.Parameter identification of Bouc-Wen model for MR fluid dampers using adaptive charged system search optimization20122682523253410.1007/s12206-012-0625-yNiknamT.GolestanehF.ShafieiM.Probabilistic energy management of a renewable microgrid with hydrogen storage using self-adaptive charge search algorithm20134925226710.1016/j.energy.2012.09.055AlatasB.Chaotic bee colony algorithms for global numerical optimization2010375682568710.1016/j.eswa.2010.02.042SchusterG. G.19882ndWeinheim, GermanyFederal Republic of Germany: Physick, GmnHCoelhoL. D. S.MarianiV. C.Use of chaotic sequences in a biologically inspired algorithm for engineering design optimization2008343190519132-s2.0-3744900751310.1016/j.eswa.2007.02.002MayR. M.Simple mathematical models with very complicated dynamics197626155604594672-s2.0-001718544310.1038/261459a0PeitgenH.-O.JürgensH.SaupeD.1992Berlin, GermanySpringerMR1185709ZBL0785.58001ZhengW. M.Kneading plane of the circle map199447122112332-s2.0-002846351310.1016/0960-0779(94)90033-7ZBL0817.58030DresslerU.FarmerJ. D.Generalized Lyapunov exponents corresponding to higher derivatives199259436537710.1016/0167-2789(92)90076-YMR1192750ZBL0778.58036ErramilliA.SinghR. P.PruthiP.1994Stockholm-Kista, SwedenRoyal Institute of TechnologyCattaniC.ChenS. Y.AldashevG.Information and modeling in complexity20122012486841310.1155/2012/868413LimS. C.EabC. H.MakK. H.LiM.ChenS. Y.Solving linear coupled fractional differential equations by direct operational method and some applications201220122865393910.1155/2012/653939MR2861883LuP.ChenS. Y.ZhengY.Artificial intelligence in civil engineering201220122214597410. 1155/2012/145974DebK.An efficient constraint handling method for genetic algorithms200018631133810.1016/S0045-7825(99)00389-8ZBL1028.90533Coello CoelloC. A.Theoretical and numerical constraint-handling techniques used with evolutionary algorithms: a survey of the state of the art200219111-121245128710.1016/S0045-7825(01)00323-1MR1877684ZBL1026.74056ChenS. Y.LiY. F.Automatic sensor placement for model-based robot vision2004341393408ZhengY.ShiH.ChenS.Fuzzy combinatorial optimization with multiple ranking criteria: a staged tabu search framework201283457472MR2978296KavehA.TalatahariS.Engineering optimization withhybrid particle swarm and ant colony optimization2009106611628MontesE. M.CoelloC. A. C.An empirical study about the usefulness of evolution strategies to solve constrained optimization problems200837444347310.1080/03081070701303470MR2440892ZBL1219.90129ZBL1219.90129CattaniC.CiancioA.Hybrid two scales mathematical tools for active particles modelling complex systems with learning hiding dynamics200717217118710.1142/S0218202507001875MR2292354ZBL1142.82019CattaniC.CiancioA.Separable transition density in the hybrid model for tumor-immune system competition201220126610124MR2874490ZBL1234.92026