SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

This article is part of the series Recent Advances in Operator Equations, Boundary Value Problems, Fixed Point Theory and Applications, and General Inequalities.

Open Access Research

Exponential and polynomial stabilization of the Kirchhoff string by nonlinear boundary control

Yuhu Wu1* and Jianmin Wang2

Author Affiliations

1 Department of Mathematics, Harbin University of Science and Technology, Harbin, 150080, P.R. China

2 Department of Electronics Science and Technology, Harbin University of Science and Technology, Harbin, 150080, P.R. China

For all author emails, please log on.

Boundary Value Problems 2013, 2013:215  doi:10.1186/1687-2770-2013-215

The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2013/1/215


Received:22 April 2013
Accepted:6 August 2013
Published:21 October 2013

© 2013 Wu and Wang; licensee Springer

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

This paper addresses the stabilization problem of the nonlinear Kirchhoff string using nonlinear boundary control. Nonlinear boundary control is the negative feedback of the transverse velocity of the string at one end, which satisfies a polynomial-type constraint. Employing the multiplier method, we establish explicit exponential and polynomial stability for the Kirchhoff string. The theoretical results are assured by numerical results of the asymptotic behavior for the system.

1 Introduction

Stabilization and vibration controllability of string or beam systems arising from different engineering backgrounds has attracted attention of many researchers [1-4]. In particular, boundary feedback stabilization of string and beam systems has become an important research area [5-7]. This is because, in a practice system, vibration is more easily controlled through a boundary point than using point sensors or actuators away from the boundaries [8,9].

There are several nonlinear mathematical models that describe the transversal vibration of stretched strings. One such model is presented in the following equation:

y t t ( x , t ) = [ a + b 0 1 y x 2 ( x , t ) d x ] y x x ( x , t ) (1)

for all x ( 0 , 1 ) and t 0 , where a > 0 , b 0 are two constants. Obviously, the above equation is a simple prototype of the classical equation

ρ h y t t ( x , t ) = [ T 0 + h E 2 l 0 l y x 2 ( x , t ) d x ] y x x ( x , t ) ,

which was proposed by Kirchhoff [10]. Here l is the length of the string; E is Young’s modulus of the material; ρ is density; h is the area of the cross section; y ( x , t ) is the transversal displacement of the point x of the string at time t. This model has been studied by researchers from the physical and mathematical points of view; see, e.g., references [11-13] and the references therein.

In this paper, we consider Kirchhoff string (1) with the following boundary conditions (see Figure 1):

y ( 0 , t ) = 0 , (2)

T ( t ) y x ( 1 , t ) = u ( t ) (3)

for all t 0 . T ( t ) denotes the tension in the string at time t. The boundary condition in equation (2) implies that the string is fixed at x = 0 . The boundary condition in equation (3) represents the balance of the transversal component of the tension in the string and the control input u which is applied transversally at x = 1 . Because the tension in the string represented by equation (1) is not constant and is given by

T ( t ) = a + b 0 1 y x 2 ( x , t ) d x (4)

for all t 0 (see [14]), the boundary condition in equation (3) can be written as

[ a + b 0 1 y x 2 ( x , t ) d x ] y x ( 1 , t ) = u ( t ) . (5)

thumbnailFigure 1. Schematic of the nonlinear Kirchhoff string with boundary control.

Shahruz and Krishna [13] investigated the stabilization of Kirchhoff string (1) with a linear negative velocity control, which means the boundary control u has a linear negative velocity feedback form

u ( t ) : = u ( y t ( 1 , t ) ) = L y t ( 1 , t ) (6)

for all t 0 , where L is a positive constant. They established exponential stability. In [15], the absolute stability of the Kirchhoff string (1) with linear sector boundary control was considered. It is well known that linear strings represented by equation (1), for which b = 0 , can be stabilized by the control law in equation (6); see e.g., references [16-18]. Moreover, Shahruz [19], Fung et al.[5], and Li and Hou [20] developed linear boundary control laws for axially moving strings. It is worth mentioning that Kobayashi [21] designed a linear parallel compensator based on boundary displacement observer and proved the string (1) can be stabilized by parallel compensator control.

In the literature mentioned above, such as [13,16] and [20], the exponential stabilization result for various string systems by linear boundary control mainly relies on the Lyapunov direct method. In this work, we investigate the stabilization of string (1) with a more general and ‘flexible’ boundary control (see hypothesis (H) in Section 2). The feedback function u is not required to satisfy a strict control law such as (6), but just satisfies some appropriate polynomial-type constraint. In this general boundary control case, it seems that the Lyapunov direct method is no more applicable. So, we need to use a more meticulous method to deal with the stabilization problem. Applying the multiplier method, we establish not only exponential stability result but also polynomial stability result for Kirchhoff string (1).

The remainder of this technical paper is arranged as follows. Section 2 describes the model of the Kirchhoff nonlinear string and introduces the control assumption. The problem of exponential and polynomial stability is addressed in Section 3. Finally, a numerical example is demonstrated where the nonlinear distributed parameter infinite-dimensional equation is solved by applying the finite element method in Section 4.

2 Problem formulation

Consider the nonlinear Kirchhoff string model as shown in Figure 1. For the sake of easy reading and later referring, the governing equation, the boundary conditions and the initial functions are put together as

for all x ( 0 , 1 ) and t 0 . Here f ( x ) and g ( x ) in equation (7d) are the initial displacement and velocity of the string, respectively. We assume that f C 1 [ 0 , 1 ] , g C [ 0 , 1 ] and that at least one of the functions f or g is not identically zero over [ 0 , 1 ] .

To obtain a precise stabilization result, we make the following hypothesis on the continuous control feedback u : R R :

(H) There exist constants L 2 L 1 > 0 and r 1 such that

