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

книги из ГПНТБ / Давыдов, В. В. Динамические расчеты прочности судовых конструкций учебник

.pdf
Скачиваний:
58
Добавлен:
22.10.2023
Размер:
12.6 Mб
Скачать

безразмерный коэффициент присоединенной массы В О Д Ы

 

ярс \ 1

(^£п\

(8.35)

 

р = 2Мг\

I

где $

по формуле Пабста равен

 

 

 

sp /2£д

 

(8.36)

 

1 + 0,8 2сп

 

Затем подсчитывают импульс силы, приведенный к месту ее при­

ложения при х = а, т. е. в центре удара

 

S„p = Sfi(a),

(8.37)

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

Если удар волны приходится в крыло, то определение силы удара следует выполнить по методу И. Т. Егорова [4] (для малопогруженных крыльев).

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

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

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

Втеории И. Т. Егорова учитывается сжимаемость воды и воздуха

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

Сила удара определяется следующей формулой

= 1 + ц

(8.38)

где М — масса судна;

[см. формулу (8.23)];

г) — коэффициент эксцентричности удара

v — скорость встречи крыла с волной в начальный момент удара,

определяющаяся по формуле (8.19),

при ф2 = ф (т. е. при

Фх = 0);

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

 

г Шх\

(8.39)

 

^

а $

по формуле Пабста

 

 

4 M f ) =

(8.40)

 

 

1 + 0 , 8 -

В формулах (8.38) — (8.39):

10 В. В. Давыдов, Н. В. Маттес

289

b — длина хорды профиля крыла (вдоль судна);

/— размах крыла (поперек судна);

£—■логарифмический декремент затухания колебаний, с-1.

На основании формул И. Т. Егорова А. В. Уткиным [10] была получена следующая зависимость логарифмического декремента ко­ лебаний £ от коэффициента присоединенной массы воды р, — хорда крыла):

 

м.

£— , м/с

м-

S— . м/с

 

 

2

 

2

0

 

41,3

1 , 0

125 .

0

, 1

46,1

1 , 1

153

0

, 2

50,4

1 , 2

182

0,3

57,3

1,3

217

0,4

64,3

1,4

252

0,5

72,3

1,5

288

0 , 6

80,4

1 , 6

329

0,7

88,4

1,7

369

0 , 8

99,2

1 , 8

413

0,9

ПО

1,9

454

 

 

 

2 , 0

496

Зависимвсть силы удара волны в крыло судна от времени носит характер, показанный на рис. 64.

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

Действительно, период встречи судна с волнами определяют по выражению (с, м, м/с)

 

г

 

^в

 

1качки —

 

Если >.в

изменяется в пределах 20—40 м, a vc — судна при ходе

в водоизмещающем состоянии

11—17 м/с, то ткачки находится в пре­

делах от 1

до 2 с.

 

 

Период свободных изгибных колебаний корпуса первого тона у

судов на подводных крыльях составляет около тх = 0,2

0,25 с,

т. е. в четыре—десять раз меньше. Поэтому и не требуется учитывать

динамичность

волнения.

Частоту

и форму свободных колебаний f1 для первого тона оп­

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

Решение уравнений (5.58) относительно X? и f\ выполняют на ЭЦВМ. Вращение учитывают поправкой по формуле (5.43), так как поправка эта невелика для судов на подводных крыльях, введение же вращения в дифференциальные уравнения сильно усложняет их решение. В результате применения метода конечных разностей полу­ чаются частота Я,х, форма колебаний первого тона /у и ее производные по формулам (5.55).

290

Дифференциальное уравнение относительно обобщенной коорди­ наты изгибных .колебаний первого тона имеет вид (8.14)

Ф1+ 2р,^1+ Я2ф1= ^ 1 - ,

(8.41)

причем ФП находят по формуле (8.33), а

 

L

 

М 1= Jm /^x « ALlmf*.

(8.42)

Для уточнения коэффициента сопротивления рг необходимо про­

вести ряд дополнительных экспериментальных исследований

судов

на подводных крыльях. Пока

рекомен­

 

 

 

дуем

приближенно полагать

по

книге

 

 

 

А. Д. Звягина и В. В. Шабарова

[5]

 

 

 

Pi

0,07

 

 

 

 

 

 

 

 

Функцию P(t), входящую в выра­

 

 

 

жение (8.33) для обобщенной силы Ф1;

 

 

 

можно взять из выражений (8.31).

 

 

 

 

Для решения

уравнения

движения

 

 

 

(8.41)

надо применить метод элементар­

 

 

 

ных импульсов,

изложенный

в § 4,

где

 

 

t

дано решение

уравнения (1.54),

совер­

Рис. 64.

Закон изменения силы

шенно такого же,

как и уравнение (8.41),.

только постоянную г следует заменить

удара волны в крыло

через

рг.

 

 

 

 

 

 

 

 

Полагая, что начальное отклонение и начальная скорость равны

нулю, получим согласно формуле (1.56)

 

 

 

 

д,

/

s=i

 

 

 

s=/

\

ф, = -------ё~р<]

sin Ы; 2

cos kts — cos ktj

2 Ф5ер sin Xts ) .

2/.,.VI

\

s=0

 

 

 

s—0

J

 

 

 

 

 

 

 

 

 

(8.43)

Вычисления,

применительные

к

формуле (8.43),* можно

вести

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

Фмакс-

Если удар приходится на плоское днище, то по формулам (8.34) — (8.37) определяют импульс силы, приведенной к центру удара 5 пр; по импульсу находят обобщенную координату по формуле (1.60).

Движение носит колебательный характер, имеет частоты колеба­

ний первого тона

и амплитуду обобщенной координаты

 

 

4 W c = | ^ .

(8.44)

* В формуле (8.43) опущены индексы, подчеркивающие номер тона колеба­

ний.

10*

291

Когда найдено максимальное значение обобщенной координаты Фх максможно определить амплитуду колебаний кор'Пуса в любой точке

И * ) = Ы * )^ 1 „ а к с -

(8.45)

Для определения изгибающего момента и срезывающей силы по формулам (8.17) и (8.18) необходимо вычислить вторые и третьи про­

изводные fi (х), Д (х), а также первую производную I (х). Это вы­ полняется методом конечных разностей по формулам (5.55).* Напря­ жения от общей вибрации в корпусе судна при ударе волны опре­ деляют по обычным зависимостям.

П р и м е р . Определить динамический изгибающий

момент в районе тео­

ретического шпангоута № 9 судна на

подводных

крыльях «Ракета» при ходе

на волнении и ударе волн размерами

1 2 X 1 , 2 м в корпус.

Главные размерения судна

 

Длина судна расчетная L, м ...........................

 

 

22,20

Ширина В, м . ..........................................................

4,40

Высота борта Н, м ......................................

...

 

3,17

Водоизмещение D, т ..................................................

 

 

25,2

Скорость судна:

 

 

55 (15,3 м/с)

на волнении о0, к м / ч ......................................

 

 

с учетом встречного волнения vc, м/с . . .

19,6

М атери ал .....................................................................

 

Дюралюминий Д16

Модуль упругости £ , кгс/см2 . . . . ...............

0,7.10е

Моменты инерции сечений I и интенсивности

масс

т для первого тона,

[с учетом присоединенных к крыльям масс воды по формуле (5.16) ] даны в столб­

цах

2 —3 табл. 2 0 .

 

 

 

 

 

 

 

 

 

 

 

 

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

изгибе от колебаний

=

60

см2.

 

 

 

 

I =

3,5 м; она захва­

 

Прежде всего была определена длина удара волны

тывает шпации 0—3.

Затем по методу А. С. Повицкого

[формулы

(8.19) —

— (8.28) ] был получен график изменения силы

удара во времени, с

которого

снята наибольшая сила удара Рынке = 6140 кге и

время

ее нарастания tH—

= 0,070 с.

 

 

 

 

 

равно

 

 

 

 

 

Полное время удара по формуле (8.24)

 

 

 

 

 

 

 

T =

- i - =

l ^ -

= 0,18с,

 

 

 

 

 

 

 

 

 

fc

19,6

 

 

 

 

 

 

а время спада нагрузки

 

*сп = 7 — fH= 0 , l l c.

 

 

 

 

 

Закон изменения силы удара во времени взят по формулам (8.31).

 

 

Значения

силы удара волны во времени

 

 

.

Время /, с .............

0,00

 

0,01

0,02

0,03

0,04

0,05

0,06

0,07 0,08 0,09

 

Сила удара Р (f),

кге

, 0

 

880

1766

2640

3520 4400

5270

6140 5550 4860

 

Время t. с . . . .

.

. 0,10

0,11

0,12

0,13

0,14

0,15

0,16

0,17 0,18

 

 

Сила удара Р (0,

кге . 4300

3600

2940

2320

1750 1270

1040

370 О

 

* Или по более точным, например,

h2fk = j^ ( fk+2+ 16/t+i —30/*+ 16*_i—fk—2);

A3/* = ~ " ( — /fc+3 + 8/ft+2— 13/*+ i+ 13/*_1— 8/fc—2+ /*!-з);

8

hi' — ( — Ik+&+ 8 / k-\-1—8 /fe—1+ Ik—2).

292

Таблица 20

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

Моменты

Интенсив­

Функция

 

ность масс

формы

m/f.

сечений *

инерции

первого

первого

 

сечений, см*

тона т ,

тона fi

кгс-с3/**3

 

 

кгс*с3/м3

 

 

0

2 , 2 0 -1 0 °

41

 

—0,8700

31,2

1

3,20-106

244

 

—0,6968

118,6

2

4,10-106

369

 

—0,5254

1 0 1 , 8

3

4,80-10°

2 0 0

 

—0,3580

25,6

4

5,30-10°

104

 

—0,1981

4,0

5

5,40-10е

105

 

—0,0513

0,3

6

5,40-10°

81

 

0,0779

0,5

7

5,40-10°

83

 

0,1820

2,7

8

5,40-10°

84

 

0,2610

5,7

9

5,40-10°

79

 

0,3188

7,7

1 0

5,33-10°

80

 

0,3262

8,5

1 1

5,33-10°

97

 

0,3062

9,1

1 2

5,29-10°

116

 

0,2532

7,4

13

5,20-10°

150

 

0,1648

4,1

14

5,10-10°

166

 

0,0428

0,3

15

5,00-10°

150

 

—0,1068

1,7

16

4,80-10°

163

 

—0,2843

13,2

17

4,60-10°

147

 

—0,4806

34,0

18

4,30-10°

81

 

—0,6931

38,8

19

3,50-10°

231

 

—0,9446

205,0

2 0

2,50-10°

204

 

—1,1396

264,5

 

 

 

 

 

2 =737,9**

* Теоретическая шпация

)авна 1 1 1

см.

 

 

** С поправкой на полусумму крайних слагаемых.

 

Расчет

частоты свободных

колебаний

первого тона

и формы свободных

колебаний f x выполнен по методу конечных разностей (см. § 2 1 ) с учетом сдвига по двадцати ординатам.

Для определения частоты и формы колебаний f t первого тона составляют систему уравнений (5.62). Частоты находят из условия равенства нулю опреде­

лителя этой системы. Низший корень соответствующего уравнения = 57,2 с- 1 подставляют в систему уравнений (5.62), из которой находят формы колебаний первого тона. Ординаты функции формы приведены в табл. 20 и показаны на

рис. 65.

Обобщенные координаты определяют из уравнения (8.41)

% + 2Р А + *i4>i = •

Для определения обобщенной массы М х по формуле (8.42) в той же табл. 20 подсчитана сумма2 mf\ — 737,9 кгс-с3 -м—2 = 0,07379 кгс-с2 -см-2 .

293

Следовательно*

М г = A L2m f\ = 111-0,07379= 8,21 кгс-с^см- 1 [ 8 т]..

Коэффициент сопротивления взят равным р х = 0,07 Хх = 0,07-57,2 « 4 с- 1 . Обобщенная сила по формуле (8.33) равна

Фxi= f ^ r - j A ( l ) d 6 *

2/i =

^ ^ - 111-1,8362 = 0,917Р(0,

21

о

21

2-35U

где I = 350 см,

а 2 /х(|) =

1,8362 найдена следующим способом:

 

 

з

/х |

 

 

Подсчет 2

 

 

о

 

 

№ с е ч

е н и й

 

 

0

 

0,8700

 

1 ..................................

'

0,6968

 

2

0,5254

 

3

 

0,3580

 

Поправка

2,4502

 

0,6140

 

Исправленная с у м м а ......................

1,8362

f(D

0 , 5 -

- 1,0 -

Рис. 65. Форма колебаний первого тона теплохода «Ракета»

С численными коэффициентами уравнение (8.41) примет вид

 

|

„ , 1

,

0 2 .

 

0,917Р (0

-

 

ф +

2-4ф +

57,22ф =

----------- —

 

 

 

 

 

 

 

8,21

 

Обобщенную координату фх определяют по формуле (8.43), взяв At = 0,01 с.

ф/ — — ------ (

sin /.tj ^

Ф $еР‘5cos kts — cos kt; 2

sin kts ) ;

2 kM

 

s=o

 

 

 

s=o

0,01 -0,917-e 4ti

sin kt;

 

Ф (ts) e s cos kts

Ф7 =

2-57,2-8,21

2

 

V

 

s=o

 

 

 

c o s ^ ^ J Ф (

i

' s

)sesin kts

;

 

 

 

 

s—0

 

 

 

 

ф/ = 9,7 • 10~6e— 41,

sin Xtj 2

P (ts) e ts cos kts — cos ktj 2 P (ts)e*ts sin kts

 

 

s =

0

 

 

 

 

s

294

По этой формуле составлена табл. 2!.* Наибольшее значение обобщенной координаты оказалось равным: ф = 0,267 см.

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица

21

 

Определение обобщенной

координаты ф по формуле (8.38)

 

 

 

 

 

 

в расчете удара

волны в корпус судна

 

 

 

 

U с

Р (/). тс

еи

 

Р (О е4/,

и

 

COSЛ. t

s i n kt

 

[4] [6], тс

 

 

 

 

 

 

ТС

 

 

 

 

 

 

 

 

 

 

1

2

 

3

 

4

 

5

 

6

 

7

 

 

8

 

0 , 0 0

0

 

 

1

 

0

 

0

 

1 , 0 0 0

0

 

 

0

 

0 , 0 1

0,880

 

1,041

0,915

0,572

0,841

0,541

 

0,770

0 , 0 2

1,766

 

1,083

1,900

1,414

0,414

0,910

 

0,785

0,03

2,640

 

1,127

2,980

1,720

—0,149

0,989

—0,445

0,04

3,520

 

1,173

4,130

2,280

—0,651

0,759

—2,700

0,05

4,400

 

1 , 2 2 1

 

5,370

2,860

—0,961

0,278

—5,160

0,06

5,270

 

1,271

6,700

3,430

—0,959

—0,284

—6,420

0,07

6,140

 

1,323

8 , 1 0 0

4,000

—0,654

—0,757

—5,300

0,08

5,550

 

1,377

7,650

4,580

—0,132

-0,991

—1,015

0,09

4,860

1,433

6,960

5,150

0,423

—0,906

 

2,955

0 , 1 0

4,300

 

1,492

6,400

5,720

0,845

—0,534

 

5,430

 

 

 

 

 

 

 

 

 

 

 

Продолжение табл.

21

t, с

2 [8], тс

[7]-[9],

тс

[4]-[7], тс

2

[11]. тс

[6].[12], тс

[юд—ИЗ].

С -

4

* 1| ) .

с м

 

 

 

 

 

 

 

 

 

 

 

 

тс

 

 

 

 

i

9

 

10

 

 

11

 

12

 

13

 

14

15

 

16

0 , 0 0

0

 

0

 

 

0

 

0

 

0.

 

0

1

 

0 , 0 0 0

0 , 0 1

0,770

 

0,417

 

0,496

 

0,496

 

0,417

 

0 , 0 0 0

0,961

0 , 0 0 0

0 , 0 2

2,325

 

2 , 1 2 0

 

1,736

 

2,727

 

1,130

 

0,990

0,923

0,008

0,03

2,665

 

2,640

 

2,350

 

7,412

—1,104

 

3,744

0,887

0,032

0,04

—0,480

—0,365

 

3,140

 

13,502

—8,800

 

8,435

0,852

0,073

0,05

—8,340

 

—2,320

 

1,495

 

18,138

17,500

16,180

0,819

0,129

0,06

19,920

5,650

1,900

 

17,733

—17,000

22,665

0,787

0,174

0,07

—30,940

23,400

—6,150

 

9,683

—6,340

29,740

0,756

0,218

0,08

—37,255

37,000

 

—7,600

—4,067

 

0,535

36,465

0,726

0,258

0,09

—35,315

32,000

 

—6,320

—17,987

—7,550

39,55

0,698

0,267

0 , 1 0

—26,930

14,350

 

—3,430

—27,737

—14,800

29,15

0,670

0,190

Определив

по

формуле

(5.55)

значение

 

второй

производной

 

на

9-м

шпангоуте

/9 = ^ - ( /8 - ^ + и = 7П(°’2610- -2-0,3188 + 0,3262)= — 4,08-Ю- 6

* Расчет может быть выполнен также на ЭЦВМ; например, по программе

ГИИВТ

Д-1.

«Интегрирование дифференциального уравнения М’ф + 7?ф +

+ Nф =

Р (0

при М, R, N постоянных» (для машины «Проминь»), 1970.

295

найдем искомый изгибающий момент. По формуле (8.17) он равен

М нзг = £ //"ф = — 7-105• 5,4 • 106.4,08-10“ 6 0,267 = — 41,3-105' кгс-см =

= — 41,3 тс-м.

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

§ 31.

КАТЕРА

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

небречь искривлением корпуса судна, учитываемым в формуле (8.1) СО

членом 2 //'Фу-

j=i

Примем, что корпус судна перемещается как твердое тело по урав­

нению

 

w = z ( t ) + ( j - - ^ 1<р(/).

(8.46)

Силу удара волны Р (t) определяем по методу А. С. Повицкого, как это было предложено Ю. А. Шиманским [24 ], т. е. по формулам (8.19) — (8.28), учитывающим обводы судна. Из графика изменения силы удара во времени (см. рис. 63) можно определить Р макс при t = t„. Закон изменения силы удара по длине примем синусоидальным, т. е.

р(х, t) = p(t) s in^ y,

где I — длина района удара,

ориентировочно равная для

катеров

трем теоретическим шпациям (0,15 L);

 

I — абсцисса, отсчитываемая от начала или конца района удара.

В дальнейшем понадобится лишь максимальное значение р (t),

поэтому проинтегрируем р (£, t)

 

по длине удара и приравняем макси­

мальной силе удара

 

 

 

 

Рмакс J0

 

1

Рмакс>

 

откуда

 

_я Р макс

 

 

 

 

(8.47)

Рмакс

2/

 

и, следовательно,

Рмакс Sin

 

 

Р (X) =

.

(8.48)

За пределами длины удара интенсивность удара р (*) равна нулю. Элемент длины корпуса dx в процессе движения будет нагружен си­ лой инерции, которая равна массе элемента mdx, умноженной на. ус­ корение, и направлена обратно ускорению

\d2w ,

т -----dx, dt2

296

и заданными ударными волновыми усилиями pdx. Интенсивность силы удара применительно к формуле (8.43) можно считать заданной функ­ цией х. Поэтому суммарная интенсивность распределенной нагрузки от удара волны без учета сил веса и уравновешивающих их сил под­ держания равна [учитывая (8.46) ]

Соответственно срезывающие силы и изгибающие моменты равны

X X

М , „ = Я \р ~ т [*“* ( т ~■ т ) ^ ] 1 '

0

0

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

Значения входящих в формулы (8.49) и ,(8.50) производных г и ф могут быть определены из условий динамического равновесия, выра­ жающих равенство нулю срезывающей силы и изгибающего момента на кормовом перпендикуляре (Vn = Ш изг]я = 0).

пп

1=0 1=0

(8.51)

Где п — номер теоретического шпангоута, на кормовом перпенди­

куляре.

Из уравнений (8.51) путем совместного их решения легко опреде­ лить максимальные ускорения z и ф, причем максимальные их значе­ ния получаются, если в формуле (8.48) рмакс взято по формуле (8.47), Зная максимальные значения г и ф и подставив их значения в фор­

*Множитель Vg при каждом суммировании введен ввиду так называемого «интегрального» суммирования, дающего удвоение.

**Третья сумма исчезает, так как предполагается, что центр тяжести судна

находится на миделе.

297.

мулы (8.49) и (8.50), можно найти срезывающую силу и изгибающий момент в любом сечении судна.

Расчет укладывается в простую табличную форму (табл. 22). Для определения перерезывающих сил и изгибающих моментов

достаточно из уравнения (8.51) найти максимальные значения уско­

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

В этом случае следует при вычислении интенсивности силы удара брать по выражению

Р (t) = P «)-$-,

(8-52)

заимствуя Р (/) из закона (8.31) или графика рис. 63, оставляя по длине удара прежнюю зависимость, т. е.

P— P(t) sin

sin •

(8.53)

Правые части уравнений (8.51) становятся функциями времени. Впрочем, учитывая кратковременность волнового удара, удобнее вер­ тикальную и угловую скорости корпуса в конце удара получить по суммарному ударному импульсу и его моменту.

§ 32

ВОДОИЗМЕЩАЮЩИЕ ГИБКИЕ СУДА

У некоторых гибких судов частота волнения близка к частоте сво­ бодных колебаний корпуса (обычно первого тона) и, следовательно, нагрузка от волнения носит динамический характер. К таким судам относятся большегрузные гибкие речные суда, например грузовые теплоходы «Волго-Дон».

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

Для расчета качки по А. Н. Крылову необходимо определить кру­

говую кажущуюся частоту волнения (с-1; + при встречном волнении,

— при попутном)

( 0 = 2 llU ± 1,25 К О '

(8.54)

где v0 — скорость судна, м/с; Яв — длина волны, м.

298

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