arxiv: v2 [math.ap] 5 Jun 2015

Transcription

1 arxiv: v2 [math.ap] 5 Jun 215 STABILIZATION OF TRANSVERSE VIBRATIONS OF AN INHOMOGENEOUS EULER-BERNOULLI BEAM WITH A THERMAL EFFECT OCTAVIO VERA, AMELIE RAMBAUD, AND ROBERTO E. ROZAS Abstract. We consier an inhomogeneous Euler-Bernoulli(EB) beam of length L clampe at both ens an subject to : an external frictional amping an a thermal effect. We prove the well-poseness of the moel an analyze the behavior of the solution as t +. The existence is prove using semigroup theory, an the exponential stabilization of solutions is obtaine consiering multiplier technique. A numerical illustration of the energy ecay is given, base on initial ata close to a real physical experiment. The results generalize the ones obtaine where there is no heat conuction taken into account (see [2, 3, 12] for EB type moels, [15, 16] for wave like equations), together with the couple systems of type ampe wave equation with heat conuction [18, 19]. Keywors: semigroup. Asymptotic behavior. Euler-Bernoulli. heat conuction. amping. Couple system. multiplier. 1. Introuction Simplifie moels of beams vibrations have intereste researchers for ecaes. One of the main issues, both from a physical an mathematical point of view for such problems is the question of stabilization of the oscillations when time goes to infinity. There exist various moels of transverse vibrations of beams. Here, we focus on a Euler-Bernoulli(EB) like moel, which is the most commonly use moel because it is simple an provies reasonable engineering approximations. See for example [6, 14, 21, 24, 17]. However, the stanar homogeneous EB equation can be generalize in orer to account for inhomogeneity of the material, heat conuction, an ampings (either viscous, or frictional). Let us also mention a moel where both transverse an axial vibrations are investigate [26]. In the present work, we investigate a moel of transverse vibrations, erive from the moel propose in [1]. Let L > be the length of a axially symmetric thin beam, along the axis x. We also suppose that the microbeam is subjecte to heat an to a non uniform frictional amping. Then the transverse isplacement u an the temperature fiel θ satisfy the following set of equations (without imension), for all t >, for all x (,L): (1.1) (1.2) u tt +(p(x)u xx ) xx +2q(x)u t +κθ xx =, θ t ηθ xx κu xxt =, wherep is a non negativefunction escribingthe non homogeneityofthe beam, an q, also non negative stans for the force of the frictional amping. The constants κ, η are assume to be strictly positive. Note that system (1.1)-(1.2) is a simplifie version of the one investigate in [1]. Inee, although we ae inhomogeneity 1

2 an amping, the authors in [1] consier a more general law of heat conuction an assume the beam is moving along the x axis at a constant velocity, an that the temperature is also ampe. Moreover the coupling constants are not the same in their case (here we only have κ > ), so their problem is more ifficult to tackle. System (1.1)-(1.2) is complemente by Dirichlet bounary conitions (1.3) u(, t) = u(l, t) =, u x (, t) = u x (L, t) =, θ(, t) = θ(l, t) =, that hol for all t >, escribing a clampe beam at both ens. An finally we have initial conitions: (1.4) u(x, ) = u (x), u t (x, ) = u 1 (x), θ(x, ) = θ (x), x [, L]. Our purpose in this work is to provie well-poseness of system (1.1)-(1.4) uner smoothness assumption on the inhomogeneity functions p, q, an to investigate the long time behavior of the energy, that is to say wether or not the energy ecays exponentially. Inee, the stuy of the total energy functional E : R + R +, efine, in our case, for all t, by [ (1.5) E(t) = 1 L ] L p(x)u 2 xx x+ u 2 t x+ θ 2 x, 2 when t, answers the question of the stabilization of the vibrations of the structure. Rates of ecay for the total energy have been establishe for various issipative moels, see for example [2, 4, 5, 7, 8, 9, 1, 11] an references therein. Up to our knowlege, the exponential stability has not been investigate for our moel, an provies a first answer to the stabilization of microbeams uner stuy in [1], with of course only the Fourier law of heat, an a non moving structure in the x irection, for a start. Moreover, our result extens the results of [15, 16] for inhomogeneous beams without heat, an the ones of [18, 19] for ampe waves couple with heat. The rest of the paper is organize as follows. In Section 2 we state the mathematical hypothesis an prove well-poseness of the initial bounary value problem (1.1)-(1.4) using a semigroup argument. In Section 3 we prove the exponential ecay of the energy (1.5) which is the main result. Finally, in Section 4 we illustrate numerically the exponential ecay of the energy by consiering a simple Finite Difference iscretization, in the spirit of [22] an initial conitions close to a physically relevant situation, inspire by [23]. 2. Well-poseness In this section, we prove existence an uniqueness of solution to (1.1)-(1.4). Let us start by introucing the main assumptions an notations Preliminaries. Throughout this paper, we always assume that q(x) an p(x) are positive efinite functions, q(x) L (, L), p(x) H 2 (, L) an that there exist constants α 1, α 2, α 3 an α 4 such that (2.1) < α 1 p(x) α 2, x [, L] < α 3 q(x) α 4, x [, L] Throughout this paper C is a generic constant, not necessarily the same at each occasion(it will change from line to line), which epens in an increasing way on the inicate quantities. 2

3 We first establish an equality escribing the issipativity property of our system. Lemma 2.1. For every smooth enough solution of the system (1.1)-(1.4) the total energy E : R + R +, efine by (1.5) satisfies (2.2) E(t) = 2 t q(x) u t 2 x η θ x 2 x. Proof. We multiply (1.1), (1.2) by u t an θ respectively. Integrating in x [, L], using the bounary conition (1.3), summing an taking the real part, the result follows straightforwar. Remark 2.2. The equations (1.5)-(2.2) can be rewritten as E(t) = 1 ( ) pu xx 2 L 2 2 (,L) + u t 2 L 2 (,L) + θ 2 L 2 (,L) an satisfies t E(t) = 2 qu t 2 L 2 (,L) η θ x 2 L 2 (,L). Let us now introuce the phase space associate to our problem. We use the stanar L 2 (, L) space, the scalar prouct an norm are enote by u, v L 2 (,L) = uv x, u 2 L 2 (,L) = u 2 x. Hence, enote H = H 2(, L) L2 (, L) L 2 (, L) enowe with the inner prouct given by (u, v, θ), (u 1, v 1, θ 1 ) H = p(x)u xx u 1xx x+ for U = (u, v, θ), Ũ = (u 1, v 1, θ 1 ) an the norm vv 1 x+ (u, v, θ) 2 H = pu xx 2 L 2 (,L) + v 2 L 2 (,L) + θ 2 L 2 (,L). θθ 1 x We can easily show that the norm H is equivalent to the usual norm in H, so that the phase space H is a Hilbert space. Finally, we have the Poincaré inequality where C p is the Poincaré constant. u 2 L 2 (,L) C p u x 2 L 2 (,L), u H1 (, L), 2.2. Existence an uniqueness. In orer to establish the well-poseness, we reformulate our problem as an abstract Cauchy problem. Taking u t = v, the initial bounary value problem (1.1)-(1.4) can be rewritten as: (2.3) t U(t) = AU(t), U() = U, t >, where U(t) = (u, v, θ) T an U = (u, u 1, θ ) T. The linear operator A : D(A) H H is given by (2.4) A u v v = 2q(x)v (p(x)u xx ) xx κθ xx. θ κv xx +ηθ xx 3