u ( w ) Ω ¯ ( w ) : = { Co { L 1 w , L 2 w } , | w | > 1 , Co { L 1 w r , L 2 w 1 r } , | w | 1 . (8)

Remark 2.1 Obviously, condition (H) is equivalent to, for all w R ,

{ u ( w ) w 0 , L 1 min { | w | , | w | r } | u ( w ) | L 2 max { | w | , | w | 1 r } .

Obviously, hypothesis (H) is a ‘flexible’ and ‘robust’ condition, which allows the feedback function u to vary in an appropriate geometric region given by a polynomial-type constraint. For example, Figure 2 illustrates a feedback control u satisfying a linear sector constraint (H) with r = 1 , L 1 = 1 , L 2 = 2 , and Figure 3 shows a feedback control u satisfying a nonlinear constraint (H) with r = 2 , L 1 = 1 , L 2 = 2 . Since u ( y t ( 1 , t ) ) y t ( 1 , t ) 0 , for all t 0 , the boundary control (7c) is the negative feedback of transversal velocity of the string at x = 1 .

thumbnailFigure 2. Linear sector constraint.

thumbnailFigure 3. Nonlinear sector constraint.

For the existence and uniqueness of the solution of the general Kirchhoff equation, we refer to [11,12] and references therein. In this work, we study the stabilization of the string in (7a) by this negative feedback boundary control u, which provides a dissipative effect.

Remark 2.2 According to boundary condition (7b) at x = 0 , we easily get

y t ( 0 , t ) = 0 for all  t 0 . (9)

We define the natural energy function of time for system (7a)-(7d) and denote it by t E ( t ) . The scalar-valued function E is defined as

E ( t ) : = 1 2 0 1 y t 2 ( x , t ) d x + 1 2 0 1 a y x 2 ( x , t ) d x + b 4 [ 0 1 y x 2 ( x , t ) d x ] 2 . (10)

Especially, from the initial displacement and velocity condition (7d), we obtain the initial energy as

E ( 0 ) = 1 2 0 1 [ g 2 ( x ) + a f x 2 ( x ) ] d x + b 4 [ 0 1 f x 2 ( x ) d x ] 2 .

Since at least one of the functions f and g is not identically zero over ( 0 , 1 ) we have E ( 0 ) > 0 .

3 Stabilization by boundary control

In this section we state and prove our main result. For this purpose we establish several lemmas.

Lemma 3.1Let y ( , ) be the solution for system (7a)-(7d). Then

0 1 y x y x t d x = y x ( 1 , t ) y t ( 1 , t ) 0 1 y x x y t d x , (11)

0 1 x y x y x x d x = 1 2 y x 2 ( 1 , t ) 1 2 0 1 y x 2 ( x , t ) d x , (12)

0 1 x y x y t t d x = 0 1 ( [ x y x y t ] t + 1 2 y t 2 ) d x 1 2 y t 2 ( 1 , t ) . (13)

Proof See the Appendix. □

Now, we give a property of the energy function E.

Proposition 3.1The time-derivative of the energy functionEin equation (10), along the solution of system (7a)-(7d) satisfies

E ( t ) = u ( t ) y t ( 1 , t ) 0 (14)

for all t 0 .

Proof Differentiating the energy function (10) with respect to t, we get

E ( t ) = 0 1 y t y t t d x + [ a + b 0 1 y x 2 d x ] 0 1 y x y x t d x . (15)

According to equation (11) and boundary control (7c), we get

[ a + b 0 1 y x 2 d x ] 0 1 y x y x t d x = u ( t ) y t ( 1 , t ) [ a + b 0 1 y x 2 d x ] 0 1 y x x y t d x . (16)

Substituting equation (16) into equation (15) and observing (7a), we obtain

E ( t ) = u ( t ) y t ( 1 , t ) + 0 1 y t t y t d x [ a + b 0 1 y x 2 d x ] 0 1 y x x y t d x = u ( t ) y t ( 1 , t ) (17)

for all t 0 . We obtain equation (14). □

Remark 3.1 From Proposition 3.1, we obtain the energy identity for system (7a)-(7d),

E ( S ) E ( T ) = T S u ( t ) y t ( 1 , t ) d t .

Therefore, the energy E is a decreasing function of time.

During the subsequent stability analysis, we utilize the following inequality.

Lemma 3.2Let y ( , ) be the solution for system (7a)-(7d). Then

0 1 x y x ( x , t ) y t ( x , t ) d x max { 1 , 1 a } E ( t ) (18)

for all t 0 .

Proof Applying the Cauchy-Schwarz inequality, we get

0 1 x y x ( x , t ) y t ( x , t ) d x 1 2 0 1 x 2 y x 2 d x + 1 2 0 1 y t 2 d x 1 2 0 1 y x 2 d x + 1 2 0 1 y t 2 d x (19)

for all t 0 . On the other hand, the definition of energy function (10) implies

a 2 0 1 y x 2 d x + 1 2 0 1 y t 2 d x E ( t )

for all t 0 . It follows from the above inequality that

1 2 0 1 y x 2 d x + 1 2 0 1 y t 2 d x { E ( t ) , a 1 , 1 a E ( t ) , 0 < a < 1 . (20)

Together with (19) and (20), we get equation (18). Hence we complete the proof of Lemma 3.2. □

Now, we present a Gronwall-type lemma (see Komornik [22], pp.124), which will play an essential role when establishing the stabilization result.

Lemma 3.3Let G : R + R + be a non-increasing function. Assume that there exists a constant ω > 0 such that

T + G q + 1 ( t ) d t 1 ω G ( T ) for all  T 0 .

Then the following estimation is true, for all t 0 ,

{ G ( t ) G ( 0 ) e 1 ω t if  q = 0 , G ( t ) G ( 0 ) ( 1 + q q ω t ) 1 q if  q > 0 .

We give a priori estimation for the energy function E ( t ) , which was established in [15]. For the sake of completeness, we give the proof here.

Lemma 3.4The energy functionEin equation (10), along the solution of system (7a)-(7d), satisfies

E ( t ) 0 1 [ x y x y t ] t d x + 1 2 y t 2 ( 1 , t ) + 1 2 a u 2 ( y t ( 1 , t ) ) (21)

for all t 0 .

Proof We multiply equation (7a) by x y x ( x , t ) and do integration over [ 0 , 1 ] , with respect to x. We obtain

0 = 0 1 x y x ( y t t [ a + b 0 1 y x 2 d x ] y x x ) d x = 0 1 ( [ x y x y t ] t + 1 2 y t 2 ) d x 1 2 y t 2 ( 1 , t ) 1 2 [ a + b 0 1 y x 2 d x ] ( y x 2 ( 1 , t ) 0 1 y x 2 d x ) , (22)

using equations (13) and (12) in Lemma 3.1. It follows from (10) that

1 2 0 1 y t 2 d x + 1 2 [ a + b 0 1 y x 2 d x ] 0 1 y x 2 d x = E ( t ) + b 4 [ 0 1 y x 2 d x ] 2 . (23)

Since a + b 0 1 y x 2 d x a , according to boundary control (7c), we have

[ a + b 0 1 y x 2 d x ] y x 2 ( 1 , t ) 1 a u 2 ( y t ( 1 , t ) ) . (24)

Hence, substituting equation (23) into equation (22) and using equation (24), one has

E ( t ) + b 4 [ 0 1 y x 2 d x ] 2 = 0 1 [ x y x y t ] t d x + 1 2 y t 2 ( 1 , t ) + 1 2 [ a + b 0 1 y x 2 d x ] y x 2 ( 1 , t ) 0 1 [ x y x y t ] t d x + 1 2 y t 2 ( 1 , t ) + 1 2 a u 2 ( y t ( 1 , t ) ) .

Since b 4 [ 0 1 y x 2 d x ] 2 0 , we complete the proof of Lemma 3.4. □

Lemma 3.5For any constant q 0 , the energy functionEalong the solution of system (7a)-(7d) satisfies the following estimation, for all S > T 0 ,

T S E q + 1 ( t ) d t C 1 E q + 1 ( T ) + 1 2 T S E q ( t ) y t 2 ( 1 , t ) d t + 1 2 a T S E q ( t ) u 2 ( y t ( 1 , t ) ) d t , (25)

where C 1 = 3 q + 2 q + 1 max { 1 , 1 a } .

Proof According to inequality (21) in Lemma 3.4, we have, for all 0 T < S ,

T S E q + 1 ( t ) d t T S E q ( t ) 0 1 [ x y x y t ] t d x d t + 1 2 T S E q ( t ) y t 2 ( 1 , t ) d t + 1 2 a T S E q ( t ) u 2 ( y t ( 1 , t ) ) d t . (26)

Moreover, using integration by parts, we get

T S E q ( t ) 0 1 [ x y x y t ] t d x d t = T S [ E q ( t ) 0 1 x y x y t d x ] t d t T S q E ( t ) E q 1 ( t ) 0 1 ( x y x y t ) d x d t .

Hence, inequality (26) becomes

T S E q + 1 ( t ) d t A 1 + A 2 + 1 2 T S E q ( t ) y t 2 ( 1 , t ) d t + 1 2 a T S E q ( t ) u 2 ( y t ( 1 , t ) ) d t , (27)

where

A 1 : = T S [ E q ( t ) 0 1 x y x y t d x ] t d t , A 2 : = T S q E ( t ) E q 1 ( t ) 0 1 ( x y x y t ) d x d t .

Firstly, we estimate A 1 and A 2 ,

A 1 = T S [ E q ( t ) 0 1 x y x y t d x ] t d t | E q ( T ) 0 1 x y x ( x , T ) y t ( x , T ) d x | + | E q ( S ) 0 1 x y x ( x , S ) y t ( x , S ) d x | ( by (18) ) max { 1 , 1 a } E q + 1 ( T ) + max { 1 , 1 a } E q + 1 ( S ) 2 max { 1 , 1 a } E q + 1 ( T ) , (28)

where the last inequality follows from the fact E ( t ) is a decreasing function. On the other hand,

A 2 = T S q E ( t ) E q 1 ( t ) 0 1 x y x y t d x d t T S | q E ( t ) E q 1 ( t ) | | 0 1 x y x y t d x | d t max { 1 , 1 a } | T S q E ( t ) E q ( t ) d t | ( by (18) ) = q q + 1 max { 1 , 1 a } | E q + 1 ( S ) E q + 1 ( T ) | q q + 1 max { 1 , 1 a } E q + 1 ( T ) . (29)

Finally, inserting the two inequalities, (28) and (29), in (27), we get inequality (25). This completes the proof of Lemma 3.5. □

We now state the main stabilization result for system (7a)-(7d).

Theorem 3.1Assume that assumption (H) holds. Then there exist three constants k 1 , k 2 , σ > 0 such that, for all t 0 ,

{ E ( t ) k 1 e σ t if  r = 1 , E ( t ) k 2 t 2 r 1 if  r > 1 . (30)

Remark 3.2 It is worth to mention that Theorem 2.4 in [13] can be viewed as a special cases of Theorem 3.1. Indeed, in the linear control case (6), the exponential stability in Theorem 3.1 coincides with the result in [13].

Proof of Theorem 3.1 We distinguish two cases related to the parameter r to establish the energy decay rate.

Case (I): r = 1 ;

Case (II): r > 1 .

In Case (I), we choose q = 0 . According to hypothesis (H), we know that

s 2 1 L 1 u ( s ) s , u 2 ( s ) L 2 u ( s ) s for all  s R .

Hence, from inequality (25) and equation (14), we deduce that, for all S > T 0 ,

T S E ( t ) d t C 1 E ( T ) + 1 2 T S [ y t 2 ( 1 , t ) + 1 a u 2 ( y t ( 1 , t ) ) ] d t C 1 E ( T ) a + L 1 L 2 2 a L 1 T S u ( t ) y t ( 1 , t ) d t = C 1 E ( T ) + a + L 1 L 2 2 a L 1 T S E ( t ) d t C 2 E ( T ) , (31)

where C 2 = C 1 + 1 2 ( 1 L 1 + L 2 a ) and C 1 is given in Lemma 3.5.

Now we deal with Case (II). In this case, we choose q = r 1 2 > 0 . We first admit the following fact (the proof is given in the Appendix).

Claim 1For any δ > 0 , we have the following estimates, for all S > T 0 ,

T S E q ( t ) y t 2 ( 1 , t ) d t r 1 r + 1 δ r + 1 r 1 T S E q r + 1 r 1 ( t ) d t T S E q ( t ) y t 2 ( 1 , t ) d t + 2 δ r + 1 2 ( r + 1 ) L 1 E ( T ) + 1 ( q + 1 ) L 1 E q + 1 ( T ) , (32)

T S E q ( t ) u 2 ( y t ( 1 , t ) ) d t r 1 r + 1 δ r + 1 r 1 T S E q r + 1 r 1 ( t ) d t T S E q ( t ) u 2 ( y t ( 1 , t ) ) d t + 2 L 2 r δ r + 1 2 ( r + 1 ) E ( T ) + L 2 ( q + 1 ) E q + 1 ( T ) . (33)

Now, inserting inequalities (32) and (33) into (25), we obtain, for all S > T 0 ,

T S E q + 1 ( t ) d t ( 1 2 + 1 2 a ) r 1 r + 1 δ r + 1 r 1 T S E q r + 1 r 1 ( t ) d t + C 3 E q + 1 ( T ) + C 4 E ( T ) , (34)

where C 3 = C 1 + a + L 1 L 2 2 a ( q + 1 ) L 1 , C 4 = 1 ( r + 1 ) L 1 δ r + 1 2 + L 2 r ( r + 1 ) a δ r + 1 2 and C 1 is given in Lemma 3.5. Now we choose δ = [ ( 1 + 1 a ) r 1 r + 1 ] r 1 r + 1 . Then it is obvious that ( 1 + 1 a ) r 1 r + 1 δ r + 1 r 1 = 1 . Hence, inequality (34) becomes

T S E q + 1 ( t ) d t 1 2 T S E q r + 1 r 1 ( t ) d t + C 3 E q + 1 ( T ) + C 4 E ( T ) .

Recalling q = r 1 2 , we get q + 1 = q r + 1 r 1 = r + 1 2 . Hence, the above inequality is rewritten as

T S E r + 1 2 ( t ) d t 2 C 3 E r + 1 2 ( T ) + 2 C 4 E ( T ) 2 ( C 3 E r 1 2 ( 0 ) + C 4 ) E ( T ) , (35)

where the last inequality follows from Remark 3.1. Finally, by letting S + in (31), (35) and using Lemma 3.3 with G ( t ) = E ( t ) , we complete the proof of Theorem 3.1.  □

Remark 3.3 According to the proof of Theorem 3.1, it is easy to see that the constants σ, k 1 and k 2 in Theorem 3.1 can be chosen as, respectively, σ 1 = 2 max { 1 , 1 a } + 1 2 ( 1 L 1 + L 2 a ) , k 1 = e E ( 0 ) , and k 2 = E ( 0 ) [ 2 r + 1 r 1 ( C 3 E r 1 2 ( 0 ) + C 4 ) ] 2 r 1 with C 3 = 3 r 1 r + 1 max { 1 , 1 a } + a + L 1 L 2 a ( r + 1 ) L 1 , C 4 = 1 r + 1 ( 1 L 1 + L 2 r a ) [ ( 1 + 1 a ) r 1 r + 1 ] r 1 2 . This means that the coefficients of the exponential or polynomial decay rate are exactly determined only by the initial tension a, the initial energy E ( 0 ) and the feedback control u. However, in the polynomial decay case, the order of decay rate is determined only by the feedback control u.

Finally, it is shown that the boundary control u stabilizes the nonlinear Kirchhoff string.

Theorem 3.2Assume that assumption (H) holds. Then there exist two constants k 1 , k 2 > 0 such that for all x ( 0 , 1 ) and t 0 ,

{ | y ( x , t ) | k 1 e σ t 2 if  r = 1 , | y ( x , t ) | k 2 t 1 r 1 if  r > 1 , (36)

where k 1 = 2 k 1 a , k 2 = 2 k 2 a and k 1 , k 2 , σare given in Theorem 3.1.

Proof According to the fact that y x ( 0 , t ) = 0 , for all t 0 , we get

| y ( x , t ) | = | 0 x y x ( z , t ) d z | 0 1 | y x ( z , t ) | d z ( 0 1 y x 2 ( x , t ) d x ) 1 2 ( by the Cauchy inequality ) 2 E ( t ) a ( by equation (14) ) (37)

for all x ( 0 , 1 ) and t 0 . By combining (37) with (30) in Theorem 3.1, we complete the proof of Theorem 3.2. □

4 Numerical results

In this section we consider a computational example for the closed-loop system (7a)-(7d). To illustrate the control performance of the boundary control law satisfying condition (H), numerical simulations by using the finite element method (FEM) are performed. We use Lagrange ‘hat’ basis with FEM equidistant meshes. The system parameters used in the simulations are a = 2 , b = 3 . The initial conditions are f ( x ) = 0.6 sin ( 5 x ) and g ( x ) = 0.5 sin ( 3 x ) . That is we consider the following Kirchhoff system:

{ y t t ( x , t ) = [ 2 + 3 0 1 y x 2 ( x , t ) d x ] y x x ( x , t ) , y ( 0 , t ) = 0 , [ a + b 0 1 y x 2 ( x , t ) d x ] y x ( 1 , t ) = u ( y t ( 1 , t ) ) , y ( x , 0 ) = 0.6 sin ( 5 x ) , y t ( x , 0 ) = 0.5 sin ( 3 x ) . (38)

The dynamic responses of the controlled Kirchhoff string were simulated under two feedback control laws:

u 1 ( x ) = 2 x sin x , x R

and

u 2 ( x ) = { | x | , | x | < 1 , x , | x | 1 .

Obviously, the feedback control function u 1 satisfies the constraint condition (H) with r = 1 , L 1 = 1 , L 2 = 3 , and u 2 satisfies the constraint condition (H) with r = 2 , L 1 = L 2 = 1 . Then, according to Theorem 3.2, the asymptotic behavior of the transverse vibration of system (38) under control law u 1 (or control law u 2 ) possesses exponential decay (or polynomial decay with degree −1, because r = 2 ). The string response y ( x , t ) of closed-loop system (38) with control law u 1 and control law u 2 are shown in Figure 4 and Figure 5, respectively. The corresponding transversal displacement at x = 0.5 is shown in Figures 6 and 7, respectively.

thumbnailFigure 4. String response of closed-loop system (38) under control law u 1 .

thumbnailFigure 5. String response of closed-loop system (38) under control law u 2 .

thumbnailFigure 6. Transverse displacement of the string at x = 0.5 under control law u 1 .

thumbnailFigure 7. Transverse displacement of the string at x = 0.5 under control law u 2 .

From Figures 4 and 5, it can bee seen that, in the case of control law u 2 , the decay of the transverse vibration relatively slow compared to the case of control law u 1 . Indeed, from Figure 6 and Figure 7, we know that in the case of control law u 1 , the transverse vibration has been suppressed exponentially ( | y ( 0.5 , t ) | K 2 e σ t ), whereas in the case of control law u 2 , the transverse vibration has been suppressed polynomially ( | y ( 0.5 , t ) | K 3 t 1 ). It is coincident with the results of Theorem 3.2.

Appendix

A.1 Proof of Lemma 3.1

Integrating by parts and (9) we get, for all t 0 ,

0 1 y x ( x , t ) y x t ( x , t ) d x = 0 1 [ y x ( x , t ) y t ( x , t ) ] x d x 0 1 y x x ( x , t ) y t ( x , t ) d x = y x ( 1 , t ) y t ( 1 , t ) 0 1 y x x ( x , t ) y t ( x , t ) d x .

So we obtain equation (11). Next, integrating by parts we compute

0 1 x y x y x x d x = 1 2 0 1 [ x y x 2 ] x d x 1 2 0 1 y x 2 d x = 1 2 y x 2 ( 1 , t ) 1 2 0 1 y x 2 d x

for all t 0 . That is, equation (12) holds. Finally, we write

0 1 x y x y t t d x = 0 1 [ x y x y t ] t d x 0 1 x y t y x t d x (39)

and

0 1 x y t y x t d x = 1 2 0 1 [ x y t 2 ] x d x 1 2 0 1 y t 2 d x (40)

for all t 0 . Substituting equation (40) into equation (39), we obtain equation (13). Thus the proof of Lemma 3.1 is complete.

A.2 Proof of Claim 1

Since r > 1 , hypothesis (H) implies

L 1 | s | r | u ( s ) | L 2 | s | 1 r if  | s | 1 , L 1 | s | | u ( s ) | L 2 | s | if  | s | > 1 .

Hence, it is true that, for all s R ,

s 2 ( 1 L 1 u ( s ) s ) 2 r + 1 1 L 1 u ( s ) s , (41)

u 2 ( s ) ( L 2 r u ( s ) s ) 2 r + 1 L 2 u ( s ) s . (42)

Then, from (41) and Proposition 3.1 we have

T S E q ( t ) y t 2 ( 1 , t ) d t T S E q ( t ) ( 1 L 1 u ( t ) y t ( 1 , t ) ) 2 r + 1 d t + T S E q ( t ) ( 1 L 1 u ( t ) y t ( 1 , t ) ) d t = T S E q ( t ) ( E ( t ) L 1 ) 2 r + 1 d t 1 L 1 T S E q ( t ) E ( t ) d t T S E q ( t ) ( E ( t ) L 1 ) 2 r + 1 d t + E q + 1 ( T ) ( q + 1 ) L 1 . (43)

From Young’s inequality, we have, for any α , β , δ > 0 ,

α β δ ρ α ρ ρ + δ ρ β ρ ρ , (44)

where ρ , ρ > 1 and 1 ρ + 1 ρ = 1 . Applying the above inequality with α = E q ( t ) , β = ( E ( t ) L 1 ) 2 r + 1 , ρ = r + 1 r 1 and ρ = r + 1 2 , we obtain, for any δ > 0 ,

T S E ( t ) q ( E ( t ) L 1 ) 2 r + 1 d t r 1 r + 1 δ r + 1 r 1 T S E q r + 1 r 1 ( t ) d t + 2 δ r + 1 2 ( r + 1 ) L 1 T S E ( t ) d t r 1 r + 1 δ r + 1 r 1 T S E q r + 1 r 1 ( t ) d t + 2 δ r + 1 2 ( r + 1 ) L 1 E ( T ) . (45)

By inserting (45) into (43), we deduce inequality (32) in Claim 1. Similarly, from (42) and Proposition 3.1 we have

T S E q ( t ) u 2 ( y t ( 1 , t ) ) d t T S E q ( t ) ( L 2 r u ( t ) y t ( 1 , t ) ) 2 r + 1 d t + L 2 T S E q ( t ) ( u ( t ) y t ( 1 , t ) ) d t = T S E q ( t ) ( L 2 r E ( t ) ) 2 r + 1 d t L 2 T S E q ( t ) E ( t ) d t T S E q ( t ) ( L 2 r E ( t ) ) 2 r + 1 d t + L 2 ( q + 1 ) E q + 1 ( T ) . (46)

Applying again inequality (44) with α = E q ( t ) , β = [ L 2 r E ( t ) ] 2 r + 1 , ρ = r + 1 r 1 and ρ = r + 1 2 , we obtain, for any δ > 0

T S E ( t ) q ( L 2 r E ( t ) ) 2 r + 1 d t r 1 r + 1 δ r + 1 r 1 T S E q r + 1 r 1 ( t ) d t + 2 L 2 r δ r + 1 2 ( r + 1 ) T S E ( t ) d t r 1 r + 1 δ r + 1 r 1 T S E q r + 1 r 1 ( t ) d t + 2 L 2 r δ r + 1 2 ( r + 1 ) E ( T ) . (47)

Inserting (47) into (46), we get inequality (33) in Claim 1.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors contributed equally to the manuscript and read and approved the final manuscript.

Acknowledgements

This work is supported by the National Natural Science Foundation of China (Grant 11226128) and the Natural Science Foundation of Heilongjiang Province of China (Grant F201113).

References

  1. Chen, G, Delfour, MC, Krall, AM, Payre, G: Modeling, stabilization and control of serially connected beams. SIAM J. Control Optim.. 25, 526–546 (1987). Publisher Full Text OpenURL

  2. Lagnese, J: Decay of solutions of wave equations in a bounded domain with boundary dissipation. J. Differ. Equ.. 50, 163–182 (1983). Publisher Full Text OpenURL

  3. Shahruz, SM: Bounded-input bounded-output stability of a damped nonlinear string. IEEE Trans. Autom. Control. 41, 1179–1182 (1996). Publisher Full Text OpenURL

  4. Guo, BZ: Riesz basis approach to the stabilization of a flexible beam with a tip mass. SIAM J. Control Optim.. 39, 1736–1747 (2001). Publisher Full Text OpenURL

  5. Fung, RF, Wu, JW, Wu, SL: Exponential stabilization of an axially moving string by linear boundary feedback. Automatica. 35, 177–181 (1999). Publisher Full Text OpenURL

  6. Morgül, Ö: Stabilization disturbance rejection for the beam equation. IEEE Trans. Autom. Control. 46(12), 1913–1918 (2001). Publisher Full Text OpenURL

  7. Smyshlyaev, A, Guo, B-Z, Krstic, M: Arbitrary decay rate for Euler-Bernoulli beam by backstepping boundary feedback. IEEE Trans. Autom. Control. 54(5), 1134–1140 (2009)

  8. Wang, JJ, Li, QH: Active vibration control methods of axially moving materials - a review. J. Vib. Control. 10, 475–491 (2004)

  9. Luo, ZH, Guo, BZ, Morgül, Ö: Stability and Stabilization of Infinite Dimensional Systems with Applications, Springer, London (1999)

  10. Kirchhoff, G: Vorlesungen über Mechanik, Teubner, Stuttgart (1883)

  11. Arosio, A, Panizzi, S: On the well-posedness of the Kirchhoff string. Trans. Am. Math. Soc.. 348(1), 305–330 (1996). Publisher Full Text OpenURL

  12. Ono, K: Global existence, decay, and blow up of solutions for some mildly degenerate nonlinear Kirchhoff strings. J. Differ. Equ.. 137, 273–301 (1997). Publisher Full Text OpenURL

  13. Shahruz, SM, Krishna, LG: Boundary control of a non-linear string. J. Sound Vib.. 195, 169–174 (1996). Publisher Full Text OpenURL

  14. Oplinger, DW: Frequency response of a nonlinear stretched string. J. Acoust. Soc. Am.. 32, 1529–1538 (1960). Publisher Full Text OpenURL

  15. Wu, YH, Xue, XP, Shen, TL: Absolute stability of the Kirchhoff string with sector boundary control (under review)

  16. Chen, G: Energy decay estimates and exact boundary value controllability for the wave equation in a bounded domain. J. Math. Pures Appl.. 58, 249–273 (1979)

  17. Chen, G: A note on boundary stabilization of the wave equation. SIAM J. Control Optim.. 19, 106–113 (1981). Publisher Full Text OpenURL

  18. Lagnese, JE: Note on boundary stabilization of wave equations. SIAM J. Control Optim.. 26, 1250–1256 (1988). Publisher Full Text OpenURL

  19. Shahruz, MS: Suppression of vibration in a nonlinear axially moving string by boundary control. ASME Des. Eng. Tech. Conf.. 106(1), 6–14 (1997)

  20. Li, TC, Hou, ZC: Exponential stabilization of an axially moving string with geometrical nonlinearity by linear boundary feedback. J. Sound Vib.. 296, 861–870 (2006). Publisher Full Text OpenURL

  21. Kobayashi, T: Boundary position feedback control of Kirchhoff’s non-linear strings. Math. Methods Appl. Sci.. 27, 79–89 (2004). Publisher Full Text OpenURL

  22. Komornik, V: Exact Controllability and Stabilization, Masson, Paris (1994)