Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 13, Article number: 19930 (2023 ) Cite this article Huge Neodymium Magnets
To research the magnetic field and mechanical characteristics of the permanent magnet governor, the static magnetic field of the sector permanent magnet is analyzed by the molecular current method in the permanent magnet governor. The magnetic flux distribution is acquired at any spatial position. Comparing the analytical value with the simulation value, the results show that they are basically consistent. Based on the analytical formula, the influence of the radial position, radial length, thickness, and pole number on the magnetic induction intensity of the permanent magnet governor is studied. Thus, it provides the theoretical reference for the structural optimized design. At the same time, a test bench was set up to measure the magnetic induction intensity. The calculation and experimental results show that the magnetic induction strength of the permanent magnet is increased by 27.5%, the axial component of the air gap flux density is increased by 14.3%, and the permanent magnet material is reduced by 7.84%.
Researchers have developed a variety of permanent magnet devices at home and abroad, the permanent magnet governor is one of them1. Permanent magnet governors can be categorized as disk governors and cylinder governors based on their structural types. Disk governors are the most popular in engineering because of their simple construction and easy installation. With high reliability, high efficiency and energy saving, vibration isolation, light load starting, and adaptable ability in harsh environments, permanent magnet governor have attracted more and more attention from energy-consuming enterprises as a national key energy-saving technology2. At present, the research of permanent magnet governors mainly focuses on performance research, structural design optimization, and energy-saving applications. By researching the transient and steady-state of axial field eddy currents, LUBIN T used the vector magnetic potential approach to obtain the analytical equation of torque, which is highly consistent with the finite element simulation results3; Liu et al.4 proposed a hybrid high-efficiency permanent magnet governor, and established the mathematical model of temperature field and conducted simulation analysis. The results showed that the casing temperature with cooling ribs was significantly lower than that without cooling ribs; Wang et al., used a bidirectional coupling method to simulate and analyze the temperature field. The results showed that the bidirectional coupling simulation was more accurate than the unidirectional coupling simulation5; Seo et al., established an analytical model for the synchronous permanent magnet coupler based on the vector magnetic potential method. The analytical model accuracy was confirmed by comparison with the FEM analysis and experimental results6; Wu et al.7 established a numerical model to study the structural parameters effect of the disk permanent magnet governor on the speed regulation performance, and the accuracy of the derived equation was proved by comparison with the FEM simulation results; Jin et al.8 obtained the influence relationship of an air gap, conductor shape, poles number and speed difference on the output torque of permanent magnet governor by using Ansoft Maxwell software; Li used the equivalent magnetic circuit method to study the magnetic field distribution, and the structural parameters influence on transmission performance of the axial radial combined permanent magnet governor. The results show that the effectiveness and model accuracy are in good agreement with the finite element method and experimental results9. Zhang and Li10 and Yao11 have studied the energy-saving principle, energy-saving benefits, and technical characteristics of permanent magnet regulators by the equivalent magnetic circuit method, which showed that it has obvious advantages compared with other speed regulation technologies. However, at present, most studies on characteristics and parameter optimization are based on finite element software or a large number of experiments and experiences. Basic research is relatively lacking, and analytical methods are rarely used for analysis, design, and optimization.
Currently, the vector magnetic potential method and the equivalent magnetic circuit method are the main numerical simulation methods in engineering. The vector magnetic potential method derives the differential equations about the magnetic vector potential from the Maxwell equations. Then, the finite element or boundary element method is used to solve the problem numerically. When calculating the external magnetic field distribution, it is necessary to know the permeability or magnetization curve of the permanent magnet in the magnetization direction. However, it is difficult to determine these parameters. The magnetic circuit method converts the magnetic field’s calculation into the magnetic circuit’s calculation by reducing the magnetic field to the magnetic circuit equivalent to a circuit in an electric field. The calculation process is simple, and the accuracy is not high. In contrast, the molecular current method concentrates the bound current on the surface of the permanent magnet. It simplifies the model to calculate the spatial distribution of the magnetic field by the equivalent current. This method has simple principles and high precision.
To solve the above problems, the molecular current method was used in this paper to analyze the magnetic induction intensity of the permanent magnet governor12. And the accuracy of the analytical solution was verified by comparing it with the finite element simulation results, and the distribution law of the magnetic field was analyzed. Based on the analytical formula, the influence of the change of permanent magnet parameters on the permanent magnet governor was analyzed, and the influence law of the parameters of the permanent magnet on the magnetic induction intensity was obtained.
The permanent magnet governor is mainly consisting of a permanent magnet rotor, conductor rotor, and regulating mechanism as shown in Fig. 1. It mainly includes the active shaft, load shaft, permanent magnet rotor, conductor rotor, yoke iron, and air gap adjustment mechanism. Its operation principle is that there is a speed difference between the permanent magnet and the conductor rotor. When the two parts move relative to each other, the conductor rotor cuts the magnetic induction lines to generate eddy currents. The induced magnetic field generated by the eddy currents and the magnetic field generated by the permanent magnets are coupled to each other to transmit torque and drive load as shown in Fig. 213.
The physical picture of permanent magnet governor.
Principle model of permanent magnet governor.
Figure 3 shows the molecular circulating current model of permanent magnets. Assuming that the sector permanent magnet is fully and uniformly magnetized in the axial direction and reaches a saturation state. According to Ampere’s law of molecular circulation, macroscopically, the permanent magnet only has surface current but no body current. The magnetic induction intensity generated by the current loop is considered the superposition of the magnetic induction generated by each current element. The single-turn coil trajectory is continuous but not smooth and cannot be integrated continuously, so the equivalent current coil is divided into a superposition of four smooth current segments.
In Fig. 3, R1 and R2 are the inner and outer diameters of the permanent magnet, respectively; h is the height of the permanent magnet; φ is the angle of the permanent magnet; \(P\left( {x,y,z} \right)\) or \(\left( {\rho ,\theta ,z} \right)\) is any point in space, which is denoted as the field point \(G\left( {x_{0} ,y_{0} ,z_{0} } \right)\) or \(\left( {\rho_{0} ,\theta_{0} ,z_{0} } \right)\) is any point on the permanent magnet, which is denoted as the source point; The current element at the source point G can be determined by Biot Savart law and the current intensity \({\text{I}} = {\text{J}}dz_{0}\) .
According to the calculation of the current element magnetic field provided by Biot Savart law, the magnetic induction of the current ring at any point P in space is following.
where, \(\mu_{0}\) is the vacuum permeability, dl is the microlinear element, l is the integral path, \(r\) is the vector diameter from the source to the field point P, and V is the distribution space volume.
The magnetic field B is generated by a permanent magnet at any point in space.
where, dBx, dBy and dBz are the magnetic field components generated by the current element at point P in three directions, respectively.
There is a relationship of \(\rho_{0} = R_{2}\) among r and \(\left| {\varvec{r}} \right|\) .
Substituting the above equations to obtain the magnetic induction produced by the current element at point P.
There is a relationship of \(\theta_{0} = \varphi\) among r and \(\left| {\varvec{r}} \right|\) .
Substituting the above equations to obtain the magnetic induction produced by the current element at point P.
By analogy, the magnetic induction intensity of CD and DA segments can be obtained. According to the current direction superposition and integration of magnetization direction z. The total magnetic fields generated by permanent magnet ABCD at any point P(x, y, z) in its external space can be obtained. Among them,
where, M is the magnetization, n is the unit vector of the outer normal on the side surface of the magnet.
From the above equations, it can be seen that the magnetic induction intensity distribution of the sector permanent magnet is related to the permanent magnet parameters, spatial location, and equivalent current density. And the expression can be expressed as,
where, R1 reflects the radial position of the permanent magnet; R2-R1 can be expressed as the radial length of the permanent magnet. The angle φ of the permanent magnet can be converted into a magnetic pole number.
In a permanent magnet governor, n permanent magnets of the same size and opposite polarity are arranged in a circle (n is an even number, n = 2k). The magnetic induction of the n permanent magnets at point P can be obtained by the following coordinate transformation.
where, m represents the m-th permanent magnet and α is the angle between the two permanent magnets.
After that, the magnetic field distribution in the external space of each permanent magnet is obtained. The analytical equation of the magnetic field generated by n permanent magnets at any point in their external space is obtained by superposition. If two adjacent permanent magnets have opposite polarity, then the three components of the magnetic field are expressed as follows.
where, B(m) is the magnetic induction intensity generated by the m-th permanent magnet.
Previous experimental studies have shown that the numerical solution of the finite element method is the more accurate method when calculating the magnetic induction intensity of permanent magnets13,14. To verify the correctness of the analytical formula derived by the molecular current method, and to analyze the effect of each parameter of the permanent magnet on the magnetic induction intensity, the analytical solution of the calculation is compared with the numerical solution of the Maxwell finite element method through an arithmetic analysis. The parameters of the sector permanent magnets are shown in Table 1.
According to the parameters in Table 1, a three-dimensional model is shown in Fig. 4. Based on the above theoretical derivation, the axial components of path one and path two and the radial component of path three have been calculated by the integral2 function in MATLAB, respectively. The results are shown in Figs. 5, 6, and 7.
Spatial coordinate system of sector permanent magnet.
Axial component variation curve of path one.
Axial component variation curve of path two.
Radial component variation curve of path three.
From Figs. 5 and 6 with the increase of the distance from the surface of the permanent magnet, the trend of the axial component of the magnetic induction intensity is as follows: firstly, the axial component of the magnetic induction intensity is maximum at the end face of the permanent magnet. It shows a decaying trend with the increase of the distance from the surface of the permanent magnet. It indicates that the air gap size is the main reason for the mechanical characteristics of the Permanent magnets15. Secondly, the distribution of the axial component of the air-gap magnetic induction intensity in the same plane transversely has a special relationship with the distance. At a distance less than a specific value, the value of magnetic induction intensity at the central region is lower than that of the edge region.
From Fig. 7, the radial direction, the radial component of the air gap flux density decreases gradually as the distance from the inner and outer sides of the permanent magnet increases.
As seen from the figures, the calculation results of the two methods are consistent, and the maximum deviation point is at the end of the permanent magnet, which is only 6.06%. The analysis method is accurate and effective.
This paper takes the above permanent magnet parameters as an example in the permanent magnet governor, and uses the analytical formula and Maxwell software to study the influence of the permanent magnet structure parameters on its air gap magnetic induction strength, and analyzes the selection principle of the main parameters of the permanent magnet of the disk permanent magnet governor, and summarizes its regular conclusions, and provides a basis for the optimal design of the permanent magnet governor. Since the axial component of the magnetic field has a direct impact on the transmission performance, the radial component is minor, and the distribution changes are more complicated, so the calculation is carried out with the axial component.
To comprehensively investigate the relationship between parameter and axial magnetic density, the structural parameters are radial position, radial length, thickness, and number of magnetic pole pairs. The middle radius of the permanent magnet and the air gap of 1 mm is chosen as the reference point. To study the variation law of each structural parameter on the axial component of the magnetic induction intensity and to provide a basis for structural optimization.
From Fig. 8, when the cross-sectional area and thickness of the permanent magnet remain unchanged, the magnetic induction intensity of the air gap does not change significantly when the radial position of the permanent magnet is changed. The reason is that changes in the shape and position of the section do not affect the magnetic potential. Therefore, the permanent magnets with different cross-sectional shapes are equal. The principle is that the parameters, such as the cross-sectional area and thickness of the permanent magnets, remain unchanged before and after equating. The transmission performance of the governor has little influence16,17.
Effect of radial position on the axial component of the air gap flux density.
From Fig. 9, when the radial length of the permanent magnet is around 9 mm, the air gap flux density reaches its maximum. As the radial length of the permanent magnet increases, the air gap flux density gradually decreases. The radial length of the permanent magnet mainly affects the radial component of the air-gap magnetic inductance, and the influence on the tangential and axial components is relatively small. Therefore, the air gap magnetic induction can be increased by reducing the radial length.
Effect of radial length on the axial component of the air gap flux density.
From Fig. 10, with the increase of permanent magnet thickness, the air gap flux density gradually increases, and the gradient of the rising curve gradually decreases. Finally, the air gap flux density tends to be stable. The reason is that the magnetoresistance and magneto dynamic potential of the permanent magnet simultaneously increases, which decreases the air gap flux, and the increase rate of the air gap flux density decreases. Therefore, the thickness needs to be selected reasonably to utilize the magnetic energy effectively.
Effect of thickness on the axial component of the air gap flux density.
The influence law of the number of magnetic pole pairs on the magnetic induction intensity is shown in Fig. 11. When the number of pole pairs is eight, the air gap flux density is maximum. Because when the duty cycle is constant, the number of pole pairs increases, the angle of each permanent magnet decreases, and the magnetic induction lines become more concentrated, which enhances the air gap flux density.
Effect of pole pairs number on the axial component of the air gap flux density.
Considering the permanent magnet number and the air gap flux density, the optimized structure parameters of permanent magnet governor are: inner diameter \(R_{1}\) = 26 mm, outer diameter \(R_{2}\) = 35 mm, height h = 12 mm, angle φ = \(\pi /12\) . The optimized structure is analyzed and calculated, and the results before and after optimization are compared (the permanent magnet performance parameters remain unchanged), the comparison results are shown in Table 2. The amount of permanent magnet before the optimization is 14,667 mm3, and the air gap flux density is 0.357 T; the amount of permanent magnet after optimization is 13,798 mm3, and the air gap flux density is 0.408 T. The amount of permanent magnet after optimization is reduced by 7.84%, and the axial component of air gap flux density is increased by 14.3%.
The sector permanent magnets are mounted for measurement tests. Figure 12 shows the test device. To measure the magnetic induction intensity of two sizes of permanent magnets in different spatial positions before and after optimization, the position of the measurement point was changed by adjusting the height of the elevating table. The magnetic induction intensity of permanent magnets was measured several times in different spatial positions, and the average value was recorded. The permanent magnet is mounted on the permanent magnet holder and clamped by an aluminum bench vice, and then the Tesla meter probe is fixed on the lifting table. In the test, a cursor universal angle ruler and a vernier caliper are used to ensure that the Tesla meter probe is aligned with the center of the permanent magnet, and a vernier caliper measures the distance between the Tesla meter probe and the upper surface of the permanent magnet. The magnitude of the magnetic induction intensity at the current position can be read from the handheld Tesla meter. By moving the platform height, the magnetic induction intensity of the sector permanent magnets is obtained in different spatial positions before and after optimization.
Measurement test of magnetic induction strength of permanent magnets. ① lifting platform ② aluminum bench vice ③ permanent magnet ④ permanent magnet holder ⑤ vernier caliper ⑥ cursor universal angle ruler ⑦ handheld Tesla meter ⑧ Tesla meter probe.
Figure 13 shows the variation curve of magnetic induction intensity measured by the sector permanent magnet at different heights obtained from the above experimental process.
Experimental results comparison before and after optimization.
The test results show that the magnetic induction intensity decreases with the in-crease in distance, and there is an intersection point at a distance of 2.5 mm. When the distance from the upper surface is less than 2.5 mm, the magnetic induction intensity after optimization is more significant than that before optimization. When the distance from the upper surface is greater than 2.5 mm, the magnetic induction intensity is greater than that after optimization. The optimization effect is most obvious on the upper surface of the permanent magnet, and the magnetic induction intensity after optimization is increased by 27.5% compared with that before optimization. Therefore, the optimized parameters using the analytical formula by the molecular current method can effectively improve the magnetic induction strength of the permanent magnet, which can guide the preliminary design and later product improvement.
This paper calculates the spatial magnetic field distribution of the sector permanent magnet in the permanent magnet governor based on the molecular current method. The correctness of the solution is verified by the numerical analysis and comparison with the finite element method. According to the analytic formula, the following conclusions can be drawn:
Without considering the number of permanent magnets and limiting the outer diameter of the disk-type permanent magnet governor, the magnetic induction intensity of the air gap does not change significantly with the radial position of permanent magnets. The radial length of the permanent magnet mainly affects the radial component of the magnetic induction, while the effect on the tangential and axial components is relatively small. The magnetic induction intensity at the air gap increases with the thickness in-creasing of permanent magnets, and stabilizes after increasing to a specific value. The magnetic density of the air gap increases with the increase in the number of pole pairs.
When considering the number of permanent magnets, it is required not to exceed the number of permanent magnets in the initial structure. After the structure optimization, the axial magnetic field strength is increased by 14.3% compared with the original one. It provides a theoretical reference for optimizing sector permanent magnets in permanent magnet governor.
The datasets analyzed during the current study available from the corresponding author on reasonable request.
Li, G. K., Jia, R. Z. & Ji, Q. S. Magnetic drive and magnetic pump. Journal of the Chinese Society of Rare Earths, vol. S1 534–537 (1994). http://qikan.cqvip.com/Qikan/Article/Detail?id=4001439122.
Davies, E. J. An experimental and theoretical study of eddy-current couplings and brakes. IEEE Trans. Power Appar. Syst. 82, 401–419. https://doi.org/10.1109/TPAS.1963.291434 (1963).
Lubin, T. & Rezzoug, A. Steady-state and transient performance of axial-field eddy-current coupling. IEEE Trans. Industr. Electron. 62, 2287–2296. https://doi.org/10.1109/TIE.2014.2351785 (2015).
Liu, Y. et al. Simulation and analysis of temperature field of new hybrid high-efficiency adjustable permanent magnetic coupler. Small Spec. Electr. Mach. 46, 69–72 (2018).
Wang, L., Jia, Z. Y. & Zhang, L. Investigation on the accurate calculation of the temperature field of permanent magnet governor and the optimization method of heat conduction. Case Stud. Therm. Eng. 13, 100360. https://doi.org/10.1016/j.csite.2018.11.011 (2019).
Seo, S. W., Kim, Y. H., Lee, J. H. & Choi, J. Y. Analytical torque calculation and experimental verification of synchronous permanent magnet couplings with halbach arrays. AIP Adv. 8, 056609. https://doi.org/10.1063/1.5006731 (2018).
Wu, M. et al. Speed governing performance of the permanent magnetic coupler under centrifugal load. J. Mech. Eng. 55, 225–232. https://doi.org/10.3901/JME.2019.14.225 (2019).
Jin, Y. W., Yang, J. W. & Su, Y. F. Analysis of permanent magnet governor magnetic field and torque of impact study. Mach. Design Manuf. 08, 5–8. https://doi.org/10.19356/j.cnki.1001-3997.2017.08.002 (2017).
Li, W. et al. Modeling and performance analysis of an axial-radial combined permanent magnet eddy current coupler. IEEE Access 8, 78367–78377. https://doi.org/10.1109/ACCESS.2020.2990204 (2020).
Zhang, Y. & Li, S. Research on the improvement scheme of permanent magnet governor. Energy Energy Conserv. 1, 166–168. https://doi.org/10.16643/j.cnki.14-1360/td.2019.01.075 (2019).
Yao, J. Energy saving analysis of Permanent Magnet Adjustable Speed technology. Electric Age. 06, 58–60(2021). http://qikan.cqvip.com/Qikan/Article/Detail?id=7104939548
Sun, Z. S., Li, X. Q. & Li, X. N. Research on permanent magnet adjustable speed driver and its application. Mach. Build. Autom. 45, 1–4. https://doi.org/10.19344/j.cnki.issn1671-5276.2016.03.001 (2016).
Zhang, H. B., Qiu, Y. J. & Jiang, S. Y. Analysis of the equivalent surface current model for the permanent magnet bearing by using the integral definition. J. Mech. Eng. 52, 54–59. https://doi.org/10.3901/JME.2016.07.054 (2016).
Choi, J. Y. & Jang, S. M. Analytical magnetic torque calculations and experimental testing of radial flux permanent magnet-type eddy current brakes. J. Appl. Phys. 111, 756–763. https://doi.org/10.1063/1.3672408 (2012).
Zhi, S. Y. & Wu, H. B. Structure parameters optimization design of the cylindrical permanent magnet governor. Electr. Drive. 47, 16–19. https://doi.org/10.19457/j.1001-2095.20171103 (2017).
Wang, J. et al. A general analytical model of permanent magnet eddy current couplings. IEEE Trans. Magn. 50, 1–9. https://doi.org/10.1109/TMAG.2013.2279073 (2013).
Dai, X. Study on the Analytical Model and Torque Characteristics of Permanent Magnet Eddy Current Couplings (Shanghai Jiao Tong University, 2016). https://doi.org/10.27307/d.cnki.gsjtu.2016.001527.
This research was funded by the National Natural Science Foundation of China (No. 12205004) and Industrial Collaborative Innovation Special Fund of Anhui Polytechnic University and Jiujiang District(No. 2022cyxtb8) and Natural Science Foundation of Anhui under Grant (No. 1908085QA25) and Anhui Key Laboratory of Mine Intelligent Equipment and Technology, Anhui University of Science & Technology(No. 201901004) and the Open Research Fund of Anhui Engineering Technology Research Center of Automotive New Technique(No. QCKJ202008) and Key Research and Development Plan of Anhui Province (No. 2022a05020006) and Major science research project of Anhui Province Colleges and Universities (No.2022AH040134).
School of Mechanical Engineering, Anhui Polytechnic University, Wuhu, 241000, China
Manman Xu, Long Shao & Yimin Lu
Wuhu HIT Robot Technology Research Institute Co., Ltd, Wuhu, 241000, China
Harbin Institute of Technology, Harbin, 150001, China
Wuhu Magnetic Wheel Transmission Technology Co., Ltd, Wuhu, 241000, China
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
M.X. and L.S. performed the data analyses, and writing the main manuscript text. Y.L. contributed to the conception of the study. X.W. offered the main onsite testing technology, and prepared Figs. 12, and 13. C.C. carried out the simulation, and prepared Figs. 5, 6 and 7. Y. F. modified Figs. 2, 3 and 4. All authors reviewed the manuscript.
Correspondence to Yimin Lu or Chuqing Cao.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Xu, M., Fu, Y., Shao, L. et al. Optimized design of permanent magnet for disk permanent magnet governor. Sci Rep 13, 19930 (2023). https://doi.org/10.1038/s41598-023-47047-2
DOI: https://doi.org/10.1038/s41598-023-47047-2
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
Linear Motor Magnet Manufacturers Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.