Материал: Прочность, устойчивость, колебания. Т. 3

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

538

С т а т и ст и ч еск и е за д а ч и к о л е б а н и й и уст о й ч и во ст и

процедуры статистической динамики. Если оператор Ь имеет вид, от­ личный от оператора (51), то решения, получаемые двумя методами, вообще говоря, не совпадают.

МЕТОДЫ ИССЛЕДОВАНИЯ НЕЛИНЕЙНЫХ И ПАРАМЕТРИЧЕСКИХ СЛУЧАЙНЫХ КОЛЕБАНИЙ

Метод .статистической линеаризации [24, 41]. Этот метод нашел ши­ рокое применение в задачах автоматического управления. В задачах колебаний механических систем наиболее распространен следующий вариант метода. Рассмотрим для простоты случай системы с одной степенью свободы. Уравнение

х + Р* + юЦх + рУ (х, х) = / (0>

(53)

где р. — малый параметр, /(/) — стационарная случайная функция, заменяют эквивалентным в некотором смысле линейным уравнением

х + Рх + (02х = / (/)•

(54)

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

>

п ,

XV (х, х)

(55)

■= ю^ +

ц ----- Ц - ! - .

 

 

X2

 

Для определения средних

значений, входящих в формулы

(55),

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

/(О — белый шум интенсивностью Р 0, а V (х, х) = х3. Тогда

XV (х, х) = 0 ;

р =

Р;

хУ (х, х) = х4 *

3 (х2)2;

ш2 = ©о +

 

2 =

“ о + Зм------—

(56)

3 ^

*

 

 

 

2р©а

 

Последнее соотношение (56) решается относительно ш2. В резуль­ тате с точностью до членов порядка ца получаем

со2 = с»о +

 

2^ ;

 

(57)

X

Ро Л 3|хР0 \

2| Ц \

2( Ц ) '

 

Нелинейные случайные колебания

539

Результаты^ исследования уравнения (53) для различных типов

функции

V (х, х) и различных возмущающих процессов / (/) приведены

в работах

[40, 41, 54]. В работах [44, 53] даны результаты исследова­

ний точности метода статистической линеаризации применительно к не­

которым простым системам. Эти исследования основаны на сравнении

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

марковских процессов. Вычисления показывают, в частности, что точ­

ность метода тем выше, чем меньше интенсивность возмущающего про­

цесса (с увеличением.последней возрастает «эффективная нелинейность»

системы).

 

 

Идея метода статистической линеаризации (априорное введение функ­

ции распределения, зависящей от конечного числа параметров, которые

далее находят из условия минимума ошибки) нашла применение и при

решении

детерминистических задач [18]. Аналогичную идею исполь­

зуют по существу и при исследовании нелинейных стохастических

уравнений методом моментов. Для того чтобы замкнуть цепочку урав­

нений для моментов, постулируют зависимость между моментами раз­

личных порядков, соответствующую выбираемому закону распределе­

ния. В рассмотренном выше примере система уравнений для моментов

второго и четвертого порядков замыкается при помощи дополнитель­

ного соотношения х4 = 3 (х2)2. В результате приходим к выражению (56)

для ш2.

 

 

Метод малого параметра. Пусть в уравнении (53) V (х, х) — поли­

ном относительно х, х. Полагая, что

 

 

х (0 = *0 (0 + рхг (0 + р3х2 (/) Н-----

 

получаем линейные стохастические дифференциальные уравнения от­

носительно хь (/), в правые части которых входят функции х/ (/) (/ <

к).

В результате функции хд {() могут быть выражены в виде квадратур

от возмущающего процесса [ ((). После этого можно найти моменты

процесса х (0 » например:

 

 

Р = ^ + 2 ц х 0 х 1 + ( р ( л ^ + 2 з д [ ) Н ------

 

Метод малого параметра применяли к системе с нелинейным демпфи­

рованием

[44] и к нелинейной системе с двумя степенями свободы [41 ].

В работе

[33] этим методом решена задача о нелинейных колебаниях

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

Метод малого параметра может быть применен к параметрическим задачам, в которых случайные функции входят в коэффициенты диффе­ ренциального уравнения. Некоторые задачи устойчивости линейных систем со случайно изменяющимися параметрами исследованы при помощи этого метода в работах [35, 38, 48]. В работе [16] рассмотрена задача, приводящая к уравнению

* + Р* + шо П + ^5 (0] * = а а х а>д(;

(58)

540 Статистические задачи колебаний и устойчивости

здесь | (0 — стационарная случайная функция, для которой

ш = 0; Е('0 Е(* + т)= /С (т).

Использование разложения (57) с последующим применением изло­ женной выше схемы вычислений дает с точностью до членов порядка р,4

а

Г

°>о

1

ф),

 

 

 

5Ш (<о0* +

где ф — постоянная

фаза, а

 

 

 

^ ± = ^ К (т)е

Рг/2 соз(о0 ± со^т Лт;

 

В работе [16] показано, что если

^ 0 при со >

0, гдеФ (со) —

спектральная плотность процесса

^+^> 0. Следовательно,

при этом условии случайные изменения собственной частоты Приводят к уменьшению «в среднем» амплитуды колебаний.

Метод усреднения. Этот метод использует известные идеи Крылова— Боголюбова в теории нелинейных колебаний. Если исследуемый колеба­ тельный процесс имеет узкополосный спектр, то уравнения движения могут быть усреднены за «период» колебаний. Затем применяют либо корреляционную теорию, либо теорию марковских процессов. Подроб­ ное изложение метода усреднения применительно к случайным функ­ циям содержится в монографии [27], где рассмотрено большое коли­ чество нелинейных и параметрических задач.

МЕТОДЫ ТЕОРИИ МАРКОВСКИХ ПРОЦЕССОВ

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

ние вероятностей в момент

зависит от распределения в предшеству­

ющий момент ^

^1 и не зависит от истории системы.

Уравнение Фоккера—Планка—Колмогорова. В качестве простей­ шего примера рассмотрим простой марковский процесс с одной пере­

менной х (/). Введем обозначения

х = х ( 0 и хх =

Можно

показать, что если пределы выражений

 

«А (0 = Нш

(*Т— *)*

(И)

*->о

т

 

Методы теории марковских процессе

541

существуют, то плотность вероятности р (х) подчиняется кинетическому уравнению

во

 

р (*) = 2 -^Г (~

)*1Х*(х)Р <*))•

(6°)

 

А= 1

 

 

Функции

(х) называют интенсивностями марковского процесса.

Из них к 1(х)

характеризует среднее

течение процесса;

х 2 (х) — его

дисперсию и т. д. Коэффициенты интенсивности связаны с корреляци­ онными функциями производной х(/) зависимостями

*2. *Л)= Х * (х ) б (* !-* .) б ( /! - /* ) + . . . (61)

(точками обозначены члены, содержащие более слабые особенности). Если все у./г (х) при /&> 2 равны нулю, то марковский процесс на­

зывают непрерывным. Уравнение (60) принимает вид

 

р (х) = — - ^ г К (*) Р (*)] + - | - •

Iх®(*) р (*)1

(62)

и носит название одномерного уравнения Фоккера—Планка. Оно опи^ сывает, например, случайное блуждание точки вдоль оси. Уравне­ ние (62) может быть обобщено на случай многомерных непрерывных

марковских

процессов

 

 

 

 

 

Р(Х1,

 

Х $)

 

 

5

 

 

 

 

 

[ха (лч,

х%.

хв)р(х1, хг.

*5> ] +

 

дха

 

 

 

 

5

$

 

 

 

 

 

а2

 

 

*$) Р (XI. *2,

*)]; (63)

 

дхадхр [*а0 (*Ь *2.

здесь

 

т

(хах ха)

 

 

 

 

 

 

->0

 

Т

 

 

 

 

 

 

(64)

Уравнение (63) называют обычно уравнением Фоккера—П ланка-

Колмогорова (уравнением

ФПК).

 

 

Пусть движение системы описывается стохастическими дифферен­

циальными уравнениями

 

5

(65)

 

 

 

 

 

 

2

 

Ха + ёа (Х1....... хз) +

Х$**^~ ^

 

 

 

Э-1

 

где

(0 — независимые белые

шумы;

 

Ь Ш

Ж

) =

6 (т) 6а/}.

Источник: https://tut-files.ru/previewfile/74208