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

книги из ГПНТБ / Варжапетян, А. Г. Готовность судовых систем управления

.pdf
Скачиваний:
3
Добавлен:
22.10.2023
Размер:
12.26 Mб
Скачать

ГЛАВА 4

тации, где М = (2, 1, 3). При такой структуре значение функции готовности системы в момент, соответствующий, например, началу 4-го этапа, было бы равно 0,785, что существенно отличается от до­ пустимого Гдоп. Функция готовности полученной системы изобра­ жена на рис. 4.1 (кривая 2).

ОПРЕДЕЛЕНИЕ ЗАВИСИМОСТЕЙ

§ 4.3

МЕЖДУ ХАРАКТЕРИСТИКАМИ ГОТОВНОСТИ

 

ИСТОИМОСТИ

Внастоящее время успешно развивается экономический подход

коценке готовности технических устройств. Повышение готовности

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

Функция

стоимости обладает следующими

свойствами:

1 )

с (kr0,

kr)

0 , так как стоимость всегда

положительна;

2 )

с (£г0,

1 ) =

о о ,

т. е. разработка аппаратуры с коэффициентом

готовности,

равным

единице, требует бесконечно больших затрат;

3)с (/гг0, kr0) = с0, где с0 — стоимость аппаратуры с коэффи­

циентом готовности kro;

4)с (£г0, kr) не убывает по kT при фиксированном kr0 и не возра­

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

Указанными свойствами могут обладать многие функции, -в част­ ности следующая зависимость [23, стр. 43; 25]:

-

<4 - " >

где kn0 — коэффициент простоя системы с первоначальным значе­ нием коэффициента готовности; kn — коэффициент простоя, соответ­ ствующий системе с повышенным значением коэффициента готов­ ности.

Параметром f можно характеризовать эффективность вложения средств для повышения готовности. Чем меньше /, тем меньше не­ обходимо затратить средств для достижения заданного уровня готов­ ности. Рассмотренная модель (4.11) выгодно отличается от других моделей тем, что она учитывает изменение стоимости системы в 'за­ висимости от изменения ее готовности, обусловленного изменением показателей безотказности и-восстанавливаемостіт-системы.

МЕТОДЫ ПОВЫШЕНИЯ готовности

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

* г = 1 — О — U " 1.

(4.12)

Стоимость резервированной системы в случае нагруженного ре­

зерва сн. р будет определяться

из выражения

 

сп. Р=

тс0.

 

(4.13)

После подстановки (4.13) в (4.12)

и преобразования получим

 

^н. р — С0

In (1 —

k r )

(4.14)

In (1 —

k TB)

 

Допустим теперь, что предложен метод, при котором увеличению коэффициента готовности до значения krl соответствует изменение

стоимости до значения сь

причем

сн р, т.

е.

 

Асн. р< А с 1,

(4.15)

где Дсн. р •— приращение

стоимости

системы

в результате приме­

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

 

/ dci \

^

(

dcK. р \

(4.16)

 

\ d k г / йг> і гоч

\

dkг ) k T> k rB

 

 

Используя выражение (4.14), получим следующий критерий

целесообразности

применения

метода:

 

( — )

 

< -

_______ £о______

(4.17)

 

(1 —

k r) ln (1 ^го).

\ d k г Д г >/,Ѵо

 

 

Рассмотрим функцию (4.11). Пусть с2 — стоимость системы с ко­

эффициентом простоя

/гп.

Тогда

 

 

 

 

С2

_ (

^ПО

__ (

1 0

^ I

(4.18)

Со

V

/

\

1— kr

J

 

Функция стоимости, определяемая выражением (4.18), и функция стоимости системы при нагруженном резервировании имеют общую исходную точку с координатами kr0 и с„. Чтобы найти другие общие

точңи, решим

уравнение

с2

= сн. р,

или

 

( 1

- W ln ( 1

-

kJ = ( 1 -

ЬаУ ln ( 1 - *rt),

(4,19)

102

' ГЛАВА 4

Графическое решение этого уравнения представлено на рис. 4.2, где построены вспомогательные функции

Уі = {1 — £ r)f l n ( l — k ry,

г/ 2 = (1 — Ar0 )M n (l— Ar0).

Функция г/i при kr = k'r имеет минимум. Из рис. 4.2 видно, что кривая у г и прямая у 2 имеют две точки пересечения. Значит, урав­ нение (4.19) имеет два решения krl и kr2.

Найдем минимум функции у х, приравняв нулю ее производную:

или после подстановки у х

X ln (1 — kT) = 0.

После упрощения получим

1

,

/ 1п (1 — *г) _

л

 

 

1 kr

\ ~ k r

Рис. 4.2. Графическое решение уравнения

откуда

 

 

(1 - kr)f ln

(1 -

kT) = (1 - kro)f X

 

£ =

1 - е ~ f .

X

ln

(1 ^то)*

 

(4.20)

 

 

Определим характер пересечения функций с2 и сн%р в точке Аг0,

для чего найдем их производные

de2

^гоУ'

f

 

dkr

 

(1 —

 

dcH, p

___________ £o_____ _

1

dkr

ln (1 ^ro)

 

^ — kp

В точке kro производные соответственно равны

 

I

dc2 N

---с

f

.

 

\

dkr Л г= * го -

1-Аро *

(

dcH, р \

 

_______________ £о________

\

dkr /Аг = * г0

0

—Аго) ln ( 1 kro)

Проверим выполнение условия (4.16):

/

dc2 \

f

dcn. р

\

\

dkr / А г = * г „

V

dkr

) kr = k ro'

После подстановки соответствующих выражений получим нера­ венство

f <

ln ( 1

kro)

 

103

МЕТОДЫ ПОВЫШЕНИЯ готовности

 

которому может удовлетворять решение

 

1

 

ér0 < l — Г Т -

(4.21)

Сопоставляя (4.20) и (4.21), приходим к выводу, что для значе­

ний Г о < 1 е 1 общей исходной точкой функций с % и сн. р

с координатами Го и с0 является точка Г і на рис. 4.2. В этом случае условие (4.21) выполняется, а для

значений

Го > 1

е f исход­

ной

точкой

будет

служить точ­

ка

Г 2, и тогда

 

 

( 'Г \

\

Л 4сц. р \

 

 

\ d k r ) kr= k r

 

\ dkr

Уаг= а-г0 '

 

 

На основании этого можно вы­

 

явить характер

взаимного распо­

 

ложения графиков

функций

с2

 

и сн р. Кривые с2 и сн. р на рис. 4.3

Рис. 4.3. Взаимное расположение

пересекаются

в точках I и 2.

Ис­

ходной точкой с координатами Г о

функций С% и сн. р.

и с0 может быть любая из этих

точек в зависимости от величины коэффициента готовности исход­

ного

образца. Точка

2 будет

исходной в

случае,

если

 

 

 

1

 

 

 

 

Г о > 1 —С0 f

 

(4.22)

При

/ = 1 имеем Го

> 0 ,6 5 7 .

Соотношение

(4.22)

на практике

всегда выполняется, и, следовательно, исходной точкой всегда будет точка 2.

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

На рис. 4.4 представлена область расположения графиков всех возможных функций стоимости, имеющих исходную точку с коорди­ натами Г о и со- Эта область ограничена линиями 2 и 3. Кривой /,

соответствующей случаю нагруженного резервирования, она разде­ ляется на две области / и II. Допустимой является область I. Гра­ фиками допустимых функций стоимости являются кривые, имеющие общую точку с координатами Г о и со. располагающиеся ниже кри­ вой 1 и не пересекающие ее во всем диапазоне значений Г > Го-

Такие графики имеют логарифмические функции вида c ' = c 0log/(l — Г)-

104

ГЛАВА 4

На рис. 4.5 изображены кривая, соответствующая случаю на­ груженного резервирования ctl p, логарифмическая кривая с' и гра­ фик искомой функции с. Кривая с может быть получена путем пере­

носа вверх

кривой с

на

величину со— с'о,

т. е.

 

 

с =

с0

log/- ( 1

К) -|- (с0 +

Со),

 

откуда

 

 

 

 

 

 

 

c = ' » ( 1 +

l0 & - r 5 f e - ) - ■

Ѵ -Щ

Основание

логарифма

/

может

изменяться

в пределах

О < / <

< 1 — kr0. Чем меньше /, тем ниже пойдет кривая и тем совершеннее

Рис. 4.4. Область расположения гра­ фиков всевозможных функций стои­ мости повышения готовности.

Рис. 4.5. Графическое построение ис­ комой функции с с использованием ло­ гарифмической кривой с' и кривой нагруженного резервирования сн. р.

метод повышения готовности. При / = 1 — krQ функция (4.23) превращается в функцию, характеризующую нагруженное резерви­ рование, а при / = 0 имеем идеальный случай увеличения надежности без затрат. Если известны исходная стоимость с0 и коэффициент

готовности, а также стоимость с и соответствующий ей коэффициент kr улучшенного образца, то согласно [22, стр. 43] параметр / можно определить по формуле

< 4 - 2 4 >

Как показано в [23], в результате анализа и обработки стати­ стических данных можно утверждать, что параметр / изменяется в пределах

l ~ k rQ< f < ( l - k r0)\

105

МЕТОДЫ ПОВЫШЕНИЯ готовности

Зависимость, аналогичная (4.23), справедлива и для вероятности безотказной работы:

С = Со( 1 + 1 о &£ £ ) ,

где Р и Р й — вероятности безотказной работы улучшенного и исход­ ного образцов соответственно.

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

стоимости

 

 

 

c =

Co( l + l o g f A ) .

(4.25)

Основание логарифма /

имеет пределы О <С / <С 1 — Р 0

или для

высоконадежных систем

0

< / < Ы.

 

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

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

примет

вид

 

С

 

 

 

 

 

 

kr = 1 — (1~ у

f С°

(4.26)

При

/ = 1 — kr0 выражение

(4.26)

превращается в

выражение

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

из т =

элементов при неограниченном

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

Ч>

 

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

повышения готовности,

можно записать

 

kr — f {X, [д,,

т, способ резервирования).

(4.27)

Функция (4.27) позволяет рассмотреть влияние на коэффициент готовности как характеристик безотказности, так и характеристик ремонтопригодности.

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

106

ГЛАВА 4

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

Для резервированных систем при нагруженном резервировании возможны два режима восстановления:

1 ) с восстановлением отказавших блоков после отказа всей

резервированной системы; 2 ) с восстановлением отказавшего блока сразу же после его

отказа.

Для первого случая среднее время безотказной работы опреде­

ляется соотношением

 

 

 

Т ср і и п г

__І _

 

(4.28)

m -j- i

І=п

 

где Tcp — среднее время безотказной

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

п — число блоков в резервированной системе; m — число исправных блоков, при котором резервированная система сохраняет свою рабо­ тоспособность.

Для второго случая

 

т \n—m+ 1

mCl

(4.29)

~Т„

где Т в — время восстановления

отказавшей системы.

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

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

 

 

 

Т а б л и ц а 4.4

Характеристики надежности и готовности резервированных систем

 

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

п * 4

Мажоритарное резервирование

45

0,860

Дублирование

80

0,920

»

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

200

0,965

107

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

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

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

АЛГОРИТМ ОЦЕНКИ ЭФФЕКТИВНОСТИ

§ 4.4

РАЗЛИЧНЫХ СПОСОБОВ ПОВЫШЕНИЯ ГОТОВНОСТИ

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

1 ) алгоритма нахождения оптимальных структур для различных

способов повышения готовности; 2 ) алгоритма мажорирования полученных оптимальных структур

с целью получения результирующей оптимальной структуры.

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

Сформулируем поставленную задачу о получении минимального коэффициента простоя при заданной стоимости с помощью методов теории динамического программирования. Для вывода основного функционального уравнения запишем рассматриваемую задачу в виде

min kn{Xj,

х2........ х„) =

 

П

 

 

0 « 2 ctxі<с0

 

 

£=1

1— П kni (*|)

 

п

(4.30)'

= min

 

ОС 2

i = l

108

ГЛАВА 4

где kn — коэффициент простоя системы; /?П(- (xt) — коэффициент простоя і-й подсистемы, когда в ней имеется х с резервных элементов; Cj- — стоимость одного элемента і-го типа; с 0 — величина, ограничи­

вающая

стоимость

системы.

 

 

 

 

Заметим,

что

 

 

 

 

 

 

 

 

min

k

n

( X ) —

 

min

П knn (хп) X

 

П

сіхі^со

 

 

 

спх „ < с „

ЛЛ

 

 

°< 2

 

 

 

 

 

 

 

і = 1

 

 

 

 

 

 

 

 

 

 

 

(1 - 1

 

 

 

 

 

 

 

х П

 

1 —

min

П k ni (Хі)

(4.31)

 

 

 

( = 1

л- i

д..

 

 

 

 

 

 

 

0 < 2

сіЧ^со~сп

 

 

 

 

 

 

 

1 = 1

 

 

 

где X =

(хІУ

х.2,

. . .,

хп).

 

 

 

 

Введем в

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

некоторую функцию

(с*), которая

определяется в нашем случае следующим образом:

f k ( с * ) = "

min

k n ( X ) =

min

1 - П

1

П k ni (Xi)

 

 

П

 

 

П

 

 

 

 

 

°=s 2 cixi<’c"

°< 2

cixi<c'

 

 

 

 

i=i

 

 

i=i

 

 

 

 

T. e. индекс при функции f означает размерность минимизируемой

функции /гп.

Тогда выражение (4.31) можно записать в виде

fn (с0) =

min

( П

knn (xn) -I-

(с0 спхп) —

 

° < Ѵ ( і < со 1 х п

 

 

 

 

П

/еп„ (х„)

(с0 — с(1лд|.

(4.32)

 

-'Vi

 

 

 

J

 

Выражение (4.32)

и

есть

функциональное уравнение,

дающее

рекуррентное решение поставленной задачи. Идея решения функцио­ нального уравнения состоит в следующем:

а) определяются оптимальные двумерные векторы состава системы для первого и второго элементов при всех значениях показателя стоимости, не превосходящих с0;

б) находятся оптимальные трехмерные векторы состава системы для третьего элемента и соответствующих пар (xlt х 2) при всех зна­ чениях показателя стоимости, не превосходящих с0;

в) процедура продолжается до тех пор, пока не будет найден оптимальный п-мерный вектор состава для /г-го элемента и соответ­ ствующего оптимального вектора (xlt х 2, . . ., хп) при значении пока­ зателя стоимости, равном с0;

г) выделяется оптимальное хп и соответствующий оптимальный вектор (хъ х 2, . . ., хп_х), которые в совокупности и дают оптималь­ ное решение.

109

МЕТОДЫ ПОВЫШЕНИЯ готовности

Алгоритм решения функционального уравнения достаточно раз­ работан и подробно описан [1, 32]. С помощью этого алгоритма может быть построена оптимальная синтезируемая структура при использовании какого-либо одного метода повышения готовности.

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

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

следовательности {(&п0, с0); (&п1, Сі); (Ігп2, с2)}, для

которой переход

в состояние с меньшим коэффициентом простоя /еп

происходит с ми­

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

Рассмотрим работу алгоритма «мажорирование» на следующем примере. Пусть даны две исходные оптимальные последовательности, соответствующие двум способам повышения готовности (табл. 4.5). Необходимо с помощью алгоритма «мажорирование» получить ре­ зультирующую мажорирующую последовательность, учитывающую оптимальное сочетание рассматриваемых способов повышения готов­ ности.

Выбираем среди всех членов заданных последовательностей член с минимальной стоимостью. Им является первый член первой после-

г 0,0571

0

 

 

довательности (

1 1 5 /• Среди оставшихся членов исходных последо­

вательностей вновь выбираем член с

минимальной стоимостью.

Это второй член первой последовательности |

1 3 g \ Он имеет коэф­

фициент простоя

kn =

0,045, меньший,

чем .у

предыдущего члена

ПО

 

 

 

 

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