Инженерия / Надёжность инженерных систем
Верхняя граница интенсивности при нуле отказов
Даже если в испытании не случилось ни одного отказа, можно оценить верхнюю доверительную границу постоянной интенсивности отказов. Чем больше суммарная наработка, тем ниже граница.
Формула
На одинаковой наработке граница при доверительном уровне 95% выше, чем при 90%.
Обозначения
- $λ_U$
- односторонняя верхняя доверительная граница интенсивности, ч⁻¹ или другая обратная единица времени
- $α$
- уровень значимости, равный 1 минус доверительный уровень, —
- $T$
- суммарная наработка всех испытуемых объектов, ч или другая единица времени
Условия применения
Наблюдалось ровно ноль отказов, суммарная плановая экспозиция положительна. Времена до отказа независимы и имеют общую постоянную интенсивность либо отказы описываются однородным пуассоновским процессом с фиксированной длительностью наблюдения.
Ограничения
Если интенсивность меняется с возрастом, эта граница не оценивает будущую интенсивность без дополнительной модели. При одном или нескольких отказах нужна формула с учётом их числа. Произвольное продление испытания после промежуточных проверок требует отдельного расчёта покрытия.
Подробное объяснение
Если время до отказа экспоненциально, вероятность проработать время t без отказа равна exp(−λt). Для независимых образцов вероятности перемножаются, а их времена в показателе складываются. Поэтому вероятность не увидеть ни одного отказа за суммарную наработку T равна exp(−λT).
Если эта вероятность для проверяемой интенсивности меньше α, наблюдённый ноль слишком редок для односторонней проверки на уровне α. Граничное значение получается из exp(−λ_UT)=α. Логарифм α отрицателен, так что λ_U положительна.
Вероятность здесь относится к результатам повторных испытаний, а интенсивность в частотной постановке считается фиксированным неизвестным параметром. Число 95% означает уровень процедуры, а не вероятность того, что уже вычисленный интервал включает параметр.
Для экспоненциальной модели средняя наработка равна 1/λ. Верхняя граница λ_U поэтому даёт нижнюю границу для средней наработки. Если тест требует ограничить λ заданным числом, целевое время находят из неравенства T≥−lnα/λ_U.
Как пользоваться формулой
- Сложите длительности фактической работы всех образцов и проверьте, что за учтённое время отказов не было.
- Выберите односторонний доверительный уровень, например 0,95, и вычтите его из единицы для получения α.
- Возьмите натуральный логарифм α со знаком минус. Разделите полученное число на суммарную наработку.
- Укажите обратную единицу времени. Для нижней границы средней наработки в этой модели обратите результат.
Историческая справка
В задачах промышленной надёжности нулевое число отказов регулярно возникает, если ресурс компонентов намного больше доступной длительности испытаний. Справочник 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
Связанные формулы
Аналитика
Оценка Каплана — Мейера при правом цензурировании
Оценка Каплана — Мейера восстанавливает долю объектов без события к заданному времени, когда некоторые наблюдения закончились раньше события. Кривая падает только в моменты событий.
Аналитика
Оценка Нельсона — Аалена для накопленного риска
Оценка Нельсона — Аалена складывает доли событий среди объектов, остававшихся под наблюдением. Она даёт накопленный риск, который может быть больше единицы.
Аналитика
Формула Гринвуда для ошибки оценки Каплана — Мейера
Формула Гринвуда оценивает выборочную дисперсию кривой Каплана — Мейера в заданный момент. Квадратный корень из дисперсии даёт стандартную ошибку этой оценки.
Инженерия
Вероятность безотказной работы при постоянной интенсивности отказов
Вероятность безотказной работы при постоянной интенсивности отказов показывает, как получить расчетную величину из проверяемых исходных данных. Формула полезна для предварительного инженерного расчета, потому что сразу связывает результат с единицами измерения и областью применимости.
Инженерия
Интенсивность отказов в экспоненциальной модели надежности
Интенсивность отказов в экспоненциальной модели надежности показывает, как получить расчетную величину из проверяемых исходных данных. Формула полезна для предварительного инженерного расчета, потому что сразу связывает результат с единицами измерения и областью применимости.
Инженерия
Распределение Вейбулла: безотказность и изменение риска
Двухпараметрическая модель Вейбулла описывает срок службы с убывающей, постоянной или растущей интенсивностью отказов. Параметр формы задаёт изменение риска с возрастом, а масштаб определяет временную шкалу.