ГИАБ. Горный информационно-аналитический бюллетень / MIAB. Mining Informational and Analytical Bulletin, 2020;(7):94-103 ОРИГИНАЛЬНАЯ СТАТЬЯ / ORIGINAL PAPER
УДК 622.23.05:622.235 DOI: 10.25018/0236-1493-2020-7-0-94-103
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ВЗАИМОДЕЙСТВИЯ БОЙКА И ИНСТРУМЕНТА ПРИ УДАРНОМ РАЗРУШЕНИИ ГОРНЫХ ПОРОД
А.Б. Жабин1, И.М. Лавит1, А.В. Поляков1, З.Э. Керимов1
1 Тульский государственный университет, Тула, Россия, e-mail: zhabin.tula@mail.ru
Аннотация: Обоснована расчетная схема, в которой ударная система машины моделируется продольно нагруженными упругими цилиндрическими стержнями при малых деформациях. Горный массив рассматривается как абсолютно твердая преграда. Боек и инструмент характеризуются длиной, площадью поперечного сечения, плотностью материала и модулем Юнга. Модель позволяет найти усилие (как функцию времени), с которым инструмент действует на массив. Считается, что перемещения и скорости стержней в поперечном направлении пренебрежимо малы по сравнению с продольным направлением, а внешние силы, отличные от реакций связей, на систему стержней не действуют. Составлено вариационное уравнение, выражающее принцип возможных перемещений. Вариации от времени не зависят. Рассмотрены начальные и граничные условия. Вариационное уравнение решалось методом прямых, при котором оператор дифференцирования по времени заменяется конечно-разностным оператором. Задача сводится к последовательно решаемым краевым задачам с переменной правой частью. Используемая конечно-разностная схема — это апробированная неявная схема Кранка-Николсон. Краевые задачи на каждом шаге интегрирования по времени решаются методом конечных элементов. В результате вариационное уравнение преобразуется к системе линейных алгебраических уравнений, приведенное решение которой определяет искомое усилие. Представлен пример расчета для гидравлической машины ударного действия G100 фирмы «Rammer» (Финляндия) в виде графической зависимости силы давления инструмента на массив от времени с шагом 0,1 мкс, а также его анализ. Относительная погрешность вычислений продолжительности удара и максимального (по модулю) значения силы не превосходит 0,1%.
Ключевые слова: горный массив, боек, ударный инструмент, математическая модель, усилие, скорость, перемещения, начальные и граничные условия, вариационное уравнение, оператор, матрица, время.
Благодарность: Результаты исследования опубликованы при финансовой поддержке ТулГУ в рамках научного проекта № НИР 2019-18.
Для цитирования: Жабин А.Б., Лавит И. М., Поляков А. В., Керимов З. Э. Математическая модель взаимодействия бойка и инструмента при ударном разрушении горных пород // Горный информационно-аналитический бюллетень. - 2020. - № 7. - С. 94-103. DOI: 10.25018/0236-1493-2020-7-0-94-103.
Mathematical model of piston-tool interaction in rock fracture by impact
A.B. Zhabin1, I.M. Lavit1, A.V. Polyakov1, Z.E. Kerimov1
1 Tula State University, 300012, Tula, Russia, e-mail: zhabin.tula@mail.ru
© А.Б. Жабин, И.М. Лавит, А.В. Поляков, З.Э. Керимов. 2020.
Abstract: The authors justify a computation model of a machine percussion system simulated by elastic cylindrical rods subjected to maximal axial load at minimal strain. Rock mass is assumed as a perfectly solid block. The piston and tool are described by the values of length, cross-section area, density and Young's modulus. The model determines the force applied by the tool on the rock as function of time. It is assumed that transverse displacements and velocities of the rods are negligeable as compared with the axial displacements and velocities, while the rods are free from the action of the external forces different from the restraining forces. The variational equation expresses the principle of possible displacements. The variations are independent of time. The initial and boundary conditions are considered. The variational equation is solved using the method of straight lines, with replacement of a time differentiation operator by the finite difference operator. The problem reduces to the successive solving of boundary value problems with variable right-hand sides. The finite difference scheme is the approved implicit scheme of Crank-Nicolson. The boundary value problems are solved using the finite element method at each step of integrating. As a result, the variational equation transforms into a system of linear algebraic equations, and the reduced solution of this system yields the wanted force. The calculations are illustrated by the tool press force-time curve plotted with a step of 0.1 ps for hydropercussion machine G100 by Rammer, Finland. The relative calculation error of the impact duration and maximal force (in absolute magnitude) is not higher than 0.1%.
Key words: rock, piston, percussive tool, mathematical model, force, velocity, displacements, initial and boundary conditions, variational equation, operator, matrix, time.
Acknowledgements: The studies were supported by the Tula State University, R&D Project No. 2019-18.
For citation: Zhabin A. B., Lavit I. M., Polyakov A. V., Kerimov Z. E. Mathematical model of piston-tool interaction in rock fracture by impact. MIAB. Mining Inf. Anal. Bull. 2020;(7):94-103. [In Russ]. DOI: 10.25018/0236-1493-2020-7-0-94-103.
Введение
Одним из перспективных способов разрушения горных пород повышенной прочности является механический удар (с использованием гидроударников), обеспечивающий высокую концентрацию нагрузки на локальном участке забоя и уступающий по энергоемкости разрушения (17 — 25 Дж/см2) только взрыву (6 Дж/см2) [1 — 8].
Ударная система технологической машины представляется цилиндрическими
стержнями, соотношения длин и сечений которых зависят от конструктивно-технологических особенностей исполнительного органа [8 — 12]. В соответствии с этим в наиболее авторитетной работе [8] для моделирования процесса энергопередачи рассмотрено соударение двух стержней, являющихся бойком 1 и волноводом (инструментом) 2. Разрушаемый массив также представлен цилиндрическим стержнем 3 со свойствами, характерными для горной по-
Рис. 1. Расчетная схема ударной системы Fig. 1. Computational model of percussion system
A_1 В С_2 Г)_3 F
. _--------x----------------------------------------------------1-►
___1 -x
роды, и имеющим массу т3 >> т1 + т2 (рис. 1).
В работе [8] приведены дифференциальные уравнения и граничные условия, а также качественный анализ решения, заключающийся в представлении волнового характера процесса, перечислении и обсуждении различных типов волн, распространяющихся в стержнях 1, 2 и 3 (см. рис. 1). Однако ни метода решения задачи, ни результатов расчетов в работе [8] нет. Отмечается только, что представленная качественная картина изменения напряжений в элементах ударной системы характеризует основные закономерности передачи энергии удара в массив.
В настоящей работе предпринята попытка дать количественный анализ процесса ударной передачи энергии в массив горной породы. Естественно принять убедительно обоснованную расчетную схему работы [8], в которой ударная система машины моделируется цилиндрическими стержнями. Однако от моделирования горного массива стержнем стоит отказаться. Это следует сделать по двум причинам.
Первая — это то, что геометрически массив даже отдаленно не напоминает стержень. Если уж его моделировать упругой или упругопластической средой, следует выбрать более реалистическую конфигурацию. Но это чрезмерно усложняет задачу.
Вторая причина — это то, что моделирование деформирования массива на основе какой-либо модели деформируемого твердого тела всегда будет неадекватным. Дело в том, что крепкие горные породы деформируются упруго вплоть до разрушения. Но, с другой стороны, в узкой зоне, прилегающей к внедряемому инструменту, образуется из частиц разрушенного материала породы так называемое ядро сжатия. Это ядро обусловливает нелинейность зависимо-
сти сила-перемещение инструмента. Получить эту зависимость можно только экспериментально, но к настоящему времени какие-либо опытные данные отсутствуют.
С учетом сказанного, ради получения более простого и однозначного решения массив следует рассматривать как абсолютно твердую преграду. Это, конечно, внесет погрешность в расчет, так как уменьшится время удара. Однако главная величина, ради которой, собственно, и следует делать расчет — импульс, передаваемый породному массиву, — останется практически неизменным.
Таким образом,исследования в этом направлении дополнят монографию [1] в решении поставленной в ней задачи.
Расчетная схема
Несмотря на разнообразие конструкций ударных инструментов [2 — 12], принцип их действия одинаков: боек разгоняется до большой скорости, затем он ударяет по скалывающему инструменту, прижатому к разрушаемому массиву. Эффективность разрушения обусловлена динамичностью процесса разрушения. Расчетная схема работы инструмента представлена на рис. 2.
При этом инструмент может отличаться конструктивными параметрами, разрушать уступ различной высоты и располагаться под любым углом к нему. Боек также может иметь конструктивные отличия. Рассматривать взаимодействие всех трех объектов, изображенных на рис. 2, нет необходимости.
Горный массив 3 (см. рис. 2.) влияет на процесс взаимодействия бойка 1 и инструмента 2 только как ограничитель перемещений инструмента, подверженного действию бойка. Эти перемещения, конечно, отличны от нуля в силу упругости массива и возможного разрушения его материала непосредственно под инструментом. Однако они
у// //////////////////////// /.
У //У /У/ ////// /////////////// ¿г,
Рис. 2. Расчетная схема: 1 — боек; 2 — ударный инструмент; 3 — горный массив; 4 — скалываемый уступ
Fig. 2. Computational model: 1 — piston; 2 — percussive tool; 3—rock mass; 4 — cut bench
малы и ясно, что без заметной погрешности их можно считать нулевыми. При этом математическая модель распадается на две модели.
Первая из них — это модель взаимодействия бойка 1 (см. рис. 2) и инструмента 2 при условии, что рабочий конец инструмента, контактирующий с массивом, неподвижен. Эта модель позволяет найти усилие (как функцию времени), с которым инструмент действует на массив.
Вторая математическая модель должна дать описание напряженного состояния массива в зависимости от времени. При включении в нее критерия прочности она позволит найти объем разрушаемого материала.
Объединение упомянутых моделей даст возможность установить влияние различных конструктивных и режимных параметров, а также прочности горного массива на объем разрушенного материала. В данной статье рассматривается первая математическая модель.
Согласно работе [8] боек и инструмент можно рассматривать как продольно нагруженные упругие стержни при малых деформациях. Будем считать, что координатная ось х совпадает с осями обоих взаимодействующих стержней (рис. 3).
Боек (стержень 1) характеризуется следующими величинами: l — длиной; F1 — площадью поперечного сечения; р1 — плотностью материала; E1 — модулем Юнга. Ударный инструмент характеризуется тем же набором параметров (они имеют индекс 2).
Для стержня перемещения и скорости в поперечном направлении пренебрежимо малы по сравнению с перемещениями и скоростями в продольном направлении. Учтем также то, что внешние силы, отличные от реакций связей, на систему стержней 1 и 2 (см. рис. 3) не действуют.
Пусть t — время. Введем в рассмотрение поле перемещений u(t, x) и поле скоростей v = U = ди/8t. В данном слу-
U
I,
Рис. 3. Расчетная схема взаимодействия бойка и инструмента: 1 — боек; 2 — инструмент Fig. 3. Computational model of piston-tool interaction: 1 — piston; 2 — tool
чае — при продольной деформации стержня — из всех компонент тензора напряжений отлична от нуля только одна. Обозначим ее через а. Согласно закону Гука связь напряжения с продольной деформацией е = ди/дх дается формулой а = Ег.
Вариационная постановка задачи
Так как задача решается численно, исходное соотношение удобно записать в вариационной формулировке. Вариационное уравнение, выражающее принцип возможных перемещений, применительно к данной задаче имеет вид
¡Г^р^-5" + ст5е^dx = 0 , (1)
где символом 5 обозначены вариации соответствующих величин. Вариации — это произвольные функции координат, ограниченные лишь одним требованием: на тех участках граничной поверхности, где заданы главные (кинематические) граничные условия (в данном случае при х = 0), вариации перемещений равны нулю. Вариации от времени не зависят.
Рассмотрим начальные и граничные условия задачи. В начальный момент Г = 0 стержни составляют один стержень длиной I = I + 12 (см. рис. 3). Перемещения, деформации и напряжения в начальный момент равны нулю:
t = 0; х е [0, I]: и = 0;
г = 0; а = 0. (2)
Скорости точек инструмента в начальный момент равны нулю, скорости точек бойка — равны (-V). Это означает, что при Г = 0 боек ударяет по инструменту:
ГX е[0,/2): V = 0;
'=0: \ (3)
Гх е(/2, /]: V = -У.
При х = 0 задается граничное усло-
Это условие является главным (кинематическим) условием. При х = I задается естественное (динамическое) граничное условие:
х = I: а = 0. (4)
В процессе удара могут быть изменения в граничных условиях в зависимости от значения напряжения в точке х = 12. Если в этой точке а < 0 (сжатие), то это значит, что контакт между стержнями имеет место. Если же получается а £ 0, то это значит, что контакт отсутствует, то есть начинается отскок бойка. С этого момента рассматривается только стержень 2, то есть х е [0, 12] (см. рис. 3) при граничном условии
I,
0.
(5)
вие:
х = 0: и = 0; у = 0.
Таким образом, входящие в уравнение (1) величины I, р, Р, а также модуль Юнга Е определяются выражениями:
1. Стержни находятся в контакте:
/ = к + ^
х е(0,/2): Р = р = р2; Е = Е2;
х б(12,1): р = р1; р = р1; Е = Е1
2. Стержень 1 отскочил от стержня 2, но при этом стержень 2 прижат к массиву:
I = 12;
X е(0,I): F = F2; р = р2; Е = Е2
Время удара определяется значением напряжения в точке х = 0. Если выполняется условие
х = 0 : а < 0, (6)
то это значит, что стержень 2 (инструмент) прижат к массиву. Как только оно нарушится, удар завершается: инструмент отскакивает от массива.
Зависимость а(Г) при х = 0 является основным результатом расчета. Величина
Р^) = )| х (7)
представляет собой усилие, с которым ударный инструмент действует на горный массив.
Метод решения задачи
Для решения уравнения (1) применяем метод прямых [13]. Суть его состоит в том, что оператор дифференцирования по времени заменяется конечно-разностным оператором. Задача сводится к последовательно решаемым краевым задачам с переменной правой частью. Используемая в данной работе конечноразностная схема — это апробированная неявная схема Кранка-Ни-колсон [14].
Пусть Дг — шаг интегрирования по времени, а п — номер шага интегрирования. Конечноразностное представление производной по времени от некоторой функции у(^х) на п-м шаге имеет вид
(8)
ду дt
Дt
|F Д( -+ 1
Д
ёх = 0,
1 п , л-1\
Г
4р .(* )2
и"5и + ап 5е
\ г А (""+^1) -
dx =
15е
dx
где у"-1 и у" — значения у на границах временного интервала. Величины, не содержащие производных по времени, представляются на п-м шаге интегрирования по времени в виде
У = 1 {у" + У"-1), (9)
то есть как среднее арифметическое граничных значений. С использованием формул (8) и (9) конечно-разностное представление уравнения (1) записывается в виде
Д 2х (10)
Величины с индексами п-1 известны из решения для предыдущего шага. Из системы (10) определяются величины с индексом п. Выразим из второго уравнения системы (10) величину у" и подставим в первое уравнение:
= — (ип - ип-1)- ; (11)
Дг4 '
Таким образом, система (10) распадается на два последовательно решаемых уравнения. Вначале из уравнения (12) определяются перемещения и", а затем из уравнения (11) — скорости у".
Краевые задачи [13] на каждом шаге интегрирования по времени решаются методом конечных элементов [14,15]. В результате вариационное уравнение
(12) преобразуется к системе линейных алгебраических уравнений:
Аип = В (13)
где и" — матрица-столбец узловых перемещений; В — матрица-столбец правых частей (напомним, что узловые перемещения и"-1 известны); А — квадратная матрица коэффициентов системы уравнений. Матрицы А и В необходимо преобразовать с целью учесть главное граничное условие. Это делается следующим образом. Пусть перемещение к-го узла (в рассматриваемом случае к = 1) равно нулю. Обнуляется к-я строка и к-й столбец матрицы А за исключением диагонального элемента. Обнуляется также к-й элемент матрицы В.
После того как система уравнений
(13) решена, необходимо найти реакцию связи в узле 1. Это сила, с которой инструмент воздействует на горный массив. Запишем первую строку матричного равенства (13)
АД" = в, + я. (14)
Фигурирующие здесь элементы матриц А и В — это элементы непреоб-разованных матриц. Поэтому равенство возможно только при добавлении к правой части величины Я. Это и есть реакция связи — сила, с которой массив действует на ударный инструмент. По
Рис. 4. Зависимость силы, действующей на массив, от времени Fig. 4. Time curve of force applied to rock mass
третьему закону Ньютона инструмент действует на массив с силой минус R.
Пример расчета
В качестве примера расчета определим зависимость (7) для гидравлической машины ударного действия G100 фирмы «Rammer» (Финляндия) [8]. И боек, и инструмент изготавливаются из стали. Поэтому E = E1 = E2 = 2 • 1011 Н/ м2; р = р1 = р2 = 7800 кг/м3. Длины и площади поперечных сечений стержней: l1 = 2,9 м, l2 = 1,45 м, F1 = 0,168 м2 и F2 = 0,0227 м2. Скорость удара — & = 2,18 м/с. Результаты расчета, представленные ниже, получены при шаге по времени At = 0,1 мкс. И боек, и инструмент разбивались на 500 конечных элементов. При этом, как показали сравнительные расчеты, относительная погрешность вычислений продолжительности удара и максимального (по модулю) значения силы P не превосходит 0,1%.
График зависимости силы давления инструмента на массив от времени представлен на рис. 4. Ступенчатый характер зависимости обусловлен распрост-
ранением упругих волн вдоль стержней. Скорость звука равна с = у]Е/ р = = 5060 м/с, время прохождения волны вдоль инструмента получается равным 12/с = 0,286 мс, что согласуется с результатами, приведенными на рис. 4.
Остроугольные «выскоки» амплитуды в зонах скачков функции, обусловленных приходом и отражением волн, представляют собой известное явление Гиббса [16], возникающее всегда при аппроксимации разрывных функций непрерывными. Они не вносят принципиальных изменений в результаты расчета.
Через 5,73 мс боек отскакивает от инструмента, но воздействие инструмента на массив продолжается. Оно закачивается при 6,01 мс. Это и есть время удара. Максимальное значение силы давления инструмента на массив можно использовать для анализа прочности инструмента. В данном случае Ртах = 12,9 МН. При этом напряжение (сжимающее) в инструменте будет равно 568 Н/мм2, что приблизительно в два раза меньше предела текучести типичной высокопрочной стали.
Заключение
При математическом моделировании взаимодействия бойка и инструмента необходимо рассматривать их как продольно нагруженные стержни, характеризующиеся длиной, площадью попе-
речного сечения, плотностью и модулем Юнга.
При неподвижной рабочей поверхности инструмента это позволяет определять усилие, с которым он действует на массив, как функцию времени.
СПИСОК ЛИТЕРАТУРЫ
1. Кравченко В.А. Обоснование гидравлических устройств ударного действия с пониженной удельной металлоемкостью для разрушения горных пород. Автореф. дис. канд. техн. наук. — Орел, 2004. — 23 с.
2. Доброборский Б. С., Овчаров А. А. Исследования машины ударного действия с подвижным ударным механизмом // Современные проблемы науки и образования. — 2014. — № 5. URL: https://science-education.ru/pdf/2014/ 5/ 468.pdf (дата обращения: 24.09.2014).
3. Абраменков Д. Э., Абраменков Э. А., Дедов А. С., Крутиков Е. И. Пневматический ударный механизм с комбинированным воздухораспределением // Известия высших учебных заведений. Строительство. —2014. — № 4. — С. 114 — 118.
4. Gorodilov L. V. Analysis of the dynamics of two-way hydropercussion systems. Part II: Influence of design factors and their interaction with rocks // Journal of Mining Science. 2013. Vol. 49. No 3. Pp. 465 — 474. DOI: 10.1134/S1062739149030143.
5. Oparin V. N., Timonin V. V., Karpov V. N., Smolyanitsky B. N. Energy-based volumetric rock destruction criterion in the rotary-percussion drilling technology improvement // Journal of Mining Science. 2017. Vol. 53. No 6. Pp. 1043 — 1064. DOI: 10.1134/S1062739117063114.
6. Redelin R.A., Kravchenko V.A., Kamanin Y. N., Panichkin A. V. and Bozhanov A. A. Study of effect of in-line hydropneumatic accumulators on output characteristics of hydraulic hammer // IOP Conference Series Earth and Environmental Science. 2017. Vol. 87. No 2. Article 022016. DOI: 10.1088/1755-1315/87/2/022016.
7. Tian J., Yang Z, Li Y., Yang L. Vibration analysis of new drill string system with hydro-oscillator in horizontal wellt // Journal of Mechanical Science and Technology. 2016. Vol. 30. No 6. Pp. 2443 — 2451. DOI 10.1007/s12206-016-0504-z.
8. Ушаков Л. С., Котылев Ю. Е., Кравченко В. А. Гидравлические машины ударного действия. — М.: Машиностроение, 2000. — 416 с.
9. Жуков И. А. Исходные основания к изучению влияния форм бойков на форму ударного импульса в машинах ударного действия // Вестник Кузбасского государственного технического университета. — 2014. — № 5 (105). — С. 25 — 27.
10. Еремьянц В. Э. К вопросу о рациональной форме бойков ударно-вращательных бурильных машин // Физико-технические проблемы разработки полезных ископаемых. — 2011. — № 5. — С. 74 — 82.
11. Авдеева А. И., Шадрина А. В., Саруев Л. А. Разработка методики расчета распространения продольных волн деформаций по буровому инструменту // Горный информационно-аналитический бюллетень. — 2014. — № 5. — С. 328 — 333.
12. Ветюков М. М., Юнгмейстер Д. А., Пивнев В. А. и др. Теоретические расчеты параметров ударных систем «поршень-боек-штанга» // Горный информационно-аналитический бюллетень. — 2010. — № 7. — С. 329 — 330.
13. Ладыженская О. А. Краевые задачи математической физики. — М.: Наука, 1973. — 408 с.
14. Стренч Г., ФиксДж. Теория метода конечных элементов. — М.: Мир, 1977. — 349 с.
15. Зенкевич О. Метод конечных элементов в технике. — М.: Мир, 1975. — 541 с.
16. Джексон Д. Ряды Фурье и ортогональные полиномы. — М.: ГИИЛ, 1948. —260 с. fj^
REFERENCES
1. Kravchenko V. A. Obosnovaniegidravlicheskikh ustroystv udarnogo deystviya s ponizhen-noy udelnoy metalloemkost'yu dlya razrusheniya gornykh porod [Feasibility study of hydraulic percussion machines of reduced steel intensity for rock fracture], Candidate's thesis, Orel, 2004, 23 p.
2. Dobroborskiy B. S., Ovcharov A. A. Analyses of percussive machines with mobile hammering mechanism. Sovremennye problemy nauki i obrazovaniya. 2014, no 5, available at: https://science-education.ru/pdf/2014/ 5/ 468.pdf (accessed 24.09.2014).
3. Abramenkov D. E., Abramenkov E. A., Dedov A. S., Krutikov E. I. Air-driven percussion mechanism with composite air distribution. Izvestiya vysshikh uchebnykh zavedeniy. Stroitelstvo. 2014, no 4, pp. 114-118. [In Russ].
4. Gorodilov L. V. Analysis of the dynamics of two-way hydropercussion systems. Part II: Influence of design factors and their interaction with rocks. Journal of Mining Science. 2013. Vol. 49. No 3. Pp. 465-474. DOI: 10.1134/S1062739149030143.
5. Oparin V. N., Timonin V. V., Karpov V. N., Smolyanitsky B. N. Energy-based volumetric rock destruction criterion in the rotary-percussion drilling technology improvement. Journal of Mining Science. 2017. Vol. 53. No 6. Pp. 1043-1064. DOI: 10.1134/S1062739117063114.
6. Redelin R. A., Kravchenko V. A., Kamanin Y. N., Panichkin A. V. and Bozhanov A. A. Study of effect of in-line hydropneumatic accumulators on output characteristics of hydraulic hammer. IOP Conference Series Earth and Environmental Science. 2017. Vol. 87. No 2. Article 022016. DOI: 10.1088/1755-1315/87/2/022016.
7. Tian J., Yang Z., Li Y., Yang L. Vibration analysis of new drill string system with hydro-oscillator in horizontal wellt. Journal of Mechanical Science and Technology. 2016. Vol. 30. No 6. Pp. 2443-2451. DOI 10.1007/s12206-016-0504-z.
8. Ushakov L. S., Kotylev Yu. E., Kravchenko V. A. Gidravlicheskie mashiny udarnogo deystviya [Hydropercussive machines], Moscow, Mashinostroenie, 2000, 416 p.
9. Zhukov I. A. Starting framework for the analysis of shape effect of pistons on impact pulse form in percussive machines. Vestnik Kuzbasskogo gosudarstvennogo tekhnicheskogo univer-siteta. 2014, no 5 (105), pp. 25-27. [In Russ].
10. Erem'yants V. E. Efficient shape of pistons in rotary-percussive machines. Fiziko-tekh-nicheskie problemy razrabotki poleznykh iskopaemykh. 2011, no 5, pp. 74-82. [In Russ].
11. Avdeeva A. I., SHadrina A. V., Saruev L. A. Calculation procedure of P-wave propagation along a drilling tool . MIAB. Mining Inf. Anal. Bull. 2014, no 5, pp. 328-333. [In Russ].
12. Vetyukov M. M., Yungmeyster D. A., Pivnev V. A. Theoretical design of percussive piston-tool-rod systems. MIAB. Mining Inf. Anal. Bull. 2010, no 7, pp. 329-330. [In Russ].
13. Ladyzhenskaya O. A. Kraevye zadachi matematicheskoy fiziki [Boundary values problems in mathematical physics], Moscow, Nauka, 1973, 408 p.
14. Strang G., Fix J. Teoriya metoda konechnykh elementov [Finite element method theory], Moscow, Mir, 1977, 349 p.
15. Zenkevich O. Metod konechnykh elementov v tekhnike [Finite element method in engineering], Moscow, Mir, 1975, 541 p.
16. Jackson D. Ryady Fur'e i ortogonal'nye polinomy [Fourier's series and orthogonal polynomials], Moscow, GIIL, 1948, 260 p.
ИНФОРМАЦИЯ ОБ АВТОРАХ
Жабин Александр Борисович1 - д-р техн. наук, профессор, e-mail: zhabin.tula@mail.ru, Лавит Игорь Михайлович1 - д-р физ.-мат. наук, профессор, e-mail: igorlavit@yandex.ru,
Поляков Андрей Вячеславович1 - д-р техн. наук, профессор, e-mail: Polyakoff-an@mail.ru,
Керимов Захит Эльдарович1 - аспирант, e-mail: k-zahit94@mail.ru,
1 Тульский государственный университет.
Для контактов: Жабин А.Б., e-mail: zhabin.tula@mail.ru.
INFORMATION ABOUT THE AUTHORS
A.B. Zhabin1, Dr. Sci. (Eng.), Professor, e-mail: zhabin.tula@mail.ru, I.M. Lavit 1, Dr. Sci. (Phys. Mathem.), Professor, e-mail: igorlavit@yandex.ru,
A.V. Polyakov1, Dr. Sci. (Eng.), Professor, e-mail: Polyakoff-an@mail.ru, Z.E. Kerimov1, Graduate Student, e-mail: k-zahit94@mail.ru, 1 Tula State University, 300012, Tula, Russia. Corresponding author: A.B. Zhabin, e-mail: zhabin.tula@mail.ru.
Получена редакцией 26.03.2020; получена после рецензии 21.04.2020; принята к печати 20.06.2020. Received by the editors 26.03.2020; received after the review 21.04.2020; accepted for printing 20.06.2020.
A_
РУКОПИСИ, ДЕПОНИРОВАННЫЕ В ИЗДАТЕЛЬСТВЕ «ГОРНАЯ КНИГА»
ИНСТРУМЕНТАРИЙ ПРОГРАММНО-АНАЛИТИЧЕСКОГО УПРАВЛЕНИЯ ЭЛЕКТРОПОТРЕБЛЕНИЕМ НА ГОРНЫХ ПРЕДПРИЯТИЯХ
(№ 1222/07-20 от 07.05.2020; 8 с.) Карпенко Сергей Михайлович1 — канд. техн. наук, доцент, Безгинов Глеб Юзович1 — аспирант, 1 ГИ НИТУ «МИСиС».
Приведен анализ существующих способов контроля энергоресурсов, а также математических методов прогнозирования электропотребления на горных предприятиях. Проведен обзор современных программных продуктов и диспетчерских систем, используемых для обеспечения сбора и обработки информации, автоматизации, контроля и мониторинга различных параметров обследуемого объекта. Предложено создание на этой основе программно-аналитических комплексов для совершенствования управления электропотреблением и энергосбережением на горных предприятиях.
Ключевые слова: горные предприятия, электропотребление, контроль, моделирование, прогнозирование, инструментарий, программный пакет, программно-аналитическое управление.
TOOLS FOR SOFTWARE AND ANALYTICAL MANAGEMENT OF POWER CONSUMPTION AT MINING ENTERPRISES
S.M. Karpenko1, Cand. Sci. (Eng.), Assistant Professor, G.Yu. Bezginov1, Graduate Student,
1 Mining Institute, National University of Science and Technology «MISiS», 119049, Moscow, Russia.
The article presents an analysis of existing methods for controlling energy resources, as well as mathematical methods for predicting power consumption at mining enterprises. A review of modern software products and dispatcher systems used to ensure the collection and processing of information, automation, control and monitoring of various parameters of the object under investigation. It is proposed to create software and analytical complexes on this basis to improve the management of power consumption and energy conservation at mining enterprises.
Key words: mining enterprise, power consumption, control, modeling, forecasting, tools, software package, software and analytical management.