Mr. J. H. Michell on the Wave–Resistance of a Ship.
Philosophical Magazine, 1898, vol.45, Ser.5, pp.106–123.
Джон Генри Мичелл (1863–1940),
Австралийский математик,
профессор математики
Университета в Мельбурне.

XI. Волновое сопротивление корабля.
Дж. Мичелл.
[*] (English)


Объектом настоящего исследования является общее решение проблемы корабельного волнообразования, производимого обводами судна, движущегося с постоянной скоростью в невязкой жидкости, приводящей к последовательному построению кривой волнового сопротивления. Сделано единственное допущение о том, что наклон касательной плоскости в любой точке его поверхности к диаметральной плоскости мал. Это условие не удовлетворяется возле днища на миделе судна у современных судов, но кажется, что это не сильно влияет на вызванные волны или сопротивление, потому что волны возникают скорее возле носа и кормы, где касательная плоскость к поверхности наклонена в направлении движения судна, чем в скругленной средней части корпуса. Пренебрежение вязкостью имеет мало последствий. Завихрение жидкости возле бортов судна, без сомнения, слегка изменяет виртуальную форму корпуса, но это влияет на угол наклона виртуальной касательной плоскости, от которого зависит волновое сопротивление, очень мало. Далее, влияние вязкости на разрушение волн, производимых носом, будет влиять на взаимодействие носовой и кормовой систем волн; но, это волны, длина которых соизмерима с длиной судна, дают основной вклад в волновое сопротивление, влияние вязкости должно быть мало. Отсюда делается заключение, что можно, следуя В. Фруду, разделить сопротивление на сопротивление трения и волновое, и это будет близко к истине.

