Геометрическое моделирование поверхностных форм отгона виража автомобильной дороги на основе циклографического отображения

Бесплатный доступ

В теоретических исследованиях по проектированию автомобильных дорог ось дороги принимается в качестве базового элемента, определяющего пространственные формы других элементов дороги: проезжей части, обочины, откосов. Исходя из этого принципа формообразования элементов дороги, в настоящей работе предложены математические модели формообразования, основанные на циклографическом методе моделирования пространства и его объектов. Формообразование поверхностных форм элементов дороги на её переходном (отгонном) участке и на вираже существенно отличается от формообразования на прямолинейном участке рядом специфических особенностей, затрудняющих и усложняющих процесс теоретического и практического формообразования. В работе показана возможность построения математических моделей формообразования для отгонного участка на основе циклографического отображения оси дороги. Приведён пример, в котором выполнен подробный разбор с визуализацией каждого этапа формообразования. Результаты работы могут быть использованы для разработки современных специализированных САПР автомобильных дорог.

Еще

Геометрическое моделирование, циклографическое отображение, автомобильные дороги, виражи, отгоны виражных участков

Короткий адрес: https://sciup.org/147232128

IDR: 147232128   |   DOI: 10.14529/build190109

Текст научной статьи Геометрическое моделирование поверхностных форм отгона виража автомобильной дороги на основе циклографического отображения

Геометрические параметры автомобильных дорог играют важнейшую роль в определении их транспортно-эксплуатационных качеств. Если геометрические параметры дороги не отвечают действующим нормам и требованиям, то эти ее недостатки нельзя исправить лучшими покрытиями полотна [1]. Автомобильная дорога является сложным строительным объектом и состоит из нескольких отличающихся друг от друга участков: прямолинейного в плане и виражного, при этом виражный участок, в свою очередь, делится на собственно вираж (круговой) и отгонные участки, служащие для плавного соединения прямолинейного в плане и кругового участков (рис. 1). По статистике наиболее аварийными участками автомобильных дорог являются различные поворотные участки, поэтому данные участки должны в максимальной степени отвечать установленным требованиям и нормам [2].

В современном проектировании и строительстве автомобильных дорог применяются различные подходы к их математическому моделированию [1, 3–7]. Все они ориентированы на рекомендации, приведенные в государственных нормах и правилах [2], согласно которым прямолинейный в плане участок дороги должен содержать дискретный каркас из типовых поперечных профилей, а на круговом участке дороги (вираже) – дискретный каркас из поперечных профилей в виде прямолинейных отрезков неизменного уклона. Известны математические модели, в которых автомобильная дорога рассматривается как объект, составные элементы которого, различаясь режимными и конструктивными характеристиками, стыкуются между собой по определенному порядку гладкости. Эти модели рассмотрены в работах [6, 7].

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

Постановка задачи. Начальным этапом формирования любого участка автомобильной дороги является проектирование ее оси. Этот этап достаточно подробно рассмотрен в работе [7]. Будем считать ось дороги заданным геометрическим элементом, который представляет собой пространственную кривую не ниже порядка гладкости С 1 [7, 9].

Рис. 1. Устройство автомобильной дороги на виражных участках

Задачу моделирования проезжей части автомобильной дороги на отгонном участке можно сформулировать следующим образом: задана ось дороги в виде гладкой, не ниже порядка гладкости С 1, пространственной кривой P ( t ) , описываемой векторным уравнением:

P ( t ) = ( x ( t ), y ( t ), z ( t )); P '( t ) * 0; t e R : To t T . (1)

Требуется построить теоретическую поверхность отгонного участка автомобильной дороги.

Теория. Сущность циклографического метода заключается во взаимно однозначном отображении точек пространства R 3 на циклы в координатной плоскости ( xy ) декартовой системы координат Oxyz . Циклографический метод был предложен и развит немецкими геометрами W. Fiedler [10], E. üller и J. Krames [11]. Современный уровень развития циклографического отображения и его многочисленные теоретические применения и практические реализации изложены в работах [12–15].

Циклографическая проекция точки ( x, y, z ) пространства R3 представляет собой направленную окружность - цикл на плоскости (ху) с центром в точке (х, у) и радиусом R = |z|. Двум направлениям цикла соответствуют аппликаты z > 0 и z < 0. Следовательно, множество да3 точек пространства R3 можно поставить во взаимно однозначное соответствие с множеством да3 циклов плоскости (ху).

Рассмотрим циклографический образ оси дороги, представляющий собой пространственную кривую. Его получение основано на построении однопараметрического множества конусов, вершины которых принадлежат данной кривой, а основания (циклы) расположены в плоскости проекций (xy) и образуют однопараметрическое множество циклов [8]. Огибающая множества циклов будет являться циклографической проекцией пространственной оси дороги. Уравнения огибающей для изменяющихся значений полуугла β = β(t) при вершине конуса отображения имеют вид [8]:

x p (1,2) ( t ) = x ( t ) + z ( t ) e ( t )

- x ( t )W)+ y ( t )7^^.

x ( t ) ______________;

- x '( t ) • ц ( t ) ± x '( t ) 7 X ( t ) 2( t )

У р (1,2) ( t ) = У (t ) + z ( t ) e ( t )

X ( t )

,

где Ц ( t ) = e ( t ) z '( t ) + e '( t ) z ( t );

X ( t ) = x '( t )2 + y '( t )2; e ( t ) = tg ( P ).

Параметрические уравнения (2) описывают огибающую, которая может быть как единым объектом, так и состоять из ветвей P 3 (1) ( t ) и P 3 (2) ( t ). Центры всех циклов образуют линию P 1 ( t ) , представляющую собой горизонтальную проекцию исходной линии – оси дороги P ( t ) . Если принять (1) и (2) в качестве уравнений направляющих линий, то можно получить линейчатую поверхность, которая будет огибающей указанного однопараметрического множества конусов отображения с вершинами на пространственной кривой P ( t ) . Параметрические уравнения этой линейчатой поверхности имеют следующий вид:

X ( t , l ) = x ( t ) + l • [ x P (i, 2) ( t ) - x ( t )];

Y ( t , l ) = У ( t ) + l • [ У p (i,2) ( t ) - У ( t )];                (3)

Z ( t , l ) = z ( t ) • (1- 1 );

t , 1 e R : T 0 t T , L 0 t L .

При рассечении поверхности (3) нормальными по отношению к линии P1 (t) плоскостями будут образовываться кривые линии с вершинами на исходной линии P(t) . Полученные линии сечений отличаются от общепринятых, соответствующих известным правилам и нормам [2]. Поэтому возникает необходимость преобразования получен- ной поверхности в другую линейчатую поверхность, сечения которой указанными нормальными плоскостями соответствовали бы общепринятым типовым нормам. Предлагаемое преобразование основано на «доворачивании» общей с линейчатой поверхностью (3) пары прямолинейных образующих каждого конуса отображения в нормальную плоскость линии P1 (t) . Это преобразование подробно рассмотрено в работе [8].

Общий принцип геометрического моделирования поверхностей дорожного полотна заключается в следующем. На плоскости проекций (ху) строится горизонтальная проекция оси дороги – линия P 1 ( t ) . Затем по обеим сторонам этой линии строятся соответствующие ей линии – эквидистанты, удаленные от P 1 ( t ) на расстояние, равное сумме полуширины проезжей части и обочины. После этого на полученных линиях как на направляющих строятся вертикальные цилиндрические поверхности. Линия пересечения такой цилиндрической поверхности и линейчатой поверхности проезжей части, построение которой основано на циклографической проекции, будет являться бровкой проезжей части. Циклографическая проекция бровки проезжей части позволяет получить новую линейчатую поверхность, которая при пересечении с вертикальной цилиндрической поверхностью обочины даёт бровку обочины дорожного полотна. Таким образом, создавая поверхность между осью дороги и бровкой проезжей части, получаем поверхность проезжей части, а между бровкой проезжей части и обочины – поверхность обочины.

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

– считаем известными граничные условия стыкуемых к рассматриваемому отгону виража и прямолинейного в плане участков дороги, а именно: координаты стыкуемых крайних точек на бровках дорожного полотна и обочин, а также значения первых и вторых производных их уравнений в стыкуемых точках;

– уклоны образующих прямых линий поверхностей проезжей части и обочины с целью наглядности заведомо увеличены относительно требуемых норм и правил [2];

– при моделировании отгона не рассматривается уширение проезжей части на подходе к виражу.

В качестве оси дороги рассмотрим кривую P ( t ) = { x ( t ), y ( t ), z ( t ) } с параметрическими уравнениями:

x ( t ) = 12 t 2 + 36 t + 2, y ( t ) = 20 t^ + 9 t 2 + 3 t + 3, z ( t ) = 10 + 3 t , t e [0,1].

Известно, что отгон виражного участка представляет собой элемент автомобильной дороги, на котором происходит переход от двускатного профиля к односкатному и наоборот [2, 3]. Автомобильная дорога, в свою очередь, состоит из поверхностей проезжей части и прилегающей к ней обочины. Условно отгон виражного участка можно разбить на несколько участков, схематично представленных на рис. 2. Деление на участки накладывает некоторые конструктивные ограничения, которые необходимо учитывать при проектировании отгона.

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

Таким образом, формообразование линейчатых поверхностей на участках 1–4 основано

Рис. 2. Схема изменения профилей нормальных сечений поверхностной формы отгона нейной образующей поверхности на соответствующем участке. Отрезок прямолинейной образующей, как было отмечено выше, представляет собой профиль нормального сечения поверхности отгона. Изменение положения профиля происходит на основе его вращения в плоскости нормального сечения относительно точки пересечения нормальной плоскости и оси дороги.

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

Рассмотрим более подробно принцип формообразования поверхности на участках с переменным полууглом при вершине конуса отображения γ, γ = γ(t) (рис. 3). Фактически задача заключается в получении такой поверхности, которая бы не приводила к излому поверхности дороги в начале и в конце участка. В нашем случае к данному участку отгона прилегает с одной стороны прямолинейный участок дороги с двухскатным профилем, где угол наклона прямолинейных образующих поверхности обочины равен γ1, а с другой идет продолжение отгона, где у внешней поверхности обочины угол изменяется до угла γ2. Следовательно, можно получить две огибающие PY1(1) и PY2(1) циклов - оснований конусов циклографического отображения с полууглами при вершинах γ1 и γ2 соответственно. Решение задачи переносится на плоскость проекций (ху) и сводится к нахождению некоторой кривой, соединяющей начало огибающей P (1) с концом огибающей PY2 (1). Впоследствии данная кривая будет определять закон изменения угла наклона образующих поверхности обочины. Для получения такой кривой используется кубический сплайн, который позволяет добиться гладкости C1 в точках стыковки с кривыми соседних участков, что, в свою очередь, позволит добиться плавности стыковки получаемых впоследствии поверхностей обочины. Уравнение сплайна имеет вид [16]:

P ( t ) а о + ах t + a 2 t + a 3 t ,

где 1 1 t < 1 2 , t e R .

Рис. 3. К образованию линейчатой поверхности на участке с переменным полууглом при вершине конуса отображения

В рассматриваемом примере при углах γ1 = 80° и γ2 = 85° параметрические уравнения кривой P Y ( 1 ) имеют вид:

x ( 1 ) 3,5 + 42,6 1 + 54,4 1 2 - 67,2 1 3;

y ( t ) — - 14,9 + 3,5 1 + 16,9 1 2 + 16,8 1 3;

z ( 1 ) 8,1 + 3 1 + 20,9 1 2 - 46,7 1 3, где 0< t <0,3.

Уравнения полученной кривой P Y ( 1 ) совместно с уравнениями оси дороги P ( t ) определяют уравнения линейчатой поверхности, которая и будет искомой поверхностью обочины с переменным углом y e [ y 1 , Y 2 ] наклона ее образующих. Дальнейшая обрезка данной поверхности вертикальной цилиндрической поверхностью внешней обочины позволяет получить линейчатую поверхность обочины на данном участке отгона (рис. 4).

Рис. 4. Участок отгона с поверхностью внешней обочины, имеющей переменный угол γ наклона образующих

Параметрические уравнения этой поверхности имеют вид:

X ( t , l ) = 3,2 + 42,6 t + 32,3 t 2 - 18 t 3 + (0,2 + 22,1 t 2 - 49,2 t 3) l , Y ( t , l ) = - 11,9 + 3,5 1 + 12,5 t 2 + 26,6 t 3 + ( - 3 + 4,4 t 2 - 9,8 t 3) l , Z ( t , l ) = 8,9 + 3 t + 6,640"9 t 2 - 2,2Ю 8 t 3 +

+ ( - 0,8 + 20,9 t 2 - 46,7 t 3) l , где 0 <  t < 0,3; 0 <  l < 1; t , l e R .

Оставшаяся внутренняя часть поверхности отгона, вплоть до четвертого участка, остается неизменной.

На четвертом участке при неизменном уклоне прямолинейных образующих поверхности обочины образующие прямые проезжей части начинают менять угол наклона до тех пор, пока не произойдет совпадение с углом наклона образующих внутренней поверхности обочины (см. рис. 2, участок 4). Это происходит за счет вращения образующих внутренней проезжей части вокруг ее бровки. Таким образом, полный внутренний поперечный профиль отгона преобразуется в одну прямую линию. Здесь применяется тот же метод, что и для создания поверхности обочины на первом участке при изменении ее угла наклона, схема которого представлена на рис. 3. Следует отметить, что при вращении образующих прямых вокруг бровки проезжей части происходит изменение оси дороги. Новая ось дороги может быть найдена следующим образом. Из ортогональной проекции оси дороги P 1 ( t ) строится вертикальная цилиндрическая поверхность. Затем находится линия пересечения этой цилиндрической поверхности и линейчатой поверхности проезжей части с переменным углом наклона образующих. Полученная линия и будет «новой» осью дороги на данном участке (рис. 5).

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

Выше было рассмотрено формообразование поверхности первого участка внешней части отгона (рис. 4). Для остальных участков внешней части отгона построим циклографическую проекцию оси дороги, имеющую неизменный полуугол при вершине конуса отображения до конца отгона. На полученной огибающей построим вертикальную цилиндрическую поверхность Ψ и линейчатую поверхность Ω , как показано на рис. 6. На цилиндрической поверхности необходимо определить несколько ключевых точек, по которым можно построить сплайны, определяющие подъем образующих линейчатой поверхности на конкретных участках отгона.

Каждая из ключевых точек лежит на границах между соответствующими участками отгона. Первая из ключевых точек – точка 1 – находится в плоскости проекций (xy) на первой образующей полученной линейчатой поверхности Ω . Для построения кубического сплайна между точками 1 и 2 имеем следующие условия: координаты точек T 1 и T 2 , значения первой T 1 ' и второй T 1 '' производных в точке 1, полученные из найденной ранее огибающей линии (циклографической проекции) с

Рис. 5. Подъем оси дороги в конце отгона виражного участка

Рис. 6. Схема образования поверхности проезжей части внешней половины отгона виражного участка

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

X s i ( 1 ) = 1,4 t 3 4,7 t 2 + 66,2 1 + 60,7;

yS 1 ( 1 ) = — 2,8 1 3 + 13,8 1 2 + 6,3 1 141,2;

zS 1( t) = —1013 + 21,812, где 0 < 1 < 1.

Для нахождения точки 2 проведем горизонтальную прямую из последней точки второго участка до пересечения с цилиндрической поверхностью Ψ так, чтобы ортогональная проекция этой прямой была перпендикулярна ортогональной проекции оси дороги P 1 ( t ) . Полученная точка и будет искомой точкой 2. Аналогичным образом, зная координаты точек T 2 и T 3 , а также производные T 2 ' и T 2 '' , строим кривую s 2 , определяющую изменение угла наклона образующих для данного участка. Уравнения кривой s 2 следующие:

xS 2 ( 1 ) = — 1,2 1 3 2,8 1 2 + 61 1 + 123,6;

yS 2 ( 1 ) = — 2,6 1 3 + 12,4 1 2 + 25,6 1 123,8;

Z s 2 ( t ) = — 0,13 t 3 + 2,7 t 2 + 13,6 1 11,8, где 0 1 1.

Точка 3 находится на пересечении последней образующей третьего участка внутренней поверх- ности проезжей части и цилиндрической поверхности Ψ. Для определения кривой на этом участке имеем шесть граничных условий: T3 , T4 , T3 ' , T4 ' , T3 '' и T4 '' (для точки 4 в ходе принятых допущений принимаем известными параметры виражного участка, стыкуемого с рассматриваемым участком отгона). Для достижения второго порядка гладкости между точками 3 и 4 строим эрмитов сплайн пятой степени s3 , уравнение которого имеет вид [17]:

P ( t ) = a 0 + axt + a 2 t + a 3 t + a 4 t + a 5 t , где 1 1 1 < 1 2 , 1 e R .

Полученные в ходе вычислений параметрические уравнения кривой s 3 имеют следующий вид:

x S 3 ( t ) = 48,9 1 5 136,5 1 4 + 105,4 1 3 + 4,8 t 2 + 25 1 + 180,6; y S 3 ( t ) = — 10,2 1 5 + 40,6 1 4 52,2 1 3 + 33,8 1 2 + 33,8 1 88,4; zS 3 ( t ) = 41,5 1 5 128,5 1 4 + 114,1 1 3 6,1 1 2 + 1,8 1 + 25,4, где 0 1 1.

Таким образом, полученные кривые между точками 1–4, состыкованные по второму порядку гладкости и рассмотренные совместно с соответствующими им участками оси дороги на отгоне виража, позволяют получить линейчатые поверхности. Отсекая данные поверхности вертикальной цилиндрической поверхностью E ( t,l ) внешней кромки проезжей части, как показано на рис. 7, получаем искомые поверхности проезжей части на внешней стороне дорожного полотна. Таким же образом получаем поверхности для внешней обочины.

Итоговая визуализация отгона виражного участка в рассматриваемом примере представлена на рис. 8. Для большей наглядности показаны насыпи, полученные путем построения циклографической проекции от бровок обочины.

Заключение

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

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

Рис. 7. Схема образования поверхности проезжей части внешней половины отгона виражного участка

Рис. 8. Итоговая модель поверхностных форм отгона виражного участка

Список литературы Геометрическое моделирование поверхностных форм отгона виража автомобильной дороги на основе циклографического отображения

  • Кузьмин, В.И. Автоматизированное конструирование виражей безопасных конструкций на закруглениях автомобильных дорог/В.И. Кузьмин, А.И. Левтеров//Вестн. Харьк. нац. автомоб.-дор. ун-та. -2009. -Вып. 47. -С. 29-33.
  • Автомобильные дороги: СНиП 2.05.02-85. -Введ. 1987-01-01. -Госстрой СССР. -М.: ЦИТП Госстроя СССР, 1986. -56 с.
  • Бойков, В.Н. Автоматизированное проектирование автомобильных дорог (на примере IndorCAD/Road)/В.Н. Бойков, Г.А. Федотов, В.И. Пуркин. -М.: Изд-во МАДИ, 2005. -224 с.
  • Пуркин, В.И. Основы автоматизированного проектирования автомобильных дорог/В.И. Пуркин. -М.:МАДИ, 2000. -141 с.
  • Сальков, Н.А. Математическое моделирование линейных и поверхностных форм автомобильных дорог на подходах к мостам/Н.А. Сальков//Труды МАДИ: Прикладные теоретические вопросы проектирования переходов через водотоки. -М., 1989. -С. 60-66.
  • Система автоматизированного проектирования автомобильных дорог IndorCAD/Road 5.0. Справочное руководство. -Томск: ООО «ИндорСофт», 2004. -275 с.
  • Сальков, Н.А. Моделирование автомобильных дорог /Н.А. Сальков. -М.: ИНФРА-М, 2012. -120 с. -(Научная мысль).- 978-5-16-006756-8.
  • ISBN: 9785160067568
  • Panchuk, K.L. Cyclographic Modeling of Surface Forms of Highways/K.L. Panchuk, A.S. Niteyskiy, E.V. Lyubchinov//IOP Conf. Series: Materials Science and Engineering. -2017. -262. -012108.
  • DOI: 10.1088/1757-899X/262/1/012108
  • Федотов, В.А. Теоретические, экспериментальные и практические основы автоматизированного проектирования развязок движения на пересечениях автомобильных дорог в разных уровнях: дис. … канд. техн. наук/В.А. Федотов. -М., 1979. -292 с.
  • Fiendler, W. Cyklographie oder Construction der Aufgabenüber Kreise und Kugeln und elementare Geometrie der Kreis-und Kugelsysteme/W. Fiendler. -Leipzig, Druckund Verlagvon B.G. Teubner, 1882. -284 p.
  • Emil MullerDr. Vorlesungenüber Darstellende Geometrie. II. Band: Die Zyklographie. Edited from the manuscript by Dr. Josef Leopold Krames/Dr. Emil Muller. -Leipzig and Vienna, Franz Deuticke, 1929. -476 p.
  • Pottmann, H. Computational Line Geometry/H. Pottmann, J. Wallner. -Berlin. Heidelberg: Springer Verlag, 2001. -565 p.
  • Pottmann, H. Applications of Laguerre geometry in CAGD/H. Pottmann, M. Peternell//Comput. Aided Geom. Design. -1998. -№ 15. -P. 165-186.
  • Cho, H.C. Clifford algebra, Lorentzian geometry, and rational parametrization of canal surfaces/H.C. Cho, H.I. Choi, S-H. Kwon, D.S. Lee, N-S. Wee//Computer Aided Geometric Design. -2004. -Vol. 21. -P. 327-339.
  • Панчук, К.Л. Циклографическая начертательная геометрия: моногр./К.Л. Панчук, Н.В. Кайгородцева. -Омск: Изд-во ОмГТУ, 2017. -232 с.
  • Rogers, D.F. Mathematical Elements for computer graphics/D.F. Rogers. -McGraw-Hill Publishing Company -New York, 1990. -604 p.
  • Завьялов, Ю.С. Методы сплайн-функций/Ю.С. Завьялов, Б.И. Квасов, В.Л. Мирошниченко. -М.: Наука, 1980. -355 с.
Еще
Статья научная