УДК 622.015.002
А.А.ГРАБСКИЙ, канд. техн. наук, профессор, a.a.grabsky@yandex.ru Московский государственный горный университет
А.A.GRABSKY, PhD in eng. sc., professor, a.a.grabsky@yandex.ru Moscow State Mining University
СТАТИСТИЧЕСКАЯ ОЦЕНКА ЭНЕРГЕТИЧЕСКИХ ПОКАЗАТЕЛЕЙ КОМБАЙНА MТS-250 (ГЕРМАНИЯ)
Приведена статистическая оценка и статистическое моделирование работы шнекофре-зерной динамической системы комбайна типа MTS-250 (Германия), осуществляющей управление работой шнека и вибрологического импульсатора. Выполнено статистическое моделирование с использованием разработанной автором цифровой модели всей системы и отдельных ее элементов, что позволяет варьировать рабочими параметрами комбайна и оценить расход энергии за время его работы на объем добытой горной массы, минимизировать рабочие параметры комбайна.
Ключевые слова: шнекофрезерный комбайн, гидромотор, рабочие параметры, цифровая модель, динамика, статистическая оценка, энергетические показатели, производительность.
STATISTICAL EVALUATION OF THE ENERGY RATES OF THE SURFACE MINER MTS-250 (GERMANY)
Statistical estimation and statistical modeling of screw and cutter dynamic system that controls operation of the screw and vibrological impulsator of the surface miner of the type MTS-250 (Germany) is given. The statistical modeling is carried out using a digital model developed by the author as of the whole system and its individual components. This allows to vary the operating parameters of the surface miner and to evaluate the energy input per amount of the mined rock during the time of its work. This will eventually minimize the operating parameters of the surface miner.
Key words, screw and cutter surface miner, hydraulic motor, operating parameters, digital model, dynamics, statistical estimation, energy datum, capacity
Энергопотребление шнекофрезерного комбайна является одним из его основных показателей, а снижение этого показателя -важной технической задачей. В общем случае расход энергии на процесс разрушения горного массива определяется внешним моментом сопротивления, воздействующим на шнекофрезерный исполнительный орган, и угловой скоростью гидромотора:
Е = Ы0югЫ, (1)
где k - некоторый конструктивный коэффициент; юг - угловая скорость вала гидромотора, юг = / dt, рад/с; фг - угол поворо-
та, рад; t - время, за которое определяется расход энергии, ч.
Для снижения внешнего момента сопротивления М0 в комбайне использован специальный виброреологический гидропульсатор, вызывающий вибрацию шнека и позволяющий уменьшить силы трения шнека о горный массив при их взаимодействии.
Экспериментальные и теоретические исследования, выполненные автором [1], показали, что внешний момент сопротивления является случайной величиной и может быть описан как стационарный случайный процесс, представляющий собой сумму нормального случайного процесса, характе-
384 -
ISSN 0135-3500. Записки Горного института. Т.199
Рис. 1. Фрагмент осциллограммы изменения давления в напорной магистрали насосной установки
карьерного комбайна MTS-250
Т1, Т2 - период колебания давления; 1- временная функция изменения момента сопротивления породы при вращении шнека; 2 - усилие, создающее виброреологический эффект при движении шнека
ристики которого зависят от физико-механических свойств горного массива М1, и случайного периодического процесса с равномерно распределенной в интервале [-п, п] фазой, параметры которого определяются настройкой виброреологического гидропульсатора М2 (рис.1).
Характеристики моментов следующие: математическое ожидание первого процесса тм^ , второго тщ = 0, корреляционная
функция первого процесса
Км1 (т) = ^!#1ехР(-а г ^ корреляционная функция второго процесса
км2 (т) = Мо2 адз(ю вт + фо),
где о2М - дисперсия внешнего момента М1,
связанного с напряжением разрушения горного массива £р, Па; аг - константа, характеризующая физико-механические свойства горного массива, с-1; М0, юв - амплитуда и частота виброреологического момента, Н/м, рад/с; ф0 - случайная фаза колебаний, равномерно распределенная в интервале [-п; п].
Закон распределения внешнего момента М1 - нормальный, закон распределения момента М2 - равномерный в интервале[-п; п].
Вероятностные характеристики суммарного момента М0, состоящего из суммы двух статистически независимых величин с
указанными корреляционными функциями и законами распределения, имеют вид [2]:
• корреляционная функция момента
М0 = М1 + М2:
Км0 = Км1 + Км2 = = стМ1 ехр(-ах) + м0! адз(ювТ + Ф0); (2)
• дифференциальный закон распределения (закон распределения В.И.Тихонова) записан следующим образом [2]:
р(м 0) = —^=х
п
х | ехр
0
М в
Мон -Т~ С08 Ф
Лф, (3)
где Мон = Мш / кст^; Мш - суммарный момент; ст^ - среднеквадратическое отклонение момента Мш.
Характер изменения принятого распределения при изменении соотношения между моментами М1 и М2 приведен на рис.2. Как видно, при изменении соотношения а от 0 до 10 закон распределения изменяется от нормального до равномерного, что делает задачу о принятии его устойчивого вида весьма проблематичной.
Второй сомножитель в уравнении (1) -угловая скорость гидромотора юг = Лфг / Л претерпевает различные случайные измене-
385
Санкт-Петербург. 2012
2
р(Мо)
0,5 0,4 0,3 0,2 0,1
^а = 0
= 2 а = 1
^-а = 5
Мон
Рис.2. Характер изменения закона распределения нормированного момента Мон от величины а = Мв1кл^
ния при работе комбайна, которые из-за нелинейности механической системы определить теоретически достаточно сложно.
Для определения статистических характеристик величины / dt была использована разработанная автором цифровая модель [3] механической системы комбайна MTS-250, осуществляющей управление работой шнека и виброреологического гидропульсатора (рис.3).
На этой модели было выполнено статистическое моделирование, заключающееся в исследовании поведения отдельных элементов системы и системы в целом при воздействии на нее полного момента М0, описываемого случайным процессом с характеристиками тм , DM и Км(т). В частности, была
установлена зависимость угловой скорости ф г гидромотора от полного момента М0.
Установлено, что при значительной жесткости привода передаточная функция
W (j<) =
Ф (j<)
M 0(ja)
= const представляет собой
коэффициент передачи в виде некоторого постоянного числа Кп.
В этом случае фг = юг = МшКп и формулу (1) можно записать в виде:
Е = Nt = M„ Kt = МШ Кп Kt
(4)
при этом энергопотребление и мощность N имеют такое же распределение, как и полный момент сопротивления Мш.
Более сложным является случай, когда жесткость гидромотора Сг2<<Сг1 и передаточная функция не является константой. В этом случае величина Е была определена экспериментально с использованием дополнительно разработанной модели (рис.4).
В модели для расчета потребляемой энергии имеются два входа, на которые подаются значения момента М0 (вход 1), и скорость юг (вход 2). Эти величины формируются в общей модели (см.рис.3, точки А и В). Далее эти величины поступают на интеграторы, где они интегрируются на определенном интервале времени Дt (час, смена,
Момент 1 dFd / dt
М dFd / dt
Fm 3 Fd
dFm / dt
U
Дизель
и 5
dFd /dt
dFm / dt
Рис.3. Схема общей модели рабочей системы комбайна для анализа динамических процессов
1 - блок формирования динамического момента дизеля; 2 - модель функционирования гидромотора: 3 - модель дизеля; 4 - гидропульсатор; 5 - блок управления гидропульсатором; 6 - блок формирования момента сопротивления
ISSN 0135-3500. Записки Горного института. Т.199
0
2
4
5
0
Рис.4. Модель для расчета усредненной мощности N и усредненного расхода энергии Е
сутки и т.д.) и делятся на это время, т.е. находятся усредненные среднеинтегральные на интервале Дt величины момента м0 и скорости юг = Лфг / Л :
М 0 =
'At
J M 0(t )dt
'At
Jp r(t )dt
0
/ At.
/ At.
Затем усредненные величины перемножаются, учитываются коэффициенты к и получается усредненная на интервале Лt мощность комбайна
N = M,„ ~ k .
(5)
Если не осуществлять деления на временной интервал Л^ то можно получить расход энергии (за час, смену, сутки):
Лt
Е(г) = к |мш(г)шг(г)Лг.
0
Если разделить расход энергии Е за время Т на объем горной массы V = QT, где Q - объемная производительность комбайна, м3/ч, то можно получить удельный показатель, равный
=
е = От '
т.е. расход энергии на кубометр произведенной продукции. В зависимости от того,
как получен расход энергии, этот показатель может быть усредненным или вероятностным.
Варьируя параметры комбайна, можно добиться минимизации этих показателей.
ЛИТЕРАТУРА
1. Грабский А.А. Анализ динамического нагруже-ния шнекофрезерного рабочего органа карьерного комбайна статистическим методом // Горная промышленность. 2012. № 4 (104). С.158-160.
2. Грабский А.А. Установление дифференциального закона распределения полного момента сопротивления на шнекофрезерном рабочем органе карьерного комбайна // Горная промышленность. 2012. № 5 (105). С.91-95.
3. Грабский А.А. Использование цифрового моделирования при анализе устойчивости динамической системы карьерного комбайна со шнекофрезерным рабочим органом типа MTS (Германия) // Научно-технический вестник Поволжья. 2012. № 6. С.204-211.
REFERENCES
1. Grabsky А.А. Analysis of dynamic loading of the screw and cutter working body of a surface miner by the statistical method // Mining Industry. 2012. № 4 (104). P.158-160.
2. Grabsky А.А. Establishment of the differential distribution law of the whole resistance moment on the screw and cutter working body of a surface miner // Mining Industry. 2012. № 5 (105). P.91-95.
3. Grabsky А.А. The use of digital simulation in the stability analysis of the dynamical system of a surface miner with screw and cutter working body of the type MTS (Germany) // Volga Region Scientific-technical Newsletter. 2012. № 6. P.204-211.
®r =
-387
Санкт-Петербург. 2012