Jump to content

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

(Перенаправлено с элемента «Граница» )

Метод граничных элементов ( МГЭ ) — это численный вычислительный метод решения линейных дифференциальных уравнений в частных производных , которые были сформулированы как интегральные уравнения (т.е. в граничной интегральной форме), включая механику жидкости , акустику , электромагнетизм (где этот метод известен как метод моментов или сокращенно МОМ ), [1] механика разрушения , [2] и обратитесь к механикам . [3] [4]

Математическая основа

[ редактировать ]

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

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

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

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

Вывод функций Грина в замкнутой форме представляет особый интерес в методе граничных элементов, особенно в электромагнетике. В частности, при анализе слоистых сред вывод функции Грина в пространственной области требует обращения аналитически выводимой функции Грина в спектральной области через интеграл по траектории Зоммерфельда. Этот интеграл не может быть вычислен аналитически, а его численное интегрирование является дорогостоящим из-за его колебательного и медленно сходящегося поведения. Для надежного анализа пространственные функции Грина аппроксимируются комплексными экспонентами с помощью таких методов, как метод Прони или обобщенный пучок функций , а интеграл оценивается с помощью тождества Зоммерфельда . [5] [6] [7] [8] Этот метод известен как метод дискретного комплексного изображения. [7] [8]

Сравнение с другими методами

[ редактировать ]

Метод граничных элементов часто более эффективен, чем другие методы, включая конечные элементы, с точки зрения вычислительных ресурсов для задач, где имеется небольшое соотношение поверхности/объема. [9] Концептуально, это работает путем создания « сетки » поверх моделируемой поверхности. Однако для многих задач методы граничных элементов существенно менее эффективны, чем методы объемной дискретизации ( метод конечных элементов , метод конечных разностей , метод конечных объемов ). Хорошим примером применения метода граничных элементов является эффективный расчет собственных частот в плесков жидкости резервуарах. [10] [11] [12] Метод граничных элементов является одним из наиболее эффективных методов численного моделирования контактных задач. [13] в частности для моделирования клеевых контактов. [14]

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

См. также

[ редактировать ]
  1. ^ В электромагнетике более традиционный термин «метод моментов» часто используется, хотя и не всегда, как синоним «метода граничных элементов»: см. ( Gibson 2008 ). дополнительную информацию по этому вопросу
  2. ^ Метод граничных элементов хорошо подходит для анализа трещин в твердых телах. Существует несколько подходов к решению проблем с граничными элементами. Одним из таких подходов является формулировка условий на трещины в терминах гиперсингулярных граничных интегральных уравнений, см. ( Ang 2013 ).
  3. ^ Порт, Р.; Ли, К. (01 октября 2014 г.). «Полная формулировка граничных элементов для задач нормального и тангенциального контакта». Физическая мезомеханика . 17 (4): 334–340. дои : 10.1134/S1029959914040109 . ISSN   1029-9599 . S2CID   137494525 .
  4. ^ «Учебное пособие по расчету контактного давления на основе BEM» . www.tribonet.org . 9 ноября 2017 г.
  5. ^ Чоу, ЮЛ; Ян, Джей-Джей; Фанг, Д.Г.; Ховард, GE (март 1991 г.). «Пространственная функция Грина замкнутой формы для толстой микрополосковой подложки». Транзакции IEEE по теории и технике микроволнового излучения . 39 (3): 588–592. Бибкод : 1991ITMTT..39..588C . дои : 10.1109/22.75309 .
  6. ^ Аксун, Мичиган (февраль 2003 г.). «Надежный подход к выводу функций Грина в замкнутой форме». Транзакции IEEE по теории и технике микроволнового излучения . 44 (5): 651–658. дои : 10.1109/22.493917 . hdl : 11693/10779 .
  7. ^ Jump up to: а б Тео, Суи-Энн (2000). «Метод дискретного комплексного изображения для функций Грина общих многослойных сред». IEEE СВЧ и волноводные буквы . 10 (10): 400–402. дои : 10.1109/75.877225 .
  8. ^ Jump up to: а б Тео, Суи-Энн; Чу, Сиу-Тек; Леонг, Мук-Сенг (февраль 2003 г.). «Анализ ошибок метода дискретных комплексных изображений и извлечения полюсов». Транзакции IEEE по теории и технике микроволнового излучения . 51 (2): 406–412. Бибкод : 2003ITMTT..51..406T . дои : 10.1109/TMTT.2002.807834 .
  9. ^ Это ( Кацикаделис 2002 ).
  10. ^ Колаи, Амир; Ракхеджа, Субхаш; Ричард, Марк Дж. (01 сентября 2015 г.). «Трехмерное динамическое выплескивание жидкости в частично заполненных горизонтальных резервуарах, подверженных одновременным продольным и боковым возбуждениям». Европейский журнал механики Б. 53 : 251–263. Бибкод : 2015EuJMB..53..251K . doi : 10.1016/j.eurotechflu.2015.06.001 .
  11. ^ Колаи, Амир; Ракхеджа, Субхаш; Ричард, Марк Дж. (31 января 2015 г.). «Совместный мультимодальный метод и метод граничных элементов для анализа эффективности защиты от выплескивания частичных перегородок в частично заполненном контейнере». Компьютеры и жидкости . 107 : 43–58. дои : 10.1016/j.compfluid.2014.10.013 .
  12. ^ Колаи, Амир; Ракхеджа, Субхаш; Ричард, Марк Дж. (14 ноября 2014 г.). Том 4А: Динамика, вибрация и управление . стр. V04AT04A067. дои : 10.1115/IMECE2014-37271 . ISBN  978-0-7918-4647-6 .
  13. ^ Попов, Валентин (2017). Контактная механика и трение – Физические основы и (глава 19) . Спрингер. стр. 337–341. ISBN  9783662530801 .
  14. ^ Порт, Роман; Попов, Валентин Л. (9 апреля 2015 г.). «Моделирование адгезионного контакта упругих тел с использованием локального критерия отрыва, зависящего от сетки, в методе граничных элементов» . Facta Universitatis, Серия: Машиностроение . 13 (1): 3–10.
  15. ^ Наджарзаде, Л., Мовахедян, Б. и Ажари, М., 2022. Численное решение задач распространения водных волн в переменных батиметриях с использованием модифицированного метода граничных элементов радиального интегрирования. Океанское машиностроение, 257, стр.111613.
  16. ^ Наджарзаде Л., Мовахедян Б. и Ажари М., 2019. Численное решение скалярного волнового уравнения с помощью модифицированного метода граничных элементов радиального интегрирования. Инженерный анализ с граничными элементами, 105, стр. 267–278.

Библиография

[ редактировать ]

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

[ редактировать ]
  • Констанда, Кристиан; Доти, Дейл; Хэмилл, Уильям (2016). Методы граничных интегральных уравнений и численные решения: тонкие пластины на упругом основании . Нью-Йорк: Спрингер. ISBN  978-3-319-26307-6 .
[ редактировать ]

Бесплатное программное обеспечение

[ редактировать ]
  • Bembel Трехмерное изогеометрическое программное обеспечение БЭМ более высокого порядка с открытым исходным кодом для задач Лапласа, Гельмгольца и Максвелла, использующее метод быстрых мультиполей для сжатия и снижения вычислительных затрат.
  • border-element-method.com Программное обеспечение БЭМ с открытым исходным кодом для решения задач акустики / Гельмгольца и Лапласа.
  • Puma-EM Высокопроизводительная параллельная программа «Метод моментов» с открытым исходным кодом / Многоуровневый быстрый многополюсный метод.
  • AcouSTO Acoustics Simulation Tool, бесплатный параллельный БЭМ-решатель с открытым исходным кодом для интегрального уравнения Кирхгофа-Гельмгольца (KHIE)
  • FastBEM Бесплатные программы быстрых мультипольных граничных элементов для решения 2D/3D задач потенциала, упругости, стоксова течения и акустических задач.
  • ParaFEM Включает в себя бесплатный параллельный БЭМ-решатель с открытым исходным кодом для задач упругости, описанный в Герноте Бире, Яне Смите, Кристиане Дуэнсере, « Метод граничных элементов с программированием: для инженеров и ученых» , Springer, ISBN   978-3-211-71574-1 (2008 г.)
  • Библиотека шаблонов граничных элементов (BETL). Программная библиотека C++ общего назначения для дискретизации граничных интегральных операторов.
  • Nemoh Программное обеспечение BEM для гидродинамики с открытым исходным кодом, предназначенное для расчета волновых нагрузок первого порядка на морские сооружения (добавленная масса, радиационное демпфирование, силы дифракции).
  • Bempp , программное обеспечение БЭМ с открытым исходным кодом для решения 3D-задач Лапласа, Гельмгольца и Максвелла.
  • MNPBEM , набор инструментов Matlab с открытым исходным кодом для решения уравнений Максвелла для наноструктур произвольной формы.
  • Симулятор контактной механики и трибологии , бесплатное программное обеспечение на базе BEM
  • MultiFEBE , решатель BEM-FEM для вычислительной механики, позволяющий связывать 2D и 3D вязкоупругие или пороупругие среды с элементами конструкции балки и оболочки (например, для задач динамического взаимодействия грунта и конструкции).
  • BE-STATIK Бесплатные BE-программы для двумерных задач потенциала, упругости и изгиба пластин (Кирхгофа).
Arc.Ask3.Ru: конец переведенного документа.
Arc.Ask3.Ru
Номер скриншота №: daa3588a11a1f04ab1b268c8e5f8e0a1__1714692240
URL1:https://arc.ask3.ru/arc/aa/da/a1/daa3588a11a1f04ab1b268c8e5f8e0a1.html
Заголовок, (Title) документа по адресу, URL1:
Boundary element method - Wikipedia
Данный printscreen веб страницы (снимок веб страницы, скриншот веб страницы), визуально-программная копия документа расположенного по адресу URL1 и сохраненная в файл, имеет: квалифицированную, усовершенствованную (подтверждены: метки времени, валидность сертификата), открепленную ЭЦП (приложена к данному файлу), что может быть использовано для подтверждения содержания и факта существования документа в этот момент времени. Права на данный скриншот принадлежат администрации Ask3.ru, использование в качестве доказательства только с письменного разрешения правообладателя скриншота. Администрация Ask3.ru не несет ответственности за информацию размещенную на данном скриншоте. Права на прочие зарегистрированные элементы любого права, изображенные на снимках принадлежат их владельцам. Качество перевода предоставляется как есть. Любые претензии, иски не могут быть предъявлены. Если вы не согласны с любым пунктом перечисленным выше, вы не можете использовать данный сайт и информация размещенную на нем (сайте/странице), немедленно покиньте данный сайт. В случае нарушения любого пункта перечисленного выше, штраф 55! (Пятьдесят пять факториал, Денежную единицу (имеющую самостоятельную стоимость) можете выбрать самостоятельно, выплаичвается товарами в течение 7 дней с момента нарушения.)