CC BY
9
УДК510.644.4: 656.71 В01: 10.18287/2541-7533-2021-20-3-160-170
ОПТИМИЗАЦИЯ РАСПРЕДЕЛЕНИЯ ТЕХНОЛОГИЧЕСКИХ РЕСУРСОВ АВИАТРАНСПОРТНОЙ СИСТЕМЫ ПРИ НАЛИЧИИ РАСПИСАНИЯ И НЕЧЁТКИХ ИСХОДНЫХ ДАННЫХ
© 2021
В. А. Романенко кандидат технических наук, доцент кафедры организации и управления перевозками на транспорте; Самарский национальный исследовательский университет имени академика С.П. Королёва; у!а [email protected]
Рассматривается задача оптимального распределения технологических ресурсов (операторов) технической или организационно-технической системы, предназначенных для обслуживания определённых объектов (операндов) по заданному расписанию. Учитываются необходимость включения в состав процесса обслуживания, совместно с основной, подготовительной и заключительной операций, возможность выделения для обслуживания операнда одного или нескольких операторов, а также зависимость продолжительности операций от факторов, характеризующихся неопределённостью. В связи с предполагаемым отсутствием статистики используются экспертно-заданные неопределённые величины в форме треугольных нечётких чисел. Оптимизационная задача формулируется как задача математического программирования с нечётким критерием и «чёткими» ограничениями, состоящая в поиске такого распределения заданного числа операторов для обслуживания каждого операнда из заданной совокупности, которое минимизирует целевую функцию, учитывающую отклонения от расписания (опоздания) с окончанием обслуживания. Характерными примерами систем, для которых задача актуальна, служат производственные комплексы авиатранспортных предприятий, функционирующие в условиях неопределённости при необходимости обеспечения регулярности и безопасности воздушных перевозок. Приведён модельный пример решения задачи распределения подвижных средств топливозаправки в узловом аэропорту с учётом особенностей его расписания. Показано, что для решения достаточно возможностей стандартного программного обеспечения персональной вычислительной техники.
Оптимизация; технологические ресурсы; операнд; треугольное нечёткое число; нечёткое математическое программирование; узловой аэропорт
Цитирование: Романенко В.А. Оптимизация распределения технологических ресурсов авиатранспортной системы при наличии расписания и нечётких исходных данных // Вестник Самарского университета. Аэрокосмическая техника, технологии и машиностроение. 2021. Т. 20, № 3. С. 160-170. DOI: 10.18287/2541-7533-2021-20-3-160-170
Введение
Объектом исследования является техническая или организационно-техническая система, располагающая комплектом операторов - технологических средств, предназначенных для осуществления определённого технологического процесса по целенаправленному преобразованию (обработке или обслуживанию) некоторых объектов -операндов. Предполагается, что обслуживание операндов должно производиться в соответствии с заданным расписанием. Для обслуживания операнда могут привлекаться как один, так и несколько операторов, величины производительности которых могут различаться в зависимости как от их типа, так и от типа преобразуемого операнда. Операнды характеризуются объёмами работ, необходимыми для их успешного обслуживания, и максимальной численностью операторов, которые могут быть одновременно использованы для обслуживания операнда. Процесс обслуживания операнда включает подготовительную, основную и заключительную операции. Один или несколько операторов, выделенных для обслуживания операнда, предполагаются занятыми в течение всего промежутка времени от момента начала подготовительной операции до момента окончания заключительной. Фактические продолжительности всех трёх операций
априори достоверно неизвестны. Главным фактором, влияющим на вероятностные распределения времени подготовительной и заключительной операций, предполагается тип оператора. Распределение времени основной операции зависит в первую очередь от потребного для обслуживания операнда объёма работ, численности используемых с этой целью операторов и интенсивности выполнения работ операторами. Перечисленные параметры также могут характеризоваться неопределённостью.
Характерными примерами рассматриваемой системы служат производственные комплексы предприятий воздушного транспорта, в частности топливозаправочный комплекс аэропорта, оснащённый монтируемыми на шасси автомобиля передвижными заправочными агрегатами (ПЗА), предназначенными для заправки воздушных судов (ВС) топливом с помощью систем централизованной заправки (ЦЗС). Перед указанными системами, функционирующими в условиях неопределённости, с особой остротой стоит задача обеспечения пунктуального выполнения расписания, так как его нарушение может служить прямой причиной нарушения регулярности и безопасности воздушных перевозок. Примеры решения задач оптимального планирования и распределения технологических ресурсов аэропортовых комплексов в весьма сложной и детальной постановке рассмотрены в работах [1; 2], в которых, однако, не учитывается неопределённый характер исходных данных. С другой стороны, к настоящему времени накоплен основательный методологический задел в области сетевого планирования и управления проектами с учётом влияния неопределённости различной природы, изложенный, например, в работах [3; 4], в которых отсутствует учёт специфики авиатранспортных систем. Ниже предлагается методика оптимизации распределения технологических ресурсов, в первую очередь систем наземного обслуживания перевозок, которая может быть реализована на базе стандартного программного обеспечения (например, в условиях аэропортового предприятия), учитывающая неопределённость исходных данных.
В соответствии с практикой предполагается, что статистика, позволяющая сформировать частотные распределения, необходимые для формулирования задачи оптимального распределения операторов в терминах теории вероятностей, отсутствует, а её накопление невозможно или затруднено. В этом случае все неопределённые величины определяются методом экспертного оценивания и задаются в форме нечётких величин, благодаря чему рассматриваемая оптимизационная задача сводится к разряду задач нечёткого математического программирования [5; 6].
Математическая модель
Остановимся на используемых ниже понятиях и приемах нечёткой арифметики. Под нечётким множеством А на универсальном множестве и понимается совокупность кортежей вида /а (и), и), где / (и) - степень принадлежности элемента и е и
нечёткому множеству А, которая задаётся как действительное число из интервала [0,1]. Функция, позволяющая для произвольного элемента универсального множества вычислить степень его принадлежности нечёткому множеству, называется функцией принадлежности. Под нечёткой величиной (НВ) понимается произвольное нечёткое множество, заданное на множестве действительных чисел [7]. Одним из наиболее часто используемых, простых и удобных для практических вычислений типов НВ является НВ с «треугольной» функцией принадлежности или «треугольное нечёткое число», характеризующее неопределённость типа «приблизительно равно». Треугольное нечёткое число (ТНЧ) А представим в виде кортежа А = ( ьа, м а, , включающего координаты опорных точек функции принадлежности / (и):
Я
Xй )=
u - La
Ma - La' Ra - u
Ra -M a'
a < u <Ma,
Ma < u<Ra,
0, u <La, Ma < u
где M a - наиболее возможное значение u (мода), L a, R a - соответственно наименьшее и наибольшее из возможных значений u, La<Ma<Ra.
Основываясь на принципе обобщения Заде [7], определим необходимые ниже
операции с ТНЧ. Пусть даны два ТНЧ: A1 = (La1,Ma1,RaA2 = (La2,Ma2,Ra2). Результатами их нечёткого сложения (+), перемножения (х), вычитания (-), деления (/), взятия максимума (max) и минимума (min) будут ТНЧ, определяемые соответственно как:
А1 + (х)А2 = (ьа1 + (х) 1а2, ма1 + (х)ма2, Яа1 + (х)Яа2),
А - (/)А~2 = ( Ч - (/)ЙС2, мС1 - (/)мС2, % - (/)ьа2), (1)
тах(тт)А, А2)= (тах^т^Ч ,ьа2), max(min)(мa1 ,ма2), тах^т)^^,Яа2)).
Для дефаззификации, то есть приведения к чёткой форме, ТНЧ используем метод центра тяжести [7], в соответствии с которым чёткое число аа - результат дефаззификации ТНЧ А = (ьа,ма,Яа) - определяется как
da =
,si~A La + M a + Ra def (A) =---. (2)
Будем считать, что рассматриваемая система располагает п операторами, которые предполагается задействовать для обслуживания к операндов. Обозначим г, j индексы оператора и операнда, соответственно г = 1,...,п, j = 1,...,к . Введём следующие обозначения величин, составляющих исходные данные задачи:
nj - максимальное число операторов, которые могут одновременно принимать участие в обслуживании операндаj, j = 1,..., к;
Qj - объём работ при обслуживании операндаj, j = 1,...,к;
Яг - максимальная производительность оператора г, г = 1,...,п ;
Я}- - максимально допустимый уровень производительности при выполнении операции по обслуживанию операндаj, j = 1,...,к;
^, Atf - плановые (заданные в соответствии с расписанием), соответственно момент времени начала и продолжительность основной операции обслуживания операнда j, j = 1,...,к;
At1П, At13 - нормативные (заданные в соответствии с нормативной картой) продолжительности соответственно подготовительной и заключительной операций, выполняемых г-м оператором, г = 1,...,п;
АТгп, АТг3 - фактическая продолжительность соответственно подготовительной и заключительной операций, выполняемых г-м ресурсом, г = 1,...,п.
Все входные величины являются либо ТНЧ, выражающими мнение эксперта, либо обычными числами и, таким образом, все полученные в результате вычислений величины также будут либо числами, либо ТНЧ.
Будем считать, что нумерация операндов упорядочена в порядке возрастания плановых моментов времени начала основной операции их обслуживания:
г0 = 0<^ <гу, у = 1,...,к,
где г0 = 0 - начальный момент модельного времени г, в рамках которого определяются временные характеристики технологического процесса.
Пусть для обслуживания операндау должно быть выделено не менее одного и не более пу операторов. Выделенные для обслуживания операнда у операторы одновременно начинают и одновременно завершают основную операцию.
Введём бинарную переменную Ху е (0,1}, принимающую значение 1, если оператор г выделяется для обслуживания операнда у, и значение 0 в противном случае, г = 1,..., п, у = 1,..., к. С учётом введённых допущений и правил нечёткой арифметики (1)
НВ А ТО фактической продолжительности выполнения основной операции обслуживания операнда у будем рассчитывать согласно выражению
А~0 =-%~т, У = 1,.., к, (3)
I Я
г=1 : Ху =1
где - производительность оператора г при обслуживании операнда у, ограниченная величинами Яг и Яу, определяемая посредством взятия нечёткого минимума (1):
Яу = шт(д,Я), г = 1,...,п, у = 1,...,к.
Учитывая неопределённый характер фактической продолжительности АТ° выполнения основной операции обслуживания операнда у, следует ожидать, что и время ТН фактического начала этой операции, и время ТуК её фактического окончания будут неопределёнными и, следовательно, будут отличаться: первое - от планового времени ¿у её начала, второе - от планового времени гК её окончания. ТНЧ ТуК и обычное число гК будем определять соответственно как
ТК = ТН +ато , гК = гу + Аг0, у = 1,..., к. (4)
Задача состоит в определении такого множества (ху, г = 1,...,п, у = 1,...,к}, которое минимизировало бы некоторую величину, связанную с промежутками времени опоздания с завершением выполнения основных операций обслуживания операндов. Учитывая, что в качестве опоздания следует рассматривать только такую ситуацию, когда фактическое время завершения основной операции превышает плановое время её за-
вершения, фактическое опоздание ATf с завершением основной операции определим как ТНЧ с использованием нечёткой операции взятия максимума (1):
ДГК = max(0, (j - tf)) j = 1,...,k. (5)
Как следует из (4) и (5), определение ATf возможно, если известно Tf, для установления которого предварительно должен быть определён ряд временных характеристик технологического процесса. С этой целью примем, что обслуживание операнда выделенным для этого оператором может начаться в том случае, если оператор готов к началу обслуживания. Состояние готовности оператора к началу выполнения подготовительной операции обслуживания операнда наступает в момент, когда рассматриваемый оператор завершает заключительную операцию обслуживания предыдущего операнда. Если до начала обслуживания рассматриваемого операнда оператор не был задействован в обслуживании других операндов, то моментом времени готовности такого оператора примем начальный момент времени t0 = 0 . Таким образом, фактическое
время наступления готовности Tf оператора i к началу обслуживания операнда j
(j £ {l,...,k}, Xj = 1) принимается равным либо 10 = 0, если операнд является первым по
порядку среди операндов, к обслуживанию которых назначен оператор i, либо моменту освобождения оператора от обслуживания предыдущего операнда:
~ lo j = m
Tf =l~ ю . = T i = 1,..., n, Л )£{l,..., k}, (6)
l1 i j(i) > J ~ Jm>1j
где Tf3 - фактическое время завершения конечной операции процесса обслуживания i-м оператором j-го операнда; j() - m-й элемент множества {j(l)}е{1,...,k} индексов операндов, для обслуживания которых назначен i-й оператор, j) е {j()}, 0 < m < k .
С использованием переменной xij выражения (3) и (6) могут быть представлены в форме, более удобной для программной реализации табличных вычислений:
T о = Qi
AT? , j = 1,..., k.
%'Xij
J T S Ri>x-
i=1
Величина 17 для любого г е{1,...,п\ и такого 7 е {1,..., к}, для которого выполняется условие Ху = 1, определяется с использованием рекуррентной формулы:
Tf =
T j*=1
j, xi,/-1=0, j*=j, i=n. j x- ,*-1 =1, j* = 2,.. i,
Г> ГТ1 Г • * 1 • * ГК • *
Здесь величины 17 определяются последовательно для 7 = 1, 7 = 2, ..., 7 = 7. Найденная таким образом величина 17* для 7 * = 7 и будет искомой Т7 .
Взаимосвязь временных характеристик представлена на рис. 1, где в качестве примера рассмотрен процесс обслуживания некоторого (/-го) операнда двумя выделенными с этой целью операторами с индексами I¡" е {1,..., п}, для которых справедливо равенство х,",/ = х,",/ = 1 (отметим, что для I е{1,..., п}\ {" , I"} здесь х^ = 0). Тонировкой
отмечен график функции принадлежности искомой НВ АТК на универсальном множестве А/, для оси абсцисс которого в качестве начала отсчёта А/ = 0 принимается чёткое значение /. В отличие от АТК функции принадлежности всех остальных НВ, приводимых на рис. 1, строятся на множестве
; ¡\ г
'1 Т-Г- I \ Т^1 / / /т ////\ / 1 / / \ / , с /У // // // / / г ~ V т^я \ ' * тК/ тч / / / / / / !() 1 * г \ г* / V \ \ \ \ N -—► Д/ оиг/зднс* 7 / \ уЮ I / \ „
-к 7нп А '/ \ -г А77 1 А ; / : / : / ; 1 V/ А / \ т \ ТЯ \т> \ \ \ \ \ \ \ у г/У / / / / / / \ \ \ \ \ \ N оператор V / \ 7 Р 1 / \ .
оператор V
Рис. 1. Временные характеристики процесса обслуживания операнда/
Запишем выражения для нечётких расчётов других временных характеристик технологического процесса обслуживания/-го операнда, предполагая ТНЧ ТГ подсчитанным. Поскольку подготовительная операция оператора I к обслуживанию операнда/ не может начаться раньше планового времени / и фактического времени готовности
Т/, то фактическое время начала этой операции определяется как ТНЧ:
/ = шах/,Т/Г), I = 1,...,п, / = 1,...,к : ху = 1. (7)
Входящее в (7) плановое время / начала подготовительной операции определяется по заданным плановым величинам / и :
/ = /-А^П, / = 1,..., п, / = 1,..., к : х/ = 1.
Фактическое время ТК окончания подготовительной операции определяется в результате нечёткого сложения ТНЧ Т/НП, подсчитанного по (7), и заданного АТП :
Т/КП = ТНП + АТ,П, I = 1,..., п, / = 1,..., к : х/ = 1.
Фактическое время Т/ начала основной операции обслуживания у-го операнда выделенными операторами, т.е. имеющими такие индексы {/}с{1,..., п}, для которых выполняется условие Ху = 1, определяется как наибольшая величина в результате сравнения планового срока ' начала указанной операции с наиболее поздним сроком окончания подготовительной операции оператора из числа выделенных для обслуживания рассматриваемого операнда. Указанная процедура осуществляется посредством последовательного взятия нечётких максимумов (1):
I',
Ту = тах1 тах ', у =1,-,к.
{1,...,п}: Ху =1
Величины фактического времени окончания основной операции дляу-го операнда Т}К и выполняющих операторов ТК совпадают и определяются как ТНЧ:
ТК - ТК = ТН + дТ°, г = 1,...,п, у = 1,...,к : Ху = 1.
Наконец, НВ фактического времени окончания заключительной операции обслуживания у-го операнда выделенными операторами определяется как ТНЧ в предположении, что заключительная операция выполняется сразу после основной:
ТК = ТуК +Д~З, г = 1,...,п, у = 1,...,к : Ху = 1.
Приведённые формулы позволяют по заданным {ху, г = 1,...,п, у = 1,...,к}, , Яу, у, ду (у = 1,...,к) и Яг (г = 1,...,п) оценить фактические опоздания ДТ^ (у = 1,...,к).
Логично полагать, что опоздания с завершением обслуживания операндов, требующих выполнения большего объёма работ, приводят к большим издержкам в том или ином смысле. Следовательно операторы должны быть так распределены для обслуживания операндов, чтобы обеспечить минимум показателю, учитывающему не только величины времени опоздания, но и соответствующие операндам объёмы работ. Исходя из представленных соображений, выберем в качестве критерия рассматриваемой оптимизационной задачи минимум следующей нечёткой суммы:
Т = £дТК0;->тт. (8)
у=1
Проектными являются бинарные переменные Ху, г = 1,...,п, у = 1,...,к . Ограничения накладываются на число операторов для обслуживания одного операнда:
1 <]ГХу < п;, у = 1,...,к . (9)
г=1
Поскольку оптимизационная задача с чёткими ограничениями (9) и нечётким критерием (8) в силу нечёткости последнего является задачей с бесконечным числом критериев [5], то формулировку (8), (9) следует рассматривать лишь как условную. Для придания постановке задачи определённости используем в качестве целевой функции
значение dS = def [~j, полученное в результате дефаззификации методом центра тяжести (2) ТНЧ S . Таким образом, задача оптимального распределения технологических ресурсов сводится к определению множества {xij е {0,l}, i = 1,...,n, j = 1,...,k}, обеспечивающего минимум целевой функции dS и удовлетворяющего ограничениям (9) при заданных Qj, Rj, tj, Atf ( j = 1,...,k) и Ri (i = 1,...,n). Решение задачи достигается известными «чёткими» методами математического программирования, в том числе используемыми программной надстройкой «Поиск решения» табличного процессора Microsoft Excel.
Модельный пример решения оптимизационной задачи
Оптимизационная задача решена на модельном примере с исходными данными топливозаправочного комплекса регионального узлового аэропорта, оснащённого ЦЗС. Предполагается, что аэропорт располагает семью ПЗА, из которых три отличаются большей производительностью. Характеристики ПЗА следующие:
- Яг = (0.9,1.0,1.1 м3/мин, АТ3 = (8.0,10.0,12.0) мин, АТгП = (4.0, 5.0, 6.0) мин, А^П = 5.25 мин, Аtг3 = 10.5 мин, если г е {1, 2, 3};
- Д = (0.45, 0.50, 0.55)м3/мин, А~3 = (6.4, 8.0, 9.6) мин, А~П = (3.2,4.0,4.8) мин, АггП = 4.2 мин, Аt13 = 8.4мин, если г е {4,..., 7}.
Модельное расписание поступления ВС на заправку, приведённое в табл. 1, отражает специфические черты потоков ВС в узловом аэропорту, состоящие в резких и значительных периодических колебаниях их интенсивностей [8]. Моделируется приблизительно двухчасовой промежуток времени массового поступления ВС на заправку, первая половина которого характеризуется ростом интенсивности потока ВС, а вторая половина - её снижением. Каждое из заправляемых ВС относится к одной из четырёх групп, различающихся характеристиками топливной системы и, как следствие, параметрами процесса заправки, сведёнными в табл. 2.
Таблица 1. Модельное расписание поступления ВС на заправку в региональном узловом аэропорту
j 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
tj, мин 10 20 30 40 45 50 55 60 65 70 75 80 90 100 110
Группа ВС I II II III III IV I II II III III IV II III III
Таблица 2. Характеристики ВС
Группа ВС I II III IV
j 1, 7 2, 3, 8, 9, 13 4, 5, 10, 11, 14, 15 6, 12
nj 2 2 1 1
Atf, мин 40 40 30 30
Qj, мин <31.5, 35.0, 38.5) <27.0, 30.0, 33.0) <9.0, 10.0, 11.0) <7.2, 8.0, 8.8)
Rj, м3/мин <0.95, 1.0, 1.05) <0.95, 1.0, 1.05) <0.48, 0.5, 0.52) <0.28, 0.3, 0.31)
Результаты оптимизации, представленные в табл. 3, демонстрируют возможность заметных отклонений от плана с завершением обслуживания ВС, однако обязательно задержка реализуется только для 10-го ВС. Каждому ПЗА предстоит заправить от 2 до 4 ВС, для большинства ПЗА это число вполне предсказуемо составляет 3 ВС.
Таблица 3. Результаты оптимизации
1 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Распределение ПЗА
1 1 1 0 1 0 0 0 0 0 1 0 0 0 0 1 0
2 0 0 0 1 0 0 0 1 0 0 0 0 1 0 0
3 1 0 0 0 0 0 1 0 0 1 0 0 0 0 0
4 0 1 0 0 0 1 0 0 0 0 0 1 0 0 0
5 0 1 0 0 0 0 0 1 0 0 0 0 1 0 0
6 0 0 0 0 1 0 0 0 0 0 1 0 0 0 1
7 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0
Временные характеристики процесса заправки ВС
~ Н 1 1 5 <10.0, 10.0, <20.0, 20.0, <37.0, 42.5, <40.0, 40.0, <45.0, 45.0, <54.1, 62.0, <55.0, 55.0, <69.1, 75.0, <65.9, 77.5, <97.0, 105.0, <75.0, 77.0, <86.6, 100.7, <98.0, 110.0, <100.0, 112.5, <110.0, 110.0,
мин 10.8) 20.6) 50.1) 40.8) 45.6) 71.7) 55.8) 89.1) 92.6) 116.5) 84.8) 116.9) 124.4) 135.0) 123.3)
~ К 1 1 5 <25.0, 27.5, <44.5, 50.0, <53.9, 62.5, <57.1, 60.0, <62.1, 65.0, <77.0, 88.7, <85.0, 90.0, <86.0, 95.0, <82.8, 97.5, <114.1, 125.0, <92.1, 97.0, <109.14 , 127.3, <114.9, 130.0, <117.1, 132.5, <127.1, 130.0,
мин 32.1) 57.3) 74.6) 63.9) 70.0) 102.5) 98.5) 106.4) 117.0) 139.7) 108.9) 147.8) 148.8) 158.2) 147.7)
АТ 'к <0.0, 0.0, <0.0, 0.0, <0.0, 0.0, <0.0, 0.0, <0.0, 0.0, <0.0, 8.7, <0.0, 0.0, <0.0, 0.0, <0.0, 0.0, <14.1, 25.0, <0.0, 0.0, <0.0, 17.3, <0.0, 0.0, <0.0, 2.5, <0.0, 0.0,
мин 0.0) 0.0) 4.6) 0.0) 0.0) 22.5) 3.5) 6.4) 12.0) 39.7) 3.9) 38.7) 18.8) 28.2) 7.7)
Графики функций принадлежности ) временных характеристик технологического процесса для двух выбранных ПЗА (I = 1,3) отображены сплошными линиями на рис. 2. Штриховые отрезки соединяют графики функций принадлежности начала и окончания промежутков времени занятости ПЗА обслуживанием того или иного ВС и таким образом условно обозначают промежутки времени, в течение которых ПЗА определённо или возможно занят обслуживанием.
Рис. 2. Результаты оптимизации
Первый ( i = 1) ПЗА практически непрерывно занят обслуживанием ВС в следующей последовательности: j = 1,3,9,14. Непродолжительная пауза возможна (но маловероятна) лишь после заключительной операции обслуживания ВС j = 9 , так как окончание этой операции может произойти раньше наступления планового времени начала выполнения подготовительной операции для следующего ( j = 14 ) ВС. Занятость третьего ( i = 3) ПЗА несколько ниже. Он предполагается к обслуживанию трёх ВС в следующей последовательности: j = 1,7,10. Обслуживанием ВС j = 1 оба ПЗА, первый и третий, занимаются совместно, поэтому временные рамки процесса обслуживания этого ВС для обоих ПЗА совпадают. Между окончанием процесса обслуживания ВС j = 1 и началом процесса обслуживания ВС j = 7 ПЗА будет простаивать в течение промежутка нечёткой продолжительностью порядка <6, 12, 17) минут.
Затраты машинного времени на решение с помощью программной надстройки «Поиск решения» табличного процессора Microsoft Excel при использовании современной персональной вычислительной техники для рассмотренных исходных данных составляют не более 5 мин.
Заключение
Предложенная методика позволяет оптимизировать распределение технологических ресурсов разнообразных систем, в том числе авиатранспортных, в условиях неполной определённости исходных данных. Представленные результаты подтверждают достаточный уровень возможностей стандартного программного обеспечения для решения рассмотренной задачи на характерном примере системы наземного обеспечения полётов ВС в узловых аэропортах. Проведённый сравнительный анализ позволил выявить заметное различие результатов оптимизации при нечётких и чётких исходных данных, дефаззифицированных методом центра тяжести, что говорит о целесообразности нечёткой постановки рассмотренной задачи.
Библиографический список
1. Babic O. Optimization of refueling truck fleets at an airport // Transportation Research Part B. 1987. V. 21, Iss. 6. P. 479-487. DOI: 10.1016/0191-2615(87)90042-7
2. Carotenuto P., Giordani S., Salvatore A., Biasini A. Resource planning for aircraft refueling in airport parking area // Transportation Research Procedia. 2019. V. 37. P. 250-257. DOI: 10.1016/j.trpro.2018.12.190
3. Балашов В.Г., Заложнев А.Ю., Новиков Д.А. Механизмы управления организационными проектами. М.: ИПУ РАН, 2003. 84 с.
4. Голенко-Гинзбург Д.И. Стохастические сетевые модели планирования и управления разработками. Воронеж: Научная книга, 2010. 283 с.
5. Орловский С.А. Проблемы принятия решений при нечёткой исходной информации. М.: Наука, 1981. 206 с.
6. Лю Б. Теория и практика неопределённого программирования. М.: БИНОМ. Лаборатория знаний, 2005. 416 с.
7. Пегат А. Нечёткое моделирование и управление. М.: БИНОМ. Лаборатория знаний, 2013. 798 с.
8. Романенко В.А. Оптимизация параметров системы трансферных авиаперевозок с учётом нечёткой и стохастической неопределённостей // Управление большими системами. 2015. № 41. С. 285-313.
OPTIMIZING THE ALLOCATION OF TECHNOLOGICAL RESOURCES OF AN AIR TRANSPORT SYSTEM IN THE PRESENCE OF A SCHEDULE AND FUZZY INPUT DATA
© 2021
V. A. Romanenko Candidate of Science (Engineering), Associate Professor of the Department of Transportation Management and Control; Samara National Research University, Samara, Russian Federation; vla [email protected]
The problem of optimal allocation of technological resources (operators) of a technical or organizational-technical system, designed to serve certain objects (operands) according to a given schedule, is considered. We take into account the necessity of incorporating, together with the main, preparatory and final operations, the possibility to select one or several operators for operand service into the service process, as well as the dependence of the operations duration on the factors characterized by uncertainty. Due to the supposed absence of statistics, expert-assigned indefinite values in the form of triangular fuzzy numbers are used. The optimization problem is formulated as a mathematical programming problem with a fuzzy criterion and clear-cut constraints, consisting in finding such a distribution of a given number of operators to serve each operand from a given set which minimizes the target function that takes into account deviations from the schedule (delay) with the service termination. Typical examples of systems for which the problem is relevant are the production complexes of air transport enterprises operating in conditions of uncertainty when it is necessary to ensure the regularity and safety of air transportation. A model example of solving the problem of allocating mobile refueling facilities at a hub airport, taking into account the peculiarities of its schedule, is presented. It is shown that the capabilities of standard personal computer software are sufficient for the solution.
Optimization; technological resources; operand; triangular fuzzy number; fuzzy mathematical programming; hub airport
Citation: Romanenko V.A. Optimizing the allocation of technological resources of an air transport system in the presence of a schedule and fuzzy input data. Vestnik of Samara University. Aerospace and Mechanical Engineering. 2021. V. 20, no. 3. P. 160-170. DOI: 10.18287/2541-7533-2021-20-3-160-170
References
1. Babic O. Optimization of refueling truck fleets at an airport. Transportation Research Part B. 1987. V. 21, Iss. 6. P. 479-487. DOI: 10.1016/0191-2615(87)90042-7
2. Carotenuto P., Giordani S., Salvatore A., Biasini A. Resource planning for aircraft refueling in airport parking area. Transportation Research Procedia. 2019. V. 37. P. 250-257. DOI: 10.1016/j.trpro.2018.12.190
3. Balashov V.G., Zalozhnev A.Yu., Novikov D.A. Mekhanizmy upravleniya organizatsionnymiproektami [Mechanisms of organizational project management]. Moscow: IPU RAN Publ., 2003. 84 p.
4. Golenko-Ginzburg D.I. Stokhasticheskie setevye modeli planirovaniya i upravleniya razrabotkami [Stochastic network models for development planning and management]. Voronezh: Nauchnaya Kniga Publ., 2010. 283 p.
5. Orlovskiy S.A. Problemy prinyatiya resheniy pri nechetkoy iskhodnoy informatsii [Problems of decision-making with fuzzy input information]. Moscow: Nauka Publ., 1981. 206 p.
6. Liu B. Theory and practice of uncertain programming. Heidelberg: Physica-Verlag, 2002. 388 p. DOI: 10.1007/978-3-7908-1781-2
7. Piegat A. Fuzzy Modeling and Control. Heidelberg: Physica-Verlag, 2001. 728 p. DOI: 10.1007/978-3-7908-1824-6
8. Romanenko V.A. Optimization of transfer air transportation system parameters considering fuzzy and stochastic uncertainties. Automation and Remote Control. 2015. V. 76, Iss. 8. P. 1500-1514. DOI: 10.1134/S0005117915080135
