Теория надежности: марковские модели отказов и структурное резервирование систем

Базовым понятием теории является функция надежности P(t), которая определяет вероятность того, что система проработает без единого отказа от момента включения до времени t. Для большинства электронных и механических компонентов интенсивность отказов (лямбда) описывается так называемой U-образной кривой (Бадьевой кривой). Она включает период приработки (ранние отказы из-за заводского брака), период нормальной эксплуатации (где интенсивность отказов строго постоянна и имеет пуассоновскую природу) и период старения (где риск поломки экспоненциально возрастает из-за износа материалов). В классическом исследовании операций фокус делается на периоде нормальной эксплуатации, что позволяет использовать аппарат марковских цепей с непрерывным временем для расчета глобальной устойчивости системы.

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

Резервирование бывает нескольких типов: нагруженное (горячее), когда резервный элемент работает параллельно с основным и изнашивается с той же скоростью; и ненагруженное (холодное), когда резервный элемент выключен и включается только в момент поломки основного. Расчет надежности восстанавливаемых систем (где вышедшие из строя элементы могут быть заменены или отремонтированы без остановки всего комплекса) требует построения марковских графов гибели и размножения. Переходы в сторону ухудшения состояния идут с интенсивностью отказов (лямбда), а переходы в сторону восстановления — с интенсивностью ремонта (мю). Решая систему дифференциальных уравнений Колмогорова, аналитики находят коэффициент готовности (Availability) — долю времени, в течение которой система способна выполнять свои функции (знаменитые пять девяток: 99.999%).

Венцом применения исследования операций в этой области является задача оптимального резервирования. Добавление резервных элементов бесконечно повышает надежность, но неизбежно увеличивает вес, габариты и финансовую стоимость системы (что критично, например, для спутников). Возникает классическая задача нелинейного целочисленного программирования: необходимо максимизировать глобальную надежность системы при строгих линейных ограничениях на общий бюджет, массу и энергопотребление. Эта сложнейшая комбинаторная задача виртуозно решается методами динамического программирования. Алгоритм Беллмана рекурсивно распределяет вес и деньги между различными подсистемами ракеты, безошибочно выявляя самые слабые звенья и направляя каждый дополнительный килограмм резервной массы туда, где он обеспечит максимальный математический прирост выживаемости миссии.