УДК 621.9
А.Н. Заикин, Л.И. Евельсон, Е.Г. Изюмова
Брянская государственная инженерно-технологическая академия
Евельсон Лев Игоревич родился в 1962 г., окончил в 1985 г. Брянский институт транспортного машиностроения, кандидат технических наук, профессор, заведующий кафедрой информационных технологий Брянской государственной инженерно-технологической академии. Имеет более 60 печатных работ в области математического моделирования, систем автоматизированного проектирования машин, информационных систем. E-mail: levelmoscow@mail.ru
Изюмова Елена Геннадьевна родилась в 1982 г., окончила в 2004 г. Брянский государственный университет, старший преподаватель кафедры информационных технологий Брянской государственной инженерно-технологической академии. Имеет 8 печатных работ в области дистанционных технологий, информатизации и автоматизации лесозаготовок. E-mail: e-izumova@mail.ru
ПОСТАНОВКА ЗАДАЧИ ОПТИМИЗАЦИИ ОБЪЕМОВ ЗАПАСОВ И РЕЖИМОВ РАБОТЫ ЛЕСОСЕЧНЫХ МАШИН
Приведены особенности создания межоперационных запасов в технологическом процессе лесосечных работ и необходимость повышения их эффективности. Сформулирована постановка задачи многокритериальной условной оптимизации объемов запасов и режимов работы лесосечных машин. Предложена реализующая ее имитационная компьютерная модель.
Ключевые слова: лесозаготовки, страховой запас, гарантийный запас, комплект машин, целевая функция, оптимизация, имитационная модель.
Лесозаготовки в нашей стране выполняются в различных природно-производственных условиях предприятиями разных форм собственности. Организация лесозаготовительного производства определяется его технологическими процессами, которые характеризуются составом применяемой техники (оборудования), количеством, последовательностью и местом выполнения технологических и транс-портно-переместительных операций в ходе осуществления технологического процесса.
Первой фазой лесозаготовок являются лесосечные работы. Основные из них включают ряд операций: валку, трелевку, обрезку сучьев, раскряжевку, погрузку. Лесозаготовки всегда начинаются с валки (В) деревьев, остальные операции могут выполняться в разной последовательности, быть совмещены, обрезка сучьев (С) и раскряжевка (Р) могут отсутствовать,
трелевка (Т) и погрузка (П) осуществляются в любом случае. В итоге выделяют разные технологические системы для лесозаготовок, например: В+Т+С+П, В+С+Т+П, В+Т+С, В+Т+П, ВСР+Т+П, ВТ+СР+П, В+Т, ВСР+Т, ВТ+СР, ВТ+П. Одной из особенностей лесосечных работ является создание между каждой парой смежных операций запаса ресурсов.
В нашей работе приняты следующие условные обозначения: Zс -объем страхового запаса, м3; 2тх - технологический запас между операция-
33
ми, м ; /тн - технический запас, м ; Zcp - средний ликвидный запас древесины на 1 га, м3; Рсм - сменная производительность сучкорезной или погрузочной машины, м /см.; Тдн - число дней, за которое будет выработан объем технологического запаса (как правило, один день); Р/, Р8 - сменная производительность основных машин на /-й и 5-й операциях (как показали
проведенные нами исследования, их значение необходимо принимать как среднестатистическую величину за последних три года для отдельных месяцев разработки лесосек для конкретного комплекта машин с учетом планируемых условий), м3; Pi max - наибольшая производительность из рассматриваемых операций, м ; Тсм - продолжительность смены, ч; Zr - объем гарантийного запаса, м3; Qi - объем выработки основных машин на i-й операции, м3, Qi = PjNjChsm; QS - то же на s-й операции, м3, QS = PsNschsm; Тп -продолжительность расчетного периода (месяца), дн.; QiD - уровень, до которого увеличивается объем выработки (Qi) на i-й операции после подключения дополнительной машины с объ-
3 D
емом выработки Qd/, м , Qi = Qi+ Qd/; Qmax - максимальный объем выработки
3D
машин на одной из операций, м ; QS -уровень, до которого увеличивается объем выработки (QS) на s-й операции после подключения дополнительной машины с объемом выработки QDS, м3, Qsd = Qs + Qds; Ni, Ns - число основных машин (на i-й и s-й операциях); N^ -число дополнительных машин; Р/д -сменная производительность дополнительной машины, м3; Qz - общий запас на лесосеке, м ; Lcp - длина зоны безопасности, м; chsm - число смен работы машин; chrd - число рабочих дней в месяце; Pa - сменная производительность машин на вывозке, м3; kmes -число месяцев разработки лесосеки.
Проведенные нами исследования показали, что запасы ресурсов многоуровневые. Первый уровень -страховой объем запасов (ZG), который состоит из объема технологических запасов (Z^), обеспечивающих зону безопасности, и объема технических запасов (ZpH), страхующих от простоев технические средства (ТС) на смежных операциях, когда одно из них выходит из строя, и рассчитывается по формуле
7 = 7 + 7 Технологический запас между операциями валки и трелевки составляет объем древесины, заготовленной с площади, равной по длине (Ь) и ширине (В) зоне безопасной работы, и может быть рассчитан по формуле Же р
Z -тх 10000
При L = 50.. .250 м; B = 50 м ZTX = (0,25.. ,1,25)ZCp.
Объем технологического запаса между операциями трелевки и обрезки сучьев, обрезки сучьев и погрузки должен обеспечивать маневренную работу технических средств и безопасность рабочих и будет прямо пропорционален производительности машин, например сучкорезных и погрузочных, и времени, за которое этот запас будет выработан:
7 = Р Т
тх см дн •
Объем технического запаса можно определить, рассмотрев простейшую линию из двух технических средств (ТС), производительности которых в единицу времени на предыдущей i-й и следующей 5-й операциях соответственно равны Pi и PS, а коэффициенты технической готовности Ктг1 и Ктг£. Проведя несложные преобразования, мы получили математическую зависимость для определения этого уровня запаса при условии бесперебойной работы ТС в течение смены:
Zm Pi max[2 (Ктг'^ Ктг£)]^см.
Второй уровень - гарантийный объем запасов (Zr), который слагается из страхового и оперативного, обеспечивающего компенсацию отклонений фактических объемов и времени поступления сырья на очередную операцию от расчетных значений. В целях максимального сокращения простоев из-за погодных, организационных,
Таблица 1
Формулы для расчета запасов
Расчетный период (месяц) Вариант Предельный запас с учетом максимальной выработки комплекта
Первый 2/ < 2* , [Тп (2?"2тахХ2тах" 2г) + 2г Ш г_ (2?-2тах)(2тах" 2.) + 2. 2И
2/ > 2* . _ (Тп2тах ^)(2тах " 2* )2 " 2тах) + ^2тах2Д* 2тах2оЗ
Второй -предпоследний 2/ < 2* V Тп (2тах 2: )2 "2тах) + ^ 2о, г 2*
2г > 2* гу Тп (2тах 2* )(2* - 2тах) + ^с г 2^
Последний 2/ < „ _ 2тах[Тп (2тах 21 )2 "2тах) + Za 2и, ] г (2тах 21Ж - 2тах) + 2тах2* /
2/ > 2* . (Тп 2* -¿сХ2тах " 2* )(2*~2тах) + ^ 2*2о* г 232о3
технических и других причин он рассчитывается по формулам табл. 1 [1].
Операции технологического процесса выполняются различными техническими средствами с разной производительностью. Поэтому, чтобы не пересечь границы страхового и гарантийного запасов, на операции лесозаготовительного процесса приостанавливают работу технического средства с большей производительностью в комплекте. Если на следующей технологической операции объем выработки машины окажется больше, чем на предыдущей, то возможен простой оборудования из-за отсутствия древесины. И, наоборот, если на предыдущей операции объем запаса увеличится, тогда на следующей машина не будет справляться с его потреблением. Объем выработки таких комплектов машин, как правило, равен минимальному объему выработки на одной из операций.
Как показали проведенные нами исследования, для повышения эффективности работы комплекта машин (возрастания объемов выработки и
снижения себестоимости лесозаготовок), целесообразно повысить объем выработки отстающих машин в комплекте до максимального уровня (объема выработки ведущей машины) за счет увеличения численности и (или) сменности работы машин на этих операциях (далее подключение дополнительных машин). Такая организация работы машин будет способствовать поддержанию страхового и гарантийного запасов на определенном, рассчитанном для конкретных условий, уровне [1]. В качестве дополнительного комплекта необходимо на валке и обрезке сучьев иметь бензиномоторную пилу, а на трелевке чокерный трактор. Этот комплект необходим также и для освоения объема древесины (до 5 %), который современные машины по своим техническим параметрам не могут освоить [3].
На рис. 1 представлены два варианта работы комплекта машин на несколько месяцев разработки лесосеки (в первый месяц (Т1) запасы создаются пополняются (72) и потребляются (73),
2, м3
'То
Технический уровень
Технологически{ 7*, уровень
Г
Ъ
Г
т, дн.
Рис. 1. Изменение запасов древесины в течение всего периода разработки лесосеки в зависимости от соотношения объемов их пополнения и потребления:
1 - при Qi<QmaX, 2 - при 2тах > Qs
во второй и следующие месяцы (Т) запасы пополняются (¿2) и потребляются (¿3), а в последний месяц (Т) пополняются (¿2), потребляются (¿3) и вырабатываются (¿4).
В первом варианте (при Qi<Qmax) за время ^ происходит создание гарантийного объема запаса машинами с объемом выработки Qi. В интервале времени ¿2 выполняются две операции технологического процесса (■ - предыдущая, 5 - следующая, 5 = 1). Так как Qi<Qs, то запас быстрее потребляется, чем пополняется. В итоге его уровень упадет до страховой величины. Дальнейшее снижение запаса допустить нельзя, так как нарушатся требования техники безопасности. Поэтому на «отстающей» операции Qi рекомендуется подключение дополнительной машины, тогда уровень запаса увеличивается до гарантийного в течение времени ¿3. Если лесосека разрабатывается в течение одного месяца, за период ¿4 происходит выработка запаса техническими средствами следующей операции с суточным объемом выработки Qs, когда машины предыдущей операции перебазировались на следующую лесосеку.
Во втором варианте (при Qmax>QS) за время ^ создается страхо-
вой запас техническими средствами предыдущей операции с суточным объемом выработки Qi. В интервале ¿2 выполняются две операции технологического процесса (^ 5). Так как Qmaх>QS, то запас быстрее пополняется, чем потребляется. В итоге его уровень поднимется до гарантийной величины. Дальнейшее увеличение запаса нежелательно, так как технические средства на «отстающей» операции Qs не будут справляться с его потреблением. На этой операции после подключения дополнительных машин уровень запаса снижается до страхового в течение времени ¿3. Если лесосека разрабатывается в течение одного месяца, за период ¿4 происходит выработка запаса техническими средствами с суточным объемом выработки Qs [2].
При планировании и организации процесса лесозаготовок необходимо управлять размерами межоперационных запасов, числом смен работы машин и временем их работы на операциях. Для определения этих параметров необходимо знать характеристики лесосеки и связи между ними.
Показатели лесозаготовительного процесса определяются, исходя из характеристик лесосек, отражающихся в технологической карте их разработки, особенностей организации процесса на предприятии, имеющегося парка машин, или рассчитываются.
Многие параметры лесозаготовительного процесса имеют вероятностный характер. Ярким подтверждением этого факта является производительность машин, которая зависит от климатических условий, квалификации работника, эксплутацион-ных характеристик и качества обслуживания машин [5-7]. Поэтому, как показали проведенные нами исследования, для расчета режимов работы
комплектов лесосечных машин необходимо производительность каждой из них принимать как среднестатистическую за последние три года для каждого месяца с учетом характеристики лесосеки и конкретного оператора.
Размещение технических средств по технологическим операциям определяется видом операций (В,Т,С,Р), числом машин (Ы, Ы/д), и их производительностью (Р:, Р/д). В итоге функция распределения технических средств по операциям лесозаготовок примет вид
Размер страхового запаса (2с) зависит от требований техники безопасности, характеристики лесосеки (2г, 2ср, Ьср) и может быть представлен в виде функции
2с = > 2сР 5 А;р
Гарантийный запас (2г) зависит от производительности (Р), числа машин (Ы), числа смен их работы (екш), числа рабочих дней (екгд) в расчетном периоде (месяце). Эта зависимость имеет вид
2Т = д,, ^, , сЬгй, сНят^
Для обеспечения максимальной загрузки отдельных машин и максимальной выработки всего комплекта при расчетах режимов работы необходимо оптимизировать объемы запасов, число смен работы основных и продолжительность работы дополнительных машин. Тогда задача оптимизации может быть сформулирована следующим образом: необходимо определить такое число смен работы основных машин и продолжительность работы дополнительных машин, при которых объемы гарантийных запасов, продолжительность разработки лесосеки и приведенные затраты будут минимальны. Оптимизационные методы позволяют выбрать наилучший ва-
риант без непосредственной проверки всех существующих [4].
Целевая функция задачи, определяющая параметры технологического процесса лесозаготовок, примет вид
2 = Ч(^2с, 2г) ^ тах;
Т = г (2с, 2г ) тт; 2 = z( р1,2 с, 2 г )->тт.
Задача относится к многокритериальной многомерной условной оптимизации. В качестве целевых функций могут выступать максимизация объемов выработки, минимизация времени работы машин, затрат на выполнение операций и т. п. Исследуемая функция зависит от большого числа параметров, ее аналитическое представление затруднительно. Поэтому целевую функцию следует задать процедурой, в которой рассчитываются основные параметры процесса лесозаготовок. На первом этапе в качестве критерия оптимизации выбрано время работы машин, а по другим целевым функциям задавались ограничения. Выполнена математическая постановка задачи оптимизации лесозаготовок.
Задача оптимизации заключается в минимизации Т — I , 2С,2, при ограничениях: = , 2ср, ¿ср „,
с2 = ё2(ра,рщ,рпмпмщ), д^адЛ),
2 = :(!•] ,1С,1Т) и управляющих параметрах: и = п(екчш, Ы Ы/д, кшеъ).
Оптимизация поставленой задачи вручную требует много времени и может привести к большим погрешностям. В целях сокращения времени оптимизации и повышения точности расчетов нами разработана имитационная модель функционирования лесосечных машин в течение всего периода освоения лесосеки с заданным объемом производства и с учетом подключения дополнительных машин, реализованная в программном обеспечении
Таблица 2
Формулы для расчета времени работы машин
Расчетный период (месяц)
Вариант
Время создания запаса
Время пополнения и потребления запаса
Время работы дополнительных машин
Время выработки запаса
Первый
Второй -предпоследний
Последний
а < а
а > &
а, < & & > &
а < а & > &
и =
и =
а а
и =
<2 =
и =
4 =
к =
и =
~2С
бтах -а
*0тах -О*
£?тах -а
бтах -а
V
О'Ш.Ж -е,
£?тах ~Оз
ч =
^ =
со-о ~Т> -о
-о8
^ =
Ч2та
Т с О
п яс-тах х^-г
4 =
Ог>,
_ ~ ах ~ & .
^ _ ^п^тах ¿г
О, ах Ой,
ч =
гп(<2юах-&)
а.
и =
и =
из
на ПЭВМ. Основные расчетные формулы, используемые в программе, приведены в табл. 1 и 2.
Обеспечение максимальной эффективности комплекта лесосечных машин достигается за счет объективных исходных данных. Входными данными для программы являются: число месяцев, номер первого месяца (в котором начинается разработка лесосеки), общий запас на лесосеке, средний запас древесины на 1 га, длина разрабатываемой ленты (зоны
ПРОГРАММА РАСЧЕТА РЕЖИМОВ РАБОТЫ ЛЕСОСЕЧНЫХ МАШИН
безопасности), число смен в один день и сменная производительность, число машин на вывозке, вид технологической системы, число рабочих дней в месяце, марка машин на каждой операции, их сменная производительность, число смен работы по основным машинам, марка машин, их число и производительность по дополнительным машинам. Ввод данных в программу осуществляется с помощью формы, представленной на рис. 2.
ОБЩИЕДАННЫЕ
повтор ввода!
Количество месяцев Номер первого месяца Марка
Общий запас на лесосеке (кбм) Средний запас (кбм/га) Длина зоны вырубки (м) Число смен в один день на вывозке Сменная производительность машин на вывозке (кбм/смену)
Подтвердите ввод общих данных
Технологическая система
г в»т»оп г в+Т+с+П С вср.т.п г вср-т г в»с»т»П с в.т.с г вт»ср»п г вт»ср
а в-с-т»п г в»т»п <- в+т г вт«п
Автоматический ввод исходных данных по основным и дополнительным машинам
повтор ввода I ВВОД ДАННЫХ ПО ОСНОВНЫМ МАШИНАМ
инка
нам
МЕСЯЦ ЧИСЛО РАБ ДНЕЙ ПАРАМЕТРЫ В |С |т п * в
МАРКА МАШИНЫ ЛП-19 ЛТ-157 ЛП-306 Тайга
1 Март 20 ЧИСЛО МАШИН 1 1 1 1
ЧИСЛО СМЕН 1 К 1 1
НОРМА ВЫРАБ01 210 170 165 180 -
ВВОД ДАННЫХ ПО ДОПОЛНИТЕЛЬНЫМ МАШИНАМ
месяц число раб. дней [параметры |в |С |т |П
марка машины Урал ТТ-4 Тайга Урал II
1 Март 20 число машин 1 1 1 1
норма выработ10 50 50 ро V
Подтвердите ввод данных по машинам нажатием на кнопку
ВЫПОЛНИТЬ РАСЧЕТЫ Рис. 2. Ввод данных в программу
Выходные параметры: объемы объемы производства основных и до-
страхового и гарантийного запасов, полнительных машин, суммарный
продолжительность работы машин по объем производства машин, ежеднев-
операциям на создании, пополнении, ная потребность машин на вывозке,
потреблении и выработке запасов, число сэкономленных дней (рис. 3).
Значения параметров по операциям
[месяц число раб. дней |1 вид запаса в-с |с-т | т-п
1 Март 18 1 страховой 200.00 200.00 180.00 ■
1 гарантийный 352.38 281.43 309.29
2 Апрель 20 2 вид запаса Л
ПРОСМОТРЕТЬ ИСХОДНЫЕ ДАННЫЕ
Продолжительность работы основных машин
Значения параметров по операциям
ВЫВОД РЕЗУЛЬТАТОВ НА ПЕЧАТЬ
[месяц |числ0РАБ дней [Период времени в-с с-т Т-П -
Шерт | 18 Создания запаса 0.95 0.95 0.86 а
Пополнения и потребления 3.81 1.81 4.31
Работа сдоп, машинами 15.24 16.29 12.93
2 Апоель 20 Пополнения и потребления 4.00 2.00 5.00
ВЫХОД ИЗ ПРОГРАММЫ
ОБЪЕМ ПРОИЗВОДСТВА И ЕЖЕДНЕВНАЯ ПОТРЕБНОСТЬ МАШИН ДЛЯ ВЫВОЗКИ
Объем производства машин
Месяц [Основных |Дополнительь|Всего [Потреб, маш. на бые[
Экономия дней 16
3800 00 2.92
Рис. 3. Результаты работы программы
С помощью разработанной программы можно реализовать имитационную модель процесса лесозаготовок. Варьируя по выбранному алгоритму управляющими параметрами и определяя объем выработки машин с учетом случайных факторов, можно находить значения целевой функции и проверять выполнение ограничений. Далее можно установить точку оптимума, используя известные методы оптимизации.
Решение поставленной задачи оптимизации позволит организовать лесозаготовительный процесс с учетом влияния случайных факторов, обеспечить минимальное время выполнения всего технологического процесса при максимальной выработке и минимизации денежных затрат.
СПИСОК ЛИТЕРАТУРЫ
1. Заикин А.Н. Теория, методы и модели интенсификации лесосечных работ: мо-ногр. Брянск: БГИТА, 2009. 204 с.
2. Изюмова Е.Г. Компьютерное моделирование режимов работы комплектов машин // Материалы Регион. науч. конф. студентов и аспирантов «Достижения молодых ученых Брянской области» / под. ред. И.А. Лаге-рева. Брянск: БГТУ, 2010. 301 с.
3. Макуев В.А. Научно-методологические основы формирования парка лесосечных машин предприятия: автореф. дис. ... д-ра техн. наук. М.: МГУЛ, 2010. 33 с.
4. Реклейтис Г., Рейвиндран А., Рэг-сдел ^.Оптимизация в технике. Кн. 1. М.: Мир, 1986. 350 с.
5. Тацюн М.В. Лесные технологии на фоне структурной перестройки ЛПК // Лесн. пром-сть. 1996. № 2. 22 с.
6. http://www.neosystems.ru
7. http://www.neweconomic.com
A.N. Zaikin, L.I. Evelson, E.G. Izyumova Bryansk State Academy of Engineering and Technology
Setting Optimization Task for Resources Volume and Logging Machines Operation Modes
The peculiarities of establishing in-process resources in the technological process of logging operations and necessity of increasing their efficiency are provided. The task setting of multi-criteria conditional optimization of resources volume and operation modes of the logging machines are formulated. The simulation model for its implementation is offered.
Keywords: logging operations, insurance resources, guarantee resources, set of machines, target function, optimization, simulation model.