Исходной информацией при числовом моделировании служили результаты решения системы дифференциальных уравнений, описывающих нестационарное изотермическое течение газа. К полученным в результате числового решения величинам давления добавляли центрированную адаптивную помеху с дисперсией, соответствующей погрешности измерительного прибора класса 1. Аналогично имитировались погрешности измерения расхода, которые соответствовали погрешностям измерительного прибора класса 2.
В табл. II.6 приведены массивы исходной информации и результаты моделирования.
Массив р — двумерный массив относительных величин р, 5 = 37 (37 временных слоев); массив Qo — относительных величин расхода; D—диаметр газопровода, м; Ах — протяженность элементарного__участка, х=15км; |0 — начальное значение величины £— 1/УА, для элементарного участка (выбирается произвольно); т, п — постоянные величины; B — sign{...}; Аг — начальное значение параметра А (выбирается произвольно); А/— интервал дискретизации по времени; Д£=1800 с; с — скорость
47
Таб. |
пи |
ца Ц.6 |
|||||
Ряехоп |
Давление по элементарным участкам |
||||||
Время |
Q.10-7, |
||||||
тыс. м*/сут |
1 |
2 |
3 |
4 |
|||
0 |
222 |
1,203283 |
,210403 |
1,206970 |
1,204463 |
||
1 |
222 |
1,168380 |
1 |
[,158418 |
1,165866 |
1,166906 |
|
2 |
222 |
1,125972 |
1 |
,117540 |
1,141580 |
1,148640 |
|
3 |
222 |
1,094797 |
1 |
[,094040 |
1,078330 |
1,098059 |
|
4 |
223 |
1,068330 |
1 |
,057070 |
1,056600 |
1,063440 |
|
5 |
224 |
1,007820 |
1 |
,013815 |
1,014646 |
t ,026064 |
|
6 |
226 |
0,9845051 |
0,965390 |
0,981010 |
0,969186 |
||
7 |
228 |
0,9387153 |
0,938252 |
0,915440 |
0,900970 |
||
8 |
233 |
0,896400 |
0,862295 |
0,841510 |
0,824728 |
||
9 |
233 |
0,808133 |
0,778330 |
0,737000 |
0,710309 |
Продолжение табл. II.6
Время |
Расход Q-10—?, тыс. м*/сут |
Давление |
; По элементарным |
участкам |
|
5 |
6 |
7 |
|||
0 |
222 |
1,208233 |
,196553 |
1,198063 |
|
1 |
222 |
1,164050 |
,162250 |
1,175530 |
|
2 |
222 |
1,130810 |
,132637 |
1,119620 |
|
3 |
222 |
1,099190 |
,091330 |
1,090800 |
|
4 |
223 |
1,056715 |
,043909 |
1,075590 |
|
5 |
224 |
1,017264 |
,996140 |
0,999590 |
|
6 |
226 |
0,947717 |
0,955930 |
0,939980 |
|
7 |
228 |
0,893700 |
0,884560 |
0,865727 |
|
8 |
233 |
0,810580 |
0,798950 |
0,765944 |
|
9 |
233 |
0,687550 |
0,635520 |
0,622060 |
звука в газе, с = 375 м/с; хо = 0,000084 с/м; ро = 0,0000307 с/м. Расчет 1/А приведен в табл. II.9, 11.10, 11.11. Эти таблицы отличаются выбором коэффициента ys.
Размерности и диапазоны изменения констант, входных и выходных величин алгоритма приведены ниже.
Параметр Я М At, с N D Ц |
Диапазон измерения
0—200
0—200
0—3600
0—5
0,5-1,4
0—200
ps, ps-!, кгс/см2 20 • 10«—76 • 104
Qs_i, кг/с 50—150
Т, К 220—350
Параметр
Ах, м
тп, п
/V кгс/см2
Диапазон измерения
10 000—50 000
3—5 0,1—0,3
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.