Next Article in Journal
A Method for Extracting Debye Parameters as a Tool for Monitoring Watered and Contaminated Soils
Previous Article in Journal
Design Space Exploration of Clustered Sparsely Connected MPSoC Platforms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Effective Method to Model and Simulate the Behavior of Inductive Angle Encoders

1
Robert Bosch GmbH, Robert-Bosch-Allee 1, 74232 Abstatt, Germany
2
Faculty of Electrical and Computer Engineering, Technische Universität Dresden, Helmholtzstraße 10, 01069 Dresden, Germany
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(20), 7804; https://doi.org/10.3390/s22207804
Submission received: 15 September 2022 / Revised: 27 September 2022 / Accepted: 3 October 2022 / Published: 14 October 2022
(This article belongs to the Section Electronic Sensors)

Abstract

:
Many angle or position sensors rank among the inductive encoders, although their sensing principle is different. The sensor design investigated in this paper is based on coupled coils, whereas the information about the position angle is modulated on the induced voltage, measured at the receiving coils. Unfortunately, no closed solution for most of the physical quantities exists, since this principle is based on eddy currents, which are rather complex to calculate for the given geometry. Consequently, the common way is to calculate the sensor quantities by a 3D finite-element (FE) simulation. However, this leads in most cases to a high time and computational effort. To overcome the limitations with respect to computational resources, a novel method is presented to reduce simulation effort and calculate regression models, which can even replace simulations. In the following investigations, D-optimal designs are used—a subdomain in the field of statistical design of experiments—and combined with a numerical implementation of Faraday’s law, in order to calculate the induced voltages afterwards from simulated magnetic field data. With this method, the sensor signals can be calculated for multiple angle positions from one simulated position by shifting the integration boundaries. Hence, simulation time is significantly reduced for a full period. The regression models obtained by this method, can predict the Tx-coil inductance, induced Rx-voltage amplitude and angular error in dependency of geometric design parameters.

1. Introduction

This work contributes to inductive angle encoders. Many angle encoders—independent of their operating principle—share one typical characteristic: The angle is not measured directly, but rather encoded in two orthogonal signals—commonly known as sine and cosine channels in quadrature—that form a vector in the complex plane. The actual angle is then calculated with the (unwrapped) arctangent function, essentially converting Cartesian coordinates to a polar angle.
The following investigations focus on a printed circuit board-based (PCB) sensor design which contains a planar coil system—comprising a transmitter coil (Tx-coil) and at least two receiver coils (Rx-coils)—a rotating coupling element (target) made out of conductive material and a signal-processing circuitry for measuring, demodulation and amplification of the coil-signals (Figure 1). For this sensor design the position angle can be derived from the induced voltage signals on the Rx-coils. The Tx-coil is arranged in an oscillator circuit in order to generate a time-varying magnetic field. Due to the concentric arrangement of the coils on the PCB, this magnetic field infuses the Rx-coils together with the target, which is placed directly above the coils in a certain distance (air gap).
In absence of a target the coupling between the Tx- and Rx-coils is zero and thus also the induced voltage signals, since the geometry of the Rx-coils leads to a differential signal. If a target is present, the magnetic field of the Tx-coil induces eddy currents, that—pursuant to Lenz’s law—create an opposing field. Accordingly, the magnetic field below the target almost cancels out, but yields an inhomogeneous field distribution in the angular direction above the Rx-coils. A rotation of the target therefore leads to an amplitude-modulated voltage signal on the Rx-coils, that oscillates with the carrier frequency of the Tx-oscillator. Note that within this paper dynamical effects—caused by high rotational speeds—are neglected (quasi-static).
In order to demodulate and amplify the induced voltage signals on the Rx-coils, an output circuitry is needed. Commonly, an application-specific integrated circuit (ASIC) is chosen for that use, since it includes most of the necessary circuitry in one electronic part. Based on the application and the utilized signal-processing circuitry or ASIC, the coil system has to fulfill several requirements. For example, the amplitude of the induced voltage signals needs to be in a specific range depending on the amplifier stage. Unfortunately, no closed analytic solution exists for the coil system quantities. This is mostly related to the eddy current distribution, which is rather complex to calculate for the given geometry. Analytical solutions exist mainly for simple geometries as shown by Nagel [1]. The common way to calculate the sensor quantities numerically is a finite-element (FE) simulation. Unfortunately, most simulation tools are commercial software and require high computational resources and time-effort.
Throughout the literature there are many publications about the application of statistical methods—especially considering optimal designs—on modeling of physical quantities of various sensor types. For example, Gianchandani and Crary presented a similar approach for the modeling of a microaccelerometer, by creating also optimal designs [2]. However, in terms of the inductive position sensing principle based on coupled coils, this method was not applied yet. Regarding this sensor principle, Passarotto et al. showed a time-efficient simulation approach which uses a surface integral method [3]. Hoxha et al. applied this method in combination with a non-linear least-squares solver, in order to optimize the sensor accuracy by changing the Rx-coil geometry [4]. However, the current literature does not provide an in-depth overview about the sensor behavior in dependency of geometric parameters. State of the art is the estimation of various sensor quantities by numerical calculation in form of a FE simulation. An alternative to FE calculations is not presented for this sensor design.
The novel method presented in this paper combines model simplification and numerical calculations with statistical methods, aiming at the reduction of simulation effort. Simulation software will still be used for data generation, but the computational resources are reduced to a minimum (Figure 2). At the end, regression models are calculated which are able to predict sensor quantities—depending on various parameters—without any further simulation.

2. Modeling and Simulation Method

The characteristics of inductive sensors which are based on coupled coils are mostly dependent on the geometry and the arrangement of the individual sensor components. Hence, the following studies focus on the influence of geometric design parameters. A generalized simulation model will be used for the investigations, which is dependent on the intermediate Rx-coil radius r Rx , the Rx peak-to-peak distance l Rx , the number of Tx-windings per layer n w , the number of Tx-layers n l , the air gap between coils and target h a and the periodicity of the sensor p. In order to reduce computational effort, the Rx-coils will be excluded from the simulation model. This is valid as the physical point of view will show, since the magnetic field distribution is—in good approximation—not influenced by the Rx-coils. The next step according to Figure 2 is to determine design parameter variations with a high informational output, considering a limited number of simulations. In the field of statistical design of experiments, D-optimal designs are used to minimize the covariance of estimated regression coefficients β . If a design is not D-optimal, more experiments are needed to achieve the same prediction accuracy. Design means in this context a set of experiments—or in this case simulations—with different parameter variations. Once the set of simulations—given by a so-called design matrix  X —is determined, the dataset of response values y can be generated for the following regression analysis. If the sensor quantity of interest is related to the induced voltage on the Rx-coils, a numeric integration of the simulated magnetic field data must follow. The most important benefit of this method is, that even for a full rotation only one target position needs to be simulated. The rotation can be represented by an iterative calculation with shifted integration boundaries. This is valid, due to the almost symmetric field distribution of the simplified coil system. After data are acquired by simulation and following numeric integration for the Rx-coil quantities, linear regression models can be fitted. For the calculation of the regression coefficients, the common least-squares method is applied.

2.1. Equivalent Circuit Model

As shown in Figure 3 the coil system including the target can be modeled as a transformer circuit. For simplification, the oscillator driver for the Tx-coil can be replaced by an AC-voltage source, which matches the definition in most simulation models. Furthermore, only one Rx-coil is considered, since in an ideal case the difference of both coils lies only in a phase-shift. All following considerations refer to an isolated environment without further conductive material—such as a housing—in close proximity to the sensor. At first, Ohm’s law is written down which results in a linear system for the whole network:
Z Tx Z Tx , Tar Z Tx , Rx Z Tar , Tx Z Tar Z Tar , Rx Z Rx , Tx Z Rx , Tar Z Rx · I Tx I Tar I Rx = U Tx U Tar U Rx , R Tx + j ω L Tx j ω M Tx , Tar 0 j ω M Tar , Tx R Tar + j ω L Tar j ω M Tar , Rx 0 j ω M Rx , Tar j ω L Rx · I Tx I Tar 0 = U Tx 0 U Rx .
The ohmic resistance of the Rx-coils is neglected, since it is a small value compared to the targeted very high resistance of the amplifier input. Therefore, the current through the Rx-coils is almost zero and hence neglectable. Note, that also proximity losses are not included inside the impedance matrix, since they are also considered to be comparatively small. Due to the reciprocity theorem follows for the mutual inductances [5,6,7,8]:
M Tx , Tar = M Tar , Tx , M Tar , Rx = M Rx , Tar .
According to the differential Rx-coil geometry it yields:
M Tx , Rx = M Rx , Tx = 0 ,
which means that there is no direct coupling between the Tx- and Rx-coils.

2.2. Tx-Coil Impedance

One important quantity of the coil system is the impedance of the Tx-coil in presence of the target. The imaginary part of this impedance determines—in combination with at least one capacitor—the frequency of the Tx-oscillator circuit. Since the mutual coupling between Tx- and Rx-coil is zero, the Rx-coil loop can be ignored in the network. Equation (1) is solved for the resulting impedance [9]:
Z Tx , Tar = R Tx + ω 2 M Tx , Tar 2 R Tar R Tar 2 + ω 2 L Tar 2 + j ω L Tx ω 2 M Tx , Tar 2 L Tar R Tar 2 + ω 2 L Tar 2 .
Equation (4) shows, that even without considering proximity effects, the resulting ohmic resistance and inductance—which can be measured e.g., with an impedance analyzer—is mainly influenced by the self-inductance of the target and the mutual inductance between the Tx-coil and the target. Accordingly as an example, a change of the air gap between the target and the coil system would affect the ohmic resistance measured at the Tx-coil connections and also the oscillator frequency.

2.3. Rx-Coil Induced Voltage

Since the sensor can be modeled as a transformer, it is obvious to calculate also a transformation ratio between the Tx- and Rx-coils which yields
U Tx U Rx = M Tx , Tar I Tar L Tx I Tx M Rx , Tar I Tar + j R Tx I Tx ω M Rx , Tar I Tar 0 .
The imaginary part of this expression is set to zero, due to the negligibly small numerator compared to the denominator. Following, the induced voltage of the Rx-coils can be written as:
U Rx = M Rx , Tar M Tx , Tar L Tx I Tx I Tar · U Tx .
As a result, the oscillator frequency is not directly included in Equation (6) and for the operating frequency range the self-inductance of the Tx-coil and the mutual inductances can be considered as constants. In good approximation, the induced voltage amplitude on the Rx-coils does not depend on the Tx-oscillator frequency.

2.4. Angular Error

Commonly, encoders are used with ambiguous sensor signals over 360° of mechanical rotation. This ambiguity is caused by the periodicity p of the encoder and leads to a distinction between the mechanical angle ϕ and the electrical angle ϕ p
ϕ p = p ϕ .
For the inductive coupled coils sensor the periodicity p equals the number of target wings. To calculate the encoded angle from a and b, the arctangent function is used. In case of ideal, undistorted signals, it yields:
atan a ( ϕ ) b ( ϕ ) = atan sin ( p ϕ ) cos ( p ϕ ) = atan tan ( p ϕ ) = ϕ p .
Equation (8) is only valid for a restricted range of p ϕ because of quadrant ambiguities and eventual division by zero. In practice, this problem is solved by a modified two-argument arctangent function well known as atan2 ( a , b ) (Figure 4).
Every encoder is influenced by undesired disturbances, which will cause an angular error. The electrical angular error be defined as the difference from the electrical reference angle
Δ ϕ p = atan ( a ( ϕ ) b ( ϕ ) ) p ϕ .
With respect to Equation (7), the mechanical angular error is defined as
Δ ϕ = 1 p Δ ϕ p .
This shows the main benefit of a periodicity p > 1 : The angular error over a mechanical rotation decreases for higher periodicities p. However, the absolute angle can no longer be determined unambiguously without further measures.

3. Parameter Definition and Model Simplification

The first step for building up a general applicable simulation model is the parametrization of the sensor geometry. Since the Rx-coil is the most decisive element, the definition focuses on its dimensions. Due to the concentric arrangement of the coils, the dimensions of the Tx-coil and target following implicitly. Regarding the simplification, the Rx-coils are discarded from the simulation model. This way of simplifying the model is valid, since there is only neglectable coupling between the Tx- and Rx-coils (see Equation (3)). However, with the numerical implementation of Faraday’s law, results for the induced voltages can still be obtained.

3.1. Definition of the Sensor Geometry

Figure 5 shows the parametrization of the sensor geometry. The dimensions of the model are mostly dependent on the Rx-coil radius r Rx and the peak-to-peak distance l Rx . For the Tx-coil and the target follows a implicit definition, since the inner radius of the Tx-coil is defined as
r i , Tx = r Rx + l Rx / 2 + l d , Tx ,
where l d , Tx is a constant distance between the Rx- and Tx-coil. For the inner radius of the target follows
r i , Tar = r Rx l Rx / 2 .
The outer radius of the Tx-coil depends on the number of windings n w and the constant copper-width l w , Cu :
r o , Tx = r i , Tx + ( 2 n w 1 ) · l w , Cu ,
whereby the outer Target radius equals the outer Tx-coil radius
r o , Tar = r o , Tx .
In the following observations the distance between the Rx- and Tx-coil is set to l d , Tx = 0.774 m m and the copper-width to l w , Cu = 0.254 m m .

3.2. Numeric Integration of Field Data

A numerical implementation of Faraday’s law is presented, to calculate the induced voltage of the polar-shaped Rx-coils from simulated B-field data. The simulation software of choice is Ansys Maxwell3D whereby the eddy current solver is used. This solver provides a steady-state solution which neglects wave propagation effects. All currents inside a conductor oscillate with the same user-defined frequency and there is no change in the electric phase angle along a conductor. The defined excitation is sinusoidal, which leads to a complex notation of all quantities. As a result of the simulation, the B-field data will be exported in complex notation on a defined grid in the 3D-space [10,11,12].
The first step is to find a mathematical description of the spanned Rx-surface S, because this is essential for the magnetic flux calculation. As shown in Figure 5, each Rx-coil consists of a positive and negative sine or cosine contour. It is necessary to place each Rx-coil on two PCB-layers in order to be able to use the same layers for both coils. This leads to jumps at every half period of the positive and negative contour (via-positions in Figure 1). Following this, the centerline of the Rx-coil tracks can be defined in cylindrical coordinates for every half period as
C + / : = ( φ , z 1 , z 2 ) R | 0 φ π / p and z 2 > z 1 ,
which yields for the coordinates
C + / = r Rx ± ł Rx / 2 cos ( p φ ) φ z 1 / 2 .
As shown in Figure 6, the spanned surface rotates around the intermediate radius r Rx . A parametrization of the surface is introduced, depending on the curvilinear coordinates ( t , ψ ) . These coordinates defining the position on the unwound surface A (Figure 6b)
A : = ( t , ψ ) R 2 | t min t t max and 0 φ π / p ,
where t denotes the parameter in the interpolated direction between C + and C as well as the second parameter ψ , which equals the angular direction in cylindrical coordinates. The inner and outer limit of A is given by
t min / max = | t Rx | 1 2 | t C + ( φ = 0 ) t C ( φ = 0 ) | = r Rx 2 + z Rx 2 1 2 l Rx 2 + ( z 1 z 2 ) 2 ,
while for the intermediate z-coordinate follows
z Rx = z 1 + z 2 2
and the inclination angle around the intermediate radius, as shown in Figure 6c:
ξ ( ψ ) = atan ( C + z ( ψ ) C z ( ψ ) C + r ( ψ ) C r ( ψ ) ) .
Now the complete parametrization is given with the vector v , which maps every point at ( t , ψ ) to the corresponding point in cylindrical coordinates
v ( t , ψ ) = r ( t , ψ ) φ ( t , ψ ) z ( t , ψ ) = ( t | t Rx | ) cos ( ξ ( ψ ) ) + r Rx ψ ( t | t Rx | ) sin ( ξ ( ψ ) ) + z Rx .
With a suitable parametrization of the spanned surface, the magnetic flux can be calculated
Φ = S B · d S = S ( B · n ) d S = A B ( v ( t , ψ ) ) · v t × v ψ = n t d t d ψ ,
whereby the normal vectors of each surface element are calculated by a cross-product of the partial derivatives of v in t- and ψ -direction. For the numerical implementation, a weighting function W is needed, which sets the normal vectors of all surface elements outside the boundaries to zero [13]
W ( t , ψ ) = 1 , C v ( t , ψ ) C + 0 , otherwise .
Within the definition of this function the notation ( t , ψ ) emphasizes that the relevant points for this case selection are not the grid points, but rather the centerpoints of the surface elements.
After the projection of the simulated B-field values B on the with W weightened normal vectors n , only flux values for surface elements inside the integration boundaries are obtained, as illustrated in Figure 6a. Finally, the induced voltage can be calculated by a summation of the magnetic fluxes of each element and consideration of the angular frequency ω :
U Rx = j ω Φ = j ω B · n .
Instead of rotating the target in the simulation model, the integration boundaries will be rotated in the numerical calculation of Faraday’s law. Accordingly, only one simulation is needed for the calculation of one electrical period of the induced voltages. Since the numerical integration is much faster than the simulation of the full coil system, this leads to a huge benefit regarding computational time. Note, that this approach is valid as long as the field distribution caused by the Tx-coil and target is rather symmetrical in the angular direction.

4. Applied Statistical Methods

Instead of simulating each new coil design mathematical models are used, which will be able to predict various sensor quantities with less computational and time effort. Since this approach is based on linear regression models, simulation data are still once needed for the model calculation. However, this effort will be comparably small, because only a limited number of simulations is required.
In statistics the independent variables are commonly called factors, but in the following the term design parameters will further be used. Within the study, the parameters will be varied in defined steps. Due to the different parameter ranges, the steps will be normalized from −1 to +1. The amount of needed steps per parameter is k + 1 , where k is the highest order of the regression polynomial. A trivial approach is to calculate all parameter combinations. In the field of design of experiments this approach is known as full-factorial design. The goal here will be to minimize the number of experiments with a maximum informational output [14].

4.1. Linear Regression

The classical linear regression model is defined as [14]
y = X · β + ϵ ,
where y is the dependent or response variable
y = y 1 y m ,
X denotes the design matrix of order m × ( n + 1 )
X = 1 x 11 x 1 n 1 x m 1 x m n ,
β corresponds to the regression coefficients
β = β 1 β n
and ϵ represents the disturbance terms
ϵ = ϵ 1 ϵ m .
The response values y will be obtained by simulation in dependence of the respective parameter values in X . Subsequently, the regression coefficients can be calculated by the least-squares method with
ϵ · ϵ T = ϵ i 2 = min ,
which leads to
β = ( X T X ) 1 X T y .
Additionally, the linear regression model is based on the following assumptions [14]:
1.
For the expectation value of the disturbance term follows E ( ϵ ) = 0
2.
The data need to be homoscedastic with Cov ( ϵ ) = σ 2 I
3.
The factors needs to be linearly independent shown by rk ( X ) = n
Note that for higher-order polynomials the regression will still be linear, since some matrix entries x i j will be substituted by a factor of order k.

4.2. D-Optimal Design

In the following, a short overview about the theory behind D-optimal designs will be given. The term design stands for a set of experiments—or in this case simulations—which should be reduced to a minimum. Accordingly, optimization means in this context the minimization of the design. Starting from a candidate set including all combinations of the parameter levels (full-factorial), the experiment should be limited to a number of N runs. As an criterion for optimization, D-optimal experiments minimize the covariance of the regression coefficients which is given by [15]
min ( Σ β ) = min Cov ( β ) = min σ 2 ( X T X ) 1 .
Since the determinant det ( X T X ) is inversely proportional to the covariance in Equation (32), a row-exchange algorithm is used—similar to the method presented by Fedorov (pp. 104–110)—to maximize this determinant [16]
max det ( X T X ) min σ 2 ( X T X ) 1 .
In the end, the algorithm chooses those N experiments of the candidate set, which lead to D-optimality. It should be noted, that the term X T X is known as Fisher information matrix [14].
By means of considering the normal distribution of two random variables, the relationship between the covariance of the regression coefficients and the determinant of the Fisher information matrix can be shown. In this regard, one important realization is the similarity between the multivariate normal distribution and an ellipsoid [17].
Ellipsoids can be mathematically described by a quadric, which is the solution set of a quadratic equation. Generally spoken, a quadric ranks as an algebraic variety. For an ellipsoid it yields
( x q ) T D ( x q ) = 1 ,
where x is the vector of the independent variables, q is the center-point vector and D a positive definite matrix. Additionally, the eigenvalues λ D i of the matrix D are related to the semi-axes a i of the ellipsoid described by Equation (34) as follows [17]
λ D 1 λ D i = 1 / a D 1 2 1 / a D i 2 .
Figure 7 shows a dataset following the distribution
N ( μ , Σ ) = N 0 0 , 1 1.5 1.5 3 .
An ellipsoid corresponding to this bivariate distribution can be described by [17]
( x μ ) T Σ 1 ( x μ ) = 1 .
Under the assumption of a positive definite matrix Σ the terms
det ( Σ ) = 1 det ( Σ 1 )
are identical, which leads—associated with Equation (35)—to
det ( Σ ) = i = 1 λ Σ i = i = 1 a Σ i 2 .
This shows that the determinant of the covariance—or in respect to Equation (32) the determinant det ( X T X ) —is proportional to the size of the describing ellipsoid and therefore a good criterion for optimizing the design. Note that a D-optimal design needs at least n + 1 experiments, where n is the number of regression coefficients.

5. D-Optimal Designs in Application

As mentioned in the introduction, regression models for the Tx-coil inductance, the induced voltage on the Rx-coils and the angular error are presented by generating D-optimal designs. For all models the highest order will be quadratic and only linear interaction terms are considered. With 6 parameters, this approach yields a model containing 28 terms:
y ˜ = β 0 + β 1 r Rx + β 2 l Rx + β 3 n w + β 4 n l + β 5 h a + β 6 p + β 12 r Rx l Rx + β 13 r Rx n w + β 14 r Rx n l + β 15 r Rx h a + β 16 r Rx p + β 23 l Rx n w + β 24 l Rx n l + β 25 l Rx h a + β 26 l Rx p + β 34 n w n l + β 35 n w h a + β 36 n w p + β 45 n l h a + β 46 n l p + β 56 h a p + β 11 r Rx 2 + β 22 l Rx 2 + β 33 n w 2 + β 44 n l 2 + β 55 h a 2 + β 66 p 2 .
Since the highest order of this polynomial is 2, at least three steps per parameter are needed. Following, a full-factorial design would require 3 6 = 729 experiments. In contrast, a D-optimal design requires at least n + 1 = 28 experiments.
A commonly used indicator for comparing different designs—especially with a different amount of experiments—is the D-efficiency which is calculated as [18]
D eff = det ( X T X ) 1 / ( n + 1 ) m .
For an orthogonal design matrix the determinant det ( X T X ) gets equal to m ( n + 1 ) . Hence, the D-efficiency would be 100 % for such a design.

5.1. Tx-Inductance Model

The first regression model gives an estimation for the Tx-coil inductance in the presence of a target. The polynomial will be dependent on the variables listed in Table 1. In order to find a good compromise between simulation time and informational output, the D-efficiency (according to Equation (41)) is compared for different numbers of experiments in Figure 8. A number of N = 32 runs will be chosen, since the D-efficiency is comparatively high with over 47 % and the available computational resources allow to run 16 simulations in parallel.
Table A1 contains the result of the row-exchange algorithm in normalized notation. Each experiment represents one simulation. Hence, the result will be 32 inductance values which will be fitted. According to monotonous increasing residuals in a first fit, the simulated inductance values are logarithmized. This is a common transformation and leads finally to an exponential model [19].
Figure 9 shows the values of the 13 most significant regression coefficients in descending order. Those values showing the effects of the actual terms, since the regression coefficients can also be interpreted as a weighting factor. In this case, the number of windings n w , number of layers n l and intermediate radius r Rx have the biggest effect on the Tx-coil inductance. Since the absolute values of the coefficients dropping after the 13th term rapidly, the model will be reduced to
L ˜ Tx = exp ( 5.72 + 107.83 m 1 r Rx + 43.62 m 1 l Rx + 1.28 n w + 1.31 n l + 170.54 m 1 h a + 29 . 5010 3 p 4.95 m 1 h a p 1 . 7410 3 m 2 r Rx 2 109 . 7010 3 n w 2 110 . 8010 3 n l 2 26 . 4610 3 m 2 h a 2 1 . 010 3 p 2 ) μ H .
The regression coefficients of the full model are listed in Table A2. In Equation (42) the regression coefficients have been transformed in order to allow unnormalized parameter levels (Table 1) inside the equation.
Additionally to the dataset of Table A1, a dataset with 32 random parameter levels is generated, but the values lie still inside the boundaries of Table 1. This dataset will be denoted as test data and is not considered for the calculation of the regression model (see Table A3). Figure 10 shows the relative deviation between the simulated value and the estimation given by the regression model from Equation (42). The predicted inductance values of Equation (42) are in excellent agreement with the simulated results, since the deviation range is only between approximately 2 % to + 2 % for the full model and between 5 % to + 5 % for the reduced model. The very low deviation in case of the training data shows that the chosen exponential polynomial leads to a good representation of the physical behavior.

5.2. Rx-Voltage Amplitude Model

