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

книги / Надежность и диагностика технологических систем

..pdf
Скачиваний:
9
Добавлен:
19.11.2023
Размер:
48.1 Mб
Скачать

3.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 — блок согласования

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