МОДЕЛИРОВАНИЕ ПРИРОДОПОДОБНЫХ КРИВЫХ Текст научной статьи по специальности «Физика»

УДК 514.851

DOI: 10.30987/1999-8775-2021-6-11-22

В.А. Короткий, И.Г. Витовтов

МОДЕЛИРОВАНИЕ ПРИРОДОПОДОБНЫХ КРИВЫХ

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

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

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

V.A. Korotky, I.G. Vitovtov

NATURE-LIKE CURVE MODELING

A physical spline is called an elastic rod the cross- section dimensions of which are rather small as compared with the length and radius of its axis curvature. Such a rod when passing through specified points obtains in natural way a nature-like shape characterized with minimum energy of inner stresses and minimum mean curvature. A search for the equation of elastic line is a difficult mathematical problem having no elementary solution.

The work purpose: the development of the experimental-rated procedure for modeling a nature-like elastic curve passing through complanar points specified in advance.

The investigation methods: methods of piece-cubic interpolation based on the application of polynomial splines and compound curves specified by parametric equations. In the paper there are considered polynomial and parametric methods of the geometric modeling of the physical spline passing through the points specified in advance. The elastic line of the physical spline is obtained experimentally.

The investigation results: it is shown that unlike a polynomial model a parametrized model on the basis of Fergusson curve gives high accuracy of approximation if in basic points there are specified tangents to the elastic line of the physical spline with large deflections.

Novelty: there is offered a simplified method for the computation of factors of an approximating spline allowing the substitution of the 2n system of nonlinear equations (smoothness conditions) by the successive solution of n systems of two equations.

Conclusions: for the modeling of nature-like curves with large deflections there is offered the application of Fergusson cubic spline passing through specified points and touching the specified straight lines in these points. The error of the modeling of the natural elastic line with free ends at n=5 does not exceed 0.4%.

Key words: affinating compression, cubic polynomial, vector derivative, Fergusson spline, graphic differentiation.

Введение

Природоподобной кривой называют линию, образованную осью упругого стержня, проходящего через наперед заданные точки. Предполагается, что размеры поперечного сечения такого стержня весьма малы по сравнению с длиной и радиусом кривизны его оси. Такую кривую называют "физическим сплайном".

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

ей внутренних напряжений и минимальной средней кривизной. Такие формы широко используются в современном архитектурно-строительном проектировании [1, 2, 3].

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

Для больших прогибов решение принципиально усложняется. Согласно [5], оно сводится к вариационной задаче, не имеющей элементарного решения. Поэтому моделирование природоподобной кри-

Постановка задачи

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

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

ность и практическую значимость поставленной задачи [6, 7].

Здесь рассматриваются два способа аппроксимации: 1) на основе кубических парабол [8]; 2) с помощью кубического параметрического сплайна, составленного из кривых Фергюсона [9, 10].

В качестве исходных данных использованы результаты экспериментального моделирования физического сплайна, проходящего через п = 5 компланарных точек.

Экспериментальное моделирование физического сплайна

На планшете закреплена масштабно-координатная бумага. В опорные точки А(1.2, 0.9), 5(2.1, 2.3), С(2.8, 4), П(5, 3.2), £(5.5, 2.1) вбиты тонкие булавки с нанесенной на них смазкой (для снижения трения). Масштаб планшета 70:1. На планшет уложена гибкая металлическая линейка, огибающая булавки. Координаты опорных точек назначены в зоне упругой деформации материала линейки.

Модель физического сплайна со свободными концами. Моделирование физического сплайна со свободными кон-

цами выполнено в двух вариантах. Согласно первому варианту, металлическая линейка проходит через указанные опорные точки А, В, С, П, Е. Согласно второму варианту, исходный набор опорных точек подвергается аффинному сжатию к оси х с коэффициентом сжатия 0.5. Новые опорные точки А', В', С', П', Е' занимают среднее положение между исходными точками А, В, С, П, Е и осью х (рис. 1). В обоих вариантах свободные концы линеек принимают свое естественное положение.

Рис. 1.Исходная и аффинно-сжатая упругие линейки со свободными концами (фото и чертеж)

Используя алгоритм графического дифференцирования [11], определяем значения производных (тангенсов углов

наклона касательных) исходного физического сплайна ЛВСБЕ в его опорных точках:

tgSA = 1,134, tgSв = 2,840, tgSc = 0,985, tgSD = -1,815, tgSE = -2,460.

