АНАЛИЗ ОБЕСПЕЧЕНИЯ ЭКСПЛУАТАЦИОННОЙ ГОТОВНОСТИ ГРУППЫ СУДОВ СМЕННО-ЗАПАСНЫМИ ЧАСТЯМИ С ИСПОЛЬЗОВАНИЕМ ВЕРОЯТНОСТНОГО МОДЕЛИРОВАНИЯ
Аннотация и ключевые слова
Аннотация (русский):
Разработана модель обеспечения эксплуатационной готовности однотипных изделий группы судов, использование которой позволяет скорректировать объёмы начального запаса сменно-запасных частей в местах базирования флота, ремонтных и производственных мощностей, заказа сменно-запасных частей для группового запаса. Процесс эксплуатации однотипных изделий группы судов предложено представить как совокупность процессов потребления и пополнения сменно-запасных частей и процессов вывода и последующего ввода изделий в эксплуатацию. При этом применялось последовательное рекуррентное решение полученной системы уравнений (модели) для каждого из условных этапов рассматриваемого периода эксплуатации методом Рунге - Кутта. Учёт распределений отклонений от их среднего значения реализован с использованием статистического моделирования. Показано успешное применение описанного подхода к задаче обеспечения эксплуатационной готовности конкретной группы судов. Использование разработки способствует решению обозначенных проблем при их внедрении в организациях, осуществляющих эксплуатацию флота.

Ключевые слова:
суда, эксплуатационная готовность, сменно-запасные части, модель, вероятностный подход
Текст
Введение Для периода эксплуатации судовых систем и оборудования существует проблема определения количества группового запаса сменно-запасных частей в местах базирования флота и необходимой интенсивности его пополнения, что осложняется особенностями эксплуатации судов, стратегий восстановлений и ремонтов. Запас должен перекрывать потребности флота в сменно-запасных частях, обеспечивая его готовность к использованию, и одновременно с этим быть минимальным из-за необходимости сокращения издержек на его содержание. Цель исследования - обеспечение эксплуатационной готовности группы судов путём определения необходимого на период эксплуатации начального запаса сменно-запасных частей, а также значений величин, характеризующих процессы восстановления изделий и пополнения запаса. Публикации по данной проблеме [1-5] позволяют судить о недостаточной степени её изученности и проявляемом к ней значительном научном и практическом интересе. Результатом анализа процессов эксплуатации судов и опубликованных работ по данной тематике стал вывод о том, что существующие модели в полной мере не позволяют учитывать сложные стратегии потребления и пополнения сменно-запасных частей, часто встречающиеся в практике эксплуатации флота. Например модель, приведённая в [3], хоть и позволяет учитывать число восстановлений изделий, но не может учесть ни интенсивность восстановлений, ни возможность пополнения запасных частей в рассматриваемый период за счёт их заказа или изготовления. Не позволяют сделать этого и модели, приведённые в [4, 5]. Превышение потребностями флота в деталях среднего значения учитывается введением показателя уровня достаточности φi, который по своей сути имеет смысл коэффициента готовности ЗИП (запасные части, инструменты и принадлежности - ГОСТ 2.601), но с учётом последствий отказов деталей. Некоторые из известных моделей, в силу принятых допущений, не могут быть использованы в случаях с нестационарными потоками событий. Модель процесса обеспечения эксплуатационной готовности судов сменно-запасными частями и её применение В ходе исследования нами была разработана модель процесса обеспечения эксплуатационной готовности судов сменно-запасными частями, обладающая достаточной универсальностью в применении и позволяющая учесть сложные стратегии потребления и пополнения сменно-запасных частей. В результате подробного изучения статистических данных об эксплуатации групп судов при разработке модели были приняты следующие допущения: - в процессе эксплуатации любой элемент может находиться либо в работоспособном состоянии, либо в состоянии ремонта; - потоки событий процесса обеспечения работоспособности изделий обладают свойствами ординарности и отсутствия последействия; - потоки событий процесса обеспечения работоспособности изделий при некоторых условиях стационарны; - число состояний в этих процессах дискретно, а время непрерывно. Совокупность этих свойств позволила для периодов, на которых потоки событий стационарны, реальный процесс обеспечения работоспособности изделий аппроксимировать марковским процессом типа гибели и размножения с непрерывным временем и дискретным числом состояний [6]. При отказе изделий для восстановления их работоспособности необходимо затратить сменно-запасные части из запаса. Этим обусловлено потребление сменно-запасных частей. Пополнение запаса может производиться в результате их изготовления, ремонта или заказа. Указанные состояния и переходы между ними представлены размеченными графами (рис. 1). Рис. 1. Графы процессов: а - поддержания работоспособности изделий; б - потребления и пополнения запаса сменно-запасных частей Введем обозначения: x1 - состояние процесса, в котором i-е изделие находится в работоспособном состоянии; x0 - состояние процесса, в котором i-е изделие находится в ремонте в результате отказа; y1 - состояние процесса, в котором k-я сменно-запасная часть находится в запасе; y0 - состояние процесса, в котором k-я сменно-запасная часть использована; t - текущее время рассматриваемого периода эксплуатации; λ - интенсивность выхода изделия из эксплуатации в результате отказа; μ - интенсивность ввода изделий в эксплуатацию в результате ремонта; γ - интенсивность пополнения запаса сменно-запасных частей в результате их ремонта, производства или заказа; ν - интенсивность расходования сменно-запасных частей из запаса; - среднее число работоспособных однотипных изделий; - среднее число сменно-запасных частей, находящихся в запасе. γ = φ + ξ + ψ где φ - интенсивность поставок; ξ - интенсивность ремонта; ψ - интенсивность производства. Опираясь на графы, определены функции приращения количества изделий в эксплуатации и количества сменно-запасных частей в запасе от времени. Через определение матожиданий этих приращений получена система дифференциальных уравнений, которую можно назвать моделью процесса обеспечения работоспособности изделий (обеспечения эксплуатационной готовности судов сменно-запасными частями) на время эксплуатации со стационарным потоком событий: В более общей форме, учитывающей количество судов, число изделий на них и значение параметров, модель приобретает следующий вид: где i - номер судна; h - количество судов в рассматриваемой группе; j - номер единицы изделия на судне; ni - число единиц изделий на i-м судне; I - номер судна группы, обеспечиваемой групповым запасом сменно-запасных частей; H - число судов в группе, обеспечиваемой групповым запасом сменно-запасных частей; J - номер единицы изделия на судне группы, обеспечиваемой групповым запасом сменно-запасных частей; Ni - число единиц изделий на I-м судне группы, обеспечиваемой групповым запасом сменно-запасных частей. Для упрощения записи представим как , а как . Для учёта нестационарности потоков событий рассматриваемый период эксплуатации можно представить как бесконечное количество бесконечно малых периодов времени со стационарными потоками событий. Тогда при , , описываемая модель представляется как система уравнений (1) где A - закон изменения интенсивности выхода из строя изделий рассматриваемой группы судов Λ во времени; B - закон изменения интенсивности восстановления готовности изделий к использованию μ во времени; C - закон изменения интенсивности пополнения запаса ЗИП γ от t. В условиях реальной эксплуатации A может иметь зависимость также от , и ; B от , , ; C от , , . Подтверждением корректности разработанной модели служит её успешное применение для решения практических задач. Рассмотрим ситуацию, когда для рыбопромыслового флота из h судов необходимо смоделировать процесс эксплуатации с заданными параметрами и скорректировать их для обеспечения достижения целевой функцией необходимого значения. В этом конкретном случае скорректировать необходимо начальное количество деталей редукторов для компрессоров в запасе my0 (работой этих изделий в большой степени определяется эффективность данной группы судов), а также значения μ и γ таким образом, чтобы объём выполненной флотом работы, пропорциональный времени использования судов, оказался не ниже требуемого Vmin. На рассматривавшихся судах установлены потребляющие сменно-запасные части изделия в таких количествах ni, что общее их число составило mx0. Групповой пополняемый ЗИП (ЗИП-Г) доступен для всех судов группы. На начало периода эксплуатации в ЗИП-Г содержалось my0 сменно-запасных частей, включая изделия резерва (холодного). Необходимо отметить, что если часть изделий работает в горячем резерве, то в период времени, когда основное изделие неисправно, интенсивность отказов резерва возрастает вместе с нагрузкой на него. При этом нарушается допущение об отсутствии последействия потока событий. Описать зависимость интенсивности отказов от времени можно используя аппарат теории Маркова, определив стратегию вывода из эксплуатации судов в случае нехватки сменно-запасных частей (рис. 2). Рис. 2. Граф процесса поддержания работоспособности изделий: i0 - состояние, в котором все изделия, имеющие горячее резервирование, находятся в эксплуатации; i1 - состояние, в котором первое изделие, имеющее горячий резерв, выведено из эксплуатации; in - состояние, в котором все изделия, имеющие горячий резерв, выведены из эксплуатации Затем, решая полученную с использованием графа систему уравнений Колмогорова, находятся вероятности пребывания системы в каждом из n + 1 состояний в каждый малый период времени, после чего можно получить приращения mx от времени. В данном случае группа судов - это единственный потребитель сменно-запасных частей для ЗИП-Г, следовательно Вследствие того, что при каждом выходе из строя изделия требуются сменно-запасные части из ЗИП-Г, для поддержания всех изделий в работоспособном состоянии существует жёсткая зависимость: . Когда сменно-запасные части на складе отсутствуют и интенсивность ремонта и изготовления частей не позволяет создать запас, для системы справедливо выражение при . Система уравнений (1) решается поэтапно методом приближенных вычислений (методом Рунге - Кутта). Время эксплуатации рассматриваемой группы судов для этого было условно разделено на шесть этапов, характеристики которых приведены в табл. 1. Таблица 1 Значения переменных, характеризующие периоды эксплуатации рассматривавшейся группы судов Период времени Переменная t0 ÷ t1 t1 ÷ t2 t2 ÷ t3 t3 ÷ t4 t4 ÷ t5 t5 ÷ t6 μ, ч-1 = Λ, const < Λ, const (предел) < Λ, const (предел) < Λ = γ, (потенциально может быть выше) > Λ, const > Λ, const γ, ч-1 const (0) const (0) < Λ, const < Λ, K/mx (↑) > Λ, const > Λ, const , ч-1 Норма Выше нормы Выше нормы Выше нормы Выше нормы Норма my, ед. > 0, ↓ > 0, ↓ > 0, ↓ const (0) > 0, ↑ > 0, ↑ mx, ед. const (mx0) < mx0, ↓ < mx0, ↓ < mx0, ↓ < mx0, ↑ < mx0, ↑ Примечание. ↑ - значение переменной возрастает; ↓ - значение переменной убывает. В период времени t0 ÷ t1 на складе имеется запас сменно-запасных частей, интенсивность ввода изделий в эксплуатацию μ обеспечивает работу всех изделий. Общая интенсивность отказов изделий группы - норма. Пополнение запаса не производится (γ = 0). Все вышедшие из строя изделия без ожидания заменяются исправными. В момент времени t1 = 700 ч происходит увеличение интенсивности эксплуатации судов, вследствие чего интенсивность отказов изделий Λ кратно возрастает (период времени t1 ÷ t2). Интенсивность ввода в эксплуатацию изделий μ меньше интенсивности их вывода Λ. Суда выводятся из эксплуатации в зависимости от используемого критерия оценки эффекта от их использования. В случае, когда количество сменно-запасных частей в запасе становится меньше 6 шт., начинается пополнение запаса с постоянной интенсивностью (γ > 0). Этому периоду соответствует интервал времени t2 ÷ t3. Когда запас сменно-запасных частей на складе заканчивается (my = 0), интенсивность ввода изделий в эксплуатацию μ становится равной интенсивности пополнения сменно-запасных частей γ (интервал t3 ÷ t4). Пополнение запаса сменно-запасных частей производится с интенсивностью, пропорциональной mx. По прошествии времени T2 = 3 500 ч от начала периода эксплуатации интенсивность пополнения сменно-запасных частей γ и интенсивность ввода изделий в эксплуатацию μ возрастают до фиксированных значений, превышающих значение интенсивности вывода Λ (интервал t4 ÷ t5). Количество судов в эксплуатации увеличивается, на складе появляется запас сменно-запасных частей (my > 0). При T3 = 8 000 ч интенсивность эксплуатации судов снижается до начального значения, следовательно, и величина интенсивности отказов изделий Λ снижается до начального значения (интервал t5 ÷ t6). Конец интервала - момент, когда все изделия становятся работоспособными. Построенный по результатам моделирования график изменений во времени количества эксплуатирующихся изделий mx(t) и количества сменно-запасных частей в запасе my(t) показан на рис. 3. Рис. 3. График изменения во времени mx(t) и my(t): -- - mx(t); ……. - my(t) Искомый объём выполненной работы при этом выражается формулой (2) где T - период времени функционирования флота; k(t) в подынтегральном выражении - средний объём работы, обеспечиваемый одним изделием на судах, находящихся в эксплуатации в момент времени t. Используя формулу (2), находим за t = [0; T3], вычислен момент, когда все изделия вновь становятся работоспособными (t6 ≈ 10 650 ч). В данном случае учтена возможность перестановки изделий с одного судна на другое. В случае нехватки сменно-запасных частей может складываться ситуация, когда судно вводится в эксплуатацию периодами при появлении сменно-запасных частей. В этом случае необходимо учитывать, что с судна сменно-запасные части не снимаются и одновременно с необходимостью решения системы (1) появляется необходимость решать задачу оптимизации состава флота при эксплуатации в условиях недостатка сменно-запасных частей. Для выявления наиболее эффективных способов увеличения V определён положительный вклад αi+ каждой из N переменных () в целевую функцию V для каждого из этапов эксплуатации рассматриваемой группы судов. Для этого использовалось выражение Полученные значения вкладов приведены в табл. 2. Таблица 2 Значения положительных вкладов переменных в целевую функцию Период времени Переменная t0 ÷ t1 t1 ÷ t2 t2 ÷ t3 t3 ÷ t4 t4 ÷ t5 t5 ÷ t6 μ 0 1,218·106 T ≈ 10 650 1,224·106 T ≈ 10 650 0 1,257·106 T ≈ 9 550 - T ≈ 9 150 γ 1,306·106 T ≈ 10 550 1,306·106 T ≈ 10 550 1,306·106 T ≈ 10 550 1,306·106 T ≈ 10 550 0 0 my0 1,306·106 T = 10550 μ и γ 1,174·106 T ≈ 10 650 1,174·106 T ≈ 10 650 1,198·106 T ≈ 10 650 1,563·106 T ≈ 10 350 1,58·106 T ≈ 3 500 - T ≈ 8 000 Зная дисперсии величин μ, γ, Λ и используя свойства среднеквадратических отклонений, были определены распределения количества , во времени и взаимообусловленные с ними распределения длительности интервалов времени t1 ÷ t2; t2 ÷ t3; t3 ÷ t4; t5 ÷ t6 для данной группы судов, после чего вычислены значения среднеквадратических отклонений величин V и t6. Используя полученные значения, вероятность выполнения флотом необходимого объёма работы можно определить по формуле . Результатом расчёта обоснованы значения страховых коэффициентов для получения необходимых значений my0, μ, γ, позволивших получить необходимый эффект с требуемой вероятностью и скорректировать вероятность достаточности ЗИП для рассматриваемого флота. Заключение Применение разработанной модели и подхода для принятия решений по изменению объёма начального запаса сменно-запасных частей, созданию дополнительных каналов их поставки, ремонтных и производственных мощностей, ресурсов, выводов об эффективности использования флота в различных условиях позволили обеспечить эксплуатационную готовность групп судов. Необходимо отметить, что для получения достаточной точности результатов использование предложенного подхода требует существенных вычислительных и временны́х затрат при выполнении расчётов статистического моделирования. Поэтапность решения, ведущая к накоплению погрешностей при переходах в расчётах от одного этапа к другому, также говорит о необходимости использования при статистическом моделировании больших выборок случайных величин. Видится возможным снятие некоторых ограничений на использование модели, наложенных принятыми при её разработке допущениями. Данные шаги в перспективе позволят сделать круг решаемых с её помощью задач ещё шире.
Список литературы

1. Барзилович Е. Ю. Модели технического обслуживания сложных систем: учеб. пособие / Е. Ю. Барзилович. М.: Высш. шк., 1982. 231 с.

2. Башуров Б. П. О стратегии технического обслуживания и ремонта вспомогательного оборудования энергетических установок судовых транспортных средств / Б. П. Башуров, Е. С. Носенко, В. В. Шарик // Машиностроение. 2004. № 10. С. 29-35.

3. Трунин С. Ф. Надежность судовых машин и механизмов: учеб. пособие / С. Ф. Трунин, Л. А. Промыслов, О. Р. Смирнов. Л.: Судостроение, 1980. 192 с.

4. Голинкевич Т. А. Прикладная теория надежности: учеб. для вузов / Т. А. Голинкевич. М. Высш. шк., 1977. 160 с.

5. Смелякова Н. Н. Техническое обслуживание машин, оборудования и приборов зарубежными фирмами / Н. Н. Смолякова. М.: ВО «Внешторгреклама», 1978. 408 с.

6. Рудницкий А. В. Обеспечение надёжности судовых систем и оборудования на этапах проектирования и эксплуатации с применением информационных технологий / А. В. Рудницкий, В. С. Виноградов // Морской вестник. 2012. № 3 (43). С. 45-49.


Войти или Создать
* Забыли пароль?