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

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

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

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

Важно отметить, что в этих вполне обычных выкладках предполагается, что величина р не равна нулю. Почему от­ брасывается возможность того, что р = 0? Есть ли для этого основания, в чем они заключаются и насколько они универ­ сальны? Почему отброшены нулевые корни а /= 0 , удовлет­ воряющие характеристическим уравнениям для схем 1, 2,

4,

5?

 

 

Предположим с самого начала, что р—0. Тогда вместо

уравнений (26.4) и (26.5) соответственно получится

 

 

f = 0 ,

(26.8)

 

V"v = 0 .

(26.9)

Решение уравнения (26.8) можно записать в виде

 

 

Г = Л - ) - 0 /,

(26.10)

а

решение уравнения (26.9) — в виде

 

 

У= Ct + C,z + C„z* -I- C4z*.

(26.11)

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

V (0 ) -0 , Vf {0 ) - 0 , V (/) = 0, У '(О - о,

то получится однородная

система

уравнений

С ,- 0 ,

 

Сг =

0,

с ,- г с гн

-

о

,

С, + 2С:1/ + ЗС4Р = 0,

определитель которой не равен нулю. Это означает, что все постоянные С{ равны нулю, т. е. движение предполо­ женного вида попросту н е с у щ е с т в у е т . Аналогич­ ным образом можно убедиться в том, что случай р—0 дол­ жен быть отброшен и для схемы 2. хотя а / —0 формально удовлетворяет характеристическому уравнению (см. таб­ лицу).

Если же рассмотреть схему 4, то результат получится иным. Для этой схемы граничные условия имеют вид

V (0) — 0, V"

Г ( /) - 0 , Г " (0 = 0.

2'Л

После подстановки сюда выражения (26.il) получим

V=Ctz,

и согласно (26.3) движение описывается выражением

v={A + В1)г,

в котором постоянная С2 включена в А и В. Полученное движение представляет собой равномерное вращение жест­ кой балки вокруг левого шарнира; такое движение в самом деле не исключено, хотя оно и не носит колебательного ха­ рактера.

Итак, нулевому значению собственной частоты отве­ чает движение балки без деформаций. Это — общее пра­ вило: если система допускает перемещения частей как жестких тел, то возникают нулевые собственные частоты. При этом число таких частот равно числу степеней свободы механической системы, которая получится из заданной си­ стемы, если ее звенья считать недеформируемыми. В схе­ мах 1—3 общая подвижность отсутствует; поэтому здесь нет и нулевых собственных частот. Схеме 4 соответствует одна нулевая собственная частота, а схеме 5 — две такие частоты.

Наличие нулевых корней следует и из формального рас­ смотрения характеристических уравнений для схем I, 2, 4, 5, но в первых двух случаях эти корни должны быть от­ брошены, а в двух других случаях они имеют смысл. Таким образом, совпадение характеристических уравнений не означает полного совпадения спектров корней (и собствен­ ных частот).

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

Исправим последние строки данной выше таблицы, вне­ ся в них пропущенные нулевые корпи (см. с. 203).

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

Во мно1их книгах можно найти формулу

202

Первые три корня характеристического уравнения

 

 

 

схемы

W h

«*1>,

«Л»

 

i

4,730

7,853

10,996

2

3,927

7,068

10,210

3

1,875

4,694

7,854

4

0,000

3,927

7,068

5

0,000

0,000

4.730

определяющую собственную частоту крутильных колебаний двухмас­ совой системы, показанной на рис. 26.1. В этой формуле У, н У, — мо­ менты инерции дисков относительно оси системы, I — длина упругого

Рис. 26.1. Крутильно-колебательная сис­

*771

тема

 

вала,, GJр — жесткость вала при кручении, Эта формула, конечно, верная, но в сущности она определяет не п е р в у ю собственную частоту, а в т о р у ю: первая собственная частота равна нулю и соот­ ветствует вращению системы вокруг оси как жесткого целого.

Теперь мы переходим к другому рассказу, в котором речь пойдет тоже о «потере» первой собственной частоты, но н е р а в н о й н у л ю (тема подсказана авторам О. А. Са­ виновым).

Будем рассматривать свободные продольные колебания стержня, находящегося в упругой среде, которая создает реакции, пропорциональные продольным перемещениям. Обозначив эти перемещения через w(z, t), запишем интен­ сивность продольных реакций среды в виде —т, где г — коэффициент жесткости среды. При этом вместо обычного дифференциального уравнения продольных колебаний

Е SPw dho р дгг dt2

и р — соответственно модуль упругости и плотность материала стержня) получится несколько более сложное дифференциальное уравнение

£

г _

(26 12)

р* dz2

р£ w ~ ~дЛ

 

203

( f — площадь сечения стержня). Разыскивая частное реше­ ние в виде w-=W(z)T(z) и разделяя переменные, приходим к дифференциальным уравнениям

Г + р*Т = 0,

(26.13)

W"-)-а*Г = 0,

(26.14)

где

 

<2815)

Решение уравнения (26.13)

Т = A sin (pt + у)

