Производство по чертежам Подбор аналогов Цены производителя Оригинальная продукция в короткие сроки
INNERпроизводство и поставка промышленных комплектующих и оборудования
Отзыв ★★★★★ Будем благодарны за отзыв в Яндексе — это помогает нам развиваться Оставить отзыв →
Правовая информация Условия использования технических материалов и калькуляторов Правовая информация →
INNER
Контакты

Калькулятор настройки ПИД-регуляторов

Калькулятор настройки ПИД-регулятора (онлайн, 9 методов, симуляция)

Калькулятор настройки ПИД-регулятора

Расчёт коэффициентов Kp, Ti, Td по 9 методам (Зиглера-Никольса, Cohen-Coon, CHR setpoint/disturbance, IMC, SIMC, AMIGO, Tyreus-Luyben) для моделей FOPDT, SOPDT, IPDT. Симуляция переходного процесса замкнутой системы.

Калькулятор реализует канонические формулы по [1] Ziegler-Nichols 1942, [2] Cohen-Coon 1953, [3] Chien-Hrones-Reswick 1952, [4] Rivera-Morari-Skogestad 1986, [5] Skogestad 2003 (SIMC), [6] Astrom-Hagglund 2004 (AMIGO), [8] Tyreus-Luyben 1992. Результаты — начальные настройки для тонкой подстройки на реальном объекте.
1 Модель объекта управления
Тип модели FOPDT — самая распространённая модель (нагрев, ёмкости, теплообменники). SOPDT — каскадные объекты с двумя постоянными времени. IPDT — интегрирующие (уровень при постоянном притоке).
Передаточная функция показана в описании метода ниже.
2 Параметры объекта
Коэффициент передачи K Отношение изменения выхода к изменению входа в установившемся режиме. Например, для нагрева: K = 2 °C на 1% мощности.
Ненулевое число (±)
Постоянная времени T, с Время, за которое выход достигает 63.2% установившегося значения после ступенчатого воздействия. Для IPDT — не используется.
> 0 секунд
Время запаздывания L, с Транспортное запаздывание — задержка между изменением управляющего сигнала и реакцией объекта. Для трубопроводов: L = длина / скорость потока.
≥ 0 секунд
Вторая постоянная T₂, с Вторая постоянная времени для модели SOPDT. Меньшая из двух обычно. Применяется правило половины Skogestad: T_eff = T₁ + T₂/2, L_eff = L + T₂/2.
> 0 секунд
Типовые объекты:
Подставит типовые K, T, L для соответствующего класса объектов.
3 Метод настройки
Метод
Подсказка по выбору внизу.
Тип регулятора
ПИД даёт лучшее качество, но дифференциальная составляющая усиливает шум измерений.
Параметр λ (τc), с Желаемая постоянная времени замкнутой системы. Малое λ — быстрее, но менее устойчиво. Рекомендации: λ ≈ L (агрессивно), λ ≈ T (умеренно), λ > T (мягко).
Структура регулятора 2-DoF структура. β — вес уставки в пропорциональной части (0…1). γ — вес уставки в дифференциальной (0…1). β=γ=1 — классический ПИД (по умолчанию). β=γ=0 — I-PD (мягкий старт без бросков). γ=0 — derivative-on-PV (Brett Beauregard, Rockwell PIDE).
β (P-weight): 1.00
γ (D-weight): 1.00
Зиглера-Никольса (1942): классический эмпирический метод для FOPDT. Даёт quarter-decay ratio (~25% перерегулирование). Лучший выбор для L/T = 0.1…1. Для L/T > 1 рекомендуется Cohen-Coon.
4 Параметры симуляции переходного процесса
Время симуляции, с
Обычно 5–10·T
Уставка (SP)
Ограничение u, ± Физическое ограничение управляющего сигнала с anti-windup. Например, для клапана: ±100%. 0 = без ограничений.
0 = без ограничения
5 Реалистичные эффекты: шум, мёртвая зона, квантование (опционально)
Покажет, как настройка ведёт себя в реальных условиях — шумный датчик, мёртвая зона привода, квантование АЦП. Помогает увидеть, что «идеальная» Z-N может оказаться непригодной из-за дёрганья клапана.
Шум измерения PV, СКО Среднеквадратическое отклонение гауссовского шума на измерении PV. Например, термопара даёт шум ±0.5–2 °C, расходомер — ±1–3% от шкалы.
0 = без шума
Мёртвая зона привода Управляющий сигнал |u| меньше этого порога не воздействует на объект. Моделирует люфт редуктора, трение клапана, гистерезис симистора. Типично 1–5% от полной шкалы.
0 = без мёртвой зоны
Шаг квантования датчика Дискретность АЦП/датчика. Например, термометр с разрядностью 0,1 °C → шаг 0.1. 12-bit АЦП на диапазоне 0–100 → шаг 100/4096 ≈ 0.024.
0 = непрерывный
Ожидание ввода данных
Заполните параметры объекта
Коэффициенты регулятора Kp PB%
Пропорциональный Kp
Время интегрирования Ti, с
Время дифференцирования Td, с
Ki = Kp/Ti, 1/с Альтернативная запись интегральной составляющей (используется в Arduino-библиотеках и некоторых PLC). Для интегрирующих объектов (IPDT) и медленных процессов Ki = 10⁻³…10⁻⁵ — это нормально.
Kd = Kp·Td, с Альтернативная запись дифференциальной составляющей. Если Td = 0 (П- или ПИ-регулятор), то Kd = 0.
Переходный процесс при ступенчатой уставке
Уставка (SP) Выход (PV) Управление (u) Измеряемый PV (шум/квантование)
Перерегулирование
Время регулир., с
Время нараст., с
IAE∫|e|dt — суммарное отклонение. Чем меньше, тем точнее регулирование.
ISE∫e²dt — сильнее штрафует большие отклонения. Полезно когда нельзя выходить за допуск.
ITAE∫t·|e|dt — штрафует поздние ошибки. Хорошо отражает экономический ущерб от долгого простоя вне уставки.
ITSE∫t·e²dt — комбинирует штраф за большие и за поздние ошибки. Самая строгая метрика.
Диагностика и рекомендации показать ▼
Автоанализ текущей настройки и типовых проблем. Кнопки «Применить» подставят соответствующие параметры и пересчитают.
У меня в реальности другой симптом:
Анализ робастности: устойчивость к ±20% разбросу параметров объекта показать ▼
Реальные объекты «плывут» — коэффициент передачи K и постоянная времени T меняются со временем (старение, нагрузка, износ). Таблица показывает, что произойдёт с переходным процессом при изменении параметров на ±20%.
Сценарий K T, с L, с Перерег., % tуст, с IAE Оценка
Сравнение всех методов
Метод Kp Ti, с Td, с Overshoot, % tset, с IAE
Зелёным — метод с наименьшим IAE на данном объекте. Симуляция выполняется при текущей уставке.
И Идентификация FOPDT по переходной характеристике
Величина ступени входа (Δu)
Размер ступенчатого изменения входа в эксперименте
Данные переходного процесса (t, y)
Каждая строка: время (с) и значение выхода через пробел или запятую. Можно вставить из Excel/CSV. После идентификации параметры подставятся в основную форму.
Введите данные и нажмите «Идентифицировать»

Онлайн-калькулятор настройки ПИД-регулятора (пропорционально-интегрально-дифференциального) рассчитывает коэффициенты Kp, Ti, Td по 9 каноническим методам — Зиглера-Никольса, Cohen-Coon, Chien-Hrones-Reswick, IMC, SIMC, AMIGO, Tyreus-Luyben — для трёх типов моделей объекта: FOPDT (первого порядка с запаздыванием), SOPDT (второго порядка), IPDT (интегратор с запаздыванием). Параметры подходят для настройки температуры, давления, расхода, уровня, оборотов двигателя, положения клапанов и задвижек, и большинства промышленных контуров автоматического регулирования.

В отличие от упрощённых форм, калькулятор включает симулятор переходного процесса замкнутой системы — после расчёта коэффициентов вы сразу видите, как поведёт себя объект при ступенчатом изменении целевого значения (уставки, setpoint). Показываются: перерегулирование, время регулирования (по полосе ±5%), время нарастания и интегральные ошибки IAE/ISE. В режиме «Сравнить методы» одновременно симулируются все 9 методов с подсветкой того, который даёт минимальную IAE — это надёжнее, чем выбирать метод по названию.

Полученные значения Kp, Ti, Td и альтернативная запись Ki = Kp/Ti, Kd = Kp·Td подходят для большинства промышленных контроллеров и измерителей-регуляторов: Овен ТРМ (ТРМ10, ТРМ12, ТРМ212, ТРМ251), REX C-100, Elhart, Термодат, преобразователей частоты, ПЛК на CodeSys, измерительных приборов с ПИД-функцией, а также для реализации на Arduino и микроконтроллерах с пересчётом единиц по форме записи (см. раздел про реализацию).

Принцип работы ПИД-регулятора простыми словами (для чайников)

Если кратко: ПИД-регулятор — это «умный кран», который сам подстраивает подачу управляющего сигнала, чтобы измеряемая величина (например, температура в печи) всегда оставалась на заданном уровне. Разберём как работает ПИД-регулятор шаг за шагом.

Сначала разберёмся с терминологией. Уставка (на английском setpoint, сокращённо SP) — это инженерное название того значения, которое регулятор должен удерживать. Например, если вы хотите поддерживать температуру 80 °C, то 80 °C — это уставка. Текущее измеренное значение датчика (фактическая температура в данный момент) на инженерном языке называется PV (process variable, регулируемая величина). Разница SP − PV — это ошибка регулирования (или просто ошибка, error). Задача регулятора — свести ошибку к нулю.

ПИД-регулятор — это устройство (или программа в контроллере), которое автоматически удерживает регулируемую величину (температуру, давление, расход, уровень, обороты) на заданной уставке. Расшифровка аббревиатуры: Пропорционально-Интегрально-Дифференциальный. Это значит, что управляющий сигнал на исполнительный механизм (клапан, нагреватель, частотный преобразователь) формируется как сумма трёх составляющих, каждая из которых реагирует на ошибку по-своему.

Бытовая аналогия — водитель за рулём: представьте, что вы едете по шоссе и хотите держать скорость 90 км/ч. Здесь 90 км/ч — желаемая скорость (это и есть «уставка»), фактическая скорость на спидометре — измеряемая величина, нажатие на педаль газа — управляющий сигнал. П-составляющая — вы давите на газ пропорционально тому, насколько сейчас медленнее: отстаёте на 10 км/ч → нажимаете сильнее, отстаёте на 2 км/ч → слабее. И-составляющая — вы замечаете, что машина едет медленнее желаемой скорости уже долго (например, в горку), и постепенно дожимаете педаль глубже, пока скорость не выйдет к 90 — устраняете накопившееся отставание. Д-составляющая — видите, что скорость начала быстро расти (поехали под уклон), и заранее отпускаете газ, чтобы не перерегулировать. Все три действия вы делаете одновременно — это и есть ПИД-регулирование.

Формально управляющий сигнал u(t) считается так:

u(t) = Kp · e(t) + (Kp/Ti) · ∫₀ᵗ e(τ) dτ + Kp · Td · de(t)/dt

где e(t) = SP − PV — ошибка регулирования (setpoint минус process variable), Kp — пропорциональный коэффициент, Ti — время интегрирования (за какое время интегральная часть повторит вклад пропорциональной), Td — время дифференцирования (на сколько вперёд экстраполируется тренд ошибки).

Чем П, ПИ и ПИД-регуляторы отличаются на практике

ТипСоставГлавное свойствоГде применяют
П-регулятортолько KpПростой, быстрый, но даёт статическую ошибку — PV никогда не достигнет SP, останется небольшое смещение. Чем больше Kp, тем меньше смещение, но при определённом значении система пойдёт в колебания.Когда статическая ошибка не критична (например, ручная коррекция, грубое регулирование уровня сыпучих материалов)
ПИ-регуляторKp + TiУстраняет статическую ошибку за счёт интегральной составляющей. Большинство промышленных контуров — это именно ПИ, а не полный ПИД.Регулирование расхода, давления, уровня, частоты вращения. Везде, где есть шум измерения (расходомеры, датчики давления)
ПИД-регуляторKp + Ti + TdДобавляет «опережение» через производную — система быстрее реагирует на изменение тренда ошибки, меньше перерегулирует.Регулирование температуры (медленные объекты), уровня в больших ёмкостях, позиционирование. Не применяется там, где сильный шум измерения
Почему ПД-регуляторы (без интегральной составляющей) применяют редко: без И-составляющей система всегда будет иметь статическую ошибку, как и чистый П-регулятор. ПД-регулирование используется в специальных случаях — например, в управлении электромеханическими позиционными системами, где обратная связь по положению и так компенсирует ошибку механически.

Зачем рассчитывать ПИД-настройки

ПИД-регулятор без правильной настройки — это случайный генератор. При завышенном Kp система пойдёт в автоколебания, оборудование начнёт стучать, клапаны изнашиваться втрое быстрее, продукт идти в брак. При заниженном Kp — управление будет вялым, уставка не достигается, технологический процесс плывёт. При неправильном Ti система либо перерегулирует и долго стабилизируется, либо никогда не выйдет в ноль ошибки. Дифференциальная составляющая Td при неверной настройке усилит шум измерительного канала и даст ложные броски управляющего сигнала.

Профессиональная настройка делается в три этапа: идентификация объекта (получение модели), расчёт начальных коэффициентов одним из канонических методов, тонкая подстройка на реальном объекте. Этот калькулятор закрывает второй этап и помогает с первым через встроенный модуль идентификации FOPDT по экспериментальной переходной характеристике.

Какие данные нужны

Минимум для расчёта по FOPDT — три параметра модели: коэффициент передачи K, постоянная времени T и время запаздывания L. Они извлекаются из экспериментальной переходной характеристики, снятой в разомкнутом режиме (когда регулятор переключён в ручной режим и подаётся ступенчатое изменение управляющего сигнала).

ПараметрЧто этоГде взять
K — коэффициент передачи объектаУстановившееся изменение выхода (PV) / изменение входа (CV)Δyуст / Δu из эксперимента
T — постоянная времени, сВремя достижения 63,2% от установившегося значения после окончания запаздыванияИз графика переходной характеристики
L — время запаздывания, сЗадержка между подачей входа и началом реакции выходаИз графика, либо встроенный модуль идентификации
T₂ — вторая постоянная времени (только для SOPDT)Для объектов с двумя ярко выраженными апериодическими звеньямиЧерез идентификацию по более сложным методам
Kv — интегральный коэффициент (только для IPDT)Скорость изменения выхода при единичном входеНаклон выходной кривой Δy/Δt при ступени Δu=1
Как получить переходную характеристику: переключите регулятор в ручной режим, дождитесь установившегося режима, дайте ступенчатое изменение управляющего сигнала на 5-10% (например, с 50% до 55%) и записывайте PV в течение времени, достаточного для установления нового стационара (обычно 3-5 постоянных времени). Запись через любую SCADA с интервалом 0,5-2 с — потом данные вставляются в калькулятор в режиме «Идентификация по характеристике».

Что считает калькулятор

Каждый параметр в результате отвечает за конкретный аспект качества регулирования:

ПараметрНазначение
Kp — пропорциональный коэффициентОсновной коэффициент, определяющий силу реакции на ошибку. Большой Kp → быстрая реакция, но риск колебаний.
Ti — время интегрирования, сЗа какое время интегральная составляющая повторит вклад пропорциональной. Малое Ti → быстрое устранение статической ошибки, но риск перерегулирования.
Td — время дифференцирования, сОпережение на основе скорости изменения ошибки. Подавляет перерегулирование, но усиливает шум.
Ki = Kp/Ti — интегральный коэффициентАльтернативная форма для параллельного представления (часто используется в PLC).
Kd = Kp·Td — дифференциальный коэффициентАльтернативная форма для параллельного представления.
Перерегулирование, %На сколько % выход превысил уставку. Допустимый предел зависит от процесса: для нагрева 5-20%, для расхода и давления — не более 5%, для уровня в баке — стремятся к 0%.
Время регулирования, сКогда выход окончательно вошёл в полосу ±5% от уставки и больше её не покидает.
Время нарастания, сВремя с 10% до 90% уставки. Характеризует быстродействие.
IAE — интегральная абсолютная ошибкаПлощадь под кривой |ошибки|. Чем меньше — тем лучше регулирование в среднем. Удобно для сравнения методов.

Модели объекта и формулы

FOPDT — первого порядка с запаздыванием

Базовая модель для 80% промышленных объектов: одно апериодическое звено плюс транспортное запаздывание. Описывается тремя параметрами K, T, L:

G(s) = K · e−L·s / (T·s + 1)

Применяется для: теплообменников, реакторов, сушилок, насосов, расходных линий с инерцией, температуры в печах. Если ваш объект ведёт себя «по S-образной кривой» — это FOPDT.

SOPDT — второго порядка с запаздыванием

G(s) = K · e−L·s / ((T₁·s + 1)(T₂·s + 1))

Применяется для каскадированных процессов: дистилляционные колонны, многоступенчатые теплообменники, системы с сильно различающимися инерциями. Калькулятор использует приём Skogestad half-rule: половина меньшей постоянной времени T₂ добавляется к L, а другая половина — к T₁. Это сводит SOPDT к эквивалентной FOPDT для применения тех же методов.

IPDT — интегратор с запаздыванием

G(s) = Kv · e−L·s / s

Применяется для объектов без самовыравнивания: уровень в баке при постоянном притоке, угол положения сервопривода, заряд ёмкости. Для таких объектов нельзя получить установившееся значение после ступени — выход растёт линейно. Тут требуются специальные методы: Z-N для IPDT, Tyreus-Luyben, AMIGO для IPDT и IMC с интегратором.

Все 9 методов настройки: как настроить ПИД-регулятор

Ниже подробно описаны 9 канонических методов настройки коэффициентов Kp, Ti, Td. Все они применяются по одной схеме: сначала идентифицируется модель объекта (FOPDT/SOPDT/IPDT — параметры K, T, L), затем по формулам метода вычисляются настройки регулятора. Калькулятор делает оба шага автоматически — нужно только ввести данные эксперимента или известные параметры объекта.

Метод 1: Зиглера-Никольса / Циглера-Никольса (Z-N, 1942)

Самый известный исторический метод (в литературе встречаются обе транслитерации фамилии Ziegler — «Зиглер» и «Циглер»). Даёт quarter-decay ratio — каждое следующее перерегулирование в 4 раза меньше предыдущего. Получается перерегулирование около 25-50% и довольно агрессивная настройка.

ПИД: Kp = 1,2 · T/(K·L); Ti = 2·L; Td = 0,5·L
ПИ: Kp = 0,9 · T/(K·L); Ti = L/0,3
П: Kp = T/(K·L)

Когда применять: для FOPDT с L/T = 0,1...1, когда допустимо умеренное перерегулирование ради быстродействия. Не годится для процессов, где перерегулирование опасно (давление, расход, чувствительные позиции).

Метод 2: Cohen-Coon (1953)

Улучшение Z-N для объектов с большим запаздыванием (L/T > 0,3). Чувствительнее к соотношению r = L/T.

ПИД: Kp = (1/K) · (T/L) · (4/3 + r/4)
Ti = L · (32 + 6r) / (13 + 8r); Td = L · 4/(11 + 2r)
ПИ: Kp = (1/K) · (T/L) · (0,9 + r/12)
Ti = L · (30 + 3r) / (9 + 20r)

Когда применять: для процессов с большим запаздыванием (длинные транспортёры, теплообменники со значительными трубопроводами, аналитические измерения с лабораторной задержкой).

Методы 3-6: Chien-Hrones-Reswick (CHR, 1952)

Четыре варианта: для slежения за уставкой (setpoint tracking) и для подавления возмущений (disturbance rejection), с 0% и с 20% перерегулирования. Это самая полная семья среди классических.

CHR setpoint 0% (PID): Kp = 0,6·T/(K·L); Ti = T; Td = 0,5·L
CHR setpoint 20% (PID): Kp = 0,95·T/(K·L); Ti = 1,4·T; Td = 0,47·L
CHR disturbance 0% (PID): Kp = 0,95·T/(K·L); Ti = 2,4·L; Td = 0,42·L
CHR disturbance 20% (PID): Kp = 1,2·T/(K·L); Ti = 2·L; Td = 0,42·L

Когда применять: CHR setpoint — когда задача в основном «следить за уставкой» (бэтч-процессы, изменение режимов). CHR disturbance — когда задача удерживать уставку при возмущениях (давление в трубопроводе при колебаниях расхода). Выбирайте 0% если перерегулирование недопустимо, 20% — если допустимо ради быстродействия.

Метод 7: Internal Model Control (IMC, 1986)

Модельный метод Rivera-Morari-Skogestad. В отличие от эмпирических Z-N и CC, IMC получает коэффициенты аналитически из условия желаемой формы замкнутой системы. Имеет настроечный параметр λ — желаемая постоянная времени замкнутого контура.

ПИД: Kp = (T + L/2) / (K · (λ + L/2))
Ti = T + L/2; Td = T·L / (2T + L)

Когда применять: когда нужно явно управлять компромиссом быстродействие-устойчивость через один параметр. Малое λ = быстро, но менее робастно. Большое λ = плавно, но медленно. Типичный выбор: λ = T (умеренный) или λ = max(0,1·T; L) (быстрый).

Метод 8: SIMC — Skogestad IMC (2003)

Уточнённая версия IMC, учитывающая ограничение интегрального времени.

ПИД: Kp = (1/K) · T/(τc + L)
Ti = min(T, 4·(τc + L)); Td = L/2

Когда применять: наиболее универсальный модельный метод. Рекомендуемый τc = L даёт компромисс между скоростью и робастностью. Хорошо работает для широкого диапазона L/T.

Метод 9: AMIGO (Aström-Hägglund, 2004)

Современный метод, оптимизированный по интегральной ошибке IAE при ограничении на максимальную чувствительность Ms ≤ 1,4. Даёт «робастные» настройки, устойчивые к ошибкам в модели.

ПИД: Kp = (1/K) · (0,2 + 0,45·T/L)
Ti = L · (0,4·L + 0,8·T)/(L + 0,1·T)
Td = 0,5·L·T / (0,3·L + T)

Когда применять: когда модель неточная или объект меняется со временем (старение, износ). AMIGO даёт чуть более консервативную настройку, чем Z-N, но гораздо лучше переносит ошибки идентификации.

Пример 1: система нагрева цеха

Параметры: тепловая инерция помещения, нагреватель 50 кВт

Из эксперимента получено: K = 2 °C/% (на каждый процент мощности +2°C в стационаре), T = 120 с (постоянная времени), L = 10 с (запаздывание датчика и теплопередачи). Соотношение L/T = 1/12 = 0,083 — малое запаздывание.

МетодKpTi, сTd, сПеререг., %tустан, сIAE
Z-N (классический)7,2205102,293,727,0
Cohen-Coon8,12523,783,58104,992,532,1
CHR setpoint 0%3,612050,341,617,3
CHR setpoint 20%5,71684,736,361,719,2
IMC (λ=10)4,1671254,81,139,814,8
SIMC (τс=10)3,08054,840,923,7
AMIGO2,845,454,8815,0123,728,1

Вывод: Z-N и Cohen-Coon дают перерегулирование более 100% — для системы отопления это много (потолок будет регулярно перегреваться). Лучшие варианты: IMC с λ=10 (overshoot 1,1%, IAE 14,8 — минимум среди всех методов) или CHR setpoint 0% (без перерегулирования, IAE 17,3). Для отопительного процесса, где важна плавность, рекомендуется IMC.

Что важно понять из примера: «классические» методы Z-N и Cohen-Coon были созданы в 1940-50-х для определённого класса задач (управление с quarter-decay ratio) и далеко не всегда оптимальны для современных требований к качеству. Калькулятор показывает все методы одновременно — выбирайте по тому, что подходит вашему процессу: для систем отопления и нагрева лучше IMC или CHR-0; для быстрого следования за уставкой — IMC или SIMC; для робастности — AMIGO.

Пример 2: регулирование расхода жидкости

Параметры: клапан + датчик расхода, прямая линия

Из эксперимента: K = 1,5 (т/ч)/% (каждый процент открытия даёт +1,5 т/ч), T = 8 с (инерция клапана и трубопровода), L = 2 с (датчик расхода). Соотношение L/T = 0,25 — умеренное запаздывание.

МетодKpTi, сTd, сПеререг., %tустан, сIAE
Z-N3,24,01,082,417,84,77
Cohen-Coon3,724,470,791,817,65,68
CHR setpoint 0%1,68,01,01,48,23,67
IMC (λ=2)2,09,00,891,07,92,98
SIMC (τс=2)1,338,01,01,210,04,33

Вывод: для расхода — где перерегулирование критично (импульсные нагрузки на оборудование, гидравлические удары) — Z-N и Cohen-Coon с перерегулированием 82-92% категорически не подходят. Оптимум: IMC с λ=2 (overshoot 1%, время установления 7,9 с, IAE 2,98). Для большинства расходных контуров рекомендуется использовать только ПИ-регулятор без дифференциальной составляющей — расходные датчики обычно зашумлены, и Td усилит шум.

Как читать переходный процесс

График в правой панели показывает реакцию объекта на ступенчатое изменение уставки. На нём три линии:

  • Серая (SP — setpoint) — заданная уставка, к которой система должна выйти. Обычно прямая горизонтальная линия (по умолчанию на уровне 1).
  • Тёмно-индиго (PV — process variable) — выход объекта, который реально получается с такими настройками регулятора. Это главная линия — она должна сойтись с SP.
  • Зелёная (u — управление) — управляющий сигнал, выдаваемый регулятором на исполнительный механизм. Полезно смотреть на его форму: резкие пики говорят об агрессивных настройках и износе клапанов.

Признаки правильной настройки

  • PV выходит к SP за 3-5 постоянных времени T объекта
  • Перерегулирование менее 5-20% (в зависимости от процесса)
  • 1-2 колебания и затем плавный выход в стационар
  • Управляющий сигнал u не имеет резких всплесков и не упирается в ограничения

Признаки плохой настройки

  • Незатухающие колебания — Kp слишком велик, либо Ti слишком мал. Уменьшайте Kp на 30-50%.
  • Очень медленный выход к SP — Kp слишком мал. Увеличьте в 1,5-2 раза.
  • Постоянное смещение от SP — отсутствие интегральной составляющей (Ti = ∞), либо Ti слишком велик. Уменьшите Ti.
  • Высокочастотные броски u — Td слишком велик, либо сильный шум измерения. Уменьшите Td или установите фильтр на вход.

Идентификация FOPDT по характеристике

Если параметры K, T, L неизвестны — их можно определить из экспериментальной переходной характеристики, снятой в разомкнутом режиме. Калькулятор использует метод Sundaresan-Krishnaswamy по двум точкам — времени достижения 28,3% и 63,2% от установившегося значения:

K = (y − y0) / Δu
T = 0,67 · (t63,2% − t28,3%)
L = 1,3 · t28,3% − 0,29 · t63,2%

Для использования модуля: перейдите на вкладку «Идентификация по характеристике», укажите величину ступени Δu и вставьте данные из SCADA или CSV-таблицы в формате «время значение» по одной паре на строку. Калькулятор автоматически определит K, T, L и предложит подставить их в основную форму.

Качество идентификации зависит от данных: ступенька должна быть достаточно большой (5-10% от диапазона), но не вводить объект в нелинейность; запись должна охватывать полное установление (3-5 постоянных времени); шум измерения должен быть минимален (можно усреднить по 3-5 точкам). Если переходная характеристика сильно отличается от S-образной формы — объект не является FOPDT, рассмотрите SOPDT или IPDT модели.

Какой метод выбрать

Не существует «универсально лучшего» метода — выбор зависит от характера объекта и требований к качеству:

ЗадачаРекомендуемый методПочему
Быстрое следование за уставкой, допустимо перерегулированиеZ-N или Cohen-CoonАгрессивная настройка, quarter-decay ratio
Плавный выход без перерегулирования (нагрев, уровень)CHR setpoint 0% или IMC с λ ≈ TГарантированно без overshoot
Подавление возмущений (давление, расход в магистрали)CHR disturbance или AMIGOОптимизированы под disturbance rejection
Большое запаздывание L/T > 0,5Cohen-Coon, IMC или SIMCЧувствительнее к запаздыванию
Малое запаздывание L/T < 0,1SIMC с τc > LИначе Z-N даст слишком агрессивный Kp
Неопределённость в модели, изменчивый объектAMIGO или IMC с большим λРобастные настройки
Универсальный выбор «по умолчанию»IMC с λ = max(0,1·T, L)Хороший компромисс для большинства случаев

Что делать, если результат не нравится

Перерегулирование слишком велико

  1. Сменить метод на менее агрессивный: Z-N → IMC или CHR-0
  2. В IMC — увеличить λ в 2-3 раза
  3. Уменьшить рассчитанный Kp на 20-30%
  4. Увеличить Ti в 1,5 раза

Слишком медленная реакция

  1. В IMC — уменьшить λ
  2. Сменить метод на более агрессивный: CHR-0 → CHR-20 → Z-N
  3. Увеличить Kp на 20-30%
  4. Уменьшить Ti на 30%

Колебательная неустойчивость

  1. Резко уменьшить Kp (в 2-3 раза)
  2. Увеличить Ti (в 2 раза)
  3. Если есть Td — обнулить
  4. Проверить корректность параметров модели — может быть, T идентифицирован неверно

Шум на выходе регулятора (резкие броски u)

  1. Уменьшить Td в 2-5 раз или обнулить
  2. Добавить фильтр первого порядка на входе измерения (постоянная фильтра 0,1·Td)
  3. Если используется производная — реализовать «derivative on PV» вместо «derivative on error»

Статическая ошибка (PV не приходит к SP)

  1. Проверить, что используется ПИД или ПИ, а не П-регулятор (Ti должно быть конечным)
  2. Уменьшить Ti в 2 раза
  3. Проверить ограничения на управляющий сигнал — возможно, исполнительный механизм упёрся в насыщение
  4. Проверить отсутствие нелинейностей (мёртвая зона клапана, гистерезис)

Применение: типовые объекты автоматического регулирования

Ниже приведены типовые классы промышленных объектов с типичными параметрами FOPDT-модели и рекомендуемым типом регулятора. Эти цифры — отправная точка для идентификации на реальном объекте, а не «универсальные настройки».

Объект управленияТип регулятораТиповые K, T, LРекомендуемый метод
Температура нагрева (печь, теплообменник, тэн, термокамера, экструдер, термодат)ПИДK = 1–5 °C/%, T = 60–600 с, L = 5–60 сIMC с λ ≈ T или CHR setpoint 0% (без перерегулирования, важно при перегреве продукта). Промышленные ПИД-терморегуляторы (Овен ТРМ10/12, REX C-100, Термодат) — это специализированные исполнения этой задачи
Расход жидкости/газа (магистральный трубопровод, дозирование)ПИ (без Td из-за шума)K = 0,5–2 (т/ч)/%, T = 5–20 с, L = 1–3 сIMC с малым λ или SIMC. Td не использовать — расходомеры зашумлены
Давление в трубопроводе или ёмкостиПИ или ПИДK = 0,5–3 бар/%, T = 5–30 с, L = 0,5–3 сCHR disturbance (главная задача — отрабатывать броски нагрузки)
Уровень в баке/резервуаре (открытый, с постоянным притоком — IPDT)ПИKv = 0,1–1 (м/с)/%, L = 2–10 сTyreus-Luyben или AMIGO для IPDT. Главное — без перерегулирования (перелив!)
Положение задвижки/клапана (электропривод, трёхходовой клапан)ПИ или ПИДK = 1 %/%, T = 5–30 с, L = 0,5–5 сIMC. Для трёхходового регулирующего клапана — обычная FOPDT-модель
Обороты двигателя (асинхронный двигатель с частотным преобразователем, шпиндель станка, насос с ЧРП)ПИK = 5–20 об/мин/%, T = 0,1–2 с, L = 0,01–0,1 сSIMC или Z-N PI. Без Td — энкодер дискретен и шумен
Ток нагрузки (управление электродвигателем по току, тиристорное регулирование, ИВЭП)ПИK = 0,5–2 А/%, T = 0,05–0,5 с, L = 0,01–0,05 сSIMC с малым τc — быстрый контур, ПИ-структура классическая для токовых петель
pH, концентрация (химическая дозация)ПИД с фильтромK сильно нелинеен, T = 30–300 с, L = 30–120 сCHR-0 setpoint, осторожно с настройками из-за нелинейности по диапазону pH
Скорость вентилятора (системы вентиляции, тяга в котлах)ПИK = 2–10 (м³/ч)/%, T = 3–30 с, L = 0,5–3 сSIMC. Часто каскад с регулятором температуры/давления выше
Калькулятор поддерживает быстрые пресеты для основных классов (кнопки «Нагрев», «Уровень в баке», «Расход», «Давление», «Колонна») — они подставляют типовые K, T, L. Это удобно как первое приближение и для тестирования различных методов на знакомых объектах.

Реализация ПИД-регулятора в контроллерах

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

Формы записи ПИД-регулятора

Существует несколько математически эквивалентных, но по-разному параметризованных форм:

ФормаУравнениеГде встречается
Параллельная (ideal/parallel)u = Kp·e + (Kp/Ti)·∫e dt + Kp·Td·de/dtКалькулятор, MATLAB, Simulink, большинство современных PLC, академическая литература
Параллельная по гейнам (gain form)u = Kp·e + Ki·∫e dt + Kd·de/dt, где Ki = Kp/Ti, Kd = Kp·TdArduino-библиотеки, простые микроконтроллерные реализации, некоторые PLC
Последовательная (interacting / series)u = Kp·(1 + 1/(Ti·s))·(1 + Td·s)·eСтарые пневматические регуляторы, некоторые старые PLC Honeywell, Foxboro
ISA-стандартнаяПараллельная с фильтром на Td: u = Kp·(e + 1/Ti·∫e dt + Td·de/dt/(1+Td/N·s))Промышленный стандарт ISA, Siemens PCS7, ABB 800xA

Калькулятор выдаёт настройки в параллельной (ideal) форме. Для контроллеров с другими формами нужны пересчёты — обычно они описаны в документации производителя.

Реализация в популярных платформах

ПлатформаФормаЕдиницы TiОсобенности
Овен ТРМ (ТРМ10, ТРМ12, ТРМ212, ТРМ251)ISA-параллельнаясекундыПрямой ввод Kp, Ti, Td. Поддерживает автонастройку. Для трёхпозиционного управления задвижками — ТРМ12
REX C-100 / C-700 / C-900ПараллельнаясекундыПараметры P, I, D. P задаётся как полоса пропорциональности (PB = 100/Kp в %). Имеет автонастройку (AT)
Siemens S7 (TIA Portal, PID_Compact)ISA-параллельнаясекундыФункциональный блок PID_Compact для S7-1200/1500. Поддерживает pretuning + fine-tuning
CodeSys (CODESYS 3.5)ПараллельнаясекундыГотовые библиотеки PID. На программируемых реле Овен ПР200, контроллерах WAGO, Beckhoff
Arduino / ESP32 (библиотека PID_v1 by Brett Beauregard)Гейновая (Ki, Kd)Используется Ki = Kp/Ti, Kd = Kp·Td. Сэмплирование задаётся отдельно
MATLAB / Simulink (PID block)Параллельная или ISAсекундыПрямая подстановка Kp, Ti, Td. Удобно для предварительной симуляции перед прошивкой
Преобразователи частоты (ABB ACS, Danfoss FC, Schneider Altivar, Веспер)ПараллельнаясекундыВстроенный ПИД для регулирования давления/расхода с обратной связью. Часто P задаётся как % от номинала

Подключение и интерфейсы регуляторов

На стороне выхода ПИД-регулятора (управляющий сигнал на исполнительный механизм) используются стандартные промышленные сигналы:

  • Токовый аналоговый 4–20 мА — наиболее распространённый стандарт. Помехоустойчив, работает на длинных линиях (до 1 км), позволяет обнаружить обрыв провода (ток меньше 4 мА — авария). Подключается к управляющим клапанам, регулирующим клапанам с электроприводом, тиристорным регуляторам мощности
  • Напряжение 0–10 В — для коротких линий и подключения к преобразователям частоты, диммерам, твердотельным реле
  • Релейный/транзисторный дискретный — для нагревателей через твердотельные реле (SSR), работает по принципу ШИМ (широтно-импульсная модуляция) с периодом обычно 1–10 секунд
  • Импульсный «больше/меньше» — для трёхпозиционного управления задвижками с электроприводом (Овен ТРМ12, ТРМ212)

