Optimum Design of Braced Steel Space Frames including Soil-Structure Interaction via Teaching-Learning-Based Optimization and Harmony Search Algorithms

Optimum design of braced steel space frames including soil-structure interaction is studied by using harmony search (HS) and teaching-learning-based optimization (TLBO) algorithms. A three-parameter elastic foundation model is used to incorporate the soil-structure interaction effect. A 10-storey braced steel space frame example taken from literature is investigated according to four different bracing types for the cases with/without soil-structure interaction. X, V, Z, and eccentric V-shaped bracing types are considered in the study. Optimum solutions of examples are carried out by a computer program coded in MATLAB interacting with SAP2000-OAPI for two-way data exchange. -e stress constraints according to AISC-ASD (American Institute of Steel Construction-Allowable Stress Design), maximum lateral displacement constraints, interstorey drift constraints, and beam-tocolumn connection constraints are taken into consideration in the optimum design process. -e parameters of the foundation model are calculated depending on soil surface displacements by using an iterative approach. -e results obtained in the study show that bracing types and soil-structure interaction play very important roles in the optimum design of steel space frames. Finally, the techniques used in the optimum design seem to be quite suitable for practical applications.


Introduction
Optimum design of steel structures prevents excessive consumption of the steel material.Suitable cross sections must be selected automatically from a predefined list.Moreover, selected profiles should satisfy some required constraints such as stress, displacement, and geometric size.Metaheuristic search techniques are highly preferred for problems with discrete design variables.
ere are many metaheuristic techniques developed recently.Some of them are genetic algorithm, harmony search algorithm, tabu search algorithm, particle swarm optimization, ant colony algorithm, artificial bee colony algorithm, teaching-learningbased optimization, simulated annealing algorithm, batinspired algorithm, cuckoo search algorithm, and evolutionary structural optimization.In literature, there are many studies available for the optimum design of structures using these algorithms.For example, Daloglu and Armutcu [1] used the genetic algorithm method for the optimum design of plane steel frames.Kameshki and Saka [2] carried out the optimum design of nonlinear steel frames with semirigid connections using the genetic algorithm.Lee and Geem [3] developed a new structural optimization method based on the harmony search algorithm.Hayalioglu and Degertekin [4] applied genetic optimization on minimum cost design of steel frames with semirigid connections and column bases.Kelesoglu and Ülker [5] searched for multiobjective fuzzy optimization of space trusses by MS Excel.Degertekin [6] compared simulated annealing and genetic algorithms for the optimum design of nonlinear steel space frames.Esen and Ülker [7] optimized multistorey space steel frames considering the nonlinear material and geometrical properties.Saka [8] used the harmony search algorithm method to get the optimum design of steel sway frames in accordance with BS5950.Degertekin and Hayalioglu [9] applied the harmony search algorithm for minimum cost design of steel frames with semirigid connections and column bases.Hasancebi et al. [10] investigated nondeterministic search techniques in the optimum design of real-size steel frames.Hasançebi et al. [11] used the simulated annealing algorithm in structural optimization.Hasancebi et al. [12] investigated the optimum design of high-rise steel buildings using an evolutionary strategy integrated with parallel algorithm.Togan [13] used one of the latest stochastic methods, teaching-learning-based optimization, for design of planar steel frames.Aydogdu and Saka [14] used ant colony optimization for irregular steel frames including the elemental warping effect.Dede and Ayvaz [15] studied structural optimization problems using the teachinglearning-based optimization algorithm.Dede [16] applied teaching-learning-based optimization on the optimum design of grillage structures with respect to LRFD-AISC.Hasançebi et al. [17] used a bat-inspired algorithm for structural optimization.Saka and Geem [18] prepared an extensive review study on mathematical and metaheuristic applications in design optimization of steel frame structures.Hasançebi and Çarbas ¸ [19] studied the bat-inspired algorithm for discrete-size optimization of steel frames.Dede [20] focused on the application of the teaching-learningbased optimization algorithm for the discrete optimization of truss structures.Azad and Hasancebi [21] focused on discrete-size optimization of steel trusses under multiple displacement constraints and load case using the guided stochastic search technique.Artar and Daloglu [22] obtained the optimum design of composite steel frames with semirigid connections and column bases.Artar [23] used the harmony search algorithm for the optimum design of steel space frames under earthquake loading.Artar [24] used the teaching-learning-based optimization algorithm for the optimum design of braced steel frames.Carbas [25] studied design optimization of steel frames using an enhanced firefly algorithm.Daloglu et al. [26] investigated the optimum design of steel space frames including soil-structure interaction.Saka et al. [27] researched metaheuristics in structural optimization and discussions on the harmony search algorithm.Aydogdu [28] used a biogeography-based optimization algorithm with Levy flights for cost optimization of reinforced concrete cantilever retaining walls under seismic loading.
In literature, there are several researches available for optimum structural design, as mentioned above.On the other hand, there are a few researches on the optimum design of braced steel space frames including soil-structure interaction.So, this study investigates a 10-storey braced steel space frame structure studied previously in literature, which is investigated for four different bracing types and soil-structure interaction.
ese bracing types are X, V, Z, and eccentric V-shaped bracings.Optimum design solutions are obtained using a computer program developed in MATLAB [29] interacting with SAP2000-OAPI (open application programming interface) [30].Suitable cross sections are automatically selected from a list including 128 W profiles taken from AISC (American Institute of Steel Construction). e frame model is subjected to wind loads according to ASCE7-05 [31] as well as dead, live, and snow loads.e analysis results are found to be quite consistent with the literature results.In this study, the vertical displacements on soil surfaces are also calculated.It is observed that minimum weights of space frames vary depending on the bracing type.Also, it can be concluded that incorporation of soil-structure interaction results in heavier steel weight.

Optimum Design Formulation
e optimum design problem of braced steel space frames is calculated as follows: where W is the weight of the frame, A k is the cross-sectional area of group k, ρ i and L i are the density and length of member i, ng is the total number of groups, and nk is the total number of members in group k. e stress constraints according to AISC-ASD [32] are defined as follows: If (f a /F a ) ≤ 0.15, instead of using (2), the stress constraint is calculated as follows: where nc is the total number of members subjected to both axial compression and bending stresses, f a is the computed axial stress, F a is the allowable axial stress under axial compression force alone, f bx is the computed bending stress due to bending of the member about its major (x), F bx is the allowable compressive bending stress about major, F ex ′ is the Euler stress, F y is the yield stress of the steel, and C mx is a factor.It is calculated from C mx � 0.6 − 0.4(M1/M2) for the braced frame member without transverse loading between the ends and C mx � 1 + ψ(f a /F e ′ ) for the braced frame member with transverse loading.e effective length factors of columns in braced frames are calculated as follows [33]: where G A and G B are the relative stiffness factors at the Ath and Bth ends of columns.e maximum lateral displacement and interstorey drift constraints are defined as follows: 2 Advances in Civil Engineering where δ jl is the displacement of the jth degree of freedom under load case l, δ ju is the displacement at the upper bound, m is the number of restricted displacements, and nl is the total number of loading cases.
where Δ jil is the interstorey drift of the ith column in the jth storey under load case l, Δ ju is the limit value, ns is the number of storeys, and nsc is the number of columns in a storey.e beam-to-column connection geometric constraint is determined as follows: where n bw is the number of joints where beams are connected to the web of the column, b fbk,i ′ is the flange width of the beam, d c,i is the depth of the column, t fl,i is the flange thickness of the column, n bf is the number of joints where beams are connected to the flange of the column, and b fbk,i and b fck,i are flange widths of the beam and column, respectively (Figure 1).

