COMPUTER-AIDED DESIGN SYSTEM APPLICATION FOR CALCULATION PARAMETERS OF STRENGTH OF BALL MILL CLASSIFICATION DEVICES
Abstract and keywords
Abstract (English):
Constructive improvement of the classification devices allows increasing capacity of ball mills, to reduce the specific energy consumption for the grinding process of the material. Developing methods of calculation that allows you to determine the magnitude and direction of applying loads to the design of the device is actual. The article describes the mathematical model of the process of the spatial motion of grinding media, their interaction with shell liners, and each other in the cone-shaped and cylindrical chambers of the mill. The method of calculating loads on classifying devices from the action of grinding media is considered based on a mathematical model. The results of calculation loads, strains and deformations acting on the cylindrical classifier of the industrial using mill in its design obtained using computer-aided design systems are presented

Keywords:
ball mill, classification device, grinding media, mathematical model of the process, methods, load, deformation
Text
Publication text (PDF): Read Download

Введение. Шаровые мельницы распространены в различных отраслях промышленности для грубого и тонкого помола различных материалов. Эти мельницы имеют относительную простоту конструкции и эксплуатации. Их универсальность позволяет использование в условиях сухого и мокрого способа измельчения как однокомпонентных, так и многокомпонентных шихт и шламов, обеспечивая достаточно высокую однородность распределения частиц в измельченном продукте по их крупности и составу. Нерациональность процессов движения шароматериальной загрузки, выделения из ее среды кондиционных частиц материала приводит к невысокой эффективности процесса помола, следствием чего является снижение производительности и повышение удельного расхода электроэнергии. Это ограничивает преимущества шаровых мельниц перед аналогичными помольными агрегатами [1–5]. Эволюционное совершенствование шаровых мельниц осуществляется на протяжении всей истории их существования. Одним из направлений их модернизации является изменение конструкций внутримельничных устройств, позволяющих осуществлять процесс выделения из шароматериальной загрузки частиц материала, достигших определенной крупности, и подавать на последующее доизмельчение или направлять в готовый продукт [6]. Это значительно снижает нерациональное переизмельчение мелких частиц и их демпфирующее воздействие на мелющие тела при разрушении крупных, что повышает эффективность процесса помола.

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

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

Так, в работе [7] предлагается методика расчета сил, действующих на элементы цилиндрического классифицирующего устройства, основывающаяся на графоаналитическом методе [8].

,            (1)

где L - длина k-го сектора классифицирующего устройства, м; i – номер сектора контура загрузки (i = 1,2,3), γшм - объемный вес шароматериальной загрузки, Н/м3; g - ускорение свободного падения, м/с2; ρ и ρ – радиусы, ограничивающие контур мелющей загрузки, м; ∆νк – угловой сектор, рад.

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

Методика. Заслуживает особого интереса дискретно-событийный подход к процессу движения мелющих тел [9–13]. В рамках этого подхода разработана математическая модель, адекватно описывающая протекающие в цилиндрической или конусообразной камере шаровой мельницы процессы пространственного движения мелющих тел, взаимодействия с бронефутеровкой и друг другом. В основу математической модели положено движение отдельно взятого сферического мелющего тела, которое последовательно взаимодействует с другими мелющими телами сферической формы, а также бронефутеровкой камеры и ограничивающих ее перегородок или днищ. Взаимодействия сферических мелющих тел с бронефутеровкой и друг с другом рассматриваются как мгновенные и вызывают изменения скоростей, которые изменяются также мгновенно. Мелющие тела на интервалах времени между соударениями движутся под действием силы тяжести по параболическим траекториям.

Время движения мелющего тела до бронефутерованной конусообразной поверхности камеры определяется из уравнения:

 (2)

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

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

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

Обозначим индексы параметров следующими цифрами: 0 и 2 – до и после взаимодействия; 1 – в момент, когда сила трения обращается в ноль.

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

(3)

,(4)

(5)

где

  .

 

 

 

Рис. 1. Бронефутерованная камера в абсолютной

системе координат

 

Рис. 2. Мелющее тело в локальной системе

координат

 

Нагрузка от воздействия мелющих тел с силой  на фрагмент поверхности бронефутеровки камеры площадью  вычисляется как  = . Касательная  и нормальная составляющие силы  рассчитываются как:

         (6)

Здесь – период времени, за которое устанавливается нагрузка; n – количество взаимодействий с мелющими телами  за период времени ; и – касательный и нормальный  импульсы -го взаимодействия.

Бронефутерованные поверхности цилиндрической или конусообразной (рис. 3) камер разделяются на участков в направлении оси вращения  и на секторов в поперечных сечениях. Определим номер участка на бронефутерованной поверхности конусообразной камеры, взаимодействующего с мелющим телом в точке, имеющей координаты . Конусообразная камера характеризуется: углом  наклона образующей, длиной  диаметром  у большего основания.

Длина  участка (рис. 3), при их количестве , равна:

                               (7)

Рис. 3. Схема к разделению конусообразной камеры на участки

 

Его номер (в направлении оси ):

                            (8)

Угол между ограничивающими сектор радиусами:

                           (9)

Номер этого сектора:

                       (10)

где – угол между радиус-вектором точки взаимодействия и осью .

                       (11)

Развёртка выделяемого фрагмента поверхности бронефутеровки камеры представляет форму трапеции. Её высота вычисляется как , а размеры оснований – нижнего и верхнего :

           (12)

    (13)

Площадь выделяемого фрагмента поверхности в форме трапеции:

(14)

Основная часть. Для компьютерной реализации процессов движения мелющих тел в камере цилиндрического классифицирующего устройства, их взаимодействия с бронефутеровкой и друг с другом формируются начальные положения мелющей загрузки [14]. Затем для мелющих тел осуществляется расчёт процессов их движения, взаимодействия с бронефутеровкой и друг с другом. На бронефутеровке устройства устанавливается номер фрагмента и координаты точки взаимодействия с мелющим устройством. Импульсы взаимодействия мелющих тел с бронефутеровкой суммируются по соответствующим фрагментам за назначенный промежуток времени и рассчитываются действующие на фрагменты нагрузки. 

Для шаровой мельницы D×L=2×10,5 м осуществлен расчёт напряжённо-деформированного состояния цилиндрического классифицирующего устройства. Была принята рабочая относительная частота вращения барабана мельницы ψ = 0,76 ψкр. Устройство длиной L = 0,9 м, с внутренним диаметром Dцку = 1,8м выполнено из стальных стержней диаметром dc = 0,06 м, расположенных с промежутками и образованием щелевых отверстий шириной aщо = 0,003 м и длиной lо = 0,2 м. Материалом стержней является сталь 09Г2С ГОСТ 2590-206 ([σт ] = 275 МПа). Устройство загружено мелющими телами с коэффициентом φ = 0,3. Так как оно длительное время эксплуатируется в условиях высокого абразивного износа, то рассматриваем состояние при 50 % износе стержней по их диаметру. Воспользовавшись, в соответствии с [14], компьютерной реализацией процессов движения мелющих тел в камере устройства, их взаимодействия с бронефутеровкой и друг с другом рассчитали нагрузки на фрагменты конструкции устройства. Полученные результаты, приведенные в таблице 1, применили для расчёта напряжённо-деформированного состояния цилиндрического классифицирующего устройства в CAD/CAM/CAE системе Unigraphicx NX [14]. Результаты расчета напряжений и деформаций в конструкции устройства приведены на  рисунках 4-6.

 

Таблица 1

Нагрузки на цилиндрическое классифицирующее устройство мельницы D×L = 2×10,5 м

Нагрузка на фрагменты устройства Р, Н

 

Угловой номер фрагмента

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

 

Угол, град.

23

28

32

37

42

46

51

55

60

65

69

74

78

83

88

92

97

102

106

111

115

120

125

129

134

138

143

148

152

157

162

166

171

175

180

185

189

№ фрагмента в направлении продольной оси

1

0

68

27

253

470

869

1194

1678

2060

2180

2956

3150

3238

3450

3361

3071

3072

3088

2645

3084

2948

3013

2892

2618

2303

2149

2114

1714

1457

1096

996

767

523

337

231

0

0

2

0

35

167

483

712

1638

1865

2503

3399

4047

4495

5083

5034

5187

4574

5562

5957

5549

5048

4762

4648

4587

4336

4202

3318

3311

2881

2761

2150

1814

1465

1204

708

455

0

0

290

3

15

126

162

508

1063

1386

1867

2683

3476

3076

3584

4111

4028

4368

4373

4742

4180

4090

4488

4531

4632

4533

4159

3599

2865

2902

2641

2400

1916

1525

1002

798

632

314

1095

0

0

4

15

39

174

641

1178

1747

2388

2601

2980

3572

4101

4739

4569

4856

4950

4511

4760

4426

4554

3770

3747

3411

3473

3400

3167

2694

2475

2131

1889

1471

1236

937

712

465

0

0

0

5

0

53

195

558

894

1422

1853

2521

3687

3787

3804

4177

4916

4890

4167

4256

4976

4115

4476

4692

4138

3692

3723

3601

3118

2423

2269

1864

1735

1499

1109

660

433

279

225

229

0

6

0

32

175

414

700

1505

1799

2697

3031

3265

4655

4738

5254

5270

5266

5777

5186

5585

5536

4841

4729

4753

5112

4209

3950

3882

3270

3078

2520

1887

1483

1134

755

469

255

428

0

7

0

21

109

485

457

925

1357

1549

2083

2738

3249

2904

3120

2863

3171

3343

3439

3314

3047

3075

3233

2776

2679

2270

2198

1993

1824

1690

1368

1214

846

708

524

355

272

0

366

 

Максимальные значения напряжений в цилиндрических стержнях – σ = 255,3 Мпа (меньше допускаемого [σт ] на 7 %); перемещений –
8,3
10-4 м характеризуют рабочее состояние устройства после  длительной эксплуатации.

 

 

 

а)                                                   б)

Рис. 4. Распределение напряжений и деформаций в конструкции цилиндрического

классифицирующего устройства  мельницы D×L=2×10,5 м: а) напряжения, б) деформации

 

а)                                                   б)

Рис. 5. Распределение напряжений и деформаций в стержне цилиндрического

классифицирующего устройства  мельницы D×L=2×10,5 м: а) напряжения, б) деформации

 

 

 

   

а)                                      б)

Рис. 6. Распределение напряжений и деформаций в дистанционной вставке цилиндрического

классифицирующего устройства  мельницы D×L=2×10,5 м: а) напряжения, б) деформации

 

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

Выводы. Получена математическая модель, описывающая параметры процессов движения мелющего тела в конусообразной и цилиндрической камерах шаровой мельницы, взаимодействия с их бронефутеровкой и друг другом. Использование математической модели позволяет осуществлять компьютерную реализацию рассматриваемых процессов в камерах мельницы для всей мелющей загрузки в совокупности и устанавливать характеризующие их параметры. Разработана методика расчета нагрузок от воздействия мелющих тел на фрагменты поверхности бронефутеровки камер. С использованием разработанного программного обеспечения выполнен расчет нагрузок, действующих на цилиндрическое классифицирующее устройство промышленной шаровой мельницы. В программной среде CAD/CAM/CAE системы Unigraphicx NX рассчитаны напряжения и деформации в конструкции устройства.

*Работа выполнена в рамках Программы развития опорного университета на базе БГТУ им. В.Г. Шухова.

References

1. Romanovich A.A. Energosberegayuschie agregaty dlya izmel'cheniya materialov ce-mentnogo proizvodstva s anizotropnoy strukturoy: Avtoref. dis. d-ra. tehn. nauk: Belgorod. 2014. 43 s

2. Naumenko Yu.V. Osnovi teorii re-zhimiv roboti barabannih mliniv. Rivne.: SPD Zelent O.I., 2009. 282 s.

3. Yorgensen. S.V. Pomol cementa. Sravnitel'nyy analiz vertikal'noy valko-voy i sharovoy mel'nic // Cement i ego pri-menenie. 2006. № 6. S. 39-44.

4. Sevost'yanov V.S. Energosberegayuschie pomol'nye agregaty. Belgorod. 2006. 435 s.

5. Harder Y. Razvitie odnoetapnyh pro-cessov izmel'cheniya v cementnoy promysh-lennosti // Cement. Izvest'. Gips. 2006. № 1. S. 24-38.

6. Pat. 2236298 Rossiyskaya federaciya, MPK VO2S 17/06. Trubnaya mel'nica / S.I. Hanin , V.V. Kaydash, A.V. Chalov, D.N. Solo-dovnikov, O.S. Hanina; zayavitel' i patento-obladatel' Belgor. gos. tehnolog. akad. stro-it. materialov. - №2003113249/22; zayavl. 05.05.03; opubl. 20.09.04, Byul. №26. 6 s.

7. Solodovnikov D.N. Cilindricheskoe vnutrimel'nichnoe klassificiruyuschee ustroystvo trubnoy mel'nicy: dis... kand. tehn. nauk. Belgorod, 2009. 162 s.

8. Kryukov D.K. Futerovki sharovyh mel'nic. M.: Mashinostroenie, 1965. 184s.

9. Brilliantov N.V., Spahn F., Hertzsch Jan-Martin, Thorsten Poschel Model for collisions in granular gases // Phisical review E. 1996. Vol. 53. № 5. P. 5382-5392.

10. R. Reichardt and W. Wiechert / Event driven simulation of a high energy ball mill. In Proceedings ASIM 2003, page 249. ASIM 2003.

11. Rosa Ramirez, Thorsten Pöschel, Nikolai V. Brilliantov, Thomas Schwager // Phys-ical review e. october 1999. Vol. 60. № 4. P.4464-4472.

12. Schwager T., Poschel T. Coeffi-cient of restitution for viscoelastic spheres: The effect of delayed recovery // Physical review E V.78. 2008. № 5. November. P. 1304-1316.

13. Thorsten Pöschel Thomas Schwager // Computational granular dunamics. Models and algoritms. Springer - Verlang Berlin Heidelberg 2005. P. 322.

14. Svidetel'stvo o gosudarstven-noy registracii programmy dlya EVM №2012615855. Sistema rascheta i vizualiza-cii parametrov dvizheniya melyuschih tel v korpuse sharovoy barabannoy mel'nicy / Starchenko D.N., Hanin S. I., Truhachev S.S.; zayavitel' i patentoobladatel' BGTU im. V.G. Shuhova; №2012613407, zayavleno 27.04.2012; zaregistrirovano. 27.06.2012.

15. El'cov M.Yu., Kozlov A.A., Se-doykin A.V., Shirokova L.Yu. Proektirovanie v NX pod upravleniem Teamcenter. 2010.783 s.


Login or Create
* Forgot password?