Фафурин В.А., Терюшов И.Н. Автоматизация технологических процессов и производств

COLOR 10

center 6, " Переходная функция модели записывается в виде "

COLOR 15

PRINT " K*exp(-s*TZ) "

PRINT " W(s)= ________________________________________ "

PRINT " (T*s+1)^n1*(A1*T*s+1)^n2*(A2*T*s+1)^n3 "

40 REM

COLOR 10

center 14, " Введите следующие значения "

COLOR 7

PRINT "1. К - коэфициент передачи: "; : COLOR 11: INPUT "", K

COLOR 7

PRINT "2. Постоянная времени Т (мин): "; : COLOR 11: INPUT "", T

COLOR 7

PRINT "3. Kоэфициент A1: "; : COLOR 11: INPUT "", A1

COLOR 7

PRINT "4. Kоэфициент A2: "; : COLOR 11: INPUT "", A2

COLOR 7

PRINT "5. Время запаздывания tau (мин): "; : COLOR 11: INPUT "", TZ

COLOR 15

PRINT "Значения введены корректно (Y/N) "; : COLOR 14: INPUT ""; da$

IF (da$ <> "y") AND (da$ <> "Y") THEN PRINT : GOTO prmin:

50 REM

COLOR 7

PRINT "6. Показатель степени n1: "; : COLOR 11: INPUT "", n1%

COLOR 7

PRINT "7. Показатель степени n2: "; : COLOR 11: INPUT "", n2%

COLOR 7

PRINT "8. Показатель степени n3: "; : COLOR 11: INPUT "", n3%

COLOR 7

PRINT "9. Воздействие F: "; : COLOR 11: INPUT "", F

COLOR 7

PRINT "10. Шаг интегpиpования H : "; : COLOR 11: INPUT "", H

COLOR 15

PRINT " Значения введены корректно (Y/N) "; : COLOR 14: INPUT "";

da$

IF (da$ <> "y") AND (da$ <> "Y") THEN PRINT : GOTO 50

IF T = 0 THEN

PRINT "Введите Т,не pавное нулю - ": PRINT : GOTO prmin

END IF

IF n1% = 0 THEN

PRINT "Введите n1,не pавное нулю - ": PRINT : GOTO prmin

END IF

init:

M1% = INT(TZ / H) 'Кол-во элементов в массиве B()

FOR i% = 1 TO n1% + n2% + n3%

Y(i%) = 0

NEXT i%

Z = 0: x = 0: l% = 0: i2% = 0

FOR i% = 1 TO M1%

B(i%) = 0

NEXT i%

GOTO slct1

slct1:

COLOR 7

PRINT "1-вывод на дисплей, 2-печать в файл DIG.DAT "; : COLOR 11:

INPUT "", i1%

IF i1% = 2 THEN

OPEN "DIG.DAT" FOR OUTPUT AS #1

PRINT #1, "**************************************"

PRINT #1, "** Расчет пеpеходной функции модели **"

PRINT #1, "**************************************”

PRINT #1, "Переходная функция модели записывается в виде"

PRINT #1, " K*exp(-s*TZ) "

PRINT #1, " W(s)=________________________________ "

PRINT #1, " (T*s+1)^n1*(A1*T*s+1)^n2*(A2*T*s+1)^n3 "

END IF

IF i1% <> 1 THEN PRINT : PRINT "Redo from start": GOTO slct1

GOTO interv

interv:

COLOR 7

PRINT "11. Шаг печати : "; : COLOR 11: INPUT "", W1

COLOR 7

PRINT "12. Kонечное вpемя интегpиpования: "; : COLOR 11: INPUT

"", W2

COLOR 15

PRINT " Значения введены корректно (Y/N) "; : COLOR 14: INPUT

""; da$

IF (da$ <> "y") AND (da$ <> "Y") THEN PRINT : GOTO interv:

IF W1 < H THEN

PRINT "Шаг печати не может быть. меньше шага интегpиpования"

GOTO interv

END IF

IF i2% = 1 THEN GOSUB rk2

GOTO round

prcout:

PRINT "X="; x, "Z="; Zarray(counter)

IF i1% = 2 THEN PRINT #1, " X="; x, "Z="; Zarray(counter)

IF x < W2 THEN GOSUB rk2: GOTO round

A$ = INPUT$(1)

GOTO graph

exgraph:

CLS

GOTO fin

round:

counter = counter + 1

MN = 1: Z1 = ABS(Z)

IF Z1 <> 0 THEN

WHILE Z1 < 1

Z1 = 10 * Z1

MN = 10 * MN

WEND

Z1 = (INT(100 * Z1 + .5)) / (100 * MN)

IF Z >= 0 THEN Z = Z1 ELSE Z = -Z1

END IF

Zarray(counter) = Z

Xmax = x

GOTO prcout

graph:

min = Zarray(1): max = Zarray(1)

FOR c = 1 TO counter

IF min > Zarray(c) THEN min = Zarray(c)

IF max < Zarray(c) THEN max = Zarray(c)

CLS

VIEW (110, 110)-(550, 390)

LOCATE 2, 5

PRINT "График зависимости Z(x)"

LOCATE 6, 5

PRINT USING "##### . ####"; max

LOCATE 25, 5

PRINT USING "##### . ###"; min

LOCATE 26, 10

PRINT "0.000"

LOCATE 26, 65

PRINT USING "##### . ###"; Xmax

WINDOW (1, min)-(counter, max)

FOR i = 0 TO 10

LINE (i * counter / 10 + 1, min)-(i * counter / 10 + 1, max), 4

LINE (1, i * ABS(max - min)/10)-(counter, i * ABS(max - min) / 10),4

FOR i = 2 TO counter

LINE (i - 1, Zarray(i - 1))-(i, Zarray(i)), 15

SCREEN 9

COLOR 15, 1

GOTO exgraph

fin:

CLS

END

rk2: ' Метод Рунге-Кутта втоpого поpядка

m% = INT(W1 / H)

FOR j% = 1 TO m%

GOSUB difsys

FOR i% = 1 TO n1% + n2% + n3%

G(i%) = Y(i%) + H * P(i%) / 4

Y(i%) = Y(i%) + 2 * H * P(i%) / 3

NEXT i%

GOSUB difsys

FOR i% = 1 TO n1% + n2% + n3%

Y(i%) = G(i%) + 3 * H * P(i%) / 4

NEXT i%

IF TZ = 0 THEN Z = Y(n1% + n2% + n3%) ELSE GOSUB dead

NEXT j%

x = INT((x + H * m%) / H + .5) * H

RETURN

difsys: ' Уpавнения Коши

P(1) = (K * F - Y(1)) / T

FOR i% = 2 TO n1%

P(i%) = (Y(i% - 1) - Y(i%)) / T

NEXT i%

FOR i% = n1% + 1 TO n1% + n2%

IF A1 <> 0 THEN

P(i%) = (Y(i% - 1) - Y(i%)) / (A1 * T)

END IF

NEXT i%

FOR i% = n1% + n2% + 1 TO n1% + n2% + n3%

IF A2 <> 0 THEN

P(i%) = (Y(i% - 1) - Y(i%)) / (A2 * T)

END IF

NEXT i%

RETURN

dead: ' Звено чистого запаздывания

l% = l% + 1

Z = B(l%)

B(l%) = Y(n1% + n2% + n3%)

IF l% < M1% THEN RETURN

l% = 0

RETURN

'Center:

' Center text on the given row.

SUB center (row%, text$)

LOCATE row%, 41 - LEN(text$) / 2

PRINT text$;

END SUB

SUB Pausa

ZP:

PZ$ = INKEY$

IF PZ$ = "" GOTO ZP

END SUB

Приложение 4

Таблица 3.1.1

Вихревой

расходомер

V - образный шаро-

вой клапан

Труба

ВЕНТУРИ или

сужающее

устройство

Двухпозиционный

проходной клапан с

поршневым

приводом

Трубка ПИТО Перепускной клапан

с поршневым

приводом

Элемент

измерения

расхода

(диафрагма)

Трехходовой клапан

с поршневым

приводом

Откалиброванн

ая диафрагма с

патрубками

Трехходовой

шаровой клапан с

поршневым

приводом

Продолжение таблицы 3.1.1

Датчик расхода

(трансмиттер)

со встроенной

диафрагмой

Пневматическое

реле

Массовый

расходомер

Соединение из плав-

кого нейлона на

случай пожара

Ротаметр с

регулятором

Предохранительный

клапан

Ротаметр Клапан сброса дав-

ления или вакуума

(дыхательный

клапан)

Шаровой

регулирующий

клапан

с

пневмоприводо

м

Электрический

сигнал

Этот штурвал у

того, у кого

надо штурвал

Программная связь

или передача

данных

Редукционный

клапан с

внешним

отбором

Пневматическая ли-

ния

Редукционный

клапан с

внутренним

отбором

Капиллярная линия

Продолжение таблицы 3.1.1

Дроссельный

регулирующий

клапан

Электромагнитный

или акустический

сигнал

Шаровой

клапан с

пневмоприводо

м

Полевой прибор

Двухходовой

соленоидный

клапан

Щитовой прибор

Угловой

поршневой

клапан

Прибор за щитом

Шиберный

клапан с

поршневым

приводом

Прибор на местной

панели

Трехходовой

соленоидный

клапа

Прибор позади

местной панели

Пара

трехходовых

соленоидных

клапанов по

схеме 1оо2

Лампочка состояния

на местной панели

Продолжение таблицы 3.1.1

Четырехходово

й соленоидный

клапан

"Прибор" РСУ

Четырехпроход

ной с двумя

соленоидами

Возврат к

исходному

состоянию через

РСУ

Клапан САМ

FLEX

Команды Открыть и

Закрыть от РСУ

Шаровой

клапан с

поршневым

приводом

Деблокирующий

ключ на РСУ

Блокировочное

устройство на

агрегате

Логика блокировки

через РСУ

Двойной

шаровой клапан

Логика блокировки

через ПАЗ

Оперативная панель

системы ПАЗ

Пневматически

й

переключатель

Оперативная панель

системы ПАЗ

Электроп

невмопреоб-

разователь (1/Р)

Оперативная панель

системы ПАЗ

Фафурин В.А., Терюшов И.Н. Автоматизация технологических процессов и производств

Подождите немного. Документ загружается.