Design and Non-Linearity Optimization of a Vertical Brushless Electric Power Steering Angle Sensor

This paper presents the design and the non-linearity optimization of a new vertical non-contact angle sensor based on the electromagnetic induction principle. The proposed sensor consists of a stator part (with one solenoidal excitation coil and three sinusoidal receiver coils) and a rotor part (with six rectangular metal sheets). The receiver coil was designed based on the differential principle, which eliminates the effect of the excitation coil on the induced voltage of the receiver coil, and essentially decouples the excitation field from the eddy current field. Moreover, the induced voltages in the three receiver coils are three-phase sinusoidal signals with a phase difference of 10°, which are linearized by CLARK transformation. To minimize the sensor non-linearity, the Plackett–Burman technique was used, which identified the stator radius and the rotor blade thickness as the key factors affecting the sensor linearity. Then, the particle swarm algorithm with decreasing inertia weights was utilized to optimize the sensor linearity. A sensor prototype was made and tested in the laboratory, where the experimental results showed that the sensor non-linearity was only 0.648% and 0.645% in the clockwise and counterclockwise directions, respectively. Notably, the non-linearity of the sensor was less than −0.696% at different speeds.


Introduction
Angle sensors are widely used in various applications where precise angle measurements are required, including automobiles [1], aviation servo systems [2], and industrial robots [3].Depending on the underlying angle measurement method, angle sensors can be divided into contact angle sensors [4] and non-contact angle sensors [5][6][7].Contact angle sensors are mainly the potentiometer-type angle sensors [8], where the presence of friction during measurements makes these sensors more prone to wear, electrical noise, short life, etc. Non-contact angle sensors, on the other hand, can be further divided into capacitive, photoelectric, and magnetoelectric angle sensors.Capacitive angle sensors have the advantage of low power consumption and high sensitivity [9][10][11], but this sensor has poor temperature stability and low measurement accuracy.Meanwhile, photoelectric angle sensors need an independent light source, and since environmental factors affect the light propagation, such sensors require tight sealing to ensure their use in harsh environments [12][13][14].Magnetoelectric angle sensors include Hall sensors, magnetostrictive sensors, and electromagnetic induction-type sensors.Hall-type angle sensors have the advantages of high reliability, long life, and fast response [15,16], but this sensor contains a specific structure of the internal permanent magnet, its interchangeability is relatively poor, and the output is non-linear.In addition, the signal of Hall-type sensors is affected by temperature, thus needing a temperature compensation device, and the sensor accuracy is relatively low.Alternatively, magnetostrictive angle sensors have simple installation, high sensitivity, good stability, high power, and high overload capacity [17,18].However, these sensors have a poor resistance to interference and cannot be used with magnetically conductive materials.In contrast, with the advantages of a good structural flexibility, low cost, no required temperature compensation, high environmental adaptability, high accuracy, and good electromagnetic compatibility [19][20][21], electromagnetic induction angle sensors are quite popular in many engineering applications.This paper mainly studies the angle sensor used in the Electric Power Steering (EPS) system.Measuring the angle is an important function of EPS sensors.The function of the EPS system is to detect the torque and direction generated by the steering wheel when the driver is steering.A critical component of the vehicle chassis system, the steering system directly influences the stability, driving comfort, and driving safety of the vehicle.The development of the steering system has progressed through four stages: (i) a Mechanical Steering system; (ii) a Hydraulic Power Steering system; (iii) an Electro-Hydraulic Power Steering; and (iv) an Electric Power Steering system.Among these developments, the EPS system has the advantage of energy saving and environmental protection, and is increasingly used in automobiles.
Existing electromagnetic induction type angle sensors have a flat structure [22][23][24][25], which mainly has the following defects: relatively large sensor size due to the planar windings; and the two rotors in the sensor exhibit a crosstalk and affect the sensor nonlinearity.These issues of planar induction sensors indicate the need to design a new induction angle sensor with a small size and a high linearity.The structural parameters of the electromagnetic induction type angle sensor affect the intensity of the magnetic induction inside the sensor, which in turn affects the linearity of the relationship between the measured angle and the output signal [26].Non-linearity is an important parameter to measure the performance of sensors [27].To effectively reduce the non-linearity of an angle sensor, it is necessary to optimize the structural parameters of the sensor.
Usually, an optimization algorithm, such as a genetic algorithm [28], a response surface method [29], or a particle swarm algorithm [30], is used for the sensor structure optimization.The selection of genetic algorithm parameters, such as mutation rate and crossover rate, seriously affects the quality of the corresponding optimal solution.At present, the selection of these parameters mostly depends on the experience.The response surface method is not suitable for discrete optimization, and the optimal parameters obtained by this method depend on the degree of fit of the regression equation.The Particle Swarm Optimization (PSO) algorithm was first proposed by Eberhart and Kennedy in 1995, the basic concept of which originates from the study of the foraging behavior of bird flocks.This algorithm is fast and efficient in searching, does not depend on the problem information, and has a strong generality.The PSO algorithm is a generalized swarm intelligence method, commonly used for solving global optimization problems.In this paper, we design a new vertical non-contact induction angle sensor, where we use the Plackett-Burman test to screen the key factors affecting the non-linearity of the sensor, and then use the PSO algorithm to optimize these selected key factors.

Theoretical Research and Simulation Analysis 2.1. Theoretical Research
The proposed vertical non-contact induction angle sensor consists of a stator and a rotor, as shown in Figure 1.The stator includes the receiving coil and the excitation coil.Furthermore, the receiving coil winding consists of three sinusoidally structured coils with a sinusoidal period of 60 • , where the spacing between the adjacent receiving coils is 10 • .The excitation coil is a 10-turn solenoidal coil, while the rotor consists of six rotor blades, each at 30 • .
The measurement principle of the proposed sensor can be considered as "electricmagnetic-electric".When the excitation coil is fed with a high-frequency alternating current, an alternating magnetic field ( → B e ) is generated, and according to the Faraday law of electromagnetic induction, the voltage is induced in the receiving coil.The rotor in the alternating magnetic field induces a secondary magnetic field ( → B r ) of the same frequency, and the eddy current field here also induces a voltage in the receiving coil.Since the rotor is distributed at intervals, the coupling area between the rotor and the receiving coil is related Sensors 2024, 24, 2469 3 of 18 to the rotor rotation angle, and the variation of the coupling area changes the induced voltage in the receiving coil.Therefore, the receiving coil induced voltage is due to the superposition of the excitation magnetic field → B e and the eddy current field → B r , as shown in Equation (1).Next, the magnetic induction intensity can be calculated as shown in Equation (2).
→ B e,r = µ 0 4π where N is the number of turns; S is the area; I is the excitation current; According to the law of electromagnetic induction as given in Equation ( 1), the induced voltage is related to the magnetic induction intensity and the magnetic flux through a closed loop area.In this work, the receiving coil is designed as a differential structure, as shown in Figure 2a.Specifically, the geometry of the two adjacent loops of the receiving coil is the same.The induced voltages are equal in magnitude and opposite in direction to the excitation magnetic field, so the effect of the excitation coil on the induced voltage of the receiving coil is zero.Evidently, the differential structure of the receiver coil eliminates the influence of the excitation field on the induced voltage of the receiver coil.It should be noted that the coupling area between the rotor and the receiver coil is the key factor affecting the sensor-induced voltage.According to Figure 2a, for the sake of illustration, the receiving coil can be regarded as a series connection of the primary coil and the secondary coil.The coupling area between the rotor and the receiving coil changes as follows: (a) when the rotor is at position 1, the rotor is in the middle of the loop consisting of the primary and secondary coils, and the direction of the voltage induced by the primary and secondary coils is demonstrated by the arrow in the figure.The induced voltages of sinusoidal loop 1 and sinusoidal loop 2 are of the same magnitude but with opposite directions (the direction of the induced voltage in sinusoidal loop 1 is indicated by a red arrow, while that in sinusoidal loop 2 is indicated by a yellow arrow).Hence, they cancel out each other, and the induced voltage of the receiving coil is zero, i.e., the zero position in Figure 2b; (b) as the rotor rotates, the coupling area between the rotor and sinusoidal loop 1 gradually decreases, while the coupling area between the rotor and sinusoidal loop 2 gradually increases.Likewise, the magnetic flux through sinusoidal loop 2 gradually increases, and the induced voltage also gradually increases.When the rotor rotates to position 2, the rotor and sinusoidal loop 2 are fully coupled, and at this point, the magnetic flux through the sinusoidal loop 2 reaches the maximum.Meanwhile, the coupling area between the rotor and the sinusoidal loop 1 is 0, hence the magnetic flux within sinusoidal loop 1 is also zero.Here, the induced voltage of the receiving coil is at the maximum, as shown by S max in Figure 2b; (c) as the rotor rotates further, the coupling area between the rotor and sinusoidal loop 2 gradually decreases, whereas the coupling area between the rotor and sinusoidal loop 3 gradually increases, and the induced voltage slowly decreases.When the rotor is again in the middle of the loop consisting of the primary and secondary coils, the induced voltage in the receiving coil is zero; (d) with more rotation of the rotor, the receiving coil induced voltage reaches a negative maximum, where the rotor is fully coupled to sinusoidal loop 3, corresponding to position −S max in Figure 2b.Finally, as the rotor rotates to position 3, the receiving coil induced voltage becomes zero again, completing one whole cycle of the measurement.The measurement principle of the proposed sensor can be considered as "elec magnetic-electric".When the excitation coil is fed with a high-frequency alternating rent, an alternating magnetic field ( e B ) is generated, and according to the Faraday law electromagnetic induction, the voltage is induced in the receiving coil.The rotor in alternating magnetic field induces a secondary magnetic field ( r B ) of the same freque and the eddy current field here also induces a voltage in the receiving coil.Since the r is distributed at intervals, the coupling area between the rotor and the receiving co related to the rotor rotation angle, and the variation of the coupling area changes the duced voltage in the receiving coil.Therefore, the receiving coil induced voltage is du the superposition of the excitation magnetic field According to the law of electromagnetic induction as given in Equation (1), the duced voltage is related to the magnetic induction intensity and the magnetic flux thro a closed loop area.In this work, the receiving coil is designed as a differential struct as shown in Figure 2a.Specifically, the geometry of the two adjacent loops of the receiv coil is the same.The induced voltages are equal in magnitude and opposite in directio the excitation magnetic field, so the effect of the excitation coil on the induced voltag the receiving coil is zero.Evidently, the differential structure of the receiver coil elimin the influence of the excitation field on the induced voltage of the receiver coil.It sho be noted that the coupling area between the rotor and the receiver coil is the key fa affecting the sensor-induced voltage.According to Figure 2a, for the sake of illustrat the receiving coil can be regarded as a series connection of the primary coil and the ondary coil.The coupling area between the rotor and the receiving coil changes as follo (a) when the rotor is at position 1, the rotor is in the middle of the loop consisting of primary and secondary coils, and the direction of the voltage induced by the primary secondary coils is demonstrated by the arrow in the figure.The induced voltages of si as the rotor rotates, the coupling area between the rotor and sinusoidal loop 1 gradu decreases, while the coupling area between the rotor and sinusoidal loop 2 graduall creases.Likewise, the magnetic flux through sinusoidal loop 2 gradually increases, the induced voltage also gradually increases.When the rotor rotates to position 2, th tor and sinusoidal loop 2 are fully coupled, and at this point, the magnetic flux thro the sinusoidal loop 2 reaches the maximum.Meanwhile, the coupling area between rotor and the sinusoidal loop 1 is 0, hence the magnetic flux within sinusoidal loop also zero.Here, the induced voltage of the receiving coil is at the maximum, as show Smax in Figure 2b; (c) as the rotor rotates further, the coupling area between the rotor sinusoidal loop 2 gradually decreases, whereas the coupling area between the rotor sinusoidal loop 3 gradually increases, and the induced voltage slowly decreases.W the rotor is again in the middle of the loop consisting of the primary and secondary c the induced voltage in the receiving coil is zero; (d) with more rotation of the rotor receiving coil induced voltage reaches a negative maximum, where the rotor is fully pled to sinusoidal loop 3, corresponding to position −Smax in Figure 2b.Finally, as the r rotates to position 3, the receiving coil induced voltage becomes zero again, comple one whole cycle of the measurement.According to the above analysis, the magnitude of the induced voltage of the rec ing coil is related to the magnetic induction intensity and the coupling area, S, betw the rotor and the receiving coil.The coupling area between the rotor and the receiving is elaborated in Figure 3, and accordingly, the variation of the rotor-receiver coil coup area can be expressed as Equation (3).
where rd is the stator radius; NR is the number of rotor blades; and hj is the receiving height.According to the above analysis, the magnitude of the induced voltage of the receiving coil is related to the magnetic induction intensity and the coupling area, S, between the rotor and the receiving coil.The coupling area between the rotor and the receiving coil is elaborated in Figure 3, and accordingly, the variation of the rotor-receiver coil coupling area can be expressed as Equation (3).
where r d is the stator radius; N R is the number of rotor blades; and h j is the receiving coil height.

Simulation Analysis
To analyze the magnetic induction intensity, eddy current, and induced voltage of the vertical non-contact induction angle sensor, Solidworks was used to build the sensor model based on the structural parameters given in Table 1, which was then imported into COMSOL for the simulations.In the simulations, the material for the excitation coil, the receiver coil, and the rotor was set to copper, and the grid division results are shown in Figure 4.The input voltage for the excitation coil is shown in Figure 5a, i.e., U = 5 × sin(2π × 10,000 × t) V, while the voltage induced in the receiving coil is shown in Figure 5b.Here, the induced voltage was generated using the rotor eddy current field, and according to the law of flutters, the direction of the induced voltage is opposite to that of the excitation voltage.The induced voltage was generated from the excitation coil to the rotor and then to the receiving coil.At the initial stages, the induced voltage had not yet reached the

Simulation Analysis
To analyze the magnetic induction intensity, eddy current, and induced voltage of the vertical non-contact induction angle sensor, Solidworks was used to build the sensor model based on the structural parameters given in Table 1, which was then imported into COMSOL for the simulations.In the simulations, the material for the excitation coil, the receiver coil, and the rotor was set to copper, and the grid division results are shown in Figure 4.

Simulation Analysis
To analyze the magnetic induction intensity, eddy current, and induced voltage of the vertical non-contact induction angle sensor, Solidworks was used to build the sensor model based on the structural parameters given in Table 1, which was then imported into COMSOL for the simulations.In the simulations, the material for the excitation coil, the receiver coil, and the rotor was set to copper, and the grid division results are shown in Figure 4.The input voltage for the excitation coil is shown in Figure 5a, i.e., U = 5 × sin(2π × 10,000 × t) V, while the voltage induced in the receiving coil is shown in Figure 5b.Here the induced voltage was generated using the rotor eddy current field, and according to the law of flutters, the direction of the induced voltage is opposite to that of the excitation voltage.The induced voltage was generated from the excitation coil to the rotor and then to the receiving coil.At the initial stages, the induced voltage had not yet reached the steady state and had a small amplitude.After the first cycle, the induced voltage reached The input voltage for the excitation coil is shown in Figure 5a, i.e., U = 5 × sin(2π × 10,000 × t) V, while the voltage induced in the receiving coil is shown in Figure 5b.Here, the induced voltage was generated using the rotor eddy current field, and according to the law of flutters, the direction of the induced voltage is opposite to that of the excitation voltage.The induced voltage was generated from the excitation coil to the rotor and then to the receiving coil.At the initial stages, the induced voltage had not yet reached the steady state and had a small amplitude.After the first cycle, the induced voltage reached the steady state with an amplitude of 34.5 mV and a frequency of 10 kHz, which was the same as the excitation frequency.The magnetic induction intensity of the rotor is shown in Figure 5c, where it can be seen that the magnetic induction intensity was large inside the excitation coil; its maximum is 12.6 × 10 −3 T. The rotor vortex distribution is shown in Figure 5d, and the maximum rotor vortex density was 1.41 × 10 −7 A/m 2 .From Figure 5d, it can be seen that the vortex was mainly distributed at the bottom and the sides of the rotor.
Sensors 2024, 24, x FOR PEER REVIEW 6 of 19 the steady state with an amplitude of 34.5 mV and a frequency of 10 kHz, which was the same as the excitation frequency.The magnetic induction intensity of the rotor is shown in Figure 5c, where it can be seen that the magnetic induction intensity was large inside the excitation coil; its maximum is 12.6 × 10 −3 T. The rotor vortex distribution is shown in Figure 5d, and the maximum rotor vortex density was 1.41 × 10 −7 A/m 2 .From Figure 5d, it can be seen that the vortex was mainly distributed at the bottom and the sides of the rotor.The receiving winding of the sensor consisted of three coils, which were staggered in a spatial layout with a 10° offset.The induced voltages in the three receiving coils represent sinusoidal curves of the same frequency and amplitude, but with a phase difference of 10°, i.e., a three-phase signal, as shown in Figure 6.This three-phase signal can be mathematically expressed using Equation (4).The receiving winding of the sensor consisted of three coils, which were staggered in a spatial layout with a 10 • offset.The induced voltages in the three receiving coils represent sinusoidal curves of the same frequency and amplitude, but with a phase difference of 10 • , i.e., a three-phase signal, as shown in Figure 6.This three-phase signal can be mathematically expressed using Equation (4).
where, U 1 , U 2 , and U 3 are the induced voltages of the receiving coil; U j is the amplitude of the induced voltage; and ω e is the angular frequency of the induced voltage.
where, U1, U2, and U3 are the induced voltages of the receiving coil; Uj is the amplitu the induced voltage; and ωe is the angular frequency of the induced voltage.

Non-Linearity Optimization
The structural parameters of the vertical non-contact induction angle sensor the magnetic induction intensity and the eddy current distribution, which in turn the waveform of the induced voltage in the receiving coil.Notably, unreasonable s tural parameters lead to distortion of the induced sine waveform, resulting in large linearity.Limiting the magnitude of non-linearity is critical to enable a satisfactory formance of the sensor; therefore, the structural parameters of the sensor also need be optimized for reduced non-linearity.Before optimization, the three induced vol needed to be linearized.According to Equation ( 5), the CLARK transformation was formed on the three induced voltages to transform the three-phase signal into a two-p signal.The signals before and after the transformation are shown in Figure 7a.As sh in Figure 7a, the signal was converted from three-phase to two-phase with the same s amplitude.Next, the transformed two-phase signal was linearized by Equation (6 linearized signal is shown in Figure 7b.
where C1 and C2 are the two-phase quadrature winding voltages; and L is the linea signal.

Non-Linearity Optimization
The structural parameters of the vertical non-contact induction angle sensor affect the magnetic induction intensity and the eddy current distribution, which in turn affect the waveform of the induced voltage in the receiving coil.Notably, unreasonable structural parameters lead to distortion of the induced sine waveform, resulting in large non-linearity.Limiting the magnitude of non-linearity is critical to enable a satisfactory performance of the sensor; therefore, the structural parameters of the sensor also needed to be optimized for reduced non-linearity.Before optimization, the three induced voltages needed to be linearized.According to Equation ( 5), the CLARK transformation was performed on the three induced voltages to transform the three-phase signal into a two-phase signal.The signals before and after the transformation are shown in Figure 7a.As shown in Figure 7a, the signal was converted from three-phase to two-phase with the same signal amplitude.Next, the transformed two-phase signal was linearized by Equation ( 6); the linearized signal is shown in Figure 7b.
where C 1 and C 2 are the two-phase quadrature winding voltages; and L is the linearized signal.
The induced voltages of the three receiving coils during one cycle of the rotor rotation were simulated.Next, the induced voltages were linearized and the non-linearity was calculated according to Equation (7).The results are shown in Table 2, demonstrating that the non-linearity between the simulated and theoretical values was 0.383%.
where U i is the ideal voltage; U s is the simulation voltage; and U max is the voltage range.The induced voltages of the three receiving coils during one cycle of the rotor rotation were simulated.Next, the induced voltages were linearized and the non-linearity was calculated according to Equation (7).The results are shown in Table 2, demonstrating that the non-linearity between the simulated and theoretical values was 0.383%.

(
) where i U is the ideal voltage; By analyzing the results in Table 2, it was concluded that the receiving coil induction voltage determined the non-linearity of the vertical non-contact induction angle sensor.Combining Equations ( 1)-(3), the most intuitive parameters characterizing the influence of the excitation magnetic field and the eddy current magnetic field on the induced voltage of the receiving coil include: the number of turns of the excitation coil and the excitation coil wire width (to characterize the excitation magnetic field effect); the rotor blade thickness and height (to characterize the vortex field action); the air gap between the stator and the rotor (to characterize the joint action of two magnetic fields); and the height of the receiving coil (to characterize the effect of the magnetic field lines passing through the area of receiving coil, generated by the two magnetic fields).The air gap between the stator and the rotor was determined by the stator radius and the rotor radius.In this work, the rotor radius was fixed at 14 mm; therefore, the change in the stator radius mainly characterized the air gap between the stator and the rotor.By analyzing the results in Table 2, it was concluded that the receiving coil induction voltage determined the non-linearity of the vertical non-contact induction angle sensor.Combining Equations ( 1)-( 3), the most intuitive parameters characterizing the influence of the excitation magnetic field and the eddy current magnetic field on the induced voltage of the receiving coil include: the number of turns of the excitation coil and the excitation coil wire width (to characterize the excitation magnetic field effect); the rotor blade thickness and height (to characterize the vortex field action); the air gap between the stator and the rotor (to characterize the joint action of two magnetic fields); and the height of the receiving coil (to characterize the effect of the magnetic field lines passing through the area of receiving coil, generated by the two magnetic fields).The air gap between the stator and the rotor was determined by the stator radius and the rotor radius.In this work, the rotor radius was fixed at 14 mm; therefore, the change in the stator radius mainly characterized the air gap between the stator and the rotor.
Simultaneous multi-objective optimization often results in wasted resources and low efficiency; therefore, the key factors affecting the sensor non-linearity index need to be appropriately identified before the actual optimization.The Plackett-Burman test is a test method for analytical analysis that is used to estimate the effect of a factor as accurately as possible with a minimum number of trials.This method is suitable for quickly and efficiently screening the most important factors among all the influencing factors, which can be then further investigated.The experimental process here was to code multiple factors at high and low levels, analyze the inter-subjective effects and significance levels of the factors based on the results of Plackett-Burman test, and then filter out the factors with the highest effect on the test results.This essentially reduced the number of factors involved in the optimization.The non-linearity of the vertical non-contact induction angle sensor was first evaluated in terms of the number of turns of the excitation coil, the stator radius, the excitation coil line width, the rotor blade thickness, the rotor length, and the receiver coil height, and then the main affecting factors were selected.Specifically, the above six factors were coded and the lowest (−1) and highest (1) levels were taken for each factor for 12 sets of tests.The experimental factor codes and the high and low levels taken are shown in Table 3. Progressively, the Plackett-Burman experimental design scheme is shown in Table 4, while Table 5 shows the evaluation table for the effect of each factor.Next, the data were subjected to a multi-distance regression analysis and the optimal equation for the non-linearity was obtained, as shown in Equation (8).
No. By analyzing the evaluation table for the effect of each factor (Table 5), the key factors affecting the sensor non-linearity were sorted according to the magnitude of their effect, as follows: X 2 > X 4 > X 1 > X 6 > X 5 > X 3 .Among them, the effects of the stator radius and the rotor blade thickness reached an extremely significant level (p < 0.001), while the effects of the excitation coil turns, the excitation coil wire width, the rotor length, and the receiver coil height were not very significant.Therefore, two key factors, the stator radius and the rotor blade thickness, were selected for inclusion in the subsequent optimization.
Typically, the factors screened using the Plackett-Burman test are used in response surface analysis studies.The optimal parameters obtained using the response surface optimization method depend on the degree of fit of the regression equation, where different regression equations yield different optimal values.In contrast, this paper used the PSO algorithm to optimize the screened factors for reduced non-linearity, which is a mainstream optimization algorithm with fast convergence and only a few setup parameters for the sensor.
In the PSO algorithm, the PSO search space is an n-dimensional space, where the particle population consists of N particles and each particle in the population is initialized with a randomized position ( → x i ) and velocity ( i,n ).At time, t, the position, → x i,t , can be considered as a set of coordinates of a point in the n-dimensional space.The particles fly through the virtual space looking for the candidate solutions and attracting the surrounding particles to a location that produces the best results.Moreover, i,n ) is the individual best position of the i-th particle at moment, t, and n ) is the global best position of the entire particle population at moment, t.The velocity update formula of a particle in the PSO algorithm is given in Equation ( 9).Meanwhile, the maximum distance traveled by a particle in one iteration cycle was determined using the velocity, as expressed in Equation (10).The position and the velocity of each particle were updated once for each iteration of the algorithm.It should be noted that a larger value of inertia weight, ω, is beneficial for improving the global search capability of the algorithm, while a smaller value improves the accuracy of local search.Therefore, the inertia weight, ω, used in this work decreased linearly with the number of iterations, as shown in Equation (11), which ensured the global search capability of the algorithm while avoiding falling into the local optimal solutions. x where c 1-PSO and c 2-PSO are the acceleration coefficients; r 1-PSO and r 2-PSO are the uniformly distributed random numbers on the interval [0, 1]; i = 1, 2, . . .N PSO , N PSO is the particle swarm size; j = 1, 2,. ..n, n is the number of spatial dimensions; ω is the value of the inertia weights; ω max is the maximum value of the inertia weight; ω min is the minimum value of the inertia weights; and t max is the maximum number of iteration steps.
To minimize the non-linearity of the proposed sensor, the first step was to establish the objective function of the optimization problem.The objective function needs to satisfy one or more conditions to obtain the optimal solution by adjusting the system parameters.The non-linearity optimization of the vertical non-contact induction angle sensors is a bi-objective optimization problem, and there is incommensurability between the two objectives, i.e., there exists a conflicting relationship between the two objectives.Specifically, when one of the objective functions reaches the optimality, the other objective function may become worse.On the other hand, this phenomenon is not present in single-objective opti-mization problems; therefore, the bi-objective optimization problem needs to be converted into a single-objective optimization problem.Accordingly, in this paper, the objective function for optimization was reconstructed, where the weight coefficients were introduced to reflect the importance of the two objectives in the overall objective function, as shown in Equations ( 12)- (14).
where L(r) is the objective function of the rotor thickness parameter; L(d) is the objective function of stator radius parameter; L min is the non-linear degree of reconstruction; w r is the rotor thickness weighting factor; w d is the stator radius weighting factor; U ir is the ideal voltage for rotor thickness; U sr is the simulation voltage for rotor thickness; U id is the ideal voltage for stator radius; U sd is the simulation voltage for stator radius; U rmax is the maximum voltage range of rotor thickness; U dmax is the maximum voltage range of stator radius.w r and w d are the weight coefficients introduced to transform the bi-objective optimization problem into a single-objective extremum problem.The weight coefficients were normalized according to the effect values obtained using the Plackett-Burman test.In particular, the stator radius effect value (E d ) of 0.053 and the rotor blade thickness effect value (E r ) of 0.024 in Table 5 were normalized to within the range [0, 1].Then, the weighting coefficients for the stator radius and the rotor blade thickness were expressed as: w r = E r E d + E r = 0.024 0.053 + 0.024 = 0.311 (16) Additionally, the search ranges for the rotor thickness and the stator radius are shown in Table 6, and the parameters of the PSO algorithm were set as follows: the number of particles was set to 30; the learning factor c 1 = c 2 = 1.4945; the maximum velocity was limited to 0.05; the minimum velocity was limited to −0.05; the number of iterations was set to 30; the maximum inertia weight value ω max = 0.9; the minimum inertia weight value ω min = 0.4.The relationship between the non-linearity and the number of iterations is illustrated in Figure 8, where it can be seen that the non-linearity stabilizes at 0.366% when the number of iterations exceeds 18. Progressively, as shown in Figure 9, the particle population aggregated at a non-linearity of 0.366%, where the corresponding rotor blade thickness was 0.52 mm and the stator radius was 15.1 mm.
By optimizing the rotor blade thickness and the stator radius, the optimal design parameters of the sensor were obtained, which are shown in Table 7.The simulation yielded a non-linearity of 0.366% for the vertical non-contact angle sensor, as shown in Figure 10.
The relationship between the non-linearity and the number of iterations is illustrated in Figure 8, where it can be seen that the non-linearity stabilizes at 0.366% when the number of iterations exceeds 18. Progressively, as shown in Figure 9, the particle population aggregated at a non-linearity of 0.366%, where the corresponding rotor blade thickness was 0.52 mm and the stator radius was 15.1 mm.By optimizing the rotor blade thickness and the stator radius, the optimal design parameters of the sensor were obtained, which are shown in Table 7.The simulation yielded a non-linearity of 0.366% for the vertical non-contact angle sensor, as shown in Figure 10.The relationship between the non-linearity and the number of iterations is illustrated in Figure 8, where it can be seen that the non-linearity stabilizes at 0.366% when the number of iterations exceeds 18. Progressively, as shown in Figure 9, the particle population aggregated at a non-linearity of 0.366%, where the corresponding rotor blade thickness was 0.52 mm and the stator radius was 15.1 mm.By optimizing the rotor blade thickness and the stator radius, the optimal design parameters of the sensor were obtained, which are shown in Table 7.The simulation yielded a non-linearity of 0.366% for the vertical non-contact angle sensor, as shown in Figure 10.

Fabrication and Testing of the Sensor
In this work, a customized sensor performance test bench was used as an experimental platform to evaluate the performance of the proposed vertical non-contact induction angle sensor, as shown in Figure 11a.The main components of the experimental equipment included a servo motor system, an angle encoder, a vertical non-contact induc-

Fabrication and Testing of the Sensor
In this work, a customized sensor performance test bench was used as an experimental platform to evaluate the performance of the proposed vertical non-contact induction angle sensor, as shown in Figure 11a.The main components of the experimental equipment included a servo motor system, an angle encoder, a vertical non-contact induction angle sensor, and a data acquisition system.The vertical non-contact angle sensor was mounted on a rotating shaft, as shown in Figure 11b, and Figure 11c shows the experimental data acquisition module.The servo motor in the system was controlled to rotate at 5 r/min in the counterclockwise and clockwise directions, respectively, and the angle measured by the vertical noncontact angle sensor was collected to calculate the non-linearity of the sensor, as shown in Figure 12.The experimental results were: non-linearity was 0.648% for the clockwise sensors; non-linearity of the counterclockwise direction sensor was 0.646%.The servo motor in the system was controlled to rotate at 5 r/min in the counterclockwise and clockwise directions, respectively, and the angle measured by the vertical non-contact angle sensor was collected to calculate the non-linearity of the sensor, as shown in Figure 12.The experimental results were: non-linearity was 0.648% for the clockwise sensors; non-linearity of the counterclockwise direction sensor was 0.646%.
Furthermore, to test the non-linearity of the proposed angle sensor at different rotational speeds, the servo motor was controlled to rotate at 15 r/min, 30 r/min, and 60 r/min in the clockwise and counterclockwise directions, and the corresponding results are shown in Figure 13.Moreover, the maximum values of the non-linearities are given in Table 8.
From Table 8, it can be seen that the non-linearity of the vertical non-contact angle sensor was 0.696% at the maximum for the different rotational speeds.It can be concluded that the rotational speed had little effect on the non-linearity of the sensor.The servo motor in the system was controlled to rotate at 5 r/min in the counterclockwise and clockwise directions, respectively, and the angle measured by the vertical noncontact angle sensor was collected to calculate the non-linearity of the sensor, as shown in Figure 12.The experimental results were: non-linearity was 0.648% for the clockwise sensors; non-linearity of the counterclockwise direction sensor was 0.646%.Furthermore, to test the non-linearity of the proposed angle sensor at different rotational speeds, the servo motor was controlled to rotate at 15 r/min, 30 r/min, and 60 r/min in the clockwise and counterclockwise directions, and the corresponding results are shown in Figure 13.Moreover, the maximum values of the non-linearities are given in Table 8.From Table 8, it can be seen that the non-linearity of the vertical non-contact angle

Sine Fit
The vertical non-contact angle sensor designed in this paper used a receiving coil with a sinusoidal structure, which provided a better sinusoidal fit than the existing diamond-shaped receiving coils.To elaborate further, suppose the diagonal of the diamondshaped receiving coil is 1 mm, the height of the sinusoidal receiving coil is 1 mm, and the rotor blade moves at a speed of 0.1 mm/s, as shown in Figure 14a.The variation curves of the coupling area between the rotor and the diamond-shaped coil and sinusoidal coil, respectively, are shown in Figure 14b.

Sine Fit
The vertical non-contact angle sensor designed in this paper used a receiving coil with a sinusoidal structure, which provided a better sinusoidal fit than the existing diamondshaped receiving coils.To elaborate further, suppose the diagonal of the diamond-shaped receiving coil is 1 mm, the height of the sinusoidal receiving coil is 1 mm, and the rotor blade moves at a speed of 0.1 mm/s, as shown in Figure 14a.The variation curves of the coupling area between the rotor and the diamond-shaped coil and sinusoidal coil, respectively, are shown in Figure 14b.
The curves for the variation in the rotor and diamond coil coupling area and the variation in the rotor and sinusoidal coil coupling area were fitted as the sinusoidal functions, where the degree of fit was measured using the goodness-of-fit function, R-squared, which can be calculated using Equation (17).
where SST is the square sum of the difference between the original data and the average value; and SSE is the sum square of errors of the fitting data and the original data.The curves for the variation in the rotor and diamond coil coupling area and the variation in the rotor and sinusoidal coil coupling area were fitted as the sinusoidal functions, where the degree of fit was measured using the goodness-of-fit function, R-squared, which can be calculated using Equation (17).
where SST is the square sum of the difference between the original data and the average value; and SSE is the sum square of errors of the fitting data and the original data.
The R-squared interval ranges from [0, 1], and the closer its value is to 1, the better the fit.The fit of the diamond-shaped coil was calculated as 0.998, while the fit of the sinusoidal coil was 1. Experiments have shown that the sinusoidal coil has a better sinusoidal fit.From the perspective of manufacturing processes, the diamond-shaped coils are easier to implement and save manufacturing costs.However, the sinusoidal coils offer the following advantages: (1) sinusoidal coils can reduce electromagnetic interference and crosstalk in high-frequency circuits, thus improving electromagnetic compatibility; (2) for highspeed signal lines and specific protocol signal lines, sinusoidal coils can reduce signal waveform distortion and delay distortion, thereby enhancing signal integrity; (3) sinusoidal coils can avoid the "sharp corners" introduced by right-angle traces, thus reducing the complexity of the PCB layout and the volume and weight of the circuit board.

Magnetic Induction Intensity
Notably, the vertical non-contact angle sensor designed in this paper used the excitation coil with a solenoidal structure to generate a large magnetic induction intensity.Table 9 compares the magnetic induction intensity of the planar type angle sensor and the sensor designed in this paper.Evident from Table 9, the sensor designed in this paper produced a greater magnetic induction intensity with a smaller excitation source; therefore, the sensor designed in this paper had a stronger anti-electromagnetic interference capability.The R-squared interval ranges from [0, 1], and the closer its value is to 1, the better the fit.The fit of the diamond-shaped coil was calculated as 0.998, while the fit of the sinusoidal coil was 1. Experiments have shown that the sinusoidal coil has a better sinusoidal fit.From the perspective of manufacturing processes, the diamond-shaped coils are easier to implement and save manufacturing costs.However, the sinusoidal coils offer the following advantages: (1) sinusoidal coils can reduce electromagnetic interference and crosstalk in high-frequency circuits, thus improving electromagnetic compatibility; (2) for high-speed signal lines and specific protocol signal lines, sinusoidal coils can reduce signal waveform distortion and delay distortion, thereby enhancing signal integrity; (3) sinusoidal coils can avoid the "sharp corners" introduced by right-angle traces, thus reducing the complexity of the PCB layout and the volume and weight of the circuit board.

Magnetic Induction Intensity
Notably, the vertical non-contact angle sensor designed in this paper used the excitation coil with a solenoidal structure to generate a large magnetic induction intensity.Table 9 compares the magnetic induction intensity of the planar type angle sensor and the sensor designed in this paper.Evident from Table 9, the sensor designed in this paper produced a greater magnetic induction intensity with a smaller excitation source; therefore, the sensor designed in this paper had a stronger anti-electromagnetic interference capability.

Volume
The proposed sensor was compared with the HELLA Torque Angle sensor.The size of the proposed vertical non-contact angle sensor was 35 × 35 × 40 mm 3 , whereas the HELLA Torque Angle Sensor had the dimensions of 70 × 70 × 15 mm 3 .Evidently, the sensor designed in this paper offered a 50% reduction in radial size and a 32.8% reduction in overall size, compared with the HELLA torque angle sensor shown in Figure 15.This essentially makes the vertical non-contact angle sensor more suitable for installation in spaces such as vehicles or robots.

Conclusions
In this paper, a new vertical non-contact induction angle sensor was designed, consisting of a solenoidal excitation coil and three receiver coils that form the stator of the sensor.In addition, a ring-shaped metal sheet with six rectangular blades was utilized to form the rotor of the sensor.Using the principle of differential measurement, the receiver coil was designed to have a sinusoidal structure, and the induced voltages generated by the excitation source in the two adjacent loops of the receiver coil cancelled out each other, which decoupled the excitation field from the eddy current field.In addition, the rotor angle was measured using the change in the coupling area between the rotor and the receiver coil, which improved the sinusoidal fit of the induced voltage.To study the sensor theoretically, a finite element analysis was performed using the COMSOL 6.1 software, where the induced voltage was a sinusoidal signal with frequency of 10 kHz and an amplitude of 34.5 mV.To analyze the non-linearity of the sensor, an algorithm using the CLARK transform method with an inverse tangent function was proposed that linearized the three-phase signal, which in turn reduced the non-linearity compared with the case when the segmented lines were used.The simulation results showed that the non-linearity of proposed sensor was only 0.383%.To reduce the sensor non-linearity, the stator radius and the rotor blade thickness were identified as the key factors affecting the sensor nonlinearity, using the Plackett-Burman test.Then, based on the effect values of the stator radius and the rotor blade thickness, an optimization objective function was constructed, following which the PSO algorithm with decreasing inertia weights was used to optimize the sensor non-linearity.The optimization result showed that the sensor non-linearity was minimum at 0.366%, when the rotor thickness was 0.52 mm and the stator radius was 15.1 mm.Besides, the experimental results showed that the non-linearity was 0.648% for the clockwise rotation and 0.646% for the counterclockwise rotation.The proposed vertical non-contact induction angle sensors offer the following advantages: (a) comparative studies with commercially available HELLA sensors show that the sensors in this paper were smaller in size; (b) the magnetic induction intensity of the vertical induction sensor was greater than that of the planar type induction angle sensor; (c) compared with the optical sensors and the HALL sensors, the vertical angle sensors had negligible sensitivity to moisture, dust, oil, etc.The above advantages make the proposed sensor a promising device for automotive and robotics applications.

Conclusions
In this paper, a new vertical non-contact induction angle sensor was designed, consisting of a solenoidal excitation coil and three receiver coils that form the stator of the sensor.In addition, a ring-shaped metal sheet with six rectangular blades was utilized to form the rotor of the sensor.Using the principle of differential measurement, the receiver coil was designed to have a sinusoidal structure, and the induced voltages generated by the excitation source in the two adjacent loops of the receiver coil cancelled out each other, which decoupled the excitation field from the eddy current field.In addition, the rotor angle was measured using the change in the coupling area between the rotor and the receiver coil, which improved the sinusoidal fit of the induced voltage.To study the sensor theoretically, a finite element analysis was performed using the COMSOL 6.1 software, where the induced voltage was a sinusoidal signal with frequency of 10 kHz and an amplitude of 34.5 mV.To analyze the non-linearity of the sensor, an algorithm using the CLARK transform method with an inverse tangent function was proposed that linearized the three-phase signal, which in turn reduced the non-linearity compared with the case when the segmented lines were used.The simulation results showed that the non-linearity of proposed sensor was only 0.383%.To reduce the sensor non-linearity, the stator radius and the rotor blade thickness were identified as the key factors affecting the sensor nonlinearity, using the Plackett-Burman test.Then, based on the effect values of the stator radius and the rotor blade thickness, an optimization objective function was constructed, following which the PSO algorithm with decreasing inertia weights was used to optimize the sensor non-linearity.The optimization result showed that the sensor non-linearity was minimum at 0.366%, when the rotor thickness was 0.52 mm and the stator radius was 15.1 mm.Besides, the experimental results showed that the non-linearity was 0.648% for the clockwise rotation and 0.646% for the counterclockwise rotation.The proposed vertical non-contact induction angle sensors offer the following advantages: (a) comparative studies with commercially available HELLA sensors show that the sensors in this paper were smaller in size; (b) the magnetic induction intensity of the vertical induction sensor was greater than that of the planar type induction angle sensor; (c) compared with the

→rl
is the distance between the excitation source and the induced conductor; is the current element; and µ 0 is the permeability of vacuum,

Figure 1 .
Figure 1.Structure of the vertical non-contact angle sensor.

eB
and the eddy current field B shown in Equation(1).Next, the magnetic induction intensity can be calculated as sho in Equation (2).
the number of turns; S is the area; I is the excitation current; r is the dista between the excitation source and the induced conductor;

Figure 1 .
Figure 1.Structure of the vertical non-contact angle sensor.

Figure 2 .
Figure 2. (a) Coupling between the rotor and receiving coil areas; (b) Induced voltage of the re ing coil corresponding to the different rotor positions.

Figure 2 .
Figure 2. (a) Coupling between the rotor and receiving coil areas; (b) Induced voltage of the receiving coil corresponding to the different rotor positions.

Figure 3 .
Figure 3. Schematic diagram of the coupling area between the rotor and the receiving coil.

Figure 4 .
Figure 4. Meshing of the vertical non-contact angle sensor.

Figure 3 .
Figure 3. Schematic diagram of the coupling area between the rotor and the receiving coil.

Table 1 .Figure 3 .
Figure 3. Schematic diagram of the coupling area between the rotor and the receiving coil.

Figure 4 .
Figure 4. Meshing of the vertical non-contact angle sensor.

Figure 4 .
Figure 4. Meshing of the vertical non-contact angle sensor.

Figure 6 .
Figure 6.Receiving coil winding and three-phase signal.

Figure 6 .
Figure 6.Receiving coil winding and three-phase signal.

Figure 7 .
Figure 7. (a) A three-phase signal is converted to a two-phase signal; (b) Two-phase signal linearization.

sU
is the simulation voltage; and max U is the voltage range.

Figure 7 .
Figure 7. (a) A three-phase signal is converted to a two-phase signal; (b) Two-phase signal linearization.

Figure 8 .
Figure 8. Relation between the non-linearity and the number of iterations.

Figure 8 .
Figure 8. Relation between the non-linearity and the number of iterations.

Figure 8 .
Figure 8. Relation between the non-linearity and the number of iterations.

Figure 10 .
Figure 10.Non-linearity of the vertical angle sensor under optimal structural parameters.

Figure 10 .
Figure 10.Non-linearity of the vertical angle sensor under optimal structural parameters.

Figure 14 .
Figure 14.(a) Schematic diagram of the coupling process between the rotor and the receiving coil; (b) Change of coupling area between the rotor and the receiving coils.

Figure 14 .
Figure 14.(a) Schematic diagram of the coupling process between the rotor and the receiving coil; (b) Change of coupling area between the rotor and the receiving coils.

Sensors 2024 , 19 Figure 15 .
Figure 15.Comparison of the HELLA angle sensor and the vertical non-contact angle sensor.1. HELLA angle sensor; 2. Vertical non-contact angle sensor.Precautions for use and safety: (1) Please operate within a temperature range of −40 to 80 °C; (2) Use the power supply voltage and load within specified ranges and specifications; (3) Keep the sensor away from other electrical equipment as far as possible to minimize the interference signals.

Figure 15 .
Figure 15.Comparison of the HELLA angle sensor and the vertical non-contact angle sensor.1. HELLA angle sensor; 2. Vertical non-contact angle sensor.Precautions for use and safety: (1) Please operate within a temperature range of −40 to 80 • C; (2) Use the power supply voltage and load within specified ranges and specifications; (3) Keep the sensor away from other electrical equipment as far as possible to minimize the interference signals.

Table 2 .
Non-linearity of vertical non-contact angle sensor.

Table 2 .
Non-linearity of vertical non-contact angle sensor.

Table 3 .
Non-linearity of vertical non-contact angle sensor.

Table 6 .
Particle search range.

Table 8 .
Non-linearity of the sensors at different speeds.

Table 8 .
Non-linearity of the sensors at different speeds.

Table 9 .
Comparison of the magnetic induction intensity of the planar type angle sensor and the sensor designed in this paper.

Table 9 .
Comparison of the magnetic induction intensity of the planar type angle sensor and the sensor designed in this paper.