Human-Induced Vibration Control of Floor Structures Using MTMD System Optimized by MATLAB-SAP2000 Interface

: Under human-induced excitations, a floor structure may suffer excessive vibrations due to its large span and low damping ratio. Vertical vibrations, in particular, can become intolerable during resonance events. A tuned mass damper (TMD) is a widely used single-degree-of-freedom dynamic vibration absorber. To enhance the serviceability of a floor structure, a multiple TMD (MTMD) system finds broad application. The parameters of the MTMD must be carefully designed to achieve satisfactory performance. However, existing studies often employ a simplified model of the floor structure with closely spaced modes to optimize the parameters of MTMD. Nonetheless, an oversimplified floor model can lead to a reduction in its control effect. To solve this problem, this study utilizes the OAPI facility of SAP2000 to build a connection with MATLAB. A multi-objective optimization algorithm based on the artificial fish swarm algorithm (AFSA) for MTMD is developed in MATLAB, while the finite element model of a real floor structure is built in SAP2000. The locations of the MTMD system are initially specified in SAP2000 and, through the proposed MATLAB–SAP2000 interface, data can be exchanged between them. Based on the structural dynamic responses to external excitations in SAP2000, the optimization process for the MTMD is carried out in MATLAB. Concurrently, the parameters of the MTMD in SAP2000 are iteratively adjusted until they reach their final optimal values. To underscore the enhancements brought about by the proposed interface and optimization method, a comparative case study is conducted. A group of MTMDs, optimized using a conventional method, is presented for reference. The numerical results indicate that, overall, the proposed MTMD system exhibits superior control effectiveness and robustness.


Introduction
With the acceleration of urban modernization and the rapid development of science and technology, modern architectural structure forms are developing towards larger volumes and spans [1][2][3][4].Simultaneously, heightened demands are placed on the serviceability of structures.For sports venues where the floor structure possesses a lighter weight and lower damping ratio, the natural frequency tends to be lower [5][6][7][8][9].Considering the 1.6 Hz to 2.4 Hz human-induced excitation frequency range, when the excitation frequency closely aligns with the structural natural frequency, resonance phenomena can compromise user comfort and even evoke a sense of panic [10][11][12][13].Most recently, serviceability evaluations on low-frequency floors, cross-laminated timber floors, concrete floors, lightweight timber floors, and reinforced concrete precast one-way joist floors can be found in [14][15][16][17][18], respectively.
An appealing approach to improve the serviceability of a floor structure involves the application of additional structural control devices [19].Among these, the TMD stands out as a popular SDOF dynamic vibration absorber [20][21][22][23].In the case of a TMD with constant mass, its control effectiveness relies on both its frequency and damping ratio [24][25][26][27][28]. Therefore, meticulous optimization of TMD parameters is essential to achieve satisfactory performance [29][30][31].Zhang and Zhang [32] optimized a TMD for tall intake towers, simplifying the tower to a lumped mass model.Jimenez-Alonso et al. [33] designed a TMD for a lively footbridge, employing a simplified SDOF system for the bridge.Though the designed TMD in these two references exhibited good control effectiveness, the error caused by the simplification of the model could affect the parameter optimization of the TMD.Rong and Lu [34] developed an Abaqus-Python interface to optimize a particle TMD for a tower.In TMD optimization design, a more detailed model of the main structure leads to more accurate TMD optimization results and improved control effectiveness.Inspired by Ref. [34], an interface between two commercial software programs can be devised to address issues arising from an overly simplified model.In recent years, the application of mesh-free techniques in dynamic analyses of rods, beams, and plates has gained interest among researchers, yielding more accurate results with the same number of nodes.Keivan et al. [35][36][37], for instance, explored the prediction of the dynamic response of an elastically supported thin beam traversed by a moving mass.
Typically, a floor structure exhibits closely spaced modes, and human-induced excitations span a broad frequency range.In such cases, a MTMD system, composed of several TMDs with different frequencies, is known to yield more effective control than a single TMD [38][39][40].Wang et al. [41] proposed an adaptive-passive variable mass MTMD for a floor slab.Liu et al. [42] presented a MTMD system for a cantilevered floor.Zuo et al. [43] optimized a MTMD for an offshore wind turbine.Jahangiri and Sun [44] proposed a MTMD for an offshore floating wind turbine.Gaspar et al. [45] proposed a MTMD to control the vibration of steel-concrete composite floors.Huang et al. [46,47] presented a shape memory alloy-based MTMD for a timber floor structure.From the literature review in this paragraph, it is evident that the optimization results of MTMD parameters are contingent on an accurate and precise structural model.However, the simplification of the structural model inevitably introduces errors that can diminish the control effectiveness of the MTMD.
In the literature mentioned above, the simplification of the structural model and the optimization of MTMD parameters are predominantly proposed using the business mathematics software MATLAB [48].Simultaneously, the widely acclaimed and extensively used structural analysis program SAP2000 [49] has been a focus of research, particularly in the realms of structural optimization and model updating based on the OAPI facility of SAP2000 connected to MATLAB.Through the developed MATLAB-SAP2000 interface, Dinh-Cong et al. [50] updated the model and assessed the damage to full-scale structures.Minh et al. [51] identified the damage in a 3D transmission tower.Artar et al. [52] optimized the design of steel truss bridges.Miano et al. [53] built fragility curves of existing reinforced concrete buildings.Optimization designs of 3D reinforced concrete frames, cold-formed steel frames, frame structures, RC counterfort retaining walls, tied-arch bridges, cablestayed footbridge, framed structures, and geometrically nonlinear lattice structures can be found in [54,55], respectively.
It is evident that the MATLAB-SAP2000 interface has found widespread application in the field of structural engineering.However, its utilization in the realm of structural control is notably absent.In existing studies, the optimization of MTMD parameters relies on a simplified model of the floor structure.Nevertheless, this oversimplification can result in the diminished control effectiveness of the MTMD.
To address this gap and provide a solution, this paper introduces the development of a MATLAB-SAP2000 interface in Section 2. Additionally, a multi-objective optimization algorithm based on the AFSA [56,57] for MTMD, considering both maximum and root mean square (RMS) accelerations at selected measuring points, is proposed in Section 3. To validate the accuracy and efficacy of the developed interface and optimization algorithm, an SDOF system constructed in SAP2000 is initially tested in Section 4. Subsequently, a large-span floor structure is presented as a case study in Section 5, with TMD/MTMD designed through the conventional method included for comparison in Sections 4 and 5. Finally, conclusions are drawn in Section 6.

