Next Article in Journal
Nano Minerals Applications in Pollutants Removal: A New Open Special Issue in Materials
Next Article in Special Issue
Viscoelastic and Electromagnetic Materials with Nonlinear Memory
Previous Article in Journal
Inverse Design of Valley-Like Edge States of Sound Degenerated Away from the High-Symmetry Points in a Square Lattice
Previous Article in Special Issue
Studies on the Mechanical Models and Behaviors for the Stamp/Film Interface in Microtransfer Printing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Magneto-Viscoelastic Materials: Memory Functionals and Rate Equations

1
Dipartimento di Ingegneria Civile, Architettura, Territorio, Ambiente e di Matematica, Università di Brescia, Via Valotti 9, 25133 Brescia, Italy
2
Dipartimento di Informatica, Bioingegneria, Robotica e Ingegneria dei Sistemi, Università di Genova, Via All’Opera Pia 13, 16145 Genova, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Materials 2022, 15(19), 6699; https://doi.org/10.3390/ma15196699
Submission received: 23 August 2022 / Revised: 22 September 2022 / Accepted: 23 September 2022 / Published: 27 September 2022
(This article belongs to the Special Issue Modelling of Viscoelastic Materials and Mechanical Behavior)

Abstract

:
The properties of viscoelastic solids subject to a magnetic field are modelled within two thermodynamically consistent approaches that are typical of models with a non-instantaneous response. One is based on memory functionals: the reversible changes are described by the instantaneous response, while the dissipativity is expressed by the dependence on histories. The other approach involves objective rate equations. While memory functionals lead to the difficulty of determining thermodynamically consistent free energy functionals, rate equations result in a simpler scheme. The greater simplicity allows the discovery of, in particular, models of magneto-hyperelastic materials, magneto-hypoelastic materials, and various forms of magneto-viscoelastic behaviour. The novelty of the procedure is based on two features: a representation formula, originating from the entropy inequality, and the use of the entropy production as a constitutive function. Relations with other approaches in the literature are examined in detail.

1. Introduction

The response of magnetic materials to an applied magnetic field is generally not instantaneous. In addition to showing a partial instantaneous response, the material gradually approaches equilibrium in a finite time dependent on the deformation. Magneto-viscoelastic materials are intended to be models that exhibit both instantaneous changes of magneto-mechanical properties and a variable time response when acted upon by a magnetic field. This subject is of interest for applications and requires both appropriate balance equations and constitutive equations. Balance equations are treated, e.g., in [1,2], while interesting constitutive equations are developed in, e.g., [3,4]. Updated lists of references are given in [5,6,7]. Despite the various approaches and procedures developed in the literature, the subject deserves further attention, hopefully to create simpler models.
Recently, we have developed a systematic approach to nonlinear materials with memory [8,9]. The crucial points are that the models are thermodynamically consistent; the entropy production is a constitutive quantity that characterizes the dissipative properties; the entropy inequality allows us eventually to derive a representation formula for the pertinent function; and the nonlinear rate equations are derived with an objective time derivative.
In essence, magnetic viscoelasticity involves the interaction between mechanical and magnetic fields within dissipative processes in materials with a non-instantaneous response. This scheme is realized by letting the independent variables occur through their histories or by considering rate-type constitutive equations. Further, in both cases, the thermal properties are modelled, and the restrictions placed by objectivity and thermodynamic consistency are investigated in detail in order to be able to make a comparison between the characteristics of the two approaches.
While balance equations are used as the standard in the literature (see, e.g., [6,10,11]), constitutive equations are placed in quite new settings. The modelling through memory functionals involves constitutive equations with a joint dependence on the present values and thermal, deformation, and magnetic field histories so that, at time t, the response of the material is determined by the present values of the temperature θ , the deformation gradient F , the magnetic field H , and the temperature gradient θ ,
θ ( t ) , F ( t ) , H ( t ) , θ ( t )
as well as the histories, up to time t,
F t , H t , θ t .
Indeed, the scalar functions (internal energy, entropy, free energy) are required to be in-variant under Euclidean transformations, and this implies that the dependence on F , H , and θ occurs through their invariants. As a result of the thermodynamic requirements placed by the second law inequality, explicit representations of the stress tensor, magnetization, and heat flux are established.
Next, we elaborate on differential models expressed by rate-type equations, and, hence, objective derivatives are required in the rate equations. Due to the occurrence of the time derivative, objectivity indicates that the Lagrangian description is more convenient as the starting step. A general scheme is allowed by letting an invariant stress T , the Green–Saint Venant tensor E , and a magnetic field vector be the independent variables. The generality is allowed by the joint dependence on T and E and by letting the entropy production be a constitutive function, which is quite new in the literature. Furthermore, to determine the restrictions placed by the second law, we apply a representation formula for tensors and vectors given by Equation (4). Depending on the degree of arbitrariness of the independent (magnetic) variable, the thermodynamic analysis proves that hyper-magnetoelasticity, hypo-magnetoelasticity, and dissipative magneto-viscoelasticity are fully characterized.
The relations with other approaches in the literature are given in Section 6. Here, we indicate two features of the present work. First, the modelling through memory functionals describes the instantaneous and non-instantaneous response to the deformation, magnetic field, and temperature gradient. The thermodynamic consistency is made formal in (8) and (9). Next, specific (linear) representations are established by taking the stress, magnetization, and heat flux in the forms of (10), (11), and (14). Secondly, rate equations are considered for the stress and magnetization in terms of the deformation and magnetic field. Through the representation formula, we find the results (28) for the stress T . All of the rates given by (28) are consistent with thermodynamics; for any free energy and any entropy production, we find a physically admissible rate equation. In particular, Equations (24) and (27) describe magneto-hyperelastic materials and magneto-hypoelastic materials, respectively.
Notation. 
We consider a body occupying a time-dependent region Ω E 3 . The motion is described by means of the function χ ( X , t ) , providing the position vector x Ω = χ ( 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 i K = X K χ i . The invertibility of X x = χ ( X , t ) is guaranteed by letting J : = det F > 0 . For any tensor A , we define | A | as ( A · A ) 1 / 2 .

2. Balance Equations

Let v ( x , t ) be the velocity field on Ω × R . A superposed dot denotes the time differentiation following the motion of the body, and hence, for any function f ( x , t ) , we have f ˙ = t f + v · f . We denote by L the velocity gradient, L i j = x j v i , and recall that
F ˙ = L F .
The right Cauchy–Green tensor C and the Green–Saint Venant deformation tensor E are defined by
C : = F T F , E : = 1 2 ( C 1 ) ,
where 1 is the second-order identity tensor. Moreover, D denotes the stretching tensor, D = sym L , and W the spin tensor, W = skw L .
Let ε be the internal energy density (per unit mass), T the symmetric Cauchy stress, q the heat flux vector, ρ the mass density, r the external heat supply, and b the mechanical body force per unit mass. Let m = M / ρ be the magnetization per unit mass and H the magnetic intensity. The balance equations for mass and linear momentum are written in the form (see, e.g., [11,12,13])
ρ ˙ + ρ · v = 0 , ρ v ˙ = · T + ρ b + f M ,
where f M is the force per unit volume of magnetic character,
f M = μ 0 ρ ( m · ) H μ 0 ϵ 0 ρ [ m ˙ × E + v × ( m · ) E ] ;
in stationary conditions, we can take E = 0 . Additionally, for later convenience, we let H be the magnetic intensity at the frame locally at rest with the body. The balance of angular momentum and energy can be written in the form
skw ( T + μ 0 ρ H m ) = 0 ,
ρ ε ˙ = μ 0 ρ H · m ˙ + T · L · q + ρ r
Let η be the entropy density and θ the absolute temperature. As for the statement of the second law of thermodynamics, we let the inequality
ρ η ˙ + · q θ ρ r θ = σ 0
hold for any process compatible with the balance equations. The scalar σ , or the entropy production, is non-negative and is viewed as a constitutive function. Hence, the thermodynamic process consists of η , q , r , σ , and the other functions occurring in the balance equations.
In terms of the magnetic Gibbs free energy
ϕ = ε θ η μ 0 m · H
the entropy inequality can be written in the form
ρ ( ϕ ˙ + η θ ˙ ) μ 0 M · H ˙ + T · L 1 θ q · θ = θ σ 0 .
Based on (2), next, we describe the magneto-viscoelasticity by the memory functionals or rate equations and examine the thermodynamic consistency.

Representation Formula

The thermodynamic analysis usually leads to relations of the form
Z · K + A · F = f ,
where, to fix ideas, Z , K , A , F are second-order tensors and f is a scalar. If K and F are arbitrary and independent, then it follows that Z = 0 , A = 0 (and f = 0 ). If, instead, K and F are not independent, then we can determine the relation between Z and A through a representation formula [8,9].
Let N be a unit tensor, | N | = 1 . Then
Z = ( Z · N ) N + Z , Z · N = 0 .
If it happens that Z · N is known, say Z · N = g , whereas Z is unknown, then Z can be expressed by
Z = ( I N N ) G = G ( G · N ) N ,
where I is the fourth-order unit tensor, and G is an arbitrary second-order tensor. Once Z · N = g is given, we can write the representation formula
Z = g N + ( I N N ) G .
Returning to (3), we let N = K / | K | to yield
Z = f A · F | K | 2 K + ( I K | K | K | K ) G .
A strictly analogous relation holds if the tensors are replaced with vectors.

3. Constitutive Assumptions

The constitutive assumptions are suggested for several purposes. First, we allow for interaction between the deformation and temperature fields with magnetization; this indicates that θ , F , H , θ are among the independent variables. Moreover, the time delay in response motivates the dependence on the histories F t , H t , θ t . Hence, we let
Γ = ( θ , F , H , θ , F t , H t , θ t )
be the set of independent variables.
The internal energy (density) ε , the entropy η , and the Gibbs free energy ϕ are invariant under a change of frame. The constitutive equations for ε , η , and ϕ are then required to provide invariant values. Now, both F and H are not invariant. Under a change of frame
x * = c + Q x , Q Q T = 1 ,
F and H change as vectors,
F * = Q F , H * = Q H .
Instead,
E = 1 2 ( F T F 1 ) , H = F T H
are invariant in that
E * = F * T F * 1 = F T Q T Q F 1 = F T F 1 = E , H * = F * T H * = F T Q T Q H = F T H = H .
Of course, H · H is also invariant,
H * · H * = Q H · Q H = H · Q T Q H = H · H .
Incidentally,
H · H = F T H · F T H = H · ( F F T ) H .
Hence, the dependence on the pair H , E does not include that on F · F , H · H . For definiteness, we let Λ be the pair of functions
f ( F ) = 1 2 F · F , h ( H ) = 1 2 H · H .
As for the dependence on θ , we observe that
x * θ = Q x θ or * θ = Q θ .
Hence, θ is a vector while the referential gradient
R θ = θ F ,
is invariant.
Let θ , F , H , and θ be differentiable as t R for any x Ω . Define the constant continuation of θ t 0 , F t 0 , and H t 0 as
θ ˜ t ( u ) = { θ ( t 0 ) , u [ 0 , t t 0 ) , θ t 0 ( u ( t t 0 ) ) , u [ t t 0 , ) , θ ˜ ( t ) = θ ˜ t ( 0 ) ,
and the like for F and H while θ = 0 on [ t 0 , t ] .
While Γ is the set of independent variables, objectivity and modelling purposes indicate that we let
ϕ = ϕ ( θ , E , H , R θ , Λ , E t , H t , R θ t , Λ t ) .
With a small abuse of notation, the dependence on time is denoted with the same symbol
ϕ ( t ) = ϕ ( Γ ( t ) ) .
We first show that the free energy functional ϕ ( Γ ) is required to satisfy a minimum property. If θ , F , H , and, hence, Λ , are constant, and θ = 0 in the interval [ t 0 , t ] , then (2) implies that
ϕ ˙ ( τ ) 0 , τ ( t 0 , t ) ,
and hence,
ϕ ( θ ˜ ( t ) , E ˜ ( t ) , H ˜ ( t ) , Λ ˜ ( t ) , 0 , E ˜ t , H ˜ t , R θ ˜ t , Λ ˜ t ) ϕ ( θ 0 , E 0 , H 0 , Λ 0 , 0 , E t 0 , H t 0 , R θ t 0 , Λ t 0 ) ,
where we let θ 0 = θ ( t 0 ) , E 0 = E ( t 0 ) , H 0 = H ( t 0 ) ,   and   Λ 0 = Λ ( t 0 ) . By the continuity of the functional, as t t 0 we have
ϕ ( θ 0 , E 0 , H 0 , 0 , Λ 0 , E ˜ t , H ˜ t , R θ ˜ t , Λ ˜ t ) ϕ ( θ 0 , E 0 , H 0 , 0 , Λ 0 , E 0 , H 0 , 0 , Λ 0 ) ,
where E 0 , H 0 , 0 , and Λ 0 are the constant histories with values of E 0 , H 0 , 0 , and Λ 0 . Thus, it follows that
ϕ ( θ 0 , E 0 , H 0 , 0 , Λ 0 , E 0 , H 0 , 0 , Λ 0 ) ϕ ( θ 0 , E 0 , H 0 , 0 , Λ 0 , E 0 t 0 , H 0 t 0 , R θ t 0 , Λ t 0 ) .
This is the content of the minimum property: among all histories E t , H t , R θ t , and Λ t with the given present values of E 0 , H 0 , 0 , and Λ 0 , none yields a smaller value of the free energy than that corresponding to the constant histories E 0 , H 0 , 0 and Λ 0 .

4. Thermodynamic Restrictions

We now compute the time derivative ϕ ˙ and substitute in the entropy inequality to obtain
ρ ( θ ϕ + η ) θ ˙ ρ E ϕ · E ˙ ρ H ϕ · H ˙ ρ θ ϕ · θ ¯ ˙ μ 0 M · H ˙ ρ f ϕ F · F ˙ ρ h ϕ H · H ˙ + T · L 1 θ q · θ ρ d ϕ ( Γ | E ˙ t ) ρ d ϕ ( Γ | H ˙ t ) ρ d ϕ ( Γ | θ ¯ ˙ t ) ρ d ϕ ( Γ t | f ˙ t ) ρ d ϕ ( Γ t | h ˙ t ) 0 .
Here, we have made the dependence on Λ = ( f , h ) and Λ t = ( f t , h t ) explicit. The linearity and arbitrariness of θ ¯ ˙ , θ ˙ imply
θ ϕ = 0 , η = θ ϕ .
Observe that
E ˙ = F T D F , H ˙ = F T L T H + F T H ˙ , F · F ˙ = ( F F T ) · L .
Since L = D + W , the occurrence of W in the inequality is through
( T + ρ F H ϕ H ρ f ϕ F F T ) · W
and this quantity has to vanish. Now F F T Sym and then ( F F T ) · W = 0 identically. Hence, we have
( T + ρ F H ϕ H ) · W = 0
which implies
T + ρ F H ϕ H Sym .
We now can write the remaining inequality in the form
( T ρ F E ϕ F T + ρ F H ϕ H ρ f ϕ F F T ) · D ( μ 0 M + ρ F H ϕ ρ h ϕ H ) · H ˙ 1 θ q · θ ρ d ϕ ( Γ | E ˙ t ) ρ d ϕ ( Γ | H ˙ t ) ρ d ϕ ( Γ | θ ¯ ˙ t ) ρ d ϕ ( Γ t | f ˙ t ) ρ d ϕ ( Γ t | h ˙ t ) 0 .
The linearity and arbitrariness of D , H ˙ imply
T = ρ F E ϕ F T ρ F H ϕ H , μ 0 M = ρ F H ϕ + ρ h ϕ H .
Inequality (7) then reduces to
1 θ q · θ ρ d ϕ ( Γ | E ˙ t ) ρ d ϕ ( Γ | H ˙ t ) ρ d ϕ ( Γ | θ ¯ ˙ t ) 0 .
According to (6) and (8), it follows that
skw T = skw ( μ 0 M H )
which is just the requirement placed by the balance of angular momentum. It is worth remarking that this requirement holds merely because of the dependence of ϕ on H through H = F T H . Instead, the dependence of ϕ on h leaves skw T unchanged.
Further restrictions, placed by the reduced inequality (9), follow by considering some particular cases. First, we assume the temperature is uniform at any time, θ t = 0 . Moreover, nonlinear constitutive equations for T and M are established by selecting the partial derivatives
E ϕ = G 0 ( θ , | E | ) E + g ( θ ) 0 G ( u ) E ( t u ) d u , G 0 , G Lin ,
H ϕ = A 0 ( θ , | H | ) H a ( θ ) 0 A ( u ) H ( t u ) d u , A 0 , A Lin .
Now, given G ( θ , | E | ) , a fully symmetric fourth-order tensor-valued function, we have
E 1 2 E · G E = G E + 1 2 | E | E · ( | E | G ) E E
and the like for a symmetric second-order tensor-valued function A ( θ , | H | ) . Hence, letting 1 and I be the second- and fourth-order unit tensors, we define
G 0 = G + 1 2 | E | E · ( | E | G ) E I , A 0 = A + 1 2 | H | H · ( | H | A ) H 1 ,
so that
G 0 E = E 1 2 E · G E , A 0 H = H 1 2 H · A H .
Hence, the sought functional ϕ ( θ , E , H , E t , H t ) takes the form
ϕ = 1 2 E · G E + g ( θ ) E · 0 G ( u ) E ( t u ) d u 1 2 H · A H a ( θ ) H · 0 A ( u ) H ( t u ) d u + Φ ( θ , E t , H t ) ,
and the functional Φ is thus far undetermined. According to (9), the functional ϕ has to satisfy the inequality
d ϕ ( Γ | E ˙ t ) + d ϕ ( Γ | H ˙ t ) 0
and this is eventually the requirement of the unknown functional Φ .
For definiteness, we consider Φ in the form
Φ ( θ , E t , H t ) = ϕ 0 ( θ ) + 1 2 g ( θ ) 0 E t ( u ) · G ( u ) E t ( u ) d u 1 2 a ( θ ) 0 H t ( u ) · A ( u ) H t ( u ) d u
and, hence, the functional ϕ can be written as ϕ = ϕ m ,
ϕ m : = ϕ 0 ( θ ) + 1 2 E · G E 1 2 g ( θ ) 0 [ E ( t ) E ( t u ) ] · G ( u ) [ E ( t ) E ( t u ) ] d u 1 2 H · A H + 1 2 g ( θ ) 0 [ H ( t ) H ( t u ) ] · A ( u ) [ H ( t ) H ( t u ) ] d u .
It follows that
d ϕ ( Γ | E ˙ t ) = 1 2 g ( θ ) 0 [ E ( t ) E ( t u ) ] · G ( u ) E ˙ ( t u ) ] d u = 1 2 g ( θ ) { ( E ( t ) E ( t u ) ) · G ( u ) ( E ( t ) E ( t u ) ) 0 0 ( E ( t ) E ( t u ) ) · G ( u ) ( E ( t ) E ( t u ) ) d u } .
The assumption G ( ) = 0 implies
d ϕ ( Γ | E ˙ t ) = 1 2 g ( θ ) 0 [ E ( t ) E ( t u ) ] · G ( u ) [ E ( t ) E ( t u ) ] d u .
Likewise, letting A ( ) = 0 , we find
d ϕ ( Γ | H ˙ t ) = 1 2 a ( θ ) 0 [ H ( t ) H ( t u ) ] · A ( u ) [ H ( t ) H ( t u ) ] d u .
Hence, (12) holds if and only if
g ( θ ) G ( u ) 0 , a ( θ ) A ( u ) 0
for all u [ 0 , ) .
The minimum property of ϕ m at the constant histories E , H holds if and only if
g ( θ ) G ( u ) 0 , a ( θ ) A ( u ) 0
for all u [ 0 , ) .
Hence, the functional ϕ m is thermodynamically consistent if G and G , as well as A and A , have opposite types of definiteness.

Heat Conduction

Now we let R θ t 0 , and, for simplicity, we look for models where R θ is independent of F and H so that the reduced dissipation inequality splits into (12) and
ρ θ d ϕ ( Γ | R θ ˙ t ) + 1 θ q · θ 0 .
Multiply this inequality by J = det F and observe that, using the referential heat flux q R = J F T q , we can write
ρ R θ d ϕ ( Γ | R θ ˙ t ) + q R · R θ 0 .
For definiteness, let q R be given by the constitutive functional
q R ( θ ( t ) , R θ t ) = α ( θ ) K 0 0 β ( u ) R θ ( t u ) d u ,
where K 0 is a positive-definite second-order tensor, while α and β are so far undetermined; we only assume β ( ) = 0 .
We let ϕ = ϕ m + ϕ c with ϕ c taken in the form
ϕ c = α ( θ ) | K 0 1 / 2 0 β ( u ) R θ ( t u ) d u | 2 .
The minimum property of ϕ , and hence of ϕ c , at R θ 0 holds if and only if α > 0 . Now, inequality (13) results in
ρ R θ K 0 0 β ( u ) R θ ( t u ) d u · 0 β ( u ) R θ ˙ ( t u ) d u + K 0 0 β ( u ) R θ ( t u ) d u · R θ 0 .
An integration by parts yields
0 β ( u ) R θ ˙ ( t u ) d u = 0 β ( u ) u R θ ( t u ) d u = β ( u ) R θ ( t u ) 0 + 0 β ( u ) R θ ( t u ) d u = β ( 0 ) R θ ( t ) + 0 β ( u ) R θ ( t u ) d u .
Hence, inequality (15) can be written in the form
( ρ R θ β ( 0 ) + 1 ) K 0 0 β ( u ) R θ ( t u ) d u · R θ ( t ) + ρ R θ K 0 0 β ( u ) R θ ( t u ) d u · 0 β ( u ) R θ ( t u ) d u 0 .
The linearity and arbitrariness of R θ ( t ) imply that
β ( 0 ) = 1 ρ R θ , β ( u ) β ( u ) 0 .
Since β ( 0 ) < 0 and ( β 2 ) 0 , then β ( u ) 0 for all u [ 0 , ) . Hence,
α ( θ ) K 0 β ( u ) 0
and, consistently, the kernel of the functional q R is negative definite.
Observe that
β ^ ( u ) = β ( u ) β ( 0 ) 0 .
Hence, we can write q R in the form
q R = α ( θ ) ρ R θ K 0 0 β ^ ( u ) R θ ( t u ) d u .
At the limit of short memory, we have
q R = α ( θ ) ρ R θ K 0 0 β ^ ( u ) d u R θ .
In the spatial description, the constitutive equation reads
q = α ( θ ) J ρ R θ F K 0 0 β ( u ) ( θ F ) ( t u ) d u .

5. Rate Equations in the Eulerian Description

It is a crucial point of magnetoelasticity, as well as of magneto-viscoelasticity, that the stress tensor need not be symmetric. Hence, the mechanical power T · L need not equal T · D , and, moreover,
T · L = J 1 ( F T R R F T ) · L = J 1 T R R · ( F T D F ) + J 1 T R R · ( F T W F ) .
Since E ˙ = F T D F , then
T · L = J 1 T R R · E ˙ + J 1 T R R · ( F T W F ) .
We start with the Eulerian description and write the Clausius–Duhem inequality in the form
ρ ( ϕ ˙ + η θ ˙ ) + T · L μ 0 M · H ˙ 1 θ q · θ = θ σ 0 .
Since we look for rate equations, objectivity indicates that the independent variables are invariant so that their time derivatives are invariant too. Hence, we assume that
ϕ = ϕ ( θ , T , E , H ) ,
where T is a stress-like variable to be identified. Consequently, the Clausius–Duhem inequality (17) can be written as
ρ ( θ ϕ + η ) θ ˙ ρ T ϕ · T ˙ ρ E ϕ · E ˙ ρ H ϕ · H ˙ + T · L μ 0 M · H ˙ 0 .
Since H ˙ = F T L T H + F T H ˙ , then we have
ρ ( θ ϕ + η ) θ ˙ ρ T ϕ · T ˙ ρ E ϕ · E ˙ + ( T ρ F E ϕ F T ρ H F H ϕ ) · ( F T E ˙ F 1 ) + ( T ρ H F H ϕ ) · W ( ρ F H ϕ + μ 0 M ) · H ˙ 0 .
Hence, it follows
η = θ ϕ ,
μ 0 M = ρ F H ϕ ,
skw ( T + μ 0 H M ) = 0 .
Thus, we replace ρ F H ϕ with μ 0 M and write the remaining inequality in the form
ρ T ϕ · T ˙ + ( ρ E ϕ + J 1 T R R + μ 0 F 1 H F 1 M ) · E ˙ 0
This result indicates that we let
T : = T R R + J μ 0 F 1 H F 1 M
so that we have
ρ R T ϕ · T ˙ + ( T ρ R E ϕ ) · E ˙ = J θ σ .
Observe that F 1 H and F 1 M are invariant vectors. To yield this
( F 1 H ) * = ( Q F ) 1 Q H = F 1 Q 1 Q H = F 1 H
and the like for F 1 M . Hence, T is an invariant tensor. Moreover, by letting M = J F 1 M we can write (19) in the form
μ 0 M = H ϕ .
If T ˙ and E ˙ are independent, then we obtain
T ϕ = 0 , T = ρ R E ϕ
and then σ = 0 . The equations in (24) are said to characterize magneto-hyperelastic materials. The results (23) and (24) allow us to write the incremental relations
μ 0 M ˙ = ρ R E H ϕ E ˙ ρ R H H ϕ H ˙ ,
T ˙ = ρ R E E ϕ E ˙ + ρ R H E ϕ H ˙ .
If, again, σ = 0 , but
T ϕ 0 , E ϕ 0 ,
then we have
ρ R T ϕ · T ˙ + ( T ρ R E ϕ ) · E ˙ = 0 ;
Equation (27) is said to characterize magneto-hypoelastic materials. In general, we can express T ˙ via the representation formula,
T ˙ = ( T ˙ · N ) N + ( I N N ) G ,
where G is any second-order tensor. Let N = T ϕ / | T ϕ | . According to (27), we have
T ˙ = ( T ρ R E ϕ ) · E ˙ ρ R | T ϕ | 2 T ϕ + ( I T ϕ | T ϕ | T ϕ | T ϕ | ) G
or
T ˙ = 1 ρ R | T ϕ | 2 [ T ϕ ( T ρ R E ϕ ) ] E ˙ + ( I T ϕ | T ϕ | T ϕ | T ϕ | ) G .
If σ 0 , then we have
ρ R T ϕ · T ˙ + ( T ρ R E ϕ ) , · E ˙ = J θ σ
and the representation formula for T ˙ generalizes to
T ˙ = ( T ρ R E ϕ ) · E ˙ J θ σ ρ R | T ϕ | 2 T ϕ + ( I T ϕ | T ϕ | T ϕ | T ϕ | ) G .
Definite forms of (28) are now established by having in mind fluid or solid behaviours.

5.1. Fluids

The interaction between deformation and magnetization can be modelled by letting both ϕ and σ depend on T and H . For definiteness, we let
ρ R ϕ = ρ R ϕ 0 ( θ ) + 1 2 α | T | 2 1 2 ξ | H | 2 , α , ξ > 0 ,
and
J θ σ = β | T | 2 + γ | H | 2 , β , γ > 0 .
Equations (18) and (19) hold. According to (19), the dependence of ϕ on H results in the magnetization relation,
μ 0 M : = μ 0 J F 1 M = ξ H .
Since
E ϕ = 0 , ρ R T ϕ = α T
then (28) can be written in the form
T ˙ = T T α | T | 2 E ˙ β α T γ | H | 2 α | T | 2 T + ( I T | T | T | T | ) G .
If, for definiteness, we let G = E ˙ / α , then it follows
T ˙ = 1 α E ˙ β α T γ | H | 2 α | T | 2 T .
In particular, the choice γ = γ 0 | T | 2 yields
T ˙ = a E ˙ b T , a = 1 / α > 0 , b = β / α + γ 0 | H | 2 / α > 0 ,
that is a generalization of the Maxwell (fluid) model. Given T 0 = T ( 0 ) , if | H | is considered known, then the linear ODE (30) with non-constant coefficients can be solved on ( 0 , T ) to obtain
T ( t ) = T 0 exp 0 t b ( λ ) d λ + 0 t a ( τ ) exp τ t b ( λ ) d λ E ˙ ( τ ) d τ .
If E ˜ is a given past history on ( , 0 ] , we assume
T 0 = 0 a ( τ ) exp τ 0 b ( λ ) d λ E ˜ ˙ ( τ ) d τ .
By substituting T 0 in (31), we have
T ( t ) = t a ( τ ) exp τ t b ( λ ) d λ E ˙ ( τ ) d τ = 0 a ( t s ) exp 0 s b ( ξ ) d ξ E ˙ ( t s ) d s ,
where E ( τ ) = E ˜ ( τ ) , τ 0 .
We observe that the memory kernel of (31) has the form
G ( τ , t τ ) : = a ( τ ) exp τ t b ( λ ) d λ = a ( τ ) exp 0 t τ b ( ξ ) d ξ
which describes an aging effect as the G function changes over time due to the presence of the a ( τ ) factor. We obtain a standard memory kernel
G ( t τ ) = a 0 exp τ t b ( λ ) d λ
if a = a 0 is constant.
Things are different if σ is viewed as a viscous term in the form
J θ σ = ν | E ˙ | 2 + ζ | H ˙ | 2 , ν , ζ > 0 .
It follows that
T ˙ = T T α | T | 2 E ˙ ν | E ˙ | 2 + ζ | H ˙ | 2 α | T | 2 T + ( I T | T | T | T | ) G .
If G = E ˙ / α , we have
T ˙ = E ˙ α ν | E ˙ | 2 + ζ | H ˙ | 2 α | T | 2 T .
Observe that
H ˙ = F ˙ T H + F T H ˙ = F T ( L T H + H ˙ )
Now,
L T H + H ˙ = H + D H ,
where H is the corotational derivative, H = H ˙ W H . Hence, we have
| H ˙ | 2 = ( H + D H ) · F F T ( H + D H ) .
Equation (32) can then be written in the form
T ˙ = T T α | T | 2 E ˙ ν | E ˙ | 2 α | T | 2 T ζ ( H + D H ) · F F T ( H + D H ) α | T | 2 T + I T | T | T | T | G .

5.2. Solids

Solids are characterized by a stress dependence such that, asymptotically, T = G E . Hence, we formally replace the T of the fluid model with T G E . Define
ρ R ϕ = ρ R ϕ 0 ( θ ) + 1 2 E G E + 1 2 ( T G E ) · A ( T G E ) 1 2 H · Ξ H
and
J θ σ = ( T G E ) · [ ( β + H · Γ H ) 1 A ( T G E ) ] , β > 0 ,
where Ξ , Γ Sym + and A , G are positive-definite fourth-order tensors. Observe that
ρ R H ϕ = Ξ H , ρ R T ϕ = A ( T G E ) , ρ R E ϕ G E G T A ( T G E ) .
Moreover, according to (19), it follows that M : = J F 1 M = μ 0 1 Ξ H .
We now apply the representation (28) by letting
N = T ϕ | T ϕ | = A ( T G E ) | A ( T G E ) | .
Within ( I N N ) G , the representation formula yields
T ˙ [ T G E + G T A ( T G E ) ] E ˙ ( T G E ) · [ β + H · Γ H ] 1 A ( T G E ) | A ( T G E ) | N A ( T G E ) · [ A 1 + G E ) ] E ˙ A ( T G E ) · [ β + H · Γ H ] 1 ( T G E ) | A ( T G E ) | N .
Hence, we have
T ˙ = ( N N ) [ A 1 + G ] E ˙ ( N N ) [ β + H · Γ H ] 1 ( T G E ) + ( I N N ) G .
Choosing, e.g.,
G = [ A 1 + G ] E ˙ [ β + H · Γ H ] 1 ( T G E )
we find
T ˙ = [ A 1 + G ] E ˙ [ β + H · Γ H ] 1 ( T G E ) ,
whence
T ˙ G E ˙ + [ β + H · Γ H ] 1 ( T G E ) = A 1 E ˙ .
Equation (33) shows that T G E evolves with a relaxation time
τ = β + H · Γ H .
Moreover, if E ˙ = 0 , then, asymptotically, we have
T = G E ,
as we expected for a solid model. We finally note that (33) takes the usual form
T ˙ = G 0 E ˙ 1 τ ( T G E )
after letting G 0 = A 1 + G .

5.3. A One-Dimensional Example

Restrict attention to one-dimensional models associated with strain, applied traction, and magnetic field in the direction e such that
E = E e e , T = T e e , H = H e .
The symbol T for the component of T is consistent with the engineering stress considered in the literature to be the ratio of the axial force over the reference area. Moreover, let G = G I , G 0 = G 0 I , and γ = e · Γ e ; Equation (34) can be written as
T ˙ = G 0 E ˙ 1 τ ( T G E ) , τ = β + γ H 2 ,
where G 0 , G , β , and γ > 0 . Assume G = G 0 τ / 4 , β = 0.1 , and γ = 1 . Then, the traction response T for a given sinusoidal strain E ( t ) = λ sin ( π 20 t ) is plotted in Figure 1 under different values of the magnetic fields H. These results are in agreement with [6] Figures 4 and 5 as they predict that the increase in the magnetic field changes the orientation of the loops and widens the hysteresis due to greater energy dissipation.

6. Relation to Other Approaches

Magneto-viscoelasticity is a broad subject that accounts for the interaction between a magnetic field and deformation while both elastic and dissipative effects are allowed. In this framework, various approaches have been developed. The great majority of them can be characterized according to the modelling of dissipation.
In [11,14,15,16,17], mechanical viscous effects are described by assuming the existence of an intermediate configuration that is related to the current configuration by an elastic deformation and to the initial configuration by a purely viscous deformation. Hence, the deformation gradient is given a multiplicative decomposition
F = F e F v .
Instead, the magnetic induction is assumed in the form
B = B e + B v .
The Cauchy–Green tensor C is considered in the form
C = F T F = F v T F e T F e F v = F v T C e F v ;
hence, both C e and C v are defined, but C C e C v . As is often the case in non-equilibrium processes, the distinction between B e and B v is based on the observation that, upon the sudden application of a constant magnetic induction, the magnetic field generated inside the material starts from an initial non-equilibrium value and then evolves to approach an equilibrium value. To model these effects, the existence of a dissipation mechanism is assumed for the magnetic induction as well. The additive decomposition of B is further motivated by the vector character B .
The entropy inequality in [11], Equation (15),
ρ Ψ ˙ + T · D M · B ˙ 0
is consistent with (2) in that
M · B ˙ = M · μ 0 ( M ˙ + H ˙ )
and, hence, we have
ρ ψ ˙ + T · D μ 0 M · H ˙ 0 , ρ ψ = ρ Ψ + 1 2 μ 0 M 2 .
Next, a free energy Ω is considered to be a function of the invariant fields C , C v , B , B v ; we observe that F B of [11] is just the magnetic induction here denoted by B . A reduced dissipation inequality then follows in the form
C v Ω · C ˙ v + B v Ω · B ˙ v 0 .
Though the approach is deeply different from ours, bearing in mind equation (12), it seems natural to view the roles of C v and B v as the analogues of the dependence on the histories of E and H .

6.1. Incremental Magnetoelastic Equations

A simpler model is established in [18] in terms of the Lagrangian counterparts of B and H , i.e., B = J F 1 B and H = F T H , denoted by B l and H l in [18]. Differently from the present approach, Otténio et al. start with a “modified free energy function” Ω ( F , B ) such that
T R = F Ω , H = B Ω ;
in components
T i K R = F i K Ω , H P = B P Ω .
The magnetic field and the deformation are then supposed to undergo changes. In the linear approximation, the changes Δ F , Δ B , and Δ T R , Δ H are related in the form
Δ T R = A Δ F + Γ Δ B , Δ H = Γ Δ F + K Δ B ,
where
A = F F Ω , Γ = F B Ω , K = B B Ω .
in components,
A i K j P = F i K F j P Ω , Γ i K P = F i K B P Ω = Γ P i K , K P Q = B P B Q Ω .
Now, we look at the increments as occurring smoothly in time ( C 1 functions), in the linear approximation
Δ T R = T ˙ R Δ t ,
and the like for the other terms, so that we can write
T ˙ R = A F ˙ + Γ B ˙ , H ˙ = Γ F ˙ + K B ˙ .
Equations (37) are formally equal to the relations in (23) of [18] where T ˙ R , F ˙ , B ˙ , and H ˙ are said to be infinitesimal increments. There is no abuse in regarding T ˙ R , F ˙ , B ˙ , and H ˙ as time derivatives. With this view the formulas in (37) have some similarity with Equations (25) and (26). Indeed, the similarity is easily clarified once we observe that the starting assumption (12) of [18] can be derived here from the entropy inequality as the elastic part of the constitutive equations.

6.2. Visco-Hyperelastic Constitutive Modelling

Mainly in connection with fluid-structure coupling problems, the stress tensor is often expressed in terms of both the right and the left Cauchy–Green tensors F T F , F F T (see, e.g., [19,20] and the references therein). The corresponding approaches lead to neither memory functionals nor rate equations, yet application of the idea underlying the visco-hyperelastic models within a simple version of the present setting is of interest.
For simplicity, we neglect the heat conduction ( q = 0 ) and the dependence on the temperature gradient ( θ ). Let θ , F , H , and D be the independent variables. Since the free energy ϕ is invariant, we assume the particular dependence ϕ = ϕ ( θ , E , H , D ) . We then prove the thermodynamic consistency of a stress tensor in the form
T = T ^ ( θ , F , H ) + T d ( D )
where
T d Sym , T d ( D ) 0 as D 0 .
Consider the entropy inequality (2). Upon computation of ϕ ˙ , we have
ρ ( θ ϕ + η ) θ ˙ ρ E ϕ · E ˙ ρ H ϕ · H ˙ ρ D ϕ · D ˙ + T ^ · L + T d · D μ 0 M · H ˙ = θ σ 0 .
The arbitrariness (and linearity) of θ ˙ and D ˙ implies
η = θ ϕ , D ϕ = 0 .
In view of (5), since L = D + W , we can write the inequality in the form
( T ^ ρ F E ϕ F T ρ H F H ϕ ) · D + T d · D + ( T ^ + ρ F H ϕ H ) · W ( μ 0 M + ρ F H ϕ ) · H ˙ = θ σ 0 .
The arbitrariness of W Skw and H ˙ implies
T ^ + ρ F H ϕ H Sym
μ 0 M = ρ F H ϕ ,
where the relation in (41) is formally equal to (6). Similar to the arbitrariness of D in (39), we conclude that the linear part is required to vanish, whence
sym [ T ^ ρ F E ϕ F T ρ H F H ϕ ] = 0 ,
while the nonlinear part is non-negative, T d · D 0 . According to (40)–(42), we can write the stress T ^ in the form
T ^ = ρ F E ϕ F T μ 0 H M , T ^ + μ 0 H M Sym .
If
T d = 2 μ D + λ ( tr D ) 1
then
0 θ σ = T d · D = 2 μ D · D + λ ( tr D ) 2
whence we obtain the classical inequalities μ 0 , 2 μ + 3 λ 0 of the viscosity coefficients.
Note that the tensor T ^ + μ 0 H M of (43) is the Eulerian analogue of the tensor T introduced in Section 4.

6.3. Rheological Equations and Relaxing Media

The use of E and H in the present models is motivated by the invariance character so that E ˙ and H ˙ are also invariant and then objective. The Eulerian descriptions might involve the left Cauchy–Green tensor F F T , rather than the right one F T F = 1 + 2 E . In this connection, we observe that, if ϕ depends on E through scalar invariants, then E ϕ in (43) is replaced by one or more scalars, say g, and
F E ϕ F T = g F F T ,
thus providing the dependence of T ^ on the left Cauchy–Green tensor F F T .
The elastic properties of nonlinear materials are often modelled in terms of F F T , while dissipativity is mainly described through the stretching tensor D . This is the case, e.g., in [21] where the Cauchy stress T is given the form
T = T E ( F F T ) + T V ,
where T V is governed by the rate equation
T V + λ D T V = ν D ,
in which D denotes the objective derivative; ref. [21] considers the Oldroyd derivative, T V = T ˙ V L T V T V L T , the Jaumann derivative T V = T ˙ V W T V + T V W , and the Cotter–Rivlin derivative T V = T ˙ V + L T T V + T V L . Equation (44) is the analogue of (38).
It is worth remarking that in Oldroyd’s model of rheological equations for fluids, the viscous stress T V is subject to Equation ([22], Equation (59))
T V + λ T V = 2 μ D + 2 μ τ D .
Viscoelastic fluid theories have also been established in terms of the left relative Cauchy–Green tensor, defined by
B t ( τ ) = F t ( τ ) F t ( τ ) T , F t ( τ ) = F ( τ ) F 1 ( t ) .
The constitutive equation is taken in the form [23]
T = η B 1 + β B 1 2 + ν B 2 ,
where
B 1 ( t ) = d d τ B t ( τ ) | τ = t , B 2 ( t ) = d 2 d τ 2 B t ( τ ) | τ = t .

7. Conclusions

The interaction between the deformation and magnetic field leads to a broad spectrum of phenomena. The most common scheme is merely that of stress and magnetization induced by the deformation gradient and the magnetic field, as is described by the linearized Equations (37) [18] or (25) and (26). However, the broad spectrum motivates an investigation of magneto-viscoelasticity, where both the equilibrium and dissipative properties of the interaction are modelled. This paper develops an extensive thermodynamically consistent framework and establishes new models by letting the constitutive properties involve memory functionals or rate equations.
The modelling through memory functionals involve both a dependence on the present value and the history of deformation, magnetic field, and temperature gradient. We then find that the relations for stress and magnetization are given by derivatives of the free energy with respect to the present values of the strain and magnetic field. Instead, the history dependence allows the description of the dissipative effects of the stress, magnetization, and heat conduction. Some relevant examples of the functionals are determined; Equation (11) is a linearized functional for the magnetization in terms of the history of the magnetic field.
Materials with a non-instantaneous response are often described through fractional-order derivatives [24,25,26] in place of the memory functional. This involves an unbounded kernel of the power-law form, which makes the thermodynamic consistency questionable. However, we are unaware of genuine models with memory functionals for the magnetic field through fractional-order derivatives.
The approach through rate equations turn out to be advantageous thanks as well to some ideas and results established in [8], in which the entropy production is viewed as a constitutive function, while a representation formula yields the general result (28) for solids and the analogue for fluids. The completeness of the scheme allows us to find the stress rate in magneto-hypoelasticity, magneto-hyperelasticity, and (dissipative) magneto-viscoelasticity. Equation (28) yields the rate T ˙ of the stress T in terms of the rate E ˙ ; the rate is characterized by the free energy ϕ ( θ , T , E , H ) and the entropy dissipation σ . The example devised in Section 5.3 predicts that the increase in the magnetic field broadens the hysteretic loops (see Figure 1). This is consistent with the view that the dissipation increases proportionally to the friction caused by inter-particle magnetic attraction [27]. The dual constitutive assumption ϕ = ϕ ( θ , E , M , H ) models the hysteresis effects in ferromagnetic materials, parameterized by the temperature θ and the strain E , as well as with hysteresis in ferroelectrics [28].
By comparing the two approaches, we can say that the dependence on histories proves convenient for linear models (non-linearities involve multiple integrals and, hence, are technically quite difficult to manage). Rate-type models are more flexible and easily allow non-linear models. However, both approaches have proved to be thermodynamically consistent.
The analysis in Section 6 of some of the approaches developed in the literature allows us to find analogies to, and differences from, the present scheme. First, there are approaches in which the strain is split into elastic and viscous parts; the present description of deformation is decomposition-free. Secondly, as expected, the instantaneous response of memory functionals or the non-dissipative version of the rate-type equations yields the analogue of magnetoelastic models (incremental equations). Thirdly, with the restriction to the instantaneous response of the magnetization, it is shown that a model holds true in which the additive terms of the stress describe the dissipative effects of the deformation within an Eulerian framework.

Author Contributions

Investigation, A.M. and C.G. All authors have contributed substantially and equally to the work reported. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The research leading to this work has been developed under the auspices of INDAM-GNFM.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Eringen, A.C.; Maugin, G.A. Electrodynamics of Continua I—Foundations and Solid Media; Springer: New York, NY, USA, 1990. [Google Scholar] [CrossRef]
  2. Ogden, R.W.; Steigmann, D.J. Mechanics and Electrodynamics of Magneto- and Electro-Elastic Materials; CISM Courses and Lectures 527; Springer: Vienna, Austria, 2011; ISBN 9783709107003. [Google Scholar]
  3. Bustamante, R.; Shariff, M.; Hossain, M. Mathematical formulations for elastic magneto-electrically coupled soft materials at finite strains; time-independent processes. Int. J. Eng. Sci. 2021, 159, 103429. [Google Scholar] [CrossRef]
  4. Hossain, M.; Steinmann, P. Modelling electro-active polymers with a dispersion-type anisotropy. Smart Mater. Struct. 2018, 27, 025010. [Google Scholar] [CrossRef]
  5. Agirre-Olabide, I.; Lion, A.; Elejabarrieta, M.J. A new three-dimensional magneto-viscoelastic model for isotropic magnetorheological elastomers. Smart Mater. Struct. 2017, 26, 035021. [Google Scholar] [CrossRef]
  6. Haldar, K. Constitutive modeling of magneto-viscoelastic polymers, demagnetization correction, and field induced Poynting effect. Int. J. Eng. Sci. 2021, 165, 103488. [Google Scholar] [CrossRef]
  7. Wang, B.; Kari, L. Constitutive model of isotropic magneto-sensitive rubber with amplitude, frequency, magnetic and temperature dependence under a continuum mechanics basis. Polymers 2021, 13, 472. [Google Scholar] [CrossRef] [PubMed]
  8. Giorgi, C.; Morro, A. Nonlinear models of thermo-viscoelastic materials. Materials 2021, 14, 7617. [Google Scholar] [CrossRef]
  9. Giorgi, C.; Morro, A. A thermodynamic approach to rate-type models of elastic-plastic materials. J. Elast. 2021, 147, 113–148. [Google Scholar] [CrossRef]
  10. Rajagopal, K.R.; Ružička, M. Mathematical modeling of electrorheological materials. Cont. Mech. Thermodyn. 2001, 13, 59–78. [Google Scholar] [CrossRef]
  11. Saxena, P.; Hossain, M.; Steinmann, P. A theory of finite deformation magneto-viscoelasticity. Int. J. Solids Struct. 2013, 50, 3886–3897. [Google Scholar] [CrossRef] [Green Version]
  12. Dorfmann, A.; Ogden, R.W. Magnetoelastic modelling of elastomers. Eur. J. Mech. A Solids 2003, 22, 497–507. [Google Scholar] [CrossRef]
  13. Garcia-Gonzalez, D. Magneto-visco-hyperelasticity for hard-magnetic soft materials: Theory and numerical applications. Smart Mater. Struct. 2019, 28, 085020. [Google Scholar] [CrossRef]
  14. Garcia-Gonzalez, D.; Hossain, M. A microstructural-based approach to model magneto-viscoelastic materials at finite strain. Int. J. Solids Struct. 2021, 208–209, 119–132. [Google Scholar] [CrossRef]
  15. Lubliner, J. A model of rubber viscoelasticity. Mech. Res. Comm. 1985, 12, 93–99. [Google Scholar] [CrossRef]
  16. Reese, S. A micromechanically motivated material model for the thermo-viscoelastic material behaviour of rubber-like polymers. Int. J. Plast. 2003, 19, 909–940. [Google Scholar] [CrossRef]
  17. Reese, S.; Govindjee, S. A theory of finite viscoelasticity and numerical aspects. Int. J. Solids Struct. 1998, 35, 3455–3482. [Google Scholar] [CrossRef]
  18. Otténio, M.; Destrade, M.; Ogden, R.W. Incremental magnetoelastic deformations, with application to surface instability. J. Elast. 2008, 90, 19–42. [Google Scholar] [CrossRef]
  19. Sugiyama, K.; Li, S.; Takeuchi, S.; Takagi, S.; Matsumoto, Y. A full Eulerian finite difference approach for solving fluid-structure coupling problems. J. Comp. Phys. 2011, 230, 596–627. [Google Scholar] [CrossRef]
  20. Upadhyay, K.; Subhash, G.; Spearot, D. Visco-hyperelastic constitutive modeling of strain rate sensitive soft materials. J. Mech. Phys. Solids 2020, 135, 103777. [Google Scholar] [CrossRef]
  21. Saccomandi, G.; Vianello, M. Shear waves in nonlinear relaxing media: A three-dimensional perspective. J. Acoust. Soc. Am. 2021, 149, 1589–1595. [Google Scholar] [CrossRef]
  22. Oldroyd, J.G. On the formulation of rheological equations of state. Proc. R. Soc. Lond. 1950, 200, 523–541. [Google Scholar] [CrossRef]
  23. Huilgol, R.R. Viscoelastic fluid theories based on the left Cauchy-Green tensor history. Rheol. Acta 1979, 18, 451–455. [Google Scholar] [CrossRef]
  24. Deseri, L.; Di Paola, M.; Zingales, M. Free energy and states of fractional-order hereditariness. Int. J. Solids Struct. 2014, 51, 3156–3167. [Google Scholar] [CrossRef]
  25. Jaishankar, A.; McKinley, G.H. Power-law rheology in the bulk and at the interface: Quasi-properties and fractional constitutive equations. Proc. R. Soc. A 2013, 469, 20120284. [Google Scholar] [CrossRef]
  26. Podlubny, I. Fractional Differential Equations; Academic Press: London, UK, 1998; ISBN 9780125588409. [Google Scholar]
  27. Hiptmair, F.; Major, Z.; Hasslacher, R.; Hild, S. Design and application of permanent magnet flux sources for mechanical testing of magnetoactive elastomers at variable field directions. Rev. Sci. Instrum. 2015, 86, 085107. [Google Scholar] [CrossRef]
  28. Giorgi, C.; Morro, A. A thermodynamic approach to rate-type models in deformable ferroelectrics. Cont. Mech. Thermodyn. 2021, 33, 727–747. [Google Scholar] [CrossRef]
Figure 1. Cyclic responses in the ( E , T ) -plane under sinusoidal stretching with λ = 3 , G 0 = 4 and H = 0.1 (dashed), H = 0.4 (solid), H = 0.7 (dotted).
Figure 1. Cyclic responses in the ( E , T ) -plane under sinusoidal stretching with λ = 3 , G 0 = 4 and H = 0.1 (dashed), H = 0.4 (solid), H = 0.7 (dotted).
Materials 15 06699 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Giorgi, C.; Morro, A. Magneto-Viscoelastic Materials: Memory Functionals and Rate Equations. Materials 2022, 15, 6699. https://doi.org/10.3390/ma15196699

AMA Style

Giorgi C, Morro A. Magneto-Viscoelastic Materials: Memory Functionals and Rate Equations. Materials. 2022; 15(19):6699. https://doi.org/10.3390/ma15196699

Chicago/Turabian Style

Giorgi, Claudio, and Angelo Morro. 2022. "Magneto-Viscoelastic Materials: Memory Functionals and Rate Equations" Materials 15, no. 19: 6699. https://doi.org/10.3390/ma15196699

APA Style

Giorgi, C., & Morro, A. (2022). Magneto-Viscoelastic Materials: Memory Functionals and Rate Equations. Materials, 15(19), 6699. https://doi.org/10.3390/ma15196699

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