STROKA-01 +15+20+01
KX KT KRK
KX - число узлов
разностной сетки по координате;
KT - число узлов разностной сетки по времени;
KRK - число вариантов счета в режимах MR
= 0,2,3 (в остальных случаях KRK =
1). Программа считает данный пример KRK раз, требуя ввода СТРОК с данными,
характеризующими водно - физические и
гидрохимические свойства грунта;
STR-02 TEKCTОВАЯ
ИНФОРМАЦИЯ О РЕШАЕМОЙ ЗАДАЧЕ
STR-03 06 03 01 01 01 01 02 00
K KS IP MR MH MC MD MGP
K - число точек, в которых определены исходные
эпюры. Исходные эпюры могут быть первоначально определены не во всех узлах
сетки, а аппроксимируются программой по законам ступенчатой или ломаной
функций,K<99;
KS - число слоев с различными свойствами в рассматриваемой области
KS £
9;
IP - признак
эпюры. При IP = 0 исходные эпюры аппроксимируются на сетке ступенчатой
функцией, а при IР = 1 - ломаной.
MR - переменная, задающая
режим счета:
MR = 0. Расчет
перераспределения компонент по глубине почвогрунта на один заданный период времени ТО . В режиме MR =
О могут подобраться параметры уравнений путем
сопоставления теоретических и экспериментальных эпюp
влажности или солей. При
этом для ряда значений параметров уравнений рассчитываются теоретические эпюры
(исходные эпюры остаются прежними для каждых новых вводимых данных).
MR = 1. Расчет режима. В ЭВМ вводится
последовательность временных интервалов ТОi
с заданным законом изменения краевых условий. На каждый новый расчетный период
времени за исходные принимаются распределения компонент, полученные на конец
предыдущего периода. На печать выдаются эпюры:
MR = 2. Расчет выходных
кривых с шагом DТ по времени до времени DP . Служит
для подбора параметров по выходам кривым. На печать
выдается изменение компонент на любой заданной глубине х(kf) с заданным
шагом. Расчет может повторяться для других
параметров.
MR = 3. Определение
времени и нормы промывки.
Расчет выполняется при МС ¹
0 и задании исходящих
потоков влаги. Счет продолжается до тех пор, пока минерализация
порового раствора в заданной
точке X(KF) не станет меньше допустимой СKF < DP.
На печать выдаются значения компонент в этой точке с шагом DТ по времени. После
выполнения условия СKF < DP
происходит
возврат назад на один шаг DТ, DТ уменьшается в 10 раз и расчет
повторяется.
MR = 4. Расчет режима. Аналогичен случаю MR = 1, но на печать выдаются значения компонент только на заданной глубинe X(KF).
KF – фиксированный узел
сетки.
MH, MC – переменные, задающие компоненты, участвующие в расчете,
могут принимать следующие значения:
MH = 1 – влагоперенос
рассчитывается;
MH = 0 – влагоперенос не
рассчитывается (задается скорость фильтрации постоянной по профилю
почвогрунта);
MC = 0 – солеперенос не
рассчитывается;
MC = 1 - солеперенос
рассчитывается по уравнению (3);
MC = 2 - солеперенос
рассчитывается по уравнениям (4) –
(5);
MD = 1 имеется
горизонтальный дренаж в однородном пласте подстилаемом водоупором (строение
зоны аэрации выше глубины заложения дренажа, можем быть неоднородным);
MD
= 2 - имеется горизонтальный дренаж в слоистой толще;
MD
= 0 – параметры дренажа не
вводятся;
MGP = 0, если
исходные концентрации вводятся в г/л;
MGP = 1, если
исходные концентрации вводятся в %
STR-04 000.000 000.500 001.000 001.500 002.500
003.500
HJ(J)
HJ(J) - координаты, в которых
определены исходные эпюры, м; HJ(1) = 0, НJ(К)= L -задает
нижнюю границу области расчета;
STR-05 000.000 000.450 001.400
HS PO G
000.300 000.450 001.400
001.000 000.440 001.490
HS(J), J=1,..,KS - верхняя
граница каждого литологического слоя, м;
НS(1) = 0 (в области расчета L);
PO(J) , J=1,..,KS - пористость (m) каждого литологического
слоя, доли 1;
G(J), J=1,..,KS - плотность ( g )
каждого литологического слоя
(объемный вес грунта), кг/м3;
STR-06 002015
KDX/KDP ИЛИ
DT/DP/KF
Управление печатью
результатов имеет две различные формы: при MR = 0 или 1 вводятся
KDX – число узлов
координатной сетки, через которые производится вывод эпюр на печать (если KDX = 1, вывод производится в каждом узле);
KDP – последний узел,
в котором производится вывод на печать (KDP £ KX).
При MR = 2,3,4 вводятся:
DT – интервал
времени, через который производится вывод на печать результатов счета, сут.
(при MR = 4: DT = DP = 0);
DP (при MR = 2 )- конечное
время, до которого производится расчет; DP (при
MR = 3) – допустимая
концентрация, до которой должно произойти опреснение раствора в слое до
глубины X(KF)
при этом г/л при MGP = 0, а % при MGP = 1;
KF – узел сетки, в котором производится вывод на печать.
STR-07 01 04 01
MHW MPS MKW
Исходные данные для
уравнения влагопереноса (вводятся при MH = 1)
MHW – переменная, определяющая вводимую эпюру, может принимать 2
значения: 1 – в качестве исходных вводятся значения потенциала P, м; 2 – в качестве исходных вводятся значения влажности W;
MPS – переменная, задающая зависимость P ( W ): MPS = 2 [ ]; MPS = 3 [ ]; MPS = 4 [ ]
MKW - переменная, задающая зависимость K ( W ): MKW = 1 [ ]; MKW = 2 [ ];
STR-08 -002.000-001.500-001.000-000.500 000.500
001.500
PJ
Имеет 2 формы в
зависимости от значения MHW:
P(J), J=1,..,K – значения
капиллярного потенциала P в тпчках HJ(J), м (вводится при MHW = 1);
W(J), J=1,..,K – значения
объемной влажности W в точках HJ(J), м (доли 1 вводится при MHW = 2, в зоне полного насыщения
полагать W = WP);
STR-09 650.000 000.500 003.000 010.000 000.650
160.000
B DDR
HDR TB KB PH
При MD
= 1 вводятся:
B - междренное расстояние, м;
DDR - диаметр дрены,
м;
HDR - глубина заложения дренажа, м;
TB - мощность
покровного слоя, м;
KB - коэффициент фильтрации покровного слоя, м/сут (дренаж
заложен в покровном слое);
PH – проводимость нижележащих слоев. При
MD = 1 PH = 0, При MD = 2, PH = åKiTi, м2/сут.; Ki, Ti - коэффициенты фильтрации и мощности
слоев ниже покровного.
STR-10 001.500 000.900 000.000 000.050 000.440
000.190 000.120
KO HK WMM WO PV PPV WZ
001.260 000.900 000.000 000.072 000.440 000.190 000.120
000.650 001.350 000.000 000.072
000.430 000.240 000.120
Для каждого
литологического слоя в области расчета вводятся:
KO - коэффициент
фильтрации, м/сут,
HK - hK или h*K в формулах (20)-(24), м;
WMM - W* в формулах (20)-(21), м3/м3;
WO - W0 , WM в формулах (21)-(24), м3/м3;
PV - полная
влагоемкость, м3/м3;
PPV -
предельно-полевая влагоемкость, м3/м3;
WZ - влажность
завядания, м3/м3;
STR-11 -500.000-500.000-500.000
PM(J)
PM(J): J=1,..,KS. Имеет две
формы: PM=n в зависимости P(W); (n<1) при MPS = 2;
PM = -500 м при MPS = 3 или 4.
STR-12 003.500 003.500 004.000
NKW(J)
NKW(J), J=1,..,KS – показатель степени nkw в зависимости K(W).
STR-13 000.458 000.707 000.931 001.081 005.000 005
000
C(J), J=1,..,K – минерализация порового раствора («сквозных» пор при MC =2) в точках HJ(J), г/л при MGP = 0, %
при MGP = 1
STR-14 N(J), J=1,..,K – минерализация порового раствора «тупиковых» пор в точках HJ(J),г/л при MGP = 0, % при MGP = 1
STR-15 001.000 000.830 000.100
CA CLD BE
001.000 000.830 000.100
001.000 000.830 000.100
Вводится на каждый
литологический слой (KS штук) -
CA – доля «сквозных» пор ( k ) в общей пористости, СА = 1 при MH =1
CLD – параметр
гидродисперсии ( l*
), м;
BE – параметр обмена
( b*
), 1/сут (BE =0 при MC = 1)
STR-16 031.000
-0.00150 0.00020 0 2 3 0 0 0 0
0
TO H1 H2 KE/G1/G2/LR/LE/CP/BC2/C2
021.000 -0.00260 0.000201 2 3 1 0 0 0 0 0
Текущие
исходные данные:
TO
H1
H2
KE
G1
G2
LR
LE
CP
BC2
C2
STR-17 000.00210 000.4780 001.4830 001.4830
000.750
EK SL
X1 XL HKOR
EK
SL
X1
XL
HKOR
STR-18 000.240 000.140 001.000 000.740
001.000-000.200 000.200
WMAX WMIN CMAX CMIN
XMAX H1 CP
WMAX
WMIN
CMAX
CMIN
XMAX
H1
CP
STR-19 EO
W1 W2 HE ИЛИ EO
ALFA
EO
W1
W2
HE
EO
ALFA
0
0