SAP2000-API
An excessively simplified model for the optimization of MTMD inevitably introduces errors, particularly for floor structures with closely spaced modes, thereby diminishing the control effectiveness of MTMD.To address this challenge, a MATLAB-SAP2000 interface is developed in this section.Leveraging this interface, the MTMD system can be optimized directly within the finite element model, resulting in higher precision and improved control effectiveness.
Utilizing the OAPI, users can employ various programming languages to manipulate the SAP2000 database.This significantly streamlines the modeling, analysis, optimization, extraction, and processing of SAP2000 calculation results, enabling seamless data connection and exchange with other software.The OAPI allows for two-way direct data conversion, facilitating the rapid and comprehensive import of large models.Furthermore, it permits multiple data exchange loops between different applications without the need to create a new model for each loop.Therefore, the OAPI module extends the innovative use of SAP2000 functionality, empowering developers to invoke the SAP2000 analysis and design engine through their own applications.
The main steps and codes of SAP2000-OAPI are shown in Figure 1.
Buildings 2024, 14, x FOR PEER REVIEW 3 of 24 TMD/MTMD designed through the conventional method included for comparison in Sections 4 and 5. Finally, conclusions are drawn in Section 6.

SAP2000-API
An excessively simplified model for the optimization of MTMD inevitably introduces errors, particularly for floor structures with closely spaced modes, thereby diminishing the control effectiveness of MTMD.To address this challenge, a MATLAB-SAP2000 interface is developed in this section.Leveraging this interface, the MTMD system can be optimized directly within the finite element model, resulting in higher precision and improved control effectiveness.
Utilizing the OAPI, users can employ various programming languages to manipulate the SAP2000 database.This significantly streamlines the modeling, analysis, optimization, extraction, and processing of SAP2000 calculation results, enabling seamless data connection and exchange with other software.The OAPI allows for two-way direct data conversion, facilitating the rapid and comprehensive import of large models.Furthermore, it permits multiple data exchange loops between different applications without the need to create a new model for each loop.Therefore, the OAPI module extends the innovative use of SAP2000 functionality, empowering developers to invoke the SAP2000 analysis and design engine through their own applications.
The main steps and codes of SAP2000-OAPI are shown in Figure 1.

MATLAB-SAP2000 Interface
It is mentioned in [44][45][46] that at present, the current research indicates a lack of a universally accepted design method for the optimization of MTMD applied to floor structures.It is usually simplified into an SDOF system firstly according to the floor mode of interest, and then the traditional optimization design formula [58][59][60][61][62][63] is used for MTMD optimization.However, this widely employed simplified MTMD design method exhibits three main shortcomings: firstly, the modes of the floor structure are often closely spaced, and the simplification of multiple degrees of freedom to a single degree of freedom often leads to deviations; secondly, there is a lack of a widely accepted method to design the frequency and damping ratio of each TMD in the MTMD system; and thirdly, the actual vibration modes at the measuring points on the floor (i.e., the locations of each TMD in the MTMD) differ, and the simplified model fails to capture these differences between measuring points.

MATLAB-SAP2000 Interface
It is mentioned in [44][45][46] that at present, the current research indicates a lack of a universally accepted design method for the optimization of MTMD applied to floor structures.It is usually simplified into an SDOF system firstly according to the floor mode of interest, and then the traditional optimization design formula [58][59][60][61][62][63] is used for MTMD optimization.However, this widely employed simplified MTMD design method exhibits three main shortcomings: firstly, the modes of the floor structure are often closely spaced, and the simplification of multiple degrees of freedom to a single degree of freedom often leads to deviations; secondly, there is a lack of a widely accepted method to design the frequency and damping ratio of each TMD in the MTMD system; and thirdly, the actual vibration modes at the measuring points on the floor (i.e., the locations of each TMD in the MTMD) differ, and the simplified model fails to capture these differences between measuring points.
Hence, to mitigate the errors stemming from the simplification of floor structures in existing studies, this paper aims to directly utilize the floor finite element model in SAP2000 for the optimization of MTMD parameters through MATLAB.This approach is adopted to enhance the control effectiveness of MTMD in mitigating human-induced vibrations.The MATLAB-SAP2000 interface is developed upon the SAP2000-OAPI.

Parameter Optimization Algorithm of MTMD
MATLAB excels in conducting intelligent algorithm optimization for an objective function within specified parameters and a given range.However, it faces limitations in modeling and simplifying complex structures.On the other hand, SAP2000 is proficient in constructing a 3D finite element model, but it can only analyze under predefined models and parameter conditions, making direct optimization of MTMD parameters challenging.Therefore, by integrating these two software packages, a comprehensive approach can be achieved.After structural modeling, intelligent algorithm-based MTMD parameter optimization and structure analysis can be conducted, enabling applications across a broader spectrum of engineering scenarios.
In this paper, we have chosen MATLAB version 2020 [48] and SAP2000 version 2022 [49] software to develop the "interface function", introducing a methodology that leverages the joint utilization of these two software tools.Structural modeling is conducted in SAP2000 and MTMD parameter optimization is performed in MATLAB, with seamless data transfer between the two platforms.Subsequently, MATLAB can automatically invoke the structural model in SAP2000 for MTMD parameter optimization and structural analysis.
The optimization outcomes and the effectiveness of vibration mitigation through MTMD are intricately tied to the chosen optimization target.In a bid to enhance the control effectiveness and robustness of the MTMD system for human-induced vibration control, this section introduces a multi-objective optimization algorithm.The algorithm integrates the AFSA [56,57] with the MATLAB-SAP2000 interface.The proposed method addresses the limitation of the direct application of intelligent optimization algorithms in actual engineering calculations.It significantly reduces the computational load, streamlining the MTMD parameter optimization process to be quicker, more convenient, and more accurate.

Artificial Fish Warm Algorithms
To expedite and enhance the effectiveness of MTMD parameter optimization, this section employs the AFSA for searching optimal parameters.AFSA is crafted based on the foraging, swarming, following, and random behavior of fish, presenting advantages such as strong flexibility, rapid convergence, high accuracy, and robust fault tolerance.
Four key parameters in the algorithm significantly influence convergence, namely visual range, step size, the number of artificial fish, and the number of attempts.Therefore, careful selection of these parameters is essential based on the specific situation and function.
The visual range parameter directly impacts algorithm convergence.Following the assessment of the fish swarm algorithm case and multiple attempts, this study determines a visual range of 0.5.
The step size parameter, influencing convergence speed and accuracy, requires a balance.A large value accelerates convergence but risks missing the optimal solution, leading to longer convergence times.On the other hand, a small value improves accuracy but may hinder algorithm convergence.Considering both speed and accuracy, along with the visual range, multiple trial calculations establish a step size of 1 in this paper.
The number of artificial fish directly affects the algorithm's convergence speed.A higher number accelerates convergence but increases the iterative calculation load, impacting computational efficiency.After several trial calculations, the number of artificial fish is set at 80 in this section.
Lastly, the number of attempts influences foraging ability and convergence efficiency.A higher number strengthens foraging ability but raises the risk of local extreme values.After considering the case and the factors mentioned above, the number of attempts is determined to be 50 in Section 4.2 and 300 in Section 5.3 in the subsequent 3D model analysis.

Parameter Optimization Based on the Interface and AFSA
Regarding the serviceability problem, the typical focus is on the structural acceleration response.In this context, the maximum and RMS accelerations stand out as two crucial indices.In the pursuit of an MTMD system with heightened robustness and superior control effectiveness, both of these indices are taken into consideration.The objective function Y i for the ith TMD in the system is shown in Equation (1).
where α 1 and α 2 are the two weight coefficients for the maximum and RMS accelerations, respectively, g 1 and g 2 are the mean maximum and RMS acceleration objective functions, respectively, and f ti and ξ ti represent the frequency and damping ratio of the ith TMD in the MTMD system, respectively.In Equation ( 1), under several external excitations for optimization, normalized objective functions for the ith TMD are: Then, from Equations ( 2) and ( 3), it can be found that when the objective functional value is 1.0, it means that the ith TMD can completely kill the vibration at the corresponding structural measuring point.Therefore, the ith TMD will have a better control effect when the value is larger.
In Equation ( 1), the sum of α 1 and α 2 is 1.0, which means that Y i for the ith TMD is also a normalized objective function.When α 1 is 1 and α 2 is 0, it means that only the maximum acceleration is considered, and it becomes a single-objective optimization problem.In this paper, they are both set to 0.5 as an example.
The suggested placement of each TMD in the MTMD system on the floor should align with the structural modal shape.Following the determination of the mass for each TMD in the MTMD system, the stiffness and damping coefficients can be optimized using the proposed algorithm, as depicted in Figure 2.
Through the multi-objective optimization algorithm proposed in Figure 2, the optimal frequency ratio f opti and the damping ratio ξ opti for each TMD can be searched, and an optimized MTMD system is obtained.
Hence, the AFSA-based multi-objective optimization algorithm for MTMD using the developed MATLAB-SAP2000 interface is outlined in Figure 3.The floor structure and the initial MTMD system are initially constructed in SAP2000.Subsequently, under various human-induced excitations in SAP2000 for optimization purposes, structural acceleration responses can be seamlessly transferred to MATLAB.For the ith step, swarming and following behaviors of the artificial fish are carried out to search for a larger objective functional value in Equation (1).In correspondence with this step, adjustments are made to the parameters of the MTMD in SAP2000.The adaptability evaluation aims to assess whether a globally optimal solution has been reached.If affirmative, the entire joint optimization process concludes.Additionally, the optimization process concludes when a predetermined maximum number of steps has been attempted.Through the multi-objective optimization algorithm proposed in Figure 2, the optimal frequency ratio  and the damping ratio  for each TMD can be searched, and an optimized MTMD system is obtained.
Hence, the AFSA-based multi-objective optimization algorithm for MTMD using the developed MATLAB-SAP2000 interface is outlined in Figure 3.The floor structure and the initial MTMD system are initially constructed in SAP2000.Subsequently, under various human-induced excitations in SAP2000 for optimization purposes, structural acceleration responses can be seamlessly transferred to MATLAB.For the ith step, swarming and following behaviors of the artificial fish are carried out to search for a larger objective functional value in Equation ( 1).In correspondence with this step, adjustments are made to the parameters of the MTMD in SAP2000.The adaptability evaluation aims to assess whether a globally optimal solution has been reached.If affirmative, the entire joint optimization process concludes.Additionally, the optimization process concludes when a predetermined maximum number of steps has been attempted.

Model Introduction
To validate the accuracy of the developed interface and optimization algorithm, in this section, a mass-spring-dashpot SDOF system built in SAP2000 is tested.The mass is set to 100 t, stiffness is 15,791 kN/m, indicating a natural frequency of 2.0 Hz, and a damping ratio of 1.0% is used, which is a typical model for a slender pedestrian structure.

Model Introduction
To validate the accuracy of the developed interface and optimization algorithm, in this section, a mass-spring-dashpot SDOF system built in SAP2000 is tested.The mass is set to 100 t, stiffness is 15,791 kN/m, indicating a natural frequency of 2.0 Hz, and a damping ratio of 1.0% is used, which is a typical model for a slender pedestrian structure.
The comfort design is based on [64,65].According to [64,65], for pedestrian structures, the vertical load per unit area can be expressed as: where P b is the vertical load of a single pedestrian and the value is 0.28 kN, f 1 means the frequency of the human-induced excitation, and r ′ represents the equivalent pedestrian density, and can be expressed as: where A is the area of the pedestrian structure, ξ and N are the structural damping ratio and the number of total pedestrians, respectively, and Ψ is the vertical load scale factor, which is shown in Figure 4.  To meticulously assess the performance of the pedestrian structure under walking induced excitation, the model is subjected to harmonic excitations ranging from 1.25 to 4.60 Hz, with an equal interval of 0.05 Hz, resulting in a total of 68 cases based on Equation (4).It is important to note that this section considers the scenario with only one pedestrian

Parameter Optimization of TMD
A TMD with a mass ratio of 1% is proposed as an example.The frequency searching range of TMD is set to be [0, 4 Hz], while the range of the damping ratio is set to be [0 0.1].Employing the optimization method outlined in Figure 4, the objective is to optimize the TMD to minimize both the maximum and RMS acceleration responses of the main structure under various excitations with different frequencies.The AFSA-based parame ter searching is carried out in MATLAB while the SDOF main structure with a TMD is built in SAP2000, and they are connected through the developed interface.The optimiza tion process can be seen in Figure 5.To meticulously assess the performance of the pedestrian structure under walkinginduced excitation, the model is subjected to harmonic excitations ranging from 1.25 to 4.60 Hz, with an equal interval of 0.05 Hz, resulting in a total of 68 cases based on Equation (4).It is important to note that this section considers the scenario with only one pedestrian.

Parameter Optimization of TMD
A TMD with a mass ratio of 1% is proposed as an example.The frequency searching range of TMD is set to be [0, 4 Hz], while the range of the damping ratio is set to be [0, 0.1].Employing the optimization method outlined in Figure 4, the objective is to optimize the TMD to minimize both the maximum and RMS acceleration responses of the main structure under various excitations with different frequencies.The AFSA-based parameter searching is carried out in MATLAB while the SDOF main structure with a TMD is built in SAP2000, and they are connected through the developed interface.The optimization process can be seen in Figure 5.
the TMD to minimize both the maximum and RMS acceleration responses of the main structure under various excitations with different frequencies.The AFSA-based parameter searching is carried out in MATLAB while the SDOF main structure with a TMD is built in SAP2000, and they are connected through the developed interface.The optimization process can be seen in Figure 5.In Figure 5, the proposed MATLAB-SAP2000 interface facilitates seamless data transfer between the two platforms.The optimization process for MTMD is conducted in MATLAB based on structural dynamic responses under external excitations in SAP2000.The parameters of MTMD in SAP2000 are adjusted iteratively until reaching the final optimal values.In each iterative step, the ASFA depicts the convergence trajectory of the optimization process.The blue dots represent the distribution of TMD parameters for the current iteration, while the red dots represent the optimal TMD parameters for that iteration.Initially, the distribution of fish is scattered, but as the iterations progress, the fish swarm gradually converges towards the center, particularly around the frequency band of 2 Hz.Analyzing the position of the red dot in each subgraph and its slight variations reveals that the target function in Equation ( 1) is more sensitive to the frequency ratio than the damping ratio in this case.It is evident that the target function reaches its optimal value at approximately the 10th step due to the consistent location, with a frequency of 1.98 Hz and a damping ratio of 6.06%.
From Figure 6, with only a few iterations, the optimization method proposed in this study can quickly find the optimal frequency and damping ratio of TMD.Den Hartog [58][59][60][61] provided optimal formulas for the frequency ratio and damping ratio in an undamped SDOF system, but this paper focus on damped systems.In this regard, Warburton [62,63] conducted numerical studies and provided results for reference.
From Figure 6, with only a few iterations, the optimization method proposed in this study can quickly find the optimal frequency and damping ratio of TMD.Den Hartog [58][59][60][61] provided optimal formulas for the frequency ratio and damping ratio in an undamped SDOF system, but this paper focus on damped systems.In this regard, Warburton [62,63] conducted numerical studies and provided results for reference.For a given TMD mass ratio , structural damping ratio , TMD structure frequency ratio , excitation-structure frequency ratio , or TMD damping ratio , the optimization objective is to find the minimum acceleration response  of the system under harmonic excitation.The response  is defined by Equations ( 6)- (10).
Objective Function Value For a given TMD mass ratio µ, structural damping ratio β, TMD structure frequency ratio f , excitation-structure frequency ratio r, or TMD damping ratio ξ, the optimization objective is to find the minimum acceleration response R of the system under harmonic excitation.The response R is defined by Equations ( 6)- (10).
The optimal frequency ratio f opt and damping ratio ξ opt should be chosen to minimize the peak value of R. Substituting the known variables into the equations, we can obtain the optimal f opt and damping ratio ξ opt through numerical analysis methods.
According to Equations ( 6)- (10), when β = 1%, we can obtain f opt = 0.9891 and ξ opt = 0.062, respectively.The optimized results are 1.98 Hz and 6.06%, respectively, which are nearly the same.Therefore, through a SDOF main structure, it has been verified in this section that the proposed AFSA-based multi-objective optimization algorithm for MTMD using the developed MATLAB-SAP2000 interface has the advantages of good convergence, fast convergence speed, and high precision.

Case Introduction
To further illustrate the improvement of the proposed method in the optimization design of MTMD parameters used for vibration control of floor structures, this section takes a large-span floor structure as an example.
The proposed floor structure has a plane size of 67.6 m × 43.0 m.It is a steel frame structure.The thickness of the passageway floor is 120 mm, the thickness of the roof floor is 150 mm, and the floor material is concrete C30.The structure's frame column adopts a 700 mm × 700 mm steel-reinforced concrete column, the material of the steel pipe is Q355, and material of the concrete is C50.The inclined strut mainly adopts 600 × 600 × 35 × 35 mm, 500 × 500 × 25 × 25 mm, BRB-120 × 120, BRB-180 × 180, and other box section steel, and the material is Q355.The frame beam is mainly made of 1100 × 600 × 35 × 35 mm, 900 × 600 × 30 × 30 mm, 600 × 450 × 14 × 30 mm H-beams and other sections, whose material is Q355.The total mass of the structure is 10,164 t.The finite element model of the floor structure is built in SAP2000 and shown in Figure 7.

Conventional Design Method
In the conventional design approach for MTMD in a floor structure, the typical procedure involves conducting an initial modal analysis.Subsequently, the mode that is particularly sensitive to human-induced excitation is treated as a SDOF damped system, and the TMD can be optimized using formulas by Warburton.However, when dealing with MDOF floor structures characterized by densely spaced modal features, using the abovementioned formulas cannot yield optimal results.The reason is that utilizing these formulas requires transforming the MDOF floor system into an equivalent SDOF system, which involves accurately assessing modal mass for the modes of interest.For a complex floor structure, achieving a precise equivalent SDOF system is challenging and the dynamic characteristics of the simplified SDOF model are different from those of the 3D model to

Conventional Design Method
In the conventional design approach for MTMD in a floor structure, the typical procedure involves conducting an initial modal analysis.Subsequently, the mode that is particularly sensitive to human-induced excitation is treated as a SDOF damped system, and the TMD can be optimized using formulas by Warburton.However, when dealing with MDOF floor structures characterized by densely spaced modal features, using the above-mentioned formulas cannot yield optimal results.The reason is that utilizing these formulas requires transforming the MDOF floor system into an equivalent SDOF system, which involves accurately assessing modal mass for the modes of interest.For a complex floor structure, achieving a precise equivalent SDOF system is challenging and the dynamic characteristics of the simplified SDOF model are different from those of the 3D model to some extent.This equivalence process introduces deviations, preventing the optimal vibration reduction effects when applying these formulas to 3D floor structures.
The structural mode shapes, sensitive to human-induced excitations, are depicted in Figure 8, with detailed information listed in Table 1. Figure 8 illustrates the vertical mode shape deformation of the entire floor structure.It is important to note that even though pedestrians exert force solely on the first floor, the connected inclined strut causes the second floor to influence the modal shape and frequency of the first floor.Consequently, a joint analysis of both floors is necessary.From Figure 8 and Table 1, the first natural frequency is close to the fourth natural frequency, while the fifteenth natural frequency is close to the sixteenth natural frequency, which increases the difficulty of conventional design.
Based on the mode shapes presented in Figure 8 and considering that TMDs should be implemented where the modal response is maximum, a MTMD system comprising ten TMDs is proposed in Figure 9.It is worth noting that, for the ease of optimizing MTMD parameters and facilitating production processing, these TMDs are divided into three groups.TMDs within the same group share identical parameters.
The structural mode shapes, sensitive to human-induced excitations, are depicted in Figure 8, with detailed information listed in Table 1. Figure 8 illustrates the vertical mode shape deformation of the entire floor structure.It is important to note that even though pedestrians exert force solely on the first floor, the connected inclined strut causes the second floor to influence the modal shape and frequency of the first floor.Consequently, a joint analysis of both floors is necessary.which increases the difficulty of conventional design.
Based on the mode shapes presented in Figure 8 and considering that TMDs should be implemented where the modal response is maximum, a MTMD system comprising ten TMDs is proposed in Figure 9.It is worth noting that, for the ease of optimizing MTMD parameters and facilitating production processing, these TMDs are divided into three groups.TMDs within the same group share identical parameters.By comparing Figures 8 and 9, TMD-I is used to control the first mode, but the location will be affected by the fifteenth and sixteenth modes as well.TMD-II is aimed at the fourth mode.TMD-III is mainly applied to control the fifteenth and sixteenth modes.By comparing Figures 8 and 9, TMD-I is used to control the first mode, but the location will be affected by the fifteenth and sixteenth modes as well.TMD-II is aimed at the fourth mode.TMD-III is mainly applied to control the fifteenth and sixteenth modes.
According to Table 1, the MTMD system is designed based on Equations ( 6)-( 10), while the results are shown in Table 2.Note that the sixteenth mode is merged into the fifteenth mode while TMD-III is aimed at controlling it in the conventional design.

Proposed Design Method
As discussed in Section 5.2, the conventional design method for MTMD relies on a simplified model of a floor structure with closely spaced modes, and an overly simplified floor model may result in a diminished control effect.In this section, the efficacy of the proposed AFSA-based multi-objective optimization algorithm for MTMD, utilizing the developed MATLAB-SAP2000 interface, is evaluated in the context of this floor structure.
In Equation ( 5), the area of the first floor is 1079 m 2 , as indicated in [46], and the total number of pedestrians is half of this area, amounting to 540.To assess the performance of the pedestrian structure under human-induced excitation, the model undergoes harmonic excitations ranging from 1.25 Hz to 4.60 Hz, with an incremental difference of 0.05 Hz, resulting in a total of 68 cases based on Equation ( 4).
The parameter searching using the AFSA is conducted in MATLAB.Simultaneously, the floor structure with the MTMD system is constructed in SAP2000, and the two components are interconnected through the developed interface.The optimization process is illustrated in Figure 10.The frequency searching range for each TMD is set between 1 Hz and 3.5 Hz, while the range of the damping ratio is established as [0, 0.15].
resulting in a total of 68 cases based on Equation ( 4).
The parameter searching using the AFSA is conducted in MATLAB.Simultaneously, the floor structure with the MTMD system is constructed in SAP2000, and the two components are interconnected through the developed interface.The optimization process is illustrated in Figure 10.The frequency searching range for each TMD is set between 1 Hz and 3.5 Hz, while the range of the damping ratio is established as [0, 0.15].Through the established MATLAB-SAP2000 interface, data seamlessly transfers between the two platforms.The optimization process for the MTMD is executed in MATLAB based on structural dynamic responses under external excitations in SAP2000.Parameters of the MTMD in SAP2000 are adjusted iteratively until reaching the final optimal values.As depicted in Figure 10, with the progression of iterations, the initially scattered fish swarm gradually converges into three distinct groups, corresponding to the optimal frequency and damping ratio for the three types of TMDs illustrated in Figure 9.In each iteration step, the three red dots symbolize the optimal TMD parameters for the three types of TMDs, as presented in Table 3.
Table 3. Optimal parameters of the three types of TMD in each iteration step.

TMD Type Parameters Initial Step
Step 10 Step 50 Step 100 Step 200 Step 300 Observing Table 3, it becomes apparent that with the advancement of iterations, the optimal parameters of the MTMD gradually stabilize.The results are summarized in Table 4, and the entire iterative process is depicted in Figure 11.A comparison between Table 2 and Table 4 reveals that the frequencies of the proposed optimized MTMD are consistently smaller than the values obtained through the conventional design method.Notably, there is a substantial difference in their optimal damping ratios.A comparative analysis of their human-induced vibration control effects will be conducted in the subsequent section.Through the established MATLAB-SAP2000 interface, data seamlessly transfers between the two platforms.The optimization process for the MTMD is executed in MATLAB based on structural dynamic responses under external excitations in SAP2000.Parameters of the MTMD in SAP2000 are adjusted iteratively until reaching the final optimal values.As depicted in Figure 10, with the progression of iterations, the initially scattered fish swarm gradually converges into three distinct groups, corresponding to the optimal frequency and damping ratio for the three types of TMDs illustrated in Figure 9.In each iteration step, the three red dots symbolize the optimal TMD parameters for the three types of TMDs, as presented in Table 3.
Table 3. Optimal parameters of the three types of TMD in each iteration step.

TMD Type Parameters Initial Step
Step 10 Step 50 Step 100 Step 200 Step 300 Observing Table 3, it becomes apparent that with the advancement of iterations, the optimal parameters of the MTMD gradually stabilize.The results are summarized in Table 4, and the entire iterative process is depicted in Figure 11.A comparison between Tables 2 and 4 reveals that the frequencies of the proposed optimized MTMD are consistently smaller than the values obtained through the conventional design method.Notably, there is a substantial difference in their optimal damping ratios.A comparative analysis of their human-induced vibration control effects will be conducted in the subsequent section.Analyzing Figure 11 reveals that after step 100, the value of the objective function remains constant, signifying that the iteration has converged.Hence, for a full-scale floor structure, this section provides verification that the proposed AFSA-based multi-objective optimization algorithm for MTMD, employing the developed MATLAB-SAP2000 interface, exhibits robust convergence characteristics.
For the real-scale floor structure model in Figure 9, when optimizing the MTMD system consisting of 10 TMDs, the CPU of the computer used in this paper is 3.50 GHz, the RAM is 64.0 GB, and the calculation time is about 20 h.

Comparison of the Control Effect
In this section, the human-induced vibration control effects of the MTMD from Sections 5.2 and 5.3 are compared to emphasize the enhancements brought about by the proposed method.The comparison focuses on structural acceleration spectra under various human-induced excitations with different frequencies, as illustrated by the ten measuring points in Figure 9, now presented in Figure 12.Analyzing Figure 11 reveals that after step 100, the value of the objective function remains constant, signifying that the iteration has converged.Hence, for a full-scale floor structure, this section provides verification that the proposed AFSA-based multi-objective optimization algorithm for MTMD, employing the developed MATLAB-SAP2000 interface, exhibits robust convergence characteristics.
For the real-scale floor structure model in Figure 9, when optimizing the MTMD system consisting of 10 TMDs, the CPU of the computer used in this paper is 3.50 GHz, the RAM is 64.0 GB, and the calculation time is about 20 h.

