Spin spacecraft dynamics modeling method with mass mutation characteristic

Abstract

The invention discloses a spinning spacecraft dynamics modeling method with mass mutation characteristics. The spinning spacecraft dynamics modeling method comprises the following steps that 100, a whole satellite system is decomposed into a satellite platform body, a three-axis superconducting magnet system, a rotation separation device, a momentum exchange device and a separable load according to the structural characteristics of the whole satellite system; and step 200, based on the orbit and structural characteristics of the satellite, proposing a derivation hypothesis about a whole-satellite kinetic equation so as to construct an ideal modeling environment. According to the dynamic modeling method for the partial spinning spacecrafts, mass mutation is considered, the modeling method can be used for on-orbit delivery spacecrafts and other spacecrafts with mass mutation characteristics, a vector mechanics method is adopted in the modeling process, the dynamic state of a system can be reflected more visually, and the modeling accuracy is improved. And in-depth research of the satellite nutation mechanism from the overall perspective is facilitated.

Description

Technical Field

[0001] The invention relates to the field of aerospace technology, and in particular to a method for dynamic modeling of a spinning spacecraft with mass mutation characteristics. Background Art

[0002] With the huge increase in human's demand for exploring and utilizing near-Earth space, traditional single-spin and dual-spin spacecraft are difficult to achieve established service goals such as continuous and high-speed in-orbit delivery of multiple payloads. Therefore, some spinning spacecraft came into being.

[0003] Some spinning spacecraft have important advantages such as flexibility, controllability and stability, and can carry out more complex or special on-orbit missions, making their position in space resource utilization increasingly prominent. Among them, the space target transportation based on the space rotation momentum exchange tether proposed by NASA and the low-orbit in-situ geomagnetic energy storage delivery based on the space-based rotation acceleration principle proposed by the Institute of Mechanics of the Chinese Academy of Sciences (patents 201900773631.X and 201910774225.5) are typical representatives. In most cases, such partially spinning spacecraft also have mass mutation characteristics, which brings great challenges to the dynamic modeling of the entire system.

[0004] At present, the available literature and authorized patents (represented by the patent with publication number CN109033604A) mainly provide system dynamics modeling methods for satellites with rotating loads or elastic / flexible attachments. The dynamic equations can only study the stress conditions of local single-machine components, and the modeling does not consider the mass mutation process. It is impossible to study the nutation instability mechanism and characteristics of the satellite platform caused by this feature. Without mastering the nutation law, it is impossible to propose effective active / passive suppression control laws and design controllers. Therefore, the dynamic equations of partially spinning spacecraft with mass mutation characteristics are crucial, and the modeling method needs to be studied urgently. This work will lay a theoretical foundation for the development of special spacecraft dynamics and control. Summary of the invention

[0005] The purpose of the present invention is to provide a method for modeling the dynamics of a spinning spacecraft with a mass mutation characteristic, so as to solve the technical problem that the modeling method of the dynamic equations of some spinning spacecraft with a mass mutation characteristic in the prior art needs to be studied urgently.

[0006] In order to solve the above technical problems, the present invention specifically provides the following technical solutions:

[0007] A method for modeling the dynamics of a spinning spacecraft with mass mutation characteristics comprises the following steps:

[0008] Step 100: Decompose the entire satellite system into a satellite platform body, a three-axis superconducting magnet system, a rotation separation device, a momentum exchange device, and a detachable payload according to the structural characteristics of the entire satellite system;

[0009] Step 200: Based on the orbit and structural characteristics of the satellite, a hypothesis for the derivation of the dynamic equation of the entire satellite is proposed to construct an ideal modeling environment;

[0010] Step 300: To facilitate the description of each component, the following coordinate system is established:

[0011] The system center-of-mass orbit coordinate system OXYZ with the system's center of mass as its origin;

[0012] ObXbYbZb with the center of mass of the satellite platform as the origin;

[0013] OcXcYcZc with the center of mass of the rotating separation device as the origin;

[0014] OdXdYdZd with the centroid of the detachable load as the origin;

[0015] Step 400, combined with the derivation assumptions of step 200, the system after load separation is regarded as a two-rigid body system, and the system angular momentum is expressed as four parts, namely the angular momentum generated by the rotating separation device around its own center of mass, the angular momentum generated by the rotating separation device center of mass around the system center of mass, the angular momentum generated by the satellite platform body rotating around its own center of mass, and the angular momentum generated by the satellite platform body center of mass around the system center of mass, which are defined as And respectively calculate;

[0016] Calculate the combined angular momentum of the above four systems;

[0017] Transforming the angular momentum based on the attitude conversion matrix between the reference system of the rotating separation device and the reference system of the satellite platform body to obtain the angular momentum of the system relative to the center of mass;

[0018] Step 500: Use the angular momentum theorem to calculate the system's relative center of mass angular momentum obtained in step 400, and combine the derived assumptions to obtain the overall dynamic equation of the satellite system.

[0019] As a preferred solution of the present invention, in step 100;

[0020] In step 100;

[0021] The satellite platform body is composed of a rigid shell, internal rigid components and a momentum flywheel;

[0022] The momentum exchange device is composed of a stator and a rotor, and the stator is fixedly connected to the satellite platform body;

[0023] The three-axis superconducting magnet system is fixedly connected to the satellite platform body through the stator, and the three-axis superconducting magnet system can adjust the attitude of the satellite platform body through external torque;

[0024] The rotating separation device is fixedly mounted on the rotor, a slide rail for continuously conveying the separable load is arranged on the rotating separation device, and a controllable instantaneous release mechanism is arranged at the end of the slide rail;

[0025] The detachable load is a single rigid body slidably mounted on the slide rail, and the detachable load is launched through the controllable instantaneous release mechanism.

[0026] As a preferred solution of the present invention, the following assumptions are derived in step 200:

[0027] In the dynamic modeling and analysis, the satellite platform body, mass separation and rotation separation device, and superconducting magnet system are all assumed to be rigid bodies;

[0028] Ignore the mass of the momentum exchange device, the gap between the stator and rotor, the damping, etc.

[0029] Assume that the mass separation rotating separation device and the load to be separated are separated instantaneously without impact between them;

[0030] It is assumed that the satellite is in a stable state before separation.

[0031] As a preferred solution of the present invention, the specific steps of establishing multiple coordinate systems in step 300 are as follows:

[0032] System center of mass orbit coordinate system OXYZ: The origin O is at the center of mass of the system, the OZ axis points to the center of the earth, the OX axis is in the satellite orbit plane, perpendicular to the OZ axis and points to the direction of satellite motion, and the OY axis, OX and OZ axes form a right-handed coordinate system;

[0033] The satellite platform body coordinate system ObXbYbZb: origin O b Located at the mass center of the satellite platform, when the three azimuth angles of the platform relative to the orbital coordinate system are zero, the directions of the axes are consistent with the directions of the corresponding axes of the orbital coordinate system OXYZ;

[0034] The rotating separation device body coordinate system OcXcYcZc: the origin of the coordinate system O c Located at the intersection of the rotating axis of the rotating separation device and the rotating plane, it is assumed that at the initial moment, the axes of the rotating separation device and the coordinate axes of the platform body coordinate system are in the same direction;

[0035] The detachable load body coordinate system OdXdYdZd: the origin of the coordinate system O d Located at the centroid of the detachable load 5, the three coordinate axes before separation are respectively parallel to the corresponding coordinate axes of the coordinate system of the rotating separation device body.

[0036] As a preferred solution of the present invention, in step 400, after the detachable load is separated, the angular momentum generated by the rotating separation device around its own center of mass is:

[0037]

[0038] in; The moment of inertia of the rotating separation device and the delivery load relative to the combined center of mass before separation;

[0039] ω is the relative angular velocity between the satellite platform body reference system and the center of mass orbit reference system;

[0040] Ω is the relative angular velocity between the optional device reference system and the satellite platform system.

[0041] As a preferred embodiment of the present invention, in step 400, the centroid of the rotating separation device is rotated around the system

[0042] The angular momentum generated by the combined center of mass is:

[0043]

[0044] Among them, m c is the mass of the rotating separation device;

[0045] r c is the radius vector of the center of mass of the rotating separation device relative to the total center of mass of the system;

[0046] ω is the relative angular velocity between the satellite platform body reference system and the center of mass orbit reference system;

[0047] Ω is the relative angular velocity between the optional device reference system and the satellite platform system.

[0048] As a preferred solution of the present invention, in step 400, the angular momentum generated by the satellite platform body rotating around its own center of mass is:

[0049]

[0050] in, The moment of inertia of the satellite platform body and the superconducting device relative to the combined center of mass;

