Open access peer-reviewed chapter

General Stability in Viscoelasticity

Written By

Salim A. Messaoudi

Submitted: December 4th, 2015 Reviewed: May 12th, 2016 Published: September 21st, 2016

DOI: 10.5772/64217

Chapter metrics overview

1,358 Chapter Downloads

View Full Metrics


In this chapter, we consider a problem which describes the motion of a viscoelastic body and investigate the effect of the dissipation induced by the viscoelastic (integral) term on the solution. Precisely, we show that, under reasonable conditions on the relaxation function, the system stabilizes to a stationary state. We also obtain a general decay estimate from which the usual exponential and polynomial decay rates are only special cases.


  • general decay
  • memory
  • relaxation function
  • stability
  • viscoelasticity

1. Introduction

Elastic materials, when subjected to a suddenly applied loading state held constant thereafter, respond instantaneously with a state of deformation which remains constant. On the other hand, Newtonian viscous fluids respond to a suddenly applied state of uniform shear stress by a steady flow process. However, there exist materials for which any suddenly applied and maintained state of uniform stress produces an instantaneous deformation followed by a flow process which might or might not be limited in magnitude as time grows. Such materials exhibit both instantaneous elasticity effects and creep characteristics. Obviously, such a behavior cannot be described by either an elasticity theory or a viscosity theory only but it combines features of each. The most interesting examples of such materials are polymers, which can display all the intermediate range of properties (glassy, brittle solid or an elastic rubber or a viscous liquid) depending on temperature and the experimentally chosen time scale. Such materials are said to possess memories.

Many scientists, such as Maxwell, Kelvin, Voigt, and Boltzmann, have contributed in modeling these phenomena. Boltzmann, in 1874, supplied the first formulation of a three-dimensional theory of isotropic viscoelasticity. He elaborated the model of a “linear” viscoelastic solid on a basic assumption which states that at any (fixed) point x of the body, the stress at any time t depends on the strain at all the proceeding times. In addition, if the strain at all preceding times is in the same direction, then the effect is to decrease the corresponding stress. The influence of a previous strain on the stress depends on the time elapsed since that strain occurred and is weaker than those strains that occurred long ago. Such properties make the model of solid, elaborated by Boltzmann, a material with (fading) memory. These memory effects are expressed by the dependence on the deformation gradient. Therefore, for these “viscoelastic” materials, the stress at each point and at each instant does not depend only on the present value of the deformation gradient but on the entire temporal prehistory of the motion. In addition, Boltzmann made the assumption that a superposition of the influence of previous strains holds, which means that the stress-strain relation is linear. Mathematically, this is interpreted by the time convolution of a “relaxation” function with the Laplacian of the solution. As a consequence, a subtle damping effect is produced. The types of equations we intend to discuss in this chapter are of the form:


where is a bounded domain with regular boundary, g is a nonincreasing positive function, referred to as the relaxation function which describes the viscoelastic material in consideration, f is an external force, and u(xt) is the position of a point x “in the reference configuration” at a time t.

In early 1970s, Dafermos [1, 2] discussed a one-dimensional viscoelastic model, where he proved, for smooth monotonically decreasing relaxation functions, various existence and asymptotic stability results. However, no rate of decay has been given. After that, viscoelastic problems have attracted the attention of many researchers and many results of existence and long-time behavior have been established. To the best of our knowledge, the first work that studied the uniform decay of solutions was presented by Dassios and Zafirapoulos [3]. In their work, Dassios and Zafirapoulos presented a viscoelastic problem in ℝ3 and proved a polynomial decay for exponentially decaying kernels. In 1994, Muñoz Rivera [4] considered, in ℝn and in bounded domains, equations for linear isotropic homogeneous viscoelastic solids, with exponentially decaying memory kernels and showed that, in the absence of body forces, solutions decay exponentially for the bounded-domain case, whereas, for the whole space case, the decay is of a polynomial rate. After that, Cabanillas and Muñoz Rivera [5] studied problems, where the kernels are of algebraic (but not exponential) decay rates and showed that the decay of solutions is algebraic at a rate which can be determined by the rate of the decay of the relaxation function and the regularity of solutions. This result was later improved by Barreto et al. [6], where equations related to linear viscoelastic plates were treated. For viscoelastic systems with localized frictional dampings, Cavalcanti et al. [7] considered the following problem:


where is a bounded domain of ℝn (n ≥ 1) with a smooth boundary ∂g is a positive nonincreasing function satisfying, for two positive constants, the conditions:

- ξ 1 g t g t - ξ 2 g t , t 0 , E2

and a(x) ≥ a0 > 0 in a subdomain ω ⊂ , with meas(ω) > 0 and satisfying some geometry restrictions. They established an exponential rate of decay. Berrimi and Messaoudi [8] improved Cavalcanti’s result by weakening the conditions on both a and g. In particular, the function a can vanish on the whole domain and consequently the geometry condition is no longer needed. This result has been later extended to a situation, where a source is competing with the viscoelastic dissipation, by Berrimi and Messaoudi [9]. Also, Cavalcanti et al. [10] have studied a quasilinear equation, in a bounded domain, of the form:

| u t | ρ u t t - u - u t t + 0 t g t - τ u τ d τ - γ u t = 0 , E3

with ρ > 0, and a global existence result for γ ≥ 0, as well as an exponential decay for γ > 0, have been established. Messaoudi and Tatar [11,12] discussed the situation when γ = 0 and established polynomial and exponential decay results in the presence, as well as in the absence, of a nonlinear source term. Fabrizio and Polidoro [13] studied a homogeneous viscoelastic equation in the presence of a linear frictional damping (aut, a > 0) and showed that the exponential decay of the relaxation function g is a necessary condition for the exponential decay of the solution energy of the solution. In other words, the presence of the memory term, with a non-exponentially decaying relaxation function, may prevent the exponential decay even if the frictional damping is linear. He also obtained a similar result for the polynomial decay case.

For more general decaying kernels, Messaoudi [14,15] considered

u t t Δ u + 0 t g t τ Δ u τ d τ = b u m 2 u , in × 0 , u x t = 0 , x , t 0 u x , 0 = u 0 x , u t x 0 = u 1 x , x E1.1

with b = 0 and b = 1 and for relaxation functions satisfying

g ( t )  ξ ( t )g( t ), t  0,   E1.2

where ξ: ℝ+→ℝ+ is a nonincreasing differentiable function. He showed that the rate of the decay of the energy is exactly the rate of decay of g, which is not necessarily of exponential or polynomial decay type. After that, a series of papers using Eq. (1.2) have appeared. See for instance, Han and Wang [16], Liu [17,18], Park and Park [19], and Xiaosen and Mingxing [20].

In this work, we intend to study the following problem:

u t t Δ u + 0 t g t τ Δ u τ d τ + a u t m 2 u t = 0 , i n × 0 , u x t = 0 , x , t 0 u x , 0 = u 0 x , u t x 0 = u 1 x , x E1.3

where is a bounded and regular domain of ℝn, a > 0 is a constant, and g is a positive nonincreasing function satisfying Eq. (1.2). We will establish some general decay results depending on the behavior of g and m.


2. Preliminary

In this section, we present some material needed in the proof of our result and state a global existence result which can be proved using the well-known Galerkin method. See, for example, [2,3]. In order to prove our main result, we make the following assumptions:

(A1) g: ℝ+ → ℝ+ is a bounded differentiable function such that


where γ(t) is a differentiable function satisfying


(A2) Concerning the nonlinearity in the damping, we assume that


Remark 2.1. Examples of functions satisfying (A1) are


for a and b constants to be chosen properly.

Proposition 2.1. Let u0u1H01×L2 be given. Assume that (A1), (A2) hold. Then problem (1.3) has a unique global solution:


Proposition 2.2. [21] Let E: ℝ+ → ℝ+ be a non-increasing function and φ: ℝ+ → ℝ+ be an increasing C2-function such that


Assume that there exist q ≥ 0 and A > 0 such that


then we have

E( t )cE( 0 ) ( 1+φ( t ) ) 1 q , t0, if q>0,E16
E( t )cE( 0 ) e ωφ( t ) , t0, if q=0,E17

where c and ω are positive constants independent of the initial energy E(0).

Next, we introduce the “modified energy”:




Remark 2.2. By multiplying Eq. (1.3) by ut and integrating over Ω, using integration by parts and hypotheses (A1), (A2), we get, after some manipulations, as in [3,20],


3. Decay of solutions

In order to state and prove our main result, we set


where ε1 and ε2 are positive constants to be specified later and


Lemma 3.1. For ε1 and ε2 so small, we have


holds for two positive constants α1 and α2.

Proof. It is straightforward to see that


where Cp is the Poincaré constant. In the other hand,

F(t)E(t) ε1 2| ut |2dx ε1 2|u|2dx ε2 2| ut |2dx ε2 2Cp( 1l)( gou)(t)12lu( t)22+12ut22+12( gou)(t) +1 γ+2u γ+2 γ+2[ ε1+ε2 2]| ut |2dx ε1 2Cp| u|2dx ε2 2Cp( 1l)( gou)(t)α1E(t), E25

for ε1 and ε2 small enough. Thus, Eq. (3.2) is established.

Lemma 3.2. Assume that m ≥ 2 and assumptions (A1), (A2) hold. Then, the functional Ψ(t) satisfies, along the solution of Eq. (1.3), the estimate:


