Design and Parameter Optimization of Double-Mosquito Combination Coils for Enhanced Anti-Misalignment Capability in Inductive Wireless Power Transfer Systems

: This paper proposes a novel double-mosquito combination (DMC) coil for inductive wireless power transfer (IPT) systems to improve their anti-misalignment capability. The DMC coil consists of a mosquito coil with single-turn spacing and a tightly wound close-wound coil. By superimposing the magnetic fields generated by both coils, a relatively uniform magnetic field distribution is achieved on the receiving coil plane. This approach addresses the challenges of significant output voltage fluctuations and reduced transmission efficiencies caused by coupling coil misalignments in conventional IPT systems. To further optimize the DMC coil, an interaction law between its parameters and the mutual inductance is established, setting the coil mutual inductance fluctuation rate as the optimization objective, and using the coil turn spacing, number of turns, and outer diameter as constraint conditions. The beetle antennae search algorithm (BAS) is employed to enhance the whale optimization algorithm (WOA), facilitating the adaptive optimization of the coil parameters. An experimental IPT system platform with a 50 mm transmission distance is developed to validate the robust anti-misalignment capability of the proposed coil. The results demonstrate that within a horizontal misalignment range of 50 mm, the system’s output voltage fluctuation rate stays below 7.4%, and the transmission efficiency remains above 83%.


Introduction
Grounded in the principles of electromagnetic induction, the inductive wireless power transfer (IPT) technology facilitates wireless power transfer by establishing electromagnetic coupling between transmitting and receiving coils [1][2][3].This innovative technology offers numerous benefits such as safety, flexibility, and convenience, making it a promising solution for a wide range of applications.These applications include electric vehicles [4], underwater operations [5], automatic guided vehicles (AGVs) [6,7], and implantable medical devices [8].
The AGV, a handling tool that follows a set path, can significantly cut labor costs, boost productivity, and is commonly employed in various industries, transportation, and other sectors.The use of IPT systems to transfer power to AGVs in real-world scenarios inevitably results in vehicular vibrations during operation, causing the position of the receiving coil, mounted on the vehicle chassis, to shift in relation to the stationary transmitting coil embedded in the road surface.This displacement leads to changes in the IPT system's mutual inductance [9], affecting both the output power and transmission efficiency.Consequently, it is crucial to prioritize the investigation of the coil structure design and parameter optimization for the coupling mechanism in order to enhance the IPT system's anti-misalignment performance [10,11].
Enhancing the anti-misalignment capabilities of IPT systems through coil structure design hinges on the construction of a uniform magnetic field [12].A transmitting coil with a uniform magnetic field distribution enables the receiving coil to be positioned freely, maintaining the stability of the output power and transmission efficiency.Numerous transmitting coil design methods have been proposed by researchers worldwide to achieve a uniform magnetic field distribution.The literature [13] employs partial magnetic field cancellation and introduces a transmitting coil structure that combines forward and reverse winding, resulting in a uniform longitudinal magnetic field intensity distribution on the receiving coil plane.The literature [14] minimizes the mutual inductance fluctuation rate within the misalignment range by connecting a reverse-wound coil in series at the transmitting end.However, both [13,14] rely on the enumeration method, which entails lengthy simulation calculations and complex optimization processes.
The literature [15] suggests a group series winding-transmitting coil structure that involves winding two close-wound coil groups in series.This design achieves a uniform transverse magnetic flux generated by the transmitting coil when its outer diameter is larger than that of the receiving coil, but it only considers a single close-wound coil combination.The literature [16] proposes a combined series winding hexagonal coil structure that merges tightly and loosely wound coils.Although this approach improves the coupling coefficients and anti-misalignment abilities of wireless power transfer systems, it merely combines mosquito and close-wound coils without further examining the impact of different turnspacing coil arrangement methods on the coil's magnetic induction intensity distribution.Furthermore, it does not consider the optimal combination method for achieving the most uniform magnetic field distribution.
This paper aims to achieve a relatively uniform magnetic field distribution by combining close-wound and mosquito coils and determining the best combination method through a comparison of the coils' magnetic induction intensity distributions generated under various arrangement methods.In addition, to minimize mutual inductance fluctuations during the coupling coil misalignment process and to enhance the output power and transmission efficiency, the paper incorporates the tentacle algorithm to improve the whale algorithm, thus deriving the optimal coil turn spacing and the number of turns under given constraint conditions.The main contributions of this paper are listed as follows: (1) A proposed double-mosquito combination (DMC) coil structure combines the magnetic induction intensity distribution analysis of the mosquito coil and the closewound coil.This coil exhibits a relatively uniform magnetic field distribution, which can reduce the mutual inductance and output voltage fluctuation rate of the IPT system after deviation.(2) The integration of the beetle antennae search algorithm into the whale algorithm is presented to calculate the novel optimal whale position using the position update formula of the beetle antennae search algorithm.Compared with the optimal whale position calculated by the conventional whale algorithm, the smaller one is selected as the optimal position within the current whale population, thereby further enhancing the convergence speed and accuracy of the whale algorithm.(3) With the minimization of the coil's mutual inductance fluctuation rate as the optimization target, the parameters of the DMC coil are optimized using the improved whale algorithm, which enhances the resistance to misalignment of the DMC coil and ensures the efficient and stable operation of the IPT system.
The rest of the paper is organized as follows: Section 2 presents an equivalent circuit diagram of the IPT system, combining the coil, and establishes the relationships among system output power, transmission efficiency, and mutual inductance.In Section 3, the magnetic field distribution of the mosquito and close-wound coils is obtained through simulation.A model of two single-turn square coils coaxially aligned is developed, and the interaction law between the coil parameters and mutual inductance of transmitting and receiving coils is analyzed.Section 4 details the DMC coil design method, optimizes the DMC coil parameters using a specific algorithm, and examines the misalignment characteristics of the optimized DMC coil.Section 5 describes the construction of an experimental platform for the IPT system, providing verification of the robust anti-misalignment capability of the DMC coil after parameter optimization within a horizontal misalignment range of ±50 mm at a transmission distance of 50 mm.Finally, Section 6 concludes the paper.

Transmission Characteristics Analysis of IPT System with Combined Coils
At present, there are mainly two types of traditional planar square coils: mosquito coils with single-turn spacing and close-wound coils with tight winding.Figure 1 illustrates the comparison between the combined coil suggested in this study and the conventional square coil structure.
simulation.A model of two single-turn square coils coaxially aligned is developed, and the interaction law between the coil parameters and mutual inductance of transmitting and receiving coils is analyzed.Section 4 details the DMC coil design method, optimizes the DMC coil parameters using a specific algorithm, and examines the misalignment characteristics of the optimized DMC coil.Section 5 describes the construction of an experimental platform for the IPT system, providing verification of the robust anti-misalignment capability of the DMC coil after parameter optimization within a horizontal misalignment range of ±50 mm at a transmission distance of 50 mm.Finally, Section 6 concludes the paper.

Transmission Characteristics Analysis of IPT System with Combined Coils
At present, there are mainly two types of traditional planar square coils: mosquito coils with single-turn spacing and close-wound coils with tight winding.Figure 1  The combined coil integrates the winding features of both the mosquito coil and the close-wound coil, resulting in a coil with varying turn spacings from the inner to the outer section using a single wire.Based on its winding characteristics, the circuit diagram for the IPT system is illustrated in Figure 2. The IPT system primarily consists of a high-frequency inverter (I), compensation network (II), magnetic coupling mechanism (III), rectifier filter circuit (IV), and load (V) [17,18].Furthermore, the system employs a straightforward structure and a more stable power factor series-series (SS) compensation network.The combined coil integrates the winding features of both the mosquito coil and the close-wound coil, resulting in a coil with varying turn spacings from the inner to the outer section using a single wire.Based on its winding characteristics, the circuit diagram for the IPT system is illustrated in Figure 2. The IPT system primarily consists of a high-frequency inverter (I), compensation network (II), magnetic coupling mechanism (III), rectifier filter circuit (IV), and load (V) [17,18].Furthermore, the system employs a straightforward structure and a more stable power factor series-series (SS) compensation network.
simulation.A model of two single-turn square coils coaxially aligned is developed, and the interaction law between the coil parameters and mutual inductance of transmitting and receiving coils is analyzed.Section 4 details the DMC coil design method, optimizes the DMC coil parameters using a specific algorithm, and examines the misalignment characteristics of the optimized DMC coil.Section 5 describes the construction of an experimental platform for the IPT system, providing verification of the robust anti-misalignment capability of the DMC coil after parameter optimization within a horizontal misalignment range of ±50 mm at a transmission distance of 50 mm.Finally, Section 6 concludes the paper.

Transmission Characteristics Analysis of IPT System with Combined Coils
At present, there are mainly two types of traditional planar square coils: mosquito coils with single-turn spacing and close-wound coils with tight winding.Figure 1 illustrates the comparison between the combined coil suggested in this study and the conventional square coil structure.The combined coil integrates the winding features of both the mosquito coil and the close-wound coil, resulting in a coil with varying turn spacings from the inner to the outer section using a single wire.Based on its winding characteristics, the circuit diagram for the IPT system is illustrated in Figure 2. The IPT system primarily consists of a high-frequency inverter (I), compensation network (II), magnetic coupling mechanism (III), rectifier filter circuit (IV), and load (V) [17,18].Furthermore, the system employs a straightforward structure and a more stable power factor series-series (SS) compensation network.In Figure 2, Uin represents the DC input voltage, while Ut signifies the inverter output voltage.It and Ir denote the currents flowing through the transmitting and receiving coils, respectively.Ct stands for the compensation capacitor at the transmitter end, and Cr corresponds to the compensation capacitor at the receiver end.Rt is the internal resistance of In Figure 2, U in represents the DC input voltage, while U t signifies the inverter output voltage.I t and I r denote the currents flowing through the transmitting and receiving coils, respectively.C t stands for the compensation capacitor at the transmitter end, and C r corresponds to the compensation capacitor at the receiver end.R t is the internal resistance of the transmitting coil, and U r indicates the rectifier input voltage.L 1 , L 2 , and L n are the coil self-inductances for the first, second, and nth turn spacings, respectively, with the sum of the coil self-inductances of n turn spacings constituting the transmitting coil self-inductance.R r refers to the internal resistance of the receiving coil, while L r represents the receiving coil self-inductance.M 1r , M 2r , and M nr are the mutual inductances between the coil and the receiving coil for the first, second, and nth turn spacings, respectively.C F is the filter capacitor of the rectifier circuit, R L denotes the load resistance, U O is the voltage across the load, and R eq represents the equivalent load resistance.
According to Kirchhoff's law, we can derive: When the system reaches the resonant state, it meets the following equation: The transmitting circuit current, denoted by I t , is given as follows: The receiving circuit current, denoted by I r , is given as follows: From this, the system output power and transmission efficiency can be obtained as follows: The influential relationship of the mutual inductance M tr on the system output power and transmission efficiency is plotted as shown in Figure 3, based on Equations ( 5) and (6), which represent the superposition of mutual inductance between the transmitting coil and the receiving coil, that is, M tr = M 1r + M 2r + . . .+ M nr .
Figure 3 demonstrates that as the mutual inductance increases, the system output power initially rises and then declines.Concomitantly, the system transmission efficiency increases with the growth of mutual inductance, eventually reaching a plateau after surpassing a specific value.It is important to note that M tr is linked to both the magnetic flux (ϕ) and the magnetic induction intensity (B).Therefore, through the design and optimization of the magnetic coupling mechanism, achieving a uniform magnetic field is possible, leading to an enhancement in the system's transmission performance.
Electronics 2024, 13, x FOR PEER REVIEW 5 of 19 (6), which represent the superposition of mutual inductance between the transmitting coil and the receiving coil, that is, Mtr = M1r + M2r + …+Mnr.Figure 3 demonstrates that as the mutual inductance increases, the system output power initially rises and then declines.Concomitantly, the system transmission efficiency increases with the growth of mutual inductance, eventually reaching a plateau after surpassing a specific value.It is important to note that Mtr is linked to both the magnetic flux (ϕ) and the magnetic induction intensity (B).Therefore, through the design and optimization of the magnetic coupling mechanism, achieving a uniform magnetic field is possible, leading to an enhancement in the system's transmission performance.

Magnetic Induction Intensity Distribution of Combined Coils
The magnetic induction intensity distribution characteristics of both mosquito coils and close-wound coils were investigated by obtaining the magnetic field distribution on the receiving coil plane at varying heights (h) for these coil types via simulation, as shown in Figures 4 and 5

Magnetic Induction Intensity Distribution of Combined Coils
The magnetic induction intensity distribution characteristics of both mosquito coils and close-wound coils were investigated by obtaining the magnetic field distribution on the receiving coil plane at varying heights (h) for these coil types via simulation, as shown in Figures 4 and 5.
(6), which represent the superposition of mutual inductance between the transmitting coil and the receiving coil, that is, Mtr = M1r + M2r + …+Mnr.Figure 3 demonstrates that as the mutual inductance increases, the system output power initially rises and then declines.Concomitantly, the system transmission efficiency increases with the growth of mutual inductance, eventually reaching a plateau after surpassing a specific value.It is important to note that Mtr is linked to both the magnetic flux (ϕ) and the magnetic induction intensity (B).Therefore, through the design and optimization of the magnetic coupling mechanism, achieving a uniform magnetic field is possible, leading to an enhancement in the system's transmission performance.

Magnetic Induction Intensity Distribution of Combined Coils
The magnetic induction intensity distribution characteristics of both mosquito coils and close-wound coils were investigated by obtaining the magnetic field distribution on the receiving coil plane at varying heights (h) for these coil types via simulation, as shown in Figures 4 and 5   Using the center of the mosquito coil and close-wound coil as the origin, a straight line was drawn through the points (−150, 0, h) mm and (150, 0, h) mm.The normalized curve of the magnetic induction intensity along this line for different heights of the mosquito coil and close-wound coil is presented in Figure 6.Using the center of the mosquito coil and close-wound coil as the origin, a straight line was drawn through the points (−150, 0, h) mm and (150, 0, h) mm.The normalized curve of the magnetic induction intensity along this line for different heights of the mosquito coil and close-wound coil is presented in Figure 6.Using the center of the mosquito coil and close-wound coil as the origin, a straight line was drawn through the points (−150, 0, h) mm and (150, 0, h) mm.The normalized curve of the magnetic induction intensity along this line for different heights of the mosquito coil and close-wound coil is presented in Figure 6.The magnetic induction intensity is most prominent in the central region of the mosquito coil, as depicted in Figures 4-6.The magnetic field gradually diminishes from the coil center outward.In contrast, the magnetic induction intensity is concentrated on both sides of the close-wound coil, resulting in a decrease in the magnetic field in the central area.This is evident from the comparative analysis of the two coil types as illustrated in the figures.
This paper proposes a combined coil that merges the characteristics of the mosquito coil and close-wound coil, with the center of the combined coil comprising the mosquito coil and the outer portion made up of the close-wound coil.The inner and outer coils are connected in the same direction to form a combined coil.This design aims to address the issues of a weak magnetic field on the outer part of the mosquito coil and a low magnetic field intensity at the center of the close-wound coil, ultimately achieving a more uniform magnetic field-intensity distribution based on the magnetic field distribution characteristics of both coil types.The magnetic induction intensity is most prominent in the central region of the mosquito coil, as depicted in Figures 4-6.The magnetic field gradually diminishes from the coil center outward.In contrast, the magnetic induction intensity is concentrated on both sides of the close-wound coil, resulting in a decrease in the magnetic field in the central area.This is evident from the comparative analysis of the two coil types as illustrated in the figures.
This paper proposes a combined coil that merges the characteristics of the mosquito coil and close-wound coil, with the center of the combined coil comprising the mosquito coil and the outer portion made up of the close-wound coil.The inner and outer coils are connected in the same direction to form a combined coil.This design aims to address the issues of a weak magnetic field on the outer part of the mosquito coil and a low magnetic field intensity at the center of the close-wound coil, ultimately achieving a more uniform magnetic field-intensity distribution based on the magnetic field distribution characteristics of both coil types.

Mutual Inductance Calculation of Combined Coils
A schematic diagram showing the horizontal misalignment of a single-turn square coil is depicted in Figure 7.

Mutual Inductance Calculation of Combined Coils
A schematic diagram showing the horizontal misalignment of a single-turn square coil is depicted in Figure 7.The transmitting coil and the receiving coil consist of two coaxially aligned coils.In this arrangement, h represents the distance between the transmitting and receiving coils, l is the side length of the transmitting coil, r is the side length of the receiving coil, and Δ is the distance the receiving coil shifts in the positive direction along the x-axis.The re- The transmitting coil and the receiving coil consist of two coaxially aligned coils.In this arrangement, h represents the distance between the transmitting and receiving coils, l is the side length of the transmitting coil, r is the side length of the receiving coil, and ∆ is the distance the receiving coil shifts in the positive direction along the x-axis.The receiving coil's sides are labeled as b 1 , b 2 , b 3 , and b 4 , while the transmitting coil's four sides are a 1 , a 2 , a 3 , and a 4 .The mutual inductance M of the square coil can be considered as the superposition of the mutual inductance between two parallel sides [19].This setup defines the key parameters and relationships essential for understanding the interaction between the transmitting and receiving coils.
After defining the coordinates of any point on sides a 1 and b 1 as (l/2, y 1 , 0) and (r/2, y 2 , h), respectively, and the coordinates of any point on sides a 2 and b 2 as (x 1 , l/2, 0) and (x 2 , r/2, h), respectively, the M a1b1 and M a2b2 values following a horizontal misalignment ∆ undergone by the receiving coil can be determined in accordance with the Neumann formula as follows [20]: In the formula, µ 0 represents the magnetic permeability of a vacuum, where Similarly, the M a1b3 , M a2b4 , M a3b3 , M a3b1 , M a4b4 , and M a4b2 can be obtainedfrom a horizontal misalignment ∆ between the transmitting and receiving coils.
The M a1b1 , M a1b3 , M a2b2 , M a2b4 , M a3b3 , M a3b1 , M a4b4 , and M a4b2 can be obtained when the horizontal misalignment ∆ is 0 and the transmitting coil and the receiving coil are aligned coaxially.
In light of Equations ( 8) and ( 9), the mutual inductance for the multi-turn mosquito coil can be derived under both alignment and misalignment scenarios.Here, we denote the number of turns of the mosquito coil as N m and the coil turn spacing as d.Given that the turn spacing d of the mosquito coil is considerably larger than the diameter D of the single-turn coil, D can be disregarded when calculating the side length of the mosquito coil.Consequently, the side length for the N m -th turn of the mosquito coil can be determined as follows: By substituting Equation (10) into Equation ( 8), the M a1b1 of the N m -th turn of the mosquito coil can be obtained as follows: The total M a1b1(total) for N m turns can be obtained as follows: The mutual inductance M mr between the multi-turn mosquito coil and the square coil comprising N r turns can be expressed as follows: The diameter of the single-turn coil is denoted by D, and the number of turns for the close-wound coil is denoted by N c .Since the turn spacing d of the close-wound coil is approximately equal to 0, it can be disregarded when computing the side length of the close-wound coil.Thus, the side length for the N c -th turn of the close-wound coil can be described as follows: Similarly, drawing upon Equations ( 8) and ( 9), the M a1b1(Nc) for the N c -th turn of the close-wound coil, the M a1b1(Nc) for the N c -th turn following horizontal misalignment, and the mutual inductance M cr between the multi-turn, close-wound coil and the square coil with N r turns can be derived.
Therefore, the mutual inductance between the combined coil and the receiving coil can be represented as the superposition of the mutual inductances among the mosquito coil, close-wound coil, and receiving coil, respectively.The corresponding equation is given as follows: In Equation (15), n m denotes the total number of mosquito coils, k signifies the k-th mosquito coil counted from the innermost layer outward, n c represents the overall quantity of close-wound coils, and p refers to the p-th close-wound coil counted from the innermost layer outward.
In conclusion, a relatively uniform magnetic induction intensity distribution can be achieved in the combined coil by the rational design of the quantity and coil parameters for both mosquito and close-wound coils.This design effectively reduces the fluctuations in mutual inductance during the misalignment process.Such an optimization strategy for coil design ensures a consistent magnetic induction intensity distribution, thereby enhancing the overall performance and stability of the combined coil system.

Design of Double-Mosquito Combination (DMC) Coils
The combined coil, as discussed in Section 3.1, is comprised of a mosquito coil at the center and a close-wound coil at the outer layer.This DMC coil exhibits three different turn spacings and can be arranged in two possible configurations: (1) with the first and third layers being mosquito coils, and the second layer being a close-wound coil; and (2) with the first and second layers as mosquito coils, and the third layer as a close-wound coil.The magnetic field distribution of the combined coil under these two conditions was simulated using COMSOL Multiphysics 5.6 software.As the turn spacing of the close-wound coil is approximately zero, the diameter D of the single-turn coil is used instead of its turn spacing.The turn spacing and width of each layer of the coil are denoted by d 1 , W 1 ; d 2 , W 2 ; and d 3 , W 3 , respectively.The specific parameters are detailed in Table 1, in which No.  Using the center of the combined coil as the origin, a straight line passing through the points (−150 mm, 0 mm, h mm) and (150 mm, 0 mm, h mm) is drawn.Figure 8 illustrates the distribution curve of the magnetic flux density along this line under different arrangement methods.
wound coil is approximately zero, the diameter D of the single-turn coil is used instead of its turn spacing.The turn spacing and width of each layer of the coil are denoted by d1, W1; d2, W2; and d3, W3, respectively.The specific parameters are detailed in Table 1 Using the center of the combined coil as the origin, a straight line passing through the points (−150 mm, 0 mm, h mm) and (150 mm, 0 mm, h mm) is drawn.Figure 8   As shown in Figure 8, the magnetic flux density distribution of the combined coil under the first arrangement method is high in the middle and low on both sides.The second arrangement method results in a more uniform distribution.Hence, to achieve a uniform magnetic flux density on the receiving coil plane of the combined coil, the second arrangement method is adopted, in which the turn spacing decreases gradually from the inside to the outside.The final schematic diagram of the DMC coil is illustrated in Figure 9.As shown in Figure 8, the magnetic flux density distribution of the combined coil under the first arrangement method is high in the middle and low on both sides.The second arrangement method results in a more uniform distribution.Hence, to achieve a uniform magnetic flux density on the receiving coil plane of the combined coil, the second arrangement method is adopted, in which the turn spacing decreases gradually from the inside to the outside.The final schematic diagram of the DMC coil is illustrated in Figure 9.As shown in Figure 9, the first and second layers are mosquito coils, and the third layer is a close-wound coil.The inter-turn distances of the three layers of coils are denoted by d1, d2, and d3, respectively.The widths of the three layers of coils are denoted by W1, W2, and W3, respectively.The number of turns of the three layers of coils are denoted by N1, N2, and N3, respectively.The total width of the coil is denoted by S. As shown in Figure 9, the first and second layers are mosquito coils, and the third layer is a close-wound coil.The inter-turn distances of the three layers of coils are denoted by d 1 , d 2 , and d 3 , respectively.The widths of the three layers of coils are denoted by W 1 , W 2 , and W 3 , respectively.The number of turns of the three layers of coils are denoted by N 1 , N 2 , and N 3 , respectively.The total width of the coil is denoted by S.
In this paper, a DMC coil is fabricated using Litz wire with a single-turn coil diameter D of 2 mm.When the inter-turn spacing d of the coil is greater than D, the diameter of the Litz wire is negligible compared to its inter-turn spacing.Conversely, when d is less than or equal to D, the inter-turn spacing is negligible compared to the diameter of the Litz wire.Hence, the calculation formula for the DMC coil parameters is given by:

Optimization of Double-Mosquito Combination (DMC) Coil Parameters
The improved whale optimization algorithm (WOA) is used to optimize the parameters of the DMC coil.Given the total width S of the coil and the transmission distance h, the optimization determines the inter-turn distance d, width W, and number of turns N of the three-layer coil.This ensures a relatively uniform magnetic field distribution and a relatively stable mutual inductance for the DMC coil.
The whale optimization algorithm (WOA) is inspired by the hunting behavior of humpback whales in nature [21].It has a simple and easy-to-implement principle, and it is widely applied to various optimization problems.In the WOA, the position of the optimal whale affects the position update of other whale individuals.To enhance the convergence speed and accuracy of the WOA, the beetle antennae search algorithm (BAS) [22] is employed to improve the position of the optimal whale.BAS searches on both sides of the optimal whale obtained using the WOA, in the domain of the side with a lower fitness value, calculates the new optimal whale position using the position update formula of the BAS, compares the optimal whale positions of the two algorithms, and selects the one with a lower fitness value as the leading whale.
Based on Equation ( 13), assuming that the total width of the coil S, the transmission distance h, and the horizontal misalignment distance ∆ are constant, the objective function is the minimum mutual inductance fluctuation rate before and after the misalignment ∆, and the inter-turn spacing, width, and number of turns of the three-layer coil are the key parameters affecting the mutual inductance fluctuation rate.From Equation ( 15), d 1 , d 2 , N 1 , and N 2 are known, and other coil parameters can be derived.Hence, the objective function of the DMC coil is: where M represents the mutual inductance of aligned coils, and M ∆ represents the mutual inductance after the coils are displaced.S is the length of the outermost side of the coil, which is 300 mm.To prevent the DMC coil from exceeding S due to the large inter-turn distance of the first and second layer coils, the inter-turn distance of the first layer coil is limited to 25 mm or less.The inter-turn distance of the second layer coil is smaller than that of the first layer coil and is limited to 15 mm or less.To reduce the cost, the number of turns of the third layer coil is limited to 29 or less.The final optimization constraints for the coil parameters are: The DMC coil parameter optimization process is illustrated in Figure 10.The DMC coil parameter optimization process is illustrated in Figure 10.The optimization steps are as follows: • Initialize the DMC coil parameters (d 1 , d 2 , N 1 , and N 2 ) and randomly generate a whale population (d 1 , d 2 , N 1 , and N 2 ).

•
Calculate the objective function f (x) for each whale in the population, and mark the whale with the minimum objective function value as X*.

•
Update the position parameters of the whales and check whether they exceed the search boundary.If so, adjust them to the boundary value.

•
Calculate the objective function value of the whales after the position update, and obtain the optimal whale position X* WOA .

•
Based on the position update mechanism of the beetle antennae search algorithm (BAS), calculate the new optimal whale position X* BAS , and compare it with X* WOA .Choose the one with the lower objective function value as the global optimal whale position X*.

•
Update the beetle antennae parameters.

•
If the number of iterations reaches or exceeds T max , output the optimal whale position, which is the optimal solution (d 1 , d 2 , N 1 , and N 2 ).Otherwise, update the whale position parameters and continue the process.
To verify the effectiveness of the improved algorithm, this paper compares the improved whale algorithm and WOA for parameter optimization.Figure 11 shows the comparison of the algorithm convergence results.
Figure 11 shows that the improved WOA and unimproved WOA have a similar convergence accuracy, but the improved WOA has a faster convergence speed and smaller objective function value.This implies that the beetle antennae search algorithm can effectively improve the WOA performance.
Table 2 shows the optimized coil parameters.
• If the number of iterations reaches or exceeds Tmax, output the optimal whale position, which is the optimal solution (d1, d2, N1, and N2).Otherwise, update the whale position parameters and continue the process.
To verify the effectiveness of the improved algorithm, this paper compares the improved whale algorithm and WOA for parameter optimization.Figure 11 shows the comparison of the algorithm convergence results.Figure 11 shows that the improved WOA and unimproved WOA have a similar convergence accuracy, but the improved WOA has a faster convergence speed and smaller objective function value.This implies that the beetle antennae search algorithm can effectively improve the WOA performance.
Table 2 shows the optimized coil parameters.

Analyzing the Misalignment Features of the Optimized Double-Mosquito Combination (DMC) Coil
This paper studies the misalignment characteristics of the optimized DMC coil, using COMSOL Multiphysics to verify its uniform magnetic field distribution and small mutual inductance fluctuation.Figure 12 shows the distribution of the magnetic field intensity of

Analyzing the Misalignment Features of the Optimized Double-Mosquito Combination (DMC) Coil
This paper studies the misalignment characteristics of the optimized DMC coil, using COMSOL Multiphysics to verify its uniform magnetic field distribution and small mutual inductance fluctuation.Figure 12 shows the distribution of the magnetic field intensity of the DMC coil (as the transmitting coil) on different height planes by changing the height h between the receiving coil and the transmitting coil.As shown in Figure 12, the magnetic induction intensity of the DMC coil is approximately uniform within the square region of [−80 mm, 80 mm] for h = 40 mm and h = 70 mm, within the square region of [−100 mm, 100 mm] for h = 50 mm, and within the square region of [−90 mm, 90 mm] for h = 60 mm.Although the areas of approximately equal magnetic field strengths vary at different heights for the DMC coil, it is evident that the magnetic field fluctuates minimally over a wide range.
Figure 13 shows the law of mutual inductance change under the horizontal misalignment of the transmitting coil at different h values.The specific data are shown in Table 3.As shown in Figure 12, the magnetic induction intensity of the DMC coil is approximately uniform within the square region of [−80 mm, 80 mm] for h = 40 mm and h = 70 mm, within the square region of [−100 mm, 100 mm] for h = 50 mm, and within the square region of [−90 mm, 90 mm] for h = 60 mm.Although the areas of approximately equal magnetic field strengths vary at different heights for the DMC coil, it is evident that the magnetic field fluctuates minimally over a wide range.
Figure 13 shows the law of mutual inductance change under the horizontal misalignment of the transmitting coil at different h values.The specific data are shown in Table 3.
mm, within the square region of [−100 mm, 100 mm] for h = 50 mm, and within the square region of [−90 mm, 90 mm] for h = 60 mm.Although the areas of approximately equal magnetic field strengths vary at different heights for the DMC coil, it is evident that the magnetic field fluctuates minimally over a wide range.
Figure 13 shows the law of mutual inductance change under the horizontal misalignment of the transmitting coil at different h values.The specific data are shown in Table 3.Based on Figure 13 and Table 3, the mutual inductance fluctuation of the DMC coil during the horizontal misalignment process is small for different transmission distances,  Based on Figure 13 and Table 3, the mutual inductance fluctuation of the DMC coil during the horizontal misalignment process is small for different transmission distances, and the maximum fluctuation rate is less than 8.2%.This demonstrates the excellent anti-misalignment performance of the DMC coil at different transmission distances.

Experimental Verification
An experimental platform based on the IPT system shown in Figure 2a was built, as depicted in Figure 14.The coupled coils and the receiving coil have an outer diameter of 142 mm × 142 mm, while the DMC coil, the mosquito coil, and the close-wound coil have an outer diameter of 300 mm × 300 mm.A was used as the controller in this experiment to collect the signals and generate the PWM waves.The transmitting and receiving coils of the magnetic coupling mechanism were wound with 0.1 mm × 200 strands of Litz wire.The coil mutual inductance was measured using a Chroma 11050-5M high-frequency LCR tester.The voltage at both ends of the load was measured using a Tektronix MSO 2024B oscilloscope, and the load current was measured using an RP1001C current probe.
The system input voltage was set to 45 V, the operating frequency f was 85 kHz, the self-inductance L r of the receiving coil was 60.3 µH, the compensation capacitance C r of the receiving coil was 58.2 nF, and the load R L was 10 Ω.According to the actual coil parameters, the compensation circuit parameters were adjusted.The mutual inductances, system output voltages, output powers, and efficiencies of the three transmitting coils with the same receiving coil were compared.The self-inductance measurements of the three transmitting coils and the system compensation capacitance values are shown in Table 4.
an outer diameter of 300 mm × 300 mm.A was used as the controller in this experiment to collect the signals and generate the PWM waves.The transmitting and receiving coils of the magnetic coupling mechanism were wound with 0.1 mm × 200 strands of Litz wire.The coil mutual inductance was measured using a Chroma 11050-5M high-frequency LCR tester.The voltage at both ends of the load was measured using a Tektronix MSO 2024B oscilloscope, and the load current was measured using an RP1001C current probe.The system input voltage was set to 45 V, the operating frequency f was 85 kHz, the self-inductance Lr of the receiving coil was 60.3 µH, the compensation capacitance Cr of the receiving coil was 58.2 nF, and the load RL was 10 Ω.According to the actual coil parameters, the compensation circuit parameters were adjusted.The mutual inductances, system output voltages, output powers, and efficiencies of the three transmitting coils with the same receiving coil were compared.The self-inductance measurements of the three transmitting coils and the system compensation capacitance values are shown in Table 4. Using the position where the receiving coil and the transmitting coil are coaxially aligned and the transmission distance is 50 mm as the reference base, we measured the  Using the position where the receiving coil and the transmitting coil are coaxially aligned and the transmission distance is 50 mm as the reference base, we measured the mutual inductance of the DMC, mosquito-coil, and close-wound coil transmitters during a 50mm horizontal misalignment of the receiving coil.The results are shown in Figure 15.Based on the experimental results, it can be observed that the mutual inductance fluctuation rate of the DMC and the receiving coil is 7.3% when the receiving coil undergoes a horizontal misalignment.This rate is significantly lower than the mutual inductance fluctuation rates of the close-wound coil (12%) and the mosquito coil (20.3%).Therefore, it can be concluded that the DMC coil is capable of ensuring the stable operation of the system under misalignment conditions.
The waveforms of the inverter output voltage Ut, inverter output current It, load terminal voltage Ul, and load current Il for the DMC, mosquito, and close-wound coils are shown in Figure 16 under conditions of alignment with the receiving coil or horizontal misalignment.Based on the experimental results, it can be observed that the mutual inductance fluctuation rate of the DMC and the receiving coil is 7.3% when the receiving coil undergoes a horizontal misalignment.This rate is significantly lower than the mutual inductance fluctuation rates of the close-wound coil (12%) and the mosquito coil (20.3%).Therefore, it can be concluded that the DMC coil is capable of ensuring the stable operation of the system under misalignment conditions.
The waveforms of the inverter output voltage U t , inverter output current I t , load terminal voltage U l , and load current I l for the DMC, mosquito, and close-wound coils are shown in Figure 16 under conditions of alignment with the receiving coil or horizontal misalignment.
Based on the experimental results, it can be observed that the mutual inductance fluctuation rate of the DMC and the receiving coil is 7.3% when the receiving coil undergoes a horizontal misalignment.This rate is significantly lower than the mutual inductance fluctuation rates of the close-wound coil (12%) and the mosquito coil (20.3%).Therefore, it can be concluded that the DMC coil is capable of ensuring the stable operation of the system under misalignment conditions.
The waveforms of the inverter output voltage Ut, inverter output current It, load terminal voltage Ul, and load current Il for the DMC, mosquito, and close-wound coils are shown in Figure 16  As shown in Figure 16, when RL = 10 Ω, the output voltage of the IPT system with the DMC coil was 29.89 V and the efficiency was 83.3% when the coils were aligned; the output voltage was 32.09 V and the efficiency was 83.13% when the coils were displaced by 50 mm; the output voltage fluctuation rate of the system during the displacement process was 7.4%.The output voltage of the mosquito coil was 38.45 V, and the efficiency was 84.45% when the coils were aligned; the output voltage was 47.47 V, and the efficiency was 84.1% when the coils were displaced by 50 mm; the output voltage fluctuation rate of the system during the displacement process was 23.46%.The output voltage of the closewound coil was 26.77 V, and the efficiency was 83.28% when the coils were aligned; the output voltage was 30.18V, and the efficiency was 82.86% when the coils were displaced by 50 mm; the output voltage fluctuation rate of the system during the displacement process was 12.7%.It can be seen that the output voltage fluctuation rates of the three coils are significantly different, while the efficiencies are similar.Therefore, under the condition of comparable efficiency, the output voltage fluctuation rate is the main criterion for comparison.Compared with the mosquito coil and the close-wound coil, the output voltage fluctuation rate of the DMC coil is significantly lower, and the system's robustness to displacement is greatly improved.
To verify the strong anti-misalignment capability of the DMC coil under different load conditions, RL was set to 20 Ω, 30 Ω, and 40 Ω respectively, and the transmission distance between the transmitting coil and the receiving coil was 50 mm.The output voltages and transmission efficiencies of the DMC, mosquito, and close-wound transmitting As shown in Figure 16, when R L = 10 Ω, the output voltage of the IPT system with the DMC coil was 29.89 V and the efficiency was 83.3% when the coils were aligned; the output voltage was 32.09 V and the efficiency was 83.13% when the coils were displaced by 50 mm; the output voltage fluctuation rate of the system during the displacement process was 7.4%.The output voltage of the mosquito coil was 38.45 V, and the efficiency was 84.45% when the coils were aligned; the output voltage was 47.47 V, and the efficiency was 84.1% when the coils were displaced by 50 mm; the output voltage fluctuation rate of the system during the displacement process was 23.46%.The output voltage of the close-wound coil was 26.77 V, and the efficiency was 83.28% when the coils were aligned; the output voltage was 30.18V, and the efficiency was 82.86% when the coils were displaced by 50 mm; the output voltage fluctuation rate of the system during the displacement process was 12.7%.It can be seen that the output voltage fluctuation rates of the three coils are significantly different, while the efficiencies are similar.Therefore, under the condition of comparable efficiency, the output voltage fluctuation rate is the main criterion for comparison.Compared with the mosquito coil and the close-wound coil, the output voltage fluctuation rate of the DMC coil is significantly lower, and the system's robustness to displacement is greatly improved.
pared with other coils in the literature.To ensure a consistent comparison, the outer diameter of the transmitting coil was set to be the same, and the receiving coil was the same closely wound coil.The coil structure schematic is shown in Figure 18.During the process of the receiving coil horizontally misaligning by 50 mm relative to the transmitting coil, the robustness to displacement of the IPT system with different transmitting coils was measured, as shown in Table 5.As can be seen from the Table 5, under the same horizontal misalignment distance, the mutual inductance, output voltage, and output power fluctuation rates of the DMC coil proposed in this paper were all the smallest during the misalignment process.Compared with the existing coils, the coil proposed in this paper exhibited a better resistance to the misalignment distance, thereby ensuring that an IPT system with this coil can work efficiently and stably when the coupling mechanism is misaligned.

Conclusions
This paper proposes a DMC coil, which combines the magnetic field distribution characteristics of a mosquito coil and a close-wound coil to address the high anti-misalignment performance problem of the IPT system.From the perspective of coil structure optimization and enhancing the uniformity of magnetic field distribution, the coil aims to achieve a relatively even magnetic field distribution and minimize the fluctuations in system mutual inductance and output power.In this study, an equivalent circuit model of the IPT system is employed to derive analytical expressions for the system's transmission power and efficiency with respect to mutual inductance.Furthermore, the paper investigates the interaction between DMC coil parameters and mutual inductance.To optimize the DMC coil parameters, we propose an improved whale optimization algorithm, which incorporates the beetle antennae search algorithm.This improved algorithm is anticipated to yield more effective optimization results for the DMC coil design.This study involves a 3D simulation of the DMC coil and an investigation into its misalignment characteristics.To assess its antimisalignment capability in comparison to other coil types, we established an experimental platform with a transmission distance of 50 mm.The experimental results reveal that the mutual inductance fluctuation rate of the DMC coil during the misalignment process is 13% lower than that of the mosquito coil and 4.7% lower than that of the close-wound coil with the same load and misalignment distance.Similarly, the output voltage fluctuation rate of the DMC coil is 16.06% lower than that of the mosquito coil and 5.3% lower than that of the close-wound coil.These findings illustrate a substantial decrease in both the mutual inductance and output voltage fluctuation rates of the DMC coil compared to other coil types.This suggests that an IPT system featuring the DMC coil can effectively provide stable power to an AGV even when there is a misalignment in the coupling coil.Based on this, to further enhance the coupling between the coils, ferrite cores can be added, and the impact of factors such as core layout and core material on the mutual inductance between the coils can be comprehensively considered in order to further refine the analysis model of the coupling mechanism and lay a theoretical foundation for the global optimization of the IPT system.

Figure 1 .
Figure 1.Structure comparison between traditional square coil and combined coil: (a) mosquito coil; (b) close-wound coil; and (c) combined coil.

Figure 2 .Figure 1 .
Figure 2. Circuit diagram for the IPT system: (a) system structure diagram and (b) equivalent circuit diagram.In Figure2, Uin represents the DC input voltage, while Ut signifies the inverter output voltage.It and Ir denote the currents flowing through the transmitting and receiving coils, respectively.Ct stands for the compensation capacitor at the transmitter end, and Cr corresponds to the compensation capacitor at the receiver end.Rt is the internal resistance of

Figure 1 .
Figure 1.Structure comparison between traditional square coil and combined coil: (a) mosquito coil; (b) close-wound coil; and (c) combined coil.

Figure 2 .
Figure 2. Circuit diagram for the IPT system: (a) system structure diagram and (b) equivalent circuit diagram.

Figure 2 .
Figure 2. Circuit diagram for the IPT system: (a) system structure diagram and (b) equivalent circuit diagram.

Figure 3 .
Figure 3.The relationship between output power, transmission efficiency, and mutual inductance. .

Figure 3 .
Figure 3.The relationship between output power, transmission efficiency, and mutual inductance.

Figure 3 .
Figure 3.The relationship between output power, transmission efficiency, and mutual inductance. .

Figure 7 .
Figure 7. Schematic diagram of the horizontal misalignment of a single-turn square coil.

Figure 7 .
Figure 7. Schematic diagram of the horizontal misalignment of a single-turn square coil.
1 and No.2 share the same turn spacing and width, as do No.3 and No.4 as well as No.5 and No.6.No.1, No.3, and No.5 correspond to the first type of arrangement, while No.2, No.4, and No.6 correspond to the second type of arrangement.
illustrates the distribution curve of the magnetic flux density along this line under different arrangement methods.

Figure 8 .
Figure 8. curve of magnetic induction intensity distribution on line.

Electronics 2024 , 19 Figure 8 .
Figure 8. Normalized curve of magnetic induction intensity distribution on line.

Figure 10 .Figure 10 .
Figure 10.Flow chart for optimization of coil parameters.The optimization steps are as follows:

Figure 11 .
Figure 11.The iterative curve of optimization.

Figure 11 .
Figure 11.The iterative curve of optimization.

Figure 13 .
Figure 13.Mutual inductance variation patterns in the case of horizontal misalignment.

Figure 13 .
Figure 13.Mutual inductance variation patterns in the case of horizontal misalignment.

Electronics 2024 ,
13, x FOR PEER REVIEW 15 of 19 mutual inductance of the DMC, mosquito-coil, and close-wound coil transmitters during a 50mm horizontal misalignment of the receiving coil.The results are shown in Figure 15.

Figure 15 .
Figure 15.Mutual inductance of coils with misalignment in the horizontal plane.

Figure 15 .
Figure 15.Mutual inductance of coils with misalignment in the horizontal plane.

Table 1 .
Parameters of mosquito coils and close-wound coils.

Table 1 .
, in which No.1 and No.2 share the same turn spacing and width, as do No.3 and No.4 as well as No.5 and No.6.No.1, No.3, and No.5 correspond to the first type of arrangement, while No.2, No.4, and No.6 correspond to the second type of arrangement.Parameters of mosquito coils and close-wound coils.

Table 3 .
Results of mutual sensing at different transmission distances.

Table 3 .
Results of mutual sensing at different transmission distances.

Table 4 .
Design values of system parameters.

Table 4 .
Design values of system parameters.

Table 5 .
Comparison of misalignment resistance in DMC systems.
ND* represents the turns in a single D coil.

Table 5 .
Comparison of misalignment resistance in DMC systems.
N D * represents the turns in a single D coil.
Author Contributions: Conceptualization, W.H. and J.H.; methodology, J.H.; software, Y.H. and J.H.; validation, J.H., Y.H. and Y.Z.; formal analysis, W.H. and J.H.; investigation, Y.C. and Y.Z.; resources, Y.C. and Y.Z.; writing-original draft preparation, J.H.; writing-review and editing, W.H. and J.H.; visualization, J.H.; supervision, Y.C.All authors have read and agreed to the published version of the manuscript.Funding: This research was funded by the National Natural Science Foundation of China under the grant number 61903129 and the Hubei University of Technology Doctoral Research Project under the grant number BSQD2020012.