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

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

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

ной на рис. 128, б,— два (абсциссы ах и а2), а для прямой, пока­ занной на рис. 128, а,— три корня (абсциссы аь а2 и а3). Изме­ нение значений этих корней при постепенном увеличении часто­ ты о показано на рис. 129, а. Конечно, полученные кривые i, II и III представляют собой только приближенное решение задачи, но дают качественно верное представление об изменении ампли­ туды вынужденных колебаний при различных частотах <о.

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

при достаточно высоких частотах ш решение становится неоднозначным; одному и тому же значению частоты со могут соответствовать три значения амплитуды;

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

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

Для построения этой кривой может быть использовано графи­ ческое решение, 'подобное показанному на рис. 128; в данном случае нужно положить PQ= 0 и искать точки пересечения кри­ вой F(a) и прямой z = тр2а, проходящей через начало коорди­ нат (рис. 130). Каждому значению частоты соответствует одно значение амплитуды. Точка р0 (рис. 129, б) соответствует случаю исчезающе малых амплитуд колебаний. Кривая IV амплитуд свободных колебаний является скелетной линией, разделяющей две ветви амплитуд вынужденных колебаний.

Кривые, данные на рис. 129, типичны для систем с жесткой характеристикой (т. е. систем с постепенным увеличением жест­ кости) . Если система имеет мягкую характеристику (когда жест­ кость уменьшается при росте х), то скелетная кривая оказы-

242

вается искривленной влево; соответственно искажаются и кри­ вые амплитуд вынужденных колебаний.

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

Сначала, при постепенном увеличении частоты to, амплитуды также увеличиваются, следуя ветви /. Если, начиная с какой-ли­ бо частоты мь частота уменьшается, то происходит «срыв» ам­ плитуды на ветвь 111 (точки k и

k') и затем медленное возраста­ ние амплитуд до точки k". При последующем уменьшении часто­ ты в точке k" амплитуда резко увеличивается (точка к?" на вет­ ви /) и затем постепенно умень­ шается по кривой /.

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

III. Таким образом, малые колебания частоты около любого зна­ чения toi обязательно приведут изображающую точку в положе­ ние k'\ в этом смысле ветвь III отвечает более устойчивым реше­ ниям, чем ветвь I.

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

Метод Бубнова Галеркина. Согласно этому методу необхо­ димо потребовать обращения в нуль интеграла, подобного (96а):

2~

0

j [mx0-{-F(x0) — P0 sin <at]x0dt—О,

о

подставляя сюда подходящую форму решения xQ{t). Например, приняв за такое решение выражение (379), получим

Зг.

0

J [— таш2 sin ш/ + F (a sin utf) — Р0sin atf] sin wtdt = 0.

о

Выполняя интегрирование, получим нелинейное алгебраи­ ческое уравнение для амплитуды колебаний а; в развернутой форме оно может быть записано после конкретизации вида ха­ рактеристики F(x). Общий характер решения соответствует уста­ новленному выше (см. рис. 129).

243

Метод прямой линеаризации. Согласно этому

методу нели­

нейное уравнение (378) заменяется линейным

 

х + pH = — sin a>t,

(382)

tn

 

в котором величина р2 определяется так, как было пояснено вы­

ше (см. стр. 81).

 

 

 

 

(382)

Стационарная часть решения линейного уравнения

имеет вид

 

 

 

 

 

 

д; =

 

Р0sin art

 

 

-----_______

 

 

 

ш{рг — О)2) ’

 

так что для амплитуды можно записать

 

 

а =

 

Л>

 

(383)

 

m (р г

(О2)

Если, например,

 

 

 

 

 

 

 

 

m

=

р ^

+ « з,

 

 

 

 

 

 

то, как было найдено выше

 

 

 

 

 

2 _

n2 J_ _2_ aaS

 

 

р2=

р1+ —

 

 

и соотношение (383) принимает вид

 

 

а =

_____ Ро

 

 

 

(n ip y

+

maa2 — mco2j

 

т. е.

 

 

 

 

 

~-maa3+

m (p2 — o>2) a = PQ;

(384)

далее из этого нелинейного уравнения нужно определить

неиз­

вестную амплитуду а.

 

 

 

 

 

Метод Дуффинга. В основе этого метода лежит прием исклю­ чения вековых членов, указанный М. В. Остроградским для за­ дачи о свободных колебаниях нелинейной системы. Следуя Дуф-

фингу, ограничимся рассмотрением

кубической

характеристики

tn

= р\х + ах3,

 

w

 

 

когда уравнение вынужденных колебаний принимает вид

х + рЬс +

осх3 =

— sin <оt.

(385)

 

 

tn

 

Решая далее это уравнение, несколько изменим ход рассуж­ дений, использованный Дуффингом; это изменение сделано в со­ ответствии с критическими замечаниями А. Н. Крылова по пово-

244

Постоянные С\ и С2 должны быть определены на данных вы­ ше начальных условиях; используя их, находим

C i = CL +

32w2 9

с2=о

 

аа°

 

и

 

 

х = a sin <at -)—

(sin ut — sin3co/).

1 32ю2 v

'

Таким образом, гармоническая сила Ро sin Ы вызывает в не­ линейной системе не только колебания основного тона с частотой ш, но и колебания с более высокой частотой.

Для построения следующего приближения нужно подставить выражение (384) в правую часть уравнения (387). Далее, для исключения возможности появления векового члена необходимо вновь положить коэффициент при sin Ы равным нулю; это даст уточненное уравнение типа (388), связывающее частоту о и ам­ плитуду а и т. д.

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

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

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

силы, так что если сила меняется

по закону Р — Р0sin со/, то

ко­

лебания в первом приближении

описываются

уравнением х =

= asin(©£— а). Можно описать

колебания

уравнением

х =

=a sin ast, принимая закон изменения силы в виде Р = A) sin (со^+

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

mx + kx-{- F (х) — Р0sin (atf + a),

(389)

а приближенное решение в виде

 

х == asinatf.

(390)

Потребуем, как и выше, чтобы уравнение (389) удовлетворя­ лось в моменты наибольших отклонений (a>t = nf2), когда sin©/ = 1. При этом

х а; х = 0; х = — аш2;

sin + а) = cos а

246

и дифференциальное уравнение (389) заменяется соотношением

 

 

— maco2 + F (а) = Р0cos ос.

 

(391)

Кроме того, потребуем, чтобы решение (390)

удовлетворяло

уравнению

(389) в моменты перехода через положение равнове­

сия со/

= 0, когда sin (со/ + а) =

sin а. Подставляя

в

уравнение

(389)

х =

0; F(x) —0; х —асо;

х = 0, sin (cot +

а)

= sin а, по­

лучим

 

£асо =

Р0 sinoc.

 

(392)

 

 

 

Таким образом, для определения неизвестных а и а служат уравнения (391) и (392). Возводя эти уравнения в квадрат и складывая, исключим фазовый угол а и получим уравнение для искомой амплитуды

 

 

 

 

 

F (а) = } / Р 2 — (£дш)2 + тао)2.

(393)

Отсюда,

при отсутствии сил вязкого сопротивления,

когда

k =

0, следует полученное ранее уравнение (380). При

малых

амплитудах

а

второе

 

слагаемое

под

радика­

 

лом мало

сравнитель­

 

но

с

первым,

 

так

что

 

решение,

полученное

 

выше

(см.

рис. 129),

 

приближенно

верно в

 

рассматриваемом слу­

 

чае. Значительная раз­

 

ница

появится

при

 

больших

амплитудах,

 

когда

второе

слагае­

 

мое под знаком

корня

 

станет близким

к пер­

 

вому

слагаемому.

Та­

 

кое

сближение

двух

 

слагаемых

 

означает,

 

что эффективное значе­

 

ние силы Р0 мало, т. е.

 

решение приближается

 

к решению

задачи

о

 

свободных

колебаниях.

 

 

Наконец, при ka®=

 

— PQ радикал

обраща­

 

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

а = -^ - k(D

является максимально возможным.

247

При этом уравнение (393) принимает вид

 

F (а) = та<л2г

(394)

т. е. определяет скелетную кривую; на этой кривой и следует от* метить наибольшую амплитуду, определяемую уравнением (394).

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

21. ЛИНЕЙНЫЕ СИСТЕМЫ С НЕСКОЛЬКИМИ СТЕПЕНЯМИ СВОБОДЫ

Принципы решения

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

Рис. 133

дельности. При таком под­

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

возмущения Pi sin со^ (или Pi cos сot), здесь Pi — амплитуда воз­ мущающей силы, действующей по i-му направлению; ю — часто­ та возмущения, общая для всех сил, приложенных к различным точкам системы.

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

Особенности каждого из этих способов проследим на простей­ шем примере двухмассовой системы (рис. 133).

Непосредственное решение. Предполагая, что внешняя нагруз­ ка разложена в тригонометрический ряд, исследуем движение системы, вызванное одной гармоникой возмущения. Силы, дейст-

248

вующие на каждую массу, обозначим через Рхsin со/ и P2sinco/. Случай, когда обе силы имеют одинаковую частоту, но разные фазы {Pi sin (со/ + ai), P2sin(co/ -f- a2)], рассмотрим ниже.

Записывая уравнения движения каждой из масс, получим по­ добно выражениям (109) и (ПО)

Рх sin ш/ — сгХх + с2 (х2— лгх) — тхХх,

(395)

Р2sin со/ — с2 (х2Хх) = т2х2,

эти уравнения принято записывать таким образом, чтобы в пра­ вых частях были только возмущающие силы:

ШхХх+ схХх— с2 (х2хг)= Рх sin со/;

(396)

т2х2+ с2 (х2— xa) = Р2 sin со/.

 

Решение этой системы, как и решение одного уравнения, со­ ставляется из двух частей: решения соответствующей однородной системы, имеющего вид выражения (120), и частного решения неоднородной системы (396).

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

Приняв частное решение в виде

 

 

 

 

 

Хх = ахsin ш/;

х2 = а2sin «/

 

(397)

и подставив его в выражение (396), получим

два уравнения G

двумя неизвестными амплитудами колебаний ахи а2:

 

Шхерах + сгах— с2 (а2 — ах)= Рх, \

(398)

 

т2о)2а2+ с2(а2 — ах) = Р2.

\

 

 

Решение (397) означает, что колебания происходят с той же

частотой, с которой происходит изменение самих сил.

 

Решая систему уравнений (398), находим

 

 

 

ах

______ Рх (с3 — т гю2) 4- Р2с2_______.

 

 

1-|—С2— //11СО) (с2— /712®

^2

 

 

 

 

(399)

 

 

 

 

 

а2

Ре (^1 “f” ^ 2 ----/ « ! © - ) -|- Р jCj

(Cj + с 2 — mjto2)(c2 — т 2о)2) —

 

Знаменатели полученных выражений совпадают с левой ча­ стью частотного уравнения (117а), если заменить в нем букву р буквой to. Следовательно, если частота возмущения о совпадает с любой из двух собственных частот рхили /?2,то знаменатели фор­

249

мул (399), согласно уравнению (117а), обратятся в нуль, а ам­

плитуды а\ и а2 станут бесконечно большими

(резонанс).

откло­

При со = 0 формулы (399) определяют

статические

нения обеих масс, вызванные силами Pi и Р2:

 

Рх + Р» ,

 

 

ch =

 

 

 

 

С1

 

Р2

 

_ ^1+^2

г

U.о

 

 

 

 

 

 

С1

 

 

с 2

 

При со—»- оо решения

(399) стремятся

к ну­

лю.

Зависимость

ам­

плитуды

й] от частоты

показана

на рис. 134;

график

построен

для

случая Pj =

1,

Р2 = О,

с\ = с2=

1,

тл = т2

= 1. В этом случае чис­

ло

резонансов

 

равно

двум, что соответствует числу степеней свободы системы и числу ее собственных частот.

При со = 1 амплитуда ct\ = 0; этот случай антирезонанса рас­ смотрен на стр. 263—264.

При помощи выражений (399) можно найти форму вынуж­ денных колебаний, определяемую отношением а2: а\. В общем случае эта форма не совпадает ни с одной из собственных форм; лишь при резонансах форма вынужденных колебаний совпадает с формой свободных колебаний.

Рассмотрим действие двух сил одинаковой частоты, но сдви­ нутых по фазе Pi sin (со/ + ai) и P2 sin (со/ + аг). Эти силы мож­ но представить в виде Pi cos cci sin со/ + Pi sin ai cos со/; P2 cos a2X X sin со/ + P2 sin a2cos со/ и затем решать две задачи:

1)действуют только «синусные» составляющие Pi cos ai sin со/

иР2 cos a2 sin со/;

2)действуют только «косинусные» составляющие Р i sin ai cos со/

ИР2 Sin CC2 cos со/.

Впервой задаче получатся уравнения

 

— m1« 2a1 -f- схаг — е2 (а2 — 0i) = Pi cos ах\

 

— m2о)2 a2 + с2 (a2 — ati) =

P2 cos a2,

а во

второй (когда частное решение

имеет вид Х\ = а\cos со/,

х2 =

а2cos со/)

 

 

m ^ai + сгах— с2 (а2 — ах) = Рг sin ах;

 

— m2oj2a2 -f с2 (а2 — ах) =

Р2 sin а2.

Каждая из этих задач может быть решена до конца; затем ре­ зультаты решения должны быть соответственно суммированы.

250

Если возмущающие силы имеют полигармоническую струк-

туру

Р%= Р11 sin со^ 4" Рla sin о)2/ -f- ... 4~ Р\пsin font‘,

Р 2 = Р21 ИП ®11 + Р22 Sin <°2^ 4" • • • + Pin ®nt>

то резонанс становится возможным в следующих случаях:

% =

л ;

«>2 =

Pi;

10з =

Р ъ

<*>„ =

Pi;

tOj =

р2;

ш2 =

р2;

ш3 =

р2;

=

р2,

т. е. при совпадении любой из /г частот возмущающей силы с любой из двух собственных частот системы.

Разложение решения по собственным формам колебаний.

Вернемся к простейшей системе, испытывающей действие возму­ щающих сил Pi sin соt и Р%sin азt

Для вывода основных зависимостей рассматриваемого спо­ соба предварительно образуем вспомогательные соотношения ис­

ходя из уравнений

(109)

и (ПО). Эти уравнения удовлетворяют­

ся как решениями

 

«и Sin (Pit 4 - «l); 1

 

 

 

 

Хх =

 

(400)

 

 

х 2 =

a2i sin ( p xt

4- #1 ), f

 

так и решениями

 

 

 

 

Xi — 0^12 sin ( p 2^

 

 

 

 

 

^ 2) 5 1

 

(401)

 

 

 

# 2 2 sin (Р г* 4- a 2) * )

 

 

 

* a =

 

 

Подставляем

в уравнения

(109)

и (110)

сначала решения

(400):

 

 

 

 

 

 

 

— «i«u 4 -c2(a2i

«ц) =

— mipian;

1

(402)

 

 

с 2 — «ц) = —m2pia21 ,

|

 

 

а затем решения (401):

 

 

 

 

 

 

4- с%(а22 — «i2) =

— #ziP2 #i2 >

(403)

с 2 (а22 — а12) = — т 2р1а22.

 

 

 

Эти соотношения необходимы для дальнейшего.

В дифференциальных уравнениях (396) неизвестными явля­ ются функции Xi и #2. Основная идея рассматриваемого способа состоит в замене этих функций двумя новыми функциями ti(t), /j}(£) согласно равенствам

Хх —

axxfx 4 "

# 1 2 / 2 ;

1

(404)

*2 =

# 2 l/l +

 

 

a 2 z h >

f

 

здесь «п и « 1 2

произвольные числа (можно, например, принять

#п == й12 = 1), с

которыми « 2 1 и « 2 2 связаны известными соотно-

251