International Scholarly Research Notices

International Scholarly Research Notices / 2011 / Article

Research Article | Open Access

Volume 2011 |Article ID 154806 |

Shanshan Wang, Yukun An, "Multiple Periodic Solutions to a Suspension Bridge Wave Equation with Damping", International Scholarly Research Notices, vol. 2011, Article ID 154806, 9 pages, 2011.

Multiple Periodic Solutions to a Suspension Bridge Wave Equation with Damping

Academic Editor: T. Tran
Received03 Sep 2011
Accepted11 Oct 2011
Published27 Dec 2011


This paper is concerned with the existence of multiple periodic solutions for a suspension bridge wave equation with damping. By using Leray-Schauder degree theory, the authors prove that the damped wave equation has multiple periodic solutions.

1. Introduction

In [1], also see [2–6], the author considered a horizontal cross-section of the center span of a suspension bridge and proposed a partial differential equation model for the torsional motion of the cross-section and treated the center span of the bridge as a beam of length 𝐿 and width 2𝑙 suspended by cables. Consider the horizontal cross-section of mass 𝑚 located at position 𝑥 along the length of the span. She treated this cross-section as a rod of length 2𝑙 and mass 𝑚 suspended by cables. Let 𝑦(𝑥,𝑡) denote the downward distance of the center of gravity of the rod from the unloaded state and let 𝜃(𝑥,𝑡) denote the angle of the rod from horizontal at time 𝑡. Assume that the cables do not resist compression, but resist elongation according to Hooke’s law with spring constant 𝐾. Then the torsional and vertical motion of the span satisfy 𝜃𝑡𝑡−𝜀1𝜃𝑥𝑥=3𝐾𝑚𝑙cos(𝑦−𝑙sin𝜃)+−𝛿𝜃𝑡+ℎ1𝑦(𝑥,𝑡),𝑡𝑡+𝜀2𝑦𝑥𝑥𝑥𝑥𝐾=−𝑚(𝑦−𝑙sin𝜃)++(𝑦+𝑙sin𝜃)+−𝛿𝑦𝑡+𝑔+ℎ2(𝑥𝑡),𝜃(0,𝑡)=𝜃(𝐿,𝑡)=𝑦(0,𝑡)=𝑦(𝐿,𝑡)=𝑦𝑥𝑥(0,𝑡)=𝑦𝑥𝑥(𝐿,𝑡)=0,(1.1) where 𝑢+=max{𝑢,0}, 𝜀1, 𝜀2 are physical constants related to the flexibility of the beam, 𝛿 is the damping constant, ℎ1 and ℎ2 are external forcing terms, and 𝑔 is the acceleration due to gravity. The spatial derivatives describe the restoring force that the beam exerts, and the time derivatives 𝜃𝑡 and 𝑦𝑡 represent the force due to friction. The boundary conditions reflect the fact that the ends of the span are hinged.

Throughout the paper [1] the author assumes that the cables never lose tension; that is, it is assumed that (𝑦±sin𝜃)≥0. In this case, we see that (1.1) becomes uncoupled, and the torsional and vertical motions satisfy, respectively, 𝜃𝑡𝑡−𝜀1𝜃𝑥𝑥=−6𝐾𝑚cos𝜃sin𝜃−𝛿𝜃𝑡+ℎ1𝑦(𝑥,𝑡),𝜃(0,𝑡)=𝜃(𝐿,𝑡)=0,(1.2)𝑡𝑡+𝜀2𝑦𝑥𝑥𝑥𝑥=−2𝐾𝑚𝑦−𝛿𝑦𝑡+𝑔+ℎ2(𝑥,𝑡),𝑦(0,𝑡)=𝑦(𝐿,𝑡)=𝑦𝑥𝑥(0,𝑡)=𝑦𝑥𝑥(𝐿,𝑡)=0.(1.3)

In paper [1], removing the damping term; that is, let 𝛿=0, changing variables, and imposing boundary and periodicity conditions, the author rewrites (1.2) as 𝑢𝑡𝑡−𝑢𝑥𝑥𝑢+𝑏sin𝑢=ğœ€â„Ž(𝑥,𝑡),(0,𝑡)=𝑢(𝜋,𝑡)=0,𝑢(𝑥,0)=𝑢(𝑥,𝜋),𝑢𝑡(𝑥,0)=𝑢𝑡(𝑥,𝜋),𝑢(𝑥,𝑡)=𝑢(𝜋−𝑥,𝑡),𝑢(𝑥,𝑡)=𝑢(𝑥,𝜋−𝑡).(1.4) And it proves that (1.4) has at least two solutions in the subspace 𝐻 of 𝐿2. Where 𝐻 is defined as 𝐻=𝑢∈𝐿2(Ω)𝑢(𝑥,𝑡)=𝑢(𝜋−𝑥,𝑡),𝑢(𝑥,𝑡)=𝑢(𝑥,𝜋−𝑡),𝑢𝑖𝑠𝜋−periodicin𝑡.(1.5)

Notice that (1.4) is particular in no damping and the selection of 𝐻. Hence, in [1] the author left a problem which is relevant to this case.

Problem 1. “Under appropriate hypotheses on the forcing term, does a similar result hold for the damped equation?”
Motivated by this problem, in this paper, we suppose that the damping is present, that is, 𝛿≠0, and study the following problem: 𝑢𝑡𝑡−𝑢𝑥𝑥+𝛿𝑢𝑡𝑢+𝑏sin𝑢=ğœ€â„Ž(𝑥,𝑡),(0,𝑡)=𝑢(𝜋,𝑡)=0,𝑢(𝑥,0)=𝑢(𝑥,𝜋),𝑢𝑡(𝑥,𝑡)=𝑢𝑡(𝑥,𝜋),𝑢(𝑥,𝑡)=𝑢(𝜋−𝑥,𝑡).(1.6)

2. Preliminaries

Let 𝑁={0,1,…} and 𝑍 be the set of integers, Λ=𝑁×𝑁. Let Ω=(0,𝜋)×(0,𝜋) and 𝐿2(Ω) be usual space of square integrable functions with usual inner product (⋅,⋅) and corresponding norm ‖⋅‖. For the Sobolev space 𝐻1(Ω), we denote the standard inner product by (𝑢,𝑣)1=(𝑢,𝑣)+(𝑢𝑥,𝑣𝑥)+(𝑢𝑡,𝑣𝑡) and norm by ‖𝑢‖1.

Define the operator 𝐿𝛿𝑢=𝑢𝑡𝑡−𝑢𝑥𝑥+𝛿𝑢𝑡∶𝐻→𝐻 by 𝐷𝐿𝛿=𝑢∈𝐻∣𝑢(𝑥,𝑡)=Λ𝑢𝑚𝑛Φ𝑚𝑛,Λ(2𝑛+1)2−4𝑚22+4𝑚2𝛿2||𝑢𝑚𝑛||2,𝐿<âˆžğ›¿î“ğ‘¢=Λ(2𝑛+1)2−4𝑚2𝑢+2𝑚𝛿𝑚𝑛Φ𝑚𝑛,forall𝑢∈𝐻.(2.1)

We know that the eigenvalues and corresponding eigenfunctions of 𝐿𝛿 are 𝜆𝑚𝑛=(2𝑛+1)2−4𝑚2Φ+2𝑚𝛿𝑖,(𝑚,𝑛)∈Λ,𝑚𝑛(𝑥,𝑡)=𝑒2𝑚𝑡𝑖sin(2𝑛+1)𝑥,(𝑚,𝑛)∈Λ.(2.2)

In order to seek the solutions of (1.6), we first investigate the properties of operator 𝐿𝛿. We have the following Lemma.

Lemma 2.1. 𝐿𝛿−1 exists, 𝐿𝛿−1∶𝐻→𝐻 is compact, and ‖𝐿𝛿−1‖=1.

