1. Introduction
High-resolution satellites and scientific experiments in space demand high environment quality; however, micro-vibrations in the space environment seriously affects the accuracy of these scientific activities [
1,
2,
3]. A maglev vibration isolation platform with six degrees of freedom (DOF), which has been successfully used in the field of low frequency vibration isolation, can significantly reduce the impact of micro-vibrations on experimental equipment and thereby increase the accuracy of scientific research [
4,
5,
6].
The maglev vibration isolation platform with six DOF is composed of a stator and a floater, as shown in
Figure 1. The stator is fixed on the spacecraft, where various disturbances exist that are caused by power generation systems, maintenance systems, thermal control systems,
etc. The floater is used to fix payloads, such as ultra-precise equipment, sensitive loads, high-resolution telescopes
etc. [
7,
8]. In order to provide an acceptable acceleration level for the payload, sensing of the absolute acceleration of the vibration level is necessary. Appropriate control algorithms can be designed to produce the corresponding forces using maglev actuators to counteract the measured vibration [
9,
10]. As the acceleration magnitude of the disturbances is small in a micro-vibration environment, high measurement precision is especially important for the improvement of isolation performance.
Figure 1.
Schematic diagram of the active vibration isolation platform.
Figure 1.
Schematic diagram of the active vibration isolation platform.
The method used to measure angular acceleration by linear accelerometers was first proposed as early as 1965 [
11]. The angular acceleration measurement precision depends closely on the accelerometer configuration. Two typical configurations that use the minimum number of accelerometers to measure six degrees of freedom accelerations have been developed in past research, as shown in
Figure 2. A cubic configuration presented by Chen in 1994 [
12] is shown in
Figure 2a. Six accelerometers are fixed at the center of each cube’s six faces to detect six DOF motion, and the sensitive directions are set along the diagonal directions [
12,
13]. The method has been proven to have high accuracy and reliability in some studies [
14]. However, the cubic configuration is suitable only when the space available is cubic and it requires high installation precision. In addition, a planar triangle configuration of the accelerometers to detect six DOF accelerations has been applied to an active vibration isolation platform [
15]. As shown in
Figure 2b, six accelerometers are fixed on three vertices of an equilateral triangle. Three of the sensitive directions are vertical. Others are along the tangential direction of a circumcircle of the equilateral triangle. Even on the same plane, the planar triangle configuration needs three vertexes of an equilateral triangle. In an electrical system, the appropriate space is rare and the sensitive directions need high installation precision.
Figure 2.
Two typical accelerometer configurations. (a) Cubic configuration; (b) planar triangle configuration.
Figure 2.
Two typical accelerometer configurations. (a) Cubic configuration; (b) planar triangle configuration.
The aim of this work is to propose a measurement model that uses a minimum number of accelerometers but has high measurement precision for six DOF accelerations. In order to overcome the shortcomings of existing configurations mentioned above, such as space constraints, installation precision demands, and so on, a new accelerometer configuration for measuring six DOF accelerations was proposed in this paper. The conditions needed for minimizing angular acceleration calculation noise are presented. Moreover, a comparison was carried out between the proposed configuration and the planar triangle configuration. The effect of the installation position and orientation errors on the output values of accelerometers for these two configurations were analyzed by simulation. The calculation precision of angular acceleration was compared for the two configurations by practical experiments.
2. Analysis of Angular Acceleration Calculation
As shown in
Figure 3, O-XYZ is the inertial frame and o-xyz is the body frame. Point
m and point
n are selected arbitrarily on a rigid body. According to the dynamics theory of rigid bodies, the acceleration of any point on a rigid body contains inertial acceleration, centripetal acceleration and tangential acceleration.
Figure 3.
Schematic diagram of a rigid body.
Figure 3.
Schematic diagram of a rigid body.
The accelerations of point
m and point
n in the inertial frame can be expressed by Equations (1) and (2), respectively:
where
is the inertial acceleration of the coordinate origin;
is the angular acceleration vector relative to the body frame;
and
are the displacement vector from the coordinate origin to the point
m and
n, respectively;
is the angular velocity relative to the body frame.
The difference of the acceleration between
m and
n can be written as Equation (3):
where
is the displacement vector from point
m to point
n.
The matrix form of Equation (3) is written as follows:
where each item in the equations contains plus or minus. Equation (4) can be solved as:
As the magnitude of micro-vibrations is tiny, almost all the research about maglev vibration isolation platforms is only focused on the acceleration level [
3,
9,
16]. The squared terms and cross product terms of the angular velocity are all neglected to construct a dynamic model of the system with high precision [
17,
18]. This can be done in order to simplify theoretical analysis for the section.
Therefore, Equation (5) can be simplified as:
If an accelerometer
Ai is installed at the point
on the rigid body, and the sensitive direction cosine of the accelerometer is
, then the output value of the accelerometer can be written as:
2.1. Configuration for One DOF Angular Acceleration
Taking the floater as an example, three placements listed in
Figure 4 can be used to acquire one DOF angular acceleration of a rigid body. The sensitive directions of accelerometers
p and
q are assumed to be along the X axis. The distances between
p and
q along X, Y and Z directions are
d2,
d1, and
d3, respectively. Then, the angular acceleration could be calculated based on output values of the accelerometers.
Figure 4.
Configurations for detecting one DOF angular acceleration.
Figure 4.
Configurations for detecting one DOF angular acceleration.
As the X axis is the sensitive direction, only the acceleration along the X direction was considered. According to Equation (6), the difference of output values between accelerometers
p and
q in
Figure 4a,b can be deduced as:
According to Equation (8), angular acceleration around the Z axis in
Figure 4a,b can be obtained:
Based on
Figure 4c, the difference of output values between accelerometers
p and
q is:
Thus, the angular acceleration around the Z axis in
Figure 4c is:
Equation (11) shows that the angular acceleration around the Y axis must be known first to calculate the angular acceleration around the Z axis. The distance
d3 will bring more noise to the angular acceleration. Therefore,
Figure 4c is inadvisable, so it can be seen that a coplanar accelerometer configuration has higher measurement precision than a non-coplanar placement and the accelerometers can be placed on any point of a plane to obtain one DOF angular acceleration with high precision.
2.2. Configuration for Two DOF Angular Acceleration
The above analysis indicates that four accelerometers can measure two DOF angular accelerations. To save space and decrease cost, three accelerometers are enough to detect two DOF angular accelerations by sharing an accelerometer with another two accelerometers. As shown in
Figure 5, three accelerometers
m,
n,
k are placed on the same plane and the sensitive directions are all along the Z axis.
Figure 5.
Configuration for detecting two DOF angular acceleration.
Figure 5.
Configuration for detecting two DOF angular acceleration.
According to Equation (6):
Solving Equations (12) and (13), the angular accelerations around the X axis and Y axis are:
2.3. Configuration for Six DOF Angular Acceleration
Based on the previous analysis, at least six accelerometers are needed to measure six DOF accelerations of a platform. As shown in
Figure 6a, by adopting six single-axis accelerometers
m,
n,
k,
p,
q and
s, three linear accelerations and three angular accelerations of the floater can be calculated.
Figure 6.
Configuration of accelerometers for a floater.
Figure 6.
Configuration of accelerometers for a floater.
In general, multi-axis integrated accelerometers are used to save space, which not only improves the installation accuracy, but also reduces cost. The floater is taken as an example. A three-axis integrated accelerometer, a two-axis integrated accelerometer and a single-axis accelerometer are adopted. The configuration of the accelerometers after assembly is displayed in
Figure 6b. Based on Equation (7), the output of each accelerometer can be calculated:
Assuming the distances between any two accelerometers space can be represented by
d because all sides are equal, then the linear and the angular accelerations of the floater can be solved and simplified as Equation (17):
Based on the theory of random variables, the root mean square (RMS) noise of the angular acceleration can be calculated as follows:
where
,
and
are the RMS noise of angular acceleration about X axis, Y axis and Z axis, respectively.
According to Equation (17), RMS noise of each angular acceleration can be calculated by Equation (19):
where
,
,
,
,
are the variances of noise outputs of each accelerometer.
According to Equation (19), maximizing the distance d as much as possible in the effective installation space is the best way to minimize the calculation noise of angular acceleration and improve measurement precision.
For the planar triangle configuration, the outputs of each accelerometer are assumed to be
a1~
a6 and the radius of the circumcircle is
R. Then linear acceleration and angular acceleration can be calculated by Equation (20) [
15]. Similarly, the noise of the angular acceleration decreases with increasing radius
R:
3. Simulations of the Effect of Installation Error on the Output Values of Accelerometers
Sensor inaccuracies, such as installation error, base bending, temperature transients, cable whip and so on, will deviate the output electric signal of the accelerometer from what it would have been, but inaccuracies such as base bending, temperature transients and cable whip, are hard to describe for the analysis. As the same accelerometers and the same installation plane are assumed for comparison in the same environment, the influence of base bending, temperature transients and cable whip can be neglected for the simulation analysis. Only the installation position error and orientation error are mainly considered for the simulation. The installation error is composed of the installation position error and sensitive direction error. Position installation error refers to a small displacement that deviates from the ideal installation position. The sensitive direction error is generally caused by a small rotation angle rotating around the ideal sensitive direction.
In this paper, the accelerometers are placed on a plane. The ideal installation position is:
The ideal sensitive direction is:
Assuming position installation errors are
,
and
(
i = 1, 2, 3), respectively, then the actual installation position is:
For the maglev vibration isolation platform, the order of magnitude of the angular error Δ
α is about 0.01°. With a three-axis accelerometer taken as a case, as shown in
Figure 7, the theoretical directions are the X-axis, Y-axis and Z-axis, respectively. The angular offset is assumed to be
as the blue shading in the figure. The direction of the angular error is defined as
, which can be any direction around the sensitive axis in the range of 0° to 360°. The actual sensitive directions are:
Figure 7.
Sensitive direction error for a three-axis accelerometer.
Figure 7.
Sensitive direction error for a three-axis accelerometer.
Therefore, the actual sensitive direction cosine of the proposed configuration yields:
Then, according to Equation (7), the linear acceleration and angular acceleration can be obtained.
To acquire the influence of the installation errors on the output value of each accelerometer, simulations were conducted. It is assumed that the linear accelerations and the angular accelerations of the floater are all sinusoidal motions with 1 Hz frequency and 1 micro-g amplitude. The distance
d is about 0.4 m for the platform. The position installation errors
,
and
(
i = 1, 2, 3) are all 0.0005 m, and the magnitude of the angular error Δ
α is 0.05°. According to Equation (7) please confirm here is Equation (7), comparison results of output values between the proposed configurations with and without installation error are shown in
Figure 8.
Figure 8.
Comparison of output results of the proposed configuration.
Figure 8.
Comparison of output results of the proposed configuration.
According to the planar triangle configuration, the comparison between the output values of the ideal installation and the actual installation under the same assumptions are displayed in
Figure 9.
Figure 9.
Comparison of output results of planar triangle configuration.
Figure 9.
Comparison of output results of planar triangle configuration.
From
Figure 8, it can be seen that actual output values of other accelerometers are almost the same as ideal output values. However,
Figure 9 shows that the differences of most output values between the ideal values and the actual values of the planar triangle configuration are significant, so the total influence of the installation error on the output values of the proposed configuration is smaller than that on the planar triangle configuration. That indicates that the proposed configuration has higher measurement precision under the condition of identical installation errors.
4. Experiments
Experiments are carried out to verify the advantages of the proposed configuration by comparing angular accelerations obtained from the proposed configuration and the planar triangle configuration. Experimental setups are shown in
Figure 10. Model 2422 accelerometers from Silicon were adopted to test the acceleration. The error is about 12 μg/(root Hz), the sensitivity is 2 V/g and its range is ±2 g. Only the reduction of the RMS error is considered and compared, as different vibrations exist in the environment during the experiments. Two accelerometer configurations were installed on the floater of the vibration isolation system and the output value of each accelerometer was acquired. First, the whole experimental setup is stationary. As there are many kinds of disturbances presented in the environment, the output values of the accelerometers contain environment noise and system noise. The RMS noises of angular accelerations of the two configurations are listed in
Table 1. The data shows that the RMS noise of the angular acceleration of the proposed configuration has a 14.89% reduction compared to that of the planar triangle configuration in the same environment.
Figure 10.
Experimental setups. (a) Signal measurement and processing system; (b) The proposed configuration; (c) The planar triangle configuration.
Figure 10.
Experimental setups. (a) Signal measurement and processing system; (b) The proposed configuration; (c) The planar triangle configuration.
Table 1.
Comparison of root mean square (RMS) noise of angular acceleration.
Table 1.
Comparison of root mean square (RMS) noise of angular acceleration.
RMS Noise | The Proposed Configuration | The Planar Triangle Configuration |
---|
σαx | 0.0117 | 0.0142 |
σαy | 0.0119 | 0.0152 |
σαz | 0.0118 | 0.0120 |
Noise reduction | 14.89% | 0 |
Then, experiments were carried out on a shaker. The linear accelerations were set as sinusoidal movement with 0.05 m/s
2 amplitude and 1 Hz frequency and the angular velocity was set as a sinusoidal movement with 0.005 m/s amplitude and 1 Hz frequency. In the same time interval, the differences between the calculated angular accelerations and the set value for the two configurations are shown in
Figure 11. During the same period, the RMS error of the difference of the angular acceleration between the measurements and the actual motion for the two configurations is shown in
Table 2. Comparing to the planar triangle configuration, the RMS error of the angular acceleration has a reduction of 19.73%.
Figure 11.
Comparison of measuring results of two configurations. (a) The proposed configuration; (b) The triangle planar configuration.
Figure 11.
Comparison of measuring results of two configurations. (a) The proposed configuration; (b) The triangle planar configuration.
Table 2.
Comparison of RMS error of angular acceleration differences.
Table 2.
Comparison of RMS error of angular acceleration differences.
RMS Error | The Proposed Configuration | The Triangle Planar Configuration |
---|
σαx | 0.0119 | 0.0148 |
σαy | 0.0116 | 0.0157 |
σαz | 0.0118 | 0.0133 |
Error reduction | 19.73% | 0 |
According to the prediction step and update step of the classical Kalman filter algorithm [
19], the measurement results after filtering seen in
Figure 12 can be obtained. The comparisons indicate that the angular acceleration formulae in the paper are correct, although there are errors due to the inaccurate positioning and orientation of the accelerometers, noises in the accelerometer outputs, and so on. Comparisons of RMS error between measurement results and filtered results are listed in
Table 3. The RMS error can reduce 70.61% after filtering. All experiments prove that the configuration proposed in this paper has higher measurement precision.
Figure 12.
Comparison of true acceleration, measurement results and filtered results.
Figure 12.
Comparison of true acceleration, measurement results and filtered results.
Table 3.
Comparison of RMS error between measurement results and filtered results.
Table 3.
Comparison of RMS error between measurement results and filtered results.
RMS Error | Measurement Results | Filtered Results |
---|
αx | 0.0117 | 0.0034 |
αy | 0.0119 | 0.0034 |
αz | 0.0118 | 0.0036 |
Error reduction | 0 | 70.61% |