862
.pdf31
которое используют при анализе работы цилиндрических барабанов убо-
рочных машин [1.2], [1.3], [1.31].
Поскольку после прохождения точки 1 частицы начинают скользить относительно обода, причем угловая скорость снижается, и, как следствие, центробежные силы уменьшаются, то в некоторой точке 2 может произойти отрыв гранул от обода.
Когда угол поворота диска t находится между точками 1 и 2 ,т.е.1 t 2 , частицы скользят по касательной к ободу, который, в свою
очередь, вращается вокруг оси гранулятора. В результате такого перемещения частиц появляется сила Кориолиса (рис.1.14).
Рис. 1.14. Схема сил, действующих на частицу в фазе скольжения по ободу гранулятора
Если относительную угловую скорость частицы обозначить r, и спроектировать силы на подвижные оси координат и , то дифференциальные уравнения относительного движения можно составить на основании принципа Д`Аламбера.
|
|
m |
d 2 |
m r2 r mg cos N m 2 |
r 2m r r |
(1.32) |
|||||
|
|
|
|||||||||
|
|
|
dt 2 |
|
|
|
|
|
|||
|
d 2 |
|
|
|
d r |
|
|
|
|
(1.33) |
|
m dt 2 |
m r dt |
N tg f mg cos mg sin sin |
|||||||||
|
|||||||||||
|
Слагающая относительного ускорения |
|
d 2 |
вдоль оси |
направлена |
||||||
|
|
|
|||||||||
|
|
|
|
|
|
|
|
dt 2 |
|
вдоль радиуса к центру О и является относительным центростремительным |
||||
ускорением, которое выражается величиной 2r. |
|
|
||
r |
|
|
||
Слагающая относительного ускорения |
d 2 |
вдоль оси |
направлена |
|
dt 2 |
||||
|
|
|
по касательной к ободу и является тангенциальным относительным ускоре-
нием, выражаемым величиной r |
d |
2 |
|
r |
. |
||
|
|||
|
dt |
|
|
Из уравнения (1.32) можно найти величину силы N и поставить ее значение в (1.33).
N mg cos sin m 2r 2m r r m r2r ,
или
N m g cos sin ( r )2 r
32
Разность между переносной угловой скоростью диска и относительной скоростью r определяет величину абсолютной угловой скорости
гранул r .
В таком случае |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
N m (g cos sin 2 r) . |
|
|
(1.34) |
|||||||||||||||||||
После подстановки этого значения N в уравнение (1.33) и незначи- |
||||||||||||||||||||||||||||||||
тельных преобразований можно получить: |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
d r |
|
|
g |
|
sin( ) sin cos sin 2tg . |
(1.35) |
|||||||||||||||||||||||||
|
|
|
|
|
||||||||||||||||||||||||||||
|
dt |
|
|
r |
|
|
|
|
|
|
|
|
cos |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Если r |
и const , то |
d |
|
d r |
. |
|
|
|
||||||||||||||||||||||||
dt |
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|||||||
Так как |
d |
|
(где t) , то |
d |
|
|
d 2 |
. |
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
dt |
|
dt 2 |
|
|
|
|
|
|
|
||||
|
|
|
С учетом этого |
|
d |
r |
|
d 2 |
. |
|
|
|
|
|
|
|
|
(1.36) |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
dt 2 |
|
|
|
|
|
|
|
|
|
|
|
|||||
Учитывая значения (1.36), дифференциальное уравнение (1.35) мож- |
||||||||||||||||||||||||||||||||
но представить |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
d 2 |
|
|
g sin( ) sin cos sin |
|
d 2 |
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
tg , |
|
||
|
|
|
dt 2 |
|
|
|
|
|
|
|
|
|
cos |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|||||||||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
sin( ) sin cos sin 0 . |
|
|||||||||||||||||
|
|
|
( )2 tg |
g |
(1.37) |
|||||||||||||||||||||||||||
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r |
|
|
|
|
|
|
|
cos |
|
|
|
|
|
|
|
Выражение (1.37) представляет собой дифференциальное уравнение второго порядка, решить которое можно путем подставки значений z
и |
1 |
|
dz |
, после чего оно принимает вид: |
|
|
|
|
|
|
|||
|
2 d |
|
|
|||
|
|
|
|
z P( ) z q( ) , |
|
(1.38) |
где |
|
|
|
|
|
|
|
|
|
|
P( ) 2tg tg , |
|
|
|
|
|
q( ) 2g sin( ) sin cos sin /(r cos ) . |
|
||
|
Соотношение (1.38) является дифференциальным уравнением перво- |
|||||
го порядка, решением которого становится |
|
|
||||
|
z e P( ) d (c q( ) e P( ) d d ) |
. |
(1.39) |
После подставки значений P(α) и q(α), а также определения постоянной интегрирования С из начальных условий (когда α = α1, то ) получаем:
|
|
2 |
|
2g cos sin |
cos( ) |
cos |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
r |
cos |
|
|
|
|
|
2 |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. (1.40) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2g |
|
|
cos sin |
|
|
|
|
cos |
|
|
e( 1 ) tg |
||||||||
|
2 |
|
|
cos( |
|
) |
|
||||||||||||||
|
|
1 |
|
||||||||||||||||||
|
|
|
r |
|
cos |
|
2 |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
33
Если в уравнение (1.40) ввести значение β = 90 , то получим соотно-
шение
|
2g cos |
|
|
|
cos |
|
|
|
|
e( 1 ) tg |
, (1.41) |
2 |
|
cos( ) |
2 |
|
|
cos( |
1 |
) |
|
||
|
|
||||||||||
|
r cos |
|
|
|
cos |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
которое точно соответствует формуле, применяемой для определения действительной угловой скорости частиц, движущихся внутри цилиндрического решета [1.3], т.е. уравнение (1.41) выражает частый случай (1.40).
В момент отрыва частиц от обода в точке 2 нормальная составляющая (рис.1.14) равна нулю:
|
|
|
|
|
|
N m 2r mg sin sin( |
2 |
/ 2) 0 |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 (g / r) sin ( cos |
) |
|
|
|
|
|
|
|
(1.42) |
||||||
|
|
|
|
|
|
2 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
После подставки значения 2 с учетом того, |
что |
2r |
k , можно |
||||||||||||||||||
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
g |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
получить уравнение: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
2g cos sin |
|
|
|
|
|
|
cos |
|
|
|||||||||||
|
|
|
|
|
|
|
cos( 2 |
) |
|
|
|
|
|
|
|||||||
|
r |
|
|
cos |
2 |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
k |
|
|
cos sin |
|
|
|
|
|
|
|
cos |
|
|
|||||||
|
|
|
|
|
|
|
cos( 1 |
) |
|
|
|
|
|
|
|
(1.43) |
|||||
2 |
cos |
|
2 |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
e 2 1 tg gr sin cos 2 .
Из уравнения (1.43) можно определить значение фазы α2, при котором произойдет отрыв частицы от поверхности обода.
Угол α2 обычно больше 90 , что представляет определенные неудобства при математических операциях с тригонометрическими функциями. Положение точки 2 можно определить острым углом В2 (рис. 1.15), измеренным относительно горизонтальной оси.
В момент отрыва частицы от обода величина ее абсолютной скорости Vaбc=Ω·r , а при проекции ее на оси координат составляет:
Vx,абс r sin B2 , |
(1.44) |
Vy,абс r cos B2 . |
|
После отрыва от обода перемещение частицы по диску под действием сил инерции, трения о диск и силы тяжести происходит по некоторой траектории. В произвольный момент времени абсолютная скорость частицы может быть представлена в виде векторной суммы переносной (Vпер) и относительной (Vотн) скорости:
|
|
|
|
|
|
|
|
V |
абс Vпер Vотн |
(1.45) |
34
Рис.1.15. Схема сил, действующих
на частицу, и направление скоро-
сти ее движения по диску гранулятора
Переносная скорость Vпер направлена перпендикулярно к текущему значению радиуса R (рис.1.16), причем R х2 y2 , а
Рис. 1.16. Схема определения |
Рис. 1.17. Схема определения |
переносной скорости |
направления силы трения |
Vх,пер
Vу,пер
R cos R |
y |
y |
|
|
||
|
|
|
||||
R |
(1.46) |
|||||
|
|
x |
|
. |
||
|
|
|
|
|
||
R sin R |
R |
x |
|
|||
|
|
|
|
|
Величина относительной скорости может быть определена равен-
ством
__ |
__ |
__ |
V отн V абс Vпер |
или в проекциях на оси координат:
V х,отн V х,абс Vх,пер |
(1.47) |
|
|
Vy,отн Vy,абс Vy,пер |
|
С учетом значений составляющих переносной скорости
Vх,отн Vх,абс y |
, |
(1.48) |
|
Vу,отн Vу,абс x |
|||
|
|
полная величина относительной скорости будет равна
Vотн Vx2,отн Vy2,отн .
35
Сила трения частицы о диск направлена в сторону, противополож-
__
ную относительной скорости V отн (рис 1.17).
Если угол наклона относительной скорости к оси Х обозначить δ, то
|
cos |
Vx,отн |
, |
sin |
Vу,отн |
, |
|||
|
|
|
|||||||
|
|
|
Vотн |
|
|
|
Vотн |
||
а составляющие силы трения |
|
|
|
|
|
||||
F |
F |
cos F |
|
Vx,отн |
|
||||
|
|||||||||
x.тр |
тр |
|
|
тр |
|
Vотн |
|||
|
|
|
|
|
|
(1.49)
Fу,тр Fтр sin Fтр Vу,отн . Vотн
Величина силы трения частицы о диск с учетом наклона его к горизонту под углом β (рис.1.13)
Fтр fmg cos . |
(1.50) |
После определения величины и направлении сил, действующих на частицу, можно составить дифференциальные уравнения движения по диску с учетом принципа Д Аламбера:
|
m |
dV |
x |
fmg cos |
Vx,отн |
|
|
|
|
|
dt |
|
Vотн |
. |
(1.51) |
||||
|
|
|
|
|
|||||
|
dVy |
|
|
|
|
|
Vy,отн |
||
m |
mg sin fmg cos |
|
|
||||||
|
|
|
|
||||||
|
dt |
|
|
|
|
Vотн |
|
|
Вычисление приращений скоростей по этим зависимостям позволяет определить координаты движущейся частицы в любой момент времени и построить ее траекторию с помощью ЭВМ.
Построение математической модели движения частиц и разработка программ расчета параметров их траекторий с помощью ЭВМ обеспечивает подбор таких размеров и режима диска, при которых длина пути, пройденного гранулой, окажется максимальной, в связи с чем вероятность наклеивания на нее пылевидных частиц увеличивается.
На основании математической модели частиц по поверхности наклонного дискового гранулятора составлена компьютерная программа Drage определения основных показателей, характеризующих его работу.
Программа позволяет построить траектории движения гранул при различных условиях работы, определить границы зон относительного покоя и относительного скольжения частиц по ободу, длину траекторий, общий путь, проходимый частицами в абсолютном и относительном движении.
С помощью этой программы прежде всего изучена зависимость эффективного значения угла наклона диска от показателя режима и фрикционных свойств семян. Под эффективными значениями угла наклона диска понимаются те их значения, при которых траектории движения семян охватывают всю рабочую поверхность диска, так как в этом случае процесс
36
наслаивания защитно-питательной оболочки ускоряется, а возможность слипания гранул – снижается.
Результаты расчета приведены в табл.1.14
Таблица 1.14
Эффективные значения угла наклона диска при различных значениях показателя кинематического режима и угла трения семян
Угол |
|
|
Показатель кинематического режима |
|
|
||||
трения |
0,2 |
0,4 |
|
0,6 |
0,8 |
1,0 |
|
1,2 |
1,4 |
θ, град. |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
23 |
38 |
40 |
|
44 |
47 |
54 |
|
61 |
68 |
28 |
38 |
41 |
|
46 |
53 |
61 |
|
68 |
77 |
33 |
38 |
45 |
|
51 |
59 |
67 |
|
75 |
85 |
38 |
43 |
49 |
|
56 |
64 |
73 |
|
81 |
- |
Анализ данных табл.1.14 показывает, что ввиду изменчивости угла трения семян о рабочую поверхность дражиратора, особенно при смачивании гранул клеящим раствором (при этом угол трения возрастает до 36…38 ), у дражиратора необходимо изменять частоту вращения, либо угол наклона диска в пределах 38…85 .
Частота вращения диска может быть определена как
n |
30 |
|
|
k g |
|
, |
|
|
r |
||||||
|
|
|
|
|
где k - показатель кинематического режима дражиратора; r - радиус диска.
Расчетная длина пути, проходимого частицей по поверхности дражиратора приведена в табл.1.15 для случая, когда θ = 33 , а r = 0,31 м.
Таблица 1.15
Расчетная длина пути, проходимого частицей по поверхности дражиратора
Режим работы |
Полная, м |
Траектории по |
Участков с отно- |
|||||||
|
β |
диску, м |
сительным сколь |
|||||||
k |
|
м/с |
||||||||
|
|
м/с |
жением, мм/с |
|||||||
|
|
|
|
|
|
|||||
0,2 |
38 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
0,4 |
45 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
0,6 |
51 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
0,8 |
59 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
1,0 |
67 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
1,2 |
75 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
1,4 |
85 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
37
Данные табл.1.15 подтверждают, что одинаковой длины пути, проходимого частицей, можно достичь при различных значениях режима работы дражиратора и угла наклона диска, но продолжительность участков с относительным скольжением растет с увеличением частоты вращения (в указанных предела) и угла β.
Влияние угла наклона дражиратора на траекторию движения частиц представлено на рис.1.18. Данный расчет проведен для конкретного дражиратора (r–0,31, k – 0,18, θ – 33 ), на котором проведены экспериментальные исследования.
Рис.1.18. Влияние угла наклона дражиратора на траекторию движения частиц
Наблюдения за работой дражиратора в этом режиме подтвердили работоспособность установки в указанных пределах кинематического режима для условий лабораторной установки представлены на рис. 1.19.
Рис.1.19. Влияние режима работы дражиратора на траекторию движения частиц
Экспериментальная проверка подтвердила расчетный характер движения частиц в указанном диапазоне изменения показателей кинематического режима. Таким образом, программы расчета параметров траектории движения гранул в дражираторе могут быть использованы для определения их эффективных параметров и режимов работы при различных свойствах семян и состава наполнителя.
Характеристика семян моркови сорта «Золотой шпиль», дражированных на экспериментальной установке, приведена в табл.1.16.
|
|
|
|
|
|
|
|
|
|
|
38 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 1.16 |
||
|
|
|
|
|
|
|
|
Характеристика дражированных семян моркови |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Размеры, |
|
Количество |
Масса 1000 |
Плотность, |
|
|
|
|
|
|
|
|
с |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
семян в драже, |
|
|
|
|
|
Всхожесть, % |
Угол трения, град |
м/ |
|||||||||||
|
мм |
семян, г |
г/см3 |
|
, |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
% |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
скорость |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Фракция |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
группо- |
одиноч- |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
семян |
|
|
|
|
|
|
|
|
сушкидо |
сушкипосле |
сушкидо |
|
сушкипосле |
сушкидо |
сушкипосле |
месяцчерез |
сталь |
резина |
сталь |
|
резина |
Критическая |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
вой |
ный |
|
|||
|
M |
|
0 |
|
1 |
|
2 |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Основная |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
фракция |
4,79 |
0,51 |
2 |
|
95 |
|
3 |
- |
48,8 |
40,3 |
0,80 |
|
0,74 |
78 |
78 |
76 |
25 |
33 |
18 |
|
19 |
16,26 |
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Мелкая |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
фракция |
3,20 |
0,58 |
8 |
|
89 |
|
3 |
- |
39,2 |
31,6 |
0,86 |
|
0,78 |
77 |
76 |
76 |
20 |
27 |
18 |
|
23 |
13,35 |
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Крупная |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
фракция |
5,88 |
0,22 |
1 |
|
71 |
|
9 |
19 |
98,0 |
63,7 |
0,77 |
|
0,73 |
78 |
77 |
75 |
22 |
30 |
23 |
|
24 |
- |
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
38
39
Дифференциальные уравнения движения частиц по рабочей поверхности наклонного дискового гранулятора представляют собой детерминированную математическую модель технологического процесса дражирования, в то время когда все основные физико-механические свойства семян являются случайными величинами.
Для устранения этого недостатка разработана программа расчета параметров траектории движения частиц по ободу и диску дражиратора на ЭВМ, в которой коэффициент трения задается случайным образом перед анализом каждой траектории.
Поскольку практически все физико-механические свойства семян имеют нормальный закон распределения, то на ЭВМ необходимо формировать случайные числа в соответствии с данным законом распределения и определенными значениями математического ожидания и дисперсии.
В настоящее время все персональные ЭВМ снабжены программами, способными формировать случайные числа с равномерным распределением в интервале (0,1) по оператору RAND или RANDOM.
Для преобразования равномерно распределенных чисел в последовательность величин, имеющих нормальное распределение, моделируют условия, при которых оказываются справедливыми предельные теоремы теории вероятностей.
Известно, что в силу центральной предельной теоремы сумма большего числа случайных слагаемых (при выполнении весьма общих условий) имеет асимптотически нормальное распределение.
В таком случае коэффициент трения можно представить в виде суммы независимых одинаково распределенных величин:
R R1 R2 R3 ... Rk .
Математическое ожидание и среднее квадратическое отклонение такой суммы равны
|
|
|
|
|
a a* k; |
|
c |
* k |
|
c |
|
|
|
где а*, * - числовые характеристики Ri.
Для равномерно распределенных величин в диапазоне (0,1) a* 1/ 2 ; * 1/(2 3),
|
k |
|
|
|
1 |
|
|
; |
||
|
|
|
|
|
||||||
a |
и |
|
|
|
|
k |
. |
|||
|
с |
|
|
|||||||
c |
2 |
|
2 |
3 |
|
|
||||
|
|
|
|
|
Считают, что R приближается к нормальному распределению уже при небольших значениях k (от 4 до 10).
Чтобы пронормировать распределение этих чисел (получить распределение с числовыми характеристиками m = 0 и = 1), каждое из чисел R=ΣRi нужно преобразовать:
~ |
|
|
|
|
k |
|
|
2 |
3 |
|
|
||||
R |
(R |
) . |
|||||
k |
|
|
|||||
|
2 |
|
Для получения чисел Rз, образующих последовательность с другими заданными параметрами mз и з, необходимо провести еще одно преобразование.
40
|
|
~ |
|
|
R m |
з |
|
|
~ |
|
|
|
|||||
Поскольку R |
|
|
|
з |
|
, то R R |
|
|
m . |
||||||||
|
|
|
|
|
|
3 |
|||||||||||
|
|
|
|
|
|
|
|
з |
|
|
|
3 |
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
~ |
|
|
|
|
|
После подставки значения |
R можно получить случайные числа с |
||||||||||||||||
нормальным распределением и заданными параметрами, |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
||
|
|
|
|
|
|
|
|
Rз |
A RI |
B , |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
i 1 |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
где A 2 |
|
|
3 |
|
|
, B m |
|
A |
k |
. |
|
|
|
||||
з |
k |
з |
|
|
|
|
|||||||||||
|
|
|
|
|
|
2 |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
Для проверки степени приближения суммы равномернораспределенных чисел к нормальному закону сформированы два ряда распределения, содержащих по 10000 чисел (табл.1.16а).
В первом k = 5, а во втором k = 10.
Таблица 1.16а
Ряды распределения случайных чисел, сформированных для моделирования изменчивости угла трения семян о рабочую поверхность дражиратора при заданных числовых характеристиках θср = 33 ; = 1
|
Нижние |
k = 5 |
|
|
k = 10 |
|
|||
Номер |
границы |
|
|
|
|
|
|||
фактически |
|
ряд с нор- |
фактически |
ряд с нор- |
|||||
класса |
классов, |
сформир. |
|
мальным |
сформир. |
мальным |
|||
|
град. |
ряд |
|
распред. |
ряд |
|
распред. |
||
1 |
29,5 |
0 |
|
|
|
1 |
|
|
|
2 |
30,0 |
3 |
|
|
62,9 |
9 |
|
|
14,7 |
3 |
30,5 |
41 |
|
|
|
47 |
|
|
51,3 |
4 |
31,0 |
177 |
|
|
168,7 |
166 |
|
170,6 |
|
5 |
31,5 |
477 |
|
|
449,9 |
468 |
|
446,5 |
|
6 |
32,0 |
969 |
|
|
936,9 |
966 |
|
918,9 |
|
7 |
32,5 |
1513 |
|
|
1523,0 |
1495 |
|
1487,7 |
|
8 |
33,0 |
1961 |
|
|
1932,8 |
1826 |
|
1894,7 |
|
9 |
33,5 |
1838 |
|
|
1951,2 |
1858 |
|
1898,3 |
|
10 |
34,0 |
1473 |
|
|
1481,7 |
1477 |
|
1496,1 |
|
11 |
34,5 |
884 |
|
|
894,9 |
991 |
|
927,6 |
|
12 |
35,0 |
446 |
|
|
422,0 |
478 |
|
452,4 |
|
13 |
36,5 |
173 |
|
|
155,3 |
161 |
|
173,5 |
|
14 |
36,0 |
41 |
|
|
56,6 |
49 |
|
|
52,4 |
15 |
36,5 |
4 |
|
|
8 |
|
|
15,1 |
|
|
|
|
|
|
|||||
|
|
φср=32,98 |
φср=33,0039 |
||||||
|
|
σ=0,9945 |
|
σ=1,0104 |
|
||||
|
|
Σ 2=18,388 |
|
Σ 2=19,34 |
|
||||
|
|
Σ 2крит=16,92 |
Σ 2крит=19,68 |
||||||
Необходимо |
отметить, |
что |
только |
для |
второго |
ряда |
2 2крит , т.е. при статистическом моделировании случайных значений коэффициента трения необходимо, чтобы k > 10.