книги / Надежность и диагностика технологических систем
..pdf3.1. Законы распределения случайных величин |
91 |
На рис. 3.3 приведены кривые распределения Пуассона при различных значениях а, которые показывают, что с ростом т
укривых распределения Пуассона появляется максимум, а при
т> 6 кривые по характеру все больше становятся похожими на известную кривую нормального распределения закона Гаусса.
Рис. 3.3. Распределение Пуассона при различных значениях параметра а
Закон Гаусса выводится из биноминального закона при стрем лении п —» ©о и т —» ©о произвольным образом. В теории вероятно стей доказывается, что при х = (т - np)/(npq)0,5вероятностьрт(П) оценивается по формуле
Р т ( п ) |
1 |
(3.11) |
|
J2nnpq |
|||
|
|
Данное выражение только незначительно отличается от рас пределения Гаусса, записанного в нормализованной форме:
У = |
1 g-*72 |
(3.12) |
№ |
|
Нормальный закон относится к классу предельных. К виду этого закона приближаются при определенных условиях так на зываемые смеси распределений, описываемые другими законами при их многократном употреблении (суммировании). Имеются
92 |
3. Математические основы теории надежности |
и некоторые другие законы, которые при определенных услови ях можно рассматривать как предельные.
Плотность вероятности случайной величины, распределенной по нормальному закону, в общем случае имеет вид:
1 |
с -(х-т х )2/(2Р2) |
(3.13) |
«* ) = |
|
|
оч/27г |
|
В данном выражении параметр тх соответствует математиче скому ожиданию, а параметр а = — среднему квадратиче скому отклонению распределения случайной величины.
Вид кривой плотности распределения по нормальному закону имеет вид симметричной колоколообразной кривой с максиму мом, соответствующим математическому ожиданию при х = тх (рис. 3.4). Очевидно, что при таких условиях величина макси мума равна f(x)шаХ = 1/(а-^/2 п ). Чем меньше величина среднего квадратического отклонения а, тем больше величина f(x)max, и наоборот.
Рис. 3.4. График плотности распределения функции Лапласа
Функция распределения для нормального закона равна:
1 |
(.x-mx?l(2<?)dX' |
(3.14) |
|
а^2п к |
|||
|
|
Несмотря на то что данная функция распределена по всей числовой оси, вероятность события на интервале от до тх — За составляет всего рх = 0,00135, а за пределами интервала (тх - За < < х < тх + За) равна рх = 0,0027. Традиционно принято считать, что случайная величина, распределенная по нормальному закону,
3.1. Законы распределения случайных величин |
93 |
находится в диапазоне «трех сигм», т.е. (тх - За < х < тх + Зо). Однако можно принять и другой диапазон.
Интеграл формулы (3.14) не выражается в элементарных функциях. Поэтому вероятность попадания случайной величи ны на участок от А до В определяется выражением
Если сделать замену переменной интегрирования по формуле t = (x -m x)l<3*j2, то новые верхний и нижний пределы соответст венно будут равны: |3 = (В -тпх)/о>/2; a = (A -m x)/aV2, а величина дифференциала dx = <W2dt.
Тогда формула (3.15) примет вид:
(3.16)
Для вычисления даннрго интеграла можно провести интег рирование численным методом на компьютере, например в сис теме MathCAD, или воспользоваться таблицами интеграла Ф(*). При практических вычислениях чаще используют табулирован ный интеграл (функцию Лапласа):
Данный интеграл имеет следующие свойства:
1.Ф(0) = 0, поскольку пределы будут совпадать.
2.Ф(оо) = 0,5, так как интегрируют только на положительной полуоси.
3.Функция Ф(£) нечетная, т.е. Ф(-£) = -Ф(*)» так как при за мене у интеграла в верхнем пределе х на -х получается отрица тельное значение определенного интеграла.
Рассмотрим участок функции плотности нормального распре деления от точки п (с отклонением от математического ожида ния 1г) до точки <7 с отклонением 12 (см. рис. 3.4). Вероятность попадания на этот участок равна:
Р(п,д) = Ф[*х/а] - Ф[г2/а].
94 |
3. Математические основы теории надежности |
Вероятность попадания на участок между точками о и q опре делится из выражения
р (о, q) = Ф (ф ) -Ф (-V o ) = Ф (ZI/ G) +Ф (k/o). (3.18)
Вероятность нахождения случайной величины в диапазоне от - 00 до точки k находится по формуле
р (-оо, к) = Ф (~ ф )~ Ф (— /о) =
= ф (со)- ф (Z4/o) = 0,5-Ф (Z4/ G). |
(3.19) |
Пример 3.5. Определить вероятность появления отказа на участке от -«в до точки k, если 14= 8, а величина среднего квадратического от клонения о = 4.
Решение.
По формуле (3.19) находим: Р(-оо, k) = 0,5 - Ф(/4/о) = 0,5 - Ф(8/4) = = 0,5 - Ф(2) = 0,5 - 0,4772 = 0,0228. Значение 0,4772 — табличное зна чение величины Ф(2).
При исследовании надежности используются и другие законы распределения случайных величин. Например, экспоненциаль ное распределение, функция распределения которого описывается формулой (3.10), имеет плотность распределения
f(x) = F\x) = A,exp(-Ajc).
Математическое ожидание
©о |
оо |
ОО |
тх = Jxf(x)dx =|xX-exp(-Xx)dx = A.J xexp(-Xx)dx. |
||
о |
о |
о |
После интегрирования по частям получим |
||
|
тх = 1/А. |
|
Дисперсия данного распределения |
|
|
оо |
оо |
оо |
DX = J x?f(x)dx =fx2X-exp(-Xx)dx = Xj я2 •exp(-Xx)dx. (3.20)
о о о
При интегрировании по частям получим значение D(x) = 1/Х2 или G(X ) = 1/Х. При экспоненциальном распределении тх = а(х) = = 1/Х.
3.1. Законы распределения случайных величин |
95 |
Экспоненциальное распределение применяется для опреде ления характеристик надежности элементов электронной аппа ратуры ТС и других элементов, а также в случаях, когда ТС эксплуатируется очень долго и многократно ремонтируется.
Логарифмическое нормальное распределение имеет смещен ный максимум и очень удобно при описании законов распреде ления отказов различных элементов, например подшипников качения и скольжения. Плотность распределения в этом случае находится по формуле
/(*) = — |
(з. 21) |
хаы2п
Математическое ожидание и дисперсия данной функции равны:
тх = ехр (а + а2/^);
Dx = ехр 2 (а+ сз2)[ехр(о2) -1].
Закон Вейбулла имеет примерно такой же вид кривой плот ности вероятности, как и предыдущий:
f(x) = ЪХ•дгь_1ехр (-Хя&). |
(3.22) |
Этим законом удобно пользоваться при исследовании надеж ности изнашивающихся или стареющих элементов. Матема тическое ожидание и дисперсия данного закона находятся по формулам:
тх = 1/суаТ(1+1/а);
Dx = l/c ^ [Г(1+ 2/а)- Г2(1- 1 /а)].
В этих формулах функция
оо
Цх) = jt^e-'dt
является гамма-функцией, табулированной в специальных таб лицах. Закон Вейбулла удобно использовать при описании не симметричных распределений, сдвинутых в сторону и с одним максимумом.
Гамма-распределение применяется в аналогичных случаях. Его функция плотности вероятности (в отдельных случаях при t = 2Хх) имеет вид:
f(t) = l/[2 4Cn -1)]*4exp(-f/2).
3.2. Классификация технических систем и элементов |
97 |
•выбор рациональной структуры ТхС (схемная надежность);
•оценка надежности системы по ее структуре и количествен ным показателям составляющих элементов;
•определение методов технического обслуживания и ремон та системы с минимизацией экономических затрат.
Надежность ТхС зависит от свойств и характеристик самой системы и ее элементов. Необходимо учитывать характер отка зов, которые могут быть как внезапными, так и постепенными, связанными с процессами старения, коробления, износа и др.
Элементы системы могут быть восстанавливаемыми и невосстанавливаемыми. Невосстанавливаемый элемент (НЭ) — это такой элемент, свойства которого после отказа восстановить не удается или его восстановление экономически нецелесообразно. НЭ после отказа удаляется из системы и заменяется новым.
Вмашиностроении к невосстанавливаемым элементам относят типовые стандартные элементы, изготавливаемые в массовом производстве (подшипники, детали крепежа, агрегаты пневмо гидроаппаратуры, пружины и т.п.).
Восстанавливаемый элемент (ВЭ) после отказа может быть отремонтирован (восстановлен). К категории ВЭ обычно относят сложные дорогостоящие детали и сборочные единицы (узлы), изготовленные в условиях мелкосерийного или единичного про изводства и обладающие свойством восстанавливаемости. На пример, достаточно легко восстановить покрытие элемента или целого изделия. К ВЭ относятся сборочные единицы, которые рассматриваются при расчете надежности как отдельный, еди ничный элемент. Отказы из-за выхода из строя одной или не скольких деталей у подобных сборочных единиц устраняются за счет их замены при ремонте. Работоспособность многих дета лей удается восстановить технологическими методами (наплав кой, плазменным напылением, сваркой и др.).
На НЭ и ВЭ действуют внезапные, постепенные и комбиниро ванные отказы (рис. 3.6). Комбинированные отказы проявляются при действии субъективных факторов и под влиянием окружаю щей среды. К комбинированным отказам можно отнести пробой электроизоляции, вызванный ее старением, воздействием тем пературы, влажности и др.
При расчете показателей надежности технических систем их следует рассматривать как некоторую эквивалентную структуру
98 |
3. Математические основы теории надежности |
Рис. 3.6. Классификация отказов элементов технических систем
взаимосвязанных элементов. В свою очередь, системы могут быть невосстанавливаемые (НС) и восстанавливаемые (ВС) (рис. 3.7). К НС прежде всего следует отнести системы однократного при менения (ракеты, средства пожаротушения и т.д.). Учитывая, что на восстановление изделий затрачиваются значительные ре сурсы, а общий уровень надежности непрерывно повышается, удельная составляющая НС и НЭ постепенно возрастает. Так, многие электронные изделия, модули и блоки в настоящее время имеют такой высокий уровень надежности, что их восстановление практически нецелесообразно, поскольку большинство изделий имеют срок службы, сравнимый с периодом их морального ста рения. Поэтому общая тенденция при создании технических си стем состоит в достижении высокого уровня надежности как составляющих элементов, так и систем в целом.
Рис. 3.7. Классификация систем обеспечения надежности
3.2. Классификация технических систем и элементов |
99 |
С другой стороны, изделия и системы могут иногда оказы ваться под воздействием субъективных факторов, влиянием ок ружающей среды (температуры, влажности, тумана, гололеда
ит.д.), в аварийных ситуациях (автомобильный транспорт), ко гда доля незначительных повреждений и отказов весьма высока (например, проколы камер, разрыв покрышек). Поэтому в по добных условиях процессы восстановления объективно необхо димы.
ТС весьма дорогие, сроки их эксплуатации достаточно дли тельны, поэтому они в подавляющем большинстве случаев отно сятся к классу восстанавливаемых. Восстанавливаемые элементы
исистемы могут находиться в трех состояниях: работоспособном; состоянии ожидания восстановления; состоянии восстановления.
Наряду со способностью к восстановлению важное значение с позиций обеспечения надежности имеют структура и свойства системы, определяющие влияние отказа каждого элемента на надежность системы в целом. При этом для расчета показателей надежности любая конструкция или изделие должны быть пред ставлены в виде элементарных блоков, соединенных между собой функциональными связями. Элементарным блоком может быть модуль, сборочная единица, узел и даже отдельная деталь.
Классификация структур систем обеспечения надежности изделий достаточно сложная (см. рис. 3.7). Системы могут быть нерезервируемые и резервируемые. Нерезервируемые системы
спозиций надежности могут быть только последовательными. Последовательные системы (ПС) имеют структуру, которая пред ставляет собой последовательную цепочку элементов (1-4), каж дый из которых, за исключением первого 1 и последнего 4> соеди нен одной функциональной связью с предыдущим и последующим элементами (рис. 3.8). Отказ хотя бы одного элемента ПС приво дит к отказу всей системы в целом. К таким системам относятся большинство механических, электромеханических, оптических, оптико-механических и других изделий, у которых усложнение структуры с целью повышения надежности приводит к резкому увеличению стоимости изделия и, соответственно, к пропорцио-
1 2 3 4
Рис. 3.8. Структура системы с последовательным соединением элементов
100 |
3. Математические основы теории надежности |
нальному снижению эффективности его применения по служеб ному назначению.
Следует отметить, что схема расчета надежности может не сов падать со схемой построения (структурой) изделия. На рис. 3.9, а приведена параллельно-последовательная блок-схема диагности ческого вихретокового прибора для контроля колебаний шпин делей и роторов в подшипниках. Однако, поскольку отказ любого блока прибора приведет к отказу прибора в целом, используют последовательную схему расчета надежности (рис. 3.9, б).
Рис. 3.9. Параллельно-последовательная блок-схема (а)
и последовательная схема обеспечения расчета надежности (б) диагностического вихретокового прибора для измерения колебаний шпинделей и роторов в подш ипниках:
1 — частотный генератор; 2 — амплитудный детектор; 3 — усилитель вы сокой частоты; 4 — пиковый вольтметр; 5 — блок настройки; 6 — распре делительная схема; 7 — вихретоковые датчики; 8 — цепи экранирования; 9 — блок питания; 10 — частотный генератор; 11 — блок согласования
ПС может использоваться и при обеспечении надежности по методу резервирования для уменьшения вероятности ложного срабатывания, когда включение участка цепи происходит при