Proof. Because we are restricted to the subspace 𝐻 of 𝐿2, and 𝜆𝑚𝑛≠ 0, we easily know 𝐿𝛿−1 exists.
We prove 𝐿𝛿−1∶𝐻→𝐻 is compact below. We find that 𝐿𝛿−1𝑢=Λ1(2𝑛+1)2−4𝑚2𝑢+2𝑚𝛿𝑖𝑚𝑛Φ𝑚𝑛,(2.3) for all ∑𝑢=Λ𝑢𝑚𝑛Φ𝑚𝑛∈𝐻. For any (𝑚,𝑛)∈Λ, we have ||(2𝑛+1)2−4𝑚2||+2𝑚𝛿𝑖2≥1,(2.4) then ‖‖𝐿𝛿−1𝑢‖‖2=Λ||||1(2𝑛+1)2−4𝑚2𝑢+2𝑚𝛿𝑖𝑚𝑛||||2≤Λ||𝑢𝑚𝑛||2=‖𝑢‖2.(2.5) Hence, ‖‖𝐿𝛿−1𝑢‖‖≤‖𝑢‖.(2.6) On the other hand, ‖‖𝐿𝛿−1𝑢‖‖21=Λ|||||1+(2𝑛+1)2+4𝑚2(2𝑛+1)2−4𝑚22+4𝑚2𝛿2𝑢𝑚𝑛|||||2,(2.7) while 1+(2𝑛+1)2+4𝑚2(2𝑛+1)2−4𝑚22+4𝑚2𝛿2=1+(2𝑛+1)2+4𝑚2(2𝑛+2𝑚+1)2(2𝑛−2𝑚+1)2+4𝑚2𝛿2≤1+(2𝑛+1)2+4|𝑚|2(2𝑛+2|𝑚|+1)2+4|𝑚|2𝛿2≤(2𝑛+2|𝑚|+1)2+1(2𝑛+2|𝑚|+1)2≤2.(2.8) Hence, ‖‖𝐿𝛿−1𝑢‖‖21≤2Λ|𝑢|2=2‖𝑢‖.(2.9) By (2.6) and (2.9), we can find that the operator 𝐿𝛿−1∶𝐻→𝐻 is compact since the embedding 𝐻1→𝐿2 is compact.
Finally, we prove ‖𝐿𝛿−1‖= 1. By (2.2) and ||𝜆𝑚𝑛||2=||(2𝑛+1)2−4𝑚2||+2𝑚𝛿𝑖2≥1.(2.10) Set 𝑢=Φ00, such that, ‖𝐿𝛿−1𝑢‖/‖𝑢‖=1. Therefore, ‖‖𝐿𝛿−1‖‖=1.(2.11) Hence, we complete the proof of this lemma.

Definition 2.2. One says that 𝑢∈𝐻 is a solution to (1.6) if 𝑢=𝐿𝛿−1(ğœ€â„Žâˆ’ğ‘sin𝑢).(2.12)
To establish the existence of multiple periodic solutions to (1.6), we use Leray-Schauder degree theory to prove the existence of multiple zeros of a related operator 𝑇1. To compute the degree of 𝑇1, we continuously deform it to a linear operator 𝑇0, the Gâteaux derivative of 𝑇1, and compute its degree via a direct calculation.
It is not difficult to show that the homotopy property of Leray-Schauder degree ensures that the degree of an operator 𝑇1 is preserved as 𝑇1 is continuously deformed to its Fréchet derivative under appropriate hypotheses. However, the nonlinear term in (1.6), 𝑓(𝑢)=sin𝑢, is not Fréchet differentiable in 𝐿2 at 𝑢=0.
There is a theorem in paper [1], in which, the author shows that, under certain conditions on the nonlinear term 𝑓 and the differential operator 𝐿, Leray-Schauder degree is indeed preserved under homotopy from the operator 𝑇1 to its Gâteaux derivative 𝑇0. This result can be used to establish multiplicity of solutions to equations of the form (1.6). The result follows.

Lemma 2.3. Let 𝐼1, 𝐼2 be open, bounded intervals in 𝑅, and define 𝑄∶=𝐼1×𝐼2. Let 𝐵 be a subspace of 𝐿𝑝(𝑄), 𝑝≥1, and define ‖u‖∶=‖u‖𝐿𝑝. Consider the problem 𝐿𝑢+𝑓(𝑢)=ğœ€â„Ž(𝑥,𝑡),(2.13) where 𝐿, 𝑓, and ℎ satisfy the following:(H1)𝐿−1 is compact;(H2)‖𝐿−1‖≤ 1;(H3)𝑓(0)=0;(H4)𝑓 is Lipschitz with Lipschitz constant 𝑀;(H5)â„Žâˆˆğµ and ℎ≤1;(H6) the Gâteaux derivative 𝑑𝑓(0,𝑢) exists and satisfies 𝑑𝑓(0,𝑢)=𝜌𝑢, where 𝜌>0 and −𝜌 is not an eigenvalue of 𝐿.Define 𝑇0: 𝐵→𝐵 by 𝑇0(𝑢)=𝑢+𝜌𝐿−1(𝑢),(2.14) and 𝑇1∶𝐵→𝐵 by 𝑇1(𝑢)=𝑢−𝐿−1(ğœ€â„Žâˆ’ğ‘“(𝑢)).(2.15) Then for 𝜀 sufficiently small, there exists 𝛾>0 such that 𝑇deg1,𝐵𝛾(0),0=deg𝑇,𝐵𝛾(0),0.(2.16)

3. Result and Proof

The main result of this paper is as follows.

Theorem 3.1. Let â„Žâˆˆğ» with ‖ℎ‖≤1, and let √𝑏∈(−25+4𝛿2√,−9+4𝛿2), 0<𝛿2<14. Then there exists 𝜀0>0 such that if |𝜀|<𝜀0, (1.6) has at least two solutions in 𝐻.

