CC BY
36
УДК 519.3
АН. ХОМЧЕНКО
Черноморский государственный университет им. Петра Могилы
Е.И. ЛИТВИНЕНКО, И.А. АСТИОНЕНКО
Херсонский национальный технический университет
БАРИЦЕНТРИЧЕСКИЕ КООРДИНАТЫ МЁБИУСА И НЕЛИНЕЙНЫЕ ФУНКЦИИ-"КРЫШКИ"
Посвящается 225-летию со дня рождения А.Ф. Мёбиуса
В статье рассматриваются возможные обобщения барицентрической задачи Мёбиуса на неоднородный стержень (одномерный симплекс). Как известно, классические результаты Мёбиуса связаны с однородным стержнем, когда априори работает архимедово правило рычага. Барицентрические координаты Мёбиуса ассоциируются с линейной функцией-"крышкой", которая часто применяется в методе конечных элементов в связи с кусочно-линейной интерполяцией функций. Модель неоднородного стержня можно построить с помощью нелинейных "полукрышек", образующих симметричную пару. Так поступают в задаче интерполяции сплайнами. Рассмотрены два характерных случая неравномерного распределения массы, когда плотность в средней точке стержня достигает минимума или максимума. Приведены примеры вычисления узловых нагрузок для стержня единичной массы со смещенным барицентром. Установлен вероятностный смысл построенных моделей.
Ключевые слова: барицентрические координаты, метод конечных элементов, кусочно-линейная интерполяция, узловые нагрузки, неоднородный стержень.
А.Н. ХОМЧЕНКО
Чорноморський державний ушверситет ím. Петра Могили
О.1. ЛИТВИНЕНКО, 1.О. АСТ1ОНЕНКО
Херсонський нащональний техшчний ушверситет
БАРИЦЕНТРИЧН1 КООРДИНАТИ МЕБ1УСА ТА НЕЛ1Н1ЙН1 ФУНКЦП-"КРИШКИ"
У статтi розглядаються можливi узагальнення барицентричног задачi Мебiуса на неоднорiдний стрижень (одномiрний симплекс). Як вiдомо, класичнi результати Мебiуса пов'язат з однорiдним стрижнем, коли апрiорi працюе архiмедове правило важеля. Барицентричт координати Мебiуса асоцiюються з лiнiйною функщею-"кришкою", яка часто застосовуеться в методi сктченних елементiв у зв 'язку з кусочно^тйною iнтерполяцiею функцт. Модель неоднорiдного стрижня можливо побудувати за допомогою нелiнiйних "напiвкришок", як утворюють симетричну пару. Так поводяться у задачi ттерполяци сплайнами. Розглянутi два характерних випадки нерiвномiрного розподшу маси, коли щшьшсть у середтй точцi стрижня досягае мiнiмума або максимума. Наведено приклади обчислення вузлових навантажень для стрижня одиничног маси зi змщеним барицентром. Установлений ймовiрнiсний зм^т побудованих моделей.
Ключовi слова: барицентричн координати, метод сктченних елементiв, кусочно^тйна iнтерполяцiя, вузловi навантаження, неоднорiдний стрижень.
A.N. KHOMCHENKO
Petra Mohyla Black Sea State University
YE.I. LITVINENKO, 1А. ASTIONENKO
Kherson National Technical University
MÖBIUS' BARYCENTRIC COORDINATES AND NON-LINEAR FUNCTIONS-"CAPS"
The probable generalizations of Möbius barycentric problem onto the inhomogeneous core (one-dimensional simplex) are shown in the article. As is known, Möbius' classical results are connected with homogeneous core, when Archimedes' Law of the Lever a priori works. Möbius' barycentric coordinates are associated with lihear function-"cap", which is often used in the finite elements method in connection with piecewise linear interpolation of functions. The model of inhomogeneous core can be built with the help of non-linear "half-caps", which form symmetric pair. This way is used in the problem of spline interpolation. Two characteristic cases of nonuniform mass distribution when density in the middle point of core reaches minimum or maximum are investigated. The examples of calculation of nodal loads for core of unit mass with shifted barycenter are shown. Probability meaning of built models is ascertained.
Keywords: barycentric coordinates, finite elements method, piecewise linear interpolation, nodal loads, inhomogeneous core.
Постановка проблемы
Решая барицентрическую задачу (о принудительном смещении барицентра симплекса), Мёбиус априори предполагал, что плотность симплекса постоянна. Так были открыты (1827 г.) удобные и полезные барицентрические координаты (Б-координаты), которые находят всё больше применений в геометрии, вычислительной математике и теории вероятностей. Представляет интерес исследование возможных обобщений на неоднородные симплексы классической задачи Мёбиуса с целью установить, какие свойства Б-координат наследуют новые решения.
Анализ предшествующих публикаций
Термин "крышка", вероятно, появился во второй половине ХХ столетия в работах инженерно ориентированных энтузиастов, развивающих метод конечных элементов (МКЭ, [1, 2]). В период становления и развития МКЭ появилось немало специфических терминов: ячейка Куранта, пирамида Куранта, клетка Кунса, функция-"пагода", симплекс, мультиплекс, серендипов элемент, "кирпич" Вильсона и т.п. На самом деле о линейной функции-"крышке" знал ещё Архимед, который в III в. до н.э. сформулировал "золотое правило" механики [3]. Так рычаг Архимеда стал прообразом одномерного симплекса, а линейные зависимости правила Архимеда - прообразом "полукрышек". Линейные функции-"крышки" были известны Ньютону, который в 1687 г. использовал разделенные разности для построения интерполяционной формулы. Позже Мёбиус распространил барицентрические идеи на двумерный и трехмерный симплексы.
В данной работе на примере одномерного симплекса исследованы нелинейные "крышки" в двух типичных случаях: когда плотность в средней точке симплекса максимальна и когда плотность в указанной точке минимальна. Нетрудно предположить, что сохранение весового баланса и соблюдение симметрии не гарантируют однозначность решения задачи Мёбиуса о смещении барицентра. Каждая модель по своему реагирует на появление эксцентриситета. Ниже предпринята попытка с помощью симметричных пар кривых ввести на неоднородном симплексе локальные координаты, обладающие характерными свойствами барицентрических координат Мёбиуса.
Основная часть
Не уменьшая общности, рассмотрим симплекс единичной длины [0; 1]. Как известно, Б-координаты одномерного симплекса имеют вид:
¿0 (х ) = 1 - х; Ьх(х ) = х. (1)
Эта симметричная пара (рис. 1) обладает следующими характерными свойствами:
1 (11 = к 0 < Ц(х)< 1; ^¿(х) = 1; 11(хк) = \'^1 (2)
1=0 I0, ^ к,
где - номер функции, к - номер узла.
Следующий шаг состоит в подборе нелинейных симметричных пар, обладающих свойствами (2). Из множества подходящих пар мы выбрали две модели (рис. 1).
0 1
Рис. 1. Три модели симметричных пар
Модель Е:
Е0 (х) = 2х3 - 3х2 +1; Е1 (х) = -2х3 + 3х2 . (3)
Модель Я :
Яо (х) =
008
пх 2
. Пх ТТ
81П--+ 008-
22
&(х ) =
81П
Тх 2
Тх Тх
81П--+ 008-
22
(4)
Примечание. В теории интерполяции функция-"крышка" составляется из двух "полукрышек" путем соединения пары в наивысшей точке. В задаче Мёбиуса удобнее соединить "полукрышки" в средней точке (в центре симметрии). Понятно, что во всех случаях барицентр симплекса совпадает с геометрическим центром симплекса. Поэтому все модели Ь, Е и Я дают одинаковое поузловое распределение единичной массы симплекса: у о =71 = 0,5. Теперь сместим центр тяжести симплекса в точку х = 0,25 , вычислим нагрузки в узлах 0 и 1. Результаты таковы:
для модели Ь : Уо = 0,75, = 0,25; для модели Е : ^0 = 0,84, У1 = 0,16; для модели Я : у0 = 0,71, = 0,29.
Заметим, что при одинаковой средней плотности 1, плотность неравномерного распределения модели Е достигает максимума 1,5 в центре отрезка, а плотность распределения модели Я в той же точке достигает минимального значения 0,785.
В формулах (2) легко усмотреть четкий вероятностный смысл Б-координат Мёбиуса. Так, функция Ь (х) определяет вероятность случайного попадания точки в [0; х]. Тогда Ь (х) определяет вероятность попадания случайной точки в [х; 1]. В рамках нелинейного закона распределения вероятностей функция 70 (х) ассоциируется с вероятностью попадания случайной точки в [х; 1], ^(х) - с вероятностью попадания случайной точки в [0; х]. Точнее, 7} (х) определяет вероятность попадания случайной ординаты в интервал 7о (х); 1], а 71 (х) - это вероятность попадания случайной ординаты в [0; 71 (х)].
Анализ показывает, что нелинейные "полукрышки" наследуют многие замечательные свойства Б-координат. Кроме отмеченного выше, укажем на то, что любая кривая из симметричной пары разделяет единичный квадрат (рис. 1) на две равновеликие фигуры. Интересно, что площадь прямоугольного треугольника с искривленной гипотенузой можно вычислить точно, пользуясь любой из приближенных формул интегрирования, например, центральных прямоугольников (1 узел), трапеций (2 узла), Симпсона (3 узла), "трех восьмых" (4 узла).
Приведём ещё две симметричные пары нелинейных "полукрышек": модель Н (с ослабленным центром):
Н 0 (х)= , Я,(х )= ^ ; (5)
Ых + V1 - х ых +л/1 - х
модель О (с усиленным центром):
00 (х) = -6х5 +15 х4 -10х3 +1; О1(х)= 6х5 -15 х4 + 10х3 . (6)
Заинтересованный читатель уже понял, что множество симметричных пар не ограничивается приведенными здесь примерами. Среди семи протестированных нами моделей наиболее плотной в центре оказалась модель О , а узловые нагрузки моделей (5) и (6) таковы:
для модели Н : /0 = 0,64, /1 = 0,36; для модели О : У0 = 0,9, /1 = 0,1.
На первый взгляд может показаться, что в барицентрической задаче Мёбиуса возникает парадокс (типа парадокса Бертрана). На самом деле, в каждом конкретном случае "работает" свой закон распределения массы, позволяющий уравновесить неоднородный рычаг при смещении точки опоры.
Выводы
В классической задаче Мёбиуса плотность постоянна, поэтому Б-координаты учитывают только геометрию. Симметричные пары нелинейных функций учитывают не только геометрию, но и физику. Представляет интерес использование нелинейных "полукрышек" в барицентрической задаче Мёбиуса для неоднородной квадратной пластины.
Список использованной литературы
1. Стренг Г. Теория метода конечных элементов / Г. Стренг, Дж. Фикс. — М. : Мир, 1977. — 350 с.
2. Марчук Г. И. Введение в проекционно-сеточные методы / Г. И. Марчук, В. И. Агошков. — М. : Наука, 1981. — 416 с.
3. Балк М. Б. Геометрия масс / М. Б. Балк, В. Г .Болтянский. — М. : Наука, 1987. — 160 с.
