Анализ и прогноз для блока ШБ3Бт

Для анализа функционирования исследуемых блоков использовались два метода математического моделирования: «Временные ряды» и «Марковские процессы».

а) Анализ на основе временных рядов

1) Был проанализирован технический паспорт № 555.4433.539т ПС на блок №110115 (изделие ШБЗБт), где зафиксированы движение изделия в эксплуатации и его поломки:

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

Отсюда можно составить таблицу роста числа поломок со временем:

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

«Простое среднее»

для общего

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

«Скользящее

среднее с линейным

трендом» для общего

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

«Экспоненциальное

сглаживание

с регулярным

трендом» для общего

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

«Линейная

регрессия» для общего

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

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

3) Прогноз поломок, связанных с отказом мгновенного расхода. Составляем таблицу числа поломок, относящихся к категории «отказы мгновенного расхода топлива»:

Составляем таблицу роста числа поломок со временем:

Вводим табличные данные в программу ПЭР в разделе «Экспоненциальное сглаживание с регулярным трендом» и получаем результаты в следующем виде:

1 2 3 4 5 6 7 8 9 10

4) Выводы: Сравнивая график общих поломок и график отказов мгновенного расхода, мы видим, что в течении всей эксплуатации прибора его общие поломки растут равномерно и что спрогнозированный участок кривой не имеет тенденции к резкому нарастанию числа отказов. На графике «отказов мгновенного расхода» мы видим, что в течении 8000 летных часов поломки, связанные с отказом мгновенного расхода практически не увеличивались, однако в последующем заметен их резкий рост в течении последующих 4000 летных часов, затем их рост продолжается и на прогнозируемом участке происходит значительное увеличение числа отказов. Также из сравнения видно, что в течении следующих 4000 летных часов, для которых производится прогноз, можно полагать, что в основном будут происходить неисправности, связанные с отказом мгновенного расхода топлива.

б) Анализ блока ШБТ3Бт с помощью теории Марковских процессов.

А) Проанализировав таблицу отказов, определили, что существует ряд типичных неисправностей. Введём обозначение S_i для состояний блока ШБТ3Бт, как системы, находящейся в любом из возможных состояний (как в исправном, так и в неисправных). Составим таблицу для видов неисправностей и их числа:

Б) В соответствии с частотой нахождения в основных неисправных состояниях рассчитаем вероятности нахождения в этих состояниях:

S2 = ( 11/30 ) * 100% = 36,67%

S3 = (6/30)* 100% = 20%

S4 = (6/30)* 100% = 20%

S5 = (7/30)* 100% = 23,33%

В) Эксперты-специалисты в области эксплуатации самолётных приборов дали следующую оценку для вероятностей переходов между состояниями:

На основании экспертных данных составляем общую таблицу для вероятностей переходов:

Г) Строим размеченный граф состояний для блока ШБТ3Бт :

Отказ датчика

мгновенного

расхода

0.05 0.3 0.0167

Отказ датчика 0.2 ИСПРАВЕН 0.15 Механические суммарного поломки

расхода

0.2333 0.05

НЕ РАБОТАЕТ

Д) Составляем стохастическую матрицу вероятностей переходов между возможными состояниями прибора:

Е) Осуществление прогноза будущих состояний.

Если в момент начала эксплуатации блок ШБТ3Бт полностью исправен, то начальное состояние системы представляется вектором {1; 0; 0; 0; 0;}. С помощью вектора начальных состояний {1; 0; 0; 0; 0;}. и стохастической матрицы вероятности переходов получаем целый ряд вероятностей нахождения системы в разных состояниях в последовательные моменты времени (время одной итерации - 2000 летных часов) :

Таким образом после 168 итераций получаем конечное состояние :

{0,0000; 0,0000; 0,5257; 0,0000; 0,4743}

На основании имеющейся таблицы итераций можно построить графики изменения состояний во времени :

Изменение состояния S1 (блок полностью исправен) со временем

Изменение состояния S2 (отказ датчика мгновенного расхода) со временем

Изменение состояния S3 (отказ датчика суммарного расхода) со временем

Изменение состояния S4 (механические поломки) со временем

Изменение состояния S5 (блок не работает) со временем