Magnetostrictive Material
Magnetostriction
Magnetostriction describes the change in dimensions of a material due to a change in its magnetization. This phenomenon is a manifestation of magnetoelastic coupling, which is exhibited by all magnetic materials to some extent. The effects related to magnetoelastic coupling are described by various names. The Joule effect describes the change in length due to a change in the magnetization state of the material. This magnetostrictive effect is used in transducers for applications in sonars, acoustic devices, active vibration control, position control, and fuel injection systems.
The inverse effect accounts for the change in magnetization due to mechanical stress in the material. This effect is also known Villari effect. This effect is mostly useful in sensors.
Magnetostriction has a quantum-mechanical origin. The magneto-mechanical coupling takes place at the atomic level due to spin-orbit coupling. From a system level, the material can be assumed to consist of a number of tiny ellipsoidal magnets which rotate due to the torque produced by the externally applied magnetic field. The rotation of these elemental magnets produces a dimensional change leading to free strain in the material. Linear Magnetostriction Model
The magnetostriction has a nonlinear dependence on the magnetic field and the mechanical stress in the material. However, the effect can be modeled using linear coupled constitutive equations if the response of the material consists of small deviations around an operating point (bias point).
It is possible to express the relation between the stress S, strain ε, magnetic field H, and magnetic flux density B in either a stress-magnetization form or strain-magnetization form:
Stress-Magnetization Strain-Magnetization where μ0 is the magnetic permeability of free space; cH and sH are the stiffness and compliance matrices measured at constant magnetic field, respectively; and μrS and μrT are the relative magnetic permeabilities measured at constant strain and constant stress, respectively. The matrices dHT and eHS are called piezomagnetic coupling matrices.
In COMSOL Multiphysics, both constitutive forms can be used; simply select one, and the software makes all necessary transformations. The following equations transform strain-magnetization material data to stress-magnetization data: You find all the necessary material data inputs within the Magnetostrictive Material node under the Solid Mechanics interface, which are added automatically when you add a predefined Magnetostriction multiphysics interface. Such a node can be also added manually under any Solid Mechanics interface similar to all other material model features. The Magnetostrictive Material uses Voigt notation for the anisotropic material data. More details about the data ordering can be found in the Orthotropic and Anisotropic Materials section.
For a crystalline material with tetragonal symmetry, the strain-magnetization form of the constitutive relations is the following:  The following material data corresponds to Terfenol-D at 100 kA/m bias and 30 MPa prestress (Ref. 5):  Nonlinear Magnetostriction
All domains have magnetization of the same magnitude |M| = Ms, but the magnetization can have different orientations characterized by the corresponding direction vector m = M/Ms for each domain. The applied magnetic field changes the domain orientation. In this section, the term domain refers to a small part of magnetic material. This is typical for micromagnetics literature, and it should not be mistaken with the concept of domain as part of the model geometry, the latter is often used in COMSOL Multiphysics documentation.
For a single crystal with cubic symmetry, the magnetostrictive strain tensor can be written as the following quadratic form:
where . Note that the magnetostrictive strain is represented by a deviatoric tensor, that is, trme) = 0. This is because the deformation is related to the magnetic domain rotation, and such process should not change the material volume.
The strain in any direction given by the directional cosines βi can be written as Using Equation 3-67, one gets
When both magnetization and measurement direction are parallel to the same crystal direction , one has m1 = β1 = 1 and all other components are zero, so that  In many applications, such alignment of the domains is achieved by applying a compressive prestress. Thus, the maximum usable magnetostriction is achieved via a 90-degree rotation of the domains For an isotropic material, λ100 = λ111 = λs, and Equation 3-67 becomes  In COMSOL Multiphysics, this description of the magnetostriction is modeled using the following equation for the magnetostrictive strain:
The strain field is deviatoric, and Equation 3-69 exhibits the same properties as Equation 3-67 at saturation, that is, when |M| → Ms. Equation 3-68 is replaced by Note that the strain vanishes when |M| → 0, which makes the model applicable in the whole range from full demagnetization to saturation.
The stress in the magnetostrictive material is modeled as For isotropic materials, the stiffness matrix cH can be represented in terms of two parameters, for example, using the Young’s modulus and Poisson’s ration. Cubic materials possess only three independent components: c11, c12, and c44.
Using Equation 3-70, one can derive a linear response around a given bias state characterized by a premagnetization vector M0. Thus, where M1 is a perturbation, and If one assumes a unidirectional state, for example it will further simplify into Magnetization
where L is the Langevin function with χm being the magnetic susceptibility in the initial linear region.
Other possible choices of the L function are a hyperbolic tangent and a linear function The latter option will make it possible to find an explicit expression for the magnetization. However, such model does not have a proper saturation behavior, and thus it should be used only in the operating range far from saturation. Both the Langevin function and hyperbolic tangent models requires the magnetization vector components to be treated as extra dependent variables.
For cubic crystals, the effective field in the material is given by
where H is the applied magnetic field. The second term in Equation 3-72 represents the mechanical stress contribution to the effective magnetic field, and thus to the material magnetization, which is called the Villari effect. The deviatoric stress tensor is related to the strain as For isotropic materials, the effective magnetic field expression simplifies into COMSOL Multiphysics solves for the magnetic vector potential A whose curl yields the vector B-field. The H-field is then obtained as a function of the B-field and magnetization.
Hystersis modeling Compared to Equation 3-72 and Equation 3-73, the effective magnetic field Heff gets one more term αMan, where α is the interdomain coupling parameter.
The change in the total magnetization caused by the change on the effective magnetic field is represented as
where cr is the reversibility parameter, and kp is the pining parameter.
Equation 3-74 can be solved using either a time dependent analysis or a stationary parametric sweep. For more details, see The Jiles-Atherton Hysteresis Model in the AC/DC Module User’s Guide.