Arduino.ru

ПИ регулятор отопления

• Войдите или зарегистрируйтесь, чтобы получить возможность отправлять комментарии

Регулятор был задуман как замена штатного в электроконвектор на дачу, потому что встроенный не обеспечивал необходимой точности поддержания температуры и при перебоях в питании самостоятельно не включался в работу, что могло нехорошо отразится на заполненом водопроводе.

Данный регулятор обеспечивает поддержание заданной температуры. Проект начинался тут.

Проработал зиму на даче, никаких замечаний к работе регулятора нет, температуру держет точно.

В качестве датчика температуры используется терморезистор NTC 10кОм при 20С, выходным устройством являеться реле (в дальнейшем планирую поставить симистор) (организован медленный ШИМ частотой 2Гц).

Схема включения терморезистора стандартная, кнопки (аналоговые). Изскетча должно быть понятно что куда цеплять. У меня девайс выполнен на своей плате, использован микроконтроллер ATMega8, тактирование от встроенного генератора 8МГц. Собрано всё в монтажной коробке 100х100х40.

Схема устройства, LCD экран не показан.

• Войдите или зарегистрируйтесь, чтобы получить возможность отправлять комментарии

2 герца на реле.

2*60*60*24*365 = 63 ляма в год. реле столько не ходят, обычно от 200 до 500 тыс, я когда печку делал себе посчитал 1 переклюяение в 20 секунд и то вышло всего на 5 лет 8 часовой работы.

• Войдите или зарегистрируйтесь, чтобы получить возможность отправлять комментарии

2 герца на реле.

Вместо реле ставим SSR на нужную мощность и все.

• Войдите или зарегистрируйтесь, чтобы получить возможность отправлять комментарии

2 герца на реле.

Вместо реле ставим SSR на нужную мощность и все.

Согласен, я выше написал что планируется замена на симистор и подача напряжения по алгоритму Брезенхема. Опыт уже есть — видео.

• Войдите или зарегистрируйтесь, чтобы получить возможность отправлять комментарии

Такую же задачу решал недавно. Реализовал алгоритм «мячик на двух резинках» с поиском равновесной мощности под заданную температуру итеративным перебором. Регулирование мощности на моем преобразователе ШИМ в ФИМ

зима показала очень высокую точность поддержания температуры. Разве что помех много от симистора но у меня это никак не проявляется, знаю, что должны быть)

• Войдите или зарегистрируйтесь, чтобы получить возможность отправлять комментарии

Такую же задачу решал недавно. Реализовал алгоритм «мячик на двух резинках» с поиском равновесной мощности под заданную температуру итеративным перебором.

Ссылочку не дадите для ознакомления

• Войдите или зарегистрируйтесь, чтобы получить возможность отправлять комментарии

Не публиковал, но если что конкретно интересует поделюсь

• Войдите или зарегистрируйтесь, чтобы получить возможность отправлять комментарии

• Войдите или зарегистрируйтесь, чтобы получить возможность отправлять комментарии

Согласен, я выше написал что планируется замена на симистор и подача напряжения по алгоритму Брезенхема. Опыт уже есть — видео.

спасибо за информацию, не знал, что такая простая вещь как пропуски полупериодов так круто называются)

логика програмная пропусков почти ничем не отличается от ФИМ.

По мне управление мощностью дело совсем не сложное. Алгорит расчета %необходимой мощности вот это куда более сложная и интересная задача. ПИД ориентирован на релейное управление. А здесь как? Я использовал алгоритм резинки. В стабильных условиях это отлично работает, а вот при резко меняющихся условий резинка входит в колебательный процесс

• Войдите или зарегистрируйтесь, чтобы получить возможность отправлять комментарии

что именно? Схема, алгоритм?

• Войдите или зарегистрируйтесь, чтобы получить возможность отправлять комментарии

2 герца на реле.

2*60*60*24*365 = 63 ляма в год. реле столько не ходят, обычно от 200 до 500 тыс, я когда печку делал себе посчитал 1 переклюяение в 20 секунд и то вышло всего на 5 лет 8 часовой работы.

Я прикидывал 600000-800000 переключений за отопительный сезон 210 дней.

• Войдите или зарегистрируйтесь, чтобы получить возможность отправлять комментарии

что именно? Схема, алгоритм?

• Войдите или зарегистрируйтесь, чтобы получить возможность отправлять комментарии

2 герца на реле.

Вместо реле ставим SSR на нужную мощность и все.

Если точно то 1% — 0,6с. (100% за 60с.). Если нужна мощность 15% то реле включаеться на 9с. а остальные 51с. отключено. И так каждые 60с.

• Войдите или зарегистрируйтесь, чтобы получить возможность отправлять комментарии

хорошо. С компа пришлю примеры кода, а на словах примерно так:

— по таймеру беру равные интервалы времени, это шаг итерации

— на каждой итерации считаю три дельты — 1. Дельта между текущей и целевой температурой 2. Дельта между целевой температурой как прогноз на следующий период 3. Дельта как ошибка прогноза на текущий период на предыдущем цикле

— все три цифры суммирую с коэффициентами и получаю корректировку на этой итеррации для мощности в плюс или минус. Коэффициенты подбираются при отладке

— кроме того есть крайние условия 1. Если разница между текущей тем и целевой больше заданной константы (кажется у меня там 10 градусов) включаю полную мощность 2. Если разница меньше (у меня 0.5 градуса) то в эту итерацию мощность не меняется

— кроме того задается жесткость резинки. Запрещается слишком большое изменение мощности за несколько итерраций, это исключает автоколебания с большой амплитудой

а целом алгоритм показал очень хороший результат для меня. Его минус в необходимости подбора параметров под конкретную задачу

RadioTexnik

Блог по электронике

Необходимость построить ПИД регулятор температуры связана с желанием построить паяльную станцию для пайки безвыводных компонентов.

К сожалению, производители современной бытовой электронной аппаратуры всё чаще используют электронные компоненты на безвыводной пайке, так называемые BGA (англ. Ball grid array — массив шариков). Процесс демонтажа данных компонентов осложняется тем, что необходимо создать выдержанные значения и временные интервалы воздействия температуры в зоне пайки так называемых «шаров». Эту задачу легко решают промышленные образцы инфракрасных паяльных станций (например ТермоПро), но их стоимость не по карману обычному обывателю домашней мастерской.

К счастью, в настоящее время на рынке электроники, который пестрит готовыми платами с микроконтроллерами на борту, например ARDUINO. Вполне по силам собрать такую паяльную станцию своими силами. Для этого потребуется усидчивость, уверенное знание ПК, а также ПРЯМЫЕ РУКИ и горячий паяльник)))

На видео, помещенном в начале этой статьи, подробно показана работа простейшей системы параметрического автоматического регулирования температуры в зоне бесконтактной пайки инфракрасным излучением при помощи простейшей китайской термопары (по схеме TC1), готового модуля (опять же китайского) ARDUINO UNO (по схеме Arduino328), силового ключа на n-канальном MOSFEET транзисторе IRFZ-44N (Q1), который управляет мощностью накала автомобильной галогеновой лампы (L1). Данная система используется в качестве стенда для «обкатки» самостоятельно написанной программы в среде визуального программирования для микроконтроллеров AVR – FlowCode v5.3.0.0.

О возникших трудностях в повторении проекта, а также о своем личном опыте не стесняйтесь излагать в комментариях. Успехов! 😉

СКАЧАТЬ исходные файлы данного проекта, а также прошивку ARDUINO UNO

PID РЕГУЛЯТОР НА ARDUINO

ПИД регулятор – мощный инструмент, позволяющий удерживать заданную величину (температура, скорость вала, положение) при помощи управляющего устройства (обогреватель, контроллер мотора, линейный привод). Вот отличная статья по теории, что такое ПИД регулятор, как он работает и как его настроить. А я предлагаю свою библиотеку для работы с PID на Arduino.

ПИД регулятор выдаёт на выходе сигнал для плавного управления управляющим устройством (диммер, транзистор), если вам нужно реле – используйте библиотеку GyverRelay.

Алгоритм ПИД регулятора выглядит так, можете использовать его напрямую в скетче:

БИБЛИОТЕКА GYVERPID

GyverPID v3.0

Библиотека классического PID регулятора для Arduino

• Быстрая и лёгкая библиотека
• Время одного расчёта около 70 мкс
• Режим работы по величине или по её изменению (для интегрирующих процессов)
• На выбор целочисленная или float модель вычисления
• Возвращает результат по встроенному таймеру или в ручном режиме

Поддерживаемые платформы: все Arduino (используются стандартные Wiring-функции)

Версии

Версия 1.1 – убраны дефайны
Версия 1.2 – возвращены дефайны
Версия 1.3 – вычисления ускорены, библиотека облегчена
Версия 2.0 – логика работы чуть переосмыслена, код улучшен, упрощён и облегчён
Версия 2.1 – integral вынесен в public
Версия 2.2 – оптимизация вычислений
Версия 2.3 – добавлен режим PID_INTEGRAL_WINDOW
Версия 2.4 – реализация внесена в класс
Версия 3.0
– Добавлен режим оптимизации интегральной составляющей (см. доку)
– Добавлены автоматические калибровщики коэффициентов (см. примеры и доку)
Версия 3.1 – исправлен режиме ON_RATE, добавлено автоограничение инт. суммы

ДОКУМЕНТАЦИЯ

Документация

Логика работы

ПИД регулятор принимает на вход две величины:

• Входной сигнал input – сигнал с датчика: температура, скорость, положение, и т.д;
• Установку setpoint – величина, к которой регулятор будет стараться регулировать входной сигнал (температуру, скорость, положение…)

С ПИД регулятора выходит выходной (управляющий) сигнал output – безразмерная величина, которая подаётся на управляющее устройство. Это может быть транзистор с ШИМ сигналом, диммер, сервопривод, и т.д. Выходной сигнал должен влиять на входной сигнал: нагреватель нагревает объект с датчиком температуры, мотор крутится и даёт значения для датчика оборотов, и т.д.

Закон управления регулятора устанавливается при помощи коэффициентов Kp , Ki и Kd .

• Kp – пропорциональный коэффициент, выходная величина будет увеличиваться пропорционально разнице входного сигнала и установки.
• Ki – коэффициент интегрирующей составляющей, отвечает за накапливающуюся ошибку, позволяет сгладить пульсации и нивелировать маленькую ошибку.
• Kd – коэффициент дифференциальной составляющей, отвечает за скорость изменения величины, позволяет уменьшить раскачку системы.

Инициализация

Можно инициализировать объект тремя способами:

• GyverPID regulator; // инициализировать без настроек (всё по нулям, dt 100 мс)
• GyverPID regulator(kp, ki, kd); // инициализировать с коэффициентами. dt будет стандартно 100 мс
• GyverPID regulator(kp, ki, kd, dt); // инициализировать с коэффициентами и dt (в миллисекундах)

Режимы и настройки

Направление регулирования: зависит от того, в какую сторону направляется управляемая величина input при увеличении управляющего сигнала output . Например: охлаждение или нагрев, разгон или торможение, и т.д. По умолчанию стоит NORMAL – регулятор считает, что увеличение управляющего сигнала output увеличит входной сигнал input . Устанавливается командой

setDirection(dir); // dir – NORMAL или REVERSE

Режим работы: режим регулирования по ошибке входного сигнала ON_ERROR или по изменению входного сигнала ON_RATE . По умолчанию стоит ON_ERROR , его рекомендуется использовать в большинстве случаев, потому что большинство процессов – самоустанавливающиеся (температура нагревателя сама установится в своём максимуме, скорость мотора – тоже). Режим ON_RATE рекомендуется использовать в интегрирующих процессах, в которых выходная величина влияет на скорость изменения входной величины, например положение моторизированного слайдера, который не остановится при управляющем сигнале, отличном от нуля. Таким процессом будет проще управлять в режиме ON_RATE . Устанавливается командой

setMode(mode); // mode – ON_ERROR или ON_RATE

Подробнее про этот режим смотри в самом конце документации, в разделе оптимизации интегральной суммы.

Пределы выхода: ограничение значения выходного сигнала, по умолчанию: 0-255 (для 8 бит ШИМ). Может быть установлено 0-180 для угла сервопривода, и т.д. Устанавливается командой

