Abstract

The nonlinear dynamic behavior of rain-wind-induced vibration of inclined cable is investigated with the consideration of the equilibrium position of the moving rivulet. The partial differential governing equations of three-degree-of-freedom on the model of rain-wind-induced cable vibration are established, which are proposed for describing the nonlinear interactions among the in-plane, out-of-plane vibration of the cable and the oscillation of the moving rivulet. The Galerkin method is applied to discretize the partial differential governing equations. The approximately analytic solution is obtained by using the method of averaging. The unique correspondence between the wind and the equilibrium position of the rivulet is ascertained. The presence of rivulet at certain positions on the surface of cable is then proved to be one of the trigger for wind-rain-induced cable vibration. The nonlinear dynamic phenomena of the inclined cable subjected to wind and rain turbulence are then studied by varying the parameters including mean wind velocity, Coulomb damping force, damping ratio, the span length, and the initial tension of the inclined cable on the model. The jump phenomenon is also observed which occurs when there are multiple solutions in the system.

1. Introduction

Dynamic behavior of corresponding differential systems has been of extensive concern and investigated by many scholars [18]. One of the systems is the inclined cable, which is a structure characterized by small stiffness, low mass, and light damping. The similar structures including transmission line, suspended cable, and submarine cable are widely used in the long-span structures, so the dynamic study of the inclined cable is of great engineering significance. Rain-wind-induced vibration of inclined cable is a phenomenon of large amplitude vibration with low frequency under the conditions of rain and wind. Many wind tunnel tests have been carried on a cable model. Large amplitude vibrations of cables were first observed by Hikami and Shiraishi [9]. Upper rivulet was regarded as the origin of large amplitude vibration due to its formulation on the cable, causing the cable cross section aerodynamically unstable. Yamaguchi [10] took the lead in aerodynamic test of cable with artificial rivulet. Aerodynamic forces acting on the cable were obtained. Matsumoto et al. [11] reported that another factor named an axial flow generated at the near wake of the inclined cable may trigger the rain-wind-induced oscillations. He stated that rain-wind-induced vibration can be explained as a vortex-induced vibration, which occurs at limited high reduced wind velocity region [12]. Bosdogianni and Olivari [13] compared wind tunnel results of cables with moving rivulet and fixed rivulet. They concluded that it was the presence of rivulet at certain positions on the surface of cable and not the motion of rivulet that caused cable instability. Then the aerodynamic forces acting on the moving rivulet were first measured by Gu et al. [14], proposing a new description of friction force between the water rivulet and the cable surface.

On the basis of the measured aerodynamic coefficients, the kinetic equations of rain-wind-induced vibration have been established. However, before the measurement of the aerodynamic coefficients on the moving rivulet, the motion of rivulet was assumed to be harmonic [1517]. With the acquisition of the aerodynamic coefficients on the moving rivulet, the analytic model for cable with the kinetic equation of the rivulet was built [18]. Zhang et al. [3] presented a two-degree-of-freedom model of the rain-wind-induced vibration of a continuous stay cable. The phenomenon of bifurcation was observed in the system using the singularity theory. Zhang et al. [19] established the three-degree-of-freedom model of rain-wind-induced vibration of inclined cable. The motion characteristics of rivulet were achieved with two rivulet models through numerical calculation and some more rational explanations of the phenomenon were obtained.

The influences of equilibrium position of rivulet on cable have been studied, in which the equilibrium position of rivulet was always assumed to be independent of wind. The equilibrium position of rivulet was altered factitiously to inspect its effects on the vibration while the mean wind velocity was fixed [15, 18]. However, the equilibrium position of rivulet is uniquely corresponding to the wind once the mass of the rivulet is fixed. Meanwhile, the similar structures including beam, suspended cables, and inclined cable without wind or rain have been extensively researched analytically to discover more dynamic phenomena [48]. However, the nonlinear dynamic behavior of rain-wind-induced vibration has been rarely investigated using approximately analytic method for its complexity.

In comparison with the literature above, the equilibrium position of the moving rivulet is considered and the nonlinear dynamic characteristics of rain-wind-induced vibration are investigated analytically by using the approximate method of averaging. The present paper established the continuum model of the inclined cable subjected to wind and rain turbulence with the consideration of the equilibrium position of the rivulet. The reliance of the equilibrium position of the rivulet on wind is figured out and the effects of the equilibrium position of the rivulet on the cable vibration are observed and proved. Furthermore, the amplitude response curves of parameters including mean wind velocity, Coulomb damping force, damping ratio, the span length, and the initial tension of the inclined cable present abundant dynamic behaviors. The jump phenomenon is also observed when multivalued solutions exit in the system.

2. Model and Nonlinear Equations

Assumptions of inclined cable are made in present paper: (1) The flexural rigidity, torsional stiffness and shear stiffness are ignored; (2) The constitutive relation of cable deformation submits to Hooke's law and the points bear the stress evenly; (3) The axial motion of the cable is ignored; (4) The influences of bridge and tower are disregarded.

The analytic model of rain-wind-induced vibration of inclined cable is shown in Figures 1(a), 1(b), and 1(c), is mean wind velocity, is cable angle, is wind angle, is the span length of the cable, and and are dynamic displacement in the and direction, respectively. is deflection at the equilibrium position of cable. Define that the vibration in the plane - is in-plane vibration, while the vibration in the plane - is out-of-plane vibration. is the mass per unit length of cable, is the initial tension of the cable, and is the dynamic tension of the cable. and are arc length of undeformed and deformed cable, respectively.

Figure 1(d) presents the force analysis of rivulet. is the equilibrium position of rivulet, is dynamic displacement of rivulet which is the angle deviating from the equilibrium position, is dynamic displacement of rivulet which is the angle deviating from direction. is the mass per unit length of rivulet, is the gravity acceleration, is the radius of cable, is Coulomb damping force whose direction is opposite to its motion, is the linear viscous damping coefficient between cable and rivulet, and effective angle of attack is the angle between effective wind velocity and direction. is the angle between relative wind velocity and z direction.

Balance equations of the system are obtained from the force analysis of the cable element and rivulet (see Figures 1(c) and 1(d)): where is the elasticity modulus of cable, is the cross sectional area of cable, and is dynamic strain. , are the damping coefficients of in-plane and out-of-plane vibration, respectively. , , and are the aerodynamic forces on the in-plane, out-of-plane, and rivulet vibration, respectively, with the function of , , , and which are defined by where is air density, is the diameter of the cable, is the characteristic length of rivulet, and is equivalent wind attack angle and can be written as

In (3), lift coefficient , resistance coefficient of cable are the fitting functions of , aerodynamic coefficients of rivulet and are the fitting functions of which are obtained in the experiment [18] and can be expressed as where are constant coefficients, , . At static situation, the force analysis of rivulet is given by The equilibrium position can be attained from (6) and expressed as . The equation of rivulet on the equilibrium position is obtained by substituting (Figure 1(d)) into system (1). Galerkin method is then adopted to discretize system (1) and the first mode of vibration is reserved, assuming that the final equations are obtained: where is the natural frequency of in-plane, out-of-plane, and rivulet vibration, respectively. , are integral constants related to physical parameters of cable and mean wind velocity, . The equations are composed of geometric nonlinearity and aerodynamic nonlinearity.

3. Average Analysis

Nondimensionalize system (8) and assume that , , , , , and ; then system (8) becomes , , and are integral constants. The analytic solutions are assumed to be is amplitude and is phase position. They can be obtained according to the orthogonality of derived equations and conjugate equations which are given by is the natural frequency of derived equations. , , , , , and . The basic solutions , , , and of derived equations and conjugate equations are obtained according to system (10), , .

Given that the rivulet vibrates mainly with the frequency of cable, which is in system (8) is related to , described as making can ascertain the interval of the symbolic function. Thus , . The last average equations can be written as

Make then is arbitrary constant; the steady state solutions are obtained by substituting (18) into (12)

4. Results Analysis

Parameters are selected as follows [20]: , , , , , , , , , , , , , and is 10% of rivulet’s mass per unit length.

4.1. Numerical Solutions

Figure 2 shows the numerical solutions of inclined cable vibration calculated form system (8) utilizing Matlab. The full line in Figure 2(a) is the time-history of in-plane vibration while the dotted line represents the out-of-plane vibration of the cable. The in-plane amplitude is about 4 times as that of out-of-plane vibration and they have the same phase. The time-history of the displacement response of the rivulet is presented in Figure 2(b). The moving rivulet vibrates at its equilibrium position. Note that, the vibration is not exhibited to be harmonic and there is not one frequency component. Figure 3 provides the vibration frequencies of in-plane and rivulet vibration. Cable vibrates at its natural frequency which is equal to 1.25 Hz as shown in the Figure 3(a). However, the frequency components of rivulet vibration are more complex. As shown in Figure 3(b), its primary frequency is the same as the natural frequency of cable; in addition, some higher order frequencies also exist in its vibration.

4.2. The Comparison between Analytic and Numerical Solutions

Figure 4 represents the analytic solutions together with the numerical results. The -coordinate is the mean wind velocity and the -coordinate is the maximum displacement response amplitude of the in-plane, out-of-plane, and rivulet vibration, respectively. Analytic solution is a smooth curve, describing the trend of the amplitude more directly and visually. Numerical solutions are several discrete points. The amplitudes of in-plane, out-of-plane, and rivulet vibration present a directly proportional relationship. The response curves also confirm that rain-wind-induced inclined cable vibration is a kind of restricted amplitude vibration rather than galloping. The analytic results are in good agreement with the numerical results. Thus, the validity of the approximately analytic method is proved.

4.3. Analysis of Equilibrium Position of Rivulet

The correspondences among the equilibrium position of rivulet, mean wind velocity, and rivulet mass are calculated from (6) and shown in Figure 5(a). The equilibrium position of rivulet is determined by mean wind velocity with certain rivulet mass. The effective wind velocity is around 1 m/s~24 m/s. Greater rivulet mass needs larger mean wind velocity to form on the surface of the inclined cable. Keeping the other parameters invariant, the effects of rivulet mass on the cable vibration subjected to variational mean wind velocity are presented in Figure 5(b). Greater rivulet mass on the surface of inclined cable will cause stronger oscillation. The impact of the equilibrium position of rivulet on the inclined cable is investigated and plotted in Figure 5(c). It can be perceived that the amplitude is up to the maximum when the equilibrium position of rivulet is around 48° and almost irrelevant to the rivulet mass. It confirms that the presence of rivulet at certain positions on the surface of inclined cable is one of the trigger that stimulate the inclined cable to vibrate.

4.4. Nonlinear Dynamic Analysis

The nonlinear dynamic behavior of cable reflected by parameters is given in Figure 6. The effect of mean wind velocity on the cable is presented in Figure 6(a), in which the full line is stable solution and the dotted line represents unstable solution. The figure illustrates that the mean wind velocity modify the number of the solutions. Two stable solutions appear at the mean wind velocity range from 7 m/s to 8.7 m/s. Figure 6(b) describes the nonlinear response curve of Coulomb damping force acting on the cable amplitude subjected to various mean wind velocities. When mean wind velocity is selected as 8 m/s, two stable solutions occur for is less than 15%. As increases, the amplitude of inclined cable first arises along the smaller branch of the stable solution and it will not jump until reaches 15%. Once increases to 15%, the smaller branch of amplitude will be substituted by the larger branch. The jump phenomena will also exist while is about 9% for is 7 m/s. The variance of will make no change to the number of solution while the mean wind velocity is 6 m/s, where the system has only one stable solution and no jump phenomena occurs. Similarly, the response curve of damping ratio acting on the inclined cable depicts the nonlinear behavior of the system (Figure 6(c)). Multiple solutions lie in the certain range of when is beyond its critical mean wind velocity 7 m/s; otherwise, unique solution is obtained irrelevant to the damping ratio.

The relation curves of amplitude and the span length of cable are obtained in Figure 6(d). It also illustrates that the mean wind velocity modified the stability of the inclined cable. When the mean wind velocity is selected as 7 m/s or 8 m/s, the amplitude will jump to a solution of higher branch as the span length of lengthens. However, the mean wind velocity of 6 m/s will not cause the modification of the number of the solutions but only a slight increase. Figure 6(e) shows the response curve of cable amplitude and initial tension . Similarly, the mean wind velocity of 7 m/s and 8 m/s will cause cable amplitude to increase drastically while the lower wind velocity acting on the cable makes little sense to the variance of the initial tension.

5. Conclusions

The continuum model of inclined cable subjected to wind with moving rivulet on its surface is established considering the equilibrium position of the rivulet. The Galerkin and average methods are adopted to analyze the cable system analytically. The correspondence between mean wind velocity and the equilibrium position of the rivulet is ascertained considered the different rivulet mass. The equilibrium position of the rivulet is proved to be one of the major roles in determining the maximum amplitude of the oscillating cable.

The nonlinear dynamic behavior of the rain-wind-induced cable vibration is then investigated considering the continuous change of parameters including mean wind velocity, Coulomb damping force, damping ratio, the span length, and the initial tension of the inclined cable. Mean wind velocity is observed to be a significant factor in jump phenomena, which is discovered in the range of multiple solutions.

Acknowledgments

The authors wish to express their special gratitude to the reviewers and the editor for their valuable comments and suggestions that led to a truly significant improvement of the paper. This work was supported by the National Natural Science Foundation of China (no. 51009107), Tianjin Natural Science Foundation (no. 13JCQNJC0420 -0), and Key Project of Tianjin Natural Science Foundation (no. 13JCZDJC27100).