Parametric Study on Responses of a Self-Anchored Suspension Bridge to Sudden Breakage of a Hanger

The girder of self-anchored suspension bridge is subjected to large compression force applied by main cables. So, serious damage of the girder due to breakage of hangers may cause the collapse of the whole bridge. With the time increasing, the hangers may break suddenly for their resistance capacities decrease due to corrosion. Using nonlinear static and dynamic analysis methods and adopting 3D finite element model, the responses of an actual self-anchored suspension bridge to sudden breakage of hangers are studied in this paper. The results show that the sudden breakage of a hanger causes violent vibration and large changes in internal forces of the bridge. In the process of the vibration, the maximum tension of hanger produced by breakage of a hanger exceeds 2.22 times its initial value, and the reaction forces of the bearings increase by more than 1.86 times the tension of the broken hanger. Based on the actual bridge, the influences of some factors including flexural stiffness of girder, torsion stiffness of girder, flexural stiffness of main cable, weight of girder, weight of main cable, span to sag ratio of main cable, distance of hangers, span length, and breakage time of hanger on the dynamic responses are studied in detail, and the influencing extent of the factors is presented.


Introduction
Large-span bridges, including cable-stayed bridge, suspension bridge, and arch bridge, need cables to be the stay cables, main cables, or hangers. With the time increasing, the steel elements of bridge are exposed to corrosion, and their resistance capacities decrease. The stay cables and hangers especially are more prone to corrosion than other elements for the diameters of their steel wires are small. Under the combination action of live loads and corrosion, stay cables and hangers may break suddenly. In china, serious corrosion and breakage of cables occurred in many bridges [1]. In 2001, eight hangers of Yibin Southgate Bridge in China broke, and the deck supported by these hangers fell into the river. Two men died in this accident. In 1995, one stay cable of Guangzhou Haiyin Bridge that was only used for 7 years broke, and it hit an oil tank truck. To avoid the breakage of cables due to corrosion, cables in many bridges were replaced in China. For example, all the stays in Jiujiang Bridge that was used for 10 years were replaced, because 70% of the stays were seriously corroded, and 1/3 of wires in some stays were broken.
Because the breakage of cable usually occurs suddenly, and it can cause strong vibration and large change of internal forces of the structure, the sudden breakage may endanger the safety of the bridge. Because one hanger is broken, Mahakam II suspension Bridge in Indonesia collapsed in 2011, as shown in Figure 1. It happened very fast, only about 30 seconds. This bridge was completed in 2002, and it was less than 10 years old. At least 11 people were dead and 30 people were missing in this accident.
The effects of breakage of stay cables on cable-stayed bridge have been studied by many researchers [2][3][4][5][6][7][8]. Considering the influences of different layout of stays, number of planes of stays, and stiffness of deck, Mozos and Aparicio [5,6] studied the effects of breakage of stays on deck, tower, and stays in detail. The study showed that it was safe for stays when adopting the dynamic amplification factor (DAF) of 2.0 recommended by PTI [9], but it was unsafe for deck and tower. Because the time that the breakage occurs in (breakage time) has significant effects on the responses of structure, Mozos and Aparicio [10] studied the breakage time through experiments and found that the breakage time was 0.00375 s for damaged cables, and 0.0085 s for undamaged cables. Zhou   and Chen [8] studied dynamic responses of a cable stayed bridge caused by abrupt cable-breakage considering dynamic bridge-vehicle interactions.
Considering the geometric and material nonlinearity, Cai et al. [11] studied the nonlinear responses and progressive collapses of cable-stayed bridge due to sudden breakage of stays. When a single stay broke, the tension of other stays changed little. When two adjacent stays broke simultaneously, the bridge collapse did not occur, but the deck had large plastic deformation and the stays yielded. Qu et al. [12] studied DAFs of Sutong Bridge in China and found that the DAFs did not exceed 3%.
Ruiz-Teran and Aparicio [13] studied the effects of breakage of stays on underdeck cable-stayed bridge. The influences of the breakage time and the dynamic loads of breakage were investigated. When the breakage time was less than 10% of the fundamental period of the damaged bridge, it had no influences on the results. The study results showed that the responses produced by breakage of stays must be calculated using dynamic analysis method.
Considering the geometric and material nonlinearity, Qiu et al. [14] and Kao et al. [15] studied the static load-bearing capacity of self-anchored suspension bridge. The studies showed that the bridge did not collapse when five hangers broke. Because the dynamic effects of the breakage of hangers were not considered, the load-bearing capacities obtained by the studies are the upper limit value.
In recent years, more than 20 self-anchored suspension bridges have been built in China [16]. Their hangers are all made of parallel high strength galvanized steel wires. In Chinese codes [17], the safety factor of hanger is 3.0 during service, and it is 1.8 during replacement of hanger. No specifications are presented for sudden loss of hangers. Because the hangers are connected with the main cables, the breakage of hangers will induce the main cable to vibrate strongly, which will further induce strong vibration and large changes of internal forces of the whole bridge. Using nonlinear static and dynamic analysis methods, the responses of a self-anchored suspension bridge due to the breakage of hangers are studied in this paper. The results can be used for reference in design and maintenance of this kind of bridge.

