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

УДК 519.651

ОБ ОРТОГОНАЛЬНЫХ ФИНИТНЫХ ФУНКЦИЯХ ВТОРОЙ СТЕПЕНИ, СВЯЗАННЫХ С ТРЕУГОЛЬНЫМИ СЕТКАМИ,

И ИХ ПРИМЕНЕНИИ В МАТЕМАТИЧЕСКОМ МОДЕЛИРОВАНИИ

В. Л. Леонтьев, А. В. Кочулимов

Ульяновский государственный университет,

432970, Ульяновск, ул. Л. Толстого, 42.

E-mails: leontievvl@ulsu.ru

Рассматриваются ортогональные финитные функции второй степени, компактные носители которых представляют собой локальные совокупности треугольников сетки. Формулируется теорема об аппроксимирующих свойствах последовательностей наборов таких функций. Выделяются области применения функций: в алгоритмах смешанных численных методов решения краевых задач, в алгоритмах линейной аппроксимации поверхностей.

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

Ортогональные финитные функции второй степени строятся на сетке треугольников, созданной в некоторой области 0, на основе классических финитных функций [1, с. 125-130] следующим образом. Каждая треугольная часть Р\, Рз, Р5 конечного носителя функции [1] разбивается на четыре треугольные подобласти тремя линиями, соединяющими середины Р2, Р4, Рб сторон Р1-Р3, Р3-Р5, Р5-Р1, затем каждая из четырёх полученных подобластей делится на три треугольника линиями, соединяющими вершины этих подобластей и одну из точек М1, М2, М3, М4 пересечения медиан треугольных подобластей. На всех двенадцати треугольных подобластях строятся вспомогательные финитные функции второй степени, равные единице в узле Мг и нулю во всех остальных пяти узлах малого треугольника, имеющего вершину в точке Мг, их конечный носитель — треугольники, окружающие узел Мг. Вспомогательные финитные функции, взятые с постоянными коэффициентами, складываются с одной из шести классических финитных функций второй степени [1], связанных с узлами Р1, Р2, Рз, Р4, Р5, Рб, таким образом модифицируя её. Коэффициенты модифицирующих линейных комбинаций выбираются, исходя из условий взаимной ортогональности шести функций, связанных с узлами Р1, Р2, Рз, Р4, Р5, Рб основной сетки. Графики функций Р(х,у), полученных для ведущих узлов Р1, Р4, изображены на рисунке для одной треугольной части конечных носителей финитных функций, связанных с ведущими узлами Р1, Р4. Аналогичный вид имеют модифицированные ортогональные финитные функции, связанные с узлами Р2, Рз, Р5, Рб.

Пусть — подпространство построенных сеточных функций ^>г(х,у), где С С(0) П ^21(0), С(0) — пространство непрерывных на 0 функций с нормой ||«||с= тах|и(х)|, С(г)(0) —банахово пространство функций,

Виктор Леонтьевич Леонтьев (д.ф.-м.н., профессор), профессор, каф. информационной безопасности и теории управления. Александр Валерьевич Кочулимов, аспирант, каф. информационной безопасности и теории управления.

Ш)

Модифицированные функции, связанные с узлами Рі, Р4

имеющих непрерывные в О производные до порядка Г включительно, (О) (г = 0,1, 2,...) — гильбертовы пространства с нормами

МІ»Ї (П)

Е

V <г

Г \1/2

I \в"и\2,

где ^ = дн /(дтхдпу) — оператор частного дифференцирования т раз по х и п раз по у; V = (т, п); т и п принимают значения 0, 1, 2; \и\ = т + п.

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

Теорема. Для любой функции и(х,у) Є (О) П С(1)(О) существует

пн Є Ш1Л

и — и

І№» (П)

А1/251/2Н ^ |и|Ц2(П) + сіп й„ 2^

|Г| = 1

0

1^Г и\\с (П),

где 0о — минимальный из углов всех треугольников сетки; 5 — площадь области 0; А = (76 — 20^7 — 4^105 + 5^15)/480 ~ 0,0031; пн — линейные комбинации сеточных функций <^г(х, у); Н —максимальная длина сторон треугольников сетки; постоянный коэффициент С, не зависящий от параметров сетки и от функции п.

В первую область применения последовательностей наборов сеточных ортогональных финитных функций ^г(х,у), связанных с треугольными сетками, входят алгоритмы смешанных численных методов решения краевых задач теории математического моделирования. Такие алгоритмы основаны на смешанных вариационных принципах, в которых одновременно и независимо друг от друга варьируются как основные неизвестные функции (в механике деформируемого твердого тела — компоненты тензоров деформаций и напряжений, компоненты вектора перемещения) и их частные производные. Эта

особенность методов, как известно, приводит к получению приближённых решений одинаковой гладкости и одного порядка точности для всех кинематических и силовых факторов. Основной недостаток смешанных численных методов — повышенное число узловых неизвестных — устраняется благодаря использованию ортогональных финитных функций, которые в отличие от классических ортогональных многочленов и специальных функций позволяют решать двумерные задачи в областях с криволинейными границами. Вычислительные затраты этих методов такие же, как у численных методов «в перемещениях», основанных на вариационном принципе Лагранжа. Вторая область применения последовательностей наборов сеточных функций <^г(ж,у) —линейная аппроксимация поверхностей объектов и сглаживание локальных дефектов поверхностей. За счёт ортогональности финитных функций аппроксимация поверхностей эллиптического, параболического, гиперболического и смешанного типов проводится в несколько раз быстрее, чем в случае использования классических финитных функций. При этом относительная ошибка аппроксимации на сетках с 1-2 тысячами узлов составляет десятые и сотые доли процента, что удовлетворяет требованиям многих приложений.

Работа выполнена при поддержке ФЦП «Научные и научно-педагогические кадры инновационной России» (гос. контракты №№ П1122, П2230).

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Марчук Г. И., Агошков В. B. Введение в проекционно-сеточные методы. М.: Наука, 1981. 416 с. [Marchuk G. I., Agoshkov V. I. Introduction to projection-grid methods. Moscow: Nauka, 1981. 416 pp.]

Поступила в редакцию 01/X/2010; в окончательном варианте — 02/II/2011.

MSC: 42C40, 41A15

ON ORTHOGONAL SECOND ORDER FINITE FUNCTIONS

ASSOCIATED WITH TRIANGULAR MESHES

AND THEIR APPLICATION IN MATHEMATICAL MODELING

V. L. Leontiev, A. V. Kochulimov

Ulyanovsk State University,

42, L. Tolstoy st., Ulyanovsk, 4329170, Russia.

E-mails: leontievvl@ulsu.ru

Second order orthogonal finite functions are studied, their compact supports are local sets of grid triangles. The theorem on approximating the properties of sequences of sets of such functions is formulated. The applications in algorithms of mixed numerical methods of boundary value problems solution, in algorithms of linear approximation of surfaces are considered.

Key words: orthogonal finite functions, triangular grids, approximating properties, mathematical modeling, boundary value problems, surface approximation.

Original article submitted 01/X/2010; revision submitted 02/II/2011.

Victor L. Leontiev (Dr. Sci. (Phys. & Math.)), Professor, Dept. of Information Security & Management Theory. Alexander V. Kochulimov, Postgraduate Student, Dept. of Information Security & Management Theory.