![](/user_photo/_userpic.png)
книги из ГПНТБ / Алферов, С. А. Динамика зерноуборочного комбайна
.pdfУравнение ведущего момента рабочего органа
Внешняя скоростная харак теристика двигателя Д-60 с регулятором:
ЯП
^ 3 = 2 9 - 183,73— г кГ'М
Mzз = 0
Момент на ведомом валу ва риатора привода ходовой части:
|
М,вар • |
|
75 |
Івар |
|
0,82 X |
|
|
|
|
I |
||||
|
|
|
hap “b |
|
|
||
X |
|
-9'32 l 1- i W |
« |
. |
п) |
||
|
|
|
|
|
|
|
кГ-м |
Мвар — |
|
75 |
Чар |
|
0,82 X |
||
|
|
1 |
|||||
|
|
|
Івар + |
|
|
||
X |
— 93,32і и I(- |
|
; |
_ 1 |
|||
— g |
|
\ z |
Мар к. п |
|
|||
|
|
|
|
|
|
|
кГ-м |
Момент, |
|
развиваемый |
перед |
||||
|
ними |
|
ведущими |
колесами |
|||
|
самоходного шасси: |
|
|
||||
|
|
Мп = 2080 X |
|
|
|||
X |
0,58 • |
0,058 |
|
|
I кГ-м |
||
|
— 28,9 |
||||||
|
|
|
|
|
|||
|
|
|
|
У |
|
/ |
|
|
Мп = —2080 X |
|
|
||||
X |
I 0,58--------- ----------- \ |
кГ-м, |
|||||
|
|
|
— 28,9— 0,9 |
|
|||
|
|
|
У |
|
|
I |
|
где ідіб. р = 7,6 -3,8 = 28,9
Момент, передаваемый на зад ние ведомые колеса при вы ключенной обгонной муфте:
М3 1 = 1195 X
X ( 0,58 — ------ÇU158 |
з кГ-м |
||
|
1,1 |
»•7,6 |
) |
|
|
||
А13. х = —1195 X |
|
||
X / 0,58 • |
|
0,058 |
I кГ-м |
|
■0,9 |
||
|
|
|
|
|
0-7,6 |
|
Условия существова- |
j |
ния уравнения |
[ |
116 < z < 183,73 Мсек
1116 > |
183,73 Мсек |
\—оо<^УИг<^0
У■
оо< —- івар1к. п
<Івар1к. п *С ■С + оо
■оо < — 28,9;
У
1 < — 28,9 < + оо
У
— оо < -----=-£- : 1
0-7,6 ^
1 < 1Г7Ж< + 00
Графическая
зависимость
Г\ т73 Z 116
Уравнение ведущего момента |
Условия существова- |
Графическая |
|||
рабочего органа |
|
|
ния уравнения |
зависимость |
|
Ведущий момент, развиваемый |
|
|
|
||
задними колесами при вклю |
|
|
|
||
ченной обгонной муфте: |
|
|
|
||
М, |
1195 X |
|
— ОО < — 30,4 1 |
|
|
X / 0,58 ■ |
0,058 |
кГ-м |
|
У |
|
|
|
|
|||
|
|
|
|
||
1,1 - — 30,4 |
|
|
|
||
|
У |
) |
|
|
\ \ С к о л ь > * е н и е |
М3 = —1195 X |
|
1 < |
30,4 < + ОО |
мп V |
|
X / 0,58 |
0,058 |
N кГ-м, |
|
|
|
|
30,4 — 0,9 , |
|
|
|
где ідіпл.р=7,6-4=30,4
Т а б л и ц а 4
Логические условия для выбора структуры динамической модели машины и системы уравнений
Логические условия изменения системы уравнений
Если z > w, то М 1 = + Мм. с Если z<^w, то Mj = — Мм. с
Если z = |
w и 1М м. с 1< 1М / 1, |
||
то М г = sign М / 1М м. с 1. |
где |
||
Mj — M u |
„ + М х X |
|
|
|
l UX\ |
|
|
X — |
>■- + Bww2 -г Jww |
|
|
1хЦ |
|
|
|
z = w |
и 1Мм. c ISs 1Mj\ |
|
|
Если 5 > 7 ,6 о |
и A43^ 0 , |
то |
|
|
|
3,8 |
|
Если |
< 5 ^ |
7,би, то |
|
|
s < |
JL |
|
|
6 ^ |
3,8 |
|
Кодовое
обозначение
модели
п и
1101
1011
1001
1111
1011
Физический смысл изменения структуры динамической
' модели машинет
Начало буксования муфты сцепления привода рабочих органов
Прекращение буксования муф ты сцепления привода ра бочих органов
Начало проскальзывания обгонной муфты привода заднего моста
Если S^> 7,6ѵ и А43^>0, то |
1101 |
Прекращение проскальзыва |
|
ния обгонной муфты при |
1001 |
вода заднего моста |
|
логические условия (табл. 4), то берут систему 1111 как удовлет воряющую этим случаям.
ЭЦВМ печатает для каждого шага кодовое обозначание си стемы уравнений, по которой производится вычисление этого
шага, и |
следовательно, кодовое обозначение структуры модели |
в данный |
момент. |
При переходе от одной системы уравнений к другой конечные |
результаты предыдущего шага являются начальными условиями
для |
следующей системы уравнений. |
|
|
||||
Начальные условия системы уравнений для каждого варианта |
|||||||
или |
группы |
вариантов даются |
в виде значений |
переменных х 0, |
|||
и о, |
w0, |
z0, |
у о, v0, |
s0 при t = |
0. |
ѵХх — ѵ%и и пропорцио |
|
Во всех системах уравнений |
величина |
||||||
нальна |
линейной |
скорости комбайна, |
которая |
была секунду |
тому назад. Эта' величина хранится в оперативной памяти ма шины.
В |
начале |
расчета можно |
принять |
ѵХх = |
ѵх |
— |
ѵ0. Для |
|||||||
различных значений Мг берут разные значения |
М ж. |
|
||||||||||||
В качестве ответов решаемой здесь прямой задачи динамики |
||||||||||||||
печатаются значения следующих величин: х, и, w, |
z, у, |
s X іпл. р, |
||||||||||||
V X і'б.р |
Х |
ід |
(где іпл_р = |
4; іб_р = |
3,8; |
ід = |
7,6). |
Кроме того, |
||||||
печатаются |
значения |
Мх, |
М и, |
Мвап х |
-—-— |
, |
М п —— г-------, |
|||||||
|
|
|
|
|
* |
“ |
вар |
ЬарЦвар |
|
п |
»6. p W n .K ’ |
|||
м з1 |
т |
|
— |
(где |
Цвар |
0,96, |
ч]п.к |
|
Лз. к |
|
0,95), |
Л4г, 5 |
||
|
Ід1пл. рЛз. |
к |
|
|
|
|
|
|
|
|
|
|
||
— 0,52 J |
vdt, |
t, код |
системы и |
к. п. д. г] |
ходовых |
колес, т. е. |
о
определяются по существу скорости отдельных частей несвобод ной системы и реакции связей.
Коэффициент полезного действия движителей или ходовых колес самоходного шасси вычисляется для каждого шага по сле дующим формулам:
при расчете шага по системе уравнений для типов структур 1111 и 1011, т. е. когда ведущими являются только передние колеса шасси,
„_ у (мп— [ М3, к I) _
'U. К— |
1 |
» |
У 3,8-7,6 Мп
при расчетах по системам уравнений для типов структур 1101 и 1001, т. е. когда ведущими являются передние и задние колеса шасси,
ѵ __________ (М п -р Ma)________
( 3,8-7,6 Мп + |
Мз) |
Для выбранных значений постоянных величин внешние на грузки Fх (t), Fu (t) и Fv (t) определяют по формулам:
Fx (0 = Q 1(1 + 0,7 sin 1,63of) X
X 0,278 + 0,044 0,0013л:] — 2,1;
Fu (t) = Q(f)xuvXu 0,044 (« 0,006 + 1,7);
(III.104)
t
2л 0,52 J V dt
Fv (t) = 472 + 335 sin----------------- |
. |
Дополнительными условиями для решения систем уравнений при частных случаях движения агрегата СШ-45 и КПН-2 являются следующие:
1.Момент, развиваемый муфтой сцепления при ее включении
сразличным темпом Т, находят по выражению
|
Мм.с = М М'Сmax( l - e ~ T ). |
|
(III.105) |
||||
При этом логические условия изменения системы уравнений |
|||||||
сохраняются, |
т. е. |
|
|
|
|
|
|
Мх = -\-ММ'С = + М М'СШХ \1 — е |
т), |
если г>ге>; |
|||||
|
|
|
|
|
|
|
(III.106) |
М1 = — Мм.с = — Мм.ст:іХ (і — е |
Г), |
если z < w . |
|||||
Величина |
является, таким образом, разрывной функцией, |
||||||
скачкообразно |
изменяющей |
в |
точке |
z = w |
свой |
знак. |
|
2. При расчете разгонов |
и |
выбегов Fx (і) |
= 0 |
и Fu (t) = 0, |
так как разгоны и выбеги рабочих органов и ходовой части агре
гата |
производятся |
без внешней нагрузки. |
3. |
При расчетах |
ряда вариантов разгона и выбега ходовой |
части, получаемых при изменении передаточного числа вари атора івар, необходимо соблюдать определенные логические пра вила, из-за. ограниченности диапазона изменения передаточного
числа вариатора: при |
івар = |
b — at, |
если |
igap ^ |
b, то |
Мвар = |
||||
= |
0; |
если Ь1 < івар < |
Ь, то Мвар вычисляется по |
соответствую |
||||||
щей |
формуле |
табл. 3; если івар^ Ь ъ |
то |
принимается |
івар = br |
|||||
и |
Мвар также |
вычисляется |
по табл. |
3, |
при этом b ■= 1,766, |
|||||
bx = |
0,82, а = |
0,4. |
|
|
|
|
|
|
|
|
и |
Подготовка систем уравнений для решений на ЭЦВМ «Урал-1» |
|||||||||
«Минск-2». Задачи |
по расчету динамики |
комбайна |
СК-4 и |
агрегата СШ-45 и КПН-2 решались в вычислительном центре Ростовского-на-Дону университета. При этом отыскивали реше ния систем из г дифференциальных уравнений первого порядка
при заданных начальных условиях. Из имеющихся численных методов для решения поставленных задач был выбран метод Рунге-Кутта.
Предварительно были проанализированы способ Эйлера и ме тод Адамса. При решении пробных вариантов оказалось, что спо
соб Эйлера |
и метод Адамса |
не могут дать решений этой задачи |
с требуемой |
точностью при |
шаге вычислений Я = 0,1 сек. |
Метод Рунге-Кутта дает погрешность порядка Я5. Для по лучения нового значения решения системы приходится подсчиты вать функцию в четырех дополнительных точках. Это прибли жает вычисляемое значение к истинному и способствует сохра нению точности.
Метод Рунге-Кутта при решении данных систем уравнений заключается в следующем. Наряду с системой дифференциальных
уравнений, |
аналогичной |
системе |
(II 1.82) |
|
|
|||||
|
|
dzx |
fi(x, |
Zu |
z2, |
|
|
|
|
|
|
|
dx |
|
• - I z |
r ) \ |
|
||||
|
|
|
|
|
|
|
|
|
||
|
|
dz2 |
f2 (x, zlt |
z2, |
|
• •> zr); j |
(III.107) |
|||
|
|
dx |
|
|||||||
|
|
%T = fr(X' |
21> г2, |
|
■• > zr)> |
|
||||
интегрировали также уравнения ^ - = |
1 и |
= zm, в которых |
||||||||
vm— добавочные переменные, |
представляющие собой интегралы |
|||||||||
вида Jt |
zdt. |
|
|
|
|
|
|
|
|
|
о |
|
|
|
1 точке вычисляли по формулам: |
||||||
Значения zt и х в k + |
||||||||||
|
zi, k+i = Zk 4- -g- kuh + |
h A |
+ |
"з- h i h -f — Я4(Я; |
||||||
|
|
|
|
xk+i = xk + |
Я, |
|
|
|||
где |
|
|
|
|
|
|
|
|
|
|
|
|
&1І |
fi |
|
%iky |
|
• |
■•> |
zrk)\ |
|
Ki — ft (xk -f- 4" h>zik + |
-j- |
ku h, |
. . -, |
zrk + 4" |
klrh |
|||||
h i ~ |
fi (^Xk H |
2~ |
Д* H |
2~ hih> |
■■• ) |
zrkA— |
hrh |
kn = fi (xk + Я, Zik + k31h, . . ., zrk + k3rh).
Рунге-Кутта и приведенных |
логических |
условий. |
Оптималь |
|
ные шаги |
были h = 0,1 сек при расчете |
на машине |
«Урал-1» и |
|
h = 0,05 |
сек— на машине |
«Минск-2». |
|
|
Транспортное запаздывание на величину т при решении си стемы уравнений на ЭЦВМ учитывалось следующим образом: запоминалось значение необходимой переменной в течение 1,1 сек в оперативной памяти машины, а затем содержимое десяти по следних ячеек сдвигалось после каждого шага. Таким образом,
машина «знала» значение |
переменной секунду, полсекунды и |
||
т. д. назад. |
|
|
|
Более трудной задачей |
являлось |
программирование функций |
|
М; |
|
в, |
|
-А,— |
2У+1 |
||
|
|
Сі~ |
|
|
|
|
г.1 |
выражающих моменты, |
развиваемые |
фрикционными передачами |
и дрИз положительной и отрицательной области возможных зна
чений Mj выбирали лишь те, которые соответствовали логическим условиям (см. табл. 1 и 3) существования уравнения и, следова тельно, физическому смыслу.
Для характеристики двигателя логические условия опреде ляют Мг = 0 при максимальной угловой скорости. В ЭЦВМ были заложены логические правила и ограничения, определяю щие работу комбайна в различных условиях.
Численные значения переменных z,-, zjt Mj и др. печатались на таблограмме через каждый шаг, а затем вычерчивались как функции времени t, давая своеобразные расчетные осциллограммы работы комбайна в заданных условиях.
§ 12. РЕЗУЛЬТАТЫ РАСЧЕТОВ ДИНАМИКИ КОМБАЙНОВ
Комбайн СК-4 — система с постоянной структурой. Варианты, решенные на ЭЦВМ «Урал-1», были получены при условии, что предварительное натяжение в ветвях клиноременных передач привода равно двум значениям: 1) ар = 10ч-11 кГІсм2 и 2) а — = 15ч-16 кГІсм2.
На рис. 43, а и б даны аналогичные полученным в гл. II рас четные графики изменения угловых скоростей Ах и Au в зависи мости от частоты X внешнего возмущения Fх (t) при различных значениях величин Jx и ар. Следует отметить, что с увеличением момента инерции в пределах Jx = 0,62 ч-0,73 кГ м сек2 вели чина Ах уменьшается даже при уменьшенном натяжении ремней ор = 10 кГ/'см2. На основании этого можно считать оправданным для комбайна СК-4 увеличение момента инерции Jх до 0,73 кГ -м сек2. С увеличением натяжения клиноременных передач до ар = 15ч-16 кГІсм2 величины Ах и Au значительно умень шаются.
Таким образом, увеличение натяжения клиноременных пере дач в данных условиях имеет большее значение, чем повышение момента инерции молотильного барабана. Соответствующие экс периментальные значения Ах и Аи близки к расчетным.
Из расчетных осциллограмм было установлено, что с увели
чением Jх уменьшаются колебания угловых скоростей |
х и и, |
а крутящие моменты на валах барабана и контрпривода |
имеют |
при больших Jx меньшие средние значения и меньшие колебания. При неравномерной подаче хлебной массы до qcp = 6 кг/сек обороты контрпривода и барабана падают главным образом из-за
Рис. 43. Расчетные графики изменения угловых скоростей Ах и Дм для комбайна СК-4 в зависимости от частоты К поступления порций хлебной массы при различных Jx в кГ-м -сек3
скольжения ременной передачи к контрприводу. Установлено, что при натяжении ремней ор = 10 кГІсм2 эта передача пере гружена при всех рабочих режимах до предела, соответствующего началу вредного скольжения. Расчетами на ЭЦВМ установлено влияние колебаний Fy (t) и Fх (t) на угловые скорости х, и, у, z.
Изменения угловых скоростей при разгоне и выбеге рабочих органов комбайна СК-4 при ор = 10 кГ/см2 и Jх = 0,62 K F X Хм сек2 показаны на рис. 44, а и б. Расчеты и эксперименты дают следующие значения времени разгона и выбега рабочих ор ганов комбайна СК-4:
|
Расчетные |
Эксперимен |
|
|
тальные |
||
Ор в кГ/см2 ....................................................... |
10 |
16 |
15 |
і,раз в сек ....................................................... |
8,5 |
6 |
5.5 |
іШб в сек ............................................................... |
10,5 |
18,5 |
19.5 |
Агрегат СШ-45 и КПН-2 — система с переменной структурой. Для данного агрегата было выполнено 27 различных вариантов расчетов. По расчетным данным доля тягового усилия задних колес шасси СПІ-45 при включенной обгонной муфте (ôeK= 5%) составляет 30% от усилия, развиваемого передними колесами, что примерно совпадает с данными испытаний этого агрегата в 1960 г. в ВИСХОМе.
При расчетах на ЭЦВМ выявлено, что продолжительность выключенного состояния обгонной муфты составляет около 40%
146
при движении шасси по полю с волнистым микрорельефом. Ми нимальный к. п. д. ходовых колес г\х к = 0,92 соответствует ра боте агрегата с включенным задним мостом (типы структур 1001 или 1101), а максимальный г\х к = 0,96— работе с выключенным (типы структур 1011 или 111І).
При разгоне ходовой части агрегата на I и II передачах мо мент Мх на привод молотильного барабана становится отрица тельным, а угловые скорости рабочих органов снижаются, осо бенно при разгоне на II передаче. Эта особенность была отмечена на осциллограммах при лабораторно-полевых испытаниях зерно уборочных комбайнов.
Рис. 44. Кривые разгона и выбега рабочих органов комбайна СК-4 при расчете на ЭЦВМ (Jx = 0,62 кГ-мсек2 и сгр = 10 кГ/см2)
Расчетные значения времени разгона и выбега ходовой части агрегата на I и II передачах совпадают с экспериментальными данными:
|
|
|
|
|
|
Расчетные |
Экспериментальные |
|
|||
Передача |
................................ |
|
|
|
I |
|
II |
I |
И |
|
|
|
|
|
|
|
с вариа- |
без |
с вариа- |
без |
|
||
tpaa в |
сек |
|
|
|
|
тором |
вариатора |
тором |
вариатора |
||
........................... |
|
|
|
3 |
|
8 |
3 |
8,4 |
|
||
івыб в |
сек ........................... |
|
|
|
4,5 |
|
— |
4,5 |
— |
|
|
Разгон и выбег рабочих органов комбайна КПН-2 получены |
|||||||||||
при |
различных |
темпах включения |
дисковой |
муфты |
сцепления |
||||||
Т = |
0 |
(рис. |
45, |
а); |
Т = |
0,7 |
сек |
(рис. 45, |
б) и Т = 1,5 |
сек |
|
(рис. 45, в). |
При |
анализе |
разгонов |
рабочих |
органов |
можно |
от |
||||
метить |
следующие |
особенности: |
|
|
|
|
1. С увеличением темпа включения Т минимальная угловая скорость коленчатого вала двигателя незначительно увеличи вается, а затем опять уменьшается. Следовательно, существует оптимальная величина темпа включения, равная по расчетам
Т0,7 сек.
2.При расчетах процесса разгона с Т > 0 отмечено возник новение колебаний крутящего момента в муфте сцепления и ко
лебаний моментов М и и М х. Эти автоколебания, по существу, можно объяснить нелинейностью характеристики муфты сцепле ния. Подобные колебательные процессы в муфте сцепления от
мечены |
и при экспериментах [17]. |
3. |
Результаты расчетов и экспериментальные данные при раз |
гоне и выбеге рабочих органов комбайна КПН-2 близко совпа
дают: |
расчетные |
(на ЭЦВМ |
«Минск-2») — tpa3 я» 6-ь6,5 сек, |
івыб ^ |
20 сек\ экспериментальные — tpa3 я» 6 сек, іеыб я« 18 сек. |
||
Были получены |
расчетные |
зависимости изменения угловых |
скоростей Ах, Au, Az и Ау при работе агрегата на I передаче
z,w,x,u z,w,x,u
с івар = 0,822 при урожайности Q = 20 ц/га |
и различных вели |
|
чинах Jх и Nde. Из графиков было |
выявлено, |
что с повышением |
мощности двигателя до Nde = 60 |
л. с. вместо Nde = 50 л. с. |
значительно уменьшаются колебания угловых скоростей рабо чих органов, больше чем при увеличении момента инерции моло тильного барабана Jх.
Таким образом, на основе метода Лагранжа с неопределенными множителями в данной главе были получены системы дифферен циальных нелинейных уравнений, описывающих большие движе ния сложных сельскохозяйственных машин с голономными и неголономными связями в приводе.
При рассмотрении несвободных механических систем с неголономными связями необходимо число уравнений, равное числу обобщенных координат, определяющих поведение системы. При этом используются экспериментальные тяговые характеристики фрикционных передач (неголономных связей), однозначно опреде ляющих неизвестные множители Лагранжа.
Г л а в а IV. |
АНАЛИЗ ПЕРЕХОДНЫХ ПРОЦЕССОВ |
|
В ПРИВОДЕ КОМБАЙНА В ЛИНЕЙНОЙ |
|
ФОРМЕ |
§ 13. ПОЛУЧЕНИЕ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ДЛЯ МОДЕЛИ ПРИВОДА И ОПРЕДЕЛЕНИЕ ПЕРЕДАТОЧНОЙ МАТРИЦЫ СИСТЕМЫ
Изучение динамики сложных машин (многосвязных систем), снабженных системами автоматического управления (САУ) и регулирования (САР), представляет собой важную задачу, ана литическое решение которой возможно пока только в линейной форме.
Аналитические методы изучения динамики сложных сельско хозяйственных машин позволяют наметить общие пути решения таких вопросов, как динамическая точность поддержания выход ных координат систем, устойчивость их движения, поведение при случайных внешних возмущениях, рациональная структура си стем и др. Для механических систем с неголономными (дифферен циальными) связями использование аналитических методов было затруднено из-за сложности получения дифференциальных урав нений, отражающих большие движения таких систем, и необхо димости, как правило, последующей линеаризации этих уравне ний.
Основная трудность при анализе сложных сельскохозяйствен ных машин с неголономными связями — это разработка идеали зированной модели, воспроизводящей действительную систему. Необходима предварительная аналитическая или эксперимен тальная оценка различных факторов и свойств системы для вы бора только немногих главных факторов, определяющих поведе ние системы в заданных условиях. Движение сельскохозяйствен ной машины при этом можно представить как движение модели с ограниченным числом степеней свободы и, следовательно, обоб щенных выходных координат.
Таким образом, независимо от метода получения и вида перво начальных дифференциальных уравнений приходим к необхо димости создания линейной модели сложной машины, имеющей г выходов и, как правило, г входов. При этом динамика модели описывается г дифференциальными линейными уравнениями. Си стема дифференциальных уравнений позволяет найти передаточ