setLimits(min, max); // установить пределы

Время итерации: время итерации можно изменить в процессе работы (не знаю, зачем, но возможность есть). Время устанавливается в миллисекундах и влияет на функцию getResultTimer() , которая с этим периодом делает новый расчёт управляющего сигнала. Также это время входит в расчёт управляющего сигнала (в И и Д составляющей). Устанавливается командой

setDt(dt); // установка времени итерации в мс

Установка/чтение параметров

Основные величины регулятора можно менять в любом месте программы любым удобным способом (кнопки, энкодер, передача через UART/GSM/WiFi, как угодно). Коэффициенты регулятора Kp , Ki и Kd можно устанавливать и читать напрямую как члены класса, например

Время итерации меняется при помощи метода setDt() (см. выше).

Величины регулятора (вход, установка, выход) также являются членами класса и к ним можно обратиться напрямую для чтения и записи:

Тип вычислений

Библиотека имеет режим целочисленных вычислений. Скорость вычислений особо не меняется, но код занимает меньше места и всё-таки должен выполняться быстрее. По умолчанию стоит режим чисел с плавающей точкой, в заголовочном файле библиотеки смотрите ключевое слово datatype , datatype будет float или int в зависимости от настройки: это некоторые переменные и функции. Настройка осуществляется дефайном перед подключением библиотеки:

Как работать с библиотекой?

Нужно скормить регулятору текущее значение величины в input , нужное значение в setpoint , провести расчёт при помощи getResult() или getResultTimer() , и после этого выходную величину output подать на управляющее устройство. Делать это нужно часто для быстрых процессов (стабилизация частоты оборотов шпинделя станка под нагрузкой: dt берём около 10-50 мс), и не очень часто для медленных процессов (удержание заданной температуры бойлера: dt можно взять пару секунд, процесс очень инерционный). Функция getResult() делает расчёт в каждый свой вызов и возвращает output , а getResultTimer() делает расчёт только при срабатывании встроенного таймера. То есть getResult() нужно вызывать по своему таймеру (для продвинутых пользователей), а getResultTimer() нужно вызывать как можно чаще, он посчитает только тогда, когда это будет нужно по своему таймеру. После расчёта можно подавать управляющий сигнал (выходную величину output ) на управляющее устройство. Смотрите пример!

Как настроить коэффициенты?

Подбор коэффициентов ПИД регулятора – индивидуальная задача, зависящая от конкретных условий и “железа”. Можно почитать статьи на эту тему: например эту, вот эту попроще, и вот эту посложнее. Первым делом нужно установить dt – об этом я писал выше. Маленький dt для быстрых процессов и побольше для медленных (инертных). Dt влияет на расчёты при неизменных коэффициентах, поэтому dt лучше не менять во время настройки, чтобы не пришлось пересчитывать все остальные коэффициенты. Диапазон коэффициентов: 0.01 – 100, т.е. довольно широк и зависит напрямую от инертности системы и выбранного времени dt. Коэффициенты должны быть положительные, противоположное направление регулирования задаётся в setDirection() .

В версии 3.0 появился автоматический тюнер коэффициентов, читай ниже.

Оптимизация интегральной суммы (экспериментально)

В реальной системе интегральная сумма может перенасыщаться и стать причиной неадекватного поведения регулятора, “заклинивания” его в крайних положениях. Библиотека предлагает несколько автоматических способов оптимизации интегральной суммы, остальную теорию можно почитать здесь.

С версии 3.1 в библиотеке работает автоматическое ограничение интегральной суммы по выходным лимитам регулятора.

Ручная оптимизация

В версии 2.1 интегральная сумма вынесена в публичный доступ как член класса, к ней можно обратиться как regulator.integral (где regulator – ваше имя объекта). Интегральная составляющая суммирует ошибку по времени, и при слишком сильном накоплении может приводить к перерегулированию (например для инерционных систем, таких как обогреватель). Для наблюдения за её состоянием можно прочитать integral , и при необходимости, например, ограничить её диапазон ( regulator.integral = constrain(regulator.integral, -500, 500); каждый раз после вызова getResult() ) или даже обнулить ( regulator.integral = 0; ) по условию.

Режим интегрального окна

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

Будьте внимательны, внутри библиотеки будет создан массив указанного размера и займёт память! Используйте только в том случае, если понимаете как это работает и для чего оно нужно, а также есть возможность наблюдать за графиком и делать выводы!

Режим оптимизации интегральной суммы

В версии 3.0 появился новый режим автоматической оптимизации интегральной суммы: она автоматически ограничивается так, чтобы выходной сигнал не превышал установленные в setLimits() пределы, то есть инт. сумма не будет бесконечно расти или уменьшаться. В то же время резкие скачки значения с датчика (вход регулятора) вблизи пределов могут приводить к обнулению интегральной суммы, поэтому входной сигнал рекомендуется фильтровать. Для активации режима оптимизации интегральной суммы нужно прописать в скетче дефайн #define PID_OPTIMIZED_I до подключения библиотеки.

Режим “пропорционально скорости”

Выше была описана смена режима работы при помощи setMode(mode); , в режиме ON_RATE регулятор лучше справляется с интегрирующими процессами (например позиция вала мотора), но иногда хорошо работает и с нагревателями, уменьшая переполнение интегральной суммы и перерегулирование в целом. В этом режиме коэффициенты ведут себя иначе и их оптимальные значения отличаются от обычного режима, логика ручной настройки также отличается: коэффициент Kp теперь работает только в паре с Ki и настраивать их нужно вместе, причём в некоторых процессах Kp не даёт никакого эффекта при нулевом Ki, либо работает неадекватно. Смотри пример simulation_linear, в котором симулируется идеальный интегрирующий процесс.

Реализация режима взята отсюда, там же есть теоретическое обоснование.

Автоматический тюнер коэффициентов 1

Тюнер тип 1

Автоматический калибровщик коэффициентов ПИД регулятора, метод “реле” http://auto-controls.blogspot.com/2009/10/pid-controllers-auto-tuning-relay.html. Данный тюнер лучше настраивает коэффициенты для удержания величины и парирования внешних возмущений.

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

Примечание: тюнер выдаёт коэффициенты для ПИ и ПИД регулятора. ПИ регулятор подразумевает, что коэффициент Д будет равен 0.

Как это работает?

1. Тюнер подаёт управляющий сигнал и ждёт стабилизации значения с датчика
2. Тюнер изменяет сигнал на некоторую величину (ступеньку)
3. Ждёт заданное время, затем меняет сигнал на ту же ступеньку, но в другую сторону
4. Начинается раскачка системы: при прохождении значения с датчика через значение стабилизации сигнал снова переключается
5. Производится анализ периода раскачки и её амплитуды, на основании этих данных вычисляются рекомендуемые коэффициенты

Как пользоваться библиотекой?

1. Инициализация и настройка

tuner.setParameters(направление, сигнал, ступенька, период, точность стабилизации, продолж. импульса, период итерации);

1. Направление :
   – NORMAL : увеличение выходного сигнала увеличивает сигнал с датчика (например обогреватель, мотор)
   – REVERSE : увеличение выходного сигнала уменьшает сигнал с датчика (например холодильник, тормоз)
2. Cигнал : базовый сигнал на управляющее устройство. Система будет ждать стабилизации по величине этого сигнала, и от него будет откладываться ступенька
3. Ступенька : величина, на которую будет изменяться сигнал в обе стороны от базового
4. Период : период опроса в ожидании стабилизации
5. Точность стабилизации : скорость изменения значения с датчика, ниже которой система будет считаться стабильной
6. Продолж. импульса : время в миллисекундах на первую раскачку
7. Период итерации : dt системы в мс, желательно должно совпадать с периодом ПИД регулятора

Пример: tuner.setParameters(NORMAL, 130, 40, 2000, 1, 2000, 30);

Калибруем нормальный процесс (увеличение сигнала увеличивает значение с датчика), базовый сигнал 130, ступенька 40 (в процессе калибровки сигнал будет меняться от 130-40 до 130+40, сразу думаем о том, чтобы это не превысило допустимый максимум), на этапе стабилизации хотим, чтобы система считалась стабильной при изменении сигнала с датчика менее, чем на 1 (условная величина датчика) за 2000 миллисекунд. На этапе первоначальной раскачки мы задали длину импульса 2000 мс, период работы всей системы – 30 мс.

2. Структура цикла
Библиотека сделана универсальной для любого датчика и управляющего устройства, цикл тюнинга организуется вот так:

3. Отладка и получение значений
3.1 Во время работы тюнера можно вызвать tuner.getAccuracy() – чем ближе его значение к 100, тем стабильнее на данный момент качается система и
тем вычисляемые коэффициенты будут более близки к идеальным
3.2 Для наблюдения за тюнером через Serial есть готовые методы:
– tuner.debugText() выводит текстовые данные (смотри скриншот в папке docs библиотеки)
– tuner.debugPlot() выводит данные для построения графика через плоттер Arduino IDE (смотри скриншот в папке docs библиотеки)
3.3 Чтобы получить коэффициенты внутри программы (без Serial) желательно задать условие if (tuner.getAccuracy() > 95) и при наступлении этого условия получить коэффициенты:

Смотрите примеры в examples/autotune

Автоматический тюнер коэффициентов 2

Тюнер тип 2

Автоматический калибровщик коэффициентов ПИД регулятора, метод Cohen-Coon https://pages.mtu.edu/

tbco/cm416/cctune.html. Данный тюнер лучше настраивает коэффициенты для переходного процесса, например разогрев с одной температуры до другой.

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

Примечание: тюнер выдаёт коэффициенты для ПИ и ПИД регулятора. ПИ регулятор подразумевает, что коэффициент Д будет равен 0.

Как это работает?

1. Тюнер подаёт стартовый управляющий сигнал и ждёт стабилизации значения с датчика
2. Тюнер запоминает минимальное значение и подаёт конечный сигнал, ждёт стабилизации
3. Тюнер запоминает максимальное значение, снова подаёт начальный сигнал и ждёт стабилизации
4. Тюнер снова подаёт конечный сигнал
5. Зная полное время процесса, тюнер измеряет сигнал в определённых точках и по специальным формулам считает коэффициенты

Как пользоваться библиотекой?

1. Инициализация и настройка

tuner.setParameters(направление, начальный сигнал, конечный сигнал, период, точность, время стабилизации, период итерации)

1. Направление :
   – NORMAL : увеличение выходного сигнала увеличивает сигнал с датчика (например обогреватель, мотор)
   – REVERSE : увеличение выходного сигнала уменьшает сигнал с датчика (например холодильник, тормоз)
2. Начальный сигнал : стартовый сигнал на управляющее устройство
3. Конечный сигнал : конечный сигнал на управляющее устройство
4. Период : период опроса в ожидании стабилизации
5. Точность стабилизации : скорость изменения значения с датчика, ниже которой система будет считаться стабильной
6. Период итерации : dt системы в мс, желательно должно совпадать с периодом ПИД регулятора

Пример: tuner.setParameters(NORMAL, 150, 200, 1000, 1, 50);

Калибруем нормальный процесс (увеличение сигнала увеличивает значение с датчика), начальный сигнал 150, конечный 200, на этапе стабилизации хотим, чтобы система считалась стабильной при изменении сигнала с датчика менее, чем на 1 (условная величина датчика) за 1000 миллисекунд. Период работы всей системы – 50 мс.

2. Структура цикла
Библиотека сделана универсальной для любого датчика и управляющего устройства, цикл тюнинга организуется вот так:

3. Отладка и получение значений
3.1 Во время работы тюнера можно вызвать tuner.getState() – вернёт номер текущего этапа работы. На 7-ом этапе можно забирать коэффициенты
3.2 Для наблюдения за тюнером через Serial есть готовые методы:
– tuner.debugText() выводит текстовые данные (смотри скриншот в папке docs библиотеки)
– tuner.debugPlot() выводит данные для построения графика через плоттер Arduino IDE (смотри скриншот в папке docs библиотеки)
3.3 Чтобы получить коэффициенты внутри программы (без Serial) желательно задать условие if (tuner.getState() == 7) и при наступлении этого условия получить коэффициенты:

Смотрите примеры в examples/autotune2