COPLANAR MULTI-TURN RENDEZVOUS IN NEAR-CIRCULAR ORBIT USING A LOW-THRUST ENGINE Текст научной статьи по специальности «Физика»

Вестник РУДН. Серия: Инженерные исследования RUDN Journal of Engineering Research

2022;23(4):283-292

ISSN 2312-8143 (Print); ISSN 2312-8151 (Online) journals.rudn.ru/engineering-researches

DOI 10.22363/2312-8143-2022-23-4-283-292 UDC 629.78

Research article / Научная статья

Coplanar multi-turn rendezvous in near-circular orbit using a low-thrust engine

Andrey A. Baranov" , Adilson P. Olivio'

aKeldysh Institute of Applied Mathematics, Russian Academy of Sciences, Moscow, Russian Federation bPeoples' Friendship University of Russia (RUDN University), Moscow, Russian Federation

^ [email protected]

Article history

Received: September 27, 2022 Revised: November 25, 2022 Accepted: November 28, 2022

Keywords:

spacecraft, calculation of maneuver parameters, space object, low-thrust engine, circular orbit, velocity impulse

Abstract. The authors describe an algorithm that allows calculating the parameters of maneuvers performed on several turns by the low-thrust engine, which ensure the flight of the active spacecraft to the specified vicinity of the target space object. The movement takes place in the vicinity of a circular orbit. Linearized equations of motion are used in solving the problem. The influence of the non-centrality of the gravitational field and the atmosphere is not taken into account. The determination of maneuver parameters takes place in three stages. At the first and third stages, the parameters of the pulse transition and the transition performed by the low-thrust engine, are determined analytically. At the second stage, the distribution of maneuvering between turns, which provides a solution to the meeting problem, is carried out by iterating over one variable. This method of solving the problem provides simplicity and high reliability of determining the parameters of maneuvers, which allows it to be used on board the spacecraft. The paper investigates the dependence of the total characteristic speed of solving the meeting problem on the number of turns of the flight and the magnitude of the engine thrust.

For citation

Baranov A.A., Olivio A.P. Coplanar multi-turn rendezvous in near-circular orbit using a low-thrust engine. RUDN Journal of Engineering Research. 2022;23(4):283-292. http://doi.org/10.22363/2312-8143-2022-23-4-283-292

Компланарная многовитковая встреча на околокруговой орбите с помощью двигателей малой тяги

А.А. Баранов3 , А.П. Оливиоь

^Институт прикладной математики имени М.В. Келдыша РАН, Москва, Российская Федерация ^Российский университет дружбы народов, Москва, Российская Федерация ^ [email protected]

История статьи Аннотация. Описан алгоритм, позволяющий рассчитать параметры

Поступила в редакцию: 27 сентября 2022 г. маневров, исполняемых на нескольких витках двигателем малой тяги,

Доработана: 25 ноября 2022 г. обеспечивающих перелет активного космического аппарата в заданную

Принята к публикации: 28 ноября 2022 г. окрестность целевого космического объекта. Движение происходит в

© Baranov A.A.. Olivio A.P.. 2022

|(cc) © © I This work is licensed under a Creative Commons Attribution 4.0 International License https://creativecommons.Org/licenses/by-nc/4.0/legalcode

Ключевые слова:

космический аппарат, расчет параметров маневров, космический объект, малая тяга, круговая орбита, импульс скорости

окрестности круговой орбиты. При решении задачи используются линеаризованные уравнения движения. Влияние нецентральности гравитационного поля и атмосферы не учитываются. Определение параметров маневров происходит в три этапа. На первом и третьем этапах параметры импульсного перехода и перехода, выполняемого двигателем малой тяги, определяются аналитически. На втором этапе распределение маневрирования между витками, обеспечивающее решение задачи встречи, осуществляется перебором по одной переменной. Данный метод решения задачи прост и гарантирует высокую надежность определения параметров маневров, что позволяет использовать его на борту космического аппарата. Исследуется зависимость суммарной характеристической скорости решения задачи встречи от числа витков перелета и величины тяги двигателя.

Для цитирования

Baranov A.A., Olivio A.P. Coplanar multi-tum rendezvous in near-circular orbit using a low-thrust engine // Вестник Российского университета дружбы народов. Серия: Инженерные исследования. 2022. Т. 23. № 4. С. 283-292. http://doi.org/10.22363/2312-8143-2022-23-4-283-292

Introduction

The problem of the meeting in a near-circular orbit using a low-thrust engines is of great importance in the practice of spacecraft (SC) flights. This problem is solved during the rendezvous and docking of spacecraft, the implementation of a group flight of several spacecraft, the formation of a given configuration of satellite systems, the removal of space debris, and maintenance of the spacecraft.

Over the past sixty years, the problem of the meeting has been considered in the papers of many authors. Particular attention was paid to the problem of rendezvous in near-circular orbits, when maneuvers are performed by high-thrust engines. A problem of this type was encountered most often in practical paper.

As the first notable papers in this area, we can note the papers of J.E. Prussing [1], who studied the problem of a meeting with a duration of no more than three turns for the case of two circular coplanar orbits, and J.-P. Marec [2], who solved the classical problem of a meeting of average duration in near-circular orbits.

Currently, three main approaches to solving complex multi-impulse problems of spacecraft maneuvering are widely used. In the first case, the problem of maneuvering in the plane of the orbit and the problem of turning the plane of the orbit are solved independently of each other. A similar scheme was used, in particular, to implement the Shuttle rendezvous with the orbital station1, to control the move-

1 Shuttle Press Kit: STS-92. Available from: https://hist0ryc0llecti0n.jsc.nasa.g0v/JSCHist0ryP0rtal/hist0ry/s huttle_pk/pk/Flight_100_STS-092_Press_Kit.pdf (accessed: 17.08.2022).

ment of geostationary satellites [3], satellites included in satellite systems [4], etc. The advantage of such a scheme is its simplicity and reliability, the disadvantage is the excessive costs of the total characteristic speed for maneuvering.

In the second case, numerical methods are used to find the optimal solution to the most complex multi-impulse problems, taking into account a wide range of restrictions [5; 6]. Most often, the simplex method is used to calculate the parameters of maneuvers [7; 8].

In the third method, at the first stage, using the solution of the Lambert problem, the parameters of the two-impulse solution of the meeting problem are determined. Then the behavior of the hodograph of the basis vector corresponding to the found solution is analyzed, and, if necessary, additional velocity impulses are added to obtain the optimal solution. This approach was first applied in the papers of Lion and Handelsman [9], Jezewski and Roosendaal [10].

There are also methods that are at the junction of different approaches. For example, in papers [11; 12] proposed numerical-analytical methods for solving the multi-impulse meeting problem, combining the advantages of the first and second of the previously listed approaches. They allow using the results obtained in the early papers of T. Edelbaum [13], J.-P. Marec [2], when solving modern practical problems.

Since the 1960s, the process of using electric rocket engines (ERE) on spacecraft has begun. Due to the high specific impulse, EREs can significantly reduce fuel costs for orbital maneuvering. However, the small (compared to traditional liquid rocket engines) thrust of the ERE leads to the need to take into account their long-term operation.

Problems of this type occupy a special place among the problems of optimal spacecraft maneuvering. A significant number of papers have been devoted to them, several very interesting monographs have been published [14; 15]. Particularly noteworthy are the papers of V.G. Petukhov [16], Petukhov and Olivio [17] and [18; 19]. Due to the complexity of the problems in which it is assumed that maneuvering is carried out using low-thrust propulsion systems, they have traditionally been solved by numerical methods using the Pontryagin maximum principle or the continuation method. In recent years, to solve problems with a large extent of maneuvers, Yu.P. Ulyby-shev successfully uses the interior point method [20].

In the method considered in this paper, the meeting problem is solved both in an impulse setting and taking into account the long-term operation of a low-thrust engine.

To analyze the relative motion of a spacecraft in the vicinity of circular orbits, it is necessary to use special mathematical models of motion. The most popular mathematical model of the relative motion of spacecraft in the vicinity of circular orbits is the Hill -Clohessy - Wiltshire (HCW) model. Linearized differential equations for the relative motion of a spacecraft in the vicinity of a circular orbit for the rendezvous and docking problem were obtained by Clohessy -Wiltshire in 1960 [21], but as early as the 19th century, similar equations were used by Hill in his theory of the motion of the Moon [22]. In this mathematical model, to obtain the equations of relative motion, a rotating (orbital) coordinate system and linearization of differential equations of relative motion are used, based on the assumption that the distance between the considered spacecraft is small compared to the average radius of the orbit. In this paper, we use linearized equations obtained by P.E. Elyasberg [23].

Due to the increase in the number of maneuvering SC, the increase in the efficiency of solving problems, at present there is a tendency to transfer the process of calculating maneuvers on board the spacecraft. This leads to the need to simplify the process of calculating the parameters of maneuvers and increase.

1. Formulation of the meeting problem

The problem of calculating the parameters of flight maneuvers between close near - circular orbits is solved in an approximate impulse formulation, within the framework of the unperturbed Keplerian motion.

The flight conditions using N velocity impulse in a fixed time from the initial orbit to a given point

of the final orbit (meeting problem) in linear approximation can be written as [12]:

N

2 (AVrisin9, +2A Vti cos ^)=Aex; (1)

i= 1

N

2 (-A Vri cos 91 +2 A Visingi )=Ae^; (2)

i=1

N

2 2 A V.=A a; (3)

i=1

N

2 (2AVri(1 - cos91)+AV, x

i=1

x(-39i +4sin 9 i ))=At, (4)

where Aex = e/cosro/- eocosroo; Aey = e/sinro/ - eosinroo; Aa = (a/- ao) / ro; At = Xo(t/- to); AVti = A Vt* / Vo;

AVri = A Vr* / Vo; / o - indices corresponding to

the final and initial orbits; e/, eo - eccentricities of the orbits; a/, ao - semi-major axes of the orbits; ro/, roo - angles between the direction to the pericenter of the corresponding orbit and the direction to the point specified on the final orbit (the axis Ox is directed to this point); t/ - the necessary arriving time to this point; to - time at which, when moving in the initial orbit, the projection of the radius-vector in the plane of the final orbit falls on a beam passing through a given meeting point; Vo, ^o - orbital and angular velocities along the reference circular orbit of radius ro (ro = a/); N - number of velocity impulse; 9 - angle of application of the i-th velocity impulse, calculated from the direction to the given meeting point towards the motion of the SC; A Vt*, A Vr* - transversal and radial components of the i-th velocity impulse respectively.

It must be taken into consideration that the angles 9 - are negative, i.e., because it was assumed that 9/ = o at the given point.

The problem of finding the parameters of optimal maneuvers can be formulated as follows: it's necessary to determine AVri, AVti, 9 (i = 1, ..., N), at which the total characteristic velocity of maneuvers AV is minimal

N

AV

¿=i ¿=i

при ограничениях (1)-(4).

2. Algorithm for solving the transfer problem

When solving the problem of transfer between coplanar orbits, the first three equations of system (1)-(4) are used.

There are three types of solutions for which the necessary optimality conditions are satisfied:

a) on the ц, A plane, the hodograph of the basis vector is an ellipse, the center of which is located on the ц axis, but is shifted from the origin of the coordinate system; the ellipse is tangent to a circle of unit radius at a point on the ц axis;

b) the hodograph of the basis vector degenerates into a point coinciding with the point of intersection of a circle of unit radius with the ц axis;

c) the hodograph of the basis vector - an ellipse centered at the origin of the coordinate system, touching the circle at two points on the ц axis.

Since all possible optimal solutions have решений A = 0, and ц Ф 0, then the velocity impulses of these solutions are purely transversal.

Assuming that the velocity impulses are applied at the optimal points for correcting the eccentricity vector,

t Ae.-tg ф- ;

Ф1 = фе; Ф2 = Ф1 + п

a vlt=V\ti ;

i=1

N

av2 t=£av2 й,

(7)

(8)

i=1

where N is the number of turns on which maneuvering is allowed.

The further goal is to choose such a distribution of the magnitudes of the velocity impulses over the coils so that equation (4) is satisfied.

To significantly simplify the solution of the problem, we assume that the magnitudes of the velocity impulses along the turns change linearly:

AVlft =AVm + +(A Vm -A Vm)(i -1)/( N -1); (9)

AV2 tl = A V211 + +(A V2tN -A V211 )(i -1)/( N -1), (10)

where AVm, AVun and AV2i1, AV2tN are the value of the velocity of impulses on the first and last allowed turns of maneuvering, which is part of the first and second velocity of solving the transfer problem.

Substituting the values of the velocity impulses calculated by formulas (9), (10) into (7) and (8) we get:

A Vit=X A Vw = 0.5 N ( A Viti + a Vin ); (11)

we find the magnitude of the velocity impulses of the optimal solution:

AVt1 = 4(Aa + Ae).

AVt2 = -4(Aa-Ae).

(5)

(6)

3. Algorithm for solving the meeting problem

When solving the meeting problem, the values

of the velocity impulses AVa, AVq, determined

when solving the transition problem, are distributed

among N turns allowed for maneuvering:

AV21=XAV2ft. = 0.5 N (A V2t1 +AV2N ). (12)

i=1

Using (11) and (12), we obtain formulas for determining A Vin, AV2n:

AV

AV = 1t -AV • 1N 0.5N 1t1;

(13)

(14)

Substituting the found values AViw, AV2tN into formulas (9) and (10), we obtain:

i=1

àVlti = 2AV1(i - 1)/ /N(N-l)+AVia[l-^]; (15)

AV2ti = 2AV2(i - 1)/ /N(N-1)+AV2tl[l-2-^]. (16)

Thus, we found the values of all impulses of velocity, expressed only through AVm and AV^i. Substituting them into equation (3), we obtain a linear equation with two unknowns AVm, AV2ii. The coefficients for velocity impulses are known, since their application angles are known:

^ = 9e + N - N); (17)

impulse one and we consider that the problem has been solved. If the duration of the maneuvers is significant, then we proceed to the solution with low thrust.

4. Solving the problem with "low thrust"

For each turn, we find what changes in the eccentricity and the semi-major axis produce velocity impulses determined on this turn

Ae, =2AVltl - 2AV2i.; (20)

Aa, =2AVU + 2AV2tI. (21)

Then we determine the required duration of low-thrust maneuvers, which will produce the same change in these elements:

Ф2i = фe + 2я(Ni - N).

(18)

Going through the value of the variable AVm, within the given limits, for each value from equation (3) we find the value of the variable AVm.

Then, using (15) and (16), we find the values of all velocity impulses. Adding the modules of all velocity impulses, we find the total characteristic velocity of the next solution. The solution for which the total characteristic velocity is minimal is taken as a solution to the meeting problem. If the total characteristic velocity of the found solution coincides with the total characteristic velocity of the solution to the transition problem, then a solution was found with the minimum possible total characteristic velocity.

At the next stage, we estimate the duration of each of the found maneuvers.

The duration of each of the maneuvers is estimated using the ratio:

w

Дф. ДV w

i ■>

(19)

where wc - the centripetal acceleration of the refer-

ence circular orbit

V

2 Л

; w - the acceleration

'О у

created by the propulsion system f w = p

m

m -

the mass of the active SC; P - thrust of its engine.

If the duration of the largest velocity impulse doesn't exceed 20°, then the solution is close to the

4sin

ДФ1

■4sin

ДФ2г

w Де.

w

2Дфь. + 2Дф2, =

wc Да.

(22)

w

From system (22) one can find the quantities

Дф1., Дф2. [24]:

. w Да _

Дф1 = —--h 2arcsin-

4wn

w„ Де

8wn cos

wc Да '

8wn

Дф2 =

wc Да 4 wn

- 2arcsin-

wc Де

8wn cos

wcДа ' 8wn

(23)

Thus, turn by turn, we find the duration of all maneuvers. The problem with low thrust is correctly solved. If the argument of the arcsine is greater than 1, then there is no solution (with the available thrust and mass of the spacecraft, it is impossible to solve the rendezvous problem for a given number of turns).

The found solution with "low thrust" gives the same change in the semi-major axis and the eccentricity vector as the original impulse solution. Equation (4) is also quite accurate, since the middle of long maneuvers coincides with the moments of ap-

plication of velocity impulses, and the same change in the major axis on the turn is made, therefore, the required time of arrival at the meeting point is provided.

5. Iterative procedure

In the formulated meeting problem, linearized equations of motion are used, the non-centrality of the gravitational field, the influence of the atmosphere, etc. are not taken into account. This leads to the fact that the actual accuracy of the fulfillment of the terminal conditions in the system (1)-(6) will be insufficient. To solve the problem with a given accuracy, you can use the iterative scheme [5; 6], which consists of the following steps:

1. At the beginning of the next iteration, the "approximate" problem is solved: under the previously accepted simplifying assumptions, the parameters of the maneuvers that ensure the formation of the "target" orbit are determined (at the first iteration, the "target" orbit coincides with the final orbit).

2. Then, taking into account the calculated velocity impulses, using the models of all necessary disturbances, an "accurate" prediction of the spacecraft motion is carried out and the parameters of the formed orbit are found.

3. The deviations of the parameters of the formed orbit from the corresponding parameters of the final orbit are calculated.

4. If the deviations exceed the allowable values, then the parameters of the "target" orbit are changed by the value of the calculated deviations, and the next iteration is carried out.

5. The procedure ends when the terminal conditions are met with the specified accuracy.

6. For "accurate" forecasting, as a rule, numerical and/or high-precision numerical-analytical integration is used. It is possible to use different forecasting methods at different iterations, but the accuracy of the forecast should increase with the number of the current iteration.

7. Numerical integration takes into account the influence of non-centrality of the gravitational field, atmosphere, light pressure, etc., the operation of the spacecraft engines is accurately modeled, therefore, in spite the fact that the maneuver parameters and are found at each iteration using the simplest motion

model, but as a result of the iterative procedure, they provide access to the final orbit with the required accuracy.

6. Problem solving examples

Considering the motion of a spacecraft (SC) relative to the point O, moving along a non-perturbed near-circular orbit with a radius of 6871 km. Let's take the gravitational parameter of the Earth equal to 3.9860044-1014 m3/s2. Consider the problem of a flight with the help of N impulses of velocity in a fixed time from the initial orbit to a given point of the final orbit from a point in phase space r0 = (10, 100, -5) km, V0 = (1, -10, 3) m/s to the origin, i.e, to the point rf = (0, 0, 0) km, with a velocity vf = (0, 0, 0) m/s. For the problem, we'll take the initial mass of the SC equal to 1000 kg, the specific impulse of the SC propulsion system is 220 seconds (2157.463 m/s), and the thrust (T) will be varied in the range from 0.362 to 100 N. The flight is carried out in N = 4 and 13 turns.

Double impulse transfer. In the Table 1 shows the results of the calculation the parameters of a two-impulse transfer between coplanar orbits, that is, the magnitude of the transversal components of the velocity impulses, the angles of application of the first and second impulses, as can be seen, the minimum value of the characteristic velocity that a SC must have for the transfer maneuver, is 4.485 m/s.

Multi-impulse solution. The value of the first velocity impulse is moved within the limits from -2.785 to 0.5 m/s with a step of 0.024 m/s.

Parameters of the optimal solutions. In the Tables 2-4 show the parameters of the optimal solution, that is, the values of the velocity impulses on the turns for the cases when N = 4 and 13.

Solution of the problem with low thrust. In the Tables 5, 6 and 7 shown the results of calculating the problem with low thrusts "Thrust = 1 and Thrust = 0.362", that is, the values of the velocity impulses and the duration of maneuvers on turns for cases when N=4 and 13. In some cases, there're no solutions because the value of the arcsine argument is out of range (-1; 1). With increased thrust, the duration of the maneuvers is reduced, and the costs of the total characteristic velocity (CXC) of the low thrust solution during thrust increment, coincides with costs of the CXC of the impulse solution.

Table 1

Results of the calculation the parameters of the coplanar transfer problem

AV1,m/s AV2,m/s |AK|,m/s <f>e° °

-2.785 1. 7 4.485 6.4 186.4 366.4

Table 2

Parameters of the optimal solution for N = 4

N AVU, m/s AV2i,m/s IAVltl + IAV2il,m/s

1 -0.024 0.848 0.872

2 -0.472 0.566 1.038

3 -0.92 0.284 1.204

4 -1.369 0.002 1.371

I -2.785 1.7 4.485

Table 3 Parameters of the optimal solution for N = 13

N AVzi,m/s |ДК„| + |ДК2!|,т/8

1 -0.001 0.199 0.2

2 -0.037 0.187 0.224

12 -0.392 0.074 0.466

13 -0.427 0.063 0.49

-2.785 1.7 4.485

Table 4

Results of the calculation of the problem with low thrust, for N = 4

T, N 0.362 0.37 0.4 0.5 1 2 5 10 100

Дф° There is not solution There is not solution There is not solution There is not solution 300.137 144.199 57.082 28.499 2.849

AV,m/s There is not solution There is not solution There is not solution There is not solution 4.726 4.541 4.494 4.487 4.485

Table 5

Results of the calculation of the problem with low thrust "Thrust = 1", for N = 4

N AVu,m/s AV2i,m/s | АУ^ +^V^m/s Лф2,°

1 -0.04 0.864 0.904 -2.561 54.87 57.431

2 -0.48 0.573 1.053 -30.473 36.424 66.897

3 -0.942 0.306 1.248 -59.841 19.434 79.275

4 -1.443 0.077 1.52 -91.65 4.884 96.534

I -2.905 1.82 4.725 -184.525 115.612 300.137

Table 6

Results of the calculation of the problem with low thrust, for N = 13

T, N 0.362 0.37 0.4 0.5 1 2 5 10 100

Дф° 809.865 791.329 728.952 578.15 285.884 142.555 56.979 28.487 2.849

AV,m/s 4.616 4.61 4.591 4.551 4.501 4.489 4.486 4.485 4.485

Table 7

Results of the calculation of the problem with low thrust "Thrust = 0.362", for N = 13

N AVlt,m/s AV2i,m/s | AV^ + |AV2i|,m/s (| AV■Li| + | AV2i|)°

1 -0.003 0.2 0.203 -0.45 35.149 35.599

2 -0.038 0.189 0.227 -6.642 33.117 39.759

12 -0.405 0.087 0.492 -71.015 15.253 86.268

13 -0.444 0.08 0.524 -77.953 13.967 91.92

I -2.85 1.765 4.615 -500.114 309. 751 809.865

Conclusion

In the paper describes an algorithm for calculating the parameters of a multi-turn multi-impulse encounter. The main advantage of the proposed algorithm is its simplicity and reliability, which makes it possible to use it not only in ground control centers, but also on board the spacecraft. Simultaneously, this algorithm makes it possible to obtain the optimal solution of the problem in the case when the initial phase belongs to the optimal phase range and the total characteristic velocity of the solution of the meeting problem coincides with the total characteristic velocity of the optimal solution of the transfer problem. The algorithm makes it possible to obtain a solution even when the maneuvers are performed by low thrust engines. The examples given in the paper confirm the ope-rability of this algorithm and the high quality of the resulting solution.

References

1. Prussing JE. Optimal two- and three-impulse fixed-time rendezvous in the vicinity of a circular orbit. AIAA Journal. 1970;8(7):46-56. https://doi.org/10.2514/3.5876

2. Marec JP. Optimal space trajectories. Studies in Astronautics (vol. 1). Amsterdam, Oxford, New York: Elsevier Sci. Pub. Co.; 1979.

3. Bulynin YuL. Ballistic support for orbital motion control of geostationary spacecraft at various stages of operation. System Analysis, Control and Navigation: Abstracts of Reports. Crimea, Yevpatoria; 2008. p. 73-74. (In Russ.)

4. Rylov YuP. Control of a spacecraft entering the satellite system using electric rocket engines. Kosmiches-kie Issledovaniya. 1985;23(5):691-700.

5. Kulakov AYu. Model and algorithms of reconfiguration of the spacecraft motion control system (dissertation of the candidate of Technical Sciences). St. Petersburg; 2017. (In Russ.)

6. Baranov AA. Algorithm for calculating the parameters of four-impulse transitions between close almost-circular orbits. Cosmic Research. 1986;24(3):324-327.

7. Lidov ML. Mathematical analogy between some optimal problems of trajectory corrections and selection of measurements and algorithms of their solution. Kos-micheskie Issledovaniya. 1971;9(5):687-706. (In Russ.)

8. Gavrilov V, Obukhov E. Correction problem with fixed number of impulses. Kosmicheskie Issledovaniya. 1980;18(2):163-172. (In Russ.)

9. Lion PM, Handelsman M. Basis-vector for pulse trajectories with a given flight time. Rocket Technology and Cosmonautics. 1968;6(1):153-160. (In Russ.)

10. Jezewski DJ, Rozendaal HL. An efficient method for calculating optimal free-space n-impulse trajectories. AIAA Journal. 1968;6(11):2160-2165. (In Russ.)

11. Baranov AA. Geometric solution of the problem of a rendezvous on close nearly circular coplanar orbits.

Cosmic Research. 1989;27(6):689-697.

12. Baranov AA, Roldugin DS. Six-impulse maneuvers for rendezvous of spacecraft in near-circular nonco-planar orbits. Cosmic Research. 2012;50(6):441-448.

13. Edelbaum TN. Minimum impulse transfer in the vicinity of a circular orbit. Journal of the Astronauti-cal Sciences. 1967;XIV(2):66-73.

14. Lebedev VN. Calculation of the motion of a spacecraft with low thrust. Moscow: Publishing House of the USSR Academy of Sciences, 1968. (In Russ.)

15. Grodzovsky GL, Ivanov YuN, Tokarev VV. Mechanics of low-thrust space flight. Moscow: Nauka Publ.; 1966. (In Russ.)

16. Petukhov VG. Method of continuation for optimization of interplanetary low-thrust trajectories. Cosmic Research. 2012;50(3):249-261. (In Russ.)

17. Petukhov VG, Olivio AP. Optimization of the finite-thrust trajectory in the vicinity of a circular orbit.

Advances in the Astronautical Sciences. 2021;174:5-15.

18. Baranov AA. Maneuvering in the vicinity of a circular orbit. Moscow: Sputnik+ Publ.; 2016. (In Russ.)

19. Baranov AA. Development of methods for calculating parameters of spacecraft maneuvers in the vicinity of a circular orbit (dissertation of the Doctor of Physical and Mathematical Sciences). Moscow; 2019. (In Russ.) Available from: http://library.keldysh.ru/diss.asp?id=2019-baranov (accessed: 12.08.2022).

20. Ulybyshev YuP. Optimization of multi-mode rendezvous trajectories with constraints. Cosmic Research. 2008;46(2):133-145. (In Russ.)

21. Clohessy WH, Wiltshire RS. Terminal guidance system for satellite rendezvous. Journal of the Aerospace Sciences. 1960;27(9):653-678. https://doi.org/10.2514Z8.8704

22. Hill GW. Researches in lunar theory. American Journal of Mathematics. 1878;1:5-26.

23. Elyasberg PE. Introduction to the theory of flight of artificial Earth satellites. Moscow: Nauka Publ.; 1965. (In Russ.)

24. Baranov AA, Prado AFB, Razumny VY, Bara-nov Jr. Optimal low thrust transfers between close near-circular coplanar orbits. Cosmic Research. 2011;49(3): 269-279. https://doi.org/10.1134/S0010952511030014

Список литературы

1. Prussing J.E. Optimal two- and three-impulse fixed-time rendezvous in the vicinity of a circular orbit // AIAA Journal. 1970. Vol. 8. Issue 7. Pp. 46-56. https://doi.Org/10.2514/3.5876

2. Marec J.P. Optimal space trajectories // Studies in Astronautics. Amsterdam, Oxford, New York: Elsevier Sci. Pub. Co., 1979. Vol. 1. 329 p.

3. Булынин Ю.Л. Баллистическое обеспечение управления орбитальным движением геостационарных КА на различных этапах эксплуатации //13 Международная научная конференция «Системный анализ, управление и навигация»: тезисы докладов. Крым, Евпатория, 2008. С. 73-74.

4. Рылов Ю.П. Управление космическим аппаратом, входящим в спутниковую систему при помощи электроракетных двигателей // Космические исследования. 1985. Т. 23. № 5. C. 691-700.

5. Кулаков А.Ю. Модель и алгоритмы реконфигурации системы управления движением космического аппарата: дис. ... канд. техн. наук. СПб., 2017. 156 с.

6. Baranov A.A. Algorithm for calculating the parameters of four-impulse transitions between close almost-circular orbits // Cosmic Research. 1986. Vol. 24. Issue 3. Pp. 324-327.

7. Лидов М.Л. Математическая аналогия между некоторыми оптимальными задачами коррекции траекторий и выбора состава измерений и алгоритмы их решения // Космические исследования. 1971. Т. 9. № 5. С. 687-706.

8. Гаврилов В.П., Обухов Е.В. Задача коррекции с ограничением на число импульсов // Космические исследования. 1980;18(2):163-172.

9. Лайон П.М., Хенделсмен М. Базис-вектор для импульсных траекторий с заданным временем перелета // Ракетная техника и космонавтика. 1968;6(1):153-160.

10. Ежевски Д.Дж., РозендаалХ.Л. Эффективный метод расчета оптимальных N-импульсных траекторий полета в космическом пространстве // Ракетная техника и космонавтика. 1968. Т. 6. № 11. С. 138-145.

11. Baranov A.A. Geometric solution of the problem of a rendezvous on close nearly circular coplanar orbits // Cosmic Research. 1989. Vol. 27. Issue 6. Pp. 689-697.

12. Baranov A.A., Roldugin D.S. Six-impulse maneuvers for rendezvous of spacecraft in near-circular noncoplanar orbits // Cosmic Research. 2012. Vol. 50. Issue 6. Pp. 441-448.

13. Edelbaum T.N. Minimum impulse transfer in the vicinity of a circular orbit // Journal of the Astronau-tical Sciences. 1967. Vol. XIV. Issue 2. Pp. 66-73.

14. Лебедев В.Н. Расчет движения космического аппарата с малой тягой. М.: ВЦ АН СССР, 1968. 108 с.

15. Гродзовский Г.Л., Иванов Ю.Н., Токарев В.В. Механика космического полета с малой тягой. М.: Наука, 1966. 678 с.

16. Петухов В.Г. Метод продолжения для оптимизации межпланетных траекторий с малой тягой // Космические исследования. 2012. Т. 50. № 3. С. 258-270.

17. Petukhov V.G., Olivio A.P. Optimization of the finite-thrust trajectory in the vicinity of a circular orbit // Advances in the Astronautical Sciences. 2021. Vol. 174. Pp. 5-15.

18. Баранов А.А. Маневрирование в окрестности круговой орбиты. М.: Спутник+, 2016. 512 с.

19. Баранов А.А. Разработка методов расчета параметров маневров космических аппаратов в окрестностях круговой орбиты: дис. ... д-р физ.-мат. наук. М., 2019. 304 с. URL: http://library.keldysh.ru/diss.asp?id=2019-baranov (дата обращения: 12.08.2022).

20. Улыбышев Ю.П. Оптимизация межорбитальных перелетов с малой тягой при ограничениях // Космические исследования. 2012. Т. 50. № 5. С. 403-418.

21. Clohessy W.H., Wiltshire R.S. Terminal guidance system for satellite rendezvous // Journal of the Aerospace Sciences. 1960. Vol. 27. Issue 9. Pp. 653-678. https://doi.org/10.2514/8.8704

22. Hill G. W. Researches in lunar theory // American Journal of Mathematics. 1878. Vol. 1. Pp. 5-26.

23. Эльясберг П.Е. Введение в теорию полета искусственных спутников Земли. М.: Наука, 1965. 540 с.

24. Баранов А.А., Прадо А.Ф.Б., Разумный В.Ю., Баранов А.А. Оптимальные переходы с малой тягой между близкими околокруговыми компланарными орбитами // Космические исследования. 2011. Т. 49. № 3. С. 278-288.

About the authors

Andrey A. Baranov, Candidate of Physical and Mathematical Sciences, leading researcher, Institute for Applied Mathematics, Russian Academy of Sciences, 4 Miusskaya Ploshchad', Moscow, 125047, Russian Federation; ORCID: 0000-0003-18239354, Scopus Author ID: 14059618600, eLIBRARY SPIN-code: 6606-3690; [email protected]

Adilson Pedro Olivio, postgraduate, Department of Mechanics and Control Processes, Academy of Engineering, Peoples' Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya St, Moscow, 117198, Russian Federation; ORCID: 0000-0001-5632-3747, Scopus Author ID: 57223139496, eLIBRARY SPIN-code: 7628-6084; [email protected]

Сведения об авторах

Баранов Андрей Анатольевич, кандидат физико-математических наук, ведущий научный сотрудник, Институт прикладной математики имени М.В. Келдыша, Российская академия наук, Российская Федерация, 125047, Москва, Миусская пл., д. 4; ORCID: 0000-0003-1823-9354, Scopus Author ID: 14059618600, eLIBRARY SPIN-код: 6606-3690; [email protected]

Оливио Адильсон Педро, аспирант, департамент механики и процессов управления, Инженерная академия, Российский университет дружбы народов, Российская Федерация, 117198, Москва, ул. Миклухо-Маклая, д. 6; ORCID: 00000001-5632-3747, Scopus Author ID: 57223139496, eLIBRARY SPIN-код: 7628-6084; [email protected]