Мёртвая петля

Описание


Цель проекта: разработать программный продукт, представляющий собой виртуальную физическую лабораторию, моделирующую движение шарика в «мёртвой петле» при различных начальных условиях

Язык программирования: C++

Фреймворк: Qt4 / Qt5 / Qt6

Лицензия: GNU GPL v3

В программе моделируется поведение шарика, который совершает «Мёртвую петлю», то есть катится по круговому сектору. В модели можно изменять: угол наклона начального сектора (по которому шар набирает скорость), его длину, радиус проходимой шаром окружности и радиус шара.

В программе выводятся графики:


Скриншоты

 


Загрузка

Исходный код:
deathloop-1.2.zip
deathloop-1.2.tar.gz

Git:
https://github.com/AlienCowEatCake/deathloop.git

Бинарные сборки:
Компилятор / Qt ОС Ссылка
MinGW 3.4.2, Qt 4.4.3 Windows 98 SE и выше deathloop_qt4.4_mingw32.zip
MSVC 2017 x86, Qt 5.6.3 Windows XP SP3 и выше deathloop_qt5.6_msvc2017_x86.zip
MSVC 2017 x64, Qt 5.6.3 Windows XP x64 SP2 и выше deathloop_qt5.6_msvc2017_x64.zip
MSVC 2022 x86, Qt 5.15.5 Windows 7 x86 и выше deathloop_qt5.15_msvc2022_x86.zip
MSVC 2022 x64, Qt 5.15.5 Windows 7 x64 и выше deathloop_qt5.15_msvc2022_x64.zip
GCC 4.8.5 i486, Qt 5.6.3 GNU/Linux i486 deathloop_qt5.6_i486-linux-gnu.elf.gz
GCC 4.8.5 x86_64, Qt 5.6.3 GNU/Linux x86_64 deathloop_qt5.6_x86_64-linux-gnu.elf.gz
Clang 9.0.0 x86_64, Qt 5.15.5 GNU/Linux x86_64 (AppImage) deathloop_qt5.15_x86_64-linux-gnu.AppImage
Clang 9.0.0 x86_64, Qt 5.15.5 GNU/Linux x86_64 (deb) deathloop_1.2.0-1_amd64.deb
Apple Clang 13.1.6, Qt 5.6.3 Mac OS X 10.6 и выше (x86_64) deathloop_qt5.6_clang64_libstdcxx.dmg
Apple Clang 13.1.6, Qt 5.15.5 Mac OS X 10.10 и выше (x86_64/AArch64) deathloop_qt5.15_universal_target10.10.dmg

Архив версий:
Releases/


Список изменений

[v1.2] - 29 Aug 2022 Предыдущие версии
[v1.1] - 05 Dec 2021 [v1.0] - 31 Dec 2016 [v0.38] - 25 Aug 2015 [v0.37] - 05 Jan 2015 [v0.36] - 24 Aug 2014 [v0.35] - 20 Aug 2014 [v0.34] - 07 Aug 2014 [v0.33] - 09 Jan 2014 [v0.32] - 02 Dec 2013 [v0.31] - 02 Dec 2013 [v0.30] - 02 Dec 2013 [v0.29] - 01 Dec 2013 [v0.28] - 11 Nov 2013 Список предыдущих изменений, увы, не сохранился.


Модель

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

[image]

Траектория движения центра шарика при прохождении петли

β - угол наклона начального участка

L - длина начального участка петли

R - радиус петли

r - радиус шарика

α - текущий угол, с помощью которого задаём положение шарика

Движение шарика разбивается на четыре основных участка: движение по наклонной части желоба, движение по входному в петлю сектору β, движение в петле и выкатывание шарика по горизонтальному участку желоба.

Поскольку мы имеем дело с качением, для решения задачи нам необходимо составить уравнение динамики вращательного движения твёрдого тела относительно оси вращения.

M=I*dω/dt

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

[image]

h=r/sqrt(2)

Поперечное сечение шарика и желоба

Вращение шарика при качении происходит вокруг оси, обозначенной пунктиром.

Воспользуемся теоремой Штейнера.

I=(2/5)*m*r^2+m*h^2

Теперь рассмотрим первый участок движения шарика:

Относительно данной оси в уравнении динамики вращательного движения твёрдого тела ненулевым является только момент силы тяжести.

M_gr=m*g*h*sin(β)

Мы можем вычислить характеристики шарика:

I*ε=M_gr -> ε=M_gr/I

a=ε*h, x=x_0+g*t^2*sin(β)*cos(β)/3.6, y=y_0-g*t^2*sin(β)*sin(β)/3.6, v_x=g*t*sin(β)*cos(β)/1.8, v_y=g*t*sin(β)*sin(β)/1.8, v=sqrt(v_x^2+v_y^2)

На втором участке при расчёте углового ускорения шарика и мгновенного ускорения поступаем аналогично. Движение рассматриваем как движение по окружности (пунктир на рисунке соответствует траектории движения центра шарика).

Угол α изменяется от – β до 2π

x=(R-r*sqrt(2))*sin(α)+(L*cos(β)+R*sin(β)), y=-(R-r*sqrt(2))*cos(α)+R

a=(-5/9)*g*sin(α), ε_tr=a/(R-r*sqrt(2))

Так как α=f(ε(α),ω), найдём α при помощи метода Эйлера. Реализация его алгоритма находится непосредственно в коде.

v=ω*(R-r*sqrt(2))

Третий участок отличается от второго смещением по ещё одной оси координат (необходимо для движения именно по петле, а не по замкнутой окружности).

Движение по последнему участку:

x=L*cos(β)+R*sin(β)+v*t

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

v^2/(R-r/2)<g*cos(α+π)