автореферат диссертации по строительству, 05.23.17, диссертация на тему:Напряженно-деформированное состояние конструкций, взаимодействующих с нелинейно-деформируемой средой

На правах рукописи

ШАШКИН Всеволод Алексеевич

НАПРЯЖЕННО-ДЕФОРМИРОВАННОЕ СОСТОЯНИЕ КОНСТРУКЦИЙ, ВЗАИМОДЕЙСТВУЮЩИХ С НЕЛИНЕЙНО-ДЕФОРМИРУЕМОЙ СРЕДОЙ

Специальность: 05.23.17 - Строительная механика

АВТОРЕФЕРАТ

диссертации на соискание учёной степени кандидата технических наук

1 9 СЕН 2013

005533107

Санкт-Петербург 2013

005533107

Диссертация выполнена в ФГБОУ ВПО «Санкт-Петербургский государственный архитектурно-строительный университет» на кафедре строительной механики

Научный руководитель: кандидат технических наук, профессор

Бабанов Владимир Владимирович

Официальные оппоненты: Глаговский Вячеслав Борисович

доктор технических наук, ОАО «Всероссийский научно-исследовательский институт гидротехники имени Б. Е. Веденеева», Санкт-Петербург, 1-й заместитель генерального директора — научный руководитель;

Улитин Виктор Васильевич

доктор технических наук, профессор каф. механики и прочности Санкт-Петербургского национального исследовательского университета информационных технологий, механики и оптики

Ведущая организация: ФГБОУ ВПО «Санкт-Петербургский

государственный политехнический университет»

Защита состоится 17 октября 2013 г. в 14:00 часов на заседании диссертационного совета Д 212.223.03 при ФГБОУ ВПО «Санкт-Петербургский государственный архитектурно-строительный университет» по адресу: 190005, г. Санкт-Петербург, 2-я Красноармейская ул., д. 4, зал заседаний (ауд. 219). Тел./факс: 8 (812) 316-58-72 Email: rector@spbgasu.ru

С диссертацией можно ознакомиться в фундаментальной библиотеке ФГБОУ ВПО «Санкт-Петербургский государственный архитектурно-строительный университет».

Автореферат разослан » С^Щ^ 2013 г.

Учёный секретарь диссертационного совета д.т.н., проф. Кондратьева Лидия Никитовна

I. ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность исследования. Учет совместной работы оснований и сооружений на протяжении многих десятилетий является основополагающим принципом расчета и проектирования в отечественных нормах. С вступлением в силу в 2010 г. 384-ФЭ «Технический регламент о безопасности зданий и сооружений» этот принцип получил статус положения Федерального закона. СП «Основания зданий и сооружений» предписывает рассматривать сооружение и его основание с учетом их взаимодействия, а нагрузки и воздействия на основания (которые СП «Нагрузки и воздействия» относит к длительным), устанавливать расчетом, исходя из рассмотрения совместной работы сооружения и основания. Аналогичное положение имеется и в СП «Свайные фундаменты», который предусматривает, что расчет системы «свайное основание - изгибаемая плита - верхнее строение» следует производить в пространственной постановке с учетом совместной работы надземных и подземных конструкций зданий, свайных фундаментов и их основания.

До недавнего времени учет взаимодействия зданий и оснований представлял собой чрезвычайно сложную задачу, решение которой можно было осуществить только на основе упрощающих предпосылок. Происходящие в последнее 10-15 лет процессы бурного развития вычислительной техники позволяют сегодня реализовать высказанное еще 40 лет назад М.И. Горбуновым-Посадовым утверждение о том, что будущее методов расчета взаимодействия основания и сооружения заключается в применении теории упругости совместно с теорией пластичности, в преодолении трудностей практических расчетов при помощи программ для ЭВМ. В настоящее время появилась возможность освободиться от неизбежных ранее упрощений и моделировать трехмерные задачи расчета конструкций на нелинейно деформируемой среде.

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

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

Степень разработанности темы исследования. Учет взаимодействия оснований и сооружений осуществлялся по двум основным направлениям. В первом из них взаимодействие конструкций с основанием моделировалось упругими опорами, податливость которых описывалась «коэффициентами постели». В основе этого направления лежит гипотеза Винклера (однопараметри-ческая модель упругого основания). Для учета распределительной способности основания были разработаны двухпараметрические (П.Л.Пастернак, В. 3. Власов и др.), а также трехпараметрические модели (В. А. Барвашов, В. Г. Федоровский). Были предложены также вариации модели Винклера с переменной по площади фундамента величиной коэффициента постели (В.Г.Федоровский).

Эти модификации позволили приблизить указанные методы ко второму направлению, которое заключалось в применении к описанию работы основания теории упругости. Рассмотрению работы гибких фундаментов на упругом полупространстве посвящены работы М. И. Горбунова-Посадова, И. А. Симвулиди, Б. Н. Жемочкина, А. Н. Синицына, В. А. Флорина, С. Н. Клепикова, Ю. К. Зарец-кого и др. На западе развитие этих методов связано с трудами М. Био, Е. де Беера, X. Грасгоффа и др. Для решения практических задач был принят ряд упрощений, которые еще полвека назад позволили рассчитывать гибкие фундаментные балки и плиты на упругом основании с помощью имеющихся средств вычислительной техники.

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

Возможность применения в проектной практике расчетов, основанных на решении смешанной задачи теории упругости и пластичности, позволяющей приблизить очертание контактной эпюры к реально наблюдаемой на практике, появилась с развитием вычислительной техники. При этом преимущественное развитие получили расчетные программы, ориентированные на решение геотехнических задач. Среди них - идеально упруго-пластические модели, развитие которых связано с именами А. К. Бугрова, А. Б. Фадеева и др.; «шатровые» модели (Д. Друкер, В. Прагер, К. Роскоу, Дж. Берланд и др.); модели «с двойным упрочнением» (П. Вермеер, Т. Бенц, Ю. К. Зарецкий). К сожалению, западные разработчики компьютерных программ совместного расчета здания и основания традиционно отдавали предпочтение упрощенным моделям, построенным на использовании «коэффициентов постели».

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

Цель и задачи исследования.

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

Объектом исследований является напряженно-деформированное состояние здания с несущими стенами, формирующееся в результате его взаимодействия с нелинейно деформируемой средой основания.

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

Задачи исследования:

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

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

3. Установить влияние различных факторов на величину предельной нагрузки на штамп на упруго-пластическом полупространстве и на деформацию штампа в допредельном состоянии.

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

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

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

Область исследования соответствует требованиям паспорта научной специальности ВАК: 05.23.17 - Строительная механика, а именно: содержанию

специальности, каковым являются методы расчета сооружений и их элементов на прочность, устойчивость при силовых воздействиях, а также следующим основным направлениям: п. 1 «Общие принципы расчета сооружений и их элементов», п. 2 «Линейная и нелинейная механика конструкций и сооружений, разработка физико-математических моделей их расчета», п. 4 «Численные методы расчета сооружений и их элементов», п.8 «Исследование нагрузок на сооружения».

Научная новизна исследования

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

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

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

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

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

Практическая ценность и реализация результатов исследований.

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

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

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

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

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

Предложенный подход к учету взаимодействия надземной конструкции и нелинейно-деформируемой среды основания получил внедрение в расчетах при разработке проектов нового строительства (жилой 15-этажный дом на Крас-нопутиловской ул. 113, жилой 16-этажноый многосекционный дом в пос. Лав-рики) и проектов усиления деформированных зданий в Санкт-Петербурге (комплекс разновысотных жилых зданий на Варшавской ул. д.59; 9-этажные общежития типовой серии 1-447-С-54, к которой относится обрушившееся здание на Двинской ул.). На основании выполненных исследований разработан раздел Справочника Проектировщика «Основания, фундаменты и подземные сооружения», посвященный рекомендациям по учету эффектов взаимодействия оснований и сооружений (готовится к печати).

Достоверность научной гипотезы обеспечивается применением сертифицированной расчетной программы, в которой реализованы расчетные модели, верифицированные на соответствие аналитическим решениям и результатам наблюдений за напряженно-деформированным состоянием системы «основание-сооружение» на статистически представительной выборке реальных объектов; соответствием результатов численного анализа по предложенной методике и фактической картины развития деформаций на реальных объектах.

Апробация работы. Основные положения и результаты диссертации докладывались на научных семинарах и конференциях, проводимых ассоциацией инженерных изысканий в строительстве в г. Москва (2011), международным обществом по механике грунтов и геотехническому строительству в Санкт-Петербурге (2012) и г. Таллине (2012), на конференции молодых ученых, проведенной Санкт-Петербургским государственным архитектурно-строительным университетом (2012), на секции строительной механики Петербургского Дома ученых (2013), на семинаре в Петербургском государственном университете путей сообщения (2013).

Публикации. Материалы диссертации опубликованы в 8 печатных работах, в том числе 4 в изданиях, входящих в перечень ведущих рецензируемых научных журналов, утвержденных ВАК Российской Федерации.

