# An analytical model for the magnetostriction strain of ferromagnetic materials subjected to multiaxial stress

Laurent Daniel

► To cite this version:

Laurent Daniel. An analytical model for the magnetostriction strain of ferromagnetic materials subjected to multiaxial stress. European Physical Journal: Applied Physics, 2018, 83 (3), pp.30904. 10.1051/epjap/2018180079. hal-01943237

HAL Id: hal-01943237

https://hal.science/hal-01943237

Submitted on 16 Jul 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# An analytical model for the magnetostriction strain of ferromagnetic materials subjected to multiaxial stress

Laurent Daniel*

GeePs – Group of Electrical Engineering – Paris, UMR CNRS 8507, CentraleSupélec, Univ. Paris-Sud, Université Paris-Saclay, Sorbonne Université, 3 rue Joliot-Curie, Plateau de Moulon, 91192 Gif-sur-Yvette Cedex, France

Received: 2 February 2018 / Received in final form: 23 May 2018 / Accepted: 26 June 2018

**Abstract.** Magnetostriction is the magnetisation-induced strain in ferromagnetic materials. It highly depends on mechanical stress. Stress state in electromagnetic devices is usually multiaxial and its effect on magnetostrictive properties is not easily predicted. In this paper, an original three-parameter analytical model for the stress-dependent magnetostriction strain of ferromagnetic materials is proposed. It is based on a simplified energetic description of magneto-elastic behaviour. It follows a similar method previously adopted for the description of the effect of stress on the magnetic permeability of magnetic materials. It is applied for the first time to the magnetostriction behaviour and results in a simple formula to express the effect of multiaxial magneto-mechanical loadings on the magnetostriction strain. The approach also naturally includes the description of the so-called $\Delta E$ effect. The analytical formula is derived in the paper. It shows very satisfying agreement with experimental results on iron-cobalt alloy and pure iron specimen.

## 1 Introduction

Magnetic and mechanical behaviours are strongly coupled [1]. Magnetisation is very sensitive to the application of stress, leading to significant effects on the performance of electromagnetic devices [2–4]. Conversely, the magnetisation process is associated with a mechanical deformation called magnetostriction. Magnetostriction can be used for actuation purposes, for instance using giant magnetostrictive materials [5–7]. It is also one of the origins of the noise emitted by electromagnetic devices such as transformers [8–10]. Magnetostriction itself is known to be sensitive to stress [11–14]. Although very significant, this effect of stress on the magnetostriction strain is rarely taken into account in the modelling of electromagnetic devices or limited to uniaxial configurations where the stress is uniaxial and applied in the direction parallel to the applied magnetic field [15–17]. However, stress state in electromagnetic devices is usually multiaxial, making the 1D approaches irrelevant. Recently, 3D magneto-mechanical models for the prediction of the magnetostriction strain have been proposed. Some of them are based on the definition of the magneto-elastic free energy of the material [18–20]. These approaches can handle full 3D configurations, but the number of parameters required for the modelling can be high and their identification can be tedious. Some others are based on a multiscale description of magneto-elastic couplings [21–27]. These models provides a very useful insight into magneto-elastic coupling effects, but they are usually too complex to be directly implemented into numerical analysis tools. In order to

overcome this issue, simplified approaches have been derived from these full multiscale approaches [28]. These simplifications lead to a loss of information, notably regarding internal heterogeneity effects, but can be implemented into structural analysis tools for electromagnetic devices [29–31]. In this paper, the simplification procedures are extended further so as to obtain a fully analytical model for the magnetostriction strain as a function of applied magnetic field and stress. It follows a similar method previously adopted for the description of the effect of stress on the magnetic permeability of magnetic materials [32]. In Section 2, the simplified 3D magneto-elastic model [28] is recalled. It is shown in Section 3 how further simplifications can provide a fully analytical definition for the stress-dependent magnetostriction strain. The magnetostriction strain under particular loadings is then detailed (Sect. 4) and a parameter identification method is proposed (Sect. 5). Section 6 is finally dedicated to an illustration of the model prediction and to the comparison to experimental results.

## 2 Simplified magneto-elastic model

A magneto-elastic constitutive law can be derived from the description of a magnetic material as a set of magnetic domains with known magnetisation ($M_s$) and random orientation [28]. This approach, detailed in reference [28] is recalled hereafter. The local free energy (Eq. (1)) of a magnetic domain $k$ is expressed as the sum of three contributions:

$$W_k = W_k^{\text{mag}} + W_k^{\text{el}} + W_k^{\text{an}}. \quad (1)$$

<footnote>* e-mail: laurent.daniel@centralesupelec.fr</footnote>
The Zeeman energy $W_k^{\text{mag}}$ (Eq. (2)) introduces the effect of the applied magnetic field on the equilibrium state. $\mu_0$ is the vacuum permeability. $\mathbf{H}_k$ and $\mathbf{M}_k$ are the magnetic field and magnetisation in the magnetic domain.

$$W_k^{\text{mag}} = -\mu_0 \mathbf{H}_k \cdot \mathbf{M}_k. \qquad (2)$$

The elastic energy $W_k^{\text{el}}$ (Eq. (3)) introduces the effect of stress on the magnetic equilibrium. $\sigma$ is the applied stress and $\varepsilon_k^\mu$ is the magnetostriction strain in the magnetic domain. Like the other energy terms, this energy is defined except for a constant. The derivation of equation (3) from the classical definition of the elastic energy is detailed in reference [24].

$$W_k^{\text{el}} = -\sigma : \varepsilon_k^\mu. \qquad (3)$$

Macroscopic anisotropy can be described through an anisotropy energy term $W_k^{\text{an}}$ added to the free energy. Equation (4) gives this additional term for a uniaxial anisotropy along direction $\mathbf{v}$, with $K$ the anisotropy constant. If we assume macroscopic isotropy, this term vanishes:

$$W_k^{\text{an}} = K(\mathbf{u}_k \cdot \mathbf{v})^2. \qquad (4)$$

Such an approach, very close to Armstrong model [22], was proposed in reference [29] in the 2D case. For a three-dimensional configuration and considering isotropic and isochoric magnetostriction [28], the following definitions for the local magnetisation $\mathbf{M}_k$ and magnetostriction strain $\varepsilon_k^\mu$ are used:

$$M_k = M_s \mathbf{u}_k, \qquad (5)$$

$$\varepsilon_k^\mu = \lambda_s \left( \frac{3}{2} \mathbf{u}_k \otimes \mathbf{u}_k - \frac{1}{2} \mathbf{I} \right). \qquad (6)$$

$M_s$ is the saturation magnetisation of the material, $\mathbf{u}_k$ is the orientation of the magnetisation in the domain $k$, $\lambda_s$ is the saturation magnetostriction constant and $\mathbf{I}$ the second-order identity tensor. Note that equation (6) corresponds to isochoric magnetostriction (tr($\varepsilon_k^\mu$) = 0), so that volume magnetostriction, occurring for high magnetic field levels, is neglected.

The magneto-elastic behaviour is obtained by defining the volume fraction $f_k$ of a domain with orientation $\mathbf{u}_k$ through the use of a Boltzmann probability function [21]:

$$f_k = \frac{\exp(-A_s W_k)}{\int_k \exp(-A_s W_k)}, \qquad (7)$$

where $A_s$ is a material parameter linked to the initial anhysteretic susceptibility $\chi^0$ [24]:

$$A_s = \frac{3\chi^0}{\mu_0 M_s^2}. \qquad (8)$$

Once the probability $f_k$ is defined, the macroscopic magnetisation $\mathbf{M}$ and magnetostriction $\varepsilon_k^\mu$ are obtained,

thanks to an averaging operation over all possible directions:

$$\mathbf{M} = \langle \mathbf{M}_k \rangle = \int_k f_k \mathbf{M}_k, \qquad (9)$$

$$\varepsilon^\mu = \langle \varepsilon_k^\mu \rangle = \int_k f_k \varepsilon_k^\mu. \qquad (10)$$

This integration step can be performed numerically using a discretisation of possible orientations $\mathbf{u}_k$ [25].

Although simplified compared to the full multiscale model [25], this approach is not analytical. Integration operations are required in equations (7), (9) and (10). It can be made analytical by considering a further simplified configuration with a limited number of domains. This is expected to aggravate the limitations of the model proposed in reference [28], particularly in the case of complex anisotropies or high level of heterogeneity. However, such a strategy was already successfully applied in reference [33] for the definition of an equivalent stress for magnetic behaviour, in reference [34] to describe the effect of plasticity on the magnetic behaviour, or in reference [32] to derive an analytical expression for the stress-dependent magnetic permeability of ferromagnetic materials. It is developed here to derive an original analytical expression for the stress-dependent magnetostriction strain.

## 3 Stress-dependent magnetostriction strain

We consider a homogeneous isotropic magnetic material subjected to a magnetic field $\mathbf{H}$ in the direction $\mathbf{x}$ ($\mathbf{H} = H\mathbf{x}$) in an orthonormal coordinate system ($O,x,y,z$). This material is simultaneously subjected to a multiaxial stress state $\sigma$ given by equation (11):

$$\sigma = \begin{pmatrix} \sigma_{xx} & \sigma_{xy} & \sigma_{xz} \\ \sigma_{xy} & \sigma_{yy} & \sigma_{yz} \\ \sigma_{xz} & \sigma_{yz} & \sigma_{zz} \end{pmatrix}_{xyz}. \qquad (11)$$

Any relative orientation between the applied magnetic field and the principal stress can be described. Following the approach proposed in reference [32], the material equilibrium state can be defined using a simplified energy description. The material is assumed to be divided into six domains noted $k$ ($k$= {1, 2, 3, 4, 5, 6}) with magnetisation oriented along $\mathbf{u}_k$ ($\mathbf{u}_k$= { $\mathbf{x}$, $-\mathbf{x}$, $\mathbf{y}$, $-\mathbf{y}$, $\mathbf{z}$, $-\mathbf{z}$}). Each domain $k$ is characterised by its free energy (Eq. (12)) obtained after simplification of equation (1). For strongly anisotropic media, an anisotropy term (e.g. Eq. (4)) should be kept in the free energy (Eq. (12)), to the price of a complexification of the analytical expressions obtained.

$$W_k = - \mu_o H M_s \mathbf{x} \cdot \mathbf{u}_k - \sigma : \varepsilon_k^\mu. \qquad (12)$$

Under these assumptions, the free energy for each domain $k$ can be explicitly written. The volume fraction $f_k$ of each domain (Eq. (7)) is then estimated from equation (13)
using a discrete summation. The relative weight for a given orientation k is directly related to the relative energy of this orientation compared to the other orientations.

$$f_k = \frac{\exp(-A_s W_k)}{\sum_{k=1}^{6} \exp(-A_s W_k)}, \quad (13)$$

which leads to

$$\left\{ \begin{array}{l} f_1 = C A_x \exp(\kappa H) \\ f_2 = C A_x \exp(-\kappa H) \\ f_3 = f_4 = C A_y \\ f_5 = f_6 = C A_z \end{array} \right. \quad (14)$$

using the following notations:

$$C = \frac{\exp\left(-\frac{\alpha}{3}\operatorname{tr}(\sigma)\right)}{\sum_{k=1}^{6} \exp(-A_s W_k)}, \quad (15)$$

$$\alpha = \frac{3}{2} A_s \lambda_s, \quad (16)$$

$$\kappa = \mu_o A_s M_s, \quad (17)$$

$$A_i = \exp(\alpha \sigma_{ii}), \quad i = \{x,y,z\}. \quad (18)$$

As shown in reference [32], the material magnetisation M is obtained, thanks to the discrete summation (Eq. (19)), which leads to the analytical definition (Eq. (20)) for the magnetisation:

$$M = \sum_{k=1}^{6} f_k M_s u_k, \quad (19)$$

$$M = \frac{A_x \sinh(\kappa H)}{A_x \cosh(\kappa H) + A_y + A_z} M_s \mathbf{x}. \quad (20)$$

This expression introduces the saturation magnetisation M_s of the material and two additional constants α and κ. Due to the form of the magnetostriction tensor (Eq. (6)), only the components σ_{xx}, σ_{yy} and σ_{zz} appear in the analytical expression of the magnetisation (Eq. (20)) — see the expression of the magneto-elastic energy in equation (3).

Similarly, the magnetostriction strain is defined by the discrete summation (Eq. (21)):

$$\varepsilon^{\mu} = \sum_{k=1}^{6} f_k \varepsilon_k^{\mu}. \quad (21)$$

It is easily shown that it can be written in the form

$$\varepsilon^{\mu} = \lambda \begin{pmatrix} 1 & 0 & 0 \\ 0 & -1/2 & 0 \\ 0 & 0 & -1/2 \end{pmatrix} \quad (22)$$

with

$$\lambda = \lambda_s \left( 1 - \frac{3(A_y + A_z)}{2(A_x \cosh(\kappa H) + A_y + A_z)} \right). \quad (23)$$

Under the considered assumptions (isotropic homogeneous material), equation (23) provides an analytical expression for the magnetostriction strain of the material under any multiaxial stress state, with any orientation with respect to the magnetic field. Three material parameters are introduced: the saturation magnetostriction λ_s and two additional parameters κ and α. These parameters are related to standard material parameters by equations (16)–(18).

## 4 Particular configurations

Further simplifications can be obtained if less complex loadings are considered. The corresponding expressions are given hereafter:

- From equation (23), it is verified that under no applied stress nor magnetic field, the magnetostriction strain is zero. It is also verified that when the magnetic field is getting very high, λ tends towards λ_s.

- If no stress is applied, the magnetostriction strain reduces to equation (24):

$$\lambda(H, \mathbf{0}) = \lambda_s \left( 1 - \frac{3}{\cosh(\kappa H) + 2} \right). \quad (24)$$

- If no magnetic field is applied, the magnetostriction strain reduces to equation (25):

$$\lambda(0, \sigma) = \lambda_s \left( 1 - \frac{3(A_y + A_z)}{2(A_x + A_y + A_z)} \right). \quad (25)$$

- The magnetostriction strain under uniaxial stress and no applied field is given by equation (26). This expression describes the so-called ΔE effect [35,36]:

$$\lambda(0, \sigma_{xx}) = \lambda_s \left( 1 - \frac{3}{A_x + 2} \right). \quad (26)$$

- If we consider a uniaxial stress σ_{xx} in the direction x, equation (23) reduces to

$$\lambda(H, \sigma_{xx}) = \lambda_s \left( 1 - \frac{3}{A_x \cosh(\kappa H) + 2} \right). \quad (27)$$

## 5 Identification of parameters κ and α

The proposed model is based on three material parameters (λ_s, κ and α). λ_s is the maximum magnetostriction strain and its significance is clear. As defined in reference [32], κ can be identified from the initial anhysteretic susceptibility χ_o^o of the material under no applied stress (Eq. (28)). This relation is obtained by combining equations (17) and (18). The parameter α can be identified from equation (29), obtained by combining equations (16) and (18). α can also be extracted from a ΔE measurement curve. The ΔE effect
is the loss of linearity in the stress–strain curve of magnetic materials [35] (see Sect. 5). It is due to the magnetostriction strain that superimposes to the elastic strain. This magnetostriction strain under no applied field, and for a uniaxial stress, is given by equation (26). By considering the initial slope *p* of the Δ*E* curve, the parameter *α* can be obtained (Eq. (30)):

$$ \kappa = \frac{3\chi_o^o}{M_s}, \qquad (28) $$

$$ \alpha = \frac{9\lambda_s\chi_o^o}{2\mu_0 M_s^2}, \qquad (29) $$

$$ \alpha = \frac{3p}{\lambda_s}, \qquad (30) $$

$$ \text{with } p = \frac{\partial\lambda(0, \sigma_{xx})}{\partial\sigma_{xx}} \Big|_{\sigma_{xx}=0}. \qquad (31) $$

## 6 Model prediction

The proposed model (Eq. (23)) can provide a prediction of the magnetostriction strain of an isotropic material for any magneto-mechanical loading, i.e. for any magnetic field **H** and any stress tensor **σ**, whatever the relative orientation between them. It must be noticed that the model describes the anhysteretic behaviour, so that hysteresis effects are not included in the proposed description. As an example, Figure 1 shows the evolution of the magnetostriction strain under no applied magnetic field as a function of stress under uniaxial, equibiaxial, hydrostatic and pure shear stress states. These stress states are, respectively, defined by a diagonal stress tensor with the values (σ, 0, 0), (σ, σ, 0), (σ, σ, σ) and (σ, −σ, 0) on the diagonal. The magnetostriction strain is measured along the direction **x**.

The results show very similar evolution as those obtained for the magnetic susceptibility in reference [32]. It is worth noting that an applied hydrostatic stress has no effect on the magnetic susceptibility, which is consistent with the fact that magnetostriction strain is isochoric (trace (ε^μ) = 0). The uniaxial stress curve (u) corresponds to the description of the Δ*E* effect, as defined in reference [35].

Very few measurements are available in the literature for an experimental validation under complex loadings. Most of the available data are restricted to uniaxial configurations for which the stress is applied in the direction parallel to the magnetic field. Here, the model has been compared to magnetostriction strain measurements under uniaxial stress presented in reference [36] for an iron–cobalt alloy. They consist of anhysteretic measurements on flat samples of dimensions 2.5 × 12.5 × 110 mm³ loaded along the longest direction. The strain was measured using strain gages. Numerical simulations performed on a geometry corresponding to the experimental setup have shown that the form effect (deformation due to the magnetic forces on the specimen [37]) has little effect on the strain measurement area, so that the presented data reflect the magnetostriction strain. Figures 2 and 3 show the strain

<figure><img src="image_2.png" /><figcaption>Fig. 1. Magnetostriction strain (maximum principal strain) under no applied magnetic field as a function of stress intensity for uniaxial (u), equibiaxial (e), hydrostatic (h) and pure shear (s) stress states (modelling parameters: λ_s = 58 × 10⁻⁶, κ = 2.2 × 10⁻³ m/A, α = 10⁻⁷ Pa⁻¹).</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>Fig. 2. Longitudinal magnetostriction strain under uniaxial stress for an iron–cobalt alloy: analytical model (left) and experimental results from reference [36] (right) (modelling parameters: λ_s = 58 × 10⁻⁶, κ = 2.2 × 10⁻³ m/A, α = 10⁻⁷ Pa⁻¹).</figcaption></figure>

measurements, respectively, parallel and perpendicular to the magnetic field direction for various stress levels. The material parameters have been identified as λ_s = 58 × 10⁻⁶, κ = 2.2 × 10⁻³ m/A and α = 10⁻⁷ Pa⁻¹. The left side of the figure shows the prediction of the proposed model, and the right side shows the experimental results. Regarding the modelling results, the definition of the magnetostriction curve (Eq. (27)) incorporates naturally the Δ*E* effect, which is the magnetostriction strain as a function of stress under no applied field. In order to be consistent with the representation of experimental results, this initial strain has been removed in the representation of modelling results. The plotted numerical results are then defined as λ = λ(H, σ_{xx}) − λ(0, σ_{xx}).
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Transverse magnetostriction strain under uniaxial stress for an iron-cobalt alloy: analytical model (**left**) and experimental results from reference [36] (**right**) (modelling parameters: $\lambda_s = 58 \times 10^{-6}$, $\kappa = 2.2 \times 10^{-3}$ m/Å, $\alpha = 10^{-7}$ Pa$^{-1}$).</figcaption></figure>

Considering the simplicity of the model (only Eq. (27) is used), Figures 2 and 3 show a very satisfactory agreement between experimental and modelling results. The model tends to slightly overestimate the longitudinal strain. On the contrary, the transverse strain is underestimated by the model. Since an isotropic approach has been used for the material, the predicted transverse strain is simply half the longitudinal strain with an opposite sign. The experimental results do not show a similar ratio, suggesting an initial anisotropy is not taken into account in the model. The agreement could probably be enhanced by introducing additional material parameters, notably to describe initial anisotropy effects.

Another experimental validation can be made based on the measurement of the $\Delta E$ effect. The $\Delta E$ effect is an apparent dependency of the Young's modulus to the stress level when a tensile or compressive test is performed on a magnetic specimen [1]. It is a manifestation of the effect of stress on the magnetostriction strain (under no applied magnetic field) [35]. Indeed, the application of stress modifies the magnetic domain structure and hence generates a magnetostriction strain. This magnetostriction strain is superimposed on the elastic strain, so the stress-strain curve obtained from a tension or compression test is non-linear. A description of the $\Delta E$ effect based on this interpretation has been thoroughly discussed in reference [35]. The modelling approach presented in this paper describes the $\Delta E$ effect through equation (26). The results have been plotted in Figures 4 and 5 for an iron-cobalt alloy (different from the alloy presented in Figs. 2 and 3) and for pure iron, respectively. The experimental results are simply extracted from reference [35]. The experimental setup is the same as the one used for Figures 2 and 3. The results from the model proposed in reference [35] are also reproduced with a dotted line.

It can be seen that both models and experimental results exhibit very satisfying agreement, both in trends and values. Only two parameters are required for the proposed modelling of the $\Delta E$ effect since the parameter $\kappa$ is

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** $\Delta E$ effect for an iron-cobalt alloy: longitudinal and transverse magnetostriction strain as a function of the applied stress $\sigma$ (uniaxial), modelling (lines) and experimental (points) results from reference [35]. The dotted line reproduces the modelling results obtained in reference [35]. The plain line shows the results from equation (26) (modelling parameters: $\lambda_s = 40 \times 10^{-6}$, $\alpha = 1.25 \times 10^{-7}$ Pa$^{-1}$).</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** $\Delta E$ effect for pure iron: longitudinal and transverse magnetostriction strain as a function of the applied stress $\sigma$ (uniaxial), modelling (lines) and experimental (points) results from reference [35]. The dotted line reproduces the modelling results obtained in reference [35]. The plain line shows the results from equation (26) (modelling parameters: $\lambda_s = 8 \times 10^{-6}$, $\alpha = 2.4 \times 10^{-8}$ Pa$^{-1}$).</figcaption></figure>

used only under applied magnetic field. The modelling results from reference [35] have simply been reproduced. In particular, the material parameters, deduced from single crystal data, have not been changed. It is clear that the modelling of the $\Delta E$ effect (dotted lines) could be improved by fitting the modelling parameters to the experimental data. On the contrary, the modelling parameters used in the present approach (Eq. (26)) have been adjusted to the experimental data, which explains the slightly better agreement observed. The main difference between the two modelling approaches is that the description proposed in reference [35] is predictive (in the sense that material parameters can be deduced from the basic properties of
single crystals) but restricted to the case where no magnetic field is applied. On the contrary, the present model is more phenomenological (in the sense that the material parameters have to be fitted on macroscopic measurements) but also more general in the sense that it can combine the application of stress and magnetic field (see Figs. 2 and 3). It can again be observed that the measured transverse strain is not half the longitudinal strain in amplitude, contrarily to the hypothesis made in both models. This can be attributed to anisotropy effects not taken into account in the modelling approaches.

## 7 Conclusion

In this paper, an original analytical model for the definition of the stress-dependent magnetostriction strain of magnetic materials is proposed. It is based on a very simplified description of the energetic equilibrium underlying magnetic behaviour. The approach results in a simple formula to express the effect of multiaxial magneto-mechanical loadings on the magnetostriction strain. The multiaxiality of stress is naturally introduced in this approach, and no assumption is made on the relative orientation between stress and magnetic field. Only three materials parameters are required to describe this complex magneto-mechanical behaviour. The three parameters can be, respectively, identified from an anhysteretic curve under no applied stress, a magnetostriction curve under no applied stress and a $\Delta E$ effect measurement. They can also be fitted from experimental results under stress. This analytical model can be useful for electromagnetic design when users require the implementation of a magnetostriction model into standard structural analysis tools. Due to the strong assumptions made in the construction of the model, it is expected that limitations will be found in cases of strong anisotropy or high heterogeneities.

## References

1. B.D. Cullity, C.D. Graham, *Introduction to Magnetic Materials*, 2nd edn (John Wiley & Sons, Inc., New York, 2009)

2. N. Takahashi, H. Morimoto, Y. Yunoki, D. Miyagi, J. Magn. Magn. Mater. **320**, e925 (2008)

3. H. Ebrahimi, Y. Gao, H. Dozono, K. Muramatsu, T. Okitsu, D. Matsuhashi, IEEE Trans. Magn. **52**, 8201404 (2016)

4. K. Yamazaki, H. Takeuchi, IEEE Trans. Ind. Appl. **53**, 963 (2017)

5. F. Claeyssen, N. Lhermet, R. Le Letty, P. Bouchilloux, J. Alloys Compd. **258**, 61 (1997)

6. G. Engdahl, *Handbook of Giant Magnetostrictive Materials* (Academic Press, Cambridge, UK, 2000)

7. A.G. Olabi, A. Grunwald, Mater. Des. **29**, 469 (2008)

8. A.J. Moses, IEEE Trans. Magn. **10**, 154 (1974)

9. B. Weiser, H. Pfutzner, J. Anger, IEEE Trans. Magn. **36**, 3759 (2000)

10. Y.H. Chang, C.H. Hsu, H.L. Chu, C.P. Tseng, IEEE Trans. Magn. **47**, 2780 (2011)

11. R.M. Bozorth, *Ferromagnetism* (Van Norstand, New York, 1951)

12. P.I. Anderson, A.J. Moses, H.J. Stanbury, IEEE Trans. Magn. **43**, 3467 (2007)

13. Y. Kai, Y. Tsuchida, T. Todaka, M. Enokizono, IEEE Trans. Magn. **50**, 6786404 (2014)

14. A. Belahcen, D. Singh, P. Rasilo, F. Martin, S.G. Ghalamestani, L. Vandevelde, IEEE Trans. Magn. **51**, 2001204 (2015)

15. M.J. Sablik, D.C. Jiles, IEEE Trans. Magn. **29**, 2113 (1993)

16. O. Bottauscio, A. Lovisolo, P.E. Roccato, M. Zucca, C. Sasso, R. Bonin, IEEE Trans. Magn. **44**, 3009 (2008)

17. H. El Bidweihy, E. Della Torre, Y. Jin, L.H. Bennett, M. Ghahremani, IEEE Trans. Magn. **48**, 3360 (2012)

18. L. Hirsinger, G. Barbier, C. Gourdin, R. Billardon, Mech. Electromagn. Mater. Struct. **19**, 54 (2000)

19. K. Fonteyn, A. Belahcen, R. Kouhia, P. Rasilo, A. Arkkio, IEEE Trans. Magn. **46**, 2923 (2010)

20. P. Rasilo, D. Singh, U. Aydin, F. Martin, R. Kouhia, A. Belahcen, A. Arkkio, IEEE Trans. Magn. **52**, 7300204 (2016)

21. N. Buiron, L. Hirsinger, R. Billardon, J. Phys. IV **9**, 187 (1999)

22. W.D. Armstrong, J. Appl. Phys. **91**, 2202 (2002)

23. L. Daniel, O. Hubert, R. Billardon, Comput. Appl. Math. **23**, 285 (2004)

24. L. Daniel, O. Hubert, N. Buiron, R. Billardon, J. Mech. Phys. Solids **56**, 1018 (2008)

25. L. Daniel, N. Galopin, Eur. Phys. J. Appl. Phys. **42**, 153 (2008)

26. L. Daniel, M. Rekik, O. Hubert, Arch. Appl. Mech. **84**, 1307 (2014)

27. D. Vanoost, S. Steentjes, J. Peuteman, G. Gielen, H. De Gersem, D. Pissoort, K. Hameyer, J. Magn. Magn. Mater. **414**, 168 (2016)

28. L. Daniel, O. Hubert, M. Rekik, IEEE Trans. Magn. **51**, 7300704 (2015)

29. L. Bernard, X. Mininger, L. Daniel, G. Krebs, F. Bouillault, M. Gabsi, IEEE Trans. Magn. **47**, 2171 (2011)

30. L. Bernard, L. Daniel, IEEE Trans. Magn. **51**, 7002513 (2015)

31. M. Liu, O. Hubert, X. Mininger, F. Bouillault, L. Bernard, IEEE Trans. Magn. **52**, 8000212 (2016)

32. L. Daniel, IEEE Trans. Magn. **49**, 2037 (2013)

33. O. Hubert, L. Daniel, J. Magn. Magn. Mater. **323**, 1766 (2011)

34. O. Hubert, S. Lazreg, IEEE Trans. Magn. **48**, 1277 (2012)

35. L. Daniel, O. Hubert, Eur. Phys. J. Appl. Phys. **45**, 31101 (2009)

36. O. Hubert, L. Daniel, IEEE Trans. Magn. **46**, 401 (2010)

37. L. Hirsinger, R. Billardon, IEEE Trans. Magn. **31**, 1877 (1995)

**Open Access** This article is distributed under the terms of the Creative Commons Attribution License https://creativecommons.org/licenses/by/4.0 which permits unrestricted use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.

**Cite this article as:** Laurent Daniel, An analytical model for the magnetostriction strain of ferromagnetic materials subjected to multiaxial stress, Eur. Phys. J. Appl. Phys. **83**, 30904 (2018)