Аннотация и ключевые слова
Аннотация (русский):
Апостериорный анализ надёжности выполняется после изготовления опытной партии аппаратуры с целью определения характеристик ее надёжности. Такие испытания необходимы потому, что на стадии проектирования устройства конструктор не располагает полными априорными сведениями, которые позволили бы заранее определить показатели надёжности с достаточно высокой достоверностью. Важным источником сбора информации о надёжности является система сбора данных о работе изделий в процессе их эксплуатации. Существуют два основных вида испытаний на надёжность. Один из них - определительные испытания, задачей которых является оценка показателей надёжности. Объект таких испытаний - крупносерийные изделия. Другой вид испытаний - контрольные испытания, в ходе которых осуществляется проверка соответствия показателя надёжности системы техническим условиям. Исследуется первый вид испытаний. Показан порядок проведения статистических испытаний радиоэлектронных систем по различным процедурам. Определены оценки среднего времени безотказной работы, достоверности среднего времени безотказной работы, длительности испытаний. Показано, что при оценке среднего времени безотказной работы метод максимального правдоподобия эффективнее испытаний по процедуре [ n , Б , r ].

Ключевые слова:
радиоэлектронная система, анализ надежности, показатели надежности, время безотказной работы, длительность испытаний
Текст
Введение Развитие радиоэлектронной промышленности приводит к быстрому росту функциональности выпускаемых изделий и усложнению структуры радиоэлектронных систем (РЭС) при одновременном повышении требований к их надёжности. Используемые модели имеют ряд недостатков, главным из которых является то, что они позволяют получить точную оценку показателей безотказной работы только в отдельных (частных) случаях [1-4]. Такая оценка пригодна для подтверждения требований технического задания, но не даёт возможности провести анализ надёжности РЭС после изготовления опытной партии аппаратуры. Именно поэтому определение характеристик надёжности изготовленных образцов РЭС является актуальной задачей. Постановка задачи и решения 1. Задачи апостериорного анализа Апостериорный анализ надёжности выполняется после изготовления опытной партии аппаратуры с целью определения характеристик ее надёжности. Для этого проводятся статистические испытания РЭС по одной из нижеперечисленных процедур [5]: а) процедура [n, Б, r] предполагает, что в испытаниях участвует n РЭС до отказов без замены отказавших систем; б) процедура [n, B, r] предполагает, что в испытаниях участвует n РЭС до отказов с заменой отказавших систем (восстановление); в) процедура [n, Б, T] предполагает, что в испытаниях участвует n РЭС в течение заданного времени (длительность испытаний) без замены отказавших систем; г) процедура [n, B, Т] предполагает, что в испытаниях участвует n РЭС в течение заданного времени T с заменой отказавших систем (восстановление); д) смешанные процедуры: [n, Б, r/T] или [n, B, r/T] предполагают, что задана длительность испытаний и число отказов. Испытания прекращаются, когда либо r, либо T достигают заданного значения; при этом, если длительность испытаний до последнего отказа , то обработка результатов выполняется по процедурам [n, Б, r] или [n, B, r], если , то обработка результатов выполняется по процедурам [n, Б, T] или [n, B, Т]; е) процедура [n, Б, r] - испытания проводятся до отказа всех n РЭС, участвующих в испытаниях; эта процедура используется редко, в основном в тех случаях, когда необходимо определить статистические характеристики последовательности отказов отдельных элементов РЭС. Каждая из процедур имеет определённые достоинства и недостатки, некоторые из них будут показаны в последующих разделах. Обработка результатов испытаний имеет целью решение одной из двух задач: 1-я задача. Определение характеристик надёжности изготовленных образцов РЭС. 2-я задача. Определение степени соответствия характеристик надёжности изготовленных образцов РЭС техническим условиям (решение этой задачи выходит за рамки нашего исследования). 2. Определение характеристик надёжности по результатам испытаний (1-я задача) 2.1. Оценка среднего времени безотказной работы (процедуры [n, Б, r]). 1. В результате испытаний по процедуре [n, Б, r] без замены отказавших систем получена выборка моментов времени отказов системы , по которой определена выборка интервалов времени между отказами (рис. 1). Рис. 1. Моменты отказов и интервалы между отказами 2. Сделана оценка полученной выборки на соответствие теоретической модели потока отказов (например, путём проверки соответствия функции распределения интервалов между отказами показательному закону распределения по критерию Колмогорова или критерию ). Решение задачи 1. Так как моменты времени t1, t2, t3, …, ti образуют последовательность зависимых случайных величин, в дальнейшем будем рассматривать интервалы времени между отказами y1, y2, …, yi, которые являются взаимно независимыми случайными величинами с показательной функцией плотности вероятностей (ФПВ): w (yi) = l (n - i + 1) exp{-l (n - i + 1) yi}, (1) где - интенсивность отказов одной системы; t*- среднее время безотказной работы одной системы; (n - i + 1) - число систем, участвующих в испытаниях с учётом отказавших. Совместная плотность вероятностей r интервалов yi с учётом их независимости равна wr (yi) = w r(y1, y2, …, yr) = . 2. Для определения оценки среднего времени безотказной работы используем метод максимального правдоподобия (МП). Обозначим функцию правдоподобия как Ly(t*) = . (2) В условиях решаемой задачи функция правдоподобия представляет собой ФПВ интервалов y при данном значении параметра t*. Сущность метода МП заключается в том, что в процессе обработки статистических данных вычисляется функция правдоподобия, а значение оценки искомого параметра (- оценка параметра t*) равно значению аргумента, при котором функция правдоподобия максимальна (рис. 2). Рис. 2. Возможный вид функций правдоподобия Тогда оценка среднего времени безотказной работы будет равна = = arg max Ly(t*); t*Î (0; t). 3. Найдем . Для поиска экстремума функции правдоподобия необходимо решить уравнение . Так как любая монотонная функция от функции правдоподобия также является функцией правдоподобия, то для упрощения решения перейдём к уравнению . С учётом (1) и (2) получим: Заменим l = 1/t* и продифференцируем: . (3) Решая (3), получим оценку среднего времени безотказной работы. . (4) 4. Преобразуем выражение (4) к моментам времени t1, t2, …, tr. yi = ti - ti-1; . Тогда , (5) где (n - r)tr - суммарное время работы не отказавших систем; - суммарное время безотказной работы всех отказавших систем; - суммарное время работы всех систем до r-го отказа. 2.2. Достоверность оценки среднего времени безотказной работы Оценка полученная по экспериментальным данным, является случайной величиной, точность приближения которой к истинному значению параметра t* зависит как от объёма испытаний, так и от методики испытаний и обработки результатов. В связи с этим требуется проверка полученной оценки на достоверность. Проверка осуществляется по четырем позициям: 1. Состоятельность оценки. 2. Несмещённость оценки. 3. Эффективность оценки. 4. Достаточность оценки. Проверка на состоятельность требует выполнения неравенства или , где e - сколь угодно малая величина. Вероятность того, что отклонение оценки от истинного значения равно нулю при увеличении числа отказов. В нашем случае оценка состоятельна в силу закона больших чисел. Для несмещённости оценки необходимо, чтобы математическое ожидание оценки было равно истинному значению параметра, т. е. для любого значения r. Для полученной оценки Следовательно, оценка является несмещенной. Для проверки эффективности оценки необходимо убедиться, что удовлетворяется неравенство Рао - Крамера: , (6) где - дисперсия оценки; Jy(t*) - информация по Фишеру, которая содержится в статистике y о параметре t*. . В рассматриваемом случае Ly(t*) - функция правдоподобия (2). Оценка будет эффективной, если неравенство Рао - Крамера (6) становится равенством: Найдем дисперсию полученной оценки: где известно из (5), дисперсия интервалов yi равна , т. к. yi - интервалы между отказами - имеют показательный закон распределения (1) который определяется одним параметром: l(n - i + 1). Найдем информацию по Фишеру - Jy(t*). Из (3) известно, что . Тогда Следовательно, и полученная оценка эффективна. Для достаточности оценки необходимо выполнение равенства (7) где - информация по Фишеру, которая содержится в статистике y о параметре . В то же время известно, что любая эффективная оценка достаточна. Таким образом, полученная нами оценка является состоятельной, несмещенной, достаточной и эффективной. Вместо оценки можно было искать оценку интенсивности отказов , такую как . Однако является смещённой и неэффективной оценкой. Математическое ожидание оценки не совпадает с измеряемым параметром, т. к. . Несмещённая оценка получится, если взять , но дисперсия несмещённой оценки оказывается равной и оценка остаётся неэффективной. 2.3. Доверительный интервал среднего времени безотказной работы Полученная в разделе 2.2 оценка среднего времени безотказной работы является точечной оценкой искомого параметра . Эта оценка является случайной величиной, которая в конкретном испытании может принять любое положительное значение - от 0 до µ, поэтому в дополнение к точечной оценке обычно определяется интервальная оценка измеряемого параметра. Имеется в виду, что по одной оценке определяется доверительный интервал , в котором находится истинное значение измеряемого параметра с заданной доверительной вероятностью , где g - доверительная вероятность (или коэффициент доверия); и - соответственно нижняя и верхняя границы доверительного интервала (рис. 3). Рис. 3. Доверительный интервал Для определения доверительного интервала надо знать функцию распределения вероятностей оценок . Для этого преобразуем выражение (7) для того, чтобы в нём использовались нормированные величины: при этом - длина интервала. Тогда (7) перепишется в виде . (8) С учётом того, что , или , выражение (8) примет вид , или . (9) Таким образом, необходимо найти ФПВ величины , т. е. w(). Из (5) можно определить, что суммарное время наработки на отказ равно . (10) Функция плотности вероятностей интервалов yi известна (1): , . В этом законе необходимо заменить переменную, чтобы получить стандартное распределение вероятностей с дисперсией равной 1. Обозначим: ; . (11) Тогда - экспоненциальное распределение с единичной дисперсией. Известно, что в этом случае имеет распределение Гаусса, а имеет распределение c2(2r) c 2r степенями свободы, которое широко используется в статистике для обработки экспериментальных данных. C учётом (10) и (11) . Введём переменную , t имеет распределение c2(2r) c 2r степенями свободы; - число отказов. Распределения c2(2r) табулированы. Для большого числа степеней свободы это распределение стремится к нормальному. Пусть и , тогда выражение (9) для доверительного интервала преобразуется к виду . (12) На рис. 4 показана ФПВ c2, заштрихованная площадь под кривой представляет собой доверительную вероятность g (вся площадь под кривой ФПВ, как известно, равна единице). Как видно из рис. 4, доверительный интервал можно расположить на оси t по-разному, т. е. решение задачи неоднозначно. Рис. 4. Функция плотности вероятностей c2 Обычно поступают так, чтобы границы интервала отсекали справа и слева одинаковые площади под кривой, равные . Тогда нижняя граница доверительного интервала точка распределения, а верхняя граница доверительного интервала точка распределения, значения которых определяются по таблицам квантилей распределения. Затем, на основании (12), , или , или . (13) Из (13) очевидно, что нижняя граница доверительного интервала равна , (14) а верхняя граница, соответственно, равна . (15) Таким образом, процедура определения верхней и нижней границ доверительного интервала сводится к следующему: 1. Перед испытанием задаётся - число отказов и - доверительная вероятность. 2. В процессе испытаний определяется (суммарное время наработки на отказ ). 3. По формулам (14), (15) вычисляются и 2.4. Длительность испытаний Длительность испытаний T = tr, т. е. совпадает с моментом r-го отказа, когда испытания прекращаются. Для процедуры испытаний (n, Б, r) это случайная величина. В то же время знание сроков проведения испытаний является важным как для исполнителей, так и для руководителей испытаний. Функцию плотности вероятностей w (Т) найти трудно, т. к. , а величины разнородны (распределение зависит от (1)), поэтому определим только среднее значение времени испытаний и его дисперсию. Среднее время испытаний , (16) т. е. больше, чем оценка среднего времени безотказной работы . Запишем (16) в виде ряда где . Тогда при , а при . Обозначим , тогда . Известно, что при функция . Значит, если , то среднее время испытаний m(T) = . Если , то m(T) = . Дисперсия времени испытаний равна Обозначим: ; Тогда . Если , то . Это означает, что дисперсия времени испытаний уменьшается с увеличением r, но эта процедура не очень эффективна, т. к. дисперсия стремится не к нулю, а к некоторому постоянному числу (как будет видно в дальнейшем, более эффективной процедурой испытаний является процедура с заменой отказавших блоков). 2.5. Оценка среднего времени безотказной работы (процедуры [n, B, r]) Рассмотрим испытания системы на надёжность по процедуре , когда отказавшие системы заменяются или восстанавливаются. В остальном задача аналогична рассмотренной в разделе 2.1. В результате испытаний получена выборка моментов отказов и интервалов между отказами (рис. 5). Рис. 5. Моменты отказов и интервалы между отказами Модель отказов для одной системы , , , или , , где - среднее время наработки на отказ. Решение задачи 1. Теперь выборка интервалов однородна и каждый интервал подчиняется закону , где - общая интенсивность отказов систем, участвующих в испытаниях. 2. Для определения оценки среднего времени безотказной работы также используем метод МП. В условиях решаемой задачи функция правдоподобия представляет собой закон распределения интервалов y при данном значении параметра t*: 3. Оценка по максимуму правдоподобия определяется как параметр, соответствующий максимальному значению функции правдоподобия: . Тогда оценка , где - суммарная наработка на отказ, общая для обоих планов испытаний. Отсюда следует, что , и, значит, качество оценки такое же, как и при процедуре [n, Б, r], при одинаковых и r. 4. Полученная оценка является несмещенной, т. к. , и эффективной, т. к. . 5. Средняя продолжительность испытаний равна . Если , то , т. е. меньше, чем для процедуры [n, Б, r], когда она была равна (). Дисперсия продолжительности испытаний равна Если n = r, то и стягивается по n, чем больше участников испытаний, тем меньше . Следовательно, по сравнению с процедурой испытаний [n, Б, r] рассматриваемая процедура более эффективна. Выводы Рассмотрена процедура испытаний опытной партии РЭС и определены следующие характеристики надёжности: - оценка среднего времени безотказной работы; - достоверность оценки и доверительный интервал среднего времени безотказной работы. Показано, что при определении оценки среднего времени безотказной работы метод МП эффективнее процедуры испытаний [n, Б, r].
Список литературы

1. Жаднов В. В. Проектная оценка надёжности радиотехнических систем / В. В. Жаднов, С. Н. Полесский // Надёжность и качество: тр. Междунар. симпоз.: в 2 т. Т. 1 / под ред. Н. К. Юркова. Пенза: Изд-во Пенз. гос. ун-та, 2006. С. 24-29.

2. Жаднов В. В. Управление качеством при проектировании теплонагруженных радиоэлектронных средств / В. В. Жаднов, А. В. Сарафанов. М.: Солон-Пресс, 2004. 464 с.

3. Артюхова М. А. Метод учёта влияния системы менеджмента надёжности предприятия при расчётной оценке показателей надёжности электронных средств / М. А. Артюхова, В. В. Жаднов, С. Н. Полесский // Радiоелектронiка, iнформатика, управлiння. 2013. № 2. С. 48-53.

4. Надёжность ЭРИ: Справочник. М.: МО, 2006. 641 с.

5. Левин Б. Р. Теория надёжности радиотехнических систем / Б. Р. Левин. М.: Сов. радио, 1978. 264 с.