close

Вход

Забыли?

вход по аккаунту

?

3336.Исследование поведения динамики механической системы при помощи блочного построения структуры уравнений описывающих её движение в среде Matlab Simulink

код для вставкиСкачать
УДК 620.178.: 621.382
ИССЛЕДОВАНИЕ ПОВЕДЕНИЯ ДИНАМИКИ МЕХАНИЧЕСКОЙ СИСТЕМЫ ПРИ
ПОМОЩИ БЛОЧНОГО ПОСТРОЕНИЯ СТРУКТУРЫ УРАВНЕНИЙ,
ОПИСЫВАЮЩИХ ЕЁ ДВИЖЕНИЕ В СРЕДЕ MATLAB – SIMULINK
Юрий Артемьевич Можаев
Сибирская государственная геодезическая академия, 630108, г. Новосибирск, ул.
Плахотного, д. 10, старший преподаватель кафедры специальных устройств и технологий,
тел. 8-983-309-51-56, e-mail: yura82@mail.ru
В статье рассматривается задача динамики механической системы при помощи
блочного построения уравнений, описывающих еѐ движение, в среде Matlab - Simulink.
Ключевые слова: динамическая система, блочное построение, амплитудо - частотная
характеристика, программное движение.
RESEARCH OF THE BEHAVIOR OF DYNAMICS OF MECHANICAL SYSTEMS WITH
BLOCK BUILDING STRUCTURE EQUATIONS DESCRIBING ITS MOVEMENT IN THE
ENVIRONMENT MATLAB – SIMULINK
Yuriy A. Mozhaev
Siberian State Geodetic Academy, 10 Plahotnogo, Novosibirsk, 630108, senior teacher, department
of special devices and technologies, tel. 89833095156, e-mail: yura6810@mail.ru
In this paper we consider the dynamics of a mechanical system with block construction of the
equations describing its motion in the environment Matlab - Simulink.
Key words: dynamical system, building block, amplitude and frequency response, program
motion.
Для моделирования поведения различных динамических систем,
используются ЭВМ. Существует большое количество алгоритмических языков,
на которых может быть выполнено решение задачи. Выбор того или иного
языка программирования зависит от многих условий. Часто решающую роль
оказывает удобство программирования, наличие проверенных математических
методов, легкость представления результатов моделирования. Такими
особенностями обладает пакет Matlab [3, 4], содержащий в своем составе
инструмент визуального моделирования Simulink [1, 2 ].
Simulink - графическая среда имитационного моделирования, позволяющая
при помощи блок-диаграмм в виде направленных графов, строить
динамические модели, включая дискретные, непрерывные и гибридные,
нелинейные и разрывные системы. Simulink сочетает в себе наглядность
аналоговых машин и точность цифровых вычислительных машин, а также
обеспечивает пользователю доступ ко всем возможностям пакета Matlab, в том
числе к большой библиотеке численных методов. Подготовка задачи для
моделирования в Simulink проводится в следующей последовательности: 1 выбор расчетной схемы, 2 - составление системы уравнений описывающих
исследуемый процесс, 3 - определение начальных условий, 4 - составление
структурной схемы, 5 - моделирование возмущающих функций, 6 - определение
исходных данных, 7 - составление модели в среде Simulink, 8 - включение
средств визуализации, 8 – тестирование, 9 – решение, 10 - анализ результатов.
Пусть поведение механической системы описывается следующей системой
дифференциальных уравнений
mbl1

,
I0
  2 z1   12  2 y,
(1)
F t 
m
2 z1   12  ý .
M
M
Исследования будем проводить при следующих параметрах: m  9êã ,
b  0.245 ì , I0  0.997êã  ì 2 , 1  1350Ãö ,  z1  84.375c 1 , l1  0.4 ì , M  49.5êã
y


, hc  0.2 ì , y0  500 ì ñ2 .
Опишем структуру программирования в Simulink и то, как производится
реализация решений. Из уравнений (1) видно, что связанными являются только
два уравнения по  и y , а  получается из  . Поэтому наш анализ
исследований и структура построения системы уравнений значительно
упрощаются. При построении модели в Simulink за основу мы возьмѐм только два
параметра  и y .
Так как в системе уравнений фигурируют переменные y , y ,  ,  , то у нас
будут фигурировать 4 интегрирующих блока Integrator. Следует отметить, что y
и y получаются из y и y , а  и  получаются из  и  путѐм
интегрирования. Этим переменным соответствуют следующие блоки: Integrator
-  , Integrator1 -  , Integrator2 - y , Integrator3 - y . Далее нам необходимо
реализовать сначала правую часть уравнения y , так как она не зависит от
производной второго порядка. Для этого добавляем два блока усилителя,
которые выполняют умножение по константе: Gain усиливает входной сигнал
переменной  , превращая еѐ в сигнал 1.822 106  , Gain1 усиливает входной
сигнал  , превращая его в сигнал 168.75 . Далее два сигнала суммируются в
блоке Add, превращаясь в 168.75  1.822  106  . После этого сигнал
увеличивается
блоком
Gain2,
превращаясь
в
сигнал
0.182 168.75  1.822 106  . Далее блок Add1 суммирует три сигнала: 1 -


0.182  168.75  1.822 106   , 2 – задаваемая внешняя электродинамическая
сила в блоке Sine Wave Function, 3 – сигнал 0  y (так как в правой части уравнения
y не входит y , то его зануляем в Gain3). Далее суммарный сигнал подаѐтся на
вход блока Integrator 2 – это и будет выражение для y .
Теперь необходимо сформировать правую часть для  . Сигнал 2 y
получается при помощи усилителя Gain4, далее этот сигнал вместе с сигналом
168.75  1.822  106  складывается в блоке Sum, и результирующий сигнал
подаѐтся на вход блока Integrator.
Переменная  формируется при помощи усилителя Gain5, где сигнал 
умножается на константу mbl1 / I0  0.885 . Для отображения результатов
используется блок Scope, куда выводится необходимая нам переменная
изменяющаяся от времени.
Для построения АЧХ в Simulink используется подпрограмма LTI Viewer [5].
Инструмент Simulink LTI-Viewer входит в состав пакета прикладных программ
Control System Toolbox и предназначен для анализа линейных стационарных
систем. С помощью данного инструмента можно легко построить частотные
характеристики исследуемой системы, получить ее отклики на единичные
ступенчатое и импульсное воздействия, найти нули и полюса системы,
построить диаграмму Найквиста и Николса. Принцип измерения АЧХ такой: на
вход системы подается сигнал с заданной амплитудой, на выходе системы
снимается показание амплитуды; затем изменяется частота входного сигнала
(при этом амплитуда поддерживается такой, как была установлена в начале),
при этом снимаются показания амплитуды на выходе. Потом находится
отношение амплитуд на выходе к входной для каждой частоты, получается АЧХ
как безразмерная величина На рис.1 (б) изображены графики двух
безразмерных амплитуд ускорений. Aout - снимаемая амплитуда на выходе
системы (в нашем случае амплитуды  и y ), Ain - снимаемая амплитуда на
входе системы (амплитуда электродинамической силы F   ).
Как видно из графиков рис. 1 на частоте 1350 Гц наблюдается режим
динамического гашения параметра y , а на частоте 1750 Гц наблюдается
максимум амплитуды параметров y и  .
Далее проведѐм исследования колебаний модели перевѐрнутого
физического маятника с четырьмя попарно равноудалѐнными осцилляторами от
оси симметрии и точки подвеса подвижной части электродинамического
возбудителя рис. 2. Где c - коэффициент, характеризующий демпфирование
при угловом перемещении; cx, y,z - коэффициенты демпфирования при
перемещении вдоль осей x ,y и z; k - коэффициент угловой жѐсткости
пружины;
k x, y, z - коэффициенты жѐсткости при линейных перемещениях
вдоль осей y и z; M , m1234 массы маятника и осцилляторов; I 0 - момент
инерции физического маятника относительно оси, перпендикулярной плоскости
рисунка и проходящей через точку О. При этом m1  m2 , m3  m4 , cz1  cz2 ,
cz3  cz3 , k z1  k z2 , k z3  k z3 , l1  l2  l3  l4 , z2 b1  b2 , b3  b4 .
Рис. 1. К решению системы дифференциальных уравнений: а) – структурная
схема в Simulink; б) – зависимости параметров y и  от частоты
Рис. 2. Модель перевѐрнутого физического маятника
Уравнения, описывающие колебательные движения системы в линейном
приближении и учѐтом того, что парциальные частоты осцилляторов лежат
значительно выше парциальной частоты угловых колебаний маятника,
имеют вид
mbl1
mbl2

1 
2 ,
I0
I0
1  2 z1 1  121  2 y,
 2  2  z 2  2  2 2  2   2 y
(2)
1  z1  z2 ,
3  z1  z2 ,
y
F t 
m
m1
2 z1   12  3 2 z3   32  ý .
M
M
M




Структурная схема, составленная в пакете Simulink рис. 3 составлялась с
учѐтом следующих физических параметров исследуемой системы:
m1  m2  9êã, m3  m4  7êã, M  45êã, ñz1  cz2  175c 1 , ñz3  cz4  93.75c1 ,
1  2800 Ãö, 2  1500 Ãö, l1  l2  0.25, l3  l4  0.5 .
Как видно из аналитических исследований рис.4, в районе двух
парциальных частот осцилляторов 1500 и 2800 Гц наблюдается искажение
вибрационного режима в рабочем направлении в виде гашения колебания, а в
поперечном направлении в районе тех же частот наблюдается значительное
усиление колебаний.
Аналитические
исследования
проводились
при
помощи
дифференциальных уравнений, которые описывают колебательные движения
механической системы, в линейном приближении без учѐта нелинейных членов
при помощи пакета Simulink путѐм блочного программирования. Построены
амплитудо-частотные характеристики (АЧХ) в случае переменной амплитуды
возбуждающей силы, зависящей от частоты. Показано, что в зоне парциальных
резонансов наблюдается динамическое гашение задаваемого программного
движения в вертикальном направлении и появляются интенсивные поперечные
колебания в направлении ортогональному программному.
Рис. 3 Структурная схема, составленная в пакете Simulink для исследования
системы дифференциальных уравнений, описывающих колебания маятника с
четырьмя осцилляторами
12
10
8
6
4
2
0
1000
2000
3000
4000
5000
6000
7000
8000
9000
10000
Рис. 4. Зависимости уровня ускорений от частоты для физического
маятника взаимодействующей с осцилляторами; - в вертикальном
направлении,
- в поперечном направлении
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Герман-Галкин С.Г. MATLAB & Simulink. Проектирование мехатронных систем на
ПК.- СПб.: КОРОНА – Век, 2008.
2. Гудвин Г.К., Гребе С.Ф., Сальгадо М.Э. Проектирование систем управления // М.:
Бином. Лаборатория знаний, 2004.
3. Добонравов В.В.Основы аналитической механики. Учебное пособие для вузов М.,
«Высшая школа», 1976. 264 c.
4. Дьяконов В., Круглов В. Математические пакеты расширения MATLAB. СПб.:
Питер, 2001.
5. Черных И. В. SIMULINK- среда создания инженерных приложений / Под общ. ред.
к.т.н. В. Г. Потемкина. - М.: ДИАЛОГ-МИФИ, 2003.
© Ю.А. Можаев, 2012
1/--страниц
Пожаловаться на содержимое документа