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