SciELO - Scientific Electronic Library Online

 
vol.68 issue2Time-dependent interactions in tunnelling dynamicsThermodynamics of viscous dark energy for the late future time universe author indexsubject indexsearch form
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • Have no similar articlesSimilars in SciELO

Share


Revista mexicana de física

Print version ISSN 0035-001X

Rev. mex. fis. vol.68 n.2 México Mar./Apr. 2022  Epub Mar 27, 2023

https://doi.org/10.31349/revmexfis.68.020703 

RESEARCH

GRAVITATION, MATHEMATICAL PHYSICS AND FIELD THEORY

Fractional viscoelastic models with novel variable and constant order fractional derivative operators

Krunal Kachhiaa 

J.F. Gomez-Aguilarb   

aDepartment of Mathematical Sciences, P. D. Patel Institute of Applied Sciences, Charotar University of Science and Technology, Changa, Anand 388421, Gujarat, India.

bCONACyT-Tecnológico Nacional de México/CENIDET. Interior Internado Palmira S/N, Col. Palmira, 62490, Cuernavaca, Morelos, México.


Abstract

This paper deals with the application of a novel variable-order and constant-order fractional derivative without singular kernel of AtanganaKoca type to describe the fractional viscoelastic models, namely, fractional Maxwell model, fractional Kelvin-Voigt model, fractional Zener model and fractional Poynting-Thomson model. For each fractional viscoelastic model, the stress relaxation modulus and creep compliance are derived analytically under the variable-order and constant-order fractional derivative without singular kernel. Our results show that the relaxation modulus and creep compliance exhibit viscoelastic behaviors producing temporal fractality at different scales. For each viscoelastic model, the stress relaxation modulus and creep compliance are derived analytically under novel variable-order and constant-order fractional derivative with no singular kernel.

Keywords: Fractional viscoelastic models; variable-order derivatives; relaxation modulus; fractional derivative operators

1. Introduction

Materials are more or less dissipative and show phenomena such as creep, relaxation and damping. Viscoelastic constitutive models are often useful to model such behavior. Viscoelasticity theory describes processes in which the state of a mechanical system depends on a complete record of actions that were performed on it. Many polymers initially undergo fast relaxation followed by very slow relaxation. This behavior can be described effectively by using fractional order operators in the constitutive laws of viscoelasticity. Fractional viscoelastic model used to described the laws of deformation for modeling the viscoelastic behavior of real materials [7]. Fractional differential models in linear viscoelasticity allow the creep and relaxation processes to be described adequately by means of simple relationships between stresses and strains with a relatively small number of adjustable parameters [15]. Fractional viscoelastic models lead to the non-Debay relaxation behavior typical of numerous experimental observations [6]. Fractional differential models in linear viscoelasticity take into account stochastic (micro-Brownian) motion of molecules at the microlevel in the phenomenological description of viscoelastic media [15,16]. Viscoelastic materials exhibit both elastic and viscous characteristics when undergoing deformation [22]. For elastic response, the stressstrain relationship is rate independent. For viscous behaviors, the stress depends on deformation rate representing as a stiffer response at a larger loading rate. The viscoelastic response shows a dependence on the deformation history [12]. The mechanics of these viscoelastic materials can be modeled as arrangements of lossless elastic springs and lossy viscous dampers (dashpots). For a distributed material the stress due to the elastic element is proportional to the strain and the stress due to the viscous element is proportional to the time rate of change of strain. The commonly used viscoelastic material models are Maxwell model and Kelvin-Voigt model. These provide a useful description of relaxation, creep and stress rate dependence for some viscoelastic materials [13].

The Fractional calculus is a powerful tool to characterize the viscoelastic behaviors due to long-time memory and global correlation information. Fractional order derivatives are used in fractional calculus models, which represent a relatively simple way to describe dynamics in complex systems. In fractional viscoelastic models, a fractional dashpot is proposed to describe the dependence of viscous response on deformation history, while an elastic spring is used to represent the time-independent elastic response. Through different combinations of elastic springs and fractional dashpots, various fractional viscoelastic models have been constructed such as fractional Maxwell model, fractional Kelvin-Voigt model and fractional Zener model.

