Научная статья на тему 'Программный комплекс имитационного моделирования и расчета стационарных вероятностей и оценки надежности резервированной системы с произвольными распределениями времени безотказной работы и ремонта её элементов'

Программный комплекс имитационного моделирования и расчета стационарных вероятностей и оценки надежности резервированной системы с произвольными распределениями времени безотказной работы и ремонта её элементов Текст научной статьи по специальности «Компьютерные и информационные науки»

CC BY
65
14
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ИМИТАЦИОННОЕ МОДЕЛИРОВАНИЕ / НАДЕЖНОСТЬ РЕЗЕРВИРОВАННЫХ СИСТЕМ / ОТНОСИТЕЛЬНАЯ СКОРОСТЬ ВОССТАНОВЛЕНИЯ / ВЕРОЯТНОСТЬ БЕЗОТКАЗНОЙ РАБОТЫ СИСТЕМЫ / ЧУВСТВИТЕЛЬНОСТЬ / ОЦЕНКА НАДЁЖНОСТИ / SIMULATION / RELIABILITY OF REDUNDANT SYSTEMS / RELATIVE RECOVERY RATE / PROBABILITY OF FAILURE-FREE SYSTEM OPERATION / SENSITIVITY / RELIABILITY ASSESSMENT

Аннотация научной статьи по компьютерным и информационным наукам, автор научной работы — Уанкпо Гектор Жибсон Кинманон, Козырев Дмитрий Владимирович

С быстрым развитием и распространением компьютерных сетей и информационных технологий, перед исследователями возникают новые сложные и комплексные задачи как прикладного, так и теоретического характера по исследованию надежности и доступности сетей и систем передачи данных. Одной из них является исследование надежности резервированных систем передачи данных между устройствами. Одним из основных условий стабильной работы системы передачи данных между устройствами является нечувствительность качества и эффективности системы к изменениям исходных параметров модели. Анализ такой чувствительности, который предполагается провести в рамках данной работы, является одним из новых направлений исследований сетей передачи данных последующих поколений. Мы рассматриваем имитационную модель восстанавливаемой системы передачи данных как модель замкнутой однородной системы холодного резервирования с одним ремонтным устройством с произвольной функцией распределения времени безотказной работы и произвольной функцией распределения времени ремонта её элементов. Для анализа и сравнения результатов были выбраны следующие распределения времени безотказной работы и времени ремонта элементов: Экспоненциальное (M), Вейбулла-Гнеденко (WB) и Логнормальное (LN). Ранее в [1] было показано, что явные аналитические выражения для стационарного распределения рассматриваемой системы удается получить не всегда. Разработанная в этой работе имитационная модель позволила исследовать надежность системы, определяемую как стационарную вероятность безотказной работы системы, а также получить оценки характеристик надёжности системы. Также получены значения коэффициента ρ (относительная скорость восстановления) при котором достигается заданный уровень надежности, построены графики зависимости вероятности безотказной работы системы и графики равномерной разности результатов имитационной модели для разных распределений в зависимости от относительной скорости восстановления. Программная реализации алгоритмов имитационного моделирования была осуществлена на основе языка R.

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

Похожие темы научных работ по компьютерным и информационным наукам , автор научной работы — Уанкпо Гектор Жибсон Кинманон, Козырев Дмитрий Владимирович

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

SOFTWARE TOOL FOR SIMULATION AND CALCULATING STATIONARY PROBABILITIES AND RELIABILITY OF A REDUNDANT SYSTEM WITH ARBITRARY DISTRIBUTIONS OF UPTIME AND REPAIR TIME OF ITS ELEMENTS

With the rapid development and proliferation of computer networks and information technologies, researchers are faced with new complex problems, both applied and theoretical, of studying the reliability and availability of networks and data transmission systems. One of them is the reliability study of redundant data transmission systems between devices. One of the main conditions for the stable functioning of the data transmission systems is the insensitivity of their quality and effectiveness to changes in the initial parameters of the model. The analysis of such sensitivity, which is supposed to be carried out as part of this work, is one of the new directions in the research of the next-generation data transmission networks. We consider the simulation model of the repairable data transmission system as a model of a closed homogeneous cold standby system with one repair unit with arbitrary distribution functions of uptime and repair time of its elements. In this paper, we study the system-level reliability using the simulation approach. Also, we obtained the values of the relative recovery speed at which the desired level of reliability is achieved, presented plots of the system failure probability and plots of the uniform difference of the obtained simulation results against the relative speed of recovery. The simulation algorithm was implemented based on the R language.

Текст научной работы на тему «Программный комплекс имитационного моделирования и расчета стационарных вероятностей и оценки надежности резервированной системы с произвольными распределениями времени безотказной работы и ремонта её элементов»

ТЕОРЕТИЧЕСКИЕ ВОПРОСЫ ИНФОРМАТИКИ, ПРИКЛАДНОЙ МАТЕМАТИКИ, КОМПЬЮТЕРНЫХ НАУК И КОГНИТИВНО-ИНФОРМАЦИОННЫХ ТЕХНОЛОГИЙ / THEORETICAL QUESTIONS OF COMPUTER SCIENCE, COMPUTATIONAL MATHEMATICS, COMPUTER SCIENCE AND COGNITIVE INFORMATION TECHNOLOGIES

УДК 519.237.5

DOI: 10.25559^ГПТО.15.201903.553-562

Программный комплекс имитационного моделирования и расчета стационарных вероятностей и оценки надежности резервированной системы с произвольными распределениями времени безотказной работы и ремонта её элементов

Г. Ж. К. Уанкпо1, Д. В. Козырев12

1 Российский университет дружбы народов, г. Москва, Россия 117198, Россия, г. Москва, ул. Миклухо-Маклая, д. 6

* [email protected]

2 Институт проблем управления им. В.А. Трапезникова РАН, г. Москва, Россия 117997, Россия, г. Москва, ул. Профсоюзная, д. 65

Аннотация

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

Мы рассматриваем имитационную модель восстанавливаемой системы передачи данных как модель замкнутой однородной системы холодного резервирования с одним ремонтным устройством с произвольной функцией распределения времени безотказной работы и произвольной функцией распределения времени ремонта её элементов. Для анализа и сравнения результатов были выбраны следующие распределения времени безотказной работы и времени ремонта элементов: Экспоненциальное (М), Вейбулла-Гнеденко ^В) и Логнормальное Ранее в [1] было показано, что явные аналитические выражения для стационарного распределения рассматриваемой системы удается получить не всегда. Разработанная в этой работе имитационная модель позволила исследовать надежность системы, определяемую как стационарную вероятность безотказной работы системы, а также получить оценки характеристик надёжности системы. Также получены значения коэффициента р (относительная скорость восстановления) при котором достигается заданный уровень надежности, построены графики зависимости вероятности безотказной работы системы и графики равномерной разности результатов имитационной модели для разных распределений в зависимости от относительной скорости восстановления. Программная реализации алгоритмов имитационного моделирования была осуществлена на основе языка R.

Ключевые слова: имитационное моделирование, надежность резервированных систем, относительная скорость восстановления, вероятность безотказной работы системы, чувствительность, оценка надёжности.

Контент доступен под лицензией Creative Commons Attribution 4.0 License. The content is available under Creative Commons Attribution 4.0 License.

Vol. 15, No. 3. 2019 ISSN 2411-1473 sitito.cs.msu.ru

Modern Information Technologies and IT-Education

Финансирование: публикация подготовлена при поддержке Программы повышения конкурентоспособности РУДН «5-100» (получатель Д. В. Козырев, разработка имитационной модели). Исследование выполнено при финансовой поддержке Российского фонда фундаментальных исследований в рамках научного проекта № 17-07-00142 а «Разработка моделей и методов анализа надежности соединений в гетерогенных сетях 5в для концепции Интернета надежных вещей» (получатель Д. В. Козырев, проведение численного анализа).

Для цитирования: Уанкпо Г. Ж. К., Козырев Д. В. Программный комплекс имитационного моделирования и расчета стационарных вероятностей и оценки надежности резервированной системы с произвольными распределениями времени безотказной работы и ремонта её элементов // Современные информационные технологии и ИТ-образование. 2019. Т. 15, № 3. С. 553562. DOI: 10.25559^ГПТО.15.201903.553-562

Современные информационные технологии и ИТ-образование

Том 15, № 3. 2019 ^ 2411-1473 sitito.cs.msu.ru

THEORETICAL QUESTIONS OF COMPUTER SCIENCE, COMPUTATIONAL MATHEMATICS, COMPUTER SCIENCE AND COGNITIVE INFORMATION TECHNOLOGIES

Software Tool for Simulation and Calculating Stationary Probabilities and Reliability of a Redundant System with Arbitrary Distributions of Uptime and Repair Time of its Elements

Hector G. K. Houankpo1, Dmitry V. Kozyrev12

1 Peoples' Friendship University of Russia, Moscow, Russia 6 Miklukho-Maklaya Str., Moscow 117198, Russia

* [email protected]

2 V.A. Trapeznikov Institute of Control Sciences of RAS, Moscow, Russia 65 Profsoyuznaya Str, Moscow 117997, Russia

Abstract

With the rapid development and proliferation of computer networks and information technologies, researchers are faced with new complex problems, both applied and theoretical, of studying the reliability and availability of networks and data transmission systems. One of them is the reliability study of redundant data transmission systems between devices. One of the main conditions for the stable functioning of the data transmission systems is the insensitivity of their quality and effectiveness to changes in the initial parameters of the model. The analysis of such sensitivity, which is supposed to be carried out as part of this work, is one of the new directions in the research of the next-generation data transmission networks.

We consider the simulation model of the repairable data transmission system as a model of a closed homogeneous cold standby system with one repair unit with arbitrary distribution functions of uptime and repair time of its elements.

In this paper, we study the system-level reliability using the simulation approach. Also, we obtained the values of the relative recovery speed at which the desired level of reliability is achieved, presented plots of the system failure probability and plots of the uniform difference of the obtained simulation results against the relative speed of recovery. The simulation algorithm was implemented based on the R language.

Keywords: Simulation, reliability of redundant systems, relative recovery rate, probability of failure-free system operation, sensitivity, reliability assessment.

Funding: This publication was prepared with the support of the RUDN University Competitiveness Enhancement Program 5-100 (recipient D. V. Kozyrev, development of a simulation model). The study was conducted with financial support from the Russian Foundation for Basic Research as part of scientific project No. 17-07-00142 a "Development of Models and Methods for Analyzing the Reliability of Connections in 5G Heterogeneous Networks for the Concept of the Internet of Reliable Things" (recipient D. V. Kozyrev, numerical analysis).

For citation: Houankpo H.G.K., Kozyrev D. V. Software tool for simulation and calculating stationary probabilities and reliability of a redundant system with arbitrary distributions of uptime and repair time of its elements. Sovremennye informacionnye tehnologii i IT-obrazovanie = Modern Information Technologies and IT-Education. 2019; 15(3):553-562. DOI: 10.25559/SITIT0.15.201903.553-562

Modern Information Technologies and IT-Education

Введение и мотивация

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

В последнее время функционирование различных аспектов современного общества стало критически зависеть от сетей связи [2,3]. С миграцией критически важных объектов в сети связи стало жизненно важным обеспечить высокую надежность и доступность сетей и систем передачи данных. Ряд предыдущих исследований [4-10] был сфокусирован на анализе надежности различных сложных телекоммуникационных систем. В частности, было проведено исследование надежности резервированных систем передачи данных в холодном режиме.

В [11] было рассмотрено создание имитационных моделей и средств для поддержки автоматизированного проектирования высоконадежных распределенных компьютерных систем. В [12] был представлен анализ надежности комбинированной электростанции, работающей на газотурбинном двигателе. Целью работы [13] была разработка модели для изучения надежности системы и анализа чувствительности доступности системы. В [14] был рассмотрен метод имитационного моделирования для моделирования надежности выполнения задания сложной системой посредством моделирования циклограммы задания, моделирования профиля среды и метода динамического моделирования надежности. В [15] были представлены методы моделирования и оценки, позволяющие проводить температурную оптимизацию надежности многопроцессорной системы на кристалле для конкретных приложений. В [16-25] также были рассмотрены различные подходы к моделированию и анализу надежности сложных систем.

В текущей работе обобщаются результаты предыдущих исследований авторов на случай холодного резервирования системы 1^01N / 01 /с произвольной функцией распределения (ФР) времени безотказной работы (в.б.р.) и произвольной ФР времени ремонта её элементов. Целью работы является проведение имитационного моделирования для нахождения значения коэффициента р (относительная скорость восстановления), при котором достигается заданный уровень надежности и построение графиков зависимости вероятности безотказной работы системы от относительной скорости восстановления; приводятся результаты вычисления оценки надёжности для разных исходных распределений.

Описание модели и постановка задачи

