ISSN 2074-9414 (Печать),
ISSN 2313-1748 (Онлайн)

ИССЛЕДОВАНИЕ ФУНКЦИОНИРОВАНИЯ ЦЕНТРОБЕЖНОГО СМЕСИТЕЛЯ НЕПРЕРЫВНОГО ДЕЙСТВИЯ МЕТОДОМ МНОЖЕСТВЕННОГО РЕГРЕССИОННОГО АНАЛИЗА

Аннотация
Целью данной статьи является разработка новой конструкции центробежного смесителя, обладающего регулируемой инерционностью, для получения сухих многокомпонентных композиций с соотношением смешиваемых компонентов от 1:75 до 1:125. Исследованы две конструкции одноконусных центробежных смесителей непрерывного действия для получения смесей сухих дисперсных материалов. На основе анализа множественной регрессии определены рациональные значения некоторых технологических и конструктивных параметров, оказывающих влияние на качество исследуемых смесей в модернизированном аппарате по сравнению с базовой конструкцией.
Ключевые слова
Центробежный смеситель, сыпучий материал, множественная регрессия
ВВЕДЕНИЕ
Введение Концепции государственной политики в области здорового питания населения РФ предусматривают создание пищевой продукции, обогащенной витаминами, минеральными и биологически активными добавками. При этом главное внимание обращается на качество пищевых продуктов и их соответствие медико-биологическим требованиям. Это особенно важно для Кузбасса, где неблагоприятная экологическая обстановка осложняется несбалансированностью рациона и отсутствием в нем необходимого количества витаминов, микро- и макроэлементов, что вызывает необходимость обогащения продуктов питания биологически ценными компонентами и освоение новых технологий их получения. Аналогичную проблему приходится решать в аграрно-промышленном комплексе при производстве комбикормов, ЗЦМ сухим способом, премиксов, БВД, различных шихт для получения стекла, искусственных алмазов и композиционных материалов. Поэтому разработка эффективных, непрерывно действующих смесительных агрегатов нового поколения для переработки мелкозернистых и дисперсных материалов с большой разницей содержания их в смеси является актуальной научной задачей, представляющей практический интерес для сельскохозяйственных, пищевых, химических и строительных отраслей. Процессы смешения широко используются во многих отраслях промышленности для приготовления сыпучих композиций из дисперсных материалов, которые можно осуществлять в непрерывном или периодическом режимах. Переход на аппаратурное оформление стадии смешивания по непрерывной схеме дает возможность автоматизировать процесс, уменьшить загрязнение производственных помещений и окружающей среды пылевыми выбросами, а также увеличить производительность труда при одновременном снижении энергопотребления, металлоемкости и себестоимости смеси. Например, конструкции, разработанные с участием авторов [4-6], обладают производительностью 400 кг/ч при энергопотреблении не более 0,6 (Вт×ч)/кг, в то время как у применяемого в промышленности оборудования этот показатель 09÷1,1 (Вт×ч)/кг. Повышение эффективности и интенсивности смешивания достигается за счет проведения процесса в тонких и разреженных слоях. К такому оборудованию относятся центробежные смесители с рабочим органом в виде усеченного конуса, закрепленным на вертикальном валу. Целью данной статьи является разработка новой конструкции центробежного смесителя (ЦС), обладающего регулируемой инерционностью, для получения сухих многокомпонентных композиций с соотношением смешиваемых компонентов от 1:75 до 1:125 на основе анализа результатов теоретических и экспериментальных исследований процесса смешивания. На основе обработки результатов были получены математические модели ЦС новой конструкции, позволяющие проанализировать возможность получения смесей заданного качества с учетом его инерционных, геометрических и технологических параметров работы. Объекты и методы исследований В качестве объектов исследований взяты два одноконусных центробежных смесителя непрерывного действия. Первый аппарат состоит из элементов 1-8, второй - дополнительно содержит элементы 9-11 (рис. 1). Рис. 1. Одноконусный центробежный смеситель непрерывного действия: 1 - корпус; 2 - крышка; 3 - загрузочные патрубки; 4 - подшипниковый узел; 5 - приводной вал; 6 - коническое днище; 7 - разгрузочный патрубок; 8 - конический ротор; 9 - углообразные турбулизаторы; 10 - перфорированные направляющие; 11 - рассеивающий диск В первом смесителе [1] на вращающемся валу 5 закреплен усеченный полый гладкий конус 8, внутрь которого на малое основание подают ингредиенты через загрузочные патрубки 3. За счет центробежной силы частицы начинают спиралеобразно скользить по внутренней поверхности конуса и смешиваться в разреженном слое. Готовый продукт сбрасывается на стенки корпуса 1, оседает вдоль них, скатывается по коническому днищу 6 и удаляется из аппарата через разгрузочный патрубок 7. В модернизированном смесителе на внутренней поверхности конуса 8 дополнительно установлены углообразные турбулизаторы 9, которые усложняют траектории движения частиц сыпучего материала по его поверхности. Полученная композиция оседает по внутренней поверхности корпуса 1 на перфорированные направляющие 10. Часть потока стекает по направляющим на рассеивающий диск 11, закрепленный на внешней поверхности ротора 8, вторая - проходит через отверстия, выполненные в теле направляющих 10. С помощью диска 11 частицы отбрасываются центробежной силой к стенкам корпуса 1, где пересекаются с потоком, прошедшим через отверстия. Элементы 9-11 обусловливают дополнительное смешивание ингредиентов и тем самым повышают качество готового продукта. С целью выявления степени влияния числа турбулизаторов h, частоты вращения ротора n и соотношения ингредиентов С на качество двухкомпонентных смесей (поваренная соль - манная крупа и сахар-песок - пшено), оцениваемое коэффициентом вариации Vc: где S - выборочное среднее квадратическое отклонение содержания ключевого компонента в пробах, взятых из смеси; - выборочное среднее значение концентрации ключевого компонента в пробах. Методики планирования экспериментов и обработки их результатов изложены в [2]. Результаты и их обсуждение В ходе экспериментов варьировались следующие параметры: количество углообразных отражателей h в диапазоне (4÷12 шт.), частота вращения ротора n (750÷950 об/мин), соотношение смешиваемых компонентов C (1:75÷1:125). Базовые (нулевые) точки и шаги варьирования приведены в табл. 1. Согласно матрице плана с соблюдением рандомизации принимались режимные и конструктивные параметры аппарата. Затем устанавливали требуемую частоту вращения ротора и подавали ингредиенты. После выхода смесителя на стационарный режим работы осуществляли отбор проб на выходе из аппарата в количестве 30 шт. для каждого эксперимента. Таблица 1 Уровни и интервалы варьирования Фактор и его обозначениеВерхний уровеньНижний уровеньЦентр планаИнтервал варьирования Количество углообразных турбулизаторов h12484 Частота вращения ротора n950750850100 Соотношение смешиваемых компонентов C1:1251:751:1001:25 Далее анализировали качество проб смеси с помощью ситового анализа, по различию гранулометрического состава, химическим способом с помощью поляриметрического метода (определение сахарозы) и ионометрическим методом (определение концентрации соли). Для получения более точного значения измеряемого параметра на каждом режиме работы проводилась серия опытов от 3 до 5. Усредненные результаты коэффициента неоднородности для исследуемых смесей (соль - манная крупа, сахар - пшено) приведены в табл. 2 и 3. Таблица 2 Результаты экспериментов, полученные на первой конструкции ЦС Названия факторов и их значенияКоэффициент неоднородности Vc, %, для смесей Cnhсоль - манная крупасахар - пшено 1:75750019,718,2 1:125750023,422,8 1:100750022,523,4 1:75950018,317,6 1:125950024,223,7 1:100950023,722,9 1:75850019,218,5 1:125850024,121,4 1:100850022,119,8 При сопоставлении табл. 2 и 3 видно, что смеситель модернизированной конструкции позволяет получать смеси с коэффициентом неоднородности (Vc ≈ 10 %) почти в 2 раза меньшим, чем на базовой конструкции (Vc ≈ 19 %). Это свидетельствует о положительном влиянии углообразных турбулизаторов и рассеивающего диска на конус СНД при смешении обоих смесей. Столь большая величина Vc обусловлена конструкцией смесителя, которая имеет всего 4 углообразных турбулизатора в роторе, состоящем всего из одного конуса, а также соотношением смешиваемых компонентов 1:125. Кроме того, при больших частотах вращения ротора частицы исследуемых материалов подвергаются частичному измельчению. Таблица 3 Результаты экспериментов, полученные на второй конструкции ЦС Названия факторов и их значенияКоэффициент неоднородности Vc, %, для смесей Cnhсоль - манная крупасахар - пшено 1:75750415,2413,16 1:125750417,9321,34 1:100750416,8415,77 1:75950414,3612,35 1:125950417,1619,54 1:100950415,5214,15 1:75850414,9412,97 1:125850417,3620,35 1:100850415,8814,94 1:75750811,710,17 1:125750813,7516,83 1:100750812,6711,72 1:7595089,1610,01 1:125950811,517,07 1:100950810,1111,46 1:7585089,769,65 1:125850811,9114,84 1:100850811,2910,98 1:757501210,379,6 1:1257501213,613,59 1:1007501212,4511,46 1:75950129,879,2 1:1259501211,9313,59 1:1009501210,5311 1:75850129,479,52 1:1258501212,5513,83 1:100850129,6611,27 Под действием воздушных потоков образовавшиеся мелкие компоненты начинают витать над ротором, что приводит к сегрегации смеси и, как следствие, к ухудшению конечного продукта. Наилучшие показатели качества исследуемых смесей получились в результате смешения при частоте вращения ротора 850 об/мин, количестве углообразных турбулизаторов, равном 12, и соотношении смешиваемых компонентов 1:75. Для оценки функциональной зависимости между независимыми функциональными переменными Х1, Х2, Х3 (представленными в табл. 5 и 6) и откликом Vc воспользовались методом множественного регрессионного анализа [3]. Результаты обработки опытных данных для исследуемых смесей проведены в интегрированной системе комплексного статистического анализа Statistica и представлены в табл. 4. Таблица 4 Краткие сведения регрессионного анализа Значения анализа Смесь сахар - пшеноСмесь соль - манная крупа Множественный R0,945470,93461 Множественный R20,893910,87350 Скорректированный R20,883960,86164 F-критерий89,8754873,65557 p-значение0,000000,00000 Стандартная ошибка оценки1,556261,76015 В таблице даны показатели качества подобранной модели регрессии. R - коэффициент множественной корреляции. R2 - коэффициент детерминации - чрезвычайно важная характеристика, показывающая долю общего разброса (относительно среднезависимой переменной), которая объясняется построенной регрессией. Итак, например, для смеси сахар - пшено R2 равен 0,89391. Это значит, что построенная регрессия объясняет более 89 % разброса значений переменной относительно среднего. Далее рассмотрим значения F-критерия и уровня значимости p. F-критерий используется для проверки гипотезы о значимости регрессии. В нашем случае F = 89,9; p = 0,00. Следовательно, построенная регрессия высоко значима, так как уровень значимости для исследуемых смесей равен нулю. Стандартная ошибка оценки показывает рассеяние наблюдаемых значений относительно регрессионной прямой. Для того чтобы узнать, какие из независимых переменных (С, n или h) дают больший вклад в предсказание коэффициента неоднородности, необходимо изучить регрессионные (или В) коэффициен-ты. Для этого представим в табл. 5 и 6 результаты множественной регрессии. Таблица 5 Результаты множественной регрессии для смеси сахар - пшено BetaStd. Err.ВStd. Err.t(23)p-level Свободный член10,63,0153,530,001284 X1 - переменная С0,510,0570,1130,0128,910,000000 X2 - переменная n-0,0410,057-0,0020,003-0,710,478294 X3 - переменная h-0,7930,057-0,790,057-13,770,000000 Таблица 6 Результаты множественной регрессии для смеси соль - манная крупа BetaStd. Err.ВStd. Err.t(23)p-level Свободный член19,33,415,670,000003 X1 - переменная С0,270,060,0620,0144,320,000138 X2 - переменная n-0,100,06-0,0050,003-1,600,119085 X3 - переменная h-0,880,06-0,920,06-14,10,000000 В табл. 5 и 6 представлены стандартизованные коэффициенты Beta и нестандартизованные В регрессионные коэффициенты (веса), их стандартные ошибки Std. Err. Beta и Std. Err. B и уровни значимости p-level, а также t-критерий для проверки гипотезы о равенстве 0 свободного члена регрессии. Beta-коэффициенты оцениваются по стандартизованным данным, имеющим выборочное среднее 0 и стандартное отклонение 1. Поэтому величина Beta-коэффициентов позволяет нам сравнить соответственные вклады каждой из независимой переменной в предположение зависимой переменной Vc. Как видно из таблиц результатов, статически значимыми являются переменные С и h (их значения в таблицах представлены курсивом), значения которых самые большие, т.е. стремятся к единице. Именно они оказывают сильное влияние (в особенности h, при котором Beta = 0,8-0,9) на коэффициент неоднородности. Причем для обеих смесей регрессионный коэффициент h отрицателен; т.е. чем больше количество углообразных турбулизаторов, тем меньше Vc. Вклад в регрессию для С положителен, поэтому чем меньше соотношение смешиваемых компонентов, тем меньше коэффициент неоднородности Vc. Частота вращения ротора оказывает очень маленькое влияние на качество смешивания, по крайней мере в исследуемом диапазоне, поэтому данная величина Х2 была исключена из регрессионного анализа. Далее получим уравнения регрессии, описывающие влияние технологических и режимных параметров СНД на качество изготавливаемых смесей в кодированном и натуральном виде соответственно. Для смеси соль - манная крупа: Y = 16,7458+0,0227×X1-1,679×X3+ +0,0003×X21-0,0029×X1×X3+ +0,083×X23; (2) Vc = 16,7739+0,19987×h-0,01684×C+ +0,0192×h2-0,00023×h×C+0,000008×C2. (3) Для смеси сахар - пшено: Y = 24,4135+0,1904×X1-1,4975×X3+ +0,0014×X21-0,001×X1×X3+ +0,0552×X23; (4) Vc = 25,13753+1,61248×h-0,01497×C+ + 0,0896×h2-0,00008×h×C+0,000006×C2. (5) Для примера рассмотрим поверхности отклика (рис. 2 и 3), описываемые расчетными регрессионными уравнениями (2 и 4), за исключением малозначимого коэффициента Х2 (то есть зависимость частоты вращения ротора на Vc). Рис. 2. Поверхность отклика для смеси соль - манная крупа Рис. 3. Поверхность отклика для смеси сахар - пшено Анализируя полученные графические интерпретации, можно отметить, что количество углообразных турбулизаторов h, равное 8-12 шт., в роторе модернизированного СНД центробежного типа оказывает наибольшее влияние на коэффициент неоднородности Vс, особенно при h = 12. Наименьшее значение Vc достигается только в том случае, когда соотношение смешиваемых компонентов составляет 1:75, но и при 1:100 готовую смесь можно считать удовлетворительной, так как Vc для обеих смесей не превышал 12 %. Увеличение С ведет к ухудшению качества смеси, что объясняется малым временем пребывания частиц в роторе аппарата, так как последний состоит всего из одного конуса. Важным элементом регрессионного анализа является оценка адекватности модели, которая основывается на анализе удаленных остатков. Они определяются как стандартизованные остатки для соответствующих наблюдений, которые получились бы при исключении соответствующих наблюдений из анализа (исключение переменной Х2). При проведении анализа с помощью множественной регрессии предполагается, что значения остатков (наблюдаемые минус предсказанные значения) нормально распределены и функция регрессии (т.е. связь между коэффициентом неоднородности Vc и независимыми переменными С и h) в действительности линейна. Если данные предположения существенно нарушаются, то будут повышаться или понижаться оценки коэффициентов регрессии (В-коэффициентов), а также завышение или занижение значимости исследуемых переменных. Графики остатков приведены на рис. 4. Значения остатков нормально распределены, так как все точки на графиках не отклоняются от прямой и большинство из них не выходит за области граничных условий. Поэтому полученные модели для исследуемых смесей с более чем 95 % доверительности адекватно описывают опытные данные. Таким образом, при исследовании влияния на качество получаемых смесей конструктивных, технологических и режимных параметров СНД выявлено, что как отдельное воздействие каждой исследуемой переменной, так и их совместное влияние соответствует положительной корреляции с откликом. Выводы Для приготовления смесей сыпучих материалов с соотношением смешиваемых компонентов от 1:75 до 1:100 целесообразно применять модернизированную конструкцию одноконусного СНД центробежного типа. На основе анализа множественной регрессии определены параметры, оказывающие влияние на качество исследуемых смесей в модернизированном аппарате. Например, наибольшее влияние на коэффициент неоднородности оказывает число углообразных турбулизаторов (Beta = 0,8-0,9), меньшее - соотношение смешиваемых компонентов (Beta = 0,2-0,5) и фактически не оказывает влияния частота вращения ротора в исследуемом диапазоне. При соотношении смешиваемых компонентов 1:100 целесообразнее проводить процесс смешения в модернизированном ЦС при частоте вращения ротора от 750 до 950 об/мин, с установкой 12 углообразных турбулизаторов. Только при этих параметрах получаемые смеси можно считать удовлетворительными (Vc ≈ 10 %) по отношению к смесям, полученным на базовой конструкции (Vc ≈ 19 %). Рис. 4. Графики остатков
СПИСОК ЛИТЕРАТУРЫ
  1. А.с. 2121870 Россия, МКИ B01 F7/28. Центробежный смеситель порошкообразных материалов / В.И. Пындак, Г.Г. Соломатин, И.Ф. Рязанкин. - Опубл. в Б.И., 1998, № 32.
  2. Грачев, Ю.П. Математические методы планирования экспериментов / Ю.П. Грачев. - М.: Пищевая пром-сть, 1969. - 315 с.
  3. Елисеева, И.И. Статистика: учебник для вузов / И.И. Елисеева. - СПб.: Питер, 2010. - 360 с.
  4. Пат. 2200055 Российская Федерация, МПК В01 F7/26. Центробежный смеситель / Иванец В.Н., Ратников С.А., Иванец Г.Е., Бородулин Д.М., Белоусов Г.Н.; заявитель и патентообладатель ГОУ ВПО КемТИПП. - № 2000128252/12; заявл. 13.10.2000; опубл. 10.03.2003, Бюл. № 7. - 10 с.
  5. Пат. 2216394 Российская Федерация, МПК В01 F7/26. Центробежный смеситель / Иванец В.Н., Бакин И.А., Белоусов Г.Н., Бородулин Д.М., Волков А.С.; заявитель и патентообладатель ГОУ ВПО КемТИПП. - № 2002130106/12; заявл. 10.11.2002; опубл. 20.11.2003, Бюл. № 32. - 6 с.
  6. Пат. 2361653 Российская Федерация, МПК В01 F7/26. Центробежный смеситель / Ратников С.А., Бородулин Д.М., Селюнин А.Н., Сибиль А.В.; заявитель и патентообладатель ГОУ ВПО КемТИПП. - № 2008115038/15; заявл. 16.04.2008; опубл. 20.07.2009, Бюл. № 20. - 3 с.
Как цитировать?
О журнале