Design of a Tuned Mass Damper Inerter (TMDI) Based on an Exhaustive Search Optimization for Structural Control of Buildings under Seismic Excitations

A tunedmass damper inerter (TMDI) is a new class of passive control device based on the inclusion of an inerter mechanism into a conventional tuned mass damper (TMD). (e inerter device provides inertial resisting forces to the controlled system, through relatively small masses, converting it in a mechanism with the potential to enhance the performance of passive energy dissipating systems. (is work presents a study of an optimal TMDI design through an exhaustive search process. TMDI device design using the cited parameter selection methodology consists in the determination of the damper critical damping ratio, ζTMDI, and frequency ratio, υTMDI, which result in the minimum structural response of a multidegree of freedom structural system, considering predefined values for mass ratio (μ) and inertance ratio (β). (e used optimization process examines all possible damping device design parameter combinations to select the set of values that results in the best device performance to reduce response parameters in a structure. Four different optimization processes are performed by independently minimizing four performance indices: J1 associated to the reduction of the structure’s maximum peak displacement, J2 calculates the minimal RMS value for the structure’s peak displacement, J3 seeks by the minimal peak interstorey drift, and JP determines the lowest value for a linear weighted combination of the abovementioned three indices. A numerical example is developed with the purpose of validating the proposed optimization procedure and to evaluate the benefits of using TMDI as controlling devices for structures under seismic excitation, by carrying out a comparative analysis to contrast the performance of the optimization alternatives developed, running up to 1968192 cases.(e obtained results show that devices designed based on exhaustive search optimization produce peak displacement reductions of up to 35% and peak structure displacement RMS reductions of up to 30%.


Introduction
Structural control is, fundamentally, a set of techniques that aims to reduce risks associated with the exposure of structures to random events originated in natural phenomena such as earthquakes, wind, and the effects of coastal waves. e main motivation for structural control is the possibility of providing with high level safety, comfort, and peace of mind to the people that inhabit the controlled structures. Structural control may be achieved through the use of diverse types of devices, mechanisms, and algorithms with varied levels of structural behavior. Successful structural control methodologies have been developed based on energy to the structural system, while in passive devices, their behavior results in the use inertial forces within the structure achieving reduction in the structural response through energy dissipation mechanisms without the need for external power sources. Semiactive and hybrid devices correspond to intermediate solutions in which the working mechanisms use characteristics of both active and passive devices with low energy consumption.
Due to their simplicity, reliability, economy, and sturdiness, passive control devices stand out among structural control alternatives, especially the devices known as tuned mass dampers (TMDs). ese mechanisms, composed by a mass, a spring, and a damper, are some of the oldest vibration control devices in existence [21]; moreover, their proved efficiency for the damage reduction and collapse prevention of structural systems [22], as they are tuned to the fundamental frequencies of the structure, make them ideal mechanisms to protect the structural integrity of a building. Some of the most representative applications of TMD can be found in base isolation [23][24][25][26], seismic protection of buildings [27][28][29][30], effect in soilstructure interaction [31][32][33], and vibration control of pedestrian bridges [34,35], among others. According to their functioning mechanism, TMD are classified into conventional TMD, pendulum TMD (PTMD), bidirectional TMD (BTMD), and liquid column tuned dampers (TLCD) [36].
Searching for improvement in the performance of these devices, an important number of researchers have been working during the last five years with a promising modification of conventional TMD, consisting in the addition of a two-node linear device, called an inerter, that develops an internal resisting force proportional to the relative acceleration between the two nodes [37].
is new device composed by an inerter and a conventional TMD has been called the tuned mass damper inerter (TMDI). Several research works have explored the design [38][39][40], optimization [39][40][41][42], and performance [43][44][45] of this new passive control device with positive results, proving that these mechanisms may outperform the conventional TMD output in structures under seismic loads [46]. Likewise, inerterbased devices have been used in specific applications to different fields of civil engineering. Examples of this are the study of TMDI for the mitigation of wind-induced response of tall buildings [47], the use of a Tuned Damper Inerter (TID) for suppression of unwanted cable vibrations due to support excitations [48], the numerical implementation of a Tuned Parallel Inerter Mass System (TPIMS) for seismic mitigation of a wind turbine tower [49], the optimal design based on analytical solution for a liquid storage tank with a inerter isolation system [50], the optimal design and seismic performance of a TMDI when installed below the isolation floor of base-isolated structures in order to enhance the earthquake resilience and reduce the displacement demand [51], the development of a comfort-based optimal method for designing floors using Tuned Inerter Mass Systems (TIMS) to reduce the vertical-vibration response of floors subjected to human-induced excitations [52], the study of a Tuned Inerter Liquid System (TILS) for vibration control by utilizing the advantages of the inerter-based system and tuned liquid element to improve control performance [53], and the investigation of the effect of connecting TMDIs between nonconsecutive floors [54], among others.
Regarding the quoted paper by Giaralis and Petrini [54], it is suggested that an enhanced behavior is found when connecting the TMDI device located at the top floor with a floor up to three levels below, although reasonable, was considered not applicable to the work described herein, as a retrofit application is exemplified in this paper.
e main objective of the study presented here is to obtain an optimal TMDI design through an exhaustive search process, which consists in the determination of the damper critical damping ratio, ζ TMDI , and frequency ratio, υ TMDI , which result in the minimum structural response of a multidegree of freedom structural system, considering predefined values for mass ratio (µ) and inertance ratio (β). e proposed optimization process examines all possible damping device design parameter combinations, to select the set of values that results in the best device performance to reduce response parameters in a structure. Four different optimization processes are performed by independently minimizing four performance indices: J 1 is associated to the reduction of the structure's maximum peak displacement, J 2 calculates the minimal RMS value for the structure's peak displacement, J 3 seeks by the minimal peak interstorey drift, and J P determines the lowest value for a linear weighted combination of the abovementioned three indices. A numerical example is developed with the purpose to validate the proposed optimization procedure and to evaluate the benefits of using TMDI as controlling devices for structures under seismic excitation, by carrying out a comparative analysis to contrast the performance of the optimization alternatives developed, running up to 1968192 cases.
is research represents a thorough and rigorous, albeit simple, optimization process applied on a model with multiple degrees of freedom. e exhaustive search-based optimization technique works efficiently in limited size problems, such as the one at hand. Furthermore, exhaustive search is reliable, easy to implement, and represents a certainty, as all possible solutions are investigated.

Tuned Mass Damper Inerter (TMDI)
e inerter device ( Figure 1) is conceptually developed based on analogies between existing mechanical and electrical systems, as a true mechanical equivalent of the electric condenser [55]. is mechanism, characterized by a small mass in comparison with the storey mass, is an element connected to two nodes. ese nodes move, respectively, from each other, generating relative accelerations that are proportional to the resisting forces produced by the device, as described by equation (1), where F is the resisting force in each node, b is the proportionality constant in terms of mass which defines the inerter, and € u i and € u j are the i and j nodal accelerations, respectively: Most proposed inerters are devices based on mechanisms developed from gears, racks, pinions, and flywheels 2 Shock and Vibration [38,39,45,56], although there are also inerter proposals based on pendulum systems [43], on ball screws [56], and on hydraulic devices [57][58][59].

Equation of Motion for an n Degree of Freedom System
Equipped with a TMDI. To pose the equations of motion that govern an n degree of freedom system equipped with a TMDI, consider the structure shown in Figure 2. Each storey stiffness, damping, and mass is represented with a spring with stiffness constant k i , viscous damper with damping constant c i , and mass concentration elements m i , respectively, where i symbolizes each level of the structure. Likewise, the TMDI device is idealized in this figure as the combination of a conventional TMD and an inerter. In this case, the conventional TMD is installed in the n th level of the building; the inerter connects the conventional TMD with the structure's level n − 1. e TMDI mass is denoted by m TMDI , while the stiffness and damping parameters associated with the damper are denoted as k TMDI and c TMDI , respectively. Constant b is the parameter quantifying the TMDI inertance.
In this way, the dynamic behavior of the n+1 degrees of freedom TMDI equipped system subjected to a ground acceleration at its base may be modelled with equation (2), where x(t), _ x(t), and € x(t) are the n + 1 displacement, velocity, and acceleration system vectors, respectively, which depend on the variation of time, 1 is a unitary vector of the order n + 1, and € x g (t) is an scalar representing the ground acceleration which also varies with time.
M is the system's mass matrix of the order (n + 1) x (n + 1), shown with the following expression: C refers to the system's damping matrix, of the order of (n + 1) x (n + 1), expressed as Finally, the system's stiffness matrix K, of the order (n + 1) x (n + 1), may be calculated as e motion described with equation (2) may be transformed into a space-state representation in order to more readily solve the dynamic system through first-degree differential equations. Consequently, in order to represent the n Shock and Vibration 3 degree of freedom system's dynamic behavior as a spacestate form, the state vector z(t) is defined as where z(t) is a vector of the order 2(n + 1) x 1 containing the displacements and velocities from the structural system response. us, equation (2) may be transformed into equation (7), where 0 and I are the null and the identity matrices, respectively, of the order (n + 1) x (n + 1).
Now, in order to simplify equation (7), the matrices A and B are introduced as per equations (8) and (9), where A is the transition state matrix, of the order 2(n + 1) x 2(n + 1), and B is the structural system's external excitation action location matrix: us, equation (7) becomes e solution for this differential linear equation of the first-degree gives the dynamic behavior of the structure equipped with the TMDI device, when subjected to an acceleration record at its base.

TMDI Design Optimal Parameters
Selection Methodology e optimal parameters for the design of the TMDI device are determined in this work through an optimization process based on the so-called exhaustive search technique. is optimization methodology is a general optimal parameter search technique that does not require specific knowledge of the dominion, as only a general state description, an operator set adjusted to the system's typical behavior, an initial state, and a goal or objective for the final state are needed [60]. e exhaustive search-based optimization methodologies work simply and directly, resulting in adequate and precise solutions, but require an important amount of processing time and power. Allowing for this, an exhaustive search algorithm examines all possible solutions of each particular state, continually comparing which of these possible solutions is most efficient in complying with the problem's restrictions with the lowest cost, which is usually defined with the system's own characteristics and the imposed conditions by the problem at hand. e main motivation for the use in this paper of the exhaustive search based optimization technique is that, in limited size problems, such as the one at hand, for the selection of TMDI device parameters, the used methodology works efficiently. is must be added to the fact that exhaustive search is reliable and easy to implement when it is compared with others optimization techniques such as evolutionary algorithms, gradient-based optimization algorithms, and nature-inspired metaheuristics; furthermore, trying all possible solutions, determining an optimal state through the comparison between different final states of the problem under study, is a certainty.
TMDI devices design, using the described parameter selection procedure, consists in determining the values for the damper critical damping ratio, ζ TMDI , and frequency ratio, υ TMDI , which are defined according to equations (11) and (12), respectively, where c TMDI , m TMDI and ω TMDI are the damping, mass and frequency TMDI device parameters, respectively, b is the inertance, and ω 1 is the structural system's natural fundamental frequency: e methodology starts with closed continuous sets of possible optimal values for parameters ζ TMDI and υ TMDI . To begin with, the process generates a division of these sets into a finite number of values, comprising the complete set of optimal values for each parameter. e optimization method uses all possible combinations of finite values generated in the process of division of the closed sets, so each one of these discrete values and its respective combination are used to calculate the TMDI damper coefficients of mass (m TMDI ), damping (c TMDI ), and stiffness (k TMDI ), from the mathematical manipulation of expressions (11) and (12), assuming a device mass ratio (µ) and an inertance ratio (β) whose values are constant through the whole optimization process, where µ and β are defined according to expressions (13) and (14) and m TMDI and m 1 are the TMDI device mass and the generalized system mass that allows to obtain the structural system natural fundamental frequency, respectively: Once the TMDI parameters are determined, the optimization procedure calculates the controlled structural system response when subjected to a specific acceleration record. is response is stored within a database where it is processed to determine the performance parameters that allowed establishing comparison criteria between the diverse examined device cases. After this step, the optimization process selects a new discrete value for ζ TMDI and υ TMDI and repeats the process anew until it obtains updated comparison parameters according to the new used characteristics. When the exhaustive search algorithm has spent all possible solutions to the optimization problem, a comparison routine defines which are the TMDI design parameters which result in the optimized structural response, measured through the minimum values of three response parameters indices, namely, the structure's maximum peak displacement, the RMS value for the structure's peak displacement, and the peak storey drift: Furthermore, a fourth performance index is included, consisting in a linear combination of the first three indices, with specific weights, according to the relevance given to each response parameter; this linear combination is labeled as performance index, J P , and is described with the following equation: ese particular weights for the combination of the performance indices were defined so as to amplify the reduction of the peak displacement, although all three performance indices show reduced responses. Figure 3 shows a logic tree diagram describing the exhaustive search algorithm parameter definition used in this study.
is is the procedure used to calculate the optimal TMDI device design parameters, which controls the response of a structure subjected to seismic excitation.

Numerical Example for a TMDI-Equipped n Degree of Freedom System Optimally Designed Using the Exhaustive Search Optimization Procedure
To study the performance of a TMDI controlled system designed with the proposed exhaustive search procedure, a numerical model based on a real structure is implemented. is model is subjected to four acceleration records of natural origin, and a comparative analysis of the control device is developed. e model used in this example is based on the building that houses Beneficencia de Antioquia, which is the governmental agency in charge of the city of Medellin's lottery, Colombia. is building has 11 levels, each with a height of 2.72 m. e structure's total height is 32.14 m above the ground level, counting the elevators' engine room and the roof. According to the Colombian earthquake resistant design code, Normas Colombianas de Diseño y Construcción Sismo Resistente, NSR-10 [61], the structure is comprised by two structural systems. In the north-south direction, parallel to the building's main facade, the structure is a reinforced concrete momentresisting frame, while in the east-west direction, the structure is a combination of a moment-resisting frame and structural walls. Figure 4 presents an image of the modelled building; Figure 5 shows the fifth floor plan, which is a typical floor configuration, and the plane frame modelled.
e H-shaped structure was used because it is part of a larger study for seismic retrofit being conducted by the research group; however, a plane frame was derived from the original structural system; the degrees of freedom of this plane frame were reduced to one degree of freedom per storey, assuming an infinitely rigid floor diaphragm so as to obtain equal horizontal degrees of freedom for all nodes in each storey, assuming axially infinitely rigid column so as to eliminate all vertical degrees of freedom for all nodes and, finally, condensing the rotational degrees of freedom as to be expressed in terms of the remaining horizontal degrees of freedom.
Expressions (17) and (18) represent the stiffness (K) and mass (M) matrices of the reduced model, respectively. e damping matrix (C) is calculated as a Rayleigh damping matrix proportional to the stiffness and mass matrices, assuming a critical damping ratio of 5% for the first and last modes of vibration. Table 1 shows the modal properties of the plane frame modelled.
Once the structure's model characteristics are defined, the TMDI parameter optimization procedure is initialized.
is algorithm employs two fixed mass ratios, µ, of 5% and 2% and uses two close sets of optimal values for parameters ζ TMDI and υ TMDI , defined as ese value intervals have been selected because they are coherent with conventional TMD devices, serving as starting points as TMDI has the dynamic behavior same as that of conventional TMD [62].
Several inertance ratios values, β, are tried in the optimal TMDI design in order to analyze the impact of that property on the structural system response performance. With that purpose, β values of 5, 10, 20, 30, 40, and 50% were used. Both the structural model and the optimization tool were  implemented in MATLAB using routines built by the authors. Table 2 shows the acceleration records used to excite the controlled system in order to evaluate the performance of the different TMDI configurations in the optimization process. Multiple records commonly found in the literature were used to encompass diverse acceleration inputs to show that the device works independently of the time-history demand. Tables 3-6 show the data obtained from the structure's response under the seismic records employed in this work for the system uncontrolled; interesting values are worth pointing out. e first one has to do with the maximum displacement magnitude, which reached 0.572 m, corresponding to the lateral movement of the building's 11 th floor under the "Kobe" earthquake record. Secondly, the maximum displacement RMS corresponds also to the same floor, but under the acceleration record from "El Centro," which reached 0.162 m, representing 165% greater displacement RMS of the 11 th floor under the "Kobe" record, which reached a value of 0.097 m. With respect to the storey drift data, the peak value was reached under the "Kobe" record for the horizontal relative displacement between floors 8 and 9, of 0.107 m, constituting a storey drift peak value of 3.92%. Under the "Algarrobo," "El Centro," and "Christchurch" records, the storey drift peak values were 0.083 m, 0.078 m, and 0.061 m, respectively, values that represent a storey drift of 3.05%, 2.86%, and 2.23%, respectively.

Response of the Reference System (without Control).
Only elastic analysis was considered in this study as in most codes, the principle of equivalent displacements is invoked to allow for elastic analysis for design, using a reduction factor for loads to consider the effect of plastification. e purpose of TMD and TMDI practical applications is to reduce the structural response to minimize plastification on the building frame and concentrate energy dissipation on the controlling device. In this paper, we seek to try a thorough parameter optimization process to produce reduced response; whether the final response lies inside or outside the linear behavior is out of the scope of this paper.

TMDI Optimization Process.
e TMDI design parameters are obtained for four different optimization alternatives, corresponding to the evaluation of indices J 1 , J 2 , J 3 , and J P . To make sure the device design parameter selection corresponds to optimal parameters, a comparative analysis of the controlled structure's behavior is made for each of the TMDI parameter optimization alternatives. In this way, the comparative analysis shows the general control device performance main tendencies found with each of the optimization alternative studied.
Tables 7-10 present a summary of the parameters ζ TMDI and υ TMDI obtained with the optimization processes conducted for TMDI mass ratios (µ) of 2 and 5% and inertance ratios (β) of 5, 10, 20, 30, 40, and 50%. Figure 6 presents the optimal values distribution dispersion for parameters ζ TMDI and υ TMDI obtained through the linear combination optimization proposal (J P ), using the mass ratios and the inertance ratios values mentioned above, for the acceleration records considered. Table 11 presents the response parameters for the controlled structure subjected to the "El Centro" earthquake record for the four optimization alternatives proposed in this study. Positive and negative values in the percentage reduction columns in this table represent, respectively, a response reduction and response increase. Based on the analysis of the results data, it may be said that the TMDI controlled system performance with the greatest mass ratio (µ � 5%) is always superior to that of the lesser mass ratio (µ � 2%). e greatest peak displacement reduction of the structure, corresponding to the 11 th floor displacement, is 6.76%, for the TMDI with µ � 2% and β � 5% and optimizations alternatives J P and J 1 ; this same parameter for the same optimizations alternatives but with device values of µ � 5% and β � 5% reaches a reduction of 25.84%. Something similar occurs with the RMS value for the 11 th floor displacement, where the greatest reduction for this parameter is obtained for the TMDI designed with values µ and β of 5%, with reductions of 46.64% for the optimization alternative J 2 , while the best control alternative for TMDI designed with µ � 2%, corresponding also to the optimization alternative J 2 , is that with β � 5%, with RMS displacement reductions for the 11 th floor of 24.47%. Finally, the peak interstorey drift reductions show that the greater the mass ratio, the greater the response parameter reduction, as the TMDI-controlled system with µ and β of 5% reduces the maximum peak storey drift in 26.22%, while the TMDI with µ � 2% and β � 5% results in reductions of 12.87%.

Controlled System Response Using the Optimization Alternative JP.
us, it can be seen that the device with mass and inertance ratio values of 5% is the TMDI with the best results in terms of response parameter reduction. Likewise, the linear combination optimization alternative, J P is, in the case of the controlled structure subjected to the "El Centro" earthquake, the most balanced and interesting alternative, as it results in important reductions of peak displacements and of peak storey drifts. e results clearly show that the linear combination optimization surpasses the single objective alternatives, J 2 and J 3 , characteristics, generally reaching considerable reduction for all analyzed response parameters. Figure 7 presents a graphic comparison of the 11 th floor displacement evolution of the structure controlled with the TMDI configured with design parameters selected by the optimization alternative J P for µ � 5% and for the different values of inertance ratio used in this work. Figure 8 shows the peak storey drift ratio for each level of the building controlled by a TMDI with µ � 5% and β between 5 and 50%.
TMDI performance optimized for the structure subjected to the "Christchurch" record confirms that inertance ratio (β) and mass ratio (µ) values of 5% are the parameter combination that achieves the best response reduction results when optimized using the linear combination alternative proposed in this work. For this specific load case, the peak displacement reductions are of up to 10%, while the 8 Shock and Vibration      Table 7: TMDI optimization parameters for µ � 0.02 and µ � 0.05 for the structure subjected to the "El Centro" record.  Table 8: TMDI optimization parameters for µ � 0.02 and µ � 0.05 for the structure subjected to the "Christchurch" record.  Table 9: TMDI optimization parameters for µ � 0.02 and µ � 0.05 for the structure subjected to the "Kobe" record.  Table 10: TMDI optimization parameters for µ � 0.02 and µ � 0.05 for the structure subjected to the "Algarrobo" record.   displacement RMS achieves up to 28% reductions. In general terms, the displacement reductions, for both the peak and the RMS values, present a somewhat uniform performance tendency, where small response reduction percentages are obtained for the structure's lower levels and much greater response reduction values for the last floors of the building, which is where the TMDI device exerts the greatest influence over the structural system. Figures 9 and 10 show these tendencies, for the peak displacements and the RMS displacements, respectively, for each one of the levels of the building, for the most efficient TMDI-controlled case, which is the one corresponding to the device with β and µ of 5% optimized with the linear combination alternative J P . Tables 12 and 13 present the structure's response under the "Christchurch" and "Kobe" records, respectively, for the TMDI device parameters obtained with the different optimization alternatives used in this work. Figure 11 presents the storey peak acceleration for the controlled structure subjected to the "Kobe" record for the optimization alternative J p and for a value of µ � 0.05. Although all cases show important reduction for this response parameter, this time, the most efficient alternative does not correspond to the one with β of 5%, but for the optimized   Figure 8: Interstorey drift ratio of the controlled structure subjected to "El Centro" acceleration record for TMDI designed with µ � 5% and the J P optimization alternative.    Figure 9: Peak displacement of each storey of the controlled structure subjected to "Christchurch" acceleration record for TMDI designed with µ � 5% and the J P optimization alternative. Uncontrolled β = 5% β = 10% β = 20% β = 30% β = 40% β = 50% Figure 10: RMS value of the displacement of each storey of the controlled structure subjected to "Christchurch" acceleration record for TMDI designed with µ � 5% and the J P optimization alternative.        Figure 11: Peak acceleration of each storey of the controlled structure subjected to "Kobe" acceleration record for TMDI designed with µ � 5% and the J P optimization alternative. Uncontrolled β = 5% β = 10% β = 20% β = 30% β = 40% β = 50% Figure 12: RMS value of the acceleration of each storey of the controlled structure subjected to "Kobe" acceleration record for TMDI designed with µ � 5% and the J P optimization alternative.  TMDI with β of 50%. In fact, for the peak storey acceleration, the linear combination optimized, J P , TMDI with µ and β of 5% is the case with smaller reduction values of all the optimization alternatives analyzed; the same situation may be seen in Figure 12 for storey acceleration RMS values. Figure 8 shows trends characteristic of drift variation with respect to height in a moment-resistant frame, with maximum drifts occurring at the third lower part of the structure's total height. Similar trends are shown in Figures 11  and 12 for floor accelerations, which are also characteristic of acceleration distribution for ultimate lateral loading, where accelerations tend to decrease with height, with respect to ground acceleration, except for the floors of the upper quarter of the structure, which tend to increase markedly. By contrast, Figures 9 and 10 show the absolute displacement distribution, which characteristically increase with height, with a higher rate of increase in the middle third of the structure's total height. Table 14 presents the controlled structure's response values obtained for the different optimization alternatives developed, when subjected to the "Algarrobo" record.
For this case, the controlled system's response reduction, under the Algarrobo record, shows that the inerter device with mass ratio and inertance ratio of 5% achieves the best performance. Figure 13 shows the greatest peak displacement reduction for µ and β of 5%; however, as the inertance ratio, β, increases, the reduction percentage diminishes, that is, the peak displacements increase, almost linearly for a mass ratio, μ, of 5%, and asymptotically for µ of 2% Figure 14 presents the reduction of the 11 th storey RMS value of displacement of the structure subjected to "Algarrobo" acceleration record. For this case, results also show that the optimized device using the linear combination alternative achieves the best results for µ of 5%, but for this input seismic record, the greatest reduction is achieved in conjunction with β of 20%, although the variation of the reduction percentage for µ of 5% does not vary beyond 20% of the maximum as β varies. For TMDI designed with µ of 2%, reductions are highest for β of 5%, and reduction decreases asymptotically as β increases.
Finally, Figure 15 shows the peak storey drift reduction for the controlled structure under the "Algarrobo" record. For this case, results also show that the optimized device using the linear combination alternative achieves the best results for µ and β of 5%, with reductions decreasing significantly for β of 10%, and with no significant variation as β increases further. For µ of 2%, reductions are the highest for β of 5%, and reduction decreases asymptotically as β increases. e response results shown in Figure 11 through 15 show that the TMDIs effectively reduces the structural system response for all combinations of µ and β; moreover, the graphs show that the highest reductions are obtained for peak interstorey drift with µ and β of 5%. at is, the drift is effectively reduced for all inertance ratio values, which implies a reduction of damages under seismic excitations.

Conclusions
is paper deals with the optimal design of TMDI using exhaustive search. e used methodology focuses on determining the best combination of parameters, so the optimal device is able to reduce the response of a structural system, through the single and combined control of three different response parameters: structure's maximum peak displacement, RMS value for the structure's peak displacement, and peak interstorey drift. e tuned mass damper inerter designed optimally have a significant potential for the response reduction of multiple degree of freedom systems under seismic excitations. e obtained results show that those devices with fixed mass ratio (μ) and  inertance ratio (β) of 5% achieved the greatest last floor reductions of peak displacements and RMS displacement values, for mean values of critical damping ratio (ζ TMDI ) and frequency ratio (υ TMDI ), of 0.03 and 0.89, respectively, with reductions of up to 25%, for peak displacement, and up to 46% for RMS displacement values. For these design parameters, the maximum peak storey drift ratio reductions are of the order of 26%. Based on these results, it may be said that the linear combination optimization alternative is a valuable tool of easy implementation for the determination of optimal design parameters for TMDI devices in structural systems subjected to seismic actions. Nonlinear analysis inclusion on this topic should be the subject of future developments on this field to ensure practical applications of TMDI devices.

Data Availability
All data included in this study are available upon request by contact with the corresponding author.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this article.  Figure 14: Comparison of the reduction of the 11 th storey RMS value of displacement of the structure subjected to "Algarrobo" acceleration record for TMDI designed with µ � 2% and µ � 5% using the J P optimization alternative.  Figure 15: Comparison of the reduction of the peak interstorey drift of the structure subjected to "Algarrobo" acceleration record for TMDI designed with µ � 2% and µ � 5% using the J P optimization alternative.  Figure 13: Comparison of the reduction of the 11 th storey peak displacement of the structure subjected to "Algarrobo" acceleration record for TMDI designed with µ � 2% and µ � 5% using the J P optimization alternative.