Модель масса-пружина-демпфер

Динамическая система, включающая массу, пружину и демпфер, очень часто приводится для изучения природы колебательных процессов, кинематики и эффектов хаотического движения (если объектов несколько) или методов численного интегрирования.

Мы рассмотрим направленную модель msd_causal.engee, которая осуществляет численное интегрирование уравнения .

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

Вторая модель msd_acausal.engee создана из ненаправленных (акаузальных) блоков библиотеки физического моделирования Engee. В верхней части находится датчик, который измеряет перемещение массы относительно точки опоры, а также блок SolverConfiguration, который указывает глобальному решателю модели, что подсоединенные к нему сигналы и блоки будет иметь собственный интегратор с собственными настройками.

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

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

Обе модели используют идентичные переменные параметры: масса груза m, жесткость пружины k и коэффициент поглощения демпфера b.

Когда нужные переменные созданы, можно будет управлять моделями, меняя значения переменных в Окне переменных. При значениях b=1, k=2, m=3 мы должны наблюдать довольно медленный, инерционный процесс с заметным затуханием.

При помощи однострочной команды мы либо загрузим (load), либо откроем уже загруженную модель (open). Здесь мы пользуемся синтаксисом однострочных условных операторов x = условие ? если_истинно : если_ложно.

Посмотрим, какие параметры интегратора прописаны в обеих моделях.

Мы видим, что обе модели имеют одинаковые настройки интегратора на верхнем уровне. Интегратор физической модели задан в блоке Solver Configuration, на верхнем уровне опрос модели производится с постоянным шагом 0.01 с.

Выполним модели и сохраним полученные результаты в переменные rc и ra для направленной и физической моделей соответственно.

Сравнение результатов выполнения:

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

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

Теперь мы получим решение направленной модели методом с переменным шагом интегрирования и посмотрим, как изменится ошибка относительно физической модели.

Чтобы сравнить графики с разным разрешением нам потребуется интерполяция.

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

В Engee встроен Инспектор данных, который позволит сравнить разные результаты запусков между собой.

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

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

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