На стороне входа — измерительные сигналы от датчиков: токовые 4–20 мА (давление, расход, уровень), термосопротивления Pt100/100М/50М, термопары типов K, J, L, S, потенциометры (для обратной связи положения), импульсные сигналы (расходомеры, тахометры).

Для связи с верхним уровнем АСУ ТП (SCADA, MES) современные регуляторы поддерживают цифровые интерфейсы — RS-485 с протоколом Modbus RTU, реже Ethernet (Modbus TCP, Profinet). Это позволяет менять уставки и коэффициенты регулирования удалённо, собирать архив переходных процессов и реализовывать каскадные схемы между несколькими контроллерами.

Структурная схема замкнутого ПИД-контура (упрощённо): уставка SP → сумматор (вычисление ошибки e = SP − PV) → ПИД-регулятор (Kp, Ti, Td) → ограничитель управляющего сигнала u → исполнительный механизм (клапан / нагреватель / ПЧ) → объект управления → датчик PV → обратная связь на сумматор. Эта структура одинакова и для одиночного контура, и для каждого контура в каскадных схемах, отличаются только величины K, T, L объекта и настройки регулятора.
Дискретная (цифровая) реализация ПИД-регулятора: в PLC и микроконтроллерах ПИД работает не непрерывно, а с дискретным шагом ΔT (sample time). Для устойчивости и точности соблюдайте правило ΔT ≤ T/10 — частота опроса должна быть минимум в 10 раз выше постоянной времени объекта. Для медленных тепловых процессов это легко (T = 120 с → ΔT = 12 с допустимо). Для быстрых контуров (ток двигателя, T = 0,05 с) нужен ΔT ≤ 5 мс. Также проверьте, что ваш контроллер поддерживает анти-windup (защита от насыщения интегратора при упоре исполнительного механизма в ограничения) — без неё регулятор может «зависнуть» на максимуме после возврата процесса.

Передаточная функция ПИД-регулятора

В операторной форме (s-плоскость Лапласа) передаточная функция параллельного ПИД-регулятора записывается так:

WR(s) = Kp + Kp/(Ti·s) + Kp·Td·s = Kp · (1 + 1/(Ti·s) + Td·s)

В замкнутом контуре управления с объектом WO(s) передаточная функция от уставки к выходу:

Φ(s) = WR(s)·WO(s) / (1 + WR(s)·WO(s))

Эта функция показывает, как система реагирует на изменения уставки. Если хочется отдельно анализировать реакцию на возмущение, используется передаточная функция от возмущения к выходу. Полюса замкнутой системы определяют устойчивость и характер переходного процесса. Калькулятор не вычисляет полюса аналитически, а получает переходный процесс прямой численной симуляцией — это наглядно и сразу даёт метрики качества (overshoot, IAE, время регулирования).

Частые вопросы

Когда нужна дифференциальная составляющая Td

Td полезна для процессов с большой инерцией и низким уровнем шума: температура, уровень, положение. Td даёт «опережение» — регулятор реагирует не только на ошибку, но и на скорость её изменения. Для быстрых и зашумлённых процессов (расход, давление в трубе) Td обычно не используют — она усилит шум измерительного канала. Правило простое: если измерение зашумлено (виден разброс ±2% и больше) — отключайте Td.

Что делать если L/T > 1 (большое запаздывание)

При очень больших запаздываниях классические методы дают неоптимальные настройки. Варианты: применить Cohen-Coon (он чувствительнее к L), использовать IMC или SIMC с большим λ или τc, рассмотреть предиктор Смита (Smith predictor) — это специальная структура, выносящая запаздывание из контура. Для очень больших L (десятки минут) часто проще принять более медленное регулирование и применить CHR setpoint 0% — главное чтобы система не была колебательной.

Почему мои реальные настройки в PLC отличаются от расчётных

Калькулятор использует стандартную параллельную форму ПИД: u = Kp·e + Kp/Ti·∫e dt + Kp·Td·de/dt. Некоторые PLC и контроллеры используют последовательную (interacting) или ISA-форму с другим определением коэффициентов. Также часто используют другие единицы — Ti в минутах вместо секунд, Ki как «repeats per minute». Сверьтесь с документацией вашего контроллера и применяйте соответствующие пересчёты. Если параметры просто не работают на реальном объекте — проверьте, что модель идентифицирована корректно (правильно ли определили K, T, L).

Почему методы Z-N и Cohen-Coon такие агрессивные

Они были разработаны в 1940-50-х годах для оптимизации по критерию quarter-decay ratio (каждое следующее перерегулирование в 4 раза меньше предыдущего), который тогда считался приемлемым. С появлением современных требований к качеству регулирования (минимум перерегулирования, точное соблюдение уставки) появились более совершенные методы IMC, SIMC, AMIGO. Z-N и Cohen-Coon сегодня — это в основном «исторические» методы для быстрой грубой настройки.

Можно ли использовать результаты сразу на реальном объекте

Результаты калькулятора — это начальные настройки, которые гарантируют устойчивость и приближают к оптимуму. Для финальной настройки практически всегда требуется подстройка на реальном объекте (так называемая «fine tuning»). Алгоритм: сначала подставьте рассчитанные коэффициенты, переведите регулятор в автоматический режим в безопасном режиме (например, при низкой нагрузке), дайте небольшое возмущение и наблюдайте отклик. Если перерегулирование слишком велико — снижайте Kp на 10-20%; если медленно — повышайте Kp на 10-20%; колебания — увеличивайте Ti. Итерация 3-5 циклов обычно достаточно.

В каком регистре уравнения PID я получу настройки

Калькулятор выдаёт коэффициенты в стандартной параллельной форме (ideal/parallel): u(t) = Kp·e(t) + Kp/Ti·∫e(t)dt + Kp·Td·de/dt. Это так называемая «академическая» форма, используемая в большинстве учебников и в Simulink/MATLAB. Если ваш контроллер использует последовательную (series, interacting) форму или ISA-стандарт — пересчёт делается по простым формулам, обычно описанным в документации контроллера. Для удобства в результатах также показываются Ki = Kp/Ti и Kd = Kp·Td — это «гейновая» форма, часто используемая в PLC.

Что такое IAE и зачем оно нужно

