лаб. №3 - форма 51 (1087400), страница 2
Текст из файла (страница 2)
В рассматриваемом случае эта функция приобретает такой вид:
причем безразмерные величины ζ и υ определяются следующими выражениями:
Такая безразмерная форма представления уравнений является предпочтительной, поскольку позволяет сократить количество параметров, а также получить решение в более общем виде. Вынесение вычисления моментов внешних действий в отдельную процедуру также позволяет сделать задачу нахождения решения более общей.
Исходными (задаваемыми) параметрами моделирования будем считать:
-
параметры самого маятника, к ним в данном случае относится коэффициент затухания ζ;
-
параметры, характеризующие внешнее воздействие, сюда входят: амплитуды виброперегрузок в вертикальном nmy и горизонтальном nmx направлениях; относительная (по отношению к частоте собственных малых колебаний маятника) частота вибрации опоры υ; начальные фазы εу, εх вибрации точки опоры;
-
начальные условия движения маятника, то есть начальное отклонение φ от вертикали и начальная безразмерная угловая скорость маятника φ'.
К выходным (моделируемым) величинам будем относить текущий угол отклонения маятника от вертикали φ(τ) и его безразмерную угловую скорость φ'(τ).
Запишем уравнение (2) в удобной для интегрирования форме:
5.2.2.2. В среде Simulink создадим подпрограмму, позволяющую по известному значению ускорения определять значения безразмерной угловой скорости и угла отклонения маятника. Схема интегрирования выглядит следующим образом:
Рисунок 3.2. Схема подсистемы интегрирования дифференциального уравнения движения маятника
При настройке первого блока Integrator в качестве значения его параметра Initial condition указать fit0, а второго fi0.
5.2.2.3. В среде Simulink создадим подпрограмму, позволяющую определять значение выражения 6. назовем ее внешние моменты сил. Схема подпрограммы представлена на следующем рисунке.
Рисунок 3.3. Схема подсистемы определения внешних моментов сил
Величина dz (относительный коэффициент затухания) использована в блоке Gain. При формировании перегрузок по горизонтальной и вертикальной осям в параметрах настройки использованы переменные nmx, nmy, nu, ex, ey. Пример параметров настройки блока вибрации приведен ниже.
Рисунок 3.4. Параметры настройки блока вибрации
5.2.2.4. На основании разработанных подпрограмм соберем основную схему, реализующую выражение 9. Схема приведена ниже.
Рисунок 3.5 Общая схема модели маятника
Перед началом интегрирования зададим параметры численного интегрирования уравнения маятника.
Рисунок 3.6. Параметры численного интегрирования.
5.2.2.5.Обеспечим присвоение значений переменных в параметрах блоков. Для этого создадим m-файл и напишем программу следующего содержания.
Рисунок 3.6. Программа присвоения значений переменным
Получим графики колебаний маятника и фазовый портрет колебаний.
5.2.2.6. Изменить входные параметры модели в соответствии с вариантами, предложенными в таблице 1.
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | |
nmx | 0 | 0.2 | 0 | 1.1 | 1 | 0.7 | 1.3 | 0.4 | 0.1 | 0 |
nmy | 1 | 1.1 | 1.2 | 0 | 0.3 | 0.6 | 0 | 0.9 | 1.2 | 1.3 |
ey | 130 | 120 | 110 | 100 | 90 | 80 | 70 | 60 | 50 | 30 |
ex | 10 | 20 | 30 | 40 | 50 | 60 | 70 | 80 | 90 | 100 |
nu | 2,3 | 2,4 | 2,5 | 2,6 | 2,7 | 2,8 | 2,9 | 3 | 3,1 | 3,2 |
dz | 0,1 | 0,2 | 0,05 | 0,06 | 0,08 | 0,09 | 0,11 | 0,12 | 0,13 | 0,14 |
fi0 | 0 | 5 | 0 | 6 | 0 | 8 | 0 | 9 | 0 | 2 |
fit0 | 1 | 0 | 12 | 0 | 11 | 0 | 3 | 1 | 6 | 2 |
11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 | 19 | 20 | |
nmx | 1.2 | 0.4 | 0.5 | 0.7 | 1.2 | 1.1 | 0.9 | 0.7 | 1.5 | 1 |
nmy | 0 | 0.6 | 0.8 | 0.5 | 0.1 | 0.3 | 0.2 | 0.7 | 0.6 | 1.2 |
ey | 20 | 10 | 20 | 30 | 40 | 50 | 60 | 70 | 80 | 90 |
ex | 120 | 130 | 110 | 100 | 90 | 80 | 70 | 60 | 40 | 10 |
nu | 3,3 | 4 | 3.4 | 3.2 | 3 | 2.5 | 2 | 1.5 | 1 | 0.5 |
dz | 0,15 | 0,17 | 0.19 | 0.16 | 0.14 | 0.12 | 0.13 | 0.17 | 0.2 | 0.21 |
fi0 | 0 | 4 | 0 | 5 | 0 | 7 | 0 | 9 | 0 | 11 |
fit0 | 8 | 1 | 6 | 0 | 4 | 0 | 2 | 0 | 1 | 0 |
Получить графики колебаний маятника и фазовый портрет колебаний.
5.2.3. Методические указания к моделированию и рекомендации к содержанию отчета
5.2.3.1. Кратко описать принципы функционирования разделов программных пакетов Matlab и Simulink.
5.2.3.2. Привести распечатки, подтверждающие выполнение всех пунктов экспериментальной части и текстовые пояснения к ним. Записать дифференциальные уравнения движения маятника, получить и описать его колебания под действием гармонической возбуждающей силы.
5.2.3.3. Привести описание полученных подсистем, их схемы. Описать входные и выходные параметры.
Отчет должен содержать следующие разделы:
1. Цель работы.
2. Программа работы.
3. Выводы. Сравнительный анализ результатов моделирования.