Изогеометрический анализ (IGA): синтез систем САПР и методов конечных элементов

Конец пропасти между чертежом и расчетом

В современной инженерной практике существует болезненный исторический разрыв между дизайном и физическим расчетом. Инженеры-конструкторы создают идеальные, математически гладкие геометрические модели деталей (кузова автомобилей, корпуса кораблей) в системах компьютерного проектирования (САПР/CAD), используя аппарат рациональных B-сплайнов (NURBS). Однако, когда эти идеальные модели передаются инженерам-расчетчикам для симуляции напряжений или аэродинамики (САЕ), геометрическая красота уничтожается. Программы для метода конечных элементов (МКЭ) принудительно аппроксимируют эти плавные кривые сеткой из плоских треугольников или прямых тетраэдров.

Процесс генерации такой конечно-элементной сетки — это самая дорогая, трудоемкая и чреватая ошибками часть любого сложного анализа, отнимающая до 80% времени инженера. Любое изменение в CAD-модели требует полного перезапуска генератора сеток. Более того, граненые аппроксимации криволинейных поверхностей в МКЭ вносят значительные погрешности в расчеты поверхностных эффектов (например, при расчете акустического рассеяния или трения обшивки самолета). В 2005 году профессор Томас Хьюз предложил радикальную идею, призванную навсегда стереть эту границу — Изогеометрический анализ (Isogeometric Analysis, IGA).

Геометрия NURBS: идеальные круги и цилиндры

Фундаментальная идея изогеометрического анализа невероятно элегантна: почему бы нам не использовать для решения дифференциальных уравнений (поиска напряжений и деформаций) те же самые базисные функции, которые использовались для рисования самой геометрии детали в CAD-системе?

CAD-системы построены на базе сплайнов NURBS (Non-Uniform Rational B-Splines). В отличие от обычных полиномов, применяемых в МКЭ, NURBS являются рациональными функциями (дробями). Это дает им математическую суперспособность — они могут абсолютно точно, без малейшей аппроксимации описывать конические сечения (идеальные круги, эллипсы, сферы, цилиндры). Переход к IGA означает, что в методе Галеркина вместо обычных полиномов Лагранжа в качестве пробных функций используются функции NURBS. В результате расчетная сетка МКЭ просто исчезает как сущность. Уравнения механики сплошной среды решаются непосредственно на «контрольных точках» CAD-сплайна, гарантируя абсолютную, 100% геометрическую точность расчетной области без граненых артефактов.

Преимущества гладкости и математические вызовы IGA

Изогеометрический анализ произвел настоящую революцию в вычислительной механике. Базисные функции NURBS обладают высокой степенью непрерывности не только самой функции, но и ее производных сквозь границы элементов (что почти недостижимо в стандартном МКЭ). Это потрясающее свойство гладкости делает IGA идеальным инструментом для решения уравнений в частных производных высоких порядков, таких как уравнение Кана-Хилларда (моделирование фазового распада сплавов) или расчеты сверхтонких оболочек (Kirchhoff-Love shells), где требуются непрерывные вторые производные.

Однако внедрение IGA в массовую промышленность столкнулось с серьезными математическими препятствиями. CAD-системы исторически разрабатывались только для описания внешней поверхности (оболочки) детали. Для решения физической задачи (например, упругости) нам нужна сплошная объемная (3D) параметризация объекта сплайнами (Trivariate NURBS), что является сложнейшей задачей геометрического моделирования. Кроме того, тензорная структура сплайнов NURBS затрудняет локальное адаптивное измельчение сетки (нельзя просто добавить точку в одном месте, не протянув линию через всю деталь). Для решения этих проблем сегодня активно разрабатываются новые, более гибкие классы сплайнов: T-сплайны, LR-сплайны и иерархические B-сплайны, которые уверенно ведут нас к полной интеграции проектирования и анализа.