Multi-Parameter Optimization of Heat Dissipation Structure of Double Disk Magnetic Coupler Based on Orthogonal Experimental Design

: The existing heat dissipation research on double disk magnetic couplers ignores the coupling inﬂuence of electromagnetic temperature–stress and other multiphysics ﬁelds, and the error between the calculation and analysis results and the measured values is large. Therefore, a multi-parameter optimization method for heat dissipation structures of double disk magnetic couplers based on orthogonal experimental design is proposed. Based on the double disk magnetic coupler model, a three-dimensional ﬁnite element model based on ﬂuid–solid–heat coupling is established, with the axial air gap length, input motor speed, the thickness of the permanent magnet in the magnetizing direction, the thickness of the copper plate, the number of ﬁns of the heat dissipation plate and the length of the ﬁns of the heat dissipation plate as design variables. Six-factor and three-level simulation experiments are designed with the minimum temperature of the heat dissipation plate as the objective function, and additionally, orthogonal experiments were designed according to the actual working conditions by selecting the optimal combination of parameters and modifying the model to perform physical tests. The results show that the variables that have the most signiﬁcant impact on heat dissipation performance from high to low are as follows: axial air gap length, input motor speed, the length of the ﬁns of the heat dissipation plate, the thickness of the permanent magnet in the magnetizing direction, the number of ﬁns of the heat dissipation plate and the thickness of the copper plate. The increase in axial air gap length can effectively reduce the temperature rise, and the maximum decrease can reach 9.76%. Under the same conditions, the input motor speeds are set to 300 r/min, 400 r/min, 500 r/min, 600 r/min and 700 r/min, respectively, and the simulation results are in good agreement with the physical test results, with a maximum error of 4.8%. The error between the simulation result and the physical test result is only 1.9% under the optimal combination of parameters obtained by the orthogonal experiment, which veriﬁes the correctness of the optimization model. In conclusion, the study is of reference signiﬁcance for the parameter optimization of the heat dissipation structure of the double disk magnetic coupler.


Introduction
Double disk magnetic couplers, a revolutionary energy-saving transmission device, boast the advantages of high efficiency speed regulation, strong adaptability and low running costs. They have become one of the key research objects of magnetic transmission systems in recent years [1,2]. During the operation of the coupler, the magnetic induction eddy current loss will generate a significant amount of heat, leading to an increase in the temperature of all parts of the transmission system. Too high a temperature will affect the working characteristics of permanent magnet materials [3,4], which in turn affects the reliability of the transmission system. Electromagnetic temperature-stress coupling energy is a rather complicated problem, which has become one of the hot topics in the field of magnetic coupling [5].
As far as international studies are concerned, R Bujakiewicz-korońska and J Koroński carried out a study on the principle of minimum entropy production for diffusion and heat transfer in an open system, deriving general formulas and typical values of the entropy production for diffusion and heat transfer for any substance that involves specific boundary conditions [6]; Andrew T. Morrision et al. explored the temperature field of a heat sink in a natural convection state by using finite element software, and calculated the thermodynamic parameters of heat sink fins in a thermal steady state using an iterative method [7]; Fabrizio Madgnetti et al. performed finite element analysis by means of a dynamic coupling method of temperature and fluid, which provided guidance for the temperature rise characteristics of a disk magnetic coupler [8]. A study was also carried out by Avram Bar-cohen in which the analytical formula of the Nusselt number was applied to design the size of a heat sink and derive the optimal size based on the "least material" theory by taking the total heat dissipation and space heat transfer coefficient of the device as the objective function [9]. As far as domestic studies are concerned, Yang Chaojun simulated and analyzed the mechanical characteristics and speed-regulating characteristics of a speedregulating magnetic coupler by using the finite element method, contributing to a better realization of the speed regulation function for the coupler [10]; Wang Lei et al. further put forward a method to calculate the heat dissipation coefficient based on the fluid-solid coupling velocity field and then analyzed the coupler temperature field, demonstrating the applicability of the coupler temperature field distribution to the heat dissipation surfaces of different couplers [11]; GE Yanjun et al. designed a magnetic coupler with the outer rotor being a permanent magnet rotor and the inner rotor being a copper sleeve rotor, and verified the heat dissipation effect of the heat dissipation structure on the magnetic coupler by simulation analysis [12]. To sum up, a plurality of studies on magnetic couplers have been carried out at home and abroad, but none of these studies have employed the fluid-solid-heat coupling analysis method to optimize the multi-parameters for the heat dissipation problem of the double disk magnetic coupler. Its structural parameters, the size of the copper plate, the air gap of the permanent magnet disk, the input speed of the motor and the load are the key factors that affect the heat dissipation performance.
The existing heat dissipation research on double disk magnetic couplers ignores the coupling influence of electromagnetic temperature-stress and other multiphysics fields, and the error between the simulation calculation and analysis results and the measured values is large. Therefore, this paper builds a finite element simulation platform for permanent magnetic field-temperature field coupling based on the thermodynamic model for the heat dissipation structure parameters of the double disk magnetic coupler. Taking the minimum temperature of the heat dissipation plate as the objective function and the mass change and average output speed as the constraint function, an orthogonal experimental method was used to optimize the heat dissipation structure with multi-parameters to obtain the optimal structure layout and parameters. Under the optimal parameter combination conditions, physical experiments are carried out for verification.

Structure
The double disk magnetic coupler is mainly composed of an input shaft, heat dissipation plate, conductor yoke, copper plate, permanent magnet plate, permanent magnet yoke, gap-adjusting device and output shaft. Considering the symmetry of the left and right sides of the double disk magnetic coupler, only the left side structural model is analyzed in this paper. Specifically, the input shaft is connected to the input motor, and a large amount of heat is generated by the magnetic induction eddy current loss, which is dissipated by the annular fins (fixed on the left side of the conductor yoke) on the radiating plate. The conductor yoke can not only fix the copper plate, but also guide the magnetic lines generated by the permanent magnets, so that more magnetic lines enter the copper plate to improve Energies 2022, 15, 8801 3 of 18 the transmission efficiency, increase the local magnetic induction intensity, and prevent or reduce the magnetic flux leakage between the permanent magnets [13]. The copper plate is fixed on the right side of the conductor yoke by hexagon bolts. The distance between the copper plate and the permanent magnet disk is equal, with air as the gas medium. The permanent magnet disk is composed of a permanent magnet, an aluminum disk and a flange, where the aluminum disk is fixed on the permanent magnet yoke. The permanent magnet yoke is fixed on the right side of the permanent magnet disk, which plays a role in controlling the magnetic flux direction and increasing the magnetic flux density in the magnetic circuit. The two permanent magnet disks are connected by a gap-adjusting device and are symmetrically distributed. The conductor yoke is connected to the output shaft, as shown in Figure 1.
large amount of heat is generated by the magnetic induction eddy current loss, which is dissipated by the annular fins (fixed on the left side of the conductor yoke) on the radiating plate. The conductor yoke can not only fix the copper plate, but also guide the magnetic lines generated by the permanent magnets, so that more magnetic lines enter the copper plate to improve the transmission efficiency, increase the local magnetic induction inten sity, and prevent or reduce the magnetic flux leakage between the permanent magnets [13]. The copper plate is fixed on the right side of the conductor yoke by hexagon bolts The distance between the copper plate and the permanent magnet disk is equal, with air as the gas medium. The permanent magnet disk is composed of a permanent magnet, an aluminum disk and a flange, where the aluminum disk is fixed on the permanent magne yoke. The permanent magnet yoke is fixed on the right side of the permanent magnet disk which plays a role in controlling the magnetic flux direction and increasing the magnetic flux density in the magnetic circuit. The two permanent magnet disks are connected by a gap-adjusting device and are symmetrically distributed. The conductor yoke is connected to the output shaft, as shown in Figure 1.

Working Principle
When the input motor is energized, the input shaft drives the copper plate on the right side of the conductor yoke to rotate, and the copper plate and the permanent magne disk move, resulting in slip. According to Faraday's electromagnetic induction principle there is an eddy current between the copper plate and the permanent magnet disk, and the induced magnetic field generated by the eddy current interacts with the magnetic field generated by the permanent magnet to generate torque, which generates heat on the sur face of the copper plate. The generated torque drives the permanent magnet disk to rotate in the same direction as the copper plate and drives the output shaft to rotate, thus trans mitting power to the load [14].

Gas Heat Transfer Control Equation
Fluid-solid-heat coupling calculations are mainly used to solve the problem that gas flow, heat transfer and structural deformation need to be considered at the same time. The key is the heat transfer at the fluid-solid boundary. The heat transfer process can be ex pressed as follows:

Working Principle
When the input motor is energized, the input shaft drives the copper plate on the right side of the conductor yoke to rotate, and the copper plate and the permanent magnet disk move, resulting in slip. According to Faraday's electromagnetic induction principle, there is an eddy current between the copper plate and the permanent magnet disk, and the induced magnetic field generated by the eddy current interacts with the magnetic field generated by the permanent magnet to generate torque, which generates heat on the surface of the copper plate. The generated torque drives the permanent magnet disk to rotate in the same direction as the copper plate and drives the output shaft to rotate, thus transmitting power to the load [14].

Gas Heat Transfer Control Equation
Fluid-solid-heat coupling calculations are mainly used to solve the problem that gas flow, heat transfer and structural deformation need to be considered at the same time.
The key is the heat transfer at the fluid-solid boundary. The heat transfer process can be expressed as follows: where λ is the coefficient of thermal conductivity of the air medium (0.023 W/m·k); T is the temperature detected on the surface of the heat dissipation plate; ∂T ∂n indicates the temperature gradient of the heat dissipation plate in the normal direction of the fluid-solid interface; α is the convective heat transfer coefficient of air; T f is the temperature of the air. The heat conduction equation between solids based on the three-dimensional heat conduction differential equation is as follows [15]: where Φ is the heating power of the input motor. If air is regarded as an incompressible gas, then in the calculation of fluid-solid-heat coupling, the mass conservation equation, momentum conservation equation and energy conservation equation should be taken into account. Its three-dimensional flow and heat transfer control equation is as follows [16]: where v is the velocity vector; ρ, f, P, µ, c P , T 1 and k represent the density, volume force, pressure, dynamic viscosity, specific heat capacity, temperature and thermal conductivity of air, respectively.

Boundary Parameter Setting
In this paper, the heat dissipation structure of a double disk magnetic coupler is taken as the research object. Figure 2 shows a simplified model of the heat dissipation structure of a double disk magnetic coupler, including a heat source, copper plate, conductor yoke and heat dissipation plate. The end face of the copper plate close to the permanent magnet plate is simplified as a heat source-generating surface, and the heat flow boundary is indicated by an arrow. The conductor yoke is connected to the fins, so that the heat generated on the copper plate can be uniformly and sequentially transmitted to the yoke plate, the heat dissipation plate and the fins of the heat dissipation plate. In this paper, the gas heat transfer control equation is calculated by the steady-state solution [17]. The heat generated by the eddy current loss is transferred from the surface of the copper plate to the surface of the fins of the heat dissipation plate, and finally, the heat is transferred away by the heat exchange between the fin surface and the environment, so as to avoid the accumulation of a large amount of heat on the surface of the copper plate. In addition, all the joints are filled with thermal grease to reduce the contact thermal resistance between components.      Table 2 displays the material properties of each component.

Fluid-Solid-Heat Coupling Finite Element Simulation
In most of the current studies on magnetic couplers, the coupling effects of multiple physical fields, such as temperature field and thermal field, are neglected, resulting in large errors between the calculated analysis results and the measured values. Therefore, in view of the influence of eddy current heat on the temperature of the heat dissipation plate in the three-dimensional magnetic field simulation of a double disk magnetic coupler, a three-dimensional finite element analysis method based on fluid-solid-heat coupling is proposed to calculate the temperature field at the heat sink and solve the temperature of the heat dissipation plate.

Three-Dimensional Model of Heat Dissipation Structure
The axial air gap length of the double disk magnetic coupler was set to 5 mm, the input motor speed to 400 r/min, the number of heat dissipation plate fins to 6, and the length of the heat dissipation plate fins to 8 mm. Its three-dimensional model is shown in Figure 3.  Based on the three-dimensional model of the heat dissipation structure of the double disk magnetic coupler, a steady-state thermal simulation is carried out. The solid domain part and the fluid domain part [18] of the heat dissipation disk are established using a bidirectional coupling solution, in which the solid domain adopts the finite element method and the fluid domain adopts the finite volume method. The rotation domain is selected as the fluid domain, with the speed assigned to 400 r/min, and the grid is set as the MRF model. The meshing is shown in the following Figure 4. Based on the three-dimensional model of the heat dissipation struc disk magnetic coupler, a steady-state thermal simulation is carried out. part and the fluid domain part [18] of the heat dissipation disk are es bidirectional coupling solution, in which the solid domain adopts t method and the fluid domain adopts the finite volume method. The r selected as the fluid domain, with the speed assigned to 400 r/min, and the MRF model. The meshing is shown in the following Figure 4. When meshing, the global mesh size is set as 8 mm, the number o the number of meshes is 4929, and the simulation time is 20 s. After gr judged that the grid meets the calculation requirements.

Permanent Magnetic Field and Temperature Field Coupling Platform
Combining Ansoft Maxwell and Ansys Workbench software, a per field-temperature field coupling platform is built. Considering the sym and right sides of the double disk magnetic coupler, only the left side he is analyzed in order to simplify the calculation volume. To analyze the e rent heat on the performance of the copper plates and a permanent ma power generated by the eddy current is introduced into the temperatur sient thermal module in turn as the heat source load. The heat generated nal heat generation stage is applied to the volume of the copper plate, a rent thermal power on the copper plate is simulated in the magnetic fie ulation results of the permanent magnetic field are combined with the t and the temperature of the heat dissipation plate is calculated by simula tion diagram is shown in Figure 5. When meshing, the global mesh size is set as 8 mm, the number of nodes is 23,313, the number of meshes is 4929, and the simulation time is 20 s. After grid inspection, it is judged that the grid meets the calculation requirements.

Permanent Magnetic Field and Temperature Field Coupling Platform
Combining Ansoft Maxwell and Ansys Workbench software, a permanent magnetic field-temperature field coupling platform is built. Considering the symmetry of the left and right sides of the double disk magnetic coupler, only the left side heat transfer model is analyzed in order to simplify the calculation volume. To analyze the effect of eddy current heat on the performance of the copper plates and a permanent magnet, the thermal power generated by the eddy current is introduced into the temperature field of the transient thermal module in turn as the heat source load. The heat generated during the internal heat generation stage is applied to the volume of the copper plate, and the eddy current thermal power on the copper plate is simulated in the magnetic field. Then, the simulation results of the permanent magnetic field are combined with the temperature field, and the temperature of the heat dissipation plate is calculated by simulation. The simulation diagram is shown in Figure 5. When the ambient temperature is 22 °C, the air gap length is 5 mm, the input motor speed is 400 r/min, and the number of fins of the heat dissipation plate is 6, 8, 10 and 12 in turn, the minimum temperatures are 38.530 °C, 38.232 °C, 37.951 °C and 37.736 °C, respectively, and the maximum temperatures are 42.915 °C and 42.66 °C, respectively. It can be observed that when the number of fins of the heat dissipation plate gradually increases, the temperature range of the heat dissipation plate gradually decreases, and the minimum temperature of the model is concentrated in the center of the fins.

Structural Parameter Model
In view of the large number of structural parameters of the double disk magnetic coupler, the orthogonal experimental design of the structural parameters of the double disk magnetic coupler is carried out using the software Design Expect 8.0.6. Orthogonal experimental design, compared to other optimization methods, is more widely used by virtue of the fact that it requires fewer combinations of experiments and is generally used for the analysis of the nonlinear influence of factors [19]. In performing this experiment, an appropriate number of highly representative test points are selected from a large number of test points, and a normalized orthogonal table is used to arrange the experiment When the ambient temperature is 22 • C, the air gap length is 5 mm, the input motor speed is 400 r/min, and the number of fins of the heat dissipation plate is 6, 8, 10 and 12 in turn, the minimum temperatures are 38.530 • C, 38.232 • C, 37.951 • C and 37.736 • C, respectively, and the maximum temperatures are 42.915 • C and 42.66 • C, respectively. It can be observed that when the number of fins of the heat dissipation plate gradually increases, the temperature range of the heat dissipation plate gradually decreases, and the minimum temperature of the model is concentrated in the center of the fins.

Structural Parameter Model
In view of the large number of structural parameters of the double disk magnetic coupler, the orthogonal experimental design of the structural parameters of the double disk magnetic coupler is carried out using the software Design Expect 8.0.6. Orthogonal experimental design, compared to other optimization methods, is more widely used by virtue of the fact that it requires fewer combinations of experiments and is generally used for the analysis of the nonlinear influence of factors [19]. In performing this experiment, an appropriate number of highly representative test points are selected from a large number of test points, and a normalized orthogonal table is used to arrange the experiment scientifically. Furthermore, the sensitivity of each factor to the influence of the optimization target is analyzed to seek the optimal structure layout and parameters of the double disk magnetic coupler.
The orthogonal experimental design adopts the three-dimensional magnetic field simulation analysis method to simulate the response value that corresponds to the preselected experimental design point of the temperature of the heat dissipation plate of the double disk magnetic coupler, and the mathematical model is as follows: where T is the objective function, g j is the constraint function; the design variable is the analytic function of the objective function and the constraint function, the number of constraint functions is m, x = (x 1 , x 2 , . . . x n ) is the design variable, the number of design variables is n, and the upper and lower values of the ith design variable are x iL and x iU , respectively.

Improved Response Surface Design Principle
Using the improved Box-Behnken (response surface) design principle, the average output speed and the temperature of the heat dissipation plate of the double disk magnetic coupler in the stable operation stage are taken as the response values.
In this paper, a six-factor and three-level simulation experiment is conducted to select sample points in the space composed of six factors. The axial air gap length is taken as 5-15 mm, the input motor speed as 400-800 r/min, the thickness of the permanent magnet in the magnetizing direction as 12.7-16.7 mm, the thickness of the copper plate as 6-10 mm, the number of fins of the heat dissipation plate as 6-10, and the length of the fins of the heat dissipation plate as 8-12 mm. The values of the six factors have not yet been determined, so the axial air gap length, input motor speed, the thickness of the permanent magnet in the magnetizing direction, the thickness of the copper plate, the number of fins of the heat dissipation plate and the length of fins of the heat dissipation plate are selected as the design variables, which are recorded as A, B, C, D, E, and F, respectively. For the convenience of analysis, the dimensionless design variables are converted into dimensionless variables, whose three levels are −1, 0, 1. Dimensional transformation is defined as follows: The response surface test factors and levels are shown in Table 3 below.

Orthogonal Experimental Design and Simulation
The experimental software Design Expert 8.0.6 is used to simulate, calculate and obtain the experimental sample data, and the orthogonal experimental simulation results are shown in Table 4.

Analysis of Variance in Orthogonal Experiment
In order to test whether the established orthogonal experimental model can be further optimized, it is necessary to analyze the variance in the model and confirm the fitting accuracy. Table 5 displays the results of the variance analysis of the orthogonal experimental model of the temperature of the heat dissipation plate. After calculation, the negative correlation coefficient R 2 of the temperature of the heat dissipation plate is 0.9248, which meets the requirements of engineering accuracy.

Significant Impact
The F-value of the variance analysis of the orthogonal experimental results can be used as a parameter to respond to significant differences. The larger the F-value of the variance analysis, the higher the significant influence [20]. For the six factors, the axial air gap length is 219.53, the input motor speed is 14.29, the thickness of the permanent magnet in the magnetizing direction is 2.95, the thickness of the copper plate is 0.23, the number of fins of the heat dissipation plate is 2.27 and the length of the fins of the heat dissipation plate is 4.55. Therefore, the order of significant influence from high to low is A > B > F > C > E > D. It can be observed that the axial air gap length is the most significant factor affecting the heat dissipation performance, followed by the input motor speed, and its value should be reasonably selected.

Single Factor Analysis
The single factor that affects the temperature change in the heat dissipation plate is shown in Figure 6.
With the increase in single factor A (axial air gap length), the temperature of the heat dissipation plate decreases from about 47 • C to about 43 • C, and then shows a slight upward trend. With the increase in single factor B (input motor speed) and C (the thickness of the permanent magnet in the magnetizing direction), the temperature of the heat dissipation plate gradually rises, with a small increase. With the increase in single factor D (the thickness of the copper plate), the temperature change in the heat dissipation plate is small. With the increase in single factor E (the number of fins of the heat dissipation plate) and F (the length of the fins of the heat dissipation plate), the temperature of the heat dissipation plate gradually decreases, with a small decrease.
In the case of single factor analysis, the percentage of the six factors' rise and fall that affect the temperature change in the heat dissipation plate is shown in Table 6 below.
It can be observed from the above table that the axial air gap length decreases the most, that is, increasing the axial air gap length to a certain extent allows the temperature of the heat dissipation plate to reach its minimum level, which is closer to the objective function T min . With the increase in single factor A (axial air gap length), the temperature of the heat dissipation plate decreases from about 47 °C to about 43 °C, and then shows a slight upward trend. With the increase in single factor B (input motor speed) and C (the thickness of the permanent magnet in the magnetizing direction), the temperature of the heat dissipation plate gradually rises, with a small increase. With the increase in single factor D (the thickness of the copper plate), the temperature change in the heat dissipation plate is small. With the increase in single factor E (the number of fins of the heat dissipation plate) and F (the length of the fins of the heat dissipation plate), the temperature of the heat dissipation plate gradually decreases, with a small decrease.
In the case of single factor analysis, the percentage of the six factors' rise and fall that affect the temperature change in the heat dissipation plate is shown in Table 6 below. It can be observed from the above table that the axial air gap length decreases the most, that is, increasing the axial air gap length to a certain extent allows the temperature of the heat dissipation plate to reach its minimum level, which is closer to the objective function Tmin.

Two-Factor Analysis
To explain the orthogonal experimental model more intuitively, the three-dimensional graphic analysis of its two-factor influence is carried out by using the software of Design Expert 8.0.6. Figure 7 shows the influence of four groups of representative three-dimensional figures on the temperature of the heat dissipation plate when six factors are combined in pairs.

Two-Factor Analysis
To explain the orthogonal experimental model more intuitively, the three-dimensional graphic analysis of its two-factor influence is carried out by using the software of Design Expert 8.0.6. Figure 7 shows the influence of four groups of representative threedimensional figures on the temperature of the heat dissipation plate when six factors are combined in pairs.    From the analysis, it can be observed that the larger the value of A, the smaller the value of B and the lower the temperature of the heat dissipation plate, and the influence of AB on the objective function is inversely proportional. When A = 15.00 mm and B = 400 r/min, the temperature of the heat dissipation plate reaches its lowest value, about 42.7 • C. That is, with the increase in the axial air gap length and the decrease in the input motor speed, the temperature of the heat dissipation plate gradually decreases. Figure 7b shows the change in the temperature response surface of the heat dissipation plate at the response point with the influence factors A and F, when the influence factors are as follows: B = 600 r/min, C = 14.70 mm, D = 8.00 mm and E = 8.00. It can be observed from the analysis that the larger the value of A, the larger the value of F and the lower the temperature of the heat dissipation plate, and the influence of AF on the objective function is proportional. When A = 15.00 mm and F = 12.00 mm, the temperature of the heat dissipation plate reaches its lowest value, about 42.7 • C. That is, with the increase in the axial air gap length and the fin length of the heat dissipation plate, the temperature of the heat dissipation plate gradually decreases. Figure 7c shows the change in the temperature response surface of the heat dissipation plate at the response point with the influence factors B and E, when the influence factors are as follows: A = 10 mm, C = 14.70 mm, D = 8.00 mm and F = 10.00 mm. It can be observed that the smaller the value of B, the smaller the value of E and the lower the temperature of the heat dissipation plate, and the influence of BE on the objective function is directly proportional. When B = 400 r/min and E = 6, the temperature of the heat dissipation plate reaches its lowest value, about 43.0 • C. That is, with the decrease in the input motor speed, the number of fins of the heat dissipation plate decreases, and the temperature of the heat dissipation plate gradually decreases. Figure 7d shows the change in the temperature response surface of the heat dissipation plate at the response point with the influence factors E and F, when the influence factors are as follows: A = 10 mm, B = 600 r/min, C = 14.70 mm and D = 8.00 mm. It can be observed that the larger the value of E, the larger the value of F and the lower the temperature of the heat dissipation plate, and the influence of EF on the objective function is proportional. When E = 10 and F = 12 mm, the temperature of the heat dissipation plate reaches its lowest value, about 43 • C. That is, with the increase in the fin number and fin length of the heat dissipation plate, the temperature of the heat dissipation plate gradually decreases.
According to the two-factor influence change diagram of the temperature of the heat dissipation plate, the response and optimization direction of different influence factors can be intuitively expressed.

Simulation Optimal Conditions and Results
In practical problems, in order to ensure better transmission performance and control the cost of the double disk magnetic coupler, the mass of the optimized double disk magnetic coupler should not be higher than 5% of the original mass, and the initial value of the influencing factors is x 1 = x 2 = x 3 = x 4 = x 5 = x 6 = 0, namely x 1 = 10, x 2 = 600, x 3 = 14.7, x 4 = 8, x 5 = 8 and x 6 = 10. The structural parameter optimization model of the double disk magnetic coupler is as follows: Equation (8) is the multi-parameter optimization model established in this paper, with the minimum temperature T of the heat dissipation plate as the objective function and the mass change ∆m and average output v u speed as the constraint functions. m 0 means the initial mass of the double disk magnetic coupler, ∆m indicates the change in its mass during optimization, and v u is the lower limit of the average output speed. The optimal results of the orthogonal experiment are as follows: A = 14.87 mm, B = 501.42 r/min, C = 13.05 mm, D = 6.80 mm, E = 9.85 and F = 11.89 mm. Considering the difficulty of processing technology and the technical limitations of the actual installation of the double disk magnetic coupler experimental device, we calculated the above approximate values. The optimal solution is that the axial air gap length is 15.00 mm, the input motor speed is 500.00 r/min, the thickness of the permanent magnet in the magnetizing direction is 13.00 mm, the thickness of the copper plate is 7.00 mm, and the number of fins of the heat dissipation plate is 10, the length of the fins of the heat dissipation plate is 12.00 mm, and the minimum temperature of the heat dissipation plate is 41.8409 • C.

Physical Test
To verify the results of the theoretical analysis and finite element simulation mentioned above, the same parameters as those used in the simulation are selected for physical tests on a double disk magnetic coupler. The field test platform is set up as shown in Figure 8, and the infrared temperature sensor with a heat dissipation structure is shown in Figure 9.
The optimal solution is that the axial air gap length is 15.00 mm, the input motor speed is 500.00 r/min, the thickness of the permanent magnet in the magnetizing direction is 13.00 mm, the thickness of the copper plate is 7.00 mm, and the number of fins of the heat dissipation plate is 10, the length of the fins of the heat dissipation plate is 12.00 mm, and the minimum temperature of the heat dissipation plate is 41.8409 °C.

Physical Test
To verify the results of the theoretical analysis and finite element simulation mentioned above, the same parameters as those used in the simulation are selected for physical tests on a double disk magnetic coupler. The field test platform is set up as shown in Figure 8, and the infrared temperature sensor with a heat dissipation structure is shown in Figure 9.  The input motor is a 380 V three-phase low-speed motor, which is controlled by a frequency converter. The rated power of this test-bed is 55 kW, the maximum speed is 3000 r/min, and the recommended speed range is 0-2000 r/min. The measuring circuit is equipped with NCTES3000 torque and a speed measuring instrument (measuring accuracy: 0.2%, measuring range: 500 N.m) and an MIK-2740 temperature inspection instrument, which are connected to a computer to complete signal acquisition and measurement tasks, and Siemens S7-200 series PLC is used for interlocking control.

Comparison between Simulation and Physical Tests
According to the optimization results of the improved response surface method, the experimental time is 225 s, the axial air gap length is 15 mm, the thickness of the The optimal solution is that the axial air gap length is 15.00 mm, the input motor speed i 500.00 r/min, the thickness of the permanent magnet in the magnetizing direction is 13.0 mm, the thickness of the copper plate is 7.00 mm, and the number of fins of the heat dis sipation plate is 10, the length of the fins of the heat dissipation plate is 12.00 mm, and th minimum temperature of the heat dissipation plate is 41.8409 °C.

Physical Test
To verify the results of the theoretical analysis and finite element simulation men tioned above, the same parameters as those used in the simulation are selected for physi cal tests on a double disk magnetic coupler. The field test platform is set up as shown in Figure 8, and the infrared temperature sensor with a heat dissipation structure is shown in Figure 9.  The input motor is a 380 V three-phase low-speed motor, which is controlled by frequency converter. The rated power of this test-bed is 55 kW, the maximum speed i 3000 r/min, and the recommended speed range is 0-2000 r/min. The measuring circuit i equipped with NCTES3000 torque and a speed measuring instrument (measuring accu racy: 0.2%, measuring range: 500 N.m) and an MIK-2740 temperature inspection instru ment, which are connected to a computer to complete signal acquisition and measuremen tasks, and Siemens S7-200 series PLC is used for interlocking control.

Comparison between Simulation and Physical Tests
According to the optimization results of the improved response surface method, th experimental time is 225 s, the axial air gap length is 15 mm, the thickness of th The input motor is a 380 V three-phase low-speed motor, which is controlled by a frequency converter. The rated power of this test-bed is 55 kW, the maximum speed is 3000 r/min, and the recommended speed range is 0-2000 r/min. The measuring circuit is equipped with NCTES3000 torque and a speed measuring instrument (measuring accuracy: 0.2%, measuring range: 500 N.m) and an MIK-2740 temperature inspection instrument, which are connected to a computer to complete signal acquisition and measurement tasks, and Siemens S7-200 series PLC is used for interlocking control.

Comparison between Simulation and Physical Tests
According to the optimization results of the improved response surface method, the experimental time is 225 s, the axial air gap length is 15 mm, the thickness of the permanent magnet in the magnetizing direction is 13.00 mm, the thickness of the copper plate is 8.00 mm, the number of fins of the heat dissipation plate is 10, and the length of the fins of the heat dissipation plate is 12.00 mm, and the input motor speed is changed to 300 r/min, 400 r/min, 500 r/min and 600 r/min in turn. The physical test results are shown in Figure 10. The simulation value is compared with the test value, as shown in Table 7. permanent magnet in the magnetizing direction is 13.00 mm, the thickness of the copper plate is 8.00 mm, the number of fins of the heat dissipation plate is 10, and the length of the fins of the heat dissipation plate is 12.00 mm, and the input motor speed is changed to 300 r/min, 400 r/min, 500 r/min and 600 r/min in turn. The physical test results are shown in Figure 10. The simulation value is compared with the test value, as shown in Table 7.  It can be observed from Figure 10 that when the double disk magnetic coupler starts up, the temperature of the heat dissipation plate gradually rises, and it reaches a stable running state after 30 s. The five curves in the figure show that the temperature of the heat dissipation plate increases with the increase in the input motor speed. By comparing the simulation results in Table 7 with the physical test results, it can be observed that the simulation results are in good agreement with the physical test results, and the maximum error is 4.8%. The optimal parameter combination obtained after the correction of No.3 in Table 7 is selected, and the error between the simulation result and the physical test result is 1.9%, which indicates that the optimized parameters meet the expected effect of the actual transmission, and can be used to guide the design of the heat dissipation structure of the double disk magnetic coupler in the permanent magnet eddy current drive system of a coal mine.  It can be observed from Figure 10 that when the double disk magnetic coupler starts up, the temperature of the heat dissipation plate gradually rises, and it reaches a stable running state after 30 s. The five curves in the figure show that the temperature of the heat dissipation plate increases with the increase in the input motor speed. By comparing the simulation results in Table 7 with the physical test results, it can be observed that the simulation results are in good agreement with the physical test results, and the maximum error is 4.8%. The optimal parameter combination obtained after the correction of No.3 in Table 7 is selected, and the error between the simulation result and the physical test result is 1.9%, which indicates that the optimized parameters meet the expected effect of the actual transmission, and can be used to guide the design of the heat dissipation structure of the double disk magnetic coupler in the permanent magnet eddy current drive system of a coal mine.

Conclusions
Based on an orthogonal experimental design method, multi-parameter optimization of a heat dissipation structure of a double disk magnetic coupler is carried out, and the temperature of the heat dissipation plate is obtained by the finite element method. The influences of the axial air gap length of the double disk magnetic coupler, the input motor speed, the thickness of the permanent magnet in the magnetizing direction, the thickness of the copper plate, the number of fins of the heat dissipation plate and the length of the fins of the heat dissipation plate on the temperature of the heat dissipation plate are analyzed. Furthermore, the optimal calculation results are selected by the orthogonal experimental design, and the structural parameters of the model are improved and optimized; thus. the heat dissipation performance of the model is improved. The following conclusions are drawn:

1.
A six-factor and three-level experiment based on the improved response surface method is designed, and it is found that different variable factors have a certain influence on the heat dissipation effect of the heat dissipation structure model of the double disk magnetic coupler. After analyzing its heat dissipation performance, it can be observed that the order of significant influence from high to low is as follows: axial air gap length > input motor speed > the length of the fins of the heat dissipation plate > the thickness of the permanent magnet in the magnetizing direction > the number of fins of the heat dissipation plate > the thickness of the copper plate. The increase in the axial air gap length can effectively reduce the temperature rise, and the maximum decrease can reach 9.76%, which has an obvious effect. According to the single-factor and double-factor influence change diagram of the temperature of the heat dissipation plate, the response of different influence factors can be intuitively expressed, which is helpful to find the optimal direction of the heat dissipation structure.

2.
The parameters of the optimized heat dissipation structure obtained by orthogonal experimental design are as follows: the axial air gap length is 15.00 mm, the input motor speed is 500.00 r/min, the thickness of the permanent magnet in the magnetizing direction is 13.00 mm, the thickness of the copper plate is 7.00 mm, the number of fins of the heat dissipation plate is 10, and the length of the fins of the heat dissipation plate is 12.00 mm. The corresponding value of the objective function, namely the minimum temperature of the heat dissipation plate, is 41.8409 • C. Under this heat dissipation structure, the optimum heat dissipation performance of the double disk magnetic coupler can be achieved.

3.
According to the influence law and optimal parameters of the temperature characteristic index of the heat dissipation plate, the physical experiment was designed. The heat dissipation performance of the theoretical and practical heat dissipation structure before and after optimization under actual working conditions was verified.
The simulation results are in good agreement with the physical test results, with a maximum error of 4.8%. The error between the simulation results of the optimal parameter combination and the physical test results after correction is only 1.9%, with high accuracy, which verifies the correctness of the numerical model. The research results have reference significance for guiding the optimization of the heat dissipation structure parameters of the double disk magnetic coupler.

Data Availability Statement:
The data that support the findings of this study are available from the corresponding author.