UP_nadezhnost_i_diagnostika

Нечеткие методы ФД

Рис.6.12. Классификация методов, использующих нечеткий анализ.

где x – n-мерный вектор состояния, u – m-мерный вектор входа, y – p-мерный вектор выхода, F – n x n-матрица динамики, G – входная n x m-матрица, H

– выходная p x n-матрица, – функция динамики.

При диагностировании линейной системы уравнения наблюдателя, настроенного на i-е техническое состояние Si , имеют вид:

K D s

Далее для простоты полагается, что в рассматриваемой системе отсутствуют возмущения. С известными решениями при компенсации возмущений можно ознакомиться, например, по известным работам [11].

Как следует из вышеизложенного в рассматриваемой задаче синтеза средств ФД можно выделить два основных вопроса – синтез банка

101

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

Нечеткое правило принятия решения об отказе. Важнейшей особенностью рассматриваемого метода является использование нечеткого анализа (см. Приложение 3) при принятии решения об отказе [26, 30 – 32], который привлекает аппарат нечетких множеств при описании невязок и отказов. Далее при описании нечетких множеств ограничимся применением трапециидальных функций принадлежности, которые для параметров, принимающих, например, только неотрицательные значения, описываются выражениями:

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

характеристики, отражающей степень уверенности разработчика в том, что в данный момент система находится в состоянии Si . Вычисление Ki

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

Казалось бы, что решение об отказе можно принимать по минимуму невязки, однако использование коэффициентов уверенности дает определенное преимущество. В целом ряде случаев оно открывает путь для применения более эффективных вариантов организации банков наблюдателей, а именно, банков взаимодействующих наблюдателей. Кроме того, использование коэффициентов уверенности позволяет сделать анализ технического состояния ДС более детальным, отделяя более очевидные ситуации, когда присутствие отказа в ДС практически не вызывает сомнения, от менее очевидных. Классификация ситуаций осуществляется путем соотнесения максимального из коэффициентов уверенности Ki с некоторым порогом A, т.е. должно выполняться

Предполагается, что невязка νi i 0,M, формируемая при сопоставлении

выходов системы и i-го наблюдателя, может быть представлена лингвистической переменной (рис. 6.14), например, с двумя термами –

102

«малая» и «большая», для которых заданы функции принадлежности μ 0νi и

объясняется переходными процессами, сопровождающими оценивание, отсутствием на практике полной адекватности используемой при синтезе наблюдателя модели диагностируемой системы, неучтенными возмущениями ее динамики или выхода. Терм «большая» соответствует ситуации, когда модель, использованная при синтезе наблюдателя, существенно неадекватна текущему техническому состоянию диагностируемой системы. Так бывает, если, например, диагностируемая система находится в i-м техническом состояние, а наблюдатель настроен на j–е техническое состояние. При этом

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

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

103

степеней для всех состояний:

~

K i Nμi .

~

μ j

j 0

Подытожим настоящий раздел, сформулировав алгоритм ФД произвольных отказов в рассматриваемой постановке.

Алгоритм 6.1.

1.Формирование перечня отказов.

2.Синтез наблюдателя для каждого из отказов.

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

4.Принятие решения об отказах путем формирования коэффициентов уверенности.

Для иллюстрации предложенного алгоритма снова обратимся к примеру на рис. 6.13.

Пример 6.1. Синтезируем средства диагностирования для обрыва обратной связи по скорости. Номинальное поведение системы описывается уравнением:

104

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

6.3.Функциональное диагностирование при поиске отказов в пространстве сигналов

При диагностировании в пространстве сигналов отказ моделируется как дополнительное слагаемое δ в уравнении динамики, т.е. для исходной линейной системы

При этом число типов M однократных отказов равно размерности n вектора состояния диагностируемой системы. Первый тип моделируется дополнительным слагаемым в первом уравнении динамики для первой компоненты вектора x, второй тип – во втором уравнении для первой компоненты вектора x и т.д. Отказы внутри типа различаются уровнем слагаемого δ.

Как следует из рис.6.11 в задаче синтеза средств диагностирования, кроме формирования правила принятия решений, необходимо определить структуру банка и самих наблюдателей. При диагностировании линейной системы уравнения наблюдателя имеют вид [1]:

для нелинейной:

получаемого добавлением к вектору xk состояния диагностируемой системы переменной i , моделирующей отказ. Причем с учетом предположения о постоянстве значения i уравнение для нее имеет вид:

Врезультате матрицы диагностируемой линейной системы при наличии

вней i-го отказа принимают вид:

105

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

этого составного вектора, в том числе формируется и оценка ˆi значения переменной i . Далее, пользуясь материалом предыдущего параграфа, нетрудно синтезировать средства ФД с банком из независимых наблюдателей. Однако оказывается, что при диагностировании в пространстве сигналов можно применить несколько более сложную, но зато и более эффективную организацию банка, когда оценка в каждом наблюдателе формируется с учетом оценок в других наблюдателях. Будем называть этот случай использованием банка взаимодействующих наблюдателей.

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

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

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

106

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

Проведем качественный сравнительный анализ методов ФД с независимыми и взаимодействующими наблюдателями, опираясь на выражение для невязки в случае независимых наблюдателей:

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

Рис.6.16. Структура алгоритма диагностирования с взаимодействующими наблюдателями, анализирующими состояния.

только от матрицы динамики F i наблюдателя, несущей информацию об отказе, но и от текущей оценки, формируемой в наблюдателе. Ясно, что в силу различия матриц динамики траектория движения каждого наблюдателя в пространстве состояний, а значит, и результирующая оценка в общем случае будет своя, а невязка большая. Последнее касается и наблюдателя, адекватного техническому состоянию ДС, особенно во время переходного процесса. Таким образом, в общем случае отсутствуют гарантии минимальности невязки адекватного наблюдателя в каждый момент времени.

Ситуация меняется при использовании банка взаимодействующих наблюдателей, поскольку при формировании невязок во всех наблюдателях используется одна и та же оценка состояния ДС, полученная осреднением:

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

107

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

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

Отказ

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

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

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

108

Si S j , то он использует модель системы в j-м техническом состоянии, а при вычислениях в нем используется оценка xˆ i вектора состояния, сформированная в предположении об i-м техническом состоянии:

M i

xˆ ik Ki, j xˆki, j , j 0

где M i - число технических состояний, из которых возможны переходы в Si , Ki, j - коэффициент уверенности для перехода Si S j , который вычисляется

аналогично коэффициентам уверенности для технических состояний. Коэффициент уверенности для каждого технического состояния вычисляется путем суммирования коэффициентов уверенности, соответствующих переходам в это состояние:

Mi

Ki Ki,j . j 0

На рис. 6.17 приведена структура алгоритма для простейшего случая одного отказа, вследствие чего техническое состояние может принимать только два значения, однако число рассматриваемых переходов равно трем (переход из неработоспособного состояния в работоспособное не рассматривается). В результате в средствах диагностирования применяются три наблюдателя. При другой интерпретации структурных блоков можно заметить сходство данного алгоритма с известным алгоритмом диагностирования стохастической системы, использующим полигауссову аппроксимацию апостериорной плотности вероятности [26].

Пример 6.2. Рассмотрим пример системы, заимствованный из работы [40] и характеризующийся матрицами:

109

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

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

Для этого примера в среде Simulink было проведено моделирование задачи диагностирования в пространстве сигналов. Моделировался перемежающийся (с возвратом в номинальное состояние) отказ θ в виде меандра на первом интеграторе (первая диаграмма на рис. 6.18а и 6.18б).

На вход системы подавался синусоидальный сигнал с амплитудой 0,5. На рис. 6.18 приведены временные диаграммы полученных коэффициентов уверенности для случаев использования независимых и взаимодействующих наблюдателей. Видно, что в случае независимых наблюдателей (рис. 6.18а) средства диагностирования не формируют ожидаемого результата, когда до появления отказа коэффициент К0 должен принимать устойчивое значение, близкое к единице, а коэффициент К1 – значение близкое к нулю. После появления отказа коэффициент К1 должен принимать устойчивое значение, близкое к единице, а К0 – значение близкое к нулю. При использовании взаимодействующих наблюдателей (рис. 6.18б), сопоставленных с техническим состояниями, временные диаграммы демонстрируют адекватную работу средств диагностирования. Так на интервале времени до появления отказа коэффициент уверенности К0 для работоспособного технического состояния после переходного процесса, связанного с начальным оцениванием, принимает значение, равное единице. После появления отказа соответствующий ему коэффициент уверенности устойчиво принимает значение, близкое единице.

110