Смекни!
smekni.com

Исследование математических операций 2 (стр. 26 из 28)

Считают, что интервалы времени между отказами

t1k,
t2k, ....
tpk представляют собой систему взаимно независимых случайных величин с плотностями распределения наработок между отказами f(
t1), f(
t2) , … , f(
tp) .

Моменты отказов или восстановлений образуют в каждом k-м опыте (испытании) ряд чисел по следующему правилу:

(5.4)

или

.

где tik - время работы (наработка) элемента до i-го отказа в k–м опыте, час,

,
.

tik - время работы (наработка) элемента между (i-1)-м и i-м отказами в k–й реализации, час,
,

.

Числа t1k, t2k, ... ,tpk образуют случайный поток, который назы­вается процессом восстановления. Этот процесс является различным для различных элементов и продолжается до окончания срока служ­бы системы. Изучением таких процессов занимается теория восста­новления.

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

 простой, при котором все функции распределения наработок до первого и между последующими отказами Fi (t) равны;

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

, i = 2, 3, 4,…

 сложный, при котором все функции распределения Fi (t) раз­личны.

Основной характеристикой процесса восстановления является функция восстановления

и ее дифференциальная характерис­тика - плотность восстановления
, определяемые по следую­щим формулам:

; (5.5)

; (5.6)

где Fn(t) и fn(t) - соответственно плотность и функция распределения на­работки до n-го отказа.

В случае независимости наработок между отказами функции распределения Fn(t) наработок до n-го отказа находятся путем по­следовательного применения правила свертки для суммы двух слу­чайных величин:

; (5.7)

F1(t) = F(t) .

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

и
по формулам (5.5), (5.6) состоит в том, что свертка (5.7) лишь для некоторых законов распределения вычисляется в конечном виде. Использование аналитических методов расчета плотности
и функции восстановления
ограничено из-за сложности математической формализации применяемых стратегий восстановления работоспособности технических систем и необходимости учета множества факторов, влияющих на замену элемента в системе. В этих условиях наиболее эффективным методом расчета
и
является метод Монте-Карло.

Расчет ведущей функции и параметра потока отказов этим методом в случае простого, общего или сложного процессов производится в следующем порядке.

По известным законам распределения наработок элементов с использованием формул преобразования (табл. 5.1) моделируются массивы случайных величин

tik между (i-1)-м и i-м отказами Размерность каждого массива равна N.

Далее вычисляются значения наработок до i-го отказа tik по следующим формулам:

; (5.8)

, (5.9)

где i – номер отказа,

k – номер реализации при моделировании,

p – максимальное число отказов элемента, получаемое в k-й реализации случайного процесса

Затем полученные случайные величины наработок tik группиру­ются по интервалам времени.

Номера интервалов, в которые попадают моменты возникнове­ния отказов t1k, t2k, ... ,tik, ... , tpk определяются по формуле:

, (5.10)

где

- наименьшее целое число, не меньшее
;

ti - величина интервала времени

Параметр и ведущая функция потока отказов в j-м интервале времени определяется по следующим формулам:

; (5.11)

; (5.12)

где nij - число попаданий случайной наработки до i-го отказа tik в j-й ин­тервал времени (

) за N реализаций.

; (5.13)

. (5.22)

где h - максимальное число интервалов времени.

Пример 5.5. Законы распределения наработок элемента системы до первого и второго отказов и соответствующие параметры этих законов приведены в следующей таблице:

№ отказа

Закон распределения

Параметры закона

a(

)

b

1

Вейбула

1,4

45,8

2

Экспоненциальный

0,3

-

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

ti = 1 час).

Решение

1. Выберем равномерно распределенное случайное число. До пустим

i = 0,725.

2. Вычислим случайные значения наработок на отказ элемента используя формулы табл. 5.1.

час.;

;

час.;

час.

3. Определим номер временного интервала, на котором произойдут отказы

первый отказ

;

второй отказ

.

В ходе первой реализации элемент системы первый раз откажет на 21-м временном интервале, а второй отказ произойдет на 22-м временном интервале.

В ходе первой реализации элемент системы первый раз откажет на 21-м временном интервале, а второй отказ произойдет на 22-м временном интервале.

Назад | Содержание | Далее

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

Задача 5.1