Determination of Nonlinear Creep Parameters for Hereditary Materials

This work proposes an effective algorithm for description of nonlinear deformation of hereditary materials based on Rabotnov’s method of isochronous creep curves. The notions have been introduced for experimental and model rheological parameters and similarity coefficients of isochronous curves. It has been shown how using them, one can find instantaneous strains at various stress levels for description of nonlinear deformation of hereditary materials at creep. Relevant equations have been determined from the nonlinear integral equation of Yu. N. Rabotnov for the application cases of Rabotnov’s fractional exponential kernel and Abel’s kernel for nonlinear deformation of hereditary materials at creep. The improved methods have been given for determination of creep parameters α, ε0, δ, β, and λ. By processing and using test results for material Nylon 6 and glass-reinforced plastic TC 8/3-250, the process has been shown for sequential implementation of the developed methods for description of linear and nonlinear deformation of these materials at creep. From the results of the experimental investigation performed by the authors of this paper, it has been determined that fine-grained, dense asphalt concrete at the temperature of 20± 2 ◦C and stresses up to 0.183 MPa at direct tension is deformed considerably in a nonlinear way. It has been shown in an illustrative way by construction of isochronous creep curves at various load durations and curves of experimental rheological parameter at various stresses. Nonlinear deformation of asphalt concrete at creep is adequately described by the proposed methods.


Introduction
Many natural (soils, rocks, wood, natural asphalt) and artificial (metals, their alloys, polymers, concretes, composites) materials, depending on temperature and load level, to one extent or another show their viscoelastic properties.Currently there are sufficiently developed viscoelastic theory and methods [1][2][3][4][5][6] which allow the determination and description of viscoelastic properties of the materials.One can distinguish linear and nonlinear viscoelasticity [4][5][6][7] in viscoelastic theory, as well as in elastic theory [8,9].One can say that currently the theory and methods of linear viscoelasticity have been developed sufficiently.In spite of the fact that as far back as 1913 V. Volterra [10] proposed to describe nonlinear viscoelasticity by binary integral equation, to date, nonlinear viscoelasticity theory and methods are on the stage of development and currently calculations for strength, stability, and longevity in many branches of engineering activity are made without consideration or with poor consideration of viscoelastic properties for the materials and elements of structures.
One simple but efficient means for consideration of deformation nonlinearity for hereditary materials was proposed in 1948 by Yu.N. Rabotnov [11].It was based on the similarity of isochronous creep curves of the materials.Meanwhile, the process of creep strain is described by the equation of Boltzmann-Volterra for linear viscoelastic materials, but strain in the left part of the equation has been replaced by the so-called "curve of instantaneous deformation", which is determined experimentally.Yu.N. Rabotnov, as the author of the method, considers that the curve of instantaneous strain is a kind of ideal imaginary curve, which is impossible to obtain in reality, as in real conditions the strain rate is always a finite value [12,13].It can be obtained from isochronous creep curves at finite values of time, considering them similar.
It is known that the selection of kernel of the integral equation and determination of its parameters is one of the most responsible actions in the description of mechanical behavior of viscoelastic materials.
As it has been said in our previous work [14], the most universal is creep kernel in the form of the fractional exponential function of Yu.N. Rabotnov [4][5][6]11].The fractional exponential function is well studied and a special table [15] has been developed for simplification of calculations.
This work proposes an effective algorithm for description of nonlinear deformation of hereditary materials based on Rabotnov's method of isochronous creep curves, and also the methods have been improved for determination of creep kernel parameters, developed in our previous work [14].
The term "hereditary" was introduced by V. Volterra and means delay of elastic strain in time [13].A.S. Lodge explains that the materials are considered "hereditary", described by creep kernel, depending on the difference of an argument: K(t − τ) [25].

Short Description of Method
An important matter in the analysis of mechanical behavior of any material is the clarification of its deformation character: linearity or nonlinearity, which is usually determined by construction of strain dependence on stress according to the results of experimental tests.A more efficient and illustrative method is known for determination of deformation nonlinearity degree for materials, proposed by Yu.N. Rabotnov [4][5][6]11].In accordance with this method, the so-called isochronous creep curves are constructed according to the test results of samples for the materials for creep at some constant stresses.
Figure 1 represents creep curves for a material at various constant stresses.Let us draw several, for example, four vertical lines, corresponding to time moments t 1 , t 2 , t 3 , and t 4 .Each vertical line crosses four creep curves corresponding to stresses σ 1 , σ 2 , σ 3 , and σ 4 .Drawing horizontal lines from specified cross points, one can find four strain values ε 1 (t 1 ), ε 2 (t 1 ), ε 3 (t 1 ), and ε 4 (t 1 ) on the vertical axis, corresponding to time t 1 .According to known four values of stress and four values of strain, one can construct an isochronous creep curve corresponding to time moment t 1 .Making similar actions, one can construct isochronous creep curves for the material at different load durations (Figure 2).Isochronous curves allow the visual evaluation of the character (linearity or nonlinearity) and nonlinearity degree of deformation for material.It is often in practice that isochronous curves are similar.The similarity property provides an opportunity to obtain all other isochronous creep curves of a material, if one of them is known.

Nonlinear Equation and Creep Kernel
Considering the similarity property of isochronous creep curves, Yu.N. Rabotnov proposed the following nonlinear integral equation to describe the process of nonlinear deformation of hereditary materials [4][5][6]11].

Nonlinear Equation and Creep Kernel
Considering the similarity property of isochronous creep curves, Yu.N. Rabotnov proposed the following nonlinear integral equation to describe the process of nonlinear deformation of hereditary materials [4][5][6]11].

Nonlinear Equation and Creep Kernel
Considering the similarity property of isochronous creep curves, Yu.N. Rabotnov proposed the following nonlinear integral equation to describe the process of nonlinear deformation of hereditary materials [4][5][6]11].
where ε(t) is strain at time moment t;  1) represents by itself the so-called "instantaneous deformation curve", determination of which will be given below.

Methods of Determination of Linear Creep Parameters for Hereditary Materials
Nearly always one can accept that the relationship between stress and strain in materials is a linear one at weak stresses.Therefore, one can use the creep curve of a material at weak stress and apply the linear viscoelasticity approach [14] for determination of creep parameters.

Parameters α, ε 0 , and δ
As it is known, creep curves of materials, depending on stress level and temperature, can have two or three characteristic strain sites [26][27][28]: site I with unstabilized creep, site II with stabilized creep, and site III of accelerating creep.
Then the method is proposed, according to which it will be sufficient to consider only site I of creep curve for determination of creep parameters of a material α, ε 0 , and δ.Having accepted n = 0, from Equation (3) we will find: One can see that the right part of the obtained equation contains a well-known Abel's function with unknown parameters α and δ.
Having divided both parts of the Equation (4) into instantaneous elasticity modulus E 0 for the case of linear deformed material, we will obtain: where ε 0 is conditionally instantaneous strain.Equation (5) contains three unknown parameters ε 0 , α, and δ.As it has been said above, the right part contains a known Abel's function with parameter of singularity α, which has the value within the interval (0, 1).Then parameter α will be considered as known, and unknown parameters ε 0 and δ will be determined with the use of the least squares method.According to the least squares method, the values of parameters ε 0 and δ should meet the following condition: where S(ε 0 , δ) is the sum of squares of deviations; ε ei are values of creep strain determined experimentally; m is number of creep strain.
Calculating under the formula the deviations of the calculated values for creep strain from those obtained experimentally, one can select optimum values for parameters α, ε 0 , and δ, providing the least value of ∆ε i .

Parameters β and λ
Having divided both parts of Equation ( 3) into instantaneous elasticity modulus E 0 , for the case of linear strain of material, we obtain: Let us rewrite Equation (10) in the following form, inserting notation where Using the least squares method, we write the extremum condition, similar to condition of (6): The following expressions are obtained for determination of parameters β and λ from two equations, composed on expressions Appl.Sci.2018, 8, 760 where The value of parameter β is determined from Equation ( 14).If the equation has an unambiguous solution, then it will be the target value for parameter β > 0. It is obvious that during determination of values for parameters β and λ from Equations ( 14) and ( 15), the previously calculated values are used for parameters α and ε 0 .
As it is known, one of the factors which limited the wide use of creep kernel with Rabotnov's fractional exponential function has been the fact that it was not an easy way to determine its parameters [22].A special table has been developed [15], and with its use, the fractional exponential function has been calculated under the argument x = β t (1−α) at x ≤ 4. At x > 4 it is recommended to use asymptotic formulas of B. D. Anin [29].Calculations have shown that the infinite series in the right part of Equation ( 10) are poorly converged at small values of argument x < 4.

