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

книги из ГПНТБ / Системы автоматического и директорного управления самолетом

..pdf
Скачиваний:
49
Добавлен:
20.10.2023
Размер:
8.62 Mб
Скачать

Тогда для п элементов при условии, что первым отказывает эле­ мент 1, будем иметь

P [ ( x ) P ^ n(x)dx,

(3.50)

о

 

где Р2 записывается аналогично (3. 19):

 

Рг-п(х )= П P i№-

(3-51)

г =2

 

Для случая, когда имеются четыре элемента 1, 2, 3 и 4 с плот­ ностями распределения времени безотказной работы соответ-

dm,

Обязательный

. Работа элемента I

 

 

Работа элемента Z

1

АЛ

Т 2

\/ d v z

/ /

Работа элемента /

 

 

Работа элемента 2 --------------------------

ХД><)

--------------

Необязательный отказ

Рис. 3. 11. К реализации отказа системы из четы­ рех элементов

ственно /2(0 , Ы 0 и Ы 0 >определим вероятность того, что за время работы t элемент 1 откажет первым, элемент 3 откажет после элемента 1. Состояние элемента 2 не играет роли после отказа элемента /, а состояние элемента 4 не играет роли после

отказа элемента 3.

 

1,

бу­

Приняв во внимание, что первым отказывает элемент

дем иметь

 

 

 

Q2=

— j P;(Ti)P2_4(ti)rfti.

(3.52)

 

о

 

 

Для вычисления вероятности Я2-Д ti) вновь применим

фор­

мулу (3.50). Из рис. 3.

11 легко видеть, что

 

 

P4-iixi) =

(3-53)

 

ч

 

 

120

Подставив (3.53) в (3. 52), получим

 

/

 

t

 

 

 

=

J

fTi)/^2(Ti) f P^(x2)P4(x2)dx1dx2

(3.

54)

или

0

 

 

 

 

 

 

 

 

 

 

 

t

t

 

 

 

(3.

55)

Q2= j

j

^('Г 121)Р ;(т242)а,т1 а,т2.

 

*1

 

 

 

 

 

Пусть теперь имеется 2 k

элементов; 1, 2, 3 ...

2 k. Требуется

определить вероятность того,

что за время работы

t

последова­

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

Искомую вероятность можно определить с помощью много­

кратного

использования

формулы

(3.50). Окончательно

будем

иметь

 

 

 

 

 

=

j

Р ' М Р М j Рз ( Ъ) Р*( * а) j Р ь ( * з ) Х .

 

 

0

 

тд

i2

 

 

Х р л хз)

f .....

f РшХ хк ) Р ^ к ^ х ^ х 2..Мхк.

(3.56)

Выражение (3.56) справедливо для любого закона распределе­ ния времени исправной работы элементов. Для наиболее часто используемого экспоненциального закона надежности формула (3. 56) запишется в следующем виде:

=

.....V

-1

j е-<х‘+х*>т* J “( ^3+^4) т2X

 

X Г

... f

е

(X2ft-i+x2ft>tft^r1fl'T2...r/tft.

 

(3.57)

i

 

4 - i

 

 

 

Для удобства нахождения решения применим к формуле пре­

образование Лапласа [12]:

 

 

 

<?*(«)=

____________________________

(3. 58)

s (s + Xyo)(s + Xjh)—4s +

1)1

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

П ^2:-l

 

 

или

 

Qk(s) = -----т д --------------

 

 

 

 

 

 

s Ills - Xjn)

 

 

 

 

 

 

1-0

 

 

121

(3.61)

где s — оператор,

2ft

 

'Vi ■ v v

(3.59)

/ = 2 / + 1

 

При неравных корнях знаменателя B(s) выражения (3.58) обратное преобразование Лапласа для Qh(s) будет иметь вид

к

л2<-1 ^

sil

 

П

B. {Sl)

(3. 60)

i - 1

___о

(Si)

 

1=0

 

 

где si = —‘ki, т. е. /-й корень знаменателя функции (3.58); B'(Si) — производная знаменателя в точке s,-.

Принимая во внимание, что sy— 0, окончательно имеем

 

П ^2i—1

1

 

a кУ0/

Q k ( t )

h-1

 

ft- 1

 

i-1

 

П Xy;

 

( — ^уо) П (Ху,- — Xy0)

 

-

1=0

 

i=i

 

 

 

 

ft-1

 

 

№~ху1*

+

\ у

j

1-1

 

ft-1

 

1=1

( — Х у / ) П ( Х у , - —

П (Ху,-— ^yi)

 

 

 

i'=0

 

1=1+1

-x■ai-D* ft—2

( ~ Ay ( k - l ) ) П (Xji;— X

1=0

Из этого общего выражения нетрудно получить частные слу­ чаи, соответствующие различным к. Так, при &=1

*•(/0 —

'-yi

 

и

Qi

е

Уй !

(3.62)

при k = 2 будем иметь

 

 

куо

 

 

 

 

 

 

 

 

— \

+ *2+ '-3 +

^41

Nl — ^'3 "I-

\ 1/2

 

Q2

е

1т‘

 

 

 

~ХУ1*

(—Ху0) (XV1 — Х(/0)

( — ^-</l) (Хуо —

О-УО^У!

Аналогично при k = 3 получим

 

 

 

 

 

(3. 63)

 

 

 

 

 

 

 

^уо—''•1 + h + ^3+

^-4 + Къ+ V,;

 

^ Д ^ 'З + ^ + ^ + ^в!

ХУ2= Х5 +

Хв;

Kysz

 

 

1

 

 

.

\ t

 

 

Qg-- 3 5

 

 

е ~

хуо‘

 

 

^У0^У1^У2

( — ^УО) Н у 1 —

Хуо)

( X.t/2 —

Хуо)

V

 

 

 

 

У2

(3. 64)

+ ■

 

 

 

 

 

 

(~^У l) (Ьуо byl) Ну2 Хуо

(— Ху2) (Хуо — Ху2) (Ху! — Ху2)

122

Очевидно, выражение (3.56) справедливо для любой реали­ зации отказа системы, а не только для обобщенной. В качестве примера в табл. 3. 2 приведены формулы, описывающие вероят­ ности <7ь </2, • ■•, Яи конкретных реализаций отказа системы, представленной на рис. 3. 8, и вероятности Q\—Q4 обобщенных реализаций. Воспользовавшись формулами (3. 62) — (3.64), не­ трудно убедиться, что выполняются следующие равенства:

Qi — <7i + <72 + '?3 + ‘74 + <75; Q‘i = Яъ~\~Ят^г Яв'

Q3r=<79_b 9 io _r<7ii~b*7i2~b<7i3>

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

Qc = Ql + <22+ Q3+ Qi'

Полученные интегральные формулы справедливы и для слу­ чая, когда в системе с резервированием замещением условия работы резервных элементов до их включения в контур управ­ ления и после включения различны. Математически это выра­ жается в различии характеристик надежности таких элементов до и после момента их включения в контур управления. Так, для системы, изображенной на рис. 3. 8, где элемент D до его вклю­ чения в контур управления работает в облегченном режиме и имеет вероятность безотказной работы PDo(0 , а после момента включения — PDn{t), получим (см. табл. 3.2):

 

(3. 65)

^Oo(T'i)j" Рд (Тг^ Р.л (^2)Рс ("П) dx^dx^.

( 3. 66)

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

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

V= 1

123

Реализа­

реали­ ция зации отказа

1А

2СА

3

DA

 

4

CDA

 

5

DCA

 

1

1

 

6

C B

q%=

 

 

7

CD B

q7 =

8

D CB

q%=

9

BD

<D =

 

10

B C D

<7io =

11

BAD

q n =

12

B C A D

qn =

13

B A C D

q 13 =

14

DB

q 14 =

Вероятность реализации отказа

t

Ч \ ~

 

J

Р 'а ( t i )

p B (Тц) P c

( Ti)

p D (Ti) d x 1

 

 

 

0

 

 

 

 

 

 

 

 

Чч =•• J P C (Tl)

f

Р 'л

(T2) P В ( t 2) P D ( t 2) d x l d X 2

 

0

 

 

Xj

 

 

 

 

 

 

 

t

 

 

 

t

 

 

 

 

 

 

J

P D O i )

j

P

'a

( t 2> p B ( t 2) P c (T2) d x x d x 2

 

0

 

 

Xj

 

 

 

 

 

 

 

 

t

t

 

 

 

 

 

 

 

?4 =

j

P C (Xl) J

P D ( T2) j

P 'a

з ) P b

(T3) d x l d x 2 rfT3

 

 

 

0

 

 

 

X!

xa

 

 

 

 

 

*

 

 

 

t

<

 

 

 

95

=

J' f f l

(Til) J P c (T2) J

Р д

( t 3) Р л

(T3) d x x d x 2 d x 3

 

 

 

0

 

 

 

Tt

Tj

 

 

 

tt

/

P'B ( x 2)

P A (t 2) P d ( t 2) dxx dx2

j P'c ( t ,) j

0

X ,

 

 

 

 

 

t

 

t

f

 

t

 

J

P c ( x x) J

P D ( t 2)

[ P B ( т 3) P A ( t 3) dxx dx2 dx3

0

 

xL

 

 

x2

t

 

t

 

 

t

 

— !' P d(T i ) j

P c (T2)

j

р'в (X3) P a (T3) dxx dx2 dx3

0

 

xi

 

 

x2

t

t

 

 

 

 

 

.1 p b

( Ti) ,f P'd ( t 2)

p A (T2) p c (T2) dxx dx2 ■

0

'ti

i

 

 

t

 

t

 

 

 

 

j

P'B ( t i ) j

P'c ( t 2)

j

P ’d ( t 3) P a ( t 3) dxx dx2 dx3

0

 

Xj

 

 

x2

к

 

f

 

 

1

P D ( т 3) P c (x3)d x xdx2dx3

— \ P ’B ( T i ) j P ^

( t 2) f

0

 

 

 

 

T2

,

К

P'c (T2 )

l.

,

f

\ Р'в ( t i ) j

.1 P a

( т з) j p 'd (T4 ) dxx dx2 dx3rft4

0

X ,

'

 

X 2

 

X 3

t_

t

 

 

t

 

t

]' Р'в ( t i ) J' P'a (t 2) f P c

(т з) j' p 'd (T4 ) dxx dx2 dx3rft4

0

T ,

 

 

T ,

 

T 3

t

t

 

 

 

 

 

j' p 'd (T i ) j

P B ( т 2)

P A ( x 2) P c (t 2) dxx dx2

0

*4

 

 

 

 

Т а б л и ц а 3 . 2

Вероятность обобщенной реализации отказа

t

Q i = — j P 'a ( t i ) р в (T i ) d x x 0

/

t

Q2 = j

P c (xi) j P в (т2) P a ( t 2) dxxdx2

0

X,

t

Qs = f p's ( T i ) p A (Ti ) X

b

t

X p ( T , ) § P D’ (x2)d xxdx2

О

1!

£

и сохранении нескольких первых членов ряда. При k = \, ограни­ чившись двумя первыми членами ряда Маклорена (v = l) вме­ сто (3.62), будем иметь

 

 

 

=

 

 

 

(3.68)

При k = 2, ограничившись тремя*

первыми членами ряда

(v = 2),

вместо (3.63) получим

 

 

 

 

 

 

 

^ = y V

/ .

 

(3.69)

При k — З, сохранив четыре первых

члена ряда

(v = 3),

вместо

(3. 64)

будем иметь

 

 

 

 

 

 

 

Qs= у

 

W

3-

 

(3.70)

И,

наконец,

применив метод

 

математической индукции, при

произвольном

k, ограничившись

 

(х + 1)-ми**

членами

(v = £)

разложения (3.67), вместо (3.61) получим

 

 

 

 

*

 

*

 

 

 

 

 

 

П ''2i~!'

 

(3-71)

 

 

 

/-1

 

 

 

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

Относительную погрешность формулы (3.71) можно опреде­ лить, если точное решение выражения (3.61) представить в сле­ дующем виде:

Q n = y i Aj t - xi ^ y i AJyj ,

U ) U )

где Aj — коэффициенты, являющиеся функциями от интенсивно­

 

стей отказав элементов;

 

 

 

 

Как известно,

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

погрешность А

величины

у

А (у) определяется

по формуле

 

 

 

 

 

 

 

юоНо,

 

 

 

 

У

 

 

 

где

у — точное, а

у — приближенное

значение

величины

у.

Нетрудно показать, что

 

 

 

 

 

 

A(Ajyj)= A(yj\

 

 

 

*

При v = 1 будем иметь Q2 в 0.

 

иметь Qh= 0.

 

 

**

При числе членов, меньшем v -И ,

будем

 

 

126

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

Из сказанного следует, что

l(Q k)= max[l{y1), Д (г/2) ,- , Д (г/,-),...].

Таким образом, погрешность при приближенном определении вероятности отказа системы не превосходит погрешности от за­ мены экспоненциальной функции с наибольшим показателем несколькими первыми членами ряда (3.67). При v = l эта по­ грешность не превосходит значений, приведенных в табл. 3. 3. При v ^ 2 эти погрешности уменьшаются.

Таблица 3. 3

л-, не более

0,045

0,135

0,19

 

Д (е х), не более

0,1%

1%

а

 

X

4,5 -1 0 -з 1/ч 1,35-10-2 1 / ч

1,9-10-2

1/ ч

— , не более

t

 

 

 

 

В этой же таблице даны предельные

значения

интенсивно­

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

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

Минимальный, самый низший порядок отказа системы, пред­ ставляет собой отказность системы. Применительно к системе на рис. 3. 8 отказ «переключение управления на исправный ре­ зервный подканал с последующим отказом последнего» имеет второй порядок по числу определяющих элементов В и D (см. табл. 3. 1) и состоит из пяти реализаций отказа: 9-я — второго порядка; 10-я и 11-я — третьего, а 12 и 13-я — четвертого по­

127

рядка. Однако рассматриваемая система является одноотказной, так как неконтролируемый отказ ее основного подканала (отказ элемента А) имеет первый порядок.

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

Q c^V Q * .

(3.72)

( f t )

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

Q*+i«Q *-

(3-73)

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

Qc~Q*mi . = Qr-

(3.74)

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

системы, изображенной

на рис. 3.8, можно приближенно определить с помощью только

первой обобщенной реализации отказа (см. табл.

3.2)

и поло­

жить ее равной

 

 

 

Q c~ - f P A{*l)PB^l)dXV

 

 

о

 

 

 

Приближенными формулами

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

пользоваться

при вычислении различного рода

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

системы

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

P '= T lP i( t y ,

(3-75)

/= 1

 

 

Q c = ? i( 0 + 2 ? '( * )

j=1

(3.76)

1=2

 

где п — число элементов соединения.

 

 

При малых величинах qu пренебрегая их произведениями,

нетрудно получить весьма простое выражение:

 

П

 

 

=

 

(3.77)

/=1

 

 

128

Отказность резервированной системы определяется протя­ женностью критического участка, характеристиками надежности элементов критической части САУ и требуемым уровнем безо­ пасности полета. Теоретически возможно, что при высокой на­ дежности элементов САУ и сравнительно малой протяженности критического участка требуемый уровень безопасности будет достигнут и без применения резервирования элементов, т. е. 'САУ может быть одноотказной. При существующих же характе­ ристиках надежности исходных элементов системы для дости­ жения требуемого уровня безопасности даже при сравнительно малой протяженности критического участка необходима как минимум двухотказная система. Это подтверждается данными табл. 3.4 (см. результаты под чертой), в которой приведено

 

 

 

 

Таблица 3.4

 

 

Минимально-допустимое 7 сР (в часах) при

 

 

^кр

г= 1

г —2

г = 3

г = 4

 

 

10

ч

100 000 000

31600

2150

562

5

ч

50 000 000

15 800

1075

281

1

ч

10 000 000

3160

215

56,2

30

мин

5 000 000

ф580

108

28,1

3

мин

500 000

158

10,8

2,81

30

с

83300

26,4

1,8

0,47

18

с

50 000

15,8

1,08

0,28

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

(3.78)

причем (Эдоп было принято равным Ю-7. Хотя результаты табл. 3.4 носят ориентировочный характер, они позволяют утверждать, что при увеличении протяженности критического участка при прочих равных условиях для обеспечения безопасности полета требуется большая отказность системы. Например, если крити­ ческий участок по времени равен продолжительности всего по­

5

132

129

лета, то требуемый уровень безопасности может быть достигнут применением как минимум четырехотказной САУ (табл. 3.4).

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

Контроль готовности САУ, его значение и виды

Для выяснения возможности снижения отказности САУ про­ анализируем приближенное (выражение (3.47), связывающее уровень риска с вероятностью отказа критической части САУ. Из этого выражения следует, что величины составляющих, опре­ деляющих значение уровня риска управления самолетом посред­ ством конкретной САУ, зависят от момента времени tQ. Первое, основное слагаемое, определяется при у с л о в и и исправности критической части САУ непосредственно перед выходом на кри­ тический участок полета. Исправность или неисправность кри­ тической части САУ контролируется с помощью специальных средств. Целью контроля является определение соответствую­ щей готовности критической части САУ к работе на критическом участке полета. Такой контроль называется контролем готов­ ности САУ к работе на критическом участке трассы, или, сокра­ щенно, контролем готовности САУ.

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

Положительный результат является основанием для исполь­ зования основного варианта полета, отрицательный — запасного.

В наиболее простом случае сам факт автоматического выдер­

живания

заданной траектории

полета с требуемой точностью

в момент

контроля готовности

tK,v может свидетельствовать об

исправности САУ при условии,

что именно проверяемая часть

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

130

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