Структура и объем работы: диссертация состоит из введения, 4 глав, основных выводов, списка литературы и 2 приложений. Она имеет объем 172

страницы печатного текста, включая 91 рис. и 11 таблиц. Список литературы включает 124 наименования, в том числе 32 на иностранных языках.

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

II. ОСНОВНЫЕ ПОЛОЖЕНИЯ И РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЙ ДИССЕРТАЦИИ, ВЫНОСИМЫЕ НА ЗАЩИТУ

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

На пути перехода задачи расчета деформаций полупространства в задачу предельного равновесия имеются теоретические и практические проблемы. Одна из них, обозначенная С.Слоаном, М.Рандольфом, А.Б.Фадеевым и др., связана с так называемым «запиранием» упруго-пластического решения. Для того, чтобы избежать «запирания», используют элементы более высокого порядка (например, в программе РЬАХ18).

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

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

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

КХ!

где К - матрица жесткости системы, Рех - вектор внешних сил, АР = КХкЛ - Рт ) - добавка в вектор правых частей, вычисляемая по разнице приращений «упругих» и допускаемых напряжений в элементе в соответствии с моделью материала.

Определим напряжения в элементе в соответствии с физической моделью ат (х) = оу + исостоящие из объемных напряжений и девиатора напряжений. Вычислим «упругие» напряжения о" = ДВХ , которые также состоят из объемных напряжений и девиатора сге = а'у + <тео. Вычислим вектор «упругих» правых частей для элемента Ре = КХ, где К-матрица жесткости элемента. Вектор упругих правых частей разделим на «объемные» и «девиаторные» компоненты:

Р£=ВТауеУ, Ре0=Ре-Р£.

Вектор сил, соответствующий физической модели, можно выразить, разделив силы на «объемные» и «девиаторные» составляющие:

Рт(Х) = Ру + Рв = ВТауУ + Вта°У • Вектор сил выражается следующим образом:

^ IIо II

При использовании выражения (1) вектор сил в соответствии с физической моделью вычисляется непосредственно по вектору «упругих» сил, в результате чего информация о виде деформации не теряется.

Решение тестовой задачи (рис. 1), демонстрирует, что условие (1) позволяет избавиться от проблем с «запиранием» решения при использовании элементов 1-го порядка и получить корректное решение, соответствующее расчету по различным программам с использованием элементов 2-го порядка. Как видно из рисунка 1, при решении задачи без использования выражения (1) некорректной является не только предельная нагрузка на фундамент, но и весь график «нагрузка-осадка» в допредельном состоянии.

В различных конечно-элементных программах используются модели с разными поверхностями текучести, близкими к классическому критерию Кулона-Мора (рис. 2, а). Например, в модели Друкера-Прагера вместо шестигранной пирамиды Кулона-Мора используется коническая поверхность.

По результатам решения конечно-элементных задач оказывается, что коническая поверхность, описанная вокруг пирамиды Кулона-Мора, приводит к существенно большей несущей способности, чем по классическому критерию

Кулона-Мора (рис. 2, б). С другой стороны, коническая поверхность, вписанная в пирамиду Кулона-Мора, не соответствует критерию прочности при сг2 = <х3, т. е. при обычном напряженном состоянии в трехосных испытаниях.

Р. кПа

Теория

" 1 1 порядок

—Z-Soil 2 порядок

—— PLAXIS 2 порядок

9 Midas 2 пооядок

-FEM models 1 порядок,

условие (1)

Рис. 1. Графики «нагрузка-осадка» для задачи о вдавливании штампа при использовании элементов 1-го и 2-го порядка.

Р. кПа

900 1100 —Усло1и«ч=Мр

—— Мрдрл»

Друквр1-Лр»гвр1 Усложи« Кулона. Море

£.11

Рис. 2. Положения предельных поверхностей текучести (а) и графики «нагрузка-осадка» (б): 1 - по условию Кулона-Мора, 2 - по выражению <у||т = Мр, где р = сг, + <уг + <х3, М - константа, 3 - по модели Друкера-Прагера.

Таким образом, при решении задач о взаимодействии конструкций с упруго-пластическим полупространством использование модели с предельной поверхностью, описываемой критерием Мизеса-Шлейхера-Боткина (Друкера-Пра-гера), может приводить к завышению предельной нагрузки и занижению деформаций в допредельном состоянии.

Устранение недостатков использования элементов 1 -го порядка позволяет получить коэффициенты Ny, рассчитанные с помощью программы FEM models,

достаточно близкие к решениям С.Слоана, С.Мартина и В.Г.Федоровского. При этом следует отметить, что результаты для малозаглубленного фундамента существенно отличаются в меньшую сторону от предельной нагрузки, рассчитанной по действующим нормам, что необходимо учитывать при расчетах и проектировании.

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

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

Решение данных задач выполнялось для идеально связной среды (с нулевым углом внутреннего трения). В соответствии с классическими решениями учет заглубления штампа выполняется через добавление по поверхности нагрузки q, равной давлению грунта выше подошвы заглубленного фундамента. Однако в реальности при увеличении заглубления начинает играть роль не только при-грузка основания рядом с фундаментом, но и прочность грунта выше подошвы штампа, что приводит к увеличению предельной нагрузки по сравнению с теоретическим решением (рис. 3, а), а) б)

р . кПа I', кПа

Рис. 3. Зависимость предельного давления от глубины заложения фундамента сI (а) и от эксцентриситета е приложения нагрузки (б) для фундамента шириной 2 м (с = 20 кПа. <р= 0°)

Как показывает численный эксперимент, при небольшой глубине фундамента происходит выпор по форме, близкой к незаглубленному фундаменту. При увеличении заглубления фундамента область выпора становится менее определенной (рис. 4, а), а зона линейного поведения на графике «нагрузка-осадка» увеличивается.

Учет эксцентриситета приложения нагрузки в отечественных нормах выполняется путем уменьшения ширины рассчитываемого фундамента Ь'= Ь- 2е, где е -эксцентриситет приложения нагрузки. Это вполне оправдано при небольшом эксцентриситете (рис. 3, б). Однако при приближении точки приложения нагрузки к краю фундамента такой подход приводит к снижению до нуля несущей способности.

е=0.1м

шшш

1..... е=0.6м

ЩИ^М^^^^^Р^^ -----1

1 е=1м

Рис. 4. Характер областей сдвига (максимальные перемещения на стадии потери устойчивости) при различном заглублении И (а) и эксцентриситете е (6) штампа шириной Ь = 2 м

(с = 20 кПа, <р=0°).

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

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

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

Рис. 5. Верификация упруго-пластической модели с независимым упрочнением и критерием прочности Кулона-Мора, реализованной в программе FEM models, на соответствие результатам натурных испытаний свай (а) и штампов (б) статическими нагрузками: 1 - численное моделирование испытаний.

Как показывают результаты численного эксперимента, выполненного для параметров упруго-пластического полупространства, соответствующих, например, инженерно-геологическим условиям южных районов Санкт-Петербурга, проявление эффекта вовлечения в работу плиты в составе свайного фундамента с низким плитным ростверком заметно в интервале расстояний между сваями в его составе от 1 Od до 16d, где d - размер поперечного сечения сваи (рис. 6).

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

При расстоянии между сваями, равном или менее 8d, несущая способность свайного фундамента с низким ростверком не отличается от несущей способности соответствующего свайного фундамента с высоким ростверком (рис. 7), а при расстоянии, равном или превышающем 24 d - от несущей способности

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

Взаимодействие нелинейно деформируемого полупространства и свай, проявляется по-разному в зависимости от шага свай. Существует такой малый шаг свай (в численном эксперименте он составлял Зс1), при котором сваи и среда в межсвайном пространстве работают как единый массив, близкий к жесткому штампу, образуя условный фундамент с подошвой на уровне нижнего конца свай. Такой фундамент обладает наибольшей несущей способностью в ряду свайных фундаментов с низким плитным ростверком и различным шагом свай.

(а) - с низким ростверком, (б) - с высоким ростверком при различном шаге свай

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

п < _ сч_

где п - шаг свай, Ысе - несущая способность одиночной сваи, И - длина сваи, с1-диаметр сваи, си - прочностная характеристика среды основания, удельный вес грунта выше подошвы условного фундамента.

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

а)

S.M

б)

В)

Нагрузка, кН

Нагрузка, кН

0 Vv™^"300 ivm ото гио о зооо 4000

m\ \ "ч, -o.í

\\ \ i

и\\ \ -0.4 ^V^V \ °

8 \ м и\\ \ \ -0.5 -Об «\Л аа \\ \ __о \\ \ VI п

О—ип \ 5 \\\\ \ ос ез • -0.7 s \ \ s \\ \ I

I I о -О.в 3 \\\ -а—\\\ I

........ее Я £ 1 — 1 09 = \\\\

5 II & , & \\\|

S.M

S.M

Рис. 8. Распределение усилий в сваях под штампом с шагом свай-, a) 3d, б) 6d, e)12d

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

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

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

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

При этом рассматривались только те сочетания, которые могут встречаться в реальности. Моделировались конкретные конструктивные схемы протяженных зданий и «точечных» зданий, охватывающие основные типы, встречающиеся на практике: с продольными несущими стенами; поперечными несущими стенами; с сеткой продольных и поперечных несущих стен. Поскольку эффекты совместного расчета, как известно, проявляются в пределах нескольких нижних этажей, рассчитывались здания высотой 10 этажей с их поэтажным возведением.

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

Для выявления эффектов взаимодействия конструкций и упруго-пластического полупространства результаты расчетов для каждого случая сравнивались с результатами соответствующего расчета на абсолютно неподатливом основании (на жестких опорах), что эквивалентно сбору нагрузок по грузовым площадям, практикуемому при проектировании. Такое сравнение представляется вполне показательным, поскольку наиболее полно выявляет исследуемый эффект, совершенно отсутствующий при рассмотрении конструкции на жестких опорах. Отношение усилия в некоторой точке конструкции здания при решении на податливом основании, к усилию в той же точке при решении на жестких опорах, назовем показателем эффекта взаимодействия конструкции и упруго-пластического полупространства К№[.

Для каждой конструктивной схемы построены графики зависимости К от соотношения модулей деформации Етнстр/Есекущ (модуля деформации материала надземных конструкций к секущему модулю деформации нелинейно-деформируемого полупространства) (рис. 9).

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

Характерно, что на всех рассматриваемых графиках отсутствует какой-либо начальный горизонтальный участок. Следовательно, не существует таких соотношений Етжтр/Есегущ, при которых эффект взаимодействия не имеет значения в аспекте концентрации усилий в контактной зоне. Во всех случаях, кроме расчета зданий на недеформируемом скальном основании в конструкциях здания происходит концентрация усилий.

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

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

Еконстр Естщ

Рис. 9. График зависимости показателя эффекта взаимодействия конструкций и нелинейно деформируемой среды основания от соотношения их жесткостей для различных конструктивных схем зданий

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

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

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

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

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

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

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

ными стенами - соответственно к 4-х к 3-кратной концентрации по сравнению с расчетом на жестких опорах).

узлы КЭ сетки уиыКЭсспш узлыКЭссткн

Рис. 11. Эпюры распределения вертикальных усилий по низу торцевой стены: а - при поэтажном возведении здания с поперечными несущими стенами, б — при возведении того же здания на свайном фундаменте (1 - при опирании свай на среду малой степени деформативности, 2 - на среду средней степени деформативности); в, г - то же, для здания с продольными и поперечными несущими стенами.

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

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

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

Рис. 12. Проявление эффекта взаимодействия реального здания и среды основания: а - вид кирпичного здания серии ш-5733 на пр.Тореза в Санкт-Петербурге; 6- КЭ модель здания; в - усиление стен обоймами; г - эпюра распределения вертикальных усилий по низу торцевой стены; д - распределение вертикальных усилий в стенах при расчете на упруго-пластическом полупространстве; г - то же, на жестких опорах.

ЗАКЛЮЧЕНИЕ

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

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

3. Установлено, что существенное влияние на характер контактного взаимодействия в задаче о штампе на упруго-пластическом полупространстве в допредельном и предельном состоянии оказывает вид предельной поверхности. Модели с предельной поверхностью, описываемой критерием Мизеса-Шлейхе-ра-Боткина (Друкера-Прагера), отличаются от классического критерия Кулона-Мора и приводят к более высокой несущей способности фундамента, чем по решению теории предельного равновесия.

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

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

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

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

Список работ, опубликованных по теме диссертации В изданиях, рекомендованных ВАК РФ:

1. Шашкин, В.А. Эффекты концентрации напряжений в конструкциях здания при взаимодействии с основанием / В.А. Шашкин // Жилищное строительство. - 2012, №9. - С. 9-14.

2. Шашкин, В.А. Расчетный анализ работы свайных фундаментов с низким и высоким ростверками с учетом нелинейной работы основания / В.А. Шашкин,

B.В. Бабанов // Основания, фундаменты и механика грунтов . - 2012, №2 -

C. 2-7.

3. Шашкин, В.А. Взаимодействие зданий на разнотипных фундаментах и основания / В.А. Шашкин // Жилищное строительство. - 2012, №5. - С. 38-42.

4. Шашкин, В.А. Решение задач предельного равновесия с использованием метода конечных элементов / В.А. Шашкин, К.Г. Шашкин // Жилищное строительство. - 2013, №4. - С. 47-50.

В других изданиях:

5. Шашкин, В.А. Конечно-элементный анализ неравномерных осадок здания, взаимодействующего с нелинейно-деформируемой средой / В.А. Шашкин // Актуальный проблемы строительства и архитектуры, Часть II. - 2012,- С. 90-92.

6. Шашкин, В.А. Численное моделирование задач предельного равновесия с помощью упруго-пластической модели / В.А. Шашкин, К.Г. Шашкин, М.В. Дунаева // Геотехника. - 2011, №4 - С.10-23.

7. Шашкин, В.А. Эффекты взаимодействия оснований и сооружений / В.А. Шашкин // Развитие городов и геотехническое строительство. - 2012, №14-С. 141-167.

8. Шашкин, В.А. Решение задач предельного равновесия с использованием метода конечных элементов / В.А. Шашкин, К.Г. Шашкин, М.В. Дунаева // Развитие городов и геотехническое строительство.-2011,№ 13.-С. 81-95.

Компьютерная верстка И. А. Я блоковой

Подписано к печати 29.08.13. Формат 60x84 1/16. Бум. офсетная. Усл. печ. л. 1,4. Тираж 120 экз. Заказ 93. Санкт-Петербургский государственный архитектурно-строительный университет.

190005, Санкт-Петербург, 2-я Красноармейская ул., 4. Отпечатано на ризографе. 190005, Санкт-Петербург, 2-я Красноармейская ул., 5.

Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования Санкт-Петербургский государственный архитектурно-строительный университет

на правах рукописи

04201363260

'пси*—

Шашкин Всеволод Алексеевич

НАПРЯЖЕННО-ДЕФОРМИРОВАННОЕ СОСТОЯНИЕ КОНСТРУКЦИЙ, ВЗАИМОДЕЙСТВУЮЩИХ С НЕЛИНЕЙНО-ДЕФОРМИРУЕМОЙ

СРЕДОЙ

Специальность: 05.23.17 - Строительная механика

Диссертация на соискание ученой степени кандидата технических наук

Научный руководитель: к.т.н., доцент, проф. каф. «Строительная механика» Ба-банов В.В.

Санкт-Петербург 2013

Введение............................................................................................................5

1 МЕТОДЫ УЧЕТА ВЗАИМОДЕЙСТВИЯ КОНСТРУКЦИЙ И ДЕФОРМИРУЕМОГО ПОЛУПРОСТРАНСТВА.......................................................................................12

1.1 Применение метода конечных элементов к расчету конструкций зданий, сооружений и оснований....................................................................................................................13

1.2 Развитие методов расчета конструкций на упругом основании.............................14

1.3 Модели работы деформируемой среды, применяемые при решении задач взаимодействия конструкций и основания...............................................................................18

1.3.1 Модель линейно деформируемой среды..................................................19

1.3.2 Идеально упруго-пластические модели....................................................20

1.3.3 Шатровые модели...............................................................................21

1.3.4 Модели с двойным упрочнением...........................................................23

1.4 Краткий обзор существующих методов решения задач предельного равновесия.............................................................................................................25

Выводы по главе 1.................................................................................................30

2 РЕШЕНИЕ ЗАДАЧ ПРЕДЕЛЬНОГО РАВНОВЕСИЯ С ИСПОЛЬЗОВАНИЕМ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕТОВ...................................................................................31

2.1. Особенности моделирования задач совместного расчета методом конечных элементов....................................................................................................31

2.2. Упруго-пластическая модель среды основания................................................34

2.3 Решение задач предельного равновесия о вдавливании штампа в упруго-пластическую среду полупространства..................................................................................40

2.3.1 Проблема «запирания» решения задач предельного равновесия методом конечных элементов с использованием упругопластических моделей...................40

2.3.2 Выявление факторов, влияющих на вид графика «нагрузка-осадка» при

решении задач о вдавливании штампа...........................................................45

2.3.3. Расчетная модель и основные результаты расчета предельной несущей способности для штампа на поверхности.......................................................48

2.3.4 Сравнение коэффициента по результатам решения упруго-пластических задач и по теории предельного равновесия....................................................49

2.3.5 Оценка несущей способности заглубленного фундамента и фундамента с эксцентрично приложенной нагрузкой.........................................................52

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

полупространстве.........................................................................................55

2.4.1 Постановка численного эксперимента...................................................55

2.4.2 Расчетная оценка несущей способности сваи, штампа на поверхности грунта и

заглубленного штампа..............................................................................57

2.4.3 Оценка несущей способности свайного фундамента.................................64

2.4.4 Анализ результатов численного моделирования штампа на сваях на упруго-пластическом полупространстве..........................................................69

2.4.5 Распределение нагрузок между сваями в составе свайного фундамента.........71

Выводы по главе 2..............................................................................................77

3 АНАЛИЗ СОВМЕСТНОЙ РАБОТЫ КОНСТРУКЦИЙ И УПРУГО-ПЛАСТИЧЕСКОГО ПОЛУПРОСТРАНСТВА..........................................................................................79

3.1 Постановка серии численных экспериментов................................................79

3.2 Проявление эффектов взаимодействия здания на фундаментной плите и упруго-пластического полупространства..........................................................85

3.2.1 Рассмотрение различных расчетных схем зданий...................................85

3.2.2 Учет поэтажного возведения здания...................................................90

3.2.3 Влияние жесткостных характеристик материала конструкции..................92

3.2.4 Учет различной жесткости нелинейно деформируемой среды..................93

3.2.5 Анализ напряженно-дсформированного состояния плитного фундамента ...94

3.3 Численное моделирование взаимодействия здания на свайном фундаменте и упруго-пластического полупространства.........................................................104

3.3.1 Влияние конструктивной схемы здания на изменение контактной эпюры

при наличии свайного фундамента............................................................105

3.3.2 Учет жесткостных характеристик материала конструкции здания на свайном фундаменте.......................................................................................108

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

3.4 Закономерности проявления эффектов взаимодействия конструкции и нелинейно деформируемой среды в зависимости от соотношения их жесткостей ....110

Выводы по главе 3.........................................................................................115

4 СЛУЧАИ ИЗ ПРАКТИКИ И РЕКОМЕНДАЦИИ ПО УЧЕТУ ЭФФЕКТА ВЗАИМОДЕЙСТВИЯ КОНСТРУКЦИЙ И НЕЛИНЕЙНО ДЕФОРМИРУЕМОЙ СРЕДЫ.......................................................................................................116

4.1 Эффекты взаимодействия конструкций и нелинейно деформируемой среды основания на примере реальных объектов......................................................116

4.1.1. Пример протяженного здания с несущими стенами..............................116

4.1.2. Пример точечного здания с несущими стенами....................................123

4.1.3. Примеры проявления эффектов взаимодействия основания и сооружения на

примере зданий серии 1-447-С-54...............................................................134

4.2 Рекомендации по учету эффектов взаимодействия здания и нелинейно

деформируемого полупространства в практическом проектировании......................149

ОСНОВНЫЕ ВЫВОДЫ....................................................................................154

Литература.....................................................................................................156

Приложение 1.................................................................................................165

Приложение 2.................................................................................................180

ВВЕДЕНИЕ

Актуальность темы диссертационного исследования

Учет совместной работы оснований и сооружений на протяжении многих десятилетий является основополагающим принципом расчета и проектирования в отечественных нормах. С вступлением в силу Э84-ФЗ «Технический регламент о безопасности зданий и сооружений» этот принцип получил статус положения Федерального закона. Закон обязывает учитывать «особенности взаимодействия строительных конструкций между собой и с основанием», «пластические и реологические свойства материалов и грунтов», а также «пространственную работу строительных конструкций» [124].

В актуализированной редакции Свода Правил «Основания зданий и сооружений» отмечается, что «сооружение и его основание должны рассматриваться в единстве, т.е. должно учитываться взаимодействие сооружения с основанием»; «нагрузки и воздействия на основания, передаваемые фундаментами сооружений, должны устанавливаться расчетом, как правило, исходя из рассмотрения совместной работы сооружения и основания» [120]. «Воздействия, обусловленные деформациями основания, не сопровождающимися коренным изменением структуры грунта» Свод Правил «Нагрузки и воздействия» относит к длительным нагрузкам [119].

Аналогичное положение имеется и в Своде Правил «Свайные фундаменты»: «сооружение и его основание должны рассматриваться совместно, т.е. должно учитываться взаимодействие сооружения со сжимаемым основанием»; «расчет конструктивной системы «свайное основание - изгибаемая плита - верхнее строение» в общем случае следует производить в пространственной постановке с учетом совместной работы надземных и подземных конструкций зданий, свайных фундаментов и их основания» [121].

В 1973 г. М.И. Горбунов-Посадов отмечал, что «учет жесткости верхнего строения представляет собой чрезвычайно сложную задачу, которая в настоящее время решается большей частью только на основе упрощающих предпосылок» [15, с.69].

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

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

появилась возможность реализовать высказанное 40 лет назад М.И. Горбуновым-Посадовым

5

утверждение о том, что будущее методов расчета взаимодействия оснований и сооружений заключается в применении теории упругости совместно с теорией пластичности грунтов. Он отмечал, что «трудности, состоящие на пути учета этих факторов в практических расчетах будут преодолены при помощи программ для ЭВМ» [16, с.5].

Действительно, уже в 80-90-х годах прошлого века упрощенные методы расчета гибких фундаментных балок и плит на упругом основании были успешно реализованы в вычислительных программах. В настоящее время, с увеличением скорости вычислений и объема памяти компьютеров появилась возможность освободиться от неизбежных ранее упрощений и моделировать трехмерные задачи расчета конструкций зданий на упруго-пластическом полупространстве. Применяя смешанную задачу теории упругости и пластичности для моделирования работы полупространства, необходимо находить такие методы численного решения задачи, которые позволяли бы рассматривать процесс перехода допредельного напряженно-деформированного состояния полупространства в предельное состояние. Численное решение такой задачи сталкивается с определенными трудностями, рассмотренными и преодоленными в настоящей работе.

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

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

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

6

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

Степень разработанности темы диссертации

Решение задач о взаимодействии оснований и сооружений развивалось по двум основным направлениям. Первое из них заключалось в том, что взаимодействие конструкций с основанием моделировалось упругими опорами, податливость которых описывалась различными моделями, использующими понятие «коэффициент постели». В основе этого направления лежит гипотеза Винклера (однопараметрическая модель упругого основания). Для учета распределительной способности основания были разработаны двухпараметрические модели (П.Л. Пастернак [46], В.З. Власов и др. [8]), а также трехпараметрические модели (В.А.Барвашов, В.Г. Федоровский [66]). Были разработаны также вариации модели Винклера с переменной по площади фундамента величиной коэффициента постели [66]. Эти модификации позволили приблизить методы с использованием «коэффициентов постели» ко второму основному направлению развития методов совместных расчетов конструкций и оснований.

Второе направление заключалось в применении к описанию работы основания теории упругости. Рассмотрению работы гибких фундаментов на упругом полупространстве посвящены работы М.И. Горбунова-Посадова [15-17], И.А. Симвулиди [52], Б.Н. Жемочкина [26], А.Н. Синицына [53], С.Н. Клепикова [37], Ю.К. Зарецкого и др. [28].

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

На западе развитие методов расчета конструкций на упругом основании связано с трудами таких известных исследователей, как М. Био [86], Е. Де Беер [92], X. Грасгофф [95] и др., хотя в практической деятельности предпочтение отдавалось методам, основанным на гипотезе Винклера.

Оба основных направления развития методов совместного расчета зданий и оснований, в силу заложенных в них допущений, не лишены известных недостатков [15]. При использовании гипотезы Винклера основание не имеет распределительной способности, т.е. деформации соседствующих с нагрузкой элементов поверхности отсутствуют. Гипотеза упругого полупространства преувеличивает распределительную способность грунта и предполагает возникновение бесконечно больших активных усилий в краевых зонах.

Возможность применения в проектной практике расчетов, основанных на решении смешанной задачи теории упругости и пластичности, позволяющей приблизить очертание контактной эпюры к реально наблюдаемой на практике, появилась в последнее время с развитием вычислительной техники. Были разработаны расчетные модели нелинейной работы среды основания, наибольшую известность среди которых получила модель Cam Clay, разработанная К. Роскоу и Дж. Берландом [87,88,112] и, позднее, модель HSM программы «Plaxis», предложенная П.Вермеером [115-117]. К сожалению, западные разработчики компьютерных программ, позволяющих рассматривать совместную работу здания и основания, традиционно отдавали предпочтение упрощенным моделям, построенным на использовании «коэффициентов постели».

Отечественная традиция применения теории деформируемого тела к описанию поведения основания получила в последнее десятилетие развитие в трудах В.М.Улицкого и его школы [61,62]. Ими была разработана программа FEM models, позволяющая рассматривать взаимодействие зданий и нелинейно деформируемой среды основания [63,82,83]. Настоящая работа является развитием предложенных ими подходов.

Объектом исследований является напряженно-деформированное состояние здания с несущими стенами, формирующееся в результате его взаимодействия с нелинейно деформируемой средой основания.

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

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

Для достижения данной цели были поставлены и решены следующие задачи:

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

• Провести верификацию численного ре