Next Article in Journal
Asymmetric Twisting of Coronal Loops
Next Article in Special Issue
A Non-Second-Gradient Model for Nonlinear Electroelastic Bodies with Fibre Stiffness
Previous Article in Journal
Dynamical Structures of Multi-Solitons and Interaction of Solitons to the Higher-Order KdV-5 Equation
Previous Article in Special Issue
Mathematical Modeling of Spherical Shell-Type Pattern of Tumor Invasion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Generalised Time-Dependent Mathematical Formulation for Magnetoelectrically Coupled Soft Solids at Finite Strains

by
Mohd Halim Bin Mohd Shariff
1,*,
Roger Bustamante
2 and
Mokarram Hossain
3
1
Department of Mathematics, Khalifa University of Science and Technology, Abu Dhabi 127788, United Arab Emirates
2
Departamento de Ingeniería Mecánica, Universidad de Chile, Santiago 837-0456, Chile
3
Zienkiewicz Centre for Computational Engineering, College of Engineering, Swansea University, Swansea SA1 8EN, UK
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(3), 628; https://doi.org/10.3390/sym15030628
Submission received: 30 January 2023 / Revised: 18 February 2023 / Accepted: 20 February 2023 / Published: 2 March 2023
(This article belongs to the Special Issue Symmetry in Finite Element Modeling and Mechanics)

Abstract

:
To date, the mechanical models of magnetoelectric couplings at finite strains have mainly been limited to time-independent constitutive equations. This paper enhances the literature by developing a time-dependent electromagnetic constitutive equation to characterise the mechanical behaviour of soft solids at finite strains and take into account the full form of the Maxwell equations. Our formulation introduces a symmetrical total stress and uses recently developed spectral invariants in the amended energy function; as a result, the proposed constitutive equation is relatively simple and is amenable to a finite-element formulation.

1. Introduction

Research in the areas of so-called smart and multifunctional materials has grown exponentially in recent years. Smart materials can be defined as the innovative substances that can alter their physical and mechanical attributes when exposed to one or more external stimuli. Some widely used external stimuli are temperature, humidity, light, pH, and acoustic, electric, and magnetic fields [1,2]. Among other smart materials, electroactive polymers (EAPs) and magnetoactive polymers (MAPs) have received unprecedented attention in recent years thanks to their myriad of potential applications. These ever-expanding applications have reached many areas, including soft robotics for targeted drug delivery and cancer therapy, flexible stretch-based sensors for wearable devices, materials for morphing and shape-shifting structures, key ingredients for rapidly expanding metamaterials, to mention a few [1,2,3,4,5]. For EAPs, while working as the so-called dielectric elastomers (DEs), an electric voltage applied along the thickness direction of a thin structure will create Coulomb forces as a result of the attractions of two opposite charges resulting in expansions in the lateral directions [5]. This mechanism of converting electric input to a mechanical output is the key in actuators for soft robotics. Furthermore, EAPs can be used in energy harvesting using ambient motions, in which a mechanical input creates electric outputs that could be the essential ingredients for creating clean and green sustainable energy.
In MAPs, a polymeric composite filled with magnetisable particles (micro or nanosized) can be activated upon the application of a remotely controlled magnetic field. Depending on the filler particles, MAPs can be decomposed into two groups, i.e., soft-magnetic MAPs, where particles have less residual magnetisation and coercivity, while in hard-magnetic MAPs, particles retain a high coercivity with a significant remanence magnetisation after the removal of the external magnetic field, see, for example, reference [6]. Note that soft-magnetic MAPs are mainly used in areas where mechanical properties such shear and loss moduli are to be tuned, while for applications in largely deformable and flexible shape-shifting structures, hard-magnetic MAPs are used. For various matrix materials, fillers, manufacturing techniques, experimental characterisations, and the computational modelling of both hard and soft MAPs, a few recent review papers can be consulted, see, for example, [6].
The experimental study, mathematical modelling and simulations, and the search for novel applications of EAPs and MAPs have so far been performed separately. However, very recently, following the so-called multiferroic materials in which a magnetoelectric coupling occurs simultaneously, the quest for largely deformable similar materials has increased manifolds where magnetoelectric coupling is a common phenomenon. Note that ceramic and metal-based magneto-electroactive hard materials have applications in producing magnetoelectric random access memories [7,8]. An interesting counterpart to these hard materials could be a soft polymeric material (elastomer or hydrogel), in which a polymeric matrix is filled with magnetisable (e.g., iron, neodymium–iron–boron) particles [9,10]. A typical application of such material is in hyperthermia for cancer patients, in which a remotely controlled dynamic magnetic field will create localised heat as a result of vibrating particles in a conductive polymeric matrix. All these nascent applications require the development of appropriate mathematical frameworks to simulate magnetoelectric coupling phenomena at finite strains.
Until recently, the mathematical modelling of magnetoelectric couplings at finite strains has been limited, in the sense that models have discarded the time-dependence of the electro-magneto-mechanically coupled problems in formulating the governing equations. In this contribution, we aim to present a generalised mathematical formulation for magnetoelectrically coupled soft materials at finite strains that takes into account time-dependent dynamical phenomena. Following our recent works [11,12] on electromagnetic–elastic materials including the modelling of residual stresses appearing in magnetoelectric elasticity, we briefly present the relevant and key governing equations that are essential to represent a magneto-electro-mechanically coupled system at large strains. Afterwards, the framework is extended to incorporate the time-dependent dynamical aspects. Finally, we use an attractive set of spectral invariants [12] to mechanically characterise the proposed amended energy function.

1.1. Remark

There are many coexisting theories and results on the subject of electromagnetic continuum mechanics theory. The physics and mathematics of electromagnetic forces in continuum mechanics are complex and tedious, and it is beyond the scope of this paper to regurgitate them here. Readers are encouraged to read references [13,14,15] for the basic concepts of electromagnetic forces in deformable continua. In this paper, we only focus on a mathematical construction (which uses some fundamental equations and concepts for a solid continuum in the presence of electromagnetic fields based on the work of references [13,14,15]) to provide a relatively simple form of constitutive equation, which is beneficial to obtain boundary value problem results (especially with the aid of the finite element method) and for design purposes. We believe the mathematical construction and proofs developed in this article are not found in the literature and hence could be employed for future research.
An attractive set of spectral invariants and their derivatives are used in our proposed model; readers may not be familiar with them and their physical concept, since their development are only recent; see for example, references [16,17,18]. The concept and derivation of spectral models are rather involved, and hence it is beyond the scope of this paper to repeat them here. Readers are encouraged to read references [16,17,18] and references therein, for details on spectral modelling.

1.2. Kinematics

The descriptions of the deformation tensor F , the right-stretch tensor U , right Cauchy–Green tensor C = U 2 , and scalar J = det F > 0 , where det denotes the determinant of a tensor, can be found in Ogden [19]. Unless stated otherwise, in this paper, all subscripts i , j , and k take the values one, two, and three.

1.3. Electromagnetic Field Equations

Here, we briefly state some fundamental equations and concepts for a solid continuum in the presence of electromagnetic fields based on the work of references [13,14,15].

1.3.1. Maxwell Equations

We use the notation q e for the density of electric charge per unit of volume and J for the volume electric current. In the case q e = 0 and J = 0 , the Maxwell equations are (see [13,15]):
div B = 0 , curl E + B ˙ = 0 , div D = 0 , curl H = D ˙ ,
where B is the magnetic induction, E is the electric field, D is the electric displacement, H is the magnetic field, and ( ) ˙ is the time derivative keeping x constant, i.e., for example B ˙ = B ( x , t ) t , d i v and c u r l are the divergence and curl operators in the current configuration.
For a vacuum, the following relations between the electric variables and the magnetic variables are:
D = ε o E , B = μ o H ,
where ε o is the electric permittivity in vacuum, and μ o is the magnetic permeability in vacuum. For a condensed matter two extra fields can be defined, which are the electric polarization P and the magnetization M , i.e.,
D = ε o E + P , B = μ o ( H + M ) ,

1.3.2. Lagrangian Electric and Magnetic Variables

From [14,20,21] the expressions for the electric and the magnetic variables in the reference configuration are
E L = F T E , D L = J F 1 D , B L = J F 1 B , H L = F T H .

1.3.3. Continuity Conditions

The imposition of boundary conditions in electromagnetic boundary value problems is not trivial, due to the requirement of Maxwell Equation (1) to be satisfied in the body and also for the exterior surrounding, which is assumed to be a vacuum space. To make the boundary conditions clear (see also Figure 1), we introduce the notation
[ [ a ] ] = a ( o ) a ( i ) ,
where (o) is associated with the outside of the body and (i) is associated with the inside of the body, which is very close to the surface of the body. The boundary conditions for the electromagnetic fields require the following continuity conditions [14]
n · [ [ B ] ] = 0 , n × [ [ E ] ] = v n [ [ B ] ] , n · [ [ D ] ] = q e S , n × [ [ H ] ] + v n [ [ D ] ] = J S ,
where v n = v · n is the normal component of the velocity field, n is the outward unitary normal vector to B t (the surface of the body in the current configuration), q e S represents a surface distribution of free charges, and J S is a surface distribution of the electric current. From now on, we assume that q e S = 0 and J S = 0 .
At the boundary, the generally nonsymmetric Cauchy stress tensor σ satisfies the continuity condition
[ [ σ σ M ] ] n = t ^ ,
where t ^ is the external mechanical traction and the Maxwell stress tensor σ M has the relation [13,14]
σ M = 1 2 ε o E · E + 1 μ o B · B 2 M · B I ε o E E + P E + 1 μ o B B B M .

1.4. Balance of Mass and Equation of Motion

The mass balance equation is [19]
ρ ˙ + div ( ρ v ) = 0 ,
where v is the velocity and ρ is the density of the body. The density in the reference configuration is denoted ρ r and ρ = J 1 ρ r .
The first law of movement is [13]
ρ D v D t = d i v σ + f e + f ,
where D v D t is the acceleration and D D t is the total time derivative, f e is the body force due to the electromagnetic interactions and f is the body force that is independent of the electromagnetic fields. In the case q e = 0 and J = 0 , the body force f e is given by [13]
f e = ( g r a d E ) T P + ( g r a d B ) T M + P × B ¯ ˙ + d i v [ v ( P × B ) ] .
The second law of movement is stated as
ϵ σ = L ,
where σ is the Cauchy stress tensor. It is clear from Equation (12) that if the vector field L 0 , then the Cauchy stress σ is a nonsymmetric second-order tensor. We note that ϵ is the permutation tensor and, in this paper, we use the operator
ϵ ( b c ) = c × b ,
where b and c are vectors.

2. Constitutive Equations

2.1. The Amended Energy Function and the Total Stress Tensor

If ψ is the Helmholtz potential and ψ = ψ ( F , E e , B ) [13,14]), we have, for an isothermal problem, the Clausius–Duhem inequality, which takes the form
σ : D ρ ψ ˙ M e · B ˙ E ˙ e · P 0 ,
where the superposed dot represents the time derivative, : denotes the inner product of two second-order tensors, and D = g r a d v . In view of the relation σ : D = F 1 σ T : F ˙ , (14) takes the form
F 1 σ ρ ψ F T : F ˙ ρ ψ B + M e · B ˙ ρ ψ E e + P · E ˙ e 0 .
Since F ˙ , B ˙ , and E ˙ e are arbitrary, we obtain the relations
σ = ρ F ψ F , P = ρ ψ E e , M e = ρ ψ B ,
where
E e = E + E ¯ , E ¯ = v × B , M e = M + M ¯ , M ¯ = v × P .
and we use the convention ψ F α i = ψ F i α [19].
In order to define a “total symmetric stress tensor”, let us consider the following relations. First, notice that from Equation(1) and (3) we have
div D E = grad E D = ε o grad E E + grad E P , div E · E 2 I = g r a d E T E .
Using the relations
grad E grad E T E = curl E × E , grad E g r a d E T P = curl E × P ,
we have
div D E 1 2 ε o E · E I = grad E T P + curl E × P + ε o curl E × E .
In view of (1), we have the relation curl E = B ˙ and from (18), we obtain
div D E 1 2 ε o E · E I = grad E T P B ˙ × P ε o B ˙ × E = grad E T P B ˙ × D ,
thus
grad E T P = div D E ε o 2 ( E · E ) I + B ˙ × D .
On the other hand, using the relations
div B B = grad B B , div B · B 2 I = grad B T B ,
grad B grad B T B = curl B × B ,
we have, in view of (1) and (3):
1 μ o div B B B · B 2 I = 1 μ o curl B × B = D ˙ × B + c u r l M × B .
Therefore, from (1), we have the relations
div M · B I = grad B T M + grad M T B , div B × M = grad M B .
Using the relation
grad M g r a d M T B = curl M × B ,
we obtain
div M · B I B M = grad B T M c u r l M × B .
and from (25) and (28), we have
grad B T M = div 1 μ o B B 1 2 B · B I + M · B I B M D ˙ × B .
To facilitate our formulation, we let
ψ ( F , E e , B ) = ψ F , F T E L e , 1 J F B = ϕ ( U , B L , E L e ) ,
where
E L e = F T E e .
We note that ϕ is objective [19]. Consider the relations
E e F = F T E e , B F = 1 J F B L F 1 + i = 1 3 e i B L e i ,
where { e i } is the unit orthogonal vector basis in Cartesian coordinates in the current configuration.
Using the previous expressions in (16), we prove that
σ = 2 ρ F ϕ C F T P E e M e · B I + B M e .
Proof. 
From (30) (see (31) and (4)) and using the chain rule, we obtain,
F ϕ F = F ψ F + ψ B B F + ψ E e E e F .
The mathematical derivations below require the operator
q ( r s w ) = ( q · r ) ( s w ) ,
where q , r , s , and w are vectors. From (16), we have
ρ F ψ F = 2 ρ F ϕ C F T = σ ( F M e ) B F ( F P ) E e F ,
and using (35), we obtain,
2 ρ F ϕ C F T = σ ( F M e ) 1 J F B L F 1 + i = 1 3 e i B L e i + ( F P ) F T E e .
We then obtain the simplified form
2 ρ F ϕ C F T = σ + M e · B I B M e + P E e ,
which implies that σ = 2 ρ F ϕ C F T P E e M e · B I + B M e . □
We prove that P v × B = B M ¯ M ¯ · B I + v P × B . Considering Cartesian coordinates, let
A = P v × B A i j = P i ϵ i r s v r b s ,
and
G = B M ¯ M ¯ · B I + v P × B G i j = B i ϵ j r s v r P s ϵ k r s v r P s B k δ i j + v i ϵ j r s P r B s .
For the case i j , in particular when i = 1 and j = 2 , we have
A 12 = P 1 ϵ 213 v 1 B 3 + P 1 ϵ 231 v 3 B 1 = G 12 ,
and the same results can be obtained for ( i , j ) = ( 1 , 3 ) , ( i , j ) = ( 2 , 1 ) , etc., then for i j , we have that A i j = G i j . For the case i = j , for example when i = j = 1 , we have
A 11 = P 1 ϵ 123 v 2 B 3 + P 1 ϵ 132 v 3 B 2 = G 11 ,
and similar results are obtained for the cases i = j = 2 and i = j = 3 .
Using the relation
P ( v × B ) = B M ¯ M ¯ · B I + v P × B ,
from (21) and (29), we have
grad E T P = div D E 1 2 ε o ( E · E ) I + B ˙ × D ,
grad B T M = div 1 μ o B B 1 2 ( B · B ) I + ( M · B ) I B M D ˙ × B .
For clarity, we restate (10) and (11), i.e.,
ρ a = div σ + f e + f ,
f e = ( grad E ) T P + ( grad B ) T M + P × B ¯ ˙ + div [ v ( P × B ) ] .
Noticing that
σ = 2 ρ F ϕ C F T P E e M e · B I + B M e ,
then, from (42) and (43), considering (40), (41), and P = D ε o E , we obtain
ρ a = div 2 ρ F ϕ C F T + ε o E E 1 2 ( E · E ) I + 1 μ o B B 1 2 ( B · B ) I + div B M ¯ M ¯ · B + v ( P × B ) P ( v × B ) + P × B ¯ ˙ + B ˙ × D D ˙ × B + f ,
and using (39), the above equation becomes
ρ a = div T + P × B ¯ ˙ + B ˙ × D D ˙ × B + f ,
where
T = 2 ρ F ϕ C F T + 1 μ o B B 1 2 B · B I + ε o E E 1 2 E · E I .
Using P × B ¯ ˙ = P ˙ × B + B × P ˙ and (3), we have that
P × B ¯ ˙ + B ˙ × D D ˙ × B = ε o E × B ¯ ˙ ,
then, from (46), we obtain
ρ a = div T ε o E × B ¯ ˙ + f .
The expression ϵ σ = L given in (12) is given as follows: From (17) and (33), we have
ϵ σ = ϵ P E + v × B + ϵ B M e ,
therefore
ϵ σ = ϵ : P E e + ϵ B M e , = P × E + M × B + P × ( v × B ) + ( v × P ) × B ,
then, using the identity v × ( P × B ) = P × ( v × B ) + ( v × P ) × B in (51), we obtain
ϵ σ = P × E + M × B + v × ( P × B ) = L
and (12) is satisfied, while T T = T .
Using (30) in (16), we obtain
P = ρ F ϕ E L e , M e = ρ J F T ϕ B L ,
and from (3), considering that the independent electromagnetic variables are B and E , we have,
D = ε o E ρ F ϕ E L e , H = 1 μ o B + ρ J F T ϕ B L v × ρ F ϕ E L e .
We can further simplify the relation in (47) by letting
Ω = ρ 0 ϕ + 1 2 μ 0 J B · B + ε 0 2 J E L · C 1 E L
and following the work of Shariff [12], with some algebra, we have the simplified version
T = 2 J 1 F Ω C F T .
In the case of incompressible bodies, J = 1 and (56) becomes
T = 2 F Ω C F T p I ,
where p is a Lagrange multiplier associated with the incompressible constraint J = 1 .

2.2. Continuity Condition in Terms of the Total Stress

Using (57) in (33) and replacing that in (7), after some manipulations, we obtain
σ σ M = T + P ( B × v ) + B ( v × P ) + [ ( P × v ) · B ] I .
If we recall the use of the superscripts (o) and (i) for outside and inside the body, respectively, (7) becomes σ ( o ) σ M ( o ) σ ( i ) + σ M ( i ) n = t ^ . Considering that outside (in vacuum) P = 0 and using (58) in the above equation, we obtain T ( i ) + P ( i ) ( B ( i ) × v ) + B ( i ) ( v × P ( i ) ) + ( P ( i ) × v ) · B ( i ) I n = T ( o ) n + t ^ . Outside the body, for the sake of simplicity, we add the contribution of T ( o ) n to t ^ , and dropping the superscript (i), the continuity condition (7) becomes
T + P ( B × v ) + B ( v × P ) + ( P × v ) · B I n = t ^ ,
and for incompressible bodies, we have,
τ p I + P ( B × v ) + B ( v × P ) + ( P × v ) · B I n = t ^ , τ = 2 F Ω C F T .

3. Spectral Invariants

We are interested in applications concerning the behaviour of electromagnetic, active, incompressible, transversely isotropic, soft solids. We assume that the Helmholtz potential to have the form
ψ ( F , E e , B , a 0 ) = ψ F , F T E L e , 1 J F B , a 0 = Φ ( U , B L , E L e , a 0 ) ,
where in the absence of electromagnetic fields and deformations, a 0 is the preferred unit vector of a transversely isotropic solid. In this paper, we prefer to express Φ in terms of the following unit vectors [16],
e = 1 E L e E L e , b = 1 B L B L , E L e = E L e · E L e , B L = B L · B L
and express the energy function as
Φ ( U , b , e , a 0 , B L , E L e ) .
The invariance under rigid rotations of the body in the reference configuration then requires that
Φ ( U , b , e , a 0 , B L , E L e ) = Φ ( QUQ T , Q b , Q e , Q a 0 , B L , E L e )
for any orthogonal tensor Q ; hence, Φ can be expressed as an isotropic invariant of U , a 0 . Following the work of Shariff [12,16] (see also references therein), the spectral basis for the scalar-valued isotropic function Ω consists of only 14 spectral invariants
λ i ( U ) = λ i ( QUQ T ) , i = a 0 · u i = Q a 0 · Q u i , b i = b · u i = Q b · Q u i ,
e i = e · u i = Q e · Q u i , B L , E L e ,
where λ i are the principal stretches and u i are the principal directions of the right-stretch tensor U , i.e.,
U = i = 1 3 λ i u i u i .
Since a 0 , b and e are unit vectors, we have the relations
i = 1 3 a 0 · u i = 1 , i = 1 3 b · u i = 1 , i = 1 3 e · u i = 1 .
In view of the relations in (67), only 11 of the 14 invariants in (65) are independent. We can then express
Φ = Φ ( λ i , i , b i , e i , B L , E L e ) .
We strongly emphasize that the function Φ must satisfy the P-property as described in [16]. The construction of a P-property Φ , where it is independent of the signs of a 0 , b , and e , can be facilitated using the following list of invariants [17]:
λ i , α i = ( a 0 · u i ) 2 , β i = ( b · u i ) 2 , γ i = ( e · u i ) 2 , ζ i = ( a 0 · b ) i b i , κ i = ( a 0 · e ) i e i , ϱ i = ( b · e ) b i e i , η i = ( a 0 · b ) ( a 0 · e ) b i e i , χ i = ( a 0 · b ) ( b · e ) i e i , ι i = ( b · e ) ( a 0 · e ) i b i , B L , E L e .
In (69), only 11 of the 14 spectral invariants are independent [16]. The advantages of spectral invariants over traditional (classical) invariants [22] are given in Appendix A. The expressions for the derivatives of Φ in C , E L e , and B L are almost exactly the same as in [11] (see Section 3 therein), with the only difference in the use of E L e instead E L . We present them here for completeness.
In the case of considering the list of invariants in (69), we have (see [11,18]):
Φ C i i = 1 2 λ i Φ λ i , Φ C i j = 1 ( λ i 2 λ j 2 ) Φ α i Φ α j i j + Φ β i Φ β j b i b j + Φ γ i Φ γ j e i e j + 1 2 ( a 0 · b ) Φ ζ i Φ ζ j ( i b j + j b i ) + 1 2 ( a 0 · e ) Φ κ i Φ κ j ( i e j + j e i ) + 1 2 ( b · e ) Φ ϱ i Φ ϱ j ( b i e j + b j e i ) + 1 2 ( a 0 · b ) ( a 0 · e ) Φ η i Φ η j ( b i e j + b j e i ) + 1 2 ( a 0 · b ) ( b · e ) Φ χ i Φ χ j ( i e j + j e i )
+ 1 2 ( a 0 · e ) ( b · e ) Φ ι i Φ ι j ( i b j + j b i ) , i j ,
where in the above expressions there is no sum in i and j. These expressions can be used in (47) and (57).
On the other hand
Φ B L i = Φ B L b i + 1 B L ( I b b ) T Φ b · u i ,
Φ E L e i = Φ E L e e i + 1 E L e ( I e e ) T Φ e · u i .
The relations (72) and (73) can be used in (53).

4. Final Remarks

In this article, a finite strain-based time-dependent mathematical framework for magnetoelectrically coupled materials is proposed, where a simple constitution is presented, as shown in Equations (49) and (56). Magnetoelectrically activated soft materials such as smart polymers and hydrogels have become increasingly essential multifunctional materials with a myriad of potential applications ranging from drug delivery with the help of tiny robots to developing energy-harvesting systems to meet the expanding demands of clean and green renewable energy. In practical applications, such as robotic devices or energy-harvesting systems, time-dependent dynamical effects must not be discarded to appropriately model a smart system activated by a magnetoelectrically coupled field; hence, the present model will fill the gap in the current literature, in the sense that, to the best of our knowledge, existing finite-strain constitutive equations can only model time-independent effects (see, for example, references [12,20]). We are currently developing the corresponding time-dependent finite-element formulation and will use it to solve practical boundary value problems and to aid the development of specific expressions for the amended energy function Ω via experimental data. Our time-dependent finite-element formulation will aid smart material designs, such as designing robotic devices, etc.

Author Contributions

Conceptualization, M.H.B.M.S., R.B. and M.H.; methodology, M.H.B.M.S., R.B. and M.H.; formal analysis, M.H.B.M.S., R.B. and M.H.; investigation, M.H.B.M.S., R.B. and M.H.; writing—original draft preparation, M.H.B.M.S., R.B. and M.H.; writing—review and editing, M.H.B.M.S., R.B. and M.H. 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.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Physical modelling often requires representations for isotropic functions [22]. A concise and efficient modelling of complex material behaviours requires the representations of constitutive equations to be as compact as possible and, in view of this, irreducible complete representations [22,23] play an important role in facilitating constitutive equations. However, Shariff [16] (see also references therein) showed that, in general, the numbers of traditional (classical) isotropic invariants [22,23,24] in irreducible bases (or minimal integrity bases [22]) were not minimal and, in general, some of them were not independent. Minimal irreducible representations are attractive, in the sense that they reduce modelling complexity. Established classical irreducible isotropic bases are generally not minimal and most of them have no direct physical interpretation; hence, the resulting constitutive equations are too complicated for practical applications (see for example reference [25]). In view of the unclear physical interpretation of classical isotropic invariants, it is not clear how to select an appropriate (or optimum) subset from the corresponding full set (which generally contains numerous elements) of irreducible isotropic bases to represent a physical model. However, using the spectral invariants developed in Shariff [16], not only is our irreducible spectral basis minimal and the spectral invariants independent, but our spectral invariants have a clear physical interpretation, which is attractive, as exemplified below.
For simplicity, we consider a strain energy for a transversely isotropic elastic solid with the preferred direction a 0 in the reference configuration. The classical invariants [22] required to describe the strain energy function are
K 1 = tr C , K 2 = 1 2 ( tr C ) 2 tr C 2 ,
K 3 = det C , K 4 = a 0 · C a 0 , K 5 = a 0 · C 2 a 0 ,
where t r denotes the trace of a tensor. Any effort to construct a reasonable constitutive equation using the invariants in (A1) is hindered due to the constraint that the invariants in (A1) depend explicitly on C and the strain energy function generally depends explicitly on these invariants; hence, it does not allow us to use more general invariants (that cannot be explicitly expressed in terms of C ) which could ease the construction of a superior strain energy function. Furthermore, except for K 3 and K 4 , the invariants in (A1) do not have a direct physical interpretation; hence, they are not experimentally attractive in the sense that they are not suitable for a rigorous experiment to obtain a specific form of strain energy function that requires to vary one invariant while keeping the remaining invariants fixed [26] (denoted by an R-experiment). However, using the spectral invariants λ i and a 0 · u i , Shariff [27] showed that it was possible to construct an R-experiment.
To summarize, it is clear that, for a model that requires a large number of classical invariants to characterize its constitutive equation, considering that most classical invariants do not have a direct physical interpretation and some of them are not independent, it is far more difficult to construct an R-experiment using classical invariants than using spectral invariants. In addition to this, as mentioned above, it is not clear how to select an optimum subset from a corresponding full set of classical-invariant irreducible bases, which generally contain numerous elements, in order to represent a physical model.
We strongly emphasize that spectral formulations are more general, since all classical isotropic invariants [22,23] can be explicitly written in terms of spectral invariants [16].

References

  1. Li, S.; Bai, H.; Shepherd, R.F.; Zhao, H. Bio-inspired Design and Additive Manufacturing of Soft Materials, Machines, Robots, and Haptic Interfaces. Angew. Chem. Int. Ed. 2019, 58, 11182–11204. [Google Scholar] [CrossRef] [PubMed]
  2. Wan, X.; Luo, L.; Liu, Y.; Leng, J. Direct Ink Writing Based 4D Printing of Materials and Their Applications. Adv. Sci. 2020, 7, 2001000. [Google Scholar] [CrossRef] [PubMed]
  3. Amjadi, M.; Pichitpajongkit, A.; Lee, S.; Ryu, S.; Park, I. Highly stretchable and sensitive strain sensor based on silver nanowire-elastomer nanocomposite. ACS Nano 2012, 8, 5154–5163. [Google Scholar] [CrossRef] [PubMed]
  4. Bednarek, S. The giant magnetostriction in ferromagnetic composites within an elastomer matrix. Appl. Phys. A 1999, 68, 63–67. [Google Scholar] [CrossRef]
  5. Hu, W.; Lum, G.Z.; Mastrangeli, M.; Sitti, M. Small-scale soft-bodied robot with multimodal locomotion. Nature 2018, 554, 81–85. [Google Scholar] [CrossRef]
  6. Bastola, A.K.; Hossain, M. A review on magneto-mechanical characterizations of magnetorheological elastomers. Compos. Part B Eng. 2020, 200, 108348. [Google Scholar] [CrossRef]
  7. Eerenstein, W.; Mathur, N.D.; Scott, J.F. Multiferroic and magnetoelectric materials. Nature 2006, 442, 759–765. [Google Scholar] [CrossRef]
  8. Fetisov, Y.K.; Bush, A.A.; Kamentsev, K.E.; Ostashchenko, Y.A.; Srinivasan, G. Ferritepiezoelectric multilayers for magnetic field sensors. IEEE Sensors J. 2006, 6, 935–938. [Google Scholar] [CrossRef]
  9. Liu, L.; Sharma, P. Giant and universal magnetoelectric coupling in soft materials and concomitant ramifications for materials science and biology. Phys. Rev. E 2013, 88, 040601. [Google Scholar] [CrossRef] [Green Version]
  10. Liu, L. An energy formulation of continuum magneto-electro-elasticity with applications. J. Mech. Phys. Solids 2014, 63, 451–480. [Google Scholar] [CrossRef]
  11. Bustamante, R.; Shariff, M.H.B.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]
  12. Shariff, M.H.B.M. Anisotropic stress softening of electromagnetic Mullins materials. Math. Mech. Solids 2023, 28, 154–173. [Google Scholar] [CrossRef]
  13. Pao, Y.H. Electromagnetic Forces in Deformable Continua. Mech. Today 1978, 4, 209–306. [Google Scholar]
  14. Eringen, A.C.; Maugin, G.A. Electrodynamics of Continua; Springer: Berlin/Heidelberg, Germany, 1990. [Google Scholar]
  15. Kovetz, A. Electromagnetic Theory; University Press: Oxford, UK, 2000. [Google Scholar]
  16. Shariff, M.H.B.M. On the Smallest Number of Functions Representing Isotropic Functions of Scalars, Vectors and Tensors. arXiv 2023, arXiv:2207.09617. [Google Scholar] [CrossRef]
  17. Bustamante, R.; Shariff, M.H.B.M. New sets of spectral invariants for electro-elastic bodies with one and two families of fibres. Eur. J. Mech.-A/Solids 2016, 58, 42–53. [Google Scholar] [CrossRef]
  18. Shariff, M.H.B.M. Spectral derivatives in continuum mechanics. Q. J. Mech. Appl. Math. 2017, 70, 479–496. [Google Scholar] [CrossRef]
  19. Ogden, R.W. Non-Linear Elastic Deformations; Dover: New York, NY, USA, 1997. [Google Scholar]
  20. Dorfmann, A.; Ogden, R.W. Magnetoelastic modelling of elastomers. Eur. J. Mech. A/Solids 2003, 22, 497–507. [Google Scholar] [CrossRef]
  21. Dorfmann, A.; Ogden, R.W. Nonlinear electroelasticity. Acta Mech. 2005, 174, 167–183. [Google Scholar] [CrossRef]
  22. Spencer, A.J.M. Theory of invariants. In Continuum Physics; Eringen, A.C., Ed.; Academic Press: New York, NY, USA, 1971; Volume 1, pp. 239–353. [Google Scholar]
  23. Zheng, Q.S. Theory of representation for tensor function. A unified invariant approach to constitutive equations. Appl. Mech. Rev. 1994, 47, 545–587. [Google Scholar] [CrossRef]
  24. Smith, G.F. On isotropic functions of symmetric tensors, skew-symmetric tensors and vectors. Int. J. Eng.Sci. 1971, 9, 899–916. [Google Scholar] [CrossRef]
  25. Merodio, J.; Rajagopal, K.R. On Constitutive Equations For Anisotropic Nonlinearly Viscoelastic Solids. Math. Mech. Solids 2007, 12, 131–147. [Google Scholar] [CrossRef]
  26. Holzapfel, G.A.; Ogden, R.W. On planar biaxial tests for anisotropic nonlinearly elastic solids: A continuum mechanical framework. Math. Mech. Solids 2009, 14, 474–489. [Google Scholar] [CrossRef]
  27. Shariff, M.H.B.M. Nonlinear transversely isotropic elastic solids: An alternative representation. Q. J. Mech. Appl. Math. 2008, 61, 129–149. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Deformation due to the application of electromagnetic fields B and E , boundary displacement, and boundary traction t ^ . B r is the reference (undeformed) configuration, B t is the current configuration, and X and x are, respectively, the position vectors of X in the reference and current configurations, where X represents a generic particle of the solid body.
Figure 1. Deformation due to the application of electromagnetic fields B and E , boundary displacement, and boundary traction t ^ . B r is the reference (undeformed) configuration, B t is the current configuration, and X and x are, respectively, the position vectors of X in the reference and current configurations, where X represents a generic particle of the solid body.
Symmetry 15 00628 g001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shariff, M.H.B.M.; Bustamante, R.; Hossain, M. A Generalised Time-Dependent Mathematical Formulation for Magnetoelectrically Coupled Soft Solids at Finite Strains. Symmetry 2023, 15, 628. https://doi.org/10.3390/sym15030628

AMA Style

Shariff MHBM, Bustamante R, Hossain M. A Generalised Time-Dependent Mathematical Formulation for Magnetoelectrically Coupled Soft Solids at Finite Strains. Symmetry. 2023; 15(3):628. https://doi.org/10.3390/sym15030628

Chicago/Turabian Style

Shariff, Mohd Halim Bin Mohd, Roger Bustamante, and Mokarram Hossain. 2023. "A Generalised Time-Dependent Mathematical Formulation for Magnetoelectrically Coupled Soft Solids at Finite Strains" Symmetry 15, no. 3: 628. https://doi.org/10.3390/sym15030628

APA Style

Shariff, M. H. B. M., Bustamante, R., & Hossain, M. (2023). A Generalised Time-Dependent Mathematical Formulation for Magnetoelectrically Coupled Soft Solids at Finite Strains. Symmetry, 15(3), 628. https://doi.org/10.3390/sym15030628

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