Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Липкович, Э. И. Процессы обмолота и сепарации в молотильных аппаратах зерноуборочных комбайнов (пособие для конструкторов зерноуборочных машин)

.pdf
Скачиваний:
23
Добавлен:
19.10.2023
Размер:
6.6 Mб
Скачать

imя интенсивности (скорости) процесса сепарации для При­ веденных значений количества обмолоченных зерен па входе.

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

Вероятностная модель процессов обмолота и сепарации: стационарный случай

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

Если процесс обмолота1зерен в подбарабанье представить как вероятностный процесс, а относительное количество об­ молоченных зерен интерпретировать как вероятность обмоло­ та, то, очевидно, вероятность Y(t) пребывания свободных зе­ рен в подбарабанье представится произведением:

Y (t)= x (t)y (t),

(146)

где y(t) — вероятность задержки зерен в потоке вороха в под­ барабанье.

Для раскрытия (146) необходимо определить вид функции y(t), что можно сделать на основе исследования процесса сепарации.

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

70

Для всех непересекаюіцихся временных отрезков интервала (о, і), независимо от того, где эти интервалы расположены, II независимо от прошлого состояния системы. Это значит, что процесс имеет два свойства: стационарность и отсутствие последействия.

Течение вероятностного процесса описывается некоторой функцией времени f(t). Пусть

y (t)= f(t).

(147)

Тогда вероятность Z (1) выделения зерен из подбарабанья за время t определится разностью:

Z ( t) = l - f ( t) .

(148)

Определяя вероятность cp(t) выделения зерна за малый промежуток времени на основании условий стационарности

иотсутствия последействия, сделаем следующие допущения:

1.cp(t) —A,t-)-0(t), А,>о (/\.=const).

2.Вероятность того, что за малый промежуток времени произойдет больше, чем одно выделение зерна, есть o(t). Вто­ рое допущение порождает требование ординарности — прак­ тической невозможности появления двух или нескольких со­ бытий за малый промежуток времени.

Тогда

 

 

f (t) = 1 —«p (t) = 1 — [Xt+o(t)].

(149)

Здесь o (t) — величина,

убывающая быстрее, чем t.

На основании теоремы

о произведении

вероятностей по­

лучим:

f(t+Ät) =f(t)f(A t) = f(t) [1—Ф(At)] = = f(t) { 1 -[Ш + о (Д ф ] )•

Следовательно,

Ä = - M ( t ) , dt

откуда f(t)=Ce~~At.

Так как f(о) = 1, то С= 1 и

f ( t) = e ~ Xt.

(150)

(151)

(152)

(153)

71

Подставляя выражение (І53) в зависимости (147) и (І48), получаем:

Y (t)= e ~ Xt,

(154)

Z (t)= l — е—Xl

(155)

Используя формулу (146) и ранее установленную зависи­ мость x(t), находим вероятность содержания свободных зе­ рен в подбарабанье:

Y(t) = (1—ае— )е—Xt.

(156)

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

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

Полученные уравнения (156) и (157) отражают в общем виде изменение частных совокупностей свободных зерен в подбарабанье и выделившихся нз него. Эти закономерности объединяют в себе оба. процесса: процесс обмолота, т. е. на­

копления свободных

зерен в подбарабанье, и одновременно

с ним протекающий

процесс сепарации.

Условие ординарности, которое заложено в выводе этих уравнений, неприемлемо для подбарабанья в целом, так как последнее может выделить в фиксированный момент времени несколько зерен нз нескольких участков по ширине и из не­ скольких слоев по толщине потока. Поэтому следует выбрать участок подбарабанья, для которого условие ординарности можно считать более пли менее выполнимым как по ширине подбарабанья, так и по толщине соломистого потока. Таким участком по ширине подбарабанья является одна ячейка пространственной решетки шириной ас (назовем этот участок минимальным). Число h минимальных участков определится выражением:

( 158)

72

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

Хз=1000 3 ^ ,

(159)

С7)і

 

где q3 — секундная подача зерна в подбарабанье, кг/сек;

Е; — вес одного зерна, г.

 

Если Хзі — количество зерен над і-ым слоем минимально­

го участка, то интенсивность Я.і выделения зерна из і-ого слоя можно определить соотношением:

Хі= ріХзі,

(160)

где рі — вероятность выделения отдельного зерна.

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

из потока зерно проникает, далее, под решетку

деки.

Если

Ро — вероятность прохождения зерна сквозь один

слой,

рд —

вероятность прохождения зерна сквозь решетку деки, і — чис­ ло слоев между зерном и решеткой деки, то полная вероят­ ность р выделения зерна под деку определится выражением:

Р= Ро‘Рд-

(161)

Каждый слой потока представляет решетку из стеблей. Рассматривая зерно как тело выпуклого контура с приведен­ ным диаметром с13> можно вычислить вероятность р0 по фор­ муле:

d:i + d

acsi п Ѳ

(162)

 

где d —диаметр стебля;

 

Ѳ — угол падения зерна;

 

ас — расстояние между осями соседних стеблей.

 

Количество зерна в і-ом слое потока вороха зависит от расположения колосьев по сечению валка. Опыты по изуче­ нию структуры валка, проведенные во ВНИПТИМЭСХ, пока­ зали, что закономерность этого распределения с хорошим приближением может быть представлена параболической зависимостью:

73

''Г^к.ѵЧ-кгѵ+кз. (163)

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

Гг:=^ (Аѵѵ2)

(°<ѵ<Д).

(164)

Для конечного числа с, слоев

имеем (164)

в дискретном

виде:

 

 

Ѵ (с і-іД

(165)

і=1

 

 

Уравнение (165) записано для симметричного случая рас­ пределения зерен относительно середины толщины валка. Все другие случаи распределения могут быть учтены коэффици­ ентом смещения вершины (или оси) параболы или отражены другими уравнениями.

Принимая распределение зерен по сечению потока в подбарабанье таким же, как и в валке, определим количество зерен в і-ом слое потока вороха:

ѵ1000q3ac(H—i2)

Х

м

i

 

 

<;r(|V (tj_.j2)

(166)

На основе выражения (160)

с

учетом формул

(162) и

(166) получаем:

 

 

 

 

^ __ЮООдзЗс р Л _

с1з4-с1 У

 

 

 

Л

а с /

£

 

 

 

 

S

(? і-і2)

(167)

Для отдельного слоя пространственной решетки уравне­ ния (156) и (157) примут вид:

Y, (t)

= (1 -а е - И )е - Ѵ ,

(168)

Z, (t) =

(1 — а е —P*)(l—e —V).

(169)

74

Ма основании представления 6 математическом ожиДанніі количество Х'3і обмолоченных зерен, находящихся в каждом слое потока в подбарабанье, составит:

Х'ві= Х зг Y,(t) =Хзі(1—ae“ ^l)e~>|1.

(170)

Следовательно,

Х'зі=Хзі (1—ае —i5l)e'“ ?'|l;

X/32= X 32( l - a e - ß 1) e - ^ t;

Х ^ = Х3с (1 -,е -Р 0 е -Ч с1

(171)

S X ,3.= X,3 = ( l - « e - P t) S X3ie-Xlt

 

i=l

 

i=l

 

и с учетом

(166)

 

 

 

x , _1000qA

> -« - ? ■

44 (. , _ iaje_ M .

 

^ l

V(Èi-j2)

 

(172)

 

 

 

 

 

i = l

 

 

Вероятность нахождения свободного зерна во всех слоях

потока вороха в подбарабанье составит:

 

Y(t) =

ХС

l - a e - ß 1

S

-X,t

=

 

S ( « - i 3)e

 

 

Хз

е

(173)

 

 

 

 

V (? i- i2) 1=

i=l

или в развернутом виде с использованием зависимости (167):

.

1 —яе—ßf J

/

Y (t)=

^ -----------^(?і—і2)ехр[ —Рд^ 1—

 

i=i

(174)

ds+dV 1000q3ac(:i—i2)

^ S C : i - i 2)

i = l

75

Аналогично можно записать и выражение для вероятно­ сти выделения свободного зерна из подбарабанья, т. е. сепа­ рирующее действие молотильного аппарата:

 

„ ...

1—ае

 

"

Л,

 

 

 

 

 

 

Z(t) =

------------- N] (Si— i2)(l — e Alf).

 

 

 

 

V (ü i-i2)

i=1

 

 

 

(175)

 

 

i=l

 

 

 

 

 

 

 

 

Преобразовав зависимость

(175), получаем:

 

 

 

 

'

 

5

U - i2

5

Si-12

_x,tl

.

Z(t) =

(1—а е -? 1)

V ----------------- V --------------- e

 

 

 

 

 

t

( 5i - i s)

 

 

 

 

-i=1 у ;

(Si— i2)

i==1y ;

 

 

 

 

 

 

i=i

 

 

i=i

 

(176)

Здесь

 

 

 

 

 

 

 

 

 

 

 

;i — i2

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=l V (H - i3j

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

Z(t) =

(l—а е “ ?1)

,

£

 

? i - i3

,

,

 

 

1— у

------------- e~x,t

 

 

 

 

 

 

 

 

s

(Si—i2)

 

(177)

 

 

 

 

 

i= 1

 

 

 

 

Нетрудно проверить, что такой же результат достигается при простой операции вычитания в соответствии с (157). Отыщем скорость процесса сепарации. Дифференцируя урав­ нение (175), находим:

ш

dt

1 ----------- E ( 5 i- i2)[(aße ?1-Х,е

- X , t

 

 

V (si —і2) i=1

 

 

 

 

1= 1-

 

 

 

При t —О

— а(ф + Х,)е_ ^ + Хі>1.

(178)

 

 

 

 

'd_Z(t)‘

_

1

5

 

 

dt t

0 “ 1

 

S ( 6 i - i 2)(l- a)X ,

(179)

 

 

V

( s i- i2) i=1

 

 

i=l

 

 

 

76

При а = 1 уравнение

(179)

обращается в нуль.

 

Исследуем структуру уравнения (179). Введем знамена­

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

 

 

 

dZit)

 

51 — 1“

(1

а)Л|

(180)

dt

t= 0

 

 

1 = 1

(?i--i2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Установим

относительное

содержание

зерна,

поступаю­

щего на минимальный участок:

 

 

 

 

Х з і

5 1 1

 

 

(181)

 

хГ

i=l

і - і 2)

 

 

 

 

 

 

 

 

Сравнивая (181) с правой частью уравнения (180), ви­ дим, что величина

$ і- і2

(1-а)

г

У (^і- і2)

і=і

представляет собой относительное содержание свободных (обмолоченных) зерен в і-ом минимальном участке подбарабанья в начальный момент. Следовательно, выражение (180) показывает, что скорость сепарации зерна в начальный мо­ мент пропорциональна относительной величине совокупности свободных зерен.

Для вычисления коэффициента пропорциональности h в соответствии с разработанными здесь аналитическими конст­ рукциями необходимо определить расстояние ас между ося­ ми стеблей. Е го величина связана с количеством стеблей в потоке вороха. Известно, что в подбарабанье стебли подвер­ гаются разрыву (перебиванию). Степень перебивания состав­ ляет 1,5—2,5 и даже более, в зависимости от секундной по­ дачи, технологического режима молотильного устройства, фи­ зико-механических свойств растительной массы. Однако, пренебрегая выделением соломы под деку, будем рассматри­ вать многослойную пространственную решетку состоящей из соломы первоначальной длины. Далее предположим, что вы­ тягивание потока растительной массы в подбарабанье про­ исходит путем уменьшения его толщины, а параметры от­

/ /

дельного слоя соломистой решетки остаются неизменными. Указанные допущения соответствуют общей структуре выве­ денных уравнений, которая предусматривает постоянство на­ чальных условий во всем интервале (o,t).

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

 

 

 

 

TiLd

 

 

(182)

 

 

 

 

 

 

 

 

где

р — средний вес одного стебля;

 

 

 

 

L — средняя длина стебля;

 

 

 

 

 

уі — объемный вес соломы.

На рис. 20 представле­

x ft),y (t),z(t)

 

 

 

 

 

]

ны графики функций x(t)>

 

 

7ІГ)

0,8

 

y(t)

и Z(t) для молотиль­

 

г п

’/ft)

 

ного аппарата типа СК-4

0,5

 

при допустимой величине

0,4

 

 

j(tl

 

подачи

растительной мас­

 

 

 

сы

в

молотилку

qm=

0,2

.

 

 

 

 

 

 

= 6,5 кг/сек. Из графиков

0

 

 

-

 

видно, что характер функ­

0

 

 

 

ЬЖгсек' ций

сепарации, построен­

 

 

 

 

 

ных на основе вероятнос-

Рис.

20. Графики

функций

x (l), Y (i),

Тной стационарной

моде-

ZSen™ вероятностион стац"онарной ли, такой же, как у де­ терминированной модели.

Поэтому эти кривые анализировать не будем.

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

78

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

Вероятностная модель процессов обмолота и сепарации: нестационарный случай

Построение модели с переменной интенсивностью выпол­ ним с помощью теории ветвящихся процессов — процесса «рождения» и «гибели» [46, 52, 58, 59].

Для производящей функции F (t,s;x) при непрерывном времени записываются следующие уравнения в частных про­ изводных:

dF(t,s;x)

-f(s,F);

(183)

ds

 

 

 

 

 

 

dF(t,s;x)

_ V

dF

 

ч——

(184)

<?t

k=l

UXK

 

где ( x j < 1, k = l, 2, .... n.

Пусть А(т) — плотность вероятности сепарации зерен («гибель»), в общем случае это функция времени; ц(т) — плотность вероятности обмолота («рождение») — также пе­ ременная во времени величина.

Построим функцию

Г(х,х) = Х(т)-[Х(т:) + іх(т)]х-Ьіх(т)х2 = (1-х)[л(-)-іі('г)].

(185)

Тогда дифференциальное уравнение (183) для F (T ,S ;X )

— (T,S—* = - X (s ) + [X (s)+ n (s)]F -p (s)F 2

(186)

есть уравнение Риккати. Отсюда получаем F ( S ,T ;X )

в дробно-

линейном виде:

 

 

F(T,S;X) = s‘+ ( i —V Л1)*

(187)

 

1—ті‘х

 

 

Для определения функции

и ц' применим второе диф­

ференциальное уравнение (184), которое перепишется с ис­ пользованием соотношения (185) в следующем виде:

79

Соседние файлы в папке книги из ГПНТБ