稀有金属(英文版) 2020,39(10),1145-1150
Numerical simulation of single roller melt spinning for NdFeB alloy based on finite element method
Xu-Chao Wang Ming Yue Dong-Tao Zhang Wei-Qiang Liu Ming-Gang Zhu
College of Materials Science and Engineering,Beijing University of Technology
Central Iron and Steel Research Institute
作者简介:*Xu-Chao Wang,e-mail:362577891@qq.com;
收稿日期:26 November 2017
基金:financially supported by the National Natural Science Foundation of China(No.51571064);the National Basic Research Program of China (No.2014CB643701);
Numerical simulation of single roller melt spinning for NdFeB alloy based on finite element method
Xu-Chao Wang Ming Yue Dong-Tao Zhang Wei-Qiang Liu Ming-Gang Zhu
College of Materials Science and Engineering,Beijing University of Technology
Central Iron and Steel Research Institute
Abstract:
The numerical simulation model of single roller melt spinning for rapid quenching process of NdFeB alloy was built,and the vacuum chamber,cooling roller and sample were taken into account as a system.The existing mature technology was in order to verify the correctness of simulation.The rapid quenching ribbons with different roll speeds were used as the simulation objects.The results of the numerical simulation and experiments show that the validity of the model has been testified and the reasons of the formation of complete quenching ribbons and by-product have been explained.The experimental thickness of the ribbons is proportional to the theoretical thickness.In the same spray condition,with the roll speed increasing,the thickness decreases linearly.At the speed range of25-30 m·s-1,the simulated calculation date is close to the experimental date,which can be considered as an ideal technological parameter.
Keyword:
NdFeB; Melt spinning; Finite element method; Numerical simulation;
Received: 26 November 2017
1 Introduction
Since the outstanding magnetic properties of NdFeB were found in 1983,it has been studied in many ways
[
1,
2,
3,
4,
5]
.Fully dense anisotropic NdFeB magnets can be produced by the rapid quenching powders
[
6,
7,
8,
9]
,followed by the hot-pressing process
[
10,
11]
.Now,vacuum rapidly quenched technology is a significant method of preparing amorphous and nanocrystalline alloys,which is widely used to prepare NdFeB powders
[
12,
13]
.The magnetic fine powders with particle size in an appropriate size have many applications
[
14]
.The powders can be used to prepare the bonded magnets
[
15]
and the hot-deformed magnets
[
16,
17,
18,
19,
20]
.The NdFeB magnets,which were prepared by hot-pressing heat deformation with powders,can have excellent magnetic property of about 432.5 kJ·m-3
[
21]
.
The process parameters have a direct effect on the magnetic powder performance,and the preparation process of the rapid quenching technique is relatively complex
[
22,
23,
24,
25,
26,
27]
,so we need to accumulate a large amount of experimental data.Therefore,the simulation of the rapid quenching process is necessary and can be helpful to the experiment.Owing to the complexity of the process,it is necessary to take different degrees of simplification
[
28,
29,
30]
.But the simulation results are still complicated and the research on the rapid quenching mechanism is not thorough enough.Yang et al.
[
31,
32,
33]
have studied the cooling speed of the process with numerical simulation,in which the change of cooling speed along the thickness of the melt-spun ribbons is not considered.In this paper,based on the simplification of the NdFeB rapid quenching process,a preliminary study on the fluid motion of the single-roll rapid quenching process was carried out,and the influence of the rapid quenching process parameters on the preparation of the melt-spun ribbons was discussed.
2 Numerical simulation of rapid quenching process
2.1 Experimental design and numerical model establishment
Figure 1 shows the schematic diagram of the melt spinning.According to the rapid quenching technology and the experimental data,the appropriate process parameter has been defined.The cast ingots were melt-spun to form ribbons,about 40 g,which was conducted by using a vacuum melt spinner with a quartz tube with a nozzle diameter of1 mm,a nozzle height of 3 mm,a roller speed of15-45 m·s-1 and a vacuum of less than 1×10-2 Pa at the temperature of 1310℃
[
34]
.The diameter of quartz tube is 21.4 mm,and the ejection pressure is 0.05 MPa.The ejection pressure determines the injection speed of the molten fluid.Some property parameters of the NdFeB alloy were investigated upon experimentation for the numerical simulation of rapid quenching.The density of the NdFeB melt is 6.94 kg2m-3,and the viscosity is 1.4 Pa-s at the temperature of 1310℃
[
35]
.
When NdFeB was smelted,the melt is forming a curved surface at the nozzle,due to the surface tension.As shown in Fig.2,the nozzle is required to get accurate injection velocity field,and the arc is pided into 16 parts,and the closer to the bottom is,the more compact it is,and in order to ensure the accuracy,FLULID141 and FLULID142 unit were used.Finally,the entire mesh model is pided into1882 nodes.At the top edge of the model,the pressure of0.05 MPa was applied.The velocity of injection at the remaining edges is 0 m·s-1.
![](/web/fileInfo/upload/magazine/14843/370770/XYJS202010006_00800.jpg)
Fig.1 Schematic diagram of rapid quenching
![](/web/fileInfo/upload/magazine/14843/370770/XYJS202010006_00900.jpg)
Fig.2 Injection process finite element model
3 Results and discussions
To run the program,the equivalent of NdFeB alloy was the molten state for the liquid jet.Figure 3 clearly shows that speed is focused on the nozzle mouth and the maximum value is 2.2585 m·s-1.The edge velocity of the nozzle decreases because the edge parts obviously gain resistance by quartz tube.
3.1 Simulation analysis of injection time
In this simulation,the law of conservation of mass is adopted.According to the quality of NdFeB,Eq.(1) can be obtained.
![](/web/fileInfo/upload/magazine/14843/370770/XYJS202010006_13100.jpg)
where Rout is the nozzle radius;vout is the velocity of injection;ρis the density of NdFeB;t is the ending time of the injection;M is the quality of NdFeB.
At last,the simulation calculated the completion time of the injection of 40 g melt;it is 3.24 s,consistent with the time in the actual experiment.
![](/web/fileInfo/upload/magazine/14843/370770/XYJS202010006_01700.jpg)
Fig.3 Distribution of velocity at nozzle
3.2 Simulating width of NdFeB ribbons
In the ejection process of the rapid quenching experiments,not all of the molten liquid will shoot vertically.It will also spray on both sides.And the width of the NdFeB ribbons is related to the distance of the nozzle to the rollers.Figure 4is the model for the broadening of the ribbons.
As shown in Fig.5,the speed of 16 nodes is got in the xor y-direction by simulation.And the maximum width of each point is calculated by the following three formulas.
![](/web/fileInfo/upload/magazine/14843/370770/XYJS202010006_02100.jpg)
where RM is the width radius of the ribbons;Hy is the distance of the roller;Vy is the initial velocity of the nozzle at each point of the nozzle;Vx is the initial velocity of the injection horizontal direction of the nozzle;Lx is thedistance of the nozzle at each point of the nozzle;H is theheight of the roller wheel at the nozzle center;n is the nthcurve of the jet.
![](/web/fileInfo/upload/magazine/14843/370770/XYJS202010006_02300.jpg)
Fig.4 A model for broadening of ribbons
![](/web/fileInfo/upload/magazine/14843/370770/XYJS202010006_02700.jpg)
Fig.6 A maximum width of ribbons at nozzle
![](/web/fileInfo/upload/magazine/14843/370770/XYJS202010006_02800.jpg)
Fig.7 Comparison of theoretical data and experimental data of ribbons thickness
![](/web/fileInfo/upload/magazine/14843/370770/XYJS202010006_02900.jpg)
Fig.5 Velocity curves at nozzle.a y vector and b x vector
As shown in Fig.6,the maximum width of the rapid quenching band is obtained at the 13rd node,2.16 mm.In Fig.3,liquid and spray center vertical velocity is high and consistent,but the edge speed is small,leading to that the center part of melt first reaches the cooling roller and the edge parts of the melt have not reached the roller with the high-temperature condition.The first melt reached the cooling roller,and then condensed quickly.And then the edges melt reached and cooled later.It is called the discontinuous rapid quenching process.The results explain why in the rapid quenching,we get not only the complete rapid quenching ribbons but also the formation of some byproducts,such as round balls or irregular ribbons.
![](/web/fileInfo/upload/magazine/14843/370770/XYJS202010006_03100.jpg)
Fig.8 Temperature field cloud picture of ribbons
![](/web/fileInfo/upload/magazine/14843/370770/XYJS202010006_03200.jpg)
Fig.9 Curve versus temperature of model at marginal area and central zone
3.3 Simulating quenching speed and ribbon thiclkness
According to the law of conservation of mass,Eq.(5) can be obtained.
![](/web/fileInfo/upload/magazine/14843/370770/XYJS202010006_03500.jpg)
where h is the ribbon thickness;vout is the average velocity of vertical direction of the nozzle;y is the actual width of the belt;vrotate is the roll speed.
From Fig.7,it can be seen that the experimental thickness of the ribbons,which is prepared by the different quenching speed,is proportional to the theoretical thickness.In the same spray condition,with the roll speed increasing,the thickness decreases linearly.The experimental thickness variation is consistent with Eq.(5).At the speed range of 25-30 m·s-1,the simulated calculation date is close to the experimental date,which can be considered as an ideal technological parameter.
3.4 Temperature field cloud picture of ribbons
Figure 8 exhibits the temperature field cloud picture of the ribbons at 0.04 s-1 departure time.With the elongation of the departure time,the temperature in the ribbons is stable,the high temperature distribution of the edge part is reduced,and the whole is presented as a rectangular ring cooling.When the departure time was 0.04 s,the highest temperature of the ribbons is about 272℃
![](/web/fileInfo/upload/magazine/14843/370770/XYJS202010006_04000.jpg)
Fig.10 a Thickness and b grain size of ribbons at roll speed of 30 m·s
![](/web/fileInfo/upload/magazine/14843/370770/XYJS202010006_04100.jpg)
Fig.11 Hysteresis loops for melt-spun ribbons at roll speed of30 m.s-1
In Fig.9,with the extension of departure time,the temperature of the edge and the center part of the ribbons are linear.The cooling rate of the rapid quenching band is also changing.The calculated maximum cooling rate at the edge of the area is about 105℃·s-1.
3.5 Test of micros truc ture and magnetic property of ribbons
From Fig.10,it can be seen that the thickness of the ribbon is about 28.7μm and the average grain size of the ribbons is about 100 nm.The ribbon was prepared by the melt spinning at 30 m·s-1.In Fig.11,the coercive force of the ribbons is about 1274.3 kA.m-1,which was test by vibrating sample magnetometer (VSM).
4 Conclusion
In this study,the rapid quenching process was simulated by the Ansys.The simulation shows that the completion time of the injection of 40 g melt is 3.24 s and the maximum speed value is 2.2585 m·s-1.The results explain why in the rapid quenching,we get not only the complete rapid quenching ribbons but also the formation of some byproducts.The experimental thickness of the ribbons is proportional to the theoretical thickness.In the same spray condition,with the roll speed increasing,the thickness decreases linearly.At the speed range of 25-30 m·s-1,the simulated calculation date is close to the experimental date,which can be considered as an ideal technological parameter.The thickness of the ribbon is about 28.7μm,and the average grain size of the ribbons is about 100 nm.
参考文献
[1] Liu WQ,Chang C,Yue M,Yang JS,Zhang DT,Zhang JX,Liu YQ.Coercivity,microstructure,and thermal stability of sintered Nd-Fe-B magnets by grain boundary diffusion with TbH3nanoparticles.Rare Met 2017;36(9):718.
[2] Zheng ZG,Zhong XC,Su KP,Yu HY,Liu ZW,Zeng DC.Magnetic properties and large magnetocaloric effects in amorphous Gd-Al-Fe alloys for magnetic refrigeration.Sci China Phys Mech Astron.2011;54(7):1267.
[3] Zhao W,Liu Y,Li J,Wang RQ,Qiu YC.Microstructure and magnetic properties of hot-deformed anisotropic NdFeB magnets prepared from amorphous precursors with different crystallization proportions.Rare Met.2017;36(4):268.
[4] Liu WQ,Chang C,Yue M,Yang JS,Zhang DT,Zhang JX,Liu YQ.Coercivity,microstructure,and thermal stability of sintered Nd-Fe-B magnets by grain boundary diffusion with TbH3nanoparticles.Rare Met 2017;36(9):718.
[5] Li JJ,Li AH,Zhu MG,Pan W,Li W.Study on corrosion behaviors of sintered Nd-Fe-B magnets in different environ mental conditions.J Appl Phys.2011;109(7):07A744.
[6] Zhang X,Ma YT,Zhang B,Li Y,Lei MK,Wang FH,Zhu MG,Wang XC.Corrosion behavior of hot-pressed nanocrystalline NdFeB magnet in a simulated marine atmosphere.Corros Sci.2014;87:156.
[7] Betancourt I,Davies HA.Exchange coupled nanocomposite hard magnetic alloys.Mater Sci Technol.2010;26:5.
[8] El-Moneim AA,Gutfleisch O,Plotnikov A,Gebert A.Corrosion behaviour of hot-pressed and die-upset nanocrystalline NdFeB-based magnets.J Magn Magn Mater.2002;248(1):121.
[9] Lee RW,Brewer EG,Schaffel NA.Processing of Neodymium-iron-boron melt-spun ribbons to fully dense magnets.IEEE Trans Magn.1985;21(5):1958.
[10] Yi PP,Lin M,Chen RJ,Lee D,Yan AR.Enhanced magnetic properties and bending strength of hot deformed Nd-Fe-B magnets with Cu additions.J Alloys Compd.2010;491(1-2):605.
[11] Chen WH,Li W,Li CJ.The role of Nb addition in Nd-Fe-B sintered magnets with high performance.J Alloys Compd.2001;319(1-2):280.
[12] Dou YW.Progress of magnetic materials.J Funct Mater.2014;45(10):10001.
[13] Chang WC,Wu SH,Ma BM,Bounds CO.The effects of La-substitution on the micro structure and magnetic properties of nanocomposite NdFeB melt spun ribbons.J Magn Magn Mater.1997;167(1-2):65.
[14] Yoshida Y,Kasai Y,Watanabe T.Hot workability of melt-spun NdFe-Co-B magnets.J Appl Phys.1991;69(8):5841.
[15] Liu Y,Tu MS,Mao XZ.Effect of plastic on the properties of bonded NdFeB permanent magnet.J Funct Mater.1995;26(2):170.
[16] Brown DN,Chen Z,Guschl P,Campbell P.Developments with melt spun RE-Fe-B powder for bonded magnets.J Magn Magn Mater.2006;303:e371.
[17] Li W,Zhu MG.Research progress of high performance metallic permanent magnetic materials.Mater China.2001;28(9-10):62.
[18] Hu ZH,Qu HJ,Zhao JQ,Luo C,Li J,Liu Y.Enhanced mechanical properties in die-upset Nd-Fe-B magnets via die-upsetting process.J Rare Earths.2012;30(11):1112.
[19] Li AH,Li W,Lai B,Wang HJ,Zhu MG.Investigation on microstructure,texture,and magnetic properties of hot deformed Nd-Fe-B ring magnets.Journal of Applied Physics.2010;107(9):09A7251.
[20] Lai B,Li YF,Wang HJ,Li AH,Zhu MG,Li W.Quasiperiodic layer structure of die-upset NdFeB magnets.J Rare Earths.2013;31(7):679.
[21] Saito T,Fujita M,Kuji T,Fukuoka K,Syono Y.The development of high performance Nd-Fe-Co-Ga-B die upset magnets.J Appl Phys.1998;83(11):6390.
[22] Li W,Zhu MG.High property rare-earth permanent magnetic materials and its pivotal preparation technique.Chin J Nonferrous Met.2004;14(Z1):332.
[23] Zhu MG,Li W,Li XM.Exchange coupling and coercivity of nanocrystalline Nd_(4.5)(Fe,Ga,Co)_(77.5)B_(18)bonded magnets.Acta Metall Sin.2003;39(2):217.
[24] Yue M,Tian M,Zhang JX,Zhang DT,Niu PL,Yang F.Microstructure and magnetic properties of anisotropic Nd-Fe-B magnets produced by spark plasma sintering technique.Mater Sci Eng B.2006;131(1-3):18.
[25] Song J,Yue M,Zuo JH,Zhang ZR,Liu WQ,Zhang DT,Zhang JX,Guo ZH,Li W.Structure and magnetic properties of bulk nanocrystalline Nd-Fe-B permanent magnets prepared by hot pressing and hot deformation.J Rare Earths.2013;31(7):674.
[26] Yue M,Niu PL,Li YL,Zhang DT,Liu WQ,Zhang JX,Chen CH,Liu S,Lee D,Higgins A.Structure and magnetic properties of bulk isotropic and anisotropic Nd2Fe_(14)B/α-Fe nanocomposite permanent magnets with differentα-Fe contents.J Appl Phys.2008;103(7):07E1011.
[27] Lai B,Li YF,Wang HJ,Wang XC,Li AH,Zhu MG,Li W,Zhang Y.Effect of wheel speed and pressure difference on melt-spun NdFeB powders and die-upset magnets.J Funct Mater.2014;45(3):03083.
[28] Yoshida Y,Kasai Y,Watanabe T,Shibata S,Panchanathan V,Croat JJ.Hot workability of melt-spun NdFe-Co-B magnets.J Appl Phys.1991;69(8):5841.
[29] Brown DN,Smith B,Ma BM,Campbell P.The dependence of magnetic properties and hot workability of rare earth-ironboride magnets upon composition.IEEE Trans Magn.2004;40(4):2895.
[30] Lai B,Li YF,Wang HJ,Li AH,Zhu MG,Wei L,Zhang Y.Model of temperature field for the preparation process of melt-spun NdFeB powders.J Rare Earths.2014;32(6):514.
[31] Yang YS,Dong WH,Chen XM.Transient heat transfer model for rapid solidification of metals.Acta Metall Sin.2003;39(3):249.
[32] Andreev SV,Bartashevich MI,Pushkarsky VI,Maltsev VN,Pamyatnykh LA,Tarasov EN,Kudrevatykh NV,Goto T.Law of approach to saturation in highly anisotropic ferromagnets Application to NdFeB melt-spun ribbons.J Alloys Compd.1997;260(1-2):196.
[33] Pollard RJ,Parker SFH,Grundy PJ.The effect of quench rate on the microstructure and coercivity of some Nd-Fe-B ribbons.J Magn Magn Mater.1988;75(3):239.
[34] Wang XC,Zhu MG,Li W,Zheng LY,Zhao DL,Du X,Du A.The microstructure and magnetic properties of melt-spun CeFeB ribbons with varying Ce content.Electron Mater Lett.2015;11(1):109.
[35] Gu M,Li L,Sun H,Dong PY.Electric processing research on temperature field and thermal stress of NdFeB material.Mach Des Manuf.2014;02:106.