Skip to main content

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

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 [14]. In particular, boundary feedback stabilization of string and beam systems has become an important research area [57]. 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 t0, where a>0, b0 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 [1113] 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 t0. 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)dx
(4)

for all t0 (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)
Figure 1
figure 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 t0, 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 [1618]. 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 t0. 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], gC[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:RR:

  1. (H)

    There exist constants L 2 L 1 >0 and r1 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 wR,

{ 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 t0, the boundary control (7c) is the negative feedback of transversal velocity of the string at x=1.

Figure 2
figure 2

Linear sector constraint.

Figure 3
figure 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)=0for all t0.
(9)

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

E(t):= 1 2 0 1 y t 2 (x,t)dx+ 1 2 0 1 a y x 2 (x,t)dx+ 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 ) ] dx+ 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.1 Let y(,) be the solution for system (7a)-(7d). Then

0 1 y x y x t dx= y x (1,t) y t (1,t) 0 1 y x x y t dx,
(11)
0 1 x y x y x x dx= 1 2 y x 2 (1,t) 1 2 0 1 y x 2 (x,t)dx,
(12)
0 1 x y x y t t dx= 0 1 ( [ x y x y t ] t + 1 2 y t 2 ) dx 1 2 y t 2 (1,t).
(13)

Proof See the Appendix. □

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

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

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

for all t0.

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

E (t)= 0 1 y t y t t dx+ [ a + b 0 1 y x 2 d x ] 0 1 y x y x t dx.
(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 t0. 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)dt.

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

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

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

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

for all t0.

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 t0. On the other hand, the definition of energy function (10) implies

a 2 0 1 y x 2 dx+ 1 2 0 1 y t 2 dxE(t)

for all t0. It follows from the above inequality that

1 2 0 1 y x 2 dx+ 1 2 0 1 y t 2 dx { 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.3 Let G: R + R + be a non-increasing function. Assume that there exists a constant ω>0 such that

T + G q + 1 (t)dt 1 ω G(T)for all T0.

Then the following estimation is true, for all t0,

{ 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.4 The energy function E in equation (10), along the solution of system (7a)-(7d), satisfies

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

for all t0.

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 dxa, 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.5 For any constant q0, the energy function E along the solution of system (7a)-(7d) satisfies the following estimation, for all S>T0,

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 0T<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.1 Assume that assumption (H) holds. Then there exist three constants k 1 , k 2 ,σ>0 such that, for all t0,

{ 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)sfor all sR.

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

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 1 For any δ>0, we have the following estimates, for all S>T0,

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>T0,

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)dt 1 2 T S E q r + 1 r 1 (t)dt+ 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 =2max{1, 1 a }+ 1 2 ( 1 L 1 + L 2 a ), k 1 =eE(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.2 Assume that assumption (H) holds. Then there exist two constants k 1 , k 2 >0 such that for all x(0,1) and t0,

{ | 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 t0, 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 t0. 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.6sin(5x) and g(x)=0.5sin(3x). 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)=2xsinx,xR

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.

Figure 4
figure 4

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

Figure 5
figure 5

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

Figure 6
figure 6

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

Figure 7
figure 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 t0,

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 dx= 1 2 0 1 [ x y x 2 ] x dx 1 2 0 1 y x 2 dx= 1 2 y x 2 (1,t) 1 2 0 1 y x 2 dx

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

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

and

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

for all t0. 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 sR,

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.

References

  1. Chen G, Delfour MC, Krall AM, Payre G: Modeling, stabilization and control of serially connected beams. SIAM J. Control Optim. 1987, 25: 526-546. 10.1137/0325029

    Article  MATH  MathSciNet  Google Scholar 

  2. Lagnese J: Decay of solutions of wave equations in a bounded domain with boundary dissipation. J. Differ. Equ. 1983, 50: 163-182. 10.1016/0022-0396(83)90073-6

    Article  MATH  MathSciNet  Google Scholar 

  3. Shahruz SM: Bounded-input bounded-output stability of a damped nonlinear string. IEEE Trans. Autom. Control 1996, 41: 1179-1182. 10.1109/9.533679

    Article  MATH  MathSciNet  Google Scholar 

  4. Guo BZ: Riesz basis approach to the stabilization of a flexible beam with a tip mass. SIAM J. Control Optim. 2001, 39: 1736-1747. 10.1137/S0363012999354880

    Article  MATH  MathSciNet  Google Scholar 

  5. Fung RF, Wu JW, Wu SL: Exponential stabilization of an axially moving string by linear boundary feedback. Automatica 1999, 35: 177-181. 10.1016/S0005-1098(98)00173-3

    Article  MATH  MathSciNet  Google Scholar 

  6. Morgül Ö: Stabilization disturbance rejection for the beam equation. IEEE Trans. Autom. Control 2001, 46(12):1913-1918. 10.1109/9.975475

    Article  MATH  MathSciNet  Google Scholar 

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

    Article  MathSciNet  Google Scholar 

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

    Article  MATH  Google Scholar 

  9. Luo ZH, Guo BZ, Morgül Ö Communications and Control Engineering Series. In Stability and Stabilization of Infinite Dimensional Systems with Applications. Springer, London; 1999.

    Chapter  Google Scholar 

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

    MATH  Google Scholar 

  11. Arosio A, Panizzi S: On the well-posedness of the Kirchhoff string. Trans. Am. Math. Soc. 1996, 348(1):305-330. 10.1090/S0002-9947-96-01532-2

    Article  MATH  MathSciNet  Google Scholar 

  12. Ono K: Global existence, decay, and blow up of solutions for some mildly degenerate nonlinear Kirchhoff strings. J. Differ. Equ. 1997, 137: 273-301. 10.1006/jdeq.1997.3263

    Article  MATH  Google Scholar 

  13. Shahruz SM, Krishna LG: Boundary control of a non-linear string. J. Sound Vib. 1996, 195: 169-174. 10.1006/jsvi.1996.0414

    Article  MATH  MathSciNet  Google Scholar 

  14. Oplinger DW: Frequency response of a nonlinear stretched string. J. Acoust. Soc. Am. 1960, 32: 1529-1538. 10.1121/1.1907948

    Article  MathSciNet  Google Scholar 

  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. 1979, 58: 249-273.

    MATH  MathSciNet  Google Scholar 

  17. Chen G: A note on boundary stabilization of the wave equation. SIAM J. Control Optim. 1981, 19: 106-113. 10.1137/0319008

    Article  MATH  MathSciNet  Google Scholar 

  18. Lagnese JE: Note on boundary stabilization of wave equations. SIAM J. Control Optim. 1988, 26: 1250-1256. 10.1137/0326068

    Article  MATH  MathSciNet  Google Scholar 

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

    MathSciNet  Google Scholar 

  20. Li TC, Hou ZC: Exponential stabilization of an axially moving string with geometrical nonlinearity by linear boundary feedback. J. Sound Vib. 2006, 296: 861-870. 10.1016/j.jsv.2006.03.012

    Article  MATH  MathSciNet  Google Scholar 

  21. Kobayashi T: Boundary position feedback control of Kirchhoff’s non-linear strings. Math. Methods Appl. Sci. 2004, 27: 79-89. 10.1002/mma.440

    Article  MATH  MathSciNet  Google Scholar 

  22. Komornik V RAM: Research in Applied Mathematics. In Exact Controllability and Stabilization. Masson, Paris; 1994. Wiley, Chichester

    Google Scholar 

Download references

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).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Yuhu Wu.

Additional information

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.

Authors’ original submitted files for images

Rights and permissions

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

Reprints and permissions

About this article

Cite this article

Wu, Y., Wang, J. Exponential and polynomial stabilization of the Kirchhoff string by nonlinear boundary control. Bound Value Probl 2013, 215 (2013). https://doi.org/10.1186/1687-2770-2013-215

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2013-215

Keywords