книги / Прочность, устойчивость, колебания. Т. 3
.pdfОбласти неустойчивости для систем уравнение |
361 |
найдется хотя бы один корень, обладающий положительной действи тельной частью. Исследование уравнения (37) показывает, что области неустойчивости матричного уравнения (36) располагаются вблизи частот
|
0 = ^ - (/1= 1, 2,. |
. ); |
(38) |
здесь |
— частоты собственных колебаний системы, загруженной ста |
тической нагрузкой с параметром а. Эти частоты определяют из урав нения
|
|
|
|
\Е — аА — Й2С | = 0. |
|
|
|
|
|
|
(39) |
|||
Кроме того, уравнение (37) позволяет |
выделить две группы |
«подо |
||||||||||||
зрительных» |
частот. Первую |
группу |
частот определяют |
по |
формуле |
|||||||||
|
|
|
д _ _ 0 / + Оа |
(/1 = |
1, 2,. |
!ф к), |
|
|
|
|
(40) |
|||
|
|
|
|
п |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
вторую группу частот |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
0 |
0 / — |
(п= 1, 2 . |
|
|
|
|
|
|
(41) |
||
|
|
|
— |
п |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Из |
более |
строгого |
анализа |
выражений |
[14, |
20, |
21] |
следует, |
что |
|||||
в случае гамильтоновых систем * области неустойчивости |
|
образуются |
||||||||||||
лишь вблизи |
частот, определяемых |
по формуле (40). Если |
система |
|||||||||||
является негамильтоновой, то возможны области |
неустойчивости, |
|||||||||||||
располагающиеся вблизи второй группы частот. |
|
|
|
|
|
|
||||||||
|
Области |
неустойчивости, лежащие вблизи частот, соответствующих |
||||||||||||
формуле (38), называют основными областями |
[5]; |
остальные |
обла |
|||||||||||
сти — комбинационные. Последнее название подчеркивает |
тот |
факт, |
||||||||||||
что колебания внутри этих областей осуществляются главным |
образом |
|||||||||||||
за счет взаимодействия какой-либо пары форм |
колебаний. |
Это |
непо |
|||||||||||
средственно следует из формул (40) и (41), |
в которые |
в симметричной |
||||||||||||
форме входят две собственные частоты й ; и |
Аналогично можно го |
|||||||||||||
ворить об основных и комбинационных параметрических |
|
резонансах. |
||||||||||||
|
Метод построения границ основных областей неустойчивости [7]. |
|||||||||||||
Границы |
областей, |
которым |
в формуле (38) отвечают значения |
п = |
||||||||||
= |
1, 3, |
|
определяют из |
уравнения |
|
|
|
|
|
|
|
|
||
Е - а А ± 1 р я --\-№ С |
|
|
|
|
О |
|
|
|
|
|
||||
|
|
— |
5 -1» |
г - < м — | е * с |
- |
4 -№ |
|
|
: |
=о. |
||||
|
|
|
О |
|
— |
В |
|
Е - а А - ^ & |
‘С'. |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
' |
(42) |
* Механическую систему с конечным числом степеней свободы называю! гамильтоновой, если уравнения ее движения могут быть представлены в кано нической форме Гамильтона. К нсгамильтоновым системам приходим, напри
мер, рассматривая колебания стержня, нагруженного следящими силами [8]
362 Параметрические колебания упругих систем
Определитель, стоящий в левой части, следует понимать в том смысле, что на месте каждого его элемента стоит матрица порядка п, где п — число членов ряда (11). Например, удерживая в определителе (42) одни лишь выписанные члены, получим определитель порядка 3п. Границы главной области неустойчивости (п = 1) в первом прибли жении могут быть найдены из уравнения 17]
|е — аА ± 4-РД---^-еас|=0. |
(43) |
Границы областей, которым в формуле (38) отвечают значения п — |
|
= 2, 4, . . ., находят из уравнений, аналогичных выражению |
(42). |
Эти |
уравнения неудобны для аналитических вычислений, поскольку |
их |
решение требует развертывания определителей высокого порядка |
и отыскания корней алгебраических уравнений высоких степеней. Эти операции, однако, не представляют затруднений для электронных цифровых машин.
Применение метода малого параметра. Для аналитических вычисле
ний более удобны формулы, |
основанные на методе малого параметра |
|||||||||||
[21, |
23]. |
Преобразуем |
уравнение |
(36) |
к |
главным |
осям |
матрицы |
||||
С ^ ( Е - а А ) : |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
+ |
Р Я Л |
СЮ 0/ = 0; |
|
|
|
(44) |
||
здесь / г — преобразованный вектор / ; |
матрица С^ 1 имеет вид |
|
|
|||||||||
|
|
|
|
0 |
|
|
0 |
|
|
|
|
|
|
|
|
0 |
&\ |
|
0 |
|
|
|
|
|
|
|
|
|
0 |
0 |
|
|
|
• _ |
|
|
|
|
|
|
|
• |
- |
|
|
|
|
|
|
|
|
а матрицу И вычисляют применением преобразования |
подобия к ма* |
|||||||||||
трице С_1Д. Иначе, |
И — |
|
|
|
|
|
|
(45) |
||||
|
|
|
|
|
|
|
|
|
||||
где V — матрица, составленная по столбцам из собственных векторов |
||||||||||||
матрицы |
С-1 {Е — аА). |
|
1. Если |
параметр р достаточно мал, |
то |
|||||||
Ограничимся случаем п = |
||||||||||||
границы главных областей неустойчивости могут быть найдены |
по |
|||||||||||
формуле |
[28] |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 .= 2 О / ± |
- |- .- |Щ , |
|
|
|
(46) |
||||
где |
Л/д— элементы матрицы |
Н. Границы |
комбинационных |
областей |
||||||||
определяют по формуле |
|
|
|
|
|
|
|
|
|
|
||
|
|
0, = |
П/ + |
а * ± |
А |
|
] |
/ » . |
|
|
|
(47) |
Если |
диагональные |
элементы матрицы |
Нвелики |
по |
сравнению |
|||||||
с побочными элементами, то |
ширина |
комбинационных |
областей |
не |
||||||||
устойчивости будет мала по сравнению с шириной главных |
областей. |
Влияние демпфирования на области неустойчивости |
363 |
Примером, где это условие не выполнено, может служить задача об устойчивости плоской формы изгиба (см. 354—355). Пусть в уравне ниях (20) и (22) М (0 = рсоз 0/. Тогда уравнение (22) имеет вид (44), где С0 = С,
•Ч 9 1-а>1
И = С”1А =
Ык
Ширина главных областей имеет в этом случае порядок р2. Границы комбинационной области согласно формуле (47) определяют как
где Р* — критическое зи чеиие параметра для соответствующе;’ ческой задачи:
ВЛИЯНИЕ ДЕМПФИРОВАНИЯ НА ГРАНИЦЫ ОБЛАСТЕЙ НЕУСТОЙЧИВОСТИ
Уравнения для границ областей неустойчивости. Ограничимся про стейшим случаем уравнения Матье (26) с диссипативным членом
^ 1 |
+ 2 * 4 - + П2 (1 — 2ц со5 00 = 0; |
(48) |
л 2 |
+ г а |
|
здесь е — коэффициент демпфирования. Области неустойчивости попрежнему располагаются вблизи частот, определяемых по формуле (2/). При п = 1,3, границы областей неустойчивости находят из урав нения [7]
. |
902 |
|
— ц |
0 |
- |
А |
. Л |
|
4Й2 |
|
|
|
|
я |
2Й |
— И |
, , |
02 |
Д . |
о |
|
0 |
|
|
|
я ' 2Й |
|
||||
0 |
|
_д___ о_ |
, |
02 |
|
= 0. (49) |
|
|
_ р |
||||||
|
л |
20 |
1 11 |
4ЕГ- |
|||
А |
39 |
|
0 |
- И |
|
|
902 |
я ’ |
20 |
|
|
|
|||
|
|
|
|
|
402 |
||
|
|
|
|
|
|
|
Влияние перемещений в невозмущенном состоянии |
365 |
Влияние демпфирования в общем случае. Если параметрические ко лебания описываются системой дифференциальных уравнений с периоди ческими коэффициентами, то расчет областей неустойчивости при на личии демпфирования усложняется. Некоторые приближенные методы указаны в книге [71. Метод малого параметра дает для определения границ комбинационных областей следующую приближенную фор мулу [28]:
О = й / |
± |
6/7 + |
екк |
Е А А / |
|
(52) |
|
|
2 |
|
|
|
|
здесь ^ — элементы |
матрицы |
демпфирования |
после |
при |
||
ведения уравнения |
|
|
|
|
|
|
|
|
+ ( Е - а Д |
- Р с°8 9 ^ ) / = 0 |
|
|
|
к главным осям матрицы |
С-1 (Е — аА). Из формулы (52) видно, |
что |
при определенных условиях введение малого демпфирования может приводить к расширению комбинационных областей неустойчивости. Аналогичное явление было обнаружено в задачах устойчивости упру гих систем, находящихся под действием непотенциальных сил [8 ].
ВЛИЯНИЕ ПЕРЕМЕЩЕНИЙ В НЕВОЗМУЩЕННОМ СОСТОЯНИИ
Постановка задачи. Рассмотренные выше задачи параметрических колебаний можно трактовать как задачи об устойчивости некоторых
режимов |
установившихся вынужденных |
|
колебаний. |
Поясним |
это |
на примере задач, показанных на рис. 1. |
В случае, |
показанном |
на |
||
рис. 1, а, |
роль невозмущенного движения |
играют продольные колеба |
ния стержня, в случае рис. 1, 6 — радиальные колебания кольца, в случае 1, в — колебания пластинки в своей плоскости и т. д. Однако весь предыдущий анализ базировался на предположении, что переме щения в невозмущенном состоянии пренебрежимо малы. Рассмотрим уточненную постановку задачи для случая упругого стержня, сжимае мого периодической продольной силон (рис. 3).
Пусть и (х, I) — продольное перемещение точек; V(*, 0 — попереч ное перемещение точек, принадлежащих оси стержня; ЕР — жесткость сечения при растяжении-сжатии. С учетом наиболее существенных не
линейных |
членов |
уравнения |
совместных |
продольных и поперечных |
|||||
колебаний |
имеют |
вид |
[6, 7, |
|
13] |
|
|
|
|
|
|
|
- ^ 4 |
|
г + " ‘ 4 т г = 0; |
|
|||
|
|
|
|
дх |
1 |
д^^ |
|
|
|
|
|
д1и |
сс |
|
д |
( |
дь \ |
д2у |
(53) |
|
Е* ~ д ^ ~ ЕР дх ( 8х |
дх ) |
+ '" д(2 ~ ° |
||||||
|
|
Через гх обозн чена продольная деформация оси стержня, т. е.
ди . 1 / до \ а
366 Параметрические' колебания упругих систем
Решение уравнений (53) должно удовлетворять-граничным условиям
|
V— |
г |
и = |
0 при Хг= 0; |
|
(54) |
|
= |
|
||||
у = — |
= 0 ; |
— Ерел = |
Р0 + Р( соз0* при * = /. |
|
||
Исследование |
устойчивости. Уравнения (53) |
граничные |
усло |
|||
вия (54) будут удовлетворены, |
если положить |
|
|
_______________ |
Решение (55) |
описывает установившиеся про- |
|||
о |
м |
дольные колебания; это |
решение будем |
назы- |
|
|
рис, ю |
вать- невозмущенным-. Исследуем устойчивость |
|||
|
|
этого решения. |
Полагая |
в уравнениях |
(53) |
и = иа-\- ы,. и — о0+ Ъ, где и и о — малые возмущения, и линеари зируя уравнение (53),относительно этих возмущений, придем к следую щим уравнениям в вариациях:.
Уравнения (56) подробно исследованы в книге [7]. Важнейший ре зультат, относящийся к главной области параметрического резонанса, состоит в следующем: помимо резонанса вблизи соотношения 0 = 2Й, возможно возбуждение поперечных колебаний 0 = ш/ (со/ — частота продольных колебаний).
Физическое истолкование результата. Вблизи 0 = со/ имеет место резонанс продольных колебаний, вследствие чего резко возрастает динамическая продольная сила в стержне. Поэтому жесткость стержня по отношению к поперечным колебаниям вблизи 0 = со; является перио дической функцией времени с большой амплитудой изменения. Глав ные области параметрического возбуждения при отсутствии демпфиро вания показаны на рис. 10. При больших значениях коэффициента возбуждения р области сливаются.
Дальнейшие подробности можно найти в книге [7]. Там же рас смотрен в уточненной постановке ряд других задач (параметрические колебания криволинейных стержней, балок, рам, пластин и т. д.).
Учет нелинейнык факторов |
367 |
УЧЕТ НЕЛИНЕЙНЫХ ФАКТОРОВ ПРИ ПАРАМЕТРИЧЕСКИХ КОЛЕБАНИЯХ
УПРУГИХ СИСТЕМ
Общие замечания. Рассмотрение параметрических колебаний в ли нейной постановке позволяет найти границы областей неустойчивости и описать поведение упругих систем в течение начального периода возбуждения параметрических колебаний. Согласно линейной теории амплитуды параметрических колебаний возрастают со временем по экспоненциальному закону. Для того чтобы найти амплитуды устано вившихся колебаний, необходимо рассмотреть задачу в нелинейной постановке, удерживая в уравнениях члены, которые обычно (напри мер, при изучении вынужденных колебаний) игнори руются.
Простейшая нелинейная задача. Нелинейные пара метрические колебания упругих систем рассматривались впервые в работах [4, 101 и позднее в статьях [6, 16, 25, 28, 29]. Приведем некоторые сведения, относя щиеся к простейшей задаче о колебаниях опертого стержня, сжатого периодической продольной силой. Пусть на подвижной опоре стержня (рлс. 11) имеется линейная упругая связь, препятствующая продольным перемещениям, продольный линейный демпфер и сосре доточенная масса. Полагая, что
V {х, 0 = / (0 51П - р
*(х,и
получим, что колебания системы приближенно описы ваются нелинейным дифференциальным уравнением
• ж + 2е ■§■+ 22 (* ~ 2* с05 е'> / +
+ 1> (/■ |
Й . |
Щ = о |
(57) |
<И' |
АР ) |
|
В уравнении (57) использованы те же обозначения, что и в урав нении (48); кроме того, в это уравнение входит нелинейная функция
* = V/3 + 2е„^ 1 + 2 ,/ [ / - $ + ( 4 ) 2] ! |
(58) |
368 Параметрические колебания упругих систем
здесь у. ро и х — некоторые постоянные коэффициенты. Первый член в формуле (58) учитывает влияние продольной упругой связи, второй — влияние продольного демпфера, третий — влияние сосредоточенной массы, закрепленной на подвижном конце.
Графики зависимости амплитуды а0 от частоты возбуждения 9 схематически показаны на рис. 12, а—в в окрестности главного параметрического резонанса.
Рассмотрены три случая:
а) у > 0, е0 = к = 0; б) е0 > 0, у —х = 0, в) и > 0, у = еп — 0*
Штриховыми линиями обозначены неустойчивые решения (подробнее см. работы [7,10]). Там же дано решение различных нелинейных задач о параметрических колебаниях стержней, стержневых систем, пласти нок. Параметрические колебания упругих систем в нелинейной поста
новке рассмотрены также |
в работах [16, |
25, 27]. |
|
|
ЛИТЕРАТУРА |
|
|
1. Б е й л и н Е. А., |
Д ж а н е л и д з е |
Г. Ю. |
Обзор работ по дина |
мической устойчивости упругих систем. Прикладная |
математика и механика. |
Т.6. Вып. 5. 1952.
2.Б е л я е в Н. М. Устойчивость призматических стержней под дей
ствием переменных продольных сил. Сб. «Инженерные сооружения и строи
тельная механика». |
Изд-во |
«Путь», 1924. |
|
3. Б о д н с р |
В. А. Устойчивость пластин под действием продольных |
||
периодических сил. |
Прикладная математика и механика. Т. 2, |
№ 1, 1938. |
|
1 Б о л о т и н |
В. В. |
О поперечных колебаниях стержней, |
вызываемых |
периодическими продольными силами. Сб. «Поперечные колебания и крити ческие скорости». Вып. I. Изд-во АН СССР, 1951.
5. Б о л о т н н |
В. В. |
О |
параметрическом возбуждении поперечных |
колебаний. Сб. «Поперечные |
колебания и критические скорости». Вып. 2. |
||
Изд-во АН СССР. 1953. |
|
|
|
6. Б о л о т и н |
В. В. О взаимодействии вынужденных и параметрически |
||
возбуждаемых колебаний. «Изв. |
АН СССР*. ОТН, № 4, 1954. |
7.Б о л о т и н В. В. Динамическая устойчивость упругих систем. Гостехиздат, 1956.
8.Б о л о т и н В. В. Неконсервативные задачи теории упругой у с т о й чивости. Физматгиз, 1961.
9. |
В о л ь м и р А. С. |
Устойчивость упругих систем. Физматгиз, |
1963. |
|||||
10. |
Г о л ь де |
н б л а т |
И. |
И. |
Динамическая устойчивость сооружений. |
|||
Стройнздат, 1948. |
и б л а т |
И. И.. |
С и з о в |
А. М. |
Справочник по расчету |
|||
11. |
Г о л ь д е |
|||||||
строительных конструкций |
на |
устойчивость и |
колебания. Стройнздат, |
1952. |
||||
12. |
Д ж а н е л и д з е |
Г. |
Ю., |
Р а д ц и г |
М. А. |
Динамическая устой |
||
чивость |
кольца под действием |
нормальных периодических сил. П р и к л а д н а я |
||||||
математика и механика. Т. |
4. Вып. 5—6. 1940. |
|
|
|
13.Д ж а и е л и д з с Г. 10. Устойчивость упругих систем при динами
ческих нагрузках. Сб. «Проблемы устойчивости в строительной механике*. Стройнздат, 1965.
14.К р е й н М. Г. Основные положения теории Л — зон устойчивости канонической системы линейных дифференциальных уравнений с периоди
ческими коэффициентами. Сб. памяти А. А. Андронова. Изд-во АН СССР, |
1955. |
||||||
|
15. М а к - Л а х л а и |
Н. В. Теория и приложения функций |
Матье. |
||||
Изд-во ИЛ, 1953. |
|
Г. |
В. О динамической устойчивости пологих |
упру |
|||
гих |
16. М и ш е н к о в |
||||||
оболочек. «Инженерный |
журнал». 1961. Т. I. Вып. 2. |
|
|
||||
|
17. Ст р е тт М. Д. |
Функции Ляме, Матье и родственные им в физике |
|||||
и технике. ДНТВУ, |
1935. |
Динамическая устойчивость авиационных |
кон |
||||
|
18. Ч е л о м е й |
В. И. |
|||||
струкций. Изд-по «Аэрофлот», 1939. |
устойчивости |
упру |
|||||
гих |
19. Ч е л о м е й |
В. Н. |
О возможности повышения |
||||
систем при помощи |
вибраций. Доклады АН СССР. |
Т. ПО, |
3, |
1956, |
Литература |
369 |
20. Я к у б о в и ч В. А. Замечание к некоторым работам |
по системам |
линейных дифференциальных уравнений с периодическими коэффициентами. Прикладная математика и механика. Т. 21. Вып. 5, 1957.
21. Я к у б о в и ч В. А. Метод малого параметра для канонических систем с периодическими коэффициентами. Прикладная математика и меха ника. Т. 23. Вып. I, 1959.
22.Е V а п • I V а п о V з к I К. М. Оп 1ке рагатсЫс гезропзс о! 5(гис1игез. Арр1. МесЬ. КеУ., уо1. 18, п. 9. 1965.
23.М е Ш е г Е. А11дсте1пе ТНеопо с!ег 8(аЫШа1 еггугилдепег 5сН\у1п-
дипдеп ектзИзсЬег Кбгрег. 1пдеп1еиг—А г с Ы у , Вб. 17, п. 6, 1949.
24.М е I 11 е г Е. Оупапйс ЬискИпд. НапбЬоок о! Елд1пеегшд МесЬа- п!сз (еб. Ьу \у. Р1йдде), Мс Ога\у-Н111, №\у-Уогк, 1962.
25.Р 1з г с г е к К. Рагатс1пс сотЫлаНоп гсзопапсе (о! 1Нс зесолб
к1пб) 1п попПпеаг зуз1ешз. Рогргалуу 1пгк., 1. 8, 1960.
26. |
5 с к ш 1<11 |
О. ЙЬег сПе |
В1еде$сЬумпдипдеп без де1епк!д де1адег1еп |
|||||
ах1а! ри1з1сгелс1 |
Ье1аз1е1еп |
31аЬс&. |
Ма1ет. №скпсН1еп. Вб. 23, п. 2, 1961. |
|||||
27. |
3 с Н т |
1 б 1 |
О. ПЬсг б1е |
(2исгзс1г\у1пдипдеп вск^асН уогдскг0тт1ег |
||||
51аЬе Ье1 ехгсп1п5ск |
апдгеИелбег |
ри|51сгспбег |
ЬапдзЬе1а$1ипд. Ма1еш. Иаск- |
|||||
г1сЫел., Вб. 23, п. 4—5, 1961. |
|
е г |
Р. |
1пз1аЬШШеп дебатрГ- |
||||
28. |
5 с Нт |
1 б 1 |
О., |
V/ е I б е п к я т т |
||||
1ег гкеоПпеагег 5сЬ\ушдипдеп. Ма111ет. №скг1сЫеп, |
Вб. 23, N 4—5, 1961. |
|||||||
29. |
Ш е I б е п На т т |
е г Р. |
ЭДсМПлеаге В1едезс1т1пдипдеп без ах1а1— |
|||||
риЫегепб Ье1аз1е1сп |
51аЬез. 1пдсп1еиг—АгсЫу, |
Вб. 20, N 5, 1952. |
Глава 7
КОЛЕБАНИЯ ПЛАСТИНОК
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ КОЛЕБАНИЙ ПЛАСТИНОК
Уравнение изгиба пластинки в ортогональной системе координат для динамического случая. Пусть срединная поверхность изотропной пластинки постоянной толщины Н отнесена к ортогональной криво линейной системе координат хгх2 (рис. 1). Дифференциальное уравнение малых поперечных колеба ний в рамках гипотез Кирхго
фа-Лява будет
ОДДш + рЛ -щ т = <?; (1)
здесь ш—прогиб пластинки; 0 =
ЕН.3
— 12 (1 _ ^2) ~ Цилиндрическая
жесткость; р — плотность мате Рнс. 1 риала; ц (х1( х2, 0 — интенсив ность поперечной нагрузки; Е —
модуль упругости; V — коэффициент Пуассона;
Н гН Л дхЛ Н ! дх,)
обобщенный оператор Лапласа; Яц Я 2 — параметры Ламе.
Моменты Мп, М22. М12 и |
перерезывающие силы |
выража |
||
ются через прогиб |
пластины |
ш следующими |
формулами: |
|
ч[ |
1 а |
дш |
дНх |
дьи |