Open access peer-reviewed chapter

Non-Classical Memory Kernels in Linear Viscoelasticity

By Sandra Carillo and Claudio Giorgi

Submitted: December 6th 2015Reviewed: May 16th 2016Published: September 21st 2016

DOI: 10.5772/64251

Downloaded: 1326


In linear viscoelasticity, a large variety of regular kernels have been classically employed, depending on the mechanical properties of the materials to be modeled. Nevertheless, new viscoelastic materials, such as viscoelastic gels, have been recently discovered and their mechanical behavior requires convolution integral with singular kernels to be described. On the other hand, when the natural/artificial aging of the viscoelastic material has to be taken into account, time-dependent kernels are needed. The aim of this chapter is to present a collection of nonstandard viscoelastic kernels, with special emphasis on singular and time-dependent kernels, and discuss their ability to reproduce experimental behavior when applied to real materials. As an application, we study some magneto-rheological elastomers, where viscoelastic and magnetic effects are coupled.


  • Materials with memory
  • Viscoelasticity
  • Unbounded memory kernels
  • Existence of solutions
  • Asymptotic solutions’
  • Behavior

1. Introduction

The stress-strain relation in linear viscoelasticity involves a convolution integral with a memory kernel. The fading memory principle requires that the memory kernel decays quickly as the elapsed time goes to infinity, but no limitation is imposed to its behavior near zero. So, a wide range of kernels may be used depending on the nature of the materials to be modeled. Starting from the rheological model of a standard viscoelastic solid, whose kernel involves a single exponential, a large variety of regular kernels have been classically employed: discrete and continuous Prony series, completely monotonic functions, etc. Recently, new viscoelastic materials, such as viscoelastic gels, have been described by virtue of convolution integral with singular kernels: for instance, fractional and hypergeometric kernels [1]. On the other hand, when the natural/artificial aging of the viscoelastic material has to be taken into account, time-dependent kernels are needed. Furthermore, the behavior of some new materials, for instance, ferrogel and magneto-rheological elastomers, can be determined by coupling viscoelastic and magnetic effects.

The material of this chapter is organized as follows. First we present the model of a viscoelastic body which represents the basis for our study. It is assumed to be homogeneous and isotropic, and its crucial feature is that the stress response at time tlinearly depends on the whole past historyof the strain up to t. Then, we look for the modeling of aging isothermal viscoelasticity, assuming that the viscoelastic structural parameters are time dependent while the material is subject to chemical or physical agents at constant temperature. Finally, singular kernel problems are addressed to, at first, in the case of a viscoelastic body and, later, when the viscoelastic behavior is coupled with magnetization. In particular, the case of magneto-viscoelastic bodies is considered. Indeed, the idea of coupling the viscoelastic behavior with magnetic effects is suggested by new materials which are obtained by inserting magnetic defects into a solid body to have the opportunity to influence the mechanical properties of the body when a magnetic field is applied.


2. Preliminary notions and notations

This section is devoted to provide the key notions concerning the model of isothermal viscoelastic body with memory. For sake of simplicity, the body is supposed homogeneous. In order to briefly introduce the subject, at the beginning, we restrict our attention to one-dimensional processes. Let εdenote the uniaxial strain and σthe corresponding tensile stress at every point xof the reference configuration of the sample. According to Boltzmann’s formulation of hereditary elasticity [2], a linear viscoelastic solid may be described by a stress-strain relation in the Riemann-Stieltjes integral form


where Gis named Boltzmann function(or memory kernel) and ε(·) is a fading strain history, namely


In particular, when the strain history vanishes from −∞ to 0, then Eq. (1) reduces to


A peculiar behavior of viscoelastic solid materials is named relaxation property: if the solid is held at a constant strain starting from a given time t0 ≥ 0, the stress tends (as t→ ∞) to a constant value which is “proportional” to the applied constant strain. Indeed, if ε(x, ·) is continuous on (−∞, t0] and


it follows that


where the relaxation modulus


is assumed to be positive. Then, using Eq. (2) and letting


the stress-strain relation (1) may be rewritten as


Of course, the choice of Gis required to satisfy some basic principles, like the fading memory principleand the dissipation principle, a thermostatic version of the second law of thermodynamics (see [3], for instance). In general, these conditions allow the memory kernel to be unbounded at the origin.

In the terminology of Dautray and Lions [4], hereditary effects with long memory rangeare represented by a convolution integral, where


whereas a short memory rangeis related to singular kernels of the Dirac delta type. In the latter case, letting Ĝ= Γδ0, where δ0denotes the Dirac mass at 0+, from Eq. (1) it follows:


where t denotes partial derviation with respect to ‘t’.

which is named the Kelvin-Voigt model. On the other hand, assumption (6) may be strengthened by letting Gbe bounded along with its derivatives


If this is the case,


and an integration by parts changes Eq. (1) into the alternate forms


where the relaxation function G'(τ)is the derivative with respect to τof the Boltzmann function G. This constitutive stress-strain relation is based on the Lebesgue representation of linear functionals in the history space theory devised by Volterra [5]. Provided that Eq. (2) holds true, the Boltzmann and the Volterra constitutive relations are equivalent. The latter approach, however, can be applied to a wilder class of strain histories (uniformly bounded, for instance), in which Eq. (2) is no longer needed.

In the three-dimensional case, all fields depend on the space-time pair (x, t) ∈ Ω × ℝ, where Ω ⊂ ℝ3 is the reference configuration. The displacement vector u(x, t) is given by


where μ(x, ·) is the motion of x, and


is the infinitesimal strain tensor. Borrowing from Eq. (9), the viscoelastic Cauchy stress tensor Tis given by



:ℝ+ → Lin(Sym) stands for the relaxation function and


A simple manipulation of Eq. (10) yields an alternate stress-strain relation:


For fading strain histories obeying Eq. (2), an integration by parts allows Eq. (10) to be rewritten as



is defined as Ĝin Eq. (4). The material is said to enjoy the fading memory principlewhen, for every ε > 0there exists a positive time shift S0(ε), possibly dependent on the strain history, such that

|0G(τ+s)E(x,tτ) dτ |=| 0G^(τ+s)dτE(x,tτ) |ε ,ss0.E13

Note that Eq. (13) does work even if G(t) is allowed to be singular and non-integrable at the origin. Indeed, the fading memory property requires that the memory kernel decays quickly as the elapsed time τgo to infinity, but no limitation is imposed to its behavior near zero.

3. Aging models in linear viscoelasticity

Aging is a gradual process in which the properties of a material change, over time or with use, due to chemical or physical agents. Corrosion, obsolescence, and weathering are examples of aging. In metallurgical processes, aging may be induced by a heat treatment (age hardening). Consequences of aging are of various types. For instance, the damages caused by melting or time-deteriorating processes are examples for decreasing stiffness in elastic springs. Instead, solidification of concrete is an irreversible transition process where the system increases its stiffness and releases a large amount of energy per volume. As pointed out in the sequel, the former type of aging is compatible with thermodynamics under isothermal conditions, while the latter involves a latent heat and then requires a non-isothermal framework. For definiteness, in this section, we investigate viscoelastic solids and assume that the viscoelastic model holds while the material is subject to chemical or physical agents at constant temperature. It is then understood that we look for the modeling of aging isothermal viscoelasticity.

In modeling aging effects, we might think that in Eq. (1), the dependence of the Boltzmann function Gon tand sis not merely through the difference t − sbut involves tand sseparately. It is a central problem to understand how to model Gand we would like to argue as far as possible on physical grounds. The recourse to physical arguments to model aging properties is not new in the literature (see, e.g., [6, 7]). Quite naturally one may refer to the classical rheological models with regular kernels (8) and hence to express the aging properties in terms of time-dependent elasticity and viscosity coefficients.

To this purpose, we first address attention to rheological models and, in particular, we consider the standard solid and the Wiechert-Maxwell model [8, 9]. Hence, we establish the functional providing the stress in terms of the strain. This procedure has the advantage of showing how the dependence on the present value and that on the history of ε are influenced by the rheological parameters. Next we generalize the model and look for the corresponding three-dimensional version. For a generic time-dependent relaxation function, a free energy is found to hold for the stress functional as a suitable Graffi-Volterra functional [3, 5, 10]. As a consequence, the stress functional is found to be compatible with thermodynamics subject to weak restrictions on the relaxation function.

3.1. Insights from a rheological model

To get some insights about the modeling of aging viscoelastic solids, we start from the classical standard linear solid where a Maxwell unit, consisting of a spring and a dashpot connected in series, is set in parallel with a lone spring. While we have in mind the behavior of the model in terms of elongation and forces, we extend the formulae to the continuum framework by the standard analogies stress-force and strain-elongation. It is understood that the model is framed within a one-dimensional picture, so that both strains and stresses are scalar fields depending on (x,t). Since the elastic and Maxwell elements are in parallel, the strain is the same for every element and the applied stress is the sum of the stress in each element (see Figure 1).

Figure 1.

Mechanical scheme of an aging standard viscoelastic solid.

Hereafter, the dependence on xof all the fields involved is understood and not written. For the Maxwell element, let εsand εdbe the strain of the spring and that of the dashpot. Hence, denoting by ε the common strain we have


Let σe be the stress on the isolated spring while σm the stress on the Maxwell element. Then, the total applied stress is given by


Moreover let kand ke be the elastic modulus (or rigidity) of the spring of the Maxwell element and of the spring in parallel, respectively, and γthe viscosity of the dashpot (Figure 1). It is the essential feature of the aging effect that k, ke and γare positive functions of the time t. In the Maxwell unit, the spring and the dashpot are in series and hence they are subject to the same stress so that, according to the Hook’s law,


where tdenotes partial differentiation with respect to time t. Using the last equality, from Eq. (14) we have


Incidentally, if ke and kare time independent then time differentiation of Eq. (14) and use of Eq. (15) give


which holds for any viscoelastic standard element. Letting g0 = ke + kand g=ke, this differential equation is equivalent to


which is commonly used in the literature. So far it is only assumed that

(A1) ke, k, γ∈ℂ1 (ℝ)  and  ke(t), k(t)0,γ(t)γ0>0for every t∈ℝ.

(A2) tα(ξ)dξ=

The last condition is fulfilled when α = k/γis a constant function, for instance.

We may regard Eq. (16) as a differential equation in the unknown εd(t). Then, integration over [t0,t] yields


It is convenient to let t0→−. By assuming that εd is uniformly bounded on (−,t], assumption (A2) allows us to take


Hence, we have


and from the representation


we obtain the stress-strain relation


which involves both the present value ε(t)and the past history ε(s),s[,t). Since


an integration by parts allows (17) to be rewritten as


provided that εis uniformly bounded on [−,t). A change of variables τ = t−swithin Eq. (17) leads to the alternate form


Finally, after introducing the so-called relative history,


the stress-strain relation may be rewritten as


3.2. Some remarks on the aging effect

To give some evidence to the aging effects, we fix a time t0 < tand we let


This statement holds even if t0 = −provided that we identify the constant values with the limits as t → −. If no aging affects the material, then


Otherwise, remembering that α = k/γ, we introduce the functions


In particular, κ, ϰand wequal unity for non-aging materials. This approach leads to identify κand wwith the aging factors of the elastic and the Maxwell elements, respectively. Moreover, Eq. (18) becomes




This suggests that aging effects can be modeled by means of two functions: κand H. In our notation, the present value ε(t)is affected by the factor κ(t), whereas the history of ε is affected by the function H(t,s). Letting


the stress-strain relation (21) may be rewritten as


For non-aging materials, κ(t) = H(t,s) ≡ 1, and this relation reduces to Eq. (5).

We end by observing that in [11] fatigue effects are modeled by using the convolution form (3) modified by the occurrence of a reduced time tr in place of time t, that is




aT being named the time-temperature shift factor. A similar approach may equally well model other aging effects. If we denote by f(t)the function associated with the aging process applied to the body, then we may introduce a rescaled time tr which is given by


To our mind the use of a rescaled time tr is an operative way of accounting for aging effects. Hereafter, we show that Eq. (18) may be represented as a linear convolution integral after introducing a suitable rescaled time. We start by letting


where, for every fixed t,


is positive and nondecreasing because of (A1). Moreover, from (A2) limsA(s)=, so that tr=A(t)plays the role of a rescaled time. Letting ε^(A(s))=ε(s), we have


and the stress-strain relation (18) may be rewritten as


This expression suggests that aging effects may be partly represented by a suitable change of the time scale within the memory integral. Indeed,


where tr=A(t), sr=A(s), and σ^(tr)=σ(t). This expression completely matches with Eq. (23) only if ke and kare constants. For non-aging materials, the scaling turns out to be linear, tr=A(t)=αt, and Eq. (24) becomes


3.3. From long to short memory: a possible aging effect

Instead of (A1) and (A2), we assume here

(B1) keγ > 0 are constants.

(B2) k ∈ ℂ1(ℝ) is positive, nondecreasing and such that


From (B1), the viscosity of the damper and the rigidity of the lone spring are constants, whereas (B2) translates the fact that the spring in the Maxwell element becomes completely rigid in the longtime. Under the additional very mild assumption

It is easily seen that (B3) always holds, for instance, when kis eventually concave down as t.

(B3) limtk(t)[k(t)]2=0,

we can prove that the Kelvin-Voigt viscoelastic model (7) is recovered when t →. Namely, letting


within (B1)-(B3) the distributional convergence


occurs as t →, so that Eq. (18) collapses into the Kelvin-Voigt stress-strain relation


The rigorous proof can be found in [12]. Since the function kt(·) is nonnegative for every t, Eq. (25) follows by showing that, for every fixed v≥ 0,


Assumptions (B1)-(B3) comply with the dissipation principle, as proved by Example 2 in Section 3.7.

3.4. The Wiechert-Maxwell model with aging

The Wiechert-Maxwell model (or Generalized Maxwell model) is composed by a bunch of (say N) Maxwell elements, assembled in parallel, and a further spring in parallel with the whole array. Since all elements are in parallel the strain is the same for every element and the applied stress is the sum of the stress in each element.

εdenotes the common strain and σe denotes the stress on the isolated spring, while σ1,…,σN are the stresses on the Maxwell pairs. Moreover, let ke, k1,…, kNbe the elastic modulus (or rigidity) of the N+ 1 springs and γ1,…,γNthe viscosity coefficients of the dashpots. It is the essential feature of the aging effect that ke, k1,…, kNand γ1,…,γNare functions of the time t(see Figure 2).

Figure 2.

Mechanical scheme of a Wiechert-Maxwell model.

The dependence of kand γon time requires that we review the elementary arguments to determine the relations among σe, σ1,…,σNand ε. For each jth Maxwell element, let εsjand εdjbe the strain of the spring and that of the dashpot. Hence, we have


As a consequence, γjtεdj+ kjεdj= kjε, and then


Previous assumptions are generalized for any j= 1,…,N, as follows:

(C1) ke, kj, γj∈ ℂ1(ℝ) and ke(t) ≥ 0, kj(t) ≥ 0, γj(t) ≥ γ0 > 0, for every t∈ ℝ,

(C2) tαj(ξ)dξ=, for every t∈ ℝ.

By regarding Eq. (26) as a differential equation in the unknown εdj(t), for the jth Maxwell element we have


The whole stress on the Wiechert-Maxwell model is then given by


To give some evidence to the aging effects as in Section 3.2, we assume that all springs and dashpots within the Maxwell elements have common aging factors.

(C3) There exist a time t0 ∈ ℝ (possibly, t0 = -) and two functions ϰ, w: ℝ → ℝ such that ϰ(t0) = w(t0) = 1 and for every j= 1,…,N


In particular, j=1Nαj(y)=j=1Nαjw(y). Defining the aging factors as follows:


the stress-strain relation (27) may be rewritten in the form (22) by letting


As in the standard solid model, the present value ε(t) is affected by the factor κ(t) only, whereas the history of εis affected by the function H(t, s).

3.5. Time-dependent linear viscoelasticity

Borrowing from the Wiechert-Maxwell solid developed above, we now state the uniaxial stress-strain constitutive equation that allows for time-dependent properties. If we introduce the function


which is defined on the half plane = {(t, s) ∊ ℝ2 :st}, the Wiechert-Maxwell constitutive law (27) may be rewritten into the general form


where ∂sdenotes partial differentiation with respect to the variable sand


for all t∈ ℝ. In addition, from (C1)–(C2) we have


For further convenience, we define Ğ: ℝ × ℝ+ → ℝ as Ğ(t, τ) = G(t, t- τ) so that


Finally, remembering that G0(t) = G(t, t), an integration by part of (28) yields


provided that Eq. (2) holds and


The classical expressions (5) and (9) are recovered from Eqs. (29) and (28), respectively, by simply assuming that G(t, s) = G(t-s), st. If this is the case, G0 and Gturn out to be constants.

We now look for a general, though linear, time-dependent three-dimensional model. According to Eq. (28), the Cauchy stress tensor Tis given by



stands for the t-dependent relaxation function and



: ×+Lin(Sym) such that


a change of the integration variable into Eq. (30) yields an alternate stress-strain relation


or equivalently




For non-aging materials, (t, s) and

(t, τ) reduce to (t-s) and (τ), respectively.

Hereafter, for ease in writing, we introduce the function


which is assumed to satisfy the following properties.

  1. (M1) 

    for every compact set ⊂ ℝ × ℝ+.

  2. (M2) For every fixed t∈ ℝ, the map τG(t,τ)is positive semi-definite, absolutely continuous and summable on ℝ+. Then, for every t∈ ℝ


    Besides, it is differentiable for all τ∈ ℝ+ and


    for every compact set

    ⊂ ℝ × ℝ+.

  3. (M3) For every fixed τ> 0, the map tG(t,τ)is differentiable for all t∈ ℝ. Besides,


    for every compact set

    ⊂ ℝ × ℝ+.

  4. (M4) There exists a nonnegative scalar function M: ℝ → ℝ+, bounded on bounded intervals, such that


    for every (t, τ) ∈ ℝ × ℝ+.

According to (M2), the t-dependent relaxation function

may be represented as


Borrowing from the scalar case,

(t) is assumed to be positive definite for every t∈ ℝ, namely


Finally, an integration by part of Eq. (30) yields


provided that Eq. (2) holds for Eand


As an advantage, within Eq. (34),

may be unbounded at the origin.

In order to stress the aging effects, we might assume the following factorization of the memory kernel G.

(M5) There exist three functions,

such that ℍ is uniformly bounded,
and, for every t∈ ℝ and s< t,


Accordingly, the stress-strain relation (34) may be rewritten into the form (22). The aging factors κand ℍ reduces to unit when non-aging materials are considered.

So far, we restrict our attention to scrutinize stress-strain relations in the form (30). In particular, for isotropic materials Gtakes the special form


where 1is the unit second-order tensor, Iis the symmetric fourth-order identity tensor, and λ, μ: ℝ × ℝ+ are named Lamé relaxation functions. Accordingly,


3.6. A Wiechert-type three-dimensional model

In the sequel, we scrutinize the special isotropic vector-valued kernel G=G111+G2I, where 1and 2are given by


as in the rheological Wiechert-Maxwell model devised in Section 3.1. We first prove that properties (M1)-(M4) hold provided that some additional restrictions are imposed on the material functions kjiand αji. Finally, we give some examples of these functions that fulfill these conditions.

  • (M1) Starting from (C1), it is quite trivial to prove this property.

  • (M2) By virtue of Eq. (35) and (C1), Gi, i=1,2, are positive and continuously differentiable with respect to tand τ. Moreover,



(t), is summable and vanishing at infinity for every t∈ ℝ. In addition, we have


Hence (M2) is fulfilled.

  • (M3) It is obviously true as αji∈ ℂ1(ℝ) by virtue of (C1)-(C2). In particular,


  • (M4) In order to prove this property we need more restrictive conditions. Since


    a sufficient condition to ensure (M4) is given by


  • In order to prove these inequalities, we now assume


    and for every t∈ ℝ we let


    It is apparent that M(t) ≥ 0 and then from (35) it follows:


    When non-aging material parameters are involved, Eq. (37) reduces to τi+ Mi≤ 0, i= 1, 2, which implies the exponential decay of the kernels.

    3.7. Some examples

    We present here some special expressions of material functions kjiand αji, j= 1, 2,…, N, i= 1, 2, which fulfill properties (M1)-(M4).

    • Example 1.

      For simplicity, we restrict our attention to a single Maxwell element. Letting j= 1. and i= 1, 2, we choose


    and then


    so that (A1) and (A2) hold true. Condition (38) is fulfilled for all t∈ ℝ and for every choice of the parameters, so that


  • Example 2.

    Otherwise, for j= 1. and i= 1,2, we can choose


    where ωi, ki, ηi> 0. Accordingly,


    so that (A1) and (A2) hold true. On the other hand, condition (38) reduces to


    which is equivalent to


    and is fulfilled for all t∈ ℝ provided that ωi≤ 3 Ki/ ηi. If this is the case,


  • 3.8. Motion, free energies, and thermodynamics

    We now derive the motion equation related to the time-dependent viscoelastic stress-strain relation (32) and we examine its compatibility with thermodynamics. The displacement field u: Ω×ℝ → ℝ3, relative to the reference configuration Ω ⊂ ℝ3, is subject to the equation of motion


    where fis the body force, per unit volume. Hence, from Eqs. (32)–(33), we obtain


    In order to introduce the initial boundary value problem for this equation, we have to take in mind that it is not invariant under time shift.

    Consistent with linear viscoelasticity, we restrict attention to isothermal processes, namely those where the temperature is constant and uniform. Hence, the local form of the second law inequality reduces to the dissipation inequality


    where ρis the mass density,

    is the Helmholtz free energy density per unit volume, and Dis the stretching tensor. Again for consistency with the linearity of the model, we let the mass density ρbe constant and take the approximation


    Accordingly, we take the dissipation inequality in the form


    In materials with memory, the motion equations are required to rule both the displacement instantaneous value u(t) and its history up to t. Letting t0 ∈ ℝ be arbitrarily fixed, we define the relative displacement historyζt(x,s), with (t,s) ∈ [t0,T] × ℝ+, by


    where ζt0 is the prescribed initial (relative) past history of uup to t0,


    Accordingly, ζt(x,0) = 0 and the motion equation (39) becomes a system


    where u: Ω×[t0, +∞) and ζt: Ω×ℝ+ → ℝ3, T∈ [t0, +∞) are the unknown variables. Their initial conditions are prescribed at t0 ∈ ℝ as follows


    Let H0= [L2 (Ω)]3 and H1= [H10(Ω)]3, and let ,jdenote the usual inner product in Hj, j= 0,1. For every tt0, we introduce the family of memory spaces



    denotes the t-dependent weighted L2 inner product equipping each
    . In this functional framework, the motion equation admits a unique regular solution. The proof of this result can be found in [12, Th. 4.5].

    Theorem 1.Lett = H1 × H0 ×

    and fH0. Under assumptions (M1)-(M4), for every T> t0 and every initial datum zt0=(ut0,vt0,ζt0)Ht0, problem (42)-(43)admits a unique solution z(t) = (u(t), ∂tu(t), ζt) on the interval [t0,T] such that




    for some C> 0 depending only on T, t0 and the size of the initial datumzt0.

    Now, we introduce a time-dependent free energy density borrowing its expression from the Graffi’s single-integral quadratic form (see [13] and references therein). Let


    For ease in writing, hereafter the dependence on xis understood and not written. In addition, we assumeρ= 1. After integrating over Ω, we end up with the total free energy functional


    Theorem 2.For an aging viscoelastic material, the dissipation inequality (40) is fulfilled provided that (M4) holds and


    Proof.First we observe that


    Then, by virtue of (42), and some integration by parts, we obtain


    and, taking into account (41),


    In summary, we end up with


    Owing to (M4), this yields


    where M(t) ≥ 0, and (44) finally implies the dissipation inequality


    4. Singular kernel models in linear viscoelasticity

    The study of singular kernel problems is motivated by the modeling of new materials and, in particular, of the mechanical behavior of some new viscoelastic polymers and bio-inspired materials. As noticed in [14], the appropriate way to handle the response of certain time-dependent systems exhibiting long tail memories is to account for power laws, both for creep and relaxation, leading to the occurrence of fractional hereditariness. Another example encountered in natural materials is mineralized tissues as bones, ligaments, and tendons. They exhibit a marked power-law time-dependent behavior under applied loads (see e.g. [15]), since the high stiffness of the crystals in such tissues is combined with the exceptional hereditariness of the collagen protein-based matrix. In all these cases, we are forced to abandon the regularity assumptions (8) and assume the memory kernels obey Eq. (6) and are unbounded at the origin.

    The idea of singular kernels to model particular cases of viscoelastic behaviors was introduced by Boltzmann [2] in the nineteenth century. The fast growth of polymer science motivated further developments of viscoelasticity in the middle of the twentieth century [16, 17], but a Volterra-type integro-differential equation with a regular kernel (typically, a finite sum of exponentials) was preferred to the Boltzmann approach in the modeling of the mechanical response [5, 18]. Later, however, many authors addressed their interest to singular kernel problems, both under the analytical as well as the model point of view [1924], and their thermodynamical admissibility was analyzed in [25]. In modern viscoelasticity, it is a central problem to understand how to model the memory kernels, and it should be argued as far as possible on physical grounds. So, the first question to answer to is why do we consider singular kernel models. More recently, new viscoelastic materials, such as viscoelastic gels, have been discovered and their mechanical properties are well described by virtue of convolution integral with singular kernels: for instance, fractional and hypergeometric kernels [1]. This applicative interest gave rise to a wide research activity concerning singular kernel problems, both in rigid thermodynamics with memory as well as in viscoelasticity (see, for instance, [2631], and especially concerning applications of fractional calculus to the theory of viscoelasticity and the study of new bio-inspired materials [15, 3235]. A recent book [36] provides an overview on this subject. In this framework, Fabrizio [37] analyzes the connection between Volterra and fractional derivatives models and shows how experimental results motivate us to adopt, as in this present article, less restrictive functional requirements on the kernel representing the relaxation modulus.

    4.1. Singular isothermal viscoelastic body with memory

    To start with, the one-dimensional classical viscoelasticity problem is recalled. It reads

    u(,0)=u0,ut(,0)=u1inΩ;  u=0onΣ=Ω×(0,T)E46

    where Ω = (0,1). When, to model the physical behavior of new materials or polymers, the regularity assumptions on the relaxation modulus are relaxed, Gis assumed to satisfy the following functional requirements

    GL1(0,T)C2(0,T), GL1(0,T), TE47

    that is, now, the relaxation function G(t) is not required to be finite at t= 0 and then Eq. (45) loses its meaning and, hence, needs to be replaced by a different one. The method to overcome this difficulty, devised in [28], consists in the introduction of a suitable sequence of regular problems, depending on a small parameter0 < ε≪ 1 which, in the limit ε→ 0 reduce to the singular problem under investigation. The key steps of the approximation strategy can be sketched as follows.

    • Let K, termed integrated relaxation function, denote

      it is well defined, since GL1(0,T),T+.

  • Then, introduce the regular problems:


  • together with the initial and boundary conditions


    • For each ε, the problem Pεis a regularapproximated problem since Gε(0) is finite and, therefore, the initial boundary value problem (49)-(50) admits a unique solution:

    • then, find approximated solutions uε, 0 < ε≪ 1,

    • show the existence of the limit solution u:=limε0uεu:=limε0uε

    • prove the uniqueness of the limit solution uwhich represents a weak solution admitted by the singularproblem.

    Note that, corresponding to each value of ε, the problem Pεis equivalent to the integralequation:

    Pε: uε(t)=0tKε(tτ)uxxε(τ)dτ+u1t+u0+0tdτ0τf(ξ)dξ,E51

    Partial derivation w.r. to t, twice, of Eq. (51) delivers Eq. (49) together with initial and boundary conditions (50). Furthermore, when ε= 0, we obtain the well-defined problem

    Pε: uε(t)=0tKε(tτ)uxxε(τ)dτ+u1t+u0+0tdτ0τf(ξ)dξ,E52

    where the superscripts, in the case ε= 0, are omitted for notational simplicity. Hence, the following theorems can be proved. Here only the outlines of the proofs are given; the details are comprised in [28] when homogeneous Dirichlet b.c.s (50) are imposed and in [27] when homogeneous Neumann b.c.s are considered.

    Theorem 1Given uε solution to the integral problem Pε (51), thenu(t)=limε0uε(t) inL2(Q), Q=Ω×(0,T).

    Proof’s outline:

    • weak formulation, on introduction of test functions φ∈ ϵ H1(Ω × (0,T) s.t. φx= 0, on ∂Ω,

    • consider separately the terms without ε,

    • the terms with uεand Kε,

    • prove convergence via Lebesgue’s theorem.

    Furthermore, the weak solution, as stated in the following theorem, is unique.

    Theorem 2The integral problem (52) admits a unique weak solution.

    Proof’s outline: The result is proved by contradiction, see [28] for details, assuming there are two different solution and, then, showing that such an assumption leads to a contradiction.

    As a final remark, we wish to emphasize that, since the isothermal rigid viscoelasticity model exhibits remarkable analogies, under the analytical point of view [38], with rigid thermodynamics with memory, then, analogous results can be obtained also in the study of singular kernel problems in such a framework [29].

    4.2. Magneto-viscoelasticity problems

    This section is concerned about a problem in magneto-viscoelasticity, again under the assumption of a memory kernel singular at the origin. The interest in magneto-viscoelastic material finds its motivation in the growing interest in new materials such as magneto-rheological elastomers or, in general, magneto-sensitive polymeric composites (see [3941] and references therein). The model adopted here to describe the magneto-elastic interaction is introduced in [42]. Evolution problems in magneto-elasticity are studied in [43] and, later magneto-viscoelasticity problems are considered in [44, 45]. Notably, under the analytical viewpoint, when the coupling with magnetization is considered, the problem to study is modeled via a nonlinear integro-differential system while the purely viscoelastic problem is linear.

    To understand the model equations, a brief introduction on the model magnetization here adopted, based on [46], who revisited the Gilbert magnetization model. Accordingly, when Ω ⊂ ℝ3 denotes the body configuration, the related magnetization changes according to the Landau Lifshitz equation, which, in Gilbert form, where mrepresents the magnetization vector reads

    γ1mtm×(aΔmmt)=0,|m|=1, γ,a+.E54

    The quantities of interest, in the general three-dimensional case, are the following ones:

    where the coefficients λklmnare subject to the condition


    Then, the following constitutive assumptions are assumed. Thus, the exchange magnetization energy is given by



    Then, the magneto-elastic energy is given by


    The viscoelastic energy is given by


    where the tensor’s entries of


    Then, the total energy of the system is given by


    taking into account, further to the single magnetic and viscoelastic contribution, of the exchange energy.

    4.2.1. A regular magneto-viscoelasticity problem

    The problem we are concerned about is the behavior of a viscoelastic body subject also to the presence of a magnetic field; in the one-dimensional case, it is modeled by the nonlinear system

    {uttG(0)uxx0tG(tτ)uxx(τ)dτλ2(Λ(m)m)x=f,mt+m|m|21ε+λΛ(m)uxmxx=0, inQE60

    where Ω = (0,1), Q:= Ω ×(0, T) and

    ≡ (0,m, where m= (m1, m2), denotes the magnetization vector, orthogonal to the conductor, since u≡ (u, 0, 0), when both quantities are written in ℝ3 in addition, vis the outer unit normal at the boundary ∂Ω, Λ is a linear operator defined by Λ(m) = (m1, m2) the scalar function uis the displacement in the direction of the conductor itself, here identified with the xaxis and λis a positive parameter. In addition, the term frepresents an external force which also includes the deformation history.

    In [44], the existence and uniqueness of the solution to the problem given by (60), together with the following initial and boundary conditions, is proved


    under the assumptions


    Then, the following existence and uniqueness result [44] holds.

    Theorem 3Given the problem (60)-(63), it admits a unique solution for any given T > 0 and ε small enough (i.e., ε<λ2G(T)), s.t.

    The proof, is based on the a priori estimate on the viscoelastic term:


    A result of existence, in a three-dimensional regular magneto-viscoelasticity problem, is given in [45].

    4.2.2. A singular magneto-viscoelasticity problem

    Now, as in the purely viscoelastic case, when the requirement G'L1(0,T)is removed, the magneto-viscoelasticity problem cannot be written under the form (60); however, since G'L1(0,T)via integration with respect to time of the integro-differential equation, it can be formulated in the following equivalent form


    The strategy to prove the existence result [47], relies on the fact that the classical problem (60) as soon as the initial time is t0 = ε, for any arbitrary ε > 0, the relaxation modulus satisfies the classical regularity requirements, namely, as in subSection 4.0.1, Gε():=G(ε+)implies that GεC2[0,T]Hence, each time-translatedapproximated problems


    with the assigned initial and boundary conditions


    is regular. Then, according to [44], the problem Pε admits a unique strong solution. According to [47], where all the needed proofs are given, the following existence result can be stated.

    Theorem 3.1For allT> 0, there exists a weak solution(u, m) to the problem(65)-(61)-(62), that is a vector function(u, m) s.t.

    which satisfies


    arbitrarily chosen test functions in Q.

    The proof, not included here, is provided in [47].

    Proof’s Outline:

    • consider the viscoelastic energy associated to the problem to obtain a suitable a priori estimate

    • consider the energy connected to interaction between magnetic and viscoelastic effects to obtain further suitable estimates

    • consider the total energy together with smooth enough initial data to estimate the energy at the generic time t

    • introduce an appropiate weak formulation and suitable test functions

    • consider separately the limit process when ε → 0

    As a closing remark, we can note that, under the applicative point of view as well as under the analytical one, the free energy associated to the model plays a crucial role. Indeed, the proof relies on estimates which are based on the free energies connected to the model here adopted. Specifically, the viscoelastic energy allows [47], also in the magneto-viscoelastic case, to prove an a priori estimate on which the subsequent results are based. This is not surprising since the connection relating free energies and evolution problems is well known; see for instance [48] and references therein.


    S. Carillo wishes to acknowledge the partial financial support of GNFM-INDAM, INFN, and SAPIENZA Università di Roma.


    • It is easily seen that (B3) always holds, for instance, when k is eventually concave down as t → ∞.

    © 2016 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

    How to cite and reference

    Link to this chapter Copy to clipboard

    Cite this chapter Copy to clipboard

    Sandra Carillo and Claudio Giorgi (September 21st 2016). Non-Classical Memory Kernels in Linear Viscoelasticity, Viscoelastic and Viscoplastic Materials, Mohamed Fathy El-Amin, IntechOpen, DOI: 10.5772/64251. Available from:

    chapter statistics

    1326total chapter downloads

    5Crossref citations

    More statistics for editors and authors

    Login to your personal dashboard for more detailed statistics on your publications.

    Access personal reporting

    Related Content

    This Book

    Next chapter

    Finite Element Modeling and Experiments of Systems with Viscoelastic Materials for Vibration Attenuation

    By Antonio Marcos G. de Lima, Luiz Fernando F. Rodovalho and Romes A. Borges

    Related Book

    First chapter

    Mass-Transfer Characteristics of a Double-Y-Type Microchannel Device

    By Susumu Nii

    We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

    More About Us