Generalised viscoelastic fibre at small strain

A straight elastic fibre is usually perceived as a one-dimensional structural component, and its similarity with a cylindrical rod makes its concept analogous, if not equivalent with the concept of an elastic spring. This analogy enables this communication to match the one-dimensional response of a relevant viscoelastic fibre with that of a viscoelastic spring and, hence, to describe its one-dimensional behaviour in the light of a new, generalised viscoelastic spring model. The model shares simultaneously properties of an elastic spring and an inelastic damper (dashpot) and this communication is interested on its applicability at small strain only. However, the form of its constitutive equation, which is based on the combined action of an internal energy function and a viscous flow potential, is non-linear as well as differential and, also, implicit in the stress. The model enables a posteriori determination of (i) the manner that the elastic and the inelastic parts of the fibre strain are assembled and form the observed total deformation, (ii) the part of stress that creates recoverable work and the part of stress wasted in energy dissipation, and (iii) the amount of work stored in the material as well as the amount of energy dissipation during the fibre deformation. A detailed analysis is presented for the case that small-strain, steady viscoelastic deformation takes place in a spatially homogeneous manner. This includes a complete relevant solution of the problem of interest and is accompanied by an adequate set of corresponding qualitative numerical results.

The known analogy between a a linearly elastic spring of elastic modulus K, and b a linearly elastic cylindrical rod of Young's modulus E, cross-sectional area A, and length L; f , σ and e represent force, axial stress and elastic strain, respectively ical tissue [12,13] were followed by substantial relevant progress, and also led to identification of fibre response that diverges considerably from that of the strong fibres considered in [1][2][3][4][5][6][7][8][9][10][11]. The existence and behavioural modelling of biological fibres that resist extension but do not support compression is thus already attracting considerable attention (see [14][15][16] and references therein). The same though does not seem to be true in the case that soft fibres exhibit viscoelastic behaviour (e.g. [17][18][19]). This communication thus aims to initiate a relevant modelling endeavour.
It is initially recalled that a straight elastic fibre is usually perceived as having the form of a flexible cord (e.g. [1,2]) or that of a slender cylindrical rod (e.g. [8]), or even the form of an Euler-Bernoulli beam in cases that the fibre can resist bending (e.g. [6]). In all cases though, it is macroscopically regarded as some kind of a onedimensional structure and, in this regard, its similarity with a straight cylindrical rod makes its concept analogous, if not equivalent to that of an elastic spring (e.g. [20,21] and Fig. 1). That analogy enables this communication to connect the one-dimensional response of a relevant viscoelastic fibre with that of the generalised viscoelastic spring model presented in [21]. It also enables adoption of a simplicity kind of notion employed in [21], according which the term "spring" is considered interchangeable/replaceable with the term "rod" and, hence, "fibre" in what follows, without further consequences. This notion is in fact unquestionable in elasticity, at least within the small strain/deformation regime ( Fig. 1) which will become the range of applicability of the present analysis.
It is then further recalled that the generalised viscoelastic spring/rod/fibre model presented in [21] shares simultaneously properties of an elastic spring and an inelastic damper (dashpot). At its most general form, that model describes the behaviour of one-dimensional deformable solids subjected to large viscoelastic deformation but, through a series of progressively simplifying considerations, it produced a simpler relevant model that is applicable at small strain and deformation. That small-strain viscoelastic model is a very appropriate, essentially ideal point of departure for the more general research endeavour mentioned above, and thus becomes the focal point of this communication.
Under these considerations, and with the help of an Appendix, Sect. 2 presents the small-strain viscoelastic spring/fibre model in a manner that enables it to stand by itself, independently of its large deformation counterpart(s) developed in [21]. Nevertheless, the model is still underpinned by the combined action of a viscous flow potential that controls the inelastic part of its response and an internal energy function that dictates the manner that energy is stored in the material due to the elastic part of its total deformation. Although the model is considered applicable at small strain only, its constitutive equation still obtains a non-linear form, which is also differential as well as implicit in the stress, regardless of the form of either of the involved potentials. However, influenced by linear elasticity considerations, Sect. 3 combines next a quadratic form of the internal energy function with a quadratic viscous flow potential and, apart from presenting the relevant form obtained by the constitutive equation, it also discusses some fundamental energy and dissipation matters and describes in detail the manner that the elastic and the inelastic parts of the spring/fibre deformation are determined. Section 4 achieves further analytical simplification by confining attention to the type of homogeneous deformation underpinned by the solution of the corresponding linearly elastic spring/fibre problem. In this manner, it introduces in a rational manner the case of steady or strictly quasi-static viscoelastic spring/fibre deformation. It is noted in passing that that point may alternatively be reached, independently of the modelling route followed in Sects. 2 and 3, through a relevant linearised procedure followed and detailed in [21]. With the help of a second Appendix, Sect. 4 then presents the exact solution of the considered linear problem in the form of a fast converging series. Appropriate use of that series solution thus underpins the presentation of a set of qualitative numerical results that reveal the manner that the total deformation splits into elastic and inelastic parts.
Section 5 continuous the analysis and demonstrates the manner that (i) a part of the created energy is stored in the material while another part is dissipated, and (ii) the total stress acting on the spring/fibre splits into work-recoverable and energy-dissipative parts. Discussions of the outlined analytical results and corresponding computational progress continue into Sect. 6, which also highlights some important research steps and directions that may follow within well-known relevant research areas.
2 The generalised viscoelastic spring/rod/fibre at small strain 2.1 Preliminaries: the linearly elastic spring/fibre The linearly elastic spring and rod/fibre depicted in Fig. 1 are considered equivalent one-dimensional elastic components, in the sense that they have the same length, L, they are both clamped at one end and, when the same external force of magnitude f is applied at their free end, they undergo identical one-dimensional deformation. The spring has elastic modulus K while its cylindrical rod/fibre equivalent has Young's modulus E and cross-sectional area A. Through this well-known equivalence which, with the help of a single co-ordinate parameter, x, is adequately explained in the figure, these two linearly elastic models and associated terminology are considered completely interchangeable in what follows, without further notice or consequence.
It is accordingly assumed that, at time t t 0 , the spring/fibre is unloaded and unstressed at length L. Loading from that initial configuration takes off through application of a force, and, for t > t 0 , gives rise to small strain e u ,x , (2.2) where u(x, t) represents the displacement and a comma denotes partial differentiation with respect to x. As is further well-known, the strain energy stored internally in the elastically deforming spring, is quadratic in the strain and positive semi-definite; the equality sign applies only in the case of zero strain. The quadratic form of this strain energy function thus gives rise to the linear constitutive equation where σ is the stress. Moreover, equilibrium requires where a dot denotes differentiation with respect to time, v u represents the velocity and ρ is the material density. It is recalled, for later use, that changes of material density are generally considered negligible in linear elasticity and, at this point, ρ is accordingly considered identical to its initial value, ρ 0 , attained at t t 0 .
It is pointed out that Fig. 1 indicates that the boundary traction, T (t), is tensile and, therefore, positive. For convenience, and unless is stated otherwise, this notion of tensile and positive T (t) will be considered valid everywhere in what follows. It is left understood though that T (t) may also be considered compressive, by appropriately reversing the sign of the affected analytical results, wherever this is necessary.

2.2
The generalised, small-strain viscoelastic spring/fibre It is now postulated that the implied one-dimensional deformation comprises simultaneous action and/or interaction of elastic and inelastic deformation, the later being due to viscous fluid-type flow. The manner that those different parts of material response are merged and compose the observed total strain, e, is initially considered unknown.
It is, instead, anticipated that, due to inevitable energy dissipation, the elastic energy, W , stored internally in the viscoelastically deforming spring changes with, and therefore depends not only on the total strain, e, but also on the rate of the inelastic strain,ė v , namely whereė v ė u ,x v ,x . Moreover, such a dismissal of the purely elastic nature of the spring implies that potential changes of its material density may not anymore be negligible, even within the small-strain regime. The mass density evolution of the implied small-strain viscoelastic spring is thus required to also obey the standard continuity equatioṅ It is emphasised that the total strain involved in (2.6) is still defined in (2.2) but is not anymore elastic. On the other hand, the appearing rate of inelastic strain is defined as follows: where a prime denotes ordinary differentiation with respect to σ , and, as is shown in [21], the existence of a viscous flow potential φ(σ ) guarantees that the rate of the dissipative work, is stationary.
Here, the stress is decomposed into two parts as follows: where σ R is the part that produces recoverable/elastic work while σ D is wasted in creation of dissipative work with rate ψ. Nevertheless, (2.5) still serves as the spring equation of motion (see also [21] as well as Appendix A). The most general form of a constitutive equation associated with the outlined small-strain viscoelastic model is as follows: This is obtained in [21] through appropriate reduction of its large deformation counterpart into the small-strain deformation regime. As is though shown in Appendix A for self-sufficiency of this communication, it is also alternatively obtainable through direct use of appropriate small-deformation postulates/considerations. Despite its small-strain regime of applicability, this is generally a non-linear differential equation for σ , even in cases that its solution might be sought and found independently of the solution of the equilibrium and the continuity equations (2.5) and (2.7a), respectively.
In summary, the small-strain viscoelastic spring/rod/fibre model can take the form of the following onedimensional initial/boundary value problem: given the form of both the strain energy function (2.6) and the flow potential introduced in (2.8), determine the three principal unknowns, u, ρ and σ , in a manner that satisfies the three principal differential equations, (2.5), (2.7a) and (2.11), subject to the boundary conditions (2.12) and the initial conditions (2.7b) and After the implied principal unknowns thus are determined, and (2.2) is used for the calculation of the total strain e, the inelastic and elastic parts of the latter, as well as the elastic part of the displacement, can, respectively, be obtained as follows: (2.14)

Quadratic potentials
It is natural for someone to search for simpler, possibly linearised forms of (2.11) by reducing as much as possible the influence that viscous flow exerts on the strain energy function, W . Influenced by the quadratic form that strain energy functions, such as (2.3), meet in linear elasticity, the strain energy function, W , as well as the viscous potential, φ(σ ), are now both considered quadratic in their arguments. It thus is initially assumed that where η is an appropriate, viscosity related material modulus. An additional quadratic term that might involve (ė v ) 2 is not included in (3.1), as this would represent purely viscous material flow that is assumed of negligible influence on the stored internal energy. The form (3.1) of W clearly satisfies the initial condition (2.6b) and enables the constitutive equation (2.11) to acquire the more specific form and, by virtue of (2.8), It is pointed out that, if η 0, this constitutive equation resembles closely its linearly elastic spring counterpart (2.4) which, however, assumes that the material density is not affected by the deformation.
The second term within the parenthesis appearing in right-hand side of either (3.2) or (3.3) reveals that this constitutive equation still preserves the non-linear differential equation features observed already in (2.11). However, further simplification becomes possible by next confining attention to the lowest polynomial form that the flow potential may attend without violation the initial condition (2.6b). This is the special case of a quadratic viscous flow potential, where G is some appropriate material modulus.
In that case, the viscous flow rule (2.8) yieldṡ and enables (3.1) to obtain the special form where use is made of the stress decomposition rule (2.10), along with the fact that, by definition, σ D contributes to dissipative work only. It is also observed that the material parameters η and G, introduced in (3.1) and (3.4), respectively, are here combined to form In what follows, this combined material parameter,η, will thus be referred to as the generalised viscosity of the model associated with the quadratic flow potential (3.4). Introduction of the latter thus enables the constitutive equation (3.3) to obtain the form which, due to involvement of the ratioσ /ė , is still differential with respect to time as well as non-linear in both σ and e.

Energy and dissipation considerations
By virtue of definition (2.9), the energy dissipation encountered during the described inelastic deformation is expected to be the following non-decreasing function of time: where σ D is to be determined, and the equality sign holds in the particular case that the spring behaves elastically (ė v σ D 0). The following connection is then established between the internal energy (3.6) of the viscoelastic spring and the strain function (2.3) of its elastic counterpart: (3.10) By virtue of (3.9), this requires which, through direct use of the internal energy expressions (2.3) and (3.6), as well the notion (3.7), leads to On the other hand, (3.10) also implies that the developed viscoelastic model may stop to represent a solid spring/fibre at some time instant, say t t F > t 0 , that While the undeformed state, e 0, serves as an essentially trivial solution of this equation, the alternative solution estimates the ratio of the total strain and the work-recoverable part of stress at the time, t t F , that viscous effects may overwhelm the observed deformation and, hence, deprive from the spring its ability to store elastic energy.

Determination of the elastic and inelastic parts of deformation and stress
Solution of the simultaneous differential equations (2.5), (2.7) and (3.8), subject to the boundary and initial conditions (2.12) and (2.13), respectively, will enable determination of the displacement, the mass density and the stress distribution within the framework outlined so far in this section. With the stress, σ , and the total strain, e, being thus determined, and by virtue of (3.4), (2.14) will then provide the inelastic and elastic strain parts as well as the elastic part of the displacement as respectively. It is pointed out that expressions (3.15) do not guarantee that the total strain, e, and its resulting elastic and inelastic parts, e e and e v , are necessarily all of the same sign. The parts σ D and σ R of the stress that create dissipative and recoverable work, respectively, may then also be estimated in a couple of different ways. The direct and more accurate route followed later (see Sect. 5) is through time differentiation and subsequent comparison of (3.9) and (3.12a). This leads to the first of the relations the second being a result of the stress decomposition (2.10). It is anticipated that at least one of σ D and σ R should be of the same sign with σ . At this stage of the analysis though, no apparent reason seems to suggest that all three σ , σ D and σ R should be of the same sign throughout the spring span.
There is also provided in [21] an alternative approach for the estimation of σ D and σ R . However, this will not be applied in the present problem of interest, as it is regarded as an approximate approach that may become applicable to more advanced/complicated generalised viscoelastic spring cases.

Spatially homogeneous viscoelastic deformation
It is now recalled that the deformation of a corresponding elastic spring is spatially homogeneous in the sense that, at any given time instance, t, the strain, e, is not only elastic, but also constant throughout the spring span (e.g. [20]). The same may not necessarily be true in the present viscoelastic spring case but, even if it is not, it can still be considered as an adequate approximation of the problem solution, at least for sufficiently small viscoelastic deformations.
Attention is accordingly now directed into the search for potential solutions of the outlined problem that are consistent with a displacement field of the form where α(t) is some suitable function of time that needs to be determined. This displacement pattern satisfies the boundary condition (2.12a) and underpins the homogeneous deformation field Moreover, it enables the continuity equation to obtain the form of the separable ordinary differential equation whose solution yields the material density of the spring as follows: Introduction of these results into the constitutive Eq. (3.8) converts the latter into the following: Differentiation of this equation with respect to the spatial co-ordinate parameter, x, followed by appropriate use of the equation of motion (2.5), leads next tô and, by virtue of (4.1) and (4.4), tô Potential solutions of this highly non-linear ODE will provide forms of the function α(t) that are consistent with the search for spatially homogeneous viscoelastic spring deformations.

Steady (strictly quasi-static) viscoelastic deformation
Attention is next confined into the simplest time-dependent solution of (4.6), namely the linear function This imposes conditions of steady viscoelastic deformation, e t/t 0 − 1, (4.8) that takes place at constant speed, v 1/t 0 , and thus enables the equation of motion (2.5) to attain naturally its quasi-static form The stress thus is independent of position but, of course, is still dependent on time. By virtue of the boundary condition (2.12b), σ is thus equal to the currently unknown boundary traction, T (t), that drives the assumed deformation, (4.1) and (4.7). As it is considered unphysical for the magnitude of either σ R or σ D to exceed the magnitude of T (t) at the boundary, it is accordingly concluded that both σ R and σ D have necessarily the same sign with σ throughout the span of the spring.
Upon adopting at this point the reasonable assumption that the directions of total stress and total strain coincide throughout the spring span, at least within the small-strain regime, one concludes that this is also the case with the direction of either of σ R or σ D . It is thus seen that and as (3.12) thus requireŝ η ηG < 0, (4.11) the material moduli η and G should be having opposite signs. It follows by virtue of (3.5) that, if then a spring element in extension (σ > 0) experiences inelastic strain of an increasing rate, while in compression (σ < 0) the rate of inelastic strain will be decreasing. This effect is evidently reversed it in the alternative case, where However, as σ equals the boundary traction throughout the spring span, the whole spring is either in extension or in compression, and this thus clearly depends on the direction of T (t). For simplicity, and by virtue of the T -direction shown in Fig. 1, only the tensile spring problem will be considered in what follows, although most, if not all of the outlined analytical results hold regardless of whether the spring is in tension or in compression. Moreover, wherever detailed specification of G and η is required for numerical applications, only the case (4.12) will be considered. This is because, as is shown later in Sect. 5 with the help of an Appendix, consideration of the alternative case (4.13) leads to values of σ R that violate the necessary condition (4.10). where Solution of this equation, subject to the initial condition (2.13b), provides the stress in the closed form where, however, the noted integrations may be accomplished only numerically. Nevertheless, an alternative form of this solution is possible, and this is provided in Appendix B, in the form (B.1) of an infinite power series expansion in e. The coefficients of that series are given by the recurrence formula (B.5), which reveals features of very fast convergence within the regime of small strain. As a matter of fact, the first term of that series is proportional to e and, as the recurrence formula shows that the second term is of the order of e 2 , the order of each subsequent term decreases very rapidly.
A single term truncation of (B.1), namely then provides an excellent approximation of the exact solution (4.16), in the sense that it produces numerical results that are practically unaffected by the contribution of the truncated higher-order terms. Indeed, the straight lines depicted in Fig. 2 are all based on (4.17) and remain unaffected if re-plotted by retaining more terms in (B.1). As will be seen later in Sect. 5 though, the corresponding work-recoverable and energy-dissipative stress parts vary with the strain in a non-linear manner.  Fig. 1. Moreover, care is taken for the total strain, e, to vary within reasonably small bounds. It is pointed out in this regard that the so-called infinitesimal strain regime may be narrower to the small-strain range employed in this communication, as is sometimes subjected to an upper-bound restriction that does not exceed 5%. Nevertheless, the infinitesimal strain regime is still part of the small-strain range employed in what follows and, if necessary, the associated discussion and relevant observations can be restricted accordingly. It should be emphasised though that a decrease or even an increase of the small strain upper limit implied in Figs. 2, 3, 4 and 5 leaves completely unaffected the qualitative features of the presented numerical results.
Under these considerations, the straight lines depicted in Fig. 2 exhibit the practically linear manner that the non-dimensional stress parameter σ σ/E e/ 1 + 2 η (4.18) Fig. 5 Elastic strain e e 1 − e v (1 − λe)e against total strain e, for the three higher values of λ employed in Fig. 3 varies with increasing the total strain e, for several values of the viscosity parameter η . Figure 2 thus clearly demonstrates the decreasing effect that increasing viscosity exerts on the stress value. It is emphasised that the non-dimensionalisation employed in (4.18) is based on the Young modulus, E, of the corresponding elastic spring/fibre, rather than its effective viscoelastic spring counterpart,Ê, defined in (4.17b). As a matter of fact, if E were replaced in (4.18) byÊ, then all straight lines shown in Fig. 2 would appropriately rotate and become identical with their depicted elastic spring counterpart, obtained by setting in (4.18) η 0. It becomes thus understood that any set of numerical results involving below the non-dimensional stress parameter σ/Ê is similarly a contracted representation of a wider set of relevant results.
In view of the stress distributions depicted in Fig. 2, it may erroneously be felt at this point that the effect of viscosity is confined within a simple replacement of the actual spring modulus, E, by its effective counterpart,Ê, or, equivalently, the replacement of a linearly elastic spring/fibre with another that possesses some appropriately reduced elastic modulus. However, as will be seen in what follows, the present viscoelastic model provides also a considerable amount of further interesting information regarding the manner that (i) the total strain splits into elastic and inelastic parts, (ii) the total stress splits into work-recoverable and energy-dissipative parts, and (iii) the spring internal energy is accumulated or dissipated during the assumed viscoelastic deformation.
Under these considerations, the inelastic and the elastic parts of the total strain e, as well as the elastic part of the total displacement, are next obtained by inserting (4.18) into (3.15) and performing the noted integrations. These are as follows: , (4.20) appears naturally at this stage of the solution as a dimensionless, inelastic strain normalisation parameter. This parameter has necessarily the same sign with G and may thus be either positive or negative. However, it has already been anticipated that only choices of G and η that comply with (4.12) enable satisfaction of the necessary condition (4.10) and, hence, someone can safely consider that λ > 0, η < 0, (4.21) in what follows.
For several values of the normalisation parameter λ, Fig. 3 thus demonstrates the parabolic manner that the inelastic part of the strain, e v , varies with increasing the total strain, e; see also (4.19a) and (4.8), respectively. It is observed that, for relatively small values of λ, the inelastic part of strain remains smaller than e. This is the case for the three lower λ-values, for each one of which Fig. 4 depicts further the corresponding variation of the elastic strain part; the e e variation that corresponds to each of the three higher values of λ is similarly depicted in Fig. 5. For the low λ-values, Figs. 3 and 4 thus show that both e v and e e remain positive and, therefore, tensile within the employed region of total strain variation.
However, for the three higher values of λ, there comes a time instance that the inelastic strain reaches the value of e and, hence, the elastic strain becomes zero. This evidently happens when the relevant parabolic curve reaches the bisector of the co-ordinate axes depicted in Fig. 3. Beyond that crossing point, the positive value of the inelastic strain tends to exceed its total strain counterpart (e v > e), as is demonstrated in Fig. 3 for the two higher values of λ. However, despite the temporary zero value that e e attains at that point, the elastic part of the deformation is still active due to the non-zero value of the strain energy stored already in the material. Figure 5 thus shows that when the spring is loaded beyond the implied crossing point, the stored strain energy enables the elastic strain to become compressive (change of sign) and to thus counteract the tensile part of e v that tends to exceed the total e-value.

Strain energy and energy dissipation: work-recoverable and energy-dissipative parts of stress
The form (3.6) of the strain energy function implies that determination of the energy stored in the deforming viscoelastic spring requires previous evaluation of the stress part σ R that produces recoverable/elastic work. It is accordingly observed that a combination of (3.5) with (3.9) leads tȯ where use is also made of (2.10). On the other hand, time differentiation of (3.12a) yieldṡ Equating the right-hand sides of (5.1) and (5.2), and making also use of (4.11), one obtains With use of (4.8) and (4.17a), this equation acquires the alternative forṁ This is essentially in the standard form (4.14) of a first-order linear ODE and, hence, (5.4) it also susceptible to the standard closed form solution (4.16).
Nevertheless, an equivalent power series solution of (5.4) is obtained in Appendix C. This is practically again more useful than (4.16), because (i) it makes it clear that only the case (4.12) needs to be considered in this investigation, (ii) enables formation of the following closed form relationship between the total stress and its work-recoverable counterpart: and (iii) thus makes evident the natural appearance of the new, non-dimensional normalisation parameter It is pointed out that, by virtue of (3.11) and (4.12), the larger (smaller) is the positive value ofλ , the stronger is the influence of the elastic (inelastic) effects of the deformation on the work-recoverable part of the stress.
The fact that both σ R and σ D have already been anticipated tensile and, therefore, positive enables (2.10) to provide the condition: If connected with (5.5), this inequality provides from the range of values ofλe to satisfy the following condition: 1 λe It can readily be verified that this condition is indeed satisfied in the present case of interest, where bothλ and e are considered positive.
By virtue of (2.10) and (5.5), the energy-dissipative part of the stress relates with the total stress as follows: Use of (5.5) and (5.9) thus enables Fig. 6 to demonstrate the highly non-linear manner that the ratios σ R /σ and σ D /σ vary within a substantially large region ofλe-values.
The curves drawn in that figure make it evident that inelastic effects dominate in regions of relatively smallλevalues. However, their influence is continually degreasing with increasingλe and is dying out withλe approaching infinity. Elastic effects thus take gradually over and become dominant for larger values of theλe-parameter.
It is next seen that, with use of (5.5), (4.7) and (4.11), the internal energy function (3.6) may take the form which, by virtue of (5.8), is clearly a decreasing function of the strain, e. This expression demonstrates the manner that the viscoelastic spring modulus, E v , changes with the deformation. Based on this formula, Fig. 7 demonstrates that manner that the ratio W (e,ė v )/W e decreases with increasing theλe-parameter, for several values of the generalised viscosity parameter,η. It is pointed out that, at the limiting case thatλe approaches zero, a double application of L' Hospital's rule reveals that and, hence, that the value of the ratio W (e,ė v )/W e is everywhere finite and well defined. However, a direct computational evaluation of the large ratio-expression appearing in (5.11) creates the wrong impressions that the left-hand side of (5.12) is indeterminate. The resulting numerical instability thus prevents the graphical user interface employed for creation of the plots depicted in Fig. 7 from producing reliable numerical results in the vicinity ofλe 0. However, the numerical results depicted in Fig. 8 make it clear that nothing unusual, unexpected or, more generally, singular takes place at the vicinity ofλe 0. These numerical results (Fig. 8) hold within in the chosen, relatively small-value region ofλe (0 ≤λe ≤ 1) for identical material parameters to those employed in Fig. 7. However, the depicted curves of the ratio W (e,ė v )/W e are drawn through use of a four-term truncation of the series solution (C.6) of (5.4), rather than its mathematically equivalent closed form (5.5).
Under these considerations, it is safe for someone to consider that all curves depicted in Fig. 7 depart from the top left corner of that figure. Hence, either of (5.10) or (5.11) can finally reveal that the value of the corresponding dissipated energy is (5.13) as can alternatively be verified by combining (5.5) and (3.12a).

Further discussion and conclusions
It is emphasised that the steady viscoelastic spring/fibre deformation considered in Sect. 4, and discussed afterwards extensively is Sect. 5, becomes possible only when the boundary condition (2.12b) is applied in consistence with (4.17) and, therefore, when the externally applied tensile traction evolves in time as follows: The relevant solution detailed in Sect. 4 then represents an exceptional viscoelastic spring/fibre type of deformation that is followed neither by stress relaxation nor by strain creep. It is, however, recalled that the stress relaxation or strain creep behaviour of a standard spring-dashpot viscoelastic constitutive model is usually studied by subjecting the model into some known, usually constant amount of strain or stress, respectively (e.g. [23][24][25]).
If such an essentially arbitrary loading condition is similarly applied externally on the present viscoelastic model, the constitutive law (3.8) will evidently fail to respond in a manner consistent with the equation of motion (2.5) or the continuity equation (2.7). It is in fact anticipated that, even if the spatially homogenous deformation conditions (4.1) were potentially applicable in such a case, the constitutive law (3.8) would probably fail to satisfy the alternative relevant form (4.6) attained in that case by the equation of motion. Most of these comments and observations hold evidently true in relevant mathematical models that do not introduce quasi-static conditions leading to (4.9) in the strict manner considered in this communication, but still employ it, in the usual manner, as a reliable approximation of corresponding boundary value problems (see also [26]).
The thus emerging equilibrium unbalance would then set up a stress relaxation and/or a strain creep process that would continue until mass density balance and stress equilibrium were both completely restored. It is accordingly concluded that the exact viscoelasticity solution presented in Sect. 4 is further expected to represent the final stage of stress relaxation and/or strain creep that follow any relevant small-strain viscoelastic deformations that takes earlier place in some unsteady (not quasi-static) manner.
Alternatively, an initially unsteady small-strain viscoelastic deformation that conforms with the spatially homogeneous deformation assumptions considered at the beginning of Sect. 4 should be followed by some appropriate stress relaxation and/or strain creep process that would enable it to approach as close as possible either (i) the relevant steady equilibrium state detailed in Sect. 4.1 or (ii) some alternative, nearer to it equilibrium state that, if available, will also satisfy the continuity equation (2.7) and the reduced form (4.6) of the equation of motion.
It is further observed with interest that, despite its non-linear form, the constitutive equation (3.8) can be rearranged in the following alternative form: This resembles closely either its well-known linear elastic spring/rod counterpart (2.4) where, however, the appearing elastic modulus, E, is constant, or its homogeneous, steady deformation viscoelastic approximation (4.17), where the effective (or reduced) constant modulusÊ emerges as an excellent approximation of the essentially unknown time-dependent modulusẼ implied in (6.2).
In the present case though, the value of the effective Young's modulusÊ depends on the generalised viscosity parameterη. In this regard, it is remarkable that the deformation-dependent modulusẼ appearing in (6.2) is so accurately approximated in (4.17) by the its constant counterpartÊ. It is in fact remarkable that, in the small-strain regime of interest, (4.17) replaces so accurately the essentially non-linear deformation response of a viscoelastic spring with what, at a first view, might be thought of as the response of a linearly elastic spring (see also the relevant straight lines depicted in Fig. 2 for η  0).
Nevertheless, the more thorough relevant investigation performed in Sects. 4 and 5 revealed that this seemingly linear material behaviour is composed by well-specified elastic and inelastic parts whose relevant strain and stress states vary non-linearly with the implied small-strain viscoelastic deformation. Moreover, the analysis detailed in Sect. 5 made possible a mathematically exact determination of the dissipated energy and the work stored in the material during the implied deformation process.
It may still be pointed out though that the analytical solution detailed in Sect. 4, as well as all relevant numerical results presented and discussed afterwards, refer only to continuous loading conditions, starting from an initial state at which the spring/fibre is both unstressed and unstrained. The fact that no relevant unloading process was considered is due to that such a process may be defined, in advance, in one of several different ways. Indeed, even in the simple case that such a spring unloading process is anticipated elastic, someone still needs to specify the relevant constant value of the spring/fibre elastic modulus.
Choosing, for instance, any of the straight lines depicted in Fig. 2 for η 0 and assuming that unloading from some point of that line takes place elastically, one has still to assign a value to the corresponding elastic (Young's) modulus. If this is assumed equal to the valueÊ reached already during the outlined viscoelastic loading process, then elastic unloading will evidently take place by tracing the same line "downwards", in the opposite direction. If, on the other hand, there is some reason (such as experimental evidence) suggesting that elastic unloading takes place under the influence of the modulus E of the corresponding "virgin" elastic spring, then unloading would trace downwards a straight line parallel to that depicted in Fig. 2 for η 0. The choice of any other elastic or viscoelastic, spatially homogeneous, steady unloading process would end up with use of some different effective or reduced Young's modulus and would thus trace downwards some different straight line. In case that such an unloading process is considered viscoelastic, the slope of that line would evidently depend not only on the anticipated value of the Young's modulus but also on the value of its associated generalised viscosity parameter,η. It is anticipated that, while the fibre will still be in tension during unloading, either of these material moduli may not necessarily be identical to their counterparts employed during the corresponding loading process.
Under these considerations, it is mentioned in passing that by (i) keeping E constant during several successive loading-unloading cycles of a relevant deformable spring/fibre and, as an effect of each of those cycles, (ii) decreasing successively the influence of the generalised viscosity parameterη, a constitutive equation of the type (3.8) may succeed to model successfully the Mullins hardening effect observed in the behaviour of several kinds of polymers and standard rubber-like materials (e.g. [27][28][29]). However, at these early modelling stages, this observation is only regarded as a possibility that requires thorough investigation and further study. This is at least because such a further study/investigation will certainly require use of a relevant constitutive equation that is valid within the large rather than the small strain/deformation regime (see also [21]).
It is noted in this context that the linearly inelastic spring model described in Sect. 4, and the corresponding analytical solution detailed in Sect. 5, may well guide the search for relevant large deformation solutions associated with, and stemming from a proper and detailed consideration of the fully non-linear inelastic spring model detailed in [21]. Such solutions are of paramount importance for a better understanding of the large inelastic deformation of one-dimensional structures, though may still be considered inadequate for relevant understanding of more realistic, three-dimensional inelastic material behaviour. Nevertheless, three-dimensional models that correspond to the nonlinear inelastic spring modelled in [21] are already available in the literature [30][31][32] and are expected to receive considerable benefit and boost from any relevant one-dimensional progress, such as the present one, that becomes gradually available. This is because the non-linear spring model developed in [21] and, hence, its small stain version considered in the present communication, may alternatively be obtained by specialising in one-dimension its three-dimensional counterpart detailed in [32]. It is recalled in this regard [21] that, by considering that the new plasticity-type postulations detailed in [32] are relevant and applicable to classical viscoelasticity situations, someone can clearly observe that all three subjects of Hyperelasticity, Plasticity and Viscoelasticity share a common theoretical background. The here presented small-strain one-dimensional model, analysis and relevant solution may thus be considered as the simplest possible example that fits within this larger, and still potential theoretical unification picture.