Passive Control of Ultra-Span Twin-Box Girder Suspension Bridges under Vortex-Induced Vibration Using Tuned Mass Dampers: A Sensitivity Analysis

: Suspension bridges’ in-plane extended conﬁguration makes them vulnerable to wind-induced vibrations. Vortex shedding is a kind of aerodynamic phenomenon causing a bridge to vibrate in vertical and torsional modes. Vortex-induced vibrations disturb the bridge’s serviceability limit, which is not favorable, and in the long run, they can cause fatigue damage. In this condition, vibration control strategies seem to be essential. In this paper, the performance of a tuned mass damper (TMD) is investigated under the torsional vortex phenomenon for an ultra-span streamlined twin-box girder suspension bridge. In this regard, the sensitivity of TMD parameters was addressed according to the torsional responses of the suspension bridge, and the reached appropriate ranges are compared with the outputs provided by genetic algorithm. The results indicated that the installation of three TMDs could control all the vulnerable modes and reduce the torsional rotation by up to 34%.


Introduction
Within the variety of bridge designs, suspension bridges are known for their capacity to span vast distances, often being the top choice for connecting distant points through extended structural spans.Consequently, the significance of wind-driven oscillations for these bridges increases as their span lengthens.These bridges can vibrate independently in four primary modes-vertical, longitudinal, lateral, and torsional-or as a blend of these modes [1][2][3].Aerodynamic phenomena, such as vortex-shedding, galloping, buffeting, static divergence, and flutter, are caused by the interaction between airflow and a broad structure like a bridge deck.Vortex shedding is the primary cause of structural vibrations in the direction perpendicular to the airflow, although other velocity-dependent forces become relevant when the induced motion is substantial [4].
Recently, vortex-induced vibration (VIV) has been detected on numerous bridges globally, including Japan's Trans-Tokyo Bay Bridge [5], Russia's Volgograd Bridge [6], China's Xihoumen Bridge, Yingwuzhou Yangtze River Bridge, and Humen Bridge [7], as well as Denmark's Great Belt Bridge.It is essential to acknowledge that vortex sheddinginduced oscillations are predominantly resonant, self-regulating, narrowly focused, and usually occur at low wind speeds, impacting driving safety, contributing to long-term fatigue damage, and reducing the average bridge life span [8][9][10].
If it is determined that a structure does not meet aerodynamic requirements, there are three strategies for bridge control under VIV as follows [11][12][13]:

•
Prevent vortex shedding by refining the bridge deck's aerodynamic design, utilizing distinct streamlined box girders.

•
Modify the bridge's dynamic properties by changing its structural layout.

•
Enhance stability by incorporating additional attachments such as spoilers, guide vanes, or flapping plates.
The tuned mass damper (TMD) system, a passive mechanical control approach, falls under the third strategy and is frequently employed in suspension bridges to manage vibrations resulting from various load types.Three primary parameters make up a TMD, specifically, mass ratio, tuning frequency, and damping ratio [14,15].Originally conceived by Frahm [16], the TMD was a spring-mass system that could only dampen vibrations of a single frequency because of its inability to retain excess energy.Den Hartog [17] later proposed a set of optimal TMD design formulas based on fixed-point theory and included a viscous damper into the system, making it the most commonly used setup.
A novel passive system was tested in a wind tunnel and shown to be effective by Kwon et al. [18], who used a TMD mechanism to activate a plate and adjust the airflow on the deck.Gu et al. [19] proposed a new lever-type TMD to counteract wind-induced vibrations, testing their suggested control system on the Yichang suspension bridge and demonstrating greater efficiency than a passive TMD.Pourzeynali and Datta [20] examined the effects of TMD parameters on accelerating the rate of flutter, finding that TMD did not introduce any instability.Chen and Kareem [21] studied TMD's effectiveness in controlling a bridge's flutter, including an optimal TMD design strategy based on the negative damping of systems.
A unique control strategy for dampening modal coupling effects and reducing resonant vibrations via TMDs was proposed by Chen and Cai [22].The effectiveness of TMD in increasing the gallop speed of flexible suspension bridges without modifying their forms was studied by Abdel-Rohman and John [23].Domaneschi et al. [24] studied how TMD can be used to manage buffeting on suspension bridges.For the Vincent Thomas suspension bridge, Alizadeh et al. [25] conducted a sensitivity analysis of flutter velocity concerning the gyration radius and placement of TMDs.Kontoni and Farghaly [26] used tuned mass dampers to mitigate soil-structure interaction effects on a cable-stayed bridge's seismic response.Finally, Mansouri et al. [27] investigated how far-fault earthquake duration, severity, and size affected the seismic response of RC bridges retrofitted with seismic bearings.
Patil and Jangid [28] and Bandivadekar and Jangid [29] reported that the performance of TMD depends on tuning frequency and optimum damping ratio.Researchers have made numerous efforts to address the mistuning issue of single TMD (STMD), such as implementing nonlinear TMD [30], adaptive-passive TMD [31], and semi-active TMD [32].Additionally, a new idea involving multiple-tuned mass dampers (MTMD), which consists of an array of sub-TMDs with differing frequencies, was put forth by researchers [33,34].
In all types of vibrations, the mass ratio was introduced as the most important factor, causing the large static displacement in the single and multi-TMD configuration [35].While in torsional motion, the distribution of mass block should be accounted for instead of its weight, increasing the mass does not necessarily improve the efficacy.In fact, the optimum distribution of mass should simultaneously provide the best performance and lightest mass.
In this paper, the effects of TMD parameters on response reduction under torsional VIV are addressed by sensitivity analysis.In this regard, an ultra-span streamlined twin-box girder suspension bridge is chosen for the case study.The VIV analysis will be done on certain vulnerable modes.Then, TMDs are placed along the span according to the shape of the considered modes.The variation of four main parameters called mass ratio, distance factor, damping ratio, and frequency ratio are investigated in the reduction in VIV.The appropriate range is provided by interpreting the related curvatures, and then the results are compared with the outputs of the genetic algorithm.Finally, the most important points will be represented.

