Next Article in Journal
Effects of Wave Orbital Velocity Parameterization on Nearshore Sediment Transport and Decadal Morphodynamics
Previous Article in Journal
Nearshore Wave Transformation Domains from Video Imagery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Performance Analysis of Conical Permanent Magnet Couplings for Underwater Propulsion

1
School of Marine Science and Technology, Northwestern Polytechnical University, Xi’an 710072, China
2
Key Laboratory for Unmanned Underwater Vehicle, Northwestern Polytechnical University, Xi’an 710072, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2019, 7(6), 187; https://doi.org/10.3390/jmse7060187
Submission received: 27 May 2019 / Revised: 14 June 2019 / Accepted: 15 June 2019 / Published: 18 June 2019

Abstract

:
Permanent magnet couplings (PMCs) are widely used in underwater propulsion because it can solve the deep-sea sealing problem effectively. In this paper, a new type of conical permanent magnet coupling (CPMC) is proposed, which is able to match the tail shape of the underwater vehicle and make full use of the tail space to increase pull-out torque capability. Based on the three-dimensional finite element method (3D-FEM), the electromagnetic characteristics of an initial model for CPMC are analyzed. In order to facilitate the design and optimization of CPMC, an equivalent three-dimensional (3D) analytical method for the pull-out torque calculation is presented, and its accuracy is verified by comparison with the 3D finite element results. Finally, the influence of design parameters such as half-cone angle, pole pair, pole arc coefficient and permanent magnet thickness on maximum pull-out torque and torque density of CPMC is analyzed, and a preliminary optimization model is obtained.

1. Introduction

The deep sea is a huge treasure that still remains to be explored. With the rapid development of marine technology, the interest in deep sea exploration is increasing both in the scientific and business community [1]. There are various types of deep-sea exploration equipment, of which the most widely used are underwater vehicles with streamlined revolved body, such as UUV, AUV, ROV, etc. [2,3,4,5].
With the increasing working depth of underwater vehicle, the deep-sea sealing problem of the power output shaft has become one of the key technologies that need to be broken in the marine science and technology field. The deep sea environment is characterized by high corrosive, high hydrostatic pressure. Traditional mechanical seals have complex structures and serious mechanical losses. At the same time, it is easy to leak grease, which will damage the sea environment.
In recent years, more and more underwater vehicles have adopted a permanent magnet coupling (PMC) on the power output shaft to achieve non-contact torque transmission, which converts the dynamic seal into the static seal, and greatly improves the stability of the underwater propulsion unit [6,7]. The existing permanent magnet coupling is mainly divided into two forms: axial PMC and radial PMC, of which the radial PMC is generally chosen for underwater propulsion [8,9]. The maximum pull-out torque is one of the most important properties of PMC because it directly determines the power output capability of the propulsion unit. Therefore, how to improve the maximum pull-out torque has always been a research focus of PMC. In paper [10], the effect of pole pairs, inner and outer permanent magnets thickness and outer rotor yoke thickness on the maximum pull-out torque of radial PMC is analyzed by three-dimensional finite element method (3D-FEM). In paper [11], the pull-out torque and magnet volume for three types of iron less cylindrical permanent magnet couplings are optimized using a nonlinear optimization method. Paper [12] optimizes the design of an axial-field magnetic coupling by using multi-objective genetic algorithms. In addition, new magnetic circuit structures are constantly being used, such as the use of single-sided or bilateral Halbach arrays to increase the maximum pull-out torque of PMC [13,14].
The existing PMC with traditional structure does not take into consideration the special needs of the underwater propulsion, and thus has many limitations in performance improvement. For the PMC used in underwater propulsion, its design should be coupled with the geometric structure of the vehicle, so as to maximize its performance. The traditional radial PMC is cylindrical and placed in the tail section, which does not match the shape of the tail and cannot make full use of the space. The line type of underwater vehicle is shown in Figure 1a, and there is a conical segment of the tail. In this paper, a new class of conical permanent magnet coupling (CPMC) for underwater vehicles is proposed, which is mounted in the tail, as shown in Figure 1b. The CPMC can make full use of the tail space to improve pull-out torque capability.
Considering the special structure of the CPMC, the pull-out torque can be calculated by 3D-FEM. However, the 3D-FEM is quite time consuming, not suitable for design and preliminary optimization of CPMC. The pull-out torque calculation methods for axial and radial PMCs have been extensively studied. Paper [15] proposes a three-dimensional analytical method for calculating the pull-out torque of radial PMC. In paper [16], simple analytical formulas are derived for computing the pull-out torque, and its accuracy is verified by 3D-FEM and measurements. Paper [17] presents torque analysis of radial PMC based on analytical magnetic field calculations and a Maxwell stress tensor method. In paper [18], the authors developed a 3D analytical method for calculating the pull-out torque of radial PMC based on the expression of one pair magnet.
In this paper, an equivalent three-dimensional analytical method is proposed, and the effects of design parameters such as half-cone angle, pole pair, pole arc coefficient and permanent magnet thickness on the maximum pull-out torque and pull-out torque density of CPMC are analyzed, which lays the foundation for the design and optimization of CPMC.

2. Configuration of the CPMC

The structure of the proposed CPMC is shown in Figure 2a. Like a traditional PMC, the CPMC is composed of an inner rotor and an outer rotor. The inner rotor is connected with motor output shaft, and the outer rotor drives the propeller. Non-contact torque transmission between the inner and outer rotors is achieved by the magnetic force. The dimensional parameters of CPMC include the front end diameter, axial length, half-cone angle, outer rotor yoke thickness, inner and outer permanent magnet thickness, air gap length and so on. Every parameter is illustrated in Figure 2b.
The CPMC is tapered as a whole and the specific size is determined by the shape of the tail section of the underwater vehicle. The main parameters of an initially designed CPMC for an autonomous underwater vehicle (AUV) are shown in Table 1. Both the inner and outer rotors use high-performance permanent magnet Nd-Fe-B, the remnant flux density is 1.27 T and the relative permeability is 1.1045.

3. Performance Analysis of Proposed CPMC

Numerical simulation is a precise and high-efficient investigation approach of the PMC. Considering the special structure of the CPMC, the 3D-FEM is usually employed to compute its electromagnetic characteristics. All of the numerical simulations are carried out by using electromagnetic field analysis software Infolytica/MagNet. A cylindrical area of 200 mm in diameter and 150 mm in length is used as the computational domain for the CPMC. The CPMC is placed at the center of the computational domain. In order to ensure the calculation accuracy, mesh refinement is carried out on the air gap, permanent magnet and the yoke of CPMC. The adaptive mesh is generated using the meshing tool in Infolytica/MagNet, as shown in Figure 3. The maximum element size of the air gap, permanent magnets and yokes are 0.5 mm, 2 mm and 2 mm, respectively, and the total number of elements in the volume mesh is 337,100.
By establishing a suitable finite element analysis model, the magnetic field distribution and electromagnetic force of CPMC can be accurately calculated. However, the 3D-FEM is quite time consuming, it is suitable for analyzing a given prototype rather than for design and optimization. As a new type of permanent magnet coupling, it is very important to develop an analytical method to calculate its pull-out torque transmitted by the CPMC.

3.1. Electromagnetic Characteristics of CPMC

By utilizing electromagnetic field analysis software Infolytica/MagNet, the flux distribution of the whole CPMC at the aligned position and the maximum pull-out torque position, can be observed from Figure 4a,b, respectively. It can be seen that the magnetic field saturation mainly occurs at the rotor yoke and near the front end face of CPMC. The main reason is that when the pole arc coefficient is 1, there is a very serious magnetic leakage between poles.
Figure 5 displays the pull-out torque and axial force of the initially designed CPMC. Obviously, the CPMC generates axial force while transmitting torque, which is one of the characteristics of CPMC. When the inner and outer rotor permanent magnets with the same magnetization direction are aligned, the torque is zero, and there is a maximum axial repulsive force. When the inner and outer rotor deflection angle is at half pole angle, the pull-out torque is maximum, and there is still a repulsive force at the axial direction.

3.2. Equivalent Three-Dimensional Analytical Method

3.2.1. Model Assumptions

To obtain an equivalent three-dimensional analytical solution of the torque transmitted by CPMC, the following assumptions have to be made:
  • The permanent magnet is uniformly magnetized.
  • The yoke material (steel 10) is infinite permeable ( μ r ), and no magnetization saturation occurred.
  • End-effect and magnetic flux leakage are neglected.
  • Ignoring the influence of isolation shell.
  • The permeability of all materials is assumed to be isotropic and homogenous.

3.2.2. Equivalent Model

The CPMC can’t directly use a two-dimensional model to simplify the calculation as a traditional cylindrical permanent magnet coupling. However, the three-dimensional analytical model is very complicated because the permanent magnets of the CPMC are irregular in shape. In the calculation of some performances of permanent magnet motors, a multi-section two-dimensional model is used to approximate the three-dimensional model [19,20,21]. Similarly, we establish an equivalent 3D analysis model here by axial segmentation, as shown in Figure 6. The CPMC is divided into N modules along the axial direction, and each module is a common cylindrical permanent magnet coupling. The pull-out torque of the CPMC is approximately equal to the sum of the torque generated by the N modules.

3.2.3. Torque Calculation

The analytical solution of the torque of permanent magnet couplings has been extensively studied For a common cylindrical permanent magnet coupling, the torque can be calculated by Maxwell equation [8,14]:
T n = L n μ 0 r n 2 B r n B θ n d θ n = 1 , 2 , 3 , N
where, n means the nth module, θ is the deflection angle of internal and external rotors, μ 0 is the permeability of vacuum, r m is the air gap radius, L n is the axial length, B r n and B θ n are the radial and tangential components of the air gap flux density of the nth module, respectively.
When N is large enough, the torque of CPMC can be expressed as follows:
T a l l n = 1 N T n
The dimension parameters of cylindrical permanent magnet coupling are shown in Figure 7. The air gap flux density is jointly excited by the inner and outer rotor permanent magnets, so it can be expressed as:
B r n r , θ = B r i r , θ + B r o r , θ
B θ n r , θ = B θ i r , θ + B θ o r , θ
where, B r i r , θ , B r o r , θ , B θ i r , θ , B θ o r , θ are the radial and tangential components of flux density of the air gap flux density when the inner and outer rotors are separately excited. Each flux density component can be calculated through an analytical method [22]. When the magnetization mode of the permanent magnets is different, the analytical expression of the magnetization components M r and M θ are different, the waveform is shown in Figure 8.
The complete solution for the flux density components in the air gap can be deduced from the general solution of Laplacian/quasi-Poissonian equations with specified boundary conditions.
B r i r , θ = k = 1 , 3 , 5 K B i k f B r i r cos k p θ
B r o r , θ = k = 1 , 3 , 5 K B o k f B r o r cos k p θ + ϕ
B θ i r , θ = k = 1 , 3 , 5 K B i k f B θ i r cos k p θ
B θ o r , θ = k = 1 , 3 , 5 K B o k f B θ o r cos k p θ + ϕ
where,
K B i k = μ 0 M k i μ r i k p k p 2 1 × A 3 k i 1 + 2 R r i R m i k p + 1 A 3 k i + 1 R r i R m i 2 k p μ r i + 1 μ r i 1 R r i R s i 2 k p μ r i 1 μ r i R m i R s i 2 k p R r i R m i 2 k p
K B o k = μ 0 M k o μ r o k p k p 2 1 × A 3 k o 1 + 2 R r o R m o k p + 1 A 3 k o + 1 R r o R m o 2 k p μ r o + 1 μ r o 1 R r o R s o 2 k p μ r o 1 μ r o R m o R s o 2 k p R r o R m o 2 k p
f B r i r = r R s i k p 1 R m i R s i k p + 1 + R m i r k p + 1
f B θ i r = r R s i k p 1 R m i R s i k p + 1 + R m i r k p + 1
f B r o r = r R s o k p 1 R m o R s o k p + 1 + R m o r k p + 1
f B θ o r = r R s o k p 1 R m o R s o k p + 1 + R m o r k p + 1
For radial magnetized,
M r n = 2 B r μ 0 α p sin k π α p 2 / k π α p 2 M θ n = 0 M n i = M n o = M r n + n p M θ A 3 n i = A 3 n o = n p
The size parameters are as follows:
R r i = R 1 , R m i = R 2 , R s i = R 4         inner   rotor R r o = R 4 , R m o = R 3 , R s o = R 1     outer   rotor

3.2.4. Comparison with 3D-FEM

