Abstract
According to the conditions of today's world, design of resistant structures against blast loading is an important subject that requires special attention. Thus, given the benefits of optimization in engineering, development and assessment of optimization methods for optimum design of structures against blast is of great importance. In this research, the optimum design of steel frame structures against blast loading is investigated. For this purpose first an optimization methodology is proposed. In the proposed method the structural analysis is performed using nonlinear explicit finite element analysis. Based on the proposed method a framework is developed and three numerical examples are investigated using different numerical optimization techniques. Results of this study show that by using nonlinear explicit FE analysis as the structural analysis method and NLPQLP optimization technique as the optimization method, the current optimization problem can be performed effectively, because the procedure is relatively accurate and computationally inexpensive.
Keywords
Optimization; blast; nonlinear analysis; steel frames; SQP; NLPQLP; PSO; MIGA
1. Introduction
In today’s world, terrorist attacks and wars are phenomena which threaten human security all over the world. Meanwhile, structures play an important role in either increasing or decreasing damages or losses. Also as it's true for other types of loads, if a structure is to be designed to have an appropriate performance against potential events and their subsequent loadings, it will bring psychological comfort for residents and the whole society even though they never occur. Thus, it is very important to conduct more studies on the behavior of structures against blast and the resistant design of structures against this type of loading, particularly under current conditions and events in the world.
Preliminary researches on blast loading, behavior and resistant design of structures against blast date back to the years of World War II (Taylor 1940Taylor, G.I. (1940), “Notes on the Dynamics of Shock Waves from Bar Explosive Charges”, UK Ministry of Home Security, Civil Defense Research Committee Paper., ^{1941a}Taylor, G.I. (1941a), “The Propagation and Decay of Blast Waves”, UK Home Office, ARP department., ^{1941b}Taylor, G.I. (1941b), “The Propagation of Blast Waves over the Ground”, UK Ministry of Home Security, Civil Defence Research Committee paper.). But most studies in this field have been conducted in about past 15 years. Some of the most important researches in recent years carried out on steel structures are as follows:
Hadianfard et al. (2012Hadianfard, M. A., Farahani, A., Jahromi, A. (2012). “On the effect of steel columns cross sectional properties on the behaviours when subjected to blast loading”. Structural Engineering and Mechanics, 44(4), 449463.) studied the effect of steel columns cross sectional properties on their behaviors when subjected to blast. Using ANSYS software, they analyzed some steel columns with different shapes of crosssection and different boundary conditions, subjected to blast loading. They concluded that shape and elasticplastic properties of sections and also boundary conditions of columns play important roles on the response of steel columns subjected to blast.
Nassr et al. (2012Nassr, A.A., Razaqpur, A.Gh., Tait, M.J., Campidelli, M., Foo, S. (2012), “Single and multidegree of freedom analysis of steel beams under blast loading”, Nuclear Engineering and Design, 242, 6377.) modeled steel beam and beamcolumns against blast load using single and multidegree of freedom models. First, they conducted an experimental study to evaluate responses of some wide flange beams subjected to blast. Then, they compared the results of the SDOF and MDOF models with the experimental tests. Based on the results, it was shown that both proposed single and multidegree of freedom models could predict the history of responses with a good accuracy. Also, they concluded that the use of a constant strain rate to calculate Dynamic Increase Factor (DIF) would lead to a conservative design.
Using passive unidirectional dampers, Saeed Monir (2013Saeed Monir, H. (2013), “Flexible blast resisitant steel structures by using unidirectional passive dampers”, Journal of Constructional Steel Research, 90, 98107.) conducted a research on the resistant steel structures against blast. He presented a new unidirectional passive damper which shows different performance against compression and tension. He concluded that using this type of damper in a ductile frame, a resistant structure would be obtained which it could absorb most of the blast energy.
Nassr et al. (2013Nassr, A.A., Razaqpur, A.Gh., Tait, M.J., Campidelli, M., Foo, S. (2013), “Strength and and stability of steel beam columns under blast load”, International Journal of Impact Engineering, 55, 3448.) evaluated the resistance and stability of steel beamcolumns subjected to blast loading. They used a onedegree of freedom model for studying the effects of axial load on the strength and stability of columns subjected to blast. The model was validated by comparing the results with experimental results and also finite element analysis. Comparing the results of the onedegree of freedom model and results from UFC 334002 code, showed that regardless of axial load to the Euler critical load ratio, the UFC method overestimates the column’s capacity for ductility coefficients greater than one.
Coffield and Adeli (2014Coffield, A., Adeli, H. (2014), “An investigation of the effectiveness of the framing systems in steel structures subjected to blast loading”, Journal of Civil Engineering and Management, 20(6), 767777.) examined the performance of earthquakeresistant bracing systems against blast. They studied three systems including Moments Resistant Frame (MRF), a Centrally Braced Frame (CBF) and an Eccentrically Bracing Frame (EBF), designed for earthquake. The results showed that the CBF system had a better resistance level in the blast scenarios considered in the study.
Elsanadedy et al. (2014Elsanadedy, H.M., Almusallam, T.H., Alharbi, Y.R., AlSalloum, Y.A., Abbas, H. (2014), “Progressive collapse potential of a typical steel building due to blast attacks”, Journal of Constructional Steel Research, 101, 143157.) studied the potential of progressive collapse in steel structures subjected to blast attacks. They analyzed a conventional multistory steel frame against blast to evaluate its vulnerability in accidental or terroristic blast scenarios. Based on the results of finite element analysis, they proposed strategies for reducing and controlling potential progressive collapse in steel structures.
Habibi and Khaledy (2015Habibi, A.R., Khaledy, N. (2015), “Development of an Exact Method to Analyze BeamColumns Subjected to Blast”, Journal of passive deffence, 5(4), 2128.) proposed an analytical method for analyzing beamcolumns against lateral blast load. Assuming elastic behavior for beamcolumn, the governing fourthorder partial differential equation was formulated, and then an analytical method for solving this differential equation was proposed. By studying a numerical example, they concluded that the influence of axial force on lateral displacement of the beamcolumns is very significant and practically this effect should not be neglected. Also, it was concluded that for different values of axial force, the maximum lateral deformation does not occur at a specified time and the time is a function of the applied axial force.
On the other hand, optimum design of structures given numerous advantages such as cost and time savings is of great importance in structural engineering. Due to advances in the fields of computers and processors, optimization science has undergone a great progress. It is always updated along with new situations and needs of the day. Some of the important researches have carried out on the optimization of steel structures in recent years are as follows:
Degertekin (2012Degertekin, S. O. (2012). “Optimum design of geometrically nonlinear steel frames using artificial bee colony algorithm”. Steel and Composite Structures, 12(6), 505522.) studied the optimum design of geometrically nonlinear steel frames using artificial bee colony algorithm (ABC). Structural weight minimization was the objective of the research while strength, displacement and size constraints were taken into account. The optimization was performed considering the geometrical nonlinearity. Three numerical examples where studied in the research using the ABC algorithm. He concluded the ABC algorithm could find better results comparing to the other metaheuristic algorithms.
Kaveh et al. (2010Kaveh, A., Farahmand Azar, B., Hadidi, A., Rezazadeh Sorochi, F., Talatahari, S. (2010), “Performance based Seismic design of steel frames using ant colony optimization”, Journal of Constructional Steel Reseach, 66, 566574.) studied the optimum seismic performance based design of steel frames using Ant Colony algorithm. To get the response of the structure at various performance levels, the nonlinear pushover analysis was performed using a computer. The structural weight was chosen as the objective function, and story drifts were selected as the constraints. In two numerical examples, they showed that the used method is superior to the traditional Genetic Algorithm.
Based on the history of nonlinear responses, Gong et al. (2012Gong, Y., Xue, Y., Xu, L., Grierson, D.E. (2012), “EnergyBased design optimization of steel building frameworks using nonlinear response history analysis”, Journal of Constructional Steel Research, 68, 4350.) presented a method for optimum design of steel frames under seismic loads. Minimization of weight, minimization of input earthquake energy and maximization of energy absorption, were considered as the three objective functions. Also, story drift and plastic hinge rotation in members were considered as design constraints. A three story building was studied as a numerical example. They concluded that the proposed method is an efficient method for designing steel frames under seismic loadings.
Habibi and Rostami (2013Habibi, A.R, Rostami, H. (2013), “Optimum design of plane Steel moment frames by using consistent approximation method”, MSc. Thesis, University of Kurdistan, Sanandaj, Iran.) developed an optimization method for designing steel frames based on Consistent Approximation (CONAP) method. Structural weight was considered as the objective function, and design constraints were considered based on AISC and Iranian steel structures design Code. Having studied multiple numerical examples, they concluded that the proposed method could easily be used to achieve economic and reliable designs. Also, it was shown that the proposed algorithm would converge in a limited number of iterations.
Kaveh and Nasrollahi (2014Kaveh, A., Nasrollahi, A. (2014), “Performancebased seismic design of steel frames utilizing charged system search optimization”, Applied Soft Computing, 22, 213221.) conducted a research on optimum performance based on seismic design of steel frames using charged system search (CSS) algorithm. The analysis of structure was performed using pushover approach assuming semirigid connections. Based on the results, they concluded that the proposed method would significantly reduce the weight of the structure compared to a traditional design.
Kaveh et al. (2015Kaveh, A., FahimiFarzam, M., & KalatehAhani, M. (2015). “Optimum design of steel frame structures considering construction cost and seismic damage”. Smart Structures and Systems, 16(1), 126.) studied the optimum design of steel frames considering construction cost and seismic damage. They used ParkAng damage index for the seismic damage and constraints were considered based on FEMA350 seismic design specifications. They adapted the nondominated sorting genetic algorithm (NSGA II) as the multiobjective optimization technique. They concluded that the proposed framework is effective for achieving convenient pareto front of possible optimal solutions.
Artar (2016Artar, M. (2016). “Optimum design of braced steel frames via teaching learning based optimization”. STEEL AND COMPOSITE STRUCTURES, 22(4), 733744.) studied the optimum design of braced steel frames using teaching learning based optimization (TLBO). The optimization problems were performed based on AISCASD design specifications. Structural weight minimization was the design objective. Design constraints were considered as stress, geometrical size, displacement and interstory drifts. By studying two numerical examples they concluded that the TLBO is a robust and applicable method for optimum design of multielement structures.
Gholizadeh et al. (2017Gholizadeh, S., Davoudi, H., Fattahi, F. (2017) “Design of steel frames by an enhanced moth flame optimization algorithm”, Steel and Composite Structures, 24(1), 129140.) studied optimum design of steel frames using mothflame optimization (MFO) and enhanced mothflame optimization (EMFO) algorithms. By solving some benchmark problems using MFO and EMFO, and comparing the results with results of the other metaheuristic methods, they concluded that the proposed EMFO method had better optimal results while it needed less computational efforts.
In the field of optimum design of structures against blast, researches are just limited to few researches on nonbuilding and nonframed structures (Taha et al. 2009Taha, M.R, ColakAltunc, A.B., AlHalik, M. (2009), “A multiobjective optimization approach for design of blastresistant composite laminates using carbon nanotubes”, Composites: Part B, 40, 522599.; Sun 2011Sun, W. (2011), “Dynamic Response Analysis and Optimal Design of a RC Slab to Blast Loads”, Advanced Materials research, 163167, 23902396., Salimi et al. 2012Salimi, H., Saranjam, B., Hoseini Fard, A., Ahmadzadeh, M. (2012), “Use of Genetic Algorithms for Optimal Design of Sandwich Panels Subjected to Underwater Shock Loading”, Journal of Mechanical Engineering, 58 (3), 156164.; Qi et al. 2013Qi, Ch., Yang, Sh., Yang, L.J., Wei, Zh.Y., Lu, Zh. H. (2013), “Blast resistance and multiobjective optimization of aluminum foamcored sandwich panels, Composite Structures, 105, 4457.; Xia et al. 2015Xia, Y., Wu, Ch., Li, Zh. X. (2015), “Optimized Design of Foam Cladding for Protection of Reinforced Concrete Members under Blast Loading”, ASCE Journal of Structural Engineering, ASCE J. Struct. Eng, 141(9), 17.). Thus, there is a gap on the field to be filled by more researches on this important subject and related topics. The main purpose of the present study is to compare the efficiency of the three optimization techniques for optimum design of steel moment frames under blast, based on a proposed procedure.
2. Blast Loading
To design structures against blast, it is required to be specified the characteristics of blast loading. Blast load is a timehistory loading which occurs in a very short period. Generally, its timehistory diagram is as shown in Fig. 1.
General timehistory diagram for blast loading (UFC 334002, 2008UFC 334002 (2008), “Structures to resist the effects of accidental explosions”, US Department of Defense, Washington (DC), USA)
When an explosion occurs on the earth surface or in front of a structure where blast waves can be reflected, this reflection amplifies the blast loading, and in this case, the effective blast pressure would be
Blast load timehistory with reflection of waves (UFC 334002, 2008UFC 334002 (2008), “Structures to resist the effects of accidental explosions”, US Department of Defense, Washington (DC), USA)
Explosions which occur on the surface of the earth are called Surface bursts (blasts), and explosions which occur in the air are called Air bursts (blasts). The most important parameters associated with blast loading, include maximum effective pressure (
A diagram for calculating blast loading parameters in surface blasts (UFC 334002, 2008UFC 334002 (2008), “Structures to resist the effects of accidental explosions”, US Department of Defense, Washington (DC), USA)
In this graph Z is the scaled distance of blast:
Where R is the distance between blast source and structure, called standoff distance and W is the mass of equivalent TNT of explosive. All of the graphs and equations for calculating blastrelated parameters are set based on TNT. For other types of explosives, the mass must be multiplied by some coefficients in order to be able to use the graphs. For some explosive types, these coefficients are shown in Table 1.
Equivalent TNT mass for some explosives (Cormie et al. 2009Cormie, D., Mays, G., Smith, P. (2009), “Blast Effects on Buildings (Second Edition)”, Thomas Telford Publishing, London, UK.)
In blast events, the effective pressure decreases as a nonlinear function by moving away from the source of the explosion, so that the local effects should be considered in nearrange explosions. Based on ASCE 5911 (2011)ASCE 5911 (2011), “Blast Protection of Buildings”, American Society of Civil Engineers, Reston, VA, USA. if the scaled distance Z be greater than
The quantity and blast characteristics a given structure is designed for, depends on various factors such as the history of terrorist attacks, the importance of the building, ease of access to the building or terroristic target, the number of occupants of the building, distance from the structure and ease of access to the threatening materials (FEMA 452, 2005FEMA 452 (2005), “A How to Guide to Mitigate Potential Terrorist Attacks Against Buildings”, Federal Management Agency, Washington (DC), USA.). Based on these factors, a simplified method is proposed in FEMA 452 for terroristic threats risk assessment. Also, there are some diagrams which can be used to estimate the potential blast capacity of some blast threats. One of these diagrams is shown in Fig. 4. For example, based on this diagram, a sedan car has a capacity of carrying a mass between 100500 lbs. (45226 kg) of TNT.
Explosive environments blast range to effect (FEMA 452, 2005FEMA 452 (2005), “A How to Guide to Mitigate Potential Terrorist Attacks Against Buildings”, Federal Management Agency, Washington (DC), USA.)
3. Nonlinear dynamic analysis
In this research, the structural analysis has been performed by finite element analysis using Abaqus software (Hibbitt et al., 2010Hibbitt, H., Karlsson, B., & Sorensen, P. (2010). “Abaqus analysis user’s manual version 6.10”. Dassault Systèmes Simulia Corp.: Providence, RI, USA.). In the finite element analysis, the nonlinear dynamic analysis is performed by direct integration methods. Direct integration methods can be performed by two different procedures: Implicit approach and explicit approach. In the “implicit” approach, it is required that in every step the structural stiffness matrix be inverted and nonlinear equations be solved. Thus, when the degrees of freedom are high, this method will be computationally expensive as it is required to calculate the inverse of the stiffness matrix and solve the nonlinear equations. In the “explicit” approach, velocity and displacement are calculated based on the known values at the beginning of each time step. Therefore, calculating the inverse of the stiffness matrix is not required. In other words, compared to the implicit method, this method requires less computational effort at each time step. However, the “implicit” method is numerically stable, while in “explicit” approach, time steps should be considered small enough to ensure the stability of the method. Thus, in dynamic problems that occur in a short time such as impact and explosion problems, the explicit method is better and requires less computational effort. Because of the mentioned computational advantages, the explicit method is used in this research.
One of the most popular explicit methods is the central difference method. In this method velocity and accelerations are written as follows (De Borst et al., 2012De Borst, R., Crisfield, M. A., Remmers, J. J., Verhoosel, C. V. (2012) “Nonlinear Finite Element Analysis of Solids and Structures (Second Edition)”, John Wiley & Sons.):
Where
Where, M is the global mass matrix, and
This equation can be solved for displacement in
Then the displacement increment can be calculated as
And the internal forces vector can be derived as follows:
Where
Where,
In the present study, we have modeled the frame structure using B21 Timoshenko beam element (Hibbitt et al., 2010Hibbitt, H., Karlsson, B., & Sorensen, P. (2010). “Abaqus analysis user’s manual version 6.10”. Dassault Systèmes Simulia Corp.: Providence, RI, USA.). Also, material nonlinearity is considered as bilinear elastoplastic steel with isotropic hardening. Also, geometrical nonlinearity is considered in the analysis. Elements sizes (mesh size) are chosen by 1/10 of member's length. Also, strain rate effects are considered in the analysis according to UFC 334002 (2008)UFC 334002 (2008), “Structures to resist the effects of accidental explosions”, US Department of Defense, Washington (DC), USA as shown in Fig.5. To ensure the accuracy of the finite element model, this model has been validated by another study. Nassr (2012Nassr, A.A. (2012), “Experimental and Analytical Study of the Dynamic Response of Steel Beams and Columns to Blast Loading”, PhD Thesis, McMaster University, Canada.) experimentally and analytically studied the response of some beam and beamcolumns against blast. Here, two beams and two beamcolumns which had been studied by Nassr are modeled. In the modeled samples, each member length is 2413 mm and section profiles are W150X24. The shape and dimensions of W150X24 section profile is shown in Fig 6. Yield stress of the steel material is 470 MPa as reported by Nassr. Since the yield stress is 470 MPa, and there is no specific curve for steel with Fy=470 MPa in Fig 5, thus, the DIF data have been interpolated between the curves for ASTM A36 (Fy=250 MPa) and ASTM A514 (Fy=700 MPa) to be used in FE modeling of the experimental tests. The experimental test setup is shown in Fig 7. Also, a schematic drawing of the location of the charge and the samples is shown in Fig 8. Summary of the modeled tests in the present study is shown in Table 2. Complete details of the experimental tests are described by Nassr (2012).
Strain Rate Effect (SRE) on Dynamic Increase Factor (DIF) (UFC 334002, 2008UFC 334002 (2008), “Structures to resist the effects of accidental explosions”, US Department of Defense, Washington (DC), USA)
The experimental test setup (Nassr, 2012Nassr, A.A. (2012), “Experimental and Analytical Study of the Dynamic Response of Steel Beams and Columns to Blast Loading”, PhD Thesis, McMaster University, Canada.)
Schematic drawing of the charge and the samples location (Nassr, 2012Nassr, A.A. (2012), “Experimental and Analytical Study of the Dynamic Response of Steel Beams and Columns to Blast Loading”, PhD Thesis, McMaster University, Canada.)
In Figs 9 to 12, history of midspan deflection of modeled samples in the current study and the results reported by Nassr are compared. As it is obvious, the used finite element model properly predicted the history of midspan deflection of the members. In Figs 912, differences between the maximum midspan deflection in the current study and the experimental study of Nassr (2012)Nassr, A.A. (2012), “Experimental and Analytical Study of the Dynamic Response of Steel Beams and Columns to Blast Loading”, PhD Thesis, McMaster University, Canada. are 1.4%, 2.8%, 0 and 3.1% respectively. Also in these Figs, it is evident that the history of responses in the current study and the multidegree of freedom model used by Nassr have a high compatibility with each other. In both studies, beam elements are used. In the current study, each member is divided into 10 elements, and the Timoshenko beam elements are used. Nassr had used Bernoulli beam elements, and each member was divided into 24 elements. It should be noted that in the process of verification it was observed that the sensitivity of responses was high to changes in the maximum blast pressure and the blast duration. Therefore, it seems that the very small differences in graphs, is due to the very small differences in the blast loading in the experimental test and the numerical model. In the current study, the blast load applied to the members by a uniform distributed timehistory load, which in the experiments, based on the pressures measured by gauges, the pressure distribution had not been completely uniform over the member's length. In addition, in real conditions, due to shape and type of obstacles in the field, and the effect of reflection of blast waves, the actual diagram of blast loading may not be exactly as the diagrams obtained by the equations and curves presented in manuals and design codes, which probably are derived in some specific and simplified conditions. It should be emphasized that the horizontal axes of the verification diagrams are in millisecond, and therefore differences between the compared graphs are very small.
4. Optimum Design problem formulation
In general, optimization is a procedure to find the best solution that it satisfies certain conditions. The general form of an optimization problem is as follows:
Where
Where
Where,
The main convergence criteria in NLPQLP are based on numerically satisfying the KKT (KarushKuhnTucker) optimality conditions. When the KKT conditions are satisfied with a desired tolerance, the algorithm will stop. More details about this technique are available in NLPQLP User Guide (Schittkowski, 2006Schittkowski, K. (2006). “NLPQLP: A Fortran implementation of a sequential quadratic programming algorithm with distributed and nonmonotone line searchuser’s guide”.). It should be noted that in such engineering problems the structural responses (functions) are not explicitly available. Thus, for obtaining derivatives required in the optimization algorithm, sensitivity analysis can be performed. This can be done easily using finite different procedure, as forward or central difference methods. The derivative of a structural response
Also, the central difference formula can be written as follows:
If there are
The Particle Swarm Optimization (PSO) algorithm was firstly developed by Kennedy and Eberhart (1995Kennedy, J, and R Eberhart. (1995). “Particle Swarm Optimization.” In Proceedings of ICNN’95  International Conference on Neural Networks, 4:194248. IEEE.). This algorithm is a population based optimization technique, which is inspired by social behavior of bird flocking or fish schooling. This algorithm has many similarities with evolutionary algorithms such as Genetic Algorithms (GA). But unlike the GA, no evolution operators such as crossover and mutation are used in PSO. In this algorithm, the potential solutions, called particles, fly through the problem space by following the current optimum particles. The best position found by a particle and its neighbors is used to decide the next position of the particle in the next iteration. The condition of a particle is the value for all design variables in the optimization problem and its velocity in the design space, and each move produces a new generation. The particle swarm algorithm starts by generating a set of random particles that make the initial population as the initial design. Also initial velocities are given to the particles. In each iteration of the algorithm, the particle is moved using two best values. The first one is the best position found by the particle (particlebest or
We assume that
Multi Island Genetic Algorithm (MIGA) (Miki et al, 1999Miki, M., Hiroyasu, T., Kaneko, M., & Hatanaka, K. (1999). “A parallel genetic algorithm with distributed environment scheme”. Systems, Man, and Cybernetics, 1999. IEEE SMC'99 Conference Proceedings. 1999 IEEE International Conference on (Vol. 1, pp. 695700). IEEE.) is an improved version of genetic algorithm (GA). The main difference of MIGA from the traditional GA is that in (MIGA) each population of individuals is divided into several subpopulations called “islands.” All of the traditional GA operations including selection, crossover, and mutation are performed on each subpopulation. Some individuals are selected from each island and migrated to different islands periodically during an operation called “migration”. The migration is controlled by two parameters which are “migration interval” and “migration rate”. Migration interval is the number of generations between each migration, and migration rate is the percentage of individuals migrated from each island. More details about MIGA algorithm can be reviewed in (Miki et al., 1999).
The purpose of a structural optimization is to find the design variables for a structural system in order to maximize or minimize the objective function, and to satisfy the design requirements. Since there is a direct correlation between the weight of a structure and the cost of materials and construction, in many structural optimization problems the structural weight is considered as the objective function which must be minimized to achieve the lowest cost. In the current study the structural weight is considered as the objective function which is being minimized during the optimization process. Also in structural problems, depending on the problem, various constraints such as stress, frequency, and deformation may be considered. As blast loads usually are such that the structural responses may go in inelastic zone, appropriate consideration of nonlinear analysis should be taken into account. In an inelastic design it is better not to consider the strength constraints (Gong et al. 2012Gong, Y., Xue, Y., Xu, L., Grierson, D.E. (2012), “EnergyBased design optimization of steel building frameworks using nonlinear response history analysis”, Journal of Constructional Steel Research, 68, 4350.). Zieman et al. (1992Zieman, R.D., Mc Guire, W., Deierlein, G. (1992), “Inelastic limit states design, part I: planar frame studies”, ASCE Journal of Structural Engineering, 118 (9), 25322549. ) showed that an inelastic design could not be used to full advantage if a design is required to satisfy both strength and deformation constraints, as the strength constraints generally prevent the structural member from yielding. Accordingly, in the present study, only the deformation constraints are taken into account. According to the UFC 334002 (2008)UFC 334002 (2008), “Structures to resist the effects of accidental explosions”, US Department of Defense, Washington (DC), USA criteria, design constraints are considered as story drifts and relative support rotation in beams and columns. Based on these criteria, maximum allowable story drift is
Relative support rotation in beam members (UFC 334002, 2008UFC 334002 (2008), “Structures to resist the effects of accidental explosions”, US Department of Defense, Washington (DC), USA)
Relative support rotation in frame members (UFC 334002, 2008UFC 334002 (2008), “Structures to resist the effects of accidental explosions”, US Department of Defense, Washington (DC), USA)
Based on the above considerations the design problem formulation is considered as follows:
Where
Here, based on the discussed assumptions an optimization methodology is proposed and assessed for optimum design of steel frames against blast loading. The flowchart of the proposed method is shown in Fig 15. As described in Fig 15, after numerical modeling of the structure, the model is initialized. This includes applying gravitational and blast loads on the structure and initializing the design variables. In the present research, the crosssectional areas of members are considered as the design variables. The other geometrical properties of the sections such as section depth, width, flange and web thicknesses are computed based on some equations of crosssectional areas. These equations can be easily derived by performing regression analysis on available steel profile sections such as American AISC or European DIN 1025 Ishape sections, or other available sections. In the next step the structural analysis is performed to derive the structural responses. In the present study, the structural analysis is carried out by explicit nonlinear dynamic finite element analysis of the structure, using Abaqus FE analysis software. As previously mentioned, using explicit nonlinear dynamic FE analysis in blast analysis problems, in addition of its high accuracy, also is relatively computationally inexpensive. After deriving the nonlinear responses the optimization problem can be formulated based on UFC334002 criteria. In the next step, by using numerical optimization techniques and based on the design formulation, an optimization step will be performed. In the case of using gradient based optimization techniques, the sensitivity analysis is required to be able to use the technique. As the structural analysis is performed using nonlinear FE analysis, the structural sensitivities can be calculated with a suitable accuracy. More accurate sensitivities result in more accurate and stable optimization procedure in the case of gradient based techniques. The design procedure is repeated until stopping criteria are satisfied. At the end, cross sectional areas and the weight of the structure at the optimum point are reported as the optimum results. According to the proposed methodology a framework is developed and some numerical examples are studied. In the next section three numerical examples are presented by assessing three numerical optimization techniques of NLPQLP, PSO, and MIGA, using the proposed method and the developed framework.
5. Numerical Examples
The following assumptions are made throughout the studied examples:
Three types of loading are considered which are dead, live and blast loads. In all stories dead and live loads are assumed to be
Where DL is the dead load, LL is the live load, and B is the blast load acting on the structure.
Design constraints are considered based on UFC 334002 (2008)UFC 334002 (2008), “Structures to resist the effects of accidental explosions”, US Department of Defense, Washington (DC), USA criteria. These constraints are the maximum story drift and relative support rotation in beams and columns. Based on the UFC criteria, allowable drift is limited to H/25, and the allowable relative support rotation is limited to 2 degrees (0.035 rad). The weight of the structure is considered as the objective function which is minimized during the optimization process. For all members Yield stress is 240
51. Onestory onebay portal frame
In this example, the optimum design of a onestory portal frame is performed. The frame topology is shown in Fig 16. Two groups of sections are considered in the design: one group for columns and another one for beam. These section groups are shown in Fig. 16 using numbers 1 and 2. The frame is optimized using the developed framework by three techniques of NLPQLP, PSO, and MIGA. The initial configuration of the optimization techniques are summarized in Tables 35. Table 6 shows the initial values of the variables and the objective function. Table 7 shows a comparison between the results of the three techniques at the optimum point. Also Convergence history is shown in Fig. 17 for 100 design evaluations. It can be seen that the optimum design obtained by NLPQLP has a weight 26% lighter than the initial design. Also optimum weight obtained by NLPQLP is 3.8% lighter than PSO and 6.12% lighter than MIGA optimum solutions. Total execution time for performing the optimization process using NLPQLP technique is 79% lower than the execution time required for PSO and 77% lower than the execution time required for MIGA. It is obvious that in this example, NLPQLP technique has a faster convergence and needs less computational efforts comparing to MIGA and PSO techniques.
52. Threestory twobay frame
The second example is a threestory twobay frame shown in Fig. 18. Design section groups are categorized into six groups as shown in Fig. 18 by numbers 1 to 6. The optimization procedure is performed using the developed framework by three optimization techniques of NLPQLP, PSO, and MIGA. The initial configuration of the optimization techniques are summarized in Tables 810. Table 11 shows the initial values of the variables and the objective function. Table 12 shows a comparison between the results of the three techniques at the optimum point. Convergence history is shown in Fig.19 for 250 design evaluations. It is evident that the optimum design obtained by NLPQLP has a weight 14% lighter than the initial design. Also optimum weight obtained by NLPQLP is 2.3% lighter than PSO. In this example no feasible design found using MIGA technique through the 250 runs. Total execution time for performing the optimization process using NLPQLP technique is 62% lower than the execution time required for PSO and 61% lower than the execution time required for MIGA. Like the previous example, NLPQLP technique shows faster convergence and better results comparing to the MIGA and PSO techniques.
53. Sixstory threebay frame
The third example is a sixstory threebay frame shown in Fig 20. Design sections are categorized into nine groups as shown in Fig. 20 by numbers 1 to 9. Using the developed framework this frame is optimized by three optimization algorithms of NLPQLP, PSO, and MIGA. The initial configuration of the optimization techniques are summarized in Tables 1315. Table 16 shows the initial values of the design variables and the structural weight. Table 17 shows a comparison between the results of the three techniques at the optimum point. Convergence history is shown in Fig.21 for 500 design evaluations. It is evident that the optimum design obtained by NLPQLP has a weight 30% lighter than the initial design. Also optimum weight obtained by NLPQLP is 8.5% lighter than PSO. For MIGA method, no feasible design was found through the 500 runs. Also, total execution time for performing the optimization process using NLPQLP technique is 67% lower than the execution time required PSO and 65% lower than the execution time required for MIGA. Just like the previous examples, the NLPQLP technique has a faster convergence and better results comparing to MIGA and PSO techniques.
6. Conclusions
In this study, optimum design of steel moment frames under blast loading was studied. For this purpose, first an optimization algorithm was proposed. The structural analysis in the proposed method is performed using nonlinear explicit dynamic finite element analysis. The nonlinear dynamic finite element analysis results in relatively accurate and realistic evaluations of structural responses. Unlike the seismic analysis problems, In the case of blast loading, the explicit dynamic analysis is very computationally effective. Based on the proposed method a framework was developed and three numerical optimization techniques were investigated through the numerical examples. Results of the numerical examples show that the NLPQLP technique was superior to the other studied optimization techniques, in the current blast optimum design problem. The NLPQLP technique had better results than the two other optimization techniques and found better designs in lower number of runs. Also, between the two methods of PSO and MIGA, the PSO method had better results. Also, the computational time of the PSO and the MIGA methods were approximately the same. Results of this research show that by using the nonlinear explicit FE analysis as the structural analysis method and NLPQLP method as the optimization technique, the optimum design of steel frames under blast loading, can be performed very effectively, because the procedure is relatively precise and computationally cost effective. Also, since the NLPQLP resulted in better optimum points comparing to the PSO and MIGA methods, the NLPQLP method is capable of converging to the global optima in the current problem. In the future researches, the effectiveness of other numerical optimization techniques can be studied for the problem of optimization of steel frames under blast loading.
References
 Artar, M. (2016). “Optimum design of braced steel frames via teaching learning based optimization”. STEEL AND COMPOSITE STRUCTURES, 22(4), 733744.
 ASCE 5911 (2011), “Blast Protection of Buildings”, American Society of Civil Engineers, Reston, VA, USA.
 Coffield, A., Adeli, H. (2014), “An investigation of the effectiveness of the framing systems in steel structures subjected to blast loading”, Journal of Civil Engineering and Management, 20(6), 767777.
 Cormie, D., Mays, G., Smith, P. (2009), “Blast Effects on Buildings (Second Edition)”, Thomas Telford Publishing, London, UK.
 De Borst, R., Crisfield, M. A., Remmers, J. J., Verhoosel, C. V. (2012) “Nonlinear Finite Element Analysis of Solids and Structures (Second Edition)”, John Wiley & Sons.
 Degertekin, S. O. (2012). “Optimum design of geometrically nonlinear steel frames using artificial bee colony algorithm”. Steel and Composite Structures, 12(6), 505522.
 Dusenberry, D.O. (2010). “Handbook for BlastResistant Design of Buildings”, Jon Wiley & Sons, Hoboken, NJ.
 Elsanadedy, H.M., Almusallam, T.H., Alharbi, Y.R., AlSalloum, Y.A., Abbas, H. (2014), “Progressive collapse potential of a typical steel building due to blast attacks”, Journal of Constructional Steel Research, 101, 143157.
 FEMA 452 (2005), “A How to Guide to Mitigate Potential Terrorist Attacks Against Buildings”, Federal Management Agency, Washington (DC), USA.
 Gholizadeh, S., Davoudi, H., Fattahi, F. (2017) “Design of steel frames by an enhanced moth flame optimization algorithm”, Steel and Composite Structures, 24(1), 129140.
 Gong, Y., Xue, Y., Xu, L., Grierson, D.E. (2012), “EnergyBased design optimization of steel building frameworks using nonlinear response history analysis”, Journal of Constructional Steel Research, 68, 4350.
 Hibbitt, H., Karlsson, B., & Sorensen, P. (2010). “Abaqus analysis user’s manual version 6.10”. Dassault Systèmes Simulia Corp.: Providence, RI, USA.
 Habibi, A.R, Rostami, H. (2013), “Optimum design of plane Steel moment frames by using consistent approximation method”, MSc. Thesis, University of Kurdistan, Sanandaj, Iran.
 Habibi, A.R., Khaledy, N. (2015), “Development of an Exact Method to Analyze BeamColumns Subjected to Blast”, Journal of passive deffence, 5(4), 2128.
 Hadianfard, M. A., Farahani, A., Jahromi, A. (2012). “On the effect of steel columns cross sectional properties on the behaviours when subjected to blast loading”. Structural Engineering and Mechanics, 44(4), 449463.
 Kaveh, A., FahimiFarzam, M., & KalatehAhani, M. (2015). “Optimum design of steel frame structures considering construction cost and seismic damage”. Smart Structures and Systems, 16(1), 126.
 Kaveh, A., Farahmand Azar, B., Hadidi, A., Rezazadeh Sorochi, F., Talatahari, S. (2010), “Performance based Seismic design of steel frames using ant colony optimization”, Journal of Constructional Steel Reseach, 66, 566574.
 Kaveh, A., Nasrollahi, A. (2014), “Performancebased seismic design of steel frames utilizing charged system search optimization”, Applied Soft Computing, 22, 213221.
 Kennedy, J. (2011). “Particle swarm optimization”. Encyclopedia of machine learning (pp. 760766). Springer US.
 Kennedy, J, and R Eberhart. (1995). “Particle Swarm Optimization.” In Proceedings of ICNN’95  International Conference on Neural Networks, 4:194248. IEEE.
 Miki, M., Hiroyasu, T., Kaneko, M., & Hatanaka, K. (1999). “A parallel genetic algorithm with distributed environment scheme”. Systems, Man, and Cybernetics, 1999. IEEE SMC'99 Conference Proceedings. 1999 IEEE International Conference on (Vol. 1, pp. 695700). IEEE.
 Nassr, A.A., Razaqpur, A.Gh., Tait, M.J., Campidelli, M., Foo, S. (2012), “Single and multidegree of freedom analysis of steel beams under blast loading”, Nuclear Engineering and Design, 242, 6377.
 Nassr, A.A., Razaqpur, A.Gh., Tait, M.J., Campidelli, M., Foo, S. (2013), “Strength and and stability of steel beam columns under blast load”, International Journal of Impact Engineering, 55, 3448.
 Nassr, A.A. (2012), “Experimental and Analytical Study of the Dynamic Response of Steel Beams and Columns to Blast Loading”, PhD Thesis, McMaster University, Canada.
 Qi, Ch., Yang, Sh., Yang, L.J., Wei, Zh.Y., Lu, Zh. H. (2013), “Blast resistance and multiobjective optimization of aluminum foamcored sandwich panels, Composite Structures, 105, 4457.
 Saeed Monir, H. (2013), “Flexible blast resisitant steel structures by using unidirectional passive dampers”, Journal of Constructional Steel Research, 90, 98107.
 Salimi, H., Saranjam, B., Hoseini Fard, A., Ahmadzadeh, M. (2012), “Use of Genetic Algorithms for Optimal Design of Sandwich Panels Subjected to Underwater Shock Loading”, Journal of Mechanical Engineering, 58 (3), 156164.
 Schittkowski, K. (2006). “NLPQLP: A Fortran implementation of a sequential quadratic programming algorithm with distributed and nonmonotone line searchuser’s guide”.
 Shi, Y., Eberhart, R. (1998). “A modified particle swarm optimizer. Evolutionary Computation Proceedings”. IEEE World Congress on Computational Intelligence. The 1998 IEEE International Conference, 6973.
 Sun, W. (2011), “Dynamic Response Analysis and Optimal Design of a RC Slab to Blast Loads”, Advanced Materials research, 163167, 23902396.
 Taha, M.R, ColakAltunc, A.B., AlHalik, M. (2009), “A multiobjective optimization approach for design of blastresistant composite laminates using carbon nanotubes”, Composites: Part B, 40, 522599.
 Taylor, G.I. (1940), “Notes on the Dynamics of Shock Waves from Bar Explosive Charges”, UK Ministry of Home Security, Civil Defense Research Committee Paper.
 Taylor, G.I. (1941a), “The Propagation and Decay of Blast Waves”, UK Home Office, ARP department.
 Taylor, G.I. (1941b), “The Propagation of Blast Waves over the Ground”, UK Ministry of Home Security, Civil Defence Research Committee paper.
 UFC 334002 (2008), “Structures to resist the effects of accidental explosions”, US Department of Defense, Washington (DC), USA
 Xia, Y., Wu, Ch., Li, Zh. X. (2015), “Optimized Design of Foam Cladding for Protection of Reinforced Concrete Members under Blast Loading”, ASCE Journal of Structural Engineering, ASCE J. Struct. Eng, 141(9), 17.
 Zieman, R.D., Mc Guire, W., Deierlein, G. (1992), “Inelastic limit states design, part I: planar frame studies”, ASCE Journal of Structural Engineering, 118 (9), 25322549.

Available Online: July 10, 2018
Publication Dates

Publication in this collection
2018
History

Received
04 Mar 2018 
Reviewed
18 June 2018 
Accepted
29 June 2018