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

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

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

Зк" % ФА+ ( -

3* л / ^

+

3*/‘ ^ ) ф / -

 

^ lki t

+ ^Л/^Л

^ ^ Л/

^

^

= —M•

(3-88)

Уравнение типа

(3.87) составляется для

тех сечений, в

которых

неизвестны прогибы,

а (3.88) — в которых неизвестны углы поворота.

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

Уравнениями (3.87), (3.88) можно пользоваться и для крайних

опор,

полагая жесткости k hj несуществующих элементов (участков

балки)

равными нулю.

Особенно простой вид уравнения имеют для призматических балок при разделении их на элементы одинаковой длины

I 0 l = I l 2 = . . . = 1 0 = I И I q i = / j 2 == • • • = l o — I*

а именно:

 

—6фЛ+ (12—Ф) Ф/—6фА—3{5ft+ 3pft— Р1\

(3.87а)

Зфд—Зф* + Р/, + 4Р; + рй = —М,

(3.88а)

где через Ф обозначено

l3M а 2

(3.89)

2EI

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

П р и м е р .

Определить вынужденные поперечные

колебания

призмати­

 

ческой свободно

опертой

балки

под

действием

возмущающей

силы Р cos at

 

= 0,225 кН,

ш =

18 с-1 ), приложенной в середине пролета (рис.

33). Масса

 

балки равна 0,086 т, длина 1 0 м, момент инерции поперечного сечения 180-10~8м4,

 

модуль Е = 1,96 -108

кН/м2.

 

 

 

 

 

 

 

 

 

Разбиваем балку на десять элементов длиной I — 1 м каждый и сосредото­

 

чиваем в разделительных сечениях массы Л1 =

0,0086 т.

 

 

 

 

 

Входящий в группу уравнений (3.87а) отвлеченный коэффициент Ф по фор­

 

муле (3.89) равен

 

13-0,0086-1 8 2

 

 

 

 

 

 

 

 

 

Ф =

 

= 0,004.

 

 

 

 

 

 

2-1,96-Ю8- 180-10~8

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнения типа (3.87а) составляем для 1, 2, 3, 4 и 5-го сечений, учитывая,

 

что ф0 = Р6 = 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

1 1 , 9 9 6 ф ! —

6 ф 2

З Р 0

+

З Р 2

=

0 ,

 

 

'—

6 ф

2 - ( -

1 1 ,9 б9 ф6 зф2 З

р

-4j- З р з0 :

 

 

 

 

6 ф

2 +

1 1 , 9 9 6 ф 3 —

6 ф 4

З Р 2

+

З Р 4 =

0 ,

 

 

 

 

 

 

 

 

 

 

 

155

 

6<p3 + 11,996ф4 — 6ф5 — Зрз = О,

6ф4+ 11,996ф5 — 6фв — ЗР4 + Зрв = 0,225-1.

Далее составляем пять уравнений типа (3.88а) для 0, 1, 2, 3 а 4 сечений

- З ф1 + 2Ро + Рт = 0;

— Зф2 + Ро + 4Pi + Рг = 0;

Зф| — Зфз —|—Pi -f- 4Р2 -f- Рз = 0;

Зф2 — Зф4 + Р2 + 4рз + Р4 = 0;

Зф3 — Зф5 + Рз + 4р4 = 0.

P c o s w t

Рис. 33. К применению метода конечных элементов

Если в пятом уравнении положить Рв = — Р4 и фв = ф4 и поделить это уравнение на два, то обе группы уравнений составят систему десяти канониче­ ских уравнений относительно десяти неизвестных.

Корни системы * и амплитуды перемещений, найденные по формулам (3.85), равны

Ф1 =

14,995

I2-14,995

2-1,96-108-180-10~1

 

1

Ф2 =

28,569

f2 = 40 мм;

Фз =

39,419

=

56 мм;

ф4 =

46,462

/ 4

=

6 6 мм;

ф5 =

48,923

/ 5

=

69 мм;

Ро =

15,234;

Р4 =

14,515;

Р2=

12,413;

Рз =

9,103;

Р4 =

4,852.

 

 

 

Форма изгибных колебаний показана на рис. 33 внизу. Расчет этой же балки точным методом (с. 131) дал такой же результат.

* Найденные на ЭЦВМ.

156

Г Л А В А

4

СИЛЫ, ВЫЗЫВАЮЩИЕ ВИБРАЦИЮ

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

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

§ 16

УСИЛИЯ, ВЫЗЫВАЕМЫЕ ПОРШНЕВЫМИ МЕХАНИЗМАМИ

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

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

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

М г — массу прямолинейно движущихся частей — поршня, пол­ зуна, штока-поршня;

М 2 l ~ s---- часть массы шатуна;* s — расстояние центра тяжести

1шатуна от верхнего (головного) подшипника, т. е. от точки А, I — длина шатуна.

Вточке В будем считать сосредоточенными величины:

М 2 ----- оставшуюся часть массы шатуна;

*

Масса

шатуна

М 2 распределена между точками

А к В таким

образом,

чтобы

общий

центр

тяжести масс,

сосредоточенных в

этих точках,

совпадал

с действительным центром тяжести

шатуна.

 

 

157

М 3 — массу мотылевого пальца и приведенную к радиусу мо­ тыля массу щек коленчатого вала.*

Даламберова сила инерции массы, сосредоточенной в точке В, как центробежная сила направлена по линии ОВ и равна

Р = г О)2

( м а - ^ + М 8) ,

(4.1)

где и = 2пп --— угловая

скорость вращения

вала, с 1;

радиус мотыля, м.

Массы, сосредоточенные в точке А и движущиеся только поступа­ тельно по оси z, дадут вертикальную силу инерции, равную взятому со знаком минус произведению массы на ускорение

Q=—[ м г - \ - М г —у- db

(а)

 

 

dt2

 

где z — координата точки А, равная

 

 

г = ОС+ СА = ОС+ ] / > — ВС2;

 

2 = Г COS at 4- ]/"

Рr2sin2to/.

 

Учитывая, что второе слагаемое под корнем много меньше пер­

вого,** можно, разлагая корень в биноминальный ряд

 

2= г cos at + 11 1---- — sin2 at

 

 

 

V

12

 

 

z = r cos со/ -4-/(1

---- — sin2 at — ^—sin4 a t — . . . ) ,

(б)

\

2/2

81*

)

 

ограничиться лишь членами, содержащими rtl в степени не выше вто­ рой

 

2= г cos at + l ( \ ---- — sin2

.

 

V

2/2

 

)

Произведя замену

 

 

 

 

sin2 at = —----- cos.2co/,

 

 

2

2

 

 

получим

r2

 

 

 

 

 

-— cos2o>/.

 

2= / -------- 1- r cos at

 

 

41

 

41

 

* Приведение выполняется на основании

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

масс относительно оси вала. Так, если масса щек равна Мщ и их центр тяжести

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

М щ, где г —ра­

диус мотыля.

 

^•2

 

** Обычно rtl меньше 1/4, а следовательно, — < 1/16.

 

158

Тогда ускорение

равно

 

 

 

d 2z

2

/

/"2а>2

о ,

— = —гаг cos ш

------- cos 2ш.

d t 2

 

I

 

Подставив это ускорение в (а), получим

Q = п о2 [ м г+ М2 ~~y~J ^cos + -у- cos 2co/j .

(4.2)

Инерционные силы Р и Q, определяемые формулами (4.1) и (4.2), вызовут ряд внутренних усилий взаимодействия между движущимися деталями кривошипного механизма, а также между движущимися

Рис. 34. Инерционные и внутренние усилия в одноцилиндровом дви­ гателе

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

Сила инерции Р растягивает мотыль и вызывает силу Р х, прикла­ дываемую рамовым подшипником к мотылю; точно такая же сила Р 2, но обратного направления, действует на раму двигателя (рис. 34, в).

Вертикальная сила инерции поршня * и прямолинейно движу­ щихся частей Q уравновешивается горизонтальной реакцией стенки цилиндра Нг и усилием R lt прикладываемым к поршню шатуном. Растягивающие шатун усилия R 2, R3 передаются мотылю (сила /?4) и уравновешиваются реакцией рамового подшипника Rb и моментом сил инерции вращающихся масс т (вала, маховика, других мотылей...) Этот момент сил инерции появляется вследствие того, что сила Rt со­ здает в данном положении кривошипно-шатунного механизма тормо­ зящий момент, т. е. слегка замедляет вращение вала. Сила Re, обрат­ ная реакции R3, приложена к раме двигателя.

* При малых углах ш/ (вблизи верхней мертвой точки) ускорение отрица­ тельно, т. е. направлено вниз, а сила инерции Q — вверх.

159

Итак, в рассматриваемом положении кривошипный механизм при­ лагает к станине и цилиндру двигателя, а через них и к корпусу судна, три силы (см. рис. 34, в)

H2 = Qtgp;

Р2 — га2 (^М2 -у- + M3j ;

d _

Q

6

co s р

где р — угол шатуна с вертикалью (с осью г).

Разлагая силы Р 2 и Re по осям координат (рис. 34, г) и группируя

силы, получим

 

 

 

Y — Р2sin соt = ra2[M2 -у- +

M3j sin at;

(4.3)

Z — P2cos at = rco2 \m 2 y - +

Mgj cos at;

 

Q = Recos P = rco2

+ M2

j cos <°t +

 

+ -r f - -+- M2 l~ s j cos 2at;

H3 = Resin p = Qtg p.

Горизонтальное усилие Y совершает свой полный цикл за один оборот вала двигателя; будем называть такие силы силами первого порядка.

Сложив вертикальные силы, получим

Z + Q = гсо2 (Мх-f М 2+ М9) cos at +

( X + М2- t z l j cos 2 a t.

Они состоят из сил первого и второго порядков

 

Z1 = rco2(M1 + i'H2 +

i,H3)cos(oC,

(4.4)

Z2 = ~~~~

+ М2 ~

~ J cos 2со/,

(4.5)

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

Горизонтальные усилия Н 2 и Н3 образуют опрокидывающую пару с моментом (см. рис. 34)

т = Н -0А.

Можно показать, * разлагая т в ряд, что главная его составляющая имеет второй порядок.

* См. В. В. Д а в ы д о в , Н. В. М а т т е с. Динамические расчеты проч­ ности судовых конструкций. М., «Транспорт», 1965, с. 243.

160

Итак, движущиеся части поршневого двигателя прилагают к не­ подвижным его частям (цилиндру, раме), а через них и к фундаменту двигателя, т. е. корпусу судна, горизонтальные инерционные силы первого порядка (4.3), вертикальные силы первых двух порядков формулы (4.4) и (4.5), а также опрокидывающие моменты.

Есть в машине инерционные усилия и высших порядков, но ампли­ туды их незначительны.

Инерционные усилия в многоцилиндровом двигателе. Инерционные усилия в многоцилиндровом двигателе, определяемые формулами (4.3) — (4.5), можно считать перенесенными на ось вала. Инерцион­ ные усилия в многоцилиндровом двигателе действуют в различных плоскостях и изменяются во времени с некоторыми сдвигами по фа­ зам, зависящими от углов установки мотылей коленчатого вала от­ носительно друг друга. Условимся отсчитывать углы установки мо­

(г,+г2)1

Рис. 35. Схема инерционных усилий в многоцилин­ дровом двигателе

тылей от мотыля одного из крайних цилиндров в направлении враще­ ния вала (угол а£ на рис. 35) и введем для краткости следующие обо­ значения:

а — г(о2 \М2-y -+ M 3j ;

Ь = п о 2 ( М х + М г + М з ) ;

(4.6)

Будем ставить у величин а, Ь, с индексы соответственно номеру цилиндра. Таким образом, будем считать, что в районе первого ци­ линдра (в точке О на рис. 35) приложены силы

агsin со/, bicosat,

в районе второго цилиндра — силы

a2sin(<o/ + a 2),' b2cos + a 2), с2cos 2 (со/ + а 2)

и, вообще, в районе любого цилиндра (рис. 35)

a/ sin(d)/ + a I), 6,cos(co/ + a (), с, cos 2 (toA+a,).

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

161

Если в результате подсчета выяснится, что неуравновешены силы

и моменты, друг другу соответствующие,

т. е. действующие в одних

и тех же

плоскостях с одинаковыми частотами (например, Y и Мг,

Zx и М

и т. д.) и фазами (уу = уМг,

= yMi . . .), то можно попы­

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

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

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

ax = аг= . . . = а\ b1 = b2= .. . = Ъ\

С\ = с2 = • ■• ~ с

при одинаковом расстоянии h между осями цилиндров приведенные формулы значительно упрощаются. В этом случае можно руководст­ воваться готовыми решениями (табл. 12).

В этой таблице приведены коэффициенты х и ц в выражениях ам­ плитуд соответствующих усилий

Y -

Мг= p,jah\

Zl — Klb

M yi =

Z2 = x2c

My = ^ch .

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

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

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

164

Т а б л и ц а 12

Коэффициенты неуравновешенных инерционных сил и моментов в двигателях [коэффициенты и и р к формулам (4.12)]

Расположение мотылей, число цилиндров и тактов

1 -цилиндровый

Неуравнове­

Неуравнове­

шенные

шенные силы

моменты

 

 

 

!

первого порядка

второго порядка х3

первого порядка Hi

второго порядка Из

1

1

0

0

2 -цилиндровый двухтактный

0

2

1

0

i

3-цилиндровый

двух- и четырех­

 

 

 

 

0

0

1,73

1,73

тактный

3 '

2*3т^

 

 

 

 

4 -цилиндровый четырехтактный

0

4

0

0

 

 

 

 

1

 

 

 

 

 

4 -цилиндровый двухтактный

4------- J

0

0

1,41

4

 

 

 

 

2

 

 

 

 

 

4 -цилиндровый двухтактный

 

 

1

0

0

2,83

2

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

1

к

 

 

 

 

 

6 -цилиндровый

четырехтактный

3 ,4 '

 

 

 

0

0

0

0

 

 

 

 

 

 

 

 

 

6 -цилиндровый двухтактный

к

 

и

 

0

0

2

6,93

 

 

 

 

6 -цилиндровый двухтактный

 

 

 

 

0

0

0

6,93

6 -цилиндровый двухтактный

К

 

I

л

0

0

0

3,46

 

 

 

 

 

 

"

I

T

'

 

 

 

 

165

Расположение мотылей, число цилиндров и тактов

8-цилиндровый четырехтактный

 

А /

4,512,7

 

 

'3,5

8-цилиндровый двухтактный

А

 

т —X4~3

 

 

Г '

8-цилиндровый двухтактный

Л

в

4 -3

 

5 '

,'а

 

 

г

8-цилиндровый двухтактный

ч

/ 7

 

4 -*

 

 

 

 

4 s

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

Неуравнове­

Неуравнове­

шенные

шенные силы

моменты

первого порядка

1 порядка

я«

первого порядка 1-М

второго порядка Ма

 

второго

 

 

 

0

0

 

0

0

0

0

 

2,61

5,66

0

0

 

1,08

5,66

0

0

 

2,16

2,83

нием газа в двигателе внутреннего сгорания или пара в паровой ма­ шине.

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

На рис. 36 показана схема передачи усилий, возникающих во время рабочего хода в одноцилиндровом двигателе, который установлен в корпусе винтового судна. Давление газа (или пара) в цилиндре че­ рез поршень передается шатуну АВ, мотылю ВО и создает крутящий момент на коленчатом и гребном вале OD. При вращении вала по ча­ совой стрелке, если смотреть с носа в корму (т. е. отточки О к точке D), со стороны воды будет прилагаться к винту помимо толкающего про­ дольного усилия момент М, направленный в сторону, противополож­ ную вращению винта. Этот момент, внешний по отношению к судну с работающим в нем двигателем, вызовет крен судна * и сопровождаю­ щий его восстанавливающий момент со стороны воды М х.

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

* Ничтожный по величине и никогда не учитываемый.

166

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