Aerodynamic Phenomena
Airflow, approaching the deck, causes two independent static and dynamic responses.The evaluation of static response is simple, while conditions are too sophisticated in computing the dynamic responses.Totally, the following three main reasons make dynamic responses: 1.
Oncoming airflow inherently includes turbulence, i.e., it is fluctuating in time and space.

2.
The flow separates at the body's sharp edges because of the friction that causes further turbulence and vortex development on the surface.As a result, there is vortex shedding in the body's wake since the flow past the body is unstable, with a variable fraction alternating from side to side.

3.
The flexible body will be vibrated by fluctuating forces.In this condition, the interaction between the oscillating body and alternation flow will cause additional forces.
In wind engineering, the first item is named buffeting, the second one is known as vortex shedding, and finally, the third part expresses the motion inducted or self-excited forces.Commonly, VIV will have occurred in low velocities, the buffeting will be seen in higher ones, and the motion-inducted forces affect the structure in too-high velocities.It is worth noting that this division is only for convenience.There are some other phenomena disturbing the stability of the bridge, such as static divergence, galloping, motion stability limit in torsion, and flutter [36].As per mentioned contents, the aerodynamic phenomena can be categorized according to the instability and serviceability limits, as shown in Figure 1.

Vortex Shedding Vibration
Vortex shedding causes bridges to vibrate in vertical and torsion degrees of freedom.VIV will be remarkable around the natural frequencies of bridge modes; it can be considered as a narrow-band frequency process.In fact, when the corresponding frequency of vortices is equal (or nearly equal) to the bridge modes frequency, the vibrations will be seen in larger amplitude.Thus, there is a certain velocity for each mode that resonance occurs in that.
where f i , D, and St are the mode frequency, depth of deck, and Strouhal number, respectively.
As VIV is the narrow band process, it is rational to evaluate the bridge response in certain modes.Thus, lower modes should be more accurately studied.As mentioned, the total response is divided into the static and dynamic parts, such that the mean value of the dynamic response is zero.In this regard, the auto spectral density of response (S r (ω)) and its variance (σ 2 r ) can be written as follows: where z and θ signify the vertical and torsional degrees of freedom, respectively.In addition, ω, φ n and x r denote the frequency, mode shape, and position in which the computation is done in there, respectively.The normalized frequency response function ( Ĥη n (ω)) and normalized auto spectral density of wind load (S Qn (ω)) are as follows: where L exp and λ are the length of the bridge experiencing the vortex inducted loads and the non-dimensional coherence length scale of vortices, respectively.In addition, M n is modal mass.The aerodynamic damping (ζ ae ) can be evaluated by the following relation: where m n is the modally equivalent and evenly distributed mass: In addition, the auto spectral density of the load may be evaluated by the following relation: where ρ, V, B, and b n are the air density, wind velocity, width of the deck, and bandwidth parameters.In addition, the standard deviation of load can be written as follows:

The Considered Bridge
The Halsafjorden suspension bridge, planned by the Norwegian Public Roads Administration (NPRA), aims to connect the west side of the fjord to Akvik on the east.Concrete towers, standing 265 m tall, support the main span on both sides, with cables extending down to rock anchor blocks across a 410-m side span.The cables have a sag of 205 m, resulting in a sag ratio of 0.1.The main cables possess a tensile strength of 1860 MPa, and each has a constructive steel area of A = 375 cm 2 .Hangers are spaced 30 m apart.
The deck features a twin box girder design, selected for its advantageous aerodynamic characteristics.Each box girder has a height of 2.5 m, a width of 11 m, and an area of 0.4430 m 2 , while the gap between them measures 10 m.Transverse stiffened steel girders, with a height of 2.5 m, a width of 1.5 m, and a constructive steel area of 0.132 m 2 , connect the two box girders.These transverse girders align with the hangers and are positioned at 30-m intervals along the bridge span.The bridge's configuration and cross-section can be seen in Figures 2 and 3, while Table 1 presents the bridge's structural and geometrical property specifications.

9825
To assess the bridge's dynamic properties, the FEM is employed to calculate the structural properties matrices.The bridge is divided into 100 specific finite elements for this purpose.Each finite element consists of two nodes at its endpoints, with multi degrees of freedom including bending rotation, torsional rotation, vertical displacement, and warping as illustrated in Figure 4. Lavasani et al. [14] provide a detailed explanation of the computation process, and the dynamic specifications of the bridge are presented in Table 2.In Table 2, TS and TA denote the torsional symmetric and torsional antisymmetric mode shape of the bridge, respectively.To investigate the VIV, the flutter derivatives of the bridge are needed, as depicted in Figure 5.

