Электронная библиотека диссертаций и авторефератов России
dslib.net
Библиотека диссертаций
Навигация
Каталог диссертаций России
Англоязычные диссертации
Диссертации бесплатно
Предстоящие защиты
Рецензии на автореферат
Отчисления авторам
Мой кабинет
Заказы: забрать, оплатить
Мой личный счет
Мой профиль
Мой авторский профиль
Подписки на рассылки



расширенный поиск

Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня Мд Салим бин Камил Мд Салим бин Камил

Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня
<
Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня
>

Диссертация - 480 руб., доставка 10 минут, круглосуточно, без выходных и праздников

Автореферат - бесплатно, доставка 10 минут, круглосуточно, без выходных и праздников

Мд Салим бин Камил Мд Салим бин Камил . Решение задач определения волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов методом конечного корня: диссертация ... кандидата технических наук: 01.02.05 / Мд Салим бин Камил Мд Салим бин Камил ;[Место защиты: Санкт-Петербургский государственный морской технический университет].- Санкт-Петербург, 2015.- 277 с.

Содержание к диссертации

Введение

ГЛАВА 1 - Волновое сопротивление однокорпусных судов 50 - 76

1.1. Предыстория

1.2. Основное уравнение – выражение Мичеля-Виглея для определения волнового сопротивления 50 - 55

1.3. Основное уравнение – Формула Ньюмена для вычисления волнового сопротивления 55 - 57

1.4. Ревизия основных уравнений 57 - 58

1.5. Техника вычислений .59 - 60

1.6. Экспериментальный анализ 60 - 61

1.7. Теоретический анализ 61 - 62

1.8. Результаты экспериментального и теоретического анализа и их апробация .62 - 76

ГЛАВА 2 - Волновое сопротивление двухкорпусных судов (катамаранов) .77 - 117

2.1. Предыстория 77

2.2. Основное уравнение Така-Лазаускаса для определения волнового сопротивления двухкорпусных судов .77 - 79

2.3. Усовершенствование основного уравнения .79 - 80

2.4. Техника вычислений 80

2.5. Экспериментальный анализ 80

2.6. Теоретический анализ 80 - 81

2.7. Результаты экспериментального и теоретического анализа и их апробация 81 - 117

ГЛАВА 3 - Волновое сопротивление трехкорпусных судов (тримаранов) 118 - 133

3.1. Предыстория 118 - 119

3.2. Основное уравнение – Формула Така-Лазаускаса вычисления волнового сопротивления тримаранов 119 - 120

3.3. Усовершенствование основного уравнения 120 - 121

3.4. Методика вычислений 121

3.5. Экспериментальный анализ 121

3.6. Теоретический анализ 121 - 123

з 3.7. Результаты экспериментального и теоретического анализа и

апробация 123 - 133

ГЛАВА 4 – Параметрическое исследование 134 - 147

4.1. Основное уравнение 134 - 135

4.2. Параметрические исследования 135

4.3. Параметрическое исследование для однокорпусных судов .135 - 137

4.4. Параметрические исследования для катамаранов .137 - 143

4.5. Параметрические исследования для тримаранов 144 - 147

Заключение .148 – 152

Литература

Введение к работе

Актуальность работы

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

Решение задачи определения потребной мощности силовой установки судна непосредственно связано с необходимостью оценки его полного сопротивления, в котором основными составляющими являются сопротивление трения и остаточное сопротивление. Расчет сопротивления трения, имеющего вязкостную природу, обычно производится с применением корреляционной кривой МКОБ ' 1957, позволяющей производить пересчет с модели на натуру. Волновое сопротивление составляет наибольшую часть остаточного сопротивления и вносит наиболее существенный вклад в общее сопротивление судна. Обычно, для определения волнового сопротивления судна проводят два типа экспериментов, связанных с двумя методами: методом форм-фактора и методом измерения волнового профиля. Применение обоих методов требует использования дорогостоящих экспериментальных установок для проведения экспериментов, традиционно включающих испытания в опытовых бассейнах. При этом точность результатов существенно зависит от метода измерений и способа обработки экспериментальных данных. Необходимо отметить, что оба вышеупомянутых метода определения волнового сопротивления, ставшие общепризнанными, имеют ряд недостатков и погрешностей, связанных со средствами измерений, человеческим

1 Международная конференция опытовых бассейнов.

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

К классическим теоретическим методам определения волнового сопротивления
однокорпусного судна можно отнести пионерскую работу Мичеля (1898), выполненную
в рамках линейной теории тонкого судна и в дальнейшем развитую Виглеем (1926 -
1942), Ньюменом (1977), Таком-Лазаускасом (1998, для многокорпусных судов) и
многими другими авторами. Однако, при использовании этой теории обнаружился ряд
трудностей, сопровождающих применение соответствующих интегральных формул
определения волнового сопротивления, в частности, связанных с необходимостью
вычисления несобственных интегралов, содержащих высокочастотные

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

