Продолжение поста «Репортаж из сортира ЭП2Д, стандарт+»
После того, как я написала данный пост, мне подумалось: какова судьба пустельги, залетевшей в туалет электричке, и зависшей над унитазом в воздухе? Не будет ли ее силами инерции поезда относить к стенкам? Понятное дело, что воздушные потоки здесь решающие. Если окно в туалете закрыто, то вуаля! Все норм. Но мы все же хотим рассчитать влияние рывков поезда точно? Поэтому я при помощи ИИ и родного С++ сделала эти расчеты. И так:
Физико-математическая постановка задачи для C++ симуляции
1. Ввод системы координат и кинематика поезда
Вводим неинерциальную систему координат (X, Y, Z), жестко связанную с унитазом вагона ЭП2Д.
Ось X — направлена продольно по ходу движения поезда.
Ось Y — направлена поперечно (в левую стенку кабины).
Ось Z — вертикально вверх.
Зададим траекторию поезда на плоскости земли с учетом крутого поворота радиуса R и постоянного продольного рывка
Угловая скорость вращения поезда на повороте:
Тогда вектор мгновенного ускорения унитаза
на плоскости земли состоит из тангенциального (из-за рывка) и центростремительного (из-за поворота) ускорений:
2. Силы, действующие на пустельгу относительно унитаза
В нашей системе координат на пустельгу массой m действуют следующие силы:
Сила тяжести:
2. Подъемная сила крыльев:
(в режиме идеального зависания на месте пустельга полностью компенсирует гравитацию:
).
3. Переносная сила инерции (главный враг птицы):
Швыряет птицу продольно при рывках и размазывает по стенке центробежной силой на поворотах.
4. Сила Кориолиса: Возникает, если птица пытается скорректировать полет и начинает двигаться относительно унитаза со скоростью
во время поворота поезда:
5. Аэродинамическое сопротивление увлекаемого воздуха:
Поскольку окно закрыто, воздух движется с поездом. Но при изменении скорости поезда возникает сила лобового сопротивления, пропорциональная относительной скорости птицы и плотности воздуха
3. Дифференциальные уравнения движения птицы для C++ (double x, y, z;)
Чтобы пустельга оставалась строго над унитазом, её результирующее ускорение должно быть равно нулю. Но биологический отклик птицы (время реакции нейронной сети мозга на взмах крыла) имеет задержку
Уравнения движения системы для численного интегрирования (метод Эйлера или Рунге-Кутты в коде):
Где
— это судорожные попытки пустельги сманеврировать крыльями, чтобы вернуться в координаты (0,0).
4. Граничные условия баг-трекера (Критерий чистоты перьев)
Кабина туалета имеет размеры: длина L, ширина W. Условие успешного завершения потока симуляции:
Если координата выходит за рамки — std::cout << "CRASH: Пустельга стукнулась в стенку вагона!" << std::endl;
Значения всех констант для кода
Параметры пустельги:
m (масса) = 0.2 кг (средний вес взрослой птицы).
S (площадь крыльев) = 0.06 м² (при размахе около 70 см).
Cx, Cy (коэффициенты лобового сопротивления) = 0.4 (в режиме зависания птица уязвима продольным потокам).
tau (время реакции нервной системы птицы) = 0.1 секунды (время на взмах-корректировку).
Параметры поезда (ЭП2Д) и геометрии:
W (ширина туалета) = 1.2 м.
L (длина туалета) = 1.4 м.
v_p (начальная скорость поезда) = 22.2 м/с (около 80 км/ч).
j_x (рывок при экстренном торможении) = -1.5 м/с³.
R (радиус крутого поворота) = 500 метров.
rho (плотность воздуха внутри) = 1.2 кг/м³.
Код программы pustelga
В итоге логи программы:
ЗАКЛЮЧЕНИЕ
В результате проведенного компьютерного моделирования динамики полета зависающей пустельги в замкнутом неинерциальном объеме сортира подвижного состава ЭП2Д были сделаны следующие выводы:
Устойчивость к продольным возмущениям (Ось X):
Пустельга (с задержкой отклика tau = 0.1 с) успешно справляется с линейным продольным торможением или разгоном и сопутствующими рывками поезда. Птица совершает затухающие колебания вдоль продольной оси X, не достигая критических барьерных зон. Максимальное приближение к фронтальной стене составило 0.568 м на временной отметке 2.8 с при лимите в 0.7 м.Потеря устойчивости в поперечной плоскости (Ось Y) — Главный деструктивный фактор:
Критический сбой системы стабилизации происходит на отметке 4.0 с симуляции из-за эффекта центробежной силы и силы Кориолиса при повороте поезда по радиусу R = 500 м.
Поскольку поперечные воздушные массы внутри заблокированного контура кабины неподвижны относительно унитаза, аэродинамическое сопротивление воздуха оказывается недостаточным для компенсации нарастающего центробежного ускорения. Из-за физиологического запаздывания нейронного ответа tau птица не успевает заложить адекватный контр-вираж, в результате чего происходит аварийный бум крыла о боковую стенку вагона на координате Y = -0.610 м (при критическом допуске 0.6 м).
Заключение экспертной комиссии
Мойте после себя унитазы, суки!!!