[0051] ω is the relative angular velocity between the satellite platform body reference system and the center of mass orbit reference system.

[0052] As a preferred embodiment of the present invention, in step 400, the platform mass center is generated around the system mass center.

[0053] The resulting angular momentum:

[0054]

[0055] Among them, mb Combine the mass of the satellite platform and superconducting device;

[0056] r b is the radius vector of the satellite platform body center of mass relative to the system center of mass;

[0057] ω is the relative angular velocity between the satellite platform body reference system and the center of mass orbit reference system;

[0058] Ω is the relative angular velocity between the optional device reference system and the satellite platform system.

[0059] As a preferred solution of the present invention, in step 400, the calculation time of the rotational inertia of the center of mass of the rotating separation device relative to the total center of mass of the system is the moment when the detachable load is separated from the rotating separation device, and it is assumed that there is no impact between the detachable load and the rotating separation device.

[0060] Compared with the prior art, the present invention has the following beneficial effects:

[0061] The dynamic modeling method of a partially spinning spacecraft in the present invention takes mass mutation into consideration. The modeling method can be used not only for in-orbit delivery spacecraft, but also for other spacecraft with mass mutation characteristics. In addition, the present invention adopts vector mechanics in the modeling process, which can more intuitively reflect the dynamic state of the system and is conducive to in-depth research on the satellite oscillation mechanism from a holistic perspective. BRIEF DESCRIPTION OF THE DRAWINGS

[0062] In order to more clearly illustrate the implementation methods of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for the implementation methods or the description of the prior art. Obviously, the drawings in the following description are only exemplary, and for ordinary technicians in this field, other implementation drawings can be derived from the provided drawings without creative work.

[0063] Figure 1 A structural schematic diagram of a method for modeling the dynamics of a spinning spacecraft with mass mutation characteristics is provided for an embodiment of the present invention;

[0064] The numbers in the figure represent the following:

[0065] 1 Satellite platform body, 2 Three-axis superconducting magnet system;

[0066] 3 rotating separation device, 301 slide rail;

[0067] 4 momentum exchange device, 401 stator, 402 rotor;

[0068] 5 Separable load. DETAILED DESCRIPTION

[0069] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0070] like Figure 1 As shown, the present invention provides a method for dynamic modeling of a spinning spacecraft with mass mutation characteristics, comprising the following steps:

[0071] Step 100: Decompose the entire satellite system into a satellite platform body 1, a three-axis superconducting magnet system 2, a rotation separation device 3, a momentum exchange device 4 and a detachable payload 5 according to the structural characteristics of the entire satellite system.

[0072] Wherein: the satellite platform body 1 is composed of a rigid outer shell, internal rigid components and a momentum flywheel.

[0073] The momentum exchange device 4 is composed of a stator 401 and a rotor 402 , and the stator 401 is fixedly connected to the satellite platform body 1 , and the relative rotation direction and speed of the stator 401 and the rotor 402 are controlled by the system.

[0074] The three-axis superconducting magnet system is fixedly connected to the satellite platform body 1 through the stator 401, and the three-axis superconducting magnet system can obtain momentum by interacting with the geomagnetic field to inject energy into the system, and can also obtain external torque for adjusting the platform attitude or flywheel unloading.

[0075] The rotating separation device 3 is fixedly mounted on the rotor 402 . A slide rail 301 for continuously conveying the separable load 5 is arranged on the rotating separation device 3 , and a controllable instantaneous release mechanism is arranged at the end of the slide rail 301 .

[0076] The detachable load 5 is a single rigid body slidably mounted on the slide rail 301 , and the detachable load 5 is launched through a controllable instantaneous release mechanism.

[0077] Step 200: Based on the orbit and structural characteristics of the satellite, a hypothesis for the derivation of the dynamic equation of the whole satellite is proposed;

[0078] 1. In the dynamic modeling and analysis, the satellite platform body 1, the rotating separation device 3, and the three-axis superconducting magnet system 2 are all assumed to be rigid bodies.

[0079] 2. Ignore the mass of the momentum exchange device 4, and the gap between the stator 401 and the rotor 402, damping and other effects.

[0080] 3. It is assumed that the rotating separation device 3 and the separable load 5 are separated instantaneously without any impact between them.

[0081] 4. Assume that the satellite is in a stable state before separation.

[0082] Step 300: In order to facilitate the description of each component, the following coordinate system is introduced:

[0083] System center of mass orbit coordinate system OXYZ: The origin O is at the center of mass of the system, the OZ axis points to the center of the earth, the OX axis is in the satellite orbit plane, perpendicular to the OZ axis and points to the direction of satellite motion, and the OY axis, OX and OZ axes form a right-handed coordinate system;

[0084] Satellite platform body 1 coordinate system ObXbYbZb: The origin Ob is located at the mass center of the satellite platform body 1. When the three azimuth angles of the platform relative to the orbital coordinate system are zero, the directions of the axes are consistent with the directions of the corresponding axes of the orbital coordinate system OXYZ.

[0085] Coordinate system OcXcYcZc of the rotating separation device 3: the origin Oc of the coordinate system is located at the intersection of the rotating axis of the rotating separation device 3 and the rotating plane, and it is assumed that at the initial moment, the axes of the rotating separation device 3 point in the same direction as the coordinate axes of the platform body coordinate system;

[0086] Coordinate system of the detachable load 5 body OdXdYdZd: The origin Od of the coordinate system is located at the center of mass of the detachable load 5, and the three coordinate axes before separation are respectively parallel to the corresponding coordinate axes of the coordinate system of the rotating separation device 3 body.

[0087] Step 400, based on the assumptions derived in step 200, the system after load separation is regarded as a two-rigid body system (basic assumption 1 in step 200), and the angular momentum of the system is expressed as four parts, namely, the angular momentum generated by the rotating separation device around its own center of mass, the angular momentum generated by the rotating separation device center of mass around the system center of mass, the angular momentum generated by the platform rotating around its own center of mass, and the angular momentum generated by the platform center of mass around the system center of mass, which are defined as

[0088] The following are the meanings of the symbols in the following formulas:

[0089] Symbols

[0090] d Distance between the load and the center of mass of the rotating separation device

[0091] m c Rotating separation device mass

[0092] m b Satellite platform and superconducting device quality

[0093] Δm Mass of load to be separated

[0094] h Distance between the rotation plane of the optional device and the center of mass of the platform

[0095] m1 Mass of load to be separated

[0096] m2 Mass of load to be separated

[0097] e Mass of load to be separated

[0098] u Mass of load to be separated

[0099] ω Relative angular velocity between the satellite platform body reference system and the center of mass orbit reference system

[0100] Ω Relative angular velocity between the optional device reference system and the satellite platform system r c Radius vector of the center of mass of the rotating separation device relative to the total center of mass of the system

[0101] r b The radius vector of the satellite platform's center of mass relative to the system's combined center of mass

[0102] Moment of inertia of load and rotating separation device relative to combined center of mass before separation

[0103] The moment of inertia of the rotating separation device relative to its own center of mass

[0104] The moment of inertia of the satellite platform and superconducting device relative to the combined center of mass

[0105] The moment of inertia of the load relative to the center of mass of the optional device

[0106] Moment of inertia of the center of mass of the rotating separation device relative to the total center of mass of the system

[0107] C The attitude transformation matrix from the reference system of the rotating separation device to the reference system of the platform body

[0108] L Control torque generated by superconducting magnetic moment and geomagnetic interaction

[0109] Among them, 1. The calculation process is as follows:

[0110] Assume that the inertia matrix of the delivered load relative to its own center of mass is The distance from the center of mass c of the rotating separation device is d. Then, in the reference system of the rotating separation device, the moment of inertia of the delivery load relative to the center of mass c of the rotating separation device can be expressed as:

[0111] in,

[0112] Formula (1)

[0113] For a rotating separation device, the center of mass offset due to load separation is:

[0114] Formula (2)

[0115] The inertia of the rotating separation device after mass separation can be expressed as:

[0116] Formula (3)

[0117] in, It is the moment of inertia of the rotating separation device and the delivery load relative to the combined center of mass before separation.

[0118] After the load is separated, the angular momentum generated by the rotating separation device around its own center of mass is:

[0119] Formula (4)

[0120] Similarly, the calculation results of the angular momentum generated by the center of mass of the rotating separation device around the combined center of mass of the system, the angular momentum generated by the platform rotating around its own center of mass, and the angular momentum generated by the center of mass of the platform around the combined center of mass of the system are as follows:

[0121] 2. Calculation

[0122] The angular momentum generated by the center of mass of the rotating separation device around the combined center of mass of the system:

[0123] Formula (5)

[0124] 3. Calculation

[0125] The angular momentum generated by the platform rotating around its own center of mass:

[0126] Formula (6)

[0127] 4. Calculation

[0128] The angular momentum generated by the center of mass of the platform around the combined center of mass of the system:

[0129] Formula (7)

[0130] 5. Calculation and simplification of angular momentum

[0131] The combined angular momentum can be expressed as: Formula (8)

[0132] in,

[0133]

[0134] is the moment of inertia of the center of mass of the rotating part relative to the combined center of mass at t = 0;

[0135] C is the attitude transformation matrix between the reference system of the rotating separation device and the reference system of the platform body.

[0136] Since the load separation occurs instantaneously, there is no impact between the load and the rotating separation device (basic assumption 3 of step 200), and:

[0137] Formula (9)

[0138] Where: ΔH is the angular momentum loss due to load separation, expressed as:

[0139]

[0140] The initial kinetic energy of the system after load separation is:

[0141] Formula (10) E(t=0)=E0-ΔE;

[0142] where is the mechanical energy loss due to load separation.

[0143]

[0144] Step 500: Use the angular momentum theorem to calculate the system relative mass center angular momentum obtained by formula (9), and obtain:

[0145] Formula (11)

[0146]

[0147] Where: M is the sum of the interference torque and control torque acting on the satellite.

[0148] set up Simplifying, we get:

[0149] Formula (12)

[0150] This equation is the overall dynamic equation of the satellite system.

[0151] The above embodiments are only exemplary embodiments of the present application and are not intended to limit the present application. The protection scope of the present application is defined by the claims. Those skilled in the art may make various modifications or equivalent substitutions to the present application within the essence and protection scope of the present application, and such modifications or equivalent substitutions shall also be deemed to fall within the protection scope of the present application.

Claims

1. A method for dynamic modeling of a spinning spacecraft with mass mutation characteristics, characterized in that: The following steps are involved: Step 100: Decompose the entire satellite system into a satellite platform body (1), a three-axis superconducting magnet system (2), a rotation separation device (3), a momentum exchange device (4), and a detachable payload (5) according to the structural characteristics of the entire satellite system; Step 200: Based on the orbit and structural characteristics of the satellite, a hypothesis for the derivation of the dynamic equation of the entire satellite is proposed to construct an ideal modeling environment; Step 300: To facilitate the description of each component, the following coordinate system is established: The system center-of-mass orbit coordinate system OXYZ with the system's center of mass as its origin; ObXbYbZb with the center of mass of the satellite platform body (1) as the origin; OcXcYcZc with the center of mass of the rotating separation device (3) as the origin; OdXdYdZd with the center of mass of the detachable load (5) as the origin; Step 400, combining the derivation assumptions of step 200, the system after load separation is regarded as a two-rigid body system, and the angular momentum of the system is expressed as four parts, namely, the angular momentum generated by the rotating separation device (3) around its own center of mass, the angular momentum generated by the center of mass of the rotating separation device (3) around the combined center of mass of the system, the angular momentum generated by the satellite platform body (1) rotating around its own center of mass, and the angular momentum generated by the center of mass of the satellite platform body (1) around the combined center of mass of the system, which are respectively defined as And respectively calculate; Calculate the combined angular momentum of the above four systems; Based on the attitude conversion matrix between the reference system of the rotating separation device (3) and the reference system of the satellite platform body (1), the combined angular momentum is transformed to obtain the combined angular momentum of the system relative to the center of mass; Step 500: Use the angular momentum theorem to calculate the system's relative center of mass angular momentum obtained in step 400, and combine the derived assumptions to obtain the overall dynamic equation of the satellite system.

2. A method for dynamic modeling of a spinning spacecraft with mass mutation characteristics according to claim 1, characterized in that: In step 100; The satellite platform body (1) is composed of a rigid outer shell, internal rigid components and a momentum flywheel; The momentum exchange device (4) is composed of a stator (401) and a rotor (402), and the stator (401) is fixedly connected to the satellite platform body (1); The three-axis superconducting magnet system is fixedly connected to the satellite platform body (1) via the stator (401), and the three-axis superconducting magnet system can adjust the attitude of the satellite platform body (1) through an external torque; The rotating separation device (3) is fixedly mounted on the rotor (402), a slide rail (301) for continuously conveying the separable load (5) is arranged on the rotating separation device (3), and a controllable instantaneous release mechanism is arranged at the end of the slide rail (301); The detachable load (5) is a single rigid body slidably mounted on the slide rail (301), and the detachable load (5) is launched through the controllable instantaneous release mechanism.

3. A method for dynamic modeling of a spinning spacecraft with mass mutation characteristics according to claim 2, characterized in that: The derivation assumptions in step 200 are as follows: In dynamic modeling and analysis, the satellite platform body (1), the rotating separation device (3), and the superconducting magnet system (2) are all assumed to be rigid bodies; Ignore the mass of the momentum exchange device (4), and the effects of the gap and damping between the stator (401) and the rotor (402); Assume that the mass rotating separation device (3) and the load to be separated are separated instantaneously without any impact between them; It is assumed that the satellite is in a stable state before separation.

4. The method for dynamic modeling of a spinning spacecraft with mass mutation characteristics according to claim 2, characterized in that: The specific steps of establishing multiple coordinate systems in step 300 are as follows: System center of mass orbit coordinate system OXYZ: The origin O is at the center of mass of the system, the OZ axis points to the center of the earth, the OX axis is in the satellite orbit plane, perpendicular to the OZ axis and points to the direction of satellite motion, and the OY axis, OX and OZ axes form a right-handed coordinate system; The satellite platform body (1) has a coordinate system ObXbYbZb: origin O b Located at the mass center of the satellite platform body (1), when the three azimuth angles of the platform relative to the orbital coordinate system are zero, the directions of the axes are consistent with the directions of the corresponding axes of the orbital coordinate system OXYZ; The rotating separation device (3) has a body coordinate system OcXcYcZc: the origin of the coordinate system O c Located at the intersection of the rotation axis of the rotating separation device (3) and the rotating plane, it is assumed that at the initial moment, the directions of the axes of the rotating separation device (3) and the coordinate axes of the platform body coordinate system are consistent; The detachable load (5) has a body coordinate system OdXdYdZd: the origin of the coordinate system O d Located at the centroid of the detachable load 5, the three coordinate axes before separation are respectively parallel to the corresponding coordinate axes of the main body coordinate system of the rotating separation device (3).

5. A method for dynamic modeling of a spinning spacecraft with mass mutation characteristics according to claim 4, characterized in that: In step 400, after the detachable load is separated (5), the angular momentum generated by the rotating separation device (3) around its own center of mass is: in; The moment of inertia of the rotating separation device (3) and the delivery load relative to the combined center of mass before separation; ω is the relative angular velocity between the satellite platform body reference system and the center of mass orbit reference system; Ω is the relative angular velocity between the optional device reference system and the satellite platform system.

6. The method for dynamic modeling of a spinning spacecraft with mass mutation characteristics according to claim 4, characterized in that: In step 400, the angular momentum generated by the center of mass of the rotating separation device (3) around the combined center of mass of the system is: Among them, m c is the mass of the rotating separation device; r c is the radius vector of the center of mass of the rotating separation device (3) relative to the total center of mass of the system; ω is the relative angular velocity between the satellite platform body (1) reference system and the center of mass orbit reference system; Ω is the relative angular velocity between the optional device reference system and the satellite platform body (1) system.

7. The method for dynamic modeling of a spinning spacecraft with mass mutation characteristics according to claim 4, characterized in that: In step 400, the angular momentum generated by the satellite platform body (1) rotating around its own center of mass is: in, The moment of inertia of the satellite platform body (1) and the superconducting device relative to the combined center of mass; ω is the relative angular velocity between the satellite platform body (1) reference frame and the center of mass orbit reference frame.

8. The method for dynamic modeling of a spinning spacecraft with mass mutation characteristics according to claim 4, characterized in that: In step 400, the angular momentum generated by the center of mass of the platform around the combined center of mass of the system is: Among them, m b The satellite platform body (1) and the superconducting device have a combined mass; r b is the radius vector of the mass center of the satellite platform body (1) relative to the system mass center; ω is the relative angular velocity between the satellite platform body (1) reference system and the center of mass orbit reference system; Ω is the relative angular velocity between the optional device reference system and the satellite platform body (1) system.

9. The method for dynamic modeling of a spinning spacecraft with mass mutation characteristics according to claim 1, characterized in that: In step 400, the moment of inertia of the center of mass of the rotating separation device (3) relative to the total center of mass of the system is calculated at the moment when the detachable load (5) and the rotating separation device (3) are separated, and it is assumed that there is no impact between the detachable load (5) and the rotating separation device (3).

Images