Microsoft Word - numero_59_art_22_3280.docx J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 326 Optimization of a reinforced concrete structure subjected to dynamic wind action Jherbyson Williams Silva Brito PPGEC/UFRGS, Brazil Jherbyson.willians@gmail.com, https://orcid.org/0000-0002-1481-1803 Letícia Fleck Fadel Miguel Department of Mechanical Engineering, Federal University of Rio Grande do Sul, Brazil letffm@ufrgs.br, https://orcid.org/0000-0001-9165-4306 ABSTRACT. This work proposes a methodology to optimize a reinforced concrete structure. For this, the Whale Optimization Algorithm (WOA) algorithm was used, an algorithm from the group of metaheuristic algorithms, which presents an easy computational implementation. As a study object, a frame structure adapted from a real reinforced concrete building was used, subjected to the dynamic action of artificially generated synoptic wind. The objective function is to reduce the volume of concrete of the structure. For that, the dimensions of the cross-sections were used as design variables, and the maximum displacement at the top imposed by the ASCE / SEI 7-10 standard as a lateral constraint, as well as the maximum story drift between floors. In addition to this structural optimization, it was also proposed the use and optimization of Tuned Mass Dampers (TMD), in different quantities, positions and parameters, improving the dynamic response of the reinforced concrete building. The results show that for this situation it was possible to reduce the concrete volume of the structure by approximately 24%, respecting the maximum limit of displacement at the top required by the standard. KEYWORDS. Optimization; Concrete structure; Metaheuristic algorithms Citation: J. W. S. Brito and L. F. F. Miguel., Optimization of a reinforced concrete structure subjected to dynamic wind action, Frattura ed Integrità Strutturale, 59 (2022) 326-343. Received: 27.09.2021 Accepted: 09.11.2021 Published: 01.01.2022 Copyright: © 2022 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. INTRODUCTION t is notable that buildings are getting thinner, mainly due to technological innovations in relation to constructive solutions and quality of new materials. Previously low and robust constructions now make room for increasingly taller and more flexible structures, because of reduction in cross-sections due to the gain in strength of the new materials. In relation to the structural design of tall buildings, the dynamic effects due to the wind are important, because if the proper care is not taken in the design, excessive vibrations may arise, causing discomfort to residents and/or users and consequent economic losses, and in an extreme case, the structure collapses, as has already happened in several cases around the world. Due to constant development and the search for efficiency and competitiveness in cost and weight reduction in structural projects, optimization methods are increasingly being applied in search of the best project, in which applications in the I https://youtu.be/Re9zARD1d3s J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 327 aerospace areas, engineering projects are being increasingly applied, equipment, automotive, among other areas of engineering. One of the first articles on optimization in tall buildings was that by [1], in which the authors describe the evolution of a group of optimization algorithms created by them, with the objective of optimizing the areas of the structural model bars, submitted a design maximum lateral constraint. Venanzi and Materazzi [2] performed a multi-objective optimization process of a communication tower for mobile devices subjected to turbulent wind loading. Yang et al [3] analyzed the control of the dynamic response of a tall building using a tuned mass damper on the top of the structure. It is a Benchmark structure, an example used as a study reference by several authors. Many classical or conventional algorithms for structural optimization are deterministic and most of them use gradient information, that is, they use function values and their derivatives. They usually work extremely well for smooth unimodal functions, however, if there is any discontinuity in the objective function, they may not converge [4]. Due to these difficulties, Heuristic and Metaheuristic algorithms have been increasingly used, which are based on randomization and local search from candidate populations, and not from a single point as happens with classical methods [5]. These algorithms are classified into four categories: evolution-based, physics-based, swarm-based, and human-behavior- based. Some examples of algorithms are the Genetic Algorithm (GA) [6], which is an evolution-based algorithm, the physics- based Gravitational Search Algorithm (GSA) [7] and the Harmony Search (HS) [8] based on human behavior. Swarm-based optimization methods reproduce the social behavior of groups of living beings. Among the best-known algorithms for this method are Particle Swarm Optimization (PSO) [9] and Ant Colony Optimization [10]. In this paper will be used the Whale Optimization Algorithm (WOA) [11] developed by Mirjalili and Lewis. Optimization algorithms can be used to optimize parameters, positions and quantities of passive energy dissipation devices, with the objective of reducing the dynamic response of the structure, respecting normative limit values. This is used to reduce the amplitude of oscillations caused by dynamic excitation in different types of structures. These systems aim to reduce the dynamic stresses that can occur in buildings, avoiding problems such as fatigue and ensuring structural safety. The purpose of structural control in civil engineering structures is to reduce vibrations produced by external forces such as winds and earthquakes, by different techniques such as modification of stiffness, mass, damping and geometric shape [12]. Vibration control systems can be divided into two main groups, active and passive systems. Passive control systems are devices composed of a mass connected to a spring and a damper to the structure, so that this device absorbs part of the energy from dynamic loading, dissipating the energy in the members of the structure. This system has advantages because it does not have the need for energy or control external to the structure, working by the vibration of the building itself, thus being of low cost compared to the active system, as well as installation, maintenance and ease of control. It’s optimally tuned to protect the structure from dynamic load at a specified frequency. However, it’s not efficient for other cases and other types of dynamic loads. This is one of the negative points in relation to active systems [13]. Active systems adjust to different vibration frequencies through sensors requiring information on structural behavior and external energy. This system is operated by hydraulic or electromechanical actuators that provide control forces to the structure from monitoring with sensors that measure excitation and/or response due to dynamic load [14]. A classic device widely studied by structural engineering researchers is the Tuned Mass Damper (TMD). It’s an efficient mechanical device that features low cost and low maintenance. Furthermore, these devices can be installed, for example in a building, without interrupting the building's operational activities. Another advantage these attenuators have over other control devices is their versatility, as they can be designed in many different shapes and sizes, which facilitates adaptations to architectural aspects and space limitations [15]. Several researchers use this device, and prove the reduction of displacement peaks at the top of structures ([3],[16] and [17]). So, the main objective of the present work is to optimize the volume of a reinforced concrete structure according to maximum displacement restrictions established by international reference standards. For this, the optimization algorithm proposed by [11] is used. The cross-sections of the members are considered as design variables. After the optimization of the structure, TMD’s are installed and optimized, in different quantities, positions and parameters. APPLICATION OF TMDS IN THE STRUCTURE here are several TMD’s design methods in the literature, among which the [18] and [19] can be highlighted, which are used for the design of simple TMD installed on top of the buildings studied. From these methods, the optimal values of stiffness, damping and tuning frequencies of the device can be found. However, in this paper these values T J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 328 will be obtained from the optimization process performed by the algorithm used, both for simulations with a single TMD at the top and for TMDs. To analyze the response of the structure with a TMD on top, as shown in Fig. 1, it is necessary to make changes in the stiffness, mass and damping matrices, in order to take into account the Degree of Freedom (DOF) of TMD. In general, the matrices can be modified according to Eqns. 1, 2 and 3, taking into account that the building is of the shear building type and only the horizontal DOF is considered, therefore, the global mass matrix is of the type diagonal, where the mass of each floor and each TMD occupies a position on the main diagonal.                       1 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 n TMD m m m m M (1)                        1 2 2 2 2 3 0 0 0 0 0 0 0 0 0 0 n TMD TMD TMD TMD k k k k k k k k k k k K (2) Figure 1: n-degree-of-freedom system (n-DOF) structure equipped with one TMD at the top floor J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 329                        1 2 2 2 2 3 0 0 0 0 0 0 0 0 0 0 n TMD TMD TMD TMD c c c c c c c c c c c C (3) In many cases, the use of a single TMD can be enough to solve the problems of excessive vibration of structures. However, for some situations, the use of only one TMD may not be efficient, as in the case of structures subjected to wind and earthquake, when more than one mode is excited. Another difficulty would be with the mass of the TMD being too large, requiring a very large space for installation. In this sense, there is the possibility of using multiple TMD’s (MTMD), each tuned to a specific frequency of the structure, thus being more efficient than a single TMD. While a single TMD is usually installed on top, MTMD can be installed in different positions and configurations, in series, in parallel, on a single floor or multiple floors. In relation to the MTMD, they can be arranged horizontally, vertically, one per floor or several per floor, as shown in Fig. 2. Thus, for TMD’s arranged horizontally and vertically in a building and in cases where the mass matrix is considered diagonal, the equation that indicates the modification of the matrix M is similar to the case of a single TMD, as already presented in Eqn. 1, and for the stiffness matrix, which is analogous to the damping matrix, we have Eqn. 4. Some examples of the use of TMD in buildings are Citicorp Center, opened in 1977, in New York, USA. This was the first application of a TMD for wind response control, with the metallic structure being 279 meters high. The attenuator was built in concrete, weighing 400 tons and installed on top of the building. The TMD mass is calculated from the mass ratio between the TMD and the structure. These values, according to researchers, can vary between 0.01 and 0.2, that is, up to about 20% of the structure's mass. It is known that the greater the mass ratio between the TMD and the structure, the greater the TMD mass value. However, the value of the mass of the TMD is not interesting because it is a very high value, as it will add additional load to the structure under analysis, making it unfeasible due to the cost and difficulty of implementation. Therefore, the optimization of the attenuator parameters becomes of fundamental importance, as well as their quantity and position. However, this problem is not so simple, requiring the use of optimization algorithms that are widespread in engineering problems. In this paper, the stiffness and damping parameters of each attenuator are obtained through structural optimization. Figure 2: n-degree-of-freedom system (n-DOF) structure equipped with TMDs possibly located in all floors of the structure J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 330                                                   1 2 1 2 ,1 2 1,1 2,1 ,1 1 2 2 3 ,2 1,2 2,2 ,2 1 , , 1 1,1 1,1 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 stFloor TMD stFloor TMD stFloorn TMD m N di d d dm i m N di d d dm i m N n di n dm n i d d d k k k k k k k k k k k k k k k k k k k k K                                                                                            ,1 2,1 ,1 ,1 1,2 1,2 2,2 2,2 ,2 ,2 , , 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 d dm dm d d d d dm dm dm n dm n k k k k k k k k k k k         (4) PROPOSED METHODOLOGY Analyzed Structure he structural model is a 42-story reinforced concrete buildong, ν = 0.2, p = 2500 kg/m³, and modulus of elasticity calculated in accordance with [20]. The total height of the structure is 105.38 meters, with 324 bars, modeled in Octave using the 2D frame element with 2 nodes and 3 degrees of freedom (DOFs) for each node. There are 185 nodes and therefore 555 DOFs. The sections of the bars are rectangular, with dimensions as explained in Tab. 1. Bars Dimensions (cm) Areas (m²) Inertia Moment (m4) 1 to 36 / 145 to 180 100x25 0.25 1.30e-4 37 to 72 / 109 to 144 25x120 0.30 0.036 73 to 108 40x157 0.628 0.129 181 to 324 12x65 0.078 2.75e-3 Table 1: Geometrical properties of the structure The structure is shown in Fig. 3. The dimensions are in meters (m). The Rayleigh Damping Matrix is used, where the critical damping ratio (ζ) was specified as 1% for the first two modes of vibration. The mass matrix of the structure is consistent where for each element, the mass and stiffness matrix in the local system, is represented by ML and KL                     2 2 140 0 0 70 0 0 0 156 22 0 54 13 0 22 4 0 13 3 70 0 0 140 0 0420 0 54 13 0 156 22 0 13 3 0 22 4 L L L L L LAL L L L L L L Lm (5) T J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 331                                  3 2 3 2 2 2 2 3 2 3 2 2 2 2 0 0 0 0 12 6 12 6 0 0 6 4 6 2 0 0 0 0 0 0 12 6 12 6 0 0 6 2 6 4 0 0 L EA EA L L EI EI EI EI L L L L EI EI EI EI LL L L EA EA L L EI EI EI EI L L L L EI EI EI EI LL L L k (6) A code was developed to calculate the natural frequencies of the structure, and later it was compared with the commercial software Ansys Mechanical APDL. The results of the first 3 natural frequencies are shown in Tab. 2: Vibration Mode Octave Frequency (Hz) Ansys Frequency (Hz) Error (%) 1 0.3989 0.3989 0 2 1.2829 1.2829 0 3 2.3999 2.3994 0.03 Table 2: Comparison between natural frequencies Wind Excitation For the generation of the synoptic wind, the wind action is considered as a stochastic process, that is, defined through an unlimited sequence of random variables. The fluctuating component of the wind is given through atmospheric turbulence, considering only the effect of the kinetic energy of the gusts, disregarding the neighborhood effects. According to the [21], the component of the global wind force, called drag force, is calculated according to Eqn. 7:  20.613 a zF C A V (7) where Ca is the drag coefficient of the structure, Az is the area of application of the dynamic wind pressure, in m² and V is the wind speed component. The drag coefficient is a function of the building geometry and the wind turbulence level, using the coefficients in Fig. 4 from [21] for low turbulence winds and Fig. 5 for high turbulence winds. The longitudinal direction is the only one that contains significant mean velocities, being considered the predominant direction. The characteristic wind velocity is given by the sum of the static and dynamic plots, can be written as Eqn. 8:     ( )x xV t V v t (8) where V is the component of the mean velocity in the longitudinal direction of the wind, ( )xv t is the floating portion of the velocity in longitudinal direction. Initially, the design speed is defined through Eqn. 9, which corresponds to the average speed in a 10-minute interval at a height of 10 meters from the ground:  0 1 2 30.69PV V S S S (9) where 0V is the basic reference wind speed where the structure is located, of a 3 s gust, exceeded on average once in 50 years, at 10 m above ground, in flat open ground. 0V must be adequate to the location of the construction, and can be obtained by the isopleths of the basic wind speed of the [21]. J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 332 Figure 3: structure studied in this paper J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 333 The factor 1S is the topographic factor that takes into account variations in terrain relief. It has a value of 1 for slightly rough terrain. 3S is a statistical factor based on the type of use of the building, with a value of 1 for the case studied, residential building. The factor 2S considers the combined effect of the terrain roughness, the variation in wind speed with the height above the terrain and the dimensions of the building or part of the building under consideration, being calculated according to Eqn. 10, with the help of Tab. 1 of [21]: 2 ( 10) p rS bF z (10) The floating component of the wind can be calculated according to the procedures of chapter 9 of [21] (Dynamic effects due to atmospheric turbulence). However, there are other procedures that, despite being more complex, are more efficient and better suited to different structures. One of these procedures is the analysis using the spectral representation method. Considering the fluctuating wind speed as a normal random process with a mean of zero, it is possible to reach the fluctuating component of the wind from the superposition of harmonic waves [22] by Eqn. 11:         1 ( ) 2 ( ) cos(2 ) N W j j j j j V t S f f f t (11) where ϕj is a random variable with a uniform probability distribution defined between 0 and 2π.  jf is obtained by  1j jf f , that is, from the division of the frequency range of interest. The calculation of the power spectral density is given by the model proposed by Harris (Eqn. 12):  2 2 5/6 * ( ) 4 (2 ) wfS f n u n (12) where:  0 WSf L n V (13)         * 0 0.4 ln refz ref V u z z (14) where V0 is the reference speed at a height of 10 meters where the structure is located, fsw is the frequency band used in the spectrum, u* is the wind flow cutoff speed and z0 is the roughness length, which is related to the height of the obstacles making up the roughness of the surface, it can be roughly estimated as one-twentieth of the average height of the obstacles. The roughness length is a parameter sensitive to changes in dimensions and density of obstacles; and, therefore, it is recommended to pay attention to the values obtained, and the use of small values favors security [23]. After calculating the velocities, the next step is calculating the correlation length. The length of vertical correlation between two points will be calculated according to the expression given by [24]. The studied structure will be inserted in the correlation plane, and the velocity for each node of the structure will be obtained through linear interpolation. The results are shown in Figs. 4 and 5. From the length of correlation, the fluctuating velocities at all nodes of the structure were obtained and then the vector of applied forces in time for each node of the structure was obtained. Once the force vector is calculated, the dynamic response itself is started (displacements, velocities and accelerations) through the Newmark Integration Method, using a time step dt = 0.1s and 999 integration points. From the analysis, an average maximum displacement at the top of approximately 0.38m was obtained, for situations where there is no energy dissipator. This value is considered high for design situations. J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 334 The Whale Optimization Algorithm (WOA) Optimization algorithms have been undergoing a constant evolution process, and this process includes metaheuristic algorithms. The Whale Optimization Algorithm (WOA) is a metaheuristic algorithm that is based on the hunting behavior of humpback whales. In general, humpback whales hunt prey through the formation of bubbles while directing them through spiral movements. As with other evolutionary algorithms, WOA starts its search process through the random generation of an initial population of candidates (whales) delimited by the design space defined by the user of the technique. The new population is formed from the previous population and new candidates generated from three operators. Once the population is created, the optimal solution is given as the best candidate of this population, generating a new population around the optimal value. This procedure is repeated until a defined stopping criterion, usually the number of repetitions of this procedure, which is Figure 4: Harris power spectrum Figure 5: Floating longitudinal component in EPS wind J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 335 defined by the programmer before the simulations. The pseudo-code of WOA is shown in Fig. 6 and according to [24] the following input parameters are necessary: Dim (number of design variables); fobj (objective function); Nsa (number of search agents, that is, the whale population); Ngen (maximum number of generations, that is, maximum iteration number); Lb (lower bound, where Lbn the lower bound of variable n, for example: Lb = [Lb1, Lb2,..., Lbn]); Ub (upper bound, where Ubn the upper bound of variable n, for example: Ub = [Ub1, Ub2,..., Ubn]). Fig. 7 shows the parameters which must be updated, where a is decreased from 2 to 0 in order to provide exploration and exploitation, respectively; A and C, are coefficients utilized to calculate the best current solution; l, is a random number in [-1,1]; and p, is a random number in [0,1] [25]. Formulation of the dynamic ptimization problem Faced with high displacement values, it is essential to resize the structure so that it respects the maximum displacement values according to appropriate codes and standards, whether national or international. In this work, we chose to use as a basis the American standard ASCE/SEI 7-10 [26], which according to its appendix C, comments that the maximum displacement value of a structure is calculated as dMáx = H/400, where H is the total height of the building. Therefore, for this problem, the maximum displacement constraint is approximately 0.26m. In addition to the maximum displacement at the top of the structure, there are also restrictions on relative displacements between floors (story drift). According to the American standard, the story drift of each floor cannot exceed 10mm (approximately 3/8 of an inch). The design variables are the heights of the cross-sections of the beams and columns (Fig. 7), keeping the base values fixed. In addition to member heights, the damping and stiffness values of each TMD are also inserted as design variables in the optimization problem. The lateral limits are in Tab. 3: Design variables Lower bound Upper bound x1/x3/x5/x7/x9 (cm) 80 120 x2/x4/x6/x8/x10 (cm) 20 80 x11 to x46 (cm) 12 100 TMD’s stiffness (N/m) 0 1.955e4 TMD’s damping (Ns/m) 0 391.08 Table 3: Lateral limits of design variables Figure 6: Pseudo-code of the WOA [25]. Initialize the whales population Xi (i = 1, 2, ..., n) Calculate the fitness of each search agent X*=the best search agent while (j < maximum number of iterations) for each search agent Update a, A, C, l, and p if1 (p < 0.5) if2 (|A| < 1) Update the position of the current search agent else if2 (|A| ≥ 1) Select a random search agent (Xrand) Update the position of the current search agent end if2 else if1 (p ≥ 0.5) Update the position of the current search end if1 end for Check if any search agent goes beyond the search space and amend it Calculate the fitness of each search agent Update X* if there is a better solution j = j+1 end while return X* J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 336 Figure 7: Design variables of the optimization problem. J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 337 The objective function consists in the minimization of the volume of the concrete of the structure, calculated as:   1 barsn structure i i i V L A (15) where L and A are the lengths and area of each member, respectively. Therefore, the optimization problem is formulated as follows: Find:  i iX h where  1, 2, , 46i Find:   ,i i iT k c where  1, 2, , AMSi N Minimize    1 barsn structure i i i V L A Subject to  0.26Maxd m e    1 10h i iv v mm Once the optimization problem was defined, 4 scenarios to control the structure were proposed, seeking the best performance of the structure with the smallest possible volume of concrete: Scenario 1: Optimization of the structure, without the presence of TMD; Scenario 2: Optimization with a TMD on top of the building; Scenario 3: Optimization with two TMDs on top of the building; and Scenario 4: Optimization with two TMD’s, one on top of the building and one on the next floor (penultimate floor). Model Concrete Volume (m³) Average Volume (m³) Original 260.72 260.72 Otim. 1 226.14 223.70 Otim. 2 219.22 Otim. 3 232.26 Otim. 4 230.61 Otim. 5 229.10 Otim. 6 217.98 Otim. 7 225.39 Otim. 8 215.72 Otim. 9 213.49 Otim. 10 223.75 Table 4: Comparison between concrete volume values from scenario 1 simulations. RESULTS AND DISCUSSIONS Scenario 1 he first scenario consists of optimizing the structure from the height of the cross-sections of the bars, without the presence of TMD, under dynamic excitation due to the wind. For this, the WOA optimization algorithm is used. 30 research agents and 900 iterations were used. Due to the probabilistic nature of the optimization algorithm, 10 independent simulations were performed and the results averaged. The number of iterations needed was evaluated using the iteration curve versus concrete volume (Fig. 8), noting that from iteration 4 or 5, the convergence of the value already occurs. However, the 30 iterations are used to ensure that all simulations are at the optimal value. 10 simulations were performed and the average of the results was calculated, where values are listed in Tab. 4. It is noticed that from the optimization process it was possible to reduce the volume of concrete in the original structure by approximately 14.33%. To verify whether the optimization algorithm was respecting the design maximum displacement T J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 338 constraints, the heights of the bars given by the algorithm in one of the simulations were used, and the dynamic analysis was performed again, comparing the optimized top displacement with the displacement in the top without optimization. The results are shown in Fig. 9. It is important to emphasize that in addition to having managed to bring the building to a maximum displacement value within the acceptable limit in the standard (it reduced from 0.38m to 0.24m), the concrete volume of the structure was also reduced. Figure 8: Convergence curve of Scenario 1. Figure 9: Displacement on top of optimized structure Scenario 2 In this scenario, a TMD was inserted at the top of the structure, connected to node 183 of the structural model. This scenario presents the classic vibration control solution with TMD, by adding one more DOF to the structure. The definition of the location of the energy dissipator is due to the fact that the building's first vibration mode generates greater movements on higher floors. Different from scenario 1, scenario 2 has two more design variables, which are TMD stiffness and damping. All other considerations are identical to scenario 1. The TMD’s mass was adopted as 3% of the structure mass as J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 339 recommended in the bibliography, resulting in about 20 tons. The results are shown in Tab. 5. As in the previous case, the number of iterations was evaluated through the convergence of the optimal response to the problems, and it is clear that from the 12th iteration on, the result has already been stabilized. It is noticed that the simulation with TMD reached values approximately 11% smaller in concrete volume than the simulation of scenario 1, and approximately 23.59% smaller than the original volume of the structure, respecting the maximum values of displacement at the top and displacement relative between floors (story drift). This behavior is based on the existing literature, showing that the insertion of TMD at the top of buildings, with its appropriate optimized parameters, reduces the amplitude of displacements due to the dynamic excitation of the wind, and thus, the volume of structure without exceeding the normative displacement limit. The average TMD properties calculated by the algorithm are mentioned in Tab. 6. Model Concrete Volume (m³) Average Volume (m³) Original 260.72 260.72 Otim. 1 187.88 199.38 Otim. 2 187.36 Otim. 3 201.16 Otim. 4 205.38 Otim. 5 212.44 Otim. 6 199.11 Otim. 7 201.5 Otim. 8 206.42 Otim. 9 203.4 Otim. 10 189.1 Table 5: Comparison between concrete volume values from scenario 2 simulations. Parameters Average Value Stiffness (N/m) 8315.12 Damping (Ns/m) 97.93 Table 6: TMD parameters obtained by optimization Scenario 3 In scenario 3, two TMDs were inserted into the structure, equivalent to scenario 2, on the top floor of the building, at nodes 182 and 184. The optimal dissipator’s design is performed for the mass ratio defined in scenario 2, that is, a total mass of 3% (1.5% for each of the 2 TMD). That is, the resulting value of approximately 20 tons is divided equally between each TMD considered in each current analysis. WOA parameters remained the same as in the previous scenario for all analyzes performed in this third scenario. As in the previous case, the number of iterations was evaluated through the convergence of the optimal answer to the problems, and it is possible to verify that convergence was reached in the 13th iteration, but it is important to perform all iterations. The concrete volume values obtained in the simulations are shown in Tab. 7. It can be seen from the results in the table that the average value of concrete volume of the optimization process is slightly higher than the average value of the previous scenario, with a single TMD at the top (approximately 4%). However, this J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 340 optimization allowed a reduction in the volume of concrete of 20.48% compared to the original structure, a good reduction in construction. The stiffness and damping values of the TMDs are in Tab. 8. There is a considerable difference between the characteristics of the two TMDs. Installing devices with different properties can be difficult or more costly compared to a solution that has the same specifications for each TMD, so for this simulation, the scenario with only one TMD on top becomes more efficient. Model Concrete Volume (m³) Average Volume (m³) Original 260.72 260.72 Otim. 1 187.88 207.32 Otim. 2 187.36 Otim. 3 201.16 Otim. 4 205.38 Otim. 5 212.44 Otim. 6 199.11 Otim. 7 201.5 Otim. 8 206.42 Otim. 9 203.4 Otim. 10 189.1 Table 7: Comparison between concrete volume values from scenario 3 simulations. Parameters Average Value Stiffness of TMD 1 (N/m) 13798.07 Damping of TMD 1 (Ns/m) 210.92 Stiffness of TMD 2 (N/m) 5418.04 Damping of TMD 2 (Ns/m) 169.61 Table 8: TMDs parameters obtained by optimization Scenario 4 In the last scenario, two TMDs were inserted in the structure, one on the top floor and another on the penultimate one, located at nodes 183 and 178. Again, the optimal dissipator’s design is performed for the mass ratio defined in scenario 2, with a value of 3%, this mass being divided equally between the two TMDs. WOA parameters remained the same as in the previous scenario for all analyzes performed in this fourth scenario. As in the previous case, the number of iterations was evaluated through the convergence of the optimal answer to the problems, being possible to verify that the convergence was reached in the 25th iteration. The convergence curve was analyzed in all simulations carried out in this research in order to guarantee the optimal design for each situation. The concrete volume values obtained in the simulations are shown in Tab. 9. From the table, it can be seen that the concrete volume results are very similar to the results of scenario 3, which also uses 2 TMD, with a difference of only 2% in the average volumes. Consequently, there was a reduction in concrete volume in relation to the original system of approximately 19.8%, an interesting value for the engineering practice. The stiffness and damping values of the TMD’s are in Tab. 10. J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 341 Model Concrete Volume (m³) Average Volume (m³) Original 260.72 260.72 Otim. 1 214.63 209.21 Otim. 2 219.86 Otim. 3 202.25 Otim. 4 213.93 Otim. 5 204.63 Otim. 6 202.3 Otim. 7 205.62 Otim. 8 210.13 Otim. 9 216.94 Otim. 10 201.78 Table 9: Comparison between concrete volume values from scenario 4 simulations. Parameters Average Value Stiffness of TMD 1 (N/m) 10008.48 Damping of TMD 1 (Ns/m) 191.52 Stiffness of TMD 2 (N/m) 5052.47 Damping of TMD 2 (Ns/m) 190.18 Table 10: TMDs parameters obtained by optimization Simulation Without TMD 1 TMD 2 TMD’s (One Floor) 2 TMD’s (Two Floors) 1 226.14 187.88 192.63 201.78 2 219.22 187.36 214.17 214.63 3 232.26 201.16 201.96 219.86 4 230.61 205.38 199.67 202.25 5 229.1 212.44 207.76 213.93 6 217.98 199.11 216.62 204.63 7 225.39 201.5 202.3 202.3 8 215.72 206.42 214.37 205.62 9 213.49 203.4 203.4 210.13 10 223.75 189.1 205.6 216.94 Average 223.37 199.38 205.85 209.21 Initial Volume: 260.72 Difference (%) 14.33 23.53 20.48 19.76 Table 11: Overall results of optimization processes J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 342 Again, there is a considerable difference between the stiffness values of the two TMDs. However, the damping values are very close, an interesting point for installing the attenuators. It should be noted that, for all simulations performed with TMDs, all restrictions on displacement between floors (story drift) and maximum displacement at the top were respected. Among the 3 scenarios with the presence of TMD’s, this scenario was the one with the highest average volume value, and scenario 1 with the lowest concrete volume values. In a way, scenario 1 was expected to be the best scenario, since the solution with 1 TMD on top is a classic solution and widely used by researchers. A summary of all results is shown in Tab. 11. CONCLUSIONS n this work, procedures for the optimization of a concrete building through a metaheuristic algorithm were presented, from a 2D frame structure of 42 floors, totaling approximately 105 meters high and 4 spans of beams, susceptible to vibrations caused by wind. Computational routines for dynamic analysis and simulation of wind force were created (comprising the action of atmospheric turbulence in the wind direction and without considering neighborhood effects), as well as a routine to obtain the dynamic response from the method of integration of Newmark. These routines were developed in Octave language and integrated with the WOA optimization algorithm. Several simulations were performed, and in all of them, smaller values of concrete volume were obtained in relation to the initial structure, which shows the effectiveness of the method proposed in this work. Given the simulations carried out in this work and the results obtained, it is clear that the procedure proposed in this paper to optimize concrete structures through metaheuristic algorithms is satisfactory, bringing good results in relation to the volume of concrete (reduction of up to 24% in relation to original structure), and also in relation to the maximum displacement required by the standard, respecting the limit value. That is, the original project, in addition to having a larger concrete volume, did not respect the maximum displacement requirements of standards, however, after the proposed optimizations, in all scenarios, the concrete volume was reduced and the structure began to respect the code limits. Regarding the insertion of TMDs, all insertion scenarios are positive and reduced the volume of concrete in a greater quantity than the optimization without TMD. The greater efficiency of scenario 2 is highlighted, which is the type of control most commonly found in the literature to reduce building displacement amplitude. Scenarios 3 and 4 are also efficient, despite the different parameters of each attenuator, which can make the TMD installation process more expensive. ACKNOWLEDGMENTS he authors acknowledge the financial support of CNPq and CAPES, Brazil. REFERENCES [1] Thornton, C.H., Joseph, L., Scarangello, T. (1990). Optimization of tall structures for wind loading, J. Wind Eng. Ind. Aerodyn., 36(1–3), pp. 235–244, DOI: 10.1016/0167-6105(90)90308-Y. [2] Venanzi, I., Materazzi, A.L. (2007). Multi-objective optimization of wind-excited structures, Eng. Struct., 29(6), pp. 983– 990, DOI: 10.1016/j.engstruct.2006.07.007. [3] Yang, J.N., Agrawal, A.K., Samali, B., Wu, J.C. (2004). A benchmark problem for response control of wind-excited tall buildings, Proc. Int. Modal Anal. Conf. - IMAC, 1(April), pp. 151–157, DOI: 10.1061/(asce)0733-9399(2004)130:4(437). [4] Miguel, L.F.F., Miguel, L.F.F. (2012). Novel metaheuristic algorithms applied to optimization of structures, WSEAS Trans. Appl. Theor. Mech., 7(3), pp. 210–220. [5] Miguel, L.F.F., Miguel, L.F.F. (2013). Assessment of modern metaheuristic algorithms - HS, ABC and FA - in shape and size optimisation of structures with different types of constraints, Int. J. Metaheuristics, 2(3), pp. 256, DOI: 10.1504/ijmheur.2013.056404. [6] Holland, J.H. (1992). Adaptation in Natural and Artificial Systems: An Introductory Analysis with Applications to Biology, Control, and Artificial Intelligence, DOI: 10.7551/mitpress/1090.001.0001. I T J. W. S. Brito et alii, Frattura ed Integrità Strutturale, 59 (2022) 326-343; DOI: 10.3221/IGF-ESIS.59.22 343 [7] Lei, Z., Gao, S., Gupta, S., Cheng, J., Yang, G. (2020). An aggregative learning gravitational search algorithm with self- adaptive gravitational constants, Expert Syst. Appl., 152, pp. 113396, DOI: 10.1016/j.eswa.2020.113396. [8] Yang, X.S. (2009). Harmony search as a metaheuristic algorithm, Stud. Comput. Intell., 191, pp. 1–14, DOI: 10.1007/978-3-642-00185-7_1. [9] Marini, F., Walczak, B. (2015). Particle swarm optimization (PSO). A tutorial, Chemom. Intell. Lab. Syst., 149, pp. 153– 165, DOI: 10.1016/j.chemolab.2015.08.020. [10] Dorigo, M., Birattari, M., Stutzle, T. (2006). Ant colony optimization, IEEE Comput. Intell. Mag., 1(4), pp. 28–39, DOI: 10.1109/MCI.2006.329691. [11] Mirjalili, S., Lewis, A. (2016). The Whale Optimization Algorithm, Adv. Eng. Softw., 95, pp. 51–67, DOI: 10.1016/j.advengsoft.2016.01.008. [12] Longarini, N., Cabras, L., Zucca, M., Chapain, S., Aly, A.M. (2017). Structural Improvements for Tall Buildings under Wind Loads: Comparative Study, Shock Vib., DOI: 10.1155/2017/2031248. [13] Saaed, T.E., Nikolakopoulos, G., Jonasson, J.E., Hedlund, H. (2015). A state-of-the-art review of structural control systems, JVC/Journal Vib. Control, 21(5), pp. 919–937, DOI: 10.1177/1077546313478294. [14] Symans, M.D., Constantinou, M.C. (1999). Semi-active control systems for seismic protection of structures: A state-of- the-art review, Eng. Struct., 21(6), pp. 469–487, DOI: 10.1016/S0141-0296(97)00225-3. [15] Varela, W.D., Battista, R.C. (2011). Control of vibrations induced by people walking on large span composite floor decks, Eng. Struct., 33(9), pp. 2485–1494, DOI: 10.1016/j.engstruct.2011.04.021. [16] Rossato, B.B., Miguel, L.F.F. (2019). Optimization of parameters of tuned mass dampers for use in tall buildings subjected to the wind action, XL CILAMCE Ibero-Latin Am. Congr. Comput. Methods Eng., pp. 12. [17] Liu, Y., Wang, K., Mercan, O., Chen, H., Tan, P. (2020). Experimental and numerical studies on the optimal design of tuned mass dampers for vibration control of high-rise structures, Eng. Struct., 211, pp. 110486, DOI: 10.1016/j.engstruct.2020.110486. [18] Den Hartog, J.P. (1956). Mechanical Vibrations, New York: McGraw-Hill. [19] Warburton, G.B. (1982). Optimum absorber parameters for various combinations of response and excitation parameters, Earthq. Eng. Struct. Dyn., 10(3), pp. 381–401, DOI: 10.1002/eqe.4290100304. [20] ABNT. (2014). ABNT NBR 6118:2014 - Projeto de estruturas de concreto, Proj. Estruturas Concreto, Procedimento. Norma Bras., pp. 238. [21] Associação Brasileira de Normas Técnicas. (1988). Forças devidas ao vento em edificações: NBR 6123, Rio de Janeiro. [22] Shinozuka, M., Jan, C.M. (1972). Digital simulation of random processes and its applications, J. Sound Vib., 25(1), pp. 111–128, DOI: 10.1016/0022-460X(72)90600-1. [23] Blessmann, J. (2013). O Vento na Engenharia Estrutural, Porto Alegre: Editora UFRGS. [24] Miguel, L.F.F., Fadel Miguel, L.F., Riera, J.D., Kaminski, J., Ramos de Menezes, R.C. (2012). Assessment of code recommendations through simulation of EPS wind loads along a segment of a transmission line, Eng. Struct., 43, pp. 1–11, DOI: 10.1016/j.engstruct.2012.05.004. [25] Brandão, F. da S., Miguel, L.F.F. (2020). Vibration control in buildings under seismic excitation using optimized tuned mass dampers, Frat. Ed Integrita Strutt., 14(54), pp. 66–87, DOI: 10.3221/IGF-ESIS.54.05. [26] American Society od Civil Engineers. (2010). Minimum Design Loads for Buildings and Other Structures: ASCE/SEI 7-10, Reston, Virginia. << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice