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

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

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

ГЛАВА 3

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

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

со

\ т

[1- R ( t ) ] .

 

J [1- R ( u ) ] d u \

 

{t

J

 

 

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

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 всегда меньше значения функции готовности, и потому коэффициент готовности применим для оценки готовности системы.

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

91

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

со

со

СЗ

Cf

я

ч

Ч О

*0

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

 

 

системы

 

 

готовности

• он

Ч

 

S О я X га

 

§5§2 =

 

^с> g я

о.

Функция

—■у Л

о

н £ âg к

 

« Ö Si §

=

 

сь о га р

 

* § g.= “

+

+

еч +

4 ^

+

*к а.

+

а. а.

со СО

+-I-

<<

с<

ся

CS

а. +

СЧ с^ + сч

“к +

сч

а_

ICQ„+ +1

оз «

*7;+ ° й

i" +

« + » « e +

-

d

^— %

 

Ci.«

rC1

^_|_co i +

c-e -h a. d

+ а. Q. CX

+

+

I «

101 r-

»а. а

C5. СЧ Cä.

1 °

« t i' I s s -

£ + + £2 M«-< 2<^«, P'co a. >->> 'CO

S-+5 * i<o

' Т+

5S+ Cl« Z

Д^А «-He

V f

«г +

а N

СО

+ СО ^ "5

£ о-

ечa.

со а.

+ СО

«a. +

Г5

А

со

н-

92

со

 

du

 

«3

=L

 

Cf

+

Tf

 

Я

+

 

 

<<

«г

* а

+

Е£

+

+

<з3

S

055

zL

d.

со

 

 

ч. у

+

 

 

_ о

 

 

■ос

со

СО

 

ч

 

>>

 

 

 

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

ч

еС

О ь е£

а

я =х я

НО«- t*

5

о2л -

 

<U* =

н f- 3® t

я о ® о :

322 £•:

*я с. d

СО

 

d.

соd. +

СО+=L

+ CO

+ C')

 

++

+ °

«All iS S t

ÜJ . CO

£Я +sJI

S = 4 - o + ^

I d .c o

CO

со \ Г

Q. T «

^3-

Sco .

Q .~ +

CJ 1

CL_j_ d .

CO

K +

CO

+

M

со*

o< +

CO

d.

со

+

I I

Q.

N

CL

ГЛАВА 3

+

ë +

g i o

I

& + + -- .- .

, <-£ d-

CD—■

w-1I -<£>

* 1 +

SCO

- CL

Ч -

 

 

Q.

d .

 

n

CO '

Ci.

 

Ci.

 

5

I

 

 

COd_ d. 'i

3

 

CO

 

Q.

 

Q.

 

r-(

 

CL

а©j

 

а

 

О

 

 

+

93

 

Глава

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

4

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

 

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

 

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

§ 4.1

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

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

ции. Разработка поставленных оптимальных задач

излагается

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

 

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

§ 4.2

ИЗБЫТОЧНОЙ СИСТЕМЫ

 

ПО КРИТЕРИЮ Г(<)

 

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

94

ГЛАВА 4

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

В формализованном виде задача может быть представлена сле­

дующим образом: максимизировать функционал

 

Р = шах ( f l

Ль ('«*)}>

(4.1)

......ІПп U=1

)

 

где тк — количество параллельно включенных элементов или узлов /г-го типа, образующих ступень резервирования; Рк (тк) — вероят­ ность безотказной работы k-й ступени резервирования; п — коли­ чество последовательно включенных ступеней, образующих основное соединение, при некотором ограничении, например ограничении на суммарный вес системы Wc

 

П

 

 

 

 

w c — £

wkinkSa 0,

 

(4.2)

 

*=i

 

 

 

или при совокупности

ограничении

 

 

п

 

 

 

 

Wcj — 2)

wjkinkЗг 0,

/ = 1, 2, . . . , т.

(4.3)

/г—1

 

 

 

Здесь wk —: вес элемента или узла /г-го

типа; Wcj — множество,

состоящее из т параметров, учитываемых

при оптимизации.

 

Таким образом, задача сводится к определению таких целочислен­ ных компонент /г-мерного вектора М = (т^ т 2, . . ., тп), которые максимизируют функционал (4.1) при одновременном выполнении условия (4.2) либо условия (4.3). Такого рода вопросы подробно исследованы в литературе, в частности в [1, 3, 6].

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

95

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

коэффициент готовности /гг или функция готовности в установившемся режиме Г (s). В первом случае максимизируется функционал

Т =

шах \ l f \ P k (tnk) dt],

 

(4.4)

 

"'l.............................."'n І О

/ . ' =

1

J

 

 

во втором —

 

 

 

 

 

 

 

 

 

K =

mv

шах

jП Kk(mk)];

 

 

 

 

,mn U =

1

)

 

 

Г (s) =

шах

/erft (m*) 4 - J p k (0 dt ,

(4.5)

 

mv - ’m n\k=\

 

 

1 о

J

 

где kVk — коэффициент

готовности

m

параллельно

включенных

элементов k-го типа

при ограничениях

(4.2) или

(4.3).

 

Следует отметить, что задачи оптимизации структуры по крите­ риям Р (t), Т и kr не адекватны. Прежде всего при оптимизации по Р (t) структура системы является функцией аргумента t, чего нет в случае критериев Т и kr. Кроме того, оптимальные структуры, полученные по данным критериям, различны. Проиллюстрируем сказанное примером.

Пример 4.1. Пусть имеется нерезервированная система, состоящая из трех

подсистем,

которые характеризуются

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

А. =

1 ч"1; А2

= 0,5

ч-1 ;

А3 =

0,33 ч_1 и весами (в

относительных единицах)

Ші =

1, ша =

3,

=

2. Интенсивности восстановлений каждой из подсистем равны соответственно

р і =

5 ч '1;

ра =

3,3 ч_1; р3 = 2ч '1. Требуется осуществить поэлементное нагру­

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

решения представлены в табл.

4.1.

 

 

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

 

 

Т а б л и ц а 4.1

 

 

 

полученные по различным критериям оптимизации

 

 

 

 

 

Кратность

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

подсистем

Критерий оптимизации

первой

второй

третьей

 

 

 

при

P (t ):

0,1

 

 

 

t =

2

2

2

»

t =

0,8

4

2

1

 

 

 

5

1

2

 

kr

 

1

3

1

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

96

ГЛАВА 4

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

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

Если оптимизацию производить для моментов времени, в которые

функция готовности k-то элемента системы Г* (t),

k =

1 , 2 , . . ., я,

t t (0 , оо) принимает минимальное значение, то в

силу

изученного

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

необходимо разбить

на отдельные этапы длительностью т;, 1 = 1 ,

2 , . . ., и структуру

оптимизировать поэтапно.

Теперь задачу можно сформулировать следующим образом.

Найти такие компонент

истемы я-мерных векторов

М (п = (т{‘\ т\п, • •

■,

1 = 1 , 2 ,

,

которые максимизируют

функционал

вида

 

Г(/) =

max

( П

(яі//1)! >

(4.6)

m('>..... >

 

 

 

где Г)/' (m[l)) — функция

готовности

k-и ступени

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

на /-м этапе оптимизации структуры системы, в области, заданной соотношениями

We — S тк>піПЗа 0, 1 = 1 , 2 , : . .

(4.7)

к= 1

 

Представляет интерес и другая постановка этой же задачи: минимизировать суммарный вес системы Wz на 1-м этапе ее эксплуа-

7 А. Г . Варжапетян

97

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

 

 

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

системы при

t 6 (0, оо) и всех

I будет не ниже заданного Г0, т. е. найти

W{!)=

min ( І и Ѵ

и Н 1 = 1 , 2 ..........

(4.8)

 

«}'>........m(0 l*=1

J

 

где wkm{l) — вес k-k резервированной ступени на l-ы этапе оптими­ зации структуры системы в области, определяемой соотношениями

Г0 -

П Г І 'Ч т П ^ О , 1 = 1 , 2 , . . .

(4.9)

А= 1

 

Следует отметить,

что оптимизация структуры

по критерию

Г (t, s) не имеет принципиальных особенностей по сравнению с опти­ мизацией по критерию Г (t).

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

особенно эффективен

при небольшом количестве компонентов, т. е.

в тех случаях, когда

резервирование производится на уровне круп­

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

Рассмотрим более детально максимизацию функционала (4.6), так как к этой задаче могут быть сведены многие часто встречающиеся

на практике ситуации. Обозначим функционал Г(/) в (4.6) через /дЯ. Тогда функциональное уравнение для решения данной задачи мето­ дом динамического программирования при ограничениях (4.7) будет иметь вид

fW =

max

[ftihiw — wNm W ) r U) (wNmW)},

(4.10)

0

г£

< Ш

 

 

где /дг*— доминирующая последовательность значений

функции

готовности и веса соответствующей части системы на

(N 1 )-м

шаге I-го этапа

оптимизации;

Г(/) (i%mjvZ)) — функция готовности

части системы, состоящей из

элементов N -го типа.

 

Для выбора

оптимальной

избыточной структуры в общем виде

при нескольких ограничениях могут быть использованы различные алгоритмы [3, 32]. Приближенный алгоритм [3, стр. 62— 67] позво­ ляет перейти от задачи оптимизации при нескольких ограничениях к задаче оптимизации с одним искусственным ограничением. Объем вычислительной работы в данном случае значительно сокращается по сравнению с точным алгоритмом [32]. Поэтому определение

98

ГЛАВА 4

членов оптимальной последовательности на каждом шаге оптимизации будем производить с помощью указанного приближенного алгоритма. Алгол-программа решения задачи оптимизации с использованием уравнения (4.10) приведена в приложении II. Программа составлена

с учетом того, что последовательность значений (Г(/) (шд/тдг')} определяется из выражения (1.35) с привлечением изложенного в § 2.5 алгоритма численного решения (см. приложение I). Значе­ ния функций Р (t), а (t) и г (t) вычисляются по формулам для резер­ вированных систем, например приведенным в [33].

Работу данного алгоритма иллюстрирует следующая модель.

Имеется

некоторая

система,

 

 

 

 

 

представляющая собой основное

^

 

 

 

 

соединение трех подсистем, от-

V.

._____________

2

 

казы и ремонты которых вза-

Г}„„

-------_Lr=r--------

 

имно независимы. Функции го­

 

 

 

 

 

товности подсистем имеют

вид,

0,8

 

 

 

 

изображенный

на

рис.

3.2,

и

 

 

 

 

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

рав­

0,7

 

 

 

 

ны

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

krl

=

0,90,

 

 

 

 

kr2 — 0,96,

kr3 =

0,86.

Функ­

 

 

 

 

 

ция

готовности

системы,

если

0,6

 

 

 

 

не приняты

специальные меры

 

 

 

 

 

по

повышению

надежности,

0,5

 

 

 

 

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

значение

 

 

10

го

Ьчч

(рис. 4.1,

кривая

1)

и

не

по­

Рис. 4.1. Функция готовности неизбыточ­

стоянна во времени. Для по­

ной (кривая 1) и избыточной

(кривая

2)

вышения

надежности

 

системы

 

 

систем.

 

 

имеется

возможность

 

исполь­

 

 

 

 

 

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

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

Рассмотренный выше алгоритм обеспечивает такое решение. В табл. 4.2 представлена стратегия оптимального изменения струк­ туры системы, полученная в результате просчета приведенной мо­ дели по указанному алгоритму. Десятый этап оптимизации доответ-

7*

99

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

 

 

 

 

Т а б л и ц а 4.2

Оптимальная стратегия изменения структуры системы

 

по критерию

(4.6)

 

 

 

 

Кратность

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

подсистем

Коэффициент

Номер этапа

 

 

второй

третьей

готовности

 

первой

 

1

3

3

2

0,958

2

3

2

3

0,922

3

3

2

3

0,934

4

3

2

3

0,944

5

3

2

3

0,947

6

3

2

3

0,951

7

3

2

3

0,953

8

3

2

3

0,955

9

3

2

3

0,957

10

2

2

4

0,958

ствует установившемуся состоянию системы. Это означает, что в даль­ нейшем структура не меняется и значение М = (2, 2, 4) является оптимальным для данного периода.

Алгоритм дает возможность решать также задачу оптимизации структуры системы с использованием уравнения (4.8) при ограни­ чениях (4.9). Для иллюстрации этой задачи в условиях рассмотрен­ ного примера выбираем такую стратегию резервирования в течение

указанного

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

и при этом

готовность ее была не ниже Г (t) — 0,91. Оптимальное

решение представлено в табл. 4.3, из которой следует, что минималь­

ный вес системы

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

периоду эксплуа-

 

 

 

 

 

 

Т а б л и ц а 4.3

Оптимальная стратегия изменения структуры системы

 

 

по критерию

(4.8)

 

 

 

 

 

Номер

Кратность

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

подсистем

 

Коэффициент

 

 

 

Суммарный

этапа

первой

второй

вес

готовности

 

третьей

 

 

1

2

 

3

2

7

0,930

2

3

 

2

3

8

0,922

3

3

 

2

3

8

0,934

4

2

 

2

3

7

0,935

5 .

2

 

2

3

7

0,925

6

2

 

2

3

7

0,928

7

2

 

2

3

7

0,932

8

2

 

2

3

7

0,934

9

2

- 2

3

7

0,936

10

2

 

1

3

6

0,938

100

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