Гебель Елена Сергеевна, Солонин Евгений Владимирович
МОДЕЛИРОВАНИЕ КАК ИНСТРУМЕНТ ИЗУЧЕНИЯ ДИСЦИПЛИНЫ «ТЕОРИЯ МЕХАНИЗМОВ И МАШИН»
Аннотация
В статье приведен опыт применения пакета Model Vision Studium для моделирования задач кинематики плоских рычажных механизмов в процессе преподавания дисциплины «Теория механизмов и машин». Рассмотрены примеры моделирования конкретных механизмов.
Ключевые слова: моделирование, теория механизмов и машин, плоские рычажные механизмы.
Перефразируя известное выражение, скажем: «Для того чтобы научиться моделировать, надо моделировать». Для лектора понимание и осмысление всех тонкостей и подходов к преподаванию дисциплины «Моделирование» приходит с опытом не одного десятилетия. Когда же студентам на «вхождение» в дисциплину отведен только один семестр, то они уже к третьей лекции тонут в новых понятиях «объектно-ориентированного численного моделирования сложных динамических систем», а вслед за этим выступают наружу огрехи в изучении базовых дисциплин: математики, физики, электротехники, теоретической механики. Более того, вдруг выясняется, что в недостаточном объеме или совсем не давались вычислительная математика и численные методы.
Поэтому при преподавании дисциплины необходимо, ориентируясь на студентов, для которых читается курс, подавать основные понятия моделирования через
© Е.С. Гебель, Е.В. Солонин, 2008
призму интересов этой конкретной специальности, стараясь добиться эффекта, когда изучение одного предмета помогает пониманию другого. Так, например, для специальности «Мехатроника» можно начать с того, что близко по духу «мехатро-никам» прежде всего как механикам, то есть использовать классические механизмы, изучаемые в рамках дисциплины «Теория механизмов и машин».
Для быстрого создания таких моделей и проведения экспериментов с ними, записи уравнений, описывающих механизм, получения временных и фазовых диаграмм, визуализации и анимации работы механизма можно использовать различные инструменты моделирования.
В статье приведен опыт применения пакета Model Vision Studium (MVS) для моделирования задач кинематики плоских рычажных механизмов. Выбор именно этого пакета определяется личным пристрастием авторов и тем, что мы начали использовать его для преподавания дисциплины «Моделирование» раньше, чем, например, пакет AnyLogic.
' i. ■ 111........■■'■hi : H'HI-"."'. ii I'i Ij-Л I. ....................................... -.11
Г. hi i tp r J VI м га i ■
ViO^-JO^I-H
. ^. "Г.
tx Л be
■ l-.. ,1
l-ll iitll ll mi l
I -.r-*-
ПШ
tiCE t * 5
H TT ■ .
i 1
Fl T' ?-4 L»" I "I
i ац hit i i in
i J -J. -
. I l—L. Г "J J 1_ Р|П "
„ iWCK ::
Рис. 1. Модель кривошипно-ползунного механизма
Рассмотрим несколько примеров, позволяющих, на наш взгляд, достаточно быстро освоить приемы моделирования таких механизмов.
Задача моделирования простейшего кривошипно-ползунного механизма полностью отражена в окнах проекта (рис. 1). Здесь X, У - координаты кривошипа, переменные Ь и р характеризуют положение коромысла и ползуна.
В примере все параметры модели заданы непосредственно в уравнениях. Длина кривошипа - 0.5, длина коромысла -
Рис. 2. Четырехзвенный механизм
0.8, длина ползуна - 0.3. Угловая скорость вращения кривошипа равна 0.2 рад/с.
При записи модели два первых уравнения записываются просто как проекции кривошипа на оси X и Y. Сложнее объяснить запись двух других уравнений. Впрочем, они становятся очевидны, когда механизм вытягивается в одну линию.
Этот простой пример позволяет быстро объяснить правила построения визуальной модели механизма.
Рассмотрим классическую задачу исследования четырехзвенного механизма (рис. 2).
В монографии [1] построена модель и описано исследование механизма в пакете AnyLogic.
Построим модель механизма в пакете MVS. Исходные параметры механизма: r1 = 180, r2= 90, r3 = 210, r4= 120.
Запишем уравнения, описывающие кинематику механизма.
Координаты точек B и C вычислим через длины звеньев AB и DC и углы.
Гхг = г2 ■ соэ^);
1У в = Г2 ■ эт^).
Гхс = Г! + г4 ■ соэ(^4);
1 Ус = г4 ■ 81П(^4).
Для звена ВС запишем два уравнения:
Гхс = хв + г3 ■ соэ((3);
1 Ус = Ув + гз ■ эт^,).
Решение системы уравнений позволяет определить координаты точек (шарниров) В, С и углы (3, (4, в зависимости от угла поворота ( входного звена АВ. Изменение угла (2 зададим уравнением (2 = Ш.
Если нас явным образом не интересует изменение угла (3, то два последних уравнения можно заменить одним уравнением
гз =уКха - хв )2 +(Уо - Ув )2 .
Уравнения и начальные условия модели показаны на рис. 3.
При заданных параметрах звеньев механизм является кривошипно-коромысло-вым. Для того чтобы исследовать поведе-
ние механизма в зависимости от длины коромысла DC, поместим параметр r4 во внутренние переменные модели. Заметим, что из условия существования механизма при заданных параметрах r1, r2, r3 длина звена DC (r4) не может быть меньше (r2 + r3) - rj =120 (j4 = 0) и больше (r2+ r3) = 300 (j4 = p) (рис. 4).
Изучение четырехзвенного механизма является важным вопросом для понимания теории плоских рычажных механизмов. Большинство плоских механизмов с одной степенью свободы имеют в своем составе подобные четырехзвенные схемы. Пакет MVS позволяет легко изменять параметры звеньев, тем самым наглядно иллюстрировать теорему Грасгофа и условия получения модификаций механизма: кривошипно-коромыслового, двухкриво-шипного, двухкоромыслового, шарнирного параллелограмма, равнобедренного че-тырехзвенного механизма.
После освоения простых моделей можно перейти к более сложным. Рассмотрим задачу моделирования кинематики ры-
Рис. 3. Запись уравнений и задание переменных модели четырехзвенного механизма
Рис. 4. Визуализация модели четырехзвенного механизма
чажного механизма с выстоем, схема которого представлена на рис. 5.
Зависимости, описывающие механизм, составляются исходя из требования разрешимости системы уравнений относительно неизвестных обобщенных координат (в данном случае углов поворота звеньев механизма) и положений шарниров звеньев.
Реализуем в пакете моделирования алгоритм решения задачи кинематического анализа многозвенного рычажного механизма с выстоем.
Рис. 5. Кинематическая схема рычажного механизма IV класса с выстоем
В качестве исходных данных приняты координаты неподвижных точек А, D и Н, кинематические параметры звеньев, найденные на этапе кинематического синтеза механизма: Ь, е, /, к, д, р, г и с -соответственно относительные длины звеньев и сторон звеньев ВЕ, ЕС, ОЕ, ЕК, вН, КС, НЬ, СЬ и CD; ла, Лс и Лн -углы треугольных звеньев ЕОЕ, ЬСК и ОНЬ (треугольное звено однозначно задается двумя сторонами и углом между ними). Длину входного кривошипа АВ принимаем равной единице.
Вычисленные на этапе кинематического синтеза координаты неподвижных шарниров определяют положение рычажного механизма на плоскости и являются «отправными» точками для составления системы уравнений.
Мысленно разъединяем звенья механизма и полагаем, что входное, выходное звено и все треугольные звенья поворачиваются в одну сторону, например против часовой стрелки, и это направление принимаем за положительное.
Обозначим углы поворота звеньев: ф - входного звена АВ, / - треуголь-
ного звена GHL, g - треугольного звена GFE, e - треугольного звена CLK, d -выходного звена CD.
Через координаты неподвижных шарниров А и D записываем координаты точки B входного звена и координаты точки C выходного:
xB = xA + cos(j); xc = xD + c ■ cos(d);
yB = xA + sin(j). yc = yD + c ■ sin(d).
Уравнения для треугольных звеньев записываем через координаты вершины, при которой задан угол. Сначала записываем уравнения для правой (если смотреть из этой вершины) стороны треугольника. Тогда при вращении заданный угол треугольника «добавляется» по направлению вращения стороны, и при записи уравнений для второй стороны берем его с плюсом. Это правило отражается в обозначении треугольных звеньев.
Через координаты точки Н выражаем координаты точек G и L треугольного звена GHL:
JxG = xH + q ■ cos(b); 1 yG = Ун + q ■sin(b);
J xi = xH +t ■cos( b+hH);
1 Уь = Ун + t ■ sin( b +hH ).
Для звена LCK:
Jxb = xc + r ■ cos(e);
1 Уь = Ус + r ■sin(e);
Jx* = xc + P ■cos(e + hc); 1 Ук = Ус + P ■ sin(e +hc).
Для звена EGF:
1xe = xg +e ■ cos(g);
1 Уе = yG +e ■ sin(g);
J= xg + f ■ cos(g+hG);
1 Уе = yG + f ■ sing + hG).
Следует обратить внимание на то, что точка L описывается «с двух сторон» и для нее уравнения записываются дважды.
На заключительном этапе составляются уравнения для звеньев KF и BE, длина которых равна к и b, соответственно:
к = л1(хр-хк)2 -Ук)2;
Ь =л1(хв - хЕ )2+ (ув - Уе )2 .
Эти уравнения называются уравнениями «невязки», хотя на самом деле они как раз «связывают» механизм в одно целое и их выполнение обеспечивает его существование в виде замкнутой кинематической цепи на всем диапазоне возможных положений звеньев.
Заметим, что уравнения для выходного звена СО можно заменить уравнением
С ^{зСс - хо )2 +ус - Уп )2 , тогда угол поворота выходного звена следует искать через одну из тригонометрических функций.
Скорость изменения угла поворота входного звена механизма по-прежнему задается уравнением (p = Юt. Значение скорости зададим в процессе моделирования.
Для решения системы требуется задать начальные условия - приближенные значения углов поворота треугольных звеньев в начальном положении механизма, когда входное звено и звено ВЕ вытянуты в одну линию.
В качестве примера с конкретными значениями параметров модели рассмотрен синтезированный в работе [2] механизм, имеющий следующие кинематические параметры: Ь = 3,75; е = 1; / = 1,5; к = 2,25; q = 1; р = 1; = 10°; г = 1,08; г = 2,48; с = 1,2343; г\н = 97,59°; т)с = 90,56°; Ло = 20°.
Координаты внешних неподвижных шарниров: А(0; 0); Я(8,2260; 2,3787); О(6,2032; -3,5702).
Для представления результатов решения задачи кинематического анализа используется трехмерная анимация работы механизма (рис. 6). Получена диаграмма перемещения выходного звена СО за один оборот входного (рис. 7). Визуализация и диаграмма перемещения звена СО позволяют легко иллюстрировать понятие «вы-стоя» или «останова» выходного звена.
Используя описанный выше поход к составлению уравнений, выполним моделирование реального рычажного меха-
Рис. 6. Визуальная модель рычажного механизма с выстоем
Рис. 7. Угол поворота выходного звена СБ в зависимости от угла поворота входного
низма игл машины ОВ-7, применяемого в текстильной промышленности. Кинематическая схема механизма представлена на рис. 8.
Конструктивные размеры и положение неподвижных шарниров 0, 01, 02 и 03 механизма (рис. 7) взяты из работы [3] .
Обозначим: а - угол поворота входного звена ОА, /3 - угол поворота треугольного звена Е02Б, у - угол поворота звена 01В, 8 - угол поворота звена 03Е, Л - угол при вершине 02 треугольного звена Е02Б.
Система координат проходит через центр О вращения кривошипа 0А, ось абсцисс совпадает с положением кривошипа ОА в момент мгновенной остановки выходного звена механизма.
Через координаты неподвижных точек записываем уравнения для шарниров Е, В, А. Соответственно:
\хР = хо - 03Р • 008(8);
1 Ур = Уо3 -0зр• 81П(8)-
\хв = Хо1 - 0, В • оо8(у);
1 Ув = Уо, -01в • 8шу).
\хА = х0 + 0А • оо8(а);
1 уА = х0 + 0А • 8т(а).
Для треугольного звена Е02Б через координаты неподвижного шарнира 02 вычислим координаты точек Е и Б:
\хв = х02 + 02 Е • 008(3); 1 Уе = У0 - 02Е • 81П(/).
Рис. 8. Механизм игл
\хв = х02 + 02 В • 008(3+Л); ¡Ув = У О, - 02 в • 81И(/ + Л).
Рис. 9. Модель механизма игл
Функциями «невязки» в этом случае будут длины звеньев АВ, АЕ, ВС, ОС и ЕС:
Ав = ЛхА -хв)2 +уА -Ув)2; АЕ = л/(У - хЕ )2 +УА - УЕ )2 ; вс =л1(хв - хс )2 +ув - Ус )2 ; пс = л](хп - хс {ув - Ус У ;
Рс =л1(хР - хс )2 +ур - Ус )2 .
Визуальная модель механизма представлена на рис. 9. Как и в предыдущем варианте, численное решение полученной системы уравнений позволяет построить диаграмму перемещения выходного звена.
Аналогичным образом получена модель восьмизвенного исполнительного механизма пресса (рис. 10) [4]. На рис. 11 показана визуальная модель пресса, закон изменения перемещения ползуна пресса и перемещение ползуна на этапе выстоя.
Рис. 10. Восьмизвенный механизм пресса
Рассмотренные примеры используются на практических и лабораторных занятиях. По окончании курса студентам предлагается самостоятельно промоделировать другие реальные механизмы.
Таким образом, применение пакета моделирования в рамках задачи кинематического анализа механизмов дает эффект быстрого начального вхождения в предмет «Моделирование», значительно активизирует процесс изучения кинематики работы многозвенных рычажных механизмов высоких классов и подвигает студентов к использованию возможностей пакета для решения других прикладных задач.
Литература
1. Карпов Ю. Имитационное моделирование систем. Введение в моделирование с AnyLogic 5. СПб.: БХВ-Петербург, 2005.
Рис 11. Визуальная модель и диаграммы работы пресса
2. Хомченко В.Г., Гебель Е.С., Солонин Е.В. Модификации и кинематический синтез рычажного механизма IV класса с выстоем выходного звена в крайнем положении по заданной циклограмме. Проблемы механики современных машин // Материалы 3-й международной конференции. Т. 1. Улан-Удэ: Изд. ВСГТУ, 2006. С. 7-10.
3. Гурвич П.И. Быстроходные основовязальные машины ОВ-7 и ОВ-8. М.: Легкая индустрия, 1970.
4. Власов В.И., Борзыкин А.Я., Букин-Батырев И.К. и др. Кривошипные кузнечно-прессовые машины. Под ред. В.И. Власова. М.: Машиностроение, 1982.
Abstract
In the article the application of Model Vision Studium software for simulation of kinematics of plane linkages within the subject Theory of Mechanics and Machines is taken into consideration. Some examples of plane linkage models are analysed.
Гебель Елена Сергеевна, старший преподаватель кафедры «Автоматизация и робототехника» Омского государственного технического университета, [email protected]
Солонин Евгений Владимирович, кандидат педагогических наук, доцент кафедры «Автоматизация и робототехника» Омского государственного технического университета, [email protected]
© Н а шм а вторы, 2008. Our authors, 2008