The nonlinear vibrations of cable in a multispan cable-stayed bridge subjected to transverse excitation are investigated. The MECS (multielements cable system) model, where multielements per cable stay are used, is built up and used to analyze the model properties of the multispan cable-stayed bridges. Then, a simplified two-degrees-of-freedom (2-DOFs) model, where the tower or the deck is reduced to a beam, is proposed to analyze the nonlinear dynamic behaviors of the beam and cable. The results of MECS model analysis show that the main tower in the multispan cable-stayed bridge is prone to the transverse vibration, and the local vibration of cables only has a little impact on the frequency values of the global modes. The results of simplified model analysis show that the energy can be transformed between the modes of the beam and cable when the nature frequencies of them are very close. On the other hand, with the transverse excitation changing, the cable can exhibit richer quasi-periodic or chaotic motions due to the nonlinear terms caused by the coupled mode between the beam and cable.
1. Introduction
Recently, as use of multispan cable-stayed bridges increases continuously, the real local vibration features of inclined cables and coupled oscillation characteristics between the cable and the bridge is becoming a new topic in the design process. The stabilization of the central towers under extreme wind or seismic vibration is a key issue, since they cannot be anchored to an outer fixed support. One solution is to use the stabilizing cables which run from the top of the central towers to a location on the deck near the side towers, such as the Ting-Kau Bridge in Hong Kong. Another solution is to increase the stiffness of the central towers itself, and most multispan cable-stayed bridges adopt this way, such as the Maracaibo Bridge in Venezuela, the Millau Bridge in France, the Mezcala Bridge in Mexico, the Dong-Ting Lake Bridge, the Yi-Ling Bridge, and the Bin-Zhou Bridge in China.
Ni et al. [1] investigated the effect of stabilizing cables on the seismic response of Ting Kau Bridge based on a validated 3D finite-element model. He et al. [2] investigated the influences of cable local modes and stabilizing cables for the Dong-Ting Lake Bridge on seismic excitation. Liu [3] investigated the seismic performance of rigid system, floating system, and passive energy dissipation system for the Bin-Zhou Bridge under two different earthquake records. In order to evaluate the stability of the main towers of Millau Bridge, Okamoto and Nakamura [4] proposed a new hybrid high tower and investigated their static and seismic behaviors. Most of the above studies focused on the stabilization of central towers, and little attention has been devoted to investigating the effects of cable vibration on the properties of tower or deck under transverse excitation.
On the other hand, all the above researchers used the finite-element (FE) model to investigate the stabilization of towers by taking the specific project as a paradigm. From the perspective of the coupled vibration between cable and tower (or deck), Fujino et al. [5] presented a 3-DOFs analytical model to investigate the auto-parametric interaction behavior of cable and beam structure, the case which under a random excitation is firstly studied by Xia et al. [6, 7]. The results show that the horizontal motions of the cable and beam are excited due to the autoparametric nonlinear coupling behavior under some cases. Caetano et al. [8, 9] studied the dynamic interactions between the cable and deck or tower in cable-stayed bridges by physical modelling and experimental testing. Georgakis and Taylor [10, 11] presented an alternative cable-deck model to investigate the nonlinear dynamics of an inclined cable which both induced by sinusoidal and stochastic support excitations. The different cable-deck interactions of the Guadiana Bridge under environmental excitations were investigated by Caetano et al. [12], using the vibration data acquisitions and a refined finite element model. However, few studies focus on the nonlinear vibration of the cable-tower (or beam) coupled system subjected to transverse excitation.
The objective of this paper is to study the nonlinear vibrations of cable in a multispan cable-stayed bridge subjected to transverse excitation. Taking the Bin-Zhou cable-stayed bridge (BZB) as a paradigm, the MECS model, where multi-elements per cable stay are used, is firstly built up and used to investigate the effects of the cable’s vibration on the properties of the BZB. Then, the possibility of the transverse resonance between the local model (cable) and global model (tower or deck) is analyzed. Based on the results of the above studies, a simplified 2-DOFs model, where the tower or the deck is reduced to a tower (or beam), is proposed to analyze the nonlinear coupled vibration of the cable and tower (or beam). The coupled relationship is completely from the dynamic interactions between the cable and tower (or beam). After that, the equations of motion are solved by using Galerkin’s method for the spatial problems and the method of multiple time scales for temporal problems. The stability of the steady-state solution is examined. Finally, the nonlinear behavior of the cable and tower (or beam) is analyzed by using the time histories, phase portraits, and Poincare maps with the transverse excitation changing.
2. Vibration Properties of the Multispan Cable-Stayed Bridge2.1. Outline of the Bin-Zhou Cable-Stayed Bridge
The BZB, as shown in Figure 1, is a three-tower cable-stayed bridge with two main spans of 300 m and two-side spans of 84 m each. The bridge deck is separated into two carriageways with a width of 13.75 m each. There are two minor towers with heights of 75.78 m in both sides and a main tower with a height of 125.28 m in the center. It has a unique feature that there is no horizontal component connecting the towers. The carriageways consist of two longitudinal pre-stressed concrete box girders with a length of 767.1 m and a height of 3 m. Along the deck edges, there are cross-girders at 7 m and 6 m intervals for the main span and two side spans, respectively. The locations of cables are consistent with the cross-girders. The bridge deck is fixed at main towers, but linked at the minor towers and two end abutments. There are 200 main stay cables anchored to the towers with intervals of 2 m (partial 2.5 m). Shown as Figure 1, from left to right the cables anchored to the left minor tower, the main tower and the right minor tower are named as A1~A12, N1~N26 and J1~J12 respectively.
View of the Bin-Zhou Bridge (cm).
The 3-D FE model of the BZB is established in ANSYS, as shown in Figure 2. The bridge members, girders, towers, and deck are simulated by spatial beam or shell elements with six DOFs at each node, respectively. In order to reflect the influence of cable vibration, the cables are simulated by spatial spar elements with 3-DOFs at each node; each cable is divided by a length of 6 m. As a result, the FE model involves 7696 nodes and 7514 elements.
FE model of the Bin-Zhou Bridge.
2.2. Modal Analysis
Based on the aforementioned criterions, the natural vibration frequencies and global modes are examined by using OECS (one-element cable system) and MECS models. The frequencies and mode shapes of the first 10 modes of OECS model and the first 300 modes of MESC model are computed. Among the latter ones, the maximum frequency is up to 1.598 Hz, and over 85% of modes are pure cable local vibration modes. This property is similar to the Ting Kau Bridge [1, 13].
According to the same global mode shapes, Table 1 shows a comparison of the natural frequencies of the OECS and MESC models. From Table 1, the first mode of the BZB is predominantly vertical bending, which differs from most of the cable-stayed bridges. This phenomenon can be attributed to the semifloating system adopted in the BZB. The vibration frequencies of the first 10 global modes of MECS model have minor difference to the OECS model and have a max difference of 3.24% in 5th mode and −0.476% in 2nd mode. Moreover, the increase in frequencies is related to the deck and the decrease is related to the towers. These results differ from the conclusion given by Wang et al. [14]. The tower modes also have a larger proportion in the first 10 global modes of the BZB, and most of them are the lateral vibrations, as seen in Figure 3. It may be explained by the fact that the lack of the horizontal component between the towers makes the overall performance of the BZB very poor in the transversal direction. Compared to the OECS model, as seen in Figure 3, the MECS model can offer all the vibration modes of the BZB including the global coupled modes and the local modes of cable stays.
Comparison of natural frequencies for the first ten global modes.
Mode order
OESC (Hz)
MESC (Hz)
Structural component
Description (MESC)
1
0.2725
0.2756
Deck
Antisymmetric vertical bending
2
0.2939
0.2925
Main Tower
Lateral bending
3
0.2961
0.2947
Main Tower
Longitudinal bending
4
0.4429
0.4470
Deck
Torsional bending
5
0.5424
0.5600
Deck
Symmetric vertical bending
6 (41)
0.6592
0.6686
Deck
Antisymmetric vertical bending
7 (42)
0.6737
0.6714
Minor tower 1, 2
Lateral + longitudinal bending
8 (56)
0.6748
0.6728
Minor tower 1, 2
Longitudinal bending
9 (57)
0.6996
0.6972
Minor tower 1, 2
Lateral bending
10 (69)
0.6996
0.6972
Minor tower 1, 2
Longitudinal + lateral bending
Note. The numbers in brackets are the mode order of the MESC model.
The lateral mode shapes of the MESC model.
The 2nd-order mode shape
The 42th-order mode shape
The 57th-order mode shape
The 69th-order mode shape
2.3. Possibility of the Transverse Resonance between the Local Model and Global Model
Table 2 shows the relationship between the natural frequencies of the global modes and part of the cables in the BZB. It is seen that the transverse coupled vibration between the tower and cable due to the 1 : 2 or 1 : 1 internal resonance mechanism are possibly to occur. For example, the natural frequencies of cables N23, N24, N25, N26, and N13, N14, N15 are in the one times or two times vicinity of the natural frequency of the 2nd lateral bending mode, respectively. On the other hand, the minor tower and cables also occur the transverse coupled vibration, such as the 7th, 8th, 9th, and 10th lateral bending modes for the cables A9 and J8.
Cables which may take place of the transverse coupled vibration between the global mode and cable.
Mode number (OECS)
Natural frequency (Hz)
Cable number
Natural frequency (Hz)
Ratio
2nd
0.2939
N23
0.313
0.939 (1 : 1)
N24
0.310
0.948 (1 : 1)
N25
0.289
1.017 (1 : 1)
N26
0.271
1.085 (1 : 1)
N13
0.648
0.454 (1 : 2)
N14
0.554
0.531 (1 : 2)
N15
0.536
0.548 (1 : 2)
7th, 8th
0.6737
A9
1.392
0.484 (1 : 2)
J8
1.373
0.491 (1 : 2)
9th, 10th
0.6996
A9
1.392
0.503 (1 : 2)
J8
1.373
0.510 (1 : 2)
3. Simplified Model3.1. Equations of Motion
In this section, the coupled vibrations between the cable and tower of the BZB under transverse excitation have been investigated by a simplified 2-DOFs model, in which the tower is reduced to a beam. The simplified model consists of a beam and a cable, each of them fixed at one end and attached to the other end as considered in Figure 4. The beam is considered as an Euler-Bernoulli beam and the cable is simulated ignoring the bending, torsional, and shear rigidities behavior. In addition, the beam and cable are assumed to be homogeneous and oscillate transversely only in outplane.
The simplified model.
Static model
Dynamic model
Considering the previous assumptions, the equations of motion for the simplified model are obtained by using the extended Hamilton principle:(1a)mbw¨b+ξbw˙b+EbIbwb′′′′+Nwb′′=Pbcos(Ωbt),(1b)mcw¨c+ξcw˙c-[Hwc′+EcAcwc′e(t)]′=Pccos(Ωct),(2)N=(H+EcAce(t))sin(θ),e(t)=1lc∫0lcwc′2dxc,
where the symbols mb and mc are the mass per unit length of the beam and cable, respectively; Ib and Ic are the lengths of the beam and cable, respectively; EbIb and EcAc are the bending and axial stiffness of the beam and cable, respectively; Eb, Ec are the Young’s modulus of the material; Ib and Ac are the moment of area and cross-sectional area, respectively; N is the axial compressive load; H is the initial tension of the cable; ξb and ξc are the damping coefficient of the beam and cable, respectively; Ωb, Ωc, Pb, and Pc are the frequency and amplitude of the transverse loads, respectively; wb is the beam transverse displacement at location xb; and wc is the transverse displacement of the cable at location xc. The overdot indicates the differentiation with respect to the time t; the prime indicates the differentiation with respect to the coordinate x.
The associated geometric and relevant mechanical boundary conditions of the beam and cable are, respectively, written as
(3)wb(0,t)=wc(0,t),wb(lb,t)=0,wc(lc,t)=0,wb′(lb,t)=0,wb′′(0,t)=0,EbIbwb′′′(0,t)=λwb(0,t),
where the λ=EcAc/lc denotes the out-of-plane supported stiffness of the beam caused by the cable.
For convenience, a set of new variables and parameters are defined as
(4)τ=ωt,x-c=xclc,x-b=xblb,w-c=wclc,w-b=wblb,ξ-b=ξbmbω,ξ-c=ξcmcω,P-b=Pbmblbω2,P-c=Pcmclcω2,ρ=mcmb,μ=EcAcH,χ=EbIblb2EcAc,βb4=ω2mblb4EbIb,βc2=ω2mclc2H,
where ω is the natural frequency of the simplified model out-plane. In the nondimensional form, (1a)–(3) become(5a)w¨b+ξbw˙b+1βb4wb′′′′+ρwb′′βc2sin(θ)+ρμe(t)βc2sin(θ)wb′′=Pbcos(Ωbt),(5b)w¨c+ξcw˙c-1βc2wc′′-μe(t)βc2wc′′=Pccos(Ωct),(6)e(t)=∫0112wc′2dx,wb(0,t)sin(θ)=wc(0,t),wb(1,t)=0,wc(1,t)=0,wb′(1,t)=0,wb′′(0,t)=0,χwb′′′(0,t)=(1μ+e(t))wc′(0).
In this paper, only the first-order modes of the beam and cable are considered. Therefore, based on the research [15, 16], the transverse displacements wb(x,t) and wc(x,t) are approximated by the transverse modes of the first order as follows:(7a)wb(x,t)=ψ1(x)p1(t),(7b)wc(x,t)=ψ1(x)p1(t)sin(θ)+ψ2(x)p2(t),where φ1(x) and φ2(x) are the mode shapes of the beam and cable, respectively. Both of them have the following form:(8a)ψ1(x)=A1sin(βbx)+A2cos(βbx)+A3sinh(βbx)+A4cosh(βbx),(8b)ψ2(x)=sin(πx).
Using the Galerkin approach, substituting (7a) and (7b) into (5a) and (5b), the nonlinear governing equations of motion with 2-DOFs for the simplified model are obtained as follows:(9a)p¨1(t)+ξ1p˙1(t)+a1p1(t)+a122p1(t)p2(t)2+a112p1(t)2p2(t)+a111p1(t)3=fbcos(Ωbt),(9b)p¨2(t)+ξ2p˙2(t)+d1p¨1(t)+d2p˙1(t)+b2p2(t)+b112p1(t)2p2(t)+b122p1(t)p2(t)2+b222p2(t)3+b1p1(t)+b111p1(t)3=fccos(Ωct),(10)e(t)=12p2(t)2I4+p1(t)p2(t)sin(θ)I5+12p1(t)2sin(θ)2I6,
where fb and fc are the amplitude of the harmonic functions and ai, bi, di, aij, bij, aijk, and bijk are the Galerkin coefficients of the simplified model, respectively. All the coefficients are defined in the Appendix.
3.2. Perturbation Analysis
The multiple scales perturbation method [17] is applied to (9a) and (9b) to obtain an approximation solution of the model. To make the nonlinear terms weak, one can substitute aijk, bijk, and cijk with εaijk, εbijk, and εcijk. Then, (9a) and (9b) can be rewritten as(11a)p¨1(t)+εξ1p˙1(t)+ω12p1(t)+εa122p1(t)p2(t)2+εa112p1(t)2p2(t)+εa111p1(t)3-εfbcos(Ωbt)=0,(11b)p¨2(t)+εξ2p˙2(t)+εd1p¨1(t)+εd2p˙1(t)+ω22p2(t)+εb112p1(t)2p2(t)+εb122p1(t)p2(t)2+εb222p2(t)3+εb1p1(t)+εb111p1(t)3=0.
Considering the fact there is 1 : 1 internal resonance between the modes of the beam and cable, primary resonance for the beam and autoparametric resonance for the cable, simultaneously, therefore, the resonant relations are represented as
(12)Ω=ω1+εσ1,ω2=ω1.
Assuming the first-order approximation solution of (11a) and (11b) in the form(13a)p1(t)=p10(T0,T1)+εp11(T0,T1),(13b)p2(t)=p20(T0,T1)+εp21(T0,T1),the time derivatives become(14a)ddt=∂∂T0∂T0∂t+∂∂T1∂T1∂t+⋯=D0+εD1+⋯,(14b)d2dt2=(D0+εD1+⋯)2=D02+2εD0D1+⋯,where Tn=εnt(n=0,1) are the fast and slow time scales, respectively.
Substituting (13a), (13b), (14a), and (14b) into (11a) and (11b) and equating the powers of ε(15a)D1,1(p10)+ω12p10=0,(15b)D1,1(p20)+ω22p20=0,(16a)D1,1(p11)+ω12p11+a112p102p20+2D1,2(p10)+ξ1D1(p10)+a122p10p202+a111p103-cos(Ωt)fb=0,(16b)D1,1(p21)+ω22p21+2D1,2(p20)+ξ2D1(p20)+b112p102p20+b122p10p202+b111p103+d1D1,1(p10)+d2D1(p10)+b1p10+b222p203=0.
The general solution of (15a) and (15b) can be expressed in the form(17a)p10=A1(T1)exp(ω1T0)+A-1(T1)exp(-ω1T0),(17b)p20=A2(T1)exp(ω2T0)+A-2(T1)exp(-ω2T0),where A1 and A2 are complex functions and A1- and A2- denote complex conjugate terms, respectively.
Substituting (17a) and (17b) and (12), into (16a) and (16b) and setting the coefficients of the secular terms to zero yield the solvability conditions as(18a)2ID1(A1)ω1+a112A2¯A12+a122A1¯A22+3a111A1¯A12-12fbexp(IT0εσ1)+2a112A1A1¯A2+Iξ1ω1A1+2a122A1A2A2¯=0,(18b)ID2ω1A1+2b122A1A2A2¯+2ID2(A2)ω1+b112A12A2¯+b1A1+Iξ2ω1A2+2b112A1A1¯A2+b122A1A22+3b222A22A2+3b111A12A1-d1ω12A1=0.Let
(19)Ak=12(xk+iyk)(k=1,2).
Substituting (19) into (18a) and (18b) and then separating the real and imaginary parts, the modulation equations obtained in the Cartesian form are as follows:(20a)ddtx1(t)=-3y13+3y1x128ω1a111+(x1x2-3y1y2)x12+y128ω1a112+(2x1x2y2-y1x22-3y1y22)8ω1a122+fb+2y1σ1ω1-x1ξ1ω12ω1,(20b)ddty1(t)=-3x13+3x1y128ω1a111+(y1x2+3x1y2)x12+y128ω1a112+(2y1x2y2+x1x22+3x1y22)8ω1a122+-2x1σ1ω1-y1ξ1ω12ω1,(20c)ddtx1(t)=-(3y12+3x12)x12+y128ω1b111-(3x12y2+3y12y2)8ω1b112-(x22+3y22)x12+y128ω1b122-3x22y2+3y238ω1b222-(b1-d1ω12)x12+y12+x2ξ2ω1-2y2σ1ω12ω1-y2x12+y12(x1ddty1(t)+y1ddtx1(t)),(20d)ddty1(t)=(x12x2+y12x2)8ω1b112+x2y2x12+y124ω1b122+3x2y22+3x238ω1b222-d2ω1x12+y12+y2ξ2ω1+2x2σ1ω12ω1+x2x12+y12(x1ddty1(t)+y1ddtx1(t)).
It is seen that the (20a), (20b), (20c), and (20d) describe a nonlinear dynamic system, indicating that the transverse vibrations of the cable and tower are nonlinear, even though the beam model considered is linear and neglects the cable’s sag effect. The nonlinear terms are only caused by the coupled behaviors between the modes of the beam and cable. It can be demonstrated that the simplified model shows the chaotic motion and period motion with the perturbation force fb changing.
3.3. Numerical Results and Discussion
In fact, there are many internal resonance forms between the modes of a beam and a cable when they act as an overall structure, such as two-to-one, one-to-one, and one-to-two. This study only focuses on the nonlinear vibrations of the simplified model in the one-to-one internal resonance case, taking into consideration the primary resonance to the beam. The parameters of the simplified model are selected as follows: mc=67.37 kg/m, Ec=2.10e11 N/m2, Ac=8.58e-3m2, lc=152.70 m, H=609.00 kN, mb=1.34e5 kg/m, lb=76.81 m, θ=30.20°, and EbIb=8.29e11 N·m2. Therefore, the nondimension parameters are calculated as follows: ρ=0.001, χ=0.068, μ=2959.3, ξ1=0.02, ξ2=0.001, and fc=0. The amplitude of the beam (fb) is chosen as a controlling parameter. The time histories, phase portraits, and Poincare maps are plotted to analyze the nonlinear dynamical motion of the simplified model. From the numerical simulations, the coupling motion between the beam and cable can be clearly found.
Figure 5 shows the nonlinear behavior of the simplified model at the force amplitude, fb=0.005, involving time histories, phase portraits, and Poincare maps. As seen in Figures 5(a) and 5(b), the amplitude history of the beam varies with a particular period, while the cable is disorders and with many subharmonics. However, there are still energy transformation occurred. In order to identify the characteristic of these motions in the time histories, the phase portraits have been calculated. As seen in Figures 5(c) and 5(d), the phase portraits of the beam and cable exhibit periodicity and nonperiodicity, respectively, which are consistent with the description in Figures 5(a) and 5(b). Since the phase portraits cannot provide enough information to determine the onset for chaotic motion, they only are used to distinguish whether the model is periodic or non-periodic. Therefore, the Poincare maps have been further calculated. As seen in Figure 5(e), the Poincare map of the beam exhibits many irregular points, confirming that the motion is chaos. On the other hand, Figure 5(f) exhibits that the return points in the Poincare map form a closed curve. Generally, this shape indicates that this motion is periodic [18]. However, since the motion of the simplified model is coupled, the beam’s motion takes on chaos characteristic, and the motion of the cable can be confirmed as chaos motion.
Nonlinear response of the simplified model when the fb=0.005 (98.3 KN): (a) Time history; (b) Local magnification of (a); (c) Phase portrait of the beam; (d) Phase portrait of the cable; (e) Poincare map of the beam; (f) Poincare map of the cable; (continued line beam, dashed line cable; p10=0.001, p20=0.032, Dp10=0.001, Dp20=-0.01).
Figure 6 shows the nonlinear behavior of the simplified model at the force amplitude, fb=0.01. As shown in Figures 6(a) and 6(b), the varying amplitudes of time history of beam and cable exhibit particular period while the cable has a few subharmonics. The phenomenon of energy transformation between them can be clearly found. As seen in Figures 6(c) and 6(d), the phase portraits of the beam and cable both exhibit periodicity and are further confirmed by Figures 6(e) and 6(f). These phenomena indicate that the motions of the beam and cable are both quasi-periodic.
Nonlinear response of the simplified model when the fb=0.01 (197.46 kN): (a) Time history; (b) Local magnification of (a); (c) Phase portrait of the beam; (d) Phase portrait of the cable; (e) Poincare map of the beam; (f) Poincare map of the cable; (continued line beam, dashed line cable; p10=0.001, p20=0.032, Dp10=0.001, Dp20=-0.01).
When the forcing amplitude is increased to fb=0.1, the nonlinear behavior of the simplified model changes to chaotic motions, as shown in Figure 7. The Poincare maps given in Figures 7(e) and 7(f) demonstrate clearly that chaotic motions exist in the simplified model again. It can be observed that the time-histories and phase portraits represented by Figures 7(a), 7(b), 7(c), and 7(d) are very similar to that of Figure 6; only the complexity of these graphics increased. On the other hand, compared to Figures 5(a), 6(a), and 7(a), the soften phenomenon in the simplified model tends to be more apparent with the force amplitude (fb) increasing.
Nonlinear response of the simplified model when the fb=0.1 (1974.58 kN): (a) Time history; (b) Local magnification of (a); (c) Phase portrait of the beam; (d) Phase portrait of the cable; (e) Poincare map of the beam; (f) Poincare map of the cable; (continued line beam,dashed line cable; p10=0.001, p20=0.032, Dp10=0.001, Dp20=-0.01).
4. Conclusions
The frequency values of the Bin-Zhou cable-stayed bridge (BZB) influenced by the cable vibration have been investigated by two FE models in this study. One is the OECS model in which one single element per cable stay is used and the other is MECS model, where multi-elements per cable stay are used. The nonlinear behaviors of the cable vibration of BZB have also been examined by a simplified model, where the tower is simplified as a beam. The motions of the simplified model are utilized by Galerkin’s method to truncate a two DOFs nonlinear coupled model. Based on above numerical experiments, some conclusions are summarized as follows.
The local vibration of the cables only has a small impact on the frequency values of the BZB but increase for the deck and decrease for the towers.
The results from the modal analysis also show that the modes of the towers are main component of the BZB in the first 10 modes, and most of them are the lateral vibrations. Compared to the OESC model, the MECS model not only offers global modes of the bridge, but also exhibits the local vibration of the cables.
The results from the simplified model analysis show that the coupled system exhibits quasi-periodic and chaotic motion with the forcing amplitude changing, even though the beam model is linear and the cable model neglects the cable’s sag effect. The simplified model also exhibits soften behavior with the forcing amplitude increasing.
The energy transform can be found due to the 1 : 1 internal resonant between the modes of the beam and cable.
Appendix
One has
(A.1)I1=∫01D(ψ2)(x)2dx,I2=∫01D(ψ2)(x)D(ψ1)(x)dx,I3=∫01D(ψ1)(x)2dx,I4=∫01ψ1(x)2dx,I5=∫01ψ1(x)D(4)(ψ1)(x)dx,I6=∫01ψ1(x)D(2)(ψ1)(x)dx,I7=∫01ψ1(x)Pwb(x)dx,I8=∫01ψ2(x)2dx,I9=∫01ψ1(x)ψ2(x)dx,I10=∫01ψ2(x)D(2)(ψ2)(x)dx,I11=∫01ψ2(x)D(2)(ψ1)(x)dx,I12=∫01ψ2(x)Pwc(x)dx,a111=12ρμsin(θ)I6I3βc2I4,a112=ρμI6I2βc2I4,a122=12ρμI6I1βc2sin(θ)I4,a1=I5βb4I4+ρI6βc2sin(θ)I4,fb=I7I4,b222=-12μI10I1I8βc2,b122=-μsin(θ)I10I2I8βc2-12μsin(θ)I11I1I8βc2,b111=-12μsin(θ)3I11I3I8βc2,b112=-12μsin(θ)2I10I3I8βc2-μsin(θ)2I11I2I8βc2,b2=-I10I8βc2,b1=-sin(θ)I11I8βc2,d1=sin(θ)I9I8,d2=ξ2sin(θ)I9I8,fc=I12I8.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The present study was sponsored by National Natural Science Foundation of China (Grants nos. 51308322, 51078119), China Postdoctoral Science Foundation (Grant no. 2013M541387, 2013M530049), Fund of State Key Laboratory of Bridge Engineering Structural Dynamics (Grant no. 201301), and Fund of National Engineering and Research Center for Highways in Mountain Area (Grant no. gsgzj-2013-03).
NiY. Q.WangJ. Y.LoL. C.Influence of stabilizing cables on seismic response of a multispan cable-stayed bridge20052021421532-s2.0-1324426311710.1111/j.1467-8667.2005.00383.xHeG. J.ZouZ. Q.NiY. Q.KoJ. M.Seismic response analysis of multi-span cable-stayed bridge2009400–4027377422-s2.0-5884913961910.4028/www.scientific.net/KEM.400-402.737LiuJ. L.2009Harbin, ChinaHarbin Institute of Technology(Chinese)OkamotoY.NakamuraS.Static and seismic studies on steel/concrete hybrid towers for multi-span cable-stayed bridges20116722032102-s2.0-7795728241810.1016/j.jcsr.2010.08.008FujinoY.WarnitchaiP.PachecoB. M.An experimental and analytical study of autoparametric resonance in a 3DOF model of cable-stayed-beam1993421111382-s2.0-002758576110.1007/BF00045250XiaY.FujinoY.Auto-parametric vibration of a cable-stayed-beam structure under random excitation200613232792862-s2.0-3294447024110.1061/(ASCE)0733-9399(2006)132:3(279)XiaY.WuQ.-X.XuY.-L.FujinoY.ZhouX.-Q.Verification of a cable element for cable parametric vibration of one-cable-beam system subject to harmonic excitation and random excitation20111435895952-s2.0-7995923803810.1260/1369-4332.14.3.589CaetanoE.CunhaA.TaylorC. A.Investigation of dynamic cable-deck interaction in a physical model of a cable-stayed bridge. Part I: modal analysis20002944814982-s2.0-0034175537CaetanoE. .CunhaA.TaylorC. A.Investigation of dynamic cable-deck interaction in a physical model of a cable-stayed bridge. Part II: seismic response2000294499521GeorgakisC. T.TaylorC. A.Nonlinear dynamics of cable stays. Part 1: sinusoidal cable support excitation20052813–55375642-s2.0-1344425941910.1016/j.jsv.2004.01.022GeorgakisC. T.TaylorC. A.Nonlinear dynamics of cable stays. Part 2: stochastic cable support excitation20052813–55655912-s2.0-1344429292710.1016/j.jsv.2004.01.023CaetanoE.CunhaA.GattulliV.LepidiM.Cable-deck dynamic interactions at the International Guadiana Bridge: on-site measurements and finite element modelling20081532372642-s2.0-4244909599010.1002/stc.241NiY. Q.ZhouH. F.ChanK. C.KoJ. M.Modal flexibility analysis of cable-stayed Ting Kau bridge for damage identification20082332232362-s2.0-3934909562210.1111/j.1467-8667.2008.00521.xWangP. H.LiuM. Y.HuangY. T.LinL. C.Influence of lateral motion of cable stays on cable-stayed bridges20103467197382-s2.0-7795026009010.12989/sem.2010.34.6.719ChengG.ZuJ. W.Dynamic analysis of an optical fiber coupler in telecommunications2003268115312-s2.0-024258396410.1016/S0022-460X(02)01575-4WuQ.TakahashiK.OkabayashiT.NakamuraS.Response characteristics of local vibrations in stay cables on an existing cable-stayed bridge200326134034202-s2.0-003746875210.1016/S0022-460X(02)01088-XNayfehA. H.MookD. T.1979New York, NY, USAJohn Wiley & SonsMR549322Chang-JianC.-W.ChenC.-K.Nonlinear analysis of a rub-impact rotor supported by turbulent couple stress fluid film journal bearings under quadratic damping20095632973142-s2.0-6394908485410.1007/s11071-008-9400-0ZBL1170.70356