Диагностика отказов системы регулирования уровня в баке (стр. 9 из 18)
Входные образцы называют измерениями или вектором особенностей (свойств). Функция, выполняемая системой распознавания образцов, состоит в отображении входного вектора особенностей в один из классов решений. В диагностике отказов, этими классами решений являются различные типы отказов, возникающих в системе. Основное преимущество нейронных сетей – их способность разбиения пространства образцов с целью классификации. Следовательно, нейронные сети могут быть использованы как классификаторы для разделения образцов рассогласований и формирования сигналов тревог. Таким образом, они могут выявлять и изолировать отказы.
3. Диагностика отказов системы регулирования уровня жидкости в баке
3.1. Постановка задачи
Реализацию описанного выше метода диагностики отказов, основанного на моделях будем выполнять применительно к системе регулирования процессом экстракции (рисунок 3.1).
Рис.3.1. Система регулирования
Рассматриваемая система состоит из регулятора, исполнительного механизма, объекта и датчиков.
В качестве объекта автоматизации рассматриваем процесс жидкостной экстракции, осуществляемый в смесителях-отстойниках. Смесители-отстойники представим упрощенно в виде двух соединенных между собой баков, изображенных на рисунке 3.2.
В бак 1 поступает жидкость (вода) с известным расходом Q1. Уровень во втором баке необходимо регулировать в соответствии с заданием. Регулирование осуществляется за счет изменения вытекающего из второго бака потока Q3. Поток Q3 изменяется задвижкой, управляемой электроприводом на базе асинхронного двигателя. Величина уровня h2 измеряется датчиком.
Рис. 3.2. Система двух баков
Необходимо выполнить диагностику отказов элементов этой системы регулирования.
Следовательно, необходимо, на основе рассмотренного выше метода аналитической избыточности, спроектировать систему диагностики отказов, которая бы позволяла выявить и изолировать отказы датчиков, исполнительных механизмов и объекта управления.
3.2. Моделирование элементов системы и отказов
Для решения поставленной задачи в первую очередь необходимо создать модель системы, позволяющую имитировать поведение рассматриваемого процесса. Данная модель так же должна включать модели возможных отказов элементов системы регулирования.
3.2.1. Модель объекта управления
Рассматриваемая система баков (рисунок 3.2) может быть описана следующими уравнениями.
Уравнения материального баланса для баков 1 и 2 имеют вид:
, (3.1) 
, (3.2)где S1 = S2 = S = 0.049 м2 – площадь основания цилиндрических баков.
Расход через трубу, соединяющую баки, в соответствии с законом Торичелли определяется по формуле:
, (3.3)где
, 
- давление воды в 1 и 2 баках соответственно,γ = 9800 Н/м3 - удельный вес воды,
К1 – коэффициент пропускной способности трубы, соединяющей баки равный K1 = 0.05 м3/час.
Расход через вентиль:
, (3.4)где К2 – коэффициент пропускной способности вентиля. К2 регулируется задвижкой и зависит от ее положения – х:
К2(х) = 10∙Кmax∙х . (3.5)
Максимальная пропускная способность вентиля принимается равной
Кmax=0.1 м3/час.
По полученным уравнениям составим модель системы в переменных состояния:
, (3.6) 
, 
(3.7)В качестве переменных состояния будем рассматривать уровни жидкости в баках h1 и h2, входами будем считать расходы
и положение задвижки x: 
, 
. (3.8)Таким образом, получим:
(3.9)Данная модель является нелинейной.
3.2.2. Модель исполнительного механизма
В качестве исполнительного механизма рассматриваем трехфазный асинхронный двигатель, передаточная функция которого при частотном управлении имеет вид:
(3.10)где Ким = 1 – коэффициент усиления , T = TМ – электромеханическая постоянная времени двигателя, определяемая следующим образом:
=0,0396 сек, (3.11)где J=0.0081 кг ∙ м2 – момент инерции ротора,
ω0 = 2 ∙ 3,14 ∙ 50 = 314 рад/сек - синхронная частота,
SМ = 0.4371 – максимальное скольжение,
Ммах = 2,8025 Н ∙ м – максимальный электромагнитный момент.
3.2.3. Моделирование датчиков
Для измерения уровня жидкости в баке 2 используется датчик уровня. По причине того, что при его моделировании не будем учитывать динамические свойства датчика, его модель может быть описана в качестве передаточной функции вида:
Ws2(s) = 1. (3.12)
Для измерения положения задвижки так же используем датчик. Он описывается так же передаточной функцией вида:
Wsх(s) = 1. (3.13)
В соответствии с техническим заданием, уровень жидкости в первом баке датчиком не измеряется. Однако, в дальнейших исследованиях необходима информация об этом сигнале. Поэтому, для восстановления этой не измеряемой величины будем использовать виртуальный датчик – наблюдатель состояния.
Спроектируем наблюдатель состояния для системы двух баков. Для этого воспользуемся описанием системы баков в переменных состояния (3.9). Для линейной динамической модели:
(3.14)где
, 
, 
для воссоздания переменных системы на основе измерений входов и выходов используется наблюдатель состояния: 
. (3.15)Система (3.9) является нелинейной, поэтому для оценки ее состояния можно построить нелинейный наблюдатель следующего вида:
. (3.16)Для выбора коэффициентов обратной связи наблюдателя H необходимо выполнить линеаризацию в некоторой рабочей точке. Например, для точки
h1-h2 = 0.16357,м линеаризованная модель для (3.9) будет иметь следующий вид:
. (3.17)Коэффициенты матрицы Н выберем исходя из условия обеспечения устойчивости наблюдателя и с учетом, того, что наблюдатель должен обладать более высоким быстродействием, чем система.
.2.одель регулятора
3.2.5. Моделирование отказов
1. Отказ исполнительного механизма
Данный отказ будем моделировать как внезапное изменение коэффициента усиления исполнительного механизма в соответствии с уравнением:
КИМ = (1+δИМ)КИМ0, (3.20)
где δИМ - величина отказа, изменяемая в пределах {-1…1}. δИМ = 0 соответствует безотказному режиму работы.
2. Отказы датчиков
Данный тип отказов проявляется во внезапном возникновении отклонений показаний датчика от действительных значений, измеряемой величины. Моделирование осуществляется в соответствии с уравнением:
, (3.21)где yR– действительное значение измеряемой величины, y – значение величины, полученное датчиком, δS- величина отказа датчика.
Будем рассматривать следующие отказы датчиков: - датчик уровня, измеряющий h2; - датчик положения задвижки.
Величины отказов датчиков определяются следующим образом:
δS1 = {-1…1}, (3.22)
δS3 = α∙t. (3.23)
3. Отказы объекта управления
Будем рассматривать следующие отказы объекта управления:
1) отверстие в баке 1
Отказ моделируется как дополнительный расход, за счет которого уменьшается уровень в баке 1 в соответствии с уравнениями:
, 
(3.24)где Kfmax=0.5 м3/час – максимально возможный поток через отверстие.
Данный отказ может быть рассмотрен как внезапный и зарождающийся. В первом случае величина отказа dО1 измеряется в пределах {0…1}. dО1 = 0 соответствует безотказному режиму работы, dО1 = 1 – максимальной величине отверстия. Во втором случае величина отказа изменяется со скоростью αО1: