Применение метода исключения варьируемого параметра при решении задач диагностирования
Портнягин Н. Н., Пюкке Г. А.,
Применение детерминированной модели диагностирования (2. 3) при решении задачи построения области работоспособных состояний ОД позволяет контролировать появление только одиночных дефектов. При этом необходимым условием проведения процедуры диагностирования является наличие предварительной информации о возникновении в системе только одиночных отказов. При появлении двойных, тройных и т. д. отказов область работоспособности на основе детерминированной модели (2.3) не определена, в силу потери однозначности соответствия между множеством значений диагностируемых параметров gi и совокупностью диагностических признаков {K1, K2 }. Недостатки такого рода при решении многомерных задач характерны для детерминированных моделей. Одним из способов их преодоления является введение в рассмотрение вероятностных характеристик при оценке состояния системы. Действительно, при рассмотрении различных вариантов сочетания вариаций параметров компонент в соотношениях (2. 3) количество попаданий в произвольную точку пространства диагностических признаков {K1, K2 } зависит от положения этой точки. Например, для трехизоварной картины (изовары 1, 2, 3; рис. 3. 6, гл. 3 ) при заданных интервалах изменения параметров составляющих компонент gi, i = 1, 2, 3, попадание в точку 1 можно выполнить посредством одновременной вариации второй и третьей компонент. Для попадания в точку 2 – первой и второй или второй и третьей. Для попадания в точку 3 – первой и второй или второй и третьей или первой и третьей. А для попадания в точку 4 необходима одновременная вариация уже трех компонент. Вариацией параметров трех компонент также можно переместиться в точки 1, 2, 3.
То есть каждой точке плоскости К1 , К2 может быть поставлено в соответствие счетное множество векторов параметров компонент, а в поле пространства диагностирования имеет место статистическая устойчивость распределения событий. Для количественной оценки полученной закономерности необходимо посредством модели (2. 3) вычислить количество векторов параметров компонент, принадлежащих каждой точке пространства диагностирования {K1, K2 } (рассматриваются все векторы с различной кратностью вариаций). Введем множество конечного числа событий А, перечисляющее все возможные состояния системы N = 2m . Это множество состоит из конечного числа подмножеств Аj , j = 0,1,2,….,m; каждое из которых включает в себя состояния нулевых, одиночных, двойных и т. д., m – кратных дефектов. Их количество определяется комбинаторно вышеописанным способом (Сm0 – количество нулевых; Cm1 – одиночных; Сm2 – двойных и т. д., Сmm – m – кратных дефектов). Под дефектом будем понимать выход параметра компоненты gi за установленные пределы [gi min , gi max], и нарушение одного или нескольких соотношений наблюдаемых признаков K1, K2. Для каждой точки пространства наблюдаемых диагностических признаков {K1, K2 } введем счетное множество событий Вк,определяющих количество реализаций вектора параметров компонент, приводящие в данную точку). Введением множества Вк формализуется оценка многозначности функции состояния ОД при его наблюдении в пространстве диагностирования. Среди совокупности реализаций, принадлежащих данной точке, отберем подмножества векторов соответствующих различной кратности вариаций. Тогда для каждой точки пространства диагностических признаков можно определить условную вероятность P(Aj | Bk) наступления события Аj при условии нахождения системы в выбранной точке. Т. к. множество Вк состоит из непересекающихся элементов, то справедлива формула полной вероятности: P(Aj) = P(Aj ¦Bk), где P(Aj) – вероятность наступления j – го события, удовлетворяющая следующим условиям 0 ≤ P(Aj) ≤ 1, P(A) =
P(Aj). Таким образом, точкам плоскостиК1 , К2 может быть поставлено в соответствие конечное множество векторов состояний ОД, в которых он может находиться, и вычислены условные вероятности P0(A0| Bk); P1(A1| Bk); P2(A2| Bk);………..; Pm(Am| Bk); возникновения нулевых, одиночных, двойных, и т. д., m – кратных отклонений параметров компонент в каждой точке пространства диагностирования {K1 , K2}. Это дает возможность каждой точке плоскости K1 , K2 поставить в соответствие вектор вероятностей всех 2m возможных состояний: Р = [Р0 Р1(1) P1(2) ……..P1(m) P2(1) P2(2)……..P2(Cm2)……..Pm(1)Pm(2)……..Pm(Cmm)] (индекс в круглых скобках обозначает номер состояния в подмножестве данной кратности; нижний индекс обозначает кратность вариации). Количество попаданий в любую фиксированную точку пространства {K1 , K2} при всех сочетаниях вариаций параметрами различных компонент ОД образуют полную группу событий. Для вектора вероятностей состояний можно написать: P0 +
PiN = 1, где R = Cmi. Решается задача распределения в пространстве диагностирования вероятностей отказов различной кратности j, (j = 0,1,2,….,m) и определении положений их максимумов. Как уже отмечалось, решить аналитически такую задачу не представляется возможным. При решении таких сложных задач, как правило, применяют метод статистических испытаний. Для его реализации можно использовать современные среды прикладного программирования, такие какMATLAB, VISUAL BASIC, EXEL и др. Метод позволяет оценить все функции условных вероятностей и решить основные задачи диагностирования. Действительно, по максимуму условной вероятности можно определить наиболее вероятные состояния ОД, если в каждой точке плоскости К1, К2 известны значения компонент вектора условных вероятностей. Как было отмечено выше, в пространстве диагностирования, при различных сочетаниях вариаций параметров компонент gi наблюдается статистическая устойчивость распределения вероятностей конечных положений точки состояния системы. Такая закономерность обусловлена наложением в каждой точке плоскости К1, К2 на точки данной кратности вариации точек вариации более высокого порядка кратности. Это дает возможность построения в пространстве диагностирования совокупности областей работоспособности ОД с заданным порогом вероятности.
.
Для построения функции вероятности работоспособности ОД всю совокупность точек плоскости К1, К2 разделяем на два непересекающихся подмножества работоспособных и неработоспособных состояний. Для выполнения процедуры деления используются условия работоспособности gimin ≥gi ≥ gi max , i = 1, m и исходные соотношения для диагностических признаков: К1 = К1 (g1, g2,…..,gm); K2 = K2 (g1, g2,…..,gm). Задается шаг изменения параметров компонент Δgi по двум направлениям Δgi+ и Δgi- каждой изовары от состояния равновесия gi ном , и вычисляются условные вероятности P(Aj) = P(Aj¦Bk), только тех точек, которые удовлетворяют системе неравенств: gi min ≥gi ≥ gi max , i = 1, m.. Тогда в каждой точке области работоспособности будут накапливаться условные вероятности, соответствующие различным сочетаниям вариаций gi и удовлетворяющие системе gi min ≥gi ≥ gi max. При машинной обработке информации выполняется вариация величинами всех m параметров составляющих компонент. Программа производит регистрацию общего количества вариаций различной кратности, приводящих в выбранную точку плоскости К1, К2 . Процедура выполняется для каждой точки пространства диагностических признаков. При таких условиях каждая кратность вариации будет иметь свой максимум вероятности в определенной подгруппе точек плоскости К1, К2. Расположение максимумов различной кратности вариаций на плоскости К1, К2 зависит от топологии изоварной картины, которая в свою очередь зависит от топологии ОД и варианта картины, выбранного по оптимизационным критериям (чувствительности, равномерности чувствительности, эквидистантности). Устойчивую закономерность распределения максимумов можно объяснить отсутствием в удаленных от центра изоварной картины точках вариаций низкой кратности. Вместе с тем вариации высокой кратности приводят к точкам на любых расстояниях от центра (в том числе и в центре). Поэтому на определенных расстояниях от центра будут образовываться зоны устойчивых максимумов вариаций различной кратности.
На рис. 4.3а приводится программа построения функции вероятности работоспособных состояний ОД “Nakoplenie”, предназначенная для оценки вероятности работоспособных состояний устройств судовых ЭСА и построения поверхности вероятностей в пространстве диагностических признаков, которая может быть использована при диагностировании сложной судовой аппаратуры.
Программа “Nakoplenie” выполняет следующие функции: накопление в численной форме функции вероятности работоспособных состояний с использованием метода Монте – Карло при заданной модели и выбранном пространстве наблюдаемых диагностических признаков; отображение в виде поверхности, полученной функции вероятности работоспособных состояний.
Далее рассмотрим построение апостериорной модели деградации системы с использованием изоварной модели диагностирования.