Результат интеллектуальной деятельности: СПОСОБ ОДНООСНОЙ ОРИЕНТАЦИИ КОСМИЧЕСКОГО АППАРАТА ВЫТЯНУТОЙ ФОРМЫ

Изобретение относится к космической технике и может быть использовано для ориентации космического аппарата (КА) при выполнение экспериментов и исследований.

Известен способ ориентации КА, включающий выставку осей КА и поддержание углового положения КА с помощью двигателей ориентации (Алексеев К.Б., Бебенин Г.Г. «Управление космическими летательными аппаратами», М.: Машиностроение, 1974).

Однако для использования данного способа необходимо расходовать рабочее тело, что приводит, кроме того, к загрязнению оптических поверхностей КА и вызывает микроускорения на борту КА.

Известен способ, включающий выставку оси КА, соответствующую минимальному моменту инерции, на центр Земли и орбитальное смещение КА (Беляев М.Ю. «Научные эксперименты на космических кораблях и орбитальных станциях», М.: Машиностроение, 1984). Данный способ используется для КА, имеющих вытянутую форму, т.е. когда момент инерции относительно продольной оси значительно меньше момента инерции относительно поперечных осей.

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

Однако вследствие неточной выставки оси КА на центр Земли появляются угловые скорости вокруг всех осей КА. Наличие угловых скоростей вокруг поперечных осей КА приводит к отклонению продольной оси аппарата от направления к центру Земли, вследствие чего ухудшается точность гравитационной ориентации КА.

Наиболее близким к предлагаемому является способ одноосной ориентации КА вытянутой формы (Патент РФ №2457159, приоритет от 30.08.2010, МПК (2006.01) B64G 1/34 - прототип), включающий выставку оси КА, соответствующей минимальному моменту инерции, на центр Земли и орбитальное смещение КА, при этом после выставки оси КА на центр Земли и орбитального смещения КА производят закрутку КА вокруг выставленной на центр Земли оси КА до требуемого момента с угловой скоростью , где I_yz - среднее значение близких по величине моментов инерции КА вокруг поперечных осей КА; I_x - момент инерции КА вокруг продольной оси; ω₀ - модуль абсолютной угловой скорости орбитальной системы координат.

Способ-прототип позволяет повысить точность одноосной ориентации конкретно рассмотренного типа КА и тем самым снизить также микроперегрузки на КА, возникающие при раскачке и переходе КА в режим неуправляемого вращения. В общем случае вращение КА с указанной скоростью не является устойчивым для всех типов КА вытянутой формы - в общем случае со временем отклонение продольной оси КА от направления к центру Земли становится все более существенным, что приводит к «кувырканию» КА и разрушению гравитационной ориентации. Это в том числе ограничивает возможности проведения экспериментов, требующих наведения научной аппаратуры на Землю и/или низкого уровня микроускорений.

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

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

Технический результат достигается тем, что в способе одноосной ориентации КА вытянутой формы, включающем развороты КА и закрутку вокруг его продольной оси, соответствующей минимальному моменту инерции, дополнительно перед выполнением закрутки космический аппарат разворачивают до совмещения продольной оси КА с плоскостью, образованной нормалью к плоскости орбиты и радиус вектором КА, и достижения углом между продольной осью КА и плоскостью орбиты максимального допустимого значения угла отклонения продольной оси КА от местной вертикали при значении угла между радиус-вектором КА и вектором, направленным из центра масс КА в центр аэродинамического давления солнечных батарей КА, <90°, после чего выполняют закрутку КА вокруг продольной оси, при этом величину угловой скорости закрутки определяют по формуле , где α≤1,

λ - отношение минимального главного центрального момента инерции к среднему значению поперечных главных центральных моментов инерции КА,

ω₀ - угловая скорость орбитального движения КА,

β₀ - максимальное допустимое значение угла отклонения продольной оси КА от местной вертикали.

Поясним предложенные в способе действия.

Запишем уравнения вращательного движения КА. КА считается твердым телом, геоцентрическое движение его центра масс - кеплеровым эллиптическим. Элементы этого движения находятся по данным радиоконтроля орбиты. Для записи уравнений введем две правые декартовы системы координат - орбитальную ОХ₁Х₂Х₃ и образованную главными центральными осями инерции КА Ox₁x₂x₃. Точка О - центр масс КА, оси ОХ₃ и OX₁ направлены, соответственно, вдоль геоцентрического радиуса-вектора точки О и по трансверсали к орбите в этой точке. Ось Ox₁ направлена вдоль продольной оси КА.