Выполнив аффинное сжатие {.Л, В, С, Б, Е} ^ {Л', В', С', Б', Е'} с коэффициентом сжатия 0.5, замечаем, что касательные к исходной и сжатой линейкам в конечных точках Л, Л' и Е, Е' пересекаются на оси х (рис. 1), следовательно, первые производные на концах сжатого физического сплайна уменьшились ровно в два раза.

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

tgSA = 0,567, tgSв = 1,205; tg5c = 0,811; tgSD = -1,00; tg5E = -1,230.

Производные в опорных точках исходной и сжатой кривой изменились не в 2 раза, а в 2,36 (точки В~В) в 1,21 (точки С~СГ), в 1,815 (точки Б~БГ), то есть деформация физического сплайна { Л, В, С, Б, Е} ^ {Л', В', С', Б', Е'} имеет в целом

Модель физического сплайна с защемленными концами. На рис. 2 представлены макеты физических сплайнов с защемленными концами. Исходный сплайн ЛВСБЕ характеризуется касательными в опорных точках:

нелинейный характер.

tgдA = 0,2; tgдв = 5,16; tgдc = 0,83; tgдD = -4,36; tgдE = -0,3. Аффинно-сжатый сплайн Л'В'С'Б'Е' с коэффициентом сжатия 0,5:

tg8A = 0,1; tg8в = 1,528; tg8c = 0,721; tg8D =-1,481; tg8E =-0,15.

Как и для сплайна со свободными концами, в случае защемленных концов наблю-

дается нелинейная деформация упругой линии физического сплайна.

Рис. 2. Исходная и аффинно-сжатая упругие линейки с защемленными концами (фото и чертеж)

Полиномиальная аппроксимация

Будем искать уравнения сегментов ЛВ, ВС, СБ, БЕ в виде кубических полиномов:

Бав(х) = а + ь1(х -хА) + С1(х-хА)2 + ^(х-хА)3, х е [хА, хв];

Яве (х) = а2 + Ь2(х - хв ) + с2(х - хв )2 + а2(х - хв )3. х е [хв > хс ];

(х) — а^ ^ (х хс ) ^ с^ (х хс ) ^ (х хс ) , х е [ х^, х^ ];

(х) = а4 + Ъ4(х-xD) + с4(х-xD)2 + й4(х-xD)3, х е [xD, хЕ]. Здесь хл = 1,2; хв = 2,1; хс = 2,8; хл = 5; хя = 5,5.

Составная кривая (1) - кубический полиномиальный сплайн. В отличие от реального (физического) сплайна его называют математическим сплайном. Как и физический сплайн, сплайн (1) проходит через заданные точки А, В, С, П, Е, удовлетворяет заданным граничным условиям и имеет второй порядок гладкости (непрерывно изменяющуюся кривизну). Тем не менее, как будет показано ниже, форма физического и математического сплайнов может существенно различаться.

а1 = Уа = 0,9 а2 = Ув = Таким образом, в системе уравнений (1) содержится не 16, а всего 12 подлежащих определению коэффициентов.

Уравнения (1) содержат 16 неизвестных коэффициентов а, Ъ, С, й = 1.4) ,

которые определяются из трех групп условий:

- инцидентности сегментов точкам А, В, С, П, Е;

- отсутствия изломов (равенство первых производных на стыках сегментов);

- непрерывного изменения кривизны (равенство вторых производных в стыковых точках).

Заметим, что из (1) непосредственно следует:

2,3; а3 = ус = 4; а4 = ул = 3,2. (2)

Условия инцидентности имеют вид:

У в = Уа + ъхК+ С1 К + ¿1 К, Ус = У в + ъ2к2 + с2к2 + й2к2 , у0 = ус + ЪъИъ + сК + ¿ЪИ1,

Уе = У о + КК + С4К + йК, где — х£> хд, — х^ х', ^^ — х^ х, — х^ х^. Условия отсутствия изломов имеют вид:

Ъ + 2ск + 3йхК = Ъ2, Ъ2 + 2с2и2 + 3 й2Ъ2 = Ъ3,

Ъ3 + 2с3к3 + 3йъЬ2 = Ъ4. Условия непрерывного изменения кривизны имеют вид:

С + 3йхК = с2 ,

С 2, ^ 3^2 — С ,

^ + — ^.

Десять уравнений (3), (4), (5) содержат двенадцать неизвестных коэффициентов Ъг, с, й (* = 1,.., 4). Два дополнительных условия назначаются либо приравниванием нулю вторых производных в концевых точках А, Е (свободные концы), либо указанием первых производных в точках А, Е (защемленные концы).

Полиномиальная модель физического сплайна со свободными концами. Вторая производная

$ав (х) = 2с + 6й (х - ^) сегмента АВ в точке А должна быть равна нулю. Отсюда следует, что С1 = 0.

Вторая

(3)

(4)

(5)

производная

S'de (х) — 2с4 + 6d4 (х - хд ) сегмента DE в точке E также должна быть равна нулю. Отсюда следует, что

S'De (хЕ ) — c4 + 3d4h4 — 0. (6)

Решая (с помощью системы компьютерной алгебры SMath Studio Desktop) систему линейных уравнений (3), (4), (5), (6) с учетом ci = 0, находим коэффициенты Ь, ci, d • Получаем математический сплайн со свободными концами, проходящий через наперед заданные точки A, B, C, D, E (рис. 3а):

ЯАВ (х) = 0,9 +1,2174(х -1,2) + 0,4175(х -1,2)3, х е [1,2; 2,1]; Явс (х) = 2,3 + 2,2319(х - 2,1) +1,1272(х - 2,1)2 -1,2089(х - 2,1)3, х е [2,1; 2.8]; Ясо (х) = 4 + 2,0329(х - 2,8) -1,4114(х - 2,8)2 + 0,1464(х - 2,8)3, х е [2,8; 5]; Я (х) = 3,2 - 2,0516(х - 5) - 0,4452(х - 5)2 + 0,2968(х - 5)3, х е [5; 5,5].

(7)

Сравнивая физический (поз. 1) и математический (поз. 1') сплайны (рис. 3а), замечаем существенное различие на интервалах ВС и СБ. Это различие уменьшается для аффинно-сжатых кривых (поз. 2 и 2'), что подтверждает известное положение: математический сплайн, составленный из кубических полиномов, тем ближе

Б'лв (ха ) = ^^^^ = 0,2 = Ъ Решая систему линейных уравнений (3), (4), (5), (8), находим коэффициенты

к физическому сплайну, чем меньше прогибы физического сплайна.

Полиномиальная модель физического сплайна с защемленными концами.

Тангенсы углов заделки равны: tgSA= 0,2, tg5E =-0,3. Граничные условия приобретают вид:

1, ^ (хе) = ^^^^ = -0,3 = ¿4 + 2с,к, + 3й4к2. (8) Ъ, С, di. Получаем уравнение кубического сплайна с защемленными концами:

Рис. 3. Полиномиальная модель физического сплайна: а - свободные концы; б - защемленные концы (1, 2 - эксперимент, 1', 2' - расчет)

^(х) = 0,9 + 0,2(х -1,2) + 2,0181(х -1,2)2 - 0,5688(х -1,2)3, х е [1,2; 2,1]; Явс(х) = 2,3 + 2,4504(х - 2,1) + 0,4824(х - 2,1)2 - 0,7336(х - 2,1)3, х е [2,1; 2,8]; Ясо (х) = 4 + 2,0473(х - 2,8) -1,0583(х - 2,8)2 + 0,0171(х - 2,8)3, х е [2,8; 5]; ЯОЕ(х) = 3,2 - 2,8572(х-5) -1,1710(х-5)2 + 4,971(х- 5)3, х е [5; 5,5].

(9)

Полиномиальный сплайн (9) проходит через заданные точки Л, В, С, Б, Е и имеет заданные производные в граничных точках Л, Е. Тем не менее, результаты моделирования (рис. 3б) показывают заметное различие физического (поз. 1, 2) и ма-

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

Параметрическая аппроксимация

Векторное уравнение сегмента параметрически заданной кубической кривой, проходящей через точки Ко, К и касаю-

щейся в этих точках прямых то, Т1, имеет вид:

г(г) = г(0)(1 - 3г2 + 2г3) + г(1)(3г2 - 2г3) + г(0)г(1 - г)2 + г(1)(-г2 + г3),

(10)

где параметр t изменяется в диапазоне [0, 1]. Векторы г(0), г(1) указывают положение опорных точек Ко, Я1, а направления векторных производных Г(0), Г(1) при t=0 и t=1 совпадают с направлениями касательных Т0, Т1 в опорных точках.

Кубическая кривая (10), предложенная Фергюсоном в 60-х годах 20 века, описывает двухпараметрическое семейство кубических кривых, инцидентных точкам К0, К и касающихся в этих точках прямых Т0, Т1. Форма любой кривой этого семейства определяется величинами | Г(0) |, | Г(1) |. При увеличении | Г(0) | кривая г(г) теснее прилегает к касательной Т0, а при увеличении | Г(1) | теснее прилегает к Т1 [12, 13].

Условие гладкости в форме Фергюсона. Согласно выражению для расчета вектора кривизны

г(0хг(0

К(г) = ■

(11)

|г(0|3 '

условие гладкого соединения кубических сегментов (10) имеет вид

'(1) х Г(1)(1) _ Г(2)(0) х Г(2)(0)

г(1)/

(12)

ги (1) = 'Т, г где Т - единичный вектор общей касательной, а коэффициенты ^(1) , ^(2) - модули первых производных в точке стыка. Подставляя (13) в (12), получаем условие гладкости:

;(2)/

(0) = (щ(2)/ щи))2Г(1)(1) .

аьг-а).

(14)

- п.,,(1)

| Г(1)(1)|3 |г(2)(0)|3

где надстрочный символ в скобках означает номер сегмента. Параметр I в уравнениях сегментов г(1)(г), г(2)(г) изменяется в диапазоне [0, 1] независимо от номера сегмента.

Учитывая, что в стыковой точке оба сегмента должны иметь общую касательную, представим первые производные, входящие в (12), в виде (2)(0) = м?(2)Т, (13)

Пусть требуется провести гладкую кривую через опорные точки К0, К1, К2, заданные радиус-векторами Г0, Г1, Г2. Будем искать решение в виде кривой, состоящей из сегментов кубических кривых вида (10). Записывая векторные производные в виде

Г (0) = щ,Т>, Г(1)(1) = Г(2)(0) = Г(2)(1) = w2T2

и подставляя в (10), получаем уравнения сегментов г(1)(^) и г(2)(^):

г(() (г) = г (1 - 3г2 + 2г3) + г (3г2 - 2г3) + щТг(1- г)2 + щ(1)Т (-2 + г3), г е [0,1] г(2) (г) = г (1 - 3г2 + 2г3) + г2(3г2 - 2г3) + щ(2)Тг(1 - г)2 + щТ2(-г2 + г3), г е [0,1]. Здесь Т0, Т1, Т2 - единичные векторы касательных Т0, Т1, Т2, а коэффициенты - модули первых производных. В общем случае ^1(1) Ф ^1(2).

Дважды дифференцируя выражение (16) и подставляя t=1, получаем значение второй производной сегмента К0-К1 в его конечной точке К1:

Г(1)(1) = 6г0 - 6г + 2щТ + 4ща)Т1.

(15)

(16) (17)

Дважды дифференцируя выражение (17) и подставляя t=0, находим вторую производную сегмента К1-К2 в его начальной точке К1:

г(2) (0) = -6г + 6г2 - 4щ(2)Т - 2щТ2.

,(2)п

Подставляя Г()(1) и Г( )(0) в (14), получаем условие гладкого соединения сегментов:

где в общем случае Л = щ(2) / щ(1).

Конструирование сплайна

Фергюсона. Через указанные на плоскости точки А, В, С, П, Е требуется провести сплайн Фергюсона - гладкую кривую, со-

3г2 + 3г (Л -1) - 3гЛ = щТЛ + 2щ (1 + Л )Т + щТ, (18)

ставленную из сегментов вида (10). В концевых точках А, Е конструируемой кривой заданы касательные та, те. Сплайн имеет вид:

Гав (г) = Га (1 - 3t2 + 2t*) + Гв (3t2 - 2tл) + ^(1 -1)2 + rB1)(-t2 +13), t е [0,1],

г вс ^) = гв (1 - 3t2 + 2t3) + гс (3t2 - 2t3) + г(,2)t (1 -1 )2 + Г£2) (-t2 +13), t е [0,1],

