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

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

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

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

Фх, Фг. • • • . фт ;

бхфх, M b . . . . £т Фт>

где kv k 2, . . . , km — любые числа.

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

Вопрос о разыскании главных координат, физически присущих системе, математически сводится к замене выбранных более или ме­ нее произвольно координат ф1; <р2, . . . , <рт другими координатами, превращающими два однородных полинома второй степени (2.5) и (2.8), т. е. кинетическую и потенциальную энергии, в суммы, не содер­ жащие произведений. Эта операция всегда возможна, однако по тру­ доемкости эквивалентна решению задачи в обобщенных, а не в глав­ ных координатах. Поэтому, если не удастся сразу подметить, что, должны представлять собой главные координаты, решение выполняют в обобщенных. При этом решение в обобщенных координатах неиз­ бежно приведет к решению в главных координатах, из которых только и могут быть составлены всякие другие координаты, как это видно из выражений (2.12а), в которых первые слагаемые в правых частях пропорциональны первой главной координате, вторые — второй и т. д. Вскрыть физический смысл какой-нибудь главной координаты можно, полагая все главные координаты, кроме одной, равными нулю или выбирая начальные условия так, чтобы в выражении (2.12а) все сла­ гаемые, кроме соответствующих интересующей нас главной коорди­ нате, отсутствовали.

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

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

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

можны

лишь

вертикальные перемещения масс. Модули упругости Е —

=2-10®

кН/м2

[~ 2 -10е кгс/см2] и G — 8,1-10? кН/м2 [— 8,1 -10б. кгс/см2].

Массами

балки

и пружины пренебречь.

80

Н а х о д и м к о эф ф и ц и е н т ж е с т к о с т и б а л к и ( к Н , м , с)

 

 

Nt =

48EI

48-2-108-2,92-10—4

1500 кН/м

150 тс/м]

Р

 

12,3?

 

 

 

=

и пружины

 

 

 

 

 

 

 

N2 =

d*G

0,034-8,1 • 107 50

кН/м

5

тс/м].

 

 

64nr3

64-28-0,09s

 

 

 

Если принять за обобщенные координаты величину опускания большой массы (фх) и величину удлинения пружины (<р2), то потенциальная энергия будет равна

u = - L n ^ + ± n24 .

Рис.

19. Система с двумя степенями свободы (к примерам на с. 80, 85, 149 и 152)

Кинетическая энергия системы

 

 

 

 

к = \

М1Ф1 + - L Л43(ф1 + ф2)2 =

- L

(Л4, + м 2) ф2 +

Л42ф1ф2 + - L м 2ф2.

Уравнение Лагранжа (2.7)

 

 

 

 

 

 

(Mi + Л4а) ф! + Л4афа"+ Ntф! = 0;

 

 

 

 

Л4аф1 + МаФа +

Л?2фа = 0.

 

 

Подстановка решений в виде гармонических функций

 

ф1 =

at cos (Xt + а)

и

ф2 = а2 cos (Xt +

а)

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

дает

 

 

 

 

 

[У, — (Мх + М г) Я2] в! — М2Я2а2 = 0;

 

 

 

M 2X2at + (N2 — М 2Х2) а2 = 0

 

 

или в численном виде

 

 

 

 

 

 

 

 

(1500 — 12Я2) at — 2Х2а2= 0;

 

 

 

— 2X2at + (50 —2Я2) аг = 0.

 

(з )

Составив определитель

частот (2.11)

 

 

 

 

1500 — 12А.2

— 2Я2

20Я4 — 3600Я2 +

75 000 = 0

 

— 2Я2

 

50— 2Я2

 

 

 

 

 

 

и решив его, найдем Я2 =

24,0455 с

2 и Я2 = 155,9545

с

2.

81

Д а л е е с л е д у е т п о д с т а в и т ь

Хр а за те м Х | в у р а в н е н и я (з) и н а й т и все a-k

(1500 — 12-24,0455)

аи — 2-24,0455-ап =

0

 

— 2 - 24,0455 • ап + (50

— 2 ■24,0455) а21= 0

 

 

(1500 — 12-155,9545) а12 — 2 ■155,9545 • а22 =

0

 

— 2 • 155,9545 • ап + (50 —

2 • 155,9545) а22 = 0

 

 

Следовательно, согласно формуле (2.12)

 

 

Ф1

ап cos (М + a i) — 0,8397а22 cos (Х2/ +

a 2);

ф2 =

25,19an cos (X]/ + a i) + агг cos (M +

а г)*

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

амплитуды

колебаний малой нижней массы в ее абсолютном движении в 25 +

+ 1 = 26

раз больше амплитуд колебаний большей массы,* либо колебания

высокой частоты Х2, когда массы сдвигаются в противоположные стороны от равновесного состояния, причем амплитуды колебаний малой массы в ее абсо­ лютном движении 1 — 0,84 = 0,16 много меньше амплитуд колебаний большой массы 0,84.**

Каждое из описанных колебаний является главным.

§ 7

ВЫНУЖДЕННЫЕ КОЛЕБАНИЯ

Общее решение. Пусть в системе с т степенями свободы в какойлибо точке ее приложена гармоническая возмущающая сила F cos со/.

Давая согласно формуле (2.4) последовательно каждой из обоб­ щенных координат приращение бф;- и определяя работу, произведен­ ную силой F cos со/ на соответствующем перемещении, можно опреде­

лить все обобщенные силы Ф 1(

Ф2, . . . , Фт . Каждая из этих обоб­

щенных сил, очевидно, должна быть пропорциональна cos со/

 

Фх =

Р! cos со/;

 

0 2=

P 2cos со/;

 

Ф;=

Р,-cos со/;

( 2 .21)

 

Фт =

Ртcos со/.

 

Величины P j зависят от амплитуды F возмущающей силы, точки

ееприложения и выбора обобщенных координат.

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

движна

%? — N 2 : М 2 = 25. Наличие большой

массы лишь незначительно сни­

жает частоту.

 

**

Квадрат частоты колебаний большой массы при отсутствии малой равен

N x : М]

= 150. Наличие малой массы лишь

незначительно увеличивает эту

частоту.

82

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

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

При учете сил сопротивления колебаниям удобно ввести функцию

рассеивания

т

 

т

 

У = 7Г 2

2 Я / аФ/Фа,

(2-22)

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

 

[ Ф / 1 с о п р = ---- (2.23)

 

 

(7ф/

 

В этом случае можно считать, что обобщенная сила,

входящая

в уравнение Лагранжа

(2.6), состоит из трех слагаемых

 

~

= [Ф/]упр + [Ф/lconp + Ф/.

(а)

Первая часть обобщенной силы обусловлена силами упругости и

учитывается по потенциальной энергии

 

[Ф /]упр

^ 7 >

(2.24)

вторая есть сила сопротивления

(2.23) и, наконец,

третья — часть

обобщенной силы, обусловленная возмущающими силами (2.21). Подставив (2.24) и (2.23) в (а), получим

_d_

д К .

д У

дП

(2.25)

dt

дф/

(Эфу

«Эфу

 

Подставив в это уравнение выражения для кинетической энергии К, функции рассеивания V и потенциальной энергии П по формулам

(2.5), (2.8) и (2.22), получим

т

 

2 (М/аФ* + Rjkф* + Nlk<pk) = Фу

(2 .26)

k = \

или, при гармонических возмущающих силах (2.21),

т

2 (М/аФа+ Rj№k + N/аФа) = Pycos at.

(2.27)

k —\

83

Уравнение (2.26) или (2.27) есть /-е уравнение системы совокуп­ ных дифференциальных уравнений второго порядка, в каждое из ко­ торых входят в общем случае все неизвестные функции ц>к и все их первые и вторые производные.

Входящая в уравнения (2.26) или (2.27) обобщенная сила сопро­ тивления (2.23) при решении конкретных задач обычно получается непосредственно, а не через производную от функции рассеивания. Функция рассеивания является определенно положительной и, как правило, не содержит слагаемых с произведениями скоростей (т. е. все Rik = 0 при / ф k).

Применение уравнений Лагранжа не является единственным спо­ собом составления дифференциальных уравнений движения. В неко­ торых случаях (особенно когда вычисление коэффициентов жесткости Njk и потенциальной энергии оказывается громоздким) с успехом мо­ жет быть использован принцип Даламбера (см. §8 и приведенный ниже на с. 92 пример).

Что касается характера сил сопротивления, то, как мы видели при рассмотрении системы с одной степенью свободы, внутреннее со­

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

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

к

сопротивлению,

пропорциональному скорости [формулы (г) на

с.

30 и (1.23)]. К коэффициентам внешнего сопротивления Rjk в этом

случае следует добавить величину, учитывающую внутреннее (гисте­ резисное) сопротивление

n I ^ikNjk

Решение уравнения (2.27) будем искать в виде частных интегралов

Ф* = akcos Ы + bksin со(,

(2.28)

где ak, bk — постоянные, подбираемые так, чтобы система уравнений (2.27) удовлетворялась тождественно. Подставив (2.28) в (2.27), по­ лучим, учитывая и гистерезисное сопротивление

т

 

2 \{Nik— М1к®2) ak+

(Rjk®+ KikNjk) bk]cos соt -f

m

 

+ 2 [—(Rjk®+ KjkNjk) ak+

(Njk— M,kCO2) bk]sin (at = P ;- cos соt.

k=i

 

Чтобы последнее выражение было справедливо для любого момента времени, необходимо

2 i(RjkMjk®2) akJr(RjkJr'XjkRjk) bk\ — Pj\

k=\

(2.29)

m

2 [( R j k ® + KjkRjk) a k + ( R j k

м jk со2) bk]= o.

k=i

 

Система 2m неоднородных алгебраических уравнений с неиз­ вестными * alt а 2, . . . , ат, Ь1г Ьг, . . . , Ьт позволяет найти все эти

Уравнения (2.29) представляют /-ю пару уравнений системы.

84

параметры, а следовательно, по формулам (2.28) все обобщенные ко­ ординаты.

П р и м е р . Определить

вынужденные колебания системы,

изображенной

на рис. 19, в предположении,

что на малую массу М 2 действует гармоническая

сила F cos a t с амплитудой F = 1 кН [102 кге) и частотой со =

5 с- 1 . Внешнее

сопротивление считать отсутствующим, а коэффициент неупругого сопротив­

ления взять равным

=

х = 0,2.

За обобщенные

координаты,

как и ранее (см. с. 81), возьмем q>x — опу­

скание массы М± и ср2 — увеличение расстояния между массами.

Давая последовательно координатам приращения 6q>i и 6<р2, найдем обоб­

щенные силы Фх — F cos at

и Ф2 = F cos at,

т. е. Рг = Р2— F.

Функция рассеивания в данной задаче, очевидно, равна

V<= _1_

xNj

kN2 ф|

2

а

а

и, следовательно, все 12 величин, входящих в систему (2.29), равны

Мц

М 12

 

М =

М 22

 

М 21

 

II

N 12

 

 

j

1 N2i

N 22

Я ц +

а

Ri2 ■

R =

 

 

 

^21 *^21

R 22 '

 

СО

 

М1 -)- м 2 м 2

12 2

 

 

м ,

 

м ,

2

■матрица масс*;

 

 

2

 

Ni

0

(I

1500

0

 

 

0

N 2

In

0

50

матрица жесткостей;

 

 

 

 

 

%N12

 

V.N-,

0

300

0

а

 

 

 

со

 

 

 

 

 

kN22

 

 

V.N.

 

10

а

 

 

 

а

 

со

- матрица сопротивлений.

 

 

сначала для j = 1, а затем для j = 2,

Составляя систему уравнений (2.29)

получим

 

 

 

 

 

. . . в общем виде

 

 

 

 

(Nn Nlaсо2) %

Rnabi

(Nj2 — Afj2cos) д2 -f- R i2ab2 — Pi',

(Nn — МцС0 2) bi

 

Ri2aci2

(Ni2M i2a 2) b2— 0;

(Nsi — Alj^co2) Q-i -f- R 2iabx

(N22 — Л122со2) q2-|- R22ab2 — P 2,

— /?2lcoai +

(N21 — M21co2) bx R 22aa2-f (N22— Afa2co2) b2 = 0,

. . . подставив буквенные символы

 

 

[TVi — (Mi +

M 2) со2] cii + xNibi M2a 2a2 = F;

k N i Qi +

[ iV (Mi

+

M 2) с о 2 ]b i

— Mc o2 26 2 = 0 ;

M 2a 2ai -f- (N2— M2co2) a2-(- y.N2b2 — F ;

•—• M2a 2bi кЫ2а2+

(N2M 2a 2) b2 — 0,

См. пример на с. 81.

85

. . . п о д с т а в и в ч и сл а

 

 

 

 

 

1200^

+ 300&1 — 50а2

 

 

=

1

300а!

— 1200Й!

+

5062

=

— 50а!

 

+

Ю62

=

1

50&1 -J- 10а2

- 0*

Решая полученную систему, найдем

 

аг = — 0,000889 м; bt = 0,003759 м; а2 = — 0,018769 м

и Ь2= 0,095553 м

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

 

= — 0,089 cos at + 0,376 sin at = 0,387 cos (at -f- 103°);

Ф2 = — 1,877 cos at + 9,555 sin at = 9,73 cos (at + 101°).

Итак, большая масса колеблется с амплитудой 0,387 см, а малая—с ампли­ тудой 9,73 см. Обе массы колеблются примерно в одной фазе (обе вниз, затем обе вверх и т. д.), т. е. «форма» колебаний соответствует колебаниям низкой частоты

(с. 82), что вполне естественно, так как частота возмущающей силы (со = 5 с-1 )

близка к низкой частоте собственных колебаний (Я-i = 1^24).

Ввиду близости к резонансу с низкой частотой сопротивление довольно сильно (в пять раз) снижает амплитуду. Без учета сил сопротивления мы полу­ чили бы

фх == — 2 cos at и ф2= — 50 cosat (см)**

§ 8,

ИЗГИБНЫЕ КОЛЕБАНИЯ БАЛОК С СОСРЕДОТОЧЕННЫМИ МАССАМИ

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

Выберем за обобщенные координаты перемещения масс вниз от положения равновесия фа, ф2, . . . , <р/, <-. . , фт .

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

денных колебаний, возмущающие силы

 

p k= p'k cos at + Pi sin со/.

(2.30)

*У второго и четвертого уравнений сменены знаки, чтобы система стала канонической.

**Амплитуда колебаний верхней массы равна 2 см, а нижней (поскольку

через ф2 было обозначено удлинение нижней пружины) —52 см.

86

Перемещение j-й массы под действием указанных сил равно

 

ф/ = Цв/* (Pk — Mkifk— RkVd-

(2.31)

 

*=1

 

 

 

В матричной форме эти уравнения имеют вид

 

или «

= * А (Р—Л1<р—-Р?),

(2.31а)

 

 

 

 

где

^4Жер + Л/?ср-f ^= .ДР,

(2.316)

 

 

 

 

8ц

612 .

• • 81т

 

 

^21

б22

■ • &2т

матрица податливости;

 

 

 

 

 

8mi

6m2

8тт

 

 

|

И,

•••

 

Mj

. /У*

 

 

w i

 

 

 

 

 

 

 

 

Рис. 20. Балка

с сосредоточенными массами

 

 

Ф1

 

 

Ф1.

 

Ф1

 

Ф =

Фа

ф =

Фа

;

Фа

 

;

 

9 =

 

 

Ф т

 

 

Ф т

 

Ф т

— матрицы-столбцы перемещений,

ускорений

и скоростей;

 

 

Pi

 

Р 1cos (0/ +

Р 1sin со^

 

 

Ра

 

Рг cos

+ P^sin со^

 

 

Рт

 

РтCOS Cirf + Pm sin (£>t

— матрица-столбец

возмущающих

сил,

из

которых в конкретных задачах

только одну обычно можно считать отличной от нуля;

87

М г

0 . . .

0

 

0

м й . . .

0

— диагональная матрица масс;*

м =

 

 

 

0

0 . . .

М т

 

R i

0 . . .

0

 

0

Я 3 . . .

0

— диагональная матрица

/? =

 

 

сопротивлений.

0

0 • • •

я т

 

В коэффициенты сопротивления Rk можно включить также и внут­ реннее сопротивление, и, как было установлено, внутреннее сопротив­ ление можно также при гармонических колебаниях считать пропор­ циональным скорости перемещения (см. с. 84).

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

ничных сил. Так,

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

под действием силы,

приложенной в k-й точке, равен

 

 

пн — Pt3

t

xk

xi

1 (.l-X k )2

если

хк^>х**

 

6EI

 

l

l

E

l2 .

 

 

 

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

 

 

 

 

 

 

 

-

В

l —xk

xi

1 {l-X k )2

где

В =

Е -

6ik ~

6 '

 

i

l

E

l2 .

 

 

EI

В случае,

если стержень имеет

переменное

сечение,

можно bod

 

 

 

 

 

 

М ]ЗГМ*ЗГ

 

пользоваться интегралами Мора ***

*/* =

£7

dx,

 

 

 

 

 

 

о

 

 

 

где Л4;изг, М Г — изгибающие моменты по длине балки от единичных усилий, приложенных в точках / и k.

Решение уравнения (2.31) в случае гармонических возмущающих сил (2.30) можно искать в форме

 

 

? = a cos со^ +

b sin at

(2.32)

*

В виде

матрицы-столбца ее представить нельзя, так как в этом случае

произведение

Му теряет смысл.

 

 

 

** Подробнее см. В. В.

Д а в ы д о в ,

Н. В.

М а т т е с, И.

Н. С и в е р -

ц е в.

Учебный справочник

по прочности судов

внутреннего

плавания. М.,

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

*** То же, с. 166.

88

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

Фх = ах cos соt + Ьгsin a t; ф2= а2cos со^ + b2sin со/;

(2.32a)

Фа= я* cos at + bk sin cot;

ц>т= amcos (at + bmsin (at.

Подставив (2.30) и (2.32) в (2.316), приравняв выражения с одина­ ковыми тригонометрическими функциями и сократив на скалярные множители cos (at и sin сat, получим систему 2т неоднородных алге­ браических уравнений относительно аъ а2, . . . , ат, blt Ь2, . . . , Ьт

(Е—АМа 2) а + A R ( a b = А Р ' ,

(2.33)

A R ( a a + ( £ — АМа 2) b = А Р ",

которая в развернутом виде может быть представлена следующим об­ разом:

(1 - б пМ 12) а 1- б 12М22а2- • • • - М * « ® 4 , +

тт

+

2 Д Я , Л

= fc2 V V -

 

 

 

 

 

62iAfiC02ai + ( l Ь22М2а2)а2— . . . — 82тМта 2ат +

 

 

А = 1

 

 

 

А=1

 

 

 

 

 

 

—6mxMxCo2a i— бт2М2(о2а2— . . . + ( 1 —6ттМт со2)а т +

 

 

т

 

 

 

т

 

 

 

 

 

 

+

2 bmkRk(abk =

2 бгаЛ

 

;

 

 

 

 

 

fe=l

 

 

 

k—\

 

 

 

 

 

 

т

 

(бцМх(о2— 1) Ь1 + 812М 2(а% +

 

2

\

 

(2.33а)

 

 

 

 

 

 

т

 

 

 

 

 

 

 

 

 

 

 

 

 

+ ' " + 8,a

 

» \ = - кs= 1 v ::

 

 

т

82kRkaak+ б^МхСО2^ + (622М2со2— 1)Ь2 +

 

2

 

a =

i

 

 

 

 

 

 

 

 

 

 

+

• • • + б2A

®

 

!&* = — 2

62*р а ;

 

 

 

 

 

 

 

 

А—О

 

 

2 b m kR k® a k +

b m iM 1(a2b 1 +

 

8 m2M 2(a2b 2 +

 

 

 

 

 

 

аз2— 1\ь

 

 

 

т

 

 

+

. . . +

м

т

<п

=

2 s

bp: .

 

1

1

т т

 

;

 

 

т А А

 

89

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