The rheological representative of fractional Maxwell model which is consisted of a spring in series with a fractional dashpot as shown in Fig. 1a). The relationship between stress σ and strain є of fractional Maxwell model can be expressed as

dνσdtν+στν=Edνϵdtν, (1)

where E is the elastic modulus of the string and τ is the relaxation time of the fractional dashpot.

FIGURE 1 Fractional viscoelastic models: a) fractional Maxwell model, b) fractional Kelvin-Voigt model, c) fractional Zener model and d) fractional Poynting-Thomson model. 

The fractional Kelvin-Voigt model is composed by a spring in parallel with a fractional dashpot as shown in Fig. 1b). The constitutive relationship has the following form

dνϵdtν+ϵτν=σEτν. (2)

The fractional three-element models include the fractional Zener model and the fractional Poynting-Thomson model. The fractional Zener model is composed of an elastic spring and a fractional Maxwell model in parallel as shown in Fig. 1c). The stress is related with the strains as

dνσdtν+ντν=(E1+E2)dνϵdtν+E1ϵτν, (3)

where E 1 is the elastic modulus of the spring in the first branch, E 2 is the elastic modulus of the spring in the branch of the fractional maxwell element and τ is the relaxation time of the fractional dashpot. While the fractional Poynting-Thomson model is consisted of a spring and fractional Kelvin-Voigt element, which are arranged in series as shown in Fig. 1d). The corresponding constitutive relationship is represented

dνσdtν+ϵτ^ν=1E^1dνσdtν+1E^1+1E^2στν, (4)

where E^1 and E^2 are the corresponding elastic moduli of the two springs, t^au is the relaxation time of the fractional dashpot. Though the fractional Zener model and fractional Poynting-Thomson model are constructed from different combinations of elastic spring and fractional dashpot, it is found that the two models are actually equivalent. If the materials constants E 1, E 2 and τ^ν of fractional Zener model in Eq. (3) are replaced by (E^1E^2/[E^1+E^2]), (E^12/[E^1+E^2]) and (E^22τ^ν/[E^1+E^2]) respectively, the obtained constitutive relationship is identical to (4).

According to the general theory of linear viscoelasticity [25], given linear model can be expressed in terms of two equivalent forms, namely the creep and the relaxation representations. To each of these representations, two material functions G(t) and J(t) respectively called relaxation modulus and creep compliance can be associate. It is also important to stress that each of these two functions contains all the physical information about the viscoelastic model. In the Laplace domain, a given (linear) constitutive equation takes the following two equivalent forms [18,25]

σ-(s)=sG-(s)ϵ-(s),      ϵ-(s)=sJ-(s)σ-(s). (5)

These fractional viscoelastic models have been successfully applied to describe the viscoelastic behaviors of various material systems [10,11,20,24,25].

1.1. Fractional Calculus

In last decades, the concept of fractional calculus has been attracting the attention around the world due to its applications in complex systems. The model involving the fractional order derivatives gives a better prediction in comparison to model using integer order derivatives. Recently, the concepts of fractional differentiation and integration have both witnessed an historical development where the power law kernel was reviewed and replaced by exponential and Mittag-Leffler kernels [3,8,9,19]. These developments come from some arguments that were raised about the singularity of the power law around the origin zero. Generally speaking, a finite-time singularity will be observed when one input parameter or variable is time and the result of the output variable blow up toward infinity at a finite time. The notion of fractional differentiation with non singular kernel was then suggested to handle such natural occurrence and also to handle those physical problem where the singularities do not occur at the initial conditions. They have been applied in many fields of science with great success [2, 5, 21, 23]. While the concept of fractional derivative with constant order have been in fashion in the last decades, they have been found to be inefficient when dealing with anomalous diffusion therefore, the concept of variable orders was suggested. The variableorder fractional derivatives are very useful when investigating the memory properties which change with time and spatial location. Therefore, variable-order fractional derivative can be used to characterize variable memory effect of the system. Constant-order fractional operators are capable of describing some very relevant physical phenomena, it cannot capture important classes of physical phenomena where the order itself is a function of either dependent or independent variables. For example, the reaction kinetics of proteins has been found to exhibit relaxation mechanisms that are properly described by a temperature-dependent fractional order [17]. Thus, the underlying physics of the reaction kinetics (captured by the order of the relaxation mechanism) changes with temperature. Hence, it is reasonable to think that a differential equation with operators that update their order as a function of temperature will better describe the protein kinetics. This study suggests that there exist classes of physical problems that would be better described by variable-order fractional operators. The differential operators are excellent suggestion for anomalous diffusion, one will realizes that they cannot be used for analytical purpose [1,4].

