способ поиска неисправного блока в дискретной динамической системе

Формула изобретения

1. Способ поиска неисправного блока в дискретной динамической системе, основанный на том, что фиксируют число k контрольных точек системы, фиксируют число m динамических элементов, входящих в состав системы, определяют время контроля Т_К>Т _ПП, используют входной сигнал x(t) на интервале t [0, T_K], регистрируют реакцию объекта диагностирования и реакцию заведомо исправной системы F_jном(t), j=1, ,k для N дискретных тактов диагностирования с дискретным постоянным шагом T_s на интервале наблюдения [0, T _k] (где T_k=T_s·N) в k контрольных точках, фиксируют число различных пробных отклонений m, определяют интегральные оценки сигналов модели для каждой из k контрольных точек, полученных в результате пробных отклонений для m одиночных дефектов блоков, для чего поочередно в каждый блок дискретной динамической системы вводят пробное отклонение параметра дискретной передаточной функции и находят интегральные оценки выходных сигналов системы для входного сигнала x(t), полученные в результате дискретного интегрирования оценки выходных сигналов для каждой из k контрольных точек и каждого из m пробных отклонений регистрируют, определяют отклонения интегральных оценок сигналов дискретной модели, полученные в результате пробных отклонений параметров, определяют нормированные значения отклонений интегральных оценок сигналов дискретной модели, полученные в результате пробных отклонений параметров, определяют отклонения интегральных оценок сигналов контролируемой дискретной системы для k контрольных точек от номинальных значений, вычисляют нормированные значения отклонений интегральных оценок сигналов контролируемой дискретной системы, по минимуму диагностического признака определяют дефектный блок дискретной системы, отличающийся тем, что одновременно подают тестовый или рабочий сигнал x(t) на вход системы с номинальными характеристиками, на вход контролируемой системы, на входы m моделей с номинальными характеристиками, в каждую из которых введены пробные отклонения параметров одного блока так, что в i-ю систему вводят пробные отклонения в i-й блок, в качестве динамических характеристик системы используют интегральные оценки, полученные для весовой функции, равной среднему арифметическому модулей производных по времени от выходных сигналов системы в различных контрольных точках, из соотношения , определяют интегральные оценки выходных сигналов F _jном(d), j=1, ,k системы с номинальными характеристиками, для чего в момент подачи тестового или рабочего сигнала на вход системы с номинальными характеристиками одновременно начинают интегрирование сигналов системы в каждой из k контрольных точек для весовой функции, путем подачи на первые входы k блоков перемножения сигналов системы на вторые входы блоков перемножения подают среднее арифметическое модулей производных по времени от выходных сигналов системы с номинальными характеристиками, выходные сигналы k блоков перемножения подают на входы k блоков интегрирования, интегрирование завершают в момент времени Т_К, полученные в результате интегрирования оценки выходных сигналов F_jном(d), j=1, ,k регистрируют, аналогично определяют интегральные оценки сигналов m моделей для каждой из k контрольных точек, полученные в результате пробных отклонений параметров каждого из m блоков, полученные в результате интегрирования оценки выходных сигналов для каждой из k контрольных точек, каждого из m пробных отклонений P_ji(d), j=1, ,k; i=1, ,m регистрируют, определяют отклонения интегральных оценок сигналов дискретной модели, полученные в результате пробных отклонений параметров одного структурного блока P_ji(d)=P_ji(d)-F_jном(d), j=1, ,k; i=1 ,m., определяют нормированные значения отклонений интегральных оценок сигналов дискретной модели, полученные в результате пробных отклонений параметров одного блока по формуле



определяют отклонения интегральных оценок сигналов контролируемой дискретной системы для k контрольных точек от номинальных значений F_j(d)=F_j(d)-F_jном(d), j=1, ,k.

вычисляют нормированные значения отклонений интегральных оценок сигналов контролируемой дискретной системы по формуле:

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

, i=1, ,m.

Описание изобретения к патенту

Изобретение относится к области контроля и диагностирования систем автоматического управления и их элементов.

Известен способ поиска неисправного блока в динамической системе (Патент на изобретение № 2451319 от 20.05.2012 по заявке № 2011129533/08(043690), МКИ 6 G05B 23/02, 2011).

Недостатком этого способа является то, что он обеспечивает определение дефектов только в непрерывной динамической системе.

Наиболее близким техническим решением (прототипом) является способ поиска неисправного блока в дискретной динамической системе (Патент на изобретение № 2444774 от 10.03.2012 по заявке № 2011101271/08(001575), MKH 6G05B 23/02, 2011).

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

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

Поставленная задача достигается тем, что предварительно регистрируют реакцию заведомо исправной дискретной во времени системы f_jном (t),j=1, ,k для N дискретных тактов диагностирования t [1, N] с дискретным постоянным шагом T_s на интервале наблюдения [0,T_k] (где Т_k=T_s ·N) в k контрольных точках, и определяют интегральные оценки выходных сигналов , j=1, ,k дискретной системы, для чего в момент подачи тестового или рабочего сигнала на вход дискретной системы с номинальными характеристиками одновременно начинают дискретное интегрирование сигналов системы управления с шагом T_s секунд в каждой из k контрольных точек с весовой функцией, равной среднему арифметическому значению модулей производных ее сигналов в контрольных точках, где усреднение производится по числу контрольных точек k. Для этого на первые входы k блоков перемножения подают сигналы системы, на вторые входы блоков перемножения подают среднее арифметическое значение модулей производных по времени сигналов, выходные сигналы k блоков перемножения подают на входы k блоков дискретного интегрирования с шагом T_s секунд, дискретное интегрирование завершают в момент времени Т_k, полученные в результате интегрирования оценки выходных сигналов f_jном(d), j=1, ,k регистрируют, фиксируют число m рассматриваемых одиночных дефектов блоков, определяют интегральные оценки сигналов модели для каждой из k контрольных точек, полученных в результате пробных отклонений для m одиночных дефектов блоков, для чего поочередно в каждый блок дискретной динамической системы вводят пробное отклонение параметра дискретной передаточной функции и находят интегральные оценки выходных сигналов систем с пробными отклонениями при том же тестовом или рабочем сигнале x(t), полученные в результате дискретного интегрирования оценки выходных сигналов для каждой из k контрольных точек и каждого из m пробных отклонений , j=1, ,k; i=1, ,m регистрируют, определяют отклонения интегральных оценок сигналов дискретной модели, полученных в результате пробных отклонений параметров разных структурных блоков P_ji(d)=P_ji(d)-F_jнорм(d), j=1, ,k; i=1, ,m, определяют нормированные значения отклонений интегральных оценок сигналов дискретной модели, полученных в результате пробных отклонений для одиночных дефектов из соотношения

в момент начала контроля t=1 на вход контролируемой системы подают аналогичный тестовый или рабочий сигнал x(t), определяют интегральные оценки сигналов контролируемой дискретной системы для k контрольных точек F_j(d), j=1, ,k, полученные значения регистрируют, определяют отклонения интегральных оценок сигналов контролируемой дискретной системы для, k контрольных точек от номинальных значений F_j(d)=F_j(d)-F_jном(d), j=1, ,k определяют нормированные значения отклонений интегральных оценок сигналов контролируемой дискретной системы из соотношения

определяют диагностические признаки из соотношения

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

Сущность предлагаемого способа заключается в следующем.

Способ основан на использовании пробных отклонений параметров модели дискретной динамической системы. Для получения диагностических признаков динамических элементов используются интегральные оценки на временном интервале T_k в k контрольных точках

Весовая функция в формуле (4) в виде среднего значения модулей производных сигналов в контрольных точках несет информацию о важности момента времени с точки зрения скорости изменения сигналов во всех контрольных точках. Чем больше средняя скорость изменения сигналов, тем с большим весом интегрируется выходной сигнал.

Используя векторную интерпретацию выражения (3), запишем его в следующем виде

где _i(d) - угол между нормированным вектором (вектором единичной длины) отклонений интегральных оценок сигналов дискретного объекта и нормированным вектором (единичной длины) отклонений интегральных оценок сигналов дискретной модели полученных в результате пробного отклонения i-ого параметра соответствующего структурного блока.

Таким образом, нормированный диагностический признак (3) представляет собой значение квадрата синуса угла, образованного в k-мерном пространстве (где k - число контрольных точек) нормированными векторами интегральных оценок пробных отклонений сигналов дискретной модели и отклонений интегральных оценок сигналов дискретного объекта диагностирования.

Пробное отклонение параметра соответствующего структурного блока, минимизирующее значение диагностического признака (3), указывает на наличие дефекта блока. Область возможных значений диагностического признака лежит в интервале [0, 1].

Таким образом, предлагаемый способ поиска неисправного блока сводится к выполнению следующих операций:

1. В качестве дискретной динамической системы рассматривают систему, например с дискретной интерполяцией нулевого порядка, с шагом дискретизации T_s, состоящую из произвольно соединенных динамических блоков, с количеством рассматриваемых одиночных дефектов блоков m.

2. Предварительно определяют время контроля Т_К>Т_ПП, где

Т_ПП время переходного процесса дискретной системы. Время переходного процесса оценивают для номинальных значений параметров динамической системы.

3. Фиксируют число контрольных точек k.

4. Одновременно подают входной сигнал x(t) (единичный ступенчатый, линейно возрастающий, прямоугольный импульсный и т.д.) на вход системы управления с номинальными характеристиками, на вход контролируемой системы, на входы m моделей с номинальными параметрами, в каждую из которых введены пробные отклонения параметров одного блока так, что в i-ю систему введены пробные отклонения в i-й блок. Принципиальных ограничений на вид входного воздействия предлагаемый способ не предусматривает.

5. Одновременно регистрируют реакцию системы f_jнорм(t), j=1, ,k на интервале t [1,N] с дискретным шагом T_s секунд на интервале наблюдения [0, T_k] (где Т_k=T_s ·N) в k контрольных точках и определяют дискретные интегральные оценки выходных сигналов , j=1, ,k системы с номинальными характеристиками, контролируемой системы F_j(d), j=1, ,k моделей с пробными отклонениями в i-м блоке P_ji (d), j=1, ,k; i=1, ,m (формула 4). Для этого в момент подачи входного сигнала одновременно начинают дискретное интегрирование сигналов системы управления с шагом Ts секунд в каждой из k контрольных точек системы с номинальными характеристиками, контролируемой системы, моделей с пробными отклонениями параметров блоков с весовой функцией, равной среднему арифметическому значению модулей производных сигналов в контрольных точках, где усреднение производится по числу контрольных точек, для чего выходные сигналы каждой системы подают на первые входы k блоков перемножения, на вторые входы блоков перемножения подают среднее арифметическое значение модулей производных сигналов системы в контрольных точках, где усреднение производится по числу контрольных точек выходных сигналов системы, выходные сигналы k блоков- перемножения подают на входы k блоков дискретного интегрирования с шагом T_s секунд, дискретное интегрирование завершают в момент времени Т_k, полученные в результате дискретного интегрирования с шагом T_s секунд оценки выходных сигналов F_jном(d) , j=1, ,k, F_j(d), j=1, ,k, P_ji(d), j=1, ,k; i=1, ,m регистрируют.

6. Определяют отклонения интегральных оценок сигналов дискретной модели, полученные в результате пробных отклонений параметров одного структурного блока P_ji(d)=P_ji(d)-F_jном(d), j=1, ,k; i=1 ,m.

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

8. Определяют отклонения интегральных оценок сигналов контролируемой дискретной системы для k контрольных точек от номинальных значений F_j(d)=F_j(d)-F_jном(d), j=1, ,k,

9. Вычисляют нормированные значения отклонений интегральных оценок сигналов контролируемой дискретной системы по формуле:

10. Вычисляют диагностические признаки наличия неисправного структурного блока по формуле (3).

11. По минимуму значения диагностического признака определяют дефектный блок.

Рассмотрим реализацию предлагаемого способа поиска дефекта для дискретной системы, структурная схема которой представлена на фиг.1.

Дискретные передаточные функции блоков:

; ,

где номинальные значения параметров: k ₁=5; Z₁=0.98; k₂=0.09516; Q₂ =0.9048; k₃=0.0198; Q₃=0.9802.

При моделировании в качестве входного сигнала будем использовать псевдослучайный сигнал (при моделировании использовался блок Band-Limited White Noise в среде Matlab). Время контроля выберем Т_k равным 10 с.

Величину пробных отклонений параметров модели выбираем равной 10%.

Моделирование процессов поиска дефектов в блоке 1 (в виде уменьшения параметра k₁ на 20%) приводит к вычислению диагностических признаков по формуле (3): J₁=0, J₂=0.184, J₃ =0.2361.

Различимость дефекта: J₂-J₁=0.184.

Для сравнения приведем диагностические признаки наличия неисправного блока с использованием экспоненциального веса при одном параметре интегрирования =0.5 (Патент на изобретение № 2444774 от 10.03.2012 по заявке № 2011101271/08(001575), МКИ⁶ G05В 23/02, 2011):

J₁=0, J₂=0.3587, J₃ =0.1605. Различимость дефекта

J=J₃-J₁=0.1605.

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

Моделирование процессов поиска дефектов в блоке 2 (в виде уменьшения параметра k₂ на 20%) для данного объекта диагностирования с использованием дифференциального веса и при таком же входном сигнале дает следующие значения диагностических признаков:

J₁=0.173, J₂=0, J₃=0.6335.

Различимость дефекта: J=J₁-J₂=0.173.

Для сравнения приведем диагностические признаки наличия неисправного блока с использованием экспоненциального веса при одном

параметре интегрирования =0.5: J₁=0.3557, J₂=0, J₃ =0.6732.

Различимость дефекта: J=J₁-J₂=0.3557.

Моделирование процессов поиска дефектов в блоке 3 (в виде уменьшения параметра k₃ на 20%) для данного объекта диагностирования при тех же состояниях дает следующие значения:

J ₁=0.2457, J₂=0.6634, J₃=0.

Различимость дефекта: J=J₁-J₃=0.2457.

Для сравнения приведем диагностические признаки наличия неисправного блока при одном параметре интегрирования =0.5:

J₁=0.1652, J₂ ,=0.668, J₃=o.

Различимость дефекта J=J₁-J₃=0.1652.

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