Научная статья на тему 'Исследование поведения коэффициентов взаимодействия системы корпус-движитель по данным натурных испытаний речных грузовых судов'

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

CC BY
12
0
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
коэффициенты взаимодействия / коэффициент попутного потока / коэффициент засасывания / натурные испытания / теория взаимодействия / судовые тяговые расчѐты / нормирование времени следования и расхода топлива / ship propulsion / wake fraction / thrust deduction / full-scale trials of ships / ship hydrodynamics / rate of ship velocity and fuel consumption

Аннотация научной статьи по физике, автор научной работы — Платов Александр Юрьевич

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Платов Александр Юрьевич

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Study of the ship-propeller interaction by trials of the river cargo ships

The calculation of the ship interaction is necessary for prediction of ships power, which are used in the problems of computation of the ships velocity and fuel consumption, as well as for the preliminary design of new ships. In the scientific and educational literature there are descriptions of several empirical methods for calculating wake fraction and thrust deduction. However, model trials data are extremely scarce, and ships trials data are completely absent. Based on the existing data of full-scale trials of ships, ship interaction coefficients were calculated and these data were compared with the available empirical methods, as well as with the hydrodynamic theory of interaction. It is shown that empirical methods do not correspond to the data of full-scale trials, even at a qualitative level. At the same time, it is shown that the data of full-scale trials qualitatively correspond to the theory of interaction.

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

УДК 623.828:532.5 DOI: 10.37890/jwt.vi76.402

Исследование поведения коэффициентов взаимодействия системы корпус-движитель по данным натурных испытаний речных грузовых судов

А.Ю. Платов

https://orcid.org/0000-0002-4589-0348

Нижегородский архитектурно-строительный университет, г. Н.Новгород, Россия

Аннотация. Вычисление коэффициентов взаимодействия необходимо для проведения судовых тяговых расчётов, которые применяются при решении задач нормирования времени следования и расхода топлива, а также для эксплуатационного обоснования новых судов. В научной и учебной литературе можно найти описание нескольких эмпирических методов для расчёта коэффициентов взаимодействия, однако данные модельных испытаний представлены крайне скудно, а данные натурных испытаний полностью отсутствуют. На основе существующих данных натурных теплотехнических испытаний судов были проведены расчёты коэффициентов взаимодействия и проведено сравнение этих данных с имеющимися методами расчёта, а также положениями гидродинамической теории взаимодействия. Показано, что эмпирические методы расчёта не соответствуют данным натурных испытаний даже на качественном уровне. Вместе с тем показано, что данные натурных испытаний качественно соответствуют положениям теории взаимодействия.

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

Study of the ship-propeller interaction by trials of the river cargo

ships

Alexander J. Platov

https://orcid.org/0000-0002-4589-0348

Nizhny Novgorod State University of Architecture and Civil Engineering, Nizhny Novgorod,

Russia

Abstract. The calculation of the ship interaction is necessary for prediction of ships power, which are used in the problems of computation of the ships velocity and fuel consumption, as well as for the preliminary design of new ships. In the scientific and educational literature there are descriptions of several empirical methods for calculating wake fraction and thrust deduction. However, model trials data are extremely scarce, and ships trials data are completely absent. Based on the existing data of full-scale trials of ships, ship interaction coefficients were calculated and these data were compared with the available empirical methods, as well as with the hydrodynamic theory of interaction. It is shown that empirical methods do not correspond to the data of full-scale trials, even at a qualitative level. At the same time, it is shown that the data of full-scale trials qualitatively correspond to the theory of interaction.

Keywords: ship propulsion, wake fraction, thrust deduction, full-scale trials of ships, ship hydrodynamics, rate of ship velocity and fuel consumption

Введение

Определение коэффициентов взаимодействия (КВ) системы корпус-движитель (СКД) является одной из важнейших задач гидродинамики судна и имеет большое практическое значение для прогнозирования ходовых характеристик грузовых судов.

Для определения КВ СКД по методу Фруда, разработанному в середине 19-го века, проводятся модельные испытания судов в гидродинамических бассейнах. Развитие вычислительной техники сделало возможным с начала 2000-х годов практическое применение методов вычислительной гидродинамики для оптимизации ходовых характеристик. Для решения задач предварительного проектирования и технического обоснования судов также используются упрощённые эмпирические методы вычисления КВ СКД, основанные, как правило, на модельных испытаниях в гидродинамических бассейнах.

Кроме указанных выше задач вычисление КВ СКД требуется для расчётов винтовых и ограничительных характеристик при нормировании времени следования, а также для рейсового планирования работы грузовых судов. Для таких задач применимы только эмпирические (инженерные) методы вычисления КВ СКД, причём конечные результаты таких вычислений должны сопоставляться с теми характеристиками режима движения грузовых судов, которые требуются экипажам судов: частотой вращения гребных винтов (ГВ), скоростью судна относительно берега и эффективной мощностью главных двигателей (ГВ). Многолетний практический опыт автора, связанный с нормированием скоростей грузовых судов, говорит о высокой погрешности инженерных методов тяговых расчётов при определении упомянутых характеристик, возникающей из-за погрешностей вычисления сопротивления суда, гидродинамических характеристик ГВ, а также КВ СКД.

Между тем в научной и учебной литературе отсутствуют данные натурных испытаний по измерению КВ СКД, вероятно, в силу нетривиальности подобных испытаний. Такое положение дел мешает как оценке адекватности существующих эмпирических методов, так и возможному совершенствованию инженерных методов расчёта ходкости.

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

Существующие методы расчёта коэффициентов взаимодействия

В классической схеме для оценки ходовых качеств судна вводится пропульсивный коэффициент г/, который определяется по выражению [1]:

1 йг?

г=Тр2^П = гРг к1 , (1)

где гр - число движителей; Д - сопротивление воды движению судна; М' -момент на валу ГВ при работе за корпусом судна; п - частота вращения ГВ; гр -к.п.д. ГВ в свободной воде; - коэффициент влияния корпуса, - коэффициент неравномерности поля скоростей.

Коэффициенты в (1) определяются по выражениям:

— Рт?р . _ Ку . • _ ПЛ гр ~ 2пМп ' гк ~ грР'ур' 1 ~ М' Р' ( 2

где и - упор ГВ за корпусом судна и в свободной воде соответственно; -момент на валу ГВ в свободной воде; ур - скорость потока через сечение ГВ.

Выражения коэффициентов через безразмерные характеристики имеют вид: _ _ к1(лр )лр ; _ _ ^ ч

Ъ-ЪЪ)2п ' Ч1 - И/ 1_1 2 (3)

где ^ (Яр) — Р/(рВ4п2) - коэффициент упора ГВ в свободной воде; К2(Яр) _ М/(рВ5 п2) - коэффициент момента ГВ в свободной воде; Яр _ рр/пВ -относительная поступь ГВ; В - диаметр ГВ; и - коэффициента попутного потока; £ -коэффициента засасывания; и -

коэффициенты неравномерности поля скоростей по упору и по моменту соответственно; Я _ V/пВ - видимая (кажущаяся) относительная поступь; К' (Я) и К'(Я) - коэффициенты упора и момента при работе ГВ за корпусом судна, рассчитываемые аналогично и К2, но для значений Р' и М' соответственно.

Взаимодействие СКД описывается в этой схеме с помощью коэффициентов и, ¿1 и ¿2-

Коэффициенты взаимодействия были введены в 1883 г. Р. Фрудом [2], Д. Тейлор в [3] предложил иное определение и, а законченный вид схема (1)-(3) приобрела уже в середине 30-х годах 20 века в цикле статей Э.Э. Папмеля [4].

Все функции в (1)-(3) определяются с помощью модельных испытаний и пересчитываются на натуру. Принципиальный порядок определения КВ следующий.

Отдельно для ГВ определяются функции К1 (Яр) и К2 (Яр), а также отдельно для корпуса функция сопротивления . Далее на испытаниях модели самоходом в жёсткой упряжке определяются функции К' (Я) и К' (Я) .

После чего коэффициент попутного потока рассчитывается из

решения уравнения

К1(Яр) _ К'(Я). (4)

При таком подходе, как можно видеть, * 1 _ 1. Иначе говоря, коэффициент попутного потока определяется из равенства .

Коэффициент неравномерности определяется затем из уравнения

*2 _ (5)

2 К2^р) ( )

Коэффициент засасывания рассчитывается из выражения

Я Ке{Х) _ Ке{Х) _ Я

_ К[(Х)_ 1~1^р)_1 ~Р' (6)

где .

Можно заметить, что в методе Э.Э. Папмеля коэффициент и определялся через отношение поступей и при нулевом упоре ГВ за корпусом судна и в свободной воде и принимался постоянным. Поэтому величина была отличной от 1 и переменной.

Очевидно, что в описанной схеме можно изменить определение коэффициента и вычислять его, исходя из соотношений моментов. Впервые такой вариант вычисления был предложен в 1937 г. Ф. Хорном на 4-ой Международной конференции опытовых бассейнов (МКОБ) [5], который показал, что при использовании равенства моментов получаются другие результаты.

Поэтому в настоящее время различают два коэффициента попутного потока: коэффициент по упору , вычисляемый из (4), и коэффициент по моменту , который вычисляется из решения следующего уравнения

К2(Яр) _ К'(Я). (7)

При таком подходе получается, что 12 = 1 , а ^ будет вычисляться из соотношений коэффициентов упора. На практике модельных испытаний, как правило, применяется метод с коэффициентом , в то время как натурные испытания дают возможность оперировать с шМ. Как утверждается в [6], разница между двумя описанными способами незначительная, однако судя по данным испытаний модели морского грузового судна, приведённым в материалах 9-ой МКОБ [7], разница между шР и шМ может достигать 11%.

На базе модельных испытаний, проведённых по вышеописанной принципиальной схеме, в разное время были предложены эмпирические зависимости для вычисления КВ. Метод Э.Э. Папмеля, предназначенный для морских судов, был модифицирован для речных судов в ЦНИИРФ (1951) [8]. Кроме этого, в разное время для расчёта коэффициентов № и t были предложены методы Е.И. Степанюка (1967) [9], М.П. Зорина (1977) [10] и В.А. Лесюкова (1982) [11]. Отечественные методы разработаны в вариантах как для открытых ГВ, так и для ГВ в насадках.

Методы для речных судов автору из иностранной литературы неизвестны. Для морских судов на основе испытаний в 1949-1950х гг. был разработан метод С. Харвальда [12] (1983), в современном варианте излагаемый в [13]. Наиболее часто в современной, в том числе и отечественной литературе рекомендуется метод Холтропа-Меннена, который был разработан в 1977 г. тоже для морских судов и усовершенствован в 1984 г. [14].

Материалы и методы исследования

В настоящем исследовании анализ поведения КВ проводился по данным натурных испытаний из [15; 16]. Эти испытания для грузовых судов были выполнены в конце 70-х-начале 80-х годов под руководством начальников теплопартий ВОРПа Е.М. Тумаринсона и «Волготанкера» - А.И. Есина. Кроме этого, привлекались неопубликованные данные испытаний судов новых проектов, выполненных в начале 1990-х сотрудником ГИИВТа В.В. Бажанкиным и в 2010-х годах начальником теплопартии Волжского пароходства С.Г. Смирновым. Всего из указанных источников автором были собраны данные по 98 судам разных проектов.

В результате теплотехнических и ходовых натурных испытаний измеряются винтовые характеристики для скорости судна и эффективной мощности главных двигателей (ГД) при разных значениях частоты вращения вала ГД. Этих данных достаточно, чтобы вычислить коэффициент попутного потока в зависимости от скорости судна или частоты вращения.

Для этого требуется использование какого-нибудь метода для вычисления коэффициентов и . В настоящем исследовании использовалась известная аппроксимация кривых действий по методу Х.В. Лербса (1975) [17].

Вычисление коэффициента засасывания может быть получено только при использовании какого-нибудь метода расчёта сопротивления, который может иметь в общем случае большую погрешность. Для данного исследования применялась аппроксимация остаточного сопротивления по формуле [18]:

Иц = 2 59,1 (1 + 114,8 (Рг - 0 , 1 3 9) 2) (В/ Т/8,42 5) 06 82 дУРг3 , (8)

где - число Фруда по длине корпуса, - ширина корпуса, - осадка, -водоизмещение, - ускорение свободного падения. Формула (9) получена по данным модельных испытаний речных грузовых судов с коэффициентом полноты водоизмещения 5 > 0 ,8.

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

Корреляция между натурными данными и методами расчёта коэффициентов

взаимодействия

Первый вопрос исследования: насколько верно эмпирические методы отражают зависимость КВ от характеристик судна. Для ответа на этот вопрос были определены коэффициенты корреляции по Пирсону. За х принято значение, вычисленное по эмпирическим методам (расч.), за у - результат расчёта по натурным данным (эксп.), то есть рассчитывалась теснота связи зависимости и/Мэксп. = / (1иМрасч) .

Для судов без насадок (всего 37 судов) был рассчитан им по натурным данным, который сравнивался с расчётным иР, определённым по методам ЦНИИРФ, Лесюкова, Харвальда и Холтропа. Коэффициент корреляции Пирсона для этих четырёх методов приведён в табл. 1.

Таблица 1

Значения коэффициента корреляции по Пирсону

№ Метод ^Мэксп . _ /(^Мрасч .) ^эксп . f(^расчJ

1 ЦНИИРФ -0,016 0,005

2 Лесюкова -0,365 -0,208

3 Харвальда -0,247 0,048

4 Холтропа -0,054 -0,141

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

зависимости £:эксп . = / ( £расч ) будет точно такой же, как в табл. 1. Вместо этого был рассчитан коэффициент корреляции между натурным значением и расчётным, но вычисленным по эмпирическим формулам, в которых подставлялось натурное значение п . Вычисленные коэффициенты Пирсона приведены также в табл. 1. Нужно заметить, что при вычислении £ предполагалось, что ^(Яр) = (Я) . На рис. 1 можно видеть значения расчётных по методу ЦНИИРФ коэффициентов в сравнении с натурными и полное отсутствие зависимости между этими величинами.

Рис. 1. Отклонение расчётного значения коэффициента попутного потока от данных натурных испытаний

Можно также видеть, что один из коэффициентов имеет отрицательное значение (-0,01 для судна проекта №795). Расчётное значение, напротив, велико (0,461).

Зависимость коэффициентов взаимодействия от скорости движения судна

Второй вопрос исследования: какова вариация значений КВ на винтовых характеристиках.

Во всех вышеперечисленных эмпирических методах расчёта КВ коэффициент попутного потока считается постоянным на обычных режимах движения судов. Это прямо следует из формул для вычисления иР. Например, формула метода ЦНИИРФ имеет вид:

иР _ 0,11 + 0,0852+ АиР, (9)

где величина растёт со скоростью судна, но равна нулю для скоростей при числе Фруда меньше 0,2. А под эту категорию попадают все обычные речные грузовые суда. Сложная аппроксимация Холтропа [14] для судов с двумя винтами вообще не содержит зависимости от скорости.

Согласно модельным испытаниям, проведённым Харвальдом в 70-е годы [12], изменение и Р в диапазоне чисел Фруда 0,15-0,24 достигает 25%, причём и Р уменьшается с ростом числа Фруда. При этом испытывались одновинтовые модели с коэффициентом полноты 5 _ 0,75 5 — 0,805 .

В исследовании В.Граффа (1961) [19] испытывалась модель двухвинтового речного судна типа «Густав Кёнигс» с коэффициентом полноты 5 _ 0,866 при разных относительных глубинах. При наибольшей относительной глубине ( )

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

В исследовании Г. Лутра (1976) [20] исследовались три модели речных двухвинтовых судов с коэффициентом полноты 87 87 . Качественные

результаты совпадают с результатами В. Граффа.

Подобные результаты можно видеть также в исследовании [21], в котором испытывались модели мелкосидящих судов с туннельными обводами кормы.

Расчёт коэффициента на основе натурных испытаний находится в полном противоречии с перечисленными выше результатами. Это можно видеть из рис. 2.

Рис. 2. Зависимость КВ от скорости судна, слева: коэффициент попутного потока, справа: коэффициент засасывания

Для большинства судов коэффициент возрастает при увеличении скорости на винтовой характеристике. У некоторых судов коэффициент им не меняется, и есть суда, у которых функция не монотонна.

Объяснить возрастающий характер им влиянием мелководья нельзя. Для судна проекта №507А относительная глубина Я/ Т = 8,3 , но при этом его коэффициент им возрастает. Для судна проекта №11 8 - и это минимальная величина из

представленных на рис. 2 судов. Однако, для этого проекта коэффициент им примерно постоянный.

Нельзя объяснить и нефизичными натурными данными. Скольжение винта, как известно находится в пределах от 0 до 1. На рис. 3 приведены рассчитанные видимое (кажущееся) и действительное скольжение для тех же проектов.

Рис. 3. Зависимость скольжения винта от скорости судна, слева: видимое, справа: действительное

Для прояснения закономерности поведения получим его приближённое аналитическое выражение. Для этого предположим, что винтовая характеристика мощности на валу ГВ Лр выражается через степенную функцию:

Лр = л (у/ у,) (10)

где и - мощность на валу ГВ и скорость при номинальной частоте вращения ГД. В «идеальном» случае степень а = 3 , но в реальности она может принимать разные значения.

Винтовую характеристику частоты вращения вала ГВ также описываем степенной функцией:

п = п о (у/ Уо) (11)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

где п 0 - номинальная частота вращения вала ГВ. Б В «идеальном» случае = 1 , но обычно немного больше 1.

Можно показать, что кривые действия ГВ хорошо аппроксимируются функцией

7=7 о ( 1-(Яр/Я2 У2 ) , (12)

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

Используя аппроксимации (10)-(12), можно после некоторых преобразований получить из уравнения (4) явное выражение для им:

wM = 1

¿0

v/(n0D) \v0

\vJ

P Y 2

M

wn

2npD5nlK02 \v0

\vJ

a-3/3

(13)

Из (13) можно видеть, что если и , то будет монотонно

возрастающей функцией скорости .

Для простоты будем аппроксимировать натурные винтовые характеристики, проводя кривые (10) и (11) через начало координат, первую и последнюю точки натурных кривых. Результаты такой аппроксимации можно видеть для двух проектов судов из рис. 4.

Рис. 4. Аппроксимация натурных винтовых характеристик Слева: проект №1577, справа: Р77 Можно видеть, что показатели степеней и должны быть переменными возрастающими величинами, но для целей анализа достаточно такой грубой аппроксимации.

Ниже в табл. 2 приведены расчёты показатели степеней и для нескольких проектов, а также показано выполнение условия .

Таблица 2

Значения показателей степеней аппроксимации винтовых характеристик

№ Номер проекта а Р а>3р

1 11 3,363 0,978 Да

2 507А 4,52 1,06 Да

3 576 4,036 1,046 Да

4 1577 4,19 1,13 Да

5 Р77 3,03 1,03 Нет

6 558 3,855 1,078 Да

7 1553 3,958 1,248 Да

8 1570 3,974 1,227 Да

9 11 4,05 1,16 Да

10 11 2,74 1,00 Нет

11 Р32 3,56 1,19 Нет

12 507Б 3,17 0,99 Да

13 576 3,63 1,15 Да

14 576 2,83 0,93 Да

15 573 4,64 1,28 Да

16 573 3,43 1,23 Нет

17 621 3,00 1,00 Да

Примечание. Первые 8 проектов соответствуют изображённым на рис. 2. Суда с одинаковыми номерами проектов имеют разные ГВ и ГД.

Из табл. 2 можно видеть, что вариация показателей степеней в аппроксимациях (10) и (11) очень велика, в том числе у судов одного номера проекта, имеющих одинаковые размерения корпусов. При этом большое значение показателя а не гарантирует выполнение условия возрастания а > 3/?, как верно и обратное: значение не означает, что функция будет убывающей. И поскольку, это условие приближённое, то не для всех проектов оно гарантирует возрастающее поведение. Например, для проекта №1553.

Рис. 5. Аппроксимация функций = / (г>) и t = / (г>)

На рис. 5. можно видеть вычисленные по (13) аппроксимации коэффициентов попутного потока и засасывания на винтовой характеристике для ряда проектов судов. Хотя аппроксимация груба и, например, для проекта №Р77 она качественно меняет поведение функций, тем не менее, в целом рис. 5 выглядит значительно яснее. Из всего 44 судов, для которых были выполнены подобные расчёты, у 33 из них условие возрастания функции выполняется.

Можно также, видеть, что поведение коэффициентов для судов проектов №1535 и №1570 отличаются относительной стабильностью. Это суда с 3-х лопастными ГВ. Чтобы получить отдельную картину, были проведены аналогичные расчёты для судов с такими винтами, которые можно видеть на рис. 6 слева.

и р Н к. а 1Г Н <| X* а н и н 1Г н а х я ¡г

Рис. 6. Аппроксимация функций V/м _ / ( V) Слева: для судов с трёхлопастными винтами, справа: для судов проекта №573

Видимо, можно предполагать, что, действительно, трёхлопастные СКВ ведут себя в целом иначе, чем четырёхлопастные.

Практически все проекты с насадками имеют либо отрицательный, либо очень малый ( ) коэффициент при движении с минимальной скоростью.

Исключением являются суда трёх проектов: Р86М, № 276 (он показан на рис. 6) и два из проектов № 573.

Расчёты по всем проектам №573 показаны на рис. 6 справа. Можно видеть, что поведение коэффициента им чрезвычайно разнообразное и не объясняется ни наличием насадки, ни параметрами корпуса.

О постоянстве парали засасывания

Третий вопрос исследования: является ли параль засасывания константой.

Как известно, на основе данных модельных испытаний Э.Э Папмель в 30-х годах вывел следующую зависимость, которая считается справедливой для всех скольжений:

£51_с о п б I , (14)

где - относительное скольжение ГВ, - шаговое

отношение нулевого упора.

Поскольку при работе на швартовах то можно видеть, что константа в (14)

равна коэффициенту засасывания на швартовах . Исходя из оснований, которые здесь не излагаются, этот коэффициент Э.Э. Папмель назвал паралью засасывания (Само слово «параль» придумал Ф.А.Брикс - создатель паральной теории гребного винта).

Изложенное эмпирическое положение встречается только в отечественной литературе, причём практически во всех книгах, посвящённым расчётам ходкости судна. На основании (14) рекомендуется пересчитывать коэффициент засасывания с «расчётного» режима, определённого по эмпирическим формулам, на иной режим движения.

На рис. 7 приведены графики зависимости (14) при изменении скорости на винтовой характеристике.

Рис. 7. Зависимость {5! = /(г>) Слева: табличное представление зависимости, справа: аппроксимация

Можно видеть, что положение (14) не выполняется, причём даже в узком диапазоне изменения относительной поступи, которое характерно для винтовой характеристики.

Согласие натурных испытаний с теорией взаимодействия

Четвёртый вопрос исследования: насколько хорошо изменение «натурных» КВ соответствует положениям ТВ. Этот вопрос разделим на два.

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

Согласно теории взаимодействия идеального движителя с корпусом суда (ТВ) коэффициент засасывания связан с коэффициентом попутного потока функциональной зависимостью. Это положение было теоретически сформулировано впервые В. Ранкином в 1865 г., см. [22], а также прямо следует из эмпирических зависимостей для вычисления .

На рис. 8. приведены зависимости между «натурными» коэффициентами и .

Рис. 8. Зависимость { = /(и'до) Слева: табличное представление зависимости, справа: аппроксимация

Можно видеть, что связь между коэффициентами довольно близкая. Коэффициент корреляции по Пирсону приведён в табл. 3.

Таблица 3

Значения коэффициента корреляции по Пирсону

№ Номер проекта tэ ксп . f ( WМ э ксп J

1 11 0,647

2 507А 0,999

3 576 0,989

4 1577 0,994

5 Р77 0,905

6 558 0,992

7 1553 0,675

8 1570 0,652

Таким образом, теоретическое положение, «поддержанное» эмпирическими зависимостями, о тесной связи и в целом подтверждается.

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

Согласно теории взаимодействия коэффициенты взаимодействия идеального движителя с корпусом судна связаны друг с другом следующими выражениями:

(12)

где - осевая вызванная скорость в струе движителя на бесконечности; -коэффициент нагрузки движителя по упору. Заметим, что £ в (15), строго говоря, является только потенциальной частью коэффициента засасывания. То же самое касается и .

Для выполнения расчётов по натурным данным для каждого набора натурных данных определялся коэффициент засасывания, показанный на рис. 2 справа, затем по (15) определялся коэффициент попутного потока . Также был построен график зависимости , аналогичный рис. 6 слева. Результаты вычислений

показаны на рис. 9.

Рис. 9. Зависимости теоретических КВ, слева: зависимость w = f (v) , справа: зависимость t = f (w)

Сравнивая рис. 9 слева и рис. 2 справа, можно видеть, что в целом коэффициенты взаимодействия связаны почти линейной зависимостью, что соответствует как теории, так и эмпирическим зависимостям.

Выводы

В результате проведённого исследования можно сделать несколько выводов.

Во-первых, вариация значений КВ СКД при изменении параметров корпуса и ГВ согласно эмпирическим методам расчёта КВ совершенно не соответствует натурным испытаниям.

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

В-третьих, широко известная формула Папмеля для пересчёта коэффициента засасывания не соответствует натурным данным.

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

Список литературы

1. Жученко М. М., Иванов В. М. Расчёты гребных винтов. М.: Машгиз, 1953. 276 с.

2. Froude R.E. Description of a method of investigation of screw-propeller efficiency. Trans. Inst. Naval Archit. 24 (1883), P. 231.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

3. Taylor D. W. Resistance of ships and screw propulsion. - Whittaker & Co., 1893. - 234 p.

4. Папмель Э.Э. Практический расчёт гребного винта. М.-Л.: изд. журн. «Торговый флот», 1926-1936. Вып.1 и 2.

5. Horrn F. Bestimmung des Mitstroms durch Versuch «Model mit Schraube». Internationale Tagung der Leiter der Schleppversuchsanstalten. Berlin, 1937. S. 54-62.

6. Molland, A. F., Tuirnock S.R., Hudson D.A. Ship Resistance and Propulsion. Practical Estimation of Ship Propulsive Power. - Cambridge: Cambridge University Press, 2011. - p. 568. - ISBN 978-0-521-76052-2.

7. Standard Procedure for Resistance and Propulsion Experiments with Ship Models // Report of Resistance Comitee. Proceedings of 9th ITTC, 1960. pp. 237-258.

8. Басин А.М. Расчёт гребных винтов // Труды ЦНИИРФ. Вып. 13. Судовые движители. Речиздат, 1951.

9. Степанюк Е.И. Расчёт гидродинамического комплекса винт-насадка. Отчёт ЛИВТ. 1967.

10. Зорин М.П., Чудновский А.М. Выбор коэффициентов взаимодействия движителя с корпусом судна по данным натурных испытаний // Труды ЛИВТ. Теория и проектирование судов речного флота. Вып. 157. Л.: Транспорт, 1976. С. 25-29.

11. Лесюков В.А. Теория и устройство судов внутреннего плавания. Учебник для вузов водн. трансп. 4-е изд., перераб. и доп. М.: Транспорт, 1982. 303 с.

12. Harvald Sv. Aa. Resistance and Propulsion of Ships. John Willey & Sons, 1983. P. 362.

13. Lothar Birk. Fundamentals of Ship Hydrodynamics. Fluid Mechanics, Ship Resistance and Propulsion, John Wiley & Sons Ltd, 2019. P. 693.

14. Holtrop J. A Statistical Reanalysis of Resistance and Propulsion Data // International Shipbuilding Progress, 1984. Vol. 31. pp. 272-276.

15. Руководство по теплотехническому контролю серийных теплоходов. - Москва: Транспорт, 1980. - 424 с.

16. Руководство по теплотехническому контролю серийных теплоходов. - Москва: Транспорт, 1986. - 205 с.

17. Van Lammeren, W.P.A., van Manen, J.D. & Oosterveld, M.W.C. The Wageningen B-screw Series. Trans. S.N.A.M.E. 1969, 77, 269-317.

18. Васильева О.Ю. Обоснование эксплуатационно-технических характеристик грузовых судов внутреннего плавания на основе критерия энергоэффективности: автореф. дис.

канд. техн. наук: 05.22.19 / Васильева Оксана Юрьевна. - Нижний Новгород, 2022. -27 с.

19. Graff W. Untersuchungen über Änderungen von Sog und Nachstrom auf beschrankter Wassertiefe in stehenden und strömendem Wasser // Schiffstechnik, Bd. 8, H.44, 1961. S. 235-242.

20. Luthra G. Untersuchung der Nachstromverteilung an einem 2-Schrauben-Binnengutermotorschiff // Versuchsanstalt für Binnenschiffbau E.V., Duisburg, Bericht Nr. 788. H.11, 1976. S. 394-400.

21. Kulczyk J., Tabaczek T. Coefficients of Propeller-hull Interaction in Propulsion System of Inland Waterway Vessels with Stern Tunnels // The International Journal on Marine Navigation and Safety of Sea Transportation. N. 3, Vol. 8, 2014. P. 377-384.

22. Lammeren W. P. A. Resistance, Propulsion and Steering of Ships: A Manual for Designing Hull Forms, Propellers and Rudders, Band 1, H. Stam-Haarlem-Holland, 1948. P. 366.

References

1. Zhuchenko M. M., Ivanov V. M. Raschety grebnykh vintov. M.: Mashgiz, 1953. 276 s.

2. Froude R.E. Description of a method of investigation of screw-propeller efficiency. Trans. Inst. Naval Archit. 24 (1883), P. 231.

3. Taylor D. W. Resistance of ships and screw propulsion. - Whittaker & Co., 1893. - 234 p.

4. Papmel' EH.EH. Prakticheskii raschet grebnogo vinta. M.-L.: izd. zhurn. «Torgovyi floT», 1926-1936. Vyp.1 i 2.

5. Horrn F. Bestimmung des Mitstroms durch Versuch «Model mit Schraube». Internationale Tagung der Leiter der Schleppversuchsanstalten. Berlin, 1937. S. 54-62.

6. Molland, A. F., Tuirnock S.R., Hudson D.A. Ship Resistance and Propulsion. Practical Estimation of Ship Propulsive Power. - Cambridge: Cambridge University Press, 2011. - p. 568. - ISBN 978-0-521-76052-2.

7. Standard Procedure for Resistance and Propulsion Experiments with Ship Models // Report of Resistance Comitee. Proceedings of 9th ITTC, 1960. pp. 237-258.

8. Basin A.M. Raschet grebnykh vintov // Trudy TSNIIRF. Vyp. 13. Sudovye dvizhiteli. Rechizdat, 1951.

9. Stepanyuk E.I. Raschet gidrodinamicheskogo kompleksa vint-nasadka. Otchet LIVT. 1967.

10. Zorin M.P., Chudnovskii A.M. Vybor koehffitsientov vzaimodeistviya dvizhitelya s korpusom sudna po dannym naturnykh ispytanii // Trudy LIVT. Teoriya i proektirovanie sudov rechnogo flota. Vyp. 157. L.: Transport, 1976. S. 25-29.

11. Lesyukov V.A. Teoriya i ustroistvo sudov vnutrennego plavaniya. Uchebnik dlya vuzov vodn. transp. 4-e izd., pererab. i dop. M.: Transport, 1982. 303 s.

12. Harvald Sv. Aa. Resistance and Propulsion of Ships. John Willey & Sons, 1983. P. 362.

13. Lothar Birk. Fundamentals of Ship Hydrodynamics. Fluid Mechanics, Ship Resistance and Propulsion, John Wiley & Sons Ltd, 2019. P. 693.

14. Holtrop J. A Statistical Reanalysis of Resistance and Propulsion Data // International Shipbuilding Progress, 1984. Vol. 31. pp. 272-276.

15. Rukovodstvo po teplotekhnicheskomu kontrolyu seriinykh teplokhodov. - Moskva: Transport, 1980. - 424 s.

16. Rukovodstvo po teplotekhnicheskomu kontrolyu seriinykh teplokhodov. - Moskva: Transport, 1986. - 205 s.

17. Van Lammeren, W.P.A., van Manen, J.D. & Oosterveld, M.W.C. The Wageningen B-screw Series. Trans. S.N.A.M.E. 1969, 77, 269-317.

18. Vasil'eva O.YU. Obosnovanie ehkspluatatsionno-tekhnicheskikh kharakteristik gruzovykh sudov vnutrennego plavaniya na osnove kriteriya ehnergoehffektivnosti: avtoref. dis. kand. tekhn. nauk: 05.22.19 / Vasil'eva Oksana Yur'evna. - Nizhnii Novgorod, 2022. - 27 s.

19. Graff W. Untersuchungen über Änderungen von Sog und Nachstrom auf beschrankter Wassertiefe in stehenden und stromendern Wasser // Schiffstechnik, Bd. 8, H.44, 1961. S. 235-242.

20. Luthra G. Untersuchung der Nachstromverteilung an einem 2-Schrauben-Binnengutermotorschiff // Versuchsanstalt für Binnenschiffbau E.V., Duisburg, Bericht Nr. 788. H.11, 1976. S. 394-400.

21. Kulczyk J., Tabaczek T. Coefficients of Propeller-hull Interaction in Propulsion System of Inland Waterway Vessels with Stern Tunnels // The International Journal on Marine Navigation and Safety of Sea Transportation. N. 3, Vol. 8, 2014. P. 377-384.

22. Lammeren W. P. A. Resistance, Propulsion and Steering of Ships: A Manual for Designing Hull Forms, Propellers and Rudders, Band 1, H. Stam-Haarlem-Holland, 1948. P. 366.

ИНФОРМАЦИЯ ОБ АВТОРАХ / INFORMATION ABOUT THE AUTHORS

Платов Александр Юрьевич, д.т.н., доцент, зав. кафедрой прикладной информатики и статистики, Нижегородский архитектурно-строительный университет (ФГОУ ВО «ННГАСУ»), 603950, Н.Новгород, ул. Ильинская, 65, e-mail: platoff@mail.ru

Alexander J. Platov, Dr. Sci. Tech, head of Applied Informatics and Statistic Chair, Nizhny Novgorod State University of Architecture and Civil Engineering, 603950, Nizhny Novgorod, Ilyinskaya st., 65

Статья поступила в редакцию 31.07.2023; опубликована онлайн 20.09.2023. Received 31.07.2023; published online 20.09.2023.

i Надоели баннеры? Вы всегда можете отключить рекламу.