База ГОСТовallgosts.ru » 13. ОХРАНА ОКРУЖАЮЩЕЙ СРЕДЫ, ЗАЩИТА ЧЕЛОВЕКА ОТ ВОЗДЕЙСТВИЯ ОКРУЖАЮЩЕЙ СРЕДЫ. БЕЗОПАСНОСТЬ » 13.180. Эргономика

ГОСТ Р 51901.16-2017 Менеджмент риска. Повышение надежности. Статистические критерии и методы оценки

Обозначение: ГОСТ Р 51901.16-2017
Наименование: Менеджмент риска. Повышение надежности. Статистические критерии и методы оценки
Статус: Принят

Дата введения: 12/01/2018
Дата отмены: -
Заменен на: -
Код ОКС: 13.180
Скачать PDF: ГОСТ Р 51901.16-2017 Менеджмент риска. Повышение надежности. Статистические критерии и методы оценки.pdf
Скачать Word:ГОСТ Р 51901.16-2017 Менеджмент риска. Повышение надежности. Статистические критерии и методы оценки.doc


Текст ГОСТ Р 51901.16-2017 Менеджмент риска. Повышение надежности. Статистические критерии и методы оценки



ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ТЕХНИЧЕСКОМУ РЕГУЛИРОВАНИЮ И МЕТРОЛОГИИ

НАЦИОНАЛЬНЫЙ

СТАНДАРТ

РОССИЙСКОЙ

ФЕДЕРАЦИИ

ГОСТР

51901.16-

2017

(МЭК 61164:2004)

МЕНЕДЖМЕНТ РИСКА

Повышение надежности. Статистические критерии и методы оценки

(IEC 61164:2004, Reliability growth — Statistical test and estimation methods,

MOD)

Издание официальное

Москва

Стандартмиформ

2017

ГОСТ Р 51901.16—2017

Предисловие

1    ПОДГОТОВЛЕН Открытым акционерным обществом «Научно-исследовательский центр контроля и диагностики технических систем» (АО «НИЦ КД») на основе собственного перевода на русский язык англоязычной версии стандарта, указанного в пункте 4

2    ВНЕСЕН Техническим комитетом ло стандартизации ТК10 «Менеджмент риска»

3    УТ8ЕРЖДЕН И ВВЕДЕН В ДЕЙСТВИЕ Приказом Федерального агентства по техническому регулированию и метрологии от 12 сентября 2017 г. No 1059-ст

4    Настоящий стандарт является модифицированным по отношению к международному стандарту МЭК 61164:2004 «Повышение надежности. Статистические критерии и методы оценки» (IEC 61164:2004 «Reliability growth —Statistical test and estimation methods». MOD) путем внесения отклонений, объяснение которых приведено во введении к настоящему стандарту.

Международный стандарт разработан IEC 56.

Наименование настоящего стандарта изменено относительно наименования указанного международного стандарта для приведения в соответствие с ГОСТ Р 1.5—2012 (подраздел 3.5).

Сведения о соответствии ссылочных национальных стандартов международным стандартам, использованным в качестве ссылочных в примененном международном стандарте, приведены в дополнительном приложении ДА

5    ВЗАМЕН ГОСТ Р 51901.16—2005 (МЭК 61164:1995)

Правила применения настоящего стандарта установлены в статье 26 Федерального закона от 29 июня 2015 г. № 162-ФЗ «О стандартизации в Российской Федерации». Информация об изменениях к настоящему стандарту публикуется в ежегодном (ло состоянию на 1 января текущего года) информационном указателе «Национальные стандартыда. а официальный текст изменений и поправок — в ежемесячном информационном указателе «Национальные стандарты». В случае пересмотра (замены) или отмены настоящего стандарта соответствующее уведомление будет опубликовано в ближайшем выпуске ежемесячного информационного указателя «Национальные стандарты». Соответствующая информация, уведомление и тексты размещаются также в информационной системе общего пользования — на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет ()

© Стандартинформ. 2017

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

II

ГОСТ Р 51901.16—2017

Содержание

1    Область применения................................................................ 1

2    Нормативные ссылки................................................................ 1

3    Термины и определения............................................................. 1

4    Обозначения....................................................................... 2

5    Модели повышения надежности при проектировании и испытаниях......................... 5

6    Модели повышения надежности, используемые для систем/объектов на стадии

проектирования ................................................................... 6

7    Планирование повышения надежности и мониторинг испытаний на повышение

надежности объекта................................................................ 10

8    Использование степенной модели для планирования программ испытаний

повышения надежности............................................................. 14

9    Статистические критерии и методы определения оценок для непрерывной степенной модели ... 14 Приложение А (справочное) Примеры планирования и анализа моделей, используемых

при проектировании и испытаниях объекта...................................25

Приложение В (справочное) Степенная модель повышения надежности. Общая информация.....34

Приложение ДА (справочное) Сведения о соответствии ссылочных национальных стандартов международным стандартам, использованным в качестве ссылочных

в примененном международном стандарте...................................38

Библиография.......................................................................39

ГОСТ Р 51901.16—2017

Введение

В настоящем стандарте описаны в виде степенной функции модель повышения надежности времени и соответствующий метод прогнозирования и приведено руководство для их использования. Существует несколько моделей повышения надежности, однако степенная модель является одной из наиболее широко используемых. В настоящем стандарте установлены процедуры определения оценок некоторых или всех характеристик, перечисленных в ГОСТ Р 51901.6—2005 (разделы 4. 6 и 7).

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

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

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

В отличие от IEC 61164:2004 е раздел 2 настоящего стандарта не включен документ1), который нецелесообразно применять в соответствии с требованиями национальной системы стандартизации.

IEC 60050—191:1990 International Electrotechnical Vocabulary (IEV) — Chapter 191: Dependability and quality of service.

IV

ГОСТ Р 51901.16—2017 (МЭК 61164:2004)

НАЦИОНАЛЬНЫЙ СТАНДАРТ РОССИЙСКОЙ ФЕДЕРАЦИИ

МЕНЕДЖМЕНТ РИСКА

Повышение надежности.

Статистические критерии и методы оценки

Risk management Reliability growth. Statistical test and estimation methods

Дата введения — 2018—12—01

1    Область применения

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

2    Нормативные ссылки

В настоящем стандарте использованы нормативные ссылки на следующие стандарты:

ГОСТ Р 50779.26—2007 Статистические методы. Точечные оценки, доверительные, предик-ционные и толерантные интервалы для экспоненциального распределения

ГОСТ Р МЭК 60605-6—2007 Надежность в технике. Критерии проверки постоянства интенсивности отказов и параметра потока отказов

ГОСТ Р S1901.6—2005 (МЭК 61014:2003) Менеджмент риска. Программа повышения надежности

Примечание — При пользовании настоящим стандартом целесообразно проверить действие ссылочных стандартов в информационной системе общего пользования — на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет игы по ежегодному информационному указателю «Национальные стандарты», который опубликован по состоянию на 1 января текущего года, и по выпускам ежемесячного информационного указателя «Национальные стандарты» за текущий гад. Если заменен ссылочный стандарт, на который дана недатированная ссылка, то рекомендуется использовать действующую версию этого стандарта с учетом всех внесенных в данную версию изменений. Если заменен ссылочный стандарт, на который дана датированная ссылка, то рекомендуется использовать версию этого стандарта с указанным выше годом утверждения (принятия). Если после утверждения настоящего стандарта в ссылочный стандарт, на который дана датированная ссылка, внесено изменение, затрагивающее положение, на которое дана ссылка, то эго положение рекомендуется применять без учета данного положения. Если ссылочный стандарт отменен без замены, то положение, в котором дана ссылка на него, рекомендуется применять в части, не затрагивающей эту осыгжу.

3 Термины и определения

8 настоящем стандарте применены термины по ГОСТ Р 51901.6, а также следующие термины с соответствующими определениями:

3.1 целевая надежность (reliability goal): Желаемый уровень надежности, которым должна обладать продукция после завершения программы повышения надежности.

Издание официальное

1

ГОСТ Р 51901.16—2017

3.2    первоначальная надежность (initial reliability): Надежность продукции на ранних стадиях лро-актирования до сокращения потенциальных режимов отказа и причин их возникновения посредством улучшения конструкции объекта.

3.3    модель повышения надежности для стадии проектирования (reliability growth model for the design phase): Математическая модель, учитывающая возможные улучшения конструкции и их значение. позволяющая математически описать повышение надежности объекта от начала до конца этапа проектирования.

3.4    средняя интенсивность отказов (объекта) (average product failure rate): Средняя интенсивность отказов объекта, рассчитанная исходя из безотказности объекта за заданный период времени.

Примечание — Изменение интенсивности отказов во времени является результатом модификации конструкции объекта.

3.5    отсроченная модификация (delayed modification): Корректирующие изменения, которые введены в объект после окончания испытаний.

Примечание — Отсроченную модификацию не проводят в процессе испытаний.

3.6    коэффициент эффективности улучшения (improvement effectiveness factor): Коэффициент, на который уменьшилась интенсивность систематических отказов в результате корректирующей модификации.

3.7    испытания типа I, испытания, ограниченные по времени (type I test, time-terminated test): Испытания на повышение надежности, заканчивающиеся в установленное время, или испытания, для которых данные могут быть получены в течение времени, не зависящего от количества отказов.

3.8    испытания типа II, испытания, ограниченные по количеству отказов (type II test, failure-terminated test): Испытания на повышение надежности, которые завершаются после получения данных об установленном количестве отказов, или испытания, для которых данные могут быть получены в течение времени, зависящего от количества отказов.

4 Обозначения

В настоящем стандарте применены следующие обозначения,

а) Для 6.1. А.1 и В.З:

Т

*а0

d(f)

“О

D

h

t

к®

Rg(T)

R(t.T)

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

начальная средняя интенсивность отказов объекта за период проектирования: количество модификаций конструкции за время ( от начала процесса проектирования;

коэффициент повышения безотказности в результате уменьшения количества отказов;

общее количество внедренных улучшений конструкции;

общая продолжительность периода проектирования, в процессе которого внедрены улучшения конструкции;

переменная времени в процессе периода проектирования от 0 до f0:

средняя интенсивность отказов объекта как функция времени в течение периода

проектирования;

заданная средняя интенсивность отказов при завершении периода проектирования t0:

целевое значение вероятности безотказной работы объекта, которого необходимо достичь за период проектирования;

вероятность безотказной работы как функция времени и улучшений конструкции;

Ь) Для 6.2. А.2 и В.4:

RgT — целевое значение вероятности безотказной работы объекта, которого необходимо достичь за период проектирования;

2

ГОСТ Р 51901.16—2017

{о

а0

^'NS

О

К

J. *.>

ч*

D*

h

RAT)

РЦТ) кs

общая продолжительность периода проектирования;

коэффициент повышения надежности в течение периода проектирования;

интенсивность несистематических (или остаточных) отказов;

общее количество прогнозируемых или внедренных улучшений конструкции в про*

цессе проектирования, выполненных для устранения слабых мест:

общее количество различных категорий отказов:

индексы;

вероятность возникновения /-го слабого места конструкции в Л-й категории отказов за установленную наработку объекта;

ожидаемое количество слабых мест конструкции кчл категории отказов за установ* ленную наработку объекта:

общее количество прогнозируемых или внедренных улучшений конструкции за пе* риод проектирования для устранения слабых мест Л-й категории отказов; интенсивность отказов, относящихся к Л-й категории отказов; начальная вероятность безотказной работы объекта за время 7; вероятность безотказной работы объекта как функция времени 7; ожидаемое время достижения целевого значения показателя надежности;

с) Для 7.1.1. 7.1.2.9. АЛ В.1 и В.2:

tr

t

d(t)

ао

\«0

КМ)

«0<7)

m

«6

р

CV

d

I

U КА

К

м

—    общее количество улучшений конструкции, внедренных за период проекгирова* ния, для уменьшения количества выявленных слабых мест:

—    общая продолжительность периода проектирования, в процессе которого енедре* ны модификации конструкций;

—    переменная времени (в процессе проектирования 0 s tz tp);

—    количество модификаций конструкции за заданное время Г0, (0£fS ^);

—    коэффициент повышения безотказности в процессе периода проектирования;

—    начальная средняя интенсивность отказов объекта при проектировании;

—    переменная средняя интенсивность отказов объекта как функция времени в процессе проектирования (0 s t s (D);

—    начальная безотказность объекта за время 7 (назначенный ресурс или установленная наработка);

—    целевое значение вероятности безотказной работы объекта за время Г. которого необходимо достичь посредством улучшений конструкции;

—    повышение надежности объекта как функция времени и улучшения конструкции:

—    целевое значение средней интенсивности отказов:

—    установленное время в период срока службы объекта (назначенный ресурс, гарантийный период, срок службы);

—    параметр масштаба для степенной модели;

—    параметр формы для степенной модели;

—    критическое значение для проверки гипотез;

—    количество интервалов для анализа сгруппированных данных;

—    среднее и индивидуальные коэффициенты эффективности улучшений;

—    количество различных типов наблюдаемых отказов категории В;

—    индексы;

—    количество отказов категории А;

—    количество отказов категории В;

—    количество наблюдаемых отказов /-го типа категории В: Кв = К,;

—    параметр статистического критерия Крамера-Мизеса;

3

ГОСТ Р 51901.16—2017

N

Ч

т

Е[ЩТ)\

W-1):

КО

т

т

t

Ты

Г

X?<v>

р

АЪ

0(7)

flP

Р,

m

Q

Е

^NS

ок

Рч

(d

R,(T)

а

К

К;г

*,

Q

Я(Т)

0(7)

количество отказов:

количество отказов в r-м интервале.

накопленное количество отказов за время испытаний Г:

математическое ожидание суммарного количества отказов за время испытаний Т: границы /-го интервала времени испытаний для анализа сгруппированных данных:

время испытаний, суммарная наработка в процессе испытаний; суммарная наработка до /-го отказа;

общая суммарная наработка в процессе испытаний по отказам типа II; общая суммарная наработка в процессе испытаний по отказам типа I; квантиль х2 — распределения с v степенями свободы уровня у; общее обозначение параметра потока отказов; квантиль стандартного нормального распределения уровня; прогнозируемый параметр потока отказов;

параметр потока отказов в момент времени Т (соответствующий времени завершения испытаний);

текущее значение среднего времени между отказами: прогнозируемая средняя наработка между отказами:

вероятность успеха на стадии /-й модификации объекта при проектировании, описанная дискретной моделью повышения надежности Барлоу. Прошана и Шойера:

количество отказов неслучайного типа, оставшихся в момент времени Г > 0 в непрерывной модели повышения надежности IBM-Рознера;

доля доработанных отказов объекта/оборудования в соответствии с моделью повышения надежности IBM-Рознера; время воздействия на объект;

интенсивность несистематических (или остаточных) отказов; интенсивность систематических отказов; интенсивность отказов А-й категории:

количество потенциальных слабых мест конструкции в А-й категории отказов; вероятность того, чтоу-е потенциальное слабое место конструкции в к-й категории отказов вызовет отказ;

ожидаемое время проектирования для достижения целевого показателя безотказности. /?е(Т);

продолжительность стадии проектирования:

начальное значение вероятности безотказной работы за время Г;

параметр для представления ожидаемой скорости повышения безотказности;

ожидаемое количество слабых мест конструкции в А-й категории отказов, которые

вызовут отказ:

доля систематических (неслучайных) отказов проектируемого объекта в начале испытаний;

количество отказов объекта в начале испытаний;

доля отказов, устраненных посредством отладки в процессе испытаний на повышение надежности;

доля первичных отказов, устраненных за время t,

ожидаемое время для устранения части q систематических отказов за время испытаний;

суммарная наработка между отказами;

4

ГОСТ Р 51901.16—2017

d) Обозначения, используемые в дискретной модели повышения надежности, 7.2

Я,

7,= 1-Я, к

t, - Лу

т

тщ)\

>..р

—    безотказность или вероятность успеха г-й конфигурации:

—    вероятность отказа /• й конфигурации:

—    количество стадий и конфигураций;

—    количество испытаний на стадии /;

—    количество отказов на стадии г,

—    суммарное количество испытаний на стадии i:

—    общее количество отказов за время испытаний Г;;

—    среднее общее количество отказов за время испытаний ц

—    параметры масштаба и формы для стеленной и дискретной моделей.

5 Модели повышения надежности при проектировании и испытаниях

Основные принципы повышения надежности объекта одинаковы для процесса проектирования и испытаний. Это обусловлено тем. что оба этапа предполагают выявление и устранение слабых мест с целью улучшения объекта и позволяют оценить улучшение путем сравнения оценки показателя надежности с целевым значением. Разница заключается в методах, используемых для проведения анализа проектирования и испытаний, и моделях, используемых для оценки повышения надежности. В ГОСТ Р 51901.6 приведено руководство по разработке программ повышения надежности и методам анализа, используемым при проектировании и в испытаниях. 8 настоящем стандарте приведена подробная информация о моделях, которые могут быть использованы для оценки повышения надежности на разных стадиях жизненного цикла объекта и для различных типов объектов, как восстанавливаемых, так и невосстанавливаемых.

Математические модели повышения надежности построены для оценки достигнутого повышения и проектного уровня надежности. Модели повышения надежности направлены на поддержку планирования программ повышения надежности путем оценки количества и значимости изменений в процессе проектирования и разработки или в период испытаний, необходимых для достижения целевого значения показателя надежности.

Модели повышения надежности могут быть представлены с помощью интенсивности (или параметра потока) отказов или вероятности успешной работы за установленное время (безотказности), как показано на рисунке 1.

Интенсивность ОпсДОе    В*рдагт«ть бвЮгвванОВ работы

ЮпежпврулуниаииЯххматруящи    КЬвнесятс улучшений imvrpyiiiKM

Рисунок 1 — Планируемое увеличение средней интенсивности отказов или вероятности безотказной работы

5

ГОСТ Р 51901.16—2017

Внутри этой общей структуры существует много моделей повышения надежности. В таблице 1 приведен обзор основных категорий. Различия между проектированием и испытаниями, тип доступных данных влияют на выбор модели. К непрерывной категории относят объекты, функционирующие во времени, например ремонтируемые объекты. К дискретной категории относят данные, которые фиксируют лишь успех/отказ объекта в испытаниях, например, для невосстаиавливаемых изделий. Для оценки повышения надежности используют классические или байесовские процедуры. Первая категория использует только наблюдаемые данные, в то время как последняя использует как эмпирические данные проектирования и испытаний, так и инженерные знания, например в отношении ожидаемого количества рассматриваемых режимов отказов.

Таблица 1 — Категории моделей повышения надежности со ссылками на статьи

Тип модели

время

Непрерывное (время}

Дискретное (количество испытаний)

Проектирование

Классическая

6.1

Байесовская

6.2

Испытания

Классическая

7.1

7.2

Байесовская

Многие модели надежности разработаны для анализа данных испытаний. В настоящем стандарте представлена одна из самых популярных моделей повышения надежности, степенная модель {также известная как AMSAA или модель Кроу (Crow) [1]) в двух формах, непрерывной и дискретной. Эта модель является обобщением модели повышения надежности Дуэйна (Duane) согласно Кроу [1]11. Хотя существуют и байесовские варианты модели Кроу, они в стандарте не представлены. Обзор различных моделей повышения надежности, пригодных для анализа данных испытаний, можно найти у Джуэлла (Jewell) [21.(3] и Ce(Xie) (4].

Документальных отчетов о моделях повышения надежности, используемых при проектировании, очень мало. Поэтому были применены модель планирования повышения надежности, являющаяся модификацией степенной модели, для использования при проектировании и байесовский вариант модели IBM-Рознера. адаптированной для проектирования. Тем не менее они применимы только для объектов, непрерывно функционирующих во времени.

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

6 Модели повышения надежности, используемые для систем/объектов на стадии проектирования

6.1    Модифицированная степенная модель планирования повышения надежности на стадии проектирования объекта

6.1.1    Общая информация

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

Данную модель используют для планирования (а не для анализа данных), оценки количества или значимости улучшений первоначальной конструкции для повышения надежности объекта до целевого

6

Число в квадратных скобках указывает на номер литературного источника в библиографют.

ГОСТ Р 51901.16—2017

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

6.1.2 Модель планирования повышения надежности на стадии проектирования объекта

Модель планирования повышения надежности Красин (Krastch) [5] может быть выведена приведенным ниже образом. Пример такой модели, а также сравнительная таблица для ее упрощенного расчета и нанесения на график, заданные начальные и целевые показатели надежности объекта приведены в приложении А.

Если начальная вероятность безотказной работы объекта за установленную наработку Г была определена на основании анализа или испытаний и составляет R0(T). тогда, предполагая, что средняя интенсивность отказов является постоянной, начальная средняя интенсивность отказов объекта за наработку Травна

л _ И*о(П) Ло0--f-

(1)

Предположение о применимости степенной зависимости обосновано тем. что отказы, наиболее неблагоприятно влияющие на надежность, устраняют в первую очередь и благодаря модификации конструкции (сокращение количества отказов) параметр потока отказов объекта постоянно улучшается (снижается) с увеличением функции 4(0- Параметр потока отказов конструкции объекта в произвольный момент времени в течение периода проектирования имеет вид:

= (2)

где 4(0 — количество модификаций конструкции в процессе проектирования:

а0 — коэффициент повышения надежности в результате сокращения количества отказов;

D — общее количество внедренных улучшений конструкции;

t0 — общая продолжительность периода проектирования, пригодная для улучшения конструкции. С помощью линейного приближения количество модификаций конструкции, как функция времени, может быть линейно распределено в процессе периода проектирования:

«0 = 0-г-.    (3)

Средняя интенсивность отказов, как функция времени, в этом случае имеет вид:

KM-KoV+D

(4)

Если целевая средняя интенсивность отказов объекта задана целевым значением вероятности безотказной работы RG(T), то целевая средняя интенсивность отказов объекта в конце периода проектирования (0 приближенно равна:

7

(5)

в тоже время:

КМ аКо ■ (■1 ♦ о ■ -rf*ttD = Чо • <1

\мз(*0> =

„-iniRomi

(1 + D)-°e .

(6)

7

ГОСТ Р 51901.16—2017

Заменяя XaG(f0) выражением, содержащим целевое значение вероятности безотказной работы, и решая (6) относительно D. получаем:

D = ехр

-1п|

pnfWl

ИКо(7)1

-1.

(7)

Решая то же уравнение относительно скорости повышения надежности, выраженного в виде функции количества модификаций конструкции, начального и целевого значений вероятности безотказной работы, можно получить

а г. =

ln(1 + D)

(8)

В процессе проектирования непрерывное улучшение безотказности объекта, которую определяют за время Т. является функцией времени t (модель повышения надежности в период времени 0 s t s f0) и может быть записано в виде:

R(t. Т) - ехр(-Х .(f) • Г).

(9)

Подставляя выражение средней интенсивности отказов, может быть получена модель повышения надежности Красин для периода проектирования 0 < I < t0

R{t. Т) = exp|-Xa0 • |1 +

R(t Т) - expl

И*о(7)1

1 + 0 irn

(10)

(11)

В приведенном выше уравнении, выражая D через начальное и целевое значения вероятности безотказной работы, повышение безотказности, как функция времени для улучшения конструкции, имеет вид:

Ж*. Л = Я0(Г)

-о©

(12)

6.1.3 Отслеживание достигнутого повышения надежности

Отслеживание достигнутого повышения надежности подразумевает простой пересчет полученной ранее функции показателя надежности объекта в момент времени улучшения конструкции для учета изменений. Значение вероятности безотказной работы, рассчитанное для одного и того же установленного ресурса или гарантийного периода, просто наносят на график модели повышения надежности для соответствующего времени проектирования.

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

Использование анализа дерева отказов с помощью доступного программного обеспечения позволяет просто и быстро проводить и отслеживать автоматическое определение оценки безотказной работы объекта после изменений.

После завершения проектирования объекта и выполнения стадии валидации объекта плановые испытания на повышение надежности могут дальше повысить безотказность объекта или выявить режимы отказов, которые не были учтены при аналитическом определении оценок. Окончательная оценка вероятности безотказной работы завершенного проекта может служить целевым показателем безотказности для испытаний на повышение надежности.

8

ГОСТ Р 51901.16—2017

Пример практического вывода и применения модели планируемого повышения для улучшения безотказности на стадии проектирования приведен в приложении А. А.1.1. Этот реальный пример поэтапно показывает, как модель построена и как ее используют.

6.2 Модифицированная байесовская модель IBM-Рознера для планирования повышения

надежности на стадии проектирования

6.2.1    Общие положения

Модель представляет собой описание повышения надежности на стадии проектирования ремонтируемого объекта, подготовлена Куигли (Quigley) и Уоллсом (Walls) (6]—(8) и основана на адаптированной байесовской модели IBM-Рознера [9]. разработанной для анализа данных испытаний и описанной в 7.1.2.

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

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

Модель может быть использована двумя способами:

a)    для прогнозирования периода времени до достижения целевого значения показателя надежности посредством прогнозирования показателя надежности конструкции — это предполагает, что средняя скорость повышения надежности может быть рассчитана:

b)    для оценки скорости повышения надежности, необходимой для достижения целевого значения показателя надежности от начальной оценки в процессе установленного периода проектирования — это предполагает, что продолжительность стадии проектирования является фиксированной.

Детали, касающиеся математического описания модели, приведены в приложении В.

6.2.2    Требования к данным

В соответствии с целью применения необходимы данные, касающиеся целевого значения вероятности безотказной работы (ЯС(Г)). а также продолжительности стадии проектирования (г0) или средней скорости повышения надежности на стадии проектирования (а0), обеспечивающей достижение целевого значения.

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

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

Необходимо оценить вероятность того, что каждое слабое место в конструкции в каждой категории отказов может привести к отказу в течение установленного ресурса объекта. Оценка может быть основана на инженерных выводах. Вероятность возникновения отказа для /-го слабого места конструкции в категории отказов к записывают как рк{.

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

= 2$ Р*,-    (13)

где О* — ожидаемое общее количество слабых мест конструкции в категории к.

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

9

ГОСТ Р 51901.16—2017

6.2.3 Оценка повышения надежности и соответствующих параметров

Ниже приведены уравнения для расчета ключевых параметров модели повышения надежности. Начальное значение показателя вероятности безотказной работы конструкции за время Т рас* считывают по формуле

Я,(7) = exp{-[>.wsr+ if,, л*0 - e-Vffi.    (14)

Повышение безотказности конструкции за время Т рассчитывают по формуле

R,(T) = exp{-[/.NSr♦ if., лкв-“^(1 -е-^)]}.    (15)

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

а0-

ln

Sicin<((1-C'**r)

-lntffG(7)]-Xwsr

to

(16)

Если скорость повышения безотказности известна, время, необходимое для достижения целевого значения показателя безотказности, может быть рассчитано по формуле

£*«1^(1-в-*»7) -1п[Яс(Т)1->.^Г

(17)

6.2.4 Отслеживание повышения надежности на стадии проектирования

Отслеживание повышения безотказности с помощью модифицированной байесовской модели IBM-Рознера выполняют также, как это описано в 6.1.3.

7 Планирование повышения надежности и мониторинг испытаний

на повышение надежности объекта

7.1    Модели непрерывного повышения надежности

7.1.1    Степенная модель

В статистических процедурах для степенной модели повышения надежности в качестве исходных данных используют отказы и наработки в процессе испытаний. За исключением метода прогнозирования (см. 9.6). модель применяют к общему набору данных об отказах (см. ГОСТ Р 51901.6—2005. рисунки 2 и 4. характеристику 3) без деления на категории.

Основные уравнения для степенной модели приведены в настоящем разделе. Справочная информация о модели приведена в приложении В.

Математическое ожидание накопленного количества отказов за время Т вычисляют по формуле

£|№(Т)] = /.71* при X > 0. р > 0. Г> 0.    (18)

где X — параметр масштаба;

6 — параметр формы (функция общей эффективности улучшений соответствует повышению надежности. если 0 < р < 1: соответствует сохранению надежности, если 6 = 1: соответствует снижению надежности, если 6 > 1).

Параметр потока отказов до испытаний в момент завершения Г (час) имеет вид:

2(Г)= — Е{ЩТ))*к$Т*-' при Г>0.    (19)

10

ГОСТ Р 51901.16—2017

Таким образом, оба параметра. X и р. влияют на параметр потока отказов в заданное время. Уравнение представляет собой прямую с углом наклона касательной N(T) коси времени в точке 7. как показано на рисунке 9 в ГОСТ Р S1901.6—2005.

Среднее время между отказами за время испытаний Т имеет вид:

0(Г> = -!-.    (20)

*Т)

8 9.1 и 9.2 приведены оценки максимального правдоподобия для параметров X и р. В 9.3 приведены критерии согласия для модели, а в 9.4 и 9.5 рассмотрены процедуры построения доверительного интервала. Распространение модели на прогнозирование повышения надежности описано в 9.6.

Модель имеет следующие характерные особенности:

- Модель проста для определения оценок.

•    Если параметры оценены по прошлым программам, это удобный способ планирования будущих программ, использующих аналогичные условия испытаний и такую же результативность улучшения (см. 7 и ГОСТР S1901.6—200S. пункты 6.4.1—6.4.7).

•    Иногда модель дает нереальные значения; например. 2(7) = ж для 7 = 0, т. е. повышение может быть бесконечным и 2(7) стремится к нулю при 7. стремящемся к бесконечности. Однако эти ограничения не влияют на практическое использование модели.

•    Модель является относительно инертной и нечувствительной к повышению надежности сразу после корректирующей модификации и может давать заниженную (пессимистическую) оценку 6(7). если только она не используется для прогнозирования (см. 9.6).

•    Обычный метод оценки предполагает, что наблюдаемые наработки известны точно, но возможен альтернативный подход, когда отказы сгруппированы в пределах известного интервала времени (см. 9.2.2).

7.1.2 Модель для фиксированного количества отказов

Данная модель, также известная как модель IBM-Рознера {9], использует следующие предположения:

•    существуют случайные отказы (функция интенсивности отказов постоянна), возникающие с интенсивностью 2:

•    существует фиксированное, но неизвестное количество неслучайных конструктивных, производственных и других отказов объекта в начале испытаний.

Ограничение модели состоит в предположении о том. что коэффициент эффективности сокращения количества отказов равен единице.

Скорость изменения МЛ во времени пропорциональна количеству неслучайных отказов, оставшихся на момент времени 7:

ЩТ)

йТ

- -К2 ■ N(T),

(21)

следовательно.

МЛге-*»т<с.

Если количество отказов в точке 7s 0 равно К,. то

N{T) s К, • е-«гт если Т> 0. Kv К2 > 0.    (22)

Если £[М Л] — математическое ожидание совокупного количества отказов за время 7. тогда

Е1М Л) а 2 • 7 + К, • (1 - е-*> Л-    (23)

Это уравнение подразумевает, что к моменту времени 7 общее количество отказов равно сумме случайных и неслучайных отказов. Здесь £[М0))= 0.

При условии, что время стремится к бесконечности

ЩМЛ) - *-    (24)

11

ГОСТ Р 51901.16—2017

Так как модель является нелинейной, расчет параметров >.. К, и Kj следует выполнять итеративным методом.

Модель позволяет прогнозировать время устранения доли отказов «<?» (доли устраненных первоначальных неслучайных отказов 0 < <j < 1).

Количество неслучайных отказов, устраненных к моменту времени Г, равно

Л/(0) - N(t) = К, - К, • ег*гт.

(25)

Следовательно, часть первоначальных отказов, устраненных к моменту времени Г. имеет вид:

0(Л

1гег»г.1

(26)

Таким образом, получив оценку К2 (Х^). время, необходимое для достижения желаемого значения q неслучайных отказов, которые следует устранить, можно рассчитать по формуле

. tn(l-g)

*2 '

(27)

Количество неслучайных отказов, оставшихся к моменту времени Г. можно рассчитать по формуле К, • в~кг т.

Как и в других непрерывных моделях зависимой переменной является общее среднее время между отказами У(Г). где

0(Г)а---.    (28)

Общее количество отхаэов в (0.0

7.2 Дискретная модель повышения надежности

7.2.1 Описание модели

Данная модель разработана Я. Кроу (L.Crow) {1] и является дискретной версией степенной модели повышения надежности.

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

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

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

Предполагается, что конфигурация системы модифицируется на к этапах, а надежность остается постоянной в течение каждого этапа испытаний. На этапе / вероятность безотказной работы r-й конфигурации равна Rjt i - 1.... к и /?, > 0. Каждое испытание на стадии / имеет одинаковую вероятность Rуспеха. Данная модель обеспечивает определение оценки максимального правдоподобия Rk для каждой конфигурации / = 1.... к. Оценкой для финальной испытанной конфигурации является ftk. Если в

12

ГОСТ Р 51901.16—2017

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

Имеется к упорядоченных этапов, где этап 1 является первым, а этап к — последним. На этапе У данные состоят из т, отказов в испытаниях продолжительности п„ где У- 1.... к.

Пусть Г, = Х^, п- — общая продолжительность испытаний на всех этапах. Значит, после завершения Аго этапа данные содержат последовательность ?, дихотомических результатов испытаний, представляющих собой успех или отказ в последовательности испытаний объекта.

8 течение Аго этапа показатель надежности R, остается постоянным. Следовательно, для каждого испытания в процессе Аго этапа испытаний вероятность успеха равна /?,. а вероятность отказа равна

ft - 1 - Rt,    (29)

8 соответствии с данной моделью вероятность отказа для конфигурации г имеет вид:

(30)

где ^ = 0, л, * f, -и к > 0. р > 0 — параметры; л — параметр масштаба; р — параметр формы.

Пусть N(/,) — общее количество отказов за время испытаний У,, также пусть E[N(/,)] — математическое ожидание накопленного количества отказов за время испытаний tr Тогда согласно данной дискретной модели

E[W(f,)] = 2^ л,

(31)

Это характеристика повышения надежности, соответствующая степенной модели для непрерывных данных.

7.2.2 Оценка

Приведенные методы подходят, если испытания включают дискретные испытания в к последовательных этапах. Корректирующие действия для повышения надежности выполняют по окончании этапов, в результате чего повышается надежность следующей конфигурации. На стадии У вероятность безотказной работы Ай конфигурации равна Rt, /= 1,... к.

Данные содержат продолжительность испытаний л, для каждой конфигурации и соответствующее количество отказов т-. У = 1,... к. Рекомендуется выбирать mr> 0.

Показатель надежности Ri рассчитывают по формуле

где 7, =

_ >4-ktU

, ). > 0. э > 0.

(32)

Оценками максимального правдоподобия параметров >. и р являются величины, удовлетворяющие следующим уравнениям

(33)

К, 141 [SJ = 0.

(34)

H,*[l?ln(0-*kln{f,-1)],

(35)

(36)

П mi Л| — ГП/ ~\

(37)

",-4+4i-T

13

ГОСТ Р 51901.16—2017

Для определения оценки максимального правдоподобия параметров, к и 0. оценки вероятности отказа и безотказной работы для *-й конфигурации рассчитывают по формулам

и ^ - 1 - f? соответственно г = 1.....к.

Оценка надежности конфигурации системы на последнем этапе к равна Rk.

8    Использование степенной модели для планирования программ
испытаний повышения надежности

В качестве исходных данных для процедур, описанных в 6.4.2.3 и 6.4.4 ГОСТ Р 51901.6—2005, используют две величины, прогнозируемые с помощью моделей повышения надежности:

•    общая суммарная наработка в часах, необходимая для достижения целей программы;

•    математическое ожидание количества отказов за период испытаний.

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

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

9    Статистические критерии и методы определения оценок
для непрерывной степенной модели

9.1    Краткий обзор

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

Для испытаний типа I с ограниченной продолжительностью испытаний Г и испытаний типа II с ограниченным количеством отказов N используют разные формулы (см. 9.2.1}.

После завершения испытаний, описанных в 9.2.1 и 9.2.2. должны быть применены соответствующие критерии согласия (см. 9.3).

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

9.2    Проверка гипотез о повышении надежности и оценка параметров

9.2.1 Случай 1 — Известны наработки до каждого отказа

Метод применим только в том случае, если наработку регистрируют для каждого отказа.

Этап 1. Исключение посторонних отказов в соответствии с 6.4.5 ГОСТ Р 51901.6—2005 и/или другой соответствующей документацией.

Этап 2. Составление набора данных о наработках до отказа (в соответствии с 7.4 ГОСТ Р 50779.26—2007). Для испытаний типа I необходимо также учитывать время завершения испытаний.

Этап 3. Вычисление тестовой статистики 39

п,

(38)

(для испытаний типа 1);

(39)

14

ГОСТ Р 51901.16—2017

(для испытаний типа 2).

(40)

где N — общее количество отказов;

Т — общая продолжительность испытаний типа I;

TN — общая продолжительность испытаний типа II;

Г- — наработка до ко отказа.

В соответствии с гипотезой о нулевом повышении надежности (наработки между отказами соответствуют гомогенному процессу Пуассона) статистика U подчиняется стандартному нормальному распределению со средним 0 и стандартным отклонением 1. Статистика U может быть использована для проверки гипотезы о наличии положительного или отрицательного изменения надежности независимо от модели повышения надежности.

Двусторонний критерий для положительного или отрицательного изменения с уровнем значимости а имеет критические значения uUa/2 и -Ut-att’ ГДО — квантиль стандартного нормального распределения уровня (1-а/2) -100%.

Если справедливо одно из неравенств;

принимают решение о положительном или отрицательном изменении надежности соответственно и анализ продолжают с этапа 4.

Если справедливо неравенство

то принимают решение об отсутствии положительного или отрицательного изменения надежности с уровнем значимости а и анализ завершают. В этом случае гипотезу об экспоненциальном распределении наработок между последовательными отказами (гомогенности процесса Пуассона) принимают с уровнем значимости а. Критические значения и,.а/2 и ~и\~ „а соответствуют одностороннему критерию для положительного или отрицательного изменения надежности с уровнем значимости а/2.

Критические значения для двустороннего критерия с уровнем значимости 0.20 составляют 1.28 и (-1.28). Критическое значение 1,28 соответствует одностороннему критерию для положительного изменения с уровнем значимости 0.10. Для других уровней значимости можно выбрать критические значения по таблицам квантилей стандартного нормального распределения.

Этап 4. Вычисление суммы

(41)

(42)

S, =    1п(Т77^) (для испытаний типа 1);

S, = X*, HTNlTt) (для испытаний типа 2).

(43)

(44)

Этап 5. Вычисление несмещенной оценки параметра р по формуле

й W-1 .    „

р = —=—(для испытании типа 1);

(45)

*    N-2

Р = -g— (для испытаний типа 2).

(46)

Этап 6. Вычисление оценки параметра X по формуле

к - N1 (Г*)1’ (для испытаний типа 1); ).-NI (TNf (для испытаний типа 1).

(47)

(48)

15

ГОСТ Р 51901.16—2017

Этап 7. Вычисление оценки параметра потока отказов 2(7) и средней наработки между отказами 3(7*) при продолжительности испытаний 7* > 0 по формулам:

i(7) = 1р7*“';

(49)

б(Г)* 1/2(7).

(50)

Примечание 1 — Оценки ЦТ) и 3(7) являются оценками параметра потока отказов и средней наработки между отказами в момент времени 7 > 0 для 7 из диапазона представленных данных. «Экстраполированные» параметры для будущего времени 7 испытаний или до его завершения могут быть получены точно так же. но ими следует пользоваться осторожно из-за свойств экстраполяции. Экстраполяция не должна превышать продолжительность испытаний.

Примечание 2 — Если программа испытаний завершена, тогда оценка 3(7) в момент времени 7s 7 или 7 « 7Ы (в зависимости от типа испытаний) является значением средней наработки между отказами испытываемой конфигурации системы в конце программы испытаний.

9.2.2 Случай 2 — Наработки объединены в группы

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

Продолжительность испытаний соответствует интервалу (0; Г), она разделена на d интервалов

0 < 1(1) <■ t{2) < ... < ((d). нл интервал представляет собой период времени между t{i - 1) и    2.....d.

?(0) s О, ((d) = Т. величины /(/) могут принимать любые значения от 0 до Т.

Этап 1. Исключение посторонних отказов в соответствии с ГОС7 Р 51901.6—2005, пункт 6.4.5, и/или другой соответствующей документацией.

Этап 2. Составление набора данных о количестве рассматриваемых отказов Nзафиксированных в ьм интервале (f(r-1); ((/)). г = 1.... d.

Общее количество рассматриваемых отказов

(51)

Для каждого интервала произведение р, N не должно быть менее 5 (при необходимости смежные интервалы объединяют до начала испытаний), где

D

ДО)

(52)

Этап 3. Вычисление статистики U для d интервалов (после объединения при необходимости) и соответствующих значений

м-МО2
рА

(53)

В соответствии с гипотезой о нулевых изменениях надежности (когда наработки между отказами подчиняются гомогенному процессу Пуассона) статистика X2 распределена в соответствии с •/^распределением с (d -1) степенями свободы. Статистика X2 может быть использована для проверки наличия положительных или отрицательных изменений надежности независимо от модели изменений надежности.

Двусторонний критерий для положительного или отрицательного изменения с уровнем значимости а имеет критическое значение

CV=x^(d-D-    (54)

Если X2 i CV. принимают решение о положительном или отрицательном изменении надежности соответственно и анализ продолжают с этапа 4.

Если X2 < CV. то принимают решение об отсутствии положительного или отрицательного изменения надежности с уровнем значимости а и анализ завершают. В этом случае гипотезу об экспоненци

16

ГОСТ Р 51901.16—2017

альном распределении наработок между последовательными отказами (гомогенности процесса Пуассона) принимают с уровнем значимости а.

Критические значения (d -1) для различных уровней значимости а и числа степеней свободы {d-^) можно найти в таблицахквантилей /^распределения, например в стандартах ГОСТ Р 50779.26 и ГОСТ Р МЭК 60605-6.

Этап 4. Для первоначального набора данных, собранных на этапе 2. вычисляют оценки максимального правдоподобия для параметра формы р. Оценка максимального правдоподобия параметра р (р) удовлетворяет следующему уравнению:

Следует учитывать, что Г(0)=0. а также Г(0)1пГ(0) - 0. Все члены f() могут быть нормированы относительно Г(Р). тогда последний член ln(f(d)) исчезнет. Для решения этого уравнения относительно р необходимо использовать итеративный метод.

Этап 5. Определение оценки параметра X (X):

Этап 6. Определение оценок параметра потока отказов ЦТ) и среднего времени между отказами 6(7*) для продолжительности испытаний 7*> 0:

Примечание 1 — Оценки z(7) и д(7*) представляют собой оценки параметра потока отказов и средней наработки между отказами в момент времени 7 > 0 для Т из диапазона представленных данных. «Экстраполированные» параметры для большей продолжительности времени Т в процессе программы испытаний или до завершения испытаний могут быть получены точно так же. но ими следует пользоваться осторожно из-за свойств экстраполяции. Экстраполяция не должна превышать продолжительности испытаний.

Примечание 2 — Если программа испытаний завершена, тогда оценка 6(Г) для момента времени Г = f(d) является значением средней наработки между отказами испытываемой конфигурации системы в конце программы испытаний.

9.3 Критерий согласия

9.3.1    Общие положения

Если известны точные значения наработок, необходимо использовать случай 1. в противном случае необходимо использовать случай 2.

9.3.2    Случай 1 — Известны данные о наработках для каждого отказа

Для оценки параметра формы р сначала следует использовать метод, описанный в 9.2.1. затем вычислить статистику Крамера-Миэиса:

где М * Ы\лТ-Т — для испытаний типа I;

М=Л7-1иГ=7*^ — для испытаний типа II;

Г, < Г2 < ... <Ти.

8 таблице 2 приведены критические значения этой статистики с уровнем значимости 0.10. Если статистика С?(М) превышает критическое значение М. приведенное в таблице, то гипотеза о том. что степенная модель соответствует данным, должна быть отклонена. В противном случае гипотеза должна быть принята.

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

(55)

X = Ni t{df.

(56)

ЦТ) = хрг'И Ъ(Т)=иЦТ).

(57)

(58)

(59)

17

ГОСТ Р 51901.16—2017

При использовании графической процедуры оценку математического ожиданий наработки до /-го отказа £[7^] наносят на график в зависимости от наблюдаемой наработки до /-го отказа 7. 8 соответствии с приложением В оценку Е[Т^ можно рассчитать по формуле

тичном линейном масштабе, как показано на рисунке А.З. Визуальная согласованность этих точек с линией под углом 45* к оси ординат является субъективным критерием применимости модели.

9.3.3 Случай 2 — Наработки объединены в группы

Данный критерий применим только в случае, когда параметр ji оценивают на основе сгруппированных данных (см. 9.2.2). Среднее количество отказов в интервале времени (ф‘- 1); ф)) приближенно равно

Для каждою интервала значение е, не должно быть меньше 5. При необходимости смежные интервалы можно объединять до начала испытаний. Для d интервалов (после объединения при необходимости) и том же значении Nt, что и в 9.2.2. вычисляют статистику

Критические значения этой статистики с числом степеней свободы (d—2) можно найти в таблице ^•распределения, например в ГОС7 Р 50779.26 и ГОС7 Р МЭК 60605-6. Если критическое значение с уровнем значимости 0.10 превышено, то гипотеза о том. что степенная модель адекватно описывает сгруппированные данные, должна быть отклонена.

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

Для каждою интервала, ограниченного точкой 1(0 (от 0 до <(/)), количество наблюдаемых отказов составляет

(60)

/1

Значения 1-^-1 затем переносят на график в соответствии с наблюдаемой наработкой 7; в иден-

(61)

(62)

myi = tpyNf

Оценку математическою ожидания количества отказов E[N{t(i))) рассчитывают по формуле

E(W(K0)] =

(63)

(64)

Из этого следует

(65)

Графическая процедура состоит в построении графика

(66)

а также линии

In). + (р - 1)1п7. 7>0,

(67)

как это показано в примере (см. рисунок А.4). 18

ГОСТ Р 51901.16—2017

Зависимость между 1пл и 6, (р - 1) и « приведена в приложении В.

При р < 1 линия является убывающей. Визуальная согласованность точек с этой линией является субъективным критерием применимости модели.

9.4 Доверительные интервалы для параметра формы

9.4.1    Общие требования

Параметр формы (I в степенной модели повышения надежности характеризует наличие изменений надежности и их величину. Если 0 < р < 1. наблюдается повышение надежности, если р - 1. надежность не изменяется, а если р > 1. наблюдается снижение надежности.

При определении доверительного интервала для р, когда известны наработки до каждого отказа, используют случай 1. Для сгруппированных отказов необходимо использовать случай 2.

9.4.2    Случай 1 — Известны данные о наработках для каждого отказа

Этап 1. Вычисление оценки р в соответствии с этапом 5 в 9.2.1.

Этап 2.

а) Испытания типа I

Для определения двустороннего доверительного интервала для р с уровнем доверия 90 % необходимо вычислить величины DL и Ои:

Для двухстороннего доверительного интервала уровня 90 % для р

п _ y.o.os(2W) . Ul~ 2{W-1) '

D _ 4.95<2W> u 2<N-1) '

(68)

(69)

Квантили можно определить no таблицам /^-распределения. приведенным, например, в стандартах ГОСТР 50779.26 и ГОСТРМЭК 60605-6.

Нижняя доверительная граница для р имеет вид:

Кв = Di' Р-    <7°)

Верхняя доверительная граница для р имеет вид:

Poes<Vp.    (71)

Соответственно значения pta и рив являются также нижней и верхней односторонними доверительными границами для р с уровнем доверия 95 %.

Ь) Испытания типа I!

Для определения двустороннего доверительного интервала для р с уровнем доверия 90 %

, _xS,06(2(M-1)).

<■    2(W-2)    '

г

У.0,Зб(2<ЛГ-1))

2{N-2)

(72)

(73)

Нижняя доверительная граница для р:

Pte 5 dl ' Р-

(74)

верхняя доверительная граница для р:

Рув = 0ц' Э-

(75)

19

ГОСТ Р 51901.16—2017

Соответственно значения 0te и 0ив являются также нижней и верхней односторонними доверительными границами для 0 с уровнем доверия 95 %.

9.4.3 Случай 2 — Наработки объединены в группы

Данные процедуры построения доверительного интервала применимы в случае, когда оценка па* раметра 0 определена на основе сгруппированных данных в соответствии с 9.2.2.

Этап 1. вычисление 0 в соответствии с этапом 4 в 9.2.2.

Этап 2. вычисление P[i)

Р(0 г -нтгг. где /—1.2.....б.    (76)

Этап 3. вычисление значения величины А

_    [f#■ tnP(f>P1пЯ(/-t)p]2

\P(jf-P{i- !)Р]

А г I

Этап 4. вычисление значения С

С =

(77)

(78)

Этап 5. Определение приближенного двустороннего доверительного интервала для 0 с уровнем доверия 90 %

s-^.    (79,

где N — общее количество отказов.

Этап 6. Нижняя доверительная граница для 0:

Pie^O-S)-

(80)

верхняя доверительная граница для 0:

(WPO+S).

(81)

Соответственно значения 0te и 0ив являются также нижней и верхней односторонними доверительными границами для 0 с уровнем доверия 95 %.

9.5 Доверительные интервалы для средней наработки между отказами

9.5.1    Общие требования

В соответствии с 9.2.1 (этап 7) и 9.2.2 (этап 6) определяют оценку средней наработки между от* казами 6(Г). Для определения доверительных интервалов оценки 6(7). когда известны наработки до каждого отказа, следует использовать случай 1. Для сгруппированных отказов — случай 2.

9.5.2    Случай 1 — Известны данные о наработках для каждого отказа

Этап 1. вычисление 0(7) в соответствии с этапом 7 в 9.2.1.

Этап 2. Для определения двустороннего доверительного интервала с уровнем доверия 90 % находят значения L и U для соответствующего объема выборки N. используя таблицу 3 для испытаний типа I или таблицу 4 для испытаний типа II.

Этап 3. Нижняя доверительная граница для 0(7):

0te = L-6(7).    (82)

Верхняя доверительная граница для 0(7):

г и&(7>.

(83)

20

ГОСТ Р 51901.16—2017

Соответственно Gte и 0ив являются также нижней и верхней односторонними доверительными границами для 0(7*) с уровнем доверия 95 %.

9.5.3 Случай 2 — Наработки объединены в группы

Приведенные процедуры определения доверительного интервала применимы в случае, когда оценка р была определена по сгруппированным данным (см. 9.2.2).

Этап 1. Вычисление р в соответствии с 9.2.2 и й(Г) в соответствии с 9.2.1 (этап 7).

Этап 2. Вычисление величины P[i):

^0гтг5г. где/* 1.2.....d.    (84)

Этап 3. Вычисление значения А:

та [Ж/)Р • inP(iy-P(f- 1)Р • \npq-1 )*]2

Этап 4. Вычисление значения D:

(85)

(86)

Этап 5. Для определения приближенного двустороннего доверительного интервала для 0(7*) с уровнем доверия 90 % вычисляют значение S:

(1.64)0

(87)

где N — общее количество отказов.

Этап 6. Нижняя доверительная граница для 0(7*):

0te = б(ГК1 - S).

(88)

верхняя доверительная граница для 0(7*):

0ца = Й(7-Я1 ♦ S).

(89)

Соответственно значения 0LB и 0ив являются также нижней и верхней односторонними довери-тельными границами для 0(7*) с уровнем доверия 95 %.

9.6 Методика прогнозирования

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

Этап 1. Выделение отказов категории А и категории В (см. ГОСТ Р 57901.6—2005. определения 3.14 и 3.15).

Этап 2. Определение времени первого появления каждого типа отказов в категории В. составление отдельного набора данных. Пусть / — количество различных типов отказов категории В (I).

Этап 3. Выполнение этапов с 1 по 5 в соответствии с 9.2.1 в отношении этого набора данных для определения оценки параметра р, используя N »/ и Т или TN для полного набора данных.

Этап 4. Для каждого из / отказов различных типов категории В в наборе данных этапа 2 назначают коэффициент эффективности улучшения Er i - 1...../. Для каждого из / отказов различных типов ка

тегории В параметр Е, (0 s Е, s 1) является инженерной оценкой ожидаемого снижения потока отказов, вызванного идентифицированной корректирующей модификацией (см. 3.5).

вычисляют среднее Е этих значений или. если приемлемо, устанавливают средний коэффициент эффективности улучшения (например. 0.7) вместо индивидуальных значений Е„ i- 1.... /. как описано выше.

21

ГОСТ Р 51901.16—2017

Этап 5. Определение оценок прогнозируемого параметра потока отказов и средней наработки между отказами

Л М].    (90)

где КА — количество отказов категории А;

К, — количество наблюдаемых отказов но типа категории В,

Т-Т или T-TN6 соответствии с этапом 3. описанным выше.

Если значения Е, не назначены, а имеется только среднее Ё. то средний член в квадратных скобках становится равным

*в<1-£).    (91>

где Кв — количество отказов категории В.

В этом случае прогнозируемый параметр потока отказов

2р = т[Кл + Кв(1“^ + / $    *]■    (Э2>

Прогнозируемое значение средней наработки между отказами

0р = 1/2р.    (93)

Таблица 2 — Критические значения критерия согласия Крамера-Миэеса с уровнем значимости 0.10

м

Критическое значение статистики

3

0.154

4

0.155

5

0.160

6

0.162

7

0.165

8

0.165

9

0.167

10

0.167

11

0.169

12

0.169

13

0.169

14

0.169

15

0.169

16

0.171

17

0.171

18

0.171

19

0.171

20

0.172

30

0.172

260

0.173

Примечание —Для испытаний типа 1 М = N: для испытаний типа II М » N- 1.

22

ГОСТ Р 51901.16—2017

Таблица 3 — Двусторонние доверительные интервалы уровня доверия 90 % для средней наработки между отказами и для испытании типа I

N

1

и

N

L

и

3

0.175

6.490

21

0.570

1.738

4

0.234

4.460

22

0.578

1.714

5

0,281

3.613

23

0.566

1.692

6

0.320

3.136

24

0.593

1.641

7

0.353

2.626

25

0.600

1.653

6

0.381

2.606

26

0.606

1.635

9

0.406

2.444

27

0.612

1.619

10

0.428

2.317

28

0.618

1.604

11

0.447

2.214

29

0.623

1.590

12

0.464

2.130

30

0.629

1.576

13

0.480

2.060

35

0.652

1.520

14

0.494

1,999

40

0.672

1.477

15

0.506

1.947

45

0.689

1.443

16

0.521

1.902

50

0.703

1.414

17

0.531

1.661

60

0.726

1.369

16

0.543

1.625

70

0.745

1.336

19

0.552

1.793

80

0.759

1.311

20

0.561

1.765

100

0.783

1.273

Примечани е —Для W> 100

где Uq — квантиль стандартного нормального распределения уровня 100 - (0.5 + у/2) %.

23

ГОСТ Р 51901.16—2017

Таблице 4 — Двусторонние доверительные интервалы уровня доверия 90 % для средней наработки между отказами и для испытаний типа II

N

1

и

N

L

и

3

0.1712

4.746

21

0.6018

1.701

4

0.2587

3.825

22

0.6091

1.680

5

0.3174

3.254

23

0.6160

1.659

6

0.3614

2.892

24

0.6225

1.641

7

0.3962

2.644

25

0.6286

1.623

8

0.4251

2.463

26

0.6344

1.608

9

0.4495

2.324

27

0.6400

1.592

10

0.4706

2.216

28

0.6452

1.578

11

0.4891

2.127

29

0.6503

1.566

12

0.5055

2.053

30

0.6551

1.553

13

0.5203

1.991

35

0.6763

1.501

14

0.5337

1.937

40

0.6937

1.461

15

0.5459

1.891

45

0.7085

1.428

16

0.5571

1.876

50

0.7212

1.401

17

0.5674

1.814

60

0.7422

1.360

18

0.5769

1.781

70

0.7587

1.327

19

0.5857

1.752

во

0.7723

1.303

20

0.5940

1.726

100

0.7938

1.267

Примечани е —Для N> 100

где Uq 5*у2 — квантигь стандартного нормального распределения уровня 100 - (0.5 ♦ у/2) %.

24

ГОСТ Р 51901.16—2017

Приложение А (справочное)

Примеры планирования и анализа моделей, используемых при проектировании и испытаниях объекта

А.1 Планирование повышения надежности на стадии проектирования объекта

А.1.1 Пример планирования для степенной модели

Реальный пример, приведенный ниже, разъясняет, как построить стеленную модель и как сопоставить фактическое повышение надежности с запланированным.

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

•    целевые показатели надежности объекта:

•    оценка значимости и количества возможных модификаций конструкции:

•    продолжительность периода проектирования, в процессе которого могут быть проведены модификации.

А.1.2 Построение модели и мониторинг повышения надежности

Модель для планирования строят следующим образом.

Требуемый ресурс объекта должен составлять 15 лет:

7 =365 -24 -15 = 1.314- 10sч.

Заданная продолжительность периода проектирования составляет 140 дней. = 140 - 24 = 3,360 ч.

Вероятность безотказной работы объекта за первые 15 лет составляет R^T = 15 пет) = 0.72. Целевым показателем надежности объекта является вероятность безотказной работы за 15 пет. равная 0.95.

Исходя из значений начальной надежности и целевого показателя надежности, рассчитывают начальную и целевую среднюю интенсивность отказов:

Принимая во внимание предполагаемое количество возможных модификаций конструкции, скорость повышения надежности равна:

Фактическое повышение надежности объекта затем наносят на график и сравнивают с моделью (см. рисунок А.1).

Модель может быть подготовлена с использованием электронных таблиц со встроенными формулами (см. таблицу А.1).

Р0(П = вхр[-ла0Г] = 0.72: *5(7) =0.95:

Модель повышения надежности строят в соответствии с уравнением

Ж*. 7) = *0<7)

25

ГОСТ Р 51901.16—2017

Таблица А.1 — Расчет модели планирования повышения надежности на стадии проектирования

А

В

C

D

E

F

G

H

l

1

Символ или функция

“О

0

T.

4

ka0

'aG

*0.

4

2

Величина

0.774 42

10

0.72

0.95

131 400

2.50E-06

3.90E-07

3 360

3

Формула

SLN(LN(D2}/LN(E2))/ LN( 1+С2)

={-LN(D2)/F2)

=HN(E2)/F2)

4

Формула Я((.Г) model

=5DS2A(({S»$2+A4*24*(LN($ES2)/LN($DS2))A(-1/$B$2)-

-A4‘24)/$I$2»$B$2))

5

((days)

fittTL

model

ReW

R(lT)_

actual

6

0

0.72

0.95

7

5

0.771 577 8

0.95

8

10

0.805411 6

0.95

9

15

0.829 525 2

0.95

10

20

0.847 690 1

0.95

11

25

0.861 927 9

0.95

12

30

0.873 425 9

0.95

0.699

13

35

0.882 930 2

0.95

14

40

0.890 934 8

0.95

15

45

0.897 780 3

0.95

16

50

0.903 710 2

0.95

0.77

17

55

0.908 902 9

0.95

18

60

0.913 492 8

0.95

19

65

0.917 582 8

0.95

20

70

0,921 253 3

0.95

21

75

0.924 568

0.95

22

80

0.927 578 2

0.95

0.82

23

85

0.930 325 6

0.95

24

90

0.932 844 4

0.95

25

95

0.935 163 1

0.95

26

100

0.937 305 5

0.95

27

105

0.939 291 8

0.95

28

110

0.941 139 1

0.95

29

115

0.942 862 1

0.95

30

120

0.944 473 3

0.95

31

125

0.945 983 7

0.95

0.93

32

130

0.947 402 9

0.95

33

135

0.948 739 2

0.95

34

140

0.95

0.95

* В строках 3 и 4 приведены формулы из электронных таблиц Excel.

26

ГОСТ Р 51901.16—2017

В соответствии с примером улучшения проведены в три этапа: при этом несколько аналогичных модификаций были объединены в одно изменение. Вначале было проведено изменение типов конденсаторов на конденсаторы, диэлектрические свойства и надежность которых гораздо выше (первая модификация конструкции охватывала 106 конденсаторов с различными показателями и различным влиянием на надежность}.

Вторым изменением стало внедрение частей с более высокими показателями надежности, так как некоторые части (конденсаторы на полупроводниках) показали на испытаниях низкую надежность.

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

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

Вероятность безотказной работы

Нет необходимости е проверке соответствия фактической пинии стеленной модели. Пример является иллюстрацией тога, как фактические улучшения соответствуют стеленному закону (чем меньше крупных улучшений, тем круче кривая).

Рисунок А.1 — Планируемое и достигнутое повышение надежности — Пример

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

А.2 Пример байесовской модели повышения надежности на стадии проектирования объекта Целевой показатель надежности объекта (вероятность безотказной работы) Я^(7) на начало этапа проектирования составляет 0.95 за 2000 часов функционирования (Г). В данном случае существует 50 потенциально слабых мест конструкции (О). Каждое слабое место конструкции связано с категорией отказов. Существуют две категории отказов, а также 20 потенциально слабых мест конструкции, связанных с отказами категории 1 (Of), и 30 — связанных с отказами категории 2 (02). Для каждого потенциально слабого места существует соответствующая вероятность появления отказа (р^). Ожидаемое количество отказов, вызванных слабыми местами категории отказов 1. равно 7 (х,). ожидаемое кошчестео отказов категории 2 — 10 (Х2). По данным эксплуатации аналогичных объектов, один несистематический отказ возникает каждые 100000 часов работы, отказы категории 1 возникают один раз каждые 200 000 часов работы, а отказы категории 2 возникают один раз каждые 75 000 часов работы. Параметр повышения надежности равен 0.00633.

27

ГОСТ Р 51901.16—2017

Начальная оценка вероятности беэогхзэной работы объекта составляет ДД2000) = ехр(-|0.00001 • 2000 + 7 • (1 - e-0.0000©5 2000) + 10 • (1 _ ©-0.000013 2000))) =

= ехр<-(0.02 + 7 • (1 - е-001) + 10(1 - в"0 026666)]) = 0.702722.

Кривая надежности имеет вид:

R(() = ехрН0.02 + е’0-008333' (7 (1- е0 0’) +10 (1- в-0*026666)]} = ехр(-{0.02 + 0.332794 ■ е'0 000331]). График R(f) на стадии проектирования в днях приведен на рисунке А.2.

Я»

Рисунок А.2 — Плакируемое повышение надежности при использовании байесовской модели

повышения надежности

Этот график показывает прогнозируемое повышение надежности на протяжении года на стадии проектирования. Ожидается, что такая интенсивная стадия проектирования позволит достичь целевого значения показателя надежности через 284 дня.

Г 0.332794 "I 1л|_-1л(0.95)-0.02 J

1л *________=

0.00833

283.8079.

Если стадия проектирования должна занимать более короткий период времени, например 100 дней, необходима более интенсивная программа с параметром повышения надежности.

[    0.332794 "I

-4п(0.95)-0.02

и--— = 0.023641.

100

А.З Данные об отказах для дискретных испытаний

Ракетный комплекс проходит серию испытаний на повышение надежности на стадии проектирования, сосгояцую в общей сложности из 68 испытаний. Отложенные корректирующие действия были выполнены после 14.33 и 48-го испытания. Начиная с 49-годо 68-го испытания конфигурация не менялась. Можно выделить четыре этапа (к = 4) испытаний: Этап 1. Испытания 1—14.

Этап 2. Испытания 15—33.

Этап 3. Испытания 34—48.

Этап 4. Испытания 49—68.

Отсюда

л, =14. п2 = 19. л3 = 15. пл = 20.

1© = 0.1, = 14. t2 = 33. (3 = 48.    = 68.

28

ГОСТ Р 51901.16—2017

Данные об отказах для оценки параметров модели:

Конфигурация 1 выявила гл, = 5 отказов.

Конфигурация 2 выявила - 3 отказа.

Конфигурация 3 выявила т$ - 4 отказа.

Конфигурация 4 выявила тА - 4 отказа.

Оценки максимального правдоподобия для параметров Ji и Р составили: i = 0.595 ир® 0.780.

Эти оценки дают

Г, = 0.333.    ‘(2 = 0.234,    f3 = 0.206.    /4 = 0.190.

Показатели надежности на стадиях 1.2.3 равны: R, = 0.667, R2 = 0.766, R$ = 0,794, а показатель надежности на заключительной стадии испытаний равен RA - 0.8Ю.

А.4 Примеры повышения надежности в процессе испытаний А.4.1 Введение

Следующие числовые примеры демонстрируют использование процедур, описанных в разделе 9. В таблице А.1 приведен полный набор данных, иллюстрирующих методы повышения надежности, если известны наработки, а в таблице А.2 — если сгруппированы в интервалы. В таблицах А.З и А.4 приведены данные для методики прогнозирования, когда корректирующие модификации отлажены на конец испытаний. При необходимости в соответствии с 9.3 применяют критерии согласия. Эти примеры могут быть использованы для выбора соответствующих компьютерных программ, реализующих методы, описанные в разделе 9.

А.4.2 Определение оценок показателей надежности А.4.2.1 Общие положения

Набор данных в таблице А.2 соответствует испытаниям, заканчивающимся через 1000 часов. Эти данные использованы в примерах А.4.2.2 и А.4.2.3 для испытаний типа I и типа II соответственно. В сгруппированном виде эти данные продет зелены в таблице А.Э для примера А.4.2.4.

А.4.2.2 Пример 1 — Испытания типа I — Случай 1 — Известны наработки до каждого отказа Этот случай рассмотрен в 9.2.1. Данные таблицы А.2 соответствуют испытаниям, заканчивающимся через 1000 часов.

a)    Проверка гипотезы о повышении надежности

U = -3,713. Критические значения для двустороннего критерия с уровнем значимости 0.20 составляют 1.28 и (-1.28). Поскольку U < -1.28, это свидетельствует о повышении надежности, и анализ продолжают.

b)    Оценка параметров

Оценка параметров степенной ыодегы:

>.= 1.0694:

^ = 0,5623.

c)    Оценка средней наработки между отказами

Оценка средней наработки между отказами за 1000 часов составила 34,2 часа.

d)    Критерий согласия

С2(М) = 0.038 при М = 52. Для уровня значимости 0,10 критическое значение в соответствии с таблицей 2 составляет 0.173. Поскогъку С2(М) < 0,173. степенная модель является приемлемой (см. 9.3 и рисунок А.З).

e)    Доверительный интервал для р

Двусторонний доверительный интервал для р с уровнем доверия 90 % составляет (0.4491: 0.7101).

f)    Доверительный интервал для средней наработки между отказами

Двусторонний доверительный интервал для средней наработки между отказами за 1000 часов, соответствующий уровню доводил 90 %. составляет (24.2 ч; 48.1 ч).

А.4.2.3 Пример 2 — Испытания типа II — Случай 1 — Известны наработки до каждого отказа Этот случай рассмотрен в 9.2.1. Данные таблицы А.2 соответствуют испытаниям, заканчивающимся через 975 часов.

a)    Проверка гипотезы о повышении надежности

U = -3.764. Критические значения для двустороннего критерия с уровнем значимости 0.20 составляют 1.28 и -1,28. Поскольку U < -1,28. это свидетельствует о повышении надежности, и анализ продолжают.

b)    Оценка параметров

Оценка параметров степенной модвпч:

>. = 1.1067; р = 0,5594.

c)    Оценка средней наработки между отказами

Оценка средней наработки между отказами за 975 часов составила 33.5 часа.

d)    Критерий согласия

29

ГОСТ Р 51901.16—2017

С2(М) = 0,041 при М » 51. Для уровня значимости 0.10 критическое значение в соответствии с табгмцей 2 составляет 0.173. Поскольку С2(М) < 0.173, стеленная модель является приемлемой (см. 9.3 и рисунок А.З).

е) Доверительный интервал р

Двусторонний доверительный интервал для р с уровнем доверия 90 % составляет (0.4456:0.7080).

0 Доверительный интервал для средней наработки между отказами

Двусторонний доверительный интервал для средней наработки между отхазами за 975 часов, соответствующий уровню доверия 90 %. составляет (24.3 ч; 46.7 ч).

А.4.2.4 Пример 3 — Случай 2 — Сгруппированные данные

Этот случай описан в 9.2.2. Использованы данные таблицы А.2. В таблице А.З приведены данные об отказах, сгруппированных в интервалы по 200 ч. Анализ этого набора данных дает результаты, приведенные ниже.

a)    Проверка гипотезы о повышем** надежности

X2 = 14.7308 с четырьмя степенями свободы. Для уровня значимости 0,20 критическое значение составляет 6.0. Поскольку X2 > 6.0. это свидетельствует о том. что наблюдается положительное или отрицательное изменение надежности, и анализ продолжают.

b)    Оценка параметров

Оценки параметров степенной модели к- 0.9615;

Р = 0.5777.

c)    Оценка средней наработки между отхазами

Оценка сродней наработки между отказами за 1000 ч составила 33,3 ч.

d)    Критерий согласия

X2 = 2,175 с тремя степенями свободы. Для уровня значимости 0,10 критическое значение составило 6.25. Поскольку X2 < 6.25. степенная модель является приемлемой (см. 9.3 и рисунок А.4).

e)    Доверительный интервал для р

Двусторонний доверительный интервал для р с уровнем доверия 90 % имеет вид: (0,4083; 0.7471).

О Доверительный интервал для средней наработки между отказами

Двусторонний доверительный интервал для средней наработки между отказами, соответствующий уровню доверия 90 %. за 1000 часов имеет вид: (20.94 ч; 45.66 ч).

А.4.3 Прогнозируемые оценки показателей надежности А.4.3.1 Общие положения

Данный пример иллюстрирует расчет прогнозируемых оценок показателей надежности (см. 9.6). когда корректирующие модификации применены в коиде испытаний.

А.4.3.2 Пример 4

Основные данные, используемые в примере, приведены в таблице А.4. Имеется общее количество отказов N » 45 при количестве отказов категории А КА » 13. для которых не проводилось корректирующих модификаций. По истечении 4000 ч испытаний в объект было введено / = 16 различных корректирующих модификаций, направленных на устранение отказов Kg = 32 категории В. Категория для каждого отказа приведена в таблице А.4. В таблице А.5 представлена дополнительная информация, используемая для прогнозирования.

Процедура состоит из следующих этапов.

Этап 1. Идентификация категорий отказов А и В

Наработки для отказов категорий А и В указаны в таблице А.4. Наработки для 16 различных типов отказов категории В приведены в таблице А.5, столбец 2.

Этап 2. Идентификация первого появления отказов различных типов категории В

Время первого появления отказов категории В для 16 различных типов приведено в таблице А.5. столбец 3. Этап 3. Анализ данных первого появления отказа

Набор данных таблицы А.5 (столбец 3) проанализирован в соответствии с этапами 4—7 в 9.2.1. Результаты приведены тже.

Оценка параметров Оценка параметров степеююй модели X - 0.0326; р * 0.7472.

Оценка параметра потока первых отказов

Оценка параметра потока отказов категории В различных типов для первого появления отказов за 4000 ч составила 0,0030 ч~1.

Критерий согласия

С2(М) = 0.085 при М = 16. Для уровня значимости 0,10 критическое значение в соответствии с таблицей 2 составляет 0.171. Поскольку С2(М) < 0.171. степенная модель является приемлемой для наработок первого появления отказов различных типов категории В.

Этап 4. Наличие коэффициента эффективности

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

30

ГОСТ Р 51901.16—2017

Пример назначенных индивидуальных коэффициентов эффективности для каждой корректирующей модификации приведен в таблице А.5. столбец 5. Среднее этих 16 коэффициентов эффективности составляет 0.72. Среднее значение в диапазоне 0.65 — 0.75 является типичным {основано на опыте).

Этап 5. Определение прогнозируемого параметра потока отказов

Для прогнозирования параметра потока отказов необходимы значения следующих величин:

7 = 4000 ч;

КЛ = 13;

/= 16:

р = 0.7472;

ё =0.72;

К, — таблица А.5. столбец 4;

Е, — таблица А.5, столбец 5.

Оценка прогнозируемого параметра потока отказов в момент времени Т = 4000 ч (по завершении испытаний) составила 0.0074 чи.

Этап 6. Определение прогнозируемого значения средней наработки между отказами

Прогнозируемая оценка средней наработки между отказами составила 135.1 ч.

Примечание — Без повышения надежности в течение 4000 ч оценка средней наработки между отказами за этот период составила 68.9 ч (4000/45). Увеличение прогнозируемой оценки средней наработки между отказами вызвано введением 16 корректирующих модификаций с соответствующими коэффициентами эффективности. Необходимо учитывать чувствительность прогнозируемого значения средней наработки между отказами к назначенным коэффициентам эффективности. Если бы было назначено среднее значение коэффициентов эффективности 0,60. прогнозируемое значение средней наработки между отказами равнялось бы 121.3 ч. Среднее значение коэффициентов эффективности 0.80 дает прогнозируемое значение средней наработки между отказами 138,1 ч.

Таблица А.2 — Полные данные — Все уместные отказы и наработки для испытаний типа I

2

4

10

15

18

19

20

25

39

41

43

45

47

66

88

97

104

105

120

196

217

219

257

260

281

283

289

307

329

357

372

374

393

403

466

521

556

571

621

628

642

684

732

735

754

792

803

805

832

836

873

975

7* * 1000 ч. А/* 52.

В таблице А.З приведены данные о наработках в порядке возрастания.

Таблица А.Э — Сгруппированные данные для примера 3. полученные из таблицы А.2

Номер труппы

Количество отказов

время испытаний, соответствующее

правой точке интервала труппы, ч

1

20

200

2

13

400

3

5

600

4

8

800

5

6

1 000

Таблица А.4 — Полные данные для прогнозируемых оценок примера 4 — Все отказы и наработки

Наработки Г,и их классификация по категориям А и В. включая различные типы категории В. ч

т,

Категория

150

253

475

540

564

636

722

871

996

61

В2

ВЗ

В4

В5

А

В5

А

В6

Ъ

Категория

1003

1 025

1 120

1 209

1 255

1 334

1 647

1 774

1 927

В7

А

В8

В2

В9

В10

В9

В10

В11

31

ГОСТ Р 51901.16—2017

Окончание таблицы А.4

Наработки X, и их классификация по категориям А и в. включая различите типы категории В. ч

ъ

Категория

2 130 А

2 214 А

2 293 А

2 448 А

2490

В12

2 508 А

2 601 В1

2 635 В8

2 731 А

ъ

Категория

2 747 В6

2 650 В13

3 040 В9

3 154 В4

3171

А

3 206 А

3 245 В12

3 249 В10

3 420 В5

Ъ

Категория

3 502 ВЗ

3 646 В10

3 649 А

3 663 В2

3 730 В6

3 794 В14

3890

В15

3 949 А

3 952 В16

Г* = 4 000 ч, N = 45, КА = 13. Kg = 32.1 = 16.

Таблица А.5 — Различные типы отказов категории В. из таблицы А.4. с указанием наработок для первого появления. количества наблкдаемых отказов и коэффициентов эффективности

Столбец №

t

2

3

4

S

Тип

Время отказа, ч

Наработку для первого появления, ч

Количестве

Назначенный ко»ффициемт аффехтнвности

1

150:2 601

150

2

0.7

2

253; 1 209: 3 663

253

3

0.7

3

475: 3 502

475

2

0.8

4

540; 3 154

540

2

0.8

5

564: 722; 3 420

564

3

0.9

6

996: 2 747

996

2

0.9

7

1003

1 003

1

0.5

8

1 120:2635:3 730

1 120

3

0.8

9

1 255: 1 647; 3 040

1 255

3

0.9

10

1 334; 1 774:3 249: 3 646

1 334

4

0.7

11

1927

1 927

1

0.7

12

2 490: 3 245

2 490

2

0.6

13

2850

2 850

1

0.6

14

3 794

3 794

1

0.7

15

3890

3890

1

0.7

16

3952

3 952

1

0.5

32

ГОСТ Р 51901.16—2017

Mmwnwwow ожидании «работай. ч

Рисунок А. 3 — Диаграмма разброса для маге магического ожидания и наблюдаемых значений наработок на основе данных таблицы А.2 при использовании степенной модели

Коянаатво сшаол Ореме иопштннй

Рисунок А.4 — Отношение наблюдаемого и оцененного количества отказе» к наработке на основе данных таблицы А.2 при использовании степенной модели

33

ГОСТ Р 51901.16—2017

Приложение В (справочное)

Степенная модель повышения надежности. Общая информация

В.1 Постулат Дуайна

Наиболее полно используемая модель повышения надежности была описана Дж.Т. Дуэйном в 1964 году [10]. В этой работе Дуэйн анализировал данные об отказах для ряда объектов в процессе испытаний на стадии проектирования. Он заметил, что накопленное количество отказов Щ Г), деленное на накопленную за время испытаний наработку 7, уменьшается и стремится к прямой линии при построении графика в логарифмическом масштабе, т.е. приблизительно:

1л(М Т)/Т) = б - Ып7. где & > 0. а > 0.

Дуайн исследовал эти графики и пришел к выводу, что накопленное количество отказов аппроксимируется степенной функцией, т. в.

ЩТ) * Х71\ где а > 0. р = 1 -а.

Опираясь на эти наблюдения. Дуэйн представил текущее мгновенное значение параметра потока отказов в момент времени Т е следующем виде:

-N(T) =    где Г> 0.

dT

Это позволило вычислить мгновенное значение средней наработки между отказами:

Г1 I’

{XpTf-1)-1 --.где Г>0.

Хр

Показатель и = 1 - р иногда называют «скоростью роста».

Постулат Дуайна является детерминированным в том смысле, что он дает модегь повышения надежности, но не описывает изменчивость данных.

В.2 Степенная модель

ЛХ Кроу в 1974 году [1] изучил степенную модель изменений надежности и сформулировал основную вероятностную модель отказов как негомогенный процесс Пуассона (NHPP) {N(7), 7 > 0} со средним значением

щ/АГ)] = ;.т*

и функцией потока отказов z(7) = ХртР-'.

Степенная модель Кроу похожа на модель повышения надежности Дуайна. Обе эти модели используют одно и то же выражение лЛ1 для математического ожидания количества отказов в момент времени 7. Однако модель Кроу позволяет определить вероятность того, что N( 7) примет конкретное значение в следующем виде:

<i.7iy,e-*,T

Pr[W(7) = n]=---- при л = 0.1,2.....

л!

Также в соответствии с этой моделью

С(лф*Упри>=1.2.....

где Гу является наработкой до /-го отказа.

Это дает полезное приближение первого порядка

Щ7})=Ы "при/=1.2.....

для математического ожидания наработки до /-го отказа.

34

ГОСТ Р 51901.16—2017

Если р = 1. то z(7) • 7. и наработки между последовательными отказами подчиняются показательному распределению со средним 1/7. (гомогенный промесс Пуассона), надежность постоянна. Функция потока отказов г(Г) уменьшается при р < 1 (положительные изменения) и увеличивается при р > 1 (отрицательные изменения).

Степенная модель повышения надежности Кроу, являющаяся вероятностной интерпретацией постулата Дуайна. использует строгие статистические процедуры для метода оценки повышения надежности. Эти методы позволяют определить оценки максимального правдоподобия параметров модели и показателей надежности объекта. определить границы доверительных интервалов и применять критерии согласия. Степенная модель была расширена Кроу в 1983 году |11] для прогнозирования повышения надежности.

В.З Модифицированная стеленная модель для планирования повышения надежности на стадии проектирования объекта

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

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

Если начальная надежность объекта за установленное время Г равна ftyT). то начальная средняя интенсивность отказов, соответствующая моменту времени Г. имеет вид:

_    1п[Яо(Г)]

*t,0 “    у

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

>•«('> = *ц0 • И+<МГ“>

или

Ч,М-*во*И + о-у-)'вв.

Здесь </(/) — количество модификаций конструкции в любое время в период проектирования. Применяя линейную аппроксимацию, количество модификаций конструкции как функцию времени можно представить в виде:

10

Целевая средняя интенсивность потока отказов объекта при известном целевом показателе надежности RG{7) имеет вид:

-1п[Яе(7)]

А«6 --т-•

В процессе периода проектирования непрерывное повышение надежности объекта как функцию времени 1 (модель повышения надежности. (О S / S <0)) можно записать в виде:

Я(() = вхр(-Х(0Г).

Применив дифференцирование и подстановки, модель повышения надежности на стадии проектирования (О < t < t0) принимает вид:

*0» W>

|_° Ч 1п[ЯЬ(7)] I J

35

ГОСТ Р 51901.16—2017

В.4 Модифицированная байесовская модель планирования повышения надежности IBM-Рознера на стадии проектирования

Эта модификация модели повышения надежности IBM-Рознера изначально была предназначена для анализа данных испытаний (8) и [6]. Версия, представленная в данном приложении, адаптирована для решений по планированию на стадии проектирования объекта.

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

Как и модель IBM-Рознера [9]. данная модель предполагает, что для объекта характерно фиксированное количество слабых мест или потенциальных отказов и в период между внедрениями модификаций конструкции интенсивность отказов постоянна. Кроме того, предполагается, что изменения конструкции являются лучшим способом устранения слабых мест.

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

Предполагается, что несистематические отказы возникают с постоянной интенсивностью (X^g) и могут быть рассчитаны с помощью данных об объектах с аналогичной конструкцией или на основании инженерных расчетов.

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

Для оценки влияния систематических отказов на надежность конструкции объекта необходимо идентифицировать по К категориям отказов все (О) возможные слабые места конструкции. Вероятность того, что каждое слабое место конструкции в каждой категории отказов приводит к отказу в процессе эксплуатации, необходимо оценить с помощью, например, тжвнерных заключений. Могут потребоваться применение методов выявления слабых мест и оценка вероятности возникновения соответствующего отказа с помощью инженерных выводов (см. например. [7]).

Прогнозируемое количество слабых мест в категории отказов Л{цк). вероятнее всего приводящих к отказу при отсутствии модификаций конструкции, может быть рассчитано по формуле

где Dk — общее количество слабых мест конструкции, ожидаемое в категории отказов к и р^, — вероятность того, что у-е слабое место конструкции приведет к отказу категории отказов к. Расчет основывается на предположении о том. что количество слабых мест для каждой категории отказов является случайной величиной, подчиняющейся распределению Пуассона.

Значение интенсивности систематических отказов также должно быть известно для каждой категории отказов. Этот параметр можно определить с помощью эмпирических или общих данных об аналогичных существующих объектах.

После того как исходные данные определены, может быть определена апостериорная оценка надежности первоначальной конструкции объекта. Эго произведение надежности несистематических отказов и надежности систематических отказов. Показатель первого — это произведение (априорного) распределения количества слабых мест и (эмпирических) данных о систематических отказах. Таким образом, надежность первоначальной конструкции имеет вид:

Rf<r) = expHXwsr+ П*(1 - «г'*1)])-

Предполагая, что модификации внедряют для устранения слабых мест, надежность конструкции объекта возрастет. Таким образом, учитывая скорость повышения надежности (ад), надежность модифицированной конструкции объекта в момент времени Г равна

R(T) = exp{-[Xws7+

Для определения оценки скорости повьаиения целевого показателя надежности необходимо заменить (Rq(T)) и установленное время периода проектирования (10) на Я(Г) и Те формуле а0 (в предыдущем уравнении). Перестановка дает

36

ГОСТ Р 51901.16—2017

Если скорость повышения надежности установлена или рассчитана, аналогично определяют оценку среднего времени до достижения целевого показателя надежности:

37

ГОСТ Р 51901.16—2017

Приложение ДА (справочное)

Сведения о соответствии ссылочных национальных стандартов международным стандартам, использованным в качестве ссылочных в примененном международном стандарте

Таблица ДА.1

Обозначение ссылочного национальною стандарта

Степень

соответствия

Обозначение и наименование соответствующего международного стандарта

ГОСТ Р 50779.26—2007

МОО

IEC 60605-4:2001 «Испытания оборудования на надежность. Часть 4. Статистические процедуры для экспоненциального распределения. Точечные оценки, доверительные, толерантные и предикционные интервалы»

ГОСТ Р МЭК 60605-6—2007

ЮТ

IEC 60605-6:2007 «Испытания оборудования на надежность. Часть 6. Критерии проверки постоянства интенсивности отказов и параметра потока отказов и методы их оценки»

ГОСТ Р 51901.6—2005 (МЭК 61014:2003)

MOD

IEC 61014:2003 «Программа повышения надежности»

Примечание — В настоящей таблице использованы следующие условные обозначения степени соответствия стандартов:

-    IDT — идентичный стандарт;

-    МОО — модифицированный стандарт.

38

ГОСТ Р 51901.16—2017

Библиография

[1]    CROW, L.H.. 1974. Reliability Analysis for Complex Repairable Systems. Reliability and Biometry, ed. F. Proschan and R.J. Serfling. pp. 379—410. Philadelphia. PA: SIAM

[2]    JEWELL. W.A. (1984) General framework for learning curve reliability growth models. Operations Research. 32. 547—558

[3]    JEWELL. W. (1985) Bayesian Extensions to a Basic Model of Software Reliability. IEEE Trans, on Software Engineering. 11. 1465—1471

[4]    XIE. M (1993) Software Reliability Models a Selected Annotated Btbliography. Software Testing. Verification and Reliability. Vol. 3. 3*28

[5]    KRASICH. Milena (1999) Analysts Approach to Reliability Improvement. Proceedings of the Annual Technical Meeting of the Institute of Environmental Sciences and Technology, pp. 180—188. Ontario. CA

[6]    QUIGLEY. J.L. and WALLS. LA.. (2003) Robust Estimation Procedures for Order Statistic Reliability Growth Models. IEEE Transactions in Reliability

[7]    WALLS. L. and QUIGLEY. J. (2001) Building Prior Distributions to Support Bayesian Reliability Growth Modelling Using Expert Judgement. Reliability Engineering and System Safety. 74. 117—128

[8]    QUIGLEY. J. and WALLS L.A. (1999) Measuring the Effectiveness of Reliability Growth Testing. Quality and Reliability Engineering International. 15. 87—93

[9]    MIL-HDBK-189, Управление ростом надежности. 13 Февраль 1981 (MIL-HDBK-189. Reliability Growth Management. 13 February 1981)

(10]    DUANE. J.T. (1964) Learning Curve Approach to Reliability Monitoring. IEEE Transactions on Aerospace 2: pp. 563—566

(11]    CROW. L.H.. 1983. Reliability Growth Projection From Delayed Fixes. Proceedings of the 1983 Annual Reliability and Maintainability Symposium, pp. 84—89. Orlando. FL

39

ГОСТ Р 51901.16—2017

УДК 658.562.014:006.354    ОКС 13.180    Т59

Ключевые слова: менеджмент риска, анализ надежности, показатели надежности, параметр потока отказов, модель повышения надежности

40

БЗ 9—2017/126

Редактор Л.В. Коретмиковэ Технический редактор В.Н. Прусакова Корректор П.С. Лысенко Компьютерная верстка А.А. Ворониной

Сдано в набор 13.09.201? Подписано о лечат» 16.t0.2017. Формат 60*84 Vj. Гарнитура Ариап. Уел. печ. п. 6.12. Уч.-иэд. л. 4.63 Тираж 22 экэ. За» 1969.

Подготовлено на основа электронной версии, предоставленной разработчиком стандарта

Издано и отпечатано во ФГУП «СТАНДАРТИНФОРМ». 123001 Моема. Гранатный пар., 4.