Shock and Vibration

Volume 2014, Article ID 497670, 14 pages

http://dx.doi.org/10.1155/2014/497670

## Closed-Form Formula of the Transverse Dynamic Stiffness of a Shallowly Inclined Taut Cable

^{1}Department of Bridge Engineering, Room 711, Bridge Building, Tongji University, 1239 Siping Road, Shanghai 200092, China^{2}Shanghai Urban Construction Design Research Institute, Shanghai 200125, China

Received 1 December 2013; Revised 30 April 2014; Accepted 6 May 2014; Published 28 May 2014

Academic Editor: Mickaël Lallart

Copyright © 2014 Dan-hui Dan 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

The segmented vibration-governed equations and their general solutions for cables acted upon by intermediate transverse forces are derived by applying Hamilton’s principle. Including the effects of sagging, flexible stiffness, clamped boundary conditions, and inclination angle of the cable, the element-wise dynamic stiffness for each cable segment, split into segments having unique transverse forces, is derived. By using methods from the global stiffness assembly process of FEM, the global level of the cables’ dynamic equilibrium equation is obtained, and, as a result, the final closed-form formula of transverse dynamic stiffness is derived. Additionally, the corresponding analytic form, without considering sagging effects, is also obtained. Case studies are conducted on the aspects of accuracy, rationality of the distribution on the spatial field, and frequency domains of dynamic stiffness calculations. By comparison with the Guyan-based static FEM reduction method, it is shown that the result obtained from the proposed closed-form solution, which includes sagging effects, is exact and rational, thus creating a powerful tool in transverse vibration analysis.

#### 1. Introduction

A current trend in structure development is exceeding existing spanning limits, for which cable structures are playing an increasingly important role. With the extensive use of cable structures, the use of lateral bracing is becoming common, including the installation of intermediate turn supports, lateral dampers, and the use of cross ties, chain bar between cables, and hanger rods [1, 2]. These elements provide transverse static forces for the main body of the cable structures and are often specially designed for improving the dynamic characteristics of the cable structures [3, 4]. The analysis and calculation of characteristics and behaviors of cables acted upon by transverse forces exerted by these elements are also important factors in the design of a cable system. A direct relationship between the excitation and vibration response can be determined by an analysis of the dynamic stiffness, and, therefore, this particular analysis method has become an important tool in cable vibration analysis [5, 6].

Due to the nonlinear behavior of cables, the study on their dynamic stiffness is limited. Currently, most studies on the dynamic stiffness of cables are focused on the development of a computation method using a 4 × 4 matrix that describes the two-dimensional dynamic stiffness of one end of a cable, while the other end is kept fixed. Kolousek first provided the series solution to the dynamic stiffness at the free end of such cables, neglecting internal damping [7]. Davenport [8] (1959) found a closed-form dynamic stiffness matrix solution, taking damping into account and later, in 1965, developed the series solution of dynamic stiffness, which also considered damping effects [9]. In 1981, on the basis of the closed-form solution given by Davenport, Irvine proposed a simplified expression for dynamic stiffness [1]. All of these studies assumed small vibration amplitudes and a small inclination angle for the cables. In 1983, Veletsos and Darbrel studied a more accurate closed-form solution to the dynamic stiffness in the upper part of a cable under harmonic excitation, at all angles, and under a parabolic initial configuration. When damping effects were not taken into account, this closed-form solution was consistent with the series solution proposed by Veletsos and Darbre [10, 11]. In 1991, Starossek derived a dynamic stiffness matrix for cable structures with inclined hinges, which took sagging into consideration, and analyzed the features of its elements [12, 13]. Sarkar and Manohar proposed a calculation scheme for the dynamic stiffness matrix of cables using finite element analysis [14]. In 2001, Kim and Chang gave an approximate analytical expression of the dynamic stiffness matrix of an inclined cable with a catenary configuration [15], which surpassed the limitations of the parabolic static configuration.

Each of these studies on dynamic stiffness considered various factors based on the taut string theory, while the terms corresponding to the bending stiffness of the cables were all neglected in the passive vibration control equations. The actual cable structures are usually inclined cables with each end clamped and with relatively high bending stiffness [16, 17]. With increased length, the effect of sag cannot be neglected [18, 19]. A calculation algorithm of the dynamic stiffness of the cables, with all of these factors taken into account, would be very valuable to engineering practices. Furthermore, dynamic stiffness is used to describe the dynamical behavior at the ends of cables and, therefore, cannot be used in the case of a cable acted upon by transverse forces. However, in engineering practices, a majority of cable systems are acted upon by transverse forces, and, therefore, the study of transverse dynamic stiffness of cables is crucial and will contribute substantially to studies involving cable damping systems [20, 21].

In this study, a cable is modeled as being divided into various segments, each acted upon by a unique transverse force, while factors including sagging, bending stiffness, clamped boundary conditions, and inclination angle are all considered. The dynamic stiffness method is used to establish an equation of motion, and its general solution corresponds to a cable system acted upon by transverse forces. The dynamic displacement of each segment is used to characterize the dimensionless vibration mode function. Eventually, an analytical expression for the transverse dynamic stiffness matrix of a cable is given in a clear, simple form. The results of the study lay a foundation for further research in system vibration.

#### 2. General Solution for Governing Equation of Cable

The focus of this study is on a commonly observed situation in engineering, namely, when one end of a cable is acted upon by transverse forces. As shown in Figure 1, due to the presence of transverse forces, the cable system can be discretized into cable segment 1 (from end a to end c), cable segment 2 (from end c to end b), and transverse force elements (such as dampers).

As shown in Figure 1, by assuming that transverse forces do not affect the static configuration of the cable, the two cable segments can be described by the same quadratic parabolic function ; when the cable vibrates, segments 1 and 2 follow different configurations, namely, and . When considering the effects of sagging, bending stiffness, clamped boundary conditions, and inclination angle and ignoring the internal damping and shear forces, following governing equation [18, 22] can be used to describe the free vibration of the two segments: where , is the flexural rigidity, is the initial tension force, is the mass per unit length, is the configuration of the cable under static equilibrium (assuming that dampers have no impact on the static configuration), is the coordinate of each segment, and represents the additional segmented force, which is an additional cable tension generated due to the difference from the static configuration during vibration, and is equal to the product of the additional strain , generated during difference from the static equilibrium position, and the axial stiffness. That is,

The additional dynamic strain can be calculated using the ratio between the change in cable length during vibration and the static cable length . Consider

In (2) and (3), the superscript represents the dynamic configuration and represents static equilibrium. A similar rule is followed below.

As shown in Figure 1, in bridge projects, the angle between the tangential direction and the chord of a cable is very small for both the static and dynamic configurations. Therefore, when ignoring the higher order terms of the static and dynamic configurations, the amount of elongation of a differential segment during vibration relative to that of the static configuration can be approximated by

The total elongation of segment can be written as

With the assumption of a parabolic static equilibrium, one can use the following expression for the configuration: where is the sag-span ratio and .

Starting from (6) and noting the relationship of , we can easily get the coordinate derivative of the static configuration function , as follows: And then, by substituting it into (5), we have Similarly, we can get the expression of :

In the static configuration, the effective length of cable segments 1 and 2 is given by the following equation [23]: where is the arch length of the cable.

By integral operation, the sagging expression for the first cable segment can be written as follows, as shown in Figure 2: Thus, where is the relative position of transverse forces and .

From (9) and (12), the following formula is obtained:

Limited to the case when , the configuration of the dynamic displacement can be written as follows, using separation of variables: where is the exponential function based on natural logarithms, is the angular frequency of vibration, and is the vibration mode function on segment . By substituting (14) into (13), the following formula is obtained: where

It can be observed, by comparing the above equation with the equation given in literature [17, 18], that when the effects of transverse forces are considered, the formula for calculating the coefficients of additional cable forces is relevant to both the vibrational mode and the transverse forces.

By substituting (14) into (1), the ordinary differential equation for the vibrational mode function is obtained as follows:

Let The dimensionless vibration equation for the two cable segments is listed as follows: where is the ratio of the axial force to bending stiffness of the cable, is the dimensionless frequency of vibration, is the corresponding fundamental frequency of a taut cable, and is the corresponding fundamental frequency of a simple beam.

Equation (19) corresponds to the free vibration equation of cable segment , and its corresponding homogeneous equation can be written as

The general solution of this ordinary differential equation is well solved in references [18, 22], which can be written as the following expression

The special solution of (19), with the nonhomogeneous item being constant, can be written as

Finally, the complete solution of (19) can be written as where In addition, it can be obtained from (24) that

Equation (23) can be transformed to a matrix form as follows: where

This same method is adopted to make (16) dimensionless, so that

In the above equation, indicates a parameter determined entirely by geometric parameters of the cable and transverse force elements and is given by

By substituting (23) into (28), subsequent arrangement, and transposition, a particular solution that is independent of the vibration mode equation is obtained and can be subsequently combined with the general solution to obtain the particular solution of the vibration equation on cable segment , which is shown as follows: where matrix shows the impact of sagging on the dynamic system, which can be calculated using the following formula:

When the effects of sagging are considered, the solution for cable segment can be given in the form of the general solution and the particular solution; namely,

The amplitude of the additional cable force can be expressed as

#### 3. Dimensionless Vibrational Mode Function

The undetermined coefficients in the solutions of the vibrational mode function solution for two cable segments given by formula (32) can be determined through the boundary conditions.

As shown in Figure 3, the cable system is divided into two independent segments and a transverse force element. The dynamic nodal force and dynamic nodal displacements act on the boundaries of these models. After discretization, the nodal forces for each element include the amplitude of the dynamic nodal shear forces , , and and the amplitude of dynamic nodal bending moments , , and . The nodal displacements include the amplitudes of dynamic displacement , , and and the dynamic swing amplitudes , , and . The time-varying part of these quantities is represented by a complex exponent , where is the complex unit.

The boundary conditions of the displacement of the two cable segments are given by the following formula:

For cable segments , each of the nodal displacements and is used as boundary condition, respectively, to obtain the expressions for the undetermined coefficients ~. Equation (32) is substituted into (34) and (35), respectively. Consider where , , and are given by the following equation:

#### 4. Dynamic Stiffness Matrix of the Cable Damper System

All internal forces on cable segment can be expressed as follows:

The equilibrium conditions for the two cable segments are given by the following set of equations:

Equation (39) is transformed to matrix form according to the expressions for cable segments. Consider where the matrix is defined as follows:

The dynamic stiffness matrix of the cable segment is defined as follows:

For each cable segment, the relationship between the dynamic nodal force and the dynamic displacement is written in matrix form as follows:

After combining the two equations, the equation for the relationship between the total nodal force and the total displacement is established. Furthermore, the dynamic equilibrium equation for a cable acted upon by transverse forces, with clamped boundary conditions at both ends, is obtained: where the total dynamic stiffness matrix of the cable system with two degrees of freedom can be obtained by combining and . The assembly process can be written in matrix operator form as follows:

The matrix operator in the equation is as follows:

can be written in the form of matrix elements as follows: where , ; , ; ; , , and and , , and corresponds to , respectively.

It can be seen from (44) and its derivation that the equilibrium equation takes into account the impact of bending stiffness, clamped boundary conditions, and cable sagging. The factors considered here comprise a model that is most closely related to the cables used in actual projects. Therefore, in theory, this equation is able to describe more accurately the cable damper system in actual projects.

Similarly, by letting the above matrix be zero, the dynamic stiffness of each cable segment, not considering sagging, is obtained as follows:

The overall dynamic stiffness matrix of the cable can be obtained using (45). The dynamic equilibrium equation of the cable can be obtained by using (44).

The horizontal dynamic stiffness of the cable with a single degree of freedom can be obtained by rearranging (47) to the following form:

#### 5. Verification

