© Ю.В. Дмитрак, Т.Д. Зиновьева, Н.Н. Сычёв, 2007
УДК 6221.926.55
Ю.В. Дмитрак, Т.А. Зиновьева, Н.Н. Сычёв
РАЗРАБОТКА РАСЧЕТНОЙ ДИНАМИЧЕСКОЙ МОДЕЛИ РАСПРЕДЕЛЕНИЯ НАГРУЗОК В ПОМОЛЬНОЙ КАМЕРЕ ВИБРАЦИОННОЙ МЕЛЬНИЦЫ
Семинар № 20
Опыт эксплуатации вибромельниц показывает, что наиболее слабыми местами с точки зрения прочности и долговечности конструкции являются подшипниковые узлы, несущая рама и помольная камера. В настоящей работе мы остановились на рассмотрении вопроса о прочности помольной камеры, т. к. именно на неё первоначально передаётся действие ударных сил со стороны мелющей загрузки. Кроме того, в МГГУ накоплен большой опыт в исследовании динамики мелющей загрузки и созданы предпосылки для начала исследований по определению влияния динамических параметров загрузки на прочностные характеристики помольной камеры. В частности группой исследователей было разработано устройство по определению динамических характеристик мелющей загрузки. Впервые стало возможным получить картину распределения сил, действующих на стенку помольной камеры. Эти данные легли в основу настоящей работы.
Характер нагружения камеры является стохастическим, при котором различные части ударных шаров воздействуют в разные времена неравномерно на стенки цилиндрической оболочки. Тем не менее, из-за циклического изменения нагрузки, обусловленной вращением возбуждающего
двигателя пространственную неравномерность нагружения можно считать периодической по времени. Для статического анализа конструкции мельницы можно рассмотреть квази-статические нагрузки, действующие в различные времена.
Г еометрическая модель конструкции вибромельницы показана на рис. 1. Основными ее силовыми элементами являются помольная камера, несущая рама, приводной вал с подшипниковыми узлами и дебалансами и пружины, обеспечивающие колебательные движения мельницы при работе двигателя. Рассмотрим нагружение конструкции в начальный момент времени. Будем считать, что как радиальная, так и касательная составляющие нагрузки имеют максимальные значения в нижней части оболочки, а минимальная амплитуда нагружения соответствует верхней части оболочки. Пусть амплитуда радиальной составляющей в нижней части оболочки равняется ?1е, а в верхней части - Р2е.
На шар действуют:
= т ■ апе = т ■ю1 ■ Н - (1)
нормальная переносная составляющая силы инерции, где апе - нормальное переносное ускорение, сое - угловая скорость вращения камеры;
Помольная камера
Рама
Пружины
Подшипниковый узел
Рис. 1. Геометрическая модель конструкции вибрационной мельницы
Рт = т ■ аТ
(2)
тангенциальная переносная составляющая силы инерции, при установившемся движении РТ = 0 ;
Рг" = т • а" = 0 - (3)
нормальная относительная составляющая силы инерции, при установившемся движении Р" = 0 .
На рис. 2 показана эпюра распределения нагрузок, действующих на стенки помольной камеры, для этого положения на примере переносной составляющей силы инерции. Нагрузки состоят из радиальной составляющей, которая пропорциональна массе ударных шаров, квадрату возбуждающей частоты колебаний и радиусу
цилиндрической оболочки и силы Кориолиса, направленной по касательной к оболочке. Опыт исследования динамики мелющей загрузки показывает, что нормальная переносная составляющая силы инерции и сила Кориолиса являются величинами большего порядка по сравнению с остальными силами, действующими на мелющую загрузку. Поэтому именно эти силы в дальнейшем учитываются при проведении статических и динамических расчетов конструкции вибромельницы. Следует отметить, что максимальное отношение значений радиальных составляющих сил, действующих в верхней и нижней частях помольной камеры равно:
F
F- = 2,8 F,
(4)
На рис. 3 показаны компоненты радиальных сил по направлениям X и У.
Характер нагружения контейнера является стохастическим, при котором различные части ударных шаров воздействуют в разные времена неравномерно на стенки цилиндрической оболочки. Тем не менее, из-за циклического изменения нагрузки, обусловленной вращением возбуждающего двигателя пространственную неравномерность нагружения можно считать периодической по времени. Для статического анализа конструкции мельницы рассматриваются ква-зистатические нагрузки.
При задании нагрузок в расчетной программе удобно пользоваться заданием компонент нагружения Рх и Ру. Компоненты нагружения могут быть легко получены из следующих соотношений:
(5)
Fy = А(р) sin(p)
Fx = А(р) cos(p) ’
где А(р)- амплитуда силы, действующей на оболочку.
Исходя из того, что значение амплитуды силы в начальный момент
/ \
времени в нижней точке ( р = — )
равно Fie, а в верхней точке (р = П) -Fie, используя линейную интерполяцию
Рис. 3. Базовые компоненты сил, действующих со стороны загрузки на помольную камеру
с помощью интерполяционного полинома Лагранжа можно записать выражения для нахождения амплитуд ,п 3л.
нагрузок на дуге ре :
А(р) = П ^ Fie (Р-П) - F,е (р-3П) j . (6)
Тогда компоненты нагрузки в рассматриваемой части оболочки могут быть записаны следующим образом:
Fv (Р) = П (Fie (Р -1) - F,e (Р - 3П) j sin(p)
(7)
Fy (р) = П ^ Fie (Р -1) + F,е (Р +1) ) sin(p)
(10)
Fx (Р) = П ^ (Р -П>') + F2e (Р + П
(ii)
Поскольку в системе расчета отсчет углов может начинаться со значения углар = п , то удобным будет воспользоваться вышеприведенными формулами вычисления амплитуды и компонент нагрузки на отрезе
Г П1
Ре [-П - —] :
Fx (Р) =П f Fie (Р-f) - F,e (р-П cos^) (8) А(р) = - f Fie (р + п) - Fie 2 F,e (р+П)
\ J п V 2 2
Аналогично можно записать выражения для вычисления амплитуды и компонент радиальной нагрузки на РУЫ = ~|Р1в(Р + п)
, п п, дуге ре [- — ^ :
(i2)
Fie 2 Fe р + п) I sin(р)
(i3)
А(Р) = П f-Fie (р-|) + F2e (р + |)| (9)
Fx (р)=П f Fie (р+п) -Fie 2 F"e (р+п i
(i4)
а)
б)
в)
1.2Е+06 -
^^2.0Е+05 -
2 -1 ) 1 2 3 4 Є
-8.0Е+05 --1.0Е+06 -
-Дтрі
РУ
РХ
Рис. 4. Зависимости составляющих радиальной нагрузки от угла ф а) - амплитуда и компоненты радиальной составляющей нагрузки; б) - амплитуда и компоненты сил Кориоли-са; в) - суммарные амплитуда и компоненты сил
Для статического анализа при действии нагрузки со смешенными на определенный угол экстремальными значениями нагрузок достаточно
только пересчитать компоненты нагружения в соответствуюшей цилиндрической системе координат, повернутой на этот угол вокруг оси 7.
На рис. 4 показаны зависимости составляюших радиальной нагрузки от угла ф. Амплитуды радиальной составляющей нагрузки в нижней и верхней частях цилиндрической поверхности соответственно равнялись
/1е = 652050 Н и /2е = 360450 Н. Амплитуды сил Кориолиса в нижней и верхней частях цилиндрической поверхности соответственно равнялись /1к = 434700 Н и /2к = 160200 Н.
В результате полученной расчётной схемы распределения нагрузок в помольной камере вибрационной мельницы стало возможным с помощью метода конечных элементов произвести прочностной расчёт отдельных узлов вибрационной мельницы. гггга
— Коротко об авторах-----------------------------------------------------
Дмитрак Ю.В., Зиновьева Т.А. - Московский государственный горный университет. Сычёв Н.Н. - ООО «Даки».
ДИССЕРТАЦИИ
ТЕКУЩАЯ ИНФОРМАЦИЯ О ЗАЩИТАХ ДИССЕРТАЦИИ ПО ГОРНОМУ ДЕЛУ И СМЕЖНЫМ ВОПРОСАМ
Автор Название работы Специальность Ученая степень
ТЮМЕНСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ
БЕКМАН Александр Дмитриевич Некоторые подходы к построению 30 геологических моделей с использованием разнородной априорной информации 05.13.18 к.ф.-мат.н.
ЮЖНО-РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ
ИСАКОВ Владимир Семенович Принципы построения и синтез функциональных механизмов строительных, подъемно-транспортных и горных машин с напряженными замкнутыми кинематическими контурами 05.05.04 05.05.06 д.т.н.