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

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

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

140 г л . III. М Н О Г О О Б Р А З И Я В Б Л И З И П Е Р И О Д И Ч Е С К И Х Р Е Ш Е Н И Й

где hx (т, Ѳ, ф, а),..., hm (т, Ѳ, ф, а) сопряжены с

(х, Ѳ,

ф, а),

..., hm (х, Ѳ, ф, а) и могут

быть определены форму­

лами

типа

 

F^m(x, а)е^п'в+т^

 

 

(т, Ѳ, ф, а)

 

 

 

I І і

[пѵ (т) +

т о (т, а)] — Н (т, а )

 

 

 

 

3. Доказательство существования точного двупарамет­ рического семейства решений уравнения (5.1). Согласно формуле (1.65), параметрическое представление точного двупараметрического семейства решений уравнения (5.1) имеет вид

X = х° (т, ф (т, t, г), а (т, t, е)) -f

+ {Ѳ (т, ф (т, t, е), а (т, t, в)) f (т, 0, ф (т, і, г), а (т, t, г), е) +

+ Ѳ (т, ф (т, t, г), а (т, t, е)) / (т, Ѳ, ф (т, t, е), а (т, t, е), е)}, (5.16)

где ф (т, і, е), а (т, і, е), /г = / (т, Ѳ,ф, (т, t,e), а (х, ва­ точные решения системы (5.3). Таким образом, для дока­ зательства существования семейства решений (5.16) нужно доказать существование точных решений ф, a, h системы (5.3).

Это доказывается по аналогии с доказательством су­ ществования решений фл а(, ht уравнений (2.6).

Аналогом интегрального уравнения (2.31) здесь будет уравнение

F (t, ф, а, е) =

с о

=j G(z)R1{t + z-, ф2; а2; F (t + г\ ф2; аг\ е); е} dz. (5.17)

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

необходимости сделать следующие

замечения.

а) Так как функции F (т, Ѳ, ф, а,

в

) , R (х, Ѳ, ф, а, /г, е)

зависят от т, Ѳ, ф, а, г, где т = et, Ѳ =

^ ѵ (т) dx = Ѳ (t, e),

то можем написать

 

о

 

 

F (et, Ѳ(t, e), ф, а, e) = Fl (t, ф, а, e),

R (et, Ѳ(t, e), ф, a, h,e) =■

 

 

= Ri(t, Ф, a, h, e)

$ Б. У Р А В Н Е Н И Я С М Е Д Л Е Н Н О М ЕН Я Ю Щ И М И С Я П АР А М Е Т Р А МИ 141

и полагать, что Fx, Rx определены для ££ [0,

£/е];анало-

гично

 

(t, ф, a, h, e),

 

Р (et, Ѳ(t, е), ф, а, h, e) =

 

 

S (et, Ѳ(t, e), ф, a, h, s) =

 

(t, ф, a, h, e).

 

 

б) В преобразовании (5.17) t

изменяется в

интервале

(О, L/e 1, а z — в интервале (—о о

,

о Функциио ) .

Рг и Rlt

стоящие в правой части преобразования, являются функ­ циями t + z. Следовательно, чтобы к функциям Flt Rx мож­ но было применить преобразование (5.17), их надо доопре­ делить на интервалах (—оо, оо). Обозначим эти функции, определенные на интервале (—оо, оо), через F2, R2. При этом полагаем, что в расширенной области функции F2, R2 об­ ладают теми же свойствами, что и функции Flt R± в исход­ ной, а в исходной области совпадают с ними.

в) При изменении t на отрезке [0, L/e]а (t) изменяется в области 21 = [а0, аД.

Являясь функцией г, аг при изменении z от—оо доф-оо может выйти из области своего определения 21. Поэтому для аг будем рассматривать расширенную область, которую

обозначим через 2t.

Таким образом, функции F2, R2 будем рассматривать

в области

 

t £ R , а£Ш, ф £ П, h £ U a, 0 < е < е о

(5.18)

и обозначать их через Fs, R3. При этом полагаем, что в области (5.18) они обладают теми же свойствами, что и

функции F (т, Ѳ, ф, а, е), R (т, Ѳ, ф, а,

г) в области (5.6),

а в области (5.6) совпадают с ними.

ф, а, h, в) также

Для функций Рг (t, ф, a, h, в), Sx (t,

будем рассматривать

расширенную область (5.18), так как

аг, ф2 определяются

из системы

 

-^ r = (t, ф, а, h, е),

=со (et, а) + Рг (t, ф, а, h, в),

иобозначать их в этой области через Р3, S 3, со (т, а) рас­ сматриваем также в расширенной области и обозначаем ее

вэтой области через ац (et, а), причем полагаем, что в расширенной области coj (et, а) обладает теми же свойства­ ми, что и со (т, а) в исходной области, а в исходной — сов­ падает с со (т, а).

142 ГЛ. III. М Н О Г О О Б Р А З И Я В Б Л И З И П Е Р И О Д И Ч Е С К И Х Р Е Ш Е Н И Й

Матрица Н (et, а) у нас определена для t £ Ю, L/e]. Вместо матрицы Н (et, а) вводим в рассмотрение матрицу Нг (t, а), определенную на интервале (—оо, оо), при этом полагаем, что на интервале (—оо, оо) она обладает тем

свойством, что для решения уравнения

= H1 (t, a)U

справедлива оценка

 

 

\U (t) | < Ке~У\‘\

 

а на промежутке

[0, L/e] совпадает с Я (т,

о).

Таким образом, вместо системы уравнений (5.3) pac-

сматривается следующая система:

 

"2 Р =

(et, а) + Р3 (t, ф, а, h, е),

 

=

S3 (t, ф, а, h, е),

(5.19)

=

H1 (t, a)h + R3 (t, if, а, h, e),

 

которая в области (5.6) эквивалентна системе (5.3). Доказав существование и единственность решений

г|) (т, t, е), а (т, t, е), h = f (т, Ѳ, if (т, t, е), а (т, t, е), е)

уравнений (5.3) и подставив их значения в формулу (1.65), получим следующее представление точного двупараметри­ ческого семейства решений уравнения (5.1):

X = х° (т, ф (т, t, е), а (т, /, е)) +

+ -у- {Ѳ(т, ф(т, t, е), а (т, t, е)) f (т, Ѳ, г|з (т, /, е), а(т, t,e), e)-f

+ Ѳ (т, ф

(т, t, е), а (т, t, e))J (т, Ѳ, ф

(т, t, е), а (т, t, е), е)} =

 

 

=

Фт(т, Ѳ, ф (т, t, е), а (т, і, е), е), (5.20)

где ф (т,

t,

е), а (т, t,

е) определяются из системы

^1'

 

© ( т , а) +

Р (т, Ѳ, ф, а, /

(т, Ѳ, ф, а, е), е),

— L

=

 

 

 

 

(5.21)

~= S (T , Ѳ, ф, а, / (т, Ѳ, ф, а, е), е)

как функции времени, зависящие от двух произвольных постоянных.

§ 5. У Р А В Н Е Н И Я С М Е Д Л Е Н Н О М Е Н Я Ю Щ ИМ ИС Я ПАР А МЕ ТР А МИ 143

4. Оценка разности между точным семейством решений и его т-м приближением. Семейство решений (5.20) пред­ ставляет собой совокупность интегральных кривых, за­

полняющих некоторую поверхность в пространстве R'1, которую обозначим через ST.

Параметрическое представление этой поверхности ана­ литически может быть задано либо в виде (5.20), где пара­ метрами являются произвольные постоянные, от которых

зависят ф (т, t,

е), а (т, t, е), либо в виде

X = X (т, % а) +

 

+ -гг {Ѳ (т, 1 ]), а) f (т, Ѳ, ф, а, е) + Ѳ (т, я|з, а) / (т, Ѳ, -ф, а, е)} =

= Фт(т, Ѳ, ф, а, е), (5.22)

где ф, а — произвольные параметры, изменяющиеся в той же области: ф £ Q, а £ Ш.

Из существования представления (5.20) с очевидностью вытекает существование представления (5.22).

Аналогичные рассуждения имеют место и относительно приближенного двупараметрического семейства решений (5.15) . Это семейство решений также целиком заполняет

некоторую поверхность в пространстве Rn, которую обо­ значим через Snp. Аналитически эту поверхность можно представить либо в виде (5.15), либо в виде

X = х° (т, ф, а) +

+ -гг {Ѳ (т, ф, а) hm(т, Ѳ, ф, а, е) + Ѳ (т, ф, а) h'n(т, Ѳ, ф, а, е)} = = Фпр (т, Ѳ, ф, а, г), (5.23)

где ф, а — произвольные параметры.

Из такой идентичности представлений (5.20) и (5.22), (5.15) и (5.23) вытекает следующий факт.

Если требуется получить оценку между точным двупа­ раметрическим семейством решений (5.20) и его m-м при­ ближением (5.15), то в аналитических выкладках вместо выражений (5.20), (5.15), представляющих собой функции многообразий (соответственно точного и приближенного), можем пользоваться выражениями (5.22) и (5.23), которые являются параметрическими представлениями соответ­ ствующих многообразий, так как, как кривые, определяе­ мые функциями (5.20) и (5.15), так и кривые, определяемые

144 ГЛ. Ш . М Н О Г О О Б Р А З И Я В Б Л И З И П Е Р И О Д И Ч Е С К И Х Р Е Ш Е Н И Й

функциями (5.22) и (5.23), лежат на одних и тех же поверх­ ностях и полностью их заполняют.

Таким образом, задача сводится к оценке близости поверхностей 5Т и 5пр, параметрическое представление которых имеет вид (5.22) и соответственно (5.23). Имеем

I Фт (т, Ѳ, ф, а, в) — Фпр (т, Ѳ, ф, а, г) \ = \х° (т, ф, а) +

+ -гг {Ѳ (ч 'Ф. а) / (т, Ѳ, ф, а, е) + Ѳ (т, ф, а)} (т, Ѳ, ф, а, е)} —

х° (т, ф, а) (Ѳ (т, ф, а) hm(т, Ѳ, ф, о, е) +

+ Ѳ (т, ф, а) hm(т, Ѳ, ф, а, е)} <

< -J- ( I Ѳ (т, ф, а) 11 / (т, Ѳ, ф, а, е) — hm (т, Ѳ, ф, а,е)| +

+ IѲ (т, ф, а) 11 / (т, 0, ф, а, е) — hm(т, 0, ф, а, г) |)

или, принимая во внимание ограниченность матриц Ѳ, Ѳ некоторой константой М\

I Фт (т, 0, ф, а, е) — Фпр (т, 0, ф, а, г) | <

< М ( I / (т, Ѳ, ф, с, е) — hm (т, Ѳ, ф, о, е) | +

+

I / (т, 0, ф, а, е) hm (т, Ѳ, ф, а, е ) | ). (5.24)

Таким образом,

требуется получить оценку

для

I / (т, 0, ф, а, е) — hm(т, 0, ф, а, е) |

(5.25)

(оценка для |/ (т, Ѳ, ф, а, е) — hm (т, 0, ф, а, в)| полу­ чается аналогично).

Чтобы оценить выражение (5.25), необходимо вместо дифференциального уравнения относительно h рассматри­

вать

интегральное

 

 

 

 

СО

 

 

 

 

SF3= j U (z) R3 (t + г; а; ф;

F3(t + г; а;

ф; е); е) dz. (5.26)

Рассмотрим

разность

 

 

 

/ (t, ф, а, е) — hm (t, ф, а, е) | =

 

 

 

 

)

U (г) [R3 (t + z;

ф, а, f (t, ф, а, е), е) —

 

 

R3 {t +

г; ф; а; hm(t, ф, а, е), е) +

em+Vm] dz

, (5.27)

где

em+i гт означает совокупность членов порядка

ет +‘.

§5. У Р А В Н Е Н И Я С М Е Д Л Е Н Н О М Е Н Я Ю Щ ИМ ИС Я ПАР А МЕ ТР А МИ 145

Принимая во внимание, что функции R 3 удовлетворяют условию Липшица по ф, а, h с функцией X (е, р) в качестве постоянной Липшица, а также ограниченность матрицы U (z) (I U (z) I С /Се~ѴІг|). можем написать

I / (t, ф, а, г) — hm (t, ф, а, е) | <

< - у

(е, р) I f (t, ф, а, е) — hm(/, ф, а, е) |] + ет+І | гт|. (5.28)

Выберем теперь такое постоянное р и такое s', чтобы для

всех р <

_

е' выполнялось неравенство

TS

 

1.

Р и е <

X(в, р) <

Такой выбор е', р всегда можно сделать, так как X (е, р)

О

при е -> 0, р -> 0.

 

 

в виде

 

Тогда неравенство (5.28) можем записать

 

I / (t, ф, а, е) — hm (t, ф, а, е) |

1 — К

X (е, р)

рт+і I г

 

 

 

 

 

 

<>£

\'rt

 

откуда

следует

 

 

 

 

 

 

 

I / (t, ф, а, е) — /Г (t, ф, а, е) | < етГт+

, (5.29)

где через гт обозначено \гт\.

 

 

 

 

Аналогичное

неравенство

имеет

место

относительно

I / ((, ф, а, е) — hm (t, ф, а, е) |.

Возвращаясь в системе (5.19) к исходному интервалу времени, мы тем самым совершим переход к исходному интервалу времени и в оценке (5.29) и, таким образом, будем иметь следующую оценку:

I / (т, Ѳ, ф, а, е) — hm(т, Ѳ, ф, а, е) | < sm+1 гт+ Аналогично

I / (т, Ѳ, ф, а, е) — hm(т, Ѳ, ф, а, е) | < em+1 гт+ • —

Возвращаясь к неравенству (5.24), можем написать

I Фт (т, Ѳ, ф, а, е) — Фпр (т, Ѳ, ф, а, е) | < 2Мет+'гт

=

- Lem+1. (5.30)

146 г л . III. М Н О Г О О Б Р А З И Я В Б Л И З И П Е Р И О Д И Ч Е С К И Х Р Е Ш Е Н И Я

§ 6. Системы уравнений, описывающие «быстрые»

и«медленные» движения

1.Основные предположения. Рассмотрим систему п +

+т уравнений с «быстро» и «медленно» изменяющимися

переменными следующего вида:

-ld xr = X (y)x + еХ* (t, X, у),

(6.1)

■7Г = еѴ V’ У)>

где X, X* и у, Y — Соответственно п - и m-векторы, X ) "

— (я X п)-матрица, е — малый положительный

параметр.

Предположим, что вектор-функция Y (/, х, у) допускает

существование среднего

по

і:

 

 

 

т

 

Y, {х, у) =

lim

) Y (it, л:, у) dt,

(6.2)

7->oo

1 0

 

и наряду с системой (6.1) рассмотрим вспомогательную систему

- ^ = Х(У)х,

(6.3)

■^- = вУ0(х, у).

Пусть выполняются следующие предположения.

1°. Система уравнений (6.3) допускает существование

двупараметрического семейства

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

решений

X = 0, у = у0 (at + ф, a)

(cut + ф = ф),

(6.4)

где у° (ф, а) — периодическая функция ф, при этом а> в общем случае является функцией а и удовлетворяет усло­ вию Липшица

 

I со (a') — со (а") | с С | а' —• а" |.

(6.5)

 

Полагаем, что С является достаточно малой постоян­

ной порядка 8, так что со медленно изменяется

с изменени­

ем

а.

 

в

2°. Функции в правой части системы (6.1) непрерывны и

области

 

4 X Од X D0 X Е £о

( 6 . 6 )

§ 6. СИСТЕМЫ У Р А В Н Е Н И Й , О П И С Ы В А Ю Щ И Е Д В И Ж Е Н И Я 147

(Dд — A-окрестность х = 0;

D„ — ст-окрестность у =

= г/° (ф, о); L — конечное число) удовлетворяют условиям

X*(/,X, у) £Сіу, {X (у) X, У (t,X, у)} £С%у.

При этих предположениях докажем существование дву-

параметрического локального

интегрального многообразия

St системы (6.1) в окрестности семейства решений (6.4), установим его свойства и сведем рассмотрение исходной системы на многообразии к двум уравнениям относительно

переменных

tf и

а.

 

 

Заметим,

что

в аналогичной постановке Дж. Хейл

[214]

исследовал

частный случай

системы

(6.1), когда

X (у) = А,

где А — постоянная матрица.

 

2. Уравнения

специального вида. Представим исходную

систему уравнений в виде

 

 

= X (у0)х + Х г (t,X, у, г),

 

 

=

еУ0 (0, у0) +

еѴ'ох (0, у0) х + гУ'0у (0, у0) y + eY1 (t, х, у),

где

 

 

 

 

(6.7)

Х г (i, X, у, е) = X (у) X — X (у0) X + еХ* (t, х, у),

 

EY ! (t, X, у) =

eY (t, х,у) еУ0 (х , у) +

еУ0 (х, у)

еУ0 (0, у0)

 

 

 

гУ0х(0, if ) X eY0y(0, у0) у.

Совершим теперь в полученных уравнениях замену пере­ менных согласно формулам

X — X1,

У = У1+ еУ'олг (0, у0) X-1 (у0) х \

где по-прежнему принадлежит области DA, а относитель­ но у1 предположим, что она принадлежит D-, при этом

о — такое, что а + Л4Д < а, {Уо* (0, у0) X "1 (у0) | < М, и, следовательно, у не выходит из области своего опре­ деления D,J.

В результате такой замены система уравнений (6.1) примет вид

d x1

= X (у0) X 1 + Х2 (t, X 1 , у1, г),

dt

dyl

(6.9>

= еУ0 (°> У°) + еУ’оу(0, у0) у1+ еУ2 (/, х \ у1),

t dt

148 ГЛ. III. М Н О Г О О Б Р А З И Я В Б Л И З И П Е Р И О Д И Ч Е С К И Х Р Е Ш Е Н И Й

где

 

Х2 (i, х \ у1, е) =

Хх (t, X 1, у1 + гУОх (0, у0) Х~1(у0), х1, г),

е У 2 (t, X 1, у 1, е )

= г У ! (t, X1, у1+ і'У'ох(0, г/°) X-1 (г /°), х 1)

- еУ'о* (0, г/°) X“ 1(г/>) Х2 (/, д \ у \ е) +

4- Е2Уоу (0, у0) Уо* (О, у0) X-1 (у0) X1.

Рассмотрим соответствующие системе (6.9) уравнения в вариациях

dbx

X(y°)öx,

dt

(6 .10)

dby

sy'oy(0,y°) Sy.

dt

Для уравнений (6.10) по аналогии с § 1 введем в рассмот­ рение интегральные матрицы X (4, а), У (ф, а, е), матрицы монодромии X (Т , а), У (Т , а, г), а также ln X (Т , а),

In У(Т, а,Е), С (а) = ln X (Т , а), В (а, в) = In У (Т , а, е),

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

е< ех разложение

В{а, е) = В0 (а) + еВх (а) +

причем в рассматриваемом случае В0 (а) = 0.

Из предположения 1 ° вытекает существование двух ну­ левых точек спектра матрицы В (а, е). Обозначим их рх = = р2 “ 0 и предположим, что они являются изолирован* ными точками спектра. Полагаем, что остальной спектр <j0 (В) матрицы В (а, е), а также спектр а (С) матрицы С (а) не пересекаются с мнимой осью и расположены слева от нее.

Рассмотрим функции

Ф (ф, а) ~ X (ф, а) е~^с <а>, Ѳ (ф, а, е) = У (ф, а, е) е~^в (а'е>

(6.11)

и соответствующие им комплексно-сопряженные функции Ф (ф, а), Ѳ (ф, а, е), для которых справедливы следующие

§ 6. СИСТЕМЫ У Р А В Н Е Н И Й , О П И С Ы В А Ю Щ И Е Д В И Ж Е Н И Я 149

соотношения:

дФ ay а)■® (а) -f Ф (ф, а) со (а) С (а) = X (у0 (ф, а)) Ф (ф, а),

дФ (Ф, а) со (а) -f Ф (ф, а) со (а) С (а) = X (у0(ф, а)) Ф (ф, а),

д \ р

 

 

 

 

 

 

 

(6.12)

"ѳ (у

’ е)'~ю ^

+

ѳ СФ’ а’е) “ (а) в (а>е) =

 

-

 

 

_

=

еУо^(0, г/° (Ф, а)) Ѳ (ф, а, е), *

 

 

 

 

.

_

с (6.13)

5Ѳ (у

’в'• м (а) +

Ѳ

Ф е) ® И

5 (а- е) =

 

 

 

 

 

=

eF0;, (0, у0(ф, а)) Ѳ (ф, а, е).

 

Введем теперь

в

системе (6.9)

вместо х1 (х\.......

х'п),

У1 (У\..... У]) новые переменные g (gb

gn), ф, а, h (hu ..,

h j

посредством следующей

замены:

 

 

 

 

X1~

\

[ф (Ф, а) g +

Ф (ф, а) g],

 

 

 

 

 

 

 

 

(6.14)

 

if =

-гг [Ѳ (Ф, а, г) h + Ѳ (ф, а, е) К].

 

Подставляя выражения (6.14) в уравнения (6.9) и прини­ мая при этом во внимание соотношения (6.12) и (6.13), по­ лучаем

1

/ ЗФ

 

.

ЗФ

-

dip

со (а)

1

2

\ Зф

&

'

Зф

^

ЧГ

 

 

 

 

 

 

ЗФ

 

!

ЗФ

-

 

 

 

 

,

д а

£

Г

д а

£ ,

 

 

 

 

+

4"ф ( т

— НН“)г) =

, a , g , h , t ) ,

 

 

 

 

 

 

 

 

 

|(6.15)

 

 

 

+

 

 

 

<«>) +

 

+

X (-35- h +

! И

% +

4 - ѳ ( f

,Н г (а).к) +

+ 4 ~ ѳ ( " § ------

s H 2 (a)h) — e Y 3 (t, ф, a , g , h, е),

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