Department of Civil Engineering, Faculty of Engineering, University of Zanjan, Zanjan 45371-38791, Iran
mbabaei@znu.ac.ir
Show less
History+
Received
Accepted
Published
2015-05-14
2016-08-20
2016-11-29
Issue Date
Revised Date
2016-11-02
PDF
(685KB)
Abstract
In this article, multi-objective optimization of braced frames is investigated using a novel hybrid algorithm. Initially, the applied evolutionary algorithms, ant colony optimization (ACO) and genetic algorithm (GA) are reviewed, followed by developing the hybrid method. A dynamic hybridization of GA and ACO is proposed as a novel hybrid method which does not appear in the literature for optimal design of steel braced frames. Not only the cross section of the beams, columns and braces are considered to be the design variables, but also the topologies of the braces are taken into account as additional design variables. The hybrid algorithm explores the whole design space for optimum solutions. Weight and maximum displacement of the structure are employed as the objective functions for multi-objective optimal design. Subsequently, using the weighted sum method (WSM), the two objective problem are converted to a single objective optimization problem and the proposed hybrid genetic ant colony algorithm (HGAC) is developed for optimal design. Assuming different combination for weight coefficients, a trade-off between the two objectives are obtained in the numerical example section. To make the final decision easier for designers, related constraint is applied to obtain practical topologies. The achieved results show the capability of HGAC to find optimal topologies and sections for the elements.
Mehdi BABAEI, Ebrahim SANAEI.
Multi-objective optimal design of braced frames using hybrid genetic and ant colony optimization.
Front. Struct. Civ. Eng., 2016, 10(4): 472-480 DOI:10.1007/s11709-016-0368-4
The aim of optimal design is to obtain the best solution based on the objective function in the given condition. Structural optimization is a complex, not simply solved, problem. If the layout of the structure is assumed, under the loading and boundary conditions, the structural responses could be obtained using the finite element analysis. The purpose of the optimization, however, is to attain an optimal structural system under the constraints. Since, the structural optimization problems usually include nonlinear and implicit functions, the analytical methods are unable to solve these practical problems. Therefore, numerical methods remain to be the only option in solving practical real optimization problems.
Three types of optimal designs in structural optimization include topology, shape and size optimization. The aim of topology optimization is to determine the optimal connection and relation of structural elements, which define the existence or non-existence of the elements. Shape optimization is the process of looking for an optimum geometry for the structure, and optimal sizing is focused on the determination of optimal cross sections for the structural elements. The proposed hybrid genetic ant colony algorithm (HGAC) in this paper is to find optimal topology and sizing for steel braced frames.
During the recent decades, there has been a tendency to develop heuristic methods to solve discrete optimization problems, which have been inspired from the nature. Although in these methods, there is no guarantee for feasibility and optimality of the solutions, however, they are able to obtain near optimal solutions, which are not possible to attain by conventional gradient-based methods. The developed heuristic methods for optimal design of structures consist of genetic algorithm (GA) [ 1, 2], evolutionary programming (EP) [ 3], simulated annealing (SA) [ 4], evolution strategies (ESs) [ 5], ant colonies (AC) [ 6, 7], particle swarm optimization (PSO) [ 8], tabu search (TS) [ 9], and cellular automata (CA) [ 10]. Elbeltagi et al. [ 11] presented a comparative study on five evolutionary heuristic optimization algorithms in 2005. GAs, among the heuristic algorithms, are most famous for structural optimal design and have a great potential in producing optimal solutions for practical problems [ 12, 13].
In 2005, Espinoza et al. demonstrated that hybrid GAs (HGAs) with local search algorithms required significantly fewer function evaluations to achieve convergence when compared with the simple GA [ 14]. As stated in the literature, GA has strong ability to exploit using crossover and mutation operators, while ant colony (ACO) is more famous to explore providing diverse population [ 15, 16]. Chan and Wong [ 17] proposed a hybrid algorithm based on GA and optimality criteria (OC) for structural topology and size optimization of tall steel frames. Chen et al. proposed a co-evolutionary model based on dynamic combination of GA and ACO [ 18]. Kicinger and Arciszewski [ 19] reported that HGAs have formed a rapidly growing research area with great potential in solving structural design problems. They studied empirical analysis of memetic algorithms for conceptual design of tall buildings. Kareem et al. reported a strong work by utilizing computational fluid dynamics to optimize tall building design [ 20]. Spence and Kareem [ 21] investigated data-enabled design and optimization for tall buildings. Performance-based topology optimization for wind-sensitive tall buildings has been studied and successive results reported [ 22].
Practical optimal topology for reinforced concrete moment resisting frames (RCMRFs) is reported in the literature [ 23]. Similar studies are reported in the literature for steel moment resisting frames with and without bracings [ 24, 25]. Kaveh et al. worked on plastic analyses of frames using genetic algorithm and ant colony algorithm [ 26]. Babaei et al. studied multi-objective optimization of RCMRFs using a non-dominated sorting genetic algorithm II (NSGA-II) [ 27]. To solve quadratic assignment problems (QAP), Xu et al. [ 28] proposed a novel adaptive GA-ACO-local search hybrid algorithm, which is a synergy between the main algorithms. More recently, optimization of tall buildings with outrigger belt trusses were investigated and optimum location and number for trusses were proposed based on two objective functions [ 30].
The following section defines the optimal design problem formulation, to review the application of GA and ACO algorithms, and to develop the proposed novel hybrid algorithm. Section 3 includes case studies to illustrate the capability of the HGAC method for topology and size optimization of steel braced frames. Significant conclusions are drawn in Section 4.
Formulation and optimization method
A minimization problem obtains the design variables so that the objective function is minimized by satisfying all the defined constraints. In general, constrained minimization problems are formulated as follows:
where is objective function; is the design variable vector of N members; and M is the number of constraints.
In multi-objective minimization problems, the design variables are obtained based on the different combinations of the objective functions. Typically, a multi-objective minimization problem could be defined as follows:
where are objective functions; n is the number of objective functions; is the design variable vector of N members; and M is the number of constraints.
There are many different methods to convert a multi-objective optimization problem into a single objective optimization problem [ 31], the most popular one of which is the weighted sum method (WSM) [ 32]. In this method, the main objective function is obtained by a linear combination of all weighted objective functions as follows:
where is a nonzero weight vector for objective functions, defined so that . If all the weights are positive, then the minimum value of the main objective function is Pareto front [ 33].
The objective function for multi-objective optimal design of steel braced frames using weighted sum method could be written as:
The stress constraints are as:
In addition, the displacements constraints are:
where and are the total weight and displacement objective function of the structure, respectively; and are the minimum weight and displacement objective functions obtained by assigning the weakest and the strongest cross sections to the elements, respectively; and are the existing and allowable bending stresses, respectively; and are the story and allowable deflections; and is the density of the structural elements.
The penalized objective function is as follows:
where and are the penalty value and exponent for weight and displacement objective functions, respectively.
Assuming and , (since ), the objective function may be written as:
Since the two objective functions are not in the same scale, they must be normalized properly. Therefore, the normalized compound objective function could be written as follows:
where and are the normalized total weight and displacement objective functions of the structure, respectively.
Genetic algorithm is one of the stochastic search methods, which is inspired by the nature. GAs have many applications in solving optimization problems. In the nature, better generations are achieved by combination of better chromosomes. Sometimes, during the combinations, mutation in chromosomes occurs, which may result in a better chromosome. There are, usually, some operators for selection, combination and mutation of chromosomes, however, GAs utilizes four operators of fitness, selection, crossover and mutation to solve the problems.
GAs, initially, create a number of solutions, renowned as first generation/population, stochastically or systematically. Each solution is named a chromosome. Subsequently, new generations are created, using selection, crossover and mutation operators. Finally, the old population is combined and replaced by some or all of the new population.
Ant system-rank (AS-rank) algorithm is the most famous method of ant colony algorithms, introduced by Bullnheimer et al. [ 34]. In this algorithm, each ant trails a pheromone base on its rank and the best ant contains the largest value of pheromone, similar to elite ant system (EAS). The pheromone updating formula is can be written as:
where is the intensity of trial; is the evaporation constant between 0 and 1; is the trial change for the ranked ants;
is the extra trial for best, elite, solution.
Since better performance of GA and ACO is represented in the literature [ 34], they are employed as the main algorithms of this research. Many different combinations of these methods are developed [ 35]. In this paper, a dynamic hybridization of GA and ACO is developed and applied to solve multi-objective optimal design of steel braced frames, which has not applied in the literature for steel braced frames optimization. The phenomenon of long time convergence, precocity, and stagnation sometimes emerge from GA and ACO. The status of execution process in each algorithm is different. In the other words, in the executing procedure stage, best chromosome, in GA, is different from best ant. Therefore, a hybrid algorithm is presented to run ACO and then switches to GA by evaluating the running state.
According to the code requirements and engineering experiences, in a frame structure the elements of lower levels are stronger than the upper levels. Although it depends on the loading condition and the topology of the structure, but in a building frame structure, beam, column or even brace sections in the first floor need to be the heaviest, while the sections at the top level of the structure need to be the lightest profiles. In this article, sorting of the element sections based on their section number is proposed and applied in the developed algorithm, as a new operator. Many tested examples and their results have shown the effectiveness of this new operator, as illustrated in experimental section.
The stagnation and slow convergence phenomenon emerges in GA and ACO. To overcome this phenomenon during the optimization process, the following definition is applied in the hybrid algorithm:
where is the generation of GA or ACO; is the average fitness in GA or ACO in the th generation; is the fitness of the th chromosome in GA or ant in ACO; is the number of chromosomes in GA or ants in ACO. The defined convergence will become small when the speed convergence tends to slow or stop.
HGAC starts by running ACO. By checking the status of the process using Eqs. (10) and 11, the algorithm switches to the other main approach and continues the process until the termination criteria met. During the process, the new operator searches for the better solution alongside the other operators. The penalized objective function is introduced and employed during the optimization process. However, not only the penalized total weight of the structure is obtained, but the total weight of the frame is also obtained to compare the results.
Numerical results and discussion
Ten-storey braced frame
In this case study, multi-objective topology and size optimization of a ten-storey three-span braced frame is considered to evaluate the developed hybrid algorithm, as shown in Fig. 1. Braces may be added to the middle span. Span length and storey height of the frame are set to 5 and 3.5 m, respectively. Uniform distributed gravity loads on the beams are considered to be 50kN/m and a lateral point load of 30 kN is applied at each story. Topology of the braces and cross section of the columns, beams and braces are defined to be design variables.
Two objective functions of the total weight and the maximum roof displacement are converted to a single objective function using Eq. (8). Columns and beams sections may be selected between W4 × 13 and W44 × 335 AISC W-shaped sections. Element groupings are illustrated in Fig. 2, to consider the symmetry for the frame. There are 10 topological and 10 cross sectional design variables for braces, and 40 cross sectional variables for beams and columns. As a whole, a total of 60 design variables are considered in this example.
The termination criterion is applied when the number of iteration is 100 or when the results of 10 consecutive iterations are similar. Since two objective functions are employed, there is not one specific optimal solution for this problem and trade-off between two objectives are obtained using WSM, as illustrated in Fig. 5. As shown in this figure, a large value of results in for the structural weight, the first objective function, to be important throughout the optimization procedure and decrease the number of braces, but it leads to an increase in the displacement of the best solution, and vice versa. When , the displacement of the frame is very small, however the weight of the structure increases drastically to control the displacement.
The results of the three different scenarios, assuming , are given in Table 1. The effect of this coefficient on the trade-off between objective functions is clear. Since all constraints are satisfied, both the penalized weight and the net weight of the structure are the same. Figure 3 shows the convergence history for HGAC, GA, and ACO, and Fig. 4 illustrates the effect of sorting elements cross section on the HGAC convergence. Moreover, Fig. 5 displays Pareto front, illustrating the topology of different scenarios. This figure shows that to optimize the weight of the structure, the frame does not need more braces, while to minimize the displacement more floors (at least 60 percent) need bracing elements. The best final topology of the structure, when the total weight of the frame is minimized by satisfying the displacement constraint, is shown in Fig. 6.
Forty-storey braced frame
In this example, multi-objective topology and sizing optimization of a 40-storey three-span braced frame is considered. The topology and loading conditions are given in Fig. 7. Span length, storey heights, bracing locations, elements grouping and uniform distributed gravity loads on the beams are the same as example 3.1. Topology of the braces and cross section of the columns, beams and braces are defined to be the design variables. Two objective functions of the total weight and the maximum roof displacement are converted to a single objective function using Eq. (8). Columns and beams sections may be selected between W4 × 13 and W44 × 335, AISC W-shaped sections. There are 40 topological and 40 cross sectional design variables for braces, 160 cross sectional variables for beams and columns, adding up to a total of 240 design variables, to be considered in this example.
A number of 100 chromosomes are set as the initial population. Mutation rate for topology and elements cross sections are set to be 20 and 5 percent, respectively. To explore design space efficiently in order to find better structural forms, a large value of topology mutation rate is assigned.
Figure 8 shows optimal topologies obtained by HGAC for the five different runs of the algorithm. These results are achieved when displacements are set as constraints during the optimization procedure. Since HGAC is a stochastic evolutionary method, different results may be obtained; however, the results have a small variance. The optimal penalized weight and the required bracings (in percentage) for every run are given in Table 2. At least 42 percent of the stories are needed to be strengthened by bracings to control displacements. The convergence history of HGAC for the five different runs and the average of these runs are illustrated in Fig. 9.
It is difficult for designers to make a decision on selecting the best topology based on the obtained forms in Fig. 8, which are not practical. To overcome this problem, another constraint is imposed in the optimization problem to find practical topologies. In this case, the algorithm does not allow a chromosome/ant with a separated bracing layout to be selected as an individual during the optimization procedure. Only ordered bracing patterns remain in the population. The total weight and the maximum displacement are converted to a single objective function using different weights of 0.50 and 0.00. When is set to be 0.50, the optimal topologies are obtained assuming the same importance for two objective functions, and both objectives are minimized. In the case of , only the maximum displacement is minimized. The optimal topologies for the second scenario are drawn in Fig. 10 for the five runs. The optimal penalized weight and the number of the required bracings for every run of these two scenarios are given in Table 3. As shown in this table, almost similar results are obtained for a specific scenario.
Conclusions
In this paper a novel hybrid genetic ant colony (HGAC) algorithm is developed for optimal design of steel braced frames. The topology and sizing of the braces along with sizing of the frame elements are considered as design variables and they are formulated to apply optimization method to the numerical examples. For multi-objective optimization, the total weight and the maximum lateral displacement of the frames are taken as the objectives to be minimized. Formulation of the problem is converted to a single-objective optimization procedure, using weighted sum method (WSM). A trade-off between the objective functions are then obtained and shown. The following results were achieved by applying the algorithm to multi-objective optimization of 10 and 40-storey braced frames: a) the algorithm is capable of finding the optimal topologies and sizing; b) the greater the weight coefficient for displacement objective function, the more braces it needs to control the displacement; c) similarity in the convergence history results for different runs of HGAC; d) proven capability of the algorithm in fast convergence and achieving solutions better than those obtained by GA and ACO.
HollandJ H. Adaptation in Natural and Artificial Systems. Ann Arbor: University of Michigan Press, 1975
[2]
GoldbergD E. Genetic Slgorithms in Search, Optimization, and Machine Learning. Boston, MA: Addison-Wesley Longman Publication Co, 1989
[3]
FogelG L, OwensA J, WalshM J. Artificial Intelligence through Simulated Evolution. New York: John Wiley, 1996
[4]
KirkpatrickS, GelattC D, VecchiM P. Optimization by Simulated Annealing. Science, 1983, 220(4598): 671–680
[5]
RechenbergI. Evolution Strategy: Optimization of Technical Systems according to the Principle of Biological Evolution. Stuttgart: Frommann-Holzboog, 1973
[6]
DorigoM. Optimization, Learning and Natural Algorithms. Pissertation for the Doctoral Degree. Politecnico di Milano, Italy, 1992
[7]
DorigoM, ManiezzoV, ColorniA. Ant system: Optimization by a colony of cooperating agents. IEEE Transactions on System, Man, and Cybernetics-Part B, 1996, 26(1): 29–41
[8]
KennedyJ, EberhartR. Particle swarm optimization. Proceedings of the International Conference on Neural Networks. 1995, 1942– 1948
[9]
GloverF. Tabu search-Part I. ORSA Journal on Computing, 1989, 1(3): 190–206
[10]
SanaeiE, BabaeiM. Topology optimization of structures using cellular automata with constant strain triangles. International Journal of Civil Engineering, 2012, 10(3): 179–188
[11]
ElbeltagiE, HegazyT, GriersonD E. Comparison among five evolutionary-based optimization algorithms. Advanced Engineering Informatics, 2005, 19(1): 43–53
[12]
KrishnamoorthyC S. Structural optimization in practice: potential applications of genetic algorithms. Structural Engineering and Mechanics, 2001, 11(2): 151–170
[13]
PezeshkS, CampC V. State of the Art on the Use of Genetic Algorithms in Design of Steel Structures. In Recent Advanced in Optimal Structural Design, ASCE (edited by S. Burns), 2002, 1–31
[14]
EspinozaF P, MinskerB S, GoldbergD E. Adaptive hybrid genetic algorithm for groundwater remediation design. Journal of Water Resources Planning and Management, 2005, 131(1): 14–24
[15]
TakahashiR. A hybrid Method of Genetic Algorithms and Ant Colony Optimization to Solve the Traveling Salesman Problem. IEEE International Conference on Machine Learning and Applications, Miami Beach, FL, 2009
[16]
ZukhriZ, PaputunganI V. A hybrid optimization algorithm based on genetic algorithm and ant colony optimization. International Journal of Artificial Intelligence & Applications, 2013, 4(5): 63–75 (IJAIA)
[17]
ChanC M, WongK M. Structural topology and element sizing design optimization of tall steel frameworks using a hybrid OC–GA method. Structural and Multidisciplinary Optimization, 2008, 35(5): 473–488
[18]
ChenM, LuQ. A co-evolutionary model based on dynamic combination of genetic algorithm and ant colony algorithm. In: Proceedings of the 6th International Conference on Parallel and Distributed Computing. Applications and Technologies, 2005, 941–944
[19]
KicingerR, ArciszewskiT. Empirical analysis of memetic algorithms for conceptual design of steel structural systems in tall building. In: Advances in Engineering Structures, Mechanics & Construction, (edited by Pandey M, Wie W C, and Xu L). Netherlands: Springer, 2006, 277–288
[20]
KareemA, SpenceS M J, BernardiniE, BobbyS, WeiD. Using computational fluid dynamics to optimize tall building design. CTBUH Journal Issue III, 2013, 3: 38–43
[21]
SpenceS M J, KareemA. Data-enabled design and optimization (DEDOpt): Tall steel building frameworks. Computers & Structures, 2013, 129: 134–147
[22]
BobbyS, SpenceS M J, BernardiniE, KareemA. Performance-based topology optimization for wind-excited tall buildings: A framework. Engineering Structures, 2014, 74: 242–255
[23]
BabaeiM. Exploring practical optimal topology for reinforced concrete moment resisting frame structures. American Journal of Civil Engineering, 2015, 3(4): 102–106
[24]
BabaeiM, Dadash AmiriJ. Determining the optimal topology for intermediate steel moment resisting frames with eccentric braces in hybrid system. International Journal of Structural Engineering (in press)
[25]
BabaeiM, JabbarM. Optimal intermediate steel moment resisting frames with different spans and story numbers. Journal of Engineering and Applied Sciences (Asian Research Publishing Network), 2016, 10(8): 223–226
[26]
KavehA, JahanshahiM, KhanzadiM. Plastic analysis of frames using genetic algorithm and ant colony algorithm. Asian Journal of Civil Engineering, 2008, 3: 9227–9246
[27]
BabaeiM, MollayiM. Multi-objective optimization of reinforced concrete frames using NSGA-II algorithm, under review
[28]
XuY L, LimM H, OngY S, TangJ. A GA-ACO local search hybrid algorithm for solving quadratic assignment problem. Genetic and Evolutionary Computation Conference, 2006, 599–605
[29]
BabaeiM, AsemaniR, KazemiF. Exploring for optimal number and location of trusses in core and outrigger belt truss system In proceedings of the 1st International & 5th National Conference of Steel and Structure, Iranian Association of Steel Structures, Iran, 2015
[30]
AthanT W, PapalambrosP Y. A note on weighted criteria methods for compromise solution multi-objective optimization. Engineering Optimization, 1996, 27(2): 155–176
[31]
KoskiJ, SilvennoinenR. Norm methods and partial weighting in multi-criterion optimization of structures. International Journal for Numerical Methods in Engineering, 1987, 24(6): 1101–1121
[32]
MarlerR T, AroraJ S. Survey of multi-objective optimization methods for engineering. Structural and Multidisciplinary Optimization, 2004, 26(6): 369–395
[33]
BullnheimerB, HartlR F, StraussC. An improved ant system algorithm for the vehicle routing problem. European Journal of Operational Research, 1999, 89: 319–328
[34]
BabaeiM. Multi-Objective Optimization of Mid and High Rise Steel Structures using an Efficient Hybrid Evolutionary Algorithm. Dissertation for the Doctoral Degree. Iran University of Science and Technology, 2012
RIGHTS & PERMISSIONS
Higher Education Press and Springer-Verlag Berlin Heidelberg
AI Summary 中Eng×
Note: Please be aware that the following content is generated by artificial intelligence. This website is not responsible for any consequences arising from the use of this content.