Автоматизация процесса флотации сульфидной золотосодержащей руды, страница 32

t,c

Рисунок 18 – Переходной процесс с оптимальными настройками ПИД-регулятора

3.4  Проверка АСР расхода сжатого воздуха на устойчивость

После определения настроек регулятора исследуем данную АСР на устойчивость, а также определим запас устойчивости системы по амплитуде  и по фазе, используя частотный критерий Найквиста. Найдем амплитудно-фазовую характеристику (АФХ) объекта. Ее получают подстановкой р= jw в передаточную функцию разомкнутой системы W(p).

Передаточная функция объекта имеет вид:

тогда

Определяем сначала передаточную функцию объекта без учета запаздывания  W0(jω).

Выделим вещественную и мнимую части. Для этого домножим числитель и знаменатель  Wo(jω) на сопряженное знаменателю выражение:

Выражения для вещественной и мнимой частей, соответственно равны:

          

Учтем запаздывание в системе:

Тогда

(14)

Получаем

Рассчитаем передаточную функцию регулятора:

(3.15)

Выделим вещественную и мнимую части и подставим значения оптимальных настроек регулятора Кп  = 0,148 м3/ч,  Ти = 2,34 c,  Тд = 0,659 c:

,

Передаточную функцию разомкнутой системы получим как произведение передаточной функции объекта и регулятора

Выделим вещественную и мнимую части:

Pраз(w) = Pоб (w)PР(w) - Qоб (w)Qp(w),

Qраз(w ) = Qоб (w)Pp(w) + P об (w)Qp(w).

Задаваясь значениями частоты ωот 0 до ∞, используя ПП Excel находим числовые значения Pоб(ω), Qоб(ω), Pр(ω), Qр(ω), Рраз(ω) и Qраз(ω), данные сводим в таблицу 10. По результатам таблицы на комплексной плоскости строим АФХ системы (рисунок 19).

Таблица 10 – АФХ объекта, регулятора и разомкнутой АСР

ω, рад/с

0

0,45

0,6

0,75

0,9

1,05

1,2

1,35

Pоб(ω)

0

-0,7801

-3,4089

-4,5943

-4,892

-4,6463

-4,0702

-3,3036

Qоб(ω)

-∞

-12,587

-10,205

-7,7414

-5,561

-3,7345

-2,25

-1,0756

Pp(ω)

0,148

0,148

0,148

0,148

0,148

0,148

0,148

0,148

Qp(ω)

-∞

-0,0967

-0,0469

-0,0112

0,0175

0,04217

0,06433

0,08482

Pраз(ω)

-∞

-1,3322

-0,9831

-0,7665

-0,6267

-0,5302

-0,4576

-0,3977

Qраз(ω)

-∞

-1,7874

-1,3504

-1,0943

-0,9086

-0,7486

-0,5948

-0,4394

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

ω, рад/с

1,5

1,65

1,8

1,95

2,1

2,25

2,4

2,55

Pоб(ω)

-2,4452

-1,5685

-0,7295

0,02914

0,67594

1,18979

1,559

1,7806

Qоб(ω)

-0,177

0,47851

0,92215

1,18392

1,29305

1,27812

1,16689

0,9859

Pp(ω)

0,148

0,148

0,148

0,148

0,148

0,148

0,148

0,148

Qp(ω)

0,10413

0,12259

0,14042

0,15775

0,1747

0,19134

0,20772

0,2239

Pраз(ω)

-0,3435

-0,2908

-0,2374

-0,1825

-0,1259

-0,0685

-0,0117

0,04278

Qраз(ω)

-0,2808

-0,1215

0,03405

0,17982

0,30946

0,41681

0,49654

0,54459

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