Acessibilidade / Reportar erro

Asymptotic Behaviour of a Viscoelastic Transmission Problem with a Tip Load

ABSTRACT

We consider a transmission problem for a string composed by two components: one of them is a viscoelastic material (with viscoelasticity of memory type), and the other is an elastic material (without dissipation effective over this component). Additionally, we consider that in one end is attached a tip load. The main result is that the model is exponentially stable if and only if the memory effect is effective over the string. When there is no memory effect, then there is a lack of exponential stability, but the tip load produces a polynomial rate of decay. That is, the tip load is not strong enough to stabilize exponentially the system, but produces a polynomial rate of decay.

Keywords:
transmission problems; memory effect; lack of exponential stability; tip load; hybrid system

1 INTRODUCTION

We consider the transmission problem for the damped vibrations of a string, whose left end is rigidly attached and in the other end has an attached hollow-tip body that contains granular material. The string is composed by two components: one of them is a viscoelastic material (with viscoelasticity of memory type) and the other is an elastic material (without dissipation effective over this component).

More precisely, let us denote by U the displacement of the string. That is

U ( x ) = u ( x ) , x ] 0 , l 0 [ v ( x ) , x ] l 0 , l [

where l is the length of the string and l 0 is the transmission point. The model that we consider in this paper is written as follows.

ρ 1 u l t - α 1 u x x + 0 t g t - s u x x · , s u x x · , s d s = 0 i n ] 0 , l 0 [ × ] 0 , + [ (1.1)

ρ 2 v t t - α 2 v x x = 0 i n ] l 0 , l [ × ] 0 , + [ . (1.2)

Figure 1:
String with Tip Load.

Here, g:[0,+) be the relaxation function, and α 1, α 2, ρ 1, ρ 2 are positive constants that reflect physical properties of the string. The boundary conditions are given by

u ( 0 , t ) = 0 , v ( l , t ) = w ( t ) , t 0 , (1.3)

and the transmission conditions are given by

u ( l 0 , t ) = v ( l 0 , t ) , α 1 u x ( l 0 , t ) - 0 t g ( t - s ) u x ( l 0 , s ) d s = α 2 v x ( l 0 , t ) , t 0 . (1.4)

We turn to model the motion of the right end with the attached tip body. We assume that the container is rigidly attached to the end x = l, and that the container and its contents have mass m and a center of mass O′ located at distance d from the end of the string. We assume that the damping effect of the internal granular material can be represented by damping coefficient γ 1, whose precise contributions are described below

ρ 3 w t t + γ 1 w t + γ 2 w .

Here, the first term is the contribution of the inertia of the container, and the second term represents the damping that the granular material provides, which is assumed to be proportional to the velocity, and so γ 1 is the damping coefficient. Thus, the force balance at the end x = l is

ρ 3 w t t + γ 1 w t + γ 2 w + α 2 v x ( l , . ) = 0 in ] 0 , + [ , (1.5)

where the parameters γ 1 and γ 2 are non-negative constants. Finally, the initial data are given by

u ( 0 ) = u 0 , u t ( 0 ) = u 1 in ] 0 , l 0 [ , (1.6)

v ( 0 ) = v 0 , v t ( 0 ) = v 1 in ] l 0 , l [ , (1.7)

w ( 0 ) = w 0 , w t ( 0 ) = w 1 . (1.8)

Here, we assume the following hypotheses on the relaxation function g:

g ( t ) 0 , t 0 , and g > 0 almost everywhere in [ 0 , + [ ; (1.9)

k 1 , k 2 > 0 : - k 1 g ( t ) g ' ( t ) - k 2 g ( t ) , t 0 ; (1.10)

0 < α : = α 1 - 0 g ( s ) d s . (1.11)

Concerning models of motion with the attached tip body, Andrews and Shillor 11 K. Andrews & M. Shillor. Vibrations of a Beam With a Damping Tip Body. Mathematical and Computer Modelling, 35 (2002), 1033-1042. establish the existence and uniqueness of the model and showed the exponential energy decay of the solution provided and extra damping term is present. In 1111 N. Zietsman, L Van Rensburg & A. Van der Merwe. A Timoshenko Beam With Tip Body and Boundary Damping. Wave Motion, 39 (2004), 199-211. Zietsman, Rensburg and Merwe consider the effect of boundary damping on a cantilevered Timoshenko beam with a rigid body attached to the free end. The authors establish the efficiency and accuracy of the finite element method for calculating the eigenvalues and eigenmodes. Although no conclusion is showed with regard to the stabilization of the system, the authors showed interesting phenomena concerning the damped vibration spectrum and the associated eigenmodes. See also the work of Feireisl and O’Dowd 77 E. Feireisl & G. O’Dowd. Stabilisation d’un syste`me hybride par un feedback non linéaire, non monotone. Comptes Rendus de l’Académie des Sciences, 326 (1998), 323-327. where is showed, for an hybrid system composed of a cable with masses at both end, the strong stability for a nonlinear and nonmonotone feedback law applied at one end.

The main result of this paper is to show that the system (1.1)-(1.8) is exponentially stable if and only if the memory effect is effective over the viscoelastic part of the material. This means that the dissipative properties given by the tip load is not enough to produce exponencial rate of decay when the memory effect is not effective. Finally, when g = 0, we prove that the system is not exponentially stable but the dissipation given by the tip load produce polynomial stability. The method we use is based on Prüss Theorem to show exponential stability. The proof of the lack of exponential stability is based on the Weyl invariance Theorem and the proof of the polinomial stability is based on the Borichev and Tomilov result.

2 EXISTENCE AND UNIQUENESS OF SOLUTIONS

To use the semigroup approach we need to rewrite the problem as an autonomous system. For this reason we introduce the history problem, obtained by replacing the equation (1.1) by the following history equation

ρ 1 u t t - α 1 u x x + - t g ( t - s ) u x x ( . , s ) d s = 0 in ] 0 , l 0 [ × ] 0 , + [ .

Following the ideas of Dafermos 44 C.M. Dafermos. On Abstract Volterra Equation with Applications to Linear Viscoleasticity. Differential and Integral Equations, 7 (1970), 554-569.), (33 C.M. Dafermos. Asymptotic Stability in Viscoleasticity. Arch. Rat. Mech. Anal., 37 (1970), 297-308. and Fabrizio 66 M. Fabrizio & A. Morro. “Asymptotic Stability in Viscoleasticity”, volume 12 of SIAM Studies in Applied Mathematics. Society for Industrial and Applied Mathematics, Philadelphia (1992)., we introduce the notation

η ( x , t , s ) : = u ( x , t ) - u ( x , t - s ) ,

with s[0,+); whence we consider the system

ρ 1 u t t - α u x x - 0 g s η x x s d s = 0 in ] 0 , l 0 [ × ] 0 , + [ (2.1)

ρ 2 v t t - α 2 v x x = 0 in ] l 0 , l [ × ] 0 , + [ (2.2)

η t + η s - u t = 0 in ] 0 , l 0 [ × ] 0 , + [ × ] 0 , + [ . (2.3)

with u, v and w, satisfying (1.5) and the initial conditions (1.6), (1.7), (1.8) and η verifying

η ( x , 0 , s ) = η 0 ( x , s ) = : u 0 ( x ) - u 0 ( x , - s ) , ( x , s ) ] 0 , l 0 [ × ] 0 , + [ , (2.4)

with boundary conditions are given by

η ( x , t , 0 ) = 0 , ( x , t ) ] 0 , l 0 [ × ] 0 , + [ , (2.5)

η ( 0 , t , s ) = 0 , ( t , s ) ] 0 , + [ × ] 0 , + [ . (2.6)

The transmission conditions now are given by

u ( l 0 , t ) = v ( l 0 , t ) , α u x ( l 0 , t ) + 0 g ( s ) η x ( l 0 , t , s ) d s = α 2 v x ( l 0 , t ) , t 0 . (2.7)

We define the total energy of the system as

E ( t ) = 1 2 0 l 0 ρ 1 | u t | 2 + α | u x | 2 + 0 g ( s ) | η x ( s ) | 2 d s d x + + l 0 l ρ 2 | v t | 2 + α 2 | v x | 2 d x + ρ 3 | w t | 2 + γ 2 | w | 2 .

Let us introduce the following spaces:

m : = H m ( 0 , l 0 ) × H m ( l 0 , l ) , m ; * m : = ( u , v ) m ; u ( 0 ) = 0 , u ( l 0 ) = v ( l 0 ) , m ; L 2 : = L 2 ( 0 , l 0 ) × L 2 ( l 0 , l ) ; H * m ( 0 , l 0 ) : = f H m ( 0 , l 0 ) ; f ( 0 ) = 0 , m ; L g 2 : = φ : + H * 1 ( 0 , l 0 ) ; 0 g ( s ) 0 l 0 | φ x ( s ) | 2 d x d s < .

We recall that Lg2 is a Hilbert space when endowed with the inner product given by

φ , ψ L g 2 = 0 g ( s ) 0 l 0 φ x ( s ) ψ x ( s ) ¯ d x d s .

With this notations, we consider the phase space

H : = ( u , v , U , V , η , w , W ) T * 1 × L 2 × L g 2 × 2 ; v ( l ) = w .

Note that the space is a Hilbert space with the norm

U H 2 = α u x L 2 ( 0 , l 0 ) 2 + α 2 v x L 2 ( 0 , l 0 ) 2 + ρ 1 U L 2 ( 0 , l 0 ) 2 + + ρ 2 V L 2 ( 0 , l 0 ) 2 + η L g 2 2 + γ 2 | w | 2 + ρ 3 | W | 2 .

where U=(u,v,U,V,η,w,W)TH.

Let us introduce the linear unbounded operator 𝒜 in as follows

A U = U V α ρ 1 u x x + 1 ρ 1 0 g ( s ) η x x ( s ) d s α 2 ρ 2 v x x U - η s W - γ 1 ρ 3 W - γ 2 ρ 3 w - α 2 ρ 3 v x ( l ) ,

with domain

D ( A ) = U = ( u , v , U , V , η , w , W ) T H ; α u + 0 g ( s ) η ( s ) d s , v 2 , ( U , V ) * 1 , V ( l ) = W , η | s = 0 = 0 , η s L g 2 , α u x ( l 0 ) + 0 g ( s ) η x ( l 0 , s ) d s = α 2 v x ( l 0 ) .

Using the hypotheses on g, a direct computation yields

Re A U , U = - γ 1 | W | 2 + 1 2 0 l 0 0 g ' ( s ) | η x ( s ) | 2 d s d x 0 ,

which means that 𝒜 is a dissipative operator. The system (2.1)-(2.7), (1.5)-(1.8) is equivalent to

U t = A U , U ( 0 ) = U 0 ; (2.8)

where

𝒰(t) = (u(t), v(t),U(t),V(t), η(t), w(t),W(t))T and 𝒰0 = (u 0, v 0, u 1, v 1, η0, w 0, w 1)T .

Under this conditions, we have

Theorem 2.1. The operator 𝒜 is the infinitesimal generator of a C 0 -semigroup of contractions (S(t)) t ≥ 0 on ℋ. Thus, for any initial data 𝒰 0 ∈ ℋ , the problem (2.8) has a unique weak (mild) solution

U C 0 ( [ 0 , [ , H ) .

Moreover, if 𝒰0D(𝒜 ), then 𝒰 is a strong solution of (2.8), that is

U C 1 ( [ 0 , [ , H ) C 0 ( [ 0 , [ , D ( A ) ) .

Proof. It easy to see that 𝒟(𝒜) is dense in ; and, since 𝒜 is a dissipative operator, it is enought to show that 0ρ(A). To do that, we will show that for F=(f1,f2,,f7)TH, there exists only one U=(u,v,U,V,η,w,W)TD(A) such that 𝒜𝒰 = F. In terms of the components, we have

U = f 1 (2.9)

V = f 2 (2.10)

α ρ 1 u x x + 1 ρ 1 0 g ( s ) η x x ( s ) d s = f 3 (2.11)

α 2 ρ 2 v x x = f 4 (2.12)

U - η s = f 5 (2.13)

W = f 6 = f 2 ( l ) (2.14)

- γ 1 ρ 3 W - γ 2 ρ 3 w - α 2 ρ 3 v x ( l ) = f 7 (2.15)

Indeed, from the equations (2.9) and (2.13), we get that ηsLg2 and that

η x , s = s f 1 x - 0 s f 5 x , τ d τ

which means that η is uniquely determined. Moreover, using (1.10) and (2.5), we can write, for each T > 0:

0 T g ( s ) 0 l 0 | η x ( s ) | 2 d x d s 2 k 2 0 T g ( s ) 0 l 0 η x ( s ) η s x ( s ) ¯ d x d s 1 2 0 T g ( s ) 0 l 0 | η x ( s ) | 2 d x d s + 2 k 2 2 0 T g ( s ) 0 l 0 | η s x ( s ) | 2 d x d s

from we obtain

η L g 2 2 k 2 η s L g 2

which enables us to conclude that ηLg2. Thus, it remains only to establish the existence and uniqueness of solution for the system

( P ) u x x = ρ 1 α f 3 - 1 α 0 g ( s ) η x x ( s ) d s v x x = ρ 2 α 2 f 4 u ( 0 ) = 0 , u ( l 0 ) = v ( l 0 ) , α 2 v x ( l ) + γ 2 v ( l ) = - ρ 3 f 7 - γ 1 f 6 α u x ( l 0 ) - α 2 v x ( l 0 ) = - 0 g ( s ) η x ( l 0 , s ) d s .

Let us consider the functional T: X → ℂ given by

T ( φ , ψ ) : = - ρ 1 0 l 0 f 3 ¯ φ d x - 0 l 0 0 g ( s ) η x ( s ) d s ¯ φ x d x + + ρ 2 l 0 l l 0 x f 4 ( τ ) d τ ¯ ψ x d x + G ~ ¯ ψ ( l )

for all (φ, ψ) ∈ X , where G~:=G-ρ2l0lf4dx and X:=*1 is a Hilbert space, endowed with the inner product

( φ , ψ ) , ( u , v ) X = α 0 l 0 u x ¯ φ x d x + α 2 l 0 l v x ¯ ψ x d x + γ 2 v ( l ) ¯ ψ ( l ) .

It’s clear that T ∈ X′; therefore, by the Riesz representation theorem we conclude that there exists only one weak solution to system (P). So we have that 0 ∈ ρ(𝒜 ). □

3 EXPONENTIAL STABILITY

In this section, we show that if hypothesis (1.9)-(1.11) hold, then the corresponding semigroup is exponentially stable. The main tool we use is Prüss’s results 99 J. Prüss. On the Spectrum of C0-semigroups. Trans. AMS, 284 (1984), 847-857., which is summarized in the following theorem.

Theorem 3.2.Let (𝒮 (t))t≥0 be a C 0 -semigroup on a Hilbert space ℋ generated by 𝒜 . Then the semigroup is exponentially stable if and only if

i ρ ( A ) , and ( i λ I - A ) - 1 L ( H ) C , λ .

In the next Lemma we show that the imaginary axis is contained in the resolvent set.

Lemma 3.1. Under the hypotheses (1.9)-(1.11), the operator 𝒜 verify

i ρ ( A ) . (3.1)

Proof. In the Theorem 2.1, we have already shown that 0 ∈ ρ(𝒜). Moreover, note that we can’t conclude that the spectrum of 𝒜 is formed only by eigenvalues, since 𝒜 −1 is not compact. So, if (3.1) is not true, then there exists λ 0 ℝ with A-1-1|λ0|, such that {iλ;|λ|<|λ0|}ρ(A) and sup{(iλ-A)-1;|λ|<|λ0|}=. Follow that there exist sequences (λ n )n ⊂ ℝ and (Un)n=((un,vn,Un,Vn,ηn,wn,Wn)T)nD(A), such that

λ n | λ 0 | , (3.2)

U n H = 1 , n , (3.3)

( i λ n - A ) U n = F n = ( f n 1 , . . . , f n 7 ) 0 in H . (3.4)

But, (3.4) is equivalent to

i λ n u n - U n = f n 1 0 in H 1 ( 0 , l 0 ) (3.5)

i λ n v n - V n = f n 2 in H 1 ( l 0 , l ) (3.6)

i λ n U n - α ρ 1 u n , x x - 1 ρ 1 0 g ( s ) η n , x x ( s ) d s = f n 3 0 in L 2 ( 0 , l 0 ) (3.7)

i λ n V n - α 2 ρ 2 v n , x x = f n 4 0 in L 2 ( l 0 , l ) (3.8)

i λ n η n - U n + η n , s = f n 5 0 in L g 2 (3.9)

i λ n w n - W n = f n 6 0 in (3.10)

i λ n W n + γ 1 ρ 3 W n + γ 2 ρ 3 w n + α 2 ρ 3 v n , x ( l ) = f n 7 0 in (3.11)

Taking the inner product of (3.4) with 𝒰 n in , we get

Re A U n , U n = - γ 1 | W n | 2 + 1 2 0 g ' ( s ) | η n , x ( s ) | 2 d s d x 0 . (3.12)

Follows from condition (1.10) and (3.12) that

γ 1 | W n | 0 and η n L g 2 0 . (3.13)

Now, we use (3.3) and (3.5) for conclude that there exist U, uL 2(0, l 0) and subsequences still denoted by (U n )n and (u n )n , such that

U n U in L 2 ( 0 , l 0 ) and u n u in L 2 ( 0 , l 0 ) . (3.14)

Moreover, using (3.3), (3.7), (3.13) and (3.14), we get

u n , x u x in L 2 ( 0 , l 0 ) , (3.15)

α u n , x x + 0 g ( s ) η n , x x ( s ) d s α u x x in L 2 ( 0 , l 0 ) .

So, from (3.14), (3.15) and (3.5), we can conclude that

u n u in H 1 ( 0 , l 0 ) and U n U in H 1 ( 0 , l 0 ) .

Proceeding analogously, we find

v n v in H 2 ( l 0 , l ) and V n V in H 1 ( l 0 , l ) .

From convergences above, remembering w n = v n (l) and W n = V n (l), follow that

w n w : = v ( l ) , W n W : = V ( l ) , and v n , x ( l ) v x ( l ) .

The convergences obtained above allow us to pass to the limit in (3.5)-(3.11), obtaining the following system

i | λ 0 | u - U = 0 i | λ 0 | v - V = 0 i | λ 0 | ρ 1 U - α u x x = 0 i | λ 0 | ρ 2 V - α 2 v x x = 0 i | λ 0 | w - W = 0 i | λ 0 | W + γ 1 ρ 3 W + γ 2 ρ 3 w + α 2 ρ 3 v x ( l ) = 0

We conclude that there exists U=(u,v,U,V,0,v(l),V(l))TD(A), such that

U n U in H , (3.16)

where (u, v) is precisely the solution of the system given by

u x x + ρ 1 | λ 0 | 2 α u = 0 in ] 0 , l 0 [ v x x + ρ 2 | λ 0 | 2 α 2 v = 0 in ] l 0 , l [ u ( 0 ) = 0 , u ( l 0 ) = v ( l 0 ) α u x ( l 0 ) = α 2 v x ( l 0 ) , v ( l ) = 0

when γ 1 > 0; or, in the case that γ 1 = 0, of the system obtained in the above system, replacing the boundary condition v(l) = 0 by v x (l) + α 3 v(l) = 0, where α3=α2-1(γ2-|λ0|2ρ3).

However, each one of these two systems has a unique solution, namely, the null solution; from which it follows that 𝒰 = 0. With this, we rewrite (3.16) as

U n 0 in H , (3.17)

which contradicts (3.3), completing the proof. □

From now, on C will denote a generic constant, that can be different in different places. Let us denote by

b : = 0 g ( s ) d s . (3.18)

We will prove that the solution 𝒰 of the resolvent equation

( i λ I - A ) U = F (3.19)

is uniformly bounded for any take F = ( f 1 , f 2 , ..., f 7)T, that is,

U H C F H F H .

In fact, in terms of the components we have

i λ u - U = f 1 (3.20)

i λ v - V = f 2 (3.21)

i λ ρ 1 U - α u x x - 0 g ( s ) η x x ( s ) d s = ρ 1 f 3 (3.22)

i λ ρ 2 V - α 2 v x x = ρ 2 f 4 (3.23)

i λ η - U + η s = f 5 (3.24)

i λ w - W = f 6 (3.25)

i λ ρ 3 W + γ 1 W + γ 2 w + α 2 v x ( l ) = ρ 3 f 7 (3.26)

The dissipative properties of 𝒜 implies that there exists a positive constant C such that

γ 1 | W | 2 + 0 l 0 0 g ( s ) | η x ( s ) | 2 d s d x C U H F H . (3.27)

The following Lemma will play an important role in the sequel.

Lemma 3.2.Under the above notations, for any ε > 0 sufficiently small, there exist a constant C ε > 0 such that, for |λ| large enough, hold

0 l 0 | U | 2 + | u x | 2 d x C ϵ U H F H + C ϵ F H 2 + ϵ α u x ( l 0 ) + 0 g ( s ) η x ( l 0 , s ) d s 2 .

Proof. Multiplying the equation (3.22) by 0g(s)η(s)¯ds and using (3.24), we get

b ρ 1 0 l 0 | U | 2 d x = ρ 1 0 l 0 0 g ( s ) η s ( s ) ¯ U d s d x - ρ 1 0 l 0 0 g ( s ) f 5 ¯ U d s d x + + α 0 l 0 0 g ( s ) η x ( s ) ¯ u x d s d x + 0 l 0 0 g ( s ) η x ( s ) d s 2 d x - ρ 1 0 l 0 0 g ( s ) η ( s ) ¯ f 3 d s d x + - α u x ( l 0 ) + 0 g ( s ) η x ( l 0 , s ) d s 0 g ( s ) η ( l 0 , s ) ¯ d s : = R 1 .

For each ε > 0, we use (1.10) and (3.27) for obtain

Re 0 l 0 0 g ( s ) η s ( s ) ¯ U d s d x ϵ U L 2 ( 0 , l 0 ) 2 + C ϵ 0 l 0 0 g ( s ) | η x ( s ) | 2 d s ϵ U L 2 ( 0 , l 0 ) 2 + C ϵ U H F H .

Using (3.27) once more we get

Re α 0 l 0 0 g ( s ) η x ( s ) ¯ u x d s d x ϵ u x L 2 ( 0 , l 0 ) 2 + C ϵ U H F H

and

| R 1 | C ϵ U H F H + ϵ α u x ( l 0 ) + 0 g ( s ) η x ( l 0 , s ) d s 2 .

Therefore, for ε > 0 sufficiently small, we have

0 l 0 | U | 2 d x C ϵ U H F H + ϵ 0 l 0 | u x | 2 d x + ϵ α u x ( l 0 ) + 0 g ( s ) η x ( l 0 , s ) d s 2 . (3.28)

On the other hand, multiplying (3.22) by u¯ and using (3.20), we get

α 0 l 0 | u x | 2 d x = ρ 1 0 l 0 | U | 2 d x + ρ 1 0 l 0 U f 1 ¯ d x + ρ 1 0 l 0 f 3 u ¯ d x + - 0 l 0 0 g ( s ) η x ( s ) u x ¯ d s d x + α u x ( l 0 ) + 0 g ( s ) η x ( l 0 , s ) d s u ( l 0 ) ¯ : = R 2 .

Since

| u ( l 0 ) | u L C u L 2 1 / 2 u x L 2 1 / 2 ,

using (3.20) we get, for each ε > 0 and for λ ≠ 0:

| R 2 | ϵ α u x ( l 0 ) + 0 g ( s ) η x ( l 0 , s ) d s 2 + C ϵ | u ( l 0 ) | 2 ϵ α u x ( l 0 ) + 0 g ( s ) η x ( l 0 , s ) d s 2 + ϵ u x L 2 2 + C ϵ u L 2 2 ϵ α u x ( l 0 ) + 0 g ( s ) η x ( l 0 , s ) d s 2 + ϵ u x L 2 2 + C ϵ | λ | U L 2 2 + C ϵ | λ | F H 2 .

Therefore, for ε > 0 sufficiently small, we have

0 l 0 | u x | 2 d x C ϵ , λ 0 l 0 | U | 2 d x + C ϵ U H F H + C ϵ , λ F H 2 + + ϵ α u x ( l 0 ) + 0 g ( s ) η x ( l 0 , s ) d s 2 .

From above inequality and (3.28), our conclusion follows. □

The next lemma is crucial to ensure that the exponential decay occurs in the case where γ 1 = 0.

Indeed, it provides an estimate for the term involving |W| 2 that can be obtained from (3.27) only when γ 1 is positive.

Lemma 3.3.There exist C > 0 such that

ρ 2 | W | 2 C U H F H + C 0 l 0 α | v x | 2 d x .

Proof. Multiplying (3.23) by (x-l0)vx¯, using (3.21), and remembering that V (l) = W , we get

ρ 2 | W | 2 + α 2 | v x ( l ) | 2 = - ρ 2 l - l 0 l 0 l | V | 2 d x - 2 ρ 2 l - l 0 l 0 l ( x - l 0 ) V f x 2 ¯ d x + + α 2 l - l 0 l 0 l | v x | 2 d x - 2 ρ 2 l - l 0 l 0 l ( x - l 0 ) f 4 v x ¯ d x .

Taking the real part, we get the desired inequality. □

Now we are in condition to show the main result of this section.

Theorem 3.3. Let us suppose that (1.9)-(1.11) hold. Then the semigroup e 𝒜t is exponentially stable.

Proof. In view of Proposition 1, we only need to show that there exist C > 0 such that:

( i λ I - A ) - 1 L ( H ) C , λ .

Since the resolvent operator is holomorphic, it is enough to prove the above inequality for |λ| large enough. In fact, multiplying (3.23) by (l-x)vx¯ and using (3.21), we get

l 0 l ρ 2 | V | 2 + α 2 | v x | 2 d x = ( l - l 0 ) ρ 2 | V ( l 0 ) | 2 + α 2 | v x ( l 0 ) | 2 + - 2 ρ 2 l 0 l ( l - x ) V f x 2 ¯ d x - 2 ρ 2 l 0 l ( l - x ) f 4 v x ¯ d x .

Taking the real part, we get

l 0 l ρ 2 | V | 2 + α 2 | v x | 2 d x C α 2 | v x ( l 0 ) | 2 + ρ 2 | V ( l 0 ) | 2 + C U H F H . (3.29)

On the other hand, multiplying (3.22) by xαux+0g(s)ηx(s)ds¯, and using (3.20) and (3.24), we get

ρ 1 | U ( l 0 ) | 2 + α u x ( l 0 ) + 0 g ( s ) η x ( l 0 , s ) d s 2 C 0 l 0 α | u x | 2 + ρ 1 | U | 2 d x + C U H F H .

From Lemma 3.2, we get, for ε small enough and for |λ| large enough, that

ρ 1 | U ( l 0 ) | 2 + α u x ( l 0 ) + 0 g ( s ) η x ( l 0 , s ) d s 2 C U H F H + C F H 2 . (3.30)

Using the transmission conditions, inequality (3.29) can be estimated by (3.30), that is

l 0 l α 2 | v x | 2 + ρ 2 | V | 2 d x C U H F H + C F H 2 . (3.31)

Moreover from Lemma 3.2 and inequality (3.30), we get

0 l 0 ρ 1 | U | 2 + α | u x | 2 d x C U H F H + C F H 2 .

Therefore, Lemma 3.3, equation (3.25), and inequality (3.31) implies

| W | 2 + | w | 2 C U H F H + C F H 2 .

From the three last inequalities and (3.27), we get

U H 2 C U H F H + C F H 2 ,

which implies in

( i λ I - A ) - 1 F H = U H C F H .

Then, the semigroup is exponentially stable. □

4 THE LACK OF EXPONENTIAL STABILITY

Now we shall prove that the dissipation given by the memory effect is necessary for exponential stability of the system. Let us consider the problem without memory effect; namely

ρ 1 u t t - α 1 u x x = 0 in ] 0 , l 0 [ × ] 0 , + [ (4.1)

ρ 2 v t t - α 2 v x x = 0 in ] l 0 , l [ × ] 0 , + [ (4.2)

ρ 3 w t t + γ 1 w t + γ 2 w + α 2 v x ( l ) = 0 in ] 0 , + [ (4.3)

with boundary conditions

u ( 0 ) = 0 , v ( l ) = w in ] 0 , + [ (4.4)

and with transmission conditions

u ( l 0 ) = v ( l 0 ) , α 1 u x ( l 0 ) = α 2 v x ( l 0 ) in ] 0 , + [ (4.5)

and initial data

( u ( 0 ) , v ( 0 ) , u t ( 0 ) , v t ( 0 ) , w ( 0 ) , w t ( 0 ) ) = ( u 0 , v 0 , u 1 , v 1 , w 0 , w 1 ) H ˘ , (4.6)

where α 1, α 2, ρ 1, ρ 2, ρ 3, γ 2 are as before, and γ 1, now, is a positive constant. Moreover, for this problem, we consider the phase space

H ˘ = U = ( u , v , U , V , w , W ) T * 1 × L 2 × 2 ; v ( l ) = w .

The total energy associed with the system is

E ˘ ( t ) = 1 2 0 l 0 ρ 1 | u t | 2 + α 1 | u x | 2 d x + l 0 l ρ 2 | v t | 2 + α 2 | v x | 2 d x + ρ 3 | w t | 2 + γ 2 | w | 2

and it is not difficult to see that, for all 𝒰H˘, we have

d d t E ˘ ( t ) = - γ 1 | w t | 2 .

Let us denote by the unbounded operator of H˘ given by

B U = U V α 1 ρ 1 u x x α 2 ρ 2 v x x W - γ 1 ρ 3 W - γ 2 ρ 3 w - α 2 ρ 3 v x ( l )

with domain

D ( B ) = U = ( u , v , U , V , w , W ) T ( * 1 2 ) × * 1 × 2 ; V ( l ) = W , α 1 u x ( l 0 ) = α 2 v x ( l 0 ) .

It is easy to see that

Re ( B U , U ) H ˘ = - γ 1 | W | 2 . (4.7)

It is not difficult to see that the operator is the infinitesimal generator of a C 0-semigroup of contractions over H˘, which we will denote by T(t). This shows that the problem (4.1)-(4.6) is well-posed.

To prove that the system (4.1)-(4.6) is not exponentially stable, the main tool we use is the Weyl’s theorem about the invariance of the essential spectral radius by compact perturbations. To do that, let us consider the following conservative system

ρ 1 u ~ t t - α 1 u ~ x x = 0 in ] 0 , l 0 [ × ] 0 , + [ (4.8)

ρ 2 v ~ t t - α 2 v ~ x x = 0 in ] l 0 , l [ × ] 0 , + [ (4.9)

ρ 3 w ~ t t + γ 2 w ~ + α 2 v ~ x ( l ) = 0 in ] 0 , + [ (4.10)

verifying the same boundary and transmission conditions and with the same initial data, where α 1, α 2, ρ 1, ρ 2, ρ 3 and γ 2 are as before. That is, with boundary conditions

u ~ ( 0 ) = 0 , v ~ ( l ) = w ~ in ] 0 , + [ (4.11)

and transmission conditions

u ~ ( l 0 ) = v ~ ( l 0 ) , α 1 u ~ x ( l 0 ) = α 2 v ~ x ( l 0 ) in ] 0 , + [ (4.12)

and initial data

( u ~ ( 0 ) , v ~ ( 0 ) , u ~ t ( 0 ) , v ~ t ( 0 ) , w ~ ( 0 ) , w ~ t ( 0 ) ) = ( u 0 , v 0 , u 1 , v 1 , w 0 , w 1 ) H ˘ . (4.13)

The total energy associed with this system is

E ~ ( t ) = 1 2 0 l 0 ρ 1 | u ~ t | 2 + α 1 | u ~ x | 2 d x + l 0 l ρ 2 | v ~ t | 2 + α 2 | v ~ x | 2 d x + ρ 3 | w ~ t | 2 + γ 2 | w ~ | 2

and it is not difficult to see that

d d t E ~ ( t ) = 0 .

Therefore the system is conservative and there is no decay. Now we are in conditions to show the main result of this section.

Theorem 4.4.The semigroup T(t) associated to system (4.1)-(4.6) is not exponentially stable.

Proof. The main ideia is to prove that T(t) have the same essential spectral radius of the semi- group associated to conservative system (4.8)-(4.13), that we denote as T 0(t). Here, we use the Weyl’s Theorem (Theorem XIII.14, 1010 M. Reed & B. Simon. “Methods of Modern Mathematical Physics”, volume 04. Academic Press Inc., California (1978).; see also Kato’s book 88 T. Kato. “Perturbation Theory for Linear Operators”. Springer-Verlag, New York (1980)., Theorem 5.35, p. 244 for details of the proof), which stablish that if the difference of two operators is compact, then the your essential spectrum radii are equals. More precisely

Theorem 4.5. Let S and T two continuous operator over a Banach space X. If S−T is a compact operator, then S and T have the same essential spectrum radius.

So, we will show that the difference T (t) −T 0(t) is a compact operator; from which we obtain

ω e s s ( T ) = ω e s s ( T 0 ) .

But since T 0(t) is unitary, then ω ess (T 0) = 0. Denoting by ω(T) and ω σ () the type of semi- group T(t) and the spectral upper bound of spectrum σ() respectively, we have that (see 55 K.J. Engel & R. Nagel. “One-parameter Semigroups for Linear Evolution Equations”, volume 194 of Graduate Texts in Mathematics. Springer-Verlag, New York (2000)., Corollary 2.11, p. 258):

ω ( T ) = max ω σ ( B ) , ω e s s ( T ) = 0 .

This imply that T(t) is not exponentially stable. In fact. Let (u, v, w) and (u~,v~,w~) be the solutions of the systems (4.1)-(4.6) and (4.8)-(4.13), respectively. Denoting by

U : = u - u ~ , V : = v - v ~ , W : = w - w ~ ,

we have that (U,V,W) is solution of the system

ρ 1 U t t - α 1 U x x = 0 in ] 0 , l 0 [ × ] 0 , + [ (4.14)

ρ 2 V t t - α 2 V x x = 0 in ] l 0 , l [ × ] 0 , + [ (4.15)

ρ 3 W t t + γ 1 w t + γ 2 W + α 2 V x ( l ) = 0 in ] 0 , + [ (4.16)

with boundary conditions

U ( 0 ) = 0 , V ( l ) = W in ] 0 , + [ , (4.17)

and transmission conditions

U ( l 0 ) = V ( l 0 ) , α 1 U x ( l 0 ) = α 2 V x ( l 0 ) in ] 0 , + [ , (4.18)

and initial data

( U ( 0 ) , V ( 0 ) , U t ( 0 ) , V t ( 0 ) , W ( 0 ) , W t ( 0 ) ) = ( 0 , 0 , 0 , 0 , 0 , 0 ) H ˘ . (4.19)

The energy associated with the system (4.14)-(4.19) is given by

E ( t ) = 1 2 0 l 0 ρ 1 | U t | 2 + α 1 | U x | 2 d x + l 0 l ρ 2 | V t | 2 + α 2 | V x | 2 d x + ρ 3 | W t | 2 + γ 2 | W | 2 .

It is easy to verify that

d d t E ( t ) + γ 1 | W t | 2 = - γ 1 w ~ t W ¯ t ,

from where follows

E ( t ) + γ 1 0 t | W t | 2 d s = - γ 1 0 t w ~ t W ¯ t d s . (4.20)

Now, let us denote by U0,n:=(u0,n,v0,n,u1,n,v1,n,w0,n,w1,n)T a bounded sequence of initial data in the phase space H˘. We will show that the corresponding sequence of solutions Un:=(Un,Vn,Un,t,Vn,t,Wn,Wn,t)T has a subsequence that converges strongly in H˘.

To show this, note that T(t)?? 0,n and T 0(t)𝒰 0,n are bounded in H˘. This implies that v~n,x(l) is bounded in L 2(0, T), for all T > 0. Therefore (4.16) implies that W n,t is bounded in H 1(0, T). Since H 1(0, T) has compact embedding in L 2(0, T), it follows that there exist subsequences, we still denote as (W n )n and (w~n)n such that

W n , t W t strongly in L 2 ( 0 , T ) ,

and similarly

w ~ n , t w ~ t strongly in L 2 ( 0 , T ) .

From the above convergences we have

0 T w ~ n , t W ¯ n , t d t 0 T w ~ t W ¯ t d t . (4.21)

Using this convergence in (4.20) it follows that [T(t)-T0(t)]U0,nH˘ converges, which implies that [T(t) − T 0(t)]𝒰 0,n converges strongly in H˘. This means that T(t) − T 0(t) is a compact operator in H˘, and therefore the proof is complete. □

5 POLYNOMIAL DECAY

In this section we show that the solutions of the system (4.1)-(4.6) decays polynomially to zero as t −1/2 . To show this, we use the Borichev and Tomilov’s Theorem (see 22 A. Borichev & Y. Tomilov. Optimal Polynomial Decay of Functions and Operator Semigroups. Math. Ann., 347 (2009), 455-478.):

Theorem 5.6.Let S(t) be a bounded C0-semigroup on a Hilbert space ℋ with generator 𝒜 such that iℝ ⊂ ρ(𝒜). Then

1 | λ | β ( i λ I - A ) - 1 L ( H ) C , λ S ( t ) A - 1 D ( A ) C t 1 / β , t > 0 .

Our starting point is to study the solution of the resolvent equation

i λ U - B U = F

which is rewritten, in terms of the components, in

i λ u - U = f 1 (5.1)

i λ v - V = f 2 (5.2)

i λ U - α 1 ρ 1 u x x = f 3 (5.3)

i λ V - α 2 ρ 2 v x x = f 4 (5.4)

i λ w - W = f 5 (5.5)

i λ W + γ 1 ρ 3 W + γ 2 ρ 3 w + α 2 ρ 3 v x ( l ) = f 6 . (5.6)

Taking the inner product with 𝒰 and using (4.7) it follows that

γ 1 | W | 2 C U F . (5.7)

Lemma 5.4.For |λ| large enough, there exist C > 0 such that

l 0 l [ ρ 2 | V | 2 + α 2 | v x | 2 ] d x C | λ | 2 U F + C F 2 .

Proof. Multiplying (5.4) by (x-l0)v¯x and using (5.2), we have

1 2 l 0 l [ ρ 2 | V | 2 + α 2 | v x | 2 ] d x = ( l - l 0 ) 2 ρ 2 | W | 2 + α 2 | v x ( l ) | 2 + ρ 2 l 0 l ( x - l 0 ) f 4 v ¯ x + V f x 2 ¯ d x . (5.8)

On the other hand, using (5.6) we get

α 2 | v x ( l ) | 2 C | λ | 2 U F + C F 2 . (5.9)

Finally, from (5.7)-(5.9) our conclusion follows. □

Lemma 5.5.For |λ| large enough, there exist C > 0 such that

0 l 0 [ ρ 1 | U | 2 + α 1 | u x | 2 ] d x C | λ | 2 U F + C F 2 .

Proof. Multiplying equation (5.3) by xu¯x, using (5.1) and the transmission conditions, we get

1 2 0 l 0 [ ρ 1 | U | 2 + α 1 | u x | 2 ] d x = l 0 2 ρ 1 | V ( l 0 ) | 2 + α 2 2 α 1 | v x ( l 0 ) | 2 + ρ 1 0 l 0 x f 3 u ¯ x + U f x 1 ¯ d x . (5.10)

On the other hand, multiplying equation (5.4) by (l-x)v¯x and using (5.2) we get

( l - l 0 ) 2 ρ 2 | V ( l 0 ) | 2 + α 2 | v x ( l 0 ) | 2 = 1 2 l 0 l [ ρ 2 | V | 2 + α 2 | v x | 2 ] d x + + ρ 2 l 0 l ( l - x ) f 4 v ¯ x + V f x 2 ¯ d x .

From this, and from Lemma 5.4 there exist C > 0, such that

ρ 2 | V ( l 0 ) | 2 + α 2 | v x ( l 0 ) | 2 C | λ | 2 U F + C F 2 . (5.11)

Combining (5.10) and (5.11), our conclusion follows. □

Now we are able to stablish the main result of this section.

Theorem 5.7.The semigroup T(t) associated to the transmission problem with load tip (4.1)-(4.6) decays polynomially as t−1/2. Moreover, if 𝒰0D( k ), then

T ( t ) U 0 H ˘ C k t k / 2 U 0 D ( B k ) .

Proof. The proof of iℝ ⊂ ρ() is analogous to the proof of Lemma 3.1. And, from Lemmas 5.4 and 5.5, follows that, for |λ| large enough, we have

U H ˘ 2 C | λ | 2 U H ˘ F H ˘ + C F H ˘ 2 ,

and this implies that

U H ˘ 2 C | λ | 4 F H ˘ 2 ,

which gives us, for |λ| large enough

( i λ I - B ) - 1 F H ˘ = U H ˘ C | λ | 2 F H ˘ .

Therefore, from Theorem 5.6 our conclusion follows. □

REFERENCES

  • 1
    K. Andrews & M. Shillor. Vibrations of a Beam With a Damping Tip Body. Mathematical and Computer Modelling, 35 (2002), 1033-1042.
  • 2
    A. Borichev & Y. Tomilov. Optimal Polynomial Decay of Functions and Operator Semigroups. Math. Ann., 347 (2009), 455-478.
  • 3
    C.M. Dafermos. Asymptotic Stability in Viscoleasticity. Arch. Rat. Mech. Anal., 37 (1970), 297-308.
  • 4
    C.M. Dafermos. On Abstract Volterra Equation with Applications to Linear Viscoleasticity. Differential and Integral Equations, 7 (1970), 554-569.
  • 5
    K.J. Engel & R. Nagel. “One-parameter Semigroups for Linear Evolution Equations”, volume 194 of Graduate Texts in Mathematics. Springer-Verlag, New York (2000).
  • 6
    M. Fabrizio & A. Morro. “Asymptotic Stability in Viscoleasticity”, volume 12 of SIAM Studies in Applied Mathematics. Society for Industrial and Applied Mathematics, Philadelphia (1992).
  • 7
    E. Feireisl & G. O’Dowd. Stabilisation d’un syste`me hybride par un feedback non linéaire, non monotone. Comptes Rendus de l’Académie des Sciences, 326 (1998), 323-327.
  • 8
    T. Kato. “Perturbation Theory for Linear Operators”. Springer-Verlag, New York (1980).
  • 9
    J. Prüss. On the Spectrum of C0-semigroups. Trans. AMS, 284 (1984), 847-857.
  • 10
    M. Reed & B. Simon. “Methods of Modern Mathematical Physics”, volume 04. Academic Press Inc., California (1978).
  • 11
    N. Zietsman, L Van Rensburg & A. Van der Merwe. A Timoshenko Beam With Tip Body and Boundary Damping. Wave Motion, 39 (2004), 199-211.

Publication Dates

  • Publication in this collection
    29 May 2023
  • Date of issue
    Apr-Jun 2023

History

  • Received
    17 Dec 2021
  • Accepted
    11 Oct 2022
Sociedade Brasileira de Matemática Aplicada e Computacional - SBMAC Rua Maestro João Seppe, nº. 900, 16º. andar - Sala 163, Cep: 13561-120 - SP / São Carlos - Brasil, +55 (16) 3412-9752 - São Carlos - SP - Brazil
E-mail: sbmac@sbmac.org.br