Structure of a Self-Anchored Suspension
Bridge and Analysis Methods used in the girder. The tower is reinforced concrete with box cross-section. The bridge has two main cables and each cable is made of 3937 paralleled, 5 mm in diameter, high strength galvanized steel wires. The span to sag ratio of main cable in mid-span is 5.5. There are 65 hangers on each side of the bridge, and the distance of the hangers along the girder is 5 m. The hangers are made of 97 paralleled, 7 mm in diameter, steel wires. The hangers are numbered from 1 to 65.

Loads.
Because this study is only used to investigate the responses of bridge caused by sudden breakage of hangers, the analysis of the structure under dead loads is carried out. The self-weight of girder, towers, main cables, and hangers is obtained from multiplying the areas of their cross-sections by their densities. The distributed dead load on the deck except self-weight is 108 kN/m. The distributed dead load on the main cables except self-weight is 0.37 kN/m. Because the weight of clamps affects the vibration of the main cable, it is also considered in dynamic analysis. The weight of each clamp from 1 to 9 is 12 kN, the weight of each clamp from 10 to 17 is 16 kN, and the weight of each clamp from 18 to 33 is 10 kN.

Analysis Model.
The bridge is modeled with threedimensional (3D) finite element (FE) model, as shown in Figure 4. In the FE model, the girder, towers, and transverse beams are modeled using 3D beam elements, and their torsion stiffness and torsion mass are considered. In order to study the influence of flexural stiffness, the main cables are modeled using beam elements. The hangers are modeled using truss elements. The influences of initial internal forces under self-weight on the stiffness of the structure are considered. The main cables between two hangers are modeled using 5 elements to consider geometrical configuration and local vibration of the main cables more precisely. The masses of clamps are concentrated on the nodes where the clamps are located. Considering the piles of the bridge had little effects on the analysis of sudden breakage of hanger, they are not modeled in the FE model. The girder is supported by two bearings at each tower as well as at each abutment. Using this FE model, the modal properties of Zhuanghe Jianshe Bridge are calculated. The lowest natural frequency of the vertical vibration mode of the bridge is 0.41 Hz, and the lowest natural frequency of the torsional vibration mode of the girder coupled with transversal vibration is 1.00 Hz.

Analysis Methods.
Responses of the bridge due to sudden breakage of hangers on the bridge are analyzed by means of nonlinear static and dynamic analysis using the finite element software ABAQUS V.6.8. For either static analysis or dynamic analysis, the geometric nonlinearity of the structure and effects of axial forces of the structure on the stiffness are considered, and the nonlinear procedures are carried out using iteration method. In dynamic analysis, the direct time integration method is used. The time history of the tension of the broken cable is shown in Figure 5. Before the broken hanger breaks, the value of its initial tension is 0 . When the hanger breaks at time 1 , its tension drops from 0 to zero during time interval Δt.
Here, Δt is called breakage time of hanger, which takes a value of 0.005 s [5] except when studying the influences of the breakage time on the dynamic responses.
The following analysis processes are adopted for analysis, as shown in Figure 6. bending moments of the girder and towers are very small.
(b) A hanger element is removed, and its tension is unloaded in time interval Δt.
(c) The structure vibrates strongly due to the sudden breakage of the hanger, and the dynamic analysis is carried out to calculate the vibration of the structure. Because of damping, the vibration attenuates with time increasing. Until the maximum node displacement amplitude of the structure is less than 0.1 mm, the state of bridge is taken as the final equilibrium state under dead loads after breakage of a hanger.
The static state before breakage of the hanger is identified by 0 . After the sudden breakage of a hanger, the structure vibrates strongly, and the dynamic state of the bridge is identified by . It is especially mentioned that the final static state obtained by dynamic analysis is just the same as the static states of bridge without the broken hanger obtained by static analysis. The final static state after breakage is identified by .
After a preliminary numerical study on the dynamic responses of the bridges to pulse loads, the time step of 0.001 s is used to calculate the dynamic process during 10 seconds after the sudden breakage of hanger, and the time step of 0.005 s is used during the other time. The time steps allow us to achieve an important reduction in computing time and to maintain adequate accuracy in the results. A Rayleigh damping of 2% is used in the dynamic analysis.

