CN112660422A - Split type satellite on-orbit platform cabin quality identification method and system - Google Patents

Split type satellite on-orbit platform cabin quality identification method and system Download PDF

Info

Publication number
CN112660422A
CN112660422A CN202011475974.7A CN202011475974A CN112660422A CN 112660422 A CN112660422 A CN 112660422A CN 202011475974 A CN202011475974 A CN 202011475974A CN 112660422 A CN112660422 A CN 112660422A
Authority
CN
China
Prior art keywords
cabins
split
satellite
platform cabin
mass
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011475974.7A
Other languages
Chinese (zh)
Other versions
CN112660422B (en
Inventor
洪振强
边志强
赵艳彬
沈毅力
张健
李文峰
陈昌亚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Institute of Satellite Engineering
Original Assignee
Shanghai Institute of Satellite Engineering
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Institute of Satellite Engineering filed Critical Shanghai Institute of Satellite Engineering
Priority to CN202011475974.7A priority Critical patent/CN112660422B/en
Publication of CN112660422A publication Critical patent/CN112660422A/en
Application granted granted Critical
Publication of CN112660422B publication Critical patent/CN112660422B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

The invention provides a quality identification method and a system for an on-orbit platform cabin of a split satellite, which are used for identifying the quality of the platform cabin after two cabins of the satellite are separated. The identification method specifically comprises the following steps: powering off a repeated locking mechanism between two cabins of the split satellite, electrifying a magnetic floating actuator between the two cabins, generating magnetic control force with known magnitude, and driving the two cabins to move; the control system keeps open loop uncontrolled, collects the information measured by the displacement sensor of the two cabins, and obtains the relative movement information of the mass center of the two cabins according to the resolving model; carrying out secondary fitting on the relative movement information of the mass center, and differentiating the fitting result twice to obtain the acceleration information of the mass center of the two cabins; and resolving to obtain the platform cabin mass by combining the known inter-cabin acting force, the acceleration information and the load cabin mass information. The method can provide technical support for the on-orbit propellant residual quantity evaluation and relative position control of the split satellite.

Description

Split type satellite on-orbit platform cabin quality identification method and system
Technical Field
The invention relates to a composite control technology of a loading cabin of a satellite platform with ultrahigh pointing accuracy and ultrahigh stability (double-super), in particular to a quality identification method and system of a split type satellite on-orbit platform cabin.
Background
The requirements of the advanced spacecraft on the attitude pointing accuracy and the stability are two orders of magnitude higher than the current level, the traditional method adopts a load and platform fixed connection type design, the dynamic characteristics of the two are deeply coupled, so that the load double super index is difficult to realize, and although the methods of active and passive micro vibration suppression and the like are adopted to obtain a certain effect, the defect of the limited fixed connection type design is that the double super index is difficult to realize.
The 'double-super' satellite platform breaks through the traditional fixed connection design, adopts a non-contact, high-precision and time-delay-free displacement sensor to realize the separation of a load cabin only provided with a quiet component and a platform cabin provided with a movable component, and thoroughly eliminates the micro-vibration influence. The traditional control logic mainly based on a satellite platform is changed, and a brand new method of 'load cabin driving, platform cabin driven and two cabin relative positions cooperative decoupling control' is adopted for the first time, so that the double super-precision of the load cabin can be realized.
The accurate identification of the quality of the two cabins is a precondition for realizing high-precision position control between the two cabins. The split satellite comprises a load cabin and a platform cabin, wherein the load cabin mainly comprises a load and a supporting structure thereof. The on-orbit quality is not changed, and the on-orbit quality can be accurately measured on the ground. The pod contains propellant and will be constantly consumed on-track, so that the pod mass is variable and needs to be identified on-track. The two-chamber mass is an important input in the relative position feed-forward control. It is therefore necessary to obtain accurate two-compartment mass position information by means of on-track calibration.
Patent document CN106249749A discloses a master-slave non-contact dual-super satellite platform variable-centroid and variable-inertia attitude control system, which does not solve the technical problem of accurate identification of two-compartment rotational inertia position.
Patent document CN109870272A discloses a method for identifying the mass of a spacecraft in orbit based on momentum conservation, but the object of the method is an integrated satellite.
In the document [1] [2] [3] [4], a method for identifying the quality characteristics of the on-orbit spacecraft based on different algorithms is proposed, but all research objects of the method are all integrated satellites, and the method is not related to and is not suitable for identifying the quality characteristics of the split satellites in a separated state.
[1]Bergmann EV,Walker BK,Levy D R.Mass property estimation for control of asymmetrical satellites[J].Journal of Guidance,Control and Dynamics,1987,10(2):483-492.
[2]Bergmann E V,Dzielski J.Spacecraft mass property identification with torque-generating control[J].Journal of Guidance,Control,and Dynamics,1990,13(2):99-103.
[3]Wilson E,Lages C,Mah R.On-line,gyro-based,mass-property identification for thruster-controlled spacecraft using recursive least squares[C]//Proceedings of the 45th Midwest Symposium on Circuits and Systems.Moffett Field,California,Ames Research Center,Aug.4-7,2002
[4]Tanygin S,Williams T.Mass property estimation using coasting maneuvers[J].Journal of Guidance,Control,and Dynamics,1997,20(4):625-632
Disclosure of Invention
Aiming at the requirement of high-precision attitude control of a load cabin of a two-cabin non-contact type double-super satellite platform, the invention aims to provide a novel split type satellite on-orbit platform cabin quality identification method and system.
According to one aspect of the invention, a novel split type satellite on-orbit platform cabin quality identification method is provided, and the method comprises the following steps:
a driving step: powering off a repeated locking mechanism between two cabins of the split satellite, electrifying a magnetic floating actuator between the two cabins, generating magnetic control force with known magnitude, and driving the two cabins of the split satellite to move;
the collection step comprises: the control system of the split satellite is kept open-loop and uncontrolled, the measurement information of the displacement sensors of the two cabins of the split satellite is collected, and the relative motion information of the mass center of the two cabins of the split satellite is obtained according to the resolving model;
and fitting differentiation: carrying out secondary fitting on the relative motion information of the mass centers of the two cabins of the split satellite, and differentiating the fitting result twice to obtain the acceleration information of the mass centers of the two cabins of the split satellite;
a platform cabin quality resolving step: and resolving to obtain the platform cabin quality by combining the known acting force and acceleration information between the two cabins of the split satellite and the load cabin quality information.
Preferably, in the driving step, a plurality of magnetic actuators are arranged between the two cabins and used for finishing the three-axis attitude control of the load cabin and the relative position control between the two cabins, an open-loop instruction is sent to enable the magnetic actuators to be combined to output magnetic control force between the cabins, and the magnetic control force acting on the platform cabin is set to be FcbThe magnetic control force acting on the load compartment is-FcbUnder the action of magnetic control force between the cabins, the acceleration of the platform cabin is abAcceleration of the load compartment of apLet the mass of the load compartment be mpMass m of the platform cabinbSo, according to Newton's second law of motion, there is
Figure BDA0002837349670000031
The relative distance between the interstitial centers of the two cabins is PCJAnd if the magnetic control force between the two cabins drives the distance between the two cabins to be larger, the distance between the two cabins is larger
Figure BDA0002837349670000032
Wherein t is time.
Preferably, in the acquiring step, the three-axis relative attitude between the two cabins is set as
Figure BDA0002837349670000034
Three-axis relative centroid displacement of [ Delta x Delta y Delta z]Two compartments were provided with 9 displacement sensors A1, A2, A3, B1, B2, B3, C1, C2 and C3 for relative position measurement, and the measured values were [ Δ z ] respectivelyA1 ΔzA2 ΔzA3 ΔyB1 ΔyB2 ΔyB3 ΔxC1 ΔxC2 ΔxC3]The installation position of 9 displacement sensors in the mechanical coordinate system is set as [ x ]i yi zi],i=A1,…C3The center of mass of the load compartment is set in the mechanical coordinate system as xpc ypc zpc]Then, the relationship between the relative position and attitude and the measured value of the displacement sensor is as follows:
Figure BDA0002837349670000033
is provided with
Figure BDA0002837349670000041
B=[ΔzA1 ΔzA2 ΔzA3 ΔyB1 ΔyB2 ΔyB3 ΔxC1 ΔxC2 ΔxC3]T Equation 5
Figure BDA0002837349670000043
Then there is
A, X ═ B formula 7
Obtaining a pseudo-inverse solution according to a pseudo-inverse algorithm
X=AT·(A·AT)-1B equation 8
Thereby obtaining the relative attitude and the relative position of the three axes and obtaining a resolving matrix of
M=AT·(A·AT)-1Equation 9
Obtaining the relative distance P between the two interstitial centers according to the calculation modelCJThen there is
PCJ=[Δx Δy Δz]TEquation 10
Preferably, in the step of fitting differentiation, the relative distance between the two capsule centers is set as PCJ_NHObtained by fitting
PCJ_NH=p1·t2+p2·t+p3Equation 11
Wherein p is1As second order fitting coefficient, p2As first-order fitting coefficient, p3Is a constant coefficient, and the linear acceleration of the relative motion between the two cabins is set as aCJThen there is
aCJ=2p1Equation 12.
Preferably, in the step of calculating the quality of the platform cabin, the magnetic control force F is generated between the cabinscbUnder the action of (2), the linear acceleration of the relative motion between the two compartments is expressed again as
Figure BDA0002837349670000042
In the above formula, aCJ、FcbAnd mpAre all known quantities, thus obtaining a platform cabin mass of
Figure BDA0002837349670000051
According to another aspect of the invention, a novel split type satellite on-orbit platform cabin quality identification system is provided, which comprises the following modules:
a driving module: powering off a repeated locking mechanism between two cabins of the split satellite, electrifying a magnetic floating actuator between the two cabins, generating magnetic control force with known magnitude, and driving the two cabins of the split satellite to move;
an acquisition module: the control system of the split satellite is kept open-loop and uncontrolled, the measurement information of the displacement sensors of the two cabins of the split satellite is collected, and the relative motion information of the mass center of the two cabins of the split satellite is obtained according to the resolving model;
fitting a differential module: carrying out secondary fitting on the relative motion information of the mass centers of the two cabins of the split satellite, and differentiating the fitting result twice to obtain the acceleration information of the mass centers of the two cabins of the split satellite;
the platform cabin quality resolving module: and resolving to obtain the platform cabin quality by combining the known acting force and acceleration information between the two cabins of the split satellite and the load cabin quality information.
Preferably, in the driving module, a plurality of magnetic suspension actuators are arranged between two cabins and used for finishing three-axis attitude control of the load cabin and relative position control between the two cabins, an open-loop instruction is sent to enable the magnetic suspension actuators to be combined to output magnetic control force between the cabins, and the magnetic control force acting on the platform cabin is set to be FcbThe magnetic control force acting on the load compartment is-FcbUnder the action of magnetic control force between the cabins, the acceleration of the platform cabin is abAcceleration of the load compartment of apLet the mass of the load compartment be mpMass m of the platform cabinbSo, according to Newton's second law of motion, there is
Figure BDA0002837349670000052
The relative distance between the interstitial centers of the two cabins is PCJAnd if the magnetic control force between the two cabins drives the distance between the two cabins to be larger, the distance between the two cabins is larger
Figure BDA0002837349670000053
Wherein t is time.
Preferably, in the acquisition module, the three-axis relative attitude between the two cabins is set as
Figure BDA0002837349670000054
Three-axis relative centroid displacement of [ Delta x Delta y Delta z]Two compartments are provided with 9 displacement sensors A1, A2, A3, B1, B2, B3, C1, C2 and C3 forRelative position measurement, respectively
[ΔzA1ΔzA2ΔzA3ΔyB1ΔyB2ΔyB3ΔxC1ΔxC2ΔxC3]The installation position of 9 displacement sensors in the mechanical coordinate system is set as [ x ]i yi zi],i=A1,…C3The center of mass of the load compartment is arranged in a mechanical coordinate system
[xpc ypc zpc]Then, the relationship between the relative position and attitude and the measured value of the displacement sensor is as follows:
Figure BDA0002837349670000061
is provided with
Figure BDA0002837349670000062
B=[ΔzA1ΔzA2ΔzA3ΔyB1ΔyB2ΔyB3ΔxC1ΔxC2ΔxC3]T Equation 5
Figure BDA0002837349670000063
Then there is
A, X ═ B formula 7
Obtaining a pseudo-inverse solution according to a pseudo-inverse algorithm
X=AT·(A·AT)-1B formula 8
Thereby obtaining the relative attitude and the relative position of the three axes and obtaining a resolving matrix of
M=AT·(A·AT)-1Equation 9
Obtaining the relative distance P between the two interstitial centers according to the calculation modelCJThen there is
PCJ=[ΔxΔyΔz]TEquation 10
Preferably, theIn the fitting differential module, the relative distance between the centers of two cabins is set as PCJ_NHObtained by fitting
PCJ_NH=p1·t2+p2·t+p3Equation 11
Wherein p is1As second order fitting coefficient, p2As first-order fitting coefficient, p3Is a constant coefficient, and the linear acceleration of the relative motion between the two cabins is set as aCJThen there is
aCJ=2p1Equation 12
Preferably, in the platform cabin quality resolving module, the magnetic control force F between the cabinscbUnder the action of (2), the linear acceleration of the relative motion between the two compartments is expressed again as
Figure BDA0002837349670000071
In the above formula, aCJ、FcbAnd mpAre all known quantities, thus obtaining a platform cabin mass of
Figure BDA0002837349670000072
Compared with the prior art, the invention has the following beneficial effects:
1. the method can provide technical support for the evaluation of the residual quantity of the on-orbit propellant and the control of the relative position of the split satellite;
2. the accurate identification of the rotational inertia of the two cabins is realized through the design of the on-orbit identification method of the rotational inertia of the two cabins;
3. the technical problem that accurate identification is difficult due to the influence of gravity and an unfolding component on the ground is solved.
Drawings
Other features, objects and advantages of the invention will become more apparent upon reading of the detailed description of non-limiting embodiments with reference to the following drawings:
FIG. 1 is a schematic diagram of the mass identification principle of a novel split type satellite in-orbit platform cabin;
FIG. 2 is a schematic diagram of a quality identification method for a novel split satellite in-orbit platform cabin;
FIG. 3 is a schematic view of a measurement curve of a displacement sensor between two compartments;
FIG. 4 is a schematic diagram of relative distance measurement and fitting between two bays.
Detailed Description
The present invention will be described in detail with reference to specific examples. The following examples will assist those skilled in the art in further understanding the invention, but are not intended to limit the invention in any way. It should be noted that various changes and modifications can be made by those skilled in the art without departing from the spirit of the invention, such as changes in the magnitude or direction of the thrust of the magnetic levitation actuator, changes in the installation position and orientation of the magnetic levitation actuator, etc. All falling within the scope of the present invention.
The invention provides a novel split type satellite on-orbit platform cabin quality identification method, which comprises the following four steps:
a driving step: powering off a repeated locking mechanism between two cabins of the split satellite, electrifying a magnetic floating actuator between the two cabins, generating magnetic control force with known magnitude, and driving the two cabins of the split satellite to move;
the collection step comprises: the control system of the split satellite is kept open-loop and uncontrolled, the measurement information of the displacement sensors of the two cabins of the split satellite is collected, and the relative motion information of the mass center of the two cabins of the split satellite is obtained according to the resolving model;
and fitting differentiation: carrying out secondary fitting on the relative motion information of the mass centers of the two cabins of the split satellite, and differentiating the fitting result twice to obtain the acceleration information of the mass centers of the two cabins of the split satellite;
a platform cabin quality resolving step: and resolving to obtain the platform cabin quality by combining the known acting force and acceleration information between the two cabins of the split satellite and the load cabin quality information.
In the driving step, a plurality of magnetic suspension actuators are arranged between the two cabins and used for finishing the three-axis attitude control of the load cabin and the relative position control between the two cabins, and an open-loop instruction is sent to enable the magnetic suspension actuators to be combined and output between the cabinsThe magnetic control force is set as FcbThe magnetic control force acting on the load compartment is-FcbUnder the action of magnetic control force between the cabins, the acceleration of the platform cabin is abAcceleration of the load compartment of apLet the mass of the load compartment be mpMass m of the platform cabinbSo, according to Newton's second law of motion, there is
Figure BDA0002837349670000081
The relative distance between the interstitial centers of the two cabins is PCJAnd if the magnetic control force between the two cabins drives the distance between the two cabins to be larger, the distance between the two cabins is larger
Figure BDA0002837349670000082
Wherein t is time.
In the collecting step, the three-axis relative attitude between the two cabins is set as
Figure BDA0002837349670000083
Three-axis relative centroid displacement of [ Delta x Delta y Delta z]Two compartments were provided with 9 displacement sensors A1, A2, A3, B1, B2, B3, C1, C2 and C3 for relative position measurement, and the measured values were [ Δ z ] respectivelyA1 ΔzA2 ΔzA3 ΔyB1 ΔyB2 ΔyB3 ΔxC1 ΔxC2 ΔxC3]The installation position of 9 displacement sensors in the mechanical coordinate system is set as [ x ]i yi zi],i=A1,…C3The center of mass of the load compartment is set in the mechanical coordinate system as xpc ypc zpc]Then, the relationship between the relative position and attitude and the measured value of the displacement sensor is as follows:
Figure BDA0002837349670000091
is provided with
Figure BDA0002837349670000092
B=[ΔzA1 ΔzA2 ΔzA3 ΔyB1 ΔyB2 ΔyB3 ΔxC1 ΔxC2 ΔxC3]T Equation 5
Figure BDA0002837349670000093
Then there is
A, X ═ B formula 7
Obtaining a pseudo-inverse solution according to a pseudo-inverse algorithm
X=AT·(A·AT)-1B equation 8
Thereby obtaining the relative attitude and the relative position of the three axes and obtaining a resolving matrix of
M=AT·(A·AT)-1Equation 9
Obtaining the relative distance P between the two interstitial centers according to the calculation modelCJThen there is
PCJ=[Δx Δy Δz]TEquation 10
In the step of fitting differential, the relative distance between the two interstitial hearts is set as PCJ_NHObtained by fitting
PCJ_NH=p1·t2+p2·t+p3Equation 11
Wherein p is1As second order fitting coefficient, p2As first-order fitting coefficient, p3Is a constant coefficient, and the linear acceleration of the relative motion between the two cabins is set as aCJThen there is
aCJ=2p1Equation 12.
In the step of resolving the quality of the platform cabin, the magnetic control force F between the cabinscbUnder the action of (2), the linear acceleration of the relative motion between the two compartments is expressed again as
Figure BDA0002837349670000101
In the above formula, aCJ、FcbAnd mpAre all known quantities, thus obtaining a platform cabin mass of
Figure BDA0002837349670000102
In the present embodiment, the parameter setting rule is as shown in fig. 2. The magnetic suspension force between two cabins is set to be 0.02N, namely
Fcb=0.02N
The system is in an open-loop uncontrolled state, under the action of magnetic control force between two cabins, the fitting result of the relative distance between the centers of the two cabins is shown in figure 4, and the linear acceleration of relative motion is obtained according to the fitting result
aCJ=2.5304×10-4m/s2
If the weight of the load cabin is 110kg, the weight of the load cabin is set to be 110kg
Figure BDA0002837349670000103
The difference between the theoretical value of 280kg and the theoretical value of 280kg is only 0.83kg, and the identification precision is better than 0.3 percent.
According to another aspect of the invention, a novel split type satellite on-orbit platform cabin quality identification system is provided, which comprises the following modules:
a driving module: powering off a repeated locking mechanism between two cabins of the split satellite, electrifying a magnetic floating actuator between the two cabins, generating magnetic control force with known magnitude, and driving the two cabins to move;
an acquisition module: the control system keeps open loop uncontrolled, collects the information measured by the displacement sensor of the two cabins, and obtains the relative movement information of the mass center of the two cabins according to the resolving model;
fitting a differential module: carrying out secondary fitting on the relative movement information of the mass center, and differentiating the fitting result twice to obtain the acceleration information of the mass center of the two cabins;
the platform cabin quality resolving module: and resolving to obtain the platform cabin mass by combining the known inter-cabin acting force, the acceleration information and the load cabin mass information.
Those skilled in the art will appreciate that, in addition to implementing the system and its various devices, modules, units provided by the present invention as pure computer readable program code, the system and its various devices, modules, units provided by the present invention can be fully implemented by logically programming method steps in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers and the like. Therefore, the system and various devices, modules and units thereof provided by the invention can be regarded as a hardware component, and the devices, modules and units included in the system for realizing various functions can also be regarded as structures in the hardware component; means, modules, units for performing the respective functions may also be regarded as structures within both software modules and hardware components for performing the methods
The foregoing description of specific embodiments of the present invention has been presented. It is to be understood that the present invention is not limited to the specific embodiments described above, and that various changes or modifications may be made by one skilled in the art within the scope of the appended claims without departing from the spirit of the invention. The embodiments and features of the embodiments of the present application may be combined with each other arbitrarily without conflict.

Claims (10)

1. A split type satellite in-orbit platform cabin quality identification method is characterized by comprising the following steps:
a driving step: powering off a repeated locking mechanism between two cabins of the split satellite, electrifying a magnetic floating actuator between the two cabins, generating magnetic control force with known magnitude, and driving the two cabins of the split satellite to move;
the collection step comprises: the control system of the split satellite is kept open-loop and uncontrolled, the measurement information of the displacement sensors of the two cabins of the split satellite is collected, and the relative motion information of the mass center of the two cabins of the split satellite is obtained according to the resolving model;
and fitting differentiation: carrying out secondary fitting on the relative motion information of the mass centers of the two cabins of the split satellite, and differentiating the fitting result twice to obtain the acceleration information of the mass centers of the two cabins of the split satellite;
a platform cabin quality resolving step: and resolving to obtain the platform cabin quality by combining the known acting force and acceleration information between the two cabins of the split satellite and the load cabin quality information.
2. The split-type satellite on-orbit platform cabin quality identification method according to claim 1, wherein in the driving step, a plurality of magnetic suspension actuators are configured between two cabins and used for completing three-axis attitude control of the load cabin and relative position control between the two cabins, an open-loop instruction is sent to enable the magnetic suspension actuators to output magnetic control force between the cabins in a combined manner, and the magnetic control force acting on the platform cabin is set as FcbThe magnetic control force acting on the load compartment is-FcbUnder the action of magnetic control force between the cabins, the acceleration of the platform cabin is abAcceleration of the load compartment of apLet the mass of the load compartment be mpMass m of the platform cabinbSo, according to Newton's second law of motion, there is
Figure FDA0002837349660000011
The relative distance between the interstitial centers of the two cabins is PCJAnd if the magnetic control force between the two cabins drives the distance between the two cabins to be larger, the distance between the two cabins is larger
Figure FDA0002837349660000012
Wherein t is time.
3. The split-type satellite in-orbit platform cabin quality identification method according to claim 2, wherein in the acquisition step, the three-axis relative attitude between two cabins is set as
Figure FDA0002837349660000013
Three-axis relative centroid displacement of [ Delta x Delta y Delta z]Between two cabins9 displacement sensors a1, a2, A3, B1, B2, B3, C1, C2, and C3 are arranged for relative position measurement, and the measured values are each set to [ Δ z [ ]A1 ΔzA2 ΔzA3 ΔyB1 ΔyB2 ΔyB3 ΔxC1 ΔxC2 ΔxC3]The installation position of 9 displacement sensors in the mechanical coordinate system is set as [ x ]i yi zi],i=A1,…C3The center of mass of the load compartment is set in the mechanical coordinate system as xpcypc zpc]Then, the relationship between the relative position and attitude and the measured value of the displacement sensor is as follows:
Figure FDA0002837349660000021
is provided with
Figure FDA0002837349660000022
B=[ΔzA1 ΔzA2 ΔzA3 ΔyB1 ΔyB2 ΔyB3 ΔxC1 ΔxC2 ΔxC3]TEquation 5
Figure FDA0002837349660000023
Then there is
A, X ═ B formula 7
Obtaining a pseudo-inverse solution according to a pseudo-inverse algorithm
X=AT·(A·AT)-1B equation 8
Thereby obtaining the relative attitude and the relative position of the three axes and obtaining a resolving matrix of
M=AT·(A·AT)-1Equation 9
Obtaining the relative distance P between the two interstitial centers according to the calculation modelCJThen, thenIs provided with
PCJ=[Δx Δy Δz]TEquation 10.
4. The split-type satellite on-orbit platform cabin quality identification method according to claim 3, wherein in the fitting differentiation step, the relative distance between the centers of two cabins is set as PCJ_NHObtained by fitting
PCJ_NH=p1·t2+p2·t+p3Equation 11
Wherein p is1As second order fitting coefficient, p2As first-order fitting coefficient, p3Is a constant coefficient, and the linear acceleration of the relative motion between the two cabins is set as aCJThen there is
aCJ=2p1Equation 12.
5. The split-type satellite in-orbit platform cabin quality identification method according to claim 4, wherein in the platform cabin quality resolving step, the magnetic control force F between the cabinscbUnder the action of (2), the linear acceleration of the relative motion between the two compartments is expressed again as
Figure FDA0002837349660000031
In the above formula, aCJ、FcbAnd mpAre all known quantities, thus obtaining a platform cabin mass of
Figure FDA0002837349660000032
6. A split type satellite in-orbit platform cabin quality identification system is characterized by comprising:
a driving module: powering off a repeated locking mechanism between two cabins of the split satellite, electrifying a magnetic floating actuator between the two cabins, generating magnetic control force with known magnitude, and driving the two cabins of the split satellite to move;
an acquisition module: the control system of the split satellite is kept open-loop and uncontrolled, the measurement information of the displacement sensors of the two cabins of the split satellite is collected, and the relative motion information of the mass center of the two cabins of the split satellite is obtained according to the resolving model;
fitting a differential module: carrying out secondary fitting on the relative motion information of the mass centers of the two cabins of the split satellite, and differentiating the fitting result twice to obtain the acceleration information of the mass centers of the two cabins of the split satellite;
the platform cabin quality resolving module: and resolving to obtain the platform cabin quality by combining the known acting force and acceleration information between the two cabins of the split satellite and the load cabin quality information.
7. The split-type satellite on-orbit platform cabin quality identification system of claim 6, wherein in the driving module, a plurality of magnetic suspension actuators are arranged between two cabins and used for completing three-axis attitude control of the load cabin and relative position control between the two cabins, an open-loop instruction is sent to enable the magnetic suspension actuators to output magnetic control force between the cabins in a combined manner, and the magnetic control force acting on the platform cabin is set as FcbThe magnetic control force acting on the load compartment is-FcbUnder the action of magnetic control force between the cabins, the acceleration of the platform cabin is abAcceleration of the load compartment of apLet the mass of the load compartment be mpMass m of the platform cabinbSo, according to Newton's second law of motion, there is
Figure FDA0002837349660000041
The relative distance between the interstitial centers of the two cabins is PCJAnd if the magnetic control force between the two cabins drives the distance between the two cabins to be larger, the distance between the two cabins is larger
Figure FDA0002837349660000042
Wherein t is time.
8. The split-type satellite in-orbit platform cabin quality identification system according to claim 7, wherein the acquisition module is configured to set a three-axis relative attitude between two cabins as
Figure FDA0002837349660000043
Three-axis relative centroid displacement of [ Delta x Delta y Delta z]Two compartments were provided with 9 displacement sensors A1, A2, A3, B1, B2, B3, C1, C2 and C3 for relative position measurement, and the measured values were [ Δ z ] respectivelyA1 ΔzA2 ΔzA3 ΔyB1 ΔyB2 ΔyB3 ΔxC1 ΔxC2 ΔxC3]The installation position of 9 displacement sensors in the mechanical coordinate system is set as [ x ]i yi zi],i=A1,…C3The center of mass of the load compartment is set in the mechanical coordinate system as xpcypc zpc]Then, the relationship between the relative position and attitude and the measured value of the displacement sensor is as follows:
Figure FDA0002837349660000044
is provided with
Figure FDA0002837349660000051
B=[ΔzA1 ΔzA2 ΔzA3 ΔyB1 ΔyB2 ΔyB3 ΔxC1 ΔxC2 ΔxC3]TEquation 5
Figure FDA0002837349660000052
Then there is
A, X ═ B formula 7
Obtaining a pseudo-inverse solution according to a pseudo-inverse algorithm
X=AT·(A·AT)-1B equation 8
Thereby obtaining the relative attitude and the relative position of the three axes and obtaining a resolving matrix of
M=AT·(A·AT)-1Equation 9
Obtaining the relative distance P between the two interstitial centers according to the calculation modelCJThen there is
PCJ=[Δx Δy Δz]TEquation 10.
9. The split-type satellite in-orbit platform cabin quality identification system according to claim 8, wherein in the fitting differential module, the relative distance between the centers of two cabins is PCJ_NHObtained by fitting
PCJ_NH=p1·t2+p2·t+p3Equation 11
Wherein p is1As second order fitting coefficient, p2As first-order fitting coefficient, p3Is a constant coefficient, and the linear acceleration of the relative motion between the two cabins is set as aCJThen there is
aCJ=2p1Equation 12.
10. The split satellite in-orbit platform cabin quality identification system according to claim 9, wherein in the platform cabin quality calculation module, the magnetic control force F between cabinscbUnder the action of (2), the linear acceleration of the relative motion between the two compartments is expressed again as
Figure FDA0002837349660000061
In the above formula, aCJ、FcbAnd mpAre all known quantities, thus obtaining a platform cabin mass of
Figure FDA0002837349660000062
CN202011475974.7A 2020-12-15 2020-12-15 Split type satellite on-orbit platform cabin quality identification method and system Active CN112660422B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011475974.7A CN112660422B (en) 2020-12-15 2020-12-15 Split type satellite on-orbit platform cabin quality identification method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011475974.7A CN112660422B (en) 2020-12-15 2020-12-15 Split type satellite on-orbit platform cabin quality identification method and system

