Содержание к диссертации
Введение
Теоретические основы исследования напряженно-деформированного и предельного состояний железобетонных конструкций в трехмерной постановке 35
1.1. Вариационная постановка задачи 35
1.2. Построение матриц жесткости трехмерных изопараметри-ческих конечных элементов 37
1.2.1. 20-ти узловой конечный элемент 38
1.2.2. Конечный элемент треугольной призмы 41
1.3. Определение узловых сил 43
1.3.1. Нормальное давление 44
1.3.2. Касательная нагрузка 46
1.3.3. Массовые силы 47
1.4. Вычисление напряжений в декартовых и местных осях 48
1.5. Модель железобетона с дискретно расположенной арматурой 49
1.5.1. Вводные замечания 49
1.5.2. Трехмерный конечный элемент с дискретными стержневыми подкреплениями 51
1.6. Конечно-элементная математическая модель трехмерноготела с дискретными мембранными подкреплениями 56
1.6.1. Вводные замечания 56
1.6.2. Мембранный конечный элемент - 58
1.7. Нелинейные задачи строительной механики железо-бетонных конструкций 61
1.7.1. Метод и алгоритм определения предельного состояния 61
1.8. Краткие сведения о программной реализации 70
1.9. Решение модельных задач 79
1.9.1. Модельная задача 79
1.9.2. Решение модельной задачи 81
1.9.3. Модельная задача 82
Исследование напряженного и предельного состояний опор внеклассного моста 84
2.1. Исследование напряженно-деформированного состояния стойки опоры моста в линейной и нелинейной постанов
ках от действия расчетных нагрузок 84
2.1.1. Постановка задачи 84
2.1.2. Анализ результатов расчетов в линейной постановке задачи 86
2.1.3. Анализ результатов расчетов в физически нелинейной постановке 88
2.2. Исследование напряженно-деформированного состояния системы буросваи-ростверк-стойка-ледозащитная оболочка в линейной постановке задачи с учетом работы металлических труб буросвай 92
2.2.1. Конструктивные особенности элементов опоры и особенности их силовых схем 92
2.2.2. Анализ результатов расчетов без включения металлической оболочки в силовую схему 94
2.2.3. Анализ результатов расчетов с учетом включения металлической оболочки в силовую схему 97
2.2.4. Анализ напряженно-деформированного состояния опоры в трехмерной физической нелинейной постановке 99
Моделирование напряженно-деформированно го и предельного состояний многослойных грунтов, взаимодействующих с деформируемыми элементами 102
3.1. Исследование напряженно-деформированного и предельного состояний склона у мостового перехода через реку Кама у села Сорочьи Горы 102
3.1.1. О постановках задач устойчивости грунтовых откосов по теории предельного состояния 102
3.1.2. Вариационная постановка задачи о НДС грунтовых откосов и ее численная реализация 106
3.1.3. Результаты расчетов откоса до и после разгрузки 111
3.1.4. Результаты расчета откоса после установки свай-шпонок 121
3.2. Исследование напряженно-деформированного и предель
ного состояний опоры моста через р. Архаровка с учетом
взаимодействия с грунтом 132
3.2.1. Вводные замечания 132
3.2.2. Геометрическая и конечно-элементная модели 134
3.2.3. Результаты расчетов 142
Моделирование реконструкции автодорожно-го городского моста через р. казанка на iii транс-портной дамбе в г. Казани 150
4.1. Особенности конструкции пролетных строений из сборных железобетонных блоков и их монтажа путем создания предварительных напряжений сжатия 150
4.2. Последовательность силовых схем, отвечающая конструкции и технологии монтажа пролетных строений 152
4.3. Последовательность расчетных схем, отвечающая этап-ности их монтажа 157
4.4. Анализ напряженно-деформированного состояния одного пролета балки от действия собственного веса и сил натяжения арматурных пучков 161
4.4.1. Основные определения и понятия 161
4.4.2. Постановка модельной задачи 163
4.4.3. Решение задачи по нетрансформирующейся схеме... 165
4.4.4. Решение задачи по трансформирующейся схеме 167
4.4.5. Сопоставление полей напряжений в сечениях балки мостового перехода через р. Казанку, найденных по двум рассматриваемым схемам 171
4.5. Управление полем напряжений в поперечных сечениях балки путем натяжения дополнительных арматурных пучков 172
4.5.1. Постановка задачи и ее решение 172
4.6. Исследование процесса формирования внутренних усилий и моментов 177
4.6.1. Решение задачи по нетрансформирующейся схеме... 177
4.6.2. Решение задачи по трансформирующейся схеме 184
4.6.3. Вычисление коэффициентов к , k\j) Ш
4.6.4. Результаты расчетов для шестиопорной балки моста через р. Казанку в г. Казани 190
4.7. Формирование дополнительных полей напряжений в сечениях балки 196
4.7.1. Постановка и решение задачи 196
4.7.2. Результаты расчетов 199
4.8. Трехмерное напряженное состояние пролетов моста после строительства 205
4.8.1. Постановка задачи 205
4.8.2. Трехмерная конечно-элементная модель 206
4.9. Методика и результаты расчетов 211
4.10. Математическое моделирование напряженного состояния балки при потере натяжения арматурных пучков в локальных зонах и процесса его восстановления при использовании балочной модели 217
4.10.1. Определение напряжений в сечениях шестиопорной балки от неравномерной потери натяжения арматурных пучков по длине балки 218
4.10.2. Результаты расчетов 223
4.10.3. Задача восстановления поля напряжений в сечениях балки после неравномерной потери натяжения арматурных пучков путем приложения дополнительной системы сосредоточенных сил 230
4.10.4. Результаты расчетов 236
4.11. Определение напряженно-деформированного состояния пролетных строений с учетом потери несущей способности и определение уровня дефектов 244
4.11.1. Результаты экспериментальных исследований 244
4.11.2. Численное исследование по определению дефектных пролетов 248
4.11.3. Определение потери натяжения пучков 253
4.12. Анализ напряженно-деформированного состояния пролетных строений после натяжения пучков по проекту НИЦ «Мосты» ЦНИИСа 257
4.12.1. Качественный анализ формирующихся полей напряжений по балочной схеме 257
4.12.2. Анализ НДС пролетов после реконструкции в трехмерной постановке 260
4.12.3. Выводы 262
Заключение 264
Список использованной литературы
- Построение матриц жесткости трехмерных изопараметри-ческих конечных элементов
- Анализ результатов расчетов в линейной постановке задачи
- Вариационная постановка задачи о НДС грунтовых откосов и ее численная реализация
- Последовательность силовых схем, отвечающая конструкции и технологии монтажа пролетных строений
Введение к работе
0.1. Объекты исследований
В процессе рассмотрения вариантов прохождения трассы автомобильных дорог и проектирования мостовых переходов проектировщики сталкиваются с задачами выбора наиболее экономичного решения при пересечении поймы и русел рек. Сложившаяся практика проектирования и расчета искусственных сооружений основана на соблюдении требований строительных норм и правил (СНиП), которые в настоящий момент не всегда учитывают новых возможностей по строительным материалам и новым технологиям, возможности проведения комплексных расчетов конструкций с учетом их взаимодействия и взаимного влияния друг на друга в строительных конструкциях, системах, поведения под нагрузкой, более достоверного определения их НДС с помощью новых возможностей ЭВМ.
Внеклассные мостовые переходы сооружаются на реках, часто имеющих глубину воды до 10 м и более, высокие скорости течения вызывающие большие общие и местные размывы. Несущие слои грунтов располагаются на глубине до 40 м от дна реки, при этом высота опор от уровня воды достигает 30 м и более.
Затраты на устройство массивных опор на глубоководных реках в сложных гидрогеологических условиях занимают до 70 процентов от общих затрат труда и времени сооружения мостов, до 60 процентов общей стоимости мостовых переходов.
Как правило, решение задач по проектированию опор, пролетных строений и других конструктивных элементов мостов основано на простейших схемах, эмпирических формулах с большими коэффициентами за-
паса, которые требуют значительных финансовых, материальных и трудовых затрат при сооружении мостовых переходов.
Создание нового конструкторско-технологического комплекса по сооружению легкой, гибкой опоры мостов с фундаментами глубокого заложения, с минимальными материальными, трудовыми и финансовыми затратами потребовало необходимости комплексного расчета опоры, в том числе с учетом несущей способности металлической оболочки буро-набивных свай.
Многоводные реки Европейской части Российской Федерации, протекающие в южном направлении, имеют крутой правый берег, часто склонный к оползню. Решение задач по стабилизации оползневых участков берега в районе прохождения трасс решаются разгрузкой склона, а также их укреплением путем забуривания, армирования и бетонирования свай или другими методами.
Решения по укреплению берегов, склонных к оползню, принимаются на основании геологических данных, а также данных, полученных за наблюдением режимных скважин.
Методики расчетов на оползневых участках, как правило, сводятся к определению кривых линий скольжения и определению коэффициентов запаса устойчивости грунтовых масс.
При принятии решения по крутым, потенциально склонным к оползню грунтовым массивам, очень важно знать распределение напряжений в грунте при разгрузке и укреплении склонов:
при каких условиях возникает предельное напряженное состояние, в результате которого происходит разрушение скелета грунта и переукладка частиц;
прогноз поведения грунтовых массивов под влиянием внешних и внутренних воздействий, изменений условия равновесия в силу различных природных и техногенных причин.
С целью достоверной оценки несущей способности опоры, определения ее трехмерного напряженно-деформированного состояния с учетом дискретного расположения в бетоне арматуры, работы стальной трубы стоек опоры и взаимодействия с прилегающим к опоре грунтом, потребовалась разработка в трехмерной постановке метода и программного обеспечения, проведения анализа напряженно-деформированного и предельного состояний системы опора-грунт.
Особый класс задач, требующих разработки методик - задачи по определению остаточной несущей способности пролетного строения железобетонного моста, состоящего из отдельных блоков, расчет и схема усиления пролетных строений таких мостов. Во второй половине XX века в СССР было построено более 20 мостов с большими пролетами такого типа, при этом все они потребовали ремонта, а некоторые из них обрушились на стадии строительства или эксплуатации.
Основным недостатком составных пролетных строений мостов является расстройство поперечных стыков между блоками и проникновение влаги в стыки между блоками.
Неравномерность потери предварительного напряжения арматурными пучками в пролетных строениях, интенсивная коррозия арматуры при раскрытии стыков, требует особого подхода к вопросам определения остаточной несущей способности всех балок в каждом пролете. Недостоверные результаты могут привести к невозможности полного восстановления несущей способности пролетного строения за счет низких усилий в напрягаемой арматуре или разрушения бетона блоков в случае больших усилий. Сочетание эксперимента и разработанных методик расчета позволяет успешно справиться с этой задачей.
0.2. Обзор исследований по применению метода конечных элементов для решения трехмерных задач теории упругости и пластичности
Основные теоретические положения теории упругости и пластичности, описание применяемых методов решения, а также решение целого ряда частных задач приведены во многих трудах отечественных и зарубежных исследователей, в частности, в монографиях [6,23,60,115,126, 130,132,149,165,190,215,217,223,228].
В связи с бурным развитием вычислительной техники в решении задач механики деформируемого твердого тела, особенно для расчета конструкций и сооружений сложной структуры и их элементов сложной геометрии, широкое развитие получили численные методы. Среди них особое место занимает метод конечных элементов (МКЭ), благодаря своей универсальности в программной реализации и возможности создания полностью автоматизированного цикла расчета. Он основан на замене исследуемого объекта совокупностью конечного числа дискретных элементов, связанных между собою в узлах. Непосредственный переход к расчетной схеме из соображений механики дает возможность естественно сформулировать граничные условия, произвольно располагать узлы сетки, сгущая ее в местах ожидаемого большого градиента искомых величин, применять метод для исследования областей, состоящих из фрагментов различной физической природы и т.д. Важно отметить естественность механической природы МКЭ.
Ход развития метода отражен в работах зарубежных исследователей Дж. Аргириса, Э. Вилсона, М. Айронса, Р.У. Клафа, O.K. Зенкевича, Дж.Т. Одена и других. Значительный вклад в теорию метода конечных элементов содержится в отечественных работах В.А. Постнова, И.Я. Хархурима, А.С. Сахарова, Л.А. Розина, И.Ф. Образцова и других.
Литература, посвященная теории и реализации метода конечных элементов, весьма обширна. Среди целого ряда монографий следует отметить работы [22,61,107,138,166,174,175,193,195,197,205,206,213]. История метода, его современное состояние и его сравнение с другими широко
используемыми численными методами отражены в обзорах [49,106, 177,302].
Одной из распространенных областей применения метода конечных элементов является решение задач прочности трехмерных тел сложной геометрии. Для расчета пространственных структур используются конечные элементы, созданные на основе различных вариационных принципов трехмерной теории упругости. Можно отметить ряд работ, описывающих построение гибридных конечных элементов [244,256,284], трехмерных конечных элементов, построенных на основе смешанного вариационного принципа [287], а в работе [218] реализован конечно-элементный вариант метода сил.
Однако наиболее популярным и часто используемым на практике является метод конечных элементов в форме метода перемещений. Во многих работах отечественных [20,68,69,76,81,87,95,114,118,139,143,145, 155,211] и зарубежных [12,108,200,248,256,305,307,328] авторов приведены схемы построения матриц жесткости конечных элементов, реализован алгоритм решения, а также решен ряд задач, причем точность решения проверялась на многочисленных тестовых примерах.
Широко используются конечные элементы различных форм, степени и вида аппроксимации перемещений. Приведены схемы построения пирамидальных [10,156], цилиндрических [207,231], тороидальных [84,88], призматических (причем с различной степенью аппроксимации перемещений по высоте призмы) [77,113,134], а также эрмитовых конечных элементов [3,117]. Характеристики перекрывающихся конечных элементов описаны в работе [259].
В большинстве случаев достижение приемлемой точности при расчете пространственных конструкций сложной геометрии получается за счет глобального сгущения конечно-элементной сетки. Сходимость результатов, полученных таким образом, называется h - сходимостью. Такой метод повышения точности решения для трехмерных задач теории
упругости описан в работах [50,134,210,262]. Кроме глобального сгущения конечно-элементной сетки иногда используется метод, когда конечно-элементная сетка остается неизменной, а повышается порядок аппроксимации элементов (сходимость решения при таком подходе называется р -сходимостью), или при неизменном количестве элементов и узлов определяется их оптимальное (по точности) расположение в конструкции (г-методика). Сюда же можно отнести группу работ [99,111,114,128,129, 131], описывающих адаптируемый подход в методе конечных элементов.
В некоторых работах [52,58,127] трехмерная конструкция рассчитывается как двумерная. Это достигается за счет введения некоторых допущений на характер деформирования конструкции при решении части задач определенного вида. Причем некоторые задачи могут сводиться к осесимметричным.
В работе [258] описывается совместное использование осесиммет-ричных и пространственных конечных элементов. Стыковка производится за счет введения множителей Лагранжа.
Для повышения точности расчета в местах локальных нагрузок, а также при резком изменении градиента кривизны поверхности конструкции производится локальное сгущение конечно-элементной сетки. Иногда для этих целей используют специальные конечные элементы с различным числом узлов на гранях (кроме того, эти элементы служат для состыковки густой и редкой конечно-элементной сеток). Такие элементы и методы построения конечно-элементных сеток описываются в работах [145,230, 275,277].
Для аппроксимации лицевых поверхностей конструкций в ряде работ [7,144,251,315] приведена методика построения вырожденных конечных элементов на основе трехмерных.
В некоторых работах [87,99,143,182,199,208] применяются различные модификации методики вывода определяющих соотношений метода,
таких как полуаналитический вариант метода конечных элементов, различные уточненные схемы, в том числе и моментная.
В [282] исследуется возможность повышения точности конечно-элементного расчета за счет введения специальных функций перемещений, совместных внутри области и совместных в среднем на границе. А в [323] для вычисления с высокой точностью производных от перемещений (в том числе и напряжений) используют несовместный конечный элемент и вводят специальные условия для обеспечения непрерывности перемещений.
Контроль ложных мод проводится для 20-ти узловых изопара-метрических конечных элементов в работах [123,286]. Исследуются варианты полного и сокращенного интегрирования матрицы жесткости, а также предложены свои комбинированные схемы интегрирования.
В последнее время появились работы, в которых описывается построение трехмерных моделей, сочетающих конечно-элементный расчет с гранично-элементным [123,261,271,314,331].
Имеется целый ряд работ [24,111,131,142,146,233,245], в которых исследуется сходимость трехмерных конечных элементов, а также производится сравнение для трехмерных задач теории упругости метода конечных элементов, метода граничных интегральных уравнений и метода конечных разностей. Отмечаются достоинства и недостатки, присущие каждому из этих методов.
В настоящее время разработаны и используются на практике пакеты прикладных программ, реализующие метод конечных элементов для решения трехмерных задач прочности для реальных конструкций. Можно отметить целый ряд работ [57,79,99,112,192,208,219,227], описывающих схему функционирования таких прикладных программ.
В работах [290,272,281,282,295,312,317,321] приводятся основные схемы построения семейств переходных элементов, построенных на основе одного конечного элемента, в [318] описывается создание и исследова-
ниє переходных трехмерных элементов, используемых при решении термоупругих задач, в [200,201,202,262] для расчета напряженно-деформированного состояния пространственных конструкций применяется комбинированная модель, где кинематические условия упругого сопряжения с обо-лочечными элементами учитываются при помощи метода штрафа. В работах [1,85,121] необходимые условия сопряжения по границе между трехмерными и оболочечными элементами реализуются путем введения в исходный функционал задачи множителей Лагранжа, параметры которых исключаются из числа варьируемых величин на уровне сборки конструкции.
Имеются работы, описывающие переходные элементы одной структуры, когда к одной грани одного конечного элемента пристыковываются несколько других [158]. В работе [4] описан элемент, состыковывающий трехмерные элементы к элементам толстой оболочки, в [319] описывается переходной элемент трехмерных и осесимметричных структур. Обзор попыток создания переходных элементов приведен в [304].
Можно отметить ряд работ [47,48,81,151,182,183], в которых на основе метода конечных элементов реализованы различные методики решения упругопластических задач. Некоторые итоги и перспективы конечно-элементных исследований физически нелинейных задач приведены в [80].
В работах [1,111,120,121,122] приводятся описания численной методики, основанной на моделировании процессов деформирования и разрушения оболочечных конструкций с учетом взаимного влияния этих эффектов при квазистатических термосиловых нагружениях.
По реализации метода начальных напряжений для упругопластических задач с использованием теории течения на основе метода конечных элементов существует много работ, среди которых можно отметить монографии [50,134,156].
0.3. Модели грунтовых сред
Взаимосвязанные процессы деформирования и фильтрации в насыщенных пористых средах составляют сущность многих явлений в природе и служат основой разнообразных технологических воздействий. Важную роль в развитии теории таких процессов играет математическое моделирование, позволяющее прогнозировать и оптимизировать технологические воздействия, интерпретировать и обрабатывать опытные данные. Как правило, они выполняются на основе модельных представлений теории фильтрационной консолидации, берущей начало с пионерской работы Терцаги К. [322] 1925 года. В ней он впервые ввел понятие эффективных напряжений и решил одномерную задачу уплотнения водонасыщенного пористого грунта в виде слоя конечной толщины.
Эта теория получила дальнейшее развитие в трудах Герсеванова Н.М. [65,66], Флорина В.А. [224,225], Цытовича Н.А. [232], Зарецкого Ю.К. [103,104,105], Био М.А. [40,41], Николаевского В.Н. [161,162,163, 164] и других. Общая математическая модель фильтрационной консолидации на базе вариационно-термодинамического подхода создана Био М.А. [40,41]. Ее глубокий анализ с позиций механики сплошной среды проведен Николаевским В.Н. [162,164]. Костериным А.В. на основе вариационных формулировок задач фильтрационной консолидации исследована их корректность и предложены обоснованные численные методы их решения [93,94]. В последние годы интерес к теории фильтрации усиливается, о чем свидетельствует быстрый рост числа публикаций по этой тематике [116, 148,247,261,287,306,313,320].
Одним из наиболее характерных свойств мягких грунтовых сред является необратимость его объемных деформаций. Чтобы учесть этот фактор в уравнениях состояния среды, то есть в зависимости между шаровыми компонентами тензоров напряжений и деформаций или между давлением р и плотностью среды р, в модели А.Ю. Ишлинского и П.В. Зволинского грунт полагается идеальной несжимаемой жидкостью, плотность которой при давлениях p
дит мгновенная «упаковка» частиц до плотности рг После этого при p>Ps грунт ведет себя как несжимаемая жидкость с плотностью р^. Разгрузка не
изменяет плотности среды.
В модели идеальной пластической сжимаемой среды Х.А.Рахматул-лина, Н.Я. Сагомяна грунт считается идеальной слшмаемой жидкостью с необратимой объемной деформацией є. Нагружение и разгрузка происходят по несовпадающим кривым. В частном случае линия разгрузки молсет быть параллельной оси давления, что означает несжимаемость при разгрузке.
Ляховым Г.М. была предложена модель водонасыщенного грунта как многокомпонентной сжимаемой среды (твердые частицы различной природы + вода + воздух) с баротропным уравнением состояния, где сжимаемость определяется относительным объемным содержанием и сжимаемостью каждого компонента. Нагрузка и разгрузка происходят по одной кривой, как в нелинейно упругой среде.
Деформационная идеальная упругопластическая модель является обобщением упругой и жесткопластической среды с внутренним трением. Реализация модели в деформационной постановке дает единственность получаемых решений, соосность напряжений и деформаций. Модель сочленяет две теории, на которых базируется современная механика грунтов: теорию упругости и теорию предельного состояния. Для описания модели достаточно обычного набора механических характеристик.
В работе [292] предлагается теоретическая модель, описывающая упругое поведение грунтов под нагрузкой, основанная на использовании принципа сохранения энергии. Соотношения модели используют закон Гука, коэффициент Пуассона и модуль Юнга, зависящие от первого инварианта тензора напряжений и второго инварианта девиатора тензора напрялсений. Параметры, необходимые для предложенных соотношений, определяются при экспериментальных исследованиях на приборах трехосного сжатия по траекториям циклического нагрулсения - нагрузка,
разгрузка. Аналогичная модель для описания нелинейно упругого поведения несвязных грунтов предложена в работе [293], но в ней пуассо-новское отношение предполагается постоянным, а модуль Юнга выражается в виде степенной функции.
Континуальная модель пористой среды с упруго деформирующимся матричным материалом предложена Феденко В.И. [222]. Учет возможностей анизотропии в пространственном распределении пор осуществляется путем использования симметричного тензора второго ранга в качестве меры пористости. Данная модель позволяет учесть диссипативные явления, возникающие в пористой среде даже при упругом деформировании матричного материала, а также получить оценку эффективных модулей в зависимости от пористости. В статье [ПО] описана модель изотропной упругой среды, содержащей анизотропное включение, симметрия упругих свойств которого может изменяться. В рамках теории упругости решена задача о перераспределении напряженно-деформированного состояния в среде.
Во всех рассмотренных моделях пренебрегается сдвиговыми напряжениями и деформациями среды. Такой подход дает удовлетворительные результаты в задачах динамики водонасыщенных грунтов, а также для неводонасыщенных сред при больших давлениях.
Для неводонасыщенных грунтов получили развитие модели, учитывающие сопротивление среды сдвиговым напряжениям. При малых нагрузках удовлетворительные результаты дают модели грунта, как линейно или нелинейно упругой среды.
В работе [263] приводится описание упругопластической модели, отражающей поведение неводонасыщенных грунтов. Данная модель позволяет дать качественное описание их следующих особенностей: объемные деформации при изменении степени водонасыщения (как набухание, так и усадка), поведение грунта при первичной консолидации при различном во-донасыщении, предельное состояние грунта в момент разрушения. В гиперболической модели деформационного типа [279] уплотненных грун-
тов не полное водонасыщение учитывается с помощью аппроксимирующей функции для касательного объемного модуля деформации. Остальные численные значения параметров модели берутся из результатов трехосных испытаний и изотропной компрессии.
Повышение уровня нагрузок вызывает появление пластических свойств, поэтому в этих случаях используются более сложные модели.
В модели несвязного грунта, основанной на механизме скольжения, зернистая среда рассматривается как набор частиц, пересекаемый большим числом потенциальных плоскостей скольжения [235]. На каждой такой плоскости, проходящей через контакты частиц, даются выражения для действующих на ней напряжений и связанных с ней составляющих разрывов поля напряжений. Суммарная деформация получается интегрированием по всевозможным ориентациям плоскостей скольжения. Связь напряжений с деформациями выводится в матричном виде сначала для отдельной плоскости, а затем вновь интегрированием дается переход к общей матрице податливости элемента грунта. Помимо угла трения минеральных частиц друг по другу остальные параметры определяются из испытаний на гидростатическое сжатие и трехосное напряжение. Работа [136] на основе системного подхода излагает методологию построения имитационной модели распределения давления в зернистых средах. На этой основе описывается характер распределения давлений и деформаций сжатия в безраспорной зернистой среде и однородном или слоистом грунтовом основаниях.
Более общая модель грунтовой среды была предложена С.С. Григоряном. В ней учтены основные свойства грунтов, существенные при кратковременных волновых процессах - нелинейность и необратимость диаграммы объемного сжатия с участком упругих деформаций при малых давлениях, упругопластический сдвиг, зависимость предела упругости при сдвиге от давления. Сдвиговая деформируемость в допредельном состоянии соответствует линейно упругой среде, а в предельной схеме Прандтля - Рейса с условием пластичности Мизеса - Шлейхера.
Олисовым В.А. при решении задачи о распространении в грунте плоской одномерной волны применена модель упругопластическои среды с линейными диаграммами нагрузки и разгрузки. В модели Зволинского А.В. деформирование грунта происходит упруго при р < ps и пластически
с постоянной плотностью при р > ps и условием Кулона.
В двухповерхностной упругопластическои модели Болдырева Г.Г. [44] параметр упрочнения заменен переменной, описывающей рост и образование пор и микротрещин. Модель включает две поверхности, поверхность нагружения и поверхность текучести. Очертания последней приняты неизменными в процессе деформирования грунта. Поверхность нагружения расширяется или сжимается изотропно от начальной поверхности текучести. Пластическая деформация определяется из ассоциированного закона течения. В работе [125] предложен упрощенный вариант упруго-пластической модели анизотропных материалов, определяемый пятью пластическими параметрами.
Сопоставление экспериментальных и расчетных данных свидетельствует о том, что дважды структурно-неоднородная расчетная модель достаточно адекватно отражает основные особенности деформирования крупноскелетных грунтов в условиях сложного напряженного состояния [62]. В предлагаемой расчетной модели используются традиционные в механике грунтов расчетные характеристики прочности и деформируемости, которые могут быть определены по результатам испытаний грунтов известными способами.
Модели грунтовых сред с учетом дилантансионных эффектов, то есть зависимости объемных деформаций (или плотности среды) от величины пластических сдвиговых деформаций, рассматривались в работах В.Н. Николаевского [160,161]. С позиции трехмерной линеаризованной теории устойчивости деформируемых сред при малых докритических деформациях развита модель упруговязкопластического грунта, учитывающая сжимаемость и дилантансию [213].
Для описания трехмерного нелинейного напряженно-деформированного состояния и дилантансии песка предназначена гипопластическая модель [237], разработанная без обращения к представлениям упруго-пластической теории, таким как поверхности скольжения, пластический потенциал и разделение деформаций на упругую и пластическую части. Модель содержит четыре параметра, определяемые с помощью трехосных компрессионных испытаний.
В работе [243] представлена упругопластическая модель песка при статическом нагружении в условиях трехосного сжатия и растяжения. Принято изотропное поведение материала в широком диапазоне действующих напряжений и анизотропных эффектов. Модель позволяет учесть предельные напряжения, дилантансию и прочностные свойства с учетом действующих главных напряжений и их поворот. В работе А.А.Вовк и Б.В. Замышляева [54] на задачах о распространении взрывных волн в грунтовой среде сделаны выводы о слабом влиянии дилантансии на параметры волн по сравнению с эффектами пластического объемного деформирования.
Для исследования поведения грунтовой среды под действием циклических нагрузок предназначена упрощенная модель уплотнения сыпучей среды [301]. Модель описана при помощи двух определяющих уравнений. Первое, дифференциального типа, описывает уплотнение грунта в зависимости от амплитуды деформации и актуальные состояния уплотнения. Второе уравнение связывает амплитуды деформации и циклические напряжения.
В микромеханической статистической модели для описания процесса уплотнения сухих песков под действием циклических деформаций в отличии от предшествующих моделей, где учитывалась только средняя пористость песков, количественным параметром выступает пространственное распределение пористости [295].
В работе [238] предложена упругопластическая модель состояния песка при циклическом нагружении, в которой принят неассоциированный
закон течения. Функция текучести и пластический потенциал являются обобщенными формами из модифицированной модели кемлей. Параметр изотропного упрочнения определяется пластической работой, связанной с разными частями объемной и девиаторной деформации. Поверхность текучести имеет коническую форму с поперечным сечением в виде шестиугольника со скругленными углами. Независимость условия текучести от направления при трехосных испытаниях на сжатие и растяжение обобщена на общий случай напряженного состояния. Модель способна прогнозировать дренированные и недренированные условия, а также учитывать влияние мембранной проницаемости.
В модели динамического деформирования упругопластической пористой среды применяется замкнутая система уравнений Прандтля-Рейса [25]. Гетерогенный подход позволил описать свойства таких сред в широком диапазоне скоростей нагружения в рамках теории пластического течения с кинематическим упрочнением. Гидродинамическая теория деформирования пористых сред обобщена на случай учета девиаторных составляющих тензора напряжений среды. Неизвестные функции модели определяются из анализа одноосных деформаций соответствующей сферической ячейки.
Математическая модель фильтрационной консолидации насыщенной пористой среды под действием внешних поверхностных сил включает в себя суммарное уравнение движения (квазиравновесия) фаз, условия неразрывности (баланса масс), закон фильтрации, реологическое соотношение для пористого скелета, граничные и начальные условия.
Вариационный подход применялся к исследованию фильтрации в деформируемой пористой среде: изучался процесс консолидации (уплотнение насыщенной пористой средой под действием внешней нагрузки). Вариационная формулировка (принцип) служит при этом основой для построения приближенного решения задачи, например, методом конечных элементов. В работе [296] построен функционал, экстремум которого дости-
гается на решении задачи консолидации нелинейно деформируемой насыщенной пористой среды с учетом конечных деформаций. Этот весьма общий результат получен на основе использования плодотворной идеи -перехода от исходных уравнений и граничных условий задачи к их производным по времени («скоростям»). Относительно «скоростей» задача становится линейной, причем коэффициенты уравнений параметрически зависят от текущих значений самих параметров состояния. В [296] освещено также современное состояние вариационной теории консолидации (дан обзор соответствующих вариационных принципов).
Контактным задачам теории упругости и вязкоупругости посвящена обширная литература [59,83,91]. Прогресс здесь связан, главным образом, с возможностью использования классического аппарата теории функций комплексного переменного для получения аналитического решения соответствующих задач.
Аналогичные по постановке задачи для насыщенных пористых сред, формулируемые в рамках схемы фильтрационной консолидации, менее изучены. Число аналитических решений краевых задач теории консолидации невелико. Основным задачам посвящены работы [119,128,252, 260,263,290]. В [128] получено приближенное решение контактной задачи о давлении штампа на полуплоскость, в [67,133] построено решение задачи о консолидации в тонком слое и в полосе.
Одним из наиболее простых и эффективных численных методов является метод конечных разностей (МКР), который применяется при решении как линейных, так и нелинейных дифференциальных уравнений с различными условиями на границе области [14,199]. К достоинству его относится простой вид аппроксимации дифференциальных уравнений, который удобен для программирования. Но метод конечных разностей является эффективным для областей достаточно простой формы.
В последнее время при решении задач фильтрации широкое применения находит более универсальный метод - метод конечных элементов
(МКЭ) [61,90,107,153,206,213]. Его преимущество - высокая устойчивость в применении к областям сложной геометрии и неоднородной структуры. МКЭ основан на замене исследуемого объекта совокупностью конечного числа дискретных элементов, связанных между собой в узлах. Непосредственный переход к расчетной схеме из соображений механики дает возможность естественно формулировать граничные условия, произвольно располагать узлы сетки, сгущая ее в местах ожидаемого большого градиента искомых величин, применять метод для исследования областей, состоящих из фрагментов различной физической природы и т.д.
Разработанные в последнее время новые вычислительные схемы, реализующие МКЭ и МКР значительно расширили класс задач, решаемых численными методами.
Использование механики грунтов в инженерной практике с каждым годом становится все более широким. Так, на основе решения ряда конкретных задач механики грунтов, а также проверки результатов в экспериментах оказалось возможным разработать весьма прогрессивный, дающий значительную экономию средств метод проектирования фундаментов по предельным состояниям грунтовых оснований [94].
В [289] для расчета грунтовых оснований и фундаментов используется метод конечных элементов, для которого записаны определяющие соотношения упруговязкопластического поведения грунта в связанной постановке, когда учитывается фильтрация жидкости в грунте (закон Дарси). Из принципа возможных приращений в скоростях записана соответствующая конечно-элементная формулировка задачи. Описана процедура пересчета входных параметров задачи, определяемых при стандартных испытаниях поведения грунта под нагрузкой, в материальные константы и функции, присутствующие в используемых определяющих соотношениях.
В работе [270] исследуется проблема уплотнения насыщенных пористых сред. Особенностью предлагаемого подхода является более детальный учет зависимости проницаемости породы от напряжений в скелете и
давления флюида, который в общем случае может быть сжимаемым (например, газ). Проницаемость принимается нелинейно зависящей от перечисленных величин. Такой учет считается существенным на больших глубинах (в геодинамических задачах) или при значительных нагрузках на поверхности. Разработана конечно-элементная модель деформации насыщенной породы. Рассматриваемая нелинейная система уравнений использована в проблеме динамики грунта при нагружении на поверхности полупространства.
Орехов В.В. в работе [176] приводит описание комплекса вычислительных программ, предназначенного для решения задач взаимодействия фундаментов с грунтовыми основаниями при статических и динамических воздействиях на основе метода конечных элементов.
В работе [220] Фадеева А.Б., Матвеенко Г.А., решение трехмерной задачи сведено к решению ряда осесимметричных задач разложением узловых нагрузок и перемещений по окружной координате в ряды Фурье. Грунт рассматривается как идеально упругопластическая среда с поверхностью текучести, описываемой критерием Боткина в октаэдрических напряжениях.
Миховой Л. [153] на основании пространственной теории Био М.А. консолидации грунта решена осесимметричная задача с применением метода конечных элементов. Грунт принят как двухфазная система, состоящая из твердой фазы (скелета) и жидкой фазы (жидкости в порах скелета). Принято, что скелет линейно деформируемый материал. Жидкость неде-формируема при полной водонасыщенности грунта и деформируема при наличии газа.
В работе [219] Фадеева А. Б., Репиной П. И., Глыбина Л.А программа обеспечивает получение серии упругопластических решений для заданной последовательности нагружения гравитационными силами, пошагового приложения строительных нагрузок, постадийной выемки котлованов для подземных выработок, введения на любом этапе конструктивных эле-
ментов (фундаментов и т.п.) Отличительной особенностью является возможность введения на любом этапе заданных перемещений узлам. Модель среды - билинейная, упругопластическая, с критерием текучести Кулона.
Stematiu D., Paunescu D. [315] предлагают модель поведения грунта с неполным насыщением под действием внешней нагрузки. Модель учитывает взаимодействие между тремя составными фазами грунта: твердым скелетом, водой и воздухом. Система дифференциальных уравнений неразрывности и баланса решена численно методом конечных элементов.
В статье Стояновича Г.М. [212] расчеты выполнялись по разработанной автором аналитически-экспериментальной методике учета вибродинамического воздействия и упругопластического напряженно-деформированного состояния земляного полотна на основе метода конечных элементов и эмпирических зависимостей снижения прочностных свойств грунтов в условиях Кулона.
В [310] выполнено численное моделирование локализации неупругой деформации в насыщенных песчаных образцах в условиях динамического нагружения в отсутствии дренажа. Использован метод конечных элементов для совместного решения уравнения баланса масс и уравнения движения.
Ng A.K.L., Small J.C. [295] методом конечных элементов исследовали консолидированное поведение ненасыщенных грунтов.
В работе [273] алгоритм адаптивного улучшения сетки разработан для нелинейных расчетов в геомеханике и основан на сглаженном представлении диаграммы напряжений в методе конечных элементов. Использована оценка ошибки в относящемся к приращениям инварианте деформаций сдвига для преобразования сетки в процессе нагружения. Алгоритм разработан в результате анализа задачи пассивного давления грунта с использованием идеальной упругопластической модели Кулона-Мора. Использован смешанный гидромеханический анализ поведения грунта в
процессе дренирования. Во всех случаях преобразование сетки признается успешным в областях с высоким градиентом деформаций.
В работе [324] представлена трехмерная численная модель, деформации которой описываются согласно нелинейной теории упругости. Математическая формулировка связанных задач представлена четырьмя уравнениями на основании принципа сохранения массы и энергии, а также уравнением равновесия. Для описания движения жидкости и воздуха в пористой среде используется закон Дарси А. В модели используются трехмерные изопараметрические двадцати узловые элементы. Метод позволяет моделировать естественно нелинейные параметры грунта.
Власюк А.П., Мартинюк П.М. [53] исследовали численное решение двумерной задачи фильтрационной консолидации глинистых грунтов. Решение получено методом конечных элементов.
В работе [222] основным недостатком приемов, рассматривающих условия предельного равновесия на некоторых кинематически возможных поверхностях скольжения - обычно круглоцилиндрических, является упрощенная картина напряженного состояния. Обычно предполагается, что в грунте действуют только вертикальные напряжения, пропорциональные глубине рассматриваемого участка поверхности скольжения от дневной поверхности. Кроме того, для определения наиболее опасного сочетания сдвигающих и удерживающих сил необходимо проведение множества Расчета по многим возможным поверхностям скольжения; оползневые тела при этом подразбиваются на достаточно крупные блоки, что вносит в результаты расчетов дополнительные погрешности. Достаточно эффективным является сочетание методов конечных элементов и предельного равновесия.
В работе Бережного Д.В., Голованова А.И., Паймушина В.Н., Сидорова И. Н. [29] разрабатывается конечно-элементная методика расчета водонасыщенной пористой среды, взаимодействующей с деформируемыми конструкциями.
Пшеничкиным А.П. [194] рассматривается деформирование во времени двухфазного грунта, который включает в себя два процесса, протекающих одновременно. Это - процесс формоизменения и объемного изменения во времени скелета грунта, происходящий в результате деформирования вязких связей между частицами грунта. Принимается, что сначала происходит выдавливание из пор воды (первичная консолидация), а затем деформирование во времени идет за счет ползучести скелета грунта (вторичная консолидация). По методу эквивалентного слоя грунта Цытовича Н.А. по теории фильтрационной консолидации, получено решение задачи уплотнения грунтов водонасыщенного основания.
В работе [43] Бойко И.П., Сахарова В.А. приведены результаты решения двумерных и трехмерных линейных и нелинейных задач взаимодействия фундаментов соседних зданий с применением численных методов на базе системы «VESNA». Используется теория пластического течения, не-ассоциированный закон деформирования грунтов основания и модифицированный критерий Мизеса-Губера-Боткина, учитывается конструктивная нелинейность системы «основание-фундамент-конструкции». Дано сравнение результатов решения задач моделей с коэффициентом жесткости основания и модели нелинейно-деформируемого слоистого грунтового массива.
В работе [233] получено точное решение пространственной задачи теории фильтрационной консолидации при осевой симметрии, которое отличается от известных приближенных решений учетом в расчетных формулах коэффициента Пуассона грунтового скелета. Это позволяет более достоверно прогнозировать развитие во времени деформаций и напряжений водонасыщенных оснований.
Основными целями работы являются постановка и решение следующих задач:
разработка высокоточной конечно-элементной методики определения напряженно - деформированного и предельного состояний железобетонных массивных тел сложной геометрии с учетом произвольности ориентации и мест расположения в теле бетона армирующих стержней, работающих на растяжение - сжатие; определения напряженно - деформированного и предельного состояний элементов конструкций указанного выше класса при наружном или произвольном внутреннем подкреплении тела бетона тонкой мембраной;
создание уточненных математических моделей, методов и программного обеспечения для определения напряженно-деформированного и предельного состояния:
железобетонных опор внеклассных мостов с учетом и без учета взаимодействия с окружающим грунтом;
крутых склонов, состоящих из многослойных грунтов, потенциально склонных к оползню.
разработка уточненных методов решения прямых и обратных задач по определению остаточной несущей способности и ее восстановления в предварительно напряженных балках пролетных строений мостов из сборных железобетонных блоков;
решение с помощью разработанных методик ряда сложных практически важных задач.
Научная новизна диссертации состоит:
1) в дальнейшем развитии метода конечно-элементного анализа
напряженно-деформированного и предельного состояний массивных трех
мерных тел сложной геометрии и неоднородной структуры;
2) в создании научных основ и программного обеспечения для мате
матического моделирования механического поведения и разрушения слож
ных строительных конструкций и сооружений на основе использования
современных достижений в области механики деформируемого твердого
тела, механики грунтов, вычислительной математики и информационных технологий;
в разработке методики и проведении исследования напряженно-деформированного и предельного состояний грунтовых массивов для оценки опасности оползневых явлений крутых склонов и способов их подкрепления и разгрузки;
в создании научных основ, программных средств и методики расчетно-экспериментальных исследований для определения остаточной несущей способности и уровня дефектов пролетных строений мостов из сборных железобетонных блоков, использование которых позволяет разработать наиболее рациональные и надежные проекты реконструкции и восстановления несущей способности вышедших из строя мостов указанного класса.
в разработке уточненных методов решения прямых и обратных задач по определению остаточной несущей способности и ее восстановления в предварительно напряженных балках пролетных строений мостов из сборных железобетонных блоков.
Достоверность основных научных результатов обеспечивается применением строгих математических методов для построения основных соотношений, сравнением полученных результатов решения некоторых тестовых задач с результатами их решения другими авторами и приближенных постановках на основе более простых моделей, анализом сходимости решений рассмотренных задач, полученных на разных конечно-элементных сетках.
Личный вклад автора. Все основные научные результаты, изложенные в диссертации принадлежат автору.
Практическая ценность результатов диссертации:
1) проведенный численный анализ напряженно-деформированного состояния опор, в том числе при их взаимодействии с грунтом, исследованных мостовых переходов позволил выявить недостатки их конструк-
ций, разработать новые конструктивные варианты, обладающие необходимыми жесткостными и прочностными характеристиками;
исследование напряженно-деформированного и предельного состояний грунтовых массивов позволило создать методику оценки опасности оползня крутых склонов рек и способов их разгрузки и укрепления;
разработанные научные основы, программные средства и методики расчетно-экспериментальных исследований для определения остаточной несущей способности и уровня дефектов пролетных строений мостов из сборных железобетонных блоков позволяют принять эффективное решение по их восстановлению и реконструкции;
разработанные методы нашли применение в практике проектирования элементов конструкций рассматриваемого класса, что осуществлено заинтересованными организациями и подтверждено соответствующими актами внедрения.
Результаты работы докладывались:
на III, VI, XI, XII международных симпозиумах «Динамические и технологические проблемы механики конструкций и сплошных сред» (г. Москва, 1997 г., 2000 г., 2005 г., 2006 г.)
на I международной конференции «Модели механики сплошной среды, вычислительные технологии и автоматизированное проектирование в авиа- и машиностроении», г. Казань, 1997 г.;
на международном конгрессе МКПК-98 «Пространственные конструкции в новом строительстве и при реконструкции зданий и сооружений (теория, исследования, проектирование, возведение)», г. Москва, 1998 г.;
на международной конференции «Численные и аналитические методы расчета конструкций», г. Самара, 1998 г.;
на республиканской научно-практической конференции «Интеллектуальные системы и информационные технологии», г. Казань, 2001 г.;
- на республиканских научных конференциях, проводимых в КГ АСУ
(г. Казань, 2004 г., 2005 г., 2006 г.)
Основные результаты диссертации изложены в монографии и работах, в которых соавторы принимали участие в постановке задач, создании программного обеспечения, проведении расчетов реальных конструкций, обсуждении полученных результатов. По материалам диссертации опубликовано 25 работ, в том числе одна монография, из них 4 единоличных публикаций.
Структура и объем работы.
Диссертация состоит из введения, четырех глав, заключения, списка литературы и изложена на 304 страницах, содержит 2 таблицы, 185 рисунков. Список литературы - 331 наименование.
Во введении проведен обзор исследований по применению метода конечных элементов для решения трехмерных задач по теории упругости и пластичности, отмечены работы, в которых излагаются существующие модели грунтовых сред и решения задач механики грунтов.
В первой главе выводятся основные соотношения метода конечных элементов для трехмерных задач теории упругости. Описывается построение матриц жесткости трехмерных конечных элементов в виде параллелепипеда и треугольной призмы с криволинейными гранями, рассматриваются методические вопросы, связанные с заданием внешних нагрузок, вычислением напряжений, а также строятся новые высокоточные трехмерные конечные элементы с учетом дискретного расположения в теле конструкции армирующих стержней с произвольной ориентацией и наличия внутри или на наружной поверхности тела подкрепляющих мембран. Излагается методика и алгоритм решения задач по исследованию предельного состояния железобетонных конструкций. Дается описание созданного прикладного программного обеспечения. Приведены решения ряда тестовых задач.
Во второй главе исследовано напряженно-деформированное состояние в стойке 11-ой опоры моста через р. Каму в составе мостового перехода через реку Кама у с. Сорочьи Горы, а также изучено влияние работы
/
металлических труб буронабивных свай на распределение полей напряжений в ростверке и буронабивных сваях в местах их соединений, где опора представляет собой сложную по геометрическим очертаниям и структуре конструкцию, относящуюся к классу комбинированных механических систем.
Третья глава посвящена моделированию напряженно-деформированного и предельного состояний многослойных грунтов, взаимодействующих с деформируемыми элементами строительных сооружений. В этой же главе даны основные сведения о физико-механических характеристиках грунтов, их модели, описывающие упругие и пластические деформации, а также соотношения теории прочности, проведено детальное исследование напряженно-деформированного и предельного состояний крутого склона у мостового перехода через р. Каму у с. Сорочьи Горы.
Проведение таких исследований было продиктовано тем, что разработанные к настоящему времени и используемые в практике проектирования сооружений методы расчета устойчивости грунтовых откосов и величины оползневого давления на проектируемую подпорную стенку основываются на априорном задании линии скольжения, ограничивающей область с достигшими предельного значения напряжениями. Использование таких методов не позволяет найти распределение полей напряжений в слоях грунтового откоса, и, как следствие, выделить места концентрации напряжений. В связи с этим для расчетного определения координат размещения удерживающих возможный оползень свай-шпонок требуется предварительное вычисление полей напряжений во всех слоях откоса.
Результаты такого определения полей напряжений во всех слоях откоса, полученные на основе разработанных методов и идеально упруго-пластической модели мягких слоев грунтов, приведены в третьей главе. Она содержит также результаты расчетов откоса до и после
его разгрузки за счет террасирования склона при действии нормативных нагрузок, а также аналогичные результаты, отвечающие действию перегрузок. Сформулирован весьма практически важный вывод о том, что крутой склон у рассматриваемого моста является устойчивым и имеет значительный запас устойчивости в отношении образования оползневых разрушений.
В этой же главе разработанные методы использованы для исследования напряженно-деформированного и предельного состояний опор моста через р. Архаровку, являющегося составной частью мостового перехода через р. Каму у с. Сорочьи Горы. Отличием представленных здесь результатов от результатов, рассмотренных ранее, является детальный анализ механизма взаимодействия опоры моста с окружающим грунтом.
Геометрическая модель рассчитываемой конструкции включает в себя армированную бетонную опору моста в стальной трубе, а также некоторый объем грунта, прилегающий к буронабивным сваям.
В четвертой главе даны постановка и решение задач механики статического деформирования трансформирующихся систем на примере пролетных строений мостов из сборного железобетона. Исследования основаны на использовании балочной расчетной схемы для пролетных строений. Изложены особенности конструкций пролетных строений из сборных железобетонных блоков и их монтажа путем создания предварительных напряжений сжатия. Установлена последовательность силовых и расчетных схем, отвечающих конструкции и технологии монтажа пролетных строений. Сформулированы и решены задачи по управлению полем напряжений в поперечных сечениях балки путем натяжения дополнительных арматурных пучков.
Аналогичное математическое моделирование пролетных строений мостов из сборного железобетона проведено методом конечных элементов по трехмерной модели, в которой учитываются все особенности геометрии
конструкции. При этом практически точно воспроизводится подкрепляя-ющее действие арматуры и предварительно натянутых пучков по месту их расположения.
В рамках использования балочной модели проведено математическое моделирование напряженного состояния балки моста при потере натяжения арматурных пучков в локальных зонах, а также сформулирована и решена задача восстановления поля напряжений в сечениях балки путем приложения дополнительной системы сосредоточенных сил.
В четвертой главе приведены постановка и решения прямых и обратных задач по определению в трехмерной постановке напряженно-деформированного состояния пролетных строений с учетом потери несущей способности и уровня дефектов. При этом задача оценки поврежденности конструкции ставится как определение величины потери предварительного натяжения пучков в зонах, где возникают максимальные растягивающие напряжения от дополнительной нагрузки, и наблюдается раскрытие блоков при эксплуатации моста.
Первым этапом в исследовании остаточной несущей способности (текущего напряженного состояния) пролетов моста было проведение экспериментального исследования жесткости балок на изгиб при действии дополнительной поперечной нагрузки. В качестве этой нагрузки использовался грузовой автомобиль, последовательно размещенный в середине всех пролетов. В результате геодезических измерений определялись величины прогибов в базовых точках, расположенных вдоль осевой линии нижней плиты балок. Точность этих измерений была на границе значений прогибов в ненагруженных пролетах, что привело к осцилляции данных. Если строить упругую линию балок моста прямо по этим результатам, то соответствующие эпюры получаются нереальными. Поэтому проводилось сглаживание экспериментальных кривых, в результате которых строились сплайновые кривые необходимой гладкости.
Далее расчетным путем подбирались величины потерь натяжения пучков, чтобы упругие линии при всех нагружениях максимально совпадали с экспериментальными. Раскрытие блоков моделируется путем «размазывания» этих дефектов по части объема (по тем зонам, в которых наблюдается реальное трещинообразование и раскрытие блоков) и использования модели идеально-пластического тела с условием на продольные растягивающие напряжения. Упругопластическая задача численно решается с применением итерационной процедуры типа метода начальных напряжений.
Основные результаты, выносимые на защиту:
разработка высокоточной конечно-элементной методики определения напряженно - деформированного и предельного состояний железобетонных массивных тел сложной геометрии с учетом произвольности ориентации и мест расположения в теле бетона армирующих стержней, работающих на растяжение - сжатие;
разработка высокоточной конечно-элементной методики определения напряженно - деформированного и предельного состояний элементов конструкций указанного выше класса при наружном или произвольном внутреннем подкреплении тела бетона тонкой мембраной;
создание уточненных математических моделей, методов, программного обеспечения для определения напряженно-деформированного состояния опор внеклассных мостов, больших массивов грунта, потенциально склонных к оползню, строительных конструкций и сооружений при их статическом нагружении;
постановка и решение прямых и обратных задач по определению остаточной несущей способности и ее восстановления предварительно напряженных пролетных строений мостов из сборных железобетонных блоков;
решение с помощью разработанных методик сложных, практически важных задач, на основе использования трехмерных уравнений теории
упругости и пластичности, прикладных вариантов теории тонкостенных
конструкции, современных достижении механики грунтов.
Построение матриц жесткости трехмерных изопараметри-ческих конечных элементов
Основным расчетным элементом, используемым в работе, является 20-ти узловой квадратичный конечный элемент, представляющий собой искривленный параллелепипед, приведенный на рис. 1.1. Для его построения будем пользоваться соотношениями трехмерной теории упругости в сочетании с изопараметрическим подходом в определении перемещений и геометрии.
Для моделирования зон перехода геометрии будем использовать конечный элемент искривленной треугольной призмы, изображенной на рис. 1.2 Для аппроксимации перемещений и геометрии в плоскости основания призмы используются L - координаты. 20-ти узловой конечный элемент
Для аппроксимации геометрии и перемещений изопараметрического пространственного элемента используются безразмерные координаты ЛХ , в которых исходный искривленный параллелепипед преобразовывается в единичный куб с гранями - 1 ,л,С 1- Это достигается за счет определения функций связи исходных координат x,y,z с , г,С,, которые представляются в виде X N У = Уп \ Z и-1 .V члД,лД (1.17) где xnyn,zn - декартовы координаты п - го узла элемента, \j/w( , л,С) функции формы элемента, которые имеют вполне определенный вид для каждого из элементов с рис. 1.1, 1.2.
Для 20-ти узлового элемента для последовательности узлов, приведенной на рис 1.1, эти функции имеют вид: -Й?(1-т12)(1 + )(1 + )- 2л;2(1Ч2)(1ЦД)(1 + л(41 = 1Д (1.18) где 4/ Л/ ЛІ - локальные координаты узловых точек. При такой аппроксимации х, у, z обладают свойством непрерывности при переходе межэлементных границ.
В соответствии с изопараметрической концепцией вектора переме щений у = и і +v +wfc определяется в виде, аналогичном (1.17), то есть (1
В качестве безразмерных координат выбираются координаты С(-1 С, і) по высоте призмы и L - координаты (0 Lx, L2, L3 і) по поперечному сечению призмы (рис 1.2). L - координатой точки треугольника является отношение величины перпендикуляра, опущенного на одну из сторон треугольника, к высоте треугольника, опущенной на эту же сторону. Компоненты радиус - вектора геометрии и вектора перемещений связаны с их узловыми значениями посредством функций формы (p„(Z,1}Z,2 - 3 C) Для треугольной призмы
В случае задания сосредоточенных сил, приложенных в узлах, достаточно определить соответствующие проекции вектора нагрузки на декартовые оси.
Если задаются распределенные по поверхности нагрузки или массовые силы, то необходимо привести их к узловым силам. Для грани 20-ти узлового элемента вводится система поверхностных координат ос,3, в которых изображенная грань будет квадратом с —1 ос,Р 1 (рис. 1.3).
Тогда поверхность этой грани в трехмерном Евклидовом пространстве с ортами і ,j, к будет описываться векторным уравнением ?(a,$=J:rHN„(a,p), (1.37) где гп- радиус-векторы узлов грани, Nn(a,fi)- функции формы для грани. Для последовательности узлов, приведенной на рис. 1.3, эти функции имеют вид / = 3ggJ-25?Ti?(l + w)(1 + i](ti)_j Работа внешних сил для нормального давления р будет определяться в виде интеграла A = \\p-vndS, (1.39) s где vn - проекция перемещений v на нормаль к поверхности грани Я, которая в каждой точке определяется как векторное произведение дг дг да ар дг дг п = — х — да ЭР у (1.40) Элемент площади поверхности dS в общем случае определится следующим образом dS = дг дг да ЭР dad (1.41) После подстановки (1.40), (1.41) в (1.39) выходит Я/ W дг дгЛ s чЭа ар/ (1.42) Аппроксимации для радиус-вектора г уже определены в (1.37). Для вектора перемещений, вследствие изопараметрического подхода, задается аналогичная аппроксимация, то есть v(a,p)=Zv„ATw(a,p). (1.43) I Распределение давления по поверхности грани удобно представить аналогично (1.37) и (1.43) Xa,p)=X#A(a,p), (1.44) 8 и-1 и задавать его можно в виде узловых значений.
После подстановки (1.37), (1.43), в (1.44) получается 8- 8 88 11 ЯЛ/" f)N Л А=1У\ ZpjT, 1гмхгиПАГ,.(а,р) .(а,р) а# . (1.45) у-1 Vy-i m-in-i -1-і да dp J Если обозначить Л . = П Да,Р)ЛгДа,р)- .- р, (1.46) то для соответствующих узловых сил в і-ом узле (нумерация локальная) будут справедливы выражения 7=1 т,п=\ Py=IlPj1L (Zn,Xr, XmZn К, , (1 -47) j=\ m,n-l 8 =ЕЛ Т,(ХтУп ymZn)lvmn 7=1 яі,п-1 Интегралы (1.46) не зависят от геометрии конкретного элемента, что позволяет отделить задачу их вычисления от задачи вычисления соответствующих узловых сил по формулам (1.47).
Анализ результатов расчетов в линейной постановке задачи
В качестве реального примера применения разработанных методов и программного обеспечения, описанных в предыдущих главах, в данном пункте рассматривается задача определения полей напряжений в стойке 11-ой опоры мостового перехода через реку Кама, построенного у села Сорочьи Горы в Республике Татарстан. Соотношения геометрических размеров стойки (размеров ее поперечного сечения и длины) позволяет отнести ее к классу стержней, работающих, главным образом, на сжатие и косой изгиб. Проектные расчеты для нее были выполнены в ОАО «Гипро-трансмост» (г. Москва), исходя из которых, выбраны ее геометрические размеры, конфигурация поперечных сечений по высоте, а также продольное и поперечное армирование. Схема армирования стойки, заданная проектной организацией в виде конструкторской документации, достаточно громоздкая и сложная. Поэтому она в данной работе не приводится.
С учетом конфигурации стойки для ее расчетов использовались изо-параметрические конечные элементы в виде параллелепипедов, описанные в первой главе. Учет дискретного расположения арматуры в теле бетона производился на основе методики, описанной в первой главе. Разбивка тела стойки на конечные элементы показана на рис.2.1, которая привела к алгебраической задаче с числом неизвестных 7500, оказавшейся достаточной для определения полей напряжений с необходимой для практики точностью.
На рис.2.1 в) показаны действующие на стойку опоры усилия и моменты, приведенные к верхнему сечению стойки. Из них направление Rz соответствует направлению надвижки пролетного строения на опоры. По величине указанные силы и моменты, заданные проектной организацией, принимались следующими Rz = 33,67", = \026,3 Т,МХ = 1067М.
При проведении расчетов эти силы «размазывались» по верхнему сечению и задавались в виде нормальных и касательных напряжений, направленных соответствующим образом. Кроме того, учитывался собственный вес стойки путем задания массовых сил.
Граничные условия задавались как жесткое защемление нижнего сечения, что и принимается на практике при проведении расчетов в соответствии со СНиПами. Следует по этому поводу указать, что при таком моделировании стойки опоры, не учитывающей податливость ростверка опоры, результаты расчетов в окрестности нижнего сечения должны оказаться завышенными по сравнению с реальным полем напряжений.
Результаты расчетов в линейном варианте представлены на рис.2.2-2.7. В частности, на рис.2.2-2.7 даны картины распределения напряжений в сечении 0-0 (плоскость Х=0) (все сечения показаны на рис.2.1 а, б, в). На рис.2.3 дана картина напряжений в сечении 1-1 (плоскость Y=4M). На рис.2.4, 2.5 даны напряжения в сечениях 2-2 и 3-3 (плоскости Y=8M И Y=12M). На рис.2.6 дано распределение напряжений в сечении 4-4 (плоскость Y=23,2M).
Вариационная постановка задачи о НДС грунтовых откосов и ее численная реализация
Исследуемый откос представляет собой протяженный вдоль берега реки склон с достаточно плавным изменением своего профиля. Исходя из этого, представляется возможным и целесообразным описать его НДС в рамках плоской задачи двумерной теории упругости. В этом случае вводятся в рассмотрение следующие функции, объединенные в векторы: і и} =\ U \ " перемещения вдоль координат х ,у; м= і Є_у - деформации, линейные и сдвиговая; ху) {о-Ь Gy Тку) - напряжения, нормальные и касательное. В общем случае напряжения GZ = 0 и определяются после расчетов в виде Gz=ll(Gx +Gy). Разрешающие уравнения берутся в вариационной форме из принципа возможных перемещений в виде #{а}Г {6Б}«4 = ШГ {bu}dA +\\{Р}Т {8u}dS, (3-15) A A S5 где [gj - вектор массовых сил, {Р} - вектор контурных сил, А - площадь расчетной области, sG - граница с заданными усилиями. Предполагается, что кинематические условия выполняются априорно на этапе задания аппроксимаций для перемещений.
Сформулированные в предыдущем пункте условия предельного состояния формально можно записать как условия для вычисления либо касательных напряжений на площадках скольжения для глинистых и песчаных грунтов т genu— т. т т = Т ,ЄСЛи—Т Т (O.10J пр — лр где т н - напряжения, вычисленные по соотношениям линейной теории упругости; либо главных растягивающих напряжений а і для скальных пород стХн,еслис7Хн апр тпр,еслисти а Прт І.Н Пр 1-і. , (3-17) где 5\н - главное напряжение по линейной теории упругости.
Указанная постановка аналогична расчету пластических материалов в рамках идеальной пластичности с отличием в направлении площадок скольжения (соотношение (3.16)) и хрупких материалов с растрескиванием поперек наибольших растягивающих напряжений (соотношение (3.17)).
При достижении предельного состояния в некоторых зонах расчетной области А применение выражений (3.16), (3.17) нарушает уравнение (3.15), т.е. получаемая система напряжений не является уравновешенной. Чтобы избавиться от этого дефекта, используется итерационная процедура типа метода начальных напряжений, широко применяемая при решении упруго-пластических задач механики деформируемых тел. Уравнение для определения {к +1) -состояния при известном А:-том получается из (3.15) добавлением в него неуравновешенной работы от разницы предельных и уравновешенных напряжений, совершаемой на соответствующих деформациях. Структурно его можно записать в виде: Я { }Т {Ы}М = /J {Q}T HdA + tt{P)T MdS + A A SS (3.18) +X tf(l - Tnp/xLk))TLk)6YdA + fl(l - anp/crS) 258idA .AT где у - деформация сдвига на площадке сдвига в направлении напряжений н\ Si - деформация растяжения на главной площадке, где действуют о\ \ Ах - площадь зоны предельного состояния скольжения, Аа - площадь зоны предельного состояния от растрескивания, А, - параметр, регулирующий скорость сходимости.
При достижении сходимости полученные компоненты вектора напряжений {а} требуют коррекции, так как они являются суммой фактических напряжений а [ и «начальных» напряжений -Сн=(і-т,ф/тн)тн в зоне Ат, стГ„ = (і- аир /ст1н)ст1н в зоне Аа. (3.19)
Для этого используются соотношения, обратные тем, которые применялись ДЛЯ ВЫЧИСЛеНИЯ (Tl,CJ2 ПО СГХ , О у, 1ху И Т, 7W ПО 3\, 52 Условие предельного состояния для а з практически не работает, так как всестороннего растяжения в решаемых задачах не встречалось, а на сжатие скальные грунты имеют значительный запас прочности по действующим напряжениям.
Последовательность силовых схем, отвечающая конструкции и технологии монтажа пролетных строений
Седьмой этап - соединение двух балок пролетного строения путем омоноличивания стыка между ними (соединительный элемент 3 на рис.4.1). Силовая схема после окончания этапа - шестиопорная неразрезная балка с поперечным сечением, изображенным на рис.4.1, в котором за-инъектированы все продольные каналы с размещенными и натянутыми арматурными пучками.
Исходя из установленной последовательности этапов строительства и монтажа пролетных строений следует сформулировать следующие основные выводы.
1. Представление пролетного строения моста в виде неразрезной бал ки, лежащей на шести шарнирных опорах и загруженной одновременно всеми приложенными к нему внешними силами, не соответствует реаль ному процессу накопления полей напряжений, деформаций и переме щений.
2. Для описания полей напряжений, деформаций и перемещений в пролетных строениях моста необходимо построить последовательность ра счетных схем, соответствующей этапности строительства и монтажа коробок.
3. Для обеспечения практической степени точности описания процесса накопления полей напряжений, деформаций и перемещений достаточно сформулировать в конце каждого этапа строительства и монтажа пролетных строений статическую задачу о равновесии сформировавшейся к концу этапа конструкции.
4. Окончательное НДС в пролетных строениях представляет собой алгебраическую сумму НДС.
5. Во всех элементах балок, появившихся на n-ом этапе строительства, НДС формируется только начиная с n-ого этапа с добавлением к нему НДС последующих этапов.
Чтобы оценить погрешность представления пролетных строений в виде многоопорной неразрезной балки без учета этапности формирования полей НДС, рассмотрим модельную задачу об НДС балки с бесконечным числом опор, находящейся в равновесии под действием собственного веса. В силу периодичности структуры конструкции достаточно рассмотреть задачу об НДС балки длиной L, у которой в точках опирання на прогиб w наложены условия:
Сопоставляя формулы (4.9) и (4.10), можно увидеть, что использование расчетной схемы в виде многоопорной неразрезной балки, принятой в расчетной практике, приводит к погрешности в 50% при определении нормальных напряжений и в 300% при определении прогибов от действия собственного веса. Отсюда следует вывод о том, что при проектировании и реконструкции пролетных строений мостовых переходов из сборного железобетона расчеты на прочность должны проводиться с использованием расчетных схем, базирующихся на учете этапности монтажа коробок пролетных строений.
При составлении последовательности расчетных схем, основанных на использовании балочной модели, в дальнейшем введем следующие предположения, упрощающие формулируемые в дальнейшем задачи:
1) в силу малости зазоров между консолями по сравнению с длиной пролета (A«L) с целью проведения качественного анализа механики деформирования рассматриваемых конструкций полагаем А « 0;
2) площади поперечных сечений пролетных строений на всех этапах их монтажа считаем постоянными (неизменяемость силовой схемы при переходе от одного этапа к последующему);
3) с целью проведения качественного анализа процесса накопления полей НДС усилия натяжения арматурных пучков, приложенных в дискретных точках Х{ по длине пролета, заменяем распределенными усилиями qHx,qex.
4) изменением положения нейтральной оси ох и расстояний от нее до арматурных тросов при переходе от этапа к этапу пренебрегаем, полагая h « h h « h .
При сформулированных упрощениях приходим к следующей упрощенной последовательности расчетных схем для определения НДС пролетных строений после завершения монтажа одной балки:
1. После завершения первого и второго этапов: четыре балки длиной L и с изгибной жесткостью EI, симметрично лежащие на опорах и нагру женные равномерно распределенной погонной поперечной нагрузкой интенсивностью q2 (действие собственного веса) и погонной осевой силой qex, приложенной на расстоянии Ив от нейтральной оси ох (рис. 4.8а).
2. После завершения третьего и четвертого этапов: две неразрезные однажды статически неопределимые балки, каждая из которых лежит на трех опорах и нагружена в пролетах 2-3 и 3-4 погонной осевой силой q", приложенной на расстоянии hH от нейтральной ох (рис. 4.8Ь).
После завершения пятого и шестого этапов: неразрезная четырежды статически неопределимая балка, лежащая на шести опорах и нагруженная погонной осевой силой q", приложенной на расстоянии hH от нейтральной оси ох (рис. 4.8с).
Завершающей сформулированную последовательность является четвертая расчетная схема, отвечающая завершению седьмого этапа строительства, которую целесообразно составить в рамках более точной математической модели, чем балочная.
