Сработал бы руль высоты на скорости 0,9 Маха? Вычислительная гидродинамика (CFD) в анализе крушен...
CFD Toolbox
0:00 / 0:00
Сработал бы руль высоты на скорости 0,9 Маха? Вычислительная гидродинамика (CFD) в анализе крушен...
25 578 просмотров · 12 ч назад
CFD Toolbox
479 подписчиков
25 578 просмотров · 12 ч назад
3 ноября NTSB сообщило, что обнаружило большую часть левого руля высоты самолета Learjet 60, потерпевшего крушение в округе Палм-Бич: три фрагмента, расположенные на расстоянии от 1,3 до 3,8 миль от кратера. Хвостовая часть разваливалась еще до удара. Расследование NTSB продолжается, причина крушения пока не установлена; в этом видео она также не названа.
Это канал, посвященный вычислительной гидродинамике (CFD), поэтому я использовал физику для ответа на четыре вопроса: где именно на траектории полета оторвался руль высоты, насколько сильно нагружался фюзеляж, работает ли руль высоты на скорости 0,9 Маха и может ли обледенение вызвать вибрацию хвостовой части. Весь анализ основан на моих собственных оценках, полученных из общедоступных данных, а модели CFD, FEA и флаттера используют ОБЩУЮ модель хвостовой части, а не данные Learjet.
Что я обнаружил (мои оценки):
Дрейф обломков (1500 сбросов, слабый ветер 10–20 узлов): фрагмент длиной 3,8 мили помещается в сброс на высоте около 26 000 футов, примерно через 37–41 секунду после начала снижения. Это до достижения расчетной скорости пикирования. Подъем на высоте 8000 футов приводит к тому, что все обломки оказываются в пределах ~1,5 мили, поэтому они не помещаются в дальний фрагмент.
Траектория транспондера: почти нулевая перегрузка в верхней точке (самолет начал пикирование), после Mmo 0,81 на +18 с, Mach 0,87–0,90 в течение последних 45 с, и динамическое давление выше расчетной скорости пикирования на высоте 22 000 футов, достигающее 1,9-кратного предела к последней точке.
SU2 RANS (k-ω SST) на типовом участке 9 %: первые 4° подъемной силы руля высоты теряют около 85 % своего эффекта в диапазоне чисел Маха от 0,80 до 0,90. При 8° подъемная сила восстанавливается, примерно на 30 % ниже числа Маха 0,80. Это мертвая зона.
Невязкий решатель этого не учитывает. Потеря происходит из-за отрыва потока, вызванного ударной волной.
Смещение центра тяжести из-за изменения числа Маха: центр давления смещается с 0,26 c до 0,39 c в диапазоне чисел Маха от 0,5 до 0,87.
Типичный FEA-анализ пластины: тот же входной параметр, который дает предельную нагрузку при расчетной скорости пикирования, дает около 97 % предельной нагрузки при зарегистрированном динамическом давлении.
Модель флаттера с 2 степенями свободы: около 1 мм льда на задней кромке снижает скорость флаттера ниже 300 м/с. Механизм подходит; образование льда при −48 °C — нет.
Предыдущее видео (четыре аэродинамических утверждения): [ССЫЛКА НА ВИДЕО_1]
Временные метки:
0:00 Отвалился хвост
0:57 Что заявило NTSB
1:38 Что говорит транспондер
2:25 Где отвалился руль высоты?
3:48 Одна временная шкала
4:25 Настройка CFD
5:09 Руль высоты перестал работать
6:29 Почему теория невязкого движения не верна
7:06 Сход с траектории Маха
7:39 Что несло хвост
8:33 Обледенение и флаттер
9:31 Доска
Данные и метод:
Траектория и скорости: гранулированный файл ADS-B Flightradar24 (неофициальный), сглаживание 6 секунд. Нагрузки усредненные, поэтому короткие пики невидимы.
Ветер и температура: радиозонд KMFL (Майами), 12Z 1 ноября 2026 г., данные предоставлены Iowa Environmental Mesonet.
Обломки: модель дрейфа Монте-Карло, 1500 выборок с изменяющимся сопротивлением, ветром и скоростью выброса. Используются только расстояния до обнаруженных обломков; их направление не разглашается.
CFD: SU2 v8, сжимаемый RANS, k-ω SST. Двумерная O-сетка из 320 × 150 ячеек, первая ячейка 2e-6 хорда, y+ ниже 1,2, Re 1,4e7. Типичный симметричный участок 9 % с элероном на кормовых 30 %. Число Маха 0,70–0,93, элерон под углом 0°, 4° и 8°.
Сравнение в невязком состоянии и эффект Маха: трансзвуковой решатель для малых возмущений (схема Мурмана и Коула, 1971 г.), угол атаки 1°.
Конечно-элементный анализ: модель пластины Кирхгофа (элементы Морли) типового конического алюминиевого стабилизатора, закрепленного на киле, откалиброванная таким образом, чтобы предельная нагрузка давала значение, равное предельному значению / 1,5. Это масштабное изображение, а не утверждение о прочности.
Флаттер: двухстепенная модель с использованием аэродинамики Теодорссена (метод p-k), несжимаемая. Крутящий момент настроен таким образом, чтобы хвостовое оперение флаттерировало при 1,2 × V_D.
Источники:
Брифинги NTSB для СМИ, 2 и 3 ноября 2026 г. (по сообщению AP, WPTV, CBS12)
Flightradar24
Murman & Cole, AIAA J. 9(1), 1971
Theodorsen, NACA Report 496 (1935)
14 CFR Part 25 (предельные и предельные нагрузки, запас по флаттеру)
Источники изображений:
Learjet 60 N208BH (не тот самолет, который попал в аварию): Wikimedia Commons, CC BY-SA 4.0. Вырезано и адаптировано.
#Learjet #CFD #Аэродинамика
---
Отказ от ответственности: Это видео предназначено только для образовательных и информационных целей и не является профессиональной инженерной консультацией или анализом аварии. Расследование NTSB продолжается, и причина еще не установлена. Весь анализ является моей собственной оценкой на основе общедоступных данных по типовому участк...