TMD Device
Tuned mass dampers have been widely employed in numerous studies for vibration control [26,[36][37][38][39][40][41][42][43], with their straightforward design process and low maintenance costs being their key attributes.On the other hand, their downsides include high sensitivity to mistuning, significant space requirements for installation, and an inability to adapt their parameters to dynamic applied loads.Nonetheless, TMDs are often the first choice for mitigating vibrations in structures, both in research and practical applications.Typically, TMDs are configured to target a specific mode (usually the lower one) to dissipate incoming dynamic energy through their dampers [37].Adjusting parameters such as mass ratio, damping ratio, and tuning frequency is crucial during the design process.Consequently, these parameters must be set to optimal values, which can be determined through an optimization process.
A too-heavy mass ratio causes considerable static deformations and changes the natural dynamic specifications.On the other hand, a light mass block does not have enough efficacy.The damping ratio should be adjusted to values, allowing free movement of the mass block.The tuning frequency may be named as the most sensitive parameter of it.As a simple way to overcome this problem, instead of tuning the TMD to a certain mode, it is rational to tune it to a ratio of a certain value that may involve the other vulnerable close modes.In this regard, the prone modes should be determined, and then an optimal ratio can be provided.The demonstration of attachment of TMD to the considered bridge is demonstrated in Figure 6.According to Figure 6, two masses are attached to the bottom of the section on the outside.

Numerical Analysis
In this paper, the torsional VIV of the considered bridge is controlled by TMD.In this regard, some parameters are used for evaluating the wind-corresponding functions, provided in Table 3.To reach a comparative view, the uncontrolled responses should be evaluated.Regarding the narrow band frequency process, the resonance velocity of each mode and its corresponding response are required.In terms of the considered bridge, computations indicate that the first six modes expose high amplitude vibrations compared to other ones, shown in Figures 7-9.In fact, the higher modes are not vulnerable to VIV. ones, shown in Figures 7-9.In fact, the higher modes are not vulnerable to VIV.     Figure 9 shows that the first and second modes experience higher amplitude vibrations, disturbing the serviceability limit.Thus, the mode shapes of these modes determine the placement of TMDs.According to Figure 7, L 4 , L 2 and 3L 4 are the most important points.Mentioned points also undergo almost high amplitudes in the other modes.Thus, the analysis is conducted according to two following conditions:   Figure 9 shows that the first and second modes experience higher amplitude vibrations, disturbing the serviceability limit.Thus, the mode shapes of these modes determine the placement of TMDs.According to Figure 7, Mentioned points also undergo almost high amplitudes in the other modes.Thus, the analysis is conducted according to two following conditions: In each case, its parameters behavior is addressed by sensitivity analysis, and the related curvatures will be interpreted.

Mass Polar Moment of Inertia
The mass polar moment of inertia involves two primary parameters: mass ratio and the distance between the mass and the rotation axis, which determine the distribution of mass around the rotation axis.In this context, this distance is referred to as the distance factor (r).The optimization process aims to achieve the lightest mass with an appropriate distribution around the rotation axis.For this purpose, r is assumed to be 1, and the sensitivity of the mass ratio is examined based on the standard deviation of the midpoint (L/2) and the point with the maximum torsional response across the entire span.It is important to note that adding a TMD may alter the location of the maximum response.Figure 10 depicts the relevant curvatures for the mass ratio.In each case, its parameters behavior is addressed by sensitivity analysis, and the related curvatures will be interpreted.

Mass Polar Moment of Inertia
The mass polar moment of inertia involves two primary parameters: mass ratio and the distance between the mass and the rotation axis, which determine the distribution of mass around the rotation axis.In this context, this distance is referred to as the distance factor (r).The optimization process aims to achieve the lightest mass with an appropriate distribution around the rotation axis.For this purpose, r is assumed to be 1, and the sensitivity of the mass ratio is examined based on the standard deviation of the midpoint (L/2) and the point with the maximum torsional response across the entire span.It is important to note that adding a TMD may alter the location of the maximum response.According to Figure 10, the maximum corresponding response of uncontrolled VIV in the middle point and whole span is the same as and equal to 6.72 degrees.In Figure 10a, by increasing the mass ratio, the middle point responses of the second mode were According to Figure 10, the maximum corresponding response of uncontrolled VIV in the middle point and whole span is the same as and equal to 6.72 degrees.In Figure 10a, by increasing the mass ratio, the middle point responses of the second mode were decreased by an appropriate slope, while the other modes' corresponding responses have stayed unchangeable.In Figure 10b, the maximum response of the whole bridge in the second mode has been decreased by an increment of mass ratio.Comparing Figure 10a,b indicates that the local and total performance of TMD is 25.5 (6.72 was decreased to 5 degrees) and 18 (6.72 was decreased to 5.5 degrees) percent, respectively.Hence, the local performance of TMD is more suitable than its total performance.In addition, another distinguishing feature of the Figure is that although TMD cannot reduce VIV-related responses of other modes, it does not increase them.Now, the computations are carried on by investigating the effects of the distance factor.In this regard, Figure 11 demonstrates the related curvatures of it.In the horizontal axis of Figure 11,

𝑟 𝐵
, signifying distance factor to deck width ratio, is utilized rather than .This ratio may provide more practical considerations.Figure 11a shows the maximum response of the middle point in the second mode, 6 degrees, is decreasing by an increment of   .The gradient of the second mode curvatures has been decreased by an increment of   , revealing there is a certain value that after it the increment of   is not too effective.Here, the mentioned value is 0.25.Mode three curvature is similar to mode five, and a little ascending has been seen in the initial part of it.Figure 11b shows that the appropriate range of   is 0.1 up to 0.2.Like the mass ratio parameter, the local performance of TMD (6 → 0.4 ⇒ 93%) is too better than the total performance of it (6 → 5 ⇒ 16%).
Comparing  In the horizontal axis of Figure 11, r B , signifying distance factor to deck width ratio, is utilized rather than r.This ratio may provide more practical considerations.Figure 11a shows the maximum response of the middle point in the second mode, 6 degrees, is decreasing by an increment of r B .The gradient of the second mode curvatures has been decreased by an increment of r B , revealing there is a certain value that after it the increment of r B is not too effective.Here, the mentioned value is 0.25.Mode three curvature is similar to mode five, and a little ascending has been seen in the initial part of it.Figure 11b shows that the appropriate range of r B is 0.1 up to 0.2.Like the mass ratio parameter, the local performance of TMD ( 6 → 0.4 ⇒ 93% ) is too better than the total performance of it (6 → 5 ⇒ 16% ).
Comparing  Figure 12a displays the main role of mass ratio in response reduction of first and second modes corresponding responses.Figure 12b shows that by increasing the mass ratio, the second-mode responses are intensified, while first-mode ones are decreased.Hence, two TMDs, located at the L 4 and 3L 4 points, provide positive local performance in all considered modes; however, their total one is not positive for all modes.Figure 13a indicates that the increment of the distance factor considerably improves the local performance of TMD for all modes.Figure 13b shows that increasing the distance factor first reduces the responses; however, after   = 0.2, the responses of other modes are intensified.

Damping Ratio
In order to investigate the sensitivity of VIV against the damping ratio, it changes in the 0 up to 20% range, demonstrated in Figure 14.In addition, the mass ratio and distance factor are considered 1% and 0.2, respectively, for both one and two TMD conditions.Figure 13a indicates that the increment of the distance factor considerably improves the local performance of TMD for all modes.Figure 13b shows that increasing the distance factor first reduces the responses; however, after r B = 0.2, the responses of other modes are intensified.

Damping Ratio
In order to investigate the sensitivity of VIV against the damping ratio, it changes in the 0 up to 20% range, demonstrated in Figure 14.In addition, the mass ratio and distance factor are considered 1% and 0.2, respectively, for both one and two TMD conditions.are intensified.

Damping Ratio
In order to investigate the sensitivity of VIV against the damping ratio, it changes in the 0 up to 20% range, demonstrated in Figure 14.In addition, the mass ratio and distance factor are considered 1% and 0.2, respectively, for both one and two TMD conditions.14a discloses that by increasing the damping ratio the L 2 point corresponding responses in modes two, three, and five are considerably reduced.Of course, the second mode response reduction gradient is too high compared to other ones in the 0 up to 10 percent range.Figure 14b indicates that the effect of the damping ratio on the maximum response of the whole bridge is too low compared to any certain point like L 2 .Especially, for mode two, the local performance of TMD is 90% ( 4.3 → 0.4) , while the total one is 7.6% ( 5.2 → 4.8 ).Thus, the damping ratio parameter is effective in the local performance of TMD, and its total efficacy is nuance.
Figure 15 demonstrates the sensitivity of responses in relation to the variation of damping ratio for the two TMDs condition.
Like the one TMD condition, the damping ratio renders suitable local performance, and its appropriate range is 0 up to 10 percent.point corresponding responses in modes two, three, and five are considerably reduced.Of course, the second mode response reduction gradient is too high compared to other ones in the 0 up to 10 percent range.Figure 14b indicates that the effect of the damping ratio on the maximum response of the whole bridge is too low compared to any certain point like  2 .Especially, for mode two, the local performance of TMD is 90% (4.3 → 0.4), while the total one is 7.6% (5.2 → 4.8 ).Thus, the damping ratio parameter is effective in the local performance of TMD, and its total efficacy is nuance.
Figure 15 demonstrates the sensitivity of responses in relation to the variation of damping ratio for the two TMDs condition.
Like the one TMD condition, the damping ratio renders suitable local performance, and its appropriate range is 0 up to 10 percent.

Frequency Ratio
According to mentioned contents of the former sections, the mistuning of TMD can considerably decline the efficacy of it.In this regard, the responses should be addressed in a range around a certain mode frequency.Figure 16 shows the VIV corresponding response in the 0.7 up to 1.4 variations of the frequency ratio.The mass ratio, the r/B, and the damping ratio are 1%, 0.2, and 10%, respectively.
Figure 16a reveals that the corresponding responses of

Frequency Ratio
According to mentioned contents of the former sections, the mistuning of TMD can considerably decline the efficacy of it.In this regard, the responses should be addressed in a range around a certain mode frequency.Figure 16 shows the VIV corresponding response in the 0.7 up to 1.4 variations of the frequency ratio.The mass ratio, the r/B, and the damping ratio are 1%, 0.2, and 10%, respectively.
Figure 16a reveals that the corresponding responses of L 2 point are declined by the increment of the frequency ratio.The final part of the second mode's corresponding curvature indicates that when the tuning frequency deviates considerably from the frequency of the modes, unfavorable effects will emerge.Figure 16a,b display that after f = 1.3ω 2 , the maximum response of the whole span and L 2 point will be increased.A similar statement of the one TMD condition can be expressed for the two TMDs.It is worth noting that one TMD is tuned to the second mode, and two TMDs are tuned to the first mode.

Combined Condition
There, TMDs whose appropriate parameters were achieved using genetic algorithm and sensitivity analysis were attached to  4 and 5.A similar statement of the one TMD condition can be expressed for the two TMDs.It is worth noting that one TMD is tuned to the second mode, and two TMDs are tuned to the first mode.

Combined Condition
There, TMDs whose appropriate parameters were achieved using genetic algorithm and sensitivity analysis were attached to  4 and 5.A similar statement of the one TMD condition can be expressed for the two TMDs.It is worth noting that one TMD is tuned to the second mode, and two TMDs are tuned to the first mode.

Combined Condition
There, TMDs whose appropriate parameters were achieved using genetic algorithm and sensitivity analysis were attached to L 4 , L 2 and 3L 4 points to reduce the VIV responses are provided in Tables 4 and 5.The middle point TMD parameters are adjusted according to the second mode, and two other ones obey the first mode specifications.Figures 18 and 19 show the standard deviation of the critical point of each mode and the total response of the bridge in uncontrolled and controlled conditions, respectively.Figure 18 indicates that TMD can successfully reduce the standard deviation of critical points.There is a little increment in mode four, which is negligible.Thus, the local performance of TMD is better in lower modes and has not had a negative role in the other modes.In addition, the parameters achieved by sensitivity analysis provide better performance than genetic algorithm ones.
For addressing the performance of TMD in response to reduction during the time, the following relation can be utilized: x(t) = N ∑ j=1 c j cos ω j t + ψ j (9) c j = 2S x ω j ∆ω j 1/2 (10) where Ψ j is arbitrary phase angles between zero and 2π, ∆ω is the frequency segment, and ω is the middle point of each segment.

Conclusions
In this paper, the torsional VIV of an ultra-span twin box girder suspension bridge, with a 2050 m span length, was controlled by TMDs.Firstly, VIV analysis of the bridge was done to determine the vulnerable modes that were the six initial ones.The parameters of TMD named mass ratio, distance factor, damping ratio, and frequency ratio were investigated under VIV to achieve their optimum ranges.A sensitivity analysis was conducted according to each mode critical point and the maximum response along the span.After computing the appropriate parameters by sensitivity analysis and genetic algorithm, the torsional VIV was controlled and comparisons were made.The most important results are as follows: The mass polar moment of inertia should be optimized in two independent phases, including mass ratio and distribution of mass around the torsion axis.Increment of mass ratio and r B improve the local and total performance of TMD, and increment of r B has a stronger impression on the local performance in comparison to the total one.The r B has the most integral role in reducing the vulnerable modes in comparison to the mass ratio.The appropriate range of mass ratio and r B are 0.5 up to 1.5 percent and 0.15 up to 0.2, respectively.
For improving the local performance, by moving away from the middle point, the mass blocks become lighter, and their distance increase.The TMD, located at the middle point, has the most weight, and the distance of its mass block is the minimum.The better total performance needs heavier mass blocks placed closer to each other.
The damping ratio has perfectly local performance, and its total efficacy is nuance.Its optimum range is 0 up to 10 percent.
The increment of frequency ratio increases the local and total efficacy of TMD such that the local one is more sensible.The optimum range of it is 1 up to 1.3 times to the considered mode.A higher value of 1.3 causes unfavorable results.
For the considered bridge, attaching three TMDs to L 4 , L 2 and 3L 4 points, comprising a 1% mass ratio, 0.15 r B , 10% damping ratio, and 1.25 frequency ratio, decreased the maximum response, 6.72 • , to 4.42 • , indicating a 34% response reduction.The considered TMDs suitably mitigated the vibration amplitude during the resonance velocity of vulnerable modes preventing the occurrence of fatigue damage in the long term.

Figure 2 .
Figure 2. In-plane view of the considered bridge (Dimensions in m).

Figure 3 .
Figure 3. Cross-section of the considered bridge, (a) Geometry of the deck (Dimensions in m), and (b) Configuration of the deck.

Figure 4 .
Figure 4. Finite element of the suspension bridge model.

Figure 4 .
Figure 4. Finite element of the suspension bridge model.

Figure 5 .
Figure 5. Flutter derivatives of the considered bridge.

Figure 5 .
Figure 5. Flutter derivatives of the considered bridge.

Figure 4 .
Figure 4. Finite element of the suspension bridge model.

Figure 5 .
Figure 5. Flutter derivatives of the considered bridge.

Figure 7 .
Figure 7.The standard deviation of six initial modes.Figure 7. The standard deviation of six initial modes.

Figure 7 .
Figure 7.The standard deviation of six initial modes.Figure 7. The standard deviation of six initial modes.

Figure 8 .
Figure 8.Total response of the bridge.Figure 8. Total response of the bridge.

Figure 8 .
Figure 8.Total response of the bridge.Figure 8. Total response of the bridge.
(a) One TMD is placed at the middle point.(b) Two TMDs are placed at the L 4 and 3L 4 points.

Figure 8 .
Figure 8.Total response of the bridge.

Figure 9 .
Figure 9. Order and magnitude of six initial modes along the wind velocity.

Figure 9 .
Figure 9. Order and magnitude of six initial modes along the wind velocity.
(a) One TMD is placed at the middle point.(b) Two TMDs are placed at the

Figure 10 .
Figure 10.Variation of the maximum torsion of (a)  2 point and (b) whole span, in relation to the mass ratio, in the one TMD condition.

Figure 10 .
Figure 10.Variation of the maximum torsion of (a) L 2 point and (b) whole span, in relation to the mass ratio, in the one TMD condition.

Figure 11 .
Figure 11.Variation of the maximum torsion of (a) L/2 point and (b) whole span, in relation to the distance factor, in the one TMD condition.
Figures 10 and 11 discloses that the   parameter can reduce the other modes' responses, or vice versa to mass ratio that cannot change the other modes' responses.Figures 12 and 13 demonstrate the sensitivity of responses in relation to the variation of mass ratio and distance factor for the two TMDs condition.

Figure 11 .
Figure 11.Variation of the maximum torsion of (a) L/2 point and (b) whole span, in relation to the distance factor, in the one TMD condition.

Figure 12 .
Figure 12.Variation of the maximum torsion of (a)  4 and 3 4 points and (b) whole span, in relation to the mass ratio, in the two TMDs condition.

Figure 12 .
Figure 12.Variation of the maximum torsion of (a) L 4 and 3L 4 points and (b) whole span, in relation to the mass ratio, in the two TMDs condition.

Figure
Figure12adisplays the main role of mass ratio in response reduction of first and second modes corresponding responses.Figure12bshows that by increasing the mass ratio, the second-mode responses are intensified, while first-mode ones are decreased.Hence, two TMDs, located at the Figure 12a displays the main role of mass ratio in response reduction of first and second modes corresponding responses.Figure 12b shows that by increasing the mass ratio, the second-mode responses are intensified, while first-mode ones are decreased.Hence, two TMDs, located at the  4 and 3 4points, provide positive local performance in all considered modes; however, their total one is not positive for all modes.

Figure 13 .
Figure 13.Variation of the maximum torsion of (a)  4 and 3 4 points and (b) whole span, in relation to the distance factor, in the two TMDs condition.

Figure 13 .
Figure 13.Variation of the maximum torsion of (a) L 4 and points and (b) whole span, in relation to the distance factor, in the two TMDs condition.

Figure 14 .
Figure 14.Variation of the maximum torsion of (a)  2 point and (b) whole span, in relation to the damping ratio, in the one TMD condition.Figure 14.Variation of the maximum torsion of (a) L 2 point and (b) whole span, in relation to the damping ratio, in the one TMD condition.

Figure 14 .
Variation of the maximum torsion of (a) L 2 point and (b) whole span, in relation to the damping ratio, in the one TMD condition.

Figure
Figure14adiscloses that by increasing the damping ratio the L 2 point corresponding responses in modes two, three, and five are considerably reduced.Of course, the second mode response reduction gradient is too high compared to other ones in the 0 up to 10 percent range.Figure14bindicates that the effect of the damping ratio on the maximum response of the whole bridge is too low compared to any certain point like L 2 .Especially, for mode two, the local performance of TMD is 90% ( 4.3 → 0.4) , while the total one is 7.6% ( 5.2 → 4.8 ).Thus, the damping ratio parameter is effective in the local performance of TMD, and its total efficacy is nuance.Figure15demonstrates the sensitivity of responses in relation to the variation of damping ratio for the two TMDs condition.Like the one TMD condition, the damping ratio renders suitable local performance, and its appropriate range is 0 up to 10 percent.

Figure
Figure 14a discloses that by increasing the damping ratio the  2

Figure 15 .
Figure 15.Variation of the maximum torsion of (a)  4 and 3 4 points and (b) whole span, in relation to the damping ratio, in the two TMDs condition.

𝐿 2 point
are declined by the increment of the frequency ratio.The final part of the second mode's corresponding curvature indicates that when the tuning frequency deviates considerably from the frequency

Figure 15 .
Figure 15.Variation of the maximum torsion of (a) L 4 and 3L 4 points and (b) whole span, in relation to the damping ratio, in the two TMDs condition.

Figure 16 .
Figure 16.Variation of the maximum torsion of (a)  2 point and (b) whole span, in relation to the frequency ratio, in the one TMD condition.

Figure 17
Figure 17 demonstrates the sensitivity of responses in relation to the variation of frequency ratio for the two TMDs condition.It is worth noting that the mass ratio,   , and damping ratio are considered 1%, 0.2, and 10%, respectively.

Figure 17 .
Figure 17.Variation of the maximum torsion of (a)  4 and 3 4 points and (b) whole span, in relation to the frequency ratio, in the two TMDs condition.
the VIV responses are provided in Tables

Figure 16 .
Figure 16.Variation of the maximum torsion of (a) L 2 point and (b) whole span, in relation to the frequency ratio, in the one TMD condition.

Figure 17
Figure 17 demonstrates the sensitivity of responses in relation to the variation of frequency ratio for the two TMDs condition.It is worth noting that the mass ratio, r B , and damping ratio are considered 1%, 0.2, and 10%, respectively.

Figure 16 .
Figure 16.Variation of the maximum torsion of (a)  2 point and (b) whole span, in relation to the frequency ratio, in the one TMD condition.

Figure 17
Figure 17 demonstrates the sensitivity of responses in relation to the variation of frequency ratio for the two TMDs condition.It is worth noting that the mass ratio,   , and damping ratio are considered 1%, 0.2, and 10%, respectively.

Figure 17 .
Figure 17.Variation of the maximum torsion of (a)  4 and 3 4 points and (b) whole span, in relation to the frequency ratio, in the two TMDs condition.
the VIV responses are provided in Tables

Figure 17 .
Figure 17.Variation of the maximum torsion of (a) L 4 3L 4 points and (b) whole span, in relation to the frequency ratio, in the two TMDs condition.

Figure 18 .
Figure 18.Standard deviation in uncontrolled and controlled conditions.Figure 18.Standard deviation in uncontrolled and controlled conditions.

Figure 18 .
Figure 18.Standard deviation in uncontrolled and controlled conditions.Figure 18.Standard deviation in uncontrolled and controlled conditions.

Figure 18 .
Figure 18.Standard deviation in uncontrolled and controlled conditions.

Figure 19 .
Figure 19.Total response of bridge in uncontrolled and controlled conditions.

Figures 20 and 21
Figures 20 and 21 show the time history of torsional responses.Figures 20 and 21 show that TMDs successfully mitigate the VIV in considered modes.

Figure 19 .
Figure 19.Total response of bridge in uncontrolled and controlled conditions.

Figures 20 and 21
show the time history of torsional responses.Figures20 and 21show that TMDs successfully mitigate the VIV in considered modes.

Figure 20 .
Figure 20.Time history of torsional response at the 4   velocity at each mode critical point.

Figure 21 .
Figure 21.Time history of torsional response at the resonance velocity of each mode at its critical point.

Figure 20 .
Figure 20.Time history of torsional response at the 4 m s velocity at each mode critical point.

Figure 20 .
Figure 20.Time history of torsional response at the 4   velocity at each mode critical point.

Figure 21 .
Figure 21.Time history of torsional response at the resonance velocity of each mode at its critical point.Figure 21.Time history of torsional response at the resonance velocity of each mode at its critical point.

Figure 21 .
Figure 21.Time history of torsional response at the resonance velocity of each mode at its critical point.Figure 21.Time history of torsional response at the resonance velocity of each mode at its critical point.

Table 1 .
Structural and geometrical specifications of considered bridge.

Table 2 .
Mode shapes, frequencies, and periods of considered bridge.

Table 2 .
Mode shapes, frequencies, and periods of considered bridge.

Table 2 .
Mode shapes, frequencies, and periods of considered bridge.
4.1.TMD DeviceTuned mass dampers have been widely employed in numerous studies for vibration control

Table 4 .
Optimum parameter of TMD reached by genetic algorithm.

Table 5 .
Optimum parameter of TMD reached by sensitivity analysis.