Next Article in Journal
Effect of Time and Frequency of Magnetic Field Application on MRF Pressure Performance
Next Article in Special Issue
Experimental Research and Multi-Physical Field Coupling Simulation of Electrochemical Machining Based on Gas–Liquid Two-Phase Flow
Previous Article in Journal
Editorial for the Special Issue on “Frontiers of Ultra-Precision Machining”
Previous Article in Special Issue
Integration of Multifocal Microlens Array on Silicon Microcantilever via Femtosecond-Laser-Assisted Etching Technology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Main Error Sources for the Error Motion Measurement of a Precision Shafting Using a T-Type Capacitive Sensor

1
Key Laboratory of Advanced Manufacturing Technology of Zhejiang Province, School of Mechanical Engineering, Zhejiang University, Hangzhou 310027, China
2
School of Mechanical Engineering, Hangzhou Dianzi University, Hangzhou 310018, China
*
Author to whom correspondence should be addressed.
Micromachines 2022, 13(2), 221; https://doi.org/10.3390/mi13020221
Submission received: 7 December 2021 / Revised: 22 January 2022 / Accepted: 26 January 2022 / Published: 29 January 2022
(This article belongs to the Special Issue Advances in Ultra-Precision Machining Technology and Applications)

Abstract

:
As a key indicator reflecting the working accuracy of rotary functional units, the error motions of the precision shafting are very necessary to be measured. In this paper, the main error sources for the error motion measurement of a precision shafting using a T-type capacitive sensor were investigated. The theoretical modeling error due to the approximate simplification for the output capacitance expressions was firstly analyzed. By means of the 3D-FEA method, the influence of fringe effects was subsequently investigated. Finally, the analysis of electrode installation errors was emphasized on the tilt error of the cylindrical electrode and coaxiality error of the fan-shaped electrode by establishing mathematical models and numerical simulation. Based on the theoretical analysis and simulation results, the methods of decreasing the approximate error and the nonlinear error caused by fringe effects were subsequently proposed; for the installation errors, the tilt error of cylindrical electrode only makes the solution of phase angle have a certain deviation and has almost no effect on solving the radial displacement, especially for the measurement range less than 0.1 mm; the measurement of the rotor tilt displacement was basically not affected by the coaxiality error of the fan-shaped electrode.

1. Introduction

The rotary functional units with precision shafting as the core are widely used in high-end equipment such as high-precision CNC machine tools, EUV lithography equipment, axial compressors and small agile satellites [1,2,3,4,5,6,7]. As a key indicator reflecting the working accuracy of rotary functional units, the error motions of the precision shafting are very necessary to be measured. This measurement provides valuable insight into the error motions and performance [8], which is favorable to optimize these rotary units.
Capacitive sensors with good stability, high speed and adaptability of extreme environments are widely adopted to measure the displacement of motion targets [9,10,11,12]. Chapman et al. carried out research on the measurement of rotary axis error motions using a cylindrical capacitive sensor (CCS) with multiple electrodes of large sense area [13,14,15]. In order to measure the multi-degree-of-freedom (multi-DOF) error motions, several separate electrodes need to be manufactured and orderly mounted along with the guard ring, which is difficult.
Compared to the CCS used in the above research, a novel T-type capacitive sensor (T-type CS) with an integrated structure was developed to achieve the simultaneous online measurement of the 5-DOF error motions of a spindle, as shown in Figure 1 and Figure 2 [16]. Both the radial electrode group (REG) and the end part electrode group (EPEG) were fabricated based on a flexible printed circuit board, which could produce a quite thinner electrode (0.035 mm) and a tiny gap (0.15 mm) between the electrode and guard ring.
This paper focuses on the analysis of main error sources for the error motion measurement of a precision shafting using a T-type CS, as shown in Figure 3. Firstly, the theoretical modeling error due to the approximate simplification for the output capacitance expression of sensing electrodes was analyzed. Then, the influence of fringe effects was investigated using the 3D-FEA method. Finally, mathematical modeling and numerical simulation were conducted to analyze electrode installation errors, emphasizing the tilt error of the cylindrical electrode and coaxiality error of the fan-shaped electrodes.

2. Basic Principle

2.1. Radial Error Motion Measurement

This measurement is performed by the REG, which is constituted by four cylindrical electrodes Er1~Er4. The four cylindrical electrodes and the radial measuring surface (CSR) of the rotor make up four cylindrical capacitors, respectively; the corresponding capacitance are denoted by C1~C4, as shown in Figure 4 [16].
Due to the arc-shaped boundary of these electrodes, the conformal mapping method was adopted to solve the Poisson Equation in electrostatic fields [17]. When the rotor has radial displacement with eccentricity α at phase angle β, caused by the generated radial run-out of a spindle, the resulting variation Δu of the equivalent spacing at any point in the complex plane t can be approximately derived as [16]:
Δ u = α cos ( θ β ) R
where R is the inside radius of the cylindrical electrode, θ is the argument of the complex number corresponding to any point on the arc-shaped boundary in the plane z.
Thereby, the differential expression of the output capacitance of a cylindrical electrode with axial width of w can be illustrated as [16]:
Δ C r = ε 0 ε r R Δ θ w g α cos ( θ β )
where ε0 is the electric permittivity of vacuum, εr is the relative dielectric constant; g is the initial radial spacing and Δθ is an equivalent length of a micro-unit that is taken on the equivalent boundary in the plane t.
Then the output capacitance of the cylindrical electrode in each quadrant is approximately expressed as follows [16]:
C 1 = ε 0 ε r R w ζ 1 ζ 2 1 g α cos ( θ β ) d θ
C 2 = ε 0 ε r R w π ζ 2 π ζ 1 1 g α cos ( θ β ) d θ
C 3 = ε 0 ε r R w π + ζ 1 π + ζ 2 1 g α cos ( θ β ) d θ
C 4 = ε 0 ε r R w 2 π ζ 2 2 π ζ 1 1 g α cos ( θ β ) d θ
where (ζ2ζ1) denotes the angular size of the cylindrical electrode.
Based on Equations (3)–(6), the expressions of the differential output capacitance of the REG can be expressed by:
CrX = C 1 + C 4 C 2 C 3
CrY = C 1 + C 2 C 3 C 4
where CrX and CrY, respectively, correspond to the components δx and δy of the displacement α in the X and Y directions.
Then, the radial error motions of the spindle can be estimated by the following equations:
δ x = f x ( C 1 + C 4 C 2 C 3 )
δ y = f y ( C 1 + C 2 C 3 C 4 )
where fx and fy are denoted the transition functions between the measured capacitance of the cylindrical electrodes and the radial run-out along the x, y directions, respectively.

2.2. Axial Error Motion Measurement

As shown in Figure 2, the fan-shaped electrodes (Ee1~Ee4) of the EPEG are also distributed in the same configuration as radial electrodes. They, along with the inside annular end face of the rotor flange (FR), make up four capacitors, respectively; the corresponding capacitance of them are denoted by C5~C8, as shown in Figure 5 [16].
When the rotor has an axial error motion along the z-axis, the variation quantity of the sum of the output capacitance of the fan-shaped electrode is given by [16]:
Δ CeZ = CeZ CeZ = 4 ε 0 ε r S e ( 1 g a 0 + Δ z 1 g a 0 )
where Se is the area of the fan-shaped electrode, ga0 is the initial value of the axial spacing, the variation of the axial spacing Δz is defined by Δz = gaga0, CeZ is the sum of the output capacitance of the fan-shaped electrodes under the axial spacing ga0 and CeZ′ is that of the fan-shaped electrodes under the axial spacing ga.
According to Equation (11), the expressions for estimating the axial error motion of the spindle can be derived as:
δ z = f z ( C 5 + C 6 + C 7 + C 8 )

2.3. Tilt Error Motion Measurement

The calculation model for the measurement of tilt error motions by EPEG is shown in Figure 6 [16]. The line segment MN represents the axial measuring surface of the rotor (i.e., the inside annular end face of the rotor flange), and its perpendicular bisector is the medial axis of the rotor.
In the boundary region of the fan-shaped electrode, a micro-plane element ΔBP at any point P is taken as the solving unit (the area of it is denoted as ΔSP). The micro-plane element ΔBP and its opposite counterpart ΔAQ compose a plane parallel capacitor (ΔAQ is in the axial measuring surface of the rotor), and the spacing between them is derived as [16]:
t = a k ρ cos ( φ γ )
where a is the distance from the intersection point Or to the EPEG; k = tan(ε), ε and γ are the amplitude and yaw angle of the tilt displacement of the rotor, respectively.
Thereby, the capacitance of the plane parallel capacitor is [16]:
Δ C e = ε 0 ε r Δ S P t = ε 0 ε r ρ d φ d ρ a k ρ cos ( φ γ )
Then, the output capacitance of the fan-shaped electrode Ee4 can be approximately obtained by integrating Equation (14) over the area of the fan-shaped electrode [16]:
C 8 = ε 0 ε r 2 π φ 2 2 π φ 1 ρ 1 ρ 2 ρ d ρ a k ρ cos ( φ + γ ) d φ
where (φ2φ1) are the angular size of each fan-shaped electrode, and ρ1 and ρ2 are the inside and outside radius of the electrode, respectively.
By considering the minor tilt displacement of a precision spindle, the power series expansion by ignoring higher-order terms was applied to express the Equation (15) (k << a):
C 8 = ε 0 ε r a 2 π φ 2 2 π φ 1 [ b + k c a cos ( φ + γ ) ] d φ
where the factors b and c are given by b = ρ 2 2 ρ 1 2 2 , c = ρ 2 3 ρ 1 3 3 .
By analogy with the electrode Ee4, the output capacitance of the other electrodes can be approximately expressed as follows:
C 7 = ε 0 ε r a π + φ 1 π + φ 2 [ b + k c a cos ( φ + γ ) ] d φ
C 6 = ε 0 ε r a π φ 2 π φ 1 [ b + k c a cos ( φ + γ ) ] d φ
C 5 = ε 0 ε r a φ 1 φ 2 [ b + k c a cos ( φ + γ ) ] d φ
By making an arrangement, Equations (16)–(19) are given as:
C 8 = ε 0 ε r a φ 1 φ 2 [ b + k c a cos ( φ γ ) ] d φ
C 7 = ε 0 ε r a φ 1 φ 2 [ b k c a cos ( φ + γ ) ] d φ
C 6 = ε 0 ε r a φ 1 φ 2 [ b k c a cos ( φ γ ) ] d φ
C 5 = ε 0 ε r a φ 1 φ 2 [ b + k c a cos ( φ + γ ) ] d φ
By combining Equations (20)–(23), the distance a can be derived as:
a = 4 ε 0 ε r b ( φ 2 φ 1 ) C 8 + C 7 + C 6 + C 5
Based on Equations (16)–(19), the expressions of the differential output capacitance of the EPEG can be expressed as follows:
CeX = C 5 + C 6 C 7 C 8
CeY = C 5 + C 8 C 6 C 7
where CeX and CeY, respectively, correspond to the components εx (around the x-axis), εy (around the y-axis) of the tilt displacement ε.
Then, the tilt error motions of the spindle can be estimated by the following equations:
ε x = f e x ( C 5 + C 6 C 7 C 8 )
ε y = f e y ( C 5 + C 8 C 6 C 7 )
where fex and fey are denoted as the transition functions between the measured capacitance of the fan-shaped electrodes and the tilt displacements about the x- and y-axis, respectively.

3. Main Error Sources Analysis of the T-Type CS

3.1. Theoretical Modeling Error Analysis

In order to theoretically analyze the radial error motion measurement with the T-type CS, we consider a small radial displacement of the precision spindle, i.e., |α/(Rg)| < 1. The power series expansion with first-order approximation (Equation (1)) was applied to express the equivalent spacing changing Δu. The exact derived expression of it is:
Δ u o = 1 2 ln { ( R g ) 2 + 2 α ( R g ) cos ( θ β ) + α 2 ( R g ) 2 }
Then, the approximate error of Equation (1) can be assessed by:
Δ u a p p . e r r o r = | Δ u o Δ u Δ u o |
Figure 7 shows the variation situation of error Δuapp.error with the changing of parameter terms α/(Rg) and cos(θβ). As mentioned above, the radial displacement α of a precision spindle is generally very small. In order to achieve relatively high sensitivity, the cylindrical electrode radius R is often designed to be tens of millimeters or even larger. Thus, the variation range of the term α/(Rg) is set to be [0, 0.1]. It can be observed that, when the ratio α/(Rg) is less than 0.01, the error Δuapp.error basically keeps lower than 2% with the variation in the term cos(θβ) in the entire phase region (the value range of cos(θβ) is [0, 1] correspondingly), except for a few phase angles. Therefore, a highly accurate value of the changing Δu calculated by Equation (1) can be obtained if the ratio α/(Rg) is less than 0.01, which can be achieved by adjusting the structural parameters of the sensor.
Since a highly accurate value of Δu can be obtained by setting the ratio α/(R − g) ≤ 0.01, the expressions (3)–(6) established based on Equation (1) could be more accurate to reflect the capacitance changing caused by the radial displacement, and is also for the expressions of the differential output capacitance of the REG derived as:
CrX = 2 ε 0 ε r R w g 1 ( α g ) 2 ( arctan 2 α g 1 ( α g ) 2 sin ζ 2 cos β 1 ( α g ) 2 ( cos 2 β + sin 2 ζ 2 ) arctan 2 α g 1 ( α g ) 2 sin ζ 1 cos β 1 ( α g ) 2 ( cos 2 β + sin 2 ζ 1 ) )
CrY = 2 ε 0 ε r R w g 1 ( α g ) 2 ( arctan 2 α g 1 ( α g ) 2 cos ζ 1 sin β 1 + ( α g ) 2 ( cos 2 β + sin 2 ζ 1 2 ) arctan 2 α g 1 ( α g ) 2 cos ζ 2 sin β 1 + ( α g ) 2 ( cos 2 β + sin 2 ζ 2 2 ) )
In a manner similar to the simplification process on the equivalent spacing changing Δu, the first-order approximation of the Equations (31) and (32) could intuitively reflect the variation’s relation of the radial displacement and the differential capacitance (CrX and CrY). Thus, the expressions of the displacement α and phase angle β can be deduced in the following form:
CrX tl = 4 ε 0 ε r R w ( sin ζ 2 sin ζ 1 ) g 2 α cos β
CrY tl = 4 ε 0 ε r R w ( cos ζ 1 cos ζ 2 ) g 2 α sin β
Therefore, the expression used to assess the approximate error (taking the CrX as an example) is given as:
CrX a p p . e r r o r = | CrX CrX tl CrX |
As shown in Figure 8, at each phase angle β, the approximate error CrXapp.error all increases with the increase in the ratio α/g. When the ratio α/g becomes larger than 0.2, the variation tendency is more significant. Figure 9 shows the error of the calculated capacitance CrXtl relative to the simulation value CrXf as the rotor has a displacement along the direction of y = x. The dependence of the relative error on the parameter α is basically consistent with the variation tendency shown in Figure 8. Note that the value CrXf was obtained based on the simulation model of the T-type CS previously established with the design parameters in [16]; these parameters were also adopted to calculate the theoretical capacitance.
From the analysis in Figure 8, it can be known that diminishing the ratio α/g could make the error CrXapp.error reduce further, but the measurement range and sensitivity of the sensor will be affected to some extent. By considering the inherent nonlinearity of the capacitor output with spacing change, the power series expansion with a third-order approximation was derived to improve the calculation accuracy [18]:
CrX tl 3 = 4 ε 0 ε r R w ( sin ζ 2 sin ζ 1 ) g ( α g + 0.95 ( α g ) 3 ) cos β
CrY tl 3 = 4 ε 0 ε r R w ( cos ζ 1 cos ζ 2 ) g ( α g + 0.95 ( α g ) 3 ) sin β
Compared to the error distribution shown in Figure 8, the error CrXapp.error is significantly decreased overall for the calculation expression of the CrX with third-order approximation, as shown in Figure 10. As the phase angle β varies from 0° to 90°, the error CrXapp.error remains less than 0.4% for the ratio α/g = 0.1 and no more than 1.2% for α/g = 0.2. Figure 11 compares the simulated differential capacitance CrXf and theoretical counterparts CrXtl3 calculated with Equation (36). The phase angle β equals 45°, that is, the rotor has a displacement along the direction of y = x. It can be observed that the deviation between the value CrXf and the CrXtl3 is closer to the actual value.
In order to theoretically analyze the tilt error motion during the measurement, we considered tiny tilt displacement of a precision spindle and reasonable structural parameters of the sensor, that is, |/a| < 1 for the parameter term /a in Equation (15). The power series expansion with first-order approximation (Equation (38)) was also applied to express the integrand term in Equation (15):
IGe tl = ρ a [ 1 + k a ρ cos ( φ + γ ) ]
Correspondingly, the expression used to assess the approximate error is given as:
IGe a p p . e r r o r = | IGe IGe tl IGe |
By considering the relative tiny tilt displacement of a precision spindle and reasonability of the sensor structural parameters, the variation range of the term /a was set to be [0, 0.2]. As shown in Figure 12, the approximate error IGeapp.error is relatively smaller over the range of the parameter terms cos(φ + γ) and /a. According to the simulation model parameters of the sensor [16], the ratio of the term /a is about 0.04 with the tilt displacement of 200 arc-sec; thus, the error IGeapp.error remains no more than 0.2% within the whole range of the term cos(φ + γ), i.e., the entire phase region. Figure 13a compares the theoretical and simulated capacitance of the end part electrode. The theoretical capacitance is integrated by the integrand term with a second-order approximation (C8tl2). It can be observed that there is little difference between C8tl2 and the counterpart using first-order approximation (C8tl), which is also the same for the theoretical value of the CeY (the differential output capacitance of the EPEG) shown in Figure 13b. Besides, similar to the theoretical value of the CrX calculated by the expression with third-order approximation (see in Figure 11), the difference between the theoretical capacitance CeYtl and the simulation one CeYf remains roughly unchanged relative to the simulation value, which indicates that the theoretical capacitance of the end part electrode solved by adopting IGetl as integrand can meet the design needs for relatively accurate results.

3.2. Influence of Fringe Effects

Due to the existence of diverging electric fields in the plate edges in Figure 14, i.e., the “fringe effects” [19], the measured capacitance of a sensor contains some additional capacitance introduced by the fringe effects. As for the sensor electrode with finite dimensions, the variation in its output capacitance is relatively small, especially for the measurement of tiny displacement. As such, the influence of the fringe effects cannot be neglected.
Figure 15 shows the distribution of the electric field between the cylindrical excitation electrode (Ed) and curved sensing electrodes (Es) for the fundamental configuration of the radial curved plate capacitor. It can be seen that the electric field concentration appears in the junction of the inner circle face of the curved plate and its side face. As shown in Figure 15a, a different intensity distribution of the electric field between the Ed and the side face of the Es and its adjacent outer circular face formed along the radial direction. Once the equipotential guard ring (Egr) is employed outside the Es, the electric field distribution between the Ed and the outer circular face of the Es and between the Ed and local region of the side face near the outer circular face can be eliminated, as shown in Figure 15b. The Egr is also classically named Kelvin guard-ring.
The variation in the electric field distribution is also reflected in the output capacitance of the Es, as illustrated in Figure 16. The use of the Egr can effectively reduce the additional capacitance introduced by the fringe effects (see Figure 16a,b). Moreover, the guard ring could decrease the variation in the additional capacitance with the reduction in radial spacing (see Figure 16d), which is helpful to improve the measurement accuracy and linearity. For the T-type CS, the Egr is distributed around the periphery of the sensing electrodes with a tiny gap (λ = 0.1 mm), forming a coplanar configuration between the Egr and the electrodes. Under the effect of the Egr, the amount of variation in the additional capacitance significantly decreased (about 10 times), and the corresponding nonlinear error reduced from 4.8% to 1.7%, as indicated in Figure 17.

3.3. Installation Errors of the Sensing Electrodes

The sensing electrodes of the T-type CS include the cylindrical and fan-shaped electrodes integrated into the REG and EPEG, respectively. These electrodes are fabricated by the FPCB process in one procedure, along with the equipotential guard ring, and the sensing electrodes have a fixed position relationship with the guard ring.
The fabricated REG with planar form (as shown in Figure 18a, solid line box) is installed in the annular groove of the stator cylindrical bore by the surface mounting method. The coaxiality error, polar angle position error and tilt error (the axial boundary of the electrode relative to the axial reference of the stator) relative to the mounting reference were produced during the process. By utilizing the electrode measurement and sensor calibration method, the coaxiality and polar angle position error can be modified and removed, respectively. Thus, the analysis of the tilt error is the main focus for the cylindrical electrodes.
By spreading out the installed REG, as shown in Figure 18a, it can be seen that a certain tilt angle εe is produced relative to the ideal position of the REG. The point Om is assumed to be the tilt point. The dash line box represents the position of the cylindrical electrodes, and the dash line segment represents the axial boundary of the REG. As mentioned above, the in-plane relationship between the cylindrical electrodes is fixed during the mounting process, and thus the tilt angle εe is the tilt error of the electrodes.
As shown in Figure 18b, by taking a micro element with arc length on the boundary D′C′ of the electrode Er4, the axial length wJ of the electrode at this location can be considered as the function of the polar angle θJ′ at point J. The solution region is vertically divided into three parts according to the constraint boundary of the length wJ, and the corresponding calculation function can be expressed as follows:
w I = k I R θ + b I
w I I = b I I
w I I I = k I I I R θ + b I I I
where kI = k1 + 1/k1, kIII = −kI, bI = b3b4, bII = b3b1, bIII = b2b1; k1 is the slope of the boundary D′C′; b1~b4 are the intercepts of the lines including D′C′, C′B′, B′A′, D′A′, respectively.
Referring to Equation (6), the approximate expression of the output capacitance of the electrode Er4 under the tilt angle εe can be obtained:
C 4 = ε 0 ε r R [ θ D θ A w I g α cos ( θ β ) d θ + θ B θ D w I I g α cos ( θ β ) d θ + θ C θ B w I I I g α cos ( θ β ) d θ ]
By considering the minor radial displacement of a precision spindle, at the condition |α/g| < 1, the power series expansion with neglection of the high-order term is applied to express the integrand term in Equation (43), and thus we have:
C 4 = ε 0 ε r R g { [ k I R θ 2 2 + b I θ k I R ( α g ) ( cos ( θ β ) + ( θ b I k I R ) sin ( θ β ) ) k I R ( α g ) 2 ( θ 2 4 b I k I R ( θ β ) 2 + cos ( 2 ( θ β ) ) 8 + ( θ b I k I R ) sin ( 2 ( θ β ) ) 4 ) ] | θ D θ A + [ 2 b I I g α g α g + α arctan ( g + α g α tan ( θ β 2 ) ) ] | θ B θ D + [ k I I I R θ 2 2 + b I I I θ k I I I R ( α g ) ( cos ( θ β ) + ( θ b I I I k I I I R ) sin ( θ β ) ) k I I I R ( α g ) 2 ( θ 2 4 b I I I k I I I R ( θ β ) 2 + cos ( 2 ( θ β ) ) 8 + ( θ b I I I k I I I R ) sin ( 2 ( θ β ) ) 4 ) ] | θ C θ B }
Similarly, the approximate expression of the output capacitance of the electrodes Er3~Er1 are derived as follows:
C 3 = ε 0 ε r R g { [ k 2 I R θ 2 2 + b 2 I θ k 2 I R ( α g ) ( cos ( θ β ) + ( θ b 2 I k 2 I R ) sin ( θ β ) ) k 2 I R ( α g ) 2 ( θ 2 4 b 2 I k 2 I R ( θ β ) 2 + cos ( 2 ( θ β ) ) 8 + ( θ b 2 I k 2 I R ) sin ( 2 ( θ β ) ) 4 ) ] | θ D 2 θ A 2 + [ 2 b 2 I I g α g α g + α arctan ( g + α g α tan ( θ β 2 ) ) ] | θ B 2 θ D 2 + [ k 2 I I I R θ 2 2 + b 2 I I I θ k 2 I I I R ( α g ) ( cos ( θ β ) + ( θ b 2 I I I k 2 I I I R ) sin ( θ β ) ) k 2 I I I R ( α g ) 2 ( θ 2 4 b 2 I I I k 2 I I I R ( θ β ) 2 + cos ( 2 ( θ β ) ) 8 + ( θ b 2 I I I k 2 I I I R ) sin ( 2 ( θ β ) ) 4 ) ] | θ C 2 θ B 2 }
C 2 = ε 0 ε r R g { [ k 3 I R θ 2 2 + b 3 I θ k 3 I R ( α g ) ( cos ( θ β ) + ( θ b 3 I k 3 I R ) sin ( θ β ) ) k 3 I R ( α g ) 2 ( θ 2 4 b 3 I k 3 I R ( θ β ) 2 + cos ( 2 ( θ β ) ) 8 + ( θ b 3 I k 3 I R ) sin ( 2 ( θ β ) ) 4 ) ] | θ D 3 θ A 3 + [ 2 b 3 I I g α g α g + α arctan ( g + α g α tan ( θ β 2 ) ) ] | θ B 3 θ D 3 + [ k 3 I I I R θ 2 2 + b 3 I I I θ k 3 I I I R ( α g ) ( cos ( θ β ) + ( θ b 3 I I I k 3 I I I R ) sin ( θ β ) ) k 3 I I I R ( α g ) 2 ( θ 2 4 b 3 I I I k 3 I I I R ( θ β ) 2 + cos ( 2 ( θ β ) ) 8 + ( θ b 3 I I I k 3 I I I R ) sin ( 2 ( θ β ) ) 4 ) ] | θ C 3 θ B 3 }
C 1 = ε 0 ε r R g { [ k 4 I R θ 2 2 + b 4 I θ k 4 I R ( α g ) ( cos ( θ β ) + ( θ b 4 I k 4 I R ) sin ( θ β ) ) k 4 I R ( α g ) 2 ( θ 2 4 b 4 I k 4 I R ( θ β ) 2 + cos ( 2 ( θ β ) ) 8 + ( θ b 4 I k 4 I R ) sin ( 2 ( θ β ) ) 4 ) ] | θ D 4 θ A 4 + [ 2 b 4 I I g α g α g + α arctan ( g + α g α tan ( θ β 2 ) ) ] | θ B 4 θ D 4 + [ k 4 I I I R θ 2 2 + b 4 I I I θ k 4 I I I R ( α g ) ( cos ( θ β ) + ( θ b 4 I I I k 4 I I I R ) sin ( θ β ) ) k 4 I I I R ( α g ) 2 ( θ 2 4 b 4 I I I k 4 I I I R ( θ β ) 2 + cos ( 2 ( θ β ) ) 8 + ( θ b 4 I I I k 4 I I I R ) sin ( 2 ( θ β ) ) 4 ) ] | θ C 4 θ B 4 }
where, k4I = k3I = k2I = kI, k4III = k3III = k2III = kIII, b4II = b3II = b2II = bII, biI = b3bi4, biIII = bi2b1(i = 4, 3, 2); bi2, bi4 are the intercepts of the lines, including Ci′Bi′, Di′Ai′ (i = 4, 3, 2), respectively; i represents the electrodes Er1~Er3 in turn.
Further, the influence of the tilt error εe on the output capacitance of the cylindrical electrode (taking Er4 as an example) can be reflected by:
C 4 t i l t . e r r o r = C 4 C 4 C 4
By referring to the structural design of the T-type CS, the range of the tilt angle εe was set to be [0.006, 0.06] deg (about 200 arc-sec) in the numerical simulation using Matlab. A tiny circumferential displacement of the cylindrical electrode was produced due to the tilt angle εe, which produces a relative error of the output capacitance under the radial displacement along the same direction. The dependence of the relative error on the angle εe is shown in Figure 19. For the single electrode, the influence of the tilt angle εe on the output capacitance is very small, with respect to the amplitude of the relative error.
Figure 20 presents the dependence of the differential capacitance on the tilt angle εe. It can be known that the influence of the tilt angle εe on the differential output capacitance of the REG is slightly larger. By considering the rotor of run-out in the given directions (β = 5°, 45°, 85°), the differential output capacitance (CrXt-εe and CrYt-εe) can be calculated by the capacitance of single electrodes under the angle εe. The difference between the differential output capacitance under the angle εe and the counterpart under an ideal position (CrXt and CrYt) could be transformed to the difference of displacement parameters, utilizing Equations (33) and (34). The difference in the displacement parameters (magnitude α and phase angle β) are shown in Figure 21. It can be known that the produced tilt error εe of the cylindrical electrode only generates a certain deviation of solving the phase angle β, which is consistent in different displacement directions (brings about the overall phase advance or lag of the displacement trajectory); the effect of tilt error εe on solving the radial displacement α can be neglected, especially for the measurement range of less than 0.1 mm.
The fabricated EPEG was also in planar form (as shown in Figure 22a, solid line box) and installed in the annular groove located at the end part of the stator by the surface mounting method. The coaxiality error, polar angle position error and parallelism error (the electrode plane relative to the axial reference of the stator) relative to the mounting reference were produced during the process. For the parallelism error, it can be modified by the electrode measurement. Thus, the analysis of the fan-shaped electrodes is emphasized on the coaxiality error and polar angle position error.
Due to the position error between the shape boundary of EPEG and the sensing electrodes and the manipulation precision, there may be an eccentricity between the geometric center Oe of the fan-shaped electrodes in the EPEG and the medial axis of the stator (the origin position of the coordinate system). As shown in Figure 22a, δe denotes the magnitude of the eccentricity, and φe is the phase angle. δe is the coaxiality error of the fan-shaped electrode.
By taking the electrode Ee1 as an example, the region confined by the points C′D′S′T′ (dash line) represents the position of Ee1 under the coaxiality error δe. As shown in Figure 22b, the micro-plane element ΔBP′ (oblique line filled) at any point P′ is taken as the solving unit, and the area of this element is denoted as ΔSP′. Then, the ΔSP′ is approximately given by:
Δ S P = d x d y
The projection point of P′ on the x′-axis is the point P′x′. By referring to Equation (13), the expression for calculating the spacing t′ at this point could be derived as:
t = a k ( x cos γ y sin γ )
Correspondingly, the capacitance of the plane parallel capacitor composed of the micro-planes ΔBP′ and ΔAQ′ is:
Δ C e = ε 0 ε r Δ S P t = ε 0 ε r d x d y a k ( x cos γ y sin γ )
Through the integration of Equation (51) over the area of the fan-shaped electrode, the output capacitance of the electrode Ee1 is approximately expressed as:
C 5 = D Δ C e = ε 0 ε r D d x d y a k ( x cos γ y sin γ )
where the integral domain D is a closed region confined by the boundary curves of the fan-shaped electrode, which can be divided into three sub-regions, i.e., DI, DII and DIII. These three sub-regions can be expressed as follows:
D I I I = { ( x , y ) | ρ 1 2 ( x x B ) 2 + y B y k C x + b C , D x x C x } , D I I = { ( x , y ) | ρ 1 2 ( x x B ) 2 + y B y ρ 2 2 ( x x B ) 2 + y B , C x x E x } , D I = { ( x , y ) | k S x + b S y ρ 2 2 ( x x B ) 2 + y B , E x x T x } .
where xB and yB are the coordinates of the geometric center Oe; D′x, C′x, E′x and T′x represent the abscissas of the corresponding points, respectively.
Equation (52) can be rewritten as:
C 5 = ε 0 ε r a D d x d y 1 k a ( x cos γ y sin γ )
As for multivariate function, it can be expanded to power series under certain conditions by analogy with the univariate power series expansion [20]. By considering the tiny tilt displacement of a precision spindle and reasonable structural parameters of the sensor, there is generally |(k/a)∙(xcos(γ) − ysin(γ))| << 1, the power series expansion with neglecting of a high-order term is applied to express the integrand term in Equation (53). Then, it can be expressed as:
C 5 = ε 0 ε r a D [ 1 + k a ( x cos γ y sin γ ) ] d x d y
By calculating Equation (54) within the integral sub-domains DI, DII and DIII, we have:
C 5 = ε 0 ε r a { D x C x d x ρ 1 2 ( x x B ) 2 + y B k C x + b C [ 1 + k a ( x cos γ y sin γ ) ] d y + C x E x d x ρ 1 2 ( x x B ) 2 + y B ρ 2 2 ( x x B ) 2 + y B [ 1 + k a ( x cos γ y sin γ ) ] d y + E x T x d x k S x + b S ρ 2 2 ( x x B ) 2 + y B [ 1 + k a ( x cos γ y sin γ ) ] d y } = ε 0 ε r a { [ ( k C + μ b C ν k C b C μ y B + ν x B ) x 2 2 + ( μ k C ν 2 k C 2 ν 2 ) x 3 3 ( 1 ν y B + μ x B ) ( x x B 2 ρ 1 2 ( x x B ) 2 + ρ 1 2 2 arcsin x x B ρ 1 ) + μ 3 ( ρ 1 2 ( x x B ) 2 ) 3 + ( ( b C y B ) ν 2 ( b C 2 ρ 1 2 ) ν 2 ( x B 2 y B 2 ) ) x ] | D x C x + [ ( 1 ν y B + μ x B ) ( x x B 2 ( ρ 2 2 ( x x B ) 2 ρ 1 2 ( x x B ) 2 ) + ρ 2 2 2 arcsin x x B ρ 2 ρ 1 2 2 arcsin x x B ρ 1 ) + μ 3 ( ( ρ 1 2 ( x x B ) 2 ) 3 ( ρ 2 2 ( x x B ) 2 ) 3 ) ν 2 ( ρ 2 2 ρ 1 2 ) x ] | C x E x + [ ( μ y B ν x B + ν k S b S μ b S k S ) x 2 2 + ( ν 2 + ν 2 k S 2 μ k S ) x 3 3 + ( 1 ν y B + μ x B ) ( x x B 2 ρ 2 2 ( x x B ) 2 + ρ 2 2 2 arcsin x x B ρ 2 ) μ 3 ( ρ 2 2 ( x x B ) 2 ) 3 + ( ( y B b S ) ν 2 ( ρ 2 2 b S 2 ) ν 2 ( y B 2 x B 2 ) ) x ] | E x T x }
Similarly, the approximate expressions of the output capacitance of the electrodes Ee2~Ee4 can be derived as follows:
C 6 = ε 0 ε r a { [ ( k C 1 + μ b C 1 ν k C 1 b C 1 μ y B + ν x B ) x 2 2 + ( μ k C 1 ν 2 k C 1 2 ν 2 ) x 3 3 ( 1 ν y B + μ x B ) ( x x B 2 ρ 1 2 ( x x B ) 2 + ρ 1 2 2 arcsin x x B ρ 1 ) + μ 3 ( ρ 1 2 ( x x B ) 2 ) 3 + ( ( b C 1 y B ) ν 2 ( b C 1 2 ρ 1 2 ) ν 2 ( x B 2 y B 2 ) ) x ] | C 1 x D 1 x + [ ( 1 ν y B + μ x B ) ( x x B 2 ( ρ 2 2 ( x x B ) 2 ρ 1 2 ( x x B ) 2 ) + ρ 2 2 2 arcsin x x B ρ 2 ρ 1 2 2 arcsin x x B ρ 1 ) + μ 3 ( ( ρ 1 2 ( x x B ) 2 ) 3 ( ρ 2 2 ( x x B ) 2 ) 3 ) ν 2 ( ρ 2 2 ρ 1 2 ) x ] | E 1 x C 1 x + [ ( μ y B ν x B + ν k S 1 b S 1 μ b S 1 k S 1 ) x 2 2 + ( ν 2 + ν 2 k S 1 2 μ k S 1 ) x 3 3 + ( 1 ν y B + μ x B ) ( x x B 2 ρ 2 2 ( x x B ) 2 + ρ 2 2 2 arcsin x x B ρ 2 ) μ 3 ( ρ 2 2 ( x x B ) 2 ) 3 + ( ( y B b S 1 ) ν 2 ( ρ 2 2 b S 1 2 ) ν 2 ( y B 2 x B 2 ) ) x ] | T 1 x E 1 x }
C 7 = ε 0 ε r a { [ ( μ y B ν x B + ν k C b C μ b C k C ) x 2 2 + ( ν 2 + ν 2 k C 2 μ k C ) x 3 3 ( 1 ν y B + μ x B ) ( x x B 2 ρ 1 2 ( x x B ) 2 + ρ 1 2 2 arcsin x x B ρ 1 ) + μ 3 ( ρ 1 2 ( x x B ) 2 ) 3 + ( ( y B b C ) ν 2 ( ρ 1 2 b C 2 ) ν 2 ( y B 2 x B 2 ) ) x ] | C 2 x D 2 x + [ ( 1 ν y B + μ x B ) ( x x B 2 ( ρ 2 2 ( x x B ) 2 ρ 1 2 ( x x B ) 2 ) + ρ 2 2 2 arcsin x x B ρ 2 ρ 1 2 2 arcsin x x B ρ 1 ) + μ 3 ( ( ρ 1 2 ( x x B ) 2 ) 3 ( ρ 2 2 ( x x B ) 2 ) 3 ) ν 2 ( ρ 1 2 ρ 2 2 ) x ] | E 2 x C 2 x + [ ( k S + μ b S ν k S b S μ y B + ν x B ) x 2 2 + ( μ k S ν 2 k S 2 ν 2 ) x 3 3 + ( 1 ν y B + μ x B ) ( x x B 2 ρ 2 2 ( x x B ) 2 + ρ 2 2 2 arcsin x x B ρ 2 ) μ 3 ( ρ 2 2 ( x x B ) 2 ) 3 + ( ( b S y B ) ν 2 ( b S 2 ρ 2 2 ) ν 2 ( x B 2 y B 2 ) ) x ] | T 2 x E 2 x }
C 8 = ε 0 ε r a { [ ( μ y B ν x B + ν k C 1 b C 1 μ b C 1 k C 1 ) x 2 2 + ( ν 2 + ν 2 k C 1 2 μ k C 1 ) x 3 3 ( 1 ν y B + μ x B ) ( x x B 2 ρ 1 2 ( x x B ) 2 + ρ 1 2 2 arcsin x x B ρ 1 ) + μ 3 ( ρ 1 2 ( x x B ) 2 ) 3 + ( ( y B b C 1 ) ν 2 ( ρ 1 2 b C 1 2 ) ν 2 ( y B 2 x B 2 ) ) x ] | D 3 x C 3 x + [ ( 1 ν y B + μ x B ) ( x x B 2 ( ρ 2 2 ( x x B ) 2 ρ 1 2 ( x x B ) 2 ) + ρ 2 2 2 arcsin x x B ρ 2 ρ 1 2 2 arcsin x x B ρ 1 ) + μ 3 ( ( ρ 1 2 ( x x B ) 2 ) 3 ( ρ 2 2 ( x x B ) 2 ) 3 ) ν 2 ( ρ 1 2 ρ 2 2 ) x ] | C 3 x E 3 x + [ ( k S 1 + μ b S 1 ν k S 1 b S 1 μ y B + ν x B ) x 2 2 + ( μ k S 1 ν 2 k S 1 2 ν 2 ) x 3 3 + ( 1 ν y B + μ x B ) ( x x B 2 ρ 2 2 ( x x B ) 2 + ρ 2 2 2 arcsin x x B ρ 2 ) μ 3 ( ρ 2 2 ( x x B ) 2 ) 3 + ( ( b S 1 y B ) ν 2 ( b S 1 2 ρ 2 2 ) ν 2 ( x B 2 y B 2 ) ) x ] | E 3 x T 3 x }
where D′ix, C′ix, E′ix and T′ix(i=1,2,3) represent the abscissas of boundary points of the electrodes in the 2nd to 4th quadrant, respectively; k C 1 = k C ,   k S 1 = k S ; μ = kcos(γ)/a, ν = ksin(γ)/a.
Further, the expression used to assess the influence of the coaxiality error on the output capacitance of the fan-shaped electrode (taking Ee1 as an example) is given as:
C 5 e c c . e r r o r = C 5 C 5 C 5
In the numerical simulation using Matlab, the range of the coaxiality error δe is set to be [−0.2, 0.2] mm, and the variation quantity of the phase angle φe is 8 degrees referring to the sensor structural design. Figure 23 shows the relative error variation in the output capacitance of fan-shaped electrodes due to the error δe under different phase angles φe. For each fan-shaped electrode, the relative error of output capacitance increases with the increase in the error δe. It should be pointed out that the magnitude of the relative error is quite small within the whole value range of the error δe. Moreover, the influence of the phase angle φe variation in the relative error is also very limited.
From Equations (27) and (28), it can be known that the tilt displacement ε and yaw angle γ are obtained based on the calculated differential output capacitance of the EPEG (CeY, CeX). As shown in Figure 24, under different yaw angle γ, the produced coaxiality error δe and phase angle φe do not cause the changing of the capacitance CeY and CeX, which indicates that the measurement of the rotor tilt displacement is basically not affected by the coaxiality error of fan-shaped electrode.
The polar angle position error of the fan-shaped electrodes equals each other and does not change the increment of the electrode polar angle. Thus, the effect of the polar angle position error on the yaw angle γ measurement can be eliminated by the sensor calibration method.

4. Conclusions

The main error sources for the error motion measurement of a precision shafting using a T-type capacitive sensor (CS) were investigated in this paper. Firstly, we analyzed the theoretical modeling error due to the approximate simplification for the capacitance expressions of each sensing electrode. The approximate error of the theoretical calculations for the radial error motion measurement could be reduced under small ratios α/(Rg) and α/g. The introduction of the third-order term in the approximate expressions of the CrX and CrY could significantly reduce the approximate error from 5.5% to 1.2% without sacrificing the measuring range and sensitivity. For the theoretical model for the tilt error motion measurement, the theoretical capacitance of the end part electrodes was solved by adopting IGetl (first-order approximation) as integrand can achieve relatively accurate results with an approximate error lower than 0.2%.
Subsequently, the 3D-FEA results indicate that for the T-type CS, the coplanar configuration of the sensing electrodes and Egr with a tiny gap (0.1 mm) could significantly decrease the amount of variation in the additional capacitance introduced by fringe effects and reduce the corresponding nonlinear error from 4.8% to 1.7%.
Finally, among the electrode installation errors, we emphatically analyzed the tilt error of the cylindrical electrode and coaxiality error of the fan-shaped electrode. Numerical simulation was carried out based on the measurement model with the installation errors. The simulation results show that the tilt error only generates a certain deviation of solving the phase angle β (consistent in different displacement directions). The effect of the tilt error on solving the radial displacement α can be neglected, especially for the measurement range of less than 0.1 mm. Additionally, the measurement of the rotor tilt displacement is basically not affected by the coaxiality error of the fan-shaped electrode.

Author Contributions

Conceptualization, W.W. and Z.C.; methodology, W.W. and K.X.; software, K.X.; data and results analysis, W.W.and K.X.; original draft preparation, K.X.; review and editing, W.W. and Z.C.; supervision, W.W. and Z.C.; Funding Acquisition, W.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (Grant No. U1709206, No. 51275465), Zhejiang Provincial Natural Science Foundation of China (Grant No. LZ16E050001) and the Key Laboratory of Advanced Manufacturing Technology of Zhejiang Province (Grant No. 2015KF02).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, G.; Chen, Y.; Lu, Q.; Wu, Q.; Wang, M. Multi-Physics Fields Based Nonlinear Dynamic Behavior Analysis of Air Bearing Motorized Spindle. Micromachines 2020, 11, 723. [Google Scholar] [CrossRef] [PubMed]
  2. IBAG. High Frequency Motor Spindles. Available online: www.ibag.ch (accessed on 20 April 2021).
  3. Nanotech. Ultra-precision Machining Systems. Our-Automotive-Revolution-Rev.0717.pdf. Available online: www.nanotechsys.com (accessed on 21 April 2021).
  4. Aerotech. ABRT400 Air-Bearing Direct-Drive Rotary Stage. Available online: www.aerotech.com/product/stages-actuators (accessed on 21 April 2021).
  5. Gerlach, B.; Ehinger, M.; Raue, H.K.; Seiler, R. Gimballing Magnetic Bearing Reaction Wheel with Digital Controller. In Proceedings of the 11th European Space Mechanisms and Tribology Symposium, Lucerne, Switzerland, 21–23 September 2005; pp. 35–40. [Google Scholar]
  6. MWI. Magnetic Bearing Momentum and Reaction Wheels with Internal Wheel Drive Electronics. Available online: www.electronicnote.com/RCG/MWI_A4.pdf (accessed on 22 April 2021).
  7. Smithanik, J.; Paul, Y. Applying API 617 to Expander-Compressors with Active Magnetic Bearings. In Proceedings of the 44th Turbomachinery Symposium. Turbomachinery Laboratories, Texas A&M Engineering Experiment Station, Houston, TX, USA, 14–17 September 2015. [Google Scholar]
  8. Marsh, E.; Grejda, R. Experiences with the Master Axis Method for Measuring Spindle Error Motions. Precis. Eng. 2000, 24, 50–57. [Google Scholar] [CrossRef]
  9. Baxter, L.K. Capacitive Sensors: Design and Applications; IEEE Press: New York, NY, USA, 1997. [Google Scholar]
  10. Nouira, H.; Vissiere, A.; Damak, M.; David, J.M. Investigation of the Influence of the Main Error Sources on the Capacitive Displacement Measurements with Cylindrical Artefacts. Precis. Eng. 2013, 37, 721–737. [Google Scholar] [CrossRef]
  11. Lin, C.; Zheng, S.; Li, P.; Shen, Z.; Wang, S. Positioning Error Analysis and Control of a Piezo-Driven 6-DOF Micro-Positioner. Micromachines 2019, 10, 542. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Bukhari, S.A.R.; Saleem, M.M.; Khan, U.S.; Hamza, A.; Iqbal, J.; Shakoor, R.I. Microfabrication Process-Driven Design, FEM Analysis and System Modeling of 3-DoF Drive Mode and 2-DoF Sense Mode Thermally Stable Non-Resonant MEMS Gyroscope. Micromachines 2020, 11, 862. [Google Scholar] [CrossRef] [PubMed]
  13. Chapman, P.D. A Capacitance Based Ultra-Precision Spindle Error Analyser. Precis. Eng. 1985, 7, 129–137. [Google Scholar] [CrossRef]
  14. Ahn, H.J.; Han, D.C.; Hwang, I.S. A Built-In Bearing Sensor to Measure the Shaft Motion of a Small Rotary Compressor for Air Conditioning. Tribol. Int. 2003, 36, 561–572. [Google Scholar] [CrossRef]
  15. Kim, J.H.; Chang, H.K.; Han, D.C.; Jang, D.Y.; Oh, S.I. Cutting Force Estimation by Measuring Spindle Displacement in Milling Process. CIRP Ann.-Manuf. Technol. 2005, 54, 67–70. [Google Scholar] [CrossRef]
  16. Xiang, K.; Wang, W.; Qiu, R.; Mei, D.; Chen, Z. A T-Type Capacitive Sensor Capable of Measuring 5-DOF Error Motions of Precision Spindles. Sensors 2017, 17, 1975. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Ahn, H.J.; Jeon, S.; Han, D.C. Error Analysis of the Cylindrical Capacitive Sensor for Active Magnetic Bearing Spindles. J. Dyn. Sys. Meas. Control. 2000, 122, 102–107. [Google Scholar] [CrossRef]
  18. Trusov, A.A.; Shkel, A.M. Capacitive Detection in Resonant MEMS with Arbitrary Amplitude of Motion. J. Micromech. Microeng. 2007, 17, 1583. [Google Scholar] [CrossRef]
  19. Lu, Q.L. Engineering Mathematics: Complex Function; Higher Education Press: Beijing, China, 1996. [Google Scholar]
  20. Gao, C. On Power Series of Functions in Several Variables. Coll. Maths. 2007, 3, 125–129. [Google Scholar]
Figure 1. Schematic of a T-type CS. The error motions δx and δy are detected using the REG, which is in the annular groove of the stator cylindrical bore; δz and εx and εy are detected using the EPEG, which is in the annular groove of the stator end part.
Figure 1. Schematic of a T-type CS. The error motions δx and δy are detected using the REG, which is in the annular groove of the stator cylindrical bore; δz and εx and εy are detected using the EPEG, which is in the annular groove of the stator end part.
Micromachines 13 00221 g001
Figure 2. Schematic of the stator: (a) The exploded view of the structure; (b,c) partial cross-sectional view of the stator. η is the thickness of the electrode, and λ is the gap between the electrode and guard ring.
Figure 2. Schematic of the stator: (a) The exploded view of the structure; (b,c) partial cross-sectional view of the stator. η is the thickness of the electrode, and λ is the gap between the electrode and guard ring.
Micromachines 13 00221 g002
Figure 3. The general sketch of the main research content.
Figure 3. The general sketch of the main research content.
Micromachines 13 00221 g003
Figure 4. The calculation model of the radial error motions: (a) Boundary before transformation; (b) boundary after transformation.
Figure 4. The calculation model of the radial error motions: (a) Boundary before transformation; (b) boundary after transformation.
Micromachines 13 00221 g004
Figure 5. The calculation model of the axial error motion.
Figure 5. The calculation model of the axial error motion.
Micromachines 13 00221 g005
Figure 6. The calculation model of the tilt error motions.
Figure 6. The calculation model of the tilt error motions.
Micromachines 13 00221 g006
Figure 7. Error distribution of the approximate calculation of the changing Δu in the complex plane t.
Figure 7. Error distribution of the approximate calculation of the changing Δu in the complex plane t.
Micromachines 13 00221 g007
Figure 8. Calculation error distribution of CrX (expansion with first-order approximation).
Figure 8. Calculation error distribution of CrX (expansion with first-order approximation).
Micromachines 13 00221 g008
Figure 9. Simulation values and theoretical one calculated by Equation (33) comparison of radial differential output capacitance.
Figure 9. Simulation values and theoretical one calculated by Equation (33) comparison of radial differential output capacitance.
Micromachines 13 00221 g009
Figure 10. Calculation error distribution of CrX (expansion with third-order approximation).
Figure 10. Calculation error distribution of CrX (expansion with third-order approximation).
Micromachines 13 00221 g010
Figure 11. Simulation values and theoretical one calculated by Equation (36) comparison of radial differential output capacitance.
Figure 11. Simulation values and theoretical one calculated by Equation (36) comparison of radial differential output capacitance.
Micromachines 13 00221 g011
Figure 12. Error distribution of the first-order approximation of the integrand term in the capacitance expression of the end part electrode.
Figure 12. Error distribution of the first-order approximation of the integrand term in the capacitance expression of the end part electrode.
Micromachines 13 00221 g012
Figure 13. Comparison of the simulation values of the end part electrode capacitance and the theoretical one calculated by the integrand term with first-order and second-order approximation respectively: (a) Single electrode capacitance; (b) Differential output capacitance.
Figure 13. Comparison of the simulation values of the end part electrode capacitance and the theoretical one calculated by the integrand term with first-order and second-order approximation respectively: (a) Single electrode capacitance; (b) Differential output capacitance.
Micromachines 13 00221 g013
Figure 14. The fringe field distribution of a parallel plate capacitor with the gap 2π (the field lines between the two plates are symmetric, and thus only consider the distribution above the centerline).
Figure 14. The fringe field distribution of a parallel plate capacitor with the gap 2π (the field lines between the two plates are symmetric, and thus only consider the distribution above the centerline).
Micromachines 13 00221 g014
Figure 15. Electric field distribution between the electrode plates of the curved plate capacitor: (a) Without the Egr; (b) using the Egr.
Figure 15. Electric field distribution between the electrode plates of the curved plate capacitor: (a) Without the Egr; (b) using the Egr.
Micromachines 13 00221 g015
Figure 16. Simulation and theoretical value comparison for the output capacitance of the curved sensing electrodes Es: (a,b) The 1st, 2nd quadrant; (c,d) variation tendency of the differential output capacitance of the electrode and their relative difference value. The C*fw and C*f in the curve diagram denote the simulation values under without and with the Egr, respectively; the C*t is theoretical calculated value. The dx is defined as a component of the displacement α in the x-axis direction.
Figure 16. Simulation and theoretical value comparison for the output capacitance of the curved sensing electrodes Es: (a,b) The 1st, 2nd quadrant; (c,d) variation tendency of the differential output capacitance of the electrode and their relative difference value. The C*fw and C*f in the curve diagram denote the simulation values under without and with the Egr, respectively; the C*t is theoretical calculated value. The dx is defined as a component of the displacement α in the x-axis direction.
Micromachines 13 00221 g016
Figure 17. Simulation and theoretical values comparison of radial differential output capacitance of the T-type CS.
Figure 17. Simulation and theoretical values comparison of radial differential output capacitance of the T-type CS.
Micromachines 13 00221 g017
Figure 18. Tilt error analysis of the radial cylindrical electrode: (a) Local view of the REG structural model and schematic of the tilt error of cylindrical electrode (planar status); (b) geometric analysis model of single electrode. The z-axis is parallel to the medial axis of the stator, and the x-axis is parallel to axial reference of it.
Figure 18. Tilt error analysis of the radial cylindrical electrode: (a) Local view of the REG structural model and schematic of the tilt error of cylindrical electrode (planar status); (b) geometric analysis model of single electrode. The z-axis is parallel to the medial axis of the stator, and the x-axis is parallel to axial reference of it.
Micromachines 13 00221 g018aMicromachines 13 00221 g018b
Figure 19. Relative error variations in the output capacitance of cylindrical electrodes under the different tilt angle εe: (a,c) The 1st quadrant; (b,d) the 4th quadrant.
Figure 19. Relative error variations in the output capacitance of cylindrical electrodes under the different tilt angle εe: (a,c) The 1st quadrant; (b,d) the 4th quadrant.
Micromachines 13 00221 g019aMicromachines 13 00221 g019b
Figure 20. The changing of the radial differential output capacitance under the different tilt angle εe: (a,c) The X-direction; (b,d) the Y-direction.
Figure 20. The changing of the radial differential output capacitance under the different tilt angle εe: (a,c) The X-direction; (b,d) the Y-direction.
Micromachines 13 00221 g020aMicromachines 13 00221 g020b
Figure 21. Comparison of the solved displacement parameters of the rotor under the different electrode poses: (a) Displacement α; (b) phase angle β.
Figure 21. Comparison of the solved displacement parameters of the rotor under the different electrode poses: (a) Displacement α; (b) phase angle β.
Micromachines 13 00221 g021
Figure 22. Coaxiality error analysis of the end part fan-shaped electrode: (a) Local view of the EPEG structural model and schematic of the coaxiality error of fan-shaped electrode (planar status); (b) geometric analysis model of single electrode.
Figure 22. Coaxiality error analysis of the end part fan-shaped electrode: (a) Local view of the EPEG structural model and schematic of the coaxiality error of fan-shaped electrode (planar status); (b) geometric analysis model of single electrode.
Micromachines 13 00221 g022
Figure 23. Relative error variations in the output capacitance of fan-shaped electrodes under the different coaxiality error δe: (a,c) The 1st quadrant; (b,d) the 2nd quadrant.
Figure 23. Relative error variations in the output capacitance of fan-shaped electrodes under the different coaxiality error δe: (a,c) The 1st quadrant; (b,d) the 2nd quadrant.
Micromachines 13 00221 g023
Figure 24. The changing of the end part differential output capacitance under the different coaxiality error δe: (a,c) The component around the y-axis; (b,d) the component around the x-axis.
Figure 24. The changing of the end part differential output capacitance under the different coaxiality error δe: (a,c) The component around the y-axis; (b,d) the component around the x-axis.
Micromachines 13 00221 g024aMicromachines 13 00221 g024b
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xiang, K.; Wang, W.; Chen, Z. Analysis of Main Error Sources for the Error Motion Measurement of a Precision Shafting Using a T-Type Capacitive Sensor. Micromachines 2022, 13, 221. https://doi.org/10.3390/mi13020221

AMA Style

Xiang K, Wang W, Chen Z. Analysis of Main Error Sources for the Error Motion Measurement of a Precision Shafting Using a T-Type Capacitive Sensor. Micromachines. 2022; 13(2):221. https://doi.org/10.3390/mi13020221

Chicago/Turabian Style

Xiang, Kui, Wen Wang, and Zichen Chen. 2022. "Analysis of Main Error Sources for the Error Motion Measurement of a Precision Shafting Using a T-Type Capacitive Sensor" Micromachines 13, no. 2: 221. https://doi.org/10.3390/mi13020221

APA Style

Xiang, K., Wang, W., & Chen, Z. (2022). Analysis of Main Error Sources for the Error Motion Measurement of a Precision Shafting Using a T-Type Capacitive Sensor. Micromachines, 13(2), 221. https://doi.org/10.3390/mi13020221

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