Файл: Липкович, Э. И. Процессы обмолота и сепарации в молотильных аппаратах зерноуборочных комбайнов (пособие для конструкторов зерноуборочных машин).pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 24.10.2024

Просмотров: 100

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

дельного слоя соломистой решетки остаются неизменными. Указанные допущения соответствуют общей структуре выве­ денных уравнений, которая предусматривает постоянство на­ чальных условий во всем интервале (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



dF = (1 -х )[Л (т)-д (т)х ]

ox

.

( 188)

"д-

 

 

 

He проводя дальнейших выкладок, которые могут пред­ ставлять лишь технический интерес, окончательно получаем вероятность yv (t) нахождения зерен в подбарабанье (сво­

бодных вместе с необмолоченными):

Уѵ(*) = е Р-

(189)

t

 

Р= f Ж - i ^ l d T .

(190)

о

 

Для того, чтобы «материализовать» полученное уравнение, следует определить вид функций /,(т) и ц(т).

Плотности вероятностей обмолота («рождения») и сепарации («гибели»)

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

x ( t ) = l —аеН 31,

(191)

где x ( t) — вероятность обмолота;

 

а — вероятность

нахождения необмолоченных зерен

в ворохе после первого удара бича.

 

Тогда плотность вероятности обмолота p,(t) определится

известным образом:

 

 

Kt) =

dx(t)

(192)

е-ß t.

 

dt

 

Плотность вероятности сепарации отыщем путем исследо­ вания характера взаимодействия зерен с соломинами на ос­ нове механики столкновений [37].

При столкновении зерен с соломинами в процессе сепара­

ции зерна рассеиваются. Если

т 3

— масса зерна

(шз=

(0,03—0,04) - ІО-4 кг.сек21м), М

масса соломины

(М =

80


(0,6—0,8) ■ІО“ 1 кг.сек2Ім), т о как показал анализ (см. рабо­ ту автора «Об интенсивности выделения зерна из грубого вороха», «Вопросы механизации и электрификации сельско­ хозяйственного производства», изд. РГУ, 1969), при гпз<^М рассеивание носит изотропный характер.

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

В общем виде диффузионный процесс образуется с по­ мощью производящего дифференциального оператора [16]:

Dnf(z) = a[( z ) ^ ^ + a2( z ) ® ^

- a 3(z)F(z),

(]93j

где ai(z)

и а? (z) — коэффициенты

диффузии и сноса;

 

a3(z )— плотность вероятности обрыва:

 

Если

К (z,z0;t) плотность вероятности перехода

системы

из состояния z0 в состояние z за время t, то с помощью (193) можно записать:

дЦг, z0; t) = Dn).(z, z0; t).

(194)

dt

 

Уравнение (194) представляет собой параболическое диф­ ференциальное уравнение, решение которого имеет вид:

0 ai)t

(

exp

(z—z(,+ a2t)2

Mz,z„;t)—

 

4a,t

V 4ica,t

[

 

—exp

(z—z0+a.,t)2—4a,tzti

(195)

4a (t

 

 

 

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

М і)= а е “ Т1,

(196)

где

а,у — постоянные.

81