Описание и анализ объекта автоматизации. Анализ технологического процесса как объекта автоматизации, страница 16

           Δh(P) = W11(P) ΔQ42(P)-W12(P) ΔQ43(P)-W13(P) ΔQ44(P)+W17(P) ΔG41(P);

Δρ43(P) = -W21(P) ΔQ42(P)-W23(P) ΔQ44(P)+W24(P) Δρ42(P)-W26(P) Δρ44(P)+W27(P) ΔG41(P)-W28(P) Δh(P),

(3.1.24.)

где W11(P) = K11/K25ThP; W12(P) = K12/K25ThP; W13(P) = K13/K25ThP; W17(P) = K17/K25ThP; W21(P) = K21/K25(Tρ43P+1); W23(P) = K23/K25(Tρ43P+1); W24(P) = K24/K25(Tρ43P+1); W26(P) = K26/K25(Tρ43P+1); W27(P) = K27/K25(Tρ43P+1); W28(P) = K28/K25(Tρ43P+1); Tρ43 = 1/ K25; Th = 1/ K25.

Экспериментально-аналитический метод построения математической модели.

Первый этап - получение экспериментальных данных (переходных характеристик), (кривые разгона) – характеристика объекта – снимаются при отключенном регуляторе.

Второй этап – обработка экспериментальных данных.

Основные каналы: 1. – по управлению (степень открытия клапана белого щёлока).    2. – по возмущению (производительностькотла(расход плава в бак растворитель плава)).   

По второму этапу существует два метода: аналитический (метод площадей или метод симаю) и графический метод (графо-аналитический метод).

I.По виду переходной характеристики определяют вид дифференциального уравнения (порядок дифференциального уравнения): если экспонента, то первый порядок, если есть перегиб, то порядок выше первого.

II. Аппроксимация или приближённая замена объектов порядка более высокого, чем первый объектом первого порядка с запаздыванием:

T*(dθ(t)/dt)+ θ(t)=КуSупр(t-Ţ)

III. Определение коэффициентов управления (Т, Ку).

Т - определяется по графику, а Ку по формуле ∆D/∆m (г/см3/%), где ∆D и ∆m соответственно изменение плотности зелёного щёлока и изменение степени открытия клапана белого щёлока.

IV. Записываем передаточную функцию, согласно имеющегося уравнения:                     -рŢ 

Wу(Р) = Ку/(Тр+1)*е

Общий вид модели получили аналитическим методом, а коэффициенты модели рассчитываются по экспериментальным данным из которых мы получаем исходные данные для программы расчёта коэффициентов модели объекта.

По разгонным характеристикам, снятым на СРК во время прохождения  преддипломной практики, строим график зависимости плотности зелёного щёлока от степени закрытия клапана подачи слабого белого щёлока:

«Таблица 3.1.1.» Разгонные характеристики.

№ п/п

t усл. сек

ρ I

ρ i- ρo

hi=( ρ i- ρo)/ (ρ ∞- ρo)

1- hi

1

2

3

4

5

6

1

0

1,1

0

0

1

2

118

1,11

0,01

0,213

0,787

3

139

1,12

0,02

0,426

0,574

4

182

1,13

0,03

0,638

0,362

5

300

1,137

0,037

0,787

0,213

6

600

1,149

0,049

1,043

-0,043

7

900

1,151

0,051

1,085

-0,085

8

1200

1,157

0,057

1,212

-0,212

9

1500

1,153

0,053

1,127

-0,127

10

1800

1,15

0,05

1,063

-0,063

11

2100

1,147

0,047

1

0

12

2400

1,147

0,047

1

0

13

2700

1,147

0,047

1

0

14

3000

1,147

0,047

1

0

Рисунок 3.1.1. Изменение степени открытия клапана белого щёлока во времени.