Файл: Липкович, Э. И. Процессы обмолота и сепарации в молотильных аппаратах зерноуборочных комбайнов (пособие для конструкторов зерноуборочных машин).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