Положение системы Ox₁x₂x₃ относительно системы OX₁X₂X₃ будем задавать углами γ, δ и β, которые введем следующим образом. Система OX₁X₂X₃ может быть переведена в систему Ox₁x₂x₃ тремя последовательными поворотами: 1) на угол δ+π/2 вокруг оси ОХ₂, 2) на угол β вокруг новой оси ОХ₃, 3) на угол γ вокруг новой оси ОХ₁, совпадающей с осью Ох₁. Матрицу перехода от системы Ox₁x₂x₃ к системе ОХ₁Х₂Х₃ обозначим , =1, где a _i - косинус угла между осями OX_i и Ox_j. Элементы этой матрицы выражаются через введенные углы с помощью формул

а ₁₁=-sinδcosβ,

а ₁₂=cosδsinγ+sinδsinβcosγ,

a ₁₃=cosδcosγ-sinδsinβsinγ,

a ₂₁=sinβ, a ₂₂=cosβcosγ, a ₂₃=-cosβsinγ,

a ₃₁=-cosδcosβ,

a ₃₂=-sinδsinγ+cosδsinβcosγ,

a ₃₃=-sinδcosγ-cosδsinβsinγ.

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

,

,

,

,

Здесь ω_i (i=1, 2, 3) - компоненты абсолютной угловой скорости КА в системе координат Ox₁x₂x₃, ω₀ - угловая скорость орбитального движения КА, I_i -моменты инерции КА относительно осей Ox_i (главные центральные моменты инерции КА), М_ℓ - вычисленный относительно точки О восстанавливающий аэродинамический момент, приложенный к КА.

Уравнения (1) позволяют оценить вращательные движения КА при различных условиях (Белецкий В.В. Движение искусственного спутника относительно центра масс. М., Наука, 1965; Белецкий В.В. Движение спутника относительно центра масс в гравитационном поле. М.: Издательство МГУ, 1975; Черноусько Ф.Л. Об устойчивости регулярной прецессии спутника. Прикладная математика и механика, 1963, т. 28, вып. 1, с. 155-157).

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

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

Здесь λ - отношение минимального главного центрального момента инерции (момент инерции КА вокруг продольной оси I_x) к среднему значению близких по величине поперечных главных центральных моментов инерции (среднее значение моментов инерции КА вокруг поперечных осей КА I_yz), ω₁ - абсолютная угловая скорость закрутки КА вокруг продольной оси.

В случае λ<<1 даже при сравнительно большом отношении |ω₁/ω₀| угол β₀ будет мал. Продольная ось спутника с неравными моментами инерции, но удовлетворяющего перечисленным выше условиям, будет совершать движение в окрестности номинального положения малые колебания.

В этом случае номинальными невозмущенными движением спутника следует считать его движения, принадлежащие двумерному интегральному многообразию уравнений движения. Это многообразие строится в виде формальных рядов по целым степеням малых параметров, характеризующих возмущения. В частности, возмущения, вызванные отличием спутника от осесимметричного, характеризуются параметром µ=(I₂-I₃)/I₁, где I₂≈I₃; возмущения, вызванные эллиптичностью орбиты, характеризуются ее эксцентриситетом и т.п. Указанное интегральное многообразие параметризуется двумя перемененными. Одна из них близка угловой скорости ω₁, вторая близка углу поворота спутника вокруг продольной оси. При малых параметрах, равных нулю, интегральное многообразие переходит в номинальный режим осесимметричного спутника, а его параметры в точности совпадают с ω₁ и указанным углом.

Формальные ряды, представляющие интегральное многообразие, можно построить не при всех значениях ω₁. Чтобы построение было возможно, величины ω₀ и ω₁ должны быть нерезонансными. Имеется в виду отсутствие резонансов между колебаниями продольной оси спутника в окрестности указанного выше номинального положения и вращением спутника вокруг этой оси. Вблизи резонанса амплитуда колебаний продольной оси спутника в окрестности номинального положения быстро возрастает, и режим разрушается.

Можно показать (Сарычев В.А. Вопросы ориентации искусственных спутников. - Итоги науки техники. Серия Исследования космического пространства. Т. 11. М.: ВИНИТИ, 1978; Хейл Дж. Колебания в нелинейных системах. М.: Мир, 1966; Черноусько Ф.Л. Об устойчивости регулярной прецессии спутника. Прикладная математика и механика, 1963, т. 27, №3, с. 474; Боголюбов Н.Н., Митропольский Ю.А. Асимптотические методы и теория нелинейных колебаний. М.: Физматгиз, 1963; Моисеев Н.Н. Асимптотические методы нелинейной механики. М.: Наука, 1969; Хайнбокел Дж., Страбл Р.А. Периодические решения систем дифференциальных уравнений, обладающих симметрией. - Механика. 1966. №1, с. 3; Пуанкаре А. Новые методы небесной механики. Избранные труды. Т. 1 и 2. М.: Наука, 1971-1972), что при λ<<1 «резонансные» значения угловой скорости вращения КА определяемых равенствами:

.

Здесь величина Ω - скорость изменения угла поворота КА вокруг продольной оси, отсчитываемого от плоскости орбиты.

С учетом (3) условие отсутствия резонансов выражается неравенствами

или

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

Условия (4) означают, что средняя угловая скорость вращения КА вокруг продольной оси не должна находиться в резонансе с собственными частотами колебаний этой оси по углам δ и β - углам, задающим положение продольной оси КА в орбитальной системе координат (β - угол между продольной осью КА и плоскостью орбиты, δ - угол между проекцией продольной оси КА на плоскость орбиты и направлением на Землю). Как указывалось выше такие частоты v определяются уравнением . При колебания продольной оси КА в плоскости орбиты и относительно этой плоскости с точностью до членов порядка O(β₀) независимы и соответствующие им частоты приближенно равны и . Тогда условия (4) принимают вид

При нарушении одного из неравенств первой группы (6) имеет место резонанс между вращением КА вокруг продольной оси и колебаниями этой оси в плоскости орбиты (δ-резонанс), при нарушении одного из неравенств второй группы (6) имеет место резонанс между вращением КА вокруг продольной оси и колебаниями этой оси в направлении, перпендикулярном плоскости орбиты (β-резонанс). С учетом и (6) соотношения для соответствующих «резонансных» значений угловой скорости вращения спутника ω₁ записываются в виде

Например, для КА с характеристикой λ≈0.15 (например, транспортный грузовой корабль (ТГК) «Прогресс») величины ν_δ/ω₀, ν_β/ω₀ составляют значения соответственно ≈1.6 и ≈1.88, а резонанс первого порядка (n=1) реализуется при значениях ω_δ,1/ω₀ и ω_β,1/ω₀ соответственно ≈0.83 и ≈0.98.

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

Рассматриваем ориентацию КА, в которой центр аэродинамического давления СБ КА расположен вне отрезка между центрами масс КА и Земли - в данной ориентации значение угла между радиус-вектором КА и вектором, направленным из центра масс КА в центр аэродинамического давления СБ КА, <90°. Тогда в процессе упомянутой выше эволюции убывает и КА со временем попадет в один из резонансов, в результате чего режим ориентированного движения быстро разрушается.

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

где β₀ - максимальное допустимое значение угла отклонения продольной оси КА от местной вертикали; а задаваемый коэффициент α≤1 определяет допустимое уменьшение угловой скорости закрутки КА относительно указанного максимально допустимого значения. Величина данного коэффициента определяется, в том числе точностями построения требуемой ориентации и отработки импульса для закрутки КА, реализуемых системой управления КА, и требованием к длительности поддержания гравитационной ориентации.

В качестве примера приведем результаты анализа возможности реализации длительного (сотни витков) поддержания гравитационной ориентации в случае, когда центр аэродинамического давления СБ КА расположен вне отрезка между центрами масс КА и Земли, для КА с характеристиками λ≈0.15, ω₀≈0.065°/сек при высоте орбиты порядка 400 км и длительности витка порядка 1,5 часа. Например, при требовании отклонения продольной оси КА от местной вертикали не более величины 15÷10° в (6) можно принять α≈0.6÷1, что соответствует начальному значению угловой скорости вращения ΚΑ ω₁≈0.25÷0.39°/сек, при которой начальное отклонение продольной оси КА от местной вертикали составляет величину β₀≈15÷40°, при этом численное моделирование показывает, что длительность поддержания гравитационной ориентации до момента наступления резонанса (т.е. до нарушения заданного требования к точности поддержания гравитационной ориентации - разрушения гравитационной ориентации) составляет не менее, чем сотни витков.

Опишем технический эффект предлагаемого изобретения.

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

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

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

В настоящее время технически все готово для реализации предложенного способа таком КА, как ТГК «Прогресс». Для реализации разворотов, закрутки и вычислений могут использоваться штатные средства системы управления ТГК «Прогресс» - система управления движением и навигацией, включая систему автономной навигации, солнечные датчики, датчики угловой скорости, двигатели ориентации, бортовой вычислитель и т.д.