CC BY
31
программы линейно возрастает с увеличением количества элементов, на обработку матрицы в направлении от очага пожара с применением теории графов.
список литературы
1. Шахраманьян М. А., Нигметов Г. М. Методика оперативной оценки последствий лесных пожаров. М.: ВНИИ ГОиЧС, 2001. 32 с.
2. Тарасевич Ю. Ю. Перколяция: теория, приложения, алгоритмы: Учебное пособие. М.: Едиториал УРСС, 2002. 112 с.
Сергей Алексеевич Астафьев
Дмитрий Юрьевич Лысенко Анатолий Сергеевич Широков
Рекомендована кафедрой мехатроники
Сведения об авторах аспирант; Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кафедра мехатроники; E-mail: [email protected] ЗАО „Транзас", Санкт-Петербург; руководитель группы Менеджеров проектов; E-mail: [email protected] ЗАО „Транзас", Санкт-Петербург; ведущий инженер-конструктор; E-mail: [email protected]
Поступила в редакцию 29.02.12 г.
УДК 681.518.3
И. И. Калапышина МОДЕЛИРОВАНИЕ ДИНАМИКИ СТАНКА ДЛЯ ОБРАБОТКИ СТЕКОЛ
Рассматривается динамика механизма для обработки оптических стекол. Представлены дифференциальные уравнения движения механизма. Для моделирования уравнений узла шлифовки оптических элементов в пакете MatLab применен редактор дифференциальных уравнений DEE.
Ключевые слова: оптические стекла, дифференциальное уравнение, автоматизированный анализ, имитационное моделирование.
Введение. Для изготовления оптических деталей наблюдательных приборов, фотографических и проекционных аппаратов используются разные типы и марки стекла, включая стекла для светофильтров и кварцевое стекло. Различные классы оптических систем требуют разной точности обработки поверхностей стекол. Оптические системы с дифракционно-ограниченным качеством изображения, которые используются, например, для космических и прецизионных телескопов, биноклей, микроскопов, требуют более высокого класса обработки, чем оптические системы с аберрационно-ограниченным качеством изображений [1]. Так, например, при обработке оптических кристаллов, таких как германий, кремний, флюорит, оптическое производство обязано выдерживать уровень шероховатости поверхности не более 2,5 нм, а в некоторых кристаллах, например фтористом кальции, требования по шероховатости еще выше — 1,1 нм. Однако при изготовлении оптических стекол автоматизированный анализ качества поверхности не производится. В условиях производства проверка поверхностей осуществляется при помощи цеховых образцов различного квалитета, для этого требуется прервать процесс обработки, разорвать кинематическую схему, а затем ее восстановить. Такая технология увеличивает время изготовления одного элемента и повышает вероятность брака. Поэтому целью настоящей работы является поиск метода, который позволит в процессе обработки оценивать текущее качество поверхности при помощи аналитических функций.
Кинематическая схема. На рис. 1, а приведено схематичное изображение шлифовального станка, б — кинематическая схема узла шлифовки оптического стекла.
Ведущее звено механизма — нижняя планшайба (ведущая) 1, на ней размещается обрабатываемая деталь 2. Верхнее звено — ведомая планшайба 3, или абразивный инструмент, он вращается относительно оси 02 за счет двигателя Д2. Также планшайба совершает дополнительный поворот вокруг собственной оси 0121 за счет сил смешанного трения, которые возникают в области контакта. Поводок 4 состоит из Г-образного стержня, на его вход подается вращательный момент от электродвигателя Д2, на выходе жестко закреплен полый цилиндр, вертикальная ось которого совмещена с валом инструмента. Электропривод Д1 вращает ведущее звено и закрепленную на нем жестко, при помощи смол, обрабатываемую деталь (она может находиться в закрепленном положении, а также получать дополнительный вращательный момент от ведущей планшайбы 1). Под действием электропривода Д2 верхняя планшайба вместе с поступающей к ней эмульсией, содержащей абразив, вращается вокруг оси поводка 02 [2]. В результате процессов трения между обрабатываемой поверхностью и ведомой планшайбой обрабатываемая деталь 3 вращается вокруг собственной оси 0121. Также она получает дополнительное вращение за счет ведущей планшайбы 1 и вязкой эмульсии между ними.
Структура дифференциальных уравнений движения узла шлифовки. При составлении дифференциальных уравнений смешанное трение имитировалось с помощью демпфирующих и упругих элементов [3]. Кинематическая схема была интерпретирована в модель, ведущая и ведомая планшайбы, а также обрабатываемая деталь которой обладают свойствами идеальных тел, инерцией и массой (для упрощения они названы: первое, третье и второе тело соответственно). Исключены из расчета движущий момент от Д2 и вращение ведомой планшайбы относительно оси 02. За обобщенные координаты в уравнениях приняты возможные перемещения относительно координатных осей для каждого тела. Между телами размещены элементы, содержащие упругую составляющую и вязкое сопротивление, которое характеризуется коэффициентом демпфирования [3]. Дифференциальные уравнения, описывающие динамику узла шлифовки, станка для обработки оптических стекол, представляют собой систему уравнений первого порядка:
Поводок
2 Д2
Привод Д2
Рис. 1
(1)
Jiqi = f(t) + Mdv - aq1 - к41, J2 q2 = f (t) - a(q2 - 41) - к (2 - 41), J343 = f (t) - a(43 - Ч2) - к(43 - 42 ),
m44 = f (t) + ^трв - ^тр ,
где Jn — момент инерции n-го тела; Mdv — движущий момент двигателя Д1; 41 = ф®2 — вращение ведущей планшайбы относительно оси OZ; 42 = ф22 — вращение обрабатываемой детали относительно оси OZ; 43 = ф®121 — вращение ведомой планшайбы относительно оси O1Z1; 44 = Х2 — перемещение обрабатываемой детали; 4n — обобщенная координата n-го тела, f (t) — внешнее воздействие, F^ — сила вязкого трения, Frp — сила сухого трения.
Для моделирования уравнений движения узла шлифовки оптических элементов в пакете MatLab применен редактор дифференциальных уравнений DEE, схема моделирования представлена на рис. 2 [4], для ее построения были введены соответствующие „машинные" переменные:
41 = x(1); 41 = x(2); ¿12 = x(3); 42 = x(4);
43 = x(5); 43 = x(6);
44 = x(7); 44 = x(8); u (1) = f (t) = ot + 5).
В результате замены переменных система уравнений (1) приобрела следующий вид:
41 = (u(1) + Mdv - ax(1) - kx(2)) / J1,
42 = (u (1) - a(x(3) - x(1)) - к (x(4) - x(2))) / J2, J3<43 = (u(1) - a(x(5) - x(3)) - к (x(6) - x(4))) / J3,
m44 = u(1) + ^рв - FjP •
(2)
(3)
Рис. 2
Результаты решения уравнений в виде фазовых пространств представлены в таблице.
Заключение. Решение системы дифференциальных уравнений (1) позволило наглядно представить эволюцию поведения скоростей ведущей планшайбы, обрабатываемой детали и ведомой планшайбы при изменении коэффициента демпфирования и упругой составляющей, а следовательно, и параметров контактного взаимодействия между телами. В процессе моделирования динамика системы изменялась от устойчивого состояния, когда фазовые пространства подсистем представляли собой центр, до неустойчивого, когда фазовые пространства подсистем видоизменялись в устойчивый, а затем в неустойчивый фокусы под воздействием внешней синусоидальной нагрузки. Таким образом формируется база данных теоретических и практических результатов, использование которой позволит прогнозировать состояния поверхности в режиме реального времени.
список литературы
1. Кирилловский В. К. Современные оптические исследования и измерения. СПб: Лань, 2010. 303 с.
2. Горбач В. Л. Кинематика рабочих органов оптических шлифовально-полировальных станков. М., 1958. 110 с.
3. Мусалимов В. М., Валетов В. А. Динамика фрикционного взаимодействия. СПб: СПбГУ ИТМО, 2006. 191 с.
4. Лазарев Ю. В. Моделирование в системе Matlab. М., 2000. 511 с.
Сведения об авторе
Ирина Ивановна Калапышина — аспирант; Санкт-Петербургский национальный исследовательский
университет информационных технологий, механики и оптики, кафедра мехатроники; E-mail: [email protected]
Рекомендована кафедрой Поступила в редакцию
мехатроники 29.02.12 г.
