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

Учебное пособие 999

.pdf
Скачиваний:
4
Добавлен:
30.04.2022
Размер:
707.4 Кб
Скачать

холодным резервом. При отказе машины М1 включается машина М2 , а при отказе М2 включается М3 . До включения

каждая из резервных машин находится в холодном резерве и отказать не может. Пустьλ1 — интенсивность потока отказов

основной машины, λ2 — остальных машин. Представим процесс, как Марковский случайный процесс (рис. 12).

Рис. 12

Введем обозначения: S1 — работает основная машина;

S2 — работает первая из резерва машина; S3 — работает вторая из резерва машина; S4 — не работает ни одна машина.

Система уравнений Колмогорова для вероятностей состояний будет

 

 

 

 

 

dp1

= −λ p ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

1

1

 

 

 

dp2

 

 

 

 

 

 

 

= −λ

2

p

2

+ λ p ;

 

 

 

 

 

 

dt

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

dp3

 

= −λ2 p3 + λ2 p2 ;

(20)

dt

 

 

 

 

 

 

 

 

 

 

dp4 = λ2 p3 . dt

Добавим к системе нормировочное условие p1 + p2 + p3 + p4 =1.

Интегрируя первое уравнение при P1 (0) =1, получим

p1 (t) = eλ1t .

Подставляя во второе уравнение, будем иметь dpdt2 = −λ2 p2 + λ1eλ1t .

40

Интегрируя с начальным условием P2 (0) = 0, получим

p2 (t) = λ2 λ1 λ1 eλ1t λ2 λ1 λ1 eλ2t .

Подставляя в третье уравнение, будем иметь

 

 

 

 

 

dp

3

= −λ

 

p

 

+

 

λ λ

2

 

 

e

λ t

 

λ λ

2

e

λ

t

.

 

 

 

 

 

 

 

 

2

3

 

 

 

 

1

 

 

 

 

1

 

 

 

 

1

 

2

 

 

 

 

 

 

 

dt

λ

2

λ

λ

2

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

При p3 (0) = 0 , после интегрирования

 

 

 

 

 

 

 

 

 

 

 

 

 

p

3

(t) =

 

 

 

λ1λ2

 

 

eλ1t

 

 

 

λ1λ2

 

 

eλ2t

 

 

λ1λ2t

eλ2t .

 

(λ

2

λ )2

 

(λ

2

λ )2

 

 

(λ

2

λ )2

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

Для

 

 

нахождения

 

 

 

 

функции p4 (t)

 

 

 

 

 

воспользуемся

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p4 =1 p(t) =1( p1 + p2 + p3 ) =

 

 

 

 

 

 

 

 

 

 

=1 eλ1t

 

 

λ1

 

 

 

 

eλ1t

+

 

 

 

λ1

 

 

 

 

eλ2t

 

 

 

 

 

 

 

 

λ

2

λ

 

 

λ

2

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

λ1λ2

eλ1t +

 

 

 

λ1λ2

 

eλ2t

+

 

 

 

λ1λ2t

 

 

eλ2t .

 

 

(λ

2

λ )2

(λ

2

 

λ )2

 

(λ

2

λ )2

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

8.1. Интенсивность отказов элемента задана графически (рис. 13). Найти закон надежности p(t) и среднее время

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

Рис. 13 Решение. Из рисунка находим, что интенсивность отказов

на участке (0,1) меняется по закону λ(t) = 2 t а на участке

t >1 постоянна и равнаλ =1 .

Надежность на участке (0,1) по формуле (8) будет

41

t

 

 

t 2

 

 

(2t )dt

= e

(2t

)

.

p (t) = e 0

2

 

1

 

 

 

 

 

Найдем p(t) на участке t >1

 

 

 

 

 

t λ(t)dt =t dt = t 1;

1

1

тогда p2 (t) = et +1.

Рис.14

График закона надежности показан на рис. 31.14. Среднее время безотказной работы равно площади, ограниченной

кривыми p1 (t) , p2 (t)

и осями координат

 

1

1

 

 

t 2

 

(2t

 

 

) dt + et +1dt =1,05 +1 = 2,05.

 

 

2

t

= p1 (t)dt + p2 (t)dt = e

0

1

0

 

 

 

 

 

1

 

 

 

 

1

(2t

t 2

 

Здесь первый интеграл е

 

) dt =1,05 вычислен

2

 

 

 

 

0

 

 

 

 

 

 

приближенно.

8.2. Плотность распределения времени безотказной работы элемента на участке (t1 ,t2 ) постоянна и равна нулю вне этого

участка (рис. 15). Найти интенсивность отказов и построить график.

42

Рис. 15

Решение. Интенсивность отказов по формуле (7) равна

λ(t) =

f (t)

=

f (t)

 

 

 

 

 

(t < t < t

2

).

 

 

 

 

 

 

 

 

 

 

p(t)

1 q(t)

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

t

 

dt

 

t t1

 

 

 

q(t) = f (t)dt =

 

 

 

=

,

 

 

t

2

t

t

2

t

1

 

 

 

 

t

t

1

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

λ(t) = t2 1t .

График интенсивности отказов приведен на рис. 16.

Рис. 16 8.3. Простая система состоит из 7 независимых элементов,

надежность каждого из которых p = 0,9 .Найти надежность

системы.

Решение. Воспользуемся формулой (10), тогда

P = 0,97 0,478.

8.4. Простая система состоит из 100 одинаково надежных

43

независимых элементов. Какова должна быть надежность отдельного элемента, чтобы надежность системы была бы не менее 0,85?

Решение. Представим формулу (10) в виде p = n P . Тогда

надежность отдельного элемента

 

p = 100 0,85; ln p =

 

1

ln 0,85; p = 0,9995.

100

 

 

8.5. Простая система состоит из трех независимых элементов. Найти интенсивность отказов системы, если плотности распределения времени безотказной работы заданы выражениями:

f1 (t) =1,

 

 

f2 (t) = 3t 2 ,

 

f3 (t) =1 2t

 

 

при 0 < t <1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Поскольку f (t) = q (t) то ненадежность каждого

элемента будет

 

 

 

 

q

 

(t) = t 3

 

 

 

q

 

(t) = t t 2

 

 

 

 

 

 

 

 

 

 

 

 

q (t) = t;

 

2

;

 

3

при 0 < t <1.

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда надежность элементов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p (t) =1 t;

p

2

(t) =1 t 3 ;

 

 

p

3

(t) =1 t +t 2

при 0 < t

<1.

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интенсивность отказа каждого элемента по формуле (7)

равна

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ (t) =

 

 

1

 

; λ

2

(t) =

 

3t

2

; λ

3

(t) =

 

 

1

2t

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

t

 

 

 

 

 

 

 

 

1 t 3

 

 

 

 

1

t +t 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интенсивность отказов системы находим по формуле (11)

Λ(t) = λ

(t)

+ λ

2

(t) + λ

3

 

(t) =

1

 

+

3t 2

 

 

+

 

1 2t

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1t

 

1 t 3

 

1 t +t 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

2(1 t +

2t 2 2t 3 +3t 4 )

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 t

3 )(1 t +t 2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8.6. Найти надежность системы, состоящей из семи элементов с надежностями p1 , p2 ,..., p7 (рис. 17).

44

Рис. 17 Решение. Поскольку в системе применяются как

«последовательное», так и «параллельное» соединение элементов, то при оценке надежности расчленяем ее на ряд подсистем. Рассматривая подсистемы как условные элементы, находим надежность системы в целом.

Подсистема I — «параллельно» включенные элементы Э2 и Э3 ; надежность p1 =1 (1 p2 )(1 p3 ).

Подсистема II — «последовательно» соединенные

элементы Э1 и I; надежность pII = pI p1.

Подсистема III — «параллельно» включенные элементы Э4 и Э5 ; надежность pIII =1 (1p4 )(1 p5 ).

Подсистема IV — «последовательное» соединение элементов Э6 и Э7 , надежность pIV = p6 p7 .

Подсистема V — «последовательно» соединенные элементы III и IV; надежность pV = pIII pIV .

Вся система — «параллельно» включенные II и V элементы; надежность P =1 (1 pII )(1 pV ).

8.7. Рассмотрим технологическую линию, состоящую из одной основной машины и трех резервных. Пусть основная машина подвергается простейшему потоку отказов с интенсивностью λ1 . Найти надежность линии с облегченным

45

резервом, т. е. когда резервные машины до включения

подвергаются

простейшему

потоку

отказов

с

интенсивностьюλ2 ,

а после

включения

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

повышается до величины λ2

0 .

 

 

 

Решение.

При

определении

надежности

линии введем

следующие обозначения состояний. Если основная машина работает, то первый индекс равен нулю, если основная машина отказала, то первый индекс равен единице, а второй индекс определяет число исправных резервных машин. Так, S0k , —

основная машина исправна, из резервных исправны k машин ( k =О, 1, 2, 3); S1k — основная машина отказала, из резервных

исправны k машин, причем одна из них работает.

Для составления уравнений Колмогорова представим граф состояний (рис. 18).

Рис. 18

Система уравнений вероятностей состояний в этом случае будет

dp03 = −(λ1 +3λ2 ) p03 ;

dt

dp02 = −(λ1 + 2λ2 ) p02 +3λ2 p03 ;

dt

dp01 = −(λ1 + λ2 ) p01 + 2λ2 p02 ;

dt

dp00 = −λ1 p00 + λ2 p01 ;

dt

46

dpdt13 = −(λ2 0 + 2λ2 ) p13 + λ1 p03 ;

dpdt12 = −(λ2 0 + λ2 ) p12 + λ1 p02 + (λ2 0 + 2λ2 ) p13 ; dpdt11 = −λ2 0 p11 + λ1 p01 + (λ2 0 + λ2 ) p12 ;

dpdt10 = λ1 p00 + λ2 0 p11.

Из интегрирования первого уравнения системы имеем

p03 = e(λ1 +3λ2 )t .

Подставляя это решение во второе уравнение системы и интегрируя, находим p02 .Продолжая последовательно процесс

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

p10 (t) =1 ( p03 + p02 + p01 + p00 + p13 + p12 + p11 ).

Отсюда надежность линии, как обратное событие, будет

P(t) =1p10 (t).

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

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

9. ЗАДАЧИ ДЛЯ САМОСТОЯТЕЛЬНОГО РЕШЕНИЯ

Найти глобальные экстремумы следующих функций:

1.1.

y = ax2 +bx +c, (a 0), x R ;

1.2.

y =

 

 

x

, x R ;

1

+ x2

 

 

 

1.3.

 

y = −3x4 +6x2 , x [2, 2] ;

47

1.4.

y = xe

x2

,

x R ;

2

 

 

 

1

 

 

e

 

( xa)2

1.5.

y =

 

 

 

2

, x R , a = const .

2π

 

 

 

 

 

 

 

 

1.6.Определить размеры открытого бассейна с квадратным дном объемом 32 м3 так, чтобы на облицовку его стен и дна пошло наименьшее количество материала.

1.7.Для хранения и охраны военного имущества требуется огородить забором прямоугольный участок площадью S=1,5 га и затем таким же забором разделить этот участок на две равные части. Определить размеры участка, при которых расход материалов будет наименьшим.

1.8.Штаб инженерной части планирует строительство

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

1.9.Два самолета летят в одной плоскости и

прямолинейно под углом 1200 с одинаковой скоростью υ км/час. В некоторый момент один самолет пришел в точку пересечения линий движения, а второй не дошел до неё d км.

Через сколько времени расстояние между самолетами будет наименьшим? Чему равно это расстояние?

1.10. Из трех досок одинаковой ширины сколачивается желоб. При каком угле наклона боковых стенок площадь поперечного сечения желоба будет наибольшей?

2.1. Найти экстремумы функции:

a) u = x2

- xy + y2 - 2x +y;

б) z = ( x2 + y2 ) e-(x2 +y2 ) ;

в)

z = ( x - 1 )2 + 4 - 2y2 ;

г)

z = ex-y

( x2 - 2y2 );

д) z = 2 - xy - 4x - 2y;

48

е) z = x3 + 8y3 - 6xy + 1;

ж) z = 3x2 - 2x y + y - 8x + 8; з) z = x2 - xy + y2 + 9x - 6y + 20; и) z = y x - y2 - x + 6y;

к) z = 3xy - x2 y - xy2 .

2.2.Найти экстремумы функций трех переменных:

а) u = x2

+ x2

+ x2

- x

x

2

+ x - 2x

3

;

1

 

2

 

3

1

 

1

 

б) u = x +

y2

 

+

z2

+

2

,

( x > 0, y > 0, z > 0);

4x

y

z

 

 

 

 

 

 

 

 

в) u = x - y2z3 ( 1 - x - 2y - 3z ), (x > 0, y > 0, z > 0).

2.3.Найти экстремум функции

z = x + 2 y

при условии x2 + y2 = 5

2.4.x2 + y2 extr; 3x + 4y = 1.

2.5.exy extr, x + y = 1.

2.6.

5x2

+ 4xy

+ y2

extr; x + y = 1.

2.7.

3x2

+ 4xy

+ y2

x extr;

x + y = 1.

2.8.

а)

1

+

1

extr;

 

1

+

 

1

 

=

1

.

x

y

 

x2

 

y2

 

 

 

 

 

 

 

 

 

 

a2

 

 

 

n

 

 

 

n

 

 

 

 

 

 

 

 

б)

xi2 extr;

xi4

1.

 

 

 

 

 

i=1

 

 

 

i=1

1

+ 1

 

1 = 1.

2.9.

x

+ y

+ z

extr;

+

 

 

 

 

 

 

 

 

 

x

 

y

 

z

2.10. a cos2 x + b cos2 y extr; y - z = π4 .

49