с1о1: 10.5862/МСЕ.55.7
Воздействие морских волн на судно, ошвартованное у причала с камерой гашения
К.т.н., доцент К.К. Семенов,
Санкт-Петербургский политехнический университет Петра Великого
старший преподаватель В.А. Леонтьев, д.т.н., профессор И.С. Нуднер,
Балтийский государственный технический университет "ВОЕНМЕХ" им. Д.Ф. Устинова
Аннотация. При проектировании причальных сооружений и систем раскрепления судов у причалов требуется обеспечить выполнение условий на возможные перемещения судна. Для обоснованного выбора конструктивного решения требуется оценить характеристики воздействия морских волн на судно с учетом возникающего волнового режима. Для инженерных изысканий удобно использовать результаты численных расчетов.
В работе представлена математическая модель для расчета воздействия волн на судно, ошвартованное у причала, содержащего камеру гашения. Для построения модели применялся аналитический метод, использующий разложение по собственным функциям.
Представлены результаты расчетов, демонстрирующие влияние характеристик причала на нагрузки, испытываемые судном.
Ключевые слова: причал; камера гашения; швартовка судна; воздействие морских волн
Введение
В настоящее время в морском гидротехническом строительстве все чаще применяются частично проницаемые конструкции берегозащитных и причальных сооружений [1]. Их волногасящие свойства обеспечиваются использованием камер гашения, лицевая часть которых представляет собой решетчатые конструкции из свайных рядов. Безопасная стоянка судов у таких причалов и безаварийное проведение погрузочно-разгрузочных работ достигаются надежной системой раскрепления, схема и параметры которой назначаются с учетом действующих на судно гидродинамических сил. Наличие причала может оказывать существенное влияние на значения действующих со стороны жидкости сил. Данная работа посвящена определению взаимосвязи гидродинамических сил (возмущающих и инерционно-демпфирующих) с геометрическими параметрами и физическими характеристиками причала, включающего камеру гашения.
Обзор литературы
Воздействие морских волн на судно, расположенное у причала, определяет выбор схемы и элементов систем раскрепления, обеспечивающих выполнение условий по допускаемым перемещениям судна при выполнении на нем погрузочно-разгрузочных работ. При проектировании возникает необходимость в определении характеристик гидродинамических воздействий на судно [2]. При этом следует учитывать конструкцию причала. Искомые значения гидродинамических характеристик, как правило, получают с помощью численного или физического моделирования.
Экспериментальное исследование воздействия морских волн на судно, ошвартованное у причала, достаточно трудоемко и сопряжено с определенными сложностями в обеспечении подобия между результатами, достигнутыми в рамках физической модели, и их натурными аналогами [3-5].
К настоящему моменту основными методами математического описания подобных задач являются [6]: аналитические методы разложения по собственным функциям или построения функций Грина, численно-аналитические методы систем дискретных источников и численные методы на разностных схемах. Перечисленные методы применяются при решении задач [7-11], близких к рассматриваемой в данной статье, и позволяют достичь приемлемой точности результатов расчета.
Известны работы, посвященные воздействию волн на одиночное судно [12, 13], на судно, расположенное у причала с непроницаемой вертикальной стенкой [14], у проницаемой преграды [15, 16], у сооружения откосного профиля [17, 18]. Публикации о нагрузках на судно,
расположенное у причала с камерой гашения, авторам неизвестны. В то же время есть практическая необходимость в рассмотрении таких вопросов. Настоящая работа посвящена решению задачи о воздействии волн на судно, расположенное у причала с камерой гашения. Представленные результаты могут быть использованы для получения инженерных рекомендаций при проектном рассмотрении систем раскрепления судов у причальных сооружений.
Постановка задачи
Взаимодействие волн с ошвартованным у причала судном исследуется при следующих предположениях: дно акватории считается горизонтальным, жидкость - идеальной и несжимаемой, ее движение - потенциальным и установившимся. Рассматривается плоская линейная задача. Судно схематизируется цилиндрическим телом прямоугольного сечения неограниченной протяженности. Причал представляет собой камеру гашения с вертикальными лицевой тонкой (пористой) и тыльной (непроницаемой) стенками.
Вводится декартова система координат Oxz, ось Ох которой совпадает с направлением распространения волн. Ось Оz проходит через центр тяжести судна вертикально вверх. Расположение судна у причала и обозначения геометрических размеров приведены на рисунке 1.
Рисунок 1. Схема расположения судна у причала
В силу линейности задача о воздействии морских волн на плавающее судно распадается на две [19]: задачу определения гидродинамических (возмущающих и инерционно-демпфирующих) сил, действующих на судно, и задачу определения параметров возникающей качки. Первая задача сводится к решению краевой задачи для уравнения Лапласа при соответствующих граничных условиях, вторая - к интегрированию системы обыкновенных дифференциальных уравнений, описывающих движение судна.
где А =
Набегающие из бесконечности волны будем описывать потенциалом скорости
М-сИ [о • €+</Зг,Чо0:~0'Ч
g■h■i
(1)
2 • о • сЬ Со
- комплексная амплитуда; И, а и X - высота, частота и длина волн
соответственно;
к0 =
2 • к
- волновое число; g- ускорение свободного падения; /' = V—Т -
мнимая единица.
Частота волн, волновое число дисперсионным соотношением
2
и местная глубина жидкости связаны между собой
О
= £•¿0 -Л Со
(2)
При взаимодействии набегающих волн с ошвартованным у причала судном часть волн отражается от преград, распространяясь в обратном направлении. Набегающие волны вызывают колебания судна, которое является дополнительным источником возникновения волн. Суммарное волновое поле будет описываться функцией [19, 20]:
ФС, z, tj= ф z, t~y Ф С*, z, 0+ Ф z, f,
(3)
где Ф г, / - дифракционный потенциал скорости, описывающий движение жидкости, вызванное присутствием судна как неподвижной преграды при набегании на него исходных волн;
z, / - потенциал скорости, описывающий волновое движение, обусловленное вынужденными колебаниями судна при отсутствии набегающих на судно волн. В этом случае исходную задачу можно представить в виде двух задач: дифракционной задачи для функции
Ф^С,г,/ и гидродинамической задачи для функции Ф^-^г.г,? , которые можно решать независимо одна от другой.
Функция Ф^, г, I удовлетворяет уравнению Лапласа
д2Ф д2Ф
дх2 dz
+ -
= 0
(4)
и следующим граничным условиям:
• кинематическому и динамическому на свободной поверхности, объединенным в одно условие
д2Ф
дг dz
непротекания через непроницаемые поверхности дна и причала
8Ф п
-= 0, z =-а , x<a+b+w\
дп
дФ
—- = 0, х = а + Ь + ч>, ге [■ 0 , дп
где п - нормаль к непроницаемой поверхности;
• непротекания через поверхность £ колеблющегося судна
(5)
(6) (7)
ЭФ ж. ^ а
(8)
где un4(,z,t - нормальная составляющая скорости судна:
ип С, z,t~j= [т • cos С, х~У- иz • cos 4}, z j)- со • 4i- cos С, z~y~ z - cos С, о • t ,
(9)
где ux , uz - проекции скорости судна на оси Ох и Оz; со - угловая скорость колебаний судна относительно оси О у, перпендикулярной плоскости Oxz; cos , cosi) - направляющие косинусы нормали n;
• на пористой стенке причала
5Ф_ 5Ф,
= ß■ф+-Ф.
(10)
дх дх
где Ф , Ф+ - значения потенциала скорости на лицевой и тыльной гранях пористой стенки; (3 -
коэффициент, зависящий от коэффициентов сопротивления f пористой стенки, пористости в и присоединенной массы Ст, а также от ее условной толщины Ь„ [21];
• на бесконечном удалении от судна (условие Зоммерфельда)
ф=A-
e
-i-k0-x
Ao+ Co]eik0x2chko- * + d^e
J^i. О-г e
x —> -00 ,
(11)
где А0 - амплитуда отраженной волны от неподвижного судна и причала; С0 - амплитуда волны, излучаемой колеблющимся судном при отсутствии набегающих волн.
С учетом равенства (3) потенциал скорости Фопределяется из соотношений (4)-(11) в предположении неподвижности судна (то есть в (9) нормальная составляющая скорости ип = 0 и судно не излучает волн, С0 = 0). Потенциал определяется из задачи (4)—(11) в
предположении, что отсутствуют набегающие и отраженные волны (А = А0 = 0).
Метод решения
Для решения задачи область О., занятая жидкостью, представляется в виде прямоугольных подобластей Оу (/'=1,2,3, 4) (рис. 2).
Рисунок 2. Разбиение области, занятой жидкостью
В каждой из областей Оу функции Ф*- и Ф^- определяются методом разделения переменных. Для функции 3= Ф^^.г,?^-Ф^^г.г,?^ч-Ф^^г,^ в областях Ц получим
следующие выражения:
ф, =
л -К,уХ
А-е " • сое
I -> --п -
Ь • « + О" X *«С+Ф-еК»-х ■ сов Ь < сС ■ е™ ,
и=0
(12)
ф2 =
о +
<Р ^ • п -ж - ^ + , п-ж-^ + с]
■ 81П ■
п=1
2-а
-сЬ-
2-а
+
I
п=1
1
_ 11-К-Х 11-К-Х
И И
I
И 11
соэ-
• 71 • С + ^
(Л-Т
„г-о/
(13)
ф, =
00 м-
е!
п п
-К„-х е п •сов
п= о
Ф4 =
и=о
где вп =
л-тИ2 . , лтг- Ц-Т^ 1 этИ ^
"7?
(14)
(15)
; К0-1-к0] значения Кп - положительные корни
2 а
уравнения G2-g^K■tg^■d , решенного относительно К, упорядоченные по возрастанию
(л = 1,2, ...), КХ<К2< ...; комплексные коэффициенты С^", (/=1.2,3,4;
п = 0, 1, 2, ...) - неизвестные.
2
Пары коэффициентов А^ и определяют потенциал Ф^^г.г,? , а С^- и -
потенциал . Соотношения (12)—(15) описывают функцию Ф, если положить
= 0 (/ = 1, 2, 3, 4; п = 0, 1, 2, ...) и их = и2 = ю = 0, и функцию Ф^-^.г,? при условии, что А$-= = 0 (/'= 1, 2, 3, 4, п = 0, 1, 2, ...) иЛ = 0.
Для определения неизвестных коэффициентов используются условия равенства давлений р/х, г, 0 и скоростей г, /)(/'= 1, 2, 3) на жидких границах подобластей и 02, Ог и Оз
Pi^a-O, z, / j= Р2 к а + 0, z, Vj Кг а - 0, z, / j= v2 k a + 0, z,
(16) (17)
условия непротекания жидкости через стенки контура судна и условие (10) на пористой стенке.
Используя свойство ортогональности собственных функций и условия (10), (16), (17), можно получить следующую систему алгебраических уравнений:
x^e~Kma . J -г xm e J m x0
n=1
x^- e d~T ed~T
'1nm
= r1n
A-eK°a +Ux-J 0+X-
n1
1^1 ^ ю- a 2-|+k
n-n . n-n• d-T
--sh-=
2 2-a
I'no~ юJ0
(18)
V P~Ko'a . T 4- V /7 - V V V .Г - - /1 . « . т
х0 е 100+х0 а У о +ZjX» е ''»О - Ае J00>
н=1
(19)
т-к-а
X' -А'„-а г _}_ , €
Хп 'е ' пт 2~Хт г)Ут
>1=О
d-T =-А.ек^-1
0m
(20)
и= 1
n п a \
x^- e d-T - y n2^- ed~T
+ ^ d om - rim
I x С e
nm m
Kma _ -, „ ~Km a
>
I. eKm a — У e Km J +
nm m m m
1*1 ^ ю-a2
. ^u7-a- |- (-1 , --ю-a
■J o"Z—
n1
n• n , n-к- d~T
--sh-=
2 2- a
I no ю •J o
(21)
¿Vi 'e + .Ун ' ^ nm _ '
n=1 2
xm^'r4m 'e
d-T j. .л d-T
Ут }Ът e
= 0,
/ >>
she, • 0,
(22)
(23)
(24)
где m = 0, 1, 2, ...,
Г1, да = 0
r1m
o, ^^ o
m
\ — 2- a, m = o
I 1, m* o ,
12, ^ = o
11, m * o ,
r4„, =
12 a, m = o
11, ^ o,
n-к- a
n-к- a
1
т-к-а
n n a
U
x
т-ж-а
m-Ti-a
\Kj-d
(
2
1 + -
вт (-К„ ■(!
п d _
2-К -(Л
■!'„ = {сое ¡Г„ • С + б/ >г =
вт
-г о
К„
Зп = Ь-сов С + =---+-----,
Кп'а Сп-а1
-т
= ——■ I СОБ-- — С ОБ Ъп-Ь + Л-
а -Т 1
а 2 • а
^т
Ъп-Щ-Т^-Ь.к.
4
— ■ т ■ л ■ бЬ
2
2-а
т-к- 4} -Т 2■ а
■ соб §ГИ • - а-сЬ
т-к- 4} -Т
2-а
п = 0, 1, 2, ...
Система уравнений (18)-(24) построена относительно неизвестных и 1, 2, 3, 4.
Система уравнений позволяет получить решения дифракционной задачи и гидродинамической задачи о колебаниях тела. Для определения неизвестных коэффициентов и В^"
необходимо в правых частях уравнений положить их = и2 = со = 0 и тогда Ар, ур= Вр.
Если же положить А = 0, то хр= Ср, Ур= В случае, когда задано ненулевое значение амплитуды А и хотя бы одно из значений их, и2, со не равно нулю, система определяет потенциал Ф: ур=Вр+В$-.
Система уравнений (18)-(24) является бесконечномерной. Для практического ее применения количество членов рядов ограничивалось сверху величиной Л/<с», то есть во всех соотношениях полагалось, что п и т принимали значения 0, 1, ... N. Таким образом, поиск решений полученной системы осуществляется ее редукцией к конечной размерности N и, соответственно, заменой функциональных рядов для потенциалов скоростей в различных расчетных областях их частичными суммами Л/-го порядка.
После определения неизвестных коэффициентов Ар, Вр, Ср и Ир М0ГУТ быть
вычислены основные гидродинамические характеристики воздействия волны на судно, ошвартованное у причала.
Величина гидродинамического давления в произвольный момент времени t может быть определена как
( > > дФ ^ ^ . *
р ^Ч, ! = -р--= -р-О-7 1
д!
\ = 1, 2
(25)
Определив давление, испытываемое поверхностью судна, можем вычислить значения действующей на судно силы [19]:
X<- |\2, СгР*^а, г, .
-т
7=1,2.
(26)
(27)
тп
d
—а
Значения и определяют нагрузку на судно, вызванную воздействием на него
морских волн. Величины и определяют гидродинамическую (инерционно-
демпфирующую) силу, действующую на судно при его качке.
Значения Х^-Х, и
связаны с присоединенными массами и коэффициентами сопротивления судна соотношениями [19]:
х «Ч>
Z
^ 22--* 22 Uxi
^ 33
dt dUz< dt
^33 -Uzi
(28) (29)
где £/,<>, V0', и,<Уи,.е™, ^22 и ц3з -
присоединенные массы качки судна; X 22 и X 33 - коэффициенты волнового сопротивления.
Из соотношений (28) и (29) легко вывести уравнения, разрешив которые, можно получить значения \х22, Л22, Н33 и Л33.
■ G ■ [122 + ^22
Z^=-Wz • 1-0-Ц22 +^22 _■
(30)
(31)
После определения возмущающих и инерционно-демпфирующих сил по уравнениям движения судна могут быть определены его предельные перемещения, необходимые для расчета усилий в швартовных связях.
Численные результаты
Представленная математическая модель была использована для расчетов, выполненных с целью оценки границ влияния причала на испытываемые судном нагрузки. Ниже приведены примеры расчетов зависимостей горизонтальной составляющей силы, воздействующей на судно, от различных параметров задачи.
Для выполненных расчетов были приняты следующие значения параметров (кроме варьируемого).
• геометрия задачи: а=15м;Ь = 5м;с( = 20м;7=4м;£) = 5м;и/=5м;
• параметры пористой стенки: Ь„= 1,0 м; 8 = 0,5; f= 2,0; Ст = 1,0;
• количество членов ряда: N = 25.
0,43 0,42
0,405 -0,400 -
0,370 -
....................I.....................J.....................L....................
х X • результаты расчета -аппроксимация
р g h a
-
N
10
15
20
0,39 0,38 0,37 0,36
*
/•
Z ■■-' р g h a • результаты расчета -аппроксимация
* N
0
10
15
20
Рисунок 3. Зависимость значения составляющих X и Т возмущающей силы
от числа N слагаемых в рядах
Рисунок 3 демонстрирует уточнение результата расчета с ростом числа членов рядов в представленных выше соотношениях для горизонтальной и вертикальной составляющих силового воздействия морских волн на судно. Видим, что уже при сравнительно небольших значениях N результат расчета перестает значимо изменяться. Приемлемыми являются значения Л/> 5. Таким образом, взятое значение N = 25 обеспечивает пренебрежимо малую погрешность от замены рядов с бесконечным числом слагаемых на их конечные суммы.
3,5
3,0
2,5
2,0
1,5
1,0
0,5
),0
- X р -ё-н-т
\ е = 0.0, 0.2,... 1.0
— без камеры гашения -одиночное судно | 1 к-Т
\0
1,0
1,5
2,0
2,5
3,0
Рисунок 4. Зависимость горизонтальной составляющей возмущающей силы от ширины камеры гашения
Рисунок 5. Зависимость горизонтальной составляющей силы от пористости лицевой грани причала
Рисунок 6. Зависимость горизонтальной составляющей возмущающей силы от осадки судна
Рисунок 7. Зависимость горизонтальной составляющей возмущающей силы от расстояния между судном и лицевой гранью причала
X р-g ■ h ■ Т
Г\
a = S, 12, 16
-1
1—( — >ез камеры гашения >диночное судно i k-T
"0,0 0,5 1,0 1,5 2,0 2,5 3,0
Рисунок 8. Зависимость горизонтальной составляющей возмущающей силы от
полуширины судна
На рисунках 4-8 представлены зависимости амплитуды X горизонтальной составляющей силы, воздействующей на судно, от воздействия морских волн, от различных геометрических параметров задачи и характеристик причала в широком диапазоне значений волнового числа к. Для удобства сопоставления на рисунки нанесены также зависимости величины X от к для одиночно стоящего судна (в отсутствие причала) и в ситуации, когда причал представляет собой только пористую проницаемую вертикальную стенку (без камеры гашения за ней). Оба эти варианта предполагают, что берег бесконечно далеко удален от судна и стенки. Расчеты данных случаев велись в рамках математических моделей, представленных в работах [22-24].
Видим, что наличие камеры гашения существенно сказывается в диапазоне значений к-Тот 0,1 до 0,5. Из зависимостей также видно влияние волногасящих свойств камеры. В случае, когда лицевая грань причала представляет собой непроницаемую вертикальную стенку (когда коэффициент пористости лицевой грани s = 0,0), судно испытывает большие нагрузки (для приведенного примера - приблизительно в 1,5 раза). Видим, что правильный выбор проектных размеров причала позволяет уменьшить нагрузки, испытываемые судном.
Заключение
В работе представлена математическая модель, описывающая воздействие морских волн на судно, стоящее у причала с камерой гашения. Модель учитывает геометрические размеры камеры гашения и физические характеристики лицевой пористой стенки (коэффициенты сопротивления и пористости, присоединенную массу). Представленные соотношения позволяют выполнить расчет нагрузок (возмущающих сил) на судно и его гидродинамических характеристик (присоединенные массы и коэффициенты сопротивления). В качестве примера представлены результаты расчетов по приведенной модели.
Литература
1. Макаров К. Н. Математическое моделирование в морской гидротехнике. Сочи: Изд-во Сочинского государственного университета туризма и курортного дела, 2008. 397 с.
2. РД 31.33.03-88. Рекомендации по определению допустимых ветро-волновых условий при проектировании морских портов. М., 1989. 80 с.
3. Frostick L.E., McLelland S.J., Mercer T.G. User guide to physical modelling and experimentation: Experience of the HYDRALAB Network. Leiden: CRC Press/Balkema, 2011. 245 p.
4. Fonseca N., Guedes Soares C. Experimental Investigation of the Nonlinear Effects on the Vertical Motions and Loads of a Containership in Regular Waves // Journal of Ship Research. 2004. Vol. 48. No. 2. Pp. 118-147.
5. Faltinsen O. Sea loads on ships and offshore structures. Vol. 1. Cambridge: Cambridge university press, 1993. 332 p.
6. Linton C.M., McIver P. Handbook of Mathematical Techniques for Wave/Structure Interactions // Chapman & Hall / CRC, 2001. 624 p.
7. Zheng Y.H., Shen Y.M., Tang J. Radiation and diffraction of linear water waves by an infinitely long submerged rectangular structure parallel to a vertical wall // Ocean Engineering. 2007. No. 34. Pp. 69-82.
8. Chong W., Watanabe E., Utsunomiya T. An eigenfunction expansion-matching method for analyzing the wave-induced responses of an elastic floating plate // Applied Ocean Research. 1995. No. 17. Pp. 301-310.
9. Vasyuk S.M., Maximov V.V., Nudner I.S., Puchkov V.S. To the problem of calculation of pitch of buoyant object moored to stationary and floating berth // International Symposium on Ship Hydrodynamics (ISSH) devoted to 85th anniversary of birthday of Abram M. Basin, St. Petersburg, 1995. Pp. 238-244.
10. Prpic-Orsic J., Slapnicar V., Turk A. Berth Operability Estimation Related to Ship Motion // Transactions of FAMENA. 2014. Vol. 38(1). Pp. 13-24.
11. Van Oortmerssen G. Forces related to motions of moored ships/analytical methods of moored ship motions // Advances in Berthing and Mooring of Ships and Offshore Structures. Springer, Netherlands. 1988. Pp. 265-281.
12. Ohkusu M. Ship motion in vicinity of a structure // Proceedings of the 1st International Conference on Behavior of Offshore Structures // The Norwegian Institute of Technology. 1974. Vol. 1. Pp. 284-306.
13. Lee Jaw-Fang. On the heave radiation of a rectangular structure // Ocean Engineering. 1995. Vol. 22(1). Pp. 19-34.
14. Hsu H.-H., Wu Y.-C. The hydrodynamic coefficients for an oscillating rectangular structure on a free surface with sidewall // Ocean Engineering. 1997. Vol. 24(2). Pp. 177-199.
15. Yang H.T.O., Huang L.H., Hwang W.S. The interference of a semi-submerged obstacle on the porous breakwater // Applied Ocean Research. 1997. No. 19. Pp. 263-273.
16. Hsu H.-H., Wu Y.-C. Scattering of water wave by a submerged horizontal plate and a submerged permeable breakwater // Ocean Engineering. 1999. No. 26. Pp. 325-341.
17. Ferreira M.D., Newman J.N. Diffraction effects and ship motions on an artificial seabed // Proceedings of the 24th International Workshop on Water Waves and Floating Bodies, 19-22 April, Zelenogorsk, Russia. 2009. 216 p.
18. Sakakibara S., Kubo M. Effect of structure types of quay wall on moored ship motions // Coastal Engineering Journal. 2008. No. 50(1). Pp. 101-122.
19. Хаскинд М. Д. Гидродинамическая теория качки корабля. М.: Наука, 1973. 328 с.
20. Ньюмен Дж. Морская гидродинамика / Перевод с английского. Л.: Судостроение, 1985. 368 с.
21. Yu X. Diffraction of water waves by porous breakwaters // Journal of Waterway, Port and Coastal Ocean Engineering. 1995. Vol. 121(6). Pp. 275-282.
22. Леонтьев В.А., Нуднер И.С., Семенов К.К. Взаимодействие волн с плавающим телом, расположенным у пористой стенки // Труды XI Всероссийской конференции «Прикладные технологии гидроакустики и гидрофизики». СПб.: Наука, 2012. С. 204-207.
23. Леонтьев В.А., Нуднер И.С., Семенов К.К. Взаимодействие морских волн с судном, раскрепленным у причала // Труды Международной конференции «Математические и информационные технологии, MIT-2013», Врнячка Баня, Сербия - Будва, Черногория. 2013. С. 386-396.
24. Shchemelinin L.G., Leontiev V.A., Nudner I.S., Semenov K.K. Berth influence on wave loads on the moored ship near it // Proceedings of the 25th (2015) International Ocean and Polar Engineering Conference ISOPE-2015, to be published.
Константин Константинович Семенов, г. Санкт-Петербург, Россия Тел. моб.: +7(911)7320171; эл. почта: semenov.k.k@gmail.com
Виктор Александрович Леонтьев, г. Санкт-Петербург, Россия Тел. моб.: +7(911)7057292; эл. почта: leontiev.vict@yandex.ru
Игорь Сергеевич Нуднер, г. Санкт-Петербург, Россия Тел. моб.: +7(911)7528499; эл. почта: igor_nudner@mail.ru
© Семенов К.К., Леонтьев В.А., Нуднер И.С., 2015