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

книги из ГПНТБ / Оптимизация процессов грузовой работы

..pdf
Скачиваний:
10
Добавлен:
27.10.2023
Размер:
26.66 Mб
Скачать

ГЛАВА 3

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

При старении и износе интенсивность отказов является возра­ стающей функцией времени.

Условию возрастания интенсивности отказов отвечают закон Релея, нормальный закон и частные случаи законов Вейбулла и

гамма-распределения. В последнее

 

 

время все чаще используется нор­

 

 

мальный закон распределения дли­

 

 

тельности безотказной работы. Это

 

 

объясняется

тем,

что

большинство

 

 

систем проходит тренировку на за-

 

 

воде-изготовителе

и поступает в эк с­

 

 

плуатацию

без

явно

 

выраженного

 

 

участка приработки. При установив­

 

 

шемся режиме работы поток отка­

 

 

зов

таких систем

близок

к простей­

 

 

шему.

Отклонения от

него обнару­

 

 

ж иваю тся только

по истечении опре­

 

 

деленного

промежутка

 

времени

в

 

 

связи

с появлением

отказов из-за

Рис. 3.9. Функция готовности си­

износа или

старения

элементов.

Т а ­

стемы

с нормальным распределе­

кая

модель

описывается

нормаль­

нием времени безотказной работы.

ным законом распределения. Кроме

 

 

того,

нормальный

закон

удобно использовать при анализе надеж ­

ности

слож ны х

систем,

для которых характерен уход параметров

элементов за допустимые пределы.

 

 

 

Н а

рис. 3 .9

показано поведение функции

готовности двух систем

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

время

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

одни и те

ж е, а среднеквадратичные

откло­

нения

различны. Как

видно из

рисунка, системы с малыми

значе­

ниями среднеквадратичных отклонений имеют более выраженный провал' функции готовности, что является результатом концентра­ ции плотности отказов в небольших интервалах [кривая Г х (£)].

Усистем с большим значением среднеквадратичного отклонения

понижение

функции готовности выражено менее

четко [кривая

Г 2 (і)], так

как отказы распределены на больших

промежутках и

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

Интенсивность отказа является возрастающ ей функцией времени в случае гамма-распределения при k > 1 (см. табл. 1.4). Следует отметить, что гамма-распределение при k > 1 наиболее удобно для описания процессов старения и износа. Этому ж е закону подчинено

85

АНАЛИЗ ПОВЕДЕНИЯ ФУНКЦИИ ГОТОВНОСТИ

распределение времени безотказной работы избыточных систем, для которых характерна схема накапливаю щ ихся повреждений, а такж е резервированных систем с включением резерва по способу замещ е­ ния при условии, что потоки отказов основной и резервных систем простейшие [3 3 ]. И, наконец, простота преобразования Л апласа плотности распределения этого закона в сочетании с возможностью ис­ пользования его для представления эмпирических распределений обу­

словливаю т

широкое применение гамма-распределения на практике.

Х арактер

поведения функции готовности в случае гамма-распре­

деления при

k = 2 и k = 3 показан на рис. 3 .10 [кривые Г 2 ( t)

14t)

ПО

 

 

 

Рис. j|3.11.

Функция

готовности

 

 

в случае распределения

времени

Рис. 3.10. Функция

готовности

безотказной

работы

по

закону

в случае Фамма-распределения вре­

 

Вейбулла.

 

 

мени безотказной

работы.

3 >

k2 >

К“ =

const.

и Г\ ( t) соответственно]. На начальном участке функция готовности близка к вероятности безотказной работы, а затем без заметного провала асимптотически приближается к установивш емуся значению.

В последнее время наметилась тенденция во многих случаях использовать не экспоненциальный закон распределения времени безотказной работы слож ны х систем, который иногда употребляется недостаточно обоснованно, а закон Вейбулла при k > 1 как наиболее точный и учитывающий процесс старения [2 2 ]. В частности, распре­ деление Вейбулла может быть использовано при ускоренных испы­ таниях элементов навигационной аппаратуры в форсированных режи­ мах. Кроме того, данное распределение с параметром k = 1,4 т-1,9 наблюдается у некоторых типов электронных ламп [3 3 ]. Распреде­ ление Вейбулла интересно еще и потому, что им можно аппроксими­ ровать, изменяя параметры, процессы приработки, нормальной

работы и старения аппаратуры . Поведение функции

готовности

при этом законе в зависимости от значений параметра

k при К =

= const показано на рис. 3 .11 .

Закон Релея удобен для описания распределения времени без­ отказной работы при интенсивном старении аппаратуры, когда отказы не удовлетворяю т условиям стационарного случайного про­ цесса. В частности, закону распределения, близкому к закону Релея, подчиняется время между отказами некоторых типов электровакуум ­ ных приборов [3 3 ].

86

ГЛАВА 3

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

отказов

с

учетом

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

от

значения

параметра а (см.

табл. 1.3) при экспоненциальном

а)

 

 

 

 

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

показана

на

рис.

 

 

 

 

3 .12 . Более глубокий, но менее

 

 

 

 

длительный

провал

функции

Г (і)

 

 

 

 

характерен для малых значений а.

 

 

 

 

 

Таким образом, из анализа пове­

 

 

 

 

дения функции готовности аппара­

 

 

 

 

туры при возрастающ ей со време­

 

 

 

 

нем интенсивности ее отказов мо­

 

 

 

 

жно

определить

последовательность

 

 

 

 

мероприятий,

которые

необходимо

 

 

 

 

осущ ествить

для

устранения

про­

 

 

 

 

вала функции готовности, что имеет

 

 

 

 

сущ ественное значение для

уникаль­

 

 

 

 

ных

судовых

систем,

эксплуатируе­

 

 

 

 

мых и в период старения.

 

 

 

 

 

 

 

ВЛИЯНИЕ ВОССТАНОВЛЕНИЯ

§ 3.5

 

 

 

 

НА ПОВЕДЕНИЕ ФУНКЦИИ

 

 

 

 

 

 

 

ГОТОВНОСТИ

 

 

 

 

 

 

 

 

 

 

 

 

При установивш емся режиме работы

Рис. 3.12. Функция готовности (а)

влияние

свойства

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

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

мости на готовность системы оце­

восстановления (б) в случае рас­

нивается

 

посредством

выражения

пределения

времени

безотказной

(1 .51). При этом необходимо лишь

 

 

работы по закону Релея.

 

 

 

Оз < o z .

знать

среднее

время

восстановле­

 

 

 

ния

Тв.

В

случае

 

неустановивш егося

режима

работы

системы

появляется необходимость в оценке зависимости поведения

функции

готовности

от

вида

закона

распределения

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

и его

параметров. Х арактер

распределения

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

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

ления времени восстановления. Т а к , исследования

некоторых авто­

ров

показы ваю т [2 0 ], что функция распределения

времени ремонта

для

многих видов радиоэлектронной аппаратуры близка к экспонен­

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

R(t) = 1 — e-n<“ ,

где параметр а мало отличается от единицы.

87

АНАЛИЗ ПОВЕДЕНИЯ ФУНКЦИИ ГОТОВНОСТИ

Многие системы обладают тем свойством, что большинство отка­

зов устраняется за

время,

близкое к среднему времени Тв, и лишь

в редких случаях

время

восстановления сущ ественно отличается

от Тв в ту или иную сторону. При такой ситуации распределение времени восстановления хорошо аппроксимируется частным случаем гамма-распределения. Данному закону, например, подчиняется время восстановления радиолокационной аппаратуры, аппаратуры связи, телевизоров и пр.

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

Проиллюстрируем влияние интенсивности экспоненциального вос­ становления на функцию готовности. В табл. 3 .2 приведены значе­ ния параметров функции Г ( t) при распределении времени безотказ­ ной работы системы по закону суперпозиции двух экспонент со зн а­ чением параметров = 0,02 ч -1 ; Х2 = 0,24 ч~1; сг = 0 ,4 ; с 2 = 0,6 в случае различных значений интенсивности, а следовательно, и среднего времени восстановления. Данные таблицы рассчитаны по формуле (2.24). Из таблицы видно, что уменьшение среднего вре­ мени восстановления влечет за собой уменьшение относительного отклонения функции готовности от стационарного значения и сокра­ щение длины провала, т. е. промежутка (tv t3) (см. рис. 3 .2).

 

 

 

 

 

Т а б л и ц а 3.2

Параметры

функции готовности

 

 

 

при различных значениях

интенсивности

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

 

 

 

Параметры функции

готовности

 

Г п■ 4

д

дг

Д н

/ 3. *

G, %

 

10

0,6925

0,1086

10

69

15,6

5

0,8182

0,1271

7

51

15,5

2,5

0,9000

0,1048

5

37

12,2

1,25

0,9474

0,0788

3

26

8,3

1

0,9574

0,0630

3

23

6,6

8 8

ГЛАВА 3

На рис. 3 .1 3 изображены кривые функции готовности одной

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

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

исключения участка приработки. В противном случае частые отка­ зы , связанны е с дефектами произ-

Рис. 3.13.

Функция готовности при

Рис. 3.14. Зависимость параметров

различных

интенсивностях восстано­

провала функции готовности от сред­

вления, подчиняющегося экспоненци­

него времени восстановления.

 

альному закону.

 

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

готовности в течение продолжительного времени.

параметров t3

Н а

рис. 3 .1 4 представлены графики зависимости

и Д Г

от величины среднего времени восстановления

Тв для рассмо­

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

саны в

§

3

.3 ,

при

уменьшении среднего времени восстановления

с 2 до

1

ч продолжительность провала функции готовности умень­

ш ается

с

80

до

30

ч.

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

89

АНАЛИЗ ПОВЕДЕНИЯ ФУНКЦИИ ГОТОВНОСТИ

ВЛИЯНИЕ РЕЗЕРВИРОВАНИЯ

§ 3.6

НА ПОВЕДЕНИЕ ФУНКЦИИ ГОТОВНОСТИ

 

Зависимость вероятности безотказной работы и других характеристик надежности от вида и кратности резервирования в теории надежности

изучена

довольно

обстоятельно [3 3 ].

П редставляет интерес такж е

вопрос

о влиянии

резервирования на

функцию готовн ости ..

При установивш емся режиме работы все изучение сводится согласно (1.51) к определению среднего времени безотказной работы резервированной системы и среднего времени ее восстановления. В случае неустановивш егося режима работы расчет функции готов­ ности резервированной системы можно производить по формуле (1 .35).

Влияние резервирования на поведение функции готовности

впределах участка приработки рассмотрим на примере общего

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

внастоящее время.

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

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

резервированной системы Рр (t)

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

ар (t) вы раж аю тся в общем виде

формулами

Р р ( 0 = 1 — [1 — Я ( / ) ] '" + 1

аР(0 = - Н )а (* )[1 - Я (0 Г ,

где Р (t) и

а (t) —•соответственно

вероятность безотказной работы

и плотность

вероятности

отказа

одной системы.

 

С учетом выражений

(2.12)

и (2.14)

имеем

 

 

Pp(t) =

1 _

( 1

— cje- м

— Съ<г-ыу»+1-

(3.7)

ар (t) =

(in 1) ( с ^ е - ы +

с2Х2е~х-() (1 — суе~ — сге~%**)т.

(3.8)

Считаем далее, что возможности ремонта не ограничены, ремонт вышедшей из строя системы начинается немедленно и распределение длительности восстановления каждого образца подчиняется экспо­

ненциальному

закону

с

параметром

ц. Д л я определения

функции

распределения

времени

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

резервированной

системы

Rp {t) воспользуемся

полученной в

[42,

стр. 83— 121]

формулой

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

со

1 т

(J

[1 - Я (и )]4 и | [ 1 - Я (01,

где R (t) — функция распределения времени восстановления одной подсистемы; tp — случайное время восстановления резервированной системы.

90

 

 

 

 

ГЛАВА 3

Очевидно, что

 

 

 

 

Яр( 0 = ^ К р < * } = 1 — р \* р > {\ =

 

с о

\ т

 

 

 

{tJ [1- R ( u ) ] d uJ\ [1- R ( t ) ] .

 

Т ак как по принятому выше

условию R

(t) =

1 — е~ ^ , то

из по­

следнего равенства получаем

 

 

 

 

# р (f) =

1

<.

 

(3.9)

И спользуя выражения (3.7),

(3.8) и (3.9),

по

изложенному

в § 2.5

алгоритму находим последовательности значений функции готов­

ности

резервированной

системы при

 

различных

т,

причем длительность

 

безотказной

работы основной и ре­

 

зервных

систем

распределена

по за ­

 

кону

суперпозиции двух

экспонент.

 

Н а рис. 3 .15 приведены графики

 

функции готовности описанной выше

 

модели

резервированной

системы

 

при т = 0 , 1 , 2 , из которых видно,

 

что общее нагруженное резервирова­

 

ние сущ ественно уменьш ает провал

 

функции готовности на участке при­

 

работки.

 

 

 

 

 

 

 

Аналогично

можно

 

рассчитать

Рис. 3.15. Функция готовности ре­

функцию

готовности

для

других

зервированной системы при различ­

видов

резервирования

и

стратегий

ных кратностях резервирования.

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

Следует

отметить,

 

что влияние других видов резерви­ рования на поведение функции готовности качественно не отлича­

ется от представленного на рис. 3 .15, но степень этого влияния зависит от эффективности применяемых схем резервирования.

Если время безотказной работы и время восстановления основной и резервных систем подчинены экспоненциальному закону, то функ­ ция готовности резервированной системы изменяется от 1 до kr. При этом kT всегда меньше значения функции готовности, и потому коэффициент готовности применим для оценки готовности системы.

Выражение для функции готовности в данном случае получается на основе математического аппарата теории простых марковских цепей [4 0 ]. В [22] получены формулы Г ( t) для различных видов и кратности резервирования, удобные при инженерных расчетах. Выраж ения для функции готовности при нагруженном резервирова­ нии различной кратности и различных стратегиях восстановления приведены в табл. 3.3, а для ненагруженного резервирования — в табл. 3 .4 .

91

чэ

Isa

Функция готовности системы при ненагруженном резервировании

Крат­

Коли­

 

ность

чество

Выражение дли функции готовности

резер­

ремонт­

вирова­

ных

 

ния

бригад

 

1

2

2

1

3

г m _

Р2 + М^

 

(РгеРі< —Ріер-‘)

У '

 

F2 +

+

*-2

Р1Р2 (Рі — Рз)

 

 

2р2 + 2Яр

 

х Ң р 2ер ' 1- р , е р*1)

 

2р 2

2Яр -f- Я"

РіРг

(Рі — Р2 )

г

 

, «

р 3 +

р 2Я + рЯ2

,

 

 

1

р3 +

р2Я + рЯ2 -+- X3 1

Я3 [р2Рз (Рг Ра)

ер'1РіРз (Рі — Рз) еРз/

,

 

 

Р1 Р2 (Рі Рз) еР2‘\

1

РіРгРз (Рі — Рз) (Рі — Ра)

(Рг — Рз)

 

 

р з + 3р 2Я + З р Я 2

,

 

1

6р3 + 6р2Я -Г ЗрЯ2 + Я3

1

3

Я3 [р2Рз (Рг Рз) еРіі ргр3 (рг — р3) epW —

 

 

— Р1Р2 (Рі — Рз) ерЭ<]

 

 

РіРгРз (Рі — Рз) (Рі — Рг) (Рг — Рз)

Таблица 3.3

Формулы для вычисления постоянн ых

 

Рі = — (А, +

р) — У Ц

 

 

Рг =

— (^ +

р) +

К ^ Р

 

Рі =

— -у

[(2Я + Зр) +

]/^р2 +

4Яр ]

Рг =

- у

[ ( 2 Я - Ь З р ) - 1 /'р 2 +

4Яр]

Рь Рг» Рз — корни

уравнения

р3 “h р2 (ЗА -j-

Зр,) -|- р

(ЗА,2 -J- 4uA -f-

+ 3fx2) + А3 +

(ХА2 +

Ар,2 + р 3 = 0

/?3

Ри Рг» Рз — корни

уравнения

/?2 (ЗА -)- 3(х) -f- р

(ЗА2 +

9Ар +

+

11ц2) + А3 + ЗцА2 +

6ц2А+

6ц3 = 0

ГОТОВНОСТИ ФУНКЦИИ ПОВЕДЕНИЯ АНАЛИЗ

Ф у н к ц и я

го то в н о с ти системы

при на гр уж е нно м резервировании

 

Крат­

Коли­

 

 

 

ность

чество

Выражение для функции готовности

р езер-

ремонт­

вирова*

ных

 

 

 

пня

бригад

 

 

 

 

1

Р2 + 2Ар

2 х Ч р 2ер ‘1 -

Ріер >‘)

 

(X2 + 2Ар -|- 2А2

РіРг (Р і

— Рг)

 

 

2

 

 

 

 

2

1

3

3

г m

...

(.I2 -f" 2Ар

2Х2 (p2ePl‘ ~

РіеР2/)

ѵ

р 2 + 2Ар + X2

Рір 2 (Рі р а)

 

 

р3 +

ЗА,р2 4- 6А2р

 

 

U

р3 + Зр2А + 6рА2 + 6 А,3

6А.3 [р2р 3 (Рі Ра) еР:1 — Рір 3 (Рі — р3) ep ’J —

:— РіРг ( Р і ~ Ра) ер*‘]

РіРгРз (Рг

Ра)

(Рі — Рз) (Рі — Рі )

г , , , ,

р3 +

З р Ч

+ ЗрА2

,

( }

 

(Р -Г

А)3

1

6А-3 [раРз (Ра — Рв) еРіі — РіРз (Рі — Рз) еРгі

,— РіРг (Рі — Рг) ерз<]

РіРгРз (Рг — Рз) (Рі — Рз) (Рі — Рг)

Т а б л и ц а 3.4

Формулы для вычисления постоянных

Ра = - 4

t(3 x +

+ 4 р А ]

Р » = - 4 -

[(ЗА. -ь 2 ц ) - У А * + 4 ц А ]

 

Рі = 2 (р +

А)

 

Ра — — (Р +

Ц

Рі» Рг» Рз — корни уравнения

р3 "Ь р2 (6Ä, -)-

Зр.) -f-

р

(I ІА* + 9Я(д, -f-

-J- Зр,2) -f- (I3

3|д.2Я.

-|-

6рА -f- 6Л3 = 0

Ри Pit Рз — корни уравнения

р3 + р2 (6р +

6А) +

р [ 11 (р + А)2] +

+

6 (р +

3ГЛАВА

А)3 = 0

 

Глава

М Е Т О Д Ы П О В Ы Ш Е Н И Я

4

Г О Т О В Н О С Т И С У Д О В Ы Х

 

С И С Т Е М У П Р А В Л Е Н И Я

 

ПОСТАНОВКА ЗАДАЧИ

§ 4.1

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

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

впоследующих параграфах настоящей главы .

АЛГОРИТМ ОПТИМИЗАЦИИ СТРУКТУРЫ

§ 4.2

ИЗБЫТОЧНОЙ СИСТЕМЫ ПО КРИТЕРИЮ Г(<)

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

94

Соседние файлы в папке книги из ГПНТБ