Метод граничного элемента

Материал из Википедии — свободной энциклопедии
Перейти к: навигация, поиск

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

Применялся изначально при решении задач Дирихле, Неймана — уравнение Лапласа[1].

Потом получил обобщение для уравнений теории упругости. Одним аналогом формул Грина в теории упругости являются формулы Бетти (упругие потенциалы на основе решения Кельвина-Сомилианы)[2] . Другой использовал Вейль (антенный потенциал)[3].

В. Д. Купрадзе обобщил постановку для граничных задач теории колебаний и других.[4][5][6]

Достоинства[править | править вики-текст]

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

Недостатки[править | править вики-текст]

Недостатками традиционной постановки метода являются:

  • Рассматриваются граничные условия одного типа либо Дирихле, либо Неймана, смешанная задача не рассмотрена. (Не составляет труда записать уравнения смешанной задачи, но они не имеют теории решения.)
  • Граница должна быть гладкой. (Полученные при решении задачи Неймана сингулярные интегралы не существуют в угловых точках кусочно гладкой границы.)
  • Матрица результирующей системы линейных алгебраических уравнений (СЛАУ), заменяющей интегральное, — полностью заполненная, в отличие от МКЭ, в котором она содержит большое количество нулей (хотя в МКЭ матрица на единицу размерности больше, так как сетка элементов наносится на всю область, а не только на границу).

Сложности[править | править вики-текст]

Так же к недостаткам можно отнести техническую сложность МГЭ:

  • Вычисление сингулярных интегралов представляет трудность. Они могут быть вычислены, например, с помощью формулы Стокса, после замены границы набором плоских элементов. Либо с помощью их регулярного представления (Перлин П. И.).
  • Разрешающие уравнения (обобщенные) Фредгольма второго рода, находятся на границе круга сходимости. То есть либо само уравнение, либо союзное к нему имеют собственные решения (отличные от нуля решения при нулевой правой части) . Что в частности, не позволяет искать решение внешней задачи Дирихле на основе потенциала двойного слоя, так как условие разрешимости нельзя сформулировать, — в общем случае неизвестна собственная функция союзного уравнения. (Хотя исходная задача имеет единственное решение, — потенциал двойного слоя не удовлетворяет «условию излучения»[1].) Способ перехода к модифицированным уравнениям известен. (Если к ним не переходить, то например, при решении внутренней задачи Неймана определитель матрицы СЛАУ при уменьшении характерного размера сети граничных элементов стремится к нулю.)

Трудности метода можно оценить прочитав предисловие Шермана Д. И. к[2].

В целом[править | править вики-текст]

  • Можно сказать, что в рамках традиционной постановки задачи Дирихле, Неймана (и соответствующие им теории упругости) для гладкой границы успешно решаются. Можно использовать аналитическое интегрирование (не всегда это рационально с точки зрения расхода машинных ресурсов) и метод последовательных приближений решения СЛАУ (на основе модифицированных уравнений), для доказательства сходимости которых используется теория интегральных уравнений Фредгольма второго рода.
  • Из-за сложности реализации и ограниченной сферы применения интерес к методу уменьшился. По крайней мере, заменой МКЭ, как ожидалось, он не стал.
  • Существует большое число постановок, отличных от традиционной. В том числе в тех случаях, когда математическая теория отсутствует, а уравнения записать можно. Например, решение на основе уравнения Фредгольма первого рода, для чего необходимо производить регуляризацию, иначе задача является некорректно поставленной (при небольшом изменении правой части, решение изменяется значительно). Смешанной задачи, где необходимо учесть возможное появление неограниченной производной искомой функции вблизи точки смены граничных условий даже при гладкой границе. Обобщение для кусочно гладкой границы (в плоском случае) можно производить с помощью уравнений для гладкой границы, путем введения весовых функций, полученных на основе исследования асимптотик решений для клина.
  • За рубежом есть сообщество исследователей МГЭ, — см.: «boundary element method»; переводную книгу:[7] Издаются тематические журналы.

Примечания[править | править вики-текст]

  1. 1 2 Сретенский Л. Н. Теория Ньютоновского потенциала.- М.: Гос. изд-во технико-теоретической литературы, 1946, 318 с.
  2. 1 2 Партон В. З., Перлин П. И. Интегральные уравнения теории упругости. — М.: Наука, 1977, 312 c.
  3. Вейль Г. Математика. Теоретическая физика. М.: Наука, 1984. −510 с.
  4. Купрадзе В. Д. Граничные задачи теории колебаний и интегральные уравнения. — М.: Гос. изд-во технико-теоретической литературы, 1950, 280 с.
  5. Купрадзе В. Д. Методы потенциала в теории упругости, М.: Гос. изд-во технико-теоретической литературы, 1963, 472 с.
  6. Купрадзе В. Д. Трехмерные задачи математической теории упругости и термоупругости, М.: Наука, 1976, 664 с.
  7. Кацикаделис Джон Т. Граничные элементы: теория и приложения. — М: Издательство АСВ, 2007 (Перевод книги: John T. Katsikadelis Boundary elements: Theory and applications, Oxford: Elsever, 2002, 336 c.)