Comparison of the Control Effect
In this section, the human-induced vibration control effects of the MTMD from Sections 5.2 and 5.3 are compared to emphasize the enhancements brought about by the proposed method.The comparison focuses on structural acceleration spectra under various human-induced excitations with different frequencies, as illustrated by the ten measuring points in Figure 9, now presented in Figure 12.
As can be seen from Figure 12, although the 10 measuring points are all situated on the same span of the floor structure, their acceleration spectrum shapes and response values are different.For measuring points 1, 7, 8, and 9, there are obviously two peaks close to the same height, that is, they are affected by multiple local vertical modes at the same time.Meanwhile, the MTMD optimized in this study has a better control effect than the MTMD optimized by the conventional method, which is especially reflected in the attenuation of the first peak at around 2.0 Hz at points 2~4 and the second peak at around 2.8 Hz at point 6 and point 10.It is worth noting in Figure 12e that for point 5, in the frequency range of 2.2 Hz to 2.6 Hz, the response of the uncontrolled structure is smaller than the response after adding the MTMDs.When MTMDs are used to control vibration at a specific frequency, they can achieve better control effects just at the frequency at the peak point.However, in a frequency range near that frequency, there may be a certain degree of response amplification.This is because the added MTMDs alter the whole structure's modes, causing a shift in the structural modal frequencies.As a result, the resonance effect near the central frequency is enhanced, leading to the possibility of response amplification, or even leading to the appearance of two new peaks.These can explain why the response of the controlled structure between the frequencies at two peak points, around 2.2 Hz to 2.6 Hz, is larger than the uncontrolled structure.What is more, in Figure 12e, the peak responses of the conventional designed controlled structure and optimized designed controlled structure are both smaller than that of the uncontrolled structure.
Figure 12 shows that the optimized parameters outperform the conventional method.Since the optimized method employs the 3D model to run an iterative computation, the results can be regarded as optimal.The conventional method, although using the optimal formula, is designed for a SDOF system.This requires us to simplify the MDOF system and find an equivalent SDOF system instead.The equivalent simplification introduces computational deviations, preventing the parameters in this formula from achieving optimal solutions.Additionally, the vibration reduction effects of MTMD in the actual MDOF model interact with each other, and this cannot be reflected in these formulas.
In Figure 12, the maximum accelerations of the case without MTMD, with the conventional designed MTMD, and with the proposed MTMD of each subgraph are compared in Table 5, while reductions of the two MTMD systems are presented in Table 6.Analyzing Figure 11 reveals that after step 100, the value of the objective function remains constant, signifying that the iteration has converged.Hence, for a full-scale floor structure, this section provides verification that the proposed AFSA-based multi-objective optimization algorithm for MTMD, employing the developed MATLAB-SAP2000 interface, exhibits robust convergence characteristics.
For the real-scale floor structure model in Figure 9, when optimizing the MTMD system consisting of 10 TMDs, the CPU of the computer used in this paper is 3.50 GHz, the RAM is 64.0 GB, and the calculation time is about 20 h.

Comparison of the Control Effect
In this section, the human-induced vibration control effects of the MTMD from Sections 5.2 and 5.3 are compared to emphasize the enhancements brought about by the proposed method.The comparison focuses on structural acceleration spectra under various human-induced excitations with different frequencies, as illustrated by the ten measuring points in Figure 9, now presented in Figure 12.As can be seen from Figure 12, although the 10 measuring points are all situated on the same span of the floor structure, their acceleration spectrum shapes and response values are different.For measuring points 1, 7, 8, and 9, there are obviously two peaks close to the same height, that is, they are affected by multiple local vertical modes at the same time.Meanwhile, the MTMD optimized in this study has a better control effect than the MTMD optimized by the conventional method, which is especially reflected in the attenuation of the first peak at around 2.0 Hz at points 2~4 and the second peak at around 2.8 Hz at point 6 and point 10.It is worth noting in Figure 12e that for point 5, in the frequency range of 2.2 Hz to 2.6 Hz, the response of the uncontrolled structure is smaller than the response after adding the MTMDs.When MTMDs are used to control vibration at a specific frequency, they can achieve better control effects just at the frequency at the peak point.However, in a frequency range near that frequency, there may be a certain degree of response amplification.This is because the added MTMDs alter the whole structureʹs modes, causing a shift in the structural modal frequencies.As a result, the resonance effect near the central frequency is enhanced, leading to the possibility of response amplification, or even leading to the appearance of two new peaks.These can explain why the response of the controlled structure between the frequencies at two peak points, around 2.2 Hz to 2.6 Hz, is larger than the uncontrolled structure.What is more, in Figure 12e, the peak responses of the conventional designed controlled structure and optimized designed controlled structure are both smaller than that of the uncontrolled structure.
Figure 12 shows that the optimized parameters outperform the conventional method.Since the optimized method employs the 3D model to run an iterative computation, the results can be regarded as optimal.The conventional method, although using the optimal formula, is designed for a SDOF system.This requires us to simplify the MDOF system and find an equivalent SDOF system instead.The equivalent simplification introduces computational deviations, preventing the parameters in this formula from achieving optimal solutions.Additionally, the vibration reduction effects of MTMD in the actual MDOF model interact with each other, and this cannot be reflected in these formulas.
In Figure 12, the maximum accelerations of the case without MTMD, with the conventional designed MTMD, and with the proposed MTMD of each subgraph are compared in Table 5, while reductions of the two MTMD systems are presented in Table 6.It is evident from Tables 5 and 6 that the proposed optimization method significantly enhances the control effectiveness of the MTMD compared to the conventional design method.Notably, the maximum reduction achieved by the conventionally designed MTMD is 50.10%, while the proposed MTMD achieves a remarkable maximum reduction of 65.51%.Specifically, for measuring points 7-9, the reduction provided by the proposed MTMD is more than twice that achieved by the conventional design method.
When the structure resonates, its vibration response is maximized, and the vibration reduction of the TMD with optimized parameters is at its maximum.When the external excitation frequency deviates from the resonance interval, the response of the structure is significantly reduced, leading to a reduction in the control effect of the TMD.Upon comparing Tables 5 and 6, it becomes apparent that for point 2, point 3, and point 4, where the acceleration response is most pronounced, the control effect of the MTMD is more significant and the optimized MTMD in this paper demonstrates superior control effectiveness.For other measurement points where the acceleration response of the structure is smaller, and the vibration mitigation effect of the MTMD is less pronounced, the MTMD system proposed in this paper still exhibits better vibration reduction performance than the MTMD system obtained through the conventional design method.
Finally, acceleration time history is compared in Figure 13.As for points 1-5 and 7-9, the excitation frequency is 2.10 Hz, while it is 2.80 Hz for points 6 and 10.
In Figure 13, it is obvious that the proposed MTMD system has a better human-induced vibration control effect at points 2-4 and points 7-9.
In Section 4, for a SDOF damped structure, the proposed method obtains an optimized TMD almost the same as the TMD from the work of G.B. Warburton [62,63].However, in this section, for a full-scale floor structure, the proposed method obtains an optimized MTMD system which is more effective than the conventional designed MTMD system based on the Warburton formula.This indicates that the Warburton formulas have certain limitations when applied to MDOF systems.One of the reasons for the deviations is the process of transforming a MDOF system into a SDOF equivalent.Additionally, these formulas cannot consider the mutual influence between MTMDs at different locations.In practice, to obtain optimal parameters for MTMD in a MDOF system, the use of interfacebased algorithm optimization is a good example.In Figure 13, it is obvious that the proposed MTMD system has a better human-induced vibration control effect at points 2-4 and points 7-9.
In Section 4, for a SDOF damped structure, the proposed method obtains an optimized TMD almost the same as the TMD from the work of G.B. Warburton [62,63].However, in this section, for a full-scale floor structure, the proposed method obtains an optimized MTMD system which is more effective than the conventional designed MTMD system based on the Warburton formula.This indicates that the Warburton formulas have certain limitations when applied to MDOF systems.One of the reasons for the deviations is the process of transforming a MDOF system into a SDOF equivalent.Additionally, these formulas cannot consider the mutual influence between MTMDs at different locations.In practice, to obtain optimal parameters for MTMD in a MDOF system, the use of interfacebased algorithm optimization is a good example.
In this study, the commonly used industrial software, SAP2000, was utilized as an interface for optimization design.This did not significantly increase time or manpower costs, yet it facilitated the optimal design of MTMD.Further, the proposed AFSA-based multi-objective optimization algorithm for MTMD can improve its control effect as well.

