CC BY
41
|Выпуск 2
Список литературы
1. Михальский В. А. Метрология в кораблевождении и решение задач навигации / В. А. Михальский,
B. А. Катенин. — СПб.: Элмор, 2009. — 288 с.
2. Груздев Н. М. Теория навигационных погрешностей / Н. М. Груздев. — СПб.: СПбВМИ, 2002. —
325 с.
3. Лаврентьев А. В. Технические средства кораблевождения / А. В. Лаврентьев, С. Ф. Лабутин, А. Б. Парамонов, А. И. Тузов: учебник. — М.: Морские средства навигации надводных кораблей, 2007. —
C. 101-114.
4. Михалъскии В. А. Теория и модели погрешностей навигационных измерений / В. А. Михальский,
A. Г. Кушнарев, М. В. Тюлькин [и др.] — СПб.: ВСОК ВМФ, 2006.
5. Пешехонов В. Г. Единая система инерциальной навигации и стабилизации «Ладога-М» / В. Г. Пе-шехонов, Б. Л. Шарыгин, Ю. В. Миронов // Морская радиоэлектроника. — 2003. — № 1. — С. 26-30.
6. Бабиченко А. В. Инженерная физика / А. В. Бабиченко, В. К. Шкред // Научно-технический журнал. — 2011. — № 11. — С. 33-52.
7. Иванов Б. Е. К вопросу о навигационной безопасности плавания / Б. Е. Иванов / Пятая Российская науч.-практ. конференция «Современное состояние и проблемы навигации и океанографии» (НО-2004). — СПб.: ГНИНГИ МО РФ, 2004. — С. 74-86.
8. Соколов А.И. Навигация и управление движением/А. И. Соколов, А. С. Фирса/Материалы XV конференции молодых ученых «Обработка сигнала в электромагнитном лаге». — СПб., 2013. — С. 88-94.
9. Смирнов Е. Л. Технические средства судовождения. — Т. 2. / Е. Л. Смирнов, А. В. Яловенко,
B. К. Перфильев. — СПб.: Элмор, 2000.
10. Борисов Д. В. Возможности приёмоиндикаторов СРНС ГЛОНАСС и NAVSTAR, установленных на кораблях и судах ВМФ / Д. В. Борисов // Научно-методический сборник. — СПб.: СПбВМИ, 2006. —
C. 54-60.
11. Алексеев С.П. Об интеграции информационного обеспечения задач навигации, стабилизации и управления движением морских подвижных объектов / С.П. Алексеев, Г.И. Емельянцев // Навигация и гидрография. — 1996. — № 2. — С. 73-76.
12. Горкавенко Ю. Н. Радиотехнические и астронавигационные средства кораблевождения. — Ч. 1 / Ю. Н. Горкавенко, А. В. Лаврентьев, Н. Б. Лукашевич и др. — СПб.: СПбВМИ, 2001. — С. 84-99.
13. Горкавенко Ю. Н. Основы устройства и использования навигационной аппаратуры потребителя «Бриз» среднеорбитных спутниковых систем ГЛОНАСС и NAVSTAR / Ю. Н. Горкавенко, А. А. Опалев, М. Ю. Шелякин. — СПб.: СПбВМИ, 2001. — С. 63-72.
14. Соловьев Ю. А. Спутниковая навигация и ее приложения / Ю. А. Соловьев. — М.: Эко-Трендз, 2003. — 326 с.
УДК 532.516.5:532.582.32
А. В. Коптев
канд. физ.-мат. наук, доц.
ИССЛЕДОВАНИЕ ВЛИЯНИЯ ПОДВОДНОГО ТЕЧЕНИЯ НА ДИНАМИКУ КОРАБЛЯ
REAEARCH OF INFLUENCE OF UNDERWATER CURRENT ON THE DYNAMICS OF THE SHIP
Предложен аналитический метод расчета воздействия подводного течения на погруженную часть корабля. Метод основан на решении краевой задачи для уравнений установившегося движения вязкой несжимаемой жидкости при граничных условиях прилипания и условии для скорости водного потока в задан-
ВЕСТИ IT
ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВ^
ной внешней точке. В качестве определяющих уравнений использовался интеграл уравнений Навье - Стокса для 2D установившегося движения. В результате решения задачи получены новые формулы, определяющие вертикальную и горизонтальную составляющие результирующей силы как функции исходных параметров, задающих геометрические размеры обтекаемого тела, а также интенсивность и направление подводного течения. Произведены вычисления с помощью пакета стандартных программ Maple и найдены оценки для значений гидродинамических сил.
On the paper under consideration we offer the analytical computational method for estimation of influence of underwater current upon the submerged part of the ship. Method is based on the decision of boundary value problem for equations of viscous incompressible steady state fluidflow with adhesion condition together with the condition for water-course velocity in the set external point. We take the integral of 2D steady state Navier - Stokes equations as the base equalizations. As a result of problem decision the new formulas for vertical and horizontal constituents of resulting force were obtained. The latter determines lift force and front resistance as functions of geometric dimensions of streamlined body together with the intensity and direction of underwater currents. Calculations are produced by means of routine package of Maple programs and values of hydrodynamic forces were found out.
Ключевые слова: судно, обтекание, вязкая несжимаемая жидкость, уравнение, интеграл, подъемная сила, лобовое сопротивление.
Key words: ship, overflow, viscous incompressible fluid, equation, integral, lift force,frontal resistance.
сил [1]. Так, основными силами, действующими на судно в вертикальном направлении, являются сила тяжести, направленная вниз, и выталкивающая архимедова сила, направленная вверх. При движении судна эти силы в среднем уравновешивают друг друга. Однако в зоне действия подводного течения, когда происходит обтекание скоростным водным потоком погруженной части судна, баланс сил изменяется, возникают дополнительные силы, имеющие гидродинамическую природу, такие как подъемная сила {Lift Force), направленная вертикально вверх, и так называемая сила лобового сопротивления {Frontal Resistance), направленная горизонтально. В общем балансе сил их влияние может быть довольно значительным. С одной стороны, судно получает дополнительную положительную плавучесть и возникает тенденция всплытия, а с другой — возникает дополнительный крен и дифферент. Указанные факторы особенно важны для высокосидящих судов, таких как круизные лайнеры и морские паромы. В истории мореплавания известны случаи, когда влияние подводного течения оказывалось решающим. Судно получало крен больше критического, теряло остойчивость, что в результате приводило к катастрофе. Даже в недалеком прошлом такие катастрофы имели место. Достаточно вспомнить гибель паромов «Эстония» на Балтике (1994 г.), «Коста Конкордия» в Тирренском море (2012 г.), «Севоль» в Желтом море (2014 г.).
Чтобы учесть, в полной мере, влияние подводного течения на динамику корабля, необходимо на основе уравнений гидромеханики исследовать обтекание погруженной части корабля водным потоком и рассчитать возникающие при этом гидродинамические силы. Для решения этой задачи предлагается пренебречь некоторыми конструктивными особенностями днища корабля и в качестве простейшей модели погруженной части рассмотреть часть кругового цилиндра радиуса R и углом раствора 2а. Предлагается рассматривать вариант плоской задачи, когда подводное течение действует в поперечном направлении и картина течения одинакова во всех плоскостях, перпендикулярных главной оси судна. В сечении погруженной части цилиндра перпендикулярной плоскостью имеем сегмент круга радиуса R с углом раствора 2а, где а — некоторое заданное значение угла из диапазона 0 < а < к. Контуром обтекаемой поверхности является нижняя круговая часть этого сегмента, а верхняя его часть соответствует уровню свободной поверхности. Считаем, что происходит безотрывное ламинарное обтекание погруженной части судна внешним потоком жидкости, которую предполагаем вязкой и несжимаемой [2] - [4]. В качестве жидкости
ПРАКТИКЕ мореплавания нередко встречаются ситуации, когда приходится осуществлять навигацию в районах действия интенсивных подводных течений. Подводные течения воздействуют на погруженную часть судна и влияют на баланс действующих
Выпуск 2
|Выпуск 2
рассматриваем морскую воду с заданными и неизменными физическими свойствами. Считаем, что р — плотность, a v — коэффициент кинематической вязкости.
Задачу рассматриваем в декартовых координатах, оси которых направим через центр окружности, получающейся в результате сечения цилиндра поперечной плоскостью. Начало координат, таким образом, отстоит от свободной поверхности на величину R cos а. Ясно, что в выбранной системе координат эта величина, в общем случае, может быть и положительной, и отрицательной,
и нулевой. Первая возможность реализуется при ^<а<п, вторая — при 0 < а< ^ и третья — при п 2 2
а--.
Условия для набегающего водного потока зададим в произвольно выбранной точке М расположенной вне цилиндра, на заданной глубине. Пусть абсцисса точкиМ0 равна Z, а ордината — И, где L иЯ представляют некоторые положительные параметры. Так что точка М0 находится на расстоянии L с левой стороны от вертикальной оси на глубине И - R cos а, считая от свободной поверхности. В точке М0 полагаем заданными и продольную, и поперечную скорость движения водного потока.
Задачу рассматриваем в безразмерных переменных, выбрав удобные для данного случая масштабы. Пусть U0 — модуль скорости жидкости в точкеМ0. Эту величину выбираем в качестве
U2
масштаба скорости. Заданную величину Z полагаем масштабом длины,-^ — масштабом ускорения, а произведение рU0 — масштабом давления.
В безразмерных переменных координаты точки М0 определяются равенствами хм = -1,
ум = —h, где h = L, Контур обтекаемого тела представляется частью окружности с уравнением
х2 + у2 = г2, где г = l — безразмерный радиус цилиндра, 0 < г < 1. Свободная поверхность соответствует плоскости, определяемой уравнением у=у0, гдеу0 = -гcos а. Интересующая нас область течения расположена ниже этой плоскости, так что вертикальная координата должна удовлетворять неравенству
J<J0- (!)
Задачу предлагается рассматривать на основе 2D уравнений Навье - Стокса для установившегося движения вязкой несжимаемой жидкости [2] - [4]. Записанные в безразмерных переменных, они имеют вид:
ди ди д( р + Ф) 1
U ~дх + V ~ду ~ дх + RE
^ д 2u
д 2u ^
дх2 + ду2
(2)
дх дх д( р + Ф) 1 U дх + V ду ~ ду + RE
2я2
д х
д V
дх2 + ду2
(3)
ди дх _ „
дх + ду '
(4)
Здесь них — соответственно продольная и поперечная скорости движения водного потока; р — давление; Re — число Рейнольдса, определяемое формулой Re ~~~~Ф — потенциал внешних сил. Для рассматриваемого случая Ф = g(y -у(), где g — безразмерное ускорение свободного падения,^ — вертикальная координата.
Уравнения (2) - (4) следует дополнить граничными условиями, которые предлагается задать следующим образом. Первое граничное условие состоит в том, что вдоль контура обтекаемого тела продольная и поперечная компоненты скорости обращаются в нуль. В этом состоит так называемое условие прилипания при обтекании тела вязкой жидкостью [2] - [4]. В принятых ранее обозначениях оно сводится к двум равенствам:
U X2+у2=r2
= 0;
V X2+y2=r2
= 0.
(5)
Предполагается, что эти равенства выполнены лишь для нижней обтекаемой части сегмента, для которой справедливо неравенство -г < у <у0. Второе граничное условие задает скорость набегающего водного потока вне обтекаемого тела и задается в точкеМ0(-1; -И). В этой точке должен быть задан вектор скорости. В соответствии с выбором масштабов модуль этого вектора (величина U0) задает масштаб скорости. Поэтому достаточно задать лишь угол наклона данного вектора или его тангенс. Последняя величина равна отношению поперечной и продольной скорости в точ-кеМ0. Будем считать, что тангенс угла наклона вектора скорости равен заданной величине k = tgQ, где 0 — значение угла наклона. Таким образом, второе граничное условие соответствует равенству
у(-!;-&) _ , и (-1; -h)
(6)
где Иик — заданные безразмерные параметры.
Чтобы упростить задачу, предлагается исходить не из уравнений Навье - Стокса (2) - (4) непосредственно, а из первого их интеграла. Эти соотношения, предложенные в работах [5] и [6], выведены с учетом сохранения нелинейных членов и представляют уравнения более низкого порядка относительно основных неизвестных. Кроме того, они определяют общую структуру решений [10], [11] и открывают новые возможности для построения решений [7], [12]. Указанные обстоятельства делают данные соотношения более простыми и удобными в математическом плане, чем исходные уравнения (2) - (4). Для рассматриваемого случая 2D установившегося движения первый интеграл уравнений Навье - Стокса сводится к трем равенствам:
U 2
p + g(y - Уо) + -у+d = 0; (7)
2 2 2 ( du dv
u -v + rE+ду
д 2^!2 ,3^12 . дл2 ду2 ’
(8)
1 (ды 5у) _ д2^j2 UV Re^ду + дх) длду '
(9)
В уравнениях (7) - (9) использованы следующие условные обозначения: ¥ — новое ассоциированное неизвестное; Uи4 — соответственно модуль скорости точки водного потока и диссипация вследствие трения, вычисляемые по формулам:
U = Vи2 + v2 ;
d =
1 (д2^12 . д2^12
2 ^ дх2 ду2
(10)
Произвольные функции одной переменной, возникающие в общем случае при интегрировании в уравнениях (7) и (8), полагаем равными нулю. К уравнениям (7) - (9) следует добавить уравнение неразрывности (4) и граничные условия (5) и (6). Таким образом, в математическом плане постановка задачи сводится к решению четырех дифференциальных уравнений: (4), (7) - (9), относительно неизвестных м, v,p, ¥ при граничных условиях (5) и (6). Предлагаемая постановка требует также задания пяти исходных параметров: угла полураствора а, числа Рейнольдса Re, безразмерных радиуса цилиндра г и глубины И для выбранной характерной точки М а также числа к, определяющего направление вектора скорости в точке М
Силы гидродинамического воздействия на погруженную часть корабля можно найти, если решена краевая задача для уравнений (4), (7) - (9) при условиях (5) и (6). Имеем задачу обтекания части цилиндра установившимся потоком вязкой несжимаемой жидкости. Задачи этого типа представляют известный класс задач теоретической гидромеханики и имеют практическое значение. Начиная с Дж. Стокса, задачи такого типа рассматривались разными авторами и для их решения применялись различные подходы. Однако и на сегодняшний день общепринятого варианта решения нет. Основные сложности возникают из-за нелинейности уравнений и из-за задания гранич-
Выпуск 2
|Выпуск 2
ного условия для набегающего потока на бесконечности [2] - [4]. При граничных условиях, аналогичных условиям (5) и (6), задача рассматривалась в работе [8]. Было предложено использовать первый интеграл уравнений Навье - Стокса, тем самым понижая порядок нелинейных членов и задавать граничное условие для набегающего потока в конечной точке. В указанной работе имело место полное обтекание цилиндра, когда роль контура обтекаемой поверхности играла целая окружность х2 + у2 = г2. ТочкаМ0 задавалась координатами х |Мо = -1, у |Мо = 0, т. е. она была расположена в плоскости симметрии обтекаемого тела.
В настоящей работе постановка задачи содержит два существенных отличия. Во-первых, контур обтекаемой поверхности представляет лишь часть окружности в соответствии с ограничением -г<у< -г cos а. Имеем не полное, а лишь частичное обтекание контура, задаваемое углом полураствора а. Во-вторых, точка М0 расположена на произвольной глубине и в общем случае несимметрично по отношению к обтекаемой поверхности. В соответствии с условием y |м = —h вводится дополнительный параметр глубины h, что дает возможность рассматривать подводные течения. Таким образом, в постановку задачи по сравнению с [8] внесено изменение во втором граничном условии (условии (6)) и вводятся два дополнительных параметра: h и а. Однако, несмотря на эти различия, для рассматриваемого случая начальные этапы решения совпадают с соответствующими этапами работы [8], и некоторыми её результатами можно воспользоваться. В частности, можно использовать выражения для поля скорости водного потока [8, формулы (13) - (14)] и для давления [8, формула (21)]:
= — -(г2 - х1 - у1 )\_а01г2 + апг 2л + 2а02г 2у - а01х2 - 4а10 ху - 5а01 у 2~ап х3 +
+х2у ■ l(a22r2 + а02)- 5апху2 + у3 ■ l(a22r2 - а02 + 2а20)];
(11)
и
г = —4 -{г2 - х2 - у2 )-[а10г2 + 2а20г 2х + апг 2у - 5а10х2
4а01ху - а10у2 +
+х3 ■ 2(а22г2 - а20 + 2а02)- 5апх2у + ху2 • 2(a22r2 + а20)- а11 у3]; (12)
16a°i. j + у (2а20а01 - апаю + j +
Р = -g(y - Уо) +1 2b20 - а120 - ^ I - *1 2a2QaiQ + aiia01 +
Re
Re r2
21 о 2 Aii 2 f 2 2 6a,,44 8ху f a2n4 2 f2a,o 2ao, 2 a,, 12a,,
+x 2| -2a2n+ — I a0i - a2 -^ I I-fl an,a,n|+ у 2| —^-^ - 2al -^ +
Re
I a0lal0- re I+у
22 r2 r
2 Re r
+
+
4x3 (a20a10 + ^ + 2a01 ) 2x2у („ + 12a10 4 2xy2 („_________ + 12a01
2 l 3
- + ana0i +
Re r2
8a2oaoi + апаю +
Re r2
+ -
8a2oaio anaoi +
Re r2
r
+
2
2
r
r
r
r
+^T (aiiaio - Ще0- - REr2 ) + Rf7ai2°- _a°2i +2r2a“ +4r^ +it) + ('
11 1 + 8x_y 7 Re )+ r4 l_
8a20
Re
4x2 у2
'у 2 2\ 8xy3 (7 2
+ -R^~(«01 + "10 l 3 A10"01 + Г a2°an -
$a2o) + y_(7a2i -3"2 + 4r2a20 + 2r2a121 - |. (13)
Re
Re
Приведенные выражения являются результатом решения определяющих уравнений с учетом лишь одного граничного условия - условия прилипания. Они справедливы и в рассматриваемом случае, с той лишь разницей, что в правой части выражения для р следует полагать у0 = -г cos а. Выбор аддитивной постоянной давления в (13) можно осуществить аналогично ра-
1 А
боте [8], приняв Ь20 = ^ «т + iT1 > а также используя полученные в работе [8, формулы (16) и (17)]
Re
выражения для коэффициентов разложения а.:
А10 —
24 г Re
Сп, Ао
— ±tk'l1 - с2; a22 — 0;
"г Re
2
24
г2 Re
(2со l); a02 — ± 2r Сол!1 ^о"; a20 —h 2r солД ^о"■
(14)
В выражении (14) с. представляет некоторую величину, удовлетворяющую неравенству
-1<с0<1.
(15)
Однако на этапе вычисления подъемной силы и силы лобового сопротивления имеем существенные отличия от работы [8]. Для рассматриваемого случая обтекаемый контур представляет лишь часть окружности, и выражение (13) следует интегрировать вдоль интервала углов ф, соответствующего неравенству -г < у < -гcosa. В качестве такого интервала использован
-П-а; _n + aj- С учетом этого обстоятельства исходные формулы для подъемной силы и силы
лобового сопротивления [4], возникающих вследствие обтекания водным потоком, приобретают вид:
п
~2 +а
п
~2+а
Fy =- | p(rcosф; rsinф)-rsinфdф; Fx = -J p(rcosф; rsinф)-rcosфdф.
(16)
Интегрирование по этим формулам с учетом выражений (13) и (14) приводит к следующим результатам:
576
Fy = gr2 (a- 0,5sin2a)---2 (-ll,66sin a + 6,llsin3 a-l,6sin5 a + c0 (1,67a-l,67sin2a-l,25sin 4a)
y r Re
+
+c0 (17,33 sin a -12,22 sin3 a + 3,2 sin5 a) + c0 (1,67 sin 4a+3,34sin 2a)
+co
(-3,2sin5 a + 5,34sin3 a - 8sin a));
Fx =
480^/T-
co
r Re2
(2a + 0,2 sin 2a + 0,5 sin 4a + c0 (l6,8sin3 a - 7,68sin5 a) - 2c2 sin 4a +c0 (30,72sin5 a - 25,6sin3 a)) •
(17)
(18)
Формулы (17) и (18), которые определяют выражения для гидродинамических сил при произвольных значениях угла полураствора обтекаемой части 0 <а<п, являются новыми. В частном случае (а = п) эти формулы приводят к выражениям для Fx, F совпадающим с выражениями для случая полного обтекания кругового контура [8, формула (25)].
Проанализируем правые части формул (17) и (18). Ясно, что за исключением особо редких случаев, F и F отличны от нуля, а значит, парадоксов типа Стокса или Даламбера [2] - [4] в данном случае нет. В правых частях рассматриваемых формул присутствуют величины, задаваемые в качестве исходных, а также другие величины, требующие уточнения. Из заданных изначально величин там присутствует число Рейнольдса Re, радиус обтекаемого цилиндра г и угол полураствора обтекаемой части а. Величина с0 нуждается в уточнении. Ее значение следует определять с помощью граничного условия (6). Заметим также, что первая группа членов в правой части выражения для F дает безразмерное значение архимедовой силы, действующей на погруженное в жидкость тело. Архимедова сила пропорциональна объему погруженной части тела и направлена вверх. В данном случае, когда рассматривается плоская задача, эта величина будет пропорциональна площади поперечного сечения погруженной части, т. е. площади сегмента окружности радиуса г и углу раствора 2а. Площадь этой фигуры равна r2(a - 0,5sin2a). Именно такой множитель и наблюдается в первом положительном слагаемом правой части формулы (17).
Известно, что архимедова сила и сила тяжести являются статичными и никак не связаны с движением тела. В противоположность этому сила воздействия водного потока представляет динамическую составляющую, которая целиком определяется движением водного потока. Ука-
а
а
а
2
2
+
Выпуск 2
|Выпуск 2
занная динамическая составляющая вносит свой вклад в правую часть формулы (17) и целиком определяет правую часть (18). Каждое из этих слагаемых зависит от с Определим значение этой величины, исходя из условия (6). С учетом выражений (11) и (12) для м, v при х = -1 ,y = —h составляем дробьУ( 1; и приравниваем ее заданной величине Уходя от знаменателя и используя и (-1; -п)
выражения (14) и (15), получаем равенство, равносильное (6):
е1со + ЕСо + Е = 41 - Со2 (eCo + Е) • (19)
Коэффициенты е. в (19) определяются через заданные параметры согласно формулам:
, , (к + Н + 5kh2 + 5h) 2 2 п t\ (Н3 + 5kh2 + 5h + к)
el = 2r(h + к)-2------------; e2 = 5 -r2 + h2 + 4kh; e3 = -r(к + h) + --------;
e4 = 2r (1 - kh) + 2
(3h3k + h2 - kh - 3).
; e5 =-5kh2 - 4h - к (1 - r2).
(20)
Уравнение (19) представляет иррациональное уравнение относительно с0. Чтобы его разрешить, необходимо избавиться от радикалов путем возведения в квадрат. Однако впоследствии следует произвести отбор корней, исходя из двух положений. Первое состоит в выполнимости ограничения (15), а второе — в том, что левая и правая части равенства (19) должны иметь одинаковые знаки. Для этого должно выполняться неравенство
Е1С0 + е2С0 + Е3 > о Е4С0 + Е5
(21)
Таким образом, в качестве критериев отбора корней при решении уравнения (19) имеем неравенства (15) и (21). Далее возводим равенство (19) в квадрат. После перегруппировки слагаемых приходим к уравнению четвертой степени относительно с которое можно представить в виде
Со4 + ?А3 + ^о2 + ?з<4) + ?4 = 0 •
(22)
Коэффициенты q. последнего уравнения определены через коэффициенты е. уравнения (19)
по формулам:
= 2-
ег + е;
?2 =
2е,е3 + е,2 + е, - е42 , _ 0 е2еъ - е4е5 ^ ^ - е5
’ ?3 _ 2----2 , 2 ’ ^4 2 , „2 •
^ ^ с4 ej2 + е4 + е4 ь4
Уравнение (22) есть приведенное уравнение четвертой степени, которое может быть сведено к уравнениям более низких степеней, если воспользоваться методом Феррари и формулами Кор-дано [9]. Таким образом, данное уравнение можно разрешить аналитически, тем самым определив с0 через исходные величины. В результате задача определения гидродинамических сил в общем виде будет решена.
Для практических целей удобно использовать численное решение уравнения (22). Этого можно достигнуть с помощью стандартных программ пакета Maple. В данной работе использована именно эта возможность. Приведем результаты расчетов. При некоторых значениях исходных параметров было произведено вычисление величины с0 и всех других интересующих величин, включая F и F Величины, определяющие масштабы, были выбраны как L= 10 м и U0 = 2 м/с.
В качестве исходных параметров были взяты а = — , к = 1. Параметр г был принят равным 0,5,
что в размерных единицах соответствует R = 5 м. Этому значению радиуса соответствует судно среднего водоизмещения с поперечным размером порядка 10 м. Параметр h полагался равным 1,5, так что точкаМ0, в которой задается второе граничное условие, находится на глубине, равной
1,5-10 - 5cosП = 12,46 м.
4
Безразмерные коэффициенты уравнения (22), найденные по формулам (21) и (23), равны = -0,4261, q2 = -0,9385, = 0,2746, = 0,1877. Численное решение уравнения (22) с помощью
программ Maple дало четыре возможных вещественных значения: с® = 0,8297, с02) = 0,7673,
е, + е;
ег + е;
(23)
r
с® = -0,3665, С04) = -0,8044, однако критериям отбора (15), (21) удовлетворяют не все из них. Если неравенству (15) удовлетворяют все четыре значения, то неравенству (21) — только Таким образом, в качестве решения уравнения (19) имеем только значение с0 = -0,8044.
Подъемная сила и лобовое сопротивление, вычисленные по формулам (17), (18), в безразмерных единицах получились следующими: F = 1,748 + 1,469 • 1041; Fx = 4,383 • 1042. Расчеты показали, что Fx является сравнительно малой величиной и что для Fy основной вклад дает первое слагаемое, определяемое архимедовой силой, тогда как вклад от динамической составляющей, вызванной подводным течением, сравнительно невелик. Однако при уменьшении величины а и значения параметра h эта тенденция начинает изменяться и вклады от гидродинамических составляющих становятся все более заметными. Они также становятся более значимыми и при увеличении интенсивности подводного течения, определяемого величиной U Обращает на себя внимание также отличный от нуля момент силы F который оказывает явное влияние на остойчивость судна в поперечном направлении.
Вывод. Предлагаемая модель и методика расчета позволяют на практике определять гидродинамические силы, возникающие вследствие подводного течения при обтекании подводной части корабля, и оценивать риски, связанные с их возникновением.
Список литературы
1. Анфимов В. Н. Устройство и гидромеханика корабля / В. Н. Анфимов. — Л.: Судостроение, 1974. —
368 с.
2. Кочин Н. Е. Теоретическая гидромеханика / Н. Е. Кочни, И. А. Кибель, Н. В. Розе. — М.: Изд-во физ.-мат. лит., 1963. -4.2. — 727 с.
3. Лойцянский Л. Г. Механика жидкости и газа / Л. Г. Лойцянский. — М.: Наука, 1987. — 833 с.
4. Валландер С. В. Лекции по гидроаэромеханике / С. В. Валландер. — Л.: Изд-во Ленингр. госуд. унта им. А. А. Жданова, 1978. — 294 с.
5. Коптев А. В. Первый интеграл и пути дальнейшего интегрирования уравнений Навье - Стокса / А. В. Коптев // Известия Российского государственного педагогического университета им. А. И. Герцена. Серия: Естественные и точные науки. — 2012. — № 147. — С. 7-17.
6. Коптев А. В. Как проинтегрировать уравнения Навье - Стокса / А. В. Коптев // Физическая механика. — Вып. 8. — СПб.: Изд-во СПбГУ, 2004. — С. 218-226.
7. Коптев А. В. Решение начально-краевой задачи для 3D уравнений Навье - Стокса и его особенности /А. В. Коптев // Изв. Российского государственного педагогического университета им. А. И. Герцена. Серия: Естественные и точные науки. — 2014. — № 165. — С. 7-18.
8. Коптев А. В. Динамические реакции подводного трубопровода на морские течения / А. В. Коптев // Вестник государственного ун-та морского и речного флота имени адмирала С. О. Макарова. — 2014. — № 4. — С. 107-114.
9. Кочаргин И. Ф. Алгебраические уравнения / И. Ф. Кочаргин. — М.: Изд-во Физматкнига, 2006. —
160 с.
10. KoptevA. V. Structure of solutions of 3D Navier- Stokes equations / A. V. Koptev // Дифференциальные уравнения и динамические системы: материалы Междунар. конф. / Математический ин-т им. В. А. Стеклова РАН. — 2014. — С. 212-214.
11. KoptevA. V. The Structure of Solution of the Navier - Stokes Equations / A. V. Koptev // Вестник национального исследовательского ядерного университета «Московский инженерно-физический институт». — Вып. 3 (6). — 2014. — С. 656 660.
12. Koptev А. V. Generator of solutions for 2D Navier - Stokes equations /А. V. Koptev // Journal of Siberian Federal University. Mathematics & Physics, Issue 7(3). — 2014. — P. 324-330.
Выпуск 2
