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

книги из ГПНТБ / Митропольский Ю.А. Интегральные многообразия в нелинейной механике

.pdf
Скачиваний:
16
Добавлен:
27.10.2023
Размер:
22.84 Mб
Скачать

20

ГЛ. I.

В В Е Д Е Н И Е

 

 

можно ввести нормальные координаты гь г%.......ZN , для ко­

торых

 

 

 

 

 

Z É=1

1/ - 4 -

Е юМ

(1.18)

 

Z

A=1

 

f \

®k — соответственно нормальные функции и собствен­

ные частоты, соответствующие квадратичным формам (1.16)), и уравнения Лагранжа для невозмущенного движе­

ния примут следующий

вид:

 

 

=

0

(ft = 1, 2, . . . , N).

(1.19)

Допустим теперь, что на рассматриваемую колебатель­

ную систему действуют малые возмущающие силы

 

zQk = е {Qi0) (q{, qc) + 2

[Qw? (<7,-, ?;) cos QJ -f

 

a

 

 

 

+ Q*2)(</i, ^)sinQ efl}

( t = l , 2 ..........N),

(1.20)

где Qa — частоты возмущающих сил, e — малый положи­ тельный параметр. Тогда, переходя в выражениях (1.20) также к нормальным координатам, получаем следующую систему нелинейных дифференциальных уравнений:

-~pL' +

- ?zdU zh *i) (ft = 1, 2..........

N), (1.21)

где Zk определяются из условия эквивалентности работ согласно формулам

N

 

 

( 1.22)

2 * = = S Q /ФІЛ

(ft =

1,2

/=і

 

 

 

Уравнения (1.21) путем замены переменных

+

x-bfi-1"**,

J

Д =

mkx_ke~1^

(1.23)

j

легко приводятся к уравнениям в стандартной форме

-~Ж~ — EX k(t, Xi) ( f t = ± l ,

± 2 , . . . ,

± N), (1.24)

в которых

—шJ

 

 

 

 

 

 

 

 

(1.25)

и, кроме того, —со-А = со*,

= Zk.

§ 1. И Н Т Е Г Р А Л Ь Н Ы Е М Н О Г О О Б Р А З И Я

21

В дальнейшем для упрощения систему дифференциаль­ ных уравнений (1.24) целесообразно записывать в вектор­

ной форме

 

- § - = е * (/,*),

(1.26)

где X , X — n-векторы (в нашем случае п — 2N). Правые части системы (1.26) являются периодическими функциями или почти периодическими в зависимости от того, как зави­ сели от t функции eQk.

К уравнениям типа (1.26) могут быть приведены также задачи о воздействии на колебательные системы внешних сил высокой частоты. Так, например, пусть движение та­ кой системы характеризуется уравнениями вида

Ц

Чь 4 г )

(s =

1, 2,

... . N), (1.27)

где Fs (wf, qlt

— периодические

(или

почти периоди­

ческие) функции t и о — «большой» параметр.

Тогда, вводя новые неизвестные согласно формулам

* =

qs = xs,

-jf- = y s (s=

1, 2,

. . . , N), (1.28)

вместо системы (1.27) получаем систему уравнений

" $ Г

= еу*’

= s f ^ x’ Хі’ Уі>

(s =

Ь 2...........N)’

 

1

 

 

(1.29)

где e =

параметр.

 

 

—---- малый

 

 

К уравнениям вида (1.26) может быть приведен и более общий случай, чем указанный выше, когда состояние дина­ мической системы характеризуется угловой переменной а и п переменными xlf х2.......хп и описывается следующей

системой

уравнений:

 

 

 

 

dx

xs(об, х^у

. • • t хп)

(s

1> 2,

. . . j

dt

da

 

 

 

 

(1.30)

= Ы (xlf . . .

, хп) +

А (а,

х1г

, хп),

dt

где X — большой параметр, X s (а,

хіу

хп), А (а, хи ...

...,хп) — периодические (или почти-периодические) функции

22

ГЛ. I. В В Е Д Е Н И Е

угловой

переменной а с периодом 2л. (Системы типа (1.30)

встречаются при изучении гироскопических систем, в тео­ рии ускорителей и т. д.)

Вводя новую независимую переменную т = Kt, можем

представить систему (1.30) в виде

 

 

 

 

 

 

 

 

-fö- =

eXs{a,

xlt

, хп)

(s =

1, 2......... п),

 

(1.31)

~

=

со(х1,

... , хп) + еЛ(а,

 

. . .

,

хп),

 

 

Ху,

 

 

 

где е =

1---- малый

параметр,

или,

исключая

т,

в

виде

dx s

=

eX*s (а,

Ху,

. . . , хп,

е)

(s —

1,

2,

...

, п),

da.

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.32)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X s (a,

Ху,

. . .

,

х п)

 

 

 

К (а, Ху,

. . . , хп, е)

 

 

 

 

Ö) (хх,

. . . .

х п) + 8Д (а,

Ху,

. . . ,

х п) '

 

 

 

 

 

К уравнениям вида (1.26) могут быть сведены также уравнения более общего вида, чем (1.21).

Как известно, нестационарные процессы в нелинейных колебательных системах во многих случаях приводятся к рассмотрению систем дифференциальных уравнений типа

-4 " { 2 ац 00 <7;} +

2

Ьц (т) <7 ( =

 

 

= eQj (т, Ѳ, qy, ...

, qN, qt, ... , qj4,

e)

( / = 1 , 2 , . . . , M),

 

 

 

 

(1.33)

где qt — обобщенные

координаты,

eQj

— внешние возму­

щающие силы, т — et — медленное

время (медленное по

сравнению с естественной единицей времени — порядка пе­ риода собственных колебаний), е — малый положительный параметр, а;/ (т) = щі (т), Ьц (т) = Ьц (т) — инерционный и квазиупругий коэффициенты, медленно изменяющиеся со

временем,

dQ

,

,

 

 

 

ѵ (т).

 

 

Система типа (1.33) с помощью замены

переменных со­

гласно формулам

 

 

 

 

 

 

N

 

 

 

 

<7і =

2

Ф

(t' = 1, 2, ...

, N),

(1.34)

k—\

§ ! И Н Т Е Г Р А Л Ь Н Ы Е М Н О Г О О Б Р А З И Я

23

где tp£fe>(т) — нормальные функции, зависящие от парамет­ ра т и удовлетворяющие условию ортогональности, может быть приведена к следующему виду:

-ZjjT +

Zk(т’ ѳ- 2i> • ■■» г‘Ѵ, Zj......... zN, e)

 

( * = 1 , 2 ......... N),

(1.35)

после чего, произведя замену согласно формулам

 

г* = Хке1Ы

^ ‘ +

х . ке-^

(1.36)

г, = «о*(т) Ѵ Ь

(Т,Л -

 

приходим к системе уравнений вида

 

 

 

= еХ (т, Ѳ, X, е)

(1.37)

|- ^ - = ѵ(т), т = eij , в которой правые части зависят от

медленного времени х и являются периодическими, почтипериодическими или квазипериодическими по Ѳ в зависи­ мости от того, как зависят от Ѳфункции Q,-.

Заметим, что к системе уравнений вида (1.37) может быть приведена и более общая система дифференциальных урав­ нений, описывающая колебания системы со многими степе­ нями свободы, характеризующиеся в невозмущенном со­ стоянии (т. е. при е = 0, х — const) лагранжевой функ­

цией

j Г N N N 1

£ = -Ö- | 2 ■аіі(х)йіЯі + 2

g iiW q â — 2

bij(x)qiqi\,

1

U ',/= i

> ./= 1

i , j = i

)

где gij

(г) =

gji (T).

 

(1.38)

 

 

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

В связи с этим мы будем рассматривать также системы, близкие к автономным, вида

~ = Х(х) +

X, е),

(1 .39)

24

ГЛ. 1. В В Е Д Е Н И Е

 

и системы, близкие к неавтономным:

 

 

= X (/, х ) + еХ* (yt, X , е).

(1.40)

К системам типа (1.39) и (1.40) при ряде предположений приводятся многочисленные задачи теории колебаний (од­ ночастотные колебания систем со многими степенями сво­ боды, не близкие к гармоническим, колебания релаксацион­ ного типа и т. д.).

Будем также рассматривать и более общие системы диф­

ференциальных уравнений вида

 

 

=

X (т, х) -f еХ* (т, 0, X,

е),

(1.41)

=

X (т, Ѳ, х ) + еХ* (т, 0,

X, г ).

(1.42)

Многие задачи приводят к рассмотрению колебательных систем, описываемых дифференциальными уравнениями, в коіорых одни переменные изменяются быстро, а другие медленно. К уравнениям такого типа относятся, например, указанные выше уравнения (1.30).

В общем виде такие уравнения можно записать следую­ щим образом:

~ ~ = X (х , у) +

(t, X, у),

(1.43)

=(*’ х’ У)-

где X , у — соответственно п- и m-векторы, е — малый по­ ложительный параметр.

Проблемы, приводящие к рассмотрению уравнений типа (1.43), встречаются в динамике спутников и искусственных небесных тел. Частным случаем системы (1.43), когда пе­ ременная X отсутствует, являются уравнения в стандартной форме. При Y (t, X , у) £= 1 приходим к рассмотрению урав­ нений типа (1.41).

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

§ 1. И Н Т Е Г Р А Л Ь Н Ы Е М Н О Г О О Б Р А З И Я

25

ных нерегулярно-возмущенных дифференциальных уравне­ ний вида

= fix, г, i), е

= F (х, z, t),

(1.44)

где е — малый параметр, х, / — «-векторы, г, F — т-век- торы, а также системы линейных нерегулярно-возмущенных дифференциальных уравнений

■^r = A(t)x + B(t)y + h(t),

(1.45)

е= С (t) у + eF (t) X + еН (t),

где Ху h — m-векторы, у, Н — п-векторы, А, В, С, F — со­ ответственно (m X т)-,X п)-, (п X «)-, (п X т)-мат- рицы.

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

 

j £ - = ALt)x(t) + A1( f ) x ( t - A )

+

 

 

+ X(t,g(t),

x(t),

X (t — A), у (t),

y(t — Д), e),

 

(ef) {/(() +

(eO y(t — b) +

(1.46)

 

+ eY (t, g (t),

X (t),

X (t — A), у (t), у (t A), e),

 

= © (t) +

G (t,

g (t),

X (t),

у (t), e),

 

где X ,

у, g — соответственно n - , m - ,

/г-векторы, А

(/), A t (/),

В (et),

Bi (et) — ограниченные

квадратные

матрицы,

© (t) — fe-вектор-функция, е — малый положительный па­ раметр.

Будут рассмотрены также нерегулярно-возмущенные уравнения с запаздыванием

~= f(t, х,у, е),

8 ЧГ = А V' х>Ха) у + В У>х>Ха) У* + § ((>х’ А'Д' У' Уь>е)’

(1.47)

26

;

 

ГЛ. I. В В Е Д Е Н И Е

где

х т-вектор,

у — п-вектор, хд — х (t — еД), //д =

=

у

(t — еД), е >

0 — малый параметр.

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

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

§2. Вспомогательные сведения из линейной алгебры

ианализа

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

1. Матричные обозначения. Рассмотрим квадратную матрицу А порядка п с элементами аи- (і, / = 1, .... п), ко­ торые, вообще говоря, являются комплексными. Матрица (кі А), где к — независимое переменное, / — единичная матрица порядка п, называется характеристической мат­ рицей для А. Ее определитель

Ф (Я) ~ det (kl — А)

(2.1)

является многочленом относительно к и называется харак­ теристическим многочленом матрицы А. Корни характери­ стического многочлена матрицы называются ее характери­ стическими числами или собственными значениями. Обоз­

начив эти корни kt

= 1, ..., л),

запишем: det (kl А)

=

П

 

 

 

= п—к).

 

 

 

1

собственных

значений) матрицы

А

Спектр (спектр

есть множество всех ее собственных значений. Он состоит

из тех значений Я,, для которых матрица (kl А)

необрати­

ма. Те значения к, для которых матрица (kJ А)

обратима,

называются регулярными, а матрица (kJ — Л)-1 называется

§ 2. С В Е Д Е Н И Я ИЗ Л И Н Е Й Н О Й А Л Г Е Б Р Ы И А Н А Л И З А

27

резольвентой матрицы Л. Спектр п собственных значений матрицы А порядка п совпадает с множеством корней ал­ гебраического уравнения

det (Я/ — А) =-- 0

(2.2)

(характеристического уравнения матрицы А), или в развер­

нутом виде:

а,

 

 

'X ап

Мп

 

 

М2

 

det

 

— Я,

(2.3)

 

= 0.

С1цЛ

Сіп

 

 

Кратность каждого корня Я/ этого уравнения равна его

алгебраической кратности т,- как собственного значения,

так что т\ + т-2 + ... = п.

Для любой квадратной матрицы А с собственными зна­ чениями Я/ матрица аА имеет собственные значения аЯ;;

матрица А р — собственные значения Я/ (р — 0, 1, 2,...); каждый многочлен или аналитическая функция / (Л) име­ ет собственные значения / (Я/).

Если две квадратные

матрицы Л

и ß имеют

собствен­

ные значения Я/ и

то

спектром

прямого произведения

Л 0

В является

множество всевозможных

произведе­

ний

Я/рй.

 

 

 

 

Имеет место следующий фундаментальный результат о

канонической форме матрицы.

 

 

Каждая квадратная матрица Л порядка п с помощью

преобразования подобия

 

 

 

 

 

J = 5Л5~'

 

 

(det S ф 0) может быть приведена к нормальной жордановой форме

У = diag [Ух (ЯА), . . . , Jrn(Ят)] < п),

(2.4)

где

 

К

1

0

•• •

0

 

о

яр 1

.. •

0

 

(Яр)

 

 

 

(р)

(2.5)

 

 

 

 

0

0

0

1

 

0

0

0

.. •

я„

 

28

ГЛ.

I. В В Е Д Е Н И Е

(р ■= 1,

m; li > 1,

lm > 1), N\p) — нильпотентная

матрица, причем каждому характеристическому корню Хр

кратности а р соответствует

одна или

несколько клеток

Жордана размерами 1{р\

.... Ір],

таких,

что

4 *4"

• • •

+

Ір*—

 

Каждой клетке Жордана

J р (Хр) порядка /р с точностью

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

Степень вырождения г собственного значения Хр пред­ ставляет собой максимальное число линейно-независимых собственных векторов матрицы А , соответствующихХр. В об­ щем случае г С ар. Если степень вырождения г собствен­ ного значения равна его кратности, т. е. г = ар, то, оче­

видно, /р*= ... = t'p = 1. В этом случае все соответствую­ щие клетки Жордана будут содержать по одному элементу (простые клетки).

Пусть Хі,..., Хт— характеристические числа матрицы Л, соответствующие различным клеткам Жордана (не обяза­ тельно различные между собой). Тогда, так как XI —

= S~lXS, А — S~lJS,

то

характеристический полином

Д (X) матрицы Л можем записать в виде

 

Д (X) = det (XI — J) = (X XJ1* . . . (X — X jm

 

(h +

’ ' •

+

Ал “ n)‘

(2.6)

Множители (X XP)1P (p =

1, .... m) называются

эле­

ментарными делителями матрицы Л, а натуральные числа Ір, т. е. размеры (порядки) клеток Жордана,— показателя­ ми элементарных делителей, соответствующих числу Хр.

Если все характеристические числа Хр имеют простые

элементарные делители (Ір — 1), то матрица Жордана J бу­

дет чисто диагонального вида

Хг

О

J =

К

= diag (Xv . . . , Ха), (2.7)

ОХп

причем числа Хр (р = 1.......п) не обязательно различны.

§ 2. С В Е Д Е Н И Я ИЗ Л И Н Е Й Н О Й А Л Г Е Б Р Ы И А Н А Л И З А

29

2. Линейные преобразования. Каждой квадратной мат­ рице А порядка я соответствует линейное преобразование А

векторного координатного пространства Сп размерности п,

т. е. вектору х = (хи ..., хп) пространства Сп ставится в соответствие вектор

А х = у = ( У ! ..........

У п ) ,

(2.8)

определяемый соотношением

п

У і = 2 a ikx k.

k=i

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

ствующего в Сп.

Говорят, что подпространство Сп линейного простран­ ства Сп инвариантно относительно преобразования А , если каждый вектор из Сп преобразование А снова переводит в некоторый вектор из Сп , т. е. А С п сг Сп .

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

Характеристическим многочленом преобразования А на­ зывается характеристический многочлен <р (Я) = det (XI

А) матрицы А.

Степень характеристического многочлена преобразова­ ния Л равна порядку матрицы А, а порядок матрицы равен

размерности пространства Сп. Поэтому степень характерис­ тического многочлена преобразования А равна размернос­ ти пространства, в котором действует это преобразование.

Так же, как и квадратная матрица, каждое линейное преобразование является корнем своего характеристиче­ ского многочлена, т. е. q> (А ) = 0.

Число % называется собственным значением линейного преобразования А , если в рассматриваемом пространстве

Соседние файлы в папке книги из ГПНТБ