Introduction of special notation
and consideration in series ( 12) and ( 16) of not time t i , but relative time (t i /t k ) have essentially increased the convergence of calculations and allowed the quick obtaining of more accurate values of the fractional exponential function practically at any values of argument.

Algorithm for Description of Nonlinear Deformation of Hereditary Materials
This section of the paper proposes a new method for description of nonlinear deformation of hereditary materials.This method is realized in several steps and performed in the following sequence: a.A new parameter is introduced called the experimental rheological parameter: where ε e (t) is the value of creep strain at time moment t, determined experimentally; ε e 0 is instantaneous strain at time moment t = 0, determined experimentally.Experimental rheological parameter k e (t) represents by itself a normalized time function in relation to experimental instantaneous strain.It has the value equal to 1 at t = 0 and more than 1 at time values of t > 0. It shows in how many times the experimental values of creep strain are more at different time values compared with instantaneous strain, which has been also obtained experimentally.
b.The values have been calculated for the experimental rheological parameter k e (t) at different time values t and stresses σ.The graphs for k e (t) at various stresses σ are constructed according to the calculation results.All curves k e (t) at various stresses coincide (for example, as in the work [30]) for physically linear material, i.e., we have only one mean curve k e (t) for all stresses (Figure 3a).There are separate curves for various stresses k e (t) (Figure 3b) for physically nonlinear material.Coincidence for some of them is not excluded.
according to the calculation results.All curves () e k t at various stresses coincide (for example, as in the work [30]) for physically linear material, i.e., we have only one mean curve As it has been said above, the construction of isochronous creep curves is an illustrative way for the evaluation of physical nonlinearity (linearity) of the material.Therefore, we have two ways of visual evaluation for nonlinearity (linearity) of deformation for materials: first-through experimental rheological parameter () e k t ; second-by construction of isochronous creep curves.c.Using Rabotnov's fractional exponential function (3) and Abel's function (4), we approximate the experimental values of creep strain for material, obtained at minimum stress σ1.Meanwhile, we find corresponding values of creep parameters α, ε0, δ, β, and λ according to the abovementioned methods.
d.In case of creep, i.e., at σ = const from integral Equation ( 1) we obtain: The right part of this integral equation is denoted through km(t): As it has been said above, the construction of isochronous creep curves is an illustrative way for the evaluation of physical nonlinearity (linearity) of the material.Therefore, we have two ways of visual evaluation for nonlinearity (linearity) of deformation for materials: first-through experimental rheological parameter k e (t); second-by construction of isochronous creep curves.
c. Using Rabotnov's fractional exponential function ( 3) and Abel's function ( 4), we approximate the experimental values of creep strain for material, obtained at minimum stress σ 1 .Meanwhile, we find corresponding values of creep parameters α, ε 0 , δ, β, and λ according to the abovementioned methods.
d.In case of creep, i.e., at σ = const from integral Equation ( 1) we obtain: The right part of this integral equation is denoted through k m (t): and it is named as model (theoretical or calculation) rheological parameter.
As it is seen from Equation (19), model rheological parameter k m (t) depends only on time t.Similar to experimental rheological parameter k e (t), the model rheological parameter has the value equal to 1 at t = 0 and more than 1 at t > 0. It shows in how many times the calculated values of creep strain are bigger compared with the instantaneous strain, determined by calculations.Whereas the experimental rheological parameter can be determined for all stresses, the model rheological parameter is determined only for minimum stress σ 1 .
e.According to Yu. N. Rabotnov's isochronous creep curves method, we find the calculated values of instantaneous strains at various stresses under the following formula: where ε e (t s ) is the experimental value of creep for material at fixed time t s at stress σ; k m (t s ) is the value of model rheological parameter at fixed time t s (similarity coefficient of isochronous creep curves according to Yu. N. Rabotnov).
f.We calculate theoretical values of creep strain at various stresses under the formula: where ε m 0 (σ) is the calculated value of instantaneous strain at stress σ, obtained under the formula (20).g.We compare calculated ε m (t) and experimental ε e (t) values of creep strain at various stresses σ.
If the comparison has shown unsatisfactory accuracy, we determine again the values of instantaneous strains at various stresses under the following formula: where ε m 0 (σ) is a mean value of instantaneous strain at stress σ, obtained at some fixed times t s ; N is total number of fixed times t s .Then the calculated values of creep strain at various stresses are made under the formula:

Material Nylon 6
The works [5,20] contain test results of material Nylon 6 at stresses 5, 10, and 15 MPa.The duration of the experiment was 100 h for all stresses.Values of creep strain for material Nylon 6 at specified stresses, obtained by processing of experimental results, are represented in Table 1.It has been also determined from the data of mentioned works that the instantaneous deformation curve of material Nylon 6 has been satisfactorily approximated by power function: 1.3334 (24) Approximation of creep strain values at minimum stress (σ = 5 MPa) using Rabotnov's fractional exponential kernel gave the following values of creep parameters: α = 0.85, ε 0 = 0.1682, β = 0.18, and λ = 1.6682.
Then the creep curve of material Nylon 6 at stress σ = 5 MPa will be described by equation From Equation ( 24) at σ = 5 MPa we find that ε m 0 = 0.1682%.Considering this fact, we determine the model rheological parameter: Appl.Sci.2018, 8, 760 9 of 17 Now, one can calculate the value of creep strain at stresses 10 and 15 MPa under the formula (21).For these stresses, the values of instantaneous strain, obtained under formula (24), are equal to 0.4238% and 0.7277%, respectively.
Experimental and calculated values of creep strain for material Nylon 6 at all three stresses are shown in Figure 4.As it is seen, convergence of calculated strains to experimental ones is good.
From Equation ( 24 Now, one can calculate the value of creep strain at stresses 10 and 15 MPa under the formula (21).For these stresses, the values of instantaneous strain, obtained under formula (24), are equal to 0.4238% and 0.7277%, respectively.
Experimental and calculated values of creep strain for material Nylon 6 at all three stresses are shown in Figure 4.As it is seen, convergence of calculated strains to experimental ones is good.

Glass-Reinforced Plastic TC 8/3-250
Samples of glass-reinforced plastic TC 8/3-250 have been tested for creep at the temperature of 23.5 ± 2 • C in the works [21,31].As the glass-reinforced plastic is an anisotropic material, the samples have been cut along textile warp (Θ = 90 • ) and at an angle of Θ = 45 • to them.

Θ = 90 •
Values of creep strain for glass-reinforced plastic (Θ = 90 • ) at four stresses are represented in Table 2. Test results have shown that the sample of glass-reinforced plastic, cut at an angle of Θ = 90 • to textile warp, is deformed linearly at a specific temperature and up to stress 349 MPa.Instantaneous strain has been satisfactorily approximated by equation: Creep strains at minimum stress, equal to 104.7 MPa, have been approximated by Rabotnov's fractional exponential kernel.Creep parameters have the following values: α = 0.8, ε 0 = 0.3501, β = 0.3, and λ = 0.0462.
The experimental and calculated values of creep strain for glass-reinforced plastic (Θ = 90 • ) at all four stresses are shown in Figure 5 for comparison.As it is seen, convergence of calculated creep strains to experimental ones is good.

Θ = 45 •
Values of creep strain for glass-reinforced plastic (Θ = 45 • ) at six stresses are represented in Table 3.According to test results, it has been determined that the samples of glass-reinforced plastic, cut at the angle of Θ = 45 • to textile warp, have been deformed nonlinearly at the tested temperature and applied stresses.
The model rheological parameter has a form of: Using the abovementioned values of instantaneous strain, one can calculate values of creep strain at all other stresses under formula (21).
Figure 6 shows experimental and calculated values of creep strain for glass-reinforced plastic (Θ = 45 • ) at all considered stresses.It is seen that convergence of experimental strains to the calculated ones is good.Processing of test results allowed determination of the following values of instantaneous strain for stresses 20.3 MPa, 40.6 MPa, 60.9 MPa, 81.2 MPa, 101.5 MPa, and 121.8 MPa: 0.1093%, 0.1982%, 0.5911%, 1.3872%, 2.4874%, and 3.9975%, respectively.
The model rheological parameter has a form of: Using abovementioned values of instantaneous strain, one can calculate values of creep strain at all other stresses under formula (21).
Figure 6 shows experimental and calculated values of creep strain for glass-reinforced plastic (Ѳ = 45°) at all considered stresses.It is seen that convergence of experimental strains to the calculated ones is good.
In this paper, bitumen of grade 100-130 has been used which meets the requirements of the Kazakhstan standard ST RK 1373-2013 [34].The bitumen grade on Superpave is PG 64-40 [35].
In this paper, bitumen of grade 100-130 has been used which meets the requirements of the Kazakhstan standard ST RK 1373-2013 [34].The bitumen grade on Superpave is PG 64-40 [35].Bitumen has been produced by Pavlodar processing plant from crude oil of Western Siberia (Russia) by the direct oxidation method.Bitumen content of grade 100-130 in the asphalt concrete is 4.8% by weight of dry mineral material.
Samples of the hot asphalt concrete in the form of a rectangular prism with the length of 150 mm, width of 50 mm, and height of 50 mm were manufactured in the following way.Firstly, the asphalt concrete samples were prepared in the form of a square slab by means of the Cooper compactor (UK, model CRT-RC2S) according to the standard EN 12697-33 [36].Then the samples were cut from the asphalt concrete slabs in the form of a prism.Deviations in sizes of the beams didn't exceed 2 mm.
Tests of the asphalt concrete samples on creep at the temperature of 22 ± 2 • C were carried out according to the direct tensile scheme in specially contracted equipment.The detailed information regarding characteristics of bitumen, asphalt concrete, and equipment can be found in the work [28].
The stresses were constant during the tests and equal to 0.041, 0.074, 0.111, 0.148, and 0.183 MPa.Durations for tests have been accepted as similar ones and equal to 570 s. 10 parallel asphalt concrete samples have been tested at each constant stress.
Mean values of creep strain for the tested 10 asphalt concrete samples are represented in Table 4. Creep curves of asphalt concrete at various stresses, constructed according to data of Table 4, are represented in Figure 7.As it is seen, the creep curves represent by themselves the lines which change monotonously with the increase of stress duration.The value of instantaneous strain (at t = 0) and slope of creep curves increase with the increase of stress.Creep curves of asphalt concrete at various stresses, constructed according to data of Table 4, are represented in Figure 7.As it is seen, the creep curves represent by themselves the lines which change monotonously with the increase of stress duration.The value of instantaneous strain (at t = 0) and slope of creep curves increase with the increase of stress.Then, using the values of stresses and strains, which correspond to various load durations, one can construct the isochronous creep curves (Figure 8).As it is seen, the asphalt concrete is deformed Then, using the values of stresses and strains, which correspond to various load durations, one can construct the isochronous creep curves (Figure 8).As it is seen, the asphalt concrete is deformed considerably in a nonlinear way at the given temperature, and the nonlinearity degree of relationship between stress and strain has been increased with the stress and strain duration increase.Nonlinearity of strain for asphalt concrete has been noted by other researchers [37,38].It should be also noted that the isochronous creep curves are similar ones.The property of similarity gives an opportunity to obtain all other isochronous creep curves, if any of them is known.
At present, it has been accepted in road science and practice that mechanical properties of asphalt concrete depend on the value and duration of the applied load and temperature, and this dependence can be described in a sufficiently accurate way with the use of linear viscoelasticity theory [39][40][41][42].The essential nonlinearity of asphalt concrete strain determined above shows the unacceptability of linear viscoelasticity theory and necessitates the use of nonlinear viscoelasticity theory.considerably in a nonlinear way at the given temperature, and the nonlinearity degree of relationship between stress and strain has been increased with the stress and strain duration increase.Nonlinearity of strain for asphalt concrete has been noted by other researchers [37,38].It should be also noted that the isochronous creep curves are similar ones.The property of similarity gives an opportunity to obtain all other isochronous creep curves, if any of them is known.At present, it has been accepted in road science and practice that mechanical properties of asphalt concrete depend on the value and duration of the applied load and temperature, and this dependence can be described in a sufficiently accurate way with the use of linear viscoelasticity theory [39][40][41][42].The essential nonlinearity of asphalt concrete strain determined above shows the unacceptability of linear viscoelasticity theory and necessitates the use of nonlinear viscoelasticity theory.Graphs of experimental rheological parameter at various stresses, constructed according to the calculation results under formula (17) with the use of data of Table 4, are shown in Figure 9. Graphs of experimental rheological parameter at various stresses, constructed according to the calculation results under formula (17) with the use of data of Table 4, are shown in Figure 9. Graphs of experimental rheological parameter at various stresses, constructed according to the calculation results under formula (17) with the use of data of Table 4, are shown in Figure 9.It is clearly seen that experimental rheological parameters of the asphalt concrete at various stresses do not merge into one curve; most of them diverge considerably from each other, which also shows the essential nonlinearity of the asphalt concrete strain at the considered temperatures and stresses.
Creep curves of the asphalt concrete at minimum stress, equal in our case to 0.041 MPa, have been approximated by Equation ( 5) with the use of Abel's kernel and the following values have been obtained for creep parameters: α = 0.3, ε 0 = 0.058, δ = 0.0138.Experimental and calculated creep curves of the asphalt concrete at the stress are represented in Figure 10 for comparison, which shows clearly that experimental creep curve of the asphalt concrete has been satisfactorily approximated by the equation of linear viscoelasticity with the use of Abel's kernel.It is clearly seen that experimental rheological parameters of the asphalt concrete at various stresses do not merge into one curve; most of them diverge considerably from each other, which also shows the essential nonlinearity of the asphalt concrete strain at the considered temperatures and stresses.
Creep curves of the asphalt concrete at minimum stress, equal in our case to 0.041 MPa, have been approximated by Equation ( 5  Then the values of instantaneous strain ε m 0 are required to calculate under formula (20) at other stresses.For this purpose, we select one value of load duration, for example, 570 s.Under the following expression, we obtain the value of calculated rheological parameter: According to the expression (32), we find that k m (570) = 2.6745 for the selected load duration and at known values of creep parameters α and δ.We have values of creep strain, equal to 0.1759%, 0.2349%, 0.3632%, and 0.5741%, respectively from Table 4 at t = 570 for stresses 0.074 MPa, 0.111 MPa, 0.148 MPa, and 0.183 MPa.According to formula (20) we find the following values of instantaneous strain of the asphalt concrete equal to 0.0658%, 0.0878%, 0.1358%, and 0.2147% at the abovementioned stresses.
Now the calculated values of creep strain for the asphalt concrete at various stresses, except for minimum ones, can be determined under formula (21).It is found out that calculated values of creep strain in initial time moments at maximum stress equal to 0.183 MPa do not sufficiently converge with experimental ones.Therefore, we redetermine the value of instantaneous strain at this stress under formula (22), having included into the calculation the values of experimental creep strain and model rheological parameter at t = 90,150, 270, and 570 s.The new value of instantaneous strain at stress 0.183 MPa has been obtained, which is equal to 0.2037.Calculated and experimental creep curves of the asphalt concrete at all stresses in comparison are represented in Figure 11.As it is seen, the calculated curves converge with sufficient accuracy with experimental ones at all stresses.

Conclusions
• Using the schematic creep curves and isochronous curves, Yu.N. Rabotnov's isochronous creep curve method has been visually explained.The nonlinear integral equation has been shown for mathematical description of the nonlinear deformation process for hereditary materials, proposed by Yu.N. Rabotnov.

•
Relevant equations have been determined from the nonlinear integral equation of Yu.N. Rabotnov for the application cases of Rabotnov's fractional exponential kernel and Abel's kernel for nonlinear deformation of hereditary materials at creep.The improved methods have been given for determination of creep parameters α, ε 0 , δ, β, and λ.

•
The detailed method has been developed for description of the nonlinear deformation process for hereditary materials.The notions have been introduced for experimental and model rheological parameters and similarity coefficients of isochronous curves.It has been shown how using them, one can find instantaneous strains at various stress levels for description of nonlinear deformation of hereditary materials at creep.

•
By processing and using test results for material Nylon 6 and glass-reinforced plastic TC 8/3-250, the process has been shown for sequential implementation of the developed methods for description of linear and nonlinear deformation of these materials at creep.The accuracy of the proposed methods is high.

•
By the results of experimental investigation, performed by the authors of this paper, it has been determined that fine-grained dense asphalt concrete at the temperature of 20 ± 2 • C and stresses up to 0.183 MPa at direct tension is deformed considerably in a nonlinear way.It has been shown in an illustrative way by construction of isochronous creep curves at various load durations and curves of experimental rheological parameter at various stresses.Nonlinear deformation of the asphalt concrete at creep is adequately described by the proposed methods.

Figure 1 .
Figure 1.Creep curves of a material at various constant stresses.

Figure 2 .
Figure 2. Isochronous deformation curves of a material.

Figure 1 . 18 Figure 1 .
Figure 1.Creep curves of a material at various constant stresses.

Figure 2 .
Figure 2. Isochronous deformation curves of a material.

Figure 2 .
Figure 2. Isochronous deformation curves of a material.

(Figure 3 .
Figure 3. Curves of experimental rheological parameter Ke(t) at various stresses σ: (a) for physically linear material; (b) for physically nonlinear material.

Figure 3 .
Figure 3. Curves of experimental rheological parameter K e (t) at various stresses σ: (a) for physically linear material; (b) for physically nonlinear material.
) at σ = 5 MPa we find that 0 m  = 0.1682%.Considering this fact, we determine the model rheological parameter:

Figure 7 .
Figure 7. Creep curves of the asphalt concrete at various stresses.

Figure 7 .
Figure 7. Creep curves of the asphalt concrete at various stresses.

Figure 8 .
Figure 8. Isochronous creep curves of the asphalt concrete.

Figure 8 .
Figure 8. Isochronous creep curves of the asphalt concrete.

Figure 8 .
Figure 8. Isochronous creep curves of the asphalt concrete.

Figure 9 .
Figure 9. Graphs of experimental rheological parameter of the asphalt concrete at various stresses.Figure 9. Graphs of experimental rheological parameter of the asphalt concrete at various stresses.

Figure 9 .
Figure 9. Graphs of experimental rheological parameter of the asphalt concrete at various stresses.Figure 9. Graphs of experimental rheological parameter of the asphalt concrete at various stresses.
) with the use of Abel's kernel and the following values have been obtained for creep parameters: α = 0.3, ε0 = 0.058, δ = 0.0138.Experimental and calculated creep curves of the asphalt concrete at the stress are represented in Figure10for comparison, which shows clearly that experimental creep curve of the asphalt concrete has been satisfactorily approximated by the equation of linear viscoelasticity with the use of Abel's kernel.

Figure 10 .Figure 10 .
Figure 10.Experimental and calculated creep curves of the asphalt concrete at stress 0.041 MPa.
Appl.Sci.2018, 8, x 16 of 18 Calculated and experimental creep curves of the asphalt concrete at all stresses in comparison are represented in Figure 11.As it is seen, the calculated curves converge with sufficient accuracy with experimental ones at all stresses.

Figure 11 .
Figure 11.Experimental and calculated creep curves of the asphalt concrete at various stresses.


Using the schematic creep curves and isochronous curves, Yu.N. Rabotnov's isochronous creep curve method has been visually explained.The nonlinear integral equation has been shown for mathematical description of the nonlinear deformation process for hereditary materials, proposed by Yu.N. Rabotnov. Relevant equations have been determined from the nonlinear integral equation of Yu.N. Rabotnov for the application cases of Rabotnov's fractional exponential kernel and Abel's kernel for nonlinear deformation of hereditary materials at creep.The improved methods have been given for determination of creep parameters α, ε0, δ, β, and λ. The detailed method has been developed for description of the nonlinear deformation process for hereditary materials.The notions have been introduced for experimental and model rheological parameters and similarity coefficients of isochronous curves.It has been shown how using them, one can find instantaneous strains at various stress levels for description of

Figure 11 .
Figure 11.Experimental and calculated creep curves of the asphalt concrete at various stresses.

Table 1 .
Creep strain values of material Nylon 6.

Table 4 .
Values of creep strain for asphalt concrete.