Joint position keeping control method and device in case of centroid deviation

Abstract

The present application relates to the technical field of spacecraft attitude orbit control, and particularly relates to a joint position keeping control method and device under the condition of mass center deviation. The method comprises the following steps: determining the direction of satellite drift caused by the tangential thrust of coplanarly installed electric thrusters when the mass center of the whole satellite deviates by 1mm, and the semi-major axis variation of the satellite in the position keeping control period, determining the target deviation direction and target deviation amount of the mass center of the whole satellite based on the drift direction of the satellite prime meridian caused by the perturbation acceleration and the semi-major axis variation of the satellite in the position keeping control period under the action of the earth perturbation force; adjusting the mass center position of the satellite based on the target deviation direction and target deviation amount; after the satellite is fixed in orbit, the coplanarly installed electric thrusters are used to autonomously perform north-south position keeping control every operation day, and since the deviation of the mass center of the whole satellite can simultaneously generate tangential thrust, the north-south position keeping control is completed, and the east-west position keeping control task is also completed.

Description

Technical Field

[0001] This invention relates to the field of spacecraft attitude and orbit control technology, and in particular to a combined position holding control method and apparatus under the condition of center of mass deviation. Background Technology

[0002] Electric propulsion, a type of low-thrust propulsion system, is characterized by its light weight, high specific impulse, and low thrust. After the satellite is positioned, coplanar electric thrusters can be used to maintain its north-south position throughout its entire lifespan. Then, chemical thrusters can be used to maintain its east-west position, thereby achieving orbital control of the satellite.

[0003] This method of decoupling the north-south position-keeping mission from the east-west position-keeping mission requires a large amount of chemical propellant for the chemical thruster. To reduce the consumption of chemical propellant, the north-south position-keeping mission and the east-west position-keeping mission can be coupled together, using coplanar electric thrusters to complete the north-south position-keeping mission while simultaneously completing the east-west position-keeping mission.

[0004] Therefore, there is an urgent need for a new joint position holding control method based on coplanarly mounted electric thrusters. Summary of the Invention

[0005] To address the problem that existing decoupling methods are not only time-consuming but also require a large amount of chemical propellant, embodiments of the present invention provide a joint position holding control method and apparatus under centroid deviation conditions.

[0006] In a first aspect, embodiments of the present invention provide a joint position-maintaining control method under centroid deviation conditions, comprising:

[0007] The perturbation acceleration caused by the satellite's fixed position determines the direction of the satellite's mean longitude drift and the change in the satellite's semi-major axis during the position-maintaining control period under the influence of Earth's perturbation force; wherein, the position-maintaining control period is 7 days and the operation day is 5 days;

[0008] Determine the direction of satellite drift caused by the tangential thrust generated by the coplanar electric thrusters when the center of mass of the entire satellite deviates by 1 mm in the +X or -X direction, as well as the amount of change of the semi-major axis of the satellite during the position-holding control period. Based on the direction of satellite longitude drift caused by the perturbation acceleration and the amount of change of the semi-major axis of the satellite during the position-holding control period under the action of the Earth's perturbation force, determine the target offset direction and target offset of the center of mass of the entire satellite.

[0009] Based on the target offset direction and target offset amount, adjust the position of the satellite's center of mass.

[0010] After the satellite is positioned in orbit, the ignition duration and right ascension of the ignition point are determined for each operational day. The coplanar electric thrusters are then used to autonomously perform north-south position-keeping control on each operational day to complete the joint position-keeping control task.

[0011] Secondly, embodiments of the present invention also provide a combined position-holding control device based on the method described in any embodiment of this specification, comprising:

[0012] The first determining unit is used to determine the direction of the satellite's mean longitude drift caused by the perturbation acceleration and the change in the satellite's semi-major axis during the position-holding control period based on the satellite's fixed position; wherein, the position-holding control period is 7 days and the operation day is 5 days;

[0013] The second determining unit is used to determine the direction of satellite drift caused by the tangential thrust generated by the coplanar electric thrusters when the center of mass of the entire satellite deviates by 1 mm in the +X or -X direction, as well as the amount of change of the semi-major axis of the satellite during the position holding control period. Based on the direction of satellite longitude drift caused by the perturbation acceleration and the amount of change of the semi-major axis of the satellite during the position holding control period under the action of the Earth's perturbation force, the target offset direction and target offset of the center of mass of the entire satellite are determined.

[0014] An adjustment unit is used to adjust the position of the satellite's center of mass based on the target offset direction and the target offset amount.

[0015] The control unit is used to determine the ignition duration and right ascension of the ignition point for each operational day after the satellite has reached its designated orbit, and to autonomously perform north-south position-keeping control on each operational day using the coplanar-mounted electric thrusters to complete the joint position-keeping control task.

[0016] Thirdly, embodiments of the present invention also provide a computing device, including a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, it implements the method described in any embodiment of this specification.

[0017] Fourthly, embodiments of the present invention also provide a computer-readable storage medium having a computer program stored thereon, which, when executed in a computer, causes the computer to perform the methods described in any embodiment of this specification.

[0018] This invention provides a combined position-keeping control method and apparatus under conditions of center-of-mass deviation. By adjusting the position of the satellite's center of mass based on the target offset direction and the target offset amount, when the satellite is performing north-south position-keeping control in orbit, the offset of the satellite's center of mass will generate tangential thrust, thereby raising or lowering the semi-major axis. This can precisely offset the change in the satellite's semi-major axis during the position-keeping control period under the influence of Earth's perturbation. Therefore, while completing north-south position-keeping control, east-west position-keeping control is also completed. This not only eliminates the need for extra time spent on east-west position-keeping but also reduces the consumption of chemical propellant by using the tangential thrust generated during north-south position-keeping control to adjust the east-west position. Attached Figure Description

[0019] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0020] Figure 1 This is a flowchart of a joint position holding control method under centroid deviation conditions provided by an embodiment of the present invention;

[0021] Figure 2 This is a schematic diagram illustrating the north-south position maintenance according to an embodiment of the present invention;

[0022] Figure 3 This is a hardware architecture diagram of a computing device provided in an embodiment of the present invention;

[0023] Figure 4 This is a structural diagram of a combined position holding control device under centroid deviation conditions provided in an embodiment of the present invention. Detailed Implementation

[0024] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are some embodiments of the present invention, but not all embodiments. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without creative effort are within the scope of protection of the present invention.

[0025] The following describes the specific implementation of the above concept.

[0026] Please refer to Figure 1 This invention provides a joint position holding control method under centroid deviation conditions, the method comprising:

[0027] Step 100: Based on the satellite's fixed position, determine the direction of the satellite's mean longitude drift caused by the perturbation acceleration and the change in the satellite's semi-major axis during the position-holding control period under the influence of Earth's perturbation force; wherein, the position-holding control period is 7 days and the operation day is 5 days;

[0028] Step 102: Determine the direction of satellite drift caused by the tangential thrust generated by the coplanar electric thrusters when the center of mass of the entire satellite deviates by 1 mm in the +X or -X direction, as well as the amount of change of the semi-major axis of the satellite during the position-holding control period. Based on the direction of satellite drift in longitude caused by perturbation acceleration and the amount of change of the semi-major axis of the satellite during the position-holding control period under the action of Earth perturbation, determine the target offset direction and target offset of the center of mass of the entire satellite.

[0029] Step 104: Adjust the position of the satellite's center of mass based on the target offset direction and target offset amount;

[0030] Step 106: After the satellite is in orbit and positioned, determine the ignition duration and right ascension of the ignition point for each operational day, and use coplanar electric thrusters to autonomously perform north-south position-keeping control on each operational day to complete the joint position-keeping control task.

[0031] In this embodiment of the invention, the position of the satellite's center of mass is adjusted based on the target offset direction and the target offset amount. When the satellite performs north-south position-keeping control in orbit, the offset of the satellite's center of mass causes the electric thrusters to generate tangential thrust in the required direction, thereby raising or lowering the semi-major axis. This can precisely offset the change in the satellite's semi-major axis during the position-keeping control period under the influence of Earth's perturbation. Therefore, the east-west position-keeping control task is completed at the same time as the north-south position-keeping control. Not only is there no need to waste extra time on east-west position-keeping, but the east-west position adjustment using the tangential thrust generated during north-south position-keeping control can also reduce the consumption of chemical propellant.

[0032] For step 100:

[0033] In this embodiment of the invention, in order to determine the target offset direction and target offset amount of the whole satellite's center of mass shift, it is necessary to determine the direction of the satellite's mean longitude drift caused by the perturbation acceleration and the amount of change of the satellite's semi-major axis during the position-maintaining control period under the action of the Earth's perturbation force, based on the satellite's fixed position.

[0034] For example, when the satellite's fixed position is at 134 degrees east longitude, the local drift acceleration is -0.00163 degrees / day^2. Since the drift acceleration is negative, the perturbation acceleration causes the satellite's mean longitude to drift westward. Under the influence of the Earth's perturbation force, the satellite's semi-major axis rises by about 127m during the position-maintaining control period.