Three-Parameter Vlasov Elastic Foundation Model
e soil reaction exerted on a structure resting on a twoparameter elastic soil is expressed in e reaction depends on the soil surface vertical displacement w, soil reaction modulus k, and soil shear parameter 2t.
ese two soil parameters, k and 2t, can be defined by in which H, υ s , and G s are the depth, Poisson's ratio, and shear modulus of the soil, respectively.In most of the classical two-parameter soil foundation models such as Pasternak, Hetenyi, and Vlasov models, the soil parameters are constants obtained by experimental tests or arbitrarily defined.However, it is highly difficult to determine these parameters experimentally.erefore, Vallabhan and Daloglu [34] developed an additional parameter c to characterize the vertical displacement profile within subsoil.
ey called this model including the third parameter c as a three-parameter Vlasov model. is model eliminates the necessity of experimental tests to determine soil parameters since these values are determined iteratively in terms of the new parameter, c. e vertical deformation profile of the subsoil is described via a mode shape function as given in e boundary values of ϕ(z) are assumed to be ϕ(0) � 1 and ϕ(H) � 0, as shown in Figure 2. e c parameter can be calculated using Equation (9) indicates that the soil parameters (k and 2t) are calculated based on the material properties and mode shape function (φ(z)).Also, it is necessary to compute the c parameter to calculate the mode shape function.It is necessary to know the soil vertical surface displacements obtained from the structural analysis to calculate the c parameter.So, it can be stated that k, 2t, ϕ, c, and w are interdependent.at is why the analysis requires an iterative procedure.For this purpose, a computer program is coded in MATLAB interacting with SAP2000 structural analysis program via OAPI (open application programming interface) to perform this iterative procedure in the threeparameter foundation model.

Advances in Civil Engineering
Using the coded program, a soil model is generated such that the soil reaction modulus k is represented by elastic area springs.e interaction between springs is taken into account using shell elements connecting the top of springs.e soil shell element with one degree of freedom at each node re ects only shear behavior of the soil.e c parameter is computed numerically in the coded program using the vertical displacements of soil shell elements.To determine the soil parameters iteratively, c 1 is assumed initially and k and 2t values are calculated.en, the structural model is analyzed using SAP2000, and the soil surface vertical displacements are retrieved to compute a new c value.e di erence between successive values of c is calculated and checked whether it is within a prescribed tolerance or not.If it is smaller than the tolerance, the iteration is terminated.Otherwise, the next iteration is performed, and the procedure is repeated until the convergence is ful lled.[3] and mimics improvisation procedures of musical harmony.It consists of three basic procedures.Operations are conducted by the harmony memory (HM) matrix.In the rst step, HM is randomly and automatically lled by the program coded in MATLAB.e form of harmony memory matrix is shown as follows:

Harmony Search Algorithm. Harmony search (HS) algorithm method is developed by Lee and Geem
where x j i is the ith design variable of the jth solution vector, n is the total number of design variables, φ(x j ) is the jth objective function value, and HMS (harmony memory size) indicates a speci ed number of solutions.In the harmony memory matrix, each row presents design variables.
In the second step, the objective function values )) of solution vectors in the harmony memory matrix are determined.In the third step, a new solution vector (x nh [x nh 1 , x nh 2 , . . ., x nh n ]) is prepared by selecting each design variable from either the harmony memory matrix or the entire section list X sl .Harmony memory consideration rate (HMCR) is applied as follows: Also, the new value of the design variable selected from the harmony memory matrix is checked whether this value should be pitch adjusted or not depending on the pitch adjustment ratio (PAR).
is decision is determined as follows: Yes, with probability of PAR No, with probability of 1 − PAR.
Detailed information about the HS algorithm can be found in the literature [3,8,9,23].
is method mimics teaching and learning processes between a teacher and students in classroom.
e person having the highest information in the class is selected as a teacher.e teacher gives his/her information to the other  ese procedures provide suitable solutions in structural optimizations.e teachinglearning-based optimization method consists of two basic steps such as teaching and learning.In the rst step, the teaching step, the rst population (class) is randomly lled in the matrix form presented as follows: class(population) where each row represents a student and gives a design solution, S is the population size (the number of students), n is the number of design variables, and f(x 1,2,...,S ) is the unconstrained objective function value of each student in the class.e student in a class having the best information is selected as a teacher of the class.His or her objective function value is the minimum in the class.e information update of students in the class is carried out with the help of the teacher as follows: where x new,i is the new student, x i is the current student, r is a random number in the range [0,1], and T F , a teaching factor, is either 1 or 2. x mean is the mean of the class de ned as x mean (mean(x 1 )....mean(x S )).If the new student has better information (f(x new,i )), the new student is replaced with the current student.In the second step, the learning step, information is shared between students.is step is similar to the rst step.If the new student presents a better information, he/she is replaced with the current student.e information update of students in the class is carried out as follows: e detailed information about the TLBO algorithm can be obtained from [13,15,16,20,24,35]. e owchart of processes in MATLAB-SAP2000 OAPI developed to get optimum solutions is presented in Figure 3.