В качестве имитационной модели резервированной системы передачи данных, состоящей из п разнотипных каналов передачи данных, рассмотрим восстанавливаемую систему многократного холодного резервирования (0!„ / 01 /1^ с одним ремонтным устройством, с произвольным законом распределения в.б.р. её элементов и произвольным законом распределения времени их ремонта.

В данной работе будет рассмотрена зависимость вероятности безотказной работы системы (01п / 01 /1^ от относительной скорости восстановления. Ставится задача нахождения стационарной вероятности безотказной работы системы

для некоторых частных случаев распределении и оценка надежности системы, при n = 3.

Модель расчета стационарных вероятностей состояний системы

(GIn /GI/1)

Определим следующие состояния моделируемой системы:

• состояние 0: один (главный) элемент работает, (N-1)- в холодном резерве;

• состояние 1: один элемент отказал и находится в ремонте, один - работает, (N - 2) в холодном резерве;

• состояние 2: два элемента отказали, один находится в ремонте, второй ждёт своей очереди на ремонт, один - работает, (N -3)- в холодном резерве;

• состояние N -: все элементы отказали, один находится в ремонте, остальные ждут своей очередь на ремонт.

• Для описания алгоритма моделирования надёжности системы (OIN / GI /1) введём следующие переменные:

• double t - часы модельного времени, меняются при отказе или восстановлении элементов системы;

• int i, j - переменные состояния системы; при наступлении события осуществляется переход из i в j;

• double tnextfail - служебная переменная, в которой хранится время до следующего отказа элемента;

• double t „ - служебная переменная, в которой

nextrepair J г ' г

хранится время до следующего окончания ремонта отказавшего элемента;

• int k - счетчик числа итераций основного цикла.

Для наглядности имитационная модель представлена графически на рисунке 1 в виде блок-схемы. Критерием останова основного цикла модели является достижение максимального модельного времени выполнения T.

Исходные данные: 1 -число отказавших элементов N —число элементов в системе X —момент отказа элемента У —момент ремонта элемента tтeк = 0 -текущее время Т —максимальное модельное время

Р и с. 1. Блок-схема имитационной модели для оценки стационарных вероятностей

F i g. 1. Block diagram of a simulation model to assess the stationary probabilities

Современные информационные технологии и ИТ-образование

Том 15, № 3. 2019 ISSN 2411-1473 sitito.cs.msu.ru

THEORETICAL QUESTIONS OF COMPUTER SCIENCE, COMPUTATIONAL MATHEMATICS, COMPUTER SCIENCE AND COGNITIVE INFORMATION TECHNOLOGIES

Для лучшего понимания и воспроизводимости имитационной модели помимо блок-схемы также приводится алгоритм дискретно-событийного процесса имитационного моделирования в виде псевдокода с комментариями (Алгоритм 1).

else i = N ; t ., := ж; j := j -1; t := t

' nextfail ' J J ' m

end

if t > T then

Алгоритм 1. Псевдокод процесса имитационного моделирования системы / 01 Ввод: а1, Ь1, N Т, N0,"^".

а1 - Среднее время между отказами элементов, Ь1 - Среднее время ремонта, N - Число элементов в системе, Т - Максимальное модельное время прогона, N0 - Число графиков траекторий, "01" - Функция распределения.

Вывод: стационарные вероятности состояний

Р Р Р Р

1 О'111 2>'">1 М'

t = T

end

if i = v then

s, := rf _ GI (X) ; ^ := rf _ GI("S(x)") [ t •= t + s ■ t •= t + s ■

nextfail ' ; lnextrepair • 1 ^ J2 ;

if t ,, , < t „ then

nextfail nextrepan

j . j + 1, t . tnex,faii t else := j -1; t := t „ . ;

J ' nextreprnr

end

r[,k]:=[,i, j];i := j;k := к +1 [;

Begin

array r []:=[0,0,0j; // многомерный массив, содержащий результаты, к-го шага основного цикла

double t := 0.0; // инициализация часов модельного времени

int i := 0; j := 0; // переменные состояния системы

double tnexfail := 0.0; // переменная, в которой хранится время до следующего отказа элемента

double tnexllepail := 0.0; // переменная, в которой хранится время до следующего окончания ремонта

int k := 1; // счетчик числа итераций основного цикла

s := rf _ Gl (Xi); // генерация произвольной случайной величины s - время до первого события (отказа)

ss := rf _ GI ((х); // генерация произвольной случайной величины ss - время ремонта отказавшего элемента)

tnextfail t + s ;

t '= t + ss ■

Lnextrepah • L ^ 1515 ;

while t < ж do

if i = 0 then

^nextiepau ' ^; J ' J + ^ ' ^nextfail;

else for (vinl : (N -l))

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

end do

Оценка длительности пребывания в каждом состоянии i, i = 0,1,2, • • • N; вычисление стационарных вероятностей

| NG

Р =_V (длительность пребывания в состоянии i/T).

' NGji 1

end

Алгоритм 1 был программно реализован на языке R.

В Таблице 1 приведены значения коэффициента р = —

bi

относительной скорости восстановлении (отношение среднего времени безотказной работы (в.б.р.) основного элемента к среднему времени ремонта отказавшего элемента), при которым достигается заданный уровень стационарной надёжности 1 -п3 = {0.9; 0.99; 0.999} .Для анализа и сравнения результатов были выбраны следующие распределения: Экспоненциальное (M), Вейбулла-Гнеденко (WB), Логнормаль-ное (LN).

Рассмотрим частные случаи модели при р = 25;N = 3;NG = 100;Г = 1000; где р = a u b = i;

bi 1

T1 - время безотказной работы системы; T2 - время ремонта отказавшего элемента

Таблица 1. Значения относительной скорости восстановлении, при которых достигается заданный уровень стационарной надёжности

Table 1. Relative recovery rate values at which a given level of stationary reliability is achieved

T M (1/a,) WB (0.5) LN (1

0.9 0.99 0.999 0.9 0.99 0.999 0.9 0.99 0.999

M (1/b) 1.6 4.2 9.1 1.5 4.7 11.3 1.6 4.4 9.7

WB (0.5) 3.2 12.2 >25 2.9 11.6 >25 3.3 11.9 >25

LN (1) 1.6 3.9 7.6 1.6 4.5 9.7 1.6 4 7.2

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

Рисунке 2, 3, 4 представлены графики вероятности безотказной работы системы; а на рисунке 5 графики равномерной разности результатов имитационной модели.

Modern Information Technologies and IT-Education

<М_3|М|1> имитац. <М_3 WB|1> шитац. <M_3|LN|1> имитац.

10 15 20

р = еа/ев

Р и с. 2. Графики зависимости вероятности безотказной работы системы {MN / GI / от относительной скорости восстановления F i g. 2. The plot of the system uptime probability (MN / GI / versus the relative recovery rate

/

/

- iWB 3|М|1> имитац. - <WB_3|WB|1> имитац. --<WB 3|LN|1=-имитац.

I

10 15

p = ea/eb

Р и с. 3. Графики зависимости вероятности безотказной работы системы (WBN / GI /1 от относительной скорости восстановления F i g. 3. The plot of the system uptime probability l^fflB / GI / versus the relative recovery rate

- •¿LN_3|M|1^ имитац. <LN_3|WB|1> имитац. <LN 3|LN|1> имитац.

10 15 20

p = ea/eb

Р и с. 4. Графики зависимости вероятности безотказной работы системы {LNN ЙGI / от относительной скорости восстановления F i g. 4. The plot of the system uptime probability (LNN ЙGI / versus the relative recovery rate

Полученные результаты демонстрируют высокую асимптотическую нечувствительность стационарной надежности системы. Видно, что различия между кривыми при «быстром» восстановлении становятся исчезающе малыми для всех рассматриваемых распределений времени ремонта элементов системы. Например, уже начиная со значения р = 10 все кривые почти неотличимы.

<M_3|M|1>&&iWB_3|M|1s-<M_3|M|1>&&=LN_3|M|1» ■;М_3|М| 1 >&&<WB_3|WB| 1 > <М_3|М| 1 >&&<LN_3|LN|1 >

I-1—

10 15

p = ea/eb

20

25

Р и с. 5. Графики равномерной разности результатов имитационной модели в зависимости от р F i g. 5. The plot of the uniform difference in simulation model results versus the р

Графические результаты из Рис. 5 показывают, что равномерная разность между моделью /M/l) , (MN /LN/l}( и (LNn / LN/1 стремится к нулю с ростом р .

Модель расчета для оценки надёжности системы

{GIn / GI /1)

Определим следующие состояния моделируемой системы:

• состояние 0: один (главный) элемент работает, (N -1)- в холодном резерве;

• состояние 1: один элемент отказал и находится в ремонте, один - работает, (N - 2)- в холодном резерве;

• состояние 2: два элемента отказали, один находится в ремонте, второй ждёт его очереди на ремонт, один - работает, (N - 3)- в холодном резерве;

• состояние N -: все элементы отказали, то есть система останавливается.

Для описания алгоритма моделирования оценки надёжности системы {GIn / GI /l) потребуются следующие переменные:

• double t - часы модельного времени, меняются при отказе или восстановлении элементов системы;

• int i, j - переменные состояния системы; при наступлении события осуществляется переход из i, в j;

• double tnextfail - служебная переменная, в которой хранится время до следующего отказа элемента;

• double t „ - служебная переменная, в которой

nextrepair J г ' г

хранится время до следующего окончания ремонта отказавшего элемента;

• int k - счетчик числа итераций основного цикла.

Современные информационные технологии и ИТ-образование

Том 15, № 3. 2019 ISSN 2411-1473 sitito.cs.msu.ru

THEORETICAL QUESTIONS OF COMPUTER SCIENCE, COMPUTATIONAL MATHEMATICS, COMPUTER SCIENCE AND COGNITIVE INFORMATION TECHNOLOGIES

Для наглядности имитационная модель представлена графически на рисунке 6 в виде блок-схемы.

if i = 0 then

^nextrepair ' ^ ; J ' J + t ' tnextfail ;

else for (vinl: (N -1)) if г = v then

* := rf _GI(Я,);^ := rf _GI((x)") [;

nextfail

.— t + Sj ; tnextrepnii ~ ^ + S2 '

Р и с. 6. Блок-схема имитационной модели для оценки надёжности системы F i g. 6. Block diagram of a simulation model to assess the system reliability

Алгоритм 1. Псевдокод процесса имитационного моделирования системы (GIn / GI /1) Ввод: a1, b1, N, NG,"GI".

a1 - Среднее время между отказами элементов,

b1- Среднее время ремонта,

N - Число элементов в системе,

NG - Число Графиков траекторий,

"GI" - функция распределения.

Вывод: Оценка надёжности системы ET

Begin

array r []:=[0,0,0 J' // многомерный массив, содержащий результаты, k-го шага основного цикла

double t := 0.0; // инициализация часов модельного времени

int i := 0; j := 0; // переменные состояния системы

double tnextffail := 0.0; // переменная, в которой хранится время до следующего отказа элемента

double tnextiepA := 0.0; // переменная, в которой хранится время до следующего окончания ремонта

int k := 1; // счетчик числа итераций основного цикла

s := rf _ GI); // генерация произвольной случайной величины s - время до первого события (отказа)

ss := rf _ GI (<5(*))' // генерация произвольной случайной величины ss - время ремонта отказавшего элемента)

tnextfail t + S '

tnextrepair • t + SS ;

if tnextfail < ^nextrepair then

J •= J + t •= tnextfail '

else , j t , tnextrepair '

end

else i = N; then Break'

end

end

r [,,к] := [t,i, j];i := j;к := k +1;

end do

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

Оценка длительности пребывания в каждом состоянии i, i = 0,1,2, • • • N ' оценка надёжности

___ 1 NG

ET = ттттУ (длительность пребывания в состоянии N);

NG i=i

End.

Алгоритм 2 был программно реализован на языке R.

В Таблице 2 приведены значения оценки надёжности системы (оценки среднего в.б.р. системы) с временем, затраченным на моделирование. Для анализа и сравнения результатов были выбраны следующие распределения: Экспоненциальное (M), Вейбулла-Гнеденко (WB), и Логнормальное (LN). Однако разработанная имитационная модель не ограничена выбором распределений.

Рассмотрим частные случаи модели при a1 = 25 = 1;N = 3; \NG = 10000 ' Tl - Время безотказной работы системы' T2 - Время ремонта отказавшего элемента.

Таблица 2. Значения оценки надёжности системы ^GIn / GI /1 Table 2. System Reliability Assessment Values (^GI^ / GI /1

M (1/a, ) WB (0.5) LN (1)

M (1/b ) 16530.34 19566.77 25.18033

WB (0.5) 28.57675 927.8087 564.099

LN (1) 249458.5 71212.42 190780.8

while t < <x do

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

Modern Information Technologies and IT-Education

Заключение

Для восстанавливаемой системы (GIn / GI /l) холодного резервирования с одним восстанавливающим элементом, с произвольной ФР в.б.р. её элементов и произвольным законом распределения времени их ремонта были получены методом имитационного моделирования уровень надежности системы достигается при сравнительно небольшом превышении средних значении в.б.р. к времени ремонта и Графические и численные результаты показывают высокую асимптотическую нечувствительность стационарной надежности системы. Различия между кривыми при «быстром» восстановлении становятся исчезающе малыми для всех рассматриваемых распределений. Было показано, что самой надежной моделью является модель с логнормальным распределением времени безотказной работы и экспоненциальным распределением времени ремонта отказавшего элемента.

Список сокращений

ФР - функция распределения в.б.р. - время безотказной работы

Список использованных источников

[1] Уанкпо Г Ж. К., Козырев Д. В. Аналитическое и имитационное моделирование надежности замкнутой однородной системы с произвольным числом источников данных и ограниченными ресурсами для их обработки // Современные информационные технологии и ИТ-образование. 2018. Т. 14, № 3. С. 552-559. DOI: 10.25559/ SITITO.14.201803.552-559

[2] Ahmed W., Hasan O., Pervez U., Qadir J. Reliability Modeling and Analysis of Communication Networks // Journal of Network and Computer Applications. 2017. Vol. 78. Pp. 191215. DOI: 10.1016/j.jnca.2016.11.008

[3] Ometov A., Kozyrev D. V., Rykov V. V., Andreev S., Gaidamaka Y. V, Koucheryavy Y. Reliability-Centric Analysis of Offloaded Computation in Cooperative Wearable Applications // Wireless Communications and Mobile Computing. 2017. Vol. 2017. Article ID 9625687, 15 pages. DOI: 10.1155/2017/9625687

[4] Rykov V. V., Kozyrev D. V, Zaripova E. Modeling and Simulation of Reliability Function of a Homogeneous Hot Double Redundant Repairable System // Proceedings of the 31st European Conference on Modelling and Simulation ECMS2017 (May 23-26, 2017, Budapest, Hungary) / Z. Z. Paprika, P. Hora'k, K. Va'radi, P. T. Zwierczyk, 'A. Vi-dovics-Dancs, J. P. R'adics (eds.). Germany, Digitaldruck Pir-rot GmbH., 2017. Pp. 701-705. DOI: 10.7148/2017-0701

[5] Houankpo H. G. K., Kozyrev D. V. Sensitivity Analysis of Steady State Reliability Characteristics of a Repairable Cold Standby Data Transmission System to the Shapes of Lifetime and Repair Time Distributions of its Elements // CEUR Workshop Proceedings. 2017. Vol. 1995. Pp. 107-113. Selected Papers of the VII Conference "Information and Telecommunication Technologies and Mathematical Modeling of HighTech Systems", Moscow, Russia, 24-Apr-2017 / K. E. Samoui-lov, L. A. Sevastianov, D. S. Kulyabov (eds.). URL: http:// ceur-ws.org/Vol-1995/paper-15-970.pdf (дата обраще-

ния: 26.07.2019).

[6] Рыков В. В., Чан А. Н. О чувствительности характеристик надежности систем к виду функций распределения времени безотказной работы и восстановления их элементов // Вестник РУДН. Серия Математика. Информатика. Физика. 2014. № 3. С. 65-77. URL: https://elibrary. ru/item.asp?id=21757061 (дата обращения: 26.07.2019).

[7] Efrosinin D., Rykov V. Sensitivity Analysis of Reliability Characteristics to the Shape of the Life and Repair Time Distributions // Information Technologies and Mathematical Modelling. ITMM 2014. Communications in Computer and Information Science. Vol. 487 / A. Dudin, A. Nazarov, R. Yak-upov, A. Gortsev A. (eds). Springer, Cham, 2014. Pp. 101112. DOI: 10.1007/978-3-319-13671-4_13

[8] Efrosinin D., Rykov V V., Vishnevsky V On Sensitivity of Reliability Models to the Shape of Life and Repair Time Distributions // 2014 Ninth International Conference on Availability, Reliability and Security. Fribourg, 2014. Pp. 430-437. DOI: 10.1109/ARES.2014.65

[9] Rykov V V, Kozyrev D. V Analysis of Renewable Reliability Systems by Markovization Method //Analytical and Computational Methods in Probability Theory. ACMPT 2017. Lecture Notes in Computer Science. Vol. 10684 / V. V. Rykov, N. Singpurwalla, A. Zubkov (eds). Springer, Cham, 2017. Pp. 210-220. DOI: 10.1007/978-3-319-71504-9_19

[10] Rykov V. V, Kozyrev D. On Sensitivity of Steady-State Probabilities of a Cold Redundant System to the Shapes of Life and Repair Time Distributions of Its Elements // Statistics and Simulation. IWS 2015. Springer Proceedings in Mathematics & Statistics. Vol. 231 / J. Pilz, D. Rasch, V. Melas, K. Moder K. (eds). Springer, Cham, 2018. Pp. 391402. DOI: 10.1007/978-3-319-76035-3_28

[11] Parshutina S. A., Bogatyrev V A. Models to support design of highly reliable distributed computer systems with redundant processes of data transmission and handling // 2017 International Conference «Quality Management,Transport and Information Security Information Technologies» (1T&QM&1S). St. Petersburg, 2017. Pp. 96-99. DOI: 10.1109/1TMQ1S.2017.8085772

[12] Lisnianski A, Laredo D., Haim H. B. Multi-state Markov Model for Reliability Analysis of a Combined Cycle Gas Turbine Power Plant // 2016 Second 1nternational Symposium on Stochastic Models in Reliability Engineering, Life Science and Operations Management (SMRLO). Beer-Sheva, 2016. Pp. 131-135. DO1: 10.1109/SMRLO.2016.31

[13] Tourgoutian B., Yanushkevich A., Marshall R. Reliability and availability model of offshore and onshore VSC-HVDC transmission systems // 11th 1ET International Conference on AC and DC Power Transmission. Birmingham, 2015. Pp. 1-8. DO1: 10.1049/cp.2015.0101

[14] Cao J., Wang Q., Shen Y Research on modeling method of complex system mission reliability simulation // 2012 1nternational Conference on Quality, Reliability, Risk, Maintenance, and Safety Engineering. Chengdu, 2012. Pp. 307-311. DO1: 10.1109/1CQR2MSE.2012.6246242

[15] Gu Z., Zhu C., Shang L., Dick R. Application-Specific MPSoC Reliability Optimization // IEEE Transactions on Very Large Scale Integration (VLSI) Systems. 2008. Vol. 16, No. 5. Pp. 603-608. DO1: 10.1109/TVLS1.2008.917574

Современные информационные технологии и ИТ-образование

Том 15, № 3. 2019 ISSN 2411-1473 sitito.cs.msu.ru

[16] Huang W., Loman J., Song T. Reliability modeling of A warm standby redundancy configuration with active ^ standby ^ active units // 2014 Reliability and Maintainability Symposium. Colorado Springs, CO, 2014. Pp. 1-5. DOI: 10.1109/RAMS.2014.6798473

[17] Korolyuk V. S., Korolyuk D. Heuristic principles of phase merging in reliability analysis // Reliability: Theory and Applications. 2017. Vol. 12, No. 1. Pp. 66-71. URL: http:// www.gnedenko.net/Journal/2017/012017/RTA_1_2017-07.pdf (дата обращения: 26.07.2019).

[18] Гнеденко Б. В. О дублировании с восстановлением // «Известия АН СССР». Техническая кибернетика. 1964. № 5. С. 111-118.

[19] Rykov V. Multidimensional Alternative Processes Reliability Models // Modern Probabilistic Methods for Analysis of Telecommunication Networks. BWWQT 2013. Communications in Computer and Information Science. Vol. 356 / A. Dudin, V. Klimenok, G. Tsarenkov, S. Dudin (eds). Springer, Berlin, Heidelberg, 2013. Pp. 147-156. DOI: 10.1007/978-3-642-35980-4_17

[20] Петровский И. Г. Лекции по теории обыкновенных дифференциальных уравнений. 4-е изд. М. Л.: ГИТТЛ, 1952. 232 с.

[21] Гнеденко Б. В. О ненагруженном дублировании // «Известия АН СССР». Техническая кибернетика. 1964. № 4. С. 3-12.

[22] Соловьев А. Д. Резервирование с быстрым восстановлением // «Известия АН СССР». Техническая кибернетика. 1970. № 1. С. 56-71.

[23] Гнеденко Д. Б., Соловьев А. Д. Оценка надежности сложных восстанавливаемых систем // «Известия АН СССР». Техническая кибернетика. 1975. № 3. С. 121-128.

[24] Севастьянов Б. А. Эргодическая теорема для марковских процессов и ее приложение к телефонным системам с отказами // Теория вероятностей и ее применения. 1957. Т. 2, вып. 1. С. 106-116. URL: http://www. mathnet.ru/links/5 70f74a0f6c43 193 96 9a963341108352/tvp4960.pdf (дата обращения: 26.07.2019).

[25] Kalashnikov V. V. Geometric Sums: Bounds for Rare Events with Applications: Risk Analysis, Reliability, Queueing // Mathematics and Its Applications. Vol. 413. Dordrecht, Boston, London: Kluwer Academic Publishers, 1997. 256 p. DOI: 10.1007/978-94-017-1693-2

Поступила 26.07.2019; принята в печать 18.08.2018; опубликована онлайн 30.09.2019.

|об авторах:|

Уанкпо Гектор Жибсон Кинманон, аспирант кафедры прикладной информатики и теории вероятностей, факультет физико-математических и естественных наук, Российский университет дружбы народов (117198, Россия, г. Москва, ул. Миклухо-Маклая, д. 6), ORCID: http://orcid.org/0000-0002-5725-0313, [email protected]

Козырев Дмитрий Владимирович, доцент кафедры прикладной информатики и теории вероятностей, факультет физико-математических и естественных наук, Российский университет дружбы народов (117198, Россия, г. Москва, ул. Миклухо-Маклая, д. 6); Институт проблем управления им. В.А.

Трапезникова РАН (117997, Россия, г. Москва, ул. Профсоюзная,

д. 65), кандидат физико-математических наук, ORCID: http://

orcid.org/0000-0003-0538-8430, [email protected]

Все авторы прочитали и одобрили окончательный вариант рукописи.

References

[1] Houankpo H.G.K., Kozyrev D.V. Analytical Modeling and Simulation of Reliability of a Closed Homogeneous System with an Arbitrary Number of Data Sources and Limited Resources for their Processing. Sovremennye informacionnye tehnologii i IT-obrazovanie = Modern Information Technologies and IT-Education. 2018; 14(3):552-559. (In Russ., abstract in Eng.) DOI: 10.25559/SITITO.14.201803.552-559

[2] Ahmed W., Hasan O., Pervez U., Qadir J. Reliability Modeling and Analysis of Communication Networks. Journal of Network and Computer Applications. 2017; 78:191-215. (In Eng.) DOI: 10.1016/j.jnca.2016.11.008

[3] Ometov A., Kozyrev D.V., Rykov V.V., Andreev S., Gaidamaka Y.V., Koucheryavy Y. Reliability-Centric Analysis of Offloaded Computation in Cooperative Wearable Applications. Wireless Communications and Mobile Computing. 2017; 2017:9625687,15pp.(InEng.)DOI:10.1155/2017/9625687

[4] Rykov V.V., Kozyrev D.V., Zaripova E. Modeling and Simulation of Reliability Function of a Homogeneous Hot Double Redundant Repairable System. In: Paprika Z. Z., Hora'k P., Va'radi K., Zwierczyk P. T., Vidovics-Dancs 'A., R'adics J.P. (eds.). Proceedings of the 31st European Conference on Modelling and Simulation ECMS2017 (May 23-26, 2017, Budapest, Hungary) Germany, Digitaldruck Pirrot GmbH., 2017, pp. 701-705. (In Eng.) DOI: 10.7148/2017-0701

[5] Houankpo H.G.K., Kozyrev D.V. Sensitivity Analysis of Steady State Reliability Characteristics of a Repairable Cold Standby Data Transmission System to the Shapes of Lifetime and Repair Time Distributions of its Elements. In: Samouilov K.E, Sevastianov L.A., Kulyabov D.S. (eds.). CEUR Workshop Proceedings. 2017; 1995:107-113. Selected Papers of the VII Conference "Information and Telecommunication Technologies and Mathematical Modeling of High-Tech Systems", Moscow, Russia, 24-Apr-2017. Available at: http://ceur-ws. org/Vol-1995/paper-15-970.pdf (accessed 26.07.2019). (In Eng.)

[6] Rykov V.V., Ngia T.A. On sensitivity of systems reliability characteristics to the shape of their elements life and repair time distributions. Bulletin of the Peoples' Friendship University of Russia. Series "Mathematics. Information sciences. Physics". 2014; 3:65-77. Available at: https://elibrary.ru/ item.asp?id=21757061 (accessed 26.07.2019). (In Russ., abstract in Eng.)

[7] Efrosinin D., Rykov V. Sensitivity Analysis of Reliability Characteristics to the Shape of the Life and Repair Time Distributions. In: Dudin A., Nazarov A., Yakupov R., Gortsev A. (eds). Information Technologies and Mathematical Modelling. ITMM 2014. Communications in Computer and Information Science, vol. 487. Springer, Cham, 2014, pp. 101-112. (In Eng.) DOI: 10.1007/978-3-319-13671-4_13

[8] Efrosinin D., Rykov V.V., Vishnevskiy V. On Sensitivity of Reliability Models to the Shape of Life and Repair Time

Modern Information Technologies and IT-Education

Distributions. In: 2014 Ninth International Conference on Availability, Reliability and Security. Fribourg, 2014, pp. 430-437. (In Eng.) DOI: 10.1109/ARES.2014.65

[9] Rykov V.V., Kozyrev D.V. Analysis of Renewable Reliability Systems by Markovization Method. In: Rykov V., Singpurwalla N., Zubkov A. (eds). Analytical and Computational Methods in Probability Theory. ACMPT2017. Lecture Notes in Computer Science, vol. 10684. Springer, Cham, 2017, pp. 210-220. (In Eng.) DOI: 10.1007/978-3-319-71504-9_19

[10] Rykov V., Kozyrev D. On Sensitivity of Steady-State Probabilities of a Cold Redundant System to the Shapes of Life and Repair Time Distributions of Its Elements. In: Pilz J., Rasch D., Melas V., Moder K. (eds). Statistics and Simulation. IWS 2015. Springer Proceedings in Mathematics & Statistics, vol. 231. Springer, Cham, 2018, pp. 391-402. (In Eng.) DOI: 10.1007/978-3-319-76035-3_28

[11] Parshutina S.A., Bogatyrev V.A. Models to support design of highly reliable distributed computer systems with redundant processes of data transmission and handling. In: 2017 International Conference «Quality Management,Transport and Information Security, Information Technologies» (IT&QM&IS). St. Petersburg, 2017, pp. 96-99. (In Eng.) DOI: 10.1109/ ITMQIS.2017.8085772

[12] Lisnianski A., Laredo D., Haim H. B. Multi-state Markov Model for Reliability Analysis of a Combined Cycle Gas Turbine Power Plant. In: 2016 Second International Symposium on Stochastic Models in Reliability Engineering, Life Science and Operations Management (SMRLO). Beer-Sheva, 2016, pp. 131-135. (In Eng.) DOI: 10.1109/SMRLO.2016.31

[13] Tourgoutian B., Yanushkevich A., Marshall R. Reliability and availability model of offshore and onshore VSC-HVDC transmission systems. In: 11th IETInternational Conference on AC and DC Power Transmission. Birmingham, 2015, pp. 1-8. (In Eng.) DOI: 10.1049/cp.2015.0101

[14] Cao J., Wang Q., Shen Y. Research on modeling method of complex system mission reliability simulation. In: 2012 International Conference on Quality, Reliability, Risk, Maintenance, and Safety Engineering. Chengdu, 2012, pp. 307-311. (In Eng.) DOI: 10.1109/ICQR2MSE.2012.6246242

[15] Gu Z., Zhu C., Shang L., Dick R. Application-Specific MPSoC Reliability Optimization. IEEE Transactions on Very Large Scale Integration (VLSI) Systems. 2008; 16(5):603-608. (In Eng.) DOI: 10.1109/TVLSI.2008.917574

[16] Huang W., Loman J., Song T. Reliability modeling of A warm standby redundancy configuration with active ^ standby ^ active units. In: 2014 Reliability and Maintainability Symposium. Colorado Springs, CO, 2014, pp. 1-5. (In Eng.) DOI: 10.1109/RAMS.2014.6798473

[17] Korolyuk V.S., Korolyuk D. Heuristic principles of phase merging in reliability analysis. Reliability: Theory and Applications. 2017; 12(1):66-71. Available at: http://www. gnedenko.net/Journal/2017/012017/RTA_1_2017-07.pdf (accessed 26.07.2019). (In Eng.)

[18] Gnedenko B.V. On cold double redundant system with restoration. Izv. ANSSSR. Tekhnicheskaya Kibernetika. 1964; 5:111-118. (In Russ.)

[19] Rykov V. Multidimensional Alternative Processes Reliability

Models. In: Dudin A., Klimenok V., Tsarenkov G., Dudin S. (eds). Modern Probabilistic Methods for Analysis of Telecommunication Networks. BWWQT 2013. Communications in Computer and Information Science, vol. 356. Springer, Berlin, Heidelberg, 2013, pp. 147-156. (In Eng.) DOI: 10.1007/978-3-642-35980-4_17

[20] Petrovsky I.G. Lectures on the theory of ordinary differential equations. Moscow, GITTL, 1952. 232 pp. (In Russ.)

[21] Gnedenko B.V. On cold double redundant system. Izv. AN SSSR. Tekhnicheskaya Kibernetika. 1964; 4:3-12. (In Russ.)

[22] Solovev A.D. On reservation with quick restoration. Izv. AN SSSR. Tekhnicheskaya Kibernetika. 1970; 1:56-71. (In Russ.)

[23] Gnedenko D.B., Solovev A.D. Reliability evaluation of complex renewable systems. Izv. AN SSSR. Tekhnicheskaya Kibernetika. 1975; 3:121-128. (In Russ.)

[24] Sevast'yanov B.A. An Ergodic Theorem for Markov Processes and Its Application to Telephone Systems with Refusals. Theory of Probability & its Applications. 1957; 2(1):104-112. (In Eng.) DOI: 10.1137/1102005

[25] Kalashnikov V.V. Geometric Sums: Bounds for Rare Events with Applications: Risk Analysis, Reliability, Queueing. Mathematics and Its Applications, vol. 413. Dordrecht, Boston, London: Kluwer Academic Publishers, 1997, 256 pp. (In Eng.) DOI: 10.1007/978-94-017-1693-2

Submitted 26.07.2019; revised 18.08.2019; published online 30.09.2019.

Hector G.K. Houankpo, Postgraduate Student of the Department of Applied Probability and Informatics, Faculty of Science, Peoples' Friendship University of Russia (6 Miklukho-Maklaya Str., Moscow 117198, Russia), ORCID: http://orcid.org/0000-0002-5725-0313, [email protected]

Dmitry V. Kozyrev, Associate Professor of the Department of Applied Probability and Informatics, Faculty of Science, Peoples' Friendship University of Russia (6 Miklukho-Maklaya Str., Moscow 117198, Russia); V. A. Trapeznikov Institute of Control Sciences of RAS (65 Profsoyuznaya Str., Moscow 117997, Russia), Ph.D. (Phys.-Math.), ORCID: http://orcid.org/0000-0003-0538-8430, [email protected]

All authors have read and approved the final manuscript.

Современные информационные технологии и ИТ-образование

Том 15, № 3. 2019 ISSN 2411-1473 sitito.cs.msu.ru

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