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

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

Корреляционные методы исследования линейных систем

529

Колебания систем с конечным числом степеней свободы под дей­ ствием случайных сил. Допустим, что движение механической системы описывается уравнениями

Оа+ 2гаыаиа + “аиа= СО(а= 1 , 2 , .. .),

(26)

где еа — коэффициенты демпфирования; юа — частоты

свободных

колебаний; точками обозначено дифференцирование по времени. Фор­ мула (25) принимает вид

Ф„ ’Э (©)

 

Ф,

(27)

(са2 + 2*еасоюа — (О2) (о)р 2ш&рй>р — со2)

 

По санкции

^

(со) можно определить вторые моменты обобщен­

ных координат.

Учитывая,

 

получим, применяя

обратное преобразование Фурье,

 

 

 

оо

 

 

 

^ = 4 1 ®к,в.в (т ) л в -

(28)

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

нат имеет вид

 

 

 

 

ис1^а1/ц ас11ип{

|

г

. . .

(29)

-------------=

-1

ф.(

(СО) (ш)*+/ Асо.

А1к Л11

2

^

 

 

Вычисление интегралов (28) и (29) для случаев, когда Ф -,^ (ш)

определяется соотношением (27), а спектральные плотности обобщен­ ных сил Ф0 0 (со) являются дробно-рациональными функциями,

может быть выполнено при помощи теоремы вычетов. В приложении к книгам [24, 26) содержатся таблицы интегралов вида

1 - 1 &.(»*>) ^

где кп (х) и дп (*') — полиномы. При помощи этих таблиц можно опре­ делять средние квадраты обобщенных координат и их производных.

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

3.

Вторые моменты обобщенных координат и обобщенных скоростей в системе с двумя степенями

 

свободы, находящейся под действием дельта-коррелированных случайных сил

Вели­

 

 

чина

°а

 

 

 

п аа

 

 

 

 

4ео“ о

 

 

 

 

0

 

 

 

 

Лаз

 

(

^

) 2

 

 

2 (еаша +

е^евр)

^ а^Р (Баа)Р + еРма)

 

 

 

°ар

 

..2

..2 , “ ““В

еа“ 0 + ев“ а

“ “

Ш» +

 

еаша + еВше

 

0

 

 

° а а

 

 

4еа“ а

 

° а Р

 

» * - * +

 

е““ в+Ев“ а

0 “

Ш2 _ т2

гаш0+ г вшр

 

°аВ

 

( <«>2 —0)о) “

 

(еа й)Р+ е вй)а)

+ 2 ^ а<°а+ЕЗс,)Р)

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

Корреляционные методы исследования распределенных систем 531

В табл. 3 приведены формулы для вторых моментов обобщенных координат и скоростей системы (26) для случая, когда 0а (() — белые шумы с интенсивностями Оаа и взаимными интенсивностями Оар. При

этом Оаа = Фаа/л, ОаЗ = Фар/л. гДе Фаа» Фа0 — спектральные плот­ ности и взаимные спектральные плотности соответственно. Из фор­

мул табл. 3 следует, что при малых значениях коэффициентов затуха­ ния, характерных для упругих систем (еа < 1), взаимная корреляция различных обобщенных координат и скоростей оказывается заметной лишь при весьма близких значениях частот 0)а и ©р. Например, при еа = ер = е

— ___________2Раре (ша + сар)________

(30)

К- °>р) + 4е2й>а“ р ( “ а +

Следовательно, если

V < | соа — свв |,

то значения уао$ оказываются пренебрежимо малыми по сравнению с главными элементами матрицы вторых моментов.

Если демпфирование мало, то для средних квадратов обобщенны координат может быть дана приближенная формула

2 __ пФУа<1а

(31)

Vа

4еа»а

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

Уравнения колебаний линейных распределенных систем. Эти урав­ нения могут быть записаны в виде

= у;

здесь <и>(.х, у, г, /) — вектор перемещений; щ (х, у, г, /) — вектор на­ грузок; Ь — некоторый линейный оператор. Пусть компоненты век­ тора у являются случайными функциями координат и времени с извест­ ными вероятностными характеристиками. Задача состоит в определении вероятностных характеристик вектора о» и его производных по коорди­ натам и времени. Для простоты в дальнейшем будем рассматривать уравнение

= у,

(32)

где ю и у — скалярные функции координат и времени. Уравнение (32) можно трактовать как «разрешающее» уравнение в теории колебаний стержней, пластинок и оболочек.

Сведение к системе с конечным числом степеней свободы. В соответ­ ствии с этим методом перемещение ни представляется в виде разложения в ряд по формам собственных колебаний — по собственным функциям Фа (*» У, 2) оператора

со

 

ш = 2 «О (0 Фа (*. У. г).

(33)

а=1

 

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

После подстановки ряда (33) в уравнение колебаний (32) получаем бесконечную систему обыкновенных линейных дифференциальных уравнении. Корреляционные функции обобщенных сил Фа (*) опреде­ ляются соотношениями

^ =

Л

^

^ 1’ Уи

У*•

 

Хфа(*1> у

2^)

фр (л*2, Г/а» 2Г г)^1 ^

(34)

здесь Кп — корреляционная функция нагрузки, т. е.

К(! = Ч(х1. Уи *1, <|) Ч (хг> У*. ?2 >У ;

(35>

\*а — нормы собственных функций фц (х, у, г). Интегрирование ведут по загруженной поверхности тела 5.

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

чек дано в работах [5, 7, 14, 39].

можно

После нахождения моментов обобщешшх координат К..

С Г 0

 

определить моменты других параметров, выражаемых через обобщенные

координаты. Например, средний квадрат 5*2 некоторого фактора 5 опре­ деляют по формуле

= Е Ц

Саср/С „ ,

(36)

Ц=1 0=1

а Р

 

где са — коэффициент влияния, входящий в формулу для фактора 5:

оо

5 = 2 СаУа, и—1

Нагрузка, которая является стационарным временным случайным процессом. В этом случае корреляционные функции (34) и (35) зависят от временного интервала т = /2 — Если выполняется условие эргодичности, то можно ввести преобразование Фурье по времени от К

Ф<7(*ь /А,

Уъ, г2; (0) =

 

оо

 

 

= -^- | К я ( *1, У и

Уг.

т ) е ~ ‘ш х<1т(37).

----ОО

 

 

Функцию Ф^ будем называть спектром пространственных корреля­ ций. Взаимные спектральные плотности Фд ^ обобщенных сил выра­

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

<Ч«в(с,,)

у"

ш)ху»

X Фа (а*11 У и г1 ) Фр (х и У и 2,) с!3± (132.

(38)

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