книги / Статистические методы в строительной механике
..pdfРассмотрим вначале одну случайную функцию -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, вообще