CC BY
271
НАУЧНОЕ ИЗДАНИЕ МЕТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 - 48211. Государственная регистрация №0421200025. КБМ 1994-0408
электронный научно-технический журнал
Методика преобразования накопительной функции седиментометра в гранулометрический состав горной породы
# 06, июнь 2013
Б01: 10.7463/0613.0576586
Дозморов П. С., Росляк А. Т.
УДК 519.688
Россия, Национальный исследовательский Томский Политехнический Университет
Эффективность целого ряда процессов порошковой технологии в различных отраслях промышленности, а также качество конечной продукции в значительной степени зависят от точности определения размеров частиц твердых компонентов, что вызывает необходимость совершенствования методов и устройств для определения гранулометрического состава [4].
Для определения гранулометрического состава дисперсных материалов используются множество методов, причем обработка результатов измерений проводится с помощью ряда аппроксимационных зависимостей [2]. В настоящей работе используется один из подходов применения метода Розина - Раммлера - Беннета путем преобразования получаемой информации в аппаратной части с целью создания алгоритма для машинной обработки данных о гранулометрическом составе. Рассмотрим эту задачу на примере метода слоевой седиментации частиц.
Способ седиментации частиц из стартового слоя [1] обеспечивает осаждение с одной высоты всех частиц анализируемой пробы порошка. В результате фиксируются все, даже самые крупные частицы, которые при обычных методах седиментационного анализа успевают достигнуть дна кюветы до начала измерений. Для реализации данного метода используется прибор «Весовой седиментометр ВС-4».
Принципиальная схема весового седиментометра представлена на рис. 1. Частицы анализируемой пробы оседают в жидкости с постоянной скоростью, которая в свою очередь зависит от размера частицы, плотности ее материала, плотности и вязкости жидкости. В процессе осаждения на приемную чашку частицы непрерывно взвешиваются высокочувствительной системой. В каждую секунду времени персональный компьютер фиксирует положение чашки. На рис. 2 приведен пример показаний прибора от времени.
В идеальном случае с каждой осевшей частицей показания прибора должны увеличиваться, однако, как видно из рисунка 2, данная функция не является монотонно возрастающей, что объясняется искажением входной информации различными шумами. Поэтому первым этапом преобразования полученных данных является интерполяция [5].
Рис. 1. Принципиальная схема седиментометра ВС-4
В идеальном случае с каждой осевшей частицей показания прибора должны увеличиваться, однако, как видно из рисунка 2, данная функция не является монотонно возрастающей, что объясняется искажением входной информации различными шумами. Поэтому первым этапом преобразования полученных данных является интерполяция [5].
Рис. 2. Функция накопления данных по времени с прибора «Весовой седиментометр ВС-4»
Р(8) - функция распределения частиц по размерам (8). Данную функцию можно получить разными методами [2].
Для нашего исследования наиболее применимым является метод Розина - Раммлера
- Беннета [6]. По этому методу кривые распределения размера частиц могут быть получены уравнением:
Тогда уравнение (2) примет вид у = ах + Ь - линейное уравнение. В этом уравнении известны пары (х{, у), 1 = 1, п, п - количество экспериментальных данных. Для нахождения а и Ь воспользуемся методом наименьших квадратов.
Функция двух переменных принимает наименьшее значение
1=1
При данных а и Ь сумма квадратов отклонений экспериментальных данных от найденной прямой будет наименьшей [3].
Чтобы найти коэффициенты, находим частные производные функции по переменным а и Ь и приравниваем эти производные к нулю. Полученная система решается методом Краммера и получаем формулы для нахождения коэффициентов по методу наименьших квадратов
Р(8)
(1)
Для нахождения параметра а уравнения (1) дважды логарифмируем:
(2)
Обозначим:
х = 1п8
(3)
п
Г(а, Ь) = ^ (у - (ах1 + Ь))2 ^ шт
(4)
п Е хгУг -Ех Е у
г=1
г=1 г=1
Ап Л2
п
Ех> - Е х-
V г=1 У
г=1
Е у- а Е х
Ь = -*=—
г=1
п
При данных а и Ь функция (4) принимает наименьшее значение.
Как видно из уравнений (3), для нахождения (, у) необходимо знать размер частиц.
Воспользуемся уравнением в условиях динамического равновесия для движущейся частицы [1]:
dw
(6)
РВ - РС - РП
т
dт ’
где
с 2
=
т
Л?
(7)
(8)
(9)
тв
В приведенных выше уравнениях g - ускорение свободного падения частицы под действием силы тяжести; w - скорость осаждения частицы относительно среды; р и ртв-
плотность среды и твердого материала (частицы) соответственно; ^ - проекция поперечного сечения частицы на направление ее движения (площадь миделева сечения); С, - коэффициент сопротивления; т - масса частицы.
Из уравнений (6 - 9) получим:
( \ Р
dw
dт
= g
1 -
V Рте У
2т
I
(10)
Уравнение (10) характеризует взаимодействие сил, в поле которых находится твердая частица или тело.
Если частица имеет шарообразную форму с диаметром 5 и осаждается в поле силы тяжести, то для определения скорости осаждения необходимо ввести условия:
• среда, в которой происходит осаждение, неограниченна;
• осаждению частицы не мешают другие частицы;
• скорость осаждения постоянна.
В соответствии с последним условием ё’/ёт = 0. При введении массы шарообразной частицы с т= (р53 /б)ртв и площади поперечного сечения частицы
/ = лб24 в
уравнение (10) получим:
ё
1
Р
Р
тв /
Эфм2
4ЪРТВ
= 0
(11)
или
%
г \
1 --Р
V Рте J
2
ос
45Рте
(12)
откуда
’ =
ос
У
4(рТ,- р)ё
(13)
3 Р?
В приведенных выше уравнениях 8 - диаметр осаждающейся частицы, м; ртв и р - плотности частицы и среды соответственно, кг/м .
Таким образом, уравнение (13) показывает зависимости скорости осаждения от коэффициента сопротивления среды, который, в свою очередь, зависит от числа
Рейнольдса Яе = ^8 / V, где V - кинематическая вязкость среды осаждения.
Рассматриваемые нами режимы осаждения частиц являются или ламинарным или переходным. Отсюда изменяется коэффициент сопротивления среды.
В случае ламинарного режима осаждения (Яе < 2 ) С = 24/Яе.
Промежуточный режим обтекания в пределах изменения 2 < Яе < 500 характеризуется меньшей зависимостью сопротивления от критерия Рейнольдса:
с = 185/Яеаб .
Применим данную теорию для нахождения размера частиц. Нам известно М = Н / ? , где Н - высота осаждения частиц, ^ - время осаждения частиц одного размера.
Подставим данное равенство в (13) и выразим время ^ . Получим
г =
3Н 2р£
(14)
V 4%8(Ртв - Р)
Подбирая размер частицы (с учетом коэффициента сопротивления среды), найдем такое время осаждения, которое совпадало бы со временем, полученным экспериментально.
Следуя данному методу, получим максимальный размер частиц в данном эксперименте. Поскольку показания прибора содержат скорость осаждения частиц, воспользуемся нормировкой униполярного показателя [7], выражающего только степень наличия некоторого качества (скорости), которое зависит от максимального и минимального показаний прибора.
Далее находим размер частиц для каждого значения показания прибора и времени его получения с учетом нормировки показания прибора. Каждое нормированное показание прибора дважды логарифмируем (у из уравнения (3)). По вышеописанному
методу находим размер частицы, логарифмируем, получаем xt из формулы (3).
Найдя, таким образом, все пары xt у, подставим их в формулы (3) для нахождения а и b .
Таким образом, мы нашли коэффициенты для уравнения Розина-Раммлера-Беннета. Зная коэффициенты уравнения Розина-Раммлера-Беннета, можно найти следующие величины:
- медиану полученного распределения;
- удельную поверхность частиц;
- дифференциальное распределение частиц по размерам;
- интегральную функцию распределения частиц по размерам;
- массовую долю частиц по фракциям в процентном соотношении;
- количество частиц в процентном соотношении.
Представленный метод анализа гранулометрического состава с позиции преобразования получаемой информации с аппаратной части прибора является, на наш взгляд, наиболее применимым в качестве метода для алгоритмизации процедуры и аппаратной обработки данных для данного прибора.
Список литературы
1. Квеско Н.Г., Росляк А.Т. Весовой седиментометр для автоматизированного измерения гранулометрического состава порошков // Заводская лаборатория. Диагностика материалов. 2000. № 7. С. 37-40.
2. Коузов П.А. Основы анализа дисперсного состава промышленных пылей и измельченных материалов. 3-е изд., перераб. Л.: Химия, 1987. 264 с.
3. Линник Ю.В. Метод наименьших квадратов и основы математико-статистической теории обработки наблюдений. Л.: Физматгиз, 1962. 352 с.
4. Росляк А.Т., Зятиков П.Н. Воздушно-центробежная классификация микропорошков. Томск: ТМЛ-Пресс, 2010. 224 с.
5. Щиголев Б.М. Математическая обработка наблюдений. М.: Наука, 1969. 344 с.
6. Rosin P., Rammler E. The Laws Governing the Fineness of Powdered Coal // Journal of the Institute of Fuel. 1933. Vol. 7, no. 31. P. 29-36.
7. Шенк Х. Теория инженерного эксперимента: пер. с англ. М.: Мир, 1972. 381 с.
SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-0408
electronic scientific and technical journal
Transformation technique for accumulative function of a sedimentometer into a grain-size composition of geological material
# 06, June 2013
DOI: 10.7463/0613.0576586
Dozmorov P.S., Roslyak A.T.
Russia, National ResearchTomsk Polytechnic University
The authors consider application sedimentation from the starting layer, which allows one to record large grains in a test sample, unlike other sedimentation methods. The analyzing method for grain-size composition of powders and particles of geological material was considered from a perspective of transformation of obtained information in a hardware component of a device as the most suitable for algorithmization of processing measured data when using the sedimentation analysis from the starting layer. Special attention was paid to the analysis of discrete values of the accumulative function by the linear method with a possibility of implementing this method as software.
Publications with keywords: Reynolds, least squares method, sedimentation, Rosin-Rammler-Bennett, normalizing, drag coefficient, Stokes
Publications with words: Reynolds, least squares method, sedimentation, Rosin-Rammler-Bennett, normalizing, drag coefficient, Stokes
References
1. Kvesko N.G., Rosliak A.T. Vesovoi sedimentometr dlia avtomatizirovannogo izmereniia granulometricheskogo sostava poroshkov [Weight sedimentometer for automated measurement of particle size distribution of powders]. Zavodskaia laboratoriia. Diagnostika materialov, 2000, no. 7, pp. 37-40.
2. Kouzov P.A. Osnovy analiza dispersnogo sostavapromyshlennykhpylei i izmel'chennykh materialov [Basis of the analysis of the dispersed composition of industrial dusts and grinded materials]. Leningrad, Khimiia, 1987. 264 p.
3. Linnik Iu.V. Metodnaimen'shikh kvadratov i osnovy matematiko-statisticheskoi teorii obrabotki nabliudenii [The method of least squares and the fundamentals mathematical-statistical theory of observation processing]. Leningrad, Fizmatgiz, 1962. 352 p.
4. Rosliak A.T., Ziatikov P.N. Vozdushno-tsentrobezhnaia klassifikatsiia mikroporoshkov [Air-centrifugal classification of micropowders]. Tomsk, TML-Press, 2010. 224 p.
5. Shchigolev B.M. Matematicheskaia obrabotka nabliudenii [Mathematical treatment of observations]. Moscow, Nauka, 1969. 344 p.
6. Rosin P., Rammler E. The Laws Governing the Fineness of Powdered Coal. Journal of the Institute of Fuel, 1933, vol. 7, no. 31, pp. 29-36.
7. Schenk H., Jr. Theories of Engineering Experimentation. New York, McGraw-Hill, 1961. (Russ. ed.: Shenk Kh. Teoriia inzhenernogo eksperimenta. Moscow, Mir, 1972. 381 p.).