To verify the accuracy of the two analytical methods for lateral dynamic stiffness, a finite element method is used to model the cable, and then a model reduction method is used to obtain the condensed stiffness and condensed mass corresponding to the degree of freedom. At this point, the corresponding transverse dynamic stiffness of cables is obtained.

Corresponding to the analytical method presented in this paper, the FEM tool-box, calfem-3.4, which is running at* Matlab* environment, is used to establish the finite element model. The cable and the damper are modeled by the two-dimensional plane beam element and a linear viscous damper model, respectively. The FEM model takes into account the effects of bending stiffness, inclination angle, clamped boundary conditions, and sagging, as well as the contributions of the cable forces and bending moments of the geometric stiffness matrix of cable. Reduction of mass and stiffness is conducted, for the two corresponding degrees of freedom, on the element acted upon by transverse forces using Guyan’s static condensation. Linear interpolation is then used to obtain the transverse dynamic stiffness matrix for the two degrees of freedom of the subordinate target point.

The basic parameters of the cable with both ends clamped which is used in the study are shown in Table 1.

##### 5.1. Transverse Dynamic Stiffness Matrix of Two Degrees of Freedom

According to the deduction process in the last section, the corresponding transverse static stiffness of the cable can be obtained by merely letting the circular frequency be . At this point, the dynamic stiffness at the estimated value of the modal frequency of the cable (provided by the formula of taut string theory in (50)) is investigated, which has a vital significance for the actual project.

Figure 4 shows the variations in static stiffness along the cable length using the three methods. The transverse axis represents the relative chord-wise length from the position of the transverse force to the anchor point of the inclined cable. The vertical axis represents the stiffness term. Due to the treatment of nodal displacements and nodal forces as described in the above section, the four terms have the same dimension. Table 2 shows a comparison of the values of the stiffness terms at the position when is 0.05, 0.1, 0.25, 0.5, 0.75, 0.9, and 0.95, respectively. Comparison of the four terms at the midpoint of the span is carried out.

As can be seen in Figure 4, the distributions of various terms of the two-dimensional dynamic stiffness matrix have similar patterns, in that the stiffness coefficient of a cable along the main diagonal of the matrix is extremely large at both ends of the cable. Moving the application point of the transverse force toward the midpoint of the cable causes this coefficient to decline rapidly. However, when varies within [0.1, 0.9], it does so smoothly, in the shape of a pan. The secondary diagonal terms transform rapidly from negative infinity to the smooth range in the middle of cable and then to positive infinity along the cable length.

The results obtained via the three presented methods for calculating the transverse static stiffness of the cable are consistent. The relationship between the values obtained in each method can be evaluated according to the following definition of error: where the subscript is . The superscripts , , , and , noSag represent the total stiffness terms obtained using the finite element method, by considering and neglecting sagging effects, respectively. As can be seen from Table 2 and Figure 5, the numerical differences between the stiffness terms obtained via the three presented methods are very small. Along the cable length, the error between the results of the first term on the main diagonal and the finite element method is less than 5.2%; is also small, within 11%. Though the relative error is given by Guyan’s static condensation using finite element analysis, the rigidity of the dynamic condensation itself also has errors. However, since the error of static stiffness is very small when , this case can be used as a reference of accuracy in the study. It is noted that, in this study, the obtained variation pattern of the transverse static stiffness of the cables along the cable length is reliable and has good computational accuracy. The dynamic stiffness method, taking sagging effects into account, is even closer to the static condensation results using finite element analysis.

Before analyzing the dynamic stiffness obtained using the three presented methods, the coefficient of dynamic stiffness is defined as follows: where is the estimated value of the first-order modal frequency of the cable according to the taut string equation, indicates the dynamic stiffness for a certain frequency, and represents the static stiffness.

Figure 6 shows the variations of the dynamic stiffness of the cable along its length, based on the estimated value of the first-order modal frequency using the three presented methods. The coefficient of the dynamic stiffness of corresponding to shows a distribution along the cable length with a shape of a positively laid pan, where the midpoint of the cable is a minimum, having a value of approximately 10% of the static stiffness. Conversely, shows a distribution with a shape of an inverted pan, where the midpoint of cable has a maximum value, getting closer to the static stiffness as it moves from the ends to the midpoint. The nearer one is to the center point, the closer the value is to the static stiffness. and , which correspond to the secondary diagonal, show discontinuities at the midpoint of cable. Their values approach positive and negative infinity, respectively.

The results of a comparison of the three presented methods can be drawn from Figures 6(a) and 6(b). The dynamic stiffness obtained using the method that considers sagging and obtained by the condensation method using finite element analysis yields results that are close to each other, whereas the method that does not take sagging effects into account produces results that are quite different. This shows that sagging effects have a significant impact on the dynamic stiffness of the cable near the first-order vibrational frequency and thus should be considered in calculations.

In some excitation tests of cables, the excitation points of even-order modes are often chosen to be at the midpoints of the cables. Therefore, it is crucial to consider the dynamic stiffness of the cables at the midpoint in engineering vibration tests. Figures 7(a) and 7(b) show the variation curve of the main diagonal terms of the dynamic stiffness matrix at the midpoint of a cable within the frequency domain (i.e., within ). As can be seen from the figure, the dynamic stiffness at the midpoint of the cable has an infinite discontinuity for even order. There is a zero crossing point in each continuous interval, namely, the odd-order modal frequency. The variation curves of the dynamic stiffness, with frequencies obtained using the three presented methods, are basically consistent. This shows that the calculated dynamic stiffness using the three methods is basically consistent within the investigated frequency range, with very small numerical differences. Figures 6 and 7 show that the calculated results given in this work for the dynamic stiffness at any frequency, using the analytical algorithm for dynamic stiffness, while considering sagging effects, are accurate and reliable.

##### 5.2. Transverse Dynamic Stiffness for a Single Degree of Freedom

The same conclusion can also be drawn by comparing the transverse dynamic stiffness of the cable for a single degree of freedom using different methods. Figure 8 shows the variations of the static stiffness along the cable length for a single degree of freedom using three methods. Figures 9(a)–9(d) show the variations of the coefficients of the dynamic stiffness along the cable length, based on the estimated value of the first-, second-, third-, and fourth-order modal frequencies, respectively (with a reference static stiffness value obtained via a finite element method). As can be seen from Figures 8 and 9, the analytical solution of the dynamic stiffness is accurate when considering sagging effects. The analytical solution when sagging effects are neglected is also accurate for calculating static stiffness, but a large error occurs when using this solution to calculate the first-order dynamic stiffness. The static stiffness obtained by Guyan’s static condensation using a finite element method gives results that are very close to those of the former two methods. The accuracy of the Guyan approach is acceptable for calculating the dynamic stiffness at fundamental frequencies. However, as the order increases, the results of the Guyan approach become quite different from those obtained from the two methods mentioned previously.

This shows that static condensation method using finite element analysis can be appropriately used as a standard to judge the accuracy of the calculated static stiffness, but not as the standard for dynamic stiffness. In the existing dynamic condensation method using finite element analysis, the stiffness matrix of the model and the mass matrix are reduced to a specified degree of freedom at a certain frequency point (i.e., the modal frequency). The question of whether it is suitable to be considered as a reference for the dynamic stiffness matrix at any frequency still needs to be explored. A new topic is thus proposed, namely, how to realize the reduction of dynamic stiffness at any frequency, a result which would play an important role in the research of the universal dynamic stiffness of structures.

##### 5.3. Vibration Tests of Full-Scale Cables

To further verify the accuracy of the analytical algorithm for the transverse dynamic stiffness matrix provided in this paper, on-site excitation data for a full-scale cable in a cable-stayed bridge are used.

The stay-cable used in the test is 267.3 m in length and 9.84 cm in diameter, with an elastic modulus of 2e11 Pa and a mass per unit length of 76.4 kg/m. A horizontal tension is applied, and the cable is fixed in the trough, with both ends clamped and secured. A TST-1000 electromagnetic vibration exciter is used at the midpoint of the cable, with power range of 0–1000 N and a maximum stroke of ±12.5 mm. In the experiment, a KD9200 resistance displacement meter is used to measure the displacement at the excitation point of the cable. Figure 10 shows the layout and plan of the on-site excitation test.

To begin, the first-order natural frequency is estimated to be using the formula for a taut string. The frequency obtained is multiplied by 0.8, 0.9, 1, 1.1, and 1.2, respectively, to be used as five different work conditions in excitation. Figure 11(a) is the time history of the displacement of the excitation point at the midpoint of cable, where the data is acquired by the on-site acquisition device, and the time history of the excitation force is recorded by the excitation device. Figure 11(b) is the curve composed of the maximum amplitude points of stable segments selected from the time histories of displacement and force. The average values are taken to estimate dynamic stiffness.

Equations (42) and (48) are used to calculate the dynamic stiffness matrix of the cable considering (method 1) and neglecting (method 2) sagging effects. Equation (45) is used to calculate the total stiffness of each matrix. Then, (49) is used to obtain the transverse dynamic stiffness with a single degree of freedom, which corresponds to the excitation point of the test. The calculation results are given in Table 3. The data in the table show, due to the consideration of sagging effects, that method 1 has better precision and yields a result that is closer to the observed value. The maximum error is less than 4%. The method that neglects sagging effects and the finite element method based on static condensation both have larger errors. This result proves that the analytical algorithm of the transverse dynamic stiffness of a cable given in this paper has better precision and reliability.

#### 6. Conclusion

The transverse vibration of cables is a primary research focus in the field of cable dynamics. The transverse force element is the most common element in cable structures. The transverse dynamic stiffness of cables is a bridge between transverse vibration and transverse force (excitation). Thus, it is of great significance for the vibration analysis of cable structures and the vibration control of an entire project.

This paper has proposed a dynamic stiffness matrix with two degrees of freedom at any point on a cable (transverse displacement and cross section rotation). On this basis, the dynamic stiffness corresponding to the transverse displacement is derived. Though several factors are simultaneously taken into account, including sagging effects, inclination, flexural stiffness, and boundary conditions, the matrix has a simple form and can be easily realized in programming applications. The existing studies have provided algorithms of dynamic stiffness at the ends of cables, which is incomparable to the results in this study. Compared with the finite element method based on static condensation, our method not only achieves a reasonable and accurate result for the static stiffness but also precisely calculates the dynamic stiffness. The variation pattern of the obtained dynamic stiffness is consistent with the findings of the existing studies. Therefore, the proposed method can be used as an effective tool in studying the transverse vibration of cables.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This paper was sponsored by the Project of National Key Technology R&D Program in the 12th Five Year Plan of China (Grant no. 2012BAJ11B01), the National Nature Science Foundation of China (Grant no. 50978196), the Key State Laboratories Freedom Research Project (Grant no. SLDRCE09-D-01), the Fundamental Research Funds for the Central Universities, and the State Meteorological Administration Special Funds of Meteorological Industry Research (Grant no. 201306102).

#### References

- E. D. S. Caetano,
*Cable Vibrations in Cable-Stayed Bridges*, vol. 9, IABSE, Zürich, Switzerland, 2007. - N. J. Gimsing and C. T. Georgakis,
*Cable Supported Bridges: Concept and Design*, John Wiley & Sons, New York, NY, USA, 2011. - E. Rasmussen, “Dampers hold sway,”
*Civil Engineering—ASCE*, vol. 67, no. 3, pp. 40–43, 1997. View at Google Scholar - P. Warnitchai, Y. Fujino, and T. Susumpow, “A non-linear dynamic model for cables and its application to a cable-structure system,”
*Journal of Sound and Vibration*, vol. 187, no. 4, pp. 695–712, 1995. View at Publisher · View at Google Scholar · View at Scopus - R. W. Clough and J. Penzien,
*Dynamics of Structures*, McGraw-Hill, New York, NY, USA, 1975. - E. de sa Caetano,
*Cable Vibration in Cable-Stayed Bridges*, IABSE-AIPC-IVBH, Zurich, Switzerland, 2007. - V. Kolousek,
*Vibrations of Systems with Curved Members*, vol. 23, International Association for Bridge and Structural Engineering, Zurich, Switzerland, 1963. - A. G. DavenPort, “The wind induced vibration of guyed and self-supporting cylindrical columns,”
*Transactions of the Engineering Institute of Canada*, pp. 119–141, 1959. View at Google Scholar - A. G. DavenPort and G. N. Steels, “Dynamic behavior of massive guyed cables,”
*Journal of the Structural Division*, vol. 91, no. 2, pp. 43–70, 1965. View at Google Scholar - A. S. Veletsos and G. R. Darbre, “Free vibration of parabolic cables,”
*Journal of the Structural Engineering*, vol. 109, no. 2, pp. 503–519, 1983. View at Google Scholar · View at Scopus - A. S. Veletsos and G. R. Darbre, “Dynamic stiffness of parabolic cables,”
*Earthquake Engineering & Structural Dynamics*, vol. 11, no. 3, pp. 367–401, 1983. View at Google Scholar · View at Scopus - U. Starossek, “Dynamic stiffness matrix of sagging cable,”
*Journal of Engineering Mechanics*, vol. 117, no. 12, pp. 2815–2829, 1991. View at Google Scholar · View at Scopus - U. Starossek, “Cable dynamics—a review,”
*Structural Engineering International*, vol. 4, no. 3, pp. 171–176, 1994. View at Publisher · View at Google Scholar - A. Sarkar and C. S. Manohar, “Dynamic stiffness matrix of a general cable element,”
*Archive of Applied Mechanics*, vol. 66, no. 5, pp. 315–325, 1996. View at Publisher · View at Google Scholar - J. Kim and S. P. Chang, “Dynamic stiffness matrix of an inclined cable,”
*Engineering Structures*, vol. 23, no. 12, pp. 1614–1621, 2001. View at Publisher · View at Google Scholar · View at Scopus - W. Lacarbonara, “The nonlinear theory of curved beams and fexurally stiff cables,” in
*Nonlinear Structural Mechanics*, pp. 433–496, Springer, New York, NY, USA, 2013. View at Publisher · View at Google Scholar - A. Luongo and D. Zulli, “Dynamic instability of inclined cables under combined wind flow and support motion,”
*Nonlinear Dynamics*, vol. 67, no. 1, pp. 71–87, 2012. View at Publisher · View at Google Scholar · View at Scopus - J. A. Main and N. P. Jones, “Vibration of tensioned beams with intermediate damper. I: formulation, influence of damper location,”
*Journal of Engineering Mechanics*, vol. 133, no. 4, pp. 369–378, 2007. View at Publisher · View at Google Scholar · View at Scopus - J. A. Main and N. P. Jones, “Vibration of tensioned beams with intermediate damper. II: damper near a support,”
*Journal of Engineering Mechanics*, vol. 133, no. 4, pp. 379–388, 2007. View at Publisher · View at Google Scholar · View at Scopus - N. Hoang and Y. Fujino, “Analytical study on bending effects in a stay cable with a damper,”
*Journal of Engineering Mechanics*, vol. 133, no. 11, pp. 1241–1246, 2007. View at Publisher · View at Google Scholar · View at Scopus - Y. Fujino and N. Hoang, “Design formulas for damping of a stay cable with a damper,”
*Journal of Structural Engineering*, vol. 134, no. 2, pp. 269–278, 2008. View at Publisher · View at Google Scholar · View at Scopus - G. Ricciardi and F. Saitta, “A continuous vibration analysis model for cables with sag and bending stiffness,”
*Engineering Structures*, vol. 30, no. 5, pp. 1459–1472, 2008. View at Publisher · View at Google Scholar · View at Scopus - H. M. Irvine, “Free vibrations of inclined cables,”
*Journal of the Structural Division*, vol. 104, no. 2, pp. 343–347, 1978. View at Google Scholar