Лабораторная работа №2 моделирование колебательных систем
Цель работы:
Составить математические модели колебательных систем с одной степенью свободы и блок-схемы их набора в пакете Simulink для различных законов демпфирования.
Исследовать характер собственных и вынужденных колебаний систем под действием различных внешних сил и ненулевых начальных условий.
Теоретическая часть
Простейшая колебательная система (например, грузик на пружине, LC ‑ колебательный контур) имеет одну степень свободы и описывается обыкновенным дифференциальным уравнением первого порядка. Так простейшая одномассовая поступательная система описывается уравнением равновесия сил, действующих на подвижную массу m:
, (2.1)
где - перемещение массы;
- сила инерции;
- демпфирующая сила сопротивления движению, зависящая от скорости ;
- усилие пружины (восстанавливающая сила);
- внешняя сила.
В зависимости от характера демпфирующей силы уравнение (2.1) движения принимает следующие формы:
В случае вязкого трения, характерного для механических крутильных систем, демпфирующая сила является линейной функцией скорости:
, (2.2)
где - линейная сила вязкого трения.
В случае квадратичного закона трения, характерного для гидравлических систем, демпфирующая сила является квадратичной функцией скорости:
, (2.3)
где - квадратичная демпфирующая сила.
а б в
Рисунок 2.1. Колебательная система:
а. с вязким линейным трением,
б. с вязким квадратичным трением,
в. с сухим трением.
В случае сухого трения, характерного для механических поступательных систем, демпфирующая сила является функцией знака скорости, а величина ее постоянна:
, (2.4)
где - сила сухого трения.
На рис. 2.1 приведены примеры колебательных систем с вязким линейным, вязким квадратичным и сухим трением.
Свободные колебания системы возникают под действием ненулевых начальных условий (, ), либо при скачкообразном приложении силы F(t). Свободные колебания затухают под действием демпфирующей силы, а характер затухания зависит от вида демпфирования (трения).
Порядок составления структурной модели
Используя уравнения (2.2), (2.3), (2.4) построить блок-схему колебательной системы с различными демпфирующими силами. Для этого необходимо:
Разрешить дифференциальное уравнение относительно слагаемого со старшей производной.
С помощью цепочки сумматоров собрать модель правой части уравнения.
Используя цепочки из двух последовательно включенных интеграторов, получить переменные и .
С помощью усилителей с регулируемым коэффициентом передачи или соответствующих нелинейных блоков получить слагаемые правой части уравнения по п. 1 и подключить их к цепочке сумматоров.
Таблица 2.1
№ вар | m | b | c | № вар | m | b | c | ||||
1 | 10 | 0.5 | 10 | 10 | 1 | 13 | 8 | 0.48 | 32 | 10 | 2 |
2 | 20 | 1.2 | 0.8 | 20 | 3 | 14 | 22 | 1.32 | 1.98 | 20 | 3 |
3 | 15 | 1.2 | 60 | 30 | 2 | 15 | 24 | 1 | 24 | 30 | 1 |
4 | 25 | 1.25 | 0.25 | 10 | 1 | 16 | 30 | 2.7 | 1.2 | 10 | 2 |
5 | 12 | 1.08 | 60 | 20 | 3 | 17 | 10 | 0.7 | 40 | 20 | 3 |
6 | 40 | 2.8 | 0.4 | 30 | 2 | 18 | 20 | 1.6 | 0.2 | 30 | 1 |
7 | 15 | 1.35 | 15 | 10 | 2 | 19 | 36 | 1.8 | 1.44 | 10 | 2 |
8 | 4 | 0.28 | 36 | 20 | 1 | 20 | 16 | 0.8 | 64 | 20 | 3 |
9 | 20 | 1.4 | 80 | 30 | 3 | 21 | 50 | 4 | 4.5 | 30 | 1 |
10 | 18 | 1.8 | 18 | 10 | 3 | 22 | 25 | 2.5 | 100 | 10 | 1 |
11 | 10 | 0.6 | 0.9 | 20 | 2 | 23 | 32 | 1.6 | 0.32 | 20 | 3 |
12 | 25 | 1.75 | 25 | 30 | 1 | 24 | 14 | 1.12 | 56 | 30 | 2 |
Проверить блок-схему - в каждом замкнутом контуре число перемен знака должно быть нечетным.
Внешнюю силу F(t) моделировать соответствующим генератором.
В полученной схеме установить коэффициенты в соответствии с вариантом (табл. 2.1).
Набрать модель в рабочем окне SIMULINK.
Порядок выполнения работы
Составить и набрать структурную модель колебательной системы с вязким трением (по уравнению (2.2)). Задать:
возмущение - единичный скачок;
начальные условия - нулевые;
метод моделирования - Рунге-Кутта;
точность - 0.01;
длительность моделирования Tk = 50 с.
Запустить моделирование командой “Start” или “Счет” (функциональная клавиша F3), по окончании вычисления вывести на экран график процесса x(t) командой “Грф” (F6) и зарисовать.
Изменяя коэффициент демпфирования (уменьшая в 2, 4, 10 раз), получить графики соответствующих переходных процессов и зарисовать. Сделать вывод о влиянии коэффициента демпфирования на затухание колебаний.
Для заданного коэффициента демпфирования исследовать влияние коэффициента передачи во внешнем контуре на частоту собственных колебаний системы. Для этого увеличивать и уменьшать указанный коэффициент в 2 и 4 раза. Сделать выводы.
Исследовать влияние ненулевых начальных условий на движение системы. Для этого установить исходные значения коэффициентов, а сигнал на выходе генератора установить равным 0. Задать последовательно начальные условия и . Зарисовать графики переходных процессов при ненулевых начальных условиях. Сделать выводы.
Исследовать поведение системы с вязким трением при подключении в виде возмущающего воздействия синусоидального сигнала при нулевых начальных условиях в районе резонанса. Для этого необходимо найти резонансную частоту системы .
(2.5)
Определить амплитуду установившихся колебаний при частоте возмущающего сигнала 0.5·, , 2·. Сделать выводы.
Исследовать систему с вязким трением при подключении в виде возмущающего воздействия меандра (прямоугольных периодических колебаний) при нулевых начальных условиях. Зарисовать процесс x(t) при частоте возмущающего сигнала 0.25·, , 4·. Описать характер движений системы. Сделать выводы.
Исследовать колебательные системы с вязким линейным, с вязким квадратичным и с сухим трением (уравнения (2.2), (2.3), (2.4) соответственно) с нулевыми начальными условиями при единичном возмущающем воздействии. Зарисовать полученные выходные характеристики и сделать вывод о влиянии различных демпфирующих сил на вид переходной характеристики.
Контрольные вопросы
1. По каким законам уменьшается амплитуда колебаний при вязком линейном трении, при вязком квадратичном трении и при сухом трении?
2. Как изменяется амплитуда вынужденных колебаний в зависимости от частоты возмущающей силы?
3. Сколько собственных частот имеют колебательные системы, описываемые обыкновенными уравнениями 3-го, 4-го, 5-го, 6-го порядка?
4. Как изменяется затухание колебаний при изменении коэффициента демпфирования?
5. Как изменится собственная частота колебательной системы при изменении массы?
Литература
1. Пановко Я.Г. Введение в теорию колебаний. - М.: Наука, 1991.
Тимошенко С.П. Колебания в инженерном деле. - М.: Физматгиз, 1959.
- Лабораторная работа №1 введение в Simulink
- Краткие сведения о пакете
- Лабораторная работа №2 моделирование колебательных систем
- Лабораторная работа №3 моделирование нелинейных и дискретных систем
- Уравнение Ван дер Поля
- Уравнение Рэлея
- Бомбометание с малых высот с учетом сопротивления воздуха
- 10. Модель логической системы
- 11. Моделирование случайных событий
- Лабораторная работа №4 моделирование и оптимизация электромеханической системы привода прокатных валков
- Лабораторная работа №5 оценивание случайных параметров и регрессия
- Лабораторная работа №6 выявление скрытых периодичностей в случайном процессе
- Лабораторная работа №7 генерация случайных процессов с заданной спектральной плотностью
- Лабораторная работа № 8
- Идентификация динамических объектов
- По переходным функциям
- Теоретические сведения
- Зарегистрированный график изменения выходной величины при скачкообразном возмущении на входе – переходную функцию технологического объекта управления (тоу) можно использовать для его идентификации.
- Идентификация с помощью настраиваемой модели
- Лабораторная работа № 9 моделирование объектов с распределенными параметрами
- Решение
- 3. ВpdeToolboxможно задавать граничные условияДирихлеиНеймана. Зададим сначала граничные условия Неймана, которые задаются следующим образом:
- После этого на правой и на левой границах пластины зададим условия Дирихле.H– весовой коэффициент, аr– заданная температура.
- 7.Для улучшения качества отображения решения можно сделать некоторые настройки. Нажать и в открывшемся диалоговом окне указать:
- Задания для самостоятельного выполнения