Abstract
In order to provide an ultraquiet environment for spacecraft payload, a six-degree-of-freedom microvibration isolation device for satellite control moment gyro (CMG) is proposed in this paper. The dynamic characteristics of the microvibration isolation device are analyzed theoretically and experimentally. The dynamic equations of the microvibration suppression device are established by using the Newton–Euler method. The dynamic responses are numerically solved and the frequency-domain characteristics of the microvibration isolation device under base excitation are analyzed. The analytical results are first verified numerically, and the two results are in good accordance. The experimental apparatus is built, and the vibration isolation performance is investigated. The acceleration transfer function is measured and the influence of the excitation amplitude on the vibration isolation performance is performed. It is shown that the amplification factor at the vicinity of the resonance frequency is within 10 dB, and the vibration isolation performance is significant at higher frequencies. The vibration attenuation performance at the main frequency of the CMG (100 Hz) is more than 30 dB. The microvibration suppression device can effectively suppress the microvibration generated by CMG during orbital operation.
1. Introduction
Remote sensing or observation of the Earth is one of the most important tasks of satellite. In order to get an ultra-high-resolution image, many disturbances that affect the performance of the sensitive payload need to be controlled on board. The most annoying disturbance that degrades the performance of payload is thermal deflection [1, 2] of satellite structure and microvibration ejected by moving parts such as cryocooler [3–5], flywheel [6, 7], and control moment gyroscope (CMG) [8]. As a matter of fact, CMGs are widely used in fast maneuver satellites because of their significant angular momentum. Nevertheless, the high-speed rotors in CMGs also make them become one of the largest microvibration sources. In a recently published paper [9], the on-orbit microvibration measurement in a remote-sensing satellite shows that flywheel and CMG are the main vibration source.
Therefore, it is urgent to isolate the microvibration generated by flywheel and CMGs in order to provide quiet environment for sensitive payload. Not surprisingly, many microvibration suppression methods including vibration source control and payload isolation [10] are employed to achieve the purpose. For example, the famous Hubble Telescope used the viscous damper developed by Honeywell to isolate flywheel vibration [11]. In Chandra X-ray Observatory [12], the Stewart-type passive vibration isolation device for flywheel is also considered, and the vibration isolation performance of which at the launch phase is also designed and tested. Li et al. [13] conducted research on the microvibration induced by flywheel and the effect on space camera; the experimental results on ground show that the camera can work normally after a vibration isolator is applied to the flywheel. Cobb et al. [14] proposed a passive-active vibration isolation system to protect the optical payload from the unwanted microvibration on board. Spanos et al. [15] at JPL invented a flexible active vibration isolation and pointing system which was successfully applied in James Webb Space Telescope. Zhang et al. [16, 17] proposed a new method to achieve better image quality of optical payload. In their research, both the control strategy and the passive vibration source isolation are implemented. Li et al. [18] also considered a similar system and the difference of which laying on the vibration isolation system is active. Zhang et al. [19] focuse on the design of internal vibration isolation structures for the CMG to reduce the vibration generated by the imbalance of the rotor. A fully passive two-stage isolation system based on viscoelastic materials has been developed in order to minimize the CMG microvibrations by Kawak [20]. This paper proposed a new type of vibration isolation device to isolate the CMG vibration transmitted to satellite. The device not only can attenuate vibrations at the main frequency of the CMG but also can suppress resonance amplification factor to a considerable low level. Some useful design guidelines and concrete engineering based experimental results can enrich the vibration isolation field for CMGs.
This paper is organized as follows. Firstly, the dynamic model of the proposed vibration isolation device based on a relaxation-type damper is deduced using Euler–Newton method. The dynamic equation is then solved numerically. Secondly, the vibration isolation performances for all six degrees of freedom are analyzed and the coupling features in different directions are discussed. Thirdly, the experimental setup is built, and the vibration isolation performance is evaluated in three translational directions. Some valuable conclusions are drawn in the last section.
2. Dynamic Modeling
2.1. The Microvibration Isolation Device
The Stewart-type structure is used very popularly to get vibration isolation performance in six degrees of freedom [21]. The schematic of the proposed microvibration isolation device for single CMG in this paper is shown in Figure 1, which is also based on a Stewart-type structure. The platform includes upper platform to connect with the CMG and the lower platform to be installed on the satellite. In between, six relaxation-type isolators are configured in a cubic pattern, which provides proper stiffness and damping to isolate the vibration generated by the CMG transmitting to the satellite. Figure 2 shows a simplified representation of the vibration isolation device in which some useful coordinates are defined. The local coordinates on the upper and lower platforms are and , respectively. The global inertial coordinate system is . The hinges of the upper and lower platforms with the isolator are defined as and . is the inertial principal axis coordinate system of the CMG.


(a)

