CC BY-NC
15Метод экспресс-анализа жидких фасованных продуктов
Белозеров Валерий Владимирович
доктор технических наук
профессор, Донской государственный технический университет, Генеральный директор, ООО "НПТ
Центр ОКТАЭДР"
344091, Россия, Ростовская область, г. Ростов-на-Дону, ул. Капирская, 22
Статья из рубрики "Физика"
Аннотация.
Принимая во внимание тот факт, что проблемам защиты потребителей 17 апреля 2017 года было посвящено заседание президиума Государственного совета (по вопросу «О национальной системе защиты прав потребителей»), и статья соавторов,
опубликованная в № 1 за 2016 год, активно цитируется и дискутируется в научном сообществе, редакционный совет принял решение опубликовать первоначальную работу автора, т.к. именно в ней было найдено решение, которое может составить техническую основу национальной системы. В статье расширена библиография, что позволяет убедиться в успешной адаптации метода экспресс-контроля жидких фасованных продуктов без вскрытия тары. В статье описан экспресс-метод весовой электрометрии и компаративного анализа с «образами эталонов», для идентификации качества и выявления контрафактной продукции.Объектом исследования явились эталоны моторного масла Castrol Magnatec и образцы контрафактной продукции, на предмет регистрации их электрических и механических характеристик, а также параметры и конструкции тары, в которую они расфасованы. Результаты исследования показали, что предлагаемый метод и средства его реализации могут осуществить абсолютную защиту потребителя от некачественной и поддельной жидкой фасованной продукции.
Ключевые слова: динамическая вязкость, кинематическая вязкость, плотность, качество продукции, сравнительный анализ, электрометрия материалов, экспресс-метод, диэлектрическая проницаемость, измеритель иммитанса, емкостный датчик
DOI:
10.7256/2453-8884.2018.2.25998
Дата направления в редакцию:
21-04-2018
Дата рецензирования:
25-04-2018
Введение
Задачи любых моторных масел, независимо от типа - смазывание сопряженных деталей и удаление продуктов износа. Очевидно, что нельзя создать универсальное масло, способное сохранять набор эксплуатационных свойств в широком диапазоне температур. При низких температура оно будет становиться более густым, а при высоких температурах его текучесть увеличивается. Диапазон температуры смазки при эксплуатации может достигать 150 градусов Цельсия, в зависимости от места, скорости и интенсивности циркуляции.
Учитывая это, автопроизводители вычисляют оптимальные характеристики моторных масел, которые должны обеспечить максимально возможный КПД силового агрегата при его минимальном износе в нормальных для данного двигателя условиях эксплуатации. При этом наиболее значимым параметром принято считать класс вязкости моторного масла, иными словами его способность оставаться на поверхности деталей и сохранять при этом текучесть. Поскольку с изменением температуры вязкость меняется, ассоциацией автомобильных инженеров США (SAE) разработана и принята
классификация по вязкости. Согласно данной классификации автомасла делятся на два класса вязкости: зимние (цифра с буквой W, маловязкие) - 0W, 5W, 10W, 15W, 20W, 25W, и летние (без буквенного кода, высоковязкие) - 20, 30, 40, 50, 60, и на три категории: для бензиновых двигателей («S»), для дизелей («С») и энергосберегающие («ЕС»), При этом, качественные характеристики масла определяет также классификация асеа (европейский аналог арі), пришедшая в 1995 году на смену устаревшей системе. Стандарты данной классификации составляют три категории и 12 классов :
А/В - бензиновые и дизельные двигатели легковых автомобилей, фургонов, микроавтобусов (А1/В1-12, АЗ/ВЗ-12, АЗ/В4-12, А5/В5-12);
С - бензиновые и дизельные двигатели с катализатором отработавших газов (С1-12, С2-12, СЗ-12, С4-12);
Е - тяжело нагруженные дизельные двигатели (Е4-12,Е6-12,Е7-12, Е9-12).
После буквы выставляется арабская цифра, обозначающая величину рабочих характеристик - чем выше цифра, тем выше требования класса
Различают кинематическую и динамическую вязкости. Кинематическая характеризует текучесть моторного масла в условиях низких и высоких температур.
Динамическая или абсолютная вязкость — это физическая величина, характеризующая силу сопротивления, возникающую при перемещении со скоростью 1 см/с двух слоев жидкости площадью в 1 квадратный сантиметр, находящихся на расстоянии в 1
сантиметр друг от друга, которую рассчитывают по формуле :
П = v-p-10-1 (1)
где р - плотность при той же температуре, при которой определялась кинематическая вязкость, кг/мЗ , V - кинематическая вязкость, м2/с.
Кинематическую вязкость определяют по общепринятым стандартам температурах и для каждого класса масел имеются определенные значения в сантистоксах (cST или сСт).
Если речь идет о парафиновых маслах, то кинематическая вязкость больше на 16 - 22%,
при разных (таблица 1)
а у нафтеновых масел эта разница меньше - от 9 до 15%. При этом, вязкость по стандарту SAE имеет цифровое или цифробуквенное обозначение, из которого можно определить сезонность смазочного материала и температуру окружающей среды, при которой его можно использовать.
Большинство производителей моторных масел защищают свою продукцию от подделки, расфасовывая их в фирменные канистры с пломбированием горлышек канистр. Канистры, пробки и наклейки периодически изменяют, применяют голографические знаки и т.д. Однако, несмотря на это, объемы контрафактной продукции не уменьшаются. Так в настоящее время третья часть всех моторных масел, которые продаются в России, являются поддельными. С прочей автохимией еще хуже: почти половина охлаждающих и тормозных жидкостей являются контрафактом -Ш-]
Таблица 1 - Сравнительные данные отечественной и зарубежной классификации
Классы рнэсскпц, Вязкость Пределы вязкости Соответствующий
принятый в России т* -1ГС, (Ст. при 1<ЮХ, cCf КЛАСС
(ГОСТ) ЯЙ бвПв* Н§ НЯН99 не белее noSAE
--- --- зд --- 0W
% 1250 зд --- 5W
2600 4.1 --- 10 w
_ ^ 6000 5.6 --- 15 W
(Ь 10400 55 --- 20 W
в --- 6.6 7,0 20
в --- 75 95 20
10 - 9.5 И5 зо
12 --- 11,5 13.0 30
14 -■ 13,0 15,0 40
16 --- 15,0 15.0 40
20 18.0 23,0 50
зуе 1250 7.0 9.5 SWfflQ
ys 2600 5.6 7.0 1QVW2Q
чуа 2600 7.0 9.5 10VW20
4/10 2600 95 115 10W/30
VI0 6000 95 (15 10W/30
Via 6000 115 13.0 10W/30
V« 6000 13.0 15.0 15W40
УЮ 10400 95 115 20W/30
У14 10400 13,0 16,0 20W/30
У16 10400 15,0 13.0 20W/30
Существующие методы и средства определения параметров моторных масел
К методам оценки физико-химических свойств моторных масел относятся способы и средства определения вязкостных характеристик, щелочности, зольности, температуры вспышки и застывания, а также определение степени чистоты и содержания в них механических примесей и воды. Кроме того, для базового масла (до введения в него присадок) определяют коксуемость и цвет. Все перечисленные методы испытаний стандартизованы и входят в соответствующие стандарты на масла, требуют
специализированных лабораторных приборов и длительного времени . В связи с чем, была поставлена задача - разработать экспресс-метод для сравнительного анализа моторных масел по их эталонам, поэтому «моторные методы» не обсуждаются
и не используются
В соответствии с определениями стандарта под кинематической вязкостью понимается сопротивление жидкости течению под действием гравитации, а под динамической - отношение сдвига к скорости сдвига жидкости .
Поэтому сущность метода определения кинематической вязкости заключается в
измерении калиброванным стеклянным вискозиметром времени истечения, в секундах, определенного объема испытуемой жидкости под влиянием силы тяжести при постоянной температуре, т.е. кинематическая вязкость моторного масла является произведением измеренного времени его истечения на постоянную вискозиметра
(некоторые из рекомендуемых стандартом вискозиметров, приведены на рис.2).
Определение вязкости выполняют для 2-х образцов жидкости (п. 9.2.1). При этом, как следует из примечаний (п.9.1.2 и к таблице А.1 обязательного приложения «А») к
указанному стандарту ^ - «Наполненный вискозиметр выдерживают в бане до тех пор, пока он не прогреется до температуры испытания... Обычно достаточно 30 мин, кроме определений очень высоких значений кинематических вязкостей. . рассчитывают среднее арифметическое значение двух измерений времени истечения. Если же два измерения не согласуются, следует повторить определение после тщательной очистки и сушки вискозиметра и фильтрации образца... Каждый диапазон кинематической вязкости требует ряда вискозиметров ... Время истечения вискозиметров должно быть не менее 200 с ... Постоянные вискозиметра проверяют ...с помощью эталонных масел».
а) - вискозиметр б) - вискозиметр в) - вискозиметр типа Канон-Фенске типа Пинкевича типа Убеллоде
Рисунок 2 - Вискозиметры
Так как, при движении жидкости под действием силы тяжести, давление жидкости пропорционально ее плотности, то динамическую вязкость вычисляют, как произведение кинематической вязкости на плотность .
На основании п.4.1. ГОСТ 8.595-2004-^- «Для измерений массы продукта___ применяют
(Изм. редакция, Изм. № 1):
- прямой метод динамических измерений;
- косвенный метод динамических измерений.
Для измерений массы продукта в мерах вместимости и мерах полной вместимости применяют:
- прямой метод статических измерений;
- косвенный метод статических измерений;
- косвенный метод, основанный на гидростатическом принципе.»
При этом, «4.5. При косвенном методе статических измерений массу продукта определяют по результатам измерений:
а) в мерах вместимости:
- уровня продукта - стационарным уровнемером или другими средствами измерений уровня жидкости;
- плотности продукта - переносным или стационарным средством измерений плотности или ареометром по ГОСТ 3900, ГОСТ Р 51069 или лабораторным плотномером в объединенной пробе, составленной из точечных проб, отобранных по ГОСТ 2517;
- температуры продукта - термометром в точечных пробах или с помощью переносного или стационарного преобразователя температуры;
- объема продукта - по градуировочной таблице меры вместимости с использованием результата измерений уровня продукта
б) в мерах полной вместимости:
- плотности продукта - переносным средством измерений плотности или ареометром в лаборатории по ГОСТ 3900, ГОСТ Р 51069 или лабораторным плотномером в точечной пробе продукта, отобранной по ГОСТ 2517;
- температуры продукта - переносным преобразователем температуры или термометром в точечной пробе продукта, отобранной по ГОСТ 2517;
- объема продукта, принятого равным действительной вместимости меры, значение которой нанесено на маркировочную табличку ...
Результаты измерений плотности и объема продукта приводят к стандартным условиям по температуре 15 °С или 20 °С, или результат измерений плотности продукта приводят к условиям измерений его объема в мерах вместимости и мерах полной вместимости.
Коэффициент объемного расширения продукта определяют в соответствии с МИ 2632
ТИ *Ж4
АНТ-2
а) 6)
Рисунок 3 - Ареометр («а») и пикнометр («б») для измерения плотности
Как следует из вышеизложенного, важнейшими свойствами масел, которые характеризуют их эксплуатационные свойства, являются изменения вязкости масел в зависимости от температуры, которые, как правило, определяются или отношением
вязкости при двух крайних температурах vmhh/vmako или п0 индексу вязкости
Расчет индекса вязкости выполняется по ГОСТ 25371 — 97, который определяет индекс вязкости (VI ), как расчетную величину, «которая характеризует изменение вязкости нефтепродуктов в зависимости от температуры» . Если индекс вязкости находится в пределах от 0 до 100, то его рассчитывают, как отношение вязкостей, определяемых при 40 °С и 100 °С по формулам
L-U L-U
Ш=------ 100: VI =---- 100,
L-H D
где U — кинематическая вязкость таблице ГОСТа, опираясь на величи
Если индекс вязкости более 100; логарифмов и таблицы стандарта (те
VI = {[(anti log N )-1]/0,0075} +100;
где U и Y - кинематические вязкости при 40 и 100°С соответственно для испытуемых нефтепродуктов;
Н = 0,1684 Y 2 +11,85 Y - 97 - полином вычисления кинематической вязкости.
Таблица 2 - Измеренные значения L , D , Н для кинематической вязкости
(2)
масла при 40 °С; значения L , Н и D находят по ну кинематической вязкости масла при 100 °С.
, то его находят по формулам с использованием зблице 2) по формулам:
N = (log Н - log U )/log Y (3)
Ккисмат- Кшсмаги- Кииематн»
'ІССХНЛ чсСкая •КСкМ
ИЗКОСТЬ L D =<£ -И) Н asıкость L D =(L - Я) И ВЭ?КрГП> L D ={L ■ Н) И
при 1О0*С. при 100 °С. при 10О
шґіс ЫкЛс *С, шАс
1 2 3 4 5 6 7 8 9 10 11 12
2 7,994 1.6 6.394 5,7 52.64 17X6 35X9 9,4 132.8 57.1 75,73
2. і 8.64 1,746 6.894 5.8 54.42 18.16 36X6 9,5 I33J 58.36 76.91
2Л 9X09 1.898 7,41 5,9 56,2 18,97 37X3 9,6 137,7 59,6 78,08
гл ю 2,056 7.944 6 57.97 19.78 38.19 9,7 140.1 60.87 79X7
2.4 10.71 2X19 8.496 6,1 59,74 20,57 39.17 9,8 N2.7 62X2 £0,46
2.5 11,45 2X9 0,063 6.2 61,52 21,28 40.15 9,9 I45X 63.54 81,67
2.6 12.21 2,567 9.647 м 63X2 22,19 41,13 10 147.7 64.86 82X7
2.7 13 7,748 ЮД5 М 65,18 23,03 42,14 10.1 150.3 66,22 84,08
2Л 13,8 2,937 10,87 6Д 67.12 23,94 43.1В 10,2 152.9 67.56 S5J
2$ 14,63 3,132 ИД 6.6 69,16 24,92 44Д4 10,3 155.4 68,9 86,51
3 15.49 3X34 12,15 6,7 71X9 25.96 45J3 10/ 158 70X5 87,72
XI E6J6 3,54 12.82 6.8 73,48 27,04 46.44 10,5 160.6 71,63 £8.95
ЗД Е7Х6 3,753 13,5 lı 6,9 75,72 28,21 47X1 10,6 I63X 71 90.39
W ES, IS 3,97 i 14,21 7 78 29,43 48 X 7 10,7 165.8 74,42 »1.4
3.4 19,12 4.196 14.03 7,1 80X5 30,63 49.61 I0.B І68Д 75,86 92.65
зх 20,09 4,428 15.66 7Х 82,39 31,7 50,69 10,9 I7IX 77X3 93,92
3.6 23.0S 4,665 16.42 7J 84X3 3X74 51,78 И 173.9 78.75 95,19
3.7 22.09 4.909 17.19 7,4 86.66 33,79 52.8В 11,1 176,6 BOX 96,45
І» 23,13 5,157 17.97 V 88.85 34,87 53.98 11,2 179.4 81,65 97,71
3.9 24.19 5,415 18,77 7,6 91,04 35,94 55,09 11,3 182,1 8JJ3 98,97
4 25Д2 5,756 19,56 7,7 93,2 37,01 ÎM IM 184,9 84,63 100Д
4,1 26,5 6,129 20,37 7,8 95.43 38.12 57 J1 11,5 187.6 86.1 1013
4 2 27.75 6,546 21.21 7.9 97,72 39.27 58.45 Iİ.6 190.4 87.61 102.8
4J 29,07 7,017 22,05 8 100 40.4 39,6 11,7 l»3J »,İS 104,1
4.4 30.48 7.56 22.92 8.1 102,3 41J7 60.74 11,8 I96X 90.75 105.4
4.Î 33.96 8,156 23,81 8Х 104.6 42,72 61.89 11,9 199 92.3 106,7
33,52 8.806 24,7 L 8,3 106,9 43,85 63.05 12 201.9 93,87 108
4.7 35.13 9,499 25.63 М тл 45,01 64.1В 12,1 204.8 95.47 109.4
4.» 36,79 10,22 26.57 8,5 111.5 46,19 65J2 I2X 207.8 97,07 110,7
4.9 İM 10,97 27,53 8,6 И 3,9 47,4 66.4В 12,3 210,7 98,66 m
5 40ДЗ 11,74 ît АЧ 8.7 пкі 67,64 IM 213,6 ТООД IlJ.Î
S.1 41.99 12,53 29/6 8,8 USX 49,75 68.79 12,5 216.6 101,9 114.7
5 2 43,76 1X32 30,43 8,9 120,9 50,96 69.94 12,6 219,6 10X6 116
SJ 45,53 14,33 зм 9 І2ЭД 522 71.1 12,7 222.6 105,3 117,4
5,4 47ДІ 14,W 32*37 91 125,7 НА 72X7 12,8 225,7 107 118.7
5Л 49,09 15,75 33.34 9Х 128 54,61 73.42 12.9 228.8 108,7 120.1
İA J№ 16M 34.32 « им 55,84 74X7 13 ası.* 110,4 121,5
продолжение Таблицы 2
1 2 3 4 5 6 7 8 9 10 11 12
13,3 235 192,1 322,9 19.3 4S3 257.* шл 37,5 1538 946,8 5ЧІ.Ї
13,2 23В,1 113,8 3242 19,9 483,6 260,9 227,7 3S L575 9725 603,1
13,3 241,2 N5,6 325,6 20 493,2 263,7 229 т5 385 1613 9983 615
13 А 244J 117,3 327 20,2 501,5 26В.5 233 39 1*53 1024 627,1
13,5 247,4 119 12В,4 20,4 510,8 274,4 236,4 39,5 1693 1052 639Д
13,6 250,6 120,8 329.8 20,6 519,9 279.8 240,1 40 1730 1079 651,S
13,7 253,8 122,6 131,2 20,8 528,8 285,3 243,5 405 1770 1106 664.2
13,8 257 124,4 332,6 21 538,4 291,3 247,1 41 1810 L133 676,6
13,9 260, t 126,2 134 21,2 547,5 296,« 250,7 41,5 1851 1162 6В9,І
14 263,3 E 28 135,4 21,4 556,7 302,6 254,2 42 1892 1191 701,9
14,1 266,6 129,8 136,8 21,6 566,4 308.6 257,8 42,5 1935 1220 714.9
К2 269,8 131,6 138,2 21,8 575,6 314,1 261,5 43 1978 1250 728,2
14,3 273 133,5 339.6 22 585,2 320 2 264,9 43,5 2021 1280 741,3
14,4 276,3 135,3 141 22,2 595 326,4 268,6 44 2064 1310 754,4
14,5 279,6 137,2 142,4 22,4 604,3 332 272.3 44,3 2108 1340 767,6
14,6 283 139,1 143.9 22,6 614,2 338,4 275.8 43 2152 1371 780,9
14,7 286,4 141,1 345,3 22,8 624,1 344,5 279,6 455 2197 1403 794,5
14,8 289.7 142,9 146,8 23 633,6 350,3 283.3 46 2243 1434 Ш2
14,9 293 144,8 148,2 23,2 356,6 286.8 465 2288 1466 821,9
15 296.5 146,8 149,7 23,4 653,8 363.3 290.5 47 2333 1498 835.5
15,1 300 N8,8 35 ІД 23,6 663,3 369 2S4.4 47,5 2380 1530 849Д
15,2 303.4 150,8 152,6 23,3 673,7 375.7 297.9 48 2426 1563 863
15,3 306.9 152,8 354,1 24 683,9 382,1 301.8 485 2473 1596 876,9
15,4 310,3 154,8 155.6 24,2 694,5 388,9 305,6 49 2523 1630 890.9
15,5 313,9 156,9 157 24,4 704,2 394.8 309,4 495 2570 1665 9053
15,6 317,5 158,9 358.6 24,6 714,9 401,9 313 50 2618 1699 919,6
15,7 321,E 161 160,1 24,8 725,7 408.8 3L7 50.5 2667 1733 933,6
15,3 324,6 E63 161.6 25 736,5 415.6 320.9 51 2717 1769 94S.2
15,9 325,3 165,2 363,1 25,2 747,2 422,4 324,9 51,5 2767 1804 962,9
16 331,9 167,3 364,6 25,4 758,2 429,5 328,8 52 2S17 IS39 977,5
16,1 335,5 169,4 366,1 25,6 769,3 436,6 332,7 52,5 2867 1875 992,1
16,2 339,2 171,5 167,7 25,8 779,7 443 336,7 53 2918 1911 1007
16,3 342,9 173,7 169,2 26 790,4 .149,8 340,5 535 2969 1947 1021
16,4 346,6 175,8 370,7 26,2 801,6 457,2 344,4 54 3020 1984 1036
16,5 зад 178,1 І72Д 26,4 812,8 464,4 348,4 54,5 3073 2022 105!
16,6 354Л 180,3 373,8 26,6 824,1 471.« 352.3 55 3126 2060 1066
16,7 358 182,5 375,4 26,8 835,5 479,1 356.4 55,5 3180 209« 1082
16,8 361.7 184,7 177 27 847 486,6 360,5 56 3233 2136 1097
16,9 365.6 187 378,6 27,2 857,5 492,9 364.6 56,5 3286 2174 1112
17 369.4 189,2 Ш2 27,4 869 500,6 368.3 57 3340 2213 1127
17,1 373.3 191,5 mj 27,6 880,6 508,3 372.3 57 5 3396 2253 1143
17,2 377,1 193,8 3833 27,8 892,3 515,9 376.4 58 3452 2293 1159
17,3 381 196,1 184.9 28 904,1 523,5 380,6 585 3507 2332 1175
17,4 384.9 198,4 186,5 23,2 915,8 53! Д 384.6 59 3563 2372 1190
17,5 388.9 200,8 188,1 28,4 927,6 538,8 388.8 w 3619 2413 1206
17,6 392,7 203 389,7 28,6 938,6 545,7 393 60 3676 2454 1222
17,7 3W,7 205,3 !9l,3 28,3 951,2 5543 i4S.fl 60,5 37Я 2496 123*
17,3 400,7 207,7 192.9 29 963,4 562J 401.1 61 3792 2538 1254
17,9 404,6 210 ІВД.6 29,2 975,4 570,1 405j 615 ЗВ50 2579 1270
13 40В,6 212,4 196,2 29,4 957,1 577.6 409.5 62 3903 2621 1286
18,1 412,6 214,8 397.8 29,6 998,9 585,3 413,5 625 3966 2664 1303
18,2 416,7 217,3 199,4 29,8 30Lİ 593.4 417,6 63 4026 27Q7 1319
18,3 420,7 219,7 201 зо 3023 601.6 421,7 63,5 40S7 2751 1336
18,4 424,9 222,2 202,6 30,5 1055 622,3 432,4 64 4147 2795 1352
18,5 429 224,7 204,3 ЗІ ЗОЙ6 643,2 443,2 645 4207 285« 1369
18,6 433,2 227,2 205,9 31,5 1119 664.5 454 65 4268 2382 1386
18,7 437,3 229,7 207.6 32 9151 686 464,9 655 4329 2927 N02
їм 441,5 232,3 209,3 32,5 1184 70« 475,9 ¥> 4392 2973 1419
18,9 445,7 234,7 211 33 1217 730Д 487 66,3 4455 301« М36
19 449.9 237,3 213,7 33,5 9251 752,8 498,1 67 4$І7 3064 N54
19,1 4542 239,8 214,4 34 32Вб 776,8 509.6 675 4580 зио N7!
19,2 458.4 242,3 2I6.J 34,5 1321 799,9 521.1 68 4645 3157 MSB
І 9,3 462.7 245 2I7J 35 1356 823,4 532,5 685 4709 3204 1506
19,4 467 247,6 219.4 35,5 1391 847Д -544 69 4773 3250 1523
19,5 471.3 250,2 221.7 36 9427 87ІД 555,6 69,5 4839 329« 1541
19,6 475.7 252,9 222,8 36,5 1464 896.5 567,1 70 4905 3346 1558
19,7 479,7 255,2 224.5 37 9509 921,8 579.3
Одним из основных свойств моторных масел являются их моющие свойства, которые характеризуют способность масел предотвращать образование высокотемпературных
углеродистых отложений на поршне двигателя. Поэтому, в соответствии со стандартом именно моющие свойства положены в основу их разделения по эксплуатационным группам. Моющие свойства моторного масла оцениваются непосредственно чистотой двигателя и косвенно щелочным числом , которое выражается через количество гидроокиси калия в мг, эквивалентное количеству всех щелочных компонентов, находящихся в 1 г масла , представляющие собой моющие присадки:
алкилсалицилаты, сульфонаты и алкилфеноляты -Ш11.
Общепринятым способом определения щелочного числа считается потенциометрическое титрование пробы спиртовым раствором соляной кислоты (установка и графики приведены на рис.З) J-Ш.
Как следует из стандарта - время определения щелочного числа одного типа моторного масла составляет от одного часа до двух.
а) - установка б) - графики титрования
Рисунок 3 - Установка и графики для определения щелочного числа
Остальные характеристики стандарта ■■-, определение каждого из которых, составляет несколько часов, не являются показателями назначения и относятся:
к показателям сохраняемости - гигроскопичность, температура помутнения, растворяемость с хладагентами, срок сохраняемости, гидролитическая стабильность, стабильность вязкости, цвет;
к эргономическим показателями - ПДК, класс токсичности;
к показателям безопасности - температура вспышки, температура самовоспламенения.
Таким образом, стандартные методы и средства не могут обеспечить экспресс-контроль моторного масла, т.к. каждый из методов требует по нескольку часов лабораторной работы.
Теоретические основы экспресс-метода весовой электрометрии моторного масла
Головной завод фирмы Castrol занимает площадь двух футбольных полей. Здесь же расположен научный центр по разработке масел. Тысячи километров труб, сотни хранилищ, уходящие в небо башни ректификационных колонн, свой порт и внутренний автотранспорт. Вместо полноценного склада готовой продукции в Гамбурге - лишь небольшой буферный перевалочный пункт (рис.4). Но будучи в 200 раз меньше по площади, чем французский Mobil, немецкий Castrol производит масла всего в 2 раза меньше - более 200 миллионов литров в год -ІШ1.
Рисунок 4 - Головной завод по производству масел фирмы Castrol
На головном заводе Castrol не занимается переработкой нефти - здесь только смешивают масла. Обходиться без значительных запасов помогает соседство с гамбургским портом - завод постоянно пополняет небольшие, на один-два дня работы, объемы ингредиентов и оперативно отгружает готовую продукцию.
Рисунок 5 - Автоматическая линия смешивания и розлива масла Castrol
Автоматизация - полная: несколько операторов с помощью компьютеров управляют смешением масел. Людей нет ни на одной из 9-ти линий розлива масла (рис.5).
Разрабатывают эти масла в трехэтажном здании, где работают 60 специалистов. В десятке лабораторий сосредоточены сотни лабораторных средств: химическая посуда, вискозиметры, спектрометры и хроматографы для анализа состава. Масла проверяют на пенообразование, умение защищать от коррозии и собственную агрессивность к металлам и резиновым уплотнениям
Линия автоматизированного розлива масел RLFMS, которая, например, используется ЛУКОЙЛом, и включена в Государственный реестр средств измерений (per. № 41083-09), имеет 11 дозаторов на 3 диапазона объема и массы дозирования канистр и
обеспечивает следующую точность :
для тары объемом 1 литр (масса от 0,82 до 0,91 кг.) - точность от 2,87 г до 3,19 г.; для тары объемом 4 литра (масса от 3,36 до 3,63 кг.)- точность от 11,76 г до 12,71 г.; для тары объемом 5 литров (масса от 4,20 до 4,54 кг.) - точность от 14,7 г до 15,89 г.
Стандарт на полимерную тару-’ требует, чтобы значения и предельные отклонения от номинальных размеров и массы тары не превышали (приложение К):
п. 9.2.1 - 0,1 мм. для геометрических размеров;
п. 9.3.1 - 0,05 мм. для толщины стенки;
п. 9.4.1 - 10% для номинальной вместимости;
п. 9.5.1 - 10% для массы тары.
Если предположить, что производители отечественной тары, укладываются в указанные допуски, то получим следующие разбросы по массе:
для тары объемом 1 литр (масса 0,07 кг.) - точность 7,0 г;
для тары объемом 4 литра (масса 0,26 кг.) - точность 26 г.
Следовательно, зная массы пустых канистр объемом 1 и 4 литра, можно взвешиванием, не вскрывая пробки и без отбора пробы, определить плотность жидкости в нём по
формуле:
Р І = (Р I - Pj )/Vij (4)
где Р і - измеренный вес / -го образца в фасованной таре; Р j - вес j -той эталонной тары; Vij = 1л, 2л...Ыл - эталонный объем і -той жидкости, заливаемой в j - объем тары.
То есть, измеряя вес образца (Р / ) на электронных весах в фирменной таре, имеющей эталонный вес (Р j ), можно вычислить плотность образца при эталонном объеме заливки в тару (Vij - 0,5 л, 1л, ...14л) со следующей точностью:
для объема 1 литр - Д% = [(820 + 2,87 + 7,0)/1000 - (820-2,87-7,0)/1000]/0,820 =
(0,01974/0,820) = 0,02407*100 = 2,4%;
для объема 4 литра - Д% = [(3360 + 11,76+26,0)/1000 - (3360-11,76-26,0)/1000]/3,360 = (0,07552/3,360)*100 = 2,24%.
Зарубежные производители, Castrol в частности, используют оборудование динамического весового дозирования CRANDALL International, которое обеспечивает
точность (при розливе в тару до 20 кг) -0,05% L2'13\ т.е. в 7 раз точнее, чем оборудование ЛУКОЙЛА, а предельные отклонения от номинальных размеров и массы тары - не более 1%. В этом случае точность определения плотности составит:
для объема 1 литр - Д% = [(820 + 0,41 + 0,7)/1000 - (820-0,41-0,7)/1000]/0,820 =
(0,00222/0,820) = 0,0027*100 = 0,27%;
для объема 4 литра - Д% = [(3360+1,68 + 2,6)/1000 - (3360-1,68-2,6)/1000]/3,360 = (0,00856/3,360)*100 = 0,25%.
Таким образом, предлагаемый в экспресс-методе способ определения массы и плотности продукции в канистре без отбора пробы (без открытия крышки и нарушения защиты канистры), сокращая время их определения на несколько порядков, обеспечивает такую же погрешность, как стандартные средства
При этом значения плотностей при различных температурах вычисляются по формуле Д.И. Менделеева !
р (t)*» Р 20°С - At ■(t - 20°С) (5)
где At - (18,310 - 13,233' р 20°С )'10-4 - температурная поправка к плотности на один градус; t - искомая температура, °С
Экспрессный метод определения кинематической вязкости нефтепродуктов через
их плотность
В России запатентован способ экспресс-определения кинематической вязкости авиационных керосинов и дизельных топлив путем регистрации относительной плотности
топлив при 20°С (Р420) в г/см3, и последующего расчета кинематической вязкости (V20) при 20°С (рис.6), по следующей зависимости J :
v20 - k1'p420 - k2, (6)
где v2o - кинематическая вязкость при 20°С, мм2/с; k^ =23,1, k2 =16,77 - для топлив с относительной плотностью от 0,780 до 0,820 г/см3; k^ =200, k2 =161,83 - для топлив с относительной плотностью от 0,820 до 0,842 г/см3.
о
В
7
6 5 4 3 г
і
о
0,78 0,79 0,80 О.Ш 0,8Й 0,83 0,84 0,85 0,86
їкотаость при 30°С {їг/см3
і . і - і ■ .і___________________________________________________________і__________і—_________t---------1---------і----------L
Рисунок б - Зависимость кинематической вязкости авиационных керосинов и
дизельных топлив от плотности по патенту РФ № 2263301.
Суть экспрессного способа заключается в том, что используется эмпирическая формула (6) зависимости кинематической вязкости топлив от их плотности (рис.б) в диапазоне от
0,780 < р420 > 0,842 г/см3^Д
При этом, когда плотность горючих определяется в полевых внелабораторных условиях в интервале температур окружающей среды от минус 10°С до плюс 30°С, ее значения не переводятся в величину при 20°С согласно ГОСТ 3900, а определяется сразу кинематическая вязкость горючих при 20°С по впервые разработанной номограмме (рис.7).
RU 2 263 301 С1
Lull И і ь 1 tİLİlili
Температура измерения плотности? °С { ~t }
Схама пальнемиид
Ь £ » w.
„ _ = р »t к S, . ^
® Л Л, 1,
Вікнсоеть нипвіл&тичєокия при 20^1 мм''/iî
О О О О ср О ifi о ю о ІҐІ
Щ t'- IÛ іО Ч- (л 14 СУІ м м о
Рисунок 7- Номограмма для экспрессного определения кинематической вязкости топлив
по патенту РФ № 2263301
Таким образом, вместо измерения скорости истечения горючего через капиллярный
стеклянный вискозиметр определяют плотность горючего с последующим переводом ее в единицы кинематической вязкости.
Запатентованный способ прошел проверку в процессе испытаний в сравнении со стандартными методами на образцах авиа керосинов и дизельных топлив производства основных нефтеперерабатывающих заводов (НПЗ) страны за последние 5 лет, результаты которой приведены в таблицах (3 и 4). Из полученных данных следует, что относительные расхождения между результатами анализов в среднем составили 5,3%, свидетельствующие о надежности заявляемого способа J— .
Таблица 3. Результаты сравнительных испытаний
Плоти4СТВ Канели та чеосав вяткогта, нм-А
ІДарка топі а на Наименование НИ? Годовал ига ври 2|GSC, гЛаа^ ГОСТ 39 Тапия ев ны а са осо 6 Росхпждениа
ГОСТ 3*09
аа^ % пт*.
ГС-1 Мрскрекхий 18*9 0,792 1,29 1.21 0,07 5,5
ГС-1 Мрокрвский 20« 0,797 1,49 1.40 0,09 5,4
РТ ЫошрпаІ 19*9 0,790 им 1.2« 0,05 3,9
Т-1 Краги ода рока І. 19*5 0,946 1,62 2,92 0,1« 5,0
РТ Воіго град оса 1 19*5 0,79« 1,49 1.98 0,04 2.0
ГС-1 Волгоград ска ■ 19*5 0,799 М2 1.4« 0,04 2,0
ТС-1 Кв райска а 1995 0,792 М§ 1,21 0,07 5,4
ГС-1 Смекай 19*5 0,799 1,99 1.4* 0,09 5,5
ТС-1 Амана а 19*5 0,799 1,94 1.4* 0.12 9,2
РТ Напкуйбниеккнй 19*5 0,791 us 1.20 0,05 4,0
Т-6 Авпаргагай 19*5 0,924 2,99 3.20 0,30 M
РТ .Ал га рек ай 19*5 0,7*0 1,44 1.50 0,06 4,2
РТ Сигровсмв А 1995 0,791 1,92 1.20 0.12 *,1
ТС-1 Гарак.овсквй 19*4 0,799 1Д4 1.2« 0,09 6.0
ТС-1 О река й 19*5 0,799 U* 1.2« 0,02 1,9
ТС-1 Нсвауфаискай 1995 0,799 1,9* 1.2* 0.13 *.4
ТС-1 Ха Барр»™ і 19*5 0,7*5 1,56 1,68 0.10 6.0
ТС-1 Рнтапскаа 2091 0,795 1,95 1.92 0,09 2.2
Всего 1? обратное Виго 15 НПЗ дим 0.07 5,3
Межлабораторные образцах топлив
испытания в семи лабораториях различных ведомств на шести для реактивных и дизельных двигателей подтвердили полученные
результаты (таблица 4) с
погрешностью в среднем
2,6% -U^1.
Таблица 4. Результаты межлабораторных испытаний
Марка тпілгм Наиер ирздаем Кв нематические вз лпкт. мн^/с РосхФждание
ГОСТ ЗЭ-20СО Эка пресс-мет ад, (инвера ла6*рот*рв1|
121 321 331 257 1.7 İİ77* E3S34 мм2А *1 ьти.
■-0.2-82 1 4.79; 4.09 4,75 4,91 4.75 5.25 4.90 4.3* 5.99 9.10 3.6
2 4,04; 4,05; 4.9« 4,69 4.90 4.90 5.20 4.« 4.29 4,8* 9,116 3.1
1-0.5-82 1 4.74; 4.09 4,75 4.92 4.73 5.20 4.91 4,47 5.İ9 «,1В 3.3
2 4.85; 4,89 4,69 4.90 4.90 5,25 4.90 4J5 - 0,18 3.7
9-0.2 (-35j) 1 2.3* 2,41 2.40 2,42 2.40 2.43 2.49 2.46 9.92 9.1
2 2v47 2,49 2.42 2.40 2.» 2.40 2.35 2.45 9,97 2.3
т-в 1 2^5 2,25 2.26 2,27 2.2* 2.31 2.29 2.34 9,U и
2 2,31 2,69 2.25 2,23 2.30! 2.20 2J5 2.3* 9,95 2.5
ТС-1 1 1.44 1,46 1,45 1.47 1.45 1.49 1.45 1.56 С, 02 2.9
2 1.45 1,39 1.40 1.43 1,40 1.43 1.49 1.45 С, 0-6 4,3
РТ 1 1.45 1,46 1.45 1.47 1.45 1.49 1.49 1.46 9.91 1,9
2 м« 1,44 1.42 1.42 1.45 1.43 1.44 1.41 9,93 3,9
Среднее 9,98 2.6
Проведенные испытания показали, что разработанный экспрессный способ позволяет надежно определять кинематическую вязкость авиационных керосинов и дизельных топлив различного происхождения, как в современных контрольных лабораториях ГСМ, так и на местах производства и применения топлив.
Технико-экономическая оценка предлагаемого экспрессного способа по сравнению со стандартными методами показала, что экспрессный способ позволяет исключить применение дорогостоящей аппаратуры, сократить продолжительность анализа
нефтепродуктов в 40 раз, а стоимость анализа в 20 раз
Экспрессный метод определения вязкости по электрическим параметрам
жидкостей
При исследовании любой электрохимической системы, ее возбуждение синусоидальным сигналом и фиксация ответа на это возмущение в состоянии равновесия, является наиболее простым методом определения её транспортных функций что позволяет
LULL;
- выполнять высокоточные измерения в большом временном (или частотном) диапазоне (104-10‘6 сек. или 10~4-10б Гц), т.к. отклик может быть длительным и устойчивым, а, следовательно, может быть усреднен по большому периоду времени;
- обрабатывать сигнал с помощью линеаризации (или другими методами).
Метод, при котором импеданс электрохимической системы измеряется как функция от частоты, называется, спектроскопией электрохимического импеданса (СЭИ, EIS), а сама система может быть описана на основе теоретического рассмотрения процессов и создания физических моделей, или путем конструирования электрической эквивалентной схемы из простейших элементов, в первую очередь, сопротивлений (резисторов) и емкостей (конденсаторов).
Для оценки параметров и получения информации о процессах транспорта носителей заряда проводится сопоставление моделей.
Импедансная спектроскопия играет все более значительную роль в фундаментальных и прикладных исследованиях. Ее можно использовать для исследования любого типа твердых и жидких материалов: ионных, смешанных, полупроводниковых и изоляторов. При этом используются соотношения проводимости в жидкостях, которые описываются законами Ома (Z, G, R, I), и определяются анализаторами импеданса и измерителями иммитанса 1 У
Помимо определения удельных параметров образцов и эталонов, может быть использовано основное уравнение характеристического импеданса :
Zİ = V(Ri + і со Li )/(G і + і со С і) (7)
где: Zi - комплексное сопротивление; щ = 2nf (циклическая частота); і - мнимая единица; Ri - сопротивление; Li - индуктивность; Gi - проводимость; Сі - емкость.
Дополнительные параметры для идентификации могут быть получены, при решении уравнения импеданса (7) и применении полученных результатов для определения критериев подобия жидкостей. Так, для определения магнитной проницаемости (р) электролитов и полупроводниковых жидкостей, имеем :
р = Z2-2İ /S-G-(l + i)2-p0-O) (8)
где і - мнимая единица; рО - магнитная постоянная; f - расстояние между измерительными электродами; S- площадь поверхности электрода; Z - значение комплексного сопротивления; G - значение проводимости; w - частота измерения.
Определение магнитной проницаемости жидкостей-диэлектриков осуществляется по формуле:
р = Z2-(l- i tgö) (9)
Циклическое измерение указанных выше параметров, дает возможность вычисления следующих критериев гомохронности (критериев подобия):
Hol = cot; (электродинамического), (10)
Но2 = р 12/ о t; (электромагнитного), (11)
НоЗ = e о / t; (диэлектрического), (12)
Но4 = L/R t; (электроиндуктивного), (13)
Но5 = C/G t; (электроемкостного), (14)
где t - время; L - измеренное значение индуктивности, а остальные обозначения такие же, как в предыдущих формулах.
Различают два основных вида зависимости электропроводности жидких сред от концентрации компонентов в них: удельную электропроводность - □ и эквивалентную - Л,
определением которых «занимается» кондуктометрия Удельная электропроводность является обратной величиной удельного сопротивления (г ) □ = 1/ , а эквивалентная и удельная электрические проводимости связаны соотношениями
Л = И0 /С и А = кем/ (15)
где С - концентрация компонентов жидкости (моль/мл); V - объем жидкости (мл), содержащий при данной концентрации 1 грамм-моль компонента.
Частным случаем анализа импеданса является диэлькометрия J-221 в которой определяются соотношения между потерями проводимости и диэлектрическими потерями, в жидких и вязких средах с диэлектрическими свойствами в частности, которые выражаются, диэлектрической проницаемостью £ , и определяются через измерения тангенса угла потерь tg5 и проводимости о:
£ = £'-/ е" —> tg5 = £"/ Е'= 0/(0 £' єО (16)
г д е £ - относительная диэлектрическая проницаемость; е" - комплексная часть
диэлектрической проницаемости; £' - действительная часть диэлектрической
проницаемости; о - активная проводимость, учитывающая оба вида потерь; (о - круговая частота; еО - диэлектрическая проницаемость вакуума (8.85-10-12 Ф/м).
Комплексная часть диэлектрической проницаемости е" стремится к нулю, как при малых, так и при больших значениях сот и достигает максимума при сот т = 1 (рис.8 «б»), в котором
в іє~£~)
=
/ (в + О
£_ =
(17)
а комплексная диэлектрическая проницаемость равна
£ = £_+-
1 + гел
£ Ч- 2 +2
Рисунок 8. Диэлектрическая проницаемость от температуры (а) и частоты (6)
После чего, используя формулы Дебая (19) и Паулса (20), определяются и усредняются времена релаксации Ї221: макроскопическое (г ) и микроскопическое (т^ ):
t = tm -(є + 2)/(е00 + 2) (19)
t = 3etm/(2£ + Єоо) (20)
Зная микроскопические времена эталонной и исследуемой жидкостей, можно использовать отношение указанных времен, которые по Дебаю связаны с динамическими вязкостями - П \, простым соотношением:
г
fil
(21)
Таким образом, определяя предлагаемым способом плотность - р, диэлектрическую проницаемость и тангенс угла потерь - tg5, и вычисляя динамическую вязкость - п через времена релаксации - t и tm (вместо определения кинематической вязкости - v с
помощью вискозиметра), вычисляем кинематическую вязкость масла по «обратной» (1) формуле стандарта (1), в т. ч. с учетом температуры и частоты
V = п/Р (22)
Инструментальные и алгоритмические особенности экспресс-контроля
Диэлектрические материалы и жидкости в большинстве случаев являются неоднородными диэлектриками, состоящими из разнородных по электрическим свойствам частиц, слоев или молекул (текстолит, пенопласты, масла и т.д.).
Для расчета и определения диэлектрических потерь и диэлектрической проницаемости таких неоднородных диэлектриков используют эквивалентные схемы замещения, как отдельных компонентов, так и материала или жидкости в целом. Возможные варианты упорядоченного расположения компонентов могут быть представлены в виде их параллельного и последовательного включения. При параллельном включении -плоскость раздела компонентов параллельна вектору напряженности электрического поля (рис. 9) и емкость такого конденсатора определится по формуле (23).
Рисунок 9 - Плоский конденсатор с 2-мя параллельно включенными диэлектриками
С = с1+са
fI
+
«о ' «а ’ ^
fj г-2
h
(23)
где S = S 1+S 2 - площадь электродов конденсатора (м2); є 1 и є 2 - диэлектрические проницаемости компонентов; С 1 и С 2 - емкости компонентов (Ф); є - эффективная диэлектрическая проницаемость; є 0 = 8,85'10-12 - электрическая постоянная (Ф/м); h -толщина диэлектрика (м).
При одинаковой толщине диэлектриков их объемные концентрации 0 пропорциональны их площадям и составляют соответственно для первого 0 1 и второго 0 2 диэлектрика:
01 =
Sı
sl+s3
Sı
s1+s2
91+92=1 (24)
Подставив уравнение (24) в (23), получим:
go gı &ı '(^i Sji @2 (gi + £ (*^i + ^а)
й h k
(25)
откуда - £ = £ !'0 ! + £ 2'0 2 (26)
Если неоднородный диэлектрик состоит не из двух, а из нескольких компонентов, то
уравнение (26) имеет вид:
2$=i
!-Я
1-І
(27)
Диэлектрические потери в таком неоднородном диэлектрике складываются из потерь в компонентах:
Р = Р і + Р 2 = U 2'ы 'С 'tgS (28)
где Р - активная мощность, рассеиваемая во всем диэлектрике (Вт); Р 1 и Р 2 - активная мощность, рассеиваемая в компонентах (Вт); ср - круговая частота (рад/с); С - емкость конденсатора (Ф); tg5 - тангенс угла диэлектрических потерь неоднородного диэлектрика; U - напряжение на электродах (В).
При параллельном включении диэлектриков (рис.9) к каждому компоненту приложено одинаковое напряжение, поэтому
Р \ — U 2 -о) 'С ^'tgö і и Р i - U 2 "со 'С 2'tg5 2 (29)
где tg£ 1 и tg£ 2 ' тангенсы углов диэлектрических потерь компонентов.
Из уравнений (28) и (29) следует, что тангенс угла диэлектрических потерь двухкомпонентного диэлектрика при параллельном соединении равен
tgâ -
с, tgâx + Cj tgöş
с
(ЗО)
Полагая, что
q _ ~ ^ і _ д g Fı ^ _ ga ^ ga
1 ~ h h 2 ~ h 3 h (31)
уравнение (30) можно привести к виду 9Х ■ £) £göl + $7 £, ■ tgâ^
(32)
tgâ = ■
Если неоднородный диэлектрик состоит из нескольких параллельно соединенных слоев, тогда уравнение (32) может быть записано в виде
---------
1 (33)
При последовательном соединении компонентов плоскость раздела их также перпендикулярна вектору напряженности электрического поля (рис.10).
Рисунок 10 - Плоский конденсатор с 2-мя последовательными диэлектриками
Однако объемные концентрации компонентов будут пропорциональны их толщинам:
Ö, = ■
*1
Aj +h2
А, + Aj
0 ! + 0 2 = 1 (34)
gı -g3 1 _ 3 ( 8%
9X - £\±92 є Fj ^25)
Формулы (26) и (35) являются частными случаями формулы Лихтенеккера: £■* - £х -$х + ^26)
гдех - константа, характеризующая пространственное расположение компонентов и принимающая значение от х = +1 (для параллельного включения компонентов) до х = -1 (последовательное расположение компонентов).
Формула Лихтенеккера (36) широко применяется для расчета диэлектрической проницаемости мелкодисперсных смесей. Она дает результаты расчета, достаточно хорошо совпадающие с измеренными величинами, если є ^ и є 2'0 2 не очень сильно отличаются друг от друга.
Практическое применение находят мелкодисперсные хаотические смеси (пенопласты, керамика и др.), имеющие предельно разупорядоченное строение и жидкие диэлектрики. Для них можно принять, что константа х , в формуле (36) стремится к нулю. Чтобы избавиться от неопределенности, продифференцируем это уравнение по х :
= £■*. ]û є = в1 Єу ю #i +в2 ■ İn f3
(37)
Примем X =0 и тогда получим, так называемый,логарифмический закон смешения:
1л Є — ^ ■ In 4"^2 "In Є2 (38)
Если неоднородный диэлектрик содержит более чем две компоненты, то формулы (37) и
(38) приводятся к виду:
1-й i-й Ж
23=1 = 23 !п Щ
і-l і-l 1
(39)
Для пенопластов, пенокерамики и других пористых материалов, состоящих из твердого и газообразного диэлектрика, а также для жидких диэлектриков, удобнее пользоваться не объемными концентрациями компонентов, а их массовым содержанием в смеси. Плотность смеси можно рассчитать на основании арифметического закона смешения:
р = 0 1'Р 1+ 0 1 'Р 2 (40)
где р - ПЛОТНОСТЬ смеси (КГ/мЗ); Р і И р 2 - плотности компонентов (кг/мЗ).
Тогда из формул (38) и (40) следует:
ву-~ Ъе = — ю «i
А А Т
(41)
гдер і - плотность монолитного (сплошного) диэлектрика; р - плотность смеси (пенопласта и т.д.).
Диэлектрические потери в многослойном неоднородном диэлектрике складываются из потерь в каждом из последовательно соединенных компонентов диэлектрика согласно уравнению (28):
Р = Ру + Р2 = U2 ■ Ф- С ■ tgâ = Uy Ф-Су- tgду л~и\-Ф С2 tgâ2 (42)
где U і и U 2 - падение напряжения на конденсаторах С і и С і - соответственно (рис.10).
Учитывая, что
Г_С1 с, и, _ с щ с_
Cj + Са ‘ 7/ С,’ ^ Ч
(43)
уравнение (42) можно записать в виде:
j _ 1 + C'ı'
g г +С
і з (44)
Принимая во внимание уравнения (31) и (40), в окончательном виде запишем:
5 Л ' + eş - ву - ügü,
«і -З
(45)
При числе слоев более двух получим:
tgS =
(46)
У всех диэлектриков величина диэлектрической проницаемости в большей степени зависит от температуры. Согласно определению, температурный диэлектрической проницаемости равен:
или меньшей коэффициент
ТКе =
1 de
~ё'Ит
град-1 (47)
где £ - относительная диэлектрическая проницаемость; Т - температура (град).
ТК є так же, как и є смесей, зависит не только от свойств и объемной концентрации компонентов смеси, но и от их взаимного пространственного расположения в электрическом поле. Поэтому формулу для расчета ТК є можно получить дифференцированием уравнения (38) по температуре:
1. — = в—.^ + в3- — -^1 = в1- ТКе-у + 9. e dT 1 e, dT 3 en dT 1 1 '
a TKe2
(48)
t.e., для мелкодисперсной хаотической смеси справедлив арифметический закон смешения.
При параллельном включении компонентов дифференцируем по температуре уравнение (26):
de ІІ Ey
de
1 d Ey
1 den
-9-----+Sn----- = e, ■ 9,-- + En ■ ft,
dT ^ dT 3 dT 1 1 Ey dT 3 ^ £a dT
(49)
Разделив левую и правую часть на £ , и, применяя формулу (47), получим:
ТКє=Єх ТКє1 + &2 — TKs2
є є (50)
При последовательном включении компонентов после дифференцирования по температуре уравнения (35) и умножения на є смеси правой и левой частей уравнения получим:
1 ds_ 1 de1 1 ds2
1 dE Е E
ТКє=-— = —в1-ТКє2 + Є2 — ■ ТЩ
и соответственно - s eh £) (52)
Для пенопластов, других пористых материалов и жидкостей, после дифференцирования уравнения (41), получим:
ТКг = -Е- ■ ТКє1 Pi
(53)
(а) - при определении р V (1 - измерительный электрод; 2 - охранное кольцо; 3 -питающий электрод; 4 - образец), (б) - при определении р S (1 - измерительный электрод; 2 - питающий электрод; 3 - нижний свободный электрод; 4 -образец)
Рисунок 11 - Схема из 3-х электродов и плоского образца
Измерение диэлектрической проницаемости и тангенса угла диэлектрических потерь на разных частотах производят, как правило, с помощью трех электродных систем (рис.11 и рис.12).
В первом случае (рис.11) емкость плоского образца рассчитывается следующим образом:
с _є-єв-$ _Є-SgTT jd^)2 _ g. Є0л-■ (ûfı+ûfa)2
* h Ah \6h ГЧ4Л
(A+d2)2
где 4 - эффективная площадь электрода (м2); dİ- диаметр
измерительного электрода (м); d 2 - внутренний диаметр охранного кольца (м).
Это позволяет получить из уравнения (54) значение относительной диэлектрической проницаемости плоского образца:
є= 0.144 — С, d
* (55)
Емкость трубчатого (коаксиального) или жидкого образца (в пФ) (рис.12) определяется по формуле:
(а) - трубчатый, (б) - стержневой, (1 - образец, 2 - измерительный электрод, 3 -
охранное кольцо)
Рисунок 12 -Система из 3-х электродов коаксиального или стержневого образца
На частотах 1000 Гц и выше охранные электроды уже не дают требуемого эффекта, и поэтому при испытаниях диэлектриков применяют датчики без охранных электродов, но в этом случае, при измерении емкости образцов, вводят поправку на краевую емкость, которая учитывается следующим образом:
Сх = Сизм - Скр - Сэл (57)
где Сизм - измеренная емкость образца (пФ); Скр - краевая емкость (пФ); Сэл - емкость измерительного электрода относительно земли, (пФ).
Значение емкости Сэл по отношению к земле определяется путем непосредственных измерений или приближенно, для плоского круглого электрода диаметром (D ), оно может быть рассчитано по формуле:
Сэл =17,7-10-12-D (58)
При этом, относительную диэлектрическую проницаемость материала определяют, как отношение емкости Сх конденсатора, в котором пространство между электродами заполнено испытуемым диэлектриком, к емкости СО при таким же образом расположенных электродах в вакууме (воздухе):
£ = С к/С q (59)
При испытании жидкостей или образцов трубчатой формы (рис.9) толщина образца h должна быть меньше 0,1 D
Скр = 3,14 (D !+Л ) (1,9£ - 5,81дЛ - 10,6) (60)
При испытании жидких диэлектриков краевую емкость можно учесть с помощью воздуха (е = 1,0) и калибровочной жидкости, для которой точно известна величина
диэлектрической проницаемости, а величина tgö <0,01 (например, криоскопический бензол при Т =20°С имеете = 2,29, a tgö = 0,005).
В этом случае вводится понятие «постоянной электродов», которая определяется по формуле:
С0 = (Ск - Св)/( e - 1) (61)
где С к- емкость с калибровочной жидкостью; Св -емкость с воздухом; є -относительная диэлектрическая проницаемость калибровочной жидкости.
Величина С q учитывает дополнительную емкость за счет изменения распределения поля у краев электродов при замене воздушного диэлектрика жидким, в связи с чем имеем следующую уточненную формулу расчета диэлектрической проницаемости жидкости с учетом краевой емкости:
£ = (Ск - Св)/С0 + 1 (62)
Таким образом, имея эталонные данные определенной совокупности из указанных выше параметров и соответствующие приборы (электронные весы и измеритель иммитанса), можно в течение нескольких минут идентифицировать образец, на предмет соответствия (или несоответствия) эталону.
Принимая во внимание вышеизложенное, были разработаны и изготовлены два емкостных датчика: плоский 2-х электродный и коаксиальный 3-х электродный.
Высота медной пластины датчика («д»), приклеенной с внутренней стороны к стеклотекстолиту толщиной 2 мм,- 200 мм., ширина - 20 мм., а зазор между ними (для образца масла) - 2 мм. (рис.13).
Из формулы (23) следует, что емкость датчика (на воздухе) должна составить:
С д = 1,005* 8,85-10-12 *0,2*0,02/0,002 = 1,77885Е-11 Ф или 17,8 пФ.
Измерениями на частоте 1 кГц зафиксированы следующие значения: С = 26,9 пФ при tg5 = 0,0154; R = 8,972 Мом; Н = 999,1 при Q = 53,72.
При калибровке в криобензоле получены следующие значения: С = 50,1 пФ при tg5 = 0,044. Тогда, применяя формулу (63) получим краевую емкость электродов:
Со- (50,1 - 2б,9)/(2,9 -1,0) = 12,21 пФ
Таким образом, получим результирующую формулу расчета относительной диэлектрической проницаемости масла:
£ = (С изм - 2б,9)/12,21 +1 (63)
Для повышения точности измерений был изготовлен 3-х электродный коаксиальный
датчик (рис.14).
Рисунок 14 - Емкостный 3-х электродный коаксиальный датчик
Высота всех медных трубок и стержня датчика («г») - 150 мм. Диаметр внешней трубки («в») 22 мм с толщиной стенки 1 мм, а внутренней - 12 мм. («б») с толщиной стенки 1 мм., диаметр стержня («а») - 4 мм. То есть зазор между внешней трубкой и внутренней (для образца масла) - 4 мм. (рис.12), а зазор между внутренней трубкой и стержнем (также для образца масла) - 3 мм.
Из формулы (56) следует, что емкость наружной секции датчика (на воздухе и без охранного кольца) должна составить:
Сдн = 6,28*1,005* 8,85-10-12 *0,15/1п(0,20/0,12) = 1,64017Е-11 Ф или 16,4 пФ,
а емкость внутренней секции датчика (между внутренней трубкой и стержнем -
Сдв = 6,28*1,005* 8,85-10-12 *0,15/1 п(0,10/0,04) = 9Д438Е-12 Ф или 9,1 пФ.
При параллельном соединении секций получается - 25,5 пФ, а при последовательном -16,4*9,1/(16,4 + 9,1) = 5,9 пФ.
Измерениями на частоте 1 кГц зарегистрированы следующие значения:
С дн = 25,7 пФ при tgd> = 0,0168; R дн = 0,111 Мом; Н дн = 999,5 Гн при Q = 58,35;
С дв = 15,8 пФ при tgö = 0,0445; R дв = 0,279 Мом; Н дв = 999,9 Гн при Q = 39,47;
При параллельном включении секций получено - 35,1 пФ при tgS = 0,0127, а при последовательном включении - 7,8 пФ при tgc5 = 0,0104.
При калибровке в криобензоле были получены следующие значения: С дн = 47,3 пФ при tgS = 0,0111, а С дв = 26,3 пФ при tg5 = 0,0221.
При параллельном включении секций получено - 73,6 пФ при tgS = 0,0231, а при последовательном включении - 16,9 пФ при tg<5 = 0,0087.
Тогда, применяя формулу (63) вычислим краевую емкость электродов:
С одн = (47,3 - 25,7)/(2,9 -1,0) = 11,37 пФ С од в - (26,3 - 15,8)/(2,9 -1,0) = 5,53 пФ
Следовательно, результирующая формула расчета относительной диэлектрической проницаемости масла «внешним и внутренним» датчиками:
£ = (С из мдн - 25,7)/11,37 +1 (64)
£ = (С измдв - 15,8)/5,53 +1 (65)
С точки зрения точности измерений и диэлектроспектроскопии
наиболее предпочтительнее использовать измерители иммитанса
Г251
компонентов масла Е7-20ІШ или Е7-25
Измерители иммитанса Е7-20 (рис.15 «а») и Е7-25 (рис.15 «б») - прецизионные приборы класса точности 0,1 с широким диапазоном рабочих частот 25 Гц - 1 МГц и высокой скоростью измерений до 25 измерений/сек., т.к. использованы микропроцессоры. Измеряемые параметры: индуктивность (Ls, Lp), емкость (Cs, Ср), сопротивление (Rs, Rp), проводимость (Gp), фактор потерь (D), добротность (Q), модуль комплексного сопротивления Z, реактивное сопротивление (Xs), угол фазового сдвига (ф); ток утечки (I). Оба прибора внесены в Государственный реестр средств измерений республики Беларусь (Е7-20 per. № 017436, а Е7-25 per. № РБ 0316 3593 07). Однако, с точки зрения экспресс-контроля, указанные приборы имеют небольшой диапазон усреднения измеряемых параметров: 2 предела усреднения - по 10 и по 100 измерениям.
а) б)
Рисунок 15 - Внешний вид Е7-20 («а») и Е7-25 («б»)
Измеритель иммитанса Е7-22J-2Ğİ, помимо автономного режима работы, единственный из приборов этого класса имеет режим измерения среднего значения по 3000 измерений, обеспечивая измерение параметров иммитанса (сопротивления, емкости, индуктивности) с базовой погрешностью 0,7%, внесен в ГСИ РФ (per. № 24969-08). Е7-22 с
коаксиальным датчиком успешно применяется в промышленности (рис.17).
OGıutıc характеристики:
лНСИЛСВ...,............................. жядкокрнстаглнчеекнД, + 1Д разрядный
выбор длдгштонд намерения............................. автематнЧесХнА/ручнпЦ
старость чтмерсния...........-.....................................2.1 ніиУс
питание........................................элемент питания 9 В гмпа «Крона»
габаритные размеры, им...................................... 192 в 62,6 *91
масса, кг..........................................................-.0.36S
Условна хранения н эксплуатации;
температуро хранения, °С ......................................... -20.. 00
рабочая тем пература. °С...............................— ............--О.. SO
относительная влажность. %..................... нс более 8S, бет конденсации влаги;
Рисунок 16 - Внешний вид и общие характеристики Е7-22
Рисунок 17 - Измеритель уровня и диэлектрических свойств жидкости Исследование эталонов и разработка экспресс-методологии
Все канистры с эталонами были взвешены, после чего вскрыты и измерены
диэлектрические параметры масел, расфасованных в них ^1. Далее эталонные масла были перелиты в другую тару, а канистры были вымыты, высушены и взвешены, после
чего вычислены их плотности и параметры вязкости по следующим алгоритмам İ2S1.
Как следует из приведенных выше формул и уравнений, измерения с помощью датчика (рис.14) двух емкостей (Сдн и Сдв) и двух тангенсов угла потерь (tg5 1,2), позволяют вычислить среднее значение относительной диэлектрической проницаемости (£ )
образца масла, и при наличии данных эталона (тэ,г|э), вязкость образца масла (по), через вычисление его макроскопического времён релаксации (т):
определить динамическую микроскопического (то) и
Сдв 'Сдн + Сдн - Сов
По =
tgS =
2C„-C„
г,
ъж
2є + \ Т°
{£ - ) СОТ
Є + Є_й)іТІ
(бб)
где С он - емкость наружного датчика в воздухе, пФ; С ов - емкость внутреннего
датчика в воздухе, пФ.
Принимая во внимание, что измерения проводятся на фиксированных частотах (100 Гц, 1000 Гц и т.д. до 1 МГц), а значениег со —> 1 (рис.1 «а») при высоких температурах (например, при парообразовании), получим:
tg<5 ■(£ + £ со со 2т 2) = (£ -£ оо)-а;г или tgS -ш 2т 2-(£ - 1)шг + tg<5 ■£ = 0 (67) Заменяя о; на 2nf и разделив каждый член уравнения на множитель при г 2, получим:
з (£-1) S ,,
Г -т—-----— +---- = 0
2 ırflgS Аж1/1
(68)
Подставляя частоту измерения (1000 Гц) и измеренные значения tgS и £ образцов, найдем макроскопические времена релаксации (г ) исследуемых образцов масла при температуре окружающей среды (22°С):
ги =
(*-1)
Г
АШж/igS \\610*4яг/г {б9)
После этого по 3-му уравнению системы (бб) определяется микроскопическое время (то) образца, а по второму - его динамическая вязкость (q о).
Для сравнения полученных результатов, с эталонными значениями при трех стандартных температурах (15°С, 40°С, 100°С) и двух отрицательных -застывания и предельной температуры эксплуатации соответствующего типа масла, использовано ограничение значения потерь в нефтепродуктах при 100°С (tg<5 <0,02) и следующие формулы температурной зависимости tg5 и е :
tgö (T ) = tgö 20 ехр[/с- (Г - 20°С)] (70)
ТКЕ — —
' е^-тУр (Га-Г,)
, 1/град. (71)
При этом изменение плотности от температуры вычисляется по формуле Менделеева :
р -____£ш_____
Ит 1ч-Д(Т-293)
или р {t ) = р 20°С - At -(t - 20°С) (72)
где рГ и р293 - плотность нефтепродуктов соответственно при температурах Г и 293 К; Зр - коэффициент объёмного расширения; At - (18,310 - 13,233' р 20°С )'10"4 -температурная поправка к плотности на один градус; t - искомая температура, °С,
и по формулам Вальтера, выражающие зависимости кинематической вязкости от температуры:
Iglg (V + 0,8) = а + b \дТ (73)
где эмпирические коэффициенты, а и d определяются по известным парам значений v и
Т, по следующим формулам:
а = Iglg {v х + 0,8) - b \дТ і и b = lg[lg (v г + 0,8)/lg (v 2 + 0,8)]/lg(7 і /T 2) (74)
Для идентификации образцов масла, использованы данные кинематической и динамической вязкости эталонов при 3-х значениях положительных температур (15°С, 40°С, 100°С) и 2-х отрицательных (застывания и эксплуатационной, зависящей от типа масла), которые были найдены по формулам Вальтера для эталонов:
5W-30 АЗ/В4 -► lglg(v +0,8) = 7,46501- 2,885031258 IgT (75)
5W-30 А5 -► lglg(v +0,8) = 8,28018- 3,21б8б8291-/дГ (76)
5W-40 АЗ/В4 -> lglg(v +0,8) = 7,23652- 2,790862353-IgT (77)
10W-40 A3/B4 —>- lglg(v +0,8) = 7,22062- 2,779909558-/g7 (78)
5W-40 DPF ->• lglg(v +0,8) = 7,29566 - 2,81431707-/gT (79)
5W-30 AP ->■ /g/gfu +0,S; = 7,71344 -2,987622646>/gT (80)
І0И/-40 64 lglg(v +0,8) = 7,22080 -2,779978337-IgT (81)
Аналогично, для идентификации образцов масла по температуре застывания, использованы формулы Вальтера (77-83), при значении кинематической вязкости 10000 мм2/с ■1321:
lg( 10000 + 0,8) = а і - b ilgT і -> r/ = 1°" ' (82)
где 10000 - кинематическая вязкость і-го нефтепродукта при температуре застывания T і, °К; Аі и Ві - константы для і -того жидкого нефтепродукта.
Подставляя полученные значения кинематических вязкостей из формул (75-82), были найдены температуры застывания эталонов, а из формул (2,3) их индексы вязкости с погрешностью менее 1%.
Результаты экспресс-испытаний партии моторного масла
По разработанной методологии были испытаны и идентифицированы три партии масла, расфасованные в 4-х и 1-литровые канистры (3068 шт., 156 шт. и 31 шт.)
Основным доказательством наличия контрафактной продукции в этих партиях масла послужили, во-первых, факты «колебаний» массы продукта в указанных емкостях (от 3551 до 3730 грамм в 4-х литровых канистрах), что свидетельствует о «ручном дозировании» при заполнении канистр маслом, а во-вторых, и это главное несоответствие плотности, динамической вязкости, кинематической вязкости при 40°С и 100°С, индексу вязкости и температуре застывания.
Таким образом, предлагаемая методология позволяет идентифицировать соответствие или отличие масел эталонам по б-ти важнейшим параметрам, а также синтезировать модель автоматизированной системы экспресс-контроля жидкой фасованной продукции, если её производитель разместит коаксиальный емкостный датчик в крышке тары, а на своем сайте «выложит образ» продукта и тары, в которую он расфасован.
Реализация модели переносного автоматизированного комплекса экспресс-контроля для производителей, потребителей и органов надзора, включая конструкцию датчика-крышки
тары и расширение метода на другие продукты подробно описана
Заключение
Представлен анализ существующих проблем в данной предметной области, а так же описаны существующие методы и средства определения параметров моторных масел, в т.
ч. со средними временами их выполнения, которые измеряются часами, из чего делается вывод, что стандартные методы и средства не могут обеспечить экспресс-контроля моторного масла.
На основе патентного поиска, а также анализа технологий и средств автоматизации при производстве моторных масел, изложены теоретические основы экспресс-метода весовой электрометрии моторного масла.
Представлена методология и результаты исследования эталонов для сравнительного анализа с образцами, по которым осуществлена экспертиза нескольких сотен канистр с маслом.
Библиография
1. Зарубежные масла, смазки, присадки, технические жидкости: ассортимент, свойства /Резников В.Д. и др. ISBN: 978-5-89551-016-2 - М.: Техинформ" МАИ,
2005. - 385с.
2. ГОСТ 4.24-84 Масла смазочные. Номенклатура показателей /ИУС 4-94 - М.:
Стандартинформ, 1994. - 14с.
3. ГОСТ 17479.1-85 Масла моторные. Классификация и обозначение. - М.: Стандартинформ, 2006.-42с.
4. ГОСТ 33-2000 (ИСОЗЮ4-94) Нефтепродукты. Прозрачные и непрозрачные жидкости. Определение кинематической вязкости и расчет динамической вязкости. - М.: Стандартинформ, 2000. - 23с.
5. ТЕХНИЧЕСКИЙ РЕГЛАМЕНТ ТАМОЖЕННОГО СОЮЗА ТР ТС 030/2012 «О требованиях к смазочным материалам, маслам и специальным жидкостям»-http://protrts.ru/reglamenty-tr-ts/tr-ts-030-2012/ .
6. ГОСТ Р 8.595-2004 Масса нефти и нефтепродуктов. Общие требования к методикам выполнения измерений-М.: Изд-во СТАНДАРТОВ, 2005. - 14с.
7. ГОСТ 3900-85 Нефть и нефтепродукты. Методы определения плотности. - М.: ИПК Изд-во СТАНДАРТОВ, 2003. - 140с.
8. МИ 2632-2001 Государственная система обеспечения единства измерений.
Плотность нефти и нефтепродуктов и коэффициенты объемного расширения и сжимаемости. Методы и программа расчета. С.-Пб.:ВНИИМ, 2001.
9. ГОСТ 25371-97 (ИСО 2909-81) Нефтепродукты. Расчет индекса вязкости по кинематической вязкости - Минск: Межгосударственный совет по стандартизации, метрологии и сертификации, 1997.-9с.
10. Кулиев А.М. Химия и технология присадок к маслам и топливам -Л.: Химия, 1985. -312с.
11. ГОСТ 11362-96 (ИСО 6619-88) Нефтепродукты и смазочные материалы. Число нейтрализации. Метод потенциометрического титрования-Минск: Межгосударственный совет по стандартизации, метрологии и сертификации, 1996.-18с.
12. Хлебушкин И.Н. Как делают масло Castrol//Авторевю.-2014.-№ 13-http ://ww w. а utoreview.ru/archive/section/detail.php?
ELEMENT_ID = 139734&SECTION_ID = 7837
13. Keith A. Crandall, Чернышев А.В. Весовое дозирование жидкостей: Новый взгляд на природу вещей-http://www.potomac.ru/news/art_gf-new.htm
14. ГОСТ Р 52620-2006 Тара транспортная полимерная. Общие технические условия. -М.: Стандартинформ, 2008. - ббс.
15. Золотов В.А., Бартко Р.В., Кузнецов А.В. Определение эксплуатационных групп моторных масел //Труды 25 ГосНИИ МО РФ.-200б.-вып.53. - 234с.
16. Зрелов В.Н., Алаторцев Е.И., Шаталов К.В., Зрелова Л.В., Бордюговская Л.Н.
Способ экспрессного определения кинематической вязкости авиационных керосинов и дизельных топлив - патент РФ 2263301, 27.10.2005, Бюл. № 30.
17. Кондуктометрия - Химическая энциклопедия-бИр://епс-dic.com/enc_chemistry/Konduktometrija-1233/ .
18. Афанасьев А.В., Москвичев А.Н., Москвичев А.А., Односевцев А.В., Орлов И.Я. Низкочастотный комплекс импедансных измерений характеристик проводящих сред // Вестник Нижегородского университета им. Н.И. Лобачевского.-2008.-№ 2, с.60-64.
19. Шабловский Я.О. Спектральный анализ негармонического и плигармонического откликов электрохимического объекта при импедансных исследованиях //Электрохимическая энергетика.-2012.-Т.2 № 12, с.111-116.
20. Белозеров В.В., Босый С.И., Кальченко И.Е., Нестеров А.А., Прус Ю. В. О
термоэлектроакустическом методе определения характеристик пожароопасности твердых и жидких веществ и материалов //Технологии техносферной безопасности.-
2010.-№ б (34).-5 с. - http://ipb.mos.ru/ttb/
21. Белозеров В.В., Бойко В.А., Голубов А.Н., Коган В.А., Прус Ю.В. О спектрометрии
жидких и вязких сред //«Спектроскопия координационных соединений»: мат-лы XI междунар. конф.-Краснодар: Куб ГУ, 2014, с. 145-146.
22. Диэлькометрия - Химическая энциклопедия-ЬИр://епс-dic.com/enc_chemistry/Dijelkometrija-527.html
23. Гусев Ю.А. Основы диэлектрической спектроскопии/уч. пособие/-Казань: КГУ,-112с.
24. Измеритель иммитанса Е7-20. Руководство по эксплуатации /УІ1ІЯИ.411218.012 РЭ -Минск: ОАО «МНИПИ», 2004.-40с.
25. Измеритель иммитанса Е7-25. Руководство по эксплуатации /УІ1ІЯИ.411218.015 РЭ-Минск: ОАО «МНИПИ», 2007.-30с.
26. Цифровой измеритель иммитанса Е7-22. Руководство по эксплуатации - Москва:
ФГУ «Ростест-Москва, 2003. - 44с.
27. Мастепаненко М.А. Информационно-измерительная система непрерывного контроля уровня топлива в емкостях: дисс.... канд. тех. наук - Ставрополь: ФГУ СГАУ, 2014,-225с.
28. Castrol-http://www.castrol.com/ru_ru/russia/products/cars/engine-oils/castrol-magnatec.html
29. Коробейников С.М., Свириденко М.В. Определение диэлектрической проницаемости и тангенса угла диэлектрических потерь кремнийорганической жидкости ПМС-20 в широком диапазоне температур: сб. науч. тр. НГТУ. - 2011. - № 2(64), с. 135-142
30. Черникин А.В. О вязкостно-температурной зависимости Филонова - Рейнольдса //Трубопроводный транспорт.-2010.-№ 6(22), с.35-37.
31. Новиков А.А. Введение в информатику процессов первичной переработки нефти: уч. пособие / ISBN 5-98298-264-4/- Томск: Изд-во ТПТУ, 2008.-200с.
32. Шадрин С.В. Эксплуатационные материалы: метод, указания - Хабаровск: Изд-во ДВГУПС, 2014. - 35 с.
33. ГОСТ Р 51634-2000 Масла моторные автотракторные. Общие технические требования - М.: Изд. Стандартинформ. 2000. - 5 с.
34. Белозеров В.В., Кудрявцев Ю.А., Плахотников Ю.Г. ИССЛЕДОВАНИЕ ПАРТИИ ФАСОВАННОГО МОТОРНОГО МАСЛА CASTROL MAGNATEC НА ПРЕДМЕТ ОПРЕДЕЛЕНИЯ КОНТРАФАКТНОЙ ПРОДУКЦИИ /отчет о НИР № 2015/12 от 23.03.2015 (ООО
"ПОЛИЭКСПEPT")-URL: https://elibrary.ru/download/elibrary_24555474_15117928.pdf
35. Троицкий В.М., Белозеров В.В. МОДЕЛЬ СИСТЕМЫ ЗАЩИТЫ ОТ ПОДДЕЛКИ ЖИДКИХ ПИЩЕВЫХ ФАСОВАННЫХ ПРОДУКТОВ // Материалы VIII Международной студенческой электронной научной конференции «Студенческий научный форум 2016»-URL: https ://w w w. scienceforum.ru/2016/pdf/26255.pdf
36. Батшев А.С., Белозеров В.В. МЕТОДОЛОГИЯ АВТОМАТИЗАЦИИ ЭКСПРЕСС-КОНТРОЛЯ ФАСОВАННЫХ МОЛОЧНЫХ ПРОДУКТОВ// Материалы IX Международной студенческой электронной научной конференции «Студенческий научный форум 2017»-URL: https ://w w w. scienceforum.ru/2017/pdf/30467.pdf
37. Белозеров В.В., Троицкий В.М., Белозеров Вл.В О МОДЕЛИ ИДЕНТИФИКАЦИИ КОНТРАФАКТА ЖИДКИХ ПИЩЕВЫХ ФАСОВАННЫХ ПРОДУКТОВ // Рациональное питание, пищевые добавки и биостимуляторы. - 2016. - № 1. - С. 26-36;-URL: http://www.journal-nutrition.rU/pdf/2016/l/35710.pdf
38. Белозеров В.В., Троицкий В.М. МОДЕЛЬ АВТОМАТИЗИРОВАННОЙ СИСТЕМЫ ЗАЩИТЫ ОТ ПОДДЕЛКИ ЖИДКИХ ФАСОВАННЫХ ПРОДУКТОВ //Международная научно-практическая конференция «Приоритетные задачи и стратегии развития технических наук»: сб. материалов - Н. Новгород: «Эвенсис», 2016, с.56-60.
39. Белозеров В.В., Батшев А.С., Любавский А.Ю. Об автоматизации идентификации жидких фасованных продуктов // Электроника и электротехника. — 2016.-№ 1,-С.135-145. DOI: 10.7256/2453-8884.2016.1.20924. URL: http://е-notabene.ru/elektronika/article_20924.html
