Файл: Баясанов, Д. Б. Автоматизированные системы управления трубопроводными объектами коммунального хозяйства.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 17.10.2024
Просмотров: 135
Скачиваний: 0
= |
t ( f / i ) i , A + i + |
( y i ) b f e l / 2 - [ ( t / i ) o , f t + x + |
(t/i)o, ftl/2 |
„ _ |
|
- |
— , |
(3.97) |
получим конечно-разностное уравнение на границе х = О,
которое можно привести к виду:
|
(^i)o, a+ i = ^ i (^ i)i , ft+i + |
v i, |
(3.98) |
||
где щ = 1; Vi = |
( U j ) о,л,— |
|
|
|
|
Используя 'представление первой производной в гра |
|||||
ничном условии (3.94) в виде: |
|
|
|
||
|
д и 2 ( х , |
<)[ |
= |
|
|
1(^2)„,, k + ( U |
дх |
U=R |
|
|
|
2 ) n , . k V 2 - W |
i n r l |
. k + l |
+ m « , - l . k \ V |
/0 nm |
получим конечно-разностное уравнение на границе стержня
х = R, которое можно привести к виду:
W d k , , k + \ = K z ( U 2 ) n , - 1 , ft-ь 1 -bv2, |
(3.100) |
гм « " а л Т ь ’
— {U 2) ni_ [ _k\ + 2 a h ( ^ f e + i) .
Аналогично получим конечно-разностное уравнение в точ
ке соприкосновения двух сред х = г, соответствующее гра
ничным условиям (3.92) и (3.93)
(^1)л, , ft+1 = |
(^2)о, fe + 1= |
х3 (^l)rtj — 1 , ft+l + |
|
|
+ X4(£/2)I f * + i + V „ |
(3.101) |
|
где и3 - |
%1 + %2 : |
|
+ Лз ; |
где |
|
|
|
{^1 [(Ul)ni — \, |
*] +^2 |
[(^2)l, ft— (^2)0, ft]}. |
Использованная выше конечно-разностная аппроксимация
на границах дает второй порядок аппроксимации получен ной схемы. Системы уравнений (3.95) и (3.96) решают мето
дом прогонки, о котором шла речь выше, а систему (3.95) решаем в виде:
(Ui)i, ft+i = «i+ i(f/i)i+ i, ft+x + P m , (г = 0 , 1 ,2 ,..., tii — !)• (З.Ю2)
С помощью (3.95) для прогоночных коэффициентов
получаем следующие рекуррентные соотношения:
В |
о |
Л6 i + F i |
« - 1 ' 2 ...... |
<злоз> |
1S8
Используя граничное соотношение (3.98), получим
|
«1 = xy, P!=Vi. |
(3.104) |
|
Систему уравнений (3.96) решаем в виде: |
|
||
(^2)j, ft+i= а / + 1(£/2)^+1, ft+i + Р /+ 1. |
|
||
|
(/ = 0, 1 , |
2....... л2 - 1 ) . |
(3.105) |
Из выражений .(3.105) и (3.96) для коэффициентов а /, |
(Д по |
||
лучаем соотношения: |
|
|
|
«У+1 ; |
С '—Л 'а )/> Р/ + 11' С ' — А ' а ) ’ |
|
|
|
(/ = 1 . |
2 ,..., n2— 1). |
(3.106) |
Величины а[, |3{ определяют из соотношений . |
|
||
—1 , ft- 1 = |
а Я1 (^l)n,. ft+l +Pn, |
(3.107) |
|
(результат прогонки через |
область 0 ^ х ^ г); |
|
|
(^2)0 , k+ 1 — |
(^1)щ — i , ft-ц i + к 4 (^2) 1 1 * 4 -1 + v3 |
(3.108) |
(соотношение на внутренней границе соприкосновения двух
сред). Исключая из соотношения (3.107) и (3.108) неизвест
ные величины {U^)nt, k+\ |
и |
(U i)n,, *4-1 |
с учетом равенст |
||||||
ва |
(U^m, ь+ 1 = |
(Д2)0, ft+i> |
являющегося следствием |
усло |
|||||
вия |
(3.92), |
получим: |
|
|
|
|
|
||
|
|
|
|
*4 |
|
|
ЯзРл,+Уз |
(3.109) |
|
|
(^2)0, й+ 1 = ! —из a |
(^2 )1 , ft+i + |
I 4 -K3 ani |
||||||
|
|
||||||||
Отсюда |
|
|
|
|
|
|
|
|
|
|
|
|
, |
к4 |
|
|
Ч-'Мз |
(3.110) |
|
|
|
ai = |
-----------> Pi = - ----------- |
||||||
|
|
|
|
! — x 3 a„ |
1— « з % |
|
|||
Для начала |
обратной прогонки необходимо знать величину |
||||||||
(U 2)„2, ft+1, |
которую определяют |
из |
соотношений: |
|
|||||
|
|
m |
a t - |
i M k + l = |
« |
,n. m a „ |
k + |
l + K t \ |
(3.111) |
|
|
|
|||||||
|
|
№ |
) „ „ |
ft4 - 1 — х 2 (^г)л2 —1 , ft4-i + У 2 - |
(3.112) |
||||
Из |
уравнений |
(3.111) и (3.112) получим |
|
||||||
|
|
|
.... |
|
И2 Р/!,+У 2 |
|
(3.113) |
||
|
|
|
(^ „ ..ft- H - 1 _ И2а'г |
• |
|||||
|
|
|
|
189
Согласно рассмотренному выше алгоритму можно сформули ровать следующие основные этапы вычислительного про
цесса:
1) |
по формуле |
(3.93),Ч используя |
соотношения (3.94) |
||||||||
в качестве |
начальных, |
у определяем |
а г+1, |
|Зг+1 |
(i = |
||||||
= 1 . 2 , ......... |
п х— 1 ); |
|
|
используя |
соотношение (3.110) |
||||||
2) |
из уравнения (3.96), |
||||||||||
в качестве |
начального соотношения, определяем a/_pi |
(5)+i |
|||||||||
(i = |
1 , 2 , ...... , п 2 — 1); |
|
|
|
|
|
|
||||
3) |
по формуле |
(3.105), |
используя |
(3.113) |
в |
качестве |
|||||
начального соотношения,^определяем (U 2)^ h+1 (/ |
= |
п 2 — 1 , |
|||||||||
п2— |
2, ... |
2, |
1, |
0); |
|
|
|
|
|
|
|
4) используя {U x)nt, fe+1 |
= |
(U 2)0i ft+1 |
в качестве |
началь |
|||||||
ного соотношения, |
из |
уравнения (3.102) находим (U x)it h+1 |
|||||||||
(i = |
п х — |
1 , |
— 2 , |
..., |
2 , |
1 , 0 ). |
|
|
|
|
|
|
§ 5. ПРОБЛЕМА ФУНКЦИОНАЛА В ЗАДАЧЕ |
|
|
|
|||||||
|
АВТОМАТИЗИРОВАННОГО УПРАВЛЕНИЯ |
|
|
|
|||||||
|
Для |
структуры |
типа |
автоматизированной |
системы |
управления трубопроводными предприятиями коммуналь
ных хозяйств характерны многомерность, многосвязанность,
взаимообусловленность, сочетание подструктур автомати
зированного и автоматического действий (как во всякой
человеко-машинной системе), наличие неформализованных
подструктур и элементов с остаточной неопределенностью, иерархический принцип организации подструктур и т. п.
Будучи многомерными, анализируемая система АСУ и
процессы в ней могут быть описаны при помощи векторных
временных числовых последовательностей £1( £2, |
..... , l h или |
|
же векторных функций времени £ (t) — |
(t), |
g 2 (t), ..., |
•••,} и операторно-матричных соотношений между ними.
Это значит, что если на входе i-й подструктуры по I-му виду
работ в /-й структуре 6 -й подсистемы описываемой системы
АСУ (в дальнейшем такой элемент |
будем |
называть [ i — |
|
— l — j — & 1-элемент) действует конечной |
длины векторная |
||
функция времени 1 (t) = {gx (i), |
(t), |
..., |
lN (t)} и в ре |
зультате функционирования элемента возникает вектор
выходных координат ф (t) = {фх (t), ф2 (i), ...... , |
ф^ (t)} |
|
длины |
М, совпадающей и отличающейся от размерности |
|
N, то |
оператор-матрица [i — I — / — 61-элемента |
будет: |
|
q ( t ) = A t ¥ (t), |
(3.114) |
190
где A t= \ A t t, A u , ... , A , J ,
Ы * ) = А и 1т(1);
4>* |
( t ) = A tt | r (0; |
(3.115) |
|
|
yM(t)=Ath l T (/).
Конкретный вид оператора ф (t) отражает свойства многосвязанности и взаимообусловленности элементов системы,
а также наличие в ее составе независимых, иначе, авто
номных структур, описываемых блочными операторами (рис. 28). Так, например, рассматривая комплекс вычисли
тельных работ, обеспечивающих формирование многоднев
ных прогнозов в подсистеме оперативного управления АСУ
трубопроводным хозяйством, вектор входных координат
подсистемы формируется на основе анализируемых отклоне
ний от плана тех или иных процессов и учитывает действие
на эти процессы многих факторов, а оператор-матрица А является здесь набором операторов экстраполяции. Век
тор на выходе U — I — j — k ]-элемента имеет составля
ющими прогнозируемые значения отклонений от плана про
цесса на определяемом интервале времени в пределах теку
щего и планового периодов.
Наличие в описываемой системе АСУ неформализован
ных подструктур обусловливает постановку задач иденти
фикации, иначе — математического моделирования та ких подструктур и их элементов. При решении этого ком
плекса задач не только определяются или уточняются па
раметры системы, но и моделируются на основе принятия
некоторых гипотез относительно применяемого класса мо
делей самих подструктур системы АСУ. На самом деле, ана лиз объектов городских коммунальных хозяйств в реаль ных условиях достаточно часто сталкивается с ситуацией, когда требуется определить динамические характеристики конкретного производственного объекта, составить его
математическую модель, а структура взаимосвязей этого
объекта неизвестна или труднодоступна. О ней лишь можно делать в большей или меньшей степени достоверные пред
положения. Доступными для наблюдений и регистрации
данных являются некоторые отдельные точки структуры, например входы и выходы скалярные или векторные. Тогда
требуется идентифицировать изучаемую производственную
191