Многокритериальная оптимизация процесса отделения корки от мякоти дыни методами экспериментально-статистического анализа Текст научной статьи по специальности «Прочие технологии»

BecmHunJBry^T/Proceedings of VSUET, № 2, 2016L

Оригинальная статья/Original article_

УДК 635.611, 519.863

DOI: http://doi.org/10.20914/2310-1202-2016-2-28-36

Многокритериальная оптимизация процесса отделения корки от мякоти дыни методами экспериментально-статистического анализа

Евгений Б. Медведков, 1 [email protected]

Маржан Е. Кизатова, 1 [email protected]

Александр А. Шевцов, 2 [email protected]

Александр С. Муравьев 2 [email protected]

1 кафедра механизации и автоматизации технологических процессов, Алматинский техн. ун-т, ул. Толе би, 100, г. Алматы, Казахстан

2 кафедра технологии хлебопекарного, кондитерского, макаронного и зерноперерабатывающего производств, Воронеж. гос. ун-т. инж. техн., пр-т Революции, 19, г. Воронеж, Россия

Реферат. Научная и инновационная деятельность по созданию энергоэффективных процессов при переработке дыни, является актуальной задачей. Одной из трудоёмких операций в данной технологии является отделение корки от мякоти с их последующим целевым применением в создании новых продуктов питания. Недостаточная научно-экспериментальная база этой операции сдерживает разработку высокоэффективных машин для её реализации. В связи с этим предлагается методика проведения эксперимента по отделению корки от дыни на опытной установке и поиск оптимальных режимов её работы методами статистического моделирования. В качестве объектов исследования были взяты дыни позднеспелых сортов: Калайсан, Торлами, Гуляби-сары. Для исследования взаимодействия факторов, влияющих на процесс отделения корки от мякоти дыни, применяли центральное композиционное ротатабельное униформпланирование и дробный факторный эксперимент. Использование метода планирования эксперимента с обработкой матрицы планирования в программном модуле Design Expert v.10, позволило получить уравнения регрессии, адекватно описывающие реальный процесс. Установлены рациональные интервалы значений отношения частоты вращения истирающего барабана к частоте вращения подающего барабана; величины зазора между подающим барабаном и срезающим ножом; угла заточки срезающего ножа; количества шипов подающего барабана; диаметра отверстий истирающего барабана. Среднеквадратичная ошибка не превышала 12,4% Графическая интерпретация уравнений регрессии представлена кривыми равных значений и инженерными номограммами, которые позволяют прогнозировать выбор рациональных значений входных факторов по трём критериям: минимальным значениям удельных энергозатрат на процесс срезания, максимальным значениям удельной производительности по мякоти и коэффициенту извлечения мякоти. Полученные данные могут быть использованы при оперативном управлении технологическими параметрами процесса с учётом геометрических размеров дыни, и её неоднородного состава.

Ключевые слова: дыня, оптимизация, мякоть дыни, многофакторный анализ_

Multi-criteria optimization of the flesh melons skin separation process _by experimental and statistical analysis methods_

Yevgenij B. Medvedkov, 1 [email protected]

Marjane E. Kizatova, 1 [email protected]

Aleksandr A. Shevtsov, 2 [email protected]

Aleksandr S. Muravev 2 [email protected]

1 Department of mechanization and automation of manufacturing processes, Almaty technological university, Tole bi, 100, Almaty, Kazakhstan

2 Department of bakery technology, confectionery, pasta and grain processing industries, Voronezh state university of engineering technologies, Revolution Av., 19 Voronezh, Russia_

Summary. Research and innovation activity to create energy-efficient processes in the melon processing, is a significant task. Separation skin from the melon flesh with their subsequent destination application in the creation of new food products is one of the time-consuming operations in this technology. Lack of scientific and experimental base of this operation holding back the development of high-performance machines for its implementation. In this connection, the technique of the experiment on the separation of the skins of melons in the pilot plant and the search for optimal regimes of its work methods by statistical modeling is offered. The late-ripening species of melon: Kalaysan, Thorlami, Gulab-sary are objects of study. Interaction of factors influencing on separating the melon skins process is carried out. A central composite rotatable design and fractional factorial experiment was used. Using the method of experimental design with treatment planning template in Design Expert v.10 software yielded a regression equations that adequately describe the actual process. Rational intervals input factors values are established: the ratio of the rotational speed of the drum to the abrasive supply roll rotational frequency; the gap between the supply drum and the shearing knife; shearing blade sharpening angle; the number of feed drum spikes; abrading drum orifices diameter. The mean square error does not exceed 12.4% Regression equations graphic interpretation is presented by scatter plots and engineering nomograms that can be predictive of a choice of rational values of the input factors for three optimization criteria: minimal specific energy consumption in the process of cutting values, maximal specific performance by the pulp and pulp extraction ratio values. Obtained data can be used for the operational management of

the process technological parameters, taking into account the geometrical dimensions of the melon and its inhomogeneous structure._

Keywords: melon, optimization, melon pulp, response surface analysis_

Для цитирования

Медведков Е. Б., Кизатова М. Е., Шевцов А. А., Муравьев А. С. Многокритериальная оптимизация процесса отделения корки от мякоти дыни методами экспериментально-статистического анализа // Вестник ВГУИТ. 2016. № 2. С. 28-36. (М:10.20914/2310-1202-2016-2-28-36

For citation

Medvedkov Y. B., Kizatova M. E., Shevtsov A. A., Muravev A. S. Multi-criteria optimization of the flesh melons skin separation process by experimental and statistical analysis methods. Vestnik VSUET [Proceedings of VSUET]. 2016. no 2 pp. 28-36. (in Russ.). doi:10.20914/2310-1202-2016-2-28-36

Введение

Структурная перестройка экономики и активизация инновационной и инвестиционной деятельности являются стратегическими факторами роста и играют фундаментальную роль в разработке прорывных технологий и оборудования в перерабатывающих отраслях промышленности. Реализация программ по переводу пищевой промышленности на наукоёмкий тип технико-технологического развития опирается на новые научные достижения и их практическую реализацию. Научная и инновационная деятельность по созданию энергоэффективных процессов при переработке бахчевых культур, в том числе дыни, является актуальной задачей [1-2]. Одной из трудоёмких операций в разветвлённой технологии обработки дынь является отделение корки от мякоти с их последующим целевым применением в создании новых продуктов питания. Недостаточная научно-экспериментальная база этой технологической операции сдерживает разработку высокоэффективных машин для её реализации. В связи с этим предлагается методика проведения эксперимента по отделению корки от дыни на опытной установке и поиск оптимальных режимов её работы методами статистического моделирования.

1.1 Материалы и методы

Экспериментальная установка для отделения корки от дыни (рисунок 1) содержит: истирающий барабан, подающий барабан с шипами, срезающий нож приёмный жёлоб, жёлоб для отвода мякоти, приводы истирающего и падающего барабанов, станину.

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

В качестве объектов исследования были взяты дыни следующих позднеспелых сортов: Калайсан, Торлами, Гуляби-сары, так как именно позднеспелые сорта дыни используются для переработки и выработки из них различных продуктов длительного хранения согласно классификации сортов дыни, предложенной А. М. Адмаевой [3-4].

пДиющии j

барабан | /' с шчпишхЧу^

kSupply drum \r

kwilh spikes /

^^>лектродвигат&ть жёлоб для

отвода шжетТК motor Щ/

К1 pulp guner

Рисунок 1. Экспериментальная установка для отделения корки от дыни

Figure 1. Expérimental installation to separate the skin from melon

Кроме того, эти сорта дыни имеют высокую степень лёжкости, по сравнению с раннеспелыми и позднеспелыми сортами, что позволяет продлить сроки переработки сырья, а, следовательно, повысить эффективность производительности линии по переработке дыни. Именно такие сорта имеют овально-удлинённую форму, позволяющую в максимальной степени механизировать процесс очистки дыни с одновременной её резкой на куски заданного размера.

1.2 Оптимизация процесса отделения мякоти дыни от корки

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

Математическая модель изучаемого процесса представляется в виде полинома второй степени:

Y=ъ0 + îb1X1 +:tb,Xf + i bjXiXj (1)

i=1 i=1 i<j

где bo - свободный член уравнения, равный средней величине отклика при условии, что рассматриваемые факторы находятся на средних, «нулевых» уровнях; X - масштабированные значения факторов, которые определяют функцию отклика и поддаются варьированию; bij -коэффициенты двухфакторных взаимодействий,

BeemHunJBTy^T/Proceedings of VSUET, № 2, 20Î6L

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

1.2.1 Обоснование выбора пределов изменения входных факторов

В качестве основных факторов, влияющих на процесс срезания корки дыни, были выбраны:

Х1 - отношение частоты вращения истирающего барабана к частоте вращения подающего барабана, с-1/с-1;

Х2 - величина рабочего зазора между подающим барабаном и срезающим ножом, мм; Х3 - угол заточки срезающего ножа, град; Х4 - количество шипов подающего

барабана в зоне захвата, шт.;

Х5 - диаметр отверстий истирающего барабана, мм.

Все эти факторы совместимы и некоррелированны между собой. Выбор интервалов изменения факторов обусловлен технологическими условиями. Пределы изменения исследуемых факторов приведены в таблице 1.

Пределы изменения входных факторов Range of input factors variation

Таблица 1.

Table 1.

Условия планирования Planning conditions Кодированное значение Coded value Пределы изменения факторов Range of input factors variation

Х1, с-1/с-1 s-1/s-1 Х2, мм mm Х3, град deg Х4, шт. psc Х5, мм mm

Основной уровень Center point 0 10,0 5 40 6 4

Верхний уровень Low level +1 12,5 7 55 7 5

Нижний уровень High level -1 7,5 3 25 5 3

Верхняя «звёздная точка» Top axial point +2 15,0 9 70 8 6

Нижняя «звёздная точка» Bottom axial point -2 5,0 1 10 4 2

Выбор интервалов изменения факторов обусловлен технологическими условиями процесса отделения корки от дыни [5], а также технико-экономическими показателями.

Критериями оценки влияния различных факторов на процесс срезания были выбраны:

Г1 - удельные энергозатраты на процесс срезания, Вт/кг;

У2 - удельная производительность по мякоти, г/с;

У3 - коэффициент извлечения мякоти (отношение количества извлекаемой мякоти к общей массе дольки), %.

Программа исследования была заложена в матрицу планирования эксперимента (таблица 2).

Для исследования применили центральное композиционное ротатабельное униформ-планирование и дробный факторный эксперимент ДФЭ 25"1. Число опытов в матрице планирования пяти входных параметров равно 32.

1.2.2 Уравнения регрессии и их анализ

После обработки матрицы планирования в программном модуле Design Expert v.10 [5] были получены следующие уравнения регрессии:

Y = 1,45 + 0,046X1 - 0,041X2

- 0,043X3 + 0,035X4 - 0,09X5 - 5,00E-03 • X12 - 0,113X22 -

-0,150X32 - 0,125X42 - 0,127X52

(2)

Y2 = 69,67 + 5,54X1 - 2,21X2 - 2,29X3 + 2,04X4 +1,96X5 - 0,69X2X3 + 0,56X2X4 --1,06X3X5 -0,81X4X5 +0,08X12 +0,20X22 +0,20X32 -0,80X42 +0,58X52

(3)

Y3 = 89,86 +1,25X1 - 2,67X2 -1,33X3

1,17X4 -1,58X5 -0,25XX -1,25XX -1,00XX

+1,25X2X3 - 0,62X2X4 + 0,75X2X5 - 0,87X3X4 - 0,24X22 - 0,25X32 - 0,86X52

Графическая интерпретация уравнений (2)-(4) в виде кривых равных значений и поверхности отклика приведена на рисунках 2-4. Установлено влияние входных факторов Х1,7 = (1,5)

на выходные Уу, у = (1,3) (рисунки 5-7). Номограммы, представляющие практический интерес, представлены на рисунках 8-9.

Y„ Вт/кг W/kg

11-5 х1( cVc1

s'/s'

(a) (b)

Рисунок 2. Кривые равных значений (a) и поверхность отклика (b) зависимости удельных энергозатрат на процесс срезания Y1, Вт/кг от: отношения частоты вращения истирающего барабана к частоте вращения подающего барабана Х1, с-1/с-1 и угла заточки срезающего ножа Х3, град

Figure 2. Contour plot (a) and surface response (b) depending on the specific energy consumption in the process of cutting Y1, W/kg by: relationships abrasive drum speed to speed the delivery drum Х1, s-1/s-1 and the sharpening angle of the shearing knife Х3, deg

у2, г/С g/s

(a) (b)

Рисунок 3. Кривые равных значений (a) и поверхность отклика (b) зависимости удельной производительность по мякоти Y2, кг/с от: величины рабочего зазора между подающим барабаном и срезающим ножом Х2, мм и диаметром отверстий истирающего барабана Х5, мм

Figure 3. Contour plot (a) and the response surface (b) depending on the specific performance of the pulp Y2, kg/s of: working gap between the supply drum and the shearing knife Х2, mm and abrading drum orifices diameter Х5, mm

Y3,%

XJ( мм mm

36 Х,,град deg

(a)

(b)

Рисунок 4. Кривые равных значений (a) и поверхность отклика (b) зависимости коэффициента извлечения мякоти Y3, % от: величины рабочего зазора между подающим барабаном и срезающим ножом Х2, мм и угла заточки срезающего ножа Х3, град

Figure 4. Contour plot (a) and the response surface (b) depending on the pulp extraction ratio Y3, % of: working gap between the supply drum and the shearing knife Х2, mm and shearing knife sharpening angle Х3, deg

(a)

(b)

X^uiT 4 Xs.mm 2

12.5 15 s-'/S"1

7 9 mm

70

б 7 8 P«

4 5 6 mm

(c)

Рисунок 5. Влияние удельных энергозатрат на процесс срезания Y1, Вт/кг (a), удельной производительность по мякоти Y2, кг/с (b), коэффициента извлечения мякоти Y3, % (c) от входных факторов X

Figure 5. Impact of the input factors Xt on specific energy consumption in the process of cutting, Y1 W/kg (a), pulp density performance Y2, kg/s (b) and pulp recovery ratio Y3 % (c)

3

Л{ К ! к \ ; \ и]_ LJy£J ТДз1 V

\ \ м \ 4 \ '' ' [91] 1

\ ^Wi \ [f

\ \ . \ \ '[88 Ш ГШ] т гид! я \ .\ \ \ \

4 'Ш ' ГСТП \ \ -. ■

10

Ж,, cVc1 s'Vs"1 -Y2-

Рисунок 6. Номограмма для определения зависимости отношения частоты вращения истирающего барабана к частоте вращения подающего барабана (Х1), с-1/с-1 и количества шипов подающего барабана в зоне захвата (Х4), шт. от выходных факторов Y'

Figure 6. The nomogram for the determination of the abrasive drum speed to the feed speed of the drum ratio ( Х1 ) , s-1/s-1 and the amount of feed drum spikes in grip zone (Х4), psc by output factors Y j

1.2.3 Поиск оптимальных режимных параметров процесса

Нахождение оптимума и решения задачи оптимизации проводили с использованием функции желательности d [...], при этом каждый Y переводили в dt е{0...1} . Если Y определяет оптимум целевой функции (максимальное или минимальное значение), или же оптимальное значение является заданным, обусловленным экономической целесообразностью, то di = 1. В случае когда оптимальное значение не соответствует интервалу заданных значений Yj £ {0.1} , то di = 0 . После определения d по каждому из Yj , определяли максимальное значение функции желательности по формуле:

D = (dd ...dj )(1/j), (5)

По сравнительной оценке множества полученных значений функции желательности D выбирали оптимальное (единственное) решение

Рисунок 7. Номограмма для определения зависимости количества шипов подающего барабана в зоне захвата (Х4), шт. и диаметр отверстий истирающего барабана (Х5), мм от выходных факторов Yj

Figure 7. The nomogram for determining of the amount of feed drum spikes in grip zone (Х4), psc and abrasive drum orifices diameter (Х5), mm by output factors Y'

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

Функцию желательности ё при решении задачи оптимизации делили на три частных случая:

1. Выходной фактор стремится к максимальному значению:

d =

0

У - L T - L 1

У < L L < y < T, У > T

(6)

где Т - целевое значение; V, Ь - нижняя и верхняя границы интервала выходного фактора соответственно; г - показатель «влияния» функции желательности (весовой коэффициент), при г = 1 - ё является линейной функцией; г > 1 оказывает большее влияние на смещение ё в сторону целевого значения; 0 < г < 1 - менее важная степень влияния.

2. Выходной фактор стремится к минимальному значению:

d =

1

U - y

y < T T < y < U , (7)

U - T

0 y > U

3. Выходной фактор стремится к заданному числовому значению:

d =

1

y - L U - T

U - y U - T 0

y < L L > y > T

T < y < U y > U

(8)

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

Так как выбранные критерии оптимизации одинаково важны и носят характер равнозначности, то показатели их «влияния» на функцию желательности принимали равными единице, г = 1.

Из массива решений Б выбраны значения В ^ 1 (рисунок 8), по которым установлены рациональные интервалы значений входных факторов: Х1 = 8...10 с-1/с-1; Х2 = 1_3 мм; X3 = 11.13 град; Х4 = 6.7 шт.; Х5 = 5.6 мм.

Рисунок 8. Определение функции желательности D для заданных критериев оптимизации Figure 8. Determination of the desirability function D for the specified criteria optimization 34

Для проверки правильности полученных результатов был поставлен ряд параллельных экспериментов. Полученные результаты попадали в рассчитанные доверительные интервалы по всем критериям качества. Среднеквадратичная ошибка не превышала 12,4%.

Заключение

Использование метода планирования эксперимента с обработкой матрицы планирования в программном модуле Design Expert v.10 позволило получить уравнения регрессии, адекватно описывающие реальный процесс.

ЛИТЕРАТУРА

1 Sannikova T. A. et al. Getting candied fruits from melon crops // Technologies for the food and processing industry of AlC-healty food. 2015. №26 P. 45-49.

2 Villanueva M. J., Tenorio M. D., Esteban M. A., Mendoza M. C. Compositional changes during ripening of two cultivars of muskmelon fruits // Food Chemistry. 2004.V. 87. № 2. P. 179-185. URL: http://dx.doi.org/10.1016/j .foodchem.2003.11.009

3 Медведков E. Б., Кизатова M. E., Ад-маева А. М., Азимова С. Т. и др. Содержание и структура пектина в плодах дыни // Новое слово в науке и практике: гипотезы и апробация результатов исследований. 2015. № 21. С. 103-107.

4 Медведков Е. Б., Кизатова М. Е., Ад-маева А. М., Байболова Л. К. и др. Упруго-пластические характеристики мякоти и мезги дыни при сжатии // Достижения вузовской науки. 2015. № 19 С. 99-104.

5 Пеленко В. В., Зуев Н. А., Ольшевский Р. Г., Иваненко В. П. и др. Оценка зависимости производительности измельчителей мяса от их конструкции и физико-механических свойств сырья. // Вестник Международной академии холода. 2015. № 1. С. 9-15.

6 Myers R. Н., Montgomery D. С., Ап-derson-Cook C. M. Response Surface Methodology: Process and Product Optimization Using Designed Experiments. Wiley (Wiley Series in Probability and Statistics), 2009. 680 p.

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

REFERENCES

1 Sannikova T. A. et al. Getting candied fruits from melon crops. Technologies for the food and processing industry of AlC-healty food. 2015, no 26, pp.45-49.

2 Villanueva M. J., Tenorio M. D., Esteban M. A., Mendoza M. C. Compositional changes during ripening of two cultivars of muskmelon fruits. Food Chemistry. 2004, vol. 87, no 2, pp. 179-185. Available at: http://dx.doi.org/10.1016/j .foodchem.2003.11.009

3 Medvedkov E. B., Kizatova M. E., Ad-maeva A. M., Azimova S. T. et al. The content and structure of the pectin in the fruit melon. Novoe slovo v nauke i praktike: gipotezy i aprobatsiya re-zultatov issledovanii. [New word in the science and practice: testing hypotheses and research results], 2015, no 21, pp. 103-107. (in Russian).

4 Medvedkov E. B., Kizatova M. E., Ad-maeva A. M., Baibolova L. K. et al., Elasto-plastic characteristics of pulp and melon mash compressive. Dostizheniya vuzovskoji nauki. [Achievements university science], 2015, no 19, pp. 99-104. (in Russian).

5 Pelenko V. V., Zuev N. A., Olszewsky R. G., Ivanenko V. P. et al. Evaluation of meat grinder performance depending on their design and physical and mechanical properties of raw-material. Vestnik Mezhdunarodnoi akademii kholoda. [Proceedings of the International Academy of Refrigeration] 2015, no. 1, pp. 9-15. (in Russian).

6 Myers R. H., Montgomery D. C., Anderson-Cook C. M. Response Surface Methodology: Process and Product Optimization Using Designed Experiments. Wiley (Wiley Series in Probability and Statistics), 2009, 680 p.

СВЕДЕНИЯ ОБ АВТОРАХ

Евгений Б. Медведков д. т. н., профессор, кафедра механизации и автоматизации технологических процессов, Алматинский техн. ун-т, ул. Толе би, 100, г. Алматы, Казахстан, [email protected]

Маржан Е. Кизатова докторант PhD, кафедра механизации и автоматизации технологических процессов, Алматинский техн. ун-т, ул. Толе би, 100, г. Алматы, Казахстан, [email protected]

Александр А. Шевцов д. т. н., профессор, кафедра технологии хлебопекарного, кондитерского, макаронного и зерноперерабатывающего производств, Воронеж. гос. ун-т. инж. техн., пр-т Революции, 19, г. Воронеж, Россия, [email protected] Александр С. Муравьев аспирант, кафедра технологии хлебопекарного, кондитерского, макаронного и зернопе-рерабатывающего производств, Воронеж. гос. ун-т. инж. техн., пр-т Революции, 19, г. Воронеж, Россия, [email protected]

КРИТЕРИЙ АВТОРСТВА

Евгений Б. Медведков предложил методику проведения эксперимента и организовал производственные испытания Маржан Е. Кизатова обзор литературных источников по исследуемой проблеме, провела эксперимент, выполнила расчеты

Александр А. Шевцов консультация в ходе исследования Александр С. Муравьев написал рукопись, корректировал её до подачи в редакцию и несет ответственность за плагиат

КОНФЛИКТ ИНТЕРЕСОВ

Авторы заявляют об отсутствии конфликта интересов. ПОСТУПИЛА 15.04.2016

ПРИНЯТА В ПЕЧАТЬ 16.05.2016

INFORMATION ABOUT AUTHORS

Yevgeniy B. Medvedkov D. t. sc., professor, mechanization and automation of manufacturing processes department, Almaty technological university, Tole bi, 100, Almaty, Kazakhstan, [email protected] Marzhan E. Kizatova Ph. D. student , mechanization and automation of manufacturing processes department, Almaty technological university, Tole bi, 100, Almaty, Kazakhstan, [email protected] Aleksandr A. Shevtsov D. t. sc., professor, bakery technology, confectionery, pasta and grain processing industries de-partment, Voronezh state university of engineering technologies, Revolution Av., 19 Voronezh, Russia, she-valol@rambler. ru

Aleksandr S. Muravev graduate student, bakery technology, confectionery, pasta and grain processing industries de-partment, Voronezh state university of engineering technologies, Revolution Av., 19 Voronezh, Russia, [email protected]

CONTRIBUTION

Evgenij B. Medvedkov proposed a scheme of the experiment and organized production trials Marjane E. Kizatova review of the literature on an investigated problem, conducted an experiment, performed computations

Aleksandr A. Shevtsov consultation during the study Aleksandr S. Muravev wrote the manuscript, correct it before filing in editing and is responsible for plagiarism

CONFLICT OF INTEREST

The authors declare no conflict of interest.

RECEIVED 4.15.2016

ACCEPTED 5.16.2016