`Journal of Applied MathematicsVolume 2010, Article ID 848620, 17 pageshttp://dx.doi.org/10.1155/2010/848620`
Research Article

## Numerical Exponential Decay to Dissipative Bresse System

Department of Mathematics, Federal University of Pará, Augusto Corrêa Street no. 01, Belém, CEP 66075-110, Pará, Brazil

Received 18 February 2010; Revised 25 May 2010; Accepted 24 June 2010

Copyright © 2010 M. L. Santos and Dilberto da S. Almeida Júnior. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We consider the Bresse system with frictional dissipative terms acting in all the equations. We show the exponential decay of the solution by using a method developed by Z. Liu and S. Zheng and their collaborators in past years. The numerical computations were made by using the finite difference method to prove the theoretical results. In particular, the finite difference method in our case is locking free.

#### 1. Introduction and Main Results

The circular arch problem is also known as the Bresse system (see [1] for details). Elastic structures of the arches type are objects of study widely explored in engineering, architecture, marine engineering, aeronautics and others. In particular, the free vibrations about elastic structures is a function of their natural properties and are an important subject of investigation in engineering and also in mathematics. In the field of mathematical analysis is interesting to know properties which relate the behavior of the energy associated with the respective dynamic model. For feedback laws, for example, we can ask what conditions about dynamic model must be ensured to obtain the decay of the energy of solutions in the time . In this sense, the property of stabilization has been studied for dynamic problems in elastic structures translated in terms of partial differential equations. A interesting property determines that the exponential decay with few feedback laws occurs only in an particular situation (see [2, 3]).

Regarding the numerical aspects related to elastic structures, special attention was dedicated for a numerical pathology know as locking phenomenon about transverse shear force or simply shear locking (see [47]). The locking phenomenon is a numerical anomaly that affect some numerical methods when applied in elastic structures, such as finite element methods for linear shape functions. Basically, this obstacle numerical consists in an over-estimation about of the coefficient of rigidity and that obviously does not correspond to the real case. Therefore, the locking phenomenon is a numerical inconsistent. In this sense, any numerical method applied to elastic structures must be to avoid this numerical anomaly.

Following the aim idea about deformation in elastic structures, we consider the Bresse system given by equations of motion where we use , and to denote the axial force, the shear force and the bending moment, respectively. These forces are stress-strain relations for elastic behavior and given by Here where is the density of material, is the modulus of elasticity, is the shear modulus, is the shear factor, is the cross-sectional area, is the second moment of area of the cross-section and is the radius of curvature. The functions , , and are the longitudinal, vertical and shear angle displacements, respectively (see Figure 1).

Figure 1: The circular arch.

From coupled equations (1.1) and (1.2) we obtain

In the literature (1.3) are the equations for the theory of circular arch. For more details see [1].

In this paper we will examine the issues concerning the asymptotic stabilization of Bresse system with frictional dissipative terms given by where and are positive constants. The initial conditions are and we assume the Dirichlet boundary conditions

Remark 1.1. If , then , and this model reduces to the well-known Timoshenko beam equations (see [1] for details).

Related to the objectives of this paper there are few results in the literature. Recently, Liu and Rao [8] studied the asymptotic behavior of the circular arch in the context of linear thermoelasticity and they proved that the exponential decay occurs if and only if the speed of wave propagation that occur in the model are equal.

Our main result is to prove that the Bresse system is exponentially stable in the presence of feedback laws and in the sequence we evidence our result by using the finite difference method of locking-free nature. This result was not considered for circular arch problem.

The paper is organized as follows. In Section 2 we establish the existence, regularity, and uniqueness of global solutions of (1.4)–(1.8). We use the semigroup techniques. In Section 3 we study the exponential decay of the strong solutions to system (1.4)–(1.8). We show the exponential decay of the solution by using a method developed by Liu and Zheng [9] and their collaborators in past years. Finally, in Section 4 some numerical aspects were considered. We use a particular discretization in finite difference method that is locking free (see [10, 11]) and our objective in this case is only to verify the numerical exponential decay for system (1.4)–(1.8). Questions about numerical analysis of numerical exponential decay as well the criterion stability for explicit time method in finite difference requires a more elaborate analysis.

