Электронная библиотека Веда
Цели библиотеки
Скачать бесплатно
Доставка литературы
Доставка диссертаций
Размещение литературы
Контактные данные
Я ищу:
Библиотечный каталог российских и украинских диссертаций

Вы находитесь:
Диссертационные работы России
Технические науки
Математическое моделирование, численные методы и комплексы программ

Диссертационная работа:

Гребенников Дмитрий Юрьевич. Математическое моделирование и исследование аналитическими и численными методами процессов распространения нелинейных волн в трубопроводах : Дис. ... канд. физ.-мат. наук : 05.13.18 Воронеж, 2006 141 с. РГБ ОД, 61:06-1/948

смотреть содержание
смотреть введение
Содержание к работе:

Введение ^ 5

Глава L Подходы к решению задач о нестационарном взаимодействии
жидкости и конструкции 13

1 Л.Различные задачи гидроупругости 13
1.2.Общая постановка совместной задачи гидродинамики и теории

упругости 19

1.3.Моделирование волновых процессов в сжимаемой жидкости 27

1.4.Численное моделирование задач гидроупругости 31
1А1. Применение разностных схем для решения задач

гидродинамики 34

1.4.2, Метод Галеркина и методы взвешенных невязок 36

1.5.Выводы к главе 1 39

Глава 2. Аналитические методы решения нестационарных задач

распространения волн давления в каналах с упругими стенками 40

  1. Одномерная постановка задачи распространения волн в канале 40

  2. Лучевой метод разложения решения волновой задачи 43

  3. Исследование структуры переднего фронта волны давления 48

  4. Основные закономерности распространения волн в трубопроводах 52

  5. Распространение волн в трубопроводах с тонкими стенками из пористого материала 58

  1. Одномерная квазистационарная постановка задачи 59

  2. Нестационарная модель течения жидкости в трубопроводе с тонкими пористыми стенками 64

2.53. Анализ переходной зоны на фронте волн сильного разрыва

с учетом поперечного расхода 69

2.6. Выводы к главе 2 73

Глава 3- Энергетический подход к исследованию напряженно-деформированного соЬтояния стенок круглой трубы из повреждаемого

материала 75

3.1. Математическая модель деформирования повреждаемого
материала 75

  1. Энергетическая модель образования микродефектов 76

  2. Плоскость реологических состояний 78

  3. Изменение параметров материала 80

3.2. Задача о деформировании стенок круглой трубы 82

  1. Решение упругой задачи Ламе 82

  2. Постановка задачи с учетом изменения параметров материала стенки трубы 85

3.3. Решение задачи о начале активного процесса накопления

повреждаемости на внутренней поверхности стенки трубы 88
і
3-4. Решение нелийейной задачи о накоплении повреждаемости методом

последовательных приближений 92

3.4.1 - Нулевое приближение решения 92

3.4.2. Первое приближение решения 95

3.5. Выводы к главе 3 98

Глава 4, Решение неЬтационарных задач с подвижной границей с помощью

вычислительных программных комплексов 99

4.1. Возможности моделирования стандартных программных
комплексов .' 99
4.1Л. Уравнения нестационарного движения жидкости 100
4.1.2. Методы дискретизации по пространству и времени 103

  1. Расширение возможностей моделирования с помощью подпрограмм 110

  2. Численные расчеты с учетом деформаций стенок 113

4,4- Влияние параметров дискретизации на представление волнового

фронта 119
4,5. Численное моделирование нелинейных особенностей в задачах

распространения волн 124

4.6- Выводы к главе 4 128

Заключение 129

Библиографический список 131

Введение к работе:

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

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

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

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

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

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

Достижение указанной цели включало решение следующих задач:

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

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

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

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

Объекты и методы исследования. Исследовался процесс распространения

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

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

і ч

Научная новизна полученных в работе результатов определяется тем, что в ней:

  1. Построена математическая модель нестационарного потока в гидролинии с распределенными по ее длине стоками или источниками массы.

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

  3. Получен критерий прочности трубопровода под действием динамического давления за водной гидроудара с учетом повреждаемости материала стенок.

  4. Разработан итерационный метод решения задачи Ламе о напряженном состоянии стенок трубы под действием внутреннего давления с учетом накопления микродефектов в структуре материала.

  5. Гидродинамический программный комплекс дополнен подпрограммами, учитывающими линейную сжимаемость жидкости и деформирование границы расчетной области.

