ОРДЕНА ТРУДОВОГО КРАСНОГО ЗНАМЕНИ
ЦЕНТРАЛЬНЫЙ НАУЧНО-ИССЛЕДОВАТЕЛЬСКИЙ
И ПРОЕКТНО-ЭКСПЕРИМЕНТАЛЬНЫЙ ИНСТИТУТ
КОМПЛЕКСНЫХ ПРОБЛЕМ СТРОИТЕЛЬНЫХ КОНСТРУКЦИЙ
И СООРУЖЕНИЙ им. В. А. КУЧЕРЕНКО
ГОССТРОЯ СССР

РЕКОМЕНДАЦИИ ПО ОПРЕДЕЛЕНИЮ
РАСЧЕТНОЙ СЕЙСМИЧЕСКОЙ
НАГРУЗКИ ДЛЯ СООРУЖЕНИЙ
С УЧЕТОМ ПРОСТРАНСТВЕННОГО
ХАРАКТЕРА ВОЗДЕЙСТВИЯ
И РАБОТЫ КОНСТРУКЦИЙ

Утверждены
Директором ЦНИИСК им. В. А. Кучеренко

11 марта 1988 г.

МОСКВА - 1989

Рекомендованы к изданию решением секции «Сейсмостойкость сооружений» Научно-технического совета ЦНИИСК им. В. А. Кучеренко Госстроя СССР.

Содержат алгоритм определения расчетной сейсмической нагрузки для сооружений с учетом пространственного характера внешнего воздействия и работы конструкции.

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

Предназначены для инженерно-технических работников проектных и научно-исследовательских институтов.

ПРЕДИСЛОВИЕ

Рекомендации составлены в развитие СНиП II-7-81 с целью практического использования методов расчета сооружений с учетом пространственного характера работы конструкций и сейсмического воздействия.

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

В Рекомендациях приведен алгоритм расчета по определению сейсмической нагрузки на сооружения с учетом пространственного характера воздействия и работы самой конструкции. Даны рекомендации по выбору расчетных динамических моделей сооружений и определению их параметров. Приводится алгоритм расчета режима свободных колебаний сооружений. Даны рекомендации по определению расчетных моделей сейсмического воздействия и выбору их численных параметров.

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

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

По методике настоящих Рекомендаций разработана автоматизированная система расчета сооружений на сейсмические воздействия «Полифем-S87».

Рекомендации разработаны в Отделе сейсмостойкости сооружений ЦНИИСК им. В. А. Кучеренко (основной текст - канд. техн. наук Ю. П. Назаров при научном руководстве д-ра техн. наук проф. Н. А. Николаенко, приложения 1, 2 - кандидаты техн. наук А. Т. Штоль и Г. В. Мамаева; приложение 3 - инж. А. Ш. Ревишвили).

ОБЩИЕ ПОЛОЖЕНИЯ

1.1. Настоящие Рекомендации составлены в развитие СНиП II-7-81 [17] с целью практического использования методов расчета сооружений с учетом пространственного характера работы конструкций и сейсмического воздействия.

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

Анализ регистрации движений грунта при землетрясениях показывает, что они являются сложными многомерными (пространственными) процессами. Сложность пространственного характера сейсмического воздействия состоит в том, что оно представляет собой поле движений грунта, определенное сейсмическими волнами с конечной скоростью их распространения. Сейсмические волны не являются монохроматическими, а представляют спектр одноименных волн различной длины, доминирующих в различные интервалы времени землетрясения. В начальные моменты времени приходят, имеющие большую скорость распространения, продольные Р-волны, затем - поперечные S-волны с меньшей скоростью распространения, но с большей интенсивностью воздействия. После S-фазы генерируются поверхностные L-волны (Рэлея-Лэмба и Лява), которые имеют еще меньшую скорость распространения. Наиболее сильные повреждения и разрушения сооружений наблюдаются при более интенсивных фазах движения грунта, определяемых спектрами S- и L-волн. Такое поле движений приводит к дилатации и ротации массива грунта в основании сооружений. Дилатация (изменение объема за счет его расширения-сжатия) массива грунта в основании сооружения определяется всеми видами волн P, S и L и приводит только к поступательным движениям основания сооружения. Эти движения происходят в трех направлениях. Ротация (изменение формы объема за счет его искажения) определяется только вихревыми составляющими поля S- и L-волн и приводит к вращению массива грунтового основания сооружения. Эти вращения происходят относительно трех осей.

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

1.3. В Рекомендациях рассматривается упругая стадия работы конструкций при колебаниях, соответствующих малым перемещениям и углам вращения масс сооружения.

1.4. Применение метода рекомендаций связано с обработкой громоздких массивов чисел. Использование рекомендуемых методов требует привлечения ЭВМ. Рекомендации составлены таким образом, что применение ЭВМ в проектной практике возможно в двух формах:

а) применение ЭВМ для выполнения отдельных этапов расчета:

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

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

- определение параметров сейсмического воздействия;

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

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

2. АЛГОРИТМ РАСЧЕТА ПО ОПРЕДЕЛЕНИЮ СЕЙСМИЧЕСКОЙ НАГРУЗКИ НА СООРУЖЕНИЯ С УЧЕТОМ ПРОСТРАНСТВЕННОГО ХАРАКТЕРА ВОЗДЕЙСТВИЯ И РАБОТЫ КОНСТРУКЦИЙ

2.1. В Рекомендациях приводится алгоритм определения сейсмической нагрузки на сооружения в наиболее простом варианте, т.е. без учета всех корреляционных связей при компонентной форме представления векторов сейсмического воздействия  и ^*). Этот алгоритм состоит из следующих разделов:

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

- формирование матриц жесткости или податливости и определение единичны усилий в элементах конструкций;

- расчет режима свободных колебаний;

- определение расчетных параметров векторов сейсмического воздействия  и ;

- расчет режима вынужденных колебаний;

- определение расчетных усилий в поперечных сечениях несущих конструкций сооружений.

^*) Структура алгоритма расчета сооружений как единых пространственных систем на сейсмические воздействия с учетом пространственного характера возмущения зависит от ряда факторов:

способа представления информации о векторах сейсмического воздействия;

учета корреляции между отдельными формами колебаний;

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

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

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

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

2.2. В основу расчета положены дискретные расчетные динамические модели, которые формируются исходя из конкретного типа сооружения и следующих положений:

- конструкции сооружения (перекрытия, отдельные конструктивные части и т.д.), жесткость которых на порядок больше жесткости других конструкций, считаются абсолютно твердыми;

- другие конструкции сооружения (элементы каркаса, ядра и диафрагмы жесткости, диафрагмы лестничных клеток и т.п.) считаются упругими связями с линейными диаграммами деформирования материала;

- в динамических моделях массы сооружения распределяются по абсолютно твердым телам, а упругие конструкции считаются невесомыми и их масса относится к соответствующим твердым элементам, т.е. динамическая модель сооружения представляется в виде отдельных (дискретных) абсолютно твердых тел, соединенных несущими невесомыми упругими элементами;

- фундамент сооружения в расчетной динамической модели принимается жестким и моделируется твердым телом.

2.3. Дискретно расположенные в динамической модели массы могут быть представлены:

- материальной точкой, которая может иметь максимум три степени свободы и характеризуется величиной массы m_k (k = 1, ..., n);

- твердым диском, имеющим три степени свободы и учитывающим геометрию распределения масс и упругих связей в плоскости диска. Диск движется в своей плоскости и характеризуется величиной массы m_k и центральным осевым моментом инерции θ_k относительно оси, перпендикулярной к плоскости диска;

- массой в виде твердого тела, имеющего шесть степеней свободы и учитывающего геометрические параметры распределения масс и упругих связей в пространстве. Инерционные параметры твердого материального тела характеризуются величиной массы m_k и тензором инерции масс относительно каких-либо осей [θ_k] (k = I, II, ..., n); в расчетах в качестве таких осей следует принимать главные и центральные оси инерции тела, для которых [θ_k] имеет вид диагональной матрицы с отличными от нуля только осевыми моментами инерции масс θ_ik (i = 1, 2, 3).

Величины масс m_k определяются в соответствии с требованиями учета нагрузок, изложенными в п. 2.1 СНиП II-7-81; моменты инерции масс θ_ik вычисляются по справочным данным.

2.4. Материал несущих конструкций сооружения (деформируемые связи) принимается упругим с линейной диаграммой деформирования. Модели этих конструкций (колонны, панели, диафрагмы и т.п.) определяются:

- видом представления масс (точка, диск, тело) в расчетной динамической, модели сооружения;

- моделями упругих связей в статической расчетной модели при формировании матриц жесткости (или податливости).

2.5. При геометрическом описании расчетных динамических моделей сооружения вводятся системы отчета и устанавливается связь между ними:

- вводится инерциальная система осей ОХ₁X₂Х₃ без привязки начала отсчета. В качестве этой системы осей может быть принята ориентация по странам света или какая-либо другая местная система трех ортогональных направлений;

- вводится система отсчета Ox₀₁x₀₂x₀₃, связанная с подвижным основанием (фундаментом), движение которого определяется моделью сейсмического воздействия; оси Ox₀₁x₀₂x₀₃ рекомендуется совмещать с основными осями здания;

- взаимосвязь между введенными системами осей определяется матрицей [φ₀], имеющей следующий вид:

,                              (1)

где φ⁰_ij - скалярные произведения единичных векторов j_x0i и j_xi (орт систем осей Ox₀₁x₀₂x₀₃ и ОХ₁Х₂Х₃); точкой отмечен знак скалярного произведения, поэтому φ⁰_ij являются тригонометрическими функциями углов между осями ОХ₁Х₂Х₃ и Ox₀₁x₀₂x₀₃; если системы отсчета ОХ₁Х₂Х₃ и Ox₀₁x₀₂x₀₃ одинаково ориентированы в пространстве (оси параллельны), то оператор (1) сводится к единичной матрице:

;

- однозначно и произвольно устанавливается нумерация инерционных элементов (материальных точек, дисков, тел) расчетной динамической модели сооружения I, II, III, ..., k, ..., n;

- с каждым из инерционных элементов неизменно связываются системы осей kx_k1x_k2x_k3. Для дисков и твердых тел эти оси должны быть главными и центральными осями инерции масс;

- определяется взаимосвязь между введенными системами осей kx_k1x_k2x_k3 и Ox₀₁x₀₂x₀₃. Эта взаимосвязь устанавливается с помощью операторов [φ_k], имеющих вид:

,                          (2)

где φ^k_ij - скалярные произведения орт j_x0i и j_xkj системы осей Ox₀₁x₀₂x₀₃ и kx_k1x_k2x_k3, поэтому φ^k_ij вычисляются аналогично (1); в частности, если оси kx_k1x_k2x_k3 и Ox₀₁x₀₂x₀₃ ориентированы в пространстве одинаково, то [φ_k] = Е - единичная матрица;

- векторами  определяются положения начала отсчета систем kx_k1x_k2x_k3 относительно осей Ox₀₁x₀₂x₀₃, где x⁰_ik - координаты центра масс тела в системе отсчета Ox₀₁x₀₂x₀₃ (i = 1, 2, 3).

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

2.6. Упругие свойства динамической модели сооружения определяются матрицей жесткости [r] или податливости [δ], которые рассчитываются по любым известным методам строительной механики.

2.7. Если расчет выполняется по методу сил (определение [δ]), то в центрах масс по направлению осей системы отсчета (или относительно их) последовательно прикладываются единичные силы Px_ik = 1 (или моменты Мx_ik = 1). От последовательного действия каждой единичной силы или момента выполняется статический расчет упругой системы. При этом для поперечных сечений упругих элементов определяются внутренние усилия, которые сводятся в табл. 1. Определяются также перемещения начал всех систем отсчета по направлению принятых осей и углы поворота данных систем отсчета относительно этих осей. В соответствии с принятой нумерацией масс из этих перемещений и углов поворота формируется матрица податливости системы [δ]. Результаты сводятся в табл. 2.

Таблица 1

Таблица 2

Таблица 3

Таблица 4

2.8. Если расчет упругой системы выполняется по методу перемещений, то на начало выбранных систем отсчета по направлению принятых осей (или относительно их) накладываются связи, которым последовательно сообщаются единичные перемещения Δx_ik = 1 (или углы поворота Δφ_ik = 1). Последовательно от каждого из этих перемещений или углов поворота выполняется статический расчет принятой упругой системы. При этом для поперечных сечений упругих элементов определяются внутренние усилия, которые сводятся в табл. 3. Определяются также реакции во всех наложенных связях. Из этих реакций, в соответствии с принятой нумерацией масс, формируется полная матрица жесткости системы [r]. Результаты формирования матрицы [r] сводятся в табл. 4.

Расчет режима свободных колебаний

2.9. Выполняется вспомогательный анализ полных матриц [r] [δ^-1]. Для выявления особенностей составляющих

 или ,                                              (3)

где

.                                                                  (4)

Таблица 5

Анализ (3) и (4) является необходимым, но недостаточным для определения частот и форм колебаний системы.

2.10. Выполняется анализ полных матриц [А] = [M^-1][r] или [B] = [А^-1] = [δ][M] для выделения обособленных составляющих

                                                  или

                                                   (5)

где                                                                                                                      (6)

Здесь [M_i] - полная матрица инерционных параметров системы, структура которой показана в табл. 5, а [M_i] - ее составляющие, соответствующие обособленным составляющим [r_i] или [δ_i].

2.11. Отдельно для каждой из обособленных составляющих [А_i] или [B_i] (если они имеются) решается задача на собственные значения (собственные числа и векторы) по частотным уравнениям.

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

                                                      (7)

                                                     (8)

где

 -

блочный вектор коэффициентов j-ой формы колебаний.

Если обособленные составляющие [A_i] или [B_i] отсутствуют, то все колебания в пространстве линейно взаимосвязаны и задача решается по полной матрице [А] или [B].

При решении прямой задачи (7) по [A] определяются собственные числа в виде Ω²_j (Ω_j - частота j-ой формы колебаний), а при решении обратной задачи (8) по [B] - собственные значения - 1/Ω²_j.

Каждая из обособленных составляющих [A_i] или [B_i] описывает свою независимую группу колебаний.

Результаты решения задачи по определению форм и частот сводятся в соответствующие таблицы (табл. 6). Такие таблицы составляются для каждой обособленной матрицы [А_i] или [B_i].

2.12. Коэффициенты форм колебаний системы нормируются по отношению к одному какому-либо коэффициенту, принимая его значение за единицу. В качестве такого нормирующего множителя удобно принимать значение коэффициента перемещения I массы в направлении 1-ой оси - Z^(j)_x11 = 1.

2.13. Критерием точности вычислений коэффициентов форм колебаний является удовлетворение их условиям ортогональности, имеющим вид

 (i ≠ j).                                             (9)

При проверке условий ортогональности строится матрица [С], компоненты которой C_ij определяются по (9) и при i = j, C_ii = M_i (М_i - приведенная к i-ой форме колебаний масса сооружений). Результаты проверки ортогональности форм колебаний сводятся в табл. 7.

Таблица 6

Таблица 7

Проверка условий ортогональности выполняется для каждой самостоятельной группы колебаний, определяемой соответствующими обособленными составляющими [A_i] или [В_i].

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

Рис. 1. Структура программы расчета режима собственных колебаний динамической системы

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

2.14. В соответствии с расчетной балльностью площадки строительства по СНиП II-7-81, определяется независимая от пространственной ориентации (инвариантная) расчетная интенсивность ускорений поступательного движения грунта в основании сооружения.

Расчетная инвариантная интенсивность ускорений поступательного движения грунта I является функцией от меньшего размера сооружения в плане В и вычисляется по формуле

I(В) = gАχ₁(В) (м/с²),                                              (10)

где g ≈ 10 м/с² - ускорение силы тяжести; А - коэффициент, принимаемый равным 0,1; 0,2; 0,4 для расчетной сейсмичности 7, 8 и 9 баллов по СНиП II-7-81.

Значение нормированной аппроксимирующей функции χ₁(В) определяется по ее графикам рис. 2 или вычисляется по формуле

χ₁(В) = e^{α(B - 25)},                                                         (11)

где В(м) - меньший размер сооружения в плане; α - аппроксимирующий коэффициент, значение которого равно -8 · 10^-4; -4,8 · 10^-3 и -1,2 · 10^-2 (1/м) для I, II и III категории грунтов по СНиП II-7-81.

Рис. 2. Графики нормированных функций χ₁(В) = e^{α(B - 25)} при α = -8 · 10^-4; -4,8 · 10^-3 и -1,2 · 10^-2 (1/м) для I, II и III категорий грунтов, соответственно

2.15. В зависимости от категории грунтов по СНиП II-7-81 определяется независимая от пространственной ориентации (инвариантная) интенсивность угловых ускорений вращательных движений грунта в основании сооружений (ускорения ротации поля сейсмических движений). Этот инвариант интенсивности определяется относительной (по отношению к инварианту ускорения поступательного движения) величиной W.

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

 (м^-1),                                                  (12)

где  = 2 · 10^-2; 6 · 10^-2 и 9 · 10^-2 (м^-1) для грунтов I, II и III категорий по СНиП II-7-81.

Значение нормированной аппроксимирующей функции χ₂(В) определяется по графикам рис. 3 или вычисляется по формуле

χ₂(В) = e^{b(B - 25)},                                                         (13)

где В(м) - меньший размер сооружения в плане; b - аппроксимирующий коэффициент, значение которого равно -7,2 · 10^-3; -1 · 10^-2 и -1,6 · 10^-2 (1/м) для I, II и III категорий грунтов по СНиП II-7-81.

Рис. 3. Графики нормированных функций χ₂(В) = e^{b(B - 25)} при b = -7,2 · 10^-3; -1 · 10^-2 и -1,6 · 10^-2 (1/м) для I, II и III категорий грунтов, соответственно

2.16. Спектральный состав компонент векторов сейсмического воздействия (ускорения поступательного движения  и углового ускорения ) определяется соответствующими коэффициентами динамичности, в зависимости от категории грунтов по СНиП II-7-81 - k_гр, коэффициентов потерь энергии при колебаниях рассматриваемой системы γ = δ/π (δ - логарифмический декремент колебаний), периодов собственных колебаний - Т = 2π/Ω (Ω - частота собственных колебаний, определяемая в соответствии с п.п. 8-12 Приложения 1).

Для вектора ускорения поступательного движения грунта значения коэффициента динамичности определяются по СНиП II-7-81 и вычисляются по формуле

,                                                 (14)

где k_гр - коэффициент, учитывающий категорию грунтов по СНиП II-7-81, значения которого принимаются по данным табл. 8.

Таблица 8

Таблица 9

Рис. 4. Графики нормированных коэффициентов динамичности вектора ускорения поступательного движения сейсмического воздействия  для грунтов I, II и III категорий

Потери энергии при колебаниях системы учитываются коэффициентом α_γ, значения которого определены для (14) в соответствии со СНиП II-7-81 и приведены в табл. 9.

Значения нормированных по соответствующим максимумам коэффициентам динамичности  определяются в зависимости от категории грунтов по СНиП II-7-81 по графикам рис. 4 или вычисляются по формулам:

для I категории грунтов

                                         (15)

для II категории грунтов

                                    (16)

для III категории грунтов

                                    (17)

Для вектора углового ускорения вращения сейсмического воздействия значения коэффициента динамичности вычисляются по формуле

,                                                 (18)

где k_гр - коэффициент, учитывающий категорию грунтов по СНиП II-7-81, значения которого определяются по табл. 8.

Таблица 10

Примечание. Для других значений γ допускается коэффициент b_γ определять по интерполяции.

Значения коэффициента b_γ, учитывающего потери энергии при колебаниях системы, определяются по табл. 10.

Нормированные значения коэффициентов динамичности  определяются по графику рис. 5 или вычисляются по формуле

                                           (19)

Рис. 5. График нормированного коэффициента динамичности вектора углового ускорения вращения сейсмического воздействия

2.17. Сейсмическое воздействие имеет произвольное направление в пространстве (см. п. 2.3 СНиП II-7-81). Ориентация в пространстве векторов ускорения поступательного движения и углового ускорения вращения грунтового основания сооружения определяется направляющими косинусами.

Вектор ускорения поступательного движения грунта в основании сооружения задается в инерциальных осях ОХ₁Х₂Х₃, в качестве которых, в соответствии с п. 2.5, может быть принята ориентация по сторонам света или какая-либо другая местная система трех ортогональных направлений. Ориентация в пространстве вектора  определяется направляющими косинусами  в этих инерциальных осях ОХ₁X₂Х₃. Значения направляющих косинусов  должны удовлетворять условиям нормировки, имеющим вид

,                                                        (20)

Вектор углового ускорения вращения грунта задается в осях Оx₀₁,x₀₂x₀₃, неизменно связанных с основанием сооружения. В качестве осей Оx₀₁x₀₂x₀₃, согласно п. 2.5, следует принимать систему главных осей сооружения. Ориентация в пространстве вектора  определяется направляющими косинусами , которые отсчитываются относительно этих же осей Оx₀₁x₀₂x₀₃. Значения направляющих косинусов  должны удовлетворять условиям нормировки в виде

.                                                       (21)

Расчетные параметры ориентации векторов сейсмического воздействия  и  определяются значениями направляющих косинусов  и  (i = 1, 2, 3) в указанных осях. В расчетах следует исходить из двух групп значений направляющих косинусов:

проверочные (наиболее вероятные для заданного региона строительства) и основные расчетные.

Расчет на наиболее вероятные значения направляющих косинусов является проверочным и выполняется в тех случаях, когда они известны для данного региона строительства. Эти значения  и  определяются по данным конкретной сейсмологической обстановки: геология региона местности; возможные очаги землетрясений и направления на них; опыт повреждений и разрушений во время прошедших землетрясений; данные обработки зарегистрированных процессов движения грунта и т.д. При этом принимаемые проверочные значения направляющих косинусов должны удовлетворять условиям нормировки (20) и (21).

Расчет на основные значения направляющих косинусов  и  является обязательным во всех случаях.

Основные расчетные значения направляющих косинусов вычисляются отдельно для каждой из учитываемых форм колебаний, исходя из условий максимума динамической реакции по рассматриваемой форме колебаний. Эти значения направляющих косинусов для j-ой формы колебаний определяются по следующим формулам:

 и ,                                      (22)

где

                                      (23)

Результаты вычислений расчетных значений направляющих косинусов сводятся в табл. 11.

Таблица 11

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

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

,                                                                                                                                       (24)

где

Таблица 12

Максимальное значение коэффициента динамичности (24) имеет место при i = j, т.е. β⁽ⁱ⁾_j < β^(j)_j при i ≠ j.

Результаты вычислений коэффициентов динамичности сводятся в табл. 12.

2.19. Если расчет упругой системы выполняется по методу сил и в соответствии с требованиями п. 2.7 выполнялось построение матрицы податливости [δ] (см. табл. 2), а также вычислялись внутренние усилия в поперечных сечениях элементов от единичной нагрузки (см. табл. 1), то основными расчетными параметрами являются векторы сейсмических сил и моментов, которые определяются по следующим формулам:

                                    (25)