Proof. Let 𝐿=𝐿𝛿 and 𝑓(𝑢)=𝑏sin𝑢, it is easy to know that 𝐿 and 𝑓 satisfy the conditions (H1–H5) in Lemma 2.3.
Reply Lemma 2.3, we define 𝑇0: 𝐻→𝐻 by 𝑇0(𝑢)=𝑢+𝑏𝐿𝛿−1(𝑢),(3.1) and 𝑇1: 𝐻→𝐻 by 𝑇1(𝑢)=𝑢−𝐿𝛿−1(ğœ€â„Žâˆ’ğ‘sin(𝑢)).(3.2) And note that zeros of 𝑇1 correspond to solutions of (1.6). To prove the theorem, we will show the following:(C1) there exists 𝑅0>0 such that for 𝑅>𝑅0, deg(𝑇1,𝐵𝑅(0),0)=1;(C2) there exists 𝛾∈(0,𝑅0) such that deg(𝑇1,𝐵𝛾(0),0)=−1.Then, since deg(𝑇1,𝐵𝛾(0),0)≠0, there exists a zero of 𝑇1 (i.e., a solution of (1.6)) in 𝐵𝛾(0). Moreover, by the additivity property of degree, deg(𝑇1,𝐵𝑅(0)⧵𝐵𝛾(0),0)≠0 and hence (1.6) has a second solution in the annulus 𝐵𝑅(0)⧵𝐵𝛾(0).
To establish (C1), define 𝑇𝛽𝑢=𝑢−𝛽𝐿𝛿−1(ğœ€â„Žâˆ’ğ‘sin(𝑢)),(3.3) or 𝛽∈[0,1], and note that this definition of 𝑇1 is consistent with our previous definition. Note also that 𝑇0 is simply the identity map; hence, for any 𝑅>0 we have deg(𝑇0,𝐵𝑅(0),0)=1. The homotopy property of degree ensures that deg(𝑇𝛽,𝐵𝑅(0),0) is constant provided that 0∈𝑇𝛽(𝜕𝐵𝑅(0)) for all 𝛽∈[0,1].
Fix 𝛽∈[0,1] and suppose 𝑢∈𝐻 solves 𝑇𝛽𝑢=0. We will show that u is bounded above by some 𝑅0>0 and that this bound is independent of 𝛽.
Since 𝑇𝛽𝑢=0, we have ‖‖𝐿‖𝑢‖=𝛽𝛿−1‖‖𝜀(ğœ€â„Žâˆ’ğ‘sin𝑢)≤𝛽0≤𝜀+𝑏‖sin𝑢‖0+𝑏𝑚(Ω)1/2<𝜀0√+𝑏2𝜋<𝑅0,(3.4) if we choose 𝑅0>𝜀0√+𝑏2𝜋.
Thus, for 𝑅>𝑅0, we have 𝑇deg1,𝐵𝑅𝑇(0),0=deg0,𝐵𝑅(0),0=1,(3.5) and (C1) above holds.
To establish (C2), let 𝜀<𝜀0; we will determine the value of 𝜀0 later. For 𝜇∈[0,1] define 𝑇𝜇𝑢=𝑢+(1−𝜇)𝐿𝛿−1(𝑏𝑢)−𝜇𝐿𝛿−1(ğœ€â„Žâˆ’ğ‘sin𝑢),(3.6) and note again that this definition of 𝑇1 is consistent with our previous definitions. We will again apply the homotopy property of degree (via Lemma 2.3) and a standard degree calculation to show that for some 𝛾>0𝑇deg1,𝐵𝛾𝑇(0),0=deg0,𝐵𝛾(0),0=−1.(3.7)
Observe that for 𝐿=𝐿𝛿 and 𝑓(𝑢)∶=𝑏sin𝑢, hypotheses (H1)–(H5) of Lemma 2.3 are satisfied. To verify hypothesis (H6), we need to show that 𝑑𝑓(0,𝑢)=𝑏𝑢.(3.8) By definition of the Gâteaux derivative, 𝑑𝑑𝑓(0,𝑢)=𝑑𝑡𝑓(0+𝑡𝑢)|𝑡=0=limℎ⟶0𝑓((𝑡+ℎ)𝑢)−𝑓(𝑡𝑢)ℎ|𝑡=0=limℎ⟶0𝑏sin(â„Žğ‘¢)ℎ.(3.9) We will show that the limit above (in 𝐻) is 𝑏𝑢.
Note first that in 𝑅 we have limℎ→0=sin(â„Žğ‘¢)ℎ=limℎ→0sin(â„Žğ‘¢)â„Žğ‘¢ğ‘¢=𝑢,(3.10) and hence |||sin(â„Žğ‘¢)ℎ|||−𝑢2⟶0,(3.11) as ℎ→0. Invoking the convexity of 𝜔2, we have |||sin(â„Žğ‘¢)ℎ|||−𝑢21≤42|||sin(â„Žğ‘¢)ℎ|||2+12|𝑢|2≤4𝑢2.(3.12)
Since 𝑢∈𝐿2, |(sin(â„Žğ‘¢)/ℎ)−𝑢|2 is dominated in 𝐿1; thus by the dominated convergence theorem, ‖‖‖𝑏sin(â„Žğ‘¢)â„Žâ€–â€–â€–âˆ’ğ‘ğ‘¢âŸ¶0,(3.13) as ℎ→0; therefore (3.8) holds. Moreover, by the form of eigenvalue of 𝐿𝛿 and our choice of 𝑏, −𝑏 is not an eigenvalue of 𝐿𝛿; therefore hypothesis (H6) of Lemma 2.3 holds. Thus, by Lemma 2.3, for sufficiently small 𝛾,𝜀>0, we have 𝑇deg0,𝐵𝛾𝑇(0),0=deg1,𝐵𝛾(0),0.(3.14)
Finally, we will show that 𝑇deg0,𝐵𝛾(0),0=deg𝐼+𝑏𝐿𝛿−1,𝐵𝛾(0),0=−1.(3.15)
Consider the finite dimensional subspace 𝑀𝑁=span{Φ𝑚𝑛}𝑁0 of 𝐻 and recall that, by compactness, 𝑏𝐿𝛿−1 can be approximated in operator norm by the operators 𝐵𝑁∶𝑀𝑁→𝑀𝑁 given by 𝐵𝑁(𝑢)=𝑏𝑁𝑁−𝑁0𝑐𝑚𝑛𝜆𝑚𝑛Φ𝑚𝑛.(3.16) By definition of Leray-Schauder degree, for 𝑁 sufficientlylarge, 𝑇deg0,𝐵𝛾(0),0=deg𝐼+𝐵𝑁,𝐵𝛾𝑀(0)𝑁=,0𝑢∈𝐼+𝐵𝑁−1(0)sgn𝐽𝐼+𝐵𝑁(𝑢),(3.17) where 𝐽Φ(𝑢) is the Jacobian determinant of Φ at 𝑢.
Since 𝐼+𝐵𝑁 can be identified with an (2𝑁+1)2×(2𝑁+1)2 diagonal matrix whose entries are 1+𝑏/𝜆𝑚𝑛, we have deg𝐼+𝐵𝑁,𝐵𝛾𝑀(0)𝑁,0=sgn𝑁𝑁𝑚=−𝑁𝑛=0𝑏1+𝜆𝑚𝑛.(3.18) Now we consider the following two cases.(D1)If 𝜆𝑚𝑛 contains imaginary part, suppose a pair of conjugate complex numbers are ğ‘ŽÂ±ğ‘ğ‘–(𝑐≠0), then, 𝑏1+ğ‘ğ‘Ž+𝑐𝑖1+=ğ‘Žâˆ’ğ‘ğ‘–(ğ‘Ž+𝑏)2+𝑐2ğ‘Ž2+𝑐2>0.(3.19)(D2)If 𝜆𝑚𝑛 is real, then 𝑚=0, here 𝜆𝑚𝑛=(2𝑛+1)2.
Since √𝑏∈(−25+4𝛿2√,−9+4𝛿2), and 0<𝛿2<14, the only negative value of 1+𝑏/(2𝑛+1)2 occurs at 𝜆00. Therefore, deg𝐼+𝐵𝑁,𝐵𝛾𝑀(0)𝑁,0=−1.(3.20) The proof of the theorem is complete.


  1. K. S. Moore, “Large torsional oscillations in a suspension bridge: multiple periodic solutions to a nonlinear wave equation,” SIAM Journal on Mathematical Analysis, vol. 33, no. 6, pp. 1411–1429, 2002. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  2. P. J. McKenna, “Large torsional oscillations in suspension bridges revisited: fixing an old approximation,” The American Mathematical Monthly, vol. 106, no. 1, pp. 1–18, 1999. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  3. K. S. Moore, Large-amplitude torsional oscillations in a nonlinearly suspended beam: a theoretical and numerical investigation, dissertation, University of Connecticut, Storrs, Conn, USA, 1999.
  4. P. J. McKenna and K. S. Moore, “Multiple periodic solutions to a suspension bridge ordinary differential equation,” in Proceedings of the Conference on Nonlinear Differential Equations (Coral Gables, FL, 1999), vol. 5 of Electron. J. Differ. Equ. Conf., pp. 183–199, Southwest Texas State Univ., San Marcos, Tex, USA. View at: Google Scholar
  5. A. C. Lazer and P. J. McKenna, “Large-amplitude periodic oscillations in suspension bridges: some new connections with nonlinear analysis,” SIAM Review, vol. 32, no. 4, pp. 537–578, 1990. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  6. P. J. McKenna, “Large-amplitude periodic oscillations in simple and complex mechanical systems: outgrowths from nonlinear analysis,” Milan Journal of Mathematics, vol. 74, pp. 79–115, 2006. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet

Copyright © 2011 Shanshan Wang and Yukun An. 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.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.