- •А.Ю. Крюков, Б.Ф. Потапов
- •Крюков, А.Ю.
- •СОДЕРЖАНИЕ
- •ПРЕДИСЛОВИЕ
- •ВВЕДЕНИЕ
- •1.1. СУЩНОСТЬ И НАЗНАЧЕНИЕ МОДЕЛИРОВАНИЯ
- •1.1.1. Понятие о моделировании и свойства моделей
- •1.1.2. Цели моделирования
- •1.1.3. Виды моделирования и классификация моделей
- •1.2.1. Общие положения и основные определения
- •1.2.2. Структура математической модели и ее построение
- •1.2.3. Иерархия математических моделей
- •1.2.5. Классификация математических моделей
- •1.2.6. Геометрическое представление математических моделей
- •1.3. ЗАКЛЮЧИТЕЛЬНЫЕ ЗАМЕЧАНИЯ
- •2.1. МОДЕЛИ МИКРОУРОВНЯ
- •2.2. МОДЕЛИ МАКРОУРОВНЯ
- •2.2.1. Общая характеристика моделей, их структура и сущность
- •1. Компонентные уравнения
- •2. Топологические уравнения
- •Компонентные и топологические уравнения электрической системы
- •2.2.5. Задания для самостоятельной работы
- •2.3. ЗАКЛЮЧИТЕЛЬНЫЕ ЗАМЕЧАНИЯ
- •3.2.1. Общие принципы моделирования полей
- •3.2.2. Особенности построения моделей
- •3.2.3. Модели стационарных полей
- •3.2.4. Модели нестационарных полей
- •3.2.7. Задание для самостоятельной работы
- •3.3. ЗАКЛЮЧИТЕЛЬНЫЕ ЗАМЕЧАНИЯ
- •4.1.1. Надежность объектов как комплексное свойство
- •4.1.2. Классификация отказов и временные понятия
- •4.3.1. Основные понятия, определения и положения
- •4.3.2. Основные характеристики случайных величин
- •4.4. ПОКАЗАТЕЛИ НАДЕЖНОСТИ
- •4.4.1. Показатели безотказности
- •4.6.1. Общая характеристика и виды моделей
- •4.6.2. Распределения, используемые в теории надежности
- •4.6.3. Композиция законов распределения
- •4.6.4. Задание для самостоятельной работы
- •4.7 ПОТОКИ ОТКАЗОВ И ВОССТАНОВЛЕНИЙ
- •4.8. РАСЧЕТ НАДЕЖНОСТИ ТЕХНИЧЕСКИХ СИСТЕМ
- •4.8.1. Общие положения
- •4.8.2. Модели параметрических отказов
- •5.1. ОСНОВЫ ТЕОРИИ ГРАФОВ И СЕТЕВЫЕ МОДЕЛИ
- •5.1.1. Общие замечания, основные понятия и определения
- •5.1.2. Решение задачи о максимальном потоке
- •5.1.3. Потоки минимальной стоимости
- •5.1.4. Практические примеры сетевых задач
- •5.1.5. Некоторые обобщения по сетевым задачам
- •5.1.6. Альтернативные методы решения сетевых задач
- •5.2.1. Основные положения
- •5.2.2. Структура принятия решения
- •5.2.4. Пример применения классических критериев
- •6.1.1. Постановка задач оптимизации и критерии
- •оптимальности
- •6.1.2. Многокритериальные задачи оптимизации
- •6.1.3. Классификация методов оптимизации
- •6.2. РЕГУЛЯРНЫЕ МЕТОДЫ ОПТИМИЗАЦИИ
- •6.2.1. Методы математического анализа
- •6.2.2. Понятие о вариационном исчислении
- •6.2.3. Принцип максимума Понтрягина
- •6.2.4. Метод динамического программирования
- •6.3. ПРЯМЫЕ МЕТОДЫ ОПТИМИЗАЦИИ
- •6.3.2. Минимаксные стратегии одномерного поиска
- •6.3.4. Многомерные методы безусловной оптимизации
- •6.3.7. Заключение по прямым методам оптимизации
- •7 1. ОБЩИЕ ПОЛОЖЕНИЯ
- •ЗАКЛЮЧЕНИЕ
- •Учебное издание
- •Учебное пособие
U - система подмножеств, выделенная в Q; р - числовая функ
ция, определенная на U и |
обладающая свойствами р(А) > О, |
p(fl) = 1 для произвольной |
последовательности попарно несо |
вместимых событий А„ 1,2,... (А[ГлА) = 0 при / * j, А,е U). Функ ция р называется вероятностью.
Читателю рекомендуется самостоятельно ознакомиться с представленным подходом к описанию природы вероятности, который оказался исключительно плодотворным для развития новых разделов теории вероятностей, в частности - теории слу чайных процессов.
4.4. ПОКАЗАТЕЛИ НАДЕЖНОСТИ
Показатели надежности - это количественная характе ристика одного или нескольких свойств, составляющих надеж ность объекта. Если показатель надежности характеризует одно из свойств надежности, то он называется единичным, если же несколько свойств —комплексным показателем надежности.
Рассмотрим единичные показатели (рис. 4.8), которые иг рают важную роль в проектных расчетах и отражают основные характеристики надежности технических систем.
4.4.1.Показатели безотказности
1.Вероятность безотказной работы (ВБР) - это вероят ность того, что в пределах заданной наработки отказ объекта не возникнет. ВБР является основной количественной характе ристикой безотказности объекта на заданном временном интер вале и наиболее полно определяет его надежность. Если обозна чить через Т время непрерывной исправной работы объекта от начала работы до первого отказа, а через t - время, за которое необходимо определить ВБР, то ВБР записывается в виде
P(t) = P{T> t},t> 0. |
(4.21) |
Рис. 4.8. Единичные показатели надежности
Функция ВБР обладает следующими очевидными свойст
вами:
1)0 </>(/) < 1;
2) функция P(t) есть невозрастающая функция своего ар
гумента, т.е. если h> t\, то P(h) < P(t\);
3 ) Р(0)=1,Д<ю) = 0.
Из третьего свойства следует, что вероятность безотказной работы определяется в предположении, что в начальный момент времени исчисления наработки объект был работоспособен.
Статистически ВБР равна
|
|
IIAI |
|
|
P(t) = |
lim - |
/=1__^ |
m |
(4.22) |
v 7 |
Д/—>0 |
дг |
Nn ‘ |
|
|
Nt)->со |
0 |
|
|
где N0 - число объектов в начале испытаний; и, - |
число отка |
завших объектов в интервале времени A/,; t - время, для которо
го определяется ВБР; N(t) - число объектов, исправно работаю щих на интервале [0, О-
Вероятность того, что отказ объекта произойдет за время,
не превышающее заданной величины t, т.е. что Т < t, равна |
|
6(0 = Q{T <t} = \ - Р(0, t > о. |
(423) |
Функция 6 (0 представляет собой интегральную функцию |
|
распределения случайной величины - времени отказа |
Если |
функция 6 (0 дифференцируема, то ее производная есть диффе ренциальный закон (плотность) распределения случайной вели чины t - времени исправной работы. Тогда можно записать:
« W - Л О и М - - Ш . |
(4.24) |
|
dt |
at |
|
Статистически вероятность отказа равна: |
|
|
1Ш |
|
|
2 > . |
N0 ~N(t) |
|
|
(4.25) |
|
г |
Nn |
|
м,->=о Jvo |
|
|
Графически функции P(t), Q(t) и /0 представлены на рис. 4.9.
Рис. 4.9. Графики функций:
1 - вероятности безотказной работы;
2 - вероятности отказа;
3 - плотности вероятности отказа
Отметим некоторые общие свойства функции распределеиия Q(t), которые, как следует из выражений (4.24), противопо ложны свойствам функции распределения ВБР P{t):
1)0 < 2 ( 0 ^ 1;
2)функция распределения Q{t) есть неубьтающая функЦия своего аргумента, т.е. при /2 > t\ П Ь ) > Ш ,
3)в начальный момент времени функция распределения
равна нулю: (2(0) = 0, |
и в |
бесконечности |
она равна единице |
|||
1. |
|
|
|
|
|
|
Проиллюстрируем |
пред |
о |
|
|
|
|
ставленные свойства. |
Для |
этого |
т |
t0 |
t |
|
будем рассматривать величину tQ |
|
|
|
|
||
как случайную точку |
Т на оси / |
|
m |
|
|
|
(рис. 4.10), которая в |
результате |
Рис. 4.10. Иллюстрация |
||||
опыта может занять то или иное |
свойств функции Q(t) |
|
||||
положение. |
|
|
|
|
|
|
Тогда функция |
распределе |
|
|
|
|
ния Q(t) есть вероятность того, что случайная точка Т в результа те опыта попадает левее точки to.
Будем увеличивать t, т.е. перемещать точку to вправо по оси абсцисс. Очевидно, при этом вероятность того, что случайная точка Т попадает левее to, не может уменьшиться, следовательно, функция распределения Q(t) с возрастанием t убывать не может.
Из введенного в определение вероятности безотказной ра боты предположения, что в начальный момент времени исчисле ния наработки объект был работоспособен, следует, что Q(0) = 0.
Плотность вероятности fit) статистически определяется по
формуле |
|
|
|
Д Л ') = |
An(At) |
(4.26) |
|
N0At ’ |
|||
|
где An(At) —число отказов за интервал времени At.
Очевидно, что |
|
( |
со |
6(0 = J/COrfT И P(t)= |
(4-27) |
где т - момент времени, предшествующий t.
Поскольку основной параметр функционирования объекта в теории надежности - время t - изменяется в пределах [0; со], то из условия нормировки (4.10) и первого свойства функции Q(t)
следует, что для функцииД/) имеет место равенство |
|
оо |
|
J/(0dt = 1 |
(4.28) |
2. Средняя наработка до отказа - это момент первого порядка (математическое ожидание) наработки объекта до пер вого отказа. Эту величину обозначают Гер и называют также средним временем безотказной работы:
(4.29)
Статистически средняя наработка до отказа однотипных объектов равна
|
I |
ср |
(4-30) |
*'0 7=1
где tj - время исправной работыJ-го объекта.
3. Гамма-процентная наработка до отказа (у - это нара ботка, в течение которой отказ объекта не возникнет с вероятно стью у, выраженной в процентах. Гамма-процентная наработка определяется из уравнения
(4-31)
При у =100% гамма-процентная наработка называется установленной безотказной наработкой, при у = 50 % гамма процентная наработка называется медианной наработкой.
4. Средняя наработка на отказ —отношение наработки восстанавливаемого объекта к математическому ожиданию числа его отказов в течение этой наработки. Для интервала наработки (<i; /2) средняя наработка на отказ определяется выражением
Т — |
*2 |
(4.31) |
*л — |
|
Q (/,;0
где То - средняя наработка на отказ; Q(/i; /2) - математическое ожидание числа отказов на интервале наработки /2 - Л.
То характеризует наработку восстанавливаемого объекта, приходящуюся в среднем на один отказ, в рассматриваемом ин тервале суммарной наработки или определенной продолжитель ности эксплуатации.
Статистически средняя наработка на отказ определяется так:
< « 2)
П /=1
где tcp- время исправной работы между ( / - 1)-м и i-м отказами объекта; п —число отказов объекта.
При достаточно большом числе отказов /ср будет стре миться к среднему времени между двумя соседними отказами. Если испытания проводятся не с одним, а с несколькими одно типными объектами, то среднее время между отказами можно определить из выражения
1 м |
(4.33) |
|
Т° = м Р с ’ |
||
|
где М - число объектов.
5. Интенсивность отказов - это отношение числа отка завших объектов в единицу времени к среднему числу объектов, продолжающих исправно работать в данный интервал времени.
Другими словами, интенсивность отказов - это условная плотность вероятности возникновения отказа невосстанавливаемого объекта, определяемая для рассматриваемого момента времени при условии, что до этого момента отказ еще не возник.
Математически интенсивность отказов определяется сле дующим образом:
4 0 = |
g(t;t + At) |
(4.34) |
|
|
At |
где q(t; t + At) - вероятность^ отказа невосстанавливаемого объ екта на интервале времени от t до / + At при условии, что до мо мента времени t отказ не возникает.
Определение интенсивности отказов на основе статисти ческих данных осуществляют по формуле
& , |
(4.35) |
N(t)At
где An(At) - число отказов объекта за период времени ( —— ;
t + At
2
Интенсивность отказов часто называют А,-характеристи- кой, она показывает, какая часть объектов выходит из строя в единицу времени по отношению к среднему числу исправно работающих объектов.
Функция Х,(/) для многих объектов имеет характерный вид (рис. 4.11) для разных периодов эксплуатации.
Период I (О < / < Л ) называют периодом приработки, в процессе которого устраняются систематические отказы эле ментов (конструктивные, производственные, эксплуатацион ные), заменяются элементы, имеющие дефекты. Все это приво дит к повышению вероятности безотказной работы и, следова тельно, снижению интенсивности отказов.
Рис. 4.11. Интенсивность отказов
Встречаются случаи, когда на начальной стадии прира ботки интенсивность отказов возрастает (участок Г). Такие зави симости могут быть описаны в следующем виде:
= |
1. |
(4-36) |
кк
Эта модель описывает смесь элементов с вероятностями безотказной работы Pk(t) и относительными долями р Группа дефектных элементов имеет малую среднюю наработку до отка за и быстрое старение.
В периоде II (t\ < t < /2), часто называемом периодом нор мальной эксплуатации, интенсивность отказов остается пример но постоянной. Для периода нормальной эксплуатации харак терно наименьшее количество отказов.
Наконец, в периоде III, называемом периодом старения и износа, интенсивность отказов увеличивается из-за необрати мых физико-химических процессов в элементах объекта, свя занных с его длительным использованием, т.е. когда вследствие естественных процессов объект или его составные части при ближаются к предельному состоянию по условиям физического
изнашивания. Отказы III периода иногда называют деградации отыми отказами.
6. Параметр потока отказов - отношение среднего числа отказов восстанавливаемого.объекта за произвольно малую его наработку к значению этой наработки.
Параметр потока отказов со(/) характеризует среднее число отказов, ожидаемых на малом интервале времени, и равен
dCljt)
(4.37)
dt ’
где Q(t) - математическое ожидание числа отказов за время /. Отметим, что средняя наработка до отказа, гамма
процентная наработка до отказа, интенсивность отказов харак теризуют надежность невосстанавливаемых объектов. Эти пока затели могут быть использованы и для характеристик надежно сти восстанавливаемых объектов, но при условии, что наработка восстанавливаемых объектов рассматривается только лишь до первого отказа. Средняя наработка на отказ и параметр потока отказов используются в качестве показателей для восстанавли ваемых систем.
Для того чтобы определить показатели, характеризующие надежность восстанавливаемого элемента, необходимо определить модель его функционирования. Простейшей является такая мо дель, при которой объект какое-то случайное время Т| работает до первого отказа, фиксируемого достоверно, затем следует мгновен ное и полное восстановление свойств объекта в момент Х\ = 1\, по сле чего элемент снова работает случайное время х2до второго от каза, определяемого достоверно, затем мгновенно восстанавлива ется до начального состояния в момент t\ = TJ + х2и т.д.
Для таких объектов моменты отказов на оси суммарной наработки или на оси непрерывного времени образуют поток отказов. В качестве характеристики потока отказов используют «ведущую функцию» Q(T) данного потока - математическое ожидание числа отказов за время t.
Статистически параметр потока отказов можно опреде лить по формуле
©(/) = An(At) |
(4.38) |
N0At |
’ |
t —At
где An(At) - число отказов объекта за период времени
2 *
t + At
~7 ~ }'
4.4.2.Показатели долговечности, ремонтопригодности
исохраняемости
Данная группа показателей вводится аналогично показа телям безотказности.
Средний ресурс Тр- это математическое ожидание ресурса.
Гамма-процентный ресурс Т^% - это наработка, в течение которой объект не достигнет предельного состояния с заданной вероятностью у, выраженной в процентах. Гамма-процентный ресурс определяется по формуле, аналогичной (4.31).
Назначенный ресурс Гр„ определяется как суммарная на работка объекта, при достижении которой применение по на значению должно быть прекращено.
Средний срок службы Тс - это математическое ожидание срока службы.
Гамма-процентный срок службы Тс% характеризуется ка лендарной продолжительностью от начала эксплуатации объек та, в течение которой он не достигнет предельного состояния с заданной вероятностью у, выраженной в процентах.
Назначенный срок службы Тсли - это календарная про должительность эксплуатации объекта, при достижении которой применение по назначению должно быть прекращено.
Назначенный ресурс (срок службы) - это характеристики, устанавливаемые на основании субъективных или организаци
онных принципов и поэтому являющиеся косвенными показате лями надежности.
Основным показателем ремонтопригодности является ве роятность восстановления. Момент восстановления работо способности объекта после отказа является случайным событи ем. Поэтому интервал времени от момента отказа до момента восстановления является случайной величиной и для характери стики ремонтопригодности может быть использована функция распределения этой случайной величины ©. Вероятностью вос становления Ря называется вероятность того, что время восста новления работоспособного состояния объекта не превысит за данного:
PB(t) = P {® < l},0< t. |
(4.39) |
Функция Рв(0 представляет собой интегральную функцию распределения случайной величины 0. Вероятность невосста-
новления на заданном интервале (, т.е. вероятность |
того, что |
0 > /, равна |
|
6 .(0 = Р{® > t) = 1 - Р Ш о <1. |
(4.40) |
Плотность вероятности момента восстановления равна: |
|
dP |
(4.50) |
/«(0 = - ^ - |
По аналогии со средней наработкой до отказа момент пер вого порядка m\{Q} (математическое ожидание) времени вос становления работоспособного состояния объекта называется
средним временем восстановления:
(4.51)
|
0 |
|
|
II |
ft? 1 |
/""'Ч |
(4.52) |
|
|
|
о
Важным показателем ремонтопригодности объекта явля ется интенсивность восстановления р(/), которая, следуя общей
методологии, аналогична показателю безотказности —интенсив ности отказов.
Показатели сохраняемости - средний срок сохраняемо сти и гамма-процентный срок сохраняемости - определяются аналогично соответствующим показателям безотказности и дол говечности. Средний срок сохраняемости - это математическое ожидание срока сохраняемости, а гамма-процентный срок со храняемости это срок сохраняемости, достигаемый объектом с заданной вероятностью у, выраженной в процентах.
4.5. АНАЛИТИЧЕСКИЕ ЗАВИСИМОСТИ МЕЖДУ
ПОКАЗАТЕЛЯМИ НАДЕЖНОСТИ
Между основными показателями надежности существуют аналитические зависимости. Приведем некоторые из них.
1.Зависимость между вероятностью безотказной работы
исредней наработкой до отказа. По выражению (4.29)
TCf = \tf(t)dt- |
(4.53) |
О |
|
Так как/ 0 = Q\t) и Q(f) = 1 -P{t), то, выполнив интегри
рование по частям, получим
Тср = -\tdP (t) = -tP (tf0 +)p(t)dt |
(4.55) |
|
О |
О |
|
Учитывая свойства функции P(t), получаем равенство ну лю первого слагаемого в правой части полученного равенства.
Тогда
Tcp = \P{t)dt. |
(4.55) |
2. Связь между ВБР и интенсивностью отказов. Согласно свойствам функций P(t) и Q(t) получим:
q(t, t + ДО = 1 - P(t-1 + At). |
(4.56) |
С другой стороны вероятность безотказной |
работы |
P(t + At) на интервале времени от 0 до t+At в соответствии с тео ремой умножения вероятностей рассчитывается по формуле
|
Pit + At) = P(t)P(t; t + At). |
(4.57) |
|||
Графическая запись приведенных рассуждений представ |
|||||
лена на рис. 4.12. |
|
|
|
|
|
P(t) |
Pit,t+At) |
|
С учетом (4.34), (4.56) и (4.57) |
||
имеем: |
|
|
|||
|
|
|
|
||
0 |
t t+At t |
|
|
j Pit + At) |
|
P it+ At) |
|
40 = |
Pit) |
(4.58) |
|
Puc. 4.12. Иллюстрация |
|
|
At |
||
|
P it)-P it + At) |
|
|||
теоремы умножения |
|
|
|||
|
|
PiOAt |
|
||
вероятностей |
|
|
|
||
|
|
|
|
||
|
|
|
Переходя к пределу при At -> О, |
||
получаем |
|
|
|
|
|
|
|
1 |
dPjt) |
Л О |
(4.59) |
|
4 0 = Pit) |
dt |
Pit) |
||
После интегрирования выражения (4.59) при начальном |
|||||
условии Pi0) = 1 получаем: |
|
|
|
||
|
|
Pit) =e " |
|
(4.60) |
где т - момент времени, предшествующий t.
Различие между интенсивностью отказов \(t) и плотно стью распределения отказов fit) заключается в том, что эле мент условной вероятности Xit)dt характеризует вероятность отказа объекта на элементарном участке (/; t + dt), если из вестно, что к моменту времени / он не отказал. Элемент веро-
ятностиХО** характеризует вероятность отказа в течение эле ментарной наработки на участке (t; t + dt), но при этом неиз вестно, исправен ли объект к моменту времени t. Кроме того, из (4.59) следует, что fit) = X(t)P(t), то есть при малых t значе ние P{t) близко к единице и, следовательно, fit) « X(t). При
t =0 P(t) - 1 |
и fit) = Ц/). Во всех остальных случаях, когда |
p { t)< \,fit)< m - |
|
3. |
Связь между вероятностью безотказной работы, интен |
сивностью отказов и средней наработкой до отказа.
Среднюю наработку до отказа можно вычислить через ин
тенсивность отказов. Подставив в выражение (4.55) значение (4.60), получим
(4.61)
4. Зависимость между плотностью вероятности времени безотказной работы и параметром потока отказов.
Пусть в момент t = 0 на испытании находится N0объектов. По мере выхода из строя отказавшие объекты заменяются но выми (выборка с возмещением), если объекты невосстанавливаемые. Тогда на основании расчета среднего числа отказавших объектов в зависимости от плотности распределения отказов^/)* параметра потока отказов со(/) и длительности рассматриваемого интервала времени можно получить уравнение Вольтерра, свя зывающее со(/) иfit):
(4.62)
О
где т - момент времени, предшествующий /, а интеграл в правой части есть число отказавших объектов из числа замененных в процессе испытаний за время от 0 до /.
Решение уравнения (4.62) показывает, независимо от вида функции f(t), т.е. независимо от закона распределения времени безотказной работы, что параметр потока отказов объектов стремится к постоянной величине - обратной средней наработке на отказ, т.е. справедливо выражение
lim©(0 = -Jr. |
(4.63) |
т0 |
|
5. Связь между вероятностью восстановления и интенсив ностью восстановления.
Рассмотрим условную вероятность PB(t/т) того, что вос становление работоспособности объекта произойдет на интервале времени следующим за интервалом х после отказа, на ко тором еще не удалось восстановить работоспособность:
р, « / *)= |
РЛ‘+ *)-РщЮ |
(4.64) |
|
1 - а д |
|||
|
|
Предел отношения PB(tlx) при t —» 0 есть не что иное, как дифференциальная плотность вероятности восстановления в момент т при условии, что объект не был восстановлен до мо мента х:
PB(tk )_ d P B(x) |
1___ |
|
р(т) = lirn |
|
(4.65) |
f->О l |
dl 1 - P B( T ) |
|
Представляя (4.65) в виде |
|
|
р(т) = - — 1п[1-Р(т)] |
(4.66) |
и учитывая, что Рв(0) = 0, получаем, аналогично (4.60),
/
- | и ( т ) А
PB(t) = \ - e " |
(4.67) |