Применение предобуславливания для ускорения вычислений в задачах динамической аэроупругости

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

Задачи, связанные с динамической аэроупругостью, возникают при проектировании различных объектов, таких как летательные аппараты и мосты. Часто их решение требует сложного численного моделирования, потребляющего большое количество вычислительных ресурсов. Эти задачи имеют решающее значение в аэрокосмической отрасли, поскольку они позволяют анализировать аэродинамические нагрузки на конструкции и предсказать их поведение [1–3]. Тем не менее высокая вычислительная сложность ограничивает их использование на практике. Поэтому актуально оптимизировать вычислительные процессы при решении этих задач.

Применение предобуславливания для ускорения вычислений в задачах...

Раскатова Марина Викторовна кандидат технических наук, доцент, Национальный исследовательский университет «МЭИ», Москва. Сфера научных интересов: анализ данных, разработка программных приложений. Автор более 100 опубликованных научных работ. ORCID: 0000-0001-7671-3312, SPIN-код: 80535041, AuthorID: 609945.

Одним из эффективных способов повышения скорости вычислений является использование метода предобуславливания [4–5]. Этот метод предназначен для ускорения процесса сходимости итерационных подходов к решению систем линейных уравнений. Основная концепция заключается в преобразовании исходной системы в эквивалентную, но с более благоприятным числовым условием, что делает ее более эффективной для итеративных техник.

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

Методика исследования

Изучим математическую модель динамической аэроупругости самолета [6], которая представлена в виде системы дифференциальных уравнений, записанных в матричной форме:

Mu’’ + Cu’ + Ku = F(u, u’, t) , (1) где M , C , K – матрицы, соответствующие массе, демпфированию и жесткости соответственно. Вектор обобщенных перемещений представим как u . Вектор аэродинамических нагрузок, который зависит от обобщенных перемещений u , их производной u’ и времени t , обозначим как F .

Для дискретизации рассматриваемой системы применяется метод конечных элементов, что приводит к образованию системы линейных алгебраических уравнений следующего вида:

A·x = b , (2) где A – матрица системы; x – вектор значений, которые необходимо найти; b – вектор свободных членов.

Для повышения производительности решения этой системы уравнений рекомендуется использовать итерационные методы с предобуславливанием. В качестве предобуславлива-теля рассматривается неполное LU-разложение матрицы A. Выбор данного типа предо-буславливателя определен его эффективностью и относительной простотой реализации.

Для оценки точности решения системы линейных уравнений (2) использовались следующие критерии.

• 1.    Относительная невязка:

• 2.    Относительная ошибка:

||Ax – b|| / ||b|| , (3) где ||·|| – евклидова норма; данный критерий показывает, насколько точно решение x удовлетворяет исходной системе уравнений (2).

|| x – x_ref| | / || x_ref ||, (4) где x_ref – референсное решение, полученное с высокой точностью; данный критерий характеризует отклонение полученного решения от точного.

Обоснование выбора метода предобуславливания

Среди рассмотренных методов предобуславливания наиболее предпочтительным является неполное LU -разложение ( ILU ). Этот выбор определен следующими преимуществами:

• 1)    высокая эффективность ускорения вычислений. Применение ILU -разложения в качестве предобуславливателя позволило сократить время решения системы (2) на 25…30 % по сравнению с решением без предобуславливания;

• 2)    относительная простота реализации. ILU -разложение является более простым в реализации по сравнению с другими методами предобуславливания, такими как алгебраический многосеточный метод;

• 3)    оптимальное сочетание скорости и точности. По результатам проведенных испытаний ILU-разложение предлагает удачный баланс между увеличением скорости вычислений и сохранением точности решения в отличие от более примитивного диагонального предобуславливания.

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

Результаты исследования и их обсуждение

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

В качестве тестового примера выбрана модель самолета, для которой была рассчитана динамическая устойчивость. Исследовались три варианта:

• 1.    Применение метода сопряженных градиентов без использования предобуслав-ливания.

• 2.    Неполное LU -разложение – использование метода сопряженных градиентов с пре-добуславливанием, основанным на неполном LU -разложении.

• 3.    ILU (0) – использование метода сопряженных градиентов с предобуславливанием, основанным на разложении ILU (0).

Результаты показали, что использование метода сопряженных градиентов с предобус-лавливанием позволило сократить время вычислений на 35…40 % по сравнению с решением системы без предобуславливания.

Применение предобуславливания для ускорения вычислений в задачах...

Применение неполного LU -разложения в качестве предобуславливателя дало прирост в 25…30 %. Метод ILU (0) также продемонстрировал положительный эффект, сократив время вычислений на 20…25 %.

Применение предобуславливания не только ускоряло вычисления, но и не приводило к значительному снижению точности решения. Эффективность предобуславливания зависит от характеристик матрицы системы (размер, степень заполнения) и параметров предобуславливания (например, степень заполнения в ILU ( k )).

Результаты времени вычислений и точности представлены в Таблице.

Таблица

Результаты времени вычислений и точности

Метод	Время вычислений, с	Относительная невязка	Относительная ошибка
Без предобуславливания	100	1e-06	1e-05
Неполное LU -разложение	65	1e-05	1e-04
ILU (0)	75	1e-05	1e-04

Источник: составлено авторами.

Перспективы дальнейших исследований

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

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

Заключение

Разработанная методика оптимизации вычислений при решении задачи динамической аэроупругости самолета, основанная на применении предобуславливания, продемонстрировала высокую эффективность. Использование неполного LU -разложения ( ILU ) в качестве предобуславливателя позволило сократить время вычислений на 25…30 %, сохраняя при этом приемлемую точность решения.

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

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

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