Собственные значения и формы колебаний арочных плотин и их напряженно-деформированное состояние от действия сейсмических сил

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

АВТОРЕФЕРАТ

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

Собственные значения и формы колебаний арочных плотин и их напряженно-деформированное состояние от действия сейсмических сил

Специальность : 05.23.07 - Гидротехническое строительство.

кандидата технических наук

Нгуен Куанг Кыонг

Москва - 2006

Работа выполнена в Государственном образовательном учреждении высшего профессионального образования

Московском государственном строительном университете

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

Официальные оппоненты доктор технических наук, профессор Гольдин Александр Львович

кандидат технических наук Доронин Федор Леонидович

Ведущая организация Центр службы геодинамических наблюдений в энергетической отрасли (ЦСГНЭО) - филиал ОАО «Инженерный центр ЕЭС»

Защита состоится « » _________ 2007 года в часов минут на заседании диссертационного совета Д212.138.03 при ГОУВПО Московском государством строительном университете, по адресу: 107066, Москва, ул. Спартаковская, д.2/1, ауд. _____

С диссертацией можно ознакомиться в научной библиотеке ГОУ ВПО Московского государственного строительного университета.

Автореферат разослан «____»_________ 200 г.

Ученый секретарь

диссертационного совета Орехов Г.В.

арочный плотина сейсмический колебание

Общая характеристика работы

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

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

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

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

Цель и задачи работы.

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

· разработка алгоритма программы автоматизированного расчета на ЭВМ арочной плотины с основанием в пространственной постановке с учетом сейсмических воздействии;

· определение собственных колебаний по предлагаемой методике комплекса: арочная плотина и основание.

· построение номограмм периодов собственных колебаний от различных наиболее существенно влияющих факторов.

· определения напряженно-деформированного состояния арочных плотин вместе с основанием от сейсмических воздействий, задаваемых характерной трехкомпонентной акселерограммой;

· исследование зависимости максимальных напряжений арочной плотины от направления сейсмических волн.

Достовернось работы определяется глубиной проверки МКЭ в течении 40 лет проверкой созданной программы решением текстовых задач и сопоставлением собственных значений, арочной плотины Ингурской ГЭС с данными натурных наблюдений, которые подтвердили достоверность полученных результатов.

Научная новизна работы данной диссертационной работы заключается в:

· исследовании 100 первых собственных значений и форм колебаний арочных плотин совместно с основанием в пространственной постановке методом конечных элементов.

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

· исследовании напряженно-деформированного состояния арочных плотин с основанием при землетрясении 9 баллов.

· анализе влияния угла подхода сейсмической волны к створу плотины на напряженно-деформированного состояние плотины.

Практическая значимость и реализация результатов исследований определяется расширение строительства арочных плотин в условиях СРВ и вообще в Юго-западной Азии. Весь этот регион сейсмически активен и строительных арочных плотин требует глубоких исследований НДС плотин при действии сейсмически сил.

Внедрение результатов исследований:

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

Апробация работы:

Основные результаты были доложены в 2-х статьях:

1) д.т.н., проф. Рассказов Л. Н., Нгуен Куанг Кыонг - Вестник МГСУ №2, 2006 - стр. 28 ч 42;

2) д.т.н., проф. Рассказов Л. Н., Нгуен Куанг Кыонг - журнал Гидротехническое строительство №9, 2006 - стр. 32 ч 36;

и будут использованы при создании учебника «Гидротехнические сооружения» - 2-я редакция.

На защиту выносятся:

· алгоритм и программа расчета собственных значений и форм колебаний.

· результаты численных исследований собственных периодов и форм колебаний арочных плотин.

· построение номограмм по определению собственных периодов различных арочных плотин.

· результаты численных исследований НДС арочных плотин при землетрясении 9 баллов.

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

Структура и объем работы:

Диссертация состоит из введения, пяти глав, основных выводов, списка используемой литературы из 78 наименований. Работа изложена на 116 страницах машинописного текста, содержит 142 рисунка и 24 таблицы.

Содержание работы

Во введении обосновывается актуальность темы диссертационной работы.

Первой глава посвящена обзору исследование по применению МКЭ для расчета плотин на сейсмические воздействия.

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

Динамический метод, основанный на решении полных динамических уравнений совместно с физическими уравнениями состояния материалов, открывает более широкие пути для исследования сложных конструкций гидротехнических сооружений. Этот метод дает возможности расчета сооружений динамическим методом не только в плоской, но и в пространственной постановках. Метод основывается на решении уравнения свободных колебаний системы с n-степенями свободы без затухания, которое записывается в виде: л{Х)_i=[D]{Х}_i . К численным методам решения этого уравнения относятся методы: метод Стодолы, метод Данилевского, метод Хольцера, метод Якоби, метод Шварца-Пиккарда, "степенной метод" и др.

На кафедре Гидротехнических Сооружений Московского инженерно-строительного института (ныне МГСУ) под руководством Л.Н. Рассказова была решена задача сейсмостойкости грунтовых плотин в плоской постановке динамическим методом. Были найдены формы собственных колебаний грунтовых плотин итерационным методом, используя характеристическое уравнение. В работе определены пульсации динамических напряжений в сооружении во время землетрясения. Используя энергетическую модель грунта, предложенную Л.Н. Рассказовым, были найдены остаточные перемещения, возникающие в плотине после прохождения сейсмической волны.

В работе Л.Н. Рассказова приводится исследование влияние конечной скорости распространения сейсмических волн в основании сооружения на его сейсмонапряженное состояние. Скорости распространения сейсмических волн изменялись от х= (платформенный эффект) до х=1500м/с.

Во второй главе описываются теоретические основы и методика решения задачи собственных значений, собственных форм колебаний и НДС арочных плотин при сейсмических воздействиях, показывается решение тестовых задач.

Решение задачи о собственных колебаниях сводится к решению динамического уравнения свободных колебаний, которое записывается в виде:

- характеристическое уравнение. После преобразования уравнения получается вид: л{Х} = [Н] {X}

где л =1/щ²; [H] =[L]^-1 [M] [L^T]^-1; {X} =[L^T] {r_o}; [L][L]^Т =[К];

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

Метод Якоб дает всё собственные значения и вся формы, но требует большего объема времени.

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

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

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

где [М], [С], К] - соответственно матрицы масс, затухания, жесткости; {г} - вектор смещения узловых точек относительно неподвижного основания; {Мх}, {Му}, {Мz} - столбцы масс, определяющие инерционные нагрузки по координатным осям X, У, Z; U_ox(t), U_ox(t), U_ox(t), - компоненты сейсмических ускорений по координатным осям X, У, Z.

После преобразования получим:

где щ - частота колебания; о - коэффициент затухания;

[М]* - диагональная матрица, значения которой равны:

M_ij^*={Х}_i^T[М]{Х}_j ={Х}_i^T[М]{Х}_i, при i=j

Для решения этого уравнения могут применяться различные численные методы, в частности метод Рунге-Кутта второго порядка, который и используется в представленной работе.

Решение тестовых задач:

Для проверки работы программы были проведены тестовые расчёты. Расчетная схема первой задачи является «бетонный» столб высотой 100м квадратного поперечного сечения (20х20м). Модуль упругости 200000кг/см² при коэффициенте Пуассона 0,2. Объемный вес был принят 2т/м³. Решалась тестовая задача двумя методами: Якоби и Шварца. В результате получили хорошее совпадение. По 1-ой и 2-ой формам столб колеблется по осям X и Y с периодами Т₁=Т₂=0,83с, по 3-ей форме столб вращается вокруг свой оси с периодом Т₃=0,19с.

рис. 1 Колебание столба в плане по 3-ей форме. Т₃=0,19с.

Для того, чтобы проверить алгоритм расчета программы по напряжениям была приведена специальная акселерограмма: горизонтальный составляющий , в момент с 2 до 14 секунд =-g. Можно считать, что при столб свободно падает. При свободном падении столб переходит в невесомое состояние, и напряжения в нем должны быть равны нулю. Из рис. 2 видно, что к моменту Т=2,4с напряжения в столбе у_z>0.

рис. 2 Пульсация напряжения в элементах при свободном падений.

В третьей главе решается задача собственных значений и форм колебаний арочных плотин различных типов и были построены номограммы для определения периоды собственных колебаний по 6-ти формам различных арочных плотин. Задача собственных значений и форм колебаний решалась в пространственной постановке с использованием шестигранных элементов. Характеристическое уравнение решалось методом Шварца-Пиккарда с использованием метода Якоби для сопоставления с методом Шварца-Пиккарда. Для получения функциональной зависимости в виде полинома от выбранных факторов использовался полнофакторный вариант из метода планирования эксперимента.

Цилиндрическая плотина Плотина в треугольном створе

рис. 3 Расчетные схемы.

В качестве факторов на основе анализа выбраны:

1) коэффициент стройности - от 0,15 до 0,25;

2) модуль деформируемости основания Е -от 50000 до 500000 кг/см²;

3) коэффициент створа ( L/H ) - от 1 до 3; где L -ширина створа

4) высота плотины Н= - от 70 до 200 м.

Таким образом, для решения поставленной задачи требовалось решить 2⁴ задач, т.е. 16 задач и одну задачу в центре плана при = 0,2, =275000 кг/см², = 2 и =135,0 м. Нормированные значения факторов в центре плана равны х₁= х₂ = х₃ = х₄ = 0.

Соответственно план для оценки собственных форм и собственных значений предоставлен в виде таб. 1. Решалась также задача на нулевом уровне для проверки адекватности. В нижней таблице приведены коэффициенты полиномов типа :

T_i=b_0i+b_1ix₁+b_2ix₂+b_2ix₃+b_4ix₄+b_12ix₁x₂+b_13ix₁x₃+b_14ix₁x₄+b_23ix₂x₃+b_24ix₂x₄+b_34ix₃x₄+b_123ix₁x₂x₃+b_124ix₁x₂x₄+b_134ix₁x₃x₄+b_234ix₂x₃x₄+b_1234ix₁x₂x₃x₄ (1)

Функции типа (1) были построены для шести собственных значений. Проверка адекватности (правый нижний угол 0) показала, что в центре плана наименьшее расхождение с расчетом (3%) имеет полином Т₁ =f[х₁, х₂ …]. Наибольшее расхождение при прогнозе - 9% (Т₆ =f[х₁, х₂ …] - строка test). При этом закономерности расхождения от номера формы нет.

рис. 4 Номограмма для определения периода колебаний цилиндрических арочных плотин

а) Периоды 1-ой, 2-ой, 3-ей форм

б) Периоды 4-ой, 5-ой, 6-ой форм

рис. 5 - Формы колебаний цилиндрической плотины

Вариант 1 (t=0.25, E=500 000 кг/см², L/H=3, H=200 м)

а) 1-ая формаб) 2-ая форма

в) 3-яя формаг) 4-ая форма

таб. 1 План расчетов цилиндрических арочных плотин

x1: коэффициент стройности 0.15ч0.25x2: модуль основания 50.000ч500.000 kg/cm2

x3=L/H : коэффицент створа 1ч3x4: высота 70ч200 m

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

а) Периоды 1-ой, 2-ой, 3-ей форм;

б) Периоды 4-ой, 5-ой, 6-ой форм.

таб. 2 План расчетов арочных плотин в треугольном створе.

x1: коэффициент стройности 0.15ч0.25x2: модуль основания 50.000ч500.000 kg/cm2

x3=L/H : коэффицент створа 1ч3x4: высота 70ч200 m

Размещено на http://www.allbest.ru/

Видимо, эти данные можно считать приемлемыми для оценки периода собственных колебаний. Коэффициент b₁ очень мал, как, впрочем, и коэффициенты взаимодействия b₁₂, b₁₃, b₁₄, b₁₂₄, b₂₃₄, b₁₂₃₄, т.е. все коэффициенты взаимодействия первого фактора (коэффициент стройности) с остальными. Из этого следует, что коэффициент стройности в выбранном диапазоне вариации влияет на периоды собственных колебаний в пределах первых шести форм очень мало.

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

На рис. 5, а и б показаны четыре первых формы колебаний в трехмерной постановке для варианта 1 (таб. 1).

Вторая форма симметрична. Формы три и четыре представляют несомненный интерес. Форма три симметричная, а форма четыре несимметричная, но характер колебаний напоминает симметричную. Третья форма (рис. 5, в) характеризуется вертикальными колебаниями и встречными колебаниями в плоскости от борта к борту. Четвертая форма носит сложный характер (рис. 3, г): имеют место существенные вертикальные колебания и переход через “0” гребня плотины в плане.

Аналогичные исследования были проведены и для треугольного створа. Те же факторы, Расчетная схема представлена на рис. 3. План расчетов арочных плотин в треугольном створе (плотина с постоянным углом) представлен в таб. 2.

Результаты во многом аналогичны: еще более мало влияние коэффициента стройности и всех взаимодействий с его участием. В связи с этим варианты плотин, отличающиеся только коэффициентом стройности, имеют очень близкие собственные значения по первой и второй, а иногда и третьей форме (см. отклики в таблице 2). Улучшилось условие адекватности. В основном расхождение 34%, и только шестая форма имеет отклонение 7%.

Наибольший период собственных колебаний равен Т₁ = 1,127 с в третьем и четвертом вариантах, а минимальное значение 0,124 с в 13-м варианте. По результатам этих исследований построена номограмма определения Т₁f(x_i)Т₆ f( x_i) (рис. 6, а и б). Номограммы упростились, т.к. выпал из рассмотрения коэффициент стройности. Теперь имеют место две шкалы и три поля (для Т₁, Т₂ и Т₃ на рис. 6, а и аналогично на рис. 6, б). Очень велико влияние модуля упругости основания (левая шкала).

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

В данной работе рассмотрим только арочные плотины.

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

Решение поставленной задачи осуществлялось на основе рассмотрения двух типов арочных плотин: с постоянным центральным углом и радиусом (цилиндрическая арка) и с постоянным центральным углом (в треугольном створе).Задачи решались в пространственной постановке совместно с основанием. Обе решенные задачи представлены на рис. 3.

Контрольные элементы с напорной стороны, в которых рассматривались напряжения, обозначены как A, B, C, D в случае цилиндрической плотины и A, B, C в треугольном створе (рис. 8). Аналогичные элементы с низовой грани обозначены A, B, C, D и соответственно A, B, C.

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

а) горизонтальный состав б) вертикальный состав

рис. 7 акселерограмма Паркфилд

Направления сейсмического воздействия представлены на рис. 8. Угол менялся с интервалом 15 от 0 (направление сейсмических волн вдоль русла реки) и до 90 (направление вдоль створа).

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

- 1-1: цилиндрическая плотина - Н=200м; t=0,25; l=3; E_осн.=500 000кг/см²;

- 1-2: цилиндрическая плотина - Н=200м; t=0,25; l=1; E_осн.= 50 000кг/см²;

- 2-1:плотина в треугольном створе-Н=200м; t=0,25; l=3; E_осн.=500 000кг/см²;

- 2-2:плотина в треугольном створе-Н= 70м; t=0,25; l=3; E_осн.=500 000кг/см²;

- 2-3:плотина в треугольном створе-Н= 70м; t=0,25; l=3; E_осн.= 50 000кг/см²;

Где Н- высота; t- коэффициент стройности ; l- коэффициент створа; Е_осн.- модуль упругости основания.

Размер автореферата не позволяет показать результаты всех вариантов, поэтому представлены лишь результаты вариантов 1-1 и 2-1.

- При угле = 0 - пульсации напряжений в контрольных элементах представлены на рис. 9. На верхних отметках в пятах (В и В) пульсация арочных и консольных напряжений существенна, арочные достигают растягивающих напряжений 15кг/см² на верховой грани и 38кг/см² на низовой грани. В ключе в нижней точке верховой грани (С) имеем арочные растягивающие напряжения до 40кг/см². Арочные напряжения как со стороны верховой грани (D), так и со стороны низовой грани (D) сжимающие, достигают -110кг/см². Максимальные напряжения в контрольных элементах представлены в таб. 3.



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

Решая задачу по варианту 1-1 получим следующие результаты:



у_A - арочное напряжение; у_K - консольное напряжение (кгс/см²)

рис. 9 Напряжения в арочной плотине по варианту 1-1 при угле б=0⁰.

а) в пятах - у основания б) в пятах - у гребня

в) в ключе - у основания г) в ключе - у гребня

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

- при подходе сейсмической волны под углом =90 (т.е. сейсмическая волна действует вдоль створа и поперек русла) картина напряжений по сравнению с первым случаем меняется мало. В элементе А растягивающие нормальные напряжения несколько ниже 5 кг/см², сжимающие консольные напряжения достигают - 20 кг/см². В верхнем поясе арочной плотины в пятах (элементе В') растягивающие арочные напряжения достигают 37кг/см² и сжимающие -110кг/см².

таб. 3 Расчетные моменты времени по пульсации напряжений в контрольных элементах (вариант 1-1; угол б=0 ).

Т(с)	Мах. Напряжение (кг/см2)	В элементе	Т(с)	Мин. Напряжение (кг/см2)	В элементе
6.08	2.39	A	5.39	-13.20	A
6.09	0.67	A'	5.26	-28.89	A'
5.78	23.44	B	5.35	-45.42	B
5.75	15.45	B'	5.32	-101.17	B'
4.59	44.44	C	5.85	-12.44	C
5.20	0.00	C'	4.58	-93.87	C'
5.82	1.41	D	5.39	-108.68	D
5.80	9.04	D'	5.37	-128.12	D'

таб. 4 Расчетные моменты времени по пульсации напряжений в контрольных элементах (вариант 1-1; угол б=45 ).

Т(с)	Мах. Напряжение (кг/см2)	В элементе	Т(с)	Мин. Напряжение (кг/см2)	В элементе
5.13	4.50	A	4.80	-18.45	A
5.14	5.84	A'	5.33	-29.56	A'
5.08	30.28	B	5.44	-51.03	B
5.03	26.08	B'	5.38	-123.20	B'
4.59	40.45	C	5.86	-10.94	C
6.54	0.00	C'	4.57	-91.37	C'
5.10	0.97	D	5.41	-95.29	D
5.07	4.71	D'	5.38	-115.52	D'

таб. 5 Расчетные моменты времени по пульсации напряжений в контрольных элементах (вариант 1-1; угол б=90 ).

Т(с)	Мах. Напряжение (кг/см2)	В элементе	Т(с)	Мин. Напряжение (кг/см2)	В элементе
5.14	7.02	A	4.82	-21.11	A
5.14	8.69	A'	5.39	-27.52	A'
5.12	40.61	B	4.80	-52.99	B
5.07	34.04	B'	5.43	-110.73	B'
4.56	25.07	C	5.20	-9.26	C
7.26	0.00	C'	4.54	-79.03	C'
5.07	0.62	D	5.59	-65.34	D
5.11	4.10	D'	4.78	-72.62	D'

таб. 6 Сопоставление максимальных напряжений в арках прямоугольных и треугольных створах - варианты 1-1 и 2-1.

вариант расчетов	максимальное напряжение (кг/см2)
No	? (градус)	арочная плотина с постоянным радиусом	в контр. Элементе	арочная плотина с постоянным углом	в контр. Элементе
1	0	32.80	D'	24.50	C
2	15	31.59	D'	23.86	C
3	30	27.61	B'	22.19	C
4	45	32.41	B'	20.25	C
5	60	34.84	B'	19.10	B'
6	75	33.79	B'	18.76	B'
7	90	30.36	B'	18.22	B'



вариант1-1: цилиндрическая плотина - Н=200м; t=0,25; l=3; E_осн.=500 000кг/см²

вариант1-2: цилиндрическая плотина - Н=200м; t=0,25; l=1; E_осн.= 50 000кг/см²

вариант2-1:плотина в треугольном створе - Н=200м; t=0,25; l=3; E_осн.=500 000кг/см²

вариант2-2:плотина в треугольном створе - Н= 70м; t=0,25; l=3; E_осн.=500 000кг/см²

вариант2-3:плотина в треугольном створе - Н= 70м; t=0,25; l=3; E_осн.= 50 000кг/см²

(Н-высота; t- коэф. стройности ; l- коэф. створа; Е_осн.- модуль упругости основания)

рис. 10 Максимальное напряжение в арочных плотинах в зависимости от направления сейсмических волн.

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

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

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

На рис. 10 показаны значения максимальных растягивающих напряжений при различных углах подхода сейсмической волны. Видно, что при цилиндрической арке максимальные напряжения при угле подхода = 60, а минимальное при = 30, но стоит уменьшить коэффициент створа до 1 и уменьшить модуль упругости основания в 10 раз и картина меняется - пик растяжения при = 75, но при 60 и 90 значения напряжений растяжения близки, а минимум достигается при = 0.

В треугольном створе напряжения достигают максимума растяжения при = 0, а минимума при = 90. При изменении модуля упругости основания в 10 раз и коэффициента створа в пределах от 1 до 3 картина не меняется: пик растяжения при = 0, а при других углах подхода значения близкие.

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

Арочная плотина Нам-Нган имеет высоту 50м, толщина у основания 8 м, поверху - 2 м; центральный угол 120.



свойства материалов для статических задач
Вид материала	E, МПа		, тс/м3
бетон плотины	30000	0,18	2,40
основание в зоне IB	4000	0,28	0,00
основание в зоне IIА	9000	0,26	0,00
основание в зоне IIB	14000	0,25	0,00
свойства материалов для динамических задач
Вид материала	E, МПа		, тс/м3
бетон плотины	30 000	0,18	2,40
основание в зоне IB	8 000	0,28	2,61
основание в зоне IIА	18 000	0,26	2,67
основание в зоне IIB	28 000	0,25	2,67

рис. 11 Пространственная сетка МКЭ и свойства материалов.

В результате получились период основного тона составляет около 0,24 с, второго и третьего - около 0,21 с, десятого - 0,15 с. На всех представленные формах собственных колебаний на гребне плотины наблюдаются максимальные горизонтальные смещения вдоль русла.

В зоне контакта с основанием растягивающие напряжение достигает 17кгс/см², и время нахождения элементов в растянутой зоне составляет около 0,2с.



у₃ на напорной граниу₁ на низовой грани

рис. 12 Напряжения плотины при сейсмическом воздействии вдоль русла.

у₃ на напорной граниу₁ на низовой грани

рис. 13 Напряжения плотины при сейсмическом воздействии вдоль створа.

Основные выводы

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

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

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

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

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

6. Не столь сильное влияние угла подхода сейсмической волны к плотине можно объяснить большим влиянием на напряжения в плотине вертикальной составляющей акселерограммы, а она не зависит от угла подхода. Горизонтальные составляющие действуют на ее фоне и потому их изменения от угла подхода не столь значительны. Видимо, в каждом конкретном случае требуется свой анализ влияния угла подхода сейсмической волны, хотя для треугольного створа во всех случаях наиболее опасен угол = 0.

7. Исследование сейсмической задачи о НДС арочной плотины на фоне её пространственного статического решения выявило целый ряд недостатков этого варианта: неустойчивость плотины в строительный период даже от собственного веса. При сейсмических воздействиях положение существенно ухудшается. Такие исследования арочной плотины дают очень большую информацию для проектирования. Исследовался также угол подхода сейсмического воздействия к оси плотины. Наихудший вариант - подход сейсмической волны вдоль русла реки.

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

Основные положения диссертации опубликованы в следующих работах автора

1) Рассказов Л. Н., Нгуен Куанг Кыонг - Собственные формы и собственные значения арочных плотин, «Вестник МГСУ», №2, 2006 - стр. 28 ч 42

2) Рассказов Л. Н., Нгуен Куанг Кыонг - Влияние направления сейсмического воздействия на напряжения в арочных плотинах, «Гидротехническое строительство», №9, 2006 - стр. 32 ч 36

Размещено на Allbest.ru

...