Моделирование системы управления с регулятором дробного порядка и исследование её устойчивости

Серия 4. Химическое машиностроение и инженерная экология шеным продуктом, как например выделяемым дрожжами Phaffia rhodozyma (они же -Xanthophyllomyces dendrorhous) каратиноидом астаксантином, который характеризуется насыщенным красным цветом, - целесообразно дополнить схему ещё одним фотоэлементом, снабжённым цветным светофильтром, пропускающим только свет с частотой, соответствующей продукту. В результате представляется возможным оценивать как концентрацию продукта, так и клеток продуцента в режиме on-line. Разумеется, погрешность такого измерения существенно выше, чем у лабораторного анализа, но с учётом его корректировок может быть полезным инструментом.

Прошедший через культу рал ьную среду Поток света и стенку поток света

Рисунок 4. Схема управления освещённостью в лабораторном фотобиореакторе

Работа выполнена в рамках выполнения государственного задания высшим учебным заведениям в части проведения научно-исследовательских работ, проект № 7.6102.2011 по теме: «Разработка систем энергоэффективного управления процессами биосинтеза».

Литература

1. Захаров З.В., Герман Л.С., Петрищева O.A., Жарко М.Ю. Культивирование дрожжей Phaffia Rhodozyma при постоянном и периодическом освещении // Известия Московского государственного технического университета МАМИ. 2012. Т. 4. № 2.

2. Мальцевская Н.В., Бирюков В.В. Влияние преры-вистого освещения на процесс роста фо-тотрофного микроорганизма Chlorella Sp. // Биотехнология, 2011. № 1. c. 47.

3. Зубов Д.В., Строков С.С., Ефремов Д. А. Автоматизация лабораторного стенда для исследования процессов культивирования фотозависимых микроорганизмов // Труды XIX Всероссийской научно-методической конференции Телематика'2012, т.1, с.68-69

Моделирование системы управления с регулятором дробного порядка и исследование её устойчивости

к.т.н. доц. Зубов Д.В., аспирант Студеникин Г.И.

Университет машиностроения zubov@msuie.ru

Аннотация. Рассмотрена возможность создания систем управления с регуляторами дробного порядка, их преимущества и возникающие при реализации слож-

Серия 4. Химическое машиностроение и инженерная экология

ности.

Ключевые слова: дифферинтеграл Римана — Лиувилля, устойчивость систем автоматизации.

Задачами исследования являются:

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

- анализ устойчивости системы при помощи метода расширенных частотных характеристик;

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

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

обеспечения адекватности математических моделей исследуемым объектам. Динамические системы, как объект моделирования, традиционно изучались путем использования классического математического анализа, в частности аппарата интегро-дифференциальных уравнений в обыкновенных и частных производных. Классический анализ предполагает, что интегралы и производные имеют порядки, выражаемые целыми числами.

Между тем, поведение целого ряда объектов и процессов не соответствует в полной мере используемым математическим моделям, необходимо разрабатывать и использовать уточненные модели, в том числе с использованием производной и интеграла нецелых порядков (нецелые порядки могут быть дробными, иррациональными и комплексными числами). Хотя история возникновения и развития дробного исчисления насчитывает уже более трех столетий, прикладное применение усложнялось сложностью вычислений, но при появлении компьютерных программ для математических расчетов данная проблема отпадает.

Важное место в общей теории обработки сигналов, моделировании и автоматическом управлении занимают операционные методы анализа. Преобразование Лапласа является одним из наиболее распространенных операционных методов, который позволяет анализировать динамические системы в переходном режиме. В рамках этого преобразования рассматриваются два пространства: пространство оригиналов (сигнальное пространство) и пространство изображений сигналов (изображение по Лапласу). Математической моделью переходного процесса динамической системы в первом пространстве являются интегро-дифференциальные уравнения. Во втором (преобразованном) пространстве математической моделью переходного процесса являются алгебраические уравнения.

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

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

апБапу{1) + Б^уЦ) + • • • + а0 Ба° у{1) =

= ЪпБРпх

где: а;, Ь - коэффициенты уравнения, а;, Р] - дробные порядки дифференциальных операторов, у(1;) - функция выхода динамической системы (функция состояния), х(1;)- функция входа динамической системы (функция управления). В случае нулевых начальных условий передаточная характеристика динамической системы в области преобразования по Лапласу принимает вид:

Серия 4. Химическое машиностроение и инженерная экология

^

а„ря" + а^р""-1 +••• + а0ра°

Для вычисления производных и интегралов дробных порядков в системах управления широко используются различные апроксимационные зависимости, базирующиеся на классической теории дробномерных дифференциальных операторов, построенных на основании частотных методов теории автоматического управления и представляющих собой приближенную динамическую модель звена с дробномерным дифференциальным оператором. Такой подход является приближенными и не имеет строгого математического обоснования. Поэтому для вычисления производных и интегралов дробных порядков чаще используется формула Грюнвальда-Летникова [1], в соответствии с которой дробная производная р аг)1(^) определяется следующим образом:

N-1 1 N-1

1 N -1 1 N -1

рЧ (*) = Шп—X ^ = Шп—X D(a, г)п1 (* - т) =

,=о Т Т 1=°

1 N=1, ,4! Г(« +1) /

тТ Г

Также может использоваться дифферинтеграл Римана — Лиувилля:

Iа (/(х)) = \/(*)(х - *)-1 а*. г(а):

Рисунок 1. Кривые равной колебательности

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

Wo (Р) = ,

Т  р +1

с параметрами объекта равными: к° = 135, т = 15, ^ = 180 .

Серия 4. Химическое машиностроение и инженерная экология

В качестве регулятора использовались: обычный ПИ-регулятор, регуляторы дробного

порядка ПИ0'5 , ПИио" , ПИ'

0,65

т-0'75

ПИ

0,95

с соответствующими передаточными функциями:

С л \

R(P) = ~kv

1 + -

1

Tu P0

где а - дробный показатель степени (0,5; 0,65; 0,75; 0,95);

кр - настройки пропорциональной составляющей регулятора;

Ти - настройка интегральной составляющей регулятора.

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

сти для систем с регуляторами ПИ0'5 , ПИ0'65 , ПИ0'75 , ПИ0'95 (см. рисунок 1) и соответствующие кривые разгона (рисунок 2), найдены оптимальные настройки ПИ-регулятора для различных степеней дробной производной по модульному интегральному критерию.

г \ 1 1 1 1' 1 1 1 1 - ь k и0*75 1 1 1 1 1

| 0,65 /0,5

1 Л 7 " — - я -------------

V ¡г V ' т v/'- • \ V/ 1 : L: - * од/ 1 1 1 1 1 1 1 1 1

Время

Рисунок 2. Кривые разгона

В результате получено, что при использовании ПИ-регуляторов дробного порядка с объектом целого порялка можно получить больший запас устойчивости, чем с обычным ПИ-регулятором, но при этом неизбежно появление статической ошибки.

Литература

1. ВасильевВ.В., СимакЛ.А. Дробное исчисление и аппроксимационные методы в моделировании динамических систем. -Киев, HAH Украины, 2008. — 256 с.

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

к.т.н. доц. Поляков А.Н., Смирнова A.C., Щелканова О.Н., Киселева A.C.

Университет машиностроения nauka@msuie.ru

Аннотация. Исследованы протекающие при очистке сточных вод процессы,