According to the chosen parametrization of the sensor geometry, the same simulation models and design matrix can be used for the regression model of the Rx-coil voltage amplitude. The proposed numerical integration method is applied to the exported field data of the simulated models. It is necessary to calculate only the induced voltage amplitude for one coil due to symmetry. However, the mesh resolution of the FE simulation has to be high enough in the vacuum region of the Rx-coils. A coarse mesh could lead to numerical inaccuracy, especially in further analysis whereby the integration boundaries are shifted for the representation of a target rotation.
Compared to the Tx-coil inductance model, the coefficient values in Figure 11 are not decreasing as rapidly after the four most significant terms. Besides the identity term, the air gap h a has the highest effect to the Rx-voltage amplitude. This is expected because the magnetic field strength decreases exponentially with the distance from the source. Interestingly, the next significant terms are not the parameters r Rx and l Rx who determine the radial dimension, but rather the number of windings n w and layers n l together with the periodicity p of the sensor. As shown with Figure 9, n w and n l have a high influence on the Tx-coil inductance. Correspondingly, the Tx-coil inductance has also an high impact on the Rx-voltage amplitude, since the Tx-current will decrease with an increasing inductance value. An explanation for the effect of the parameter p is the impact on the spanned Rx-coil surface, which can decrease for high periodicities. This is related to the polar sine geometry of the Rx-coils and the increasing slope for higher periodicities, which leads to less coverage of the available area.
Analogously to the Tx-coil inductance model, the model for the Rx-voltage amplitude is reduced to the significant terms displayed in Figure 11:
U ˜ Rx = exp ( 6.37 22.25 m 1 r Rx + 313.41 m 1 l Rx 882.00 × 10 3 n w 1.10 n l 1.02 × 10 3 m 1 h a 30.70 × 10 3 p + 21.40 × 10 3 m 2 r Rx h a + 5.86 m 1 r Rx p 20.11 m 1 l Rx n w + 38.00 × 10 3 m 2 l Rx h a 57.07 m 1 h a p 16.60 × 10 3 m 2 l Rx 2 + 125.40 × 10 3 n w 2 + 127.20 × 10 3 n l 2 4.70 × 10 3 p 2 ) m V .
In Equation (43) the regression coefficients have been transformed in order to allow unnormalized parameter levels (Table 1) inside the equation.
The prediction is less accurate as for the Tx-coil inductance model, but gives still an acceptable result (Figure 12). The relative deviations for the full and reduced models lie between 17 % to + 15 % . Note, that this model relates to a driving voltage amplitude of U Tx = 1 V , which is linearly related to the Rx-voltage amplitude as Equation (6) shows. For different driving voltages the result of Equation (43) can be scaled.

5.3. Angular Error Model

The angular error is a special case, since not only the maximum value is of high interest, but the harmonic composition. In order to predict this behavior a description in form of a Fourier series is introduced with
Δ ϕ p ˜ = H 0 ˜ + k = 1 5 H k ˜ cos ( k ϕ p + δ k ˜ ) ,
where H k ˜ denotes the amplitude and δ k ˜ the phase of the harmonic of order k. A regression polynomial will be calculated for each of those characteristic values, which leads to 11 models in total for a fit up to order k = 5 . Similarly to the other models, data will still be obtained by the numeric integration method, whereby now the signals over one electrical period are needed. Previously, this method was only applied for one target position or reference angle. As mentioned before, only one position is simulated per experiment and the B-field data are exported. The induced voltage signals over one electrical period can be obtained by iterative numeric integration with shifted integration boundaries in the angular direction. With the complete Rx-voltage signals, the electrical angle (Equation (8)) and electrical angular error (Equation (9)) can be calculated. Afterwards, a discrete Fourier transform (DFT) is applied to analyse the harmonic composition of the angular error and obtain the data for the regression analysis.
It turns out that the model approach of Equation (40) leads to an insufficient prediction accuracy. Therefore, the air gap h a is eliminated from the study, whereby the remaining parameters and levels still equals Table 1. The air gap has a high influence on the angular error, since the magnetic field distribution over the Rx-coils changes significantly with this parameter and in turn, the field distribution determines finally the angular error composition. Furthermore, the design needs to be recalculated due to the different number of parameters. The new design contains 25 experiments and needs to be evaluated for each air gap separately (Table A4).
The model approach equals Equation (40). Only the terms containing h a will be set to zero. The target values of the harmonic amplitudes H k will again be logarithmized due to exponential behavior, but the phase information δ k and offset H 0 will still be linear. For readability reasons the whole model equations are not written down. The coefficient values can be obtained from the appendix.
Figure 13 shows nine randomly chosen examples of the test dataset for an air gap h a = 0.5 m m . The left side contains the composed angular error from the regression models, as written in Equation (44). Accordingly, the right side shows only the agreement between the values obtained by numeric integration and the regression models for the harmonic amplitudes. Observing the results, only 0th and 4th order harmonics are present in the angular error. Harmonics of other orders are numerically equal to zero. This is related to the simplified simulation model, which contains only the idealized coil system without any further copper-tracks of a PCB. Considering the prediction accuracy, the mostly dominant 4th harmonic is well fitted with a relative deviation in the range of 23 % to + 11 % . In some cases, the deviation of the offset is huge with up to ± 260 % for very small absolute values. For most applications this is acceptable, since an offset of the angular signal can be compensated.
Table 2 shows exemplarily the comparison for model no. 20 of Table A4. A step size of 1 deg is chosen, which leads in summary to 72 individual simulations. In the parallelized case, 16 target positions were calculated simultaneously. If the available computational resources are limited, only a sequential calculation might be possible, where only one simulation runs at a time. With the simplified model and following numeric integration, only one target position needs to be simulated. To reduce discretization errors, the mesh resolution in the area where the Rx-coils would lie (vacuum) was increased. A lower mesh density could also lead to an acceptable solution accuracy and would require less memory. In the numeric integration, the Rx-surface was discretised in 10 steps in the radial and 180 steps in the angular direction. Regarding computational time, the simplified model is the fastest way of calculation. The additional effort for the following numeric integration is thereby negligible, since the calculation is done in under 10 s.

5.4. Validation with Measurements

As a last step, a proof of the presented models will be given by a measurement of a real sensor-system (see Figure 14). The measurement results will be compared with a simulation model containing the whole PCB-layout and the presented prediction models. Similar to the simulation, one angle is measured at a time without continuous movement of the target (quasi-static). The results of Table 3 are in good agreement with the previous test datasets. For the angular error in Figure 15 follows a good prediction of the dominant 4th harmonic, but a high deviation in the offset. Furthermore, also harmonics up to the 3rd order are present in the measurement and simulation, but cannot be obtained from the regression models. An reasonable explanation for the wrong prediction of these low order harmonics are the not-considered copper tracks of the PCB, since the simulation regarding the full layout is in good agreement with the measurement. However, as the right plots of Figure 15 shows, those harmonics can be eliminated by standard compensation methods—as proposed by Kuntz et al.—and therefore are not highly relevant [20].

6. Conclusions

A novel method for the modeling and simulation of inductive coupled-coil-based angle encoders was introduced to reduce computational effort and calculate regression models for the prediction of sensor quantities without even a FE simulation. The method contains four major steps, starting with the simplification and parametrization of the sensor geometry in order to obtain a generally applicable simulation model. By virtue of the physical sensor principle, the observations focus on six geometric parameters: the intermediate Rx-coil radius r Rx , the Rx peak-to-peak distance l Rx , the number of Tx-windings per layer n w , the number of Tx-layers n l , the air gap between coils and target h a and the periodicity of the sensor p. For the calculation of a regression model, each parameter needs to be varied in defined steps. For instance, 729 simulations would be required if three steps are considered per parameter. In order to reduce this effort, the next step in the proposed method is to find a D-optimal design, which reduces the number of simulations to N = 32 or N = 25 if the air gap is eliminated from the model.
At this point, data can be generated by simulation for the selected parameter variations. Due to the chosen parametrization, the same simulations can be used for the calculation of the Tx-coil inductance and the Rx-voltage amplitudes. Although the Rx-coils were excluded from the simulation model, the resulting induced voltages can still be obtained by a numerical implementation of Faraday’s law. Therefore, an integration method is presented, which is applied to the exported B-field data from the simulation. This method has its biggest benefit when it comes to the calculation of a full electrical period. Normally, this leads to multiple simulations, since every target position needs to be calculated separately. With this method only one simulation is required, since the induced voltages for every angle can be calculated from one field data file by shifting the integration boundaries in the angular direction.
Once the data are acquired, linear regression models can be fitted by the common least-squares method. The deviation between the predicted values from the calculated models and the simulated results were in the range of 2 % to + 2 % for the Tx-inductance and from 17 % to + 15 % for the Rx-voltage amplitudes. A model for the electrical angular error was also calculated, although this is a special case, since a definition of a Fourier series is needed in order to predict the complete signal behavior. As the results show, only 0th and 4th order harmonics appear in the electrical spectrum. The in most cases dominant 4th order harmonic can be predicted with a deviation in the range of 23 % to + 11 % . Unfortunately, the prediction of the offset is less accurate with up to ± 260 % . However according to the measurement, the 0th harmonic can be compensated by simple correction methods. Correspondingly, the 4th is of higher relevance, since simple corrections do not have an impact.
As an outlook, further parameters could be included in the regression models or different parameter ranges might be investigated. Additionally, the applied methods in combination with the calculated regression polynomials leave space for optimization. It is conceivable that the application of optimization algorithms on those models also lead to new insights.

Author Contributions

Conceptualization, R.A.D., G.G. and S.F.; methodology, R.A.D.; software, R.A.D.; validation, R.A.D.; formal analysis, R.A.D.; investigation, R.A.D.; resources, R.A.D. and S.F.; data curation, R.A.D.; writing—original draft preparation, R.A.D.; writing—review and editing, R.A.D., G.G. and S.F.; visualization, R.A.D.; supervision, G.G. and S.F.; project administration, G.G. and S.F.; funding acquisition, S.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Robert Bosch GmbH—Chassis Systems Control, Germany.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. D-Optimal Design for the Tx-Inductance and the Rx-Voltage Amplitude

Table A1 shows the full design matrix in normalized notation and the response values used for the calculation of the regression models of the Tx-coil inductance and the Rx-voltage amplitude. The regression coefficients for the full models are listed in Table A2. The randomly generated test parameter levels with the corresponding response values can be obtained from Table A3.
Table A1. D-optimal design with 32 experiments. The parameter levels are normalized from −1 to +1. In the experiments those levels correspond to the values in Table 1.
Table A1. D-optimal design with 32 experiments. The parameter levels are normalized from −1 to +1. In the experiments those levels correspond to the values in Table 1.
No.Design MatrixTraining Data
I r Rx l Rx n w n l h a p r Rx l Rx r Rx n w r Rx n l r Rx h a r Rx p l Rx n w l Rx n l l Rx h a l Rx p n w n l n w h a n w p n l h a n l p h a p r Rx 2 l Rx 2 n w 2 n l 2 h a 2 p 2 log ( L Tx / μ H ) log ( U Rx / mV )
1+1+1+1−1−1−1−1+1−1−1−1−1−1−1−1−1+1+1+1+1+1+1+1+1+1+1+1+10.5714.011
2+1−1+100+1+1−100−1−100+1+100000+1+1+100+1+11.7390.805
3+1+1+1−1−1+1+1+1−1−1+1+1−1−1+1+1+1−1−1−1−1+1+1+1+1+1+1+10.8072.571
4+1−1−1+1−1−1+1+1−1+1+1−1−1+1+1−1−1−1+1+1−1−1+1+1+1+1+1+11.1622.906
5+1−1−1−1+1−1+1+1+1−1+1−1+1−1+1−1−1+1−1−1+1−1+1+1+1+1+1+11.1872.878
6+1−1+1+1−1−1−1−1−1+1+1+1+1−1−1−1−1−1−1+1+1+1+1+1+1+1+1+11.2633.383
7+1+1−1−1−1−1+1−1−1−1−1+1+1+1+1−1+1+1−1+1−1−1+1+1+1+1+1+10.4993.776
8+1+1−1−1+1−1−1−1−1+1−1−1+1−1+1+1−1+1+1−1−1+1+1+1+1+1+1+11.6573.198
9+1−1−1−1−1−1−1+1+1+1+1+1+1+1+1+1+1+1+1+1+1+1+1+1+1+1+1+1−0.1883.780
10+100+10−1+10000000000−1+100−100+10+1+12.3832.825
11+1+1−1+1+1+1+1−1+1+1+1+1−1−1−1−1+1+1+1+1+1+1+1+1+1+1+1+13.1590.979
12+10+1−1+1+1000000−1+1+10−1−10+1000+1+1+1+101.8222.153
13+1−1+1−1−1−1+1−1+1+1+1−1−1−1−1+1+1+1−1+1−1−1+1+1+1+1+1+10.1843.660
14+1−1−1−1−1+1+1+1+1+1−1−1+1+1−1−1+1−1−1−1−1+1+1+1+1+1+1+10.0321.030
15+1+1+10+1−1+1+10+1−1+10+1−1+1000−1+1−1+1+10+1+1+12.7392.951
16+1+1+1+1+1+1−1+1+1+1+1−1+1+1+1−1+1+1−1+1−1−1+1+1+1+1+1+13.2722.026
17+1−1−1+1+1+10+1−1−1−10−1−1−10+1+10+100+1+1+1+1+102.5480.780
18+1+100−1+1000−1+100000000−100+100+1+101.4322.501
19+1+1−1+1−1−1−1−1+1−1−1−1−1+1+1+1−1−1−1+1+1+1+1+1+1+1+1+11.5803.281
20+1+1+1+1−1−1+1+1+1−1−1+1+1−1−1+1−1−1+1+1−1−1+1+1+1+1+1+11.9503.352
21+1−1−1+1−1+1−1+1−1+1−1+1−1+1−1+1−1+1−1−1+1−1+1+1+1+1+1+11.2162.063
22+1−10+1+1−1−10−1−1+1+10000+1−1−1−1−1+1+10+1+1+1+12.4652.666
23+1−1+1+1+10+1−1−1−10−1+1+10+1+10+10+10+1+1+1+10+12.7861.340
24+1−1+1−1+1−1−1−1+1−1+1+1−1+1−1−1−1+1+1−1−1+1+1+1+1+1+1+11.3243.322
25+1+1−1−10+1−1−1−10+1−1+10−1+10−1+100−1+1+1+10+1+11.3072.283
26+1−10−1+1+1−10+1−1−1+10000−1−1+1+1−1−1+10+1+1+1+11.3602.247
27+1−1+1−1−1+1−1−1+1+1−1+1−1−1+1−1+1−1+1−1+1−1+1+1+1+1+1+10.2153.031
28+1+1−1+1+1−10−1+1+1−10−1−1+10+1−10−100+1+1+1+1+102.9942.591
29+10−10+10−1000000−10+10000−100+10+10+12.2612.227
30+10+1+1−1+1+100000+1−1+1+1−1+1+1−1−1+10+1+1+1+1+11.8221.654
31+1+10−1+10+10−1+10+10000−10−10+10+10+1+10+11.9572.513
32+1+1+1+1000+1+1000+1000000000+1+1+10002.7482.652
Table A2. Regression coefficients for the full model of the Tx-inductance and Rx-voltage amplitude. The coefficients are given in normalized notation. Hence, the parameter levels must be normalized from −1 to +1, according to the boundaries listed in Table 1. Note that those values were calculated for logarithmized data (exponential model).
Table A2. Regression coefficients for the full model of the Tx-inductance and Rx-voltage amplitude. The coefficients are given in normalized notation. Hence, the parameter levels must be normalized from −1 to +1, according to the boundaries listed in Table 1. Note that those values were calculated for logarithmized data (exponential model).
TermRegression Coefficient β i
L ˜ Tx U ˜ Rx
I1.9052.465
r Rx 278.9 × 10 3 195.8 × 10 3
l Rx 109.1 × 10 3 156.2 × 10 3
n w 624.5 × 10 3 280.5 × 10 3
n l 642.4 × 10 3 336.1 × 10 3
h a 66.5 × 10 3 689.1 × 10 3
p 48.6 × 10 3 301.1 × 10 3
r Rx l Rx 10.2 × 10 3 7.6 × 10 3
r Rx n w 76.2 × 10 6 5.1 × 10 3
r Rx n l 3.9 × 10 3 5.7 × 10 3
r Rx h a 3.0 × 10 3 107.0 × 10 3
r Rx p 3.3 × 10 3 117.2 × 10 3
l Rx n w 3.4 × 10 6 50.3 × 10 3
l Rx n l 643.6 × 10 6 29.2 × 10 3
l Rx h a 5.7 × 10 3 97.2 × 10 3
l Rx p 3.6 × 10 3 3.2 × 10 3
n w n l 10.9 × 10 3 37.1 × 10 3
n w h a 2.5 × 10 3 42.8 × 10 3
n w p 6.2 × 10 3 6.7 × 10 3
n l h a 5.2 × 10 3 4.5 × 10 3
n l p 2.4 × 10 3 21.1 × 10 3
h a p 19.8 × 10 3 228.3 × 10 3
r Rx 2 43.4 × 10 3 45.8 × 10 3
l Rx 2 6.8 × 10 3 103.7 × 10 3
n w 2 109.7 × 10 3 125.4 × 10 3
n l 2 110.8 × 10 3 127.2 × 10 3
h a 2 26.5 × 10 3 1.2 × 10 3
p 2 16.0 × 10 3 75.1 × 10 3
Table A3. 32 randomly generated parameter levels inside the boundaries of Table 1. Those data are not used for calculating the regression models and are therefore denoted as test data.
Table A3. 32 randomly generated parameter levels inside the boundaries of Table 1. Those data are not used for calculating the regression models and are therefore denoted as test data.
No.ParametersTest Data
r Rx / mm l Rx / mm n w n l h a / mm p log ( L Tx / μ H ) log ( U Rx / mV )
118.156.38431.0272.5132.539
219.065.23432.1042.5371.936
311.275.49321.3630.7972.889
419.139.12232.3241.4432.716
516.328.47220.8620.4463.781
610.986.59341.0322.0712.580
712.789.75430.7992.3792.536
815.475.17330.7791.7862.750
919.587.19422.2461.9122.307
1019.656.91231.6611.3152.958
1111.588.83431.6032.2562.373
1219.718.98330.7942.1143.142
1319.575.93232.2181.3622.200
1414.857.45331.7411.8382.623
1518.007.23331.2011.9432.885
1611.428.23321.5320.9343.051
1714.228.55321.3061.1542.998
1819.168.77440.6573.1832.605
1917.926.38320.9861.2223.219
2019.598.40340.7552.6442.847
2116.568.28430.8752.4912.765
2210.365.81340.9831.9992.499
2318.495.59421.3371.7542.651
2419.347.49430.6032.5132.958
2516.799.80322.3161.3702.529
2617.586.70342.3922.5421.991
2717.437.93221.4840.5593.393
2813.926.12241.4861.5722.417
2916.558.76341.1872.5722.457
3011.716.28442.3022.7091.655
3117.067.53421.2481.7812.730
3210.328.50230.7270.8813.178

Appendix B. D-Optimal Design for the Angular Error

Analogously to Appendix A, Table A4 shows the full design matrix in normalized notation and the response values used for the calculation of the regression models of the 0th and 4th harmonic amplitude of the angular error. The regression coefficients for the full models are listed in Table A5 and the test data can be obtained from Table A6.
Table A4. D-optimal design with 25 experiments. The parameter levels are normalized from −1 to +1. In the experiments those levels correspond to the values in Table 1, although h a is not considered in this design.
Table A4. D-optimal design with 25 experiments. The parameter levels are normalized from −1 to +1. In the experiments those levels correspond to the values in Table 1, although h a is not considered in this design.
No.Design MatrixTraining Data
I r Rx l Rx n w n l p r Rx l Rx r Rx n w r Rx n l r Rx p l Rx n w l Rx n l l Rx p n w n l n w p n l p r Rx 2 l Rx 2 n w 2 n l 2 p 2 H 0 / deg log ( H 4 / deg )
1+1−1+1−1+1+1−1+1−1−1−1+1+1−1−1+1+1+1+1+1+1−1.110−0.541
2+1+1−1+1−1+1−1+1−1+1−1+1−1−1+1−1+1+1+1+1+1−0.519−2.018
3+1−1−1+1−1−1+1−1+1+1−1+1+1−1−1+1+1+1+1+1+10.005−1.764
4+1+1+1−1−1+1+1−1−1+1−1−1+1+1−1−1+1+1+1+1+1−0.349−0.659
5+1−1−1−1−1+1+1+1+1−1+1+1−1+1−1−1+1+1+1+1+1−0.849−1.607
6+1+1+1+1+1+1+1+1+1+1+1+1+1+1+1+1+1+1+1+1+1−0.316−0.958
7+1+1+1+1−1−1+1+1−1−1+1−1−1−1−1+1+1+1+1+1+10.083−0.296
8+1+1−1+1+1−1−1+1+1−1−1−1+1+1−1−1+1+1+1+1+10.303−2.387
9+10+1+1−1+10000+1−1+1−1+1−10+1+1+1+1−0.647−0.692
10+1+1+1−1+1−1+1−1+1−1−1+1−1−1+1−1+1+1+1+1+10.288−0.529
11+1+1−1−1−1−1−1−1−1−1+1+1+1+1+1+1+1+1+1+1+10.457−1.915
12+1−1+1−1−1−1−1+1+1+1−1−1−1+1+1+1+1+1+1+1+10.129−0.323
13+1−1−1+1+1+1+1−1−1−1−1−1−1+1+1+1+1+1+1+1+1−0.807−1.940
14+1−100+1−100−1+100000−1+100+1+10.178−1.120
15+1+1+10+10+10+100+10000+1+10+10−0.108−0.832
16+1−1+1+1+1−1−1−1−1+1+1+1−1+1−1−1+1+1+1+1+10.071−0.615
17+1−1−1−10−1+1+10+1+10+10+10+1+1+10+10.390−2.359
18+1−1+1+100−1−100+100000+1+1+100−0.484−0.620
19+1−1−1−1+10+1+1−10+1−10−100+1+1+1+10−0.223−2.419
20+100−1−100000000+10000+1+10−0.207−1.002
21+10−100+1000000−10000+100+1−0.473−2.135
22+1−10+1−1+10−1+1−1000−1+1−1+10+1+1+1−1.147−0.916
23+1+1−1−1+1+1−1−1+1+1+1−1−1−1−1+1+1+1+1+1+1−0.158−2.786
24+1−1+10−1+1−10+1−10−1+100−1+1+10+1+1−1.370−0.407
25+1+10+10−10+10−10000−10+10+10+10.173−0.887
Table A5. Regression coefficients for the model of the 0th and 4th harmonic amplitude. The coefficients are given in normalized notation. Hence, the parameter levels must be normalized from −1 to +1, according to the boundaries listed in Table 1. Note that the values for the model H 4 ˜ were calculated for logarithmized data (exponential model).
Table A5. Regression coefficients for the model of the 0th and 4th harmonic amplitude. The coefficients are given in normalized notation. Hence, the parameter levels must be normalized from −1 to +1, according to the boundaries listed in Table 1. Note that the values for the model H 4 ˜ were calculated for logarithmized data (exponential model).
TermRegression Coefficient β i
H 0 ˜ H 4 ˜
I 201.5 × 10 3 1.156
r Rx 198.7 × 10 3 96.9 × 10 3
l Rx 74.2 × 10 3 788.1 × 10 3
n w 61.3 × 10 3 30.9 × 10 3
n l 69.7 × 10 3 191.4 × 10 3
p 442.8 × 10 3 49.7 × 10 3
r Rx l Rx 34.8 × 10 3 43.7 × 10 3
r Rx n w 18.4 × 10 3 11.0 × 10 3
r Rx n l 8.5 × 10 3 37.7 × 10 3
r Rx p 139.5 × 10 3 118.9 × 10 3
l Rx n w 39.2 × 10 3 70.7 × 10 3
l Rx n l 19.1 × 10 3 90.5 × 10 3
l Rx p 7.2 × 10 3 78.2 × 10 3
n w n l 42.5 × 10 3 42.0 × 10 3
n w p 32.7 × 10 3 23.8 × 10 3
n l p 4.5 × 10 3 26.3 × 10 3
r Rx 2 127.9 × 10 3 28.7 × 10 3
l Rx 2 33.7 × 10 3 276.3 × 10 3
n w 2 81.1 × 10 3 47.1 × 10 3
n l 2 52.2 × 10 3 23.6 × 10 3
p 2 45.1 × 10 3 83.4 × 10 3
Table A6. 25 randomly generated parameter levels inside the boundaries of Table 1, but without considering the air gap h a . Those data are not used for calculating the regression models and are therefore denoted as test data.
Table A6. 25 randomly generated parameter levels inside the boundaries of Table 1, but without considering the air gap h a . Those data are not used for calculating the regression models and are therefore denoted as test data.
No.ParametersTest Data
r Rx / mm l Rx / mm n w n l p H 0 / deg log ( H 4 / deg )
118.158.793320.041−0.703
219.068.72327−0.303−0.815
311.276.96333−0.157−1.290
419.138.28235−0.078−1.006
516.325.862320.208−1.743
610.988.53336−0.462−0.872
712.785.16443−0.207−2.259
815.476.38336−0.276−1.581
919.585.23337−0.278−2.045
1019.655.49247−0.085−2.266
1111.589.12435−0.366−0.858
1219.718.473420.101−0.881
1319.576.593430.031−1.546
1414.859.75338−0.452−0.764
1518.005.17432−0.005−1.883
1611.427.19428−0.761−1.059
1714.226.91325−0.315−1.141
1819.168.83239−0.280−0.991
1917.928.982420.154−0.832
2019.595.93345−0.107−1.935
2116.567.45422−0.085−0.867
2210.367.23339−0.905−1.074
2318.498.234310.148−0.717
2419.348.55227−0.234−0.813
2516.798.77438−0.394−0.972

References

  1. Nagel, J.R. Induced eddy currents in simple conductive geometries: Mathematical formalism describes the excitation of electrical eddy currents in a time-varying magnetic field. IEEE Antennas Propag. Mag. 2017, 60, 81–88. [Google Scholar] [CrossRef]
  2. Gianchandani, Y.B.; Crary, S.B. Parametric modeling of a microaccelerometer: Comparing I-and D-optimal design of experiments for finite-element analysis. J. Microelectromechanical Syst. 1998, 7, 274–282. [Google Scholar] [CrossRef]
  3. Passarotto, M.; Qama, G.; Specogna, R. A fast and efficient simulation method for inductive position sensors design. In Proceedings of the 2019 IEEE SENSORS, Montreal, QC, Canada, 27–30 October 2019; pp. 1–4. [Google Scholar]
  4. Hoxha, A.; Passarotto, M.; Qama, G.; Specogna, R. Design Optimization of PCB-Based Rotary-Inductive Position Sensors. Sensors 2022, 22, 4683. [Google Scholar] [CrossRef] [PubMed]
  5. Schwab, A.J. Begriffswelt der Feldtheorie; Springer: Berlin/Heidelberg, Germany, 2019; Volume 8. [Google Scholar]
  6. Jiao, D.; Ni, L.; Zhu, X.; Zhe, J.; Zhao, Z.; Lyu, Y.; Liu, Z. Measuring gaps using planar inductive sensors based on calculating mutual inductance. Sens. Actuators A Phys. 2019, 295, 59–69. [Google Scholar] [CrossRef]
  7. Weber, H.; Baran, H.; Utermöhlen, F.; Schuster, C. Macromodeling of Mutual Inductance for Displaced Coils Based on Laplace’s Equation. IEEE Trans. Instrum. Meas. 2021, 70, 1–11. [Google Scholar] [CrossRef]
  8. Ye, L. Characteristics analysis of self-sensing double-coil eddy current damper. In Proceedings of the 2017 IEEE 2nd Information Technology, Networking, Electronic and Automation Control Conference (ITNEC), Chengdu, China, 15–17 December 2017; pp. 1581–1585. [Google Scholar]
  9. Cauer, E.; Mathis, W.; Pauli, R. Life and work of wilhelm cauer (1900 1945). In Proceedings of the 14th International Symposium of Mathematical Theory of Networks and Systems (MTNS 2000), Perpignan, France, 19–23 June 2000; pp. 1–10. [Google Scholar]
  10. Li, Y.; Berthiau, G.; Feliachi, M.; Cheriet, A. 3D finite volume modeling of ende using electromagnetic T-formulation. J. Sens. 2012, 2012, 785271. [Google Scholar] [CrossRef] [Green Version]
  11. Granula, B.; Porzig, K.; Toepfer, H.; Gacanovic, M. A Comparison between an AV and V Formulation in Transcranial Magnetic Stimulation. In Proceedings of the Comsol Conference, Rotterdam, The Netherlands, 23–25 October 2013; pp. 23–25. [Google Scholar]
  12. Theodoulidis, T.P.; Kantartzis, N.V.; Tsiboukis, T.D.; Kriezis, E.E. Analytical and numerical solution of the eddy-current problem in spherical coordinates based on the second-order vector potential formulation. IEEE Trans. Magn. 1997, 33, 2461–2472. [Google Scholar] [CrossRef] [Green Version]
  13. Davis, P.J.; Rabinowitz, P. Methods of Numerical Integration; Dover Publications, Inc.: Mineola, NY, USA, 2007. [Google Scholar]
  14. Fahrmeir, L.; Kneib, T.; Lang, S.; Marx, B.D. Regression Models; Springer: Berlin/Heidelberg, Germany, 2021; pp. 23–84. [Google Scholar]
  15. de Aguiar, P.F.; Bourguignon, B.; Khots, M.; Massart, D.; Phan-Than-Luu, R. D-optimal designs. Chemom. Intell. Lab. Syst. 1995, 30, 199–210. [Google Scholar] [CrossRef]
  16. Fedorov, V.V. Theory of Optimal Experiments; Academic Press: Cambridge, MA, USA, 1972. [Google Scholar]
  17. Friendly, M.; Monette, G.; Fox, J. Elliptical insights: Understanding statistical methods through elliptical geometry. Stat. Sci. 2013, 28, 1–39. [Google Scholar] [CrossRef]
  18. Özdemir, A.; Turkoz, M. Development of a D-optimal design-based 0–1 mixed-integer nonlinear robust parameter design optimization model for finding optimum design factor level settings. Comput. Ind. Eng. 2020, 149, 106742. [Google Scholar] [CrossRef]
  19. Altun, E. The log-weighted exponential regression model: Alternative to the beta regression model. Commun. Stat.—Theory Methods 2021, 50, 2306–2321. [Google Scholar] [CrossRef]
  20. Kuntz, S.; Dauth, R.; Gerlach, G.; Ott, P.; Fella, S. Harmonic analysis of the arctangent function regarding the angular error introduced by superimposed Fourier series for application in sine/cosine angle encoders. Sens. Actuators A Phys. 2022, 344, 113585. [Google Scholar] [CrossRef]
Figure 1. Overview of the sensor coil system (a)—including a target—and the signal flow from encoder signals (b) to the calculated angle and error (d) by calculation of the unwrapped arctangent function (c). Note that the signals a and b are already demodulated. The figure does not consider the signal-processing circuitry.
Figure 1. Overview of the sensor coil system (a)—including a target—and the signal flow from encoder signals (b) to the calculated angle and error (d) by calculation of the unwrapped arctangent function (c). Note that the signals a and b are already demodulated. The figure does not consider the signal-processing circuitry.
Sensors 22 07804 g001
Figure 2. Modeling approach starting with a parametrization and simplification of the sensor geometry (a), followed by application of statistical methods (b), data generation (c) and calculation of regression models (d). The geometric parameters of the sensor model are the intermediate Rx-coil radius r Rx , the Rx peak-to-peak distance l Rx , the number of Tx-windings per layer n w , the number of Tx-layers n l , the air gap between coils and target h a and the periodicity of the sensor p. For the statistical methods the design matrix X , the response variable y and the regression coefficients β must be introduced.
Figure 2. Modeling approach starting with a parametrization and simplification of the sensor geometry (a), followed by application of statistical methods (b), data generation (c) and calculation of regression models (d). The geometric parameters of the sensor model are the intermediate Rx-coil radius r Rx , the Rx peak-to-peak distance l Rx , the number of Tx-windings per layer n w , the number of Tx-layers n l , the air gap between coils and target h a and the periodicity of the sensor p. For the statistical methods the design matrix X , the response variable y and the regression coefficients β must be introduced.
Sensors 22 07804 g002
Figure 3. Equivalent circuit diagram of the coil system of the inductive coupled coil sensor. The oscillator driver for the Tx-coil is replaced by an AC-voltage source. This is a valid approximation which is also done in most simulation models. Note, that for simplicity just one Rx-coil circuit is drawn. The ohmic resistance of the Rx-coil is neglected since the amplifier stage input for the voltage measurement has a high impedance to avoid current flow.
Figure 3. Equivalent circuit diagram of the coil system of the inductive coupled coil sensor. The oscillator driver for the Tx-coil is replaced by an AC-voltage source. This is a valid approximation which is also done in most simulation models. Note, that for simplicity just one Rx-coil circuit is drawn. The ohmic resistance of the Rx-coil is neglected since the amplifier stage input for the voltage measurement has a high impedance to avoid current flow.
Sensors 22 07804 g003
Figure 4. Calculated encoder angle in dependence of the mechanical reference angle ϕ for a sensor with periodicity p = 4 . For comparison, the angle is calculated by, atan2, and unwrapped atan2 functions. An exaggerated angular error deviating from the ideal output is shown as an example.
Figure 4. Calculated encoder angle in dependence of the mechanical reference angle ϕ for a sensor with periodicity p = 4 . For comparison, the angle is calculated by, atan2, and unwrapped atan2 functions. An exaggerated angular error deviating from the ideal output is shown as an example.
Sensors 22 07804 g004
Figure 5. Sensor geometry shown as an example with periodicity p = 5 . In the radial dimension, the intermediate radius r Rx and the peak-to-peak distance l Rx of the Rx-coil define primarily the overall dimensions. The radii of the Tx-coil and target follow implicitly. Note that this is a simplified sketch, containing neither coil connections or vias.
Figure 5. Sensor geometry shown as an example with periodicity p = 5 . In the radial dimension, the intermediate radius r Rx and the peak-to-peak distance l Rx of the Rx-coil define primarily the overall dimensions. The radii of the Tx-coil and target follow implicitly. Note that this is a simplified sketch, containing neither coil connections or vias.
Sensors 22 07804 g005
Figure 6. Surface definition for numeric integration of simulated B-field data for the Rx-coils. A surface-grid is defined, whereby only elements inside the boundaries—determined by the centerlines C + and C of the Rx-coil tracks—contribute to the overall magnetic flux. The 3D-view of a cutout of the Rx-coil with the surface normals for the relevant elements is displayed in (a). Due to the complex geometry of the coil with crossing tracks on different layers, the surface rotates around the intermediate radius t Rx . Accordingly, a parametrization of the unwound surface A ( t , ψ ) is introduced. (b) shows the discretization of this surface with low resolution for better illustration. The crossed surface-elements lie inside the integration boundaries and contribute therefore to the total magnetic flux. For the transformation in ( t , ψ ) coordinates, the inclination angle ξ for every position at angle φ needs to be known. (c) shows the cross-section in the rz-plane. The calculation of ξ and the following parametrization requires only simple trigonometric expressions as shown in Equations (21) and (20).
Figure 6. Surface definition for numeric integration of simulated B-field data for the Rx-coils. A surface-grid is defined, whereby only elements inside the boundaries—determined by the centerlines C + and C of the Rx-coil tracks—contribute to the overall magnetic flux. The 3D-view of a cutout of the Rx-coil with the surface normals for the relevant elements is displayed in (a). Due to the complex geometry of the coil with crossing tracks on different layers, the surface rotates around the intermediate radius t Rx . Accordingly, a parametrization of the unwound surface A ( t , ψ ) is introduced. (b) shows the discretization of this surface with low resolution for better illustration. The crossed surface-elements lie inside the integration boundaries and contribute therefore to the total magnetic flux. For the transformation in ( t , ψ ) coordinates, the inclination angle ξ for every position at angle φ needs to be known. (c) shows the cross-section in the rz-plane. The calculation of ξ and the following parametrization requires only simple trigonometric expressions as shown in Equations (21) and (20).
Sensors 22 07804 g006
Figure 7. Example of two correlated random variables y 1 and y 2 . The random variables are jointly normally distributed which is shown by the contours (a) of the underlying bivariate normal distribution (b). Accordingly, this plot elucidates the similarity to an ellipsoid, whereat the eigenvectors lying in the direction of the semi-axes.
Figure 7. Example of two correlated random variables y 1 and y 2 . The random variables are jointly normally distributed which is shown by the contours (a) of the underlying bivariate normal distribution (b). Accordingly, this plot elucidates the similarity to an ellipsoid, whereat the eigenvectors lying in the direction of the semi-axes.
Sensors 22 07804 g007
Figure 8. Comparision of the D-efficiency of designs with different number of experiments. The compared designs containing the terms of Equation (40).
Figure 8. Comparision of the D-efficiency of designs with different number of experiments. The compared designs containing the terms of Equation (40).
Sensors 22 07804 g008
Figure 9. Pareto-diagram related to the normalized regression coefficients of the Tx-coil inductance model. It shows the 13 most significant coefficients in descending order. The term denoted by I represents the identity term, which is related to β 0 . Please note that those values correspond to the exponential model as shown in Equation (42).
Figure 9. Pareto-diagram related to the normalized regression coefficients of the Tx-coil inductance model. It shows the 13 most significant coefficients in descending order. The term denoted by I represents the identity term, which is related to β 0 . Please note that those values correspond to the exponential model as shown in Equation (42).
Sensors 22 07804 g009
Figure 10. Relative deviation between the training or test data and the regression model of the Tx-coil inductance. The full-model (a)—considering all coefficients from Equation (40)—is compared with the reduced model (b) from Equation (42). The training data respond to the experiments of Table A1. The test dataset contains also 32 experiments with randomly generated parameter levels inside the boundaries listed in Table 1.
Figure 10. Relative deviation between the training or test data and the regression model of the Tx-coil inductance. The full-model (a)—considering all coefficients from Equation (40)—is compared with the reduced model (b) from Equation (42). The training data respond to the experiments of Table A1. The test dataset contains also 32 experiments with randomly generated parameter levels inside the boundaries listed in Table 1.
Sensors 22 07804 g010
Figure 11. Pareto-diagram related to the normalized regression coefficients of the Rx-voltage amplitude model. It shows the 16 most significant coefficients in descending order. The term denoted by I represents the identity term, which is related to β 0 . Please note that those values correspond to the exponential model as shown in Equation (43).
Figure 11. Pareto-diagram related to the normalized regression coefficients of the Rx-voltage amplitude model. It shows the 16 most significant coefficients in descending order. The term denoted by I represents the identity term, which is related to β 0 . Please note that those values correspond to the exponential model as shown in Equation (43).
Sensors 22 07804 g011
Figure 12. Relative deviation between the training or test data and the regression model of the Rx-voltage amplitude. The full-model (a)—considering all coefficients from Equation (40)—is compared with the reduced model (b) from Equation (43). The training data respond to the experiments of Table A1. The test dataset contains also 32 experiments with randomly generated parameter levels inside the boundaries listed in Table 1.
Figure 12. Relative deviation between the training or test data and the regression model of the Rx-voltage amplitude. The full-model (a)—considering all coefficients from Equation (40)—is compared with the reduced model (b) from Equation (43). The training data respond to the experiments of Table A1. The test dataset contains also 32 experiments with randomly generated parameter levels inside the boundaries listed in Table 1.
Sensors 22 07804 g012
Figure 13. Comparison of the electrical angular error signal and harmonics of nine randomly chosen examples of the test dataset. For each harmonic of the orders 0 to 5 a regression model is calculated for the harmonic amplitude and phase. The right plots show the harmonic amplitudes (b), whereby the left plots show the total error curve (a) according to Equation (44). The data were calculated for an air gap of h a = 0.5 m m and the index of each plot title (in brackets) refers to Table A6.
Figure 13. Comparison of the electrical angular error signal and harmonics of nine randomly chosen examples of the test dataset. For each harmonic of the orders 0 to 5 a regression model is calculated for the harmonic amplitude and phase. The right plots show the harmonic amplitudes (b), whereby the left plots show the total error curve (a) according to Equation (44). The data were calculated for an air gap of h a = 0.5 m m and the index of each plot title (in brackets) refers to Table A6.
Sensors 22 07804 g013
Figure 14. Overview about the test setup. The sensor PCB is mounted on a linear axis for adjusting the air gap. On the opposite side is the target. The rotation is performed by a stepper motor. This motor also includes an encoder system, which gives the reference angle for the measurement. The measurement procedure is quasi-static. The rotation will stop at each angle, while the output voltages are measured with a digital multimeter.
Figure 14. Overview about the test setup. The sensor PCB is mounted on a linear axis for adjusting the air gap. On the opposite side is the target. The rotation is performed by a stepper motor. This motor also includes an encoder system, which gives the reference angle for the measurement. The measurement procedure is quasi-static. The rotation will stop at each angle, while the output voltages are measured with a digital multimeter.
Sensors 22 07804 g014
Figure 15. Comparison of the electrical angular error signals (a,c) and harmonics (b,d) between measurement, simulation containing the full PCB-layout and prediction of the presented regression models for experiment no. 4 of the training dataset for an air gap of 0.5 m m . The left side (a,b) shows the raw results. On the right side (c,d) an amplitude-, offset- and orthogonality-correction was applied to the induced voltage signals and a mean compensation to the predicted angular error in case of the regression model.
Figure 15. Comparison of the electrical angular error signals (a,c) and harmonics (b,d) between measurement, simulation containing the full PCB-layout and prediction of the presented regression models for experiment no. 4 of the training dataset for an air gap of 0.5 m m . The left side (a,b) shows the raw results. On the right side (c,d) an amplitude-, offset- and orthogonality-correction was applied to the induced voltage signals and a mean compensation to the predicted angular error in case of the regression model.
Sensors 22 07804 g015
Table 1. Parameters describing the sensor geometry. For each parameter three levels are required, since the maximum order of the regression polynomial is two.
Table 1. Parameters describing the sensor geometry. For each parameter three levels are required, since the maximum order of the regression polynomial is two.
VariableLevelsQuantity
r Rx [ 10 , 15 , 20 ] m m Intermediate Rx-coil radius
l Rx [ 5 , 7.5 , 10 ] m m Rx-coil peak-to-peak distance
n w [ 2 , 3 , 4 ] No. of Tx-windings per layer
n l [ 2 , 3 , 4 ] No. of Tx-layers
h a [ 0.5 , 1.5 , 2.5 ] m m Air gap
p [ 1 , 5 , 9 ] Periodicity
Table 2. Computational resources for no. 20 of Table A4 evaluated by a sequential and parallelized simulation of the full model over one electrical period and one simulation of the simplified model with following numeric integration.
Table 2. Computational resources for no. 20 of Table A4 evaluated by a sequential and parallelized simulation of the full model over one electrical period and one simulation of the simplified model with following numeric integration.
SequentialParallelizedSimplified
Comp. time 138 h 9 h 35 min 4 h 01 min
Max. RAM 26.8 GB 428.8 GB 77.7 GB
CPU cores1161
Table 3. Comparison between measurement results and prediction of the models according to Equations (42) and (43). A scaling factor of 1.2 was introduced to the model of the Rx-voltage amplitude, since the Tx driving-voltage amplitude equals 1.2 V for the tested sample.
Table 3. Comparison between measurement results and prediction of the models according to Equations (42) and (43). A scaling factor of 1.2 was introduced to the model of the Rx-voltage amplitude, since the Tx driving-voltage amplitude equals 1.2 V for the tested sample.
MeasurementModelRelative Deviation
L Tx 2.12 μ H 2.05 μ H 3.30 %
U Rx 67.61 m V 71.13 m V + 5.21 %
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dauth, R.A.; Gerlach, G.; Fella, S. An Effective Method to Model and Simulate the Behavior of Inductive Angle Encoders. Sensors 2022, 22, 7804. https://doi.org/10.3390/s22207804

AMA Style

Dauth RA, Gerlach G, Fella S. An Effective Method to Model and Simulate the Behavior of Inductive Angle Encoders. Sensors. 2022; 22(20):7804. https://doi.org/10.3390/s22207804

Chicago/Turabian Style

Dauth, Robert Alexander, Gerald Gerlach, and Sina Fella. 2022. "An Effective Method to Model and Simulate the Behavior of Inductive Angle Encoders" Sensors 22, no. 20: 7804. https://doi.org/10.3390/s22207804

APA Style

Dauth, R. A., Gerlach, G., & Fella, S. (2022). An Effective Method to Model and Simulate the Behavior of Inductive Angle Encoders. Sensors, 22(20), 7804. https://doi.org/10.3390/s22207804

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

Article Metrics

Back to TopTop