#### 2. The Semigroup Setting

In this section we will study the existence and uniqueness of strong and global solutions for the system (1.4)–(1.8) using the semigroup techniques. To give an accurate formulation of the evolution problem we are introducing the product Hilbert space with norm given by Let , and we define the operator given by with domain

Therefore, system (1.4)–(1.8) is equivalent to where . It is not difficult to see that is a dissipative operator in the phase space . More precisely we have Under the above notation we can establish the following theorem.

Theorem 2.1. The operator is the infinitesimal generator of -semigroup of contraction in .

Proof. Since is dense in and is a dissipative operator, to prove Theorem 2.1, it is sufficient to prove that belongs to the resolvent set of , that is, (see [12, Theorem ]). Let us take , then there is an only such that It easily follows from (2.7) or the standard result on the linear elliptic equations that (2.7) has a unique solution . Therefore, and . The proof is now complete.

#### 3. Exponential Decay

In this section we will prove that the semigroup on is exponentially stable. Here we will use necessary and sufficient conditions for -semigroups being exponentially stable in a Hilbert space. This result was obtained by Gearhart [13] and Huang [14], independently (see also [15]).

Theorem 3.1. Let be a -semigroup of contractions on Hilbert space . Then is exponentially stable if and only if hold, where is the resolvent set of .

Theorem 3.2. The -semigroup of contractions , , generated by , is exponentially stable.

Proof. To prove the exponential stability of , it remains to verify the properties (3.1) and (3.2) of Theorem 3.1. First we prove that Suppose that the conclusion of (3.1) is not true. Then there is a such that, , (spectrum of and since is compact, must be an eigenvalue of . Let be with , such that . Using the definition of it follows that if and only if Multiplying equations (3.4) (3.5), and (3.6) by and , respectively, and performing integration by parts on , we arrive at which implies that and consequently from where we can conclude that , which contradicts . This completes the proof of (3.1).
We now prove (3.2) by a contradiction argument again. Suppose that (3.2) is not true. Then there are a sequence and a sequence such that From we have that there exists a unique sequence such that Denoting by it follows that Let us denote by From we have the following equations in : Taking the inner product of with in and then taking the real part, we obtain from where we can conclude that Considering the equation then from (2.6) which implies that , and then From we can conclude that From equations (3.14), (3.16) and (3.18) we obtain, after some calculations, In summary, we have proved that which contradicts . Thus, the proof is complete.

#### 4. Unlocking for Spatial Finite Difference Scheme

For numerical verification of the exponential decay of system (1.4)–(1.8) we use the total discretization in finite difference method. With respect to numerical schemes in this case, some aspects should be taking into account. The first one concerns to numerical phenomenon known as shear locking, which affects some numerical models applied to vibration problems in structures as shell, plates and beams. For the locking phenomenon (see [47]) for ample discussion for plane beams. For fast discussion about this numerical problem, we considerer the case of plane beams described by theory of Timoshenko and given by the following equations: are resulting of the equations (1.4)–(1.5) for .

By using the finite element method standard with linear shape functions, the rigidity coefficient is modified for where considering a rectangular geometry with width a and thickness . In particular, for plane beams must be have where is the spatial division (see ([10, 11] for details). Consequently, and clearly this value for rigidity coefficient does not correspond to the real case for fixed. In reality, with the equation for angle rotation in (4.2) is rewrite as and the energy respective is and obviously this energy is different of the real case with = . In general the numerical schemes in finite element method or finite difference method are locking free/over-estimation if not exist any additional term about the coefficients of the system, because the shear locking problem is basically a numerical anomaly characterized by a over-estimation about rigidity coefficients and dependent of parameter .

To avoid this numerical anomaly the finite element method can be used, however, special care must be taking in to account for choose of basis functions [6, 7]. Naturally, this shear locking/over-estimation can be to affect equations (1.4)–(1.6). See the studies by Loula et al. in [1618] for numerical treatment in finite element for circular arch problem.

In our case, we use the total discretization in finite difference method and to avoid shear locking/over-estimation we make a particular discretization for the functions of zero derivative such as and for a numerical operator of second order in relation to .

For our purposes, we use the space-time explicit method applied to equations (1.4)–(1.6) and we define , for and nets where and for and . The numerical scheme consists in to find and such that for , with the following numerical operators of second order for a function :

For and we denote the average of on the points and , respectively. This approximation avoid any over-estimation about the coefficients of equations (1.4)–(1.6). Then, we have,

The boundary conditions in the numerical context are given by and the discretizations to initial conditions are given by

The numerical energy associated for (4.7)–(4.13) is given by

Equations (4.7)–(4.9), in fact, are locking free because the numerical energy (4.14) is compatible with the continuous energy in (2.2), because the coefficients and are exactly those in (2.2), without any dependence with . To verify this, we have the following proposition.

Proposition 4.1. For all the energy (4.14) of solutions of the discrete equations (4.7)–(4.9), with initial conditions (4.13) and any boundary conditions (4.12) is such that

Proof. The proof is more extensive and we have omitted it here. For a idea of the proof, we use the multiplicative techniques such as performed in [19, 20], that is we multiply the equations (4.7)–(4.9) by and , respectively, and after we applied the discrete sum for . Then after some properly simplifications and taking in to account the boundary conditions given in (4.12) we get (4.14).

It is clear that discrete equations (4.7)–(4.9) are all consistent with for truncation error. Then, by Lax equivalence theorem the equations (4.7)–(4.9) are convergent if, and only if, they are stable. But another numerical limitation occurs for explicit time methods when applied to vibrations problem in elastic structures, with, say, respect to numerical stability. In particular, for explicit time method in finite difference applied to vibrations in Timoshenko beams, the prevailing numerical stability is given by where . The limitation occurs when is small. The understanding and overcoming for this limitation was studied by Joseph P. Wright in [10, 11]. Naturally, the same numerical problem affect (4.7)–(4.9), in function of the dependence of thickness . Then, this problem should be the objective of study in another opportunity, and for our purpose of numerical verification to exponential decay, we consider fixed.

For numerical example, we consider , thickness , width , , and the following initial conditions:

First we reproduce the conservative case, and for dissipative case, we consider . For the computational domain chosen 32 spatial points and 312 points in the time domain were given by seconds (see Figures 1, 2, 3, 4, 5, 6, and 7).

Figure 2: Conservative case, .
Figure 3: Conservative case, .
Figure 4: Total dissipation, .
Figure 5: Total dissipation, .
Figure 6: Total dissipation, .
Figure 7: Total dissipation, .

With Figures 8, 9, 10, 11, 12, and 13 we illustrate the simulation results for vibration in angle rotation and the behavior of numerical energy in the cases conservative and dissipative. Of course the numerical energy (4.14) must be conservative in the absence of damping.

Figure 8: Conservative case, .
Figure 9: Conservative case, .
Figure 10: Total dissipation case, .
Figure 11: Total dissipation, .
Figure 12: Conservative case.
Figure 13: Dissipative case.

#### Acknowledgments

The first author acknowledges the support of the CNPq (306338/2008-4) and FAPESPA-PA-Brazil. The authors are thankful to the referees of this paper for the valuable suggestions which improved the paper.

#### References

1. J. E. Lagnese, G. Leugering, and E. J. P. G. Schmidt, Modeling, Analysis and Control of Dynamic Elastic Multi-Link Structures, Systems & Control: Foundations & Applications, Birkhäuser, Boston, Mass, USA, 1994.
2. J. E. Muñoz Rivera and R. Racke, “Mildly dissipative nonlinear Timoshenko systems—global existence and exponential stability,” Journal of Mathematical Analysis and Applications, vol. 276, no. 1, pp. 248–278, 2002.
3. J. E. Muñoz Rivera and R. Racke, “Global stability for damped Timoshenko systems,” Discrete and Continuous Dynamical Systems. Series A, vol. 9, no. 6, pp. 1625–1639, 2003.
4. A. F. D. Loula, T. J. R. Hughes, and L. P. Franca, “Petrov-Galerkin formulations of the Timoshenko beam problem,” Computer Methods in Applied Mechanics and Engineering, vol. 63, no. 2, pp. 115–132, 1987.
5. A. F. D. Loula, T. J. R. Hughes, L. P. Franca, and I. Miranda, “Mixed Petrov-Galerkin methods for the Timoshenko beam problem,” Computer Methods in Applied Mechanics and Engineering, vol. 63, no. 2, pp. 133–154, 1987.
6. G. Prathap and G. R. Bhashyam, “Reduced integration and the shear-flexible beam element,” International Journal for Numerical Methods in Engineering, vol. 18, no. 2, pp. 195–210, 1982.
7. T. J. R. Hughes, R. L. Taylor, and W. Kanoknukulchai, “A simple and efficient finite element method for plate bending,” International Journal for Numerical Methods in Engineering, vol. 11, no. 10, pp. 1529–1543, 1977.
8. Z. Liu and B. Rao, “Energy decay rate of the thermoelastic Bresse system,” Zeitschrift für Angewandte Mathematik und Physik, vol. 60, no. 1, pp. 54–69, 2009.
9. Z. Liu and S. Zheng, Semigroups Associated with Dissipative Systems, vol. 398 of Chapman & Hall/CRC Research Notes in Mathematics, Chapman & Hall/CRC, Boca Raton, Fla, USA, 1999.
10. J. P. Wright, “A mixed time integration method for Timoshenko and Mindlin type elements,” Communications in Applied Numerical Methods, vol. 3, no. 3, pp. 181–185, 1987.
11. J. P. Wright, “Numerical stability of a variable time step explicit method for Timoshenko and Mindlin type structures,” Communications in Numerical Methods in Engineering, vol. 14, no. 2, pp. 81–86, 1998.
12. A. Pazy, Semigroups of Linear Operators and Applications to Partial Differential Equations, vol. 44 of Applied Mathematical Sciences, Springer, New York, NY, USA, 1983.
13. L. Gearhart, “Spectral theory for contraction semigroups on Hilbert space,” Transactions of the American Mathematical Society, vol. 236, pp. 385–394, 1978.
14. F. L. Huang, “Characteristic conditions for exponential stability of linear dynamical systems in Hilbert spaces,” Annals of Differential Equations, vol. 1, no. 1, pp. 43–56, 1985.
15. J. Prüss, “On the spectrum of ${C}_{0}$-semigroups,” Transactions of the American Mathematical Society, vol. 284, no. 2, pp. 847–857, 1984.
16. A. F. D. Loula, L. P. Franca, T. J. R. Hughes, and I. Miranda, “Stability, convergence and accuracy of a new finite element method for the circular arch problem,” Computer Methods in Applied Mechanics and Engineering, vol. 63, no. 3, pp. 281–303, 1987.
17. T. J. R. Hughes, R. L. Taylor, and W. Kanoknukulchai, “A simple and efficient finite element for plate bending,” International Journal for Numerical Methods in Engineering, vol. 11, no. 10, pp. 1529–1543, 1977.
18. B. D. Reddy and M. B. Volpi, “Mixed finite element methods for the circular arch problem,” Computer Methods in Applied Mechanics and Engineering, vol. 97, no. 1, pp. 125–145, 1992.
19. M. Negreanu and E. Zuazua, “Uniform boundary controllability of a discrete 1-D wave equation,” Systems & Control Letters, vol. 48, no. 3-4, pp. 261–279, 2003.
20. W. Strauss and L. Vazquez, “Numerical solution of a nonlinear Klein-Gordon equation,” Journal of Computational Physics, vol. 28, no. 2, pp. 271–278, 1978.