книги из ГПНТБ / Варжапетян, А. Г. Готовность судовых систем управления
.pdfГЛАВА 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. |
|
со |
|
|
|
ч. у |
+ |
|
|
_ о |
|
|
|
■ос |
со |
СО |
|
ч |
|
||
>> |
|
|
|
Ф ункция готовности системы при нагруженном резервировании
ч
еС
• |
О ь е£ |
а |
я =х я |
S« |
НО«- 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) = |
шах |
IП |
/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