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

книги / Статистические методы в строительной механике

..pdf
Скачиваний:
6
Добавлен:
12.11.2023
Размер:
18.02 Mб
Скачать

Рассмотрим вначале одну случайную функцию -X(tf). Ее кор­ реляционная функция

7/2

^ x x (T) = lim

~r

J x(t)x(t + x)dt

(5.10)

Т-+00

—TJ2

 

 

 

является согласно формуле

(5.7)

четной функцией х. Допустим,

что она может быть представлена при помощи косинус-преобра­ зования Фурье

 

 

(5.11)

Обратное преобразование имеет вид

 

 

оо

 

Фхх (ш) =

J Кхх (т) cos о)т dt.

(5-12)

 

о

 

Вновь введенная функция

Фхх (со), зависящая от частоты со,

называется спектральной плотностью стационарного случайного процесса X(t). Учитывая (5.8) и (5.10), получим, что

00

Ф хх(«)cto.

(5.13)

о

Для установления смысла спекгральной плотности заметим, что вели­

' X X

чина X‘2(t)

пропорциональна

средней

 

 

 

мощности

случайного процесса.

Сле­

 

 

 

довательно, произведение

Фхх

(to)rfa

 

 

Ф(Ш)

соответствует той доле мощности,

ко­

 

 

 

 

 

торая заключена

в интервале частот от

 

 

 

(о до (D+ dœ. Поэтому функция

Фхх (ш)

 

 

I

называется

иногда также

спектраль­

 

 

ной мощностью, спектральной плотно­

 

 

I

стью энергии и т. п.

 

 

 

-Ш,

о

-1— CL)

 

 

 

Ь)г

Рассмотрим

некоторые типичные

 

 

 

примеры

спектральных

плотностей.

 

Рис.

65

Если случайный процесс не коррелиро­ ван, то корреляционная функция обращается тождественно в

нуль всюду, за исключением точки т=0. Ее можно выразить, та­ ким образом, через дельта-функцию 6(т):

К ххМ = КаЬ(т),

Ф хх(®) = &

(5.14)

Здесь Ко— некоторая константа.

ТЪ

 

 

Из второй

формулы

видно,

что каждая частота вносит одинаковый вклад в мощность слу­ чайного процесса. Эта наиболее хаотичная случайная функция называется «белым» шумом. «Белый» шум, являющийся физи-

ческим аналогом «белого» света в оптике, представляет собой все же чрезмерную абстракцию, поскольку его полная мощность, пропорциональная интегралу

00

J Ф*х(ш)сК

о

оказывается бесконечно большой. Чтобы мощность процесса оставалась конечной, вводится ограничивающая частота шс (рис. 65).

Противоположным случаем является процесс с детермини­ рованной частотой 0. В этом случае

Фхх (ш) = Со(в — | ш|),

т. е. вся мощность процесса сосредоточена в бесконечно малой окрестности частоты 0 (рис. 66).

Рассмотрим, наконец, весьма важный случай, когда подав­ ляющая часть мощности процесса заключена в некотором ин­ тервале вокруг преобладающей частоты процесса 0. Корреля­ ционная функция должна иметь вид затухающей кривой с ча-

Ф(и)

Рис. 66

стотой 0 ; затухание этой кривой соответствует ослаблению кор­ реляционной связи по мере увеличения интервала т (рис. 67). Одно из возможных выражений для Кхх (•) имеет вид

Кхх С1) = Ко

х 1cos 0т,

(5.15)

где /Со и а — константы. Соответствующая-

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

ность

 

 

 

Фхх ("> =

к

 

 

g/Cp

+

1

____ «*о

 

(5.16)

 

(ш — 0)» + «а

тс

(а) + 0)® -J- а*

 

имеет максимумы вблизи частот в ='±<о.

случай нескольких

кор­

Результаты

легко

обобщаются

на

релирующих между

собой

случайных

процессов Хх (t), Х2 (£),

Xa(t). Соответствующие спектральные плотности вводятся

формулами:

 

 

 

 

00

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К x . x k (т)

=

J

ФXjXh (ш) COS ют du,

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

(5.17)

 

 

 

 

 

 

00

 

 

 

 

 

 

ф XjXk (<0) =¥

J

/СхухА(т) cos шт dï.

 

 

 

Практическое

 

определение

 

 

 

 

 

 

спектральных плотностей

по эм­

 

 

 

 

 

 

пирическим

данным

сводится

0 «?'.*♦

 

 

 

обычно

к численному

определе­

 

 

 

нию корреляционных

функций и

 

 

\ о

 

 

 

выполнению

 

преобразования

В Ю'-

 

^ Плои.{ода

Фурье. Операции могут быть ме­

 

 

 

ханизированы и, в частности, мо­

 

-4

 

 

41 Уо

 

гут выполняться общими и специ­

5 Ю

 

 

 

 

ализированными

электронными

 

 

 

 

 

 

вычислительными машинами. Не­

 

1-4

 

 

 

 

которые сведения об этом мож­

 

 

 

 

 

 

 

 

 

 

 

но найти в книгах (85, 106].

 

 

 

 

 

 

 

Для некоторых реальных про­

210-4

 

 

 

 

цессов

были

построены

кривые

 

 

 

 

 

 

спектральных

плотностей.

На

 

 

 

 

 

 

рис. 68

показана

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

 

 

2000

то

6000

8Û0\

плотность для

пульсаций

давле­

 

 

ния в аэродинамическом следу за

 

 

 

 

C J '/с е к

плохо

обтекаемым

телом

[70].

 

 

Рис.

68

 

 

Процесс состоит из «фона» с бес­

 

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

составляю­

порядочными пульсациями и

чисто

 

щей, заключающей в себе около 41°/о мощности всего процесса.

51. Случайное воздействие на простейшую линейную систему

Прежде чем переходить к анализу поведения упругих си­ стем под действием случайных сил, остановимся на хорошо изу­ ченном случае линейной колебательной системы с одной сте­ пенью свободы. Пусть / — обобщенная координата, ыа — собст­ венная частота, р — характеристика затухания, равная отноше­ нию логарифмического декремента затухания к числу 2я. Обо­

значив обобщенную случайную силу через Q (t), получим урав­ нение движения

/ + 2рш0/ + шо / = Q (0-

(5-18)

Его решение может быть записано при помощи импульсной пе­ реходной функции h(t—т), равной величине обобщенной ко­ ординаты ] в момент времени t от единичного импульса, прило­ женного в момент времени т:

t

f { t ) = j h {t T) Q (X) d z .

(5.19)

Вычислим при помощи формулы (5.6) корреляционную функ­ цию обобщенной координаты К// (£ц h) :

 

*/,(<!. ЧЬ = № № .

 

(5-20)

Подставляя в формулу

(5.20)

выражение

(5.19),

найдем

 

K f / (t„ y =

j

I f t p . — ч ) А ( / „ — <>)Х

 

 

 

—00 —00

 

 

 

 

X K QQ ы

z2) d z l d z2,

 

(5.21)

где K QQ

(ть T2) — корреляционная функция

для

обобщенной

силы.

Если внешняя

нагрузка

начала действовать в момент

времени t=0, а при £<0 система находилась в покое, то в ка­ честве нижних пределов интегрирования в формуле (5.21) сле­ дует взять t\= /2—0.

Рассмотрим подробнее стационарное случайное воздействие. Производя в формуле (5.21) замену переменных

 

tx

Т, =

г

г

Х»

 

 

ti, t%-- = Х2,

-- ^1 =

 

получим

оо

оо

 

 

 

 

 

 

 

 

 

 

/С//(х) = j

j* h

(xj) h (x2) K QQ (x — x2 + xi)

dxj d z 2.

(5.22)

оô

Вдальнейшем штрихи при tî и x2 опускаем.

Составим преобразование Фурье от левой и правой частей равенства (5.22). Замечая, что1

ОО

/CQ Q (T) = - ^ - J Ф<?<2 (ш)е'шт do),

*>/W

Л ® ,

 

1 Здесь используется то обстоятельство, что спектральная плотность - четная функция частоты <•>.

и изменяя порядок интегрирования, придем к соотношению

во

j =

00

00

= J

dm I

j h (tj) h (т2)

®QQ (ш) еш(T_Tj+Tl) dxt dx2.

(5.23)

00

о

о

 

 

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

передаточную функцию линейной

системы

F (tco), равную преобразованию Фурье от импульсной

переход­

ной функции:

 

 

 

 

 

 

 

00 "

 

 

 

F (im) =

J Л ('с) é~iandx.

(5.24)

Передаточная функция проще всего определяется как отноше­ ние установившегося решения линейной системы при гармони­ ческом воздействии exp(fW) к этому воздействию. Так, для уравнения (5.18)

F (!») =

--------!-----

;-----.

(S.25)

 

wjj— оГ

2(V<»a>0

 

Разумеется, можно выразить передаточную функцию и через импульсную переходную функцию, подставляя в формулу (5.24) выражение

Цх) =

 

О,

если т < 0,

 

(5.26)

1

е |о,°т sin (w0 Y

1— р2х),

если х >

“О/ Г н ?

0.

 

 

 

 

 

Принимая

во внимание

формулу

(5.24),

вместо

(5.23) по­

лучим соотношение

 

 

 

 

 

00

 

 

00

 

 

 

J Ф/7 («О еш dm = —00

j F (im) F ( — im) 0 QQ (ш) eimxdm. (5.27) —00

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

00

^ h (х) etwx* dx = F ( im).

о

Соотношение (5.27) должно выполняться при любом т; от­ сюда

Ф( ш) = F (im) F ( — im) ФQQ (m).

(5.28)

Формула (5.28) устанавливает весьма простую связь между спектральной плотностью обобщенной силы Q(t) и спектраль­ ной плотностью обобщенной координаты f(t) (между спек­ тральными плотностями «входа» и «выхода» линейной си­

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

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

~ мулам (5.11) и (5.13) подсчитать значения кор­ Рис. 69 реляционной функции и средний квадрат «выхо­

да».

Для системы, описываемой уравнением (5.18) с передаточ­ ной функцией (5.25), формула (5.28) примет вид

Ф// (“>)

Фрр (<■>)

(5.29)

(„ * _ в«)« + (2р»од»

 

 

Если затухание не очень велико, то система выполняет роль «фильтра», пропуская главным образом пульсации с частота­ ми, близкими к woЗдесь имеет место явление, похожее на ре­ зонанс вынужденных колебаний при детерминированном перио­ дическом воздействии (рис. 69). _

Средний квадрат обобщенной координаты f2 найдем по фор­ муле (5.13)

фд<? (<■>)

(5.30)

(<»0 — 4- (2fkü<o0)2

52. Колебания упругой системы при случайном воздействии

Исследование поведения упругой системы представляет не­ которые трудности по сравнению с изложенным выше простей­ шим случаем. Это происходит главным образом потому, что колебания упругих систем описываются дифференциальными уравнениями в частных производных. Так, изгибные колебания прямолинейного стержня, загруженного поперечной нагрузкой д(х, /), в рамках обычных предположений сопротивления мате­ риалов описываются уравнением

 

J L

(El —

) 4 - [F diw

— Я(х>0 -

(5.31)

 

dx*

\

dxi )

dt*

 

Здесь w (x,

t) — прогиб

 

стержня,

El — изгибная

жесткость,

F— площадь

поперечного

сечения,

 

р — плотность

материала.

Изгибные колебания пластины

постоянной

толщины в рамках

обычных предположений описываются уравнением

 

 

 

ШДш + р/г ^

=

 

ÿt

t),

(5.32)

где w(x, у, t) — нормальный прогиб, D — цилиндрическая жест­

кость, h — толщина пластины. В обоих

случаях

затухание не

учитывается. Более сложные примеры

получим,

рассматривая

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

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

Пусть w(x, у, 2, t) — динамическое

перемещение

точек си­

стемы, wK{x, у, z) — формы собственных

незатухающих колеба­

ний системы. Они удовлетворяют уравнению

 

в случае изгибных колебаний стержня, уравнению

 

£>Д Д wk—p/îi4 wfc= 0

 

в случае изгибных колебаний пластины и т. д. Здесь

шА— соот­

ветствующая частота собственных колебаний.

Представив перемещение w(x, у, z, t) .в виде ряда

 

w(x, у, z , 0 = 2

У* 2).

(5-33)

к=Л

 

 

легко придем к уравнениям, описывающим поведение обобщен­ ных координат fh{t):

fk + 2j3ft cuAfk +

щ fl{ — Qk(0

(5.34)

(A = 1, 2,

.

. я).

 

Здесь Qk (t) — обобщенная сила, которая в общем случае трех­

мерного тела объемом V, загруженного

на

поверхности S на­

грузкой q(x, у, z, t), определяется по формуле

 

 

 

J J q(x,

у, г,

i) wk {x,

у, z) dS

 

QAt) =

J J

J

 

Z)d V

(5.35)

 

 

 

y '

 

 

 

 

V

 

 

 

 

Члены, содержащие обобщенные скорости fk , введены уже

после перехода к

обыкновенным

дифференциальным уравне­

ниям;

— характеристики

затухания,

равные логарифмиче­

ским декрементам

 

соответствующих форм

собственных коле­

баний, разделенным на 2я. Следовательно,

здесь использовано

предположение, что затухание — линейное и

что диссипативные

связи между отдельными формами колебаний отсутствуют. Для приближенного анализа это предположение является вполне приемлемым.

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

вообще говоря, нужно учитывать и

 

корреляцию

между раз­

личными степенями свободы. Рассмотрим этот

вопрос подроб­

нее.

 

 

нас

интересует

 

средний

квадрат

прогиба

Допустим, что

 

w(x, у, 2, t). Из формулы (5.33) следует, что

 

 

 

я»2

(х,

У,

г, t) = w

(х, у, z,

/) w

(.х,

у,

z, t) =

 

=

2 2

K fjffl'

*)тЛх>У>

 

У’ *).

(5.36)

 

 

/=1 А=1

 

 

 

 

 

 

 

 

где Kfjfn

(tu

h) — корреляционные

 

функции второго

порядка

для обобщенных координат

f j ( t ) .

 

Аналогично

для

среднего

квадрата изгибающего момента в сечениях

стержня

получаем

формулу

 

 

 

 

________________

 

 

 

МЦх, t) = (El)2 дга> <х,

t)

&ш<х, t)

 

 

 

 

 

 

 

 

дх2

 

дх2

 

 

 

= {Е1Г 2

 

 

d2Wj (х)

d*wk (х) .

(5.37)

2

 

*> ~ d x 2

dxa

 

 

 

 

 

 

/=1 *=1

 

 

 

 

 

 

 

 

Таким образом, необходимо уметь

вычислять

корреляционные

функции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

&

= h

('.) h

('»)

 

 

(5.38)

 

 

 

 

 

 

 

 

 

 

по крайней мере для значений аргументов ti= t2=t.

К необходимости изучения взаимных корреляционных функ­ ций обобщенных «координат можно подойти также иначе. До­ пустим, что интересующие нас случайные функции Vj(t) (на­ пример, напряжения, перемещения и т. п.) выражаются в виде линейных функций обобщенных координат:

<*/.('>

 

 

k - i

 

 

 

(/ = 1,

2,

. m),

 

где Cjk— некоторые не

зависящие от

времени

коэффициенты.

Тогда корреляционные

функции

для

v} (/) определяются по

формуле

 

 

 

 

 

к , , к (<„ t,) = V

2

с'*“* КЩ (,‘>

<5-39>

 

0 = 1

Р =1

 

 

 

Если внешнее воздействие распределено по нормальному за­ кону, то согласно свойствам нормального распределения пара­ метры Vj, полученные в конечном счете путем некоторого ли­ нейного преобразования внешнего воздействия, тоже будут распределены по нормальному закону. Отсюда совместная я-мерная плотность вероятности для параметров Vj (/) в каж­ дый фиксированный момент времени t согласно формуле (5.3) будет

 

 

Р К

V 2 ,

 

vm) =

 

 

 

 

1

■ер Г— -i- 2

2 {*0 ' (<■ О }

/

*

• (5.40)

 

 

V (2ч)" \Kv (t, t)\

L

^ j=l A=1

 

J

 

Здесь K7X

0 — матрица, обратная матрице Kv

{t, t)

с эле­

ментами KvjvJt’

*)•

 

 

 

 

 

 

Формулы для

корреляционных

функций

(^* ^2) легко

получаются

аналогично

тому, «ак

были выведены

формулы

для системы с одной степенью свободы. Решение

уравнений

(5.34) представляется

через

импульсные переходные

функции

hk {t—т) следующим образом:

 

 

 

 

 

*) Q* W dx.

Подставляя это выражение в формулу (5.38), получим:

и

и

K fjfk (к, ( ,) =

— xj)hk (tt — TJ)KQ/ К та) rfxjdT,. (5.41)

J’

f

00—ôo

 

Если внешняя нагрузка является стационарной случайной функцией времени, то, вводя спектральные плотности нагрузки <pQjQk (to) и спектральные плотности обобщенных координат

Ф///Л (со), получим

= T J ф/ ; Л «

<542)

— 00

 

где

 

Ф/1fk (ш) = Fj( — *“) Fk

®QjQk H -

 

(5.43)

Нетрудно показать,

что

если

корреляционные

функции

KQJQJ Ï) являются

действительными

функциями т, то

корреля­

ционные функции

Kfffk (т),

определяемые

по формулам

(5.42)

и (5.43), тоже являются

действительными

функциями

т. При

этом

 

 

 

 

 

 

 

 

 

 

 

( - ’)•

 

 

(5 44 )

Как уже указывалось, наибольший интерес представляют значения корреляционных функций при t\ = t2=t и, в частности, при стационарном случайном процессе — величины К/^к (0).

Учитывая, что

 

Шу— «2 4"2il3yù)yù>

легко найдем:

ф,QjQk H d(i

К

(5.45)

V/ fk

(“/ — “2— 2/pytfly ш) (wA— ш2 + 2i'$ku ku>)

J »

53. Условия, при которых взаимной корреляцией обобщенных координат можно пренебречь

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

то все Ф0 .Q (ш) = 0 при } ф k Формула (5.45) принимает V/ я

вид

Kfkfk

(“>А— ш3)2+

й>)2

(5.46)

 

K/jfb (0) = 0 (у ф k).

При этом вычисления допускают значительные упрощения. Формулы типа (5.46) были положены в основу работ (70, 118, 119] и др.; при этом анализ условий, при которых допустимо пренебрежение взаимной корреляцией, оставлен без внимания.

То, что спектральные плотности Ф///Л (о) при H=k, вообще

Соседние файлы в папке книги