Jie Tang^{1} , Dengqing Cao^{2} , Zhaohong Qin^{3} , Haibo Li^{4} , Dengshuo Chen^{5}
^{1, 2, 5}School of Astronautics, Harbin Institute of Technology, Harbin, 150001, P. R. China
^{3, 4}Beijing Institute of Structure and Environment Engineering, Beijing, 100076, P. R. China
^{2}Corresponding author
Journal of Vibroengineering, Vol. 20, Issue 2, 2018, p. 10351050.
https://doi.org/10.21595/jve.2017.18494
Received 15 April 2017; received in revised form 23 July 2017; accepted 9 August 2017; published 31 March 2018
Copyright © 2018 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
In the launching process, the spacecraft situated in an extreme dynamic environment may withstand by various dynamic loads such as noise loading in the fairing, motor excited vibration, shock of the separation devices, etc. To achieve a successful launching, the device called wholespacecraft vibration isolator is usually installed between the adapter and the spacecraft to prevent the later from being damage. A novel WSVI device, which composed of supporting leaf springs, voice coil motors (VCM) and actuator supports, is designed to suppress the structural vibration of the spacecraft in this paper. The novel WSVI device features small space footprint and light weight, and satisfy the design requirement of vibration isolation without changing of the Payload Adaptor Fitting structure. A dynamic model of the WSVI device is established to evaluate the performance of the system. The dynamic characteristics and responses subjected to external excitation are studied for the spacecraft installed with WSVI. The vibration isolation performance is analyzed after turning the VCM into passive dampers. The results show that the novel WSVI device, which can reduce the amplitude of vibration response of the spacecraft significantly, is valid for vibration suppression of the spacecraft.
Keywords: wholespacecraft vibration isolation, voice coil motors, dynamic model, vibration suppression.
During the launching process of the rocket, the spacecraft may suffer an extreme dynamic environment. The launching dynamic loads include low frequency dynamic loads such as liftoff, motor excitation, buffet, motor starts and shutoffs. Spacecraft are also subjected to shock loads due to staging, spacecraft separation, and fairing separation [1]. These dynamic loads may cause damage to the structure of the spacecraft, which affect the normal function of the devices carried by the spacecraft. The wholespacecraft vibration isolation (WSVI) technique using a vibration isolation device to connect the rocket and spacecraft, instead of traditional rigid connection by Payload Adaptor Fitting (PAF) structure. By using this technique, the dynamic loads transferred to the spacecraft can be isolated, and the launch environment will be improved significantly.
As an important technology affecting the development of space technology, wholespacecraft isolation has drawn attention of researchers for a long time. On the initiative of United States Air Force Research Laboratory (AFRL) and Space Vehicles Directorate, CSA Engineering Inc [25] began the research about wholespacecraft isolation and shock resistance test of the spacecraft. They have proposed the SoftRide UniFlex and SoftRide MultiFlex WSVI systems. In order to solve the problem on the dimensional limit of the SoftRide MultiFlex and further enhance the performance of the dampers, a SoftRide OmniFlex isolation system was developed by CSA Engineering Inc in 2008 [68]. And now the SoftRide OmniFlex platform has become the main solution to WSVI problem of CSA Engineering Inc. In the frame of the ESA GPS program European Aeronautic Defense & Space (EADS) [9, 10] has developed an active payload adaptor for Ariane 5. And the proposed system is a truss structure with 24 active struts in the load path using piezoceramics as actuator. The active payload adaptor is thought to replace the conventional passive payload adaptors, for instance, the conical Carbon Fiber Reinforced Plastics (CFRP) shell structures. In 2004, the AFRL cooperated with ATA Engineering Inc developed a MultiPayload Adapter (MPA) for the Peacekeeper Space Launch Vehicle [11]. ATA’s MPA consists of an annular flat plate that has top and bottom face sheets separated by radial ribs and closeout rings. Also, this vibration isolator has the advantage of small footprint and strong commonality but has barely lateral vibration isolation ability. Besides, Scholars like Griffin, Khorrami [12, 13], etc. also have proposed different WSVI ideas.
The WSVI schemes reported in the literature were mainly focused on the passive vibration isolation, the isolation performance of which may not be very effective in the low frequency region. A novel WSVI device, composed of supporting leaf springs, voice coil motors and actuator supports, is designed to realize active vibration isolation of the spacecraft. The proposed device is not only designed to have the feature of small space footprint and light weight, but also satisfy the design requirement of vibration isolation without changing of the PAF structure. The isolation system is elegant in its simplicity, which can be used universally for different type of spacecraft.
Compared with other types of actuators, electromagnetic actuator (also called VCM) has the advantages of compact structure, fast response, wide control frequency band, large actuating displacement, easy assembly, large ratio of output force between quality and volume of the actuator, etc. It is a new type of actuator integration of mechanical, electrical and magnetic components. What is more, only electric energy provided for its normal operation, which has a highly engineering application value in aerospace engineering where the electric power is widely used. The structure of the electromagnetic actuator adopted in this thesis is as shown in Fig. 1. The cutaway view depicted in Fig. 1 illustrates a unique configuration inside the VCM, where the soft iron and magnet are interval arranged along the central linear bearing to ensure a maximum performance of the magnet.
Fig. 1. Cutaway view of the electromagnetic actuator
In order to derivate the driving force expression of the VCM, Maxwell 16.0 is utilized to simulate and analyze the magnetic field of the actuator, and the magnetic field lines and intensity distribution are shown in Fig. 2. It is shown that magnet steel with a high magnetic permeability can constrain the distribution of the magnetic field line, prevent it from diverging. The special shape and configuration of the permanent magnet steel contribute to a regular magnetic field in the actuator, where magnetic field line mainly along the area I and II when it crosses the air gap. Ideally, make a hypothesis that all the magnetic field line pass through the air gap with a height of $h$ (the thickness of the soft iron), and perpendicular to the coils, as shown in Fig. 3.
Fig. 2. The magnetic field lines and intensity distribution of the actuator
a) Magnetic field lines
b) Magnetic intensity distribution
According to the Gauss's law for magnetism, the total magnetic flux through a closed surface $S$ is equal to zero. The equation is defined by:
Hence in area I and II we can get:
Then:
where ${B}_{r}$ is the remanence intensity of the permanent magnet, ${B}_{s1}$ and ${B}_{s2}$ are the magnetic flux density of area I and II, respectively.
Fig. 3. Magnetic induction line in the actuator
By the Ampere theorem, Electric conductor in magnetic field will be affected by force, which contribute to the constitutive equation of the VCM:
where $n$ is the number of the coils, $B$ is the remanence intensity, $L$ is the length of per coil, $I$ is the current of the coil. Hence, the driving force of the actuator can be obtained from Eq. (4). For convenient we define $\lambda =nBL$ (N/A) as the force constant of the actuator, which is the main parameter that reflecting the performance of actuator comprehensively.
PAFs are usually adopted to connect the launcher and the payload. The WSVI system (Fig. 4) is proposed to replace the PAF with passive or active adaptor so that the dynamic load transmit to the payload can be attenuated. The novel WSVI device proposed in this paper is mainly composed of supporting leaf springs, electromagnetic actuators, actuator supports, etc. The supporting leaf spring is a Utype connector, the upper limb of which is connected to the bottom of the satellite and the lower limb is connected to the top of the PAF. All the leaf springs provide the support stiffness for the isolation platform. The electromagnetic actuators installed in the interior of the PAF, the mover connect to cantilever end of the leaf spring, and the stator fixed on the support. The supports also fixed in the interior of the PAF by bolted connection, which are used for supporting of the electromagnetic actuators. The supporting leaf spring, electromagnetic actuator and support form an actuation unit. The whole WSVI platform composed of eight actuation units, and uniformly installed in circumference direction of the PAF. The assembled WSVI platform is as shown in Fig. 5.
Fig. 4. The WSVI system
Fig. 5. The proposed WSVI platform
The visual representation of the mathematical model of the WSVI system is shown in Fig. 6. The spacecraft is assumed to be a rigid body (with the consideration that the solar panels are folded during the launching process) supported from below by 8 suspension legs which connect the top of the PAF to the bottom of the spacecraft. A single suspension leg will be defined here as one leaf spring and one electromagnetic actuator. The global coordinate frame $O$ is set with the orthonormal ${e}_{1}$ and ${e}_{2}$ axes parallel to the bottom of the PAF and the ${e}_{3}$ axis normal to the bottom of the PAF.
Fig. 6. Mathematical model of the system
The origin of the body coordinate frame $U$ is located at the center of gravity (CG) of the spacecraft. The generalized displacement vector of the body coordinate frame $U$ in the global coordinate frame $O$ is defined as:
where ${\mathbf{x}}_{t}^{O}={\left[\begin{array}{lll}{x}_{x}& {y}_{x}& {z}_{x}\end{array}\right]}^{T}$ is the translational displacement vector correspond to the displacements about the ${u}_{1}$, ${u}_{2}$, ${u}_{3}$ axes. ${\mathbf{x}}_{r}^{O}={\left[\begin{array}{lll}{\varphi}_{x}& {\theta}_{y}& {\psi}_{z}\end{array}\right]}^{T}$ is the rotational displacement vector corresponding to the rotations about the ${u}_{1}$, ${u}_{2}$, ${u}_{3}$ axes.
The generalized displacement vector of the body coordinate frame $B$ in the global coordinate frame $O$ is defined as:
where ${\mathbf{w}}_{t}^{O}={\left[\begin{array}{lll}{x}_{w}& {y}_{w}& {z}_{w}\end{array}\right]}^{T}$ is the translational displacement vector correspond to the displacements about the ${b}_{1}$, ${b}_{2}$, ${b}_{3}$ axes, ${\mathbf{w}}_{r}^{O}={\left[\begin{array}{lll}{\varphi}_{w}& {\theta}_{w}& {\psi}_{w}\end{array}\right]}^{T}$ is the rotational displacement vector corresponding to the rotations about the ${b}_{1}$, ${b}_{2}$, ${b}_{3}$ axes.
Let ${K}_{j}$ be the $j$th leaf spring coordinate frame and ${A}_{j}$ be the $j$th actuator coordinate frame. Transformation between different coordinates can be described by transformation matrix, where ${\mathbf{T}}_{{K}_{j}}^{O}$ is the transformation matrix from the local fixed frame ${K}_{j}$ to the global frame $O$. Euler angles${\varphi}_{{K}_{j}}$, ${\theta}_{{K}_{j}}$, ${\psi}_{{K}_{j}}$ can be used to describe the transformation of the $j$th leaf spring, and they are corresponding to the rotations about the ${k}_{{1}_{j}}$, ${k}_{{2}_{j}}$,${k}_{{3}_{j}}$ axes, respectively. Let ${K}_{j}$ rotation about ${k}_{{3}_{j}}$${k}_{{1}_{j}}$${k}_{{2}_{j}}$ in turn, here, $\varphi ={\varphi}_{{K}_{j}}$, $\theta ={\theta}_{{K}_{j}}$ and $\psi ={\psi}_{{K}_{j}}$, then the elementary rotation matrices can be written as:
And then the expression of the transformation matrix is:
where $s=\mathrm{s}\mathrm{i}\mathrm{n}$, $c=\mathrm{c}\mathrm{o}\mathrm{s}$. ${\mathbf{T}}_{{K}_{j}}^{O}$ is an orthogonal matrix, thus ${\mathbf{T}}_{{K}_{j}}^{O}=({\mathbf{T}}_{{K}_{j}}^{O}{)}^{T}$. The transformation matrices ${\mathbf{T}}_{U}^{O}$ and ${\mathbf{T}}_{{A}_{j}}^{O}$ are similarly defined.
In fact, the rotational displacement of the spacecraft is very small, namely the rotational displacements of body frame with respect to the global frame are assumed to be small, besides the body frame $U$ and global frame $O$ share the same initial orientation, then the position vector of any point n fixed in the body frame $U$ can be expressed in the global frame $O$ as:
where ${\mathbf{p}}_{nU}^{U}={\left[\begin{array}{lll}{p}_{x}& {p}_{y}& {p}_{z}\end{array}\right]}^{T}$ is the position vector from bode frame $U$ to point $n$, and ${\stackrel{~}{\mathbf{P}}}_{nU}^{U}$ is the skewsymmetric matrix defined for vector cross product operation given as:
The translational displacement of any point $n$ fixed in the body frame $U$ can be expressed in the global frame $O$ as:
The dynamic equation of the system can be written as:
where ${\overline{\mathbf{M}}}_{}^{O}$ is the generalized mass matrix, and ${\overline{\mathbf{f}}}_{K}^{O}$, ${\overline{\mathbf{f}}}_{A}^{O}$, ${\overline{\mathbf{f}}}_{g}^{O}$ are the generalized global forces on the body frame due to the leaf springs, actuators, and gravity, respectively. The generalized mass matrix can be expressed in the global frame as:
where ${\stackrel{}{\mathbf{J}}}^{O}={\mathbf{T}}_{U}^{O}{\mathbf{J}}^{O}({\mathbf{T}}_{U}^{O}{)}^{T}$, and ${\mathbf{J}}^{O}={\mathbf{J}}_{c}^{U}+m({\mathbf{r}}_{OU}^{T}{\mathbf{r}}_{OU}{\mathbf{I}}_{{}_{3\times 3}}{\mathbf{r}}_{OU}{\mathbf{r}}_{OU}^{T})$, in which $m$ is the mass of the payload, ${\mathbf{J}}_{c}^{U}$ is the moment of inertia matrix that the spacecraft with respect to its barycenter, ${\mathbf{J}}^{O}$ is the transition of moment of inertia matrix, ${\stackrel{}{\mathbf{J}}}^{O}$ is the moment of inertia matrix that the spacecraft with respect to global frame $O$.
The generalized force vector ${\overline{\mathbf{f}}}_{g}^{O}$ due to gravity applied to the body frame $U$ can be written as:
The stiffness of $j$th leaf spring defined in the local frame can be written as:
where ${k}_{1}$, ${k}_{2}$, ${k}_{3}$ are the lateral and longitudinal stiffness of the leaf spring, along ${k}_{{1}_{j}}$, ${k}_{{2}_{j}}$, ${k}_{{3}_{j}}$, respectively, as shown in Fig. 7(a) and:
The generalized force vector due to all $n$ leaf springs applied to the body frame $U$ can be expressed as:
where ${\stackrel{}{\mathbf{f}}}_{{K}_{j}}^{O}={\left[\begin{array}{ll}{\left({\mathbf{f}}_{{K}_{j}}^{O}\right)}^{T}& {\left({\mathbf{t}}_{{K}_{j}}^{O}\right)}^{T}\end{array}\right]}^{T}$, and ${\mathbf{f}}_{{K}_{j}}^{O}$, ${\mathbf{t}}_{{K}_{j}}^{O}$ are the force and torque vectors applied to the body frame $U$ due to the $j$th leaf spring, respectively. Then, the force vector for$j$th leaf spring can be written as:
And the torque vector ${\mathbf{f}}_{{K}_{j}}^{O}$ is:
Hence, the generalized force vector ${\stackrel{}{\mathbf{f}}}_{K}^{O}$ can be written as:
where ${\stackrel{}{\mathbf{K}}}^{O}=\sum _{j}^{n}{\stackrel{}{\mathbf{K}}}_{j}^{O}$ is the generalized global stiffness matrix, and:
Fig. 7. Coordinate frame of leaf spring and actuator
a) Coordinate frame of $j$th leaf spring
b) Coordinate frame of $j$th actuator
In this WSVI scheme, the electromagnetic actuator vertically mounted on the support by universal spherical joints which allows only axial force to be transmitted through the actuator. Hence, the force of the $j$th electromagnetic actuator defined in its local coordinate frame ${A}_{j}$ (see Fig. 7(b)) can be expressed as:
where ${f}_{{a}_{j}}=\lambda {i}_{j}$, $\lambda $ is the force constant of the electromagnetic actuator, ${i}_{j}$ is the drive current in $j$th electromagnetic actuator. The generalized global force vector due to all $n$ actuators applied to the body frame $U$ can be described as:
where the generalized force vector due to $j$th electromagnetic actuator is expressed as:
Let $\mathbf{H}=\lambda \left[{\mathbf{H}}_{1}\begin{array}{ll}\cdots & {\mathbf{H}}_{n}\end{array}\right]$, and $\mathbf{i}=[{i}_{1}\begin{array}{ll}\cdots & {i}_{n}\end{array}{]}^{T}$. Then Eq. (25) can be written as:
Hence, Eq. (14) can be written as:
The equation of motion of the system can be expressed in state space form, namely:
where $\mathbf{x}={\left[\begin{array}{ll}({\stackrel{}{\mathbf{x}}}^{O}{)}^{T}& ({\dot{\stackrel{}{\mathbf{x}}}}^{O}{)}^{T}\end{array}\right]}^{T}$ is the state vector, $\mathbf{w}={\left[\begin{array}{ll}({\stackrel{}{\mathbf{w}}}^{O}{)}^{T}& 1\end{array}\right]}^{T}$ is the external disturbances, and $\mathbf{u}=\mathbf{i}$ is the system input. Then, the system matrix is:
the input matrix is:
and the disturbance matrix is:
Up to now, the kinematical equations of the WSVI system are described mathematically and the dynamic equations of the proposed WSVI system are established. Moreover, the corresponding statespace expression of the system is developed for vibration analysis and controller design in the future.
A scaleddown model of the spacecraft with WSVI device has been designed for further investigation. The physical parameters of the spacecraft model and the WSVI device are listed in Table 1.
Substituting the parameters listed in Table 1 into Eq. (28), the natural frequencies of the system are calculated and listed in Table 2, in which the first order frequency correspond to the rolling modal around the lateral axis, the second order frequency is the longitudinal vibration along the rocket axis, the third order frequency is the lateral vibration normal to the rocket axis, and the fourth frequency is the torsional vibration about the vertical axis.
Table 1. Physics parameters of the system
Physics parameter

Value

Mass of satellite

23.46 kg

Moment inertia of the satellite

Diag (0.94 0.94 1.02) kg·m^{2}

Distance between body frame $U$ to global frame $O$

0.788 m

Distance between body frame $U$ to the bottom of the satellite

0.456 m

Leaf spring stiffness

Diag (3.3e5 3e5 2e5) N/m

Force constant of the actuator

41 N/A

Gravitational acceleration

9.8 m/s²

Passive damping coefficient

10 N·s/m

Table 2. Natural frequencies of the system
Items

First order

Second order

Third order

Fourth order

Frequency (Hz)

51.35

261.15

376.90

418.39

The dynamic loads from the rocket are mainly considered here for the dynamic response analysis of the system. Assume that $\mathbf{w}={\left[\begin{array}{lll}0.05\mathrm{s}\mathrm{i}\mathrm{n}\left(5t\right)& zeros\left(\mathrm{5,1}\right)& 1\end{array}\right]}^{T}$, namely, only consider the lateral disturbance, the response of the system then can be obtained as shown in Fig. 8. The result shows that the lateral disturbance not only excite the lateral response (Fig. 8(a)), but also excite the rotational response (Fig. 8(b)) and longitudinal response (Fig. 8(c)) of the system.
Fig. 9 shows the response of the system due to the longitudinal excitation. It is noted that only the vertical response be excited (red line in Fig. 9), as a result of the assumption that the system is a symmetrical one.
Fig. 8. Response due to lateral excitation: a) lateral response, b) rotational response, c) longitudinal response, d) responses in other directions)
Fig. 9. Response due to vertical excitation
According to Faraday’s Law of Electromagnetic Induction, the electromotive force ($E=BLv$) and current ($I=BLv/R$) are induced when the closed coils moves in magnetic field. After substituting them into Eq. (4) the Ampere’s force then can be written as:
where $R$ is the resistance of the coils, $v$ is relative velocity of coils.
Thus, the electromagnetic actuator can equivalent to a viscous damper, when shorts out the coil of the actuator. And the damping coefficient:
is a constant only related to the physic parameters of the actuator. Then, the damping force due to $j$th damper can be written in the global coordinate frame as:
and the total damping force then be described as:
Substituting Eq. (34) and Eq. (35) into Eq. (28), the dynamic equation of the damped system can be written as:
Fig. 10 shows the response of the system subject to a fixed frequency excitation $0.05\mathrm{s}\mathrm{i}\mathrm{n}\left(150t\right)$ (N) in the vertical direction, the dotted line represents the response of the original system and the solid line represent the response of the damped system, respectively. It is noted that the response of the damped system has a lower magnitude compared to that of the original system. The free responses with an initial displacement 0.05 m and initial velocity 0.05 m/s, are shown in Fig. 11. It can be easily observed from Fig. 11 that the damped system attenuated rapidly.
Fig. 10. Response under vertical excitation
Fig. 11. Free response in vertical direction
Fig. 12 to Fig. 14 shows that the magnitudes of the response of the system with a fixed frequency excitation $0.05\mathrm{s}\mathrm{i}\mathrm{n}\left(150t\right)$ (N) in the lateral direction, and a lower magnitude obtained after considering the electromagnetic damping.
Fig. 12. Lateral response under lateral excitation
Fig. 13. Rotational response under lateral excitation
_{}
Fig. 14. Vertical response under lateral excitation
The free responses with an initial displacement 0.05 m and initial velocity 0.05 m/s in the lateral direction are shown in Fig. 15 to Fig. 17, respectively. The damped system illustrated the same attenuation trend compared to the original system.
_{}
Fig. 15. Lateral free response (with an initial displacement and velocity in lateral direction)
Fig. 16. Rotational free response (with an initial displacement and velocity in lateral direction)
Fig. 17. Vertical free response (with an initial displacement and velocity in lateral direction)
It is observed that the WSVI has a good isolation performance when considering the VCMs as viscous dampers. The response amplitude of the system can be suppressed effectively, and it is also noted that the system has a higher decay speed in vertical direction than lateral direction, this is because of the actuators are all mounted vertically, which lead to a larger damping force in the vertical direction.
A scaling model is built, and an experimental setup is constructed (see Fig. 18), sinusoidal sweepingfrequency vibration tests were carried out on the experimental device in both lateral and longitudinal directions for further verification the isolation performance of the proposed WSVI platform. In the two experiment cases, the shaking table provides a sweeping acceleration signals from 5 Hz to 100 Hz with an amplitude of 1 m/s^{2} and 5 m/s^{2} respectively as the disturbances. The undamped and damped systems are tested in each case, where the undamped and damped system corresponding to the open and short circuit of the VCM, respectively. The amplitudefrequency responses captured from the two vibration experiments are presented in Fig. 19. It should be noted that, for better demonstrate the isolation perform of the WSVI platform the stiffness of the leaf springs employed in this experiment is lower than that of the simulation.
Fig. 18. Photo of the WSVI platform and experiment setup
It can be observed form the experiment results that in the lateral experiment case the amplitude at the resonance peak is significantly decreased from 29.15 m/s^{2} to 11.23 m/s^{2}, and in the longitudinal experiment case the response at the resonance frequency of the system is attenuated from 88.88 m/s^{2} to 27.46 m/s^{2} due to the contribution of the electromagnetic damping.
Fig. 19. Frequency response of the payload
a) Lateral
b) Longitudinal
In this work, a novel wholespacecraft vibration isolation device has been developed based on the VCM. An analytical model for the generated force of the VCM is derived by applying the electromagnetic equations. Then, the structural features of the proposed WSVI device were illustrated and the kinematics equation for the system were formulated. The dynamic model and the state space equation of the system were developed for the study of the dynamic characteristics. In addition, the physics parameters of a scaleddown WSVI device were given and then, the natural frequency, the dynamic response under exogenous disturbances were obtained. Finally, based on the electromagnetism theorem, the electromagnetic actuator was turned into viscous damper, and the isolation performance of the WSVI system was verified though the simulation and experiment.
The newly proposed WSVI device features small space footprint, light weight. It can satisfy the design requirement of WSVI without changing of the PAF structure with its novel design. From the passive vibration isolation results, it is clear that the new WSVI device performed very well to reduce structure vibration transmitted to the spacecraft. What should be noted is that, this paper mainly discussed the passive isolation performance of the proposed device, the future research will focus on the investigation of active vibration isolation performance.
The authors are grateful to the Foundation of Civil Astronautics and the Preresearch Foundation of General Armament Department of China for financial support in this study.