ОЦЕНКА ВОЗМУЩАЮЩЕГО МОМЕНТА ПРИ МОДЕЛИРОВАНИИ ДВИЖЕНИЯ СИСТЕМЫ НАВЕДЕНИЯ АНТЕНН

Полный текст

Этим обусловлен переход на использование крупногабаритных антенн и высоких связных частот, что приводит к чрезвычайно жестким требованиями к управлению наведением. Ввиду жесткого лимита веса и объема космического аппарата созданные большие конструкции антенн являются гибкими в силу выбора легкого материала, а упаковка такой конструкции получается более плотной. Также большие космические конструкции имеют большое количество форм колебаний, которые обладают низкой частотой и близко расположены друг к другу. Современные космические средства связи предъявляют всё более и более жесткие требования к функциональным возможностям предоставляемых услуг связи, её непрерывности и готовности, независимо от местонахождения потребителя и времени суток. Данные функции реализуемы лишь при наличии возможностей оперативного и точного наведения антенн на наземные пункты, международную космическую станцию, космические аппараты, находящиеся на низкой и круговой орбите. 22 Математика, механика, информатика В результате антенная система становится более сложным многомерным объектом управления, и это диктует новые требования к системам управления положением антенн. Необходимость ограничения возмущающего момента на аппарат при соблюдении заявленной точности и скорости сопровождения объекта создала задачу построения системы наведения антенн, не влияющей на ориентацию и стабилизацию космического аппарата при выполнении основной функции. Конструкция антенны, для управления которой проектировалась система наведения антенн, выполнена из 16 углепластиковых спиц, на которые натянуто радиоотражающее сетеполотно. Антенна в рабочем положении (рис. 1) установлена на штангу. Рис. 1. Общий вид антенны в рабочем положении На этапе создания системы наведения антенн было уделено внимание оценке возмущающего момента при математическом моделировании работы системы при управлении положением нежесткой антенной. Моделирование даёт возможность увидеть, как будет вести себя система в действительности. Для получения реальной картины поведения системы была разработана математическая модель, которая учитывает воздействие факторов, оказывающих наибольшее влияние на работу системы. Электромеханический модуль (ЭММ) блока механического системы наведения антенн (БМ СНА) моделируемой системы наведения антенны состоит из двух параллельных ветвей по пять ступеней, последняя ступень является волновой передачей, связанной с выходным валом. В данной модели представлена только одна из ветвей. По результатам анализа структуры ЭММ БМ СНА разработана расчётная эквивалентная схема модели. Данная схема представлена в виде кинематической схемы, с учётом нежёсткости и люфтов редуктора и упругости антенны, на рис. 2. Все ступени редуктора содержат нелинейность типа люфт. Математическая модель редуктора (рис. 2) с учётом нежёсткости и люфтов представляет собой систему дифференциальных уравнений (2)-(11), представленных ниже. В уравнениях обозначено: J1 - момент инерции вала синхронного двигателя с электромагнитной редукцией (СДЭР) и первой шестерни редуктора; AI1, AI2, AI3 - амплитуда тока в соответствующей фазе; CM1, CM2, CM3 - коэффициенты момента СДЭР для соответствующей фазы; 01 - угол поворота ротора СДЭР и первой шестерни редуктора; SIN(2-n-F-t), SIN(2-n-F-t+120°), SIN(2-n-F-t-120°) -синусоидальная форма тока в соответствующей фазе; SIN(2-n-0i-t), SIN(2-n-0i-t+120°), SIN(2-n-01-t-120°) -синусоидальная форма момента СДЭР в соответствующей фазе; COF1, COF2 - коэффициенты, учитывающие вторую и третью гармоники в форме момента СДЭР; PF1, PF2, PF3 - погрешность сдвига фаз для соответствующей фазы СДЭР; F - частота тока в фазах двигателя; KD - коэффициент вязкого трения (коэффициент демпфирования) в СДЭР; MTP1 - момент сухого трения на валу СДЭР; K - коэффициент электромагнитной редукции СДЭР; 02, v1, а, рь р2, и1, уь у2, р - углы поворота соответствующих шестерён в редукторе (рис. 2); д - угол поворота антенны; J2, J3, J4, J5, J6, J7, J8, J9 - моменты инерции соответствующих шестерён в редукторе (рис. 2); J10 - момент инерции антенны (рис. 2); R1, R2, R3, R4, R5, R6, R7, R8, R9, R10 - радиусы соответствующих шестерён в редукторе; C2, C3, C4, C6, C8 - жесткости на кручение соответствующих осей в редукторе (рис. 2); C1, C5, C7 - упругое сопротивление («твердость») зубцов (первой, четвертой и пятой ступеней редуктора) шестерён при соприкосновении зубцов после выборки люфта; C9 - жесткость на кручение выходного вала, в месте крепления антенны к БМ, для моделирования первого тона упругих колебаний антенны; L1, L2, L3 - люфты на первой, четвертой и пятой ступенях редуктора (рис. 2); KD2, KD4, KD6, KD8 - коэффициенты вязкого трения при демпфировании упругих колебаний в редукторе; KD9 - коэффициент вязкого трения при демпфировании упругих колебаний антенны; MTP2 - момент сухого трения на последней (выходной) ступени редуктора; DX, DY, DZ - относительное линейное перемещение зубьев шестерён соответствующих радиусов; DXL, DYL, DZL - «проникновение» одного зуба в другой после того, как выберется люфт DX = 01 -R1+ 02 -R2, (1) если DX < 0, то DXL= DX + L1, при этом, если DXL > 0, то DXL=0, если DX > 0, то DXL = DX - L1, при этом, если DXL < 0, то, DXL = 0 то есть люфт не выбран. 23 Вестник СибГАУ. № 3(49). 2013 r6,c4, ь, Rs J5 J8 R10 R7 С5 - RaCe, 91 С7 МТР2 С3,92 ' МТРЗ с,, m Рис. 2. Кинематическая расчётная схема модели ЭММ БМ СНА Математическая модель синхронного двигателя с электромагнитной редукцией в виде дифференциальных уравнений имеет следующий вид: J • 01 = AI1 • CM 1 • SIN(2-п-F•t)X х [ SIN(K • 01 + PF1) + COF1 • SIN(2 • K •01 + PF1) + + COF2 • SIN(3 • K •01 + PF1)] + AI2 • CM2 х x SIN(2 • п • F • t +120°) • [ SIN(K •01 +120° + PF2) + + COF1 • SIN(2 • K •01 +120° + PF2) + + COF2 • SIN(3 • K •01 +120° + PF2)] + AI3 • CM3 x (2) x SIN (2 •nF • t -120°) •[ SIN (K •01 -120° + PF 3) + + COF1 • SIN(2 • K •01 -120° + PF3) + COF2 • SIN x x (3 • K •01 -120° + PF 3) ]-KD ^ - - DXL • C1 • R1 - MTP1 • SIGN(0), J2 •02 =-(02 -C1)•C2-(02-u1)x x KD2 - DXL • C1 • R2, • • • • J3 •u1 = (02-u1) •C2-(02-u1) KD2 - - (N2•u1 +a)• N2• C3 -(N2•u1 + a)• N2• KD3, J4 •a = -(N2•u1 +a) • C3-(N2 •u1 + a) • KD3 - - (N3^a + P1) • N3 •C4 -(N3 ^p’^N3^KD4, N2 = R3/R4, N3 = R5/R6, DY = p1 -R7+ p2 -R8, J5 •p1 =-(N3•a + p1)C4 + (N3•a+P1)> x KD4 - DYL•C5• R7; J6 ^2 = -DYL • C5 • R8 - (6) (7) - (P2 -Y1)• C6-(P2-y>KD6; DZ = y1 -R9 + y2 R10, если DZ < 0, то DZL = DZ + L3, при этом, если DZL > 0, то DZL = 0, если DZ > 0, то DZL = DZ - L3, при этом, если DZL < 0, то, DZL = 0. (3) (4) (5) J7-Y1 = (в2-Y1) C 6 + (в2-Y1) x x KD6 - DZL •C7 • R9, J8 • Y 2 = (P-Y 2) C8 + (P-Y 2^ KD8 - - DZL •C7 • R10 - MTP2 • SIGN(y2), J9 • P = -(p-Y2)C8-(p-Y2)^KD8 + + (ц-p) C9 + (ц-p)^ KD9, •• • • Jw ^ = -(^-p) C9 - (ц-p) x x KD9 -MTP3 • SIGN(ц). (8) (9) (10) (11) В блоке управления системы наведения антенны реализована импульсная форма. При импульсной форме тока в фазах двигатель работает в шаговом режиме. Период тока T разделен на шесть частей, средняя величина в радианах определяется по формуле: если DY < 0, то DYL = DY + L2, при этом, если DYL > 0, то DYL = 0; если DY > 0, то DYL = DY - L2, при этом, если DYL < 0, то, DYL = 0. ^1Ш 2 •п 6K (12) Средняя угловая скорость вычисляется по формуле: • 2^ п F K (13) 24 Математика, механика, информатика Для построения амплитудно-частотных характеристик колебаний скорости и возмущающих моментов использовалась программная среда Matlab, функция fft [2, с. 138], реализующая алгоритм быстрого преобразования Фурье. Для оценки возмущающего момента, действующего на космический аппарат (в данном случае при работе СНА в режиме программного наведения), рассматривался диапазон частот 0,1-1 Гц. Для вариантов расчётов были получены частотные спектры возмущающего момента. По результатам анализа проведённых расчётов и полученных частотных спектров были определены максимальные амплитуды возмущающего момента и соответствующие им частоты, которые приведены в таблице. Возмущающий момент немного более 0,02 Н • м (200 гс см) присутствует в вариантах расчётов 7, 8, 10. Проведённая оценка по результатам моделирования работы системы показала: - величина полученных возмущающих моментов не превышает 0,02 Нм в диапазоне частот 0,1-1,0 Гц во всей временной области; - наблюдаемые возмущающие моменты не корел-лируются с частотами управления, а определяются параметрами зубчатых пар редуктора и выходной ступени. В результате проведенной работы была получена математическая модель системы наведения антенны, позволяющая моделировать динамические процессы, возникающие в ней, а так же производить оценку возможных возмущающих моментов, которые могут создаваться при её работе. Данные результаты могут быть использованы для оценки динамических процессов на космическом аппарате в целом и оценки влияния на его точность ориентации и ориентацию в пространстве антенны при наведении на наземные пункты, международную космическую станцию, космические аппараты, находящиеся на низкой и круговой орбите. Максимальные амплитуды возмущающего момента и соответствующие им частоты № варианта аЦУ,' а Ц У, c Частота 1, Гц Момент 1, Нхм Частота 2, Гц Момент 2, Нхм 1 50 0,2 0,409 0,0142 0,247 0,0104 2 50 0,5 0,1698 0,0184 0,3196 0,0135 3 50 -0,2 0,368 0,01 0,667 0,00856 4 50 -0,5 0,0499 0,0190 0,29 0,0143 5 150 0,2 0,4489 0,0133 0,92 0,0159 6 150 0,5 0,066 0,0191 0,332 0,0105 7 150 -0,2 0,227 0,0293 0,399 0,0194 8 150 -0,5 0,06 0,0223 0,238 0,0209 9 400 0 0,078 0,0104 0,306 0,0103 10 500 0,5 0,829 0,0388 0,889 0,0139 11 0 0,5 0,83 0,0213 0,76 0,0148 12 10 -0,5 0,648 0,0106 0,282 0,00856 Передаточное число примененного редуктора = R2 • R4 • Ж, R8 ' R10 = ^9. Р R1 R3 R5 R7 R9 С учетом параметров редуктора были определены собственные частоты колебаний ступеней редуктора с учётом нежесткости на кручение [1, с. 305]: ю = , (14) V J 102,6 Гц, 127,11 Гц, 103,4 Гц, 12,08 Гц, 20,4 Гц. На основе проведенного модального анализа антенны в рабочем положении на штанге было получено значение частоты первого тона собственных колебаний антенны 0,83 Гц. Максимально возможный момент инерции антенны составляет: J10 = 50 кг •м2. При проектировании системы были задано следующее ограничение по работе системы: во всем диапазоне скоростей не должен создаваться возмущающий момент более 0,02 Н м (200 гс см), в диапазоне частот 0,1-1 Гц амплитуда периодической составляющей возмущающего момента должна быть не более 0,001 Н •м (10 гс •см). Для оценки возмущающего момента, действующего на космический аппарат при работе системы наведения антенн, рассматривались несколько вариантов. Данные варианты приведены в таблице. Возмущающий момент определяется как произведение углового ускорения на момент инерции нагрузки, приведённый к выходному валу ЭММ БМ СНА. То есть оценка возмущающего момента при движении БМ СНА производилась относительно места закрепления на элементе конструкции спутника. По результатам проведенного моделирования движения антенны для оценки возмущающего момента от БМ СНА использовалась формула: мвм = jh . (15) dt 25 Вестник СибГАУ. № 3(49). 2013 Рис. 3.Частотный спектр возмущающего момента на выходе редуктора (вариант 8) Рис. 4. Частотный спектр возмущающего момента на выходе редуктора (вариант 10)

Об авторах

А. В. Гончарук

ОАО «Информационные спутниковые системы» имени академика М. Ф. Решетнёва»

Email: gav78@iss-reshetnev.ru
Россия, 662970, Красноярский край, г. Железногорск, ул. Ленина, 52

М. В. Лукьяненко

Сибирский государственный аэрокосмический университет имени академика М. Ф. Решетнёва

Email: info@sibsau.ru
Россия, 660014, г. Красноярск, просп. имени газеты «Красноярский рабочий», 31

В. А. Раевский

Сибирский государственный аэрокосмический университет имени академика М. Ф. Решетнёва

Email: info@sibsau.ru
Россия, 660014, г. Красноярск, просп. имени газеты «Красноярский рабочий», 31

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

Дополнительные файлы