Файл: Баясанов, Д. Б. Автоматизированные системы управления трубопроводными объектами коммунального хозяйства.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