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

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

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

приложения л-ro и (л+1)-го импульсов, имеем

9 “ ^

(sin /»/

!-sin p ((— Т)-\

 

 

+ sin p (f— 27)

+ sin p (( — nT)\

 

 

apt?o

sin p (t—kT)

(n T < t< ( n + l) T ).

(28.5)

 

 

 

 

t—nT,

Для

начала

этого

интервала времени, т. е. при

по выражению (28.5) находим

 

 

 

Я=

П-I

 

 

 

 

£

sin (л—ft) р 7 .

(28.6)

Вывод выражений (28.5) и (28.6) прост и ясен, а их вид может показаться даже привлекательным своей внешней компактностью. Однако наши симпатии исчезнут, как только мы вдумаемся в процедуру необходимых выкладок и представим себе их объем. Пусть, например, нужно оп­ ределить движение в интервале времени (997, 1007’) и найти отклонения через каждую десятую часть этого интервала, т. е. в моменты (д= 997, ^=99,1 Г, (*=99,27, . . ., = =99,97, (10= 1007. Согласно выражению (28.5) нам придется вычислить одиннадцать сумм:

99

9 .“ -^ £ sin (99 — ft)p7,

99

 

9 i = - i £ s i n ( 9 9 , l - f t ) p 7 ,

(28.7)

У*=•

 

99

=sin (100—ft) р7,

*=ж«

вкаждую из которых входит сто слагаемых.

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

Хотя эти неудобства описанного способа очевидны, но мы пока воздержимся от окончательных оценок: нужно еще

посмотреть, каковы достоинства и недостатки других спо­ собов.

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

q« A cos pt + В sin pt,

где

Д = — ^ { Р ( т ) sinprdx,

О

t

я = —

cos px dx.

 

II

(28.8)

(28.9)

Теперь положим, что в моменты времени О, Т, 2Т, . . .

к системе прикладываются одинаковые импульсы S. Тогда для интервала времени \пТ, (л+1) Т\ по выражениям (28.9) находим

А = — ^ [ s i n рТ + sin2/> 7'+ . . . + sin npT], (28.10)

5 = ^ -[1 Ч w s p T cos2pT -{-. . . + cosnpT], (28.11)

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

дящую в

выражение (28.10),

на 2 sin

и запишем ряд

очевидных

равенств

 

 

2 sin -^ - sin рТ = cos —

c.os^~ ,

 

2sin

sin 2рТ — cos Щ -

cos

,

2 s in

sin npT =»cos

pT — cos

pT.

Сложив все эти равенства, получим вместо выражения, стоящего в скобках формулы (28.10),

cos ~ — c o s — ^— рТ

 

s in — i — рГ -shi у

рТ

- —- - f

=

-~~f

.

^*ш —

 

s*n ~~2~

 

222

Аналогично этому, выражению, стоящему в скобках фор* мулы (28.11), можно придать форму

, Sff -f~ 1

м . . рТ

sin — g—

рТ + sin

2sin рТ

Окончательно формулы (28.10), (28.11) приобретают вид

 

И+1

 

 

S %in:zY -P T-%Xn2 Рт

 

А =

 

~Грт

 

 

 

sin

 

(28.12)

 

2n

1 _

р Т

В -

sin —

2рТ

-1-sln

 

ар

sin РТ

 

 

 

 

Таким образом, для определения движения на любом ин­ тервале нет необходимости в утомительном суммировании, которое требуется в первом способе; нужно лишь предвари­ тельно найти числа А и В по формулам (28.12), а затем воспользоваться весьма простым выражением (28.8).

Тем не менее рассматриваемый способ все же содержит

один существенный

недостаток.

Поскольку решение, данное в выражениях (28.8) и

(28.12), не является

периодическим и м е н я е т с я от

одного периода к другому (в эти выражения входит номер л), то для полного описания движения придется повторять вычисления заново для каждого нового периода. Эго тем более досадно, что, по-видимому, главной частью движения несомненно является п е р и о д и ч е с к и й п р о ц е с с , «навязываемый» упругой системе импульсами, повторяю­ щимися с периодом Т.

Чем объяснить, что полученное решение не обнаружи­ вает такой периодичности? Ответ на этот вопрос звучит несколько парадоксально: это решение «слишком» точное — в нем содержится не только описание стационарного коле­ бательного процесса с периодом Т, но и отражено влияние начальных условий, которое в действительности из-за не­ упругих сопротивлений очень быстро затухает. Кажется, что эта «примесь» к стационарному процессу вообще не может быть исключена, поскольку исходное уравнение записано для строго упругой системы, лишенной неупругих сопротивлений.

223

Можно ли вообще получить чисто периодическое реше­

ние,

если

опираться

и м е н н о н а э т о

у р а в н е ­

н и е ,

т. е.

не вводить

в уравнение члены,

выражающие

демпфирование? Да, можно, как это показывают два сле­ дующих способа решения.

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

суммы

гармоник

с периодами

Т,

Т , -j

Т............ т. е.

в виде

ряда

Фурье

 

 

 

 

 

 

Р(1) = $

+

Ё

а» cos

2nnt

f Ц

К sin

2яя/

(28.13)

 

 

 

П»1

 

~ Т ~

л = |

 

Т

 

Идея решения состоит в том, чтобы сначала найти движения, вызываемые каждой из гармоник, а затем — опираясь на принцип наложения — просуммировать эти «парциальные» решения. Для вычисления коэффициентов On и Ьп следует воспользоваться формулами Фурье

г

 

o „ = |r S / , (0 co s-~ -fiK

(л = 0 . 1,2, . . . ) ,

г

(28.14)

bH= ^ \ p ( t ) s m ~ i d t

<л=1. 2. ...).

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

г

 

 

4<

 

$ Я (0 cos

dt =

Urn

J р (t) cos

dt = S,

о

 

 

о

(28.15)

 

 

 

 

J Р (*) sin

■dt =

lim

J P (t) sin

dt — 0

о

 

At-*0

 

 

и по формулам (28.14)

находим

 

224

Таким образом.

(28.16)

Я*1

и нам придется иметь дело не с мгновенными импульсами, а с бесконечной суммой, входящей в выражение (28.16);

,

2S

2 п я (

общий

член суммы имеет вид - у cos - у —.

Рассмотрим движение, которое вызывается действием только одного такого слагаемого. Для этого нужно решить дифференциальное уравнение движения

dta

2 S

2пЯ?

(28.17)

а Т COS-

 

Частное решение (не связанное с начальными условиями) ищем в виде

9я = Ля с о з - ^ .

(28.18)

Амплитуду Ап можно найти путем подстановки решения (28.18) в уравнение (28.17):

( - i = ? . ) 4 + P м „ - £ .

откуда следует

(28.19)

Таким образом, решение (28.18) имеет вид

25

cos

2лл/

 

 

Т

(28.20)

Суммируя отдельные решения типа (28.20), т. е. складывая движения, вызываемые каждой из гармоник ряда (28.16), находим для одного периода

л __ s

, 25

V

(28.21)

а р Ч

'H tT

Z +

 

Здесь первый член соответствует первому слагаемому SIT в выражении (28.16). Хотя решение (28.21) требует сумми­ рования бесконечного ряда, но этот недостаток не может

заслонить важное достоинство решения

(28.21) — его п е-

р и о д и ч и о с т ь, благодаря которой

достаточно найти

225

движение в течение всего лишь одного периода. Это свойстворешения (28.21) выгодно отличает его от непериодических решений, полученных двумя первыми способами.

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

1.Каким образом удалось исключить «примеси», на­ рушающие периодичность?

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

Ответ на первый вопрос прост. Ведь решение (28.18) уравнения (28.17) представляет только чисто периодиче­ скую (с периодом Т) часть движения, так как в нашем ре­ шении заведомо «отсечены» те члены полного решения, которые могли бы нарушить эту периодичность.

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

периодическое решение,

мы к о м п е н с и р у е м

н е ­

п о л н о т у с а м о г о

у р а в н е н и я . С этой

нефор­

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

чем полученные выше два других решения.

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

Рассмотрим какой-либо одни из периодов Т, приняв за начало отсчета времени момент исчезновения последнего импульса. В течение рассматриваемого периода колебания

являются свободными и описываются решением

 

q= Acos pt | В sin pt

(0 < / < Т).

(28.22)

Если q„— начальное смещение и q„— начальная

скорость,

то постоянные А и В равны

 

 

 

A = q„

B = - - f ,

(28.23)

223

решение (28.22) можно

записать

в виде

 

q= ft cos pt +

sin pi

(0 < / < T).

(28.24)

Дифференцируя по времени l, находим скорость

 

q= —qtpslnpl + qt cospt.

(28.25)

В конце этого периода, непосредственно перед приложением очередного импульса (т. е. при (= Л * имеем

= q9cos р Т + ^ - sin рТ,

(28.26)

q, = —p ft sin pT + ft cos pT.

(28.27)

Сразу после приложения очередного импульса

смещение

q сохранит свое значение (28.26):

 

ft = ft — ft cos рт + sin рт,

(28.28)

но скорость мгновенно изменится на величину S/a и при

учете выражения

(28.27) составит

 

ft =

S

(28.29)

—Pft sin рт + q„cos pT + —.

Хотя эти выкладки просты, но возможно, что читатель еще не уловил, с какой целью они сделаны. Все станет ясным после следующего решающего шага: вследствие предпола­ гаемой п е р и о д и ч н о с т и процесса величины Яг и должны быть равны соответственно величинам ft и ft, т. е.

ft — Я• cos рТ + ijj- sin рТ,

(28.30)

f t ” —pq9s\npT + f t cos pT+ £ .

(28.31)

Таким образом, мы получили систему двух алгебраических уравнений с двумя неизвестными ft и ft; решив ее, найдем

 

« • - 4 -

<28-32)

и закон движения (28.22) принимает вид

 

? = Z ^ ( sin ^ + c0SP/,c le £ f )

< ° < * < Л -

(28.33)

И все? Да, это окончательное выражение исчерпываю­ щим образом решает задачу. Здесь нет конечных (как в спо-

227

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

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

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

соответственно составляют ql и <£. Конечно, они отлича­ ются от величин дви дв, которые относятся к чисто периоди­ ческому процессу и даны выражениями (28.32). Запишем тождества:

(28.34)

Первые слагаемые правых частей отвечают периодическому (с периодом Т) движению (28.33), а вторые слагаемые по­ служат причиной свободных колебаний (с периодом 2л(р)

9*=(^-slrd*4 ) С°М+(4 “ Я?)sin P f ' <28-35>

Здесь, в отличие от (28.33), время t отсчитывается от под­ линного начала процесса. Таким образом, движение будет описываться надлежаще образованной (с учетом различия в началах отсчета времени) суммой решений (28.33) и (28.35). Еще раз повторим, что вторая часть имеет практическое значение в течение небольшого промежутка времени вблизи начала процесса.

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

в

выражении (28.33).

 

Прежде

всего заметим, что если рТ12=пл (где п — це­

лое), то

J—* <зои амплитуды перемещений стремятся

к

бесконечности, т. е. наступает ударный резонанс. Если

228

м =2л/Т

— угловая частота

приложения

импульсов, то

условие

ударного резонанса

приобретает

вид

 

 

(28.36)

При всех иных соотношениях частот отклонения оказы­ ваются конечными. Наибольшее отклонение системы со­ гласно (28.33) равно

* - . - £ / > + ' * « • 4 - 1

H s" 4

Так как дробь S/(ap) есть максимальное отклонение, вы­ зываемое одним мгновенным импульсом 5, то выражение

Р -

I

(28.38)

 

2 1sin -5£- <•>

можно назвать коэффициентом влияния повторности. На рис. 28.1 изображена зависимость р от отношения частот

/

1

Рис. 28.1. Коэффициент влия­ ния повторности импульсов Q5

О

оi/p. Как видно, здесь возможно неограниченно большое число ударных резонансов (соответственно формуле (28.36)). Наименьшее возможное значение р равно 1/2.

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

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

Первый из способов можно найти в книге Н . К. Снитко «Динамика сооружений» (М.: Госстройиздат, 1960, с. 69). Второй способ был пред­ ложен А. Н. Крыловым (см. его книгу «Вибрация судов»,— М.: ОНтИ,

1936, § 16.2);

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

см. в

книге

Г. С. Горелика

«Колебания и волны» (М.: Физматгиз,

1959,

гл. II,

§7). В книге С. А. Бернштейна *Основы динамики сооружений» (М.:

Госстройиздат, 1938, § 51) тот же способ применен к задаче о действии периодических импульсов конечной продолжительности и дана некая

229

oneпнл суммарного эффекта. Казалось бы, что эту опенку можно путем несложного предельного перехода отнести к случаю действия мгновен­ ных импульсов, однако найденный С. А. Бернштейном результат со­ вершенно ошибочен — в выкладках принято, что при любых целых числах т и п «cos n (2т—1)л=»—1, так как n (2m—1) — нечетное чис­ ло» (с. 153),

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

Четвертый способ был предложен в 1918 г. Дуффингом для общего случая действия периодической вынуждающей силы и с учетом вязкого сопротивления излагается во всех изданиях книги Л. Г. Лойиякского и А. И. Лурье «Курс теоретической механики» (см., например, 6-е изд., М.: Наука, 1983, § 184). Опираясь на приведенное там решение, можно перейти и к случаю действия периодических мгновенных импульсов на систему без трения; впрочем, в указанном издании этот частный слу­ чай рассмотрен и независимо (см. с. 79—81). См. также «Курс строитель­ ной механики» И. М. Рабиновича, т. 2 (М.: Госстройиздат, 1954), где разъяснен и прием дополнительного учета действительно задаваемых начальных условий с помощью выражения (28.35).

«Математическую автобиографию» Н. Г. Чеботарева см. в журнале «Успехи математических наук» (1948, 1. 3, №4 (26)).

В первые два издания нашей книги в фразу, предшествующую вы­ ражению (28.37), вкралась неточность — вместо «наибольшее отклоне­ ние» было написано «амплитуда». Этой неточности, а также еще одному, пожалуй, более серьезному авторскому промаху (в 4 4) была посвящена статья Y. Рironпеан «Sut la stability des systemes elastiques et de quel* ques erreurs» (Atinales des Ecole nationale suplrieure de mecanjque, Nan* tes, I sem„ 1975).

§ 29. Об «обратной» форме дифференциальных уравнений колебаний

Уравнении Лагранжа представляют собой наиболее общую форму записи дифференциальных уравнений коле­ баний механических систем со многими степенями свободы. Пусть рассматривается консервативная система с л степеня­ ми свободы, причем qt— обобщенные координаты, aik— инерционные коэффициенты и сц,— коэффициенты жестко­ сти. Кинетическая и потенциальная энергия равны

г

U - T I I W A*

(29.1)

i.ft

i . к

 

Соответственно уравнения Лагранжа приобретают форму

п

2 (а1кЧк+ сг*<7») “

(/, fc = l, 2, . . . . п), (29.2)

*=]

 

где Qi— обобщенные силы, являющиеся заданными функ­ циями времени (вынуждающие силы).

230

Соседние файлы в папке книги