9. Решая полученную систему линейных алгебраических уравнений, получим коэффициенты ©О, для всех узлов сетки, далее по формуле (22) получим решение дифференциального уравнения (17) с граничными условиями (19)-(20). 10. Далее пункты 6-9 алгоритма циклически повторяются для последующих шагов по времени. Вышеприведенный алгоритм позволяет неустановившегося движения воды в открытых руслах. Подготовка данных для моделирования. Основными данными для моделирования являются гидравлические и морфометрические параметры дна и потока воды: — данные дна участка, которые задаются по узлам разбитой сетки области в следующей таблице с полями { х, у, 2}, где 1 — номер узла, х — координата по оси х, у — координата по оси у, 2 — координата по оси 5. осуществить моделирование Например: П х У 2 1 234201.37 441578.18 0.10 2 234471.73 441578.18 0.20 3 234742.08 441578.18 0.30 — элементы области и способы аппроксимации задаются в следующей таблице с полями {& #, 1, 1з, 4, арг}, где К — номер элемента, 11 — первый узел элемента, 1› — второй узел элемента, 13 — третий узел элемента, 14 — четвертый узел элемента, арг — способ аппроксимации при выборе базисов элемента. Например: К і] і2 і3 і4 арг 1 41 42 12 11 лИн 2 42 43 13 12 лИн 3 43 44 14 13 лИн — — граничные условия задаются с помощью таблицы с полями {& 1, 12, 13, М, где & — номер границы, 11 — первый узел границы, 1› — второй узел границы, 13 — третий узел границы, 14 — четвертый узел границы и др. Например: И 12 13 н 15 16 300 330 270 240 210 180 Ю | - | =- 301 271 231 201 171 161 Результаты моделирования. Было промоделировано двумерное неустановившееся движение воды для следующих случаев. Пример 1. Прямоугольный участок русла с сужением с помощью стенок с двух сторон. Параметры прямоугольного русла по ширине 2700,0 м, а по длине 7600,0 м. Русло разбито на 261 четырехугольный элемент, примерная площадь каждого четырехугольника равна 68—67 тыс. кв. м. Участок русла (рис. 1 а) был разбит на прямоугольные элементы, в середине участка имеется сужение участка с помощью твердых тонких стенок. Длина стенки составляет три ширины прямоугольного элемента, ширина стенки равна одну треть длины элемента. -24-