СПОСОБ ПОИСКА НЕИСПРАВНОГО БЛОКА В ДИСКРЕТНОЙ ДИНАМИЧЕСКОЙ СИСТЕМЕ
METHOD OF SEARCHING FOR FAULTY UNIT IN DISCRETE DYNAMIC SYSTEM

Изобретение относится к области контроля и диагностирования систем автоматического управления и их элементов.
Известен способ поиска неисправного блока в динамической системе (Патент на изобретение №2451319 от 20.05.2012 по заявке №2011129533/08(043690), МКИ 6 G05B 23/02, 2011).
Недостатком этого способа является то, что он обеспечивает определение дефектов только в непрерывной динамической системе.
Наиболее близким техническим решением (прототипом) является способ поиска неисправного блока в дискретной динамической системе (Патент на изобретение №2444774 от 10.03.2012 по заявке №2011101271/08(001575), MKH 6G05B 23/02, 2011).
Недостатком этого способа является то, что он предполагает интегрирование специальных тестовых сигналов с использованием экспоненциальной весовой функции и обеспечивает определение дефектов с невысокой различимостью, то есть обладает невысокой помехоустойчивостью.
Технической задачей, на решение которой направлено данное изобретение, является расширение функциональных возможностей способа путем применения рабочего диагностирования (без использования тестового воздействия), увеличение помехоустойчивости способа диагностирования дискретных систем автоматического управления путем улучшения различимости дефектов и уменьшение аппаратных затрат на вычисление весовой функции. Это достигается с помощью замены экспоненциальной весовой функции функцией, являющейся средним арифметическим модулей производных по времени сигналов системы с номинальными характеристиками контролируемой системы и моделей с пробными отклонениями.
jном
s
k
k
s
s
s
k
jном
ji
ji
jнорм
Поставленная задача достигается тем, что предварительно регистрируют реакцию заведомо исправной дискретной во времени системы f(t),j=1,…,k для N дискретных тактов диагностирования t∈[1, N] с дискретным постоянным шагом T на интервале наблюдения [0,T] (где Т=T·N) в k контрольных точках, и определяют интегральные оценки выходных сигналов
, j=1,…,k дискретной системы, для чего в момент подачи тестового или рабочего сигнала на вход дискретной системы с номинальными характеристиками одновременно начинают дискретное интегрирование сигналов системы управления с шагом T секунд в каждой из k контрольных точек с весовой функцией, равной среднему арифметическому значению модулей производных ее сигналов в контрольных точках, где усреднение производится по числу контрольных точек k. Для этого на первые входы k блоков перемножения подают сигналы системы, на вторые входы блоков перемножения подают среднее арифметическое значение модулей производных по времени сигналов, выходные сигналы k блоков перемножения подают на входы k блоков дискретного интегрирования с шагом T секунд, дискретное интегрирование завершают в момент времени Т, полученные в результате интегрирования оценки выходных сигналов f(d), j=1,…,k регистрируют, фиксируют число m рассматриваемых одиночных дефектов блоков, определяют интегральные оценки сигналов модели для каждой из k контрольных точек, полученных в результате пробных отклонений для m одиночных дефектов блоков, для чего поочередно в каждый блок дискретной динамической системы вводят пробное отклонение параметра дискретной передаточной функции и находят интегральные оценки выходных сигналов систем с пробными отклонениями при том же тестовом или рабочем сигнале x(t), полученные в результате дискретного интегрирования оценки выходных сигналов для каждой из k контрольных точек и каждого из m пробных отклонений
, j=1,…,k; i=1,…,m регистрируют, определяют отклонения интегральных оценок сигналов дискретной модели, полученных в результате пробных отклонений параметров разных структурных блоков ΔP(d)=P(d)-F(d), j=1,…,k; i=1,…,m, определяют нормированные значения отклонений интегральных оценок сигналов дискретной модели, полученных в результате пробных отклонений для одиночных дефектов из соотношения
j
j
j
jном
в момент начала контроля t=1 на вход контролируемой системы подают аналогичный тестовый или рабочий сигнал x(t), определяют интегральные оценки сигналов контролируемой дискретной системы для k контрольных точек F(d), j=1,…,k, полученные значения регистрируют, определяют отклонения интегральных оценок сигналов контролируемой дискретной системы для, k контрольных точек от номинальных значений ΔF(d)=F(d)-F(d), j=1,…,k определяют нормированные значения отклонений интегральных оценок сигналов контролируемой дискретной системы из соотношения
определяют диагностические признаки из соотношения
по минимуму значения диагностического признака определяют порядковый номер дефектного блока.
Сущность предлагаемого способа заключается в следующем.
k
Способ основан на использовании пробных отклонений параметров модели дискретной динамической системы. Для получения диагностических признаков динамических элементов используются интегральные оценки на временном интервале T в k контрольных точках
Весовая функция в формуле (4) в виде среднего значения модулей производных сигналов в контрольных точках несет информацию о важности момента времени с точки зрения скорости изменения сигналов во всех контрольных точках. Чем больше средняя скорость изменения сигналов, тем с большим весом интегрируется выходной сигнал.
Используя векторную интерпретацию выражения (3), запишем его в следующем виде
i
где φ(d) - угол между нормированным вектором (вектором единичной длины) отклонений интегральных оценок сигналов дискретного объекта
и нормированным вектором (единичной длины) отклонений интегральных оценок сигналов дискретной модели
полученных в результате пробного отклонения i-ого параметра соответствующего структурного блока.
Таким образом, нормированный диагностический признак (3) представляет собой значение квадрата синуса угла, образованного в k-мерном пространстве (где k - число контрольных точек) нормированными векторами интегральных оценок пробных отклонений сигналов дискретной модели и отклонений интегральных оценок сигналов дискретного объекта диагностирования.
Пробное отклонение параметра соответствующего структурного блока, минимизирующее значение диагностического признака (3), указывает на наличие дефекта блока. Область возможных значений диагностического признака лежит в интервале [0, 1].
Таким образом, предлагаемый способ поиска неисправного блока сводится к выполнению следующих операций:
s
1. В качестве дискретной динамической системы рассматривают систему, например с дискретной интерполяцией нулевого порядка, с шагом дискретизации T, состоящую из произвольно соединенных динамических блоков, с количеством рассматриваемых одиночных дефектов блоков m.
К
ПП
2. Предварительно определяют время контроля Т>Т, где
ПП
Т время переходного процесса дискретной системы. Время переходного процесса оценивают для номинальных значений параметров динамической системы.
3. Фиксируют число контрольных точек k.
4. Одновременно подают входной сигнал x(t) (единичный ступенчатый, линейно возрастающий, прямоугольный импульсный и т.д.) на вход системы управления с номинальными характеристиками, на вход контролируемой системы, на входы m моделей с номинальными параметрами, в каждую из которых введены пробные отклонения параметров одного блока так, что в i-ю систему введены пробные отклонения в i-й блок. Принципиальных ограничений на вид входного воздействия предлагаемый способ не предусматривает.
jнорм
s
k
k
s
j
ji
s
k
s
jном
j
ji
5. Одновременно регистрируют реакцию системы f(t), j=1,…,k на интервале t∈[1,N] с дискретным шагом T секунд на интервале наблюдения [0, T] (где Т=T·N) в k контрольных точках и определяют дискретные интегральные оценки выходных сигналов
, j=1,…,k системы с номинальными характеристиками, контролируемой системы F(d), j=1,…,k моделей с пробными отклонениями в i-м блоке P(d), j=1,…,k; i=1,…,m (формула 4). Для этого в момент подачи входного сигнала одновременно начинают дискретное интегрирование сигналов системы управления с шагом Ts секунд в каждой из k контрольных точек системы с номинальными характеристиками, контролируемой системы, моделей с пробными отклонениями параметров блоков с весовой функцией, равной среднему арифметическому значению модулей производных сигналов в контрольных точках, где усреднение производится по числу контрольных точек, для чего выходные сигналы каждой системы подают на первые входы k блоков перемножения, на вторые входы блоков перемножения подают среднее арифметическое значение модулей производных сигналов системы в контрольных точках, где усреднение производится по числу контрольных точек выходных сигналов системы, выходные сигналы k блоков- перемножения подают на входы k блоков дискретного интегрирования с шагом T секунд, дискретное интегрирование завершают в момент времени Т, полученные в результате дискретного интегрирования с шагом T секунд оценки выходных сигналов F(d) , j=1,…,k, F(d), j=1,…,k, P(d), j=1,…,k; i=1,…,m регистрируют.
ji
ji
jном
6. Определяют отклонения интегральных оценок сигналов дискретной модели, полученные в результате пробных отклонений параметров одного структурного блока ΔP(d)=P(d)-F(d), j=1,…,k; i=1…,m.
7. Определяют нормированные значения отклонений интегральных оценок сигналов дискретной модели, полученные в результате пробных отклонений параметров одного блока по формуле:
j
j
jном
8. Определяют отклонения интегральных оценок сигналов контролируемой дискретной системы для k контрольных точек от номинальных значений ΔF(d)=F(d)-F(d), j=1,…,k,
9. Вычисляют нормированные значения отклонений интегральных оценок сигналов контролируемой дискретной системы по формуле:
10. Вычисляют диагностические признаки наличия неисправного структурного блока по формуле (3).
11. По минимуму значения диагностического признака определяют дефектный блок.
Рассмотрим реализацию предлагаемого способа поиска дефекта для дискретной системы, структурная схема которой представлена на фиг.1.
Дискретные передаточные функции блоков:
;
,
1
1
2
2
3
3
где номинальные значения параметров: k=5; Z=0.98; k=0.09516; Q=0.9048; k=0.0198; Q=0.9802.
k
При моделировании в качестве входного сигнала будем использовать псевдослучайный сигнал (при моделировании использовался блок Band-Limited White Noise в среде Matlab). Время контроля выберем Т равным 10 с.
Величину пробных отклонений параметров модели выбираем равной 10%.
1
1
2
3
Моделирование процессов поиска дефектов в блоке 1 (в виде уменьшения параметра k на 20%) приводит к вычислению диагностических признаков по формуле (3): J=0, J=0.184, J=0.2361.
2
1
Различимость дефекта: ΔJ-J=0.184.
6
Для сравнения приведем диагностические признаки наличия неисправного блока с использованием экспоненциального веса при одном параметре интегрирования α=0.5 (Патент на изобретение №2444774 от 10.03.2012 по заявке №2011101271/08(001575), МКИ G05В 23/02, 2011):
1
2
3
J=0, J=0.3587, J=0.1605. Различимость дефекта
3
1
ΔJ=J-J=0.1605.
Приведенные результаты показывают, что фактическая различимость нахождения дефектов этим способом выше, следовательно, выше будет и помехоустойчивость способа.
2
Моделирование процессов поиска дефектов в блоке 2 (в виде уменьшения параметра k на 20%) для данного объекта диагностирования с использованием дифференциального веса и при таком же входном сигнале дает следующие значения диагностических признаков:
1
2
3
J=0.173, J=0, J=0.6335.
1
2
Различимость дефекта: ΔJ=J-J=0.173.
Для сравнения приведем диагностические признаки наличия неисправного блока с использованием экспоненциального веса при одном
1
2
3
параметре интегрирования α=0.5: J=0.3557, J=0, J=0.6732.
1
2
Различимость дефекта: ΔJ=J-J=0.3557.
3
Моделирование процессов поиска дефектов в блоке 3 (в виде уменьшения параметра k на 20%) для данного объекта диагностирования при тех же состояниях дает следующие значения:
1
2
3
J=0.2457, J=0.6634, J=0.
1
3
Различимость дефекта: ΔJ=J-J=0.2457.
Для сравнения приведем диагностические признаки наличия неисправного блока при одном параметре интегрирования α=0.5:
1
2
3
J=0.1652, J,=0.668, J=o.
1
3
Различимость дефекта ΔJ=J-J=0.1652.
Минимальное значение диагностического признака во всех случаях правильно указывает на дефектный блок, а данный способ в двух случаях из трех улучшает фактическую различимость дефектов, следовательно, увеличивает помехоустойчивость диагностирования. Кроме того, заявляемый способ позволяет осуществлять диагностирование в условиях реального функционирования объекта диагностирования (рабочее диагностирование).