Design and Analysis of a Compression and Separation Device for Multi-Satellite Deployment

The launch method of one arrow with multiple satellites can greatly shorten the time for constellation networking and improve the deployment efficiency. A new compression and separation device with a four-bar perimeter arrangement is proposed for multi-satellite compaction and inorbit release. A compression device with gap elimination is designed to implement the reliable compaction of stacked flat satellites. An electromagnetic separation device is proposed to achieve the fast, low-interference release of multi-satellites. The dynamic model with flexible guide bars is established. The separation characteristics of multiple satellites are analyzed by the kinematic simulation. The prototype is developed, and the related experiment is implemented. The results show that the four-guide-bar-edge arrangement scheme with a gap elimination device achieves reliable locking and fast separation under a vibration environment. The dynamic separation characteristics of satellites are investigated by the air floatation experiments. The results show that a stable separation speed and low disturbance angular velocity are achieved under 10% spring error.


Introduction
Multiple satellites working in concert can implement space tasks that cannot be accomplished by a single robust satellite, such as remote sensing, constellations, and formation [1][2][3]. The satellites are generally launched and released in a single rocket, which can greatly shorten the time of constellation networking and improve the launch efficiency. The compression and separation devices for multi-satellite launch deployment scenarios are mainly divided into three categories: box-storage type [4][5][6], wall-mounted type [7], and stacked compression type [8,9]. The box-storage type is applied to the storage and release of multiple nano-satellites and pico-satellites of the same or similar size. The wall-mounted type is suitable to the launch scenario of multiple rectangular satellites or satellites with trapezoidal cross-section. The stacked compression type is applied to compress and release the multiple flat-shaped satellites, such as Starlink and Smart Dusts [10]. In the paper, the stacked compression method is adopted to lock and release multiple flat satellites.
In a multi-satellite launch, the design of the locking connection and separation system for multiple satellites is the key technology to be solved first. The layout and installation of multiple satellites in the rocket, the anti-collision problem of separation, and the connection method to minimize the impact between satellites and ensure effective separation into orbit are issues that need to be considered when designing the locking connection and separation system. In the compression connection and separation system, traditional locking release devices and initiating explosive devices [11] are used on spacecraft, such as pin pullers, explosive bolts, initiating explosive nuts, and cutters. The unlocking impact of initiating explosive devices is large, reaching 3000~20,000 g. It also has the disadvantages of hard storage and transportation, safety risks, and can only be used once so they cannot be tested rection. The satellite group is pre-tightened by the centralized compression device b tween the satellite group and the carrier platform. The pre-tightening force is maintain by the unlocker to compress the satellite group so that the system has sufficient stiffne After the centralized device is unlocked, the controllable ejection device is responsible f the connection between satellites and between satellites and carrier platforms. When t satellites are separated, the controllable connection separation mechanism acts accordi to the predetermined deployment time to remove the connection between the satell and the carrier platform to realize controllable separation of the satellites one by one.

Compression Mechanism
The compression mechanism is the key component to achieving high load-bearin and its structural form of compression and release directly determines high load-bearin and the performance of the subsequent mechanism. The compression mechanism mainly connected to the base and the satellite bearing column at the top through t guide bar, which is connected by thread, and the unlocker unlocks and releases. The d tailed structure composition is shown in Figure 3.
The guide bar passes through the interior of the satellite load-bearing column connect the force and displacement of the base section to the top unlocking section. T upper end of the guide bar is connected to the upper joint flange. The load-bearing scre passes through the upper joint flange. The lower end is connected to the preload nut a the locknut, and the upper end is connected to the non-self-locking nut. T non-self-locking nut is arranged in the unlocker, and the lower surface of the unlock base plate is in contact with the upper surface of the uppermost satellite load-beari rection. The satellite group is pre-tightened by the centralized compression device between the satellite group and the carrier platform. The pre-tightening force is maintained by the unlocker to compress the satellite group so that the system has sufficient stiffness. After the centralized device is unlocked, the controllable ejection device is responsible for the connection between satellites and between satellites and carrier platforms. When the satellites are separated, the controllable connection separation mechanism acts according to the predetermined deployment time to remove the connection between the satellite and the carrier platform to realize controllable separation of the satellites one by one.

Compression Mechanism
The compression mechanism is the key component to achieving high load-bearing, and its structural form of compression and release directly determines high load-bearing, and the performance of the subsequent mechanism. The compression mechanism is mainly connected to the base and the satellite bearing column at the top through the guide bar, which is connected by thread, and the unlocker unlocks and releases. The detailed structure composition is shown in Figure 3.
The guide bar passes through the interior of the satellite load-bearing column to connect the force and displacement of the base section to the top unlocking section. The upper end of the guide bar is connected to the upper joint flange. The load-bearing screw passes through the upper joint flange. The lower end is connected to the preload nut and the locknut, and the upper end is connected to the non-self-locking nut. The non-self-locking nut is arranged in the unlocker, and the lower surface of the unlocker

Compression Mechanism
The compression mechanism is the key component to achieving high load-bearing, and its structural form of compression and release directly determines high load-bearing, and the performance of the subsequent mechanism. The compression mechanism is mainly connected to the base and the satellite bearing column at the top through the guide bar, which is connected by thread, and the unlocker unlocks and releases. The detailed structure composition is shown in Figure 3.
column. The lower end of the guide bar is connected to the lower joint flange. The loading nut and the conical head nut form a threaded connection with the lower joint flange, and the loading nut contacts the spherical washer, which contacts the spherical surface of the fixed base. Figure 3. Structure of compression mechanism: 1-loader post; 2-guide bar; 3-locknut; 4-preload nut; 5-upper joint flange; 6-load-bearing screw; 7-non-self-locking nut; 8-unlocker; 9-electrical connector upper interface; 10-electrical connector lower interface; 11-fixed base; 12-conical head nut; 13-retention spring; 14-lower joint flange; 15-spherical washer; 16-loading nut; 17-return spring; 18-holding base.
In the pressed state, the loading nut squeezes the spherical washer by applying a pressing force. Therefore, a vertical downward reaction force is received from the fixed base, and the reaction force is transmitted to the guide bar. The guide bar applies a vertical downward preload force to the preload nut through the upper joint flange, so as to transmit the preload force to the unlocker. Then the force is transmitted to the load-bearing columns of the satellite through the lower base plate of the unlocker to realize the compression of the satellite group.
When unlocking and releasing, the unlocker is energized to release the circumferential restriction on the non-self-locking nut. The guide bar drives the load-bearing screw to move downward, freeing it from the thread constraint of the non-self-locking nut and releasing the compression force on the satellite group. The lower end of the guide bar allows the conical head nut to be embedded in the tapered groove of the retention base by the action of the retention spring, keeping the compression assembly in a fixed position in the microgravity environment.

Separation Mechanism
In order to realize the controllable separation function of the satellites, the inter-satellite controllable connection and separation mechanism is designed by using a combination of the electromagnetic connection and the spring ejection. Since only the load-bearing column in the satellite structure can bear the load, the electromagnetic connection assembly and the ejection separation assembly are both arranged at the contact connection of the load-bearing column. To ensure the uniformity of the force, four electromagnetic connection assemblies and ejection separation assemblies are added between each group of satellites, as shown in Figure 4. . Structure of compression mechanism: 1-loader post; 2-guide bar; 3-locknut; 4-preload nut; 5-upper joint flange; 6-load-bearing screw; 7-non-self-locking nut; 8-unlocker; 9-electrical connector upper interface; 10-electrical connector lower interface; 11-fixed base; 12-conical head nut; 13-retention spring; 14-lower joint flange; 15-spherical washer; 16-loading nut; 17-return spring; 18-holding base.
The guide bar passes through the interior of the satellite load-bearing column to connect the force and displacement of the base section to the top unlocking section. The upper end of the guide bar is connected to the upper joint flange. The load-bearing screw passes through the upper joint flange. The lower end is connected to the preload nut and the locknut, and the upper end is connected to the non-self-locking nut. The non-self-locking nut is arranged in the unlocker, and the lower surface of the unlocker base plate is in contact with the upper surface of the uppermost satellite load-bearing column. The lower end of the guide bar is connected to the lower joint flange. The loading nut and the conical head nut form a threaded connection with the lower joint flange, and the loading nut contacts the spherical washer, which contacts the spherical surface of the fixed base.
In the pressed state, the loading nut squeezes the spherical washer by applying a pressing force. Therefore, a vertical downward reaction force is received from the fixed base, and the reaction force is transmitted to the guide bar. The guide bar applies a vertical downward preload force to the preload nut through the upper joint flange, so as to transmit the preload force to the unlocker. Then the force is transmitted to the load-bearing columns of the satellite through the lower base plate of the unlocker to realize the compression of the satellite group.
When unlocking and releasing, the unlocker is energized to release the circumferential restriction on the non-self-locking nut. The guide bar drives the load-bearing screw to move downward, freeing it from the thread constraint of the non-self-locking nut and releasing the compression force on the satellite group. The lower end of the guide bar allows the conical head nut to be embedded in the tapered groove of the retention base by the action of the retention spring, keeping the compression assembly in a fixed position in the microgravity environment.

Separation Mechanism
In order to realize the controllable separation function of the satellites, the inter-satellite controllable connection and separation mechanism is designed by using a combination of the electromagnetic connection and the spring ejection. Since only the load-bearing column in the satellite structure can bear the load, the electromagnetic connection assembly and the ejection separation assembly are both arranged at the contact connection of the load-bearing column. To ensure the uniformity of the force, four electromagnetic connection assemblies and ejection separation assemblies are added between each group of satellites, as shown in Figure 4.  The electromagnet in the electromagnetic connection assembly is in clearance fi with the inner hole of the lower base of the electromagnetic connection. The electro magnetic connection base has two symmetrical slots, and the magnet connection flange is designed with two symmetrical bosses that fit with its clearance. The electromagne preload spring provides thrust to ensure that the electromagnet can maintain close con tact with the iron block despite errors due to processing and assembly. The sliding sleeve guide rod in the ejection separation assembly is in clearance fit with the bore of the lowe base of the electromagnetic connection.

Natural Coordinate Method Based Rigid Body Satellite Dynamic Model
For rigid body satellites, the natural coordinate method is used to describe the pla nar unit. The natural coordinate method is used to describe the planar multi-body prob lem by selecting a series of representative "nodes", whose Cartesian coordinates consti tute the natural coordinates of the multi-body system together. Figure 5 shows the representation of a planar rigid satellite based on the NCF method, whose motion position can be determined by two nodes,  The electromagnet in the electromagnetic connection assembly is in clearance fit with the inner hole of the lower base of the electromagnetic connection. The electromagnetic connection base has two symmetrical slots, and the magnet connection flange is designed with two symmetrical bosses that fit with its clearance. The electromagnet preload spring provides thrust to ensure that the electromagnet can maintain close contact with the iron block despite errors due to processing and assembly. The sliding sleeve guide rod in the ejection separation assembly is in clearance fit with the bore of the lower base of the electromagnetic connection.

Natural Coordinate Method Based Rigid Body Satellite Dynamic Model
For rigid body satellites, the natural coordinate method is used to describe the planar unit. The natural coordinate method is used to describe the planar multi-body problem by selecting a series of representative "nodes", whose Cartesian coordinates constitute the natural coordinates of the multi-body system together. Figure 5 shows the representation of a planar rigid satellite based on the NCF method, whose motion position can be determined by two nodes, P 1 and P c and A, containing four generalized coordinates, as shown in Table 1.
The generalized coordinate vector of a planar rigid satellite can be expr The natural coordinates corresponding to the planar rigid body satell sented by the matrix as: The total number of natural coordinates of the system is 2n, expressed coordinate array as: The generalized coordinate vector of a planar rigid satellite can be expressed as: The natural coordinates corresponding to the planar rigid body satellite are represented by the matrix as: The total number of natural coordinates of the system is 2n, expressed in a natural coordinate array as: Set n i as the normal base vector of the motion plane, select vector u i = n i × r (i) 12 as the reference vector, and u i is orthogonal to r (i) 12 and has the same modulus l (i) in the motion plane. With the base point P Decompose the vector diameter 1 of any point P (i) within the rigid body relative to the base point P  Representing each vector as a coordinate array in O 0 , E (0) , we can obtain: The mass of the rigid body is m i , and the vector diameter R (i) represented by Equation (6) is integrated within the range of the rigid body. We obtain the vector diameter R (i) c of the rigid body center of mass O ci relative to O 0 as: The matrix D (i) is defined as: Let the vector diameters of the point P (i) and the point P After that 1 , E (i) as follows: The rate of change of vector r (i) 12 with respect to time t due to the rotation of a rigid satellite is: Let both sides be multiplied by the r (i) 12 vector, and express the vector r (i) 12 as a coordinate array to obtain the angular velocity expressed in natural coordinates as: G (i) is a 1 × 4 row matrix, expressed as: The rotational moment of inertia of the rigid satellite with respect to the center of mass O ci is J i . The coordinate array L (i) of the momentum of the rigid satellite with respect to J i G (i) is the inertia matrix of the rigid satellite with respect to the center of mass O ci .
Taking the derivative of Equation (5)  Q i = 0, we obtain: .  (7) and (15) into the momentum theorem and the momentum moment theorem for the center of mass, the dynamics equation of the rigid satellite plane motion is obtained as:

Dynamic Model of Flexible Guide Bar Based on Floating Coordinate System Method
Due to the relative displacement between the mass points during the motion of the flexible body, no reference system can be completely fixed with the flexible body. Taking the floating coordinate system as the reference system, the actual motion of the flexible body can be understood as the synthesis of the large-scale rigid body motion of the floating coordinate system relative to the inertial coordinate system and the deformation motion relative to the floating coordinate system, and the flexible guide bar and its floating coordinate system are shown in Figure 6.
force acting on the rigid satellite and the principal moment relative to the ( ) (7) and (15) into the momentum the momentum moment theorem for the center of mass, the dynamics equatio satellite plane motion is obtained as:

Dynamic Model of Flexible Guide Bar Based on Floating Coordinate System M
Due to the relative displacement between the mass points during the flexible body, no reference system can be completely fixed with the flexible the floating coordinate system as the reference system, the actual motion body can be understood as the synthesis of the large-scale rigid body m floating coordinate system relative to the inertial coordinate system and the motion relative to the floating coordinate system, and the flexible guid floating coordinate system are shown in Figure 6.
The mass matrix of the flexible body is: The floating coordinate system of the flexible guide bar is established at the midpoint of the fixed constraint at the bottom end, and the bottom end of the guide bar is a fixed constraint, so it does not produce a wide range of rigid body motion. The motion of the guide bar is its own flexible deformation.
The relative displacement u at the point P can be expanded by using the modal Ψ j (x, y)(j = 1, 2, . . . , n) as: where q j (t)(j = 1, 2, . . . , n) are the modal coordinates. The mass matrix of the flexible body is: Each submatrix is defined as: The stiffness matrix of the flexible body is: For a plane beam the stiffness matrix K can be simplified as: The flexible guide bar generalized coordinate array Q 2 is represented by the rigid motion coordinates R O1 , the azimuthal Euler angular coordinates θ and the modal coordinates Q f , denoted as: The Lagrange multiplier method is applied to establish the kinetic equation as: Coupled with the constraint equation, the multi-body dynamic equation of the flexible guide bar is obtained as: In order to analyze the dynamics of the long-range guide separation mechanism, a contact collision model of the rigid body satellite and the flexible guide bar is required.

Contact Force Model
The geometric description and vector model of the contact between the satellite and the guide bar are as follows.

Contact Form of Rigid Body Satellite and Guide Bars
To represent the relationship between the rigid satellite and the flexible guide bar clearly, the clearance distance between the rigid satellite guide hole and the guide bar is enlarged. There are four possible forms of contact between the rigid satellite and the guide bars, as shown in Figure 7.

1.
No contact between the two elements in the radial direction; 2.
Point contact occurring between the side of the guide bar and one endpoint of the satellite guide hole; 3.
Point contact occurring between the sides of two guide bars and two endpoints of the satellite guide hole; 4.
The sides of the two guide bars overlap with the sides of the satellite guide holes leading to line contact; 2. Point contact occurring between the side of the guide bar and one endpoint of t satellite guide hole; 3. Point contact occurring between the sides of two guide bars and two endpoints the satellite guide hole; 4. The sides of the two guide bars overlap with the sides of the satellite guide ho leading to line contact;

Judgment Basis of Contact State
The contact forms between the rigid satellite and the guide bars are complex a varied. How to discern the contact and separation state between the two elements is a k issue in the study of the guided separation process.
The prerequisite to determining whether and when the contact between the satell and the guide bar occurs is to find the first point where the contact occurs, i.e., to det mine the potential contact point, so as to establish the expression of the embedding dep between the two elements caused by the contact deformation. A simplified model solving the embedding depth is shown in Figure 8. The embedding depth is solved the left guide hole of the satellite as an example, and the situation is similar on the oth side.

Judgment Basis of Contact State
The contact forms between the rigid satellite and the guide bars are complex and varied. How to discern the contact and separation state between the two elements is a key issue in the study of the guided separation process.
The prerequisite to determining whether and when the contact between the satellite and the guide bar occurs is to find the first point where the contact occurs, i.e., to determine the potential contact point, so as to establish the expression of the embedding depth between the two elements caused by the contact deformation. A simplified model for solving the embedding depth is shown in Figure 8. The embedding depth is solved for the left guide hole of the satellite as an example, and the situation is similar on the other side. Vector r A and vector r H can be expressed as follows: Because of: We can obtain: The embedding depth can be expressed as: The embedding depth δ is a function of coordinates of point O c3 and point M, i.e., δ = δ (O c3 , M).
The modeling of the contact collision force is the key to describing the interaction force of the guide and separation mechanism. According to the different relative motion states between the rigid body satellite and the flexible body guide bar, the contact collision force also has three forms. When two elements are in contact with each other, a contact force is generated between the elements, which is F c . The contact force is F c decomposed into the collision force F n in the normal direction of the contact point and the friction force F t in the tangential direction, as shown in Figure 9, which transforms the contact force problem into the collision and friction problem during the motion of the mechanical elements.
The embedding depth  is a function of coordinates of point 3 The modeling of the contact collision force is the key to describing the interaction force of the guide and separation mechanism. According to the different relative motion states between the rigid body satellite and the flexible body guide bar, the contact collision force also has three forms. When two elements are in contact with each other, a contact force is generated between the elements, which is c F . The contact force is c F decomposed into the collision force n F in the normal direction of the contact point and the friction force t F in the tangential direction, as shown in Figure 9, which transforms the contact force problem into the collision and friction problem during the motion of the mechanical elements. Obviously, the points of action of the contact forces acting at contact points A and H are at different locations. To establish a unified expression, the contact force is converted into a contact force acting on the center of mass of the satellite and the guide bar. During Obviously, the points of action of the contact forces acting at contact points A and H are at different locations. To establish a unified expression, the contact force is converted into a contact force acting on the center of mass of the satellite and the guide bar. During the conversion of the contact force, an additional moment of action is generated for the satellite and the guide bar, which acts equivalently at the center of mass.
Since the normal collision force model and the tangential friction force model contain relative velocity terms, the relative velocities of the two elements moving in contact with each other are also determined in order to calculate the normal collision force and the tangential friction force between the two elements.
Derivation of Equation (25) yields the velocity at the point of contact as: .
The relative velocity of the contact point can be expressed as: Projecting the relative velocity to the normal and tangential directions of the contact point, the relative normal velocity and relative tangential velocity can be obtained as: By substituting the relative velocities of the contact points into the normal collision force model and the tangential friction model, the contact forces can be obtained.
Based on the Hertz contact model, a contact force model considering hysteresis damping called the L-N nonlinear spring damping model was proposed by Lankarani and Nikravesh.
The contact force is expressed using the L-N contact force model as: where K is the nonlinear spring stiffness coefficient, δ is embedding depth, D is generalized damping coefficient, . δ is collision relative velocity, and C e is recovery coefficient, is relative velocity at the initial moment of contact.
Substituting the expression for the generalized damping factor, the contact force can be obtained as: This contact force model is suitable for contact between satellites and guide bars, which has low-speed collision and large recovery coefficient. Figure 10a shows the classical model for expressing friction, which fits Coulomb friction model. Friction is proportional to the normal contact force, and when the normal contact force is constant, the friction is constant and the direction is opposite to the speed of motion. Friction can be expressed as: where µ is the friction coefficient between the surfaces of two objects, F t is tangential friction, and v is the relative velocity between the two objects. The improved Coulomb friction model introduces a dynamic modified coefficien cd, which is related to the relative velocity between objects, and its expression is as fo lows: where v1 and v0 are relative velocity thresholds.

Separation Dynamic Model with Contact Collisions
The relative motion of the satellite and the guide bar is expressed in four states: fre state, collision state, detached state, and contact state. The free state is defined as the stat in which the satellite and the guide bar are separated from each other with no collision o contact. In this state, there is no contact force between the inner hole of the satellit load-bearing column and the guide bar. The contact state is defined as the state in whic the surface deformation occurs between the inner hole of the satellite load-bearing co umn and the outer surface of the guide bar by mutual extrusion. The contact state gen erates contact force, which can be decomposed into normal collision force due to norma relative motion and tangential friction force due to tangential motion. The detached stat and the collision state are two transition stages. The detached state is the transition stag from the contact state to the free state. The collision state is the transition stage from th free state to the contact state.
When the free state ends, the motion state is changed into the collision state. Whe the collision state ends, the motion state is changed into the contact state. The next mo ment after the end of the contact state, the motion state is transferred to the detached stat The next moment after the end of the detached state, the motion state is immediatel transferred to the free state, and so on. Therefore, it is necessary to judge the relativ motion state of the moving subelements, because each change of motion state means th beginning of contact or separation between the moving subelements. The Coulomb friction model with modified coefficient introduced is shown in Figure 10b, and it can be expressed as: where c f is coefficient of friction between the two surfaces and v t is the relative velocity between the two objects.
The improved Coulomb friction model introduces a dynamic modified coefficient, c d , which is related to the relative velocity between objects, and its expression is as follows: where v 1 and v 0 are relative velocity thresholds.

Separation Dynamic Model with Contact Collisions
The relative motion of the satellite and the guide bar is expressed in four states: free state, collision state, detached state, and contact state. The free state is defined as the state in which the satellite and the guide bar are separated from each other with no collision or contact. In this state, there is no contact force between the inner hole of the satellite load-bearing column and the guide bar. The contact state is defined as the state in which the surface deformation occurs between the inner hole of the satellite load-bearing column and the outer surface of the guide bar by mutual extrusion. The contact state generates contact force, which can be decomposed into normal collision force due to normal relative motion and tangential friction force due to tangential motion. The detached state and the collision state are two transition stages. The detached state is the transition stage from the contact state to the free state. The collision state is the transition stage from the free state to the contact state.
When the free state ends, the motion state is changed into the collision state. When the collision state ends, the motion state is changed into the contact state. The next moment after the end of the contact state, the motion state is transferred to the detached state. The next moment after the end of the detached state, the motion state is immediately transferred to the free state, and so on. Therefore, it is necessary to judge the relative motion state of the moving subelements, because each change of motion state means the beginning of contact or separation between the moving subelements.
Firstly, the relative motion state at the current moment is determined by the embedding depth δ derived in the previous section.
When δ < 0, the satellite and the guide bar have not been in contact and the relative motion state is free state.
When δ = 0, the satellite and the guide bar begin to collide or begin to separate. When δ > 0, the satellite continuously contacts with the guide bar. The relationship between δ of two adjacent time steps and the relative velocity between the two in the collision normal direction at the current moment are used to further determine the switching of the motion state.
The initial relative motion state of the satellite and the guide bar is the free state. In the process of separation motion, the flow chart of the judgment method to discern the relative motion state of the satellite and the guide bar is shown in Figure 11.
The initial relative motion state of the satellite and the guide bar is the free state. In the process of separation motion, the flow chart of the judgment method to discern the relative motion state of the satellite and the guide bar is shown in Figure 11. In order to introduce the contact forces into the flexible multi-body dynamics equation and establish a unified form of the dynamics equation, introduce a step function u( ).
The contact force Fc can be expressed as: Then, the dynamic model of satellite guided separation system with contact interaction can be obtained as: From the dynamic model, it can be seen that the main factors affecting the motion parameters in the guided satellite separation process are the length of the load-bearing column, the inner hole diameter of the load-bearing column, the diameter of the guide bar, and the initial velocity and initial angular velocity of the satellite. The subsequent In order to introduce the contact forces into the flexible multi-body dynamics equation and establish a unified form of the dynamics equation, introduce a step function u(δ).
The contact force F c can be expressed as: Then, the dynamic model of satellite guided separation system with contact interaction can be obtained as: M ..
From the dynamic model, it can be seen that the main factors affecting the motion parameters in the guided satellite separation process are the length of the load-bearing column, the inner hole diameter of the load-bearing column, the diameter of the guide bar, and the initial velocity and initial angular velocity of the satellite. The subsequent content will present dynamics simulation analysis and experimental investigation of the key factors.

Separation Characteristics Simulation
Since the sixth satellite is located at the bottom and released last, it has the longest travel inside the guide bars. So, we select the sixth satellite to analyze the separation characteristics. As shown in Figure 12, it spends about 9.5 s for the satellite to be completely separated from the guide bar with a length of 3.85 m. The separation velocity and acceleration are shown in Figure 13. The separation speed fluctuates due to the collision between the satellite and the guide bars. The maximum acceleration is about 1.5 m/s 2 . Due to the friction and collision with the guide rod, several smaller negative accelerations occur on the z-axis. The angular acceleration around the x-axis of the satellite is shown in Figure 14.
Because of the stiffness error of the spring, the satellite has a negative angular acceleration around the x-axis.
Aerospace 2022, 9, x FOR PEER REVIEW between the satellite and the guide bars. The maximum acceleration is abo Due to the friction and collision with the guide rod, several smaller negative ac occur on the z-axis. The angular acceleration around the x-axis of the satellite Figure 14. Because of the stiffness error of the spring, the satellite has a negat acceleration around the x-axis.     Aerospace 2022, 9, x FOR PEER REVIEW between the satellite and the guide bars. The maximum acceleration is about 1 Due to the friction and collision with the guide rod, several smaller negative accel occur on the z-axis. The angular acceleration around the x-axis of the satellite is sh Figure 14. Because of the stiffness error of the spring, the satellite has a negative acceleration around the x-axis.     Aerospace 2022, 9, x FOR PEER REVIEW between the satellite and the guide bars. The maximum acceleration is about 1 Due to the friction and collision with the guide rod, several smaller negative accel occur on the z-axis. The angular acceleration around the x-axis of the satellite is s Figure 14. Because of the stiffness error of the spring, the satellite has a negative acceleration around the x-axis.      Figure 15 shows comparison of the separation characteristics of the six satellites. The maximum velocity of the satellites is all 0.347 m/s due to the same spring. The satellites collide with the guide bars several times, causing the velocity to decrease. Because of the different number of collisions, there is a difference in the velocity after separation. The spring difference leads to a large angular velocity of the satellite, which is well suppressed after many contact collisions with the guide bars. However, due to the short guided stroke of the first satellite, the angular velocity reached 2.47 • /s at separation, while the other satellites kept it within 1 • /s.

Experiment Investigation
To investigate the function and characteristics of the compression and release device, experiments were conducted on the compression and unlocking characteristics, and the dynamic characteristics of the satellite release, respectively.

Compression Test
In order to check whether the structures in the system can withstand the sinusoidal and random vibration environment during the rocket launch, and to investigate the defects in structural design and material selection as well as processing and assembly, the sinusoidal and random vibration tests were conducted on the system. The two test satellites were mounted on the adapter table and connected by the centralized compression and release mechanism. The centralized compression and release mechanism applied a compression force of 45 kN to keep the system in a compressed state, as shown in Figure  16. The vibration amplitude is 1.17 mm when the sine vibration frequency is 5~8 Hz. The vibration amplitude is 0.3 g when the frequency is 8~100 Hz. The random vibration input conditions are shown in Table 2.
The results of the vibration experiments are shown in Figure 17. Under the action of

Experiment Investigation
To investigate the function and characteristics of the compression and release device, experiments were conducted on the compression and unlocking characteristics, and the dynamic characteristics of the satellite release, respectively.

Compression Test
In order to check whether the structures in the system can withstand the sinusoidal and random vibration environment during the rocket launch, and to investigate the defects in structural design and material selection as well as processing and assembly, the sinusoidal and random vibration tests were conducted on the system. The two test satellites were mounted on the adapter table and connected by the centralized compression and release mechanism. The centralized compression and release mechanism applied a compression force of 45 kN to keep the system in a compressed state, as shown in Figure 16. The vibration amplitude is 1.17 mm when the sine vibration frequency is 5~8 Hz. The vibration amplitude is 0.3 g when the frequency is 8~100 Hz. The random vibration input conditions are shown in Table 2.  The successful unlocking of the centralized compression and release device is important for the smooth separation of the satellite in orbit. Especially after being subjected to complex vibration environment, it can release the compression force and realize the unlocking function. The unlocking separation experiment is shown in Figure 18.
The counterweight suspension is used to simulate the space microgravity environment. Each satellite is connected to the counterweight at four corners through wire ropes and fixed pulleys. The mass of the counterweight is equal to that of the test satellite and is proportionally distributed according to the position of the center of mass of the test satellite. When the centralized compression and release mechanism is unlocked, the testing satellite is separated along the vertical direction under the thrust of the separation spring. The unlocking process of the centralized compression and release mechanism starts at the moment of SMA wire is energized and ends at the moment when the compression  The results of the vibration experiments are shown in Figure 17. Under the action of sinusoidal vibration with the frequency range of 5 to 100 Hz, the x-direction and y-direction resonate at 29.1 Hz and 31.2 Hz, respectively, producing a maximum response of 4.53 g and 4.29 g. The z-direction resonates at 77.6 Hz with a maximum response of only 0.68 g, reflecting the good compression performance. According to the static strength analysis, it is known that the compression and release device can withstand an acceleration response of 7.5 g in the z-direction, which verifies that the compression and release device can be reliably locked under the given vibration conditions. In the frequency range from 10 to 2000 Hz, the resonance occurs at 585 Hz and 580 Hz in the x-direction and y-direction respectively by random vibration, and the maximum response is 0.95 g 2 /Hz and 0.81 g 2 /Hz respectively. The resonance occurs at 530 Hz in the z-direction, and the maximum response is only 0.09 g 2 /Hz. Therefore, the experiments reflect the good compression performance.

Unlocking Experiment
The successful unlocking of the centralized compression and release device is important for the smooth separation of the satellite in orbit. Especially after being subjected to complex vibration environment, it can release the compression force and realize the unlocking function. The unlocking separation experiment is shown in Figure 18.
force of the device on the satellite group is released to zero. The time experienced in process is defined as the unlocking time. In the practical test, it is difficult to identif moment when the compression force is released to zero. So, the moment when load-bearing screw is disengaged from the unlocker is used as the endpoint. The ave unlocking and releasing time was 312 ms after three tests with input power of 5 A plied to the unlocker from a DC power. This experiment indicates that the centra compression and release mechanism can still unlock and release the pressing force predetermined time to achieve the pressing and releasing function after experienc complex vibration environment.

Separation Characteristics Experiment
The scaled-down prototype and experimental setup are shown in Figure 19 satellite pedestal is fixed on the side of the marble platform. The four electromag actuators and the ejection separation assembly are mounted at the four corners o base attachment plate. The satellite is located above the air-floating platform to elim the effect of friction during the separation process. The gas cylinder is mounted insid satellite The satellite is kept in the initial position by the suction of the electromag actuator to overcome the thrust of the separation spring. A high-speed camera rec the whole process of the motion of the target point located on the satellite. The an velocity sensor is mounted on the satellite to record the angular velocity. Figure 19. Separation velocity test. The counterweight suspension is used to simulate the space microgravity environment. Each satellite is connected to the counterweight at four corners through wire ropes and fixed pulleys. The mass of the counterweight is equal to that of the test satellite and is proportionally distributed according to the position of the center of mass of the test satellite. When the centralized compression and release mechanism is unlocked, the testing satellite is separated along the vertical direction under the thrust of the separation spring. The unlocking process of the centralized compression and release mechanism starts at the moment of SMA wire is energized and ends at the moment when the compression force of the device on the satellite group is released to zero. The time experienced in this process is defined as the unlocking time. In the practical test, it is difficult to identify the moment when the compression force is released to zero. So, the moment when the load-bearing screw is disengaged from the unlocker is used as the endpoint. The average unlocking and releasing time was 312 ms after three tests with input power of 5 A supplied to the unlocker from a DC power. This experiment indicates that the centralized compression and release mechanism can still unlock and release the pressing force at a predetermined time to achieve the pressing and releasing function after experiencing a complex vibration environment.

Separation Characteristics Experiment
The scaled-down prototype and experimental setup are shown in Figure 19. The satellite pedestal is fixed on the side of the marble platform. The four electromagnetic actuators and the ejection separation assembly are mounted at the four corners of the base attachment plate. The satellite is located above the air-floating platform to eliminate the effect of friction during the separation process. The gas cylinder is mounted inside the satellite The satellite is kept in the initial position by the suction of the electromagnetic actuator to overcome the thrust of the separation spring. A high-speed camera records the whole process of the motion of the target point located on the satellite. The angular velocity sensor is mounted on the satellite to record the angular velocity. base attachment plate. The satellite is located above the air-floating platform to eliminat the effect of friction during the separation process. The gas cylinder is mounted inside th satellite The satellite is kept in the initial position by the suction of the electromagneti actuator to overcome the thrust of the separation spring. A high-speed camera record the whole process of the motion of the target point located on the satellite. The angula velocity sensor is mounted on the satellite to record the angular velocity.   Figure 20 illustrates the test results of the separation speed. Figure 20a shows that the maximum separation velocity of the satellite is 0.345 m/s. The satellite collides with the guide bars, resulting in a slight decrease in velocity. The velocity decreases to 0.315 m/s when the satellite is completely detached from the guide bars. The complete detachment velocity 0.315 m/s of the satellite is smaller than the simulation result 0.35 m/s. The difference of 0.035 m/s is due to the tilt of the placed spring, which leads to a shortening of the stroke. The difference in the spring stiffness coefficient causes the satellite to develop an angular velocity around the x-axis. The small bending of the guide bars caused by the gravity leads to an increase in the probability of collision occurrence, while the friction between the mechanisms also leads to energy loss. Figure 20b shows the experimental results of the separation spring under a 10% error condition. The maximum separation velocity of the satellite is 0.34 m/s. The satellite collides with the guide bars, causing a stepwise decrease in velocity. The velocity decreases to 0.252 m/s when the satellite is completely detached from the guide bars. The complete separation velocity 0.252 m/s of the satellite is smaller than the simulated value 0.278 m/s. The experimental results show that the four-guide-bar corner arrangement ensures the satellite separate smoothly and the complete separation velocity loss is small. Even under the large spring error, a reliable release can still be successfully achieved.
Aerospace 2022, 9, x FOR PEER REVIEW 20 of 22 Figure 20 illustrates the test results of the separation speed. Figure 20a shows that the maximum separation velocity of the satellite is 0.345 m/s. The satellite collides with the guide bars, resulting in a slight decrease in velocity. The velocity decreases to 0.315 m/s when the satellite is completely detached from the guide bars. The complete detachment velocity 0.315 m/s of the satellite is smaller than the simulation result 0.35 m/s. The difference of 0.035 m/s is due to the tilt of the placed spring, which leads to a shortening of the stroke. The difference in the spring stiffness coefficient causes the satellite to develop an angular velocity around the x-axis. The small bending of the guide bars caused by the gravity leads to an increase in the probability of collision occurrence, while the friction between the mechanisms also leads to energy loss. Figure 20b shows the experimental results of the separation spring under a 10% error condition. The maximum separation velocity of the satellite is 0.34 m/s. The satellite collides with the guide bars, causing a stepwise decrease in velocity. The velocity decreases to 0.252 m/s when the satellite is completely detached from the guide bars. The complete separation velocity 0.252 m/s of the satellite is smaller than the simulated value 0.278 m/s. The experimental results show that the four-guide-bar corner arrangement ensures the satellite separate smoothly and the complete separation velocity loss is small. Even under the large spring error, a reliable release can still be successfully achieved. A large angular velocity of separation causes satellite spin and brings adverse effects on the subsequent missions. Considering the case of a 10% error in the spring stiffness coefficient, the separation experiment results are shown in Figure 21. The additional disturbance moment from the spring inconsistency makes the satellite collide several times during the separation process. The angular velocity decreases after each collision. A large angular velocity of separation causes satellite spin and brings adverse effects on the subsequent missions. Considering the case of a 10% error in the spring stiffness coefficient, the separation experiment results are shown in Figure 21. The additional disturbance moment from the spring inconsistency makes the satellite collide several times during the separation process. The angular velocity decreases after each collision. The trend of the experimental results and simulation results are consistent. The angular velocity obtained by dynamics simulation is 1.88 • /s. The angular velocity at the moment of satellite detaching from the guide bars in the experimental results is 1.34 • /s. This shows that the separation angular velocity is reduced due to the contact collision between the guide bars and the satellite.

Conclusions
According to multi-satellite launch requirements and mission targets, a comp and release scheme with four guide bars arranged at corners is proposed. A ne pression and release device is designed correspondingly, which achieves good co sion and separation characteristics. In the frequency interval of 5~100 Hz, under tion of sinusoidal vibration, the maximum responses in the x, y, and z directions g, 4.29 g, and 0.68 g, respectively, indicating that the compression and release dev be reliably locked under the given vibration conditions. The maximum response x, y, and z directions were 0.95 g 2 /Hz, 0.81 g 2 /Hz, and 0.09 g 2 /Hz, respectively, un action of random vibration of 10 to 2000 Hz, reflecting good compression perfo The average time of unlocking and releasing is only about 312 ms. This verifies compression and separation device can still unlock and fast release after experien complex vibration. The maximum velocity of the satellite release process is 0.345 m the complete separation velocity is 0.315 m/s. Under 10% spring error, the ma satellite separation velocity is 0.34 m/s, the complete release velocity is reduced m/s, and the separation angle is 1.34°/s. The results show that the configuratio ranging guide bars at the four corners makes the satellite separate smoothly, complete separation velocity loss is small. Reliable release and low interference velocity are successfully achieved, even with large spring errors.  The separation experimental results show that both separation velocity loss of the satellite and the disturbance angular velocity are small under the four-guide-bar corner arrangement. Thus, the good separation characteristics are achieved successfully.

Conclusions
According to multi-satellite launch requirements and mission targets, a compression and release scheme with four guide bars arranged at corners is proposed. A new compression and release device is designed correspondingly, which achieves good compression and separation characteristics. In the frequency interval of 5~100 Hz, under the action of sinusoidal vibration, the maximum responses in the x, y, and z directions are 4.53 g, 4.29 g, and 0.68 g, respectively, indicating that the compression and release device can be reliably locked under the given vibration conditions. The maximum responses in the x, y, and z directions were 0.95 g 2 /Hz, 0.81 g 2 /Hz, and 0.09 g 2 /Hz, respectively, under the action of random vibration of 10 to 2000 Hz, reflecting good compression performance. The average time of unlocking and releasing is only about 312 ms. This verifies that the compression and separation device can still unlock and fast release after experiencing the complex vibration. The maximum velocity of the satellite release process is 0.345 m/s, and the complete separation velocity is 0.315 m/s. Under 10% spring error, the maximum satellite separation velocity is 0.34 m/s, the complete release velocity is reduced to 0.252 m/s, and the separation angle is 1.34 • /s. The results show that the configuration of arranging guide bars at the four corners makes the satellite separate smoothly, and the complete separation velocity loss is small. Reliable release and low interference angular velocity are successfully achieved, even with large spring errors.