With the development of global deepwater oil and gas and offshore wind power resources, there is an increasing demand for ultra-deep water pile hammers (UDPH). In this paper, a permanent magnet synchronous motor (PMSM) with oil fill is proposed for UDPH on the basis of underwater 2500 m depth and 450 kJ impact energy requirements. First, the rated parameters of the PMSM are calculated based on the power required by hydraulic pumps, and then, the method of combining empirical formulas and path algorithms is used to complete the electromagnetic scheme design. Second, finite element analysis of the electromagnetic field is carried out by Maxwell under both no-load and load conditions to verify the rationality of the PMSM. Furthermore, aiming at improving efficiency and reducing cogging, parameter optimization is carried out by a robust parameter optimization design method, and the optimal parameter scheme of the PMSM is obtained. In addition, loss analysis and temperature field coupling calculations of the optimized PMSM are conducted to verify that the temperature of this motor can meet the requirements under different working conditions. Finally, a comparative analysis is conducted between the proposed PMSM and the three-phase asynchronous motor used in the underwater 2,500 m-450 kJ UDPH. The results indicate that the new PMSM introduced in this paper demonstrates significant advantages in terms of weight, volume, and efficiency; moreover, the new PMSM is reasonable and more suitable for UDPH. Therefore, the PMSM proposed in this paper is reasonable and can provide a theoretical basis for improving the traditional UDPH.

Deepwater anchor piles are the foundation for ensuring the development of global offshore oil, gas, and wind energy resources, and an underwater pile hammer is the key to ensuring the installation of deepwater anchor piles. Currently, the development and production technology of pile hammers are basically mature, including the MHU series developed by MENCK Company in Germany, the S series manufactured by IHC Company in the Netherlands, and the MODEL series produced by HPSI Company in the United States.1 However, most of them are small and medium-sized hammers that are used for driving piles above water or on land. Only MENCK and IHC are capable of manufacturing hydraulic piling hammers suitable for underwater applications.2–4 The motors, which are the power source of the hammer driving system, are generally induction motors for both MENCK and IHC. For example, the motors used in underwater hydraulic piling hammers produced by MENCK are three-phase squirrel cage asynchronous motors produced by Foss Company in Germany.5,6 The three-phase asynchronous motor has strong constant speed performance, but this performance also determines the poor speed regulation of the asynchronous motor. Due to the disadvantages of the asynchronous motor, such as its low power factor, large volume, and low efficiency, it is not ideal to apply this motor to driving hydraulic pumps of underwater pile hammers in practical. In addition, the three-phase asynchronous motor used for UDPH is filled with pure water for cooling. Although the cooling effect is good, due to the strong conductivity of water, the insulation requirements of the motor are very strict, which increases the manufacturing cost and the risk of an internal short circuit. As the permanent magnet synchronous motor (PMSM) has the advantages of high power factor, small size, lightweight, high efficiency, low noise, and good dynamic performance, etc., it can replace the traditional three-phase asynchronous motor in the UDPH.7–9 Its characteristics are highly matched with the motor operation requirements under deepwater conditions, and it is the current research focus of high-performance underwater motors all over the world.10 Therefore, it is necessary to study the underwater PMSM applied for UDPH.

Recent advancements in offshore oil and gas exploration have heightened the application of deepwater pile hammer technology in marine engineering. Since the 1990s, Menck has introduced a deepwater pile hammer capable of operating at water depths of up to 300 m. In parallel, IHC has been innovating deepwater pile hammers that utilize water as a medium, replacing hydraulic oil, though there’s no documented usage yet. Research from MENCK highlights the significant role of their MHU-S and -T series of underwater hydraulic hammers in oil and gas exploration. These hammers, equipped with dual-action hydraulic technology and a modular design, achieve an energy transfer efficiency above 95%. This efficiency positions the MHU series as lightweight yet powerful tools for diverse pile configurations.11 Moreover, Ref. 12 offers a comprehensive analysis of subsea structures in deepwater projects, discussing the functionality and depth limitations of hydraulic hammers. The study also notes that, when combined with MENCK's deepwater girdle-type power pack, these hammers can facilitate installations in depths reaching 3000 m. On another note, Pop et al. highlighted the surge in offshore pile driving activities due to the growing number of offshore wind farms in European waters. Such activities emit high-pressure acoustic waves in the surrounding environment, raising concerns about marine life.13 Addressing this, Jiang et al. explored the acoustic response of column vibration, focusing on the noise produced by pile-driving hammers. They proposed a double-shell pipe pile structure, showcasing its efficacy in noise reduction.14 

At present, PMSM is widely used in industrial production, digital control centers, industrial machinery drives, and other fields. It is well known that electric vehicles have been studied since the 1980s, and the PMSM is currently used as the driving motor. Since the release of the Prius I in 1997, Toyota has newly developed the Prius IV multilayer magnetic barrier with built-in PMSM with an output power of 53 kW and a maximum efficiency of 97%. The 2016 BMW i3 high-speed PMSM developed by BMW has an output power of 126 kW and an efficiency of 94%.15 In 2017, Tesla developed the Model 3, whose driving motor is also a PMSM, with a maximum output power of 192 kW and a maximum torque of 430 Nm.16 In order to enhance the power density of high-power electric motorcycle drive motors, a novel PMSM with an optimized asymmetric eccentric air-gap structure is proposed.17 In rail transit, the Japan Railway Comprehensive Technology Research Institute produced the first RMT1 direct drive PMSM in the 1990s, whose weight was reduced by one-third and efficiency was enhanced by 2.7%18 compared with asynchronous motors of the same specification. The direct drive PMSM produced jointly by CRRC Zhuzhou and Siemens has been put into use in the subways in Shenyang, Changsha, and Shenzhen. In the field of the military with its strict requirements, PMSM is still applicable.

Deepwater PMSMs, known for their high efficiency, reliability, and compact size, are increasingly adopted in deep-sea endeavors like offshore oil and gas extraction, submarine operations, and subsea exploration. Table I gives the application of radial PMSM in submarine and destroyer. Gambhire et al. delved into underwater robotic systems, particularly highlighting the potential of sliding mode control in PMSM speed servo systems, offering fresh insights into deepwater motor control.19 Reference 20 discussed the role of PMSM in oil production for deep-well artificial lifts, with a focus on the temperature analysis of a submersible PMSM driving a progressing cavity pump (PCP). Tian et al. explored the integration of PMSM in submarines, emphasizing the dual-position feedback control and drive system during navigation.21 Furthermore, Ref. 22 detailed the use of a PMSM in designing a controlled power system for underwater transmission mechanisms, reflecting humanity's deepening engagement with marine exploration.

TABLE I.

Application of radial PMSM in submarine and destroyer.

ManufacturerRated power (MW)Speed (rpm)Model
BBC 1.5 180 Swedish submarine 
DRS 1.7 120 FFG(X) frigate 
Siemens 1.76 120 U212 submarine 
Siemens 3.9 150 U214 submarine 
DRS 36.5 127 DDG1000 destroyer 
ManufacturerRated power (MW)Speed (rpm)Model
BBC 1.5 180 Swedish submarine 
DRS 1.7 120 FFG(X) frigate 
Siemens 1.76 120 U212 submarine 
Siemens 3.9 150 U214 submarine 
DRS 36.5 127 DDG1000 destroyer 

However, despite the widespread adoption of PMSM in sectors, such as automotive, rail transport, and submarines, there is a noticeable gap in research and practical applications concerning pile hammers, especially in the realm of UDPH.

Based on the above, a novel PMSM for UDPH with an underwater depth of 2500 m and an impact energy of 450 kJ is proposed in this paper. The method of combining field calculation and magnetic circuit calculation is adopted to complete the design of the motor. Then, electromagnetic field finite element analysis and temperature field analysis are conducted to verify the rationality of the new motor.23,24 In addition, aiming at improving efficiency and reducing cogging, parameter optimization is carried out by a robust parameter optimization design method, and the optimal parameter scheme of the PMSM is obtained. Furthermore, in order to verify that the temperature of this motor can meet the requirements, loss analysis and temperature field coupling calculations of the optimized PMSM are conducted. The calculation results of the electromagnetic field and temperature field prove the rationality of the motor designed in this paper, and it meets the application requirements of UPDH. Therefore, this study can provide a new theoretical basis for improving the traditional UDPH and further assist in the development of deepwater oil and gas resources and offshore wind energy. Table II gives the comparison between traditional and new motors.

TABLE II.

Traditional and new motor for UDPH.

Traditional motorNew motor
Motor type Induction motor PMSM 
Internal circulating medium Water Oil 
Weight Heavier Lighter 
Volume Higher Smaller 
Efficiency Lower Higher 
Traditional motorNew motor
Motor type Induction motor PMSM 
Internal circulating medium Water Oil 
Weight Heavier Lighter 
Volume Higher Smaller 
Efficiency Lower Higher 

Currently, only three-phase asynchronous induction motors are utilized for UDPH, with PMSM yet to be implemented. There has not been any significant research on the use of submersible PMSM in UDPH. It is important to note that three-phase asynchronous induction motors and PMSM possess distinct rotor magnetic poles. The former uses induction coils for the rotor’s magnetic poles, while the latter employs permanent magnets. This distinction leads to fundamental differences in their electromagnetic and structural designs. PMSM offers clear benefits, including high magnetic energy, reduced size and weight, and superior efficiency and output torque.

In this study, PMSMs are used to power hydraulic pumps, which, in turn, drive the hammer. Figure 1 depicts the UDPH’s layout.

FIG. 1.

System of UDPH with PMSM.

FIG. 1.

System of UDPH with PMSM.

Close modal

According to the requirements of water depth (2500 m), impact energy(450 kJ), and frequency of the UDPH, the rated parameters of the motor are calculated and determined, as given in Table III.25,26

TABLE III.

Rated parameters of PMSM.

Parameter nameIndicators
Rated power 132 kW 
Number of phase 
Number of pole pairs 
Rated speed 1500 rpm 
Rated torque 840 Nm 
Frequency 50 Hz 
Pressure ≥25 MPa 
Rated voltage 3000 V 
Temperature 3 °C 
Parameter nameIndicators
Rated power 132 kW 
Number of phase 
Number of pole pairs 
Rated speed 1500 rpm 
Rated torque 840 Nm 
Frequency 50 Hz 
Pressure ≥25 MPa 
Rated voltage 3000 V 
Temperature 3 °C 

Compared to conventional ground motors, the PMSM examined in this research faces notable challenges due to its operation at a daunting depth of 2500 m underwater.

  1. Sealing performance: Given the deepwater conditions and seawater’s conductivity, it is crucial to implement reliable sealing measures. This prevents seawater and micro-organisms from infiltrating the motor, which could lead to short circuits or even significant motor damage.

  2. Pressure resistance: Seawater pressure escalates with increasing depth, rising by ∼1 MPa for every 100 m. At 2500 m, the motor encounters a pressure of around 25 MPa. Such immense pressure can compromise delicate structures, making it essential to factor in this pressure during motor design.

  3. Specification requirements: The motor’s capacity and weight are vital considerations, especially concerning energy transmission during UDPH’s hammering process. Enhancing UDPH efficiency entails reducing the motor’s weight and volume without compromising its power.

  4. Temperature management: PMSMs operating in deep-water settings grapple with limited heat dissipation compared to their counterparts in the air. Coupled with their higher power density, this results in more significant losses per unit volume. Consequently, various motor components undergo a pronounced temperature increase. Such elevated temperatures not only affect the motor’s performance and longevity but also threaten its insulation, risking potential motor damage.

  5. Corrosion resistance: The electrochemical nature of seawater, combined with material erosion and marine microbial activity, makes metals susceptible to corrosion. Therefore, it is imperative that the motor’s casing be crafted from corrosion-resistant materials, given its direct exposure to seawater.

Compared to ground motors, the structure of deep-water PMSM may vary based on their application in different underwater situations, primarily to address the specific requirements and challenges posed by the marine environment. These include factors such as pressure and cooling demands. The deep-water permanent magnet motors studied in this article exhibit the following specific structures. Fig. 2 shows the layout of PMSM for UDPH.

  1. Pressure compensator. In order to withstand the high pressure, a pressure compensator is installed at the end of the motor. It is a positive pressure type that maintains the internal pressure of the motor slightly higher than the external seawater pressure to prevent the strong corrosive seawater from pouring into the motor.

  2. Oil-filled systems. In contrast to typical water-filled motors, the pressure protector’s interior and the motor assembly are entirely filled with insulating motor oil. Due to water’s poor insulating properties, higher conductivity, and the potential for increased instability under high pressure, oil is filled instead of water to offer better protection for motor insulation and also lubricate rotating components.

  3. Oil circulation systems. In addition to featuring a pressure protector, at the other end of the motor, there is an impeller installed. The impeller, protector, and oil passages together form an oil circulation system. During the circulation process, the motor oil dissipates heat while also lubricating bearings and other components. The hollow motor shaft and the gap between the stator and rotor are the main passageways for the motor oil.

  4. Temperature rise estimation. Considering the harsh operating environment and limited heat dissipation conditions of deepwater motors, accurate temperature rise estimation is crucial for the construction of the PMSM studied in this study. In the calculation process, the flow of internal motor oil and external seawater should both be carefully taken into account during the calculation process. As described in more detail in the following parts, this study will employ a multi-physics two-way coupled method.

FIG. 2.

Construction of the PMSM for UDPH.

FIG. 2.

Construction of the PMSM for UDPH.

Close modal

According to the requirements of water depth (2500 m), the pressure compensator needs to provide a pressure of no less than 25 MPa and the volume is 20 l. As the vibration-absorption structures are part of UDPH, a detailed design will not be presented in this paper.

First, the basic size parameters of the motor according to the empirical formula, including the pole-slot fit, the inner diameter and length of the stator core, the structure size of the stator and the rotor, etc., are calculated according to the empirical formula.27,28

Then, according to the basic structural parameters calculated above, further calculations are conducted to get the magnetic circuit parameters, including stator tooth width, stator tooth calculated length, stator yoke calculated length and height, rotor yoke calculated length and height, etc. After obtaining the magnetic circuit parameters of the motor, the no-load magnetic circuit of the motor can be calculated.29 The process is shown in Fig. 3.

FIG. 3.

Calculate process of the magnetic circuit.

FIG. 3.

Calculate process of the magnetic circuit.

Close modal
The following is the primary formula for PMSM:
(1)
where Dil is the inner diameter of the stator core, Lef is the length of the stator core, n is the rotation speed, P′ is the calculated power of the motor, A is the electric load, Bδ is the air-gap flux density, αi is the polar arc coefficient, Kdp is the waveform coefficient of airgap magnetic field, and KNm is the armature winding coefficient.
The expression for the airgap magnetic density Bδ is
(2)
where Br is the residual magnetic density, hm is the magnetization direction length of the permanent magnet, μr is the relative permeability of permanent magnets, and σ is the leakage coefficient.
The no-load electromotive force (EMF) is generated by the fundamental magnetic field in the airgap,
(3)
where E0 is the no-load EMF, f is the current frequency, Kdp is the winding factor, N is winding turns, φδ0 is the no-load main magnetic flux, and Kφ is the magnetic flux waveform coefficient of airgap.
The slots and pole combination of PMSM is
(4)
where Q is the number of stator slots, q is the slots per pole per phase, p is the number of pole pairs, and m is the phase number.
Based on experience, the slot width Bs0 for this motor design is 2 mm, the slot height Hs0 is 1 mm, and the slot slope height Hs1 is 1.5 mm. The stator tooth width is
(5)
where t is the armature tooth pitch, Bδ the airgap magnetic flux density, KFe is the stacking coefficient of silicon steel sheets, and Bt is the magnetic flux density of the stator tooth.
The height of the stator yoke is
(6)
where τ1 is the polar pitch, Bj is the magnetic flux density of stator york, and αi is the effective pole arc coefficient.

After the previous calculation and the continuous adjustment of parameters, the basic design scheme of the PMSM is determined as given in Table IV.

TABLE IV.

Design scheme of PMSM.

ParametersValueParametersValue
Stator outer diameter 280 mm Permanent magnet N38 
Stator inner diameter 167 mm Magnet thickness 6 mm 
Rotor outer diameter 159 mm Pole pairs 
Rotor inner diameter 70 mm Slots 30 
Silicon steel sheet DW540 Length of armature 1670 mm 
ParametersValueParametersValue
Stator outer diameter 280 mm Permanent magnet N38 
Stator inner diameter 167 mm Magnet thickness 6 mm 
Rotor outer diameter 159 mm Pole pairs 
Rotor inner diameter 70 mm Slots 30 
Silicon steel sheet DW540 Length of armature 1670 mm 

The finite element simulation model is established to verify the rationality of the electromagnetic design of the motor. The basic guiding ideology of the finite element method is to combine variational principle and subdivision interpolation to solve partial differential equations. The specific implementation process is to divide the region to be solved into finite elements, discretize the continuous problem, and establish and solve the differential equations with each element node as the unknown quantity.30 The finite element analysis model and meshing model of the motor are established by ANSYS Maxwell, as shown in Figs. 4 and 5.

FIG. 4.

Finite element analysis model.

FIG. 4.

Finite element analysis model.

Close modal
FIG. 5.

Meshing of the analysis model.

FIG. 5.

Meshing of the analysis model.

Close modal

In Maxwell 2D, the no-load state boundary conditions and current excitation are set, and the no-load state finite element analysis is carried out to obtain the no-load back EMF, magnetic density cloud map, and radial airgap magnetic density.

It is shown in Fig. 6 that the no-load back EMF waveform of the PMSM is an ideal sinusoidal flat-topped wave,31 and the voltage peak value is lower than the rated voltage of the motor by 3000 V, indicating that the motor can start normally. From the Fourier spectrum analysis diagram of the no-load back EMF, it can be seen that there are fewer harmonics, which meet the design requirements.32 

FIG. 6.

No-load back EMF and Fourier decomposition. (a) Back EMF of three phases under no load. (b) Fourier decomposition of phase A no-load back EMF.

FIG. 6.

No-load back EMF and Fourier decomposition. (a) Back EMF of three phases under no load. (b) Fourier decomposition of phase A no-load back EMF.

Close modal

Figure 7 shows that the airgap flux density of the PMSM designed is an ideal flat-top wave with small fluctuation. It can be seen from the Fourier spectrum decomposition that the harmonic content is less, indicating less fluctuation. The maximum airgap magnetic density is 1.1 T, which meets the design requirements.

FIG. 7.

Flux density and Fourier decomposition. (a) No-load Radial flux density of airgap. (b) Fourier decomposition of airgap radial flux density under no load.

FIG. 7.

Flux density and Fourier decomposition. (a) No-load Radial flux density of airgap. (b) Fourier decomposition of airgap radial flux density under no load.

Close modal

The no-load magnetic flux density and magnetic density vector are shown in Fig. 8, which indicate that the maximum magnetic flux density in the stator teeth is 1.65 T and the maximum magnetic flux density in the stator yoke is 1.43 T, both within a reasonable range.

FIG. 8.

No-load magnetic flux density and magnetic density vector distribution diagram. (a) Magnetic flux density distribution at 0.02 s under no load. (b) Magnetic flux density vector distribution at 0.02 s under no load.

FIG. 8.

No-load magnetic flux density and magnetic density vector distribution diagram. (a) Magnetic flux density distribution at 0.02 s under no load. (b) Magnetic flux density vector distribution at 0.02 s under no load.

Close modal

Through the finite element analysis of the load state of the PMSM, the changes in the output torque, current, magnetic flux density, and magnetic flux density vector under the loading state can be observed to further verify the rationality of the design.33 

From the loading torque curve shown in Fig. 9(a) and the loading three-phase current curve in Fig. 9(b), it can be seen that the operation of the motor is gradually stable and the fluctuation value is small after 60 ms, which meets the design requirements.

FIG. 9.

Loading torque and current. (a) The curve of output torque under loading. (b) Three-phase current curve under loading.

FIG. 9.

Loading torque and current. (a) The curve of output torque under loading. (b) Three-phase current curve under loading.

Close modal

Figures 10(a) and 10(b), respectively, show the magnetic flux density and the magnetic flux density vector at 0.02 s under load. As shown in Fig. 10, when the motor is running under load, the overall distribution of the magnetic flux of the motor is relatively uniform, and there is no magnetic saturation, indicating that there is less magnetic flux leakage. Thus, the motor can maintain good operating performance under load.

FIG. 10.

Loading magnetic flux density and magnetic density vector distribution diagram. (a) Magnetic flux density distribution at 0.02 s under loading. (b) Magnetic flux density vector distribution at 0.02 s under loading.

FIG. 10.

Loading magnetic flux density and magnetic density vector distribution diagram. (a) Magnetic flux density distribution at 0.02 s under loading. (b) Magnetic flux density vector distribution at 0.02 s under loading.

Close modal

The Taguchi robust parameter optimization method can obtain the optimal parameter combination by establishing a reasonable test scheme and reducing test times. The sensitivity of product manufacturing processes to various types of noise factors is reduced through optimal design parameter combinations, and as a result, the quality stability of the product is improved. The idea of this method is to improve product quality through design rather than testing, which can reduce manufacturing costs and achieve ideal products. With considering noise factors, improving efficiency and reducing cogging torque of this motor are set as the optimization objections,34 and ultimately operating performance and stability of this motor are improved.

Controllable factors affecting motor efficiency and cogging torque are selected as follows: stator notch width Bs0, stator slot depth Hs2, polar arc coefficient αp, airgap length δ, and permanent magnet thickness hm, as shown in Fig. 11. Three levels of each controllable factor are taken: initial value, low-level value, and high-level value. The processing tolerance of stator inner diameter ΔDil, airgap Δδ, and permanent magnet thickness Δhm that are prone to errors in the motor manufacturing process are selected as noise factors, and two levels of each noise factor are selected.

FIG. 11.

Diagram of the stator core and slot.

FIG. 11.

Diagram of the stator core and slot.

Close modal

Each controllable factor takes three levels, namely, level 1 (initial value), level 2 (low-level value), and level 3 (high-level value). Then, the controllable factor level table is established as given in Table V.

TABLE V.

Controllable factor level table.

Controllable factorLevel 1 (initial value)Level 2 (low level)Level 3 (high level)
Bs0 (mm) 1.5 2.5 
Hs2 (mm) 18 17 19 
αp 0.75 0.7 0.8 
δ (mm) 0.9 1.1 
hm (mm) 
Controllable factorLevel 1 (initial value)Level 2 (low level)Level 3 (high level)
Bs0 (mm) 1.5 2.5 
Hs2 (mm) 18 17 19 
αp 0.75 0.7 0.8 
δ (mm) 0.9 1.1 
hm (mm) 

A robust parameter optimization orthogonal table with different levels of controllable factors and noise factors was established, and all calculations were completed according to this table. Table VI presents the efficiency optimization orthogonal table and calculation results, and Table VII demonstrates the cogging torque optimization orthogonal table and calculation results, respectively.

TABLE VI.

Robust parameter design efficiency optimization orthogonal table.

Noise factors
Control factors
ΔDil (mm)000.20.2
Δδ (mm)−0.050.05−0.050.05
Δhm (mm)−0.050.050.05−0.05
No.Bs0 (mm)Hs2 (mm)αpδ (mm)hm (mm)η1(%)η2(%)η3(%)η4(%)η̄(%)S/N (dB)
1.5 17 0.7 0.9 95.46 95.47 95.35 95.55 95.46 39.60 
1.5 17 0.7 0.9 94.44 94.46 94.31 95.54 94.69 39.53 
1.5 17 0.7 0.9 93.31 93.33 93.17 93.43 93.31 39.40 
1.5 18 0.75 95.05 95.06 94.92 95.16 95.05 39.56 
1.5 18 0.75 93.88 93.89 93.73 93.99 93.87 39.45 
1.5 18 0.75 92.59 92.61 92.43 92.72 92.58 39.33 
1.5 19 0.8 1.1 94.71 94.73 94.58 94.83 94.71 39.53 
1.5 19 0.8 1.1 93.41 93.43 93.26 93.55 93.41 39.41 
1.5 19 0.8 1.1 92.08 92.01 91.84 92.14 92.02 39.28 
10 17 0.75 1.1 95.16 95.17 95.03 95.27 95.16 39.57 
11 17 0.75 1.1 93.95 93.96 93.81 94.07 93.95 39.46 
12 17 0.75 1.1 92.62 92.64 92.46 92.76 92.62 39.33 
13 18 0.8 0.9 94.41 94.41 94.27 94.53 94.41 39.50 
14 18 0.8 0.9 93.39 93.04 92.87 93.16 93.12 39.38 
15 18 0.8 0.9 91.53 91.53 91.35 91.67 91.52 39.23 
16 19 0.7 95.45 95.47 95.34 95.55 95.45 39.60 
17 19 0.7 94.40 94.42 94.28 94.52 94.41 39.50 
18 19 0.7 93.26 93.28 93.11 93.38 93.26 39.39 
19 2.5 17 0.8 94.54 94.55 94.39 94.67 94.54 39.51 
20 2.5 17 0.8 93.15 93.15 92.98 93.27 93.14 39.38 
21 2.5 17 0.8 91.60 91.61 91.41 91.74 91.59 39.24 
22 2.5 18 0.7 1.1 95.58 95.60 95.47 95.69 95.58 39.61 
23 2.5 18 0.7 1.1 94.52 94.55 94.39 94.65 94.53 39.51 
24 2.5 18 0.7 1.1 93.36 93.39 93.21 93.48 93.36 39.40 
25 2.5 19 0.75 0.9 94.80 94.82 94.67 94.92 94.80 39.54 
26 2.5 19 0.75 0.9 93.55 93.56 93.40 93.67 93.54 39.42 
27 2.5 19 0.75 0.9 92.17 92.18 92.00 92.30 92.16 39.29 
Noise factors
Control factors
ΔDil (mm)000.20.2
Δδ (mm)−0.050.05−0.050.05
Δhm (mm)−0.050.050.05−0.05
No.Bs0 (mm)Hs2 (mm)αpδ (mm)hm (mm)η1(%)η2(%)η3(%)η4(%)η̄(%)S/N (dB)
1.5 17 0.7 0.9 95.46 95.47 95.35 95.55 95.46 39.60 
1.5 17 0.7 0.9 94.44 94.46 94.31 95.54 94.69 39.53 
1.5 17 0.7 0.9 93.31 93.33 93.17 93.43 93.31 39.40 
1.5 18 0.75 95.05 95.06 94.92 95.16 95.05 39.56 
1.5 18 0.75 93.88 93.89 93.73 93.99 93.87 39.45 
1.5 18 0.75 92.59 92.61 92.43 92.72 92.58 39.33 
1.5 19 0.8 1.1 94.71 94.73 94.58 94.83 94.71 39.53 
1.5 19 0.8 1.1 93.41 93.43 93.26 93.55 93.41 39.41 
1.5 19 0.8 1.1 92.08 92.01 91.84 92.14 92.02 39.28 
10 17 0.75 1.1 95.16 95.17 95.03 95.27 95.16 39.57 
11 17 0.75 1.1 93.95 93.96 93.81 94.07 93.95 39.46 
12 17 0.75 1.1 92.62 92.64 92.46 92.76 92.62 39.33 
13 18 0.8 0.9 94.41 94.41 94.27 94.53 94.41 39.50 
14 18 0.8 0.9 93.39 93.04 92.87 93.16 93.12 39.38 
15 18 0.8 0.9 91.53 91.53 91.35 91.67 91.52 39.23 
16 19 0.7 95.45 95.47 95.34 95.55 95.45 39.60 
17 19 0.7 94.40 94.42 94.28 94.52 94.41 39.50 
18 19 0.7 93.26 93.28 93.11 93.38 93.26 39.39 
19 2.5 17 0.8 94.54 94.55 94.39 94.67 94.54 39.51 
20 2.5 17 0.8 93.15 93.15 92.98 93.27 93.14 39.38 
21 2.5 17 0.8 91.60 91.61 91.41 91.74 91.59 39.24 
22 2.5 18 0.7 1.1 95.58 95.60 95.47 95.69 95.58 39.61 
23 2.5 18 0.7 1.1 94.52 94.55 94.39 94.65 94.53 39.51 
24 2.5 18 0.7 1.1 93.36 93.39 93.21 93.48 93.36 39.40 
25 2.5 19 0.75 0.9 94.80 94.82 94.67 94.92 94.80 39.54 
26 2.5 19 0.75 0.9 93.55 93.56 93.40 93.67 93.54 39.42 
27 2.5 19 0.75 0.9 92.17 92.18 92.00 92.30 92.16 39.29 
TABLE VII.

Robust parameter design cogging torque optimization orthogonal table.