Conclusions
To enhance the serviceability of a floor structure, the application of a MTMD system is widespread.The parameters of the MTMD should be meticulously designed to achieve satisfactory performance.However, the prevalent use of a simplified model for the floor In this study, the commonly used industrial software, SAP2000, was utilized as an interface for optimization design.This did not significantly increase time or manpower costs, yet it facilitated the optimal design of MTMD.Further, the proposed AFSA-based multi-objective optimization algorithm for MTMD can improve its control effect as well.

Conclusions
To enhance the serviceability of a floor structure, the application of a MTMD system is widespread.The parameters of the MTMD should be meticulously designed to achieve satisfactory performance.However, the prevalent use of a simplified model for the floor structure with closely spaced modes in optimizing MTMD parameters may result in a diminished control effect.
To address this challenge, the present study utilizes the OAPI feature of SAP2000 to establish a connection with MATLAB.A multi-objective optimization algorithm, based on the AFSA, is developed in MATLAB, while the finite element model of an actual floor structure is constructed in SAP2000.Initial positions for the MTMD system are designated in SAP2000.Subsequently, through the devised MATLAB-SAP2000 interface, data can be seamlessly exchanged between the two platforms.The optimization process for the MTMD is conducted in MATLAB based on structural dynamic responses obtained from SAP2000, and the parameters of the MTMD in SAP2000 are finely tuned until reaching the ultimate optimal values.Through two numerical simulations, the following conclusions can be drawn: (1) Through a SDOF main structure, it has been verified in this section that the proposed AFSA-based multi-objective optimization algorithm for TMD using the developed
are the maximum values of the maximum structural accelerations of different excitations at the measuring point where the ith TMD is located, without and with the ith TMD, respectively, while .. y rms and .. y tirms represent the maximum values of the RMS responses of different excitations, respectively.It should be noted that .. y max and .. y timax are probably not from the same external excitation case, and it is likely that the excitation frequencies of these two cases are different... y rms and .. y max are not necessarily from the same excitation case either.

Figure 3 .
Figure 3. Flow chart of the MATLAB-SAP2000 interface for MTMD optimization.

Figure 6 .
Figure 6.Objective function value with the process of iteration.

Figure 6 .
Figure 6.Objective function value with the process of iteration.

Figure 7 .
Figure 7.The 3D view of the floor structure in SAP2000.

Figure 7 .
Figure 7.The 3D view of the floor structure in SAP2000.

Figure 8 .
Figure 8. Mode shapes: (a,b) Structural and the first floor s first mode; (c,d) Structural and the first floor s fourth mode; (e,f) Structural and the first floor s fifteenth mode; (g,h) Structural and the first floor s sixteenth mode.

Figure 8 .
Figure 8. Mode shapes: (a,b) Structural and the first floor's first mode; (c,d) Structural and the first floor's fourth mode; (e,f) Structural and the first floor's fifteenth mode; (g,h) Structural and the first floor's sixteenth mode.

Figure 9 .
Figure 9.The layout of the MTMD system on the first floor.

Figure 9 .
Figure 9.The layout of the MTMD system on the first floor.

Figure 11 .
Figure 11.Objective function value with the process of iteration of the MTMD system.

Figure 11 .
Figure 11.Objective function value with the process of iteration of the MTMD system.

Figure 11 .
Figure 11.Objective function value with the process of iteration of the MTMD system.

Table 2 .
Conventional design results of the MTMD system.

Table 4 .
Optimized parameters of the MTMD system by the proposed method.

Table 4 .
Optimized parameters of the MTMD system by the proposed method.

Table 5 .
Comparison of the maximum accelerations at the measuring points.

Table 6 .
Comparison of the reductions of different MTMDs at the measuring points.