Skip to main content

Steady-state solutions for a suspension bridge with intermediate supports

Abstract

This work is focused on a system of boundary value problems whose solutions represent the equilibria of a bridge suspended by continuously distributed cables and supported by M intermediate piers. The road bed is modeled as the junction of N=M+1 extensible elastic beams which are clamped to each other and pinned at their ends to each pier. The suspending cables are modeled as one-sided springs with stiffness k. Stationary solutions of these doubly nonlinear problems are explicitly and analytically derived for arbitrary k and a general axial load p applied at the ends of the bridge. In particular, we scrutinize the occurrence of buckled solutions in connection with the length of each sub-span of the bridge.

MSC:35G30, 74G05, 74G60, 74K10.

1 Introduction

In this paper, we investigate the solutions of a system of one-dimensional nonlinear problems describing the steady-states of an extensible elastic suspension bridge with M intermediate supports (piers). In particular, we assume that the road bed of the bridge (deck) is composed of N=M+1 extensible elastic beams which are clamped to each other, pinned at their ends to each pier and suspended by a continuous distribution of flexible elastic cables. On account of the midplane stretching of the beams due to their elongations, a geometric nonlinearity appears into the bending equations of the deck.

In the case of a bridge with a single span (no intermediate pier), the problem can be recast into a non-dimensional setting, where its length is supposed to be unitary, for simplicity. Let u=u(x):[0,1]R be a non-dimensional field which accounts for the downward deflection of the bridge in the vertical plane with respect to its reference configuration, and let u + stand for its positive part, namely

u + ={ u if  u 0 , 0 if  u < 0 .

Assuming that both ends of the bridge are pinned, the equation of the bending equilibrium looks like (see [1])

{ u ′′′′ + ( p u L 2 ( 0 , 1 ) 2 ) u + k u + = 0 , u ( 0 ) = u ( 1 ) = u ( 0 ) = u ( 1 ) = 0 .
(1.1)

The term k u + accounts for the restoring force due to one-sided springs which models the supporting cables. Since we confine our attention to stationary conditions, we neglect the dynamical coupling between the deck and the main cable. The constant p represents a non-dimensional measure of the axial force acting at the ends of the span in the reference configuration. Accordingly, p is negative when the span is stretched, positive when compressed. The symbol ′ represents the derivative with respect to the argument.

Assume now that the bridge span is composed of N extensible beams whose internal ends are hinged to M=N1 intermediate piers. Let x i (0,1), i=1,,M, be the points at which these supports are located and let

0= x 0 < x 1 << x N 1 < x N =1.

We denote by I n =( x n 1 , x n ), n=1,,N, the corresponding intervals between two subsequent piers and by | I n | the corresponding length, so that

n = 1 N | I n |=1.

According to the assumptions, the deflection u of the whole span obeys

u( x )= u ( x )=0,=0,1,,N.

Furthermore, we assume that consecutive sub-spans are mutually clamped at the common end, which in turn implies the continuity of u across the supported points x i , i=1,,M (see Figure 1).

Figure 1
figure 1

The joint connecting two consecutive sub-spans and the intermediate pier.

Let u n = u n (x): I n R, n=1,,N, account for the downward deflection of the n th beam in the vertical plane, and let p n be the axial load acting on it. Of course, u n can be viewed as the restriction of u on I n , namely u n =u | I n . Then, taking advantage of the nonlinear model (1.1) for a single supported span, each u n solves the following boundary value problem:

{ u n ′′′′ + ( p n u n L 2 ( I n ) 2 ) u n + k u n + = 0 , in  I n u n ( x n 1 ) = u n ( x n ) = u n ( x n 1 ) = u n ( x n ) = 0 .
(1.2)

Such solutions must fulfill the mutual groove condition

u i ( x i )= u i + 1 ( x i ),i=1,2,,M.
(1.3)

We finally observe that the total elongation of the span equals the sum of the elongations of each beam, namely

u L 2 ( 0 , 1 ) 2 = n = 1 N u n L 2 ( I n ) 2 .

Throughout the paper, we assume a uniform distribution of the axial load, that is, p n =p, n=1,2,,N. In addition, we let

k=κ π 4 ,κ R + ,p=b π 2 ,bR.

Our aim is to scrutinize the existence of suitable buckled solutions for u, which can be obtained by joining buckled solutions for u n on each sub-span, n=1,2,,N. For later convenience, we denote such solutions by u ( N ) =( u 1 ,, u N ). We prove that they exist provided that the lengths of the sub-spans are properly chosen.

1.1 Early contributions

In recent years, an increasing attention has been payed to the analysis of buckling, vibrations and post-buckling dynamics of nonlinear beam models (see, for instance, [2, 3]). As far as we know, most of the papers in the literature deal with approximations and numerical simulations, and only few works are able to derive exact solutions (see [47]).

The investigation of solutions to BVP (1.1), in dependence on p, represents a classical nonlinear buckling problem in the literature on structural mechanics. The notion of buckling, introduced by Euler more than two centuries ago, describes a static instability of structures due to in-plane loading. In this respect, the main concern is to find the critical buckling loads, at which a bifurcation of solutions occurs, and their associated mode shapes, called postbuckling configurations. In the case k=0, a careful analysis of the corresponding buckled stationary states was performed in [7] for all values of p in the presence of a source with a general shape (see also [8]). By replacing u + with u in (1.1), we obtain a simpler model which was scrutinized in [1, 5].

In addition, it is worth noting that (1.1) represents the static counterpart of quite many different evolution problems arising both from elastic, viscoelastic and thermo-elastic theories. An example is the following quasilinear equation:

t t u+ x x x x u+ ( p x u L 2 ( 0 , 1 ) 2 ) x x u+ t u+k u + =f,
(1.4)

which is obtained by matching the modeling of the extensible elastic beam (see [9, 10]) with the well-known equation describing the motion of a damped suspension bridge (see [11, 12])

t t u+ x x x x u+ t u+k u + =f.
(1.5)

Free and forced vibrations of (1.5) were recently scrutinized in [13, 14], whereas the long-time behavior of (1.4) was described in [15] for all values of p.

Obviously, solutions to BVP (1.1) represent the steady states of a lot of models more general than (1.4), for instance, when either the rotational inertia (as in the Kirchhoff theory) or some kind of damping are taken into account. In particular, (1.1) works either when external viscous forces are added or when some structural dissipation phenomena occur in the deck, as in thermoelastic and viscoelastic beams (see, for instance, [1618]).

When the geometric nonlinear term into (1.1) is disregarded, the existence of nontrivial (positive) solutions to the corresponding system were established in [19] by the variational method. Therein, some nonlinearly perturbed versions were also scrutinized, but the set of assumptions made there no longer holds when the full model is considered.

1.2 Outline of the paper

To the best of our knowledge, this is the first paper in the literature dealing with exact solutions to the doubly nonlinear BVPs (1.2), even for N=1. As is well known, the analysis of the corresponding set S κ of their stationary solutions takes a great importance in the longterm dynamics of the corresponding evolution system, especially when its structure is nontrivial [20]. The main results of this paper concern the steady states analysis of a bridge with N=1,2,3 sub-spans and are stated in Section 2. In Section 2.1 we scrutinize the case of a single span without piers (N=1) and we prove that increasing the value of the lateral load p, first a negative u ( 1 ) , then a positive u ( 1 ) + buckled solution appear at equilibrium. In Section 2.2 a bridge with a single pier (N=2) is considered. When the position of the pier is allowed to be asymmetric (| I 1 || I 2 |), we establish a condition on ε=| I 1 | in order that buckled static solutions exist. In particular, we prove that ε1/2 as k0. Taking advantage of these results, the analysis of a bridge with two symmetrically-placed piers (N=3) is performed in Section 2.3, where buckled static solutions are proved to exist provided that | I 1 |=| I 3 || I 2 | fulfills a suitable condition. In Section 3 we deal with the general problem of a bridge with M=N1 piers, and we discuss separately the cases when the number M of the piers is either odd or even. All buckled solutions are determined in a closed form and belong to C 2 (0,1). Each of them is constructed by rescaling and suitably collecting positive and negative solutions, u ( 1 ) + and u ( 1 ) . For any given N, a general explicit formula is established to compute the bifurcation values as a function of k.

2 Stationary states I

2.1 A single span without piers (N=1)

The set S κ of the bridge equilibria in this case consists of all (weak) solutions u ( 1 ) to the following boundary value problem on the interval [0,1] (see Eq. (1.1))

{ u ′′′′ + ( b π 2 u L 2 ( 0 , 1 ) 2 ) u + κ π 4 u + = 0 , u ( 0 ) = u ( 1 ) = u ( 0 ) = u ( 1 ) = 0 .
(2.1)

It is worth noting that S κ is bounded in H 2 (0,1) H 0 1 (0,1) for all bR, κ R + (see [20]).

When κ=0, a general result was established in [7] for a class of non-vanishing sources. In [4], the same strategy with minor modifications was applied to a problem close to (2.1), where the one-sided springs are replaced by unyielding ties. We summarize here the results concerning stationary solutions in the case κ=0.

Theorem 2.1 (see [7], Th. 4.1)

When κ vanishes, the set S 0 is finite for all bR. Letting S b ={jN:b> j 2 } and j =| S b |, S 0 has exactly 2 j +1 elements:

u 0 =0and u j ± (x)= C j ± sinjπx,j S b ,

where

C j ± =± 1 j 2 ( b j 2 ) .

The general case is much more complicated. Since the scheme devised in [5, 7] does not work in the present situation, we obtain here a limited result.

Theorem 2.2 (Existence of buckled solutions)

When κ>0, besides the null solution u 0 , which exists for all bR, problem (2.1) admits

  • a negative buckled solution, u ( 1 ) , ifb>1;

  • a positive buckled solution, u ( 1 ) + , ifb>1+κ.

Proof For all values bR, problem (2.1) has at least the trivial solution, u 0 =0. For all b>1, besides the null solution, we obtain also the negative buckled solution

u ( 1 ) (x)= A 1 sin(πx),with  A 1 ( = C 1 ) = 2 ( b 1 ) ,

which solves the Woinowsky-Krieger problem (see [21])

{ u ′′′′ + ( b π 2 u L 2 ( 0 , 1 ) 2 ) u = 0 , in  [ 0 , 1 ] , u ( 0 ) = u ( 1 ) = u ( 0 ) = u ( 1 ) = 0 .
(2.2)

If b1, there exists only the null solution. Therefore, b=1 will be referred to as the smallest bifurcation value. By paralleling [5], we infer that when b>1+κ, system (2.1) admits the positive buckled solution

u ( 1 ) + (x)= A 1 + sin(πx),with  A 1 + = 2 ( b 1 κ ) .

Since κ>0, we can easily check that u ( 1 ) + (x)<| u ( 1 ) (x)| for all x(0,1) (see Figure 2).  □

Figure 2
figure 2

The buckled solutions u ( 1 ) (a) and u ( 1 ) + (b) when κ=0.1 , b=1.125 .

2.2 A bridge with a single pier (N=2)

Now we consider a bridge with a single pier at x= x 1 and two sub-spans. The whole solution u ( 2 ) =( u 1 , u 2 ) is obtained by joining their deflections, u 1 and u 2 , which solve problems (1.2) for n=1,2 and fulfill the mutual groove condition

u 1 ( x 1 )= u 2 ( x 1 ).
(2.3)

When | I 1 |=| I 2 |, it is easy to check that this condition cannot be satisfied by any buckled solution. Then we choose | I 1 || I 2 |. Let x 1 =ε be the point at which the pier is located, so that | I 1 |=ε. For the sake of definiteness, let u 1 >0 and u 2 <0.

Theorem 2.3 When κ>0, problems (1.2) for N=2 admit two buckled solutions, called u ( 2 ) + and u ( 2 ) , provided that ε= ε κ and b> b κ , where the values of ε κ and b κ are defined in (2.9) and (2.10), respectively.

Proof In order to establish the form of u 1 and u 2 , we rescale the domains [0,ε] and [ε,1] by virtue of the transformations x=εy and x=(1ε)z+ε, respectively. Hence, we get the problem

{ w ′′′′ + ε 2 ( b π 2 1 ε 2 w L 2 ( 0 , 1 ) 2 ) w + ε 4 κ π 4 w = 0 , in  [ 0 , 1 ] , v ′′′′ + ( 1 ε ) 2 ( b π 2 1 ( 1 ε ) 2 v L 2 ( 0 , 1 ) 2 ) v = 0 , in  [ 0 , 1 ] , w ( 0 ) = w ( 1 ) = w ( 0 ) = w ( 1 ) = 0 , v ( 0 ) = v ( 1 ) = v ( 0 ) = v ( 1 ) = 0 ,
(2.4)

where w(y)= u 1 (εy) and v(z)= u 2 ([1ε]z+ε). Obviously, the null solution v=w=0 exists for all bR. On the contrary, nontrivial solutions occur under special conditions.

Arguing as in Section 2.1, we exploit the explicit expression of solutions u ( 1 ) + and u ( 1 ) , respectively, and we obtain

w(y)= A 2 + (ε)sin(πy)andv(z)= A 2 (ε)sin(πz),

where, for some ε κ (0,1),

A 2 + ( ε κ )= ε κ 2 [ b 1 ε κ 2 κ ε κ 2 ] , A 2 ( ε κ )=( ε κ 1) 2 [ b 1 ( 1 ε κ ) 2 ] .

Of course, such solutions both exist provided that

b> b κ =max { 1 ( 1 ε κ ) 2 , 1 + κ ε κ 4 ε κ 2 }

and ε κ satisfies the continuity condition (2.3), namely

w ( 1 ) ε κ = v ( 0 ) 1 ε κ .
(2.5)

Then, by joining w and v, we obtain the whole solution

u ( 2 ) + (x)={ A 2 + ( ε κ ) sin ( π x / ε κ ) , x [ 0 , ε κ ) , 0 , x = ε κ , A 2 ( ε κ ) sin ( π ( x ε κ ) / ( 1 ε κ ) ) , x ( ε κ , 1 ] .
(2.6)

The function u ( 2 ) =( u 2 , u 1 ), which is obtained by symmetrizing u ( 2 ) + with respect to the vertical line x=1/2, solves (1.2) under the same conditions (see Figure 3). Its expression is given by replacing ε κ with 1 ε κ and A 2 + with A 2 into (2.6),

u ( 2 ) (x)={ A 2 ( ε κ ) sin ( π x / ( 1 ε κ ) ) , x [ 0 , 1 ε κ ) , 0 , x = 1 ε κ , A 2 + ( ε κ ) sin ( π ( x 1 + ε κ ) / ε κ ) , x ( 1 ε κ , 1 ] .
(2.7)
Figure 3
figure 3

The buckled solutions u ( 2 ) + (a) and u ( 2 ) (b) when κ=5.9 , b=5.75 .

Computation of ε κ and b κ . For each κ>0, in order to compute the unknown value ε κ , which ensures the existence of buckled solutions, we exploit the continuity condition (2.5), namely

(1 ε κ ) A 2 + ( ε κ )+ ε κ A 2 ( ε κ )=0.

Thus, the required value ε κ has to satisfy the system

{ 1 ε κ 2 + κ ε κ 2 = 1 ( 1 ε κ ) 2 , 0 < ε κ < 1 ,
(2.8)

which implies

κ= 2 ε κ 1 ε κ 4 ( 1 ε κ ) 2 .
(2.9)

It can be easily checked (see Figure 4) that ε κ (1/2,1) and

lim κ 0 ε κ ε 0 = 1 2 , lim κ + ε κ =1.

In order to compute the critical value of the axial force, b κ , we observe that by virtue of (2.9), we have

b κ = 1 ( 1 ε κ ) 2 = 1 + κ ε κ 4 ε κ 2 .
(2.10)
Figure 4
figure 4

The value of ε κ as a function of κ .

Summarizing, for any given κ>0, if b b κ no buckled solution exists, if b> b κ there exists a unique value of ε κ , which produces two buckled solutions. As a consequence, b κ represents the bifurcation value. According to (2.10), we get

A 2 + ( ε κ )= ε κ 2 ( b b κ ) , A 2 ( ε κ )=( ε κ 1) 2 ( b b κ )

and A 2 < A 2 + . In addition, for all κ>0, we can easily obtain

b κ > b 0 =4.

Then, when b4, the set S κ contains only the trivial solution u 0 . □

In the limit κ0, it follows that

A 2 , A 2 + ( b 4 ) / 2 .

This means that solutions u ( 2 ) + and u ( 2 ) tend to coincide with the second bifurcation branch of problem (2.2), as expected.

2.3 A bridge with two piers (N=3)

When the bridge has two intermediate piers and three sub-spans, we shall construct solutions u ( 3 ) =( u 1 , u 2 , u 3 ), where u 1 , u 2 and u 3 solve (1.2). It is easy to check that no buckled solution exists when | I 1 |=| I 2 |=| I 3 |. Then we choose | I 1 |=| I 3 || I 2 | and construct a buckled solution by joining three (suitably rescaled) functions which have the form of either u ( 1 ) + or u ( 1 ) .

Theorem 2.4 When κ>0, problems (1.2) for N=3 admit two buckled solutions, u ( 3 ) + and u ( 3 ) , provided that b> c κ + and b> c κ , respectively, where the values of c κ + and c κ are defined in (2.12) and (2.15).

Proof The buckled solutions are constructed as follows. The former, u ( 3 ) + , has two positive components with the same shape on I 1 and I 3 , and a negative one on I 2 . The latter, u ( 3 ) , is negative in I 1 and I 3 , but positive in I 2 . Accordingly, let

w j = u ( 3 ) + | I 2 j 1 , j = 1 , 2 and v 0 = u ( 3 ) + | I 2 , v j = u ( 3 ) | I 2 j 1 , j = 1 , 2 and w 0 = u ( 3 ) | I 2 ,

so that w k >0 and v k <0 for k=0,1,2.

In order to construct u ( 3 ) + , we exploit the same argument as in Section 2.2. First, we argue that

| I 1 |=| I 3 |=ξ,| I 2 |=12ξ,ξ(0,1/2).

Then, we apply the rescaling procedure with the scale factor ξ/(1ξ)=ε and, by virtue of the transformations x=ξy and x=(12ξ)z+ξ, we obtain

{ w j ′′′′ + ξ 2 ( b π 2 1 ξ 2 w j L 2 ( 0 , 1 ) 2 ) w j + κ ξ 4 π 4 w j = 0 , in  [ 0 , 1 ] , j = 1 , 2 , v 0 ′′′′ + ( 1 2 ξ ) 2 ( b π 2 1 ( 1 2 ξ ) 2 v 0 L 2 ( 0 , 1 ) 2 ) v 0 = 0 , in  [ 0 , 1 ] , w j ( 0 ) = w j ( 1 ) = w j ( 0 ) = w j ( 1 ) = 0 , j = 1 , 2 , v 0 ( 0 ) = v 0 ( 1 ) = v 0 ( 0 ) = v 0 ( 1 ) = 0 .
(2.11)

Taking into account that ε= ε κ is uniquely and implicitly given by (2.9), the corresponding value of ξ becomes

ξ κ = ε κ 1 + ε κ ,

so that

lim κ 0 ξ κ = ξ 0 = 1 3 , lim κ + ξ κ = 1 2 .

Accordingly, we obtain

w j (y)= A 3 + ( ξ κ )sin(πy),j=1,2, v 0 (z)= A 3 ( ξ κ )sin(πz),

where

A 3 + ( ξ κ )= ξ κ 2 [ b 1 ξ κ 2 κ ξ κ 2 ] , A 3 ( ξ κ )=(2 ξ κ 1) 2 [ b 1 ( 1 2 ξ κ ) 2 ] .

Of course, such solutions exist at the same time, provided that b> c κ + , where

c κ + =max { 1 ( 1 2 ξ κ ) 2 , 1 + κ ξ κ 4 ξ κ 2 } =max { ( 1 + ε κ ) 2 ( 1 ε κ ) 2 , ( 1 + ε κ ) 4 + κ ε κ 4 ε κ 2 ( 1 + ε κ ) 2 } ,

where c κ + represents the bifurcation value. By virtue of (2.9), we get

c κ + = ( 1 + ε κ ) 2 ( 1 ε κ ) 2 .
(2.12)

If this is the case, we havethe whole solution (see Figure 5a)

u ( 3 ) + (x)={ A 3 + ( ξ κ ) sin ( π x / ξ κ ) , x [ 0 , ξ κ ] , A 3 ( ξ κ ) sin ( π ( x ξ κ ) / ( 1 2 ξ κ ) ) , x [ ξ κ , 1 ξ κ ] , A 3 + ( ξ κ ) sin ( π ( x 1 + ξ κ ) / ξ κ ) , x [ 1 ξ κ , 1 ] .
(2.13)

In order to construct u ( 3 ) , we proceed as before by letting

| I 1 |=| I 3 |=η,| I 2 |=12η,η(0,1/2).

By means of a rescaling procedure with the scale factor η/(1η)=1ε, the transformations x=ηz and x=(12η)y+η lead to

{ v j ′′′′ + η 2 ( b π 2 1 η 2 v j L 2 ( 0 , 1 ) 2 ) v j = 0 , in  [ 0 , 1 ] , j = 1 , 2 , w 0 ′′′′ + ( 1 2 η ) 2 ( b π 2 1 ( 1 2 η ) 2 w 0 L 2 ( 0 , 1 ) 2 ) w 0 + κ ( 1 2 η ) 4 π 4 w 0 = 0 , in  [ 0 , 1 ] , v j ( 0 ) = v j ( 1 ) = v j ( 0 ) = v j ( 1 ) = 0 , j = 1 , 2 , w 0 ( 0 ) = w 0 ( 1 ) = w 0 ( 0 ) = w 0 ( 1 ) = 0 .
(2.14)

Accordingly, we find

η κ = 1 ε κ 2 ε κ ,

where ε κ is given by (2.9), and

lim κ 0 η κ = η 0 = 1 3 , lim κ + η κ =0.

So, we obtain

v j (z)= A 3 ( η κ )sin(πz),j=1,2, w 0 (y)= A 3 + ( η κ )sin(πy),

where

A 3 ( η κ ) = η κ 2 [ b 1 η κ 2 ] , A 3 + ( η κ ) = ( 1 2 η κ ) 2 [ b 1 ( 1 2 η κ ) 2 κ ( 1 2 η κ ) 2 ] .

Of course, such solutions do exist together provided that

b> c κ =max { 1 η κ 2 , 1 + κ ( 1 2 η κ ) 4 ( 1 2 η κ ) 2 } =max { ( 2 ε κ ) 2 ( 1 ε κ ) 2 , ( 2 ε κ ) 4 + κ ε κ 4 ε κ 2 ( 2 ε κ ) 2 } ,

where c κ represents the bifurcation value. By virtue of (2.9), we get

c κ = ( 2 ε κ ) 2 ( 1 ε κ ) 2 .
(2.15)

By means of (2.9) it is easy to check that c κ < c κ + . The complete solution is given by (see Figure 5b)

u ( 3 ) (x)={ A 3 ( η κ ) sin ( π x / η κ ) , x [ 0 , η κ ] , A 3 + ( η κ ) sin ( π ( x η κ ) / ( 1 2 η κ ) ) , x [ η κ , 1 η κ ] , A 3 ( η κ ) sin ( π ( x 1 + η κ ) / η κ ) , x [ 1 η κ , 1 ] .
(2.16)

 □

Figure 5
figure 5

The buckled solutions u ( 3 ) + (a) and u ( 3 ) (b) when κ=0.1 and b=1.225 .

3 Stationary states II

In this section we generalize the problem to a bridge with N sub-spans and M=N1 piers. The existence of buckled solutions is investigated in connection with the length of the sub-spans. Indeed, it is easy to check that no buckled solution exists when all of them are of the same length. Then a buckled solution may be obtained by collecting and joining N (suitably rescaled) functions of the same form as either u ( 1 ) + or u ( 1 ) . To this end, we are forced to consider separately the cases when the number M of the piers is either odd or even. In the former case, indeed, we adopt a strategy which is close to that applied in Section 2.2. In the latter, we iterate the procedure devised in Section 2.3.

Theorem 3.1 For any κ>0, NN, (1.2) admits two buckled solutions.

  • In the odd case, M=2m1, there exist u ( 2 m ) + and u ( 2 m ) provided thatb> b κ (m), where the value of b κ (m)is characterized in (3.5).

  • In the even case, M=2m, there exist u ( 2 m + 1 ) + and u ( 2 m + 1 ) provided thatb> c κ + (m)andb> c κ (m), respectively, where the values of c κ + (m)and c κ (m)are characterized in (3.6).

Proof The odd case. Solutions u ( 2 m ) + and u ( 2 m ) , m=(M+1)/2=N/2, are assumed to change the sign alternately on the sub-spans. The superscript + (−) means that the solution is positive (negative) on I 1 . In order to construct them, we join 2m rescaled functions of the same form as either u ( 1 ) + or u ( 1 ) . Arguing as in the step N=2, we let

| I 2 j 1 |= ε κ m ,j=1,,m,| I 2 h + 2 |= 1 ε κ m ,h=0,,m1,

where ε κ is given by (2.9).

Construction of u ( 2 m ) + . Let

w j = u ( 2 m ) + | I 2 j 1 ,j=1,,m, v h = u ( 2 m ) + | I 2 h + 2 ,h=0,,m1,

so that w j >0 and v h <0. Since

| I 2 j 1 I 2 j |= 1 m ,

we stress that each restriction

u ( 2 m ) + | I 2 j 1 I 2 j ,j=1,,m,

is similar to u ( 2 ) + , rescaled by a factor 1/m. In particular,

I 2 j 1 = ( j 1 m , j 1 + ε κ m ) , j = 1 , , m , I 2 h + 2 = ( h + ε κ m , h + 1 m ) , h = 0 , , m 1 .

Therefore, on I 2 j 1 we have

w j (x)= A 2 m + sin ( π ( x j 1 m ) / ε κ m ) ,j=1,,m,

while on I 2 h + 2 we obtain

v h (x)= A 2 m sin ( π ( x h + ε κ m ) / 1 ε κ m ) ,h=0,,m1,

where

A 2 m + = ε κ m 2 [ b m 2 ε κ 2 κ ε κ 2 m 2 ] , A 2 m = ε κ 1 m 2 [ b m 2 ( 1 ε κ ) 2 ] .
(3.1)

Such solutions exist provided that

b> b κ (m)=max { m 2 ε κ 2 + κ ε κ 2 m 2 , m 2 ( 1 ε κ ) 2 } .
(3.2)

By virtue of Lemma 3.2, which will be proved later, we get

b κ (m)= m 2 b κ ,m=N/2.

Construction of u ( 2 m ) . Let

v j = u ( 2 m ) | I 2 j 1 ,j=1,,m, w h = u ( 2 m ) | I 2 h + 2 ,h=0,,m1,

so that w h >0 and v j <0. Arguing as before, we obtain

| I 2 j 1 |= 1 ε κ m ,j=1,,m,| I 2 h + 2 |= ε κ m ,h=0,,m1,

where ε κ is given by (2.9) and

u ( 2 m ) | I 2 j 1 I 2 j ,j=1,,m,

is similar to u ( 2 ) , rescaled by a factor 1/m. In particular,

I 2 j 1 = ( j 1 m , j ε κ m ) , j = 1 , , m , I 2 h + 2 = ( h + 1 ε κ m , h + 1 m ) , h = 0 , , m 1 ,

so on I 2 j 1 we have

v j (x)= A 2 m sin ( π ( x j 1 m ) / 1 ε κ m ) ,j=1,,m,

while on I 2 h + 2 we obtain

w h (x)= A 2 m + sin ( π ( x h + 1 ε κ m ) / ε κ m ) ,h=0,,m1,

with A 2 m + , A 2 m defined as in (3.1).

The even case. In this case, we construct the solutions u ( 2 m + 1 ) + and u ( 2 m + 1 ) , where m=M/2=(N1)/2.

Construction of u ( 2 m + 1 ) + . Let

w j = u ( 2 m + 1 ) + | I 2 j 1 , j = 1 , , m + 1 , v h = u ( 2 m + 1 ) + | I 2 h + 2 , h = 0 , , m 1 ,

so that w j >0 and v h <0. Arguing as in the step N=3, we obtain

| I 2 j 1 | = ξ κ , j = 1 , , m + 1 , | I 2 h + 2 | = 1 ( 1 + m ) ξ κ m , h = 0 , , m 1 ,

where ξ κ is given by ξ κ = ε κ m + ε κ . Therefore

| I 2 j 1 | = ε κ m + ε κ , | I 2 h + 2 | = 1 ε κ m + ε κ

and

I 2 j 1 = ( j 1 m + ε κ , j 1 + ε κ m + ε κ ) , j = 1 , , m + 1 , I 2 h + 2 = ( h + ε κ m + ε κ , h + 1 m + ε κ ) , h = 0 , , m 1 .

So, on I 2 j 1 we have

w j (x)= A 2 m + 1 + sin ( π ( x j 1 m + ε κ ) / ε κ m + ε κ ) ,j=1,,m+1,

while on I 2 h + 2 we obtain

v h (x)= A 2 m + 1 sin ( π ( x h + ε κ m + ε κ ) / 1 ε κ m + ε κ ) ,h=0,,m1,

where

A 2 m + 1 + = ε κ m + ε κ 2 [ b ( m + ε κ ) 2 ε κ 2 κ ε κ 2 ( m + ε κ ) 2 ] , A 2 m + 1 = ε κ 1 m + ε κ 2 [ b ( m + ε κ ) 2 ( 1 ε κ ) 2 ] .

Such solutions exist provided that

b> c κ + (m)=max { ( m + ε κ ) 2 ε κ 2 + κ ε κ 2 ( m + ε κ ) 2 , ( m + ε κ ) 2 ( 1 ε κ ) 2 } .
(3.3)

By virtue of Lemma 3.2, we get

c κ + (m)= ( m + ε κ ) 2 ( 1 ε κ ) 2 ,m=(N1)/2.

Construction of u ( 2 m + 1 ) . Let

v j = u ( 2 m + 1 ) | I 2 j 1 ,j=1,,m+1, w h = u ( 2 m + 1 ) | I 2 h + 2 ,h=0,,m1,

so that w h >0 and v j <0. Arguing as in the step N=3, we obtain

| I 2 j 1 | = η κ , j = 1 , , m + 1 , | I 2 h + 2 | = 1 ( 1 + m ) η κ m , h = 0 , , m 1 ,

where η κ is given by η κ = 1 ε κ m + 1 ε κ . Therefore

| I 2 j 1 |= 1 ε κ m + 1 ε κ ,| I 2 h + 2 |= ε κ m + 1 ε κ

and

I 2 j 1 = ( j 1 m + 1 ε κ , j ε κ m + 1 ε κ ) , j = 1 , , m + 1 , I 2 h + 2 = ( h + 1 ε κ m + 1 ε κ , h + 1 m + 1 ε κ ) , h = 0 , , m 1 .

As a consequence, on I 2 j 1 we have

v j (x)= B 2 m + 1 sin ( π ( x j 1 m + 1 ε κ ) / 1 ε κ m + 1 ε κ ) ,j=1,,m+1,

while on I 2 h + 2 we obtain

w h (x)= B 2 m + 1 + sin ( π ( x h + 1 ε κ m + 1 ε κ ) / ε κ m + 1 ε κ ) ,h=0,,m1,

where

B 2 m + 1 + = ε κ m + 1 ε κ 2 [ b ( m + 1 ε κ ) 2 ε κ 2 κ ε κ 2 ( m + 1 ε κ ) 2 ] , B 2 m + 1 = ε κ 1 m + 1 ε κ 2 [ b ( m + 1 ε κ ) 2 ( 1 ε κ ) 2 ] .

Such solutions exist provided that b> c κ (m), where

c κ (m)=max { ( m + 1 ε κ ) 2 ε κ 2 + κ ε κ 2 ( m + 1 ε κ ) 2 , ( m + 1 ε κ ) 2 ( 1 ε κ ) 2 } .
(3.4)

By virtue of Lemma 3.2, we get

c κ (m)= ( m + 1 ε κ ) 2 ( 1 ε κ ) 2 ,m=(N1)/2.

 □

Lemma 3.2 (Characterization of the bifurcation values)

For any given κ>0 and mN,

b κ (m)= m 2 ( 1 ε κ ) 2 = m 2 b κ ,
(3.5)
c κ (m)= ( m + 1 ε κ ) 2 ( 1 ε κ ) 2 < c κ + (m)= ( m + ε κ ) 2 ( 1 ε κ ) 2 ,
(3.6)

where ε κ and b κ are computed by (2.9) and (2.10), respectively.

Proof First, in view of (3.2), we have to prove that

m 2 ( 1 ε κ ) 2 m 4 + κ ε κ 4 ε κ 2 ,

which is equivalent to

κ m 4 2 ε κ 1 ε κ 4 ( 1 ε κ ) 2 .

By replacing the expression of κ given by (2.9), we obtain

1 m 4 ,

which implies that (3.5) holds for all mN. Then, in view of (3.4), we have to prove

( m + 1 ε κ ) 2 ( 1 ε κ ) 2 ( m + 1 ε κ ) 4 + κ ε κ 4 ε κ ( m + 1 ε κ ) 2 ,

which is equivalent to

κ ( 2 ε κ 1 ) ( m + 1 ε κ ) 4 ε κ 4 ( 1 ε κ ) 2 .

By virtue of (2.9), we obtain

1 ( m + 1 ε κ ) 4 ,

which is identically satisfied for all admissible values ε κ and for all mN. Hence, (3.6)1 holds. Finally, in view of (3.3), we need to prove

( m + ε κ ) 2 ( 1 ε κ ) 2 ( m + ε κ ) 4 + κ ε κ 4 ε κ ( m + ε κ ) 2 ,

and this is equivalent to

κ ( 2 ε κ 1 ) ( m + ε κ ) 4 ε κ 4 ( 1 ε κ ) 2 .

Applying (2.9) as in the previous cases, we obtain 1 ( m + ε κ ) 4 , which is identically satisfied for all admissible values ε κ and for all mN so that (3.6)2 follows. Moreover, starting from (3.6) and taking into account that 1/2< ε κ <1, it is easy to check that

c κ (m)< c κ + (m).

 □

Previous results can be summarized by means of a simple sketch which highlights the main bifurcation features (see Figure 6).

Figure 6
figure 6

The bifurcation portrait: for N odd (on the left, m=N/2 ) and N even (on the right, m=(N1)/2 ).

Finally, it is worth noting that by removing the coupling between the road-bed and the cable, we recover well-known results (see, for instance, [5, 8])

Remark 3.3 In the limit κ0, from (2.9) it follows ε κ 1/2, so that buckled solutions exist even if all sub-spans are equal. Moreover,

b κ (m) ( 2 m ) 2 = N 2 , c κ (m), c κ + (m) ( 2 m + 1 ) 2 = N 2 .

References

  1. Bochicchio I, Giorgi C, Vuk E: On some nonlinear models for suspension bridges. In Evolution Equations and Materials with Memory, Proceedings. Edited by: Andreucci D, Carillo S, Fabrizio M, Loreti P, Sforza D. Università La Sapienza Ed., Roma; 2012:1-17. Rome, 2010

    Google Scholar 

  2. Choo SM, Chung SK: Finite difference approximate solutions for the strongly damped extensible beam equations. Appl. Math. Comput. 2000, 112: 11-32. 10.1016/S0096-3003(99)00005-3

    Article  MathSciNet  MATH  Google Scholar 

  3. Fang W, Wickert JA: Postbuckling of micromachined beams. J. Micromech. Microeng. 1994, 4: 116-122. 10.1088/0960-1317/4/3/004

    Article  Google Scholar 

  4. Bochicchio I, Giorgi C, Vuk E: Steady states analysis and exponential stability of an extensible thermoelastic system. Comun. SIMAI Congr. 2009, 3: 232-244.

    Google Scholar 

  5. Bochicchio I, Vuk E: Buckling and longterm dynamics of a nonlinear model for the extensible beam. Math. Comput. Model. 2010, 51: 833-846. 10.1016/j.mcm.2009.10.010

    Article  MathSciNet  MATH  Google Scholar 

  6. Chen JS, Ro WC, Lin JS: Exact static and dynamic critical loads of a sinusoidal arch under a point force at the midpoint. Int. J. Non-Linear Mech. 2009, 44: 66-70. 10.1016/j.ijnonlinmec.2008.08.006

    Article  Google Scholar 

  7. Coti Zelati M, Giorgi C, Pata V: Steady states of the hinged extensible beam with external load. Math. Models Methods Appl. Sci. 2010, 20: 43-58. 10.1142/S0218202510004143

    Article  MathSciNet  MATH  Google Scholar 

  8. Reiss EL, Matkowsky BJ: Nonlinear dynamic buckling of a compressed elastic column. Q. Appl. Math. 1971, 29: 245-260.

    MathSciNet  MATH  Google Scholar 

  9. Ball JM: Initial-boundary value problems for an extensible beam. J. Math. Anal. Appl. 1973, 42: 61-90. 10.1016/0022-247X(73)90121-2

    Article  MathSciNet  MATH  Google Scholar 

  10. Woinowsky-Krieger S: The effect of an axial force on the vibration of hinged bars. J. Appl. Mech. 1950, 17: 35-36.

    MathSciNet  MATH  Google Scholar 

  11. Lazer AC, McKenna PJ: Large-amplitude periodic oscillations in suspension bridges: some new connections with nonlinear analysis. SIAM Rev. 1990, 32: 537-578. 10.1137/1032120

    Article  MathSciNet  MATH  Google Scholar 

  12. McKenna PJ, Walter W: Nonlinear oscillations in a suspension bridge. Arch. Ration. Mech. Anal. 1987, 98: 167-177.

    Article  MathSciNet  MATH  Google Scholar 

  13. Abdel-Ghaffar AM, Rubin LI: Non linear free vibrations of suspension bridges: theory. J. Eng. Mech. 1983, 109: 313-329. 10.1061/(ASCE)0733-9399(1983)109:1(313)

    Article  Google Scholar 

  14. Ahmed NU, Harbi H: Mathematical analysis of dynamic models of suspension bridges. SIAM J. Appl. Math. 1998, 58: 853-874. 10.1137/S0036139996308698

    Article  MathSciNet  MATH  Google Scholar 

  15. Bochicchio I, Giorgi C, Vuk E: Longterm damped dynamics of the extensible suspension bridge. Int. J. Differ. Equ. 2010., 2010: Article ID 383420

    Google Scholar 

  16. Bochicchio I, Vuk E: Longtime behavior for oscillations of an extensible viscoelastic beam with elastic external supply. Int. J. Pure Appl. Math. 2010, 58: 61-67.

    MathSciNet  MATH  Google Scholar 

  17. Giorgi C, Naso MG: Modeling and steady states analysis of the extensible thermoelastic beam. Math. Comput. Model. 2011, 53: 896-908. 10.1016/j.mcm.2010.10.026

    Article  MathSciNet  MATH  Google Scholar 

  18. Giorgi C, Pata V, Vuk E: On the extensible viscoelastic beam. Nonlinearity 2008, 21: 713-733. 10.1088/0951-7715/21/4/004

    Article  MathSciNet  MATH  Google Scholar 

  19. An Y, Fan X: On the coupled systems of second and fourth order elliptic equations. Appl. Math. Comput. 2003, 140: 341-351. 10.1016/S0096-3003(02)00231-X

    Article  MathSciNet  MATH  Google Scholar 

  20. Bochicchio I, Giorgi C, Vuk E: Long-term dynamics of the coupled suspension bridge system. Math. Models Methods Appl. Sci. 2012., 22: Article ID 1250021

    Google Scholar 

  21. Dickey RW: Dynamic stability of equilibrium states of the extensible beam. Proc. Am. Math. Soc. 1973, 41: 94-102. 10.1090/S0002-9939-1973-0328290-8

    Article  MathSciNet  MATH  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Claudio Giorgi.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

CG conceived the study, participated in its design, verified all calculations and drew all figures, EV participated in the design of the study, performed the proofs of theorems and lemmas, carried out all calculations and drafted the manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Giorgi, C., Vuk, E. Steady-state solutions for a suspension bridge with intermediate supports. Bound Value Probl 2013, 204 (2013). https://doi.org/10.1186/1687-2770-2013-204

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2013-204

Keywords