[0035] Regarding steps 102 and 104:

[0036] You can refer to this. Figure 2 This scheme utilizes a two-degree-of-freedom vector adjustment mechanism to achieve coplanar mounting of the electric thrusters. The electric thrusters can provide thrust in both the radial and normal directions, enabling control of orbital inclination and eccentricity. The radial direction points towards the Earth's center, while the normal direction is perpendicular to the orbital plane. When the satellite's center of mass shifts, the electric thrusters generate tangential thrust during inclination control. When the satellite shifts westward and its semi-major axis rises due to perturbations, the required tangential thrust in the desired direction and magnitude can adjust the satellite's east-west position, ensuring it follows its predetermined orbit.

[0037] In some implementations, the change in the semi-major axis of the satellite during the position-holding control period caused by the tangential thrust generated by the coplanarly mounted electric thrusters for every 1 mm deviation of the satellite's center of mass is determined by the following steps B1-B4:

[0038] Step B1: Obtain the tangential and normal thrust of the primary and backup electric thrusters, which are coplanarly mounted, when the center of mass of the entire satellite deviates by 1 mm.

[0039] In this embodiment, the coplanarly mounted electric thrusters include a primary electric thruster and a backup electric thruster, and the thrust of the primary electric thruster and the backup electric thruster are different. It is necessary to obtain the tangential thrust and normal thrust of the primary electric thruster and the tangential thrust and normal thrust of the backup electric thruster respectively.

[0040] Step B2: Based on the tangential thrust, normal thrust, and the pre-set normal daily average velocity increment, calculate the corresponding tangential daily average velocity increment when the primary electric thruster is started and when the backup electric thruster is started.

[0041] In some implementations, the tangential daily average velocity increment is calculated using the following formula:

[0042]

[0043] In the formula, ΔV t ΔV represents the daily average velocity increment in the tangential direction. n For the daily average velocity increment in the normal direction, F t For tangential thrust, F n This is the normal thrust.

[0044] In this embodiment, the daily average normal velocity increment is determined as follows: the tilt drift velocity is approximately 0.75° to 0.95° per year, and the velocity increment Δv maintained at the north-south position is between 41 m / s and 51 m / s. The average value Δv = 46 m / s is used in the design. Dividing Δv by 365 days yields the daily average normal velocity increment ΔV. n It is 0.126 m / s.

[0045] In this embodiment, since the thrust of the primary electric thruster and the backup electric thruster are different, it is necessary to use the above formula to calculate the tangential daily average velocity increment when using the primary electric thruster alone and the tangential daily average velocity increment when using the backup electric thruster alone.

[0046] Step B3: Based on the vitality formula and the daily average tangential velocity increment, calculate the daily average semi-major axis change when the primary electric thruster is started and when the backup electric thruster is started.

[0047] In some implementations, the daily average change in the semi-major axis is calculated using the following formula:

[0048]

[0049] In the formula, ΔV t denoted as the daily average velocity increment, v as the satellite velocity, μ as the gravitational constant, r as the satellite radius vector, a as the semi-major axis, and Δa as the daily average change in the semi-major axis.

[0050] From the above formula, we can obtain the relationship Δa=27.4×ΔV t .

[0051] Similarly, by substituting the daily average tangential velocity increment corresponding to the use of the primary electric thruster alone into the formula, we can obtain the daily average semi-major axis change when the primary electric thruster is started; by substituting the daily average tangential velocity increment corresponding to the use of the backup electric thruster alone into the formula, we can obtain the daily average semi-major axis change when the backup electric thruster is started.

[0052] Step B4: Based on the daily average semi-major axis change and the number of days in the position holding control period, determine the amount of semi-major axis change caused by the tangential thrust generated by the coplanar electric thrusters when the center of mass of the satellite deviates by 1 mm during the position holding control period.

[0053] In this embodiment, multiplying the daily average semi-major axis change when the primary electric thruster is activated by 7 yields the semi-major axis change caused by the tangential thrust generated by the primary electric thruster during the position-holding control period when the satellite's center of mass deviates by 1 mm. Multiplying the daily average semi-major axis change when the backup electric thruster is activated by 7 yields the semi-major axis change caused by the tangential thrust generated by the backup electric thruster during the position-holding control period when the satellite's center of mass deviates by 1 mm.

[0054] Therefore, the effect of each 1mm deviation of the entire star's center of mass on the semi-major axis is shown in Table 1 below.

[0055] Table 1

[0056]

[0057] In some implementations, the target offset direction and target offset amount of the entire star's center of mass are determined in the following manner:

[0058] When the center of mass of the satellite is offset in the +X and -X directions respectively, the direction of satellite drift caused by the tangential thrust generated by the coplanar electric thrusters is determined. The direction of the offset of the center of mass of the satellite corresponding to the direction opposite to the direction of satellite drift in longitude caused by perturbation acceleration is determined as the target offset direction of the center of mass of the satellite.

[0059] The change in the semi-major axis of the satellite during the position-holding control period under the influence of Earth's perturbation is divided by the change in the semi-major axis of the satellite during the position-holding control period caused by the tangential thrust generated by the coplanar electric thrusters when the satellite's center of mass deviates from the target offset direction by 1 mm, thus yielding the target offset of the satellite's center of mass.

[0060] Continuing with the example of a fixed position at 134 degrees east longitude, step 100 has determined that when the satellite's fixed position is at 134 degrees east longitude, the perturbation acceleration causes the satellite's mean longitude to drift westward. Under the influence of the Earth's perturbation force, the satellite's semi-major axis rises by approximately 127 meters during the position maintenance control period.

[0061] When the center of mass of the entire satellite shifts in the +X direction, the tangential thrust generated by the coplanar electric thrusters causes the satellite to accelerate westward. When the center of mass of the entire satellite shifts in the -X direction, the tangential thrust generated by the coplanar electric thrusters causes the satellite to accelerate eastward. The perturbation acceleration causes the satellite's mean longitude to drift westward. Therefore, the target shift direction of the center of mass of the entire satellite is towards the -X direction.

[0062] As can be seen from Table 1 above, the primary electric thruster is usually used. As can be seen from the table above, the change in the semi-major axis of the satellite during the position-keeping control period of 127m under the action of Earth perturbation is divided by the change in the semi-major axis of 11.55m during the position-keeping control period corresponding to the primary electric thruster in Table 1, and the target offset of the entire satellite's center of mass is obtained as 11mm.

[0063] Therefore, when the satellite is on Earth, adjusting the satellite's center of mass to shift 11mm in the -X direction will cause the electric propulsion ignition during the position holding control cycle to accelerate the satellite eastward drift. The electric propulsion effect is opposite to the direction of Earth's natural perturbation, producing a canceling effect, which helps to reduce the frequency of east-west position holding and reduce propellant consumption.

[0064] Regarding step 106:

[0065] In some implementations, step 106 may include the following steps 200-206:

[0066] Step 200: When an autonomous position holding control command is received, determine the change in orbital inclination vector within the position holding control period based on the current satellite time and the position holding control period;

[0067] Step 202: Based on the change in the orbital tilt vector within the position holding control cycle, determine the total ignition duration within the position holding control cycle, so as to determine the ignition duration corresponding to each operating day within the position holding control cycle;

[0068] Step 204: Determine the right ascension of the ignition point based on the change in the orbital inclination vector during the position holding control cycle;

[0069] Step 206: Based on the ignition duration and right ascension of the ignition point corresponding to each operation day, the coplanar electric thrusters are used to autonomously perform north-south position holding control within each operation day to complete the joint position holding control task.

[0070] In this embodiment of the invention, the satellite autonomously calculates the ignition start time and duration for multiple operational days based on the current satellite time, a pre-determined position-holding control period, and current onboard orbital parameters, enabling the satellite to autonomously complete the multi-day position-holding process. This solution not only predicts the ignition duration and start time for multiple operational days but also simplifies the calculation, avoids orbital recursion, and is more conducive to onboard implementation.

[0071] For step 200:

[0072] In some implementations, step 200, "determining the change in orbital inclination vector within the position-holding control period based on the current satellite time and position-holding control period," may include steps S1-S4:

[0073] S1, based on the current star time, the initial star time and J2000 time, calculate the Julian century, and determine the ascending node longitude based on the Julian century;

[0074] S2, based on the ecliptic longitude of the ascending node, calculates the rate of change of the orbital inclination vector over one year;

[0075] S3, transform the zero tilt vector of the instantaneous true coordinate system to the J2000 coordinate system, and obtain the magnitude of the orbit tilt vector in the x and y directions corresponding to the zero tilt vector in the instantaneous true coordinate system in the J2000 coordinate system;

[0076] S4. Based on the rate of change of the track inclination vector within one year, the magnitude of the track inclination vector in the x and y directions corresponding to the zero inclination vector in the instantaneous true coordinate system in the J2000 coordinate system, and the position holding control period, determine the change of the track inclination vector within the position holding control period.

[0077] In some implementations, the rate of change of the orbital inclination vector over a year is determined as follows, i.e., steps S1-S2 may include:

[0078] The Julian century can be calculated using the following formula:

[0079]

[0080] in,

[0081] t J2000 =t+473342400

[0082] In the formula, T represents the Julian century, and t J2000 t represents the number of seconds between the current satellite time and the J2000 time, and t represents the number of seconds between the current satellite time and the initial satellite time.

[0083] It should be noted that J2000.0 is an epoch used in astronomy. The prefix "J" indicates that it is a Julian epoch, not a Bessel epoch. It refers to the Julian date TT2451545.0, or TT time January 1, 2000, 12:00.

[0084] Calculate the ecliptic longitude of the ascending node using the following formula:

[0085] Ω l =2.182439 - 33.757 × T

[0086] In the formula, T represents the Julian century, and Ω l The ecliptic longitude of the ascending node;

[0087] The rate of change of the orbital inclination vector over a year can be calculated using the following formula:

[0088] dΔi x =0.097cosΩ l +0.852

[0089] dΔi y =0.13sinΩ l

[0090] In the formula, dΔi x ′ and dΔi y ′ represent the rates of change of the orbital inclination vector in the x-direction and y-direction, respectively, over a year, Ω l The ecliptic longitude of the ascending node;

[0091] The zero tilt vector of the instantaneous true coordinate system is transformed to the J2000 coordinate system using the following formula, i.e., step S3 includes:

[0092] Δi xTOD =KTodcof1*C PR (2,0)

[0093] Δi yTOD =KTodcof1*C PR (2,1)

[0094] In the formula, Δi xTOD and Δi yTOD These represent the x- and y-direction magnitudes of the orbital inclination vector corresponding to the zero inclination vector in the instantaneous true coordinate system in the J2000 coordinate system, respectively. KTodcof1 is a pre-set coordinate system transformation adjustment coefficient, and C... PR (2,0) is the element in the 3rd row and 1st column of the precession nutation transformation matrix, where 2 and 0 are the row and column numbers, respectively. PR (2,1) is the element in the 3rd row and 2nd column of the precession nutation transformation matrix, where 2 and 1 are the row number and column number, respectively.

[0095] In some implementations, the change in the track inclination vector during the position holding control cycle is calculated using the following formula, i.e., step S4 may include:

[0096] Δi x =dΔi x ′ / 365*TJetPeriod*Kcomp+(i·sin(Ω)-Δi xTOD R2D

[0097] Δi y =dΔi y ′ / 365*TJetPeriod*Kcomp+(-i·cos(Ω)-Δi yTOD R2D

[0098] In the formula, Δi x and Δi y dΔi represents the change in the orbital inclination vector in the x-direction and the change in the y-direction, respectively, during the position holding control cycle. x ′ and dΔi y ′ represents the rate of change of the orbital inclination vector in the x and y directions respectively over one year; TJetPeriod is the predetermined position-keeping control period, which is 7; Kcomp is the correction value, which defaults to 1; i is the satellite's orbital inclination; Ω is the right ascension of the satellite's ascending node; Δi xTOD and Δi yTODR2D represents the magnitudes of the orbital inclination vector in the x and y directions corresponding to the zero inclination vector in the instantaneous true coordinate system in the J2000 coordinate system, respectively, and R2D is the conversion between radians and degrees.

[0099] Regarding step 202:

[0100] In some implementations, step 202 may include steps A1-A3:

[0101] A1, based on the change in track inclination vector during the position holding control cycle, calculate the magnitude of track inclination change, and determine the required velocity increment in the normal direction based on the magnitude of track inclination change.

[0102] In this embodiment of the invention, the magnitude of the change in track inclination Δi is calculated using the following formula. jet The required velocity increment Δv in the normal direction jet :

[0103]

[0104]

[0105] In the formula, Δi x and Δi y These represent the changes in the track inclination vector in the x-direction and the changes in the y-direction, respectively, during the position holding control cycle.

[0106] A2 determines the total ignition duration within the position-holding control cycle based on the required velocity increment in the normal direction, the satellite's mass, and the satellite's nominal thrust.

[0107] In this embodiment of the invention, the total ignition duration ΔT within the position holding control cycle is calculated using the following formula. jet :

[0108]

[0109] In the formula, Msa is the satellite mass, and F jetnorm The nominal thrust is calculated using the south main component as the nominal value by default, corresponding to 0.0297 N.

[0110] A3. Divide the total ignition duration within the position holding control cycle into 5 equal parts to obtain the ignition duration corresponding to each operating day.

[0111] In this embodiment, position holding control is typically required five days a week, so the number of operating days is 5. The total ignition duration ΔT is calculated as follows. jet Dividing the data into 5 equal parts, we can obtain the ignition duration ΔT for each operating day. id .

[0112] In some implementations, after dividing the total ignition duration within the position holding control cycle into five equal parts in A3, before obtaining the ignition duration corresponding to each operating day, the method further includes:

[0113] Obtain the time limit for coplanar mounted electric thrusters;

[0114] Determine whether the total ignition time, after being divided equally, exceeds the time limit.

[0115] If the limit is exceeded, the time limit will be set as the ignition time corresponding to each operation day;

[0116] If the total ignition time is not exceeded, the total ignition time will be divided equally and used as the ignition time for each operation day.

[0117] In this embodiment, to ensure the safe operation of the coplanar-mounted electric thruster, a daily working time limit is generally set for the coplanar-mounted electric thruster. Under normal circumstances, the total ignition time after being divided equally is less than the limit. However, in order to further protect the coplanar-mounted electric thruster, it is necessary to determine whether the total ignition time after being divided equally exceeds the limit. If it exceeds, the limit is determined as the ignition time corresponding to each operating day. If it does not exceed, the total ignition time after being divided equally is used as the ignition time corresponding to each operating day.

[0118] Regarding step 204:

[0119] In some implementations, the ignition point ascension is calculated using the following formula:

[0120] OmgF=atan2(Δi x ,Δi y )

[0121] In the formula, OmgF is the right ascension of the ignition point, and Δi x and Δi y These represent the changes in the track inclination vector in the x-direction and the changes in the y-direction, respectively, during the position holding control cycle.

[0122] Regarding step 206:

[0123] In this step, the ignition duration ΔT corresponding to each operating day in step 202 is obtained. id Following the ignition point right ascension determined in step 204, ignition is performed once at the ascending node and once at the descending node on each operating day. The sum of the two ignition times must equal the corresponding ignition time ΔT for each operating day. id Therefore, on each operational day, when the satellite reaches the right ascension of the two ignition points, it activates the coplanar-mounted electric thrusters for position holding control, until the total duration reaches the ignition duration ΔT. idThe satellite stops automatically during north-south position holding. Since the center of mass of the entire satellite has been pre-adjusted according to the target offset direction and amount, the electric thrusters generate tangential thrust in the required direction during north-south position holding, thereby raising or lowering the semi-major axis. This effectively counteracts the change in the satellite's semi-major axis during the position holding control period caused by Earth's perturbation. Therefore, east-west position holding is completed simultaneously with north-south position holding, eliminating the need for additional time spent on east-west position holding. Furthermore, adjusting the east-west position using the tangential thrust generated during north-south position holding reduces chemical propellant consumption.

[0124] like Figure 3 , Figure 4 As shown, embodiments of the present invention provide a combined position-holding control device based on the method described in any embodiment of this specification. The device embodiment can be implemented in software, hardware, or a combination of both. From a hardware perspective, as... Figure 3 The diagram shown is a hardware architecture diagram of a computing device containing a joint position-keeping control device under centroid deviation conditions, as provided in an embodiment of the present invention. (Except for...) Figure 3 In addition to the processor, memory, network interface, and non-volatile memory shown, the computing device in the embodiment may also include other hardware, such as a forwarding chip responsible for processing packets. Taking software implementation as an example, such as... Figure 4 As shown, a device in a logical sense is formed by the CPU of its computing device reading the corresponding computer program from non-volatile memory into memory and running it. This embodiment provides a joint position-keeping control device for centroid deviation based on any embodiment of this specification, comprising:

[0125] The first determining unit 401 is used to determine the direction of the satellite's mean longitude drift caused by the perturbation acceleration and the change in the satellite's semi-major axis during the position-holding control period based on the satellite's fixed position; wherein, the position-holding control period is 7 days and the operation day is 5 days;

[0126] The second determining unit 402 is used to determine the direction of satellite drift caused by the tangential thrust generated by the coplanar electric thrusters when the center of mass of the satellite deviates by 1 mm in the +X or -X direction, as well as the amount of change of the semi-major axis of the satellite during the position holding control period. Based on the direction of satellite drift in longitude caused by perturbation acceleration and the amount of change of the semi-major axis of the satellite during the position holding control period under the action of Earth perturbation, the target offset direction and target offset of the center of mass of the satellite are determined.

[0127] Adjustment unit 403 is used to adjust the position of the satellite's center of mass based on the target offset direction and the target offset amount;

[0128] The control unit 404 is used to determine the ignition duration and right ascension of the ignition point for each operating day after the satellite is in orbit and to autonomously perform north-south position holding control using coplanar electric thrusters on each operating day in order to complete the joint position holding control task.

[0129] It is understood that the structures illustrated in the embodiments of the present invention do not constitute a specific limitation on a combined position holding control device under centroid deviation conditions. In other embodiments of the present invention, a combined position holding control device under centroid deviation conditions may include more or fewer component units than illustrated, or combine certain component units, or split certain component units, or arrange different component units. The illustrated components may be implemented in hardware, software, or a combination of software and hardware.

[0130] The information interaction and execution process between the modules in the above-mentioned device are based on the same concept as the method embodiment of the present invention, and the specific details can be found in the description of the method embodiment of the present invention, and will not be repeated here.

[0131] This invention also provides a computing device, including a memory and a processor. The memory stores a computer program, and when the processor executes the computer program, it implements a joint position holding control method under centroid deviation conditions according to any embodiment of this invention.

[0132] This invention also provides a computer-readable storage medium storing a computer program, which, when executed by a processor, causes the processor to perform a joint position-maintaining control method under centroid deviation conditions according to any embodiment of this invention.

[0133] Specifically, a system or apparatus equipped with a storage medium may be provided, on which software program code implementing the functions of any of the embodiments described above is stored, and the computer (or CPU or MPU) of the system or apparatus may read and execute the program code stored in the storage medium.

[0134] In this case, the program code read from the storage medium can itself implement the function of any of the above embodiments, and therefore the program code and the storage medium storing the program code constitute part of the present invention.

[0135] Examples of storage media used to provide program code include floppy disks, hard disks, magneto-optical disks, optical disks (such as CD-ROM, CD-R, CD-RW, DVD-ROM, DVD-RAM, DVD-RW, DVD+RW), magnetic tapes, non-volatile memory cards, and ROMs. Alternatively, program code can be downloaded from a server computer via a communication network.

[0136] Furthermore, it should be clear that not only can the program code read by the computer be executed, but also the operating system or other components operating on the computer can be instructed based on the program code to perform some or all of the actual operations, thereby realizing the function of any of the embodiments described above.

[0137] Furthermore, it is understood that the program code read from the storage medium is written to the memory set in the expansion board inserted into the computer or to the memory set in the expansion module connected to the computer. Then, based on the instructions of the program code, the CPU or other components installed on the expansion board or expansion module execute some and all of the actual operations, thereby realizing the function of any of the above embodiments.

[0138] It should be noted that, in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus.

[0139] Those skilled in the art will understand that all or part of the steps of the above method embodiments can be implemented by hardware related to program instructions. The aforementioned program can be stored in a computer-readable storage medium. When the program is executed, it performs the steps of the above method embodiments. The aforementioned storage medium includes various media that can store program code, such as ROM, RAM, magnetic disk, or optical disk.

[0140] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A joint position holding control method under centroid deviation conditions, characterized in that, include: The perturbation acceleration caused by the satellite's fixed position determines the direction of the satellite's mean longitude drift and the change in the satellite's semi-major axis during the position-maintaining control period under the influence of Earth's perturbation force; wherein, the position-maintaining control period is 7 days and the operation day is 5 days; Determine the direction of satellite drift caused by the tangential thrust generated by the coplanar electric thrusters when the center of mass of the entire satellite deviates by 1 mm in the +X or -X direction, as well as the amount of change of the semi-major axis of the satellite during the position-holding control period. Based on the direction of satellite longitude drift caused by the perturbation acceleration and the amount of change of the semi-major axis of the satellite during the position-holding control period under the action of the Earth's perturbation force, determine the target offset direction and target offset of the center of mass of the entire satellite. Based on the target offset direction and target offset amount, adjust the position of the satellite's center of mass. After the satellite is positioned in orbit, the ignition duration and right ascension of the ignition point are determined for each operational day. The coplanar electric thrusters are used to autonomously perform north-south position-keeping control on each operational day to complete the joint position-keeping control task. The change in the semi-major axis of the satellite during the position-holding control period caused by the tangential thrust generated by the coplanar electric thrusters when the center of mass deviates by 1 mm is determined as follows: The tangential and normal thrusts of the primary and backup electric thrusters, which are coplanarly mounted, are obtained when the center of mass of the entire satellite deviates by 1 mm. Based on the tangential thrust, the normal thrust, and the preset normal daily average velocity increment, calculate the corresponding tangential daily average velocity increment when the primary electric thruster is started and when the backup electric thruster is started, respectively. Based on the vitality formula and the tangential daily average velocity increment, the daily average semi-major axis change corresponding to the start of the primary electric thruster and the start of the backup electric thruster are calculated respectively. Based on the daily average semi-major axis change and the number of days in the position holding control period, the tangential thrust generated by the coplanar electric thrusters caused by the 1mm deviation of the satellite's center of mass results in the change of the satellite's semi-major axis within the position holding control period. The daily average tangential velocity increment is calculated using the following formula: In the formula, For the tangential daily average velocity increment, This represents the daily average velocity increment in the normal direction. For tangential thrust, Normal thrust; The daily average semi-major axis change is calculated using the following formula: In the formula, For the tangential daily average velocity increment, For satellite speed, The gravitational constant, For satellite radius vector, For the semi-major axis, This represents the daily average change in the semi-major axis; The target offset direction and target offset amount of the entire star's center of mass are determined in the following manner: When the center of mass of the entire satellite is offset in the +X and -X directions respectively, the direction of satellite drift caused by the tangential thrust generated by the coplanar electric thrusters is determined. The direction of the offset of the center of mass of the entire satellite, which is opposite to the direction of the satellite's mean longitude drift caused by the perturbation acceleration, is determined as the target offset direction of the center of mass of the entire satellite. The target offset of the satellite's center of mass is obtained by dividing the change in the satellite's semi-major axis during the position-holding control period under the influence of Earth's perturbation by the tangential thrust generated by the coplanar electric thrusters when the satellite's center of mass deviates from the target offset direction by 1 mm.

2. The method according to claim 1, characterized in that, After the satellite reaches its designated orbit, the ignition duration and right ascension of the ignition point are determined for each operational day. The coplanar-mounted electric thrusters then autonomously perform north-south position-keeping control on each operational day to complete the joint position-keeping control task, including: When an autonomous position-holding control command is received, the change in orbital inclination vector within the position-holding control period is determined based on the current satellite time and the position-holding control period. Based on the change in the orbital inclination vector within the position holding control cycle, the total ignition duration within the position holding control cycle is determined, so as to determine the ignition duration corresponding to each operating day within the position holding control cycle. The right ascension of the ignition point is determined based on the change in the orbital inclination vector during the position holding control cycle. Based on the ignition duration and right ascension of the ignition point for each operating day, the coplanar electric thrusters autonomously perform north-south position holding control within each operating day to complete the joint position holding control task.

3. The method according to claim 2, characterized in that, The determination of the orbital inclination vector change within the position-holding control period based on the current satellite time and the position-holding control period includes: Based on the current star time, the initial star time, and the J2000 time, calculate the Julian century, and determine the ascending node ecliptic longitude based on the Julian century; Based on the ecliptic longitude of the ascending node, calculate the rate of change of the orbital inclination vector over one year; Transform the zero tilt vector in the instantaneous true coordinate system to the J2000 coordinate system to obtain the magnitude of the orbit tilt vector in the x and y directions corresponding to the zero tilt vector in the instantaneous true coordinate system in the J2000 coordinate system. Based on the rate of change of the track inclination vector within one year, the magnitude of the track inclination vector in the x and y directions corresponding to the zero inclination vector in the instantaneous true coordinate system in the J2000 coordinate system, and the position holding control period, the change in the track inclination vector within the position holding control period is determined.

4. A combined position holding control device based on the method of any one of claims 1-3, characterized in that, include: The first determining unit is used to determine the direction of the satellite's mean longitude drift caused by the perturbation acceleration and the change in the satellite's semi-major axis during the position-holding control period based on the satellite's fixed position; wherein, the position-holding control period is 7 days and the operation day is 5 days; The second determining unit is used to determine the direction of satellite drift caused by the tangential thrust generated by the coplanar electric thrusters when the center of mass of the entire satellite deviates by 1 mm in the +X or -X direction, as well as the amount of change of the semi-major axis of the satellite during the position holding control period. Based on the direction of satellite longitude drift caused by the perturbation acceleration and the amount of change of the semi-major axis of the satellite during the position holding control period under the action of the Earth's perturbation force, the target offset direction and target offset of the center of mass of the entire satellite are determined. An adjustment unit is used to adjust the position of the satellite's center of mass based on the target offset direction and the target offset amount. The control unit is used to determine the ignition duration and right ascension of the ignition point for each operational day after the satellite has reached its designated orbit, and to autonomously perform north-south position-keeping control on each operational day using the coplanar-mounted electric thrusters to complete the joint position-keeping control task.

5. A computing device comprising a memory and a processor, wherein the memory stores a computer program, and the processor, when executing the computer program, implements the method as described in any one of claims 1-3.

6. A computer-readable storage medium having a computer program stored thereon, which, when executed in a computer, causes the computer to perform the method of any one of claims 1-3.

Images