Научная статья на тему 'ПОЛУМАРКОВСКАЯ МОДЕЛЬ ПРОЦЕССА ФУНКЦИОНИРОВАНИЯ СРЕДСТВ ИЗМЕРЕНИЙ'

ПОЛУМАРКОВСКАЯ МОДЕЛЬ ПРОЦЕССА ФУНКЦИОНИРОВАНИЯ СРЕДСТВ ИЗМЕРЕНИЙ Текст научной статьи по специальности «Математика»

CC BY
109
28
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ПОЛУМАРКОВСКИЙ ПРОЦЕСС / СРЕДСТВА ИЗМЕРЕНИЙ / ПЕРИОДИЧНОСТЬ ПОВЕРКИ / ПРОЦЕСС ФУНКЦИОНИРОВАНИЯ

Аннотация научной статьи по математике, автор научной работы — Ефремов Валерий Александрович, Мищенко Владимир Ильич, Мищенко Илья Владимирович

Актуальность и цели. Для обоснования величины интервала между поверками и расчета показателей эффективности функционирования при использовании по назначению средств измерений (СИ) используются различные модели функционирования в ходе использования по назначению. На основе предложенной полумарковской модели процесса функционирования средств измерений военного назначения обосновано определение периодичности их поверки. Показана актуальность применения теории полумарковских процессов для моделирования функционирования СИ. Материалы и методы. Разработана полумарковская модель процесса функционирования СИ, позволяющая комплексно учесть влияние на надежность средств измерений следующих факторов: периодичности поверки СИ, интенсивности использования по назначению, вероятности возникновения неисправности рабочего эталона, вероятности неисправности СИ, вероятности отказа ЗИП, продолжительности проведения поверки, продолжительности диагностирования СИ, продолжительности замены неисправных элементов, вероятностей возникновения ошибок I и II рода. Показателем эффективности функционирования СИ выбран коэффициент готовности средств измерений, являющийся аналогом коэффициента метрологической исправности для СИ. Результаты и выводы. С помощью данной модели можно рассчитать величину оптимального интервала между поверками СИ.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Ефремов Валерий Александрович, Мищенко Владимир Ильич, Мищенко Илья Владимирович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

SEMI - MARKOV MODEL OF THE FUNCTIONING PROCESS MEASURING INSTRUMENTS

Background. To substantiate the value of the interval between verifications and to calculate the performance indicators of functioning when measuring instruments (MI) are used for their intended purpose, various models of the functioning of the means during intended use are used. On the basis of the proposed semi-Markov model of the process of functioning of measuring instruments for military purposes, the determination of the frequency of their verification is justified. The relevance of the application of the theory of semi-Markov processes for modeling the functioning of MI is shown. Materials and methods. A semi-Markov model of the MI functioning process has been developed, which allows to comprehensively take into account the influence of the following factors on the reliability of measuring instruments: the frequency of MI verification, the intensity of intended use, the probability of a malfunction of the working standard, the probability of a malfunction of the MI, the probability of a failure of the spare parts, the duration of verification, the duration of MI diagnostics, the duration of replacement of faulty elements, the probabilities of errors I and II rhoda. The indicator of the effectiveness of the functioning of the MI is the readiness coefficient of measuring instruments, which is an analogue of the coefficient of metrological serviceability for the MI. Results and conclusions . With the help of this model, it is possible to calculate the value of the optimal interval between the verifications of MI.

Текст научной работы на тему «ПОЛУМАРКОВСКАЯ МОДЕЛЬ ПРОЦЕССА ФУНКЦИОНИРОВАНИЯ СРЕДСТВ ИЗМЕРЕНИЙ»

УДК 389.14

doi:10.21685/2307-5538-2022-3-3

ПОЛУМАРКОВСКАЯ МОДЕЛЬ ПРОЦЕССА ФУНКЦИОНИРОВАНИЯ СРЕДСТВ ИЗМЕРЕНИЙ

В. А. Ефремов1, В. И. Мищенко2, И. В. Мищенко3

1 2' 3 Военно-космическая академия имени А. Ф. Можайского, Санкт-Петербург, Россия

1 2' 3 vka@mil.ru

Аннотация. Актуальность и цели. Для обоснования величины интервала между поверками и расчета показателей эффективности функционирования при использовании по назначению средств измерений (СИ) используются различные модели функционирования в ходе использования по назначению. На основе предложенной полумарковской модели процесса функционирования средств измерений военного назначения обосновано определение периодичности их поверки. Показана актуальность применения теории полумарковских процессов для моделирования функционирования СИ. Материалы и методы. Разработана полумарковская модель процесса функционирования СИ, позволяющая комплексно учесть влияние на надежность средств измерений следующих факторов: периодичности поверки СИ, интенсивности использования по назначению, вероятности возникновения неисправности рабочего эталона, вероятности неисправности СИ, вероятности отказа ЗИП, продолжительности проведения поверки, продолжительности диагностирования СИ, продолжительности замены неисправных элементов, вероятностей возникновения ошибок I и II рода. Показателем эффективности функционирования СИ выбран коэффициент готовности средств измерений, являющийся аналогом коэффициента метрологической исправности для СИ. Результаты и выводы. С помощью данной модели можно рассчитать величину оптимального интервала между поверками СИ.

Ключевые слова: полумарковский процесс, средства измерений, периодичность поверки, процесс функционирования

Для цитирования: Ефремов В. А., Мищенко В. И., Мищенко И. В. Полумарковская модель процесса функционирования средств измерений // Измерения. Мониторинг. Управление. Контроль. 2022. № 3. С. 29-38. doi:10.21685/2307-5538-2022-3-3

SEMI - MARKOV MODEL OF THE FUNCTIONING PROCESS MEASURING INSTRUMENTS

V.A. Efremov1, V.I. Mishchenko2, I.V. Mishchenko3

1 2 3 Military Space Academy named after A.F. Mozhaisky, St. Petersburg, Russia

1 2 3 vka@mil.ru

Abstract. Background. To substantiate the value of the interval between verifications and to calculate the performance indicators of functioning when measuring instruments (MI) are used for their intended purpose, various models of the functioning of the means during intended use are used. On the basis of the proposed semi-Markov model of the process of functioning of measuring instruments for military purposes, the determination of the frequency of their verification is justified. The relevance of the application of the theory of semi-Markov processes for modeling the functioning of MI is shown. Materials and methods. A semi-Markov model of the MI functioning process has been developed, which allows to comprehensively take into account the influence of the following factors on the reliability of measuring instruments: the frequency of MI verification, the intensity of intended use, the probability of a malfunction of the working standard, the probability of a malfunction of the MI, the probability of a failure of the spare parts, the duration of verification, the duration of MI diagnostics, the duration of replacement of faulty elements, the probabilities of errors I and II rhoda. The indicator of the effectiveness of the functioning of the MI is the readiness coefficient of measuring instruments, which is an analogue of the coefficient of metrological serviceability for the MI. Results and conclusions. With the help of this model, it is possible to calculate the value of the optimal interval between the verifications of MI.

Keywords: semi-Markov process, measuring instruments, verification periodicity, functioning process

For citation: Efremov V.A., Mishchenko V.I., Mishchenko I.V. Semi-Markov model of the functioning process measuring instruments. Izmereniya. Monitoring. Upravlenie. Kontrol' = Measurements. Monitoring. Management. Control. 2022;(3):29-38. (In Russ.). doi:10.21685/2307-5538-2022-3-3

© Ефремов В. А., Мищенко В. И., Мищенко И. В., 2022. Контент доступен по лицензии Creative Commons Attribution 4.0 License / This work is licensed under a Creative Commons Attribution 4.0 License.

Интервалы между поверками (ИМП) являются одной из важных характеристик СИ, поскольку они определяют не только их технические характеристики, но и существенно влияют на боеготовность вооружения. Поэтому требования к ИМП постоянно ужесточаются, что находит свое отражение в технических требованиях на вновь разрабатываемые и модернизируемые СИ. На ранних этапах проектирования для установления ИМП используются расчетные методы, приведенные в РМГ 74-2004. В начале рассчитывается первоначальное значение интервалов между поверками на этапе утверждения типа СИ согласно РМГ 74-2004. В этом документе приведен ряд методик, основанных на использовании показателей метрологической надежности СИ, таких как вероятность того, что в течение заданной наработки метрологический отказ не возникнет или средняя наработка на метрологический отказ, равная математическому ожиданию наработки до первого метрологического отказа. Согласно РМГ 74-2004 аналогом коэффициента метрологической исправности СИ является комплексный показатель надежности СИ - коэффициент готовности, т.е. Кми ~Кг.

Наступление метрологического отказа может быть обнаружено только при поверке СИ, результаты которой позволят утверждать, что отказ произошел в период времени между двумя последними поверками. Является очевидным, что существует оптимальная величина ИМП.

Действительно, частые поверки приводят к отрыву СИ от использования их по назначению, подрывают боеготовность вооружения и военной технтки (ВВТ). Последнее обусловлено тем, что СИ, часто находясь на поверке, не могут обеспечить контроль и оценку технического состояния ВВТ. Это приводит к снижению коэффициента готовности и коэффициента метрологической исправности.

Редкие поверки приводят к тому, что в СИ накапливаются метрологические отказы, которые могут быть выявлены только при их поверке. И чем реже поверки, тем этих отказов больше. Накопление метрологических отказов приводит к тому, что после поверки СИ отправляются на ремонт. Это опять же приводит к снижению Кми и Кг. Таким образом, существует оптимальное значение ИМП, обеспечивающее максимум Кми и Кг и адаптацию значений этих показателей к достигнутому уровню безотказности и метрологической исправности СИ.

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

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

Граф этой модели представлен на рис. 1.

X ---------«

----- --

Рис. 1. Граф состояний модели процесса функционирования простых СИ: X и | - интенсивности отказов и восстановлений простых СИ соответственно

Решения системы уравнений Колмогорова для такой модели в стационарном режиме, представленной графом на рис. 1, на этапе нормальной эксплуатации имеют вид

|+Х |+Х

Так как X = —, a u = —, то P = -

T' ТВ 1

1

1

To

(

T

1 T3

1 1

—+—

T Tв

Л

T ^+1

T.

T.+Tв

= Kr, где T - средняя

наработка на отказ; ТВ - среднее время восстановления.

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

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

Рис. 2. Граф состояний для модели процесса эксплуатации СИ с избыточностью

Множество состояний графа Я = {Я}, * = 1,..., г (г - количество состояний графа, для графа на рис. 2, очевидно, г = 3) включает: Я - состояние исправности; Я2- состояние, когда на исправном СИ проводится поверка; Я3 - состояние, когда на неисправном СИ проводится

поверка, отыскивается неисправность и после этого средство восстанавливается.

В соответствии с работой [1] полумарковский процесс задается начальным состоянием и двумя матрицами - матрицей переходов вложенной марковской цепи (' = {м>у},(у) = 1,..., г) и матрицей условных функций распределения продолжительности пребывания в состояниях ¥(?Ш?) = {F¡J(?)},(/,]) = 1,..., г).

Для графа, представленного на рис. 2, эти матрицы имеют вид

0

W = 1 1

1 - Pa 0 0

СИ 0 ,

0 ^(t) *1з(0

F(t ) = F21(t ) F^(t )

0 0

В матрице W

Pc

1 си '1 - ехр(-^ситимп x

где РСИ - вероятность возникновения неисправности СИ; А,си - интенсивность возникновения метрологических неисправностей; Тимп - интервал между поверками. Неравные нулю элементы матрицы Б(?) имеют вид

Fu(t ) = \

0, t < T

им

1, t > T

'; ^(t )=\

0, t < Tn

0,

t < 0,

-^XPÍ-^¿L ,0 < t < t ,

имп

1 - exp(-Àситимп )

1,

t > ти

i jn> T-* 10,t<тп + тВ, F21(t) ii, t. Tn; ^ = I1, t > Тп +Тв,

0

где ТП - математическое ожидание продолжительности поверки; ТВ ние продолжительности восстановления СИ.

- математическое ожида-

Стационарное распределение вероятностей состояний полумарковского процесса в соответствии с работой [1] определяется по формуле

Рт —

щ==г-1—, ^ ] =1 г, (1)

У Рт,

где р - стационарная вероятность пребывания вложенной марковской цепи в 7-м состоянии; т7 - математическое ожидание продолжительности пребывания СИ в 7-м состоянии.

Стационарное распределение вероятностей пребывания вложенной марковской цепи в состояниях Р = (р , Р2, Р3) отыскивается в результате решения системы уравнений, записанной в матричном виде и удовлетворяющее условию нормировки

Р = РЖ, У [=1 р = 1, 0 <Р} < 1, ] = 1,г. (2)

Решение этой системы уравнений имеет вид

1 1 - Р Р

Р = _ • Р = си . Р = -* си

Р 2; Р2 2 • Рз 2 .

Математические ожидания пребывания в состояниях определяются по традиционной формуле [1]

т, = П1 - Ъ (Г)] Я,7= 1Г, (3)

здесь Ъ7 (/) - безусловные функции распределения вероятностей ухода процесса из соответствующих состояний, отыскиваемые в соответствии с [1] по формуле

ъ с)=у 3=1 ъ (') ^ ,(7, ])=17,

где ) - элементы матрицы ); м> ^ - элементы матрицы W .

Тогда средние времена пребывания полумарковского процесса в состояниях Ъ ) определятся как

т1 = Т~[1 - еХР (-—сиТимп )] ; т2 = ТП • тз = ТП + ТВ.

—си

Таким образом, показателем эффективности функционирования СИ в процессе использования по назначению будет вероятность застать полумарковский процесс и, следовательно, объект исследования в исправном состоянии в произвольный момент времени, т.е.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

—- [1- ехр (-КитШп )]

П = П--. (4)

—[1 - ехр(--си7;мП )] + Ти + [1 - ехр(-—^)]Тв

—си

Представленное выражение является аналогом коэффициента технического использования, приведенного в ГОСТ 27.002-2015. Как видно из формулы (4), приведенная модель учитывает безотказность объекта исследования — , продолжительность поверки ТП, величину межповерочного интервала Тимп и продолжительность восстановления ТВ . Недостатком модели является отсутствие учета качества поверки, которое определяется достоверностью контроля технического состояния СИ, проводимого в ходе поверки.

Исследованиям процессов использования по назначению сложных технических систем посвящены работы профессоров Ю. К. Беляева, А. Д. Соловьева, Е. Ю. Барзиловича, В. А. Каштанова, А. В. Смагина, Л. И. Волкова, Е. И. Сычева. В. И. Мищенко и др.

Среди них по совокупности учитываемых факторов выделяются именно полумарковские модели, разработанные профессорами В. И. Мищенко, Е. И. Сычева и Л. И. Волкова [2-4].

Как показывает анализ моделей процесса функционирования технических систем [2, 3, 5-7], наиболее адекватной с точки зрения учета влияния достоверности контроля параметров

при поверке СИ является модель, предложенная профессором В. И. Мищенко [2]. Эта модель взята в качестве базовой. Ее суть заключается в том, что контроль пригодности СИ, проводимый в рамках поверки, сводится к измерению, как правило, некоторого постоянного набора независимых параметров. По результатам контроля принимается решение об отнесении состояния объекта к пригодному или непригодному. Обнаружение метрологических неисправностей, обусловивших непригодность, осуществляется только при проведении поверки. Периодичность поверки образца СИ определяется как продолжительность временного промежутка Тимп , по истечении которого СИ поверяется в течение случайного времени Тп, ограниченного сверху временем ТП, т.е. Тп < Тп. Своевременная поверка средств измерений обеспечивает требуемый уровень их пригодности к использованию по назначению. Чем выше уровень сложности СИ (соответственно, уровень избыточности), тем критичнее становятся принятые решения о годности СИ к достоверности полученных результатов. При этом большое значение приобретает метрологическая исправность рабочего эталона. Поэтому при разработке модели процесса функционирования СИ при использовании по назначению необходимо учитывать возможность снижения достоверности показаний СИ в случае метрологического отказа рабочего эталона при проведении поверки СИ. Таким образом, при использовании СИ по назначению могут возникнуть следующие ситуации:

1. Перед началом проведения поверки СИ с вероятностью 1 - Рси (Тимп) пригодно к использованию по назначению, рабочий эталон с вероятностью 1 - РЭ (ТП) пригоден к использованию по назначению. В результате проведения поверки СИ в течение ограниченного сверху времени Тп < Тп поверитель может принять следующие решения:

а) СИ пригодно (правильное решение) - с вероятностью 1 -а1. После проведения поверки СИ продолжает применяться по назначению до следующей поверки;

б) СИ непригодно (ошибочное решение, соответствующее ложной регистрации метрологической неисправности СИ, т.е. ошибке первого рода) с вероятностью ). Причем эта ошибка обусловлена в первую очередь ограниченной метрологической надежностью эталона. В этом случае проводится дополнительная поверка, позволяющая выявить достоверно отсутствие метрологического отказа. Проведение дополнительной поверки проводится в течение ограниченного сверху времени тдп < ТдП .

2. До проведения поверки СИ с вероятностью 1 - РСИ (Тимп) пригодно к использованию по назначению, рабочий эталон с вероятностью РЭ (ТП) отказал. В зависимости от исхода поверки, а также характера метрологического отказа рабочего эталона (частичного или полного) по результатам поверки СИ могут быть произведены следующие оценки его состояния:

а) СИ непригодно (ложный метрологический отказ) - с вероятностью Ъ2. В этой ситуации проводится дополнительная проверка;

б) СИ пригодно - с вероятностью 1 -а2.

3. До начала проведения поверки СИ с вероятностью РСИ (Тимп) оказалось непригодным, рабочий эталон с вероятностью 1 - РЭ (ТП) пригоден к использованию по назначению. В результате проведения поверки о состоянии СИ могут быть приняты следующие решения:

а) СИ непригодно (правильное решение) - с вероятностью 1 - в1. После обнаружения метрологического отказа СИ проходит дополнительную поверку в течение времени тдп < ТдП . Затем, в случае наличия отказавшего элемента в ЗИП, с вероятностью 1 - Рзип, проводится замена этого элемента на работоспособный за время < ТЗ (временем поиска в комплекте ЗИП пренебрегается). Если же требуемый элемент, с вероятностью Рзип, отсутствует и в комплекте ЗИП, тогда оформляется заявка на предприятие-изготовитель, которая удовлетворяется в течение случайного времени экстренной поставки tЭд < ТЭП и после замены и проведения поверки продолжает использоваться по назначению;

б) СИ пригоден (пропуск отказа) - с вероятностью в2. В этом случае СИ продолжает использоваться по назначению с метрологическим отказом, снижающим эффективность функционирования СИ при использовании по назначению, до следующей поверки.

4. К началу проведения поверки СИ с вероятностью Рси (Тимп) непригодно, рабочий эталон с вероятностью РЭ (ТП) неисправен. В зависимости от результатов поверки и характера отказа рабочего эталона производится следующая оценка состояния СИ:

а) СИ непригодно - с вероятностью 1 — в2, после чего проводится дополнительная проверка и замена отказавшего элемента с использованием комплектов ЗИП, описанным в п. 3а;

б) СИ пригодно (пропуск отказа) - с вероятностью в2, что приводит к ситуации, описанной в п. 3б.

Граф, соответствующий этой модели, представлен на рис. 3. Предложенная модель процесса функционирования СИ учитывает основные факторы, влияющие на это функционирование, и позволяет определять основные показатели эффективности функционирования СИ при использовании по назначению.

Рис. 3. Граф полумарковской модели процесса функционирования при использовании по назначению СИ

Множество возможных состояний Я = {1,2,... , г} включает 9 состояний (т.е. г = 9): Я1 - СИ пригодно и используется по назначению (исходное состояние полумарковского процесса); Я2 - СИ непригодно и функционирует с неисправностью, обусловившей непригодность, до начала поверки; Я3 - СИ поверяется при условии, что оно к началу проведения поверки было пригодно; Я4 - СИ поверяется при условии, что оно к началу проведения поверки было непригодно; Я5 - СИ проходит дополнительную проверку при условии ложной регистрации пригодности; Я6 - СИ используется по назначению с неисправностью, обусловившей заключение о непригодности, до очередной поверки; Я - СИ проходит дополнительную проверку, подтверждающую его неисправность, и восстанавливается; Я - СИ ожидает экстренной поставки элемента, отсутствующего в комплекте ЗИП; Я9 - СИ восстанавливается, производится замена отказавшего элемента из состава ЗИП, после чего производится повторная поверка, подтверждающая его исправность и используется по назначению.

Матрица ' переходных вероятностей вложенной марковской цепи имеет вид

' =

0 Рси (1 — Рси ) 0 0 0 0 0 0

0 0 0 1 0 0 0 0 0

— а 0 0 0 а 0 0 0 0

0 0 0 0 0 в 1 -в 0 0

1 0 0 0 0 0 0 0 0

0 0 0 1 0 0 0 0 0

0 0 0 0 0 0 0 Р зип (1 — Рзип )

0 0 0 0 0 0 0 0 1

1 0 0 0 0 0 0 0 0

a = al(l-P„) + a2P„; ß = F;(l-P„)+ßP

где a,ß - вероятности возникновения ошибок контроля I и II рода; P - вероятность возникновения метрологического отказа рабочего эталона.

Для нормального этапа процесса функционирования СИ интенсивности отказов элементов считаются константами и тогда вероятность возникновения метрологического отказа СИ в течение периода между поверками Pœ определится по традиционной формуле

= ¡СИ (тимр ) = 1 - еХР(-^СИТимп X

(5)

где РСИ (Тимп) - вероятность возникновения неисправности, обусловившей непригодность СИ, на интервале [0,Тимп]; XСИ - интенсивность метрологических отказов СИ.

В матрице W указана также вероятность Рзип , которая соответствует вероятности отказа ЗИП, т.е. вероятности того, что после обнаружения отказа и обращении к ЗИП требуемый элемент в ЗИП будет отсутствовать. Эта вероятность рассчитывается по традиционным методикам, приведенным в ГОСТ РВ 27.3.03-2005.

Ненулевые элементы матрицы Б(г) определяются с учетом особенностей обеспечения адекватности математического представления реального процесса и представлены ниже

Fl2(t ) = i

0,

1 - e'x~t

1 - e l,

t < 0,

-X „Т™ ,0 < t < Тимп,

1, t — Т,„п ;

Fu(t) = FM(t ) = J0, t < ^ F24(t) = J

I1, t — т имп ;

,-ХгиТи,

( -1)

1 - g-V™

l,

t < 0, 0 < t < Тим

t — T„„ ;

10, г < ТП, 10, г < Тдп,

ЗД = ^(0 = ЗД = ^(0 = г>тп; ад = ^ г>тда

[1, г > тП; [1, 1 > тДП;

Г0, г < т , Г0, г < тЭД, Г0, г < тЗАМ + ТП,

к8(0 = Р79(0<! д Р89(0 = ^ ЭД Р91(0^ ' ЗАМ П' 1 г > т . 1г > т . 1 г > т + т

В представленных формулах обозначены: тимп - периодичность поверки СИ; тП - продолжительность проведения поверки СИ; тдП - продолжительность повторной проверки в случае ошибки 1-го рода; т - продолжительность диагностирования СИ; тЭд - экстренная доставка исправного элемента; тЗАМ - продолжительность замены неисправных элементов.

В соответствии с формулами (1)-(3) были получены соответствующие параметры полумарковского процесса:

- стационарные вероятности вложенной марковской цепи - набор (Р1, Р2,..., Р9) :

P

l

1 2(1 + % + F (1 - PJ + P(1 - P3IJ

; P2

P.,,

D

2(1 + D) + F(1 - P ) + P (1 - Pw )'

Pз =-

1 - P.„

2(1 + D) + F(1 - PCи ) + Pси (1 - Pзип )

; P4 =

P.

D

2(1 + D) + F (1 - P. ) + P. (1 - P3un)

■■О"

Р = ^ — Рси) . Р Рси (1 — В) .

^5 - Р ' Р6 = "

2(1 + В) + F (1 — Рси ) + Рси (1 — Рзип ) В

2(1 + В) + F(1 — Рси ) + Рси (1 — Рзип )

Р Р Р

Р __1си_. р * --

тл ч ^ О

7 2(1 + В-) + F(1 —Рси) + Рси(1 — Рзип) 8 2(1 + В) + F(1 —Рси) + Рси (1 —Рзип)

Рси

Р =--,

9 2(1 + В) + F (1 —Рси) + Рси (1 — Рзип)'

- математические ожидания продолжительности пребывания в состояниях Fi (?) -набор (т^ т2,..., т9):

т, = Р ; т2 = -^(Т — т,); т3 = т4 = ТП;

1 л си' 2 V имп 1 /' 3 4 П'

Хси Рси

т5 = ТДП; тб = Тимп; т7 = ТД ; т8 = ТЭД; т9 = ТЗАМ.

Распределение вероятностей состояний полумарковского процесса - набор (п1, п2,..., п9):

1

Р

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

X си Р

= - си

1 (Тимп + Тп )) —Рси + Р:^ ^ + Тд^(1 — Рси )+ Рси (Тд + тзам + Тп + Р^д ) ^ ' п =ХсиТ,Мп — Рси , п = Тп (1 — Рси ) , п = РТп , п = ТДП(1 — Рси ) ,

2" Хсиб ' 3" д ' 4" вд> д '

п = Рси (1 — В)Тимп • п = РсиТД , _ = РсиРзипТЭД , _ = Рси (ТЗАМТП)

Пб = вд ' п" д ' п" д ' п= д '

Здесь п1 соответствует вероятности застать объект поверки, т.е. СИ в состоянии пригодности в произвольный момент времени, т.е. коэффициенту готовности. Но теперь этот коэффициент явно зависит от интенсивности метрологических отказов, обусловивших непригодность, величины межповерочного интервала, продолжительности поверки, ошибок первого и второго родов и продолжительностей диагностирования, замены и экстренной доставки.

Заключение

Представленная в данной статье полумарковская модель процесса функционирования СИ позволяет комплексно учесть влияние на надежность средств измерений различных факторов: периодичности поверки СИ, интенсивности использования по назначению, вероятности возникновения неисправности рабочего эталона, вероятности неисправности СИ, вероятности отказа ЗИП, продолжительности проведения поверки, продолжительности диагностирования СИ, продолжительности замены неисправных элементов, вероятностей возникновения ошибок I и II рода. С помощью данной модели можно рассчитать величину оптимального интервала между поверками СИ.

Список литературы

1. Харламов Б. П. Непрерывные полумарковские процессы. СПб. : Наука, 2001. 431 с.

2. Мищенко В. И., Храмов М. Ю. Проблематика эксплуатации сложных технических систем. СПб. : Политехника-сервис, 2016. 218 с.

3. Сычев Е. И. Метрологическое обеспечение радиоэлектронной аппаратуры. М. : Татьянин день, 1993. 274 с.

4.

5.

6.

7.

Волков Л. И. Управление эксплуатацией летательных комплексов. М. : Высш. шк., 1981. 368 с. Новиков А. Н., Кравцов А. Н., Ширямов О. А. Математическая модель прогнозирования изменения значения критической составляющей погрешности рабочего эталона единицы величины с учетом априорной информации // Труды Военно-космической академии имени А. Ф. Можайского. 2018. Вып. 664. С. 102-109.

Новиков А. Н. Алгоритм индивидуального прогнозирования предельных экономически целесообразных сроков эксплуатации измерительных комплексов // Вестник СибГУТИ. 2016. № 4. С. 19-24. Мищенко В. И., Кравцов А. Н., Мамлеев Т. Ф. Полумарковская модель функционирования резервируемых средств измерений с учетом периодичности поверки // Измерительная техника. 2021. № 4. С. 22-27.

References

1. Kharlamov B.P. Nepreryvnye polumarkovskie protsessy = Continuous semi-Markov processes. Saint Petersburg: Nauka, 2001:431. (In Russ.)

2. Mishchenko V.I., Khramov M.Yu. Problematika ekspluatatsii slozhnykh tekhnicheskikh system = Problems of operation of complex technical systems. Saint Petersburg: Politekhnika-servis, 2016:218. (In Russ.)

3. Sychev E.I. Metrologicheskoe obespechenie radioelektronnoy apparatury = Metrological support of electronic equipment. Moscow: Tat'yanin den', 1993:274. (In Russ.)

4. Volkov L.I. Upravlenie ekspluatatsiey letatel'nykh kompleksov = Management of the operation of aircraft complexes. Moscow: Vyssh. shk., 1981:368. (In Russ.)

5. Novikov A.N., Kravtsov A.N., Shiryamov O.A. Mathematical model of forecasting changes in the value of the critical component of the error of the working standard of the unit of magnitude taking into account a priori information. Trudy Voenno-kosmicheskoy akademii imeni A. F. Mozhayskogo = Proceedings of the Military Space Academy named after A.F. Mozhaisky. 2018;(664):102-109. (In Russ.)

6. Novikov A.N. Algorithm of individual forecasting of marginal economically expedient-different service life of measuring complexes. Vestnik SibGUTI = Bulletin of SibSUTI. 2016;(4):19-24. (In Russ.)

7. Mishchenko V.I., Kravtsov A.N., Mamleev T.F. Semi-Markov model of functioning of redundant measuring instruments taking into account the frequency of verification. Izmeritel'naya tekhnika = Measuring equipment. 2021;(4):22-27. (In Russ.)

Информация об авторах /Information about the authors

Валерий Александрович Ефремов

адъюнкт,

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

(Россия, г. Санкт-Петербург, ул. Ждановская, 13) E-mail: vka@mil.ru

Valeriy A. Efremov

Adjunct,

sub-department of metrological support

of weapons, military and special equipment,

Мilitary Space Academy

named after A.F. Mozhaisky

(13 Zhdanovskaya street, St. Petersburg, Russia)

Владимир Ильич Мищенко

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

(Россия, г. Санкт-Петербург, ул. Ждановская, 13) E-mail: vka@mil.ru

Илья Владимирович Мищенко

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

(Россия, г. Санкт-Петербург, ул. Ждановская, 13) E-mail: vka@mil.ru

Vladimir I. Mishchenko

Doctor of technical sciences, professor,

professor of the sub-department of metrological

support of weapons, military and special equipment,

Мilitary Space Academy

named after A.F. Mozhaisky

(13 Zhdanovskaya street, St. Petersburg, Russia)

Il'ya V. Mishchenko

Candidate of technical sciences,

senior lecturer - head of the metrological service,

Мilitary Space Academy

named after A.F. Mozhaisky

(13 Zhdanovskaya street, St. Petersburg, Russia)

Измерение. Мониторинг. Управление. Контроль. 2022. № 3 ■о-...................................................................................................

Авторы заявляют об отсутствии конфликта интересов / The authors declare no conflicts of interests.

Поступила в редакцию/Received 05.05.2022 Поступила после рецензирования/Revised 06.06.2022 Принята к публикации/Accepted 11.07.2021

i Надоели баннеры? Вы всегда можете отключить рекламу.