3. Последовательность вычислений при решении задачи n2
Все вычисления при решении задачи N2 выполняются на ЭВМ. Для выполнения расчетов составлена специальная программа, загрузочный модуль которой имеет наименование DINBAR.EXE. Работа программы организована в среде WINDOWS в диалоговом режиме.
Система исходных уравнений (4.1) - (4.3) в программе решается методом характеристик. Это позволяет вычислять приращения отметок дна в
зоне заиления верхового участка водохранилища с течением времени. Результаты расчетов выводятся на печать в табличном виде.
В качестве исходных данных для выполнения расчетов с использованием программы DINBAR задаются следующие параметры:
- ширина реки в верхнем сечении, м;
- глубина реки в верхнем сечении, м;
Q
- расход воды, м
/c;
DB - коэффициент расширения;
К - шаг по времени, сут.;
М - число моментов времени;
N - число расчетных сечений.
Программа счета
состоит из головной программы и двух
подпрограмм. Головная программа
организует ввод параметров задачи (
,
выполняет вычисление линейной координаты
для нулевого момента времени для каждого
из расчетных сечений. После этого в
головной программе вычисляются значения
глубины русла для нулевого момента
времени
в каждом сечении с учетом уклона дна
.
По величине коэффициента расширения
русла
определяется значение ширины русла
и площади живого сечения
в каждом сечении для нулевого момента
времени. Затем вычисляется угловой
наклон линий характеристик и определяются
значения глубины деформированного
русла на любой заданный момент времени.
Подпрограмма вывода организует вывод на дисплей результатов вычислений и их печать на принтере.
По результатам
вычислений строятся совмещенные профили
дна в моменты времени
и
.
Графики, изображающие результаты расчетов, оформляются в виде рис.4.2.
ЛАБОРАТОРНАЯ РАБОТА N5
Расчет деформаций дна реки методом конечных
разностей
ЦЕЛЬ РАБОТЫ: выполнить прогноз деформаций на участке реки, решая уравнение деформации методом конечных разностей.
ИСХОДНЫЕ ДАННЫЕ: расход воды, число расчетных поперечных сече- ний, шаг по времени, число расчетных моментов времени. Кроме этих данных для решения задачи задаются массив расстояний поперечных сечений от постоянного начала, массивы площадей живого сечения, ширин и средних глубин сечений.
ПОРЯДОК ВЫПОЛНЕНИЯ РАБОТЫ
1. Общие указания
Рассматривается участок реки длиной 20-30 ширин русла. Применяется одномерная модель потока. Движение воды за избранный промежуток времени считается квазиустановившимся. Влияние деформаций на положение свободной поверхности не учитывается.
Все вычисления выполняются на ЭВМ по заранее составленной программе расчета. Программа имеет сокращенное наименование DINBED и реализована в диалоговом режиме в среде WINDOWS. Уравнение деформации в программе решается методом конечных разностей по левой явной трехточечной разностной схеме. Результаты расчета приращений отметок дна реки выводятся на печать.
Исходные уравнения,
решаемые в программе, записываются в
виде 4.1 - 4.3. Берега считаются недеформируемыми
и изменение ширин русла по длине участка
задается. Уравнение деформации (4.1)
решается при заданных начальном и
граничном условиях.
Начальным условием служит распределение площадей сечений по длине участка в начальный момент времени. Верхняя граница участка назначается в сечении с пренебрежимо малыми деформациями дна.
Граничное условие
имеет вид
.
В качестве исходных данных для выполнения расчетов задаются следующие параметры.
Q - расход воды,
м
/c;
N - число поперечных сечений;
L - расстояние от постоянного начала, км;
- площадь поперечного сечения, м ;
В - ширина сечения, м;
К - шаг по времени, сут,;
М - число моментов времени.
2. Разностная запись уравнений
В двумерной области
,
вводится прямоугольная сетка координат
с переменным шагом по длине и постоянным
шагом по времени. Разностная запись
уравнения деформации имеет вид
(5.1)
При этом средняя
глубина в живом сечении и величина
приращения отметки дна для j
- го момента времени определяютя
соответственно
(5.2)
(5.3)
3. Последовательность вычислений
Программа счета состоит из головной программы и четырех подпрограмм.
Головная программа является программой-диспетчером. С помощью подпрограмм осуществляется ввод параметров задачи и массивов исходных данных и основной цикл расчетов значений площади живого сечения и приращения отметок дна за весь рассматриваемый период времени. Специальная подпрограмма формирует таблицы результатов вычислений и выводит их на дисплей, а при необходимости на печать.
Ввод исходных данных осуществляется в соответствии с командами монитора. После выполнения программы результаты расчета приращений отметок дна выводятся на экран монитора в виде таблицы.
По табличным данным строятся совмещенные профили дна для всех расчетных моментов времени и выполняется анализ полученных результатов.
СОДЕРЖАНИЕ
стр.
Введение
Лабораторная работа N1
Лабораторная работа N2
Лабораторная работа N3
Лабораторная работа N4
Лабораторная работа N5