Разработка системы ориентации и стабилизации малых космических аппаратов — различия между версиями

Материал из Department of Theoretical and Applied Mechanics
Перейти к: навигация, поиск
(Экпериментальные катушки)
м (Эксперименты)
 
(не показано 14 промежуточных версий 2 участников)
Строка 1: Строка 1:
==Руководитель==
+
'''Работу выполнил: [[Буй Ван Шань]]'''
А.С.Мурачев
+
 
 +
'''Научный руководитель: [[Андрей Мурачёв]]'''
 +
 
 +
Загрузить работу: [[Media:Диссертация Шань 2016.pdf|1.25МB]]
 +
 
 +
Загрузить презентацию: [[Media:Шань - кубсат-2016.pdf|1.27МB]]
 +
 
 
==Введение==
 
==Введение==
 
[[File:Picture1.png|200px|thumb|left|Рис.1.CubeSat 1U]]
 
[[File:Picture1.png|200px|thumb|left|Рис.1.CubeSat 1U]]
Строка 48: Строка 54:
 
*Напряженность имитатора приблизительно можно выразить следующей формулой
 
*Напряженность имитатора приблизительно можно выразить следующей формулой
 
[[File:HelmholzEquation.png|200px|thumb|center]]
 
[[File:HelmholzEquation.png|200px|thumb|center]]
 +
[https://www.comsol.com/model/download/313721/models.acdc.helmholtz_coil.pdf| Инструкция по моделированию катушек в Comsol]
  
 
===Экспериментальный образец имитатора===
 
===Экспериментальный образец имитатора===
Строка 88: Строка 95:
  
 
==Эксперименты==
 
==Эксперименты==
[[File:Motion.png|530px|thumb|left|Рис.5. ψ_желаемый=п/3]]
+
[[File:Motion.png|560px|thumb|left|Рис.5. ψ_желаемый=п/3]]
[[File:Motion1.png|530px|thumb|right|Рис.6. ψ_желаемый=2п/3]]
+
[[File:Motion1.png|560px|thumb|left|Рис.6. ψ_желаемый=2п/3]]
 
На рис.5 и рис.6 показаны графики реальных угловых положений спутника для разных заданных углов ψ_желаемый=п/3 и ψ_желаемый=2п/3 соответственно. Абсолютная ошибка ориентации спутника составляет 10 град. (~0.17 рад).
 
На рис.5 и рис.6 показаны графики реальных угловых положений спутника для разных заданных углов ψ_желаемый=п/3 и ψ_желаемый=2п/3 соответственно. Абсолютная ошибка ориентации спутника составляет 10 град. (~0.17 рад).
 +
 +
Файл для Arduino [[Media:CoilManagment.rar]]
 +
 +
Программа моделирования C# [[Media:Helmholtz_coil.rar ]]
  
 
==Выводы==
 
==Выводы==

Текущая версия на 00:46, 2 июля 2016

Работу выполнил: Буй Ван Шань

Научный руководитель: Андрей Мурачёв

Загрузить работу: 1.25МB

Загрузить презентацию: 1.27МB

Введение[править]

Рис.1.CubeSat 1U
  • Одной из самых важных задач при разработке космических аппаратов является создание системы

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

Задачи[править]

  • Создание магнитной системы управления угловым движением КА.
  • Разработка алгоритма управления при стабилизации КА по МПЗ.
  • Создание имитатора поля.

Магнитное управление, принцип работы и особенности[править]

Момент взаимодействия магнитной системы К.А с магнитным полем Земли (МПЗ):[править]

Equation1.png

Где: L – вектор магнитного момента К.А.; B – вектор индукции МПЗ.

Особенности магнитной системы управления[править]

  • Нельзя создать управляющий момент вдоль направления поля. Все возможные положения М заключены в плоскости, перпендикулярной В.
  • Управление по осям оказывается зависимым.
  • При совпадении вектора магнитного момента L и поля B управление невозможно

Принцип работы магнитной системы управления[править]

  • По принципу действия магнитная система демпфирования с катушечными и электромагнитными МИО – это по существу электродвигатель постоянного тока. В нем можно выделить все основные функциональные элементы: статор – это Земля, поле статора – МПЗ, ротор – КА, обмотки ротора – катушки МИО, коллектор – электронная часть системы, регулирующая и переключающая токи в МИЮ.
  • Отличаются эти системы только режимом работы и «конструктивными» особенностями. Так, в отличие от классического электродвигателя рассматриваемый «двигатель» работает большей частью в режиме торможения; обычный режим двигателя бывает только, когда есть необходимость в увеличении скорости вращения КА. «Конструктивная» особенность состоит в том, что здесь используемся в некотором смысле обращенная конструкция: ротор – КА с МИО – находится не внутри статора, а вне его (Земли). Кроме того, «ротор» не фиксирован относительно «статора»: он совершает линейное движение относительно «статора», а также имеет три степени свободны и соответственно три системы обмоток (катушек МИО) – это как бы двигатель со сферическим ротором.
  • Понятно, что энергетические соотношения, характерные для таких магнитных систем, будут такими же, как и для электродвигателей постоянного тока.

Создание имитатора поля[править]

Напряженность поля[править]

  • Выполним оценку максимальной требуемой напряженности магнитного поля в данном случае. Для простоты определим параметры поля, чтобы механические моменты, действующие на макет, могли бы вынудить его совершать гармонические колебания.
  • Требуемая индукция поля определяется следущим образом:
Equation2.png

где Т - период колебаний, J - момент инерции макета, A – амплитуда колебания.

  • Подставляя данные, получаем, что необходимая напряженность поля составляет H≈200 А/м. Если учитывать поправку на геомагнитное поле, то H надо увеличить еще на 50-100 А/м, чтобы надежно подавить геомагнитный фон. Таким образом получаем, что необходимая напряженность поля равна 250 А/м (~5 раз больше МПЗ).

Конструкция имитатора[править]

  • Существуют различные геометрические конфигурации токовых катушек для создания однородное поле. Наиболее распространенней конфигурацей имитатора поля является "кольца Гельмгольца".Обозначим радиус катушки через а_с, расстояние между центрами катушек - d_c (рис.4). Показано [1][2], что оптимальной конфигурацией является та, для которой выполняется соотношение d_c⁄(a_c=1.116).
Рис.2.Кольца Гельмгольца
  • Напряженность имитатора приблизительно можно выразить следующей формулой
HelmholzEquation.png

Инструкция по моделированию катушек в Comsol

Экспериментальный образец имитатора[править]

Рис.3.Имитатор поля

Создание магнитной системы управления[править]

Параметры токовых катушек[править]

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

  • материалом катушки - медь,
  • масса катушки M = 0.02 Кг,
  • диаметр витка l = 0.1 м,
  • площадь сечения материала обмотки должна быть равна S_r=0.0234мм^2, что соответствует диаметру d = 0.172 мм,
  • напряжение на катушке равно 12 В,
  • количество витков берем равным 240,
  • длина катушки (в случае однослойной намотки) ,
  • мощность каждой катушки W = 2 Вт.

Экпериментальные катушки[править]

Рис.4.Макет спутника с токовыми катушками
  • В соответствии с выбранными параметрами были созданы токовые катушки для спутника. Проволока, используемая для намотки катушек, была выбрана с диаметром 0.2 мм. Сопротивление каждой катушки R=N∙r=53.5 Ом. Измеряемое сопротивление каждой катушки составляет ≈55 Ом. То есть отклонение измеренного сопротивления от расчетного составляет ≈3 %. Такая разница может вызвана превышением числа витков в процессе ручной намотки. Используя формулы (20), получим создаваемый дипольный момент равен 0.42 Ам2. Относительное отклонение от требуемого значения момента составляет σ=5%. Разницы моментов и сопротивления катушек будут усчитаны при выводе закона управления системы. *Расположение двух катушек на макете спутника показано на рисунке справа.

Уравнение и закон управления движением спутника[править]

Уравнение движения спутника[править]

  • Уравнение движения спутника относительно вертикальной оси запишется в следующем виде:
MotionEquation.png

где: m – дипольный момент КА; ω – угловая скорость макета вокруг вертикальной оси; ψ – угол отклонения дипольного момента КА от вектора индукции МПЗ. J - момент инерции спутника относительно вертикальной оси

  • Если спутник содержит две токовых катушки, его дипольный момент складывается из моментов двух катушек. Перепишем уравнение движения:
MotionEquation2.png

где: m1,m2 – дипольный момент первой и второй катушек соответственно; ψ – угол между моментом первой катушки и индукция МПЗ.

Закон управления[править]

  • Желаемый дипольный момент катушек определяется следующим образом:

m=Kp*e+Kd*de/dt + Ki∫de где: - Kp, Kd, Ki: пропорциональный, дифференцирующий интегральный коэффициенты, соответственно; - e="ψ"_желаемый-"ψ" : разность между желаемым и измеряемым углами.

  • Для получения данного дипольного момента, необходимо приложить дипольные моменты вдоль двум остальным осям в соответствии с выражением:
Moment.png

где b = B/|B| - нормированный вектор индукции поля в связанной системе координат.

Эксперименты[править]

Рис.5. ψ_желаемый=п/3
Рис.6. ψ_желаемый=2п/3

На рис.5 и рис.6 показаны графики реальных угловых положений спутника для разных заданных углов ψ_желаемый=п/3 и ψ_желаемый=2п/3 соответственно. Абсолютная ошибка ориентации спутника составляет 10 град. (~0.17 рад).

Файл для Arduino Media:CoilManagment.rar

Программа моделирования C# Media:Helmholtz_coil.rar

Выводы[править]

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

Литература[править]

  • 1. M.Pastena and M.Grassi, Optimum Design of a three-axis Magnetic Field Simulator, IEEE Transaction on Aerospace and Electronic Systems, Vol.38, No.2, April 2002, pp.488-501.
  • 2. Коваленко А.П. Магнитные системы управления космическими летательными аппаратами. М., «Машиностроение», 1975, 248с.
  • 3. Д.В.Сивухин. Общий курс физики. Учебное пособие: для вузов. В 5 т. Т.III. Электричество. 4-е изд., стереот. М.: Физматлит; Издательство МФТИ, 2002 г, 656 с.
  • 4. Айзерман М.Ф. Классическая механика. «Наука» Изд.; М, 1980г, 368 с.
  • 5. Овчинников М.Ю., Цветков Е.А. Москва, 2005 Проектирование имитатора геомагнитного поля в составе лабораторного стенда для отработки способов управления ориентацией микроспутников.
  • 6. П.Л.Калантаров, Л.А.Цейтлин. Расчет индуктивностей: Справочная книга. 3-е изд., переработанное и дополненное. Л.: Энергоатомиздат, Ленинградское отделение, 1986 г, 488 с.
  • 7. Active magnetic control system for gravity gradient stabilized spacecraft. Francois Martel, Parimal K.Pal and Mark Psiaki, Cornell University
  • 8. Control systems, robotics and automation - Vol. II - PID Control - Araki M.
  • 9. Севастьянов Н.Н., Бранец В.Н., Панченко В.А., Казинский Н.В., Кондрации Т.В., Негодяцев С.С. Анализ современных возможностей создания малых космических аппаратов для дистанционного зондирования Земли. Труды МФТИ No3, том 1. -2009.
  • 10. Гущин В.Н. Основы устройства космических аппаратов. ¬ М.: Машиностроение, 2003.
  • 11. Ziegler, J.G & Nichols, N. B. (1942). "Optimum settings for automatic controllers". Transactions of the ASME 64.
  • 12. Euston, M.; Coote, P.; Mahony, R.; Kim, J.; Hamel, T. A complementary filter for attitude estimation of a fixed-wing UAV. In Proceedings of IEEE/RSJ International Conference on Intelligent Robots and Systems, Acropolis Convention Center, Nice, France, 22–26 September 2008.