Full Text

Введение На сегодняшний день повышение эффективности процесса измельчения материалов является одним из актуальных вопросов во многих сферах производства. Несмотря на широкий ассортимент существующих конструкций измельчительных машин, принципиально отличающихся друг от друга спецификой измельчения, потребность в совершенствовании последних остается на прежнем уровне [1]. Прежде всего это связано с современной тенденцией к миниатюризации, идея которой заключается в нанотехнологическом открытии и исследовании совершенно новых свойств наночастиц веществ. В связи с этим возросли и требования к проектируемым измельчительным машинам. С одной стороны, конструируемая технологическая машина должна обеспечивать необходимый минимальный размер так называемой «частицы-продукта» (т.е. конечный размер частицы после помола), а с другой - обладать высокой энергоэффективностью и производительностью. Данный подход поставил перед проектировщиками измельчительной техники ряд новых перспективных задач, связанных не только с обеспечением процесса измельчения как такового, но и с созданием инновационных машин, позволяющих адаптировать режим дробления к конкретному измельчаемому материалу с учетом особенностей его физико-механических свойств. Поставленная задача требует разработки соответствующего математического аппарата (метода расчета) для определения фактического значения прилагаемой нагрузки, необходимой для разрушения материала, от которой, в свою очередь, зависит энергоемкость процесса измельчения. При этом следует учесть ряд индивидуальных параметров обрабатываемого материала, связанных с исходным размером частиц, их прочностными свойствами, условиями разрушения и т.д. Одним из прогрессивных направлений совершенствования процесса измельчения является разрушение твердого вещества посредством многократного приложения силы переменной (периодической) величины [2-8]. Известно, что при циклической нагрузке разрушение материала возможно и в основном реализуется при напряжениях ниже предела текучести [9] в области упругих деформаций [10]. Можно констатировать, что эффективность применения данного принципа усиливается при использовании в процессе дробления мелющих тел, имеющих, как правило, сферическую форму. К таким машинам-измельчителям относятся шаровые, вибрационные, бисерные, планетарные и другие мельницы [1], в которых помол материала происходит при его интенсивном контактном взаимодействии с мелющими телами-шарами, выполняющими функцию передатчика потенциальной энергии разрушения частицам измельчаемого вещества. Количество энергии, требуемой для измельчения камня до определенного размера, зависит от размера, формы, взаимного расположения частиц, прочности, хрупкости, однородности, изотропности исходного материала и т.д. Поэтому установить аналитическое соотношение между энергоемкостью и механическими характеристиками измельчаемого вещества можно лишь приближенно. Следует также отметить, что движение мелющих шаров осуществляется в довольно сложном (хаотичном) режиме. Эти обстоятельства создают определенные трудности для разработки соответствующей механико-математической модели, позволяющей достаточно точно и адекватно оценивать напряженно-деформированное состояние и прочность твердого хрупкого вещества при его разрушении в процессе измельчения, зависящем от многих факторов случайного характера, трудно поддающихся расчету [10-13]. Сложность процесса помола, исключающая возможность создания единой универсальной теории, привела к появлению общепринятых гипотетических предпосылок о замене действительных кусков материала самой разнообразной формы идеализированными расчетными схемами кубического или шарообразного вида [10]. Как показал анализ экспериментально-теоретических исследований в этой проблемной области, наблюдается отсутствие конкретных физико-математических зависимостей, реалистично аппроксимирующих процесс измельчения [1, 10]. Представленная статья посвящена математическому описанию процесса помола невесомой сферической модели частицы измельчаемого хрупкого материала радиусом [10] путем ее соударения с шарообразными абсолютно твердыми телами большого радиуса и массой (рис. 1). На механическую систему, изображенную на рис. 1, действуют в противоположных направлениях некоторые внешние уравновешенные массы , имеющие заданную скорость в начальный момент времени , предшествующий присоединению к . Рис. 1. Общая расчетно-конструктивная модель взаимодействия разрушаемой частицы 1 с шарами-измельчителями 2 Fig. 1. General computational and constructive model of interaction of the destroyed particle 1 with the grinding balls 2 Исследование ударных процессов относится к числу наиболее востребованных практикой фундаментально-прикладных разделов сопротивления материалов и теоретической механики, посвященных оценке поведения различных конструкций в условиях интенсивного силового воздействия импульсного характера. Такие динамические нагрузки имеют место при эксплуатации многих современных сооружений, машин, механизмов и приборов [14-16]. Сформулированный в названии статьи вопрос является предельным частным случаем [14] конструктивно-нелинейной [17] задачи механики сплошной твердой среды об ударно-силовом контактировании двух шаров [15, 18, 19], движущихся вдоль оси, соединяющей их центры, со скоростями и после их точечного соприкосновения. В ее основе - классическая аналитическая зависимость Герца-Штаермана [18, 20] о реактивной статической нагрузке , возникающей между сферическими телами, когда (рис. 2): , (1) где , - радиусы шаров; u - сближение взаимодействующих тел в результате их локальной (контактной) деформации; , , , - соответственно коэффициенты Пуассона и модули упругости контактирующих материалов. Рис. 2. Классическая схема соударения шарообразных элементов по Герцу [18-21] Fig. 2. The classical scheme of the spherical elements collision according to Hertz [18-21] В известной фундаментальной теоретической публикации [18] Г. Герц высказал основополагающее предположение, подтвержденное экспериментальной проверкой [14, 15, 22], что «комбинируя статическое сжатие в частях тел, лежащих непосредственно у места соприкосновения, с общими уравнениями движения для остальных частей тел, мы, вероятно, могли бы получить закон для соударения тел любой формы». В то же время следует учитывать в настоящее время тот очевидный факт, что точная теория удара, основанная на местных деформациях и колебательно-волновом распространении напряжений в упругом теле, оказывается пока еще очень сложной задачей [15, 21]. Поэтому ниже излагается прикладная математическая модель ударно-силового импульса, возникающего в контактирующих шарообразных элементах (см. рис. 1). Принимая во внимание допущение о недеформируемости мелющих шаров, а также симметрию и буквенные обозначения схемы рис. 1, представляем формулу (1) в более компактном упрощенном виде: , (2) после выполнения замен , , , , , (3) где , - упругие постоянные измельчаемой частицы. На рис. 3 дана конкретизированная схематизация процесса помола, адаптированная к решению поставленной динамической контактной задачи теории упругости, базирующейся на соотношении (2) и модели рис. 1. Рис. 3. Уточненное моделирование соударения шаров Fig. 3. Detailed modeling of ball collisions В основе рассматриваемой теории измельчения используем теорему об изменении кинетической энергии движущейся со скоростью массы [15, 24], формулу (2) и следующие классические предпосылки [15, 18, 22, 24]: а) материал измельчаемой частицы - однородный, сплошной, изотропный и физически линейный, т.е. подчиняется закону Гука [17, 23, 25], а его прочностные свойства соответствуют хрупкому напряженному состоянию с разными временными сопротивлениями при сжатии и растяжении [25, 26], что характерно, например, для мрамора, гранита, песчаника и т.п. [26, 27]; б) соударяющиеся поверхности предполагаются совершенно гладкими [15]; в) удар является неупругим при нулевом коэффициенте восстановления [24, 28]; г) характер аналитических соотношений между усилиями , и перемещениями , остается неизменным при статической и динамической силах, где - время, т.е. (рис. 3) , (4) а стационарная функция выражается формулой (2); д) физико-механические характеристики , измельчаемого вещества остаются постоянными в условиях динамического нагружения, что допустимо для скоростей соударения [22]; е) установлено [14, 22], что вибрации и колебания, возникающие в движущемся шаре при ударе, не оказывают заметного влияния на его продолжительность и на величину напряжений в материале; ж) время завершения активной фазы ударного импульса (рис. 4) [22] адекватно моменту достижения наибольшей абсолютной динамической деформации (рис. 3), соответствующей, в свою очередь, скорости перемещения суммарной массы : (5) иными словами, в случае реактивная продольная сила удара , (6) а при равнодействующая эллиптической функции контактного давления [9, 10] , , (7) с экстремумом (8) в центре круговой площадки (на рис. 3 она заштрихована) радиусом , (9) становится равной нулю (, , рис. 4); Рис. 4. Обобщенный вид графика ударно-силового импульса Fig. 4. Generalized view of the pulse-force pulse graph з) пренебрегаем собственным весом частицы, по сравнению с весовой составляющей шара-измельчителя (при ); и) в процессе математической аппроксимации напряженного состояния разрушаемого материала учитываются только местные (контактные) деформации, и это допущение является обоснованным, с общепринятой в настоящее время точки зрения, применительно к сферической модели (рис. 1-3) [14, 15, 18, 22]; к) определяющим фактором при оценке общей несущей способности шаровой расчетной схемы измельчаемой частицы из хрупкого вещества согласно [21] является наибольшее главное растягивающее напряжение: , , (10) действующее в радиальном направлении на круговой границе поверхности контакта радиусом (см. рис. 3). Продолжая решение задачи, находим скорость в момент времени точечного соприкосновения шарообразного мелющего тела и сферической поверхности частицы (см. рис. 3), применяя для неупругого удара по допущению в) теорему о равенстве количества движения взаимодействующих тел до и после соударения [24, 28]: (11) откуда (12) где - искомый внешний статический нагрузочный параметр, тождественный силе измельчения , когда ; , , - соответственно плотность [26], объем [29] и вес стального шара-измельчителя (см. рис. 3), (13) - ускорение силы тяжести. Для определения максимального сближения (см. рис. 3) воспользуемся известным оригинальным способом Рейли-Юнга [14, 30, 31], основанным на законе сохранения механической энергии [24, 28] , (14) где - кинетическая составляющая выражения (14) от массы ; - потенциальная энергия локальной деформации разрушаемой частицы (см. рис. 3). В момент равновесия компонента , а скорость , и, следовательно, кинетическая энергия . При (см. рис. 4) система взаимодействующих шаров получает экстремальное перемещение и соответствующую потенциальную энергию , когда и . Поэтому уравнение (14) можно записать в такой форме для этих двух граничных значений времени и : (15) или , (16) следовательно, . (17) Учитывая нелинейный характер функции (4) в зависимости от [32], раскрываем равенство (17), используя (4) при , (12), (13) и [14, 22, 29, 30]: (18) Из энергетического баланса (18) получаем . (19) Имея в соответствии с (19) сближение , находим согласно (4), (8) и зависимости [21] (20) искомые максимальные динамические параметры (см. рис. 3): , (21) , (22) для продолжительности контакта (см. рис. 4) [15, 22] (23) где - упругая константа материала частицы , ; (24) - физико-механический функциональный коэффициент , ; (25) - геометрическая характеристика модели рис. 3, , ; (26) , , - значения гамма-функций или эллиптических интегралов второго рода [29, 33]: ; ; (27) . Для определения статического усилия измельчения используем прежде всего формулу (10), при и , (28) где - максимальное реактивное давление (см. рис. 3), (29) от равнодействующей силы [21,26], учитывая обозначения коэффициентов (24) и (26). С целью математического описания предельного состояния измельчаемого шара применяем критерий прочности Галилея-Лейбница [21] в виде равенства [25-27] (30) с наибольшим по модулю эквивалентным статическим напряжением , которое широко распространено в проектно-технических расчетах для весьма хрупких материалов, по-разному сопротивляющихся сжатию и растяжению (, см. предпосылку а)). В справочно-нормативной литературе [26] формула (30) известна как гипотеза максимальных растягивающих [21] нормальных напряжений (первая теория прочности). Подстановка выражений (28), (29) в условие несущей способности (30) приводит к алгебраическому уравнению относительно , (31) из решения которого (32) Повторяя алгоритм (28)-(31), для ударного действия наибольшей инерционной силы (см. рис. 4), соответствующей статической нагрузке (32), будем иметь, учитывая (22), (24)-(26), главное функциональное соотношение решаемой задачи, представляющее собой критерий разрушения твердого хрупкого материала при его измельчении по гипотезе Галилея-Лейбница [25-27]: . (33) Для граничного случая , (34) аналогичного (30), находим в общем виде предельное значение скорости , меньше которой частица заданного радиуса не будет измельчаться с учетом допущения д) [22]: , (35) откуда после исключения на основании (25) и (32), (36) Последующую расчетно-исследовательскую апробацию условий (33), (36) целесообразно выполнять в конкретном численном виде, зависящем от скорости и радиального размера сферической модели разрушаемого материала. С этой целью рассмотрим измельчаемое вещество из высокопрочного гранита, имеющего физико-механические константы [26, 27, 34] МПа МПа , МПа , при скорости соударения , изменяющейся в пределах , (37) радиусах поверхностей контакта мм мм (38) и удельном весе (13) мелющих шаров (см. рис. 3). Уточняем также допущение з) условиями , (39) учитывающими предполагаемые конструктивно-технологические особенности процесса помола. Воспользовавшись приведенными количественными данными, - вычисляем постоянную в соответствии с (24): ; (40) - находим зависимость , руководствуясь (32) и (40): (41) - определяем согласно (26), (38) значения функции , когда мм (табл. 1). . (42) Таблица 1 Результаты расчета геометрического параметра Table 1 The results of calculation the geometric parameter , мм 0 0,02 0,04 0,06 0,08 0,10 0,12 , 50,05 25,05 16,72 12,55 10,05 8,38 Далее, непосредственно проверяя соотношение (36) для каменных материалов с очень малыми величинами коэффициента (см. результат (40)) и предел прочности на растяжение, можно убедиться, что в подкоренном выражении (36) первое слагаемое всегда многократно - на несколько порядков меньше второго, т.е. (43) и на основании (43) формула (36) принимает упрощенный вид: (44) после подстановки физико-математических констант , и замены суммарной кривизны контактирующих поверхностей (см. рис. 3) в соответствии с (26) при выполнении условия (см. (39)). Потребовав соблюдения равенства , (45) возможно количественно оценить тот минимально допустимый размер частицы, при котором обеспечивается механико-математическая корректность и точность выведенных аналитических зависимостей (33), (44), содержащих физические постоянные , , . Заменяем предел прочности в общей формуле (33) и коэффициенты , для гранита, используя выражения (25), (40), (41). В итоге получаем два функциональных соотношения: , (46) откуда при и мм имеем (47) а из равенства (45) или находим (48) Графики функций , , , построенные по данным табл. 1, 2, 3, 4, изображены на рис. 5 и 6. Соответствующие расчеты проведены по представленным ниже формулам, в основе которых зависимости (46), (47): при мм, (см. табл. 1 и 2) (49) при мм, (см. табл. 1 и 2) (50) при (см. табл. 1 и 3) (51) при (см. табл. 1 и 3) (52) при (см. табл. 4) . (53) Таблица 2 Результаты вычислений функций (49), (50) Table 2 The results of the calculations of the functions (49), (50) , 0 200 400 600 800 1000 1200 , 0 2,968 3,916 4,606 5,168 5,65 6,078 , 0 4,092 5,399 6,35 7,125 7,79 8,379 а б Рис. 5. Графическая иллюстрация функциональных формул (49)-(52): а - график зависимости эквивалентного напряжения от скорости ; б - график зависимости эквивалентного напряжения от радиуса частицы Fig. 5. Graphs of the functional formulae (49)-(52): a - graph of the dependence of the equivalent stress on the speed ; b - graph of the dependence of the equivalent stress on the radius of the particle Таблица 3 Численная информация о функциональных зависимостях (51), (52) Table 3 The calculations of the functional dependences (51), (52) 0 0,02 0,04 0,06 0,08 0,10 0,12 , 0 2,943 3,38 3,665 3,881 4,057 4,208 , 0 2,135 2,452 2,658 2,815 2,943 3,052 Таблица 4 Расчетные данные к функции (53) Table 4 Calculation data for the function (53) , мм 0,002 0,004 0,006 0,008 0,01 0,012 0,014 0,016 0,018 0,02 , 618,9 437,6 354,3 309,5 276,8 252,7 234 218,9 206,4 195,8 Рис. 6. Гиперболический график зависимости (53) Fig.6. The hyperbolic graph of the dependence (53) Завершая данную статью, дополняем аналитические соотношения (19)-(23), (32), (33), (44) полуэмпирическими формулами, аппроксимирующими с высокой точностью ударно-силовой импульс в промежутке времени (см. рис. 4) [22]: (54) где , , - динамические характеристики (19), (21), (23). В то же время следует подчеркнуть, что параметры (54) и время соударения (23) не имеют принципиального значения с практической точки зрения при изучении процесса помола материала. Однако необходимость их использования может возникнуть в случае уточненного математического моделирования разновидности удара, рассмотренной в этой статье (см. рис. 1 и 3), базируясь на фундаментальных положениях механики упругодеформируемого твердого тела [20] и колебательно-волновых теориях [14, 15, 21]. Комплексный анализ проведенных расчетно-теоретических исследований позволяет сделать следующие выводы: 1. В пределах ограничения скорости соударения [22] получено оригинальное инженерно-техническое решение динамической задачи теории упругости о контактном взаимодействии деформируемой неподвижной шаровой модели измельчаемого хрупкого материала с двумя абсолютно твердыми сферическими телами, имеющими заданную массу и максимальную кинетическую энергию (см. рис. 1). 2. Разработанная инновационная механико-математическая модель доведена до простых функциональных соотношений (24), (26), (32), (33), (44), наглядно проиллюстрированных (см. рис. 5, 6) численным примером определения несущей способности шара из гранита. 3. Установлено существование минимально допустимых значений радиуса частицы и скорости мелющего шара (формулы (33), (44), рис. 5), в пределах которых материал не будет разрушаться. 4. Выведенные формулы и алгоритм расчета можно непосредственно использовать при физико-математическом описании и регулировании процесса дробления строительных и рудных материалов измельчительными машинами ударно-вибрационного типа. 5. Степень помола теоретически варьируется и прогнозируется с помощью уравнения (44) по заданному предельному радиальному размеру частицы , соответствующему требуемой скорости соударения (рис. 1, 3 и 6; пример - формула (53)) при условии, что максимально допустимая величина адекватна очень мелкому измельчению твердого вещества, когда порядок оценивается мм (48). 6. Результаты работы можно обобщить и модифицировать применительно к определению мощности и производительности реальной мельницы с подбором оптимальных проектно-технологических параметров: скорости удара , радиуса разрушаемой частицы, диаметра шара-измельчителя и частоты вынужденных колебаний привода-активатора.

About the authors

B M Abdeev

D. Serikbayev East Kazakhstan State Technical University

S R Baigereyev

D. Serikbayev East Kazakhstan State Technical University

G A Guryanov

D. Serikbayev East Kazakhstan State Technical University

References

Statistics

Views

Abstract - 206

PDF (Russian) - 119