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
[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 |
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)x−EALc[∫L0y2x+ℓ2θ2xξ3]s″ξ3 −2EALc[∫L0s″yξ3](s′ξ−yxξ3)x+2EAℓ2Lc[∫L0s″θξ3](θxξ3)x,x∈(0,L),t>0,(M3+2mξ)θtt=GKℓ2θxx+2H0(θxξ2+3s′yxθxξ4)x−2EALc[∫L0yxθxξ3]s″ξ3 −2EALc[∫L0s″θξ3](s′ξ−yxξ3)x+2EALc[∫L0s″yξ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.
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(GKℓ2+2H0ξ2+6H0s′yxξ4+2AEξ3Lc∫L0s″yξ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)x−EA(M+2mξ)Lc[∫L0y2x+ℓ2θ2xξ3]s″ξ3−2EA(M+2mξ)Lc[∫L0s″yξ3](s′ξ−yxξ3)x+2EAℓ2(M+2mξ)Lc[∫L0s″θξ3](θxξ3)x,h2(y,θ)=(M3+2mξ)−1[(2H0ξ2+6H0s′yxξ4+2AEξ3Lc[∫L0s″yξ3])xθx−2EALc[∫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):={u∈H2(0,L);u(x)=0 on {0,L}},H3∗=H3∗(0,L):={u∈H3(0,L);u(x)=u″(x)=0 on {0,L}},H4∗=H4∗(0,L):={u∈H4(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 s∈C4(0,L)∩C3([0,L]). Let y0∈H4∗, θ0∈H3∗ and y1,θ1∈H2∗. 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 θxx∈C0([0,T]×[0,L]), and then θtt∈C0([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 t≥0 as follows
E(t)=M2∫L0[y2t+ℓ2θ2t3]dx+m∫L0[y2t+ℓ2θ2t]ξdx+EI2∫L0y2xxdx+GK2∫L0θ2xdx−(Mg+2mg)∫L0ydx−2H0∫L0s′yxdx+H0∫L0[y2x+ℓ2θ2xξ2]dx+AELc[(∫L0s′yxξdx)2+(∫L0ℓs′θxξdx)2]+H0∫L0[s′yx(y2x+3ℓ2θ2x)ξ4]dx−AELc∫L0s′yxξdx∫L0[y2x+ℓ2θ2xξ3]dx−AELc∫L0ℓs′θxξdx∫L0ℓyxθ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 s′yx, 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.
On the Hilbert space L2(0,L), we denote the standard scalar product and norm by
(u,v)2:=∫L0uvdx;‖u‖2:=(∫L0u2dx)1/2,∀u,v∈L2(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,‖u‖H10:=‖u′‖2,∀u,v∈H10(0,L);(u,v)H2∗:=(u″,v″)2,‖u‖H2∗:=‖u″‖2,∀u,v∈H2∗(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,‖u‖H2W:=‖√Wu″‖2,∀u,v∈H2∗(0,L);(u,v)H3W:=(√Wu‴,√Wv‴)2,‖u‖H3W:=‖√Wu‴‖2,∀u,v∈H3∗(0,L);(u,v)H4W:=(√Wu⁗,√Wv⁗)2,‖u‖H4W:=‖√Wu⁗‖2,∀u,v∈H4∗(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 ‖s′‖0=maxx∈[0,L]|s′|. For any T>0, we define a bounded set S by
S:={y∈C0([0,T];H4∗(0,L)):supt∈[0,T]‖y‖2H4Φ≤Λ2}, |
where Λ satisfies
0<Λ<M+6mM+6m√1+‖s′‖20[2H01+‖s′‖0+GKℓ2][6H0‖s′‖0L5/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 V⊂L2 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,v⟩V=a(u,v),∀u,v∈V, |
which is self-adjoint. We denote by A−1 its inverse which is also self-adjoint.
Let H be another Hilbert space satisfying V⊂H and we restrict to the case where the injection of V in H is compact. In this case A−1 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}k∈N of H such that
A−1ek=μkek,ek=μkAek∀k∈N, |
where {μk}k∈N is decreasing and goes to 0. The family {ek}k∈N 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}k∈N 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}i∈N 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}k∈N 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)t∀t>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 gt∈C0([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 y0∈H4∗(0,L) and y1∈H2∗(0,L). Then for any T>0 the problem (2.3) admits a unique strong solution
y∈¯Y∩C2([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 ‖ϵk‖H4K=1, then ‖ϵk‖H2∗=√λk↓0. For any n≥1, denote En:=span{ϵ1,ϵ2,...,ϵn} and we put
(y0)n:=n∑k=1(y0,ϵk)H4Kϵk→y0 in H4∗(0,L),(y1)n:=n∑k=1(y1,ϵk)H2∗ϵk→y1 in H2∗(0,L). | (2.4) |
Then we look for
yn=n∑k=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)t‖2H2∗+‖yn‖2H4K=‖(y1)n‖2H2∗+‖(y0)n‖2H4K+2∫t0∫L0(yn)xxxxtgdxds=‖(y1)n‖2H2∗+‖(y0)n‖2H4K+2((yn)xxxx(s),g(s))2|t0−2∫t0∫L0(yn)xxxxgtdxds. |
Since g∈C0([0,T];L2(0,L)) and gt∈C0([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|t0−∫t0∫L0(yn)xxxxgtdxds|≤C+σ‖yn‖2H4K+CT(∫t0‖yn(s)‖2H4Kds)1/2. |
Hence,
‖(yn)t‖2H2∗+(1−σ)‖yn‖2H4K≤C+CT(∫t0‖yn(s)‖2H4Kds)1/2. | (2.7) |
Then by Proposition 1, (2.7) can be written as
‖(yn)t‖2H2∗+(1−σ)‖yn‖2H4K≤C+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)=n∑k=1yk(t)ϵk. |
Now we prove that {yn} admits a strongly convergent subsequence in a suitable sense. Let m>n≥1 and define
ym,n:=ym(x,t)−yn(x,t)=m∑k=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)t‖2H2∗+‖ym,n‖2H4K=‖(y1)m−(y1)n‖2H2∗+‖(y0)m−(y0)n‖2H4K→0, 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
yn→yin ¯Y as n→∞. |
From the equation in (2.3), we also have the additional regularity y∈C2([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 ft∈C0([0,T];L2(0,L)). We consider another linear problem
{θtt−F(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 θ0∈H2∗(0,L) and θ1∈H10(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 ‖ei‖H2F(t)=Fi(t)>0 for every t∈[0,T]. Then ‖ei‖H10=√μi(t)Fi(t)=αi with αi↓0 as i→+∞.
For any n≥1, denote a subset En:=span{e1,e2,...,en} and let
(θ0)n=n∑i=1(θ0,ei)H2F(0)ei→θ0 in H2∗(0,L)(θ1)n=n∑i=1(θ0,ei)H10ei→θ1 in H10(0,L). | (2.10) |
Then we look for θn in the form of
θn=n∑i=1θin(t)ei, |
which solves the following approximating problem
(θn)tt−F(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)t‖2H10+‖θn‖2H2F(t)=‖(θ1)n‖2H10+‖(θ0)n‖2H2F(0)−2(f(t),(θn)xx(t))2+2(f(0),(θn)xx(0))2+2∫t0((θn)xx,ft)2dτ+∫t0((θn)2xx,Ft)2dτ. |
Since F,Ft are bounded and f,ft∈C0([0,T];L2(0,L)), by the Hölder inequality and the Young inequality, there exists δ∈(0,1) such that
‖(θn)t‖2H10+δ‖θn‖2H2F(t)≤C+C∫t0‖θn(τ)‖2H2F(τ)dτ≤C+CTsupt∈[0,T]‖θn‖2H2F(t). |
For a smaller T>0 if necessary, it follows that
supt∈[0,T]‖(θn)t‖2H10+supt∈[0,T]‖θn‖2H2F(t)≤CT,∀t∈(0,T), |
where CT is independent of n and we can write
θn=n∑i=1θi(t)ei. |
Let m>n≥1 and define
θm,n:=θm(x,t)−θn(x,t)=m∑i=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)t‖2H10+‖θm,n‖2H2F(t)=Cm,n+∫t0(Ft,(θm,n)2xx)2ds≤Cm,n+c∫t0‖θm,n(s)‖2H2F(s)ds≤Cm,n+cTsupt∈[0,T]‖θm,n‖2H2F(t), |
where Cm,n=‖(θ1)m−(θ1)n‖2H1+‖(θ0)m−(θ0)n‖2H2F(0). Hence, for a smaller T>0 if necessary, we have
supt∈[0,T]‖(θm,n)t‖2H10+supt∈[0,T]‖θm,n‖2H2F(t)≤(1−cT)Cm,n→0, 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 θtt∈C0([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]), f∈C0([0,T];H2(0,L)), ft∈C0([0,T];L2(0,L)), then we have
Lemma 2.3. Assume that θ0∈H3∗(0,L) and θ1∈H2∗(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=n∑k=1ϑkn(t)uk. |
Let
(θ0)n:=n∑k=1(θ0,uk)H3F(0)uk,(θ1)n:=n∑k=1(θ1,uk)H2∗uk, |
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)xtt−F(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)t‖2H2∗+‖θn‖2H3F(t)=‖(θ1)n‖2H2∗+‖(θ0)n‖2H3F(0)−2∫t0(Fx(s)(θn)xx,(θn)xxxt)2ds+∫t0(Ft,(θn)2xxx)2ds−2∫t0(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)2ds≤C∫t0‖θn(s)‖2H3F(s)ds |
and
∫t0(Fx(s)(θn)xx,(θn)xxxt)2ds=−∫t0(Fxx(s)(θn)xx,(θn)xxt)2ds−∫t0(Fx(s)(θn)xxx,(θn)xxt)2ds≤C∫t0‖θn(s)‖H2∗‖(θn)t(s)‖H2∗ds+C∫t0‖θn(s)‖H3F(s)‖(θn)t(s)‖H2∗ds≤C∫t0‖θn(s)‖2H3F(s)ds+C∫t0‖(θn)t(s)‖2H2∗ds. |
The condition f∈C0([0,T];H2(0,L)) allows us to estimate
∫t0(fx,(θn)xxxt)2ds=−∫t0(fxx,(θn)xxt)2ds≤C(∫t0‖fxx(s)‖22ds)1/2(∫t0‖(θn)t(s)‖2H2∗ds)1/2≤C+∫t0‖(θn)t(s)‖2H2∗ds. |
Hence, we have
‖(θn)t‖2H2∗+‖θn‖2H3F(t)≤C+C∫t0‖θn(s)‖2H3F(s)ds+C∫t0‖(θn)t(s)‖2H2∗ds. |
For a smaller T>0 if necessary, one gets
supt∈[0,T]‖(θn)t‖2H2∗+supt∈[0,T]‖θn‖2H3F(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 θtt∈C0([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 s∈C4(0,L)∩C3([0,L]) and q(x,t)∈S. Let y0∈H4∗, θ0∈H3∗ and y1,θ1∈H2∗. 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,η))t∈C0([0,T];L2(0,L)). |
Hence, the statements of Theorem 2.4 follows by Lemmas 2.1, 2.2 and 2.3.
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 y0∈H4∗, θ0∈H3∗, y1,θ1∈H2∗ and h(0)=h1(y(x,0),θ(x,0))∈L2(0,L) satisfy
‖y0‖2H4Φ+‖y1‖2H2∗+‖θ0‖2H3Ψq(x,0)+‖θ1‖2H2∗+2‖h(0)‖2‖y0‖H4Φ+16‖h(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
‖y‖2H4Φ+‖yt‖2H2∗+‖θ‖2H3Ψq(x,t)+‖θt‖2H2∗=‖y0‖2H4Φ+‖y1‖2H2∗+‖θ0‖2H3Ψq(x,0)+‖θ1‖2H2∗+∫t0((Ψq(x,t))t,θ2xxx)2ds−2∫t0((Ψq(x,t))xθxx,θxxxt)2ds−2∫t0((h2(z,η))x,θxxxt)2ds+2∫t0∫L0(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)2ds≤c∫t0‖θ‖2H3Ψq(x,s)ds≤cTsupt∈[0,T]‖θ(t)‖2H3Ψq(x,t)≤c1T‖(y(t),θ(t))‖2¯Y×¯Θ |
and
2∫t0((Ψq(x,t))xθxx,θxxxt)2ds=−2∫t0((Ψq(x,t))xxθxx,θxxt)2ds−2∫t0((Ψq(x,t))xθxxx,θxxt)2ds≤cT(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
2∫t0((h2(z,η))x,θxxxt)2ds=−2∫t0((h2(z,η))xx,θxxt)2ds≤CT+14supt∈[0,T]‖θt(t)‖2H2∗≤CT+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,
2∫t0∫L0(y)xxxxth1(z,η)dxds=2(yxxxx(t),h1(z,η))2−2(yxxxx(0),h(0))2−2∫t0∫L0(h1(z,η))tyxxxxdxds≤2‖h1(z,η)‖2‖y(t)‖H4Φ+2‖h(0)‖2‖y0‖H4Φ+CT+18supt∈[0,T]‖y(t)‖2H4Φ≤14supt∈[0,T]‖y(t)‖2H4Φ+16(‖h(0)‖22+(ChT)2L)+CT+2‖h(0)‖2‖y0‖H4Φ≤14‖(y(t),θ(t))‖2¯YׯΘ+16‖h(0)‖22+CT2+CT+2‖h(0)‖2‖y0‖H4Φ. |
Then we deduce that
(12−(c1+c2)T)‖(y(t),θ(t))‖2¯YׯΘ≤ρΛ2+CT2+CT. |
Let T≤T0−22(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=za−zb,η=ηa−ηb,y=ya−yb,θ=θ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
‖y‖2H4Φ+‖yt‖2H2∗+‖θ‖2H3Ψq(x,t)+‖θt‖2H2∗=2∫t0∫L0(h1(za,ηa)−h1(zb,ηb))yxxxxtdxds+∫t0∫L0[(Ψq(x,t))tθ2xxx−2(Ψq(x,t))xθxxθxxxt]dxds−2∫t0∫L0[h2(za,ηa)−h2(zb,ηb)]xθxxxtdxds:=I+II−III. |
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)‖2≤CTsupt∈[0,T]‖[h1(za,ηa)−h1(zb,ηb)]t‖2≤CT‖(z(t),η(t))‖¯YׯΘ. |
Hence, from (3.2) we obtain
I=2(yxxxx(t),[h1(za,ηa)−h1(zb,ηb)]|t0)2−2∫t0∫L0(h1(za,ηa)−h1(zb,ηb))tyxxxxdxds≤CTsupt∈[0,T]‖[h1(za,ηa)−h1(zb,ηb)]t‖2‖(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=∫t0∫L0(Ψq(x,t))tθ2xxxdxds−2∫t0∫L0(Ψq(x,t))xxθxxθxxtdxds−2∫t0∫L0(Ψq(x,t))xθxxxθxxtdxds≤cT(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,η))t∈C0([0,T]×L2(0,L)), it yields that
III=−2∫t0∫L0(h2(za,ηa)−h2(zb,ηb))xxθxxtdxds≤CT‖(y(t),θ(t))‖¯YׯΘ‖(z(t),η(t))‖¯YׯΘ. |
Therefore,
(1−cT)‖(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 q0∈S, there exists a unique solution (y1,θ1)⊂Q of the problem (1.4) with Ψy(x,t)=Ψq0(x,t). Moreover, y1∈S. 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)=Ψyn−1(x,t) has a unique solution (yn,θn) (n≥1).
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)xxxx⇀yxxxx, 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
‖Ψyn−1(x,t)(θn)xx−Ψy(x,t)θxx‖2=‖(Ψyn−1(x,t)−Ψy(x,t))(θn)xx‖2+‖Ψy(x,t)((θn)xx−θxx)‖2≤C‖Ψyn−1(x,t)−Ψy(x,t)‖2+C‖(θn)xx−θxx‖2→0. |
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.
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.
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. |