Идентификация объекта управления

В данном разделе выполнен расчёт автоматической системы регулирования соотношения руда и расхода воды в мельницу, ([8], [11]).

Расчет АСР начинается с составления математического описания объекта управления, при котором определяются его статические (коэффициент усиления - К_об) и динамические (запаздывание - t, постоянные времени – Т) параметры.

Расчет ведем при внешнем возмущающем воздействии DX_вх =10 % хода регулирующего органа.

Показатели качества регулирования, определяемые технологическим процессом:

- статическая ошибка DG_cт= 0;

- максимальное динамическое отклонение DG₁ £ 2,5 м³/ч;

- время регулирования t_р ≤ 8 с;

- переходной процесс с 5% перерегулированием.

Данные кривой разгона приведены в таблице 3.11.

Таблица 3.11 – Данные кривой разгона

График кривой разгона и возмущающего воздействия DX_вх (t), % ХРО представлен на рисунке 3.20, 3.21.

Рисунок 3.20 – График возмущающего воздействия X_вх (t), % ХРО

DG(t)

Рисунок 3.21 – График кривой разгона объекта ∆G(t), м³/ч

Единичная DG⁰(t) и нормированная DG^Н(t) переходные функции определяются по следующим формулам:

DG⁰(t) = DG(t) / A, (16)

DG^Н(t) = DG⁰(t) / DG⁰(T_у), (17)

где A – скачкообразное возмущающее воздействие, при котором снята переходная характеристика (А = 3 % ХРО); T_у – время переходного процесса; DG^Н(T_у) – установившееся значение переходной характеристики; DG⁰(T_у) = 0,91 м³/ч.

Результаты расчётов единичной и нормированной переходных функций сведены в таблицу 3.12.

Таблица 3.12 – Ординаты переходных функций

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

По результатам расчетов таблицы 3.12 строю единичную и нормированную переходные характеристики, рисунок 3.22.

DG(t),

1 - ΔG(t)⁰ – еденичная; 2 - ΔG(t)^н – нормированная

Рисунок 3.22 – Переходные характеристики

Из зависимости DG⁰(t) находится величина коэффициента усиления объекта K_об:

K_об = DG⁰(T_у) = 0,91. (18)

При определении динамических свойств объекта, сам объект представляется последовательным соединением нескольких инерционных звеньев первого порядка и одного звена чистого запаздывания. Исходя из экспериментальных данных, величина запаздывания объекта составляет 0,5 секунду. Уточним это значение, рассчитав величину дополнительного запаздывания τ_д. Для этого возьмем на графике нормированной переходной функции точки t_A и t_B,

Для этого на рисунке 3.22, по нормированной переходной характеристике, графически определим время t_a и t_b, для ΔG^н_а(t)=0,1 и ΔG^н_b(t)=0,7 соответственно:

t_a = 1,2 – 0,5 = 0,7с и t_b = 4,36 – 0,5 = 3,86с.

Определим динамические характеристики объекта при аппроксимации его последовательным соединением апериодического звена и звена запаздывания.

Определим величину дополнительного запаздывания:

, (19)

;

Общее запаздывание объекта τ_об:

t_об = t + t_д = 0,5 + 0,397 = 0,898 с; (20)

Определим постоянную времени Т_об:

с. (21)

Таким образом, аппроксимирующая передаточная функция:

. (22)

Для определения точности аппроксимации экспериментальной переходной функции решением дифференциального уравнения первого порядка с запаздывающим аргументом рассчитываются ординаты аппроксимирующей кривой:

(23)

(24)

Для определения среднеквадратической ошибки аппроксимации вычисляется отношение d_Ап:

(25)

Результаты расчетов приведены в таблице 3.13.

Таблица 3.13 – Ординаты переходных функций

По данным таблицы 3.13 строю график аппроксимирующей кривой показан на рисунке 3.23.

DG(t)

1 – при аппроксимации решением дифференциального уравнения первого порядка;

2 – нормированная сглаженная

Рисунок 3.23 – Переходные характеристики

По данным таблицы 3 рассчитываем среднеквадратичную ошибку аппроксимации δ:

%. (26)

Погрешность аппроксимации удовлетворительная, так как не превышает d < 3 %.

Построим кривую переходного процесса в среде MATLAB.

Зададим передаточную функцию объекта регулирования:

>> Wо=tf(0.91,[5,732 1],’inputdelay’,0.897)

Transfer function:

0.91

exp(-0.897*s) * -----------

5.732 s + 1

строим кривую переходного процесса:

>> step(Wо)

DG(t)

Рисунок 3.24 – Кривая переходного процесса объекта регулирования

Из полученного графика видно, что переходной процесс является монотонным, оптимальным, перерегулирования нет.