CC BY
25
Наука и Образование
МГТУ им. Н.Э. Баумана
Наука и Образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2016. № 01. С. 38-51.
Б01: 10.7463/0116.0831018
Представлена в редакцию: 14.12.2015 Исправлена: 28.12.2015
© МГТУ им. Н.Э. Баумана
УДК 629.78
Радиальное равновесие протяженного объекта на круговой орбите
Коровин В. В.1' ' когоущлтШтыидц
:МГТУ им. Н.Э. Баумана, Москва, Россия
Анализ орбитального движения и равновесия протяженных объектов относится к числу актуальных проблем механики космического полета. В статье рассмотрен переход от материальной точки к объектам в виде гантели, однородного стержня, затем - типа гирлянды. Сформулированы уравнения равновесия для случая круговой орбиты. Показано, что движение центра масс подобных объектов является некеплеровым и они испытывают растяжение в радиальном направлении. Выведены аналитические выражения для силы натяжения и проанализированы свойства. Предложено рассматривать действующую на протяженный объект гравитационно-центробежную нагрузку как самоуравновешенную систему сил кажущегося веса с переменным вдоль радиуса кажущимся ускорением. Полученные результаты использованы для анализа космического лифта. Рассчитана минимальная длина лифта в виде однородного стержня. Определено распределение силы натяжения по длине и требования к прочности материала. Показан характер и величина действующего в сечениях протяженных объектов ускорения микрогравитации.
Ключевые слова: космический лифт, космическая тросовая система, равновесие на круговой орбите, орбитальная скорость, сила натяжения, орбитальный центр, ускорение микрогравитации
1.Введение
Задачи механики орбитального движения и равновесия протяженных объектов привлекают внимание отечественных и зарубежных специалистов на протяжении последних десятилетий. Это определяется интересом к таким проектам, как космические тросовые связки (системы) и космические лифты.
Космический лифт - конструкция, находящаяся в равновесии относительно вращающейся Земли или другого небесного тела [1, 2, 3, 4]. Космическая тросовая связка может находиться в относительном равновесии (гравитационно-уравновешенное положение) на круговой орбите [4, 5, 6, 7].
При исследовании космического лифта вопросы статики являются первоочередными. Особенностью лифта является известная угловая скорость вращения. Она определяется планетой, с которой связан лифт.
Исследование механики орбитальных тросовых связок во многих случаях не подразумевает анализа статического равновесия, т.к. решаются задачи динамики. Если рассматривается движение относительно центра масс, часто используется предположение о кеп-леровом движении центра масс протяженного объекта [8, 9, 10, 11]. Такое допущение справедливо в линейном приближении. Действительное гравитационное поле в этом случае заменяется его линейным аналогом в окрестности центра масс системы, что не приводит к заметным погрешностям, пока размеры объекта существенно меньше радиуса орбиты центра масс [5, 6]. Свободным от таких ограничений является подход, когда уравнения движения записываются в геоцентрической инерциальной системе отсчета [12].
В ряде работ допускается перенос результатов анализа линеаризованных моделей на объекты, немалые в сравнении с орбитальным радиусом (космический лифт, протяженные связки) [3]. В результате, например, центр масс космического лифта оказывается на геостационарной орбите, что, строго говоря, неверно.
В целом, вопросы статики протяженных орбитальных объектов многократно рассматривались отечественными и зарубежными авторами. При этом обычно анализ выполнялся в рамках решения конкретной прикладной задачи. В настоящей работе рассматриваются нелинейные уравнения равновесия одномерных объектов в центральном поле тяготения на круговой орбите. Ставится цель определения условий статического равновесия, исследования характера нагружения, получения удобных и физически наглядных расчетных соотношений.
2.. Материальная точка и гантель
На материальную точку в центральном поле тяготения действует гравитационная
сила
_ -^ГР - - — Г '
где
[I - гравитационный параметр притягивающего центра; М - масса материальной точки; Г - радиус-вектор.
Уравнения движения материальной точки в проекциях на оси инерциальной системы
отсчета, связанной с центром притяжения, имеют вид:
" - ^
Х1 ~ г3 '
где - координаты,
Г = / X2 + у2 + 22.
Введем связанную орбитальную систему отсчета . Начало отсчета совпада-
ет с материальной точкой, ось - местная вертикаль; ось
Су - трансверсаль, направлена по орбитальному движению; ось Сг - нормаль к плоскости орбиты и дополняет систему до правой.
В этой неинерциальной системе отсчета материальная точка находится в состоянии равновесия
Рг р+Ре = 0, (1)
где - переносная сила инерции.
Ограничимся случаем круговой орбиты. Тогда
Ре = М П 2 Г = Рц ,
где Рц - центробежная сила; ПП - угловая скорость орбитального движения.
Из условия равновесия (1) следует
п = $ ■ (2)
Здесь П - кеплерова угловая скорость для орбиты радиуса г, т.е. угловая скорость материальной точки в центральном поле тяготения на круговой орбите при отсутствии возмущающих сил.
Рассмотрим равновесие на круговой орбите простейшего тела, не являющегося материальной точкой. Это гантелеобразный объект: две материальные точки и М2, соединенные идеальной невесомой нерастяжимой связью. Пусть «гантель» ориентирована вдоль местной вертикали. Орбитальная система отсчета С ху г в этом случае связана с центром масс гантели (рис. 1).
Рис.1 Равновесие гантелеобразного объекта
Для равновесия необходимо равенство центробежных и гравитационных сил:
м^п - ^) = о ¿=1 4 1 ' Отсюда определяется угловая скорость орбитального движения П:
(3)
П2 = Ц Г М±Г^+М2Г£ 1
Сравним полученный результат с кеплеровой угловой скоростью (2). Положим Мх = М2 = М, расстояние между материальными точками равно 2 /. Радиус орбиты центра
масс гантелеобразного тела равен й , гх = й — / , г2 = й + /. Тогда:
п2 = [ ¡(¡^ ■ (4)
При получим , что соответствует кеплеровой угловой скорости. При
/ 0 из (4) следует: П > П ке п л ер ( й ) ■ Таким образом, движение центра масс гантелеобраз-ного объекта некеплерово.
Определим реакцию связи, действующей на материальные точки. Уравнение равновесия верхней массы (рис. 1):
М (п2Г2— I) — г = 0 , (5)
где - реакция связи.
Подставляя в (5) выражение для П 2 (4), находим
г=! М -^р- ■ (6)
Из (6) видно, что при положительных и будет положительное натяжение . При . Следовательно, гантелеобразный объект в радиальном положении на круговой орбите растянут.
3. Однородный стержень
Рассмотрим равновесие однородного стержня погонной массой и длиной , расположенного радиально на круговой орбите. Радиус орбиты центра масс стержня - й. Тогда его нижний и верхний концы движутся по орбитам радиусами и соответственно (рис. 2).
Рис.2 Радиальное равновесие однородного стержня
Условие равновесия стержня имеет интегральный вид
¡™(тП 2г-^) аг = 0. (7)
На стержень действуют центробежная сила
Я+1
Fц = тП2 I гйг = 2 т П2 Я I ;
Я-1
и гравитационная сила
Я+1
[ с1г 21
Кп = —11т — = — и т ——— .
гр И- ] г2 И- д2 _ ¿2
Я-1
Приравнивая силы, находим выражение для квадрата угловой скорости стержня
п2=1о^ ■ (8)
При I — 0 имеем: П2 = ц/ й3 - кеплерова угловая скорость.
Равновесие элементарного участка стержня йг определяется уравнением [7]:
£ = т(£ — П2г) , (9)
где - продольная сила в стержне (сила натяжения в тросе, если под стержнем понимается трос).
Интегрируя (9), получим
Т(г) = — т 0 + 1 П 2г2) + с . (10)
Подставляя в (10) выражение для угловой скорости (8) и граничное условие Т ( й — =0, найдем
3 й2 + /2
с = \хт
2Й (Я2 - 12) '
Т(г) = — |1 771
1 г^-г^ зк
г 2к (к2-г2) 2 (к2-г2). Максимальное натяжение Тт а х определяется условием
(11)
— = 0 . (12)
йг у '
Из (12) получим значение г, соответствующее Тт ах
г ( Ттах ) = [й (И2 —/ 2) ] 1 /3 . (13)
Анализ выражения (13) дает следующее.
При , или сечение с максимальным натяжением всегда ниже цен-
тра масс (середины однородного стержня).
При , т.е. сечение с максимальным натяжением приближается
к центру масс протяженного в радиальном направлении объекта снизу, если размеры объекта малы в сравнении с радиусом его орбиты.
При I ~ й (случай космического лифта) г ( Тт ах ) — 0. Это означает, что максимальное натяжение будет ближе к нижнему концу космического лифта постоянной погонной массы.
Условие радиального равновесия стержня (7) показывает, что гравитационная и центробежная силы, действующие на стержень, интегрально равны. При этом с ростом г центробежная сила растет, гравитационная - падает. Определим радиус сечения , в котором /гр и уравновешены и, следовательно, имеет место невесомость. Условие равновесия элемента стержня в этом случае:
йт П2г =
-у 2
Подставляя сюда значение угловой скорости (8), получим
г = [й (й2 — /2 ) ] 1/3 , (14)
что совпадает с формулой (13) - радиусом максимального натяжения г ( Тт ах) .
Сечение протяженного в радиальном направлении орбитального объекта, в котором существует невесомость и возникает наибольшая радиальная растягивающая сила, называют гравитационным центром или орбитальным центром (ОЦ) [5].
Из сравнения формулы для квадрата орбитальной скорости (8) с радиусом орбитального центра (13), (14) ясно, что движение ОЦ происходит с кеплеровой угловой скоростью (2):
П 2 =_£_= Л- (15)
П * Я (Я 2 г2 ) дзц . (15)
Таким образом, если задан радиус орбиты центра масс протяженного объекта, его орбитальная скорость определяется из условия равенства центробежных и гравитационных сил (7). Определяя, затем, радиус орбиты, для которой найденная угловая скорость является кеплеровой, находим положение ОЦ протяженного орбитального объекта (ПОО).
4. Общий подход к решению задачи и числовой пример
При смещении в радиальном направлении от ОЦ равенство гравитационной и центробежной сил нарушается (рис. 2). Выше ОЦ центробежная сила превышает гравитационную, ниже - наоборот. Следовательно, точки ПОО, радиус которых отличается от радиуса ОЦ, находятся в условиях микрогравитации. Величина ускорения микрогравитации а м к гр (г) определяется разностью ускорений от центробежной и гравитационной сил. Она
растет нелинейно с удалением от орбитального центра:
гз_^ з
а М КГр (г) = Д г2р з . (16)
т коц
Ускорение микрогравитации направлено от ОЦ по радиус-вектору в зенит при
и к центру притяжения (в надир) при .
Ускорение (16) можно трактовать, как кажущееся ускорение, а определяемые им силы - как силы кажущегося веса. Тогда протяженный объект в радиальном равновесии на круговой орбите растянут самоуравновешенной системой сил кажущегося веса, как показано на рис. 2.
Космическая тросовая система в общем случае имеет вид «гирлянды» из сосредоточенных масс и соединительных тросов. Радиальное положение на круговой орбите в центральном поле тяготения является одной из возможных форм равновесия подобных объектов [6]. Эта форма равновесия является устойчивой, так как при отклонении от нее силы кажущегося веса создают восстанавливающий момент. Множество способов применения тросов в космосе предполагают создание именно радиально-уравновешенных систем.
Для статического расчета такого объекта необходимо определить орбитальную угловую скорость связки по известному радиусу какой-либо точки, или наоборот, по заданной скорости найти радиус, обеспечивающий орбитальное равновесие. В случае космиче-
ского лифта угловая скорость задается небесным телом, с которым связан лифт. Тогда требуется определить минимальную длину лифта, обеспечивающую радиальное равновесие. Далее, производится расчет одномерной конструкции на действие самоуравновешенной нагрузки в виде сил кажущегося веса от ускорения (16). При этом, сила натяжения в сечении r для ПОО с орбитальными радиусами нижнего и верхнего концов R! и R2 соответственно, определяется выражением:
Д2
Т(г) = j т(ОаМкгр(0^ + ^М;амкгр(г;) =
Г i
^мкгр (f)df+ ZуМуОмкгр (Г}) ) , (17)
где Mi - сосредоточенные массы на участке гирлянды от r до R2 ;
Mj - сосредоточенные массы на участке от Rj до r .
С помощью полученных здесь аналитических выражений проанализируем, кратко, свойства простейшего земного космического лифта. Такой лифт представляет собой однородный стержень, закрепленный на экваторе Земли и простирающийся в радиальном направлении. Минимальная длина лифта определяется из условия отсутствия сжимающих напряжений. Запишем условие равновесия (7) в виде
J^2 (mil2 r - dr = 0 , (18)
где R 2 - радиус верхнего конца, определяемый из условия (18).
Раскрывая (18), получим, при m = с о ns t, уравнение относительно R 2 в виде:
l 2R ! ■ R l — ( 2 f + l 2R l) ■ R 2 + 2 fR ! = 0 (19)
Примем:
R х = 63 78 к м - экваториальный радиус Земли (радиус нижнего конца лифта);
f = 3,986 ■ 1 0 5 — - гравитационный параметр Земли;
П = 7,292 ■ 1 0 _ 5 — - угловая скорость вращения Земли.
При этих числовых значениях корень уравнения (19) R2 = 150 162 км. Два других корня соответствуют нулевой или отрицательной длине лифта. Следовательно, полная длина космического лифта 2 / = R 2 — R х=143784 км. Полученный результат соответствует данным других авторов [1, 2, 4].
Геоцентрический радиус центра масс лифта составит R = R ! + I = 782 70 км. Используя (13), (14) или (15) находим радиус орбитального центра R 0 ц = 42 164 км, что равно радиусу геостационарной орбиты. Таким образом, центр масс лифта постоянного сечения выше геостационарной орбиты на 36106 км.
Выражения (11), (17) позволяют вычислить силу натяжения в лифте Т (г) . Разделив Т ( г) на погонный вес т^ , получим результат в виде разрывной длины Ь Р ( г) (рис. 3). Эта величина показывает, какой разрывной длиной (показатель удельной прочности) должен обладать материал лифта в сечении г, чтобы выдержать действующую нагрузку. Максимум разрывной длины имеет место на геостационарной орбите и равен 4950 км.
Требования к прочности материала космического лифта существенно снижаются при создании лифта переменного поперечного сечения [2, 4, 13]. В любом случае, центр масс лифта находится намного выше орбитального центра.
5000 4000
2000 1000 о
1? ?'."' 45 (?'."' 75 9'."' 105 120 1?5 15'."' геоцентрический радиус ¥._ км * Ю-3
Рис. 3. Нагружение космического лифта
На рис. 4 показано возникающее в сечениях лифта кажущееся или микрогравитационное ускорение.
Термин «микрогравитация» оправдан для объектов, размер которых существенно меньше радиуса орбиты. В случае космического лифта на нижнем конце выражение (16) дает а ( й х) = — 9 , 8 м/с 2 , что соответствует земному ускорению свободного падения. Знак
«минус» в данном случае означает направление действия против радиуса-вектора (к центру притяжения). Следовательно, кажущийся вес на нижнем конце космического лифта совпадает с весом на Земле.
Ускорение равно нулю (невесомость) в орбитальном центре лифта, т.е. на геостационарной орбите.
г
Рис. 4 Гравитационно-центробежное ускорение вдоль космического лифта
5. Заключение
В статье рассмотрено равновесие на круговой орбите объектов в виде гантели, однородного стержня, затем - типа гирлянды. Показаны присущие им свойства:
- движение центра масс некеплерово. Кеплерову угловую скорость имеет точка протяженного в радиальном направлении объекта, называемая орбитальным центром;
- протяженный объект на орбите растянут гравитационно-центробежными силами;
- для орбитального центра характерно наличие невесомости и максимума силы натяжения. Орбитальный центр всегда ниже центра масс.
Получены следующие результаты:
- выведены аналитические выражения для орбитальной угловой скорости центра масс рассмотренных протяженных объектов;
- получены аналитические выражения силы натяжения в типовых случаях и проанализированы их свойства;
- показан характер и величина действующего в сечениях протяженных радиально-ориентированных объектов ускорения. Оно определяется разностью гравитационных и центробежных сил;
- предложено рассматривать действующую на протяженный объект гравитационно-центробежную нагрузку как самоуравновешенную систему сил кажущегося веса с переменным вдоль радиуса кажущимся ускорением. Это дает общую формулу для определения силы натяжения в сечениях радиально-ориентированных одномерных объектов (17). Выражение (17) является более простым и физически наглядным в сравнении с исполь-
зуемыми в работах [2, 4], где натяжение определяется через потенциал гравитационно-центробежного поля.
Выполнен расчет простейшего космического лифта. Определены минимальная длина лифта в виде однородного стержня, распределение силы натяжения по длине и требования к прочности материала. Показано распределение кажущегося ускорения.
Корректность используемых аналитических выражений и подхода к решению задачи в целом подтверждена совпадением ряда полученных результатов с данными других авторов [1, 2, 4] или с физически очевидными величинами.
Список литературы
1. Pearson J. The Orbital tower: a spacecraft launcher using the Eath's rotational energy // Acta Astronautica. 1975. Vol. 2, no. 10. P. 785-799. DOI:
2. Поляков Г.Г. Собрание трудов. Т. 1. Привязные спутники, космические лифты и кольца. Астрахань: Изд-во Астрахан. гос. пед. ун-та, 1999. 580 с.
3. Bolonkin А. Non Rocket Space Launch and Flight. Elsevier, 2005. 488 p.
4. Нуралиева А.Б. О динамике троса космического лифта: дис. ... канд. физ.-мат. наук. М., ИПМ им. М.В. Келдыша РАН, 2012.
5. Arnold D.A. The behavior of long tethers in space // Journal of the Astronautical Sciences. 1987. Vol. 35, no. 1. P. 3-18. [русск. перев.: Арнольд Д.А. Поведение системы из двух связанных длинным тросом тел // Астронавтика и ракетодинамика: Экспресс-информация. 1988. № 2. С. 1-9.]
6. Белецкий В.В., Левин Е.М. Динамика космических тросовых систем. М.: Наука, 1990. 336 с.
7. Усюкин В.И. Строительная механика конструкций космической техники. М.: Машиностроение, 1988. 392 с.
8. Калашников Л.М., Малышев Г.В., Свотин А.П. Тросовые системы в орбитальных транспортных операциях // Полет. 2001. № 4. С. 51-56.
9. Сидоров И.М., Веселова Г.В. О возможности реализации пассивной одноосной гравитационной стабилизации международной космической станции. М.: ИКИ РАН, 2013. 18 с. (Препринт / ИКИ РАН; № 2169).
10. Иванов В.А., Купреев С.А., Ручинский В.С. Космические тросовые системы: учебное пособие / под ред. проф. В.А. Иванова. М.: Альфа-М; ИНФРА-М, 2014. 208 с.
11. Иванов В.А., Купреев С.А., Ручинский В.С. Орбитальное функционирование связанных космических объектов: учеб. пособие / под ред. проф. В.А. Иванова. М.: ИНФРА-М, 2014. 320 с.
12. Коровин В.В., Попов А.В., Усюкин В.И. Динамика неуправляемого развертывания космической тросовой связки // Инженерный журнал: наука и инновации. 2012. № 8. DOI: 10.18698/2308-6033-2012-8-444
13. Ледков А.С., Пикалев Р.С. Исследование влияния движения подъемника на динамику космического лифта // Наука и образование. МГТУ им. Н.Э. Баумана. Электронный журнал. 2014. № 5. С. 206-216. DOI: 10.7463/0514.0710704
Science ¿Education
of the Baumail MSTU
Science and Education of the Bauman MSTU, 2016, no. 01, pp. 38-51.
DOI: 10.7463/0116.0831018
Received: 14.12.2015
Revised: 28.12.2015
© Bauman Moscow State Technical Unversity
Radial Equilibrium of Extended Object in a Circular Orbit
V.V. Korovin1' 'korovinwigSbmstuju
:Bauman Moscow State Technical University, Moscow, Russia
Keywords: space elevator, space tethered system, equilibrium in a circular orbit, the orbital velocity,
tension force, the orbital center, microgravity acceleration
Over the past decades the problems in mechanics of orbital motion and equilibrium of extended objects draw attention of domestic and foreign experts because of their interest in such projects as the space tether systems and space elevators.
The article deals with the nonlinear equations of equilibrium one-dimensional objects in a central gravitational field in a circular orbit. It considers transition from the material point to the objects in the form of dumbbells, a homogeneous rod, and a type of garland. It is shown that the motion of the mass center of these objects is non-Keplerian. In some cases the analytical expressions of the angular velocity of the center of mass are obtained.
A point of the object extended in the radial direction, called the orbital center, has a Keplerian angular velocity. The orbital center is characterized by a zero gravity and a maximum of the tension force. The orbital center is always below the center of mass.
The gravitation-centrifugal forces provide stretching the extended object in orbit. The analytic expressions for the tension force in typical cases have been received, and their properties have been analyzed. The article shows the nature and acceleration value in the sections of the radially extended objects of acceleration that is determined by difference between gravitational and centrifugal forces.
The article offers to consider a gravitation-centrifugal load, which acts on the extended object, as the self-balanced system of forces of the apparent weight with the variable apparent acceleration along a radius. This gives a simple general formula for determining the tensile force in the cross sections of the radially oriented one-dimensional objects.
The calculation of the simplest space elevator has been fulfilled. Provided that there is a lack of compressive stresses, have been determined the minimum length of an elevator as a homogeneous rod, distribution of the tensile force by length, and requirements for material strength. The distribution of the apparent acceleration is shown.
The coincidence of a number of the results with those of other authors or physically apparent values proves the correctness of used analytical expressions and approach to the problem, as a whole.
References
1. Pearson J. The Orbital tower: a spacecraft launcher using the Eath's rotational energy. Acta Astronautica, 1975, vol.2, no. 10, pp. 785-799. DOI:
2. Polyakov G.G. Sobranie trudov. T. 1. Privyaznye sputniki, kosmicheskie lifty i kol'tsa [Collected Works. Vol. 1. Tethered satellites, space elevators and rings]. Astrakhan, ASPU Publ., 1999. 580 p. (in Russian).
3. Bolonkin A. Non Rocket Space Launch and Flight. Elsevier, 2005. 488 p.
4. Nuralieva A.B. O dinamike trosa kosmicheskogo lifta. Kand. dis. [On the dynamics of a space elevator cable. Cand. dis.]. Moscow, Keldysh Institute of Applied Mathematics (Russian Academy of Sciences), 2012. (in Russian, unpublished).
5. Arnold D.A. The behavior of long tethers in space. Journal of the Astronautical Sciences, 1987, vol. 35, no. 1, pp. 3-18.
6. Beletskii V.V., Levin E.M. Dinamika kosmicheskikh trosovykh system [Dynamics of space tether systems]. Moscow, Nauka Publ., 1990. 336 p. (in Russian).
7. Usyukin V.I. Stroitel'naya mekhanika konstruktsii kosmicheskoi tekhniki [Structural mechanics of space application constructions]. Moscow, Mashinostroenie Publ., 1988. 392 p. (in Russian).
8. Kalashnikov L.M., Malyshev G.V., Svotin A.P. Tethered systems in orbital transport operations. Polet = Flight, 2001, no. 4, pp. 51-56. (in Russian).
9. Sidorov I.M., Veselova G.V. O vozmozhnosti realizatsii passivnoi odnoosnoi gravitatsionnoi stabilizatsii mezhdunarodnoi kosmicheskoi stantsii. Preprint no. 2169 [On the possibility of realization of a passive single-axis gravity stabilization of the international space station. Preprint no. 2169]. Moscow, Space Research Institute of Russian Academy of Sciences, 2013. 18 p. (in Russian, unpublished).
10. Ivanov V.A., Kupreev S.A., Ruchinskii V.S. Kosmicheskie trosovye sistemy [Space tethered systems]. Moscow, Al'fa-M Publ.; INFRA-M Publ., 2014. 208 p. (in Russian).
11. Ivanov V.A., Kupreev S.A., Ruchinskii V.S. Orbital'noe funktsionirovanie svyazannykh kosmicheskikh ob"ektov [The orbital operation of the associated space objects]. Moscow, INFRA-M Publ., 2014. 320 p. (in Russian).
12. Korvin V.V., Popov A.V., Usyukin V.I. Dynamics of Uncontrollable Deployment of a Tethered Satellite System. Inzhenernyy zhurnal: nauka i innovatsii = Engineering Journal: Science and Innovation, 2012, no. 8. DOI: 10.18698/2308-6033-2012-8-444 (in Russian).
13. Ledkov A.S., Pikalev R.S. Study of influence of climber motion on the space elevator dynamics. Nauka i obrazovanie MGTU im. N.E. Baumana = Science and Education of the BaumanMSTU, 2014, no. 5, pp. 206-216. DOI: 10.7463/0514.0710704 (in Russian).
