Lyapunov-Based Control for Suppression of Wind-Induced Galloping in Suspension Bridges
This paper investigates the suppression of galloping in a suspension bridge due to wind loads. The galloping phenomenon can be destructive due to the high-amplitude oscillations of the structure. Two controllers are proposed to generate the control force needed to suppress the vertical galloping in the suspended cables and in the bridge deck. SIMULINK software is used to simulate the controlled system. The simulation results indicate that the proposed controllers work well. In addition, the performance of the system with the proposed controllers is compared to the performance of the system controlled with a tuned mass damper.
Suspension bridges are one of the most popular types of bridges. Very long spans of suspension bridges are being planned . Long, suspended steel cables such as those used in suspension bridges are prone to vibrations induced by wind and moving loads [2–6]. When suspended cables are subjected to wind disturbance, the system behaves nonlinearly due to its flexibility and because there is coupling between the motion of the bridge deck and the motion of the suspended cable.
The wind forces exerted on the cables and the bridge deck are complex. The aerodynamic stability of the cables and the bridge deck depends on many parameters including the wind velocity, the shape, and the size of the cross-sections and the angles of attack. When the wind velocity exceeds a certain critical value, the oscillations of the cable and/or the bridge deck will increase until a steady state response is achieved. This can be quite destructive to the structure due to the large amplitude of the oscillations. The aerodynamic instability of the bridge deck could occur simultaneously in both the torsional and vertical vibrational modes; this is called flutter and occurs when the torsional stiffness is close in magnitude to the flexural stiffness . Aerodynamic instability can also occur exclusively in the vertical mode if the torsional stiffness is much larger than the flexural stiffness, and this form is called galloping. In this paper, we will consider only the galloping instability for a suspension bridge where changing the shape or size of the cross-sections is difficult.
The occurrence of galloping is difficult to predict. The effect of the galloping phenomenon can be quite destructive, due to the high-amplitude oscillations of the structure in the direction perpendicular to the mean wind direction . For structures with low damping such as suspended cables, suspension bridges, or tall buildings, galloping can easily occur at low, steady or unsteady wind speeds . The nonlinearity behavior of the suspended bridge under wind excitation can cause internal resonances which affect the safety of the structures [7, 10, 11].
Galloping of suspended cables due to resonant harmonic loading has been studied in [10, 11] while galloping due to wind effects was studied in . It has been shown by Abdel-Rohman in  that increasing the damping is a very effective factor in increasing the critical wind speed at which galloping occurs. Attempts to introduce passive or semiactive damping into stay cables were made in [13–18].
In this paper, we proposed two active control schemes to suppress the galloping due to wind loading in suspension bridges. The designs of both controllers are based on Lyapunov theory; the proposed controllers guarantee the asymptotic stability of the controlled suspension bridge.
The paper is organized as follows. The dynamic model of the suspension bridge subject to wind forces is presented in Section 2. In Section 3, a SIMULINK model is presented for each subsystem of the suspension bridge. In Section 4, two state feedback controllers are proposed. Simulation results are given and discussed in Section 5. Finally, some concluding remarks are given in Section 6.
In the sequel, we denote by the transpose of a matrix or a vector . We use to denote a positive- (negative-) definite matrix . Sometimes, the arguments of a function are omitted in the analysis when no confusion can arise.
2. Equations of Motion of the System
The derivation of the basic equations of motion of the suspended cables subject to wind forces can be found in [12, 17, 19, 20]. For ease of presentation, the dynamic model is rederived in this paper.
Using the displacement directions defined in Figure 1, the general equations of motion are where is the spatial coordinate along the cable curved length (which is ), is the time, is the horizontal coordinate along the cable span, is the equation of the cable static profile, , , and are, respectively, the displacements of the cable at location in the horizontal, vertical, and transversal directions, is the static tension, is the additional dynamic tension in the cable, is the damping coefficient, is the mass of the cable per unit length, is the gravitational force, and and are the external loading per unit length in the vertical and transverse directions, respectively.
The nonlinear strain-displacement relationship during the deformation of the cable  is given by where is the modulus of elasticity, is the cross-section area of the cable, is the deformed cable segment, and is the undeformed cable segment.
The variables and are defined such that
The two-hinged bridge deck is suspended from the cables using vertical hangers located at as shown in Figure 2. is the length of the bridge deck and is the active control force installed at the general location along the bridge.
Based on the assumption of a small curvature regime and neglecting the horizontal motion (which can be safely disregarded because there is zero longitudinal loading ), the equations of motion from the static equilibrium position of the cable and the bridge can be simplified as follows: where is the flexural rigidity of the bridge deck, and are, respectively, the damping and mass per unit length of the bridge deck, is the stiffness of the vertical hangers (each located at and ), is the vertical response of the bridge deck, and indicate the response of the bridge deck and the suspended cables at and , respectively, and the is the Dirac-Delta function which is used to introduce, into the differential equations, the active control force . The forces , , and are, respectively, the lateral and vertical wind forces on the suspended cable and the vertical wind force on the bridge deck.
Previous studies have shown that locating the control force at the most flexible location of the structure provides the most feasible and optimal performance of the controlled structure as was shown by Abel-Rohman in . Thus, we set (see Figure 2).
The displacement functions , , and are considered to be the contribution of the first modes of vibrations and are assumed to take the forms where and are, respectively, the cable first mode shapes in the transversal and vertical directions, which can be determined using linear theory of cables  and satisfy the following boundary conditions: where is a constant chosen to make at the midspan of the cable, . The parameter is computed from the following relations:
Here, the terms and denote the vertical weight on the cable and the horizontal component of the static tension , respectively. The force is obtained from the static equilibrium of the cable as follows: in which is the sag in the cable profile of which is given by the equation:
For a two-hinged bridge deck, the mode shape can be assumed to take the form
Substituting (2.5) into (2.4) and applying an integral transformation one obtains the equations of motion for the suspended cable in the transverse and vertical directions and for the vertical motion of the bridge deck as follows: in which is the damping ratio in the suspended cable, and are, respectively, the natural frequencies of the cable in and directions, is the damping ratio in the bridge deck, and is the natural frequency of the bridge deck.
Expressions of the forces , , and in the above equations are given in  as follows: where the coefficients and are functions of the basic wind speed .
3. Simulation of the Uncontrolled Suspension Bridge Model Using SIMULINK
The dynamic model of the uncontrolled suspension bridge is simulated using the SIMULINK software. At first, a SIMULINK model is built for each subsystem of the suspension bridge. Then, the SIMULINK model subsystems are combined to simulate the overall system.
3.1. SIMULINK Model for the Unsuspended Bridge
The response of the bridge alone is obtained by simulating (2.13) with , , and, set to zero. Galloping occurs when , and therefore the critical wind speed is calculated to be 10.39?(m/sec) . The SIMULINK model of the unsuspended bridge is shown in Figure 3.
3.2. SIMULINK Model for the Suspended Cable Alone
If the cable is constructed alone, then the transverse and vertical responses of the suspended cable are obtained by simulating (2.11) and (2.12) with and set to zero. The critical wind speed of the cables is calculated to be 7.30?(m/sec) . The SIMULINK model of the suspended cable is shown in Figure 4 while the wind forces are simulated in Figure 5.
3.3. Simulation Results of the Uncontrolled Suspension Bridge
In the case of a suspension bridge, the response of the uncontrolled suspension bridge is simulated using (2.11)–(2.13) with and set to zero. The SIMULINK model of the uncontrolled suspension bridge is shown in Figure 6.
Abdel-Rohman and Joseph  showed that the critical wind speed is approximately 10.4?(m/sec). Therefore, we simulated the response of the uncontrolled suspended bridge at wind speed 12?(m/sec) which causes galloping as shown in Figure 7. It can be seen that the uncontrolled suspension bridge response has galloping in the vertical direction of both the cable and the bridge deck.
In order to suppress galloping due to wind forces at mean wind speed ?(m/sec), in the cable and bridge vertical responses, a vertical cable between the bridge deck and the suspended cables is used to install a hydraulic actuator able to generate the required active control force on the bridge deck. The controller will generate two equal and opposite active control forces. In the next section, a controller design based on Lyapunov theory is proposed. Two state feedback controllers are provided: the first one is a nonlinear controller and the second one is a linear controller.
4. Controllers Design
In order to insure the controllability of the pair , the element (2,3) of the matrix is set to a constant , and, consequently, is subtracted from the second element in the vector.
Hence, the equation of the system in (4.1) can be written as
The objective of this paper is to design control schemes to improve the stability of the system by suppressing the oscillations in the suspension bridge under wind excitation.
Remark 4.1. The simulation results indicate that the nonlinear function in (4.4) is uniformly bounded, and hence it can be assumed that the nonlinear term satisfies the following cone-bounding constraint: where is a positive scalar and is the Euclidian norm of a vector.
Remark 4.2. It can be checked that the pair in (4.4) is controllable. Hence, the poles of the closed-loop system can be selected such that the response of the linear part of the system (i.e., ) is as desired.
4.1. Design of the First Controller
In this section, a nonlinear controller is used to suppress the oscillations in the suspension bridge due to wind loads. The control law is divided into two parts: a linear part and a nonlinear part. The linear part of the controller is designed by using the standard pole placement technique. The nonlinear part of the controller is designed to guarantee the asymptotic stability of the closed-loop system.
Let the matrix be such that where is the feedback gain vector obtained by the standard pole placement technique. Let the symmetric positive definite matrix be the solution of the following Lyapunov equation : where .
Proposition 4.3. The control law given by the following equations (4.8)–(4.10) when applied to the suspension bridge system described by (4.4) guarantees the asymptotic stability of the closed-loop system: with where is a positive design parameter and is the sign function:
Proof. Using (4.4), (4.8), and (4.9), it follows that
Consider the following Lyapunov function candidate:
Note that for and for . Taking the derivative of with respect to time and using (4.12), (4.10), (4.7), and (4.5), it follows that where and are the minimum and maximum eigenvalues of a matrix, respectively.
Therefore, it can be concluded that if the matrices and are selected such . Hence, the control scheme given by Proposition 4.3 guarantees the asymptotic stability of the closed-loop system.
4.2. Design of the Second Controller
To simplify the computations involved in the implementation of the first controller, a linear controller is designed for the suspension bridge under wind loads as described by (4.4). Again, the control law is divided into two parts. The first part of the controller is designed by using the pole placement technique. The second part of the controller is designed to guarantee the asymptotic stability of the closed-loop system.
Define the matrix is such that
Let the symmetric positive definite matrix be the solution of the following Lyapunov equation: where . Also, let a positive design parameter be such that
Proof. Using (4.4), (4.18), and (4.19), it follows that
Consider the following Lyapunov function candidate:
Note that for and for . Taking the derivative of with respect to time and using (4.21), (4.20), (4.16), and (4.5), it follows that The choice of guarantees that .
Therefore, it can be concluded that . Hence, the control scheme given by Proposition 4.4 guarantees the asymptotic stability of the closed-loop system.
5. Simulation Results
The simulation results of the controlled suspension bridge using the proposed controllers are presented and discussed. The system is simulated using the SIMULINK software. The data for the suspension bridge used in  is used here. The details of the bridge deck and the cable and wind parameters are given in Table 1. From this table, the parameters of the model of the suspension bridge can be calculated as given in the appendix. It is worth mentioning that the torsional mode frequency which is equal to 1.276?rad/sec is more than three times greater than the frequency of the vertical bending mode of 0.404?rad/sec. Thus, our assumption of ignoring the torsional frequency is reasonable.
In this example, it is assumed that the wind speed profile is a constant and equals to the basic wind speed . The coefficients and can then be determined as functions of the basic wind speed as given in the appendix. In all simulations, the following initial conditions are used ?m and ?m/sec.
5.1. Simulation Results Using the First Controller
The parameters of the first controller are , the feedback gain vector , and , where is the identity matrix of appropriate dimension. The feedback gain vector is computed by MATLAB using the linear part of the system: where clp is the desired closed-loop poles of the linear part of the system. The matrix which is the solution of the Lyapunov equation in (4.7) is computed by MATLAB as follows:
Any wind speed greater than 10.4?m/sec will cause galloping and can be used for simulation purposes. In this paper, a wind speed of 12?m/sec (equivalent to 43.2?Km/hr) is used. This wind speed equals the wind closed-loop poles of the linear suspension bridge as shown in Figure 7. The results are shown in Figure 8.
Figure 8(a) shows the vertical displacement of the suspended cable, . It can be seen that the oscillations in the response decay to zero in about 30 seconds. Figure 8(b) shows the vertical displacement of the bridge deck, . It can be seen that the oscillations in the response decay to zero in about 30 seconds. Figure 8(c) shows the transverse displacement of the suspended cable, . It can be seen that the response oscillates with a slow decaying envelop as time increases. Hence, it can be concluded that the first control scheme is able to greatly reduce the oscillations in and . Comparing Figures 7(c) and 8(c), we noticed that the controller did not have much of an effect on the amplitude of the transverse displacement oscillations while it reduced the frequency of the oscillations.
It is worth noting that the active control force was introduced only for the vertical displacements, and hence the oscillations of the cable and the bridge deck in this direction was suppressed, but the cable response in the transverse direction oscillates for some time before decaying to a constant amplitude. The transverse response can be improved if the bridge deck is supplemented by a horizontally sliding damper along the mean wind direction. However, for this example, the transverse response is stable irrespective of the value of the mean wind speed . Moreover, the plot of the active control force versus time is shown in Figure 9. It can be seen that the force stays within a reasonable range.
5.2. Simulation Results Using the Second Controller
The simulation results of the controlled suspension bridge using the second controller are presented in this section. The parameters of the controller are 106, used in first controller, and . The wind speed used in testing the performance of the proposed controller is 12 (m/sec). The results are shown in Figure 10.
Figure 10(a) shows the vertical displacement of the suspended cable, . It can be seen that the oscillations in the response decay to zero in about 30 seconds. Figure 10(b) shows the vertical displacement of the bridge deck, . It can be seen that the oscillations in the response decay to zero in about 30 seconds. Figure 8(c) shows the transverse displacement of the suspended cable, . It can be seen that the response oscillates with a slow decaying envelop as time increases. Hence, it can be concluded that the first control scheme is able to greatly reduce the oscillations in and . Comparing Figures 7(c) and 10(c), we noticed that the controller did not have much of an effect on the amplitude of the transverse displacement oscillations.
The plot of the active control force versus time is shown in Figure 11. It can be seen that the force stays within a reasonable range.
Although we have almost the same performance when applying the two controllers, the second controller (which is linear) is much simpler in the implementation than the first controller. Thus, one can suggest applying the second controller.
5.3. Comparison with a TMD Controller
In this section, a comparison with the TMD controller is presented. A tuned mass damper with a weight 5% of that of the first mode mass of the bridge deck, a natural frequency of 0.97, and a damping ratio of 20% is used. The tuned mass damper TMD is added to the bridge deck at . The wind speed used in testing the performance of the three controllers is 12 (m/sec). The results are shown in Figures 12 and 13.
It can be seen from these figures that the proposed controllers gave better results when compared to the passive TMD controller. The proposed controllers have smaller oscillation amplitudes and faster decay rates. For further comparison, the norms of the vertical displacements of the suspended cable and the bridge deck are computed and listed in Table 2.
It can be seen that the linear controller gives the smallest norm among the three controllers. This gives more advantage to the second controller.
The control of the nonlinear oscillations of suspension bridges due to wind loads is investigated in this paper. In order to control the vertical oscillations of the suspended cables and the bridge deck, a hydraulic actuator can be installed between the bridge deck and the suspended cables. This actuator is used to generate an active control force on the bridge deck. A linear and a nonlinear control schemes are presented to generate the active control forces. These controllers guarantee the asymptotic stability of the closed-loop system. The performance of the controlled system is investigated through simulations using the SIMULINK software. The simulation results indicate that the proposed control schemes work well. Moreover, simulation results indicated that the proposed controllers gave better results when compared to the TMD controller.
A. Parameters of the Suspension Bridge Model
B. Coefficients of Wind Force
C. Coefficients of Wind Force
M. Abdel-Rohman, “Effect of unsteady wind flow on galloping of tall prismatic structures,” Nonlinear Dynamics, vol. 26, no. 3, pp. 231–252, 2001.View at: Google Scholar
A. H. Nayfeh and M. Abdel-Rohman, “Analysis of wind excited vibrations of cantilever beams using the method of multiple scales,” Journal of Sound and Vibration, vol. 144, pp. 87–93, 1991.View at: Google Scholar
P. Yu, Y. M. Desai, A. H. Shah, and N. Popplewell, “Three-degree-of-freedom model for galloping. Part 1: formulation,” Journal of Engineering Mechanics, vol. 119, no. 12, pp. 2404–2425, 1993.View at: Google Scholar
P. Yu, Y. M. Desai, N. Popplewell, and A. H. Shah, “Three-degree-of-freedom model for galloping. Part 2: solutions,” Journal of Engineering Mechanics, vol. 119, no. 12, pp. 2426–2448, 1993.View at: Google Scholar
A. Luongo and G. Piccardo, “Non-linear galloping of sagged cables in 1:2 internal resonance,” Journal of Sound and Vibration, vol. 214, no. 5, pp. 915–940, 1998.View at: Google Scholar
M. Abdel-Rohman and H. Askar, “Control by passive TMD of wind-induced nonlinear vibrations in cable stayed bridges,” Journal of Vibration and Control, vol. 2, no. 2, pp. 251–267, 1996.View at: Google Scholar
Y. L. Xu and Z. Yu, “Vibration of inclined sag cables with oil dampers in cable-stayed bridges,” Journal of Bridge Engineering, vol. 3, no. 4, pp. 194–203, 1998.View at: Google Scholar
E. A. Johnson, B. F. Spencer Jr., and Y. Fujino, “Semi active damping in stay cables—a preliminary study,” in Proceedings of the 17th International Modal Analysis Conference, pp. 417–423, Kissimmee, Fla, USA, February 1999.View at: Google Scholar
G. V. Rao and R. N. Iyengar, “Internal resonance and nonlinear response of a cable under periodic excitation,” Journal of Sound and Vibration, vol. 149, pp. 25–41, 1991.View at: Google Scholar
M. Irvine, Cable Structures, Dover, New York, NY, USA, 1992.
M. Abdel-Rohman, “Design of optimal observers for structural control,” Journal of Control Theory and Applications, vol. 131, no. 4, pp. 158–163, 1984.View at: Google Scholar
D. S. Naidu, Optimal Control Systems, CRC Press, Boca Raton, Fla, USA, 2003.
E. Simiu and R. H. Scanlan, Wind Effects on Structures: Fundamentals and Applications to Design, Dover, Mineola, NY, USA, 3rd edition, 2008.