4 It remains to efine the omain D(A) of the operator: D(A) = { (u, v, θ) H : u H 2 (, L), p(x)u xx κθ H 2 (, L), θ H 2 (, L), κv +ηθ H2 (, L)}. We are now reay to state our existence result. Theorem 2.3. For any U H there exists a unique solution U(t) = (u, u t, θ) of (1.1)-(1.4) (or, equivalently, to the abstract Cauchy problem (2.3)), satisfying However, if U D(A) then u C([, [; H 2 (, L)) C1 ([, [; L 2 (, L)) θ C([, [; L 2 (, L)). u C 2 ([, [; H 2 (Ω)) C2 ([, [; L 2 (, L)) θ C([, [; H 2 (, L)) C1 ([, [; L 2 (, L)) pu xx +κθ, κu t +ηθ C([, [; H 2 (, L)). Proof. Firstly, we show that the operator A generates a C -semigroup of contractions on the space H. We will show that A is a issipative operator an belongs to the resolvent set of A, enote by (A). Then our conclusion will follow using the well known Lumer-Phillips theorem [2]. We compute: AU, U H = + = 2iIm p(x)v xx u xx x [ 2q(x)v (p(x)u xx ) xx κθ xx ]v x+ + 2iκIm p(x)v xx u xx x 2 Taking real parts we obtain v x θ x x η q(x) v 2 x θ x 2 x. (2.5) Re AU, U H = 2 q(x)v 2 L 2 (,L) η θ x 2 L 2 (,L). (κv x +ηθ x ) x θ x Thus A is a issipative operator. Next, the omain D(A) is clearly ense in H an the operator is close. Finally, let us show that ρ(a), the resolvent of A. In fact, given F = (f, g, h) H, we must show that there exists a unique U = (u, v, θ) in D(A) such that AU = F, that is to say: (2.6) (2.7) (2.8) v = f in H 2 (, L) Taking v = f in (2.8) we have 2q(x)v (p(x)u xx ) xx κθ xx = g in L 2 (, L) κv xx +ηθ xx = h in L 2 (, L). (2.9) θ xx = 1 η (h κf xx) L 2 (, L). It is known that there is a unique θ H 2 (, L) satisfying (2.9) an θ xx L 2 (,L) 1 η h κf xx L 2 (,L) C F H 4

5 for a positive constant C. Moreover, we have that (κv+ηθ) xx = h L 2 (, L) an we conclue that (κv +ηθ) H 2 (, L). Replacing (2.6) into (2.7) we have (2.1) Moreover (p(x)u xx ) xx = g 2q(x)f κθ xx in L 2 (, L). u 2 xx x = 1 p(x) p(x)u2 xx x 1 p(x)u 2 xx α x = 1 p(x)u xx 2 L 1 α 2 (,L). 1 It follows using the Lax-Milgram theorem that there exists a unique function u H 2(, L)suchthat(2.1)holsinL2 (, L). Itiseasytoshowthat U H C F H for a positive constant C. Therefore we conclue that ρ(a). Therefore A is the generator of a C semigroup of contractions an the proof of Theorem 2.3 is complete. 3. Asymptotic Behavior In this section we show that the energy ecays uniformly exponentially with time. We state our main result as follows. Theorem 3.1. Let u, θ be solution to system (1.1)-(1.4). Then there exist positive constants M an m such that E(t) M E()e mt, t, where E is the total energy of the system, given by (1.5). We will prove this theorem using multiplier technique an a suitable Lyapunov. We begin with preliminaries Preliminaries. We recall the energy equality from Lemma 2.1, (3.1) t E(t) = 2 q(x)u 2 t x η θx 2 x. Let us now efine an auxiliary function that will help builing our Lyapunov: (3.2) F 1 (t) = uu t x+ q(x)u 2 x. We conclue this subsection by a lemma to control the time erivatives of the auxiliary function F 1. Lemma 3.2. Let u,θ be a smooth solution to (1.1)-(1.4). Then the auxiliary function F 1 efine by (3.2) satisfies the following estimates. There exists µ > such that for all t > : (3.3) µ E(t) F 1 (t) µe(t). Moreover, the erivative of F 1 verifies: (3.4) t F 1(t) = 2E(t)+R 1 (t). with the remainer R 1 such that: there exist non negative constants C 1, C 2, C 3 >, such that, for all α >, for any t >, it hols: ( (3.5) R 1 (t) C 1 u 2 t +C ) θx 2 +C 3 α pu 2 α xx. 5

6 Proof. We begin with the first estimate. From the efinition of F 1, (3.2), On the other han, using the classical Young an Poincaré inequalities successively, together with the assumption on functions p an q, (2.1) we obtain (3.6) F 1 (t) 1 2 u t 2 x+c(c p,α 1,α 4 ) p(x)u 2 xx x, where C(C p,α 1,α 4 ) is a non negative constant, epening on Poincaré constants, an the upper an lower bouns of p an p. The estimate (3.3) follows straightforwar with µ = 1 2 +C(C p,α 1,α 4 ). Next, we ifferentiate F 1 : by using (1.1) an integrating by parts we have t F 1(t) = = with = = 2E +R 1, uu tt x+ u 2 t x+ t q(x)u 2 x u( 2q(x)u t (p(x)u xx ) κθ xx )x+ p(x)u 2 xx x κ R 1 = 2 θu xx x+ u 2 t + θ 2 κ u 2 t x u 2 t x+ t θu xx. q(x)u 2 x Hence, using the Poincaré inequality, (2.1) an the generalize Young s inequality, we have, for all α >, for all t > : R 1 (t) 2 u 2 t +C p θx 2 + C p θx α p(x)u 2 2α 2α xx. 1 Estimate(3.5)followsirectly. Wewillchoseaparticularα > inthefollowing. With this Lemma, we are now reay to prove the exponential stability Proof of Theorem 3.1. We introuce a Lyapunov functional: (3.7) G(t) = E(t)+ε (F 1 (t)+f 2 (t)), where ε > is a positive constant to be ajuste later (as the α > of the previous Lemma). Let us compute the erivative of G. From (3.1), together with (3.4) an (3.5) from Lemma 3.2, we get: (3.8) t G(t) = 2 q(x)u2 t x η θ2 x x 2 εe(t)+εr 1(t) ( ) 2 ε C1 L α 3 q(x)u2 t x ( ( η εc L α)) θ2 x x ( ) L ε 2E(t) αc 3 pu2 xx. Now we will chose the constants ε an α. We nee to control the last term in (3.8) in orer to have a goo ifferential inequality. We also nee to ensure the positivity of the Lyapunov G. 6

7 On the one han, from the efinition of the total energy, (1.5), it hols: 2E(t) C 3 α pu 2 xx = (1 αc 3 ) Hence, choosing α > small enough so that we get, injecting in (3.8): t G(t) 2KεE(t) K 1 pu 2 xx + K := 1 αc 3 >, u 2 t x+ q(x)u 2 t x K 2 θx 2 x θ 2 x. where K > is efine above, while α is now fixe to ensure K >, an: ( K 1 = 2 ε C ) ( 1, K 2 = η εc ). α 3 α On the other han, from the efinition of G, (3.7), an from (3.3) in Lemma 3.2, we have, for all t > : (3.9) (1 ε µ) E(t) G(t) (1+ε µ) E(t). Therefore, it suffices now to choose ε >, small enough so that: 1 ε µ ( >, ) (3.1) K 1 = 2 ε C1 α 3 >, ( K 2 = η εc α) >. Hence, the positivity of the K i, i = 1,2 gives: (3.11) G(t) 2KεE(t), t where ε is now fixe, so that (3.1) is satisfie. Finally combining (3.9) with (3.11) integrate in time, we get, for all t > : with E(t) M 2Kε (1 εµ) The conclusion of Theorem 3.1, that is t E(s)s, M = G() 1 εµ = 1+εµ E(). 1 εµ E(t) M E()e mt, follows immeiately from the Gronwall Lemma, with: an the proof is complete. m = 2 K ε 1 εµ >, M = 1+εµ 1 εµ, 7

8 4. Numerical simulation In this section we provie a numerical illustration of Theorem 3.1, that is the stabilization of the vibrations thanks to the thermal effect. Few numerical simulations have been on for a EB type equation, couple with a amping an/or heat conuction. See for example [22], where the authors use a Finite Element scheme in space, but o no care about the time iscretization, since they are intereste in the Optimal Control Moel. Here we want to illustrate our result of exponential stability for physically relevant initial conitions. Therefore, since we o not aim at esigning a high orer scheme, we consier a classical secon orer Finite Difference scheme in both time an space (see for example [28]), taking care with the iscretization of the bilaplacian operator ( [25, 13, 27]). We consier the following initial conitions: { u(x,) =, x (,L), θ(x,) = Ae x2 σ 2, x (,L), with A =.2, whereas the stanar eviation σ 2 changes from 1 4 to 1 2. These initial conitions represent a possible physical experiment, on to test the response of a beam (the eformation) when we apply an initial heat pulse in the center of the beam, see [23]. We plot in Figure 4 the iscrete energy (compute with iscrete integration at all iscrete times t n, of the continuous expression (1.5)), versus time in Log Log scale, an observe the exponential ecay for three values of the stanar eviation. Conclusion In this work, we show the exponential stabilization of an EB type moel of flexible inhomogeneous beam, taking into account frictional amping an heat conuction of Fourier type. We illustrate our result numerically, consiering physically relevant initial conitions, an our result generalizes the results obtaine for the wave equation an for EU moels without heat. We think that this result is a first step towar a better unerstaning of more complicate moels, such as the one erive in [1]. 8

9 Acknowlegements Octavio Vera thanks the support of the Fonecyt project We thank Centro CRHIAM Project Conicyt/Fonap an Re Doctoral REDOC.CTA. Amelie Rambau thanks the support of Fonecyt project References [1] A.E. Abouelregal, an A.M. Zenkour. Thermoelastic problem of an axially moving microbeam subjecte to an external transverse excitation. Journal of Theoretical an Applie Mechanics. vol. 53 (215) pp [2] K., Ammari, M., Dimassi, M., Zerzeri. The rate at which energy ecays in a viscously ampe hinge Euler Bernoulli beam. Journal of Differential equations. Vol 257, (9) (214) [3] K. Ammari an M. Tuesnak. Stabilization of Bernoulli-Euler beams by means of a point feeback force. SIAM Journal of Control an Optimization. Vol (2) [4] G. Chen. Energy ecay estimates an exact bounary-value controllability for the wave equation in a boune omain. Journal e Mathématiques Pures et Appliquées. Vol. 58.(1979) [5] G. Chen. A note on the bounary stabilization of the wave equation. SIAM J. Control Opt. 19(1981) [6] A. Dimarogonas. Vibration for Engineers. 2n e. Prentice-Hall. Inc. (1996). [7] V. Komornik. Rapi bounary stabilization of wave equations. SIAM J. Control an Optimization. 29(1991) [8] V. Komornik, E. Zuazua. A irect metho for bounary stabilization of the wave equation. J. Math. Pures Appl. 69(199) [9] J. Lagnese. Note on bounary stabilization of wave equations. SIAM j. Control Opt. 26(1988) [1] J. Lagnese. Decay of solutions of wave equations in a boune region with bounary issipation. J. Diff. Eqns., 5(1983) [11] J.L. Lions. Exact controllability, stabilization an perturbation for istribute system. SIAM review. Vol. 3. 1(1988)1-68. [12] K. Liu an Z. Liu. Exponential ecay of energy of the Euler-Bernoulli beam with locally istribute Kelvin-Voigt amping. SIAM Journal of Control an Optimization. Vol (1998) [13] Xu-zheng Liu, Xia Cui, Jia-guang Sun. FDM for multi-imensional nonlinear couple system of parabolic an hyperbolic equations. Journal of Computational an Applie Mathemathics, 186 (26), [14] L. Meirovitch. Funamentals of Vibrations. International Eition. McGraw-Hill. (21). [15] P.K. Nani, G. C. Gorain an S. Kar. A note on stability of longituinal vibrations of an inhomogeneous beam. Applie Mathematics. 3(212) [16] G.C. Gorain. Exponential Stabilization of Longituinal Vibrations of an Inhomogeneous Beam. Journal of mathematical sciences.198 (214) [17] S.M. Han, an H. Benaroya, an T. Wei. Dynamics of transversely vibrating beams using four engineering theories. Journal of Soun an Vibrations. 225 (1999) [18] S. Misra, an R. Ghosh, an G. C. Gorain. A note on stability of a ampe wave couple with thermal effect. Preprint (215) [19] S. Misra, an R. Ghosh, an G. C. Gorain Stabilization of transverse vibrations of an inhomogeneous beam. Preprint. [2] A. Pazy. Semigroup of linear operators an applications to partial ifferential equations. Springer-Verlag (1983). New york. [21] W. T. Thomson. Theory of Vibration with Applications. 2n e. (1981). [22] Yu, Xinan Ren, Zhigang an Zhang, Qian an Xu, Chao. A Numerical Metho of the Euler- Bernoulli Beam with Optimal Local Kelvin-Voigt Damping. Journal of Applie Mathematics. (214) [23] X. Wang, X. Xu. Thermoelastic wave inuce by pulse laser heating. Appl. Phys. A. 73, (21). [24] W. Weaver, S.P. Timoshenko, D.H. Young. Vibration Problems in Engineering. 5th e. John Wiley an Sons. Inc. (199). 9

10 [25] L.V. Kantorovitch, V.I. Krylov. Approximation Methos of Higher Analysis, Interscience, New-York, [26] E. Kaya-Cekin, E. Aulisia, A. Ibragimov. Stability analysis of inhomogeneous equilibrium for axially an transversely excite nonlinear beam. Communications on Pure an Applie Analysis, vol 1 (5). (211) [27] L.F. Richarson. The approximate arithmetical solution by finite ifference of physical problems involving ifferential equations with an application to the stress in Mansory am. Transaction of the Royal Society of Lonon, Ser.A, 21, (191). [28] G. D. Smith, Numerical Solution of Partial Differential Equations: Finite Difference Methos. Clarenon Press. (1985) σ 2 =1-2 σ 2 =1-3 σ 2 = E t Figure 1. Exponential ecay Department of Mathematics, Universia el Bío-Bío, PO Box 5C, Concepción, Chile aress: overa@ubiobio.cl Department of Mathematics, Universia el Bío-Bío, PO Box 5C, Concepción, Chile aress: arambau@ubiobio.cl Department of Physics, Universia el Bío-Bío, PO Box 5C, Concepción, Chile aress: rrozas@ubiobio.cl 1