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

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

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

fl( 0 = 2 ^ {qi' q"

qr>Ъ (fr+ь Çr+2)

. qm, /).[(7.11)

 

 

На рис. 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

л

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 ) =

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