автореферат диссертации по транспорту, 05.22.19, диссертация на тему:Совершенствование теории судовождения на внутренних водных путях

доктора технических наук
Тихонов, Вадим Иванович
город
Нижний Новгород
год
2011
специальность ВАК РФ
05.22.19
цена
450 рублей
Диссертация по транспорту на тему «Совершенствование теории судовождения на внутренних водных путях»

Автореферат диссертации по теме "Совершенствование теории судовождения на внутренних водных путях"

правах рукописи

Тихонов Вадим Иванович

СОВЕРШЕНСТВОВАНИЕ ТЕОРИИ СУДОВОЖДЕНИЯ НА ВНУТРЕННИХ ВОДНЫХ ПУТЯХ

Специальность 05.22.19 - эксплуатация водного транспорта,

судовождение

Автореферат диссертации на соискание ученой степени доктора технических наук

г

19 МАЙ 2011

Н. Новгород-2011

4847760

Работа выполнена в Федеральном государственном образовательном учреждении высшего профессионального образования «Волжская государственная академия водного транспорта» (ВГАВТ, г. Нижний Новгород).

Официальные оппоненты: доктор технических наук, профессор

Ваганов Александр Борисович;

доктор технических наук, профессор Палагушкин Борис Владимирович;

доктор технических наук, профессор Юдин Юрий Иванович.

Ведущая организация: ФГОУ ВПО «Морская государственна академия имени адмирала Ф.Ф. Ушакова».

Защита диссертации состоится 16 июня 2011 г. в 14 часов на за седании диссертационного совета Д.223.001.01 при Волжской го сударственной академии водного транспорта в ауд. 231 по адресу: 603950, г. Нижний Новгород, ул. Нестерова, д. 5А.

С диссертацией можно ознакомиться в библиотеке ВГАВТ.

Автореферат разослан « "7 » мая 2011 г.

Ученый секретарь диссертационного совета, доктор технических наук, профессор

__А.Н. Ситнов

I. ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

Исследованию действующих на судно гидродинамических усилий посвящены труды A.M. Басина, Я.И. Войткунского, В.В. Вью-гова, О.И. Гордеева, А.Д. Гофмана, В.Г. Павленко, Р.Я. Першица, J1.M. Рыжова, В.Г. Соболева, К.К. Федяевского и многих других отечественных и зарубежных ученых. Разработанные к настоящему времени эмпирические методы определения корпусных сил и их моментов, основанные на результатах систематических модельных испытаний, позволяют с достаточной точностью решать задачи теории корабля, связанные с установившимся движением судна (оценка поворотливости судов, проектирование движительно-рулевого комплекса, обеспечивающего нормированную управляемость судна, и т. п.). В то же время подстановка вычисленных по

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

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

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

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

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

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

Корректность аналитических методов определения коэффициентов присоединённых масс, а также продольных и поперечных составляющих корпусных усилий проверена путём сопоставления

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

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

предложена уточнённая модель циркуляционно-отрывного обтекания судового корпуса потоком жидкости;

обоснована общая структура гидродинамических характеристик судна;

получено уравнение, аналогичное интегралам Эйлера, Бернулли и Громеки, которое позволяет определять давления реальной жидкости на поверхность обшивки судового корпуса;

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

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

откорректированы общие уравнения плоскопараллельного движения судна;

составлена частная система уравнений неустановившегося движения судна на повороте реки.

Перечисленные исследования для судов внутреннего и смешанного (река - море) плавания выполнены впервые.

Достоверность результатов. Закономерность распределения скоростей частиц жидкости в плоском пограничном слое, на основе которой получена формула зависимости экстраполятора трения от числа Рейнольдса, подтверждается экспериментальными исследованиями Я.И. Войткунского.

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

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

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

учитывающие индивидуальные геометрические характеристики судового корпуса аналитические методы определения -коэффициентов присоединённых масс;

-продольных и поперечных составляющих усилий циркуляци-онно-отрывного характера;

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

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

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

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

преподавательского состава ВГАВТ (1982-2007 г.г.), на семинарах по управляемости судов и судовождению Волжско-Камского правления НТО им. акад. А.Н. Крылова (г. Горький, 1985-1987 г.г.), на курсах повышения квалификации инженерно-технических работников МРФ РСФСР при ГИИВТе (г. Горький, 1986-1989 г.г.), на научно-промышленном форуме «Великие реки 2008» (г. Н.Новгород, 2008 г.), на международной научной конференции «Гидродинамика. Механика. Энергетические установки», посвященной 145-летию со дня рождения академика А.Н. Крылова (г. Чебоксары, 2008 г.).

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

ковским институтом инженеров водного транспорта (ГИИВТ) - в настоящее время Волжская государственная академия водного транспорта (ВГАВТ) - при выполнении научно-исследовательской работы «Подготовить к изданию справочник маневренных характеристик судов новых проектов» (тема № XII - 5.16 плана НИР и ОКР МРФ РСФСР 1985 г.) по заданию Главной инспекции по безопасности судоходства и охране объектов Министерства речного флота РСФСР. «Справочник маневренных характеристик судов» был издан заказчиком в ЦБНТИ МРФ в 1989 г. и разослан в пароходства с целью распространения на судах для практического использования судоводителями.

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

Публикации. По теме диссертации автором опубликованы двадцать две работы общим объёмом более 30 п. л., восемь из которых - в рецензируемых ВАК РФ изданиях.

Структура и объем работы. Диссертация состоит из перечня основных условных обозначений, введения, восьми разделов, заключения, списка литературы и трёх приложений. Общий объём работы составляет 287 страниц. Основная часть рукописи содержит 264 страницы, включая 9 таблиц, 28 рисунков и библиографию (123 наименования). В приложениях представлены геометрические и скоростные характеристики тринадцати типов судов речного флота, результаты расчётов углов дрейфа этих судов на установившихся циркуляциях, а также документы, подтверждающие использование результатов исследований автора.

2. КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

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

Уравнения произвольного движения твёрдого тела в идеальной жидкости впервые были получены Г. Кирхгофом ещё в 1869 г., а впервые составлены и применены для анализа движения судна в 1949 г. A.M. Васиным. В случае плоского движения на спокойной воде эти уравнения представляются следующим образом:

{™ + Лп)~--{т + Л22)ууа)-Л26со2 =Х; (1)

(т + Л22}-— + (т + Ли \>ха> + Л26 = Y; (2)

[?2+Л6б)~ + {Л22-Лп )vxvy + dvv

+Л26Ухсо + Л2в—?- = М2, (3)

где т - масса судна;

Лп,Л22 - присоединённые массы жидкости;

vx,vy - проекции вектора линейной скорости v центра масс (ЦМ) судна на оси х и у связанной с ним системы координат; t - время;

со - угловая скорость вращения судна относительно вертикальной оси z, проходящей через его ЦМ;

Л 26 ~ присоединённый статический момент;

Jz - момент инерции судна относительно вертикальной оси z;

Л66 - момент инерции присоединённых масс жидкости;

X,Y - проекции главного вектора приложенных к судну сил неинерционной природы на оси х и у подвижной системы координат;

Mz - проекция вектора главного момента приложенных

к судну сил неинерционной природы на ось z подвижной системы координат.

Усилия, обусловленные инерционностью присоединённой жидкости, определяются через её кинетическую энергию Т2, причём

Т2 = 0,+ Z22v2y +Я66а>2)+ Z26vyсо. (4)

Отмечено, что полученная Г. Кирхгофом формула (4) без соответствующих пояснений или специальных оговорок не может считаться составленной корректно, поскольку движение неравных масс воды по направлениям х и у противоречит свойству её неразрывности.

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

A.26w2, Лпухф и момента Л26ухсо в левых частях уравнений (1)-

(3) быть не должно.

Замечено, что разные способы составления уравнений движения суДна, применённые В.Г. Павленко, JI.M. Витавером, В.В. Вьюговым, О.И. Гордеевым, В.В. Саленеком и др., привели в конечном итоге к тем же выражениям (1)-(3). Это позволило предположить, что причины одних и тех же ошибок, допущенных как Г. Кирхгофом и A.M. Васиным, так и многочисленными их последователями, кроются в некорректности формулы (4), а также в подходе к системе судно - жидкость как к некоей фиктивной механической системе, масса которой меняется в зависимости от направления её движения.

Проанализированы существующие методы определения величин Лп, Ля, Л26, Л66 и Jz. Отмечено, что для судов речного

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

ширины В, осадки Т и коэффициента общей полноты S, то есть не учитывают индивидуальных особенностей обводов судна. Кроме того, коэффициент присоединённого статического момента &26 принимается равным нулю.

Усилия неинерционного происхождения Хт, Гг и Мг, действующие на подводную часть судового корпуса, принято именовать гидродинамическими характеристиками (ГДХ) судна и представлять в следующей форме:

Xr=0,5CXrpLTv2; (5)

Yr = 0,5CyrpLTv2- (б)

Мг =0,5CmrpL2Tv2. (7)

Здесь Схг,Суг>Стг - безразмерные коэффициенты гидродинамических усилий;

Р - плотность жидкости.

Задача оценки ГДХ является одной из труднейших, поэтому, начиная с конца 40-х годов прошлого столетия, на разработку методов определения ГДХ были направлены усилия М.Я. Алферьева, Н.И. Анисимовой, A.M. Басина, Г.И. Ваганова, A.B. Васильева, Я.И. Войткунского, В.Ф. Воронина, В.В. Вьюгова, А.Д. Гофмана, В.Г. Павленко, Р.Я. Першица, Л.М. Рыжова,' Г.В. Соболева, К.К. Федяевского и целого ряда других отечественных и зарубежных исследователей.

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

Суг = cxß + с2со + c^ß\ß\ + cAß\cö\+ с5Ш\Щ; (8)

С,„г = тф + m2cö + m3ß\ß\ + mAß\a\ + т5ЩЩ, (9)

где с,-, т{ - безразмерные коэффициенты, зависящие от геометрических характеристик судового корпуса;

Р - угол дрейфа судна по центру масс; со = со 1/у - безразмерная угловая скорость вращения судна.

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

Л22уусо, Л26со2, ЛиУха>, (л22 и ^2(ухсо переносятся в

правые части уравнений (1)-(3) и включаются в состав гидродинамических усилий неинерционной природы Хг> 7Г и Мг. При этом ни в одном из существующих эмпирических методов указанные инерционные усилия в явном виде не фигурируют.

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

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

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

Рассмотрена природа возникновения так называемой присоединённой жидкости. Поскольку судовой корпус занимает в окружающей его воде вполне определённый объём, то по мере движения судна происходит процесс непрерывного вытеснения движущимся корпусом некоторого (необходимого для его «присутствия» в жидкой среде) объёма воды и заполнения ею того объёма, который высвобождается судном.

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

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

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

1) обтекание судового корпуса потоком жидкости имеет явно выраженный циркуляционно-отрывной характер;

2) циркуляция скорости потока происходит лишь в носовой и кормовой оконечностях корпуса;

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

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

В векторной форме суммарное сопротивление жидкости движению судна представлено следующим образом:

Путём проецирования главного вектора Л на оси х и у, а вектора главного момента М - на ось г связанной с судном подвижной системы координат получены общие структурные формулы для определения ГДХ судна:

Я = -Кцир + -^отр +Яу+Кцг+Кд-

(10)

(11) (12)

Mv=Mm»+Mmp+My+Mw+Me. (13)

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

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

vo6T =-(у-упр)- (14)

Доказано, что в пределах пограничного слоя средняя скорость частиц присоединённой жидкости vnp = 0,5v, а средняя скорость

обтекания судна жидкостью vo6t = -0,5v. Это обстоятельство позволило сделать вывод о том, что в пределах пограничного слоя средняя удельная кинетическая энергия Ек, приходящаяся на единицу массы как присоединённой жидкости, так и жидкости, обтекающей судовой корпус, определяется выражением:

Ек = JvnpC/v = \vo6jdv = 0,5 J\dv = 0,25v2 . (15)

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

Ё = gz + — + 0,25v2 = const. (16)

Р

Здесь е ~ средняя удельная энергия, приходящаяся на единицу массы жидкости в пределах пограничного слоя;

S - ускорение свободного падения; z - аппликата рассматриваемой точки пограничного слоя относительно начала координат;

Р — давление жидкости.

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

лишь мгновенными значениями скоростей относительного движения. Поэтому уравнение (16) позволяет определять корпусные ГДХ как при установившемся, так и неустановившемся движении судна. Данное обстоятельство подтверждает обоснованность выдвинутой К.К. Федяевским гипотезы квазистационарности.

Установлена объективная закономерность распределения скоростей частиц в пограничном слое жидкости, обтекающей судовой корпус, -

. БШ—у\, (17)

^ В т-ТГ £ )

где их - скорость жидкости, обтекающей судно в пределах пограничного слоя; £ - толщина пограничного слоя.

Отмечено, что зависимость (17) полностью соответствует результатам экспериментальных исследований распределения скоростей в турбулентных гидродинамических следах и спутных струях, создаваемых сходящим с корпуса судна или его модели пограничным слоем, которые получены Я.И. Войткунским.

Путем решения интегрального соотношения Т. Кармана с использованием уравнения (17) получена следующая формула зависимости экстраполятора трения Су0 от числа Рейнольдса:

0,7434 л/Яе

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

%=0,25 ргА1ъщ. (19)

Здесь Ацг, Ьц, - соответственно амплитуда и ширина фронта волны.

Показано, что известная в гидродинамике судна формула (19) может быть преобразована к виду

С/о^-2^- (18)

I 2 —

К = -Р»пЛ№Ъ\У = Рп$Ш ■

1

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

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

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

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

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

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

Отмечено, что для точек потока, лежащих в плоскости одной и

той же ватерлинии выражение (16) примет следующий вид:

где рь р2 - давления жидкости в рассматриваемых точках потока;

V,, у2 - скорости частиц жидкости в этих точках. Показано, что проекции на оси координат х и у элементарного гидродинамического усилия ¿/<2, возникающего в какой-либо

л+0,25/>У!2 =/>2+0,25/>У| ,

(21)

точке A[x,y,z) поверхности S обшивки судового корпуса, могут быть представлены в виде

dQx - -pdS cos q cos у; (22)

dQy = pdS sin q cos у. (23)

Здесь q - измеряемый в четвертной системе счета курсовой угол нормали к ватерлинии в рассматриваемой точке л(х,у,г); у - снижение нормали п к поверхности корпуса в этой точке относительно нормали к ватерлинии.

Произведения dS cos q cos у и dS sin q cos y представляют собой проекции элементарной площади dS соответственно на плоскость мидельшпангоута и диаметральную плоскость (ДП) судна. Следовательно, в случае плоскопараллельного движения проекции dX и dY действующих на судовой корпус элементарных гидродинамических сил будут определяться выражениями:

dX = -ÁpxdSM; (24)

dY = Ар ydS д, (25)

где Др — разность динамических давлений жидкости в равноотстоящих от ДП и лежащих в плоскости одной .и той же ватерлинии точках обшивки носовой и кормовой оконечностей судового корпуса;

Др - разность динамических давлений жидкости в симметричных относительно ДП точках обшивки внешнего и внутреннего бортов корпуса судна.

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

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

2) носовая оконечность - часть корпуса, расположенная в нос от цилиндрической вставки;

3) кормовая оконечность - часть корпуса, расположенная в корму от цилиндрической вставки.

Тогда при сохранении объёмов его составных частей значения средней длины и коэффициента полноты водоизмещения носовой оконечности будут

/н=0,5стн1-/цн; <5нЦ0,551-/Уцн)//н, (26)

а кормовой -

/к = 0,5сгк£-/цк; SK =(0,5£1-/уцк)//к. (27)

Здесь сгн - коэффициент полноты носовой половины диаметрального батокса; /цн - средняя длина цилиндрической вставки в носовой

половине корпуса судна; Ры - коэффициент полноты мидельшпангоута;

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

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

ÏÏh =arctg(/1,/0,5£/?M); (28)

-arctg(/K/0,5i?/?M). - (29)

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

qH = arctg(/H,4H/0,55&J; (30)

qK=arctg{lKAK/0,5Bf3M), (31)

где Ан, Ак - коэффициенты компенсации погрешностей в определении динамических давлений.

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

А/я - cos2 уи, АУк = cos2 fK и АГа = cos2 , которые учитывают

средние снижения нормалей к поверхности обшивки относительно нормалей к ватерлиниям соответственно в носовой и кормовой оконечностях, а также в области цилиндрической вставки корпуса.

Доказано, что усилия, действующие как на корпус реального судна, так и на его аналог, тождественны. То есть геометрически эквивалентный аналог является гидродинамически адекватным аналогом реального судового корпуса. Это позволило сделать вывод о существовании возможности для разработки аналитических методов определения гидродинамических усилий, действующих на судно. И для того чтобы воспользоваться этой возможностью, необходимо представить погруженную часть судового корпуса в виде её эквивалентного аналога, а значения дн, дк, АУа, А^ и АГа определить численными методами с помощью теоретического чертежа судна.

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

Рассмотрен применённый Г. Кирхгофом и ставший классическим подход к получению выражения

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

(32)

где

V, = ; у6=а>;

<Р\=псо$дсо$у, ф2 = л вт ^ соя у; срв =и(х8т^соз^-7соз?соз/),

(33)

г икг 6 г/

Я 26 = Р ] ^ 2-<Я, = />ДХ81П$С05?'-

- у соэ д соэ эш # соэ ^ исК";

г д<р6 ч2 Л66 = р \(ре-р Нхзт^соз^-^создсоэх) "¿Я,

5

Г. Кирхгоф и получил формулу (4).

Выражения (32)—(34) со всей очевидностью свидетельствуют о том, что на оси связанной с судном системы координат проецируется не скорость присоединённой жидкости, а её масса/ Поэтому формула (4) явно противоречит принципу неразрывности жидкости и без сопутствующих оговорок или пояснений не может быть признана корректной и пригодной для составления уравнений движения судна.

Отмечено, что в любой точке Л(х,у,г) смоченной поверхности судового корпуса р\иХдиуА -ихлиул ^¿й? = 0 (здесь иХл , иуА - составляющие линейной скорости присоединённой жидкости). По-

ченного Г. Кирхгофом, а называемого, почему-то, моментом М. Мунка) в природе не существует.

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

Однако теория движения идеальной жидкости не позволяет проинтегрировать уравнение (35), поскольку скорости частиц идеальной жидкости приравниваются нулю либо на бесконечном удалении от судна, либо на неподвижных твёрдых границах водоёма. Следовательно, величина п, входящая в подынтегральную часть выражения (35) и представляющая собой не что иное, как толщину слоя присоединённой жидкости, становится неопределённой.

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

этому никакого «инерционного» момента

5

+ со2х2 бш4 д + 2ууб)хз'т4 д)соБ4 у пй5.

(35)

жающее энергию плоскопараллельного движения реальной присоединённой жидкости,

Т2 = 0,25р JJj(v2 cos4 q + v2 sin4 q + v

+ co2x2 sin4 q + 2vyco xsin4 q)cosAy dV, (36)

где V - объём погруженной части корпуса судна.

Интегрирование уравнения (36) по объёмам составных частей судового корпуса позволило представить кинетическую энергию присоединённой воды следующим выражением:

Г2 =0,5 (тий2и +тпй2и +ткй2к +ткй2к + юцй*ц +

+ k66Jzú)2 + 2mLk26vyco). (37)

В формуле (37) обозначено: тп,тк,тп - массы соответственно носовой, кормовой оконечностей и цилиндрической вставки; йХн,йун,йХк,иук,йу[1 - средние значения составляющих линейных скоростей присоединённой жидкости соответственно в носовой и кормовой оконечностях, а также в области цилиндрической вставки корпуса судна.

Коэффициенты к66 и к26 находятся по выражениям:

ГцАА,)/2£т (38)

k66 = ("ЛА'н sin2 qn + AYySjK sin2 qK+A

Тн k26 =

ArJh

0,25a2 -l2H

sin2 qu -

~АгЛ

0,25cr2 -I2K sin2 qK + Ari¡pu

í2 -I2

'цн 'ЦК

AS.

(39)

Доказано, что формула (37) уже ни при каких условиях не позволит получить «инерционный» момент (я22 -ЛпУхуу в левой

части уравнения (3), поскольку в любой из составных частей судового корпуса

дТ2 _ дТ2 _ / _ _ \ ---"Л "аТ"*' = т'Кмл (40)

Эмт

У)

Если скорости присоединённой воды приравнять к скоростям судна, то выражение (37) преобразуется в формулу (4), в которой

Ли=тки; Л22 = тк21; Л26 = mLk16 ; Л66 = J2k66, (41) а коэффициенты и к22 находятся следующим образом:

K)/2S; (42)

(43)

Л ^

к\\ = {aYhSJh cos2 +АГк8к1к COSz q

hi = [AyH5JHsin2?H + ArKSjKsin2 qK + A^J^/IS

то есть выражения (38) и (43) свидетельствуют о том, что к66 = к22.

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

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

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

Сх = -тА' cos2 р + mAz,sin2 В +

лцир -ч) ' ХЪ ^

+ (т2АХ4 - тк22)ш sin р + {mk'2(¡ - m2Ax¡ - m3Ax¡ )óJ2 ; (44)

Сицир =-m2BXx sinj3cosj3+(m2BX2 + т3В"Х2)соcosp. (45>

В формулах (44) и (45) обозначено: т = 25B/L - безразмерная масса судна;

4) - Ри (cos2 qK - cos2 qH )¡4S ; Л3 = A. (sin2 gH - sin2 qK )/4S; Л4 = fil (sin 9H cos <7h + sin ?к cos ;

¿22 = (¿Л sin2 gH + SjK sin2 qK )jlS;

^26 -

¿J 0,25o-2 ~/2K

/ N

8b 0,25<т2 -QB eos2 qH

-

21

cos2qK-

/4S;

4*5 = Pv. fe sinqK COSqH - SKlK sinqK cosqK)/l6S2 ;

A"Xi=Pl{sm2qa-sm2qK)l\928l-, BXl=fiu (sinqH COS+ sineosqK)/l6S2 ;

5*2 = A, fe sin?H cos?H - sJk sin?K cos^K)/l6^ ;

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

С, = mAXncos2 fi + m3A'co2 \

х5

Ст =-тъВ'х7с0cosр-т2ВХ2 sin2 fi +

'отр

*2 'V

+ т BXicoúnfi-m ВХ}со .

(46)

(47)

Здесь

Л0 = А. cos2 qH + cos2 )/4д; = fil (АГи cos2 qH + cos2 qK )/l 92¿3 ; B'X2 = ^ (4H cos2 gH + 4K cos2 )/9653 ; BXl = fil (лУк sin2 qn - АГк sin2 qK )/Ъ2д2 ; = A* sin2 <?н + 5JkA/k sin2 <7к )/16<52

/ \

II 2* 5Лп 0,25и2 -/2Н sin2 qH

V >

[0,25^-/2

ЦК

Sin2?«

'32<52 .

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

Xv =0,5Cv pLTv2,

(48)

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

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

СУХ={СЛ+АСЛ1+ЧЩ:' <49)

Здесь ДСу- поправка на шероховатость обшивки корпуса; К^ -

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

0. = ЗЬВ + 2ЬТ, , (50)

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

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

Принято допущение о том, что сила Ху по ширине судна не меняется, и момент этой силы относительно вертикальной оси Мух «0.

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

С054/?; (51)

(52)

где Рг - у/- число Фруда;

, В ( ,2 4 — ¡1 4 — \

4гх = м\Аг«С0* 1«+АгкС05 ч*)-

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

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

.=1А „-.-а^а. (53)

(54)

(55)

АСг = т Ап эт В соэ В:

•*цир Х1 ' '

АСх = -т2 А' совт/З;

*0Хр г-

«О-

В выражениях (53) и (54) обозначено:

01 ( ■ - - ■ -

Аг. =-— (эш (¡н соэ + эш сое с{к);

4.='

, 16 Т8

вь

16Г<52

В итоге сделан вывод о том, что для математического моделирования неустановившегося криволинейного движения судна коэффициенты продольных составляющих гидродинамических усилий, действующих на его корпус, могут быть представлены в виде: ~ С \

т

С*г =

\ Ах0

+ Су +Гг2Л1Р сое2/?

СОв /? +

+ т2Ахх^п/?|созр + ШАХз /?-тк'22\а>ып Р\;

(56)

Стх"0. (57)

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

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

представлены в следующем виде:

С =Ах$\ПРС,О$/}-{А2 -А2 -А^)С5со5Р +

+ Ат,$т2(3-ААа>5т/3+А5й)2\ Ст =В15тРсо5Р-(В2+В'2-В^)Шсо5р+

+ Вг вш^ Р - В4а> бш /3 + В5со 2,

(58)

(59)

А, = 2(íH sin qH cos qH - /к sin qK eos gK ) ; A2 = (o,25a2 -/2H)singH cos¿fH + + (o,25cr2 ~l2KjsinqK cosqK ; k{x = (jH/H eos2 + ¿>K/K eos2 ¡7K )/2<5 ; ^ =

eos2 qH-SKlKAy¡t eos2 Л

s2

Лз = /„4„ sin2 ÍH AyK sin2 Чк + /iA, ;

Д, = (o,25o-2 -72н ]аГп sin2 q„ --(o,25a¡-Цк)АГк sin2 qK + (/£„ -/Ц2КЯЦ ¡

11|^0,125<r3 -Í3H sin2 gH +

Ac — —

5 3

/ > / \

+ 0,125a3-/3K 4к Sin2 + /3 +Г3 'цн ^'цк

V J ^ J

5, =-

0,-125(7,"] -Z3H

singH cos#H -

0,125a3-Ркк

sinqK cos qK

/

B, = = = 0,5Л4; BA = 2AS ;

В

B"2= — 2 2 L

с \

2

H " ЦН

4 7

0,25 al-L

АУи cos2 qK +

f

\

0,25a2 -/2K

4K cos2 Ík

7 \

0,0625(j„ -/цН 4Hsin2?H + V У

+

f \ f л

4

Ia -Ia

'цн 'цк

Ч

4Ц-

0,0625сга -1

ЦК

У

Аи sm2qK

Отмечено, что именно моментом Мцир1 =0,5BlpI?Tvxvy, обусловленным парой неинерционных сил циркуляционной природы, приложенных к носовой и кормовой оконечностям корпуса, объясняется причина широко известной тенденции удлинённых тел разворачиваться наибольшей стороной против направления относительного движения потока воды.

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

Yy=0,5CVypLTvzy, (60)

а коэффициент Су, исходя из предположения о том, что аффинная

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

Cvr(Cf0+ACftl + KFy)^. (61)

Здесь KF - коэффициент пропорциональности между вязкостным

У

сопротивлением формы и сопротивлением трения при поперечном движении судна.

Для анализа и приближённой оценки поперечных усилий вязкостной природы были введены следующие допущения:

1) силы сопротивления воды движению судна лагом, обусловленные её вязкостью, приложены к поверхности обшивки днища судового корпуса;

2) касательные напряжения распределены по поверхности обшивки днища корпуса равномерно.

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

Суу = А3у sin2 р- ААуasmр + А5усо2 ; (62)

Cmr = В3у sin2 р-B^msmр + В5у а2, (63)

-б.

0,25 а2к-Цк

г \

0,25сг2 -/цН

\

в (р -р л

Им\ 'цн ЦК

/

/

К =С,

В_

уу

/

\

0,125^ -/Зн

+

\

/

3 73

ЦК ' Им

/ V

0,125^-/^ +Д 6В

Г +Г

'цн т'цк

¿Зу=Суу~-> В3у=0,5А4у; ВАу=2А5у;

В5 у =Сг/

уу 4 Г

0,0625<74-/4н

ч

0,0625с74-/~4к

Л /

+ А

у

/~4 -/4 'цн ЦК

Коэффициенты поперечных составляющих усилий, обусловленных волнообразованием, представлены в виде

ст = зт4 Р-А^соэт3/3 +А5(усо2эт2/3;

ст№= В2№О) соз3 £ ■+ в3№ эт4 /?- А;вт V ■ В выражениях (64) и (65) обозначено:

Аъш = Рг2 ^(1иА2Ги эт4 +1КА£ эт4 + /цЛ2ц); V

(64)

(65)

Л ь 2 Т

0,25 ан2-/2н

^ "

л

0,25<т2-/2к

^2кзт4?к +

/2 -Р

'цн 'цк

Тц

0,125ст3 -/

+

0,125а3 -Í3K

Al sin4 qK +

B_ 4 T

7K ' f

0,25^-/1

Г +Г

'цн т'цк

Га

АГН cos4 ?н +

Г \

0,2 5cr2-/2K А2Гк cos4 qK

ч J _

2

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

Суд =А!0 зт/?соз/?-/12<?й>соз/?;

Ств =В10 ьтРеа&р-В^Шсжр.

(66) (67)

Здесь

Т

Л20 - Т

S.

1Вв {sju sin qK cos qH + SjK sin qK cos qK);

¿U 0,25(71 -I2

ЦН

singH cos#H -

0,25<t2-/«k

sin#K cosgK

% - A4

2 Вв

3 T

0,125<73 -/1

singH cos;?,, +

0,125 a3-/3,

smgK cos qK

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

СУГ = Суш,р +СУотр +Cyv +Cy¡V +СУО =

= {Ai + Ale )sinpcosp- [a2 + Ale -A'2- A2 jfflTcosP +

+ {А3 + АЪу )ИП ^!sin Р-{АА+А4у)\°>\5{пР + (А5+ ASy )|©| ® +

+ A3w |sin/?]sin3 p - A4¡v JftjJsin3 p + ASw fflv sin2 p; (68) С = С +c +c +c +c =

И Г >"ию тотр ту т>Г т6>

+

Г \ i 1 дТ^ 3 3

дй)ч

+ = ^ М.2,3), (70)

r=l А=1 , ая-,-

= Ц + В]9 jsin ficosfi- [в2 + Вч + В'2- В2 )а> cos р + + (я3 + ВЪу )|sin/?|sin/?-(в4 + ВАу )\Ш\ sin P + (bs+ В5у )\Ш| ñ ■

+ B2w a cos3 р + B3w |sin/?|sin3 Р - B4w (oTjsin3 (69)

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

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

d_

dt

где Ts - кинетическая энергия системы судно-жидкость;

со ¡ - квазискорости [со j =vx;co2 = vy;co3 =со);

- трёхиндексные символы;

я¡ ~ квазикоординаты = х; я2 = у; п3 =y/ = cot)\ F¡ - квазисилы (Fl-X;F2= У; F3 = Mz).

Показано, что применение формулы (70) для составления уравнений движения судна, строго говоря, некорректно, поскольку коэффициенты при обобщенных и квазискоростях являются функциями лишь курса судна у/ = cot, который, в свою очередь, явно зависит от времени t.

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

с»»

в которых коэффициенты при скоростях явно зависят от времени /.

Для системы, состоящей из одиночного водоизмещающего судна и присоединённой жидкости, дТ3 I дп1 - 0. Однако у более сложных систем - плотовые, буксируемые или изгибаемые составы — кинетическая энергия может функционально зависеть от квазикоординат. Поэтому уравнения (71), по мнению автора, мо1уг быть использованы для математического моделирования плоского движения в жидкости любых водоизмещающих плавсредств.

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

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

2. Жёсткой связи между корпусом судна и частицами присоединённой воды не существует. Следовательно, присоединённая жидкость не совершает вместе с судном движения по криволинейной траектории и не испытывает воздействия инерционных усилий центробежного характера.

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

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

(т + Л1})—- т\>усо = X;

(72)

(73)

(74)

При составлении системы частных уравнений были сделаны следующие допущения (не противоречащие известным допущениям В.Г. Павленко):

1) скорость течения с постоянна по величине на оси судового хода и направлена по касательной к ней;

2) в поперечном сечении русла скорость течения меняется пропорционально расстоянию до мгновенного центра кривизны оси судового хода;

3) управление судном ведется так, что его ЦМ движется по криволинейной траектории, совпадающей с осью судового хода.

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

[т + Ли\сс>5в—-\т + Ли + 2т—^эш р— + v а/ \ 11 V) сН

\

+ т\ 1+2-

V

- (ш + Л22 )эт Р~ -

/

т + Л22 +2 т—

vcosB-[- + Л

+ 1 + 2 -уа соэ р + Л26 = ¥;

^со , . Л

■--ят И--

(76)

с

./г + Л66 Л26 вт р

V

V

М л

-Л2(усо$р^- = Мг. (77)

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

На основе существующих исследований, изложенных в работах В.Г. Павленко, получена формула

ту2

9 = 22——, (78)

&

где Я - радиус кривизны траектории движения ЦМ судна, которая позволяет приближённо определять углы крена судов внутреннего и смешанного плавания при неустановившемся криволинейном движении. Отмечено, что результаты вычислений угла крена в по формуле (78) для случаев установившихся циркуляций удовлетворительно соответствуют данным теоретических расчетов, выполненных В.Г. Павленко, наблюдаемым значениям угла в при проведении натурных циркуляционных испытаний, а также экспертным оценкам опытных капитанов грузовых и пассажирских судов речного флота.

С целью общей проверки корректности результатов исследований и приближенной оценки действующих на судно гидродинамических усилий сделано допущение о том, что для речных судов Ау^ »/?м, а для определения величин Ан, Ак, /1/н и Ау^ предложены выражения, полученные на основе анализа теоретических чертежей девяти грузовых и четырех пассажирских судов речного флота с выпуклыми или слегка вогнутыми ватерлиниями: для грузовых судов

Ан= 1,5-5н//Ц

л =1,5-4/дЛ 1 ;

для пассажирских судов

Ан=1,57(1-ён/Щ

4=1,57(1-4/^);} 1 >

для грузовых и пассажирских судов

АГи =ст10р(\-21цн)ЬТ/В;

АУк =стк20^(1-2 1цк)ЬТ/В,

где егно, сгко - коэффициенты полноты диаметрального батокса в

носовой и кормовой оконечностях судового корпуса, определяемые по формулам: 0-но = 4 /(0,5 -/цн ); сгко = 1К /(0,5 - /цк).

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

(81)

сопротивлением трения при продольном KF и поперечном KF

х у

движении судна предложены следующие выражения:

* I

т

KFx = 3,09-^

8 «У(1-/Ц)Г

(82)

fF ,1298* (83)

ЬУ а V L

Здесь сг = 0,5(<тн + сгк) - коэффициент полноты диаметрального

батокса.

Формулы (82) и (83) получены на основе модельных и натурных испытаний девяти грузовых и четырех пассажирских судов. В качестве результатов модельных испытаний приняты подсчитанные по методу В.В. Вьюгова значения коэффициентов CXQ общего

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

Отмечено, что в работах A.M. Басина и Я.И. Войткунского влияние кривизны судовой поверхности на величину сопротивления трения оценивается в 2-8% от значения С^ , подсчитанного по

формуле Прандтля - Шлихтинга. Кроме того, согласно экспериментальным исследованиям Я.И. Войткунского, коэффициент силы вязкостного сопротивления воды движению круглого плоского диска С, « 0,085.

С целью проверки корректности формул (82) и (83) по выражениям

Cf =-С/"Крх 7со хТ00%; Cv =Cf (l + KF )

* 0,455(IgRe) ' y fyK y!

для тринадцати типов судов были найдены значения Срх и С,.^ .

Результаты расчётов позволили сделать вывод о том, что значения KF и KF , определяемые формулами (82) и (83), удовлетвори-

х у

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

дольном движении судна, а также экспериментальным данным, полученным Я.И. Войткунским.

Для общей оценки корректности результатов теоретических исследований были выполнены расчёты значений коэффициентов присоединённых масс кп, Аг22 и коэффициентов сопротивления воды продольному движению Сх, а также углов дрейфа Р

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

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

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

3) корректность метода аналитического определения поперечных составляющих гидродинамических усилий неинерционной природы подтверждается данными натурных циркуляционных испытаний судов.

3. ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

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

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

2. Уточнена роль присоединённой жидкости в формировании действующих на судно гидродинамических усилий. Показано, что под термином «присоединённая жидкость» следует понимать некую условную составную часть динамической системы судно -жидкость, которая позволяет выразить действующие на судовой корпус усилия как инерционного, так и неинерционного происхождения.

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

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

5. Установлена объективная закономерность распределения скоростей частиц воды в плоском пограничном слое, которая не противоречит фундаментальным законам механики и подтверждается экспериментальными исследованиями Я.И. Войткунского. Данная закономерность распространяется и на турбулентное и на ламинарное движение воды в пограничном слое.

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

7. Получено уравнение, аналогичное интегралам Эйлера, Бер-нулли и Громеки, которое позволяет определять давления реальной жидкости на поверхность обшивки судового корпуса.

8. Показано, что сила волнового сопротивления численно равна произведению скоростного напора волны на площадь её фронта. Отмечено, что известная в гидродинамике формула для определения силы волнового сопротивления и выражение, полученное авто-

ром на основе упомянутого в предыдущем пункте уравнения, идентичны.

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

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

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

- масс и моментов присоединённой жидкости;

- продольных и поперечных составляющих усилий циркуляци-онно-отрывного характера;

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

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

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

14. Откорректированы общие уравнения произвольного плоского движения системы судно - жидкость.

15. Составлены частные уравнения неустановившегося движения судна на повороте реки.

16. Установлено, что корректность выполненных автором теоретических исследований подтверждается данными модельных и натурных испытаний судов.

4. ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Справочник маневренных характеристик судов / сост. В.И. Тихонов; под ред. Д.Ф. Бирюкова. - М.: ЦБНТИ МРФ РСФСР, 1989.-319 с.

2. Тихонов, В.И. Волнообразование и волновое сопротивление / В.И. Тихонов // Судовождение и безопасность плавания, водные пути, гидротехнические сооружения и экологическая безопасность судоходства: сб. науч. тр. / ВГАВТ. - Н. Новгород: ФГОУ ВПО ВГАВТ, 2006.-С. 28-33.

3. Тихонов, В.И. Исследование управления судном при переходе с прямого курса на поворот заданного радиуса / В.И. Тихонов // Труды / Горьков. ин-т инж. водн. трансп. - Горький, 1987. - Вып. 223.-С. 27-47.

4. Тихонов, В.И. Кинетическая энергия плоскопараллельного движения системы судно-жидкость / В.И. Тихонов // Труды ( ВГАВТ. - Н. Новгород, 1998. - Вып. 291. - Ч. 1. - С. 29-38.

5. Тихонов, В.И. Обоснование возможности условной замены реального судового корпуса его эквивалентным аналогом / В.И. Тихонов // Сборник трудов X Международной научной школы «Гидродинамика больших скоростей» и Международной научной конференции «Гидродинамика. Механика. Энергетические установки» (к 145-летию со дня рождения академика А.Н. Крылова). -Чебоксары: ЧПИ МГОУ, 2008. - С. 287-294.

6. Тихонов, В.И. Общие уравнения плоского движения судна / В.И. Тихонов // Труды / ВГАВТ. - Н. Новгород, 1998.-Вып. 291. -Ч. 1.-С. 58-68.

7. Тихонов, В.И. Определение параметров пограничного слоя жидкости / В.И. Тихонов // Судовождение и безопасность плавания, водные пути, гидротехнические сооружения и экологическая безопасность судоходства: сб. науч. Тр. / ВГАВТ. - Н.Новгород: ФГОУ ВПО ВГАВТ, 2005. - С. 15-21.

8. Тихонов, В.И. Основы теории динамической системы судно - жидкость / В.И. Тихонов. - Н. Новгород: ФГОУ ВПО ВГАВТ, 2007. - 262с. (монография).

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

ческие сооружения и экологическая безопасность судоходства: сб. науч. тр. / ВГАВТ. - Н.Новгород: ФГОУ ВПО ВГАВТ, 2005. - С. 41-51.

10. Тихонов, В.И. Продольные составляющие циркуляционно-отрывных усилий, возникающих на корпусе судна / В.И. Тихонов // Судовождение и безопасность плавания, водные пути, гидротехнические сооружения и экологическая безопасность судоходства: сб. науч. тр. / ВГАВТ. - НИовгород: ФГОУ ВПО ВГАВТ, 2005. - С. 51-62.

И. Тихонов, В.И. Роль присоединенной жидкости в формировании усилий, действующих на корпус судна / В.И. Тихонов // Судовождение и безопасность плавания, водные пути, гидротехнические сооружения и экологическая безопасность судоходства: сб. науч. тр. / ВГАВТ. - Н.Новгород: ФГОУ ВПО ВГАВТ, 2005. - С. 11-15.

12. Тихонов В.И. Сопротивление воды движению судна и его составляющие / В.И. Тихонов, М.В. Небасов // Судовождение и безопасность плавания, водные пути, гидротехнические сооружения и экологическая безопасность судоходства: сб. науч. тр. / ВГАВТ. -Н.Новгород: ФГОУ ВПО ВГАВТ, 2005.-С.38-41.

13. Тихонов, В.И. Уравнения Эйлера-Лагранжа для математического моделирования движения системы судно-жидкость / В.И. Тихонов//Труды / ВГАВТ. - Н. Новгород, 1998.-Вып. 291.-Ч. 1. - С. 43-48.

14. Тихонов, В.И. Характер обтекания корпуса судна потоком жидкости / В.И. Тихонов // Судовождение и безопасность плавания, водные пути, гидротехнические сооружения и экологическая безопасность судоходства: сб. науч. тр. /ВГАВТ. - Н.Новгород: ФГОУ ВПО ВГАВТ, 2005. - С.33-38.

Публикации в рецензируемых ВАК РФ изданиях

1. Тихонов, В.И. Закономерности движения жидкости в плоском пограничном слое / В.И. Тихонов // Речной транспорт (XXI век). - 2007.-№ 2. - С. 77-79.

2. Тихонов, В.И. Моделирование кинетической энергии реальной присоединенной жидкости / В.И. Тихонов // Наука и техника транспорта. - 2008. - № 4. - С. 5-8.

3. Тихонов, В.И. Общая структура гидродинамических характеристик судового корпуса / В.И. Тихонов, М.В. Небасов // Речной транспорт (XXI век). - 2006. - № 6. - С. 58-59.

4. Тихонов, В.И. Способ аналитического определения присоединённых масс жидкости / В.И. Тихонов // Речной транспорт (XXI век).-2010.-№5.-С. 81-82.

5. Тихонов, В.И. Способ определения волнового сопротивления воды движению судна / В.И. Тихонов // Речной транспорт (XXI век). - 2008. -№ 3. - С. 80-81.

6. Тихонов, В.И. Уточнение понятия «присоединенная жидкость» / В.И. Тихонов // Речной транспорт (XXI век). - 2006. - № 6. -С. 64-65.

7. Тихонов, В.И. Уточнение системы уравнений плоскопараллельного движения судна / В.И. Тихонов // Речной транспорт (XXI век). - 2008. - № 3. - С. 81-83.

8. Тихонов, В.И. Циркуляционно-отрывная модель обтекания корпуса судна / В.И. Тихонов // Речной транспорт (XXI век).-2006.-№ 5. - С. 70-72.

с

Формат 60x84 '/го- Гарнитура «Тайме». Ризография. Усл. печ. л. 1,92. Уч.-изд. л. 1,9. Тираж 150 экз. Заказ 604.

Издательско-полиграфический комплекс ФГОУ ВПО «ВГАВТ»

603950, Нижний Новгород, ул. Нестерова, 5а

Текст работы Тихонов, Вадим Иванович, диссертация по теме Эксплуатация водного транспорта, судовождение

Волжская государственная академия водного транспорта

На правах рукописи

05201151652

ТИХОНОВ Вадим Иванович

СОВЕРШЕНСТВОВАНИЕ ТЕОРИИ СУДОВОЖДЕНИЯ НА ВНУТРЕННИХ ВОДНЫХ ПУТЯХ

Специальность 05.22.19 — эксплуатация водного транспорта,

судовояедение

Диссертация на соискание ученой степени доктора технических наук

Н. Новгород — 2011

СОДЕРЖАНИЕ

с.

Основные условные обозначения 5

Введение 9

1. Анализ исследований динамического взаимодействия водоизме-щающего судна с окружающей его жидкостью и задачи настоящей работы 14

1.1. Характеристики криволинейного движения судна 14

1.2. Методы определения инерционных усилий, действующих на корпус судна при криволинейном движении 19

1.3. Методы определения гидродинамических усилий неинерционной природы 32

1.4. Цель и задачи настоящего исследования 43

2. Основные закономерности динамического взаимодействия воды с корпусом судна 46

2.1. Присоединенная жидкость и ее роль в формировании действующих на судно усилий 46

2.2. Сравнительный анализ характера обтекания крыльев и судна потоком жидкости 50

2.3. Структура гидродинамических усилий, действующих на корпус судна при его произвольном плоском движении 57

2.4. Характеристики движения жидкости в плоском пограничном слое 61

2.5. Волнообразование и волновое сопротивление 73

3. Обоснование возможности условной замены корпуса судна его эквивалентным аналогом 79

3.1. Основные допущения и замечания 79

3.2. Динамические давления воды на поверхность корпуса судна 85

3.3. Средневзвешенная нормаль к ватерлинии 94

с.

3.4. Гидродинамические усилия, действующие на корпус судна 103

4. Массы и моменты присоединенной жидкости 109

4.1. Кинетическая энергия идеальной жидкости 109

4.2. Кинетическая энергия реальной жидкости 118

4.3. Определение масс и моментов присоединенной воды 123

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

5.1. Циркуляционные усилия 135

5.2. Отрывные усилия 148

5.3. Силы вязкостной природы 156

5.4. Усилия, обусловленные волнообразованием 158

5.5. Анализ влияния крена судна на продольные составляющие корпусных усилий 164

6. Поперечные составляющие гидродинамических усилий, действующих на корпус судна 169

6.1. Циркуляционные усилия 169

6.2. Отрывные усилия 180

6.3. Силы вязкостной природы 189

6.4. Усилия, обусловленные волнообразованием 194

6.5. Анализ влияния крена судна на поперечные составляющие корпусных усилий 199

7. Уравнения плоского движения системы судно-жидкость 205

7.1. Преобразование уравнений Эйлера - Лагранжа для математического моделирования движения системы судно-жидкость 205

7.2. Составление общих уравнений движения системы судно-жидкость 211

7.3. Частная система уравнений неустановившегося движения судна на повороте реки 219

8. Анализ корректности выполненных исследований 228

8.1. Крен судна при его криволинейном движении 228

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

8.3. Вязкостное сопротивление формы при поперечном движении судна 239

8.4. Сопоставление результатов исследований с данными модель. ных и натурных испытаний судов 250

Заключение 257 Список использованной литературы 261 Приложение 1. Геометрические и скоростные характеристики судов 273 Приложение 2. Углы дрейфа судов на установившихся циркуляциях 276 Приложение 3. Документы, подтверждающие использование результатов исследований 284

ОСНОВНЫЕ УСЛОВНЫЕ ОБОЗНАЧЕНИЯ

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

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

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

В - расчетная ширина судна, м; Р — угол дрейфа в центре масс (ЦМ) судна, рад; Рм - коэффициент полноты миделыппангоута, ед.;

С - коэффициент сопротивления трения эквивалентной гладкой пла-

стины (экстраполятор трения), ед.; Сцг — коэффициент волнового сопротивления, ед.;

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

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

АУк =соз 2ук

АУн =соь2 ун

Ауц =со з2уц

дк — коэффициент полноты водоизмещения кормовой оконечности

корпуса судна, ед.; дн - коэффициент полноты водоизмещения носовой оконечности корпуса судна, ед.; е — толщина пограничного слоя воды, м; ¥г — число Фруда, ед.;

/^о — число Фруда при продольном движении судна со скоростью у(), ед.;

л

g = 9,81 м/с — ускорение свободного падения;

ук, уп — снижения нормалей к поверхностям обшивки относительно нормалей к ватерлиниям соответственно в кормовой и носовой оконечностях корпуса, рад; у — снижение нормалей к поверхности обшивки относительно нормалей к ватерлиниям в области цилиндрической вставки корпуса, рад;

г] = 1,305 • 10"6 м2 /с - кинематический коэффициент вязкости воды (при

=10°С);

Jz — момент инерции судна относительно вертикальной оси г, проходящей через его ЦМ, кгм2; Кр ,КРу — коэффициенты пропорциональности между вязкостным сопротивлением формы и сопротивлением трения соответственно при продольном и поперечном движении судна, ед.; ки,к22 — коэффициенты присоединенных масс, ед.; к26 — коэффициент присоединенного статического момента, ед.; к66 - коэффициент присоединенного момента инерции, ед.; Ь - расчетная длина судна, м; 1К — длина кормовой оконечности корпуса, м;

/ - длина носовой оконечности корпуса, м;

1Ц — длина цилиндрической вставки, м;

1Ц К — длина цилиндрической вставки в кормовой половине корпуса, м; / — длина цилиндрической вставки в носовой половине корпуса, м; 1К=1К/Ь - относительная длина кормовой оконечности корпуса, ед.; 1Н=1Н/Ь - относительная длина носовой оконечности корпуса, ед.; 1Ц К = 1Ц К / Ь - относительная длина цилиндрической вставки в кормовой

половине корпуса, ед.; / =1Ц Н/ Ь - относительная длина цилиндрической вставки в носовой

половине корпуса, ед.; Я,,, Л22 — массы присоединенной жидкости, кг;

Л26 - присоединенный статический момент, кгм;

Л6б — присоединенный момент инерции, кгм2;

т - масса судна, кг;

т = 25В/Ь — безразмерная масса судна, ед.;

Мг - момент гидродинамических сил, действующих на подводную часть судового корпуса, относительно вертикальной оси г, проходящей через ЦМ судна, Нм; р - давление жидкости, Н/м2;

- курсовые углы нормалей к ватерлиниям соответственно в кормовой и носовой оконечностях корпуса, рад; Я.к>Ян ~ средние значения курсовых углов нормалей к ватерлиниям в

кормовой и носовой оконечностях корпуса, рад; Я — радиус кривизны траектории ЦМ судна, м; Яе — число Рейнольдса, ед.; р — массовая плотность воды, кг/м3;

сгк — коэффициент полноты кормовой половины диаметрального ба-токса, ед.;

<7Н - коэффициент полноты носовой половины диаметрального баток-са, ед.;

Т ~ расчетная осадка судна, м;

Т8 - кинетическая энергия системы судно-жидкость, Нм; Тх - кинетическая энергия судна, Нм; Т2 - кинетическая энергия присоединенной воды, Нм; / - время, с;

в - угол крена судна, рад;

V - линейная скорость движения ЦМ судна, м/с;

у0 - скорость прямолинейного движения судна, м/с (в табл. П.1 - км/ч); гх >уу ~ проекции вектора линейной скорости движения ЦМ судна на

оси подвижной системы координат, м/с;

V ,ууо - проекции вектора линейной скорости движения ЦМ судна на

оси неподвижной системы координат, м/с;

Л

Ф — потенциал скорости, м/с;

Ф — среднее в пределах пограничного слоя жидкости значение потен-

Л

циала скорости, м /с; Ф1->(р] — единичные потенциалы, ед.;

Хг - продольная составляющая гидродинамической силы, действующей на подводную часть судового корпуса, Н; Уг — поперечная составляющая гидродинамической силы, действующей на подводную часть судового корпуса, Н; у/ - курс судна, рад;

О - площадь смоченной поверхности корпуса судна, м2;

со - угловая скорость вращения судна относительно вертикальной оси

г, проходящей через его ЦМ, 1/с; со = со иV - безразмерная угловая скорость ед.

ВВЕДЕНИЕ

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

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

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

Исследованию действующих на судно гидродинамических усилий посвящены труды А.М. Басина, Я.И. Войткунского, В.В. Вьюгова, О.И. Гордее-ва, А.Д. Гофмана, В.Г. Павленко, Р.Я. Першица, Л.М. Рыжова, В.Г. Соболева, К.К. Федяевского и многих других отечественных и зарубежных ученых. Разработанные к настоящему времени эмпирические методы определения корпусных сил и их моментов, основанные на результатах систематических модельных испытаний, позволяют с достаточной точностью решать задачи теории корабля, связанные с установившимся движением судна (оценка по-

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

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

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

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

>

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

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

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

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

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

новременно является гидродинамически тождественным (адекватным) аналогом реального судового корпуса.

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

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

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

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

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

На защиту выносятся:

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

— особенности циркуляционно-отрывного обтекания судового корпуса потоком воды;

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

— характер распределения скоростей частиц воды в пограничном слое;

— уравнение потенциального движения жидкости в пограничном слое.

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