На защиту выносятся следующие основные положения:

  1. Обобщенная математическая модель нестационарного потока в гидролинии с распределенными стоками или источниками массы.

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

  1. Критерий прочности трубопровода под действием динамического давления за волной гидроудара с учетом повреждаемости материала стенок.

  2. Итерационный метод решения задачи Ламе о напряженном состоянии стенок трубы под действием внутреннего давления с учетом накопления микродефектов в структуре материала.

  3. Метод решения совместной гидроупругой задачи с помощью гидродинамического программного комплекса дополненного подпрограммами; реализующими линейную сжимаемость жидкости и деформирование "стенки канала.

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

нелинейных волновых процессов в гидролинии, учитывающее её

конструктивные особенности и возможные дефекты стенок дает аналитический

инструмент, позволяющий выявить закономерности распространения волн в

трубопроводах различной структуры, что позволяет проанализировать влияние

отдельных факторов на процесс распространения волны при различных

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

учетом накопления Люврежденности позволяет определять степень износа

і материала и определить характер поведения стенок с учетом меняющихся

упругих свойств материала.

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

Личный вклад автора, Определение направлений исследований, постановка задачи,> получение уточненных моделей и построение

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

Апробация работы. Результаты диссертационной работы докладывались и обсуждались на (четвертой и пятой всероссийской научно-технической конференции «Прикладные задачи механики и тепломассобмена в авиастроении» (Воронеж 2003, 2004 г. г.); конференции «Современные проблемы механики и прикладной математики» (Воронеж 2004 г.); научных семинарах Воронежского государственного университета, Воронежского технического университета, Воронежской государственной технологической академии в 2003-2006 г.г.

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

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

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

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

подвижной и неподвижной жидкости.

!

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

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

превышение давления внутри гидролинии. Интенсивность волн разрывов так

же зависит от наличия стоков или источников. В случае линейно

распределенного стока вдоль трубопровода скорость распространения волны

замедляется при движении вниз по потоку и тем самым время существования

волны на определенном участке трубопровода увеличивается, и интенсивность

волны уменьшается сильнее по сравнением со случаем ее распространения в

трубопроводе без распределенных стоков. В случае распределенных

источников (втока жидкости) скорость волны гидроудара вниз по потоку

і возрастает. Время существования волны на отрезке длины L уменьшается, и

интенсивность волны падает меньше чем для случая ее движения вниз по

потоку с распределенными стоками.

Рассмотрена динамическая теория, описывающая нестационарное

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

Полученные результаты показывают, что в нулевом приближении разложения в

ряд по толщине тонкой стенки, решение системы уравнений соответствует

I і

решению в квазистационарном приближении перемещений стенки трубопровода.

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

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

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

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

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

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

уделено оценке погрешностей дискретной модели и численного метода. Для

численного моделирования явления гидроупругости разработана подпрограмма

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

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

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

I I

Подобные работы
Швецов Александр Степанович
Исследование процессов токообразования в никель - металлогидридной системе методом математического моделирования
Зарипов Дамир Мунзирович
Математическое моделирование динамики трубопровода под действием волн давления в транспортируемой жидкости
Хачумов Максим Александрович
Разработка и исследование методов моделирования Интернет портала
Дехканбаев Дмитрий Саттаркулович
Исследование методов моделирования и разработка программного обеспечения для изучения фрактальных свойств системы N тел
Буланов Сергей Георгиевич
Разработка и исследование методов программного моделирования устойчивости систем линейных дифференциальных уравнений на основе матричных мультипликативных преобразований разностных схем
Давудпур Марьям
Разработка и исследование декомпозиционного метода моделирования дискретных структур
Резников Владимир Борисович
Разработка и исследование метода построения программного комплекса моделирования для распределенных систем с многоуровневым представлением сложных объектов
Михеев Сергей Александрович
Математическое моделирование и символьно-численные методы исследования гравитирующей быстровращающейся сверхплотной конфигурации в постньютоновском приближении
Минаков Сергей Владимирович
Структурное моделирование в классе задач о назначении и исследование генетического метода решения
Мамаева Саргылана Николаевна
Исследование характеристик электронно-оптических систем с полевыми катодами методами математического моделирования

© Научная электронная библиотека «Веда», 2003-2013.
info@lib.ua-ru.net