where C is a “generic” positive constant independent of t.

Proof. By using Eq. (1.3), we easily see that


We now estimate the third term of the right-hand side of Eq. (3.4), using Young’s inequality and (A2). Thus, we get


where cδ is a constant depending on δ. For the fourth term of the right-hand side of Eq. (3.4), we get


We then use Cauchy-Schwarz inequality, Young’s inequality, and the fact that


to obtain, for any η > 0,


By combining Eqs. (3.4)(3.7), we arrive at

Ψ(t)ut2dx12|u(t)|2dx+δCEm22(0)u22+12(1+η)(1l)2|u(t)|2dx+cδ|ut|mdx+12(1+η)(1l)0tg(tτ)|u(τ)u(t)|2dτdxut2dx12[1(1+η)(1l)22δCEm22(0)]|u(t)|2dx+12(1+η)(1l)(gou)(t)+cδ|ut|mdx. E32

By choosing η = l/(1 − l) and

(0), estimate (3.3) is established.

Lemma 3.3. Assume that m ≥ 2 and assumptions (A1), (A3) hold. Then, the functional Xt (t) satisfies, along the solution of Eq. (1.3) and for any δ > 0, the estimate


Proof: By using Eq. (1.3), we easily see that


Similarly to Eq. (3.4), we estimate the right-hand side terms of Eq. (3.9). So for any δ > 0, we have

Ω(0tg(ts)u(s)ds)(0tg(ts)(u(t)u(s))ds)dxδΩ(0tg(ts)|u(s)|ds)2dx+14δΩ(0tg(ts)|u(t)u(s)|ds)2dxδΩ(0tg(ts)(|u(t)u(s)|+|u(t)|)ds)2dx+14δ(0tg(ts)ds)Ω0tg(ts)|u(t)u(s)|2ds dx2δΩ(0tg(ts)|u(t)u(s)|ds)2dx+2δ(1l)2Ω|u|2dx+14δ(1l)(gu)(t)(2δ+14δ)(1l)(gu)(t)+2δ(1l)2Ω|u|2dx, E3.11
Ω|ut|m2ut(t)0tg(tτ)(u(t)u(τ))dτdxδΩ|0tg(tτ)(u(t)u(τ))dτ|mdx+cδΩ|ut|mdxδ(0tg(τ)dτ)m1Ω0tg(tτ)|u(t)u(τ)|mdτdx+cδΩ|ut|mdxδ(1l)m1C0tg(tτ)u(t)u(τ)2mdτ+cδΩ|ut|mdxδ(1l)m1Cp(2E(0)l)m220tg(tτ)u(t)u(τ)22dτ+cδΩ|ut|mdxδ(1l)m1Cp(2E(0)l)m22(gu)(t)+cδΩ|ut|mdx, E3.12

A combination of Eqs. (3.9)(3.13), then, yields Eq. (3.8).

Theorem 3.4. Let u0u1H01× L2 (Ω) be given. Assume that (A1), (A2) hold. Then, for any t0 > 0, there exist positive constants K and λ such that the solution of Eq. (1.3) satisfies


Proof: We start with the case m ≥ 2. Since g(0) > 0, then there exists t0 > 0 such that


By using Eqs. (2.3), (3.1), (3.3), and (3.8), we obtain


At this point, we choose δ so small that


Whence δ is fixed, the choice of any two positive constants ε1 and ε2 satisfying




We then pick ε1 and ε2 so small that Eqs. (3.2) and (3.16) remain valid and, further,


Therefore, we arrive at


for two constants cβ > 0. We multiply (3.17) by γ(t) and use Eqs. (1.2) and (2.3), to get


This implies that




Again, by using the fact that γ′(t) ≤ 0, letting


and noting that L~E, we arrive at


A simple integration of Eq. (3.18) over (t0t) leads to


We obtain, then, Eq. (3.14) by virtue of equivalence of E and L.

To establish Eq. (3.15), we re-estimate Eqs. (3.5) and (3.12), for m < 2, as follows


Similarly, we have


By repeating all above steps and using Eqs. (3.19), (3.20) instead of Eqs. (3.5), (3.12), we arrive at



which gives


By multiplying Eq. (3.21) by γ(t)Eq(t), for q > 0 to be specified later, and using (A1), Eq. (2.3), and Young’s inequality, we get


By choosing q = (2 − m)/(2m − 2) (hence, qm/(2 − m) = q + 1] and taking μ small enough, Eq. (3.22) yields


By recalling that γ ′ (t) ≤ 0 and integrating (3.23) over (ST), S ≥ t0, we get


for some positive constant A. Therefore, Proposition 2.2 gives (3.15). This completes the proof.

Remark 3.1. Estimates (3.14) and (3.15) also hold for t ∈ [0, t0] by virtue of continuity and boundedness of and γ.

Remark 3.2. This result generalizes and improves many results in the literature. In particular, it allows some relaxation functions which satisfy


instead of the usual assumption 1 ≤ ρ < 3/2.

Remark 3.3. Note that the exponential and the polynomial decay estimates, given in early works, are only particular cases of Eq. (3.14). More precisely, we obtain exponential decay for γ(t)≡a and polynomial decay for γ(t) = a(1 + t)− 1, where a > 0 is a constant.



The author would like to express his sincere thanks to King Fahd University of Petroleum and Minerals for its support. This work has been funded by KFUPM under project # IN151002.


  1. 1. Dafermos C.M., Asymptotic stability in viscoelasticity, Arch. Rational Mech. Anal. 37 (1970), 297–308.
  2. 2. Dafermos C.M., On abstract Volterra equations with applications to linear viscoelasticity, J. Differ. Equations 7 (1970), 554–569.
  3. 3. Dassios G. and Zafiropoulos F., Equipartition of energy in linearized 3-d viscoelasticity, Quart. Appl. Math. 48(4) (1990), 715–730.
  4. 4. Munoz Rivera J.E, Asymptotic behavior in linear viscoelasticity, Quart. Appl. Math. 52(4) (1994), 628–648.
  5. 5. Cabanillas E.L. and Munoz Rivera J.E., Decay rates of solutions of an anisotropic inhomogeneous n-dimensional viscoelastic equation with polynomial decaying kernels, Comm. Math. Phys. 177 (1996), 583–602.
  6. 6. Barreto R., Lapa E.C., and Munoz Rivera J.E., Decay rates for viscoelastic plates with memory, J. Elasticity 44(1) (1996), 61–87.
  7. 7. Cavalcanti M.M., Domingos Cavalcanti V. N., and Soriano J.A., Exponential decay for the solution of semilinear viscoelastic wave equations with localized damping, Electron. J. Differential Equations 2002(44) (2002), 1–14.
  8. 8. Berrimi S. and Messaoudi S.A., Exponential decay of solutions to a viscoelastic equation with nonlinear localized damping, Elect J. Diff. Eqns. 2004 (88) (2004), 1–10.
  9. 9. Berrimi S. and Messaoudi S.A., Existence and decay of solutions of a viscoelastic equation with a nonlinear source, Nonlinear Anal. TMA 64 (2006), 2314–2331.
  10. 10. Cavalcanti M.M., Domingos Cavalcanti V. N., and Ferreira J, Existence and uniform decay for nonlinear viscoelastic equation with strong damping, Math. Meth. Appl. Sci. 24 (2001), 1043–1053.
  11. 11. Messaoudi S.A. and Tatar N.E., Exponential and polynomial decay for a quasilinear viscoelastic equation, Nonlinear Anal. TMA 68 (2007), 785–793.
  12. 12. Messaoudi S.A. and Tatar N.E., Global existence and uniform stability of solutions for a quasilinear viscoelastic problem, Math. Meth. Appl. Sci. 30 (2007), 665–680.
  13. 13. Fabrizio M. and Polidoro S., Asymptotic decay for some differential systems with fading memory, Appl. Anal. 81(6) (2002), 1245–1264.
  14. 14. Messaoudi S.A. General decay of solutions of a viscoelastic equation, J. Math. Anal. Appl. 341 (2008), 1457–1467.
  15. 15. Messaoudi S.A., General decay of the solution energy in a viscoelastic equation with a nonlinear source, Nonlinear Anal. TMA 69 (2008), 2589–2598.
  16. 16. Han Xi. and Wang M., General decay of energy for a viscoelastic equation with nonlinear damping, Math. Meth. Appl. Sci. 32(3) (2009), 346–358.
  17. 17. Liu W.J., General decay of solutions to a viscoelastic wave equation with nonlinear localized damping, Ann. Acad. Sci. Fenn. Math. 34(1) (2009), 291–302.
  18. 18. Liu W.J., General decay rate estimate for a viscoelastic equation with weakly nonlinear time-dependent dissipation and source terms, J. Math. Phys. 50(11), art. no. 113506 (2009).
  19. 19. Park J. and Park S., General decay for quasilinear viscoelastic equations with nonlinear weak damping, J. Math. Phys. 50, art. 083505 (2009).
  20. 20. Xiaosen H. and Mingxing W., General decay of energy for a viscoelastic equation with nonlinear damping, Math. Methods Appl. Sci. 32(3) (2009), 346–335.
  21. 21. Komornik V., Exact Controllability and Stabilization. The Multiplier Method, Masson-John Wiley, Paris (1994).

Written By

Salim A. Messaoudi

Submitted: December 4th, 2015 Reviewed: May 12th, 2016 Published: September 21st, 2016