Из Википедии, свободной энциклопедии
Перейти к навигации Перейти к поиску

hp-FEM - это общая версия метода конечных элементов (FEM), численного метода решения уравнений в частных производных, основанного на кусочно-полиномиальных аппроксимациях, в котором используются элементы переменного размера (h) и полиномиальной степени (p) . Истоки HP-FEM восходят к новаторской работе Барны А. Сабо и Иво Бабушки [1] [2] [3] [4] [5] [6]], которые обнаружили, что метод конечных элементов сходится экспоненциально быстро, когда сетка уточняется с использованием подходящей комбинации h-уточнений (разделение элементов на более мелкие) иp-уточнения (увеличение их полиномиальной степени). Экспоненциальная сходимость делает этот метод очень привлекательным по сравнению с большинством других методов конечных элементов, которые сходятся только с алгебраической скоростью. Экспоненциальная сходимость hp-FEM была не только предсказана теоретически, но и наблюдалась многочисленными независимыми исследователями. [7] [8] [9]

Отличия от стандартного МКЭ [ править ]

HP-FEM отличается от стандартного FEM (низшего порядка) во многих аспектах. [10]

  • Выбор функций формы более высокого порядка [ необходим пример ] : Для начала, полиномы более высокой степени в элементах могут быть сгенерированы с использованием различных наборов функций формы. Выбор такого набора может существенно повлиять на кондиционирование матрицы жесткости и, в свою очередь, на весь процесс решения. Эта проблема была впервые задокументирована Babuska et al. [11]
  • Автоматическая адаптация к HP : в HP-FEM элемент может быть усовершенствован множеством различных способов. Один из способов - просто увеличить степень его полинома, не разделяя его в пространстве. Или элемент может быть разделен геометрически, и к подэлементам могут быть применены различные степени полинома. Число кандидатов на уточнение элементов легко достигает 100 в 2D и 1000 в 3D. Следовательно, очевидно, что одного числа, указывающего размер ошибки в элементе, недостаточно для управления автоматической HP-адаптивностью (в отличие от адаптивности в стандартном FEM). Для получения дополнительной информации о форме ошибки в каждом элементе необходимо использовать другие методы, такие как эталонные решения или соображения аналитичности . [12]
  • Соотношение процессорного времени сборки и решения : В стандартном МКЭ матрица жесткости обычно собирается быстро, но она довольно большая. Поэтому, как правило, решение дискретной задачи занимает большую часть общего вычислительного времени. Напротив, матрицы жесткости в HP-FEM обычно намного меньше, но (для того же размера матрицы) их сборка занимает больше времени, чем в стандартном FEM. В основном это связано с вычислительными затратами на числовую квадратуру, которая должна иметь более высокую точность и, следовательно, быть более высокого порядка по сравнению со стандартным МКЭ, чтобы воспользоваться преимуществами более высокой скорости сходимости.
  • Аналитические задачи : HP-FEM труднее понять с аналитической точки зрения, чем стандартный FEM. [ согласно кому? ]Это касается многих методов, таких как дискретные принципы максимума (DMP) для эллиптических задач. Эти результаты утверждают, что, обычно с некоторыми ограничивающими предположениями о сетке, кусочно-полиномиальное приближение МКЭ подчиняется аналогичным принципам максимума, что и лежащая в основе эллиптическая УЧП. Такие результаты очень важны, поскольку они гарантируют, что приближение остается физически допустимым, не оставляя возможности вычислить отрицательную плотность, отрицательную концентрацию или отрицательную абсолютную температуру. DMP довольно хорошо изучены для FEM самого низкого порядка, но совершенно неизвестны для HP-FEM в двух или более измерениях. Первые DMP в одном пространственном измерении были сформулированы недавно. [13] [14]
  • Проблемы программирования : гораздо сложнее реализовать решатель hp-FEM, чем стандартный код FEM. Множество проблем, которые необходимо преодолеть, включают (но не ограничиваются ими): квадратурные формулы высшего порядка, функции формы высшего порядка, информацию о связности и ориентации, связывающую функции формы в эталонной области с базисными функциями в физической области и т. Д. [15]

Пример: проблема Fichera [ править ]

Проблема Fichera (также называемая угловой проблемой Fichera) является стандартной задачей эталонного теста для адаптивных кодов FEM. Его можно использовать, чтобы показать резкую разницу в производительности стандартного FEM и HP-FEM. Геометрия задачи - это куб с отсутствующим углом. Точное решение имеет сингулярный градиент (аналог бесконечного напряжения) в центре. Знание точного решения позволяет точно рассчитать ошибку аппроксимации и, таким образом, сравнить различные численные методы. Для иллюстрации проблема была решена с использованием трех различных версий адаптивного FEM: с линейными элементами, квадратичными элементами и hp-FEM.

  • Проблема Фичеры: сингулярный градиент.

  • Задача Fichera: сравнение сходимости.

Графики сходимости показывают ошибку аппроксимации как функцию числа степеней свободы (DOF). Под глубиной резкости мы подразумеваем (неизвестные) параметры, которые необходимы для определения приближения. Количество степеней свободы равно размеру матрицы жесткости. Читатель может видеть на графиках, что сходимость hp-FEM происходит намного быстрее, чем сходимость обоих других методов. На самом деле разрыв в производительности настолько огромен, что линейный МКЭ может вообще не сходиться за разумное время, а квадратичному МКЭ потребуются сотни тысяч или, возможно, миллионы степеней свободы, чтобы достичь точности, которую HP-FEM достиг примерно с 17000 степеней свободы. Получение очень точных результатов с использованием относительно небольшого количества степеней свободы - основная сила HP-FEM.

Почему HP-FEM настолько эффективен? [ редактировать ]

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

  • Слева: сетка, состоящая из двух линейных элементов.
  • Справа: сетка, состоящая из одного квадратичного элемента.

Кусочно-линейная аппроксимация.Квадратичное приближение.

Хотя количество неизвестных в обоих случаях одинаково (1 степень свободы), ошибки в соответствующей норме составляют 0,68 и 0,20 соответственно. Это означает, что квадратичное приближение было примерно в 3,5 раза эффективнее кусочно-линейного. Когда мы проделаем еще один шаг и сравним (а) четыре линейных элемента с (б) одним четвертым элементом (p = 4), то обе дискретные задачи будут иметь три степени свободы, но приближение четвертой степени будет примерно в 40 раз эффективнее. Выполнив еще несколько таких шагов, читатель увидит, что разрыв в эффективности открывается очень быстро.

Напротив, маленькие элементы низкого порядка могут улавливать мелкомасштабные особенности, такие как сингулярности, намного лучше, чем большие элементы высокого порядка. HP-FEM основан на оптимальной комбинации этих двух подходов, которая приводит к экспоненциальной сходимости. Обратите внимание, что эта экспоненциальная сходимость выражается по оси ошибки по отношению к степеням свободы. В реальных приложениях мы обычно учитываем время вычислений, необходимое для достижения того же уровня точности. Для этого индикатора производительности h- и hp-уточнение может дать аналогичные результаты, например, см. Окончательный рисунок в [16] (ссылка на WebArchive [17] ). Как только становится сложнее программировать и распараллеливать hp-FEM по сравнению с h-FEM, превосходная сходимость hp-уточнения может оказаться непрактичной.

Что такое HP-адаптивность? [ редактировать ]

Некоторые сайты FEM описывают hp-адаптивность как комбинацию h-адаптивности (разделение элементов в пространстве с сохранением их полиномиальной степени фиксированной) и p-адаптивности (только увеличение их полиномиальной степени). Это не совсем так. HP-адаптивность значительно отличается как от h-, так и p-адаптивности, поскольку hp-уточнение элемента может быть выполнено множеством различных способов. Помимо p-уточнения, элемент может быть разделен в пространстве (как в h-адаптивности), но существует множество комбинаций степеней полинома на подэлементах. Это показано на рисунке справа. Например, если треугольный или четырехугольный элемент разделен на четыре подэлемента, где степени полинома могут изменяться не более чем на два, то это дает 3 ^ 4 = 81 кандидата на уточнение (без учета полиномиально анизотропных кандидатов).Аналогично, разделение шестигранника на восемь подэлементов и изменение их полиномиальных степеней не более чем на два дает 3 ^ 8 = 6,561 кандидата на уточнение. Ясно, что стандартных оценок ошибок МКЭ, обеспечивающих одно постоянное число на элемент, недостаточно для управления автоматической HP-адаптивностью.

Функции формы высшего порядка [ править ]

В стандартном МКЭ работает только с функциями формы, связанными с вершинами сетки (так называемые вершинные функции ). В отличие от этого, в hp-FEM, кроме того, рассматриваются граничные функции (связанные с ребрами элементов), функции граней (соответствующие граням элементов - только 3D) и пузырьковые функции (полиномы более высокого порядка, которые обращаются в нуль на границах элементов). На следующих изображениях показаны эти функции (ограниченные одним элементом):

  • Вершинная функция.

  • Краевая функция.

  • Функция лица.

  • Пузырьковая функция.

Примечание: все эти функции определены во всем интерьере элемента!

Коды hp-FEM с открытым исходным кодом [ править ]

  • Deal.II : deal.II - это бесплатная библиотека с открытым исходным кодом для решения уравнений в частных производных с использованием метода конечных элементов.
  • Концепции : C / C ++ hp-FEM / DGFEM / BEM библиотека для эллиптических уравнений, разработанная в SAM, ETH Zurich (Швейцария) и в группе К. Шмидта в TU Berlin (Германия).
  • 2dhp90, 3dhp90: коды Fortran для эллиптических задач и уравнений Максвелла, разработанные Л. Демковичем из ICES, UT Austin.
  • PHAML: параллельный иерархический адаптивный многоуровневый проект. Программное обеспечение с конечными элементами, разработанное в Национальном институте стандартов и технологий, США, для численного решения двумерных эллиптических уравнений в частных производных на параллельных компьютерах с распределенной памятью и многоядерных компьютерах с использованием методов адаптивного измельчения сетки и методов многосеточного решения.
  • Проект Hermes : библиотека C / C ++ / Python для быстрого прототипирования пространственно-пространственно-временных адаптивных решателей hp-FEM для большого разнообразия PDE и мультифизических систем PDE, разработанная группой hp-FEM в Университете Невады, Рино ( США), Институт термомеханики, Прага (Чешская Республика), и Университет Западной Богемии в Пльзене (Чешская Республика) - с инженерным программным обеспечением Agros2D, созданным на основе библиотеки Hermes.
  • PHG : PHG - это набор инструментов для разработки параллельных адаптивных программ конечных элементов. Подходит для h-, p- и hp-fem. PHG в настоящее время активно разрабатывается в Государственной ключевой лаборатории научных и инженерных вычислений Института вычислительной математики и научных / инженерных вычислений Китайской академии наук (LSEC, CAS, Китай). PHG работает с соответствующими тетраэдрическими сетками и использует пополам для адаптивного локального уточнения сетки и MPI для передачи сообщений. PHG имеет объектно-ориентированный дизайн, который скрывает детали распараллеливания и предоставляет общие операции с сетками и функциями конечных элементов абстрактным способом, позволяя пользователям сосредоточиться на своих численных алгоритмах.
  • MoFEM - это код анализа методом конечных элементов, предназначенный для решения мультифизических задач с произвольными уровнями приближения, различными уровнями детализации сетки и оптимизированный для высокопроизводительных вычислений. Он предназначен для управления сложностями, связанными с неоднородным порядком приближений для пространств L2, H1, H-div и H-curl.
  • Sparselizard - это мультифизическая, адаптируемая к hp, удобная для пользователя библиотека конечных элементов C ++ с открытым исходным кодом, которая в настоящее время разрабатывается в университете Тампере, Финляндия. Он сочетает в себе трехмерные тетраэдры и двумерные треугольники / четырехугольники, конформное адаптивное уточнение сетки с иерархическими функциональными пространствами H1 и H-curl произвольного порядка для общих статических и переходных HPFEM.

Коммерческое программное обеспечение hp-FEM [ править ]

  • StressCheck - это инструмент анализа методом конечных элементов с возможностями HP, ориентированными на подробный структурный анализ.

Ссылки [ править ]

  1. ^ BA Szabó , AK Mehta: p-конвергентные приближения конечных элементов в механике разрушения, Int. J. Num. Meth. Engng, Volume 12, pp. 551-560, 1978.
  2. ^ И. Бабушка , Б. А. Сабо и И. Н. Кац: p-версия метода конечных элементов, SIAM J. Numer. Anl., Volume 18, pp. 515-544, 1981.
  3. ^ И. Бабушка , Б. А. Сабо , О скоростях сходимости метода конечных элементов, Int. J. Numer. Meth. Англ., Том 18, стр. 323-341, 1982.
  4. ^ И. Бабушка : p- и hp-версии метода конечных элементов: современное состояние, конечные элементы: теория и приложения, под редакцией DL Dwoyer, MY Hussaini и RG Voigt, Нью-Йорк, Springer-Verlag, 1988 .
  5. ^ BA Szabó , I. Babuška : Анализ методом конечных элементов, John Wiley & Sons, ISBN  978-0-471-50273-9 , 1991.
  6. ^ И. Бабушка , Б.К. Го: версии h, p и hp метода конечных элементов: базовая теория и приложения, Advances in Engineering Software, Volume 15, Issue 3-4, 1992.
  7. ^ JM Melenk: HP-методы конечных элементов для сингулярных возмущений, Springer, 2002
  8. ^ К. Шваб: p- и hp- методы конечных элементов: теория и приложения в механике твердого тела и жидкости, Oxford University Press, 1998
  9. ^ П. Солин: Уравнения в частных производных и метод конечных элементов, J. Wiley & Sons, 2005
  10. ^ П. Солин, К. Сегет, И. Долезел: Методы конечных элементов высшего порядка, Chapman & Hall / CRC Press, 2003
  11. ^ И. Бабушка, М. Грибель и Дж. Питкаранта, Проблема выбора функций формы для конечного элемента p-типа, Междунар. J. Numer. Методы Engrg. (1989), стр. 1891–1908.
  12. ^ L. Demkowicz, W. Rachowicz, и Ph. Devloo: A Fully Automatic hp-Adaptivity, Journal of Scientific Computing, 17, № 1–3 (2002), 127–155
  13. ^ П. Солин, Т. Вейходский: слабый дискретный принцип максимума для hp-FEM, J. Comput. Прил. Математика. 209 (2007) 54–65
  14. ^ Т. Вейходский, П. Солин: Дискретный принцип максимума для конечных элементов высшего порядка в 1D, Math. Comput. 76 (2007), 1833–1846 гг.
  15. ^ L. Demkowicz, J. Kurtz, D. Pardo, W. Rachowicz, M. Paszynski, A. Zdunek: Computing with hp-Adaptive Finite Elements, Chapman & Hall / CRC Press, 2007
  16. ^ http://hpfem.org/wp-content/uploads/doc-web/doc-examples/src/hermes2d/examples/maxwell/microwave-oven.html
  17. ^ https://web.archive.org/web/20180807173436/http://hpfem.org/wp-content/uploads/doc-web/doc-examples/src/hermes2d/examples/maxwell/microwave-oven.html