гсв (() = гс (1 - 3t2 + 2t3) + гв (3t2 - 2t3) + г^(1 -1)2 + Г(3) (-t2 +13), t е [0,1],

(19)

гве ^) = г0 (1 - 3t2 + 2t3) + ге (3t2 - 2t3) + г^(1 -1)2 + г£ (-2 +13), t е [0,1].

Здесь Га = ^Т, Ге = ^Т - векторные производные в граничных точках Л, Е. Надстрочный индекс в скобках означает номер сегмента. Направления единичных векторов Тл, TE заданы касательными тл,

гвв1} = ^рТв, гв2-1 = , гс2 = ^^^тс,

где Хв, Хс, Хб - единичные касательные векторы в точках В, С, Б. Производные (2о) определяются из условий гладкости.

ТЕ. Модули Ж4, мЕ векторных производных определяют "степень прилегания" сплайна (19) к этим касательным. Остальные производные, входящие в (19), определяются выражениями:

3)_,,.(3)т - (3) (3)т - (4) _ ,,.(4)т (20)

= ^с Тс, ГD = WD TD, ГD = WD TD, (20)

Согласно (18), условия гладкости имеют вид:

или, с учетом (20):

3гс + 3гв(4 -1) -3Га4 = <Та4 + 2^в2)(1 + лв)Тв + ТК® /яс, + 3Гс (% -1) - 3Гв4 = Ц2)Тв4 + 2wC3) (1 + яс )Тс + TDwD4) / л 3Ге + 3Г^(% -1) -3Гс4 = ^.3)Тс4 + + лD)TD + <Те,

3гс + 3гв (% -1) - 3Га% = ГаЛ + 2гв2) (1 + Яв) + Г^ / Яс,

+ 3Гс (Яс -1) - 3гвЛс = гв2)лс + 2гсз) (1 + Яс) + Г^ / ЛD,

3Ге + 3^ (Л -1) - 3Г% = гСЛ + 2^4) (1 + ЯD ) + Ге ,

(21)

(21')

(2)

ЯD = / К

где % = w(2) / w(1), % = / стыковых точках В, С, Б).

Условия (21), (21') могут трактоваться двойственным образом. Если зафиксировать скалярные величины Яв, Яс, Яп, то

система уравнений (21') становится системой линейных уравнений относительно неизвестных векторных величин Гд2), Г<13), Г(4) . Если же задать направления единичных касательных векторов Хв, Хс,

(3)

(отношения модулей первых производных в

Хэ, то система (21) становится системой нелинейных уравнений относительно ска-

лярных величин Яв, Яс, Яв, , ,

w

(4)

Фиксируя в (21') скалярные величи-получаем условия гладко-

ны Яв, Яс, Я сти:

2(1 + Яв Уж + / Яс = 3хс + 3хв (Яв - 1) - 3хАЯв - Я^А cos 3А ,

2(1 + Яв ув2 + с / Яс = 3ус + 3 у в (Я2в - 1) - 3уаЯ2в - Я2вWA sin За , Я^вх + 2(1 + Яс )гсх + ^ / ЯD = 3XD + 3хс (Яс -1) - 3хв яс , Я2с4? + 2(1 + Яс у® + / ЯD = 3УD + 3 Ус (Л -1) - 3ув яс,

+ 2(1 + Я Уш = 3хЕ + 3XD (Я - 1) - 3хсЯЯП - WE cos 3Е ,

Я2DíCr> + 2(1 + ЯD = 3уе + 3УD Я -1) - 3усЯ% - WE sin Зе ,

(22)

где 5л и 5е - углы наклона касательных тл, тЕ к оси х.

В условиях (22) можно варьировать пятью управляющими параметрами: модулями первых производных мл, ме в граничных точках Л , Е и величинами Яв, Яс, Я . При изменении этих величин меняется форма сплайна (19).

Пример 1. Найти аппроксимирующую функцию для физического сплайна ЛВСБЕ со свободными концами (рис. 1), варьируя величинами Яв, Яс, Яв .

Положив в (22) мл = ме = 1, ^в = 1, Хс=2, Аб=0.5, находим компоненты векторных производных Г™, г®, г® , г!3, Г^3), Г^4). Раскладывая (19) по осям х, у, получаем

скалярные функции вида х = х^), у = у(), позволяющие построить сплайн Фергюсо-на (рис. 4а). Погрешность аппроксимации достигает 4 %.

Пример 2. Найти аппроксимирующую функцию для физического сплайна АВСПЕ с защемленными концами (рис. 2), варьируя величинами Лв, Лс, .

Подставляя в (22) координаты опорных точек А, В, С, П, Е и значения гgSA = 0.2, г%5Е = -0.3, положив Ж4= ^е= 1, Ав = Ап =1, Ас = 2, получаем систему линей-

ных уравнений, из которой находим компоненты векторных производных, входящих в коэффициенты сплайна (19). Раскладывая (19) по осям х, у, получаем скалярные функции вида х = х^), у = у(), позволяющие построить сплайн. На рис. 4, б представлен физический сплайн (поз. 1) и сплайны Фергюсона при Ас = 2 (поз. 2) и Ас = 3 (поз. 3). Погрешность аппроксимации более 10 %.

а)

б)

Рис. 4. Аппроксимация физического сплайна: а - свободные концы; б - защемленные концы (поз. 1 - эксперимент, поз. 2, 3 - расчет)

Примечание. В условиях гладкости (21') можно произвольным образом задавать модули WA, ^е первых производных в крайних точках А, Е, а также величины Ав, Ас, Ап. Следовательно, имеется 5 управляющих параметров, влияющих на форму сплайна (19). Задача подбора этих параметров для аппроксимации реальной упругой линии может быть решена лишь путем многочисленных "численных экспериментов". К тому же нет уверенности, что найденный таким образом вариант аппрок-

2.

Z =

-2 =

z = w(2)

yvC 5

симирующего сплайна является наилучшим.

Таким образом, условия гладкости (21'), содержащие пять варьируемых параметров, не могут быть рекомендованы для моделирования физического сплайна.

Сплайн Фергюсона с заданными касательными. Рассмотрим алгоритм построения сплайна (19), проходящего через заданные опорные точки и касающегося в этих точках заданных прямых. Введем обозначения

z

= w(3)

C

z5 = W>5

z6 = <

Подставляя в (21), получаем систему нелинейных алгебраических уравнений относительно неизвестных zi, ..., Z6:

3ä2Zj2 + (3h - wA cos ôA )z2 - 2 cos ôB (z\z2 + z2z ) - z\z3 cos ôc = 0;

3Ay2Zj + (3Ay - wA sin 5a )z2 - 2 sin ôB (Zj z2 + z2 z ) - zx z3 sin 5C = 0;

22 22 1 z2 + " '

2 , 2

(23)

3^z3 + 3h2z4 - z4z2 cosôB -2cosôc(z4z3 + z4z3) -zz3 cosôD = 0; 3 Ay3z32 + 3Ay2z^ - z2z2 sin ôB - 2 sin öc(z^ + z2z3) - zsz\ sin ôD = 0; 3h4z52 + 3hz62 - z4z2 cos ôc - 2 cos ôD (z6z52 + z2z5 ) - z\wE cos ôE = 0; 3AyAz\ + 3Ay1>z\ - z4z2 sin ôc - 2 sin ôD (z6z52 + z2z ) - z\wE sin öE = 0,

где

h = xb - xa ; h2 = xc - xb ; h3 = хо - xc ; h4 = xe - хо ; Ay = Ув- y a =; A^ = Ус- Ув ; лУэ=y о - Ус ; аУ4 = Уе - Уо •

Задавая углы дА, ^ наклона касательных и модули WA, ^е векторных производных Г л = щ ТА, ГЕ = щ Т в концевых точках А, Е и решая нелинейную систему уравнений (23), находим модули ¿ь.^б производных в стыковых точках В, С, П, что позволяет полностью определить коэффициенты сплайна (19).

Рассмотренный алгоритм предполагает итерационное решение нелинейной системы (23), что может стать серьезной вычислительной проблемой. Решение

можно упростить, используя разработанный авторами алгоритм подстановки.

Алгоритм подстановки. Зафиксируем первый сегмент АВ сплайна (19), произвольно задав модули первых производных щ и щ(1) в точках А, В. Варьируя величинами щ и щ(1) , добиваемся удовлетворительного совпадения расчетного сегмента АВ с соответствующим сегментом физического сплайна. Подставляем найденные

(1)

величины wA и w уравнений системы (23):

21 в первую пару

3h2z2 + (3h - wA cos ôA )z2 - 2 cos ôB (zx z2 + z2 z ) - zx z cos 5C = 0, 3Ay2z2 + (3 Ay - wA sin ôA)z\ - 2 sin ôB (z,2z2 + z2z) - z\z sin ôc = 0.

(24)

Решая (24), находим неизвестные ^ = щ/, 2 = щс и подставляем их во вторую пару уравнений системы (23):

3hz32 + 3^z42 - z2z2 cos ôB- 2 cos ôc(zAz\ + z2z) - zz2 cos ôD = 0,

3Ay3z32 + 3Ay2z2 - z\z2 sin ôB - 2 sin ôc(zAz\ + z2z) - zz32 sin ôD= 0. _™(3) „ _ «,(3)

(25)

Решая (25), находим неизвестные= щЛ 2 = щ> и подставляем их в третью пару уравнений системы (23):

3h4z\ + 3hz62 - z4z2 cos ôc - 2 cos ôD ( z6z2 + z I z5 ) - z2 wE cos ôE = 0, 3Ay4z2 + 3Ay3z2 - z^ sin ôc - 2sin ô (z6z2 + z,2z ) - z2wE sin ôE = 0.

(26)

wE

Решая (26), находим неизвестные 26, г^е. Все величины, входящие в коэффициенты сплайна (19), определены.

Примечание. Алгоритм подстановки позволяет заменить решение системы шести нелинейных уравнений (23) последовательным решением трех пар нелинейных уравнений (24), (25), (26).

Управление формой сплайна Фергюсона с заданными касательными

Раскладывая векторные уравнения (19) по осям х, у, получаем скалярные па-

раметрические Фергюсона:

уравнения

сплайна

xab

Уав хвс Увс

xcd

Усо

xde

Уое

= ^ (1 - 3t2 + 2t3) + xB (3t2 - 2t3) + w^ cos SAt(1 -1)2 + w™ cos SB (-t2 +13), = y^ (1 - 3t2 + 2t3) + ys (3t2 - 2t3) + w sin SAt(1 -1)2 + w(1) sin SB (-t2 +13), = xs (1 - 3t2 + 2t3) + xc (3t2 - 2t3) + wf cos ^t(1 -1)2 + w(,2) cos (-t2 +13) = yB (1 - 3t2 + 2t3) + yc (3t2 - 2t3) + w(2) sin SBt(1 -1)2 + w(2) sin Sc (-t2 +13),

(27)

= xc (1 - 3t2 + 2t3) + xD (3t2 - 2t3) + w(3) cos Sct(1 -1)2 + w(3) cos SD (-t2 +13),

„(3),

= yc (1 - 3t2 + 2t3) + yD (3t2 - 2t3) + w(3) sin Sct(1 -1)2 + sin 3D (-t2 +13),

,,(3),

= xD(1 -3t2 + 2t3) + xB(3t2 -2t3) + w(4) cosSDt(1 -1)2 + wB cosSE(-t2 +13),

= yD (1 - 3t2 + 2t3) + yB (3t2 - 2t3) + w£° sin t(1 -1)2 + wE sin SE (-t2 +13).

В уравнениях (27) можно варьировать модулями мл, ме первых производных в крайних точках Л, Е. Если используется алгоритм подстановки, то можно произвольно задавать модули и w(1) первых производных в граничных точках Л , В первого сегмента. Таким образом, имеется всего два управляющих параметра, влияющих на форму сплайна, проходящего че-tg5A = 1,134, tg5в = 2,840; tg5c = вычисляем синусы и косинусы углов наклона касательных За , 5В, Зс, Зв, Зе и подставляем их в условия гладкости (23).

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

Пример 3. Моделирование физического сплайна ЛВСБЕ со свободными концами (рис. 2).

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

0,985; tgöD = -1,8154 tgöE = -2,460,

Положив варьируемые параметры WA=WE=1, решаем систему уравнений (23) методом Ньютона (SMath Studio):

* = w(1) = 2,175; z. = w(2) = 1,784; z.

= w(2) = 1,845; z = w(3) = 2,0335;

yC

- ™(3) = 4,850; zÄ = w<4) = 1,1935.

(28)

Z5 = w D

'6 ~ "О

Подставляя (28) в (27), получаем уравнения аппроксимирующего сплайна, отличающегося от физического сплайна менее, чем на 0,4% (рис. 5а).

tgЗA = 0,2; tgЗв = 5,16; tgЗc = 0,83; ^ =-4,36; tgЗE =-0,3 .

Пример 4. Моделирование физического сплайна с защемленными концами (рис. 3). Используем экспериментальные значения:

У

3

2

Лс л

А и D

jв Vi

/ в

v

2

а)

У

\tF

и 5 х

1 <\г О

Г'с//С \2 1 О

В

1

"Я Z в

2

б)

Рис. 5. Физический сплайн (поз. 1) и сплайн Фергюсона (поз. 2): а - свободные концы; б - защемленные концы

5 х

Варьируя модулями w

и w

(1)

в

производ-

ных rA, rB1}

в граничных точках А , В первого сегмента, добиваемся совпадения формы расчетного сегмента АВ с сегментом АВ физического сплайна. В рассматривае-

мом примере требуемое совпадение достигнуто при щ = щ4 = 2. Используя алгоритм подстановки, последовательно решаем системы уравнений (24), (25), (26):

z2 =

= 1,978; z = w(2) = 1,723; z4 = w(3) = 2,493;

-5 "Э ,

Подставляя (29) в (27), получаем скалярные уравнения, позволяющие построить аппроксимирующий сплайн (рис. 5б). По-

z = w(3) = 3,672; zÄ = w(n4) = 1,292; w = 1,356.

(29)

грешность аппроксимации составила менее 1,5%.

Заключение

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

сегментов сплайна) предлагается способ подстановки, позволяющий заменить решение системы 2(п-2) нелинейных уравнений последовательным решением п-2 пар нелинейных уравнений (п - количество базисных точек физического сплайна). На ряде примеров показано, что ошибка моделирования упругой линии физического сплайна со свободными концами не превышает 0,4%, а погрешность моделирования упругой линии с защемленными концами составляет менее 1,5%.

СПИСОК ЛИТЕРАТУРЫ

1. Голованов, Н. Н. Геометрическое моделирование / Н. Н. Голованов. - Москва : Изд-во физико-математической литературы, 2012. - 472 с.

2. Иванов, Г. С. Теоретические основы начертательной геометрии / Г. С. Иванов. - Москва : Машиностроение, 1998. - 157 с.

3. Glaeser, G. Geometrie und ihre Anwendungen in Kunst, Natur und Technik / G. Glaeser. - Moscow : Springer Spektrum, 2014. - 508 pp. - DOI 10.1007/978-3-642-41852-5.

4. Завьялов, Ю. С. Сплайны в инженерной геометрии / Ю. С. Завьялов, В. А. Леус, В. А. Ско-роспелов. - Москва : Машиностроение, 1985. -224 с.

5. Фокс, А. Вычислительная геометрия. Применение в проектировании и на производстве / А. Фокс, М. Пратт. - Москва : Мир, 1982. - 304 с.

6. Конопацкий, Е. В. Вычислительные алгоритмы моделирования одномерных обводов через k наперед заданных точек / Е. В. Конопацкий, А. А. Крысько, А. И. Бумага // Геометрия и графика. - 2018. Т. 6, № 3. - С. 20-32. - DOI: 10.12737/article 5bc457ece18491.72807735.

1. Golovanov, N.N. Geometric Modeling / N.N.

Golovanov. - Moscow: Publishing House of Physi-

cal-Mathematical Literature, 2012. - pp. 472.

7. Любчинов, Е. В. О гладкости стыковки линий и поверхностей при циклографическом моделировании поверхностных форм автомобильных дорог / Е. В. Любчинов, К. Л. Панчук // Вестник Южно-Уральского государственного университета. Серия «Строительство и архитектура». -2020. - Т. 20. - № 1. - С. 52-62. - DOI: 10.14529/build200106.

8. Понтрягин, Л. С. Кубическая парабола / Л. С. Понтрягин // Квант. - 1984. - №3. - С. 10-14, 32.

9. Уокер, Р. Алгебраические кривые / Р. Уокер. -Москва : Либроком, 2009. - 240 с.

10. Шикин, Е. В. Кривые и поверхности на экране компьютера / Е. В. Шикин, Л. И. Плисс. -Москва : Диалог-МИФИ, 1996. - 240 с.

11. Прасолов, В. В. Геометрия / В. В. Прасолов, В .М. Тихомиров. - Москва : Изд-во МЦНМО, 2013. - 336 с.

12. Савелов, А. А. Плоские кривые / А. А. Савелов. - Москва : Либроком, 2009. - 296 с.

13. Короткий, В. А. Кубические кривые в инженерной геометрии / В. А. Короткий // Геометрия и графика. - 2020.- Т. 8. - № 3. - С. 3-24. -DOI: 10.12737/2308-4898-2020-3-24.

