УДК 531.25 Вестник СПбГУ. Математика. Механика. Астрономия. 2024. Т. 11 (69). Вып. 1 МБС 70С20
Устойчивость плавания судов с сечениями в форме эллиптического и гиперболического сегментов
А. С. Смирнов1'2, И. А. Кравчинский2'3
1 Санкт-Петербургский политехнический университет Петра Великого, Российская Федерация, 195251, Санкт-Петербург, ул. Политехническая, 29
2 Институт проблем машиноведения Российской академии наук, Российская Федерация, 199178, Санкт-Петербург, Большой пр. В. О., 61
3 ООО «Центр технологического консалтинга»,
Российская Федерация, 194044, Санкт-Петербург, Большой Сампсониевский пр., 28/2
Для цитирования: Смирнов А. С., Кравчинский И. А. Устойчивость плавания судов с сечениями в форме эллиптического и гиперболического сегментов // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2024. Т. 11 (69). Вып. 1. С. 171-184. https://doi.org/10.21638/spbu01.2024.112
В статье рассматриваются две задачи об устойчивости тривиального положения равновесия плавающих судов с сечениями в форме эллиптического и гиперболического сегментов. Дается обзор примеров на устойчивость плавания тел и излагаются ключевые принципы ее исследования методами аналитической статики. Для каждой из представленных задач путем достаточно серьезных математических построений получается точное выражение для потенциальной энергии в рамках принятой конфигурации и вычисляется его квадратичная аппроксимация вблизи исследуемого равновесного состояния. На ее основе устанавливаются условия устойчивости положения равновесия в терминах трех безразмерных параметров, а также анализируются предельные случаи. Осуществляется сопоставление промежуточных выражений и конечных результатов, полученных в процессе обсуждения каждой из задач, и выявляются их общие черты и отличительные особенности. Найденные решения иллюстрируются в виде семейств границ областей устойчивости на плоскости двух безразмерных параметров при различных значениях третьего параметра. Эти результаты представляют фундаментальное теоретическое значение и могут оказаться полезными для практических приложений.
Ключевые слова: плавающее судно, эллиптический и гиперболический сегменты, устойчивость, статический анализ, граница области устойчивости, плоскость безразмерных параметров.
1. Введение. Одним из ключевых вопросов в корабельном деле является исследование устойчивости плавающих тел в рамках одного из важнейших разделов механики — гидростатики, причем в морской науке для этого понятия используется специальный термин «остойчивость» [1, 2]. Как известно, основы статики корабля были заложены еще Архимедом, а существенный толчок в развитии учения об устойчивости плавающих тел произошел благодаря трудам Л. Эйлера и его многочисленных последователей. Подобные задачи имеют, помимо фундаментального теоретического значения, и широкое практическое применение, поскольку плаваю© Санкт-Петербургский государственный университет, 2024
щие в жидкости тела могут моделировать самые разнообразные водные транспортные средства, имеющие различное назначение, — понтоны, баржи, суда, лодки и т. п. [3, 4].
Особую важность при этом приобретает исследование устойчивости плавающих тел в пространстве параметров задачи, которое позволяет выявить характер равновесных положений при плавном изменении этих параметров [5]. Это обстоятельство зачастую и приводит к неожиданной потере устойчивости водного объекта с последующим его опрокидыванием или переворотом на бок [6], что требует весьма тщательного анализа.
В то же время статический анализ тела, погруженного в жидкость, даже при весьма несложной геометрии его сечения далеко не всегда удается осуществить до конца аналитическими методами, особенно если необходимо провести глобальный анализ, который подразумевает определение всех возможных положений равновесия плавающего тела и исследование их устойчивости [6, 7]. Поэтому в современных руководствах и задачниках по механике твердых тел и теории устойчивости присутствует весьма небольшое количество задач гидростатической устойчивости [8, 9], что и побуждает лишний раз обратить серьезное внимание на эту область.
Наиболее известными задачами об устойчивости плавающих тел являются задача о плавании понтона с прямоугольным сечением, решенная А. П. Фан-дер-Флитом [10, 11], и задача о плавании лодки с треугольным сечением, связанная с именем А. Ю. Давидова [9, 12]. В качестве развития первой из них можно указать статью [5], где рассматривается плавание понтона со смещенным центром масс, а также работу [6], где осуществляется глобальный статический анализ плавающей прямоугольной баржи.
Еще больший практический интерес представляют задачи об устойчивости плавания барж с трапециевидным и пятиугольным сечениями, которые сочетают геометрические особенности как задачи о плавании прямоугольного понтона, так и задачи о плавании треугольной лодки [13]. Это связано с тем, что тела с упомянутыми сечениями представляют собой несколько более удачное приближение к реальным плавающим средствам.
Однако существенно важный шаг при моделировании плавающих тел в жидкости заключается в переходе от брусьев с кусочно-прямолинейными сечениями, о которых шла речь во всех указанных выше вариантах, к брусьям с сечениями, границы которых содержат криволинейные участки [2], и от общих теоретических сведений к рассмотрению конкретных задач [7]. Одним из таких примеров является задача об устойчивости плавания судна с сечением в форме параболического сегмента, решенная авторами в работе [14], где помимо исследования устойчивости тривиального положения равновесия определяются и наклонные равновесные состояния.
В настоящей работе предлагается продолжить рассмотрение плавающих тел с сечениями, представляющими собой сегменты кривых второго порядка, а именно исследовать устойчивость плавающих судов с сечениями в форме эллиптического и гиперболического сегментов.
2. Исследование устойчивости плавающих тел методами аналитической статики. Прежде чем переходить к постановке и решению указанных задач, целесообразно указать основные положения, которые обычно принимаются при исследовании устойчивости плавающих тел методами аналитической статики. В связи
с тем, что подробное изложение этих моментов приведено в [5, 7, 9, 13], ограничимся здесь лишь их достаточно краткой сводкой. Предположим, что плавающее средство моделируется при помощи длинного бруса постоянного поперечного сечения с плотностью ръ (поверхностная плотность сечения бруса, т. е. масса на единицу площади), погруженного в жидкость плотности pi (также на единицу площади). В силу сказанного можно ограничиться рамками плоской задачи, при этом горизонтальная координата центра масс бруса в статической задаче отбрасывается как несущественная, а вертикальная координата исключается путем рассмотрения лишь равнообъемных наклонений бруса по углу поворота, что является стандартным приемом в статике плавающих тел [2, 15, 16]. В результате в рассмотрении остается лишь одна обобщенная координата — угол наклона бруса а задача сводится к одностепенной. Для исследования устойчивости необходимо записать и проанализировать выражение для потенциальной энергии плавающего бруса, которая равна П = P(уъ — yi), где P = Pi = Ръ, причем Ръ и Pi — веса бруса и вытесненной жидкости соответственно, равные друг другу, а уъ и yi — вертикальные координаты центров масс бруса и вытесненного объема, которые отсчитываются от ватерлинии [17]. При этом следует отметить, что условие равнообъемности в рамках плоской задачи приводится к виду Si = рБъ = const, где Бъ — площадь сечения бруса; Si — вытесненная площадь жидкости, а р = ръ/pi — безразмерная плотность, причем из условия плавания ръ < pi, так что 0 < р < 1.
3. Устойчивость плавания судна с сечением в форме эллиптического сегмента. В качестве первой задачи рассмотрим судно с сечением в форме эллиптического сегмента, погруженное в жидкость, и исследуем устойчивость его симметричного положения равновесия (рис. 1). Будем полагать, что полуоси эллипса равны a и b, а высота его сечения — H.
Отметим, что простейшие выражения для задачи о плавании судна с целико-вым эллиптическим сечением приведены в [18], где упоминается, что такая форма сечения не встречается в практике судостроения, хотя и удобна для теоретических расчетов, представляющих определенный интерес. В то же время поставленная нами задача об устойчивости плавания судна с сечением в форме эллиптического сегмента имеет уже существенно большее практическое значение, хотя ее решение оказывается более трудоемким.
Рис. 1. Судно с сечением в форме эллиптиче-
ского сегмента.
Приступая к его построению, запишем уравнение эллипса в осях £п, жестко связанных с эллипсом:
£2 п2
Будем использовать следующую параметризацию эллипса:
£ = a sin ф, п = —b cos ф, (2)
которая удобна для дальнейших действий. Здесь ф — параметр, причем нижней точке эллипса £ = 0, п = —b отвечает ф = 0. Положим, что значение ф* соответствует точке эллипса, в которой эллиптическая граница сечения бруса переходит в прямолинейную, т. е. когда п* = H — b. Тогда
Н — Ъ = — bcosф*, 7 = — = 1 — cosф* = 2sin2 -у, (3)
где y — безразмерный параметр задачи, 0 < y ^ 2, а 0 < ф* ^ п, причем случай Y = 2 отвечает целиковому эллиптическому сечению (ф* = п). Нетрудно понять, что площадь эллиптического сегмента Sb дается формулой
Sb = 2 j £dr] = 2ab j sin2 ф<1ф = у (2ф,_ - sin 2ф,.), (4)
-b о
тогда как положение его центра масс на оси п определяется выражением
п* „ „
2 ( , 2ab2 f ,2 , 4b sin3 ф„
1ъ = тг V&V =--/ sin фсовф(1ф = -——-. . 5
Sb J Sb J 3(2фФ - sin 2фФ)
-b о
Далее запишем уравнение ватерлинии в осях £п в отклоненном положении бруса на угол у:
П = - tg у£ - L, (6)
где L = L(y) — пока что не известная функция (рис. 1). Определим координаты точек пересечения ватерлинии (6) и эллипса (1). Объединяя эти уравнения, приходим к квадратному уравнению:
(a2 tg2 у + b2)£2 + 2a2Ltgу£ + (L2 - b2)a2 =0, (7)
решение которого имеет вид:
б,2 = aHg^ + b2 (-aLtgy±b^a2tg2y + b2-L2^) , (8)
где под £i будем подразумевать меньший корень (со знаком «—»), а под £2 — больший (со знаком «+»), как это показано на рис. 1. Следовательно, из (6) будем иметь
Vl'2 = ~a2tgV + &2 + . (9)
b
a2 tg2 у + b2
Тогда согласно (2) можно получить и соответствующие этим точкам значения ф\22, которые будут определяться следующими формулами:
cos ф]_ 2 = ~ТТТ~-Г? ( bL ± atg<p J a2 tg2 <р + Ъ2 - L2 ) ,
а2 tg2 ф + b2 \ V J (10)
sinV>i,2 = 9, }——To ( -aL tg <p ± bJa? tg2 <p + b2 - L2 ) .
a2 tg2 ф + b2 \ v у
Вычислим теперь вытесненную площадь жидкости используя параметрические уравнения (2) и формулы (10):
У 2
= У [— tg ф(а sin ф) — L + b cos -0] a cos —d-ф = V>i (11)
ab ( , гЬ^а2 tgV + Ь2 -L2\
= ТГ I r2 — rl---—ГЪ- •
2 \ a2 tg2 ф + b2 I
Обозначая Ф2 — —i = 2a, вычислим sin2a и cos2a с помощью формул (10):
• „ • , , • , , 2L\/a2 tg2 ip + b2 — L2 sin 2a = sin Ф2 cos фi — sin фi cos Ф2 = -
a2 tg2 ф + b2 (12)
2L2
cos 2a = cos ф 1 cos Ф2 + sin sin Ф2 = —ñ-г — 1.
a2 tg2 ф + b2
Учитывая теперь, что S¡ = pSb, приходим согласно (4), (11) и (12) к следующему соотношению:
2a — sin 2a = р (2фФ — sin2—*). (13)
Таким образом, что, хотя величины ф\ и Ф2, очевидно, зависят от угла у, величина а от у не зависит. Выражение (13) позволяет установить зависимость а от ф* при любом заданном значении р из допускаемого диапазона 0 < р < 1. При этом, очевидно, 0 < а < ф*. Что же касается зависимости Ь(у), то ее проще всего установить из второго выражения (12):
L2 = i(a2tgV + 62)(l + cos2a) = (a2 tg2 <р + Ь2) cos2 a. (14)
Вестник СПбГУ. Математика. Механика. Астрономия. 2024■ Т. 11 (69). Вып. 1 175
Перейдем далее к вычислению координат центра масс вытесненной площади жидкости в осях по формулам
-2
Si Фi
4а2 tg у (a2 tg2 у + Ъ2 — L2)3/2 4a2 sin3 a tg у
ф2
= — / ÇdS = — / a sin ф [— tg y(a sin ф) — L + 6 cos ф] a cos ipdip = Si J Si J
3(a2 tg у + 62)2(2a — sin 2a) 3(2a - sin 2a) ^Ja2 tg2 (p + fJ2''
ф2
J rjdS = I I ( — tg tp(as
l 4b2 (a2 tíф2 у + Ъ2 - L2)3/2 4Ъ2 sin3 a 1
ф2
Vi = тг 'idS = —— / [(-tgy(asin - Ь)2 - (-bcos ^)2] acos V'dV'
S^ 2Si J
3(а2 tg2 у + Ь2)2(2а — йш 2а) 3(2а-зш2а) у^Н^^Т^2'
(15)
где учтено выражение (14). Тогда вертикальная координата центра масс вытесненной площади относительно ватерлинии определится формулой
4 sin3 a cos у л/a2 tg2 у + Ь2
Vi = {L + 'i]i cosy + £csmy = Leos y--—-——-, 16
3(2a — sin 2a)
а соответствующая координата центра масс бруса
,т\ т 4b sin3 ф* cos w
уь = (L + щ cosy = Leos у- ——-. 17
3(2-0* — sin
Следовательно, потенциальная энергия плавающего бруса
П = Р( - ) = -Р eos I sin3 «У7«2 tg2 у + б2 Ъ sin3 ф, \
(Уь Vi) 3 cosy^ 2a-sin 2a 2ф,_ - sin 2ф,_ J ' 1 '
а ее квадратичная аппроксимация вблизи положения равновесия у = 0 с точностью до несущественной аддитивной постоянной примет вид
П = -РЪ
3
(S2 — 1) sin3 a sin3 ф* 2
2a — sin 2a 20* — sin 20*
у2, (19)
где введен еще один безразмерный параметр S = a/b. Из (19) вытекает условие устойчивости равновесного положения у = 0:
ô2 > 1 _ sin3 ф.,{ 2a-sin 2а) ^^
sin3 а(2ф* — sin2^*)
Отметим, что зависимость а от р и ф* дается формулой (13), а связь безразмерного параметра y с ф* определяется вторым выражением (3). Все сказанное позволяет заключить, что условие устойчивости (20) получено в терминах трех ключевых безразмерных параметров задачи р, S и y .
Рассмотрим теперь предельные случаи. Сначала обсудим вариант, когда брус практически не погружен в жидкость, т. е. величина р является исчезающе малой. Переходя в (20) к пределу при р ^ 0 и учитывая, что при этом а ^ 0 согласно
(13), а также принимая во внимание, что в таком случае (2а — sin 2а)/ sin3 а ^ 4/3, получим следующее условие:
¿2 > !__4sin3^_
3(2-0* — sin 2-0*)' 1 j
которое фактически является условием устойчивости равновесия бруса на горизонтальной плоскости. Кроме того, если рассматривается целиковый эллипс, т. е. ф* = п, то как условие (20) при 0 < а < п (т. е. 0 < р < 1), так и условие (21), полученное при р ^ 0, приводятся к одному и тому же неравенству 5 > 1. Оно означает, что для устойчивости как плавающего тела с эллиптическим сечением, так и аналогичного бруса на горизонтальной плоскости должно быть а > Ь. Что же касается другого частного случая, когда брус почти полностью погружен в жидкость, то для него значение р оказывается близким к единице. Переходя к пределу при р ^ 1 и учитывая, что при этом исходя из (13) будем иметь а ^ ф*, нетрудно получить, что неравенство (20) переходит в условие устойчивости 52 > 0, которое выполняется при любых значениях параметра 5, отвечающих смыслу задачи.
Также необходимо отметить, что полученные в настоящей задаче формулы для координат центров масс бруса и вытесненной площади могут быть получены и несколько иным путем без вычисления интегралов, если принять во внимание то, что эллипс с полуосями а и Ь может быть найден из круга радиуса Ь растяжением (при а > Ь) вдоль оси £ в а/Ь раз, а для кругового сегмента все интересующие геометрические параметры хорошо известны. Тем не менее в следующей задаче о плавании судна с сечением в форме гиперболического сегмента такую простую аналогию провести не удается, тогда как ход решения по смыслу оказывается тем же самым. Именно поэтому в задаче о судне с сечением в форме эллиптического сегмента были представлены все выкладки, позволяющие далее осуществить более наглядное сопоставление двух задач.
4. Устойчивость плавания судна с сечением в форме гиперболического сегмента. Рассмотрим теперь задачу об устойчивости плавания судна с сечением в форме гиперболического сегмента высотой Н (рис. 2).
Описывая гиперболический контур сечения в связанных с ним осях £п уравнением
п2 е
б2 а2 - 1, (22)
введем следующую параметризацию гиперболы:
£ = а эИ ф, п = Ь сИ ф, (23)
где ф — параметр, причем нижней точке гиперболы £ = 0, п = Ь отвечает значение ф = 0. Будем считать, что значение ф* > 0 соответствует точке гиперболы, в которой гиперболическая граница сечения бруса переходит в прямолинейную, т. е. когда п* = Н + Ь. Тогда
Н + Ъ = ЪсЪф*, 7 = у = сЪф,. - 1 = 2эЬ2 у, (24)
где 7 — безразмерный параметр задачи. Тогда площадь гиперболического сегмента Бь определится формулой
п* ф*
Яь = 2 У £<1г) = 2аЪ J вЬ2 ф<1ф = у (вЪ2ф,_ -2ф*), (25)
ь о
тогда как положение его центра масс на оси п дается выражением
п* ф* 3
2 [ м 2аЬ [ V2 / ь Ы, 45эЬ ф-,_
ьо
Уравнение ватерлинии в осях £п в отклоненном положении бруса на угол у в этой задаче записывается в виде
п = - tg у£ + Ь, (27)
где Ь = Ь(у) — функция, подлежащая определению (рис.2). Объединяя уравнения ватерлинии (27) и гиперболы (22) в систему, приходим к квадратному уравнению:
(Ь2 - а2 tg2 у)£2 + 2а2Ь tg у£ + а2(Ь2 - Ь2) = 0, (28)
из которого определяются координаты £12 точек их пересечения:
= 62-а2Ч2у • (29)
Как и прежде, под £1 будем подразумевать меньший корень (со знаком «—»), а под £2 — больший (со знаком «+»), как это приведено на рис.2. Пользуясь формулой (27), найдем отвечающие им координаты п1,2:
т'2 = Ь2-а2^2у • (30)
ъ
Ь2 — а2 tg2 у
Соответствующие этим точкам значения ф^ согласно (23) подчиняются выражениям
с11^1'2 = Та-п 2 (ь^То^ул/а2 tg2 у + Ь2 - Ъ2) ,
Ь — а tg у \ / (31)
вЬ-01,2 = Ъ2 _ д2 ^ ^-аЬ tg^f±Ъ\J а2 tg2 у + Ь2- Ъ2^ .
Перейдем теперь к нахождению площади вытесненной жидкости 5;, используя параметрические уравнения (23) и формулы (31):
Si = [— tg у (a sh ф) + L — b ch ф] a ch фё,ф
ab ~2
' 2L Va2 tg2 y+ L2-b2
~{Ф2 ~Фl) + "
(32)
b — a tg у
Обозначая ф2 — ф1 = 2а, вычислим sh2a и ch2a с помощью формул (31):
2 Ьу/ a 2 tg 2 у + L 2 — b2
sh 2a = sh ф2 ch ф1 — sh ф1 ch ф2
b2 - a2 tg2 y
2L2
ch 2a = ch ch 02 — sh 0i sh 02 = -^—„--1.
b2 — a2 tg2 у
(33)
А тогда из соотношения Si = pSb с учетом (25), (32) и (33) будем иметь
sh a — a = р (sh 2ф* — 2ф*). (34)
Применительно к этому выражению можно сделать тот же вывод, что и по отношению к (13): хотя величины ф1 и ф2 зависят от угла у, величина a от него не зависит, и она определяется при любом заданном значении р из диапазона от 0 < р < 1 по значению параметра ф*, при этом 0 < a < ф*. Наконец, зависимость L(y) устанавливается из второго выражения (33):
L2 = \{Ь2 - a2 tg2 у) (1 + ch2a) = (b2 - a2 tg2 y) ch2 a. (35)
Вычислим далее координаты центра масс вытесненной площади жидкости в
осях
Ф2
1 Г 1 Г
= — / ÇdS = — / a sh 0 [— tg <p(a sin ф) + L — b ch 0] a ch фdф = Si Si
Si Si
Si 'i
4a2 tg y(a2 tg2 y + L2 — b2)3/2 4a2 sh3 a tg y
3(62 — a2 tg2 y)2(sh2a — 2a) 3(sh2a-2a) ^52 _ a2 tg2 ^'
Ф2
(36)
4 = ^¡1VdS= À / [(_tg^(ash^) + ¿)2 " (^chV1)2] ach0d0
Si 'i
4b2(a2 tg2 y + L2 — b2)3/2 4b2 sh3 a 1
3(62 — a2 tg y)2(sh2a — 2a) 3(sh2a-2a) ^b2 - a2 tg2 y'
Вестник СПбГУ. Математика. Механика. Астрономия. 2024. Т. 11 (69). Вып. 1 179
2
где принята во внимание формула (35). Следовательно, вертикальная координата центра масс вытесненной площади относительно ватерлинии будет
4sh3acos у^Ъ2 -a2 tg2 у У1 = (т - L) cos у + io sm у = -L cos у +-3(gh2a_2a)-, (37)
тогда как соответствующая координата центра масс бруса есть
I т\ т 4bsh3 ф* cosу
yb = (щ - L) cos у = —L cos у + 3(sh2^_2vg • (38)
Таким образом, получаем выражение для потенциальной энергии плавающего бруса:
тт Pi ^ 4р ( sh3"Vtg2^
П = P(yb-yi) =-Pcosy \ ——-—---—--- , 39
3 \sh2^ — 2ф* sin 2а — 2а I
а ее квадратичная аппроксимация вблизи положения равновесия у = 0 с точностью до несущественной аддитивной постоянной есть
П= -РЪ
3
((52 +1) sh3 a sh3V, 2
у , (40)
sh 2а — 2а sh 2ф* — 2ф*
где также введен безразмерный параметр 6 = а/Ь. Из (40) устанавливаем условие устойчивости равновесного положения у = 0:
§2 > эЬ3 Ъ2а - 2а) _ 1 sh3 а^И2ф* — 2ф*)
Здесь зависимость а от р и ф* дается формулой (34), а связь 7 с ф* — вторым выражением (24). Поэтому можно считать, что условие устойчивости (41) получено в терминах трех безразмерных параметров: р, 6 и 7.
Совершенно ясно, что при рассмотрении предельного случая р ^ 0 мы вновь фактически будем иметь задачу об устойчивости равновесия бруса с сечением в форме гиперболического сегмента на горизонтальной плоскости, причем из (34) устанавливаем, что при этом также а ^ 0. Принимая теперь во внимание, что в условии устойчивости (41) будет в указанном случае — 2а)/sh3 а ^ 4/3, приведем его
к виду
¿2 > ___1 (42)
Наконец, остается осуществить предельный переход при р ^ 1 с учетом того, что а ^ ф* согласно (34). В результате получим, что условие устойчивости (41) приводится к виду 62 > 0 и оно выполняется при любых положительных значениях параметра 6.
Таким образом, из сопоставления формул (1)-(21) предыдущего раздела и формул (22)-(42) настоящего раздела видно, что все различия между ними заключаются в некоторых знаках и используемых в решении функциях (тригонометрических и гиперболических соответственно). Кроме того, в первой задаче 0 < ^ ^ 2, 0 < ф* ^ п, тогда как во второй задаче параметр 7 может принимать любые положительные значения, как и величина ф* .
5. Графическая иллюстрация результатов. Перейдем теперь к графической интерпретации полученных решений. Обращаясь сначала к задаче об устойчивости плавания судна с сечением в форме эллиптического сегмента, зададимся различными значениями параметра 7 и построим отвечающие им границы областей устойчивости согласно (20) и с учетом (3) и (13) на плоскости параметров рб (рис.3, а). Аналогичным образом поступим и применительно к задаче об устойчивости плавания судна с сечением в форме гиперболического сегмента, опираясь на условие устойчивости (41) и принимая во внимание формулы (24) и (34) (рис.3, б). Отметим, что в обеих задачах области устойчивости будут лежать выше построенных границ. Видно, что с увеличением значения 7 область устойчивости сужается, тогда как при фиксированном 7 с увеличением параметра р диапазон значений б, отвечающих устойчивому равновесию, расширяется. Наконец, построенные графические иллюстрации находятся в полном соответствии со всем сказанным ранее относительно частных случаев, когда величина р стремится либо к 0, либо к 1.
Р Р
а б
Рис. 3. Границы областей устойчивости плавания судна с сечением в форме: а — эллиптического сегмента, б — гиперболического сегмента.
6. Заключение. В настоящей статье были решены две задачи об устойчивости плавания судов с сечениями в форме эллиптического и гиперболического сегментов, которые являются достаточно хорошим приближением к реальным плавающим средствам. Для решения этих задач были использованы методы аналитической статики, сводящиеся к построению точного выражения для потенциальной энергии плавающего тела и его исследованию вблизи тривиального положения равновесия. Проделанные построения позволили получить условия устойчивости в каждой из задач в достаточно компактной форме с учетом введения безразмерных параметров. Была дана графическая иллюстрация этих условий путем построения отвечающих им областей устойчивости на плоскости двух параметров при различных значениях третьего параметра. Полученные результаты оказались интересными в теоретическом отношении, а также они имеют важное прикладное значение в корабельном деле. В заключение следует отметить, что полученные точные выражения для потенциальной энергии в рамках рассмотренных конфигураций (т. е. принятых вариантов погружения брусьев с расположением прямолинейных частей гра-
ниц их сечений выше ватерлинии) имеют существенно важное значение. В самом деле, они оказываются полезными в случае неустойчивости исследуемого симметричного положения равновесия, поскольку позволяют при необходимости получить посредством аналитического или численного анализа смежное с ним устойчивое равновесие, если оно отвечает той же конфигурации. Данное обстоятельство лишний раз свидетельствует о том, что избранный с самого начала статический подход к анализу устойчивости и конструирование точных формул для потенциальной энергии являются вполне оправданными, а сами эти формулы могут стать основой для последующих исследований.
Литература
1. Седов Л. И. Механика сплошной среды,. Москва, Наука (1952).
2. Борисов Р. В., Луговский В. В., Мирохин Б. В., Рождественский В. В. Статика корабля. Санкт-Петербург, Судостроение (2005).
3. Черданцев С. В. Теоретические основы расчета понтонов, используемых на угольных разрезах. Физико-технические проблемы разработки полезных ископаемых 1, 61—69 (2013).
4. Черданцев Н. В., Черданцев С. В. Остойчивость понтонов в зумпфах угольных разрезов на больших углах крена. Вестник Кузбасского государственного технического университета 4 (98), 32-37 (2013).
5. Скалина М.К., Смирнов А. С. Равновесие и устойчивость прямоугольной баржи со смещенным центром тяжести в жидкости. Неделя науки СПбПУ. Материалы научной конференции с международным участием. ИПММ. Санкт-Петербург, 18-23 ноября 2019 г., 95-98 (2019).
6. Скалина М.К., Смирнов А. С. Устойчивость положений равновесия плавающей прямоугольной баржи в условиях плоской задачи. Труды семинара «Компьютерные методы в механике сплошной среды» (2020-2021), 41-55 (2021).
7. Меркин Д. Р., Смольников Б. А. Прикладные задачи динамики твердого тела. Санкт-Петербург, Изд-во С.-Петерб. ун-та (2003).
8. Гарбарук А. В. Механика жидкости и газа. Задачи и решения. Ч. 1. Санкт-Петербург, Аркуш (2003).
9. Меркин Д. Р., Бауэр С.М., Смольников Б. А., Смирнов А. Л. Теория устойчивости в примерах и задачах. Москва; Ижевск, Регулярная и хаотическая динамика (2007).
10. Фан-дер-Флит А. П. Остойчивость прямоугольного понтона при наклонении от 0 до 90°. Известия Санкт-Петербургского политехнического института 10, 281-307 (1908).
11. Фан-дер-Флит А. П. Начальная остойчивость прямоугольных параллелепипедов. Известия Санкт-Петербургского политехнического института 10, 397-406 (1908).
12. Жуковский Н. Е. Полное собрание сочинений. Т. 9. Математика. Астрономия. Речи. Доклады. Характеристики и биографии. Москва, Ленинград, ОНТИ НКТП СССР (1937).
13. Смирнов А. С., Хашба Т. Н. Устойчивость плавания барж с трапециевидным и пятиугольным сечениями. Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия 9 (67), 740-752 (2022). https://doi.org/10.21638/spbu01.2022.415
14. Кравчинский И. А., Смирнов А. С. Устойчивость плавания судна с сечением в форме параболического сегмента. Неделя науки ФизМех. Сборник статей Всероссийской научной конференции. Санкт-Петербург, 4-9 апреля 2022 г., 202-205 (2022).
15. Войткунский Я. И. (ред.) Справочник по теории корабля: в 3 т. Т. 2. Статика судов. Качка судов. Ленинград, Судостроение (1985).
16. Валле Пуссен Ш.-Ж. Лекции по теоретической механике. Т. 2, пер. с фр., Москва, ИИЛ (1949).
17. Ржаницын А. Р. Устойчивость равновесия упругих систем. Москва, ГИТТЛ (1955).
18. Постон Т., Стюарт И. Теория катастроф и ее приложения, пер. с англ. Москва, Мир (1980).
Статья поступила в редакцию 17 февраля 2023 г.;
доработана 1 мая 2023 г.; рекомендована к печати 31 августа 2023 г.
Контактная информация:
Смирнов Алексей Сергеевич — ассистент, мл. науч. сотр.; smirnov.alexey.1994@gmail.com Кравчинский Игорь Андреевич — аспирант, инженер; igor.krav2000@gmail.com
Stability of floating vessels with cross sections in the form of elliptic and hyperbolic segments
A. S. Smirnov1'2, I. A. Kravchinsky2'3
1 Peter the Great St. Petersburg Polytechnic University,
29, ul. Polytechnicheskaya, St. Petersburg, 195251, Russian Federation
2 Institute for Problems in Mechanical Engineering of the Russian Academy of Sciences, 61, Bolshoy pr. V. O., St. Petersburg, 199178, Russian Federation
3 "CMT-Engineering",
28/2, Bolshoy Sampsonievskiy pr., St. Petersburg, 194044, Russian Federation
For citation: Smirnov A. S., Kravchinsky I. A. Stability of floating vessels with cross sections in the form of elliptic and hyperbolic segments. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy, 2024, vol. 11 (69), issue 1, pp. 171-184. https://doi.org/10.21638/spbu01.2024.112 (In Russian)
The paper considers two problems on the stability of the trivial equilibrium position of floating vessels with cross sections in the form of elliptic and hyperbolic segments. A review of examples on the stability of floating bodies is given and the key principles of its study by methods of analytical statics are outlined. For each of the presented problems, by means of rather serious mathematical constructions, an exact expression for the potential energy is obtained within the framework of the accepted configuration, and its quadratic approximation is calculated near the equilibrium state under study. On its basis, stability conditions of the equilibrium position are established in terms of three dimensionless parameters, and limit cases are also analyzed. A comparison of intermediate expressions and final results obtained in the process of discussing each of the problems is carried out, and their common features and distinctive features are revealed. The found solutions are illustrated as families of boundaries of stability regions on the plane of two dimensionless parameters for different values of the third parameter. These results are of fundamental theoretical importance and may be useful for practical applications.
Keywords: floating vessel, elliptic and hyperbolic segments, stability, static analysis, boundary of the stability region, plane of dimensionless parameters.
References
1. Sedov L.I. Mekhanika sploshnoi sredy. Moscow, Nauka Publ. (1952). (In Russian) [Eng. transl.: Sedov L.I. In: 2 vols. Singapore, World Scientific Publishing Co. (1997)].
2. Borisov R. V., Lugovsky V. V., Mirokhin B. V., Rozhdestvensky V. V. Ship statics. St Petersburg, Sudostroenie Publ. (2005). (In Russian)
3. Cherdantsev S. V. Theoretical basis for calculation of pontoons used in open pit mines. Fiziko-tekhnicheskie problemy razrabotki poleznykh iskopaemykh 1, 61—69 (2013). (In Russian) [Eng. transl.: Journal of Mining Science 49 (1), 52-59 (2013) https://doi.10.1134/S1062739149010073].
4. Cherdantsev N. V., Cherdantsev S. V. Stability of pontoon at coal pit under greater corners tumbling. Vestnik Kuzbasskogo gosudarstvennogo tekhnicheskogo universiteta 4 (98), 32-37 (2013). (In Russian)
5. Skalina M.K., Smirnov A. S. Equilibrium and stability of a rectangular barge with a displaced center of gravity in a liquid. Week of Science of SPbPU. Materials of scientific conference with international participation. IAMM. St Petersburg, November 18-23, 2019, 95-98 (2019). (In Russian)
6. Skalina M.K., Smirnov A. S. Stability of equilibrium positions of a floating rectangular barge under conditions of a flat problem. In: Computer methods in continuum mechanics. 2020-2021, 41—55 (2021). (In Russian)
7. Merkin D.R., Smolnikov B.A. Applied problems of the dynamics of a rigid body. St Petersburg, St Petersburg University Press (2003). (In Russian)
8. Garbaruk A. V. Mechanics of liquid and gas. Tasks and solutions. P. 1. St Petersburg, Arkush Publ. (2003). (In Russian)
9. Merkin D.R., Bauer S.M., Smirnov A. L., Smolnikov B.A. The theory of stability in examples and problems. Moscow, Izhevsk, R&C Dynamics (2007). (In Russian)
10. Fan-der-Flit A. P. Stability of a rectangular pontoon at an inclination from 0 to 90°. Izvestia of the St Petersburg Polytechnic Institute 10, 281—307 (1908). (In Russian)
11. Fan-der-Flit A. P. Initial stability of rectangular parallelepipeds. Izvestia of the St Petersburg Polytechnic Institute 10, 397-406 (1908). (In Russian)
12. Zhukovsky N. E. Complete set of works. Vol. 9. Mathematics. Astronomy. Speeches. Reports. Characteristics and biographies. Moscow, Leningrad, ONTI NKTP USSR Publ. (1937). (In Russian)
13. Smirnov A.S., Khashba T.N. Stability of floating barges with trapezoidal and pentagonal sections. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy 9 (67), 740752 (2022). https://doi.org/10.21638/spbu01.2022.415 (In Russian) [Eng. transl.: Vestnik St Petersburg University. Mathematics 55, iss.4, 504-512 (2022). https://doi.org/10.21638/spbu01.2022.415].
14. Kravchinsky I. A., Smirnov A. S. Stability of floating vessel with a cross section in the form of a parabolic segment. Week of Science PhysMech. Collection of articles of the All-Russian scientific conference. St. Petersburg, April 4-9, 2022, 202-205 (2022). (In Russian)
15. Voytkunsky Ya. I. (ed). Handbook of ship theory in three volumes. Vol. 2. Statics of courts. The pitching of ships. Leningrad, Sudostroenie Publ. (1985). (In Russian)
16. Ch.-J. de la Vallee Poussin. Leçons de mécanique analytique. T.II. (1925) [Rus. ed.: Ch.-J. de la Vallee Poussin. Lektsii po teoreticheskoi mekhanike. T. 2. Moscow, IIL Publ. (1949)].
17. Rzhanitsyn A.R. Stability of equilibrium of elastic systems. Moscow, GITTL Publ. (1955). (In Russian)
18. Poston T., Stewart I. Catastrophe theory and its applications. London, San Francisco, Melbourne, Pitman (1978) [Rus. ed.: Poston T., Stewart I. Teoriia katastrof i ee prilozheniia. Moscow, Mir Publ. (1980)].
Received: February 17, 2023 Revised: May 1, 2023 Accepted: August 31, 2023
Authors' information:
Aleksey S. Smirnov — smirnov.alexey.1994@gmail.com Igor A. Kravchinsky — igor.krav2000@gmail.com