Publications (2)

Publication Number Publication Date
CN112660422A true CN112660422A (en) 2021-04-16
CN112660422B CN112660422B (en) 2023-03-24

Family

ID=75405990

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011475974.7A Active CN112660422B (en) 2020-12-15 2020-12-15 Split type satellite on-orbit platform cabin quality identification method and system

Country Status (1)

Country Link
CN (1) CN112660422B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114408220A (en) * 2022-01-25 2022-04-29 上海卫星工程研究所 On-track calibration method and system for force arm of magnetic suspension actuator
CN114919774A (en) * 2022-05-20 2022-08-19 南京航空航天大学 On-orbit calibration method for Lorentz force actuator of non-contact load undisturbed satellite platform

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2491101A1 (en) * 2003-12-30 2005-06-30 Canadian Space Agency Zero-g emulating testbed for spacecraft control system
CN105035361A (en) * 2015-07-31 2015-11-11 上海卫星工程研究所 Satellite with ultrahigh pointing accuracy and ultrahigh stability under dynamic-static isolation and principal-subordinate cooperative control
US20170190446A1 (en) * 2015-12-31 2017-07-06 Lawrence Ellis Williams, Sr. Centripetal aerodynamic platform spacecraft
CN111099045A (en) * 2019-12-12 2020-05-05 上海卫星工程研究所 Novel double-super-satellite dynamics and control air floatation platform full-physical simulation method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2491101A1 (en) * 2003-12-30 2005-06-30 Canadian Space Agency Zero-g emulating testbed for spacecraft control system
CN105035361A (en) * 2015-07-31 2015-11-11 上海卫星工程研究所 Satellite with ultrahigh pointing accuracy and ultrahigh stability under dynamic-static isolation and principal-subordinate cooperative control
US20170190446A1 (en) * 2015-12-31 2017-07-06 Lawrence Ellis Williams, Sr. Centripetal aerodynamic platform spacecraft
CN111099045A (en) * 2019-12-12 2020-05-05 上海卫星工程研究所 Novel double-super-satellite dynamics and control air floatation platform full-physical simulation method

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114408220A (en) * 2022-01-25 2022-04-29 上海卫星工程研究所 On-track calibration method and system for force arm of magnetic suspension actuator
CN114408220B (en) * 2022-01-25 2024-01-30 上海卫星工程研究所 On-orbit calibration method and system for force arm of magnetic suspension actuator
CN114919774A (en) * 2022-05-20 2022-08-19 南京航空航天大学 On-orbit calibration method for Lorentz force actuator of non-contact load undisturbed satellite platform

Also Published As

Publication number Publication date
CN112660422B (en) 2023-03-24

Similar Documents

Publication Publication Date Title
CN112684697B (en) Split type satellite in-orbit two-cabin rotational inertia identification method and system
CN112660422B (en) Split type satellite on-orbit platform cabin quality identification method and system
Kechichian Motion in general elliptic orbit with respect to a dragging and precessing coordinate frame
CN103076807B (en) A kind of control method of drive lacking flexible spacecraft attitude stabilization
Orr et al. Space launch system ascent flight control design
Richie et al. Simultaneous attitude control and energy storage using VSCMGs: Theory and simulation
Zareb et al. Fuzzy-PID hybrid control system to navigate an autonomous mini-Quadrotor
Kusuda et al. Feedback control with nominal inputs for agile satellites using control moment gyros
Somov et al. Health checking autonomous attitude control system of Earth-observing miniature satellite in initial orientation modes
Schwartz The distributed spacecraft attitude control system simulator: from design concept to decentralized control
Ge et al. Model-free regulation of multi-link smart materials robots
Nakaya et al. Formation deployment control for spinning tethered formation flying-simulations and ground experiments
CN112577669B (en) Split type satellite in-orbit two-cabin mass center identification method and system
Nanos et al. On parameter estimation of space manipulator systems with flexible joints using the energy balance
Aydogan et al. Adaptive LQR stabilization control of reaction wheel for satellite systems
CN114115315B (en) Method and system for controlling release and recovery of tethered satellites
Maghami et al. Laser interferometer space antenna dynamics and controls model
Xuan et al. From PID to L1 adaptive control for automatic balancing of a spacecraft three-axis simulator
Reiter et al. MEMS control moment gyroscope design and wafer-based spacecraft chassis study
Culton et al. Design and development of an unrestricted satellite motion simulator
Oland et al. Path Following Control of Quadrotors Using Controllers Developed for General Euler-Lagrange Systems
CN113247310B (en) Estimation method and system suitable for continuous attitude maneuver times of satellite
Richie Combined attitude control and energy storage for small satellites using variable speed control moment gyroscopes
Mashtakov et al. Study of the accuracy provided by small satellite attitude determination & control system
Richie et al. A practical variable-speed control moment gyroscope steering law for small satellite energy storage and attitude control

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant