Доклады Российской академии наук. Физика, технические науки, 2021, T. 499, № 1, стр. 66-72

ПОДХОД ДЛЯ ПОСТРОЕНИЯ ДИНАМИЧЕСКИХ МОДЕЛЕЙ ПРОЦЕССА РАСКРЫТИЯ ТРАНСФОРМИРУЕМЫХ КОСМИЧЕСКИХ КОНСТРУКЦИЙ

В. Н. Бакулин 1*, С. В. Борзых 2**

1 Институт прикладной механики Российской академии наук
Москва, Россия

2 ПАО Ракетно-космическая корпорация “Энергия” им. С.П. Королева
Королев, Московская область, Россия

* E-mail: vbak@yandex.ru
** E-mail: rigidbor@gmail.com

Поступила в редакцию 17.04.2021
После доработки 17.04.2021
Принята к публикации 22.04.2021

Полный текст (PDF)

Аннотация

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

Ключевые слова: универсальный подход, динамические модели, раскрываемые космические конструкции, структурно-сложная механическая система, солнечные батареи, уравнения связей

ВВЕДЕНИЕ

Создание современных летательных аппаратов (ЛА) требует решения многих актуальных проблем [1, 2]. В настоящее время существует устойчивая тенденция постоянного расширения области применения размещенных в космосе объектов – орбитальных станций, спутников, спутниковых группировок и т.д. Кроме ставших традиционными прикладных направлений – связи, мониторинга земной поверхности, космической биологии и медицины, космического материаловедения [3] и ряда других, появляются новые направления, нацеленные на решение перспективных научных и прикладных задач, таких как исследования дальнего космоса [4], глобальная энергетика и экология [5, 6], освоение внеземных ресурсов. Космические аппараты (КА) или целые комплексы (рис. 1) в связи с этим приобретают все более сложную структуру, для которой характерным является наличие крупногабаритных космических трансформируемых конструкций – различного типа антенн, выдвижных штанг, солнечных батарей и т.д. [79].

Рис. 1.

Международная космическая станция.

Такого рода трансформируемые или раскрываемые конструкции на этапе выведения ракетой-носителем, межорбитальным буксиром или разгонным блоком должны находиться в компактном сложенном транспортировочном состоянии и трансформироваться в рабочее конечное положение после окончания выведения, на целевой орбите.

Перспективные схемы раскрытия трансформируемых космических конструкций, проблемы математического моделирования динамики процесса их раскрытия рассмотрены в работе [7], моделирование безопасного отделения пилотируемых транспортных кораблей от орбитальной станции – в работе [10], исследование процессов отделения малых космических аппаратов от транспортно-пускового контейнера, установленного на грузовом корабле “Прогресс МС”, как во время автономного полета корабля, так и в период его пребывания в составе Международной космической станции, – в [11], алгоритмы отделения космического корабля “Союз МС” от нестабилизированной Международной космической станции – в [12].

Совершенно очевидно, что отказ системы раскрытия по какой-либо причине (непрохождение команды на раскрытие, отказ пиротехники, отказ механизмов и т.д.) практически неизбежно означает невыполнение целевой задачи миссии.

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

ПОСТАНОВКА ЗАДАЧИ

Наиболее распространенным подходом при моделировании движения сложных механических систем, который применялся, в том числе и авторами, является использование динамических уравнений, основанных на классических формах Эйлера, Лагранжа и Аппеля. Это позволяет ограничить размерность системы линейных уравнений, решаемых на каждом шаге интегрирования, числом обобщенных координат механической системы [13]. Однако при таком подходе неизбежен весьма трудоемкий этап составления для механической системы в целом сложных исходных выражений для момента количества движения (метод Эйлера), кинетической энергии (метод Лагранжа) и “энергии ускорений” (метод Аппеля), которые жестко привязаны к структуре механической системы, что исключает возможность оперативного формирования расчетной модели системы с меняющейся структурой и использование ранее составленных исходных выражений даже для систем с незначительно измененной структурой.

В то же время процесс раскрытия крупногабаритных трансформируемых конструкций может осуществляться в несколько последовательных фаз, и для моделирования раскрытия такого рода конструкций применение описанных выше подходов представляется нецелесообразным.

Принципиальным недостатком такого подхода также является и то, что использование классических форм динамических уравнений возможно только при отсутствии трения в механической системе. Кроме того, например, при использовании уравнений Лагранжа 2-го рода невозможно напрямую определить реакции в узлах сочленения элементов трансформируемой конструкции, а это является одной из основных целей анализа процесса.

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

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

ПОДХОД ДЛЯ ПОСТРОЕНИЯ ДИНАМИЧЕСКИХ МОДЕЛЕЙ ПРОЦЕССА РАСКРЫТИЯ ТРАНСФОРМИРУЕМЫХ КОСМИЧЕСКИХ КОНСТРУКЦИЙ

Рассматривается механическая система, включающая корпус КА, штанги и панели солнечных батарей (СБ) (рис. 2). Для каждого из тел системы записываются уравнения движения центра масс и уравнения вращения относительно него на основании законов изменения количества движения и кинетического момента, а для определения силовых факторов, обусловленных взаимодействием отдельного тела с соседними телами, записываются специальные уравнения связей, сводимые в конечном итоге к системе линейных уравнений с переменными по времени коэффициентами.

Рис. 2.

Пример общего вида космического аппарата с солнечными батареями: слева – рабочее положение; справа – транспортировочное положение.

Поскольку связь между отдельными элементами СБ выполнена в виде оси вращения, то действие на данную панель соседних “отброшенных” элементов учитывается в виде сил реакций (три компоненты вектора сил) и моментов (две компоненты, поскольку связь допускает разворот относительно общей оси x) (рис. 3).

Рис. 3.

Расчетная схема.

Уравнения поступательного перемещения центра масс каждого тела в инерциальных осях и уравнения вращательного движения относительно центра масс в связанных осях соответственно будут иметь вид

${{m}_{i}}{{a}_{{oi}}} = {{F}_{{oi}}} + {{R}_{i}},$
$\left[ {{{I}_{{oi}}}} \right]\frac{{\tilde {d}{{\omega }_{i}}}}{{dt}} + {{\omega }_{i}} \times \left( {\left[ {{{I}_{{oi}}}} \right]{{\omega }_{i}}} \right) = {{L}_{{oi}}} + {{L}_{i}} + {{L}_{{ci}}},$
где aoi, ${{\varepsilon }_{i}} = \frac{{\tilde {d}{{\omega }_{i}}}}{{dt}}$ – соответственно ускорение центра масс и угловое ускорение i-го тела системы, локальная производная $\frac{{\tilde {d}{{\omega }_{i}}}}{{dt}}$ берется в его связанной системе координат, mi, $\left[ {{{I}_{{oi}}}} \right]$ – соответственно масса и тензор инерции i-го тела, Foi – главный вектор внешних сил, включая силы, создаваемые средствами раскрытия батарей, ${{R}_{i}}$ – главный вектор сил реакций связей, действующих на i-е тело, Loi – главный момент внешних сил, ${{L}_{i}}$ – главный момент сил реакций связей, ${{L}_{{ci}}}$ – главный реактивный момент.

Отсюда можно выразить ${{a}_{{oi}}}$ и ${{\varepsilon }_{i}}$:

(1)
$\begin{gathered} {{a}_{{oi}}} = m_{i}^{{ - 1}}\left( {{{F}_{{oi}}} + {{R}_{i}}} \right), \\ {{\varepsilon }_{i}} = {{[{{I}_{{oi}}}]}^{{ - 1}}}\{ - {{\omega }_{i}} \times \left( {\left[ {{{I}_{{oi}}}} \right]{{\omega }_{i}}} \right) + {{L}_{{oi}}} + {{L}_{i}} + {{L}_{{ci}}}\} . \\ \end{gathered} $

Результатом интегрирования уравнений (1) будут скорости центров масс тел Voi, координаты центров масс roi в инерциальной системе координат и угловые скорости ${{\omega }_{i}}$ тел.

Для определения реакций связей необходимо записать дополнительно уравнения связи, вид которых определяется конкретным характером связи. Поскольку два последовательно соединенных тела имеют общие точки связи (рис. 3), то радиусы-векторы в инерциальной системе координат и полные скорости этих точек, принадлежащих двум разным телам, должны быть одинаковы. Кроме того, одинаковыми должны быть проекции угловой скорости на оси Yci, Zci системы координат связи.

С учетом изложенного можно записать

(2)
$\begin{gathered} {{{\mathbf{a}}}_{{i,i + 1}}} - {{{\mathbf{a}}}_{{i + 1,i}}} = 0 \\ {{{\mathbf{\varepsilon }}}_{i}} - {{{\mathbf{\varepsilon }}}_{{i + 1}}} + {\mathbf{\omega }} \times {{{\mathbf{\omega }}}_{{i + 1}}} = 0, \\ \end{gathered} $
где ${{a}_{{i,i + 1}}}$, ${{a}_{{i + 1,i}}}$ значения ускорений некоторой связи, расположенной на оси вращения и принадлежащей i-му и (i + 1)-му телам (рис. 3). Первый индекс в индексации ускорений точек связи означает принадлежность к данному телу (звену СБ), а второй – индекс того тела, с которым связано данное тело.

Появление последнего слагаемого во втором из выражений (2) обусловлено вращением системы координат связи OciXciYciZci относительно систем координат i-го и (i + 1)-го тел.

В выражениях (2) полные ускорения точек связи i-го и (i + 1)-го тел запишутся следующим образом:

$a_{{i,i + 1}}^{{}} = {{a}_{{oi}}} + {{\omega }_{i}} \times ({{\omega }_{i}} \times r_{{i,i + 1}}^{{}}) + {{\varepsilon }_{i}} \times r_{{i,i + 1}}^{{}},$
$a_{{i + 1,i}}^{{}} = {{a}_{{oi + 1}}} + {{\omega }_{{i + 1}}} \times ({{\omega }_{{i + 1}}} \times r_{{i + 1,i}}^{{}}) + {{\varepsilon }_{{i + 1}}} \times r_{{i + 1,i}}^{{}}.$

Индексация радиус-векторов (рис. 3) такая же, как индексация ускорений точек связи в выражении (2).

В последнее выражение входят линейные ускорения центров масс тел и угловые ускорения их вращения относительно центра масс. Они определены выражениями (1). Входящие в (1) главные вектора реакций связей ${{R}_{i}}$, моментов реакций связей ${{L}_{i}}$ и реактивных моментов ${{L}_{{ci}}}$ i-го тела (рис. 3) определяются соотношениями:

(3)
$\begin{gathered} {{R}_{i}} = \left[ {{{M}_{{ci \to i}}}} \right] \cdot R_{i}^{c} - \left[ {{{M}_{{ci - 1 \to i}}}} \right] \cdot R_{{i - 1}}^{c}, \\ {{L}_{{ci}}} = \left[ {{{M}_{{ci \to i}}}} \right] \cdot L_{i}^{c} - \left[ {{{M}_{{ci - 1 \to i}}}} \right] \cdot L_{{i - 1}}^{c}, \\ {{L}_{i}} = - [\Phi _{{i,i + 1}}^{{}}] \cdot R_{i}^{c} + [\Phi _{{i,i - 1}}^{{}}] \cdot R_{{i - 1}}^{c}, \\ \end{gathered} $
$\begin{gathered} {{R}_{{i + 1}}} = - \left[ {{{M}_{{ci \to i + 1}}}} \right] \cdot R_{i}^{c} + \left[ {{{M}_{{ci + 1 \to i + 1}}}} \right] \cdot R_{{i + 1}}^{c}, \\ {{L}_{{ci}}} = - \left[ {{{M}_{{ci \to i + 1}}}} \right] \cdot L_{i}^{c} + \left[ {{{M}_{{ci + 1 \to i + 1}}}} \right] \cdot L_{{i + 1}}^{c}, \\ {{L}_{{i + 1}}} = [\Phi _{{i + 1,i}}^{{}}] \cdot R_{i}^{c} - [\Phi _{{i + 1,i + 2}}^{{}}] \cdot R_{{i + 1}}^{c}, \\ \end{gathered} $
в которых компоненты тензоров [Ф] определяются из выражений
$\begin{gathered} \left[ {{{M}_{{ci \to i}}}} \right] \cdot R_{i}^{c} \times r_{{i,i + 1}}^{{}} = [\Phi _{{i,i + 1}}^{{}}] \cdot R_{i}^{c}, \\ \left[ {{{M}_{{ci - 1 \to i}}}} \right] \cdot R_{{i - 1}}^{c} \times r_{{i,i - 1}}^{{}} = [\Phi _{{i,i - 1}}^{{}}] \cdot R_{{i - 1}}^{c}, \\ \end{gathered} $
$\begin{gathered} \left[ {{{M}_{{ci \to i + 1}}}} \right] \cdot R_{i}^{c} \times r_{{i + 1,i}}^{{}} = {\text{[}}\Phi _{{i + 1,i}}^{{}}{\text{]}} \cdot R_{i}^{c}, \\ \left[ {{{M}_{{ci + 1 \to i + 1}}}} \right] \cdot R_{{i + 1}}^{c} \times r_{{i + 1,i + 2}}^{{}} = [\Phi _{{i + 1,i + 2}}^{{}}] \cdot R_{{i + 1}}^{c}, \\ \end{gathered} $
где необходимые матрицы перехода между системами координат определяются известными методами [1315].

Подстановки (3) в (2) дают окончательное выражение полного линейного ускорения i-й точки i-го тела. Выражение ускорения i-й точки (i + 1)-го тела получается аналогично. Тогда первое из уравнений связи (2) запишется следующим образом:

$\begin{gathered} m_{i}^{{ - 1}}({{F}_{{oi}}} + \left[ {{{M}_{{ci \to i}}}} \right] \cdot R_{i}^{c} - \left[ {{{M}_{{ci - 1 \to i}}}} \right] \cdot R_{{i - 1}}^{c}) + \\ + \;{{\omega }_{i}} \times ({{\omega }_{i}} \times r_{{i,i + 1}}^{{}}) + {{[{{I}_{{oi}}}]}^{{ - 1}}}\left\{ { - {{\omega }_{i}} \times \left( {[{{I}_{{oi}}}]{{\omega }_{i}}} \right) + } \right. \\ + \;{{L}_{{oi}}} - [\Phi _{{i,i + 1}}^{{}}] \cdot R_{i}^{c} + [\Phi _{{i,i - 1}}^{{}}] \times \\ \times \;R_{{i - 1}}^{c} + \left[ {{{M}_{{ci \to i}}}} \right] \cdot L_{i}^{c} - \left[ {{{M}_{{ci - 1 \to i}}}} \right] \cdot L_{{i - 1}}^{c}\} \times r_{{i,i + 1}}^{{}} = \\ = \;m_{{i + 1}}^{{ - 1}}({{F}_{{oi + 1}}} + - \left[ {{{M}_{{ci \to i + 1}}}} \right] \cdot R_{i}^{c} + \\ + \;\left[ {{{M}_{{ci + 1 \to i + 1}}}} \right] \cdot R_{{i + 1}}^{c}) + {{\omega }_{{i + 1}}} \times ({{\omega }_{{i + 1}}} \times r_{{i + 1,i}}^{{}}) + \\ + \;{{\left[ {{{I}_{{oi + 1}}}} \right]}^{{ - 1}}}\{ - {{\omega }_{{i + 1}}} \times ([{{I}_{{oi + 1}}}]{{\omega }_{{i + 1}}}) + {{L}_{{oi + 1}}} + [\Phi _{{i + 1,i}}^{{}}] \times \\ \times \;R_{i}^{c} - [\Phi _{{i + 1,i + 2}}^{{}}] \cdot R_{{i + 1}}^{c} - \left[ {{{M}_{{ci \to i + 1}}}} \right] \cdot L_{i}^{c} + \\ + \left[ {{{M}_{{ci + 1 \to i + 1}}}} \right] \cdot L_{{i + 1}}^{c}\} \times r_{{i + 1,i}}^{{}}. \\ \end{gathered} $

Второе из уравнений связи (2) путем аналогичных подстановок примет вид

$\begin{gathered} {{\left[ {{{I}_{{oi}}}} \right]}^{{ - 1}}}\{ - {{\omega }_{i}} \times (\left[ {{{I}_{{oi}}}} \right]{{\omega }_{i}}) + {{L}_{{oi}}} - [\Phi _{{i,i + 1}}^{{}}] \cdot R_{i}^{c} + \\ + \;[\Phi _{{i,i - 1}}^{{}}] \cdot R_{{i - 1}}^{c} + \left[ {{{M}_{{ci \to i}}}} \right] \cdot L_{i}^{c} - \left[ {{{M}_{{ci - 1 \to i}}}} \right] \cdot L_{{i - 1}}^{c}\} - \\ - \;{{\left[ {{{I}_{{oi + 1}}}} \right]}^{{ - 1}}}\{ - {{\omega }_{{i + 1}}} \times \left( {\left[ {{{I}_{{oi + 1}}}} \right]{{\omega }_{{i + 1}}}} \right) + {{L}_{{oi + 1}}} + [\Phi _{{i + 1,i}}^{{}}] \times \\ \times \;R_{i}^{c} - [\Phi _{{i + 1,i + 2}}^{{}}] \cdot R_{{i + 1}}^{c} - \left[ {{{M}_{{ci \to i + 1}}}} \right] \cdot L_{i}^{c} + \\ + \;\left[ {{{M}_{{ci + 1 \to i + 1}}}} \right] \cdot L_{{i + 1}}^{c}\} + {{\omega }_{i}} \times {{\omega }_{{i + 1}}} = 0. \\ \end{gathered} $

В общем виде полученные уравнения связи можно записать следующим образом:

(4)
$\begin{gathered} \text{[}A_{{i - 1}}^{R}]{\text{ }}R_{{i - 1}}^{c} + [A_{i}^{R}]{\text{ }}R_{i}^{c} + [A_{{i + 1}}^{R}]R_{{i + 1}}^{c} + \\ + \;[A_{{i - 1}}^{L}]{\text{ }}L_{{i - 1}}^{c} + [A_{i}^{L}]{\text{ }}L_{i}^{c} + [A_{{i + 1}}^{L}]{\text{ }}L_{{i + 1}}^{c} = {{B}_{i}}, \\ [C_{{i - 1}}^{R}]{\text{ }}R_{{i - 1}}^{c} + [A_{i}^{R}]{\text{ }}R_{i}^{c} + [A_{{i + 1}}^{R}]{\text{ }}R_{{i + 1}}^{c} + \\ + \;[A_{{i - 1}}^{L}]{\text{ }}L_{{i - 1}}^{c} + [A_{i}^{L}]{\text{ }}L_{i}^{c} + [A_{{i + 1}}^{L}]{\text{ }}L_{{i + 1}}^{c} = {{D}_{i}}. \\ \end{gathered} $

Таким образом, определение реакций связи, как это видно из структуры соотношений (4), сводится к решению системы линейных уравнений относительно неизвестных компонент векторов реакций связи $R_{i}^{c}$ и реактивных моментов $L_{i}^{c}$, которые затем подставляются в (1).

Уравнения связи в виде (4) справедливы для всех связей в “к”-й панели (крыле) (рис. 3), за исключением первой (i = 1) и последней (i = Iк).

Для i = Iк  $R_{{k,i + 1}}^{c} = 0$$L_{{k,i + 1}}^{c} = 0$.

Для первого тела (корпуса) главный вектор и главный момент первой связи (со штангами) получается суммированием по числу крыльев, т.е. для k-го крыла при индексе связи i = 1 уравнения связи будут иметь вид

$\begin{gathered} {\text{ }}\sum\limits_{k = 1}^k {[A_{{k,1}}^{R}]R_{{k,1}}^{c}} + [A_{{k,2}}^{R}]R_{{k,2}}^{c} + \\ + \;\sum\limits_{k = 1}^k {[A_{{k,1}}^{L}]L_{{k,1}}^{c}} + [A_{{k,2}}^{L}]L_{{k,2}}^{c} = B_{{k,1}}^{{}}, \\ \sum\limits_{k = 1}^k {[C_{{k,1}}^{R}]R_{{k,1}}^{c}} + [C_{{k,2}}^{R}]R_{{k,2}}^{c} + \\ + \;\sum\limits_{k = 1}^k {[C_{{k,1}}^{L}]L_{{k,1}}^{c}} + [C_{{k,2}}^{L}]L_{{k,2}}^{c} = D_{{k,1}}^{{}} \\ \end{gathered} $

В качестве иллюстрации на рис. 4 приведены эпюры изгибающего момента для аппарата, изображенного на рис. 2 (область положительных значений – 1, отрицательных – 2). Эпюры получены выборкой по времени максимальных значений в период процесса раскрытия. Они позволяют определить моментную нагрузку на вал привода раскрытия аппарата (значения, соответствующие нулевой координате), а также пиковые значения момента в “крыльях” солнечной батареи.

Рис. 4.

Огибающая эпюра изгибающего момента.

Отметим особенности моделирования многофазных схем раскрытия солнечных батарей. В процессе раскрытия общее число степеней свободы системы меняется (увеличивается), так как для каждой новой фазы характерно снятие дополнительных связей.

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

ВЫВОДЫ

Предложенный подход имеет следующие преимущества:

1. Универсальность подхода. Уравнения движения всех тел (свободных или нет) записаны одинаково. Отличие состоит в структуре главных векторов сил и моментов реакций связей, определяемой конкретным типом связи.

2. Подход дает возможность гибко формировать расчетную схему любой пространственной конфигурации. Подход легко поддается формализации и алгоритмизации.

3. Часто одной из основных целей расчета является определение реакций связей. В рамках данного подхода определение реакций связей – неотъемлемая процедура.

4. Определяется как абсолютное, так и относительное движение любого элемента рассматриваемой механической системы.

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

  1. Липанов А.М., Алиев А.В. Проектирование РДТТ. М.: Машиностроение, 1995. 397 с.

  2. Бакулин В.Н., Острик А.В. Комплексное действие излучений и частиц на тонкостенные конструкции с гетерогенными покрытиями. М.: Физматлит, 2015. 280 с.

  3. Легостаев В.П., Марков А.В., Сорокин И.В. Целевое использование Российского сегмента Международной космической станции: значимые научные результаты и перспективы // Космическая техника и технологии. 2013. № 2. С. 3–18.

  4. Рулев Д.Н., Рулев Н.Д. Планирование наблюдений астрономических объектов с космического аппарата с учетом ограничений на моменты выполнения наблюдений // Космическая техника и технологии. 2019. № 1 (24). С. 58–67.

  5. Израэль Ю.А., Лиознов Г.Л., Расновский А.А. Возможности космических и ядерных технологий в реформировании мировой энергетики XXI века // Известия РАН. Энергетика. 2008. № 3. С. 3–19.

  6. Коротеев А.С., Семенов В.Ф., Семенов Ю.П., Сизенцев Г.А. и др. Космическая техника и космонавтика в решении экологических проблем мировой энергетики XXI века // Известия РАН. Энергетика. 2006. № 1. С. 142–155.

  7. Бакулин В.Н., Борзых С.В., Щиблев Ю.Н. Перспективные схемы раскрытия трансформируемых космических конструкций и проблемы моделирования динамики процесса их раскрытия // Материалы XIII Междунар. конф. по Прикладной математике и механике в аэрокосмической отрасли (AMMAI'2020), посвященной 90-летию МАИ и 100-летию со дня рождения академика И.Ф. Образцова, сентябрь 2020. Алушта. М.: Изд-во МАИ, 2020. С. 405–407.

  8. Зимин В.Н., Бей Н.А. Трансформируемые антенны больших размеров для геостационарных космических аппаратов // Антенны. 2005. № 10. С. 24–27.

  9. Хамиц И.И., Филиппов И.М., Бурылов Л.С., Медведев Н.Г. Трансформируемые крупногабаритные конструкции для перспективных пилотируемых комплексов // Космическая техника и технологии. 2016. № 2 (13). С. 23–33.

  10. Bakulin V.N., Borzykh S.V. Modeling of the Deployment Process Dynamics for Large-Sized Transformable Space Structures. Russian Aeronautics. 2020. V. 63. P. 610–617.

  11. Бакулин Д.В. и др. Моделирование процесса раскрытия солнечных батарей // Математическое моделирование. 2004. № 6. С. 88–92.

  12. Бакулин В.Н., Борзых С.В., Щиблев Ю.Н., Ильясо-ва И.Р. Система уравнений минимальной размерности для описания процесса раскрытия солнечных батарей // Краевые задачи и математическое моделирование. Труды 10-й Всероссийской конференции. Новокузнецк: Изд-во НФИ КемГУ, 2010. Т. 2. С. 138–142.

  13. Ильясова И.Р. Динамика процесса раскрытия многозвенных солнечных батарей //Вестник Самарского гос. аэрокосмического ун-та им. С.П. Королева. 2012. Т. 35. № 4. С. 88−93.

  14. Бакулин В.Н., Борзых С.В., Ильясова И.Р. Математическое моделирование процесса раскрытия многозвенных солнечных батарей // Вестник МАИ. 2011. № 3. С. 295–302.

  15. Bakulin V.N., Bogomolov N.V., Borzykh S.V. Separation Algorithm of the Soyuz MS Spacecraft from Nonstabilized International Space Station. Russian Aeronautics. 2019. V. 62. № 4. P. 577–584.

Дополнительные материалы отсутствуют.

Инструменты

Доклады Российской академии наук. Физика, технические науки