Noise factors
Control factors
ΔDil (mm)000.20.2
Δδ (mm)−0.050.05−0.050.05
Δhm (mm)−0.050.050.05−0.05
No.Bs0 (mm)Hs2 (mm)αpδ (mm)hm (mm)Tc1 (N m)Tc2 (N m)Tc3 (N m)Tc4 (N m)Tc̄(N m)S/N (dB)
1.5 17 0.7 0.9 4.41 4.15 4.43 4.15 4.29 −12.64 
1.5 17 0.7 0.9 4.51 4.24 4.53 4.24 4.38 −12.84 
1.5 17 0.7 0.9 4.60 4.32 4.60 4.31 4.46 −12.99 
1.5 18 0.75 3.36 3.09 3.37 3.08 3.23 −10.18 
1.5 18 0.75 3.42 3.15 3.42 3.15 3.29 −10.34 
1.5 18 0.75 3.45 3.18 3.45 3.18 3.32 −10.42 
1.5 19 0.8 1.1 0.00 0.00 0.00 0.00 0.00 52.36 
1.5 19 0.8 1.1 0.00 0.01 0.00 0.01 0.00 46.12 
1.5 19 0.8 1.1 0.01 0.01 0.01 0.01 0.01 41.91 
10 17 0.75 1.1 5.24 4.86 5.26 4.85 5.05 −14.08 
11 17 0.75 1.1 5.34 4.98 5.36 4.98 5.16 −14.27 
12 17 0.75 1.1 5.40 5.03 5.42 5.03 5.22 −14.36 
13 18 0.8 0.9 0.00 0.00 0.00 0.00 0.00 51.40 
14 18 0.8 0.9 0.02 0.01 0.01 0.01 0.01 39.10 
15 18 0.8 0.9 0.01 0.01 0.01 0.01 0.01 40.36 
16 19 0.7 7.08 6.67 7.10 6.66 6.88 −16.75 
17 19 0.7 7.26 6.86 7.28 6.85 7.06 −16.99 
18 19 0.7 7.39 6.98 7.41 6.98 7.19 −17.14 
19 2.5 17 0.8 0.00 0.01 0.00 0.01 0.00 46.14 
20 2.5 17 0.8 0.01 0.01 0.01 0.01 0.01 39.76 
21 2.5 17 0.8 0.02 0.02 0.02 0.02 0.02 35.46 
22 2.5 18 0.7 1.1 10.09 9.47 10.13 9.45 9.79 −19.82 
23 2.5 18 0.7 1.1 10.38 9.75 10.42 9.74 10.07 −20.07 
24 2.5 18 0.7 1.1 10.57 9.94 10.60 9.94 10.26 −20.23 
25 2.5 19 0.75 0.9 8.94 8.43 8.97 8.42 8.69 −18.78 
26 2.5 19 0.75 0.9 9.09 8.60 9.11 8.60 8.85 −18.94 
27 2.5 19 0.75 0.9 9.18 8.68 9.20 8.69 8.94 −19.03 
Noise factors
Control factors
ΔDil (mm)000.20.2
Δδ (mm)−0.050.05−0.050.05
Δhm (mm)−0.050.050.05−0.05
No.Bs0 (mm)Hs2 (mm)αpδ (mm)hm (mm)Tc1 (N m)Tc2 (N m)Tc3 (N m)Tc4 (N m)Tc̄(N m)S/N (dB)
1.5 17 0.7 0.9 4.41 4.15 4.43 4.15 4.29 −12.64 
1.5 17 0.7 0.9 4.51 4.24 4.53 4.24 4.38 −12.84 
1.5 17 0.7 0.9 4.60 4.32 4.60 4.31 4.46 −12.99 
1.5 18 0.75 3.36 3.09 3.37 3.08 3.23 −10.18 
1.5 18 0.75 3.42 3.15 3.42 3.15 3.29 −10.34 
1.5 18 0.75 3.45 3.18 3.45 3.18 3.32 −10.42 
1.5 19 0.8 1.1 0.00 0.00 0.00 0.00 0.00 52.36 
1.5 19 0.8 1.1 0.00 0.01 0.00 0.01 0.00 46.12 
1.5 19 0.8 1.1 0.01 0.01 0.01 0.01 0.01 41.91 
10 17 0.75 1.1 5.24 4.86 5.26 4.85 5.05 −14.08 
11 17 0.75 1.1 5.34 4.98 5.36 4.98 5.16 −14.27 
12 17 0.75 1.1 5.40 5.03 5.42 5.03 5.22 −14.36 
13 18 0.8 0.9 0.00 0.00 0.00 0.00 0.00 51.40 
14 18 0.8 0.9 0.02 0.01 0.01 0.01 0.01 39.10 
15 18 0.8 0.9 0.01 0.01 0.01 0.01 0.01 40.36 
16 19 0.7 7.08 6.67 7.10 6.66 6.88 −16.75 
17 19 0.7 7.26 6.86 7.28 6.85 7.06 −16.99 
18 19 0.7 7.39 6.98 7.41 6.98 7.19 −17.14 
19 2.5 17 0.8 0.00 0.01 0.00 0.01 0.00 46.14 
20 2.5 17 0.8 0.01 0.01 0.01 0.01 0.01 39.76 
21 2.5 17 0.8 0.02 0.02 0.02 0.02 0.02 35.46 
22 2.5 18 0.7 1.1 10.09 9.47 10.13 9.45 9.79 −19.82 
23 2.5 18 0.7 1.1 10.38 9.75 10.42 9.74 10.07 −20.07 
24 2.5 18 0.7 1.1 10.57 9.94 10.60 9.94 10.26 −20.23 
25 2.5 19 0.75 0.9 8.94 8.43 8.97 8.42 8.69 −18.78 
26 2.5 19 0.75 0.9 9.09 8.60 9.11 8.60 8.85 −18.94 
27 2.5 19 0.75 0.9 9.18 8.68 9.20 8.69 8.94 −19.03 

Taguchi regards the signal-to-noise ratio (S/N) as the evaluation index of robust parameter design, and the higher the S/N, the higher the stability of the product. The optimization objective of this time is to improve the motor efficiency and reduce the cogging torque. Obviously, efficiency is usually expected to be higher, while the cogging torque is desired to be lower. Therefore, the optimization objective of efficiency is the larger-the-better characteristic, and conversely, it is the smaller-the-better characteristic for the optimization objective of cogging torque.

The S/N for smaller-the-better characteristic can be expressed as follows:
(7)
The S/N for larger-the-better characteristic can be expressed as follows:
(8)
In Formula (7) and (8), yi is the ith value of the output characteristic.

The average value of the total efficiency S/N and total cogging torque S/N can be obtained by the S/N results given in Tables VI and VII. Table VIII illustrates the average value of total S/N for both efficiency and cogging torque.

TABLE VIII.

Average value of the total efficiency and cogging torque S/N.

NameAverage efficiency S/N Sη̄Average cogging torque S/N STC̄
Value 39.4421 4.1763 
NameAverage efficiency S/N Sη̄Average cogging torque S/N STC̄
Value 39.4421 4.1763 

The mean value of efficiency S/N and efficiency for each controllable factor at different levels are given in Tables IX and X, respectively. The main effect diagrams of efficiency S/N and efficiency are demonstrated in Figs. 12 and 13.

TABLE IX.

Average value of efficiency S/N.

LevelBs0Hs2αpδhm
39.45 39.45 39.50 39.43 39.56 
39.44 39.44 39.44 39.44 39.45 
39.43 39.44 39.38 39.46 39.32 
Delta 0.02 0.01 0.12 0.02 0.23 
Rank 
LevelBs0Hs2αpδhm
39.45 39.45 39.50 39.43 39.56 
39.44 39.44 39.44 39.44 39.45 
39.43 39.44 39.38 39.46 39.32 
Delta 0.02 0.01 0.12 0.02 0.23 
Rank 
TABLE X.

Average value of efficiency.

LevelBs0Hs2αpδhm
93.90 93.83 94.45 93.67 95.02 
93.76 93.78 93.75 93.76 93.85 
93.69 93.75 93.16 93.93 92.49 
Delta 0.21 0.08 1.29 0.26 2.53 
Rank 
LevelBs0Hs2αpδhm
93.90 93.83 94.45 93.67 95.02 
93.76 93.78 93.75 93.76 93.85 
93.69 93.75 93.16 93.93 92.49 
Delta 0.21 0.08 1.29 0.26 2.53 
Rank 
FIG. 12.

Main effect diagram of efficiency S/N.

FIG. 12.

Main effect diagram of efficiency S/N.

Close modal
FIG. 13.

Main effect diagram of efficiency.

FIG. 13.

Main effect diagram of efficiency.

Close modal

The mean value of cogging torque S/N and cogging torque for each controllable factor at different levels are given in Tables XI and XII, respectively. The main effect diagrams of the cogging torque S/N and cogging torque are demonstrated in Figs. 14 and 15.

TABLE XI.

Average value of cogging torque S/N.

LevelBs0Hs2αpδhm
7.8873 4.4655 −16.6067 3.9599 6.4050 
4.1426 4.4231 −14.4881 4.3940 3.5052 
0.4989 3.6402 43.6236 4.1749 2.6186 
Delta 7.3884 0.8253 60.2303 0.4342 3.7864 
Rank 
LevelBs0Hs2αpδhm
7.8873 4.4655 −16.6067 3.9599 6.4050 
4.1426 4.4231 −14.4881 4.3940 3.5052 
0.4989 3.6402 43.6236 4.1749 2.6186 
Delta 7.3884 0.8253 60.2303 0.4342 3.7864 
Rank 
TABLE XII.

Average value of cogging torque.

LevelBs0Hs2αpδhm
2.5519 3.1771 7.1533 4.4026 4.2137 
4.0656 4.4413 5.7489 3.4433 4.3163 
6.2925 5.2913 0.0077 5.0640 4.3798 
Delta 3.7406 2.1144 7.1457 1.6206 0.1661 
Rank 
LevelBs0Hs2αpδhm
2.5519 3.1771 7.1533 4.4026 4.2137 
4.0656 4.4413 5.7489 3.4433 4.3163 
6.2925 5.2913 0.0077 5.0640 4.3798 
Delta 3.7406 2.1144 7.1457 1.6206 0.1661 
Rank 
FIG. 14.

Main effect diagram of cogging torque S/N.

FIG. 14.

Main effect diagram of cogging torque S/N.

Close modal
FIG. 15.

Main effect diagram of cogging torque.

FIG. 15.

Main effect diagram of cogging torque.

Close modal

From Figs. 12 and 13, it can be seen that the controllable factor combinations meeting with the maximum efficiency S/N and efficiency are Bs0(1.5), Hs2(17), αp(0.7), δ(1.1), and hm(5). Figures 14 and 15 show that the controllable factor combinations meeting with the maximum cogging torque S/N and minimum cogging torque are Bs0(1.5), Hs2(17), αp(0.8), δ(1), and hm(5).

Variance analysis can quantitatively analyze the relative importance of each controllable factor to the S/N, thereby determining the optimal combination. Table XIII gives the relative importance that each controllable factor to both efficiency S/N and cogging torque S/N.

TABLE XIII.

Relative importance of controllable factors to S/N.

Controllable factorS/N (η)S/N (Tcog)
SSRelative importance (%)SSRelative importance (%)
Bs0 0.0006 0.54 81.8878 1.15 
Hs2 0.0001 0.07 1.2959 0.02 
αp 0.0215 20.39 7009.1 98.50 
δ 0.0009 0.85 0.2827 0.00 
hm 0.0824 78.15 23.5320 0.33 
TSS 0.1054 100.00 7116.098 100.00 
Controllable factorS/N (η)S/N (Tcog)
SSRelative importance (%)SSRelative importance (%)
Bs0 0.0006 0.54 81.8878 1.15 
Hs2 0.0001 0.07 1.2959 0.02 
αp 0.0215 20.39 7009.1 98.50 
δ 0.0009 0.85 0.2827 0.00 
hm 0.0824 78.15 23.5320 0.33 
TSS 0.1054 100.00 7116.098 100.00 

From Table XIII, it can be seen that the relative importance of Hs2, δ, and hm to the efficiency S/N is larger than that to cogging torque S/N, while the relative importance of Bs0, αp to cogging torque S/N is larger than that to efficiency S/N. The optimal Pareto factor combination for this design is determined as Bs0(1.5), Hs2(17), αp(0.8), δ(1.1), and hm(5), as given in Table XIV. Efficiency and cogging torque are recalculated according to the optimized parameter, and the comparison of results before and after optimization is given in Table XV.

TABLE XIV.

Comparison of controllable factors before and after optimization.

Controllable factorBs0 (mm)Hs2 (mm)αpδ (mm)hm (mm)
Before optimization 18 0.8 1.1 
After optimization 0.75 
Controllable factorBs0 (mm)Hs2 (mm)αpδ (mm)hm (mm)
Before optimization 18 0.8 1.1 
After optimization 0.75 
TABLE XV.

Comparison of efficiency and cogging torque before and after optimization.

EfficiencyCogging torque
Before optimization 93.7 5.56 
After optimization 94.8 0.0024 
Improvement +1.1% −99.9% 
EfficiencyCogging torque
Before optimization 93.7 5.56 
After optimization 94.8 0.0024 
Improvement +1.1% −99.9% 

From Table XV, it can be seen that after this robust parameter optimization, the motor efficiency has been increased by 1.1%, the cogging torque has been reduced by 99.9%, and the motor performance has been greatly improved.

The winding loss, which is also called copper loss, is caused by winding resistance under the action of the current; thus, it is related to the resistance and the current. Different from traditional induction motors, the permanent magnet synchronous motor used for the ultra-deepwater pile hammer only has stator winding loss but no rotor winding loss. Winding loss accounts for a large proportion of losses and is also one of the main sources of motor heating and temperature rise.

According to Joule's law, the calculation formula for motor winding loss is as follows:
(9)
where m is the phase of motors, I1 is the current value, ρ is the resistivity, N denotes the number of conductors, L denotes the copper conductor length, and r denotes the cross-sectional radius of the conductor.

The wingding loss can be calculated by Maxwell, and the loss curve is shown in Fig. 14.

It can be seen from Fig. 16 that the winding loss of the motor tends to be stable after 50 ms, and the loss value is about 3.75 kW.

FIG. 16.

Graph of winding loss.

FIG. 16.

Graph of winding loss.

Close modal

As is known, the resistivity increases with temperature, and consequently, the winding loss will increase with temperature. The influence of temperature on winding loss is considered by the two-way coupled method introduced in detail in our early study,20 and the calculation process is iterative.

The total core loss of PMSM is equal to the sum of loss generated by the fundamental magnetic field and each harmonic wave. Considering the alternating magnetic field, rotating magnetic field, and harmonic magnetic field, the core loss calculation formula is
(10)
where n is the harmonic number, fi is the ith harmonic frequency, Brmi is the radial flux density amplitude of the ith harmonic, Btmi is the tangential flux density amplitude of the ith harmonic, Kh is the hysteresis loss coefficient, Kc is the eddy current loss coefficient, and Ke is the additional loss coefficient. Kh can be obtained by fitting the B–P curve of a silicon steel sheet, and Kc can be obtained by Eq. (11). As the frequency of this motor is low, the core loss is basically composed of hysteresis loss and eddy current loss, and the additional loss can be ignored,
(11)
where d, σ, and ρ are the thickness, conductivity, and density of the silicon steel sheet, respectively.

The core loss curve, hysteresis loss curve, and eddy current loss curve of the motor can be obtained through Maxwell simulation, as shown in Fig. 17. As the additional loss accounts for a small portion of core loss, it can usually be ignored.

FIG. 17.

Graph of core loss.

FIG. 17.

Graph of core loss.

Close modal

Figure 17 indicates that the core loss of the motor tends to be stable after 50 ms and is about 3.1 kW under stable operation, of which the hysteresis loss is about 2.5 kW, accounting for about 80%, and the eddy current loss is about 560 W, accounting for about 20% of the core loss.

Mechanical loss is usually caused by the friction inside the bearing and the friction between rotor rotation and oil. Referring to the experience, 2% of the rated power is selected as the loss value and without a detailed analysis of it. There is eddy current loss generated inside the permanent magnet due to harmonic magnetic fields in the airgap. However, it is very small compared to copper loss and stator core loss so that it has almost no impact on motor temperature.

In order to verify that the temperature rise and output performance of the optimized motor meet the design requirements in a deep water environment of 2500 m, temperature field and electromagnetic field calculations in deep water environments are conducted, respectively, in this section.

Deepwater PMSM, unlike those operating in air, is filled with oil to balance the external pressure of seawater. This method, however, leads to additional fluid friction losses when the rotor rotates. Moreover, the need for deepwater PMSM to be compact and have a high power density results in increased losses per unit volume. Consequently, these significant losses cause a rise in the temperature of the motor. High temperature rise not only affects the performance and service life of the motor but also damages the insulation layer of the motor. Therefore, the accurate temperature rise calculation for these deepwater PMSMs is crucial. Due to the oil inside and seawater outside, the flow of the fluid field has a significant impact on heat dissipation of the motor. It is necessary to analyze the internal fluid field and calculate losses of the motor to employ a two-way coupled method for the accurate temperature rise calculation. Such comprehensive analysis ensures the motor's safe and prolonged operation, while also adhering to temperature rise constraints.

Considering the skin effect of the motor, the influence of temperature on the winding loss,35–37 and the effect of the fundamental and harmonics on the iron loss, the motor loss is calculated.38–41 Through the calculated winding loss and iron loss, the heat generation rate of the motor heat source is further calculated, as given in Table XVI. The thermal conductivity of each thermally conductive component is obtained by consulting relevant materials, as given in Table XVII.15,42,43

TABLE XVI.

Heat generation rate of the motor heat source (W/m3).

Types of lossComponentHeat generation rate (W/m3)
Copper loss Winding 5.3874T + 3724.5705 
  (T is temperature) 
Core loss Stator core 3 × 104 
Mechanical loss Rotor 7 × 103 
Types of lossComponentHeat generation rate (W/m3)
Copper loss Winding 5.3874T + 3724.5705 
  (T is temperature) 
Core loss Stator core 3 × 104 
Mechanical loss Rotor 7 × 103 
TABLE XVII.

Relative thermal conductivity.

ComponentMaterialThermal conductivity (W/m °C)
Housing 45# 45.5 
Stator core DW540-50 16.27 
Stator winding Copper 390 
Lubricating oil 45 transformer oil 0.12 
Magnet N38 7.6 
Rotor core Steel-1010 45 
Shaft 40 Cr 25.1 
ComponentMaterialThermal conductivity (W/m °C)
Housing 45# 45.5 
Stator core DW540-50 16.27 
Stator winding Copper 390 
Lubricating oil 45 transformer oil 0.12 
Magnet N38 7.6 
Rotor core Steel-1010 45 
Shaft 40 Cr 25.1 

We establish the finite element analysis model of the motor and draw the mesh with ICEM, as shown in Fig. 18.44–46 

FIG. 18.

Temperature field calculation model. (a) 3D model for temperature field calculation including PMSM and seawater. (b) Discretization by meshing for the finite element analysis model.

FIG. 18.

Temperature field calculation model. (a) 3D model for temperature field calculation including PMSM and seawater. (b) Discretization by meshing for the finite element analysis model.

Close modal

The 3D model for the temperature field calculation into ANSYS Fluent was imported, and then, we set the boundary conditions.47 Assuming that the initial temperature of the motor and seawater is 3 °C, the seawater around the casing, the motor airgap oil, and the bearing lubricating oil are fluid. The heat generation rate and thermal conductivity were loaded to the corresponding components according to Tables XVI and XVII, respectively.48 

Lubricating oil is mainly in the gap and the shaft internal channel, and the oil circulation route is shown in Fig. 19. The flow rate affects the motor temperature; meanwhile, temperature plays a part in changing the viscosity of the fluid. As the viscosity only affects the mechanical loss that accounts for a small proportion of whole loss, the effect of temperature on fluid is not considered in this paper.

FIG. 19.

Lubricating oil circulation.

FIG. 19.

Lubricating oil circulation.

Close modal

According to Sec. IV of the robust parameter optimization, those factors that affect efficiency and cogging torque are optimized. The new values are given in Table XIV. In order to confirm that the output performance of the motor meets the design requirements under the optimized parameter combination, it is necessary to carry out electromagnetic field finite element calculation. Based on the electromagnetic field analysis of non-optimization, we conduct another calculation for the optimized electromagnetic design. The relevant parameter settings in Ansys Maxwell refer to Tables IV and XIV.

The PMSM designed in this paper is internally oil filled, which is different from the existing way of internal water filling of the asynchronous motor for the deep-water pile hammer. The heat generated in the motor is mainly dissipated by the seawater outside the motor and the lubricating oil inside the motor, both of which are fluids, and their velocity will affect the heat dissipation of the motor. Underwater 2500 m, the seawater is relatively stable. According to the situation that is not conducive to heat dissipation, the seawater flow rate is set to 0.3m/s. In order to study the influence of the lubricating oil flow rate on heat dissipation of the motor, the flow rate was set to 0, 0.1, 0.3, and 0.5 m/s, respectively. The temperature contours of the motor as a whole and each component under different flow rates are obtained through simulation.

Through all the simulation calculations of the PMSM temperature field with different lubricating oil flow rates, all the results of different flow rates are summarized in Table XVIII for comparison.

TABLE XVIII.

Maximum temperature of the four main components at different lubricating oil flow rates.

ComponentV = 0 m/s ( °C)V = 0.1 m/s ( °C)V = 0.3 m/s ( °C)V = 0.5 m/s ( °C)
Winding 67.01 61.99 61.85 61.81 
Stator core 43.07 35.35 35.12 35.09 
Magnet 22.21 5.19 4.70 4.56 
Rotor core 20.82 5.29 4.81 4.67 
ComponentV = 0 m/s ( °C)V = 0.1 m/s ( °C)V = 0.3 m/s ( °C)V = 0.5 m/s ( °C)
Winding 67.01 61.99 61.85 61.81 
Stator core 43.07 35.35 35.12 35.09 
Magnet 22.21 5.19 4.70 4.56 
Rotor core 20.82 5.29 4.81 4.67 

As shown in Figs. 20 and 21 and Table XVIII, the highest temperature part of the motor is the stator winding. When the lubricating oil flow rate is 0 m/s, the winding temperature is 67.01 °C, which is 64 °C higher than the temperature of the seawater around the motor; when the lubricating oil has a flow rate, the temperature is about 62 °C, which decreases by about 7%. The reason for this is that the winding loss of the motor is the largest, and so, the heat generated by the winding is the highest; when the lubricating oil flows, it takes away more heat, and the temperature of winding is lowered; in addition, due to the thermal conductivity of the insulating material around the copper wire is low, heat dissipation of the winding is slow.

FIG. 20.

Temperature contours of the overall motor at different lubricating oil flow rates. (a) Lubricating oil flow rate at 0 m/s. (b) Lubricating oil flow rate at 0.1 m/s. (c) Lubricating oil flow rate at 0.3 m/s. (d) Lubricating oil flow rate at 0.5 m/s.

FIG. 20.

Temperature contours of the overall motor at different lubricating oil flow rates. (a) Lubricating oil flow rate at 0 m/s. (b) Lubricating oil flow rate at 0.1 m/s. (c) Lubricating oil flow rate at 0.3 m/s. (d) Lubricating oil flow rate at 0.5 m/s.

Close modal
FIG. 21.

Temperature contours of the winding at different lubricating oil flow rates. (a) Lubricating oil flow rate at 0 m/s. (b) Lubricating oil flow rate at 0.1 m/s. (c) Lubricating oil flow rate at 0.3 m/s. (d) Lubricating oil flow rate at 0.5 m/s.

FIG. 21.

Temperature contours of the winding at different lubricating oil flow rates. (a) Lubricating oil flow rate at 0 m/s. (b) Lubricating oil flow rate at 0.1 m/s. (c) Lubricating oil flow rate at 0.3 m/s. (d) Lubricating oil flow rate at 0.5 m/s.

Close modal

Figure 22 and Table XVIII show that the temperature of the magnet is much lower than that of the winding and stator core. When the lubricating oil is still, the temperature of the magnet is 22.21 °C, which is the highest than that of the lubricating oil at different flow rates. When the lubricating oil has a flow rate, the temperature of the magnet decreases significantly; it is 4.56 °C with oil flow rate 0.5 m/s, which is about 80% lower than that of still. This is because the lubricating oil in the airgap is direct contact with the magnet, which can take away more heat, and when the lubricating oil has a flow rate, the heat taken away is greatly increased. It has a similar trend for rotor core temperature as shown in Fig. 23.

FIG. 22.

Temperature contours of the magnet at different lubricating oil flow rates. (a) Lubricating oil flow rate at 0 m/s. (b) Lubricating oil flow rate at 0.1 m/s. (c) Lubricating oil flow rate at 0.3 m/s. (d) Lubricating oil flow rate at 0.5 m/s.

FIG. 22.

Temperature contours of the magnet at different lubricating oil flow rates. (a) Lubricating oil flow rate at 0 m/s. (b) Lubricating oil flow rate at 0.1 m/s. (c) Lubricating oil flow rate at 0.3 m/s. (d) Lubricating oil flow rate at 0.5 m/s.

Close modal
FIG. 23.

Temperature contours of the rotor core at different lubricating oil flow rates. (a) Lubricating oil flow rate at 0 m/s. (b) Lubricating oil flow rate at 0.1 m/s. (c) Lubricating oil flow rate at 0.3 m/s. (d) ;-Lubricating oil flow rate at 0.5 m/s.

FIG. 23.

Temperature contours of the rotor core at different lubricating oil flow rates. (a) Lubricating oil flow rate at 0 m/s. (b) Lubricating oil flow rate at 0.1 m/s. (c) Lubricating oil flow rate at 0.3 m/s. (d) ;-Lubricating oil flow rate at 0.5 m/s.

Close modal

Figure 24 and Table XVIII indicate that the stator core has the highest temperature except for the windings. When the lubricating oil flow rate is 0 m/s, its maximum temperature is 43.07 °C, which is about 24 °C lower than that of winding. That is because the core loss of the motor is smaller than the copper loss when the motor is running under the rated state. In addition, the stator core is next to the seawater through the housing that has high thermal conductivity and rapid heat dissipation ability.

FIG. 24.

Temperature contours of the stator core at different lubricating oil flow rates. (a) Lubricating oil flow rate at 0 m/s. (b) Lubricating oil flow rate at 0.1 m/s. (c) Lubricating oil flow rate at 0.3 m/s. (d) Lubricating oil flow rate at 0.5 m/s.

FIG. 24.

Temperature contours of the stator core at different lubricating oil flow rates. (a) Lubricating oil flow rate at 0 m/s. (b) Lubricating oil flow rate at 0.1 m/s. (c) Lubricating oil flow rate at 0.3 m/s. (d) Lubricating oil flow rate at 0.5 m/s.

Close modal

When the lubricating oil inside the motor has a flow rate, the motor temperature decreases significantly compared with that with oil still. However, when the flow rate of lubricating oil continues to increase, the temperature of each part does not decrease significantly. Therefore, when designing the motor cooling system, it is only necessary to give a certain flow rate to the lubricating oil, and the requirements for its flow speed are not high.

In addition to the temperature field analysis of the underwater 2500 m, the electromagnetic field analysis of the optimized motor is also calculated.

Figures 25(a) and 25(b), respectively, show the magnetic flux density and the magnetic flux density vector after optimization under load. As shown in Fig. 25, when the motor is running under load, the overall distribution of magnetic flux of the motor is relatively uniform, and there is no magnetic saturation, indicating that there is less magnetic flux leakage. Therefore, the motor after optimization has reasonable electromagnetic distribution.

FIG. 25.

Loading magnetic flux density and magnetic density vector distribution diagram after optimization. (a) Magnetic flux density distribution. (b) Magnetic flux density vector distribution.

FIG. 25.

Loading magnetic flux density and magnetic density vector distribution diagram after optimization. (a) Magnetic flux density distribution. (b) Magnetic flux density vector distribution.

Close modal

From the loading torque curve shown in Fig. 26(a) and the loading three-phase current curve in Fig. 26(b), it can be seen that the stable output power is about 1.375 K N/m and the motor is gradually stable after 50 ms, which meets the design requirements.

FIG. 26.

Loading torque and current after optimization. (a) The curve of output torque under loading. (b) Three-phase current curve under loading.

FIG. 26.

Loading torque and current after optimization. (a) The curve of output torque under loading. (b) Three-phase current curve under loading.

Close modal

Through the temperature field calculation, it can be concluded that when the motor is running under rated conditions, the operating temperature distribution of the optimized motor meets the design requirements with different oil flow rates, and the magnet will not be demagnetized at high temperature. Meanwhile, the electromagnetic analysis indicates that the electromagnetic field distribution and output power are all rational and meet the design requirements. Therefore, the PMSM designed in this paper is reasonable.

Figure 27 presents the prototype of the UDPH, which has undergone validation through both land-based simulation tests and sea trials.49 As illustrated in Fig. 28, the driving system of this pile hammer employs a traditional three-phase asynchronous motor. Unlike PMSM, the rotor of this motor is composed of excitation windings, and its interior is filled with pure water instead of insulating oil. The motor has an external diameter of ∼370 mm, and when filled with water, it weighs around 1.2 tons. Furthermore, the hydraulic system of this UDPH is powered by four three-phase asynchronous motors. In all tests, including the land-based experiments and sea trials, all systems demonstrated stable performance.

FIG. 27.

The 450 kJ UDPH prototype of underwater 2500 m depth.

FIG. 27.

The 450 kJ UDPH prototype of underwater 2500 m depth.

Close modal
FIG. 28.

The original induction motor prototype for UDPH.

FIG. 28.

The original induction motor prototype for UDPH.

Close modal

In order to test the output performance of the induction motor used by the UDPH, a motor is employed as loading as shown in Fig. 29. The results are basically consistent with those of the sea trial experiments.

FIG. 29.

The original induction motor output performance test.

FIG. 29.

The original induction motor output performance test.

Close modal

To demonstrate the benefits of the PMSM for UDPH, a performance comparison between the traditional three-phase asynchronous motor used in the prototype and the PMSM introduced in this study is conducted. The results are detailed in Table XIX.

TABLE XIX.

Comparison between the conventional motor and new PMSM.

FactorsOriginal induction motorPMSM
Rated power 132 kW 132 kW 
Weight 1.3 t 0.72 t 
Efficiency 85% 90% 
Power factor 0.85 0.98 
Winding temperature 75.6 °C 61.8 °C 
Speed control Asynchronous Synchronous 
FactorsOriginal induction motorPMSM
Rated power 132 kW 132 kW 
Weight 1.3 t 0.72 t 
Efficiency 85% 90% 
Power factor 0.85 0.98 
Winding temperature 75.6 °C 61.8 °C 
Speed control Asynchronous Synchronous 

According to Table XIX, the PMSM is noticeably lighter than the asynchronous motor, with a weight reduction of 0.58 t per unit and a total decrease of 2.32 t across four units. When using the PMSM, this reduced weight results in less gravitational potential energy during the hammering process, which leads to energy conservation. When operating under stable conditions, the PMSM registers a temperature of ∼61.8 °C, in contrast to the 75.6 °C of the three-phase asynchronous motor. A significant factor to consider is that both the stator and rotor of the three-phase asynchronous induction motor are equipped with copper windings, leading to more copper loss. The PMSM, on the other hand, employs permanent magnets instead of rotor windings, effectively eliminating rotor copper loss. This design not only ensures a lower temperature rise but also enhances motor insulation. In addition, the PMSM’s capacity to adjust speed in real-time is superior to that of the three-phase asynchronous motor in terms of accuracy.

The comparison of outer diameter and axial length is given in Table XX. It is obvious that PMSM has a smaller volume compared to the traditional induction motor, which is about 90 mm smaller in outer diameter and about 250 mm shorter in axial length. This is of great help in reducing the weight of UDPH and improving energy efficiency.

TABLE XX.

Comparison of dimensions between two motors.

FactorsOriginal induction motor (mm)PMSM (mm)
Outer diameter 370 280 
Axial length 2800 2450 
FactorsOriginal induction motor (mm)PMSM (mm)
Outer diameter 370 280 
Axial length 2800 2450 

In this paper, based on the overall design parameters of the ultra-deepwater pile hammer, the underwater 2500 m depth, and the 450 kJ impact energy requirements, the calculation of rated parameters and the electromagnetic design of the PMSM are completed. Then, the finite element analysis of PMSM is carried out by Maxwell. The no-load and load analysis results show that the back EMF, magnetic density, and output torque of the motor are all reasonable. Finally, the various losses of the motor are analyzed, the temperature field simulation model is established, and the temperature field analysis under different lubricating oil flow rates is completed. The results show that the temperature distribution of the motor is reasonable. To further verify the superiority of the PMSM, we compared it with the traditional three-phase asynchronous motor based on the sea trial results of the 450 kJ and 2500 m UDPH prototype. From the analysis, it is evident that the PMSM demonstrates significant advantages in terms of weight, volume, and temperature rise. Through the simulation analysis of the electromagnetic field and temperature field, combined with the comparison of sea trial results with the traditional three-phase asynchronous motor, we further confirmed the rationality and efficiency of the PMSM.

However, the research and application of PMSM for UDPH are still in the early stages, and there are several issues that require further investigation. The current cooling structure is relatively simple, leading to a high temperature rise in the windings, which, in turn, demands higher insulation requirements. In addition, in terms of optimization, multi-objective optimization design has not been pursued, and only manufacturing errors have been considered. In the future, the motor’s cooling structure should be thoughtfully designed to reduce the temperature rise, thereby enhancing the motor’s efficiency. Furthermore, continued research in optimization is essential. Emphasis should be placed on multi-objective optimization design and the application of global optimization algorithms to identify the best solutions.

The authors would like to thank the editor and the reviewers for their helpful comments and suggestions. This research was funded by the National Key R&D Program of China (Grant Nos. 2021YFB3401400 and 2021YFB3401402), the Fundamental Research Funds for the Central Universities—the Opening Fund of National Engineering Laboratory of Offshore Geophysical and Exploration Equipment, China University of Petroleum, Qingdao Grant No. 266580, China (Grant No. 20CX02303A), and the Ministry of Industry and Information Technology of the P R China (Grant No. CH02N20) under Grant No. (2020)313.

The authors declare no conflict of interest.

Liping Tan: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Investigation (equal); Writing – original draft (equal); Writing – review & editing (equal). Yucong Wang: Data curation (equal); Formal analysis (equal); Software (equal). Teng Wang: Formal analysis (equal); Validation (equal). Changjiang Li: Methodology (equal); Validation (equal). Fenghua Miao: Investigation (equal); Software (equal). Wensheng Xiao: Data curation (equal); Formal analysis (equal); Investigation (equal). Junguo Cui: Data curation (equal); Funding acquisition (equal); Project administration (equal). Hongyan Wang: Data curation (equal); Investigation (equal).

Data sharing is not applicable to this article as no new data were created or analyzed in this study.

