скачать На правах рукописи Пешев Ярослав Иванович СТРУКТУРНОЕ МОДЕЛИРОВАНИЕ И АЛГОРИТМ УПРАВЛЕНИЯ ПОДВИЖНЫМИ ОРГАНАМИ ОБРАБАТЫВАЮЩЕГО ЦЕНТРА С ЦЕЛЬЮ ОБЕСПЕЧЕНИЯ ТРЕБУЕМОЙ ДИНАМИЧЕСКОЙ ТОЧНОСТИ ИНСТРУМЕНТА ПРИ ВОСПРОИЗВЕДЕНИИ ЗАДАННОЙ ТРАЕКТОРИИ Специальность 05.13.06 – Автоматизация и управление технологическими процессами и производствами (промышленность) Автореферат диссертации на соискание ученой степени кандидата технических наук Самара – 2010 Работа выполнена на кафедре «Электропривод и промышленная автоматика» Государственного образовательного учреждения высшего профессионального образования «Самарский государственный технический университет». Научный руководитель: доктор технических наук, профессор ^ Официальные оппоненты: доктор технических наук, профессор Галицков Станислав Яковлевич кандидат технических наук, доцент ^ Ведущая организация: ГОУ ВПО «Тольяттинский государственный университет» (г. Тольятти) Защита диссертации состоится «28» декабря 2010 г. в 09 часов 00 мин на заседании диссертационного совета Д 212.217.03 ГОУ ВПО «Самарский государственный технический университет» по адресу: г. Самара, ул. Первомайская, 18, 1 корпус, ауд. №4 (Учебный центр «СамГТУ Электрощит»). Отзывы по данной работе в двух экземплярах, заверенные печатью, просим направлять по адресу: Россия, 443100, Самара, ул. Молодогвардейская, 244, Главный корпус на имя ученого секретаря диссертационного совета Д 212.217.03; факс: (846) 278-44-00; e-mail: D21221703@list.ru С диссертацией можно ознакомиться в библиотеке Самарского государственного технического университета (ул. Первомайская, 18). Автореферат разослан __ ноября 2010 г. Ученый секретарь диссертационного совета ^ Д 212.217.03 Общая характеристика работы Актуальность темы. Требования к качеству обработки деталей для последующих технологических операций непрерывно возрастают, так, в 1990 году величина поля допуска составляла значение 5 мкм, а в настоящее время стремительно приблизилась к величине 1-3 мкм. В таких условиях крайне актуальными становятся вопросы разработки многооперационных координатно-расточных станков (обрабатывающих центров) класса точности C, способных повысить качество выпускаемой продукции и её конкурентоспособность, как на региональном, так и на мировом рынках. Реализуемые на подобных станках технологические процессы и операции (металлообработка, контрольные измерения) предъявляют высокие требования к точности позиционирования рабочих органов, а также к динамическим процессам перемещения подвижных узлов заключающиеся в отсутствии ударов в кинематических звеньях электропривода подачи при разгоне, торможении, позиционировании; монотонность позиционирования; обеспечение требуемых значений скоростной ошибки и ошибки по ускорению при одновременном обеспечении высоких скоростей перемещения подвижных узлов. Стремление к цели обеспечения требуемой точности положения инструмента в рабочем пространстве прецизионного металлорежущего станка, приводит к необходимости решения целого ряда сложных и взаимосвязанных проблем таких как: обеспечение динамики перемещения и позиционирования рабочих органов линейных интерполирующих координат; компенсация погрешностей от взаимовлияния интерполирующих координат; обеспечение требуемой жёсткости несущей системы станка (станина, стойка, механические крепления несущих элементов); компенсация разворота и всплывания подвижных узлов на направляющих; компенсация температурных деформаций. Значительный вклад в области повышения точности прецизионных координатно-расточных станков за счёт введения системы автоматического управления (САУ) подвижными узлами и несущими элементами был внесён такими учёными как: Ж.С. Равва, К.В. Вотинов, В.А. Кудинов, Д.Н. Решетов, С.Я. Галицков, А.П. Соколовский, С.Я., В.Г. Болтянский, Н.Н. Красовский, Г.М. Уланов, З.М. Левин. Вклад в решение проблем по улучшению динамических характеристик станков внесён зарубежными исследователями: Г. Шлезингером, Ф. Кенигсбергом, И. Тлусты. Динамические погрешности положения инструмента в рабочем пространстве станка имеют крайне сложную природу возникновения, и зачастую носят случайный характер. В условиях работ по развитию нанотехнологий этот вид погрешностей приобретает особое значение, поэтому данная диссертационная работа актуальна и своевременна. Работа выполнена в рамках федеральной государственной программы «Развитие точного машиностроения и станкостроения». ^ является влияние динамики САУ линейными интерполирующими координатами и САУ устройством автоматической смены инструмента (УАСИ) на точность положения инструмента в рабочем пространстве станка и эксплуатационную надёжность по точности. ^ является формирование алгоритмов управления траекториями перемещений подвижных органов линейных интерполирующих координат и руки механизма смены инструмента, которые обеспечивают требуемую величину динамической ошибки положения инструмента в рабочем пространстве станка. ^ Структурное представление процесса формирования динамической погрешности положения инструмента в рабочем пространстве станка. Разработка управляющих алгоритмов для САУ линейными интерполирующими координатами и УАСИ, которые обеспечивают заданную динамическую ошибку положения инструмента в рабочем пространстве и требуемую долговечность станка по точности. Для достижения цели, поставленной в работе, необходимо решить следующие задачи:
^ теоретические исследования базируются на методах теории автоматического управления, классической механики, линейной алгебры, теорий вероятности и математической статистики. Вопросы анализа и синтеза систем управления решались методами дифференциального и интегрального исчисления, компьютерного моделирования, численными методами. Метод экспериментального исследования на стенде-станке использовался для получения исходных данных, проведения и уточнения результатов теоретического анализа. ^ заключается в следующем: Диссертационная работа расширяет и углубляет теоретические представления о формировании статических и динамических погрешностей положения инструмента в рабочем пространстве станка, которые определяют качество обработки деталей. В диссертации получены следующие основные результаты, отличающиеся научной новизной:
^ состоит в следующем:
^ . Основные результаты теоретических и экспериментальных исследований внедрены:
Апробация работы. Основные положения и результаты исследований докладывались и обсуждались на следующих конференциях:
Публикации. По теме диссертации опубликовано 6 работ, в том числе 3 из них – в журналах, рекомендованных ВАК РФ (журнал СТИН №10-2009 и сборник серии «Технические науки» СамГТУ №2(22) – 2008, СамГТУ №1(23) – 2009), 3 публикации – в трудах и материалах международных научно-технических конференций. ^ :
^ Диссертационная работа состоит из введения, четырех глав, заключения, библиографического списка из 58 наименований и приложения. Основной текст работы изложен на 206 страницах, диссертация содержит: 95 рисунков, 11 таблиц, 1 приложение, библиографический список на 6 страницах. ^ Во введении в краткой форме изложены актуальность проведённого исследования; научная новизна; практическая значимость диссертационной работы; перечень задач, решаемых в диссертации; приводятся основные положения, выносимые на защиту. ^ показана роль прецизионных многооперационных станков (обрабатывающих центров) в автоматизации технологических процессов серийного и мелкосерийного производства с учётом особенностей развития российского станкостроения. Представлены анализ природы возникновения и классификация ошибок положения режущего инструмента в рабочем пространстве станка. Формулируются цели и задачи исследования. ^ разработана математическая модель формирования погрешности положения инструмента в рабочем пространстве станка. Первый этап синтеза математической модели формирования погрешности положения инструмента в рабочем пространстве станка состоит в представлении каждой линейной интерполирующей координаты обрабатывающего центра в виде автономного сепаратного канала, который учитывает как упругие, инерционные и диссипативные процессы в серводвигателе, а также передаточных и исполнительного механизмов, так и свойства системы автоматического управления положением исполнительного подвижного узла. Описание упруго-диссипативных процессов линейной интерполирующей координаты удобно представить в виде расчётной схемы. ![]() Рисунок 1 Для линейной интерполирующей координаты X упомянутая выше расчётная схема (см. рисунок 1) представлена в виде последовательной цепи четырёх упруго-диссипативных и трёх инерционных элементов, причём C1X, C2X и Д1X, Д2X – крутильная жёсткость и коэффициент демпфирования соединительной муфты и винта, соответственно; C3X, C4X и Д3X, Д4X – линейная жёсткость и коэффициент демпфирования шариковинтовой пары (ШВП) и упорного подшипника со стыками; mX – масса рабочего органа (совместно с массой обрабатываемой заготовки); FТР X – сила трения в направляющих; ![]() С использованием вышеуказанной расчётной схемы и основного уравнения движения электропривода получим систему операторных уравнений (1), описывающую динамические процессы механической части линейной интерполирующей координаты X: ![]() В системе операторных уравнений (1) введены следующие обозначения: ![]() ![]() ![]() ![]() ![]() ![]() где ![]() ![]() ![]() Передаточная функция (2) по управляющему воздействию синхронного серводвигателя с постоянными магнитами синтезирована с применением известной модели обобщенного двухфазного синхронного серводвигателя. В (2) введены следующие обозначения: m1 – число фаз статора серводвигателя, Zn – число пар полюсов, J – приведённый к валу серводвигателя момент инерции; Mс – момент сил сопротивления; ![]() ![]() ![]() ![]() ![]() ![]() Таким образом, математическая модель автономного сепаратного канала линейной интерполирующей координаты X синтезирована на основе методики синтеза САУ по идеологии систем подчинённого регулирования, системы операторных уравнений (1) и передаточной функции (2). Методика синтеза математической модели автономного сепаратного канала линейной интерполирующей координаты ^ имеет полную аналогию с методикой синтеза математической модели автономного сепаратного канала линейной интерполирующей координаты X. ![]() Рисунок 2 Расчётная схема линейной интерполирующей координаты Z, изображена на рисунке 2, в которой введены обозначения: C1Z, C2Z и Д1Z, Д2Z – крутильная жёсткость и коэффициент демпфирования соединительной муфты и винта, соответственно; C3Z, C4Z и Д3Z, Д4Z – линейная жёсткость и коэффициент демпфирования шариковинтовой пары и упорного подшипника со стыками; ![]() ![]() ![]() Динамические процессы механической части линейной интерполирующей координаты Z описываются системой операторных уравнений (3), в которой: ![]() ![]() ![]() ![]() ![]() ![]() Влияние противовеса, разгружающего по массе шпиндельную бабку, на динамические показатели линейной интерполирующей координаты Z, отражено системой операторных уравнений (4): ![]() Математическая модель автономного сепаратного канала линейной интерполирующей координаты Z синтезирована на основе методики представления синхронного серводвигателя с постоянными магнитами на роторе как объекта управления (2), методики синтеза САУ по идеологии систем подчинённого регулирования, а также систем операторных уравнений (3) и (4). ![]() Рисунок 3 На рисунке 3 изображена структурная схема автономного сепаратного канала линейной интерполирующей координаты Z. Указанная структурная схема включает в себя звенья, учитывающие упругие, инерционные, диссипативные процессы в серводвигателе подачи, передаточных и исполнительных механизмов линейной интерполирующей координаты Z; контур стабилизации тока серводвигателя подачи; контур стабилизации частоты вращения вала серводвигателя подачи и контур положения шпиндельной бабки. Элементы, показанные на рисунке 3 пунктиром, учитывают наличие противовеса, обеспечивающего разгрузку по массе подвижного органа (шпиндельной бабки). На рисунке 3 обозначены: ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() Реализация возможности моделирования межканальных связей линейных интерполирующих координат через несущие элементы станка, требует рассмотрения каждого автономного сепаратного канала с позиций формирования собственной помехи линейной интерполирующей координатой и реакции подвижного органа линейной интерполирующей координаты на внешние силовые возмущающие воздействия. Структурная схема канала реакции линейной интерполирующей координаты X на внешние возмущающие воздействия изображена на рисунке 4, где ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() Рисунок 4 Структурно, каналы реакции линейных интерполирующих координат Y и Z на внешние силовые возмущающие воздействия полностью аналогичны каналу реакции линейной интерполирующей координаты X. Собственная помеха линейной интерполирующей координаты возникают при перемещении подвижного узла этой координаты и обусловлена неидеальностью геометрии элементов ШВП, направляющих, наличием на направляющих посторонних твёрдых частиц (загрязнений), а также другими неучтёнными факторами. В процессе сравнения, и статистической обработки результатов экспериментов, согласно положениям теории инженерного эксперимента, было предложено представлять собственную помеху ![]() ![]() - составляющая ![]() ![]() ![]() - составляющая ![]() Данная составляющая может быть представлен в виде сигнала ![]() ![]() ![]() С целью синтеза генератора составляющей ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() Математическую модель формирования составляющих собственной помехи линейной интерполирующей координаты X иллюстрирует структурная схема на рисунке 5, на которой обозначены: ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() Рисунок 5 Математическая модель формирования составляющих собственной помехи линейной интерполирующей координаты ^ структурно полностью аналогична математической модели формирования составляющих собственной помехи линейной интерполирующей координаты X. Разгон/торможение в процессе перемещения шпиндельной бабки способствует возникновению сильных гармонических колебаний противовеса преимущественно в вертикальной плоскости. Это обусловлено тем, что груз, подвешенный на двух тросах, соединённых с подвижным органом Z, являет собой по существу пружинный маятник. Поэтому, очевидно, что имеет место быть составляющая ![]() ![]() Представляя противовес в виде пружинного маятника, получим выражение математической модели формирования составляющей ![]() ![]() ![]() Рисунок 6 Структурная схема формирования составляющих собственной помехи линейной интерполирующей координатой Z представлена на рисунке 6, на которой обозначены: ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() При синтезе многосвязанной модели формирования погрешности положения инструмента в рабочем пространстве станка учитывались конструктивные особенности взаимного механического соединения линейных интерполирующих координат станка. Силовые возмущающие воздействия от УАСИ были определены экспериментально на разных этапах технологической операции автоматической смены инструмента. В результате чего можно выделить четыре источника силовых возмущающих воздействий: ![]() ![]() ![]() ![]() Система операторных уравнений, отражающая взаимную связь автономных сепаратных каналов и УАСИ: ![]() В (5) дифференцирующие звенья ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() Многосвязанная модель формирования погрешностей положения инструмента в рабочем пространстве станка, синтезированная с применением системы операторных уравнений (5), показана на рисунке 7. Блок «X» включает в себя автономный сепаратный канал линейной интерполирующей координаты X, канал реакции линейной интерполирующей координаты X на внешние возмущающие воздействия, канал формирования составляющих собственной помехи линейной интерполирующей координаты X. Состав блоков «Y» и «Z» аналогичен блоку «X». ![]() Рисунок 7 ^ с использованием многосвязанной модели, структурная схема которой изображена на рисунке 7, осуществляется имитационное моделирование формирования погрешности положения инструмента в рабочем пространстве станка при отработке линейной интерполяции; при воздействии на подвижные органы (стол, салазки, шпиндельная бабка) внешнего динамического силового возмущающего воздействия; при технологическом процессе шлифования с осцилляцией шпиндельной бабкой; при технологической операции автоматической смены инструмента. Исследование возникновения погрешности положения инструмента в рабочем пространстве станка при отработке круговой и линейной видов интерполяций позволило вывести критерий формирования ошибки положения стола и салазок при отработке линейной и круговой интерполяции. Указанный критерий позволяет установить аналитическую зависимость между абсолютным максимальным значением составляющей собственной помехи линейной интерполирующей координаты и угловым ускорением вращения вала серводвигателя подачи подвижного органа линейной интерполирующей координаты, а также линейной скоростью перемещения подвижного органа линейной интерполирующей координаты. Результаты экспериментов и моделирования показали следующие результаты: для участка разгона/торможения вала серводвигателя подачи допустимый диапазон изменения ускорений (0-1,5 м/c2). При ускорении 1,5 м/c2 для линейной интерполирующей координаты X, максимальное отклонение салазок от заданного положения составило – 1,8 мкм, шпиндельной бабки – 0,6 мкм. При ускорении 1,5 м/c2 для линейной интерполирующей координаты Y, максимальное отклонение стола от заданного положения составило – 1,1 мкм, шпиндельной бабки – 0,8 мкм. При движении стола со скоростью 5400 мм/мин, максимальное отклонение салазок от заданного положения составило – 2,3 мкм, шпиндельной бабки – 0,4 мкм. Движение салазок со скоростью 5400 мм/мин, взывает максимальное отклонение стола от заданного положения – 1,5 мкм, шпиндельной бабки – 0,8 мкм. Проведённый анализ формирования погрешности положения инструмента в рабочем пространстве станка при воздействии на подвижные органы внешнего динамического силового возмущающего воздействия позволил найти допустимый диапазон максимальных значений динамической составляющей силы резания ![]() Результатом изучения законов формирования ошибки положения инструмента в рабочем пространстве станка при осцилляции шпиндельной бабкой является разработка способов уменьшения силовых возмущающих воздействий от противовеса в режиме осцилляции и, как следствие, обеспечение требуемых точностных показателей обрабатывающего центра. Путём исследования возмущающих воздействий от УАСИ была доказана необходимость применения САУ УАСИ, осуществляющая режим «мягкая смена инструмента». В указанном режиме управление электроприводом УАСИ осуществляется специально разработанной системой, которая обеспечивает плавное замедление перемещения руки при захвате, выеме, посадке инструмента и при сходе руки с инструмента. В режиме «мягкая смена инструмента» в десятки раз снижается максимальное отклонение рабочих органов линейных интерполирующих координат X, Y и Z по сравнению с традиционным режимом автоматической смены инструмента. В аварийном режиме автоматической смены инструмента применение САУ УАСИ позволило снизить приблизительно в 90 раз максимальное отклонение рабочего органа линейной интерполирующей координаты Z. ^ получены критерии формирования траекторий перемещения формообразующих линейных интерполирующих координат X и Y для линейной и круговой видов интерполяций. Критерии позволяют определить значения угловых ускорений серводвигателей подач и линейных скоростей перемещения подвижных узлов, при которых обеспечивается заявленный класс точности C. Максимальная ошибка положения стола при отработке линейной интерполяции линейными интерполирующими координатами Y и/или Z и установившейся частоте вращения вала серводвигателя подачи (Y и Z): ![]() Уравнение (6) справедливо для следующих диапазонов изменения линейных скоростей перемещения салазок и шпиндельной бабки: ![]() ![]() Максимальная ошибка положения стола при отработке линейной интерполяции линейными интерполирующими координатами Y и/или Z и разгоне/торможении серводвигателя подачи (Y и Z): ![]() Уравнение (7) справедливо для следующих диапазонов установившихся линейных скоростей перемещения салазок и шпиндельной бабки: ![]() ![]() Максимальная ошибка положения салазок при отработке линейной интерполяции линейными интерполирующими координатами X и/или Z и установившейся частоте вращения вала серводвигателя подачи (X и Z): ![]() Уравнение (8) справедливо для следующих диапазонов изменения линейных скоростей перемещения стола и шпиндельной бабки: ![]() ![]() Максимальная ошибка положения салазок при отработке линейной интерполяции линейными интерполирующими координатами X и/или Z и разгоне/торможении серводвигателя подачи (X и Z): ![]() Уравнение (9) справедливо для следующих диапазонов установившихся линейных скоростей перемещения стола и шпиндельной бабки: ![]() ![]() Абсолютное значение максимальной ошибки положения стола при отработке круговой интерполяции линейными интерполирующими координатами (X и Y): ![]() Абсолютное значение максимальной ошибки положения салазок при отработке круговой интерполяции линейными интерполирующими координатами (X и Y): ![]() Уравнения (10) и (11) справедливы, если в процессе отработки круговой интерполяции не выполняется условие: ![]() Анализ и моделирование возмущающих факторов в процессе отработки технологической операции автоматической смены инструмента позволили синтезировать алгоритмы перемещения руки механизма смены инструмента, при которых возмущающие воздействия от УАСИ не оказывают существенного влияния на эксплуатационную надёжность станка по точности. Работа выполнена на кафедре «Электропривод и промышленная автоматика» Самарского государственного технического университета в период обучения в аспирантуре с 2006 по 2009 г. Заключение Основные научные и практические результаты, достигнутые в работе, заключаются в следующем:
^
Автореферат отпечатан с разрешения диссертационного совета Д212.217.03 ГОУ ВПО «Самарский государственный технический университет» (протокол №15 от 22 ноября 2010 г.) Тираж 100 экз. Заказ №93 Отпечатано на ризографе. ГОУ ВПО «Самарский государственный технический университет» Отдел типографии и оперативной полиграфии 443100, г. Самара, ул. Молодогвардейская, 244
|