показывает, что движение представляет собой гармониче­ ские колебания с угловой частотой р\ решение уравнения (26.14)

W ^ C ^ in a z + C jcosaz

(26.16)

вместе с граничными условиями позволит найти значения р. Остановимся на случае, когда левый конец стержня (z=0) свободен, а на правом конце (г=1) со стержнем связана со­ средоточенная масса /И. При этом граничные условия на обоих концах будут относиться к продольной силе, кото­ рую можно записать в виде

N = EF%j~= EFW'T.

(26.17)

Так как при г—0 продольная сила тождественно равна нулю, то первым граничным условием будет служить соот­ ношение 1Г'(0) = 0; при этом согласно (26.16) С,—0 и, сле­ довательно, дальше можно исходить из более короткой записи

W =Сгcosotz.

(26.18)

На правом конце стержня продольная сила равна силе инерции груза, которую можно представить в следующем виде:

— м Щ г_ = — М \ Р ( 1 ) Т = М Ч Р ( 1 ) р * т =

=+ (26.19)

Приравнивая (26.17) и (26.19) при г=1, получим второе гра­ ничное условие

EFW ф — MW (I) + £ ) ,

(26.20)

204

которое при учете выражения (26.18) приводит к харак­ теристическому уравнению

 

 

 

 

 

 

 

 

 

 

 

(26.21)

Здесь

m—M'ipFl) — относительная

масса

груза. е=*

—rlV (EF) — безразмерный

коэффициент жесткости

среды.

Прежде всего

отметим,

 

i W .

1

 

 

что если

полученное урав­

 

 

 

нение имеет корень + (ос/)*,

8

 

1

 

 

то имеется также отрица-

 

 

 

 

тельный

корень — (а/)*;

6

J

 

 

 

обоим этим корням соглас-

 

 

 

но (26.15) отвечает одно и

t,

 

 

 

то же значение собствен-

 

 

 

ной частоты. Поэтому даль­

 

 

 

 

ше можно не расематри-

?

 

 

 

вать отрицательные значе­

Q

 

 

 

ния а 1.

 

 

 

 

 

 

 

'

Если построить графики

 

гс

 

 

обеих

частей

уравнения

 

 

 

l

 

 

(26.21),

то абсциссы точек

 

//

 

4

c

 

пересечения графиков опре­

 

 

 

делят корни уравнения. В

 

7

 

j

~a l i i

 

качестве

примера

это сде­

 

 

 

лано на рис. 26.2 для

 

i

 

ftg e/

 

т= 1 я

е= 1 . О|лошиой

и

 

i

 

 

штрихпунктирной линиями

 

 

 

 

 

показаны графики соответ­

 

i

 

 

 

 

ственно

 

левой

и

правой

 

 

 

 

 

частей

уравнения

(26.21).

Рис.

26.2.

Графическое

решение

Как видно, в интервале

 

 

уравнения

(26.21)

 

значений

0, я

существует

 

 

ее

абсцисса а / = 1,956

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

графиков;

согласно (26.15) определяет наименьшее значение соб­ ственной частоты

(26.22)

Если уменьшать относительную массу груза, то штрихпунктирная линия будет несколько изменяться, причем точка ее пересечения со сплошной линией будет переме­

щаться в направлении, указанном на рисунке

стрелкой,

т. е. собственная частота будет в о з р а с т а т ь .

Пока в

203

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

Если на правом конце нет груза, то вместо граничного

условия (26.20) нужно

принять

W(l)=0. Эго приводит

к характеристическому

уравнению

sin а /= 0 ; наименьший

корень at—0 позволяет согласно (26.15) найти первую соб­

ственную

частоту

 

 

Р -

(26-23)

которая

оказывается м е н ь ш е й ,

чем по выражению

(26.22) . Вот это уже странно: устранение концевого гру­ за должно было повлечь за собой увеличение собственной частоты.

Возникает естественное подозрение, что выражением (26.22) определяется не первая, а в т о р а я собственная частота и, следовательно, выражения (26.22) и (26.23) нельзя сравнивать. Но где же была «упущена» п е р в а я собственная частота для стержня с грузом?

Если внимательно обдумать наши выкладки, то можно заметить одно сомнительное место: при решении трансцен­ дентного уравнения (26.21) мы считали, что его корни а / должны быть в е щ е с т в е н н ы м и .

Однако из соотношения (26.15) видно, что отрицатель­ ность а* (и соответственно мнимость а), вообще говоря, мо­ жет не противоречить положительности />*, и поэтому впол­ не возможна. Поэтому нужно проверить, имеются ли кор­ ни характеристического уравнения в мнимой области.

Считая, что а г< 0 , удобно ввести

положительную

величи­

ну р*=—а*; тогда а —ф, tg a l — i

th р/ и характеристиче­

ское уравнение (26.21) приобретает вид

 

t h p / = m ( - p / - f - - p ) .

(26.24)

На рис. 26.3 представлено графическое решение этого уравнения для принятых выше значений постоянных т= 1, е= 1 . Здесь сразу обнаруживается корень р/=0,735, кото­ рому соответствует (р/)*^0,540 и (al)*——0,540. Подставляя этот результат в соотношение (26.15), находим

 

/>=0,678

Y j p -

(26.25)

Это и

есть правильное значение собственной частоты, .

именно с

ним можно сравнить

результат (26.23). Как и

2%

предполагалось выше, результат (26.22) относится ко вто­ рой собственной частоте. Читатель может проверить, что первым собственным частотам (26.23) и (26.25) соответствует несколько особенная собственная форма — перемещение всего стержня как жесткого целого вдоль своей оси.

В связи с переходом от уравнения (26.21) к уравнению (26.24) вспоминаются образные слова французского мате­ матика Ж. Адамара: «Наиболее короткий путь между двумя

истинами в вещественной области часто проходит через мни­ мую область».

Таблица, приведенная в начале параграфа, заимствована со с. 356 книги А. Н. Крылова «Вибрация судов» (Л.: М.: ОНТИ, 1936). В том же виде ее можно найти на с. 77 книги А. П. Филиппова «Колебаимя деформируемых систем» (М.: Машиностроение, 1970).

§27 . Медленно изменяющиеся силы

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

207

приводятся и далее, а понятие «медленности» так и ос­ тается чисто качественным.

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

Яркий пример указан А. Н. Крыловым *): хотя при вы­ стреле из артиллерийского орудия длительность действия давления газов на стенки орудия составляет всего 0,01 с (для рассмотренного А. Н. Крыловым случая), однако при изучении радиальных колебаний ствола такое нагруже­ ние можно считать статическим! Наряду с этим можно было бы привести противоположные примеры, когда нужно считать динамической нагрузку, которая монотонно воз­ растает и затем монотонно в течение нескольких секунд убы­ вает до нуля.

В настоящем параграфе мы дадим некоторые количе­ ственные оценки и начнем с простого случая, когда меха­ ническая система с о д н о й с т е п е н ь ю с в о б о д ы нагружена длительно действующей гармонической вынуж­ дающей силой Рфsin (о/. Обозначим через q{t) обобщенную координату, т. е. характерное перемещение системы в процессе установившихся вынужденных колебаний; тогда можно получить

 

q—A sin

,

(27.1)

где

 

 

 

 

 

 

(27 2)

1— амплитуда

вынужденных

колебаний,

с — коэффи­

циент жесткости

системы, р — ее собственная

частота.

*) Алексей Николаевич Крылов (1863— 1945) — механик, мате­ матик, кораблестроитель и историк науки. Академик (с 1916 г.). Герой Социалистического Труда (с 1943 г.). С 1890 г. в течение почти пятиде­ сяти лет преподавал в Морской академии в Петербурге (Ленинграде).

Отношение амплитуды колебаний к перемещению PJc, вызванному статически действующей силой Р» (т. е. при <о—*-0),

_ с\ А \

|

I

|

(27.3)

Р,

~ |

1—<о*/Р*

| ’

 

называется коэффициентом динамичности. При малых значениях отношения сolp коэффициент (л мало отличается от единицы. £сли, например, пренебрегать динамическими эффектами, составляющими меньше 5% от статических, то статической можно считать такую силу, при действии кото­ рой коэффициент динамичности Не превосходит значения 1,05. Из этого условия можно с помощью (27.3) найти усло­ вие «статичности» действующей силы co//?<0,218 или

Т/Тя> 4,59,

(27.4)

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

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

Pot 2р ■ лю EF ж о g 2р 1

Здесь Д= тг 1 /

— — наименьшая собственная частота, Е и

2 Г

р

р — соответственно модуль упругости и плотность мате­ риала стержня, F — площадь сечения, / — длина. Так как статически действующая сила Ра вызывает перемещение конца PJ/ (ЕЕ), то коэффициент динамичности равен

р =

Если и в данном случае исходить из допустимости пяти­ процентной погрешности, то для частоты «медленной» силы получится условие со/р<;0,239, для периода — неравенство

Т/Тс > 4Д8,

(27.5)

которое не слишком отличается от условия (27.4).

200

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

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

Рис. 27.1. а) Заданное изменение силы; б) изменение силы в первой вспомогательной задаче; в) изменение силы во второй вспомогательной задаче

сначала возрастающая по линейному закону Р=()7, а затем при £><* остающаяся постоянной. По-видимому, при достаточно больших значениях времени силу можно считать статической, а при весьма малых значениях /„ — внезапно приложенной; поэтому при постепенном убыва­ нии ** коэффициент динамичности будет увеличиваться оче­ видно, от единицы до двух.

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

движения имеет вид

 

а<7+ е<? = р/

(27.6)

(а — инерционный коэффициент). Если записать общее решение и затем подчинить его нулевым начальным услови­ ям, то в результате получится для £>0

р/

р sin pt

(27.7)

С

ср

 

Этим результатом можно воспользоваться и во второй вспомогательной задаче, когда внешняя сила меняется так, как показано на рис. 27.1, в. Для этого нужно в (27.7) изменить знак q и вместо t подставить t—1*. Тогда для /> /* получится

- Г

Р и

-<»)

Р sinр (/ /,)

(27.8)

"

I

с

~ёр

 

210

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