The CPMC is divided into N modules and the number of N affects the analytical calculation results. When N changes from 5 to 600, the maximum pull-out torque of the initially designed CPMC calculated by analytical method is shown in Figure 9. As can be seen from the figure that N of 200 and 600 gave substantially the same results. Considering the accuracy and time economy, set N equal to 200 for the following calculations.
The proposed analytical method and 3D-FEM are used to calculate the pull-out torque of CPMC at different deflection angles, and the pull-out torque curve of half cycle is obtained, as shown in Figure 10. It can be seen that the torque curve obtained by the two methods have the same trend, and the results calculated by the equivalent three-dimensional analytical method are slightly larger than the 3D-FEM results. When the inner and outer rotor deflection angle is 15°, the maximum pull-out torque is obtained. The relative error between the calculation results of the two methods is 7.4%, which proves the validity and accuracy of the analytical method. The main reason for the error is that the equivalent multi-section method ignores the effects of yoke magnetic field saturation, end effect and magnetic flux leakage.

4. Design Consideration and Optimization

In engineering applications, we most concerned about the maximum transmission torque ( T max ), because it directly determines the power output capability of the propulsion unit. Since permanent magnets are relatively expensive, it is also necessary to increase the T max generated by per unit volume of permanent magnets and define the torque density S T , as follows:
S T = T max / V p
where, V p is the volume of the permanent magnets of CPMC. Considering that the S T is very meaningful for making full use of permanent magnets and reducing the cost of CPMC. In the next study, the influence of design parameters on T max and S T of the initially designed CPMC will be analyzed by equivalent three-dimensional analytical method.

4.1. Influence of Half-Cone Angle

Keeping the diameter of the front end face and axial length unchanged, the variation of the T max and S T of the CPMC with the half-cone angle is studied by analytical method, as shown in Figure 11. It can be seen that as the half-cone angle increases, the T max and S T are continuously reduced. In practical applications, the size of the half-cone angle is often determined by the shape of the tail section of the underwater vehicle.

4.2. Influence of Pole Pairs

Keeping the size parameters constant, the T max and S T of the CPMC versus different pole pairs are calculated with the pole pair increasing from 5 to 10. As shown in Figure 12, both the T max and S T increase first and then decrease as the number of polo pairs increases, and have a maximum value when the pole pair is 8.

4.3. Influence of Pole Arc Coefficient

The change of pole arc coefficient means a change of the volume of the permanent magnet, and a larger volume of permanent magnets often means a bigger pull-out torque. Figure 13 shows the variation of T max and S T versus the pole arc coefficient. The T max increases with the increase of the pole arc coefficient, and the S T increases first and then decreases as the pole arc coefficient increases, and have a maximum value when the pole arc coefficient is 0.7.

4.4. Influence of Permanent Magnet Thickness

Assuming that the thickness of internal and external permanent magnets is the same, the thickness of rotor yoke is at least 1.5 times that of permanent magnets. Keeping other parameters unchanged, and only changing the thickness of permanent magnet, the variation of T max and S T are shown in Figure 14. The maximum value of T max and S T are obtained with different thicknesses.

4.5. Interaction Between the Pole Pair, Pole Arc Coefficient and Permanent Magnet Thickness

The half-cone angle of the CPMC designed for underwater propulsion is 12°, which is determined by the vehicle’s profile. The analytical method can be used to analyze the interaction between the pole pairs, pole arc coefficient and permanent magnet thickness (as shown in Figure 15, Figure 16, Figure 17, Figure 18, Figure 19 and Figure 20), which is the basis for the optimal design of CPMC. The pole pair number is a discrete value within a limited range and determining the pole pair in advance can simplify the optimization process.
It can be seen from the figures below that under different pole arc coefficients, the CPMC with a pole pair of 8 always has the maximum T max and S T . The CPMC with a pole pair of 8 also has good performance at different permanent magnet thickness. Hence, for the CPMC studied in this paper, the pole pairs be 8 is more suitable for further study. Then, the CPMC optimization design problem can be simplified to an optimization problem of two independent variables the permanent magnet thickness h m and the pole arc coefficient α .

4.6. Preliminary Optimization Design

In this paper, the non-dominated sorting genetic algorithm-II (NSGA-II) was utilized to solve the multi-objective optimization of CPMC. NSGA-II is one of the most widely used multi-objective optimization algorithms, which has the characteristics of high efficiency, minimum user interaction, solutions uniformly distributed and so on [23,24].
Generally, NSGA-II is a method of searching for the Pareto optimal solutions according to several certain known models (surrogate models). In this paper, these surrogate models are established by using Kriging method, such as maximal pull-out torque ( T max ) and torque density ( S T ). The constraint conditions were given based on the actual size of the CPMC. The multi-objective function could be defined as:
V : F X = M A X T max X , S T X s . t . 1.5 h m 3 0.65 α 1.0 X = h m , α T
NSGA-II could find a set of Pareto optimal solution and corresponding design parameters of the CPMC.
The h m changes from 1.5 mm to 3 mm with a 0.5 mm interval and the pole arc coefficient α changes from 0.65 to 1.0 with a 0.05 interval. 40 samples were analyzed by using the equivalent three-dimensional analytical method. The objective T max and S T is shown in Table 2 and Table 3. The response surfaces of two sub objective function are built as shown respectively in Figure 21 and Figure 22.
The optimization procedure uses 100 individuals evolving during 50 generations, the mutation probability is 0.05, and the recombination probability is 0.5. Figure 23 shows the Pareto front of the maximum pull-out torque and torque density of CPMC. The points A and C are the optimum design points which can be obtained by single-objective optimization. For the CPMC, it is necessary to achieve a given pull-out torque, while making the torque density as large as possible. Setting the pull-out torque to 14 Nm, the optimum design points is marked with the red five-pointed star B. The optimum design parameters and objective values obtained by NSGA-II are shown in Table 4. The pull-out torque of the preliminary optimized model is improved by 17.545% and the torque density increased from 5 . 101   ×   10 5   N m 2 to 5 . 997   ×   10 5   N m 2 when compared with that of the initially designed CPMC.

5. Conclusions

A special CPMC for underwater propulsion is proposed and analyzed in this paper. The electromagnetic characteristics are analyzed by 3D-FEM at first. Then, an equivalent three-dimensional analytical method used to compute the pull-out torque of CPMC has been derived based on a multi-section two-dimensional model. The accuracy of proposed analytical method is verified by comparing with the 3D finite element results. Based on the analytical method, the influence of design parameters such as half-cone angle, pole pair, pole arc coefficient and permanent magnet thickness on maximum torque and torque density of CPMC are analyzed, and a preliminary optimization model is obtained. Further work will focus on the experimental research of the CPMC, as well as the multi-level and multi-disciplinary optimization of the overall performance of CPMC. Furthermore, since the flux density of the CPMC yoke is unevenly distributed in the axial direction, the thickness of the yoke could be unequal in the axial direction or a Halbach array could be used to reduce the volume of the yoke.

Author Contributions

Conceptualization, Y.L., and Y.H.; methodology, Y.H. and Y.L.; software, Y.L. and Z.M.; validation, Y.H., B.S. and Z.M.; formal analysis, Y.L. and W.T.; writing—original draft preparation, Y.L. and Z.M.; writing—review and editing, B.S.; supervision, Z.M. and W.T.; project administration, Z.M. and W.T.; funding acquisition, Y.H., Z.M. and B.S.

Funding

This research was funded by the National Science Foundation of China (Grant Nos. 51179159, 61572404) and the Shaanxi Province Youth Science and Technology New Star Project (Grant No. 2016KJXX-57).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, Z.H.; Nie, S.L.; Zhang, L.M.; Yuan, S.H. Development of seawater hydraulic pump tester in deep-sea simulated environment. In Proceedings of the 2015 International Conference on Fluid Power and Mechatronics, Harbin, China, 5–7 August 2015; pp. 667–671. [Google Scholar]
  2. Byron, J.; Tyce, R. Designing a vertical/horizontal AUV for deep ocean sampling. In Proceedings of the OCEANS 2007, Vancouver, BC, Canada, 29 September–4 October 2007; pp. 1–10. [Google Scholar]
  3. Fletcher, B. UUV master plan: A vision for navy UUV development. In Proceedings of the OCEANS 2000 MTS/IEEE Conference and Exhibition, Providence, RI, USA, 11–14 September 2000; pp. 65–71. [Google Scholar]
  4. Zhang, B.S.; Song, B.W.; Mao, Z.Y.; Jiang, J.; Sun, C.Y. Hydrodynamic characteristics and stability simulation of four-rotor dish-shaped UUV landing on the seabed. In Proceedings of the OCEANS 2016 MTS/IEEE Monterey, Monterey, CA, USA, 19–23 September 2016; pp. 1–5. [Google Scholar]
  5. Provencher, L.; Sarraillon, S. The MASKI+underwater inspection robot: A new generation ahead. In Proceedings of the 2016 4th International Conference on Applied Robotics for the Power Industry (CARPI), Jinan, China, 11–13 October 2016; pp. 1–6. [Google Scholar]
  6. Charpentier, J.F.; Fadli, N.; Jennane, J. Study of ironless permanent magnet devices being both a coupling and an axial bearing for naval propulsion. IEEE Trans. Magn. 2003, 39, 3235–3237. [Google Scholar] [CrossRef]
  7. Qi, C.P.; Zhao, J.X.; Zheng, Y.Y. Design of Halbach array magnetic coupling for underwater propeller. Micromotors 2017, 50, 80–83. [Google Scholar]
  8. Cheng, B.; Pan, G. Analysis and structure optimization of radial Halbach permanent magnet couplings for deep sea robots. Math. Probl. Eng. 2018, 2018, 7627326. [Google Scholar] [CrossRef]
  9. Yu, Y.L.; Zhao, B. Study on eddy current loss and system transmission efficiency in insulation shell of the under-water permanent magnet shaft coupling. Trans. China Electr. Tech. Soc. 2014, 29, 15–22. [Google Scholar]
  10. Högberg, S.; Mijatovic, N.; Holbøll, J.; Jensen, B.B.; Bendixen, F.B. Parametric design optimization of a novel permanent magnet coupling using finite element analysis. In Proceedings of the 2014 IEEE Energy Conversion Congress and Exposition, Pittsburgh, PA, USA, 14–18 September 2014; pp. 1465–1471. [Google Scholar]
  11. Charpentier, J.F.; Lemarquand, G. Optimal design of cylindrical air-gap synchronous permanent magnet couplings. IEEE Trans. Magn. 1999, 35, 1037–1046. [Google Scholar] [CrossRef]
  12. Dolisy, B.; Mezani, S.; Lubin, T.; Lévêque, J. A new analytical torque formula for axial field permanent magnets coupling. IEEE Trans. Magn. 2015, 30, 1–8. [Google Scholar] [CrossRef]
  13. Yuan, D.Q.; Li, J.P.; He, Y.Q.; Zhang, H. Characteristic analysis of transmission torque of magnetic coupling with Halbach array. J. Magn. Mater. Devices 2011, 42, 36–40. [Google Scholar]
  14. 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. 2018, 8, 0566091–0566096. [Google Scholar] [CrossRef]
  15. Ravaud, R.; Remarquand, V.; Remarquand, G. Analytical design of permanent magnet radial couplings. IEEE Trans. Magn. 2010, 46, 3860–3865. [Google Scholar] [CrossRef]
  16. Lubin, T.; Mezani, S.; Rezzoug, A. Experimental and theoretical analyses of axial magnetic coupling under steady-state and transient operations. IEEE Trans. Ind. Electron. 2014, 61, 4356–4365. [Google Scholar] [CrossRef]
  17. Choi, J.Y.; Shin, H.J.; Jang, S.M.; Lee, S.H. Torque analysis and measurements of cylindrical air-gap synchronous permanent magnet couplings based on analytical magnetic field calculations. IEEE Trans. Magn. 2013, 49, 3921–3924. [Google Scholar] [CrossRef]
  18. Zhu, Z.N.; Meng, Z. 3D torque calculation of magnetic coupling and its characteristic parameters analysis. Electr. Mach. Control 2017, 21, 94–101. [Google Scholar]
  19. Gyselinck, J.J.C.; Vandevelde, L.; Melkebeek, J.A.A. Multi-slice FE modeling of electrical machines with skewed slots-the skew discretization error. IEEE Trans. Magn. 2001, 37, 3233–3237. [Google Scholar] [CrossRef]
  20. Urresty, J.C.; Riba, J.R.; Romeral, L.; Garcia, A. A simple 2-D finite-element geometry for analyzing surface-mounted synchronous machines with skewed rotor magnets. IEEE Trans. Magn. 2010, 46, 3948–3954. [Google Scholar] [CrossRef]
  21. Hong, C.; Huang, W.; Hu, Z. Performance calculation of a dual stator solid rotor axial flux induction motor using the multi-slice and multi-layer method. IEEE Trans. Magn. 2019, 55, 1–9. [Google Scholar] [CrossRef]
  22. Zhu, Z.Q.; Howe, D.; Chan, C.C. Improved analytical model for predicting the magnetic field distribution in brushless permanent-magnet machines. IEEE Trans. Magn. 2002, 38, 229–238. [Google Scholar] [CrossRef] [Green Version]
  23. Vo-Duy, T.; Duong-Gia, D.; Ho-Huu, V.; Vu-Do, H.C.; Nguyen-Thoi, T. Multi-objective optimization of laminated composite beam structures using NSGA-II algorithm. Compos. Struct. 2017, 168, 498–509. [Google Scholar] [CrossRef] [Green Version]
  24. Deb, K.; Pratap, A.; Agarwal, S.; Meyarivan, T. A fast and elitist multiobjective genetic algorithm: NSGA-II. IEEE Trans. Evol. Comput. 2002, 6, 182–197. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Underwater vehicle. (a) Parameters of line type. (b) Installation position of CPMC.
Figure 1. Underwater vehicle. (a) Parameters of line type. (b) Installation position of CPMC.
Jmse 07 00187 g001
Figure 2. Structure of CPM (a) Components (b) Geometrical parameters.
Figure 2. Structure of CPM (a) Components (b) Geometrical parameters.
Jmse 07 00187 g002
Figure 3. Mesh diagram of CPMC.
Figure 3. Mesh diagram of CPMC.
Jmse 07 00187 g003
Figure 4. Flux distribution of the initially designed CPMC obtained by 3D-FEM. (a) At the aligned position. (b) At the maximum pull-out torque position.
Figure 4. Flux distribution of the initially designed CPMC obtained by 3D-FEM. (a) At the aligned position. (b) At the maximum pull-out torque position.
Jmse 07 00187 g004
Figure 5. Pull-out torque and the axial force of the initially designed CPMC obtained by 3D-FEM.
Figure 5. Pull-out torque and the axial force of the initially designed CPMC obtained by 3D-FEM.
Jmse 07 00187 g005
Figure 6. Equivalent 3D analysis model of CPMC.
Figure 6. Equivalent 3D analysis model of CPMC.
Jmse 07 00187 g006
Figure 7. Two-dimensional model of cylindrical permanent magnet coupling.
Figure 7. Two-dimensional model of cylindrical permanent magnet coupling.
Jmse 07 00187 g007
Figure 8. Waveforms of magnetization components M r and M θ . (a) Radial magnetization. (b) Parallel magnetization.
Figure 8. Waveforms of magnetization components M r and M θ . (a) Radial magnetization. (b) Parallel magnetization.
Jmse 07 00187 g008
Figure 9. Influence of N on the analytical calculation results.
Figure 9. Influence of N on the analytical calculation results.
Jmse 07 00187 g009
Figure 10. Comparison of 3D-FEM and analytical method.
Figure 10. Comparison of 3D-FEM and analytical method.
Jmse 07 00187 g010
Figure 11. T max and S T versus half-cone angle.
Figure 11. T max and S T versus half-cone angle.
Jmse 07 00187 g011
Figure 12. T max and S T versus pole pairs.
Figure 12. T max and S T versus pole pairs.
Jmse 07 00187 g012
Figure 13. T max and S T versus pole arc coefficient.
Figure 13. T max and S T versus pole arc coefficient.
Jmse 07 00187 g013
Figure 14. T max and S T versus permanent magnet thickness.
Figure 14. T max and S T versus permanent magnet thickness.
Jmse 07 00187 g014
Figure 15. T max versus pole arc coefficient for several pole pairs.
Figure 15. T max versus pole arc coefficient for several pole pairs.
Jmse 07 00187 g015
Figure 16. S T versus pole arc coefficient for several pole pairs.
Figure 16. S T versus pole arc coefficient for several pole pairs.
Jmse 07 00187 g016
Figure 17. T max versus permanent magnet thickness for several pole pairs.
Figure 17. T max versus permanent magnet thickness for several pole pairs.
Jmse 07 00187 g017
Figure 18. S T versus permanent magnet thickness for several pole pairs.
Figure 18. S T versus permanent magnet thickness for several pole pairs.
Jmse 07 00187 g018
Figure 19. T max versus pole arc coefficient for several permanent magnet thickness.
Figure 19. T max versus pole arc coefficient for several permanent magnet thickness.
Jmse 07 00187 g019
Figure 20. S T versus pole arc coefficient for several permanent magnet thickness.
Figure 20. S T versus pole arc coefficient for several permanent magnet thickness.
Jmse 07 00187 g020
Figure 21. Response surface of sub objective T max .
Figure 21. Response surface of sub objective T max .
Jmse 07 00187 g021
Figure 22. Response surface of sub objective function S T .
Figure 22. Response surface of sub objective function S T .
Jmse 07 00187 g022
Figure 23. Pareto front of two optimization objectives.
Figure 23. Pareto front of two optimization objectives.
Jmse 07 00187 g023
Table 1. Main parameters of the initially designed CPMC.
Table 1. Main parameters of the initially designed CPMC.
Design ParametersValue
Front end diameter D i 60 mm
Axial length L 1 55 mm
Half-cone angle β 12°
Outer rotor yoke thickness h t 3 mm
Outer permanent magnet thickness h m o 2 mm
Air gap length h a 2 mm
Inner permanent magnet thickness h m i 2 mm
Pole pairs p 6
Pole arc coefficient α 1
Table 2. T max of 40 samples.
Table 2. T max of 40 samples.
T max
(Unit: Nm)
Permanent   Magnet   Thickness   h m   ( mm )
1.52.02.53.03.5
α 0.658.99210.69211.18910.7849.832
0.709.77511.64012.19011.75310.718
0.7510.40512.42613.03212.57611.474
0.8010.89413.05413.71513.24712.093
0.8511.25613.53114.24113.76812.579
0.9011.50413.86514.61314.13912.919
0.9511.64914.06314.83514.36113.125
1.011.69614.12914.90914.43413.193
Table 3. S T of 40 samples.
Table 3. S T of 40 samples.
S T   ( Unit :   10 5 N m 2 ) Permanent   Magnet   Thickness   h m   ( mm )
1.52.02.53.03.5
α 0.656.7646.4595.8215.2514.907
0.706.8276.5295.8885.3124.963
0.756.7846.5065.8755.3014.953
0.806.6586.4075.7975.2304.886
0.856.4746.2515.6655.1114.775
0.906.2496.0505.4904.9534.628
0.955.9955.8135.2804.7644.451
1.05.7185.5485.0414.5484.250
Table 4. Optimum design parameters and objective functions.
Table 4. Optimum design parameters and objective functions.
VariableValue
Permanent magnet thickness (mm)2.168
Pole arc coefficient α 0.869
T max (Nm)14.136
S T ( 10 5 N m 2 )5.997

Share and Cite

MDPI and ACS Style

Li, Y.; Hu, Y.; Song, B.; Mao, Z.; Tian, W. Performance Analysis of Conical Permanent Magnet Couplings for Underwater Propulsion. J. Mar. Sci. Eng. 2019, 7, 187. https://doi.org/10.3390/jmse7060187

AMA Style

Li Y, Hu Y, Song B, Mao Z, Tian W. Performance Analysis of Conical Permanent Magnet Couplings for Underwater Propulsion. Journal of Marine Science and Engineering. 2019; 7(6):187. https://doi.org/10.3390/jmse7060187

Chicago/Turabian Style

Li, Yukai, Yuli Hu, Baowei Song, Zhaoyong Mao, and Wenlong Tian. 2019. "Performance Analysis of Conical Permanent Magnet Couplings for Underwater Propulsion" Journal of Marine Science and Engineering 7, no. 6: 187. https://doi.org/10.3390/jmse7060187

APA Style

Li, Y., Hu, Y., Song, B., Mao, Z., & Tian, W. (2019). Performance Analysis of Conical Permanent Magnet Couplings for Underwater Propulsion. Journal of Marine Science and Engineering, 7(6), 187. https://doi.org/10.3390/jmse7060187

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop