Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

511_SHerstneva_O._G._Modelirovanie_funktsionirovanija_ehlementov__

.pdf
Скачиваний:
4
Добавлен:
12.11.2022
Размер:
1.48 Mб
Скачать

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

Используя, результаты проведенных исследований и сопоставив их с наблюдаемыми характеристиками, получим следующие системы уравнений:

Для 1-ого варианта разбиения:

РППВ

(1 )

 

РППН

(1 )

 

 

РИВ

 

 

 

2

(1 )(1 )

Н

 

 

 

РИ 2

 

Решение системы уравнений (28) имеет вид:

1

 

РППн

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

в

 

в

Ри 2

РПП

 

в

 

 

 

 

 

 

РПП

 

 

 

1

 

 

 

 

 

 

 

 

 

(28)

(29)

Для 2-го варианта разбиения система уравнений (граф наблюдаемых событий тот же) имеет вид:

РППв

РППн

РПв

РПн

 

(1 )

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

( 30)

1

 

 

 

 

 

 

 

 

 

 

 

 

(1 )

 

 

1

 

 

 

 

 

 

 

 

Решение системы уравнений (30) приведено в виде, удобном для практического применения:

 

 

 

 

 

 

 

 

 

 

 

в

 

 

 

1 РП

 

 

 

 

 

н

 

 

 

 

 

 

 

РП

 

 

 

 

(31)

 

 

 

 

 

 

 

РППн

 

 

 

 

 

 

 

1

 

н

 

 

 

РПП

 

 

 

 

 

 

 

Р

н

 

 

 

 

 

 

ПП

 

51

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

Таким образом, на примере приведенных выше исследований системы сигнализации ОКС№7, показано применение разработанного метода расчета реальных показателей надежности функционирования участка сигнальной сети с учетом данных подсистемы контроля и диагностики систем сигнализации.

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

Разработаны математические модели и получены аналитические формулы для определения показателей надежности, к которым относятся:

-- вероятность трансформации сигнала подтверждения в сигнал запроса (ошибка контроля I рода);

-- вероятность трансформации сигнала запроса в сигнал подтвер-

ждения (ошибка контроля II рода);

- - вероятность искажения сигнальной единицы.

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

4. Методика сбора и обработки статистических данных системы контроля и управления.

4.1. Постановка задачи сбора и обработки статистических данных

Методика сбора и обработки статистических данных о надежности сетевых элементов приведена при условии, что сетевой элемент используется непрерывно в течении всего периода наблюдений в условиях нормальной эксплуатации с регистрацией параметров его функционирования [21, 22, 23, 24, 25, 26].

Экспериментальному определению подлежат следующие показатели, оценивающие надежность:

1)средняя наработка на отказ;

2)среднее время восстановления;

3)интенсивности событий, являющиеся параметрами экспоненциального распределения: параметр потока отказов; интенсивность восстановления; интенсивность проведения периодического контроля;

4)вероятности событий, вычисляемые как соотношения чисел соответ-

52

ствующих событий.

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

4.2. Оценка средней наработки на отказ

Средняя наработка на отказ - это отношение наработки восстанавливаемого объекта к математическому ожиданию числа его отказов в течении этой наработки.

 

T

t

 

 

,

(32)

 

M m(t

 

)

 

 

 

 

 

 

где t - суммарная наработка объекта;

 

m(t )

- число отказов на наработку t ;

 

М – знак математического ожидания.

Статистически средняя наработка на отказ одного объекта определяется отношением суммарной наработки объекта к наблюдаемому числу отказов:

ˆ

 

t

 

 

I

 

 

 

 

 

 

 

T

m(t )

 

m(t )

 

, (33)

 

 

 

 

 

 

 

 

 

 

 

m(t ) tK

 

 

 

 

 

 

 

K 1

 

 

где m(t ) - наблюдаемое число отказов за суммарную наработку; tK - наработка на к-ый отказ;

- наработка объекта от момента устранения последнего отказа до окончания наблюдения за объектом.

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

Пусть имеются N однотипных объектов, которые непрерывно используются по назначению и наблюдаются в течении времени t. Отказавший объект мгновенно заменяется на работоспособный. Фиксируются моменты возникновения отказов. Исходные статистические данные представляют собой не полностью определенную выборку. Оценка проводится по плану испытаний NRT, согласно которому одновременно испытывают N объектов. Отказавшие во время испытаний объекты заменяются новыми, испытания прекращают по истечении наработки t.

53

ˆ

 

NT

,

(34)

T

m

где t – продолжительность наблюдения, совпадающая с наработкой на отказ одного объекта;

N – число однотипных объектов;

m – суммарное число отказов за время наблюдения.

Точечная оценка средней наработке на отказ производится согласно ISO 2836-80:

Оценка интенсивности отказов при экспоненциальном определении наработки на отказ:

ˆ

I

 

m

.

(35)

 

 

 

T

 

Nt

 

Интервальная оценка может быть выполнена в следующем варианте. Нижняя и верхняя доверительные границы для средней наработки на отказ при экспоненциальном распределении наработки на отказ определяются согласно [СТ СЭВ 2836-80].

 

 

2Nt

 

 

 

ˆ

TH

 

 

 

2mT

 

 

 

 

 

; (36)

 

2

 

2

 

 

 

1 2,2m 2

 

 

1

2,2m 2

 

2Nt

 

ˆ

 

 

TB

 

2mT

,

(37)

2

2

 

 

 

 

 

2,2m

 

2,2m

 

 

где m – число отказов за время наблюдения; 1-α – доверительная вероятность;

2p, k - квантиль 2 - распределения уровня Р с числом степеней свободы,

равным К.

Таблица квантилей 2 - распределения приведена в [9].

По найденным доверительным границам находятся доверительные границы для интенсивности отказов:

 

 

 

 

 

 

I

 

; (38)

 

H

TB

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

.(39)

 

B

TH

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Кроме интервальной оценки можно найти предельную относительную

ˆ

 

 

 

 

 

 

 

 

ошибку оценки T

зависящую от объема наблюдений (N и Т) и от принятой

доверительной вероятности при экспоненциальном распределении наработки на отказ [6].

54

4.3. Оценка среднего времени восстановления.

Среднее время восстановления - это математическое ожидание восстановления работоспособного состояния [II].

ˆ tgG(t) tg(t)dt , (40)

0 0

где G(t) - функция распределения времени восстановления; g(t) - плотность распределения времени восстановления.

Статистическое определение среднего времени восстановления объекта:

 

I

n

 

ˆ

i ,

(41)

 

 

n i 1

 

где n - число восстановлений объекта;

i - восстановления объекта после i-го отказа.

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

 

I

n

 

ˆ

i

(42)

 

 

n i 1

 

где n - число восстановлений;

i - продолжительность i-го восстановления.

При экспоненциальном распределении времени восстановления оценка интенсивности восстановления:

 

I

 

 

.

(43)

 

Интервальная оценка выполняется следующим образом. Нижняя и верхняя доверительные границы для среднего времени восстановления определяются:

H

 

2n ˆ

 

 

, (44)

2

 

 

 

 

1 / 2,2n

 

55

 

 

 

2n ˆ

 

B

2

. (45)

 

 

 

 

 

 

/ 2,2n

По найденным доверительным границам находятся доверительные границы для интенсивностей восстановления:

H

 

 

I

;

(46)

B

 

 

 

 

 

B

 

 

I

.

(47)

 

 

 

 

H

 

 

 

По объему наблюдений n и доверительной вероятности при экспоненциальном распределении времени восстановления может быть найдена предельная относительная ошибка оценки [6].

4.4. Оценка вероятностей состояний.

В соответствии с применяемыми данными оценке подлежат вероятности попадания в состояния восстановления:

-вероятность того, что устройство попадет в состояние восстановления

сотказом, обнаруживаемым непрерывным контролем РН;

-вероятность того, что устройство попадет в состояние восстановления

сотказом, приводящим к продолжительному неработоспособному состоянию

РП.

-вероятность того, что устройство попадет в состояние восстановления

сдвумя видами отказа одновременно РНП;

-вероятность того, что устройство попадет на восстановление в работо-

способном состоянии РР Кроме этого, оценке подлежат условные вероятности ошибок контроля

I и II рода и числовые характеристики случайных величин (наработки до отказа, времени восстановления и др.). В [10, 11] приведена формула для оценки вероятности диагностирования, которая вычисляется через соотношения соответствующих чисел испытания. Точечная и интервальная оценка производятся по известным формулам. Например, интервальная оценка вероятности события при 0,1<P <0,9 производится по формулам:

H

[ t 2

/ 2n t

v(1 v) / n t 2 / 4n 2

 

/1 t 2 / n ; (48)

 

 

 

 

 

 

B

[ t 2

/ 2n t

v(1 v) / n t 2 / 4n2

/1 t 2 / n , (49)

56

где v

m

-точечная оценка вероятности, nv 4,

n (1 v) n .

n

 

 

 

 

 

При nv 10,

n (1 v) 10 формулы интервальной оценки (48, 49) упрощают-

ся.

 

 

 

 

 

 

 

 

 

 

 

 

 

vH v t

v(1 v) / n, (50)

 

 

 

 

 

 

 

 

 

vH v t

v(1 v) / n, (51)

 

Коэффициент t берется из таблиц [10, 11]. Наиболее общая схема об-

работки статистических данных для оценки показателей надежности включает в себя следующие этапы:

-выбор предположения о виде распределения наработки или времени восстановления;

-анализ однородности статистического материала;

-оценка параметров распределения;

-проверка согласия эмпирического распределения с заданным теоретическим распределением;

-оценка показателей надежности.

В большинстве случаев при расчетах надежности коммутационных узлов и станций принято, что наработка до отказа, время восстановления и случайное время между проверками распределены по экспоненциальному закону. Проверка согласия эмпирического и теоретического распределений производится по [8, 13, 12].

4.5. Оценка показателей надежности

Исходными данными для оценки средней наработки на отказ в течении одного квартала являются:

-продолжительность испытаний t = 2160час.

-доверительная вероятность интервальной оценки 0,95 ;

-число однотипных объектов, находящихся в эксплуатации, N=2, 5, 10,

16, 250, 256, 500.

-числе отказов m однотипных объектов в течении квартала.

Формулы для нижней и верхней доверительной границы следуют из формул интервальной оценки средней наработки на отказ среднего восстановления:

TH TB

H

B

ˆ

aH T , (52)

ˆ

aBT ;(53)

aH ˆ, (54)

aB ˆ;(55)

57

где и - коэффициенты для нижней и верхней границы соответственно; эти коэффициенты вычисляются по формулам (для 0,95 ):

aH 2m / 20,975;2m ;

(56)

aВ 2m / 20,025;2m .

(57)

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

ˆ

 

T Nt / m; (58)

 

 

1

m

 

ˆ

i , (59)

 

 

 

 

m i 1

 

где m – число отказов для оценки Т и число восстановлений для оценки

m I .

 

В таблице 1 приложения А, приведены квантили 2

распределения, взя-

тые в [8, 12,13], а также значения коэффициентов и

, вычисленные по

приведенным формулам.

В таблице 2 приложения А приведена точечная оценка средней наработки на отказ по результатам испытаний за квартал.

Значения оценки округлены до целых чисел. Если испытания проводятся в течении шести месяцев или года, то m -число отказов за этот период испытаний; значения точечной оценки следует увеличить в 2 и в 4 раза соответственно.

4.5.1. Текущая оценка средней наработки на отказ и среднего времени восстановления.

Для решения задачи вычисления точечной оценки путем корректировки текущего значения запишем оценку средней наработки на отказ в виде:

ˆ

m

N `tm

 

T

 

, (60)

m

 

 

 

где m - суммарное число отказов;

m - момент последнего, m-го отказа, отсчитываемый от начала наблю-

дения. Вообще говоря, tm- это любой момент времени, отсчитываемый от начала наблюдения, между m-м и (m+1)-м отказами.

В дальнейшем точечную оценку будем называть текущей оценкой. Вы-

разим ее через предыдущую оценку ˆ . Путем несложных преобразо-

T m I

ваний текущая оценка может быть представлена в виде:

58

n

 

ˆ

m T m 1 T m , T 0 0 ;

(61)

 

T

 

 

ˆ

m N m

m 1 T m 1 / m ,

(62)

 

T

где

ˆ

 

 

 

 

T (m) - корректирующее слагаемое за счет m-го отказа (приращение

оценки на m-м шаге).

Аналогично, оценка среднего времени восстановления:

n

 

ˆ(n) 1/ n i

(63)

i 1

 

может быть представлена в виде:

 

(n) (n 1) (n),

(0) 0 , (64)

ˆ

 

(n) [ n )n 1)]/ n,

(65)

ˆ

 

где ˆ(n) - корректирующее слагаемое за счет n-го восстановления (или приращение оценки на n-м шаге).

Видно, что с увеличением объема выборки (как m, так и n) корректи-

ˆ

ˆ

и T (m) имеют тенденцию к уменьшению.

рующие слагаемые (n)

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

Таким образом, для вычисления текущей оценки средней наработки на отказ необходимы следующие сведения:

I) число однотипных объектов N;

2)значение предыдущей оценки Т(m-1);

3)текущий номер отказа m;

4)момент последнего отказа tm;

5)момент предыдущего отказа tm-1 ;

Замечание. Момент времени tm-1 и tm отсчитываются от нуля, являющегося началом периода испытаний.

Для вычисления текущей оценки среднего времени восстановления необходимы следующие сведения:

1)значение предыдущей оценки ˆ(n 1) ;

2)текущий номер восстановления n;

3)продолжительность последнего восстановления n

4.5.2. Объединение точечных оценок средней наработки на отказ и среднего времени восстановления

59

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

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

Т1(m1) и Т2(m2),

где m1 и m2 - количество отказов на первом и втором периоде соответственно.

После несложных преобразований оценка для объединения этих двух периодов будет иметь вид:

ˆ

m1

ˆ

ˆ

m2 . (66)

T

m2 [m1T (m1 ) m2T (m2 )] / m1

Аналогично определяются оценки среднего времени восстановления для двух периодов эксплуатации:

ˆ n1 n2 [n1 ˆ1 (n1 ) n2 ˆ2 (n2 )]/[n1 n2 ] . (66)

где ˆ1 (n1 ) и ˆ2 (n2 ) - оценки среднего времени восстановления для двух периодов эксплуатации;

n1 и n2 количество восстановлений на первом и втором периоде соответственно.

Видно, что объединенная оценка является линейной комбинацией объединяемых оценок, при этом коэффициенты линейной комбинации пропорциональны числам событий (числам отказов и числам восстановлений соответственно).

4.6. Оценка интенсивности событий для случая полностью и неполностью определенной выборки.

Несмещенная оценка x для интенсивности ω, а также нижняя и верхняя ( xH и xB соответственно) доверительные границы интенсивности при односторонней доверительной вероятности вычисляются по следующим формулам:

 

n

 

 

x

 

 

 

 

x

 

 

для n>1

x [n 1] / xi ; xH

 

,

xB

 

,

(67)

 

 

 

i 1

 

 

r5

 

 

 

r4

 

n=1

xH

x

, xB

 

x

,

 

 

(68)

 

 

 

 

 

 

r1

 

 

 

r2

 

 

 

60