Министерство образования и науки Российской Федерации

Федеральное агентство по образованию

Саратовский государственный технический университет

Балаковский институт техники, технологии и управления

Анализ динамических систем

методом математического моделирования

Методические указания

по выполнению лабораторной работы

для студентов специальности 120100

всех форм обучения

Одобрено

редакционно-издательским советом

Балаковского института техники

технологии и управления

Балаково 2008

Введение

При исследовании динамических систем они описываются дифференциальными уравнениями. Для решения даже относительно простых дифференциальных уравнений требуется большое количество вычислений. Поэтому практическое исследование динамических систем невозможно без использования вычислительной техники. В данной работе на примере дифференциального уравнения второго порядка рассматриваются решения дифференциальных уравнений методом математического моделирования в системе SimulinkMatLab.

1. Построение математической модели динамической системы

На механическую систему, содержащую пружину жесткостью k, демпфер с коэффициентом вязкостиb, тело массойM(рис. 1) действует силаF(t).Необходимо построить математическую модель, описывающую движение системы во времени.

Входной переменной является сила, действующая на систему .

За выходную переменную примем положение подвижной массы М относительно базовой поверхностиx(t). Структурная схема системы приведена на рис.2.

Составим уравнение движения на основе второго закона Ньютона, согласно которому момент инерции системы равен сумме сил, действующих на систему.

Момент инерции системы равен .

Сила упругости пружины на основании закона Гука пропорциональна деформации пружины .

Сила сопротивления движению поршня демпфера на основании закона Ньютона для идеальной жидкости пропорциональна скорости движения поршня демпфера .

Тогда уравнение равновесия Ньютона

.

Приведем уравнение к стандартному виду, когда все выходные координаты находятся в левой части, входные – в правой части, а коэффициент при выходной переменной равен единице

.

Мы получили дифференциальное уравнение второго порядка.

2. Решение дифференциального уравнения методом математического моделирования в системе Simulink

Разрешим дифференциальное уравнение относительно старшей производной .

Считая, что вторая производная выходной переменной известна, проинтегрируем ее два раза и получим.

Рассчитаем вторую производную выходной переменной путем суммирования элементов правой части уравнения. Подав выход сумматора на вход первого интегратора, получим структурную схему системы, решающую данное дифференциальное уравнения (рис. 3).

Особенностью данной схемы является необходимость для нахождения второй производной выходной переменной знания этой второй производной. Это условие компенсируется заданием начальных значений первой производной и самой выходной переменной, которые называются начальными условиями .

В общем случае решение дифференциального уравнения является суммой двух составляющих . В зависимости от начального состояния системы и наличия входного воздействия могут быть различные случаи движения системы.

1. Свободное движение системы определяет движение системы из заданного начального состояния () в установившееся состояние при отсутствии входного воздействия и описывается свободной составляющей.

Пример 1. Массу М сместили на какое-то расстояние и отпустили. Под действием пружины система придет в движение. Движение системы в данном случае описывается свободной составляющей при ненулевом начальном отклонениии нулевой начальной скорости.

Пример 2. Массу сместили на какое-то начальное отклонение и толкнули с постоянной скоростью. Движение системы в данном случае описывается свободной составляющей при заданном начальном отклонении и заданной начальной скорости.

2. Вынужденное движение системы из начального состояния (с нулевыми начальными условиями). При воздействии на систему, находящуюся в начальном состоянии внешней силы F(t)движение системы описывается вынужденной составляющей. При этом закон движения системы зависит от зависимости изменения внешней силы во времени.

Пример 1. На систему вниз стала действовать постоянная сила . Движение массы М вниз до уравновешивания действующей силысилой упругости пружиныописывается вынужденной составляющей.

Пример 2. На систему действует синусоидальное входное воздействие . Система из состояния покоя перейдет в новое установившееся состояние, которым будет синусоидальное движение системы с круговой частотойи амплитудой, зависимой от амплитуды входного воздействия и параметров системы.

3. Вынужденное движение системы с ненулевыми начальными условиями. Данное движение описывается суммой свободного и вынужденного движений и зависит от начальных условий и закона изменения во времени входного воздействия.