Published: 27 September 2024

A Zhu-Wang-Tang damage constitutive model for sintered NdFeB considering crack spacing

Yuwei Yao1
Guolai Yang2
Lei Li3
Liqun Wang4
1, 2, 3, 4School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing, People’s Republic of China
Corresponding Author:
Guolai Yang
Views 3
Reads 0
Downloads 21

Abstract

Currently, research on the dynamic damage characteristics of sintered NdFeB is still in its early stages. Previous dynamic mechanical experiments on sintered NdFeB have shown that it is a strain rate-sensitive material. So, it is necessary to establish an appropriate damage constitutive model to describe the dynamic mechanical behavior of sintered NdFeB, with the aim of expanding its practical application range. This paper first establishes a damage evolution model by combining the Weibull distribution with a wing crack propagation model that considers crack spacing. Then, the damage model is integrated with the Zhu-Wang-Tang (ZWT) constitutive model to create the ZWT damage constitutive model. The model is fitted against previous experimental data to determine the specific parameters. Finally, to verify the accuracy of the newly established damage constitutive model, it is compared with the damage constitutive model proposed by Li. The comparison results fully affirmed the accuracy of the newly established ZWT damage constitutive model.

A Zhu-Wang-Tang damage constitutive model for sintered NdFeB considering crack spacing

Highlights

  • Consider crack spacing in damage models.
  • Combining the damage model with Zhu-Wang-Tang constitutive model.
  • Developing a new constitutive model of sintered NdFeB.

1. Introduction

Sintered NdFeB rare earth permanent magnet materials are functional materials that were discovered by Sagawa [1] in the 1980s. This material possesses excellent magnetic properties and has spurred the development of industries and technologies like computer hard disk technology, nuclear magnetic resonance, and scanners [2]. Sintered NdFeB serves as energy storage devices that enable efficient conversion of energy and information. Consequently, the volume of NdFeB required to generate the same magnetic field is smaller compared to other magnets, which aligns well with the contemporary demands for lightweight and integrated production. During the manufacturing process and practical applications, these materials inevitably encounter impact loads, making it essential to investigate the dynamic response of sintered NdFeB.

Sintered NdFeB has fundamental mechanical properties that include bending strength, compressive strength, elastic modulus, fracture toughness, and hardness [3, 4]. Many scholars have conducted research on its basic mechanical properties. Szymura et al. [5] examined S-NdFeB permanent magnets with different compositions and found that impact strength, bending strength, and fracture toughness increased with the enhancement of coercivity. Wang Huanran et al. [6] used a split Hopkinson pressure bar to investigate the fracture of NdFeB magnets under uniaxial compression. While these studies primarily focus on the basic mechanical properties of NdFeB, no damage constitutive equations have yet been established to describe the dynamic behavior of sintered NdFeB. Therefore, its dynamic damage constitutive equation requires further exploration.

Li et al. [7] conducted experiments on the dynamic characteristics of sintered NdFeB and proposed a damage constitutive model. However, it did not consider the influence of crack spacing in the damage process. Therefore, this paper first fully considers the influence of crack spacing in the damage process, and constructs the damage model of sintered NdFeB based on the weibull distribution and the swing crack growth model. Then, the damage evolution process was combined with the ZWT constitutive model to establish the ZWT damage constitutive model (ZWTDCM). Based on previous dynamic mechanical characteristics tests of sintered NdFeB, ZWTDCM of sintered NdFeB was fitted. Finally, the newly established model was compared with the model proposed by Li. The comparison results fully demonstrated the rationality of the newly established model.

2. Damage constitutive model

2.1. The definition of damage

Referring to the damage evolution process established by Ren [8] based on the Weibull distribution, damage is defined as:

1
D=γNl,N=kεe,

where, γ is the crack shape parameter, and N is the crack density. k and e are material characteristic parameters, and l is the crack propagation length.

Take the derivative of Eq. (1) with respect to time:

2
D˙=γekεe-1ε˙l3+3γkεel2l˙,

where, l˙ is the crack propagation speed.

2.2. Wing crack propagation model considering crack spacing

The condition for crack propagation under impact is [9]:

3
KID=KICd,

where, KICd is the fracture toughness. The expression for dynamic stress concentration factor KID is [10]:

4
kl˙=1-l˙l˙m1-l˙2l˙m-1,l˙m0.4lR,KID=kl˙KI,

where, l˙m is the maximum crack propagation speed, and lR is the Rayleigh wave speed. From Eq. (4), l˙ can be obtained:

5
l˙=l˙m1-KICd2KI-KICd.

Fig. 1 shows the equivalent representation of the microstructure of sintered NdFeB, where σ1 and σ2 represent the far-field loads applied to a single crack, and the initial length of the crack is 2a. When the stress concentration factor at the tip of the micro-crack reaches the fracture toughness, frictional slip occurs at points Q and Q', resulting in the formation of new tensile cracks P and P'. The initial nucleation direction of the tensile cracks makes an angle (φ) with σ1. As the cracks propagate, the nucleation direction eventually becomes parallel to σ1.

Fig. 1Wing crack propagation model

Wing crack propagation model

Under the action of far-field loads, KI at the tip of an I-type crack, considering the spacing of micro-cracks, is given by [11]:

6
KI=Fsinθdsinπl+l*d-σ22dtanπl2d,F=2lτ*=2l12σ1-σ2sin2θ-12fσ1+σ2-σ1-σ2cos2θ,l*=0.27l,

where, F is the normal force on the crack surface, τ* is the driving force for crack sliding, and f is the friction coefficient.

2.3. ZWT damage constitutive model and Li’s model

The ZWT constitutive model [12] is employed to describe the stress-strain behavior. The model consists of only one spring and one Maxwell body. The expression is as follows:

7
σ=E0ε+E10tε˙exp-t-τθ1dτ,

where, E0 is a nonlinear coefficient. E10tε˙exp-t-τθ1dτ represents a Maxwell body with a relaxation time used to describe the mechanical behavior under high strain rates. At the same time, to fully consider the damage during the loading process, the damage evolution process is combined with the ZWT constitutive model, resulting as follows:

8
σ=1-DE0ε+E10tε˙exp-t-τθ1dτ,D=γkεel.

The damage constitutive model proposed by Li [7] is shown as follows:

9
σ=(1-D)×E0ε+E10tε˙exp-t-τθ1dτ,
D=Nln1+ll*×a2cos2φπ(sin22φ-4f2sin2φ-2fsin2φ+fsin4φ)+a222fsin2φ-24fsin4φ-22f2sin4φ1+Nln1+ll*×a2cos2φπsin22φ-4f2sin2φ-2fsin2φ+fsin4φ+a222fsin2φ-24fsin4φ-22f2sin4φ.

2.4. Comparative analysis of the results from the two damage constitutive models

Referring to the experimental results in [7], the parameters of the newly established model were fitted in ISTOPT. Table 1 shows the fitted parameters. From the fitting results, it can be observed that the fitting accuracies of the new model at strain rates of 700, 1200, and 1500 are 0.99, 0.94, and 0.96, respectively, while Li's model has fitting accuracies of 0.97, 0.96, and 0.95 at the same strain rates. The two models show little difference in fitting accuracy, both exhibiting good fitting performance. The comparison of the two models is shown in Fig. 2.

Fig. 2Comparison of the two models

Comparison of the two models

a) 700 s-1

Comparison of the two models

b) 1200 s-1

Comparison of the two models

c) 1500 s-1

In the initial stage, the experimental stress-strain curves are nearly linear, and during this phase, both models can fit the trend of the changes quite well. Subsequently, the stress-strain relationship exhibits nonlinear behavior, at which point both models show some deviation from the experimental data, but the errors remain within a reasonable range. Therefore, it can be seen that the trends in strain change for both models are also consistent with the expected patterns. Thus, the accuracy and rationality of ZWTDCM established in this paper, which considers the effect of crack spacing, have been thoroughly validated.

Table 1Parameters for the new model

Parameters
Values
Parameters
Values
k
4×109
2a (μm)
5
e
2.5
θ
π/4
KIC (Pa/m)
4×106
d (μm)
120
f
0.7
E0(MPa)
538330
E1 (MPa)
2000
θ1 (μs)
0.55

3. Conclusions

In this paper, the influence of crack spacing in the damage evolution process is fully considered, and the damage evolution process of sintered NdFeB is given by combining the Weibull distribution with the swing crack growth model considering crack space. In addition, the ZWT damage constitutive model was constructed by combining the damage process with the ZWT constitutive model. Based on previous experimental data on the dynamic mechanical characteristics of sintered NdFeB, the parameters of ZWTDCM are determined and ZWTDCM was compared with Li’s model. On the one hand, the two models have high fitting accuracy; on the other hand, both can reflect the stress-strain laws of sintered NdFeB under impact. However, compared with Li’s model, this model fully considers the influence of crack gaps during the damage process, so the accuracy and superiority of this model have been further verified.

References

  • M. Sagawa, S. Fujimura, N. Togawa, H. Yamamoto, and Y. Matsuura, “New material for permanent magnets on a base of Nd and Fe (invited),” Journal of Applied Physics, Vol. 55, No. 6, pp. 2083–2087, Mar. 1984, https://doi.org/10.1063/1.333572
  • B. E. Davies, R. S. Mottram, and I. R. Harris, “Recent developments in the sintering of NdFeB,” Materials Chemistry and Physics, Vol. 67, No. 1-3, pp. 272–281, Jan. 2001, https://doi.org/10.1016/s0254-0584(00)00450-8
  • M. Lin, H. J. Wang, P. P. Yi, and A. R. Yan, “Effects of excessive grain growth on the magnetic and mechanical properties of hot-deformed NdFeB magnets,” Journal of Magnetism and Magnetic Materials, Vol. 322, No. 15, pp. 2268–2271, Aug. 2010, https://doi.org/10.1016/j.jmmm.2010.02.023
  • H. J. Wang, A. H. Li, M. G. Zhu, and W. Li, “Sintered Nd-Fe-B magnets with improved impact stability,” Journal of Magnetism and Magnetic Materials, Vol. 307, No. 2, pp. 268–272, Dec. 2006, https://doi.org/10.1016/j.jmmm.2006.04.012
  • S. Szymura et al., “Domain structure, magnetic and mechanical properties of Nd-Fe-B magnets with different grain size,” Physica Status Solidi (a), Vol. 141, No. 2, pp. 435–443, 1994.
  • H. Wang, Y. Wan, D. Chen, G. Lei, and C. Ren, “Dynamic fracture of sintered Nd-Fe-B magnet under uniaxial compression,” Journal of Magnetism and Magnetic Materials, Vol. 456, pp. 358–367, Jun. 2018, https://doi.org/10.1016/j.jmmm.2018.02.058
  • L. Li, G. Yang, and Q. Sun, “Dynamic damage study of sintered NdFeB in electromagnetic buffer under intensive impact load: Experiments, numerical simulation, and interval uncertain optimization,” Mechanics of Advanced Materials and Structures, Vol. 29, No. 24, pp. 3523–3539, Sep. 2022, https://doi.org/10.1080/15376494.2021.1905116
  • H.-L. Ren and J.-G. Ning, “Mechanical characteristics and constitutive model of alumina ceramic subjected to shock loading,” Transactions of Beijing Institute of Technology, Vol. 27, No. 9, pp. 761–764, 2007.
  • G. Ravichandran and G. Subhash, “A micromechanical model for high strain rate behavior of ceramics,” International Journal of Solids and Structures, Vol. 32, No. 17-18, pp. 2627–2646, Sep. 1995, https://doi.org/10.1016/0020-7683(94)00286-6
  • H. Deng and S. Nemat-Nasser, “Dynamic damage evolution in brittle solids,” Mechanics of Materials, Vol. 14, No. 2, pp. 83–103, 1992.
  • R. Huilan and N. Jianguo, “A dynamic compressive damage model for ceramics,” Acta Mechanica Solida Sinica-Chinese Edition, Vol. 27, No. 3, p. 303, 2006.
  • L.-L. Wang, “Stress wave propagation for nonlinear viscoelastic polymeric materials at high strain rates,” Journal of Mechanics, Vol. 19, No. 1, pp. 177–183, May 2011, https://doi.org/10.1017/s1727719100004184

About this article

Received
28 August 2024
Accepted
25 September 2024
Published
27 September 2024
SUBJECTS
Mathematical models in engineering
Keywords
Weibull distribution
wing crack propagation
damage model
ZWT constitutive model
Acknowledgements

This research was financially supported by the “National Natural Science Foundation of China” (Grant No. 52105106), the “China National Postdoctoral Program for Innovative Talents” (Grant No. BX20230493), and the “Jiangsu Province Natural Science Foundation” (Grant No. BK20210342).

Data Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflict of interest

The authors declare that they have no conflict of interest.