Nonlinear dynamical analysis for the cable excited with parametric and forced excitation.
Along with the new materials, new technology continues to improve, large span cable is applied widespread, and the cable-stayed bridge becomes one of the main forms of long-span bridges in the word. Because of the flexibility, quality, and small damping characteristic, the dynamical problem of the cable structure is relatively prominent. In recent decades, the dynamics of stayed cable has gained widespread concern in the bridge engineering. Researches have shown that the breeze and light rain conditions may induce the cable having dramatic large amplitude vibrations [1-4]. Such sharp vibration may cause the fatigue fracture of cable or junction of fatigue cracks in the cable anchor, undermine the claim of anticorrosion system, and cause harm to the safety and durability of the bridge.
In the earlier 1970s, some researchers began to study the cable vibration in cable-stayed bridge . Variety induced vibration theory is proposed, such as vibration induced with rain and wind, secondary axial flow vibration, three-dimensional vortex shedding, and wake galloping. But so far stayed cable vibration mechanism was not studied clearly. It is shown that the cable vibration is a complex phenomenon affected by multiple factors. The aerodynamic characteristics caused by the wind and rain and the geometric nonlinearity of the bridge structure. The parametric vibration about cable is studied in the literature , in which a single mass is used to simulate continue cable. The parametric vibration about cable is also studied in literature , in which the cable is considered a vertical chord. Considering the initial deflection, the tension, and external excitation, the coupling dynamics of the cable about parametric vibration and forced vibration is studied .
The vertical direction vibration of the deck will be induced with the wind loads or other live loads role. When an order of the natural frequency of the deck is similar to the twice the cable's natural frequency, slight vibration of the bridge deck will inspire cable vibration sharply [9, 10]. This phenomenon is significant feature of the parametric vibration. Therefore, studying this vibration problem is meaningful to design and protection of the cable structure. In the literature , the equations governing the linearized dynamics of the cable-stayed beam are presented and the main features of the parametric vibration phenomenon are investigated. The literature  investigates the coupled nonlinear vibration behaviors of the cable/deck system by numerical analysis of a proposed three-degree-of-freedom model and using the Runge-Kutta integration method reveal some fundamental characteristics of the nonlinear cable vibration.
Parametric vibration system is the incentive depending on the time and appears as a parameter in the control equation. A feature of such a vibration system is the system parameter with the time change. When the natural frequency of the system has a multiple relationship with the excitation frequency, even a small excitation can cause complex responses [13, 14] of the system which is different from the external excitation.
Based on the literature , using bifurcation theory and nonlinear vibration theory, the coupling effect on parametric excitation and external force to the cable vibration is studied. Using numerical method, the bifurcation response is obtained, in which some nonlinear phenomenon is shown with figures.
2. Nonlinear Dynamical Equation of the Stay Cable
In order to study vibration mechanism of the stay cable bridge, the cable vibration system is modeled as Figure 1. In this mechanics model, the initial deflection and off-plane motion are not considered. Using nonlinear vibration theory , the dynamical equation of the cable in plan xy is expressed as (1). Consider
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (1)
With the boundary condition
w(0,t) = 0, w(L, t) = D(t) sin y, (2)
in which w denotes the displacement of cable in y direction, L denotes the length, A is the section area of the cable, E denotes elastic modulus, [rho] denotes the density of the material, and [gamma] is the angle between cable and vertical direction. D(t) = D cos [omega]t is the vertical displacement of the connection point of the cable and the bridge deck. If only the first order modal is considered, the solution of (1) can be denoted by
w(x, t) = [x/L]D(t)sin + y(t)sin[[pi]x/L]. (3)
For purpose of analysis, some nondimensional parameters are introduced as follows:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4)
Then (1)-(3) can be rewritten as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (5)
[??](0, [tau]) = 0, [??](1, [tau]) = d cos [tau] sin [gamma], (6)
[??]([eta], [tau]) = [eta]d cos [tau] sin [gamma] + y([tau])sin [pi][eta]. (7)
Substituting (7) into (5) and adding the viscous damping item, using Galerkin integration method, the differential equation about y(t) can be simplified as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (8)
when [gamma] = 0 and (8) has not forced vibration item, which means the system (1) has only parametric vibration and this case has been introduced in literature . In general case, [gamma] [not equal to] 0, forced vibration and parametric vibration items exist in (8) at the same time and system (1) should be considered a dynamical system with combined forced and parametric excitation. Considering [d.sup.2] is higher order small quantities than d, (1/2)[d.sup.2][cos.sup.2][tau] [sin.sup.2] [gamma] can be ignored and (8) can be simplified as
[d.sup.2]y/d[[tau].sup.2] + 2[xi][dy/d[tau]] + [[[lambda].sup.2] + [c.sup.2.sub.2]d cos [tau] cos [gamma]]y([tau])
+ [[c.sup.2.sub.2]/2][y.sup.3](t) - [2/[pi]]d sin [gamma] cos [tau] = 0. (9)
3. Bifurcation Response Equations
In order to take advantage of multiscale method analysis, ignoring the effect of deflect and using the substitution parameters,
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (10)
Then, (9) can be rewritten as
[d.sup.2]y/d[[tau].sup.2] + 2[epsilon][mu][dy/d[tau]] + [[lambda].sup.2]y + [[epsilon].sup.2][beta][y.sup.3] + [epsilon]fy cos [tau] + p cos [tau] = 0. (11)
Based on the multiple scale method, the perturbation solution of (11) can be expensed as the power series about small parameter [epsilon]:
y([epsilon], [tau]) = [y.sub.0]([T.sub.0], [T.sub.1], [T.sub.2]) + [epsilon][y.sub.1]([T.sub.0], [T.sub.1], [T.sub.2])
+ [[epsilon].sup.2][y.sub.2]([T.sub.0], [T.sub.1], [T.sub.2]) + ..., (12)
in which [T.sub.i] = [[epsilon].sup.i][tau]. Considering the differential connections,
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (13)
Substituting (12) into 11), a series perturbation equations can be obtained as follows:
[D.sup.2.sub.0][y.sub.0] + [[lambda].sup.2][y.sub.0] = -p cos [T.sub.0], (14)
[D.sup.2.sub.0][y.sub.1] + [[lambda].sup.2][y.sub.1] = -2[D.sub.0][D.sub.1][y.sub.0] - 2[mu][D.sub.0][y.sub.0] - [fy.sub.0]cos [T.sub.0]. (15)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (16)
Solving (14), y0 can be expressed as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (17)
in which G = p/2(1 - [[lambda].sup.2]). Substituting 17) into 15), then (15) can be rewritten as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (18)
Considering the 1/2 main parametric resonance condition, assuming
[lambda] = [[omega].sub.0]/[omega] = 1/2 = [epsilon][sigma], (19)
in which [sigma] is tuning parameter, substituting (19) into (18), and eliminating the long-term items, then (20) is obtained as follows:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (20)
Solving (18), [y.sub.1] can be obtained as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (21)
Substituting (17), (21) into 16) and eliminating the long-term items, then (22) is obtained as follows:
-2i[lambda][D.sub.2]A - [D.sup.2.sub.1]A - 2[mu][D.sub.1]A - [beta](3[A.sup.2][bar.A] + 6[G.sup.2]A)
- [f.sup.2]A/4(2[lambda] + 1) = 0. (22)
Associating (20), eliminating the [D.sup.2.sub.1]A and [D.sub.1]A items from (22), and considering 2[lambda] [approximately equal to] 1, then (23) can be obtained as follows:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (23)
Considering the differential connection,
dA/d[tau] = [epsilon][D.sub.1]A + [[epsilon].sup.2][D.sub.2]A. (24)
Equation (24) can be rewritten as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (25)
Assuming that A = (1/2)a([tau])[e.sup.i[phi]([tau])] and is substituted into (25), separating the real part and the imaginary part, (26) can be obtained as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (26)
[theta] = [epsilon][sigma][tau] - 2[phi]. (27)
For da/d[tau] = d[theta]/d[tau] = 0 corresponding to the steady-state solution of the system, (26) can be rewritten as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (28)
Eliminating [theta] from (28), then
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (29)
In addition to ordinary solution a = 0, the nonordinary solution of (29) meets the following equation:
[[mu].sup.2] + [([1/2][sigma] - [epsilon][[mu].sup.2] + 6[epsilon][beta][G.sup.2] + 3[epsilon][beta][a.sup.2]).sup.2] = [([1/2]f - f[sigma][epsilon]).sup.2]. (30)
Solving (30), the nonordinary solution can be obtained as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (31)
Considering different parameters, there may exist some possible solutions as follows:
(I) [((1/2)/ - f[sigma][epsilon]).sup.2] - [[mu].sup.2] = 0, (1/2)[sigma] - [epsilon][[mu].sup.2] + 6[epsilon][beta][G.sup.2] < 0, and system has one nonzero solution;
(II) [((1/2)f - f[sigma][epsilon]).sup.2] - [[mu].sup.2] > 0, (1/2)[sigma] - [epsilon][[mu].sup.2] + 6[epsilon][beta][G.sup.2] < 0, [((1/2)[sigma] - [epsilon][[mu].sup.2] + 6[epsilon][beta][G.sup.2]).sup.2] > [((1/2)f - f[sigma][epsilon]).sup.2] - [[mu].sup.2], and system has only two nonzero real solutions;
(III) [((1/2)f - f[sigma][epsilon]).sup.2] - [[mu].sup.2] > 0, (1/2)[sigma] - [epsilon][[mu].sup.2] + 6[epsilon][beta][G.sup.2] < 0, [((1/2)[sigma] - [epsilon][[mu].sup.2] + 6[epsilon][beta][G.sup.2]).sup.2] < [((1/2)f - f[sigma][epsilon]).sup.2] - [[mu].sup.2], and system has only one nonzero real solution;
(IV) [((1/2)f - f[sigma][epsilon]).sup.2] - [[mu].sup.2] > 0, (1/2)[sigma] - [epsilon][[mu].sup.2] + 6[epsilon][beta][G.sup.2] > 0, [((1/2)[sigma] - [epsilon][[mu].sup.2] + 6[epsilon][beta][G.sup.2]).sup.2] > [((1/2)f - f[sigma][epsilon]).sup.2] - [[mu].sup.2], and system has only one nonzero real solution;
(V) other conditions, system exist only one zero solution.
From the above-mentioned discussion, it is sure that there are critical excitation amplitude [f.sub.crit] and fold bifurcation point [f.sub.T]. In order to avoid saddle-node bifurcation, the critical bifurcation point [f.sub.crit] is important. Consider
[f.sub.crit] = [+ or -] [mu]/[epsilon][omega]. (32)
At this time, the system parameters meet with
[([1/2][sigma] - [epsilon][[mu].sup.2] + 6[epsilon][beta][G.sup.2]).sup.2] = [([1/2][f.sub.T] - [f.sub.T][sigma][epsilon]).sup.2] - [[mu].sup.2]. (33)
The supercritical bifurcation is decided with bifurcation point number [f.sub.T] in (33).
4. Numerical Analysis for the Cable Vibration
In this section, the numerical method and software MATHMATIC are used to solve a cable system. It aims at analyzing the vibration frequency of deck and inclination angle of the cable effect on the cable vibration. The parameters of the structure are defined as Table 1, in which the damping is not considered and the vibration amplitude of the anchorage point is defined as D = 5 mm. the first natural vibration frequency is obtained f = 2.6 Hz and [[omega].sub.0] = 16.352.
Using the numerical method, the bifurcation equation (29) and the different equation (9) are solved, in which the deck vibration amplitude is assumed as constant. The parametric resonance bifurcation curves corresponding to various inclination angles of the cables can be obtained as Figure 2.
Results of bifurcation points of parametric resonance corresponding to various cables' inclination angles are given in Figure 2. It shows that inclination angles change the dislocation of bifurcation points.
When the inclination angles of the cable are made small, two bifurcation points become closer and the amplitude of resonance becomes larger. The parametric resonance becomes a dominant factor compared to the forced oscillation. In order to observe the effect of the inclination angle of amplitude response, the amplitude-inclination angle curves are given in Figure 3; in those curves, the vibration frequencies of the anchor point are fixed.
From Figure 3, it can be found that tune frequency parameter a affects not only the bifurcation point values but also the numbers of bifurcation points. Based on the numerical analyzing, the rules are found as follows. (1) When [sigma] [less than or equal to] 0, the system has only one critical bifurcation point and the system has only one steady state response solution. (2) The nonsteady response solution occurred as a increase, and a fold bifurcation occurred. (3) When a reach [sigma] fixed value, the systems' response changed to be more complicated.
In order to observe the effect of parametric excitation, the curves time domain and phase plan result of some systems are given in Figures 4 and 5.
From Figures 4 and 5, some interesting phenomenon is observed. As the case of 0 < [gamma] < [pi]/2, the vibration responses are caused by parametric excitation and forced excitation or the corresponding joint action. As y decreases, the parametric excitation plays a major role; this conclusion meets well with . When the value of y is larger, parametric excitation and forced excitation combined effects are obvious; its dynamic behavior is more complex and will need further research. So, this paper will not repeat them.
In this paper, nonlinear dynamical characters of the cable in cable-stayed bridge are studied. First, using nonlinear structure dynamics, a nonlinear model of a cable in cable-stayed bridge excited with deck vibration is proposed. Secondly, using multiscale method, the 1/2 principle parametric resonance is studied and the bifurcation equation is obtained. And then the bifurcation characters of the dynamical system are studied. At last, the effect rules of system parameters to the dynamic behavior of the system are studied. Through study of this dynamical system, some useful conclusions are obtained.
(1) Parametric vibration of the cables in the cable-stayed bridge is induced with the vertical vibration of the deck. At the same time, the vertical vibration of the deck has the forced excitation effect on the cables.
(2) In general, the vibration response of the cable is joint effect of the parametric and forced excitation.
(3) As y decreases, the parametric excitation plays a major role. When the value of y is larger, parametric excitation and forced excitation combined effects are obvious.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This work is supported by the National Science Foundation of China (Grants nos. 51108047, 11272270, and 11102028), Science and Technology Project in Xiamen (Grant no. 3502Z20143028), and XMUT Sino-Foreign Cooperation Foundation (Grant no. E201301500).
 M. Matsumoto, N. Shiraishi, and H. Shirato, "Rain-wind induced vibration of cables of cable-stayed bridges," Journal of Wind Engineering & Industrial Aerodynamics, vol. 43, no. 1-3, pp. 2011-2022, 1992.
 Y. Hikami and N. Shiraishi, "Rain-wind induced vibrations of cables stayed bridges," Journal of Wind Engineering and Industrial Aerodynamics, vol. 29, no. 1-3, pp. 409-418, 1988.
 J. L. Lilien and A. P. Da Costa, "Vibration amplitudes caused by parametric excitation of cable stayed structures," Journal of Sound and Vibration, vol. 174, no. 1, pp. 69-90, 1994.
 H. Yamaguchi and Y. Fujino, "Stayed cable dynamics and its vibration control," in Proceedings of the International Symposium on Advances in Bridge Aerodynamics, pp. 235-253, Balkema Rotterdam, Rotterdam, The Netherlands, 1998.
 G. Tagata, "Harmonically forced, finite amplitude vibration of a string," Journal of Sound and Vibration, vol. 51, no. 4, pp. 483-492, 1977.
 Y. Fujino, P. Warnitchai, and B. M. Pacheco, "Experimental and analytical study on internal resonances in a cable-stayed beam model," Proceedings of Japan Society of Civil Engineering, vol. 432, pp. 109-118, 1991.
 Y. Fujino, P. Warnitchai, and B. M. Pacheco, "An experimental and analytical study of autoparametric resonance in a 3DOF model of cable-stayed-beam," Nonlinear Dynamics, vol. 4, no. 2, pp. 111-138, 1993.
 Z. Wang and B. Sun, "Parametric vibration of the cable in cable--stayed bridge," China Civil Engineering Journal, vol. 35, no. 5, pp. 14-22, 2002.
 A. Pinto Da Costa, J. A. C. Martins, F. Branco, and J. L. Lilien, "Oscillations of bridge stay cables induced by periodic motions of deck and/or towers," Journal of Engineering Mechanics, vol. 122, no. 7, pp. 613-621, 1996.
 N. C. Perkins, "Modal interactions in the non-linear response of elastic cables under parametric/external excitation," International Journal of Non-Linear Mechanics, vol. 27, no. 2, pp. 233-250, 1992.
 V Gattulli and M. Lepidi, "Localization and veering in the dynamics of cable-stayed bridges," Computers and Structures, vol. 85, no. 21-22, pp. 1661-1678, 2007
 Z. Kang, K. Xu, and Z. Luo, "A numerical study on nonlinear vibration of an inclined cable coupled with the deck in cable-stayed bridges," Journal of Vibration and Control, vol. 18, no. 3, pp. 404-416, 2012.
 A. H. Nayfeh and D. T Mook, Nonlinear Oscillations, John Wiley & Sons, New York, NY, USA, 1973.
 A. H. Nayfeh, Perturbation Method, John Wiley & Sons, New York, NY, USA, 1973.
C. Z. Qian, C. P. Chen, and G. W. Zhou
College of Civil Engineering and Architecture, Xiamen University of Technology, Xiamen, Fujian 361024, China
Correspondence should be addressed to C. Z. Qian; firstname.lastname@example.org
Received 10 April 2014; Revised 20 August 2014; Accepted 7 September 2014; Published 19 October 2014
Academic Editor: Ali R. Ashrafi
TABLE 1: The essential parameters of the cable. Mass per unit length ([rho]A) 0.391 kg/m Length of cable (L) 11.90 m Axial tensile force (T) 1.5 kN Section area (A) 7.28 [mm.sup.2] Elastic modulus (_E) 8.242 x 104 MPa Extensional rigidity (EA) 6.0 x 105 N
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Qian, C.Z.; Chen, C.P.; Zhou, G.W.|
|Publication:||Journal of Applied Mathematics|
|Date:||Jan 1, 2014|
|Previous Article:||Multiobjective fuzzy mixed assembly line sequencing optimization model.|
|Next Article:||Geometric lattice structure of covering and its application to attribute reduction through matroids.|