ЭЛЕКТРОННЫЕ И ЭЛЕКТРОМАГНИТНЫЕ
УСТРОЙСТВА
УДК 681.5.11
DOI: 10.17586/0021-3454-2015-58-6-458-463
АНАЛИЗ ДИНАМИЧЕСКИХ СВОЙСТВ ТРЕХМАССОВЫХ ИСПОЛНИТЕЛЬНЫХ ОСЕЙ С ОДНИМ И ДВУМЯ ВРАЩАЮЩИМИ МОМЕНТАМИ
В. А. Толмачев, М. В. Никитина, М. Е. Сергеева
Университет ИТМО, 197101, Санкт-Петербург, Россия E-mail: [email protected]
Предложена векторно-матричная математическая модель, пригодная для анализа механических переходных процессов в однодвигательных и двухдвигатель-ных трехмассовых механизмах исполнительных осей, представляемых структурными схемами неразветвленного типа. Приведены расчетные соотношения для определения резонансных частот. Достоверность последних подтверждена результатами спектрального анализа кривых механических переходных процессов, проведенных с использованием пакета MatLab. Приводятся расчетные соотношения для определения параметров эквивалентных двухмассовых моделей двухдвигательных электроприводов с трехмассовыми механизмами осей, полученные на основе преобразований исходных моделей, методом электрических аналогий. Результаты исследования закладывают основу для выбора структурных решений и для создания строгой методики синтеза систем управления следящих электроприводов телескопов с двухдвигательными трехмассовыми механизмами осей.
Ключевые слова: исполнительная ось, трехмассовый механизм, математические модели, механические переходные процессы, резонансные частоты, спектральный анализ.
Введение. В настоящее время широко применяются оптические комплексы высокоточных наблюдений. Они позволяют решать такие задачи, как автономный поиск и обнаружение космических объектов, определение расстояния до космических аппаратов, измерение угловых координат объекта и сбор информации о них.
К основным задачам, возникающим при создании следящих электроприводов комплексов высокоточных оптических измерений, относится обеспечение:
— широкого диапазона (до 18 000) скоростей слежения при малых значениях среднеквадратичных ошибок (1—2");
— плавного движения оптической оси комплекса при инфранизких скоростях слежения (до единиц угловых секунд за секунду) в условиях нежесткости конструкции опорно-поворотного устройства (ОПУ) и возможного возникновения механического резонанса на частотах от единиц до нескольких десятков герц.
Для достижения указанных условий при построении систем управления применяются безредукторные электроприводы на основе бесконтактных высокомоментных вентильных электродвигателей и прецизионных цифровых датчиков координат с разрешающей способно-
стью по углу менее сотых долей угловой секунды [1, 2]. Помимо того, важную роль в решении указанных задач играет конструктивное исполнение всех элементов измерительного телескопа и в первую очередь — осей его ОПУ. Угловые частоты резонанса механизма оси, вызванного крутильными деформациями, в конечном итоге определяют полосу пропускания частот системы управления и влияют на ее динамическую точность.
В большинстве случаев ОПУ оптических комплексов являются двухосными и построены по альт-азимутальному принципу [3]. Известные безредукторные следящие электроприводы осей ОПУ являются однодвигательными. Конструкция азимутальных осей двухосных ОПУ, как правило, не позволяет использовать безредукторные двухдвигательные приводы, однако при построении угломестных осей это вполне возможно.
Выигрыш по массогабаритным показателям можно получить, используя двигатели меньшей мощности на обеих полуосях вилки телескопа [3]. Однако при такой конструкции необходимо обеспечивать высокие динамические качества системы слежения.
Определить целесообразность использования двух двигателей в следящем электроприводе с трехмассовой исполнительной осью можно, проанализировав процессы в механической подсистеме. Именно такой анализ и является предметом настоящей статьи.
Векторно-матричные математические модели трехмассового механизма с одним и двумя вращающими моментами. Структурная схема трехмассового механизма, представляемого моделью последовательного типа с двумя вращающими моментами М1 и М2, приведена на рис. 1 (здесь й1, 02, 03 — угловые скорости первой; второй и третьей масс, М12 и М23 — моменты упругих связей; М1 и М2 — вращающие моменты на полуосях; С12 и С23 — крутильные жесткости левого и правого валов оси; •1, •3 — моменты инерции левого и правого валов; J2 — момент инерции трубы телескопа Мс1 и Мс2 — моменты нагрузок на опорах оси).
м
м
1
•Л р
м
12 р
м
1 с 23 м 23
• 2 Р о 2 р
•з р
о
Рис. 1
Механические переходные процессы можно описать векторно-матричным уравнением
лу+ви, (1)
Т Т
где У = [01М12 02 М23 03] — вектор состояния, и = [М1 М2 Мс1 Мс2 Мс3] — вектор управления, а матрицы системы А и управления В имеют вид
А =
0 1 - • 0 0 0 1 • 0 1 - • 0 0
С12 0 -Си 0 0 0 0 0 0 0
0 1 ~2 0 1 - т2 0 , В = 0 0 0 1 - 0
0 0 С23 0 -С23 0 0 0 0 0
0 0 0 1 0 0 1 • 0 0 1 - •
Векторно-матричная математическая модель (1) позволяет рассчитать процессы в системе как с одним, так и с двумя вращающими моментами. При расчете с одним вращающим моментом М1 в исходных данных программы расчета достаточно положить М2=0.
1
460 В. А. Толмачев, М. В. Никитина, М. Е. Сергеева
Аналитический расчет резонансных частот механизма. В работе [4] предложена методика синтеза системы управления электропривода с четырехконтурной структурой и двух-массовой моделью механизма оси. Методика синтеза предполагает определение частоты механического резонанса и коэффициента соотношения масс, устанавливающих предельную полосу пропускания скоростной подсистемы электропривода. В работе [5] для случая многомассового механизма оси предложено при синтезе скоростного контура системы подчиненного регулирования угла поворота ориентироваться на низшую резонансную частоту, которую, как правило, определяют экспериментально.
Для рассматриваемого случая трехмассового механизма (см. рис. 1) резонансные частоты могут быть получены расчетным путем. Следует заметить, что с помощью аналитических соотношений для расчета резонансных частот, приведенных в известном учебнике [6], не представляется возможным провести спектральный анализ.
Используя матрицу А (1), получим характеристическое уравнение ёе1;(р1-А)=0:
р (р4 + ар2 + Ь ) = 0, (2)
где а = С12• з(• + •2) + С23■/,(J2 + .3) , Ь = С12С2з(•/, + •2 + •3) .
Решив биквадратное уравнение (2), получим корни характеристического уравнения в виде р1=0, р2,3=±/юр1, р4,5=±/юр2, где юр1 и юр2 — резонансные частоты механизма оси:
Юр1 =
2 ,/р1=Юр1/(2^); Юр2 =
| ,/р2=®р2/(2^). (3)
2
Примем следующие значения параметров механизма: •1=.3=50, .2=400 кг-м , С12 = С23=8-106 Нм/рад. При таких параметрах механизма по формулам (3) получим /р1=63,7, /р2=71,2 Гц.
Математическое моделирование и спектральный анализ механических переходных процессов. На рис. 2, а приведены результаты математического моделирования механических переходных процессов для трехмассового механизма с одним вращающим моментом М1 при принятых выше параметрах. Кривые 1, 2 и 3 в масштабе 10:1 представляют собой отклонения скоростей масс АО соответственно 01, О2 и О3 от кривой 4, описывающей усредненное движение Оср(0, соответствующее движению эквивалентной одномассовой модели:
^ср (') = (4)
• £
при М^=М1=2,62 Нм и •^=.1+.2+.3 и в масштабе 1:1.
На рис. 2, б представлены результаты выполненного с помощью МаШСаё численного спектрального анализа несинусоидальных кривых 1 и 3 (рис. 2, а). Как видно из графика, в спектре кривых имеют место пики на частотах 63,7 и 71,2 Гц, значения которых совпадают с полученными расчетным путем.
Как видно из рис. 2, а, кривая 2, соответствующая скорости второй массы (труба телескопа), представляет собой практически синусоиду с частотой /р2=71,2 Гц, что подтверждают результаты спектрального анализа, представленные на рис. 2, в. Как показали исследования, такая картина имеет место при симметрии правой и левой полуосей конструкции угломест-ной оси относительно оси трубы телескопа, т.е. при •1=.3 и С12=С23. При существенной асимметрии имеют место несинусоидальные колебания второй массы, и резонансные пики в спектре кривой 2 проявляются на двух частотах.
а) QxlCT4 1/с
7
C
-7 -14
6C
б)
AQxlC-1/с
4
12C
в)
AQx1C-6, 1/с
1С 5
18С
t, мс
4C
6C
Zp1 fp2
8C
4C
6C
fp2 80 f, Гц
f, Гц Рис. 2
На рис. 3, а приведены результаты математического моделирования механических переходных процессов для трехмассового механизма с двумя вращающими моментами при принятых выше параметрах механизма. Кривые 1, 2 и 3 представляют собой отклонения скоростей масс AQ соответственно Q1, Q2, Q3 в масштабе 10:1 от графика (4) усредненного движения эквивалентной одномассовой модели в масштабе 1:1 при моменте М^=М1+М2 и М=М2=1,31 Нм.
а)
Qx10-4. 1/с
6
-6
60
б)
Aßx10-6 1/с
32
16
120
в)
Aßx10-1/с
180 t, мс
8
1
J 1
40
60
80
f, Гц
40
60
80
f Гц
Рис. 3
Видно, что кривые 1—3 представляют собой практически синусоиды с частотой_/р2=71,2 Гц, это подтверждают результаты спектрального анализа (рис. 3, б, в). Колебания первой и третьей масс имеют одинаковую амплитуду и фазу, противоположную фазе колебаний второй массы. Такие процессы свойственны механическим переходным процессам в двухмассовом механизме.
0
5
2
0
0
0
4
462 В. А. Толмачев, М. В. Никитина, М. Е. Сергеева
Трехмассовая модель угломестной оси с симметричной системой масс и двумя вращающими моментами при М1=М2 может быть приведена к двухмассовой модели (например, с использованием метода электрической аналогии [7]) с эквивалентами параметрами •*1=.1+.3, .*2=.2, С*12=С12+С23, М*=М1+М2. Коэффициент соотношения масс для нее определяется выражением [6]:
т. = • '1 + • '2 , (5)
• '1
а соотношение амплитуд колебаний скоростей масс находится как
О 'т1
О 'т2
= у' - 1. (6)
При принятых параметрах механизма оси имеем •*1=100, .*2=400 кг-м2, С*12= =1,6-105 Нм/рад и, следовательно, у'э=5, _/*р=/р2=71,2 Гц, й'т1/й'т2=4. Справедливость последних цифр легко установить по рис. 3, а.
Как показали проведенные исследования, при М1 ^ М2 колебания первой и третьей масс оси становятся несинусоидальными и в спектрах кривых 1 и 3 (рис. 2, б) появляются пики на двух частотах _/р1=63,7, _/р2=71,2 Гц. Колебания второй массы сохраняются синусоидальными и имеют место на высшей частоте — /р2. При существенной асимметрии даже для М1=М2 имеют место несинусоидальные колебания второй массы и резонансные пики в ее спектрах на двух частотах /*р1 и /*р2. Колебания первой и третьей масс при этом сохраняются синусоидальными, но имеют разные частоты.
Заключение. Результаты анализа механических переходных процессов в однодвига-тельных и двухдвигательных механических подсистемах следящих электроприводов с трех-массовыми исполнительными осями, полученные с использованием предложенных векторно-матричной математической модели и расчетных соотношений для определения резонансных частот, позволяют сформулировать некоторые рекомендации к системам управления двух-двигательных следящих электроприводов угломестных осей ОПУ телескопа.
1. Двухдвигательные механизмы целесообразно применять при построении систем управления следящих электроприводов трехмассовых угломестных осей симметричного типа.
2. При построении систем управления следует ориентироваться на структурные решения, обеспечивающие равенство в статических и динамических режимах работы вращающих моментов, действующих на правые и левые полуоси.
список литературы
1. Васильев В. Н., Томасов В. С., Шаргородский В. Д., Садовников М. А. Состояние и перспективы развития прецизионных электроприводов комплексов высокоточных наблюдений // Изв. вузов. Приборостроение. 2008. Т. 51, № 6. С. 5—12.
2. Томасов В. С., Денисов К. М., Толмачев В. А. Следящие электроприводы систем наведения оптико-механических комплексов нового поколения. Проблемы и достижения // Тр. V Междунар. (XVI Всерос.) конф. по автоматизированному электроприводу АЭП-2007. СПб, 2007. С. 175—177.
3. Исупов А. Н., Исупов К. С., Храмов С. Н. Определение амплитудно-частотных характеристик альт-азимутального опорно-поворотного устройства крупногабаритного наземного телескопа // Изв. вузов. Приборостроение. 2008. Т. 51, № 6. С. 38—44.
4. Толмачев В. А. Синтез следящего электропривода оси опорно-поворотного устройства // Изв. вузов. Приборостроение. 2008. Т. 51, № 6. С. 68—72.
5. Борцов Ю. А., Соколовский Г. Г. Автоматизированный электропривод с упругими связями. СПб: Энергоатомиздат, 1992. 288 с.
6. Ключев В. И. Теория электропривода: Учеб. для вузов. М.: Энергоатомиздат, 2001. 704 с.
7. Саушев А. В., Шошмин В. А. Моделирование многомассовых механических систем электроприводов методом электрической аналогии // Вестн. гос. ун-та морского и речного флота им. адмирала С. О. Макарова. 2010. № 4. С. 57—64.
Сведения об авторах
Валерий Александрович Толмачев — канд. техн. наук доцент; Университет ИТМО, кафедра электротехники и прецизионных электромеханических систем; E-mail: [email protected] Мария Владимировна Никитина — канд. техн. наук; Университет ИТМО, кафедра электротехники и
прецизионных электромеханических систем; E-mail: [email protected]
Мария Евгеньевна Сергеева — аспирант; Университет ИТМО, кафедра электротехники и прецизи-
онных электромеханических систем; E-mail: [email protected]
Рекомендована кафедрой Поступила в редакцию
электротехники и прецизионных 27.11.14 г.
электромеханических систем
Ссылка для цитирования: Толмачев В. А., Никитина М. В., Сергеева М. Е. Анализ динамических свойств трехмассовых исполнительных осей с одним и двумя вращающими моментами // Изв. вузов. Приборостроение. 2015. Т. 58, № 6. С. 458—463.
ANALYSIS OF DYNAMIC PROPERTIES OF THREE-MASS EXECUTIVE AXES WITH ONE AND TWO TORQUES
V. A. Tolmachev, M. V. Nikitina, M. E. Sergeeva
ITMO University, 197101, Saint Petersburg, Russia E-mail: [email protected]
Mechanical transient processes in single-engine and two-engine three-mass mechanisms of executive axes of servo drives of support-rotating systems are considered. Vector - matrix mathematical models describing the processes are proposed. Formulas for calculation of resonance frequencies are presented, the model data are confirmed by results of spectral analysis of the transient processes carried out with the use of MatLab package. Transformation of initial models with the use of electric analogy method allows calculating parameters of equivalent two-mass models of two-engine servo drives with three-mass mechanisms of the axes. The results provide a basis for selection of structure design and development of rigorous method for synthesis of control systems of telescope servo drives with two-engine three-mass mechanisms of the axes.
Keywords: executive axis, three-mass mechanism, mathematical model, mechanical transient processes, resonance frequencies, spectral analysis.
Data on authors
Valery A. Tolmachev — PhD, Associate Professor; ITMO University, Department of Electrical
Engineering and Precision Electromechanical Systems; E-mail: [email protected] Maria V. Nikitina — PhD; ITMO University, Department of Electrical Engineering and Preci-
sion Electromechanical Systems; E-mail: [email protected] Maria E. Sergeeva — Post-Graduate Student; ITMO University, Department of Electrical En-
gineering and Precision Electromechanical Systems; E-mail: [email protected]
Reference for citation: Tolmachev V.A., Nikitina M. V., Sergeeva M. E. Analysis of dynamic properties of three-mass executive axes with one and two torques // Izvestiya Vysshikh Uchebnykh Zavedeniy. Priboros-troenie. 2015. Vol. 58, N 6. P. 458—453 (in Russian).
DOI: 10.17586/0021-3454-2015-58-6-458-463