Рус Eng Cn Translate this page:
Please select your language to translate the article


You can just close the window to don't translate
Library
Your profile

Back to contents

Software systems and computational methods
Reference:

Defining Morphometric Properties of the Micro-Disperse Aluminum Oxide System by the Image Analysis Method

Ignatova Anna Mikhailovna

PhD in Technical Science

leading researcher at Perm National Research Polytechnic University

614000, Russia, Perm Krai, Perm, str. Komsomol'ski prospect, 29

iampstu@gmail.com
Zemlyanova Marina Aleksandrovna

Doctor of Medicine

professor at Perm National Research Polytechnic University

614045, Russia, Permski Krai, Perm, str. Monastyrskaya, 82, room No. 201

zem@fcrisk.ru
Stepanko Mark Sergeevich

Bachelor of Arts at Perm National Research Polytechnic University

614990, Russia, Perm Region, Perm, Bukirev's str., 15

mmm2920@mail.ru
Ignatov Mikhail Nikolaevich

Doctor of Technical Science

professor at Perm National Research Polytechnic University

614000, Russia, Perm Region, Perm, str. Komsomol'ski prospect, 29

iampstu@gmail.com

DOI:

10.7256/2454-0714.2017.3.23526

Received:

06-07-2017


Published:

06-10-2017


Abstract: The subject of the research is the morphometric properties of Al2O3 (aluminum oxide) micro-disperse system. The authors pay special attention to the fractional composition and particle sphericity factor for defined particle groups constituting the disperse system. The authors state that image analysis allows to use the sphericity factor for evaluating elements of the particle surface combined with the interpretation of a real geometrical shape of individual particles. This allows to better define the dispersion factor as an important property of a disperse system because this parameter may be a measure of substance activity in a disperse system. Digital images for analysis ere obtained using a scanning electronic microscope with a high resolution S3400N HITACHI (made in Japan). Images were analyzed using a universal software ImageJ-FiJi. The object of the research was the powder Al2O3 of Sigms-Aldrich brand (article 265497, α-phase, particle average size is 10 micron). The powder was made at mills with spherical grinding bodies. The main conclusion of the research is that image analysis is a powerful tool that can be used to define physical properties of disperse systems. For the first time in the academic literature the authors focus on analyzing angles created by particle brim peaks. By calculating the dispersion factor, the authors have demonstrated that data obtained during particle image analysis allows to increase efficiency in defining the real factor of particle dispersity with the average validity of 5%. 


Keywords:

image analysis, morphology, sphericity factor, granulometry, surface texture, disperse systems, Aluminium oxide, dispersion factor, microscopy, metrology

This article written in Russian. You can find original text of the article here .

Традиционно основным свойством дисперсных (порошковых) систем является гранулометрический состав, поскольку соотношение частиц разного размера определяет их технологические свойства (насыпная плотность, текучесть, площадь поверхности и т.д.). Исследования последних лет [1-3] обращают внимание на взаимосвязь эксплуатационных свойств дисперсных систем и изделий из них (абразивная способность, твердость, прочность, водопоглощение, пористость и т.д.) не только с гранулометрическим составом, но и с формой частиц. В связи с чем, актуальной задачей является выявление морфометрических характеристик дисперсных систем.

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

В настоящее время гранулометрический (фракционный) состав порошковых материалов определяют различными методами (ситовой анализ, анализ изображений, лазерный анализ и т.д.) [4]. Однако, точность и достоверность результатов этих методов применительно к системам разного размерного уровня и фракционного диапазона различна. По мнению большинства исследователей, для систем, состоящих из частиц, средний размер которых более 50-100 мкм, все эти методы дают хороший и достоверный результат. Снижение достоверности и точности возникает при оценке гранулометрического состава порошка, состоящего из частиц размером менее 50 мкм. В публикациях за последние пять лет [5-7] для таких целей все чаще используют анализ изображений, полученных с помощью оптических и электронных микроскопов.

Дисперсные системы Al2O3 со средним размером частиц менее 50 мкм, то есть относящиеся к микродисперсным, имеют широкое применение в промышленности. Они являются основным компонентом при изготовлении композитов, керамических материалов, защитных покрытий, абразивных, полирующих и реакционных составов, а также косметических средств и некоторых видов материалов для медицины (стоматологические имплантаты, биосовместимые протезы и т.д.) [8-10]. Микродисперсные системы Al2O3 могут быть получены разными способами, в том числе, в дробильных устройствах и мельницах разного типа, поэтому их морфологический и гранулометрический состав может варьироваться в широком диапазоне, в зависимости от способа получения.

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

Часто гранулометрическую характеристику порошковых материалов, представляют в виде гистограммы, в которой высота столбцов определяется величиной массовой доли частиц определенного размера от общего количества в процентах. Статистические данные о гранулометрическом составе принято представлять в виде графика логарифмического нормального распределение, который строится на основе Гауссовской функции, аргумент в которой представлен логарифмом диаметров частиц. Гауссовская функция не является единственным и универсальным инструментом интерпретации дисперсного состава. В публикациях [11-12] отмечается, что применимость того или иного распределения для описания состава конкретной дисперсной системы определяется главным образом способом получения этой системы. Для порошков, полученных в валковых и щековых дробилках и стержневых мельницах используют формулу распределения Годена [13] Материалы, измельченные в шаровой мельнице, хорошо поддаются описанию с помощью распределения Розина–Раммлера [14]. Для дисперсных систем, состоящих из частиц менее 30 мкм, предпочтительнее использовать логарифмическое распределение А.Н. Колмагорова [15]. Распределение А.Н. Колмагорова также широко используют при описании данных полученных при анализе изображений [16].

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

Для получения морфометрической характеристики дисперсной системы оксида алюминия (Al2O3) использовали цифровые изображения, полученные на сканирующем электронном микроскопе высокого разрешения S3400N «HITACHI» (Япония). Анализ изображений проводился универсальным программным обеспечением ImageJ-FiJi. В качестве объекта исследования использовался порошок Al2O3 марки Sigms-Aldrich (артикул 265497, α-фаза, средний размер частиц 10 мкм), полученный в мельницах с шарообразными мелющими телами.

Как известно [17], изображения РЭМ могут быть получены на двух режимах: топографического и композиционного контраста. Топографический контраст позволяет получить информацию о состояние поверхности частиц, такие изображения хорошо использовать для трехмерных реконструкций. Композиционный контраст, в особенности, когда речь идёт о дисперсных системах высокой химической частоты, наилучшим образом подходит для изучения геометрических параметров частиц. Для анализа была получена серия исходных изображений на обоих режимах в количестве 20 штук, пример, представлен на рис. 1.

а

б

Рис. 1. Изображения частиц дисперсной системы Al2O3 (РЭМ):

а – режим химического контраста, б – режим топографического контраста

Достоверность геометрических измерений при анализе изображений зависит от качества его подготовки. Основным и первоначальным комплексом действий по подготовке изображения (рис. 2, а) является удаление фона, реликтов контраста и перевод в двух цветную (бинарную) цветовую схему [18]. На полученном изображение, часто присутствуют частицы, границы которых выходят за рамки изображения, и изображения частиц, расположенных таким образом, что при переводе в бинарный вид формируются не соответствующие действительности двумерные проекции (рис. 2, б). Такие недочеты влияют на результат анализа повышая его погрешность. Часть таких реликтов может быть скорректирована программными средствами избирательно (вручную), о чем авторы указывали в предшествующих работах [18]. В результате можно получить изображение, анализ которого, позволит получить более достоверные данные как о диапазоне фракционного состава, так и о форме частиц (рис. 2, в).

а

б

в

Рис. 2. Подготовка изображений частиц дисперсной системы Al2O3 (РЭМ) для анализа: а – исходное изображение, б – черновая обработка, в – чистовая обработка

Результат анализа частиц, по двумерным проекциям, полученных при обработке результатов РЭМ были получены данные о размере и форме 2361 отдельных частиц. Совокупный анализ данных об отдельных частицах позволил установить гранулометрический (рис. 3, а) и морфологический составы (рис. 4, а), логнормальные распределения этих показателей (рис. 3,б и 4,б) и провести морфометрический анализ, показывающий как морфологический, состав изменяется в разных размерных группах (рис. 5).

Гранулометрический состав, полученный по результатам статистической обработки величин диаметров эквивалентных сфер, для изученной дисперсной системы Al2O3 показал, что основная доля частиц (46,50% от общего количества частиц) имеет размер 3-6 мкм эквивалентного диаметра. Медиана логарифмического нормального распределения составила 5,65±0,05 мкм (рис. 3, б, табл.1).

Морфологический состав, характеризующий распределение частиц в соответствие с коэффициентом сферичности Ксф [19] (1 – идеальная сфера, 0,1 – игольчатая форма), показал, что основная доля частиц (60,41 % от общего количества частиц) имеет коэффициент сферичности 0,8-0,9. Медиана логарифмического нормального распределения по коэффициенту сферичности составила 0,835±0,001 мкм (рис. 4, б, табл. 1).

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

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

а

б

Рис. 3. Гранулометрический состав дисперсной системы Al2O3: а – гистограмма, б –логнормальное распределение

а

б

Рис. 4. Морфологический состав дисперсной системы Al2O3:

а – гистограмма, б –логнормальное распределение

Рис. 5. Морфологический состав частиц в разных размерных группах дисперсной системы Al2O3

Таблица 1

Анализ данных морфометрических характеристик частиц дисперсной системы Al2O3

Статистический показатель

Размер частиц, мкм

Коэффициент сферичности частиц

Количество обработанных изображений

20

20

Медиана логнормального распределения

5,65±0,05

0,835±0,001

Среднее значение

5,65

0,83

Стандартное квадратичное отклонение

2,5

0,06

Коэффициент вариации

44,24

7,22

Минимальное значение

0,16

0,40

Максимальное значение

15,22

0,98

Коэффициент сферичности является обратной величиной коэффициента формы (Кф), количественного показателя, который широко используется для оценки формы частиц дисперсных абразивных материалов и предусмотрен ГОСТ 9206-80 [20]. Согласно принятого стандарта, коэффициент формы определяет имеет ли частица изометричную (Кф=1,0-1,4), пластинчатую (Кф=1,4-2,2) или переходную (Кф=2,2-3,0) форму. Статистическим анализом (табл. 2) установлено, что частицы дисперсной системы Al2O3 в большинстве своем (94,8% от общего количества частиц) обладают изометричной формой. Учитывая ранее выявленные морфометрические показатели коэффициента сферичности и соответствующие им геометрические формы (табл. 3) [21], можно утверждать, что основная масса частиц дисперсной системы Al2O3 имеет форму объемных выпуклых многогранников.

Таблица 2

Анализ данных о коэффициенте формы частиц дисперсной системы Al2O3

Статистический показатель

Коэффициент формы

Количество обработанных изображений

20

Количество обработанных измерений

2360

Медиана логнормального распределения

1,181±0,002

Среднее значение

1,2

Стандартное квадратичное отклонение

0,11

Коэффициент вариации

9,14

Минимальное значение

1,01

Максимальное значение

2,5

Таблица 3

Соответствие коэффициента сферичности формам геометрических тел

Тип фигуры

Кф

Ксф

Тип фигуры

Кф

Ксф

сфера

1

1

призма а/2а/3а

1,38

0,725

икосаэдр

1,07

0,937

цилиндр h=5d

1,45

0,691

додекаэдр

1,1

0,912

тетраэдр

1,49

0,67

цилиндр h=1,5d

1,16

0,86

диск h=r/3

1,68

0,594

октаэдр

1,18

0,846

цилиндр h=10d

1,72

0,58

цилиндр h=d/2

1,21

0,827

цилиндр d=h

2,28

0,438

куб

1,24

0,806

диск h=r/10

3,1

0,323

призма a/a/2a

1,3

0,767

диск h=r/15

4,55

0,22

призма а/2а/2а

1,31

0,761

Выявленная характеристика формы является объемной, что позволяет определить совокупную характеристику дисперсной системы – фактор дисперсности (fдисп) [22]. Фактор дисперсности определяется отношением площади поверхности к объему твердого тела частицы. В работе [23] отмечается, что стабильность тела увеличивается с увеличением размера, а фактор дисперсности является обратной величиной стабильности. Иными словами, чем мельче частицы вещества, тем меньше его стабильность. Причем это правило распространяется на частицы не только в нанодисперсном, но и в микродисперсном диапазоне. Разумеется, можно определить величину фактора дисперсности по данным о диаметре частиц и параметрам эквивалентных сфер. Однако, в зависимости от формы многогранника величины его объема и площади поверхности могут отличатся от эквивалентных сфер. Следовательно, необходимо получить максимально доступное количество информации о состояние поверхности и форме частиц при изучении изображений, для того чтобы повысить точность определение фактора дисперсности.

Наблюдение частиц дисперсной системы в режиме топографического контраста методом РЭМ, позволило выявить, что на поверхности частиц в размерном диапазоне 3-5 мкм, имеются скопления более мелких частиц размером 0,1-0,4 мкм (рис. 6, а). При этом текстура поверхности по данным трехмерной реконструкции [18] приобретает шероховатость, с высотой рельефа 0,5-1 мкм (рис. 6, б, в). Анализ изображений 23 граней частиц размером 5-15 мкм показал, фракционный и морфологический состав и то насколько плотно частицы размером 0,1-0,4 мкм, распределены по их поверхности (табл. 4).

а

б

в

Рис. 6. Изучение поверхности микрочастиц Al2O3: а – исходное изображение, б – трехмерная реконструкция, в – профиль участка с частицами диаметром 0,35 мкм

Таблица 4

Анализ данных о диаметре частиц формирующих микрорельеф

Статистический показатель

Величина показателя

Количество обработанных изображений

20

Количество обработанных изображений граней частиц

23

Среднее количество частиц на грани

7,2

Грануломтерический состав, мкм, %:

0,1-0,2

52

0,2-0,3

25

0,3-0,4

23

Медиана логнормального распределения, мкм

0,19±0,01

Среднее значение, мкм

0,22

Стандартное квадратичное отклонение

0,08

Коэффициент вариации

37,89

Минимальное значение, мкм

0,12

Максимальное значение, мкм

0,39

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

Характеристика, которая отличает одни геометрические тела от других по мимо формы и количества граней, является величина двугранного угла образованного при вершине многогранника [24]. Углы, образованные гранями частиц, измеряли визуально по изображениям отдельных фрагментов (рис. 7), а затем провели автоматизированный анализ с условием приведения формы частиц к эквивалентному эллипсу и с вычислением угла между радиусом и большой полуосью эквивалентного эллипса. Статистический анализ автоматизированных измерений величины угла эллептичности показал, что все измерения можно поделить на две группы, усредненная величина которых составляет 59,84 и 128,12 º. Среди Платоновых геометрических тел близкими параметрами, а именно угол при вершине 125,26º обладает кубооктаэдр (рис. 8, а), относящийся к правильным многогранникам, которому соответствует полуправильный многогранник или каталаново тело – ромбододекаэдр (рис. 8, б) с углами при вершинах 70,53 и 109,47 º [25].

В совокупности результаты оценки рельефа и определения наиболее близкого геометрического тела, вычисление объема и площади поверхности которых известно, для описания частиц предложено использовать для при расчете площади поверхности (1) и объема частиц (2) дисперсной системы Al2O3 приведенные коэффициенты 1,1 и 1,3 соответственно:

S=1,1·Sэкв (1)

V=1,3·Vэкв (2)

где Sэкв и Vэкв – площадь поверхности и объем эквивалентного ромбододекаэдра, у которого длина ребра равна радиуса эквивалентной сферы.

Рис. 7. Измерение углов кромок частиц дисперсной системы Al2O3

Uniform polyhedron-43-t1.pnghttp://www.tryunity.net/uploads/2/3/8/9/23893634/rhombicdodecahedron_orig.jpg

а б

Рис. 8. Геометрические тела близки по своим параметрам к параметрам частиц дисперсной системы Al2O3: а - кубооктаэдр, б – ромбододекаэдр

При подсчете указанных величин был вычислен фактор дисперсности частиц исследуемой системы Al2O3, который оказался выше фактора дисперсности, вычисляемого по параметрам эквивалентных сфер в 1,2 раза (рис. 9, табл. 5).

Таблица 5

Анализ данных факторе дисперсности частиц Al2O3

Статистический показатель

По диаметру эквивалентных сфер

По данным анализа изображений (эквивалентный ромбододекаэдр)

Медиана логнормального распределения

1,12±0,03

1,17±0,03

Среднее значение

1,40

1,46

Стандартное квадратичное отклонение

1,70

1,77

Коэффициент вариации

82,33

82,33

Минимальное значение

0,39

0,41

Максимальное значение

35,84

37,20

а

б

Рис. 9. Фактор дисперсности частиц: а – по величине диаметра эквивалентных сфер, б – по уточненным данным с помощью анализа изображений (эквивалентные ромбододекаэдры)

Таким образом, методом анализа цифровых изображений, были получены морфометрические характеристики микродисперсной системы Al2O3. Установлено, что основная доля частиц (46,50% от общего количества частиц) имеет размер 3-6 мкм эквивалентного диаметра, 60,41 % от общего количества частиц имеет коэффициент сферичности 0,8-0,9. На поверхности частиц в размерном диапазоне 3-5 мкм, имеются скопления более мелких частиц размером 0,1-0,4 мкм, со средней частотой 7 частиц на одну грань. Анализ углов образованных вершинами кромок частиц позволил установить, что наиболее близка геометрическое тело к реальной форме частиц – ромбододекаэдр. По данным о микрорельефе предложены коэффициенты для уточненного расчета объема и площади поверхности частиц, равные 1,3 и 1,1 соответственно. Расчет фактора дисперсности показал, что данные полученные при анализе изображений частиц, позволяют повысить эффективность в определение реального фактора дисперсности частиц в среднем на 5%.

References
1. Stachowiak, G. W. Numerical characterization of wear particle morphology and angularity of particles and surfaces. In New Directions in Tribology. Plenary and Invited Papers of the First World Tribology Congress (Hutchings, I. M., ed.). MEP Publications, London, 1997, 371–389.
2. Stachowiak, G. B. and Stachowiak, G. W. The effects of particle characteristics of three-body abrasive wear. In Proc. 9th Nordic Symposium on Tribology NORDTRIB2000. Porvoo, 2000, 2, 434–443.
3. Stone V, Nowack B, Baun A, van den Brink N, Kammer F, Dusinska M, Handy R, Hankin S, Hassellöv M, Joner E, Fernandes TF: Nanomaterials for environmental studies: Classification, reference material issues, and strategies for physico-chemical characterization. Sci Total Environ 2010, 408(7):1745-1754.
4. Naumov S.V., Kanina A.E., Ignatova A.M., Ignatov M.N.O FRAKTsIONNOM SOSTAVE SVAROChNYKh FLYuSOV Nauchno-tekhnicheskii vestnik Povolzh'ya. 2013. № 2. S. 166-169.
5. Z. Zhang, J. Yang, L. Ding, and Y. Zhao, “An improved estimation of coal particle mass using image analysis,” Powder Technology, vol. 229, pp. 178–184, 2012.
6. Zelin Zhang, Jianguo Yang, Lihua Ding and Yuemin Zhao, An improved estimation of coal particle mass using image analysis, Journal of Powder Technology, vol.229.pp 178-184, 2012.
7. TASDEMIR A., OZDAG H., ONAL G., 2011. Image analysis of narrow size fractions obtained by sieve analysis. An evaluation by log-normal distribution and shape factors. Physicochemical Problems of Mineral Processing, 46, 95–106.
8. Christel, P.; Meunier, A.; Dorlot, J-M. (1988) Biomechanical Compatibility And Design Of Ceramic Implants For Orthopedic Surgery, In: Bioceramics: Material Characteristics Versus In Vivo Behavior. Ducheyne, P. & Lemons, J. E. Annals Of New York Academy Of Sciences,Vol. 523 pp. 234-256.
9. Piconi, C.; Condo, S.C.; Kosmač, T. Chapter 11. Alumina and zirconia-based ceramics for load-bearing applications. In Advanced Ceramic for Dentistry; Shen, J., Kosmač, T., Eds.; Butterworth-Heinemann (Elsevier): Wathmann, MA, USA, 2014; pp. 219–253.
10. Denry, I.; Holloway, J.A. Ceramics for dental applications: A review. Materials 2010, 3, 351–368.
11. Kouzov P. A. Osnovy analiza dispersnogo sostava promyshlennykh pylei i izmel'chennykh materialov. – L.: Khimiya, 1987. – 264 s.
12. Shabanova N. A., Popov V. V., Sarkisov P. D. Khimiya i tekhnologiya nanodis-persnykh oksidov. – M.: IKTs Akademkniga, 2006. – 310 c.
13. Linch A.Dzh. Tsikly drobleniya i izmel'cheniya. Modelirovanie, optimizatsiya, proektirovanie i upravlenie: Per. s angl.. — M.: Nedra, 1981. – C. 343.
14. Azarov V. N. [i dr.]. Metodika mikroskopicheskogo analiza dispersnogo sostava pyli s primeneniem personal'nogo komp'yutera (PK) / // Zakonodatel'naya i prikladnaya metrologiya. – 2004. – № 1. – S. 46-48.
15. Kolmogorov, A.N. O logarifmicheski normal'nom zakone raspredeleniya razmerov chastits pri droblenii / A.N. Kolmo-gorov // Dokl. AN SSSR. – 1941. – T. 31. – № 2. – S. 99.
16. Faria, Pons, Feyo De Azevedo, Rocha, & Vivier, “Quantification of the Morphology of Sucrose Crystals by Image Analysis,” Powder Technology, 133, 2003.
17. Mironov V. L. Osnovy skaniruyushchei zondovoi mikroskopii: ucheb. posobie dlya studentov starshikh kursov vysshikh uchebnykh zavedenii Rossii. Nizhnii Novgorod: Institut fiziki mikrostruktur, 2004. – 110 s
18. Ignatova A.M., Vereshchagin V.I. PRIMENENIE METODA ANALIZA IZOBRAZhENII V ISSLEDOVANII I STATISTIChESKOI OTsENKE PARAMETROV ChASTITs TVERDOI SOSTAVLYaYuShchEI SVAROChNYKh AEROZOLEI SILIKATNOGO I OKSIDNOGO SOSTAVA Vestnik Permskogo natsional'nogo issledovatel'skogo politekhnicheskogo universiteta. Mashinostroenie, materialovedenie. 2017. T. № 1. S. 41-57.
19. Belenko O.A. Vliyanie razmerov i formy chastits na svoistva atmosfernykh aerozolei [Tekst] / O.A. Belenko // Geo-Sibir'-2006. T.3. Monitoring okruzhayushchei sredy, geoekologiya, distantsionnye metody zondirovaniya Zemli i fotogrammetriya. Ch. 1: sb. materialov mezhdunar. nauch. kongr. «Geo-Sibir'-2006».-Novosibirsk: SGGA, 2006.-S. 163-167.
20. GOST 9206-80 POROShKI ALMAZNYE. TEKhNIChESKIE USLOVIYa
21. Gavrilova N. N. Mikroskopicheskie metody opredeleniya razmerov chastits dis-persnykh materialov: ucheb. posobie / N. N. Gavrilova, V. V. Nazarov, O. V. Yarovaya. – M.: RKhTU im. D. I. Mendeleeva, 2012. – 52 s.
22. Vakalova T.V., Khabas T.A., Erdman S.V., Vereshchagin V.I. Praktikum po osnovam tekhnologii tugoplavkikh nemetallicheskikh i silikatnykh materia-lov. – Tomsk: Izd. TPU, 1999. 160 s.
23. Vereshchagin V.I., Kozik V.V., Borilo L.P., Pogrebenkov V.M., Syryamkin B.I. Polifunktsional'nye neorganicheskie materialy na osnove prirodnykh i iskusstvennykh soedinenii. Tomsk: izd-vo Tom. un-ta,-2002.-359 s.
24. I.Kh. Sabitov,. Ob''emy mnogogrannikov. — M.:MTsNMO, 2002. — 32 s.
25. I. Kh. Sabitov, Ob''em mnogogrannika kak funktsiya ego metriki, Fundam. prikl. matem. 1996. T. 2, No. 4. S. 1235—1246.
26. Greg S., Sing K. Adsorbtsiya, udel'naya poverkhnost', poristost'. 2-e izd. M.: Mir, 1984.