1
E.
Zhu
, “
Overview of the development of hydraulic pile hammers at home and abroad
,”
Fluid Transm. Control
79
(
6
),
5
11
(
2016
).
2
A.
Jain
,
T.
Gupta
, and
R. K.
Sharma
, “
Measurement and repair techniques of corroded underwater piles: An overview
,”
Int. J. Eng. Res. Appl
6
(
7
),
19
27
(
2016
).
3
H. R.
Stewart
, in
Sable Compression Facilities: Technical Overview, Offshore Technology Conference, 2007
(
OnePetro
,
2007
).
4
F.
Gozeling
and
E. G.
Van Der Velde
, “
Pile driving prediction contest
,” in
Application of Stress-Wave Theory to Piles: Test Results
(
CRC Press
,
2021
), pp.
103
140
.
5
P.
Brunning
and
N.
Zhen
, in
Some Insights Regarding the Installation of Driven Piles for Deep Water Developments, The Twenty-Fifth International Ocean and Polar Engineering Conference
(
OnePetro
,
2015
).
6
D.
Arentsen
, “
Pile driving demonstration
,”
Application of Stress-Wave Theory to Piles: Test Results
(
CRC Press
,
2021
), pp.
141
154
.
7
T.
Jokinen
,
V.
Hrabovcova
, and
J.
Pyrhonen
,
Design of Rotating Electrical Machines
(
John Wiley & Sons
,
2013
).
8
H.
Abu-Rub
,
A.
Iqbal
, and
J.
Guzinski
,
High Performance Control of AC Drives with Matlab/Simulink
(
John Wiley & Sons
,
2021
).
9
A.
Hekmati
and
I.
Sadeghi
, “
A comprehensive review on high speed permanent magnet motors and their modern applications
,”
Electromech. Energy Convers. Syst.
1
(
2
),
61
67
(
2021
).
10
S.
He
,
S.
Shi
,
Y.
Zhang
, and
W.
Chen
, “
Design and experimental research on a deep-sea resonant linear ultrasonic motor
,”
IEEE Access
6
,
57249
57256
(
2018
).
11
MENCK Subsea Pile Driving Hammers Sizes 150–300kj, http://217.199.187.193/menck.com/index.php/pilesize/150-300kj (2023/8/31).
12
J. M.
Audibert
and
J.
Won
, in
Hybrid Subsea Foundation for Deepwater Developments, ISOPE International Ocean and Polar Engineering Conference, 2018
(
ISOPE
,
2018
), pp. ISOPE-I-18-148.
13
C.
Pop
,
V.
Zania
, and
B.
Trimoreau
, “
Numerical modelling of offshore pile driving
,” in
Proceedings of XVI European Conference on Soil Mechanics and Geotechnical Engineering
, Edinburgh, UK, 13-17 September 2015 (International Society for Soil Mechanics and Geotechnical Engineering, 2015), pp.
1351
1356
.
14
Q.
Jiang
,
X.
Wang
,
M.
Yu
,
M.
Tang
,
H.
Zhang
, and
J.
Wu
, “
Theoretical study of vibro-acoustics of fluid-pile-soil coupled system and experimental research of noise reduction of small-scale pile driving
,”
Ocean Eng.
252
,
110997
(
2022
).
15
P.
Ramesh
and
N. C.
Lenin
, “
High power density electrical machines for electric vehicles—Comprehensive review based on material technology
,”
IEEE Trans. Magn.
55
(
11
),
0900121
(
2019
).
16
M.
Coutinho
,
D.
Bento
,
A.
Souza
,
R.
Cruz
,
F.
Afonso
,
F.
Lau
,
A.
Suleman
,
F. R.
Barbosa
,
R.
Gandolfi
,
W.
Affonso
, Jr.
et al, “
A review on the recent developments in thermal management systems for hybrid-electric aircraft
,”
Appl. Therm. Eng.
227
,
120427
(
2023
).
17
W.
Cui
,
B.
Hu
, and
T.
Pan
, “
A novel PMSM with optimized asymmetric eccentric air-gap for high-power electric motorcycle
,”
Electr. Power Compon. Syst.
47
(
19–20
),
1854
1863
(
2019
).
18
K.
Matsuoka
, “
Development trend of the permanent magnet synchronous motor for railway traction
,”
IEEE J. Trans. Electr. Electron. Eng.
2
(
2
),
154
161
(
2007
).
19
S. J.
Gambhire
,
D. R.
Kishore
,
P. S.
Londhe
, and
S. N.
Pawar
, “
Review of sliding mode based control techniques for control system applications
,”
Int. J. Dyn. Control
9
,
363
378
(
2021
).
20
W.
Xiao
,
L.
Tan
,
J.
Cui
,
H.
Wang
,
C.
Li
,
H.
Qin
, and
J.
Zhang
, “
A novel thermal analysis method based on a multi-physics two-way coupled method and its application to submersible permanent magnet synchronous motors
,”
Electronics
12
(
5
),
1155
(
2023
).
21
G.
Tian
,
L.
Jiang
,
L.
Yuan
,
Z.
Xu
,
Y.
Zhang
, and
C.
Li
, “
The design and implementation of control and driving system with double position feedback for submarine steering
,”
J. Phys.: Conf. Ser.
2187
,
012019
(
2022
).
22
P.
Mou
,
Z.
Xu
,
J.
Bao
,
J.
Zhu
, and
Z.
Zheng
, “
Design of controllable power system for underwater transmission mechanism
,” in
International Conference of Fluid Power and Mechatronic Control Engineering (ICFPMCE 2022)
(
Atlantis Press
,
2022
), pp.
235
240
.
23
P. T.
Luu
,
J.
Lee
,
J.
Lee
, and
J.
Park
, “
Electromagnetic and thermal analysis of permanent-magnet synchronous motors for cooperative robot applications
,”
IEEE Trans. Magn.
56
(
3
),
7512804
(
2020
).
24
A.
Sheela
,
M.
Suresh
,
V. G.
Shankar
,
H.
Panchal
,
V.
Priya
,
M.
Atshaya
,
K. K.
Sadasivuni
, and
S.
Dharaskar
, “
FEA based analysis and design of PMSM for electric vehicle applications using magnet software
,”
Int. J. Ambient Energy
43
(
1
),
2742
2747
(
2022
).
25
L.
Wang
and
D.
Chen
, “
Anvil configuration optimization of an underwater pile driving hammer
,” in
Advances in Energy and Environment Research: Proceedings of the International Conference on Advances in Energy and Environment Research (ICAEER2016), Guangzhou City, China, August 12–14, 2016, 2017
(
CRC Press
,
2017
), p.
375
.
26
P.
Pillay
and
R.
Krishnan
, “
Modeling, simulation, and analysis of permanent-magnet motor drives. I. The permanent-magnet synchronous motor drive
,”
IEEE Trans. Ind. Appl.
25
(
2
),
265
273
(
1989
).
27
J.
Xia
,
T.
Dong
,
C.
Wang
, and
J.
Zhao
, “
Low speed high torque PMSM design based on unequal teeth structure
,” in
2008 International Conference on Electrical Machines and Systems
(
IEEE
,
2008
), pp.
3274
3277
.
28
M.
Mikhak-Beyranvand
,
J.
Faiz
, and
J.
Khoshtarash
, “
Impacts of number of Poles and slots on armature reaction and performance of ironless permanent magnet motors
,”
Electr. Power Compon. Syst.
48
(
19–20
),
1979
1991
(
2021
).
29
H.
Guo
,
X.
He
,
J.
Xu
,
W.
Tian
,
G.
Sun
,
L.
Ju
, and
D.
Li
, “
Design of an Aviation dual-three-phase high-power high-speed permanent magnet assisted synchronous reluctance starter-generator with antishort-circuit ability
,”
IEEE Trans. Power Electron.
37
(
10
),
12619
12635
(
2022
).
30
J.
Jin
,
Electromagnetic Analysis and Design in Magnetic Resonance Imaging
(
Routledge
,
2018
).
31
E.
Rahimi Namaghi
and
H. R.
Izadfar
, “
Introducing Analytical formulas for calculating back EMF of the contra-rotating double-rotor axial flux permanent magnet synchronous machine
,”
Electr. Power Compon. Syst.
51
(
7
),
684
693
(
2023
).
32
G. T.
de Paula
,
J. R. B. D.
Monteiro
,
B. P.
de Alvarenga
,
T. E.
de Almeida
,
W. C.
Pereira
, and
M. P.
de Santana
, “
On-load back emf of PMSM using Maxwell stress tensor
,”
IEEE Trans. Magn.
54
(
7
),
8107511
(
2018
).
33
M. S.
Kumar
and
S. T.
Revankar
, “
Development scheme and key technology of an electric vehicle: An overview
,”
Renewable Sustainable Energy Rev.
70
,
1266
1285
(
2017
).
34
T.
Mohanarajah
,
M.
Nagrial
,
J.
Rizk
, and
A.
Hellany
, “
A novel method to optimize permanent magnet assisted synchronous reluctance machines
,”
Electr. Power Compon. Syst.
48
(
9–10
),
933
943
(
2020
).
35
Z.
Zhu
,
S.
Xue
,
W.
Chu
,
J.
Feng
,
S.
Guo
,
Z.
Chen
, and
J.
Peng
, “
Evaluation of iron loss models in electrical machines
,”
IEEE Trans. Ind. Appl.
55
(
2
),
1461
1472
(
2019
).
36
L.
Cuiping
,
C.
Pujia
,
C.
Shukang
, and
C.
Feng
, “
Research on motor iron losses and temperature field calculation for mini electric vehicle
,” in
2014 17th International Conference on Electrical Machines and Systems (ICEMS)
(
IEEE
,
2014
), pp.
2380
2383
.
37
D.
Ouamara
and
F.
Dubas
, “
Permanent-magnet eddy-current losses: A global revision of calculation and analysis
,”
Math. Comput. Appl.
24
(
3
),
67
(
2019
).
38
A.
Hemeida
,
P.
Sergeant
, and
H.
Vansompel
, “
Comparison of methods for permanent magnet eddy-current loss computations with and without reaction field considerations in axial flux PMSM
,”
IEEE Trans. Magn.
51
(
9
),
8107511
(
2015
).
39
B.
Bilgin
,
J.
Liang
,
M. V.
Terzic
,
J.
Dong
,
R.
Rodriguez
,
E.
Trickett
, and
A.
Emadi
, “
Modeling and analysis of electric motors: State-of-the-art review
,”
IEEE Trans. Transp. Electrif.
5
(
3
),
602
617
(
2019
).
40
H.
Ahn
,
H.
Park
,
C.
Kim
, and
H.
Lee
, “
A review of state-of-the-art techniques for PMSM parameter identification
,”
J. Electr. Eng. Technol.
15
(
3
),
1177
1187
(
2020
).
41
A.
Tikadar
,
N.
Kumar
,
Y.
Joshi
, and
S.
Kumar
, “
Coupled electro-thermal analysis of permanent magnet synchronous motor for electric vehicles
,” in
2020 19th IEEE Intersociety Conference on Thermal and Thermomechanical Phenomena in Electronic Systems (ITherm), 2020
(
IEEE
,
2020
), pp.
249
256
.
42
R.
Wrobel
,
P. H.
Mellor
,
M.
Popescu
, and
D. A.
Staton
, “
Power loss analysis in thermal design of permanent-magnet machines—a review
,”
IEEE Trans. Ind. Appl.
52
(
2
),
1359
1368
(
2015
).
43
S.
Sirimanna
,
T.
Balachandran
, and
K.
Haran
, “
A review on magnet loss analysis, validation, design considerations, and reduction strategies in permanent magnet synchronous motors
,”
Energies
15
(
17
),
6116
(
2022
).
44
S.
Khalesidoost
,
J.
Faiz
, and
E.
Mazaheri Tehrani
, “
An overview of thermal modelling techniques for permanent magnet machines
,”
IET Sci., Meas. Technol.
16
(
4
),
219
241
(
2022
).
45
P.
Roy
,
M.
Towhidi
,
F.
Ahmed
,
A. J.
Bourgault
,
S.
Mukandan
,
A.
Balamurali
, and
N. C.
Kar
, “
A comprehensive review of thermal design and analysis of traction motors
,” in
2019 IEEE 28th International Symposium on Industrial Electronics (ISIE), 2019
(
IEEE
,
2019
), pp.
203
208
.
46
A.
Tikadar
,
D.
Johnston
,
N.
Kumar
,
Y.
Joshi
, and
S.
Kumar
, “
Comparison of electro-thermal performance of advanced cooling techniques for electric vehicle motors
,”
Appl. Therm. Eng.
183
,
116182
(
2021
).
47
Q.
Zhang
,
Z.
Li
,
X.
Dong
,
Y.
Liu
, and
R.
Yu
, “
Comparative analysis of temperature field model and design of cooling system structure parameters of submersible motor
,”
J. Therm. Sci. Eng. Appl.
13
(
5
),
051002
(
2021
).
48
Z.
Li
,
Q.
Zhang
,
R.
Yu
,
L.
Zhang
, and
C.
Xu
, “
Analysis of influence factors of underwater motor temperature field based on magneto-thermal-flux coupling
,” in
2018 IEEE International Conference on Mechatronics and Automation (ICMA), 2018
(
IEEE
,
2018
), pp.
773
778
.
49
L.
Dan
, “
China’s first ultra-deepwater piling hammer sea trial was successful
,” http://tj.people.com.cn/n2/2023/0705/c375366-40482313.html.