RU2680928C1 - Способ поиска неисправного блока в непрерывной динамической системе на основе функции чувствительности - Google Patents

Способ поиска неисправного блока в непрерывной динамической системе на основе функции чувствительности Download PDF

Info

Publication number
RU2680928C1
RU2680928C1 RU2018111888A RU2018111888A RU2680928C1 RU 2680928 C1 RU2680928 C1 RU 2680928C1 RU 2018111888 A RU2018111888 A RU 2018111888A RU 2018111888 A RU2018111888 A RU 2018111888A RU 2680928 C1 RU2680928 C1 RU 2680928C1
Authority
RU
Russia
Prior art keywords
determined
input
output signals
estimates
signals
Prior art date
Application number
RU2018111888A
Other languages
English (en)
Inventor
Сергей Викторович Шалобанов
Сергей Сергеевич Шалобанов
Original Assignee
Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет"
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" filed Critical Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет"
Priority to RU2018111888A priority Critical patent/RU2680928C1/ru
Application granted granted Critical
Publication of RU2680928C1 publication Critical patent/RU2680928C1/ru

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B23/00Testing or monitoring of control systems or parts thereof
    • G05B23/02Electric testing or monitoring

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

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

Description

Изобретение относится к области контроля и диагностирования систем автоматического управления и их элементов.
Известен способ поиска неисправного блока в динамической системе (Способ поиска неисправного блока в динамической системе: пат. 2435189 Рос. Федерация : МПК7 G05B 23/02 (2006.01) / Шалобанов С.В., Шалобанов С.С. - №2009123999/08; заявл. 23.06.2009; опубл. 27.11.2011, Бюл. №33).
Недостатком этого способа является то, что он использует задание величин относительных отклонений параметров передаточных функций для моделей с пробными отклонениями.
Наиболее близким техническим решением (прототипом) является способ поиска неисправностей динамического блока в непрерывной системе на основе функции чувствительности (Способ поиска неисправностей динамического блока в непрерывной системе на основе функции чувствительности: пат. 2580405 Рос. Федерация : МПК G05B 23/02 (2006.01) / Шалобанов С.В., Шалобанов С.С. - №2015110545/08; заявл. 24.03.2015; опубл. 10.04.2016, Бюл. №10).
Недостатком этого способа является то, что он использует вычисление производных функции чувствительности для параметров передаточных функций моделей.
Технической задачей, на решение которой направлено данное изобретение, является уменьшение вычислительных затрат, связанных с реализацией моделей со структурной чувствительностью блоков передаточных функций.
Поставленная задача достигается тем, что предварительно регистрируют реакцию заведомо исправной системы fj ном(t), j=1,…,k на интервале t ∈ [0, TK] в к контрольных точках, и определяют интегральные оценки выходных сигналов Fj ном (α), j=1,…,k системы, для чего в момент подачи входного сигнала на вход системы с номинальными характеристиками одновременно начинают интегрирование сигналов системы управления в каждой из k контрольных точек с весами e-αt где
Figure 00000001
, путем подачи на первые входы к блоков перемножения сигналов системы управления, на вторые входы блоков перемножения подают экспоненциальный сигнал e-αt, выходные сигналы k блоков перемножения подают на входы k блоков интегрирования, интегрирование завершают в момент времени Тк, полученные в результате интегрирования оценки выходных сигналов Fj ном (α), j=1,…,k регистрируют, фиксируют число m блоков системы, определяют интегральные оценки выходных сигналов модели для каждой из k контрольных точек и каждой из m моделей с функциями структурной чувствительности для каждого из m блоков динамической системы, для чего соединяют связью две модели: на вход первой модели подают тестовый сигнал x(t), выходом первой модели определяют вход контролируемого блока, соединяют выход первой модели со входом второй, входом второй модели становится выход контролируемого блока, находят интегральные оценки выходных сигналов системы для параметра α и входного сигнала x(t), полученные в результате интегрирования оценки выходных сигналов для каждой из k контрольных точек и каждой из m совмещенных моделей со структурной функцией чувствительности
Vji(α) j=1,…k; i=1,…m регистрируют, определяют нормированные значения интегральных оценок выходных сигналов модели, полученные в результате структурной функции чувствительности каждого из соответствующих блоков из соотношения
Figure 00000002
замещают систему с номинальными характеристиками контролируемой, на вход системы подают аналогичный входной сигнал x(t), определяют интегральные оценки сигналов контролируемой системы для k контрольных точек Fj(α), j=1,…,k для параметра интегрирования α, определяют отклонения интегральных оценок сигналов контролируемой системы для к контрольных точек от номинальных значений ΔFj(α)=Fj(α)-Fj ном (α), j=1,…,k, определяют нормированные значения отклонений интегральных оценок сигналов контролируемой системы из соотношения
Figure 00000003
определяют диагностические признаки из соотношения
Figure 00000004
по минимуму значения диагностического признака определяют порядковый номер структурного дефекта.
Таким образом, предлагаемый способ поиска структурного дефекта сводится к выполнению следующих операций:
1. В качестве динамической системы рассматривают систему, состоящую из произвольно соединенных динамических блоков, с количеством рассматриваемых блоков m.
2. Предварительно определяют время контроля ТК≥ТПП, где ТПП - время переходного процесса системы. Время переходного процесса оценивают для номинальных значений параметров динамической системы.
3. Определяют параметр интегрального преобразования сигналов из соотношения
Figure 00000005
4. Фиксируют число контрольных точек k.
5. Предварительно определяют нормированные векторы
Figure 00000006
интегральных оценок выходных сигналов модели, полученные в результате интегрирования функций чувствительности i-го блока каждого из m блоков для номинальных значений параметров передаточных функций блоков и определенного выше параметра α, для чего выполняют пункты 6-10.
6. Подают входной сигнал x(t) (единичный ступенчатый, линейно возрастающий, прямоугольный импульсный и т.д.) на вход системы управления с номинальными характеристиками. Принципиальных ограничений на вид входного тестового воздействия предлагаемый способ не предусматривает.
7. Регистрируют реакцию системы fj ном(t), j=1,…,k на интервале t ∈ [0, TK] в k контрольных точках и определяют интегральные оценки выходных сигналов fj ном(α), j=1,…,k системы. Для этого в момент подачи входного сигнала на вход системы управления с номинальными характеристиками одновременно начинают интегрирование сигналов системы управления в каждой из k контрольных точек с весами e-αt, где
Figure 00000007
для чего сигналы системы управления подают на первые входы k блоков перемножения, на вторые входы блоков перемножения подают экспоненциальный сигнал e-αt, выходные сигналы k блоков перемножения подают на входы к блоков интегрирования, интегрирование завершают в момент времени Тк, полученные в результате интегрирования оценки выходных сигналов fj ном(α), j=1,…,k регистрируют.
8. Определяют интегральные оценки выходных сигналов модели чувствительности для каждой из к контрольных точек, полученные в результате использования структурной функции чувствительности каждого из m блоков, для чего поочередно для каждого блока динамической системы соединяют связью две модели: на вход первой модели подают тестовый сигнал x(t), выходом первой модели становится вход контролируемого блока, соединяют выход первой модели со входом второй, входом второй модели становится выход контролируемого блока, и выполняют пункты 6 и 7 для одного и того же входного сигнала x(t). Полученные в результате интегрирования оценки выходных сигналов для каждой из k контрольных точек и каждой из m совмещенных моделей со структурной функцией чувствительности Vj i(α), j=1,…,k; i=1,…,m регистрируют.
9. Определяют нормированные значения интегральных оценок выходных сигналов модели, полученные в результате использования структурной функции чувствительности по формуле
Figure 00000008
10. Замещают систему с номинальными характеристиками контролируемой. На вход системы подают аналогичный входной сигнал x(t).
11. Определяют интегральные оценки выходных сигналов контролируемой системы для к контрольных точек Fj(α), j=1,…k, осуществляя операции, описанные в пунктах 6 и 7 применительно к контролируемой системе.
12. Определяют отклонения интегральных оценок выходных сигналов контролируемой системы для k контрольных точек от номинальных значений ΔFj(α)=Fj(α)-Fj ном (α), j=1,…,k.
13. Вычисляют нормированные значения отклонений интегральных оценок выходных сигналов контролируемой системы по формуле
Figure 00000009
14. Вычисляют диагностические признаки наличия неисправного блока по формуле (3).
15. По минимуму значения диагностического признака определяют структурный дефект.
Рассмотрим реализацию предлагаемого способа поиска дефектного блока для системы, структурная схема которой представлена на рисунке (см. фиг. Структурная схема объекта диагностирования).
Передаточные функции блока 1:
Figure 00000010
блока 2:
Figure 00000011
блока 3:
Figure 00000012
номинальные значения параметров: T1=5 с; k1=1; k2=1; Т2=1 с; k3=1; Т3=5 с. При поиске одиночного дефекта в виде отклонения постоянной времени Т1=4 с в первом блоке путем подачи ступенчатого тестового входного сигнала единичной амплитуды и интегральных преобразований сигналов для параметра α=0.5 и Тк=10 с получены значения диагностических признаков на основе пробных отклонений параметров модели (Способ поиска неисправного блока в динамической системе: пат. 2435189 Рос. Федерация : МПК7 G05B 23/02 (2006.01) / Шалобанов С.В., Шалобанов С.С. - №2009123999/08; заявл. 23.06.2009; опубл. 27.11.2011, Бюл. №33) при использовании трех контрольных точек (КТ1, КТ2 и КТ3), расположенных на выходах блоков: J1=0; J2=0.78; J3=0.074. Минимальное значение признака J1 однозначно указывает на наличие дефекта в первом блоке, а разность между третьим и первым признаками может количественно характеризовать практическую (апостериорную) различимость этого дефекта. Тот же дефект, найденный путем применения функции чувствительности и вычислений по формуле (3) с использованием тех же контрольных точек (КТ1, КТ2 и КТ3), времени контроля Тк=10 с и параметра интегрирования α=0.5 дает следующие значения диагностических признаков: J1=0; J2=0.78; J3=0.074. Анализ значений диагностических признаков показывает, что значения второго и третьего признака, полученные при использовании функции чувствительности, практически такие же, как и при использовании пробных отклонений. Это позволяет сделать вывод, что практическая различимость дефекта первого блока практически такая же, как и при использовании заявляемого способа. Различимости дефектов второго и третьего блоков при поиске их с использованием функции чувствительности тоже практически такие же, как и при использовании пробных отклонений параметров модели.
Моделирование процессов поиска дефектов во втором и третьем блоках для данного объекта диагностирования, при том же параметре интегрирования α и при единичном ступенчатом входном сигнале дает следующие значения диагностических признаков:
При наличии дефекта в блоке №2 (в виде уменьшения параметра Т2 на 20%, дефект №2): J1=0.7829; J2=0; J3=0.746.
При наличии дефекта в блоке №3 (в виде уменьшения параметра Т3 на 20%, дефект №3) J1=0.07406; J2=0.7471; J3=0.
Минимальное значение диагностического признака во всех случаях правильно указывает на дефектный блок.

Claims (1)

  1. Способ поиска неисправного блока в непрерывной динамической системе на основе функции чувствительности, основанный на том, что фиксируют число блоков m, входящих в состав системы, определяют время контроля ТК≥ТПП, где ТПП - время переходного процесса системы, фиксируют число k контрольных точек системы, используют входной сигнал на интервале t∈[0, TК], регистрируют реакцию объекта диагностирования и реакцию заведомо исправной системы ƒj ном(t), j=1, …, k на интервале t∈[0, TK] в k контрольных точках, замещают систему с номинальными характеристиками контролируемой, на вход системы подают аналогичный входной сигнал x(t), вычисляют диагностические признаки наличия дефекта, по минимуму диагностического признака определяют дефект, отличающийся тем, что определяют параметр интегрального преобразования сигналов из соотношения
    Figure 00000013
    в качестве динамических характеристик системы используют интегральные оценки сигналов, определяют интегральные оценки выходных сигналов
    Figure 00000014
    j=1, …, k исправной системы, для чего в момент подачи входного сигнала на вход системы с номинальными характеристиками одновременно начинают интегрирование сигналов системы управления в каждой из k контрольных точек с весами
    Figure 00000015
    где
    Figure 00000016
    путем подачи на первые входы k блоков перемножения сигналов системы управления, на вторые входы блоков перемножения подают экспоненциальный сигнал
    Figure 00000017
    выходные сигналы k блоков перемножения подают на входы k блоков интегрирования, интегрирование завершают в момент времени Тк, полученные в результате интегрирования оценки выходных сигналов
    Figure 00000018
    j=1, …, k регистрируют, определяют интегральные оценки сигналов контролируемой системы для k контрольных точек
    Figure 00000019
    j=1, …, k для параметра
    Figure 00000020
    определяют отклонения интегральных оценок сигналов контролируемой системы для k контрольных точек от номинальных значений
    Figure 00000021
    j=1, …, k, определяют нормированные значения отклонений интегральных оценок сигналов контролируемой системы из соотношения
    Figure 00000022
    определяют интегральные оценки выходных сигналов модели для каждой из k контрольных точек, полученные в результате использования функции структурной чувствительности каждого из m блоков динамической системы, для чего соединяют связью две модели: на вход первой модели подают тестовый сигнал x(t), выходом первой модели становится вход контролируемого блока, соединяют выход первой модели со входом второй, входом второй модели становится выход контролируемого блока, и находят интегральные оценки выходных сигналов системы для параметра
    Figure 00000020
    и входного сигнала x(t), полученные в результате интегрирования оценки выходных сигналов для каждой из k контрольных точек и каждой из m совмещенных моделей со структурной функцией чувствительности
    Figure 00000023
    j=1, …, k; i=1, …, m регистрируют, определяют нормированные значения интегральных оценок выходных сигналов модели, полученные в результате применения структурной функции чувствительности для каждого из соответствующих блоков из соотношения
    Figure 00000024
    вычисляют диагностические признаки из соотношения
    Figure 00000025
    i=1, …, m, по минимуму диагностического признака определяют дефектный блок.
RU2018111888A 2018-04-02 2018-04-02 Способ поиска неисправного блока в непрерывной динамической системе на основе функции чувствительности RU2680928C1 (ru)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2018111888A RU2680928C1 (ru) 2018-04-02 2018-04-02 Способ поиска неисправного блока в непрерывной динамической системе на основе функции чувствительности

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2018111888A RU2680928C1 (ru) 2018-04-02 2018-04-02 Способ поиска неисправного блока в непрерывной динамической системе на основе функции чувствительности

Publications (1)

Publication Number Publication Date
RU2680928C1 true RU2680928C1 (ru) 2019-02-28

Family

ID=65632619

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2018111888A RU2680928C1 (ru) 2018-04-02 2018-04-02 Способ поиска неисправного блока в непрерывной динамической системе на основе функции чувствительности

Country Status (1)

Country Link
RU (1) RU2680928C1 (ru)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2719747C1 (ru) * 2019-05-22 2020-04-22 Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" Способ поиска неисправного блока в непрерывной динамической системе на основе функции чувствительности
RU2744823C1 (ru) * 2020-06-22 2021-03-16 Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" Способ поиска неисправных блоков в непрерывной динамической системе на основе введения пробных отклонений

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6326758B1 (en) * 1999-12-15 2001-12-04 Reliance Electric Technologies, Llc Integrated diagnostics and control systems
WO2005071582A2 (en) * 2004-01-27 2005-08-04 British Telecommunications Public Limited Company Detection of abnormal behaviour in dynamic systems
RU2580405C1 (ru) * 2015-03-24 2016-04-10 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" Способ поиска неисправностей динамического блока в непрерывной системе на основе функции чувствительности
RU2586859C1 (ru) * 2015-03-11 2016-06-10 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" Способ поиска неисправного блока в непрерывной динамической системе на основе смены позиции входного сигнала
RU2613630C1 (ru) * 2016-03-09 2017-03-21 Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" Способ поиска неисправного блока в непрерывной динамической системе на основе введения пробных отклонений

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6326758B1 (en) * 1999-12-15 2001-12-04 Reliance Electric Technologies, Llc Integrated diagnostics and control systems
WO2005071582A2 (en) * 2004-01-27 2005-08-04 British Telecommunications Public Limited Company Detection of abnormal behaviour in dynamic systems
RU2586859C1 (ru) * 2015-03-11 2016-06-10 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" Способ поиска неисправного блока в непрерывной динамической системе на основе смены позиции входного сигнала
RU2580405C1 (ru) * 2015-03-24 2016-04-10 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тихоокеанский государственный университет" Способ поиска неисправностей динамического блока в непрерывной системе на основе функции чувствительности
RU2613630C1 (ru) * 2016-03-09 2017-03-21 Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" Способ поиска неисправного блока в непрерывной динамической системе на основе введения пробных отклонений

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2719747C1 (ru) * 2019-05-22 2020-04-22 Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" Способ поиска неисправного блока в непрерывной динамической системе на основе функции чувствительности
RU2744823C1 (ru) * 2020-06-22 2021-03-16 Федеральное государственное бюджетное образовательное учреждение высшего образования "Тихоокеанский государственный университет" Способ поиска неисправных блоков в непрерывной динамической системе на основе введения пробных отклонений

Similar Documents

Publication Publication Date Title
RU2429518C1 (ru) Способ поиска неисправностей динамического блока в непрерывной системе
RU2435189C2 (ru) Способ поиска неисправного блока в динамической системе
RU2541857C1 (ru) Способ поиска неисправностей в непрерывной динамической системе на основе введения пробных отклонений
RU2439647C1 (ru) Способ поиска неисправного блока в непрерывной динамической системе
RU2450309C1 (ru) Способ поиска неисправностей динамического блока в непрерывной системе
RU2528135C1 (ru) Способ поиска неисправного блока в непрерывной динамической системе на основе смены позиции входного сигнала
RU2439648C1 (ru) Способ поиска неисправного блока в динамической системе
RU2680928C1 (ru) Способ поиска неисправного блока в непрерывной динамической системе на основе функции чувствительности
RU2444774C1 (ru) Способ поиска неисправного блока в дискретной динамической системе
RU2586859C1 (ru) Способ поиска неисправного блока в непрерывной динамической системе на основе смены позиции входного сигнала
RU2613630C1 (ru) Способ поиска неисправного блока в непрерывной динамической системе на основе введения пробных отклонений
RU2473105C1 (ru) Способ поиска неисправностей блоков в непрерывной динамической системе
RU2461861C1 (ru) Способ поиска неисправного блока в непрерывной динамической системе
RU2676365C1 (ru) Способ поиска неисправного блока в непрерывной динамической системе на основе введения пробных отклонений
RU2464616C1 (ru) Способ поиска неисправностей динамического блока в непрерывной системе
RU2506623C1 (ru) Способ поиска неисправного блока в дискретной динамической системе
RU2541896C1 (ru) Способ поиска неисправного блока в дискретной динамической системе на основе анализа знаков передач сигналов
RU2579543C1 (ru) Способ поиска неисправного блока в дискретной динамической системе на основе смены позиции входного сигнала
RU2453898C1 (ru) Способ поиска неисправных блоков в динамической системе
RU2451319C1 (ru) Способ поиска неисправного блока в динамической системе
RU2669035C1 (ru) Способ поиска неисправного блока в непрерывной динамической системе на основе смены позиции входного сигнала
RU2658547C1 (ru) Способ поиска неисправного блока в дискретной динамической системе на основе смены позиции входного сигнала
RU2661180C1 (ru) Способ поиска неисправного блока в непрерывной динамической системе на основе смены позиции входного сигнала
RU2580405C1 (ru) Способ поиска неисправностей динамического блока в непрерывной системе на основе функции чувствительности
RU2562429C1 (ru) Способ поиска неисправного блока в непрерывной динамической системе на основе смены позиции входного сигнала

Legal Events

Date Code Title Description
MM4A The patent is invalid due to non-payment of fees

Effective date: 20200403