Применение метода исключения варьируемого параметра при решении задач диагностирования
Портнягин Н. Н., Пюкке Г. А.,
Как было отмечено выше, в пространстве диагностирования, при различных сочетаниях вариаций параметров компонент gi наблюдается статистическая устойчивость распределения вероятностей конечных положений точки состояния системы. Такая закономерность обусловлена наложением в каждой точке плоскости К1, К2 на точки данной кратности вариации точек вариации более высокого порядка кратности. Это дает возможность построения в пространстве диагностирования совокупности областей работоспособности ОД с заданным порогом вероятности.
Для построения функции вероятности работоспособности ОД(Рис.3.7) всю совокупность точек плоскости К1, К2 разделяем на два непересекающихся подмножества работоспособных и неработоспособных состояний. Для выполнения процедуры деления используются условия работоспособности
gi min ≥ 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. Для каждой точки пространства диагностических признаков (К1, К2) будем определять P(К1, К2)=P(A1/Bк) вероятность нахождения ОД в работоспособном состоянии и вероятность нахождения ОД в неработоспособном состоянии:
Pн(К1, К2) = .
Построения будем выполнять методом статистических испытаний с использованием пакетов Visual Basic и EXEL. Модель (2. 3) и генератор случайных чисел реализованы в пакете Visual Basic, построение поверхности идентификации выполнены в пакете EXEL. Испытания будем выполнять по схеме Бернулли. Результаты по каждому из n испытаний будут включать следующие массивы: значения координат на плоскости диагностических признаков [K1]n , [K2]n ; значения компонент [gi ]n; значения функции принадлежности испытания к области работоспособности ([F]n = 1, если система gi min ≥gi ≥ gi max не нарушена и [F]n = 0, если нарушена). Результаты по каждому испытанию накапливаются в два двумерных массива работоспособных и неработоспособных состояний. По столбцам массивов откладываются дискретные значения [K1]n , а по строкам значения [K2]n . Далее вычисляются численные значения третьего массива [P{K1, K2}]вероятности работоспособности объекта диагностирования. По завершении выполнения алгоритма будут сформированы два массива: массив количества событий, соответствующих работоспособности ОД [RR{K1, K2}] и массив количества событий, соответствующих неработоспособности ОД [RN(К1, К2)]. Вычисление массива вероятности работоспособности ОД [P{K1, K2}] будем выполнять по соотношению условных вероятностей: для каждой точки пространства {K1, K2}.(Рис. 3.8.)
.
Алгоритм приведен на рисунке 3.7, где блоки 1 - 5 генерируют диагностико-регулировочную модель (2.3), блоки 6 - 10 моделируют генератор случайных чисел, блоки 6 - 10 выполняют процедуры накопления вероятностей и построения поверхности идентификации. Построив поверхность функции вероятности работоспособности и определив на ней линии уровней с заданными вероятностями Р = 0.95; 0.9; 0.8;……, получим границы области работоспособности ОД с заданными порогами вероятности. Функция вероятности P{K1, K2}, представленная на рисунке 3.8, построена в результате опробования разработанного алгоритма на базе рассмотренного в работе [8] примера диагностической модели мостового выпрямителя, часто используемого в схемах судовой автоматики.
Разработанный алгоритм решает поставленную задачу для случая возникновения дефектов максимально возможной для данного ОД кратности в вероятностной ее постановке. На основе полученных результатов можно сделать вывод об эффективности предлагаемого метода при рассмотрении решаемой задачи. Приведенная на рис. 3.9. поверхность получена при количестве испытаний N = 103. Дальнейшее увеличение N не приводит к заметным изменениям результата, что свидетельствует о сходимости алгоритма. Затраты времени при реализации алгоритма для рассматриваемого примера составляют 15 минут, при использовании процессора Intel Pentium c тактовой частотой 166 МГц.
Вид функции P(К1, К2) будет определяться топологией ОД и параметрами составляющих компонент. Построение области работоспособности сводится к нахождению в каждой точке значения вероятности нахождения объекта в работоспособном состоянии. Определив таким образом функциюР(К1, К2) и задавшись требуемым пороговым значением величины вероятности, получим область работоспособности объекта в пространстве основных диагностических признаков К1, К2.
В работе [102] выполнено опробование разработанного алгоритма на примере диагностической модели мостового трехфазного выпрямителя, часто использующегося в схемах возбуждения судовых синхронных генераторов и других цепях схем судовой автоматики, там же рассмотрен вопрос построения области работоспособности. Разработанный алгоритм позволяет решать задачу для дефектов любой кратности, что особенно важно при диагностировании схем судовой автоматики, содержащих многополюсные элементы. Принятый вероятностный подход оценки состояния ОД позволяет решать задачи локализации множественных дефектов, а также решать задачи оценки и регулирования запаса работоспособности судовой аппаратуры. Эффективность метода достаточно высока при решении основных задач диагностирования.
В КамчатГТУ разработан автоматизированный программно-аппаратный комплекс решения задач диагностики, на основе которого поставлен мониторинг технического состояния ОД. Полученные экспериментальные данные используются для решения задачи определения предаварийных состояний ОД на основе метода исключения варьируемого параметра. Процедура определения предаварийных состояний ОД выполняется на основе анализа компонент вектора условных вероятностей состояний ОД, полученного на основе метода статистических испытаний. Каждая компонента вектора определяется по формуле статистической вероятности:
Рn = Wn /W, где W - количество вариаций параметров компонент различной кратности, приведшие в выбранную точку пространства диагностических признаков, Wn - количество вариаций параметров компонент кратности n, приведшие в эту же точку пространства диагностических признаков. Количество компонент вектора условных вероятностей составляет 2m , где m - количество компонент, составляющих ОД. Современный персональный компьютер позволяет выполнить прямой расчет вектора Р, в каждой измеренной точке пространства диагностических признаков К1, К2, методом статистических испытаний. Если значения компонент вектора Р определены для каждой точки наблюдения, то может быть определена величины запаса работоспособности по логарифмическому критерию: R = lg (P0 / max Pn), где n = 1, 2m . Здесь: Р0 - функция вероятности работоспособных состояний ОД, с максимумом в центре изоварной картины. Функция max Pn характеризует расположение максимума вероятности возникновения дефекта заданной кратности в выбранной точке. Критерий R работает следующим образом: при R > 1, вероятность работоспособных состояний ОД, на порядок превышает вероятности других событий и процедура определения предаварийных состояний ОД малоэффективна, т. е. ОД находится в работоспособном состоянии с большим запасом работоспособности. При значениях 0< R <1 ОД находится в работоспособном состоянии, но запас работоспособности невелик, и по условию: Nd = n
max Pn можно определить номер наиболее вероятного дефекта. При значениях R<0 ОД неработоспособен и находится в состоянии дефекта, процедура определения предаварийных состояний ОД нецелесообразна, так как необходимо решать задачу локализации и определения дефектных элементов. На рис. 3.9 приведена экранная форма выполнения программы построения области предаварийного состояния электронного стабилизатора.
.
Таким образом, область, в которой процедура определения предаварийных состояний ОД эффективна, задается условием 0< R <1 . Это условие выполняется вблизи границы области работоспособности и выделяет переходную зону на карте диагностирования.