(b)
The following symbols are used to build the dynamic equation of the system. is the position vector in coordinate , is the position vector of the hinge where the upper platform and the isolator connected in coordinate , is the position vector of the hinge where the lower platform and the isolator connected in coordinate , is the length vector of each isolator between the upper and lower platforms in coordinate , is the position vector of the mass center of the CMG with the upper platform in coordinate , and are the angular velocity and angular acceleration of the upper platform, and are the angular velocity and angular acceleration of the lower platform, is the inertial matrix of the CMG with upper platform in coordinate , and and are the inertial matrix of the upper and lower legs in coordinates and . and are the position vectors of the upper and lower legs in coordinates and .
2.2. Modeling of the Strut and Upper Platform
Different from many researchers who consider the single strut as a Kevin model or an active actuator, the strut used in this paper is a relaxation-type damper. The strut model and the load applied to the strut are shown in Figure 3.

As this vibration isolation device is a Stewart-type platform, the strut is symmetric along the center axis of the upper and lower platforms. The position vector of the ith strut in coordinate can be expressed aswhere , , , and .
The unit vector for the ith strut is given as
The translational velocity and acceleration along the strut for the ith strut arewhere .
The acceleration of the mass center of the upper and lower legs can be expressed aswhere and .
Equations of moment of momentum for the upper and lower legs are given as
Equation (6) can be simplified to express the force applied to the upper platform from one single strut, which is
As shown in Figure 4, the force applied to the upper platform from one single strut is . The gravity of the upper platform and the CMG is , and the force and moment disturbance generated by the CMG are assumed as and .

By using Newton–Euler method, the dynamic equation for the upper platform is given aswhere is the absolute acceleration of the upper platform in coordinate , is the position vector of the mass center of the upper platform in coordinate , and is the inertia matrix of the upper platform in coordinate .
2.3. Response Solution
The force applied to the upper platform from each strut can be expressed aswhere is the deflection of the joint of the damping and the parallel connected spring in the ith strut.
Applying Laplace transform to equation (9), one can get the simplified equationwhere , , , and .
Combining equations (7)–(9), the dynamic equation of the vibration isolation device subjected to base excitation can be derived asw where , , , , , , and .
Equation (11) gives the dynamic responses of the upper platform when the base is excited with harmonic vibration.
3. Performance of the Microvibration Isolation Device
The microvibration isolation device is designed to isolate vibration ejected by CMG and the parameters of which are given in Table 1. The relaxation-type damper has been manufactured beforehand and the parameters of which are tested and recorded in Table 2. Meanwhile, the structure parameters of the microvibration isolation device are summarized in Table 3.
3.1. Natural Frequency Analysis
First, the natural frequencies of the microvibration isolation in six directions are analyzed. By simplifying equation (11) and neglecting damping, one can get the free vibration equation of the isolation device without damping as
In equation (12), is the equivalent stiffness matrix and is the equivalent mass matrix. The first six-order natural frequencies of the vibration isolation system can be obtained by solving equation (12). The modal analysis is also conducted by using finite element software and the model of which is shown in Figure 5. The results of the natural frequency obtained from the finite element software are compared with the theoretical results solved from equation (12). The comparison results are given in Table 4 and the two results are in good accordance. Moreover, in order to understand the natural characteristics of this system thoroughly, the corresponding mode shapes are also summarized in Figure 6.


(a)

(b)

(c)

(d)

(e)

(f)
From the results of natural frequency and the mode shape, one can conclude that the third-order mode and the fifth-order mode are the Z-direction translational and Z-axis twist pure mode, respectively. The other four modes are coupled mode in X-direction and Y-direction.
3.2. Relaxation-Type Damper Performance
As mentioned above, the flexible strut used here is a relaxation-type damper. The advantage of this kind of damper over a Kevin model damper is that it can control the resonance amplification without losing capacity of attenuating higher-frequency vibration if the parameters of the damper are chosen carefully. To clarify this point more clearly, the transmissibility of the relaxation-type damper used in this paper is plotted in Figure 7. For the sake of comparison convenience, the transmissibility curve of a system with a Kevin model is also plotted in the same figure. One can see that although the two dampers exhibit the same resonance amplification factor, the vibration attenuation rate at the main frequency of CMG (100 Hz) has more than10 dB differences, which is the advantage of the relaxation-type damping.

3.3. Coupling Analysis
The main reason that induced the coupling in different directions for the vibration isolation device considered here is that the mass center is out of the plane of the upper platform. Without loss of generality, the parameters and are chosen to illustrate the coupling characteristic of the device. The transmissibility of the two cases is plotted in Figures 8 and 9 for base excitation.


By observing the two figures, one can conclude that if the mass center of the CMG is in the plane of upper platform, there is no coupling for the device. When the mass center of the CMG is higher than the plane of upper platform, the X-direction translation, Y-direction rotation, Y-direction translation, and X-direction rotation are coupled. Two peaks in the transmissibility are observed. However, the movements in Z-direction translation and rotation are independent of the other direction in spite of the changing of the mass center of the CMG. This conclusion is very important when designing such kind of device to acquire better vibration isolation performance in certain direction.
3.4. Vibration Isolation Performance in Translational and Rotational Direction
To evaluate the vibration isolation performance of the device, equation (13) is numerically solved by the Runge-Kutta method in MATLAB software. The base is applied with an excitation amplitude of 30 mg according to satellite experiment and the frequency range is from 2 to 40 Hz. The absolute transmissibility is used to represent the vibration performance, which is defined as the ratio of acceleration of upper and lower platforms. The translational transmissibility and rotational transmissibility are plotted in Figures 10(a) and 10(b). Meanwhile, the resonance amplification factor and the vibration attenuation rate at the main frequency of the CMG are summarized in Table 5.

