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

книги из ГПНТБ / Д’Анжело, Г. Линейные системы с переменными параметрами. Анализ и синтез-1

.pdf
Скачиваний:
9
Добавлен:
19.10.2023
Размер:
6.9 Mб
Скачать

Размерность частотной характеристики для ускорений ФѴ) (ш >

будет равна нулю. Точно такие же размерности частотных ха­ рактеристик получатся и из формулы (III.75), что вытекает из структуры этой формулы. Те же размерности имеют передаточ­ ные функции Ф (s). Из выражения (III.77) следует, что весовые функции перемещений будут иметь размерность Т, а ускоре­

ний — Т~1,

что согласуется с результатами, изложенными в § 5-

настоящей

главы.

 

 

 

 

 

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

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

формулы

(III.

77)

или

(III. 74)

можно преобразовать к действительной

форме.

 

Если

система

описывается

уравнениями (III.

64), матрица

которых

состоит

из

дифференциальных операторов второго порядка

с

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

являются дробно-рациональными

функциями аргумента (со.

Поэтому они могут быть представлены в виде

 

Фй(м>) =

И

+ ІФ« (to),

(HI.81)

где Ф і— четная, а Ф2— нечетная

функция аргумента со.

Урав­

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

этому условию и удовлетворительно описывают колебания

реаль­

ных сооружений, поэтому будем считать представление (III. 81) справедливым по отношению к аналитическим и эксперименталь­ ным частотным характеристикам. Подставляя (III. 81) в (III. 75) и (III. 77), приходим к выражению весовой функции в действи­ тельной форме:

 

 

 

 

 

со

 

 

 

 

 

 

 

k(t) = -^-

I [Ф16 (ш) cos ut — Ф2Й(ш) sin <о£] du.

 

(III.82)

 

 

 

 

 

о

 

 

 

 

 

Возможность представления частотной характеристики

в

ви­

де

(III.

81)

с

Фі—четной и Ф2— нечетной функциями

со

нужно

рассматривать

как

необходимое условие существования весовой

функции. Если бы указанные условия четности

не соблюдались»

функция h(t)

имела

бы комплексные коэффициенты,

что

лишено

физического

смысла.

 

 

реакции и

 

Найдем действительную форму h(t) через спектры

воздействия. Вывод

сделаем независимо от формул

(III.

81)

и

(III. 82) для того, чтобы подчеркнуть физический смысл

усло­

вия

(III.

81).

на входе системы, заданное в

виде

 

ускорения

 

Воздействие

 

основания сооружения, и реакции на всех выходах суть действи­ тельные функции времени. Их .комплексные спектры имеют вид

со

(im) =

j w (t) e~M dt = lx(ш) - il2(о»),

(III.83)

0

 

где /j — косинус-преобразование Фурье функции w (t) и

 

 

00

 

/j

= § w (t) cos M t ,

(III.84)

 

<>

 

120

a Jj — синус-преобразование

со

/ 2 = J w (і) sin mtdt.

о

Аналогично, для спектра реакции напишем

Sx (fco) = (со) iJ2(ш),

где

Уі = f л: (£) cos соЛ#

о

D

У2 = J а:(£) sin cuA/£

(111.85)

(111.86)

(111.87)

Из выражений (III.84) и (III.85) или (III.86) вытекает следую­ щее свойство комплёксных спектров действительных функций:

 

 

 

 

 

S ( - m ) = S*(iv).

 

 

 

(III.88)

Частотную характеристику представим в виде

 

 

 

 

Ф (ш) =

Sw('“)

= Ct ( CD) -f- іС2(ш).

 

 

 

 

 

 

 

S x ( Іи:)

 

 

 

 

 

 

Возможность такого представления следует из (III.83) и

(III.86);

оно не содержит никаких предположений о

четности

функций

С, и

С ,

 

 

 

 

 

функции

теперь

можно

преобра­

Выражение для весовой

зовать:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*(*)= f

-- m(/M)

6Wüfco = f

$w (*'“)

е

ІШ , .

 

 

 

5 , И

acu -j-

 

 

 

 

J

 

5 Г(/Ш)

 

 

J

 

 

 

 

 

 

—OO

Л

 

 

 

0

 

 

 

 

 

 

OO

 

 

 

 

OO

 

 

 

 

 

 

п

~ай' ы)У

 

п

 

 

öfcu

 

 

J*

Г

 

( г’ш)

- і ш і J

 

Г , п

. . п ,

2<оі

 

 

+ J

 

 

 

 

<С“= )(С.+‘С,)«

 

 

 

ТО

 

 

 

 

 

СО

 

 

 

 

 

 

+

(С, — г'С,) e~lU3idw = 2

j(С, cos cu£ — С2 sin cozf)<іш,

'111.89)

где

 

 

 

 

 

 

АЛ + АА

 

 

 

 

 

 

 

 

 

Ct =

 

 

 

 

 

 

 

 

 

г2

I

г2

 

 

 

 

 

 

 

 

 

 

 

У 1

+

у 2

 

 

 

(III.90)

 

 

 

 

 

г

_V : + АЛ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 ""

Г)

 

Г)

 

 

 

 

Л + У£

Вычисления Л(У) по формулам (III. 89) и (III. 90) легко мо­ гут быть запрограммированы для ЭЦВМ.

121'

В технике анализа линейных систем для экспериментального определения весовых функций нередко используются соотно­ шения между корреляционными функциями на входе и выходе системы, называемые соотношениями Винера—Ли [46, 408]:

со

 

 

 

K r m =

1 щ

 

 

 

 

 

 

—оо

 

 

 

 

 

 

Kyy(t)=

I h { t - ^ ) K yxWdx

[f

(III.91)

 

 

 

или

—со

 

 

 

 

 

 

оо

 

 

 

 

 

 

 

 

 

 

где

?h/i — корреляционная

функция весовой функции, для кото­

рой

справедливы

равенства

 

 

 

 

СО=

оо

А(О А +

 

оэ

 

 

 

J

-с) dt

= JL j I Ф (гео)I2 Л

» .

(III.92)

 

 

—ОО

 

—ОО

 

 

 

Построение

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

функций для воздействий

взры­

вов и слабых землетрясений, для которых в основном и предназ­ начены излагаемые здесь методы, затрудняется тем, что записи имеют малую продолжительность, в силу чего корреляционные функции можно построить только для малых значений t. Можно, однако, определить весовые функции по формулам (III. 91),если

по численным значениям корреляционных

функций, определен­

ных для

малых t, найти их

аналитические

аппроксимирующие вы­

ражения,

аналогичные

тем,

которые

использованы в § 2 гл. I.

Тогда в

интегральных

уравнениях

(III. 91) будут известны точ­

ные значения корреляционных функций в левой части под знаком интеграла и для численного определения ординат h(t) формаль­ но окажется возможным применение метода последовательного исключения неизвестных из уравнений (III. 51) с треугольной матрицей. Полезность получаемых оценок функций h(t) должна быть дополнительно исследована.

Изучение корреляционных функций на входе и выходе сис­ тем инженерно-сейсмометрической службы представляет интерес

и в других отношениях, например,

для исследования

вероят­

ностных характеристик

слабых землетрясений.

Поэтому

вычис­

ление функций Kxx{t),

К ху (t) и Kyy(t)

должно

быть включено в

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

•122

§ 8. С тандартны е расчеты на сейсмические воздействия

Согласно строительным нормам и правилам, принятым в

СССР, а также зарубежным строительным кодам, сейсмостойкие сооружения должны удовлетворять расчету на условные статиче­ ские нагрузки, заменяющие действие землетрясений. Анализ сей­ смических движений основания сооружения, так же как и иссле­ дование динамических свойств сооружений, показывает, что ве­ личины расчетных нормативных нагрузок и характер их воздей­ ствия не эквивалентны действию землетрясений и не отражают динамических свойств сейсмического процесса. В главе I было показано, что нормативная спектральная кривая, принятая в СНиП ІІ-А 12-69, не соответствует вероятностным спектрам уско­ рений грунта. Тем же недостатком страдают и спектры, принятые в заграничных кодах. Было показано также, что путем уточнения спектральных характеристик на основе изучения статистических свойств акселерограмм нельзя получить хорошего соответствия расчетных сейсмических воздействий действительным.

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

Периодограммы сейсмических

ускорений,

приведенные в главе

I (рис. 15—19), показывают, что

спектры

отдельных акселеро­

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

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

На рис.

33

приведен спектр ускорений акселерограммы 8-8

Г-40. При изменении

периода первой формы собственных колеба­

ний здания

на

10%

сейсмическое воздействие изменится почти

в два раза.

 

 

 

При расчете по спектральным кривым применяются различные методы сглаживания спектров, к числу которых относится и пост­ роение вероятностных спектров. Эти методы не применимы к от­ дельным акселерограммам, так как сглаженные акселерограммы не имели бы никакого смысла.

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

нимается

в качестве методики, исключающей

случайные ошибки.

 

 

 

 

 

На самом

деле

такое

 

решение

 

 

 

 

 

проблемы

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

 

 

 

 

 

вильно

и не исключает

 

фактора

 

 

 

 

 

случайности.

Вполне

возможно,

 

 

 

 

 

что некоторая

частота, имеющая

 

 

 

 

 

существенное

значение для дан­

 

 

 

 

 

ной

модели, на всех выбранных

 

 

 

 

 

для расчета акселерограммах ока­

 

 

 

 

 

жется

вблизи пика или

провала

 

 

 

 

 

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

 

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

Рис. 33.

Расположение

периода

Поэтому исключения фактора слу­

первой

формы колебаний

соору­

чайности

необходимо

достигать

жения

относительно

спектра ре­

другими

средствами и для

каж­

 

 

акции.

 

дой

акселерограммы

в

отдель­

 

 

 

 

 

ности.

 

 

 

 

 

 

 

Одной из возможностей является построение моделей с пере­

менным

спектром.

 

модель

сооружения примем столбец

ве­

За математическую

совых функций hk (t ) в соответствии с выражением (III.5),

неза­

висимо

от того, каким

способом — теоретически

или

 

экспери­

ментально — они

определены.

Установим

важное

 

свойство

весовых

функций

ускорений, заключающееся

в

том,

что

если

частоты всех форм колебаний умножить на один и тот же коэф­

фициент V,

то весовая функция для новой

системы

получится

из исходной путем умножения функции и ее аргумента

на тот

же коэффициент ѵ.

Обозначим

весовую

функцию

ускорений

системы,

имеющей

некоторый

спектр

собственных

частот

ші, ш 2......, через h(t).

Весовую функцию другой системы,

имею­

щей спектр частот

ѵш2.......

обозначим

Лѵ(£). Тогда

будет

справедливо соотношение

 

 

 

 

 

 

А (t) = vA (v^).

 

 

(III.93)

Это выражение имеет и самостоятельное значение. Оно может

быть использовано при вычислении весовых

функций

или

при пе­

рерасчетах, связанных с изменением параметров жесткости или на­ грузок сооружения.

Доказательство следует из формулы (III. 10). Под знаками экспоненциальных и тригонометрических функций находится про­

124

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

Для весовых функций перемещений, согласно формуле (III.9), соотношение между Ли (t) и (t) будет

h it) — h {'it).

Обозначив весовые функции буквами, принятыми для обоз­ начения физических величин, к которым они относятся, со знач­ ком 0 наверху, напишем следующую систему соотношений:

перемещения

«V ( 0 = 7 “ М

скорости

п ѵi t ) = ѵ ( ч і )

ускорения

w w(t) = YW (yt)

I •

(III.94)

нагрузки

sv(0 = vs(v0 перерезывающие силы

Весовые функции (III.94) зависят от параметра ѵ, следователь­ но, от него будет зависеть и реакция, определяемая по форму­ ле (III.6):

X (ѵ, t) = f w0 (t) Av (t — Д di.

(III.95)

Максимальное абсолютное значение реакции

при воздействии

данной акселерограммы, соответствующее некоторому моменту времени t — t(v), обозначим х(ѵ);

t

 

X (v) = max^ |® '0(т)//ѵ

(III.96)

Изменяя значения параметра ѵ в пределах

от Vj до ѵ2

ѵі ѵ < • V,,

 

можем найти среднее по параметру значение максимума реакции

123

(III.97)

Формулой (III.96) можно воспользоваться для „сглаживания“ результатов расчета на воздействие акселерограммы. Приняв

 

V,

 

(ІІІ.97), получим результат,

и произведя расчет по

формуле

 

< 1 <

ѵ2

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

системы относительно пиков и провалов спектра

акселерограм­

мы. Для акселерограмм,

рассмотренных в главе

I, достаточно

принять смещение спектра в ту или другую сторону на 20%, т. е.

'^ = 0 ,8 ; ѵ2=1,2.

Выражения абсолютных

 

максимумов

реакции

будут иметь следующий

вид:

 

 

 

для

перемещений

 

 

 

 

 

 

для

скоростей

 

 

V

оJ

 

 

 

 

п(ѵ) =

т а х /

t

 

 

d-z ;

(III.99)

 

[щ0(т)ц° Гѵ(£ — -)]

 

 

 

 

b

1

1

 

 

для

ускорений,

нагрузок

и перерезывающих сил

 

 

x(v) =

max; V

t

«'о (т) h [v ( t — x)]of- .

(III.100)

 

J

Расчет сооружений с переменным спектром можно обосновать и с другой точки зрения. Известно, что жесткость сооружения не

остается постоянной за время его существования. Для

железо­

бетонных

конструкций в

начальной

стадии

эксплуатации

жесткость

может

возрастать, в

дальнейшем

при

накопле­

нии

деформаций,

возникающих

под

воздействием

различ­

ных факторов,—снова уменьшаться. Доказано [124], что

в про­

цессе

колебаний сооружения

жесткость

не остается постоянной,

существенно уменьшаясь в режиме больших амплитуд. Предло­ женный метод расчета с переменным спектром косвенно отража­ ет эти особенности поведения сооружений. Изменение частот на 20% соответствует изменению жесткости приблизительно на 40%.

При дискретном шаге интегрирования формулу (III.97) можно записать в виде

(1II.101)

где ѵт — значение ѵ в точке т, s — число участков подразделения отрезка [ѵь ѵ2], x(vm) — максимум модуля реакции, соответствую­ щий значению параметра ѵ в m-й точке подразделения. Для вычис­

126

ления интегралов (III.95) может быть использован алгоритм (111.50} с добавлением операции сдвига индекса k—/ путем умножения его на V и умножения или деления всей суммы «а ѵ, если это требу­

ется в соответствии іс формула;ми. Так, для

вычисления ускорении

будем иметь:

 

 

 

 

ft

 

 

Xmk = * ( Vm *ft ) =

2 Wj ^

{k_}) .

(HI. 102)

 

J=1

 

 

Усложнение алгоритма, связанное с умножением

индекса при

h на величину ѵт , заключается

в том, что

произведение '>m(k—j)

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

h (t) = hk , { k - \

 

 

где Д^— шаг интегрирования.

 

 

При достаточно малом шаге интегрирования это упрощение не

внесет существенных искажений. Можно несколько

выиграть в

точности, принимая ступенчатую функцию по формуле

k ( t ) = hk~12

hk , { k -

і ) д г < г < м / .

 

Применяя указанные

алгоритмы,

можно составить

програм­

мы для вычисления при помощи ЭЦВМ средних реакций систем с переменными спектрами «а воздействие акселерограмм.

Изложенный метод позволяет исключить один из факторов случайности, с которым 'связан расчет на воздействие акселеро­ грамм— фактор случайного расположения -спектра модели отно­ сительно опектра акселерограммы. Другой такой фактор—случай­ ный характер самой акселерограммы. Широко известный метод ос­ реднения результатов расчета по нескольким акселерограммам пред­ назначен для оглажив.ания случайных свойств каждой отдель­ ной акселерограммы. Однако для того чтобы этот прием достигал цели, необходимо соблюдение некоторых условий.

Первое условие заключается в том, что выборка акселеро­ грамм должна быть -представительной. В большинстве случаев для этого -стараются группировать акселерограммы по интенсив­ ности землетрясения, отбирая отдельно семи- и отдельно — вось- мибалльны-е акселерограммы. С точки зрения исследований, ре­ зультаты которых изложены в главе I, этот прием не достигает цели. Выборка может считаться представительной только в том случае, если она содержит акселерограммы достаточно разнооб­ разного спектрального состава (конечно, частотный диапазон мо-

127

жет быть согласован с региональными характеристиками возмож­ ных землетрясений, если таковые известны). Этому условию удовлетворяют пять акселерограмм, которые 'подробно исследо­ ваны в главе I. Покажем периоды максимумов периодограмм Та и спектров реакции Т г этих акселерограмм:

 

 

Тп

Тх

8-5

Г-20

0,13 -0,22 —0.32

0,13

8-6

Г-69

0,26

0,26

8-3

Г-52

0,48

0,30

8-8

Г-10

0,66

0,21—0,44

7-25

Г-40

1,55

1,55

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

диапазона: 8-5 Г-20; 8-3 Г-52 и 7-25 Г-40. Если

ответственность

задачи превалирует над вопросами стоимости

машинного

вре­

мени, в выборку из пяти акселерограмм можно

добавить

еще

другие, например 7-16 Г-50, что расширяет ее частотный диапазон. Второе условие заключается в сопоставимости элементов вы­ борки. Подбор акселерограмм по интенсивности не решает этой задачи, так как оценка интенсивности не определяет механичес­ ких и статистических свойств сейсмического движения. Так, де­

вять восьмибалльных акселерограмм в табл.

1 имеют дисперсии,

изменяющиеся от 7200 до 800 сж/се/с4, т. е.

в девять раз. Срав­

нимость акселерограмм должна достигаться

их нормированием

по стандарту. При определении нормирующего множителя можно

руководствоваться оценкой балльности,

приведенной

в

§ 1 гла­

вы

I.

 

 

 

 

 

 

 

Для выборки в п акселерограмм с номерами г= 1, 2,...п рас­

четная формула сейсмического воздействия будет иметь вид

 

жср-р

 

 

 

 

(ШЛОЗ')

где

л'ср — воздействие

r -й

акселерограммы,

определенное со сгла­

 

живанием по спектру собственных частот (формула

 

(111.98);

 

 

 

 

 

 

 

ср — расчетное значение стандарта;

 

 

 

 

/ — региональная

вероятность

землетрясений

с

частот­

 

ным составом,

соответствующим

r -й акселерограмме:

 

П

 

 

 

 

 

 

Если отсутствуют данные о распределении вероятностей зем­ летрясений по частотному составу, можно производить расчет по наибольшим значениям реакции*

*ср.р = т а х г Т ~ Ч р .

(Ш.103")

 

128

Де rx ^ — величина, определенная со сдвигом спектра для г-й акселерограммы по формуле (Ш. 98); ар — расчетное значение

стандарта.

Можно .полагать, что расчет на воздействие нескольких аксе­ лерограмм с двойным осреднением—іпо спектру модели и по со­ вокупности акселерограмм—будет более точным методом прог­ нозирования возможных сейсмических воздействий, чем расчеты по спектральным кривым.

§ 9. Сейсмические реакции двумерных систем

 

Приведенные в предыдущем параграфе

формулы решают

задачу о расчете одномерных систем на воздействие акселеро­ грамм. На основе обобщенного принципа суперпозиции (III. 19) можно построить решения более сложных задач. Рассмотрим воз­ действие по закону акселерограммы на двумерные системы при синхронном движении основания и с учетом конечной скорости распространения сейсмических возмущений по продольной оси здания. Последняя задача типична для проблемы сейсмостойко­ сти протяженных систем.

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

 

 

 

/ге,л-,у(0 =

/г(*- А'- У! *)•

(III.104)

Ось Л' направлена по

горизонтали,

ось

К — по

вертикали, I—

абсцисса точки

приложения импульса в основании системы, л-и

у —координаты

точки

выхода.

 

как ипульс ускорения, со­

Входной

импульс

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

общающий

единичную скорость

бесконечно малому отрезку оси

основания d\ . Реакция на единичный импульс,

распределенный

по всей длине,

будет

 

 

 

 

 

 

 

Ц х ,у ;

X, у;

t ) d \ .

(III.105)

Это выражение в задачах сейсмостойкости может быть принято

за весовую функцию сейсмической реакции двумерной

системы.

Реакция на воздействие акселерограммы запишется как

 

t

_

— x) dx ==

 

R {х, у; t) = j w 0(x) h (X, y \ t

 

0

 

 

 

t

 

 

 

= ^dk |гіу0 (т) /г($, X, у; t

x) dx.

(III.106)

L 0

 

 

 

9—248

129

Соседние файлы в папке книги из ГПНТБ