Definition 1.1. The generalized Mittag-Leffler function Eδ,νλ(z) is defined as follows [27]

Eδ,νλ(z)=n=0(λ)nznΓ(δn+ν)n!,(δ,ν,λC,R(δ)>0,R(ν)>0,R(λ)>0),

where (λ) n is the Pochhammer symbol defined by

(λ)n=Γλ+nΓλ,     (λ)0=1,

(λ)n=k=1n(λ+k-1),n1.

Definition 1.2. Let fH 1(α,b) and 0 < δ ≤ 1, then Atangana-Koca derivative of fractional order in Caputo sense is given as [4]:

 AKC0Dtδ(f(t))=1g(δ)×0tf'(τ)Eδ,δδ(-g(δ)(t-τ)δ)  dτ, (6)

where the function g(δ) is well-defined such that

limδ01g(δ)0tf'(τ)Eδ,δδ(-g(δ)(t-τ)δ)  dτ=0tdf(τ)dt=f(t)-f(0).

The Laplace transform for the Atangana-Koca fractional derivative (6) is given by

LAKC0Dtδ(f(t))}(s)=1g(δ)(sf-(s)-f(0))×s-nδ-1(1-g(δ))δ. (7)

Definition 1.3. Let g(x) ∈ C 1[α,b] and f(x) a differential function in an open interval I. The Atangana-Koca fractional variable order derivative in Caputo sense is given by [1]

 AKV0Dtg(x)(f(t))=0tf'(τ)exp(-g(x)(t-τ))  dτ. (8)

The Laplace transform for the Atangana-Koca fractional variable order derivative is given by

L{AKV0Dtgδ(f(t))}(s)=sf-(s)-f(0)s+g(δ). (9)

The m-th derivative of the function Eδ,ν(m)(atδ) with δ > 0 and ν > 0 is given by

Eδ,ν(m)(atδ)=k=0(k+m)!k!aktδkΓ(δk+δm+ν).

The Laplace transform of the function tδm+ν-1Eδ,ν(m)(atδ) is given as

L{tδm+ν-1Eδ,νm(atδ)}(s)=m!sδ-ν(sδ-a)m+1,  δ>0,  νR,  R(s)>|a|1δ. (10)

Delgado et al [26] applied a new fractional derivative with variable and constant order applied to a mass-spring-damper system. Gómez-Aguilar [14] obtained analytical solutions to electrical circuits of non-integer order via fractional derivatives with and without singular kernels.

2. Maxwell model with Atangana-Koca fractional derivative

In this section, we will investigate the analytic solutions of the relaxation modulus and creep compliance for fractional Maxwell model with fractional derivative of variable and constant order with non singular kernel.

Maxwell model with AKC fractional derivative with constant order

Let us consider (1) involving the AKC fractional derivative with constant order (6) as follows:

 AKC0Dtν(σ(t))+στν=EAKC0Dtν(ϵ(t)). (11)

We apply Laplace transform to (11) and obtain the following relationship,

s-nνg(ν)(1-g(ν))νσ-(s)+σ-(s)τν=s-nνg(ν)(1-g(ν))νϵ-(s), (12)

which gives

σ-(s)=Es-nνs-nν+g(ν)(1-g(ν))ντ-νϵ-(s). (13)

The Laplace transform of relaxation modulus can then be derived as

G-(s)=σ-(s)sϵ-(s)=Eτνg(ν)(1-g(ν))ν×s-1snν+τνg(ν)(1-g(ν))ν. (14)

By applying the inverse Laplace transform to (14) and using (10), the relaxation modulus of fractional Maxwell model with Atangana-Koca fractional order derivative (6) is obtained as

GMAKC(t)=Eτνtnνg(ν)(1-g(ν))νEnν,nν+1×-τνtnνg(ν)(1-g(ν))ν. (15)

The Laplace transform of the creep compliance can be estimated from

J-(s)=1s2G-(s)=g(ν)(1-g(ν))νEτν×1s1-nν+1sτνg(ν)(1-g(ν))ν. (16)

The creep compliance then be obtained through inverse Laplace transform of (16) as following

JMAKC(t)=1E1+g(ν)(1-g(ν))ντνt-nνΓ(1-nν). (17)

Figures 2a)-2b) shows the mechanical response of fractional Maxwell model involving the Atangana-Koca fractional derivative for several ν values, arbitrarily chosen.

FIGURE 2 Mechanical properties of fractional Maxwell model involving the Atangana-Koca fractional derivative. In a) relaxation modulus and b) creep compliance for several ν values, arbitrarily chosen. 

Maxwell model with AKC fractional derivative with variable order

Let us consider (1) involving the AKC fractional derivative with variable order (8) as follows:

 AKV0Dtg(ν)(σ(t))+στg(ν)=EAKV0Dtg(ν)(ϵ(t)). (18)

We apply Laplace transform to (18) and obtain the following relationship,

sσ-(s)s+g(ν)+σ-(s)τg(ν)=Esϵ-(s)s+g(ν). (19)

which gives

σ-(s)=Esϵ(s)τg(ν)s(τg(ν)+1)+g(ν). (20)

By applying the inverse Laplace transform to (20) and using (10), the relaxation modulus of fractional Maxwell model with Atangana-Koca fractional variable order derivative (8) is obtained as

GMAKV(t)=Eτg(ν)τg(ν)+1e-g(ν)τg(ν)+1t. (21)

The Laplace transform of the creep compliance can be estimated from

J-(s)=1s2G-(s)=s(τg(ν)+1)+g(ν)s2Eτg(ν). (22)

The creep compliance then be obtained through inverse Laplace transform of (16) as following

JMAKV(t)=τg(ν)+1Eτg(ν)1+tg(ν)τg(ν)+1. (23)

Figures 3a)-3d) shows the mechanical response of fractional Maxwell model involving the Atangana-Koca fractional derivative with variable order for several g(ν) values, arbitrarily chosen.

FIGURE 3 Mechanical properties of fractional Maxwell model involving the Atangana-Koca fractional derivative with variable order. In a)-b) relaxation modulus and c)-d) creep compliance for several g(ν) values, arbitrarily chosen. 

3. Kelvin-Voigt model with Atangana-Koca fractional derivative

In this section, we will investigate the analytic solutions of the relaxation modulus and creep compliance for fractional Kelvin-Voigt model with fractional derivative of variable and constant order with non singular kernel.

Kelvin-Voigt model with AKC fractional derivative with constant order

Let us consider (2) involving the AKC fractional derivative with constant order (6) as follows:

 AKC0Dtν(ϵ(t))+ϵτν=σEτα. (24)

We apply Laplace transform to (24) and obtain the following relationship,

s-nνg(ν)(1-g(ν))νϵ-(s)+ϵ-(s)τν=ϵ-(s)Eτν, (25)

which gives

σ-(s)=s-nντν+g(ν)(1-g(ν))νg(ν)(1-g(ν))νEϵ-(s). (26)

The Laplace transform of relaxation modulus can then be derived as

G-(s)=σ-(s)sϵ-(s)=E1s+τνg(ν)(1-g(ν))ν1snν+1. (27)

By applying the inverse Laplace transform to (27) and using (10), the relaxation modulus of fractional Kelvin-Voigt model with Atangana-Koca fractional order derivative (6) is obtained as

GKVAKC(t)=E1+τνg(ν)(1-g(ν))νtnνΓ(nν+1). (28)

The Laplace transform of the creep compliance can be estimated from

J-(s)=1s2G-(s)=1Esnν-1snν+τνg(ν)(1-g(ν))ν. (29)

The creep compliance then be obtained through inverse Laplace transform of (16) as following

JKVAKC(t)=1EEnν,1-τνtnνg(ν)(1-g(ν))ν. (30)

Figures 4a)-4b) shows the mechanical response of fractional Kelvin-Voigt model involving the Atangana-Koca fractional derivative for several ν values, arbitrarily chosen.

FIGURE 4 Mechanical properties of fractional Kelvin-Voigt model involving the Atangana-Koca fractional derivative. In a) relaxation modulus and b) creep compliance for several ν values, arbitrarily chosen. 

Kelvin-Voigt model with AKC fractional derivative with variable order

Let us consider (2) involving the AKC fractional derivative with constant order (8) as follows:

 AKV0Dtg(ν)(ϵ(t))+ϵτg(ν)=σEτg(ν). (31)

We apply Laplace transform to (31) and obtain the following relationship,

ss+g(ν)+1τg(ν)ϵ-(s)=σ-(s)Eτg(ν), (32)

which gives

σ-(s)=s(τg(ν)+1)+g(ν)s+g(ν)Eϵ-(s). (33)

The Laplace transform of relaxation modulus can then be derived as

G-(s)=σ-(s)sϵ-(s)=E1s+τg(ν)s+g(ν). (34)

By applying the inverse Laplace transform to (34) and using (10), the relaxation modulus of fractional Kelvin-Voigt model with Atangana-Koca fractional variable order derivative (8) is obtained as

GKVAKV(t)=E(1+τg(ν)e-g(ν)t). (35)

The Laplace transform of the creep compliance can be estimated from

J-(s)=1E1s-τg(ν)s+g(ν)τg(ν)+1. (36)

The creep compliance then be obtained through inverse Laplace transform of (36) as following

JKVAKV(t)=1E1-τg(ν)1+τg(ν)e-g(ν)tτg(ν)+1. (37)

Figures 5a)-5d) shows the mechanical response of fractional Maxwell model involving the Atangana-Koca fractional derivative with variable order for several g(ν) values, arbitrarily chosen.

FIGURE 5 Mechanical properties of fractional Kelvin-Voigt model involving the Atangana-Koca fractional derivative with variable order. In a)-b) relaxation modulus and c)-d) creep compliance for several g(ν) values, arbitrarily chosen. 

4. Three-element model with Atangana-Koca fractional derivative

In this section, we will investigate the analytic solutions of the relaxation modulus and creep compliance for fractional Zener model and fractional Poynting-Thomson model with fractional derivative of variable and constant order with non singular kernel.

Zener model with AKC fractional derivative with constant order

Let us consider (3) involving the AKC fractional derivative with constant order (6) as follows:

 AKC0Dtν(σ(t))+στν=(E1+E2)AKC0Dtν(ϵ(t))+E1ϵτα. (38)

We apply Laplace transform to (38) and obtain the following relationship,

s-nνg(ν)(1-g(ν))νσ-(s)+σ-(s)τν=(E1+E2)s-nνg(ν)(1-g(ν))νϵ-(s)+E1ϵ-(s)τν, (39)

which gives

σ-(s)=(E1+E2)s-nντν+E1g(ν)(1-g(ν))νs-nντν+g(ν)(1-g(ν))νϵ-(s). (40)

The Laplace transform of relaxation modulus can then be derived as

G-(s)=σ-(s)sϵ-(s)=(E1+E2)τνg(ν)(1-g(ν))ν1s(snν+τνg(ν)(1-g(ν))ν)+E1snν-1snν+τνg(ν)(1-g(ν))ν. (41)

By applying the inverse Laplace transform to (41) and using (10), the relaxation modulus of fractional Zener model with Atangana-Koca fractional order derivative (6) is obtained as

GZAKC(t)=(E1+E2)τνtnνg(ν)(1-g(ν))νEnν,nν+1-τνtnνg(ν)(1-g(ν))ν+E1Enν,nν+1-τνtnνg(ν)(1-g(ν))ν. (42)

The Laplace transform of the creep compliance can be estimated from

J-(s)=1s2G-(s)=τνE1g(ν)(1-g(ν))ν1ssnν+E1+E2E1τνg(ν)(1-g(ν))ν+1E1snν-1ssnν+E1+E2E1τνg(ν)(1-g(ν))ν. (43)

The creep compliance then be obtained through inverse Laplace transform of (43) as following

JZAKC(t)=τνtnνE1g(ν)(1-g(ν))νEnν,nν+1-E1+E2E1τνtnνg(ν)(1-g(ν))ν+1EEnν,1-E1+E2E1τνtnνg(ν)(1-g(ν))ν. (44)

Figures 6a)-6b) shows the mechanical response of fractional Zener model involving the Atangana-Koca fractional derivative for several ν values, arbitrarily chosen.

FIGURE 6 Mechanical properties of fractional Zener model involving the Atangana-Koca fractional derivative. In a) relaxation modulus and b) creep compliance for several ν values, arbitrarily chosen. 

Zener model with AKC fractional derivative with variable order

Let us consider (3) involving the AKC fractional derivative with constant order (8) as follows:

 AKV0Dtg(ν)(σ(t))+στg(ν)=(E1+E2)AKV0Dtg(ν)(ϵ(t))+E1ϵτg(ν). (45)

We apply Laplace transform to (44) and obtain the following relationship,

ss+g(ν)+1τg(ν)σ-(s)=(E1+E2)ss+g(ν)+E1τg(ν)ϵ-(s), (46)

which gives

σ-(s)=((E1+E2)τg(ν)+E1)s+E1g(ν)s(τg(ν)+1)+g(ν)ϵ-(s). (47)

The Laplace transform of relaxation modulus can then be derived as

G-(s)=σ-(s)sϵ-(s)=E1s+E2τg(ν)τg(ν)+11s+g(ν)τg(ν)+1. (48)

By applying the inverse Laplace transform to (47) and using (10), the relaxation modulus of fractional Zener model with Atangana-Koca fractional variable order derivative (8) is obtained as

GZAKV(t)=σ-(s)sϵ-(s)=E1+E2τg(ν)τg(ν)+1e-g(ν)τg(ν)+1t. (49)

The Laplace transform of the creep compliance can be estimated from

J-(s)=1s2G-(s)=s(τg(ν)+1)+g(ν)s(s((E1+E2)τg(ν)+E1)+E1g(ν)) (50)

The creep compliance then be obtained through inverse Laplace transform of (36) as following

JZAKV(t)=1E1-(E1+E2τg(ν))E1((E1+E2)τg(ν)+E1)e-E1g(ν)(E1+E2)τg(ν)+E1t. (51)

Remark 4.1 If the materials constants E 1 , E 2 and τ^ν of fractional Zener model in Eq. (3) are replaced by E^1E^2/(E^1+E^2) , E^12/(E^1+E^2) and E^22τ^ν/(E^1+E^2) respectively, the obtained constitutive relationship is identical to (4). The analytical solutions of the relaxation modulus and creep compliance for fractional Poynting-Thomson model can be obtain from fractional Zener model.

Figures 7a)-7d) shows the mechanical response of fractional Zener model involving the Atangana-Koca fractional derivative with variable order for several g(ν) values, arbitrarily chosen.

FIGURE 7 Mechanical properties of fractional Zener model involving the Atangana-Koca fractional derivative with variable order. In a)-b) relaxation modulus and c)-d) creep compliance for several g(ν) values, arbitrarily chosen. 

5. Results and discussions

In this section, we compare the performance of fractional viscoelastic models with fractional derivatives of constant and variable order of Atangana-Koca type. The mechanical response of Maxwell model with new fractional derivatives of constant and variable order are shown in Figs. 2a)-2b) and 3a)-3d) respectively. It can be seen from Figs. 2a)-2b) that creep compliances with Atangana-Koca constant order derivative decrease sharply as ν → 1. The Atangana-Koca fractional derivative of constant order ν results in slowest relaxation of viscoelastic materials as ν → 0, which can be used to characterize ultraslow relaxation of viscoelastic materials. Numerical solutions of Eq. (23) have been depicted in Figs. 3a)-3d), for different values of the fractional order g(ν). It is seen that the system exhibits behaviors such that the solution continuously depends on the time variables. The creep compliance of fractional Maxwell model with variable order derivative exponential increase with time instead of linear increase compared with fractional Maxwell-model with constant order fractional derivative. The numerical results indicate that the fractional order has an important influence on the relaxation modulus and the creep compliance and the general solution of the fractional equations depends on the parameter ν and when the order is a function rather than a constant of arbitrary order g(ν), respectively.

Figures 4a)-4b) and 5a)-5d) shows the mechanical response of fractional Kelvin-Voigt model with Atangana-Koca fractional derivative with constant and variable order. The creep compliance of fractional Kelvin-Voigt model exponential increase with time instead of linear increase compared with fractional Maxwell-model.

Figures 6a)-6b) plot the relaxation modulus and creep compliance for different ν of fractional Zener model with Atangana-Koca fractional derivative with constant, the relaxation modulus exponential decreases with increasing value of ν.

6. Conclusions

The analytical solutions of the stress relaxation modulus and creep compliance for fractional Maxwell model, fractional Kelvin-Voigt model, fractional Zener model and fractional Poynting-Thomson model are obtained. It is found that the fractional Zener model has similar performance as fractional Maxwell model. These solutions represent a new family of solutions for the viscoelastic models, which allows for the possibility of multiple solutions that are not observed in experiments (complex dynamic behaviors are observed). For different time periods we have different memory abilities. The representation of generalized viscoelastic models by variable-order and constant-order fractional derivative also allows a deeper insight into the physics behind fractional stress-strain relations. The variable-order fractional derivative can be used to characterize variable memory effect of the systems. It is shown that the constant-order model is not able to capture all the details of the variable-order solution, particularly in the areas of transition between dynamic regimes. We concluded that this new fractional variable-order operator is more suitable for modeling real-world complex problems than all existing fractional variable-order operators. The experimental verification of these models will be reported in a future paper.

Conflicts of Interest

The authors declare no conflict of interest.

Acknowledgments

José Francisco Gómez Aguilar acknowledges the support provided by CONACyT: Cátedras CONACyT para jóvenes investigadores 2014 and SNI-CONACyT.

References

1 B. T. Alkahtani, I. Koca and A. Atangana, A novel approach of variable order derivative: Theory and Methods, J. Nonlinear Sci. Appl., 9 (2016), 4867-4876. [ Links ]

2 B.S.T. Alkahtani, Chua’s circuit model with Atangana-Baleanu derivative with fractional order, Chaos Solitons and Fractals, 89 (2016), 547-551. [ Links ]

3 A. Atangana and D. Baleanu, New fractional derivatives with nonlocal and non-singular kernel: theory and application to heat transfer model, Thermal Science, 20(2) (2016), 763-769. [ Links ]

4 A. Atangana and I. Koca, New direction in fractional differentation, Math. Nat. Sci., 1 (2017), 18-25. [ Links ]

5 A. Atangana and J.F. Gomez Aguilar, Decolonisation of fractional calculus rules: Breaking commutativity and associativity to capture more natural phenomena, Eur. Phys. J. Plus, 133 (2018), 133-166. [ Links ]

6 R. L. Bagley and P. J. Torvik, A theoretical basis for the application of fractional calculus to viscoelasticity, J. Rheol., 27 (1983), 201-210. [ Links ]

7 G. S. Blair, The role of Psychophysics in Rheology, Journal of Colloid Science, 2(1) (1947), 21-32. [ Links ]

8 M Caputo and M Fabrizio, A new definition of fractional derivative without singular kernel, Prog. Fract. Differ. Appl., 1(2) (2015), 73-85. [ Links ]

9 M Caputo and M Fabrizio, On the notion of fractional derivative and applications to the hysteresis phenomena, Meccanica, 52(13) (2017), 3043-3052. [ Links ]

10 B. Carmichael, H. Babahosseini, H. Mahmoodi and S. N. Agah, The fractional viscoelastic response of human breast tissue cells, Phys. Biol., 12(4) (2015), 046001. [ Links ]

11 Y. Carrera, G. Avila-de la Rosa, E. J. Vernon-Carter and J. Alvarez-Ramirez, A fractional order Maxwell model for NonNewtonian fluids, Physica A: Stat. Mech. Appl., 482 (2017), 276-285. [ Links ]

12 J. D. Ferry, Viscoelastic properties of Polymers, 3rd edition, John Wiley and Sons, New York, 1980. [ Links ]

13 Y. C. Fung, Biomechanics: mechanical properties of living tissues., 2nd edition, Springer-Verlag, New York, 1993. [ Links ]

14 J. F. Gómez-Aguilar, Fundamental solutions to electrical circuits of non-integer order via fractional derivatives with and without singular kernels, Eur. Phys. J. Plus ., 133 (2018) 133. [ Links ]

15 W. G. Glockle and T. F. Nonnenmacher, Fractional integral operators and Fox functions in the theory of viscoelasticity, Macromolecules, 24 (1991), 6424-6434. [ Links ]

16 W. G. Glockle and T. F. Nonnenmacher, Fractional relaxation and the time-temperature superposition principle, Rheol. Acta, 33 (1994), 337-343. [ Links ]

17 W. G. Glockle, and T. F. Nonnenmacher, A fractional calculus approach to self-similar protein dynamics, Biophys. J., 68 (1995), 46-53. [ Links ]

18 M. E. Gurtin and E. Sternberg, On the linear theory of viscoelasticity, Arch. Rational Mech. Anal., 11 (1962), 291-356. [ Links ]

19 J Hristov, Derivatives with non-singular kernels from the Caputo-Fabrizio definition and beyond: appraising analysis with emphasis on diffusion models, Frontiers in fractional calculus, Bentham Science Publishers, Sharjah (2017), pp. 270-342. [ Links ]

20 L. Jianmin, X. Rui and C. Wen, Fractional viscoelastic models with non-singular kernels, Mechanics of Materials, 27 (2018), 55-64. [ Links ]

21 D.C. Labora, J.J. Nieto, R. Rodriguez-López. Is it possible to construct a fractional derivative such that the index law holds., Progr Fract Differ Appl, 4(1) (2018), 1-3. [ Links ]

22 R. Lakes, Viscoelastic Materials, Cambridge University Press, Cambridge, 2009. [ Links ]

23 K. Ervin Lenzi, A. Tateishi Angel and V. Haroldo Ribeiro, The role of fractional time-derivative operators on anomalous diffusion, Front. Phys., 5 (2017), 1-9. [ Links ]

24 D. Lei, Y. J. Liang and R. Xiao, A fractional model with parellel fractional Maxwell elements for amorphous thermoplastic, Physica A: Stat. Mech. Appl. , 490 (2018), 465-475. [ Links ]

25 F. Mainardi, Fractional Calculus and Waves in Linear Viscoelasticity, Imperial College Press, London, 2010. [ Links ]

26 V. F. Morales-Delgado, J. F. Gómez-Aguilar, M. A. Taneco-Hernández and R. F. Escobar-Jiménez, A novel fractional derivative with variable and constant order applied to a massspring-damper system, Eur. Phys. J. Plus ., 133 (2018), 78. [ Links ]

27 T. R. Prabhakar, A singular integral equation with a genearlized Mittag-Leffler function in the kernel, Yokohama Math. J., 19 (1971), 7-15. [ Links ]

Received: December 27, 2020; Accepted: June 02, 2021

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License