Reponses due to Sudden Breakage of a Single Hanger
Considering that two hangers are not possible to break at the same time, the dynamic responses caused by breakage of a single hanger is studied in this paper. From the dynamic analysis results, it can be found that the sudden breakage of a hanger produces very strong vibration of the structure; some internal forces of the structure become very large during the vibration. The bending moments of girder, bending moments of towers, and tensions of main cables caused by the sudden breakage of a hanger are not large enough to control the Number of hangers Hanger tension (kN) 1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 T 0 T s T d,max design. The tensions of hangers and reactions of bearings are so large that they may control the design. Figure 7 shows the initial value 0 , final value , and maximum value ,max of tensions of the other hangers when hanger 23 breaks. It can be seen from the figure that breakage of a hanger has large effects on tensions of the hangers near the broken hanger and has little effects on tensions of the hangers far away from the broken hanger. The maximum tension ,max of each hanger during vibration is much larger than the corresponding final value .
To obtain the maximum values of tension ratios ,max / 0 and / 0 of every hanger, maximum tension ,max and final tension of each hanger are calculated when every one of the other hangers breaks. Figure 8 shows the maximum values of tension ratios ,max / 0 and / 0 of the hangers from 1 to 33. The tension ratio ,max / 0 of hanger 30 is the largest, that is, 2.22. However, the tension ratio / 0 of hanger 30 is only 1.48, which is only 66.7% of the tension ratio ,max / 0 . So, it can be concluded that the difference of the results between static analysis and dynamic analysis cannot be neglected. Additionally, the hanger is farther from the tower, the vibration caused by its breakage is stronger. So the tension of hanger increases with the distance between the hanger and tower increasing. The dynamic effects on hanger tensions are so marked that the maximum tension reaches 2.2 times of the initial value. If the breakage of a hanger is for the reason of corrosion, the other hangers may be also exposed to corrosion, and their resistance capacities decrease too. The breakage of one hanger may induce one-by-one breakage of the other hangers, and the bridge may collapse after several hangers break. So, the breakage of a hanger can endanger the safety of the bridge seriously.
The reaction forces ,max − 0 , ,min − 0 , and − 0 of the two bearings at tower T1 produced by sudden breakage of every one of the hangers from 1 to 33 on the left side are shown in Figures 9 and 10. Because breakage of any one hanger can induce strong vertical and torsional vibrations, vibrations further induce large changes of reaction forces of bearings. The maximum increment of reaction force of the left bearing is produced by breakage of hanger 13, and it is 1.86 times the tension of the broken hanger. The maximum decrement of reaction force of the right bearing is produced by breakage of hanger 6, and it is 1.35 times the tension of the broken hanger. Both the maximum increment and the maximum decrement of reaction force of the right bearing are produced by breakage of hanger 11, and they are 1.51 and 1.54 times the tension of the broken hanger, respectively. So, the reaction forces of bearings produced by breakage of a hanger cannot be ignored. Additionally, the reaction forces of bearings due to torsion of girder are related to the distance between the two bearings, and the reaction forces increase with the distance decreasing. A relatively large distance is adopted in the bridge analyzed in this paper.

Parametric Study on Influencing Factors
In order to find out the relationship between the responses caused by breakage of hanger and the structure, based on the bridge analyzed above, breakage of the left hanger H0 at the position of 1/4 of mid-span, as shown in Figure 11, is taken as an example to study the influences of flexural stiffness of girder, torsion stiffness of girder, flexural stiffness of main cable, weight of girder, weight of main cable, span to sag ratio of main cable, distance of hangers, span length, and breakage time of hanger on the dynamic responses. The mid-span length is 200 m for all above factors except span length. For general suspension bridge, the initial tensions of all hangers under dead load are designed to be a same value. In the suspension bridge analyzed in this paper, the initial tension of each hanger is 0 . To demonstrate the extent of the influences on hanger tensions, ratio of the maximum tension ,max to initial tension 0 is adopted in the following analysis. Because the breakage of hanger H0 has marked influence only on the adjacent hangers HL and HR, only the tensions of the two adjacent hangers are presented when the influencing factors are studied parametrically. Additionally, considering that the bearing reactions change largely during the dynamic responses, they are presented for a part of the influencing factors. Figure 12 shows relationship between tension ratio ,max / 0 of hangers HL and HR and flexural stiffness of girder when hanger H0 breaks suddenly. In the figure, is ratio of the changed flexural stiffness to the actual flexural stiffness of the designed girder. When the flexural stiffness of girder increases from 0.1 to 100 times of actual value, the tension ratio of hanger HL changes in the range of 2.085∼2.099, and it changes only by 0.67%. The tension ratio of hanger HR changes in the range of 2.093∼ 2.100, and it changes only by 0.33%. It can be concluded that flexural stiffness of girder has small effects on hanger tension. Figure 13 shows reactions of the four bearings at the two towers change with the flexural stiffness ratio increasing. It can be seen from the figure that the bearing reactions increase with the flexural stiffness of girder increasing. When the flexural stiffness ratio is between 0.01 and 3.0, the reactions caused by breakage of hanger H0 are nearly same except for the reaction of the left bearing on tower T2. But when the flexural stiffness ratio is larger than 3.0, the bearing reactions change very markedly. Particularly, when the flexural stiffness ratio is about 50.0, the peak values occur in all curves at the same time. The maximum reaction is 3946.5 kN, which is 2.31 times the initial tension of the broken hanger ( 0 =    Figure 14: Effects of torsion stiffness of girder on hanger tensions. Figure 14 shows relationship between tension ratio ,max / 0 of hangers HL and HR and torsion stiffness of girder when hanger H0 breaks suddenly. In the figure, is ratio of the changed torsion stiffness to the actual torsion stiffness of the designed girder. When the torsion stiffness ratio increases from 0.1 to 0.3, the tension ratio of hanger HL increases from 1.980 to 2.093, the tension ratio of hanger HR increases from 2.005 to 2.103, and they increase 5.7% and 4.9%, respectively. When the torsion stiffness ratio increases from 0.3 to 4.0, the hanger tension ratios decrease a bit. When the torsion stiffness ratio is larger than 4.0, the hanger tension ratios remain nearly unchanged. Figure 15 shows the trend that the bearing reactions increase with the increment of torsion stiffness of girder. When the torsion stiffness ratio increases from 0.01 to about 50, the reactions are nearly in the range of 750∼2000 kN. When the torsion stiffness ratio is larger than 50, the reactions increase markedly, and the maximum value is 4814.0 kN at the stiffness ratio of 100. is reasonable when the local vibration of main cable is not considered in structural analysis of suspension bridge. Because the flexural stiffness has large effects on the dynamic characteristics of the cable, it cannot be ignored in analysis of the violent local vibration of main cable caused by sudden breakage of hanger. In this paper, the main cable is assumed as a solid cylinder with curved geometry, which has the same sectional area as the actual main cable. By changing the flexural stiffness of main cable from 0.001 to 100 times the flexural stiffness of the solid cylinder, the effects of the flexural stiffness of main cable on the dynamic responses caused by breakage of hanger H0 are studied. Figure 16 shows the relationships between tension ratios of hangers HL and HR and the flexural stiffness ratio of main cable. It can be seen that the flexural stiffness of main cable has very marked influences on hanger tensions. When the flexural stiffness ratio of main cable increases from 0.001 to 0.01, the tension ratios nearly remain a constant value of about 2.10. When it increases from 0.01 to 0.05, the tension ratio decreases to be a value of about 2.01. When it increases from 0.05 to 1.5, the tension ratio increases to be a maximum value of 2.240. When it increases from 1.5, the tension ratio decreases rapidly. The tension ratio of hanger HR decreases to be a value of 1.675 when the flexural stiffness ratio is 100.

Weight of Structure.
Weight of structure is one of the primary factors that affect the dynamic characteristics of the bridge. By changing weight of girder and main cable, the effects of weight on the dynamic responses caused by breakage of hanger H0 are studied. Figure 17 shows the changes of tension ratios of hangers HL and HR when the weight of girder increases from 0.8 to 1.2 times the actual weight of the designed girder. The tension ratios change less than 0.2%, so the effects of weight of girder can be ignored. Figure 18 shows the relationship between tension ratios of hanger and the weight of main cable. It can be seen from the figure that the weight of main cable has larger effects than the weight of girder. When the weight of main cable  Figure 19 shows the relationship between tension ratios of hanger and the ratio of span to sag of main cable. It can be seen from the figure that the tension ratio of hanger decreases with the ratio of span to sag increasing. When the span to sag ratio increases from 4 to 12, tension ratio of hanger HL decreases from 2.088 to 2.023, and tension ratio of hanger HR decreases from 2.097 to 2.056. They decrease by 3.1% and 2.0%, respectively.

Distance of Hanger.
When the span length of the bridge remains unchanged, the effects of distance of hanger on the dynamic responses caused by sudden breakage of hanger H0 are studied. Figure 20 shows the relationship between tension ratios of hanger and the distance of hanger. The tension ratio of hanger increases with the distance of hanger increasing. When the distance of hanger increases from 5 m to 15

Span
Length. The mid-span length of the actual bridge is 200 m. In order to study the effects of span length on the dynamic responses caused by sudden breakage of hanger, the mid-span length is changed from 120 m to 400 m. In the parametric study, the distances of hangers remain to be 5 m. Figure 21 shows the relationship between tension ratios of hanger and mid-span length. With the mid-span length increasing, the tension ratio of hanger increases firstly and then decreases. The minimum value of tension ratio of hanger HL is 2.077, and the maximum value is 2.113, which is larger than the minimum value by 1.7%. The minimum value of tension ratio of hanger HR is 2.087, and the maximum value is 2.119, which is larger than the minimum value by 1.5%. Figure 22 shows effects of breakage time Δt of hanger on hanger tension ratios. The breakage time of hanger affects the tension ratio of hanger very markedly. The tension ratio of hanger decreases rapidly with the breakage time increasing. When the breakage time increases from 2 ms to 70 ms, tension ratio of hanger HL The Scientific World Journal decreases from 2.118 to 1.607, and tension ratio of hanger HR decreases from 2.135 to 1.612. The tension ratios of the two hangers increase by 24.1% and 24.5%, respectively. So, the breakage time of hanger should be studied in detail through experimental studies to make the numerical study be more accurate.

Conclusions
For self-anchored suspension bridge, its main cables are anchored directly to two ends of its stiffening girder, and the stiffening girder subjected to a very large compression force is supported by hangers. So the hangers are very important for the safety of the whole bridge. With the time increasing, the resistance capacities of the hangers decrease due to corrosion, and the hangers may break suddenly. The study on the responses of a self-anchored suspension bridge to breakage of hangers reaches the following conclusions.
(1) The sudden breakage of a hanger produces a very strong vibration and large changes of internal forces of the bridge. During the vibration, the maximum tension of hanger reaches 2.22 times the initial value. The maximum increment of reaction force of bearing is 1.86 times the tension of the broken hanger, and the maximum decrement of reaction force of bearing is 1.54 times of tension of the broken hanger.
(2) Within the studied influencing factors, the flexural stiffness of main cable, distance of hangers, and breakage time of hanger have very significant effects on the dynamic responses caused by sudden breakage of a hanger. Their effects on hanger tension ratio are larger than 10%. Torsion stiffness of girder, weight of main cable, span to sag ratio of main cable, and span length have significant effects on the dynamic responses caused by sudden breakage of a hanger, and their effects on hanger tension ratio are in the range of 2∼5%. Flexural stiffness and weight of girder have very small effects on the dynamic responses caused by sudden breakage of hanger, and their effects on hanger tension ratio are less than 1%. Some of the above influencing factors have marked effects on reactions of the bearings. Particularly, when the factors take some values, the torsion vibrations of girder are very strong, and the bearing reactions occur to be very large values, which are several times of the tension of broken hanger.