Одним из важных показателей в технике является надежность, элементами которой являются: наработка до отказа, наработка между отказами. Согласно ГОСТ Р53480-2009 [1], отказ – это потеря способности изделия выполнять требуемую функцию; наработка до отказа – наработка, накопленная от первого использования изделия или от его восстановления до отказа; наработка между отказами – интервал времени между двумя последовательными отказами. На основе наработки между отказами формируется поток отказов, по параметрам которого можно судить о надежности.
При производстве штучных изделий с длительным сроком функционирования возникает проблема малой статистики. Это обусловлено тем, что для подобных изделий создание значительного числа опытных образцов для формирования статистики отказов становится нецелесообразным и экономически не выгодным [2]. В результате данные, формирующие потоки отказов таких изделий, поступают непосредственно после запуска в эксплуатацию, а количество данных получаемых на выходе очень мало и не позволяет обоснованно применить к ним вероятностные методы. Следовательно, имитация потока отказов требует специального подхода. В этом случае для восстановления распределения статистики можно воспользоваться гистограммным методом, рассматриваемым в рамках численного вероятностного анализа [3].
Гистограммой называется случайная величина, плотность распределения которой представлена кусочно-постоянной функцией. Гистограмма P определяется сеткой {xi | i = 0, …, n}, на каждом отрезке [xi−1, xi], i = 1, . . . , n гистограмма принимает постоянное значение pi.
Пусть требуется найти гистограмму pz суммы
и пусть p(x1, x2, …, xn) — плотность распределения вероятностей случайного вектора (x1, x2, …, xn). Тогда вероятность попадания z в интервал (zi , zi+1) соответственно равна
Реализация арифметических операций основана на работе с p(x, y) –совместной плотностью вероятности двух случайных величин x, y. Пусть pz – гистограмма, приближающая плотность вероятности, арифметической операции над двумя случайными величинами x ∗ y, где ∗ ∈ {+, −, ·, /, ↑}. Тогда вероятность попадания величины z в интервал [zi , zi+1] определяется по формуле
Для решения задачи восстановления распределения отказов по малой выборке используют имитационное (статистическое) моделирование потока отказов.
Функционирование оборудования характеризуется случайной величиной T, представляющей значение интервала времени между двумя последовательными отказами
Рассмотрим построение гистограммы PТ для величины Т. Пусть для величины Т известна выборка (Т1, Т2,…, ТN). Обозначим nj количество членов Ti повторной выборки, попавших в интервал [zi—1, zi], тогда
Таким образом, для каждой входной величины Ti может быть известно не только среднее значение Ti, но и гистограмма Pi. Такой подход можно интерпретировать как построение функции плотности вероятности некоторой случайной величины. Пр y в виде гистограммы, т.е. ставится задача: зная гистограммы PТi построить гистограмму Py, аппроксимирующую распределение выходной величины y. Обычно для этих целей используется метод Монте–Карло [4].
Проблематичностью применения метода Монте–Карло является длительность и ресурсоемкость проведения имитационного эксперимента. Для каждой реализации случайной величины используется генератор псевдослучайных величин, осуществляющих набор действий сначала по генерации равномерно-распределенной псевдослучайной величины, а затем по ее отображению на ту форму представления, которая принята за основу в конкретной модели. Чем сложнее представление, тем больше затрачивается процессорного времени и тем дольше идет имитационный эксперимент. Поэтому там, где можно свести просчет модели к разумному по сложности аналитическому преобразованию, это приветствуется. Конечно, важна и алгоритмическая осуществимость таких преобразований, которую система моделирования должна уметь осуществлять без участия человека [5].
Методика заключается в следующем:
- Все данные хi разбиваются на интервалы [ 0, xi], ( хi, xi+1]…, ( хj, xj+1].
- Строится гистограмма распределения Px.
- Строиться гистограмма накопленной вероятности распределения Px’.
- Гистограмма накопленной вероятности разбивается на интервалы [0, xi]’, ( хi, xi+1]’,…, ( хj, xj+1]’.
- С помощью генератора случайных чисел создается набор псевдослучайных величин mi∈[0, 1].
- Определяется к какому интервалу ( хi, xi+1] принадлежит mi..
- Строится гистограмма накопленной вероятности распределения псевдослучайной величины.
- Проводиться сравнительный анализ полученных гистограмм, с целью получить ошибку модели.
- Генерируется произвольное число значений, описывающих поток отказов, как в плане частоты сбоев, так и в плане их длительности.
С помощью данной методики, например, было осуществлено восстановление распределения сбоев в работе одного из узлов космического аппарата. Для этого на основе статистики были построены гистограммы плотности распределения частот сбоев (рисунок 1) и их длительностей (рисунок 2).
Рисунок 1. Гистограмма плотности распределения частот появления сбоев
Рисунок 2. Гистограмма плотности рапределения длиьельности сбоев
После построения этих гистограмм использя вычислительню среду MathCAD были определены уравнения сглаживающих кривых. С помощью которых было восстановлено распределение в первом приближениии.
Далее, с помощью генератора случайных чисел был получен ряд сулчайных значений di∈[0;1] [6]. После чего для каждой случайной величие di было найдено соответствующее значение xi и восстановлены гистограммы функций распределения часоты появления сбоев и их длительностей (рисунки 3 и 4 соответсвенно).
Результаты показывают, что восстановление функций распределения, опираясь на гистограммы, и их последующая имитация [5] действительно позволяет получить приближенную функцию распределения при малом объеме выборки, тогда как для аналитических методов требуется гораздо большее количество статистики [4]. Таким образом, удалось отработать методику восстановления потока отказов по малому объему статистики. Были разработаны основные алгоритмы методики, по которым в дальнейшем будет создаваться соответствующее программное обеспечение.
Список литературы:
- ГОСТ Р53480-2009 Надежность в технике. Термины и определения. – Москва: Стандартинформ, 2 — 26 с.
- Uglev V.A., Popova O.A., Dobronets B.S. The accuracy calculation control of reliability indices for equipment responsible appointment // International Siberian Conference on Control and Communications (SIBCON). – Omsk: OmGTU, 2015. – pp. 5-8. DOI: 10.1109/SIBCON.2015.7147248.
- Добронец Б.С., Попова О.А. Численный вероятностный анализ неопределенных данных – Красноярск: Сиб. федер. ун-т, 2014. – С. 167.
- Dobronets B., Popova O. Numerical probabilistic approach for data nonparametric analysis // Applied methods of statistical analysis. Nonparametric approach: Proceedings of the international workshop, – pp. 376-384.
- Углев В.А. Выбор между методом Монте-Карло и гистограммной арифметикой при реализации моделей с элементами случайности // Имитационное моделирование. Теория и практика: Материалы VI Всероссийской научно-практической конференции. В 2 т. Том 1. — Казань: ФЭН, 2013. – С. 278-281.
- Углев В.А., Устинов В.А. Имитационное моделирование: учеб. пособие // Абакан: Сиб. федер. ун-т., 2011. – 116 с.[schema type=»book» name=»МОДЕЛИРОВАНИЕ ПОТОКА ОТКАЗОВ ПО МАЛОМУ ОБЪЕМУ СТАТИСТИКИ» description=»Описание методики моделирования потока отказов по малым объемам статистики. Проведено моделирование потока отказов по малому объему статистики. Анализ полученных результатов.» author=»Копытов Иван Васильевич» publisher=»Басаранович Екатерина» pubdate=»2016-12-14″ edition=»euroasia-science_6(27)_23.06.2016″ ebook=»yes» ]