Figure 1 shows a schematic of the torsional thrust balance. An aluminum frame was fixed to a platform. The pendulum arm was attached to the frame with a flexural pivot, forming a torsional pendulum configuration with rotational axis defined by the pivot. The electrostatic combs were used to generate calibration force. The deflection displacement of the arm was measured by fiber optic displacement sensor at a certain distance from the pivot axis. To damp the oscillation of the arm to quickly stabilize to steady state deflection, an eddy current damper was used. The stopper was used to protect the displacement sensor probe from impact caused by accidental deflection of the arm.
2.1. Methodology for the Calculation of the Thrust
The dynamics of the torsional balance can be described by the harmonic oscillator equation:
where
θ is the deflection angle of the arm,
I is the moment of inertia about the rotation axis,
Ft(
t) is the force produced by the thruster,
Lt is the distance between thruster mount position and rotational axis,
c is the damping constant,
k is effective torsional spring rate,
ζ and
ωn are the damping ratio and natural frequency, respectively.
Assume that the thrust generated by the thruster is constant,
Ft(
t) =
T, the solution of Equation (1) is:
When the arm reaches a steady state in a short period of time, the new position can be found by taking
into Equation (4):
This equation shows that the thrust
T can be determined by measuring
θ if
k and
Lt are known. However, it is not convenient to directly measure the arm deflection angle
θ. For small deflections, the small angle approximation is valid and the arm deflection
θ can be determined by measuring the linear displacement at a distance
Ls away from the rotational axis:
The linear displacement can be measured more accurately than the deflection angle; for the instance of
θ < 4°, the arm position can be determined within 0.1% [
5]. The effective torsional spring rate
k is determined by the flexural pivot and thruster power supply wire. Due to the hysteresis of pivot and the wire routing method,
k is not constant and may vary with the arm deflection. Therefore, Equation (5) can be rewritten as:
We can see that determining the thrust comes from function F(x) and the measured displacement at steady state position x. The function f(x) can be fitted by exerting a sequence of known force Fe at a certain distance away from the rotational axis and monitoring the corresponding displacement of the arm.
2.2. Torsional Balance Design
To reduce the longitudinal dimension of the torsional balance, an asymmetric arm structure about the rotational axis was adopted. The arm was made of hard aluminum alloy to ensure sufficient structural strength and small moment of inertia. The thruster to be measured will be mounted on the longer end of the arm. The center of the thruster mounting table was designed to be 270 mm away from the rotational axis, and the ratio to the distance between the electrostatic comb attached position and the rotational axis was 2.5. To reduce the influence of the deviation of the center of mass of the entire arm from the rotational axis [
5], a counterweight structure that can be orthogonally adjusted in the plane perpendicular to the rotational axis was designed on the arm’s short end. By adjusting the counterweight, the arm can be balanced and make the mass center of the arm as close to the rotational axis as possible.
To make the torsional balance more compact and easier to assemble, only one C-Flex model F10 flexural pivot (torsional spring rate is 42.08 N·μm/rad) was employed to provide the torsional restoring force. It was chosen mainly because of its good compatibility with vacuum, its frictionless rotation, and its low hysteresis. The shift of geometric center of the pivot will be minimal on the order of 0.02% of its diameter for small rotation angles.
2.4. Electrostatic Comb Design and Calibration
The electrostatic force devices are preferred for calibration because it can generate a very small force in a repeatable and accurate manner both at atmosphere and in vacuum. The use of an electrostatic force requires no direct contact with the torsional balance, and no cabling is required on the stand itself, which makes electrostatic calibration particularly well suited for μN thrust balance [
8]. There are two main types of electrostatic devices, namely, electrostatic plates and electrostatic combs. However, the force produced by the electrostatic plates is inversely proportional to the square of the gap between the plates, and even a small change in the gap can lead to apparent force variation [
9]. In the torsional balance application, the gap will change during the calibration as the arm deflects; this will result in a relatively large calibration error. However, this gap dependance can be eliminated by using the electrostatic combs as shown in
Figure 2. The force produced by the engaged combs is given by [
10].
where
Nc is the number of comb pairs,
ε0 is the permittivity, 2
x0 is the engagement distance, and other structural parameters are defined and shown in
Figure 2. It can be found from Equation (8) that the electrostatic force is proportional to the engagement distance, which is beneficial for torsional balance calibration as the error resulting by the gap change is relatively small. The electrostatic combs were used in this work with
c =
d = 1 mm,
g = 0.5 mm,
Nc = 5, and
l = 14 mm.
The combs were made of aluminum and encased in Polyetheretherketones base. One comb was attached to the arm and served as the grounded part, while its counterpart was mounted on a 3-axis manual linear stage and supplied with a voltage. The grounded comb was also grounded to the vacuum facility through the torsional balance, which ensured there was no wire attached to any part influencing the balance characteristic. The engagement distance and the alignment of the electrostatic combs were determined by the stage linear micrometer (resolution is 0.02 mm).
The electrostatic force produced by the electrostatic comb may deviate from the theoretical predicted by Equation (8). Therefore, experiment calibration must be performed prior to application. The calibration progress is illustrated in
Figure 3. The grounded comb was placed vertically on the analytical mass balance (Metttler Toledo XS205) with a grounding wire attached to the pallet. The resolution of the balance is 0.01 mg which is equal to 0.098 μN (Local gravitational acceleration is 9.7915 m/s
2). The variable potential comb assembly was positioned above it through a 3-axis manual linear stage, which can ensure the engagement distance is adjustable with an accuracy of 0.01 mm. The applied voltage was provided by a high voltage power supply (Matsusada ES-5R1.2) and monitored by a differential voltage probe (Tektronix P5205A). The balance was reset to zero to avoid shift prior to each measurement. The voltage and balance reading were recorded after the balance reading was stable for more than 15 s. Each calibration condition was repeated three times to minimize the random error. To minimize the vibration and air flow influencing the calibration results, the electrostatic combs experiment was performed on a gas floating platform and covered by an acrylic box.
The force produced by the electrostatic combs at different applied voltage with engagement distance varying is shown in
Figure 4a. It can be found that the force is almost constant in the engagement distance range of 2–11 mm, which is beneficial to reduce the calibration error caused by the electrostatic comb arrangement. We set the engagement distance as 6 mm, the electrostatic force varying with the applied voltage is shown in
Figure 4b, and the maximum electrostatic force is about 135 μN at 1000 V. The fitting function for the experimental data is given as:
2.5. Eddy Current Damper Design
To suppress electrochemical degradation of the ILET, the emitter voltage polarity will alternate with a period of several seconds. When the voltage polarity alternates from positive to negative, the thrust will decrease from a value to zero and then increase to another certain value, and vice versa. Therefore, the torsional thrust balance dynamic response characteristic is very important for ILET thrust measurement; the time desired for the balance settling to steady state should be less than the one-half voltage alternation period. A damping ratio of 0.4~0.8 generally gives a good response [
11]. In order to obtain the torsional balance damping ratio without additional damping, we excited the torsional balance by an impulse. The FODS output voltage as a function of time is shown in
Figure 5. Assuming
k is constant (42.08 N·μm/rad), the moment of inertia, damping ratio, and natural frequency were estimated by fitting a damped sinusoid in the form shown as Equation (4) to the data, resulting in:
I = 0.012 kg·m
2,
ωn = 1.955 rad/s, and
ζ = 5.526 × 10
−4, which is much smaller than the desired value range. According to Equation (2), the desired damping ratio can be conveniently achieved by increasing the damping constant. We substituted
I,
k, and the desired damping ratio into Equation (2), resulting in the desired damping constant of 0.018 <
c < 0.036 N·s/m.
The damping can be achieved through an eddy current damper, which is contactless and well compatible with vacuum. The configuration of the eddy current damper used in this work is depicted in
Figure 6; it consists of a cylindrical permanent magnet and a square copper plate fixed on the arm’s short end. When the plate moves in the magnetic field generated by the permanent magnet, the damping force (Lorentz force) will generate in the plate, making the arm rapidly settle to a stable deflection position. The damping force can be written as [
11]:
where
Lp is the distance between the plate center and the rotational axis of the arm,
is the angular velocity of the arm,
dp is the gap between the magnet and plate,
rc is the radius of the plate inscribed circle,
t and
σ are, respectively, the thickness and conductivity of the plate, and
Bz is the magnetic flux density in
z direction. The integral result of the above equation can be simplified as [
12]:
where
is the average magnetic flux density within the equivalent inscribed plate and
α is the correction coefficient containing the plate geometric factor. The damping constant generated by the eddy current damper is given as:
According to the calculation method proposed by Bae [
12],
α is ~0.5 for the square plate. However, the
ce calculated by the author of [
13] referring to Bae’s method is 2.6 times bigger than the experimental result. We believe this is mainly due to the fact that the correction coefficient does not take into account the skin effect. Therefore, the correction coefficient was corrected to 0.2 in this work. The copper plate (
σ = 5.7 × 10
7 S/m) used in this design is 45 mm × 45 mm and the thickness is 2 mm, and the
Lp is 245 mm. We substituted all the parameters to Equation (12) to determine
following that a field of about 0.09~0.13 T to cover the desired damping constant range is needed. The COMSOL Multiphysics software was used to calculate the magnetic field above the permanent magnet; after varying several parameters referring to the commercial off-the-shelf permanent magnet, a 37.6 mm diameter with 7 mm height NdFeB permanent with a remanence of 1.3 T was settled. The magnet was mounted on a manual linear stage; the magnetic field in the plate can be varied by means of adjusting
dp. The simulation and experimental results of
ce as a function of
dp are shown in
Figure 7, and the blue labels are the damping ratio calculated from the experimental data. The experiment was performed in the atmospheric environment and without the thruster mock. It was found that the simulation and experimental data follow the same trend, and the simulation results are a little higher than the experimental ones when
dp > 6 mm, and the opposite is true when
dp < 6 mm. The
dp is set to 6 mm in the torsional balance calibration experiment and the corresponding damping ratio is 0.56.