CC BY
43
УДК 658.581.001.57 В.Н. Шиловский
Шиловский Вениамин Николаевич родился в 1945 г., окончил в 1970 г. Петрозаводский государственный университет, доктор технических наук, профессор кафедры технологии металлов и ремонта ПетрГУ. Имеет более 150 печатных трудов в области надежности и ремонта лесных машин.
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ОПТИМИЗАЦИИ ОБЪЕМОВ ПОСТАВОК И РАСПРЕДЕЛЕНИЯ ЗАПАСНЫХ ЧАСТЕЙ МЕЖДУ ТЕРРИТОРИАЛЬНО РАЗОБЩЕННЫМИ ПОТРЕБИТЕЛЯМИ
Представлена вероятностная модель обеспеченности и необеспеченности запасными частями лесозаготовительных машин, позволяющая определить оптимальный объем поставки и распределения запасных частей между территориально разобщенными потребителями.
Ключевые слова: лесозаготовительные машины, запасные части, оптимизация поставки, распределение, вероятностная модель.
Парк машин лесозаготовительного предприятия, как правило, состоит всего из нескольких десятков единиц и характеризуется малым числом отказов за период Тп, обычно принимаемый равным одному году. Поэтому задача определения оптимального числа запасных частей (ЗПЧ) носит чисто дискретный характер, что дает возможность решать ее с использованием схемы Бернулли (СБ) [1]. Использование СБ предполагает, что за рассматриваемый промежуток времени с машиной может произойти только одно из двух событий: объект откажет или не откажет.
Если вероятности двух, трех и более отказов одной машины за плановый период таковы, что этими событиями нельзя пренебречь, то необходимо многократное использование СБ, т. е. период разбивается на такие интервалы времени, в пределах которых вероятность более одного отказа достаточно мала для того, чтобы ее можно было не учитывать.
Разбиение будем считать удовлетворительным, если при дальнейшем удвоении числа интервалов новое значение математического ожидания (МО) числа отказов всех машин предприятия отличается от предыдущего не более чем на 1 ... 3 %. Это условие предполагает суммирование нескольких независимых простейших потоков отказов с любыми законами распределения [2].
Обозначим через д7 вероятность отказа одной детали в течение 7-го интервала времени с учетом возможных замен на предыдущих интервалах.
Вероятность отказа одной детали в течение 7-го интервала при отсутствии замен отказавших деталей на предыдущих интервалах
/, = ^ - (1)
где Е7 - интегральная функция распределения вероятности отказа детали в конце 7-го интервала.
Для первого интервала д1 = /1. Число деталей, отказавших в течение второго интервала, с учетом замен на первом, равно:
М2=(/12 +/2)М = (д1/1 +/2)М, (2)
где /1 - вероятность отказа деталей в первом интервале;
/2 - вероятность отказа во втором интервале при отсутствии замен деталей в первом;
(ЯХМ)- число отказавших деталей во втором интервале из ^М деталей, установленных взамен отказавших в первом;
/2М - число отказавших деталей во втором интервале из М деталей, установленных к началу первого интервала.
Для решения данной задачи примем следующий порядок распределения запчастей:
каждый лесопункт (ЛЗП) при отсутствии ЗПЧ вначале обращается к ближайшему их резерву - центральному складу (ЦС) или другому ЛЗП, затем к следующему по удаленности резерву и т. д.;
к резерву любого ЛЗП и ЦС вначале обращается ближайший ЛЗП, затем следующий по удаленности и т. д.
В формализованном виде эти правила обращения представляют собой множество упорядоченных пар натуральных чисел Ж7(с, О, где с - номер ЛЗП или ЦС, к которому обращается ЛЗП с номером t, причем с принимает одно из следующих значений: 1, 2, ..., Т, 5, а t = 1, 2, ..., Т. Индекс 7- номер обращения, который изменяется от единицы до Т2, где Т - число ЛЗП. Величина = Т + 1 обозначает номер ЦС.
Функционирование рассматриваемой системы парк машин - резерв ЗПЧ - внешний источник ЗПЧ исследуем с помощью следующей вероятностной модели. Для каждого ЛЗП и любого 7-го интервала времени вводятся две независимые случайные величины. Первая показывает, что на данном ^м ЛЗП в течение 7-го интервала произошло т отказов; вероятности этих событий Рц(т) задают выражением
= I (1 - ' (3)
(А', -т)\т\
где N - число машин на ^м лесопункте.
Вторая случайная величина описывает состояние резерва запасных частей на ЛЗП и показывает, что резерв данного ^го ЛЗП на интервале 7 содержит К ЗПЧ. Вероятность этих событий будет обозначаться через Q7¡t(K), К = 0, 1, ..., К Через произведение этих случайных величин задается еще одна группа событий, состоящая в том, что на данном ЛЗП не может быть устранено более п отказов. Вероятности этих событий будут обозначаться через У7л(п) и называться вероятностями п заявок, п = 0, 1, ...,
Знать вероятности этих трех групп событий достаточно для вычислений МО числа отказов на данном ЛЗП, устраняемых за счет собственного резерва ЗПЧ ЛЗП, резервов других ЛЗП, ЦС и внешнего источника; по этому числу отказов определяют МО суммарных затрат.
Нельзя установить, какой из отказов, происшедших в течение данного интервала времени, произойдет раньше, какой - позже. Считаем, что они происходят одновременно, поэтому резерв данного ЛЗП вначале используется машинами этого ЛЗП, и лишь затем происходит обращение к нему других ЛЗП. Обозначим обращение машин 7-го ЛЗП к собственному резерву запчастей К, через (обращение 7-го ЛЗП к к-му ЛЗП через / —» к.
Вероятность наличия т ЗПЧ на 1-м ЛЗП в интервале времени ] (после обращения 7 —» 7 на этом интервале к резерву К,) при предположении о независимости отказов в любых двух различных интервалах времени вычисляется по формуле
м
О] / (т) - X 0^и(л)Рл(с), (4)
' п-т '
с=п-т
где от = 1,2,
М = К1 при т+Ы(> К{, М = + т при т + < К(.
Значит,
0л,(О) = 1-£0л,(1н). (5)
т=1
Вероятность т заявок от 7-го ЛЗП в интервале ] после обращения 7 —> 7 равна:
м
У]-Лт)= I Р],{п)0]_и{с\ (6)
п=т 3 '
с=п-т
где т = 1, 2,
М = при т + К(> Щ М = К(+т при т + К(< Щ
Q]■-\,t (с) - вероятность наличия с ЗПЧ на 7-м ЛЗП по истечении (/ - 1)-го интервала.
Значит,
м,
= т). (7)
т=1
В начальный момент времени Qo,t(с) задают следующим образом: Qo,t(Kt) = 1; Qo,t(с) = 0 для с < К
Математическое ожидание числа отказов, устраненных на 7-м ЛЗП в течение ]-го интервала за счет резерва К, равно:
к,
М01= , (8)
к=1
где
К N,
А{К) = £ mPJt (т) + К X PJt (и) при Nt > К; (9)
т=1 п=К+1
К
А(К) = XmPJ t (т) при К > N,. (10)
т=1
Вероятность наличия т ЗПЧ на k-м ЛЗП на j-м интервале после обращения t —> к вычисляется по формуле, аналогичной формуле (4), но вместо Qj.\j (п) и Рр(с) будут соответственно О '¡Лп) ~ вероятность наличия на к-м ЛЗП на j-м интервале п ЗПЧ после обращения к —> к или т —» к. предшествовавшего обращению / —> и (с) - вероятность с заявок обращения t^-k:
м
QjAm)= Z Q'Jk(n)VJt{c), (11)
п-т J' J'
с=п-т
где m = 1, 2, ..., Kk;
M = Kk при m + Nt < Kk; M=Nt + m при m + Nt < Kk.
QM(0) = l~TQM(m). (12)
т=1
Вероятность m заявок от /-го ЛЗП на j-м интервале после обращения t —> £ вычисляют по формуле, аналогичной формуле (6):
м
VjAm)= Z V'jt{n)0jk{c), (13)
п=т J' J'
с—п—т
где т = 1, 2,
M=Nt при п + Кк> Nt; М=Кк + т при т + Кк< Nt.
Здесь V р(п) - вероятность п заявок от /-го ЛЗП на j-м интервале перед обращением t —> к, a Qjik(c) - вероятность наличия с ЗПЧ на к-м ЛЗП на j-м интервале в момент обращения / —> к.
N> ,
VJÁ0) = \-TV;^>)- (14)
т= 1
Математическое ожидание числа отказов, устраняемых на t-м ЛЗП в течение j-то интервала за счет резерва m-то ЛЗП при обращении t —> m, равно:
M02Jt(m)=^Qjm(K)AK), (15)
к=1
где
К N,
А{К) = ZcFjj (с) + К Z VJ t (и) при Nt > (16)
т=1 П=к+1
ы,
А{К)=Ъ^,(с)щъК>Ыь (17)
с=1
Здесь Qj,m(K) и У/,(о) - соответственно вероятность наличия к ЗПЧ на т-м ЛЗП и вероятность с заявок от /-го ЛЗП на /-м интервале при обращении / —» т.
Вычисления по формулам (11), (13), (15) для обращений / —> к производятся в последовательности, задаваемой стратегией обращений Б,(к, /). Центральный склад формально рассматривается как ЛЗП с резервом в К(Т + 1) ЗПЧ, на котором отсутствуют машины, и вероятность с заявок от ЦС (с > 0) равна нулю, а вероятность «нулевой» заявки - единице.
Математическое ожидание числа отказов, устраняемых на /-м ЛЗП в течение /-го интервала за счет резерва других ЛЗП и ЦС, равна сумме:
Т+1
Г2,., = ЕМ12,.,(|и). (18)
т=1
Математическое ожидание числа отказов, устраняемых на /-м ЛЗП в течение у-го интервала за счет ЗПЧ из внешнего источника, равно разности:
МО3^ = МО^ - МО1^ - Т2^, (19)
где МО ^ - МО числа отказов на всех машинах /-го ЛЗП в течение у-го интервала.
м,
МО = (20)
с=1
Частная задача, в которой каждый ЛЗП имеет доступ лишь к части общего резерва ЗПЧ предприятия (собственному резерву ЦС), решается аналогично, но для другой стратегии обращения. Для этого случая стратегия формулируется в виде следующих правил:
а) каждый ЛЗП вначале обращается к собственному резерву, а затем к резерву ЦС;
б) к резерву ЦС вначале обращается ближайший ЛЗП, затем следующий по удаленности и т. д.
Математическое ожидание общих затрат по всему предприятию при всех вариантах поставки и распределения запчастей за весь расчетный период равно:
2 = 21 + 22 + 23 + 24 + 25, (21)
где 21, 22, 23, 24, 25 - МО затрат, соответственно связанных с устранением отказов за счет собственных резервов ЛЗП; других ЛЗП или ЦС; ЗПЧ из внешнего источника; на создание резерва ЗПЧ; из-за недоиспользования ЗПЧ.
Минимум МО общих затрат (2) является условием оптимального объема поставки и распределения запасных частей между территориально распределенными потребителями лесозаготовительных машин.
Выводы
1. При определении потребности в ЗПЧ за критерий оптимальности необходимо принять минимум суммарных затрат, основными составляющими которых являются потери, связанные с простоями техники как из-за недостатка, так и переизбытка ЗПЧ.
2. Выбор оптимального резерва ЗПЧ невозможен без учета его распределения между отдельными группами машин.
3. Стратегия распределения запчастей между территориально разобщенными потребителями с учетом их оперативного взаимного обращения позволяет на порядок снизить простои ЛЗМ и одновременно общий потребный запас ЗПЧ.
СПИСОК ЛИТЕРАТУРЫ
1. Попов Ф.П. Обеспечение лесозаготовительной техники запасными частями. - М.: Экология, 1991. - 144 с.
2. Шиловский В.Н. Теоретические основы и стратегии организации маркетинга и менеджмента технического сервиса территориально распределенных машин и оборудования. - Петрозаводск: ПетрГУ, 2001. - 324 с.
Петрозаводский государственный университет
Поступила 10.02.02
V.N. Shilovsky
Mathematical Model of Supplies Volume Optimization and Spare Parts Distribution among Territorially Scattered Customers
The probability model of provision and non-provision of spare parts for forest-harvesting machines is provided allowing to set the optimal volume of spare parts supply and distribution among territorially scattered customers.
