книги / Механика деформируемого твердого тела.-1
.pdfся после схода груза с балки и составляет
Р т а х = V У\ + (У*/к)2 = \ l P l s/ { 3 E J ) ,
где
(21.14)
График этой функции показан на рис. 21.4, б (по оси абсцисс отложены значения v/(kl) = 1/(5). Наибольший прогиб всего на 14% превосходит статическое значение Р/3/(3 £ 7), и, конечно, ни о каком резонансе говорить невозможно.
При формальном использовании выражения (**) для
балки с |
распределенной массой у нас получилось бы, |
что при |
v/(kl) = 1/я должно возникнуть некое «резо |
нансное» состояние. Для упрощенной модели некоторый всплеск значения ц происходит при несколько меньшем значении v/(kl) = 0,20, но главное, конечно, в том, что этот всплеск очень слабо выражен. По-видимому, можно не беспокоиться не только за существующие конструкции этого типа (при прохождении одиночных грузов!), но и за будущие — независимо от того, какие скорости движе
ния грузов окажутся реально достижимыми. |
|
|||||
Из |
(21.14) |
следует, что при больших значениях |
||||
v/(kl) величина ц |
приближенно |
равна |
3(5/8 = 3&//(8г;); |
|||
к этому же результату можно |
прийти |
и независимым |
||||
путем. |
Найдем |
импульс |
вынуждающей |
силы |
||
з р А 2/ л |
vt\ |
, |
|
/0. .очч |
|
|
— 2.-^1 — J |
(см. правую часть |
(21.12)): |
|
|
l/v
о
и теперь примем, что он мгновенно, приложен к концу балки. Тогда вызываемое таким импульсом максимальное перемещение конца балки найдется в виде
_ |
S |
3kl PZ3 |
Ут&х |
тк ~ |
8i7 3E J % |
т. е. оказывается таким же, как и по выражению (21.14) при больших значениях v/(kl).
В работе В. В. Болотина [И] рассматривается действиё гру за, движущегося вдоль, балки, но при учете его инерции; в этом
случае давление груза на балку равно сумме веса и силы инерции груза, а соответствующие уравнения имеют переменные коэффи циенты (в этом состоит существенное отличие от уравнений, полу ченных А. Н. Крыловым). Хотя речь идет как будто о движении одиночного груза, а следовательно, о процессе ограниченной дли тельности, автор, убедившись в устойчивости решения соответству ющего однородного уравнения, пишет: «Нашим дальнейшим ис следованиям будут подлежать исключительно частные решения неоднородного дифференциального уравнения» и затем ищет кри тические значения скорости груза, при которых в частном реше нии возникают резонансные слагаемые. Конечно, этот прием анали за уместен только в условиях периодического воздействия, когда по балке движется регулярная последовательность грузов, а не един ственный груз. В книге В. В. Болотина [12], вышедшей в свет вско ре после публикации [И ], отчетливо оговорено, что упомянутое решение относится имепно к случаю движения серии грузов; по этому поводу см. также статью А. Б. Моргаевского [39].
§ 22. Разумная «непоследовательность» учета трения в решениях некоторых задач
теории механических колебаний
В широко известной книге Литлвуда [27] можно про читать: «При изложении математического рассуждения мастерство заключается в умении дать образованному читателю возможность сразу, не заботясь о деталях, схва тить основную идею; последовательные дозы должны быть такими, чтобы их можно было глотать «с ходу»; в случае неудачи, или если бы читатель захотел что-ни будь проверить, перед ним долита с т о я т ь четко ограни ченная маленькая задача (например, проверить тожде ство; две пропущенные тривиальности могут в совокуп ности образовать непреодолимое препятствие)».
Вряд ли кто-нибудь не согласится с этой мыслью, од нако при ее практическом осуществлении возможны ошибки. Здесь же от промахов никто не застрахован — не так-то просто правильно определить, что окажется три виальным для конкретного читателя, а что требует до полнительных разъяснений. К чему приводит пропуск тривиальностей, о котором писал Литлвуд, можно судить по следующей невымышленной истории.
Во многих книгах разбирается задача о действии пе риодических мгновенных импульсов S неизменного на правления на линейную систему с одной степенью сво боды (рис. 22.1). Как известно, здесь вовсе не обязатель но пользоваться разложением в ряды Фурье — простое решение в замкнутой форме легко получить, исходя из
свойства периодичности искомого процесса установив шихся вынужденных колебаний.
Для этого рассматривается один из типовых промежуткс/В времени (любой) длительностью Г, расположен ный между моментами действия двух последовательных
Р*
V s |
|
|
0 t |
Ît2 |
? |
с/ Г ,
Рис. 22.1. Односторонние периодические импульсы
мгновенных импульсов, причем < £ < T (£4*-> 0). Ведя отсчет времени от начала этого промежутка, можно за писать дифференциальное уравнение* в виде
ш хЛ -сх^ 0 ( 0 < £ < Г ) (22.1)
(обозначения обычные; трение для упрощения не |
уч |
|
тено). |
|
|
Решение этого уравнения можно записать в виде |
|
|
х = хх cos kt + ^ s i n (0 < £ < Г). |
(22.2) |
Здесь к = (с/т)112, а индекс «1» означает, что соответ ствующие величины относятся к началу рассматриваемого промежутка времени (t — tt).
Для скорости получается выражение
V= —kxi sin kt + Vi cos kt. |
(22.3) |
Хотя входящие сюда начальные значения координаты и скорости Xi и Vi заранее неизвестны, выразим через них значения хТ и vTв конце промежутка (при £ = Г) :
хтя= |
Ï?I |
^ 4) |
хх cos кТ + -j sin кТ, |
||
vT = |
— кххsin кТ + vxcos кТ. |
|
'Для начала следующего промежутка времени |
(момент |
|
tz-+ Т) должно быть |
|
Эти условия «сшивания» (припасовывания) очевидны — первое условие означает, что за нулевое время координа та не может получить конечного приращения, а второе непосредственно следует из теоремы об изменении коли чества движения.
Поскольку разыскивается периодическое движение с периодом импульсов Г, состояния системы в моменты
и t2 должны быть одинаковыми:
|
х2= хи |
v2= |
Vi. |
|
|
|
Подставляя |
сюда (22.4) |
и |
(22.5), получаем |
два урав |
||
нения с двумя неизвестными |
и vx: |
|
|
|||
|
ххcos кТ + |
|
sin кТ = |
xlf |
|
|
|
|
|
|
$ |
у,. |
|
—кх, sin кТ + vj cos кТ --------- |
|
|||||
1 |
* 1 |
|
' |
771 |
1 |
|
Из этих уравнений находим |
|
|
|
|
|
|
|
S . |
JcT |
’ |
|
S |
/оо а\ |
■ ^ = 2 ы Г с^ |
Х |
|
|
(226> |
||
и согласно (22.2) получаем решение |
|
|
||||
х = |
-т- j - ^sin |
+ |
clg -Ç -co sw ). |
(22.7) |
Это решение относится к промежутку времени между моментами действия двух любых последовательных им пульсов, но, конечно, всякий раз начало отсчета времени нужно совмещать с началом промежутка. Подчеркнем, что в (22.7) нельзя подставлять £ > Т\
Далее легко анализируются свойства найденпого дви жения (определяются максимальные отклонения, фор мулируются условия резонанса кТ/2 = пп и т. д.).
Приблизительно так излагается это экономное и изящ ное решение в современных книгах и, в частности, в по следнем издании курса [28]. Такое изложение, вполне ясное для подготовленного читателя, из-за отсутствия некоторых комментариев вызывает недоуменные вопросы у начинающего (пропуск тривиальностей!). Среди различ ных недоразумений, притом возникавших не только на студенческом уровне, автору пришлось встретиться со следующими:
1. |
« Н е п р и я т и е » представления о |
мгновенном им |
||
пульсе |
конечной величины, |
поскольку |
из |
такого пред |
ставления следует, что сила |
при ударе |
оказывается 6eç- |
конечно большой (нужно иметь в виду, что задача о дей ствии периодических импульсов рассматривается в курсе задолго до изучения темы «удар»). Если здесь не помо жет ссылка на понятие о дельта-функции Дирака, то для разъяснения полезно сослаться па курс сопротивления материалов и указать на аналогичную замену нагрузки, распределенной на малом участке длины балки, сосре доточенной силой (равнодействующей).
2. В о з р а ж е н и е , связанное с тем, что каждый очереднойщмпульс вызывает скачок скорости и, по-видимому, изменяет кинетическую Энергию системы. Но на проме жутке между двумя импульсами энергия остается посто янной (система консервативная) и, следовательно, на раз личных промежутках энергия неодинаковая, а этого не может быть при периодическом процессе.
Читатель, вероятно, уже заметил, в чем ошибка этого рассуждения: хотя каждый импульс действительно вызы вает скачок скорости, но кинетическая энергия при этом не м е н я е т с я , , потому что абсолютное значение ско рости остается прежним — до приложения очередного им пульса скорость равна —S/(2m), а после этого меняет
знак, |
оставаясь |
той же-по величине S/(2m). |
3. |
Н е у д о в |
л е т в о р е н н о с т ь тем, что решение |
представлено не в привычном виде аналитического выра жения, справедливого для любых значений времени t, а в «кусочном» виде, относящемся только к одному ти пичному промежутку времени. Поскольку в курсах мате матики и механики этому представлению (весьма важно му для практики) уделяется очень мало внимания, не исключено, что учащййся впервые с ним встречается именно здесь. Для разъяснений может быть уместна ссылка также на курс сопротивления материалов, где при построении эпюр постоянно приходится иметь дело с функциями, имеющими различные аналитические вы ражения на разных промежутках изменения аргумента.
4. П р е т е н з и я по поводу того, что полученное пе риодическое решение не удовлетворяет произвольно за даваемым подлинно начальным условиям (для момента, непосредственно предшествующего приложению первого импульса). Формально отвечая на эту претензию, доста точно сослаться на то, что разыскивалось лишь ч а с т н о е решение (чисто вынужденные колебания), а оно и не должно удовлетворять начальным условиям. Но еще луч ше показать, как в надшей задаче фактически достроить общее решение.
Для этого прежде всего определяются следующие из (22.7) значения х и х в подлинно начальный момент, непосредственно предшествующий первому импульсу (они такие же, как и в момент t = Т):
* ( ° > - 2 = |
* « > > - - £ • |
<22-8> |
Если подлинные начальные условия заданы в виде |
|
|
х(0) = х0, |
i(0)=i>o, |
(22.9) |
то дополнительное движение я*(т), которое нужно доба вить к периодическому движению (22.7), порождается различиями между (22.8) и (22.9), а именно начальными условиями
** (°) = |
х0 _ ^ |
Ctg -Ç-, |
X* (0) = г0 + -JU |
(22.10) |
||
Следовательно, дополнительное движение описывается |
||||||
выражением |
|
|
|
|
|
|
** = (Жо - |
-2Ш c l g - f )cos^T + |
( х |
+ |
2§fc) sin |
<22Л1> |
|
в котором |
т — время, |
отсчитываемое |
от подлинного нача |
|||
ла процесса, причем |
т = (/ — l)T + t |
(/ — помер |
проме |
жутка времени, для которого вычисляется я*). Оконча тельно общее решение можно записать в виде
* = (жо — 2§* |
ct^ -y -) cos к [(/ — 1) T + |
t] + |
+ |
(тг + 2§ft)sin А' ^ - Ц Т |
+ 1] + |
+ ~Siï[sinkt + clg^i~coskt)’ (22 12)
(Здесь хочется сразу отметить, что дополнительное дви жение в действительности будет: затухающим из-за влия ния неизбежных сил трения, но это все же другой вопрос,. к которому мы еще вернемся в п. 6.)
5.З а б в е н и е (незнание, недооценка) того, что част
ное решение линейного дифференциального уравнения с периодической правой частью всегда представляет собой периодическую функцию того же периода (возможно, описываемую «кусочным» образом); именно это, удобней шее для анализа свойство положено в основу приведен ного выше решения. Забывая о сказанном, иногда пред
лагают искать общее решение с помощью интеграла Дюамеля в виде
Т
х = xQcos кх + -~-sin кх + J(? (£) sin к(х — £) |
(22.13) |
о |
|
где Ç ( T ) — вынуждающая сила, которая в нашей задаче может быть записана через дельта-функцию Дирака
оо |
|
<2(0 = 5 2 в [ т - ( / - 1 ) Л . |
(22.14) |
3=1
Несомненно, такое предложение продиктовано «луч шими чувствами» и опирается на высокую репутацию интеграла Дюамеля*), который действительно дает уни версальное решение задачи о колебаниях системы под действием произвольной вынуждающей силы. Однако во многих конкретных случаях, когда нетрудно сразу ука зать частное решение, совершенно не обязательно обра щаться к интегралу Дюамеля. Это относится в особенно сти к задачам о периодическом возбуждении; здесь и удобно, и естественно выделить периодическое частное решение (с периодом правой части). В этих случаях ре шение в форме (22.13) потребует, бесспорно, более гро моздких выкладок.
Выделение периодической части решения всегда ради кально упрощает выкладки — не только в задачах дина мики, но и в статических задачах, где аргументом служит, пространственная координата и речь идет* о периодично сти в.пространстве (об этом см. в конце параграфа).
6. С о м н е н и е в самостоятельной ценности найденного решения, поскольку к нему нужно еще что-то добавлять, чтобы результат удовлетворял начальным условиям. Но дополнительное движение в действительности оказывает ся затухающим, и в этом мы убедимся, если учтем трение и при определении этого движения будем исходить не из уравнения (22.1) , а из уравнения
шх + Ъх+ сх = 0. |
(22.15) |
Поэтому довольно сцоро после начала процесса дополни тельное движение станет практически неощутимым и ре
*) Жан Мари Констан Дюамель (1797*—1872) — французский математик, член Парижской академии наук (с 1840 г.), член-кор респондент Петербургской академии наук (с 1859 г.). Интеграл Дюамеля был предложен в 1834 г.
шение (22.7) достаточно точно опишет полное движение. Иными словами, частное решение (22.7) нужно толковать как «асимптотически верное».
Но тогда возникает новый (пожалуй, наиболее обо снованный) вопрос: почему роль дополнительного дви жения исследуется па основе уравнения (22.15) с учетом сил трения, тогда как искомое частное решение (22.7) получено, исходя из уравнения (22.1), в котором трение не отражено? Почему трение признается существенным фактором на одном этапе решения и им пренебрегается на другом? Нет ли здесь непоследовательности?
Конечно, если речь идбт о действии относительно ред ких импульсов, когда за время Т система успевает со вершить несколько полных колебаний (с собственной частотой к), трение следует учесть и при построении ча стного решения. Отметим, что такое решение (и притом для общего случая действия периодической вынуждаю
щей силы*)) существует и опирается |
на те |
же сообра |
|||||
жения |
о |
периодичности |
(см. [28], |
§ |
184). |
При этом, |
|
в частности, вместо (22.6) можно |
получить |
следующие |
|||||
выражения: |
|
|
|
|
|
||
|
|
= _S______ g”rsinКТ |
|
|
|
||
|
|
Х° ~~~тпк*• |
е2пТ_ 2епТcos 7г*Г + |
1 1 |
|
||
|
|
s |
епТ — cosк*Т — |
sinк*т |
|
||
|
|
V° ~ ™ |
е2пТ — 2епТ соз k* T + |
1 |
|
||
(здесь |
п = |
Ъ/(2т), к* = |
(к2— п2)1/2), |
которые |
позволяют |
||
построить,частное решение, учитывающее трение: |
|||||||
|
|
|
vn+ пхп |
\ |
(0<Zt<iT). |
||
|
|
( х0 cos k*t -1— ^ — - sin k#tj |
Такое решение, в частности, дает возможность определить конечные размахи колебаний при резонансе.
Более того, если импульсы весьма редкие (Т >> 2тс/к#), то за 'время Т действие предыдущего импульса практи чески полностью затухает и речь может идти о приложе нии одиночных импульсов.
С другой стороны, если импульсы достаточно частые
(Т <С 2л//г*), то |
за время |
Т трение практически пе ус |
пеет сказаться |
и решение |
(22.7) окажется достаточно |
) Опо было предложено Дуффтшгом в 1918 г.
точным. В этих случаях оно имеет самостоятельную цеп кость, и им можно пользоваться в достаточном удалении от подлинного начала процесса.
Что касается формальной непоследовательности, то ведь «непоследовательности» такого типа можно обнару жить чуть ли не во всех прикладных теориях, когда один и тот же фактор, учитываемый на одних этапах решения, не учитывается на других. Вспомним, например, «прин цип неизменности начальных размеров», которым поль зуются почти всегда в сопротивлении материалов при за шей уравнений равновесия; здесь деформации считаются отсутствующими, хотя на других этапах решения они обязательно учитываются. Можно привести пример и яз теории колебаний. Известно, что далеко от резопанса ре жим вынужденных колебаний можно определять без уче
та трения, но т;ама |
ориентация анализа на э т о т ре |
жи м опирается на |
невысказываемое предположение о |
затухании роли начальных условий и потому косвенпо содержит своеобразный учет трения.
Рис. 22.2. Многопролетная балка периодической структуры: а) за данная схема; б) схема к определению периодической части реше ния; в) схема к определению дополнительного решения
Выделение периодической части решения удобно и при анализе статических задач о деформировании систем с периодической (в пространстве) структурой. Пусть, например, речь идет об изгибе многопролетной неразрез ной балки с равными и одинаково нагруженными про летами (рис. 22.2,а). В данном случае целесообразно4ис кать решение в виде суммы периодического решения, описывающего изгиб любого пролета по схеме, показан
ной па рис. 22.2, б, и дополнительного решения, соответ ствующего схеме нагружения, данной на рис. 22.2, в.
Моменты, показанные на рис. 22.2, б, одинаковы по величине вследствие условий периодичности. Согласно тем же условиям одинаковыми должны быть и углы поворота сечений ф па концах выделенного пролета — это и позволит найти значение М (т. е. выразить М через параметры пролетной нагрузки).
Дополнительное решение для схемы па рис. 22.2, в можно найти, например, методом фокусов. Во всяком слу чае ясно, что оно будет быстро затухать с удалением от концов заданной балки; поэтому для далеких от кон цов балки пролетов найденное раньше периодическое решение будет мало отличаться от полного -решения. Впрочем, независимо от этого в задачах с периодической структурой всегда удобно искать решение в виде суммы периодической части и дополнительной части — последняя позволяет точно удовлетворить заданным начальным (или граничным) условиям.