Acessibilidade / Reportar erro

A transmission problem for the Timoshenko system

Abstract

In this work we study a transmission problem for the model of beams developed by S.P. Timoshenko [10]. We consider the case of mixed material, that is, a part of the beam has friction and the other is purely elastic. We show that for this type of material, the dissipation produced by the frictional part is strong enough to produce exponential decay of the solution, no matter how small is its size. We use the method of energy to prove exponential decay for the solution.

transmission; Timoshenko; beams; exponential decay; frictional damping


A transmission problem for the Timoshenko system

C.A. RaposoI; W.D. BastosII; M.L. SantosIII

IDepartment of Mathematics, UFSJ, Praça Frei Orlando, 170 – 36307-352 São João del-Rei, MG

IIDepartment of Mathematics, UNESP, Rua Cristovão Colombo, 2265 – 15054-000 São José do Rio Preto, SP

IIIDepartment of Mathematics, UFPA, Rua do Una, 156 - 66113-200 Pará, PA E-mails: raposo@ufsj.edu.br / waldemar@ibilce.unesp.br / mlsantos@ufpa.br

ABSTRACT

In this work we study a transmission problem for the model of beams developed by S.P. Timoshenko [10]. We consider the case of mixed material, that is, a part of the beam has friction and the other is purely elastic. We show that for this type of material, the dissipation produced by the frictional part is strong enough to produce exponential decay of the solution, no matter how small is its size. We use the method of energy to prove exponential decay for the solution.

Mathematical subject classification: 35J55, 35J77, 93C20.

Key words: transmission, Timoshenko, beams, exponential decay, frictional damping.

1 Introduction

The transverse vibration of a beam is mathematically described by a system of two coupled differential equations given by

Here, L is the length of the beam in its equilibrium position, t is the time variable and x is the space coordinate along the beam. The function u = u(x,t) is the transverse displacement of the beam and y = y(x,t) is the rotation angle of a filament of the beam. The coefficients r, I r, E, I and K are the mass per unit length, the polar moment of inertia of a cross section, Young's modulus of elasticity, the moment of inertia of a cross section and the shear modulus respectively. We denote r1 = r, r2 = I r, b = EI, k = K and we obtain directly from (1.1) the following system

The mathematical model describing the vibrations of beam with fixed extremities is formed by the system (1.2), boundary conditions

u(0, t) = u(L, t) = y(0, t) = y(L, t) = 0, t > 0,

and initial data

u(·, 0) = f0, ut (·, 0) = f1, y(·, 0) = y0, yt (·, 0) = y1, in (0, L).

If friction is taken into account, the system (1.2) becomes

where a and b are positive constants (we assume a = b = 1). The terms a ut and b y t represent the attrition acting in the vertical vibrations and in the angle of rotation of the filaments of the beam, respectively.

Dissipative properties associated to the system (1.3) have been studied by several authors by considering dissipative mechanism of frictional or viscoelastic type. The frictional dissipation, obtained by introduction of a frictional mechanism acting on the entire domain or on the boundary, was studied in [7, 8, 9]. The viscoelastic dissipation, given by a memory effect, was considered in [2] and [6].

An interesting problem comes out when the dissipation acts only on a part of the domain. In the present paper we consider a frictional mechanism acting only on the part of the domain given by 0 < x < L0 with 0 < L0 < L. We prove that for every L0 the energy of the system decays exponentially to zero as time goes to infinity. In other words, our result states that dissipative properties of the system are transferred to the whole beam and stabilizes the system. The main result of this paper is Theorem 2 and its corollary, both in section 5.

The mathematical model which deals with this situation is called a transmission problem. From the mathematical point of view a transmission problem consist of an initial and boundary value problem for a hyperbolic equation for which the corresponding elliptic operator has discontinuous coefficients. Hence, we can not expect to have regular solutions in the role domain. In the next section we establish the transmission problem and define appropriately the notion of solution considered. We use Hm and Lp to denote the usual Sobolev and Lebesgue spaces [1].

2 The transmission problem

In this section we describe precisely the transmission problem treated in the paper and establish existence and regularity of solution. We begin by introducing the notation

Using the notation above, model (1.3) can be written in the following form:

with boundary conditions,

u(0, t) = v(L, t) = y(0, t) = f(L, t) = 0, t > 0,

transmission conditions,

and initial data

We define the notion of weak solution to the system (2.1)-(2.6) as follows:

Definition 1.Let, m and2be the spaces defined by

We say that (u,v,y,f) is a weak solution to the problem (2.1)-(2.6) if for every (,w) Î (0,T;2Ç) we have:

and

The transmission problem for a single hyperbolic equation was studied by Dautray and Lions [3], who proved the existence and regularity of solutions for the linear problem. The existence and regularity of solutions to the transmission problem for the Timoshenko system is given in the following theorem:

Theorem 1. If (u0,v0),(y0,f0) Î and (u1,v1),(y1,f1) Î 2, then there exists a unique weak solution (u,v,y,f) to the system (2.1)-(2.6) satisfying:

(u, v), (y, f) Î C (0,¥;) Ç C1(0,¥;2).

Moreover, if (u0,v0),(y0,f0) Î 2Ç and (u1,v1),(y1,f1) Î , then the weak solution is a strong solution and satisfies

(u, v), (y, f) Î C(0,¥;2Ç ) Ç C1(0,¥;) Ç C2(0,¥;2).

Proof. For the proof we proceed in a quite similar manner as in [3].

The total energy associated to the system is defined by

Next we prove that the total energy associated to the system is decreasing for every t > 0.

Lemma 1. Let (u,v,y,f) be the strong solution to the system (2.1)-(2.6), then

Proof. Multiplying (2.1) by ut and integrating by parts over the interval (0,L0), we get

Now, multiplying (2.2) by yt and integrating by parts over (0,L0) we obtain

Multiplying (2.3) by vt and integrating by parts on (L0,L), we get

Multiplying (2.4) by ft and integrating by parts on (L0,L) leads to

Now observe that

and

Summing up (2.7), (2.8), (2.9) and (2.10), and using (2.11) and (2.12) together with the hypothesis of transmission we obtain

3 Technical lemmas

Now we develop a series of technical results in order to facilitate the proof of the main result of the paper. We begin by constructing a functional (t), equivalent to the energy functional, which satisfies (t) < C (0), C < 1. In order to do so, we use some multiplier techniques (usually associated to control problems) and the following restrictions on the boundary conditions for the elastic part of the beam:

Lemma 2. Let us define

Then

Proof. Multiply (2.1) by x ux and integrate by parts over (0,L0) to get

Multiply (2.2) by x yx and integrate by parts over (0,L0) to obtain

Multiplying (2.3) by x vx and integrating by parts over (L0,L) leads to

Now multiply (2.4) by x fx and integrate over (0,L0) to obtain

Summing up (3.3), (3.4), (3.5), (3.6), and making use of the hypothesis on the transmission, the punctual terms are canceled and we get

Now using (3.1), (3.2) and the Young's inequality [4] we obtain

Now, if we define

and choose N1 > we conclude

It is worth noticing that the estimate above is important in two aspects. First, it recovers a part of energy with minus sign. Second, it will play a fundamental role in the next two lemmas controlling punctual terms which will come up in the search for other negative terms of the energy.

Lemma 3. Define

Then

Proof. Multiply (2.1) by x (ux+y) and integrate by parts over (0,L0) to obtain

Using Young's inequality, we get

Now, defining

we obtain

If we choose N2 > we conclude that there exists > 0 such that

and then

Lemma 4. Let 3 be defined as

Then

Proof. Multiply (2.3) by (x-L)(vx+f) and integrate by parts over (L0,L) to get

Using Young's inequality we obtain

Now we set

and verify that

It follows then

Observe that in the attempt to recover the total energy of the system with negative sign we introduced the lemmas 3 and 4. Now we need to control them. It will be achieved with the aid of the next section.

4 Compactness

This section is dedicated to discuss the argument of compactness employed in the proof of the main result of the paper. First we introduce a notation; the symbol is used to denote convergence in the norm of the Sobolev space L¥ as in [5]. For sake of completeness, we state the following result due to J.U. Kim.

Lemma 5. Let (uk) be a sequence of functions satisfying

uku in L¥(0, T, Hb(0, L)),

ut in L2(0, T, Ha(0, L)),

as k ® ¥, with a < b. Then

uk ® u in C([0, T ], H r (0, L)),

for some r < b.

Proof. See [5].

Lemma 6 (Lemma of compactness). If we define

then, for every h > 0 there exists a constant C h > 0 independent of the initial data, such that

for every strong solution (u,v,y,f) to the system (2.1)-(2.6) and sufficiently large T.

Proof. We use a contradiction argument. Define

Suppose that there exists a sequence of initial data

and a positive constant h0 > 0 such that the corresponding solution (un,y n), (v n,f n) of the problem

with boundary conditions,

u n(0, t) = v n(L, t) = y n(0, t) = f n(L, t) = 0, t > 0,

transmission conditions,

and initial data

satisfies

and the following inequality

Then the integral

En(t)dt is bounded for every n Î N,

and also,

and

Now we observe that En(t) > 0 and that En(t)dt is bounded. Hence En(t) is bounded and we can take a subsequence of (un,y n),(vn,f n) (for which we use the same notations) such that

Applying the lemma 5 we conclude that for r < 1

It follows from (4.1) that

We observe that the convergences (4.2) and (4.3) result in

Now, applying Poincare's inequality we obtain

This estimates implies

B(L0, t) = 0,

which is a contradiction to (4.4). This completes the proof of the lemma.

We are now ready to prove the main result of this paper, that is, the exponential decay of the energy associated to the transmission problem for the Timoshenko System with frictional dissipation. This is the content of the next section.

5 Exponential decay

Theorem 2. Let (u,y,v,f) be a strong solution to the transmission problem for the Timoshenko System defined by (2.1)-(2.6). Then there exist positive constants C and w such that

E(t) <C E(0)e-w t .

Proof. We start defining

(t) = N3

1(t) + 2(t) + 3(t).

It follows from lemmas 2, 3, and 4

Now, integrating this inequality over (0,T) and using the Lemma of Compactness we obtain

If we choose and fix h < C0 and N3 such that N3C1 > C h, we get

Since E(t) decreases, we have

Using (5.2) in (5.1) we obtain

(T) - (0) < - T C2(T).

Now observe that for sufficiently large N we have

from what follows that

(T) - (0) < - T (T),

or else

Note that a does not depend on the initial data, and hence, by using the semigroup property we have

For t > 0, there exists a natural n and a real r, 0 < r < T such that t = nT+r. This is equivalent to

Now, using the inequalities (5.3) and (5.4) n times we obtain

Observing once more that E(t) decreases we have

where w = -ln().

Finally, using (5.3) we obtain

E(t) < 4a-1E(0) e-w t ,

and conclude the proof.

We can extend the previous theorem to the weak solutions by using simple density argument and the laws of semi-continuity for the energy functional. In this direction we have the following corollary.

Corollary 1. Under the hypothesis of the previous theorem, there exist positive constants C and w, such that

E(t) <C E(0)e-w t ,

for every weak solution (u,y,v,f) of the system (2.1)-(2.6).

6 Concluding remarks

During the past several decades, many authors have studied the same physical phenomenon for the Timoshenko system formulated into different mathematical models. Our approach to this problem is important not only from mathematical but mainly from the physical point of view with applications in Mechanics, amongst other branches of science. The system studied here is a model for vibrating beams subjected to two frictional mechanisms. More precisely, we proved that the presence of two frictional damping acting in a natural way on a small part of the beam, is enough to stabilize the whole beam. Moreover, it stabilizes quickly (at exponential rate). To the best of our knowledge, our result is the first in this direction. In this sense, this work generalizes the results previously obtained for Timoshenko's system where attrition acting in the whole beam was considered.

Received: 23/I/07. Accepted: 09/II/07.

#668/06.

  • [1] R.A. Adams, Sobolev Spaces, Academic Press, New York, (1975).
  • [2] F. Ammar-Khodja, A. Benabdallah, J.E.M. Rivera and R. Racke, Energy decay for Timoshenko systems of memory type. J. Differential Equations, 194(1) (2003), 82-115.
  • [3] R. Dautray and J.L. Lions, Analyse mathématique et calcul numérique pour les sciences et les techniques, Vol. 1, Masson S.A., Paris (1984).
  • [4] G.H. Hardy, J.E. Littlewood and G. Polya, "A Theorem of W. H. Young" § 8.3 in Inequalities, 2nd ed., Cambridge University Press, Cambridge, (1988), pp. 198-200.
  • [5] J.U. Kim, A Boundary Thin Obstacle Problem for a Wave Equation. Comm. Partial Differential Equations, 14 (8 & 9) (1989), 1011-1026.
  • [6] Z. Liu and C. Peng, Exponential stability of a viscoelastic Timoshenko Beam Adv. Math. Sci. Appl., 8 (1998), 343-351.
  • [7] C.A. Raposo, J. Ferreira, M.L. Santos and N.N.O. Castro, Exponential Stability for the Timoshenko System With Two Weak Damping Appl. Math. Lett., 18 (2005), 535-541.
  • [8] Dong-Hua Shi and De-Xing Feng, Exponential decay of Timoshenko beam with locally distributed feedback IMA J. Math. Control Inform., 18 (2001), 395-403.
  • [9] Dong-Hua Shi, De-Xing Feng and Qing-Xu Yan, Feedback stabilization of rotating Timoshenko beam with adaptive gain. Int. J. Control, 74(3) (2001), 239-251.
  • [10] S.P. Timoshenko and J.M. Gere, Mechanics of Materials, D. Van Nostrand Company, Inc, New York (1972).

Publication Dates

  • Publication in this collection
    24 July 2007
  • Date of issue
    2007

History

  • Accepted
    09 Feb 2007
  • Received
    23 Jan 2007
Sociedade Brasileira de Matemática Aplicada e Computacional Sociedade Brasileira de Matemática Aplicada e Computacional - SBMAC, Rua Maestro João Seppe, nº. 900 , 16º. andar - Sala 163, 13561-120 São Carlos - SP Brasil, Tel./Fax: 55 16 3412-9752 - São Carlos - SP - Brazil
E-mail: sbmac@sbmac.org.br