Авторефераты по всем темам >>
Авторефераты по техническим специальностям
На правах рукописи
НГУЕН ДЫК НГИА
ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ НЕСУЩЕЙ СПОСОБНОСТИ БЕТОННОЙ ГРАВИТАЦИОННОЙ ПЛОТИНЫ
НА СКАЛЬНОМ ОСНОВАНИИ
Специальность:05.23.07 - Гидротехническое строительство.
АВТОРЕФЕРАТ
диссертации на соискание ученой степени
кандидата технических наук
Москва - 2012г.
Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Московский государственный строительный университет
Научный руководитель кандидат технических наук
Толстиков Виктор Васильевич
Официальные оппоненты Марчук Алексей Николаевич
профессор, доктор технических наук, главный научный сотрудник института Физики Земли РАН.
Николаев Дмитрий Валерьевич
кандидат технических наук, ведущий инженер ОАО Атомэнергопроект.
Ведущая организация ОАО НИИЭС
Защита состоится л ноября 2012 г. в часов минут на заседании диссертационного совета Да212.138.03 при Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Московский государственный строительный университет, по адресу: 129337, Москва, Ярославское шоссе, д.26, ауд. №9 Открытой сети образования в строительстве МГСУ.
С диссертацией можно ознакомиться в научной библиотеке ФГБОУ ВПО Московский государственный строительный университет.
Автореферат разослан л____ октября 2012 г.
Ученый секретарь
диссертационного совета Орехов Генрих Васильевич
ОБЩАЯ ХАРАКТЕРИСТИКИ РАБОТЫ
Актуальность.
В современном плотиностроении наблюдаются следующие тенденции: рост числа гравитационных плотин из укатанного бетона с пониженным расходом цемента, увеличение их высоты и обжатие профиля.
Практически доказано, что невозможно достигнуть абсолютную монолитность бетонной кладки, сопряжения её с основанием и избежать появления контактной трещины при допущении растягивающих напряжений под напорной гранью.
Развитие численных методов анализа (совершенствование МКЭ) позволяет разработать методику более достоверного определения напряженно-деформированного состояния (НДС) трещиноватых скальных массивов и взаимодействующих с ними элементов сооружений. Оценивать влияние отдельных крупных трещин, прогнозировать изменение фильтрационных свойств массивов в зависимости от изменения их напряженного состояния.
По существующим нормам проектирования в России, Вьетнаме и других странах оценка несущей способности плотин выполняется по линейной зависимости Кулона - Мора, связывающей средние значения нормальных и предельных касательных напряжений в контактном сечении или в субгоризонтальных трещинах вблизи контакта, что не учитывает сложный характер статической работы системы бетонная гравитационная плотина - основание.
В связи с этим, актуальным является изучение возможных схем потери несущей способности бетонной плотины на скальном основании, их математическое моделирование, прогнозирование и оценка несущей способности системы плотина основание.
Цель и задачи работы:аисследование несущей способности бетонных гравитационных плотин на скальных основаниях на основе анализа фактических предельных состояний комплекса плотина - основание с учетом возможного раскрытия контактного шва, разуплотнения основания, сдвиговых нарушений в сжатой зоне основания и бетона плотины, изменения фильтрационного режима в скальном основании при учете пропускной способности дренажа.
В соответствии с поставленной целью и на основе анализа результатов предыдущих исследований в задачи диссертационной работы входило:
1.аВыполнить анализ имеющихся аварий в бетонных плотинах, приводящих к изменению статических условий работы системы плотина - основание по сравнению с заложенными в проекте; выполнить обзор существующих аналитических и численных исследований статической работы бетонных плотин в предельном состоянии. Проанализировать критерии оценки несущей способности бетонной плотины на скальном основании.
2.аУсовершенствовать программный комплекс CRACK с учётом следующих дополнений: а) за критерий сдвиговой прочности контакта и трещин при высоком уровне сжимающих напряжений принять критерий Фэйрхерста (параболическая огибающая кругов Мора); б) для повышения точности решения физически нелинейных задач алгоритм метода начальных напряжений МКЭ дополнить пошаговой схемой решения при приложении сдвигающих нагрузок; в) при решении задачи в статико-фильтрационной постановке учесть с помощью итерационной процедуры реальную пропускную способность дренажа основания.
3.аПроанализировать методы достижения предельного состояния системы бетонная плотина - скальное основание. Провести серию численных экспериментов по возможности математического моделирования различных схем разрушения плотины на скальных основаниях. Выбрать методику доведения системы плотина - основание до разрушения и оценки несущей способности бетонной плотины на скальном основании.
4.аИсследовать механизм разрушения системы бетонная плотина - блочное скальное основание, сравнить результаты исследований на физической модели с результатами математической модели.
5.аИсследовать изменение фильтрационного режима основания по мере увеличения сдвигающей нагрузки за счет увеличения раскрытия контакта, образования трещин в основании и пропускной способности дренажа.
6.аИсследовать на основе методики научного планирования эксперимента статическую работу при расчётной нагрузке и несущую способность системы плотина - основание при доведении её до разрушения. Получить факторные зависимости, графики и номограммы для определения различных параметров, характеризующих статическую работу плотины и её несущую способность от изменения различных факторов (заложения низовой грани плотины, деформационных и прочностных свойств основания и контакта, параметров подземного контура).
7.аВыполнить расчёты несущей способности плотины ШонЛа (СРВ).
Научная новизна работы данной диссертационной работы характеризуется следующим:
1.аПоказана возможность математического моделирования в рамках решения нелинейных задач МКЭ поведения системы бетонная гравитационная плотина - скальное основание вплоть до разрушения с учётом: снижения сдвиговой прочности контакта и трещин при высоком уровне нормальных напряжений; образования зон разуплотнения и дробления в скальном основании; изменения фильтрационного режима по мере процесса разрушения с учётом изменения проницаемости скального массива и пропускной способности дренажа.
2.аПолучены эпюры фильтрационного противодавления на подошву плотины на различных стадиях продвижения контактной трещины и разуплотнения основания при учете пропускной способности дренажа. Предложены характерные эпюры в зависимости от раскрытия контакта и положения инъекционной и дренажной завес в основании.
3.аПолучены факторные зависимости, графики и номограммы для определения различных параметров, характеризующих статическую работу плотины и её несущую способность от изменения различных факторов (заложения низовой грани плотины, деформационных и прочностных свойств основания и контакта, параметров подземного контура).
4.аВ анализе напряженно-деформированного состояния и несущей способности высокой бетонной плотины Шонла (СРВ) с учетом реального фильтрационного режима и его изменения, с явным моделированием основных крупных трещин и их прочностных и деформационных свойств. В сопоставлении полученных коэффициентов запаса несущей способности с полученными при проектировании по нормам России и США.
Достоверность результатов проведенных исследований подтверждается результатами натурных исследований на ряде бетонных плотин в России, во Вьетнаме и в других странах. Математическая модель также подтверждается сопоставлением полученных результатов с выполненными ранее на кафедре гидросооружений МГСУ исследованиями на физической модели плотины на блочном основании.
Практическая значимость работы состоит в возможности использования разработанной методики и программы для определения прочности и устойчивости системы бетонная гравитационная плотина - скальное основание при наиболее полном учёте факторов, влияющих на несущую способность. Полученные номограммы и графики позволяют при проектировании бетонных плотин оперативно оценивать прочность и устойчивость сооружений на предварительной стадии проектирования.
Внедрение результатов исследований: Оценка несущей способности бетонной гравитационной плотины ШонЛа (СРВ) с использованием предложенной методики и сопоставление результатов с принятыми в проекте.
Апробация работы: Результаты диссертационной работы были доложены: на юбилейной конференции л 80 лет факультету ГСС, Москва, 2010; на Международной межвузовской научно-практической конференции молодых ученых, докторантов и аспирантов Строительство - формирование среды жизнедеятельности, Москва, 2011; на заседании кафедры Гидротехнических сооружений МГСУ, 2012.
Публикации. По теме диссертации опубликовано 3 печатных работ, в том числе 2 в научных журналах, входящих в список ВАК для публикации результатов по кандидатским диссертациям.
Структура и объем работы:
Диссертация состоит из введения, пяти глав, основных выводов, списка используемой литературы из 149 наименований. Полный объём диссертации - 204 страницы, включая 129 страниц текста, 95 рисунков и 36 таблиц.
СОДЕРЖАНИЕ РАБОТЫ
Во введении обосновывается актуальность темы диссертационной работы.
В первой главеаприведен анализ по имеющимся авариям в бетонных плотинах, приводящих к изменению статических условий работы системы бетонная гравитационная плотина - скальное основание по сравнению с заложенными в проекте. Выполнен обзор существующих аналитических и численных исследований статической работы бетонных плотин в предельном состоянии. Проанализированы критерии оценки несущей способности бетонной плотины на скальном основании.
Во второй главеаприведены описания расчетных моделей усовершенствованного программный комплексаCRACK с учётом следующих дополнений:
- Критерием сдвиговой прочности контакта и трещин при высоком уровне сжимающих напряжений принят критерий Фэйрхёрста (параболическая огибающая кругов Мора).
- Для повышения точности решения физически нелинейных задач алгоритм метода начальных напряжений дополнен пошаговой схемой приложения сдвигающих нагрузок.
Рис. 1 - Контактный элемент и схемы итерационного процесса для зависимости σn=f(v), и τ=f(u), где: Rpа-апрочность на растяжение; τр и τrаЦапиковая и остаточная прочность при сдвиге.
Для моделирования трещин и швов используются модифицированные контактные элементы Гудмана (рис.а1).
Критерий сдвиговой прочности трещин используются и линейная зависимость Кулона-Мора (1) и параболическая огибающая кругов Мора предложенная С.аФэйрхерстом (2) позволяющая учесть при высоких нормальных напряженияхσn снижение сдвиговой прочности контакта и трещин.
τпред. = С + σn.tgφ (1),
τпред. = (m-1)[Rр(Rр- σn)]1/2 (σn<0) (2), где τпред - сдвиговая прочность трещины, m=(|Rсж|/Rр+1)1/2, (Rсж<0), C - сцепление трещины, σn-нормальное напряжение, Rсж- прочность на одноосное сжатие материала стенки трещины,Rр - прочность на растяжение материала стенки трещины (контактирующих блоков).
Для описания поведения бетона и пород основания используется математическая модель, позволяющая в рамках плоской задачи учесть полную диаграмму поведения бетона в сложном напряженном состоянии. Модель является обобщением упругой и пластической среды с предельной поверхностью разрушения, определяемой критерием Мора. Зависимость σ=f(ε) принимается в виде идеализированных диаграмм, позволяющих моделировать идеально-пластическое поведение материала: (а) - разупрочняющийся материал с бесконечной хрупкостью, (б) - материал с конечнойхрупкостью (Рис. 2).
Рис. 2 - Идеализированные зависимости осевой деформации ε3 от напряжений σ3 при различных боковых обжатиях Р
В области сжатия и растяжения напряжения ограничены критерием Мора. Принята, при построении модели, параболическая огибающая кругов Мора предложенная С.аФэйрхерстом.
Процесс нагружения.аДля повышения точности решения физически нелинейных задач алгоритм метода начальных напряжений дополнен пошаговой схемой приложения воздействий и осуществления решения.
Сущность этого алгоритма в том, что: нагрузка прикладывается малыми ступенями в той последовательности, в какой происходит реальное нагружение в натуре. Решение в настоящем шаге процесса нагружения - получается суммированием НДС системы предыдущего шага и решения от межшаговой нагрузки. После этого выполняется прочностной анализ и определяются невязки напряжений для итераций внутри шага нагружения.
Существуют два метода достижения предельного состояния системы плотина - скальное основание. Первый - увеличение вызывающих разрушение нагрузок (рис.а3). Второй - уменьшение сдвиговых характеристик контакта. Тогда:
nQ = (Qрас + ΔQ)/Qрас. (3),
nψ = (tgψрас / tgψ) (4),
tgψ = tgφ + c/σсред. (5),
где: nQ - коэффициент перегрузки, nψ- коэффициент уменьшения характеристик контакта, Qрас. - расчётная нагрузка, tgφ, c - расчётные характеристики контакта, σсред.а-асреднее нормальное напряжение в контакте.
Оба метода используются при проведении исследований. Полученные нами результаты показывают, что при низких характеристиках контакта, когда реализуется схема сдвига по контакту, разрушение происходит при близких значениях nQ и nψ (рис.а3). При высоких характеристиках контакта они дают результаты со значительной разницей потому, что схема потери несущей способности системы плотина - основание более сложная. И в практических исследованиях второй подход встречает определённые трудности, связанные с необходимостью одновременного изменения характеристик прочности контакта, бетона и скального массива, прочности в различных системах трещин. В настоящей работе увеличение вызывающих разрушение нагрузок выбрано методом достижения предельного состояния системы бетонная плотина- скальное основание.
Рис. 3 - Схема достижения предельного состояния системы плотинаа-аскальное основание при различных способах доведения до разрушения.
Исследованы возможные схемы разрушения бетонной гравитационной плотины на различных скальных основаниях (табл.а1). Результаты на рис.а4 показывают полученные схемы разрушения системы плотина-основание.
Таблицаа1
Свойства различных оснований и контактов бетон-скала
№ | Физико-механические характеристики | Обозна-чение | Тип основания | Единицы измерения | ||
Жёсткое | Среднейжесткости | Мягкое | ||||
1 | Модуль деформации | Eос | 25000,00 | 15000,00 | 5000,00 | МПа |
2 | Прочность нарастяжение | Rр | 0,25 | 0,17 | 0,10 | МПа |
3 | Прочность насжатие | Rсж | 20,00 | 12,50 | 5,00 | МПа |
4 | Коэффициент Пуассона | νос | 0,24 | 0,24 | 0,24 | а |
5 | Объёмный вес | γос | 25,00 | 25,00 | 25,00 | кН/м3 |
а | а | Контакт бетон-скала | ||||
6 | а | tgφ | 2,00 | 1,50 | 1,00 | а |
7 | Уголтрения | φ | 63,43 | 56,31 | 45,00 | градусов |
8 | Сцепление | С | 3,00 | 2,00 | 1,00 | МПа |
9 | Нормальная жёсткость | Kn | 25а000,00 | 25а000,00 | 25а000,00 | МПа/м |
10 | Касательная жёсткость | Ks | 10а000,00 | 10а000,00 | 10а000,00 | МПа/м |
На мягком основании система плотина-основание разрушилась по схеме близкой к плоскому сдвигу. Однако, это не плоский сдвиг по контакту, а более сложная схема. При этом появилась зона разуплотнения (разрушение за счет главных растягивающих напряжений) основания под напорной гранью, раскрытие контакта, сдвиговое разрушение части контакта, зона дробления (разрушения) основания под низовой гранью плотины. Часть основания сдвигалась вместе с плотиной в нижний бьеф (Схема 1).
Плотина, на жёстком прочном основании, в предельном состоянии также имела зону разуплотнения основания под напорной гранью, раскрытие контакта, сдвиговое разрушение части контакта, разрушение бетона вблизи низовой грани от сжатия. Можно сказать, что реализовалась схема близкая к опрокидыванию относительно низовой кромки подошвы плотины (Схема 2).
Механизм разрушения плотины построенной на основании средней жесткости более сложен. Им является сложная комбинация отдельных разрушений: сдвиг плотины на части контакта с основанием, разрыв скального массива под напорной гранью и дробление его под низовой гранью плотины, разрушение бетона от сжатия в низовом клине плотины. В момент предельного состояния плотина с частью основания опрокинулась относительно оси в глубине основания (Схема 3).
Тип основания | Схема деформированного состояния в момент разрушения | Состояние элементов |
Мягкое | ||
Среднеей жёсткости | ||
Жёсткое |
Рис. 4 - Состояние системы бетонная плотина - скальное основание при предельном состоянии
Состояние контактных элементов: 1- упругая работа элемента; 2- элемент в условиях нормального отрыва; 3- элемент в пластическом состоянии по условиям сдвига.
Состояние сплошных элементов: 4- в элементе по одному из направлений была превышена прочность на растяжение; 5- в элементе в условиях двухосного растяжения была превышена прочность на растяжение по двум направлениям; 6- произошло разрушение от сжатия при растягивающих напряжениях по другой оси; 7- разрушение произошло в условиях двухосного сжатия.
Разрушение системы плотина-основание по схеме 3 упоминается в IV приложении СНиПаа2.02.02-85, однако формулы, определяющие положение оси опрокидывания, и схематизация работы массива в виде равномерной эпюры на участке смятия несколько условны и численный подход для конкретного случая будет необходимым.
На рис.а5 построены зависимости, иллюстрирующие возможные схемы разрушения системы плотина-основание с учетом прочности бетона плотины, массива основания и их контакта. На мягком основании (K2а<аnQпред.), эта система разрушилась по схеме плоского сдвига. На основании средней жесткости (K2а>аnQпред.), эта система разрушилась за счёт разрушения и бетона и скалы. Затем эта система разрушилась за счёт разрушения бетона в зоне низовой клина плотины (nQпред. возрастает незначительно).
Рис.а5 - Зависимость коэффициентов устойчивости от tg
Таким образом, можно предложить методику определения несущей способности системы бетонная гравитационная плотинаа-аскальное основание косвенно через коэффициент перегрузки. Сущность этой методики заключается в следующем:
- Находится момент предельного состояния системы бетонная плотина - скальное основание путём увеличения порождающих разрушение нагрузок при не изменяющих остальных нагрузках до момента разрушения этой системы.
- Коэффициент перегрузки, при котором система плотина-основание достигнет предельного состояния является общим коэффициентом устойчивости этой системы.
Для выполнения исследований была создана программа автоматического построения сеток конечных элементов для различных видов основания: основание, рассматриваемое как сплошная среда, блочное основание (основание с регулярными системами трещин рис.а6), что позволило существенно снизить трудоемкость подготовки расчетных схем и исключить ошибки в массивах исходной информации.
Рис. 6 - Слева: Интерфейс программы Block Foundation подготовки сетки конечных элементов на блочном скальном основании. Справа: Сетка конечных элементов из 2806 узлов; 840 контактных элементов; 1582 обычных элементов.
Рассмотренная выше методика и алгоритм расчета сооружений, с учетом контактных взаимодействий по имеющимся нарушениям оплошности, были использованы при исследовании статической работы облегченной бетонной плотины расположенной на блочном скальном основании. Эти результаты позволили апробировать расчётную модель.
Рис. 7 - Сравнение картин разрушения на физической и математической моделях
Результаты расчета сопоставлялись с результатами физического моделирования выполненного ранее в масштабе 1:200 на кафедре гидросооружений МГСУ (рис.а7). Результаты, полученные на физической и математической моделях достаточно хорошо совпадают. Модели разрушаются при близких коэффициентах перегрузок (nQ=1,38 - физическая модель и nQ=1,40 - математическая). Математическая модель хорошо отображает механизм разрушения системы плотина-основание в виде предельного поворота.
В третьей главеарассматривается методика, и приводятся результаты решения задач в статико-фильтрационной постановке. Эти задачи рассматривались некоторыми авторами при расчётной нагрузке, когда изменение условий работы системы плотина - основание не очень значительно, и дренаж основания считается с бесконечной пропускной способностью, т.е. напор в дренаже постоянным. При перегрузке НДС этой системы значительно изменяется с появлением раскрытия контакта, разрушения основания, что приводит к росту фильтрационного расхода в дренаже. Когда этот расход превышает пропускную способность дренажа, тогда напор в дренаже увеличивается, т.е. меняются граничные условия решения задачи. Поэтому необходимо добавить процедуру нахождения напора в дренаже с конечной пропускной способностью. Эта процедура показана на рис.а8.
- при определённом значении nQ, задаем начальный напор (Нi) в дренаже, равный полученному напору в предыдущем шаге процесса нагружения.
- на шаге i-ом фильтрационной процедуры решив статико-фильтрационную задачу с напором в дренаже Нi=const, получим расход в дренаже по формуле:
(6), где n - количество окружающих дренаж элементов,
Kj - проницаемость окружающего элемента j
Нj - действующий напор на элементе j;
Li - длина конечного элемента j;
- с таким расходом Qi по кривой пропускной способности дренажа, напор в дренаже должен быть HiТ по формуле:
(7), где m- коэффициент расхода в дренажной скважине; - площадь поперечного сечения скважины; g - ускорение свободного падения; k - количество скважин на 1 пог. м.
Точность решения до i-ого шага:
(8),
- напор для следующего шага:
(9).
Процедура закончится, когда di , где - заданная точность процедуры.
| 1 - пропускная способность дренажа 2 - зависимость расхода в дренаже от заданных напоров | |
Рис. 8 - Процедура для определения напора в дренаже | Рис. 9 - Противодавление по подошве плотины при различных значениях nQ |
Решив статико-фильтрационную задачу для плотины на сплошном скальном основании и на блочном скальном основании можно получить противодавление по подошве плотины (рис.а9) в зависимости от nQ и соответствующих нарушений подземного контура. Если рассмотрим эти эпюры вместе с раскрытием контакта, получим результат, представленный на рис.а10, где λт - относительная длина раскрытия контакта, λц - относительное расстояние от напорной грани до оси цементационной завесы, λда-аотносительное расстояние от напорной грани до оси дренажа.
| λт = 0 | 0а<аλт<аλц | λцааλт λд | λтλд |
Рис. 10 - Противодавление по подошве плотины в зависимости от раскрытия контактаи рекомендуемая нами эпюра для статических расчётов
В четвёртой главе применена теория планирования эксперимента при определении несущей способности бетонной гравитационной плотины на скальном основании. Выполнены два плана эксперимента для плотины высотой Hа=а100ам: первый план выполнен при расчётной нагрузке, а второй план для определения предельного состояния системы плотина - основание. Использована квадратичная модель с числом экспериментов Nа=а2n+2n, где n - количество варьируемых факторов. При обоих планах варьируются по три фактора, то Nа=а14 расчётов.
В первом рассмотрены следующие факторы: Уклон низовой грани m2, варьируется с 0,60,8; Модуль деформации основания: 525аГПа; Относительное расстояние от напорной грани до оси цементационной завесы λц: 0,10,3 В качестве функций отклика рассматривались: λр,та-аотносительная длина раскрытия контакта бетон-скала (Y1); σпа-амаксимальные главные сжимающие напряжения в бетоне (Y2); σо - максимальные главные сжимающие напряжения в основании (Y3); Kу - коэффициенты устойчивости сооружения на сдвиг по контакту со сдвиговыми характеристиками по полевым испытаниям (Y4); Kу,i - коэффициенты устойчивости сооружения на сдвиг по контакту с сдвиговыми характеристиками по СНиПу 2.02.02-85 (Y5Y7), где i = 1: f1а=а0,80, c1а=а0,2аМПа, i = 2: f2а=а0,85, c2а=а0,3аМПа, iа=а3:аf3а=а0,95, c3а=а0,4аМПа. Полученные в результате эксперимента значения функций отклика Y приведены в табл.а2.
Таблицаа2
Значения функций отклика
№ экспер. | Y1 | Y2 | Y3 | Y4 | Y5 | Y6 | Y7 |
λр.т. | σп. (МПа) | σо. МПа) | Kу | K1 | K2 | K3 | |
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
1 | 0.600 | 17.53 | 13.60 | 1.34 | 0.81 | 0.90 | 1.03 |
2 | 0.700 | 22.02 | 16.37 | 1.19 | 0.77 | 0.84 | 0.96 |
3 | 0.625 | 14.24 | 11.85 | 3.19 | 0.85 | 0.93 | 1.07 |
4 | 0.750 | 19.97 | 15.27 | 2.70 | 0.83 | 0.90 | 1.02 |
5 | 0.000 | 5.27 | 3.66 | 3.26 | 1.65 | 1.89 | 2.22 |
6 | 0.050 | 5.50 | 3.84 | 2.96 | 1.46 | 1.68 | 1.98 |
7 | 0.125 | 3.22 | 2.79 | 6.90 | 1.39 | 1.59 | 1.87 |
8 | 0.150 | 3.34 | 2.88 | 6.88 | 1.40 | 1.60 | 1.88 |
9 | 0.725 | 18.36 | 14.28 | 2.03 | 0.84 | 0.91 | 1.04 |
10 | 0.075 | 3.88 | 3.08 | 5.24 | 1.53 | 1.75 | 2.05 |
11 | 0.100 | 6.72 | 4.96 | 2.60 | 1.33 | 1.52 | 1.78 |
12 | 0.400 | 5.87 | 5.38 | 4.14 | 0.85 | 0.98 | 1.14 |
13 | 0.375 | 6.27 | 5.39 | 3.08 | 0.91 | 1.05 | 1.22 |
14 | 0.250 | 5.51 | 4.57 | 3.88 | 1.18 | 1.34 | 1.56 |
Функции отклика представляются в виде:
(10), где: xi (i =13) - исследуемые факторы. С использованием матричного уравнения метода наименьших квадратов были получены значения коэффициентов регрессии для выбранных функций отклика (табл.а3)
Таблицаа3
Значения коэффициентов уравнения регрессии
Функция откликаа | Y1 | Y2 | Y3 | Y4 | Y5 | Y6 | Y7 |
в центре плана | 0.275 | 5.588 | 4.686 | 3.640 | 1.090 | 1.240 | 1.450 |
b0 | 0.295 | 5.960 | 5.021 | 3.466 | 1.088 | 1.242 | 1.446 |
b1 | -0.297 | -7.090 | -5.512 | 1.479 | 0.333 | 0.403 | 0.488 |
b2 | 0.110 | -1.040 | -0.424 | 1.246 | -0.070 | -0.083 | -0.099 |
b3 | 0.020 | 0.980 | 0.563 | -0.016 | 0.003 | 0.000 | -0.001 |
b4 | 0.105 | 5.159 | 3.660 | 0.169 | 0.097 | 0.088 | 0.099 |
b5 | -0.045 | 0.335 | 0.145 | -0.096 | 0.002 | 0.008 | 0.014 |
b6 | 0.017 | -0.070 | -0.043 | 0.014 | -0.043 | -0.047 | -0.056 |
b7 | 0.022 | 0.141 | 0.128 | 0.525 | -0.053 | -0.059 | -0.069 |
b8 | -0.018 | -1.234 | -0.741 | 0.040 | -0.015 | -0.014 | -0.014 |
b9 | 0.003 | 0.142 | 0.069 | -0.007 | 0.027 | 0.031 | 0.034 |
b10 | -0.009 | -0.167 | -0.092 | 0.077 | 0.022 | 0.024 | 0.029 |
По полученным результатам построены графики и номограммы для определения различных параметров, характеризующих статическую работу плотины и её несущую способность от изменения различных факторов.
Рис.а11 - Зависимости относительного раскрытия контакта от модуля деформации основания при X3а=а0 | Рис.а12 - Номограмма для определения относительной длины раскрытияконтакта | ||
Рис.а13а-аМаксимальные главные сжимающие напряжения в бетоне в зоне низового клина плотины (H = 100 м) | Рис.а14а-аМаксимальные главные сжимающие напряжения в основании в зоне низового клина плотины (Hа=а100ам) | ||
Рис.а15 - Коэффициент устойчивости плотины против сдвига по контакту, характеристики которого по полевым экспериментам | Рис.а16 - Коэффициент устойчивости плотины против сдвига по контакту, характеристики которого по СНиПу | ||
Во втором плане также рассмотрены три фактора, однако, вместо фактора модуля деформации рассматривалась прочность на сжатие скалы, варьируемая с 5аМПа до 20аМПа. При этом можно определить предельный коэффициент перегрузки nQпред через формулу (11) или по ниже номограммам. Где xi = (-11), i=13.
(11)
Рис.17 - Зависимости коэффициента устойчивости от прочности на сжатие основания при (λда=а0,1) | Рис.18 - Номограмма для определения коэффициента устойчивости плотины, определенного косвенно через nQпред |
Вапятой главеабыла оценена несущая способность бетонной гравитационной плотины ШонЛа возведенной из укатанного бетона во Вьетнаме.
Рис. 19 - Расчётный профиль и отметки плотины ШонЛа: PMF = 218,07ам, УЗВа=134ам, уровень воды в нижнем бьефе соответственно PMF 136,95ам.
Плотина ШонЛа имеет высоту 138ам, построена на реке Да в верхнем бьефе ГЭС ХоаБинь. Расчётный профиль плотины ШонЛа приведен на рис.а19. В проекте несущая способность оценена по Вьетнамским и Российским нормам; и по нормам США. В данной работе выполнены расчёты для определения несущей способности расчётного профиля и сечения D14 по предложенной методике, описанной во второй главе, при сочетании Р3 (PMF) на основе линейной механики разрушения и МКЭ.
Рис. 20 - Зависимости параметров прочности плотина ШонЛа для сочетания нагрузок Р3: слева - по линейной механике разрушения; права -с учётом податливости основания.
Рис. 21 - Состояние системы плотина-основание при nQ=1,32 - предельном состоянии с данными по Вьетнамским и Российским нормам
Рис. 22 - Напряжённо-деформированное состояние системы плотина-основание в предельном состоянии (nQ=1,85) с данными по нормам США
Рис. 23 - Сетка конечных элементов сечения D14: 3603 узлов, 3571 обычных элементов, 140 контактных элементов
Рис. 24 - Состояние элементов сечения D14 при предельном состоянии nQа=а1,5 с данными по Российским нормам (состояния элементов см. рис.а4) | Рис. 25 - Состояние элементов сечения D14 при предельном состоянии nQа=а2,5 с данными по нормам США (состояния элементов см. рис.а4) |
Таблица1
Сравнение расчётных результатов (МКЭ) с проектируемыми параметрами
По Российским нормам | По нормам США | ||||
по проекту | по расчёту | разница% | по проекту | по расчёту | разница % |
Расчётное сечение | |||||
1,524 | 1,32 | -13,38 | 2,86 | 1,85 | -35,31 |
Сечение D14 | |||||
1,61 | 1,50 | -6,83 | 3,04 | 2,50 | -17,76 |
ЗАКЛЮЧЕНИЕ
- Анализ аварий бетонных гравитационных плотин на скальных основаниях показывает, что характерными их причинами являются, нарушение прочности на растяжение в зоне контакта со стороны напорной грани с раскрытием контактного шва и образованием трещин в массиве основания, сопровождаемое увеличением фильтрационного противодавления. При этом возможно нарушение прочности на сдвиг и сжатие в зоне низового клина плотины, приводящие к потере несущей способности. Натурными исследованиями на ряде плотин выявлено значительное раскрытие контактного шва и явление трещинообразования в глубине основания (разуплотнение основания), сопровождающиеся увеличением эпюры фильтрационного противодавления до 100% действующего напора по длине раскрытого контакта. Это может привести к существенному ухудшению схемы работы сооружения по отношению к заложенной в проекте.
- Накопленный опыт математического моделирования с использованием МКЭ позволил построить расчетный аппарат для моделирования поведения системы плотина-основание при учете почти всех особенностей их работы на уровне имеющихся экспериментальных данных, а именно:
-апредставимость в пределах расчетной схемы МКЭ с необходимой подробностью конструктивных особенностей сооружения и структуры массива горных пород основания;
-авозможность численного моделирования в явном виденарушений сплошности расчетной области (крупные трещины, технологические и конструктивные швы, зоны контактов, разломы и т.д.) при решении контактной задачи МКЭ в нелинейной постановке;
-аобеспечение расчетов на температурные и квазидинамические воздействия;
-адостижимость требуемой точности в определении деформаций и напряжений, что особенно важно при решении нелинейных задач;
-арешение связанной задачи статики и фильтрации с учётом связи между НДС скального массива и его проницаемостью, влиянием раскрытия контакта и пропускной способности дренажа на формирование фильтрационного противодавления.
Используемые при этом упругопластические модели поведения бетона и скалы основания и их программная реализация позволяют производить расчеты плотин на скальных основаниях вплоть до исчерпания несущей способности.
- Сравнение результатов экспериментальных исследований разрушения системы бетонная плотина - скальное основание с результатами расчётов показало, что результаты, полученные на физической и математической моделях достаточно хорошо совпадают. Модели разрушаются при близких коэффициентах перегрузок (nQ=1,38 - физическая модель и nQ=1,40 - математическая). Математическая модель хорошо отображает механизм разрушения системы плотина-основание в виде предельного поворота. Это даёт возможность при постановке широких исследований использовать расчетный метод и существенно сократить трудоёмкие экспериментальные исследования.
- При реализации численного эксперимента с использованием теории планирования при расчётной нагрузке получены зависимости относительной длины раскрытия контакта, главных сжимающих напряжений и коэффициентов устойчивости от основных параметров плотины - уклона низовой грани, податливости основания и расстояния от напорной грани до оси цементационной завесы с учётом изменения противодавления воды по подошве плотины в зависимости от длины контактной трещины, разуплотнения основания и пропускной способности дренажа основания. Выполнен анализ полученных зависимостей и выявлены особенности формирования напряженно-деформированного состояния и параметров раскрытия контактного шва от каждого из принятых факторов, а также от их комплексного воздействия для рассмотренных расчётных схем. Например, показано, что раскрытие контакта заметно превосходит от 1,5 до 3,5 раз глубину зоны растягивающих напряжений, которая принята в действующих нормах для оценки местной прочности. При пересечении трещиной оси дренажа это отличие может доходить до 45 раз.
- Доведение рассмотренных вариантов плотин на различных по жёсткости и прочности основаниях показало, что практически реализуются две схемы разрушения системы плотина - основание. Схема плоского сдвига (с раскрытием контакта), характерна для плотин с коэффициентами сдвига в контакте tgψ 1,1. При более высоких параметрах сдвига схемой разрушения становится схема предельного поворота с разрушением основания иногда и бетона в сжатой зоне плотины и массива основания. А сдвиг вторичен, как окончательный механизм перемещения сооружения после разрушения. Эта схема хорошо согласуется с результатами исследований других авторов выполненных как на моделях плотин, так и при натурных испытаниях.
По предлагаемой методике была исследована несущая способность бетонных гравитационных плотин различного профиля, при изменении модуля деформации скального основания, прочности бетона, скального массива и контакта бетон-скала. Для удобства инженерного применения полученных результатов приведены формулы, графики и номограммы, позволяющие простейшими арифметическими действиями решать инженерные задачи по определению несущей способности бетонных сооружений для принятых условий.
- В проекте несущая способность плотины ШонЛа была оценена по прочностным характеристикам контакта и основания, соответствующим различным нормам Российским (Вьетнамским), и Американским. На основе этих данных, оценка несущей способности плотины ШонЛа выполнена на основе расчётов МКЭ по предложенной методике с учётом изменения статических условий работы: учтено раскрытие контактного шва, образование трещин в основании и бетоне плотины, разрушение основания под низовым клином плотины, изменение фильтрационного режима. Несущая способность плотины ШонЛа по предложенной методике ниже, чем в проекте на 13-35% в зависимости от характеристик основания и схемы потери несущей способности. Меньшее значение 13% соответствует прочностным показателям, принятым по Российским нормам и схеме разрушения близкой к схеме плоского сдвига. Схема предельного поворота дает более существенное отличие до 35%.
Основные положения диссертации опубликованы
в следующих работах автора:
1)аТолстиковаВ.В., Нгуен Дык Нгиа. Численное моделирование разрушения системы бетонная гравитационная плотина - блочное скальное основание //Гидротехническое строительство №10, 2011 - стр. 6 11.
2)аТолстиковаВ.В., Нгуен Дык Нгиа. Численное исследование возможных схем разрушения системы бетоннаяагравитационная плотина - скальное основание //Вестник МГСУ №5, 2011 - стр. 41 47.
3)аТолстиковаВ.В., Нгуен Дык Нгиа. Численное исследование возможных схем разрушения системы бетоннаяагравитационная плотина - скальное основание//Сборник трудов Четырнадцатая Международная межвузовская научно-практическая конференция молодых учёных, докторантов и аспирантов, Апрель 2011аЦастр.а678аа684.
Авторефераты по всем темам >>
Авторефераты по техническим специальностям