УДК 624.044.3:699.841
DOI: 10.14529/build160102
ОБ ОСОБЕННОСТЯХ МЕТОДОЛОГИИ НЕЛИНЕЙНОГО СТАТИЧЕСКОГО АНАЛИЗА И ЕГО СОГЛАСОВАННОСТИ С БАЗОВОЙ НОРМАТИВНОЙ МЕТОДИКОЙ РАСЧЁТА ЗДАНИЙ И СООРУЖЕНИЙ НА ДЕЙСТВИЕ СЕЙСМИЧЕСКИХ СИЛ
А.В. Соснин
Московский государственный университет путей сообщения (МГУПС (МИИТ)), Смоленский филиал, г. Смоленск
В случае интеграции метода нелинейного статического анализа (далее - НСМ) в практику проектирования зданий и сооружений в сейсмических районах России, процедура его применения должна быть понятна инженерам. Рассматриваются основные допущения и преимущества методологии НСМ, а также его согласованность с линейно-спектральным методом расчёта (далее - ЛСМ). Приводятся примеры сооружений, расчёты которых выполнены по НСМ в том числе в комплексе SAP2000.
Главная задача проектирования зданий и сооружений в сейсмических районах - предотвращение обрушения системы в целом при редком (сильном) землетрясении. Линейно-спектральный метод расчёта (далее - ЛСМ) не всегда гарантирует требуемый уровень сейсмостойкости зданий, иначе их обрушения не происходили бы. При решении проблемы с помощью положений, удовлетворяющих требованиям СП 14.13330.2014 [1], возникает ряд трудностей [2-4], к основным из которых относятся: чувствительность к алгоритмам конечно-элементного моделирования, резкая избирательность акселерограмм к отдельным собственным частотам объекта и отсутствие в нормах «прозрачной» методики.
Для повышения качества расчётной оценки сейсмостойкости многоэтажных зданий массовой застройки в качестве основного инструментария автор предлагает использовать методологию нелинейного статического расчёта [5, 6]. Полное наименование метода, указываемое в зарубежной нормативной документации и научных публикациях - Nonlinear Static Analysis Procedure (Pushover Analysis). Необходимость обсуждения возможности практического использования НСМ спровоцирована своеобразным скептическим отношением некоторых отечественных практиков и учёных к процедурам в нём реализуемым, которые, кажутся им излишне упрощёнными и дискуссионными. При этом основной акцент делается на отсутствие ясного алгоритма расчёта. В статье выполнен обзор основных допущений предлагаемой методологии, которая, действительно, имеет ряд нюансов [7, 8].
В настоящее время НСМ применяется при проектировании сейсмостойких конструкций, зданий и сооружений в различных странах мира, в том числе Китае, Индии, США, Иране, Италии, Греции, Индонезии и др. Входящие в него расчётные процедуры интегрированы в такие комплексы, как, например, SAP2000, ETABS, Midas Gen, FESPA, SOFISTIK, SCADA PRO, INDYAS, DRAIN-2DX, SeismoStuct, ADAPTIC, RUAUMOKO 3D. В англоязычной терминологии полное наименование метода в некоторых его модификациях (например, [9]) часто сопровождается словом-префиксом "simplified". Дословный перевод этого слова имеет первое общепринятое значение «упрощенный». Учитывая выполнение оценки сейсмостойкости сооружения с использованием результатов в формате ADRS [5, 6], автор работы считает, что перевод слова "simplified" необходимо рассматривать в семантической интерпретации «адаптированный» или «интуитивно понятный». А в контексте сравнения НСМ с процедурой расчёта в нелинейной динамической постановке, технически корректнее применять слово «simplified» в смысловом переводе «менее сложный». Такой транслейт-оттенок хорошо согласуется с трактовкой, приведённой в Большом толковом словаре ([10, с. 1394 ]). Преимущества НСМ и большая достоверность результатов его применения в оценке сейсмостойкости (по сравнению с ЛСМ), - очевидны, и подтверждается резолютивной частью многочисленных зарубежных исследований. В случае внедрении НСМ в отечественную практику проектирования, расчётные положения метода должны быть понятны рядовым инженерам. В статье речь пойдёт о подоснове Pushover методологии - процедуре &А& метода спектра несущей способности (англ. Capacity Spectrum Method) [11]. С базовыми терминами и определениями НСМ можно ознакомиться в работах [2-6].
Что же концептуально представляет собой рассматриваемый метод? НСМ - это диаграммный метод расчёта, применяя который инженер как бы сознательно допускает, что передача сейсмических
сил сооружению произойдет настолько медленно, что за учитываемое время между внутренними напряжениями в несущих конструкциях и действующими на систему квазистатическими силами обеспечится равновесие. Принимается, что вся приходящаяся на сооружение энергия землетрясения передается системе при одностороннем параметрически изменяющемся загружении. Уравнение движения многомассовой системы в постановке рассматриваемого метода записывается в виде:
М -и + R = -М - I - а
где М - матрица масс; и - вектор перемещений; Р - вектор внутренних сил; I - единичный вектор; - ускорение грунта основания как функция времени.
Параметры диссипации в уравнении (1) отсутствуют, поскольку будут учитываться при построении расчетного спектра реакции Sa(T). Вектор внутренних сил К эквивалентен вектору внешних сил Р (вектору сдвигающей силы Vgh.fr), который определяется из выражения
Р = р - М-Ф, (2)
здесь р - параметр, контролирующий величину
горизонтальной реакции; Ф - вектор форм, не меняющийся в процессе колебаний системы.
Вектор перемещений определяется из выражения
и = Ф-щ , (3)
здесь щ - перемещение верха системы.
При этом матрица-вектор форм Ф для упрощения расчетов нормализуется таким образом, чтобы компоненты в верхней части системы равнялись единице. Уравнение (1) тогда выглядит так:
т т т
Ф -М -Ф-щ + р-Ф -М-Ф = -Ф -М-1-аё. (4)
Одним из результатов выполнения расчета по НСМ является кривая несущей способности (рис. 1).
Каждому конкретному особому сочетанию нагрузок будет соответствовать своя кривая равновесных состояний системы. Важнейшим нюансом метода является алгоритм конвертации этой кривой из формата «горизонтальная сейсмическая реакция (сдвигающая сила) в основании - перемещение верхней точки системы» в спектр несущей способности формата ADRS [2-6]. Вычисляются коэффициент преобразования (трансформации) PF1 и эффективная модальная масса ^ (в доРис. 1. Графическая визуализация зависимости «горизонтальная сейсмическая реакция - перемещение верхней точки системы» с участками, характеризующими степень повреждения основных несущих конструкций
Примечания:
лях единицы), соответствующая основной форме колебаний (коэффициент модальной массы а1). Коэффициент а1 определяется по формуле:
" N 12
I (^-Фи) _«=1 g _
I (^ + Фи) _ 1=1 g _
Коэффициент трансформации может принимать значение от 1,0 до 1,6; определяется он из следующего выражения:
V (^) 1=1 g
V (^-Ф12„) 1=1 g
здесь N - уровень верха (основного покрытия)
М>системы (здания); — - масса 1-й точки (этажа;
яруса) системы (здания); ф1г - амплитуда колебаний г-й точки (этажа; яруса) системы (здания) по первой основной форме колебаний.
Оба указанных коэффициента зависят от формы приложения сейсмических сил по высоте конструктивной системы (см. таблицу).
Принимается, что схема приложения сейсмических сил остаётся неизменной до окончания итерационного расчёта, соответствующего рассматриваемому направлению сейсмического воздействия.
Одним из основополагающих допущений №Р-процедур является преобразование системы со многими степенями свободы (MDOF), соответствующей параметрам рассматриваемой пространственной системы, в эквивалентный ей осциллятор
Значения коэффициента модальных масс и коэффициента трансформации в зависимости от принятой схемы распределения сейсмических сил
Тип схемы1 Условное наименование схемы Выражение для определения сейсмических сил2 Эффективная 3 модальная масса (a), в д.ед. Коэффициент преобразования (PF)
V w, - h 1 1
V w, - h
Примечания:
тип 1 - низшей формой собственных колебаний простейшей консольной модели;
тип 2 - эквивалентным прямолинейным очертанием низшей формы колебаний двумерной модели простой рамной системы;
тип 3 - расчётным откликом системы, определённым суммированием вкладов всех учитываемых форм по любому из правил суммирования (SRSS, CQC, метод 10 %, DSC метод Розенблюэта, метод двойных сумм NRC, метод Лендли-Яу, метод Гупты, GMC);
тип 4 - первой формой собственных колебаний системы с первым «гибким» этажом.
wx, wi - вес яруса системы (индексом х обозначен ярус, на который действует вычисляемая сейсмическая сила);
Ax - горизонтальное перемещение системы в пределах яруса, на который действует вычисляемая сила Рх;
hi - высота яруса системы;
Vsh.b - горизонтальная сейсмическая реакция системы в уровне основания.
с ^=5 % от критического (SDOF) [3]. Принимается, что при внешнем воздействии масса осциллятора Мщ отклоняется на расстояние ¿<цеЩ, что соответствует перемещению покрытия сооружения на величину А (см. рис. 1). Применяя полученный в выражении (6) коэффициент трансформации, получаем
Д = Р^ - Б^^ , (7)
отсюда
Sd,eq (8)
PF1 -Ф1,roof
Тогда уравнение движения (1) для эквивалентного одномассового осциллятора можно записать в виде:
Meq * Sd,y,eq + Feq - "Meq & ag • (9)
Преобразование MDOF в SDOF обусловлено необходимостью обеспечения корректного наложения спектра несущей способности на спектр реакции (который, как известно, строится для SDOF).
Требуемый расчётный спектр реакции определяется только с учётом положения траектории поиска точки упруго-пластических свойств системы (англ. performance (target) point), независимо от кривой несущей способности. Преобразование кривой несущей способности в спектр выполняется в соответствие со следующими выражениями:
PF1 -Ф1, roof
здесь Sa - спектральное ускорение; Sd - спектральное перемещение; Vshb - горизонтальная сейсмическая реакция в уровне низа (средней планировочной поверхности земли) системы (здания); ^ - расчетный сейсмический вес системы (здания); Агоо! - перемещение верха (основного покрытия)
системы (здания).
Построение спектра реакции ¿а(Т) в формате ADRS производится так:
Sd - Sa (T) •
В качестве оценочного критерия сейсмостойкости, среди прочих, рассматривается положение точки свойств [2, 5, 6]. Определив величину оценочного перемещения можно проанализировать наступление предельного состояния по критерию необрушения. В графическом виде способ получения значения этого перемещения представлен на рис. 2.
Итерационный расчет проводится до выполнения условия сходимости:
If < 0,05.
После определения конечного значения соответствующего расположению точки свойств, важно не забывать преобразовать его в формат MDOF (рис. 2, этап 5 ):
t arg. disp
- PFV Фи • Sd .
После прохождения всех вышеуказанных этапов расчета можно получить информацию о коэффициентах податливости системы и о коэффициентах редукции К к упругому спектру реакции (как в [5]).
Рис. 2. Графическая визуализация пошагового алгоритма применяемой расчётной процедуры метода нелинейного статического анализа
Характерной чертой NSP-процедур является возможность учёта понятия энергоёмкости. В методе CSM потенциальная энергия деформирования (фактически равная площади фигуры, расположенной под спектром несущей способности) определяется упрощённо: значение Es0 находится из площади прямоугольного треугольника (см. рис. 2 [6]):
Eso = 0,5• (api • dpl). (15)
В той же работе [6] приведена обобщённая графическая интерпретация алгоритма, применяемого для оценки энергий в указанной расчётной процедуре, и положения по определению поправочного коэффициента к [4, с. 82] для затухания, характеризующего работу системы в упруго-пластической стадии.
Автор согласен, что изложенная методология не свободна от недостатков и допущений. При том, что НСМ присущи почти все допущения теории предельного равновесия А.А. Гвоздева [12], а качество оценки сейсмостойкости объекта напрямую зависит от грамотности и инженерной интуиции расчётчика, недостатки метода перекрываются
его главным преимуществом - возможностью анализа реакции системы в момент приближения её к предельному состоянию при действии сейсмических сил [2-6]. Расчётные процедуры НСМ позволяют упразднить сразу несколько частных коэффициентов, характерных ЛСМ в редакции свода правил СП 14.13330.2014 [1]: коэффициент K учитывается непосредственной упруго-пластической работой конструкций; а значения А, ßi и Kv рассматриваются при построении спектра реакции Sa(T).
Метод нелинейного статического анализа применяется для оценки сейсмостойкости многоэтажных зданий, стальных пространственных конструкций, стадионов, мостов, и даже таких сооружений, как подземные станции метро и даже дамбы. В качестве частного примера можно привести известный национальный стадион Beijing National Stadium (рис. 3), построенный в Пекине для проведения летних Олимпийских игр 2008 года, расчёты которого также выполнялись и в комплексе SAP2000.
Сейчас, в период разработки национальных приложений к Еврокодам, важно накопить отечественный практический опыт использования НСМ
Рис. 3. Национальный Олимпийский стадион в Пекине (Beijing National Stadium, также известный как «Птичье гнездо» / Bird Nest) [13]: а - общий вид сооружения; б - расчётная модель стадиона; в - спектры несущей способности системы при направлении сейсмических нагрузок вдоль оси Х; г - спектры несущей способности системы при направлении сейсмических сил вдоль оси Y
Примечания:
в расчётах зданий и сооружений и провести широкое обсуждение проблем, связанных с его применением на фоне нормативных и обозримых альтернативных методик расчёта. Эффективным способом для реализации дискуссий могут стать доклады на интеллектуальных площадках крупных строительных вузов.
В заключение (в качестве экспертного мнения) автор предлагает рассмотреть пояснения профессора Эдварда Вильсона (E.L. Wilson; Калифорнийский университет, Беркли, США), данные им в эссе [14]. Профессор Э. Вильсон отмечает, что его практический опыт в оценках сейсмостойкости сооружений прямым динамическим методом и методом Pushover показал различные результаты, но вероятность ошибки в оценке отклика той же системы спектральным методом будет значительно больше.
Выводы
Благодарность автора
Автор выражает искреннюю благодарность доктору технических наук, профессору кафедры «Сопротивление материалов» Московского государственного строительного университета (МГСУ) Джинчвелашвили Гураму Автандиловичу за дружескую поддержку и наставления.
for the Seismic Analysis of Structures / E.L. Wilson, 2014. - 2 p. - http://www.edwilson.org/History/ Pushover.pdf
Соснин Алексей Викторович, старший преподаватель кафедры «Здания и сооружения на транспорте», Московский государственный университет путей сообщения (МГУПС (МИИТ)), Смоленский филиал, syabryauskas@mail.ru.
Поступила в редакцию 22 декабря 2015 г.
DOI: 10.14529/build160102
ON PECULIARITIES OF THE NONLINEAR STATIC ANALYSIS AND ITS COORDINATION WITH THE STANDARD CALCULATION PROCEDURE (CSM) OF BUILDINGS AND STRUCTURES UNDER SEISMIC LOADS
A.V. Sosnin, syabryauskas@mail.ru
Moscow State University of Railway Engineering, Smolensk Branch, Smolensk, Russian Federation
The nonlinear static analysis should be clear for engineers in case of its integration into building design in seismic regions of Russia. Basic assumptions and advantages of the analysis, as well as its coordiantion with the linear spectral calculation method are considered by the author. The examples of structures calculated according to NSA, in particular in SAP2000 software are given.
References
Received 22 December 2015
ОБРАЗЕЦ ЦИТИРОВАНИЯ
Соснин, А.В. Об особенностях методологии нелинейного статического анализа и его согласованности с базовой нормативной методикой расчета зданий и сооружений на действие сейсмических сил / А.В. Соснин // Вестник ЮУрГУ. Серия «Строительство и архитектура». - 2016. - Т. 16, № 1. - С. 12-19. Ш1: 10.14529/ЬшИ160102
FOR CITATION
Sosnin A.V. On Peculiarities of the Nonlinear Static Analysis and its Coordination with the Standard Calculation Procedure (CSM) of Buildings and Structures under Seismic Loads. Bulletin of the South Ural State University. Ser. Construction Engineering and Architecture. 2016, vol. 16, no. 1, pp. 12-19. (in Russ.). DOI: 10.14529/build160102