Recovered PDF 17204 PDF
Recovered PDF 17204 PDF
Recovered PDF 17204 PDF
www.actamat-journals.com
Abstract
A computational model is proposed to predict the stability of magnetic domain structures and their temporal evolution in giant
magnetostrictive materials by combining a micromagnetic model with the phase-eld microelasticity theory of Khachaturyan. The
model includes all the important energetic contributions, including the magnetocrystalline anisotropy energy, exchange energy, magnetostatic energy, external eld energy, and elastic energy. While the elastic energy of an arbitrary magnetic domain structure is
obtained analytically in Fourier space, the LandauLiftshitzGilbert equation is solved using the ecient GaussSeidel projection
method. Both Fe81.3Ga18.7 and Terfenol-D are considered as examples. The eects of elastic energy and magnetostatic energy on
domain structures are studied. The magnetostriction and associated domain structure evolution under an applied eld are modeled
under dierent pre-stress conditions. It is shown that a compressive pre-stress can eciently increase the overall magnetostrictive
eect. The results are compared with existing experiment measurements and observations.
2005 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.
Keywords: Phase-eld models; Simulation; Domain structure; Magnetostriction
1. Introduction
Magnetostriction is a phenomenon in which a material changes its physical dimensions when there is a
change in its magnetization. The change in magnetization can be induced by the application of a magnetic
eld or a stress eld. Magnetostriction is observed, to
diering degrees, in all ferromagnetic materials. While
typical magnetostrictive strains of magnetic materials
are of the order of 105106, some compounds containing rare earth elements have magnetostriction as
high as 103. Among these, Terfenol-D (TbxDy1 xFe2,
x = 0.3) has been extensively investigated for use in sensor and actuator devices because of its high magnetostriction to anisotropy ratio [1,2]. It has been recently
*
Corresponding author. Tel.: +1 814 865 0389; fax: +1 814 865 2917.
E-mail address: jzz108@psu.edu (J.X. Zhang).
1359-6454/$30.00 2005 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.
doi:10.1016/j.actamat.2005.03.002
2846
essentially no simulation of macroscopically inhomogeneous domain patterns and their evolution under external elds in giant magnetostrictive materials.
One type of approach assumed that the magnetizations of domains are constrained by the anisotropy
energy to lie along the easy magnetic axes. The ferromagnetic constitutive behavior was then described by
the distribution of a xed quantity of available magnetic
domains within materials among all the possible directions. The model assumed that each domain is magnetized independently without being aected by the
magnetization states of neighboring domains. This type
of approach has been used in many works to study the
magnetostriction under an external magnetic eld [20
24]. But this model is not suitable for describing a material with low magnetocrystalline anisotropy or the
system states close to saturation when the magnetic eld
is applied along directions far away from those of easy
axes. Zhu et al. [25] studied the inuence of external
stress on the cohesive eld in magnetic thin lms. In
their micromagnetic model, without solving the total
compatible strain, they use the stress-free magnetostrictive strain caused by the magnetization instead; thus the
intrinsic stress resulting from the elastic incompatibility
of magnetostrictive strains is not taken into account.
The intrinsic stress is critically important for materials
with large magnetostriction [26]. Fabian and Heider
[27] have calculated the magnetostrictive self-energies
of a magnetic particle imbedded in a non-magnetic matrix using the nite element method and the continuum
theory of defects. Recently, Shu et al. [26,28] developed
a modied boundary integral formalism to calculate the
intrinsic stress induced by incompatible magnetostrictive
strains, but their method is conned to some simple twodimensional (2D) cases due to the complexity of solving
the elastic equilibrium equations.
The main purpose of this paper is to describe a
computational approach to model the stability and evolution of magnetic domain structures in giant magnetostrictive materials. It combines the micromagnetic model
for magnetic domain evolution [29] and the phase-eld
microelasticity theory [30,31] for the elastic solutions
in a magnetostrictive material with arbitrary distributions of magnetic domains. The phase-eld microelasticity theory has been extensively used in computer
simulations of structural phase transformations and
microstructure evolution in both bulk systems, see, for
example, Refs. [3234], and for more extensive references [35], and thin lms [3638]. The model is able to
predict the detailed domain structure and their evolution
under an applied eld without a priori assumptions on
domain morphologies. The main limitation of the present model as well as essentially all prior phase-eld models for solid state phase transformations is the
assumption of periodic boundary conditions, and hence
the results are only valid when the simulation system size
oM
ca
c0 M Heff 0 M M Heff ;
ot
Ms
where Ms is the saturation magnetization, c0 is the gyromagnetic ratio, a is the damping constant, and He is the
eective magnetic eld, which can be represented as a
variational derivative of the total free energy of the system with respect to magnetization
1 oE
;
2
Heff
l0 oM
where l0 is the permeability of vacuum, the total free energy of magnetostricitve materials is given by
E Eanis Eexch Ems Eexternal Eelastic ;
where Eanis, Eexch, Ems, Eexternal, Eelastic are the magnetocrystalline anisotropy energy, exchange energy, magnetostatic energy, external eld energy, and elastic
energy, respectively.
The magnetocrystalline anisotropy energy of a cubic
crystal is
Z
Eanis
K 1 m21 m22 m21 m23 m22 m23 K 2 m21 m22 m23 dV ;
4
where mi are the components of the unit magnetization
vector, m = M/Ms. K1 and K2 are the anisotropy
constants.
The exchange energy is determined solely by the spatial variation of the magnetization orientation and can
be written as
Z
Eexch A m21;1 m21;2 m21;3 m22;1 m22;2 m22;3
m23;1 m23;2 m23;3 dV ;
for example, mi,j = omi/oxj, where xj is the jth component of position vector in the Cartesian coordinates.
The magnetostatic energy of a system can be written
as
Z
1
Ems l0 M s Hd m dV ;
6
2
where Hd is the stray eld that is determined by the
long-range interaction among the magnetic moments
in the system. The stray eld Hd is governed by [29]
H d1;1 H d2;2 H d3;3 M s m1;1 m2;2 m3;3 ;
The solution of the potential for Eq. (9) is given in Fourier space by,
M s m1 nn1 m2 nn2 m3 nn3
;
10
n21 n22 n23
p
where i 1; ni are the coordinates in Fourier space,
and /(n), mi(n) are the Fourier transforms of / and mi,
respectively. The value of / in real space can be obtained through an inverse Fourier transform of /(n).
The value of Hd can then be calculated by Eq. (8).
Note that the origin of the Fourier space
(n1 = n2 = n3 = 0) is a singularity point in Eq. (10) and
its contribution corresponds to the demagnetization
eld caused by the average magnetization of the system.
This can easily be seen by writing the local magnetization eld as
dMr;
Mr M
11
/n i
12
2847
15
where eij is the total strain. The corresponding elastic energy can be expressed as
Z
Z
1
1
cijkl eij ekl dV
cijkl eij e0ij ekl e0kl dV ;
Eelastic
2
2
16
where cijkl is the elastic stiness tensor. The summation convention for the repeated indices is employed
and i, j, k, l = 1, 2, 3. For a cubic material with its
three independent elastic constant c11, c12 and c44 in
the Voigts notation, the elastic energy can be rewritten as
Eelastic
Z
1
c11 e211 e222 e233
2
2
c12 e11 e011 e22 e022 e22 e022 e33 e033
e11 e011 e33 e033 2c44 e12 e012
0 2
0 2
e23 e23 e13 e13 dV
17
2848
Z
18
1
c11 e211 e222 e233
2
Eelastic2
Z ("
2c44
m21 m22
3
k111
2
m22 m23
2
K 31
19
2 #
3
c11 c12 k100
2
)
m21 m23
dV ;
21
Following Khachaturyans theory [31], the total
strain eij(r) may be represented as the sum of homogeneous and heterogeneous strains:
eij r eij gij r:
22
The homogeneous strain is dened in such a way so that
Z
23
gij r dV 0:
The homogeneous strain represents the macroscopic
shape change of a system generated due to the formation
of a domain structure. The heterogeneous strain does
not aect the macroscopic shape of a system.
The equilibrium heterogeneous strain satises the
mechanical equilibrium condition given by the Euler
equation with respect to the elastic displacement
24
25
26
K 33
Dn l2 k 2l vn6
lv2k 2l vn2 n21 n22 n21 n23 n22 n23
v2 3k 3l vn21 n22 n23 ;
29
Z
rij;j 0;
K 32
20
3
k100 c11 c12 e11 m21 e22 m22 e33 m23
2
6k111 c44 e12 m1 m2 e23 m2 m3 e13 m1 m3 dV :
Eelastic3
where X i icijkl e0kl nnj ; ui n and e0kl n are the Fourier transforms of ui and e0kl , respectively, Nij(n) are
cofactors of a 3 3 matrix K[n],
2
3
K 11 K 12 K 13
6
7
28
Kn 4 K 21 K 22 K 23 5
27
l c44 ;
and n2 ni ni :
30
2
Z
V
cijkleijekl cijkleij e0kl dV
2
Z
1
cijkl gij r e0ij gkl r e0kl dV ;
31
2
where V is the total volume of the system.
33
where
and mi mj represent the volume average of
m2i and mi mj over a system containing a domain structure, respectively.
When a system is subject to a homogeneous applied
stress raij , the total potential energy is given by the sum
of elastic energy of the system and the potential of the
mechanical loading,
Ep Eelastic V raijeij :
oeij
oeij
oeij
Z
Vcijklekl cijkl e0kl dV V raij 0:
Therefore
Z
1
ekl
e0kl dV sijkl raij ;
V
The total strain can be calculated by Eq. (22). Consequently, the elastic energy can be obtained for a magnetostrictive material with an arbitrary magnetic domain
structure by Eq. (16).
39
1 a2
s:
t
c0 M s
The equation becomes,
om
m heff am m heff :
os
35
2849
40
1
1 oEexch
1 oE0
Ms
l0 oM
l0 oM
0
1
1
oE
A Dm
Ms
l0 oM
heff
36
41
A Dm hm;
37
i 1; 2; 3:
43
m3
mn3
g1 m2
g2 m2
44
2850
Step 2:
0 1 0
1
m1 aDsA Dm
m1
1 h1 m
B C B
C
@ m2 A @ m2 aDsA Dm
2 h2 m A:
m
3
m3
aDsA
Dm
3
Step 3:
0 n1 1
0 1
m1
m1
1 B C
B n1 C
@ m2 A @ m2 A:
jm j
m
mn1
3
3
45
h3 m
46
gni n
mni n Ds~
hi mn
;
2
2
1 n1 n3 n23 A Ds
47
the magnetostatic energy. Mechanical clamped boundary conditions were used, which means that the overall
dimensions of the system do not change. At the beginning, the random initial conguration evolves to domains as shown in Fig. 1(a). All the six dierent kinds
of orientation domains determined by the cubic anisotropy energy, including the domains with magnetization
along the positive x1 direction ([1 0 0] domain), negative
x1 direction (1 0 0 domain), positive x2 direction ([0 1 0]
domain), negative x2 direction (0 1 0 domain), positive
x3 direction ([0 0 1] domain), and negative x3 direction
(0 0 1 domain), are present with almost equal fractions.
Fig. 1(b) illustrates the domain structure after 40,000
steps of evolution, indicating domain growth and coarsening. The domains are separated by 180 or 90 domain
walls. A schematic diagram showing the 180 domain
and 90 domain walls separating the domains is given
in Fig. 2. The 90 domain walls tend to align along
the [1 1 0] or 1 1 0 directions separating [1 0 0] domain
(or 1 0 0 domain) and [0 1 0] domain (or 0 1 0 domain).
The alignment becomes increasingly strong at the later
stage of evolution as shown in Fig. 1(c) and (d). Since
in this case, the magnetostatic energy is not included
and the exchange energy is isotropic, the domain wall
alignment must be entirely due to the anisotropic elastic
interactions. For 90 domain walls, the condition of
elastic compatibility between the associated domain
pairs requires that the tangential components of the
(a)
(b)
Fig. 2. Schematic diagram of the 180 domain wall and 90 domain
wall separating the domains: (a) 180 domain wall; (b) 90 domain
wall.
2851
Fig. 3. Simulated domain structure in the absence of both magnetostatic and elastic energies. Black = [1 0 0] and 1 0 0 domains, light
gray = [0 1 0] and 0 1 0 domains, and dark gray = [0 0 1] and 0 0 1
domains.
2852
It should be noted that some other material parameters of FeGa alloy beside magnetostriction constants
were also found to depend on composition and temperature, such as saturation magnetization and cubic
anisotropy constant. They should also aect the domain
structure of FeGa alloy, which was not studied in this
work. Thus, complicated behaviors of domain structures
in FeGa are expected depending on composition and
temperature.
Fig. 5 shows the temporal evolution of domain structure when the magnetostatic energy and elastic energy
are both taken into account. In this case, the domain
wall orientation is constrained by both the elastic compatibility and the magnetic compatibility. Two adjacent
domains are magnetically compatible if the normal components of the magnetization vectors of both domains
with respect to their common domain wall are equal.
Therefore, only head to tail domain walls survive. Since
the discontinuities in the strain do not occur for 180 domain walls, the 180 domain walls are only constrained
by the magnetic compatibility. It is shown that the 180
domain walls are found to be parallel to the magnetization direction of associated domain pair (Fig. 5(d)). Finally, a twin-like domain structure was obtained as
Fig. 6. Simulated domain structure with dierent applied pre-stress (along x1 axis): (a) 5 MPa (61% domain with magnetization along x1 axis); (b) 0
MPa (49%); (c) 5 MPa (32%). Dark gray = [1 0 0] domains, black = 1 0 0 domains, light gray = 0 1 0 domains, and white = [0 1 0] domains.
300
0.6
0.4
0.2
0.0
(a)
2
4
6
Applied Field (kA/m)
(b)
300
5 MPa
0 MPa
-5 MPa
250
200
150
100
50
0
2
4
6
8
Applied Field (kA/m)
Magnetostriction (x10-6)
M/Ms
0.8
Magnetostriction (x10-6)
5 MPa
0 MPa
-5 MPa
1.0
(c)
5 MPa
0 MPa
-5 MPa
250
200
150
100
50
0
0.0
0.2
0.4
0.6
M/Ms
0.8
1.0
Fig. 7. The eld dependence of magnetization and magnetostriction for the giant magnetostrictive material with dierent pre-stress applied: (a) M/Ms vs.
applied eld curves; (b) magnetostriction vs. applied eld curves; (c) magnetostriction vs. M/Ms curves.
2853
Fig. 8. The domain structure evolution under applied eld along [1 0 0] (without pre-stress): (a) 0 kA/m; (b) 0.80 kA/m; (c) 4.38 kA/m; (d) 7.96 kA/m.
Dark gray = [1 0 0] domains, black = 1 0 0 domains, light gray = 0 1 0 domains, and white = [0 1 0] domains.
2854
a compressive stress will favor domains with magnetization perpendicular to the x1 axis. Fig. 6 shows the domain structures of samples with dierent applied prestresses. As expected, while a tensile stress increases
the fraction of domains along x1 axis, a compressive
stress applied along x1 axis decreases it.
The magnetostriction and associated domain structure evolution under an applied magnetic eld for the
samples with dierent applied pre-stress were shown in
Figs. 7 and 8. In prior experiments [35], to obtain large
magnetostriction at a relatively low magnetic eld, the
samples were usually prepared in the shapes of a thin
disk or a rod to decrease the demagnetization eect.
To approximate this case, in our simulation we use the
demagnetizing factor of innite thin sheet (Nx =
Ny = 0, Nz = 1). A magnetostriction vs. M/Ms curve for
the magnetization process was given in Fig. 7(c). The
initial at stage of the curve show that the initial increase of magnetization is accompanied by a small magnetostrictive change. As Fig. (8) shows, such an initial
magnetization stage mainly consist of 180 domain wall
movement (such as region A in Fig. 8(a) and (b)), which
does not produce any magnetostrictive change in dimensions. 90 domain wall movement also occurred by
switching from 1 0 0 to [0 1 0] (or 0
1 0) domains or
from [0 1 0] (or 0
1 0) to [1 0 0] domains (such as region
B in Fig. 8(a) and (b)). As magnetostriction vs. M/Ms
curve measures the average behavior of the system, the
eect of such two distinct 90 domain wall motion is
equal to a 180 domain switching from
1 0 0 to [1 0 0]
domain, which does not produce magnetostrictive
change. With further magnetization, the fraction of the
1 0 0 domain deceases, and the associated 180 domain
switching decreases too. As shown in the magnetostriction vs. M/Ms curve, the slope increases with magnetization. Finally, a single domain with saturated
magnetization and magnetostriction along the applied
eld was obtained as shown in Fig. 8(d). It is shown that
a pre-stress can indeed have a signicant eect on the
magnetostriction. A larger saturation magnetostriction
was obtained for systems with a compressive pre-stress,
while a smaller magnetostriction was obtained when a
tensile pre-stress was applied.
c44 = 4.87 1010 N/m2. Because of the extensive computational time required for 3D simulations, in this paper
we only focused on the prediction of 3D domain structures without an applied eld. An example of 3D domain structure is shown in Fig. 9. Domains with
magnetization along the 8 possible 111 easy directions
were observed. The obtained domain wall orientations
of non-180 walls, which constrained by the elastic compatibility, agree with those predicted by theory of James
and Kinderlehrer [17,18] and with prior experimental
observations [6,11]. It should be emphasized that these
domain wall orientations were predicted without a priori
assumptions. The orientations of 180 domain wall are
not exactly along those planes determined by the magnetic compatibility, which comes from the competition
between magnetostatic and exchange energy, since the
exchange energy prefers a at domain wall with minimization of wall. Extensive studies on the magnetic behavior of Terfenol-D such as the magnetic process with
applied magnetic eld or stress will be carried out using
extensive 3D simulations.
Acknowledgment
The authors are grateful for the nancial support by
the National Science Foundation under the Grant No.
DMR 01-22638.
References
[1] Clark AE. Magnetostrictive rare earth-Fe2 compounds. In:
Wohlfarth EP, editor. Ferromagnetic materials. Amsterdam: North-Holland; 1980.
[2] Jiles DC. J Phys D: Appl Phys 1994;27:1.
[3] Clark AE, Restor JB, Wun-Fogle M, Lograsso TA, Schlagel
DL. IEEE Trans Mag 2000;36:3238.
[4] Clark AE, Hathaway KB, Wun-Fogle M, Restor JB, Lograsso
TA, Keppens VM, et al. J Appl Phys 2003;93:8621.
[5] Raque S, Cullen JR, Wuttig M, Cui J. J Appl Phys 2004;95:6939.
[6] Engdahl G. Handbook of giant magnetostrictive material. New
York, NY: Academic Press; 2000.
[7] Holden AP, Lord DG, Grundy PJ. J Appl Phys 1996;79:4650.
[8] Dooley J, Graef MD, McHenry ME. J Appl Phys 1998;83:6837.
2855
[9] Holden AP, Lord DG, Grundy PJ. J Appl Phys 1996;79:6070.
[10] Lord DG, Holden AP, Grundy PJ. J Appl Phys 1997;81:5728.
[11] Schmidt J, Tickle R, Skidmore GD, Merton C, James RD,
Dahlberg ED. J Magn Magn Mater 1998;190:98.
[12] Shull RD, Quandt E, Shapiro AT, Glasmachers S, Wuttig M. J
Appl Phys 2004;95:6948.
[13] Busbridge SC, Piercy AR. J Appl Phys 1993;73:5354.
[14] Wang BW, Busbridge SC, Li YX, Wu GH, Piercy AR. J Magn
Magn Mater 2000;218:198.
[15] Wang BW, Busbridge SC, Guo ZJ, Zhang ZD. J Appl Phys
2003;93:8489.
[16] Brown WF. Magnetoelastic interactions. Berlin: Springer; 1966.
[17] James RD, Kinderlehrer D. Philos Mag B 1993;68:237.
[18] James RD, Kinderlehrer D. J Appl Phys 1994;76:7012.
[19] DeSimone A, James RD. J Mech Phys Solids 2002;50:283.
[20] Jiles DC, Thoelke JB. J Magn Magn Mater 1994;134:143.
[21] DeSimone A, James RD. J Appl Phys 1997;81:5706.
[22] Zhao XG, Lord DG. J Magn Magn Mater 1999;195:699.
[23] Armstrong WD. J Appl Phys 2002;91:2202.
[24] Park WJ, Son DR, Lee ZH. J Magn Magn Mater 2002;248:
223.
[25] Zhu B, Lo CCH, Lee SJ, Jiles DC. J Appl Phys 2001;89:7009.
[26] Shu YC, Lin MP, Wu KC. Mech Mater 2004;36:975.
[27] Fabian K, Heider F. Geophys Res Lett 1996;23:2839.
[28] Shu YC, Yen JH. Mater Res Soc Symp Proc 2004;785:D13.4.1.
[29] Hubert A, Schafer R. Magnetic domains: the analysis of magnetic
microstructure. Berlin: Springer; 1998.
[30] Khachaturyan AG, Shatalov GA. Sov Phys JETP 1969;29:557.
[31] Khachaturyan AG. Theory of structural transfomations in
solid. New York, NY: Wiley; 1983.
[32] Chen LQ, Wang YZ, Khachaturyan AG. Philos Mag Lett
1992;65:15.
[33] Wang Y, Khachaturyan AG. Acta Mater 1997;45:759.
[34] Jin YM, Artemev A, Khachaturyan AG. Acta Mater
2001;49:2309.
[35] Chen LQ, Hu SY. Phase-eld method applied to strain-dominated
microstructure evolution during solid-state phase transformations. In: Raabe D, Barlat F, Roters F, Chen LQ, editors.
Continuum scale simulation of engineering materials fundamentals microstructures process applications. New York,
NY: John Wiley; 2004.
[36] Li YL, Hu SY, Liu ZK, Chen LQ. Appl Phys Lett 2001;78:3878.
[37] Li YL, Hu SY, Liu ZK, Chen LQ. Acta Mater 2002;50:395.
[38] Seol DJ, Hu SY, Li YL, Shen J, Oh KH, Chen LQ. Acta Mater
2003;51:5173.
[39] Wang XP, Garcia-Cervera CJ, E WN. J Comput Phys 2001;171:
357.
[40] Hartmann U, Mende HH. J Appl Phys 1986;59:4123.