Инженерия / Надёжность инженерных систем

Верхняя граница интенсивности при нуле отказов

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

Опубликовано: Обновлено:

Формула

$$\lambda_U=\frac{-\ln\alpha}{T}$$
Схема Наработка снижает верхнюю границу
500100015002000250000.0010.0020.0030.0040.0050.006λ_U, 1/чT, ч95%90%

На одинаковой наработке граница при доверительном уровне 95% выше, чем при 90%.

Обозначения

$λ_U$
односторонняя верхняя доверительная граница интенсивности, ч⁻¹ или другая обратная единица времени
$α$
уровень значимости, равный 1 минус доверительный уровень, —
$T$
суммарная наработка всех испытуемых объектов, ч или другая единица времени

Условия применения

Наблюдалось ровно ноль отказов, суммарная плановая экспозиция положительна. Времена до отказа независимы и имеют общую постоянную интенсивность либо отказы описываются однородным пуассоновским процессом с фиксированной длительностью наблюдения.

Ограничения

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

Подробное объяснение

Если время до отказа экспоненциально, вероятность проработать время t без отказа равна exp(−λt). Для независимых образцов вероятности перемножаются, а их времена в показателе складываются. Поэтому вероятность не увидеть ни одного отказа за суммарную наработку T равна exp(−λT).

Если эта вероятность для проверяемой интенсивности меньше α, наблюдённый ноль слишком редок для односторонней проверки на уровне α. Граничное значение получается из exp(−λ_UT)=α. Логарифм α отрицателен, так что λ_U положительна.

Вероятность здесь относится к результатам повторных испытаний, а интенсивность в частотной постановке считается фиксированным неизвестным параметром. Число 95% означает уровень процедуры, а не вероятность того, что уже вычисленный интервал включает параметр.

Для экспоненциальной модели средняя наработка равна 1/λ. Верхняя граница λ_U поэтому даёт нижнюю границу для средней наработки. Если тест требует ограничить λ заданным числом, целевое время находят из неравенства T≥−lnα/λ_U.

Как пользоваться формулой

  1. Сложите длительности фактической работы всех образцов и проверьте, что за учтённое время отказов не было.
  2. Выберите односторонний доверительный уровень, например 0,95, и вычтите его из единицы для получения α.
  3. Возьмите натуральный логарифм α со знаком минус. Разделите полученное число на суммарную наработку.
  4. Укажите обратную единицу времени. Для нижней границы средней наработки в этой модели обратите результат.

Историческая справка

В задачах промышленной надёжности нулевое число отказов регулярно возникает, если ресурс компонентов намного больше доступной длительности испытаний. Справочник NIST/SEMATECH по инженерной статистике рассматривает этот случай в разделе о постоянной интенсивности и показывает переход от границы интенсивности к границе средней наработки.

Та же модель используется в описанном там планировании демонстрационных испытаний на заданную наработку. В общей записи для разных чисел отказов появляются квантили хи-квадрат. При нуле отказов граница получается напрямую из вероятности пустого пуассоновского счёта, поэтому вычисляется через один натуральный логарифм.

Историческая линия формулы

Эта граница получается обращением вероятности нуля событий в модели с постоянной интенсивностью. Справочник NIST/SEMATECH приводит её для испытаний без отказов вместе с эквивалентной односторонней границей средней наработки.

Пример

Десять одинаковых изделий работали по 100 ч без отказов. Предполагаем постоянную интенсивность и независимость, программа испытания задана заранее. Найдём одностороннюю 95%-ю верхнюю границу. Наработка T=10·100=1000 ч. При α=0,05 получаем −ln(0,05)=2,995732, поэтому λ_U=2,995732/1000≈0,002996 ч⁻¹. Обратное значение даёт нижнюю границу средней наработки примерно 333,8 ч. Подстановка в вероятность нулевого результата exp(−0,002995732·1000) даёт 0,05. Показатель безразмерен, поскольку интенсивность в ч⁻¹ умножается на часы.

Частая ошибка

Ноль отказов не даёт нулевой верхней границы. Складывайте наработку всех изделий, а не только длительность эксперимента по календарю. Для 95%-й односторонней границы берите α=0,05. Десятичный логарифм вместо натурального исказит границу в модели с экспоненциальной вероятностью.

Практика

Задачи с решением

Испытание при уровне 90%

Условие. При T=500 ч отказов нет. Найдите верхнюю границу с доверительным уровнем 90%.

Решение. Берём α=0,10, поэтому числитель ln 10≈2,302585. Делим на 500 ч.

Ответ. Около 0,004605 ч⁻¹.

Сколько часов нужно

Условие. Требуется при нуле отказов получить 95%-ю границу не выше 0,001 ч⁻¹. Какова минимальная целая наработка в часах?

Решение. Делим −ln(0,05) на 0,001 и получаем 2995,732… ч. Округляем вверх.

Ответ. 2996 ч.

Дополнительные источники

  • https://www.itl.nist.gov/div898/handbook/apr/section4/apr451.htm
  • https://www.itl.nist.gov/div898/handbook/apr/section3/apr311.htm
  • https://www.itl.nist.gov/div898/handbook/apr/section1/apr132.htm

Связанные формулы

Аналитика

Оценка Каплана — Мейера при правом цензурировании

$\widehat S(t)=\prod_{t_i\leq t}\left(1-\frac{d_i}{n_i}\right)$

Оценка Каплана — Мейера восстанавливает долю объектов без события к заданному времени, когда некоторые наблюдения закончились раньше события. Кривая падает только в моменты событий.

Аналитика

Оценка Нельсона — Аалена для накопленного риска

$\widehat H(t)=\sum_{t_i\leq t}\frac{d_i}{n_i}$

Оценка Нельсона — Аалена складывает доли событий среди объектов, остававшихся под наблюдением. Она даёт накопленный риск, который может быть больше единицы.

Аналитика

Формула Гринвуда для ошибки оценки Каплана — Мейера

$\widehat{\mathrm{Var}}\!\left[\widehat S(t)\right]=\widehat S(t)^2\sum_{t_i\leq t}\frac{d_i}{n_i(n_i-d_i)}$

Формула Гринвуда оценивает выборочную дисперсию кривой Каплана — Мейера в заданный момент. Квадратный корень из дисперсии даёт стандартную ошибку этой оценки.

Инженерия

Вероятность безотказной работы при постоянной интенсивности отказов

$R(t)=e^{-\lambda t}$

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

Инженерия

Интенсивность отказов в экспоненциальной модели надежности

$\lambda=\frac{r}{T}$

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

Инженерия

Распределение Вейбулла: безотказность и изменение риска

$R(t)=e^{-(t/\eta)^\beta},\qquad h(t)=\frac{\beta}{\eta}\left(\frac{t}{\eta}\right)^{\beta-1}$

Двухпараметрическая модель Вейбулла описывает срок службы с убывающей, постоянной или растущей интенсивностью отказов. Параметр формы задаёт изменение риска с возрастом, а масштаб определяет временную шкалу.