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
, the magnetic field
, and the temperature gradient
,
as well as the histories, up to time
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 , , 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
, the Green–Saint Venant tensor
, and a magnetic field vector be the independent variables. The generality is allowed by the joint dependence on
and
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
. 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 . The motion is described by means of the function , providing the position vector . The symbols denote the gradient operator with respect to , . The function is assumed to be differentiable; hence, we can define the deformation gradient as or, in suffix notation, . The invertibility of is guaranteed by letting . For any tensor , we define as .
2. Balance Equations
Let
be the velocity field on
. A superposed dot denotes the time differentiation following the motion of the body, and hence, for any function
, we have
. We denote by
the velocity gradient,
, and recall that
The right Cauchy–Green tensor
and the Green–Saint Venant deformation tensor
are defined by
where
is the second-order identity tensor. Moreover,
denotes the stretching tensor,
, and
the spin tensor,
.
Let
be the internal energy density (per unit mass),
the symmetric Cauchy stress,
the heat flux vector,
the mass density,
r the external heat supply, and
the mechanical body force per unit mass. Let
be the magnetization per unit mass and
the magnetic intensity. The balance equations for mass and linear momentum are written in the form (see, e.g., [
11,
12,
13])
where
is the force per unit volume of magnetic character,
in stationary conditions, we can take
. Additionally, for later convenience, we let
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
Let
be the entropy density and
the absolute temperature. As for the statement of the second law of thermodynamics, we let the inequality
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
, and the other functions occurring in the balance equations.
In terms of the magnetic Gibbs free energy
the entropy inequality can be written in the form
Based on (
2), next, we describe the magneto-viscoelasticity by the memory functionals or rate equations and examine the thermodynamic consistency.
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
are among the independent variables. Moreover, the time delay in response motivates the dependence on the histories
. Hence, we let
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
and
are not invariant. Under a change of frame
and
change as vectors,
Instead,
are invariant in that
Of course,
is also invariant,
Hence, the dependence on the pair
,
does not include that on
,
. For definiteness, we let
be the pair of functions
As for the dependence on
, we observe that
Hence,
is a vector while the referential gradient
is invariant.
Let
, and
be differentiable as
for any
. Define the constant continuation of
, and
as
and the like for
and
while
on
.
While
is the set of independent variables, objectivity and modelling purposes indicate that we let
With a small abuse of notation, the dependence on time is denoted with the same symbol
We first show that the free energy functional
is required to satisfy a minimum property. If
, and, hence,
, are constant, and
in the interval
, then (
2) implies that
and hence,
where we let
. By the continuity of the functional, as
we have
where
, and
are the constant histories with values of
, and
. Thus, it follows that
This is the content of the minimum property: among all histories and with the given present values of , and , none yields a smaller value of the free energy than that corresponding to the constant histories and .
4. Thermodynamic Restrictions
We now compute the time derivative
and substitute in the entropy inequality to obtain
Here, we have made the dependence on
and
explicit. The linearity and arbitrariness of
imply
Since
, the occurrence of
in the inequality is through
and this quantity has to vanish. Now
and then
identically. Hence, we have
which implies
We now can write the remaining inequality in the form
The linearity and arbitrariness of
imply
Inequality (
7) then reduces to
According to (
6) and (
8), it follows that
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
through
. Instead, the dependence of
on
h leaves
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,
. Moreover, nonlinear constitutive equations for
and
are established by selecting the partial derivatives
Now, given
, a fully symmetric fourth-order tensor-valued function, we have
and the like for a symmetric second-order tensor-valued function
. Hence, letting
and
be the second- and fourth-order unit tensors, we define
so that
Hence, the sought functional
takes the form
and the functional
is thus far undetermined. According to (
9), the functional
has to satisfy the inequality
and this is eventually the requirement of the unknown functional
.
For definiteness, we consider
in the form
and, hence, the functional
can be written as
,
The assumption
implies
Likewise, letting
, we find
Hence, (
12) holds if and only if
for all
.
The minimum property of
at the constant histories
holds if and only if
for all
.
Hence, the functional is thermodynamically consistent if and , as well as and , have opposite types of definiteness.
Heat Conduction
Now we let
, and, for simplicity, we look for models where
is independent of
and
so that the reduced dissipation inequality splits into (
12) and
Multiply this inequality by
and observe that, using the referential heat flux
, we can write
For definiteness, let
be given by the constitutive functional
where
is a positive-definite second-order tensor, while
and
are so far undetermined; we only assume
.
We let
with
taken in the form
The minimum property of
, and hence of
, at
holds if and only if
. Now, inequality (
13) results in
An integration by parts yields
Hence, inequality (
15) can be written in the form
The linearity and arbitrariness of
imply that
Since
and
, then
for all
. Hence,
and, consistently, the kernel of the functional
is negative definite.
Hence, we can write
in the form
At the limit of short memory, we have
In the spatial description, the constitutive equation reads
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
need not equal
, and, moreover,
Since
, then
We start with the Eulerian description and write the Clausius–Duhem inequality in the form
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
where
is a stress-like variable to be identified. Consequently, the Clausius–Duhem inequality (
17) can be written as
Since
, then we have
Thus, we replace
with
and write the remaining inequality in the form
This result indicates that we let
so that we have
Observe that
and
are invariant vectors. To yield this
and the like for
. Hence,
is an invariant tensor. Moreover, by letting
we can write (
19) in the form
If
and
are independent, then we obtain
and then
. The equations in (
24) are said to characterize magneto-hyperelastic materials. The results (
23) and (
24) allow us to write the incremental relations
If, again,
, but
then we have
Equation (
27) is said to characterize magneto-hypoelastic materials. In general, we can express
via the representation formula,
where
is any second-order tensor. Let
. According to (
27), we have
or
If
, then we have
and the representation formula for
generalizes to
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
and
. For definiteness, we let
and
Equations (
18) and (
19) hold. According to (
19), the dependence of
on
results in the magnetization relation,
Since
then (
28) can be written in the form
If, for definiteness, we let
, then it follows
In particular, the choice
yields
that is a generalization of the Maxwell (fluid) model. Given
, if
is considered known, then the linear ODE (
30) with non-constant coefficients can be solved on
to obtain
If
is a given past history on
, we assume
By substituting
in (
31), we have
where
.
We observe that the memory kernel of (
31) has the form
which describes an aging effect as the
G function changes over time due to the presence of the
factor. We obtain a standard memory kernel
if
is constant.
Things are different if
is viewed as a viscous term in the form
Now,
where
is the corotational derivative,
. Hence, we have
Equation (
32) can then be written in the form
5.2. Solids
Solids are characterized by a stress dependence such that, asymptotically,
. Hence, we formally replace the
of the fluid model with
. Define
and
where
and
are positive-definite fourth-order tensors. Observe that
Moreover, according to (
19), it follows that
.
We now apply the representation (
28) by letting
Within
, the representation formula yields
Choosing, e.g.,
we find
whence
Equation (
33) shows that
evolves with a relaxation time
Moreover, if
, then, asymptotically, we have
as we expected for a solid model. We finally note that (
33) takes the usual form
after letting
.
5.3. A One-Dimensional Example
Restrict attention to one-dimensional models associated with strain, applied traction, and magnetic field in the direction
such that
The symbol
T for the component of
is consistent with the
engineering stress considered in the literature to be the ratio of the axial force over the reference area. Moreover, let
,
, and
; Equation (
34) can be written as
where
, and
. Assume
,
, and
. Then, the traction response
T for a given sinusoidal strain
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
Instead, the magnetic induction is assumed in the form
The Cauchy–Green tensor
is considered in the form
hence, both
and
are defined, but
. As is often the case in non-equilibrium processes, the distinction between
and
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
is further motivated by the vector character
.
The entropy inequality in [
11], Equation (
15),
is consistent with (
2) in that
and, hence, we have
Next, a free energy
is considered to be a function of the invariant fields
; we observe that
of [
11] is just the magnetic induction here denoted by
. A reduced dissipation inequality then follows in the form
Though the approach is deeply different from ours, bearing in mind equation (
12), it seems natural to view the roles of
and
as the analogues of the dependence on the histories of
and
.
6.1. Incremental Magnetoelastic Equations
A simpler model is established in [
18] in terms of the Lagrangian counterparts of
and
, i.e.,
and
, denoted by
and
in [
18]. Differently from the present approach, Otténio et al. start with a “modified free energy function”
such that
in components
The magnetic field and the deformation are then supposed to undergo changes. In the linear approximation, the changes
, and
are related in the form
where
in components,
Now, we look at the increments as occurring smoothly in time (
functions), in the linear approximation
and the like for the other terms, so that we can write
Equations (
37) are formally equal to the relations in (23) of [
18] where
, and
are said to be infinitesimal increments. There is no abuse in regarding
, and
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
(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 (
) and the dependence on the temperature gradient (
). Let
, and
be the independent variables. Since the free energy
is invariant, we assume the particular dependence
. We then prove the thermodynamic consistency of a stress tensor in the form
where
Consider the entropy inequality (
2). Upon computation of
, we have
The arbitrariness (and linearity) of
and
implies
In view of (
5), since
, we can write the inequality in the form
The arbitrariness of
and
implies
where the relation in (
41) is formally equal to (
6). Similar to the arbitrariness of
in (
39), we conclude that the linear part is required to vanish, whence
while the nonlinear part is non-negative,
. According to (
40)–(
42), we can write the stress
in the form
If
then
whence we obtain the classical inequalities
of the viscosity coefficients.
Note that the tensor
of (
43) is the Eulerian analogue of the tensor
introduced in
Section 4.
6.3. Rheological Equations and Relaxing Media
The use of
and
in the present models is motivated by the invariance character so that
and
are also invariant and then objective. The Eulerian descriptions might involve the left Cauchy–Green tensor
, rather than the right one
. In this connection, we observe that, if
depends on
through scalar invariants, then
in (
43) is replaced by one or more scalars, say
g, and
thus providing the dependence of
on the left Cauchy–Green tensor
.
The elastic properties of nonlinear materials are often modelled in terms of
, while dissipativity is mainly described through the stretching tensor
. This is the case, e.g., in [
21] where the Cauchy stress
is given the form
where
is governed by the rate equation
in which
denotes the objective derivative; ref. [
21] considers the Oldroyd derivative,
, the Jaumann derivative
, and the Cotter–Rivlin derivative
. Equation (
44) is the analogue of (
38).
It is worth remarking that in Oldroyd’s model of rheological equations for fluids, the viscous stress
is subject to Equation ([
22], Equation (59))
Viscoelastic fluid theories have also been established in terms of the left
relative Cauchy–Green tensor, defined by
The constitutive equation is taken in the form [
23]
where
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
of the stress
in terms of the rate
; the rate is characterized by the free energy
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
models the hysteresis effects in ferromagnetic materials, parameterized by the temperature
and the strain
, 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.