Метод конечных элементов

Численный метод решения физических или инженерных задач

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

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

МКЭ — это общий численный метод решения уравнений в частных производных с двумя или тремя пространственными переменными (т. е. некоторых краевых задач ). Чтобы решить проблему, МКЭ делит большую систему на более мелкие и простые части, называемые конечными элементами . Это достигается за счет особой дискретизации пространства по размерностям, которая реализуется путем построения сетки объекта — числовой области решения, имеющей конечное число точек. Формулировка краевой задачи методом конечных элементов в конечном итоге приводит к системе алгебраических уравнений . Метод аппроксимирует неизвестную функцию в области определения. ^[1] Простые уравнения, моделирующие эти конечные элементы, затем собираются в более крупную систему уравнений, моделирующую всю проблему. Затем FEM аппроксимирует решение, минимизируя соответствующую функцию ошибок с помощью вариационного исчисления .

Изучение или анализ явления с помощью FEM часто называют анализом конечных элементов ( FEA ).

Базовые концепты

Разделение всей области на более простые части имеет ряд преимуществ: ^[2]

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

Типичная реализация метода включает в себя:

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

Глобальная система уравнений имеет известные методы решения и может быть рассчитана на основе начальных значений исходной задачи для получения численного ответа.

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

• набор алгебраических уравнений для стационарных задач,
• набор обыкновенных дифференциальных уравнений для нестационарных задач.

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

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

Практическое применение FEM известно как анализ конечных элементов (FEA). FEA, применяемый в технике , представляет собой вычислительный инструмент для выполнения инженерного анализа . Он включает в себя использование методов создания сетки для разделения сложной задачи на мелкие элементы, а также использование программного обеспечения, закодированного с помощью алгоритма FEM. При применении FEA сложная проблема обычно представляет собой физическую систему с базовой физикой , такой как уравнение балки Эйлера-Бернулли , уравнение теплопроводности или уравнения Навье-Стокса, выраженные либо в PDE, либо в интегральных уравнениях , в то время как разделенные малые элементы Сложная задача представляет собой различные области физической системы.

FEA может использоваться для анализа задач в сложных областях (например, автомобили и нефтепроводы), когда область изменяется (например, во время твердотельной реакции с движущейся границей), когда желаемая точность варьируется по всей области или когда решение отсутствует. гладкость. Моделирование FEA представляет собой ценный ресурс, поскольку исключает необходимость многократного создания и тестирования сложных прототипов для различных ситуаций с высокой точностью. ^{[ нужна цитата ]} Например, при моделировании лобового столкновения можно повысить точность прогнозирования в «важных» областях, таких как передняя часть автомобиля, и уменьшить ее в задней части (таким образом снижая стоимость моделирования). Другим примером может служить численный прогноз погоды , где более важно иметь точные прогнозы развития сильно нелинейных явлений (таких как тропические циклоны в атмосфере или водовороты в океане), а не относительно спокойных областей.

Четкое, подробное и практическое изложение этого подхода можно найти в книге « Метод конечных элементов для инженеров» . ^[3]

История

Хотя трудно указать дату изобретения метода конечных элементов, этот метод возник из-за необходимости решения сложных задач упругости и структурного анализа в гражданской и авиационной технике . ^[4] Его развитие можно проследить до работ Александра Хренникова ^[5] и Рихарда Куранта ^[6] в начале 1940-х годов. Другим пионером был Иоаннис Аргирис . В СССР внедрение метода в практику обычно связывают с именем Леонарда Оганесяна. ^[7] Он также был независимо переоткрыт в Китае Фэн Кангом в конце 1950-х и начале 1960-х годов на основе расчетов конструкций плотин, где он был назван методом конечных разностей, основанным на вариационном принципе . Хотя подходы, используемые этими пионерами, различны, у них есть одна важная характеристика: дискретизация сетки непрерывной области на набор дискретных подобластей, обычно называемых элементами.

Работа Хренникова дискретизирует область, используя аналогию с решеткой , в то время как подход Куранта делит область на конечные треугольные подобласти для решения эллиптических уравнений в частных производных второго порядка , которые возникают из проблемы кручения цилиндра . Вклад Куранта был эволюционным, опираясь на большой объем более ранних результатов для PDE, разработанных лордом Рэлеем , Вальтером Ритцем и Борисом Галеркиным .

Метод конечных элементов получил настоящий импульс в 1960-х и 1970-х годах благодаря разработкам Дж. Х. Аргириса с коллегами из Штутгартского университета , Р. В. Клафа с коллегами из Калифорнийского университета в Беркли , О. К. Зенкевича с коллегами Эрнестом Хинтоном , Брюсом Айронсом. ^[8] и другие из Университета Суонси , Филипп Дж. Сиарле из Парижского университета 6 и Ричард Галлахер с коллегами из Корнеллского университета . Дальнейший импульс в эти годы дали доступные программы конечных элементов с открытым исходным кодом. НАСА спонсировало оригинальную версию NASTRAN . Калифорнийский университет в Беркли сделал широко доступной программу конечных элементов SAP IV ^[9] . В Норвегии общество по классификации судов Det Norske Veritas (ныне DNV GL ) разработало Sesam в 1969 году для использования при анализе судов. ^[10] Строгая математическая основа метода конечных элементов была предоставлена ​​в 1973 году публикацией Гилберта Стрэнга и Джорджа Фикса . ^[11] С тех пор метод был обобщен для численного моделирования физических систем в самых разных инженерных дисциплинах, например, электромагнетизме , теплопередаче и гидродинамике . ^{[12] [13]}

Техническое обсуждение

Структура методов конечных элементов

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

Примерами вариационной постановки являются метод Галёркина , разрывный метод Галёркина, смешанные методы и т. д.

Под стратегией дискретизации понимается четко определенный набор процедур, которые охватывают (а) создание сеток конечных элементов, (б) определение базисной функции на эталонных элементах (также называемых функциями формы) и (в) отображение опорные элементы на элементы сетки. Примерами стратегий дискретизации являются h-версия, p-версия , hp-версия , x-FEM , изогеометрический анализ и т. д. Каждая стратегия дискретизации имеет определенные преимущества и недостатки. Разумным критерием выбора стратегии дискретизации является достижение почти оптимальной производительности для самого широкого набора математических моделей в конкретном классе моделей.

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

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

Иллюстративные задачи P1 и P2

Следующие две задачи демонстрируют метод конечных элементов.

P1 — одномерная задача

$${\displaystyle {\text{ P1 }}:{\begin{cases}u''(x)=f(x){\text{ in }}(0,1),\\u(0)=u(1)=0,\end{cases}}}$$

${\displaystyle f}$ ${\displaystyle u}$ ${\displaystyle x}$ ${\displaystyle u''}$ ${\displaystyle u}$ ${\displaystyle x}$

P2 — двумерная задача ( задача Дирихле )

$${\displaystyle {\text{P2 }}:{\begin{cases}u_{xx}(x,y)+u_{yy}(x,y)=f(x,y)&{\text{ in }}\Omega ,\\u=0&{\text{ on }}\partial \Omega ,\end{cases}}}$$

где – связная открытая область на плоскости с хорошей границей (например, гладкое многообразие или многоугольник ), и обозначают вторые производные по и соответственно. ${\displaystyle \Omega }$ ${\displaystyle (x,y)}$ ${\displaystyle \partial \Omega }$ ${\displaystyle u_{xx}}$ ${\displaystyle u_{yy}}$ ${\displaystyle x}$ ${\displaystyle y}$

Задачу P1 можно решить непосредственно путем вычисления первообразных . Однако этот метод решения краевой задачи (БВП) работает только при наличии одного пространственного измерения. Это не распространяется на проблемы более высокой размерности или такие проблемы, как . По этой причине мы разработаем метод конечных элементов для P1 и наметим его обобщение на P2. ${\displaystyle u+u''=f}$

Наше объяснение будет состоять из двух этапов, которые отражают два основных шага, которые необходимо предпринять для решения краевой задачи (BVP) с использованием FEM.

• На первом этапе необходимо перефразировать исходный BVP в его слабой форме. Для этого шага обычно не требуется никаких вычислений. Трансформация выполняется вручную на бумаге.
• Второй шаг — дискретизация, при которой слабая форма дискретизируется в конечномерном пространстве.

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

Слабая формулировка

Первым шагом является преобразование P1 и P2 в их эквивалентные слабые формулировки .

Слабая форма P1

Если решает P1, то для любой гладкой функции , которая удовлетворяет граничным условиям смещения, т. е. при и , мы имеем ${\displaystyle u}$ ${\displaystyle v}$ ${\displaystyle v=0}$ ${\displaystyle x=0}$ ${\displaystyle x=1}$

И наоборот, если для каждой гладкой функции выполняется условие (1) , то можно показать, что это позволит решить проблему P1. Доказательство проще для дважды непрерывно дифференцируемого ( теорема о среднем значении ), но его также можно доказать в смысле распределения . ${\displaystyle u}$ ${\displaystyle u(0)=u(1)=0}$ ${\displaystyle v(x)}$ ${\displaystyle u}$ ${\displaystyle u}$

Мы определяем новый оператор или отображение , используя интегрирование по частям в правой части (1): ${\displaystyle \phi (u,v)}$

где мы использовали предположение, что . ${\displaystyle v(0)=v(1)=0}$

Слабая форма P2

Если мы интегрируем по частям, используя форму тождеств Грина , мы увидим, что если решает P2, то мы можем определить для любого по ${\displaystyle u}$ ${\displaystyle \phi (u,v)}$ ${\displaystyle v}$

$${\displaystyle \int _{\Omega }fv\,ds=-\int _{\Omega }\nabla u\cdot \nabla v\,ds\equiv -\phi (u,v),}$$

где обозначает градиент и обозначает скалярное произведение в двумерной плоскости. Еще раз можно превратить в скалярный продукт на подходящем пространстве однажды дифференцируемых функций, которые равны нулю на . Мы также предполагали, что (см. пространства Соболева ). Также можно показать существование и единственность решения. ${\displaystyle \nabla }$ ${\displaystyle \cdot }$ ${\displaystyle \,\!\phi }$ ${\displaystyle H_{0}^{1}(\Omega )}$ ${\displaystyle \Omega }$ ${\displaystyle \partial \Omega }$ ${\displaystyle v\in H_{0}^{1}(\Omega )}$

Схема доказательства существования и единственности решения.

Мы можем условно считать их абсолютно непрерывными функциями, которые находятся в точках и (см. Пространства Соболева ). Такие функции (слабо) однажды дифференцируемы, и оказывается, что симметричное билинейное отображение тогда определяет скалярное произведение , которое превращается в гильбертово пространство (детальное доказательство нетривиально). С другой стороны, левая часть также является скалярным произведением, на этот раз в пространстве Lp . Применение теоремы о представлении Рисса для гильбертовых пространств показывает, что существует единственное решение (2) и, следовательно, P1. Это решение априори является только членом , но с использованием эллиптической регулярности будет гладким, если оно есть. ${\displaystyle H_{0}^{1}(0,1)}$ ${\displaystyle (0,1)}$ ${\displaystyle 0}$ ${\displaystyle x=0}$ ${\displaystyle x=1}$ ${\displaystyle \!\,\phi }$ ${\displaystyle H_{0}^{1}(0,1)}$ ${\displaystyle \int _{0}^{1}f(x)v(x)dx}$ ${\displaystyle L^{2}(0,1)}$ ${\displaystyle u}$ ${\displaystyle H_{0}^{1}(0,1)}$ ${\displaystyle f}$

Дискретизация

Функция с нулевыми значениями в конечных точках (синий) и кусочно-линейной аппроксимацией (красный) ${\displaystyle H_{0}^{1},}$

P1 и P2 готовы к дискретизации, что приводит к общей подзадаче (3). Основная идея состоит в том, чтобы заменить бесконечномерную линейную задачу:

Найдите такое, что ${\displaystyle u\in H_{0}^{1}}$

${\displaystyle \forall v\in H_{0}^{1},\;-\phi (u,v)=\int fv}$

с конечномерной версией:

где – конечномерное подпространство . Существует множество возможных вариантов (один из вариантов ведет к спектральному методу ). Однако в качестве пространства кусочно-полиномиальных функций мы возьмем метод конечных элементов. ${\displaystyle V}$ ${\displaystyle H_{0}^{1}}$ ${\displaystyle V}$ ${\displaystyle V}$

Для проблемы P1

Берем интервал , выбираем значения с и определяем по: ${\displaystyle (0,1)}$ ${\displaystyle n}$ ${\displaystyle x}$ ${\displaystyle 0=x_{0}<x_{1}<\cdots <x_{n}<x_{n+1}=1}$ ${\displaystyle V}$

$${\displaystyle V=\{v:[0,1]\to \mathbb {R} \;:v{\text{ is continuous, }}v|_{[x_{k},x_{k+1}]}{\text{ is linear for }}k=0,\dots ,n{\text{, and }}v(0)=v(1)=0\}}$$

где мы определяем и . Заметим, что функции в не дифференцируемы согласно элементарному определению исчисления. Действительно, если тогда производная обычно не определена ни при каком , . Однако производная существует при любом другом значении , и эту производную можно использовать для интегрирования по частям . ${\displaystyle x_{0}=0}$ ${\displaystyle x_{n+1}=1}$ ${\displaystyle V}$ ${\displaystyle v\in V}$ ${\displaystyle x=x_{k}}$ ${\displaystyle k=1,\ldots ,n}$ ${\displaystyle x}$

Кусочно-линейная функция в двух измерениях

Для проблемы P2

Нам нужен набор функций . На рисунке справа мы проиллюстрировали триангуляцию 15 -сторонней многоугольной области на плоскости (внизу) и кусочно-линейную функцию (вверху, в цвете) этого многоугольника, линейную на каждом треугольнике триангуляции; пространство будет состоять из функций, линейных в каждом треугольнике выбранной триангуляции. ${\displaystyle V}$ ${\displaystyle \Omega }$ ${\displaystyle \Omega }$ ${\displaystyle V}$

Можно надеяться, что по мере того, как лежащая в основе треугольная сетка становится все более мелкой, решение дискретной задачи (3) в некотором смысле будет сходиться к решению исходной краевой задачи P2. Чтобы измерить тонкость сетки, триангуляция индексируется параметром с действительным значением, который считается очень маленьким. Этот параметр будет связан с наибольшим или средним размером треугольника в триангуляции. По мере уточнения триангуляции пространство кусочно-линейных функций также должно меняться с увеличением . По этой причине часто читают вместо того, чтобы читать литературу. Поскольку мы не проводим такого анализа, мы не будем использовать эти обозначения. ${\displaystyle h>0}$ ${\displaystyle V}$ ${\displaystyle h}$ ${\displaystyle V_{h}}$ ${\displaystyle V}$

Выбор основы

Интерполяция функции Бесселя

Для завершения дискретизации необходимо выбрать базис . В одномерном случае для каждой контрольной точки выберем кусочно-линейную функцию, значение которой равно at и нулю в каждой , т. е. ${\displaystyle V}$ ${\displaystyle x_{k}}$ ${\displaystyle v_{k}}$ ${\displaystyle V}$ ${\displaystyle 1}$ ${\displaystyle x_{k}}$ ${\displaystyle x_{j},\;j\neq k}$

$${\displaystyle v_{k}(x)={\begin{cases}{x-x_{k-1} \over x_{k}\,-x_{k-1}}&{\text{ if }}x\in [x_{k-1},x_{k}],\\{x_{k+1}\,-x \over x_{k+1}\,-x_{k}}&{\text{ if }}x\in [x_{k},x_{k+1}],\\0&{\text{ otherwise}},\end{cases}}}$$

для ; этот базис представляет собой сдвинутую и масштабированную палаточную функцию . Для двумерного случая мы снова выбираем одну базисную функцию на каждую вершину триангуляции плоской области . Функция — это уникальная функция, значение которой равно нулю в каждом месте . ${\displaystyle k=1,\dots ,n}$ ${\displaystyle v_{k}}$ ${\displaystyle x_{k}}$ ${\displaystyle \Omega }$ ${\displaystyle v_{k}}$ ${\displaystyle V}$ ${\displaystyle 1}$ ${\displaystyle x_{k}}$ ${\displaystyle x_{j},\;j\neq k}$

В зависимости от автора, слово «элемент» в «методе конечных элементов» относится к треугольникам области, кусочно-линейной базисной функции или к тому и другому. Так, например, автор, интересующийся изогнутыми областями, может заменить треугольники изогнутыми примитивами и таким образом описать элементы как криволинейные. С другой стороны, некоторые авторы заменяют «кусочно-линейный» на «кусочно-квадратичный» или даже «кусочно-полиномиальный». Тогда автор мог бы сказать «элемент более высокого порядка» вместо «полином более высокой степени». Метод конечных элементов не ограничивается треугольниками (тетраэдрами в трехмерных симплексах или симплексах более высокого порядка в многомерных пространствах). Тем не менее, его можно определить на четырехугольных подобластях (шестигранниках, призмах или трехмерных пирамидах и т. д.). Формы более высокого порядка (криволинейные элементы) могут определяться полиномиальными и даже неполиномиальными формами (например, эллипсом или кругом).

Примерами методов, которые используют кусочно-полиномиальные базисные функции более высокой степени, являются hp-FEM и спектральный FEM .

В более продвинутых реализациях (адаптивные методы конечных элементов) используется метод оценки качества результатов (основанный на теории оценки ошибок) и изменение сетки во время решения с целью достижения приближенного решения в некоторых пределах от точного решения непрерывной задачи. . Адаптивность сетки может использовать различные методы; наиболее популярными являются:

• движущиеся узлы (r-адаптивность)
• перерабатывающие (и неочищенные) элементы (h-адаптивность)
• изменение порядка базовых функций (p-адаптивность)
• комбинации вышеперечисленного ( hp-adaptivity ).

Небольшая поддержка основания

Решение двумерной задачи в диске с центром в начале координат и радиусом 1 с нулевыми граничными условиями. а) Триангуляция. ${\displaystyle u_{xx}+u_{yy}=-4}$

(б) Разреженная матрица L дискретизированной линейной системы

(c) Вычисленное решение, ${\displaystyle u(x,y)=1-x^{2}-y^{2}}$

Основное преимущество такого выбора базиса состоит в том, что внутренние продукты

$${\displaystyle \langle v_{j},v_{k}\rangle =\int _{0}^{1}v_{j}v_{k}\,dx}$$

$${\displaystyle \phi (v_{j},v_{k})=\int _{0}^{1}v_{j}'v_{k}'\,dx}$$

матрица Граманосителем ${\displaystyle j,k}$ ${\displaystyle \langle v_{j},v_{k}\rangle }$ ${\displaystyle (j,k)}$ ${\displaystyle v_{k}}$ ${\displaystyle [x_{k-1},x_{k+1}]}$ ${\displaystyle \langle v_{j},v_{k}\rangle }$ ${\displaystyle \phi (v_{j},v_{k})}$ ${\displaystyle |j-k|>1}$

Аналогично, в плоском случае, если и не имеют общего ребра триангуляции, то интегралы ${\displaystyle x_{j}}$ ${\displaystyle x_{k}}$

$${\displaystyle \int _{\Omega }v_{j}v_{k}\,ds}$$

$${\displaystyle \int _{\Omega }\nabla v_{j}\cdot \nabla v_{k}\,ds}$$

Матричная форма задачи

Если написать и тогда задача (3), приняв за , примет вид ${\displaystyle u(x)=\sum _{k=1}^{n}u_{k}v_{k}(x)}$ ${\displaystyle f(x)=\sum _{k=1}^{n}f_{k}v_{k}(x)}$ ${\displaystyle v(x)=v_{j}(x)}$ ${\displaystyle j=1,\dots ,n}$

Если мы обозначим через и векторы-столбцы и , и если мы позволим ${\displaystyle \mathbf {u} }$ ${\displaystyle \mathbf {f} }$ ${\displaystyle (u_{1},\dots ,u_{n})^{t}}$ ${\displaystyle (f_{1},\dots ,f_{n})^{t}}$

$${\displaystyle L=(L_{ij})}$$

$${\displaystyle M=(M_{ij})}$$

$${\displaystyle L_{ij}=\phi (v_{i},v_{j})}$$

$${\displaystyle M_{ij}=\int v_{i}v_{j}dx}$$

Не обязательно предполагать . Для общей функции задача (3) с for становится фактически проще, так как матрица не используется, ${\displaystyle f(x)=\sum _{k=1}^{n}f_{k}v_{k}(x)}$ ${\displaystyle f(x)}$ ${\displaystyle v(x)=v_{j}(x)}$ ${\displaystyle j=1,\dots ,n}$ ${\displaystyle M}$

где и для чего . ${\displaystyle \mathbf {b} =(b_{1},\dots ,b_{n})^{t}}$ ${\displaystyle b_{j}=\int fv_{j}dx}$ ${\displaystyle j=1,\dots ,n}$

Как мы уже обсуждали ранее, большинство записей и равны нулю, поскольку базисные функции имеют небольшую поддержку. Итак, теперь нам нужно решить линейную систему с неизвестным , где большинство элементов матрицы , которую нам нужно инвертировать, равны нулю. ${\displaystyle L}$ ${\displaystyle M}$ ${\displaystyle v_{k}}$ ${\displaystyle \mathbf {u} }$ ${\displaystyle L}$

Такие матрицы известны как разреженные матрицы , и для таких задач существуют эффективные средства решения (гораздо более эффективные, чем фактическое обращение матрицы). Кроме того, они симметричны и положительно определены, поэтому предпочтительным является такой метод, как метод сопряженных градиентов . Для не слишком больших задач по-прежнему хорошо работают разреженные LU-разложения и разложения Холецкого . Например, оператора обратной косой черты MATLAB (который использует разреженный LU, разреженный Холецкий и другие методы факторизации) может быть достаточно для сеток с сотней тысяч вершин. ${\displaystyle L}$

Матрицу обычно называют матрицей жесткости , а матрицу называют матрицей масс . ${\displaystyle L}$ ${\displaystyle M}$

Общая форма метода конечных элементов

В целом метод конечных элементов характеризуется следующим процессом.

• Выбирается сетка для . В предыдущем случае сетка состояла из треугольников, но можно также использовать квадраты или криволинейные многоугольники. ${\displaystyle \Omega }$
• Затем выбираются базисные функции. В нашем обсуждении мы использовали кусочно-линейные базисные функции, но обычно используются кусочно-полиномиальные базисные функции.

Отдельного рассмотрения заслуживает гладкость базисных функций. Для эллиптических краевых задач второго порядка достаточно кусочно-полиномиальной базисной функции, которая является просто непрерывной (т. е. производные разрывны). Для уравнений в частных производных более высокого порядка необходимо использовать более гладкие базисные функции. Например, для задачи четвертого порядка, такой как , можно использовать кусочно-квадратичные базисные функции, которые равны . ${\displaystyle u_{xxxx}+u_{yyyy}=f}$ ${\displaystyle C^{1}}$

Еще одним соображением является отношение конечномерного пространства к его бесконечномерному аналогу в приведенных выше примерах . Метод соответствующих элементов — это метод, в котором пространство является подпространством пространства элементов непрерывной задачи. Пример выше является именно таким методом. Если это условие не выполняется, мы получаем метод несогласных элементов, примером которого является пространство кусочно-линейных функций на сетке, непрерывных в каждой средней точке ребра. Поскольку эти функции, вообще говоря, разрывны по ребрам, это конечномерное пространство не является подпространством исходного . ${\displaystyle V}$ ${\displaystyle H_{0}^{1}}$ ${\displaystyle V}$ ${\displaystyle H_{0}^{1}}$

Обычно имеется алгоритм разделения заданной сетки. Если основным методом повышения точности является разделение сетки, есть h -метод ( h обычно представляет собой диаметр наибольшего элемента в сетке). Таким образом, если кто-то показывает, что ошибка с сеткой ограничена сверху по , для некоторых и , то имеется метод порядка p . При определенных гипотезах (например, если область выпуклая) метод кусочного полинома порядка будет иметь ошибку порядка . ${\displaystyle h}$ ${\displaystyle Ch^{p}}$ ${\displaystyle C<\infty }$ ${\displaystyle p>0}$ ${\displaystyle d}$ ${\displaystyle p=d+1}$

Если вместо уменьшения h увеличить степень полиномов, используемых в базисной функции, получится p -метод. Если объединить эти два типа уточнения, получится hp -метод ( hp-FEM ). В hp-FEM степени полинома могут варьироваться от элемента к элементу. Методы высокого порядка с большим равномерным p называются спектральными методами конечных элементов ( SFEM ). Их не следует путать со спектральными методами .

Для векторных уравнений в частных производных базисные функции могут принимать значения в . ${\displaystyle \mathbb {R} ^{n}}$

Различные типы методов конечных элементов

АЕМ

Метод прикладных элементов или AEM сочетает в себе функции FEM и метода дискретных элементов или (DEM).

А-ФЭМ

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

CutFEM

Подход Cut Finite Element был разработан в 2014 году. ^[14] Этот подход заключается в том, чтобы «сделать дискретизацию как можно более независимой от геометрического описания и минимизировать сложность построения сетки, сохраняя при этом точность и надежность стандартного метода конечных элементов. " ^[15]

Обобщенный метод конечных элементов

Обобщенный метод конечных элементов (GFEM) использует локальные пространства, состоящие из функций, не обязательно полиномов, которые отражают имеющуюся информацию о неизвестном решении и, таким образом, обеспечивают хорошую локальную аппроксимацию. Затем используется разбиение единицы , чтобы «связать» эти пространства вместе, чтобы сформировать аппроксимирующее подпространство. Эффективность GFEM была показана при применении к проблемам с областями, имеющими сложные границы, проблемам с микромасштабами и проблемам с пограничными слоями. ^[16]

Смешанный метод конечных элементов

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

Переменная – полиномиальная

HP -FEM адаптивно комбинирует элементы с переменным размером h и полиномиальной степенью p для достижения исключительно быстрой экспоненциальной скорости сходимости. ^[17]

hpk-FEM

hpk-FEM адаптивно комбинирует элементы с переменным размером h , полиномиальной степенью локальных аппроксимаций p и глобальной дифференцируемостью локальных аппроксимаций ( k -1) для достижения наилучших показателей сходимости.

КСЭМ

Расширенный метод конечных элементов (XFEM) — это численный метод, основанный на обобщенном методе конечных элементов (GFEM) и методе разделения единицы (PUM). Он расширяет классический метод конечных элементов, обогащая пространство решений решениями дифференциальных уравнений с разрывными функциями. Расширенные методы конечных элементов обогащают аппроксимационное пространство, позволяя естественным образом воспроизводить сложные особенности, связанные с интересующей задачей: разрыв, особенность, пограничный слой и т. д. Было показано, что для некоторых задач такое встраивание признака задачи в аппроксимационное пространство может значительно улучшить скорость сходимости и точность. Более того, обработка проблем с разрывами с помощью XFEM устраняет необходимость создания сетки и повторного построения сетки поверхностей разрывов, тем самым уменьшая вычислительные затраты и ошибки проецирования, связанные с традиционными методами конечных элементов, за счет ограничения разрывов краями сетки.

Некоторые исследовательские коды реализуют этот метод в разной степени:

1. Получить ФЕМ++
2. xfem++
3. openxfem++

XFEM также реализован в таких программах, как Altair Radios, ASTER, Morfeo и Abaqus. Он все чаще используется в других коммерческих программах конечных элементов, при этом доступно несколько плагинов и реальных основных реализаций (ANSYS, SAMCEF, OOFELIE и т. д.).

Метод конечных элементов с масштабированной границей (SBFEM)

Метод конечных элементов с масштабированной границей (SBFEM) был предложен Song and Wolf (1997). ^[18] SBFEM внес один из наиболее полезных вкладов в область численного анализа проблем механики разрушения. Это полуаналитический фундаментальный метод без решения, сочетающий в себе преимущества формулировок и процедур конечных элементов и дискретизации граничных элементов. Однако, в отличие от метода граничных элементов, никакого фундаментального дифференциального решения не требуется.

S-FEM

S-FEM, сглаженные методы конечных элементов, представляет собой особый класс алгоритмов численного моделирования физических явлений. Он был разработан путем объединения бессеточных методов с методом конечных элементов.

Метод спектральных элементов

Методы спектральных элементов сочетают в себе геометрическую гибкость конечных элементов и высокую точность спектральных методов. Спектральные методы представляют собой приближенное решение уравнений в частных производных слабой формы, основанное на лагранжевых интерполянтах высокого порядка и используемое только с определенными правилами квадратур. ^[19]

Бессеточные методы

Разрывные методы Галёркина

Предельный анализ методом конечных элементов

Метод растянутой сетки

Итерация Лубиньяка

Итерация Лубиньяка — это итерационный метод в методах конечных элементов.

Метод конечных элементов кристаллопластичности (CPFEM)

Метод конечных элементов кристаллопластичности (CPFEM) — это усовершенствованный численный инструмент, разработанный Францем Ротерсом. Металлы можно рассматривать как кристаллические агрегаты, которые ведут себя анизотропно при деформации, например, при аномальном напряжении и локализации деформации. CPFEM, основанный на скольжении (скорости деформации сдвига), может рассчитывать дислокацию, ориентацию кристаллов и другую информацию о текстуре, чтобы учитывать анизотропию кристалла во время процедуры. Он применялся при численном исследовании деформации материала, шероховатости поверхности, трещин и т. Д.

Метод виртуальных элементов (VEM)

Метод виртуальных элементов (VEM), предложенный Бейраном да Вейгой и др. (2013) ^[20] как расширение миметических методов конечных разностей (MFD) является обобщением стандартного метода конечных элементов для произвольной геометрии элементов. Это позволяет допускать обычные многоугольники (или многогранники в 3D), которые имеют очень неправильную и невыпуклую форму. Название «виртуальный» происходит от того факта, что знание локального базиса функции формы не требуется и фактически никогда не рассчитывается явно.

Связь с методом градиентной дискретизации

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

Сравнение с методом конечных разностей

Метод конечных разностей (FDM) является альтернативным способом аппроксимации решений УЧП. Различия между FEM и FDM:

• Наиболее привлекательной особенностью FEM является его способность относительно легко обрабатывать сложную геометрию (и границы). Хотя FDM в своей базовой форме ограничивается обработкой прямоугольных форм и простыми их изменениями, обработка геометрии в FEM теоретически проста. ^{[2] [21]}
• FDM обычно не используется для нестандартной геометрии САПР, но чаще используется для прямоугольных или блочных моделей. ^[22]
• FEM обычно обеспечивает более гибкую адаптивность сетки, чем FDM. ^[21]
• Наиболее привлекательной особенностью конечных разностей является простота реализации. ^[21]
• Можно рассматривать FDM как частный случай подхода FEM по нескольким причинам. Например, FEM первого порядка идентичен FDM для уравнения Пуассона , если задача дискретизируется с помощью регулярной прямоугольной сетки, где каждый прямоугольник разделен на два треугольника.
• Есть причины считать математическую основу аппроксимации методом конечных элементов более надежной, например, потому, что качество аппроксимации между точками сетки в FDM низкое.
• Качество аппроксимации FEM часто выше, чем в соответствующем подходе FDM, но это сильно зависит от проблемы, и можно привести несколько примеров обратного.

Как правило, FEM является методом выбора во всех типах анализа в строительной механике (т. е. при расчете деформаций и напряжений в твердых телах или динамике конструкций). Напротив, вычислительная гидродинамика (CFD), как правило, использует FDM или другие методы, такие как метод конечного объема (FVM). Проблемы CFD обычно требуют дискретизации проблемы на большое количество ячеек/точек сетки (миллионы и более). Поэтому стоимость решения благоприятствует более простому приближению более низкого порядка внутри каждой ячейки. Это особенно актуально для задач «внешнего потока», таких как поток воздуха вокруг автомобиля, самолета или симуляция погоды.

Приложение

3D-модель переноса загрязнений – поле концентрации на уровне земли

3D-модель переноса загрязнений - поле концентрации на перпендикулярной поверхности

Различные специализации в области машиностроения (например, авиационная, биомеханическая и автомобильная промышленность) обычно используют интегрированный FEM при проектировании и разработке своей продукции. Некоторые современные пакеты FEM включают в себя конкретные компоненты, такие как тепловые, электромагнитные, жидкостные и структурные рабочие среды. При структурном моделировании FEM чрезвычайно помогает в визуализации жесткости и прочности, а также в минимизации веса, материалов и затрат. ^[23]

FEM позволяет детально визуализировать места изгиба или скручивания конструкций, указывая распределение напряжений и смещений. Программное обеспечение FEM предоставляет широкий спектр возможностей моделирования для управления сложностью моделирования и системного анализа. Аналогичным образом, требуемый уровень точности и связанные с ним требования ко времени вычислений могут управляться одновременно для решения большинства инженерных приложений. FEM позволяет создавать, уточнять и оптимизировать целые конструкции еще до их изготовления. Сетка является неотъемлемой частью модели, и для получения наилучших результатов ее необходимо тщательно контролировать. Как правило, чем больше элементов в сетке, тем точнее решение дискретизированной задачи. Однако существует значение, при котором результаты сходятся, и дальнейшее измельчение сетки не увеличивает точность. ^[24]

Конечно-элементная модель коленного сустава человека ^[25]

Этот мощный инструмент проектирования значительно улучшил как стандарты инженерного проектирования, так и методологию процесса проектирования во многих промышленных приложениях. ^[26] Внедрение FEM существенно сократило время доставки продукции от концепции до производственной линии. ^[26] Тестирование и разработка были ускорены в первую очередь за счет улучшения конструкции первоначальных прототипов с использованием FEM. ^[27] Таким образом, преимущества FEM включают повышенную точность, улучшенное проектирование и лучшее понимание критических параметров проектирования, виртуальное прототипирование, меньшее количество прототипов оборудования, более быстрый и менее дорогой цикл проектирования, повышение производительности и увеличение доходов. ^[26]

В 1990-х годах метод FEM был предложен для использования в стохастическом моделировании для численного решения вероятностных моделей ^[28], а затем и для оценки надежности. ^[29]

Смотрите также

• Метод прикладного элемента
• Метод граничных элементов
• Лемма Сеа
• Компьютерный эксперимент
• Метод прямой жесткости
• Оптимизация компоновки разрывов
• Метод дискретных элементов
• Метод конечных разностей
• Машина конечных элементов
• Метод конечных элементов в строительной механике
• Метод конечного объема
• Метод конечных объемов для нестационарного течения
• Метод бесконечных элементов
• Интервальный конечный элемент
• Изогеометрический анализ
• Решеточные методы Больцмана
• Список пакетов программного обеспечения для конечных элементов
• Бессеточные методы
• Передвижной клеточный автомат
• Многопрофильная оптимизация проектирования
• Мультифизика
• Патч-тест
• Метод Рэлея-Ритца
• Картографирование пространства
• СТРЭНД7
• Тесселяция (компьютерная графика)
• Ослабленная слабая форма

Рекомендации

1. Дэрил Л. Логан (2011). Первый курс метода конечных элементов . Cengage Обучение. ISBN 9780495668275.
2. Редди, JN (2006). Введение в метод конечных элементов (Третье изд.). МакГроу-Хилл. ISBN 9780071267618.
3. Хюбнер, Кеннет Х. (2001). Метод конечных элементов для инженеров . Уайли. ISBN 978-0-471-37078-9.
4. Лю, Винг Кам; Ли, Шаофань; Парк, Гарольд С. (2022). «Восемьдесят лет метода конечных элементов: рождение, эволюция и будущее». Архив вычислительных методов в технике . 29 (6): 4431–4453. arXiv : 2107.04960 . дои : 10.1007/s11831-022-09740-9 . ISSN  1134-3060. S2CID  235794921.
5. Хренников, Александр (1941). «Решение задач упругости каркасным методом». Журнал прикладной механики . 8 (4): 169–175. Бибкод : 1941JAM.....8A.169H. дои : 10.1115/1.4009129.
6. Курант, Р. (1943). «Вариационные методы решения задач равновесия и колебаний». Бюллетень Американского математического общества . 49 : 1–23. дои : 10.1090/s0002-9904-1943-07818-4 .
7. "СПб ЭМИ РАН". emi.nw.ru. _ Архивировано из оригинала 30 сентября 2015 года . Проверено 17 марта 2018 г.
8. Хинтон, Эрнест; Айронс, Брюс (июль 1968 г.). «Сглаживание экспериментальных данных методом наименьших квадратов с использованием конечных элементов». Напряжение . 4 (3): 24–27. doi :10.1111/j.1475-1305.1968.tb01368.x.
9. «Программное обеспечение и руководства SAP-IV». Электронная библиотека NISEE, Интернет-архив сейсмической инженерии. Архивировано из оригинала 9 марта 2013 г. Проверено 24 января 2013 г.
10. Гард Полсен; Хокон с Андерсеном; Джон Петтер Коллетт; Ивер Танген Стенсруд (2014). Building Trust, История DNV 1864-2014 . Лисакер, Норвегия: Динамо Форлаг А/С. стр. 121, 436. ISBN. 978-82-8071-256-1.
11. Стрэнг, Гилберт ; Фикс, Джордж (1973). Анализ метода конечных элементов . Прентис Холл. ISBN 978-0-13-032946-2.
12. Олек С. Зенкевич; Роберт Л. Тейлор; Дж. З. Чжу (31 августа 2013 г.). Метод конечных элементов: его основы и основы. Баттерворт-Хайнеманн. ISBN 978-0-08-095135-5.
13. Бате, KJ (2006). Процедуры конечных элементов . Кембридж, Массачусетс: Клаус-Юрген Бат. ISBN 978-0979004902.
14. Челедони (27 февраля 2023 г.). «CutFEM: дискретизация уравнений в частных производных и геометрии». ЭКМИ . Проверено 13 октября 2023 г.
15. Бурман, Эрик; Клаус, Сюзанна; Хансбо, Питер; Ларсон, Матс Г.; Массинг, Андре (16 ноября 2015 г.). «CutFEM: Дискретизация геометрии и уравнений в частных производных». Международный журнал численных методов в технике . 104 (7): 472–501. дои : 10.1002/nme.4823 . ISSN  0029-5981.
16. Бабушка, Иво ; Банерджи, Удай; Осборн, Джон Э. (июнь 2004 г.). «Обобщенные методы конечных элементов: основные идеи, результаты и перспективы». Международный журнал вычислительных методов . 1 (1): 67–103. дои : 10.1142/S0219876204000083.
17. П. Солин, К. Сегет, И. Долезель: Методы конечных элементов высшего порядка, Chapman & Hall/CRC Press, 2003
18. Сон, Чхонмин; Вольф, Джон П. (5 августа 1997 г.). «Масштабированный граничный метод конечных элементов - псевдоним последовательного метода бесконечно малых ячеек конечных элементов - для эластодинамики». Компьютерные методы в прикладной механике и технике . 147 (3–4): 329–355. Бибкод : 1997CMAME.147..329S. дои : 10.1016/S0045-7825(97)00021-2.
19. «Методы спектральных элементов». Государственная ключевая лаборатория научных и инженерных вычислений . Архивировано из оригинала 10 августа 2017 г. Проверено 28 июля 2017 г.
20. Бейран да Вейга, Л.; Бреззи, Ф.; Канджиани, А.; Манзини, Дж.; Марини, LD; Руссо, А. (2013). «Основные принципы методов виртуальных элементов». Математические модели и методы в прикладных науках . 23 (1): 199–214. дои : 10.1142/S0218202512500492.
21. Топпер, Юрген (январь 2005 г.). «Ценообразование опционов с конечными элементами». Уилмотт . 2005 (1): 84–90. дои : 10.1002/wilm.42820050119. ISSN  1540-6962.
22. «В чем разница между FEM, FDM и FVM?». Дизайн машины . 18 апреля 2016 г. Архивировано из оригинала 28 июля 2017 г. Проверено 28 июля 2017 г.
23. Кирицис, Д.; Эммануилидис, Ч.; Корониос, А.; Мэтью, Дж. (2009). «Управление инженерными активами». Материалы 4-го Всемирного конгресса по управлению инженерными активами (WCEAM) : 591–592.
24. «Анализ методом конечных элементов: как создать отличную модель» . Ковентивные композиты . 18 марта 2019 г. Проверено 5 апреля 2019 г.^{[ постоянная мертвая ссылка ]}
25. Нагиби Бейдохти, Хамид; Янссен, Деннис; Хошгофтар, Мехди; Спренгерс, Андре; Пердахиоглу, Эмин Семих; Бугаард, Тон Ван ден; Вердоншот, Нико (2016). «Сравнение динамического неявного и явного моделирования конечного элемента собственного коленного сустава» (PDF) . Медицинская инженерия и физика . 38 (10): 1123–1130. doi :10.1016/j.medengphy.2016.06.001. PMID  27349493. Архивировано (PDF) из оригинала 19 июля 2018 г. Проверено 19 сентября 2019 г.
26. Гастингс, Дж. К., Джадс, М. А., Брауэр, Дж. Р., Точность и экономичность магнитного анализа методом конечных элементов , 33-я ежегодная национальная релейная конференция, апрель 1985 г.
27. Макларен-Мерседес (2006). «McLaren Mercedes: особенность — стресс, чтобы произвести впечатление». Архивировано из оригинала 30 октября 2006 г. Проверено 3 октября 2006 г.
28. Пэн Лонг; Ван Цзиньлян; Чжу Цидин (19 мая 1995 г.). «Методы высокой точности расчета вероятностей конечных элементов». Журнал вычислительной и прикладной математики . 59 (2): 181–189. дои : 10.1016/0377-0427(94)00027-X.
29. Халдар, Ачинтья; Махадеван, Шанкаран (2000). Оценка надежности с использованием стохастического анализа конечных элементов . Джон Уайли и сыновья. ISBN 978-0471369615.

дальнейшее чтение

• Г. Аллер и А. Крейг: Численный анализ и оптимизация: введение в математическое моделирование и численное моделирование .
• К.Дж. Бат: Численные методы анализа методом конечных элементов , Прентис-Холл (1976).
• Томас Дж. Р. Хьюз: Метод конечных элементов: линейный статический и динамический анализ методом конечных элементов, Прентис-Холл (1987).
• Дж. Часкалович: Методы конечных элементов для инженерных наук , Springer Verlag, (2008).
• Эндре Сюли : Методы конечных элементов для уравнений в частных производных.
• О. К. Зенкевич, Р. Л. Тейлор, Дж. З. Чжу: Метод конечных элементов: его основа и основы , Баттерворт-Хейнеманн (2005).
• Н. Оттосен, Х. Петерссон: Введение в метод конечных элементов, Прентис-Холл (1992).
• Сюзанна К. Бреннер, Л. Риджуэй Скотт: Математическая теория методов конечных элементов , Springer-Verlag, Нью-Йорк, ISBN 978-0-387-75933-3 (2008).
• Зохди, Т.И. (2018) Учебник по методу конечных элементов для начинающих, расширенная версия, включая примеры тестов и проектов. Второе издание https://link.springer.com/book/10.1007/978-3-319-70428-9.
• Лешек Ф. Демкович: Математическая теория конечных элементов , SIAM, ISBN 978-1-61197-772-1 (2024).