УДК 621.865.8-83
Н.Ф. КАРНАУХОВ, Ю.В. ПУДОВА, М.Н. ФИЛИМОНОВ
МОДЕЛИРОВАНИЕ КОЛЕБАНИЙ ЗВЕНА ТЕХНОЛОГИЧЕСКОЙ МАШИНЫ В РЕЖИМЕ ДИНАМИЧЕСКОГО ТОРМОЖЕНИЯ АСИНХРОННОГО ДВИГАТЕЛЯ ПРИ ЧАСТОТНОМ УПРАВЛЕНИИ
Моделирование колебаний звена технологической машины (ТМ) с асинхронным электроприводом частотного управления позволяет выявить причины возникновения возмущающих воздействий, определить способы гашения колебаний. Предлагается «ступенчатое» программное регулирование тока подмагничивания статора асинхронного двигателя (АД) при формировании тормозного момента на валу с учетом амплитуды и фазы колебаний, параметры которых определяются виртуальным фазочастотометром.
Ключевые слова: звено подвижности ТМ, амплитуда и фаза колебаний, двухтоковое торможение АД, виртуальный измеритель колебаний.
Введение. В режиме разгона, торможения частотного асинхронного электропривода (ЭП) исполнительного механизма (ИМ) технологической машины (ТМ) в подвижных звеньях могут возникать вынужденные колебательные процессы, отрицательно влияющие на точностные показатели технологического оборудования (ТО), долговечность и эксплуатационные возможности средств автоматизации [1, 2]. Проанализированы причины возникновения колебательных процессов и определены их виды в ТМ [3, 4]. В зависимости от характера изучаемого явления и цели исследования колебаний в ТМ число степеней свободы движения звена ИМ определяют по расчетной схеме, приближенно отражающей конструктивные особенности исследуемого узла [4, 5]. Поставлена задача гашения колебаний звена подвижности ИМ с одной степенью свободы, в частности, захватного устройства (ЗУ) промышленного робота (ПР) с частотно-управляемым электроприводом (ЧУЭП) на базе короткозамкнутого АД.
Постановка задачи. При проектировании новых высокопроизводительных ТМ на базе ЧУЭП, функционирующего по схеме автономный инвертор напряжения (АИН) - асинхронный двигатель (АД), появляется необходимость в предварительной оценке параметров колебательности упругой системы: безударные движения ИМ в начальный период разгона и более сложные процессы торможения, обеспечивающие остановку ИМ в заданном положении с учетом прогнозируемой точности. Колебательность движения ИМ (звена подвижности ТМ) зависит от ряда факторов: массы звена подвижности, скорости перемещения и интенсивности торможения системы [1, 3]. Для оценки параметров колебаний отдельных звеньев ТМ обычно используют специальную измерительную аппаратуру, которую можно заменить виртуальными моделями инструментального оснащения и расширенных возможностей, выполненных в программном пакете Matlab+Simulink [6]. Математическая модель. Для поиска метода гашения колебаний выходного звена подвижности
- захватного устройства (ЗУ) с консольным креплением в направляющей цилиндрической системы координат ПР «Электроника НЦТМ-01» была составлена упрощенная расчетная схема (рис.1).
а 1 а
I
а)
р q
А < 1111 Ї 1 1 ї Ї 1 1 ї )
О в а F ) в
*7 Г* РА ** Ь 1 X
Рис.1. Упрощенная расчетная схема взаимодействующих усилий в направляющих звена подвижности ПР
При анализе расчетной схемы принято равномерное распределение массы руки ПР, а массы ЗУ и транспортируемого груза сосредоточены в одной точке. Рассмотрены основные силовые соотношения в схеме (рис.1) с определением реакций в опорах при нулевых значениях А и В изгибающих моментов относительно точек А и В с учетом параметров [1]: а - вылет консоли направляющей (ход объекта манипулирования), I - расстояние между опорами; d - диаметр направляющей; f- коэффициент трения между направляющей и втулкой конструкции. \ - ускоре-
ние свободного падения, q - распределенная нагрузка звена подвижности; G - вес груза;
Q1, Ql - силы трения в опорах А и В соответственно, F- движущая сила.
Значения реакций опор А и В определяются соотношениями [3]
р 2G(a - /Ф) - q(l2 - а2 ) - 2^[<і(ї+а)
ка 21 '
Яп = G(l + а - /Ф) + q(l + а)(
I + а
М
(1)
(2)
Для поиска деформации руки робота по оси X в точке D (рис.1,б) использовался интеграл Мора [4], что позволило определить перемещение ЗУ под действием приложенных нагрузок
х = , (3)
где к и к2 - коэффициенты при сосредоточенной и распределенной нагрузках, определяемые выражениями:
~Л3 -I ~ л
—1
—[2(-+1) -^ (-+3(1 - /))], 6Е/ а а а
—
І I
0/1 Т7 [2/_ (_~ - (3 + 2_) + 6/Ф(_ + 1)) + (3 + 4-----------^3)] ,
24 Е/ а а а а а а
(4)
(5)
где Е- модуль упругости при прогибе.
Если под сосредоточенной силой G (рис.1) принять алгебраическую сумму силы тяжести и силы инерции, тогда:
G = —тх — mg, (6)
где т- масса перемещаемого груза.
Для перехода к анализу динамики перемещения центра масс ЗУ ось следует направить в сторону действия силы тяжести (по оси Х) и, после некоторых преобразований, можно записать неоднородное дифференциальное уравнение
X + кГ2Х = G ,
где
к
і
и О = g + q—21 , тк„ к
(7)
(8)
11 'М1
где кп, к21 - коэффициенты различия направленности силы трения при четных и нечетных колебаниях ЗУ.
В результате математического моделирования [7] были получены графические зависимости колебаний ЗУ в точке позиционирования (рис.2) для различных перемещаемых масс (груза) без изменения условий торможения.
Рис.2. Зависимость колебаний ЗУ манипулятора в точке позиционирования при нагрузке от 1 до 12 кг т=1, тг=3, т3=5, т4=7, т5=9, т6=12 кг)
4
2
3
Для измерения амплитуд и фаз вынужденных колебаний при моделировании колебательных процессов использовался виртуальный инструментарий - фазочастотометр (рис.3), позволяющий методом сравнения [6] определять искомые параметры гармонических колебаний.
В
SiiiH W ?|\'Н
□
п
Scope 2
_Г
КоЛвбаНИЯ
Un
Вход Смещение входа Выход
Смещение выхода I начала измерения Кол-во циклов на входе
Амплитуда
Eaid/s
Hz
Grad
Rad
Display
Display 1
Фачочастотометр
Scope
N
Рис. 3. Схема подключения фазочастотометра
Процесс измерения включает этапы:
- считывание с входа интервала времени, по истечении которого устанавливаются стационарные колебания;
- измерение длительности одного полупериода колебания входного сигнала и смещения выходных колебаний по отношению к входным колебаниям;
- вычисление фазы (Ф = 2п (1/20) с учетом постоянных смещений входных и0 и выходных ХО колебаний, поскольку последние могут происходить относительно ненулевого уровня.
Результаты выполненных измерений параметров колебательного процесса позволяют построить взаимные зависимости (рис.4).
а)
8 10 12
0 2 4
8 10 12
б)
В)
Рис.4. Взаимные зависимости колебаний: а - амплитуды от массы; б - фазы колебаний от массы; в - амплитуды колебаний от фазы
Полученные зависимости расчетной схемы могут быть использованы в качестве входных переменных для реализации алгоритмов нечеткой логики при управлении частотным электроприводом ИМ ТМ, а также исходными величинами для решения сложной задачи по снижению колебаний ЗУ в процессе работы ПР. Применение ЧУЭП характеризуется возможностью программного управления гашением колебаний ЗУ в процессе регулирования постоянного тока фазы АД в режиме двухтокового динамического торможения [1, 2]. При этом формируются расчетные динамические моменты на валу АД и исключаются возможные удары в механической части ЭМС [1, 3]. Такая система управления АД позволяет (по сигналам датчиков контролируемых величин) в соот-
ветствии с приведенными зависимостями (рис.4) формировать целесообразное торможение звена подвижности в режиме демпфирования колебаний амплитуды ЗУ на определенном уровне (от 0,2 до 0,7 мм), зависимом от чувствительности датчиков и коэффициента усиления сигналов в СУ. В общем случае, процесс торможения АД можно рассматривать как результат взаимодействия эквивалентного двигательного и тормозного моментов на валу, обеспечиваемых равенством намагничивающих сил постоянного тока подмагничивания статора [1, 7] и приведенного значения переменного тока ротора АД. Использование двухтокового способа формирования тормозного момента для гашения колебаний ЗУ при подаче импульсов тока (рис.5) возможно в режиме импульсного регулирования тока подмагничивания, когда интегральная площадь импульса постоянного тока приравнивается к интегральной площади соответствующего периода колебаний звена подвижности, ЗУ. Начало периода формирования импульсов тока подмагничивания статора должно совпадать по времени с периодом вынужденных колебаний звена подвижности, ЗУ, но действовать в противофазе.
Рис.5. Метод гашения колебаний ЗУ ПР в процессе позиционирования
Такое управление электроприводом возможно при программном управлении и наличии силомоментных датчиков на конструктивных элементах ИМ, датчиков контролируемых параметров информационного канала СУ. При этом решение задачи управления ЭП предполагает разработку соответствующего программного обеспечении (ПО) и предсказуемость изменений контролируемых колебаний ЗУ (кривая 1, рис.5), а также возможность формирования импульсов тока подмагничивания обмотки статора АД с оценкой интегрального значения площади каждого импульса, ограниченными экспоненциальной кривой 2 затухания колебательного процесса. Для реализации решения необходимо выбрать датчики измеряемых величин по фазочастотным характеристикам и диапазонам контролируемых параметров. При усложнении задачи исследования и необходимости расширения функциональных возможностей датчиков целесообразно применять известные решения виртуальных приборов. Такая необходимость появляется при управлении, например, ЧУЭП ТМ как сложным объектом в условиях неопределенности с использованием алгоритмов нечеткой логики.
Выводы. Моделирование вынужденных колебаний звена подвижности ТМ, обусловленных процессами «разгона-торможения» электропривода, позволяет:
- определить характерные зависимости вынужденных колебаний от внезапного изменения скорости движения, тормозного усилия и инерционных составляющих перемещаемых масс ИМ;
- разработать направления по совершенствованию конструкции звена подвижности и направляющих ТМ, выбрать рациональное сочетание механических характеристик электропривода и ИМ (ЗУ).
Библиографический список
1. Карнаухов Н.Ф. Демпфирование колебаний захватного устройства промышленного робота в режиме двухтокового динамического торможения асинхронного двигателя при частотном управлении / Н.Ф. Карнаухов, М.Н. Филимонов, Ю.В. Пудова // Вестник ДГТУ. - 2009. - Т. 9. - 308 с.
2. Карнаухов Н.Ф. Особенности формирования двухтокового динамического торможения асинхронного двигателя мехатроннной системы при частотном управлении / Н.Ф. Карнаухов, М.Н. Филимонов, Н.В. Деркачев // Мехатроника-2008: материалы IV междунар. науч.-практ. студ. коллоквиума (Новочеркасск, 18-20 июня 2008) / Южно-Рос. гос. техн. ун-т. - Новочеркасск, 2008.
- С. 17-20.
3. Робототехнические системы в сборочном производстве / под ред. Е.В. Пашкова. - Киев: Выща шк., 1987. - 272 с.
4. Левитский Н.И. Теория механизмов и машин / Н.И. Левитинский. - М.: Наука. Глав. ред. физ.-мат. литературы, 1979.
5. Карнаухов Н.Ф. Энергетические показатели электропривода при частотном способе управления асинхронным двигателем / Н.Ф. Карнаухов, В.А. Прус, М.Н. Филимонов // Динамика технологических систем: тр. VIII междунар. науч.-техн. конф. / Донск. гос. техн. ун-т. - Ростов н/Д, 2007. - Т. III. - 219 с.
6. Данилов А.И. Построение на элементах Simulink измерителей фазы, частоты и амплитуды / А.И. Данилов // Проектирование инженерных и научных приложений в среде MATLAB: тр.
II науч. конф. - 2004.
7. Шрейнер Р.Т. Математическое моделирование электроприводов переменного тока с полупроводниковыми преобразователями частоты / Р.Т. Шрейнер. - Екатеринбург: УРО РАН, 2000.
- 654 с.
References
1. Karnauhov N.F. Dempfirovanie kolebanii zahvatnogo ustroistva promyshlennogo robota v rejime dvuhtokovogo dinamicheskogo tormojeniya asinhronnogo dvigatelya pri chastotnom upravlenii / N.F. Karnauhov, M.N. Filimonov, Yu.V. Pudova // Vestnik DGTU. - 2009. - T. 9. - 308 s. - in Russian.
2. Karnauhov N.F. Osobennosti formirovaniya dvuhtokovogo dinamicheskogo tormojeniya asinhronnogo dvigatelya mehatronnnoi sistemy pri chastotnom upravlenii / N.F. Karnauhov, M.N. Filimonov, N.V. Derkachev // Mehatronika-2008: materialy IV mejdunar. nauch.-prakt. stud. kollokviuma (Novocherkassk, 18-20 iyunya 2008) / Yujno-Ros. gos. tehn. un-t. - Novocherkassk, 2008.
- S. 17-20. - in Russian.
3. Robototehnicheskie sistemy v sborochnom proizvodstve / pod red. E.V. Pashkova. - Kiev: Vyscha shk., 1987. - 272 s. - in Russian.
4. Levitskii N.I. Teoriya mehanizmov i mashin / N.I. Levitinskii. - M.: Nauka. Glav. red. fiz.-mat. literatury, 1979. - in Russian.
5. Karnauhov N.F. Energeticheskie pokazateli elektroprivoda pri chastotnom sposobe upravle-niya asinhronnym dvigatelem / N.F. Karnauhov, V.A. Prus, M.N. Filimonov // Dinamika tehnologicheskih sistem: tr. VIII mejdunar. nauch.-tehn. konf. / Donsk. gos. tehn. un-t. - Rostov n/D, 2007. - T. III. -219 s. - in Russian.
6. Danilov A.I. Postroenie na elementah Simulink izmeritelei fazy, chastoty i amplitudy / A.I. Danilov // Proektirovanie injenernyh i nauchnyh prilojenii v srede MATLAB: tr. II nauch. konf.
- 2004. - in Russian. 7. Shreiner R.T. Matematicheskoe modelirovanie elektroprivodov peremennogo toka s po-luprovodnikovymi preobrazovatelyami chastoty / R.T. Shreiner. - Ekaterinburg: URO RAN, 2000. - 654 s. - in Russian.
Материал поступил в редакцию10.06.10.
N.F. KARNAUKHOV, J.V. PUDOVA, M.N. FILIMONOV
VIBRATION SIMULATION OF TECHNOLOGICAL MACHINE BODY IN THE DYNAMIC BRAKING MODE OF ASYNCHRONOUS MACHINE WITH FREQUENCY CONTROL
Simulation of the vibration level of technological machines (TM) with an asynchronous electric frequency control allows to detect causes of interference effects, to identify oscillation suppression technics. A steplike program control stator current bias in the formation of AS torque on the shaft taking into account the amplitude and phase of the oscillation parameters is determined by the virtual phase-frequency meter.
Key words: TM mobility element, fluctuation amplitude and phase, two current inhibition AS, virtual measuring vibrations.
КАРНАУХОВ Николай Федорович (р. 1938), профессор кафедры «Робототехника и мехатрони-ка» Донского государственного технического университета, кандидат технических наук (1983), доцент (1991). Окончил Уральский электромеханический институт (1963).
Область научных интересов: электропривод и автоматизация промышленных установок и технологических комплексов, полупроводниковые преобразователи частоты промышленных установок и устройств целевого назначения.
Автор 109 публикаций, 4 учебных пособий.
Имеет 7 авторских свидетельств (СССР) и 3 патента на изобретения (РФ).
ПУДОВА Юлия Владимировна, магистрант Донского государственного технического университета. Окончила Донской государственный технический университет (2009).
Область научных интересов: управляемые электроприводы мехатронных систем и их исследование.
Автор 2 публикаций.
ФИЛИМОНОВ Максим Николаевич (р. 1976), заместитель декана факультета «Автоматизация, мехатроника и управление», старший преподаватель кафедры «Робототехника и мехатроника» Донского государственного технического университета. Окончил Донской государственный технический университет (1999).
Область научных интересов: частотно-регулируемый электропривод мехатронных систем и технологического оборудования.
Автор 16 публикаций.
Имеет 2 патента на изобретения (РФ).
Nikolay F. KARNAUKHOV (1938), Professor of the Robotics and Mechatronics Department, DSTU. Candidate of Science in Engineering (1983), Associate Professor (1991). He graduated from Ural Electromechanical Institute (1963).
Research interests: electric drive and automation of the industrial installations and technological complexes, solid-state frequency converters, industrial facilities and special purpose devices.
Author of 109 scientific publications, 7 certificates of authorship of the USSR, 3 patents and 4 manuals.
Julia V. PUDOVA, Undergraduate DSTU. She graduated from the Don State Technical University (2009).
Research interests: guided electric mechatronic systems.
Author of 2 scientific publications.
Maxim N. FILIMONOV (1976), Associate Dean of the Automation, Robotics and Management Faculty, Senior Lecturer of the Robotics and Mechatronics Department, dStu. He graduated from Don State Technical University (1999).
Research interests: frequency-regulated electric mechatronic systems and process equipment.
Author of 16 scientific publications and 2 patents of the RF.