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

книги из ГПНТБ / Алферов, С. А. Динамика зерноуборочного комбайна

.pdf
Скачиваний:
7
Добавлен:
22.10.2023
Размер:
11.7 Mб
Скачать

Уравнение ведущего момента рабочего органа

Внешняя скоростная харак­ теристика двигателя Д-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 Мп +

Мз)

Для выбранных значений постоянных величин внешние на­ грузки (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 увеличение момента инерции до 0,73 кГ сек2. С увеличением натяжения клиноременных передач до ар = 15ч-16 кГІсм2 величины Ах и Au значительно умень­ шаются.

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

Из расчетных осциллограмм было установлено, что с увели­

чением уменьшаются колебания угловых скоростей

х и и,

а крутящие моменты на валах барабана и контрпривода

имеют

при больших 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 ц/га

и различных вели­

чинах и Nde. Из графиков было

выявлено,

что с повышением

мощности двигателя до Nde = 60

л. с. вместо Nde = 50 л. с.

значительно уменьшаются колебания угловых скоростей рабо­ чих органов, больше чем при увеличении момента инерции моло­ тильного барабана Jх.

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

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

Г л а в а IV.

АНАЛИЗ ПЕРЕХОДНЫХ ПРОЦЕССОВ

 

В ПРИВОДЕ КОМБАЙНА В ЛИНЕЙНОЙ

 

ФОРМЕ

§ 13. ПОЛУЧЕНИЕ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ДЛЯ МОДЕЛИ ПРИВОДА И ОПРЕДЕЛЕНИЕ ПЕРЕДАТОЧНОЙ МАТРИЦЫ СИСТЕМЫ

Изучение динамики сложных машин (многосвязных систем), снабженных системами автоматического управления (САУ) и регулирования (САР), представляет собой важную задачу, ана­ литическое решение которой возможно пока только в линейной форме.

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

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

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

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