(a)

(b)
The transmissibilities in X-direction and Y-direction translation are basically the same because of the symmetrical feature of the device. In all six directions, the minimum natural frequency is 6.84 Hz (X-direction translational) and the maximum natural frequency is 28.65 Hz (Y-direction rotational). The vibration attenuation rate at 100 Hz is beyond 30 dB, which means that the vibration isolation efficiency is more than 94%. Moreover, the resonant amplification factor is below 3 times.
4. The Experiment
4.1. The Scheme
The experimental scheme is shown in Figure 11. It includes a microvibration excitation table, the vibration isolation device, accelerometers, LMS data acquisition system, and suspension system.

To cancel the gravity, a soft suspension system is built in Figure 12. The vibration table which simulates the excitation of sweep sine frequency applies the excitation to lower platform; both the accelerations on the vibration table and on the CMG are measured in order to calculate the absolute transmissibility. The translational transmissibilities in X-direction, Y-direction, and Z-direction are tested.

4.2. The Soft Suspension System
Four elastic strips are used to suspend the CMG. First, the deflection length and the stiffness of the strip are tested. The force-deflection test apparatus is shown in Figure 13, and the force-deflection curve is fitted by the LMS method, which is given in Figure 14. In the figure, one can calculate that the stiffness of the strip is about 0.245 N/mm, which yields a suspension frequency of about 0.3 Hz. The natural frequency of the suspension system is far lower than the first natural frequency of the vibration isolation system, which means that the suspension system can affect the experimental result a little. Moreover, in the vertical direction, the suspension length is more than 3 meters; therefore, the horizontal suspension frequency is also low enough.


4.3. Time-Domain Results
The input and output accelerations measured in the three directions are shown in Figures 15 to 17. The input signal of the vibration table is a 30 mg amplitude and 2–200 Hz sine sweep vibration. In the time domain in Figures 15 to 17, one can see clearly the amplification and attenuation of the vibration isolation device. In the Z-direction translational excitation, only one peak is observed, which means that there is no coupling in this direction. The amplification factor at the resonance frequency is less than 2 times.



4.4. Frequency-Domain Results
The transmissibilities in the three directions are shown in Figures 18 to 20. Both the experimental and theoretical results are plotted. Overall, the proposed vibration isolation device has excellent vibration attenuation performance. More than 20 dB attenuation is obtained at the main frequency of CMG, which is 100 Hz. There are slight deviations of the theoretical and experimental results in X-direction and Y-direction, though the resonance frequency is basically accurate. The reason for the deviation is possibly the coupling and the joint friction or assembly error. Nevertheless, the theoretical results can be used to predict such kind of vibration isolation device in engineering, since the trends are in good consistency and the deviations only occur at frequencies of small range.



In the Z-direction transmissibility, the theoretical and experimental results are in superb accordance. This is mainly because, as mentioned above, there is no coupling in this direction. In Figure 20, a transmissibility of a vibration isolation system with Kevin model is also plotted for comparison convenience. The curve is plotted by assuming a same resonance amplification factor for both systems. It can be seen from the figure that the device proposed in this paper can provide better vibration isolation performance at the main frequency of the CMG. The device with relaxation-type dampers may isolate 10 dB more vibration than that with a Kevin model according to the experimental results. Moreover, the amplification factor at resonance frequency of this device is below 10 dB, which is an excellent index compared with traditional vibration isolation systems on market.
4.5. Effect of the Excitation Amplitude
In fact, the model studied in this paper is a linear system. In other words, the amplitude of excitation is irrelevant to the vibration isolation performance. However, this feature is verified experimentally for the sake of completeness. The excitation amplitude is varied as 10 mg, 20 mg, and 30 mg in Z-direction and the transmissibilities are plotted in Figure 21 for different excitation amplitudes. From the figure, one can easily distinguish that the amplitude of the excitation has rare effect on the vibration isolation performance.

5. Conclusions
A microvibration suppression device is proposed based on a relaxation-type damper and cubic configured Stewart platform in this paper. Both theoretical and experimental investigations are carried out. Some conclusions can be drawn as follows:(1)The height between the mass center of the CMG and the upper platform of the device is the main factor that induces coupling effect; this parameter should be chosen carefully when designing this kind of device.(2)The proposed device can provide more than 30 dB vibration attenuation at the main frequency of the CMG and less than 10 dB amplification at resonance, which is better than that of a traditional linear one.(3)The excitation amplitude has rare effect on the vibration isolation performance and the device can function in a wide amplitude range.
Data Availability
The data are available upon request to the corresponding author.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors acknowledge the financial support from the National Natural Science Foundation of China under Grant no. 51875363.