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

УДК 622.742

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ РАЗДЕЛЕНИЯ ГОРНЫХ ПОРОД В СЕПАРАТОРЕ ПО ТРЕНИЮ И УПРУГОСТИ

Ляпцев С. А., Потапов В. Я., Потапов В. В., Семериков Л. А., Беридзе Е. Т.

Представлена конструкция сепаратора по трению и упругости. Приведены математические уравнения, описывающие поведения частиц горной массы на наклонной поверхности узла стратификации и повторяющихся ударах свободно падающих частиц о наклонные поверхности отражающих элементов.

Ключевые слова: сепаратор; дифференциальные уравнения; математическое моделирование; зависимости.

Сепаратор СПРУТ [1] предназначен для дачу материала в зону разделения, но и подразделения сыпучих материалов на основе готавливает материал с различными коэф-различий в коэффициентах трения и восста- фициентами трения к разделению. После новления горных пород. Сепаратор включает: прохождения по лотку сформированный

корпус 1 с загрузочным лотком в виде наклон- поток продуктов в виде веера подается на ной плоскости с трамплином 2 и отражатель- ярусно расположенные отражательные эле-ные элементы 3, закрепленные консольно в менты, установленные таким образом, чтобы

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

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

Процесс движения каждой частицы возможно описать математической моделью, включающей уравнения движения на каждом этапе разделения и дифференциальные уравнения движения частицы в циркулирующем потоке воздуха [2, 3, 4].

Свободный полет частицы начинается из точки В со скоростью V направленной по Разделение сыпучих материалов в се- касательной к дуге окружности трамплина параторе СПРУТ осуществляется поэтапно: в данной точке. Если частица при движении исходный продукт подается на поверхность коснется поверхности отражательных эле-лотка, который обеспечивает не только по- ментов, то произойдет удар.

раме 4 (рис. 1).

Рис. 1. Схема аппарата СПРУТ:

1 - корпус сепаратора; 2 - наклонный желоб; 3 - отражательные элементы; 4 - рама

№ 4(36), 2014

33

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

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

Множество вариантов движения частицы при различных значениях исходных параметров убеждает в необходимости проведения математического эксперимента на ЭВМ. Это позволит не только предсказать поведение в сепараторе частиц с различным содержанием полезного компонента, но и подобрать наиболее рациональные конструктивные параметры самого сепаратора.

Рассмотрим движение частицы по углу стратификации. Получив скорость, равную V частица выходит на криволинейный участок дуги AB (рис. 2).

Рис. 2. Силы, действующие на частицу при движении по дуге окружности

Для определения скорости частицы VB в точке В составим дифференциальные уравнения движения частицы в проекциях на оси естественной системы координат (т, n):

Г dV

m • — = mg • sm9 - F dt

<

V2

m----= N - mg • cos9. (1)

l r

Решая систему (1), получаем дифференциальное уравнение:

2 иф + fU = g (ыпф - f (2)

Общее решение дифференциального уравнения (2) имеет вид:

jr г -2/ф 2g(2 f2 +1)

U = Ce JT \------еоэф +

4/2 +1

2 gf

+—¥— БШф, (3)

4 /2 +1

где С = 2ge2/p[Z/r(sin в - /cos в) - [/sin в - (2/ + + 1)cos в]/(4/ + 1)].

Скорость частицы в конце криволинейного участка:

=2 g^(sin p- f„ • cosP) •e2 +

+ 42Г X X f • [sin(P - Y) - e2/ск sinP] +

+(2 fL +1) •[e2/ск cos p - cos(p - Y)]}, (4)

где у - угол наклонной касательной к данной криволинейной поверхности в точке В; R1 - её радиус, м.

На следующем участке частицы выходят на выпуклую криволинейную поверхность -дефлектор (рис. 3). На подвижную частицу действуют сила тяжести G и реакция криволинейной поверхности (нормальная составляющая N, касательная - F ).

Уравнения движения частиц по дефлектору в естественной форме (в проекциях на оси т и n) имеют вид:

Гтат = -mgcos ф - F^

[man = - N + mgsin ф, (5)

34

Известия Уральского государственного горного университета

где ат = dVT / dt = Rd2ф/dt2 - касательное ускорение, an = VT2 / R = R2 (df /dt) - нормальное ускорение; R2 - радиус дефлектора, м.

Рис. 3. Схема движения частиц по поверхности дефлектора

Возможны два типа движения частицы по дефлектору. Первый соответствует безотрывному движению по его поверхности при небольших скоростях движения. Во втором случае возможен отрыв частицы от поверхности дефлектора из-за воздействия центробежной силы инерции. Для определения момента отрыва аналитически проинтегрирована указанная система уравнений. Величина силы трения при скольжении частицы:

F = f N

тр J ск ’

где N определяется из второго уравнения системы (5) в виде:

N = mgsin ф

R2

(6)

Тогда первое уравнение этой системы может рассматриваться как линейное неоднородное уравнение:

U - 2U = -g ( cos ф + /„sin ф), (7)

где через U обозначена величина

г2

U = ■

2 R

(8)

Решение дифференциального уравнения (7) может быть представлено суммой общего

решения однородной его части:

U„- 2U = 0 ^ U„ = Се2 (9)

и частного решения

Uчаст = Acos ф + Bsin ф, (10)

неопределенные коэффициенты которого находятся из уравнения (7).

После подстановки в общее решение

U = Цод^ ^^част (11)

начальных условий ф0 = у и U0 = VB2 / R2 получим выражение, определяющее в конечном виде величину (8), входящую в выражение для нормальной составляющей силы реакции (6).

Отрыв частицы от поверхности дефлектора произойдет при условии N = 0, т. е.:

gsin ф - U = 0.

(12)

Отсюда получается трансцендентное уравнение, определяющее угол отрыва частицы от дефлектора:

R2 f (ф - ф0 )

sin ф = — e v ' х

g

2 _ 2g

R

R

ЩЦ [sin Y (v~ -1)+3f-cos Y ]j

(f +1) [sin ф(f -1)]+%кс“ ф, (13)

VB

где ®0 =-£-.

Rr,

В результате интегрирования получены значения угла отрыва частицы от дефлектора и ее скорости в этот момент.

Условие отрыва частицы от трамплина имеет вид:

при G < или V ОgR (14)

R2

Если дуга BC симметричная, то при неотрывном движении частиц условие ее схода с криволинейной поверхности:

VC > 0. (15)

По закону сохранения энергии:

Е0 = Е + Е , (16)

0 с тр’ 4 '

№ 4(36), 2014

35

тг2

mVC

2

C

+ 1 Fw dl,

B

(17)

где l — длина окружности, м.

Для обеспечения условия (13) необходимо, чтобы в верхней точке D дуги ВС кинетическая энергия частицы превышала работу сил сопротивления движению на пути DC:

а mV2

mSfскR 1 C0s а da < —f-

o 2

(18)

Используя неравенство (18) и закон сохранения энергии, можно получить условия безотрывного движения частицы по дуге BC:

mVB > 2R2

(1 — cos f) mg +2mgfZK J cos a da

(19)

Если дугу BC ограничить в точке D, то в этом случае частицы падают вниз с начальной скоростью V = 0 и минимальной скоростью по оси х.

Условие безотрывочного и безостановочного движения частицы до точки D (находится из закона сохранения энергии):

mVB > 2R2

ф

(1 —cos ф) mg + mgfск J ras a da

о

.(20)

Для обеспечения движения частицы с отрывом от дефлектора и без встречи ее с последним необходимо выполнить условие:

m (^1п a)2 > mg(1 - cos a)R22,

V > l2g(l-c°s ф)I 2gR

B \ sin2ф V (1 + cos ф)

Как указывалось выше, участков свободного полета у частицы возможно несколько. Первый из них соответствует вылету частицы с криволинейного трамплина со скоростью VB под углом (в — у) к горизонтали (рис. 4). Поэтому проекции начальной скорости на этом участке на оси координат системы хОу имеют вид:

|V = vb cos (в — Y).

[К. = VBsin (в-т). (21)

Координаты точки начала полета опреде-

ляются через параметры полки и криволинейного трамплина:

[ хВ = Lcos в + rsin в — rsin (в — у), lyB = Zsin в — rcos(в — у) — rcos в. (22)

Если сопротивление движению при свободном полете отсутствует, то частица движется с ускорением свободного падения g, направленным параллельно оси Оу, поэтому движение вдоль оси Ох равномерное, а вдоль оси Оу — равноускоренное. В этом случае имеем:

Гх = XB + VJ,

[ У = yB + Vy t + 0,5 gt 2. (23)

Следовательно, движение частицы происходит по параболе

У = yB + tg(P — y)(x - xB) + [g(x - XB)2]/

Vв2cos2(в — y). (24)

Конечными условиями для участка свободного полета являются условия соударения с одним из отражательных элементов.

Положение /-го отражательного элемента характеризуется следующими геометрическими параметрами: его длиной l и шириной боковой грани b, м; координатами y z. вершины двухгранного угла 0 в месте его крепления, м; углом наклона у элемента к горизонтали (рис. 4).

Рис. 4. Геометрические параметры отражательного элемента

В каждый момент времени полета частицы выполняются проверки: нахождение

36

Известия Уральского государственного горного университета

ближайшего элемента; сравнение координат частицы с координатами этого отражательного элемента; проверка возможности взаимодействия частицы с боковой плоскостью сепаратора.

Встреча с элементом фиксируется, в случае если расстояние между частицей и плоскостью грани /-го отражательного элемента не превосходит задаваемой погрешности е:

■ 9 ■ ,( ч 9 ( ч ■ 0

xsin—sin у ± (y - yjcos—- (z - zjsin—cos у

< е (25)

(знаки ± соответствуют правой (+) и левой (-) грани отражательного элемента). При этом в силу ограниченности длины элемента и ширины его грани, а также в силу расположения отражательного элемента выпуклостью вверх, условие контакта с плоскостями грани должно соответствовать неравенствам:

z > Zj + x tg у,

x

<------

cos у

< l,

.0 - У/)2 + (z - z

x tg y)2 < b2.

(26)

Таким образом, процесс интегрирования системы дифференциальных уравнений, описывающих свободный полет частицы, прекращается для тех координат частицы (xk, yk, zk) и той грани отражательного элемента, которые удовлетворяют системе неравенств (25)-(26). Если при этом хотя бы одно из неравенств для всех x, y, z не выполняется, частица вылетает из зоны ударного разделения, не коснувшись ни одного из отражательных элементов.

Для определения скорости частицы после соударения с плоскостью отражательного элемента применяем теорему об изменении количества движения при ударе [2, 3]:

mU - mV = SN •n - STP т, (27)

где V - вектор скорости частицы в конце участка свободного полета в точке (xk, yk, U - скорость частицы после удара; п - единичный вектор нормали к плоскости грани, с которой произошел контакт; т - единичный вектор касательной, направленный вдоль линии пересечения плоскости удара с

плоскостью грани; S , S - величины ударных импульсов нормальной реакции N и силы трения скольжения, причем STP = f SN в соответствии с гипотезой Рауса.

Для удобства векторного представления коэффициент восстановления k при ударе определяется из зависимости [2], записанной в виде скалярного произведения:

U • n = -k V • n (28)

Чтобы определить величину ударного импульса нормальной реакции грани отражательного элемента, умножим векторное уравнение (27) скалярно на единичный вектор нормали n :

mU • n - mV • n = SN • n • n - STP т • n. (29)

Учитывая, что n • n = 1, т • n = 0, а также используя соотношение (28), получим

SN =-m (1 + kj)V• n. (30)

После подстановки этого выражения в уравнение (27) получим векторную зависимость, определяющую скорость частицы после удара:

U = V - (1 + j Vn • n • (n - fj •?), (31)

из которой нетрудно определить проекции скорости U на оси координат x, y, z, являющейся начальной скоростью на новом участке свободного полета:

Ux = Vx - (1 + kj) • (V,n, + Vyny + Vznz) • (n

‘Uy = Vy - (1 + kj) • (Vxnx + Vyny + Vznz ) • (ny

Uz = Vz -(1 + kj)• (Vxnx + Vyny + Vznz)• (n,

fTx)

■fTy)

fjT) (32)

Для удобства пользования формулами (32) определены проекции входящих в нее векторов через кинематические параметры частицы и геометрические параметры грани. Единичный вектор нормали n , перпендикулярный плоскостям граней отражательного элемента (25), имеет проекции:

n

. 0 .

sin —sin y, 2

n

y

±

0

cos—, 2

n =- sin—cosy.

z 2 r

0

(33)

№ 4(36), 2014

37

Единичный вектор касательной т можно найти из выражения двойного векторного произведения:

где b

т = n х b, (34)

вектор бинормали:

г V г b к—хn, (35)

V

V N V; + V2 + V2. (36)

Из уравнения (35) следует:

bx = V {Vynz - Vzny)

< by = V (Vznx - Vxnz)

bz = V(Vxny - Vynx) (37)

а из уравнения (34) находим проекции вектора т на выбранные оси координат:

т = n b - nb ,

x y z z y ’

< т = n b - n b ,

y z x

т = n b - n b .

z x y y X’

Отскок от грани отражательного элемента может и не произойти, если угол отражения близок к 90°. В этом случае из-за неровности поверхности частиц начнется их скольжение по плоскости грани. При этом скорость частицы будет убывать до тех пор, пока частица не покинет боковую поверхность грани. Если принять угловую погрешность расчетов равной п, то критерием начала скольжения по боковой поверхности грани может служить неравенство Un < Ux ■ tg п, то есть U n + U n + U n <

x x y y z z

< (Uxтx + Uyтy + Uzтz) ■ tg П- (39)

В этих условиях происходит не свободное, а стесненное движение частицы по плоскости грани:

• 0 . л, , 0

x Sin — Sin у л (y - y.) nos —

0

(z - z.)sin —cos у = 0.

(40)

Уравнение данного движения включает уравнение связи (40) и три дифференциальных уравнения, составленных на основе второго закона Ньютона для частицы:

V

m ■ а = mg + Nn - fj N—,

где m - масса частицы; a - её ускорение; N и fjN - величины нормальной реакции и силы трения скольжения, соответственно; g - ускорение свободного падения; n и V/v - направляющие вектора нормали к плоскости и касательной к траектории движения на плоскости.

Величину нормальной реакции N нетрудно найти из уравнения (41), спроектировав его на нормаль n :

m ■ а ■ П = mg ■ П + Nn ■ П - fj N

V ■ n

~V~ ’

откуда с учетом, что a ■ n = 0 и V ■ n лучим

N

mg ■ n = mg sin — cos у.

0

0, по-

(42)

Таким образом, уравнениями движения частицы по грани отражательного элемента являются следующие:

.0 0 . х x

x = g sin ^ cos у (sin ^ sin у - fj v ,

z = g - g sin 2 cos у (sin 2 cos у - fjV ,(43)

.----------- 0

где V = yjx2 + Vy2 + z2, Vy = Л t^-(z cos у -

- x sin у).

Интегрирование системы уравнений (43) следует вести до тех пор, пока частица не дойдет до нижней кромки грани, т. е. при

(44)

x = Ux; z = Uz

Начальными условиями для системы уравнений (43) являются значения проекций скорости

|z cos у - b - x sin у < s. (45)

и координат точки контакта с гранью xk, zk. Частица покидает поверхность грани в точке

38

Известия Уральского государственного горного университета

L, координаты которой xL, zL удовлетворяют неравенству (44), аyL определяется из уравнения (40)в виде:

Q

yL = Уi ±tg2[(z-zi)cosY-xsinY]. (46)

Таким образом, выражения (32) могут служить начальными значениями скорости на следующем участке свободного полета из точки соударения с координатами (х Ук, zk).

В случае скольжения частицы по грани свободный полет начинается в точке с координатами xL, y zL со скоростью, проекции которой определяются из системы уравнений (43). Далее вновь следует обратиться к решению системы дифференциальных уравнений свободного полета частицы с новыми начальными условиями. Чередование процессов интегрирования дифференциальных уравнений с соударениями о грани отражательных элементов происходит до тех пор, пока частица не покинет зону ударного разделения.

Указанный алгоритм реализуется на ПЭВМ со случайной выборкой начальных параметров движения (точки начала движения по наклонной плоскости и ее начальной скорости), а также случайным выбором параметров самой частицы (ее размеров и физических параметров трения и упругости). Блок-схема моделирования работы разделительного аппарата приведена на рис. 5. Результаты численного моделирования позволяют оценить эффективность разделения при выбранных конструктивных параметрах аппарата, а также выбрать их наиболее рациональные значения.

Организация вычислений при моделировании разделения частиц по различию в коэффициентах трения и упругих свойств заключается в следующем:

1. Моделируется разгон частиц по загрузочному лотку с учетом сил трения. Начальное положение частицы по ширине загрузочного лотка определяется генератором случайных чисел. Параметры частиц могут задаваться двумя способами: а) перебором по порядку дискретных значений параметров: для каждого сочетания параметров (типа ча-

стиц) реализуется несколько траекторий; б) генерированием случайной последовательности типов частиц в соответствии с заданным распределением частиц по вещественному и фракционному составу.

2. Моделируется траектория свободного полета частиц до встречи с отражательным элементом. Элементы траектории рассчитываются дискретно с автоматическим выбором шага по времени.

3. Определяются момент и координаты точки встречи частицы с отражающим элементом, боковой стенкой или дном.

4. Определяются условия перехода на скольжение. В случае перехода дальше рассчитывается траектория скольжения. В противном случае определяются траектории после отскока.

5. После взаимодействия с отражательным элементом повторяются вычисления по п. 2 до следующей встречи с элементом. Вычисления прекращаются при полной потере частицей начальной высоты. Блок-схема алгоритма моделирования приведена на рис. 5.

6. В процессе моделирования накапливается статистика о координатах частицы в момент падения на дно.

7. Для каждого типа частиц рассчитывается средняя длина траектории. Различие координат позволяет судить о разделении частиц по размерам и другим признакам. Результат выдается в виде таблицы и графика.

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

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

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

№ 4(36), 2014

39

ние угла наклона загрузочного желоба. Так Такая же закономерность имеет место при скорости подачи исходного продукта V0 при V0 = 0,5; 1,5; 2,0 м/с. С учетом ограничен-= 1,0 м/с увеличение угла наклона с 40 до 55° ного расстояния падения породных частиц снижает показатель разделения в 4 раза. лучшие условия разделения обеспечиваются

Рис. 5. Алгоритм моделирования процесса движения частицы в разделительном аппарате

при углах наклона загрузочного желоба в интервале 40-45° при V0 = 0,5...1,0 м/с. При постоянном угле наклона загрузочного желоба изменение скорости подачи исходного материала в интервале от 0,5 до 2,0 м/с существенно не влияет на разделение продуктов по трению. С увеличением длины желоба показатель разделения повышается в 1,5-2 раза, но при этом увеличивается дальность падения породы от точки отсчета. В этом случае габариты аппарата могут увеличиваться как по высоте, так и по длине.

Расчеты показывают, что для обеспече-

ния контакта всех частиц потока материала с поверхностью рабочих элементов необходимо их расположить, как минимум, в два яруса, причем ширина щели между соседними элементами должна быть меньше или равна ширине элемента.

Траектория движения частицы после отражения от поверхности элемента определяется параметрами движения частицы в момент взаимодействия с поверхностью разделительного элемента, параметрами отражающей плоскости, а также параметрами частиц, которые характеризуются их упругими

40

Известия Уральского государственного горного университета

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

1. Породные частицы крупностью более 2 мм в момент первого удара имеют скорости, равные (4,5-5,0) м/с. После удара скорость уменьшается до 3,4-3,8 м/с. Направление движения частиц в плане изменяется в среднем на угол 31° от оси у (продольной оси отражающих элементов). На породные частицы менее 2 мм значительное влияние оказывает сопротивление воздуха, поэтому скорость их падения в момент первого удара меньше, чем у крупных частиц. Так, у частиц размером 1 мм скорость падения в момент первого удара равна 4,0-4,1 м/с. Соответственно снижается скорость отражения.

2. Нераспушенные агрегаты асбеста («пешка») в момент первого удара имеют угол падения несколько больший, чем у породы, равный 63-64° и скорость 4,0-4,5 м/с. После удара скорость уменьшается до 2,5-3) м/с. Изменение направления движения от оси у составляет 31-32°.

Породные частицы имеют большую упругость, поэтому при отражении получают большую скорость и движутся по более пологой траектории чем «пешка», которые движутся после удара по более крутой траектории.

3. Распушенные волокна после первого удара о поверхность элементов теряют скорость движения с 1,7—1,8 м/с до 0,4-0,5 м/с.

4. Анализ траектории движения частиц

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

5. Для всех значений двугранного угла с увеличением угла наклона элементов показатель разделения повышается и достигает максимальной величины при g = 25°. С дальнейшим увеличением угла наклона элементов показатель разделения падает.

Максимальные значения показателя разделения получены при двугранных углах 0 в интервалах 90-110°.

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

В интервале ширины элементов от 0,05 до 0,15 м показатель разделения практически не изменяется и равен 1,38-1,4 м. С дальнейшим увеличением ширины элементов показатель разделения снижается до 1,1 м для b = 0,3 м. Наиболее рациональными можно считать ширину элементов b = 0,14 - 0,15 м, что при двугранном угле 0 = 90° соответствует уголку № 100 фасонного проката.

7. Скорость подачи частиц в загрузочный желоб должна составлять 1,0.. .1,5 м/с.

8. Минимальная рабочая длина отражательных элементов составляет 1,14 м с учетом рационального угла наклона g = 25.

На основании проведенных исследований методом имитационного моделирования была предложена конструкция аппарата с учетом рассмотренных условий [2].

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Вебер Г. Э., Ляпцев С. А. Дополнительные главы механики для горных инженеров. Екатеринбург: УрГУ, 1998.

2. Ляпцев С. А., Потапов В. Я. Математическое описание поведения рудных частиц в воздушном потоке разделительных аппаратов // Современные проблемы науки и образования. 2012. № 1. URL: www.science-educa-tion.ru/101-5493.

3. Ляпцев С. А., Потапов В. Я., Афанасьев А. И. Аппараты для разделения горных пород по упругофрикционным свойствам. LAP LAMBERT Academic Publishing, Deutschland, Saarbrucken 2014. 90 c.

4. Никитин Н. Н. Курс теоретической механики. СПб., 2010. 720 с.

№ 4(36), 2014

41

5. Понтрягин Л. С. Обыкновенные дифференциальные уравнения. М.: Наука, 1982. 331 с.

6. Методика определения упругих и фрикционных характеристик сыпучих материалов / В. Я. Потапов [и др.] // Изв. вузов. Горный журнал. 1998. № 5-6. С. 103-108.

7. Сепаратор для разделения материалов по трению и упругости / В. Я. Потапов [и др.]// Патент РФ № 111780, 2011. Бюл. № 36126.

Поступила в редакцию 20 октября 2014 г.

Ляпцев Сергей Андреевич - доктор технических наук, заведующий кафедрой технической механики, профессор. 620144, г. Екатеринбург,ул.Куйбышева,30, Уральский государственный горныйуниверситет. E-mail: [email protected]

Потапов Валентин Яковлевич - доктор технических наук, профессор кафедры горной механики, доцент. 620144, г. Екатеринбург, ул. Куйбышева, 30, Уральский государственный горный университет. E-mail: [email protected]

Потапов Владимир Валентинович - кандидат технических наук, доцент кафедры горной механики. 620144, г. Екатеринбург, ул. Куйбышева, 30, Уральский государственный горный университет. Семериков Леонтий Андреевич - аспирант кафедры горной механики. 620144, г. Екатеринбург, ул. Куйбышева, 30, Уральский государственный горный университет.

Беридзе Елизавета Тимуровна - инженер кафедры горной механики. 620144, г. Екатеринбург, ул. Куйбышева, 30, Уральский государственный горный университет.

42

Известия Уральского государственного горного университета