IAE (Integral Absolute Error) — это интегральная абсолютная ошибка, площадь под графиком |SP − PV(t)| за всё время переходного процесса. Это удобный единый числовой критерий качества: чем IAE меньше — тем регулирование лучше «в среднем». В режиме «Сравнить методы» калькулятор симулирует все 9 методов на одном объекте и подсвечивает зелёным тот, что даёт минимальный IAE — это эмпирически оптимальный выбор для данного объекта. Однако IAE — не единственный критерий: иногда важнее минимизировать перерегулирование (для механики), время установления (для бэтч-процессов) или энергию управления (для долговечности приводов).

Подходит ли калькулятор для каскадного управления

Калькулятор настраивает одиночный ПИД-контур. Для каскада сначала настройте внутренний контур (обычно расход или скорость) на быстрое слежение за уставкой методом IMC или CHR-0, а затем внешний контур (температура, уровень) — на 5-10 раз более медленное регулирование, считая внутренний контур «идеальным». Для внешнего контура постоянная времени внутреннего добавляется к его собственной L. Каскад даёт лучшее качество, чем одноконтурное регулирование, но требует двух разделённых настроек — оба раза с помощью этого калькулятора.

Чем адаптивный ПИД-регулятор отличается от обычного

Адаптивный ПИД самостоятельно подстраивает коэффициенты Kp, Ti, Td в процессе работы — например, при изменении нагрузки, температуры окружающей среды, износа исполнительных механизмов. Используются такие подходы как gain scheduling (несколько наборов настроек для разных режимов с переключением между ними) или непрерывная подстройка по идентификации модели в реальном времени. Промышленные адаптивные ПИД встречаются в Siemens SIMATIC, ABB 800xA, а простые формы автонастройки (auto-tuning, AT) есть даже в недорогих регуляторах Овен ТРМ и REX C-100. Этот калькулятор даёт стационарные коэффициенты для конкретной модели — но они могут использоваться как стартовые в адаптивных системах.

В чём разница между ПИД, нечётким (fuzzy) и релейным регулированием

Релейный (двухпозиционный) регулятор включает/выключает исполнительный механизм при пересечении уставки с гистерезисом — это самый простой способ управления, годится только для медленных процессов с большой инерцией (бытовые термостаты). Нечёткий (fuzzy logic) регулятор использует лингвистические правила вида «если ошибка большая положительная, то управление сильно увеличить» — он эффективен для нелинейных объектов с трудно формализуемой моделью, но требует экспертной настройки правил. Классический ПИД остаётся базовым решением для линейных и слабо нелинейных объектов с математической моделью — даёт оптимальное по IAE/ISE регулирование при правильно подобранных коэффициентах. На практике 90%+ промышленных контуров — это ПИД или его варианты.

Что такое автонастройка (auto-tuning) и заменит ли она ручную

Автонастройка — встроенная в регулятор функция, которая сама проводит эксперимент с объектом (обычно метод релейной обратной связи Aström-Hägglund 1984) и вычисляет настройки. Поддерживается большинством современных приборов: Овен ТРМ (АТ), REX C-100 (AT), Siemens PID_Compact (pretuning + fine-tuning). Это удобно для линейных стабильных объектов, но имеет ограничения: эксперимент происходит «вслепую» от пользователя, не учитывает специфику процесса (например, недопустимость колебаний температуры выше допуска), не выбирает между методами с разными приоритетами (быстродействие vs устойчивость). Расчёт через калькулятор по идентифицированной FOPDT-модели даёт больше контроля и понимания происходящего — особенно если параметры объекта получены из исторических данных без проведения активного эксперимента.

Дисклеймер: калькулятор предназначен для расчёта начальных настроек ПИД-регулятора и образовательных целей. Реальная настройка автоматических систем регулирования на промышленных объектах требует учёта множества факторов — нелинейностей объекта, ограничений исполнительных механизмов, специфики измерений, требований безопасности и регламентов производства. Результаты калькулятора должны рассматриваться как стартовые приближения, требующие подтверждения на математической модели и последующей подстройки на реальном объекте квалифицированным специалистом. Калькулятор не предназначен для применения в системах функциональной безопасности (SIL 1-4) без дополнительной валидации.
Источники и нормативная база:
  1. Ziegler J.G., Nichols N.B. Optimum Settings for Automatic Controllers // Transactions of the ASME, 1942, Vol. 64, pp. 759-768.
  2. Cohen G.H., Coon G.A. Theoretical Consideration of Retarded Control // Transactions of the ASME, 1953, Vol. 75, pp. 827-834.
  3. Chien K.L., Hrones J.A., Reswick J.B. On the Automatic Control of Generalized Passive Systems // Transactions of the ASME, 1952, Vol. 74, pp. 175-185.
  4. Rivera D.E., Morari M., Skogestad S. Internal Model Control. 4. PID Controller Design // Industrial & Engineering Chemistry Process Design and Development, 1986, Vol. 25(1), pp. 252-265.
  5. Skogestad S. Simple analytic rules for model reduction and PID controller tuning // Journal of Process Control, 2003, Vol. 13, pp. 291-309.
  6. Aström K.J., Hägglund T. Revisiting the Ziegler-Nichols step response method for PID control // Journal of Process Control, 2004, Vol. 14, pp. 635-650.
  7. Tyreus B.D., Luyben W.L. Tuning PI Controllers for Integrator/Dead Time Processes // Industrial & Engineering Chemistry Research, 1992, Vol. 31(11), pp. 2625-2628.
  8. O'Dwyer A. Handbook of PI and PID Controller Tuning Rules. 3rd Edition. Imperial College Press, 2009. — 624 p.
  9. Aström K.J., Hägglund T. PID Controllers: Theory, Design, and Tuning. 2nd ed. ISA, 1995. — 343 p.
  10. Денисенко В.В. Компьютерное управление технологическим процессом, экспериментом, оборудованием. — М.: Горячая линия–Телеком, 2009. — 608 с.
  11. Ротач В.Я. Теория автоматического управления. 5-е изд. — М.: МЭИ, 2008. — 396 с.
  12. Vassiljeva K. ISS0065 Control Instrumentation. Lecture 8: PID Tuning. Tallinn University of Technology, 2013.

Заказать товар

ООО «Иннер Инжиниринг»