СПОСОБ ПОИСКА НЕИСПРАВНОСТЕЙ ДИНАМИЧЕСКОГО БЛОКА В НЕПРЕРЫВНОЙ СИСТЕМЕ Российский патент 2012 года по МПК G05B23/02 

Описание патента на изобретение RU2464616C1

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

Известен способ поиска неисправного блока в динамической системе (Положительное решение от 12.07.2010 г. о выдаче патента на изобретение по заявке №2009123999/08(033242), МКИ6 G05В 23/02, 2010).

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

Наиболее близким техническим решением (прототипом) является способ поиска неисправностей динамического блока в непрерывной системе (Патент РФ на изобретение №2429518 по заявке №2010128421/08(040385), МКИ6 G05В 23/02, 2011).

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

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

Поставленная задача достигается тем, что фиксируют число m параметров всех динамических блоков контролируемой системы, регистрируют реакцию заведомо исправной системы , j=1, 2, …, k на интервале в k контрольных точках, и определяют интегральные оценки выходных сигналов Fj ном(d), j=1, …, k системы, для чего в момент подачи тестового или рабочего сигнала на вход системы с номинальными характеристиками одновременно начинают интегрирование сигналов этой системы для каждой из k контрольных точек с весовой функцией, равной среднему арифметическому значению модулей производных ее сигналов в контрольных точках, где усреднение производится по числу контрольных точек. Для этого на первые входы k блоков перемножения подают сигналы системы, на вторые входы блоков перемножения подают среднее арифметическое значение модулей производных по времени сигналов, выходные сигналы k блоков перемножения подают на входы k блоков интегрирования, интегрирование завершают в момент времени Tк, полученные в результате интегрирования оценки выходных сигналов Fj нoм(d), j=1, …, k регистрируют, одновременно определяют интегральные оценки сигналов m моделей для каждой из k контрольных точек, полученные в результате пробных отклонений поочередно

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

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

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

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

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

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

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

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

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

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

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

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

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

5. Одновременно регистрируют реакцию системы с номинальными характеристиками , реакцию контролируемой системы , реакции моделей с пробным отклонением в i-ом параметре в k контрольных точках j=1, 2, …, k на интервале .

6. Одновременно определяют интегральные оценки выходных сигналов Fj ном(d), j=1, …, k системы с номинальными характеристиками, контролируемой системы Fj(d), j=1, …, k, моделей с пробными отклонениями в i-ом параметре Pji(d), j=1, …, k; i=1, …, m (формула 4). Для этого в момент подачи входного сигнала одновременно начинают интегрирование сигналов в каждой из k контрольных точек системы с номинальными характеристиками, контролируемой системы, моделей с пробными отклонениями параметров блоков с весовой функцией, равной среднему арифметическому значению модулей производных сигналов в контрольных точках, где усреднение производится по числу контрольных точек, для чего выходные сигналы каждой системы подают на первые входы k блоков перемножения, на вторые входы блоков перемножения подают среднее арифметическое значение модулей производных сигналов системы в контрольных точках, выходные сигналы k блоков перемножения подают на входы k блоков интегрирования, интегрирование завершают в момент времени Tк, полученные в результате интегрирования оценки выходных сигналов Fj ном(d), j=1, …, k, Fj(d), j=1, …, k, Pji(d), j=1, …, k; i=1, …, m регистрируют.

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

ΔPji(d)=Pji(d)-Fj ном(d), j=1, …, k; i=1, …, m.

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

, j=1, …, k; i=1, …, m.

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

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

, j=1, …, k.

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

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

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

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

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

; ; ,

номинальные значения параметров: T1=5 с (J1); k1=1 (J2); k2=1 (J3); T2=1 с (J4); k3=1 (J5) и Т3=5 с (J6).

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

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

Моделирование процессов поиска дефектов в виде уменьшения параметра T1 на 20% приводит к вычислению диагностических признаков по формуле (3): J1=5.793e-005; J2=0.1457; J3=0.266; J4=0.3127; J5=0.1405; J6=0.2147. Различимость дефекта: ΔJ=J2-J1=0.1457.

Для сравнения приведем диагностические признаки наличия неисправного блока с использованием экспоненциального веса при одном параметре интегрирования α=0.5 (Положительное решение от 11.04.2011 г. о выдаче патента на изобретение №2429518 по заявке №2010128421/08(040385), МКИ6 G05В 23/02, 2011): J1=2.293e-008; J2=1.469e-007; J3=0.7843; J4=0.7828; J5=0.07422; J6=0.07399. Различимость дефекта ΔJ=J2-J1=1.2397e-007.

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

Моделирование процессов поиска дефектов в виде уменьшения параметра Т2 на 20% для данного объекта диагностирования с использованием дифференциального веса и при таком же входном сигнале дает следующие значения диагностических признаков: J1=0.2086; J2=0.7451; J3=0.006981; J4=0.001517; J5=0.7004; J6=0.6441.

Различимость дефекта ΔJ=J3-J4=0.005464.

Для сравнения приведем диагностические признаки наличия неисправного блока с использованием экспоненциального веса при одном параметре интегрирования α=0.5: J1=0.7828; J2=0.7825; J3=3.695e-006; J4=1.52e-008; J5=0.7455; J6=0.7462. Различимость дефекта: ΔJ=J3-J4=3.6798e-006.

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

J1=0.2186; J2=0.1878; J3=0.5691; J4=0.6809; J5=0.003594; J6=2.192e-005.

Различимость дефекта: ΔJ=J5-J6=0.003594.

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

J1=0.07407; J2=0.07384; J3=0.7476; J4=0.7457; J5=7.831e-007; J6=1.531e-010.

Различимость дефекта ΔJ=J5-J6=7.829469е-007.

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

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

Похожие патенты RU2464616C1

название год авторы номер документа
СПОСОБ ПОИСКА НЕИСПРАВНОГО БЛОКА В ДИНАМИЧЕСКОЙ СИСТЕМЕ 2011
  • Воронин Владимир Викторович
  • Киселев Владислав Валерьевич
  • Шалобанов Сергей Викторович
  • Шалобанов Сергей Сергеевич
RU2451319C1
СПОСОБ ПОИСКА НЕИСПРАВНОГО БЛОКА В ДИСКРЕТНОЙ ДИНАМИЧЕСКОЙ СИСТЕМЕ 2012
  • Киселев Владислав Валерьевич
  • Шалобанов Сергей Викторович
  • Шалобанов Сергей Сергеевич
RU2506623C1
СПОСОБ ПОИСКА НЕИСПРАВНЫХ БЛОКОВ В ДИНАМИЧЕСКОЙ СИСТЕМЕ 2012
  • Киселев Владислав Валерьевич
  • Шалобанов Сергей Викторович
  • Шалобанов Сергей Сергеевич
RU2473106C1
СПОСОБ ПОИСКА НЕИСПРАВНОСТЕЙ ДИНАМИЧЕСКОГО БЛОКА В НЕПРЕРЫВНОЙ СИСТЕМЕ 2010
  • Шалобанов Сергей Сергеевич
RU2429518C1
СПОСОБ ПОИСКА НЕИСПРАВНОГО БЛОКА В ДИНАМИЧЕСКОЙ СИСТЕМЕ 2010
  • Шалобанов Сергей Викторович
  • Шалобанов Сергей Сергеевич
RU2439648C1
СПОСОБ ПОИСКА НЕИСПРАВНОГО БЛОКА В НЕПРЕРЫВНОЙ ДИНАМИЧЕСКОЙ СИСТЕМЕ 2013
  • Шалобанов Сергей Викторович
  • Шалобанов Сергей Сергеевич
RU2513504C1
Способ поиска топологического дефекта в непрерывной динамической системе на основе введения пробных отклонений 2016
  • Шалобанов Сергей Сергеевич
RU2616512C1
Способ поиска неисправного блока в непрерывной динамической системе на основе введения пробных отклонений и анализа знаков передач 2019
  • Шалобанов Сергей Викторович
  • Шалобанов Сергей Сергеевич
RU2711000C1
СПОСОБ ПОИСКА НЕИСПРАВНЫХ БЛОКОВ В НЕПРЕРЫВНОЙ ДИНАМИЧЕСКОЙ СИСТЕМЕ 2012
  • Шалобанов Сергей Викторович
  • Шалобанов Сергей Сергеевич
RU2519435C1
СПОСОБ ПОИСКА НЕИСПРАВНОГО БЛОКА В ДИСКРЕТНОЙ ДИНАМИЧЕСКОЙ СИСТЕМЕ 2012
  • Воронин Владимир Викторович
  • Шалобанов Сергей Сергеевич
RU2486569C1

Реферат патента 2012 года СПОСОБ ПОИСКА НЕИСПРАВНОСТЕЙ ДИНАМИЧЕСКОГО БЛОКА В НЕПРЕРЫВНОЙ СИСТЕМЕ

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

Формула изобретения RU 2 464 616 C1

1. Способ поиска неисправного блока в динамической системе, основанный на том, что фиксируют число m параметров, входящих в состав системы, определяют время контроля Tk≥ТПП, используют входной сигнал x(t) на интервале t∈[0, TK], фиксируют число k контрольных точек системы, регистрируют реакцию контролируемой системы fj(t), j=1, 2, …, k, регистрируют реакцию системы с номинальными характеристиками fj ном(t), j=1, 2, …, k, на интервале t∈[0, Tk] в k контрольных точках, определяют интегральные оценки выходных сигналов системы, для чего в момент подачи сигнала на вход системы с номинальными характеристиками одновременно начинают интегрирование сигналов системы в каждой из k контрольных точек путем подачи на первые входы k блоков перемножения выходных сигналов системы, на вторые входы блоков перемножения подают весовую функцию, выходные сигналы k блоков перемножения подают на входы k блоков интегрирования, интегрирование завершают в момент времени Тк, полученные в результате интегрирования оценки выходных сигналов регистрируют, определяют интегральные оценки выходных сигналов модели для каждой из k контрольных точек, полученные в результате пробных отклонений каждого из m параметров всех блоков, для чего в каждый параметр динамической системы вводят его пробное отклонение и находят интегральные оценки выходных сигналов модели для весовой функции и входного сигнала x(t), полученные в результате интегрирования оценки выходных сигналов для каждой из k каждого из m параметров всех блоков, для чего в каждую i-ю модель поочередно вводят соответствующее пробное отклонение i-го параметра передаточной функции всех блоков динамической системы и находят интегральные оценки выходных сигналов систем с пробными отклонениями при том же тестовом или рабочем сигнале x(t), полученные в результате интегрирования оценки выходных сигналов для каждой из k контрольных точек каждого из m пробных отклонений Pji(d), j=1, …, k; i=1, …, m, регистрируют, одновременно на вход контролируемой системы подают тестовый или рабочий сигнал x(t), определяют интегральные оценки сигналов контролируемой системы для k контрольных точек Fj(d), j=1, …, k, полученные значения регистрируют, определяют отклонения интегральных оценок сигналов модели, полученные в результате пробных отклонений соответствующих параметров всех блоков ΔPji(d)=Pji(d)-Fj ном(d), j=1, …, k; i=1, …, m, определяют нормированные значения отклонений интегральных оценок сигналов модели, полученные в результате пробных отклонений соответствующих параметров всех блоков

определяют отклонения интегральных оценок сигналов контролируемой системы для k контрольных точек ΔFj(d)=Fj(d)-Fj ном(d), j=1, …, k, определяют нормированные значения отклонений интегральных оценок сигналов контролируемой системы контрольных точек и каждое из m пробных отклонений регистрируют, определяют отклонения интегральных оценок сигналов модели, полученные в результате пробных отклонений соответствующих параметров ΔPji(d)=Pji(d)-Fj ном(d), j=1, …, k; i=1, …, m, определяют нормированные значения отклонений интегральных оценок сигналов модели, полученные в результате пробных отклонений соответствующих параметров из соотношения

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

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

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

определяют интегральные оценки выходных сигналов Fj ном(d), j=1, …, k, системы с номинальными характеристиками, для чего в момент подачи входного сигнала на вход системы с номинальными характеристиками одновременно начинают интегрирование сигналов системы в каждой из k контрольных точек для весовой функции, путем подачи на первые входы k блоков перемножения сигналов системы на вторые входы блоков перемножения подают среднее арифметическое модулей производных по времени от выходных сигналов системы с номинальными характеристиками, выходные сигналы k блоков перемножения подают на входы k блоков интегрирования, интегрирование завершают в момент времени Тк, полученные в результате интегрирования оценки выходных сигналов Fj ном(d), j=1, …, k регистрируют, аналогично определяют интегральные оценки сигналов m моделей для каждой из k контрольных точек, полученные в результате пробных отклонений каждого из m параметров, полученные в результате интегрирования оценки выходных сигналов для каждой из k контрольных точек, каждого из m пробных отклонений Pji(d), j=1, …, k; i=1, …, m, используют для вычисления диагностических признаков.

Документы, цитированные в отчете о поиске Патент 2012 года RU2464616C1

RU 2009123999 А, 27.12.2010
СПОСОБ ПОИСКА НЕИСПРАВНОСТЕЙ ДИНАМИЧЕСКОГО БЛОКА В НЕПРЕРЫВНОЙ СИСТЕМЕ 2010
  • Шалобанов Сергей Сергеевич
RU2429518C1
JP 2009290349 A, 10.12.2009
US 4851985, 25.07.1989
EP 1324165 A2, 02.07.2003.

RU 2 464 616 C1

Авторы

Киселев Владислав Валерьевич

Шалобанов Сергей Сергеевич

Даты

2012-10-20Публикация

2011-11-01Подача