автореферат диссертации по строительству, 05.23.17, диссертация на тему: Динамический расчет многоэтажных зданий при землетрясении
Автореферат диссертации по теме "Динамический расчет многоэтажных зданий при землетрясении"
На правах рукописи
Нгуен Куок Донг
ДИНАМИЧЕСКИЙ РАСЧЕТ МНОГОЭТАЖНЫХ ЗДАНИЙ ПРИ ЗЕМЛЕТРЯСЕНИИ
05.23.17 - Строительная механика
АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук
Работа выполнена на кафедре строительной механики в ГОУ ВПО «Санкт-Петербургский государственный архитектурно-строительный университет».
Научный руководитель: заслуженный деятель науки и техники РФ,
доктор технических наук, профессор, почетный член РААСН Масленников Александр Матвеевич (СПБГАСУ)
Официальные оппоненты: доктор технических наук, профессор
Рутман Юрий Лазаревич (ОАО КБ специального машиностроения, г. Санкт Петербург);
кандидат технических наук, ст. препод. Дмитровская Любовь Николаевна (ПГУПС)
Ведущая организация: ГОУ ВПО «Санкт-Петербургский
государственный политехнический университет»
Защита диссертации состоится 24 декабря 2010 г. в 14 ч.ЗО мин. на заседании диссертационного совета Д 212.223.03 при ГОУ ВПО «Санкт-Петербургский государственный архитеюурно-строительный университета по адресу: 190005, Санкт-Петербург, 2-я Красноармейская ул., д. 4, зал заседаний (219). Факс: (812) 316-58-72
С диссертацией можно ознакомиться в фундаментальной библиотеке ГОУ ВПО «Санкт-Петербургский государственный архитектурно-строительный университет».
Автореферат разослан ноября 2010 г.
Ученый секретарь диссертационного
совета д.т.н., проф. ^ Л.Н.Кондратьева
ОБЩАЯ ХАРАКТЕРИСТКА РАБОТЫ Актуальность темы диссертации. Ежегодно на земном шаре, включая Вьетнам, происходит свыше 300 тыс. землетрясений, большинство из которых имеет небольшую силу или проявляется в ненаселенных районах. Однако некоторые очаги сильных землетрясений располагаются близко к населенным районам. По переписи на территории Вьетнама было более 13 землетрясений свыше 7 баллов. Такие землетрясения сопровождаются разрушениями несей-смостоиких зданий и сооружений, гибелью людей и уничтожением материальных и культурных ценностей, накопленных трудом многих поколений.
В настоящее время чрезвычайно актуальными являются исследование конструкций в процессе разрушений, учет влияния податливости основания на колебания конструкций и защита гражданских высотных зданий, находящихся в сейсмически активных районах. При сильных землетрясениях в процессе колебаний происходит частичное разрушение конструкций (отслоение бетона, обрыв связей, появление пластических деформации и т. п.). В таких случаях сооружение становится конструктивно нелинейным. Возможные места разрушения устанавливаются на основании статического расчета. В большинстве случаев расположение мест образования пластических шарниров в конструкциях при действии статической и динамической нагрузок совпадают. Поэтому исследование в процессе разрушении конструкции является одной из важнейших проблем для повышения сейсмостойкости конструкции.
При расчете зданий на сейсмическое воздействие по СНиП учет основания производится через коэффициент динамичности, который зависит от категории грунта по сейсмическим свойствам и от периодов собственных форм колебаний. При расчете высоких зданий на воздействие фрагментов реальных акселерограмм непосредственное применение коэффициента динамичности не оговаривается. В диссертации приводится исследование влияния податливости основания на колебания зданий путем непосредственного расчета по акселерограммам.
Повышение сейсмостойкости зданий достигается различными способами. Один их них применение систем специальной сейсмозащиты, обеспечивающих проскальзывание зданий относительно их оснований. В настоящее время благодаря созданию новых конструкционных материалов, прочных и долговечных пластмасс, идея сейсмоизоляции зданий с помощью элементов скольжения получила логическое развитие и конструктивное оформление. В диссертации предлагается расчет систем с сейсмоизолирующим скользящим поясом на воздействие фрагментов реальных акселерограмм. Изложенное выше свидетельствует об актуальности предлагаемой темы исследования. Цели и задачи диссертационной работы
Диссертационная работа направлена на исследование поведения статически неопределимых стержневых систем при динамическом воздействии в виде колебания основания, что имеет место при землетрясении. Для достижения указанной цели предлагается осуществить решение следующих задач: - создание специального численного метода решения поставленных задач;
- исследование совместимости результатов расчета по спектральному методу с расчетом на основе акселерограмм;
- исследование влияния коэффициента поглощения энергии на колебание зданий при землетрясении;
- исследование влияния податливости основания на колебание зданий при землетрясении;
- исследование эффективности систем с сейсмоизолирующим скользящим поясом с нелинейным демпфированием;
- разработка алгоритма расчета и составление программы учета дефектов, возникающих в здании при землетрясении.
Научная новизна диссертации:
- составлена программа решения системы дифференциальных уравнений движения сосредоточенных масс;
- на основании расчетов по реальным акселерограммам установлена область достоверности расчетов по спектральной теории;
- созданы алгоритм и программа расчета влияния податливости основания на колебание зданий при землетрясении;
- созданы алгоритм и программа учета дефектов в здании при землетрясении;
- создана математическая модель, алгоритм и программа расчета систем с сейсмоизолирующим скользящим поясом с нелинейным демпфированием.
Практическая значимость. Разработаны алгоритмы и на их основе составлены программы, которые позволяют еще на стадии проектирования конструкций исследовать их работоспособность при динамических воздействиях различного характера и дать рекомендации, направленные на улучшение их работы. Приведен ряд практических рекомендаций, полученных в процессе расчетов.
Достоверность разработанных теоретических положений, математических моделей и выводов обеспечивается корректностью постановки задачи, использованием общепринятых в строительной механике допущений, тестированием программного комплекса, а также сопоставлением с работами других авторов.
На защиту выносятся:
- программы сейсмического расчета численным методом многоэтажных зданий;
- результаты оценки достоверности спектрального метода на основе расчета по акселерограммам;
- результаты влияния коэффициента поглощения на колебания при землетрясении;
- результаты расчета и выводы по влиянию податливости основания на колебание железобетонных и металлических конструкций при сейсмическом воздействии;
- математическая модель сейсмоизолирующего скользящего пояса с нелинейным демпфированием;
- результаты влияния дефектов высоких зданий, возникающих при различных землетрясениях;
- выявление динамического эффекта, возникающего при внезапном разрушении в процессе землетрясения.
Апробация работы и публикации. Основные результаты диссертационной работы доложены и обсуждены на 61-й и 62-й международных научно-технических конференциях молодых ученых, Санкт-Петербург, СПбГАСУ, 2007 г. и 2008 г.; на 66-й научной конференции профессоров, преподавателей, научных работников, инженеров и аспирантов университета, Санкт-Петербург, СПбГАСУ, 2009 г.; на ХХШ Международной конференции «Математическое моделирование в механике деформируемых тел и конструкций. Методы граничных и конечных элементов». СПб., 2009. Основное содержание диссертации опубликовано в 6 печатных работах, в том числе две статьи в изданиях из перечня ВАК РФ.
Структура и объем диссертации. Диссертация состоит из введения, 5 глав, списка использованных источников из 116 наименований и содержит 144 страницы основного текста, 81 рисунок, 19 таблиц, 5 приложений.
Во введении обосновывается актуальность рассматриваемой темы, приведены общая характеристика работы и основные положения, которые автор выносит на защиту.
В первой главе дается обзор развития и современного состояния теории сейсмостойкости. Существенный вклад в современное состояние теории сейсмостойкости внесли Я.М. Айзенберг, B.C. Беляев, А.Н. Бирбраер, М. Био, Дж. Блюм, И.В. Гольденблат, К.С. Завриев, Г.Н. Карцивадзе, И.В. Корчинский, A.M. Масленников, C.B. Медведев, А.Г. Назаров, Ш.Г. Напетваридзе, H.A. Ни-галаенга, Н. Мононобе, Ф. Омори, B.C. Поляков, Э. Розенблюэт, Ю.Л. Рутман, O.A. Савинов, A.M. Уздин, Т. Хаузнер, Э.И. Хачиян, С.Г. Шульман и другие.
Вторая глава посвящена методам расчета на сейсмическое воздействие. Выполнен вручную расчет 5-этажной железобетонной рамы по СНиП. Результаты ручного счета сравниваются с таковыми, полученными с помощью ПК ЛИРА и SCAD (табл. 1). Анализ результатов показал, что ЛИРА дает результаты более близкие к официальному расчету. Поэтому эта программа далее используется как вспомогательная для дальнейших расчетов.
Макс, перемещение 5-го этажа, м Расчет вручную Расхождение, %
Расчет по SCAD 0,0546 0,0756 27,8
Расчет по ЛИРА 0,0779 0,0756 3
Для решения поставленных задач составлена специальная программа в символах вычислительной системы Ма&АВ. Эта же система использована и для составления всех последующих программ. Программа представляет численный метод решения системы дифференциальных уравнений движения со-
средоточенных масс. В общем виде уравнения при сейсмическом воздействием имеют вид.
Щ + СУ, + КГ, = -[Л/^]Г0(4 (1)
Здесь М, С, К- матрицы масс, демпфирования и жесткости;У, , У, ,У, -векторы перемещений, скоростей, ускорений масс системы с конечным числом степеней свободы; —матрица квазистатических перемещений (при консольной схеме здания эта матрица будет единичной); У0 (/) - амплитуды заданных ускорений.
При прямом интегрировании уравнения (1) реализуется численная итерационная процедура по времени. Чтобы найти решение уравнений (1) на каком-то интервале времени временной отрезок делится на конечное число интервалов А* и равновесие рассматривается в дискетных точках данного отрезка. Значения векторов для каждого последующего интервала времени вычисляются с помощью результатов, полученных на предыдущем шаге времени.
Метод будет безусловно устойчивым, если коэффициенты интегрировании будут находиться в пределах ае
где а = а7 /А/; р = а8 /(Л/)2 совершенно очевидно, что оптимальная устойчивость и наибольшая точность зависит от величин а и р от решаемой задачи. Эти вопросы рассматривались в диссертации специально. На основе численных экспериментов Н.МНьюмарк показал, что при решении линейных задач оптимальную устойчивость и наибольшую точность обеспечивают граничные значения а и Р, т. е. а = 0,5 и р = 0,25.
Для иллюстрации предлагаемого численного метода на рис. 1 изображены расчетные схемы рам, а на рис. 2, 3, 4 приведены развернутые во времени колебания масс верхних этажей от фрагментов трех акселерограмм в 8 баллов.
Для оценки достоверности численного метода составлена программа аналитического решения этой же задачи. В предлагаемом аналитическом методе сейсмическое воздействие в виде акселерограммы аппроксимируется трапе-цоидальными импульсами (рис. 5). Задача решается в линейной постановке.
— ——— ——— Рис. 2. Перемещения в 5-этажной раме от Рис. 1. Расчетные схемы различных акселерограмм.
1 1 2/ 3 ' V/ 4 Время, с
Рис. 3. Перемещения в 10-этажной раме от различных акселерограмм.
Рис. 4. Перемещения в 18-этажной раме от различных акселерограмм.
1,2,3 - Перемещения верхнего этажа от Карпатского землетрясения, при землетрясении в США 1954, при землетрясении Эль-Центро.
В отличие от численного интегрирования дифференциальных уравнений движения сосредоточенных масс, предложенный метод требует существенно меньше вводимой исходной информации и времени счета.
Дифференциальное уравнение движения сосредоточенных масс при кинематическом воздействии имеют следующий вид
Учет затухания колебаний осуществляется только одним членом ехр(- а; 7).
Решение системы уравнений находится в виде:
где со - спектральная переменная.
Правые части уравнений с помощью преобразования Фурье заменяются спектральными функциями:
Далее вычисляются вторые производные У^ (') от (3) и подставляются вместе с (4) в (2). После приведения подобных членов и вынесения общего множителя ехр(ко/)сЛо под знаком интеграла остается система линейных уравнений относительно неизвестных функций Уу (со):
11-8,,/^со 1г,(ш)-812т2ю У2((о)-. -д1пт„(й У„(и) = Г,(со),
-З^ш^са)-!- 1-522т2со2 ]у2(со)-. -52итлсо2у„((о)= F2(<!¡)>
-8п1т1(о Г,(со)-6я2т2со Г2(со)-. + ^1-5пятлсо2у„(со)= ^„(ш)
Неизвестные функции определяются по правилу Крамера: Г,(ш) = £,//);
Для произвольного значения ] - к
, Л ^2 4 , ч„ 2(л-2) , . 2(л-1)
, X 2 4 , , 2(л-2) , ч. 2(л-|)'
Здесь 7=1,2. л, кроме (5)
В выражении (5) множители С,-,2(л-1) ПРИ и равны соответствующим алгебраическим дополнениям А^ определителя а1п. Множители ^у,2(л-2) при со2'"-2' (степени только четные) вычисляются следующим образом: А]к принимаются за исходные определители; для имеющихся в этих определителях элементов с одинаковыми индексами составляется сумма алгебраических дополнений. При этом знаки берутся с учетом уменьшенного размера определителя и учитываются знаки исходных определителей. Эта процедура продолжается вплоть до множителя при со4. При вычислении алгебраических дополнений берутся только неодинаковые члены.
Перед интегрированием удобно разделить числитель и знаменатель выражения для УДю) на а2п и ввести обозначение Ьй = \/а2п , а затем умножить
числитель и знаменатель на (-1) . Знаменатель можно представить в виде произведения, т. е.
Здесь а>ч — корни полинома, находящегося в знаменателе.
Для дальнейшего решения необходимо уточнить вид воздействия. Как уже отмечалось, акселерограмма состоит из трапеций, как правило, «перевернутых» (см. рис.5). В этом случае
. ехр(-/С0Т]) ехр(- гшт,) 1
Рис. 5. а - произвольная акселерограмма; б - фрагмент акселерограммы в форме трапецеидального импульса
Расчетная схема при сейсмическом воздействии (учитываются только горизонтальные колебания) принимается такой, что на все массы воздействие оказывается одинаковым. В этом случае все ЬкА будут равны единице. Подстановка значения (8) в (5), а затем в (6) и интегрирование по вычетам приводит к окончательным результатам. Для времени ^ < т,
(ехр(- )э1п Юу/, - ехр(- ))
Здесь¿2 -угловые частоты свободных колебаний
конструкции. В формулах введены обозначения:
g = 1,3. (2п - 3) - нечетные числа; <7=1,2. п, кроме к;
Изложенным путем вычисляются перемещения. Силы инерции находятся как обычно:
По приведенным формулам составлен алгоритм и программа для ПК. В качестве примера рассмотрена рама четырехэтажного здания с учетом основания. На рис. 6 и 7 представлены результаты сравнения перемещений аналитического решения 1-го и 4-го этажей здания с учетом основания при воздействии Карпатского землетрясении с результатами численного решения. Сравнение графиков позволяет сделать вывод о достоверности результатов численного решения, так как ординаты графиков получились одного порядка.
Рис. 6. Перемещение 1-го этажа. 1- Расчет Рис. 7. Перемещение 4-го этажа. 1- Расчет численным методом. 2- Расчет численным методом. 2- Расчет
аналитическим методом. аналитическим методом.
Глава 3. В главе 2 было показано, (впрочем, это отмечается и в технической литературе) расхождение результатов расчета по СНиП с результатами расчета с использованием фрагментов акселерограмм. Естественно возникает вопрос о соответствии расчетов по СНиП более строгим расчетам с использованием акселерограмм. Для установления соответствия расчетов по СНиП расчетам с использованием акселерограмм выбраны три рамы различной жесткости. В начале рамы рассчитаны по СНиП с помощью ПК ЛИРА; расчет по этой программе совпадет с ручным расчетом, как показано в главе 2. Затем рамы рассчитывались с использованием акселерограммы с переменным средним периодом.
За основу исследований взята структура акселерограммы Калифорнийского землетрясения 1949 г. По Российской шкале это землетрясение относится к 7 баллам. Число ординат амплитуды ускорений 231. При Д* = 0,023 с. акселерограмма за ? =5,313 с. имеет 16,6 циклов (средний период колебаний
Г = 0,32 е.). Наиболее опасны для зданий частоты колебаний от 0,05 до 10 Гц. Меняя средний период акселерограммы от 0,05 до 10 Гц, определяется Д/. Например, для 0,1 Гц время действия акселерограммы / = 0,1-16,6 = 1,66 с. Д/ = //ЛГ = 1,66/231 = 0,00719 е., для 10 Гц М = 0,719 с.В качестве оценки величины воздействия принято перемещение массы верхнего этажа.
Расчеты выполнены с помощью программы численного метода. Результаты вычислений приведены на графиках (рис. 8, 9). Графики для каждой рамы построены в относительных координатах. По горизонтальной оси отложены соотношения периодов Т/Т, а по вертикальной оси - отношение перемещений 2™" / 20. Здесь Т0- основной период колебания конструкции; 20-перемещение массы верхнего этажа, определяемый по СНиП. Т. и 2™*- соответственно средние периоды используемых акселерограмм и максимальные перемещения верхнего этажа, получаемые с их помощью.
На уровне 1 (см. рис. 8) показан результат расчета по СНиП. Квадратиками показаны результаты, получаемые с использованием акселерограмм: 1 - Калифорнийское землетрясение 1949 года; 2 - Карпатское землетрясение 1977 года; 3 - Землетрясение Эль Центро (1940 г.); 4 - Землетрясение в США (1954 г.). Положение этих точек находится следующим образом: на горизонтальной оси откладывается соотношение периодов колебаний, как указано выше, а на вертикальной оси - отношение полученного с помощью акселерограммы результата к результату по СНиП. При этом учитывается балльность акселерограммы, т. е. полученное значение делится на результат, полученный по СНиП соответствующей балльности. Поскольку в основу расчета положена акселерограмма Калифорнийского землетрясения, то точка 1 легла на кривую графика. График позволяет определить процент расхождения результатов по СНиП и по акселерограммам. Графики по сути представляют собой огибающую воздействия акселерограмм с периодами от 0,07 с до 7,0 с.
Рис. 8. Зависимость перемещений от периодов для 5-этажной рамы
Расчет по СНиП 7 баллов = 0,1136 м Г= 1,7165 с. V = 0,25
/Расчет п / 7 баллов 1А( Т-4,151 ¥=1'02 J СНиП г18= 1,01 м 4 с. . . .
О 0,4 0,8 1 1,2 1,4 Отношение периодов колебаний
2 3 4 Отношение периодов колебаний
Рис. 9. Зависимость перемещений от периодов:
а - для 10-этажной рамы, б - для 18-этажной рамы
Анализ графиков на рис. 8,9, дает основание сделать следующие выводы о расчете по спектральной теории (по СНиП).
1. Для «жестких» сооружений (рис. 8) расчет по СНиП представляет как бы средние значения результатов расчета по фрагментам реальных акселерограмм, так как значения последних могут быть как больше, так и меньше значений результатов по СНиП. Видимо, спектральный метод и разрабатывался для такого типа малоэтажных зданий. Но, в отличие от расчета по акселерограммам, СНиП не улавливает возможные для этих зданий резонансы.
2. Для зданий «средней» жесткости (рис. 9, а) спектральный метод почти полностью обеспечивает надежность зданий. Резонансные явления практически исключаются.
3. Совершенно не приспособлен спектральный метод для расчета «гибких» сооружений (рис. 9, б). Для таких зданий и сооружений метод дает большой запас прочности. Резонансные явления здесь полностью исключены. Впрочем, по СНиП расчет зданий более 16 этажей рекомендуется производить с использованием акселерограмм.
В этой же главе приведено исследование влияния коэффициента поглощения на максимальное перемещение здания при широком изменении коэффициента поглощения.
Глава 4 посвящена учету упругого основания и сейсмоизолирующему скользящему поясу. Исследование проводится на частном примере 10-этажной рамы.
Для учета этого явления между грунтом и стенкой фундамента вводятся упругие связи, учитывающие податливость грунта. Грунт представляется моделью Винклера с одним коэффициентом постели и с двусторонней непрерывной связью. Для учета массы фундамента вводится дополнительная степень свободы. В итоге расчетная схема для протяженных зданий примет вид (рис. 10). Жесткость упругих связей представляется через коэффициент по-
стели к (кН/м3). При принятой расчетной схеме при приложении к фундаменту единичной безразмерной силы все здания переместятся на величину 50 = 1/к , как твердое тело. Принимается, что фундамент столбчатый с размерами в разрезе 0,8x0,2 м. В расчет вводится одна рама, фундамент шириной 1 м, отличное основание к- 20 кг/см3 = 2.105 кН/м3, или с учетом ширины фундамента.
Рис. 11. Матрица податливости с учетом Рис. 10. Расчетная схема 10-этажного здания основания
к = 2,0-ю' -2-1 = 4-Ю5 кН/м ■ Податливость б0 = 1/4.ю' = 2,5-Ю"6 кН/м ■
Вследствие увеличения числа степеней свободы изменяются размеры матриц. Если приложить /<'= 1 к фундаменту, то он сместится на величину 80. На эту же величину сместятся все массы; если приложить Г = 1 к любой из масс, то фундамент также сместится на 80. Перемещение массы будет 8^+80. На рис. 11 представлена матрица податливости с учетом основания.
В качестве объекта принята одна 10-этажная рама (рис. 10), выделенная из здания. Фундамент принят столбчатым, недеформируемым, с площадью контакта с грунтом 2 м2. Это несколько идеализированная схема. Масса фундамента принята равной 70 т. Воздействие на конструкцию принималось в виде фрагмента акселерограммы Карпатского землетрясения (по российской шкале интенсивности - 8 баллов). Условно принимается, что акселерограмма записана на грунте рядом со зданием.
На рис. 12 представлено развернутое во времени изменение изгибающего момента, возникающего внизу средней стойки (сплошная линия - без учета основания; штриховая линия - для к = 2,0 кг/см3; пунктирная линия -для к = 8,0 кг/см3). Совершенно очевидно, что податливость грунта для принятой расчетной схемы в виде консоли не оказывает заметного влияния на усилия. На рис. 13 показан график колебаний массы 10-го этажа при тех же условиях, что и на рис. 12. Этот график также указывает на незначительное влияние податливости основания на колебания зданий, естественно, если грунт неводонасыщенный.
Время, с Рис. 12. Изменение момента в стойке
Рис. 13. Перемещение 10-го этажа от Карпатского землетрясения
В этой же главе предлагается расчет систем с сейсмоизолирующим скользящим поясом на воздействие фрагментов реальных акселерограмм. Система непосредственно соединена с нелинейным демпфером с жестокостью г