Моделирование воздушно-решетной зерноочистительной машины

Бесплатный доступ

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

Процесс, показатели сепарации, математическая модель, результаты моделирования

Короткий адрес: https://sciup.org/14249477

IDR: 14249477

Текст научной статьи Моделирование воздушно-решетной зерноочистительной машины

Еф§м = {F,A,G§[k§m(x),Т§m(х)]}^max ;(1)

AM c Am , F c F, x e Gm (x, u);(2)

be <[bc], bпз <[bпз], §p <[§p].(3)

Здесь F – вектор входных воздействий на принятую в зерноочистительной машине систему операций.

F = {q, aj, W, M (bj), ^(bj ), fQ (B), fQ (H) } ,(4)

А – вектор управляющих факторов системы, обеспечивающих ее функционирование,

^A = {B ш, hn,, Vei, fq (Bn), fv (Вп), Pcx (P), fQx (H ), /qx (B), К §m (х)} .

Векторы входных и управляющих воздействий на пневмосепаратор 1 F и A П1 , на решетный модуль – F Р и A Р , на пневмосепаратор 2 – F П2 и A П2 показаны на рис.1.

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

В формуле (1) х - набор технологических операций из GM ( x , u ) множества, реализуемых в К 5 м функциональной схеме ВРЗОМ. G 5 [ К 5 м ( х ), T 5( х ) ] - функционал, определяющий показатели технологического процесса в ВРЗОМ (полнота выделения j-х компонентов, потери зерна, содержание в очищенном зерне других j -х компонентов), для принятой К 5 м ( х ) схемы.

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

В м = { е ф , е bj ц , w , b j , b ю, b с, 5 3 i , 5 з , Q o , Q ф , Q ок , е jxi , y n jxi , y n xi , y jxi , y xi , b n jxi , b jxi } ,

где Еф - критерий эффективности реализации технологического процесса ВРЗОМ; еbjц - полнота выделения из зернового материала отделяемых j-х компонентов; 5зi, 5з - потери зерна на xi операции и общие; Q0, Qф, Qок - масса выхода очищенного зерна, фуражных отходов и отходов разных категорий; ynjxi, ynxi - масса j-х и всех компонентов, содержащихся в зер новом материале после реализации xi -й операции; yjxi,yxi - масса j-х и всех компонентов, выделенных из зернового материала после реализации xi -й операции; bnjxi - содержание j-го компонента в массе очищаемого материала и в выделяемых фракциях bjxi; е jxi -полнота прохода j-го компонента в очищаемый зерновой материал после выполнения xi -й операции.

Ъ-е,

Е ф = 1 - ^ ----,                                (7)

L a j = 2

функция цели:

E ф ^ max .                                  (8)

Обозначения элементов, входящих в выражения (1)-(8), определены в предыдущей работе автора [1].

В общем виде плотность вероятности распределения зернового материала по ширине В пневмоканала и решетного модуля - /q (B), а плотность вероятности распределения скорости воздушного потока в пневмоканале - /V (Bп). Тогда количество зернового материала Qp, попадающего на каждый р-й участок (р=1, 2, …, е) k-й ширины рабочего органа, f/Q (в )dB

Q p = BQ0- p - 1 -------,                                (9)

f /q (B )dB p а средняя скорость воздушного потока на р-м участке:

f / v ( B ) dB

Vp = BV p - 1 .                                 (10)

p ke

f /v (B) dB

p

Полнота выхода j -х компонентов зернового материала в очищенное пневмосепаратором зерно:

t Q, ^ ( Q p V )

£   = p '-------------------, nkj          Qe

где enkj(Qp,Vp) - полнота выхода J-го компонента в очищенное пневмосепаратором зерно для ус ловий сепарации, определенных подачей Qp зернового материала и скоростью V, воздушного потока на р-м участке пневмосепаратора e S<2 ,V )= fГQ ,W,a.,MM ,f (B), f (B),B ,h ,f 1,                    (12)

nkj \ p p /           p        j 0 vj v        Q         ni ni ni где MvJ - математические ожидания средних скоростей витания j-х компонентов в исходном материале; Вn, - ширина i-го пневмоканала; hni - глубина i-го пневмоканала; fni- условия ввода зернового материала в пневмоканал.

Конкретные математические модели, определяющие еnkJ для вертикальных пневмокана- лов, пневмосепараторов с активным вводом зернового материала в вертикальный или наклонный пневмоканал определены [1].

При выполнении ограничений (1)-(3) можно рассчитать показатели процесса сепарации:

- содержание сорных примесей (j' =1, 2, ..., с) в очищенном зерне bC = L «^ [^.епJ ; j=C         V j=i          )

( 13)

- полнота выделения ebc пневмосепаратором в отходы сорных (j' =1, 2, „., с) примесей е b c

= t [ aj (1 - e п Kj ) ] / ^Laj ' J = C                     j = 1

– содержание j -х компонентов в очищенном зерне

b

= a j e C j

b

/ L

V j = 1

I a j e п Kj ;

)

– полнота выделения j-го компонента зернового материала в пневмосепараторе ebj =(1-ек);                                          (16)

– содержание j-го компонента зернового материала в легких отходах bo = a J '1 e^. )

;

– масса отходов и масса очищенного в пневмосепараторе зерна b

Q O = L Qa j (1 п kj ) , кг/(м - с) ; Q k = Q - Q о , кг/(м - с).                 (18)

j = 1

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

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

F p = [ Q p , a pj ,W , M ( b ), a 2 ( b ), fQ ( B i ), fQ ( H ), f ( / i ) ] .                       (19)

Активные средства, определяющие эффективность функционирования решетного модуля, найдем вектором А р :

Л = [ a p , P i , R i , n, T, l i , B i , b i , c , K , ( x ) ] , (20) где все показатели определены [1].

Различные базовые решетные модули, которые можно использовать в ВРЗОМ, математические модели, описывающие процессы функционирования базовых и комплексных (составленных из различных структур базовых модулей) решетных модулей, определены [1].

Рассмотрим, для примера, в общем случае многоярусный решетный модуль (рис.2) с неоднородными по длине ( 1 ц,... 1 1 n ) и высоте (1, 2, „., к ) решетами и с различными кинематическими параметрами. В общем виде на ярусы решет поступает q 1 j ,q 2j ,..., q nj ,q K j количество j -го компонента сыпучего материала с их торцов, определяемое плотностью вероятности распределения f Q ( H ) случайной величины Q (подача сыпучего материала на ярусы (1, 2, „., к решет)) по высоте H расстановки ярусов решет: Q=QK (18), содержание j -х компонентов в сыпучем материале b j определяется из выражения (15). Распределение случайной величины q 5 по ширине В р сепараторов в ярусах определяется плотностью вероятности f ( B ) В качестве допущения здесь и далее будем считать поступление сыпучих материалов на ярусы решет постоянными по времени (процесс квазистатичный) и с равномерным распределением j -х компонентов в сыпучем материале ( j =1, 2, …).

а )                                                      б )

Рис.2. Многоярусный решетный модуль ( а ). Плотность f g ( H) вероятности поступления зернового материала на к - е решетные ярусы и плотность вероятности f 5 q ( B ) его распределения по их ширине Bp ( б)

Полноту просеивания j -го компонента на 5 -м решетном ярусе можно определить по из- вестному выражению [1]:

p = 1 £ 5 j =

Е q 5 p a j^ ykq 5 p )

q 5 a j l

где    q 5 p - приведенная подача сыпучего материала на p -й участок (по ширине B )

BQ

(p = 1,2,..., e) n -го решетного яруса, q5p = ""^"P5Qp ; k=B\e — ширина p -го участка яруса решет; P5qp - вероятность попадания случайной величины Q на p -й участок по ширине

5 -го решетного яруса.

P 8 Qp = P rQ P 5 pq .                                  (22)

Здесь Р Г q - вероятность попадания случайной величины Q на Г - ( Г - 1 ) участок ( 5

решетный ярус):

Г

J f Q ( H ) dH

Г - 1

P r Q    С            ;

J f Q ( H ) dH

Г = 1

J f q ( B ) dB

P 5 pq =              -                           (23)

J f 5 q ( B ) dB

p = 1

При этом подача q5 сыпучего материала на 5 -й решетный ярус q5= QPrQk;                                     (24)

85 j ( q 5 p ) - полнота просеивания j -го компонента на p -м участке 5 -го решетного яруса при подаче q 5 p сыпучего материала на этот участок, определяется из различных выражений в зависимости от многих факторов:

85j (q5p ) = f (Q, aj , W, Y, fj , L5m , dm , am , Pm , Rm , Nm , n), где fj, у - вид сыпучего материала (зерно, семена трав и др.) и его плотность; L5m - длина m - го решета 5 -го решетного яруса; d5m - рабочий размер отверстий m -го решета 5 -го яру са; а5m,р5m,R5m,N5m - угол наклона к горизонту, направленность, амплитуда и частота колебаний m -го решета в 5 -м ярусе; n - количество решет в ярусе.

l     jfQ ( H ) dH J f 5 q ( B ) dB

A Q r----------pr--------5jp (q5p)

p = 1 J f Q ( H ) dH j f 5 q ( B ) dB

r = 1                 p = 1

8 5 j —

A q 5 a j

j где 85jp    8(n-1)ojp + [1   85(n-1)ojp ] 85njp .

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

k

A q 5 a 5 j 85 j

= 5=1 __________

,

8j"   Qajk где a5j - содержание j -го компонента в сыпучем материале, поступающем на 5 -й решетный ярус.

Сход j -го компонента с p -го участка 8 m -го решета определяется из выражения

Y, -

8 mpj

m

= q 8 р а7 П a 8 1 р ( h ) exP

m

2y 2B, 202.b, V-V-L, I 8 i 8 pi p8 ij 8 ij 8 ij 8 i

i = 1

V

4 q 8 ip

У.

.

Принимая, что для каждого p -го участка и для каждого 8 m -го решета параметры сепарации а 8 mj ( h ) , b 8 mj , V 8 mj различны, но постоянны в пределах ограниченной длины 8 m -го решета, массу просеивания j -го компонента сыпучего материала на длине Х p -го участка 8 m -го решета можно определить из выражения

У

Y n 8 mj ( X ) = Y 8 ( m 1)/

1 a 8 mj ( h ) e xP

Y 8 m 8 m p8 mj 8 mj

V

{Z Y .( m —O j

У У

а полноту просеивания j -го компонента на длине l 8 m -го решета модуля - из выражения

l

^^, q 8 mp a 8 mj 88 mj ( q 8 mp , L )

_        P = 1

8 8 mJ

J f Q ( H ) dH mi

E Qa 8j --------П(1 — 88 j) aj l j       J fQ (H) dH i=0

8= 1

.

Полнота просеивания j -го компонента 8mj через все m -е решета к ярусов определяется из выражения

k

8 mj = ^ 8 8 mj

8= 1

r

J fQ (H) dH r—1

c

J f Q ( H ) dH

8= 1

.

Полнота просеивания j -го компонента 8j через все решета k -ярусного решетного моду ля

8 j = 8 ( n 1 ) oj + [ 1 — 8 ( n 1 ) 0 j ] 8 nj ,

где 8 ( n 1 ) 0 j - полнота просевания j -го компонента через все ( n — 1 ) последовательные решета в k ярусах.

Полнота схода 8 CX 8 mj j -го компонента со всего 8 m -го решета

8 СХ 8 mj =

q 8 mp a j 8 cXmj ( q 8 mp ) p = 1

q 8 m a j l

,

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

k

8 CXmj = ^ P ГQ 8 CX 8 mj .

8= 1

Полученные выражения (19)-(33) позволяют оценить проход j -х qпδmj и всех qпδm компо- нентов сыпучего материала через δm -е решета отдельных ярусов и всех k ярусов Qпδmj , Qпδm и содержание j-х компонентов в проходе каждого m -го решета δ -го яруса bδmj и всех m -х ярусов bmj [1].

Входное воздействие на второй пневмоканал определяется вектором F П2 , активные средства – вектором А П2 (рис.1).

Математическая модель оценки показателей функционирования второго пневмоканала аналогична (9)-(12) при подаче зернового материала в пневмосепаратор 2, равной выходу очищенного в решетном модуле зерна. Проведенные исследования показали, что построенные математические модели с 95-процентной (полнота выделения сорных примесей – 90%) доверительной вероятностью адекватно описывают процесс сепарации в воздушно-решетной зерноочистительной машине [1].

Параметрический синтез. Параметрический синтез функциональной схемы и задаваемой техническим заданием вариации величин аргументов векторов входных F (4) и управляющих A (5) воздействий проводится с использованием соответствующих математических моделей методом нелинейного программирования – методом сканирования с ограничениями, так как отсутствует унимодальность функции цели (8) и нелинейные ограничения (2)-(3).

На втором иерархическом уровне синтезируется функциональная подсистема решетного модуля (РМ) с делителем зернового потока (ДЗП) для ВРЗОМ.

Для синтезированных рациональных параметров рассматриваемой структуры РМ (например, схема РМ, число решетных полотен, форма и размеры отверстий решет, кинематические параметры: частота, амплитуда колебаний, угол наклона решет к горизонту, угол направленности их колебаний) проводится многомерный анализ процесса функционирования РМ при различных величинах вектора выходных F p (19) и управляющих A p (20) воздействий, определяемых выходными показателями функционирования ПС.

Эффективность сепарации зерновых (пшеница, рис) и семян подсолнечника с учетом ограничений размеров решетных полотен (ТУ 23.2.2068-89), изготавливаемых в России, а также приводов решетных модулей и систем их очистки от забивания отверстий выявлена для одноярусных трехрешетных модулей.

Определен рациональный решетный модуль с параллельной и последовательной схемами сепарации зернового материала, функциональные показатели которых соответствуют заданным, для очистки зерна с учетом проведенного параметрического и структурного синтеза комплексных решетных модулей [1]

С использованием математической модели (1)-(33) был разработан алгоритм программ для ЭВМ ПК, позволяющий проводить параметрический синтез (8) с учетом принятых ограничений (2)-(3) и значений аргументов векторов F , А , A p, Bм [1].

Для проведения многомерного анализа оценки основных показателей функционирования обоснованного варианта ВРЗОМ (рис.3) при вариации аргументов вектора F (рис.1) осуществлено моделирование процесса системной сепарации зерна.

Рис.3. Технологическая схема воздушно-решетной зерноочистительной машины: 1 – камера приемная; 2 – скельператор барабанный; 3, 4 – каналы аспирационные;

5 – камера осадочная; 6 – 4-ярусный решетный модуль

При моделировании приняты следующие допущения: плотность распределения вероятности поступления зернового материала по ширине решетного модуля подчиняется нормальному закону с коэффициентом вариации 20,26%; плотность вероятности распределения зернового материала по высоте решетных ярусов характеризуется коэффициентом вариации X q ( H ) , изменяющимся в пределах от 0 до 54,4%.

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

Полнота выделения мелких сорных примесей s„ (рис.4, 6 ) в интервале подач от 2,775 до в см

4,635 кг/с незначительно уменьшается при изменении коэффициента вариации от 0 до 54,41%. Для интервала подач 4,635-8,58 кг/с sB незначительно снижается при изменении коэффициента см вариации от 0 до 23%. При дальнейшем увеличении Xq(H) происходит существенное снижение полноты выделения мелких сорных зерновых примесей. Агропоказатель bc< 1% (bc - содержание сорных примесей в очищенном зерне) выполняется на всем рассматриваемом диапазоне вариации подач Q и Xq(H) (рис.5,а).

Существенное влияние изменения коэффициента Xq(H) и подачи оказывают на выделение из рассматриваемого материала зерновых примесей sвпз на всем интервале изменения подач. Например, для Q = 3,705 кг/с при изменении Xq(H) от 0 до 54,4% sвпз снижается на 5%, а для Q = 8,115 кг/с при изменении Xq(H) в том же интервале sвпз снижается на 7%. При этом агропоказатель bПЗ< 2% (содержание зерновых примесей в очищенном зерне) выполняется при различных значениях Xq(H) для разных подач. Так, при Xq(H) =24% максимально допустимая по- дача составляет 6,495 кг/с (23,2 т/ч). Для Q = 5,65 кг/с и более увеличение Xq(H) существенно влияет на потери зерна в фуражные и неиспользуемые отходы. Так, для Q = 6,495 кг/с макси мально допустимое значение Xq(H) составляет 35,66%, а для Q = 8,115 кг/с Xq(H)max = 18,1%.

Рис.4. Зависимость полноты выделения компонентов из зернового материала, очищенного воздушно-решетной машиной от подачи Q зернового материала в машину и величины коэффициента вариации X неравномерности поступления зер нового материала на отдельные решетные ярусы: а - крупных £вкп и б - мелких сорных £„ всм

Выявлена усредненная допустимая неравномерность поступления зернового материала на решетные ярусы в пределах Xq(H) =28-30%. При синтезе основных параметров зерноочистительной машины для различных основных величин Q и /q (H) усредненные рациональные па раметры: Vb = 6 м/с, n = 7,92 об/с, Ь1 = 1,7 мм, b2 = 2,0 мм, b3 = 3,6 мм.

Рис.5. Зависимость содержания в зерне, очищенном воздушно-решетной машиной, сорных примесей bс (а) и чистоты aпо (б) зерна от подачи Q зернового материала в машину и величины коэффициента вариации X неравномерности поступления зернового материала на отдельные решетные ярусы

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

Список литературы Моделирование воздушно-решетной зерноочистительной машины

  • Ермольев Ю.И. Проектирование технологических процессов и воздушно-решетных и решетных зерноочистительных машин/Ю.И. Ермольев, А.В. Бутовченко, М.Н. Московский, М.В. Шелков. -Ростов н/Д: Издательский центр ДГТУ, 2010. -638 с.
  • Ермольев Ю.И. Технологические основы интенсификации процесса сепарации зерна воздушно-решетными зерноочистительными машинами и агрегатами: дис. … д-ра техн. наук: 05.20.01. -Ростов н/Д, 1990. -345 с.
  • Ermol'ev Yu.I. Proektirovanie tehnologicheskih processov i vozdushno-reshetnyh i reshetnyh zernoochistitel'nyh mashin/Yu.I. Ermol'ev, A.V. Butovchenko, M.N. Moskovskii, M.V. Shelkov. -Rostov n/D: Izdatel'skii centr DGTU, 2010. -638 s. -In Russian.
  • Ermol'ev Yu.I. Tehnologicheskie osnovy intensifikacii processa separacii zerna vozdushno-reshetnymi zernoochistitel'nymi mashinami i agregatami: dis.… d-ra tehn. nauk: 05.20.01. -Rostov n/D, 1990. -345 s. -In Russian.
Статья научная