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

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

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

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

Как было указано раньше, системы из и г уравнений типа (III.81) и (III.82) описывают большие движения роторных орга­ нов сложных сельскохозяйственных машин, в частности, зерно­ уборочных комбайнов. Величины а уи Zj (/ = 1, 2, . . ., г) являются обобщенными координатами или выходами динамической системы. При анализе переходных процессов в приводах комбайнов пред­ ставляют интерес в основном изменения рабочих скоростей zjt определяющих качество технологического процесса машины. Поэтому, если считать для моделей приводов комбайна приведен­

ные

моменты

инерции

постоянными, т. е. J f (а) я« const (у =

= 1, 2, . . ., г),

то, как было показано, при анализе дина­

мики

комбайна

можно

исходить из уравнений (II 1.82).

Для аналитического исследования переходных процессов в ме­ ханической системе, описываемой уравнениями (III.82), эти урав­ нения необходимо линеаризовать, исходя из условия, что выход­ ные координаты zu z2, ■■-, zr получают малые приращения Azlt Az2, . . ., Azr, вызванные действующими на машину внешними возмущениями и регулирующими факторами со стороны органов управления. Для получения линеаризованных уравнений раз­ ложим все функции fj (/ = 1 , 2 , . . . , г) уравнений (III.82) в ряды по малым приращениям и учтем только линейные элементы при­ ращений. Таким образом, получим

d Дzx

 

dfi_

A^iT

 

àfi

AZ2 -f-

 

LAzr -|

dfi

Af;

~ d t ~

 

дгг

 

dz-

 

dt

d A22

 

dfx

Azx ■

 

df*

Az.

 

 

dh

Af;

dt

 

dzx

 

 

dz,

 

 

 

dt

(IV. 1)

d &zr

_

dfr

д I

 

dfr

д

 

Azr -f-

dt

At.

~ З Г

-

dz.

AZl +

 

dzü

dzr r 1

 

В системе уравнений

(IV. 1) сделаём следующие подстановки:

 

Æ

. -

n

A L . — п

àfi

_

 

 

 

dzx

n ’

dz2

12T ‘ ’ *’

dzr

 

 

 

 

 

Aгі =х[,

Az2 = x’2> . . ., Azr = xr\

 

 

^ M

= h (0,

 

А/ = h (0.

• • -,

- f 1 A* = M0>

тогда система линеаризованных уравнений (IV. 1) будет

dx1

 

~ d f

~

« і Л

“ Ь

«12*2 +

• ■ ■ +

n \r x r +

h

( 0 ;

 

 

dx2

 

 

 

« 22*2 -+------- b n.2rx'r +

 

 

 

 

~df

=

n 2i*l

+

/2 (0;

( I V . 2)

 

<ІЛГГ

 

ППХ'і +

 

----- Ь n rrX'r + fr (0-

 

 

“ЗГ =

«г2*:2 +

 

Видоизменяя систему

(IV.2), получаем

 

 

 

 

(Pi -

«п) *1 -

«i2*2----------- « іЛ =

/і (0;

 

 

«21*1

~Ь ( р і

 

« 2 2 ) *2

'

« 2 г* г

=

/2 ( 0 і

(ІѴ.З)

 

 

 

 

 

 

 

 

 

 

-

« ,1 * і - «г2*2 + ■■• +

( Р і -

« гг) *г =

fr ( 0 ,

 

где

 

P1 = 4t

— символ дифференцирования;

fi (0. /2

(0.

• • -, fr (0 — внешние

возмущения,

ограничен­

 

 

 

 

 

ные по модулю.

 

 

 

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

уравнений:

 

 

*і (Р)

fl (р)

 

А ■ *2(Р) =

/2 (р)

(IV.4)

*г(Р)

М р )

 

ѵгде

А — передаточная

матрица,

опреде­

 

ляющая

динамические

свойства

 

модели

механической системы;

х г (р), х 2 (р), . . ., хг (р) — изображения

выходных

коорди­

 

нат х[,

х'ѵ

. . ., х'г

по

Лапласу;

(р)> /2 (р)I • • -,‘fr (р) — изображения

внешних

возмуще­

 

ний / х (0,

/ 2 (0,

• ■

fr (0 по

 

Лапласу;

 

 

 

 

р— комплексная переменная преобра- - зования.

Передаточная матрица Л с рангом г согласно системе уравне­ ний (ІѴ.З) имеет вид

(р — «п)

«12

■ ■ — «1Г

 

-- «21

(Р — «22)

• •

— «2Г

(IV.5)

— «г1

пп

(Р — Пгг)

 

Некоторые из коэффициентов п матрицы А могут быть равны нулю.

Обозначая матрицы-столбцы изображений выходных коорди­

нат Xj (р) (/ =

1 , 2 , . . . ,

г) через

X , а

внешних

возмущений

fi (р) ( / = 1 , 2 ,

. . г)

через F,

получаем

систему

(IV.4) в виде

 

 

A X

= F.

 

(IV.6)

На основании выражения (IV.6) изображения выходных ко­ ординат по Лапласу при известных изображениях внешних воз­ мущений будут

X = A~1-F, (IV.7)

где А~г — матрица обратного преобразования для динамической модели машины.

Для определения изображения какой-либо переменной, на­

пример х х (р) при f 1 (р) =р 0 и

/ 2 (р)

=

/3 (р) =

• • ■=/ , (р) = О,

запишем

 

 

 

 

* і (Р) =

м и (Р)

fi

(P),

(IV.8)

D(P)

где М 1г(р) — минор определителя D (р) передаточной матрицы А, получающийся при вычеркивании первого столбца и первой строки.

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

ний

Д (р),

/ 2 (р),

. . ., fr (р):

 

 

 

 

 

*і (Р)

Мц(р)D(P)

h (P) +

h (P) + ■• •

+

D(p) frip), (IV.9)

где

M д (p) — миноры определителя

D (p),

получающиеся при

 

 

вычеркивании

первого

столбца

и соответственно

 

 

/ =

1, 2, . . .,

г строки.

 

 

Из выражения (IV.9) можно определить передаточные функ­

ции Wjx (р)

для х г (р) при

действии

возмущений

(р).

г »(р) =

Т ^ Г .

=

 

К .(/>) =

М п (Р)

 

D(P)

 

 

 

 

 

(IV. 10)

Матрица обратного преобразования

 

 

 

Wn(P)

Wu (p) .

Wn (p)

 

 

А ~1= Wn {p)

1Ѵ22 (p) .

■Wn(p)

(IV.11)

 

Wir (P)

Wtr(p).

w rr(p)

 

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

Амплитудно-частотная характеристика А (со), показывающая изменение модуля одной из выходных координат в зависимости

от частоты внешнего возмущения, определяется так:

 

А (ю) = I W (ш) I = V I U (®) I2 + I V О'®) I2,

(IV. 12)

фазо-частотная характеристика ср (со), определяющая фазу вы­

ходного

колебания,

 

Ѵ(т)

 

 

Ф (со) =

arctg

(IV. 13)

 

U(a>) ’

где

I W (iсо) I — модуль

амплитудно-фазовой

характери­

 

стики;

 

 

 

U (со) и V (і'со) — действительная и мнимая частотные харак­ теристики.

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

Общепринятыми динамическими характеристиками в области действительного переменного t являются реакции системы на еди­ ничную импульсную функцию (дельта-функцию или функцию

Дирака)

6 (t)- и ступенчатую функцию f (t) = 1.

 

 

 

Если

изображение переходного процесса в операторной форме

будет L

I X,- (t) I =

W (p) fj (p),

то для f; (t) = 1

имеем

(p)

=

= L I fj

(t) I = - y

и для

(t)

= fi (t) имеем ff (p)

= L\f,

(Q| =

1.

Соответствующие изображения для выходных координат будут

L \ Xl{t)\ = W{p)±- и L I kj (t)\= W {р)\,

где Xj- (t) — переходный процесс при единичном ступенчатом воз­ действии или переходная функция;

kj (t) — импульсная переходная функция,

Если считать, что для устойчивых систем переходные процессы при нулевых начальных условиях затухают достаточно быстро, то поведение системы (реакция) определяется вынужденным дви­ жением [см. выражения (IV.8) и (IV.9)].

Рассмотрим получение линейных уравнений движения дина­ мической системы и анализ их применительно к четырехмассовой модели привода комбайна СК-4. Для этого разлагаем члены урав­ нений (II 1.98) в ряды и ограничиваемся только линейными эле­ ментами приращений. При линеаризации пренебрегаем силами сухого трения в динамической системе, считая их постоянными. Тогда получим следующую систему линейных дифференциальных уравнений:

 

 

d Ах

-

(

дМх

 

дМс. g. х \

Ах 4

 

 

 

 

dt

дх

 

 

дх

)

 

 

 

 

дМхда

Аи

 

 

 

ХУ

 

 

 

d Au

 

дМи

 

дМх

44’^+^

 

 

 

 

 

 

 

дМс. а. и

Au

 

Ju dt

 

ди

 

ди

 

‘.л

ди

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дМг

 

 

+'

dz

Az -

dFu(t)

д..

I f .

 

дх

 

 

 

 

 

 

 

 

r

dAy

( àMy

 

 

dM,С- S. у

) b y

-t-

(IV. 14)

dt~

Vdu

K- n

ду

 

 

 

I

дМу .

 

»

 

 

дМу

» .

г

г-.

 

+

~дГ Ік. пА2 +

4 Г 7 Ік. П^еар +

Ту,

 

 

г

d Az

 

і

дМг

 

дМи

1

 

 

 

 

 

2

dt

і

 

dz

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz

Ік. п

 

dMc в. Z-'j

Az

dMu

 

■Au

 

 

ди

 

 

 

1у<у

 

 

dz

 

 

 

VI«

 

 

 

 

 

 

 

 

 

 

 

дМу

Ік. п

ілу

 

 

дМу

Ік. п

Ai.

 

 

 

 

ду

V *

 

 

 

dieap

1у%

 

вару

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дМх

 

 

 

и^х^х

 

 

 

 

 

 

 

 

дх

 

 

 

(и — xixГ

 

 

 

 

 

 

дМсв. X

 

д ( V 2) _

- 2 В'

 

 

 

 

 

дх

 

 

 

дх

 

 

 

 

 

 

 

 

дМX

 

 

&XXix

 

 

 

 

 

 

 

 

 

ди

 

 

ХІх)2 У

 

 

 

 

dFx (t)

[если в формуле (111.95) считать 1

a sin Xt

ày

0];

 

 

 

 

 

 

 

 

 

 

 

 

* 1, a ß'p^x

 

 

dFx (t)

 

 

 

 

 

 

 

 

 

 

 

Д/ = f

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

дМи _

гіиВи .

 

дМс. в. и

_

' „2 !

 

 

 

 

d j Buu

=

 

 

ди

(Z

иіи)2

 

 

ди

 

 

du

 

 

 

 

 

 

 

 

 

 

 

 

дМи

 

BgUiu .

 

 

 

 

 

 

 

dz

 

 

(г — Шц)2

 

 

 

■= Q$’ku [если в формуле (III.96) считать

р„ы 0]

 

 

dFu (0

д f _ г' .

 

 

 

 

‘W

 

 

 

 

0/

1 ~~ I

 

 

 

 

 

 

 

 

 

 

 

 

(Z — УІварС)г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[здесь ІварС

^

itu/c^K. ni

ішк

А

(СМ. рИС.

11)],

 

Q

 

 

 

 

ам,.

 

 

с^еарВуУ

 

 

 

 

 

 

 

ÔZ

 

(z — у/варС)а

 

 

 

 

<3-'We. в.ÿ

 

ö (ß> 2)

=

2^г/;

 

 

 

 

 

ду

 

 

 

 

ду

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дМ,

 

 

л

 

 

ß ^2

 

 

 

 

 

Д_ =

 

 

 

 

 

 

•Эг'вар

 

ÿ

(z — y iea p c f

 

 

 

 

 

 

 

dßy (О

 

 

 

 

 

 

 

 

 

 

 

 

d/

д* =

f„;

 

 

 

 

 

 

дМг

 

 

 

ß,

 

 

 

 

 

 

 

dz

-

 

(С — г)2

 

 

 

 

амс

 

 

 

ô(ô>2)

=

2ßzz;

 

 

 

 

 

dz

 

 

 

 

dz

 

 

 

 

 

 

Me. в. *, Me. в. и, Mc. e. y, Mc. в. z — моменты, учитывающие

со­

 

 

 

 

 

 

 

 

противление

 

воздуха

для

 

 

 

 

 

 

 

 

валов X, и, у, г.

 

Для упрощения уравнений вводим следующие относительные

или безразмерные величины:

 

 

 

 

 

Дг

 

 

Дх..

и

,

 

Ди

У

=

Ду_.

 

 

X ’

 

= —

ъ

Z

 

 

и принимаем

 

 

 

U

 

 

 

Y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

У*Х = 1;

JJJ

=

1;

JyY =

1; JzZ =

1,

 

 

где X, U, Y, Z — заранее выбранные значения выходных коорди­ нат, относительно которых сравнивают вели­ чины Ах, Au, Ау, Az.

Вводим также следующие

обозначения:

/дМх

дмс.в. х

-)х = пі;

V дх

дх

 

дМх T-,

=

П9:

 

дих и

дМ_и

дМх

 

дМп

ди

ди ; г,"

 

ди

 

ух U

 

 

 

дМи т,

 

« e;

 

ô 2 ^ =

dFx (і)

Y _п •

 

= п 3,

ду

 

 

 

 

 

 

дМх

T

-

'

дх

 

dFu (t)

-,

п7;

ду

 

Y =

 

 

 

1

Г-».

 

Ф * 3

 

дМу ■

J _

дМг

дМи

1

dz

dz

і ух"

 

 

um

 

 

àMu

 

 

ди

 

 

дМу

 

 

ду

 

 

дМи

 

 

ді,еар

11dMc. в ,ц

) Ÿ = n8

ду

 

дМу .

дівар

 

 

дМи ік. п

дМп

‘lit

dz

dz

■U — п12 t 1иПи

ік-п Ÿ = п13 t

1У%

 

«к. /

Пи-

Іу^у

 

После подстановки в систему (IV. 14) принятых обозначений получим

 

t l \ ) Х —

4 ~ П3У = fXI

 

п5х +

+ Л4) U -f п7у

n6z

= f u\

 

 

4“ ^в) У

^9^

 

 

(IV. 15)

 

~ fу ~\~flU)ieap>

 

4-

ПізУ 4~ (р 4 - пп ) г

= fz Пиіварі

 

где р' = --------символ дифференцирования;

 

 

 

f'z — внешнее

возмущение

на

двигатель, связанное

с изменениями барометрического давления,

тем­

пературы среды,

качества

топлива и

т. п.

(часто принимают fz =

0,

а івар = Аівар)-

 

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

156

чается в подстановке вместо р' и функций х', и', z', f комплексной переменной р и соответствующих изображений функций х, и, у, z, f. Для изображений функций система (IV. 15) в матричной форме будет

X

 

 

и

fu

(IV.16)

 

 

Уfу "Ь «lo^'sap

Z

fz

«14leap

где А — передаточная матрица модели привода комбайна СК-4, имеющая вид

Р + Пі

 

«3

0

 

Пь

P + «4

«7

— «6

(IV. 17)

0

0

P + «8

—«9

 

0

«12

«13

P +«11

 

Система уравнений (IV. 16), устанавливающая связь между зависимыми или регулируемыми выходными координатами х, и, у, z и внешними или регулирующими факторами івар, fx, fu, fy, fz и характеризующая динамические свойства машины при работе на поле, называется системой уравнений, определяющей движение комбайна как управляемого процесса. Система (IV. 16) действи­ тельна при работе дизеля со статическим малоинерционным всережимным регулятором числа оборотов, применяемым на самоход­ ных зерноуборочных комбайнах СК-3, СК-4 и др. При работе двигателя комбайна с таким регулятором система линейных урав­ нений в развернутом виде будет

 

( р

+ п х )

X — «2« - f п 3 у =

Д;

 

 

«5* +

(р +

п4) и +

«7г/ — п в г =

 

 

(Р^~па)у

 

n9z =

fy 4" nwhap\

(IV. 18)

 

«12« + «13?/ + (P +

« ІО Z —

k y =

f l — Пиівар,

 

 

kiz + (Т2Р2 + Tip +

1) y = фФ (р),

 

где

k — коэффициент

усиления в уравнении двигателя от

 

перемещения

рейки;

 

 

 

 

Y— изображение

выходной

координаты

(движения

 

рейки)

всережимного

регулятора;

 

 

k x — постоянный

 

коэффициент воздействия от двига-

 

теля на регулятор;

 

 

 

абсолютной

7+

Tl-— постоянные времени,

весьма малые по

 

величине;

 

 

 

 

 

 

ф— регулирующее воздействие со - стороны ручного управления подачей топлива;

Ф(р) — оператор управляющего воздействия.

Последнее выражение в системе (IV. 18) представляет собой уравнение движения статического регулятора оборотов двигателя.

Если считать Т\ ^ 0 и Т 1 «=* 0, то уравнение движения

регуля­

тора будет

 

Y = ftjz + Ч>Ф (р).

(IV. 19)

Подставляя значение у из выражения (IV. 19) в уравнение двигателя в системе (IV. 18), получим для члена ky следующее значение:

ky = k [kxz + г|іФ (р) ].

Подставив, в свою очередь, это значение ky в систему (IV. 18)

и обозначив

ли — kk\ = ли

и

fz +

k\\іФ (р) = f2,

получим

систему из четырех линейных

уравнений

(IV. 16).

(IV. 16)

Постоянство

коэффициентов

л,-

в системе уравнений

сохраняется только при данном режиме нагружения комбайна. При изменении режима работы или характеристики отдельных параметров необходимо снова вычислять коэффициенты линеари­ зованных дифференциальных уравнений. При значительных откло­ нениях величин X, и, у, z, превышающих 0,1 номинального зна­ чения, линеаризованные уравнения становятся недействитель­ ными. В этих случаях необходимо пользоваться приведенными раньше нелинейными уравнениями и решать их на ЭЦВМ. В табл. 5 даны приближенные значения коэффициентов лудля двух нагрузочных режимов комбайна СК-4 при уборке пшеницы

урожайностью Q — 26 ц/га,

работе на II

передаче с ік, „ = 2,79

и

івар = 0,94

и

натяжении клиноременных

передач

ар =

=

15 кГ/см2. Нагрузочные режимы для комбайна СК-4 определены

 

 

 

 

 

 

 

Т а б л и ц а 5

 

 

 

Значения

коэффициентов п/

 

 

 

Нагрузочный

 

 

 

 

 

 

1

 

Лі

п2

п4

Пъ

«в

«7

 

режим

 

Первый

19

16,4

0,27

63,38

—29,4

47,4

0,76

 

Второй

0,6

0,74

0,27

3,6

—1,4

2,9

0,76

 

Нагрузочный

«8

П9

П іо

Пх1

«12

«13

«14

 

режим

 

Первый

70

46,7

—2790

50

-1 7 ,5

—20,5

—832

Второй

4,9

3

—120

5,3

-1 ,1

—1,4

—35,6

на основе тяговых характеристик, приведенных в гл. I. Для каж­ дого варианта внешних нагрузок определен установившийся кине­ матический режим комбайна, т. е. найдены начальные значения координат х 0, и0, у 0, z0, которые использованы при определении коэффициентов гі/ линейных уравнений.

Выражения для изображений выходных координат х, и, у , 2 получаются в матричной форме из системы (IV. 16) в виде

X

ÎX

 

и

ÎU

(IV. 20)

= А-1

 

УÎ у Н~ ^Ю^вар

 

Z

fz ' ^14^вар

где Л "1— матрица

обратного преобразования для передаточной

матрицы

(IV. 17).

 

Элементы матрицы А ' 1 вычисляются по формулам (IV.9) и (IV.10). Применительно к системе (IV.20) матрица А " 1 в соответ­

ствии с формулой

(IV. 11)

имеет вид

 

 

 

 

Wxx(p)

WuÄ P)

W yx(p)

w zx{p)

 

А-1 =

Wxu(p)

Wuu(p)

W yu(p)

WAP)

(IV.21)

 

Wxy(p)

Wuy(p)

W yy(p)

Wzy{p)

 

 

WXZ(P)

W A P )

w yz{p)

wzz{p)

 

Решения для изображений выходных координат х, и, у, z

могут быть получены также на основе теоремы Крамера

[32 ] без

использования передаточных функций W (р). Однако передаточ­ ные функции WH (р) позволяют:

1) определить амплитудно-фазовые, а также отдельно ампли­ тудно-частотные и фазо-частотные характеристики линейных моде­ лей рабочих органов и частей зерноуборочного комбайна, которые можно использовать для качественной и количественной оценки его работы;

2)найти аналитические выражения переходных процессов для любой из выходных координат при любом внешнем воздействии;

3)рассчитать автоматический регулятор постоянства подачи хлебной массы и провести анализ переходных процессов в разом­ кнутом и замкнутом контуре динамической системы с регулятором;

4)провести точный анализ поведения машины при воздей­ ствии на нее внешних возмущений в виде стационарных случай­ ных функций и оценить при этом статистический характер каждой из выходных координат;

5)определить коэффициенты усиления kSi = Wп (0) для каж­ дой из передаточных функций, приняв оператор р — 0; величины kji показывают в безразмерной форме значение соответствующей

выходной координаты / при действии единичного статического возмущения ft\ при ficm ф 0 х, — Wп (0) ficm\

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