Учебная работа. Модель надежности
Модель надежности
Оглавление
Введение
Обзор литературы
ключевые понятия теории массового обслуживания
Процессы восстановления
Марковский процесс
Полумарковские процессы
Стратегия управления и её свойства
Резюме
Основные понятия теории надёжности
Надёжность восстанавливаемого элемента
Надёжность системы
Условие задачи
Обзор решения задачи
исследовательская часть
Постановка новой задачи
Решение в общем виде
Численный пример
Заключение
список литературы
Введение
Теории надёжности, как самостоятельной науке, уже более 50 лет, за которые она проделала большой путь. Первые работы в рамках данной науки содержали методологические понятия, математические модели создания высоконадёжных технических систем и задачи математической статистики. Данные научные исследования давали или системы, состоящей из различных изделий. однако в этих трудах отсутствовало описание рационального поведения технического персонала, не было практических инструкций для повышения надёжности объектов на практике.
Необходимость в подобного рода исследованиях заставила учёных обратить внимание на данное направление. Стали появляться научные работы, в которых авторы изучали оптимальную стратегию для профилактики отдельно работающего изделия. Следующим этапом стали исследования по обслуживанию сложных многоэлементных систем. Данная область в силу своей специфики является широкой и в то же время высоко востребованной. Моя работа в свою очередь также заключена в рамках этой проблемы.
Начав исследование с разбора базовых понятий и обзора накопленного материала по тематике, я приведу свою собственную задачу, а также результаты, которые мне удалось получить.
теория надежность массовое обслуживание
Обзор литературы
Ключевые понятия теории массового обслуживания
Процессы восстановления
Говоря о теории массового обслуживания, в первую очередь стоит обратить внимание на процессы восстановления. Процесс восстановления является частным случаем случайного потока однородных событий ξ (t), для которого область значений есть целые неотрицательные числа, и все траектории являются неубывающими ступенчатыми функциями. Чаще всего такой процесс задаётся с помощью совместного распределения случайной последовательности {tn= tn (ω), 1≤n<∞}, которая определяет моменты скачков, а номер n определяет номер скачка процесса ξ (t). Рекуррентным потоком с запаздыванием называется процесс восстановления, у которого функция распределения времени первого скачка отличается от функции распределения остальных скачков. Ступенчатый случайный процесс ξ (t) называется потоком с ограниченным последействием, если случайные величины, задающие периоды времени между скачками, взаимно независимы. Простой процесс восстановления можно определить как последовательность независимых положительных одинаково распределённых величин, задаваемых функцией распределения интервала времени между скачками F (t).
Функцией восстановления называют математическое ожидание функции ξ (t)
Стационарным процессом восстановления называется процесс восстановления, у которого распределение числа восстановлений, произошедших на любом интервале заданной длины, не зависит от расположения этого интервала.
Альтернирующим процессом восстановления называется процесс восстановления, у которого четные и нечётные скачки распределены по двум разным функциям распределения.» [2]
Пуассоновский процесс
простейшим (пуассоновским) потоком ξ (t) называется простой процесс восстановления, у которого интервалы между моментами восстановления распределены по экспоненциальному закону
с параметром .» [2]
Пуассоновский процесс обладает свойств:
·Стационарность
·Отсутствие последействия — марковское свойство;
·Ординарность — равенство нулю вероятности нулевой длительности интервала.
Марковский процесс
Случайный процесс ξ (t) с непрерывным временем и дискретным множеством состояний E (конечным или счётным) называется марковским, если для любого целого n>0, любого набора моментов 0
Полумарковские процессы
Полумарковский процесс ξ (t) определяется однородной двумерной марковской цепью или однородным процессом марковского восстановления:
Однородная марковская цепь определяется переходными вероятностями, которые будем в дальнейшем называть полумарковским ядром
где n≥1, и некоторым начальным распределением
кроме этого полагаем
Из неотрицательности величин следует, что при .
Для полумарковского процесса верно следующее свойство: при известном настоящем для марковского процесса с дискретным временем прошлое и будущее независимы.» [2]
Свойства полумарковского ядра:
·;
·неубывающие по x, непрерывные слева функции;
·
·
Стратегия управления и её свойства
Стратегией управления называется правило выбора моментов принятия решений и правило выбора решений из заданного множества допустимых решений (управлений).
Если решения принимаются в некоторые заранее определённые или случайные моменты времени, то такое управление (стратегия) называется дискретным управлением.
Если стратегия зависит только от истории процесса, то такая стратегия называется неупреждающей.
Стратегия управления называется марковской, если при принятии решений используется только информация о настоящем состоянии процесса.
Стратегия управления называется однородной (стационарной), если правило принятия решений не зависит от календарного времени.
Стратегия управления называется рандомизированной, если в процесс принятия решений введён случайный эксперимент, который определяется некоторой вероятностной мерой, построенной на измеримом пространстве допустимых управлений.
Стратегия управления называется детерминированной (нерандомизированной), если вероятностная мера, построенная на измеримом пространстве допустимых управлений, является вырожденной. [2]
Резюме
В данном разделе кратко изложены базовые определения и понятия теории массового обслуживания. Я счёл это необходимым для более глубокого понимания последующего материала, так как многие инструменты теории надёжности опираются на результаты, полученные в рамках этой науки. Обозначив ключевые моменты теории массового обслуживания, перейдём непосредственно к самой теории надёжности.
Основные понятия теории надёжности
В течение эксплуатации изделия внутри структуры происходят необратимые процессы, оказывающие, в том числе, влияние на свойства надёжности объекта. Под надёжностью авторы подразумевают склонность изделия к сохранению этих свойств. Иными словами, это процесс изменения качества системы во времени. Что же именно представляет из себя теория надёжности. В книге [1] приведено следующее определение: Теория надёжности — это новая научная дисциплина, изучающая общие закономерности, которых следует придерживаться при проектировании, испытаниях, изготовлении, приемке и эксплуатации изделий для получения максимальной эффективности от их использования» [1].
Первым базисным понятием в данной науке является безотказность. Под безотказностью авторы понимают способность изделия сохранять работоспособность (т.е. не иметь отказов) в течение определённого интервала времени в определённых условиях эксплуатации [1]. Отказ, в свою очередь, представляет из себя нарушение в изделии, ведущее к полной или частичной потере работоспособности. Для наглядности авторы приводят лампы накаливания в качестве изделия, максимально удобного для работы в рамках теории надёжности. Отказы подразделяются на постепенные, в результате постепенного износа изделия, и внезапные, связанные с качеством изделия. Так же авторы приводят третий вид отказов — самоустраняющиеся отказы, которые характерны преимущественно для вычислительных машин.
Другим ключевым термином теории надёжности является долговечность. Под долговечностью изделия понимают его способность к длительной эксплуатации при необходимом техническом обслуживании, в которое могут входить и различные виды ремонтов [1]. По окончанию срока, определяющего долговечность, изделие либо невозможно восстановить, либо не выгодно экономически. Измеряться долговечность может как временем работы объекта, так и числом циклов его работы или объемом произведённой работы. Авторы обращают внимание, что для некоторых систем или объектов понятия долговечности и безотказности могут быть синонимичны, однако в общем случае это две разные характеристики надёжности.
Третий показатель надёжности связан с изделиями, которые в течение своего обслуживания подвергаются ремонту и называется ремонтопригодностью. Ремонтопригодность характеризует такие аспекты как время обнаружения отказа, время и стоимость ремонтных работ и прочее. Классификация ремонтно-восстановительных работ приведена в следующей Tаблице 1.
Таблица 1.
Глубина восстановления системыНазвания восстановительных работВосстановительная работа проводится с работоспособной системойВосстановительная работа проводится с неработоспособной системойНикакого обновления в системе не проводитсяПлановый (внеплановый) осмотр или проверка работоспособности- Проводится полное обновлениеПлановая (внеплановая) предупредительная профилактикаПлановый (внеплановый) аварийно-профилактический ремонтПроводится обновление части системыПлановая (внеплановая) предупредительная профилактика части системыПлановый (внеплановый) аварийно-профилактический ремонт части системы
Таким образом, теория надёжности, как наука, базируется на трех вышеупомянутых терминах, совокупность которых и определяет то, что мы называем надёжностью изделия. В свою очередь, качество и надёжность изделия напрямую связаны с экономическими показателями изделия. Перед производителем практически всегда стоит вопрос, какое выбрать сочетание цены и качества изделия, для получения наилучших показателей экономической эффективности. По мнению авторов, необходимо разработать систему, связывающую факторы, влияющие на надёжность изделия, и его себестоимость. Приведу слова авторов о задачах теории надёжности как науки в их понимании. Теория надёжности устанавливает и изучает количественные характеристики (критерии) надёжности, исследует связь между показателями экономичности, эффективности и показателями надёжности;
разрабатывает методы проведения испытаний на надёжность и методы обработки и оценки результатов этих испытаний;
разрабатывает методы контроля надёжности, методы оптимальных режимов профилактических (регламентных) работ при эксплуатации изделий, методы обоснования норм запасных частей (элементов, деталей).
В теории надёжности разрабатываются методы установления режимов и выбора характеристик, обеспечивающих оптимальную надёжность, методы выбора оптимальных конструкций и схем, обеспечивающих заданную надёжность, оптимальные методы отыскания неисправностей в сложной аппаратуре и т.п.
При решении задач теория надёжности использует результаты исследований физических и химических процессов, лежащих в основе явлений, связанных с потерей качества [1].
Авторы приводят перечень основных математических методов, используемых в теории надежности. Это базовые методы таких дисциплин, как:
·теория вероятностей и математическая статистика
·теория информации
·теория массового обслуживания
·линейное и нелинейное программирование
·математическая ·статистическое моделирование
В завершение раздела авторы приводят различные примеры, показывающие что в различных ситуациях приоритетными оказываются разные характеристики надёжности, и решать каждую задачу нужно отталкиваясь от специфики эксплуатации изделия.
Надёжность элемента, работающего до первого отказа
Завершив раздел, связанный с основными понятиями, авторы [1] изучают характеристики надёжности элементов, работающих до первого отказа.
За t=0 обозначим момент начала работы изделия, а t=τ является моментом отказа, при этом τ будем называть временем жизни элемента. Время жизни изделия является случайной величиной, распределённой по следующему закону распределения
Величина Q (t) отражает вероятность отказа элемента до момента времени t. Авторы обращают внимание на то, что предположение о моделировании τ как случайной величины является серьезным допущением по ряду факторов, однако является лучшим на данный момент математическим средством.
Наряду с временем жизни часто употребляется другая функция
подразумевающая вероятность безотказной работы за время t. Иначе эту функцию называют функцией надёжности. Она представляет из себя монотонно убывающую функцию, равную 1 в нуле и имеющую предел
Примерный вид этой функции приведён ниже на рис.1.
Рис. 1.
далее авторы рассматривают такую характеристику надёжности как опасность отказа. фактически, она выражается следующим образом
Решив это уравнение относительно P (t), получим
.
накопленный опыт в данной области показывает, что обычно функция выглядит подобно тому, как это изображено на рис. 2.
Рис. 2.
Как мы видим на рисунке, график разделён на три области. Первая область именуется периодом приработки и большие значения функции связаны с тем, что на этом этапе происходит выявление бракованных изделий или других проблем, связанных с эксплуатацией.
После периода приработки следует область нормальной работы, а низкие значения функции объясняются тем, что бракованные изделия уже сломались, а остальные изделия имеют маленькую вероятность отказа в связи с небольшим износом.
Последний этап, называемый периодом старения, связан с постепенным возрастанием значений функции из-за необратимых процессов в изделии, связанных с износом.
Согласно этой структуре, авторы говорят о том, что первый период обычно не продолжительный и не берётся в рассмотрение, а срок службы изделия по документации заканчивается раньше периода старения. Исходя из этих соображений, можем принять, что
,
а, следовательно, функция надёжности выглядит как
Иными словами, мы пришли к экспоненциальному закону распределения с вероятностью отказа
и плотностью
По этому поводу авторы говорят следующее: Экспоненциальный закон очень популярен в теории надёжности. Эта популярность объясняется тем, что экспоненциальный закон физически очень естественный, прост и удобен для использования. почти все задачи, возникающие в теории надёжности для экспоненциальных законов распределения, оказываются на порядок проще чем для произвольных законов. почти все формулы в теории надёжности в случае экспоненциального закона резко упрощаются.
Основная причина этого состоит в том, что экспоненциальный законвероятность безотказной работы на данном интервале (t,t+ τ) не зависит от времени предшествующей работы t, а зависит только от длины интервала τ«. Из теории массового обслуживания это свойство мы знаем как марковское свойство. Математически оно выглядит следующем образом
Несмотря на участившуюся критику в сторону использования экспоненциального закона в связи с частыми случаями, когда его использование не оправдано контекстом задачи, оно остаётся самым распространённым и самым часто применяемым на практике.
Надёжность восстанавливаемого элемента
теперь предположим, что элемент имеет возможность восстанавливаться. характер восстановления может быть различным, как это было показано ранее в Таблице 1, однако в данный момент будем считать, что элемент заменяется на новый, а так же что время восстановления равно 0. В таком случае мы уже имеем не одну случайную величину τ, а набор независимых случайных величин (τ1, τ2,…,) распределённых по одинаковому закону распределения F (t):
Данная задача представляет из себя процесс восстановления из теории массового обслуживания. ключевым элементов для его изучения является считающий процесс ν (t), случайная величина, равная числу отказов, происшедших за время t. Закон распределения ν (t) выглядит следующим образом
Надёжность системы
После обзора, посвященного надёжности отдельно взятых изделий, логично рассмотреть систему из нескольких элементов, так как на практике обычно приходится иметь дело именно с различными сложными структурами. Будем предполагать, что элементы отказывают независимо друг от друга.
Самым важным и в то же время самым простым случаем является система из последовательного соединения элементов в рассмотрении до первого отказа. Под последовательным соединением подразумеваем, что при отказе одной из составляющих, отказывает вся система. В таком случае для безотказной работы всей системы до момента времени t необходима безотказная работа всех её частей на этом интервале. Тогда, если P (t) — функция надёжности системы, а p1 (t), p2 (t),…,pn (t) — функции надежности n отдельных элементов, из которых она состоит, то получаем
Это означает, что
откуда
Таким образом, если надёжности всех составляющих системы имеют экспоненциальное распределения с параметрами , то и надёжность самой системы имеет экспоненциальное распределение с параметром, равным сумме параметров подсистем.
Рассматривая параллельное соединение элементов, имеет смысл вместо функций надёжности рассматривать вероятности отказа Q (t)
где qi (t) — вероятность отказа i-го элемента. В этом случае надёжность системы уже не будет подчиняться экспоненциальному закону. Для случая равных элементов:
Задача по обслуживанию двух последовательно соединённых подсистем
Условие задачи
рассмотрим задачу, которую В.А. Каштанов и Е.Ю. Барзилович исследуют в книге [3]. Для начала, кратко приведу условие задачи:
Рассмотрим систему, состоящую из двух последовательно соединённых подсистем. (Отказ одной из подсистем означает отказ системы.) Считаем, что обе подсистемы равнонадёжны и имеют одинаковую функцию распределения времени безотказной работы F (x) и, кроме этого, что отказы в подсистемах появляются независимо и появившийся отказ обнаруживается мгновенно.
Предположим, что в исследуемой системе возможно проведение двух видов внеплановых восстановительных работ, которые начинаются только в моменты отказов системы:
·внепланового аварийно-профилактического ремонта, полностью обновляющего систему (обновляется и отказавшая подсистема и работоспособная);
·внепланового аварийно-профилактического ремонта части системы, когда обновляется только отказавшая подсистема.
таким образом, вмешиваться в работу системы и проводить восстановительные работы можно только после отказа системы, причем при восстановлении в обязательном порядке проводится обновление отказавшей подсистемы. Другое решение — решение об обновлении работоспособной части системы — с одной стороны, влечёт за собой дополнительные затраты времени или средств, с другой стороны, увеличивает время до появления последующего отказа системы, и, следовательно, это решение необходимо принимать таким образом, чтобы обеспечить оптимальные значения показателей, характеризующих качество функционирования системы. [3]
Исследователи ставят перед собой задачу нахождения оптимальной стратегии обслуживания данной сложной системы. При этом, авторы, отталкиваясь от общих соображений и теоремы о вырожденных распределениях, ищут оптимальную стратегию в классе вырожденных пороговых стратегий (без всякого обоснования). Под пороговой стратегией подразумевается следующее: если при отказе одной из подсистем наработка второй меньше некой τ, то с вероятностью единица принимается решение заменить только отказавшую подсистему, иначе так же с вероятностью достоверного события принимается решение о замене обеих подсистем.
Обзор решения задачи
Авторы решают данную задачу путём построения полумарковского процесса. В первую очередь, задаётся дифференциальный закон
Помимо того, определена функция распределения неотрицательной случайной величины для каждой пары :
Это позволяет построить ступенчатый полумарковский процесс , значения которого образуют обобщённую цепь Маркова, а длительности пребывания процесса в том или ином состоянии являются случайными величинами, распределение которых зависит от состояний, принимаемых процессом на рассматриваемом и последующем периодах. С каждой реализацией процесса можно связать некоторую случайную прибыль, определив её на каждом периоде как функцию длительности периода и состояния процесса на рассматриваемом и последующем периодах
[3]
Определяется средняя прибыль за период при условии, что процесс переходит из состояния х в у
а так же средняя Прибыль за промежуток времени, который находился в состоянии х
Определив среднюю длительность пребывания в состоянии x как
можем выразить среднюю удельную Прибыль, которая является главным показателем качества работы системы в данной задаче:
Исследовательская часть
Постановка новой задачи
Задачей моего исследования является решение задачи о двух последовательно соединённых подсистемах для дискретного времени и конечного множества состояний. Также передо мной стоит вопрос о проверке предположения об оптимальности пороговой стратегии.
Решение в общем виде
Дискретное время подразумевает, что поломка подсистем может происходить только в определённые моменты времени. Обозначим один квант времени как ∆. В таком случае, после начала работы системы, поломка может произойти через время равное ∆, 2∆, 3∆,…,N∆. Обозначаем момент времени N∆ как момент, к которому система с вероятностью единица откажет. Имеется в виду, что существуют вероятности отказа в моменты i∆, такие что
Данные вероятности зависят от характеристик подсистем и известны изначально. С помощью этих вероятностей можно построить ступенчатую функцию F (x), которая представляет из себя вероятность того, что подсистема откажет к моменту времени x
Кроме этого, для удобства введём функцию
,
которая равна вероятности противоположного события, то есть вероятность того, что подсистема не откажет к моменту времени x.
Будем рассматривать моменты времени отказа. Так как наработка может принимать только дискретные значения, в данной задаче будем иметь случайный процесс , определяемый как наработка неотказавшего элемента, который принимает значения из конечного множества состояний. особенностью задачи является то, что в отличии от непрерывного случая, возможна одновременная поломка обеих подсистем с вероятностью отличной от нуля. Это означает, что в множестве возможных состояний E, будет особое состояние, которое мы в дальнейшем будем обозначать как 0, когда отказали обе подсистемы. В остальных случаях отказывает одна из подсистем, а наработка второй выражается в квантах ∆, отлична от 0, но не превышает (N-1) ∆. Состояние, когда одна система отказала, а вторая имеет наработку i∆ будем обозначать как i. Тогда множество состояний будет выглядеть следующим образом
Решение о выборе ремонтных работ принимается во всех состояниях, за исключением 0, когда мы с вероятностью единица события обновляем обе подсистемы. Обозначим за решение u=1 ремонт одной подсистемы, а за u=2 — ремонт обеих составляющих. В таком случае в состояниях принимаются решения . Множество называется множеством допустимых управлений. Вероятностные меры на пространствах , то есть рандомизированные стратегии определяются заданием вероятностей
.
Будем предполагать, что при работе системы мы получаем Доход, а при ремонте убытки.
Известно [4], что таким образом построенный функционал — математическое ожидание дохода за единицу времени при длительном функционировании системы выражается равенством с использованием стационарных вероятностей
Средняя прибыль в состоянии i зависит от решения
где c0 — Прибыль, который приносит система за единицу времени в рабочем состоянии, с1 и c2 — затраты на единицу времени проведения ремонтных работ 1 и 2 соответственно, T1 и T2 — среднее время проведения ремонтных работ 1 и 2, а Mij — математическое ожидание времени работы системы до поломки при условии, что система находится в состоянии i и было принято решение j.
Средняя длительность периода в состоянии i выглядит следующим образом
кроме того, оптимальную стратегию можно искать в классе вырожденных, то есть либо либо .
Таким образом для каждой вырожденной стратегии необходимо найти и Mij. Для того чтобы найти математические ожидания, необходимо перебрать все варианты развития событий. В Таблице 2 приведены все возможные события для вычисления Mi2 (т.е. при принятии решения 2) и их вероятности. Так как при ремонте обеих составляющих система обновляется, то для всех состояний j Mi2 будут одинаковы.
Таблица 2.
P1P2…PN-1PNP1P1* P1P1*P2…P1*PN-1P1*PNP2P2*P1P2*P2…P2*PN-1P2*PN………………PN-1PN-1*P1PN-1*P2…PN-1*PN-1PN-1*PNPNPN*P1PN*P2…PN*PN-1PN*PN
Столбец слева — вероятности отказа одной из подсистем в момент Pi, верхняя строка — вероятности отказа второй подсистемы в момент Pi. Элементы таблицы представляют из себя одновременное выполнение двух условий. Сумма всех элементов с множителем P1 равна вероятности того, что система проработает один квант времени, после чего выйдет из строя. Сумма всех элементов с множителем P2, но не включающих множитель P1 равна вероятности, что система проработает 2 кванта времени, и т.д. На Рис.3 данная структура приведена более наглядно.
Рис. 3.
Таким образом, получив вероятности для каждого времени работы, мы находим математическое ожидание длительности работы для любого состояния при принятии решения 2.
В ситуации с принятием решения 1 такой однородности не будет, и для каждого состояния будет своя таблица. поэтому приведём таблицу в общем виде для состояния i. надо также помнить, что в связи с присутствием наработки у одной из подсистем, каждый элемент таблицы должен быть поделён на условную вероятность того, что он проработал без поломки i∆. Таблица 3 представляет из себя набор всех возможных событий при принятии решения 1 в состоянии i.
Таблица 3.
P1P2…PN-1PNPi+1Pi+1*P1Pi+1*P2…Pi+1*PN-1Pi+1*PNPi+2Pi+2*P1Pi+2*P2…Pi+2*PN-1Pi+2*PN………………PN-1PN-1*P1PN-1*P2…PN-1*PN-1PN-1*PNPNPN*P1PN*P2…PN*PN-1PN*PN
В данном случае структура суммирования остаётся такой же, однако так как это уже не квадратная, а прямоугольная таблица, то максимальная длительность будет не N∆, а (N-i) ∆. То есть при расчёте математического ожидания будут фигурировать не N различных длительностей и их вероятности, а N-i.
Посчитав все Mij, по приведённым ранее формулам вычисляются . Как уже говорилось раньше, состояние 0 является особым и решение в нём не принимается. Характеристики для нулевого состояния:
К настоящему моменту мы имеем среднюю прибыль и среднюю длительности периода для нулевого состояния, а также эти же характеристики, зависящие от решения для всех остальных состояний.
Теперь вычислим стационарные распределения . Они находятся из системы линейных уравнений
где
— элемент полумарковского ядра, вероятность перехода из состояния k в состояние i при решении u до момента времени t.
Записав все элементы полумарковского ядра, мы сможем вычислить среднюю удельную Прибыль. Заметим, что для решения задачи нам достаточно записать , что позволит существенно сократить громоздкость формул.
В первую очередь запишем элементы вида :
для осуществления этого перехода обе подсистемы должны отказать одновременно. В остальных случаях одна из подсистем должна выйти из строя на шаге j, а вторая нет.
Теперь осталось рассмотреть элементы вида . здесь мы будем иметь дело с элементами с наработкой, а, следовательно, будут присутствовать условные вероятности.
Обе подсистемы выходят из строя в один и тот же момент, с учётом наработки.
далее рассмотрим 4 случая для переходов из i в j. На данном этапе нужно обговорить, что переходы в этом случае для некоторых i и j могут осуществляться двумя способами: либо выходит из строя обновлённая система, либо поломка происходит в старой системе.
При этом переход, когда ломается новая составляющая возможен для тех случаев, когда i
Таким образом, из-за этих двух условий мы получаем четыре различные области, для которых элементы полумарковского ядра будут разные и записывать их мы будем отдельно.
Мы записали все элементы полумарковского ядра для всех возможных решений. Задача решена в общем случае, вся работа для нахождения удельной средней прибыли проделана. Теперь, используя полученные результаты, решим задачу для конкретного примера и попробуем доказать пороговость оптимальной стратегии.
Численный пример
Будем решать задачу для N=5. Остальные условия приведены в Таблице 4.
Таблица 4.
Распределение вероятностей отказаP1P2P3P4P50,20,150,150,20,3Удельные Прибыль и затратыСреднее время ремонтовс0с1с2T1T2159101,52
При выборе вероятностей отказа Pi я руководствовался общими соображениями и накопленным опытом, о котором я писал в обзоре литературы. Элемент P1 имеет завышенное P4 и P5 связаны с износом изделия, поэтому тоже имеют значения выше.
Затраты на первый и второй ремонт практически идентичны, так как основная их часть — это плата ремонтной бригаде, и только малая часть идёт на покупку необходимых комплектующих. акцент сделан на отличии среднего времени ремонтных работ, которое оказывает влияние и на итоговую стоимость ремонта.
Множество состояний E будет состоять из 5 элементов:
Чтобы вычислить математические ожидания времени пребывания нужно выписать таблицы, аналогичные тем, которые я описывал ранее.
При решении 2 для всех состояний будет одна таблица (Таблица 5).
Таблица 5.
P1P2P3P4P5P10,040,030,030,040,06P20,030,02250,02250,030,045P30,030,02250,02250,030,045P40,040,030,030,040,06P50,060,0450,0450,060,09∆2∆3∆4∆5∆0,360,21750,17250,160,09
В других таблицах приведу ситуации для разных состояний при решении 1:
Таблица 6.
P1P2P3P4P5P20,030,02250,02250,030,045P30,030,02250,02250,030,045P40,040,030,030,040,06P50,060,0450,0450,060,09∆2∆3∆4∆0,280, 1950,1750,15
Таблица 7.
P1P2P3P4P5P30,030,02250,02250,030,045P40,040,030,030,040,06P50,060,0450,0450,060,09∆2∆3∆0,250, 2050, 195
Таблица 8.
P1P2P3P4P5P40,040,030,030,040,06P50,060,0450,0450,060,09∆2∆0,260,24
Таблица 9.
P1P2P3P4P5P50,060,0450,0450,060,09∆0,3
далее с помощью полумарковского ядра считаем переходные вероятности при решениях 1 и 2 (Таблица 10). Выбор нужной для расчётов строки зависит от выбранной стратегии.
Таблица 10.
P (x,y); u=10123400000010,1781250,150,2843750,218750,1687520,16153850,184615380,20,38461540,069230830,170,320,3900,1240,20,8000P (x,y); u=20123400,2150,320, 1950,150,1210,2150,320, 1950,150,1220,2150,320, 1950,150,1230,2150,320, 1950,150,1240,2150,320, 1950,150,12
После этого задача доведена до вычисления средней удельной прибыли.
После того, как задача реализована в Excel, проверим пороговость оптимальной стратегии. Так как решение мы ищем в классе вырожденных стратегий, переберём все возможные вырожденные комбинации. Всего их 24=16 штук. результат данных вычислений представлен в Таблице 11.
Таблица 11.
Номер стратегииРешение в состоянии 1Решение в состоянии 2Решение в состоянии 3Решение в состоянии 4Средняя удельная прибыльЯвляется ли стратегия пороговой? 111224, 196да211124,150да311214,076нет412224,065да511114,023да612124,008нет712213,942нет812113,881нет921223,773нет1021123,709нет1122223,643да1221213,578нет1322123,576нет1421113,509нет1522213,449нет1622113,378нет
Из данной таблицы мы видим, что оптимальная стратегия оказалась пороговой, что подтверждает предположение о пороговости оптимальной стратегии. Так же можем заметить, что из 5 лучших стратегий 4 являются пороговыми, что также подтверждает правильность соображений авторов [3] по этому поводу.
Заключение
В данном исследовании мною в общем виде была решена задача об обслуживании двух последовательно соединённых подсистем с дискретной моделью времени и конечным множеством состояний. Помимо этого, данное решение было реализовано в Excel и оттестировано на конкретном примере. Наконец, на численном примере мною было подтверждено предположение о пороговости оптимальной стратегии, на которое опирались В.А. Каштанов и Е.Ю. Барзилович в процессе решения аналогичной задачи для непрерывного времени.
Список литературы
1. Гнеденко Б.В., Беляев Ю.К., Соловьев А.Д., Математические методы в теории надежности: Основные характеристики надежности и их статистический анализ. Изд.2-е, испр. и доп. М.: Книжный дом «ЛИБРОКОМ», 2017.
. Каштанов В.А. Элементы теории случайных процессов.
. Барзилович Е.Ю., Каштанов В.А. Некоторые математические вопросы теории массового обслуживания., М., Изд-во советское радио, 1971.
. Вопросы математической теории надежности (под ред. Б.В. Гнеденко) М., Изд-во Радио и связь, 1983.