ANALYSIS OF PROVISION OF OPERATIONAL AVAILABILITY OF THE GROUP OF THE VESSELS WITH SPARE PARTS USING PROBABILISTIC MODELING
Authors:
1.
Scientific Industrial Center "ANSTEKO"
2.
Volga Academy of Water Transport
3.
Nizhniy Novgorod State Technical University named after R. E. Alekseev
4.
Central Designing Bureau "Lazurit"
Type:
Article
Pages:
from 51
to 58
Status:
Published
Received:
20.02.2015
Accepted:
25.02.2015
Published:
25.02.2015
Subject area:
UDK 629.5.07
GRNTI 45.01 Общие вопросы электротехники
GRNTI 55.42 Двигателестроение
GRNTI 55.45 Судостроение
GRNTI 73.34 Водный транспорт
GRNTI 44.31 Теплоэнергетика. Теплотехника
GRNTI 45.01 Общие вопросы электротехники
GRNTI 55.42 Двигателестроение
GRNTI 55.45 Судостроение
GRNTI 73.34 Водный транспорт
GRNTI 44.31 Теплоэнергетика. Теплотехника
Language:
Russian
Keywords:
vessels, operational availability, spare parts, model, probabilistic approach
Abstract (English):
The paper presents the model, designed to ensure the operational availability of single-type components of the group of vessels. Its usage aids to adjust the volume of the initial stock of spare parts in the fleet bases, repair and production facilities, ordering of the spare parts for a group stock. The process of operation of single-type components of the group of ships is supposed to be presented as a set of the processes of consumption and replenishment of the spare parts and the processes of exception and subsequent introduction of the items into operation. The consistent recurrent solution of the developed system of equations (model) was applied to each of the conditional stages of the considered operational period using Runge - Kutt method. Accounting the distributions of the deviations from their mean value is implemented using statistical modeling. The successful application of this approach to the issue of providing availability of the particular group of vessels is demonstrated. Using the development contributes to solving the mentioned issues in the organizations engaged in the fleet operation.
The paper presents the model, designed to ensure the operational availability of single-type components of the group of vessels. Its usage aids to adjust the volume of the initial stock of spare parts in the fleet bases, repair and production facilities, ordering of the spare parts for a group stock. The process of operation of single-type components of the group of ships is supposed to be presented as a set of the processes of consumption and replenishment of the spare parts and the processes of exception and subsequent introduction of the items into operation. The consistent recurrent solution of the developed system of equations (model) was applied to each of the conditional stages of the considered operational period using Runge - Kutt method. Accounting the distributions of the deviations from their mean value is implemented using statistical modeling. The successful application of this approach to the issue of providing availability of the particular group of vessels is demonstrated. Using the development contributes to solving the mentioned issues in the organizations engaged in the fleet operation.
Keywords:
vessels, operational availability, spare parts, model, probabilistic approach
vessels, operational availability, spare parts, model, probabilistic approach
Введение Для периода эксплуатации судовых систем и оборудования существует проблема определения количества группового запаса сменно-запасных частей в местах базирования флота и необходимой интенсивности его пополнения, что осложняется особенностями эксплуатации судов, стратегий восстановлений и ремонтов. Запас должен перекрывать потребности флота в сменно-запасных частях, обеспечивая его готовность к использованию, и одновременно с этим быть минимальным из-за необходимости сокращения издержек на его содержание. Цель исследования - обеспечение эксплуатационной готовности группы судов путём определения необходимого на период эксплуатации начального запаса сменно-запасных частей, а также значений величин, характеризующих процессы восстановления изделий и пополнения запаса. Публикации по данной проблеме [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. Barzilovich E. Yu. Modeli tehnicheskogo obsluzhivaniya slozhnyh sistem: ucheb. posobie / E. Yu. Barzilovich. M.: Vyssh. shk., 1982. 231 s.
2. Bashurov B. P. O strategii tehnicheskogo obsluzhivaniya i remonta vspomogatel'nogo oborudovaniya energeticheskih ustanovok sudovyh transportnyh sredstv / B. P. Bashurov, E. S. Nosenko, V. V. Sharik // Mashinostroenie. 2004. № 10. S. 29-35.
3. Trunin S. F. Nadezhnost' sudovyh mashin i mehanizmov: ucheb. posobie / S. F. Trunin, L. A. Promyslov, O. R. Smirnov. L.: Sudostroenie, 1980. 192 s.
4. Golinkevich T. A. Prikladnaya teoriya nadezhnosti: ucheb. dlya vuzov / T. A. Golinkevich. M. Vyssh. shk., 1977. 160 s.
5. Smelyakova N. N. Tehnicheskoe obsluzhivanie mashin, oborudovaniya i priborov zarubezhnymi firmami / N. N. Smolyakova. M.: VO «Vneshtorgreklama», 1978. 408 s.
6. Rudnickiy A. V. Obespechenie nadezhnosti sudovyh sistem i oborudovaniya na etapah proektirovaniya i ekspluatacii s primeneniem informacionnyh tehnologiy / A. V. Rudnickiy, V. S. Vinogradov // Morskoy vestnik. 2012. № 3 (43). S. 45-49.