Работы Виглея с применение волнового интеграла Мичеля связаны с применением теории к упрощенным «математическим» корпусам и, насколько известно автору, не применялись к реальным корпусам и конфигурациям судна. Интегралы вычислялись с применением представления поверхности корпусов в виде математических выражений и использования громоздких функциональных рядов неудобных для расчета. Ньюмен и Так-Лазаускас применили аналогичный подход для определения волновых интегралов в случае однокорпусных и многокорпусных судов. Эти основополагающие методы нельзя признать удобными для применения проектантами и инженерами-кораблестроителями к судовым корпусам реальной геометрии.

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

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

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

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

Волновое сопротивление является одной из главных составляющих полного сопротивления судна. Его можно найти экспериментально с применением так называемого «метода форм-фактора» и другими методами, например, посредством измерения волнового профиля. Основная структура сопротивления судна с учетом концепции «форм-фактора» может быть представлена графически как показано на Рис. 1.2. При этот значение форм-фактора (1+к), необходимое для экспериментальной оценки сопротивления судна может быть найдено экспериментально. Соответствующий анализ производится методом Прохаски (1966), одобренным МКОБ.

Объект исследования

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

Цель диссертационного исследования

Основная цель данной диссертации состоит в вычислении интегралов волнового сопротивления линейной теории тонкого судна, полученных в работах Мичеля, Виглея, Ньюмена и Така-Лазаускаса для однокорпусных, двухкорпусных и трехкорпусных судов посредством разработанного автором метода конечного корня, что позволило найти более точное, устойчивое и сходящееся решение задачи определения волнового сопротивления. Данный метод также оказался приспособленным ко всем типам судовых корпусов, и более простым для использования в качестве инженерного инструмента прогнозирования волнового сопротивления судна, полного сопротивления и оценки потребной мощности судна в ходе проектирования. Результаты, полученные посредством предложенного автором метода конечного корня, подлежали апробации путем сравнения с существующими и полученными автором экспериментальными данными, а также с расчетными данными, полученными другими получившими признание методами.

Решаются следующие задачи

1) Провести обзор существующей литературы по исследованиям, выполненным ранее и в настоящее время с акцентом на рассмотрение различных теорий определения волнового сопротивления судна Rw а и различных составляющих полного сопротивления Rt;

  1. Проанализировать разработанные ранее и в настоящее время теории и методы теоретического и экспериментального определения волнового сопротивления;

  2. Идентифицировать и выбрать в качестве основы для дальнейших исследований одну из наиболее подходящих теорий волнового сопротивления;

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

  4. Изучить возможности усовершенствования получаемых в избранной теории математических выражений для вычисления волнового сопротивления;

  5. Разработать расчетные процедуры для волновых интегралов теории волнового сопротивления судна;

  6. Получить имеющиеся экспериментальные данные по полному сопротивлению Rt для однокорпусных и многокорпусных судов (катамаранов и тримаранов) из банков данных исследовательских центров и институтов, специализирующихся в области судостроения;

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

  8. Выполнить анализ экспериментальных данных по волновому сопротивлению Rw судна с применением метода форм-фактора и модельного коэффициента волнового сопротивления Cw на основе опытных данных по полному сопротивлению Rt;

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

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

  11. Оценить эффективность предложенных теоретических методов решения задачи определения волнового сопротивления;

  1. Продемонстрировать возможность внедрения принятой расчетной методологии для теоретического и экспериментального анализа;

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

Стелень проработанности научной проблемы

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

Эти интегралы являются несобственными определенными интегралами;

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

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

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

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

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

Научная новизна работы

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

  1. Метод конечного корня обеспечивает определенный прорыв в части расчетной реализации интегралов волнового сопротивления для однокорпусных, двухкорпусных и трехкорпусных судов в условиях высокочастотных осцилляции подынтегральных функций;

  2. Метод конечного корня приводит к устойчивому и физически обоснованному решению задачи определения волнового сопротивления с применением волновых интегралов;

  3. Метод конечного корня обеспечивает решение задач определения волнового сопротивления с удовлетворительной точностью;

  1. Метод конечного корня вносит определенный вклад в развитие научных знаний в области проектирования судов в целом и, в частности, в понимание актуальной проблемы определения волнового сопротивления судов;

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

  3. Параметрические исследования, основанные на применении метода конечного корня, могут использоваться для прогнозирования волнового сопротивления судов с геометрически подобными корпусами.

Практическая значимость работы

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

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

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

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

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

Апробация и внедрение результатов работы

осуществлялась посредством проведения экспериментальных исследований в ряде научных центров. Один из них -Малазийский технологический университет (UTM), являющийся членом Международной конференции опытовых бассейнов (МКОБ), с 1996 года располагает опытовым бассейном, оборудование которого поставлено немецкой компанией Thyssen Rheinstall Technik. Бассейн имеет длину 120 м, ширину 4 м, глубину 2.5 м, максимальную скорость движения тележки 5 м/с, ускорение 1 м/с2 при минимальном промежутке времени измерений 10 секунд. Система измерений и обработки данных (DAAS) спроектирована и поставлена нидерландским морским исследовательским центром MARIN и способна оперативно собирать данные приборов и датчиков, хранить эти данные, выполнять стандартный анализ сопротивления судна и

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

Теоретическая и методологическая основа диссертации

Соотвествующая теория, основные уравнения и методы их решения изложены в следующих параграфах.

В 1898 году Мичель на основе линейной теории разработал математический метод определения волнового сопротивления. Он сформулировал задачу определения волнового сопротивления тонкого судна, движущегося по поверхности идеальной жидкости бесконечной глубины. Используя интегральное преобразование Фурье, Мичель нашел потенциал скорости и получил связь волнового сопротивления судна и формы его корпуса в виде интеграла, широко известного сегодня как «интеграл Мичеля».

ф = (2U/n2)\\\\M,Q{[cos(nz - e)cos(n - 0 cosmtf - х)]/[ V(m2 + n2)J}expf-y<(m2 + n2)d^dCdmdn; 0 ,0 -oo < < oo

- (2U2/ng)\\\№Q[mexp(-m2U2(z + Q)/g]/[<(m2U4/^ - 1)] sin{m(x - Q + утА(т2и4/^- l)}d&Cdm; g/U2 -oo < < oo

+(2U2/Kg)\\\M,Q[mexp(-m2U2(z + 0)^]/[А(т2и4/ - 1)] cos{m(x - Q
+ упЫ(1 - m
2U4/g2)}d&Cdm; 02, 0 <С<o, -oo <

По Мичелю волновое сопротивление может быть определено посредством

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

Rw = -2р U)\ (d(f)/dx) (drj/dx) dxdz (2)

Дифференцируя потенциал скорости ф по переменной х и подставляя результат в приведенное выше выражение с учетом записанных ниже формул для интегралов I и J, получим формулу Мичеля для определения волнового сопротивления судна, широко известную как «интеграл Мичеля»;

Rw = (4pU4/ng)\(f + ^)т2/[<(т2и4^ - l)]dm = (4pU4/ng)\(f + J2) X2/[<(X2 - l)]dX;

g/U2 (3)

Где; X = mv2/g

J - \\f(x,z)exp(-X2gz/U2)sin(Xgx/U2)dxdz; -oo < x < +oo, 0 < z < oo (4)

/ = \\f(x,z)exp(-X2gz/U2)cos(Xgx/U2)dxdz; -oo < x < +oo, 0 < z < oo (5)

Дальнейшие работы, основанные на интеграле Мичеля, выполненные Виглеем (1926,

1931, 1942 и 1948), в частности, позволили получить более глубокое представление о

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

Rw = (4p^/nU2)\(P + J2)sec3ed9; 0<9<п/2 (6)

J=bd \\(Sr]/S{)exp(-dgCsec2e/U2)sm (lg&ec6/U2)d&C; ~1 <<+!, 0 (7)

I = bd \\(dr]/dtyexp(-dg(;sec26/U2)cos({g&ec6/U2)e&C; -1 <<+!, 0 (8)

Оригинальная переменная X, использованная Мичелем, была заменена на sec0,

параметры переписаны в безразмерном виде, а пределы интегрирования по % и L,

преобразованы в [-1, 1] и [0, 1] соответственно.

Затем, в 1977 году Ньюмен в своем исследовании заключил, что при постоянной плотности воды величина волнового сопротивления пропорциональна квадрату скорости судна (U2) и взвешенному интегралу квадрата волновой амплитуды (|А(0)р) с учетом куба косинуса угла распространения волн (cos30) при преобладании поперечных волн. При получении выражения для определения волнового сопротивления он связал волновое сопротивление трехмерного судна с распространением энергии волн вниз по течению. Он показал связь полученного интеграла, включающего функцию волновой амплитуды А(0), с интегралом линейной теории тонкого судна Мичеля. Ньюмен нашел непосредственную связь распространения волновой энергии в следе движущегося судна с его волновым сопротивлением в пространственном случае. Рассматривая контрольный объем,ограниченный вертикальной контрольной поверхностью х = константа и движущейся вместе с судном, можно ассоциировать поток энергии через контрольную поверхность с работой судна по преодолению волнового сопротивления. Энергия каждой волновой составляющей движется в направлении 0 с групповой скоростью Vg. Поэтому скорость прохождения энергии через контрольную поверхность, движущуюся через жидкость со скоростью U может быть найдена как (Vgcos0 - U). Полный поток энергии в направлении х может быть найдена умножением (Vgcos0 - U) на плотность энергии '/2pgA2 и интегрированием по ширине контрольной поверхности, что приводит к выражению;

dE/dt = Vipg I A2 (VgCosd - U)dz, -oo < z < oo (9)

Волновая амплитуда А может быть аппроксимирована с помощью метода стационарной фазы следующим образом;

А=(\А(в) | (2k/\RG(6) "\)0J; Vg = й Vp = 'Л Ucos9 (10)

Подставляя записанное выше выражение в интеграл потока энергии дает;

D = Trpg \{\А (в) \2/RG "(в)}(1- V2cos2e)dz; -oo < z < oo (11)

К)

dz можно заменить выражением xsec29sin6 + zsec36(l + sin26) = 0, полученное из условия d/dd{(xcos6 + zsin0)/cos20} = 0. Производя дифференцирование, получим;

х (d/d6)(sec26sine) + z (d/de){sec39(l + sin2в)} = -(dz/de){sec39(l + sin2в)}

Дважды дифференцируя фазовую функцию RG(6) = (g/U2)(xsec6 + zsec26sin6), получим;

RG "(в) = (g/U2)[x (d/de)(sec29sine) + z (d/dd){sec36(l + sin2e)}] (12)

Следовательно; \dz/d6\ = \RG"\/[(g/U2)sec30(l + sin26)]= (U2/g)cos36\RG"\/(2 - cos2в) Подставляя dz в интеграл потока энергии, получаем выражение для определения

волнового сопротивления по Ньюмену в виде;

D=RW= (жрU2/2) \\А(в)12cos36de; -ж/2 <в<ж/2 (13)

А(в) = (2vM)sec36\\пх exp[vsec2d(z+ixcose)]dzdx; d< z<0, -L/2 (14)

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

В земной системе координат (x,y,z) и связанной системе координат (хр, ур, Zjo);

Q(x,y) = Re ]Aj(Q) exp[(-ig/U2cos2в) (xpcose + y/osin9)]exp[(-ig/U2cos26)(xcos6 +

ysindjjde; -ж/2 <6<ж/2 (15)

Функция волновой амплитуды с учетом взаимодействия п корпусов;

А(в) = YAj(O) exp[(-ig/U2cos2e)(xcose + ysine)]dO;j = 1 to n (16)

Математическое выражение для волнового сопротивления двухкорпусного судна можно получить из общей формулы для волнового сопротивления многокорпусного судна при j = 1, 2, и симметричных относительно х = 0 корпусах и значениях функции волновой амплитуды, вычисленных в диаметральной плоскости при у = 0 и может быть записано в следующем виде;

Rw = (жри2/2) \\А(в)Ср\2HIFT(6)cos3ede; -ж/2 <в<ж/2 (17)

Где (npU2/2) \\A(9)cp\2cos3 Odd; -ж/2 < в < к/2 есть выражение для волнового сопротивления, полученное Ньюменом для однокорпусного судна. По Ньюмену [9], волновая функция А (в);

А(в) = (2v/n)sec36 \\(дц/дх)ехр[ь sec2в/z + ixcos 9)]dzdx; 0 (18)

По Таку-Лазаускасу фактор поперечного взаимодействия корпусов HIFt(O) для симметричного катамарана с идентичными корпусами может быть записан в виде;

HIFt/Є) = 4cos2[0.5(2p/L) sin6/(Fn2cos2d)] (19)

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

Q(x,y) = Re }Aj(6) exp[(-ig/U2cos26)(xjocos9 + yjosin6)]exp[(-ig/U2cos26)(xcose +

ysine)]d9 -ж/2 <в<ж/2 (20)

Rw = (npU2/2) l\A(e)\2cos36de; -ж/2 <в<л/2 (21)

А (в) = ZA (6)j exp[-i(g/U2cos26 (xjcosd + yjsinO)], j = 1 n; (22)

Для тримаранов, n = З

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

Rw = (npU2/2) \\A(e)jdh\2HIFt(0)jHIFl(0)jcos3ed6; -ж/2 <9<ж/2 (23)

HIFl(O) = cos2(gx/U2cos(8)); HIFT(0) = cos2[0.5(2p/L) sinO/(Fn2cos2в)] (24)

Информационная база исследования

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

  1. Исследования 1 и 2 для однокорпусных судов - экспериментальные данные по волновому сопротивлению корабля прибрежного действия (Littoral Combatant Ship - LCS) любезно предоставлены Королевским флотом Малайзии, Судостроительный завод Боустед, Малайзия и Кораблестроительным исследовательским центром HSVA, Гамбург, Германия. В HSVA контрактные испытания были проведены в октябре 1999 года с участием автора диссертации;

  2. Исследования 3, 4, 5 и 6 для однокорпусных судов - экспериментальные данные по корпусу Виглея получены из банка данных Малазийского технологического университета (Universiti Teknologi Malaysia - UTM). Эти эксперименты проводились в опытовом бассейне Морского технологического центра университета в мае 2001 года;

  1. Исследование 7 для однокорпусных судов- экспериментальные и расчетные данные по волновому сопротивлению судов, заимствованные из опубликованных работ: (і) Сансеон Цзю, Исследование полного и вязкостного сопротивления для судна Виглея с параболическими обводами, Институт гидравлических исследований Айовы, университет штата Айова, Айова Сити, Айова 52242, апрель 1983. (ii) Кампана Е., Лалли Ф., Булгарелли Ю. Численный метод решения задачи определения волнового сопротивления при нелинейных граничных условиях на свободной поверхности //Arch. Mech. 1989. том. 41. № 2-3. с. 439-447. (Ш) СВ. Сонг и Р.Е. Баддур, Расчет нелинейных корабельных волн, Международный журнал по численным методам в механике жидкости, ноябрь 1989;

  2. Исследования 1 и 2 для двухкорпусного судна - экспериментальные данные по корпусу Виглея, полученные в Малазийском технологическом университете (UTM), где эксперименты проводились в опытовом бассейне Морского технологического центра, UTM, Малайзия в мае 2001;

  3. Исследования 3 и 4 для двухкорпусного судна - экспериментальные и расчетные данные по волновому сопротивлению, опубликованные в статье: М. Инсел, А.Ф. Молланд и Дж. Ф. Велликом, Прогнозирование волнового сопротивление катамарана по линейной теории, Труды конференции САМО 94, 59-67М, 1994;

  4. Исследования 1, 2, 3 и 4 для трехкорпусных судов- экспериментальные данные по судну Виглея, полученные из Малазийского технологического университета (UTM), где эксперименты проводились автором в опутовом бассейне Морского технологического центра, UTM, Малайзия в октябре 2011;

  5. Исследование 5 трехкорпусного судна - экспериментальные данные по корпусу Виглея, полученные путем измерений волнового профилей, заимствованные из опубликованной статьи: Настия Дегиули, Андреа Вернер, Здравко Долинер, Экспериментальное определение волнового сопротивления тримарана посредством измерения волновых профилей, XVII Всемирный конгресс IMEKO, Дубровник, Хорватия, июнь 22-27, 2003;

Публикации

Прогресс в исследованиях, проводившихся диссертантом с начала 2010 года отражен в публикациях, подготовленных единолично и в со-авторстве, а также в докладах, представленных на трех международных конференциях, в том числе двух в Малайзии (ICET 2011 и MARSTEC2011) и одной в России (WMTC2012). Две статьи по применению метода конечного корня для однокорпусных и двухкорпусных судов опубликовано в журнале «Морские интеллектуальные технологии» - МИТ No 3 (25) Т.1 2014. Еще одна статья по применению метода конечного корня для трехкорпусных

судов (тримаранов) опубликована в журнале «Морской вестник» № 4 (52) декабрь 2014, ISSN 1812-3694. Более подробные ссылки на статьи, опубликованные в трудах конференций и в журналах приведены в списке публикаций по теме диссертации.

Структура работы

Основное уравнение – выражение Мичеля-Виглея для определения волнового сопротивления

Так и Лазауская выполнили подробные теоретические исследования, основанные на линейной теории, минимизируя волновое сопротивление за счет вариации продольного и поперечного расположения корпусов. При этом отдельные корпуса многокорпусного судна были приняты идентичными. Целью было минимизировать волновое сопротивление за счет оптимизации взаимного расположения корпусов и соответствующего взаимного гашения волновых систем, в частности чисто поперечных волн. Исследованы следующие конфигурации, для которых теоретические выкладки могут быть найдены в [37]: CAT - Катамаран с симметричными корпусами WEI - Катамаран не симметричный в поперечном направлении TRI - Симметричный в поперечном направлении тримаран (все корпуса идентичны) ARR - Тримаран в конфигурации стрелы (центральный корпус больше боковых) TET - Симметричный в поперечном направлении квадримаран SLI - Четырехкорпусное судно в конфигурации стрелы DIA - Четырехкорпусное судно в конфигурации бриллианта Основные теоретические результаты приведены ниже и характеристики волнового сопротивления в функции от скорости модели для различных конфигураций системы корпусов приведены на Рис. B.21 и B.22 Приложения. CAT, TRI и TET – Обычно имеют высокое волновое сопротивление при малых скоростях ARR, SLI и DIA – Минимальное волновое сопротивление, взаимное погашение поперечных волн, а DIA имеет минимальное волновое сопротивление среди всех конфигураций практически при всех скоростях B.2.5. Исследование результатов других авторов по волновому сопротивлению. В следующих параграфах рассматривается другие работы по получению формул для определения волнового сопротивления.

В 1904 год Прандтль рассматривал вопросы определения волнового сопротивления с учетом вязкости с применением результатов разработанной им теории пограничного слоя. В районе 1908 года Хавелок выполнил несколько теоретических работ по определению гидродинамических характеристик моделей судов различной конфигурации с учетом волнообразования. Виглей и Вайнблюм изготовили серию моделей судов специальной формы в 1920 году для апробации интеграла Мичеля путем сравнения теоретических результатов с экспериментальными данными.

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

В 1957 году Бхатачарья [8] рассматривал волновое сопротивление на глубокой воде ускоряющейся системы давлений. Он нашел выражение для возвышения свободной поверхности в следующем виде ; wgpy = \kV у (к) exp(ik)dk \ sin Wx exp[-p.T+ik{cr+V2ft(2t-r)}]dr;

В середине 20 века линейная теория волнового сопротивления изучалась ведущими учеными различных стран, в том числе Лунда (Норвегия), Вехаузена (США), Андерсона (Швеция), Инуи, Баба, Маруо и Бесшо (Япония), Така и Докторса (Австралия), Эггерса (Германия) и др. Эти ученые начали применять методы определения волнового профиля и анализа следа и эксперименты, проводимые в опытовом бассейнах. В 1961 В.Г. Сизов из Одесского института инженеров морского флота, внес существенный вклад в изучение волнового сопротивления, где он впервые решил трудную задачу определения волнового сопротивления с учетом нелинейных эффектов второго порядка.

В период с 1962 по 1964 Гесс и Смит разработали панельный метод определения полей скоростей и давлений вокруг полностью погруженного тела произвольной формы. Это стало важным шагом в истории развития численных методов гидродинамики. С 1962 по 1968 год Инуи провел тщательные наблюдения волновых картин. Качественное измерение волновых параметров с целью определения волнового сопротивления было важной инициативой. Изучение и интерпретация волновых структур генерированных движущимися в воде судами дало интуитивные идеи для более объективного изучения волновых систем и волнового сопротивления.

Графф, Крахт и Вайнблюм (1964) рассчитали коэффициент волнового сопротивления судна посредством интеграла Мичеля для различных моделей при различных значениях осадки и сравнили полученные результаты с коэффициентом остаточного сопротивления найденным экспериментально. Сравнение результатов теории и эксперимента также были проведены Тейлором и Гертлером. Результаты представлены на Рис. B.5. При этом можно видеть расхождения между результатами теоретических расчетов и данными эксперимента. Теоретические работы по нелинейным подходам к вычислению волнового сопротивления были выполнены Икехатой (1965), Эггерсом (1966) и Маруо (1966). Йим (1968) и Вехаузен (1969), который на писал фундаментальные работы по теории волнового сопротивления второго порядка с другими граничными условиями. В 1969 Баба предложил рассмотреть новую составляющую сопротивления по результатам теоретических и экспериментальным исследованиям полного сопротивления судов. Путем определения потери полного напора, он оценил дополнительную потерю количества движения вне следа. Предварительные теоретические исследования этой составляющей показали, что она управляется законом подобия Фруда. Изучение этой составляющей позволило ученым объяснить тот факт, что волновое сопротивление, определяемое традиционными методами представлялось слишком малым по сравнению с разностью полного и вязкостного сопротивления.

В 1976 году Гэдд выдвинул свою идею о том, что для удовлетворения нелинейных граничных условий на свободной поверхности можно применить распределение гидродинамических особенностей, то есть такой же метод как и при моделировании обтекания корпуса судна. Используя этот подход Гэдд учел нелинейные эффекты на свободной поверхности и применил точные граничные условия на поверхности тела. В 1977 году Доусон написал компьютерную программу для расчета волн и волнового сопротивления. Работа Доусона стала основой большинства вычислительных оценок волнового сопротивления посредством компьютерных программ, что сегодня является основным инструментом расчета волнового сопротивления судна. Увеличение мощности вычислительной техники стимулирует впечатляющий прогресс численных методов. Решение математических задач по оценке волнообразования стало значительно менее времяемким по сравнению с использованными ранее методами «ручного счета».

Основное уравнение Така-Лазаускаса для определения волнового сопротивления двухкорпусных судов

Предыстория. В соответствии с изложенным во Введении к диссертации, в качестве основного метода исследования, принятого в этой работе принята линейная теория Мичеля (1898), которая изучалась в целом ряде последующих работ. В качестве основного инструмента исследования проблем волнового сопротивления в рамках теории тонкого судна избрана формула Мичеля для вычисления волнового сопротивления, широко известная как интеграл Мичеля, усовершенствованная в дальнейшем Виглеем (1942), и называемая здесь формулой или интегралом Мичеля-Виглея, а также формула Ньюмена (1977), которую можно считать идентичной интегралу Мичеля. Мичель нашел потенциал скорости посредством интегральной теоремы Фурье, в то время как Ньюмен использовал метод источников и стоков с применением функции Грина. Интересно, что вывод Ньюмена приводит к некоторой форме интеграла Мичеля, которая, впрочем, несколько отличается в части использованных функций волновой амплитуды. Ниже приведены оригинальные выводы выражений для потенциалов скорости по Мичелю и Ньюмену.

Основное уравнение – выражение Мичеля-Виглея для определения волнового сопротивления. Вывод формулы Мичеля волнового сопротивления (Интеграл Мичеля) приведен ниже в соответствии с работой [3]. Линейная теория тонкого судна Мичеля развита в рамках следующих допущений и условий :

Диаметральная плоскость судна совпадает с плоскостью y = 0, а невозмущенное положение свободной поверхности – с плоскостью z = 0, ось Ox направлена в по движению судна, а ось Oz – вертикально вниз. Начало системы координат расположено на миделе судна как показано на Рис. 1.1.

Примем решение в виде acosn(z - h)cos( mx+ a )cos( py+ р ), где m2+n2+p2=0, а m вещественное число, причем поверхность воды простирается от x = - 00 до x = + оо; n и p могут быть вещественными или мнимыми, но если p является мнимым (= ip ), то последний множитель принимает вид e–p y. Это слагаемое удовлетворяет условию dф/dz = 0 при z = h и удовлетворяет уравнению; dftdz = (if/ d2 ф/dx2 Если n tann h = -U2m 2/g, это уравнение имеет бесконечное число действительных корней и один чисто мнимый корень, удовлетворяющий уравнению n tanh n h = U2m 2/g, [n = in ]. Только чисто мнимый корень непосредственно связан с волновым сопротивлением. Что касается параметра p, он всегда принимает мнимое значение для действительных корней n, а также для случая мнимых корней, если m n . Условие dф/dy = -Ud/dx = -Uf(x,z) требует разложения заданной функции f(x, z) в виде; ЦИ amncos n(z - h)cos (mx + а), где суммирование по всем значениям m принимает форму интегрирования. Предположим сначала, что функция периодична по x так что; f(x + , z) = f(x - , z),

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

Основное уравнение - Формула Ньюмена для вычисления волнового сопротивления. Поток энергии, генерируемый в результате движения судна в его следе может быть связан с волновым сопротивлением трехмерного судна [9]. Рассматривая контрольный объем, ограниченный вертикальной контрольной поверхностью x = Constant и движущийся вместе с судном, можно связать поток энергии через контрольную поверхность с работой, совершаемой судном для преодоления волнового сопротивления. Энергия каждой двумерной составляющей движется в направлении с групповой скоростью Vg. Таким образом, скорость передачи энергии через контрольную поверхность, которая движется относительно объема жидкости со скоростью U, дается выражением (Vgcos - U). Полный поток энергии в направлении x получится, если умножить выражение (Vgcos - U) на плотность энергии gA2 и проинтегрировать по ширине контрольной поверхности, что дает; Wfia dE/dt = J/2pg \Az(VsCos6 - U)dz, -co z со Волновая амплитуда А с помощью метода стационарной фазы может быть аппроксимирована следующим образом;

Ревизия основных уравнений. В данном диссертационном исследовании сделана попытка усовершенствовать основные формулы волнового сопротивления, полученные Мичелем-Виглеем и Ньюменом как показано в предыдущих параграфах на основе фундаментального предположения о том, что волновое сопротивление также зависит от конкретной величины максимального угла распространения комбинации расходящихся и поперечных волн первичной волновой системы. Этот угол, обозначаемый как pmax принимает различные значения в зависимости от формы корпуса и скорости (числа Фруда) и определяет верхний предел интегрирования вместо используемого в оригинальных формулах Мичеля-Виглея и Ньюмена значения /2.

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

Основное уравнение – Формула Така-Лазаускаса вычисления волнового сопротивления тримаранов

Методика вычислений. Здесь используется тот же метод интегрирования и методика вычислений как и в случае волнового сопротивления Rw однокорпусных и двухкорпусных судов в Главе 1 и в Главе 2. При осуществлении вычисления интеграла I (в продольном направлении) центральный корпус разбивается на 192 интервала равной длины вдоль длины модели судна (то есть 193 шпангоута или 193 узловых сечения) и 8 равных интегралов в вертикальном направлении (9 ватерлиний) в согласовании с разбивкой боковых корпусов. Боковые корпуса разбиваются на 138 равных промежутков (139 шпангоутов или узловых сечений) точно также как и в случае однокорпусных и двухкорпусных судов как сделано в Главе 1 и в Главе 2. Для целей интегрирования интервалы продольного и поперечного разбиения для центрального корпуса и боковых корпусов одинаковы. В угловом направлении интегрирование проводится при угловом интервале в один (1) градус в диапазоне 0 до 88 градусов. Кроме того, подходящий единый интервал выбирается для трех (3) последних значений угла (два равных угловых интервала) вплоть до pmax. Интегралы I для бокового и центрального корпусов вначале рассчитываются в отдельных подпрограммах при той же методике вычислений интегралов I как и в случае отдельного корпуса или для однокорпусного судна как в Главе 1. Далее вычисленные значения используются для проведения интегрирования в вертикальном направлении и при последующем вычислении волнового сопротивления Rw тримарана. Кроме коэффициента поперечного взаимодействия корпусов HIFT() как для двухкорпусного судна, в связи с наличием продольного разноса центрального и боковых корпусов катамарана, учитывается так называемый коэффициент продольного взаимодействия корпусов HIFL(). При определении Rw интегрирование проводилось систематически и последовательно с учетом диаграммы расчета, показанной в Таблицах 3.1 и 3.2 и следующей диаграммой расчета, показанной на Рис. 3.3 Приложения 5.

3.5. Экспериментальный анализ. Модельные эксперименты при различных условиях для Исследований 1 - 4 проводились в опытовом бассейне Морского технологического центра Малазийского теоретические данные по волновому сопротивлению. При этом соблюдались все необходимые подготовительные мероприятия и стандартные процедуры. Экспериментальный анализ был проведен в соответствии с рекомендациями и процедурами МКОБ 1957 и с применением метода Прохаски форм-фактора (см. Главу 1). Результаты экспериментального анализа представлены в Таблицах 3.3. - 3.7 и на графиках нанесенных на Рис. 3.4 - 3.7 и позволяют определить форм-фактор (1+k) при различных условиях и конфигураций корпуса. По завершении модельных испытаний данные анализировались по методике МКОБ с целью получения полного сопротивления и сопротивления трения с использованием метода формфактора.

Теоретический анализ. Теоретический анализ, соответствующий условиям испытаний, был выполнен в соответствии с методиками аналогичными использованным для однокорпусных и двухкорпусных судов в Главе 1 и Главе 2 с применением форматов и расчетных методик теоретического определения волнового сопротивления Rw как в Приложении. Как обычно в заключительной части расчетов итерациям и определяется конечный корень подынтегральной функции, соответствующий верхнему пределу в интеграле волнового сопротивления и представляющий собой точное значение максимального угла распространения волн pmax первичной волновой системы. Результаты теоретического анализа подытожены в приведенных ниже Таблицах, далее строятся графики зависимостей коэффициентов волнового сопротивления в функции числа Фруда, а также кривые подынтегральной функции fRw, как иллюстрация существования конечных корней. Характерные результаты расчета волнового сопротивления Rw и подынтегральные функции формы Rw [N] приведены на последующих страницах, а также в Таблицах в Приложении 5. Исследование 5 иллюстрирует апробацию метода конечного корня в сравнении с экспериментальными методами определения волнового сопротивления посредством измерения профиля волн в следе за судном (Метод Ландвебера преобразования Фурье, метод Хогбена матричных элементов, метод Баба эквивалентых особенностей). Соответствующие результаты представлены в Таблицах 3.12 - 3.13 и на Рис. 3.14.

Параметрические исследования. С учетом основного уравнения, приведенного в параграфе 6.1 и использования математически определенной формы корпуса Виглея, изменение полуширины b или B/2 приводит к изменению наклона поверхности корпуса и, следовательно, к пропорциональному изменению волновой амплитуды. Аналогичным образом, изменение глубины d или осадки T корпуса оказывает влияние на волновую амплитуду и, таки образом, на волновое сопротивление. Нижеследующие случаи параметрических исследований по однокорпусному, двухкорпусному и трехкорпусному судну рассмотрены в этой главе с целью исследования влияния параметров формы судна на волновое сопротивление с учетом результатов, полученных в предыдущих исследованиях в Главе 1, Главе 2 и Главе 3 соответственно. Результаты параметрического анализа дают полезные представления о том, как оптимизировать форму судна для минимизации волнового сопротивления. Кроме того, данные по волновому сопротивлению могут быть интерполированы для получения значений волнового сопротивления и его коэффициента при различных значениях параметров формы корпуса; полуширины b, глубины d, а также вариации этих величин. Данные параметрического исследования для однокорпусных, двухкорпусных и трехкорпусных судов, полученные в данной Главе приведены в приведенной ниже Таблице 4.1.

Параметрическое исследование для однокорпусных судов

Ревизия основных уравнений. В данном диссертационном исследовании сделана попытка усовершенствовать основные формулы волнового сопротивления, полученные Мичелем-Виглеем и Ньюменом как показано в предыдущих параграфах на основе фундаментального предположения о том, что волновое сопротивление также зависит от конкретной величины максимального угла распространения комбинации расходящихся и поперечных волн первичной волновой системы. Этот угол, обозначаемый как pmax принимает различные значения в зависимости от формы корпуса и скорости (числа Фруда) и определяет верхний предел интегрирования вместо используемого в оригинальных формулах Мичеля-Виглея и Ньюмена значения /2.

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

Исследование подынтегральной функции интеграла волнового сопротивления на предмет определения конечного корня является достаточно трудоемким. Однако найти конечный корень можно графическим путем посредством построения графика подынтегральной функции в функции от угла распространения волн . i.e. Rw [N] в функции от При этом конечный корень определяется как точка возврата, в которой имеется скачкообразный переход от отрицательного наклона кривой к положительному наклону кривой (минимум). Конечный (или последний) корень подынтегральной функции может быть определен в ходе итеративного процесса на основе некоторого начального значения с повторением вычислений вплоть до получения окончательного значения, за которым других корней при увеличении угла не имеется. Этот подход позволяет получить приемлемое согласование теоретических и экспериментальных данных. Отметим, что при дальнейшем увеличении от конечного корня до /2 подынтегральная функция интеграла волнового сопротивления, содержащая тригонометрические функции sec, sec2 и sec3 обнаруживает тенденцию к неограниченному росту вблизи /2. В дальнейшем во всех случаях, рассмотренных в данной работе, интеграл волнового сопротивления вычисляется с верхним пределом pmax. Таким образом, предложенное автором усовершенствование формул Мичеля-Виглея и Ньюмена для вычисления волнового сопротивления состоит в замене верхнего предела интегрирования /2 на pmax:

Техника вычислений. В применяемой вычислительной процедуре в всех случаях корпус судна разбивался на 138 равных частей по длине и на 8 равных частей по осадке. Интегрирование начинается в каждом шпангоутном сечении, а далее производится в направлении с целью расчета интегралов I и J при расположении начала координат на миделе как показано на Рис. 1.1. Окончательное интегрирование для вычисления волнового сопротивления Rw проводится по углу от нуля до максимального значения угла распространения волн с интервалом один градус до точного максимального значения pmax при более мелком разбиении в окрестности этого значения в направлении = /2. Вычисления проводятся в табличной форме как показано в Таблицах 1.1 - 1.6 и в соответствии с диаграммой расчетных процедур, приведенной на Рис. 1.2 в Приложении 3 и с применением следующих вычислительных процедур: