УДК 621.311
СТОХАСТИЧЕСКИЕ ЗАКОНОМЕРНОСТИ ОПЕРАЦИОННЫХ ЦИКЛОВ ЛЕСОПИЛЬНЫХ РАМ
© Агеев С.П., д-р техн. наук, доц.
Северный (Арктический) федеральный университет имени М.В. Ломоносова, наб. Северной Двины, 17, г. Архангельск, Россия, 163002 E-mail: doctor.mart11@mail.ru
Характерной особенностью процессов деревообработки является то, что в них постоянно проявляется воздействие различных случайных факторов. Поэтому для более полного изучения основных закономерностей функционирования производственных процессов на деревообрабатывающих предприятиях необходимо рассматривать их как разновидности случайных процессов. Одним из наиболее энергоемких деревообрабатывающих производств является лесопильное производство. При этом лесопильная рама, со сравнительно небольшой околостаночной механизацией составляющая отдельный участок, как наиболее мощный электроприемник в значительной степени определяет характер его электропотребления. Цель настоящей статьи -установление вероятностных закономерностей рабочих циклов распиловки древесины на лесопильных рамах.
Для проведения исследований предложена математическая модель нагружения механизма резания главного привода лесопильной рамы, согласно которой лесопильная рама рассматривается как система массового обслуживания. Для анализа продолжительности распиловки отдельных бревен использован метод этапов, разработанный А.К. Эрлангом.
С использованием предложенной модели найден закон распределения длительности рабочих циклов лесопильной рамы, показана зависимость параметров закона распределения от ее режимных параметров и геометрических характеристик распиливаемого сырья.
В результате исследований установлено, что длительность рабочих циклов лесопильных рам как случайная величина подчинена обобщенному закону распределения Эрланга (к + 1)-го порядка. Согласно центральной предельной теореме при увеличении к обобщенный закон Эрланга неограниченно приближается к нормальному. Как показала практика моделирования технологических и энергетических процессов деревообработки, вполне удовлетворительную точность можно получить при аппроксимации законов Эрланга нормальным законом уже при к > 15. При этом средняя ошибка аппроксимации не более 0,001 %.
Полученные результаты могут быть использованы при анализе производительности и оптимизации параметров оборудования поточных линий, а также при решении вопросов нормировании электропотребления процессов лесопиления.
Ключевые слова: лесопильная рама, рабочий цикл лесопильной рамы, система массового обслуживания, плотность распределения вероятностей, дифференциальные уравнения, преобразование Лапласа, теорема свертывания, геометрические характеристики распиливаемого сырья.
Как известно, энергетические и технологические режимы в большинстве производственных процессов деревообработки взаимосвязаны. Поэтому при анализе электропотребления необходимо одновременно учитывать и особенности протекания технологических процессов с учетом влияния факторов, характеризующих как сам процесс, так и отдельные его операции.
Одним из наиболее энергоемких деревообрабатывающих производств является лесопильное. При этом лесопильная рама (ЛР) со сравнительно небольшой околостаночной механизацией составляет отдельный участок и как наиболее мощный электроприемник в значительной степени определяет характер его электропотребления [3].
Характерной особенностью процесса лесопиления является то, что в нем постоянно проявляется воздействие различных случайных факторов (форма и размеры пиловочных сортиментов, изменчивость продолжительности технологических операций и др.) [1]. В результате подача сырья к технологическим агрегатам происходит не ритмично, а с разными интервалами времени. Такие процессы не могут быть детерминированными, их необходимо рассматривать как разновидности случайных процессов.
Целью данной статьи является установление вероятностных закономерностей рабочих циклов распиловки древесины на лесопильных рамах, необходимых для решения задач нормирования электропотребления процессов деревообработки.
В работе [2] была представлена математическая модель нагружения механизма резания главного привода ЛР. Согласно этой модели лесопильная рама рассматривается как система массового обслуживания (СМО). Для анализа продолжительности распиловки отдельных бревен
был использован метод этапов, разработанный А.К. Эрлангом [6]. По этому методу распиловка бревен рассматривается как случайный дискретный процесс, состоящий из ряда последовательно выполняемых этапов (стадий).
Пусть для распиловки некоторого бревна требуется выполнить к стадий процесса. Тогда моменты начала стадий будут характеризовать следующие состояния СМО: Е1 ,Е2, ■ ■ ■ ,Ек. Рассматриваемая СМО действует таким образом, что когда обслуженное требование (бревно) покинет СМО, на ее вход поступает новое требование. Это новое требование занимает этап 1 (СМО переходит в состояние и находится в нем случайное время После ухода из этапа 1 требование поступает на этап 2 (СМО переходит в состояние Е2) и будет находиться там случайное время Т 2 и т.д. После окончания к-й стадии требование покидает СМО, а сама СМО переходит в состояние Е0 , обусловленное межторцовыми разрывами между бревнами. Только после этого в СМО может поступить новое требование, и она вновь переходит в состояние
Найдем закон распределения длительности рабочих циклов Тц лесопильной рамы. Эта задача сводится к нахождению закона распределения времени пребывания СМО в состояниях Полагая, что переход требования от одного этапа к следующему осуществляется под воздействием пуассоновского потока событий с интенсивностью X, построим размеченный граф состояний СМО (рис. 1). Состояние Е11 является поглощающим. Переход в него означает окончание текущего рабочего цикла и начало следующего.
Рис.1. Размеченный граф состояний Соответствующая этому графу система дифференциальных уравнений имеет следующий
вид:
= (0; = Хр2 (0 + Ар1 (0;
йъ = + Хрй_!(0; йъ = -АоРоСО + СО-
Начальные условия, при которых нужно интегрировать эти уравнения: Р1(0) = 1; р2( 0) = Рз(0) = - = ро(0) = 0. В соответствии с [1] вероятность
Р1ОО =
Для отыскания вероятности р 2 ( 0 получим дифференциальное уравнение
+Ар2(0 = Хе АГ. Преобразуем это уравнение по Лапласу:
зР2(з) + ХР2(з) =
5 + X
где 5 — параметр преобразования; изображение функции Следовательно,
Р20)= "
(5 + А)2 '
Используя таблицы изображения основных функций [7], найдем оригинал:
Р2(0 =
Рассуждаем аналогично:
)\ ] 1 = 1,2.....к■
Для отыскания вероятности составим дифференциальное уравнение
йр^Ь)
+ ¿оРо(0 =
Преобразуя это уравнение по Лапласу, получим
Xfc-1
sP0(s) + X0P0(s) = APfe(s) = A-
PoW = ТГТТТТГТТ^ ■ W
(5 + А)/с '
откуда
А*
(5 + Ао)(5 + Х)*
Для нахождения оригинала функции р 0 ( С) рассмотрим один из возможных методов, использующих теорему свертывания [7]. Для этого представим изображение (1) как произведение двух изображений:
Ро (5 ) = Ро 1 (5 ) Ро 2 (5 ) , (2)
где
Согласно теореме свертывания произведение (2) есть ¿-изображение функции р о ( С ) , определенной интегралом с переменным верхним пределом:
Ро О) = Р0 1 О) Р0 2 О ) - | Р о 1 00 Р о 2 ( t - х) dх . ( 3 )
Приняв
tk-1
согласно (3) получим
или
Г xfe_1
Р0 (s) = Afc I --—е-е- я° (i - x)dx ,
J (к — 1)! о
р0 ( t) = Afc g ° I xfc - ^ (Я° - Я) dx.
(к — 1)! J
о
Функция распределения F( t) = р( Гц < ^равна вероятности того, что к моменту времени i СМО выйдет из состояния £0 (закончится холостой ход) и окажется в состоянии £1 (лесопильная рама приступит к распиловке очередного бревна).
Следовательно, F ( t) = р 11 ( t) , откуда плотность распределения случайной величины Тц
frt, F(t) dPll(t)
или
/ ( t) = A0Afc g °° I xfc - 1e*(A°-X)dx. (4)
(к — 1)! J о
Закон распределения, имеющий такую плотность вероятности, называется обобщенным законом Эрланга (/с + 1 ) -го порядка, так как случайная величина 7Ц является суммой (/с + 1 ) независимых случайных величин, распределенных по показательному закону. У этого закона k параметров равны X, а один - X0.
Для практического использования формулы (4) вычислим интеграл методом по частям.
Тогда
к
f(t) = X0Xke~Xot
( 1)]+1 tk J ^e(Ao_A)t+.
¿(/-/ ) !(An-
¿^ (k-jy.Oo-XV (X0 ~ W
7=1
(5)
Параметры закона распределения могут быть определены следующим образом. В работе[2] было получено
где М Гр — математическое ожидание времени распиловки одного бревна, с; ЙГр — дисперсия времени распиловки одного бревна, с2. Учитывая, что
(Ь\ М1 мТр=м(й)=—
имеем
мь
х = иЪь'
где и - скорость подачи, м/с;
Ь - длина бревен, м; МЬ, ОЬ - математическое ожидание и дисперсия Ь в сортировочной партии.
Таким образом, параметр X зависит от геометрических размеров распиливаемого сырья и скорости подачи бревен.
Как показали экспериментальные исследования, проведенные на лесозаводах г. Архангельска, длительность вспомогательного времени за операционный цикл (межторцовый разрыв) является случайной величиной, распределенной по показательному закону с параметром Известно, что для показательного закона распределения математическое ожидание
(6)
Ло
дисперсия
ОТв = -2 . (7)
Откуда
1о
и
Таким образом, параметр Ло может быть определен экспериментальным путем.
Число стадий распиловки определяется по следующему выражению [2]:
£>Гр м ■ ^
Из (9) следует, что число стадий распиловки бревен не зависит от режимных параметров ЛР, а определяется только геометрическими характеристиками бревен в сортировочной партии.
Чтобы показать влияние геометрических размеров бревен на форму кривой плотности распределения обобщенного закона Эрланга, рассмотрим следующий пример.
Пример. На лесопильных рамах 2Р75-1/2 (мощность двигателя главного привода 110 кВт, частота вращения коленчатого вала 325 об/мин) осуществляется распиловка бревен с брусовкой при выпиливании одного бруса размерами (0,60 ... 0,73)а?в (где йв - вершинный диаметр). Порода древесины - сосна, ель. Число пил в поставе - 8, шаг зубьев - 26 мм. Построим семейство кривых распределения Эрланга для следующих случаев:
1) длина бревен с есть случайная величина с математическим ожиданием
и дисперсиями
2) длины бревен с есть случайные величины с математическим ожиданием и дисперсией
Решение
1. Расчетные данные, необходимые для построения указанных кривых распределения, представлены в табл. 1. Расчетная техническая посылка, ограниченная шероховатостью поверхности пиломатериалов, заполнением впадин зубьев опилками и мощностью электродвигателя механизма резания ЛР согласно [4] Д р = 4 3 м м . В этом случае скорость подачи и = 0, 2 3 3 м/с.
Математическое ожидание длительности вспомогательного времени (межторцового разрыва) примем равным откуда параметр закона
Тогда математические ожидания времени распиловки и рабочего цикла
МЬ 5,1
мт»=-= ТТТТ^Т = 21,89 с;
р и 0,233
МГц = МГр + МТВ = 21,89 + 1,90 = 23,79 с; Дисперсия длительности рабочих циклов
Я 7Ц = Я Гр + Я Гв, (11)
где ЯГв = 3 , 5 6 с2( с м . ф ор мулу ( 8 ) ).
Таблица 1
Дисперси я Я I, м2 Число стадий распиловки к Парамет р я, С"1 Дисперсия, с2
ОГр отц
0,2 130 5,94 3,69 7,25
0,3 87 3,96 5,53 9,09
0,4 65 2,97 7,37 10,93
2. Расчетные данные, необходимые для построения кривых распределения в зависимости от вершинного диаметра (¿в бревен, представлены в табл. 2.
Таблица 2
Диаметр бревен (¿В,СМ Расчетная посылка ДР,мм Скорост ь подачи и, м/с Пар аметр X С"1 Диспер сия отц, с2 Математич еское ожидание МГЦ, с
14 46 0,249 3,18 10,01 24,27
30 34 0,184 2,35 15,37 31,49
36 27 0,146 1,86 22,32 38,67
На рис. 2 и 3 представлены семейства кривых распределения закона Эрланга, рассчитанных по формуле (5).
0,14
10 14 18 22 26 30 34 38 42 46 50 54 Операционное время цикла, с
Рис. 2. Семейство кривых распределения закона Эрланга в зависимости от вершинного диаметра (ёв) бревен: 1 - ёв = 14 см; 2 - 30; 3 - 36 см
10 14 18 22 26 30 34 38 42 Операционное время цикла, с
Рис. 3. Семейство кривых распределения закона Эрланга в зависимости от дисперсии длины бревен (БЬ) диаметром 20 см: 1 - БЬ = 0,2 м2; 2 - 0,3; 3 - 0,4 м2
Согласно центральной предельной теореме [5] при увеличении к обобщенный закон Эрланга (£+1)-го порядка неограниченно приближается к нормальному. Как показала практика моделирования технологических и энергетических процессов деревообработки, вполне
удовлетворительную точность можно получить при аппроксимации законов Эрланга нормальным законом спараметрами М7Ц и D7Ц, вычисленными по формулам (10) и (11), уже при к > 1 5 . При этом ошибка аппроксимации составила не более 0,01 %. Кривые распределения, представленные на рис. 2 и 3, также наглядно подтверждают указанный вывод.
СПИСОК ЛИТЕРАТУРЫ
1. Агеев С.П., Мелехов В.И. Вероятностная модель производственного процесса лесопильного цеха // Актуальные проблемы развития лесного комплекса: материалы Междунар. науч.-техн., конф. 8-10 декабря 2009 г. Вологда, 2010. С. 91-93.
2. Агеев С.П. Математическое моделирование процессов распиловки древесины // СПбЛТА. Вып.179. СПб: ЛТА, 2007. С. 142-152.
3. Агеев С.П. Показатели энергоиспользования лесопильных рам // Лесн. журн. 2009. № 3. С. 131 -138. (Изв. высш. учеб. заведений).
4. Варфоломеев Ю.А., Дружинин И.С., Дьячков Ю.А. Справочник по лесопилению / Под ред. А.М. Копейкина. М.: Экология, 1991. 496 с.
5. Вентцель Е.С., Овчаров Л.А. Теория вероятностей и ее инженерные приложения: учеб. пособие для втузов. М.: Высш. шк., 2000. 480 с.
6. КлейнрокЛ. Теория массового обслуживания. М.: Машиностроение, 1979. 432 с.
7. Римский-КорсаковБ.С. Операционное исчисление. М.: Высш. шк., 1960. 146 с.
Поступила 17.10.12
Stochastic Pattern of the Sash Operating Cycles Ageev S.P., Doctor of Engineering, Associate Professor
Northern (Arctic) Federal University named after M.V. Lomonosov, Naberezhnaya Severnoy Dviny, 17, Arkhangelsk, 163002, Russia E-mail: doctor.mart11@mail.ru
The feature of woodworking is the appearing of forcing of various random factors actually. Therefore, for total analysis of basic regularity of the functioning flow process at woodworking enterprises these processes must be considered as the variety of random processes. One of the most energy-intensive carpentry is the sawmill generating. The sash with relatively small nearby machine mechanization is a separate plot, and as the most powerful electric collector, largely determines the nature of power. The purpose of this article is to establish the probability of timber sawing operational cycles on the sashes.
The mathematical model of loading the cutting device of the main sash drive is offered for carrying out the researches. According to this model, the sash is considered as a system of queueing. For the sawing length analysis of specific logs the step method has been used, designed by A.K. Erlang.
Using the ordered model, the distribution law of the operational cycle period of the sash has been found, the dependence of distribution law from the sash characteristic quantity and the geometric parameters of cutting raw materials has been shown.
The research found that the operational cycles period of the sash is submitted to the generalized Erlang distribution law (£+1) degree as a random quantity. According to the Central limit theorem under the buildup k the generalized Erlang law is closed to normal non-limiting. As the practice of modeling of engineering and energy woodworking process shows, fair accuracy can be find in the approximation of the Erlang laws to the normal law already by k > 15. The average error of the approximation was not more than 0.001%. The distribution curves, presented in the article, also confirm this conclusion.
The obtained results can be used in capability investigation and optimization of the conveyer equipment parameters, as well as in addressing the issues of normalization power consumption of cutting processes.
Keywords: sash, sash operational cycle, queuing system, the density of probability distributions, differential equations, Laplace transform, collapsing theorem, geometrical adjectives of cutting raw materials.
REFERENCES
1. Ageev S.P., Melekhov V.I. Veroyatnostnaya model' proizvodstvennogo protsessa lesopil'nogo tsekha [Probabilistic Model of the Production Process of Timber Shopfloor]. Aktual'nye problemy razvitiya lesnogo kompleksa: mat. mezhdunar. nauchno-tekhn. konf. [Actual Problems of Forestry Complex Development: Math. Intern. Tech. Sci. Conf.]. 8-10 December 2009. Vologda, 2010. pp. 91-93.
2. Ageev S.P. Matematicheskoe modelirovanie protsessov raspilovki drevesiny [Mathematical Modeling of Sawing Wood Processes]. St. Petersburg, 2007, vol. 179, pp. 142-152.
3. Ageev S.P. Analiz elektropotrebleniya ramnykh lesopil'nykh potokov [Analysis of Power Consumption of Frame Sawing Threads]. Sovershenstvovanie energeticheskikh sistem i tekhnologicheskogo oborudovaniya: Sb. nauch. tr. [Development of Energy Systems and Technologic Equipment: Collected Papers]. 2012, pp. 17-23.
4. Varfolomeev Yu.A., Druzhinin I.S., D'yachkov Yu.A., et. al. Spravochnik po lesopileniyu [Sawmilling Guide]. Moscow, 1991. 496 p.
5. Venttsel' E.S., Ovcharov L.A. Teoriya veroyatnostey i ee inzhenernye prilozheniya [Probability Theory and its Engineering Applications]. Moscow, 2000. 480 p.
6. Kleynrok L. Teoriya massovogo obsluzhivaniya [Queueing theory]. Moscow, 1979. 432 p.
7. Rimskiy-Korsakov B.S. Operatsionnoe ischislenie [Operational calculus]. Moscow, 1960. 146 p.