Текущая версия страницы пока не проверялась опытными участниками и может значительно отличаться от версии, проверенной 5 марта 2021 года; проверки требуют 13 правок.
Текущая версия страницы пока не проверялась опытными участниками и может значительно отличаться от версии, проверенной 5 марта 2021 года; проверки требуют 13 правок.
Метод конечных элементов
Решение методом конечных элементов двухмерной магнитостатической задачи (линии и цвет означают направление и величину магнитной индукции)
Разбиение на конечные элементы. Размер элементов можно менять, уменьшая его вблизи интересующей области, и увеличивая — для снижения затрат процессорного времени
Область, в которой ищется решение дифференциальных уравнений, разбивается на конечное количество элементов. В каждом из элементов произвольно выбирается вид аппроксимирующей функции. В простейшем случае это полином первой степени. Вне своего элемента аппроксимирующая функция равна нулю. Значения функций на границах элементов (в узлах) являются решением задачи и заранее неизвестны. Коэффициенты аппроксимирующих функций обычно ищутся из условия равенства значения соседних функций на границах между элементами (в узлах). Затем эти коэффициенты выражаются через значения функций в узлах элементов. Составляется система линейных алгебраических уравнений. Количество уравнений равно количеству неизвестных значений в узлах, на которых ищется решение исходной системы, прямо пропорционально количеству элементов и ограничивается только возможностями ЭВМ. Так как каждый из элементов связан с ограниченным количеством соседних, система линейных алгебраических уравнений имеет разрежённый вид, что существенно упрощает её решение.
Если говорить в матричных терминах, то собираются так называемые матрицы жёсткости (или матрица Дирихле) и масс. Далее на эти матрицы накладываются граничные условия (например, при условиях Неймана в матрицах не меняется ничего, а при условиях Дирихле из матриц вычёркиваются строки и столбцы, соответствующие граничным узлам, так как в силу краевых условий значение соответствующих компонент решения известно). Затем собирается система линейных уравнений и решается одним из известных методов.
С точки зрения вычислительной математики, идея метода конечных элементов заключается в том, что минимизация функционала вариационной задачи осуществляется на совокупности функций, каждая из которых определена на своей подобласти.
Метод получил широкое применение при проектировании сооружений,[1] а также при моделировании моделей движения, к примеру, грунта.
Из недостатков метода стоит отметить влияние размера сетки на конечные результаты.
Функция на с нулевыми значениями на концах (голубая), и аппроксимация этой функции отрезками (красная).
Базисные функции vk (голубые) и линейная комбинация из них, которая аппроксимирует искомую функцию (красная).
Визуализация деформации машины при асимметричном ударе, выполненная по результатам расчета методом конечных элементов.
Пусть в одномерном пространстве Р1 необходимо решить следующее одномерное дифференциальное уравнение для нахождения функции на промежутке от 0 до 1. На границах области значение функции равно 0:
где известная функция, неизвестная функция от . вторая производная от по .
Решение поставленной задачи методом конечных элементов разобьём на 2 этапа:
Переформулируем граничную задачу в так называемую слабую (вариационную) форму. На этом этапе вычислений почти не требуется.
На втором этапе разобьём слабую форму на конечные отрезки-элементы.
После этого возникает проблема нахождения системы линейных алгебраических уравнений, решение которой аппроксимирует искомую функцию.
Если есть решение, то для любой гладкой функции , которая удовлетворяет граничным условиям в точках и , можно записать следующее выражение:
на конечные промежутки, и получим новое пространство :
(3) такое, что
где кусочная область пространства . Есть много способов для выбора базиса . Выберем в качестве базисных функций такие , чтобы они представлялись прямыми линиями (полиномами первой степени):
для (в данном примере )
Если теперь искомое приближённое решение представить виде , а функцию аппроксимировать как , то с помощью (3) можно получить следующую систему уравнений относительно искомых :
Метод конечных элементов сложнее метода конечных разностей в реализации. У МКЭ, однако, есть ряд преимуществ, проявляющихся на реальных задачах: произвольная форма обрабатываемой области; сетку можно сделать более редкой в тех местах, где особая точность не нужна.
Долгое время широкому распространению МКЭ мешало отсутствие алгоритмов автоматического разбиения области на «почти равносторонние» треугольники (погрешность, в зависимости от вариации метода, обратно пропорциональна синусу или самого острого, или самого тупого угла в разбиении). Впрочем, эту задачу удалось успешно решить (алгоритмы основаны на триангуляции Делоне), что дало возможность создавать полностью автоматические конечноэлементные САПР.
В Китае в 1950-х годах Фэн Кан предложил численный метод решения дифференциальных уравнений в частных производных для расчета конструкций плотин. Этот метод был назван методом конечных разностей на основе вариационного принципа, что может рассматриваться как еще один независимый способ реализации метода конечных элементов.
Хотя перечисленные подходы различаются между собой в деталях, они имеют одну общую черту: дискретизация непрерывной области сеткой в набор дискретных поддоменов, обычно называемых элементами.
Дальнейшее развитие метода конечных элементов связано также с решением задач космических исследований в 1950-х годах.
В СССР распространение и практическая реализация МКЭ в 1960-х годах связана с именем Леонарда Оганесяна.
Существенный толчок в своём развитии МКЭ получил в 1963 году после того, как было доказано, что его можно рассматривать как один из вариантов распространённого в строительной механике метода Рэлея — Ритца, который путём минимизации потенциальной энергии сводит задачу к системе линейных уравнений равновесия. После того, как была установлена связь МКЭ с процедурой минимизации, он стал применяться к задачам, описываемым уравнениями Лапласа или Пуассона. Область применения МКЭ значительно расширилась, когда было установлено (в 1968 году), что уравнения, определяющие элементы в задачах, могут быть легко получены с помощью вариантов метода взвешенных невязок, таких как метод Галёркина или метод наименьших квадратов. Это сыграло важную роль в теоретическом обосновании МКЭ, так как позволило применять его при решении многих типов дифференциальных уравнений. Таким образом, метод конечных элементов превратился в общий метод численного решения дифференциальных уравнений или систем дифференциальных уравнений.
С развитием вычислительных средств возможности метода постоянно расширяются, также расширяется и класс решаемых задач. В настоящее время предложено большое количество реализаций метода конечных элементов при моделировании процессов диффузии[2], теплопроводности[3], гидродинамики[4], механики[5], электродинамики[6] и др.
Галлагер Р. Метод конечных элементов. Основы: Пер. с англ. — М.: Мир, 1984
Деклу Ж. Метод конечных элементов: Пер. с франц. — М.: Мир, 1976
Зенкевич О. Метод конечных элементов в технике — М.: Мир, 1975.
Зенкевич О., Морган К. Конечные элементы и аппроксимация: Пер. с англ. — М.: Мир, 1986
Сегерлинд Л. Применение метода конечных элементов — М.: Мир, 1979. — 392 С.
P. Austrell, O. Dahlblom, J. Lindemann, A. Olsson, K.-G. Olsson, K. Persson, H. Pettersson, M. Ristinmaa, G. Sandberg, and P.-A. Wernberg. CALFEM: A Finite Element Toolbox : Version 3.4. Structural Mechanics, LTH, 2004. isbn: 91-8855823-1