# Simple Program to Investigate Hysteresis Damping Effect of Cross-Ties on Cables Vibration of Cable-Stayed Bridges

Simple Program to Investigate Hysteresis Damping Effect of Cross-Ties on Cables Vibration of Cable-Stayed Bridges Hindawi Publishing Corporation Home Journals About Us About this Journal Submit a Manuscript Table of Contents Journal Menu Abstracting and Indexing Aims and Scope Annual Issues Article Processing Charges Articles in Press Author Guidelines Bibliographic Information Contact Information Editorial Board Editorial Workflow Free eTOC Alerts Reviewers Acknowledgment Subscription Information Open Focus Issues Focus Issue Guidelines Open Special Issues Published Special Issues Special Issue Guidelines Abstract Full-Text PDF Full-Text HTML Full-Text ePUB Linked References How to Cite this Article Complete Special Issue Advances in Civil Engineering Volume 2012 (2012), Article ID 463134, 15 pages doi:10.1155/2012/463134 Research Article Simple Program to Investigate Hysteresis Damping Effect of Cross-Ties on Cables Vibration of Cable-Stayed Bridges Panagis G. Papadopoulos , Andreas Diamantopoulos , Haris Xenidis , and Panos Lazaridis Department of Civil Engineering, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece Received 11 October 2011; Revised 21 March 2012; Accepted 24 March 2012 Academic Editor: Husam Najm Copyright © 2012 Panagis G. Papadopoulos et al. This is an open access article distributed under the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract A short computer program, fully documented, is presented, for the step-by-step dynamic analysis of isolated cables or couples of parallel cables of a cable-stayed bridge, connected to each other and possibly with the deck of the bridge, by very thin pretensioned wires (cross-ties) and subjected to variation of their axial forces due to traffic or to successive pulses of a wind drag force. A simplified SDOF model, approximating the fundamental vibration mode, is adopted for every individual cable. The geometric nonlinearity of the cables is taken into account by their geometric stiffness, whereas the material nonlinearities of the cross-ties include compressive loosening, tensile yielding, and hysteresis stress-strain loops. Seven numerical experiments are performed. Based on them, it is observed that if two interconnected parallel cables have different dynamic characteristics, for example different lengths, thus different masses, weights, and geometric stiffnesses, too, or if one of them has a small additional mass, then a single pretensioned very thin wire, connecting them to each other and possibly with the deck of the bridge, proves effective in suppressing, by its hysteresis damping, the vibrations of the cables. 1. Introduction The pretensioned cables in a typical cable-stayed bridge of medium size [ 1 ], as they are very long with a length of magnitude order 100 m, and a pretension axial force of magnitude order 1000 kN, exhibit, perpendicularly to their axis, a very small geometric stiffness, corresponding to their fundamental vibration mode, of a magnitude order only 50 kN/m. Also perpendicularly to their axis, they exhibit a very small intrinsic damping, due to their material internal friction. For the previous reasons, they are often subjected to large amplitude vibrations. And, if the external excitation is approximately periodic, with a period close to a natural period of the cable, for example, the fundamental one, then resonance may happen, and vibration amplitudes increase excessively and are maintained, with no significant reduction for a long time, unless special measures are taken. Two usual reasons for the previous cable vibrations of cable-stayed bridges are the following. (1) A pretensioned cable exhibits a sag under its self-weight. Because of traffic, the ends of the cable, on pylon and deck, are subject to a variation of their displacements; thus the elongation and axial force of the cable vary, which implies variation of its geometric stiffness, too, as well as variation of the sag of the cable. This vibration, due to variation of geometric stiffness, is called parametric excitation. (2) The successive pulses of a wind pressure exert a drag force, perpendicularly to the vertical plane of cables, at one side of the bridge. The variation of this drag force causes vibration of the cables. In [ 2 ], a complete description of the problem of cable vibrations of cable-stayed bridges is presented, as well as a state of the art of various types of dampers for cable vibrations (viscous dampers, cross-ties, and others), along with case studies of dampers on real bridges. The viscous dampers, although widely mentioned in literature, present some problems: usually, they are installed at the ends of a cable, where they are not very helpful; it seems that their main role is a slight reduction of cable’s length, thus a slight increase of its geometric stiffness. Rarely, they are installed at intermediate points of a cable, where they are more helpful; however, this installation is difficult. On the other hand, the cross-ties are preferable, for the following reasons: they are light and cheap, they are easily installed and pretensioned, and they easily replaced when damaged. And a great advantage of them is that although they are very thin, with a ratio of cross-section area of the cable to that of the cross-tie of a magnitude order 1000, however, the axial elastic stiffness, of a single pretensioned cross-tie, is comparable in magnitude with the geometric stiffness of a cable, that is of magnitude order 50 kN/m, along the same direction, perpendicularly to cable axis. Also, as the cross-ties are very thin, they are almost invisible, so they do not harm the aesthetics of the bridge. For the previous reasons, recently many researchers recommend the use of cross-ties to suppress large amplitude cable variations of cable-stayed bridges. In [ 3 – 6 ], analytical studies are performed on cross-ties or hybrid systems consisting of viscous dampers and cross-ties. Here, a simplified analytical method is proposed [ 7 ], in order to investigate the hysteresis damping effect of cross-ties, where, for every individual cable, an SDOF model is adopted, approximating its fundamental vibration mode. The geometric nonlinearity of the cables is taken into account by their geometric stiffness. At same time, the proposed method is accurate, as it includes the material nonlinearities of the cross-ties, by their compressive loosening, tensile yielding, as well as hysteresis stress-strain loops. A short computer program (only about 120 Fortran instructions), fully documented, is presented for the step-by-step dynamic analysis [ 9 ] of isolated cables or couples of parallel cables, connected to each other and possibly with the deck of the bridge, by very thin pretensioned wires (cross-ties) and subjected to variation of the axial forces of cables due to traffic [ 8 ] or to successive pulses of wind drag force. Seven numerical experiments are performed. And, based on them, observations are made on the effectiveness of a single pretensioned very thin wire, connecting a couple of cables of a cable-stayed bridge, in suppressing, by its hysteresis damping, their large amplitude vibrations. 2. Equations of the Problem Figure 1 shows a part of a typical cable-stayed bridge [ 1 ], with a pylon consisting of two vertical legs, which are connected by two transverse beams, a part of the deck with a slender rectangular plate section, two couples of pretensioned parallel inclined cables at each side of the bridge, connected by their ends to the pylon and the deck and two pretensioned very thin vertical in-plane cross-ties, at two sides of bridge, which intend to suppress parametric vibration of cables due to traffic loads ๐ , as well as two out-of-plane horizontal cross-ties, which intend to suppress cable vibrations due to wind forces ๐ . Figure 1: Part of a typical cable-stayed bridge [ 1 ]. In the following, the equations of nonlinear dynamic analysis will be written, for a specific cable structure consisting of two parallel pretensioned cables (1 and 2 in Figure 1 ), connected by two very thin pretensioned cross-ties (3 and 4 in Figure 1 ) to each other and with the deck of the bridge and subjected to parametric vibration, due to traffic. In subsequent applications, by simple and obvious modifications of these equations, other cable structures (combinations of main cables and cross-ties) subjected to traffic or wind excitation will be treated, too. In the following analysis, the inclination of the cables will be ignored for the sake of simplicity. So, in Figure 2 (a), the couple of horizontal parallel pretensioned main cables represents the inclined cables 1 and 2 of Figure 1 , whereas the very thin pretensioned vertical wires represent the in-plane cross-ties 3 and 4 of Figure 1 . Figure 2: (a) Geometric, static, and dynamic parameters of a cable structure. (b) Primary σ - ε diagram describing the nonlinear axial stress-strain law of a cross-tie. (c) Given time-history of axial forces of cables. 2.1. Geometric Equations For every individual cable, a simplified SDOF model is adopted, which approximates its fundamental vibration mode. This unique DOF, for every cable, is the displacement of its center perpendicularly to its axis, that is, the vertical displacements downwards ๐ข ๐ and ๐ข ๐ข in Figure 2 (a), for the upper and lower cables, respectively. The geometric equations, relating the displacements of cables with the elongations and strains of cross-ties, are the following, according to Figure 2 (a): Δ โ ๐ = โ ๐ − ๐ข ๐ + ๐ข ๐ข − โ o ๐ , Δ โ ๐ข = โ ๐ข − ๐ข ๐ข − โ ๐ข ๐ , ( 1 ) where Δ โ ๐ and Δ โ ๐ข are elongations of upper and lower cross-ties, respectively, โ ๐ and โ ๐ข their design (nominal) lengths, and โ o ๐ and โ ๐ข ๐ are their initial undeformed lengths. And the strains ๐ o and ๐ ๐ข of upper and lower cross-ties, respectively, are ๐ o = Δ โ ๐ โ o ๐ , ๐ ๐ข = Δ โ ๐ข โ ๐ข ๐ . ( 2 ) 2.2. Constitutive Equations Figure 2 (b) is the primary stress-strain diagram, which describes the axial nonlinear stress-strain law of a cross-tie, made of the same high-strength steel as the main cables. This σ - ε law includes compressive loosening, tensile yielding, as well as hysteresis stress-strain loops, resulting from the obvious in loading-unloading rule Figure 2 (b). There is only one constitutive variable, the plastic strain ๐ p l of the cross-tie. The present stress ๐ is an obvious function of present strain ε and present value of plastic strain ๐ p l : ๎ท ๐ = ๐ ๐ , ๐ p l ๎ธ , ( 3 a ) whereas the variation of plastic strain Δ ๐ p l can be expressed as a function of present strain ε , variation of present strain Δ ๐ , as well as present value of plastic strain ๐ p l , in an obvious manner by the loading-unloading-reloading rule of Figure 2 (b): Δ ๐ p l = Δ ๐ p l ๎ท ๐ , Δ ๐ , ๐ p l ๎ธ . ( 3 b ) 2.3. Static Equations The axial forces of the cross-ties are ๐ ๐ = ๐ o ๐ด ๐ค , ๐ ๐ข = ๐ ๐ข ๐ด ๐ค , ( 4 ) for the upper and lower ties, respectively, where ๐ด ๐ค is the cross-section area of the very thin wires (cross-ties), whereas the vertical nodal forces applied at the centers of the cables, upper and lower one, respectively, are, according to Figure 2 (a) ๐น ๐ = ๐บ ๐ − ๐พ ๐บ ๐ ๐ข ๐ + ๐ ๐ , ๐น ๐ข = ๐บ ๐ข − ๐พ ๐บ ๐ข ๐ข ๐ข − ๐ ๐ + ๐ ๐ข , ( 5 ) where the downwards direction has been taken as positive, ๐บ ๐ and ๐บ ๐ข are weights at centers of upper and lower cables, respectively, and ๐พ ๐บ ๐ = 2 ๐ ( ๐ก ) ๐ ๐ , ๐พ / 2 ๐บ ๐ข = 2 ๐ ( ๐ก ) ๐ ๐ข / 2 ( 6 ) are their geometric stiffnesses, where ๐ ( ๐ก ) is the given time-history of the axial forces of the cables. 2.4. Dynamic Equations Damping is ignored, as the material internal friction of the cables is meaningless. The vertical accelerations at the centers of upper and lower cable are ฬ ๐ข ๐ = ๐น ๐ ๐ ๐ , ฬ ๐ข ๐ข = ๐น ๐ข ๐ ๐ข , ( 7 ) where ๐ ๐ and ๐ ๐ข are lumped masses at centers of cables (Figure 2 (a)), whereas the upper dots mean derivation with respect to time. 2.5. External Excitation Within the input data of the problem, the time-history of external excitation is given, which is here the variation of axial forces of cables due to traffic. The function ๐ ( ๐ก ) is assumed to be described by a piecewise linear curve, as shown in Figure 2 (c). And within each time interval, between two successive nodes, a linear interpolation is performed, in order to find, from a specific time instant ๐ก , the corresponding axial force ๐ of the cables. 2.6. Initial Value Problem A state vector is introduced: ๐ฒ = { ๐ฎ ๐ฏ ๐ } , ( 8 ) consisting of the vertical displacements ๐ฎ = { ๐ข ๐ ๐ข ๐ข } (Figure 2 (a)) and velocities ๐ฏ = { ฬ ๐ข ๐ ฬ ๐ข ๐ข } of the centers of upper and lower cables, respectively, as well as of the constitutive variables ๐ = { ๐ ๐ p l ๐ ๐ข p l } , which are the plastic strains (Figure 2 (b)) of upper and lower cross-ties, respectively. By combining all the previous equations, ( 1 ) up to ( 8 ), a system of first-order ordinary nonlinear differential equations is obtained: ฬ ๐ฒ = ๐ช ( ๐ก , ๐ฒ ) , ( 9 a ) which, along with the initial value of the state vector ๐ฒ ( 0 ) = ๐ฒ ๐ ( 9 b ) for time ๐ก = 0 , and with sought function the time-history of the state vector ๐ฒ ( ๐ก ) , constitutes an initial value problem. 2.7. Proposed Algorithm For the step-by-step dynamic analysis (direct time integration) of the previous initial value problem of ( 9a ) and ( 9b ), the algorithm of trapezoidal rule is proposed: ๐ฒ ๐ + 1 = ๐ฒ ๐ + 1 2 ๎บ ๐ช ๎ท ๐ก ๐ ๐ฒ ๐ ๎ธ ๎ท ๐ก + ๐ช ๐ + 1 ๐ฒ ๐ + 1 ๎ธ ๎ป Δ ๐ก , ( 1 0 ) where ๐ and ๐ + 1 are two successive steps of the algorithm. This coincides with the algorithm of constant average acceleration of Newmark’s group of algorithms for step-by step dynamic analysis. The aforementioned algorithm is combined with a predictor-corrector technique, with two corrections per step, PE(CE) 2 , where, in this symbol, ๐ means prediction and ๐ถ correction of the state vector y , whereas ๐ธ means evaluation of the function ๐ช ( ๐ก , ๐ฒ ) of ( 9a ). In more detail, the proposed predictor-corrector technique can be written, within any ๐ th step of the algorithm, as follows: Prediction ๐ฒ ๐ ๐ + 1 = ๐ฒ ๐ ๎ท ๐ก + ๐ช ๐ ๐ฒ ๐ ๎ธ Δ ๐ก , ( 1 1 a ) First correction ๐ฒ 1 ๐ + 1 = ๐ฒ ๐ + 1 2 ๎บ ๐ช ๎ท ๐ก ๐ ๐ฒ ๐ ๎ธ ๎ท ๐ก + ๐ช ๐ + 1 ๐ฒ ๐ ๐ + 1 ๎ธ ๎ป Δ ๐ก , ( 1 1 b ) Second and final correction ๐ฒ ๐ + 1 = ๐ฒ ๐ + 1 2 ๎บ ๐ช ๎ท ๐ก ๐ ๐ฒ ๐ ๎ธ ๎ท ๐ก + ๐ช ๐ + 1 ๐ฒ 1 ๐ + 1 ๎ธ ๎ป Δ ๐ก . ( 1 1 c ) Thanks to the aforementioned predictor-corrector technique, no solving of algebraic system is needed, within each step of the algorithm. The stability criterion of the proposed algorithm is [ 9 ] ๐ m a x Δ ๐ก < 2 . 0 r a d , ( 1 2 ) that is, Δ ๐ก < ๐ m i n / ๐ ; otherwise a divergent solution results, whereas the accuracy criterion is at least ๐ m a x Δ ๐ก < 0 . 5 r a d , ( 1 3 ) that is, Δ ๐ก < ๐ m i n / 4 ๐ = ๐ m i n / 1 2 . 5 6 ; otherwise a significant accumulated truncation error appears, which is expressed as amplitude decay of the vibration, as well as period elongation. 3. Computer Program Based on the proposed algorithm of previous Section 2.7 , a simple and very short computer program has been developed, with only 115 Fortran instructions, consisting of the MAIN program (79 instructions) which performs the algorithm of step-by-step dynamic analysis, and of three subroutines: ( 1 ) Subroutine EVAL (17 instructions) which evaluates the present strain and stress state of the cable structure under consideration, ( 2 ) subroutine SE (9 instructions) which describes the nonlinear uniaxial stress-strain law of a cross-tie, and ( 3 ) subroutine NHIST (10 instructions) describing the given time-history of the external excitation, which is, here, the variation, with respect to time, of the axial force of cables due to traffic. A full documentation of the previous computer program is presented as Appendix, consisting of the description of program line by line in Appendix A.1 , of the complete list of Fortran instructions in Appendix A.2 , and the variables explanation in Appendix A.3 . The documentation of the computer program is completed by the series of seven applications in Section 4 . The program is particularly oriented to the specific third application of Section 4.3 , as already mentioned in the equations of problem in Section 2 . However, by simple and obvious modifications of the computer program, all the other numerical experiments, in Section 4 , can be treated, too. 4. Applications (Numerical Experiments) Seven applications (numerical experiments) follow, on the dynamic analysis of isolated pretensioned cables of a cable-stayed bridge or couples of parallel cables connected to each other and possibly with the deck of the bridge by very thin pretensioned single wires (cross-ties). Three of these cable structures are subjected to variation of axial forces of cables due to traffic (parametric excitation) and four of them are subjected to successive pulses of drag force due to a strong wind. As already mentioned, the previously presented algorithm, in Section 2 , is oriented only to the specific third application of Section 4.3 . However, by simple and obvious modifications of this algorithm, all the other applications present in Section 4 , which are presented herein after, have been analysed, too. 4.1. First Application: Isolated Cable Subject to Traffic As shown in Figure 3 , an isolated cable is subjected to a periodic variation of its axial force with a period equal to its fundamental one. Up to the fifth cycle that the excitation lasts, the vibration amplitude of the cable increases up to about 1.6 m and then remains constant; only a slight algorithmic damping is observed. Figure 3: First application: isolated cable subject to traffic. (a) Given geometry and loading. (b) Axial stress-strain diagram of high-strength steel. (c) Initial static analysis. (d) Parameters of dynamic analysis. (e) Given time-history ๐ ( ๐ก ) of axial force of cable. (f) Resulting time-history ๐ข ( ๐ก ) of vertical displacement at center of cable. 4.2. Second Application: Couple of Interconnected Cables Subject to Traffic The cable of first application is, in Figure 4 , connected to another shorter parallel cable by a thin cross-tie. Both cables are subjected to a periodic variation of their axial force with a period equal to the fundamental one of the cable system. Up to the fifth cycle that the excitation lasts, the vibration amplitude of the upper cable increases up to 1.4 m and that of the shorter lower cable up to 0.7 m. Then, both amplitudes are gradually reduced, the upper one up to 0.4 m and the lower one up to 0.3, in 10 sec. Because of the different lengths of the two cables, thus different geometric stiffness, masses, and weights, too, large differences of displacements at the ends of the cross-tie are obtained, thus large stress-strain loops with a total width Δ ๐ ≈ 7 6 mm/m, which are responsible for the significant hysteresis damping which is achieved. Figure 4: Second application: couple of interconnected cables, subject to traffic. (a) Given geometry and loading. (b) Crosssections of main cables and cross-tie. (c) Initial static analysis. (d) Parameters of dynamic analysis. (e), (f) Resulting time-histories of vertical displacements of centers of upper and lower cables. (g) Resulting time-history of axial force of cross-tie. (h) Resulting hysteresis stress-strain loops of the cross-tie. 4.3. Third Application: Couple of Cables Connected to Each Other and to Deck, Subject to Traffic The cable system of the second application is, in Figure 5 , supplied by one more thin cross-tie connecting the lower cable with the deck of the bridge. Because of the stiffness of this additional cross-tie, the maximum vibration amplitudes of both cables are significantly reduced, that of the upper cable to 0.7 m and that of the lower cable to 0.3 m. However, at same time, this reduction of displacements has a consequence less wide stress-strain loops in the upper cross-tie with total width Δ ๐ ≈ 3 8 mm/m and very thin stress-strain loops in the lower cross-tie with only Δ ๐ ≈ 6 mm/m, resulting in low values of hysteretic damping. Figure 5: Third application: couple of cables, connected to each other and to deck, subject to traffic. (a) Given geometry and loading. (b) Cross sections of main cables and cross-ties. (c) Initial static analysis. (d) Parameters of dynamic analysis. (e), (f) Resulting time-histories of displacements of upper and lower cables. (g), (h) Resulting time-histories of axial forces of upper and lower cross-tie. (i), (j) Resulting hysteresis stress-strain loops of upper and lower cross-ties. 4.4. Fourth Application: Isolated Cable Subject to Wind As shown in Figure 6 , an isolated cable is subjected to a resonant periodic wind drag force. Up to the fifth cycle that the excitation lasts, the vibration amplitude increases up to about 3.0 m and then remains constant; only a slight algorithmic damping is observed. Figure 6: Fourth application: isolated cable subject to wind. (a) Given geometry and loading. (b) Parameters of dynamic analysis. (c) Given time-history of wind drag force. (d) Resulting time-history of displacement of center of cable. 4.5. Fifth Application: Couple of Interconnected Cables Subject to Wind Two identical parallel cables are, in Figure 7 , interconnected by a thin cross-tie. A wind drag force acts on one cable only; initially, this cable exhibits larger displacements, but gradually the movement is transferred to the other cable, too. So, the displacements are divided by two, compared with those of previous fourth application. During the initial stage of displacements transfer from one cable to the other, stress-strain loops of the cross-tie with a total width of medium size Δ ๐ ≈ 3 3 mm/m appear. From this point on, as the two cables are identical and perform similar movements, no more yielding of the cross-tie appears, as shown in Figure 7 (g), thus no more stress-strain loops and hysteretic damping, too. Figure 7: Fifth application: couple of interconnected cables subject to wind. (a) Given geometry and loading. (b) Cross sections of main cables and cross-tie. (c) Initial static analysis. (d) Parameters of dynamic analysis. (e), (f) Resulting time-histories of displacements of two main cables. (g) Resulting time-history of axial force of cross-tie. (h) Resulting hysteresis stress-strain loops of the cross-tie. 4.6. Sixth Application: Couple of Interconnected Cables, Additionally Connected by Diagonals to Pylon and Deck, Subject to Wind The cable system of the fifth application is, in Figure 8 , supplied with diagonal ties connected with the pylon and the deck of the bridge. These diagonal ties offer a small additional stiffness, perpendicularly to the cables, which slightly reduces their displacements. However, at same time, this restriction of displacements further reduces the total width of stress-strain loops of the cross-tie to Δ ๐ ≈ 2 7 mm/m, thus reducing the hysteretic damping. Figure 8: Sixth application: couple of cables interconnected by cross-tie, additionally connected by diagonals to pylon and deck, subject to wind. (a) Given geometry and loading. (b) Cross sections of main cables, cross-tie and diagonals. (c) Initial static analysis. (d) Parameters of dynamic analysis. (e), (f) Resulting time-histories of displacements of two main cables. (g), (h) Resulting time-histories of axial forces of diagonal bars. (i) Resulting time-history of axial force of cross-tie. (j) Resulting stress-strain loops of the cross-tie. 4.7. Seventh Application: Couple of Interconnected Cables Subject to Wind, with Small Additional Mass on One Cable The cable system of fifth application is in Figure 9 , supplied by a small additional mass to the one of the two cables. So, the two cables have now different dynamic characteristics and they perform significantly different movements. As a consequence, large differences of displacements at the ends of the cross-tie result, thus wide stress-strain loops with a total width Δ ๐ ≈ 7 3 mm/m, which implies significant hysteretic damping. The displacement amplitudes of the cables are now only about one-fourth of those of the fifth application. Figure 9: (a) Seventh application: Couple of inter-connected cables subject to wind with additional small mass on one cable. (a) Given geometry, loading and additional small mass. (b) Comparison of cross-sections of cables and cross-tie. (c) Initial static analysis. (d) Parameters of dynamic analysis. (e), (f) Resulting time-histories of displacements of two cables. (g) Resulting time-history of axial force of cross-tie. (h) Resulting hysteresis stress-strain loops of the cross-tie. 5. Conclusions Cable vibrations of cable-stayed bridges have been examined. Either isolated cables or couples of parallel cables, connected to each other and possibly with the deck of the bridge, by a very thin pretensioned wire (cross-tie), have been considered. External excitation is either traffic, which causes displacements of cable ends on deck and pylon, thus variation of axial forces, geometric stiffnesses and sags of cables, too (parametric excitation), or successive pulses of drag force due to a strong wind, perpendicularly to a vertical cables’ plane at one side of the bridge. The proposed analytical model is on the one hand simplified, as an SDOF oscillator is adopted for every individual cable, approximating its fundamental vibration mode. However, on the other hand, the proposed analytical model is accurate, as it takes into account the geometric nonlinearity of the cables by their geometric stiffness; also it includes the material nonlinearity of the cross-ties by their compressive loosening, tensile yielding, and hysteretic stress-strain loops. The equations of the problem of dynamic analysis, oriented to a specific cable structure, have been written, consisting of the geometric, constitutive, static, and dynamic ones, as well as of the given time-history of the external excitation. By combining these equations, an initial value problem is obtained. For the step-by-step dynamic analysis of this problem, the algorithm of trapezoidal rule is proposed, combined with a predictor-corrector technique, with two corrections per step. So, no solving of algebraic system is required within each step of the algorithm. Based on the proposed algorithm, a short computer program has been developed, with only 115 Fortran instructions, consisting of the main program and three subroutines. A full documentation is given for this program, which means transparency of computation. Seven numerical experiments have been performed by the aforementioned program, three with variation of axial forces of cables due to traffic (parametric excitation) and four with successive pulses of drag force due to a strong wind. On the basis of previous series of numerical experiments some observations with practical usefulness are made. (These are not strict theoretical conclusions, but simple observations based on the results of numerical experiments.) It is confirmed by the series of numerical experiments, the great advantage of pretensioned cross-ties, that although they are very thin, with ratio of cross-section area of a cable to that of a cross-tie of magnitude order 1000, however, they possess an axial elastic stiffness comparable in magnitude to the geometric stiffness of cables, with magnitude order 50 kN/m, along the same direction, that is perpendicularly to cables axes. The in-plane cross-ties (within a vertical cables plane) are intended to suppress cables’ variations from parametric excitation due to traffic, whereas the out-of-plane cross-ties (transverse ones connecting cables at two sides of bridge) are intended to suppress cables vibrations from successive pulses of drag force due to a strong wind. General observation from all numerical experiments: in a couple of parallel cables connected to each other and possibly with the deck of bridge by cross-ties, even a single cross-tie proves effective by its hysteresis damping (due to stress-strain loops) in suppressing large amplitude cable vibrations under the following circumstances: if the two cables have different dynamic characteristics, for example, different lengths which imply different masses, weights, and geometric stiffnesses, too, or if one of them has a small additional mass. Appendices A. Documentation of the Proposed Computer Program In this appendix, documentation is given for the proposed computer program, for the step-by-step dynamic analysis of the third application, that of a couple of parallel cables connected to each other and to the deck of bridge by a thin cross-tie and subject to a variation of their axial forces. A.1. Description of Program Line by Line The description refers to the complete numbered list of Fortran instructions of Algorithms 1 , 2 , 3 , 4 . Algorithm 1: MAIN program. Algorithm 2: Subroutine EVAL. Algorithm 3: Subroutine SE. Algorithm 4: Subroutine NHIST. MAIN Program The first seven lines include nonexecutable statements. Particularly, in the three first lines, the COMMON instructions connect the MAIN program with the three subroutines, by their common variables. In the next 17 lines, 8 up to 24, the input data are read: geometric data and density of steel in lines 8-9, the parameters of σ - ε law of steel along with the pretension stresses of cables and cross-ties in line 16, and the time-history of axial forces of cables given by the coordinates of nodes of piecewise linear curve ๐ - ๐ก in lines 21–24. In lines 10–15 and 17–20, some simple preliminary calculations are performed to determine cross-section areas and pretension forces of cables and cross-ties, as well as masses and weights of cables and yield strain of steel. In lines 25–35, the initial characteristic equation of the cable structure is solved, so that to find its natural frequencies and periods. In lines 36-37, the time scale of the given time history of external excitation is expanded so that to obtain a period equal to the fundamental one ๐ m a x of the structure, in order to cause resonance, whereas in line 38, the minimum natural period ๐ m i n of the structure dictates the time-step length Δ ๐ก of the algorithm, so that to assure accuracy of computation. In lines 39–49, the initial conditions are established: time ๐ก = 0 in line 39, determination of initial static displacements of cables in lines 40-41, evaluation of undeformed lengths of cross-ties in lines 42–44, zero initial velocities of cables in lines 45-46, zero initial plastic strains of cross-ties in lines 47-48, and evaluation of initial strain and stress state of structure by calling subroutine EVAL in line 49. In line 50, any step of algorithm begins by increasing time ๐ก by Δ ๐ก . In line 51, by calling subroutine NHIST, the present value of axial force of cables is determined. In lines 52–55, the prediction of values of displacements and velocities is performed, and in lines 56–58 by calling subroutine EVAL the corresponding plastic strains and accelerations are found. So in lines 52–58 the prediction of state vector within a step of algorithm is performed. In lines 59–65 the first correction of value of state vector is made by use of trapezoidal rule. And in lines 66–70, the second and final correction. In lines 71–73, the output data of present step of algorithm are written (time ๐ก , axial force ๐ of cables, displacements ๐ข ๐ and ๐ข ๐ข of cables, axial forces ๐ ๐ and ๐ ๐ข of cross-ties, strains and stresses ε o - σ ο and ε u - σ u of cross-ties). In lines 74–79 if a maximum time has been exhausted, the algorithm is interrupted. Otherwise, we continue to the next step of the algorithm. Subroutine EVAL Lines 1–3 are nonexecutable statements. In lines 4–7, from the displacements of cables, the elongation, strain, stress by calling subroutine SE, and axial force of upper cross-tie are determined. In lines 8–11, the corresponding quantities are found for the lower cross-tie. In lines 12-13, the vertical nodal forces on centers of upper and lower cables are determined and in lines 14-15 the corresponding accelerations. Subroutine SE In lines 3–5, the new plastic strain of the cross-tie is found. In lines 6-7, the stress of the cross-tie is determined. Subroutine NHIST In line 5 is found the time interval where present time ๐ก is included. In line 6, a linear interpolation is performed between the two end-nodes of the previous time interval, in order to find the axial force ๐ of cables corresponding to present time ๐ก . A.2. List of Fortran Instructions MAIN Program See Algorithms 1 , 2 , 3 , and 4 . A.3. Variables Explanation MAIN Program ๐ด = ๐ ๐ × ๐ ๐ข coefficient of characteristic equation AC: cross-section area of a cable AW: cross-section area of a wire (cross-tie) ๐ต = − ( ๐พ 1 1 ๐ ๐ข + ๐พ 2 2 ๐ ๐ ) coefficient of characteristic equation ๐ถ = ๐พ 1 1 × ๐พ 2 2 − ๐พ 2 1 2 coefficient of characteristic equation ๐ท = ( ๐ 2 − 4 ๐ ๐ ) 1 / 2 coefficient of characteristic equation DC: cross-section diameter of a cable DENS: density of steel DHO: elongation of upper tie DHU: elongation of lower tie DW: cross-section diameter of wire (cross-tie) DT = Δ ๐ก , time steplength of algorithm ELAST: initial elasticity (Young) modulus EO: strain of upper tie EOPL: plastic strain of upper tie EOPLP: prediction of EOPL EOPL1: first correction of EOPL EU: strain of lower tie EUPL: plastic strain of lower tie EUPLP: prediction of EUPL EUPL1: first correction of EUPL EVAL: subroutine for evaluation of strain and stress state of the structure EW0: pretension strain of wires (cross-ties) EY: yield strain of steel FO: vertical nodal force at center of upper cable FU: vertical nodal force at center of lower cable GO: vertical acceleration at center of upper cable GOP: prediction of GO GO1: first correction of GO GU: vertical acceleration at center of lower cable GUP: prediction of GU GU1: first correction of GU HO: design (nominal) length (height) of upper cross-tie HO0: undeformed length (height) of upper cross-tie HU: design (nominal) length (height) of lower cross-tie HU0: undeformed length (height) of lower cross-tie LO: length of upper cable LU: length of lower cable MO: mass of upper cable MU: mass of lower cable ๐ : axial force of a cable NHIST: subroutine for given time-history of ๐ NK: ordinate of a node of piecewise linear curve ๐ - ๐ NN: number of nodes of piecewise linear curve ๐ - ๐ SO: stress of upper cross-tie STIF1: ๐พ 1 1 elements of stiffness matrix STIF2: ๐พ 2 2 elements of stiffness matrix STIF12: ๐พ 1 2 elements of stiffness matrix SU: stress of lower cross-tie SW0: pretension stress of wires (cross-ties) SY: yield stress of steel ๐ : time TIEIN: input file TIEOUT: output file TK: abscissa of a node of piecewise linear curve ๐ - ๐ TMAX: ๐ก m a x , maximum time of observation TO: axial force of upper cross-tie TU: axial force of lower cross-tie ๐ 1 : ๐ m i n extreme natural periods of structure ๐ 2 : ๐ m a x extreme natural periods of structure UO: vertical displacement of center of upper cable UOP: prediction of UO UO1: first correction of UO UU: vertical displacement of center of lower cable UUP: prediction of UU UU1: first correction of UU VO: vertical velocity of center of upper cable VOP: prediction of VO VO1: first correction of VO VU: vertical displacement of center of lower cable VUP: prediction of VU VU1: first correction of VU WO: weight at center of upper cable WU: weight at center of lower cable W1: ๐ m a x extreme natural frequencies of the structure W2: ๐ m i n extreme natural frequencies of the structure. Subroutine EVAL Only the following variable is different from those of MAIN program: SE: subroutine for stress-strain law of a cross-tie. Subroutine SE Only the following variables are different from those of MAIN program: ๐ธ : strain of a cross-tie EPL: plastic strain of a cross-tie ๐ : stress of a cross-tie. Subroutine NHIST All the variables are the same as in MAIN program. References R. Walther, B. Houriet, W. Isler, P. Moia, and J. F. Klein, Cable-Stayed Bridges , Th. Telford, London, UK, 2nd edition, 1999. E. Caetano, “Cable vibrations in cable-stayed bridges,” IABSE-Structural Engineering Documents, 2007. L. Caracoglia and D. Zuo, “Effectiveness of cable networks of various configurations in suppressing stay-cable vibration,” Engineering Structures , vol. 31, no. 12, pp. 2851–2864, 2009. View at Publisher · View at Google Scholar · View at Scopus L. Caracoglia and N. P. Jones, “Passive hybrid technique for the vibration mitigation of systems of interconnected stays,” Journal of Sound and Vibration , vol. 307, no. 3–5, pp. 849–864, 2007. View at Publisher · View at Google Scholar · View at Scopus L. Caracoglia and N. P. Jones, “In-plane dynamic behavior of cable networks. Part 1: formulation and basic solutions,” Journal of Sound and Vibration , vol. 279, no. 3–5, pp. 969–991, 2005. View at Publisher · View at Google Scholar · View at Scopus L. Caracoglia and N. P. Jones, “In-plane dynamic behavior of cable networks. Part 2: prototype prediction and validation,” Journal of Sound and Vibration , vol. 279, no. 3–5, pp. 993–1014, 2005. View at Publisher · View at Google Scholar · View at Scopus P. G. Papadopoulos, A. Diamantopoulos, P. Lazaridis, H. Xenidis, and C. Karayiannis, “Simplified numerical experiments on the effect of hysteretic damping of cross-ties on cable oscillations,” in 9th International Conference on Computational Structures Technology , Athens, Greece, September 2008. P. G. Papadopoulos, J. Arethas, P. Lazaridis, E. Mitsopoulou, and J. Tegos, “A simple method using a truss model for in-plane nonlinear static analysis of a cable-stayed bridge with a plate deck section,” Engineering Structures , vol. 30, no. 1, pp. 42–53, 2008. View at Publisher · View at Google Scholar · View at Scopus P. G. Papadopoulos, “A simple algorithm for the nonlinear dynamic analysis of networks,” Computers and Structures , vol. 18, no. 1, pp. 1–8, 1984. View at Scopus var _gaq = _gaq || []; _gaq.push(['_setAccount', 'UA-8578054-2']); _gaq.push(['_setDomainName', 'hindawi.com']); _gaq.push(['_trackPageview']); (function() { var ga = document.createElement('script'); ga.type = 'text/javascript'; ga.async = true; ga.src = ('https:' == document.location.protocol ? 'https://ssl' : 'http://www') + '.google-analytics.com/ga.js'; var s = document.getElementsByTagName('script')[0]; s.parentNode.insertBefore(ga, s); })();
http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png
Advances in Civil Engineering
Hindawi Publishing Corporation
http://www.deepdyve.com/lp/hindawi-publishing-corporation/simple-program-to-investigate-hysteresis-damping-effect-of-cross-ties-bMbk3UdbrC
/lp/hindawi-publishing-corporation/simple-program-to-investigate-hysteresis-damping-effect-of-cross-ties-bMbk3UdbrC