Research Article  Open Access
Determination of the Creep Parameters of Linear Viscoelastic Materials
Abstract
Creep process of linear viscoelastic materials is described by the integral equation of BoltzmannVolterra in which creep kernel is approximated by Rabotnov’s fractional exponential function. The creep equation contains four unknown parameters: , singularity parameter; , fading parameter; , rheological parameter; and , conditionally instantaneous strain. Twostage determination method of creep parameters is offered. At the first stage, taking into account weak singularity properties of Abel’s function at the initial moment of loading, parameters and are determined. At the second stage, using already known parameters and , parameters and are determined. Analytical expressions for calculating these parameters are obtained. An accuracy evaluation of the offered method with using experimentally determined creep strains of material Nylon 6 and asphalt concrete showed its high accuracy.
1. Introduction
Many materials used in technique and production (rocks, soils, metals, concretes, polymers, composites, etc.) under loading show viscoelastic properties. Currently a number of methods for quantitative evaluation of viscoelastic properties of materials are known [1–4]. One of the most frequently used methods is creep testing under uniaxial stress [5]. In the mathematical description of creep process of linear viscoelastic material the integral equation of BoltzmannVolterra with corresponding creep kernel is widely used. A large number of analytical expressions for the creep kernel are known [1, 4, 5].
In testing samples of viscoelastic materials on creep, usually it is difficult to measure instantaneous elastic strain . This is due to the fact that in the initial time after application of load to sample the rate of creep is extremely high. Its accurate measurement is practically impossible due to the dynamic effects caused by inertia of loading devices and measuring equipment.
Creep process of many materials particularly at small loads and low temperatures proceeds for a very long time. At high times of loading the rate of creep starts to fade and strain slowly approaching asymptotic value.
Thus, we can assume that creep kernel of integral equation must have the following two properties: firstly, having a weak singularity at the initial time point and, secondly, having a property of exponential function at high times of loading.
The most universal description of viscoelastic properties of materials satisfying the above requirements is a kernel in the form of fractional exponential function of Rabotnov [4–7]. Fractional exponential function is well studied and for simplification of calculations by using it, a special table has been developed [8]. Another advantage of the fractional exponential function is that for known creep kernel parameters constructed with using the fractional exponential function a kernel of relaxation becomes definite [9, 10].
Some methods for determining the parameters of Rabotnov’s creep kernel are known: using MittagLeffler’s function [11, 12], using LaplaceCarson’s transform [9, 12–14], and a method of direct approximation [9, 12–15].
In this paper we offer a new double stage method for determining the creep kernel parameters of linear viscoelastic materials constructed with using Rabotnov’s fractional exponential function. At the first stage, taking into account weak singularity properties of Abel’s function at the initial moment of loading, parameters and are determined. At the second stage, using already known parameters and , parameters and are determined.
Fractional derivative constitutive models for finite deformation of viscoelastic materials developed with using a kernel close to Abel’s kernel are given in [16]. A new model for the linear theory of generalized thermoviscoelasticity using a relaxation kernel close to Abel’s kernel is demonstrated in [17].
2. Creep Equation and Kernel
As is known generally creep curve of the viscoelastic material has three characteristic sites: I, site with a decreasing rate of creep (unstabilized creep), II, site with a constant (minimum) rate of creep (stabilized creep), and III, site of accelerating creep, preceding failure (Figure 1).
In the present work sites I and II of creep curves of linear viscoelastic materials under uniaxial tension are described by the integral equation of BoltzmannVolterra:where is strain at time point , , are stresses at time points and , is creep kernel, is instant modulus of elasticity, is observation time, and is time, preceding observation time.
Taking into account that for creep process , from (1) we will findwhere is conditionally instantaneous strain.
In (2) it was accounted that on Hooke’s law .
As mentioned above, creep kernel recorded using Rabotnov’s fractional exponential function in the following form [6]:where is Rabotnov’s fractional exponential function, , , and are creep kernel parameters (, , and ) and is gammafunction.
Substituting expression (3) for creep kernel in the creep equation (2) we obtain the following equation:
Equation (4) represents creep equation recorded using Rabotnov’s fractional exponential function. It contains four unknown parameters: is the singularity parameter; is the attenuation parameter; is the rheological parameter; and is conditionally instantaneous strain.
3. Method of Determining the Creep Parameters
3.1. Parameters and
Taking from series (4) we find the first term:
As one can see, obtained expression contains wellknown Abel’s function with parameter . Abel’s function at has a singularity at order . Based on this property of Abel’s function, unknown parameters and are determined by use of (5).
Knowing that parameter takes value from interval (0, 1) we will consider as unknown only parameters and .
According to the least squares method the best values of parameters and are those for which the following condition is met:where is sum of squares of deviations, is values of creep strain determined experimentally, and is number of creep strains.
From two equations based on expressions and , we find expressions for determining the parameters and :
Setting values of parameter from interval (0, 1) with a certain step from expression (7) we will find values of parameter . Substituting the found values of parameter and the corresponding values of singularity parameter in expression (8) values of parameter are determined.
Further, sequentially substituting values of singularity parameter and calculated corresponding values of parameters and in expression (5), values of creep strain are calculated.
If we designate an average deviation of calculated values of creep strains from experimental values through , as a selecting criterion of the best values of parameters , , and the following condition can be accepted:
3.2. Parameters and
Rewrite equation of creep (4) in the following form:
According to the least squares method analogically to condition (6) we write an extremum condition using the equation of creep (10):
From two equations based on the expressions and , expressions for determining the parameters and are found: but series is determined by expression (11).
Value of parameter of fading is determined from (13) by substituting trial values. If (13) has a unique solution, it will be desired value of the parameter . Obviously, for determining the values of parameters and from expressions (13) and (14) earlier calculated values of parameters and are used.
4. Experimental Approbation
4.1. Material Nylon 6
In work [9] the creep curve of Nylon 6 at a stress of MPa was approximated by expressionwhere [%] is strain and [h] is time.
Firstly, for determination of parameters and we calculate some strain values in time interval = 0,01⋯0,9 h. Time values and corresponding creep strain values are given in the first and the second columns of Table 1.
Substituting values of time , strain , and various values of singularity parameter from interval 0,05⋯0,95 with step of 0,05 in expressions (7) and (8) and taking into account condition (9) it has been found that ; ; .
With consideration of found values of parameters , , and , (5) constructed using Abel’s function for describing strain in site I of creep curve takes the form
Calculated by (17) values of creep strain and corresponding approximation errors are presented in the third and the fourth columns of Table 1.
As can be seen from Table 1, (17) with parameters and found by the offered method has high accuracy of approximation.
Since the parameter characterizes rate of material creep strain at long times, for definition it and parameter consider creep process of material Nylon 6 in the time interval = 0,01⋯100 h. Time values and corresponding values of creep strain calculated according to expression (16) are given in the first and the second columns of Table 2.

Using values of time , strain , parameters , , and trial values of parameter , from expression (13) it has been found that . Further from expression (14) we will define .
After substitution values of the found parameters , , , and (10) and auxiliary function (11) will have the next form:
Calculated by (18) values of creep strain and corresponding errors of approximation are given in the last two columns of Table 2.
As can be seen from Table 2 creep (18) with parameters , , , and defined by the offered method has high accuracy of approximation.
4.2. Asphalt Concrete
The strains calculated by approximation formula (16) were taken for determination of the above creep parameters of material Nylon 6 instead of strains obtained experimentally. The reason for this was the fact that in work [9] the creep strains were given only graphically in the form of a creep curve and numerical values of the strain are not shown. It appeared almost impossible to find the strain values from the graphical creep curve with sufficient accuracy especially at short times of loading.
Possibly, in order to reduce the length of material in published works containing an approximation of creep curves of various materials, as a rule, information on creep strain is given graphically, for example, in works [9, 10, 12–14, 18–20].
Further, we will show determination of creep parameters and by the above offered method (expressions (7) and (8)) using experimentally determined values of strain at site I of creep of hot finegrained asphalt concrete. This type of asphalt is used for application on top layer of road surface in many countries, including Kazakhstan. Asphalt concrete mix that meets the requirements of the Kazakhstan standard ST RK 1225 [21] was prepared with use of viscous road bitumen of grade BND 100/130 from Pavlodar petrochemical plant, which meets the standard requirements of ST RK 1373 [22]. Samples of asphalt concrete with rectangular crosssection dimensions of 5 × 5 × 15 cm were prepared by means of Cooper compactor (UK, model CRTRC2S) in accordance with European Union standard EN 1269733 [23]. Detailed information on preparation of asphalt mixtures and samples, as well as standard properties of bitumen and asphalt concrete, is given in [24]. Samples of asphalt concrete were tested for creep according to the scheme of uniaxial tension at a temperature of °C and different loads.
Isochronous lines of deformation for the asphalt concrete constructed using creep curves within site I at durations of loadings 20, 40, and 60 s are shown in Figure 2. As it can be seen dependence between stress and strain is linear; that is, asphalt concrete under these conditions is linearly deformable material and for description of its creep process the equation of BoltzmannVolterra is applicable.
Figures 3 and 4 show site I of creep curves of asphalt concrete at stresses of 0.16 and 0.19 MPa, respectively. In these figures points designate the creep strain defined experimentally and lines designate its approximation by Abel’s kernel (5). It is seen that compliance of approximating curves to the experimental data is quite good. Average deviations of calculated strains from experimental creep curves at 0,16 and 0,19 MPa are 1.5% and 0.5%, respectively. The values of parameters and found by the offered method are equal to and at MPa and and at MPa.
Thus, the integral equation of BoltzmannVolterra with creep kernel constructed using Rabotnov’s fractional exponential function with parameters defined by the offered method describes the creep process of linear viscoelastic materials with high accuracy.
5. Conclusions
(1)Describing the creep process of linear viscoelastic materials by the integral equation of BoltzmannVolterra with creep kernel which is constructed using Rabotnov’s fractional exponential function is suggested.(2)A new method of determining the creep kernel parameters of linear viscoelastic materials constructed using Rabotnov’s fractional exponential function is offered. It is implemented in two stages. At the first stage based on the property of singularity of Abel’s function at the initial moment of time parameters and are defined. At the second stage using already known parameters and parameters and are defined. Analytical expressions for calculating these parameters are obtained.(3)Applicability of the offered method is evaluated using data on creep strain of material Nylon 6 and hot finegrained asphalt concrete obtained experimentally. The results showed high accuracy of the offered method.
Competing Interests
The authors declare that they have no competing interests.
References
 R. M. Cristensen, Theory of Viscoelasticity: An Introduction, Academic Press, New York, NY, USA, 1971.
 N. W. Tschoegl, The Phenomenological Theory of Linear Viscoelastic Behavior: An Introduction, Springer, Berlin, Germany, 1989. View at: Publisher Site  MathSciNet
 J. D. Ferry, Viscoelastic Properties of Polymers, John Wiley & Sons, New York, NY, USA, 3rd edition, 1980.
 Yu. N. Rabotnov, Mechanics of Deformed Solid Body, Nauka, Moscow, Russia, 1988.
 Yu. N. Rabotnov, Creep of Structure Elements, Nauka, Moscow, Russia, 1966.
 Y. N. Rabotnov, Elements of Hereditary Mechanics of Solids, Nauka, Moscow, Russia, 1977. View at: MathSciNet
 Yu. N. Rabotnov, “Balance of the elastic medium with a aftereffect,” Applied Mathematics and Mechanics, vol. 12, pp. 53–62, 1948. View at: Google Scholar
 Y. N. Rabotnov, L. Kh. Papernik, and E. N. Zvonov, Tables of a Fractional Exponential Function of Negative Parameter and Its Integral, Nauka, Moscow, Russia, 1969.
 Yu. V. Suvorova, “About Yu.N. Rabotnov's nonlinear hereditary equation and its applications,” News of the Russian Academy of Sciences, Mechanics of solids, no. 1, pp. 174–181, 2004. View at: Google Scholar
 Y. V. Suvorova, “Master equations for viscoelastic materials taking into account the influence of temperature and humidity,” Strength of Materials, vol. 32, no. 5, pp. 490–496, 2000. View at: Publisher Site  Google Scholar
 M. I. Rozovskii, “Some features of elasticity hereditary media,” News of the USSR Academy of Sciences, no. 2, pp. 30–36, 1961. View at: Google Scholar
 V. P. Golub, P. V. Fernati, and Y. G. Lyashenko, “Determining the parameters of the fractional exponential heredity kernels of linear viscoelastic materials,” International Applied Mechanics, vol. 44, no. 9, pp. 963–974, 2008. View at: Publisher Site  Google Scholar
 E. N. Zvonov, N. I. Malinin, L. Kh. Papernik, and B. M. Tseitlin, “Determination of the creep characteristics of linear elastic hereditary materials by using a digital computer,” Mechanics of Solid, no. 5, pp. 76–82, 1968. View at: Google Scholar
 Yu. V. Suvorova and A. V. Mosin, “Determination of parameters of the Rabotnov's fractional exponential function with use of integral transform and modern software,” Problems of Mechanical Engineering and Automation, vol. 4, pp. 54–56, 2002. View at: Google Scholar
 D. A. Gavrilov and V. A. Markov, “Numerical method of determining the rheologic parameters of composites from test results,” Mechanics of Composite Materials, vol. 22, no. 4, pp. 417–421, 1987. View at: Publisher Site  Google Scholar
 M. Fukunaga and N. Shimizu, “Fractional derivative constitutive models for finite deformation of viscoelastic materials,” Journal of Computational and Nonlinear Dynamics, vol. 10, no. 6, Article ID 061002, 8 pages, 2015. View at: Publisher Site  Google Scholar
 M. A. Ezzat, A. S. ElKaramany, and A. A. ElBary, “Thermoviscoelastic materials with fractional relaxation operators,” Applied Mathematical Modelling, vol. 39, no. 2324, pp. 7499–7512, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 R. D. Maksimov, L. Jirgens, J. Jansons, and E. Plume, “Mechanical properties of polyester polymerconcrete,” Mechanics of Composite Materials, vol. 35, no. 2, pp. 99–110, 1999. View at: Publisher Site  Google Scholar
 R. D. Maksimov and E. Plume, “Longterm creep of hybrid aramid/glassfiberreinforced plastics,” Mechanics of Composite Materials, vol. 37, no. 4, pp. 271–280, 2001. View at: Publisher Site  Google Scholar
 E. A. Sokolov and R. D. Maksimov, “Possibilities of predicting the creep of a polymeric fiberreinforced plastic from the properties of the components,” Polymer Mechanics, vol. 14, no. 6, pp. 807–813, 1978. View at: Publisher Site  Google Scholar
 ST RK 12252003, “Asphalt concrete mixes for road, airfield and asphalt concrete,” Technical specifications, Astana, Kazakhstan, 2003. View at: Google Scholar
 “Bitumen and bituminous binders. Oil road viscous bitumens,” ST RK 13732005, Technical Specifications, Astana, Kazakhstan, 2005. View at: Google Scholar
 European Committee for Standardization, “Bituminous mixtures. Test methods for hot mix asphalt —part 33: specimen prepared by roller compactor,” Tech. Rep. EN 1269733, 2003. View at: Google Scholar
 B. Teltayev, “A new failure criterion for asphalt mixtures under fatigue loading,” International Journal of Pavement Research and Technology, vol. 8, no. 4, pp. 276–282, 2015. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Alibay Iskakbayev et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.