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

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

..pdf
Скачиваний:
4
Добавлен:
13.11.2023
Размер:
23.13 Mб
Скачать

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

из них эпюра Мкбыла ограничена прямой линией. Напомним также о спра-

ведли'вости равенства = = б/it, составляющего содер­ жание теоремы Максвелла.

Основные уравнения за­ дачи и частотное уравнение.

Рассмотрим свободные ко­ лебания балки, несущей массы mb m2,...,mn. Разви­ ваемые ими силы инерции

—т\Уи —т2у2, ...» —тпуп яв­ ляются единственной нагруз­ кой на упругий «скелет» си­ стемы в процессе колебаний. Поэтому соответственно формуле (133) можно запи­ сать для любого перемеще­ ния

yrflfff.. \ТТШШтгтгг^

Рк= 1

П

 

Уг = “ 2 ЩУ^Чк

 

k = l

 

<* = 1,2, — п),

 

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

 

Ух = — Щ У х Ь и —

т.гу2б12 — . . .тпуп81п;

У* = — тхУхб21

— ... — т„упб2„;

 

(134)

Уп = — тгУгдпх — пцу»ЬП1 — ...тпупЬпп.

Входящие сюда коэффициенты т&имеют обобщенный смысл; так как в тех задачах, где одно из перемещений ун есть угол по­ ворота, под frih следует понимать не массу, а момент инерции k-ro

груза относительно оси, проходящей

через

его центр тяжести

перпендикулярно плоскости колебаний.

в простейшем случае

од­

Основная система уравнений (134)

ной степени свободы приводит к одному уравнению с одной

не­

известной функцией

 

 

 

 

У\=

п^ухви;

 

тух-f сух = 0,

 

оно эквивалентно стандартному уравнению

так

как с = 1 : 6ц.

 

 

 

 

Частное решение системы (134) может быть

принято в виде

yi = ats\n(pt + a).

(135)

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

&1 (Шхбцр2 — 1) 4- а-ъМгЬхър1“Ь ••• аптпд1пр2= 0;

a im i6 21p2 + а2 ( т 2622р2 — ! ) + • ♦ • ^ т , гб 2ггр2 = 0; ( i 3g)

ахШхЬгар 4- я2ш2бп2р 4- • • • 4~ ^ni^n^nnP^ 1)

Эта система удовлетворяется тривиальным решением

Q x == Q-ч = . . . CLn = 0 ,

которое соответствует отсутствию колебаний, но не является един­ ственным. Амплитуды «г одновременно не обращаются в нуль, если определитель, составленный из коэффициентов системы уравнений (136), равен нулю:

mifiiip2 — 1

т 2б12р2

.

. .

. mn8inp2

 

/Я1621Р2

т 2б22р2 -г- 1

.

. .

. тпЪпР2

= 0. (137)

mi6„ip2

т 26л2р2

. . . .

Мп&ппР2 — 1

 

Развернув определитель (137), получим частотное уравнение п-й степени для квадрата частоты р2; напишем это уравнение в виде 1

1 _

ЬхР*+ Ьгр1 -

63ре +

. . . + ( _ l)nbyn = 0.

(138)

Число корней уравнения

(138)

равно п\ обозначим эти корни

через р f , р|

.... р2 , расположив

их в порядке возрастания. Все

эти корни вещественны и положительны, поэтому для частоты р

получится также

п ответов

(отрицательные решения несущест­

венны) :

__

 

___

 

___

 

Pi = + V Pi

; р2 =

+ V pt ,

рп=

+ V pl .

 

Общее решение системы уравнений

(134)

составится

в виде

суммы частных решений типа (135):

 

 

 

 

П

 

 

 

 

 

У-1=

2

aiksin {pkt + *k) (« =

1 , 2 , . . . , n),

(139)

 

&e=l

 

 

 

 

 

1 При такой расстановке знаков все коэффициенты Ь\ положительны.

104

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

Уг =

0ц sin (pit +

ai) 4* ап sin (p2t + a2) ... -f 0irt sin (pnt -f- a„);

y2=

a21 sin (pit +

c^) + a22 sin (p2t 4- <*2) 4- . •. 4- 02rt sin (pnt 4- <x„);

Уп = amsin (pi* + otj) + an2sin (p„( 4- «„) 4- - 4-0„n sin (p^ 4- «„)•

Таким образом, в каждом направлении i = 1, 2, ..., /г проис­ ходят многотонные колебания, причем число слагаемых равно числу собственных частот, т. е. совпадает с числом степеней сво­ боды системы п.

Сосредоточим внимание на вычислении собственных частот, оставив 'пока в стороне вопрос об определении амплитуд а,-.

В частном случае системы с двумя степенями свободы систе­ ма уравнений (136) принимает вид

0i (0*i6np2— . 1) 4- 02m26i2p2 =

0;

(140)

0i0*i62iP2 4- 02 (W0622P2 —

1) =

0,

 

и вместо определителя (137) имеем

 

 

 

mi6np2 — 1

m26i2p2

=

0.

(141)

 

 

0*i62iPa

m2622p 2

1

 

 

Частотное уравнение (138) становится биквадратным и име­ ет следующие решения:

р 1,2 — [0*1611 4" 0*2622i

(0 *1 6 1 1 /Я2622)2— 4 (611622— 612) 0*10*2]:

• 2/7ii/02 (611622 612).

(142)

Пример 11. Найти собственные частоты колебаний консоли, показанной на рис. 54, а с учетом инерции вращения концевого груза. Эта система имеет две степени свободы и ее движение описывается двумя координатами уу (вер- тикальное перемещение груза) и у2 (угол поворота груза).

Строим эпюры изгибающих моментов от силы Pi = 1 и момента Р2 = 1 (рис. 54, б— д) и вычисляем по формуле Верещагина

S“ = 1 E T ; 6u=6*,=

Ш ; 6г,= "ЕГ-

Далее пользуемся формулой (142),

подставляя в нее Щ — m, т 2 = т р 2

(р — радиус инерции груза):

 

105

Если, как обычно, отношение I : р велико, то последний корень можно при­ ближенно преобразовать так:

/

J L

 

 

-

 

 

9р4

,

3£2_

+ 1

Зр2

^

 

Зр2

] /

 

+

1 Г

9р4

У

 

 

 

1р2 V21*

т

3PL

\

Зр2

1

 

/2

 

 

2/2

^ /

2/2

г 2

т

Зр2

 

Тогда

6EJ

 

 

 

 

 

 

 

 

 

2

[ ' +

зр»

* L + - L + - £ - M

 

Р Х , 2 ~

т Р

± ('

212

 

Зр2 )\ ’

Принимая знак минус, получим квадрат низшей (первой) частоты:

 

Pl =

6EJ /

1

_

Зр2

\

3EJ

f

Зр2

\

 

 

т Р

[

2

 

2/2

/

т Р

V

/2

)

 

 

 

Квадрат высшей (второй) частоты оказывается во

много

раз большим:

 

 

6EJ I 2Р

3 \

3EJ

.

4Р

)

 

 

Р2 =

/Я/3

V Зр» т

2

j

ml*

V т

Зр»

 

Так,

например,

если

р :1 =

0,1,

то

вторая частота

оказывается в 137 раз

больше

первой.

 

 

 

 

 

 

частдты колебаний

балки, несущей

Пример 12. Определить собственные

три сосредоточенных груза (см. рис. 23).

 

 

 

 

 

106

Для определения единичных перемещений строим эпюры моментов от еди­ ничных сил Pi = 1, Рг = 1 и Р2= 1 (рис. 55). По формуле Верещагина нахо­ дим единичные перемещения

бц — 633 = 75fe; 622 — 243й; 621 6J3 = 632 = 6JJ = 117&; 6^3 = 83^ = 51fe.

Р

В этих выражениях к =

Теперь составляем определитель (141):

75mfep2 — 1

117ткрг

51mfep2

= 0.

117mkp2

243/nfep2— 1

117mfcp2

51mftp2

117mfep2

75mfep2 — 1

 

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

77 760 {mkjPf — 12 096 (mjfep2)2 + 393/nftp2 — 1 = 0

и имеет три корня

(на операциях вычисления этих корней не останавливаемся).

Собственные формы колебаний. Если в систему уравнений (136) подставить какой-либо из корней частотного уравнения, то одно из уравнений станет следствием остальных, т. е. незави­ симых уравнений оказывается только п — 1. Эти уравнения связывают между собой п амплитуд ацуa2i, ...» ап* [второй индекс означает номер корня р ?, который подставлен в систему (136)]

и поэтому позволяют выразить все амплитуды через какую-либо одну, например, первую. Отношения ац-.ац, а3{:ац,... опреде­ ляют t-ю собственную форму колебаний. Таким образом, каж­ дой частоте отвечает вполне определенная форма колебаний; число этих форм совпадает с числом степеней свободы системы

и равно п.

установить, что любые две собственные

Как и выше, можно

формы колебаний взаимно ортогональны:

 

mxaualk + тгаиа.2к+ ... + mnaniank= 0.

 

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

собственные формы колебаний консоли

(см.

рис. 54). Так как система обладает двумя степенями свободы, то число

соб­

ственных форм также равно двум..

Подставив данные задачи

m, = m; т г = « р > ; *u =

впервое уравнение системы (140)

а2 (тхбцр2 — 1) -}- О2/Я2612Р* — 0,

107

Принимая здесь
находим первую собственную форму:
Ряс. 56

получим

/ ml3

\

ml2р2

■р2 = 0.

а‘ [ Т Ё 7 р

1) + а °-~2Ё7

 

Она иллюстрирована на верхней схе­ ме рис. 56 и характеризуется в основ­ ном линейным смещением концевого груза.

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

нять р- = р|; отсюда находим

g22 ^

2/

&12

Зр2

эта форма иллюстрирована на нижней схеме рис. 56 и характеризуется глав­ ным образом поворотом концевой массы.

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

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

Представим уравнения (136) в виде

ах = р2(ш^ ц Дх+ m26 i2 a2

+

. .. + тпд1пап);

 

а2— р2

+ m2 6 22a2

+

юп82пап)‘>

(143)

 

 

 

 

ап= Р2 (mi8niai + m26n2a2 + . . . + тпЬппап).

Задавшись произвольной системой амплитуд (fli)o, {0.2) 0, ...» (fln)o (нулевое приближение для формы колебаний), подставим их в правые части системы (143) и вычислим выражения, входя­ щие в скобки. Пусть полученные величины будут (6 1 ) 0, (&г)о, ...» (Ьп)о- Тогда

ai = Р2 Фг)0; а2= р2 (6 2)0; . . . ап= ра (&„)„•

Эти выражения можно рассматривать как новое приближе­ ние для формы колебаний. Без ущерба можно опустить общий множитель р2, так как масштаб формы колебаний не имеет зна­

108

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

(ai)i = Фх)ф (0 2 ) 1 = Ф Х ... (an)i = ф Х

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

в(k — 1) -м приближении будет отличаться от системы амплитуд

вk-м приближении только общим масштабом, т. е.

(a l)A

__{ а ъ ) к ________

 

{ ап)/г

(ai)k -i

1

* * ’

(an)k- i '

Этим завершается определение основной формы колебаний. Затем сразу находим и основную частоту: так как

Фдм = Р*Фдь- ь

то для квадрата частоты получим

2 e (g*)*

(g4 - 1

Колебания автомобиля

Будем рассматривать автомобиль как систему упруго связан­ ных жестких тел 1—5 (рис. 57, а). Здесь тело 1 схематически представляет собой кузов автомобиля, а тела 2—5 — колеса, мас­ сы которых примем сосредоточенными.

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

Ух — вертикальное перемещение центра тяжести ку­ зова;

У 2г У ъу */4» Ув— вертикальные

перемещения

центров тяж ести

колес;

 

поперечной

оси;

Ув— поворот кузова относительно

г/7 — поворот кузова

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

продольной

оси.

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

вать состоящим

из двух взаимно не связанных процессов

(рис. 57, бив):

 

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

мещением кузова

(ух), поворотом кузова вокруг поперечной оси

(Уб) и попарно равными перемещениями обоих передних колес (Уг = У а ) и обоих задних колес (у3 = Уъ) ;

109

поперечных (боковых) колебаний, характеризуемых поворо­ том кузова вокруг продольной оси (,у7 ) и попарно равными пере­ мещениями обоих левых колес 2 = Уз) и обоих правых колес

(#4 = */5).

Соответственно этому продольные колебания описывают­ ся четырьмя, а поперечные колебания — тремя дифференциаль­ ными уравнениями. Рассмотрим продольные колебания, которые имеют основное значение.

Обозначим жесткости шин через с, жесткости передних и зад­ них рессор соответственно через сп и с3, массы кузова и колеса

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

Ап = У\ + ау6 у2 (передняя рессора), Аз = У\Ьу6— г/ 3 (задняя рессора).

Уравнения движения составим в форме Лагранжа. Кинетиче­ ская энергия системы складывается из следующих частей:

т у \

дара0§

—------ 1-------------кинетической энергии кузова;

Z

Z

о

тьи\

 

*

—------- кинетической энергии передних колес;

 

о

•2

 

ткУ3

 

* —---------кинетической энергии задних колес.

 

Суммарная кинетическая энергия

 

^

# 1 + р2 1/б) + 2 т л (yl + уз) J

(144)

ПО

Потенциальная энергия состоит из энергии деформации рессор

2 ■ !£ ? .? -+ 2 -^ - = сп (Угih 4- аУв)2+ с3 (ух Уз Ьуо)2

и энергии сжатия шин

2 С + 2 ^ - = с (у1 + у1).

Суммарная потенциальная энергия

 

 

 

П = сп (У1 — У2 + аУз)2+ са (ух — уз — Ьуо)* + c(yl + у\).

(145)

В данном случае уравнения Лагранжа имеют вид

 

tnyi + 2с„ fa — у2+ ауй) + 2са {уг - у3Ьув) = 0;

 

2 ^ з — 2са (ух —уй+ аув) + 2су3 =

0;

^ 4g

2тку3 — 2с3(у1 —Уз — Ьу6) +

2суй=

0;

 

mp2i/8 -1- 2сп(уг— уъ + аув) а — 2с3 {ух — у9 — Ьул) 6 = 0..

 

Частное решение этой системы

 

 

 

у{ = atsin (pt + а)

( / = 1 , 2 , 3 , 6).

(147)

Подстановка выражения (147) в систему уравнений (146) приведет к уравнению той же структуры, что и в задачах, рас­ смотренных выше, в частности, соответственно числу степеней

свободы обнаружатся четыре значения собственной частоты ко­ лебаний.

Считая шины недеформируемыми, рассмотрим упрощенную схему продольных колебаний (рис. 58); такая система обладает двумя степенями свободы, соответствующими координатам ух и Уе. Положим в выражениях (144) и (145) уз — уз 0. Тогда со­ ответственно для кинетической и потенциальной энергии получим

п = сп {Уг + ауву + сэ (ух — Ьув)а.

111