Физическое моделирование движения судна в битых льдах Текст научной статьи по специальности «Строительство и архитектура»

ПРОБЛЕМЫ КОРАБЛЕСТРОЕНИЯ И ОКЕАНОТЕХНИКИ

УДК 629.124

В.А. Зуев, Е.М. Грамузов ФИЗИЧЕСКОЕ МОДЕЛИРОВАНИЕ ДВИЖЕНИЯ СУДНА В БИТЫХ ЛЬДАХ

Нижегородский государственный технический университет им. Р.Е. Алексеева

Рассмотрен уточненный метод пересчета сопротивления битого льда с модели на натуру с учетом плотности воды и льда.

Ключевые слова: ледовые условия, сопротивление, битые льды, судно, физико-механические характеристики.

Большое разнообразие ледовых условий, нестабильность свойств и характеристик битых льдов делают практически невозможным достаточно точно и определенно зафиксировать многочисленные параметры, характеризующие натурную ледовую обстановку. Однако это обстоятельство не снимает актуальности вопроса оценки ледовых качеств судов в битых льдах. По-прежнему судостроителей волнует вопрос ледового сопротивления в этих условиях. Для описания ледовой обстановки и ее связи с ледовым сопротивлением вводят понятие «стандартных ледовых условий», позволяющих в некоторой степени однозначно зафиксировать ледовые условия относительно наиболее важных параметров, влияющих на ледовое сопротивление. К ним следует отнести природные мелкобитые льды и битые льды в каналах, образованных после похода ледокола в сплошном ледяном покрове.

При регламентации условий ледового плавания Российским морским регистром судоходства используются следующие определения [4]:

• сплоченность - мера сплошности ледяного покрова, характеризуемая отношением площади, занимаемой льдами, к общей площади рассматриваемого участка акватории (оценивается по 10-балльной шкале);

• разреженный лед - лед сплоченностью 4-6 баллов, в котором большинство льдин не соприкасаются между собой;

• сплоченный лед - лед сплоченностью 7-8 баллов, в котором большинство льдин соприкасаются между собой;

• очень сплоченный лед - лед, сплоченность которого равна или более 9 баллов, но менее 10 баллов;

• сплошной лед - лед сплоченностью 10 баллов. Лед также характеризуется толщиной.

Российский речной регистр характеризует битые льды только толщиной [4]. Очевидно, что такое описание ледовых условий нельзя признать достаточным для адекватной оценки ледовой ходкости.

Мелкобитые льды принято характеризовать: толщиной h - средняя толщина плавающих льдин, как правило, равная толщине разрушенного ледяного покрова за счет ветра, волнения или после прохода ледокола; сплоченностью S - отношение площади, покрытой льдом, ко всей площади поверхности участка водоема. При анализе ледового сопротивления

© Зуев В.А., Грамузов Е.М., 2011.

к ним присоединяют: плотность льда рл - масса льдины, отнесенная к ее объему; коэффициент трения / пары лед - наружная обшивка корпуса судна. Для описания ледовых каналов к характеристикам ледовых условий добавляют ширину канала ВК. Заметного разрушения льдин при контактах между собой и корпусом не происходит, поэтому прочие физико-механические характеристики льда принято не включать в рассмотрение.

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

Устоявшимся считается разделение ледового сопротивления ЯЛ на чистое ледовое сопротивление ЯЧЛЛ и сопротивление воды ЯВ: ЯЛ= ЯЧЛ+ ЯВ. При этом чистое ледовое сопротивление разделяют на независящее от скорости V, так называемое прямое сопротивление ЯП, и скоростное сопротивление Яск: КЧЛ= Кп+ Кск. Сопротивление воды принято определять обычными методами и считать его независимым от чистого ледового сопротивления ^ЧЛ.

С физической точки зрения прямое сопротивление связывают с затратами энергии на притапливание льдин корпусом судна, а скоростное сопротивление - с динамическими процессами страгивания льдин при контакте с корпусом и их раздвигания в воде корпусом судна. Во всех случаях контакта льдин с корпусом возникают сопутствующие силы трения.

Картина взаимодействия корпуса судна с битым льдом существенно зависит от его сплоченности.

Так, например, для льда сплоченностью 10 баллов взаимодействие можно представить следующим образом. Льдины, взаимодействуя с корпусом, поворачиваются, притапливаются и раздвигаются, следуя обводам подводной части корпуса ледокола. Траектория движения льдин существенно зависит от формы корпуса. Практически весь сломанный лед попадает под корпус, образуя в его подводной части ледяную «рубашку». Облегание притопленным льдом подводной части корпуса происходит на значительном расстоянии от форштевня по длине судна и вниз от действующей ватерлинии. Ширина контакта обшивки с льдинами по периметру шпангоута с каждого борта равна, примерно, полуширине ледовой ватерлинии в соответствующем районе по длине судна (рис. 1). Между зонами облегания льда правого и левого борта образуется разрыв в ледяной рубашке. В носовой части судна обломки льда образуют оболочку, плотно прилегающую к наружной обшивке. Лед в данном районе корпуса прижимается к обшивке силами его плавучести и гидродинамическими силами, поскольку внешняя нормаль к наружной обшивке составляет острый угол с направлением движения. В кормовой части плотность оболочки уменьшается за счет противоположного действия гидродинамических сил. При приближении к корме возможны отрывы ледяной оболочки. Устойчивость ледяной оболочки особенно в кормовой оконечности зависит от формы корпуса подводной части и скорости движения судна. Пройдя наиболее широкую часть корпуса судна, часть льда всплывает по бортам, а часть остается под днищем и попадает в винты. После прохождения судна практически весь лед всплывает в канале, засоряя его до 9-10 баллов.

S;

Рис. 1. Схема облегания корпуса обломками льда:

1 - нижняя граница ледяной «рубашки»; 2 - траектория движения льдины)

При движении судна в разреженных битых льдах льдины поворачиваются, притапли-ваются и раздвигаются. Взаимодействие с корпусом, за счет меньшего стеснения льдин, происходит в районе действующей ватерлинии, лед практически не попадает под корпус. При промежуточных значениях сплоченности 6-9 баллов имеет место переходная картина взаимодействия, существенно зависящая от скорости хода.

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

Суть проведения испытаний, которая состоит в использовании в качестве моделей льда материалов, удовлетворяющих критериям подобия.

Анализ природы действующих на судно сил и условий их подобия позволил установить, что для моделирования движения в битых льдах необходимо выполнить следующие условия:

• форма корпуса модели должна быть геометрически подобной натурному судну. Масштаб геометрического подобия X;

• число Фруда одинаково для модели и натуры Fr=idem;

• сплоченность льда одинакова для модели и натуры S=idem;

• плотность воды и льда одинаковы для модели и натуры:

р5 = idem, рл = idem;

• коэффициент трения льда о корпус судна одинаков для модели и натуры f= idem;

При выполнении этих условий ледовое сопротивление натурного судна при движении

в битых льдах можно определить пересчетом в соответствии с законом подобия Фруда по кубу масштаба X3.

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

Геометрическое подобие модели может быть выполнено и проверено с достаточной точностью.

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

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

Соответствие плотностей воды и льда натурным условиям должны быть подвергнуты дополнительному анализу, который удобно выполнить для наиболее регулярных условий битого льда 10-балльной сплоченности.

Сопротивление битого льда можно представить в виде

Ro = Rn + Ru + R, (1)

где Rn - сопротивление, обусловленное плавучестью льдин;

Ru - инерционная составляющая сопротивления;

R - гидродинамическая часть сопротивления, связанная с сопротивлением воды движению льдин.

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

Применяя принцип суперпозиции, можно записать

dN = dN п + dN u + dN г. (2)

Если считать, что трение обшивки о лед подчиняется закону Кулона-Амонтона, то

dF = fdN, dF„ur = fdNnur. (3)

Рассмотрим последовательно каждую составляющую сопротивления льда.

Обозначим углы, характеризующие форму обводов корпуса некоторой точки судовой поверхности:

Ф1 - угол между касательной к батоксу и ОП;

Ф 2 - угол между касательной к ватерлинии и ДП;

ф3 - угол наклона шпангоута к вертикали.

Введем также направляющие косинусы нормали и некоторого касательного направления для произвольной точки судовой поверхности. Для краткости записи обозначим их:

nx = cos (n, x );ny = cos (n, y );nz = cos(n, z);

x x = cos (x, x );x y = cos (x, y );X z = cos (x, z ).

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

1Бф2 1

n = , n„ = '

x

д/l + tg^ + tg2ф3 ' у д/l + 1§2ф2 + 1ё2фз '

1§фз 2 , 2 2 Л . 1§ф2

3 • nx + nv + nx = 1; г§ф1 =- 2

"г I ; ~г, 5 "ж 1 'У ' Ч

