Исследование структурной надежности методом статистического моделирования (стр. 6 из 14)

что соответствует

.

Разница между верхней и нижней оценками надежности теперь составляет:

, что позволяет прекратить расчеты, так как заданная точность H ± 0,01 достигнута. В качестве оценки надежности рассматриваемой сети принимаем среднеарифметическое с гарантией, что . При этом из полного множества, включающего девять сечений и восемь путей, нам удалось ограничиться рассмотрением всего трех сечений и четырех путей.

Рисунок 1.6. - Пример сети для двусторонней оценки надежности

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

1.2.4 Метод статистической оценки структурной надежности

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

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

Назовем разрезом подмножество элементов системы, удаление которых приводит к потере работоспособности. Рассмотрим некоторый разрез u, в который входит ровно z элементов. Частота выпадения такого разреза при статистических испытаниях стремится по ходу испытаний к ее вероятности :

.

Если обозначить через N общее число элементов сети, то вероятность Ри можно записать в виде:

.

Изменим исходные показатели надежности системы таким образом, чтобы каждый сомножитель

второго произведения увеличился в γ раз. Другими словами, вместо элемента с надежностью pj введем элемент с надежностью p'j такой, чтобы удовлетворялось условие:

(1.19);

Тогда надежность нового элемента будет определяться формулой:

.(1.20)

Если произвести преобразование (1.19) для всех элементов сети, то вероятность выпадения разреза u в процессе испытаний изменится и составит:

;

Введем коэффициент

увеличения частоты выпадания разреза

(1.21)

тогда, подставляя в (1.21) старое и новое значения частоты выпадания разреза u, получим:

.

Если в разрез u входит ровно z элементов, то:

,

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

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

Пересчет вероятности появления разреза u из преобразованной системы в исходную производится в соответствии с (1.21):

(1.22)

Для выполнения обратного преобразования (1.22) кроме факта отказа системы необходимо фиксировать и вес z соответствующего разреза.

Из (1.22) следует, что каждый случай появления разреза u с весом z в преобразованной системе соответствует

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

,

где

– число элементов, вышедших из строя при i-м отказе системы.

При возникновении очередного отказа m в преобразованной системе оценки надежности Рm исходной системы уточняются в соответствии с выражением:

,

где

– общее число просмотренных реализации состояний преобразованной системы.

Возникает вопрос, каким следует выбирать параметр преобразования γ для максимального убыстрения процесса статистических испытаний конкретной системы? Из (1.20) следует, что при γ = 1 изменение исходной надежности не происходит и убыстрение отсутствует. Если выбрать γ слишком большим, то в преобразованной системе будут в основном возникать разрезы большого веса, не характерные для исходной системы, причем их вклад в результирующую надежность при больших z в соответствии с (1.22) будет невелик. Поэтому параметр γ следует выбирать таким образом, чтобы максимизировать вероятность возникновения наиболее “ вероятных ” разрезов.

Рассмотрим этот вопрос более подробно. Для простоты положим, что показатели надежности всех элементов системы одинаковы и равны р. Обозначим через P(z) вероятность возникновения отказовых состояний веса z. Очевидно, что вероятность потери работоспособности системы тогда будет:

.(1.23)

Обычно для реальных систем значение

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

.

Для надежных систем в выражении (1.23) можно пренебречь всеми членами, кроме первого ненулевого, соответствующего минимальному разрезу веса

, который и будет наиболее вероятным. Таким образом, задача состоит в том, чтобы максимизировать коэффициент убыстрения для наиболее вероятного разреза, т. е. найти максимум функции при , , :

.(1.24)

Из условия

получаем оптимальное значение

.(1.25)

Подставляя выражение (1.25) в (1.20), нетрудно убедиться, что оптимальное значение

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

Убыстрение темпа набора статистики отказов в соответствии с (1.24), когда

, составит: