автореферат диссертации по технологии продовольственных продуктов, 05.18.12, диссертация на тему:Исследование процессов шелушения ячменя с целью создания малогабаритного шелушителя горизонтального типа

На правах рукописи

ГАРАНИН СЕРГЕЙ АЛЕКСАНДРОВИЧ

ИССЛЕДОВАНИЕ ПРОЦЕССОВ ШЕЛУШЕНИЯ ЯЧМЕНЯ С ЦЕЛЬЮ СОЗДАНИЯ МАЛОГАБАРИТНОГО ШЕЛУШИГГЕЛЯ ГОРИЗОНТАЛЬНОГО ТИПА

Специальность 05.18.12 - "Процессы и аппараты пищевых производств"

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

Москва 2005

Работа выполнена

В Государственном образовательном учреждении высшего профессионального образования «Московский государственный университет пищевых производств»

Научный руководитель - доктор технических наук, профессор

Глебов Леонид Александрович

Официальные оппоненты - доктор технических наук, профессор

Дулаев Валериан Георгиевич - кандидат технических наук, доцент Драгилсв Абрам Иосифович

Ведущая организация- ФГУПНИИ«МИР-ПРОДМАШ»

Защита диссертации состоится <2- %> 2005 г. в « У_2» часов на заседа-

нии Диссертационного совета Д.212.148.05 при ГОУ ВПО «Московский государственный университет пищевых производств» по адресу: 125080, Москва, Волоколамское шоссе, д.11., корпус А, аудитория 302.

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

С диссертацией можно ознакомиться в библиотеке МГУПГТ.

Автореферат разослан «£-4 2005 г.

Учёный секретарь

Диссертационного совета,

кандидат технических наук, доцент

( , А. С. Максимов

1м6

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы диссертации. Шелушение зернового сырья при производстве

( ^ ! ' - .1 /( !

комбикормов для молодняка скота и птицы является одной из важнейшей технологи-

1 ,(■ 1 г .[ , / >■

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

I !|| I 1 > I

качество и усвояемость. В настоящее время имеется целый ряд технологических ма-

' \ л го. .11

шин для шелушения зерновых культур, однако они имеют достаточно сложную конст-

' У, > "I I (/

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

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

1 Изучить технологические свойства зерна, влияющие на эффективнбеть шёлуШения;

2 Разработать математические модели рабочих процессов Для малогаб&рй'гного шелушителя с горизонтальным ротором;

3 Создать экспериментальную установку на базе малогабаритного шелушителя для проведения натурных экспериментов; • -1 *

4 Создать математическую модель процесса шелушения зерна и провеетИ 'экспериментальную проверку модели;

5. Сформулировать целевую функцию оптимальности процесса шелушения, поставить и решить задачу оптимальности процесса шелушения; „

6. Разработать алгоритм и создать программу для расчета на ГЖ оптимальных технологических и конструктивных параметров работы малогабаритного шелушителя; '

7 Усовершенствовать конструкцию шелушителя и провести промышленную апробацию результатов исследований.

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

РОС. НАЦИОНАЛЬНАЯ |

библиотека { 99

екая обработка результатов исследований, а также решение задачи оптимизации процесса шелушения проводились на ПК с использованием программ компьютерной математики MathCAD, SoIidWorks и CosmosFloWorks.

Научная новизна исследований заключается в следующем:

- в результате аналитических исследований получена и экспериментально проверена математическая модель рабочих процессов шелушителя с горизонтальным ротором;

- методами статистического анализа определено влияние технологических свойств зерна на количественные и качественные показатели процесса шелушения, что позволило в дальнейшем решить задачу оптимизации процесса шелушения;

- показана возможность использования компьютерной CAE - системы COSMOS-FIoWorks для моделирования кинетики зерновой смеси в рабочей зоне шелушителя;

- в результате статистического анализа экспериментальных данных получены эмпирические формулы зависимости эффективности шелушения ячменя от различных конструктивных и технологических параметров работы шелушителя, которые также позволили решить задачу оптимизации процесса шелушения;

- доказана возможность шелушения ячменя для кормовых целей в зоне свободного воздействия на него абразивных поверхностей;

- сформулирована целевая функция оптимальности процесса шелушения и методами нелинейного программирования на ПК определены оптимальные конструктивные параметры и рациональные режимы работы шелушителя для ячменя.

Объектом исследований является технологический процесс шелушения зерна, его режимы и закономерности.

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

Практическая значимость исследований:

1. Результаты теоретических и экспериментальных исследований получили практическую реализацию в совершенствовании оборудования для шелушения ячменя. Результаты исследований используются в учебном процессе;

2. Разработана математическая модель, алгоритм и компьютерная программа «Ше-лушитель» для определения оптимальных технологических параметров работы малогабаритного шелушителя.

3 Разработана конструкторская документация на малогабаритный шелушитель горизонтального типа, изготовлен опытный образец шелушителя, который установлен на ЗАО «Ремонтно-механичсский завод». Данные промышленных испытаний подтвердили высокие технико-экономические показатели малогабаритного шелушителя.

Апробация работы. Основные положения диссертационной работы докладывались на Международной научно - практической конференции «Индустрия продуктов здорового питания - третье тысячелетие» (Москва, МГУПП, 1999 г.); на расширенном заседании кафедры «Технологическое оборудование предприятий хлебопродуктов» (Москва, МГУПП, 2005г.).

На защиту выносятся: результаты исследований физико-механических свойств исходного продукта - зерна ячменя;

- математическая модель шелушителя горизонтального типа с абразивным ротором; эмпирические зависимости влияния различных конструктивных параметров шелушителя и технологических параметров исходного сырья на эффективности и энергоемкость процесса шелушения:

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

алгоритм и программа для определения оптимальных конструктивных и технологических параметров шелушителя горизонтального типа.

Публикации. По материалам исследований опубликовано 9 печатных работ.

Объем работы. Диссертационная работа состоит из введения, 5 глав, выводов, списка литературы и приложений. Она изложена на 174 страницах машинописного текста, основного текста 158 стр, содержит 57, иллюстраций и 8 таблиц. Список использованной литературы включает 212 наименований, в том числе 21 иностранный источник.

СОДЕРЖАНИЕ ДИССЕРТАЦИИ

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

^ '1

В первой главе приведён анализ научных публикаций и патентных материалов, отражающих современное состояние вопросов теории и практики шелушения зерно-

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

Большой вклад в развитие исследований в области изучения технологических, физико-механических и биохимических свойств зерна злаковых внесли ГА Еюров, М.Е. Гинзбург, Я М. Жазлин, Е М. Мельников, И Т. Мерко, П.Л. Шумилин, A.B. Лыков, В.В. Гортинский, В.В. Красников, Л.А. Трисвягский, а также зарубежные ученые М. М. Макмастерз, Дж Шелленбергер, В. Шеффер, Д. Джонс, Э. Калинский.

Исследованию микроструктуры зерна посвятили свои труды Бабиченко Л.В Беркутова Н.С , Швецова И.Л. Грищенко А.П., Изосимов В.П Петренко Т П., Попова Е.П., Pomeranz Y, Sacks I и другие.

Вопросами влияния гидротермической обработки (ГТО) на технологические свойства зерна в различные годы знанимались Абрамов С.Ю., Белиловская A.C., Берестов А.П., Гинсбург A.C., Гончаров З.Д., Гусев П Г., Егоров Г.А, Киреев В.М., Кретов С И., Лисицина Н В., и др. Многочисленные исследования показали, что в результате ГТО микроструктура анатомических частей зерна претерпевает необратимые изменения, которые усиливаются при повышении температуры процесса, а также в случае применения пара, СВЧ - обработки и т.д. При применении рациональных способов и режимов ГТО оболочки легче отделяются от ядра, ядро меньше дробится.

Огромный вклад в составление классификаций способов шелушения зерна был сделан такими учеными как Гинсбург М.Е., Егоров Г.А., Мельников Е.М., Панфилов В.А., Демский А.Б., Борискин М.А., Соколов А.Я. и др

Исследованию технологии шелушения зерновых культур и принципам действия оборудования, применяемого в процессе шелушения, посвятили свои работы отечественные и зарубежные ученые: Аммосова А М., Бутковский В. А., Гинсбург М Е., Глебов Л А., Гринберг E.H., Демский А Б., Дулаев В.Г., Егоров Г. А , Мельников Е.М., Жизлин Я.М., Каминский Д.В., Коломиец М.Г!.. Харламов С.В и др., в трудах которых раскрыты как общие, так и некоторые частные вопросы теории и конструирования рабочих органов шелушильных машин

Наибольший интерес представляют работы, посвященные исследованиям и разработке шелушильных машин, таких ученых как Я.М. Жизлин, E.H. Гринберг, А Я. Соколов, В.П. Зайцев, Б.Н. Гойхенберг, Д.В. Секерин, Е.В. Семенов, Л. А. Глебов, В.В. Гусев и пр.

На основе проведенного анализа обоснован выбор объекта исследования, сформулированы задачи диссертационной работы '

Вторая глава посвящена методологическим аспектам экспериментальных исследований Дано описание используемых в р'аботе приборов и экспериментальных установок, методик проведения зкспериментов' Лредсгавчены результаты исследований физико-механические свойства исходного сырья - ,ерна ячменя '

г, ' /' , | _л( 1

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

Для исследования влияния различных конструктивных и технологических параметров малогабаритного шелушителя с горизонтальным абразивным ротором на эффективность процесса шелушения ячменя в лаборатории кафедры ТОПХ МГУПГТ была создана экспериментальная установка, состоящая и* электродвигателя мощностью 7,5 кВт, цепного вариатора скоростей и шелушителя, оборудованного системой контроля частоты вращения вала и системой контроля потребляемой мощности. Экспериментальный образец шелушителя с горизонтальным ротором разработан на кафедре ТОПХ и изготовлен на опытном заводе АО «Мир» Схема экспериментальной уста-

хо и про 1\1аа ^ - съемный б\нкер для сбора отходов шелушения (ра!би1 на 5 секций) 6- сито для вывода продуктов шелушбния 7 ' йыпускной пагр\бок, 8-бункер для сбора шечушенмого )ерна,9- )лек-гродвш-атсть, 10- клиноременная передача. 11 - вариатор 12-унругая вт>точно-пальцевая муфта

Прочностные характеристики исходного сырья исследуются для определения предельно-допустимых усилий воздействия рабочих органов шелушителя на обрабатываемое зерно В настоящей работе среднее усилие разрушения зерновки ячменя при статическом сжатии определялось по методике, разработанной проф. Глебовым ЛАВ соответствии с этой методикой отбираются три навески по 100 г каждая, которые рассортировываются по крупности на пять фракций. Зерно, полученное сходом с сит с отверстиями размером 3,0x20; 2,8x20; 2,5x20; 2,0x20 и 1,8x20 мм относится соответственно к 1,2, 3,4 или 5 фракции.

В результате обработки экспериментальных данных методами математической статистики получены эмпирические зависимости среднего усилия ^ разрушения ячменя для различных фракций при влажности зерна IV от 12 до 20 %. Уравнения регрессии

имеют вид:

- для 1 -ой фракции зерна Г, = 0,377И/; - 14,03 IV + 271,4, (I) достоверность аппроксимации Я2 = 0,997;

-для 2-ой фракции зерна А = 0,3370^- 14,04ИМ 261,7, (Я2 = 0,997), (2)

- для 3-ей фракции зерна Г, = 0^77 И* - 14,04^ 4- 251,3, (Я7 = 0,993); (3)

- для 4-ой фракции зерна /Г, = 0.263И7" - 11,85IV + 223,9, (Я2 = 0,997); (4)

- для 5-ой фракции зерна Г5 = 0.516И,} - 20.14Г+271,1, (/^ = 0,980) (5)

При моделировании физических процессов в рабочей зоне шелушителя используется коэффициент трения зерна по ситовой и абразивной поверхности Для определения коэффициента трения использована методика В.В. Гортинского. Результаты исследований обработаны на ПК, в результате получены эмпирические зависимости коэффициента трения ячменя / по абразивной и ситовой поверхности в зависимости от влажности IV продукта:

- для средней крупности шлифзерен абразива 80 мкм

/т = 3,5-10"3 IV2 - 0,081 IV + 0,872, (Я2 = 0,995); (6)

- для средней крупности шлифзерен абразива 100 мкм

/,оо -= 3,6-10"3 т - 0,083 IV + 0,893. (Я3 = 0,996); (7)

- для средней крупности шлифзерен абразива 125 мкм

/т = 4,3-Ю1 ^-0,102^+ 1,03, (^ = 0,996); (8)

- для решетного полотна № 20 /т = 4,3-10"3 IV1 -0.102Ж+ 1,03, (Я2 = 0,996). (9)

Третья глава посвящена математическому моделированию рабочих процессов в шелушителе горизонтального типа с абразивным ротором. Полученная математиче-

екая модель позволяет аналитическим путем определять оптимальные режимы движения рабочих органов, производительность и мощность привода шелушителя.

Движение потока зерновок в рабочей зоне шелушителя может быть приближенно моделировано напорным винтовым течением вязкой несжимаемой жидкости внутри цилиндра радиусом 11 При этом осевая составляющая и. местной скорости потока или скорости частицы (зерновки) обусловлена перепадом давления с1р/сЬ потока на входе и выходе из канала, а окружная к, и радиальная скорости иг составляющие местной скорости потока - совместным движением потока в полости между корпусом шелушителя и абразивными дисками Если /?/- радиус абразивного диска, то в целях упрощения постановки задачи можно предполагать, что радиальная и окружная скорости потока индуцируются цилиндром некоторого приведенного радиуса, например а~ ■ /?, . Тогда приходим к модели течения смеси в полости между двумя цилиндрами радиусами 11 и а, когда внешний цилиндр радиуса ¡1 (корпус шелушителя) неподвижен, а внутренний, моделирующий абразивный диск радиуса а, движется с заданной угловой скоростью со, и, кроме того, с некоторой осевой скоростью ы, обусловленной величиной производительности шелушителя П С целью математической формапиза-

I >Ъ '

ции кинетики потока зерносмеси в шелушителе отнесем анализируемый объем зер-

* J г > 1

носмеси к циликдрической системе координат г, <р, : которая жестко связана с ротором, имеющим вращающийся вал и пакет насаженных на него абразивных дисков, на-

1 , I' I

правляя оси координат так, как показано на рис. 2.

1С

Рис 2 Схема к расчету поля скоростей частиц в рабочей зоне шелушителя

Вводя цилиндрическую систему координаг (г, ц>, г) записываем систему уравнений установившегося осесиммегричного движения жидкости, моделирующей сыпучую смесь:

- •^к"'-^]; ос)

г р с г V Г - )

и и

г

р

л "" - -

г"

¿¡и. ; (12)

р о:

+ о , (13)

дг дг

где

д с 4 д2 д д2 1л=иг—\-и,—; Д, =—г + — + —-дг ' дг дг гдг дг'

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

В таком случае в качестве граничных условий могут быть приближенно приняты условия прилипания смеси к стенке канала:

мЛ=м?)=мг= 0, при г = Я (14)

иг= 0, = V, «г=1с при г = а, (15)

где V = охг - окружая скорость внутреннего цилиндра, моделирующий абразивный рабочий орган; ш - угловая скорость абразивного ротора.

Предполагая, что расход смеси (производительность) Я есть заданная и постоянная величина, естественно считать, что ди, /дг = 0. Тогда вследствие уравнения неразрывности (13) и граничных условий (14). (15) получим иг -0, т.е. в рамках принятых допущений поток смеси в поперечном сечении следует считать квазикруговым. Кроме того, рассматривая поток вдали от питающею канала (бункера), полагаем, что и ди№ 'дг - 0. Тогда, дифференцируя обе части уравнения (10) по г, получаем:

дгр!(дгдг) = (д I дг)(др! дг) = 0 , откуда др/дг = %(г). Иначе говоря, в осевом направлении давление изменяется так же, как и в случае течения Пуазейля. Вследствие уравнений (11), (12) и граничных условий (14), (15) получаем:

и(, = КЯ(г2-а2)/(К2-«г)/г , (16)

иг =С(г2~Л2) + С,1п(г/Л). (17)

На значение осевой составляющей скорости и2 накладывается условие постоянства расхода смеси

П = 2п\и/Л'. (18)

Тогда, согласуя (16) с (17) и (18), находим

С = (2II- ч>тту) /[л-(Я2 - а2 )(а2 - Я2 + г)],

Г = [а2 - К2 + 2а2 1п(Л/а)]/ 1п(о/ /?)

Производигельность шелушителя:

кк _

// = 2тг jiijdr = Ink ji R - r)rdr = k(R - a)\R - 2a)

ii и /

¡f=*uiR-a)(R + 2a). (19>

3

В реальных условиях учитываться коэффициент заполнения рабочей зоны шелушителя т:

n = *MR-aXR~2a). (20)

Давление потока:

/>= j[('42-«2H^-a2)f/'-V'- =

= [к/г/( /гэ - a2 )f [ra / 2 - 2fl2 lit»— л4 /(2/-3)]* /»„(г). ' (2|)

Дифференциальные уравнения траектории частицы

dr-'ur rd<p'u,r ct u. или, учитывая, что и, 0, dr 0, г В const,

r2(R2-a~)d(p!fR(r2-а2)]-= ± /{с[г2 - Л2 +(R2 -a2) \n(r! R)/In (а/

откуда получаем

г аср, г Д (22)

где a = r['2 ~Rl "u'/?2-a2)\n(Rlr)!\n{c,l R)\2(R2 -я2 )/[к/?(г2 -<г)}

r = (23)

Момент всех сил вязкости, распределенных по кдкой-либо окружности радиуса г, относительно оси симметрии:

, 2г , , 4xAVRa2

м=) prrr-d<p = ~---(24)

„ К -<J

гдер,,,, - составляющая тензора напряжений:

Зт К

дг г к - а г Таким образом, момент сил вя ¡кости, распределенных по любой окружности, относительно оси симметрии не зависит от радиуса этой окружности, но зависит от динамической вязкости среды X, окружной скорости сгенки ротора V и радиусов цилиндров а и R С помощью формулы (24) может быть рассчитан крутящий момент, сообщаемый ротору цлелушитсля (длинной, L) со стороны потока обрабатываемой сыпучей смеси С учетом реальной работы шелушителя в формулу введен коэффициент заполнения зоны шелушения г.

М„ = --<Л = -------— I

о

Я2~а2

(26)

Мощность привода шелушителя с горизонтальным ротором рассчитывается исходя из крутящего момента (26)

N , = М „ со, =

4 жД УЯа 2 г

Л

2

(27)

Таким образом, в соответствии с формулами (24-27) мощность Л',„ в рамках принятых допущений, зависит от размеров шелушителя, частоты вращения ротора, а также от динамической (эффективной) вязкости зерносмеси.

Для того чтобы оценить порядок таких параметров как производительность и мощность привода условно полагаем, что И "=0,108-Ю, 115 м, а -0,086 м, р =698 кг/м1, т =0,75, м>=0,016-5-0,028 м/с, ^=1200+2000 об/мин, £=0,446 м. Динамическая вязкость зерносмеси X зависит от различных параметров сырья, и может находиться в пределах 1-20 ГШ:.

Проведенные на ПК расчеты выявили возрастающий характер зависимости производительности от приведенной осевой скорости потока продукта п1 (рис. 3). При этом, как видно из рисунка, для значений /?=0,108 м и /М), 115 м получены различающиеся кривые производительности. На рис. 4 представлена зависимость мощности Ыр привода, при различных скоростях вращения абразивного ротора, от параметра вязкости X продукта. Из рисунка видно, что затрачиваемая мощность зависит от скорости вращения абразивного ротора шелушителя и возрастает с увеличением вязкости продукта П кг/с

О 12

• /!' 0,108 м, а-0.085 м - Н-а 1 I? м '« = 0 085 м

к м/с

Рис. 3 График зависимости производительности шелушителя от приведенной осевой скорости потока продукта

• и,, - 1600 об/мин

• п, 2000 об/мин

Л, ПаЪ

Рис. 4. График зависимости мощности привода шелушителя ог динамической (эффективной) вязкости продукта

Для анализа кинетики зерновой смеси в рабочей зоне шелушителя, а также для исследования влияния дополнительных конструктивных элементов на работу шелушителя в работе использован метод компьютерного моделирования с помощью системы COSMOSFloWorks Данный программный продукт является разработкой немецкой фирмы NIK Л GmbH и предназначена для гидродинамического анализа. В COSMOSFloWorks движение и текучей среды моделируется с помощью уравнений Навье-Стокса, описывающих в нестационарной постановке законы сохранения массы, импульса и энергии этой среды. В COSMOSFloWorks используется численный метод -метод конечных объемов, вся расчетная область покрывается расчетной сеткой, значения независимых переменных рассчитываются в центрах ячеек.

В результате компьютерного моделирования исследовано влияние различных конструктивных изменений на кинетику зерновой смеси в рабочей зоне шелушителя. На рис. 5.а представлены результаты расчета флектории движения продукта при производительности шелушителя 400 кг/ч и частоте вращения абразивного ротора 2000 об/мин. Из рис. видно, что траектория имеет вид винтовой линии внутри цилиндрического корпуса шелушителя. При этом цвет линий характеризует изменение скорости продукта при перемещении На рис 5,6 представлена траектория продукта при установке в рабочей зоне двух направляющих под углом 30° против направления движения зерновой смеси В этом случае продукт отбрасывается в междисковое пространство,

I !

тем самым увеличивается интенсивность взаимодействия зерна с абразивной поверхностью Как показали результаты экспериментальных исследований, установка данных лопаток в реальных условиях приводит к повышению эффективности шелушения

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

При оценке интенсивности воздействия рабочих органов шелушителя на обрабатываемый продукт особый интерес представляет абсолютная скорость и общее давление продукта на рабочие поверхности абразивных дисков. На рис. 6,а представлена карта скоростей зернового потока в сечении шелушителя и на поверхности четырех абразивных дисков при частоте вращения ротора пр = 1200 об/мин По цвету заливки можно судить о скорости продукта в различных точках рабочей зоны шелушителя. При частоте вращения ротора пр = 1200 об/мин скорость продукта на периферии дисков находится в пределах 811 м/с. При частоте вращения ротора 2000 об/мин (рис. 6,6) скорость продукта на внешнем цилиндре абразивов максимальна, и достигает 15-20 м/с, а ближе к валу составляет 6-8 м/с Сравнивая результаты математического и компьютерного моделирования с экспериментальными исследованиями, и отталкиваясь от таких контролируемых параметров как производительность, время обработки продукта, эффективность шелушения и пр. можно с уверенностью говорить об адекватность полученных в результате математического моделирования результатах.

а)

0*

г

-АГ ,:К V: . ■ л ь, мм.« 1

*ЛГ ' V*(I; ЙгЦЙ "

У -Рк .- У л

ип, «V Ч. ^¡'^

б)

Чаг~

Рис. 6. Анализ поля скоростей в рабочей зоне шелушителя при различных технологических параметрах: а) п = 1200 об/мин; <2 = 400 кг/ч, б) п = 2000 об/мин; 0 = 400 кг/ч.

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

Исследовано влияние крупности шлифзерен ¡л на коэффициент шелушения Кш ячменя при частоте вращения ротора шелушителя 2000 об/мин (влажность зерна -12 %).

150 200 230 300 350 400

о о - средняя кр^пноаь зерна 65 мкм ---- полиномиальная per wo порядка Ц к[/ч

хх - средняя крупность крна 80 мкм --полиномиачьиая per Ч-го порядка

U (3- средняя крупность зерна 100 мкм - полиномиальная pet 3-го порядка

+ +- средняя крупность терна 125 мкм - - - полиномиальная per Vro порядка

Рис. 7. Зависимость коэффициента шелушения ячменя от п$а|рводительности шелушителя при средней крупности шлифзерен абразивов ~ 65, 80, 100 и 125 мкм

Результаты исследований статистически обработаны на ПК, в результате получены уравнения полиномиальной регрессии 3-го порядка.

- для абразивов со средней крупностью шлифзерен 65 мкм

Кш №s) = -1,6 10"5 П3 + 0,01Ш - 4,10П + 460,5 (/?•' - 0,983); (28)

- для абразивов со средней крупностью шлифзерен 80 мкм

К„П1» - -1,6 10 s П3 f 0,014/7г - 4,10П 1 467,4 (R2 = 0,998); (29)

- для абразивов со средней крупностью шлифзерен 80 мкм

К, „„an = -U 2-Ю'5 П3 +0,011/7' -3.41П + 420,2{Г?= 1); (30)

- для абразивов со средней крупностью шлифзерен 125 мкм

К,„ (,ъ, =-1,510П3 - 0,013П2 3,9П + 462,4 (R2 = 0,998) (31)

Проведена серия экспериментов по определению влияния различных конструктивных элементов в конструкции шелушите тя на такие показатели процесса шелушения как коэффициент шелушения (Кш) и производительность (ГГ) Па рис 8.а представлена зависимость коэффициента шелушения ячменя от количества установленных на роторе гонков (соответственно 0 - без гонков, с 1, 2 или 3 — мя гонками), при угле установки гонков в 30° и 45° и частоте вращения ротора 1600 и 2000 об/мин Влияние

гонков на производительность показано на рис. 8,6 Из гистограмм видно, что установка гонков приводит к повышению производительности шелушителя (с 310 до 420 кг/ч, при частоте вращения ротора 2000 об/мин), но вместе с тем существенно снижается коэффициент шелушения с 63 % до 40 % .

Кш, % П, кг/ч

- пг - 1600 об/мин О - пе -2000об/мин(<р = 30°) 0 -я,. - 1600 об/мин (ф = 30°) о -я, 2000 об/мин (Ф = 45") в - я, - 1600 об/мин (ф = 45°)

Рис 8. Зависимость коэффициента шелушения ячменя (а) и производительности (б) от количества и угла установки гонков при частоте вращения ротора 1600 и 2000 об/мин

Исследовано влияние влажность зерна на эффективность процесса шелушения С целью определения математических зависимостей в экспериментах по шелушению использовался ячмень влажностью 12, 16 и 20 %. На рис. 9 представлены основные результаты исследований. Эмпирические уравнения полиномиальной регрессии 3-го порядка имеют вид:

- влажность 12 %: Кш(12%1 -1,4 10"5 Г? 0,012Г12 - 3,68Я + 439,8 (Я"' = 0,990); (32) -влажность »И 6%: Кш,,,,,, -8,4 10"6 II3 - 8,4/0° П2 3,93/7 + 403,5 (1? = 0,999); (33)

- влажность И=20 %: К,„-2,5■■ 10 ' П3 + 0,019Г12 - 5,21Г14 517,5 (Л2 = 0,999); (34)

Кривая, аппроксимирующая значения Кш при влажности 20 % проходит значительно ниже кривых, соответствующих влажности 12 % и 16 %. Одной из причин снижения эффективности шелушения зерна ячменя можно считать «засаливание» шлифзерен абразивов, которое выявлено после визуального осмотра рабочих органов шелушителя.

Рис. 9. Изменение коэффициента шелушения в зависимости от производительности шелушителя при влажности ячменя 12, 16 и 20 %

Для решения задачи оптимизации процесса шелушения необходимо иметь математическую модель, описывающую зависимость коэффициента шелушения Кш и удельной мощности затрачиваемой на шелушение зерна, от производительности П и частоты вра.щения ротора п одновременно Для нахождения эмпирических зависимостей проведена ¿ерйя научных исследований на экспериментальном шелушителе Ис-следрэания проводились пру частоте вращения ротора п = 1200, 1600, 2000 и 2400 об/мин Производительность .изменялась в пределах от 120 до 450 кг/ч

Результаты исследований для коэффициента шелушения К,„ представлены на рис. 10. График получен путем двумерной аппроксимации полиномом 3-ей степени экспериментальных данных. Уравнение двухмерной регрессии имеет вид' Кю (п, П) =3,26-108 п3 -1,29-10'4п К),46л-1,03-10"5/7^8,67-10"77:-2,62/7-- 2,1210'8 п2П+ +7,83 10"'«/Т-3,94 10Л?Я+59,03. (35)

Результаты двумерной аппроксимации зависимости удельной мощности /V,,, от производительности шелушителя и частоты вращения ротора представлены на рис 11. Уравнение двумерной регрессии имеет,¡еид:

ПМУ(пЛ) - -3,26 10 "V -2,06- 10>'г5,8 10 «-8,31 ■ 10*7/'в,24-10"7/'-0,29/7+ , +4,2' 10'' п'П-1,23-10'7л/7;+8,24-10 'гс/7+60,51. (36)

Рис. 10. Зависимость коэффициента шелушения от производительности шелушителя и частоты вращения ротора ,¥,,,. кВтч/т

Рис. 11. Зависимость удельной мощности от производительности шелушителя и частоты вращения ротора

В пятой главе, с целью определения оптимальных конструкционных параметров шелушителя поставлена и решена задача нелинейного программирования. Сформулирована целевая функция оптимизации, найдены оптимальные значения конструктивных параметров шелушителя. Разработан алгоритм и написана программа для ПК «Шелушитель», с помощью которой можно рассчитать оптимальные параметры при проектировании шелушителя.

Параметрическая схема процесса шелушения представлена на рис. 12. Основными выходными характеристиками процесса являются удельные затраты энергии производительность П и коэффициент шелушения Кш. Основными входными характе-

19

ристиками, влияющими на эти показатели, являются характеристики продукта влажность IV, его гидротермическая обработка (ПО) и нагрузка (¿\ констр>ктивныс параметры' частота вращения ротора п, номер абразива //, количество и угол установки гонков и лопаток к ¡. установка задвижки к„ которая позволяет увеличить время обработки продукта, а также технологическая схема процесса шелушения к(-\

Экспериментальные исследования пока?али. что наибольшее влияние на процесс шелушение оказывает производительность машины и частота вращения ротора Введем ряд упрощающих топушений и будем считать, что коэффициент шелушения Кш и удельный расход электроэнергии /V,,, затрачиваемый на шелушение, зависят в основном от производительности П и частоты вращения ротора шелушителя и

У

I 5

г 5

Влажность черна ^ ''ъ

иаф\'1(.э О к.

ГТО

МШ1014(ЦНП1Ш11

шелуиштель с абрашиным ротором

ПроизводительностьП кг V

Тсхнот »ффективностъ К„

Удел (исход жергмиЛм

Рис. 12. Параметрическая схема процесса шелушения зерна

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

Я„ш,е>,П) - Кш(пЛ) Ш1,,-И„) ПЦ.'т3, (37)

где П - производительность шелушителя, кг/ч; Кш(п,Г1) - коэффициент шелушения, %; Ы>0(п,П) - удельный расход электроэнергии, кВт ч/т: //,„ - цена шелушенного зерна, руб/кг; 'Цо - цена не шелушенного зерна, руб/кг; //,, цена электроэнергии, руб/кВт-ч.

Определение параметров п и П, соответствующих максимуму функции (37) осуществляется на ПК с использованием программы \1athCAD численным методом. При решении задачи оптимизации задается область допустимых решений (ограничения) В

результате решения определена оптимальная частота вращения абразивного рогора и производительность шелушителя данного типа.

По результатам аналитических и экспериментальных исследований была разработана конструкторская документация, изготовлен и испытан на ЗАО «Ремонтно-мехакический завод» (г. Челябинск) опытный образец шелушителя горизонтального типа.

Для сравнения результатов, полученных в ходе испытания экспериментального шелушителя использовался промышленный шелушитель А1-ЗШН-Э. При проведении промышленных испытаний было установлено, что при однократном пропуске через малогабаритный шелушитель такой культуры, как ячмень, коэффициент шелушения по сравнению с машиной А1-ЗШН-Э повысился в среднем на 5 %, а удельный расход энергии снизился в среднем на 20 %.

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

В результате внедрения шелушителя данной конструкции на частных и фермерских хозяйствах затраты на процесс шелушения ячменя по сравнению с А1-ЗЦ1Н-3 снизились в среднем на 12 %.

Результаты исследований с 1999 г. используются в учебном процессе при изучении дисциплины '(Технологическое оборудование предприятий хлебопродуктов» для студентов специальностей 1706 и 2301

В приложениях представлены основные результаты экспериментов, листинг расчетных модулей для системы МаЛСАО, листинг программы для ПК «Шелушитель», акты об использовании результатов исследований в учебном процессе и на производстве.

ОСНОВНЫЕ ВЫВОДЫ И РЕЗУЛЬТАТЫ РАБОТЫ

В результате анализа научно-технической литературы, проведенных теоретических и экспериментальных исследований сделаны следующие выводы: 1 Исследование пяти фракций различной крупности ячменя показали, что: - для разрушения крупных зерен требуется большее усилие, чем мелких;

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

2 Аналитически и численным путем получена математическая модель рабочих процессов шелушителя с горизонтальным ротором В частности, она позволяет рассчи-

I

тать производительность шелушителя, крутящий момент и мощность привода шелушителя.

3. Показана возможность использования современных компьютерных САЕ - систем для моделирования кинетики сыпучей зерновой среды в шелушителе с вращающимся абразивным ротором. В результате этого получены кривые траектории движения продукта в рабочей зоне шелушителя (линии тока), определены значения скоростей и давления в зерновой смеси для различных точек рабочей зоны шелушителя при различных конструктивных и технологических параметрах шелушителя.

4. Исследовано влияние крупности шлифзерен на коэффициент шелушения ячменя. Установлено, что при увеличении крупности шлифзерен с 63 мкм до 125 мкм коэффициент шелушения увеличивается, но в месте с тем возрастает и относительное содержание эндосперма в проходовой фракции продуктов шелушения. В результате статистической обработки результатов исследований на ПК получено уравнение двумерной регрессии зависимости коэффициента шелушения К„, от крупности шлифзерен // и производительности /7

5. Установлено, что при использовании абразивных кругов с крупностью зерен 125 мкм в конце рабочей зоны шелушителя имеет место повышенное выкрашивание эндосперма. Для стабилизации процесса шелушения по длине рабочей зоны предложено использовать пять первых абразивных кругов с крупностью шлифзерен 125 мкм, а трех последующих с крупностью зерен абразива - 100 мкм.

6 Изучен эффект от установки различных дополнительных конструктивных элементов на качественные показатели работы шелушителя Результаты экспериментов пока-запй, что установка двух' пластин, возвращающих зерновой поток в рабочую зону, позволяет значительно увеличить коэффициент шелушения (на 14 %), при снижении производительности на 8 %.

7. В результате экспериментальных исследований подтверждены основные закономерности аналитических исследований и уточнена математическая модель процесса шелушения В результате статистической обработки результатов экспериментов на ПК с использованием программы сМа^ОАО» получены уравнения двумерной регрессии для определения коэффициента шелушения и удельного расхода энергии при

работе шелушителя в зависимости от частоты вращения ротора п и производительности шелушителя П, с учетом влажности исходного продукта.

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

9 Разработан алгоритм и программное обеспечение для поиска оптимальных конструктивных параметров работы шелушителя. В результате решения данной задачи на ПК, в системе MathCAD были найдены оптимальные значения производительности и частоты вращения ротора шелушителя.

СПИСОК ОПУБЛИКОВАННЫХ РАБОТ ПО ТЕМЕ ДИССЕРТАЦИИ

1 Таранин С.А., Глебов Л.А., Огурцов Ю.М. Исследование процессов шелушения ячменя на машине с горизонтальным расположением ротора. Тезисы доклада. Международная научно-теоретическая конференция «Молодые ученые - пищевым и перерабатывающим отраслям АПК», г. Москва, МГУПП, 1997.

2. Таранин С.А., Глебов Л.А , Огурцов Ю.М. Малогабаритный шелушитель горизонтального типа для ячменя. Тезисы доклада. Международная научно-техническая конференция «Ресурсосберегающие технологии пищевых производств», г. Санкт-Петербург, СПГАХиПТ. 1998.

3. Таранин С.А. Усовершенствование конструкции малогабаритного шелушителя для ячменя Тезисы доклада. Международная научно-практическая конференция «Индустрия продуктов здорового питания - третье тысячелетие (человек, наука, техника, экология)», г. Москва, МГУПП, 1999.

4. Таранин С.А Эффективность использования шелушенного ячменя при производстве комбикормов. Тезисы доклада. Научно-техническая конференция «Молодые ученые - пищевым и перерабатывающим отраслям АПК (технологические аспекты производства)», - М , МГУПП, 1999.

5. Глебов Л.А., Таранин С.А. Моделирование процесса движения отдельной зерновки в непосредственном соприкосновении с поверхностью абразивного круга шелушителя для зерновых культур с горизонтальным расположением ротора. - М.: Хлебпродинформ, информационный сборник «Научно-технические достижения и передовой опыт в отрасли хлебопродуктов», 2000,- вып. 4-5.- с. 3-6.

6. Таранин С.А. Эффективность использования малогабаритного шелушителя для зерновых культур в фермерских и частных хозяйствах. - М.: Хлебпродинформ, информационный сборник «Научно-технические достижения и передовой опыт в отрасли хлебопродуктов», 2000.- вып 4-5.- с 7-8.

7. Семенов Е.В., Глебов Л.А., Таранин С.А. О движении сыпучей смеси в рабочей зоне шелушителя. - М.: Хранение и переработка сельхозсырья, 2000,- №10,- с. 69-71.

8. Таранин С.А. Яблоков А.Е Исследование малогабаритного шелушителя. - М.: Объединенный научный журнал, 2005. №20. -с. 69-70.

9. Таранин С А. Яблоков А.Е. Оптимизация процесса шелушения ячменя в малогабаритном шелушителе. - М.: Объединенный научный журнал, 2005. №22. -с. 87-89.

I Подписано в печать 16.11.05

Формат 30x4^ Бумага типографская № 1 Печать офсетная. Уч.-изд. л. 0",9. Печ. л. 1,1. Тираж 120 экз. Заказ 310. г 125080, Москва, Волоколамское ш.. 11 Издательский комплекс МГУГ1П

*

í

í

i I

I

i,

}

I

I ]

»

I

t I

T é

\ \

Р24 186

РНБ Русский фонд

2006-4

26650

Введение.

Глава 1. Обзор современного состояния теории и практики шелушения зерновых культур. Задачи исследования.

1.1 Основные сведения о физико-механических свойствах зерна.

1.2 Технологические свойства зерна, влияющие на процесс шелушения.

1.3 Влияние гидротермической обработки (ГТО) зерна крупяных культур на процесс шелушения.

1.4. Классификация методов шелушения зерновых культур.

1.5. Обзор технологического оборудования для шелушения зерновых культур.

1.5.1. Шелушение зерна сжатием и сдвигом.

1.5.2. Шелушение зерна однократным и многократным ударом.

Ч 1.5.3. Шелушение зерна постепенным снятием оболочек в результате их интенсивного истирания.

1.6. Оценка технологической эффективности шелушения.

11 1.7. Краткие выводы, цели и последовательность исследований.

Глава 2. Экспериментальные установки и методология исследований. Исследование физико-механических свойств исходного сырья.

2.1. Описание экспериментальных установок, используемых в исследованиях.

• 2.1.1. Установка для определения прочностных характеристик зерна.

2.1.2. Установка для определения коэффициентов внешнего и внутреннего трения зерна.

2.1.3. Экспериментальная модель абразивного шелушителя с горизонтальным расположением ротора для зерновых культур.

2.2. Методология экспериментальных исследований.

2.2.1. Методика определения среднего усилия разрушения частиц зерна при сжатии.

2.2.2. Методика определения коэффициентов внешнего и внутреннего трения зерна.

2.2.3. Определение числа повторностей проводимых экспериментов и статистическая обработка результатов экспериментальных исследований

2.2.4. Методика определения коэффициента шелушения зерна.

2.2.5. Методика определения производительности шелушителя.

2.2.6. Методика определения энергозатрат.

2.2.7. Стандартные методики.

2.3. Характеристики ячменя, используемого в качестве сырья при проведении экспериментов.

2.3.1. Определение прочностных характеристик исходного ячменя.

2.3.2. Определение коэффициента трения ячменя по абразивной поверхности.

2.4. Выводы по главе.

Глава 3. Аналитические исследования и разработка математических моделей рабочих процессов в шелушителе с горизонтальным расположением ротора.

3.1. Математическое моделирование рабочих процессов в шелушителе с горизонтальным абразивным ротором.

3.2. Анализ результатов математического моделирования.

3.3. Компьютероное моделирование и исследование кинетики зерновой смеси в рабочей зоне шелушителя и использованием CAE (Computer Aided

Engineering) - системы COSMOSFlo Works.

3.3.1. Анализ влияния дополнительных конструктивных элементов на кинетику зерновой смеси.

3.3.2. Анализ поля скоростей и давления в рабочей зоне шелушителя при различных технологических параметрах.

3.4. Выводы.

Глава 4. Экспериментальные исследования по определению влияния конструктивных параметров шелушителя на качественные показатели его работы.

4.1. Исследование влияния крупности зерна абразивных кругов и производительности шелушителя на эффективность шелушения ячменя.

4.2. Исследование и оптимизация эффективности процесса шелушения ячменя по длине рабочей зоны шелушителя.

4.3. Исследование влияния длительности обработки ячменя на коэффициент шелушения.

4.4. Исследование влияния количества пропусков зерна через шелушитель на эффективность шелушения.

4.5. Изучение влияния дополнительных элементов, вносимых в конструкцию шелушителя, на эффективность шелушения ячменя и производительность шелушителя.

4.5.1. Влияние гонков на эффективность работы шелушителя.

4.5.2. Влияние элементов, возвращающих поток зерна в рабочую зону шелушителя на эффективность шелушения.

4.5.3. Влияние элементов, задерживающих зерновой поток в рабочей зоне шелушителя на эффективность шелушения ячменя. ИЗ

4.6. Исследование влияния частоты вращения ротора и производительности на эффективность работы шелушителя.

4.7. Исследование влияния влажности ячменя и производительности на эффективность работы шелушителя.

4.8. Определение удельного расхода электроэнергии на процесс шелушения ячменя.

4.9. Выводы по главе.

Глава 5. Оптимизация конструктивных параметров малогабаритного шелушителя и практическое внедрение результатов исследований.

5.1. Постановка задачи оптимизации конструктивных параметров малогабаритного шелушителя абразивного типа.

5.2. Разработка алгоритма и программного обеспечения для решения задачи поиска оптимальных параметров работы шелушителя.

5.3. Пример решения задачи оптимизации.

5.4. Производственные испытания малогабаритного шелушителя горизонтального типа.

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

Как известно, одним из основных показателей качества комбикормов для птицы является обменная энергия. Для кукурузы этот показатель равен 328-330 ккал на 1 кг, а для ячменя - 305 ккал. Ограничением для использования ячменя служит значительное содержание клетчатки. Это не позволяет вводить в комбикорма не шелушенный ячмень в количестве, достаточном для того, чтобы выдержать показатель по обменной энергии.

Шелушение ячменя, то есть отделение лузги (пленчатость ячменя варьирует в пределах от 8 до 15 %), приводит к снижению содержания клетчатки до 2,2 %, т.е. в 2,5 раза, что позволяет увеличить количество вводимого в комбикорма ячменя. Это способствует резкому снижению потребности в кукурузе при сохранении требуемой обменной энергии. Более того, общее содержание сырой клетчатки в комбикормах для птицы с шелушенным ячменем меньше, чем в комбикормах с кукурузой.

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

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

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

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

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

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

Для достижения указанной цели были поставлены и решены следующие задачи:

1. Изучить технологические свойства зерна, влияющие на эффективность шелушения;

2. Разработать математические модели рабочих процессов для малогабаритного шелушителя с горизонтальным ротором;

3. Создать экспериментальную установку на базе малогабаритного шелушителя для проведения натурных экспериментов;

4. Создать математические модели процесса шелушения зерна;

5. Провести экспериментальную проверку математических моделей и аналитических исследований шелушителя;

6. Сформулировать целевую функцию оптимальности процесса шелушения, поставить и решить задачу оптимальности процесса шелушения;

7. Разработать алгоритм расчета оптимальных параметров работы малогабаритного шелушителя, создать программу для расчета шелушителя на ПК;

8. Усовершенствовать конструкцию шелушителя и провести промышленную апробацию результатов исследований.

Методы исследования. В работе использованы методы математического и натурного моделирования, метод нелинейного программирования, методы математической статистики. Математическая обработка результатов исследований, а также решение задачи оптимизации процесса шелушения проводились на ПК с использованием программы компьютерной математики MathCAD. При геометрическом моделировании шелушителя использовалась CAD - система SolidWorks, а для гидродинамического анализа - система COSMOSFloWorks.

Научная новизна исследований заключается в следующем:

- в результате аналитических исследований получена и экспериментально проверена математическая модель рабочих процессов шелушителя с горизонтальным ротором;

- методами статистического анализа определено влияние технологических свойств зерна на количественные и качественные показатели процесса шелушения, что позволило в дальнейшем решить задачу оптимизации процесса шелушения;

- показана возможность использования компьютерной САЕ — системы С08М08Р1о\УогкБ для моделирования кинетики зерновой смеси в рабочей зоне шелушителя;

- в результате статистического анализа экспериментальных данных получены эмпирические формулы зависимости эффективности шелушения ячменя от различных конструктивных и технологических параметров работы шелушителя, которые так же позволили решить задачу оптимизации процесса шелушения;

- доказана возможность шелушения ячменя для кормовых целей в зоне свободного воздействия на него абразивных поверхностей;

- сформулирована целевая функция оптимальности процесса шелушения и методами нелинейного программирования на ПК определены оптимальные конструктивные параметры и рациональные режимы работы шелушителя для ячменя.

Объектом исследований является технологический процесс шелушения зерна, его режимы и закономерности.

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

Практическая значимость исследований:

1. Результаты теоретических и экспериментальных исследований получили практическую реализацию в совершенствовании оборудования для шелушения ячменя с целью оптимизации конструкции и ресурсосбережения;

2. Разработаны математические модели, алгоритмы и компьютерная программа «Шелушитель» для определения оптимальных технологических параметров работы малогабаритного шелушителя горизонтального типа.

3. Показана возможность использования современных компьютерных САЕ -систем для моделирования кинетики сыпучей зерновой среды в шелушителе с вращающимся абразивным ротором. В результате компьютерного моделирования получены значения различных физических параметров зерновой смеси в рабочей зоне шелушителя, проанализировано влияние конструктивных изменений на кинетику продукта.

4. Разработана конструкторская документация на малогабаритный шелушитель горизонтального типа, изготовлен опытный образец шелушителя, который установлен на ЗАО «Ремонтно-механический завод». Данные промышленных испытаний подтвердили высокие технико-экономические показатели малогабаритного шелушителя.

Апробация работы. Основные положения диссертационной работы докладывались на Международной научно - практической конференции «Индустрия продуктов здорового питания - третье тысячелетие» (Москва, МГУПП, 1999 г.); на расширенном заседании кафедры «Технологическое оборудование предприятий хлебопродуктов» (Москва, МГУПП, 2005 г.).

На защиту выносятся: математическая модель шелушителя горизонтального типа с горизонтальным абразивным ротором; эмпирические зависимости влияния различных конструктивных параметров шелушителя на эффективности и энергоемкость процесса шелушения; результаты исследований физико-механических свойств исходного продукта — зерна ячменя; программные средства, разработанные на основе математической модели шелушителя для расчета оптимальных конструкционных параметров машины; многофакторные зависимости изменения степени шелушения ячменя от конструкционных параметров шелушителя, а также от технологических параметров исходного продукта;

- целевая функция оптимизации процесса шелушения по технико-экономическим показателям работы шелушителя; алгоритм и программа для определения оптимальных конструктивных и технологических параметров шелушителя горизонтального типа.

Публикации. По материалам исследований опубликовано 9 печатных работ.

Объем работы. Диссертационная работа состоит из введения, 5 глав,

4.9. Выводы по главе

1. Исследовано влияние крупности шлифзерен на коэффициент шелушения ячменя. Установлено, что при увеличении крупности шлифзерен с 65 мкм до 125 мкм коэффициент шелушения увеличивается, но в месте с тем возрастает и относительное содержание эндосперма в проходовой фракции продуктов шелушения. В результате статистической обработки результатов исследований на ПК получено уравнение двумерной регрессии зависимости коэффициента шелушения Кш от крупности шлифзерен /л и производительности П:

Кш (И,П) =-7,25-10"4 п3 +0,23л2 -20,9я+4,53 • 10'6П3 - 7,23-10"3Я2-1,92Я

- 6,73-Ю"5 п2П+3,48-10"5я772-5,08-10"3лЯ+593,13.

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

3. Исследована возможность повышения эффективности шелушения путем задержки зерна в рабочей зоне шелушителя на 30 — 90 сек. Эксперименты показали, что для достижения требуемого коэффициента шелушения в 85 % необходима задержка продукта на 60 сек. Но это приводит к снижению производительности шелушителя в 2,5 раза (с 300 кг/ч до 120 кг/ч), следовательно, таким образом повышать коэффициент шелушения не целесообразно.

4. Изучен эффект от установки различных дополнительных конструктивных элементов на качественные показатели работы шелушителя. Результаты экспериментов показали, что в целом, использование дополнительных конструктивных элементов в шелушителе не эффективно, т.к. хотя это и приводит к повышению коэффициента шелушения, но вместе с тем снижается производительность. Однако, установка двух пластин, возвращающих зерновой поток в рабочую зону позволяет значительно увеличить коэффициент шелушения (на 14 %), при снижении производительности на 25 кг/ч.

5. С целью получения математических моделей процесса шелушения проведена серия экспериментов. В результате статистической обработки результатов исследований на ПК с использованием программы «МаШСАО» получены следующие уравнения двумерной регрессии:

- для зависимости коэффициента шелушения от частоты вращения ротора п и производительности шелушителя Я:

Кш (п,П) =3,264-Ю-8 п3 -1,984-10*4и2+0,465и-1,03• 105Я3+8,67' 103Я2-2,62Я

- 2,12'Ш8 и2Я+7,83- 10"7лЯ2-3,94-юЛз/7+59,03;

- для зависимости коэффициента шелушения от влажности зерна ¡V и производительности Я шелушителя:

КШ(Ж,П) =0,033 ¡¥3-2,52 0^+53,28 ^-7,84-10^+7,24-10-3Я2-2,ЗЗЯ

- 8,76-Ю-5 1,26-10"5 ЖП2-0,01 ЦГП+10,47;

- для зависимости удельного расхода энергии ИУд одновременно от двух параметров — частоты вращения ротора п и производительности Я:

Ыуд (п,П) = -3,26-Ю'10 п -2,06- 10^-5,8-10"3л+8,31 • 10^+3,24-10'4Л2-0,29П+ +4,2-10-9 п2П-1,23-1 0"7иЯ2+8,24- 10'5пП+60,51.

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

Глава 5. Оптимизация конструктивных параметров малогабаритного шелушителя и практическое внедрение результатов исследований 5.1. Постановка задачи оптимизации конструктивных параметров малогабаритного шелушителя абразивного типа

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

Параметрическая схема процесса шелушения представлена на рис. 56. Основными выходными характеристиками процесса являются удельные затраты энергии Муд, производительность П и коэффициент шелушения Кш. Основными входными характеристиками, влияющими на эти показатели, являются характеристики продукта: влажность IV, его гидротермическая обработка (ГТО) и нагрузка <2; конструктивные параметры: частота вращения ротора п, номер абразива ¡л, количество и угол установки гонков кг и лопаток кл, установка задвижки кз, которая позволяет увеличить время обработки продукта, а также технологическая схема процесса шелушения кСх

Экспериментальные исследования показали, что наибольшее влияние на процесс шелушение оказывает производительность машины и частота вращения ротора. В связи с этим введем ряд упрощающих допущений. В дальнейшем будем считать, что коэффициент шелушения Кш и удельный расход электроэнергии Л^, затрачиваемой на шелушение, зависят только от производительности П и частоты вращения ротора шелушителя п: Кш (п,П) =3,264-Ю*8 п3 -1,984-10*4«2+0,465«-1,03• 10'3П3+8,67-10'3П2-2,62П-- 2,12-10'8 п2П+ 7,83• 10'7пП2-3,94-10"4«Я+59,03; п,П) = -3,26-КГ10 п3 -2,06- 10"бл2-5,8-10"3л+8,31 • 10"8773+3,24-10АП2-0,29П+ +4,2' Ю'9 п2П-\,1Ъ- 10"7«Я2+8,24- 10-5«Я+60,51.

Производительность П. кг/ч -^

Влажность зерна IV, %

Малогабаритный шелушитель с абразивным ротором

Нагрузка 0, кг/ч

Технол. эффективность К,» %

ГТО

Удел, расход энергии Л^, -►

Рис. 56. Параметрическая схема процесса шелушения зерна

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

Ятах(пД) = Кш (п,П) • П-(Цш-Ц0)-Щд(п,П) -П-Цэл-10'3, (5.1) где П— производительность шелушителя, т/ч; Кш(п,П) — коэффициент шелушения, %; ЫУд(п,П) — удельный расход электроэнергии, кВт-ч/т; Цш - цена шелушенного зерна, руб/т; Ц0 - цена не шелушенного зерна, руб/т; Цэл - цена электроэнергии, руб/кВт-ч.

Определение параметров п и П, соответствующих максимуму функции (5.1) осуществляется на ПК с использованием программы МаШСАБ численным методом. При решении задачи оптимизации необходимо задавать область допустимых решений (ограничения).

5.2. Разработка алгоритма и программного обеспечения для решения задачи поиска оптимальных параметров работы шелушителя

Решение данной нелинейной задачи оптимизации возможно лишь численными методами, с использованием специальных прикладных задач. Наиболее удобным программной оболочкой, позволяющей решить задачи подобного типа, является система компьютерной математики МаШСАБ. Однако прежде чем приступить к написанию специализированного программного модуля (рабочего листа) для расчета необходимо разработать алгоритм работы программы.

Структурная схема алгоритма нахождения оптимальных конструктивных параметров работы малогабаритного шелушителя представлена на рис. 57. Процедура поиска оптимальных решений включает в себя несколько этапов:

1. Ввод оператором исходных данных и диапазон поиска оптимальных решений;

2. Расчет теоретической производительности и коэффициента шелушения, удовлетворяющих исходным данным;

3. Сравнение расчетных значений с минимально допустимыми;

4. В случае не удовлетворения последнего условия выдается сообщение об ошибке, и возврат в начало программы для корректировки исходных данных;

5. Из базы данных загружаются математические модели шелушителя, при необходимости, предусмотрена возможность их корректировки;

6. Решается задача оптимизации численным методом;

7. Результаты вычислений выводятся на печать.

Рис. 57. Структурная схема алгоритма расчета оптимальных параметров шелушителя для зерна

Программа расчета оптимальных параметров шелушителя реализована в системе компьютерной математики MathCAD-2001. Рабочий лист программы представлен в прил. 8. Обработка программы проводилась на ПК с процессором Intel Pentium IV ЗГТц.

При решении задачи пользователю необходимо задать диапазон изменения переменных, входящих в расчетный блок: минимально и максимально допустимый коэффициент шелушения — Km in и Ктах\ минимальную и максимальную скорость вращения ротора шелушителя - nmin и птах, а также текущую стоимость тонны неошелушенного и шелушенного зерна (переменные Ц0 и Цш), стоимость кВт-ч электроэнергии - Цэл.

В расчетном блоке представлены математические модели параметров процесса шелушения: К(п,Р) — коэффициент шелушения, как функция от частоты вращения ротора (п) и производительности (Р) шелушителя; N(n,P) — удельные затраты электроэнергии, как функция от частоты вращения ротора (п) и производительности (Р) шелушителя. При необходимости, эти функции могут быть легко откорректированы. Критерий оптимальности (5.1) записан в виде функции R(n,P) (см. прил. 8).

Для поиска значений переменных п и Р, при которых функция R(n,P) имеет максимальное значение, в программе используется функция maximize (R,n,P). Эта функция используется в составе блока решения, открываемого директивой Given, и возвращает вектор неизвестных, при которых заданная функция имеет максимальное значение. Внутри блока располагаются ограничительные условия в виде неравенств. Перед блоком задаются начальные значения искомых переменных.

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

5.3. Пример решение задачи оптимизации

При решении задачи нахождения оптимальных параметров работы шелушите ля использованы следующие постоянные величины: Удельный расход электроэнергии Nyd = 8 — 25 кВт-ч/т; Крупность абразивных зерен ^ = 125 и 100 мкм; Влажность продукта W= 12 %; Цена не шелушенного зерна, Ц0 = 3,6 руб/кг; Цена шелушенного зерна, Цш = 3,64 руб/кг; Цена за электроэнергию, Цэл = 1,36 руб/кВт-ч. Рассмотрено три варианта работы шелушителя:

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

Коэффициент шелушения Кш— 10-90 %; Частота вращения ротора шелушителя п = 1200 — 2200 об/мин; Диапазон производительности шелушителя: П = 200 — 400 кг/ч. Результатом решения являются следующие значения.

- оптимальная частота вращения ротора шелушителя: п = 1799 об/мин;

- оптимальная производительность: П= 358 кг/ч; -оптимальный коэффициент шелушения: Кш = 46 %;

- удельные затраты электроэнергии: Nyd = 12,25 кВтч;

- значение целевой функции оптимизации: R(n,ü') = 683,9.

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

Коэффициент шелушения Кш = 58 - 65 %; Частота вращения ротора шелушителя п = 1200 - 2200 об/мин; Диапазон производительности шелушителя: П= 200 - 400 кг/ч. Результатом решения являются следующие значения.

- оптимальная частота вращения ротора шелушителя: п = 1762 об/мин;

- оптимальная производительность: П = 279,6 кг/ч;

- оптимальный коэффициент шелушения: Кш= 58 %;

- удельные затраты электроэнергии: Иуд = 15,18 кВтч;

- значение целевой функции оптимизации: Я(п,П) = 642,9.

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

Коэффициент шелушения Кш- 85 - 95 %;

Частота вращения ротора шелушителя п = 1200 - 2200 об/мин;

Диапазон производительности шелушителя: П = 200 - 400 кг/ч.

Результатом решения являются следующие значения.

- оптимальная частота вращения ротора шелушителя: п = 1940 об/мин;

- оптимальная производительность: 77 = 200 кг/ч;

- оптимальный коэффициент шелушения: Кш- 86 %;

- удельные затраты электроэнергии: Иуд = 20 кВтч;

- значение целевой функции оптимизации: Я(п,П) = 625,7.

5.4. Производственные испытания малогабаритного шелушителя горизонтального типа

По результатам аналитических и экспериментальных исследований была разработана конструкторская документация и изготовлен на ЗАО «Ремонтно-механический завод» (г. Челябинск) опытный образец малогабаритного шелушителя горизонтального типа для шелушения зерновых культур.

В 2001 г. опытный образец шелушителя был испытан да данном заводе. Для сравнения результатов, полученных в ходе испытания экспериментального шелушителя использовался промышленный шелушитель А1-31ПН-3.

По данным, предоставленными специалистами ЗАО «Ремонтно-механический завод» (прил. 9) можно засвидетельствовать следующую практическую значимость результатов исследований.

При проведении промышленных испытаний было установлено, что при однократном пропуске через малогабаритный шелушитель такой культуры, как ячмень, коэффициент шелушения по сравнению с машиной А1-3111Н-3 повысился в среднем на 5%.

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

Удельный расход энергии при применении малогабаритного шелушителя при обработке ячменя по сравнению с А1-311Ш-3 снизился в среднем на 20%.

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

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

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

В результате внедрения шелушителя данной конструкции на частных и фермерских хозяйствах затраты на процесс шелушения ячменя по сравнению с А1-3111Н-3 снизились в среднем на 12%.

Результаты исследований с 1999 г. используются в учебном процессе при изучении дисциплины «Технологическое оборудование предприятий хлебопродуктов» для студентов специальностей 1706 и 2301 (прил. 10).

Заключение и общие выводы

Целью проведенных исследований является повышение эффективности процесса шелушения ячменя (для дальнейшего использования в составе комбикормов) путем совершенствования и оптимизации конструктивно-технологических параметров шелушителя для малого хозяйства. В рамках поставленной цели определен круг актуальных задач и приведены их решения.

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

На основе теоретических и экспериментальных исследований, получены следующие новые результаты:

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

2. Экспериментально определены физико-механические свойства исходного сырья, исследованы прочностные характеристики исходного ячменя для фракций различной крупности. В результате статистической обработки результатов исследований получены эмпирические зависимости усилия разрушения от влажности продукта для пяти, различных по крупности, фракций ячменя. Установлено, что для разрушения более крупных зерен требуется большее усилие. С увеличением влажности продукта усилие разрушение снижается.

3. Определены зависимости коэффициента трения ячменя по абразивной поверхности с различной крупностью шлифзерен. Установлено, что крупность шлифзерен в пределах 63 - 125 мкм, при одной и той же влажности продукта, не оказывает большого влияния на коэффициент трения.

4. В результате аналитических исследований получены формулы:

- для определения производительности шелушителя с горизонтальным абразивным ротором;

- для нахождения момента внутренних сил вязкости распределенных по окружности, описываемой рабочую зону шелушителя;

- для определения крутящего момента и мощности привода шелушителя с горизонтальным абразивным ротором.

- установлено, что момент сил вязкости зависит от динамической вязкости среды, окружной скорости ротора, радиуса абразивного цилиндра и внутреннего радиуса шелушителя;

5. Показана возможность использования современных компьютерных САЕ -систем для моделирования кинетики сыпучей зерновой среды в шелушителе с вращающимся абразивным ротором. В результате компьютерного моделирования получены кривые траектории движения продукта в рабочей зоне шелушителя (линии тока), определены значения скоростей и давления в зерновой смеси для различных точек рабочей зоны шелушителя при различных конструкционных и технологических параметрах шелушителя

6. Установлено, что при использовании абразивных кругов с крупностью зерен 125 мкм в конце рабочей зоны шелушителя имеет место повышенное выкрашивание эндосперма. Для стабилизации процесса шелушения по длине рабочей зоны предложено использовать пять первых абразивных кругов с крупностью шлифзерен 125 мкм, а трех последующих с крупностью зерен абразива - 100 мкм.

7. Исследовано влияние крупности шлифзерен на коэффициент шелушения ячменя. Установлено, что при увеличении крупности шлифзерен с 65 мкм до 125 мкм коэффициент шелушения увеличивается, но в месте с тем возрастает и относительное содержание эндосперма в проходовой фракции продуктов шелушения. В результате статистической обработки результатов исследований на ПК получено уравнение двумерной регрессии зависимости коэффициента шелушения Кш от крупности шлифзерен ц и производительности П.

8. Исследована возможность повышения эффективности шелушения путем задержки зерна в рабочей зоне шелушителя на 30 - 90 сек. Эксперименты показали, что для достижения требуемого коэффициента шелушения (85 %) необходима задержка продукта на 60 сек. Но это приводит к снижению производительности шелушителя в 2,5 раза (с 300 кг/ч до 120 кг/ч), следовательно, таким образом повышать коэффициент шелушения не целесообразно.

9. Изучен эффект от установки различных дополнительных конструктивных элементов на качественные показатели работы шелушителя. Результаты экспериментов показали, что установка двух пластин, возвращающих зерновой поток в рабочую зону позволяет значительно увеличить коэффициент шелушения (на 14 %), при снижении производительности на 8 %.

10. По результатам экспериментальных исследований получены математические модели процесса шелушения. В результате статистической обработки результатов экспериментов на ПК с использованием программы «МаЛСАЕ)» получены следующие уравнения двумерной регрессии:

- для зависимости коэффициента шелушения от частоты вращения ротора п и производительности шелушителя Я:

- для зависимости коэффициента шелушения от влажности зерна IV и производительности Я шелушителя:

- для зависимости удельного расхода энергии Л^ одновременно от двух параметров - частоты вращения ротора п и производительности Я.

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

12. Разработан алгоритм и программное обеспечение для поиска оптимальных конструктивных параметров работы шелушителя. В результате решения данной задачи на ПК, в системе МаШСАГ) были найдены оптимальные значения производительности и частоты вращения ротора шелушителя.

1. Абрамов С.Ю. Влияние влажности и температуры зерна на эффективность его переработки: Дис. . канд. техн. наук —М., 1984. -137 с.

2. Авдусь П.Б. Определение качества зерна, муки и крупы. М.: Колос, 1976.-с. 150-167.

3. Адлер Ю.П. Введение в планирование эксперимента. — М.: Металлургия, 1969. -157 с.

4. Адлер Ю.П., Маркова Е.В. Планирование эксперимента при поиске оптимальных условий. — М.: Наука, 1971.-283 с.

5. Александров В.Г. О строении эндосперма зерновки злака // Ботанический журнал, 1939. -Т.24.

6. Алямовский A.A. и др. Компьютероне моделирование в инженерной практике / Алямовский A.A., Собачкин A.A., Одинцов Е.В., Харитонов А.И., Понамарев Н.Б. СПб.: БХВ-Петербург, 2005. - 800 с.

7. Аммосова A.M. Исследование технологического процесса шелушения хлопковых семян: Дис. . канд. техн. наук —М., 1963. —143 с.

8. Анискин В.И. Гигроскопические свойства зерна различных культур. М.: Колос, 1967. -86 с.

9. Афанасьев В.А. Исследование тепловой обработки ячменя с применением ИК-нагрева при производстве комбикормов: Дис. . канд. техн. наук М., 1979. -179 с.

10. Бабиченко JI.B. Исследование крахмала взорванных зерен с помощью сканирующего микроскопа // Известия ВУЗов СССР, Пищевая технология, 1973.-N6.

11. И. Бабиченко JI.B. Изменение ультраструктуры крахмального зерна гречихи при тепловой обработке // Известия ВУЗов СССР, Пищевая технология, 1978. -N2.

12. Бардышев Г.М. Справочник мукомола, крупянщика, комбикормщика. -М.: Колос, 1973. -225 с.

13. Баум А.Е., Юкиш А.Е. Прогрессивная технология хлебоприемных и зерноперерабатывающих предприятий. М.: Колос, 1978. -с. 67-73.

14. Баум А.Е. Производство, хранение и переработка зерна в Швеции. М.: ЦНИИТЭИ Минзага СССР, 1979. -с. 54-65.

15. Белиловская A.C., Зенкова А.Н. Об основных признаках и показателях физико-механических свойств крупяного ячменя // Труды ВНИИЗ, 1954.-Вып. 53. -с. 129-144.

16. Белиловская A.C. Изменение структурно-механических свойств зерна в результате ГТО // ВНИИЗ, Тезисы докладов, 1960. -с. 34-40.

17. Берестов А.П. Совершенствование шелушения и ГТО овса: Дис. . канд. техн. наук М., 1987. -179 с.

18. Беркутова Н.С. Некоторые особенности микроструктуры зерна пшеницы // Мукомольно-элеваторная промышленность, 1970. -N8. -с. 19-20.

19. Беркутова Н.С., Швецова И.А. Микроструктура пшеницы. М.: Колос, 1977. -187 с.

20. Браслин С.Н. Совершенствование технологии выделения ядра с целью выработки гречневой крупы улучшенного качества: Дис. . канд. техн. наук-М., 1983.-169 с.

21. Братухин А.М., Суровегин И.А. Исследование процесса обработки зерна ржи в шелушильных машинах ЗШН // Труды ВНИИЗ, 1969. -Вып. 66. -с. 77-85.

22. Бродская И.А. Комбикормовая промышленность за рубежом // ЦНИИТЭИ Минхлебопродуктов СССР, Хранение и переработка зерна, экспресс-информация, 1989. -Вып. 7. —22 с.

23. Бродская Л.Г. Комбикормовая промышленность социалистических стран // ЦНИИТЭИ Минзага СССР, Комбикормовая промышленность, экспресс-информация, 1986. -Вып. 4. -16 с.

24. Бродская И.А., Паукова М.В. Комбикормовая промышленность за рубежом // ЦНИИТЭИ Минхлебопродуктов СССР, Хранение и переработка зерна, экспресс-информация, 1988. -Вып. 2. -24 с.

25. Бутко В.П. Исследование плотности зерна пшеницы как перспективного показателя оценки его технологических свойств: Дис. . канд. техн. наук -М., 1973.-145 с.

26. Бутковский В.А., Гафнер JI.A. Эксплуатация оборудования мельниц и крупозаводов.-М.: Колос, 1974.-е. 143-152.

27. Бутковский В.А. Технология мукомольного, крупяного и комбикормового производства. М.: Колос, 1975. - 237 с.

28. Бывайков М. Е. Исследование и разработка алгоритмов и программных методов для прогнозирования трендов параметров технических объектов. Автореф. дис. канд. техн. наук. -М., 1989. 17 с.

29. Вайнберг А. А. Надёжность оборудования предприятий по хранению и переработке зерна. Учебное пособие для вузов по специальности "Машины и аппараты пищевых производств". -Киев: Одесса: Высшая школа, 1986. -406 с.

30. Вайнштейн Б.Н. Абразивы и абразивные инструменты. Конспект лекций. Пермь: Перм. политехнический ин-т, 1976. —49 с.

31. Винников Т.И., Сенаторский Б.В. Основные направления развития и совершенствования комбикормового производства за рубежом // ЦНИИТЭИ Минзага СССР, Комбикормовая промышленность, экспресс-информация, 1978. -Вып. 7. -16 с.

32. Галицкий P.P. Оборудование зерноперерабатывающих предприятий. М.: Агропромиздат, 1990.-е. 136-148.

33. Галкин A.B. Указатель действующих государственных, отраслевых стандартов и технических условий в пищевой промышленности. М., 1976. -53 с.

34. Гафиер JI.A. Основы технологии приема, хранения и переработки зерна. -М.: Колос, 1975. -с. 281-296.

35. Георги Н.В. Новые виды технологического оборудования для переработки зерна. Одесский технол. ин-т пищевой промышленности им. М.В. Ломоносова. Киев, 1991. -116 с.36.