Коэффициент k в (25) соответствует k₁ СНиП II-7-81, значения которого согласно табл. 3 СНиП II-7-81 предназначены для сооружений симметричных в плане, а так же для зданий регулярной конструктивной схемы, когда эффектами пространственной работы конструкций можно пренебречь.

Для сооружений промышленного или гражданского назначения, имеющих сложную конструктивную схему, необходим учет пространственной работы конструкций. В этих случаях значения коэффициента k отличаются от соответствующего значения k₁ по СНиП II-7-81. Поэтому значение коэффициента k следует принимать по специальному обоснованию в соответствии с принятой конструктивной схемой сооружения и ее возможностями к развитию пластических деформаций и местных повреждений, не приводящих к выходу из строя сооружения в целом. Предварительное значение коэффициента k должно задаваться в задании на проектирование и затем корректироваться в ходе выполнения проектных работ и уточняться в соответствии с выбранной конструктивной схемой с последующим окончательным утверждением заказчиком.

Значения векторов сейсмических сил и моментов вычисляются для каждой из рассчитываемых форм колебаний и по каждому варианту ориентации векторов воздействия. Для каждой формы колебаний максимальные значения векторов сил и моментов определяются при значениях направляющих косинусов по (22). Результаты расчета сводятся в табл. 13.

Таблица 13

Таблица 14

2.20. Если расчет упругой системы выполняется по методу деформаций и в соответствии с требованиями п. 2.8 выполнялось построение матрицы жесткости [r] (см. табл. 4), а также вычислялись внутренние усилия в поперечных сечениях упругих элементов от единичных перемещений, то основными расчетными параметрами являются векторы перемещений центров масс и углов вращения масс, которые определяются по следующим формулам:

         (26)

Значения векторов перемещений центров масс и углов поворота масс относительно их центра вычисляются для каждой из рассчитываемых форм колебаний и по каждому варианту ориентации векторов воздействия. Для каждой формы колебаний значения векторов (26) имеют максимальные величины при направляющих косинусах по (22). Результаты этого расчета сводятся в табл. 14.

Определение расчетных значений внутренних усилий в поперечных сечениях несущих конструкций сооружения

2.21. Если расчет упругой системы выполнялся по методу сил, то для каждой отдельной формы колебаний, по данным п. 2.19, по величинам фактических сейсмических сил и моментов, корректируются внутренние усилия от единичной нагрузки по данным п. 2.7:

                        (27)

где N^(j,i)₁, N^(j,i)₂, N^(j,i)₃ - поперечные и продольные силы по j-ой форме при i-ой ориентации векторов сейсмического воздействия; M^(j,i)₁, M^(j,i)₂, M^(j,i)₃ - изгибающие и крутящие моменты по j-ой форме при i-ой ориентации векторов сейсмического воздействия для рассматриваемого сечения  и  - единичные усилия; S^(j,i)_lk, M^(j,i)_lk - фактические сейсмические силы и моменты по j-ой форме колебаний для i-ой ориентации векторов сейсмического воздействия.

Внутренние усилия (27) вычисляются по каждой форме колебаний и для каждого из вариантов ориентации векторов сейсмического воздействия. Для каждой формы колебаний максимум усилий (27) имеет место при значениях направляющих косинусов по (22). Результаты этого расчета сводятся в табл. 15.

2.22. Если расчет упругой системы выполняется по методу деформаций, то для каждой формы колебаний внутренние усилия от единичных перемещений и углов вращения по данным п. 2.8 корректируются на величины фактических перемещений и углов вращения масс по данным п. 2.20:

                         (28)

где  и  - единичные усилия; x^(j,i)_lk и α^(j,i)_lk - фактические перемещения и углы вращения масс по j-ой форме колебаний для i-ой ориентации векторов воздействия.

Таблица 15

Таблица 16

Таблица 17

По (28) внутренние усилия вычисляются для каждой формы колебаний и для каждого варианта ориентации векторов сейсмического воздействия. Максимум усилий для каждой формы колебаний имеет место при значениях направляющих косинусов по (22). Результаты этого расчета сводятся в табл. 16.

2.23. Для каждого i-го варианта ориентации векторов сейсмического воздействия по внутренним усилиям, соответствующим отдельным формам колебаний, N^(j,i)_l и M^(j,i)_l определяются расчетные значения внутренних усилий  и

 и ,                                      (29)

где f - число учитываемых форм колебаний при определении расчетных значений внутренних усилий в l-ом направлении (l = 1, 2, 3) поперечного сечения для рассматриваемого i-го варианта ориентации векторов сейсмического воздействия.

Результаты расчета сводятся в табл. 17.

2.24. Из расчетных усилий (29), соответствующих отдельным i-ым вариантам ориентации векторов сейсмического воздействия, по данным табл. 17 принимаются максимальные значения усилий для проверки несущей способности сечений в сочетании с усилиями от других видов нагрузок. Принятый i-ый вариант значений  и  определяет наиболее опасную ориентацию векторов сейсмического воздействия для рассчитываемого сооружения, при которой в рассматриваемом сечении данной конструкции возникают максимальные значения внутренних усилий.

Приложение 1

ПРИМЕР РАСЧЕТА ПО ОПРЕДЕЛЕНИЮ СЕЙСМИЧЕСКОЙ НАГРУЗКИ НА КОЛОННЫ СИЛОСНОГО КОРПУСА

Приведенный в разделе 2 алгоритм расчета проиллюстрирован на простейшем примере определения нагрузки на колонны силосного корпуса (рис. 6). Параметры сооружения взяты из проекта железобетонного силосного корпуса, разработанного институтом Сибгипробиосинтез (г. Красноярск) для Андижанского гидролизного завода. В соответствии с проектом сечение колонн 0,6×0,6 м.

Описание расчетной динамической модели сооружения

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

2. Согласно п. 2.3 инерционные параметры при полностью заполненном резервуаре равны: m₁ = 636 · 10³ кг; θ_1I = θ_2I = 1324 · 10⁴ кг/м²; θ_3I = 286 · 10⁴ кг/м².

3. По п. 2.4 железобетонные колонны моделируем упругими стержневыми конечными элементами с геометрическими характеристиками в соответствии с рис. 6.

Рис. 6. Конструктивная схема (а) и расчетная динамическая модель (б) силосного корпуса

4. Согласно п. 2.5 определяются геометрические параметры динамической модели:

- оси ОХ₁Х₂Х₃ принимаем параллельными осям Ox₀₁x₀₂x₀₃;

- начало осей Ox₀₁x₀₂x₀₃ связываем с геометрическим центром фундамента (см. рис. 6);

- [f₀] = Е - единичная матрица;

- твердому цилиндру присваиваем номер 1;

- вводим главные центральные оси инерции цилиндра Ix₁₁x₁₂x₁₃;

- [φ_I] = Е;

-  (м).

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

5. По п. 2.6. для статического расчета рассматриваемой упругой системы воспользуемся программами, основанными на конечном элементе.

6. В соответствии с п. 2.7 от последовательного действия единичных сил и моментов выполняем статический расчет и определяем внутренние усилия в рассматриваемых сечениях, а также перемещения центра массы (табл. 1) и углы поворотов цилиндра. В качестве примера в табл. 18 сведены внутренние усилия для нижнего сечения колонны. В табл. 19 сведены единичные перемещения и углы поворотов. Табл. 19 определяет матрицу [δ].

Расчет режима свободных колебаний

7. Согласно п. 2.9 анализ матрицы [δ] (табл. 19) показывает, что она не содержит обособленных составляющих.

Таблица 18

Единичные внутренние усилия в нижнем сечении колонны