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

ИНЖЕНЕРНО-АГРОПРОМЫШЛЕННЫЕ

СПЕЦИАЛЬНОСТИ

УДК 631.362.001.573

Ю.И. ЕРМОЛЬЕВ, М.Ю. КОЧКИН МОДЕЛИРОВАНИЕ ПРОЦЕССА ФУНКЦИОНИРОВАНИЯ ЗЕРНООЧИСТИТЕЛЬНОГО АГРЕГАТА

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

Введение. Традиционная технология очистки продовольственного зерна в сельском хозяйстве предусматривает последовательный пропуск всего обрабатываемого зернового материала через комплекс зерноочистительных машин в агрегате. При наличии в обрабатываемом зерновом материале примесей, которые выделяются только на триерном блоке, возникают трудности (невысокая производительность триерного блока - 5-10т/ч), что создает диспропорцию между возможностями технологического оборудования в поточной линии. Производительность воздушно-решетных зерноочистительных машин (ВРМ) в среднем в два раза выше производительности последовательно работающих триерных блоков. Параллельное же объединение нескольких триерных блоков вызывает значительное увеличение габаритных размеров, стоимость агрегата, усложняет его конструкцию, энергоемкость, ремонт и, что особенно важно, может привести к росту приведенных затрат на очистку зернового материала. Изменение структуры агрегата, технологических свойств и подачи обрабатываемого материала в агрегат оказывает влияние на качество очищенного зерна.

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

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

Е фю = {Р, Л ^ [кш(х\ Тш(л)]}  тах;

Л1 с Д, Р с Р, х е (х, и);

в1 £ \в} I 8 з £ 8 з ] .

ВЗОА =

РЗОА = | Q,a . ,Ж,М

Рис.1. Структурная схема зерноочистительного агрегата

Здесь Р = РсМ - вектор входных воздействий на принятую в ЗОА систему операций, показанную на рис.1, где О. - подача зернового материала в ЗОА; а . >W - содержание в исходном материале j-х компонентов, их

влажность; М(Ь.) ,о 2 (Ь.) математические ожидания и дисперсии размеров признаков разделения j-х компонентов; f о (В) - плотность вероятности распределения подачи Q зернового материала по ширине В рабочих органов зерноочистительных машин в ЗОА;

А. - векторы управляющих факторов элементов системы, обеспечивающих ее функционирование (АВРМ - воздушно-решетной машины, АС -скельператора, АПС - пневмосепаратора, АРМ - решетного модуля, Атб - триерного блока), где Вш, ^ - ширина и глубина пневмосепараторов ВРМ; Vвi - рабочая скорость воздушного потока; плотность вероятностей распределения подачи fq(Bп) зернового материала и воздушного потока ЦВп) по ширине Вп пневмосепаратора; р* ( р ) - функциональные схемы

решетных модулей; плотности вероятности распределения зернового материала по решетным ярусам fQx(H) и по ширине решет fQx(В) в решетных модулях зерноочистительной машины ОЗС-50/25/10; К8м (х) - функциональная схема и параметры ВРМ;

Дс, Вс, Пс- диаметр, ширина и частота вращения скельператора; Г, Б, dп - ширина, длина отверстий и диаметр проволоки скельператора;

а, в, к, п- угол наклона, направленность, амплитуда и частота колебаний решет; Т, В^ Ьг тип решета, его длина, ширина, размеры отверстий; Км(Х)- функциональная схема решетного модуля;

Д, L,nтц - соответственно диаметр, длина и частота вращения триерных цилиндров, а - угол наклона приемного лотка, КВт(х) - функциональная схема триерных блоков;

X - набор технологических операций из Gi (х, и) множества, реализуемых в К8 м функциональной схеме ЗОА; GК 8 М ( Т, 8 МХ ( ) - математическая модель, определяющая показатели технологического процесса в ЗОА для принятой К8 м (Х) -й его схемы.

Выходные показатели функционирования воздушно-решетной машины (входные для триерных блоков) определяются вектором ВВРМ и

всего агрегата ВЗОА , независимые аргументы которого случайные в вероятностно-статистическом смысле величины (см. рис.1), где Еф - критерий эффективности реализации технологического процесса в воздушно-решетной машине или зерноочистительном агрегате [1]; еЬ]-т - полнота выделения из зернового материала отделяемых j-х компонентов, 8 з.,8 з - потери

зерна на хг операции и общие; Q0фQ „О- - масса выхода очищенного зерна, фуражных отходов и отходов разных категорий соответственно; масса Уф,,ух, >х и всех компонентов, содержащихся в зерновом материале после реализации хг -й операции; масса у',хг, ухг >х и всех компонентов, выделенных из зернового материала после реализации хг -й операции; содержание >го компонента в массе Ьп,хг очищаемого материала и

в выделяемых фракциях Ь,хг; полнота прохода е jxi >го компонента в

очищаемый зерновой материал после выполнения хг -й операции.

Обоснуем математическую модель ^ К^ мМХ) Т (X) ЗОА. В

соответствии с принципиальной схемой рассматриваемого варианта (см. рис.1), ЗОА содержит зерноочистительную машину ОЗС-50/25/10 и один или два триерных блока (ТБ) типа ЗАВ-10.90000А с известными функциональными связями между частными технологическими операциями.

Для этих условий полнота выхода е Ь, ;Ко компонента исходного

зернового материала в очищенную в агрегате фракцию зерна

Е

е Ь = П е Ь,м ,

I = 1

где еЬ,м - полнота выхода >го компонента, поступившего в соответствующую зерноочистительную машину (К=1(только ВРМ)), К=2(ВРМ+ТБ), К=3(ВРМ+ два ТБ)), - в зерно, очищенное рассматриваемой машиной. Полнота выхода е в отходы ;Ко компонента в ЗОА

Е

е ЬОТ, ~ X 0ЬОТ,М ' ЬЬОТ,М ! 0 ' а j ,

М = 1

где Оьог,-м - количество >го компонента зернового материала выделенного М-й машиной в отходы; ььот,м - содержание ;Ко компонен-

та в отходах М-й машины.

В машине ОЗС-50/25/10 последовательно реализуются четыре частные операции (см. рис.1). Для этих условий функционирования полнота выхода ;Ко компонента в очищенную машиной фракцию зерна

4

е = е

СЬ,1 СЬ,ОП

ОП = 1

Полнота выхода ;Ко компонента исходного зернового материала в сорные отходы (сход QСХС со скельператора, выход из отстойных камер

двух пневмосепараторов QЛ1 и Q Л2 аспирационной системы ОЗС-

50/25/10, сход Q СХР с решет решетного модуля )

е О, ~ ( 0СХС Ь,С + 0Л1 Ьj1 + 0Л 2 Ьj 2 + 0СХР Ь,р) 10 а,,

где Ь,с ; Ь,1; Ь,2 ; Ь,р - соответственно содержание >го компонента в сходовой фракции со скельператора, во фракциях, выделенных в легкие

отходы пневмосепараторами №1 и №2, в сходовой фракции с решет решетного модуля. Полнота выхода ]-го компонента исходного зернового материалы в фуражные отходы

е вф,\ ~ °П2, ЬП2, ! 0 а, /

где ьп2, - содержание Но элемента в фуражных отходах.

Полнота выхода ^го компонента зернового материала в фуражные отходы после триерного блока

е еф]Т = 0ФТ ЬфТ, ! бмо Ьм /

где 0ФТ, ЬФТ, - количество и содержание ] - го компонента во фракциях,

выделенных ТБ в фуражные отходы.

Для операции (распределение зернового материала по ширине

скельператора) очевидно, что е Ь,РУ = 1. Воздействие этой операции на последующую определится плотностью вероятности /о (в) распределения подачи О зернового материала по ширине В скельператора. Рассмотрим процесс пневмосепарации гетерогенной сыпучей среды в пневмоканале шириной В и глубиной S (рис. 2) с подачей Q в него гетерогенной сыпучей среды с известной или задаваемой плотностью вероятностей /о (В) распределения этой подачи по ширине бокового ввода в пневмоканале и плотностью вероятности /у (В) распределения средних по сечению пневмоканала скоростей воздушного потока.

Для этих условий полнота прохода ]-го компонента (]=1,2,...Ь) зернового материала в очищенную на этой частной операции фракцию [2]

е

0е р, оу )

п Р = 1 ,

е ЬП = ----------------т ,В т ш,

Ь,П в рк ,В

где К = В /е - количество зернового материала Qp, попадающего на каждый р-ый участок (р=1,2,...,е) к-й ширины рабочего органа:

сепаратора

_ В 0 р- 1 ''р~ * е

р

Л (в) лв

В0

О,

Л ( в ) лв

р= 1

а средняя скорость воздушного потока р-м участке

р

Л, (В) лв

в V ,-1

ур = *

Л, ( в ) лв

р= 1

где ер, (°рур) - полнота выхода в очищенное зерно >го компонента зернового материала с р-го участка ширины пневмоканала (р = 1,2, ..., е ) определится из известных выражений [1,3] при Q= Q Р и

V = Ур

Содержание сорных О = 1,2,...,с) примесей Ьс в очищенном зерне

Ь

Ьс _ а,е Ь,П ! а,е Ь,П ,

,= с ,= 1

полнота выделения пневмосепаратора в отходы сорных 0 = 1,2,...,с) примесей

Ь

е Ьс = а, (1 -е ЬП ) / а, ,

3 = с ,= 1

содержание ;Н х компонентов в очищенном зерне

Ь

Ь, _ а,е ь,п ! а,е Ь,П ,

,=1

полнота выделения j-го компонента из зернового материала в пневмосепараторе

е Ь, = (1 - е Ь,П ) ,

содержание Ь ,1 j-го компонента зернового материала в легких отходах

Ь

Ь,1 _ а,(1 - е Ь,п )/( а,(1 - е Ь,П ) ) ;

,=1

масса выделенных отходов

Ь

рл _ 0а, (1 - е Ь,П ) ;

,=1

е

масса очищенного в пневмосепараторе зерна

Qкл= Q - Q .

Входное воздействие на рассматриваемый решетный модуль (см.

рис.1, частные операции) выражается вектором Fp, независимые аргументы которого случайные в вероятностно-статистическом смысле величины Рр = Ьрм. а„. W, М (*Да 2 (Ь1), / (В,), / (Н), / (I,) ]. Активные средства, определяющие эффективность функционирования решетного модуля, определялись вектором Арм (см.рис.1), где ар1 . ар; - подача зернового материала Qpм=Qк из (1), содержание в нем ^х компонентов; /а (В,), /а (Н) - плотность вероятности распределения Но компонента по ширине В, иго решета в ярусе и по высоте ярусов в решетном модуле; /(|) - закономерность поступления зернового материала на решета.

Приняв известную гипотезу об аддитивности процесса сепарации сыпучего материала на р -х участках решетных ярусов, полнота просеивания ] -го компонента на 5 -м решетном ярусе (5=1,2,...,4) для воздушнорешетной машины ОЗС-50/25/10) при известных /а (Н) и fq (В) определится выражением [3]

г р

е /а(Н)аН / ,тв

а ^------------ £Г----------- а}5 £ 5 ^5 рР

/а (Н ^Н /5 q(B)dB ,

г=1 р=1

£5 Я 5 а3

з

где £ 5 ;р =£5 1 ] + (1 - £ 5 1; )£5 2; .

При этом подача сыпучего материала на 5 -й решетный ярус

ч$ =а рка .

Полнота просеивания ] -го компонента на р -м участке 5 -го решетного яруса при подаче 45р сыпучего материала на этот участок определится из различных выражений [1] в зависимости от многих факторов:

£5 ; ( 45 р ) = / ( а. а} .W.g . /. Ц5 т. Лт ,а т. Р т . Щт. Nт. П ) , где /; .7 - вид сыпучего материала (зерно, семена трав и др.) и его плотность; Ц5т - длина т -го решета 5 -го решетного яруса; d5m - рабочий размер отверстий т -го решета 5 -го яруса; а 5т, Р 5т, Щт, N5m - угол наклона к горизонту, направленность, амплитуда и частота колебаний т -го решета в 5 -м ярусе; П - ко-

личество решет в ярусе.

При известной полноте просеивания £ 5; ; -го компонента на отдельном ярусе решет (5= 1,2,...,к) полнота просеивания ; -го компонента на четырех ярусах решетного модуля определится из выражения [3]

е

X Я. о' Яіі £ о}

о=і___________

&а/к

где а8ъ - содержание Ъ -го компонента в сыпучем материале, поступаю-

С* V ^

щем на 8 -й решетный ярус.

Полнота схода Ъ -го компонента с 8т -го решета [3]

Х Ч8траЪе СХ8тЪ [я8тр ) р= і ,

£ СХОтЪ = ]

Я8та

а со всех т -х решет к -ярусного решетного модуля

к

£ СХтъ = X РЩ£ СХ8 тъ

8 = і

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

Для реализации поставленных задач исследований с использованием методов математического моделирования и выбранного критерия оптимизации Еф рассмотрены три схемы ЗОА (см. рис.1): 1 - одна ВРМ

ОЗС-50/25/10; 2 - одна ВРМ ОЗС-50/25/10 и последовательно функционирующий ТБ ЗАВ-10.90000А; 3 - одна ВРМ ОЗС-50/25/10 и два параллельно функционирующих ТБ ЗАВ- 10.90000А.

В качестве исходного зернового материала принято зерно пшеницы, поступающее для первичной очистки в ЗОА после зерноуборочного комбайна, с усредненным (ЮФО, Ростовская область) содержанием основных j-х компонентов: зерно пшеницы (0,886), зерновые примеси (0,040), мелкие сорные примеси (0,010), полова (0,005), соломистые примеси (0,003), крупные примеси (0,040), овсюг (0,001), дробленое зерно (0,015). Состав зернового материала варьировался в следующих интервалах: щуплое зерно {0,025; 0,03; 0,04; 0,05; 0,06; 0,07; 0,08}, дробленое зерно {0,005; 0,01; 0,02; 0,03}, овсюг {0,001;0,01;0,02;0,03}. Подача зернового материала в зерноочистительный агрегат варьировалась следующим образом: при функционировании агрегата по первой и второй схемам: {6; 9; 12; 15; 18; 21} т/ч; по третьей схеме: {12; 15; 18; 21} т/ч.

20

15

Q,т/ч 10 5 0

25 20

15

0,т/ч

10

5 0

0,5 1 1,5 2 2,5

____________АІ,%

ВРЗОМ 0ЗС-50/25-10

отделение очистки ЗОА

16 14 12 10

0,т/ч 8

6

4 2 0

0,1 1 2 3

_________________АІ,%(овсюг)___________

ВРЗОМ 0ЗС-50/25-10

отделение очистки ЗОА

а

2,5 3 4 5 6 7 8

ВРЗОМ 0ЗС-50/25-10

отделение очистки ЗОА

Рис.3. Изменение производительности Q агрегата при вариации содержания в исходном зерновом материале щуплого зерна, дроблёного зерна и овсюга при функционировании одной ВРМ

Параметрическая оптимизация зерноочистительных машин в агрегате (рабочая скорость воздушных потоков в аспирационных системах, форма и размеры отверстий на трех решетных ярусах, частота колебаний решетного модуля при известных [1] плотностях вероятности ^р(В) распределения подачи зернового материала Ю) по ширине (В) пневмоканала и решетных сепараторов, а также рабочей скорости (Мо) воздушного потока К(В) и др. для дискретно задаваемой величины подачи исходного зернового материала и его технологических свойств (влажность, размерные характеристики, скорость витания основных компонентов [1] содержание основ-

ных компонентов) и оценка показателей функционирования отдельных зерноочистительных машин и агрегата в целом проведены на ЭВМ по программе «SH1 S». Выявлены рациональные параметры машин в зерноочистительном агрегате и оценены основные показатели функционирования по различным схемам. Проведена оценка величины производительности различных рассматриваемых агрегатов при вариации подач зернового материала и его технологических свойств (различное содержание в исходном зерновом материале его основных отделяемых компонентов: щуплое зерно, дробленое зерно, овсюг) при обеспечении агропоказателей процесса сепарации. Результаты моделирования представлены на рис.3 и 4.

Qj/ч

25

20

15

0,5 1 1,5 2 2,5

Аj,%

отделение очистки ЗОА

ВРЗОМ 0ЗС-50/25-10

Qj/ч

0,1

1

2

3

________Aj,%o(QB^r)_________

отделение очистки ЗОА

ВРЗОМ 0ЗС-50/25-10

Рис.4. Изменение производительности Q агрегата при вариации содержания в исходном зерновом материале щуплого зерна, дроблёного зерна и овсюга при функционировании одной ВРМ в составе отделения очистки ЗОА

Результаты моделирования. Разработанная математическая модель ЗОА и проведенное моделирование процессов функционирования современных вариантов зерноочистительных агрегатов при первичной очистке зерна пшеницы показали, что при росте содержания в исходном зерновом материале зерновых примесей (щуплое зерно более 4%, дробленое зерно более 1,5%) производительность ВРМ ОЗС-50/25/10 существенно снижается (до 14,55 т/ч). При содержании в исходном зерновом материале овсюга

0,1% производительность ВРМ 0ЗС-50/25Д0 достигает 14,55 т/ч, а при содержании в исходном зерновом материале овсюга более 0,1% агротребования на процесс очистки зерна в воздушно-решетной машине не выполняются. Использование одного триерного блока ЗАВ 10.90000А совместно с ВРМ ОЗС-50/25/10 (схема №2) несущественно изменяет производительность отделения очистки ЗОА (6 т/ч - 11,25т/ч).

Увеличенное количество триерных блоков (до двух) (схема №3) существенно улучшает технологические показатели ЗОА: при варьировании содержания зерновых примесей (щуплое зерно, дробленое зерно) в пределах 2,5-8% и 0,5-2% соответственно и овсюга в пределах 0,1-3% производительность отделения очистки ЗОА повышается до 21 т/ч, но при содержании в исходном зерновом материале дробленого зерна более 2,5% агротребования не выполняются.

Выводы. Обоснованы и разработаны современные методы математического моделирования процесса функционирования зерноочистительного агрегата, которые позволили, с использованием методов параметрического синтеза, количественно оценивать функциональные показатели отделения очистки зерноочистительного агрегата с одним или двумя триерными блоками. Установлено, что рост эффективности триерования зернового материала значимо влияет на показатели функционирования отделения очистки ЗОА, на его адаптацию к изменению технологических свойств зернового материала и его подачу в агрегат. Целесообразность использования одного или двух триерных блоков к одной воздушно-решетной машине типа ОЗС-50/25/10 в отделении очистки зерноочистительного агрегата определяется эффективностью триерных блоков с учетом экономического анализа функционирования агрегата с различными структурами.

Библиографический список

1. Ермольев Ю.И. Технологические основы интенсификации процесса сепарации зерна воздушно-решетными зерноочистительными машинами и агрегатами: Дис.... д-ра техн. наук. / Ю.И. Ермольев. - Ростов н/Д, 1990. - С.300.

2. Ермольев Ю.И. Тенденции и перспективы развития технологий

и технических средств для семенной очистки зерна/ Ю.И.Ермо-льев, М.В.Шелков, М.Н.Московский // Известия высших учебных заведений. Северо-Кавказский регион. Серия «Технические науки»- Ростов н/Д, 2005. - С.112-119.

3. Ермольев Ю.И. Моделирование процесса сепарации зерна в

воздушно-решётной зерноочистительной машине / Ю.И. Ермольев, М.В. Шелков // Научные основы реше-

ния проблем сельскохозяйственного машиностроения. - Тула: Тул. гос. ун-т, 2003. - С. 86-95.

Материал поступил в редакцию 23.10.07.

U.I. ERMOLJEV, M.U. KOCHKIN

MODELING OF THE PROCESS OF FUNCTION A CLEANING UNIT OF GRAIN

A general approach of a mathematical model building of function a cleaning unit of grain with different structural schemes is given, a rational parameters and functional indexes are synthesis.

ЕРМОЛЬЕВ Юрий Иванович (р.1943), доктор технических наук, профессор, заведующий кафедрой «Сельскохозяйственные машины и оборудование» ДГТУ. Окончил РиСхМ (1969).

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

Имеет более 200 научных публикаций (в т.ч. три монографии).

КОЧКИН Максим Юрьевич (р.1981), аспирант кафедры «Сельскохозяйственные машины и оборудование» ДГТУ. Окончил ДГТУ в 2004 году. Специалист в области сепарации сыпучих сред.