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

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

Корреляционные методы исследования распределенных систем 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)

Корреляционные методы исследования распределенных систем 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)

кх Дку

в (кх,

ку) т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] рассмотрены другие случаи пространственной и вре­ менной корреляции, исследовано влияние инерции вращения и дефор­ маций сдвига, получены формулы для случайных колебаний оболочек, произведены подробные вычисления для сферической оболочки. Огра­ ничимся тем, что приведем формулу для среднего квадрата нормального

Источник: https://studfile.net/preview/19992756/