Обзор экспериментальных работ по этому вопросу, как и краткое изложение теоретических работ Рассела, Рэнкина, двух Фрудов дано в интересной работе Вайта «Manual of Naval Architecture» (1894), часть. xi («Справочник по корабельной архитектуре»). Математика последнего времени по этому вопросу обращена, главным образом, к объяснению интересных наблюдений волновой картины, за исключением статьи сэра Томсона – Кельвина (W. Thomson (Lord Kelvin), Phil. Mag. (1886–7)), в которой критическая скорость в канале, полученная экспериментально задолго до этого Расселом, была описана математически. Можно обратиться к «Гидродинамике» Лэмба (1895) и к «Популярным лекциям» сэра В. Томсона (to Sir W. Thomson's «Popular Lectures on Navigation» (1891), p. 450) для обсуждения волновой картины.

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

Fig. 1

Возьмем вертикальную диаметральную плоскость судна  у=0, и поверхность невозмущённой жидкости z=0, ось Ох направлена по направлению движения судна и ось Oz  вертикально вниз. Мы можем предположить, что судно находится в покое, а жидкость движется назад с постоянной скоростью v на расстоянии от вызванных волн. Движение предполагается установившимся, и потенциал скорости записывается в виде –vx + f. Поскольку наклон поверхности корпуса  к плоскости у=0 везде мал, f  мал тоже, и мы можем пренебречь квадратами вызванных скоростей по сравнению с первым порядком f. На поверхности воды пусть высота волы z  будет направлена ниже главного уровня ( х,у ). Тогда

  (1)

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

уравнение давления, которое, поскольку

,   ....    (q.p.)

дает

и, следовательно, вместе с (1)

   ....   (2)

Учитывая симметрию судна относительно диаметральной плоскости у = 0, мы имеем df/dy=0, когда у = 0, за исключением той части судна, где h является полушириной на (x,z),

       (скажем),   ...   (3)

и это условие принимается на плоскости у=0 вместо поверхности судна, подтверждая то же, что и в уравнении (1). Окончательно, df/dz=0 на дне воды z = h.

Теперь рассмотрим решение для  f  в той части жидкости, где у положительно, с данными граничными условиями на  z = 0, z = h, y = 0.

Типичный член в решении

a cos n ( z–h ) cos ( mx+ a ) cos ( py+ b ),

где: m2+n2+p2=0. Здесь mдолжно быть взято действительным, когда вода простирается от  x = –¥ to x = +¥; n и p могут быть либо действительными, либо мнимыми, но если p мнимо [=ip'] последний множитель должен быть взят в форме e–p'y.

Этот член удовлетворяет df/dz=0 на  z=h и также удовлетворяет уравнению (2), если:

n tan nh = –v 2m 2/g     .....     (4)

Это уравнение имеет бесконечное число действительных корней и один чисто мнимый, который даётся

n' tanh n'h = v 2m 2/g ,    [n = in'].

Мы видим, что мнимый корень один представляет волновое сопротивление. Так как p всегда мнимое при действительном корне n, то также он мнимый корень, если m>n’.

Условие (3) теперь потребует выражение данной функции f (x, z) в форме

åå amn cos n( z h ) cos ( mx + a ),

где суммирование по m будет взято в форме интеграла.

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

f ( x +l, z ) = f ( xl, z ),

и положим

,

где r положительное целое число.

По методу Фурье

где A0n принимается, как обычно.

Поскольку функции cos n( zh ) все сопряжённые, то легко показать, что из этого мы получим

и

где:  A00 должно быть разделено пополам; и коэффициенты членов, данных мнимыми корнями, здесь, как всюду ниже, получены, если положить n=in'; так что

и также для Brn'.

Отсюда следует теорема:

Теперь пусть l станет бесконечным, и положив:

p r / l = m
p / l = dm,

мы получим:

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

n h = r p + e ,
h dn = p ,
tan nh = tane,
cos n( z – h )= (–)rcos ( nz– e ),
2n + (sin 2n h ) / h = 2n,
n
tan e = – km 2................... [k=v2/g],

а также:

tanh n'h = 1,
n = k m2,

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

Я могу отметить, что работа, подобная настоящей, даёт теорию подавления осцилляций, благодаря волнообразованию. Это я надеюсь дать в следующей статье.

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

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

Интегрирование последних двух членов по n, viz.

и двойной интеграл становится:

Вышеупомянутый интеграл есть: f( x,z ), и затем исчезает с тройным интегралом в данной формуле.

Рассматривая теперь для простоты воду бесконечной глубины, получаем сразу требуемое решение для f как

   ....   (5)

что дает:

,      when y = 0.

В этом выражении внимание должно быть обращено на множитель:

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

Оставляя преобразование интегралов с одной стороны, мы переходим к расчёту волнового сопротивления (R).

Пусть dp представляет собой возрастающее давление из-за волнового возмущения. Тогда: 

,

двойной интеграл берётся по диаметральной плоскости судна. Теперь, отмеряя от невозмущённой поверхности:

p = p + g r z – ½ r q 2 + ½ r v 2,

и отсюда:

     ...   (q.p.)

та, что:

 .

Подставляя величину df / dx, мы видим, что первый и третий интегралы в выражении для f  ничего не добавляют к сопротивлению, потому что:

,

из-за нечётного множителя sin ( x - x ), и отсюда:

(6)

где: l = mv2/g, и:

Если судно одинаково в носу и в корме, то I = 0, а начало координат берётся на миделе.

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

Видим, что:

   ....  (7)

где F( x ) меньше, чем наибольшее значение функции f(x, z) для данной величины х, и отсюда, если мы подставим большое число t, вместо ∞, как верхний предел λ, частью R пренебрежём, потому что её порядок не больше, чем:

    или    v 2/t 2,

и это исчезает, когда v = ∞ если мы возьмём:

.

В части R сохранённая l g/v2 всюду мала, так что мы можем применить круговую функцию и написать:

и      

Теперь:

,

судно имеет конечную длину и h = 0 в обеих оконечностях. Тогда, используя формулу (7):

и

Значимый член имеет порядок:

v–2 log v 2 , v –10/3 , …

и все исчезают, когда v = ∞. Следовательно, сопротивление, в конечном счёте, исчезает. Конечно, этот результат доказан только для очень коротких по сравнению с [§] глубиной воды судов.

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

Рассмотрим две элементарных площадки s, s' в ( x,z ), ( x',z' ) на борту судна, и пусть θ, θ ' будут наклонами горизонтальных линий в этих элементах к оси x. Сопротивление, благодаря этим двум элементам, равно

или, скажем:

 

Теперь, приняв на момент z = x + iy, и взяв интеграл

по ограниченной круговой области, где x и y положительны, мы получим:

,

и, следовательно [§] , используем

(8)

где: k = log 2 – g = .11593 … .

Из которого:

,

поскольку:

J0"(s) + J0'(s)/s + J0(s) = 0
J0'(s) = J1'(s)

[**]и также для Yo(s).

Теперь взяв:

и положив:

(9)

из уравнения (8).

H может быть легко посчитано с помощью механического интегратора в случае, когда 4r небольшое, и этот случай является обычным судном. Сложные таблицы J0 и J1 теперь используемые, и таблицы [**] Y0 и Y1 были посчитаны мистером Б. A. Смитом, который любезно подготовил таблицы kJ0 – Y0 и kJ1 – Y1, и приводятся в этой статье
Теперь мы имеем:

где:

 ,

и выражение для взаимного воздействия:

,

где:

r = g ( z + z' ) / v2,
s = g ( x – x' ) / v
2.

Элементы на противоположных концах судна s будут в общем большие по сравнению с единицей и с ; и в этом случае мы можем положить:

   ...   q.p.

и также для ( s – m ) и тогда приблизительно:

и сопротивление равно:

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

v2 = g l / 2 p ,

и эта формула может быть написана в виде:

.

Это даёт максимум сопротивления, приблизительно, когда

x – x' = ( n+7/8 )·l    ..........    [n - целое]

и максимум, связанный с

x – x' = ( n+3/8 )·l.

Эта формула соответствует интерференции носовой и кормовой волн, которая была так подробно рассмотрена Р.Е. Фрудом. Когда два элемента находятся на одном и том же вертикальном сечении судна, может быть дана другая форма описания. Положив x – x' = 0, интеграл может быть рассмотрен как:

.

Положим:

 l2 = ½ ( 1+ m),

что дает:

,

и

l2 = ½ (m – 1 );

следовательно:

или, если  m = cosh f ,

,

где:  K0 – Бесселя функция, введённая таким образом Греем и Матье (pp. 67, 90).

Отсюда:

,

поскольку K0' = K1 и соответствующий член в сопротивлении:

Функции K не были посчитаны, но их общий курс известен. Будет достаточно в настоящий момент следовать Грею и Матье, часть vii.

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

Пусть поверхность судна задана в виде:

y = ±c ( 1 + cos ax )( 1 + cos bz ),

между:

 x=± p /a,
z = 0   and   p /b;

Fig. 2.

так что для положительных значений:

Здесь I = 0, и

где k = v2/g; и отсюда:

который лучше рассчитывать с помощью механического интегратора.

Предположим, для примера, в квадратных футах:

v = 20 (скорость судна),
2 p /a = 200 (длина судна),
p /b = 20 (осадка ниже ватерлинии),
8c = 32 (наибольшая ширина),

тогда интеграл найден и равен - 620, и сопротивление около: R=940 lbs. wt.

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

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

Мы, как прежде, берём установившееся движение и используем такую же систему осей. Уравнение давлений даёт на поверхности:

,   как и раньше;

в то время как уравнение неразрывности для длинных волн даёт:

где h глубина воды, и это даёт:

    ...   (q.p.)

Дифференцируя по х и подставляя для f из уравнения давления, мы получаем:

или      

где с скорость свободных длинных волн.

Судно прямостенное и тянется до дна воды, кинематическое уравнение для судна имеет вид:

,

которое при дифференцировании по х и использовании уравнения давления даёт

и, конечно же,     

по остальной плоскости y = 0.

Теперь, если v > c, т.е. если скорость судна больше, чем свободных волн, уравнение

решается в форме:

      ........    (10)

где граничное условие даёт

,

или

      .......     (11)

Форма решения (10) выбрана такого порядка, чтобы расходящиеся волны были за кормой.

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

Сопротивление (R) даётся в виде:

,

так что оно бесконечно, когда скорость судна равна скорости свободных волн, и в конце концов изменяется по скорости.
Если v < c, берутся разные уравнения для потенциальной формы z :

 ,

при принятом

.

Теперь решение:

,

где:  r' 2 = ( x' – x )2 + y' 2,  и здесь нет волнового сопротивления.

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

В таблице, приведенной здесь, может быть ошибка порядка 1 на последнем месте или, возможно, 2 в величинах x > 3 или 4.


Мельбурнский Университет, 9 августа 1897 г.

Таблица для
k Jo – Y o и  k J1 – Y1
Составлена мистером Б.А. Смитом. M.C.E.

    x.      k Jo(x)–Yo(x) k J1(x)–Y1(x)    x.     k Jo(x)–Yo(x)k J1(x)–Y1(x)
.00.41.92432.7384
.014.7200100.0261.42.89722.6822
.024.027450.0453.43.87062.6286
.033.621533.3951.44.84462.5773
.043.333125.0767.45.81902.5282
.053.109120.0903.46.79402.4813
.062.925816.7695.47.76942.4362
.072.770514.4002.48.74532.3929
.082.635912 6255.49.72162.3514
.092.516811.2470.50.69832.3111
.102.409910.1457.51.67532 2729
.112.31339.2459.52.65282.2357
.122.22458.4971.53.63062.1999
.132.14287.8645.54.60882.1653
.142.06707.3230.55.58732.1319
.151.99616.8545.56.56612.0995
.161.92976.4450.57.54532 0681
.171.86716.0843.58.52482.0377
.181.80795.7642.59.50462 0083
.191.75175.4780.60.48161.9798
.201.69825.2209.61.46501.9521
.211.64724 9888.62.44561.9251
.221.59834.7779.63.42641.8988
.231.55154.5855.04.40761.8732
.241.50664.4091.65.38901.8483
.251 .46334.2476.66.37071.8241
.261.42164 0983.67.35251.8005
.271.38133.9603.68.33461.7775
.281.34243.8323.69.31691.7550
.291.30463.7131.70.29951.7329
.301.26803.6020.71.28231.7114
.311.23263.4982.72.26531.6904
.321.19813.4007.73.24851.6699
.331.16453.3091.74.23191.6496
.341.13193.2233.75.21551.6299
.351.10003.1423.76.19931.6105
.361.06903.0656.77.18331.5914
.371.03872.9932.78.16751.5728
.381.00922.9245.79.15181.5544
.39.98032.8593.80.13631.5365
.40.95192.7973.81.12111.5188

Таблица (продолжение).

    x.      k Jo(x)–Yo(x) k J1(x)–Y1(x)    x.     k Jo(x)–Yo(x)k J1(x)–Y1(x)
.82.10601.50154.8.4277–.3354
.83.09101.49134.9.4588–.2817
.84.07631.46745.0.4846–.2323
.85.06171.45095.1.5052–.1786
.86.04731.43465.2.5203–.1214
.87.03301.41855.3.5300–.0690
.88.01891.40265.4.5343.–.0159
.89.00501.38705.5.5332+.0373
.90–.00881.37155,6.5269.0892
.91–.02241.35635.7.5155.1393
.92–.03591.31135.8.4991.1873
.93–04931.32615.9.4782.2326
.94–.06251.31176.0.4527.2749
.95–.07561.29726.1.4232.3139
.96–.08851.28296.2.3900.3492
.97–.10121.26886.3.3535.3805
.98–.11381.25486.4.3140.4077
.99–.12631.21086.5.2721.4306
1.0–.13861.22716.6.2281.4488
1.1–.25481.09666.7.1826.4625
1.2–.3583.97576.8.1358.4715
1.3–.4501.86166.9.0881.4758
1.4–.5304.75267.0..0408.4755
1.5–.6008..61777.1–.0065.4704
1.6–.6604.51607.2–.0531.4609
1.7–.7101.41737.3–.0986.4471
1.8–.7500.35137.4–.1424..4290
1.9–.7805.25837.5–.1842.4071
2.0–.8017.16827.6–.2237.3814
2.1–.8142.08127.7–.2604.3524
2.2–.8180–.00237.8–.2910.3203
2.3–.8138–.08217.9–.3243.2854
2.4–.8017–.15808.0–.3510.2484
2.5–.7824–.22938.1–.3739.2093
2.6–.7560–.29598.2–.3928.1684
2.7–.7234–.35768.3–.4076.1206
2.8–.6848–.41398.4–.4182.0840
2.9–.6407–.46198.5–.4214.0408
3.0–.5920–.51008.6–.4264–.0020
3.1–.5390–.54928.7–.4242–.0438
3.2–.4823–.58238.8–.4177–.0851
3.3–.4227–.60928.9–.4071–.1251
3.4–.3606–.62999.0–.3926–.1637
3.5–.2970–.64449.1–.3744–.2001
3.6–.2320–.65259.2–.3527–.2343
3.7–.1666–.65469.3–.3277–.2657
3.8–.1014–.65059.4–.2996–.2940
3.9–.0368–.64069.5–.2688–.3194
4.0+.0266–.62519.6–.2359–.3412
4.1+.0881–.60429.7–.2012–.3594
4.2.1472–.57819.8–.1642–.3738
4.3.2035–.54729.9–.1261–.3845
4.4.2566–.512010.0–.0874–.3912
4.5.3058–.472810.1–.0181–.3911
4.63510–.430010.2–.0087–.3923
4.7.3917–.381210.3+.0302

[*] Communicated by the Author.
[†] 'On the Leading Phenomena of the Wave–Making Resistance of Ships'. Trans.Inst.Naval Architects, 1881
[‡] White, p. 470
[§] See Gray and Mathews, 'Bessel Functions' p.65, and Ex.18, p.230. The formula (8) was given by Weber.
[**] ‘Messenger of Mathematics,’ 1896.

Перевод с английского Ады Шоломовны Готман,
Новосибирская государственная академия водного транспорта.

Philosophical Magazine, 1898, vol. 45, Ser. 5, pp. 106–123.

Phil. Mag, S, .5. Vol . 45 No. 272. Jan. 1898. K