Корреляционные методы исследования распределенных систем 531
В табл. 3 приведены формулы для вторых моментов обобщенных координат и скоростей системы (26) для случая, когда 0а (() — белые шумы с интенсивностями Оаа и взаимными интенсивностями Оар. При
этом Оаа = Фаа/л, ОаЗ = Фар/л. гДе Фаа» Фа0 — спектральные плот ности и взаимные спектральные плотности соответственно. Из фор
мул табл. 3 следует, что при малых значениях коэффициентов затуха ния, характерных для упругих систем (еа < 1), взаимная корреляция различных обобщенных координат и скоростей оказывается заметной лишь при весьма близких значениях частот 0)а и ©р. Например, при еа = ер = е
— ___________2Раре (ша + сар)________
(30)
К- °>р) + 4е2й>а“ р ( “ а +
Следовательно, если
2е V < | соа — свв |,
то значения уао$ оказываются пренебрежимо малыми по сравнению с главными элементами матрицы вторых моментов.
Если демпфирование мало, то для средних квадратов обобщенны координат может быть дана приближенная формула
2 __ пФУа<1а |
(31) |
Vа
4еа»а
ПРИМЕНЕНИЕ КОРРЕЛЯЦИОННЫХ МЕТОДОВ К ИССЛЕДОВАНИЮ КОЛЕБАНИЙ РАСПРЕДЕЛЕННЫХ СИСТЕМ
Уравнения колебаний линейных распределенных систем. Эти урав нения могут быть записаны в виде
= у;
здесь <и>(.х, у, г, /) — вектор перемещений; щ (х, у, г, /) — вектор на грузок; Ь — некоторый линейный оператор. Пусть компоненты век тора у являются случайными функциями координат и времени с извест ными вероятностными характеристиками. Задача состоит в определении вероятностных характеристик вектора о» и его производных по коорди натам и времени. Для простоты в дальнейшем будем рассматривать уравнение
= у, |
(32) |
где ю и у — скалярные функции координат и времени. Уравнение (32) можно трактовать как «разрешающее» уравнение в теории колебаний стержней, пластинок и оболочек.
Сведение к системе с конечным числом степеней свободы. В соответ ствии с этим методом перемещение ни представляется в виде разложения в ряд по формам собственных колебаний — по собственным функциям Фа (*» У, 2) оператора
со |
|
ш = 2 «О (0 Фа (*. У. г). |
(33) |
а=1 |
|
Корреляционные методы исследования распределенных систем 533
Статистические характеристики случайных нагрузок. Для решения практических задач в рамках корреляционной теории необходимо знать
математическое |
ожидание |
ц (х, |
у, г, /) и корреляционную функцию |
||
Кд(хг, Ух» гх» |
Уч» |
гг» У |
нагрузки или |
в случае эргодичес- |
|
кой стационарной нагрузки — спектр |
пространственных корреляций |
||||
С*1» У\> *Г> |
*2» У2. 22; |
со). В |
случае |
нагрузок, |
вызванных акусти |
ческими шумами, турбулентными пульсациями давления в пограничном слое, нерегулярным волнением моря, спектры давления в отдельных точках поверхности находятся из эксперимента. Соответствующие дан ные приведены в работах [17, 19, 42, 43, 57]. В задаче о движении авто мобиля по неровной дороге спектр возмущений может быть вычислен по спектру длин волн неровностей, определенному путем геодезических измерений.
Полное экспериментальное определение корреляционной функции или спектральной плотности пространственно-временного случайного процесса связано, как правило, со значительными трудностями. По этому для приближенных оценок обычно используют простейшие аналитические зависимости. Некоторые из них указаны ниже.
Д е л ь т а к о р р е л и р о в а н н а я в п р о с т р а н с т в е н а г р у з к а . Это предельный случай нагрузки с весьма малым (по сравнению с длинами волн возбуждаемых форм колебаний) масштабом
пространственной корреляции. В этом случае |
|
|
Уи *1» *2. 1/2» V. |
&>) = |
|
= V (щ) 6 (х2 — хг) Ь {у2 — ух) Ь (г2 — гг), |
(39) |
|
где У (со) — заданная функция частоты; |
б (х) — дельта-функция. |
|
При использовании такой модели весьма упрощаются вычисления. Например, для спектральных плотностей Фд д^(а) в случае ортого
нальных |
форм собственных |
колебаний получаем, что |
|
|
Ф, |
(ш) = ЧГ(со) |
(40) |
П о л н о с т ь ю к о р р е л и р о в а н н а я |
в п р о с т р а н |
||
с т в е |
н а г р у з к а . Это |
противоположный |
случай большого (по |
сравнению с размерами конструкции) масштаба пространственной кор реляции. В этом случае спектр пространственных корреляций не зави сит от координат.
Н а г р у з к а в о л н о в о г о т и п а . Эта нагрузка характери зуется тем, что ее корреляционная функция зависит от времени и коор
динат |
через |
посредство |
параметра: |
|
|
||
|
|
т |
т |
* г — *1 |
У г — У 1 |
* » — * ! - |
|
|
|
|
1 |
сх |
Су |
сг |
* |
здесь |
сх, Су, |
сг — постоянные |
(скорости |
распространения волн на |
|||
грузки). В этом случае |
|
|
|
|
|||
|
|
|
|
*Г» |
0 ) = |
|
|
534 Статистические задачи колебаний и устойчивости
Нагрузки такого типа встречаются в задачах о воздействии на кон струкцию дальнего акустического поля или турбулентных гидродина мических пульсаций струи [43, 46, 49], о движении в статистически неоднородной среде, о движении нагрузки по мосту [4] и т. п.
А к у с т и ч е с к и е н а г р у з к и . Задача о колебаниях упру гих систем под действием акустического излучения работающих двига телей приобрела в последние годы большую важность в связи с так на зываемой проблемой «акустической усталости» конструкций [6, 43]. Экспериментальные данные по частотным спектрам пульсаций давления в различных точках акустических полей работающих двигателей приведены в работах (42, 43, 49]. Пространственную корреляцию в принципе можно рассчитывать в соответствии с теорией Лайтхилла [52], исходя из решения неоднородного волнового уравнения. В некоторых случаях, однако, пространственную корреляцию можно оценивать на основании чисто геометрических соображений [32].
Метод интегральных оценок. При применении корреляционных ме тодов к решению задач о колебаниях распределенных систем в ряде слу чаев возникают трудности, связанные с необходимостью проводить суммирование по весьма большому числу возбуждаемых форм колеба ний. Тогда эффективное приближенное решение может быть получено по методу интегральных оценок [4, 7, 39]. Суть метода заключается
впереходе от суммирования в формулах типа (36) к интегрированию
впространстве некоторых параметров. В соответствии с теорией дина мического краевого эффекта [5 ] такими параметрами являются волновые
числа кх, ку, кг форм колебаний.
Рассмотрим, например, пластинку, находящуюся под действием нормальной нагрузки. С учетом формулы (31) имеем
(41)
Согласно методу интегральных оценок выражение (41) приближенно
записывают в виде |
л |
|
|
|
|
|
|
<г |
я |
Д2 (кх, |
Ы Ф (кх, |
ки) |
йкх &ку\ |
(42) |
|
4Дкх Дку |
в (кх, |
ку) т3 (кх, |
ку) |
||||
здесь Д&х'= |
-2- Д |
а |
спектральная |
плотность обобщенных |
|||
сил Ф0 0 (ш), собственные частоты ша , коэффициенты демпфирова-
ния еа и коэффициенты влияния са трактуются как функции волновых чисел, принимающих все значения из первого квадранта.
Применения метода интегральных оценок. В случае пластинки по стоянной толщины, загруженной нормальной дельта-коррелированной на срединной поверхности нагрузкой с учетом формулы (40), получим следующую оценку для прогиба во внутренней области:
I 1 / Р |
\» ГТМА-. |
(43) |
~ 80а \ рЛ ) |
^ е (г) г ’ |
|
Корреляционные методы исследования распределенных систем 535
здесь Р — цилиндрическая жесткость; рк — плотность пластинки, отнесенная к единице площади срединной поверхности;
1
(44)
Для максимальных нормальных напряжений во внутренней области находим оцеику
а2 |
(3 + 2у + Зу’-) | |
Ч' (г) (1г |
’ |
(45) |
|
|
е(г)г |
|
где V — коэффициент Пуассона. Для максимальных напряжений в за делке имеем
|
|
|
1 |
|
|
О2 |
9 |
/_ 0 _ \ |
2 Г У (г)б.г |
(46) |
|
Ж |
\ рк ) |
] е (г) г |
|||
|
|
и т. д. В табл. 4 для ориентировки приведены формулы для средних квадратов о>2 и о2 и эффективной частоты
со
1С02Фд ((О) Й<0 |
2 |
(47) |
|
П |
|
входящей в формулы для определения долговечности при циклических нагрузках (см. гл. 8, т. 1). Формулы выведены для усеченного времен ного белого шума
Ч'(со)=( С 0 П 8 1 |
(со < |
|
(48) |
|
1 |
0 |
(йн или СО> (Од). |
||
Кроме того, введены |
обозначения |
|
|
|
|
1 |
( |
Р \ |
2 V |
|
ю 'р " V\ Рлр ь )/ |
> |
||
|
9 |
|
|
V |
|
1бк4 (3+2у + Зу2)-2_; |
|||
гн и вн — значения параметров |
г и е |
при |
со = (йн. |
|
В работе [39] рассмотрены другие случаи пространственной и вре менной корреляции, исследовано влияние инерции вращения и дефор маций сдвига, получены формулы для случайных колебаний оболочек, произведены подробные вычисления для сферической оболочки. Огра ничимся тем, что приведем формулу для среднего квадрата нормального