Abstract
The study addresses the physical background and modeling of linear viscoelastic response functions and their reasonable relationships to the Caputo-Fabrizio fractional operator via the Prony (Dirichlet series) series decomposition. The problem of interconversion with power-law and exponential (single and multi-term functions) has been discussed. Special attentions have been paid on the Prony series decomposition approach, the related interconversion problems and the expression of the viscoelastic constitutive equations in terms of Caputo-Fabrizio fractional operator.
Nothing is done in a vacuum; we must all stand on our forefathers to better ourselves and the world around us.
(Sir Isaac Newton, in his letter to rival Robert Hooke in 1676)
Science is built up of facts, as a house is built of stones; but an accumulation of facts is no more a science than a heap of stones is a house.
(Henri Poincare)
1. Introduction
This article addresses the physical background of modeling of dissipative phenomena, precisely response functions such as stress-strain relationships in the framework of the linear viscoelasticity and their reasonable relationships to the Caputo-Fabrizio fractional operator [1] via the Prony (Dirichlet series) series decomposition [2].
The appearance of the new definition of fractional derivatives with non-singular kernels was provoked by needs to model dissipative transport processes in many new materials appearing in modern technologies [1]. Recently the main achievements, especially results related to diffusion problems, were analyzed Hristov in [3] and we will avoid the thorough browsing and comments of published results (see also the rich list of reference in Hristov [3]). In the context of diffusion problems it was demonstrated that the Caputo-Fabrizio operator appears naturally in diffusion models [3, 4] when the flux-gradient relationship is expressed by a Jeffrey relaxation kernel and the Maxwell-Cattaneo concept of flux.
At the same time the Caputo-Fabrizio fractional operator was criticized [5–7] with points of view based on the classical fractional calculus with singular power-law memory kernel and with examples from the signal processing [6] and to some extent touching formal rheological relationship [5]. Despite the mathematical exactness of the counterexamples they do not focus the attention on the physical basis leading to exponential memory kernels. As a matter of fact, we need a clear answer what really this new operator models and how it appears in the modeling of physical problems, despite the fact that some properties known from the classical fractional calculus, such as the index law [6] , are not satisfied.
1.1. Motivation of This Study
The analysis done here and the principle task are oriented to modeling of viscoelastic constitutive relationships in terms of Caputo-Fabrizio operator naturally appearing through Prony series decomposition [2] of stress relaxation functions and not obeying power-law behaviors. That is, we focus on viscoelastic materials with dynamics which cannot be modeled with the classical fractional derivatives of Riemann-Liouville and Caputo [8].
This introduction avoids huge citations of works on Caputo-Fabrizio operators since without understanding of the physical background and the logic of its appearance any analysis of models created by formalistic fractionalization (see the comments in Hristov et al. [3]) is unproductive. Moreover, to the author personal experience as editor many manuscripts devoted to applications of the Caputo-Fabrizio derivative in formalistically fractionalized existing models are directly rejected by the reviewers with motivations based of the opinions in the criticizing articles [5–7]. This situation resembles that in the Catch 22 movie without perspective for escape. The existing situation is a consequence of some main reasons: (1) The Caputo-Fabrizio operator does not hold some properties such as the semi-groups, which are existing with the classical fractional derivatives with singular kernels, the strange form of the associated fractional integral and these issues cast doubts when it is applied inadequately (in blind manner) to various functions in a manner known from the integer-order calculus. (2) The formalistic approach by simple replacements of integer or fractional order (with singular kernels) derivatives with the Caputo-Fabrizio operator in existing models without taking into account the physics behind. (3) Last but not least, the human factor of author's rivaling which is dividing the scientific society into competing groups rejecting the achievements of each other. All these elements of the current situation create a discouraging atmosphere and disbelieve that this new fractional operator really cannot model natural phenomena and actually stops the further research, generally among the young researchers highly aspiring publications of submitted manuscripts.
The deep physics behind the fractional operator with exponential kernel (and the author's long time experience in science) motivated this study and the efforts are oriented to show that the existing knowledge and models, as well as techniques of data treatment, in the framework of linear viscoelasticity, lead naturally to formulation of the Caputo-Fabrizio fractional operator. This is in the context of the Sir Isaac Newton quote at the beginning of the article: the steps ahead on the shoulder of existing facts and results on the road to creations of new information are natural ways and actually the exciting moments in the beautiful journey in the world of science.
1.2. Aim and Paper Organization
This article is organized as follows: section 3 presents briefly the main properties of the Caputo-Fabrizio operator that will be used further in the analysis of the viscoelastic problems. Section 3 addresses the constitutive equations of viscoelasticity based on the Boltzmann superposition principle [9] and fading memory approach incorporated the hereditary stress and strain integrals. The principle problems of the quality and choice of the response function are discussed and the main properties requited are outlined in section 3.3. The interconversion of relaxation and creep response functions is discussed in section 4: the linear and non-linear scale-invariant (power-law) (related to application of Riemann-Liouville and responses are analyzed. Section 5 focuses on decompositions of experimental response functions by Prony series leading to discrete spectra and the related interconversions by examples with single term, two-terms and multi-term responses of exponential type, as well as interconversion of relaxation and creep compliance expressed as Prony series. Section 6 demonstrates how the constitutive equations (based on the Maxwell model) can be expressed in terms of the Caputo-Fabrizio fractional operator. The subsection 6.3 demonstrates that approximation of the response function by Bessel functions (of Maxwell-like materials) of first kind and expressed as infinite Dirichlet series naturally leads to incorporation of the Caputo-Fabrizio operator in the constitutive equations. Section 7 demodulates briefly how constitutive relationship (following the idea of Bagley and Torvik) with two fractional operators of Caputo-Fabrizio (of different orders) can be formulated.
2. Caputo -Fabrizio Operator
2.1. Definition
The Caputo -Fabrizio operator is defined as [1]
In Equation (1) M(α) is a normalization function such that M(0) = M(1) = 1. This definition is of Caputo-type because there is a convolution of the derivative df(t)/dt. The explanations in [1] relate the development of Equation (1) to the classical Caputo derivative [8] by mechanistic replacements of the singular kernel (by a non-singular exponential kernel) and the normalization function (see the explanations in Hristov [3].
From Equation (1) it follows that if f(t) = C = const., then as in the classical Caputo derivative [8]. Actually, Equation (1) is a convolution of f(t) and the convolution operator K [1]
The integration by parts in Equation (1) results in an alternative form Caputo et al. [10]
Futher, the integer order differentiation of follows the rule [1]
The associated fractional integral is Hristov [3] and Losada and Nieto [11]
With the assumption that M(α) = 1 used by Caputo and Fabrizio [1, 10]. The second definition of Losada and Nieto [11] is
Hereafter we will use the definition (Equation 6).
2.2. Laplace Transform
The Laplace transform of with a = 0 has the following Laplace transform LT given with p variable [1] taking into account the general rule of Laplace transform of a convolution, namely
Hence,
2.3. Fractional Derivative of Elementary and Transcendental Functions
Linear function f(t) = Ct [1]
Power-Law function f(t) = Ctβ [3, 12].
For f(t) = Ctβ and β > 0 the fractional derivative is
For α = 1 the expression (Equation 10) reduces to the classical (integer-order) result Cβtβ−1.
Exponential function f(t) = exp(βt) [3, 12]
For α → 1 the second exponential term in the nominator of Equation (11) goes to zero and therefore .
For β < 0 we have
For α → 1 the second exponential term in the nominator of Equation (12) goes to zero and therefore .
2.4. Caputo-Fabrizio Fractional Operator: Determination of the Fractional Parameter
In the Caputo-Fabrizio operator there is formal ambiguity because the stretched time is multiplied by a dimensional factor α/(1 − α) which should have a dimension s−1. This contradicts the definition of the fractional order (parameter) because physically α is dimensionless. The answer given in [3, 13] resolved the problem by nondimesionalization of the exponential function by help of characteristic time scale of the relaxation process t0 (the maximum time of the experiment in the sense of the rheological tests discussed here) , namely
The nondimensalization does not change the meaning of the exponential relaxation function but avoid any doubts about the definition of the fractional order α as [3, 13].
and relates it to data that can be really recovered from experimental data, such the relaxation and retardation times (see the sequel)
The relationship (Equation 14) says that for τ/t0 = 1 we get α = 1/2. Further, depending on the ratio τ/t0 we may have fractional orders roughly arranged in two groups [14]: a) when 0 ≤ τ/t0 ≤ 1 we have fractional orders α ∈ [0.5, 1.0] and the relaxation time τ is shorter than the macroscopic process observation time t0, and b) 1 ≤ τ/t0 < ∞ the fractional orders are α∈(0, 0.5] since the relaxation time τ is larger than the macroscopic process observation time t0. Qualitatively, for τ/t0 < 1 the relaxations could be considered as fast (rapid) relaxations, while τ/t0>1 is related to slow relaxations (see in Hristov [14] the comments of numerical values of Prony decomposition and relaxation times)
Last to this point, but not least, as it was commented in Hristov [14], the ratio τ/t0 is not integer and expressing the memory kernel as where we get a fractional operator. Therefore, the memory kernel of the Caputo-Fabrizio operator is controlled by a non-integer parameter in the context of what is needed to say that this operator is fractional, despite the fact that it does not repeat exactly the properties of the Classical Riemann-Liouville and Caputo derivatives.
3. Constitutive Equations of Viscoelasticity: Fading Memory Approach
3.1. Boltzmann Superposition Principle and Fading Memory Concept
The fading memory concept relating the flux to its gradient, for simple materials [15–17], is expressed by the following relation relating the flux and the gradient , namely
This definition is actually the Boltzmann linear superposition functional (Equation 16))
relating the present state of the flux to its history [9, 16–18] through the influence function (memory kernel) R(t, z) during the time interval defined by τ. In Equation (16) m and λ are transport coefficients (diffusivities) with real physical meanings as it will be demonstrated in the sequel.
The memory function could be unbounded and scale invariant such as R(t, τ) = t−μ with integration singularity at or bounded and not scale invariant such as R(t, τ) = e−t/τ (see the analysis in Hristov [14] and the comments further in this article).
3.2. Stress-Strain Viscoelasticity Response and Hereditary Integral Construction and Response Functions
The linear theory of elasticity [19] (chapter 1) considers the stress in a sheared solid body as a quantity proportional to the shear, while in the liquids the shearing stresses are proportional to the rate of shear . Most solid materials, for example polymers, compromising both effects are called viscoelastic. When a slab of solid material under a shearing motion caused by a step change in the stress load applied to it (Heaviside unit step function H0(t)) exhibits a strain (in one dimension) [19] (chapter 1)
With perfect elastic behavior of the body [19] σ(t) = σ0H0(t) for t > 0. On the contrary, in an ideal viscous fluid the stress is infinite and for t > 0 and the strain is ε(t) = (σ0/η)t, thus introducing the coefficient of viscosity η. Real materials do not shear with infinite speeds that is the reason of the concept of a finite relaxation time τ [19]. Precisely, in solids the stresses attain finite values for long times. In contrast, in viscous fluids the stresses approach zero.
The task of this study addresses the functional representations of the viscoelastic material responses, that is : 1) the stress relaxation function R(k, t), that is the stress history due to a shear step of size ε, and 2) the creep function C(t) (shear history) due to unit stress σ applied. In the linear viscoelastic theory [19] the responses can be approximated as R(ε, t) = G(t)ε + O(ε3) and C(σ, t) = J(t)σ + O(σ3), thus defining the linear stress relaxation modulus G(t) and the linear creep compliance J(t). The common mechanistic models explaining the behavior of viscoelastic materials utilize linear springs and dashpots coming from the Maxwell interpretation [20–22] and modeling the pure elastic behavior and the pure viscous state, correspondingly.
Superposition of single-step material responses results in functional relationships of stress and strain in the linear viscoelasticity concept [19] incorporating the a time lag in G(t) and J(t) through hereditary stress (Equation 18) and creep integrals (Equation 19), namely
Both σ(t) a and ε(t) are causal functions and therefore the lower terminals in Equations (18, 19) are set at t = 0.
Applying the fading memory concept [9, 20] it is possible to relate the instantaneous responses G∞ and J∞ corresponding to equilibrium states (for long times) when the effects of memory terms (convolution integrals) fade out, namely
3.3. Response Function: General Properties and Requirements
For a linear isotropic viscoelastic body exerting uniaxial stress the Boltzmann superposition principle formulates the constitutive equation relating the strain responses by help of the following hereditary integral [23]
The relaxation function G(t) should decrease monotonically and account for short and long-time strains, thus the condition (Equation 23) should be obeyed [24–26]
Coleman and Noll, in their article on foundation principles of linear viscoelasticity [27], point out that there is no universal approach in definition of unique relaxation function R(s) [the fading memory R(s) in Equation (15) its is equivalent to G(s) in Equation (22) . Despite this, two principle features are required to be obeyed by R(s):
R(s) is defined for 0 ≤ s < ∞ and R(s) > 0, and
R(s) decays monotonically to zero for large s , that is .
Boltzmann (1874) in [9] (see also [14, 27]) suggested two memory functions widely applied so far, namely
a) Generalized power-law memory function R(s) = sr(s+1)−μ of order r when r < μ. For r = 0 the memory function R(s) is unbounded (singular) at t0+ and this scale invariant kernel forms the constructions of the classical fractional integrals and derivatives [8].
b) Ordinary memory kernel , bounded at t0+ as exponential function e−βs , which is not invariant with respect to the time scale. The definition of Caputo-Fabrizio fractional operator utilizes this memory kernel.
The adequate mathematical structure in the construction of the relaxation (memory) function is the main problem that should be resolved in the modeling of viscoelastic responses. Generally, the relaxation function should adequately describe the natural process and therefore its structure should be tested (defined)by a data conversion algorithm (by fitting experimental data). Despite this intuitive and to greater extent logical approach the response function should satisfy some requirement defined in Garbarski [28] and Winter [29], and summarized in Table 1; with some author's comments related to application of fractional operators in the viscoelastic constitutive equations. Further, in the response of viscoelastic material to strain excitation, the stress relaxation spectrum (relaxation modulus) provides complete information concerning the time-dependent part of the material response. In the opposite situation, with materials undergoing stress excitations the retardation spectra (the compliances) [30] provide the required information. Hence, if the spectra are known it is possible to calculate the response of any excitation [20].
Table 1
| Properties of the Response function [28] | Properties of the data conversion algorithm [29] | Comments (present work) | |
|---|---|---|---|
| 1 | Be as simple as possible due to consequent practical applications | Good fit of the experimental data | |
| 2 | Should be positive and monotonically decreasing since non-monotonically decreasing functions have no physical meanings. | Avoidance of overfitting, that is the algorithm should be able to find optimum amount of details (i.e., parameter) | |
| 3 | After integration the function should be convergent to a limited value at infinite time | The format of the relaxation function should not be predetermined. It has to be freely adjustable during the data fitting | This is a general comment but working with fractional derivatives we have a limitation of kernels (response functions) that can be use |
| 4 | To have a simple Laplace transform | The resulting material parameter should have physical meaning | |
| 5 | Flexible to be adjusted to experimental data taken from relaxation tests. | Minimization of truncation error | |
| 6 | To allow calculations of consequent parameters of the viscoelastic relaxations such as spectra, moduli, etc. | Checking the of experimental data quality | This important since the parameter determinations are ill-posed problems |
| 7 | Possibility to be tabulated if analytical representation is not possible | For practical use R(s) should be expressed by a function or a sum which can be easy integrable in various viscoelastic models and calculations | |
| 8 | The integral of the spectrum should be convergent to a finite value | For practical reasons the number of parameters should small as much as possible | This point is very important in the construction of relaxation kernels of fractional operators |
Response function desired properties.
Besides, the thermodynamic theory of linear viscoelastic materials was developed extensively in the sixties of the last century [31]. In this context, the second law of thermodynamics yields severe restrictions on the constitutive properties [32–34]. This problem is beyond the scope of the present analysis but we may say that when the constitutive equation (Equation 20) contains exponential kernel it is compatible to second law of thermodynamics [34].
The principle relaxation functions used and the applicability of the aforementioned requirements will be analyzed next.
4. Interconversion of Relaxation and Creep: Problem and Power-law Responses
4.1. Interconversion Problem
In the linear viscoelastic materials the relationships between the creep compliance J(t) and the relaxation response G(t) are expressed as convolution integrals in accordance with the Boltzmann superposition principle [9] (see Equations 20, 21)
With the Laplace transforms of Equations (20, 21) we get equivalent relationships in the p-space (p is the transform variable)
Hence, we may recast Equation (24) as
Equating the right-hand sides of Equation (25) we get an implicit fundamental relationship
with the constraints that
Explicit form of Equation (26) can be obtained if the analytical forms of either G(t) orJ(t) is known. Examples of interconversion related to power-law and exponential memory kernels are discussed in the sequel.
4.2. Scale-Invariant (Power-Law) Response
4.2.1. Physical Preliminaries Related to the Power-Law Response
Relatively short-time relaxation history modeled by time-dependent power-law can be observed in the time evolution of G(t) from zero to the end of the time of observation t0 [19]: commonly at the beginning of the relaxation process (short times) and along the long tail (long times). The power law is scale-invariant. Materials exhibiting such behavior are called power-law viscoelastic materials [35–42] and can be easily detected by the linear behaviors of J(t) ~ tp and G(t) ~ t−p in logarithmic scales [43]. Actually, we have to memorize that this type of relaxation was modeled in the article of Caputo (1967) where the construction of the Caputo derivative was conceived [44]. The power-law functional relationships of the relaxation and the compliance, actually are the reasons to recognize the integral and derivatives of the classical fractional calculus [8] as adequate modeling tools since they are based on the same type of kernels [39, 42, 45], known also as weak singular kernels. As commented by Tschoegl [20] the scaling relationship G(t) ~ t−p is a fractional version of the Trouton law . Here, we have to mention that Metzler et al. [39] discussing relaxations of filled /polymers stated that non-exponential (non-Debay) relaxation implies memory (Sic!). We may consider this declarative opinion (the results of Glockle and Nonnemacher [46]) as adequate to the situation in the 90s of the last century, when fractional calculus was only related to power-law kernels [47].
Moreover, as a valuable analysis showing when the power-law is the adequate response function we refer to Chapter 5 of Findley et al. [48] where it is clearly demonstrated that for many materials the power-law response function tμ with μ < 0.5 [48] is a good short-time approximation for materials such as plastics, metals and concrete. For short-time loading the creep of many different rigid plastics with sufficient accuracy can be presented as linear power-law , where μ is stress-independent and nearly temperature-independent, too. The same approach following from the linear approximation approach of Pipkin [19] is commented in 3.2 . In this case the corresponding retardation spectrum is Findley et al. [48]
The function (Equation 28) may be considered as continuum spectrum of material retardation times proportional to t−α . From this position, the step from Equation (28) toward modeling of the relaxation processes by the classical Riemann-Liouville derivative is straightforward. In this context, Bagley and Torvik [49] using (Equation 28) demonstrated that following constitutive equation holds
After application of the Leibniz rule to Equation (29) the result is
For ε(0) = 0, since no strain at t = 0 exists we get
Then, for G1 = G0 = const. the stress relaxation modulus G(t) and the creep compliance J(t) are
4.2.2. Example 1. Interconversion of Power-Law Relaxation and Creep: Linear Case
Here, we consider again the power-law response in the context of the response function interconversion and in the sense of the relationships (Equations 26, 27). With G(t) = At−μ , where A is a data fitting coefficient and 0 < μ < 1 , the Laplace transform in Equation (26) provides a power-law creep (Equation 33) compliance [50]
Following Hanyga, neither experimental nor theoretical reasons lead to the assumption that the creep rate function and the stress relaxation function are bounded and regular at t = 0 [26]. The same point of view is valid for the unbounded kernels singular at t = 0 [14]. In the seminal study of Boltzmann [9] the value of α = 1 was suggested. The condition 0 < α < 1 comes from the fact that for α > 1 we get an infinite propagation speed [26]. Actually, the idea to apply fractional calculus, as is demonstrated by the model developed in the foundation studies of Scott-Blair [51–53] and Bagley and Torvik [49] comes from the experimental (empirical) findings of Nutting [54, 55]: precisely, for some viscoelastic materials the power-law relationship [ε(t)/Const] ~ σntα is satisfied. Consequently, we may obtain (Equation 32).
Since the power-law relaxation is not the main problem discussed in this work, at the end of this point, we refer to the study of Carillo and Giorgi [56] (section 4 in this chapter) where it is clearly demonstrated that the response of the material may be for short-time or as long tail (for long times) modeled by time-dependent power-law function. Here, we may repeat again (see the comments in Hristov [14], too) that, if the power-law is not exhibited by the material response functions then it is inadequate to apply the power-law memory kernels.
4.2.3. Example 2.Interconversion of Power-Law Relaxation and Creep: Non-linear Case
If the material exhibits a non-linear viscoelastic behavior, for instance, as commented by Lakes et al. [50], that is
The relaxation function G(t, ε) can be expressed as a sum of products, precisely sum of products of functions of time and functions of strain, that is Lakes and Vanderby [50] G(t, ε) = Gt(t)g(ε) . Assuming no interaction between the steps (immediate and delayed Heaviside steps) strains in the summation series and may write [50]
Consequently, the creep compliance is
and for the linear case we will obtain the relationship (Equation 26).
To have explicit form of these relationships, if we assume a power-law approximation in time [50] J(σ, t) = A(σ)tμ, that is , then we get
The nonlinear material exhibits a relaxation response which contains a sum of power-law terms, as given in equation (Equation 38) (see more comments in Lakes and Vanderby [50]).
Hereafter, non-linear behavior generally related to temperature-dependent and aging viscoelastic materials will not be discussed since the principle task of the present analysis is to show the correct origin and the physical background leading to the formulation of the Caputo-Fabrizio fractional operator.
5. Interconversion of Non Power-law Relaxation and Creep
5.1. Non Power-Law Response: Relaxation Curve Approximation
The selection of approximation function which would fit the experimental points depends on the type of materials and would be established by simple trial error process. In this context, two principle issues arise when experimental data for linear viscoelastic materials should be treated [57]: (1) Appropriate relaxation curve approximation and (2) interconversion problem. The experimental data are commonly taken in the time-domain (or frequency-domain) tests. Now, we stress the attention on discrete spectrum approximation by Prony series [2] where the condition for the monotonicity of the functions (see the sequel) is satisfied.
5.2. Discrete Relaxation Spectrum by Prony Series Decompositions and Relevant Interconversions
Unlike the linear elasticity, where the material functions are related algebraically, the relationships in the linear viscoelasticity are time-dependent, as already was demonstrated in the previous section of the article. One known method is to find the unknown functions G(t) and J(t) by fitting the data points by a known function. However, when the focus is on the adequate construction of a fractional operator with a memory kernel obeying all mandatory conditions imposed on it, then the choice of approximating function is highly restricted. If the power-law is not the adequate choice then a series approximation is the more suitable solution. Since we address the construction of the Caputo-Fabrizio operator, thus it is natural to tackle data fitting by series of exponential terms (matching the function of the desired memory kernel) known as Prony series (also as Dirichlet series [58]).
5.2.1. Prony Series Decompositions
The viscoelastic relaxation function can be expressed by a discrete relaxation spectrum through a decomposition as a Prony series gP(t)with Nϕ terms [25, 59–64] with rate constants βi, namely
or through normalized weights (amplitudes or normalized relaxation moduli) λi [21, 59] as
In Equations (39, 40) the parameters g∞ and gi are equilibrium (at large times) values and the relaxation moduli (stiffness), respectively, are constrained according to Brinson and Brinson [21]
The Prony series components have spectral strength gi and relaxation time τi. This is the so-called discrete Prony series representation known also as discrete relaxation spectrum [20, 65–69]. Besides, the series approximation (Equation 39) may be substituted in formulation of the materials law such as Maxwell, Kelvin Voigt, etc. The relaxation time τi associated with the ith element is related to the characteristic time of the spring-damper (dashpot) element and can be defined as the ratio of the viscosity over the elastic modulus, that is τi = gi/λi [68].
The first derivative of the Prony series (Equation 42) for t = 0 is finite, a fact irrespective of the number of terms used in the approximation.
As commented by Bradshaw and Brinson [57], an exact solution of interconversion problem is possible if the known function is approximated by Prony series where all coefficients are positive [using the forms (Equation 39, 40) ] and the basic interconversion relationship (Equation 26) is satisfied. In this case, by applying Laplace transform solutions, the resulting Prony series are analytically exact [70–73].
The Prony series approximations lead to the generalized Maxwell viscoelastic body (known also as Maxwell-Wiechert model) with Nϕ spring-dashpot elements in parallel. The Prony series decompositions are applicable to any viscoelastic models through their time-dependent shear and bulk moduli [74–76]. For a deep thermodynamic analysis of such type of materials termed viscoelastic solids of exponential type (VESET) we refer to works of Fabrizio et al. [32, 33].
A common step in approximation by the truncated exponential sums is the definition of preliminary stipulated decay rates βi , in a logarithmic scale that is , for Nϕ = 4 , for instance β1…4 = 100, 10, 1, 0.1 [77]. This approach, is useful because the corresponding fractional parameters in the αi in this case can be easily calculated (see further Equation 65) as α1…4 = 0.009, 0.09, 0.5, 0.9, respectively. Alternatively, fixing the points (in normal or logarithmic scale) ti along the time axis we may define the corresponding βi = 1/τi.
The parameter estimation is important and the first step was done in the seminal work of Prony (1795) [2] and several algorithms have been developed among them: graphically by log−log plots [78, 79], least squares method [80–82], nonlinear optimization methods [83], genetic based algorithm [84], from the continuous relaxation time spectrum [68, 72], logarithm equidistant distribution of relaxation times (known as R- method) [68, 77, 85], quasi-linearization for multi-exponential decay curves fitting [86], etc.
The number of the exponential terms in Equation (39) depends on the accuracy of data fitting. Commonly fourth-order Prony series fit adequately the stress-relaxation data in cases of non-linear viscoelastic behaviors [62, 87–93], while long-term relaxation tests need 10-15 terms [76, 85, 94–100]. In general, the problem corresponds to identifications of the kernel of an integral Fredholm equation which, actually, is ill-posed problem [68, 101–103] requiring Tikhonov regularization [104]. Comments on the required terms in the Prony decomposition are available in [14].
5.3. Examples of Interconversions With Exponential Terms
5.3.1. Example 3. Interconversion of a Single Exponential Model [105]
For a single exponential model (the Maxwell model) the relaxation and creep functions are
The condition (Equation 27) is automatically satisfied, while for t = 0 in (Equation 43) (and taking into account ;Equation 27) the results is
The Laplace transforms of (Equations 43, 44) lead to an equation that should be solved
Since (Equation 46) should be satisfied for all p it follows that τ1 = j1λ1 and this condition with (Equation 46) leads to Anderssen et al. [105]
5.3.1.1. Interconversion from relaxation to creep
5.3.1.2. Interconversion from creep to relaxation
More details related to the sensitivities and relative errors of these two interconversions are available in Anderssen et al. [105]
5.3.2. Example 4. Interconversion of a Double Exponential Model [105]
In this case the relaxation and creep functions are
The values of the relaxation times τi are zeros of
and consequently
5.3.3. Example 5. Interconversion of Multi-Term Exponential Model [105]
By Prony decomposition the relaxation and creep functions can be generally presented through discrete spectra as Anderssen et al. [105]
The inverse times, that is the relaxation times τi and the retardation times λi satisfy
The advantage of this (Prony decomposition) approach is that the monotonicity is automatically satisfied [105]
The attempts to solve (Equation 26) address many approaches for either J(t) or G(t) depending on the experimental data available [30, 106–109]
Now we recall that the Caputo-Fabrizio fractional operator [1, 3, 4, 110, 10] uses an exponential kernel matching the basic element of the Prony series.
5.3.4. Example 6. Interconversion of Prony Series of Relaxation and Creep
With Prony decomposition of the relaxation spectrum, the next step is to determine the creep modulus (compliance) J(t). The interconversion equation (Equation 56) simple gives [22, 70, 71, 105, 111, 112].
As a rule in the rheological studies, the discrete exponential models for G(t) and J(t), corresponding to associated relaxation and creep spectra H(τ) and L(τ) [105] (as sums of delta functions ) are used that assures the monotonic behavior automatically and we have
The constraint imposed [70, 105] by (Equations 26, 27) (see also Equation 56) allows H(τ) to be defined as function of L(τ) and vice versa. We will skip this general problem solved in [105] and to some extent demonstrated by Example 3 and 4, and will focus the attention on the discrete approximation of the relaxation (and compliance) spectra by Prony series.
The Prony series approximation of the stress relaxation (in a normalized form) is
The corresponding form of J(t) is
The inverse times, that is the relaxation times τi and the retardation times λi satisfy the inequalities (Equation 55). The attempts to solve (Equation 56) address many approaches for either J(t) or G(t) depending on the type of experimental data available [30, 106–109].
6. Caputo-Fabrizio Operator in the Constitutive Viscoelastic Equations
6.1. Relaxation Function in Terms of Caputo-Fabrizio Operator
Thus, applying the Prony approximation of the relaxation curve and substituting (Equation 61) in the convolution integral (Equation 22) the following approximation is obtained [14, 113]
Since we operate with a Prony series, that is with a finite sum of exponentials, then the inversion of the summation and the integral yields [113]
Thus, the memory effect from the convolution integral can be easy incorporated in each term of the Prony series [99],namely
The finite sum of N+1 terms σi(t) leads directly to the generalized Maxwell model, that is at any time t we have
σi(t) is a product of the spring modulus Ei and its current strain ki(t), which to some extent could be considered as a hidden material variable. The clear physical meaning of this result is that the strain ki(t) at a given time t is expressed as convolution integral with an exponential kernel [14]. Therefore, the model parameters that should be defined, following the approximation of the right-hand side of (Equation 62), are [14, 113] : the single separate spring stiffness E0 and the spring stiffness Ei as well as the relaxations time τi of each i−th Maxwell element.
6.2. Relationships of the Relaxation Time Spectrum and Fractional Order Spectrum
As it was mentioned at the beginning, the fractional parameter α is related to the dimensionless relation time as α = 1/(1 − τ/t0). Since we have a spectrum of relaxation times τi, then the spectrum of the fractional orders (parameters) is
Therefore, ki(t) can be expressed in a form related to the construction of the Caputo-Fabrizio operator, namely
6.2.1. Stress Relaxation in Terms of Caputo-Fabrizio Operator
Thus, the constitutive equation of the stress relaxation can be presented as Hristov [14]
Now, we turn on the determination of the spectrum of fractional orders αi = τi/t0 . From experimental data fittings there are a limited number of numerical values of relaxation times τi. Now, the principle problem at issue is the determination of the characteristic time t0. Since the experiments last limited times then we may assume that t0 equal the elapsed time te of the experiments. The literature data concerning Prony decompositions reveal (see the analysis in [14]) that the relaxation times form two groups (see the comments at the begging when the fractional order determination of the Caputo-Fabrizio operators was commented):: (1) relaxation times less then the elapsed time τi < te ⇒ τi/te ≤ 1 and (2) relaxation time greater than the elapsed time τi > te ⇒ τi/te ≥ 1 .
Consequently, addressing the fractional orders αi we have: i) fast relaxations for αi ∈ [0.5−1) corresponding to τi/te ≤ 1 and ii) slow relaxations for αi ∈ (0−0.5] when τi/te ≥ 1. Numerical examples supporting this estimates are reported in Hristov [14]
6.2.2. Creep Compliance in Terms of Caputo-Fabrizio Operator
By analogy of the stress relaxation expression in terms of Caputo-Fabrizio derivative, we may transform the Prony series decomposed compliance (Equation 54) as
where are the scaled (dimensionless) retardation times.
Hence, with the construction of the convolution integral describing the strain history we have
Thus, the strain history is
Then, the strain can be expressed as
where
Therefore,(Equation 73) defines a Caputo-Fabrizio operators with respect to σ(t) and the fractional order βi is related to the scaled retardation time ( from the spectrum) as
6.3. Description of Maxwell-Type Viscoelastic Media With Material Responses Modeled by Bessel Functions in Terms of Caputo-Fabrizio Operators
In the last two years (2016-2017) an interesting approach was developed by the group around Professor Mainardi [114–119] which could be considered as attempts to generalize the relaxation functions in the linear viscoelastic models, an approach also investigated in Colombaro et al. [120] and Guisti and Colombaro [121]. The main idea comes from the possibility to represent the relaxation function in a viscoelastic Maxwell-type body by infinite discrete spectrum with times related to the zeros of Bessel functions of the first kind [118, 119]. Here we will present some key points of these studies,in the context of the ideas developed in this article, that is to show the incorporation of the Caputo-Fabrizio operators in the relaxation functions expressed as infinite Dirichlet series (or infinite Prony series as the authors of these studies defined them).
The analysis developed in Colombaro et al. [114], Colombaro et al. [115], Colombaro and Guisti [116], Colombaro et al. [117], Guisti and Mainardi [118], and Guisti and Mainardi [119] is based on the power series representation of the modified Bessel function of the fist kind as Abramowitz and Stegun [122]
with an asymptotic representation
Now, if the function Fν(t) has a Laplace transform [118]
It can be presented in the time domain as Guisti and Mainardi [118]
which is an expression as Dirichlet series, locally integrable, positive and increasing function for t > 0 [118], that is
because for s → ∞ from Equation (77) we have
If the is considered as a part of the Laplace transformation of the relaxation memory function Φν(t) ( here we use the original notations as that used in Guisti and Mainardi [118]), that is , we have
In the time domain the relaxation function Φν(t) is [recall equation (Equation 78)]
Now, let us turn on to the interconversion problem. From as it is suggested [118] that
In the time domain the relaxation modulus G(t) is
Skipping details in calculation, the creep memory function denoted as in the Laplace domain is Colombaro et al. [114], Colombaro and Guisti [116], and Guisti and Mainardi [119]
and
which follows from the interconversion relationships.
For further deep reading in this elegant mathematical studies we refer to Colombaro et al. [114], Colombaro and Guisti [116], Colombaro et al. [117], Guisti and Mainardi [118], and Guisti and Mainardi [119] where it is was demonstrated that the asymptotic behaviors of the viscoelastic responses in the so-called Bessel medium are (precisely in Colombaro et al. [114])
which actually follows from Equation (79).
The main idea to use a relaxation function expressed by modified Bessel functions of first kind comes from the possibility to calculate the sum of the infinite series of reciprocal positive zeros of the Bessel functions Jν [118], namely
With (Equation 89) the relaxation modulus becomes [118]
Hence, the relaxation modulus is represented as an infinite Dirichlet series which are absolutely converging, which actually represent a discrete spectrum where the stiffness of the nth element (spring-dashpot) is
(the prefix B means Bessel) and the relaxation times are . That is, in the context of the analysis done in this study, we have
Therefore, we get an expression in terms of an infinite sum of exponential functions, as in the Prony approximations, but now the relaxation times are defined by the zeros of the Bessel function. If a time scale t0 of modeled relaxation process exists, then nth dimensionless relaxation time is 0 < τn/t0 ≤ 1, as it was demonstrated several times with the Prony series. Now, let us consider that τn/t0 = (1 − βn)/βn ,where 0 < β < 1 and then we may re-write (Equation 91) as
Then, the construction the convolution integral of the stress relaxation with G(t) yields
Interchanging the orders of the summation and the integral in equation (Equation 93) , we get
Further, the step toward incorporation of the Caputo-Fabrizio operator is straightforward, namely
or in a compact form as
with the condition
Here denotes a Caputo-Fabrizio fractional operator with a fractional order βn based on the positive zeros of Jv.
Now, the naturally question coming to mind is: How this result can be applied to real data related to the stress relaxation and creep compliance of real viscoelastic materials? Unfortunately, no answers to this question exist in all works dealing with so-called Bessel media [114–119]. The first problem immediately appearing is: how the relaxation times can be related to the real data? As mentioned in several points of this article, the relaxation times corresponds to measurements taken at equidistantly distributed (in normal scale or logarithmic scale) points along the time axis. How, this real approach could be related to the zeros of the Bessel function is a question still remaining open. The second problem comes from the impossibility to work with infinite sum in (Equations 90 , 95 , and 97) . Actually, all computer simulations should use finite number of terms that immediately leads to truncated Dirichlet series obeying the condition (Equation 97). This immediately, transforms the approximation of the relaxation function to approximation through Prony series, but the first problem formulated above, still remains unanswered.
Actually, the approximation by Prony series or infinite Dirichlet series, generally speaking, is an approach considering approximations of Non-Debye responses (relaxations) by superpositions of sub-processes (as Debye relaxations) with different relaxation times (an approach widely applied in the relaxation processes in glass transitions [123], for instance).
7. Final Comments
The ideas and results developed in this work focused on the new fractional operator conceived in [1, 10] (2015) naturally appearing when we stayed on the shoulders of two classical results: the Prony series (1795) [2] and the Boltzmann superposition principle (1864) [9].
It was demonstrated that in many cases there are viscoelastic materials which experimental behaviors exhibit strong departures from the power-law. In such cases it is natural to rise the questions about the adequate modeling of the dynamic processes in such media and to ask for new fractional operators. This is, actually, the same question raised by Bagley and Torvik [49, 124] who in case of power-law media suggested the constitutive relationship
Following Bagley and Torvik [49, 124] for a homogeneous viscoelastic materials the constitutive equation is
In (Equation 99) the derivatives have power-law kernels (Riemann-Liouville derivatives), which for N = 1 (one-term series of fractional derivatives) results in the simple expression
containing two fractional derivatives with different orders.
This article does not focus on development of different viscoelastic models based on the Caputo-Fabrizio fractional operator since this is out of its scope and draws new problems to be resolved. Despite this, if the kernels in the convolution integrals departure from the power-law we may consider a similar (formally) constitutive equation, namely
where, following the main idea of the Prony decompositions of the relaxation and the compliance, the series of fractional order of both sides of the equation have equal numbers of terms. Following [20, 107] the retardation times λi are satisfying the conditions
Therefore, from αi = 1/(1+τi/t0) and βi = 1/(1+λi/t0) it follows from Equation (102) that the fractional orders should satisfy the inequalities
The example of Renardy [125] (see also [26] and [14]), related to polymer rheology, reveals that a discrete relaxation spectrum with accumulation point at zero behaves as a power-law for short times, namely
Thus, in the asymptotic case (for t → 0 ) we may expect that the model (Equation 101) would converge to the model (Equation 100).
If we suggest only, for example, that N = 1, which actually is a departure of the main idea of Prony series, we have
This equation contains two fractional derivatives of different orders as in Equation (100). For instance, fractional viscoelastic models containing fractional derivatives and operators of two different order are thoroughly analyzed in Rossikhin and Shitikova [126] and Rossikhin and Shitikova [127] in the light of power-law materials. All these formal similarities and dissimilarities provoke new ideas that should be resolved.
To recapitulate, actually, we demonstrated that the Caputo-Fabrizio fractional operator naturally appears in the constitutive viscoelastic equations based on hereditary integrals when the material responses do not match the power-law. Moreover, it plays the same role as the fractional operators (derivatives) based on power-law memory kernels when they are inapplicable. At this moment the task of this study is accomplished and a lot of future works would be developed starting from results obtained here.
Statements
Author contributions
The author confirms being the sole contributor of this work and has approved it for publication.
Conflict of interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1.
CaputoMFabrizioM. A new definition of fractional derivative without singular kernel. Progr Fract Differ Appl. (2015) 1:73–85. 10.12785/pfda/010201
2.
deProny GR. Essai experimentale at analitique. J Ecole Polytech. (1795) 1:24–76.
3.
HristovJ. Derivatives with non-singular kernels: from the Caputo-Fabrizio definition and beyond: Appraising analysis with emphasis on diffusion models. In: Frontiers in Fractional Calculus, 2017 , BhalekarS editor. Sharjah: Bentham Science Publishers. (2017). pp. 269–342
4.
HristovJ. Transient heat diffusion with a non-singular fading memory: from the cattaneo constitutive equation with Jeffrey's kernel to the caputo-fabrizio time-fractional derivative. Thermal Sci. (2016) 20:765–70. 10.2298/TSCI160112019H
5.
GuistiA. A comment on some new definitions of fractional derivative. Nonlinear Dyn. (2018) 93:1757–63. 10.1007/s11071-018-4289-8
6.
OrtigueiraMDMachadoJAT. A critical analysis of the Caputo–Fabrizio operator. Commun Nonlinear Sci Numer Simul. (2018) 59:608–11. 10.1016/j.cnsns.2017.12.001
7.
TarasovVE. No nonlocality. No fractional derivative. Commun. Nonlinear Sci. Numer. Simul. (2018) 62:157–63. 10.1016/j.cnsns.2018.02.019
8.
PodlubnyI. Fractional Differential Equations. New York, NY: Academic Press (1999).
9.
BoltzmannL. Zur Theorie der Elastischen Nachwirkung. Sitzungsber Akad Wiss Wien Mathem Naturwiss. (1874) 70:275–300.
10.
CaputoMFabrizioM. Applications of new time and spatial fractional derivatives with exponential kernels. Progr Fract Differ Appl. (2016) 2:1–11. 10.18576/pfda/020101
11.
LosadaJNietoJJ. Properties of a new fractional derivative without singular kernel. Progr Fract Differ Appl. (2015) 1:87–92. 10.12785/pfda/010202
12.
HristovJ. The non-linear Dodson diffusion equation: Approximate solutions and beyond with formalistic fractionalization. Math Natural Sci. (2017) 1:1–7. 10.22436/mns.01.01.01
13.
HristovJ. Derivation of fractional Dodson's equation and beyond: Transient mass diffusion with a non-singular memory and exponentially fading-out diffusivity. Progr Fract Differ Appl. (2017) 3: 255–70. 10.18576/pfda/030402
14.
HristovJ. Response functions in linear viscoelastic constitutive equations and related fractional operators. Math Model Natur Phenomena. 10.1051/mmnp/2018067
15.
StormML. Heat conduction in simple metals. J Appl Phys. (1951) 22:940–51. 10.1063/1.1700076
16.
GurtinME. On the thermodynamics of materials with memory. Arch Rational Mech Anal. (1968) 28:40–50. 10.1007/BF00281562
17.
MillerRK. An integrodifferential equation for rigid heat conductors with memory. J Math Anal Appl. (1978) 66:313–32. 10.1016/0022-247X(78)90234-2
18.
NunciatoJW. On heat conduction in materials with memory. Q Appl Math. (1971) 29:187–273.
19.
PipkinAC. Lectures on Viscoelasticity Theory, 2nd Edn., Berlin: Sprimger-Verlag (1972).
20.
TschoeglNW. The Phenomenological Theory of Linear Viscoelastic Behaviour: An Introduction. New York, NY: Springer–Verlag (1989).
21.
BrinsonHBrinsonLPolymer Engineering Science and Viscoelasticity. New York, NY: Springer (2008).
22.
FerryJD. Viscoelastic Properties of Polymers. New York, NY: Wiley (1980).
23.
MarkovitzH. Boltzmann and the beginning of rheology. Trans Soc Rheol. (1977) 21:384–98. 10.1122/1.549444
24.
CiambellaJPaoloneAVidoliS. A comparison of nonlinear integral-based viscoelastic models through compression tests on filled rubber. Mech Mater. (2010) 42:932–44. 10.1016/j.mechmat.2010.07.007
25.
CiambellaJPaoloneAVidoliS. Identification of the viscoelastic properties of soft materials at low frequencies: Performance, ill-conditioning and extrapolation capabilities of fractional and exponential models. J Mech Behav Biomed Mater. (2014) 37:286–98. 10.1016/j.jmbbm.2014.05.031
26.
HanygaA. Wave propagation in media with singular memory. Math Comp Modell. (2001) 34:1399–421. 10.1016/S0895-7177(01)00137-6
27.
ColemanBNollW. Foundations of linear Viscoelasticity. Rev Modern Phys. (1961) 33:239–49. 10.1103/RevModPhys.33.239
28.
GarbarskiJ. The application of an exponentially-type function for the modeling of viscoelasticity of solid polymers. Polymer Eng Sci. (1992) 32:107–14. 10.1002/pen.760320206
29.
WinterHH. Analysis of dynamic mechanical data L inversion into relaxation time spectrum and consistency check. J Nonnewtonian Fluid Mech. (1997) 68:225–39. 10.1016/S0377-0257(96)01512-1
30.
EmriITschoeglNW. Determination of mechanical spectra from experimental data. Int J Solids Struct. (1995) 32:817–26. 10.1016/0020-7683(94)00162-P
31.
DayWA. The Thermodynamic Theory of Simple Materials With Fading Memory. Berlin: Springer (1972).
32.
FabrzioMGiorgiCNasoMG. Viscoelastic Solids of Exponential type. 1. Minimal representations and Contralability. Meccanica (2004) 39:531–46. 10.1007/s11012-003-6437-5
33.
FabrzioMGiorgiCNasoMG. Viscoelastic Solids of Exponential type. 2. Free energies . Stability and attractors. Meccanica (2004) 39:547–61. 10.1007/s11012-003-9400-6
34.
FabrizioMMorroA. Thermodynamic restrictions on relaxation functions in linear viscoelasticity. Mech Res Commun. (1985) 12:101–5. 10.1016/0093-6413(85)90077-1
35.
CrookRALettonA. The viscoelastic behaviour of notched glassy polymers. Eng Fracture Mech. (1993) 44:167–73.
36.
JaishankarAMcKinleyGH. A fractional K-BKZ constitutive formulation for describing the nonlinear theology of multiscale complex fluids. J Rheol. (2014) 58:1751–88. 10.1122/1.4892114
37.
JozwiakBOrczykowskaMDziubinskiM. Fractionalgeneralization of maxell and Kelvin-Voight models for biopolymer characterization. PLoS ONE (2015) 10:e0143090. 10.1371./journal.pone.0143090
38.
MeralFCRostonTJMaginR. Fractional calculus in viscoelasticity: an experimental studyCommun Nonlinear Sci Numer Simulat. (2010) 15:939–45. 10.1016/j.cnsns.2009.05.004
39.
MetzlerRSchickWKilianH-GNonnenmacherTF. Relaxation in filled polymers: a fractional calculus approach. J Chem Phys. (1995) 103:7180–6. 10.1063/1.470346
40.
SongKWKukHYChamgGS. Rheology of concentrated xanthan gum solutions: oscillatory shear flow behaviour. Korea Austr Rheol J. (2006) 18:18–24.
41.
WyattNBLiberatoreMW. Rheology and viscosity scaling of polyelectrolyte xanthan gum. J Appl Polymer Sci. (2010) 114:4076–84. 10.1002/app.31093
42.
SchiesselHMetzlerRBlumenANonnenmacherTF. Generalized viscoelastic models: their fractional equations with solutions. J Phys A Math Gen. (1995) 28:6567–84. 10.1088/0305-4470/28/23/012
43.
PapagiannnopoulosASotiropoulosKPispasS. Particle tracking microrheology of the power-law viscoelasticity of xantan solution. Food Hydrocolloids (2016) 61:201–10. 10.1016/j.foodhyd.2016.05.020
44.
CaputoM. Linear models of dissipation whose Q is almost frequency independent-II. Geophys J R Astr Soc. (1967) 13:529–39. 10.1111/j.1365-246X.1967.tb02303.x
45.
MainardiF. Fractional Calculus and Waves in Linear Viscoelasticity. London-Singapore: Imperial College Press and World Scienific (2010).
46.
GlockleWGNonnemacherTF. Fractional relaxation and the time-temperature superposition principle. Rheol Acta (1994) 33:337–43. 10.1007/BF00366960
47.
PadovanJGuoY. General response of viscoelastic systems modelled by fractional operators. J Franklin Inst. (1988) 325:247–75. 10.1016/0016-0032(88)90086-5
48.
FindleyWNLaiJSOnaranKCreep and Relaxation of Nonlinear Viscoelastic Materials. New York, NY: North-Holland Pub (1974).
49.
BagleyRLTorvikPJ. On the fractional Calculus model of viscoelastic behavior. J Rheol. (1986) 30:133–55. 10.1122/1.549887
50.
LakesRSVanderbyR. Interrelation of creep and relaxation: a modelling approach for ligaments. J Biomech Eng. (1999) 121:612–5. 10.1115/1.2800861
51.
Scott-BlairGWReinerM. The rheological law underlying the Nutting equation. Appl Sci Res. (1951) A2:225–34. 10.1007/BF00411984
52.
Scott-BlairGW. A model to describe the flow curves of concentrated suspensions of spherical particles. Rheol Acta(1965) 6:201–3. 10.1007/BF01976435
53.
Scott-BlairGW. Mathematics and rheology. Rheol Acta (1972) 11:237–40. 10.1007/BF01993026
54.
NuttingPG. A new general law of deformation. J Franklin Inst. (1921) 191:679–85. 10.1016/S0016-0032(21)90171-6
55.
NuttingPG. Deformation in relation to time, pressure and temperature. J Franklin Inst. (1946) 242:449–58. 10.1016/0016-0032(46)90636-9
56.
CarilloSGiorgiC. Non-classical memory kernels in linear viscoelasticity. In: Viscoelastic and Viscoplastic Materials, Mohamed FathyEl-Amin editor. Split: IntechOpen (2016). p. 295–331. 10.5772/64251
57.
BradshawRDBrinsonLC. A sing control method for fitting and interconverting material functions for linear viscoelastic solids. Mech Time Depend Mater. (1997) 1:85–108. 10.1023/A:1009772018066
58.
MandelbrotS. Dirichlet Series. Principle and Methods.Dordrecht: D. Reidel Publ. Co. (1972).
59.
CzyzJASzyszkowskiW. An effective method for non-linear viscoelastic structural analysis. Comput Struct. (1990) 37:637–46. 10.1016/0045-7949(90)90096-K
60.
CristensenRM. Theory of Viscoelasticity. An introduction.New York, NY: Academic Press (1982).
61.
BilstonLELiuZPhan-ThienN. Linear viscoelastic properties of bovine brain tissue in shear. Bioreology (1997) 34:377–85.
62.
PhillipsAPankajPMayFTaylorKHowieCUsmaniA. Constitutive models for impacted morsellised cortico-cancelous bone. Biomateials (2006) 27:2162–70. 10.1016/j.biomaterials.2005.10.034
63.
RavikumarNNobleCCramphornETaylorZA. A short constitutive model for ballistic gelatine at surgical strain rates. J Mech Behav Biomed Mater. (2015) 47:87–94. 10.1016/j.jmbbm.2015.03.011
64.
LondonoJGBerger-VergioatLWaismanH. A Prony series type viscoelastic solid coupled with continuous damage law for polar ice modeling. Mech Mater. (2016) 98:81–97. 10.1016/j.mechmat.2016.04.002
65.
BabaeiBDavarianAPryseKMElsonELGeninGM. Efficient and optimized identification of generalized viscoelastic relaxation spectra. J Mech Behav Biomed Mater. (2016) 55:32–41. 10.1016/j.jmbbm.2015.10.008
66.
MauroJCMauroYZ. On the Prony representation of stretched exponential relaxation. Phys A (2018) 506:75–87. 10.1016/j.physa.2018.04.047
67.
McGrumNGBuckleyCPBucknallCB. Principles of Polymer Engineering. 2nd Edn. Oxford: Oxford Science Publications (2003).
68.
JalohaDConstantinescuANeviereR. Revisiting the identification of the generalized Maxwell models from experimental results. Int J Solids Struct. (2015) 67–68: 169–81. 10.1016/j.ijsolstr.2015.04.018
69.
SunYChenJHuangB. Characterization of asphalt concrete linear viscoelastic behaviour utilizing the Havriliak-Negami complex modulus model. Constr Build Mater. (2015) 99:226–34. 10.1016/j.conbuildmat.2015.09.016
70.
GrossB., Mathematical Structure of the Theories of Viscoelasticity. Paris: Hermann (1953)
71.
GrossB. Distribution functions in linear viscoelastic theory. J Appl Phys. (1987) 62:2763–70. 10.1063/1.339404
72.
BaumgaertelMWinterH. Interrelation between continuous and discrete relaxation time spectra. J Non Newtonian Fluid Mech. (1992) 44:15–36. 10.1016/0377-0257(92)80043-W
73.
MeadDW. Numerical interconversion of linear viscoelastic functions. J Rheol. (1994) 38: 1769–95. 10.1122/1.550526
74.
FernandezPLamelaRey RMFenandezCanteli A. Viscoelastic characterisation of the temporomandibular joint disc in bovines. Strain (2011) 47:188–93. 10.1111/j.1475-1305.2008.00502.x
75.
Prieto-MunozPYinHTestaR. Mechaniscs of an adhesive anchor system subjected to a pullout load, II: viscoelastic analysis. J Struc Eng. (2013) 140:04013053. 10.1061/(ASCE)ST.1943-541X.0000822
76.
ShouZChenFYinH. Self-heating of a polymeric particulate composite under mechanical excitations. Mech Mater. (2018) 117:116–25. 10.1016/j.mechmat.2017.11.003
77.
LaksariKShafieianMDarvishK. Constitutive model for brain tissue under finite compression. J Biomech. (2012) 45:642–6. 10.1016/j.jbiomech.2011.12.023
78.
ChoiSTLeeSREarmmeYY. Flat indentation of viscoelastic polymer film on a rigid substrate. Acta Mater. (2008) 56:5377–87. 10.1016/j.actamat.2008.07.030
79.
SchaperyRA. Viscoelastic behavior and analysis of composite materials. In: Mechanics of Composite Materials, Vol. 2, SendeckyiGP editor. New York, NY: Academic Press (1974). p. 85–169.
80.
FajmanPSeinojaJ. A simplified approach to time-dependent subsoil-structure interaction. Comput. Struct. (2007) 85:1514–23. 10.1016/j.compstruc.2007.01.024
81.
TayebAArfaouiMZineAHamdiABenabdallahJIchchouM. On the non-linear viscoelastic behaviour of rubber-like materials: constitutive description and identificationInt J Mech Sci. (2017) 130:437–47. 10.1016/j.ijmecsci.2017.06.032
82.
VandenbergheEChoucharinaSLucaSDeKetelaere BDeBaerdemaeker JClaesJ. Spatio-temporal gradients of dry matter content and fundamental material parameters of Gouda cheese. J Food Eng. (2014) 142:31–48. 10.1016/j.jfoodeng.2014.05.019
83.
HillSA. The analytic representation of viscoelastic material properties using optimization techniques. NASA TM-108394, February, 199, (1993).
84.
MitraSMitraA. A genetic algorithms based techniques for computing the nonlinear least squares estimates of parameters of sum of exponential model. Expert Syst Appl. (2012) 39: 6370–9. 10.1016/j.eswa.2011.12.033
85.
HuangAWLuCHWuSCVinciRPBrownWLLinMT. Viscoelastic mechanical properties measurement of thin AL and Al-Mg films using bulge testing. Thin Solid Films (2016) 618:2–7. 10.1016/j.tsf.2016.03.064
86.
EnderleinJErdmannR. Fast fitting of multi-exponential decay curves. Optics Commun. (1997) 134:371–8. 10.1016/S0030-4018(96)00384-7
87.
DeHoffPHBarrettAALeeRBAnusaviceKJ. Thermal compatibility of dental ceramic system using cylindrical an spherical geometries. Dental Mater. (2008) 24:744–52. 10.1016/j.dental.2007.08.008
88.
DrozdovA. A constitutive models for nonlinear viscoelastic media. Int J Solids Struct. (1997) 34:2685–707. 10.1016/S0020-7683(96)00178-3
89.
AndrewsSHJRattnerJBShriveNGRonskyJL. Swelling significantly affects the material properties of the menisci in compression. J Biomech. (2015) 48:1485–9. 10.1016/j.jbiomech.2015.02.001
90.
KosaABerezvaiS. Visco-hyperelastic characterization of polymeric foam materials. Mater Today Proc. (2016) 3:1003–8. 10.1016/j.matpr.2016.03.037
91.
LinRC. On a nonlinear viscoelastic material law at finite strain for polymers. Mech Res Commun. (2001) 28:365–72. 10.1016/S0093-6413(01)00185-9
92.
LiuFLiCLiuSGeninGMHuangGLuTJXuF. Effects of viscoelasticity on skin pain sensation. Theor Appl Mech Lett. (2012) 5:222–6. 10.1016/j.taml.2015.11.002
93.
SlikerLJCiutiGRentschlerMMenxciasiA. Frictional resistance model for tissue-capsule endoscope sliding contact in the gastrointenstinal tract. Tribol Int. (2016) 2012:472–84. 10.1016/j.triboint.2016.06.003
94.
CuiTLinCWChienCHChaoYJVanZee JW. Service life estimation of liquid silicone rubber seals in polymer electrolyte membrane fuel cell environments. J Power Sources (2011) 196:1216–21. 10.1016/j.jpowsour.2010.08.075
95.
DrozdovAFinite Elasticity and Viscoelasticity. Singapore: World Scientific (1996).
96.
LaiJBakkerA. Analysis of the non-linear creep of high-density polyethylene. Polymer(1995) 36:93–9.
97.
PelayoFLamela-ReyMMinuz-CalventeMLopez-AenlleMAlvarez-VasquezA. Study of the time-temperature behaviour of PVB: Application of laminated glass elements. Thin Walled Struct.(2017) 119:324–31. 10.1016/j.tws.2017.06.030
98.
SchneiderJHilckenJAronenAKarvinenROlesenJFNielsenJ. Stress relaxation in tempered glass caused by hear soak testing. Eng Struct. (2016) 122:42–9. 10.1016/j.engstruct.2016.04.024
99.
ShengWAlcornRLewisA. A new method for radiation forces for floating platforms in waves. Ocean Eng. (2015) 105:43–53. 10.1016/j.oceaneng.2015.06.023
100.
TroyerKLEstepDJPuttlizChM. Viscoelasticity effects during loading play an integral role in soft tissue mechanics. Acta Biomater. (2012) 8:240–3. 10.1016/j.actbio.2011.07.035
101.
LavrentievMSome ill-Posed Problems of Mathematical Physics. Berlin: Springer (1967).
102.
HansenPC. Numerical tools for analysis and solution of Fredholm integral equations of the first kind. Inverse Prob. (1992) 8:849–72. 10.1088/0266-5611/8/6/005
103.
HonerkampJ. Ill-posed problems in rheology. Rheol Acta (1989) 28:363–71. 10.1007/BF01336803
104.
HonerkampJWeeseJ. Determination of the relaxation spectrum by a regularization method. Macromolecules (1989) 22:4372–7. 10.1021/ma00201a036
105.
AnderssenRSDaviesARdeHoog FR. On the interconversion integral equation for relaxation and creep. ANZIAM J. (2007) 48:346–63. 10.21914/anziamj.v48i0.64
106.
ParkSWSchaperyRA. Methods of Interconversion between linear viscoelastic material functions from single experiment. Int J Solids Sruct. (1999) 36:1653–75. 10.1016/S0020-7683(98)00055-9
107.
Luk-CyrJCrochonTliCLevesqueM. Interconversion of linearly viscoelastic material functions expressed as Prony series: a closure. Mech Time Depend Mater (2013) 17:53–82. 10.1007/s11043-012-9176-y
108.
PlazekDJEcheverriaI. Do not cry for me Charlie Brown, or with compliance comes comprehension. J Rheol. (2000) 44:931–41. 10.1122/1.551117
109.
HopkinsILHammingRW. On creep and relaxation. J Appl Phys. (1957) 28:906–9. 10.1063/1.1722885
110.
HristovJ. Steady-state heat conduction in a medium with spatial non-singular fading memory: derivation of caputo-fabrizio space-fractional derivative with Jeffrey's kernel and analytical solutions. Thermal Sci. (2017) 21:827–39. 10.2298/TSCI160229115H
111.
LoyRJdeHoog FRAnderssenRS. Interconversion of Prony series for relaxation and creep. J Rheol. (2015) 59:1261–70. 10.1122/1.4929398
112.
LoyRJAderssenRS. Interconversion relationships for completely monotone functions. SIAM J Math Anal. (2014) 46:2008–13. 10.1137/120891988
113.
CanestratiFStimilliABahiaHAVirgiliA. Pseudo-variables method to calculate HMA relaxation modulus through low-temperature induced stress and strain. Mater Design (2015) 76:141–9. 10.1016/j.matdes.2015.03.063
114.
ColombaroIGuistiAMainardiF. A class of linear viscoelastic models based on Bessel functtionsMeccanica (2017) 52:825–32. 10.1007/s11012-016-0456-5
115.
ColombaroIGuistiAMainardiF. On the propagation of transient waves in a viscoelastic Bessel mediumZ Angew Math Phys. (2017) 62:2017. 10.1007/s00033-017-0808-6
116.
ColombaroIGuistiA. Bessel models of linear viscoelasticity. Int J Theor Appl. Mech. (2018) 3:26–31.
117.
ColombaroIGuistiAMainardiF. A one parameter class of fractionalMaxwell-like models. Appl Math Comp Sci AIP Conf Proc. (2017) 1836:020003. 10.1063/1.4981943
118.
GuistiAMainardiF. On infinite series concerning zeros of Bessel functions of the first kind. Eur Phys J Plus(2016) 131:206. 10.1140/ep.jp/i2016-16206-4
119.
GuistiAMainardiF. A dynamic viscoelastic analogy for fluid-filled elastic tubes . Meccanica (2016) 51:2321–30. 10.1007/s11012-016-0376-4
120.
ColombaroIGuistiAVitaliS. Storage and dissipation of energy in Prabhakar viscoelasticity. Mathematics 2018 6:15. 10.3390/math6020015
121.
GuistiAColombaroI. Prabhakar-like fractional viscoelasticity. Comm Nonlinear Sci Numer Simulat. (2018) 56:138–43. 10.1016/j.cnsns.2017.08.002
122.
AbramowitzMStegunIAHandbook of Mathematical Functions. New York, NY: Dover (1965).
123.
HuntA. Non-Debye relaxation and the glass transition. J Non Crystal Solids (1993) 160:183–227.
124.
BagleyRLTorvikJ. Fractional calculus – a different approach to the analysis of viscoelastically damped structures. AIAA J. (1983) 21:741–8. 10.2514/3.8142
125.
RenardyM. Some remarks on the propagation and non-propagation of discontinuities in linear viscoelastic liquids. Rheol Acta (1982) 21:251–4. 10.1007/BF01515713
126.
RossikhinYuAShitikovaM. Analysis of viscoelastic rod dynamics via models involving fractional derivatives or operators of different orders. Shock Vibration Digest (2004) 36:3–26. 10.1177/0583102404039131
127.
RossikhinYuAShitikovaM. Comparative analysis of viscoelastic models involving fractional derivatives of different orders. Frac Calc Appl Anal. (2007) 10:111–21.
Summary
Keywords
linear viscoelasticity, response functions, interconversion, power-law response, non power-law responses, prony series, caputo-fabrizio operator
Citation
Hristov JY (2018) Linear Viscoelastic Responses: The Prony Decomposition Naturally Leads Into the Caputo-Fabrizio Fractional Operator. Front. Phys. 6:135. doi: 10.3389/fphy.2018.00135
Received
27 September 2018
Accepted
07 November 2018
Published
05 December 2018
Volume
6 - 2018
Edited by
Dumitru Baleanu, University of Craiova, Romania
Reviewed by
Yilun Shang, Northumbria University, United Kingdom; Aliyu Isa Aliyu, Federal University Dutse, Nigeria
Updates
Copyright
© 2018 Hristov.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jordan Yankov Hristov jordan.hristov@mail.bg
This article was submitted to Mathematical Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.