Название: Алгоритм определения динамических характеристик гидроупругих систем для управления гидросооружениями
Вид работы: реферат
Рубрика: Информатика и программирование
Размер файла: 93.68 Kb
Скачать файл: referat.me-133789.docx
Краткое описание работы: Рассматриваются методика и алгоритм решения задачи гидроупругости для грунтовых плотин, взаимодействующих с водной средой.
Алгоритм определения динамических характеристик гидроупругих систем для управления гидросооружениями
Маткаримов П.Ж.
Рассматриваются методика и алгоритм решения задачи гидроупругости для грунтовых плотин, взаимодействующих с водной средой.
Рассматриваемая проблема представляется актуальной при проектировании гидросооружений, так как включает в себя совершенствование модели сооружения, модели взаимодействия и алгоритмизацию поиска собственных значений системы “плотина-водохранилище”. На сегодняшний день методы исследования задач гидроупругости, т.е. совместных колебаний конструкции и жидкости весьма разнообразны. Применение тех или иных методов решения рассматриваемых задач диктуется многими обстоятельствами: характером задачи, целью исследования, принятой схематизацией явления, требуемой точностью, возможностью вычислительных средств и др.
Изложим постановку и методику решения задач гидроупругости, связанных с определением динамических характеристик упругих грунтовых плотин, взаимодействующих с полубесконечным слоем жидкости. Динамические характеристики (собственные частоты, формы колебаний) являются основными регламентирующими характеристиками (паспортом) сооружений, позволяющими заранее судить о его динамических свойствах и его поведении при различных воздействиях
Опираясь на современные достижения науки в этой области и основываясь на имеющихся материалах по влиянию жидкости на напряженно-деформированное состояние гидротехнических сооружений при динамических и сейсмических воздействиях, при формулировании задачи жидкость считаем идеальной и несжимаемой, волнообразование на свободной поверхности не учитываем. Тогда потенциал скорости движения жидкости должен удовлетворять уравнению Лапласа
(1)
и граничным условиям [1, 2]:
;
и
,
; (2)
на напорной грани (где
), скорости частиц жидкости и точек грани плотины по направленнию нормали n одинаковы
(3)
Тогда выражение для потенциала скоростей, удовлетворяющего уравнению (1) и условиями (2)-(3), будет иметь вид [2]
,
,
Далее рассматривается динамическая задача гидроупругости для грунтовых плотин. При этом для постановки задачи используется вариационное уравнение Лагранжа, основанное на принципе Даламбера:
(4)
и кинематическое граничное условие в основании:
,
. (5)
Здесь ,
,
- соответственно, компоненты вектора перемещений, тензоров деформаций и напряжений;
,
- изохронные вариации перемещений и деформаций;
-плотность материала элементов рассматриваемой системы;
- гидродинамическое давление воды.
Физические свойства тела описываются соотношениями между напряжениями и деформациями вида
(6)
Величины и
являются константами Ламе (индекс n относится к телу, с соответствующими механическими характеристиками).
В соотношении Коши учитываются только линейные члены
, i, j=1,2,3 (7)
Все задачи, поставленные в данной работе, решаются на базе метода конечных элементов (МКЭ). В частном случае, когда рассматриваются гармонические колебания полное смещение и потенциал скорости
можно представить в виде
,
, (8)
где - упругие перемещения стенки плотины, зависящие только от координаты
.
Собственные колебания грунтовой плотины с учетом водной среды водохранилища представляют собой упорядочное движение грунтовой плотины, протекающее при отсутствии внешних воздействий. Решение проблемы заключается в следующем: ищется нетривиальное решение уравнения (4) при однородных кинематических условиях в виде (5).
Постановка (8) в (4) сводит данную задачу к действительной вариационной задаче о собственных значениях в виде
(9)
,
где - амплитуда напряжений,
- искомая собственная частота и форма колебаний плотины с учетом водной среды,
- гидродинамическое давление воды на стенку плотины, которое имеет вид
(10)
При этом для решения вариационной задачи (9) используется закон Гука, геометрические соотношения Коши (7) и стандартная процедура МКЭ с использованием треугольного конечного элемента с линейной аппроксимацией поля перемещений внутри элемента.
При этом для i-го узла n-го конечного элемента первые вариации работ упругих, инерционных сил и гидродинамического давления относительно ui будут иметь вид
, (11)
,
Здесь l, m - константы Ламе, r, r0 - плотность материала сооружений и воды.
После интегрирования этих выражений получим строки матриц жесткости, массы и присоединённой массы воды, соответствующие перемещениям ui для n-го элемента. Если эти операции повторить для узлов j и k, а также для перемещения v, то мы получим матрицы жесткости и массы порядка 6´6 для n-го конечного элемента и соответствующие матрицы масс от воды к узлам конечного элемента (если эти узлы соприкасаются с водой). Таким образом, сформировав для каждого конечного элемента свою матрицу [K], [M], [Mв ] и объединив их, получим алгебраическую задачу на собственные значения для рассматриваемого сооружения с учетом взаимодействия с водой:
, (12)
где [K]-матрица жесткости плотины, [Mc ]=[M]+[Mв ] - суммарная матрица массы плотины и массы воды, w и {u} - искомые собственная частота и собственный вектор плотины, взаимодействующей с водой.
Решая уравнение (12) методом Мюллера и квадратного корня определим собственную частоту и форму колебания гидроупругой системы “грунтовая плотина с водой”. Такая постановка и методика решения задачи о собственных колебанях гидроупругих систем позволяют получить, в отличие от известных работ, более реальную картину динамики системы.
С помощью изложенного выше алгоритма исследованы собственные колебания грунтовых плотин, взаимодействующих с водной соедой.
Анализ полученных результатов позволяет сделать вывод о том, что учет взаимодействия воды с плотиной приводит к уменьшению частот собственных колебаний и для первых восьми частот составляет от 2 до 10% в зависимости от вида плотины. Эти отклонения зависят от многих факторов: от высоты плотины, характера форм колебаний сооружения, а также от значения коэффициента откоса, контактирующего с водой. При этом самым значимым фактором является крутизна верхового откоса, определяемая коэффициентом откоса - с увеличением коэффициента откоса соответственно уменьшается влияние воды на динамические характеристики плотины. Для вертикальной стенки значение гидродинамического давления воды является определяющим.
Список литературы
1. Кульмач П.П. Гидродинамика гидротехнических сооружений (Основные плоские задачи). -М.: Изд-во АН СССР. –1963. 190 с.
2. Шульман С.Г. Расчеты сейсмостойкости гидросооружений с учетом влияния водной среды. М.: Энергия, 1976. 336 с.
Похожие работы
-
Особенности реализации машинно-ориентированных алгоритмов расчета частотных характеристик канала воздействия
Рассмотрены проблемы формализованного анализа динамики сложных технологических объектов на базе топологических моделей. Приведены результаты машинной реализации алгоритмов расчета частотных характеристик.
-
Методика создания программ
Разработка большой программы отличается от разработки малой в двух основных аспектах. Текст программы для решения большой проблемы занимает много больше места, чем одна страница.
-
Задача о Ханойских башнях
Анализ и характеристика рекурсивного алгоритма решения задачи о Ханойских башнях, а также порядок его временной сложности в соответствии с пошаговым алгоритмом. Методика и особенности разработки программы, печатающей последовательность действий, ее текст.
-
Определение динамических характеристик системы
Понятие идентификацией систем. Взаимная корреляционная функция между входом и выходом системы. Определение импульсной переходной функции объекта по статистическим характеристикам, а также импульсной переходной функции по детерминированным характеристикам.
-
Алгоритм «рамо»
Имея представление о механизме функционирования минимальной ТС, можно предложить несложный рекурсивный алгоритм «Машина Открытий» (РАМО). Напомним, что суть создания ТС заключается в организации двунаправленного энергетического канала (ЭК).
-
"Последняя миля" — оптические решения
В последние несколько лет все большую популярность приобретают решения по организации "последней мили", где средой передачи информации является волоконно-оптический кабель.
-
Алгоритмические языки: использование множеств
Изучение способов описания и использования множеств, разработка алгоритма и составление программы для решения задачи. Нахождение в последовательности целых чисел таких, которые встречаются в ней ровно два раза. Набор программы, ее отладка и тестирование.
-
Алгоритм формирования ключей в процессе функционирования DES
Процесс и основные этапы реализации алгоритма формирования ключей в процессе функционирования DES с помощью языка программирования C++. Особенности работы программного алгоритма и его пошаговая реализация. Листинг получившейся программы, пример ее работы.
-
Методика и алгоритмы контроля работоспособности и диагностики сейсмометрических каналов
Приводятся результаты исследований и разработок на системе автоматического телеконтроля работоспособности регистрационных комплексов в системах сейсмомониторинга на крупных особо важных сооружениях.
-
Разработка программ с использованием динамической памяти
Поиск источников ориентированного графа. Использование динамических структур при работе с графами. Способы представления графов, операции над ними, описание программной реализации. Процедуры и функции языка. Функции работы с динамической памятью, графами.