книги / Основы прикладной теории упругих колебаний
..pdfПри пользовании формулой Верещагина необходимо разде лить длину стержня на такие участки, чтобы в пределах каждого
из них эпюра Мкбыла ограничена прямой линией. Напомним также о спра-
ведли'вости равенства = = б/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
получим
/ 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