2. Ivanov, G.S. Theoretical Fundamentals of Descriptive Geometry / G.S. Ivanov. - Moscow: Mechanical Engineering, 1998. - pp. 157.

3. Glaeser, G. Geometrie und ihre Anwendungen in Kunst, Natur und Technik / G. Glaeser. - Moscow : Springer Spektrum, 2014. - 508 pp. - DOI 10.1007/978-3-642-41852-5.

4. Zaviyalov, Yu.S. Splines in Engineering Geometry / Yu.S. Zaviyalov, V.A. Leus, V.A. Skorospelov. -Moscow: Mechanical Engineering, 1985. - pp. 224.

5. Fox, A. Computational Geometry. Application in Design and Production / A.Fox, M. Pratt. - Moscow: Mir, 1982. - pp. 304.

6. Konopatsky, E.V. Computational algorithms of one-dimensional outlines through k points specified in advance / E.V. Konopatsky, A.A. Krysko, A.I. Bumaga // Geometry and Graphics. - 2018. Vol.6, No.3. - pp. 20-32. - DOI: 10.12737/article_5bc457ece18491.72807735.

7. Lyubchinov, E.V. On smoothness of line and surface joint at cyclographic modeling of highway surface forms / E.V. Lyubchinov, K.L. Panchuk // Bul-

Ссылка для цитирования:

letin of South-Urals State University. "Construction and Architecture" Series. - 2020. - Vol. 20. -No.1. - pp. 52-62. - DOI: 10.14529/build200106.

8. Pontryagin, L.S. Cubic Parabola / L.S. Pontryagin // Quantum. - 1984. - No.3. - pp. 10-14, 32.

9. Wocker, R. Algebraic Curves / R. Wocker. - Moscow: Librokom, 2009. - pp. 240.

10. Shikin, E.V. Curves and Surfaces on Computer Display / E.V. Shikin, L.I. Pliss. - Moscow: Dia-logue-MIFI, 1996. - pp. 240.

11. Prasolov, V.V. Geometry / V.V. Prasolov, V.M. Tikhomirov. - Moscow: MCNMO Publishers, 2013. - pp. 336.

12. Savelov, A.A. Flat Curves / A.A. Savelov. - Moscow: Librokom, 2009. - pp. 296.

13. Korotky, V.A. Cubic curves in engineering geometry / V.A. Korotky // Geometry and Graphics. -2020. - Vol.8. - No.3. - pp. 3-24. -DOI: 10.12737/2308-4898-2020-3-24.

Короткий, В.А. Моделирование природоподобных кривых /В.А. Короткий, И.Г. Витовтов //Вестник Брянского государственного технического университета. - 2021. - № 6. - С. 11 - 22 . БOI: 10.30987/1999-8775-2021-611-22.

Статья поступила в редакцию 16.03.21. Рецензент: д.т.н., доцент Оренбургского государственного университета, Рассоха В.И.,

член редсовета журнала «Вестник БГТУ». Статья принята к публикации 25.05.21.

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

Короткий Виктор Анатольевич, профессор кафедры «Инженерная и компьютерная графика», Южно-Уральский государственный университет, email: ospolina@mail.ru.

Korotky Victor Anatolievich, Prof. of the Dep. "Engineering and Computer Graphics", South-Urals State University, e-mail: ospolina@mail.ru.

Витовтов Игорь Георгиевич, к. физ.-мат. н., доцент кафедры «Естественнонаучные дисциплины», Челябинский институт путей сообщения, e-mail: vitovtovig@gmail.com.

Vitovtov Igor Georgievich, Can. Sc. Phys.-Math., Assistant Prof. of the Dep. "Natural Science Subjects", Chelyabinsk Institute of Railway Transport, e-mail: vitovtovig@gmail.com.