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

Диагностика и надежность автоматизированных систем. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

В первом случае для обеспечения требующейся мощности достаточно двух генераторов, следовательно, третий будет в резерве. Структурная схема первой системы из двух основных и одного резервного элементов представлена на рис. 19, а среднее время наработки до отказа в соответствии с(41) для двух основных (N=2) и одного резервного (М=1) генераторов равно

T1

cp

=

M + 1

=

1

.

 

 

 

 

N × l

 

l

 

 

 

 

Во

 

втором

случае обеспечивают требуемую мощность - че

тыре генератора, а два будут в резерве. Схема для расчета надежности дана на рис. 31.

Рис. 31. Схема для расчета надежности системы из шести генераторов

Среднее время наработки до отказа второй системы будет

равно

T2cp =

M + 1

=

3

.

 

 

 

N × l

4l

Как видно, среднее время наработки до отказа первой системы больше чем второй.

61

 

1

 

 

 

 

 

P1(t)

 

 

 

 

 

 

P2(t)

0.5

 

 

 

 

 

 

0

100

200

300

400

500

 

0

t

Рис. 32. Вероятности безотказной работы систем для λ=0,005

Графики вероятностей безотказной работы для=0,005λ (рис. 32) показывают, что вероятность безотказной работы второй системы на начальном этапе превышает вероятность первой и гамма-процентная наработка до отказа при больших гамма может быть больше у второй системы.

4.4. Расчет надежности систем с поэлементным резервированием при различной интенсивности отказов элементов

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

t

f1-2 (t )= ò f1 (t )× f 2 t(- t)dt .

0

Для экспоненциального распределения вероятности безотказной работы fi (t) = li exp(- li t) и

62

t

f1-2 (t )= òl1 exp(- l1t)l2 exp[- l2 (t - t)]dt = l1l2 exp(- l2 t )´

0

´ òt exp[- t(l1 - l2 )]dt =

l1l2

[exp(- l1t )- exp(- l2 t )].

 

0

l2 - l1

Вероятность отказа группы равна интегралу от плотности вероятности отказа группы или времени ее работы до отказа, вероятность безотказной работы равна

t

l1l2

é1

- exp(- l1t)

 

1- exp(- l2 t)ù

P1-2 (t )=1- òf1-2 (t d)t =1-

-

 

 

 

ê

 

 

 

 

 

 

ú

l

 

- l

 

l

 

l

 

0

 

2

1

ë

 

1

 

 

 

 

2

û

или

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

(t )=

 

[l

2

exp(- l

1

t)- l

1

exp(- l

2

t)].

(42)

 

 

 

 

1-2

 

l2

- l1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для большего количества элементов можно последовательно

 

выполнять свертку плотностей f1-2 (t) с f3 (t), полученное выра-

 

жение для f1-3 (t)

- с f 4 (t )

и т. д.

 

 

 

 

 

 

 

 

 

 

 

Можно также записать плотность распределения вероятности

 

отказа элемента в преобразованиях Лапласа

 

 

 

 

 

 

 

 

 

F

(s )= ¥l

i

e-lit e-st dt =

li

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

ò

 

 

 

 

 

 

 

 

 

 

 

s + li

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

затем перемножить изображения по Лапласу всех плотностей, т. к.

 

свертке оригиналов соответствует произведение изображений, разде-

 

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

 

нию оригинала, выполнить обратное преобразование Лапласа и вы-

 

честь полученную функцию вероятности отказа из единицы:

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

li

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F1-n (s )= Õ

 

 

 

 

 

,

 

 

 

 

 

 

(43)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1 s

+ li

 

 

 

 

 

 

 

 

 

 

 

Q1-n (s )=

1

 

 

n

li

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Õ

 

 

 

 

 

.

 

 

 

 

(44)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s i=1 s + li

 

 

 

 

 

 

 

 

Обратное

 

 

преобразование

 

 

 

Лапласа

для -

дробн

63

рациональной функцииY(s) = Bm (s) , где n и m - порядки мно-

An (s)

гочленов знаменателя и числителя, причем n>m, можно выполнить, используя формулу вычетов. При всех различных корнях sk уравнения An (s) = 0 она имеет вид:

 

 

 

 

 

n

 

Bm

(s k )

× esk ×t , где A 'n

 

dA n

 

 

 

 

 

y(t) = å

(s k ) =

 

s=sk .

 

 

 

 

 

 

 

 

 

 

 

 

 

k=1 A'n (s k )

 

 

 

 

ds

 

 

 

 

Для

 

 

n

изображения

по

 

Лапласу

 

вероятности

отказа

Q1-n (s )=

1

 

li

 

 

 

 

 

 

 

 

 

 

 

 

 

Õ

 

 

 

 

 

запишем

уравнение

 

для определения

 

 

 

 

 

 

 

 

 

 

 

 

 

s i=1 s + li

 

 

 

 

 

 

 

 

корней sk,

приравняв

нулю

знаменатель дробно-рациональной

 

функции:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sÕ(s + li )= 0 .

 

 

 

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Получим корни уравнения: s 0

= 0, s1

= -l1 , s 2 = -l 2,... . Они все

 

будут различными, так как различные все интенсивности отказов

 

li . Производная от знаменателя дробно-рациональной функции

 

по s равна:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d é n

 

 

 

ù

 

n

 

n

 

 

 

 

n

 

 

 

êsÕ(s + li )ú

= Õ

(s + li

)+ sÕ

(s + li )+ s

Õ(s + li )+ ...

 

 

ds

 

 

ê

 

 

 

ú

 

i=1

 

i=2

 

 

i=1,i¹2

 

 

 

ë i=1

 

û

 

 

 

 

 

+ s Õn (s + li )+ snÕ-1(s + li ).

 

 

 

 

 

 

 

 

i=1,i¹n-1

 

 

 

i=1

 

 

 

 

 

 

 

Заменяя s на sk и учитывая, что s0=0 и е0=1, получим выражение вероятности отказа

 

 

 

n

 

n

-lk t

Õli

 

Q1-n (t )= åe

i=1

+1.

 

 

n

 

k=1

 

- lk

Õ(li - lk )

 

 

 

 

i=1,i¹k

64

Разделив числитель и знаменатель наλk и вынося минус за знак суммы, можно переписать последнее выражение в виде

 

 

 

 

n

 

 

 

 

 

 

n

 

 

 

li

 

 

 

 

 

 

 

 

Q1-n (t )= -åe

-lk t

 

Õ

 

+1.

 

 

 

 

 

 

 

- lk

 

 

 

 

 

 

 

 

k =1

 

 

 

i=1,i¹k li

 

 

 

 

 

 

 

 

Вероятность безотказной работы из зависимости P(t)=1-Q(t)

 

будет равна

n

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

li

 

 

 

 

 

 

 

 

 

P1-n (t )= åe-lk t Õ

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k=1

 

 

i=1,i¹k

 

li - lk

 

 

 

 

 

 

 

 

Для двух элементов (n=2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

 

 

(t )=

e-l1tl2

+

e-l2 tl1

.

 

 

(45)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1-2

 

 

 

l2 - l1 l1 - l2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Как видно, выражения (42) и (45) совпадают.

 

 

 

 

 

Для трех элементов (n=3) вероятность безотказной работы

P1-3 (t )=

 

 

e

-l1t l

2

l

3

 

 

 

 

e-l2t l l

3

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

1

 

 

+

 

(l2

-l1 )(l3

-l1 )

(l1

-l2 )(l3

-l2 )

(46)

 

 

 

 

 

 

e

-l3t l l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

1

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(l -l

3

)(l

2

 

-l

3

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

до отказа контроллера104 ч, датчиков 2.104 ч, исполнительных механизмов: 2000, 1500 и 1200 ч соответственно. Найти среднее время наработки до отказа резервированной и нерезервированной систем.

Присвоим значения среднего времени наработки до отказа элементов и запишем выражения для интенсивностей отказов:

tk:=104

td:=2.104

t1:=2000

t2:=1500 t3:=1200 ,

lk :=

1

, ld :=

1

,

l1

:=

1

,

l2

:=

1

,

l3

:=

1

.

tk

 

 

 

 

 

 

td

 

 

t1

 

 

t2

 

 

t3

Вероятность безотказной работы для системы без резервирования:

65

P1 ( t)
e- lk×t ×e- ld×t ×e- l1×t
e- (lk+ld+l1)×t

.

Среднее время наработки до отказа нерезервированной системы:

T1cp :=

1

 

T1cp = 1.538´ 103

ч.

lk + ld + l1

 

 

 

Получим выражение для вероятности работоспособного состояния системы с использованием преобразований Лапласа, для этого используем выражения (43), (44) изображения по Лапласу плотности вероятности отказа резервированного исполнительного механизма и изображение вероятности отказа

3

li

 

 

Fim(s)

 

Fim( s) := Õ

Qim( s) :=

.

s + li

 

i = 1

 

s

 

 

 

 

 

Выполним обратное преобразование Лапласа для Qim(s)

и запишем выражение для вероятности безотказной работы - ре зервированного исполнительного механизма Pim(t):

Pim( t) := 1 -

é

(- .500e-3)×t

(- .667e-3)×t

(- .833e-3)×t

ù

( -10.) ×e

+ 15.×e

- 6.×e

+ 1.

 

ë

 

 

 

û

Вероятность безотказной работы системы равна:

P2 ( t) := e- lk×t ×e- ld×t ×Pim( t),

а среднее время наработки ее до отказа

ó¥

T2cp := ô P2 ( t) dt, T1cp = 1538.46154 .

õ0

Аналогичный результат получим при использовании фор-

мулы (46).

66

Задача 9. Система управления содержит резервированный контроллер с нагруженным резервом, пять основных датчиков температуры и к ним два резервных ненагруженных, заменяющих любой из основных, три исполнительных устройства. Интенсивности отказов элементов заданы в табл.9. Найти среднее время наработки системы до отказа.

 

 

 

 

 

 

 

Таблица 9

λk.10-4

λd.10-4

λi.10-4

λk.10-4

λd.10-4

λi.10-4

вар

вар.

.

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0,10

0,30

0,50

14

0,24

0,44

0,36

1

0,11

0,31

0,49

15

0,25

0,45

0,35

2

0,12

0,32

0,48

16

0,26

0,46

0,34

3

0,13

0,33

0,47

17

0,27

0,47

0,33

4

0,14

0,34

0,46

18

0,28

0,48

0,32

5

0,15

0,35

0,45

19

0,29

0,49

0,31

6

0,16

0,36

0,44

20

0,30

0,50

0,30

7

0,17

0,37

0,43

21

0,31

0,51

0,29

8

0,18

0,38

0,42

22

0,32

0,52

0,28

9

0,19

0,39

0,41

23

0,33

0,53

0,27

10

0,20

0,40

0,40

24

0,34

0,54

0,26

11

0,21

0,41

0,39

25

0,35

0,55

0,25

12

0,22

0,42

0,38

26

0,36

0,56

0,24

13

0,23

0,43

0,37

27

0,37

0,57

0,23

5. РАСЧЕТ НАДЕЖНОСТИ ВОССТАНАВЛИВАЕМЫХ СИСТЕМ

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

67

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

5.1. Показатели надежности восстанавливаемых систем

В качестве характеристик надежности восстанавливаемых систем, так же как и для невосстанавливаемых, используется показатели безотказности [2] такие как: вероятность безотказной работы P(t) на интервале времени, вероятность отказа Q(t), среднее время наработки до первого отказа Тср, интенсивность отказов и плотность вероятности отказов.

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

Показатели потока отказов:

-ведущая функция потока отказовW(t), равная математическому ожиданию числа отказов за время t;

-параметр потока отказов (39)

w(t )= dW , dt

J I

w(t )= å ål j,i Pj (t ),

(47)

j=1i=1

 

где Pj (t) - вероятность работоспособного состояния, из которого возможен переход в неработоспособное; lj,i - интенсивность пе-

рехода из j- го работоспособного состояния вi - е неработоспособное; J, I – количество работоспособных и неработоспособных состояний объекта.

Показатели ремонтопригодности и восстанавливаемости:

-

вероятность восстановления объекта за заданное времяt,

равная

вероятности

,тогочто

время

t

 

 

 

 

восстановления

68

не превысит заданное

R (t )= P(t в £ t )

или статистическая оценка вероятности, равная отношению восстановленных за время t объектов nв к количеству объектов, поставленных на восстановление Nов:

~ ( )=

R t n в / Nов ;

- вероятность отсутствия восстановления объекта за задан-

ное время t

S(t )= 1 - R (t );

- плотность распределения вероятности восстановления

r(t )= dR ; dt

- интенсивность восстановления

m(t )= r(t) ; S(t )

-среднее время восстановления

¥¥

Tв = ò t × r(t )dt = òS(t )dt ;

00

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

¥

D в = ò(t - Tв )2 × r t(dt) .

0

Комплексные показатели надежности [2]:

-коэффициент готовности Кг(t), называемый в технической литературе функцией готовности, равный вероятности работоспособного состояния объекта в момент времениt; он учитывает возможность как сохранения безотказности, так и восстановления объекта путем ремонта;

-стационарный коэффициент готовности, равный вероятности работоспособного состояния объекта в произвольный -мо мент времени его использования и определяемый как предел коэффициента готовности (функции готовности)

Kг = lim Кг (t);

(48)

t®¥

 

69

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

K ог = K г (t)P(t з / t ),

где P(t з / t ) - условная вероятность сохранения работоспособного состояния на интервале t з , при условии, что к началу интервала объект был работоспособен.

Среднее время наработки на отказ на интервале(0, t) определяется как отношение среднего времени работоспособного -со стояния объекта за времяt к математическому ожиданию числа отказов за время t

t

ò K г (t )dt

q

cp

(t) =

0

.

(49)

t

 

 

 

 

ò w(t )dt

0

Для бесконечно большого интервала эксплуатации (при t ®¥)

 

K г

 

qcp =

lim w(t ).

(50)

t®¥

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

Показатели долговечности:

-ресурс – суммарная наработка объекта от начала его эксплуатации или ее возобновления после ремонта до момента перехода объекта в предельное состояние, при котором дальнейшее использование его либо невозможно, либо нецелесообразно;

-средний ресурс – математическое ожидание ресурса:

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

70

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]