Next Article in Journal
Estimation of the Band Gap of Carbon Nanotube Bundles
Next Article in Special Issue
Direct Identification of the Continuous Relaxation Time and Frequency Spectra of Viscoelastic Materials
Previous Article in Journal
Utilisation of Waste Sludge from Drinking Water Treatment as a Filler Material in Hot Mix Asphalt
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sampling Points-Independent Identification of the Fractional Maxwell Model of Viscoelastic Materials Based on Stress Relaxation Experiment Data

Department of Technology Fundamentals, Faculty of Production Engineering, University of Life Sciences in Lublin, 20-612 Lublin, Poland
Materials 2024, 17(7), 1527; https://doi.org/10.3390/ma17071527
Submission received: 8 March 2024 / Revised: 21 March 2024 / Accepted: 24 March 2024 / Published: 27 March 2024

Abstract

:
Considerable development has been observed in the area of applying fractional-order rheological models to describe the viscoelastic properties of miscellaneous materials in the last few decades together with the increasingly stronger adoption of fractional calculus. The fractional Maxwell model is the best-known non-integer-order rheological model. A weighted least-square approximation problem of the relaxation modulus by the fractional Maxwell model is considered when only the time measurements of the relaxation modulus corrupted by additive noises are accessible for identification. This study was dedicated to the determination of the model, optimal in the sense of the integral square weighted model quality index, which does not depend on the particular sampling points applied in the stress relaxation experiment. It is proved that even when the real description of the material relaxation modulus is entirely unknown, the optimal fractional Maxwell model parameters can be recovered from the relaxation modulus measurements recorded for sampling time points selected randomly according to respective randomization. The identified model is a strongly consistent estimate of the desired optimal model. The exponential convergence rate is demonstrated both by the stochastic convergence analysis and by the numerical studies. A simple scheme for the optimal model identification is given. Numerical studies are presented for the materials described by the short relaxation times of the unimodal Gauss-like relaxation spectrum and the long relaxation times of the Baumgaertel, Schausberger and Winter spectrum. These studies have shown that the appropriate randomization introduced in the selection of sampling points guarantees that the sequence of the optimal fractional Maxwell model parameters asymptotically converge to parameters independent of these sampling points. The robustness of the identified model to the measurement disturbances was demonstrated by analytical analysis and numerical studies.

1. Introduction

For several decades, fractional-order rheological models have been used to describe, analyze and improve the viscoelastic properties of different materials. In addition to theoretical research dedicated to fractional-order rheological models [1,2,3,4,5], hundreds of studies have been conducted on the applicability of such models for specific materials to describe their mechanical properties. The applicability of such models to the description of different polymers is well known, for example, poly-isobutylene [4], polyurea and PET [6], shape memory polymers [7], amorphous polymers [8] and flax fiber-reinforced polymer [9]. Fractional viscoelastic models are also used for modeling laminated glass beams in the pre-crack state under explosive loads [10]; stress relaxation behavior of glassy polymers [11]; description of fiber-reinforced rubber concrete [12]; viscoelastic modeling of modified asphalt mastics [13]; and modeling rate-dependent nonlinear behaviors of rubber polymers [14]. The modeling and simulation of viscoelastic foods, for example, food gums [15], carrot root [16], fish burger baking [17], is another field of application of rheological fractional models. Due to the non-integer order of the operations of integration and differentiation, the fractional-order models have improved flexibility and better adjustment to material characteristics, both in the time and frequency domains, compared to those of the classic integer-order models.
Although over the last several decades different fractional differential models have been proposed for modeling the viscoelastic processes in materials, the fractional Maxwell model (FMM) is the best known [4,5]. The relaxation modulus of the FMM, described by the product of Mittag-Leffler and inverse power functions, allows for the modeling of a very wide range of stress relaxation processes in materials. Describing the rheological properties of polymers by the FMM [18,19] is well known. However, the FMM was also applied, for example, for computational modeling and analysis on the damping and vibrational behaviors of viscoelastic composite structures [20], viscoelastic flow in a circular pipe [21], effect of temperature on the dynamic properties of mixed surfactant adsorbed layers at the water/hexane interface [22,23] and constitutive equations of the Mn-Cu damping alloy [24]. Fractional viscoelasticity described by the Maxwell model turned out to model both exponential and non-exponential relaxation phenomena in real materials.
Different identification methods for the recovery of the parameters of the non-integral-order models, including the FMM, from both static [16,25,26,27,28] and dynamic [12,29,30,31] experiments data have been proposed so far. It is known that different identification methods in association with different experiment plans result in different identification data yield models, which may differ [32]. Generally, the identification result, i.e., the chosen model, is influenced by the three entries that are necessary for model identification: the set of models from which the best model is chosen, the rule for the optimal model selection and the measurement data obtained in the experiment [32,33]. For the selected class of models, here, the set of the fractional Maxwell models, the identified model depends on the identification rule and the experiment data. The model parameters are usually determined by guaranteeing the “best-possible” fit to the measurements. Therefore, parameters of the optimal model are dependent on the measure applied for evaluating the “best” [32]. The mean-square approximation error is the predominant selection of the model quality measure, which results in a standard least-squares identification task. For the selected identification index, the model identified is usually dependent, sometimes even very strongly, on the experiment data. This is the case with FMM identification methods known in the literature [12,16,25,26,27,28,29,30,31]. This paper deals with the problem of the FMM identification using measurement data from the stress relaxation test. Therefore, the sampling instants used in the experiment and discrete-time measurements of the relaxation modulus compose the set of the experiment data. To build the optimal fractional Maxwell model whose parameters do not depend on sampling instants applied in the stress relaxation test is the aim of this paper.
In the previous paper [33], the problem of the least-squares approximation of the relaxation modulus has been considered for an assumed wide class of relaxation modulus models. Models being continuous, differentiable and Lipschitz continuous with respect to the parameters have been assumed. The main results in [33] refer to the models that are determined asymptotically, when the number of measurements tend to infinity. Whenever some applicability conditions concerning the chosen class of models are satisfied, the asymptotically optimal FMM parameters can be determined using the measurement data obtained for sampling instants selected randomly due to the appropriate randomization, even when the true relaxation modulus description is completely unknown. For the exponential Maxwell and the exponential stretched Kohlrausch–Williams–Watts models, the applicability conditions are satisfied [33]. It should be noted that the concept of identification being measurement point-independent comes from the Ljung paper [34] and the paper of [35], in which the optimal identification problems for dynamic and static systems have been considered.
In this paper, the concept of introducing an appropriate randomization for the selection of sampling instants at which the measurements of the relaxation modulus are recorded is applied for the fractional Maxwell model identification. Following [33], the real material description is completely unknown and only the measurement data of the relaxation modulus are available for model identification. Identification consists of determining the FMM that solves the problem of an optimal least-squares approximation of a real relaxation modulus. The complicated form of the relaxation modulus of the FMM (the product of Mittag-Leffler and inverse power functions) implies that the applicability of the sampling points-independent identification for FMM identification is not obvious. It is known that the relaxation modulus of the FMM is continuous and differentiable with respect to its four parameters [36]. However, the satisfaction of the Lipschitz continuous property with the bounded Lipschitz constant is proved in this paper for the first time, to guarantee the applicability of the experiment randomization concept.
A complete identification scheme leading to the strongly consistent estimate of the optimal model was specified. Assuming that the measurements are corrupted by additive disturbances, the stochastic-type analysis of the model convergence was carried out, and the exponential rate of convergence was demonstrated both analytically and by numerical studies. For materials described by the unimodal Gauss-like spectrum of relaxation used to describe the rheological properties of the materials [37,38,39] and by the Baumgaertel, Schausberger and Winter (BSW) spectrum [40,41] successfully applied for modeling the polymers [42,43], based on the simulation experiments, both the asymptotic properties and noise robustness of the algorithm were numerically studied. To improve the clarity of this article, the proof of the new FMM Lipschitz property is moved to Appendix A. The tables with the results of the numerical studies are given in Appendix B.

2. Materials and Methods

2.1. Material

A linear viscoelastic material subjected to small deformations for which the uniaxial, non-aging and isotropic stress–strain equation is given by a Boltzmann superposition equation [44]
ς t = t G t τ d ε τ d τ d τ
is considered, where ς t and ε t are, respectively, the stress and strain and G t denotes the linear relaxation modulus. By Equation (1), the stress ς t at time t depends on the earlier history of the strain rate described by the first-order derivative d ε τ d τ via the kernel given by the relaxation modulus G t .
The modulus G t is the stress induced in the material described by constitutive Equation (1) by the unit step strain ε t imposed. It is assumed for the studied material that the mathematical description of the modulus G t is completely unknown. However, the real relaxation modulus G t is accessible by measurement with a certain accuracy for an arbitrary time t T . Here, T = t 0 , T with the initial time t 0 > 0 and T .
We make the following assumption [33]:
Assumption 1.
The relaxation modulus  G t  of the material is bounded on  T , i.e.,  sup t T G t M < .

2.2. Fractional Maxwell Model

Constitutive equation of the fractional order Maxwell model is as follows [2,4,45]:
τ r α β d α β ς t d t α β + ς t = G e τ r α d α ε t d t α ,
where G e denotes the elastic modulus, τ r means the relaxation time, α and β are non-integer positive orders of fractional derivatives of the strain ε t and stress ς t , respectively. In this paper, d α d t α f ( x ) = D t α f ( x ) means the fractional derivative operator in the sense of Caputo’s of a function f ( x ) of non-integer-order α with respect to variable t and with a starting point at t = 0 , which is defined by [1,4]
D t α f t = 1 Γ n α 0 t t 1 n α 1 d n d t n f t d t ,
where n 1 < α < n and Γ n is Euler’s gamma function [1] (Equation (A.1.1)).
The FMM (2) can be considered as a generalization of the classic viscoelastic Maxwell model being the series connection of the ideal spring with a dashpot (see Figure 1a) described by a differential equation of the first order [44,46]:
d ς t d t + 1 τ r ς t = G e d ε t d t ,
with the elastic modulus G e of the spring, the relaxation time τ r = η / G e , where η means the viscosity of the dashpot.
A series connection (see Figure 1c), analogical to the classic Maxwell model, of two elementary fractional Scott-Blair elements G e 1 , τ r 1 , α and G e 2 , τ r 2 , β , both described by the fractional differential equation of the general form [2,4,45]
ς t = G e τ r α d α ε t d t α ,
with the parameters G e , τ r , α (see Figure 1b), yields the FMM described by Equation (2), where the parameters G e 1 , τ r 1 , α and G e 2 , τ r 2 , β uniquely determine the parameters E and τ r of the FMM (2); for details, see [16]. The four parameters G e , τ r , α , β of the FMM (2), compared with only two parameters G e , η , or equivalently G e , τ r of the classic Maxwell model (3), are important for the improvement in the FMM accuracy and flexibility.
The uniaxial stress response of the FMM (2) imposed by the unit step strain ε t , i.e., the time-dependent relaxation modulus G t , for an arbitrary 0 < β < α 1 is given by the formula [2,4,5]:
G t = G e t τ r β E α β , 1 β t τ r α β ,
where E κ , μ x is the generalized two-parameter Mittag-Leffler function defined by series being convergent in the whole z-complex plane [1,2]:
E κ , μ x = n = 0 x n Γ κ n + μ .
Further, for the description of the FMM identification task, relaxation modulus model (5) is denoted as
G M t , g = G e t τ r β E α β , 1 β t τ r α β ,
to emphasize the dependence on a four-element vector of model parameters
g = α β G e τ r T ,
where the subscript ‘ M ’ means the model.
For the special case α = β , the FMM (2) reduces to the Scott-Blair model (compare (4))
2 ς t = G e τ r α d α ε t d t α ,
and the relaxation modulus is described by
G M t , g = G e 2 Γ 1 α t τ r α .
Let us consider the following set of the FMM admissible model parameters:
G = g : β 0   β α 1 ;   G e m i n G e G e m a x ; τ r m i n τ r τ r m a x
where β 0 > 0 is an arbitrarily small positive number and the maximal and minimal values of elastic modulus G e and relaxation time τ r follow from the a priori knowledge concerning the material under investigation and are such that G e m i n > 0 and τ r m i n > t 0 . G is a compact subset of the four-dimensional real space R 4 .
The properties of the two-parameter Mittag-Leffler function and the model (7) have been studied by many authors [1,2,3,4,5]. The function E κ , μ x (6) is completely monotonic on the negative real axis for 0 < κ 1 and μ κ , i.e., the function E κ , μ x is completely monotonic for x > 0 , Ref. [4] (Equation (E.32)). Whence, since t 0 > 0 , by virtue of (6), for any t T , and any g G , we have
E α β , 1 β t τ r α β < E α β , 1 β 0 = 1 Γ 1 β 1 .
Let us introduce the function [4] (Equation (E.53))
e κ , μ x ; λ = x μ 1 E κ , μ λ x κ ,
which, comparing (7) and (13), enables describing the relaxation modulus G M t , g (7) in compact form as follows
G M t , g = G e τ r β e α β , 1 β t ; τ r β α .
The function e κ , μ x ; λ (13) is known to play a crucial role in many problems of fractional calculus [4] (p. 372) because it has many excellent and useful properties; some of them were used in this paper. The function e κ , μ x ; λ is completely monotonic for x > 0 when 0 < κ μ 1 whenever the parameter λ > 0 [4] (p. 373) as the product of two completely monotonic functions, which by (14) implies the complete monotonicity of the relaxation modulus model G M t , g for t > 0 whenever 0 β < α 1 . This means, in particular, that for t > 0 and g G , such that 0 < β < α 1 , the positive definite model G M t , g (7) monotonically decreases with increasing t > 0 . Therefore, for any t > 0 and any g G , such that 0 < β < α 1 , in view of (12)–(14), we have
G M t , g G e τ r β t 0 β E α β , 1 β t 0 τ r α β G e m a x m 0 ,
where m 0 is defined below by the sequence of inequalities valid for any t T and any g G
t τ r β τ r m a x t 0 β τ r m a x t 0 = m 0 ,
where t 0 > 0 .
For the case α = β , the relaxation modulus G M t , g (10) is also a completely monotonic function of the time for t > 0 , which in view of (16) is uniformly bounded for t T and g G by G e m a x m 0 / 2 .
Therefore, there exists a positive constant M 1 = G e m a x m 0 such that
sup t T , g G G M t , g M 1 < ,
i.e., the modulus G M t , g is uniformly bounded on the set T × G .
Inequality (17) combined with Assumption 1 implies the upper bound
sup t T , g G G t G M t , g M + M 1 < .
The Lipschitz continuity of the model G M t , g with respect to parameter g , which is not obvious, in particular, with respect to non-integer orders of fractional derivatives, is fundamental to guarantee the convergence of the optimal models for the applied here experiment randomization. Therefore, before the identification concept and the respective algorithm are presented, the Lipschitz property of the mapping G M t , g (7) will be proved, as summarized in the following theorems. The quite tedious proofs are moved into Appendix A.1.

2.3. Lipschitz Continuity of FMM with Respect to Model Parameters

Due to the relation between the parameter α and β , let us consider two cases separately when (a) β < α and (b) β = α . Therefore, the set of admissible model parameters G (11) is decomposed on two disjoint subsets:
G 1 = g : β 0   β < α 1 ;   G e m i n E G e m a x ; τ r m i n τ r τ r m a x
and
G 2 = g : β 0   β = α 1 ;   G e m i n E G e m a x ; τ r m i n τ r τ r m a x ,
in which the relaxation modulus G M t , g is described by the formulas (7) and (10), respectively. The bounded set G 1 is non-closed, i.e., the compactness property of the set G (11) is lost here, while G 2 is compact.
The following spectral representation derived in [47]
G M t , g = G e τ r α π 0 τ r v α β sin π β + sin π α τ r v 2 α β + 2 τ r v α β cos π α β + 1 v α 1 e t v d v ,
which results from the known spectral representation of the two-parameter Mittag-Leffler function [1] (Theorem 4.18, Equations (4.7.17) and (4.7.15)) and is valid for 0 < β < α 1 , will be used for g G 1 . Applying the differential approach in Appendix A.1, the next result is proved.
Theorem 1.
Let  G 1  defined by (19) be the set of the fractional Maxwell model admissible parameters. Then, the relaxation modulus  G M t , g  (7) of the FMM (2) is continuous and differentiable with respect to  g  (8) for any time  t T  and
sup t T , g G 1 g G M t , g 2 M 2 < ,
where  g G M t , g  denotes the gradient of the function  G M t , g  with respect to the vector  g ; here,  · 2  is the Euclidean norm in the space  R 4 .
The above theorem means, in particular, that for an arbitrary small positive β 0 , the mapping G M : T × G 1 R defined according to Equation (7) is, uniformly with respect to the time t T , a Lipschitz continuous function of the vector of model parameters g with Lipschitz constant M 2 .
In the case (b) β = α , for the set of model parameters G 2 (20), the FMM (2) is described by the power-law relaxation modulus G M t , g (10) and the absolute boundness of the gradient g G M t , g is resolved by the next result proved in Appendix A.2.
Theorem 2.
Let  G 2 , defined by (20), be the set of the fractional Maxwell model admissible parameters with equal parameters α  and  β . Then, the relaxation modulus  G M t , g  (10) of the model (9) is continuous and differentiable with respect to  g  (8) for any time  t T  and
sup t T , g G 2 g G M t , g 2 M 3 < .
From the proofs of the above theorems, especially from the nonnegative definiteness of the derivatives G M t , g E (A5), G M t , g τ r (A6) and the two last elements of the gradient g G M t , g (A53), the following property is derived.
Property 1.
Let  G  defined by (11) be the set of the FMM (2) admissible parameters. Then, for any fixed time  t T , the relaxation modulus  G M t , g  described by (7) or (10) monotonically increases with increasing parameters  G e  and  τ r  and other parameters being fixed, i.e., the greater parameters  G e  and  τ r  are, the greater the relaxation modulus  G M t , g  is for the given  t T .

2.4. Relaxation Modulus Measurements

Following [33,35], let T 1 , , T N be independent random variables with a common probability density function ρ t ; T is the support of ρ t . Then, let G i = G T i be the related relaxation modulus of the material for i = 1 , , N . Let G ¯ i denote their measurements corrupted by additive noise Z i , i.e., G ¯ i = G i + Z i , recorded in the stress relaxation experiment [44,46,48].
The two assumptions concerning the measurement noises are taken (compare Assumptions 5 and 6 in [33]) as follows:
Assumption 2.
The measurement noise  Z i  is a time-independent, i.e., independent of the variables  T i , sequence of independent identically distributed (i.i.d.) random variables with zero mean  E Z i = 0  and a common finite variance  E Z i 2 = σ 2 < .
Assumption 3.
The measurement noises  Z i  are bounded by  δ , i.e.,  Z i δ <  for  i = 1 , , N .
Both the above assumptions and Assumption 1, concerning the real relaxation modulus, are natural in the context of the relaxation modulus identification [33].

2.5. Identification Problem

FMM identification involves selecting from a given class of models defined by (7) and (10) the model that best fits the measurement data. Suppose an identification experiment resulted in a set of measurements G ¯ T i = G T i + Z i at the sampling times T i t 0 > 0 , i = 1 , , N . The mean-squares index
Q N g = 1 N i = 1 N G ¯ T i G M T i , g 2 ,
is taken as a measure of the FMM model accuracy. Here, the lower index denotes the number of measurements. Then, the problem of the optimal model identification consists of the solution of the minimization task
m i n g G     Q N g = Q N g N ,
where g N is the optimal model parameter. Since, due to the continuity of the model G M t , g with respect to the parameter g , the index Q N g is a continuous function of g and the set of admissible parameters G (11) is compact, the existence of the solution to the optimization problem (25) immediately results from the Weierstrass theorem about the extreme of continuous function on the compact set [49]. Since the minimum g N can be not unique, let G N denote the set of vectors g N that solve the optimization task (25).
The parameters g N of the identified relaxation modulus model G M t , g N are dependent on the measurement data, in particular, on the sampling instants T i . To make the model independent of specific sampling instants T i , the optimal sampling points-independent approximation problem is stated in the following subsection.

2.6. The Optimal FMM

Let us consider the following problem of determining such an FMM that minimizes the global approximation error:
Q g = T G t G M t , g 2 ρ t d t ,
where the selected weight function, such that 0 ρ t M 0 < , is a density on the set T , i.e., T ρ t d t = 1 .
The integral (26) is absolutely integrable, uniformly on G , both for the bounded or unbounded domain T as the product of a function G t G M t , g 2 , in view of (18) bounded uniformly for t , g T × G , and absolutely integrable function ρ t . Therefore, the integral (26) is well defined for any g G .
The problem of the optimal approximation of the real modulus G t within the class of the fractional Maxwell models relies on determining the parameter g that minimizes Q g over the set G , i.e., in solving optimization task
m i n g G   Q g = Q g .
Due to continuity of G M t , g with respect to the vector g , the index Q g (5) is a continuous function of g , and thus, the existence of the solution g follows from the previously mentioned Weierstrass theorem concerning the extreme of continuous function on the compact set. Let the set of model parameters g solving (27) be denoted by G . Any parameter g G does not depend on the specific time instants applied in the experiment.

3. Results and Discussion

In this section, the analysis of the asymptotic properties of the identified fractional Maxwell model, when the number of measurements tend to infinity, is conducted. The rate of the convergence of this model to the optimal FMM, which does not depend on the experiment data, is studied. The resulting identification algorithm is outlined. Next, the analytically proven properties of the identification method are verified by numerical simulations and studies. Two example materials are simulated. In the first, the “real” material is described by a unimodal Gauss-like relaxation spectrum [37,38,39] with short relaxation times and the Baumgaertel, Schausberger and Winter (BSW) spectrum [40,41] with long relaxation times. Such models are used to describe the rheological properties of various materials, especially polymers and biopolymers. Based on the noise-corrupted data from the simulated randomized stress relaxation experiment, the optimal FMM models are determined. The asymptotic properties and noise robustness have been studied.

3.1. Convergence

The empirical index Q N g (24) can be obtained by the replacement of the integral in Q g (26) with the finite mean sum of squares, which is clear from a practical point of view. For i = 1 , , N , by Assumption 2, the expected value is
E G T i + Z i G M T i , g 2 = Q g + σ 2 ,
whence, by (24), the expected value is
E Q N g = Q g + σ 2 .
To investigate the stochastic-type asymptotic properties of the empirical identification task given by (25), some properties derived in [35] will be used. Note, that Assumptions A1–A3 from [35], concerning the compactness of the set of model admissible parameters, continuity, differentiability and Lipshitzness of the model are satisfied here. Taken above, Assumption 2 is identical with Assumption A5 in [35], while property (18) also means that Assumption A4 from [35] is satisfied, i.e., all the assumptions from [35] hold here.
By (28), Property 2 from [35] implies the next result.
Property 2.
When Assumptions 1 and 2 are satisfied, then
s u p g G     Q g + σ 2 Q N g 0     w . p . 1     a s     N ,
where  w . p . 1  means “with probability one”.
By (28) and (29), the empirical identification index Q N g (24) is arbitrarily close to its expected value, uniformly in g over the set G . In consequence, the model parameter g N solving empirical identification task (25) can be related to the parameter g that solves the sampling points-independent minimization task (27). From the uniform in g G convergence of the index Q N g in (29), we conclude immediately the main result of this subsection, c.f., Assertion in [35] or Equation (3.5) in [34].
Property 3.
Assume that Assumptions 1 and 2 hold,   T 1 , , T N  are independently and randomly selected from   T , each according to the distribution with probability density function   ρ t . If the solutions to the minimization problems (25) and (27) are unique, then
g N g   w . p . 1     a s     N
and
G M t , g N G M t , g   w . p . 1     a s     N .
for all   t T . If the minimization problems (25) and (27) do not have unique solutions, then for any convergent subsequence of the sequence   g N , where   g N G N ,
g N G   w . p . 1     a s     N
and for any   t T and some   g G ,the convergence in (31) holds.
The existence of a convergent subsequence of g N so that the asymptotic property (32) holds results directly from the compactness of G (11). Therefore, under Assumptions 1 and 2, the optimal parameter g N of the FMM is a strongly consistent estimate of some parameter g G .
Since, by Theorems 1 and 2, the model G M t , g is Lipschitz on G uniformly in t T , then the almost-sure convergence of g N to g in (30) implies that, c.f., (Ref. [35]: Remark 2):
s u p t T     G M t , g N G M t , g 0     w . p . 1     a s     N .
i.e., that G M t , g N is a strongly consistent estimate of the optimal FMM G M t , g , uniformly on T .
Concluding, when Assumptions 1 and 2 are satisfied, the arbitrarily fine approximation of the FMM with the optimal parameter g can be determined (almost everywhere) as the number of measurements N grow enough, even if the real description of the material modulus is fully unknown.

3.2. Exponential Rate of Convergence

Analyzing the convergence in (30) and (32), the question immediately arises of how fast g N tends to some g G as N grows large. As in [35], the distance between the model parameters g N and g will be evaluated by means of the integral identification index Q g (26), i.e., in the sense of the difference Q g Q g N . For this purpose, it will be checked how fast, for a given small ε > 0 , the probability P Q g Q g N ε tends to zero, as N increases. From the well-known Hoeffding’s inequality [50], the upper bound of this probability can be derived, analogous to inequality (15) in [35] or inequality (22) in [33] (for details, see Appendix A.1 in [33]):
P Q g Q g N ε 2 e x p N ε 2 8 M ^ 2 ,
for any ε > 0 , where
M ^ = 2 M + M 1 2 + σ 2 + δ 2 + 2 M + M 1 δ ,
with the constants M and M 1 defined in Assumption 1 and Equation (17), respectively, the noises’ variance σ 2 and upper bound δ are introduced by Assumptions 2 and 3.
The inequality (34) describes the influence of the number of measurements N and the noises’ ”strength” on the rate of convergence. For ε being fixed, the bounds for P Q g Q g N ε decrease exponentially to zero as N increases. The convergence rate is the higher, the lower is M ^ (35). In particular, a quick inspection of (35) shows that for stronger measurement noises, the rate of convergence is reduced. Larger δ and σ yield a greater decrease in the rate. This is as expected, since with large disturbances, the measurements are not very adequate. Simultaneously, the larger M + M 1 , i.e., in view of the estimation (18), the greater the discrepancy between the real modulus and the FMM, the worse the convergence.

3.3. Identification Algorithm

In view of the convergence properties (30), (31) the computation of the parameter g N approximation the parameter g of the optimal FMM requires the next steps:
  • Select randomly from the set T the sampling times t 1 , , t N , choosing each t i independently, according to the probability distribution of the density ρ t defined given by the weight function in the integral Q g (26).
  • Conduct the stress relaxation experiment [44,46,48], measure and store the measurements G ¯ i of the relaxation modulus for the selected time instants t i , i = 1 , , N .
  • Solve the identification optimization task (25) and compute the identified model parameter g N .
  • Put N ¯ = N and g N ¯ = g N . To extend the set of experiment data, select new N N ¯ .
  • Repeat Steps 1–3 for a new N , that is, randomly choose new sampling times, conduct the rheological experiment once more for a new sample of the material and determine the next g N .
  • Examine if g N ¯ g N 2 < ε , where ε is a small positive number, to check if g N ¯ is an adequate approximation of g . If yes, stop the scheme and take g N ¯ as the approximate value of g . Otherwise, go again to Step 4.
Remark 1.
A less restrictive testing regarding whether   Q N ¯ g N ¯ Q N g N < ε   holds can be used as an alternative for the stopping rule from Step 6. Both types of stopping rules are commonly used in numerical optimization techniques.

3.4. Numerical Studies

The results of the numerical studies are concerned with the asymptotic properties of the determined optimal FMM and the influence of the measurement noises on this model. Apart from the theoretical analysis above, these simulation studies make it possible to show the respectability and effectiveness of the method developed for FMM identification.
Firstly, it is assumed that the rheological properties of the material are described by the Gaussian-like distribution of the relaxation spectrum, which were used to represent the viscoelastic properties of numerous materials, e.g., polyacrylamide gels [48], native starch gels [38], glass [39], poly(methyl methacrylate) [37], polyethylene [51] and carboxymethylcellulose (CMC) [52]. The spectra of various biopolymers determined by many researchers are Gaussian in nature, for example, cold gel-like emulsions stabilized with bovine gelatin [53], fresh egg white-hydrocolloids [52], some (wheat, potato, corn and banana) native starch gels [38], xanthan gum water solution [52] and wood [54,55].
Next, it is assumed that the material is modeled by the Baumgaertel, Schausberger and Winter (BSW) spectrum [40,41], which was used to describe the viscoelasticity of various polymers; for example, polydisperse polymer melts [42,43], polymethylmethacrylate (PMMA) [56], polybutadiene (PBD) [56] and polymer pelts [57].
The “real” material and the FMM model were simulated in Matlab R2023b, The Mathworks, Inc., Natick, MA, USA. Functions MLFFIT2 [58] and MLF [59], provided by Podlubny, were used for the FMM simulation and numerical solution of the optimal identification tasks.

3.5. Material I

Consider the material whose relaxation spectrum is described by the unimodal Gauss-like distribution:
H τ = ϑ e 1 τ m 2 / q / τ ,
where the parameters are as follows [60]: ϑ = 31520   P a · s , m = 0.0912 s 1 and q = 3.25 × 10 3 s 2 . The related relaxation modulus is [60]
G t = π q 2 ϑ   e 1 4 t 2 q m t e r f c 1 2 t q m q ,
where the complementary error function e r f c x is given by [4] (Equation (C.2))
e r f c x = 2 π     x e z 2 d z .
Following [47], for numerical simulations, the time interval T = 0 , 200 seconds is chosen. Hence, the weighting function in Q g (26) is ρ t = 1 200 s 1 . The elements of the optimal parameter vector g solving the measurement-independent optimization task (27) are given in Table 1.
The N sampling instants t i for the simulated stress relaxation test were selected randomly according to the uniform distribution on T . A normal distribution with zero mean value and variance σ 2 was applied to the random independent generation of the additive measurement noises z i . In the noise robustness analysis, the standard deviations σ = 2 , 5 , 8   P a were used. In the analysis of the model asymptotic properties, for any σ numbers of measurements, N N have been applied, where N = 50 ; 100 ; 200 ; 500 ; 1000 ; 2000 ; 5000 ; 7000 ; 10,000 ; 12,000 ; 15,000 .

3.5.1. Asymptotic Properties

Then, for every pair N , σ , the optimal parameter g N was determined through solving the minimization task (25). The elements of the vectors g N , the mean square indices Q N g N and integral Q g N indices, and the relative percentage errors of the approximation of the measurement-independent parameter g , defined as
E R R = g N g 2 2 / g 2 2 · 100 % ,
are given in Table A1, Table A2 and Table A3 for the three standard deviations of the noises. The model approximation error was also estimated via the relative mean error defined as (compare (24))
Q N r e l g = 1 N i = 1 N G ¯ T i G M T i , g 2 G ¯ T i 2 .
The optimal model parameters g N as the functions of the number of measurements N are illustrated by Figure 2 for the noises of σ = 2 , 5 , 8   P a . In any subplot, the values of the related parameters of the sampling points-independent model g are depicted by horizontal purple lines. The asymptotic properties are also illustrated by Figure 3 juxtaposing the empirical index mean-square index Q N g N , Equation (24), and the integral quadratic sampling instants-independent index Q g N , Equation (26), as the functions of N with the index Q g , marked with horizontal lines. In Figure 2 and Figure 3, a logarithmic scale is applied for the horizontal axes. These plots confirm the asymptotic properties of the proposed identification algorithm. The convergence of g N to the parameter g is directly translated into the convergence of Q g N into Q g , especially for N 5000 . The values of the index Q N g N for N = 50 , small compared to those for N 100 (see Table A1, Table A2 and Table A3), result from the good fit of the FMM, whose four parameters are optimally selected in problem (25), to only 50 measurement points. For more measurement points, such a good fit is, generally, impossible whenever the real characteristic does not depend on the pre-assumed class of models. A comparison of Figure 2 and Figure 3b with Figure 3a shows that the impact of stronger noises on the values of the empirical index Q N g N is much stronger than the impact of the noises on the values of the FMM parameter g N and, consequently, also on the integral index Q g N , which does not directly depend on the measurement noises. Given Equation (28), this property is natural and fully justified.
The quality of the real modulus G t approximation by the FMM is illustrated in Figure 4, where the measurements G ¯ i of the real modulus G t fitted by the optimal model G M t , g N are plotted for the N = 100 and N = 10,000 measurements and the strongest disturbances; σ = 8   P a . Although, for the N = 100 measurements, the models G M t , g N and G M t , g differ slightly (see small subplot), for the N = 10,000 measurements, they are practically identical, which is confirmed by the values of E R R (37) equaling 0.52% for N = 100 and equaling only 5.58 × 10−4 % for N = 10,000 (see Table A3). Even for the strongest noises, the relative errors E R R (37) of the parameters g and g N discrepancy is smaller than 0.002% for N 200 . This almost excellent fitting of the experiment data by the model G M t , g N is confirmed by the values of the relative square model approximation index Q N r e l g N (38), which for N 200 and the weakest noises does not exceed 0.015%, while for the strongest noises, it does not exceed 0.28%. For the noises considered, the values of the model fit indices Q N g N (24) and Q N r e l g N (38) and the integral quadratic index Q g N (26) indicate an excellent fit of the model to the experiment data and the fast convergence of g N to g as N tends to infinity; compare Table A1, Table A2 and Table A3.

3.5.2. Noise Robustness

To examine the effect of the measurement noises, for every pair N , σ , the simulated experiment was repeated n = 50 times. In each experiment repetition, the measurement noises z i were generated independently and randomly with a normal distribution, with a zero mean value and variance σ 2 .
Having in mind the definition of the index Q N r e l g (38), for the n -element sample, the mean relative relaxation modulus approximation error was determined as follows:
E R R Q N r e l = 1 n j = 1 n Q N r e l g N , j ,
for any pair N , σ , where the vector of the optimal FMM parameters g N , j was computed for j-th experiment repetition, j = 1 , , n .
For the true relaxation modulus approximation, the mean optimal integral error
E R R Q = 1 n j = 1 n Q g N , j
was determined for every pair N , σ .
The distance between the vector g N , j and the measurement-independent vector g for the n element sample was estimated by the mean relative error defined as follows (compare E R R (37)):
M E R R = 1 n j = 1 n g N , j g 2 2 / g 2 2 · 100 % .
The indices E R R Q N r e l a (39) and E R R Q (40), as the functions of N and σ , are depicted in the bar in Figure 5, while the index M E R R (41) is shown in Figure 6.
From Figure 5b, it is seen that for N > 2000 , the number of measurements do not essentially affect the integral index E R R Q , either for weak or strong noises, while both the empirical index E R R Q N r e l and mean relative error M E R R decrease exponentially with the increasing number of measurements, which confirms the analytical analysis performed above. The M E R R index is of order 0.55% for N = 100 , it does not exceed 10 3   % for N 1000 and is smaller than 5 × 10 5   % even for the strongest disturbances. This, practically, means determining the sampling points-independent parameter g . The algorithm ensures the very good quality of the measurement approximation even for large noises. The values of the relative relaxation modulus approximation error E R R Q N r e l , which due to the “real” modulus model difference is lower bounded by 3.191 × 10 4 %, already for N 100 measurements do not exceed 0.35 % , and for N 1000 , fall below 0.028%. The course of the mean integral sampling points-independent index E R R Q (40) as the function of N indicates the asymptotic independence of the model from the sampling points, especially for N 5000 .

3.6. Material II

Consider the material described by the empirical spectrum of relaxation times τ introduced by Baumgaertel, Schausberger and Winter [40,41],
H τ = β 1 τ τ c ρ 1 + β 2 τ τ c ρ 2 e τ τ m a x ,
which is known to effectively describe polydisperse polymer melts [42,43], with the coefficients [43,47,61] as follows: β 1 = 6.276 × 10 2   M P a , β 2 = 0.127   M P a , τ c = 2.481   s , τ m a x = 2.564 × 10 4   s , ρ 1 = 0.25 and ρ 2 = 0.5 . The spectrum H τ uniquely defines the relaxation modulus G t by the following integral [44]:
G t = 0 H τ τ e t / τ d τ .
Following [47], the time interval T = 0 , 2000 seconds is taken for the experiment simulations; the weighting function is ρ t = 1 2000 s 1 . The elements of the optimal parameter vector g , which solve the measurement-independent optimization task (27), are given in Table 2.
As previously described, in the simulations, the sampling points t i were randomly selected according to the uniform distribution on T . The standard deviations σ = 3 , 6 , 8   k P a of the random normally distributed noises z i combined with the number of measurements N N were used for the analysis of the model asymptotic properties.

3.6.1. Asymptotic Properties

For every pair N , σ , the elements of the optimal model parameter g N , the empirical Q N g N , Q N r e l g N and integral Q g N indices and the relative percentage errors E R R (37) are given in Table A4, Table A5 and Table A6 in Appendix B. The dependence of the optimal model parameters g N on the number of measurements N for the noises of σ = 3 , 6 , 8   k P a are illustrated by Figure 7. Figure 8 illustrates the empirical Q N g N and integral Q g N indices as the functions of N ; the value of Q g is marked with purple horizontal lines. These plots confirm the asymptotic properties of the proposed identification algorithm. Figure 8a shows the impact of noises on the values of the empirical index Q N g N .
The approximation of the real modulus G t by the FMM is illustrated in Figure 9, where the measurements G ¯ i of the real modulus G t along with optimal models G M t , g N and G M t , g are plotted for the N = 100 and N = 10,000 measurements and the strongest noises σ = 8   P a . However, for N = 100 , the model parameter error E R R = 5.35 % , while for N = 10,000 , we have E R R = 0.15 % ; both for N = 100 and N = 10,000 , the models G M t , g N and G M t , g differ slightly and the respective empirical indices are Q N r e l g N = 3.33 × 10−4 % and Q N r e l g N = 2.0 × 10−7%, respectively.

3.6.2. Noise Robustness

For every pair N , σ , the simulated experiment was repeated n = 50 times. The mean relative relaxation modulus approximation error E R R Q N r e l (39), the mean optimal integral error E R R Q (40) and the mean relative error of the parameter g approximation M E R R (41) were determined. The indices E R R Q N r e l a and E R R Q are depicted in Figure 10 as the functions of N and σ . Figure 11 illustrates the dependence of the index M E R R on N and σ .
The mean integral error E R R Q for N 12,000 does not depend essentially on the number of measurements, either for small or large noises (see Figure 10b), while both the empirical index E R R Q N r e l and mean relative error M E R R decrease exponentially with the increasing number of measurements, the M E R R for N 7000 . For N 1000 , the M E R R index does not exceed 1.01%, for N 7000 it does not exceed 0.22%, while for N 10,000 , it falls below 0.05% even for the strongest disturbances. The globally optimal parameter g was determined. As is seen from Figure 9, the algorithm practically ensures an excellent quality of the relaxation modulus approximation even for the strongest noises. The values of the relative relaxation modulus approximation error E R R Q N r e l , already for the N 100 measurements, do not exceed 3.3 × 10−4% and for N 1000 , fall below 8.3 × 10−6%. From the course of the mean integral sampling points-independent index E R R Q (40), as the function of N , we can conclude that the model is practically independent on the sampling instants for N 12,000 , independently on the measurement noises. The above combined with the close to zero values of E R R Q N r e l means the determining of the globally optimal model with the parameter g . In conclusion, the courses of both the index E R R Q N r e l (38), and the indices M E R R (41) and E R R Q (40) as the functions of N , indicate the asymptotic independence of the model from the sampling points for a sufficiently large number of measurements.

4. Conclusions

The fractional Maxwell model allows for the modeling of a very wide range of stress relaxation processes in materials. The goal of the FMM identification is, generally, not to achieve a true description of the genuine relaxation modulus, but one that is “optimally accurate” in the assumed sense of the square weighted approximation error and does not depend on the particular sampling instants used in the stress relaxation experiment. The stochastic-type analytical analysis and numerical studies demonstrated that, despite the fact that the real description of the relaxation modulus is completely unknown, an arbitrarily exact approximation of the sampling points-independent optimal FMM can be identified based on the relaxation modulus data sampled randomly, according to respective randomization, when the number of the measurements applied in the experiment appropriately grow large. The four parameters of the approximate FMM are strongly consistent estimates of the parameters of the sampling points-independent model minimizing the integral square approximation error. The resulting identification scheme is simple and useful in application. It requires only the a priori, before the experiment is performed, independent random choice of the time instants at which the relaxation modulus is recorded from the assumed set according to a stationary rule.
Although this article is about modeling the relaxation modulus, the proposed identification scheme can also be successfully applied to the identification of the fractional-order models of creep compliance using the measurements obtained in the retardation test, whenever the respective set of sampling instants is open to manipulation during experimental data collecting. Therefore, the applicability of the identification asymptotically independent of the time instants used in the rheological experiment, to other fractional-order models determination, in particular, Kelvin–Voight, Zener and anti-Zener models, can be the subject of future research.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The author declares no conflicts of interest.

Appendix A

Appendix A.1. Proof of Theorem 1

Note, firstly, that the assumption g G 1 means, in particular, that 0 < β 0   β < α 1 . The assumption t T , by definition of the set T = t 0 , T , where t 0 > 0 , implies inequality t > 0 .
Differentiability and whence continuity of the relaxation modulus model G M t , g (7) with respect to parameter G e is obvious. The two-parameter Mitteg-Leffler function E κ , μ x (6) is known to be differentiable and continuous with respect to real argument x and parameters κ and μ [36]. Therefore, by (13) and (14), the differentiability of the function G M t , g with respect to the positive relaxation time parameter τ r and parameters α and β directly results. To show that the condition (22) is satisfied, it is enough to prove that the partial derivatives of G M t , g with respect to the four model parameters are bounded uniformly on T × G 1 .
To prove the boundness of the partial derivatives with respect to the relaxation time parameter τ r and the orders α and β of the stress and strain derivatives, the spectral representation (21) of the FMM will be applied together with the property concerning the absolute boundness, uniform on the set T × G 1 , of some definite integrals, which results from the following known [62] property concerning the absolute integrability of the product of absolutely integrable and bounded functions.
Property A1
([62]). If the function  f x  is absolutely integrable in the interval  a ,  and the function  g x  is bounded in  a , , then the product  f x g x  is absolutely integrable in  a , .
Let us consider definite integral
I 0 t , g = 0 r 0 v , t , g f 0 v , t , g d v ,
where the function f 0 v , t , g R + × T × G 0 R is absolutely integrable with respect to the variable of integration v in R + , uniformly on the set T × G 0 , where R + = 0 , , G 0 R k , k 1 and the function r 0 v , t , g : R + × T × G 0 R is absolutely bounded, uniformly on R + × T × G 0 . The first assumption means that there exists a positive constant m ¯ such that
0 f 0 v , t , g d v m ¯ < ,
for any t T and any g G 0 , while the second assumption yields
r 0 v , t , g m ̿ < ,
for any v , t , g R + × T × G 0 . From Property A1, the convergence of the integral
0 r 0 v , t , g f 0 v , t , g d v ,
and, in consequence, of the integral I 0 t , g (A1) for any t T and any g G 0 follows. In view of (A3) and (A4), we have
I 0 t , g 0 r 0 v , t , g f 0 v , t , g d v m ̿ 0 f 0 v , t , g d v m ̿ m ¯ = M ¯ .
Therefore, the next result holds.
Property A2.
If the function  f 0 v , t , g  is absolutely integrable with respect to  v  in the interval  R + = 0 , , uniformly on the set  T × G 0  of the rest arguments  t , g , and the function  r 0 v , t , g  is bounded uniformly on  R + × T × G 0 , then the product  r 0 v , t , g f 0 v , t , g  is absolutely integrable in  R +  for any  t , g T × G 0  and the integral  I 0 t , g  (A1) is absolutely bounded uniformly on  T × G 0 .
Below, the proof is divided into four parts related to the four model parameters.

Appendix A.1.1. Uniform on T × G 1 Boundness of the FMM Derivative with Respect to G e

From (7), we have
G M t , g G e = t τ r β E α β , 1 β t τ r α β ,
whence, by (12) and (16), the uniform boundness of the above derivative on the set T × G 1 follows.

Appendix A.1.2. Uniform on T × G 1 Boundness of the FMM Derivative with Respect to τ r

By (7) and (13), we can express G M t , g as follows
G M t , g = G e e α β , 1 β t τ r ; 1 .
Whence, the partial derivative with respect to the relaxation time, in the respective notation, is given by
G M t , g τ r = G e 1 t τ r 2 d d x e α β , 1 β x ; 1 x = t τ r ,
where, due to the complete monotonicity of the function e α β , 1 β x ; 1 , the negative derivative d d x e α β , 1 β x ; 1 monotonically increases to zero for x > 0 .
From (A5) and (A6), the nonnegative definiteness of G M t , g G e and G M t , g τ r , being positive for any g G 1 and t < , follows; therefore, Property 1 is formulated.
To examine the asymptotic properties of G M t , g τ r as t and as t t 0 , let us express (A6), applying the known differentiation formula [4] (Equation (E.55))
e κ , μ x ; λ = d d x e κ , μ + 1 x ; λ ,
in the form
G M t , g τ r = G e 1 t τ r 2 e α β , β t τ r ; 1 ,
or, having in mind definition (13), directly in terms of the Mittag-Leffler function
G M t , g τ r = 1 G e τ r t τ r β E α β , β t τ r α β .
The following asymptotic approximation [63] (Equation (12)), see also [4] (Equation (E.30)):
E κ , μ t κ n = 1 1 n 1 t κ n Γ μ κ n ,
which holds for t , applied to (A7), yields
G M t , g τ r 1 G e τ r t τ r β t τ r α β Γ α + n = 2 1 n 1 t τ r α β n Γ β α β n ,
whence, for large times, especially for t τ r α β 1 , we obtain the asymptotic long-time approximation
G M t , g τ r 1 G e τ r Γ α t τ r α ,
where means “approximately equal”. Therefore, derivative G M t , g τ r tends to zero as t for any admissible parameter g G 1 .
To estimate the value of G M t , g τ r (A7) for t = t 0 , the series representation
G M t , g τ r = G e τ r t τ r β n = 0 1 n + 1 t τ r α β n Γ α β n β
resulting directly from (6) and (A7), is used. The first summand of the series is positive, while the next elements are positive or negative, depending on the index n and the relation between parameters α and β . Since
t τ r β n = 0 1 n + 1 t τ r α β n Γ α β n β t τ r β n = 0 t τ r α β n Γ α β n β ,
and the argument of the gamma function is such that
α β n β β > 1 ,
which in view of the monotonicity of the gamma function implies
Γ α β n β Γ x m i n = 0.8856032 > 0.8856 ,
where x m i n 1.4616321 is the real nonnegative argument at which a minimum of the function Γ x occurs [64]. In view of (A8) and (A9), having in mind the nonnegative definiteness of G M t , g τ r , we obtain the following estimation
G M t , g τ r < G e 1.1292 τ r t τ r β n = 0 t τ r α β n ,
which, by the inequalities τ r m i n τ r τ r m a x , for t = t 0 , implies the next estimation
G M t , g τ r t = t 0 < G e 1.1292 τ r m i n τ r m a x τ r m i n β t 0 τ r m i n β n = 0 t 0 τ r m i n α β n .
By the assumption t 0 < τ r m i n , the above estimation can be rewritten in compact form as
G M t , g τ r t = t 0 < G e 1.1292 τ r m i n τ r m a x τ r m i n β t 0 τ r m i n β 1 t 0 τ r m i n α β .
Since, for an arbitrary β 0   β < α 1 , we have
t 0 τ r m i n β 1 t 0 τ r m i n α β = 1 t 0 τ r m i n β t 0 τ r m i n α < 1 t 0 τ r m i n β ,
inequality (A10) for any g G 1 implies
G M t , g τ r t = t 0 < G e 1.1292 τ r m i n τ r m a x t 0 β ,
which means that
G M t , g τ r t = t 0 < G e m a x   1.1292 τ r m i n τ r m a x t 0 = G e m a x   1.1292 τ r m i n m 0 ,
i.e., the derivative for t = t 0 is bounded, uniformly on the set G 1 , where positive parameter m 0 is defined in Equation (16).
Since the continuous function G M t , g τ r of the time t is bounded for t = t 0 and for t for any fixed g G 1 , derivative G M t , g τ r as a function of the time is bounded both for the bounded and not-bounded set T . However, due to the non-compactness of the set G 1 , from which α = β is excluded, the uniform on T × G 1 boundness of G M t , g τ r is not obvious. Therefore, it should be examined if the maximum of G M t , g τ r (with respect to the time) is bounded, as α β + . To this end, an alternative to (A6) and (A7), the representation of G M t , g τ r is derived based on the spectral representation given by Equation (21).
Differentiation of (21) on both sides with respect to τ r yields
G M t , g τ r = G e α τ r α 1 π 0 r v , g v α 1 e t v d v + G e τ r α π 0 α β τ r v α β 1 sin π β τ r v 2 α β + 2 τ r v α β cos π α β + 1 v α e t v d v G e τ r α π 0 2 α β τ r v 2 α β 1 + 2 α β τ r v α β 1 cos π α β τ r v 2 α β + 2 τ r v α β cos π α β + 1 r v , g v α e t v d v ,
where the function
r v , g = τ r v α β sin π β + sin π α τ r v 2 α β + 2 τ r v α β cos π α β + 1 ,
by (21), is such that
G M t , g = G e τ r α π 0 r v , g v α 1 e t v d v .
Whence, introducing the notations
r 1 v , g = τ r v α β sin π β τ r v 2 α β + 2 τ r v α β cos π α β + 1 ,
r 2 v , g = τ r v 2 α β + τ r v α β cos π α β τ r v 2 α β + 2 τ r v α β cos π α β + 1 ,
Equation (A11) can be rewritten as
G M t , g τ r = α τ r G M t , g + α β G e τ r α 1 π 0 r 1 v , g v α 1 e t v d v 2 α β G e τ r α 1 π 0 r 2 v , g r v , g v α 1 e t v d v ,
or in a more compact form as a linear combination:
G M t , g τ r = α τ r G M t , g + α β G e τ r α 1 π I 1 t , g 2 α β G e τ r α 1 π I 2 t , g ,
where the integrals:
I 1 t , g = 0 r 1 v , g v α 1 e t v d v ,
I 2 t , g = 0   r 2 v , g r v , g v α 1 e t v d v
The denominator
q v , g = τ r v 2 α β + 2 τ r v α β cos π α β + 1
of the fractions r v , g , r 1 v , g and r 2 v , g is positive for any v 0 , whenever α β 1 , i.e., for any admissible parameter g G 1 , which satisfies the following inequalities
0 < α β 1 β 1 β 0 < 1 .
For α β + , the denominator q v , g τ r v α β + 1 2 1 . By (A13), (17) and nonnegative definiteness of r v , g (A12) on the set R + × G 1 , we have
0 r v , g v α 1 e t v d v = 0 r v , g v α 1 e t v d v = π G e τ r α G M t , g π G e τ r α M 1
for any t , g T × G 1 , i.e., the function r v , g v α 1 e t v as the function of the variable v is absolutely integrable, uniformly on the set T × G 1 , with the constant m ¯ (compare (A2)) given by
0 r v , g v α 1 e t v d v m ¯ = π G e m i n γ 2 M 1 ,
where the parameter
0 < γ 2 = min β 0 α 1 τ r m i n α < .
In view of Property A2, bearing in mind inequality (17), or (15), to prove the absolute uniform boundness of the derivative G M t , g τ r (A16) on the set T × G 1 , it is enough to demonstrate that the integrals I 1 t , g (A17) and I 2 t , g (A18) are convergent and absolutely bounded, uniformly on the set T × G 1 . For this purpose, we express the two integrals as definite integrals of the product of some absolutely integrable function and a bounded function.
The continuous, nonnegative definite for any v , g R + × G 1 , function r 1 v , g (A14) is equal to zero for v = 0 , tends to zero for v and takes the maximal value for v = 1 / τ r , independently on the values of α and β . Whence, for any v , g R + × G 1 , the inequalities hold
r 1 v , g sin π β 2 cos π α β + 2 1 2 cos π 1 β 0 + 2 = 1 2 2 cos π β 0 = m 1 < ,
i.e., r 1 v , g is absolutely bounded uniformly on R + × G 1 . By the following notable integral [4] (Equation (A.21):
0 v α 1 e t v d v = 0 v α 1 e t v d v = Γ α t α ,
which holds for any α > 0 and t > 0 , v α 1 e t v is the absolutely integrable function of the variable v 0 for any t T . Recalling the definitions of the sets T and G 1 (19), and the monotonicity of the gamma function Γ α for 0 < α 1 , we immediately obtain the estimation
0 v α 1 e t v d v Γ β 0 γ 1 ,
valid for any t T and any g G 1 , where
0 < γ 1 = min β 0 α 1 t 0 α < .
Therefore, according to Property A2, the integral I 1 t , g (A17) is convergent for any t > 0 and any g G 1 and absolutely bounded by the upper bound equal to m 1 Γ β 0 / γ 1 , uniformly on the set T × G 1 .
It is easy to check that continuous function r 2 v , g (A15) is equal to zero for v = 0 and tends to 1, as v , independently on the values of α and β from the set G 1 , i.e., for α > β . Function r 2 v , g can be expressed as
r 2 v , g = 1 τ r v α β cos π α β q v , g 1 q v , g ,
with q v , g defined by (A19), where the absolute value of the second summand takes the maximal value for v = 1 / τ r , while the third summand takes the maximal value whenever τ r v α β = cos π α β , if cos π α β < 0 , and for v = 0 in the opposite case. Therefore, for any v , g R + × G 1 , the next estimation holds
r 2 v , g 1 + cos π α β 2 cos π α β + 2 + m α , β ,
where
m α , β = 1 1 cos 2 π α β             i f   cos π α β < 0   1                                                       i f   cos π α β 0 .
Since, for g G 1 , the inequality α > β holds, for α β + , we have cos π α β 1 and m α , β 1 . Simultaneously, if α 1 and β β 0 , then m α , β 1 1 cos 2 π β 0 > 1 , whenever β 0 < 1 2 and m α , β 1 for β 0 1 2 . Therefore, for any v , g R + × G 1 , by (A26) and (A27), we have
r 2 v , g 1 + 1 2 2 cos π β 0 + 1 1 cos 2 π β 0 = 1 + m 1 + m 2 < ,
where
m 2 = 1 1 cos 2 π β 0 < ,
and according to Property A2, the integral I 2 t , g (A18) is convergent for any t T and any g G 1 and absolutely bounded uniformly on T × G 1 , with the upper bound π 1 + m 1 + m 2 G e m i n γ 2 M 1 resulting from (A20) and (A28).
Combining the absolute boundness of the three summands of the right-hand side of (A16), uniform on the set T × G 1 , the respective uniform boundness of G M t , g τ r is proved.

Appendix A.1.3. Uniform on T × G 1 Boundness of the FMM Derivative with Respect to β

Differentiation of Equation (21) on both sides with respect to β yields
G M t , g β = G e τ r α π 0 l n τ r v τ r v α β sin π β + π τ r v α β cos π β τ r v 2 α β + 2 τ r v α β cos π α β + 1 v α 1 e t v d v + 2 G e τ r α π 0 r v , g l n τ r v τ r v 2 α β + l n τ r v τ r v α β cos π α β π τ r v α β s i n π α β τ r v 2 α β + 2 τ r v α β cos π α β + 1 v α 1 e t v d v
where r v , g is given by (A12), which using the notations r 1 v , g (A14) and r 2 v , g (A15) and introducing functions
r 3 v , g = τ r v α β cos π β τ r v 2 α β + 2 τ r v α β cos π α β + 1 ,
r 4 v , g = τ r v α β s i n π α β τ r v 2 α β + 2 τ r v α β cos π α β + 1 ,
can be rewritten in compact form as a linear combination:
G M t , g β = G e τ r α I 3 t , g 2 G e τ r α I 4 t , g G e τ r α π I 5 t , g + 2 G e τ r α π I 6 t , g ,
of four integrals:
I 3 t , g = 0 r 3 v , g v α 1 e t v d v ,
I 4 t , g = 0   r 4 v , g r v , g v α 1 e t v d v ,
I 5 t , g = 0 r 1 v , g l n τ r v v α 1 e t v d v ,
I 6 t , g = 0 r 2 v , g r v , g l n τ r v v α 1 e t v d v .
To prove the absolute boundness of the derivative G M t , g β (A32), uniform on the set T × G 1 , it is enough to demonstrate that the four above integrals are convergent and absolutely bounded, uniformly on T × G 1 .
For any v , g R + × G 1 , the continuous function r 3 v , g (A30) satisfies the following inequalities (compare (A14) and (A22))
r 3 v , g cos π β 2 cos π α β + 2 1 2 cos π 1 β 0 + 2 = 1 2 2 cos π β 0 = m 1 ,
where m 1 is defined in (A22), i.e., r 3 v , g is absolutely bounded, uniformly on R + × G 1 , which combined with the absolutely integrability of the function v α 1 e t v for any t T and any g G 1 , according to Property A2, implies the convergence of the integral I 3 t , g (A33) for any t T and any g G 1 and its absolute boundness by the upper bound m 1 Γ β 0 / γ 1 resulting from and (A24), (A25) and (A37), valid uniformly on the set T × G 1 .
The nonnegative function r 4 v , g (A31) is for any v , g R + × G 1 bounded by
r 4 v , g s i n π α β 2 cos π α β + 2 1 2 cos π 1 β 0 + 2 = 1 2 2 cos π β 0 = m 1 ,
similarly to r 1 v , g (A14) and r 3 v , g (A30), which combined with the absolutely integrability of r v , g v α 1 e t v for any g G 1 and any t T , implies the convergence of the integral I 4 t , g (A34) for any t , g T × G 1 and the absolute boundness of I 4 t , g by the upper bound π   m 1 G e m i n γ 2 M 1 , derived from (A20) and (A38), uniformly on the set T × G 1 .
To demonstrate the convergence and uniform boundness of the integrals I 5 t , g (A35) and I 6 t , g (A36), let us express them as follows
I 5 t , g = 0 r 12 v , g r 5 v , t , g v α 1 e 1 2 t v d v ,
I 6 t , g = 0 r 22 v , g r 5 v , t , g r v , g v α 1 e 1 2 t v d v .
where, by (A14) and (A15)
r 12 v , g = sin π β τ r v 2 α β + 2 τ r v α β cos π α β + 1 ,
r 22 v , g = τ r v α β + cos π α β τ r v 2 α β + 2 τ r v α β cos π α β + 1 ,
and nonnegative continuous function
r 5 v , t , g = τ r v α β l n τ r v e 1 2 t v = l n τ r v τ r v α β e 1 2 t v .
Since (compare (A13)),
G M t 2 , g = G e τ r α π 0 r v , g v α 1 e 1 2 t v d v ,
the nonnegative function r v , g v α 1 e 1 2 t v is absolutely integrable for any t , g T × G 1 with the upper bound π G e m i n γ 2 M 1 ; compare (A20). In turn, by (A23), v α 1 e 1 2 t v is an absolutely integrable function of the variable v 0 , uniformly on the set T × G 1 , with the upper bound 2 Γ β 0 / γ 1 ; compare (A24) and (A25).
Functions r 12 v , g (A41) and r 22 v , g (A42) are absolutely bounded uniformly on R + × G 1 , since the following estimations hold for any v , g R + × G 1 :
r 12 v , g 1 1 cos 2 π β 0 = m 2 ,
r 22 v , g 1 2 2 cos π β 0 + 1 1 cos 2 π β 0 = m 1 + m 2 < ,
where constants m 1 and m 2 are defined in (A22) and (A29), respectively.
To examine the properties of r 5 v , t , g (A43), the asymptotic properties as v 0 + and v are studied. This function is expressed as
r 5 v , t , g = l n τ r v τ r v α β e 1 2 t v ,
where the nominator tends to and the denominator tends to + , as the variable v 0 + . Therefore, by applying the L’Hospital’s rule, in view of α > β > 0 , we obtain
lim v 0 + r 5 v , t , g = l i m v 0 + 1 τ r α β + 1 2 t   τ r v τ r v α β e 1 2 t v = 0 .
Since τ r v α β l n τ r v tends to + , while e 1 2 t v tends to zero, as the variable v tends to infinity, using the L’Hospital’s rule double times to the right expression in (A43), we have
lim v r 5 v , t , g = l i m v τ r + α β l n τ r v τ r 1 2 t e 1 2 t v τ r v 1 α β ,
and next
lim v r 5 v , t , g = l i m v τ r + α β l n τ r v τ r 1 2 t e 1 2 t v τ r v 1 α β = l i m v α β 1 2 t 1 2 t v + 1 α + β τ r v 1 α + β e 1 2 t v = 0 + .
Applying two known inequalities e x < k ! k ! + x k and [65]
ln x n x 1 n 1 ,
being valid for any integer n , k > 0 and real x > 0 , by putting n = 1 and k = 2 function r 5 v , t , g (A43) for any v , t , g R + × T × G 1 can be estimated by
r 5 v , t , g 2 τ r v α β τ r v 1 2 + 1 4 t τ r 2 τ r v 2 τ r v α β τ r v + 1 1 + 1 8 t 0 τ r 2 τ r v 2 τ r m a x v α β τ r m a x v + 1 1 + 1 8 t 0 τ r 2 τ r v 2 .
For τ r m a x v < 1 , the right inequality in (A47) implies
r 5 v , t , g < 2 1 + 1 8 t 0 τ r 2 τ r v 2 2 ,
while for τ r m a x v 1 , the middle inequality in (A47) yields
r 5 v , t , g τ r v α β 1 + τ r v α β 2 1 8 t 0 τ r 2 2 1 8 t 0 τ r 2 16 t 0 τ r m a x 2 = 16 m 0 2 .
where positive m 0 is defined in Equation (16). Combining (A48) and (A49), we obtain the inequality
r 5 v , t , g < m a x 16 m 0 2 , 2 = m 4 ,
valid for any v , t , g R + × T × G 1 , which together with (A44) and (A45) means an absolute boundness of continuous functions r 12 v , g r 5 v , t , g and r 22 v , g r 5 v , t , g , respectively, and in view of the absolute integrability of v α 1 e 1 2 t v and r v , g v α 1 e 1 2 t v , imply the convergence of the integrals I 5 t , g (A39) and I 6 t , g (A40) for any t T and any g G 1 . The absolute boundness of I 5 t , g and I 6 t , g , uniform on the set T × G 1 , with upper bounds estimations 2 Γ β 0 m 2 m 4 γ 1 and π M 1 m 1 + m 2 m 4 G e m i n γ 2 , respectively, follows from Property A2. Therefore, the absolute boundness of G M t , g β (A32), uniform on T × G 1 , is proved.

Appendix A.1.4. Uniform on T × G 1 Boundness of the FMM Derivative with Respect to α

The same integral properties and the spectral representation (21) will be applied to prove the boundness of the partial derivative with respect to parameter α . Differentiation (21) on both sides with respect to α yields
d G M t , g d α = ln τ r G M t , g + G e τ r α π 0 ln τ r v τ r v α β sin π β + π cos π α τ r v 2 α β + 2 τ r v α β cos π α β + 1 v α 1 e t v d v 2 G e τ r α π 0 r v , g ln τ r v τ r v 2 α β + ln τ r v τ r v α β cos π α β π τ r v α β sin π α β τ r v 2 α β + 2 τ r v α β cos π α β + 1 v α 1 e t v d v + G e τ r α π 0 r v , g v α 1 ln v e t v d v ,
where r v , g is described by (A12). Having in mind (A13), recalling the notations r 1 v , g (A14), q v , g (A19), r 2 v , g (A15), r 4 v , g (A31), r 5 v , t , g (A43) and introducing the integrals
I 7 t , g = 0 cos π α q v , g v α 1 e t v d v ,
I 8 t , g = 0   r v , g r 5 v , t , g v β 1 e 1 2 t v d v ,
we can express d G M t , g d α (A51) in a compact form as
d G M t , g d α = G e τ r α π I 5 t , g + G e τ r α I 7 t , g 2 G e τ r α π I 6 t , g + 2 G e τ r α I 4 t , g + G e τ r α π I 8 t , g ,
where the integrals I 5 t , g (A35), I 6 t , g (A36) and I 4 t , g (A34) are absolutely bounded uniformly on T × G 1 . Therefore, only the convergence and boundness of the two new integrals I 7 t , g and I 8 t , g must be proved.
The boundness of cos π α q v , g , uniform on R + × G 1 , with the upper bound
cos π α q v , g 1 1 cos 2 π β 0 = m 2 ,
combined with the absolute integrability of v α 1 e t v for any t , g T × G 1 , yields the convergence of the integral I 7 t , g (A52) and its absolute boundness, uniform on T × G 1 .
From (A12), the upper bound of the nonnegative function r v , g follows
r v , g 1 2 2 cos π β 0 + 1 1 cos 2 π β 0 = m 1 + m 2 ,
where constants m 1 and m 2 are defined in (A22) and (A29), respectively; therefore, the absolute boundness of r 5 v , t , g , uniform on R + × T × G 1 (c.f., (A50)), and integrability of v β 1 e 1 2 t v imply both the convergence and the absolute boundness of I 8 t , g , uniformly on the set T × G 1 .
The partial derivatives of the FMM with respect to the four model parameters are proved to be absolutely bounded uniformly on the set T × G 1 ; therefore, the uniform boundness (22) of the gradient g G M t , g follows. The theorem is proved.

Appendix A.2. Proof of Theorem 2

Since, for any t T function, G M t , g (10) is differentiable with respect to g and in this case β = α , the four-element vector of model parameters g (8) is as follows
g = α α G e τ r T ,
the gradient is given by
g G M t , g = G M t , g α G M t , g α 1 2 Γ 1 α t τ r α G e α 2 Γ 1 α t t τ r 1 α T ,
with the partial derivative
G M t , g α = G e 2 Γ 1 α ψ 1 α + l n τ r t τ r t α ,
where
ψ x = Γ x Γ x
is the digamma (or psi) function [66] defined as the logarithmic derivative of the gamma function [66]:
ψ x = d d x l n Γ x .
For 0 < β 0 α < 1 , the nonnegative digamma function ψ 1 α strictly decreases from finite ψ 1 β 0 < 0 to , while the positive gamma function Γ 1 α strictly increases from Γ 1 β 0 > 0 to + . To evaluate the first summand of G M t , g α (A54), the following result proved by Mező and Hoffman [67] is helpful in providing an infinite product representation of the entire function ψ z Γ z .
Property A3
([67] (Theorem 2.1)). For all  z C
ψ z Γ z = e 2 γ z k = 0 1 z α k e z / α k ,
where  C  is the set of complex numbers,  γ  is the Euler’s constant and  α k  are the zeros of the digamma function  ψ z .
It is known [67] that the zeros α k are real, and all but one are negative; here, α 0 is the positive zero, and α 1 ,   α 2 , are the negative ones in decreasing order. The course of ψ z Γ z for the real 0 z 1 is illustrated by Figure A1. For α = 1 , by (A56), the quotient ψ 1 α Γ 1 α is as follows ψ 0 Γ 0 = 1 , while for α = β 0 , the quotient ψ 1 α Γ 1 α tends to ψ 1 Γ 1 = γ . Therefore, for any t , g T × G 2
E ψ 1 α 2 Γ 1 α τ r t α G e m a x 2 τ r m a x t 0 ,
i.e., the first summand of G M t , g α (A54) is bounded.
Figure A1. The quotient ψ z Γ z of the psi function ψ z (A55) by the gamma function Γ z for real argument 0 z 1 .
Figure A1. The quotient ψ z Γ z of the psi function ψ z (A55) by the gamma function Γ z for real argument 0 z 1 .
Materials 17 01527 g0a1
Using the inequality (A46), the second summand of G M t , g α (A54) can be estimated as follows
G e 2 Γ 1 α l n τ r t τ r t α m G e 2 Γ 1 α τ r t 1 m 1 τ r t α ,
where an integer m > 0 . From the above, for t , g T × G 2 , the next inequality follows
G e 2 Γ 1 α l n τ r t τ r t α m   G e m a x 2 τ r m a x t 0 1 m 1 τ r m a x t 0 ,
i.e., the second summand of G M t , g α (A54) is bounded, too.
Since, for 0 < β 0   α 1 gamma function Γ 1 α Γ 1 β 0 > 1 , the last two elements of the gradient (A53) are nonnegative definite (c.f., Property 1) and bounded by 1 2 τ r m a x t 0 and G e m a x 2 t 0 τ r m a x t 0 , respectively, for any t , g T × G 2 . The theorem is proved.

Appendix B

Appendix B.1. The Results of the Numerical Studies for Material I

Table A1. The elements α N , β N , G e N and τ r N of the FMM parameter vector g N solving identification task (25) for real relaxation modulus (36) of the material described by the unimodal Gauss-like distribution, the mean-square identification indices Q N g N , Equation (24), the mean relative square model approximation index Q N r e l g N , Equation (38), the sampling points-independent integral indices Q g N defined by the optimization task (26), and the relative errors E R R (37) of the FMM parameter g approximation for N relaxation modulus measurements independently disturbed by additive, zero mean, normally distributed noises with standard deviation σ = 2   k P a .
Table A1. The elements α N , β N , G e N and τ r N of the FMM parameter vector g N solving identification task (25) for real relaxation modulus (36) of the material described by the unimodal Gauss-like distribution, the mean-square identification indices Q N g N , Equation (24), the mean relative square model approximation index Q N r e l g N , Equation (38), the sampling points-independent integral indices Q g N defined by the optimization task (26), and the relative errors E R R (37) of the FMM parameter g approximation for N relaxation modulus measurements independently disturbed by additive, zero mean, normally distributed noises with standard deviation σ = 2   k P a .
N Q N g N   k P a 2 Q N r e l g N Q g N   k P a 2 E R R   [ % ] α N β N G e N   k P a τ r N s
509.295186 × 10−52.102093 × 10−32.863765 × 10−30.5270.9637888.644275 × 10−22.86257315.170858
1005.614324 × 10−45.907274 × 10−51.350707 × 10−31.2240.9656147.207174 × 10−22.74521916.225154
2002.185904 × 10−41.109724 × 10−44.959376 × 10−34.537 × 10−30.9547258.815551 × 10−23.06597614.289093
5005.559271 × 10−45.053629 × 10−55.262647 × 10−47.740 × 10−30.9152321.261495 × 10−23.11387912.68090
10004.594139 × 10−41.476206 × 10−45.268387 × 10−43.104 × 10−40.9209371.341017 × 10−23.08128512.926718
20004.668059 × 10−43.243456 × 10−65.396114 × 10−45.315 × 10−20.9290992.375304 × 10−23.01556313.557098
50005.289533 × 10−43.282759 × 10−75.213213 × 10−42.876 × 10−40.9205021.445854 × 10−23.08148912.984136
70005.224292 × 10−43.978418 × 10−65.215398 × 10−41.654 × 10−40.9203641.419721 × 10−23.08275412.979817
10,0005.082251 × 10−43.739909 × 10−65.219846 × 10−49.301 × 10−50.9203271.391317 × 10−23.08374612.942422
12,0005.150186 × 10−41.363350 × 10−55.208965 × 10−41.689 × 10−50.9201871.504288 × 10−23.08545512.942422
15,0005.247854 × 10−43.687418 × 10−65.205439 × 10−41.938 × 10−80.9200141.472034 × 10−23.08668012.949456
Table A2. The elements α N , β N , G e N and τ r N of the FMM parameter vector g N solving identification task (25) for real relaxation modulus (36) of the material described by the unimodal Gauss-like distribution, the mean-square identification indices Q N g N , Equation (24), the mean relative square model approximation index Q N r e l g N , Equation (38), the sampling points-independent integral indices Q g N defined by the optimization task (26), and the relative errors E R R (37) of the FMM parameter g approximation for N measurements independently disturbed by additive, zero mean, normally distributed noises with standard deviation σ = 5   k P a .
Table A2. The elements α N , β N , G e N and τ r N of the FMM parameter vector g N solving identification task (25) for real relaxation modulus (36) of the material described by the unimodal Gauss-like distribution, the mean-square identification indices Q N g N , Equation (24), the mean relative square model approximation index Q N r e l g N , Equation (38), the sampling points-independent integral indices Q g N defined by the optimization task (26), and the relative errors E R R (37) of the FMM parameter g approximation for N measurements independently disturbed by additive, zero mean, normally distributed noises with standard deviation σ = 5   k P a .
N Q N g N   k P a 2 Q N r e l g N Q g N   k P a 2 E R R   [ % ] α N β N G e N   k P a τ r N s
501.246174 × 10−43.1182072 × 10−32.953417 × 10−30.5010.9620598.720788 × 10−22.86828715.170858
1006.042548 × 10−49.788569 × 10−39.505292 × 10−40.5420.9507115.528005 × 10−22.85944115.193338
2002.469407 × 10−42.622478 × 10−44.888504 × 10−33.999 × 10−30.9545638.747687 × 10−23.06725014.289093
5005.754599 × 10−41.136117 × 10−75.288353 × 10−41.1543 × 10−40.9177991.594113 × 10−23.09003712.837417
10004.8724914 × 10−46.769619 × 10−45.261736 × 10−42.129 × 10−40.9211421.374989 × 10−23.08221912.926718
20004.8687842 × 10−43.040833 × 10−75.400918 × 10−45.345 × 10−20.9292342.394509 × 10−23.01536513.557098
50005.4973567 × 10−43.940686 × 10−65.211862 × 10−43.026 × 10−40.9204721.459525 × 10−23.08135412.984135
70005.401956 × 10−41.559930 × 10−65.214093 × 10−41.573 × 10−40.9203501.428889 × 10−23.08285212.979817
10,0005.261243 × 10−42.744246 × 10−65.221499 × 10−45.516 × 10−40.9209051.446157 × 10−23.07947312.979817
12,0005.356062 × 10−45.782985 × 10−55.208668 × 10−42.386 × 10−50.9201221.495904 × 10−23.08521612.942422
15,0005.457918 × 10−44.535517 × 10−65.205517 × 10−42.513 × 10−70.9199801.477095 × 10−23.08656912.949456
Table A3. The elements α N , β N , G e N and τ r N of the FMM parameter vector g N solving identification task (25) for real relaxation modulus (36) of the material described by the unimodal Gauss-like distribution, the mean relative square model approximation index Q N r e l g N , Equation (38), the sampling points-independent integral indices Q g N defined by the optimization task (26), and the relative errors E R R (37) of the parameter g approximation for N measurements independently disturbed by additive, zero mean, normally distributed noises with standard deviation σ = 8   k P a .
Table A3. The elements α N , β N , G e N and τ r N of the FMM parameter vector g N solving identification task (25) for real relaxation modulus (36) of the material described by the unimodal Gauss-like distribution, the mean relative square model approximation index Q N r e l g N , Equation (38), the sampling points-independent integral indices Q g N defined by the optimization task (26), and the relative errors E R R (37) of the parameter g approximation for N measurements independently disturbed by additive, zero mean, normally distributed noises with standard deviation σ = 8   k P a .
N Q N g N   k P a 2 Q N r e l g N Q g N   k P a 2 E R R   [ % ] α N β N G e N   k P a τ r N s
501.823164 × 10−44.314509 × 10−32.993805 × 10−30.5040.9615658.769949 × 10−22.86764115.1708579
1006.647523 × 10−40.3535919.411604 × 10−40.5170.9498575.428199 × 10−22.86482215.193338
2002.961081 × 10−44.334879 × 10−44.686224 × 10−34.277 × 10−30.9542038.597868 × 10−23.06658214.289093
5006.078210 × 10−42.946700 × 10−55.271033 × 10−43.193 × 10−40.9180841.579814 × 10−23.09223812.837417
10005.328032 × 10−42.789884 × 10−35.267889 × 10−41.300 × 10−40.9213251.356372 × 10−23.08320312.926718
20005.248407 × 10−43.026647 × 10−75.403471 × 10−45.202 × 10−20.9292972.403479 × 10−23.01632013.557098
50005.883023 × 10−41.431097 × 10−55.210845 × 10−43.179 × 10−40.9204421.473303 × 10−23.08121912.984135
70005.759879 × 10−41.780341 × 10−75.220610 × 10−42.087 × 10−30.9216581.538542 × 10−23.07262213.066791
10,0005.616302 × 10−42.127802 × 10−65.220589 × 10−45.578 × 10−40.9208551.448935 × 10−23.07943312.979817
12,0005.739756 × 10−42.163690 × 10−45.211071 × 10−48.574 × 10−40.9211091.594595 × 10−23.07768513.006321
15,0005.843841 × 10−45.116492 × 10−65.207537 × 10−44.857 × 10−40.9208231.552819 × 10−23.07992113.006321

Appendix B.2. The Results of the Numerical Studies for Material II

Table A4. For the optimal FMM approximating the relaxation modulus (43) of the material described by the BSW spectrum (42): the elements α N β N , G e N and τ r N of the vector g N solving identification task (25), the mean-square identification indices Q N g N , Equation (24), the mean relative square model approximation index Q N r e l g N , Equation (38), the sampling points-independent integral indices Q g N defined by the optimization task (26), and the relative errors E R R (37) of the parameter g approximation for N relaxation modulus measurements independently disturbed by additive normally distributed noises with standard deviation σ = 3   k P a .
Table A4. For the optimal FMM approximating the relaxation modulus (43) of the material described by the BSW spectrum (42): the elements α N β N , G e N and τ r N of the vector g N solving identification task (25), the mean-square identification indices Q N g N , Equation (24), the mean relative square model approximation index Q N r e l g N , Equation (38), the sampling points-independent integral indices Q g N defined by the optimization task (26), and the relative errors E R R (37) of the parameter g approximation for N relaxation modulus measurements independently disturbed by additive normally distributed noises with standard deviation σ = 3   k P a .
N Q N g N   M P a 2 Q N r e l g N Q g N   M P a 2 E R R   [ % ] α N β N G e N   M P a τ r N s
501.315712 × 10−53.597318 × 10−75.638335 × 10−52.3420.6776826.870627 × 10−21.3705855.418555 × 103
1001.161929 × 10−55.298049 × 10−74.939762 × 10−54.1510.6561346.493121 × 10−21.4098995.094185 × 103
2001.03161 × 10−51.124491 × 10−85.287861 × 10−54.7590.650046.373051 × 10−21.4214095.002011 × 103
5009.772230 × 10−62.140316 × 10−82.966257 × 10−50.8970.6864757.329572 × 10−21.3362615.791700 × 103
10001.216504 × 10−51.472521 × 10−82.964808 × 10−50.8470.6917237.332173 × 10−21.3316435.808836 × 103
20009.462435 × 10−62.397636 × 10−93.703687 × 10−51.9280.6747096.955942 × 10−21.3649645.509183 × 103
50003.372717 × 10−51.236704 × 10−82.439945 × 10−50.0610.7503368.307524 × 10−21.2450236.555588 × 103
70003.499392 × 10−53.0361484 × 10−92.578132 × 10−50.1790.7611538.483605 × 10−21.2312096.669047 × 103
10,0003.974638 × 10−51.327619 × 10−92.524376 × 10−50.1360.7579998.429268 × 10−21.2353416.633616 × 103
12,0003.289041 × 10−52.540537 × 10−102.384835 × 10−51.384 × 10−30.7358398.062010 × 10−21.2654476.373837 × 103
15,0003.259757 × 10−52.789747 × 10−102.383509 × 10−56.536 × 10−50.7374898.098389 × 10−21.2625616.402808 × 103
Table A5. For the optimal FMM approximating the relaxation modulus (43) of the material described by the BSW spectrum (42): the elements α N , β N , G e N and τ r N of the parameter vector g N , the mean-square identification indices Q N g N , Equation (24), the mean relative square model approximation index Q N r e l g N , Equation (38), the sampling points-independent integral indices Q g N defined by the optimization task (26), and the relative errors E R R (37) of the parameter g for N measurements corrupted by the noises with standard deviation σ = 6   k P a .
Table A5. For the optimal FMM approximating the relaxation modulus (43) of the material described by the BSW spectrum (42): the elements α N , β N , G e N and τ r N of the parameter vector g N , the mean-square identification indices Q N g N , Equation (24), the mean relative square model approximation index Q N r e l g N , Equation (38), the sampling points-independent integral indices Q g N defined by the optimization task (26), and the relative errors E R R (37) of the parameter g for N measurements corrupted by the noises with standard deviation σ = 6   k P a .
N Q N g N   M P a 2 Q N r e l g N Q g N   M P a 2 E R R   [ % ] α N β N G e N   M P a τ r N s
505.334837 × 10−51.607199 × 10−65.383509 × 10−51.7830.6915837.062629 × 10−21.3526475.543339 × 103
1004.625811 × 10−51.912785 × 10−65.226179 × 10−54.8530.6523636.407151 × 10−21.4212374.988252 × 103
2004.104462 × 10−56.4067158 × 10−86.964808 × 10−57.8410.6322935.905374 × 10−21.4679874.606199 × 103
5003.664178 × 10−57.819189 × 10−83.074233 × 10−51.0390.6810717.285674 × 10−21.3428915.745449 × 103
10003.879589 × 10−56.729167 × 10−82.961505 × 10−50.8690.6898127.328929 × 10−21.3332335.801105 × 103
20003.618881 × 10−51.385717 × 10−93.604767 × 10−51.7860.6752996.992477 × 10−21.3620995.542684 × 103
50005.999448 × 10−53.264396 × 10−82.456599 × 10−58.109 × 10−20.7520938.335898 × 10−21.2425586.579817 × 103
70006.1067925.065819 × 10−92.598612 × 10−50.1950.7632588.506666 × 10−21.2292296.680137 × 103
10,0006.592012 × 10−51.717204 × 10−92.536338 × 10−50.1440.7594018.444666 × 10−21.2340446.640191 × 103
12,0005.916358 × 10−52.831521 × 10−92.384641 × 10−55.855 × 10−40.7371508.077941 × 10−21.2640566.382156 × 103
15,0005.893401 × 10−51.984301 × 10−102.384002 × 10−52.644 × 10−40.7382838.108685 × 10−21.2617016.408038 × 103
Table A6. For the optimal FMM approximating the relaxation modulus (43) of the material described by the BSW spectrum (42): the elements α N , β N , G e N and τ r N of the parameter vector g N , the mean-square identification indices Q N g N , Equation (24), the mean relative square model approximation index Q N r e l g N , Equation (38), the sampling points-independent integral indices Q g N defined by the optimization task (26), and the relative errors E R R (37) of the parameter g for N measurements corrupted by the noises with standard deviation σ = 8   k P a .
Table A6. For the optimal FMM approximating the relaxation modulus (43) of the material described by the BSW spectrum (42): the elements α N , β N , G e N and τ r N of the parameter vector g N , the mean-square identification indices Q N g N , Equation (24), the mean relative square model approximation index Q N r e l g N , Equation (38), the sampling points-independent integral indices Q g N defined by the optimization task (26), and the relative errors E R R (37) of the parameter g for N measurements corrupted by the noises with standard deviation σ = 8   k P a .
N Q N g N   M P a 2 Q N r e l g N Q g N   M P a 2 E R R   [ % ] α N β N G e N   M P a τ r N s
509.520744 × 10−52.941612 × 10−66.384002 × 10−51.4480.7015867.198997 × 10−21.3402515.627729 × 103
1008.218302 × 10−53.329357 × 10−65.485285 × 10−55.3480.6498886.349515 × 10−21.4288674.918087 × 103
2007.290485 × 10−51.233057 × 10−78.330015 × 10−510.3520.6207455.578770 × 10−21.5011444.339259 × 103
5006.446603 × 10−51.355064 × 10−73.166061 × 10−51.1450.6774847.255754 × 10−21.3474285.713136 × 103
10006.646278 × 10−51.240244 × 10−72.966288 × 10−50.8850.6885287.326539 × 10−21.3343195.795691 × 103
20006.395292 × 10−52.673029 × 10−93.542823 × 10−51.6930.6756997.016833 × 10−21.3601855.565162 × 103
50008.738813 × 10−55.134265 × 10−82.469399 × 10−59.607 × 10−20.7532788.354811 × 10−21.2409156.595934 × 103
70008.835062 × 10−56.705427 × 10−92.613160 × 10−50.2050.7646778.522041 × 10−21.2279136.687353 × 103
10,0009.319389 × 10−52.003354 × 10−92.544649 × 10−50.1490.7603318.454841 × 10−21.2331886.644512 × 103
12,0008.640779 × 10−56.119895 × 10−92.384921 × 10−52.494 × 10−40.7380148.088431 × 10−21.2631456.387532 × 103
15,0008.625405 × 10−51.518245 × 10−102.384484 × 10−54.544 × 10−40.7387928.115159 × 10−21.2611576.411273 × 103

References

  1. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications; Springer: Berlin/Heidelberg, Germany, 2014. [Google Scholar] [CrossRef]
  2. Schiessel, H.; Metzler, R.; Blumen, A.; Nonnejunacher, T.F. Generalized viscoelastic models: Their fractional equations with solutions. J. Phys. A Math. Gen. 1995, 28, 6567–6584. [Google Scholar] [CrossRef]
  3. Metzler, R.; Nonnenmacher, T.F. Fractional relaxation processes and fractional rheological models for the description of a class of viscoelastic materials. Int. J. Plast. 2003, 19, 941–959. [Google Scholar] [CrossRef]
  4. Mainardi, F. Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models; World Scientific: Hackensack, NJ, USA, 2022. [Google Scholar]
  5. Mainardi, F.; Spada, G. Creep, Relaxation and Viscosity Properties for Basic Fractional Models in Rheology. Eur. Phys. J. Spec. Top. 2011, 193, 133–160. [Google Scholar] [CrossRef]
  6. Meng, R.; Yin, D.; Drapaca, C.S. A variable order fractional constitutive model of the viscoelastic behavior of polymers. Int. J. Non-Linear Mech. 2019, 113, 171–177. [Google Scholar] [CrossRef]
  7. Shiru, G.; Hai, P.; Mengsen, Y.; Ziheng, S.; Dejun, L.; Jinyong, X.; Yiying, F. New perspective on the creep characteristic of fiber–dependent shape memory polymers: Variable–order fractional constitutive model. J. Mater. Res. Technol. 2023, 24, 5177–5185. [Google Scholar] [CrossRef]
  8. Liang, Z.; Li, J.; Zhang, X.; Kan, Q. A viscoelastic-viscoplastic constitutive model and its finite element implementation of amorphous polymers. Polym. Test. 2023, 117, 107831. [Google Scholar] [CrossRef]
  9. Xu, B.; Blok, R.; Teuffel, P. An investigation of the effect of relative humidity on viscoelastic properties of flax fiber reinforced polymer by fractional-order viscoelastic model. Compos. Commun. 2023, 37, 101406. [Google Scholar] [CrossRef]
  10. Viviani, L.; Di Paola, M.; Royer-Carfagni, G. Fractional viscoelastic modeling of laminated glass beams in the pre-crack state under explosive loads. Int. J. Solids Struct. 2022, 248, 111617. [Google Scholar] [CrossRef]
  11. Xiang, G.; Yin, D.; Meng, R.; Cao, C. Predictive model for stress relaxation behavior of glassy polymers based on variable-order fractional calculus. Polym. Adv. Technol. 2021, 32, 703–713. [Google Scholar] [CrossRef]
  12. Popov, I.I.; Shitikova, M.V.; Levchenko, A.V.; Zhukov, A.D. Experimental identification of the fractional parameter of the fractional derivative standard linear solid model for fiber-reinforced rubber concrete. Mech. Adv. Mater. Struct. 2023, 1–9. [Google Scholar] [CrossRef]
  13. Hajikarimi, P.; Ehsani, M.; EL Haloui, Y.; Tehrani, F.F.; Absi, J.; Nejad, F.M. Fractional viscoelastic modeling of modified asphalt mastics using response surface method. Constr. Build. Mater. 2022, 317, 125958. [Google Scholar] [CrossRef]
  14. Cai, W.; Wang, P. Rate-dependent fractional constitutive model for nonlinear behaviors of rubber polymers. Eur. J. Mech. A/Solids 2024, 103, 105186. [Google Scholar] [CrossRef]
  15. David, S.A.; Katayama, A.H. Fractional Order for Food Gums: Modeling and Simulation. Appl. Math. 2013, 4, 28194. [Google Scholar] [CrossRef]
  16. Stankiewicz, A. Fractional Maxwell model of viscoelastic biological materials. BIO Web Conf. 2018, 10, 02032. [Google Scholar] [CrossRef]
  17. Bainy, E.M.; Lenzi, E.K.; Corazza, M.L.; Lenzi, M.K. Mathematical Modeling of Fish Burger Baking Using Fractional Calculus. Therm. Sci. 2017, 21, 41–50. [Google Scholar] [CrossRef]
  18. Lei, D.; Liang, Y.; Xiao, R. A fractional model with parallel fractional Maxwell elements for amorphous thermoplastics. Phys. A Stat. Mech. Its Appl. 2018, 490, 465–475. [Google Scholar] [CrossRef]
  19. Hou, T.; Chen, H. Isothermal physical aging of PEEK and PPS investigated by fractional Maxwell model. Polymer 2012, 53, 2509–2518. [Google Scholar] [CrossRef]
  20. Elmoghazy, Y.H.; Safaei, B.; Asmael, M.; Sahmani, S.; Zeeshan, Q.; Qin, Z. Computational Modelling and Analysis of Effect of Viscoelastic Materials on Damping and Vibrational Behaviors of Composite Structures—An Extensive Review. Arch. Computat. Methods Eng. 2024, 1–52. [Google Scholar] [CrossRef]
  21. Gritsenko, D.; Paoli, R. Theoretical Analysis of Fractional Viscoelastic Flow in Circular Pipes: General Solutions. Appl. Sci. 2020, 10, 9093. [Google Scholar] [CrossRef]
  22. Loglio, G.; Kovalchuk, V.I.; Bykov, A.G.; Ferrari, M.; Krägel, J.; Liggieri, L.; Miller, R.; Noskov, B.A.; Pandolfini, P.; Ravera, F.; et al. Interfacial Dilational Viscoelasticity of Adsorption Layers at the Hydrocarbon/Water Interface: The Fractional Maxwell Model. Colloids Interfaces 2019, 3, 66. [Google Scholar] [CrossRef]
  23. Kovalchuk, V.I.; Loglio, G.; Bykov, A.G.; Ferrari, M.; Krägel, J.; Liggieri, L.; Miller, R.; Milyaeva, O.Y.; Noskov, B.A.; Ravera, F.; et al. Effect of Temperature on the Dynamic Properties of Mixed Surfactant Adsorbed Layers at the Water/Hexane Interface under Low-Gravity Conditions. Colloids Interfaces 2020, 4, 27. [Google Scholar] [CrossRef]
  24. Mao, B.; Zhu, R.; Wang, Z.; Yang, Y.; Han, X.; Zhao, Q. A Modified Fractional Maxwell Numerical Model for Constitutive Equation of Mn-Cu Damping Alloy. Materials 2020, 13, 2020. [Google Scholar] [CrossRef]
  25. Schmidt, A.; Gaul, L. Parameter Identification and FE Implementation of a Viscoelastic Constitutive Equation Using Fractional Derivatives. Proc. Appl. Math. Mech. 2002, 1, 153–154. [Google Scholar] [CrossRef]
  26. Xiao, Z.; Haitian, Y.; Yiqian, F. Identification of constitutive parameters for fractional viscoelasticity. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 311–322. [Google Scholar] [CrossRef]
  27. Zhou, H.W.; Yi, H.Y.; Mishnaevsky, L.; Wang, R.; Duan, Z.Q.; Chen, Q. Deformation analysis of polymers composites: Rheological model involving time-based fractional derivative. Mech. Time-Depend. Mater. 2017, 21, 151–161. [Google Scholar] [CrossRef]
  28. Zhou, S.; Cao, J.; Chen, Y. Genetic Algorithm-Based Identification of Fractional-Order Systems. Entropy 2013, 15, 1624–1642. [Google Scholar] [CrossRef]
  29. Lewandowski, R.; Chorążyczewski, B. Identification of the parameters of the Kelvin–Voigt and the Maxwell fractional models, used to modeling of viscoelastic dampers. Comput. Struct. 2010, 88, 1–17. [Google Scholar] [CrossRef]
  30. Beda, T.; Chevalier, Y. New Methods for Identifying Rheological Parameter for Fractional Derivative Modeling of Viscoelastic Behavior. Mech. Time-Depend. Mater. 2004, 8, 105–118. [Google Scholar] [CrossRef]
  31. Yin, Y.; Yang, Z.; Shi, M. Circular arc rules of complex plane plot for model parameters determination of viscoelastic material. Mech. Time-Depend. Mater. 2021, 25, 631–643. [Google Scholar] [CrossRef]
  32. Ljung, L. System Identification: Theory for the User, 1st ed.; Prentice Hall: Englewood Cliffs, NJ, USA, 1999. [Google Scholar]
  33. Stankiewicz, A.; Juściński, S. How to Make the Stress Relaxation Experiment for Polymers More Informative. Polymers 2023, 15, 4605. [Google Scholar] [CrossRef]
  34. Ljung, L. Convergence analysis of parametric identification methods. IEEE Trans. Autom. Control 1978, 23, 770–783. [Google Scholar] [CrossRef]
  35. Hasiewicz, Z.; Stankiewicz, A. On input-dependent system identification by Monte Carlo approach. IEEE Trans. Autom. Control 1986, 31, 192, Erratum in IEEE Trans. Autom. Control 1985, 30, 480–483. [Google Scholar] [CrossRef]
  36. Apelblat, A. Differentiation of the Mittag-Leffler Functions with Respect to Parameters in the Laplace Transform Approach. Mathematics 2020, 8, 657. [Google Scholar] [CrossRef]
  37. Muzeau, E.; Perez, J.; Johari, G.P. Mechanical spectrometry of the beta-relaxation in poly(methyl methacrylate). Macromolecules 1991, 24, 4713–4723. [Google Scholar] [CrossRef]
  38. Carrillo-Navas, H.; Hernández-Jaimes, C.; Utrilla-Coello, R.G.; Meraz, M.; Vernon-Carter, E.J.; Alvarez-Ramirez, J. Viscoelastic relaxation spectra of some native starch gels. Food Hydrocoll. 2014, 37, 25–33. [Google Scholar] [CrossRef]
  39. Wang, J.; Wang, X.; Ruan, H. On the mechanical β relaxation in glass and its relation to the double-peak phenomenon in impulse excited vibration at high temperatures. J. Non-Cryst. Solids 2020, 533, 119939. [Google Scholar] [CrossRef]
  40. Baumgaertel, M.; Schausberger, A.; Winter, H.H. The relaxation of polymers with linear flexible chains of uniform length. Rheol. Acta 1990, 29, 400–408. [Google Scholar] [CrossRef]
  41. Baumgaertel, M.; Winter, H.H. Interrelation between continuous and discrete relaxation time spectra. J. Non-Newton. Fluid Mech. 1992, 44, 15–36. [Google Scholar] [CrossRef]
  42. Cho, K.S. Power series approximations of dynamic moduli and relaxation spectrum. J. Rheol. 2013, 57, 679–697. [Google Scholar] [CrossRef]
  43. Lee, S.H.; Bae, J.-E.; Cho, K.S. Determination of continuous relaxation spectrum based on the Fuoss-Kirkwood relation and logarithmic orthogonal power-series approximation. Korea-Aust. Rheol. J. 2017, 29, 115–127. [Google Scholar] [CrossRef]
  44. Ferry, J.D. Viscoelastic Properties of Polymers, 3rd ed.; John Wiley & Sons: New York, NY, USA, 1980. [Google Scholar]
  45. Gorenflo, R.; Mainardi, F. Fractional Calculus: Integral and Differential Equations of Fractional Order. In Fractals and Fractional Calculus in Continuum Mechanics; Carpinteri, A., Mainardi, F., Eds.; Springer: Wien, Austria; New York, NY, USA, 1997; pp. 223–276. [Google Scholar]
  46. Malkin, A.I.A.; Malkin, A.Y.; Isayev, A.I. Rheology: Concepts, Methods and Applications; ChemTec: Deerfield Beach, FL, USA, 2006; Available online: https://books.google.pl/books?id=8rGafjhgz-UC (accessed on 5 April 2023).
  47. Stankiewicz, A. On Applicability of the Relaxation Spectrum of Fractional Maxwell Model to Description of Unimodal Relaxation Spectra of Polymers. Polymers 2023, 15, 3552. [Google Scholar] [CrossRef] [PubMed]
  48. Pérez-Calixto, D.; Amat-Shapiro, S.; Zamarrón-Hernández, D.; Vázquez-Victorio, G.; Puech, P.-H.; Hautefeuille, M. Determination by Relaxation Tests of the Mechanical Properties of Soft Polyacrylamide Gels Made for Mechanobiology Studies. Polymers 2021, 13, 629. [Google Scholar] [CrossRef] [PubMed]
  49. Andrei, N. Modern Numerical Nonlinear Optimization; Springer: Cham, Switzerland, 2022. [Google Scholar] [CrossRef]
  50. Hoeffding, W. Probability inequalities for sums of bounded random variables. J. Am. Stat. Assoc. 1963, 58, 13–30. [Google Scholar] [CrossRef]
  51. Kwakye-Nimo, S.; Inn, Y.; Yu, Y.; Wood-Adams, P.M. Linear viscoelastic behavior of bimodal polyethylene. Rheol. Acta 2022, 61, 373–386. [Google Scholar] [CrossRef]
  52. Cirillo, G.; Spizzirri, U.G.; Iemma, F. Functional Polymers in Food Science: From Technology to Biology, Volume 1: Food Packaging; Wiley: New York, NY, USA, 2015; Available online: https://onlinelibrary.wiley.com/doi/book/10.1002/9781119108580 (accessed on 20 June 2023).
  53. Lorenzo, G.; Checmarev, G.; Zaritzky, N.; Califano, A. Linear viscoelastic assessment of cold gel-like emulsions stabilized with bovine gelatin. LWT Food Sci. Technol. 2011, 44, 457–464. [Google Scholar] [CrossRef]
  54. Bardet, S.; Gril, J. Modelling the transverse viscoelasticity of green wood using a combination of two parabolic elements. C. R. Mec. 2002, 330, 549–556. [Google Scholar] [CrossRef]
  55. Kurenuma, Y.; Nakano, T. Analysis of stress relaxation on the basis of isolated relaxation spectrum for wet wood. J. Mater. Sci. 2012, 47, 4673–4679. [Google Scholar] [CrossRef]
  56. Choi, J.; Cho, K.S.; Kwon, M.K. Self-Similarity and Power-Law Spectra of Polymer Melts and Solutions. Polymers 2022, 14, 3924. [Google Scholar] [CrossRef]
  57. Unidad, H.J.; Goad, M.A.; Bras, A.R.; Zamponi, M.; Faust, R.; Allgaier, J.; Pyckhout-Hintzen, W.; Wischnewski, A.; Richter, D.; Fetters, L.J. Consequences of Increasing Packing Length on the Dynamics of Polymer Melts. Macromolecules 2015, 48, 6638–6645. [Google Scholar] [CrossRef]
  58. Podlubny, I. Fitting Data Using the Mittag-Leffler Function. MATLAB Central File Exchange. Available online: https://www.mathworks.com/matlabcentral/fileexchange/32170-fitting-data-using-the-mittag-leffler-function (accessed on 23 July 2023).
  59. Podlubny, I. Mittag-Leffler Function with Two Parameters. MATLAB Central File Exchange. Available online: https://www.mathworks.com/matlabcentral/fileexchange/8738-mittag-leffler-function (accessed on 23 July 2023).
  60. Stankiewicz, A. Two-Level Scheme for Identification of the Relaxation Time Spectrum Using Stress Relaxation Test Data with the Optimal Choice of the Time-ScaleFactor. Materials 2023, 16, 3565. [Google Scholar] [CrossRef]
  61. Stankiewicz, A. A Class of Algorithms for Recovery of Continuous Relaxation Spectrum from Stress Relaxation Test Data Using Orthonormal Functions. Polymers 2023, 15, 958. [Google Scholar] [CrossRef] [PubMed]
  62. Fikhtengol’ts, G.M. The Fundamentals of Mathematical Analysis; International Series of Monographs in Pure and Applied Mathematics; Pergamon Press: Oxford, UK; London, UK, 2013; Volume 2. [Google Scholar]
  63. Lemes, N.H.T.; dos Santos, J.P.C.; Braga, J.P. A generalized Mittag-Leffler function to describe nonexponential chemical effects. Appl. Math. Model. 2016, 40, 7971–7976. [Google Scholar] [CrossRef]
  64. Deming, W.; Colcord, C. The Minimum in the Gamma Function. Nature 1935, 135, 917. [Google Scholar] [CrossRef]
  65. Topsøe, F. Some bounds for the logarithmic function. In Inequality Theory and Applications; Cho, Y.J., Kim, J.K., Dragomir, S.S., Eds.; Nova Science Publishers: New York, NY, USA, 2007; p. 137. [Google Scholar]
  66. Alzer, H. On Some Inequalities for the Gamma and Psi Functions. Math. Comput. 1997, 66, 373–389. [Google Scholar] [CrossRef]
  67. Mező, I.; Hoffman, M.E. Zeros of the digamma function and its Barnes G-function analogue. Integral Transform. Spec. Funct. 2017, 28, 846–858. [Google Scholar] [CrossRef]
Figure 1. Viscoelastic models: (a) classic Maxwell model; (b) fractional Scott-Blair model of an order α ; (c) fractional Maxwell model; elastic modulus G e , G e 1 , G e 2 , viscosity η , relaxation times τ r , τ 1 , τ 2 .
Figure 1. Viscoelastic models: (a) classic Maxwell model; (b) fractional Scott-Blair model of an order α ; (c) fractional Maxwell model; elastic modulus G e , G e 1 , G e 2 , viscosity η , relaxation times τ r , τ 1 , τ 2 .
Materials 17 01527 g001
Figure 2. Dependence of the parameters of the FMM approximating the “real” relaxation modulus (36): (a) α N ; (b) β N ; (c) G e N ; and (d) τ r N on the number of measurements N for disturbances σ = 2 , 5 , 8   P a ; the horizontal purple lines are related to the optimal parameters α , β , G e and τ r independent on the sampling instants used in the rheological experiment.
Figure 2. Dependence of the parameters of the FMM approximating the “real” relaxation modulus (36): (a) α N ; (b) β N ; (c) G e N ; and (d) τ r N on the number of measurements N for disturbances σ = 2 , 5 , 8   P a ; the horizontal purple lines are related to the optimal parameters α , β , G e and τ r independent on the sampling instants used in the rheological experiment.
Materials 17 01527 g002
Figure 3. The indices of the “real” relaxation modulus (36) approximation by the FMM: (a) the mean-square empirical index Q N g N , Equation (24), (b) the integral quadratic sampling instants-independent index Q g N , Equation (26), as the functions of the number of measurements N and noises σ = 2 , 5 , 8   P a ; the horizontal purple lines correspond to the optimal integral index Q g defined in Equation (27).
Figure 3. The indices of the “real” relaxation modulus (36) approximation by the FMM: (a) the mean-square empirical index Q N g N , Equation (24), (b) the integral quadratic sampling instants-independent index Q g N , Equation (26), as the functions of the number of measurements N and noises σ = 2 , 5 , 8   P a ; the horizontal purple lines correspond to the optimal integral index Q g defined in Equation (27).
Materials 17 01527 g003
Figure 4. The measurements G ¯ i (red points) of the “real” relaxation modulus (36) and optimal FMM models: sampling points-independent G M t , g and empirical G M t , g N for N measurements and normal distribution noises with the standard deviation σ = 8   P a : (a) N = 100 ; (b) N = 10,000 .
Figure 4. The measurements G ¯ i (red points) of the “real” relaxation modulus (36) and optimal FMM models: sampling points-independent G M t , g and empirical G M t , g N for N measurements and normal distribution noises with the standard deviation σ = 8   P a : (a) N = 100 ; (b) N = 10,000 .
Materials 17 01527 g004
Figure 5. Dependence of the mean indices of the “real” relaxation modulus G t (36) optimal approximation by the FMM: (a) relative empirical error E R R Q N r e l (39), (b) integral error E R R Q (40) on the number of measurements N and the noises’ standard deviations σ .
Figure 5. Dependence of the mean indices of the “real” relaxation modulus G t (36) optimal approximation by the FMM: (a) relative empirical error E R R Q N r e l (39), (b) integral error E R R Q (40) on the number of measurements N and the noises’ standard deviations σ .
Materials 17 01527 g005
Figure 6. Dependence of the mean relative error M E R R (41) between the optimal parameters g N and g of the FMM approximating the “real” relaxation modulus G t (36) on the number of measurements N and the noises’ standard deviation σ .
Figure 6. Dependence of the mean relative error M E R R (41) between the optimal parameters g N and g of the FMM approximating the “real” relaxation modulus G t (36) on the number of measurements N and the noises’ standard deviation σ .
Materials 17 01527 g006
Figure 7. The parameters of the FMM approximating the relaxation modulus (43) of material described by the BSW relaxation spectrum (42): (a) α N ; (b) β N ; (c) G e N ; and (d) τ r N as the functions of the number of measurements N for noises σ = 3 , 6 , 8   k P a ; the horizontal purple lines correspond to optimal model parameters α , β , G e and τ r being independent on the sampling instants used in the experiment.
Figure 7. The parameters of the FMM approximating the relaxation modulus (43) of material described by the BSW relaxation spectrum (42): (a) α N ; (b) β N ; (c) G e N ; and (d) τ r N as the functions of the number of measurements N for noises σ = 3 , 6 , 8   k P a ; the horizontal purple lines correspond to optimal model parameters α , β , G e and τ r being independent on the sampling instants used in the experiment.
Materials 17 01527 g007
Figure 8. The indices of the BSW relaxation modulus (42), (43) approximation by the FMM: (a) the mean-square empirical index Q N g N , Equation (24), (b) the integral quadratic sampling instants-independent index Q g N , Equation (26), as the functions of the number of measurements N for the noises σ = 2 , 5 , 8   P a ; the horizontal purple lines correspond to the optimal integral index Q g defined in Equation (27).
Figure 8. The indices of the BSW relaxation modulus (42), (43) approximation by the FMM: (a) the mean-square empirical index Q N g N , Equation (24), (b) the integral quadratic sampling instants-independent index Q g N , Equation (26), as the functions of the number of measurements N for the noises σ = 2 , 5 , 8   P a ; the horizontal purple lines correspond to the optimal integral index Q g defined in Equation (27).
Materials 17 01527 g008
Figure 9. The measurements G ¯ i (red points) of the real relaxation modulus (43) of the material described by the BSW spectrum (42) and the fractional Maxwell optimal models: sampling points-independent G M t , g and empirical G M t , g N for N measurements and additive random normally distributed noises with standard deviation = 8   P a and zero mean value: (a) N = 100 ; (b) N = 10,000 .
Figure 9. The measurements G ¯ i (red points) of the real relaxation modulus (43) of the material described by the BSW spectrum (42) and the fractional Maxwell optimal models: sampling points-independent G M t , g and empirical G M t , g N for N measurements and additive random normally distributed noises with standard deviation = 8   P a and zero mean value: (a) N = 100 ; (b) N = 10,000 .
Materials 17 01527 g009
Figure 10. Dependence of the mean indices of the “real” BSW relaxation modulus (42), (43) approximation by the FMM: (a) the mean relative empirical error E R R Q N r e l (39), (b) the mean optimal sampling points-independent integral error E R R Q (40) on the number of measurements N and the noises’ standard deviation σ .
Figure 10. Dependence of the mean indices of the “real” BSW relaxation modulus (42), (43) approximation by the FMM: (a) the mean relative empirical error E R R Q N r e l (39), (b) the mean optimal sampling points-independent integral error E R R Q (40) on the number of measurements N and the noises’ standard deviation σ .
Materials 17 01527 g010
Figure 11. Dependence of the mean relative error M E R R (41) between the parameters g and g N of the FMM approximating the “real” BSW relaxation modulus (42), (43) on the number of measurements N and the noises’ standard deviation σ .
Figure 11. Dependence of the mean relative error M E R R (41) between the parameters g and g N of the FMM approximating the “real” BSW relaxation modulus (42), (43) on the number of measurements N and the noises’ standard deviation σ .
Materials 17 01527 g011
Table 1. The components α , β , G e and τ r of the FMM parameter g solving the optimal identification problem (27) and the optimal integral quadratic indices Q g defined by (27) for the “real” relaxation modulus G t (36).
Table 1. The components α , β , G e and τ r of the FMM parameter g solving the optimal identification problem (27) and the optimal integral quadratic indices Q g defined by (27) for the “real” relaxation modulus G t (36).
Q g   k P a 2 α β G e   k P a τ r s
5.2054279 × 10−40.9200291.469033 × 10−23.08672312.949456
Table 2. The elements α , β , G e and τ r of the FMM parameter g solving the optimal identification task (27) and the optimal integral quadratic indices Q g defined by (27) for the “real” relaxation modulus G t (42), (43).
Table 2. The elements α , β , G e and τ r of the FMM parameter g solving the optimal identification task (27) and the optimal integral quadratic indices Q g defined by (27) for the “real” relaxation modulus G t (42), (43).
Q g   M P a 2 α β G e   M P a τ r s
2.383349 × 10−50.7367068.088257 × 10−21.26341256.397636 × 103
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Stankiewicz, A. Sampling Points-Independent Identification of the Fractional Maxwell Model of Viscoelastic Materials Based on Stress Relaxation Experiment Data. Materials 2024, 17, 1527. https://doi.org/10.3390/ma17071527

AMA Style

Stankiewicz A. Sampling Points-Independent Identification of the Fractional Maxwell Model of Viscoelastic Materials Based on Stress Relaxation Experiment Data. Materials. 2024; 17(7):1527. https://doi.org/10.3390/ma17071527

Chicago/Turabian Style

Stankiewicz, Anna. 2024. "Sampling Points-Independent Identification of the Fractional Maxwell Model of Viscoelastic Materials Based on Stress Relaxation Experiment Data" Materials 17, no. 7: 1527. https://doi.org/10.3390/ma17071527

APA Style

Stankiewicz, A. (2024). Sampling Points-Independent Identification of the Fractional Maxwell Model of Viscoelastic Materials Based on Stress Relaxation Experiment Data. Materials, 17(7), 1527. https://doi.org/10.3390/ma17071527

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop