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

книги / Устойчивость и колебания упругих систем. Современные концепции, парадоксы и ошибки

.pdf
Скачиваний:
8
Добавлен:
12.11.2023
Размер:
13.57 Mб
Скачать

из

2п начальных условий — по два условия для каждой

из

обобщенных координат.

Лагранж кратко коснулся того особенного случая, ког* да среди корней частотного уравнения (25.8) имеются рав­ ные корни. В подобных случаях, пишет Лагранж, в реше­

нии возникают члены, содержащие время / вне знака сину­ са (вековые или секулярные члены); например при равенстве

двух корней появится решение вида

Як - (Ak + sin (pi + ф ). (25.9)

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

Однако свободные колебания консервативной системы не могут неограниченно возрастать — это противоречило бы очевидным энергетическим соображениям; при таком про­ цессе энергия системы увеличивалась бы до бесконечности и система представляла бы собой своеобразный вечный дви­ гатель.

В чем же дело? Может быть, в задачах о свободных колебаниях консервативных механических систем случай равных корней невозможен?

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

Рис. 26.1. Система с двумя сте­ пенями свободы. При равных жесткостях в направлениях к н у собственные частоты равны одна другой

У

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

J91

Эти два примера непосредственно свидетельствуют о том. что случай равных корней физически легко осуществим, причем появление вековых членов в решении противоречит здравому смыслу и закону сохранения энергии. Однако Лагранж уклонился от обсуждения этого парадокса: «...так как изложение этих случаев не представляет интереса для рассматриваемого нами вопроса, то мы на нем не будем ос­ танавливаться».

Утверждение Лагранжа о появлении вековых членов было опубликовано в 1788. г., когда вышло в свет первое издание его книги; оно было оставлено Лагранжей и во втором издании, которое вышло из печати в 1812 г.— за год до смерти Лагранжа. В 1853 г. «Аналитическая меха­ ника» вышла третьим изданием, которое редактировал зна­ менитый математик Жозеф Бертран; и здесь утверждение Лагранжа было сохранено.

Лишь в 1858— 1859 гг. К. Вейерштрасс *) и О. И. Со­ мов **) установили, что утверждение Лагранжа ошибочно, и д а ж е в с л у ч а е р а в н ы х к о р н е й частотного уравнения решения сохраняют чисто тригонометрическую форму, а вековые члены вида (25.9) появиться не могут.

Дело заключается в следующем. В случае кратных кор­ ней система вида (25.1) имеет обычно решение с множите­ лями /, Р и т. д. Но мы имеем дело с о с о б ы м случаем системы линейных дифференциальных уравнений, которая обладает определенными свойствами симметрии (25.2). В этом особом случае она и при наличии равных корней для рг имеет только решения вида (25.3) ***).

Остановимся подробнее иа механической системе, в ко­

торой случай равных корней

не только возможен,

но и

ж е л а т е л е н . Рассмотрим

задачу о свободных

верти

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

жен сам автомобиль, а на

рис. 25.2,6 — его упрощенная

•)

Карл

Вейерштрасс

(1815— 1897) — немецкий

математик,

с 1856 г.— профессор

Берлинского университета,

иностранный

член-

корреспондент

(1864 г.) и иностранный почетный член (1895 г.)

Петер­

бургской

Академии

наук.

 

(1815—1876) — с

1841 [.

профессор

**) Осип

Иванович Сомов

Петербургского

университета,

с

1857 г— член-корреспондент

Петер­

бургской

Академии наук. Автор

ряда крупных работ по математике

ианалитической механике.

***) Симметрия матрицы коэффициентов не является необходимым

условием для отсутствия вековых

решений в случае

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

См. об этом

в книге

А. И. Лурье «Операционное

исчисление! (М.:

Гостехиэдат,

1950, §9,

10, с.

106— 113).

 

19?

динамическая схема (в более полном анализе учитываются также массы колес и связанных с ними частей). Приведен­ ные коэффициенты жесткости задней и передней подвески автомобиля (учитывающие жесткости рессор и пневматиков) обозначены через сх и са. Центр тяжести расположен на расстояниях а и р соответственно от задней и передней опо­ ры. Масса кузова обозначена через т, а момент инерции

кузова

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

 

перечной

оси.

проходя­

 

щей через

центр тяжес­

 

ти,

через

дар* • (р — ра­

 

диус инерции).

 

 

что

 

 

Предполагается,

 

деформации самого кузо­

 

ва

пренебрежимо

малы

 

по

сравнению

с осадка­

 

ми

опор;

поэтому

 

в ди­

 

намической

схеме гори­

 

зонтальный стержень бу­

 

дем считать

совершенно

 

жестким.

Кроме

 

того,

 

положим,

что

горизон­

 

тальные

колебания

си­

 

стемы невозможны.

об­

Рис. 25.2. Автомобиль и его упрощен­

 

Такая

система

 

ная динамическая схема для анализа

ладает

двумя

степеня­

колебаний в вертикальной плоскости

ми свободы

и

за

обоб­

 

щенные координаты <?1 и qt примем вертикальное пере­ мещение центра тяжести (положительное направление вниз) и угол поворота стержня вокруг центра тяжести (по­ ложительное направление — против хода часовой стрелки). При этом осадки упругих опор равны соответственно Qt+aqt (задняя опора) и <?, ~bqt (передняя опора).

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

 

т _ Я,Ч1 .22*2?

(25.10)

 

1

 

*"

2

 

 

 

П _cl

( f i - h o f i t 2

1 Ъ

(<?1—А?*)*

(25.11)

1

2

 

 

1

2

Соответственно уравнения

Лагранжа

 

d

( дТ

\ _

 

ОТ

<И1

(25.12)

di

\dcjj )

&!/

дсц

 

193

принимают вид

m qi + (С, -1 с,) </, +

(c,fl^c,/>) q , = 0 ,

{3

mp*qt + (cta — сф) ql + (с,а2 f cj^) q2 = 0.

' ' 1

После подстановки

частного

решения

(25.3) приходим

к однородной системе уравнений

 

 

 

—тр'Ах+ (с, +

с,) Аг + (с,й —саЬ) Аг= О,

(25.14)

—/пр>М а + {сха—cjb) Аг + (сга*+ сай*) Ла = 0

 

и получаем определитель частот

 

 

 

)с,+ с,—тр*

с/х—сф

I

п

(25.15)

\с\а—сф

cta2 -(-са64—т*рар*1

'

из которого следует частотное уравнение

 

 

 

 

 

 

<*■«»

Как правило, это уравнение имеет р а з л и ч н ы е

корни.

Каждому из них соответствует свое определенное соотно­ шение между амплитудами, определяющее собственную форму колебаний. Это соотношение можно найти, напри­ мер, из первого уравнения системы (25.14):

Аг

тр*с-t —с%

(25.17)

А,

cto - ctb

 

То же получится и из второго уравнения системы (25.14). Остановимся на частном случае, когда выполняется ра­

венство

 

р* *= ай,

(25.18)

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

Если условие (25.18) выполняется, то частотное урав­ нение приобретает вид

р- -

«■+!>)<>■ +

- о

<*•!#>

и имеет следующие неравные корни:

 

 

, _

Ci (о+о)

_ fi (о+ 6)

(25.20)

Р*

mb

*

то

 

194

Для определения собственных форм подставляем эти корни поочередно в соотношение (25.17). Для первой формы полу­ чим, полагая р*=р\:

__

(25.21)

л» ь

 

Здесь вторые индексы обозначают номер рассматриваемой собственной формы. После подстановки p*=pl в соотноше­ ние (25.17) найдем вторую собственную форму

 

____ J_

(25.22)

Alt

а

 

Эти формы изображены на рис. 25.3. Их особенностью является неподвижность одной подвески автомобиля прн колебаниях другой подвески. Практически это означает.

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

Рис

25.3. Собственные фор­

Рис. 25.4. Динамическая

схема

мы

колебаний

автомобиля

автомобиля

с независимыми

 

при р*=а6

колебаниями

обеих

подвесок

ством рассматриваемой

схемы.

Конечно, свободные коле­

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

р а з ­

л и ч н ы м и

частотами.

 

 

 

 

 

 

В другом

частном случае,

когда

 

 

 

 

 

 

cto = гф,

 

'(25.23)

уравнения (25.14) становятся

независимыми:

 

 

 

 

/4, (—трг+ с, + с,) =

0.

 

(25.24)

 

 

Л ,(—

+<;,«*+ с,6*)==0.

(25.25)

Независимость уравнений означает возможность чисто вертикальных колебаний без поворотов (подпрыгивание)

195

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

 

^ # 0 ,

 

4 , =

0

(25.26)

и из уравнения

(25.24)

 

 

 

 

 

тр*-'Гсг

С> =

0

(25.27)

получаем квадрат

первой частоты

 

 

 

Pi

«1+ с»

 

(25.28)

 

т

 

 

Для случая галопирования автомобиля, когда

 

 

Л, = 0,

АгФ 0,

(25.29)

из уравнения (25.25) находим квадрат второй частоты

 

с,а*-|

с**2

 

(25.30)

 

 

трг

 

 

 

 

 

Таким образом, частоты подпрыгивания и галопирования Pi и piy вообще говоря, различны.

 

 

 

 

Возможно ли равенство обеих

 

 

 

 

частот? Схема автомобиля, обла-

_

__^ __________ -

дающего таким свойством, приз-

~

,-^fj

 

-5 ё-

нается

особенно

рациональной

 

 

а)

ш

с точки

зрения

комфортабель­

 

 

 

ности езды. Ясно, что в этом

 

 

 

 

случае сильные колебания авто­

 

 

 

 

мобиля

возникают

лишь

при

 

 

 

 

некоторой

е д и н с т в е н н о й

 

 

 

 

частоте

толчков,

тогда как в

 

 

 

 

случае Р\Фр* сильные колебания

Рис. 25.5. Собственные фор­

смогут

возникнуть

при

двух

мы

колебаний

 

автомобиля

значениях этой частоты *).

вид­

 

при cta=ctb

Из

уравнения

 

(25.16)

но, что его корни оказываются

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

+

1* 4CiC« (g-j-Ь)г = 0.

(25.31)

*) Другим примерим систем, для которых равенство собственных частот считается достоинством, могут служить центрифуги с жестким ротором; в некоторых патентах таких центрифуг равенство частот специально оговаривается.

Ш

Легко проверить, что оно выполняется, если о д н о в р е ­ м е н н о удовлетворены оба условия (25.18) и (25.23). Более того, можно доказать, что эти условия не только достаточны, но и необходимы для выполнения равенства (25.31).

Если динамическая схема автомобиля удовлетворяет

условиям (25.18) и (25.23), то

е д и н с т в е н н а я

собст­

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

 

Р - У

.

(25.32)

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

4i= Ai sin (pt -гф|), <?г= Ла51п(/7/-|-фг). (25.33)

Постоянные Ai, А,, ф! и ф, определяются начальными

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

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

Вот пример, который приводит в своей книге А. И. Лу­ рье *) (см. ссылку на с. 192). Дифференциальные уравне­ ния движения системы с двумя степенями свободы имеют вид

Qi + Я» + 9i + 9* =

(25.34)

9*4* 3^* + 9i + 3$, = 0.

 

Разыскивая частное решение в виде

9[ = aie,J,

== агеи .

(25.35)

*) Анатолия Исакович Лурье (1901—1981). член-корреспоидент АН СССР, в течение пятидесяти лет профессор Ленинградского поли­ технического института, автор многих трудов по теории упругости, устойчивости систем автоматического управления, аналитической механике.

197

мы придем к следующему характеристическому уравнению:

 

 

 

(fc* +

+ 2) (Я,Н- 1)* =

0.

(25.36)

Среди его

корней

 

 

 

А* = — 1,

К =

>», = — 1 - И ,

^ =

— 1— » (25.37)

один является

кратным, и вследствие этого решение с о ­

д е р ж и т

вековые члены. Например, при начальных ус­

ловиях

 

 

 

 

 

 

 

<7i (0) = Яг(0) = <?,(0) = 0, qt(0 )^ 0

можно найти,

что решение имеет вид

 

 

 

?i =

?i (0)(te- t + e -<—е~х cost),

 

<7, =

qx(0) (te_< —е~* sin t)

 

 

и свидетельствует о

з а т у х а ю щ е м

д в и ж е н и и

(так как ~*->0 при t-*-оо), хотя и содержит вековые члены.

Перевод «Аналитической механики» Лагранжа на русский язык вышел в 1950 г. (М.: Гостехиздат, тт. I, 2). Об ошибке Лагранжа см.

в «Лекциях по

колебаниям» Л. И. Мандельштама (М.: Наука, 1972)

и в упомянутой

на с. 207 книге А. Н. Крылова.

§ 26. Пропущенные частоты

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

Л = «1 V ~ .

(26.1)

в котором EJ — жесткость при изгибе, т — интенсивность распределенной массы.

По значениям корней характеристических уравнений можно заключить, что собственные частоты убывают при переходах от схемы 1 к схеме 3, а затем вновь начинают увеличиваться (можно заметить даже строгую симметрию результатов относительно «средней» схемы 3: они совпадают для схем 1 и 5, а также для схем 2 и 4). Это может вызвать

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

198

л«

сх емы

I

2

3

4

Схема балки

Я

ш

i

ш

1---------------------------

А

1------------------------------

 

4 4 /,

Перпые три корня характеристического

уравнения

Характеристическое

уравнение

№1),

«*/),

(«0,

c h a /c o s a /= l

4,730

7 ,8 5 3

10,996

tb a l = tg al

3,927

7,068

10,210

ch al cos al = — 1

1,875

4,694

7,854

th a / = t g al

3,927

7 .0 6 8

10,210

5

ch al cos a l« 1

4,730

7,853

10,996

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

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

Прежде всего запишем дифференциальное уравнение свободных колебаний

E Jl F + m -W = 0’

(26.2)

в котором v(2, /) — вертикальное отклонение произвольной точки оси балки от положения равновесия, г — горизон­ тальная координата точки, / — время. Будем искать част­ ное решение в виде произведения двух функций

v= V(z)T(t).

(26.3)

Тогда уравнение (26.2) примет вид

или, иначе,

EJV'vT-\ mV

О,

 

 

 

 

 

...

 

 

£ J V l v

Т

 

 

 

mV

Т ~ р '

 

где р — пока

неизвестная

постоянная («постоянная разде­

ления). Таким образом,

вместо

 

уравнения

(26.2) получа­

ются два следующих уравнения:

 

 

 

 

t + p*T = 0,

(26.4)

 

V'v a*V= 0

(а* =-•

(26.5)

Из решения

уравнения

(26.4)

 

 

 

 

Г =

/4 sin (pt + y)

(26.6)

непосредственно видно, что постоянная р представляет со­ бой собственную частоту. Для ее определения нужно об­ ратиться к решению уравнения (26.5):

V= С, sin ал + Сгcos az -f Ся$h «г 4- С, ch аг, (26.7)

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

200