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

книги / Механика деформируемого твердого тела.-1

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

ся после схода груза с балки и составляет

Р т а х = 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, в можно найти, например, методом фокусов. Во всяком слу­ чае ясно, что оно будет быстро затухать с удалением от концов заданной балки; поэтому для далеких от кон­ цов балки пролетов найденное раньше периодическое решение будет мало отличаться от полного -решения. Впрочем, независимо от этого в задачах с периодической структурой всегда удобно искать решение в виде суммы периодической части и дополнительной части — последняя позволяет точно удовлетворить заданным начальным (или граничным) условиям.