книги / Статистические методы в строительной механике
..pdffl( 0 = 2 ^ {qi' q" • |
• qr>Ъ (fr+ь Çr+2) |
. qm, /).[(7.11) |
Iх |
|
|
На рис. 119 показано, как гори помощи двух членов ряда мо
жет быть |
описано |
усиление |
высокочастотных |
'составляющих |
||||||||
спектра |
на |
заключительной |
стадии |
землетрясения. |
Здесь |
|||||||
<Di(o>) |
и Ф2(ш) —г спектральные |
плотности |
для |
функций |
<fr(0 |
|||||||
и с?2(0 |
в формуле (7.11). |
|
|
|
|
|
|
|
||||
Вычислим корреляционные функции для обобщенных ко |
||||||||||||
ординат fit (0 . Принимая, что |
|
|
|
|
|
|
||||||
ак(0 = |
2 |
|
q2’ |
’ qn ^ |
(^+ЬР9г+2) |
. qm>t), |
||||||
и используя формулу'(7.4), найдем |
|
|
и и |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
\ л . |
|
|
|
|
|
« - < - » • |
2 |
|
И |
|
|
|
|
|
|
|
|
|
|
|
!*■»• N..... |*у 0. 0 |
|
|
||
|
|
х |
• • • \ |
s t a — xs) Л . .p-i ta ) — |
|
^ ta) x |
|
|||||
|
|
Pl.P-l |
«Pï. N |
...»Ps . V-s ta. ъ |
**)<ta*i |
• dxs' |
||||||
Здесь JK |
|
|
|
— кор |
|
|
|
|
|
|
||
реляционные |
функции |
для |
|
|
|
|
|
|
||||
<fv (*)• |
|
|
|
|
|
|
|
|
|
|
|
|
Поскольку ускорения ak (/) |
|
|
|
|
|
|
||||||
имеют распределение, |
практи |
|
|
|
|
|
|
|||||
чески являющееся |
симметрич |
|
|
|
|
|
|
|||||
ным, то средние значения |
fh{t) |
|
|
|
|
|
|
|||||
и вообще |
все корреляционные |
|
|
|
|
|
|
|||||
функции с |
нечетным |
s |
могут |
|
|
|
|
|
|
|||
быть положены равными нулю. |
|
|
|
|
|
|
||||||
При этом |
сильно |
возрастает |
|
|
|
|
|
|
||||
роль корреляционных |
функций |
|
|
|
|
|
|
|||||
второго порядка: |
|
|
|
|
|
|
|
|
|
|||
|
Kftfkta. *2) = |
|
|
|
|
|
|
|
|
|||
tx и |
hj ta — "Ohk(*a — x2) Aj ta ) Akta ) K |
|
t a — *2) à'i dx2. (7.12) |
|||||||||
|
J |
|||||||||||
|
|
|
|
|
|
|
|
J |
|
* |
|
|
Заметим, что корреляционные функции (7.12) могут быть найдены и другим путем — по методу канонических разложе ний, на что было указано в работе [21].
79. Пример
Рассмотрим задачу о поведении простейшей одномассовой системы при .сейсмическом воздействии, которое описывается выражением
|
а(НлА |
. |
(7.13) |
||
Здесь А |
и с — (параметры, |
|
если(<0 |
|
|
характеризующие (величину макси |
|||||
мального |
ускорения |
и |
продолжительность |
землетрясения, |
|
q>(/)— стационарная |
случайная |
функция. Для |
автокорреля |
ционной функции второго порядка К 99 (ti—т^) — К W(T) возьмем (выражение
К99(х) = К0е~ “1т 1cos 6 т, |
(7.14) |
где Ко — некоторая константа, которая выбирается так, чтобы выполнялось условие (7.10), б — 'преобладающая частота зем летрясения, а — величина, характеризующая степень корреля ции. Чем меньше значение а, тем сильнее корреляция, т. е. тем ;£же диапазон частот, в котором заключена большая часть энер-
гии сейсмического воздействия. Последнее видно из формулы для спектральной плотности энергии
а также из ,рис. 120.
Обработка акселерограмм землетрясений дает следующие пределы изменения параметров: с=0,3-*-0,7 \}сек, 0 = 10-î- -î-40 1/сек.
Наибольший интерес представляют значения корреляцион ной функции t), поскольку они равны среднему квадра ту обобщенной координаты конструкции при сейсмическом -воз действии. Выведем формулы для вычисления #//(** 0- Подста вим выражения (7.13) в формулу (7.12) при j —k и U—tz (ин декс k всюду опущен). Учитывая, что
получим
tt
лr»
j Jexp IPK + XJS) — <х|тд — TjJJX
X sin Q (h — “ÎI) sin 2 (t2— т2) cos 0(xi — x2) dxt dx2. (7.15)
Здесь p=ie— с. Формулу (7.15) можно также представить после элементарных преобразований в виде
где
где
Форм-улу (7Л8) .можно также записать в виде
(7.20)
где
АУКо
(7.21)
Q Y2Q | р |
Чтобы оценить влияние расположения максимума спектраль ной плотности на поведение среднего квадрата обобщенной ко ординаты f, были произведены вычисления для трех случаев: п= 1,1, п=1,2 и л=1,3. При этом было принято, что а='2,0 1/сек (следует отметить, что, по предварительным данным, реальные акселерограммы имеют более слабую корреляцию). Соответ ствующие кривые, сглаженные относительно малых колебаний высокой частоты, показаны на рис. 121. Приближение «преоб ладающей» частоты 6 к собственной частоте системы й, как и следовало ожидать, вызывает увеличение максимальных значе ний среднего квадрата обобщенного перемещения. Вместе с тем, максимальные значения достигаются в более поздние мо менты времени и, чтоособенно важно, удерживаются в течение более продолжительного времени. Последнее обстоятельство имеет, как мы увидим позднее, существенное значение для оцен ки вероятности разрушения при сейсмическом воздействии..
На рис. 122 представлены кривые для трех различных значе ний константы а, характеризующей корреляцию ускорений в различные моменты времени. При этом взято постоянное отно шение частотn=il,il.
80. Случай многокомпонентного сейсмического воздействия
Будем называть сейсмичеокое воздействие многокомпонент ным, если оно описывается в общем случае тремя линейными смещениями и тремя вращениями основания по отношению к трем ортогональным осям. В этом разделе будет показана воз можность приложения развитой выше теории к многокомпо нентным сейсмическим воздействиям [135].
Рассмотрим сооружение, отнесенное к декартовой системе координат хи х2, х3 (рис. 123). Допустим, что сооружение яв ляется линейно-упругим телом. Обозначим компоненты объем
ных сил через Xj(xu х г , Х з ) , а компоненты вектора упругих .пе ремещений— через Uj- (хи х2, x3t t). Пусть перемещения точек тела под действием статических внешних сил определяются из системы трех линейных дифференциальных уравнений
Lj (ult и2, и3) = Xj |
(7.22) |
( /= 1., 2, 3).
Явную форму этих уравнений здесь не выписываем.
Предположим, что основание .сооружения является абсолют но твердым телом, движение которого задано в стохастическом смысле. Свяжем систему координат с основанием сооружения. Пусть a j(t) — линейные ускорения начала координат, $j{t) — угловые ускорения при вращении относительно координатных осей (/='!, 2, 3). Пренебрегая перемещениями Uj по сравнению с характерными размерами сооружения, найдем переносные ус корения a*j (хи Хз, хз, t) в произвольной точке:
O l = + Х 3 р2 |
* 2 Ре» |
|
Û2= «2+ Х1Рз-ХЭPl» |
(7.23) |
|
Дз= «з+ *2 Pi— |
xiРг- |
|
Уравнения движения могут быть получены из уравнений (7.23) подстановкой сил инерции .согласно принципу Далайнора
/ * |
3*и/ \ |
(7.24) |
Xj = — p у а/ + |
-gp-J |
|
U = 1. 2, 3). |
(722) |
|
Здесь р — .массовая 'плотность. |
Решение уравнений |
с правыми частями, задаваемыми согласно (7.23) и (7.24), оря однородных граничных условиях ищем в виде
i> ^2»%3, 0 = |
fk (О *?jkC^i» |
*^2J |
(7.25) |
(/ = |
1*. 2, 3). |
|
|
Здесь fk (t) — обобщенные координаты, |
<руЛ |
(х\, х2, х3) — |
«тройки» собственных функций для задачи о собственных коле баниях при неподвижном основании. Эти собственные функции удовлетворяют уравнениям
L j (<р1А, <р2А» <Рзк ) — 9 Ш1 «Ру* = 0 |
(7.26) |
(/ = 1, 2, 3), однородным граничным условиям на поверхности тела и усло виям ортогональности
J JJ Р (<?и 9и + |
|
Ъг + <Рз* <?з/) <W = |
0, (k Ф /), |
(7.27) |
|||||
где V — объем упруго,по тела. В уравнениях |
(7.26), как и ранее, |
||||||||
через сùh обозначена одна |
|
из частот собственных |
колебаний. |
||||||
Полнота системы собственных функций |
|
и сходимость разло |
|||||||
жения (7.25) следуют из известных общих теорем. |
|
будем счи |
|||||||
В дальнейшем число членов в разложении |
(7.25) |
||||||||
тать конечным и равным |
п. Тогда |
решение |
в |
форме |
(725) |
||||
должно трактоваться как .приближенное. |
|
с правыми |
частя |
||||||
Подставим ряд (7.25) |
в уравнения |
(7.22) |
|||||||
ми в форме (7.23) и (7.24). Принимая |
во внимание соотноше |
||||||||
ния (7.26) и (7.27), получим систему |
уравнений в |
обобщенных |
|||||||
координатах |
|
|
|
|
|
|
|
|
|
/ |
* + « » / . - в . |
|
|
|
|
|
(7-28) |
||
( 4 = 1 , 2, . . . л ) . |
|
|
|
|
|
|
|||
Обобщенные силы Qk(t) |
определяются .по формуле |
|
|
||||||
|
|
Qk (t) = |
|
|
|
|
|
|
|
J J j* Р [(«Ч+*. Ра—*гРз) 9ik + |
(аа+*хРз—*з Рх) ср2*Ч^(аз—Н-*2pi— |
dV |
î î l P ( + Т2* + |
dV |
Эта формула может быть записана в форме
(7.29)
о»Ю “ - 2 * / » « у ( о - 2 х'* р'<о. /=х /“1
где
*Jh = I l J
P ( *?îft + <?2ft + “fft ) dV
(7.30)
Ш |
P ^ Z + l Ÿ/+2. ft — -*/+2 ^/+ 1 . fe) |
|
bfk =■ |
IJ I P ( ?lft + ?!* + Ûk ) dV |
(<р« =?u . ы ■ |
|
|
|
Рассмотрим теперь пространственную раму или другую про |
||
странственную систему, состоящую |
из стержневых элементов. |
Чтобы упорядочить множество точек системы, выберем путь об хода элементов и введем координату s, измеряемую вдоль оси элементов от некоторой начальной точки. Вместо собственных
функций <p/ft(*i, Х2,*з) |
получим |
собственные |
функции Фу* (s). |
|||||
Пусть т (s)— погонная |
масса |
|
стержня, М-,— массы, сосредо |
|||||
точенные в точках |
системы |
с координатами |
s= sv Формулы |
|||||
(7.30) |
могут быть переписаны следующим образам: |
|||||||
|
Jm(s) Vk ds +2 Мт Ц |
|
|
|||||
,к J |
т (s) ( ?ik + |
?2k+ ?!fe) * + |
2 |
MÏ ( |
|
I,., |
||
|
|
|
|
|
7 |
|
|
T |
4k= |
Jm (s) (#y+1 9f+2,k |
|
xj+2 fj+i.k ) |
(7.31) |
||||
|
|
|
2 |
MT( Tl* + |
fift + |
|||
J «*(*) ( Ÿ1A 4- ?!fc+ ?3ft) ds+ |
Ÿ3ft)| |
|||||||
|
|
|
|
|
T |
|
|
|
|
2 |
(*/+1 Ÿ/+2,ft— |
xi+2 Ÿ/+1.A ) I J„J |
|
||||
+ |
|
|
|
|
|
|
|
|
Jw (s) ( ??*+ |
+ Фза) ds + 2 |
Tf |
M i ( <fXA+ |
+ |
)| s_s |
|||
|
|
|
|
|
|
|
ï |
Силы внутреннего и внешнего трения до сих пор не учиты вались. Эти силы, как обычно, можно формально учесть путем добавления в уравнения (7.28) членов, содержащих обобщен ные скорости. Пренебрегая диссипативными связями между формами колебаний, придем к рассмотренным ранее уравнени ям в форме (7.2). При этом вместо я*(0 в правые части под ставляются обобщенные силы, определяемые по формулам
(7.29), (7.30) и (7.31). Дальнейшие выкладки остаются без су щественных изменений. Так, для получения корреляционной функции обобщенных координат имеем формулу
_________ it ft
МММ**) = J JA yft-M M **— с*) Qj('i)Qb{4) d4 d.Xb(7.32)
|
О о |
|
|
|
|
|
|
|
|
|
( М |
= |
1,2, |
.л). |
|
||
Согласно формуле |
(7.29), корреляционная функция |
обоб |
||||||
щенных сил |
Qj(ti)Qk (М |
может быть |
выражена через корре |
|||||
ляционные функции линейных и угловых ускорений |
|
|||||||
Q, |
|
|
|
3 |
3 |
|
|
|
й ) |
w |
= |
2 |
2 |
[»«/* -■ .№ > » .(« + |
|
||
|
|
|
|
Ш—1Л=1 |
|
|
||
+ * m j K k |
«m (*l) |
Ря(М + |
Xmy*л* |
(**) Рл(М + |
|
|||
|
|
+ Х„,Ч«М Щ Ж )]- |
(7-33) |
Способ апециализации сейсмического воздействия в форме (7.7) легко может быть распространен на многокомпонентные воздействия:
®у(0 = А/(*71>Цъ |
• ?г>0?у (Яг+v Qr+v ' |
' |
I |
|
Ру (О = |
{QI, QÎ, |
. (jri t) фу(7г+^>Çr+V |
• Qs*t)‘ | |
|
Здесь Лy (0 |
и ДДО — детерминированные функции, <р,- (0 и |
фу (^)— стационарные случайные функции. Простейшая форма аат1роксима1Ц'ии имеет вид
A J = A f ) e - cj t ; 9 .(t,)ь (t2) = К}e*J1,1 " ,г 1cosду(tг- /2), (7.35)
где Л|0), Су, /Су, ау и 0у— случайные числа ((параметры) земле
трясения. Аналогичное представление может быть введено для угловых уокорений р/(<)• Подсгавим формулы (7.33), (7.34) » (7.35) в формулу (7.32). Тогда интеграл, стоящий в правой час ти, может быть выражен через элементарные функции. Даль нейшие вычисления производятся методом, развитым в п. 79.
81. Оценка вероятности разрушения и ожидаемого срока
службы
Основной задачей теории являются вычисление ожидаемого срока службы сооружения и разработка методов проектирова ния таких сооружений, вероятность разрушения которых в те чение установленного срока службы будет не больше установ ленного значения. Определение системы корреляционных функ ций служит лишь первым шагом на пути решения проблемы.
Правда, их знание позволяет легко найти закон изменения во времени математического ожидания и среднего квадрэтического значения обобщенных координат; однако Для практических це лей этого оказывается недостаточно.
•Вторым шагом по пути решения проблемы является опре деление плотности вероятности обобщенных координат в фик сированный момент времени и.яри фиксированных значениях случайных параметров Ц\, Яг, - •. Ят• Эту плотность вероятности будем обозначать через p(fu /2, ...fjt) . Чтобы найти эту плот ность вероятности, нужно знать полную систему моментов для
fi» Î2, fn• |
Эти моменты равны корреляционным функциям |
.... fps |
(iu к, . . . t s) при tl = k = - ~ = t s= t - Дело оводит- |
ся к решению некоторой обобщенной задачи теории моментов. Основываясь на свойствах землетрясений, можно внести не которые упрощения. Как уж указывалось, компоненты ускоре ния трунта, особенно горизонтальные компоненты, имеют рас пределение, близкое к симметричному, — оба направления прак тически равновероятны. Поэтому математическое ожидание ус корения, • а также все моменты нечетного порядка можно при нять равными нулю. »Пр« этом важнейшее значение приобрета ют моменты второго порядка и соответствующие корреляцион ные функции. Для случая, когда плотность вероятности подчи няется нормальному закону, достаточно найти лишь моменты второго порядка. Вопрос о том, насколько истинное распределе ние приближается к нормальному, требует дополнительного ис следования. Известно, что нормальное распределение имеет место в тех случаях, когда рассматриваемая величина форми руется под влиянием большого количества статистически неза висимых случайных причин. Возможно, что такую ситуацию мы имеем при сильных землетрясениях в пунктах, не слишком близких к эпицентру. Если даже распределение и не является нормальным, оно, по-видимому, все же может быть приближен но принято таковым для оценки не ’слишком малых вероятно
стей.
Важнейшее свойство нормального распределения состоит в том, что оно остается нормальным после линейного преобразо вания случайных величин (п. 11). Поэтому для линейно-упру гих конструкций .при малых перемещениях из нормального распределения обобщенных ускорений вытекает нормальное распределение обобщенных координат, а также любых их ли нейных функций (напряжений, смещений и т. п.). Для плот ности вероятности p(fu fi, fn\ к имеем формулу (5.3)
P ( f i , h , |
f n \ t ) = |