Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Скворцов, М. И. Теория и практика решения задач кораблевождения с учетом влияния систематических ошибок учебник

.pdf
Скачиваний:
15
Добавлен:
19.10.2023
Размер:
6.54 Mб
Скачать

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

противления

воды

от угла

дрейфа, т. е. считать

 

R\{Vx,a)~

= R\(VX).

Возникающая

при

этом

 

погрешность

будет

пренебрежимо

мала [6], [3]. Ограничиваясь

также

первыми

членами

разложения

функций P(VX)

 

и

fli(V«)

в

степен­

ные

ряды,

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W

A V V в = -Ж W v , в + Я 3 (qw,

W) -

 

 

 

 

 

 

 

- Я 2

( 0 ,

 

1/ш) +

Я 3 ( в ) ;

 

 

 

(2.50)

 

 

WW.b

 

 

 

 

д Р

d

R t

 

 

 

 

V.oi)

 

 

 

 

 

 

 

 

 

dV

(W

 

 

 

 

 

 

С

вполне

удовлетворительной

точностью

можно

счи­

тать,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R2{qw,W)

 

 

 

= cxf2(qw)W>;

 

 

 

 

(2.52)

 

 

 

 

Я а ( 0 ,

 

VJ

 

=

0^(0^1;

 

 

 

 

 

(2.53)

где с\, с2

— коэффициенты,

 

которые

при неизменной

осад­

 

 

 

ке

корабля

практически

постоянны;

 

 

 

Cx(qw)

 

значение

коэффициента

продольной

состав­

 

 

 

 

ляющей

 

полной

аэродинамической

силы

 

 

 

 

при

курсовом

угле

относительного

вет-

 

 

 

 

pa

qw;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С,.(0) — то

же,

при

<?w = 0.

 

 

 

 

 

 

 

Выражая

величины

 

AVW

и Vm

в узлах,

a

W — в

ме­

трах

в секунду, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

AVw

 

=

dkw,

 

 

 

 

 

 

(2.56)

где

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

_ * . - w . l * , ) - .

 

 

 

 

( М 7 )

k w

— коэффициент,

 

зависящий

только

 

от

осадки ко­

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

70

Подобно функции f\(qw), функция ^ ( ^ w ) должна опре­ деляться из результатов аэродинамических продувок мо­ дели надводной части корпуса корабля. У современных кораблей ее значения, соответствующие одинаковым кур­ совым углам относительного ветра, примерно одинаковы. Поэтому в случаях, когда данные аэродинамических ис­ следований модели корабля отсутствуют, можно пользо­ ваться приведенными в § 3.8 значениями этой функции, полученными осреднением опубликованных в печати ре­ зультатов продувок моделей кораблей нескольких типов. Возникающая при этом ошибка в учете влияния ветра будет заведомо меньшей, нежели при полном пренебреже­ нии этим влиянием.

Описание движения корабля по взволнованной поверх­ ности моря, как известно, связано с целым рядом суще­ ственных трудностей. Удовлетворительной теоретической модели, которая могла бы стать основой для предвы­ числения или. хотя бы косвенного учета величины R3 (u), входящей в выражение (2.51), пока не создано. Поэтому рассматривать вопрос об учете изменений скорости хода корабля, вызываемых волнением моря, мы не будем.

§ 2.3. О П Р Е Д Е Л Е Н И Е СООТВЕТСТВИЯ СКОРОСТИ ХОДА ЧИСЛУ ОБОРОТОВ Д В И Ж И Т Е Л Е Й ,

ПОПРАВКИ ЛАГА, РАСХОДА Т О П Л И В А ' И ДАЛЬНОСТИ

 

 

ПЛАВАНИЯ

КОРАБЛЯ

 

 

Применяемые

ныне

способы обработки наблюдений,

выполняемых

при

испытаниях

ходовых

качеств кораблей,

не свободны-

от

ряда

существенных недостатков.

Поль­

зуясь этими способами,

приходится для

исключения

влия­

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

71

[76], [85], но коренного улучшения методов обработки на­ блюдении при этом достигнуто не' было.

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

Чтобы применить к обработке наблюдений способ наи­ меньших квадратов, надо в первую очередь выразить за­ висимость между измеренными и искомыми величинами в виде уравнений поправок. Условимся о правиле знаков. Обозначим символом i номер пробега. Всем пробегам, вы­ полненным (хотя бы и на разных режимах) в одном на­ правлении (например, в направлении, совпадающем с на-

хправлением первого пробега на первом режиме), будем приписывать нечетные номера; пробегам, выполненным в обратном направлении,— четные. Продольную проекцию

вектора скорости течения, имеющую то

же направление,

что и нечетные

пробеги, будем считать

положительной,

противоположно

направленную — отрицательной.

Номер режима будем обозначать символом /. Для каж­ дого режима зададимся приближенным значением VjC ско­ рости хода (например, тем, которое соответствовало на­ значенному на том же режиме числу оборотов движите­ лей во время предыдущих испытаний). Скорость хода, со­ ответствующую при штиле заданному на режиме числу оборотов движителей, будем отыскивать в виде

Vj = VJC + AV..

(2.58)

Пренебрегая изменением скорости хода от

волнения

моря (при значительном волнении, когда это

изменение

заметно, испытания не проводятся), выразим зависимость

между

искомыми величинами

и результатами

наблюдений,

выполненных

на i-м пробеге

для определения

аванса ':

 

Vjc

+ AV} + Т , (Л/, - Np)

+ AVWi +

 

 

 

+ ( - l ) ' +

4 / =

V , i

(2.59)

 

/ dV \

 

 

 

где

 

—значение

частной производной скоро­

 

 

сти хода по числу оборотов движите- -

 

 

лей при N = Ni3

уз • (об/мин)-1 ;

1 Авансом называется отношение V: N, т. е. расстояние в морских

милях, проходимое кораблем за 60 оборотов гребного винта.

72

JV;среднее

число оборотов

движителей

на

i-м

пробеге,

об/мин;

 

 

NJ3—заданное

на /-м режиме число оборо­

тов

движителей,

об/мин;

 

 

AVwl—средняя

 

за время /-го пробега величи­

на изменения скорости хода от влия­

ния

ветра,

уз;

 

 

 

 

 

vTl—средняя

 

за

время пробега

величина

продольной

проекции

скорости

тече­

ния,

уз;

 

 

 

 

 

 

Vi — наблюденная на

i-м пробеге средняя

величина продольной

проекции

абсо­

лютной

скорости

(скорости

перемеще­

ния корабля относительно земной по­

верхности),

уз.

 

 

 

 

 

Воспользовавшись выражением (2.56)

и

представив ве­

личину коэффициента kw

как

сумму

его

значения

kwc,

найденного на предыдущих испытаниях однотипных кораб­

лей ', и поправки Akw, определяемой теперь,

получим

уравнение поправок

в виде

 

 

atWj + dtAkw

+ ( - 1 ) ' + Ч ~ h = ^

( 2 - 6 0 )

где

 

at = + l ;

 

( 2 . 6 1 )

 

 

 

d _ V % ~ s ^ w d

w f ;

( 2 б 2 )

h = Vi~

Vjc ~

T; Щ-

- dfaz.

(2.63)

Величину YJ, входящую в выражение (2.63), можно вы­ числить, пользуясь паспортной диаграммой ходовых ха­ рактеристик корабля, таблицей соответствия скорости хода оборотам движителей, составленной по результатам пре­ дыдущих испытаний, или по приближенной формуле, вы­ текающей из известной формулы В. И. Афанасьева:

Т ; = 0 , 9 ~ ^ .

(2.64)

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

1 Если

величина kw

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

73

наблюдений, выполненных

для

определения

поправки

лага:

 

 

 

 

 

 

! +

W

^

+

( - 1 ) '

+ 1 ^ = ^ ,

(2-65)

где

 

 

 

 

 

 

 

 

Л. =

3 6 0 0 ^ - ;

(2.66)

 

 

 

 

'/

 

рол1—разность

отсчетов

лага,

зафиксированная на i-м

пробеге,

мили;

 

 

 

 

Г.промежуток времени, в течение которого произо­ шло это приращение отсчета лага, секунды;

Л1—средняя

скорость1

нарастания

отсчета

лага на

i-м пробеге, уз;

 

 

 

Ал.—искомая

поправка

лага на /-м режиме,

%.

Отсюда легко получить следующее уравнение

попра­

вок:

 

 

 

 

 

д ; д л 7 + ( _ 1 / + ' ^ - / ; =

*;, -

(2.67)

где

 

 

 

 

 

а\=

Л г. 100;

 

(2.68)

 

1\ =

У-Л,

 

(2.69)

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

— при двух пробегах на режиме

д ^ = 4 - ( / 1 + /2 ); д.*, = (/; + д : ( « ; + «;); (2 . 70 )

— при трех пробегах на режиме

АЛ; = (1[ + 2Г2 + Г3):(а,1 + 2а'2 + а3),

(2.71)

1 Не следует отождествлять величину Jit со скоростью по лагу,

которая равна скорости нарастания отсчета лага, умноженной на коэффициент лага.

74

где индексами 1, 2, 3 обозначены номера пробегов на ре­ жиме.

Если к тому же пренебречь влиянием ветра, то оценки аванса и поправки лага совпадут с оценками, доставляе­

мыми традиционными

способами обработки наблюдений,

но вычисления будут

проще, поскольку их придется вести

с величинами, имеющими меньшее число значащих цифр. Представив зависимости между искомыми и измерен­ ными величинами в виде уравнений поправок (2.60) и (2.67), мы получим возможность аппроксимировать про­ дольную проекцию скорости течения не только степенной

функцией

номера пробега, как в

традиционных

способах,

но и любыми другими функциями

номера

пробега

или вре­

мени, в

частности тригонометрическими

суммами (это

важно при испытаниях в районе, где наблюдаются зна­ чительные приливо-отливные течения). Например, можно считать

 

VT! = A-J +

х2

sin шТ/ + л- 3 cos юГ;,

 

 

(2.72)

где х ь х2,

хл— искомые коэффициенты;

 

 

 

 

 

Т; — промежуток

времени

между

произволь­

 

 

ным моментом, принятым за начало счета

 

 

времени,

и

средним

 

моментом

1-го про­

 

 

бега.

 

 

 

 

 

 

 

 

 

 

Тогда,

подставив выражение

(2.72)

в уравнение

(2.60),

получим это уравнение в виде

 

 

 

 

 

 

 

afiV;

+ bnxx

+

bl2x2

+

bi3x3

+ dAkw

— lt =

v.,

 

(2.73)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^ = ( - 1 ) ' + ' s i n

cor.;

 

 

 

 

 

^•3 =

( - l) ( ' + 1 cos<o7Y

 

 

 

(2.74)

Аналогичный

вид примет

уравнение (2.67). Одни и те

же коэффициенты

х и х2,

xs

могут

включаться

в

число

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

составляемых по

результатам пробегов на разных режимах

(если

только

интервалы

времени

между

режимами

не

слишком

вели­

ки) . Это позволит

ограничиваться

меньшим,

чем

в тра­

диционных

способах, числом

пробегов

на каждом режиме,

без ущерба для

точности

аппроксимации

течения

и той

75

точности, с какой будут найдены оценки искомых вели­ чин.

При испытаниях ходовых качеств кораблей обычно производится измерение крутящего момента на гребных валах. Если рассматривать конкретный гребной винт фик­ сированного шага, установленный на конкретном корабле, то при небольших углах дрейфа крутящий момент на греб­ ном валу является функцией двух аргументов: скорости хода VXi н числа Л/,- оборотов гребного вала в минуту. Бу­ дем отыскивать крутящий момент, соответствующий за­ данному на /-м режиме числу оборотов движителей Nj3 и той скорости хода V}, которую будет иметь корабль при

этом числе оборотов при штиле, в виде

 

 

 

 

 

 

 

 

 

 

 

M j = /VIJ c

-f Л/И,,

 

 

 

 

 

(2.75)

где М)(. — счислимое

(полученное

на

предыдущих

ис­

 

 

пытаниях)

значение

крутящего

момента;

 

AAlj

— искомая

 

поправка.

 

 

 

 

 

 

 

 

 

Зависимость между искомыми и измеренными на г'-м

пробеге величинами

может

быть

выражена

в виде

 

М. =

MjC

+

Ш,

+

v., (К,, -

V,)

+

v; (Nt

-

NJS),

(2.76)

где

 

Mt

— средняя

величина

 

крутящего

момента

 

 

 

 

по измерениям на i-м пробеге;

 

 

\Х/—[~ЩГ)

—частная

производная

 

крутящего

мо-

 

 

 

1

мента по скорости

хода;

 

 

 

 

 

/ дМ \

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vy. =

 

 

—частная

производная

 

крутящего

мо­

 

 

 

мента

по

числу

оборотов

 

гребного

 

 

 

 

вала

в

минуту.

 

;

 

 

 

 

 

 

Заметим,,

что Vxl

= Vt

(— 1 ) ' + I vTl.

Учитывая

также

выражения

(2.58) и

(2.59),

получим

еще одно

уравнение

поправок

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а'ЛМ.

+

^ Д k w -

/; =

v\,

 

 

 

(2.77)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а/ =

+ 1 ;

 

 

 

 

 

 

 

(2.78)

/; =

Л4. -

Mjz -

( , V

; +

vy ) (N, -

Nj3)

-

p-jclfac.

(2.79)

Если величина Mjc принимается по результатам пре­ дыдущих испытаний испытуемого или однотипного с ним

76

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

 

•/.Nl3H

(2.80)

V

=

-;

 

ко)

30,8(1 —w)

 

'

 

Мус

 

(2.81)

 

1 +

 

(2.82)

где Vj{u)

— скорость хода, на

которой .при заданном

чис­

 

ле оборотов гребнрго винта в минуту крутя­

 

щий

момент

становится

 

равным

нулю, уз;

И-—конструктивный

шаг винта, м;

 

 

 

w — коэффициент

попутного

 

потока.

 

 

 

Для устанавливаемых на быстроходных военных ко-'

раблях

широколопастных

винтов

величина

коэффициен­

та х может выбираться из табл. 2.2

по отношению

H:D

конструктивного

шага

винта

к его

диаметру.

 

 

 

 

 

 

 

 

 

Т а б л п ц а

2.2

Значение коэффициента

•/. для широколопастных

винтов

Н: D

0,8

1,0 .

1.2

1.-1

1.6

 

1,28

1,22

1,18

1,17

1.16

Коэффициент

попутного

потока,

если

его

значение,

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

муле

Тейлора:

 

 

 

 

для винтов в

диаметральной плоскости

 

 

 

w = 0,55

-

0,05;

(2.83)

для - бортевых

винтов

 

 

 

 

 

w = 0,553 -

0,20,

(2.84)

где 3 — коэффициент

полноты

водоизмещения.

 

77

Если на предыдущих

испытаниях

крутящий

момент не

измерялся, то величины

Mjc,

можно

расцчитать,

пользуясь паспортной диаграммой ходовых характеристик корабля. Для этого на диаграмме (рис. 2.3) следует найти точку А, соответствующую заданному числу оборотов дви­ жителей Nj3 и скорости хода Vic\ заметить, какая мощ-

Рис. 2.3. Паспортная диаграмма хо­

довых

характеристик

корабля:

/ — график

N

зависимости

потребляемой

мощности

от скорости

хода

при от­

носительной

поступи Х = const;

г — г р а ­

фик зависимости потребной

мощности от

скорости хода; 3 график

зависимости по­

требляемой

мощности от

скорости

хода

при числе

оборотов движителен/У • =

const

ность Ne^ соответствует этой точке; разделив ее на чис­ ло движителей, найти мощность /V,( 1 ) , приходящуюся на

один движитель; рассчитать крутящий момент (согласно Международной системе единиц (СИ) крутящий момент должен выражаться в ньютон-метрах, но пока на кораб­ лях приходится пользоваться паспортными диаграммами, в которых за единицу измерения крутящего момента при­ нят килограмм силы на метр):

М ; с = 7 1 6 , 2 - ^ - ;

(2.85)

78

затем найти приращение ординат двух произвольных то­ чек В и С, лежащих на проходящей через точку А линии /Vj3 = const; учитывая соотношение (2.85), найти соответст­ вующее приращение крутящего момента; разделив его на приращение скорости, вычислить величину p.j. Подобным

образом по приращению ординат точек А и D,

соответст­

вующих одной и той же скорости, но разным

значениям

числа оборотов, рассчитать величину

 

V J C .

 

 

 

 

 

Если

и

паспортная

диаграмма

 

отсутствует,

величи­

на

M j U может быть'рассчитана

 

по

 

формулам:

 

 

 

 

 

 

X =

30,8 (1

 

w)

 

 

>

 

 

(2.86)

 

 

 

M

JC

= k pn

D'»,

 

 

 

 

(2.87)

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

где

X — относительная

поступь гребного

винта;

 

 

 

w — коэффициент

попутного

потока;

 

 

 

 

 

D — диаметр

винта,

м;

 

 

 

 

 

 

 

 

 

 

k2—коэффициент

 

упора,

 

снимаемый с кривой

дей­

 

 

 

ствия гребного

винта [6];

 

 

 

 

 

 

 

Р—массовая

плотность,

для

морской

воды

равная

 

 

 

104 к г с - с 2 - м - 4 ;

 

 

 

 

 

 

 

 

 

 

 

и — число оборотов

гребного

 

винта в

секунду.

 

 

Если

(на двух-

или трехвинтовом

корабле)

измерения

крутящего

момента

выполнялись

на двух или трех

валах,

то для каждого из них на каждом

 

пробеге

надо

состав­

лять свое

 

уравнение

вида

(2.77). Существенным

является

то

обстоятельство,

что

в

уравнение

(2.77)

величины

Vi

и u T i не

входят. Следовательно,

 

если производятся изме­

рения крутящего момента на валу,

 

коэффициент

измене­

ния

скорости хода

от

влияния

ветра

может

быть

опреде­

лен из уравнений, мало зависимых от тех, по которым определяется аванс. Это, -в свою очередь, позволит повы­ сить точности, с какими оценки остальных искомых вели­ чин отыскиваются из систем уравнений (2.60) и (2.67). Подобно тому, как это сделано в отношении ветра, в урав­ нениях (2.60) и (2.77) может быть учтено и то влияние, какое на скорость хода и среднюю величину крутящего момента оказывает волнение моря. Но, как уже было от­ мечено, для этого надо располагать удовлетворительным математическим описанием явления.

Повышению точности определения . искомых величин способствует учет дополнительных соотношений между ними. К их числу относится, например, формула В.И.Афа-

7?

Соседние файлы в папке книги из ГПНТБ