Научная статья на тему 'Разработка и верификация метода коллокации и наименьших квадратов с полиномами седьмой степени для решения бигармонического уравнения'

Разработка и верификация метода коллокации и наименьших квадратов с полиномами седьмой степени для решения бигармонического уравнения Текст научной статьи по специальности «Математика»

CC BY
12
3
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
метод коллокации и наименьших квадратов / кусочно-полиномиальный базис / автоматическая склейка решения / бигармоническое уравнение / изгиб пластины / least-squares collocation method / piecewise polynomials / automatic solution continuity / biharmonic equation / plate bending

Аннотация научной статьи по математике, автор научной работы — Лука Сергеевич Брындин

Предложен и реализован новый сеточный вариант метода коллокации и наименьших квадратов с полиномами седьмой степени (C3-МКНК7), обладающий непрерывностью вплоть до третьих производных кусочно-полиномиального решения в смысле наименьших квадратов. C3-МКНК7 принципиально отличается от предыдущих версий МКНК отсутствием условий согласования, в которых в явном виде в нескольких точках на границах между соседними ячейками требовалась непрерывность решения и его производных. Приведены значения точности и порядка сходимости C3-МКНК7 при численном решении двумерных краевых задач для бигармонического уравнения в квадратной области и области с криволинейной границей. Показаны преимущества C3-МКНК7 над предыдущими вариантами МКНК.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Лука Сергеевич Брындин

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Development and verification of the least-squares collocation method with seventh degree polynomials for the biharmonic equation

A new version of the leastsquares collocation method (C3-LSCM7) with seventh degree polynomials is proposed and implemented. The method has continuity up to the third derivatives of the piecewise polynomial solution in the sense of least squares. This is achieved by using the values of the solution and its derivatives at the vertices of the grid cells as unknown terms. The C3-LSCM7 is fundamentally different from previous versions of the LSCM in the absence of matching conditions. They explicitly require continuity of the solution and its derivatives at several points on the boundaries between neighboring cells. To solve the problem in a domain with a curved boundary a grid is constructed with rectangular cells. The solution of small irregular cells is continued from neighboring independent ones. Verification of the C3-LSCM7 is carried out by solving two-dimensional boundary value problems for a biharmonic equation in a square and in domain with the curvilinear boundary. The condition numbers of a global matrix and transition matrices from values at nodes to coefficients of polynomial expansion are studied. The advantages of the C3-LSCM7 over previous versions of the LSCM are shown.

Текст научной работы на тему «Разработка и верификация метода коллокации и наименьших квадратов с полиномами седьмой степени для решения бигармонического уравнения»

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

2024 Математика и механика № 91

Tomsk State University Journal of Mathematics and Mechanics

МАТЕМАТИКА

М ATHEMATICS

Научная статья

УДК 519.635.1 MSC: 65N35, 65N50

doi: 10.17223/19988621/91/1

Разработка и верификация метода коллокации и наименьших квадратов с полиномами седьмой степени для решения бигармонического уравнения

Лука Сергеевич Брындин

Институт теоретической и прикладной механики им. С.А. Христиановича Сибирского отделения РАН, Новосибирск, Россия, [email protected]

Аннотация. Предложен и реализован новый сеточный вариант метода коллокации и наименьших квадратов с полиномами седьмой степени (С3-МКНК7), обладающий непрерывностью вплоть до третьих производных кусочно-полиномиального решения в смысле наименьших квадратов. С3-МКНК7 принципиально отличается от предыдущих версий МКНК отсутствием условий согласования, в которых в явном виде в нескольких точках на границах между соседними ячейками требовалась непрерывность решения и его производных. Приведены значения точности и порядка сходимости С3-МКНК7 при численном решении двумерных краевых задач для бигармонического уравнения в квадратной области и области с криволинейной границей. Показаны преимущества С3-МКНК7 над предыдущими вариантами МКНК.

Ключевые слова: метод коллокации и наименьших квадратов, кусочно-полиномиальный базис, автоматическая склейка решения, бигармоническое уравнение, изгиб пластины

Благодарности: Работа выполнена в рамках государственного задания ИТПМ СО РАН.

Для цитирования: Брындин Л.С. Разработка и верификация метода коллокации и наименьших квадратов с полиномами седьмой степени для решения бигармонического уравнения // Вестник Томского государственного университета. Математика и механика. 2024. № 91. С. 5-17. аог 10.17223/19988621/91/1

© Л.С. Брындин, 2024

Original article

Development and verification of the least-squares collocation method with seventh degree polynomials for the biharmonic equation

Luka S. Bryndin

Khristianovich Institute of Theoretical and Applied Mechanics of the Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russian Federation, [email protected]

Abstract. A new version of the least-squares collocation method (C3-LSCM7) with seventh degree polynomials is proposed and implemented. The method has continuity up to the third derivatives of the piecewise polynomial solution in the sense of least squares. This is achieved by using the values of the solution and its derivatives at the vertices of the grid cells as unknown terms. The C3-LSCM7 is fundamentally different from previous versions of the LSCM in the absence of matching conditions. They explicitly require continuity of the solution and its derivatives at several points on the boundaries between neighboring cells. To solve the problem in a domain with a curved boundary a grid is constructed with rectangular cells. The solution of small irregular cells is continued from neighboring independent ones. Verification of the C3-LSCM7 is carried out by solving two-dimensional boundary value problems for a biharmonic equation in a square and in domain with the curvilinear boundary. The condition numbers of a global matrix and transition matrices from values at nodes to coefficients of polynomial expansion are studied. The advantages of the C3-LSCM7 over previous versions of the LSCM are shown. Keywords: least-squares collocation method, piecewise polynomials, automatic solution continuity, biharmonic equation, plate bending

Acknowledgments: The research was carried out within the state assignment of Ministry of Science and Higher Education of the Russian Federation.

For citation: Bryndin, L.S. (2024) Development and verification of the least-squares collocation method with seventh degree polynomials for the biharmonic equation. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika - Tomsk State University Journal of Mathematics and Mechanics. 91. pp. 5-17. doi: 10.17223/19988621/91/1

Введение

Краевые задачи для бигармонического уравнения четвертого порядка встречаются в различных разделах механики сплошных сред, например в теории Кирхгофа-Лява изгиба пластин [1] и гидродинамике при малых числах Рейнольд-са [2]. Главной особенностью и одновременно сложностью решения таких задач является наличие высоких производных, входящих как в само уравнение (четвертый порядок), так и в граничные условия (вплоть до третьего порядка) [1], что приводит к плохой обусловленности [2-4].

Известно несколько основных способов численного решения бигармонического уравнения. Первый из них связан с непосредственной дискретизацией исходной задачи в сильной постановке. Так поступают, например, в коллокационных (псевдо) спектральных методах (КМ) [3-5] и методах конечной разности (МКР) [6].

Другой подход основан на сведении бигармонического уравнения к системе, состоящей из двух уравнений с частными производными (УЧП), в которых порядок старших производных равняется двум [1, 2, 7]. Однако в этой ситуации может возникнуть дополнительная работа с краевыми условиями [2]. Наконец, в методе конечных элементов (МКЭ) ставится слабая постановка задачи, что приводит к понижению порядка старшей производной исходного решения в обновленном УЧП с возможностью использовать элементы второго или даже первого порядка [8].

Отметим, что, несмотря на хорошо проработанную технику решения линейных УЧП четвертого порядка МКР [2, 6] и МКЭ [8] существует не так много опубликованных результатов интегрирования бигармонического уравнения в нерегулярных областях и со сложными граничными условиями с помощью этих методов [3]. Кроме того, в КМ по сравнению с МКЭ уменьшаются временные затраты на формирование системы линейных алгебраических уравнений (СЛАУ) за счет отсутствия процедуры интегрирования, а также в некоторых случаях количество ненулевых элементов матрицы СЛАУ у КМ меньше [9]. Таким образом, КМ обладают своими достоинствами, выраженными прежде всего простотой применения и вычислительными аспектами. Однако решение УЧП высокого порядка в сильной постановке с помощью КМ может привести к матрице СЛАУ с достаточно большим числом обусловленности [3, 4], особенно в случае использования полиномов высоких степеней. Поэтому в таких ситуациях целесообразно использовать и исследовать влияние различных предобусловливателей [10], а также решать СЛАУ численными методами линейной алгебры, которые не ухудшают число обусловленности исходной матрицы, например ортогональными.

Метод коллокации и наименьших квадратов (МКНК) [4, 5, 11, 12] является одним из представителей КМ. Его принципиальное отличие от многих других заключается в необходимости решения переопределенных СЛАУ, состоящих из уравнений коллокации, условий согласования и краевых условий. При ее решении большую роль играют значения весовых множителей, на которые домножа-ются уравнения приближенной задачи. Они влияют на точность решения, скорость сходимости метода итераций по подобластям решения глобальной СЛАУ и ее обусловленность [11, 12]. В настоящий момент в МКНК поиск оптимальных значений весовых коэффициентов проводится экспериментальным путем. Для каждой задачи и вида используемых полиномиальных пространств это своя отдельная и непростая история. Про подобные проблемы выбора весовых множителей также указано и в других работах, например в [13, 14].

Для решения бигармонического уравнения сеточным КМ требуется склеивать производные приближенного решения вплоть до третьего порядка на границах между соседними ячейками [4, 15]. При этом в МКНК возникает четыре весовых множителя [4], а в сумме с весовыми коэффициентами уравнений коллокации и краевых их количество равняется семи. В работе [11] показана сильная чувствительность числа итераций сходимости решения глобальной СЛАУ в зависимости от значений весовых множителей условий согласования в случае, когда в МКНК применялись полиномы четвертой степени для аппроксимации бигармонического уравнения. Очевидно, что построение варианта МКНК без условий согласования заметно снизит зависимость численного решения от значений весовых множителей, которых станет еще и меньше.

Цель данной работы - построение МКНК с автоматической склейкой решения и его производных вплоть до третьего порядка, который позволит решить двумерное бигармоническое уравнение в сильной постановке. На сегодняшний день в ряде случаев уже разработаны различные варианты КМ, основанные на сплайн-аппроксимации, для УЧП, содержащих старшие производные второго порядка. Так, в работе [16] построен КМ с использованием полиномов второй степени для решения уравнений Пуассона и Гельмгольца. В [7] КМ с полиномами второй степени применялся для решения бигармонического уравнения, которое было сведено к системе, состоящей из двух УЧП второго порядка. В настоящей статье разработан новый вариант МКНК решения краевых задач для бигармонического уравнения без перехода к системе, принципиально отличающийся от предыдущих (см.: [4, 5, 11] и цитируемые там работы):

- отсутствием условий согласования;

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

Отметим, что в предыдущих работах по применению МКНК для решения бигармонического уравнения (см.: [4, 5, 11] и цитируемые там работы) неизвестными при решении глобальной переопределенной СЛАУ выступали как раз коэффициенты полинома, а в данной работе - отмеченные выше значения решения и его производных в вершинах ячеек сетки. Численно установлено: чтобы получить прямоугольную матрицу глобальной СЛАУ и прямоугольную матрицу перехода с конечным числом обусловленности и склейкой решения вплоть до третьих производных на сетке с прямоугольными ячейками, необходимо использовать пространства полиномов седьмой степени. Поэтому такой вариант МКНК назовем ^-МКНК Для решения задачи в области с криволинейной границей здесь применялся сеткопостроитель из статьи [5].

Постановка задачи и описание численного метода

Рассмотрим бигармоническое уравнение на примере его приложения к теории Кирхгова-Лява изгиба пластин [1] в двумерной области имеющее следующий вид:

д w „ д w д w q „

-Г + 2—5—7 +-Г= , (*,У) еО, (1)

Sx4 сх2ду2 dy4D ( )

где w(x, y) - прогиб изотропной пластины при изгибе поперечной нагрузкой

Et3

q(x, y), D =-— - жесткость пластины при изгибе, t = const - толщина,

12(1 -v2)

E = const и v = const - модуль Юнга и коэффициент Пуассона материала соответственно.

Дополним (1) краевыми условиями защемления [1]

w = 0, — = 0, (х, у) еда (2)

дп

или шарнирного закрепления

w = 0, Ми = Mxtf + Myn\ + 2Mxynln1 = 0, (х,у) е дО, (3)

где n = (n, n2) - внешняя

tfw ду

M = - D

(д2 w

-r + v—

dx2 э-2

единичная

(

My = - D

нормаль к границе области dQ,

д2 w

д w д w .

v^r + ^r\, Mxy = -D(i + v)

ду )' х 4 ' дхду

Задача данного исследования заключается в разработке нового сеточного варианта МКНК без условий согласования, позволяющего численно решить краевые задачи (2) и (3) для (1) в квадратной области и области с криволинейной границей областях.

Для построения расчетной сетки в случае областей с криволинейной границей используется подход из [5]. Кратко приведем основные его моменты. Поместим область ^ в прямоугольник и покроем его сеткой, состоящей из N х N прямоугольных ячеек. Граничными нерегулярными ячейками будем называть ячейки, имеющие непустое пересечение с дО и не совпадающие с исходными материнскими прямоугольными ячейками, которые содержат их внутри себя. Если центр материнской ячейки находится внутри то такие ячейки считаются самостоятельными, иначе - несамостоятельными (см. ячейки с номерами 1-5 на рис. 1). В качестве решения в каждой фиксированной точке несамостоятельных ячеек принимается решение, продолженное из самостоятельной ячейки, центр которой находится ближе всего к рассматриваемой точке. Так, например, на рис. 1 решение из ячейки с номером 6 продолжается в светло-серую область несамостоятельной ячейки 4, из ячейки номер 7 - в темно-серую ее область.

3° я О OJB> о о о о о о О О О О •е-9 о о о о о о о о 0 0 о о е^йо о о о^ао о о о о о о о о 4 та °

о <яГО о о Ъ о о О О 0 0 О О О О 0 о о о 0 0 О 0 О О О 0 О 0 0 0 OjO ЦО 0 0 о_8. о 7 "L

<А о о о Ф о о о о о о о о о о о о о о о о о о о о о о о о о о о о о о ш> о о о 9

о о о i> о о о о о о о О О О О г \ 5° ч / о о о о 0 о о о о о о ^ о о о <\

Ь О О О Ьо о о О О О О ❖ ООО 0 О О 0 О 0 О 0 ооо 4 0 о о i

& О 0 0 <Я О О О 0 0 0 0 О О 0 0 0 ❖ О 0 0 0 О 0 0 0 0 0 О 0 0 0 0 0 0 ф 0 о о и>

° ы ° ° о ofa о о о опо о о о о о о о о о о о о о о о о о о о о о О О gf о о оjao о

г о о о о о о о о OsB> О 0 о о^з-в. о о о о о о о о 0 0 О 0 о о о о о о о о О 0 -¿ЯГО о у>дП °2

Рис. 1. Пример расчетной сетки в кольце в МКНК при N х N = 4 х 4, где символ «О» означает точки коллокации, «□» - точки записи краевых условий, «о» - центры материнских ячеек, включающих нерегулярные несамостоятельные ячейки Fig. 1. Example of a calculation grid in a ring in the LSCM at jV x A? = 4 x4 where the symbol □ means collocation points, □ are points of boundary condition records, and о are centers of mother cells including irregular nonautonomous cells

В каждой j -й самостоятельной ячейке, j = 1,...,Ncdb, где NceIb - их количество, введем локальную систему координат

=

x - x„

hi

^2 =

y - yc h

(4)

где (xc, yc) - центрj-й прямоугольной ячейки, 2h, 2й2 - ее размеры.

В каждой_/-й ячейке приближенное решение представляется в виде:

7 7-4

^(Ьж,= ЦЬ = (с),V], (5)

2=о 2=о

где С^,] - 36 неизвестных коэффиЦиентов, с] = (Соо ],С0Ж]С^,],...,С70,])Т ,

V ] = (ж, £„..., ЬЖ2 Ь2,..., Ь7)Т.

Будем искать приближенное решение в виде совокупности его значений и значений всех его производных вплоть до третьего порядка в каждом узле (вершине) самостоятельных прямоугольных ячеек. Прямоугольная матрица перехода А размера 40 х 36 позволяет перейти от коэффициентов разложения с.г в (5)

к необходимым значениям в узлах посредством решения следующей переопределенной СЛАУ:

Ac. = w. = (wh1,w*-1,wh1,...,wh1 ,wh'2,...,wh'2 .,...,wh'4,...,wM .)T,

t J J v J ' X,J' y,J' ' УУУ,Р ])■■■■> yyy,J' ' J ' ' УУУ,]' '

(6)

где цифра в верхнем индексе после h обозначает значение в одном из четырех узлов, буквенный нижний индекс - частную производную, например

я3 ь

ьж д ^ /гЖ *-Ж\ /5-Ж 5-Ж Л

/ =—5— (ьп Чг), (Ь ,ь2) - локальные координаты первого из четырех узлов

& ду

ячейки. Вектор с. из (6) определяется в смысле наименьших квадратов и равен с . = , где ОЯ - QR-разложение матрицы А : А = ОЯ, О - матрица

с ортогональными столбцами размера 40 х 36, Я - верхнетреугольная матрица размера 36 х 36.

Отметим, что наличие локальной системы координат (4) влечет за собой одинаковые координаты узлов каждой ячейки и одинаковый вид матриц перехода А . При попытке использовать для аппроксимации (5) полиномы более низкой степени оказывается, что итоговая СЛАУ имеет неполный ранг, а в случае применения полиномов более высокой степени матрица А становится недоопределен-ной. В обоих этих вариантах решение адекватно не удавалось найти.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Для определения ^ в узлах сетки в каждой ячейке выписываются уравнения коллокации и краевые условия, если ячейка является граничной. Пусть точка коллокации в ]-й самостоятельной ячейке имеет координаты (ЬЖс, Ь2с). Тогда подстановка приближенного решения (5) в уравнение (1), умноженное на весовой множитель рс, в этой точке дает следующее выражение:

+ д4 1 _д_4 Л х д^ + ** д^2 + h4у

7 7-1.

Ц^, ] ьУс = . (7)

¿,=012=0 ^

Обозначим вектор размера 36, элементами которого являются Ь^, через V, а векторы, составленные из четвертых производных в точке (Ьк, ) по Ь

д4 д4 д4

четырежды, по Ь и по Ь2 дважды, по Ь2 четырежды, через —-V, —-—-V и —-V соответственно. Тогда уравнение (7) с учетом с = можно переписать в виде:

Г 1 т 1 д4 1 д4 1 д4 ) с

Рс К Q w } >77^77 у + 2Ъ2,2 V + 74774 у 1 = Р 77 . (8)

Преобразуем левую часть уравнения (8), перенося операцию умножения на матрицу во второй аргумент скалярного произведения, в результате чего

получим следующее уравнение на неизвестные значения решения и его производных в узлах у-й ячейки:

( ( 1 а4 1 д4 1 д4 ЛЛ

w ,, Q(R-1)т

=p. (9)

v v h v + h hi v+h 8^4v

В каждой самостоятельной ячейке выписываются уравнения коллокации (9) в точках, которые являются прямым произведением множеств, состоящих из корней полиномов Лежандра четвертой степени (см. рис. 1).

Аналогичным образом выписываются уравнения в граничных ячейках. Так, например, если точка границы области 8Q имеет локальные координаты , ),

вектор v состоит из элементов ^ j6 ^ , то краевое условие для первого уравнения (2) будет иметь вид:

pK (w,, Q(R-1)T v) = 0, (10)

где p^ - весовой множитель первого уравнения (2) или (3).

Для записи краевых условий в граничной ячейке по алгоритму из [5] выбирается часть границы области, на которой в зависимости от количества соседних самостоятельных ячеек Nnei^ равномерно расставляется 4(4 - Nnel^) точек.

Объединяя уравнения коллокации и краевые условия в каждой ячейке, получим глобальную переопределенную СЛАУ, для решения которой используется диагональный предобусловливатель [10] и прямой ортогональный метод, реализованный в библиотеке SuiteSparse [17] с возможностью распараллеливания на графических процессорах с помощью технологии CUDA.

Результаты численных экспериментов

В представленных ниже результатах численных экспериментов погрешность вычислялась по следующим формулам:

max max | wh (x,, y) - (x,, y )|

" r max max | we (x., y. )| ,

где wh, wex - приближенное и точное решение задачи соответственно, (xt, yt) -равномерно распределенные в у-й ячейке точки. В точках, не принадлежащих области погрешность не вычислялась.

Порядок сходимости определялся по формуле

И»

logi E

R = 1

logi

где Ш

- погрешности (11) на сетках с общим количеством неизвестных

N и N2 соответственно.

Пример 1

Рассмотрим квадратную шарнирно закрепленную (3) пластину О = (0,1)2, на которую действует нагрузка д(х, у) = 0.18т(ях)8т(я:у). В таком случае прогиб пластины определяется следующим выражением [1]:

. q(x, у)

4 *=■ 4я4 D

(13)

На рис. 2 приведены результаты сходимости С3-МКНК7 и hp-варианта (hp-МКНК) из работы [5], в котором также использовались полиномы седьмой степени (5) для аппроксимации бигармонического уравнения (1), где n = 4dof , DOF- общее количество степеней свободы, т.е. неизвестных, используемых в вариантах МКНК для определения решения, C, C2 = const. В работе [4] были продемонстрированы достоинства h-, p- и hp-МКНК по сравнению с МКР [2], МКЭ [8] и спектральными методами [3, 15, 19] на примерах решения различных краевых задач для бигармонического уравнения в нерегулярных областях. Параметры С3-МКНК7 были взяты следующими: pc = h2h2 , весовые множители перед

h h

первым и вторым уравнениями (3) равнялись p^ = 1, p^ = . В hp-МКНК [5] h h

pc = hj2h2 , pbo = 1, p42 =~jy, весовые коэффициенты условий согласования перед решением, первой, второй и третьей производными по нормали между соседними ячейками полагались равными рщ = 1, рщ = 0.12, рщ = h1h2, рщ = 0.125h1h2

соответственно. Глобальная СЛАУ в hp-МКНК решалась с помощью метода итераций по подобластям [11].

0.0 J-2.5

3 s-

кГ

м -5.0 -7.5

N ь N 4 \

4 \ \ > \ N \ N \

Ci 4 A \ ¿J \ S 4 N

1.5 2.0

п

b

Рис. 2. Сходимость С3-МКНК7 (А) и hp-МКНК [5] (о) в примере 1 (а) и С3-МКНК7 (А) в примере 2 (b) Fig. 2. Convergence of the (a) C3-LSCM7 (А) and hp-LSCM [5] (о) in example 1 and (b) C3-LSCM7 (А) in example 2

а

Из рисунка видно, что в С3-МКНК7 достигаются лучшие значения погрешности по сравнению с йр-МКНК, а также порядок сходимости С3-МКНК здесь

выше и равен приблизительно 8. Отношения количества уравнений к количеству неизвестных в СЛАУ в двух методах отличаются совсем немного: степень переопределения СЛАУ в С3-МКНК равнялась приблизительно 1.6, а в йр-МКНК - 1.8.

Исследуем число обусловленности матриц A глобальных СЛАУ в С3-МКНК7 и влияние диагонального предобусловливателя [10]. В этом примере на сетке 10 х 10 (БОЕ = 1 210) число обусловленности в спектральной норме еопй2 (А) = 9.97е + 8,

а с применением предобусловливания сопй2 (А) = 9.54е + 8. Аналогичная картина с точки зрения отсутствия разницы между значениями сопй2 (А) также наблюдалась в примере 2. Числа обусловленности матриц перехода сопй2 (А) варьировали от 104 до 106.

Для оценки погрешности склейки проведены расчеты значений приближенного решения и его частных производных вплоть до третьего порядка во всех узлах области. В каждом узле указанные значения рассчитывались во всех ячейках, содержащих рассматриваемый узел, и находилась максимальная относительная разница между ними. Установлено, что все отклонения находятся на уровне погрешности приближенного решения. Например, на сетке размера 10 х 10 погрешность Ц-Е^Ц = 3.07e-7, а отклонения рассматриваемых величин находятся в пределах от 1.42е-10 до 6.97е-7.

Пример 2

Рассмотрим круглую защемленную по краям пластину с центральным отверстием О = {(х, у): 1 < х2 + у2 < 52}, на которую действует постоянная нагрузка д. В этом случае решение задачи имеет вид [18]:

4 С 2

м,(х, у) = + —г2 (1оё2 г -1) + —4— + — ^ г + —

(14)

64В 3 2 4

где г х2 + у2 . Константы С3-С6 определяются из решения СЛАУ размера 4 х 4,

которая получается подстановкой (14) в краевые условия на внешнем и внутреннем контурах пластины. В силу громоздкости выражения для них здесь не приводятся.

В таблице и на рис. 2, Ь приведены результаты сходимости С3-МКНК7. В таблице число означает количество самостоятельных ячеек. В данном случае порядок сходимости метода колебался от 4 до 7.

Результаты сходимости C3-МКНК7 в примере 2

N х N БОЕ — да Я

5 х 5 20 320 1.30е+0 -

10 х 10 76 1 000 1.30е-1 4.04

20 х 20 304 3 520 4.11е-3 5.49

40 х 40 1 212 13 080 2.21е-4 4.45

80 х 80 4 816 50 080 1.46е-6 7.48

Несмотря на то, что второй пример может быть решен в одномерной постановке, здесь было важно, во-первых, показать возможность удобного применения сеткопостроителя [5] для областей с криволинейной границей при использовании прямоугольных ячеек в С3-МКНК7, а во-вторых, провести верификацию метода на задаче с известным точным решением. Более того, такую задачу не удастся решить р-вариантом различных КМ [3, 4, 19], в которых используются глобальные аппроксимации во всем прямоугольнике, включающем исходную область, без перехода, например, в полярную систему координат, поскольку решение (14) содержит слагаемые logr, для которых (0, 0) является точкой ветвления, где решения обращаются в бесконечность. Отметим также, что в более сложных областях перейти в какую-то специальную удобную систему координат далеко не всегда возможно. Таким образом, это подчеркивает актуальность развития сеточных численных методов.

Заключение

Впервые разработан и исследован вариант МКНК без условий согласования, позволяющий решать краевые задачи для бигармонического уравнения в сильной постановке. Его верификация проведена на решении задач изгиба пластин в рамках теории Кирхгофа-Лява. Сравнение реализованного здесь С3-МКНК7 с hp-МКНК из работы [5] при одинаковых степенях полиномов и DOF показало преимущества нового варианта. На тестовом бесконечно гладком решении в квадратной области показан высокий порядок сходимости, в среднем равный восьми. Продемонстрирована возможность применения С3-МКНК7 для решения бигармониче-ского уравнения в области с криволинейной границей. Однако вопрос уменьшения числа обусловленности cond2 (A) в МКНК с применением полиномов высоких степеней остается важным и открытым.

Список источников

1. Тимошенко С.П., Войновский-Кригер С. Пластины и оболочки. М. : Физматгиз, 1996. 625 с.

2. Chen G., Li Z., Lin P. A fast finite difference method for biharmonic equations on irregular

domains and its application to an incompressible Stokes flow // Advances in Computational Mathematics. 2007. V. 29, № 2. P. 113-133. doi: 10.1007/s10444-007-9043-6.

3. Shao W., Wu X., Chen S. Chebyshev tau meshless method based on the integration-differen-

tiation for biharmonic-type equations on irregular domain // Engineering Analysis with Boundary Elements. 2012. V. 36, № 12. P. 1787-1798. doi: 10.1016/j.enganabound.2012.06.005

4. Беляев ВА., Брындин Л.С., Голушко С.К., Семисалов Б.В., Шапеев В.П. H-, р- и hp-варианты

метода коллокации и наименьших квадратов для решения краевых задач для бигармо-нического уравнения в нерегулярных областях и их приложения // Журнал вычислительной математики и математической физики. 2022. Т. 62, № 4. С. 531-552. doi: 10.31857/S0044466922040020

5. Брындин Л.С., Беляев В.А., Шапеев В.П. Разработка и верификация упрощенного hp-

варианта метода коллокации и наименьших квадратов для нерегулярных областей // Вестник Южно-Уральского государственного университета. Сер. Математическое моделирование и программирование. 2023. Т. 16, № 3. С. 35-50. doi: 10.14529/mmp230303

6. Ben-Artzi M., Croisille J.P., Fishelov D. An Embedded Compact Scheme for Biharmonic

Problems in Irregular Domains // Advanced Computing in Industrial Mathematics: 11th Annual Meeting of the Bulgarian Section of SIAM. Cham : Springer, 2018. V. 728. P. 11-23. doi: 10.1007/978-3-319-65530-7 2

7. Bialecki B., Fairweather G., Karageorghis A., Maack J. A quadratic spline collocation method

for the Dirichlet biharmonic problem // Numerical Algorithms. 2019. V. 83. P. 165-199. doi: 10.1007/s11075-019-00676-z

8. Guo H., Zhang Z., Zou Q. A C0 Linear Finite Element Method for Biharmonic Problems //

Journal of Scientific Computing. 2018. V. 74, № 3. P. 1397-1422. doi: 10.1007/s10915-017-0501-0

9. Schillinger D., Evans J.A., Reali A., Scott M.A., Hughes T.J.R. Isogeometric collocation: Cost

comparison with Galerkin methods and extension to adaptive hierarchical NURBS discretizations // Computer Methods in Applied Mechanics and Engineering. 2013. V. 267. P. 170232. doi: 10.1016/j.cma.2013.07.017

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

10. Ramsak M., Skerget L. A subdomain boundary element method for high-Reynolds laminar flow using stream function-vorticity formulation // International Journal for Numerical Methods in Fluids. 2004. V. 46, № 8. P. 815-847. doi: 10.1002/fld.776

11. Голушко С.К., Идимешев С.В., Шапеев В.П. Метод коллокаций и наименьших невязок в приложении к задачам механики изотропных пластин // Вычислительные технологии. 2013. Т. 18, № 6. С. 31-43.

12. Исаев В.И., Шапеев В.П., Еремин С.А. Исследование свойств метода коллокации и наименьших квадратов решения краевых задач для уравнения Пуассона и уравнений Навье-Стокса // Вычислительные технологии. 2007. Т. 12, № 3. С. 53-70.

13. Ильин В.П., Кныш Д.В. Параллельные методы декомпозиции в пространствах следов // Вычислительные методы и программирование. 2011. Т. 12, № 1. С. 110-119.

14. Алексидзе М.А. Решение граничных задач методом разложения по неортогональным функциям. М. : Наука, 1978. 352 с.

15. Shao W., Wu X. An effective Chebyshev tau meshless domain decomposition method based on the integration-differentiation for solving fourth order equations // Applied Mathematical Modelling. 2015. V. 39, № 9. P. 2554-2569. doi: 10.1016/j.apm.2014.10.048

16. Fairweather G., Karageorghis A., Maack J. Compact optimal quadratic spline collocation methods for Poisson and Helmholtz problems: formulation and numerical verification // Journal of Computational Physics. 2010. V. 230, № 8. P. 2880-2895. doi: 10.1016/j.jcp.2010.12.041

17. Davis T.A. Algorithm 915, SuiteSparseQR: Multifrontal multithreaded rank-revealing sparse QR factorization // ACM Transactions on Mathematical Software. 2011. V. 38, № 1. P. 8:18:22. doi: 10.1145/2049662.2049670

18. Ike C.C. Mathematical solutions for the flexural analysis of Mindlin's first order shear de-formable circular plate // Mathematical Models in Engineering. 2018. V. 4, № 2. P. 50-72. doi: 10.21595/mme.2018.19825

19. Mai-Duy N., See H., Tran-Cong T. A spectral collocation technique based on integrated Che-byshev polynomials for biharmonic problems in irregular domains // Applied Mathematical Modelling. 2009. V. 33, № 1. P. 284-299. doi: 10.1016/j.apm.2007.11.002

References

1. Timoshenko S., Woinowsky-Krieger S. (1959) Theory of Plates and Shells. New York,

St. Louis: McGraw-Hill Book Comp.

2. Chen G., Li Z., Lin P. (2007) A fast finite difference method for biharmonic equations on

irregular domains and its application to an incompressible Stokes flow. Advances in Compu-tationalMathematics. 29(2). pp. 113-133. DOI: 10.1007/s10444-007-9043-6.

3. Shao W., Wu X., Chen S. (2012) Chebyshev tau meshless method based on the integration-

differentiation for biharmonic-type equations on irregular domain. Engineering Analysis with BoundaryElements. 36(12). pp. 1787-1798. DOI: 10.1016/j.enganabound.2012.06.005.

4. Belyaev V.A., Bryndin L.S., Golushko S.K., Semisalov B.V., ShapeevV.P. (2022) H-, p-, and

hp-versions of the least-squares collocation method for solving boundary value problems for biharmonic equation in irregular domains and their applications. Computational Mathematics and Mathematical Physics. 62(4). pp. 517-537. DOI: 10.1134/S0965542522040029.

5. Bryndin L.S., Belyaev V.A., Shapeev V.P. (2023) Razrabotka i verifikatsiya uproshchennogo

hp-varianta metoda kollokatsii i naimen'shikh kvadratov dlya neregulyarnykh oblastey [Development and verification of a simplified hp-version of the least-squares collocation method for irregular domains]. Vestnik Yuzhno-Ural'skogo gosudarstvennogo universiteta. Seriya Matematicheskoe modelirovanie i programmirovanie - Bulletin of the South Ural State University. Series Mathematical Modelling, Programming & Computer Software. 16(3). pp. 3550. DOI: 10.14529/mmp230303.

6. Ben-Artzi M., Croisille J.P., Fishelov D. (2018) An Embedded Compact Scheme for Bi-

harmonic Problems in Irregular Domains. Advanced Computing in Industrial Mathematics: 11th Annual Meeting of the Bulgarian Section of SIAM. Cham: Springer. 728. pp. 11-23. DOI: 10.1007/978-3-319-65530-7_2.

7. Bialecki B., Fairweather G., Karageorghis A., Maack J. (2019) A quadratic spline collocation

method for the Dirichlet biharmonic problem. Numerical Algorithms. 83. pp. 165-199. DOI: 10.1007/s11075-019-00676-z.

8. Guo H., Zhang Z., Zou Q. (2018) A C0 Linear Finite Element Method for Biharmonic Problems.

Journal of Scientific Computing. 74(3). pp. 1397-1422. DOI: 10.1007/s10915-017-0501-0.

9. Schillinger D., Evans J.A., Reali A., Scott M.A., Hughes T.J.R. (2013) Isogeometric colloca-

tion: Cost comparison with Galerkin methods and extension to adaptive hierarchical NURBS discretizations. Computer Methods in Applied Mechanics and Engineering. 267. pp. 170-232. DOI: 10.1016/j.cma.2013.07.017.

10. Ramsak M., Skerget L. (2004) A subdomain boundary element method for high-Reynolds laminar flow using stream function-vorticity formulation. International Journal for Numerical Methods in Fluids. 46(8). pp. 815-847. DOI: 10.1002/fld.776.

11. Golushko S.K., Idimeshev S.V., Shapeev V.P. (2013) Metod kollokatsiy i naimen'shikh nevyazok v prilozhenii k zadacham mekhaniki izotropnykh plastin [Application of collocations and least residuals method to problems of the isotropic plates theory]. Vychislitel'nye tekhnologii - Computational Technologies. 18(6). pp. 31-43.

12. Isaev V.I., Shapeev V.P., Eremin S.A. (2007) Issledovaniye svoystv metoda kollokatsii i naimen'shikh kvadratov resheniya krayevykh zadach dlya uravneniya Puassona i uravneniy Nav'ye-Stoksa [An investigation of properties of the least squares collocation method for solution of boundary value problems for the Navier-Stokes and Poisson equations]. Vychislitel'nye tekhnologii - Computational Technologies. 12(3). pp. 53-70.

13. Il'in V. P., Knysh D. V. (2011) Parallel'nyye metody dekompozitsii v prostranstvakh sledov [Parallel decomposition methods in the trace spaces]. Vychislitel'nye metody iprogrammirovanie - Numerical Methods and Programming. 12(1). pp. 110-119.

14. Aleksidze M. A. (1978) Resheniye granichnykh zadach metodom razlozheniya po neortogo-nal'nym funktsiyam [Solution of boundary value problems by the method of decomposition by non-orthogonal functions]. Moscow: Nauka.

15. Shao W., Wu X. (2015) An effective Chebyshev tau meshless domain decomposition method based on the integration-differentiation for solving fourth order equations. Applied Mathe-maticalModelling. 39(9). pp. 2554-2569. DOI: 10.1016/j.apm.2014.10.048.

16. Fairweather G., Karageorghis A., Maack J. (2010) Compact optimal quadratic spline collocation methods for Poisson and Helmholtz problems: formulation and numerical verification. Journal of Computational Physics. 230(8). pp. 2880-2895. DOI: 10.1016/j.jcp.2010.12.041.

17. Davis T.A. (2011) Algorithm 915, SuiteSparseQR: Multifrontal multithreaded rank-revealing sparse QR factorization. ACM Transactions on Mathematical Software. 38(1). pp. 8:1-8:22. DOI: 10.1145/2049662.2049670.

18. Ike C.C. (2018) Mathematical solutions for the flexural analysis of Mindlin's first order shear deformable circular plate. Mathematical Models in Engineering. 4(2). pp. 50-72. DOI: 10.21595/mme.2018.19825.

19. Mai-Duy N., See H., Tran-Cong T. (2009) A spectral collocation technique based on integrated Chebyshev polynomials for biharmonic problems in irregular domains. Applied Mathematical Modelling. 33(1). pp. 284-299. DOI: 10.1016/j.apm.2007.11.002.

Сведения об авторе:

Брындин Лука Сергеевич - младший научный сотрудник Института теоретической и прикладной механики им. С.А. Христиановича Сибирского отделения РАН (Новосибирск, Россия). E-mail: [email protected]

Information about the author:

Bryndin Luka S. (Junior Researcher, Khristianovich Institute of Theoretical and Applied Mechanics of the Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russian Federation). E-mail: [email protected]

The article was submitted 29.02.2024; accepted for publication 03.10.2024

Статья поступила в редакцию 29.02.2024; принята к публикации 03.10.2024

i Надоели баннеры? Вы всегда можете отключить рекламу.