Математическая модель работы рубильной машины с учетом ее технических отказов Текст научной статьи по специальности «Прочие сельскохозяйственные науки»

40

Труды БГТУ, 2016, № 2, с. 40-44

УДК 630*363.7:519.85

Е. А. Леонов, В. В. Игнатенко, Д. В. Клоков

Белорусский государственный технологический университет

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ РАБОТЫ РУБИЛЬНОЙ МАШИНЫ С УЧЕТОМ ЕЕ ТЕХНИЧЕСКИХ ОТКАЗОВ

В статье представлена математическая модель работы рубильной машины с учетом ее технических отказов на различных стадиях работы. Разработанная модель базируется на применении теории массового обслуживания, критериев вероятностей состояний.

Работа рубильной машины характеризуется следующими состояниями: S0 - рубильная машина исправна, но не производит сбор, подачу, измельчение древесного сырья и выброс щепы; S11 - рубильная машина исправна, осуществляется сбор и подача древесного сырья к механизму измельчения; S12 - рубильная машина исправна, осуществляется измельчение древесного сырья; S13 - рубильная машина исправна, осуществляется выброс щепы; S21 — отказ механизма подачи; S22 - отказ механизма измельчения; S23 - отказ механизма выброса щепы. В такой модели имеют место два типа потоков: древесины и отказов оборудования. Приоритетом пользуется поток отказов, так как при их наступлении они «обрабатываются» (осуществляется ремонт машин) в первую очередь.

Целью теоретических исследований являлось установление рациональных режимов работы оборудования для измельчения древесного сырья на щепу, что позволит обеспечить в дальнейшем эффективность применения уже созданных рубильных машин. Разработанная методология позволяет определять: оптимальную продолжительность подачи рабочего органа к измельчаемому древесному сырью (штабелю дров или отходов), а также оптимальные сроки восстановления работоспособности механизмов подачи, измельчения древесного сырья и выброса щепы рубильной машины в условиях ее производственной эксплуатации.

Ключевые слова: модель, рубильная машина, вероятность, технический отказ, параметры.

E. A. Leonov, V. V. Ignatenko, D. V. Klokov

Belarusian State Technological University

THE MATHEMATICAL MODEL OF CHIPPER WORK GIVEN ITS TECHNICAL FAILURES

The article presents a mathematical model of the chipper based on its technical failures at different stages of work. The model is based on the application of queuing theory, probability criteria states.

Work chipper has the following states: S0 - the chipper is working, but does not produce the collection, supply, wood chopping and chips ejection; S11 - the chipper is working, produce the collecting and wood supply; S12 - the chipper is working, produce the wood chopping; S13 - the chipper is working, produce the chips ejection; S21 - the wood feeding mechanism is out of order; S22 - the wood chopping mechanism is out of order; S23 - the chips ejection mechanism is out of order. In this model, there are two types of flows: wood and equipment failures. Priority is given to the flow of failures, since they occur they are "processed" (carried out repair of machinery) in the first place.

The purpose of theoretical research was to establish rational modes equipment for grinding of raw wood on the wood chips, which will provide in the future effectiveness of established chippers. The developed methodology allows to determine: the optimal duration of the working body to chopped wood raw material (pile of logs or wood residues) as well as the optimal timing recovery performance of wood feeding, chopping and chips ejection mechanisms in terms of its manufacturing operation.

Key words: model, chipper, probability, technical denial, parameters.

Введение. С 2011 г. в Беларуси в различных министерствах и ведомствах, а также на частных предприятиях на древесном топливе работало более 3000 котлов мощностью от 0,012 до 20 МВт, а также 11 мини-ТЭЦ с установленной электрической мощностью от 1,2 до 4,23 МВт и тепловой мощностью от 6,5 до 16,4 МВт. Ввод в действие данных

объектов, работающих на биотопливе, потребовал решения задачи гарантированного обеспечения их сырьем. С этой целью в системе предприятий Минлесхоза созданы мощности по производству топливной щепы в 47 лесхозах с объемом производства около 500 тыс. т условного топлива, или около 2 млн. пл. м3 [1].

Е. А. Леонов, В. В. Игнатенко, Д. В. Клоков

41

Для обеспечения выполнения показателей, заложенных в целевых программах, у потребителей должны аккумулироваться потоки древесного сырья из различных источников. Принципиальное отличие работы мини-ТЭЦ, как потребителей древесного энергетического сырья, заключается в сезонной аритмии. Помимо этого, поступающее из различных источников сырье имеет широкий диапазон размерно-качественных характеристик.

Обеспечение устойчивого снабжения мини-ТЭЦ энергосырьем путем разрешения технологических и организационных противоречий требовало инновационного подхода к решению задачи создания его запасов.

В условиях сезонного характера производства и потребления энергетического сырья решение данного вопроса возможно при создании гибких лесоэнергетических терминалов (ЛЭТ)[1]. Именно ЛЭТ по организационной структуре и технологии наилучшим образом отвечает требованию переработки древесного сырья в широком диапазоне размерно-качественных характеристик, дает возможность хранить его и подготавливать к использованию в соответствии с запросами мини-ТЭЦ.

Отличительными особенностями ЛЭТ от традиционных лесных складов являются:

- применение мобильной системы специализированных машин;

- гибкий технологический процесс измельчения сырья, допускающий изменение мест и зон работы машин, хранения сырья и древесного топлива;

- переработка древесного сырья в широком диапазоне размерно-качественных характеристик;

- возможность:

выбора и изменения места расположения ЛЭТ в транспортно-технологической схеме

освоения ресурсов сырья в зависимости от конкретных производственных условий;

- разделения ЛЭТ на несколько составных частей;

- функциональной и территориальной интеграции с другими структурными образованиями (лесными складами, биржами сырья, деревообрабатывающими производствами, энергообъектами и др.) [1].

Примеры разработанных и внедренных на предприятиях системы Минлесхоза технологических схем функционирования ЛЭТ приведены на рис. 1.

Реализация разработанных технологических процессов производства топливной щепы в условиях гибких ЛЭТ требует рациональной организации работы применяемого для этих целей оборудования, в том числе передвижных рубильных машин - основного оборудования в технологической цепочке производства древесного топлива [2-5].

Основная часть. Математическая модель работы рубильной машины разработана с учетом возможных технических отказов (рис. 2). Для нее характерны следующие состояния: 50 -рубильная машина исправна, но не производит сбор, подачу, измельчение древесного сырья и выброс щепы; - рубильная машина исправна, осуществляется сбор и подача древесного сырья к механизму измельчения; ^12 - рубильная машина исправна, осуществляется измельчение древесного сырья; 513 - рубильная машина исправна, осуществляется выброс щепы; <£21 - отказ механизма подачи; 822 - отказ механизма измельчения; Б23 - отказ механизма выброса щепы.

В такой модели имеют место два типа потоков: древесного сырья и отказов отдельных узлов оборудования. Приоритетом пользуется поток отказов, т. к. при их наступлении они «обрабатываются» (ремонт) в первую очередь.

а б

Рис. 1. Технологические схемы функционирования ЛЭТ при производстве топливной щепы: а - на промежуточном складе; б - на автономном терминале

Из свободного состояния S0 в рабочее S11 система переходит с интенсивностью X подачи рабочего органа к древесному сырью (штабелю дров или валу древесных отходов). В дальнейшем осуществляется обработка древесины посредством ее перемежения через подающий механизм (состояние Sn), механизм измельчения (состояние Si2) и механизм выброса щепы (состояние Si3) с интенсивно-стями обработки древесного сырья ц и ц2 соответственно.

Обратный переход системы в свободное состояние осуществляется посредством выброса полученной щепы в кузов автощеповоза или бурт с интенсивностью ц3.

При наступлении отказа механизма подачи (состояние S2i), измельчения (состояние S22) или выброса щепы (состояние S23) система с соответствующей интенсивностью (Xu, Х12 или X13) перейдет из рабочего состояния (S11, S12 или S13) в соответствующее ему состояние отказа (S21, S22 или S23). После выполнения ремонтных работ вышедшего из строя узла рубильной машины с темпом Ц11, Ц12, Ц13 система вернется в одно из рабочих состояний (S11, S12 или S13). Неизвестные параметры X, Ху и цгу устанавливаются следующим образом: X = 1/tH, где tH - продолжительность времени подачи рабочего органа к древесному сырью (штабелю дровяной древесины или валу отходов), мин; ц = Шпод, Ц2 = 1/tизм, Цз = 1/tBH5, где ^од, соответственно продолжительности

¿изм

времени перемещения древесного сырья, его измельчения и выброса щепы, мин; Хц = Шподот, где ¿тдот - продолжительность времени между отказами механизма подачи древесного сырья, мин; ц:2 = 1/^одв, где ^^ -продолжительность времени восстановления работоспособности механизма подачи, мин; Х12 = 1/^измот, где ¿измот - продолжительность времени между отказами механизма измельчения древесного сырья, мин; ц12 = 1/^измв, где ¿измв - продолжительность времени восстановления работоспособности механизма измельчения, мин; Х13 = 1/^выбОТ, где ¿выб™ - продолжительность времени между отказами механизма выброса щепы, мин; Ц13 = 1/?выб , где ¿выб -продолжительность времени восстановления работоспособности механизма выброса щепы, мин.

Функционирование рассматриваемой рубильной машины описывается системой дифференциальных уравнений Колмогорова (1), полученных на основании размеченного графа состояний (рис. 2), где Р0(0, РЮ, г = 1,2, I = 1,2,3 - вероятности соответствующих состояний рубильной машины в момент времени ¿:

~Р~ = ЦзР13 -^0'

dB

12

dB.

= ХР0 - (Х„ + Ц 1)P11 + Цпр21, = Ц1 P11 - (Х12 + Ц2)P12 + Ц12Р22' ■= Ц 2 Р12 - (Х 13 + Ц3 )р13 + Ц13 Р23 '

(1)

dP21

21

dP

dP2

23

dt

= Х 11Р11 - Ц 11Р2 = Х 12 Р12 - Ц 12 Р22 ,

= Х 13Р13 - Ц 13Р23 *

2 3

Условие нормировки ^^Рг] + Р0 = 1.

г=1 I=1

Начальные условия для решения системы (рубильная машина исправна и свободна): Ро(0) = 1, ру(0) = 0.

So Х S11 Ц1 S12 Ц2 S13 Ц3

Х1 L111 Х12 i -112 Х13 J Ц13

S21 S22 S23

Рис. 2. Размеченный граф состояний рубильной машины

С учетом того, что моделирование и исследование работы рубильной машины производится на протяжении длительного промежутка времени (месяц, год и т. д.), значения вероятностей ро(0, Pi,j(t) будут колебаться относительно определенного среднего значения. В этом случае можно допустить, что P0 « const, Pu « const, P12 « const, P13 « const, P21 « const, P22 « const, P23 « const. Ошибка при принятии данного допущения не превышает 8% [2]. Система дифференциальных уравнений (1) примет вид:

0 = цзP3 -XPo,

0 = XPo - (X11 + ц 1)Pi + H-iiP>i, 0 = Ц1 P11 - (X12 + ц 2)P12 + ^2^2 ,

0 = ц2P12 - (X 13 + ц3)P13 + Ц13(2)

0 = X 11P11 - ц nP21,

0 = X 12P12 - ц 12P22 ,

0 = X 13 P13 - ц 13 P23 .

Е. А. Леонов, В. В. Игнатенко, Л. В. Клоков

43

Решив систему уравнений относительно вероятностей состояний Po, Pli, P12, P13, P21, P22, P23, получим выражения для расчета режимов работы рубильной машины:

P -АР.

1 il 1 0'

М1

А

М2

P =_P •

-42 _ 0 ?

А

М3

P =-P •

J13 0 '

я,

P -р AP•

1 21 Fil 1 0 ' М1

P -р AP•

J 22 r12 0 '

М2

P -р AP

1 23 K13 0 '

М3

Я12 . _ Я13

(3)

(4)

(5)

(6)

(7)

(8)

где рц =—; Р12 =—; Р1з =■

Ми М12 М13

С учетом условия нормировки и выражений (3)-(8), получим уравнение для расчета вероятности состояния Р0:

P0 -

1+ Х

(Р11 + 1) , (Р12 + 1) , (Р13 + 1)

М1

М2

М3

. (9)

Подставив в выражения (3)-(8) значения вероятности Р0, найдем значения вероятностей

P1, P11, P12, P13, Р21, P22, Р23.

Полученные зависимости вероятностей состояний рубильной машины позволяют установить рациональные значения ее параметров. Технология работы с зависимостями следующая. На основе технических характеристик принимается ряд параметров, например, ц1, ц2,

ц3, Х11, Х12, Х13 и из построенных зависимостей устанавливаются искомые параметры. Например, X, Цп, Ц12 либо Ц13.

На рис. 3 приведен пример установления одного из названных параметров.

Вероятность, Рр

1 1

1 1

1 1

1 1

1 I

1 1

1 |

1 1

1 т

Д21

Д21

Рис. 3. Зависимость вероятностей состояний рубильной машины

Установленное рациональное значение ц21* позволяет определить рациональную, в данном случае, продолжительность восстановления механизма измельчения:

t " -

изм

1

*

Ml 2

При этом обеспечивается надлежащая производительность машины, т. к. достигается практически максимальная ее величина Pv (вероятность работы).

Заключение. Разработанная модель позволяет при заданных характеристиках работы исследуемой рубильной машины получать рациональные режимы подачи сырья на обработку и ремонта технологического оборудования в случае технических отказов. Это приведет к росту производительности оборудования без существенных финансовых затрат.

Литература

1. Леонов Е. А. Устойчивое снабжение энергообъектов древесным топливом с созданием минимально необходимых запасов // Труды БГТУ. 2014. № 2: Лесная и деревообраб. пром.-сть. С. 17-19.

2. Клоков Д. В., Леонов Е. А., Турлай И. В. Модель работы форвардера с учетом надежности // Труды БГТУ. 2015. № 2: Лесная и деревообраб. пром-сть. С. 23-26.

3. Клоков Д. В., Турлай И. В., Леонов Е. А. Оборудование лесопромышленных предприятий. Лабораторный практикум. Минск: БГТУ, 2015. 200 с.

4. Технология и оборудование комплексного использования древесного сырья. Практикум / А. С. Федоренчик [и др.]. Минск: БГТУ, 2014. 274 с.

5. Федоренчик А. С., Клоков Д. В., Леонов Е. А. Энергетическое использование древесной биомассы. Практикум. Минск: БГТУ, 2015. 212 с.

References

1. Leonov E. A. Sustainable supply of wood fuel energy facilities with the establishment of minimum required reserves. Trudy BGTU [Proceedings of BSTU], 2014, no. 2: Forest and Woodworking Industry, pp. 17-19 (In Russian).

2. Klokov D. V., Leonov E. A., Turlay I. V. Model of forwarder work in respect of reliability. Trudy BGTU [Proceedings of BSTU], 2015, no. 2: Forest and Woodworking Industry, pp. 23-26 (In Russian).

3. Klokov D. V., Turlay I. V., Leonov E. A. Oborudovanie lesopromyshlennykh predpriyatiy. Labora-tornyypraktikum [Equipment timber companies. Laboratory practicum]. Minsk, BGTU Publ., 2015. 200 p.

4. Fedorenchik A. S., Mokhov S. P., Klokov D. V., Leonov E. A. Tekhnologiya i oborudovanie kompleksnogo ispol'zovaniya drevesnogo syr'ya. Praktikum [Technology and equipment for the integrated use of wood raw material. Practical work]. Minsk, BGTU Publ., 2014. 274 p.

5. Fedorenchik A. S., Klokov D. V., Leonov E. A. Energeticheskoe ispol'zovanie drevesnoy biomassy. Praktikum [Energy use of wood biomass. Practical work]. Minsk, BGTU Publ., 2015. 212 p.

Информация об авторах

Леонов Евгений Анатольевич - кандидат технических наук, старший преподаватель кафедры лесных машин и технологии лесозаготовок. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). E-mail: debager13@rambler.ru

Игнатенко Василий Васильевич - кандидат физико-математических наук, доцент, доцент кафедры высшей математики. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). E-mail: turlai@belstu.by

Клоков Дмитрий Викторович - кандидат технических наук, доцент кафедры лесных машин и технологии лесозаготовок. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). E-mail: klokoff_dm@belstu.by

Information about the authors

Leonov Evgeniy Anatol'evich - PhD (Engineering), Senior Lecturer of the Department of Logging Machinery and Technology. Belarusian State Technological University (13а, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: debager13@rambler.ru

Ignatenko Vasiliy Vasil'evich - PhD (Physical and Mathematical), Assistant Professor, Assistant Professor of the Department of Higher Mathematics. Belarusian State Technological University (13а, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: ignatenko@belstu.by

Klokov Dmitriy Viktorovich - PhD (Engineering), Assistant Professor of the Department of Logging Machinery and Technology. Belarusian State Technological University (13а, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: klokoff_dm@belstu.by

Поступила 15.02.2016