V1 + гв2Ф2 + ^§2фз 1§Фз

Нормальное усилие, действующее на элемент обшивки йП, обусловленное плавучестью льда

dNn = п2йХп = пг (р - рл ^НйП. (4)

Проекция этих усилий dNn и dFn на направление движения (ось Ох):

Х = dNnnx + dFnхж . (5)

Интегрируя (5) с учетом (4) и (3) по всей зоне облегания корпуса льдом а , получим выражение для составляющей сопротивления, связанной с плавучестью льда

К = (р - р л \(ппх + А X х )йП. (6)

п л

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

Изменение кинетической энергии элемента льда в момент контакта под действием элементарной силы, приложенной со стороны обшивки,

dNu йц = К2 + Кх, (7)

где dNu - нормальная инерционная сила, действующая на элемент ватерлинии; dц - перемещение вдоль нормали;

К + К - изменение кинетической энергии элемента льда.

Деля уравнение (7) на йЦ и подставляя выражение для кинетической энергии, получим

ы Си Н- V)

dNu =— • 2-и., (8)

йц 2

где т - масса элемента льда;

Ущ и Ущ - скорости элемента льда до и после контакта;

си - безразмерный коэффициент, учитывающий присоединенные массы воды.

Пусть приращение скорости льда происходит скачкообразно в районе КВЛ от нуля до

скорости, обусловленной скоростью судна, тогда

^ =0; = ™х • (9)

Массу элемента льда представим в виде

т = рл ИЫхЫу. (10)

Выразим нормальное перемещение обшивки Ыц через поступательное перемещение судна Ых:

Ыц = Ыхпх. (11)

Подставляя (9), (10) и (11) в (8), получим

р V2

^ = 0иН^-пхйу • (12)

Проецируя инерционные и сопутствующие силы трения на направление движения и интегрируя по ширине носовой ветви КВЛ,

В/2

К =р л^2 / Си К2 + /ПхТ х Ыу • (13)

0

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

2 2 = сг = сг р^пх2ЫО• (14)

Касательную компоненту и суммарную проекцию гидродинамических сил, действующих на элемент льда, найдем по формулам (3) и (5).

Интегрируя по площади зоны контакта, определим

■"г Г";2 / Сг(Пх3 + /Пх2Т х] О л

В интегралы (6), (13) и (15) входит косинус угла между касательной к траектории движения льда в данной точке и направлением движения Т х . Из наблюдений за движением льда относительно обшивки следует, что лед движется в направлении, составляющем минимальный угол с направлением движения какой только допускают обводы корпуса в данном районе. В частности, по плоскому днищу движение происходит приблизительно по батоксам, а по вертикальному борту по ватерлиниям.

Определим Т х, используя направляющие косинусы внешней нормали. В силу орто-

Кг =PV2 { Сг(Пх3 + /Пх2Т х)ЫО • (15)

гональности нормали п и вектора Т , совпадающего по направлению с движением льда Кроме того,

ПхТ х + ПуТ у + ПТ Г = 0 • (16)

Т х2 + Т у2 + Т г 2 = 1 • (17)

Подставляя Т х, полученный из (16) в (17), получим уравнение, связывающее в неявной форме углы (Т , х) и (Т , у), определяющие направление касательной к траектории льдин в данной точке

х)'(Т' х)] = Тх2 +Т 2 +

Пх Пу

Тх + ^Ту

2

-1 = 0 • (18)

У

В силу гладкости функции (т, х) = /(т, х) при повороте вектора Т на касательной плоскости ее экстремум однозначно определится условием

_ а^а(т, y)

d(ï,y) dF/d(z,x) '

которое в развернутом виде приводит к уравнению

х 2 + -У Xy -

—х + — х

vx "у

V — — У

_ 0. (20)

Совместное решение уравнений (20) и (18) определяет искомый x x зависимости от направляющих косинусов внешней нормали

x x = ny2 + nz2 = ±^ 1 - nx2 . (21)

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

R = (p-p л )F (h, nx, nz, x x, f, Пл )+p л v2 F (h, nx, x x, f, 5) + pv2 F3 (nx, x x, f, ). (22)

Очевидно, что подобие функций F^ F2, F3 определяется геометрическим подобием и равенством коэффициентов трения модели и натуры. Соответствие плотностей воды и льда и их влияние на прогнозирование сопротивления натуры является предметом настоящего анализа.

При испытаниях моделей судов в битых льдах в настоящее время используются обычные (не ледовые) опытовые бассейны [1, 2, 3], наполненные пресной водой. Как показали исследования [1], наиболее подходящим материалом для имитации битого льда, получившим распространение в большинстве опытовых бассейнов, в настоящее время является полиэтилен высокого давления. Плотность этого материала 0,92 т/м3. Таким образом, плотность воды и льда с большой точностью соответствует льду пресных водоемов. Для морского льда необходимо вводить коррективы. Во-первых, морская вода имеет большую плотность, а во-вторых, морской лед - меньшую плотность. Огромное противоречие прогнозов ледового сопротивления, по нашему мнению, является результатом игнорирования этих фактов. При этом нетрудно заметить, что составляющие прямого и скоростного сопротивления подлежат непосредственному экспериментальному разделению. Этот факт позволяет ввести коррективы в схему пересчета ледового сопротивления с учетом фактических плотностей воды и льда.

Прямое сопротивление пропорционально разности плотностей воды и льда (p-p^). При плотностях морской воды и льда 1,025 т/м3, 0,89 т/м3 и пресноводных 1,000 т/м3, 0,92 т/м3 эта разница равна 0,175 т/м3 и 0,080 т/м3 соответственно. То есть прямое сопротивление при пересчете на морской лед необходимо увеличивать в 1,7 раза.

Скоростное сопротивление, разделяемое на инерционное и гидродинамическое, не позволяет экспериментальное разделение. Как видно из формулы (22), инерционное сопротивление пропорционально плотности льда, а гидродинамическое - плотности воды. Учитывая большую определенность плотности воды, по аналогии с гипотезой Фруда для остаточного сопротивления, будем полагать (до получения более достоверных данных) скоростное сопротивление пропорциональным плотности воды.

В заключение приведем пример пересчета ледового сопротивления судна-снабженца проекта 22450 ледового класса Arc 5. Модель в масштабе 1:60 испытывалась в плитках полиэтилена толщиной 20 мм при сплоченности 10 баллов в канале шириной 0,48 м.

Последовательность пересчета ледового сопротивления с модели на натуру выглядит следующим образом:

1. Скорость модели и м.

2. Сопротивление модели RM.

3. Число Фруда.

4. Скорость судна он.

5. Сопротивление воды модели Д^.

6. Чистое ледовое сопротивление модели Д^ = Дм - Д^.

7. Прямое ледовое сопротивление модели (и ^ 0)

8. Скоростное сопротивление модели = д^ - д^ .

9. Прямое сопротивление судна Д = Д^

о р_-рл

р -р.

х3.

л

с

10. Скоростное ледовое сопротивление судна дн = дм ^ я3.

ск ск п

р

11. Сопротивление воды судна Д.

12. Полное сопротивление судна Д = Д^ + Дснк + Д.

Индексами «н» и «м» обозначены натура и модель, а индексами «п» и «с» пресный и соленый лед соответственно.

Результаты пересчета ледового сопротивления морского судна без учета поправок на плотность воды и льда и с учетом соответствующих поправок приведены на рис. 2.

3000 2500 2000 1500 1000 500 0

9 > ,

г»

4

/

*

У г .

к*

т-

■т

Р2 Р3

v,м/с

0

1

2

3

4

Рис. 2. Сопротивление натуры:

Я1 - без учета влияния плотности воды и льда;

Я2 - с учетом влияния плотности воды; Я3 - с учетом влияния плотности воды и льда

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

Библиографический список

1. Зуев, В.А. Средства продления навигации на внутренних водных путях / В.А. Зуев. - Л.: Судостроение, 1986. - 208 с.

2. Рывлин, А.Я. Испытания судов во льдах / А.Я. Рывлин, Д.Е. Хейсин. - Л.: Судостроение, 1980. - 207 с.

3. Сазонов, К.Е. Ледовая управляемость судов / К.Е. Сазонов. - СПб.: ЦНИИ им. академика А Н. Крылова, 2006. - 251 с.

4. Российский Морской Регистр Судоходства. Правила. - СПб., 2008. Т. I.

Дата поступления в редакцию 01.02.2011

V.A. Zuev, E.M. Gramuzov PHYSICAL MODELING OF VESSEL MOTION IN BROKEN ICE

Refined methods of recalculation the resistance broken ice from the model to nature are considered. These methods are taking into account the density of water and ice.

Key words: ice conditions, resistance, broken ice, vessel, physical and mechanical characteristics.