Table of Contents

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

Уникальная аэродинамическая сложность лопастей ротора

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

Трансонский поток на переднем лезвии

На наступающей стороне диска ротора наконечник лопасти движется в том же направлении, что и свободный поток, в результате чего локальные относительные числа Маха могут превышать 0,9. На этих скоростях на поверхности лопасти образуются сильные ударные волны, что приводит к резкому увеличению сопротивления (волнового сопротивления) и потенциальному разделению пограничного слоя. Это условие, известное как наступающая сжимаемость лопасти, ограничивает максимальную скорость наклона вертолета. CFD позволяет инженерам детально исследовать структуру удара, анализируя контуры давления и распределения чисел Маха для разработки форм лопастей, которые ослабляют или задерживают эти удары. Пропущенные кончики и тонкие секции аэродинамической пленки часто оцениваются и оптимизируются с помощью CFD для управления этим трансонным потоком без ущерба для производительности наведения.

Обратный поток и динамическая остановка на отступающем лезвии

На отступающих сторонах лопасть движется против свободного течения. Для поддержания сбалансированного момента качения отступающих лопастей необходимо оперировать под гораздо более высокими углами атаки. Возле корня отступающих лопастей скорость вращения ниже скорости перелета вперед, вызывая область обратного потока, где воздух перемещается от задней кромки к передней. Далее под бортом высокий угол атаки приводит к динамическому сваливанию. При этом явление большой вихрь формируется по передней кромке, растет по хорде и сбрасывается в будку. Этот вихрь вызывает временный скачок подъема, за которым следует внезапный и насильственный скачок. Динамический скачок накладывает сильные торсионные нагрузки и вибрации. Высокоточный CFD с использованием таких методов, как Detached Eddy Simulation (DES) является теперь стандартным методом прогнозирования наступления и интенсивности динамического стека, позволяющим конструкторам выбирать аэродинамические фольги и планформы, задерживающие сваливания и уменьшающие скачки нагрузки.

Blade-Vortex Interaction (BVI) и шумопотоки

Вертолетный шум, особенно при спуске или низкоскоростном маневрировании, преобладает во взаимодействии лезвия-Вортекса. В этом сценарии лезвие проходит непосредственно через задающий кончик вихря, генерируемый предыдущим лезвием. Это быстрое изменение локального угла атаки производит резкий импульс давления, который излучается как шум. Шум BVI является основной проблемой шума сообщества и ограничением на полетных операциях. Анализ CFD в сочетании с кодами акустической аналогии является основным инструментом для понимания BVI. Путем моделирования разворота, конвекции и траектории кончика кончика инженеры могут исследовать геометрии кончика лезвия (например, сметенные, суженные или ангедральные кончики), которые рассеивают вихрь или сдвигают его путь, чтобы уменьшить силу взаимодействия.

Вычислительные рамки для CFD Ротари-Винг

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

Перезагрузка сетей и движущиеся границы

Наиболее широко применяемый подход к CFD ротора — метод перезагрузки или сетки Химеры. Сетки, облегающие корпус, генерируются вокруг каждого отдельного лезвия и фюзеляжа. Эти сетки накладываются на более крупную фоновую декартовую сетку, которая покрывает вычислительную область. Сетки общаются посредством интерполяции на своих границах. По мере вращения или шага лезвий их сетки жестко перемещаются через фоновую сетку. Такой подход позволяет избежать сильного искажения сетки, которое могло бы возникнуть в результате деформации одной структурированной сетки. Генерирование высококачественных сеток перезагрузки для винтокрылых машин — это специализированный навык, требующий тщательного внимания к вырезанию дыр и разрешению зазора между компонентами сетки. Точность решения сильно зависит от качества и разрешения сетки в области лезвия.

Моделирование турбулентности для отдельных и вращающихся потоков

Точная модель турбулентности остается основным источником неопределенности в CFD лопасти ротора. Для прикрепленных потоков в наведении и умеренной скорости вперед модели Рейнольдса-Усредненного Навье-Стокса (RANS) такие как Spalart-Allmaras (SA) и Shear Stress Transport (SST) k-ω обеспечивают надежные прогнозы интегрированных нагрузок. Однако модели RANS плохо работают для массивного разделения и очень неустойчивых потоков характеристик динамических ситуационных и вихревых волновых взаимодействий. Для этих условий необходимы симуляции масштабирования и вихревых потоков. Для этих условий требуется симуляция с отложенным Эдди (DES), особенно вариант с отложенным DES (DDES), действует как RANS в прикрепленном пограничном слое и переключается в режим большого Эдди-симуляции (LES) в отдельных регионах. Этот гибридный подход обеспечивает разрешение, необходимое для захвата переходных вихревых структур при высокой вычислительной стоимости, но осуществимой для инженерного проектирования.

Численные схемы и временное разрешение

Потоки лопастей ротора требуют точных пространственных схем высокого порядка для минимизации численного рассеивания. Распределение искусственно рассеивает верховой вихрь, заставляя его слишком быстро ослабевать и приводя к недооценке нагрузок BVI и индуцированного сопротивления. Схемы с низким рассеиванием, такие как Weighted Essentially Non-Oscillatory (WENO) или методы восходящего ветра более высокого порядка, имеют решающее значение для сохранения силы вихря в течение нескольких оборотов ротора. Временно моделирование должно быть точным по времени, с использованием двойных по времени шагов или неявных схем для разрешения неустойчивого потока. Типичная революция ротора разделена на сотни или тысячи физических шагов, каждый из которых требует нескольких внутренних итераций для сходки нелинейных уравнений.

Методологии оптимизации для геометрии клинков

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

Проектирование параметризации пространства

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

Оптимизация на основе градиента и смежные реле

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

Суррогатная оптимизация и многоцелевой подход

Когда пространство проектирования является высоко нелинейным или когда информация о градиенте ненадежна, суррогатная оптимизация (SBO) является привлекательной альтернативой. SBO включает в себя построение модели поверхности отклика, такой как модель кригинга или модель функции радиальной основы, из набора точек проектирования, оцененных CFD. Суррогат дешев для оценки, позволяя оптимизатору быстро исследовать пространство проектирования. Подходы многоточности усиливают SBO, комбинируя несколько высокоточных оценок CFD (DES) со многими оценками низкой точности (RANS или методы вихревой решетки). Этот подход уравновешивает точность и вычислительные затраты, позволяя оптимизировать сложные явления, такие как динамический ларь или шум BVI в практические сроки.

Тематические исследования в CFD-управляемой оптимизации клинков

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

Оптимизация формы подсказки для акустических и эксплуатационных прибылей

Наконечник лопасти является наиболее важной областью для определения как производительности, так и шума. Сильно свернутый наконечник вихря увеличивает индуцированное сопротивление висению и приводит к шуму BVI в прямом полете. Оптимизация на основе CFD произвела расширенные формы наконечника, такие как наконечник BERP (Британская экспериментальная программа ротора) и наконечники с захлестнутыми застегнутыми анхедралями. Эти формы рассеивают процесс образования вихря, распространяя вихрь на большую площадь и уменьшая пиковые скорости вихря. Исследования показали, что оптимизированные наконечники могут уменьшить шум BVI на несколько децибел, одновременно улучшая показатель зависания на 1-3 пункта. Компромисс часто включает в себя повышенную сложность и стоимость изготовления, но прирост производительности, подтвержденный CFD, обеспечивает обоснование этих передовых конструкций.

Распределение фольги и планформа Twist

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

Проблемы с Rotor Blade CFD

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

Вычислительные затраты и время оборота

Высокоточные DES или LES-моделирования полного ротора в прямом полете чрезвычайно дорогостоящи с точки зрения вычислений, часто требующие дней или недель на больших высокопроизводительных вычислительных кластерах (HPC). Это ограничивает количество вариантов проектирования, которые могут быть оценены в типичной временной шкале промышленного проекта. Стоимость обусловлена необходимостью тонкого разрешения сетки (особенно в верхнем вихре и пограничном слое), небольшими временными шагами для акустического разрешения и итеративной конвергенцией, необходимой для точных нагрузок. Снижение стоимости этих симуляций за счет лучших решателей, ускорения GPU и эффективной выборки является основным направлением текущих исследований.

Численность диссипации и сохранение вихрей

Поддержание прочности и структуры верхнего вихря на нескольких оборотах ротора, возможно, является единственной величайшей технической проблемой для CFD вихря ротора. Численное рассеивание, присущее всем схемам дискретизации, имеет тенденцию размазывать ядро вихря и уменьшать его пиковую скорость вихря. Это приводит к неточному прогнозированию шума BVI и индуцированного сопротивления на лопастях нисходящего потока. В то время как схемы более высокого порядка помогают, они увеличивают вычислительную стоимость. Адаптивное сетчатое уточнение (AMR), где сетка динамически совершенствуется в области вихревого следа, предлагает мощное решение, но добавляет значительную сложность решателю и управлению сеткой.

Будущие направления и новые технологии

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

Машинное обучение для моделирования и проектирования турбулентности

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

GPU-ускоренные CFD и цифровые близнецы

Переход к массово параллельным решениям на основе GPU демократизирует высокоточные CFD. Решения, разработанные для архитектур GPU, могут достигать ускорений от 10 до 100 раз по сравнению с традиционными кластерами CPU для того же размера задачи. Это сокращение времени оборота позволяет использовать высокоточные CFD ранее в цикле проектирования и открывает дверь для создания цифровых двойников роторных систем. Цифровой двойник включает в себя виртуальную модель, которая постоянно обновляется с данными о нагрузке и производительности в полете, что позволяет осуществлять мониторинг в режиме реального времени, прогнозное обслуживание и оптимизацию здоровья и производительности лопастей в течение всего срока эксплуатации.

Многодисциплинарная оптимизация (MDO)

Будущее оптимизации лопастей заключается в тесной связи аэродинамики со структурной динамикой и акустикой. Аэроупругая пошивка, где структура лопасти предназначена для изгиба и скручивания таким образом, чтобы уменьшить нагрузки или улучшить производительность, требует одновременного анализа CFD и вычислительной структурной механики (CSM). Совместные фреймворки MDO, которые интегрируют CFD, CSM и акустические коды, позволят проектировать действительно интегрированные роторные системы, которые легче, тише и эффективнее, чем те, которые разработаны с помощью последовательных, однодисциплинарных процессов.

Заключение

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