УДК 519.632.4+519.635.1
DOI: 10.14529/ mmp190312
РЕШЕНИЕ ЭЛЛИПТИЧЕСКИХ УРАВНЕНИЙ В ПОЛИГОНАЛЬНЫХ ОБЛАСТЯХ МЕТОДОМ КОЛЛОКАЦИИ И НАИМЕНЬШИХ КВАДРАТОВ
В.П. Шапеее1'2, Л.С. Брындин1'2, В.А. Беляев1
1 Институт теоретической и прикладной механики им. С.А. Христиановича СО РАН, г. Новосибирск, Российская Федерация
2 Новосибирский государственный университет, г. Новосибирск, Российская Федерация
В данной работе рассматривается новый вариант метода коллокации и наименьших квадратов (КНК) для численного решения краевых задач для эллиптических уравнений в полигональных областях, в том числе в многосвязных. Возможности этого варианта и численные эксперименты рассмотрены на примерах решения уравнения Пуассона и неоднородного бигармонического уравнения. В качестве приложения решение неоднородного бигармонического уравнения использовано для моделирования напряженно-деформированного состояния (НДС) изотропной упругой тонкой пластинки полигональной формы, находящейся под действием поперечной нагрузки. Новый вариант метода КНК основан на триангуляции исходной области, чем принципиально отличается от предложенных ранее более сложных вариантов метода КНК решения краевых задач для уравнений с частными производными (УЧП) в нерегулярных областях. Установлено, что приближенные решения рассмотренных задач на последовательности измельчающихся сеток сходятся с повышением порядка и с высокой точностью совпадают с тестовыми решениями.
Ключевые слова: метод коллокации и наименьших квадратов; полигональная многосвязная область; уравнение Пуассона; неоднородное бигармоническое уравнение; напряженно-деформированное состояние.
Введение
Важным направлением вычислительной математики и прикладного математического моделирования является разработка новых численных методов решения краевых задач для УЧП в нерегулярных областях. Актуальность данного направления очевидна, поскольку многие явления в реальном мире происходят в областях со сложной геометрической формой. Поэтому востребована возможность моделирования таких явлений в областях с формой, соответствующей действительности, чтобы на практике в расчетах моделируемых явлений и физических процессов получались достоверные результаты.
Бигармонические уравнения имеют много приложений в разных областях науки и техники [1,2]. В данной работе в качестве приложения решение неоднородного бигармонического уравнения используется для моделирования НДС изотропной упругой тонкой пластинки полигональной формы в рамках теории тонких пластин [3]. На данный момент существует много разных численных методов решения бигармонического уравнения, сведения о которых можно найти в [1,2]. Однако, многие из них применимы только в канонических областях. Здесь также следует отметить, что популярна
идея решения краевых задач для УЧП в нерегулярных областях путем отображения их в каноническую область, для которой проще применить численный метод. О недостатках такого подхода указано в работе [4].
Также известно, что численное решение бигармонического уравнения, несмотря на его линейный характер, вызывает трудности из-за наличия производных четвертого порядка в дифференциальном уравнении. В этом случае имеет место плохая обусловленность исходной задачи и соответственно приближенной задачи, полученной тем или иным способом ее аппроксимации. Тем не менее, в случае шарнирного опирания пластинки многоугольной формы решение задачи сводится к последовательному решению системы двух уравнений Пуассона с краевыми условиями Дирихле [3,5]. В этом случае система линейных алгебраических уравнений (СЛАУ), к решению которой сводится численное решение задачи, оказывается обусловлена лучше, чем в первом случае. В случае других краевых условий различными авторами предлагались нетривиальные алгоритмы, позволяющее свести исходное бигармони-ческое уравнение к решению системы двух уравнений Пуассона [1,5]. Важным обстоятельством является то, что метод КНК способен численно решить с повышенной точностью исходное бигармоническое уравнение без сведения его к системе уравнений более низкого порядка.
За последнее время проекционно-сеточный метод КНК показал свою универсальность и эффективность при решении краевых задач для УЧП в областях с криволинейной границей [6], в выпуклых четырехугольных областях [7, 8], в произвольных треугольных областях [4], а также в областях с дискретно заданной границей [2,9]. В этих работах было показано преимущество применения регулярной сетки с прямоугольными ячейками внутри области и <одинарного> слоя нерегулярных ячеек (н-ячеек) около границы области. Для достижения аппроксимации повышенной точности дифференциальной задачи использовались «законтурные>, т.е. не принадлежащие исходной области решения дифференциальной задачи части граничных прямоугольных ячеек, отсеченных границей. Кроме того, «малые> по размеру (в сравнении с остальными) и/или вытянутые н-ячейки полагались «несамостоятельными> и присоединялись к соседним «самостоятельным> ячейкам.
В настоящей работе предлагается новый вариант метода КНК для решения уравнения Пуассона и неоднородного бигармонического уравнения в полигональных областях. Здесь в методе КНК для построения расчетной сетки применяется триангуляция исходной области. Излагаемый подход обладает рядом преимуществ и достоинств, среди которых в первую очередь следует выделить его применимость в случае невыпуклой и многосвязной полигональной области; относительную простоту реализации в сравнении с предложенными ранее вариантами метода КНК для нерегулярных областей.
1. Задача Дирихле для уравнения Пуассона 1.1. Постановка задачи и численный метод ее решения
Рассмотрим краевую задачу Дирихле для уравнения Пуассона
Ап = /(Х1,Х2), (Х1 ,Х2) е П, п|йп = д(х1,Х2) (1)
в полигональной области П с границей 8П (рис. 1 а), где А = + -щ, и(х\,х2) -искомое решение, /(х1,х2) и д(х1,х2) - заданные функции.
A гр g □.□ в
B
D
Рис.
C
а) б)
1. Полигональная двусвязная область решения задачи. Жирные линии -
стороны 12 «больших> треугольников. Пунктирные линии - стороны ячеек регулярных сеток в «больших> треугольниках (а). Фрагмент расчетной области, где • - точки коллокации; □ - точки согласования; х - точки записи краевых
условий (б)
Разобьем изначально исходную полигональную область П на М «больших> треугольников (на рис. 1 а) 12 «больших> треугольников). Далее в каждом «большом> треугольнике зададим регулярную сетку с М2 треугольными ячейками одинакового размера и подобными содержащему их «большому> треугольнику (рис. 1 а). Обозначим через Мсец3 = NN2 - количество расчетных ячеек в области. Те ячейки, ни одна сторона которых не лежит на границе области 8П, назовем внутренними, остальные - граничными. Введем локальную систему координат в каждой ячейке
°бЛаСТП (Л - -П,) (.Г2 - Щ) ,,,,
!/1 = —.!/2 = . (2)
где (х\з, х2з-) - центр пересечения медиан ]-ой треугольной ячейки (] = 1,..., Мсец3), 2Нц - длина стороны ячейки, к которой проведена высота размера 2^, £ =1,..., ь(уьу2) = м(х!(у1),х2(у2)). Задача (1) после замены (2) в локальных переменных примет вид
1
д2 v 1
+
д2 v
hit dVi2 h^t дУ22 у = #(xi(yi),x2(yi)),
f (xi (У1) , Х2(У2)) , (У1,У2) G
(У1,У2) G П, t =1,
Nt.
(3)
15
В ]-ой ячейке приближенное решение ьщ(уъу2) = Ф% задачи (3) ищем в ви-
г=1
де линейной комбинации с неопределенными коэффициентами базисных элементов пространства полиномов от двух переменных. В качестве базисных элементов взяты Фг(УъУ2) ^ {уГ1 УТ1 0 ^ а1 ^ 4, 0 ^ а2 ^ 4 — а1}, I = 1,..., 15. Глобальное приближенное решение состоит из «локальных> аналитических «кусков>, связанных с центром соответствующей ячейки. Для определения в каждой ячейке выписывается переопределенная «локальная> СЛАУ. В нее входят уравнения коллокации,
условия согласования и краевые условия на если ячейка является граничной. Совокупность СЛАУ во всех ячейках - глобальная СЛАУ определяет глобальное решение. Здесь число уравнений коллокации, условий согласования и краевых условий в каждой ячейке в сумме равно 32.
Уравнения коллокации, умноженные на h1th2t, в каждой j-ой ячейке выписываются в 16 точках коллокации (рис. 1 б) и имеют вид:
ЙЛ& + S50 = (4)
где t = 1,..., Nt, kc - весовой множитель уравнения коллокации в переопределенной СЛАУ при применении взвешенного метода наименьших квадратов.
В качестве условия согласования решения в нескольких точках согласования (рис. 1 б) на общих сторонах между соседними ячейками требуется непрерывность линейной комбинации с весовыми множителями значений решения и его производной по нормали dv
k-mo'Vhj кт 1 ^^ kmoVh кт 1 ^^ , (5)
где Uj - внешняя нормаль к границе j-ой ячейки; v^j и v^ - пределы приближенного решения при стремлении изнутри и извне к границе j-ой ячейки, соответственно; km0, kmi - весовые множители условий согласования.
Если сторона граничной j-ой ячейки принадлежит границе области, то в нескольких точках (рис. 1 б) на ней выпишем краевые условия
kbo vhj = kbo g (xi,X2), (6)
где kb0 - весовой множитель краевого условия Дирихле.
В каждой ячейке, объединяя уравнения (4)-(6), относительно brj получим переопределенную локальную СЛАУ [2,4,6-9,12]
15
Y^Birbrj = F, l = 1,..., 32, j = 1,...,Nceiis. (7)
r=1
Решение системы (7) на каждой итерации находится из условия минимума функционала Ф = J^1=1(Btrbrj — Ft)2 по коэффициентам brj, r =1,..., 15. Глобальная СЛАУ решается в итерационном процессе Гаусса - Зейделя, который продолжается до тех пор, пока не выполнится условие
max j — bn | <e, (8)
где b™ - r-й (r = 1,..., 15) коэффициент полинома, аппроксимирующего решение в ячейке с номером j на n-й итерации, величина е - заданная малая константа [2,4,69,12].
1.2. Результаты численных экспериментов
Во всех представленных ниже таблицах в пункте 1.2 приведены значения абсолютной погрешности приближенного решения в равномерной норме
Ec = || U — Uh ||с = max ( max |u(x1m,X2m) — Uhj(x1m,X2m)|), (9)
j=1,...NceUs m=1,...Lj
где Lj - количество равномерно распределенных контрольных точек (x1m,x2m), взятых в j-ой ячейке для подсчета в них погрешности, u - точное решение задачи (1),
Uh - приближенное решение в области П, u^j - приближенное решение в j-ой ячейке. Порядок сходимости погрешности приближенного решения вычисляется по формуле
E
R = bg2^, (10)
где Ec - значение погрешности на текущей сетке размера N2, Ep - значение погрешности на сетке размера N2/4. В численных экспериментах были взяты следующие значения параметров: Lj = 100, е = 10-14 в случае примеров 1-3 и е = 0, 5• 10-10 в случае примера 4, кс = 1, кто = 1, kmi = by/huh2t (t = 1 ,...,Nt), кь0 = 1,5. В качестве начального приближения решения взяты brj = 0, 4, где r = 1,..., 15, j = 1,..., Nceiis. Пусть Niter - количество итераций. Здесь в расчетах комбинируется операция продолжения (prolongation) в методе Федоренко и на всех сетках ускорение итерационного процесса по Крылову [10,11].
Пример 1. В табл. 1 приведены результаты численного эксперимента с тестовым решением u(x1,x2) = ex1 + ex2 в четырехугольной области с вершинами (0;0), (1;0), (0,8; 0,3), (-0,3; 0,9). Исходная область разбивалась на Nt = 2 «больших> треугольника.
Таблица 1
Результаты численных экспериментов
Пример 1 Пример 2
N cells Ec R Niter Ncells Ec R Niter
98 2,24E-8 - 81 392 7,40E-10 4,91 74 1568 3,00E-11 4,62 85 6272 1,85E-12 4,01 121 192 4,44E-8 - 78 768 l,46E-9 4,92 70 3072 4,70E-11 4,95 70 12288 2,38E-12 4,84 91
Пример 3 Пример 4
Ncells Ec R Niter Ncells Ec R Niter
400 3,73E-8 - 95 1600 1Д7Е-9 4,99 128 6400 5,51E-11 4,40 77 25600 3,49E-12 4,01 139 1536 9,89E-2 - 243 6144 5,98E-3 4,03 445 24576 3,69E-4 4,01 1085 98304 2,30E-5 4,00 2010
Пример 2. В табл. 1 приведены результаты численного эксперимента с тестовым решением u(xi, x2) = cos(xi) + cos(x2) в двусвязной области (рис. 1 a), вершины внешней границы: (cos sin (cos sin (cos Щ-] sin Щ-), (cos sin (eos sin
/ 5n • 5п \ ... /1 7n 1 7п \ /1 1 \
(eos sin вершины внутренней границы: cos j sin (j cos j sin (±cosff;±sinff), (I cos ff; I sin (I cos ff; I sin (I cos ^; I sin if). Здесь Nt = 12.
Пример 3. В табл. 1 приведены результаты численного эксперимента с тестовым решением u(x1,x2) = sin(x1) + sin(x2) в невыпуклой шестиугольной области с вершинами (0; 0), (0,5; 1,25), (1,5; 1,5), (1; 2,25), (2,5; 2), (2; 0,5). Здесь Nt = 4.
Пример 4. В табл. 1 приведены результаты численного эксперимента с тестовым решением u(x1,x2) = e10xi + e10x2 + e5xi+5x2, имеющим весьма большие градиенты в шестиугольной области с вершинами, которые совпадают с вершинами внешней границы из примера 2. Здесь Nt = 6.
Из приведенных результатов видно, что порядок сходимости погрешности приближенного решения не хуже четвертого, что обеспечено высокой степенью базисных полиномов в рассматриваемом функциональном пространстве. Такой порядок сходимости ожидался и он не был понижен несмотря на то, что в коэффициент остаточного члена аппроксимации решения, естественно, входят большие значения производных. Однако в этом примере величина погрешности больше по сравнению со случаем решений задач без особенностей, что и следовало ожидать.
Пример 5. Рассмотрим пример задачи Дирихле для уравнения Пуассона (1) с f (хьх2) = 1, д(хьх2)|гп = 0 [9] в шестиугольной области из примера 4. В этом примере из-за наличия разрыва вторых производных в точках излома границы имеет место особенность в решении задачи (см. раб. [9] и ссылки внутри). Это является причиной увеличения погрешности численного решения задачи и снижения порядка его сходимости, значение которого существенно зависит от рассматриваемой подобласти, взятой в П (табл. 2). Для оценки порядка сходимости численного решения в равномерной норме вычислялась величина
Е = I ^сМв — и^еш/4 ||с = 1имсе11, (х1? ,х2?) - имсеШ/4(ХЧ ,х2?)1, (11)
где Q - количество равномерно распределенных контрольных точек (х^ , х2^) в области; им1,и2 - численное решение, полученное на сетке, содержащей Ысец3 ячеек, имсе11а/4 - численное решение на сетке, содержащей Мсец3/4 ячеек. Порядок сходимости решения формально определялся по формуле
^ т Е„ 1 имсМз — имсе11з/4 Нс
К- = 1°§2 = 1°§2 7Г- --> (12)
Ес Н UNceЛls/4 — иМсе11в/16 Н|С
где ЕР =Н иМсеш — и-^е113/4 Н|С, Ес =Н иМсеш/4 — иМсец,/16 Н|С и иМсец,/16 - численное
решение на сетке, содержащей N^/16 ячеек. В табл. 2 приведены полученные в численном эксперименте значения нормы (11) и соответствующие порядки ее сходимости (12) при Q = 200Мсец3 во всей исходной области П ив подобласти П1: х2 + х2 ^ 0, 52, которая не содержит точек излома границы.
Таблица 2
Результаты численных экспериментов в задаче с разрывом второй производной
II - и^е11е/2 IIс Щп К Е пх К
II «1536 - «384 ||с 6,84Е-4 II «6144-«1536 ||с 2,41Е-4 1,50 «24576 - «6144 с 8,55Е-5 1,49 «98304 - «24576 С 3,02Е-5 1,50 8Д9Е-6 9,93Е-7 3,04 1,22Е-7 3,02 1,52Е-8 3,00
2. Неоднородное бигармоническое уравнение 2.1. Постановка задачи и описание метода
Пусть требуется определить НДС упругой изотропной пластинки полигональной формы, закрепленной на краях и находящейся под действием поперечной нагрузки произвольного вида. Будем считать, что нагрузки достаточно малы, и деформирование пластинки является упругим. Пластинка предполагается тонкой, что позволяет свести трехмерную задачу теории упругости к двумерной задаче теории тонких
пластин. Согласно [3], разрешающую систему уравнений запишем в перемещениях, отыскав которые можно полностью рассчитать НДС пластинки. Величины прогибов пластинки определяются из решения неоднородного бигармонического уравнения
А2и> = (Жьх2)бО, (13)
где ^(х1,х2) - прогиб срединной поверхности пластинки, д(х1 ,х2) - поперечная нагрузка, Д = Ек3/(12(1 — V2)) - жесткость пластинки при изгибе, к - толщина пластинки, Е, V - модуль Юнга и коэффициент Пуассона изотропного материала пластинки, полигональная область П - проекция срединной поверхности пластинки на плоскость (х1, х2). В этой задаче на каждом куске границы пластинки можно задать в частности одно из следующих условий ее закрепления:
дт
и> = О, -г— = 0 - защемленный край, (14)
дщ
т = 0, Мпт = 0 - шарнирно закрепленный край, (15)
Мпт = 0, УПт = 0 - свободный край, (16)
где - внешняя нормаль к границе области дП, Мп - изгибающий момент, Уп - обобщенная поперечная сила [3]. К решению поставленной задаче применим алгоритм, аналогичный описанному выше. Отличия будут в записи уравнений приближенной задачи. Аналогично уравнениям (4), (5) и (6) в каждой ]-ой ячейке выписываются 1) уравнения коллокации, умноженные на к^к^, в 16 точках коллокации (рис. 1 б)
= ь*2), (17)
е ду\ ду>дуI ад ) '
где кс - весовой множитель уравнений коллокации, £ = 1,...,
2) условия согласования (умноженные на к14к24 во втором случае) на каждой общей стороне между соседними ячейками в нескольких точках согласования (рис. 1 б)
ь 1)и. + Ь - к ■,%+]<■ ^ Г1Я1
^то ип] I "тех ^^ "тоо ип I "тг ■> у-1-")
ПиЫг (кт2^| + ктз^) = НиЬъ {кш^ + , (19)
где щ - внешняя нормаль к границе ]-ой ячейки, £ = 1,..., г^ и г^ - пределы приближенного решения задачи при стремлении изнутри и извне к границе ]-й ячейки, соответственно, кто, кт1, кт2, ктз - весовые множители условий согласования;
3) краевые условия в нескольких точках записи краевых условий (рис. 1 б) на границе области, если ячейка является граничной.
В данном случае каждая <локальная> СЛАУ состоит из 48 уравнений: 16 уравнений коллокации и по два уравнения в каждой точке для записи условий согласования или краевых условий в случае граничной ячейки.
2.2. Результаты численных экспериментов
Пример 6. Рассмотрим шарнирно закрепленную на краях (15) прямоугольную пластинку размера ^ х , находящуюся под действием поперечной нагрузки д = 105 в1п(пх1/^1) вт(пх2/^2) (Па). В этом случае задача имеет точное решение [3]
w(x 1,х2) =-:L-LJ-(20)
V ^ nW (dj + dlf К '
Ранее эта задача была решена методом КНК в работе [12] с использованием пространства полиномов четвертой степени. В случае краевых условий (15) решение бигармонического уравнения (13) в полигональной области можно свести к последовательному решению системы двух уравнений Пуассона с краевыми условиями Дирихле [3]. Как отмечено во введении, задача Дирихле для уравнения Пуассона -хорошо обусловленная задача по сравнению с краевой задачей для бигармоническо-го уравнения, что, очевидно, сказывается на точности численного решения задачи. Пусть теперь Ec - относительная погрешность численного решения, порядок сходимости R которой вычислялся по формуле (10). В табл. 3 приведено сравнение результатов численного решения уравнения [12] и результатов решения системы двух уравнений Пуассона. Здесь прямоугольная область покрывалась регулярной сеткой с прямоугольными ячейками. В этом случае был использован вариант 1 метода КНК из работы [8]. В расчетах использовались следующие параметры: d1 = 10 м, d2 = 10 м, h = 0,1 м, E = 200 ГПа, v = 0, 28, е = 10-14. В табл. 3 в Niter в скобках указано количество итераций для первого уравнения Пуассона (аналогично в табл. 4). В работе [12] в случае Nceas = 12800 нет информации о количестве итераций.
Таблица 3
Сравнение результатов численных экспериментов
Пример 1, работа [12] Пример 1, эта работа
Ncells Ее R Niter Ес R Niter
100 400 1600 6400 12800 1,82Е-2 - 129 2,83Е-3 2,68 282 4,56Е-4 2,63 593 7,57Е-5 2,58 2090 1,32Е-5 2,51 1,00Е-4 - 68 (58) 6,32Е-6 3,98 81 (74) 3,95Е-7 4,00 112 (123) 2,47Е-8 3,99 162 (135) 1,57Е-9 3,97 251 (161)
Таблица 4
Сравнение результатов численных экспериментов
Пример 2, случай 1 Пример 2, случай 2
Ncells Ес R Nuer Ес R Nuer
48 192 768 3072 12288 3,97Е-4 - 45 3,74Е-5 3,40 43 6,24Е-6 2,58 81 1Д5Е-6 2,43 241 1,29Е-4 - 81 (81) 4,43Е-6 4,86 81 (81) 1,37Е-7 5,01 86 (81) 4,63Е-9 4,88 97 (85) 2,48Е-10 4,22 174 (161)
Пример 7. Рассмотрим полигональную пластинку с вырезом из примера 2 пункта 1.2, при этом исходная область растягивалась в 5 раз и сдвигалась на вектор (5; 5). Зададим на границе 80 значения т и Мп/ш, вычисленные из точного решения (20). В табл. 4 приведено сравнение результатов численного решения исходного бигармонического уравнения (случай 1) и результатов численного решения системы двух уравнений Пуассона (случай 2) в пластинке с вырезом. В расчетах использовались следующие параметры: = 9, 65 м, = 9, 65 м, к = 0,1 м, Е = 200 ГПа, V = 0, 28, е = 10-10 для численного решения бигармонического уравнения, е = 10-14 для численного решения системы двух уравнений Пуассона. Здесь N = 12. Здесь
для краткости значения d1 и d2 округлены до сотых. Значения весовых параметров в первом случае взяты из работы [2], а во втором случае из пункта 1.2. В первом случае при Nce11s = 12288 значение погрешности численного решения не изменялось по сравнению с предыдущим случаем.
Пример 8. Рассмотрим два примера плохообусловленных задач, точное решение которых не известно. Пусть в первом случае шестиугольная пластинка, а во втором случае эта же пластинка, но с вырезом, находятся под действием равномерной поперечной нагрузки q = const. В первом случае два края пластинки защемлены (14), два - шарнирно закреплены (15), два - свободны (16). Во втором случае внутренние края пластинки шарнирно закреплены (15), внешние четыре края - защемлены (14), два - свободны (16). Напряженное состояние пластины описывается функцией распределения интенсивности напряжений
I=^(ax-ayy + al + ^y + 6axy^ (21)
где <jx = -E(wxx + vwyy), (Jy = —E(wyy + vwxx), аху = -E(wxy), Е = Ех3/( 1 - и2), x3 Е [-h/2, h/2]. Отметим, что здесь относительно легко вычислить любые производные приближенного решения, поскольку в методе КНК оно является кусочно-полиномиальным. Это одно из достоинств используемого метода. В обоих случаях в расчетах использовались следующие параметры задачи: h = 0,1 м, E = 200 ГПа, V = 0, 28, x3 = h/2, б = 10-10. В первом примере q = 200 КПа, NceUs = 1536, Nt = 6, а во втором q = 1 МПа, Nce11s = 3072, Nt = 12. На рис. 2, 3 показаны форма прогиба w a) и интенсивность напряжений I (21) б).
Рис. 2. Величина прогиба w (a) деформированной пластинки (два края защемлены, два - шарнирно закреплены, два - свободны) и интенсивность напряжений I (б)
Заключение
Предложен и реализован новый вариант метода КНК для численного решения эллиптических уравнений в полигональных областях, в том числе в многосвязных. Проведен ряд численных экспериментов, в которых показано, что для задач с достаточно гладкими решениями, в том числе с большими градиентами, имеет место
повышенный порядок сходимости к нулю погрешности решения. В этом случае при решении уравнения Пуассона наблюдался порядок сходимости не ниже четвертого, а при решении бигармонического уравнения не меньше второго.
10 8 6 4 2 0
Ж
%
э
I, ГПа
1.48 1.11 0.74 0.37
46
xi, м
8 10
а) б)
Рис. 3. Величина прогиба эд (а) деформированной пластинки с вырезом (внутренние края пластинки шарнирно закреплены, внешние четыре края -защемлены, два - свободны) и интенсивность напряжений I (б)
Здесь также рассмотрен пример задачи Дирихле для уравнения Пуассона, когда на границе области в точках ее излома имеет место особенность - разрыв вторых производных искомого решения. В этом случае порядок сходимости численного глобального решения задачи понижается и был равен полтора в отличие от случаев, когда краевые условия и правая часть уравнения «согласованы» и достаточно гладкие. Но за исключением окрестности особых точек на границе, порядок сходимости решения остается повышенным и равен трем. Это хорошо согласуется с теорией краевых задач для эллиптических уравнений.
Кроме того, в работе приведены результаты решения задач расчета НДС изотропных пластин сложной формы. Решение таких задач имеет самостоятельное прикладное значение.
Работа выполнена в рамках Программы фундаментальных научных исследований государственных академий наук на 2013-2020 годы (проект № АААА-А17-117030610136-3).
2
Литература
1. Guo Chen. A Fast Finite Difference Method for Biharmonic Equations on Irregular Domains and Its Application to an Incompressible Stokes Flow / Guo Chen, Zhilin Li, Ping Lin // Advances in Computational Mathematics. - 2008. - V. 29, № 2. - P. 113-133.
2. Шапеев, В.П. Решение с повышенной точностью бигармонического уравнения в нерегулярных областях методом коллокации и наименьших квадратов /В.П. Шапеев, В.А. Беляев // Вычислительные методы и программирование. -2018. - Т. 19, №4. - С. 340-355.
3. Тимошенко, С.П. Пластины и оболочки / С.П. Тимошенко, С. Войновский-Кригер. -М.: Физматгиз, 1963.
4. Шапеев, В.П. Решение краевых задач для уравнений с частными производными в треугольных областях методом коллокации и наименьших квадратов / В.П. Шапеев,
B.А. Беляев // Вычислительные методы и программирование. - 2018. - Т. 19, № 1. -
C. 96-111.
5. Сорокин, С.Б. Переобусловливание при численном решении задачи Дирихле для бигар-монического уравнения / С.Б. Сорокин // Сибирский журнал вычислительной математики. - 2011. - Т. 14, № 2. - С. 205-213.
6. Шапеев, В.П. Варианты метода коллокации и наименьших невязок повышенной точности в области с криволинейной границей / В.П. Шапеев, В.А. Беляев // Вычислительные технологии. - 2016. - Т. 21, № 5. - С. 95-110.
7. Беляев, В.А. Варианты метода коллокации и наименьших невязок для решения задач математической физики в трапециевидных областях / В.А. Беляев, В.П. Шапеев // Вычислительные технологии. - 2017. - Т. 22, № 4. - С. 22-42.
8. Беляев, В.А. Варианты метода коллокации и наименьших невязок для решения задач математической физики в выпуклых четырехугольных областях / В.А. Беляев, В.П. Шапеев // Моделирование и анализ информационных систем. - 2017. - Т. 24, № 5. - С. 629-648.
9. Беляев, В.А. Решение задачи Дирихле для уравнения Пуассона методом коллокации и наименьших квадратов в области с дискретно заданной границей / В.А. Беляев, В.П. Шапеев // Вычислительные технологии. - 2018. - Т. 23, № 3. - С. 15-30.
10. Федоренко, Р.П. О скорости сходимости одного итерационного процесса / Р.П. Федорен-ко // Журнал вычислительной математики и математической физики. - 1964. - Т. 4, № 3. - С. 559-564.
11. Saad, Y. Numerical Methods for Large Eigenvalue Problems / Y. Saad. - Manchester: Manchester University Press, 1991.
12. Голушко, С.К. Метод коллокаций и наименьших невязок в приложении к задачам механики изотропных пластин / С.К. Голушко, С.В. Идимешев, В.П. Шапеев // Вычислительные технологии. - 2013. - Т. 18, № 6. - С. 31-43.
Василий Павлович Шапеев, доктор физико-математических наук, главный научный сотрудник, Институт теоретической и прикладной механики им. С.А. Хри-стиановича СО РАН (г. Новосибирск, Российская Федерация); профессор, кафедра «Математическое моделирование:», Новосибирский государственный университет (г. Новосибирск, Российская Федерация), shapeev.vasily@mail.ru.
Лука Сергеевич Брындин, лаборант, Институт теоретической и прикладной механики им. С.А. Христиановича СО РАН (г. Новосибирск, Российская Федерация); магистрант механико-математического факультета, кафедра «Математическое моделирование», Новосибирский государственный университет (г. Новосибирск, Российская Федерация), bryndin-1996@mail.ru.
Василий Алексеевич Беляев, аспирант, младший научный сотрудник, Институт теоретической и прикладной механики им. С.А. Христиановича СО РАН (г. Новосибирск, Российская Федерация), belyaevasily@mail.ru.
Поступила в редакцию 23 ноября 2018 г.
MSC 35J15, 35J40, 41A10, 65N50 DOI: 10.14529/mmp190312
SOLVING ELLIPTIC EQUATIONS IN POLYGONAL DOMAINS BY THE LEAST SQUARES COLLOCATION METHOD
V.P. Shapeev1'2, L.S. Bryndin1'2, V.A. Belyaev1
1 Khristianovich Institute of Theoretical and Applied Mechanics, Siberian Branch
of Russian Academy of Sciences, Novosibirsk, Russian Federation
2 Novosibirsk State University, Novosibirsk, Russian Federation
E-mail: shapeev.vasily@mail.ru, bryndin-1996@mail.ru, belyaevasily@mail.ru
The paper considers a new version of the least squares collocation (LSC) method for the numerical solution of boundary value problems for elliptic equations in polygonal domains, in particular, in multiply connected domains. The implementation of this approach and numerical experiments are performed on the examples of the inhomogeneous biharmonic and Poisson equations. As an application, we use the nonhomogeneous biharmonic equation to simulate the stress-strain state of isotropic elastic thin plate of polygonal form under the action of transverse load. The new version of the LSC method is based on the triangulation of the original domain. Therefore, this approach is fundamentally different from the previous more complicated versions of the LSC method proposed to solve the boundary value problems for partial derivative equations in irregular domains. We make the numerical experiments on the convergence of the approximate solution to various problems on a sequence of grids. The experiments show that the solution to the problems converges with high order and, in the case of the known analytical solution, matches with high accuracy with the analytical solution to the test problems.
Keywords: least squares collocation method; polygonal multiply connected domain; Poisson's equation; nonhomogeneous biharmonic equation; stress-strain state.
References
1. Guo Chen, Zhilin Li, Ping Lin. A Fast Finite Difference Method for Biharmonic Equations on Irregular Domains and its Application to an Incompressible Stokes Flow. Advances in Computational Mathematics, 2008, vol. 29, no. 2, pp. 113-133. DOI: 10.1007/s10444-007-9043-6
2. Shapeev V.P., Belyaev V.A. [Solving the Biharmonic Equation with High Order Accuracy in Irregular Domains by the Least Squares Collocation Method]. Numerical Methods and Programming, 2018, vol. 19, no. 4, pp. 340-355. (in Russian)
3. Timoshenko S.P., Woinowsky-Krieger S. Theory of Plates and Shells. N.-Y., Toronto, London, McGraw-Hill, 1959.
4. Shapeev V.P., Belyaev V.A. Solving Boundary Value Problems for Partial Differential Equations in Triangular Domains by the Least Squares Collocation Method. Numerical Methods and Programming, 2018, vol. 19, no. 1, pp. 96-111. (in Russian)
5. Sorokin, S.B. Preconditioning in the Numerical Solution to a Dirichlet Problem for the Biharmonic Equation. Numerical Analysis and Applications, 2011, vol. 4, no. 2, pp. 167-174. DOI: 10.1134/S1995423911020078
6. Shapeev V.P., Belyaev V.A. Versions of High Order Accuracy Collocation and Least Residuals Method in the Domain with a Curvilinear Boundary. Computational Technologies, 2016, vol. 21, no. 5, pp. 95-110. (in Russian)
B.n. fflaneeB, ^.C. BpuHflHH, B.A. Be^HeB
7. Belyaev V.A., Shapeev V.P. The Versions of Collocation and Least Residuals Method for Solving Problems of Mathematical Physics in the Trapezoidal Domains. Computational Technologies, 2017, vol. 22, no. 4, pp. 22-42. (in Russian)
8. Belyaev V.A., Shapeev V.P. Versions of the Collocation and Least Residuals Method for Solving Problems of Mathematical Physics in the Convex Quadrangular Domains. Modeling and Analysis of Information Systems, 2017, vol. 24, no. 5, pp. 629-648. DOI: 10.18255/18181015-2017-5-629-648 (in Russian)
9. Belyaev V.A., Shapeev V.P. Solving the Dirichlet Problem for the Poisson Equation by the Least Squares Collocation Method with Given Discrete Boundary Domain. Computational Technologies, 2018, vol. 23, no. 3, pp. 15-30. DOI: 10.25743/ICT.2018.3.15956 (in Russian)
10. Fedorenko R.P. The Speed of Convergence of One Iterative Process. Computational Mathematics and Mathematical Physics, 1964, vol. 4, no. 3, pp. 227-235. DOI: 10.1016/0041-5553(64)90253-8
11. Saad, Y. Numerical Methods for Large Eigenvalue Problems. Manchester, Manchester University Press, 1991.
12. Golushko S.K., Idimeshev S.V., Shapeev V.P. Application of Collocations and Least Residuals Method to Problems of the Isotropic Plates Theory. Computational Technologies, 2013, vol. 18, no. 6, pp. 31-43. (in Russian)
Received November 23, 2018