Design Example
A 10-storey braced steel space frame example taken from literature [36] is studied considering four di erent types of Advances in Civil Engineering bracing such as X, V, Z, and eccentric V. e behavior of the frame is investigated with and without considering the e ect of soil-structure interaction.e frame example is exposed to wind loads according to ASCE7-05 [31] in addition to dead, live, and snow loads.Optimum cross sections are practically selected from a prede ned list of 128 W pro les taken from AISC. e stress constraints according to AISC-ASD [32], maximum lateral displacement constraint (H/400), interstorey drift constraint (h/400), and beam-tocolumn geometric constraints are subjected to the optimum design of the braced steel space frames.In the analyses, the steel modulus of elasticity, E, and yield stress, F y , are taken as     Figure 4 shows the typical plane view of a 10storey steel frame.Also, Figure 5 represents the three dimensional view of a V-braced frame.Each storey has a height of 3.66 m (12 ft).Modulus of elasticity for the concrete is taken as 32,000,000 kN/m 2 , Poisson's ratio is 0.2, and weight per unit volume is 25 kN/m 3 .All floors excluding the roof are exposed to a dead load of 2.88 kN/m 2 and a live load of 2.39 kN/m 2 .e roof floor is exposed to a dead load of 2.88 kN/m 2 and a snow load of 0.75 kN/m 2 .e total gravity loading on the beams of roof and floors is tabulated in Table 1 [36].Wind loads are applied to the frame according to ASCE7-05 [31].Wind loads calculated for the 10-storey braced frame are presented in Table 2 [36].Modulus of elasticity of the soil, E s , is taken to be equal to 80,000 kN/m 2 .e depth of the soil stratum to the rigid base is taken as H s � 20 m, and Poisson's ratio of the soil is equal to 0.25.

X-Braced Steel Space
Frame. Figure 6 shows the typical 2D view of an X-braced steel space frame without and with soil-structure interaction.Optimum results of the X-braced space frame are given in Table 3. Soil parameters for the  8 Advances in Civil Engineering space frame on elastic foundation are presented in Table 4.Moreover, Figures 7 and 8 display settlements of the soil surface for two different algorithm methods and design histories of optimum solutions, respectively.

V-Braced Steel Space
Frame. Figure 9 shows the typical 2D view of a V-braced steel space frame without and with soil-structure interaction.Optimum results of the V-braced space frame are shown in Table 5. Soil parameters for the space frame on elastic foundation are given in Table 6.Also, Figures 10 and 11 show the soil surface settlements obtained using two different algorithms and design histories of optimum solutions, respectively.

Z-Braced Steel Space
Frame. Figure 12 represents the 2D side view of a Z-braced steel space frame without and with Advances in Civil Engineering soil-structure interaction.Optimum results of the Z-braced space frame are presented in Table 7. Soil parameters for the space frame on elastic foundation are given in Table 8.Moreover, Figures 13 and 14 represent settlements of the soil surface carried out with two different algorithm methods and design histories of optimum solutions, respectively.

Eccentric V-Braced Steel Space
Frame. Figure 15 shows the 2D side view of an eccentric V-braced steel space frame without and with soil-structure interaction.e braces are connected to the beam from one-third of the beam length.Optimum results of the eccentric V-braced space frame are presented in Table 9. Soil parameters for the space frame on elastic foundation are presented in Table 10.Moreover,   It is observed from Tables 3, 5, and 7 that the minimum weights of the braced frames for the case without soilstructure interaction are very similar to the ones available in literature research [36].In this study, the V-braced type provides the lowest steel weight of 1083.80 kN by using teaching-learning-based optimization.Z-braced and X-braced types provide the second and third low weights, 1095.27kN and 1170.03,respectively.Moreover, the minimum weight of the eccentric V-braced frame, 1275.01 kN, is nearly 15%, 14.1%, and 8.2% heavier than the minimum steel weights of the V-, Z-, and X-braced frames, respectively.On the other hand, harmony search algorithm presents 2.6-4.7%heavier minimum steel weights than the ones obtained from teaching-learning-based optimization for the case without soil-structure interaction.e tables including optimum results also show that the value of interstorey drift is very close to the limit value (h/400).erefore, the displacement constraints play very crucial roles in the optimum design of the braced frames.Five Advances in Civil Engineering independent runs are performed for each braced type for the case without soil-structure interaction.
In the case with soil-structure interaction, the minimum weights of all braced frames increased depending on settlements on the soil surfaces.It is observed from Tables 4, 6, 8, and 10 that the soil parameters of 10-storey braced steel frames on elastic foundation are similar for all cases.e minimum steel weights are mostly obtained by teachinglearning-based optimization.For the X-braced frame, the minimum weight obtained by TLBO for the case with soilstructure interaction is 4.01% heavier than the weight of the frame excluding soil-structure interaction.
is ratio is 5.32% for the harmony search algorithm.Moreover, settlement values on the soil surfaces are nearly −0.66 cm as seen in Figure 7.For the V-braced frame including soilstructure interaction, TLBO and HS present 1.66 and 5.9% heavier weights, respectively.e settlements in this braced  12 Advances in Civil Engineering frame are nearly −0.65 cm as given in Figure 10.For the Z-braced frame with soil-structure interaction, the minimum weights obtained are 5.11% and 3.73% heavier by using TLBO and HS, respectively.e settlements in this braced frame are similar to the ones of the other braced frames.For the eccentric V-braced frame with soil-structure interaction, the minimum steel weights obtained are 6.27 and 7.94% heavier by using TLBO and HS, respectively.Moreover, the convergences of optimum solutions with iteration steps are seen in Figures 8,11,13, and 17 in detail.

Summary and Conclusions
In this study, the optimum design of a 10-storey steel space frame braced with X, V, Z, and eccentric V-shaped bracings including soil-structure interaction is investigated.Optimum solutions are obtained using two different metaheuristic algorithm methods: teaching-learning-based optimization (TLBO) and harmony search (HS).For this purpose, a code is developed in MATLAB computer program incorporated with SAP2000-OAPI (open application programming interface).Required cross sections are automatically selected from a list of 128 W profiles taken from AISC (American Institute of Steel Construction). e Advances in Civil Engineering 13 frame model is exposed to wind loads according to ASCE7-05 in addition to dead, live, and snow loads.e stress constraints in accordance with AISC-ASD (American Institute of Steel Construction-Allowable Stress Design), maximum lateral displacement constraints, interstorey drift constraints, and beam-to-column connection constraints are applied in analyses.A three-parameter Vlasov elastic foundation model is used to consider the soil-structure interaction effect.e summary of the results obtained in this study are briefly listed below: (i) It is observed from analyses that the minimum weight of the space frame varies by the types of bracing.e lowest steel weight, 1083.80 kN, is obtained for the V-braced steel frame by using TLBO.Z-braced and X-braced types provide the second and third low weights, 1095.27kN and 1170.03,respectively.ese results are similar to the ones available in literature [36].e heaviest among them is the minimum weight of the eccentric V-braced frame, 1275.01 kN.(ii) Harmony search algorithm presents 2.6-4.7%heavier steel weights than the ones obtained from teaching-learning-based optimization for the frames without soil-structure interaction.Although the lighter analysis results are obtained in TLBO, a representative structure model in TLBO is analyzed twice in an iteration step by SAP2000 programming.On the other hand, it is enough to analyze the system once in HS. is situation requires longer time for the analysis in TLBO.(iii) Interstorey drift values are very close to its limit value of 0.915 cm (h/400).erefore, the constraints are important determinants of the optimum design of the braced frames.(iv) Consideration of soil-structure interaction results in heavier steel weight.For the X-braced frame including soil-structure interaction, the minimum weights are obtained to be 4.01 and 5.32% heavier by using TLBO and HS, respectively.For the V-braced frame, these values are calculated to be 1.66 and 5.9% heavier, respectively.For the Z-braced frame, these values are obtained to be 5.11% and 3.73% heavier, respectively.Moreover, for the eccentric V-braced frame, the minimum weights are obtained to be 6.27 and 7.94% heavier by these algorithm methods.(v) Settlement values on the soil surfaces are nearly 0.61-0.67cm for all braced frames.(vi) Finally, the techniques used in optimizations seem to be quite suitable for practical applications.An adaptive setting for the parameters will be very useful and user-friendly especially for the structures with a large number of members as in the case here.is will be considered in future studies.Advances in Civil Engineering

Figure 2 :
Figure 2: A space frame on three-parameter elastic foundation.(a) 3D frame on elastic subsoil.(b) Mathematical model.

Figure 3 :
Figure 3: Flowchart for the optimum design algorithm by HS and TLBO for space frames on elastic foundation.

Figure 4 :
Figure 4: Typical plane view of a 10-storey steel frame.

Figure 5 :
Figure 5: ree dimensional view of a V-braced frame.

Figure 6 :
Figure 6: 2D view of the X-braced steel space frame without and with soil-structure interaction.(a) e case without soil-structure interaction.(b) e case with soil-structure interaction.

Figure 7 :Figure 8 :Figure 9 :
Figure 7: Settlements of soil surface for the X-braced steel space frame with two different algorithm methods (cm).(a) TLBO.(b) HS.

-Figure 10 :Figure 11 :
Figure 10: Settlements of soil surface for the V-braced steel space frame with two different algorithm methods (cm).(a) TLBO.(b) HS.

Figure 12 :
Figure 12: 2D view of the Z-braced steel space frame without and with soil-structure interaction.(a) e case without soil-structure interaction.(b) e case with soil-structure interaction.

Figures 16 and 17
Figures 16 and 17  show settlements of the soil surface carried out with two different algorithm methods and design histories of optimum solutions, respectively.It is observed from Tables 3, 5, and 7 that the minimum weights of the braced frames for the case without soilstructure interaction are very similar to the ones available in literature research[36].In this study, the V-braced

-Figure 13 :Figure 14 :
Figure 13: Settlements of soil surface for the Z-braced steel space frame with two different algorithm methods (cm).(a) TLBO.(b) HS.

Figure 15 :
Figure 15: 2D view of the eccentric V-braced steel space frame without and with soil-structure interaction.(a) e case without soil-structure interaction.(b) e case with soil-structure interaction.

Figure 16 :Figure 17 :
Figure 16: Settlements of soil surface for the eccentric V-braced steel space frame with two different algorithm methods (cm).(a) TLBO.(b) HS.

Table 4 :
Soil parameters for 10-storey X-braced steel space frame on elastic foundation.

Table 3 :
Optimum results of the X-braced space frame.

Table 5 :
Optimum results of the V-braced space frame.

Table 6 :
Soil parameters for the 10-storey V-braced steel space frame on elastic foundation.

Table 7 :
Optimum results of the Z-braced space frame.

Table 8 :
Soil parameters for the 10-storey Z-braced steel space frame on elastic foundation.

Table 9 :
Optimum results of the eccentric V-braced space frame.

Table 10 :
Soil parameters for the 10-storey eccentric V-braced steel space frame on elastic foundation.