Research article

Existence result for a nonlinear nonlocal system modeling suspension bridges

  • Received: 07 September 2018 Accepted: 28 November 2018 Published: 04 December 2018
  • MSC : 35G61, 74B20

  • A nonlinear nonlocal partial di erential system modeling suspension bridge is considered. We analyze the well-posedness of the "hyperbolic" type system through a Galerkin procedure. A correspond linear problem admits a unique solution, which makes us find that the original system also has a solution with high regularity.

    Citation: Yongda Wang. Existence result for a nonlinear nonlocal system modeling suspension bridges[J]. AIMS Mathematics, 2018, 3(4): 608-624. doi: 10.3934/Math.2018.4.608

    Related Papers:

    [1] Zayd Hajjej . A suspension bridges with a fractional time delay: Asymptotic behavior and Blow-up in finite time. AIMS Mathematics, 2024, 9(8): 22022-22040. doi: 10.3934/math.20241070
    [2] Xing Zhang, Lei Liu, Yan-Jun Liu . Adaptive NN control based on Butterworth low-pass filter for quarter active suspension systems with actuator failure. AIMS Mathematics, 2021, 6(1): 754-771. doi: 10.3934/math.2021046
    [3] Zayd Hajjej . On the exponential decay of a Balakrishnan-Taylor plate with strong damping. AIMS Mathematics, 2024, 9(6): 14026-14042. doi: 10.3934/math.2024682
    [4] Tamer Nabil . Ulam stabilities of nonlinear coupled system of fractional differential equations including generalized Caputo fractional derivative. AIMS Mathematics, 2021, 6(5): 5088-5105. doi: 10.3934/math.2021301
    [5] Yaqiong Liu, Yunting Li, Qiuping Liao, Yunhui Yi . Classification of nonnegative solutions to fractional Schrödinger-Hatree-Maxwell type system. AIMS Mathematics, 2021, 6(12): 13665-13688. doi: 10.3934/math.2021794
    [6] Bashir Ahmad, Ahmed Alsaedi, Areej S. Aljahdali, Sotiris K. Ntouyas . A study of coupled nonlinear generalized fractional differential equations with coupled nonlocal multipoint Riemann-Stieltjes and generalized fractional integral boundary conditions. AIMS Mathematics, 2024, 9(1): 1576-1594. doi: 10.3934/math.2024078
    [7] Bashir Ahmad, Ahmed Alsaedi, Mona Alsulami, Sotiris K. Ntouyas . Existence theory for coupled nonlinear third-order ordinary differential equations with nonlocal multi-point anti-periodic type boundary conditions on an arbitrary domain. AIMS Mathematics, 2019, 4(6): 1634-1663. doi: 10.3934/math.2019.6.1634
    [8] Ahmed Alsaedi, Fawziah M. Alotaibi, Bashir Ahmad . Analysis of nonlinear coupled Caputo fractional differential equations with boundary conditions in terms of sum and difference of the governing functions. AIMS Mathematics, 2022, 7(5): 8314-8329. doi: 10.3934/math.2022463
    [9] Madeaha Alghanmi, Shahad Alqurayqiri . Existence results for a coupled system of nonlinear fractional functional differential equations with infinite delay and nonlocal integral boundary conditions. AIMS Mathematics, 2024, 9(6): 15040-15059. doi: 10.3934/math.2024729
    [10] Ahmed M.A. El-Sayed, Eman M.A. Hamdallah, Hameda M. A. Alama . Multiple solutions of a Sturm-Liouville boundary value problem of nonlinear differential inclusion with nonlocal integral conditions. AIMS Mathematics, 2022, 7(6): 11150-11164. doi: 10.3934/math.2022624
  • A nonlinear nonlocal partial di erential system modeling suspension bridge is considered. We analyze the well-posedness of the "hyperbolic" type system through a Galerkin procedure. A correspond linear problem admits a unique solution, which makes us find that the original system also has a solution with high regularity.


    1. Introduction

    The present paper is devoted to the following differential system, which is composed by a fourth order differential equation and a second order one. Namely,

    {(M+2mξ)ytt=EIyxxxx+H0(2yxξ2+3s(y2x+2θ2x)ξ4)xEALc[L0y2x+2θ2xξ3]sξ3                             2EALc[L0syξ3](sξyxξ3)x+2EA2Lc[L0sθξ3](θxξ3)x,x(0,L),t>0,(M3+2mξ)θtt=GK2θxx+2H0(θxξ2+3syxθxξ4)x2EALc[L0yxθxξ3]sξ3                             2EALc[L0sθξ3](sξyxξ3)x+2EALc[L0syξ3](θxξ3)x,   x(0,L),t>0. (1.1)

    The system (1.1), suggested by Arioli-Gazzola [2], aims to model the main span of a suspension bridge. In [2], the authors considered torsional instability through several numerical results, see [8,9] for some more recent results on torsional instability. Here we give a simple sketch for suspension bridges in Figure 1. The deck, which is viewed as a thin rectangular plate, is hooked to a large number of hangers. At their upper endpoints, the hangers are connected with the main cables, which are sustained by the towers. We assume that the origin of the coordinate system is at the point O and that the positive vertical axis is oriental downwards, see Fig. 1.

    Figure 1. Sketch of a suspension bridge.

    In system (1.1), L is the length of the deck between the two towers and 2 is its width, y=y(x,t) is the vertical (downwards) displacement of the longitudinal midline of the deck, θ=θ(x,t) is the torsional angle of the deck. The position of the cables at rest is denoted by s=s(x)<0 (the positive direction is vertical downwards). Hence, the function

    ξ=ξ(x):=1+s(x)2

    represents the length of the cable at rest. For the other parameters M,m,E,I,H0,Lc,A,G,K which are positive constants, one can find their physical meanings in [2,Section 4.4].

    As in Arioli-Gazzola [1], the degenerate plate (modeling the deck of the main span of the bridge) is assumed to be hinged at its two short edges where the cross sections between the towers cannot rotate. Hence, the boundary conditions associated to (1.1) are

    y(0,t)=y(L,t)=yxx(0,t)=yxx(L,t)=θ(0,t)=θ(L,t)=0,t>0. (1.2)

    In order to consider an initial value problem, we introduce the following conditions:

    y(x,0)=y0,yt(x,0)=y1,θ(x,0)=θ0,θt(x,0)=θ1,x(0,L). (1.3)

    The system (1.1) is suggested, recently, by Arioli-Gazzola [2] for describing the dynamics of a suspension bridge. In that article, they made several numerical experiments for the nonlinear model (1.1) and gave a new explanation for the torsional instability in suspension bridge, but they did not fully study the system (1.1) from the theoretical point of view. Hence, the aim of this paper is to solve the initial-boundary value problem (1.1)-(1.2)-(1.3) theoretically. For simplicity, we write them as follows

    {ytt+Φyxxxx=h1(y,θ),x(0,L),t>0,θttΨy(x,t)θxx=h2(y,θ),x(0,L),t>0,y(0,t)=y(L,t)=yxx(0,t)=yxx(L,t)=θ(0,t)=θ(L,t)=0,t>0,y(x,0)=y0,yt(x,0)=y1,θ(x,0)=θ0,θt(x,0)=θ1,x(0,L), (1.4)

    where Φ=Φ(x)=EIM+2mξ, Ψy(x,t)=(M3+2mξ)1(GK2+2H0ξ2+6H0syxξ4+2AEξ3LcL0syξ3) and the two operators h1(,), h2(,): H4(0,L)×H3(0,L)L2(0,L) are given by

    h1(y,θ)=H0M+2mξ(2yxξ2+3s(y2x+2θ2x)ξ4)xEA(M+2mξ)Lc[L0y2x+2θ2xξ3]sξ32EA(M+2mξ)Lc[L0syξ3](sξyxξ3)x+2EA2(M+2mξ)Lc[L0sθξ3](θxξ3)x,h2(y,θ)=(M3+2mξ)1[(2H0ξ2+6H0syxξ4+2AEξ3Lc[L0syξ3])xθx2EALc[L0yxθxξ3]sξ3](M3+2mξ)12EALc[L0sθξ3](sξyxξ3)x.

    The second order differential equation in (1.4) is of mixed type and the unpleasant terms h1(y,θ), h2(y,θ) are nonlinear and nonlocal. Moreover, the term h1(y,θ) contains the second order derivatives of θ in x whereas h2(y,θ) includes the couplings such as yxθx and the second order derivatives of y in x. These peculiarities show that (1.4) is a nonstandard problem and it appears quite challenging to prove existence of a solution. Here we restrict to the case where the second order equation in (1.4) is a "hyperbolic" type one, i.e. Ψy(x,t)>0 and then we investigate the system (1.4) by the Galerkin method as the researchers study this type of (nonlinear) hyperbolic equations or systems, see for instance [1,3,10,11,12] and references therein. However, note that the complicated terms h1(y,θ) and h2(y,θ) include the second order derivatives of y,θ, we may obtain solutions by the Galerkin method only if the initial data in higher regularity functional spaces such as H4×H3. Moreover, the terms h1(y,θ) and h2(y,θ) are somewhat different from the cases, such as Ball [3] or Temam [11] and references therein we met before, it seems out of reach to obtain existence results by applying their methods directly.

    In this paper we follow a different path and tackle the problem (1.4) in some high regularity functional spaces. We first solve a correspond linear problem and we obtain existence and uniqueness results. Then we show that the original system admits a local solution by applying the Contraction Mapping principle and the iterative method.

    Let

    H2=H2(0,L):={uH2(0,L);u(x)=0 on {0,L}},H3=H3(0,L):={uH3(0,L);u(x)=u(x)=0 on {0,L}},H4=H4(0,L):={uH4(0,L);u(x)=u(x)=0 on {0,L}},

    which are complete subspaces of the Hilbert spaces H2(0,L), H3(0,L) and H4(0,L) respectively. Denote

    Y:=C0([0,T];H4(0,L))C1([0,T];H2(0,L))C2([0,T];L2(0,L));
    Θ:=C0([0,T];H3(0,L))C1([0,T];H2(0,L))C2([0,T];L2(0,L)).

    We say that (y,θ)Y×Θ is a strong solution of (1.4) if it fulfills the initial conditions and if the following equations are satisfied for every t[0,T]:

    ytt+Φyxxxx=h1(y,θ),θttΨy(x,t)θxx=h2(y,θ).

    Then we prove

    Theorem 1.1. Assume that sC4(0,L)C3([0,L]). Let y0H4, θ0H3 and y1,θ1H2. Then for a small T>0 there exists a strong solution (y,θ) of (1.4). Moreover,

    θC2([0,T];H10(0,L)).

    Remark 1.2. The component θ of (y,θ) is a classical solution to the second order equation.

    By the Sobolev embedding theorem, we have θxxC0([0,T]×[0,L]), and then θttC0([0,T]×[0,L]). Hence, θC2([0,L]×[0,T]).

    Remark 1.3. The existence of global solutions of (1.4) is still open.

    According to [2,Section 2.2], where the authors analyzed all the energies of the bridge, we can write the total energy of the structure for t0 as follows

    E(t)=M2L0[y2t+2θ2t3]dx+mL0[y2t+2θ2t]ξdx+EI2L0y2xxdx+GK2L0θ2xdx(Mg+2mg)L0ydx2H0L0syxdx+H0L0[y2x+2θ2xξ2]dx+AELc[(L0syxξdx)2+(L0sθxξdx)2]+H0L0[syx(y2x+32θ2x)ξ4]dxAELcL0syxξdxL0[y2x+2θ2xξ3]dxAELcL0sθxξdxL0yxθxξ3dx. (1.5)

    Although the energy is conserved, we cannot guarantee that the solution of (1.4) is global in time. The reason is that the functions syx, sθx and yxθx may have indefinite sign, one does not know if the sum of the last three terms in (1.5) is positive or negative. If the sum has the negative sign, then it can happen that it goes to at some time whereas the sum of the remainder terms in (1.5) goes to + even if the energy still remains constant. In that case the solution may blow up in finite time.


    2. Preliminaries

    On the Hilbert space L2(0,L), we denote the standard scalar product and norm by

    (u,v)2:=L0uvdx;u2:=(L0u2dx)1/2,u,vL2(0,L).

    Then the scalar products, the norms of the Hilbert spaces H10(0,L) and H2(0,L) can be defined, respectively, by

    (u,v)H10:=(u,v)2,uH10:=u2,u,vH10(0,L);(u,v)H2:=(u,v)2,uH2:=u2,u,vH2(0,L).

    Recalling that the spaces H2, H3, H4 in Section 1, we denote the duality pairings between H2, H3, H4 and their dual spaces, respectively, by ,H2, ,H3 and ,H4. Given a bounded function W=W(x)>0, we define the weighted scalar products and the corresponding weighted norms on these spaces, respectively, by

    (u,v)H2W:=(Wu,Wv)2,uH2W:=Wu2,u,vH2(0,L);(u,v)H3W:=(Wu,Wv)2,uH3W:=Wu2,u,vH3(0,L);(u,v)H4W:=(Wu,Wv)2,uH4W:=Wu2,u,vH4(0,L).

    We remark that W=Φ,Ψy(x,t),K and F(t) in the sequel.

    Assume that s=s(x)C4(0,L)C3([0,L]) and let s0=maxx[0,L]|s|. For any T>0, we define a bounded set S by

    S:={yC0([0,T];H4(0,L)):supt[0,T]y2H4ΦΛ2},

    where Λ satisfies

    0<Λ<M+6mM+6m1+s20[2H01+s0+GK2][6H0s0L5/2+2AELcL7/2]1. (2.1)

    For any T>0, let

    ¯Y:=C0([0,T];H4(0,L))C1([0,T];H2(0,L)),¯Θ:=C0([0,T];H3(0,L))C1([0,T];H2(0,L)), (2.2)

    then ¯Y×¯Θ is a Banach space when it endowed with the norm

    (y,θ)¯YׯΘ:=(supt[0,T]y(t)2H4Φ+supt[0,T]yt(t)2H2+supt[0,T]θ(t)2H3Ψy(x,t)+supt[0,T]θt(t)2H2)1/2.

    Assume that VL2 is a Hilbert space endowed with the scalar product (,)V. We denote its dual space by V and the duality pairing by ,V. Given a bilinear continuously symmetric form a(,) on V, one can associate with a(,) a linear operator A from V into V defined by

    Au,vV=a(u,v),u,vV,

    which is self-adjoint. We denote by A1 its inverse which is also self-adjoint.

    Let H be another Hilbert space satisfying VH and we restrict to the case where the injection of V in H is compact. In this case A1 can be considered as a self-adjoint compact operator in H. Then we are able to use the elementary spectral theory of self-adjoint compact operators in Hilbert space (see for instance [7]) and we infer that there exists a complete orthonormal family {ek}kN of H such that

    A1ek=μkek,ek=μkAekkN,

    where {μk}kN is decreasing and goes to 0. The family {ek}kN is also orthogonal for a(,) in V.

    Problem 1. Assume that 0<K=K(x) is a bounded function. We consider the following linear problem in H4(0,L)

    {y=λKyxxxx,x(0,L),y(x)=yxx(x)=0,x{0,L}.

    Let A=Kd4dx4. Then from the above argument there exist a basis {ϵk}kN of H4(0,L) (which ia also a basis of H2(0,L)) and a positive sequence {λk} such that

    (ϵk,ϵk)H2=λk(ϵk,ϵk)H4K.

    Problem 2. Given T>0, for every t[0,T], let 0<F(t)=F(x,t) be a bounded function in x. For every t[0,T], we consider the problem

    {θ=μF(t)θxx,x(0,L),θ(x)=0,x{0,L}.

    Assume that A=F(t)d2dx2 for every t[0,T]. Then there exists a sequence of functions {μi(t)} and a complete orthonormal family {ei}iN of H2(0,L) and H10(0,L) such that

    (ei,ei)H10=μi(t)(ei,ei)H2F(t).

    Problem 3. For every t[0,T], we consider the linear problem in H3(0,L)

    {θx=νF(t)θxxx,x(0,L),θ(x)=0,x{0,L}.

    In fact, one can similarly obtain that a basis denoted by {uk}kN of the Hilbert spaces H3(0,L) and H2(0,L) satisfies

    (uk,uk)H2=νi(t)(uk,uk)H3F(t).

    In the sequel, we will use the Gronwall-type inequality, which can be deduced from [4] or [5,Lemma A.5/p.157].

    Proposition 1. Let ψC1(R+) be such that ψ(0)=0, 0ψ(t)C1+C2ψ(t) (with C1,C2>0) for all t>0. Then

    ψ(t)(C1+C2)24t2+(C1+C2)tt>0.

    Now we deal with several linear problems. For any T>0, assume that g=g(x,t)C0([0,T];L2(0,L)) and gtC0([0,T];L2(0,L)), we consider

    {ytt+Kyxxxx=g,x(0,L), t>0,y(x,t)=yxx(x,t)=0,x{0,L}, t>0,y(x,0)=y0,yt(x,0)=y1,x(0,L), (2.3)

    where K>0 is a bounded function as in Problem 1. Then we prove

    Lemma 2.1. Assume that y0H4(0,L) and y1H2(0,L). Then for any T>0 the problem (2.3) admits a unique strong solution

    y¯YC2([0,T];L2(0,L)),

    where ¯Y is defined in (2.2).

    Proof. Assume that {ϵk} is a basis of H4(0,L), and also is a basis of H2(0,L), see Problem 1. Let ϵkH4K=1, then ϵkH2=λk0. For any n1, denote En:=span{ϵ1,ϵ2,...,ϵn} and we put

    (y0)n:=nk=1(y0,ϵk)H4Kϵky0 in H4(0,L),(y1)n:=nk=1(y1,ϵk)H2ϵky1 in H2(0,L). (2.4)

    Then we look for

    yn=nk=1ykn(t)ϵk,

    which solves the following variational problem

    ((yn)tt,φ)H2+(yn,φ)H4K=(g,φxxxx)2,φEn,t>0. (2.5)

    Taking φ=ϵk (k=1,2,...,n) in the equation (2.5), we obtain n equations

    λk(ykn)tt+ykn(t)=(g,(ϵk)xxxx)2. (2.6)

    According to the theory of linear ODE's, one can find a unique solution ykn(t) to (2.6) for all t[0,T] and the initial conditions

    ykn(0)=(y0,ϵk)H4K,(ykn)(0)=(y1,ϵk)H2.

    Hence, there exists a unique solution yn(x,t) of (2.5).

    For any fixed T>0, let φ=(yn)t in the equation (2.5) and integrating on (0,t) with t<T, then

    (yn)t2H2+yn2H4K=(y1)n2H2+(y0)n2H4K+2t0L0(yn)xxxxtgdxds=(y1)n2H2+(y0)n2H4K+2((yn)xxxx(s),g(s))2|t02t0L0(yn)xxxxgtdxds.

    Since gC0([0,T];L2(0,L)) and gtC0([0,T];L2(0,L)), there exists a constant σ(0,1) such that (by the Hölder inequality and the Young inequality)

    2|((yn)xxxx(s),g(s))2|t0t0L0(yn)xxxxgtdxds|C+σyn2H4K+CT(t0yn(s)2H4Kds)1/2.

    Hence,

    (yn)t2H2+(1σ)yn2H4KC+CT(t0yn(s)2H4Kds)1/2. (2.7)

    Then by Proposition 1, (2.7) can be written as

    (yn)t2H2+(1σ)yn2H4KC+CT,t(0,T), (2.8)

    where C and CT are independent of n. Hence, the components ykn do not depend on n, that is,

    yn=yn(x,t)=nk=1yk(t)ϵk.

    Now we prove that {yn} admits a strongly convergent subsequence in a suitable sense. Let m>n1 and define

    ym,n:=ym(x,t)yn(x,t)=mk=n+1yk(t)ϵk.

    Hence,

    ym,n(x,0)=(y0)m(y0)n,(ym,n)t(x,0)=(y1)m(y1)n.

    Subtracting the two equations (2.5) (with n and m) and taking φ=(ym,n)t, we obtain by using the orthogonality of the {ϵk} and by integrating over (0,t) with t<T

    (ym,n)t2H2+ym,n2H4K=(y1)m(y1)n2H2+(y0)m(y0)n2H4K0, by(2.4).

    Therefore, {yn} is a Cauchy sequence in ¯Y. By completeness of these spaces we conclude that there exists y¯Y such that

    ynyin ¯Y as n.

    From the equation in (2.3), we also have the additional regularity yC2([0,T];L2(0,L)).

    Finally, arguing by contradiction and assuming that there are two solutions, we subtract the two linear equations for the two solutions and we obtain a homogeneous linear problem; then one can easily show that the two solutions are identical.

    For any T>0, let f=f(x,t)C0([0,T];L2(0,L)) and ftC0([0,T];L2(0,L)). We consider another linear problem

    {θttF(t)θxx=f,x(0,L), t(0,T),θ(x,t)=0,x{0,L}, t(0,T),θ(x,0)=θ0,θt(x,0)=θ1,x(0,L), (2.9)

    where F(t)=F(x,t)>0 satisfies F(t) and Ft(t) are bounded. Then we prove

    Lemma 2.2. Assume that θ0H2(0,L) and θ1H10(0,L). Then for a small T>0 the problem (2.9) admits a unique solution

    θC0([0,T];H2(0,L))C1([0,T];H10(0,L))C2([0,T];L2(0,L)).

    Proof. Assume that {ei} is a basis of the spaces H2(0,L) and H10(0,L) with eiH2F(t)=Fi(t)>0 for every t[0,T]. Then eiH10=μi(t)Fi(t)=αi with αi0 as i+.

    For any n1, denote a subset En:=span{e1,e2,...,en} and let

    (θ0)n=ni=1(θ0,ei)H2F(0)eiθ0 in H2(0,L)(θ1)n=ni=1(θ0,ei)H10eiθ1 in H10(0,L). (2.10)

    Then we look for θn in the form of

    θn=ni=1θin(t)ei,

    which solves the following approximating problem

    (θn)ttF(t)(θn)xx=f,x(0,L),t>0. (2.11)

    Testing the equation (2.11) with φ=(ei)xx (i=1,2,...,n) and integrating over (0,L), we obtain n equations

    α2i(θin(t))tt+F2i(t)θin(t)=(f,(ei)xx)2. (2.12)

    According to the method of successive approximations (see Coddington-Levinson [6]), one can find a unique solution θin(t) to (2.12) in [0,T] with T>0 small, and hence there exists a unique solution θn(x,t) of (2.11) with the initial conditions

    θin(0)=(θ0,ei)H2F(0),(θin)(0)=(θ1,ei)H10.

    Then for this fixed T>0, testing the equation (2.11) with φ=(θn)xxt and integrating on (0,L)×(0,t) with t<T, we have

    (θn)t2H10+θn2H2F(t)=(θ1)n2H10+(θ0)n2H2F(0)2(f(t),(θn)xx(t))2+2(f(0),(θn)xx(0))2+2t0((θn)xx,ft)2dτ+t0((θn)2xx,Ft)2dτ.

    Since F,Ft are bounded and f,ftC0([0,T];L2(0,L)), by the Hölder inequality and the Young inequality, there exists δ(0,1) such that

    (θn)t2H10+δθn2H2F(t)C+Ct0θn(τ)2H2F(τ)dτC+CTsupt[0,T]θn2H2F(t).

    For a smaller T>0 if necessary, it follows that

    supt[0,T](θn)t2H10+supt[0,T]θn2H2F(t)CT,t(0,T),

    where CT is independent of n and we can write

    θn=ni=1θi(t)ei.

    Let m>n1 and define

    θm,n:=θm(x,t)θn(x,t)=mi=n+1θi(t)ei,

    and then

    θm,n(x,0)=(θ0)m(θ0)n,(θm,n)t(x,0)=(θ1)m(θ1)n.

    Subtracting the two equations (2.11) (with n and m) and testing the difference equation with φ=(θm,n)xxt, we obtain by integrating over (0,L)×(0,t) with t<T

    (θm,n)t2H10+θm,n2H2F(t)=Cm,n+t0(Ft,(θm,n)2xx)2dsCm,n+ct0θm,n(s)2H2F(s)dsCm,n+cTsupt[0,T]θm,n2H2F(t),

    where Cm,n=(θ1)m(θ1)n2H1+(θ0)m(θ0)n2H2F(0). Hence, for a smaller T>0 if necessary, we have

    supt[0,T](θm,n)t2H10+supt[0,T]θm,n2H2F(t)(1cT)Cm,n0, by(2.10).

    Therefore, {θn} is a Cauchy sequence in C0([0,T];H2(0,L))C1([0,T];H10(0,L)). By completeness of these spaces we conclude that there exists θC0([0,T];H2(0,L))C1([0,T];H10(0,L)) such that

    θnθin C0([0,T];H2(0,L))C1([0,T];H10(0,L)) as n.

    From the equation in (2.9), one can prove that θttC0([0,T];L2(0,L)) and hence θC2([0,T];L2(0,L)). By contradiction one can get the uniqueness result and we finish the proof.

    Furthermore, if F(x,t)C2([0,L])×C1([0,T]), fC0([0,T];H2(0,L)), ftC0([0,T];L2(0,L)), then we have

    Lemma 2.3. Assume that θ0H3(0,L) and θ1H2(0,L). Then for a small T>0 the unique solution θ of the problem (2.9) satisfies

    θ¯ΘC2([0,T];H10(0,L)),

    where ¯Θ is defined in (2.2).

    Proof. Assume that {uk} is the basis of H2(0,L) and H3(0,L). Then the solutions θn of the problem (2.11) can be written in the form of

    θn=nk=1ϑkn(t)uk.

    Let

    (θ0)n:=nk=1(θ0,uk)H3F(0)uk,(θ1)n:=nk=1(θ1,uk)H2uk,

    so that

    (θ0)nθ0 in H3(0,L),(θ1)nθ1 in H2(0,L) as n. (2.13)

    Differentiating (2.11) with respect to x, we have

    (θn)xttF(t)(θn)xxx=Fx(t)(θn)xx+fx,x(0,L),t>0. (2.14)

    For the fixed T>0 in Lemma 2.2, testing (2.14) with (θn)xxxt and integrating over (0,L)×(0,t) with t<T, we are led to

    (θn)t2H2+θn2H3F(t)=(θ1)n2H2+(θ0)n2H3F(0)2t0(Fx(s)(θn)xx,(θn)xxxt)2ds+t0(Ft,(θn)2xxx)2ds2t0(fx,(θn)xxxt)2ds.

    Now we estimate the three nonlocal terms. Since F(x,t)C2([0,L])×C1([0,T]), we have

    t0(Ft,(θn)2xxx)2dsCt0θn(s)2H3F(s)ds

    and

    t0(Fx(s)(θn)xx,(θn)xxxt)2ds=t0(Fxx(s)(θn)xx,(θn)xxt)2dst0(Fx(s)(θn)xxx,(θn)xxt)2dsCt0θn(s)H2(θn)t(s)H2ds+Ct0θn(s)H3F(s)(θn)t(s)H2dsCt0θn(s)2H3F(s)ds+Ct0(θn)t(s)2H2ds.

    The condition fC0([0,T];H2(0,L)) allows us to estimate

    t0(fx,(θn)xxxt)2ds=t0(fxx,(θn)xxt)2dsC(t0fxx(s)22ds)1/2(t0(θn)t(s)2H2ds)1/2C+t0(θn)t(s)2H2ds.

    Hence, we have

    (θn)t2H2+θn2H3F(t)C+Ct0θn(s)2H3F(s)ds+Ct0(θn)t(s)2H2ds.

    For a smaller T>0 if necessary, one gets

    supt[0,T](θn)t2H2+supt[0,T]θn2H3F(t)C,C is independent of n.

    Then there exists a unique θ¯Θ (by following the procedure in the proof of Lemma 2.2) such that

    θnθin ¯Θ as n.

    From the equation in (2.9), one can prove that θttC0([0,T];H10(0,L)) and then θC2([0,T];H10(0,L)). The proof is finished.

    With the results above and according to Ball [3], where the author considered some nonlinear (nonlocal) hyperbolic problems, we are able to prove

    Theorem 2.4. Assume that sC4(0,L)C3([0,L]) and q(x,t)S. Let y0H4, θ0H3 and y1,θ1H2. If T>0 is small enough, then for any (z,η)¯Y×¯Θ there exists a unique solution (y,θ)¯Y×¯Θ of the problem

    {ytt+Φyxxxx=h1(z,η),x(0,L),t>0,θttΨq(x,t)θxx=h2(z,η),x(0,L),t>0,y(0,t)=y(L,t)=yxx(0,t)=yxx(L,t)=θ(0,t)=θ(L,t)=0,t>0,y(x,0)=y0,yt(x,0)=y1,θ(x,0)=θ0,θt(x,0)=θ1,x(0,L). (2.15)

    Proof. The condition q(x,t)S leads to the function Ψq(x,t)>0 and

    Ψq(x,t)C2([0,L])×C0([0,T]).

    Moreover, (z,η)¯YׯΘ, then we have

    h1(z,η)C0([0,T];L2(0,L)),h2(z,η)C0([0,T];H2(0,L))

    and

    (h1(z,η))t,(h2(z,η))tC0([0,T];L2(0,L)).

    Hence, the statements of Theorem 2.4 follows by Lemmas 2.1, 2.2 and 2.3.


    3. Proof of Theorem 1.1

    This section is devoted to the proof of Theorem 1.1. First we define a subset of ¯Y×¯Θ (see (2.2)) by

    Q:={(y,θ)¯YׯΘ:(y,θ)2¯YׯΘγΛ2 and (1.3) holds},

    where Λ satisfies (2.1) and γ(0,1). For any (z,η)Q¯YׯΘ, the unique solution (y,θ) of (2.15) allows us to define a map Γ:Q¯Y×¯Θ by

    (y,θ)=Γ(z,η),(z,η)Q.

    Let y0H4, θ0H3, y1,θ1H2 and h(0)=h1(y(x,0),θ(x,0))L2(0,L) satisfy

    y02H4Φ+y12H2+θ02H3Ψq(x,0)+θ12H2+2h(0)2y0H4Φ+16h(0)22ρΛ2, (3.1)

    where 0<ρ<γ/T0 with T0>2. Then we prove

    Lemma 3.1. Assume that (3.1) holds. Then if T>0 is sufficiently small, Γ(Q)Q.

    Proof. Given any (z,η)Q, the solution (y,θ) satisfies

    y2H4Φ+yt2H2+θ2H3Ψq(x,t)+θt2H2=y02H4Φ+y12H2+θ02H3Ψq(x,0)+θ12H2+t0((Ψq(x,t))t,θ2xxx)2ds2t0((Ψq(x,t))xθxx,θxxxt)2ds2t0((h2(z,η))x,θxxxt)2ds+2t0L0(y)xxxxth1(z,η)dxds.

    Now we estimate the nonlocal terms. Since Ψq(x,t)C2([0,L])×C1([0,T]), we have

    t0((Ψq(x,t))t,θ2xxx)2dsct0θ2H3Ψq(x,s)dscTsupt[0,T]θ(t)2H3Ψq(x,t)c1T(y(t),θ(t))2¯YׯΘ

    and

    2t0((Ψq(x,t))xθxx,θxxxt)2ds=2t0((Ψq(x,t))xxθxx,θxxt)2ds2t0((Ψq(x,t))xθxxx,θxxt)2dscT(supt[0,T]θ(t)2H3Ψq(x,t)+supt[0,T]θt(t)2H2)c2T(y(t),θ(t))2¯YׯΘ.

    The condition h2(z,η)C0([0,T];H2(0,L)) allows us to obtain by the Young inequality

    2t0((h2(z,η))x,θxxxt)2ds=2t0((h2(z,η))xx,θxxt)2dsCT+14supt[0,T]θt(t)2H2CT+14(y(t),θ(t))2¯YׯΘ;

    Since (z,η)Q implies that h1(z,η)C1([0,T]) for every x(0,L), it follows that

    |h1(z(x,t),η(x,t))|=|h(t)||h(0)|+ChT, for every t[0,T].

    Hence,

    2t0L0(y)xxxxth1(z,η)dxds=2(yxxxx(t),h1(z,η))22(yxxxx(0),h(0))22t0L0(h1(z,η))tyxxxxdxds2h1(z,η)2y(t)H4Φ+2h(0)2y0H4Φ+CT+18supt[0,T]y(t)2H4Φ14supt[0,T]y(t)2H4Φ+16(h(0)22+(ChT)2L)+CT+2h(0)2y0H4Φ14(y(t),θ(t))2¯YׯΘ+16h(0)22+CT2+CT+2h(0)2y0H4Φ.

    Then we deduce that

    (12(c1+c2)T)(y(t),θ(t))2¯YׯΘρΛ2+CT2+CT.

    Let TT022(c1+c2)T0 with T0>2, then

    (y(t),θ(t))2¯YׯΘρT0Λ2+T0(CT2+CT).

    If T>0 is smaller enough, then

    (y(t),θ(t))2¯YׯΘγΛ2,

    which proves that Γ(Q)Q.

    Next, we show that the map Γ is contractive for sufficiently small T.

    Lemma 3.2. Assume that (3.1) holds. If T>0 is sufficiently small, then Γ is a contractive map.

    Proof. Let (za,ηa), (zb,ηb) be two different elements in Q and let (ya,θa)=Γ(za,ηa), (yb,θb)=Γ(zb,ηb) and we denote

    z=zazb,η=ηaηb,y=yayb,θ=θaθb,

    which yields that z(x,0)=η(x,0)=y(x,0)=θ(x,0)=0, and hence,

    h1(za(x,0),ηa(x,0))h1(zb(x,0),ηb(x,0))=0. (3.2)

    Differentiating the second equation in (2.15) with respect to x both side and then subtracting the two problems (2.15) for (ya,θa) and (yb,θb), we get a system

    {ytt+Φyxxxx=h1(za,ηa)h1(zb,ηb),x(0,L),t>0,θxttΨq(x,t)θxxx=(Ψq(x,t))xθxx+(h2(za,ηa))x(h2(zb,ηb))x,x(0,L),t>0,y(0,t)=y(L,t)=yxx(0,t)=yxx(L,t)=θ(0,t)=θ(L,t)=0,t>0,y(x,0)=0,yt(x,0)=0,θ(x,0)=0,θt(x,0)=0,x(0,L). (3.3)

    Multiplying (3.3) by (yxxxxt,θxxxt) and integrating over (0,L)×(0,t) with t<T, one has

    y2H4Φ+yt2H2+θ2H3Ψq(x,t)+θt2H2=2t0L0(h1(za,ηa)h1(zb,ηb))yxxxxtdxds+t0L0[(Ψq(x,t))tθ2xxx2(Ψq(x,t))xθxxθxxxt]dxds2t0L0[h2(za,ηa)h2(zb,ηb)]xθxxxtdxds:=I+IIIII.

    Since h1(z,η)C1([0,T];L2(0,L)) and h1(z(,t),η(,t))C1([0,T]), we have

    h1(za,ηa)h1(zb,ηb)2CTsupt[0,T][h1(za,ηa)h1(zb,ηb)]t2CT(z(t),η(t))¯YׯΘ.

    Hence, from (3.2) we obtain

    I=2(yxxxx(t),[h1(za,ηa)h1(zb,ηb)]|t0)22t0L0(h1(za,ηa)h1(zb,ηb))tyxxxxdxdsCTsupt[0,T][h1(za,ηa)h1(zb,ηb)]t2(y(t)H4Φ+CT(y(t)H4Φ(z(t),η(t))Y×ΘCT(y(t),θ(t))¯YׯΘ(z(t),η(t))¯YׯΘ.

    Since Ψq(x,t)C2([0,L])×C1([0,T]), we have

    II=t0L0(Ψq(x,t))tθ2xxxdxds2t0L0(Ψq(x,t))xxθxxθxxtdxds2t0L0(Ψq(x,t))xθxxxθxxtdxdscT(supt[0,T]θ(t)2H3Ψq(x,t)+supt[0,T]θt(t)2H2)cT(y(t),θ(t))2¯YׯΘ.

    Since h2(z,η)C0([0,T];H2(0,L)),(h2(z,η))tC0([0,T]×L2(0,L)), it yields that

    III=2t0L0(h2(za,ηa)h2(zb,ηb))xxθxxtdxdsCT(y(t),θ(t))¯YׯΘ(z(t),η(t))¯YׯΘ.

    Therefore,

    (1cT)(y(t),θ(t))2¯YׯΘCT(y(t),θ(t))¯YׯΘ(z(t),η(t))¯YׯΘ.

    Then if T>0 is small enough, we have

    (y(t),θ(t))¯YׯΘ<α(z(t),η(t))¯YׯΘ, with 0<α<1

    and we finish the proof.

    Then by the Contraction Mapping principle, there exists a unique local solution to (2.15) with h1(y,θ) and h2(y,θ) in place of h1(z,η) and h2(z,η).

    Now we prove the existence result of (1.4). For any given q0S, there exists a unique solution (y1,θ1)Q of the problem (1.4) with Ψy(x,t)=Ψq0(x,t). Moreover, y1S. Let Ψy(x,t)=Ψy1(x,t) in (1.4), then the problem (1.4) with Ψy(x,t)=Ψy1(x,t) has a unique solution which is denoted by (y2,θ2). Repeating this process, we obtain that the problem (1.4) with Ψy(x,t)=Ψyn1(x,t) has a unique solution (yn,θn) (n1).

    Since (yn,θn)Q, we have that (yn,θn) is bounded in C2([0,T];L2(0,L))×C2([0,T];H10(0,L)) from the system (1.4). Hence, there exists (y,θ)Q(C2([0,T];L2(0,L))×C2([0,T];H10(0,L))) such that for every t[0,T]

    (yn,θn)(y,θ), weakly in H4(0,L)×H3(0,L),
    ((yn)tt,(θn)tt)(ytt,θtt), weakly in L2(0,L)×H10(0,L),

    which implies that for every t[0,T]

    (yn)xxxxyxxxx, weakly in L2(0,L)

    and for every t[0,T]

    h1(yn,θn)h1(y,θ),h2(yn,θn)h2(y,θ),Ψyn(x,t)Ψy(x,t), strongly in L2(0,L).

    Moreover, for every t[0,T], we have (θn)xxθxx strongly in L2(0,L), then it yields that

    Ψyn1(x,t)(θn)xxΨy(x,t)θxx2=(Ψyn1(x,t)Ψy(x,t))(θn)xx2+Ψy(x,t)((θn)xxθxx)2CΨyn1(x,t)Ψy(x,t)2+C(θn)xxθxx20.

    Then (y,θ)¯Y×¯Θ satisfies the two equations of (1.4) in L2(0,L) for every t[0,T] and we claim that (y,θ) is a solution of the problem (1.4) and we finish the proof.


    Acknowledgments

    The author is grateful to his colleagues for several remarks which lead to an improvement of the present paper. This work is supported by Gaocengci Rencai Zizhu Xiangmu of Hebei Province: CL201701 and partially supported by National Natural Science Foundation of China: 11701162.


    Conflict of interest

    The author declares no conflict of interest.


    [1] G. Arioli, F. Gazzola, On a nonlinear nonlocal hyperbolic system modeling suspension bridges, Milan J. Math., 83 (2015), 211-236.
    [2] G. Arioli, F. Gazzola, Torsional instability in suspension bridges: the Tacoma Narrows Bridge case, Commun. Nonlinear Sci., 42 (2017), 342-357.
    [3] J. M. Ball, Initial-boundary value problems for an extensible beam, J. Math. Anal. Appl., 42 (1973), 61-90.
    [4] C. Baiocchi, Soluzioni ordinarie e generalizzate del problema di Cauchy per equazioni differenziali astratte non lineari del secondo ordine in spazi di Hilbert, Ricerche Mat., 16 (1967), 27-95.
    [5] H. Brezis, Operateurs maximaux monotones, North-Holland, Amsterdam, 1973.
    [6] E. A. Coddington, N. Levinson, Theory of ordinary differential equations, McGraw-Hill, New York, 1955.
    [7] R. Courant, D. Hilbert, Methods of mathematical physics, Intersciences Publishers, New York, 1953.
    [8] A. Falocchi, Structural instability in a simplified nonlinear model for suspension bridges, AIMETA 2017-Proceedings of the 23rd Conference of the Italian Association of Theoretical and Applied Mechanics, 4 (2017), 746-759.
    [9] A. Falocchi, Torsional instability and sensitivity analysis in a suspension bridge model related to the Melan equation, Commun. Nonlinear Sci., 67 (2019), 60-75.
    [10] A. Ferrero, F. Gazzola, A partially hinged rectangular plate as a model for suspension bridges, Discrete Cont. Dyn-A, 35 (2015), 5879-5908.
    [11] R. Temam, Infinite-dimensional dynamical systems in mechanics and physics, Applied Mathematical Sciences, 68, Springer, 1997.
    [12] Y. Wang, Finite time blow-up and global solutions for fourth order damped wave equations, J. Math. Anal. Appl., 418 (2014), 713-733.
  • Reader Comments
  • © 2018 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4029) PDF downloads(600) Cited by(0)

Article outline

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog