CC BY
65
УДК 519.72
РЕШЕНИЕ ЗАДАЧИ ОПТИМИЗАЦИИ НАДЕЖНОСТИ В ДУАЛЬНОЙ ВЫЧИСЛИТЕЛЬНОЙ СРЕДЕ
И.Л. Каширина
В статье рассматривается математическая модель многокритериальной задачи оптимального резервирования. Предлагаются метод ветвей и границ отыскания полной совокупности Парето-оптимальных решений данной задачи, а также генетический алгоритм, позволяющий получить представительную аппроксимацию множества Парето
Ключевые слова: оптимальное резервирование, многокритериальная задача, множество Парето, метод ветвей и границ, генетический алгоритм
Постановка задачи
Рассматривается сложная техническая система', состоящая из п элементов. Универсальным принципом обеспечения надёжности сложной системы является резервирование ее элементов, которое может быть аппаратным, информационным, временным, программным. Наибольшее распространение получили две схемы резервирования: резервирование замещением «один из двух» (1оо2 - 1 ои1^ 2) и резервирование по схеме мажоритарного голосования «два из трех» (2оо3). Системы без резервирования обозначаются 1оо1.
Предположим, что каждый элемент исследуемой системы (а, значит, и система в целом) характеризуется тремя основными параметрами: вероятностью безотказной работы р^, стоимостью Sj и среднем временем безотказной работы г^, У е {1,2,...,и} .Зависимость стоимости, вероятности безотказной работы и среднего времени работы до отказа от способа резервирования представлена в таблице.
Зависимость параметров системы от способа резервирования
Критерий 1оо1 1оо2 2оо3
Вероятность безотказной работы P 2P-P2 3p2-2p3
Стоимость S 2SG 4S
Среднее время ра-ботыдо отказа т 3Т/2 5Т/6
Через О обозначен коэффициент, увеличивающий стоимость схемы 1 оо2, если для данного резервируемого элемента существует надежный блок переключения на резерв. Если такой блок отсутствует, для элемента вводится запрет на возможность быть зарезервированным методом 1 оо2. Вводится ограничение на среднее время безотказной работы резервируемой системы: оно должно быть не меньше заданного значения.
Задача состоит в выборе способа резервирования для каждого элемента системы, который позволил бы получать наибольшую возможную вероятность безотказной работы при наименьшей стоимости системы с учетом выполнения указанного ограничения. Математическая модель бикритериальной задачи оптимального резервирования может быть
сформулирована в виде:
P = П(Fjxu + (2Pj - pj)xy + (3p) - 2pj)x3j ) ® max (1) j=i n
S = X (SjXij + 2SjGjX)j + 4SjX3j) ® min (2)
j=i
X1 j + X2 j + X3 j = 1 ,
Y (— X +—X +—X ) < — j=iTj Xlj + 3Tj X)J + 5Т: Xy) < Tc
j = 1, n (3)
(4)
Обозначения:
Г1, если элементу ] назначается способ Ху = < резервиров ания 1 ;
[ 0, иначе.
Здесь I е {1,2,3} - номер способа резервирования (соответственно, 1оо1, 1оо2, 2оо3). Р — вероятность работы без отказа всей системы, общая стоимость системы. Тс — минимальное среднее время работы до отказа всей системы. Ограничение (4) имеет соответствующий вид потому, что среднее время работы до отказа всей системы связано со средним временем наработки до отказа ее элементов
следующим соотношением: jl =
Y т
V j=
Каширина Ирина Леонидовна - ВГУ, канд. техн. наук, доцент, тел. 8 903 653 92 93, e-mail: kash.irina@mail.ru
Отметим, что с целью получения более простой однокритериальной модели задачи можно было бы отказаться от целевой функции (2), добавив соответствующее ограничение на итоговую стоимость системы. Однако, в данном случае существует несколько практических доводов в пользу рассмотрения именно бикритериальной модели. Во-первых, стоимость сложной системы является довольно высокой и может варьироваться в достаточно широких пределах. В связи с этим, имеет смысл отыскания всего множества парето-оптимальных решений бик-ритериальной задачи (или хотя бы его представительной аппроксимации), чтобы получить представление о возможных диапазонах изменения стоимости и надежности и иметь возможность выбрать в итоге компромиссный вариант. Во-вторых, задача конструирования такого рода систем является технически сложной и содержит ряд трудно формализуемых требований. Поэтому получение не одного, а сразу множества решений позволит выбрать из них наиболее приемлемое с конструкторской точки зрения.
Методы решения
Для отыскания решения данной задачи в рабо-
те разработано и исследовано с помощью вычислительно эксперимента два метода: алгоритм ветвей и границ, организованный таким образом, чтобы находить все множество Парето-оптимальных решений и генетический алгоритм, ориентированный на
получение некоторой представительной аппроксимации Парето-оптимального множества.
На рис. 1 приведена блок-схема метода ветвей и границ.
Рис. 1. Блок- сх
Рассмотрим подробнее основные этапы этого метода.
Ветвление.В данной задаче на каждом этапе используется следующее правило ветвления: Qк = Q1k и ^2 и ^3 , где Qк - это подмножество множества допустимых решений, в котором для резервирования элемента к выбран способ I, 1=1,2,3. Если способ резервирования 1оо2 для соответствующего элемента невозможен, то полагается Qк = Як и ^к .
Построение множества рекордных реше-ний.Множество рекордных решений для данной задачи представляет собой совокупность точек, которые являются недоминируемыми в множестве всех рассмотренных алгоритмом на данный момент вариантов. Процесс формирования такого множества следующий. Пусть в какой-то момент работы алгоритма имеется некоторое множество рекордных ре-
метода ветвей и границ
шений Rec и получено новое допустимое решение
new -г^ г» , new new
x . Если $x е Rec: x f x , то решение x в рекордное множество не добавляется. Если же в
new new
Rec нет решений, доминирующих x , то x включается в множество R, причем из R исключаются все решения, которые являются доминируемыми решением xnew. Следует отметить, что если
new new new
$xе Rec: x ~ x и x Ф x , то x также включается в множество R .
Оценка подмножеств допустимых решений и сокращение дерева поиска.На этом этапе для каждого подмножества решаются по две специально построенные в работе оценочные задачи, а также определяется, содержит ли данное подмножество допустимые точки (те, в которых выполняется ограничение на среднее время безотказной работы). Если допустимых точек нет, подмножество отбрасывается.
Стратегия построения и обхода дерева вариантов. Как показал вычислительный эксперимент, наиболее существенное сокращение перебора дает рассмотрение координат в порядке убывания значения отношения Sj / Pj . Предварительное упорядочивание элементов по убыванию коэффициента Sj / Pj и последующее ветвление дерева вариантов в
соответствии с полученной очередностью является мощным способом увеличения скорости поиска решения задачи повышения надежности систем прогнозирования граничных состояний с помощью предложенного метода ветвей и границ.
Несмотря на достаточно высокую эффектив-
ность предложенного алгоритма, его практическое использование ограничено для задач с размерностью более 30-40 резервируемых элементов. Поэтому для решения поставленной задачи в случае, если она имеет большую размерность, с помощью численного эксперимента были разработаны и исследованы различные модификации приближенных генетических алгоритмов и выбран наиболее эффективный алгоритм для решения данной задачи, схема которого представлена на рис. 2, позволяющий получать хорошую аппроксимацию множества Парето за приемлемые временные затраты [1].
Рис. 2. Блок-схема
Для этого алгоритма разработана специальная методика назначения приспособленности, включающая следующие шаги.
1) Для каждого решения в популяции вычисляется вектор целей (5, Р).
2) Из текущей популяции выбирается множество недоминируемых внутри этой популяции решений, они запоминаются и временно исключаются из рассмотрения.
3) Далее ищется множество недоминируемых
генетического алгоритма
решений в полученном усеченном множестве, и они тоже исключаются. Эта процедура проделывается до тех пор, пока все решения не будут исключены из популяции.
4) Затем все решения ранжируются: принадлежащие последнему исключённому множеству получают ранг 1, предпоследнему ранг 2. Решения, первыми исключенные из рассмотрения, получают максимальный ранг. Внутри каждого исключенного множества решения имеют одинаковый ранг.
5) Далее, в отдельности для каждой группы
решений одного ранга, происходит назначение скалярных оценок приспособленности. Предположим, ранг k имеет m решений. Тогда решение, сумма евклидовых расстояний от которого до остальных решений данного ранга максимальна, получит оценку k+(m-1)/m. Решение со второй по величине суммой расстояний получит оценку k+(m-2)/m. Решение с минимальной суммой расстояний будет иметь оценку к.Такой подход к оцениванию экземпляров популяции настраивает алгоритм не только на поиск недоминируемых решений (любая недоминируемая строка будет иметь оценку выше любой доминируемой), но и на поддержание разнообразия популяции (удаленные точки получают более высокие оценки), что в результате обеспечивает лучшую аппроксимацию Парето-оптимального множества.
С целью дальнейшего повышения надежности системы предлагается использовать дуальную вычислительную среду, в которую, кроме оптимизационной модели, включается модуль построения нейросетевых моделей некоторых элементов этой системы с целью осуществления дополнительного резервирования в случае, если по результатам решения оптимизационной задачи для данного элемента не предусмотрен резерв (например, элемент слишком дорогостоящий). То есть предлагается в качестве дополнительного резервного канала рассматривать вместо физической реализации некоторый виртуальный образ в виде нейросетевой модели, имитирующей преобразование входных воздействий от других компонентов в выходные величины, поскольку именно нейросетевая модель способна решать слабо формализованные задачи, оперативно аппроксимироватьпроизвольные непрерывные многопараметрические зависимости, обрабатывать данные, представленные в разнотипных шкалах.
Вычислительный эксперимент
Результаты численных экспериментов подтвердили эффективность предлагаемых методов и целесообразность их объединения в дуальную вычислительную среду.
В частности, для алгоритма метода ветвей и границбыл организован численный эксперимент с целью определения наилучшей стратегии обхода дерева вариантов. Анализировались стратегии, основанные на упорядочивании резервируемых элементов по убыванию характеристик, вычисляемых на основе значений параметров pj, и tj . Коэффициенты целевых функций определялись с помощью специальной рандомизированной процедуры. Рассматривались системы, содержащие от 15 до 50 элементов с возможностью резервирования. Упорядочивание элементов по убыванию характеристик Sj / Pj дало наилучшие результаты, при этом время работы алгоритма сокращается в несколько раз.
Если бы решение задачи осуществлялось методом полного перебора, то, например, в каждой задаче размерности 15 необходимо было рассмотреть 3 5 степени вариантов, то есть перебрать 14348907 решений. По результатам численного эксперимента, метод ветвей и границ, не использующий стратегию упорядочивания элементов, для этой размерности просматривает в среднем 40 000 решений, а алгоритм, использующий эту стратегию, просматривает в среднем 8000 решений (что, примерно, в 1800 раз меньше, чем метод полного перебора), при этом средний объем Парето-оптимального множества, составил 65 решений.
С помощью построенной многокритериальной оптимизационной модели и разработанного для нее алгоритма решена задача повышения надежности управляющих трактов стендовой информационно-управляющей системы для диагностики и испытаний жидкостных ракетных двигателей(ЖРД).
Для элементов, зарезервированных по схеме 1оо2, в рамках дуальной вычислительной среды далее были построены нейросетевые модели и сформированы дополнительные виртуальные каналы измерений, позволяющие при физическом резервировании по дублированной схеме 1 оо2 использовать более наглядный метод мажоритарного контроля 2оо3[2-3]. Возможность использования нейро-сетевого резервирования в системах управления и регулирования основных параметров ЖРД подтверждена на более чем 60-и огневых испытаниях в ОАО КБХА.
Литература
1. Львович Я.Е. Генетический алгоритм решения многокритериальной задачи повышения надежности резервирования/Я.Е.Львович, И.Л. Каширина, А.А. Тузи-ков//Информационные технологии, 2012.-№6- С.56-60.
2. Каширина, И.Л.Нейросетевая модель датчика давления в камере сгорания для наземных огневых испытаний ЖРД[Текст] / И.Л. Каширина, Я.Е Львович, А.А. Тузиков // Вестник Воронежского государственного технического университета. - 2010. - Т. 6. -№ 11. - С. 4-7.
3. Каширина И.Л. Нейросетевое резервирование дублированных измерений параметров при наземных огневых испытаниях ЖРД/ И.Л. Каширина, Я.Е. Львович, А.А. Тузиков// Информационные технологии. 2011. № 9. С. 74-78.
4. Львович, Я. Е. Оптимизация построения решающих правил при управлении испытаниями [Текст] / Я. Е.Львович, И. Л. Каширина, А. А. Шостак // Вестник Воронежского государственного технического университета. - 2012.-Т. 8. - № 5.-С. 22-24.
5. Подвальный, С.Л. Многоальтернативные системы: обзор и классификация [Текст] / С.Л Подвальный // Системы управления иинформационные технологии. - 2012. - Т. 48. - № 2. - С. 4-13.
6. Фролов В.Н., Львович Я.Е., Повальный С.Л. Проблема оптимального выбора в прикладных задачах. Воронеж: ВГУ. 1980. - 140 с.
Воронежский государственный университет
SOLUTION OF THE TASK OF OPTIMIZATION OF RELIABILITY IN THE DUAL
COMPUTING ENVIRONMENT I.L. Kashirina
In the article the mathematical model of a multicriteria problem of optimal reservation. Is offered the method of branches and borders find a full set of Pareto-optimal solutions, as well as the genetic algorithm to obtain a representative approximation of the Pareto set
Key words: optimal redundancy, multicriteria problem, Pareto set, the method of branches and borders, genetic algorithm
