книги / Устойчивость и колебания упругих систем. Современные концепции, парадоксы и ошибки
.pdfприложения л-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 + —. |
Хотя эти выкладки просты, но возможно, что читатель еще не уловил, с какой целью они сделаны. Все станет ясным после следующего решающего шага: вследствие предпола гаемой п е р и о д и ч н о с т и процесса величины Яг и q» должны быть равны соответственно величинам 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