Nonlinear Models of Thermo-Viscoelastic Materials

The paper develops a general scheme for viscoelastic materials, where the constitutive properties are described by means of measures of strain, stress, heat flux, and their time derivatives. The constitutive functions are required to be consistent with the second law of thermodynamics. Indeed, a new view is associated with the second law: the non-negative expression of the entropy production is set equal to a further constitutive function. The introduction of the entropy production as a constitutive function allows for a much wider range of models. Within this range, a scheme to obtain nonlinear models of thermo-viscoelastic materials subject to large deformations is established. Notably, the Kelvin–Voigt, Maxwell, Burgers, and Oldroyd-B viscoelastic models, along with the Maxwell–Cattaneo heat conduction, are obtained as special cases. The scheme allows also for modelling the visco-plastic materials, such as the Prandtl–Reuss work-hardening function and the Bingham–Norton fluid.


Introduction
Viscoelasticity, especially the model of linear viscoelastic solids, traces back to Boltzmann [1], who, in essence, considered an elastic material with memory. The model elaborated by Boltzmann was based on the following assumptions.
At any point of the body, the stress at any time t depends upon the strain at all preceding times. If the strain up to time t is in the same direction, then the effect is to reduce the corresponding stress. The influence of a previous strain on the stress depends on the time elapsed since that strain occurred and is weaker for strains that occurred long ago. In addition, a superposition of the influence of previous strains holds, which means that the stress-strain relation is linear. Consequently, in the linear (infinitesimal) approximation, the Cauchy stress T is given by the infinitesimal strain ε in the form where K 0 and K(s) take values in the space of fourth-order tensors for any s ≥ 0 while T, K 0 , K, and ε are considered at fixed points of the body. The function K on [0, ∞) is called the Boltzmann function. To adhere to the standard notation, let G on [0, ∞) be defined by G is called the relaxation function, G 0 is the instantaneous elastic modulus, and G ∞ is the equilibrium elastic modulus. We can then write [2,3] It is assumed that, for solids, both G 0 and G ∞ and their difference G 0 − G ∞ are positive definite [2,4].
Hence time differentiation -denoted by a superposed dot -of (1) yieldṡ whence we have the Maxwell modelṪ The purpose of this paper is to establish thermodynamically consistent models of viscoelastic materials characterized by rate (or differential) equations for stress and strain. The approach is quite general in that the constitutive equations are nonlinear and threedimensional and involve finite deformations. Furthermore, the thermodynamic restrictions are determined so that coupling effects are allowed between stress rate, strain rate, and entropy production; formally, this is realized through a novel representation formula.
The memory effects of viscoelasticity are modeled jointly with thermal properties. In this respect, we observe that dissipative properties are ascribed to the materials models; there are approaches or models where dissipative effects are associated with external damping or boundary conditions [7].
Our approach is inherently thermodynamic. We derive constitutive vector-and tensorvalued rate equations consistent with thermodynamics in that they obey the second law (or Clausius-Duhem) inequality. This is the way we distinguish physically admissible models from purely mathematical ones, as is standard in rational thermodynamics. The pertinent variables of the constitutive equations comprise stress, strain, heat flux, and their time derivative.
It is a key point of our approach that the entropy production enters as a non-negative constitutive function. This in turn improves the generality of the materials models (see, for instance, [8]). Depending on the type of model, we proceed in the referential (Lagrangian) description or in the spatial (Eulerian) description. By arguing in the referential configuration, we obtain objective descriptions in the corresponding Eulerian description [9]. The main advantage of the present approach is the possibility of establishing nonlinear, thermodynamically consistent, objective rate-type relations describing memory and dissipative effects.
It is worth mentioning References [10,11], where viscoelastic models are developed in the relativistic Landau-Lifshitz frame on the basis of Onsager's linear non-equilibrium thermodynamics. The models are based on the energy momentum tensor; the particlenumber current; and where appropriate, the rheology equations. The constitutive equations are required to be consistent with the maximum property of the entropy functional in the process of linear regression. The maximum property of the entropy seems to be the analog of our requirement of non-negative entropy production. Instead, our scheme involves properties and processes such as nonlinearities and hysteretic effects [8,12,13].

Balance Laws and Entropy Inequality
We consider a body occupying the time-dependent region Ω ⊂ E 3 . The motion is described by means of the function χ(X, t), providing the position vector x ∈ Ω in terms of the position X, in a reference configuration R, and the time t, so that Ω = χ(R, t). The symbols ∇, ∇ R denote the gradient operator with respect to x ∈ Ω, X ∈ R. The function χ is assumed to be differentiable; hence, we can define the deformation gradient as F = ∇ R χ or, in suffix notation, F iK = ∂ X K χ i . The invertibility of X → x = χ(X, t) is guaranteed by letting J := det F > 0. Let v(x, t) be the velocity field, on Ω × R + . A superposed dot denotes time differentiation following the motion of the body, and hence, for any function f (x, t), we haveḟ = ∂ t f + v · ∇ f . We denote by L the velocity gradient, L ij = ∂ x j v i and recall thatḞ = LF.
Moreover D denotes the stretching tensor, D = SymL, and W the spin tensor W = SkewL. In terms of F, the (right) Cauchy-Green tensor C and the Green-Saint Venant deformation tensor E are defined by Let ε be the internal energy density (per unit mass), T be the symmetric Cauchy stress, q be the heat flux vector, ρ be the mass density, r be the (external) heat supply, and b be the mechanical body force per unit mass. The balance equations for mass, linear momentum, and energy are taken in the forṁ Let η be the entropy density and θ the absolute temperature. According the second law of thermodynamics, we take it that the inequality holds for any process compatible with the balance equations; σ is said to be the entropy production (per unit volume) [14]. Consequently, admissible constitutive equations are required to satisfy inequality (4). Multiplying by the absolute temperature θ and substituting ∇ · q − ρr from the energy Equation (3) 3 , we have Letting ψ = ε − θη (Helmholtz free energy density), we can write this inequality in the form Within the isothermal setting, the term θσ is usually referred to as the rate of (mechanical) dissipation [15]. Hence, non-dissipative models are characterized by σ = 0, the vanishing of the entropy production.
The modelling of the constitutive properties is made simpler by using referential, Euclidean invariant quantities. Let T RR is called the second Piola-Kirchhoff stress. Both T RR and q R are Euclidean invariant. Under any change of frame with rotation matrix Q, x * = c(t) + Q(t)x, det Q = 1, Lemma 1. Given a second-order tensor A, let N = A/|A|, |A| = √ A · A. If Z is a second-order tensor such that only Z · N is known, then Z = (Z · N)N + (I I I − N ⊗ N)G (14) for any second-order tensor G. A representation formula such as (14) also holds when A, Z and G are vectors, provided that I I I is replaced by 1, the second-order identity tensor.
Proof. If Z is a tensor such that Z is given while Z ⊥ is undetermined, then the representation of Z allows for any tensor Z ⊥ subject to Z ⊥ · N = 0. Now, since [(I I I − N ⊗ N)G] · N = 0 for any second-order tensor G, then [(I I I − N ⊗ N)G] is any possible value of Z ⊥ . Hence the relation Z = Z + Z ⊥ results in the representation (14).

Constitutive Relations
To describe thermo-viscoelastic effects, we let be the set of variables; in light of the model under consideration, the single variablesĖ,Ṫ RR , q R , ∇ R θ need not be mutually independent. Hence, we let ψ, η depend on Ξ R . In addition, we assume that η is continuous while ψ is continuously differentiable. Upon evaluation ofψ and substitution in (7), we obtain The linearity and arbitrariness ofθ, ∇ Rθ ,Ë,T RR ,q R imply that ψ is independent of ∇ R θ,Ė,Ṫ RR ,q R and hence so that the entropy inequality reduces to Further restrictions depend on the arbitrariness or possible constraints onĖ,Ṫ RR ,q R , ∇ R θ. In particular, since ψ is independent of ∇ R θ, ifĖ andṪ RR are related to each other but independent ofq R , ∇ R θ, then lettingq R = ∇ R θ = 0 we write (15) in the form where σ ET R is the entropy production density σ R whenq R = ∇ R θ = 0. In terms of the Gibbs free energy density, inequality (16) becomes Likewise, if, as is usually the case,q R , q R and ∇ R θ are independent ofĖ,Ṫ RR , then Apparently, σ R = σ ET R + σ Q R . Hereafter, the entropy productions σ ET R and σ Q R are assumed to be non-negative constitutive functions of Ξ R to be determined according to the constitutive model. (18) is common to many approaches where both q R and ∇ R θ are independent variables [16][17][18]. Since it is framed in the referential description then inequality (16) has the advantage that the material time derivative is objective [9,19,20]. Moreover, it makes consistency with thermodynamics much easier than it happens with histories [21] or summed histories [22].

Remark 1. Inequality
IfĖ,Ṫ RR andq R are independent of each other, then we have the following: Accordingly, ψ depends only on θ, E. In addition, T RR is no longer an independent variable but equals ρ R ∂ E ψ, which can be viewed as the constitutive relation for T RR , as it happens with thermo-hyperelastic materials.
It is of interest to look for the corresponding Eulerian descriptions. To this purpose, we observe that T (see (12)) and likewise q satisfy with q being the Truesdell rate of q.
Since ψ = ψ(θ, E, T RR , q R ), then we consider the dependence on E E E , T, q through E, T RR , q R . The definitions of E E E , T RR , and q R allow us to write Hence, in view of (12), we have and, in view of (10), Likewise, we evaluateq R , ∂ q R ψ and q R · ∇ R θ. Consequently, inequalities (16) and (18) can be written in the Eulerian form where σ ET and σ Q are non-negative constitutive functions of Ξ and σ ET + σ Q = σ.
In terms of the Gibbs free energy φ =φ(θ, E E E , T, q), inequality (17) can be given the form It is worth remarking that the application of (5) in connection withψ(θ, E E E , T, q, F) would result in inequalities different from (19) though the entropy production is invariant, namely

Hypo-Thermoelastic Solids
Assume that σ ET R = 0, σ Q R = 0 and Hence, Equations (16) and (18) become Equations (21) and (22) are said to describe hypo-thermoelastic solids; owing to the condition σ R = 0, such models are viewed as non-dissipative. More explicit relations foṙ T RR andq R follow by appealing to Lemma 1.
First, let N = ∂ T RR ψ/|∂ T RR ψ| and Z =ṪRR. The representation formula (14) is applicable in that Z · N is known, Accordingly, we select G = J J J RRĖ , where J J J RR is an arbitrary fourth-order tensor, and let J J J RR depend on θ, E, T RR , q R . It follows from (21) thaṫ This equation can be written in a compact form by letting We can then write (21) in the forṁ Likewise, we can also apply (14) with the vectors N = ∂ q R ψ/|∂ q R ψ| , Z =q R and G = J R ∇ R θ, where J R is is an arbitrary second-order tensor. Using (22), we have and thenq We can write (22) in the compact form, Relations (24) and (26) describe the constitutive properties of hypo-thermoelastic solids and C C C RR and K R are referred to as hypo-elastic and hypo-thermal tensors, respectively. By the arbitrariness of J J J RR and J R , it follows that there are infinitely many tensors C C C RR and K R compatible with a given free energy ψ. Moreover, C C C RR and K R need not be positive definite and this is a remarkable difference from the positive definiteness of the elastic tensor and the conductivity tensor.

Remark 2.
Relative to rigid heat conductors, if we assume ε = c v θ + ε 0 , c v > 0, the balance of internal energy takes the form, ρ R c vθ = −∇ R · q R + ρ R r and hence, it follows from (26) that This equation is the same as that resulting from Green and Naghdi's type II theory (see [23][24][25][26]). It is hyperbolic if K R is positive-definite, but this is not implied by thermodynamics.
Otherwise, if hypo-thermoelastic constitutive relations are given in advance, for in-stanceṪ then from (21) and (22), we obtain and the arbitrariness ofĖ and ∇ R θ implies In [27], the existence of a thermodynamic potential ψ satisfying (21) or (22) is investigated by exploiting the overdetermined systems (28) with assigned tensorsĈ C C RR andK R . Proposition 1. If we replace (28) into (23) and (25), we obtain Hence, by properly choosing the arbitrary terms (namely, if we let J J J RR =Ĉ RR and J R =K R ), we obtain the expected identities C RR =Ĉ RR and K R =K R .
A similar approach can be developed if the Gibbs free energy φ is given instead of ψ.
and assuming ∂ E φ = 0, it followṡ Hence, by paralleling previous arguments, we obtaiṅ where with X RR being an arbitrary (possibly vanishing) fourth-order tensor. Two further examples are now given by starting with the potential, φ or ψ.

Example 1.
Let φ be the Gibbs free energy defined by (8), and let with α and β being real-valued functions. Consequently, Let X RR = µI. From (31), we have the following: Apparently, H RR enjoys the major symmetry. Choosing µ = α, we obtain The resulting constitutive equation in the form (30) is generated by the well known Prandtl-Reuss work hardening plasticity theory provided that E and T RR are replaced by their deviatoric parts (see [4] Sect. C).

Example 2.
Let M be a non-singular, fully symmetric, fourth-order tensor; M be a non-singular, symmetric, second-order tensor; and G G G be a smooth function from Sym to Sym. Consider the free energy ψ defined by As a consequence, we infer that (28) holds witĥ Then, in light of Remark 3, we let J J J RR =Ĉ RR and J R =K R into (23) and (25) to obtaiṅ In the linear case, G G G(E) = LE, L being a (not necessarily symmetrical) fourth-order tensor, we have the following:Ṫ RR = C RRĖ , C RR := L + M −1 .

Hypo-Elastic Models with Thermal Dissipation
Unlike hypo-thermoelastic models, which exhibit no entropy production, Maxwell-Cattaneo and Fourier-like hypoelastic models allow for a non-negative value of σ R . A quite general class of hypo-elastic models with thermal dissipation is characterized by σ ET R = 0 but σ R ≡ σ Q R ≥ 0, so that (21) still holds along with (18). As in the previous section, the assumption ∂ T RR ψ = 0 leads to (24), whereas (18) can be satisfied in two different ways, depending on whether it is ∂ q R ψ = 0 or ∂ q R ψ = 0.

Fourier-Like Models
Assuming ∂ q R ψ = 0, it follows that ψ = ψ(θ, E, T RR ) and (18) reduces to the wellknown Fourier inequality This condition is satisfied, for example, by letting q R be given by Fourier's law, q R = −κ∇ R θ, where κ is a positive-definite second-order tensor. In that case, we have The same result follows also by using the representation formula (14). Let A be any symmetric, non-singular, second-order tensor, possibly parameterized by the temperature θ. Assuming A∇ R θ = 0 and applying (14) with N = A∇ R θ/|A∇ R θ|, Z = A −1 q R , and G = 0, it follows from (34) that Hence, letting κ = A 2 and σ Q R = ∇ R θ · κ∇ R θ/θ 2 , we recover Fourier's law q R = −κ∇ R θ.

Maxwell-Cattaneo-like Models
Assume ∂ q R ψ = 0. Equation (18) can then be written in the form Applying (14) with where K R = K R (θ, E, T RR , q R ) is given by (25) and For definiteness, we now show that a class of models for heat conduction, of the Maxwell-Cattaneo type [28], follows from (35). Let ψ depend on q R via ξ = |q R | n , n ≥ 2. Hence, ∂ q R ψ = n ∂ ξ ψ |q R | n−2 q R , ∂ ξ ψ = 0 and inequality (18) becomes Let σ Q R = |q R | 2 /θκ, where κ is a positive-valued function. In view of (35), we have the following: Consequently, κnρ R θ∂ ξ ψ|q R | n−2q R + q R = −κ∇ R θ Can be viewed as a Maxwell-Cattaneo equation with Playing the role of relaxation time and κ representing the heat conductivity. If n = 2, then τ = 2κρ R θ∂ ξ ψ.
In this case, τ reduces to a function of the temperature alone when κ = κ(θ) and ψ is a quadratic function of q R .

Thermo-Viscoelastic and Viscoplastic Models
Models of thermo-viscoelastic materials are characterized by non-negative entropy productions σ ET R and σ Q R independent ofĖ,Ṫ RR ,q R and ∇ R θ. In general, both the free energy and the entropy productions are functions of (θ, E, T RR , q R ). We recognize two different classes of models.
Assuming ∂ T RR ψ, ∂ q R ψ = 0 and applying the representation formula (14), we can rewrite (16) and (18) in the forṁ G RR and g R being respectively an arbitrary tensor and an arbitrary vector.
In terms of the free enthalpy φ, Equation (37) 1 takes the forṁ where Incremental relations of this type describe a thermo-viscoelastic behaviour when σ ET R is smooth. If instead σ ET R is only piecewise smooth and vanishes in a suitable open region then a rate-dependent thermo-viscoplastic behaviour is concerned; for the one-dimensional isothermal case, see for instance [8,29,30].
Otherwise, when ∂ T RR ψ = ∂ q R ψ = 0 relations (16) and (18) reduce to The latter equation has been scrutinized in Section 5.1. Concerning T RR , let M be any symmetric, non-singular fourth-order tensor, possibly parameterized by θ. Assuming MĖ = 0 and applying (14) with N = MĖ/|MĖ|, Z = M −1 [T RR − ρ R ∂ E ψ] and G = 0, we can express T RR in the form We now describe some simple examples of (37)-(39) via special choices of the free energy and the non-negative functions σ ET R and σ Q R .

Thermo-Viscoelastic Behaviour
Let M be a fully symmetric, positive-definite, fourth-order tensor; M be a symmetric, positive-definite, second-order tensor; G G G be a smooth function from Sym to Sym; and ψ be given by (32) as in Example 2. Then, we have and we can apply (37), where Since [X ⊗ Y]Z = (Y · Z)X for tensors or vectors X, Y, Z, we have from (37), we have the following: simply by choosing G RR and g R equal to the right-hand members of the equalities in (40). When G G G is linear, namely G G G(E) = G ∞ E, from (40), we recover the standard linear solid (or Zener) model with Mawell-Cattaneo heat conduction, where G 0 = G ∞ + M −1 and κ = M −1 /θ. Here, G 0 and G ∞ stand for the usual elastic and relaxation moduli, respectively, whereas κ stands for the conductivity tensor. The corresponding free energy is For solids, G ∞ > 0. In the special case G ∞ = 0, we obtain the Maxwell model for fluids.Ṫ It is worth looking for the possible relation between (41) 1 (1). Under the small strain assumption, we observe that

and the classical linear equation in Equation
The approximation being motivated by linearity. Hence, we might replace (1) with As an example, let Hence, upon time differentiation, we finḋ whence (41) 1 follows.

The Bingham-Norton Model
This model describes the behaviour of an elastic-perfectly viscoplastic solid. In our setting, it can be derived by assuming ∂ E ψ = ∂ q R ψ = 0. Accordingly, ρ R ∂ E φ = −T RR , and then from (38), we obtainĖ where (for simplicity, we let G RR = 0) Choosing λ, S y > 0 and letting We recover the Bingham-Norton model [31] for n > 1. If, instead, n = 1 and S y = 0, the model reduces toĖ which may be viewed as a generalized Maxwell-Wiechert model. In any case, the free energy is a generic function ψ = ψ(θ, T RR ).
The Bingham-Norton model (42) and (43) allows for a special decomposition of the deformation rate,Ė =Ė e +Ė vp . The elastic strain rate is represented bẏ whereas the rate of visco-plastic strain is only a function of the stress and depends on the initial yield stress S y (there is no influence of hardening) This happens by letting where H denotes the Heaviside step function.

The Kelvin-Voigt Model
Then, (39) becomes When large deformations are involved, this is the three-dimensional continuummechanics version of the well-known rheological model realized by a spring and a dashpot in parallel, where the stress is the sum of terms proportional to the strain and the strain rate. Under the small strain assumption, the linear constitutive Equation (44) is approximated by replacing E andĖ with ε andε, respectively.

Thermo-Viscoelastic and Viscoplastic Models in the Spatial Description
The difficulty in formulating incremental, or rate-type, models in the spatial description of finite deformations is related to the appropriate selection of objective time derivatives. This difficulty is overcome here by translating the visco-thermoelastic models previously formulated by means of (19) in the spatial description.
Taking into account that E E E = D, Equation (19) are written in the form where both the free energy and the non-negative entropy productions are functions of (θ, E E E , T, q). Assuming ∂ T ψ, ∂ q ψ = 0 are nonzero and applying the same arguments of the previous sections, we can rewrite (45) in the form where with G and g being, respectively, an arbitrary tensor and an arbitrary vector.
In terms of the free enthalpy φ, Equation (38) takes the form where G being an arbitrary second-order tensor.

Thermo-Viscoelastic Behaviour
Let A be a fully symmetric, positive-definite, fourth-order tensor and A be a symmetric, positive-definite, second-order tensor. Optionally, both A and A are parameterized by θ. Moreover, let F F F be a monotone function from Sym to Sym. Given Both are non-negative, and from (45), we have the following: Now, letting ψ such that We obtain These equations are certainly satisfied if where

Kelvin-Voigt Model in the Spatial Description
When ∂ T ψ, ∂ q ψ = 0, relations (46) reduce to The latter equation can be handled as in Section 5.1. Applying the procedure exhibited in Section 6.3, the former relation yields where A is a symmetric, non-singular fourth-order tensor, possibly parameterized by θ. Let Then, (51) becomes In the isotropic case, F 0 = λ 0 1 ⊗ 1 + 2µ 0 I, A 2 = λ1 ⊗ 1 + 2µI, we obtain

Bingham-Norton Visco-Plastic Fluid
Assuming ∂ E E E ψ = ∂ q ψ = 0 and X = 0 from (47), we have the following: Choosing λ, τ y > 0, n ∈ N ∪ {0} and We recover the Bingham-Norton model in the spatial description. The free energy is given by a generic function ψ = ψ(θ, T). In particular, when We have the following: so that H enjoys both the minor and major symmetries. The special case α ≡ 0 leads to the well-known model of a Bingham plastic fluid: is referred to as Norton-Hoff nonlinear viscosity; it reduces to a constant value, µ = λ, when n = 1.

Higher-Order Rate Models
Non-Newtonian phenomena exhibited, e.g., by asphalt and some biomaterials show properties that can be associated to different relaxation times or rather to higher-order differential equations. This view is often described by having recourse to the Burgers model [32] or the best known Oldroyd-B fluid model [33].

Burgers Material
While sometimes the model is motivated by rheological analogs [34] and framed within a scheme with intermediate reference configurations [35], here, the model is investigated within the thermodynamic approach developed so far. Let The stress S in the Burgers model describes an incompressible fluid and is given by the second-order rate equation We expect that, for a compressible fluid, the same equation should hold with the Oldroyd derivative S replaced with the Truesdell derivative S. The thermodynamic consistency of (52) can be proven by having recourse to the counterpart of (52) in the reference configuration. Let S RR = JF −1 SF −T be the corresponding stress in the reference configuration. We let be the analog of the Burgers model in the referential version. To investigate the thermodynamic consistency of (53), we assume the pressure p, the entropy η, and the free energy ψ as functions of the set of variables Ξ R = (θ, E, S RR ,θ,Ė,Ṡ RR ,θ,Ë,S RR ).
The Clausius-Duhem inequality becomes The linearity and arbitrariness of ... θ , Assume that p and η are independent ofθ andθ. Hence, the linearity and arbitrariness ofθ, and next ofθ, imply Consequently, the inequality simplifies to Further restrictions follow depending on the degrees of arbitrariness. The compressibility of the fluid is described by letting ψ depend on E through J = [det(1 + 2E)] 1/2 . Since J = ρ R /ρ, the classical relation p = ρ 2 ∂ ρ ψ becomes Consequently, (54) becomes We now investigate the consistency of the Burgers-like evolution Equation (53) with the reduced inequality (55), where the free energy ψ has the form ψ = ψ(θ, J, S RR ,Ė,Ṡ RR ).
For definiteness, we let β = 0 and considerS RR as a function of S RR ,Ṡ RR ,Ė,Ë. Upon substitution ofS RR from (53) into (55), we have The linearity and arbitrariness ofË imply and then It is worth contrasting the restrictions obtained with what arises from an application of the Lemma 1 about a second-order tensor Z, with G being an arbitrary second-order tensor. Here, we observe that, in view of (55), equationS is consistent with thermodynamics if it holds for any pair of functions ψ, σ of θ,Ė, S RR ,Ṡ RR subject to σ ≥ 0. The selection Z =S RR and N = ∂Ṡ RR ψ/|∂Ṡ RR ψ| into (60) together with the exploitation of (61) providë For definiteness we consider the case ν = 0, and then we let Accordingly, so that (62) becomes Letting βG = 2µĖ − λṠ RR , we obtain exactly Equation (53) with ν = 0. It is of interest to determine the Eulerian version of the rate-type Equation (53). We know thatṠ To evaluateS RR , we observe thaẗ Moreover, sinceĖ = F T DF, then with D being the Cotter-Rivlin derivative. SubstitutingṠ RR ,S RR ,Ė,Ë in (53), we obtain where B = FF T is the Finger strain tensor. Alternatively, in light of (13), If the fluid is incompressible, then The power of the pressure vanishes. Consistently, we let ∂ E ψ = 0. The results (56) and (57) hold unchanged while (63) changes to Hence, in the linear approximation (J 1, B 1, D D • D), the Burgers model (52) is recovered.

Oldroyd-B Model
This model, often used to describe the flow of viscoelastic fluids, can be regarded as an extension of the upper-convected Maxwell model and represents an incompressible fluid filled with elastic bead and spring dumbbells. In this model, named after J. G. Oldroyd [33], the stress is given by the second-order rate equation Letting λ = λ p , µ = µ s + µ p , ν = 2µ s λ p , we can split (64) into a polymeric (viscoelastic) part separately from the solvent (viscous) part, where µ s and µ p denote the viscosity coefficients of the solvent and of the solute, respectively. To prove the thermodynamic consistency of (64), we start from its counterpart in the reference configuration: where T RR = JF −1 TF −T , and assume the entropy η and the free energy ψ as functions of the set of variables Ξ R = (θ, E, T RR ,θ,Ė,Ë).
To determine the Eulerian version of the rate-type Equation (65), we follow the same procedure as in the previous subsection, and taking into account the incompressibility condition, we obtain In the linear approximation (J 1, B 1, D D • D), Equations (64) and (68) coincide. This proves the thermodynamic consistency of Oldroyd-B model provided only that λ = 0, µ > 0.

Conclusions
Thermodynamically consistent viscoelastic models are established on the basis of the following points. The constitutive equations are essentially represented in the Lagrangian description in terms of the Piola-Kirchhoff stress T RR and the Green-St.Venant strain E. For first-order rate type models, the second-law inequality is then stated in the form The entropy production σ R = Jσ is taken as a constitutive function of the chosen independent variables; a similar treatment is developed for the heat conduction effects. The corresponding second-law inequality in the Eulerian description is found to be (ρ∂ E E Eψ − T) · D + ρ∂ Tψ · T= −θσ ≤ 0, with E E E being the Eulerian Almansi tensor. This formulation proves profitable in the analysis of the thermodynamic consistency of rate-type equations in the Eulerian description. A more complex second law expression than (69) is involved for higher-order rate type models such as Burgers and Oldroyd-B (see § 8).
We now comment on the main novelty of our approach. IfĖ andṪ RR are assumed to be independent, then one finds that If insteadĖ andṪ RR are not independent, then we determine the mathematical consequence via Lemma 1. For instance, the selection N = ∂ T RR ψ/|∂ T RR ψ| yields a class of rate-type constitutive equations for the stress in the forṁ where G is an arbitrary second-order tensor. Equation (70) shows that first-order rate equations, consistent with thermodynamics, are determined by a pair of the free energy ψ and the entropy production σ. As a familiar example, the Maxwell model (2), follows by letting where G 0 is a fully symmetric, positive-definite, fourth-order tensor. Indeed, after replacing σ and ∂ T RR ψ in (70), we obtaiṅ The first and last fractional terms cancel each other out, and we haveṪ Throughout this paper, by a proper selection of the independent variables, we determined constitutive equations in the rate-type form, possibly of higher order than the first. To our mind, this has the technical advantage that the recourse to histories is avoided. Moreover, the role of the entropy production as a constitutive quantity allows for the modelling of further dissipation effects as is the case in hysteretic phenomena [8].
Following the lines developed in this paper, future work will be devoted to phase transitions and hysteresis in shape memory alloys. Data Availability Statement: The study did not report any data.