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

книги из ГПНТБ / Коптев, В. В. Вопросы динамики сложных сельскохозяйственных агрегатов

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

1 ^,

(к + 1 )-й массами;

к + 1 — параметры приводимой цепи;

со и

соК— угловые скорости вала приведения и приводимой

 

массы.

Приведение поступательного движущихся масс трактора к принятому при расчете участку приведения осуществлялось из­ вестными методами [1 1 2 ].

Iпос

G г2

(59)

 

где 1 пос— приведенный момент инерции поступательно движу­ щихся масс трактора;

G — вес трактора;

г — радиус ведущих колес трактора;

iT — передаточное отношение трансмиссии; g — ускорение свободного падения.

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

Движение системы, изображенной на рис. 8 , а, может быть описано следующей системой дифференциальных уравнений:

Mi — Ii <?i — Р12 (?2 — ?i)

- С12 (ср2 — ?!) =

0;

 

Ii ? i —

Pi- x,i ( Ti—x — ? i ) —

C i _ u

(<pi - i —

tpil

+-

 

+ Pi,i+i(?i—•?i+ i)+

Ci,i+i (®i —

+x) =

Mj (t);

/gn\

In ?n

Pn—l.n ( f n — l

?n) — Cn—l,n (<pn-l

 

 

Tn) = Mc(t),

Всистеме уравнений (60):

Mi — вращающий момент двигателя;

В — момент инерции .вращающихся частей двига­

теля;

сопротивления на

pu +i — коэффициент неупругого

участке между i-й и (1 +

1 )-й массами;

См-и — жесткость валопровада на том же участке;

Мс (t), Mi (t) — моменты сопротивлений, приложенные к конеч­ ной и промежуточным массам.

79

Для ^перевода системы уравнений (60) в относительные ве­ личины производим заменау переменных с помощью выражения

'1

 

 

(61)

Величина S| показывает относительные отклонения текущих

значений угловой скорости

вала

двигателя

от скорости

иде­

ального холостого хода соо-

 

в системе

уравнений (60)

про­

Преобразование переменных

изводилось следующим образом:

 

 

 

 

d <pf

d Ю|

 

 

 

dt2 ~

dt

 

 

d ?i + i di

По условиям замены переменных

( ю1+1 — “ i) =

( Sj — S^i);

(?i+i — ? i)= .f ( « 1 + 1 +

co,)dt = co0 j (S, — Si ti dt

Произведя замену переменных с помощью выражения (61) и разделив правые и левые части уравнений (60) на соо, получим следующую систему интегро-дифференциальных уравнений, за­ писанных в относительных величинах:

+ Ii S,

+

p12 (S, -

S2) + C12 f ( S, -

S2)dt = 0;

I , - 3

,

-

Si_i)

4- С,_! J (^ — S,_i)dt +

+ J 1,1+1

(S| —■S,+i)

+

C,,1+1 j (Sj — Si + i)dt = тп;(t );

 

 

 

 

 

 

(62)

InSn +

Pn-i,n(Sn —

Sn-i) + Cn_i,n J (Sn -

Sn-i) dt =

 

 

 

 

 

= mc(t),

 

80

где

mt

Mi

— относительный вращающий момент

co0

 

 

 

двигателя;

mc(t) = — ,

т ((г )=

— —------ относительные моменты сопро-

“ о

 

 

“о

 

 

 

тивления, приложенные к i-й

 

 

 

массе и к ведущему звену.

Так как в качестве вала приведения принят коленчатый вал двигателя, все введенные в систему уравнений (62) параметры считаются приведенными к коленчатому валу.

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

тическую характеристтику в виде

 

Mi =А—Всо

(63)

или при оперировании переменной S!

М, = — S,,

К

где А и В — постоянные, 1завиеящие от протекания статической характеристики;

со— текущее значение угловой скорости коленчатого ва­ ла на линейном участке регуляторной ветви;

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

в самом двигателе.

Характерной’ особенностью переходных процессов как в элек­ тродвигателях, так и в двигателях внутреннего сгорания' являет­ ся зависимость падения (или возрастания) угловой скорости двигателя от скорости нарастания (сброса) момента двигателя, а также от величины постоянной времени двигателя. Для двига­ теля внутреннего сгорания эта величина включает в себя период задержки воспламенения, запаздывание регулирующего органа, запаздывание подачи топлива ит. д. Используя электромехани­ ческие аналогии, для двигателя исследуемого типа может быть принято следующее обобщенное уравнение [17]:

81

Si =

kco0(M ,

+■ T n p M j),

(64)

где Si — относительная

угловая

скорость двигателя,

определяе­

мая 1по формуле (4—30)

 

Тпр — постоянная времени двигателя; <±>о— угловая скорость идеального холостого хода.

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

тивления целесообразно задавать в виде функции

 

Мс t — Мс т] (t>,

(65)

где Мс = const, а г) (t ) — функциональный множитель.

Задавая различные значения функциональному

множителю,

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

• 1. При Mcr](t) = 1 функция нагружения изображает так на­ зываемый единичный скачок (рис. 9, а). В качестве нагрузочно­ го режима единичный скачок момента сопротивления представ­ ляет собою мгновенный наброс момента сопротивления, малове­ роятный в реальных случаях нагружения.

Рис, 9. Типы возмущений

82

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

2. Функция нагружения, график протекания которой изобра­ жен на рис. 9, б, представляет замедленный яаброс момента со­ противления. В интервале времени Т S® t ^ 0 значения функцио­ нального множителя исчисляются по формуле

 

rl(t)

=

,

где Т — время нарастания нагрузки.

В

интервале времени t >

Т

функциональный множитель

т] (t) =

1 .

 

 

Функция нагружения указанного типа соответствует случаям включения машины-орудия на ходу трактора, изменению нагру­ зочного уровня трактора, торможению и друпим случаям, на­ гружения.

3. При значении ri (t) = sin cojt функция нагружения пред­ ставляет собою гармонически изменяющийся момент сопротив­ ления заданной частоты (см. рис. 9, в). Использование указан­ ного способа нагружения необходимо для аналитического опре­ деления частотных характеристик исследуемой динамической системы.

5.На рис. 9, г изображен график функции нагружения, для

которой p(t) = 1 + S sin oil t, где S — характеристика амплитуды

колебаний в долях Мс . Функция изображает гармонические ко­ лебания момента сопротивления вокруг заданного нагрузочного уровня.

В соответствии с теорией приближенных вычислений спектр функций указанного вида может изображать изменения момента сопротивления при работе трактора на большинстве сельскохо­ зяйственных операций [1 2 ].

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

Система уравнений (62) с учетом

характеристики двигателя

может быть окончательно записана в следующем виде:

S

=

km0 (m1 -

Тпр mt);

 

mt

-j-

I ( Sj + ^12

(Si — So) -f- C (2 |

(Si ■ So)dt = 0;

83

h S ^ P - M 'S , — S ,- i ) +

I S, - Si-!) d t +

+ 31,1+1 (Si Si_i) + Ci ,i+ i

J ( Sj i Si+, ) d t ~ m i ( t ) ;

......................................................................

(66)

In S n + P n - i , n ( S n — S n --i) + C n - i , n J" ( S n — S n - i ) d t —

=mc (t)

i= 1 , 2 . . . , n.

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

принять нулевыми, т. е. при t = 0 S( = 0 , Sj = 0 .

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

' [ к ш0 (1 + T n p p ; j щ , ( р ) — S j ( р )

0;

mi(p)

+

OlP + Pl2 + -

C12)S!(P) -

(P12

+

 

 

 

P

 

 

 

 

1----- (-4 2 ) S2 (p) — 0 ;

 

 

 

 

 

P

 

 

 

 

 

 

 

(Pi -

1,1

"1-----C i _ i , i ) S i _ i . i ( p )

+ ( I i P + P i , i - H ------ С.,, -j +

 

 

P

 

 

 

P

(67)

+ Pi,1 + 1

4--------C i , i + i ) S j ( p ) — (Pi,i+i 4-

 

 

 

 

 

 

 

P

 

 

 

 

 

H-------- Ci,i + i) S i + i ( p ) —

m i ( p ) ;

 

 

 

p

 

 

 

 

 

 

 

( P n - l , n 4 -----Cn-l,n)Sn-l(p)

(In P

+

Pn-I,n

+

 

 

p

 

 

 

 

 

+ —

C n—i,n) sn( p ) =

m c ( p ) ;

 

 

 

P

 

 

i =

2,

, n.

 

 

 

 

 

 

 

84

Если обозначить дифференциальный оператор при оператор­ ных изображениях переменных величин, стоящий в ш-й строке и k-м столбце, через атК(р), а индексы ш и к — изменяющимися в пределах 0 *£ m < n, 0 *£ к s£n, то детерминант системы, опи­ сываемой уравнениями (67),

о■о —«Ю

U(p) =

0

 

0

 

0

а,л

 

0

0

0

 

«и

«12

0

0

0

 

^21

«22

0

0

0

( 68)

0

0

 

1,п—2 #п—1,п--1 «г

 

0

0

0

,п—1

а,

 

Выражения для передаточных функций системы могут быть получены из детерминанта системы.

W mJ p)

Uo(Р) .

Ws,(p)

Ui(p>

(69)

 

 

 

U(p) ’

 

 

U(p)

 

где Wm,(p), Wsj (р) — передаточные

функции

для

относитель­

 

 

 

ного момента двигателя и относительной

U0 (p),U

 

 

скорости i-й массы;

 

 

1 (р) — соответствующие

им алгебраические до­

 

 

 

полнения.

 

 

 

 

 

Передаточная функция для момента сил упругости на участ­

ке между i-й и

(i +

1) -й массами из [19]

 

 

 

 

 

 

 

 

П

 

 

 

 

 

 

1

-

S

IkPWSl(p)f

W « u + 1 (р)

Uu +i(p>

k=i-rl__________

(70)

 

 

 

Uip) -

1 + Tu+i (Р)

 

 

где Т — 1,1+1- — постоянная времени. ^и-н

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

85

 

соких порядков

возника­

 

ют

некоторые

затрудне­

 

ния,

справиться

с

кото­

 

рыми помогают вычисли­

 

тельные машины.

 

 

Вместе с тем. имеется

 

возможность

 

установить

 

ряд

важных

выводов и

 

положений

при

анализе

 

двухмассовой

 

системы,

 

состоящей из двигателя и

 

приведенных

соответству­

 

ющим образом остальных

Рис. 10. Переход от п-массовой

подвижных

частей

трак­

к двухмассовой

тора

(рис.

10).

Переход

 

от простой крутильной си­

стемы, изображенной на рис. 8, к рассматриваемой .двухмас­ совой системе осуществляется с помощью следующих зависимо­ стей [18]:

пп

i= i

i = i

где 1[ — приведенная податливость двухмассовой системы;

1\ — податливость i-го участка приводимой системы;

1,' и I,- ■— моменты инерций масс двухмассовой системы и

для участка приводимой крутильной системы. Определение искомых переходных процессов, а также значе­

ний изучаемых величин с помощью двухмассовых систем являет­ ся в известной мере приближенным. Погрешность при проведе­ нии таких исследований для рассматриваемого типа приводов составляет 3—5% [19]. Указанная точность вполне приемлема для проведения инженерных расчетов.

Система уравнений, описывающих движение двухмассовой крутильной системы: .

k

o)0 (m

-I- Т пр 1Щ ) — S =

0;

 

m , + I j S t + M S ! — S 2) + С 12 J ( S , — S 2) d t = 0;

( 7 1 )

h

S 2 - f

p12( S 2 — S , ) f S a2

j ( S 2 — S O d t = m c ( t ) .

 

m

Значение детерминанта системы уравнений (71), записанного в изображениях,

к « 0(1

ТпрР)

— 1

О

U(p) =

1

II р ф Р12Ч— "

- ( p,s+ cf )

 

р

(72)

 

0

- ( 3” + т )

I2 Р 4 Pi 2 Ч-----—

 

 

р

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

WM(p )= —

-------- P-± d..iP_+ jo------

.

(73)

)2

Р ^+ Дз Р3+ й 2 р 2+йх р + д 0

 

 

Значения коэффициентов а, и d, в выражении (73) определены экспериментально-теоретически.

Из физической сущности колебательных явлений в приводе, связанной с учетом диссипации энергии в материале, корни ха­ рактеристического полинома р4 -)-а3р3 + агр2 + аф + йо в общем случае можно искать в виде:

P i 2 = — Yi

+ i k i ;

 

Рз4 = —Y2

+ ik2 -

(74)

Для определения переходной функции движущего момента можно пользоваться известной формулой разложения, применяе­ мой при обратном преобразовании Лапласа:

к

Ь м

( t )

=

| W m ( P ) | p = o 4-

Ul2 (Р>

е Рк4

(75)

У

 

 

 

 

 

Р U ( Р )

Р = Рк

 

 

 

 

 

i = l

 

 

Определим

смысл

математических преобразований, связанных

с подстановкой

в формулу (75)

пары комплексно-сопряженных

87

корней (например,

pk =

—Y i+ ik ). Определение первого члена

в выражении (75),

т. е.

|W m ( p ) | p= o, не вызывает затруднений.

Внашем случае это выражение равно единице.

Врезультате подстановки пары комплексно-сопряженных

корней во второе слагаемое имеем

£

-iCPL

ер„; = (Д,

iB,) е~(—

+

1=1

Р u ^

pk = — 7i±ik

 

 

+ (А; — iB i)e-(T i-'ki)t= 2 А ,е -н ‘

+ (А, + 1В ,)е -‘ -*‘ +

+ (А, — iB ,)e 'V .

Принимая во внимание формулы Эйлера [98] для перехода

от показательных функций к тригонометрическим:

 

. е1 1 ‘ =

cos k1 1 + i sin kit;

 

e - ’

i f =

cos кj t —

i sin ki t,

имеем

 

 

 

 

 

 

(A, +

B,) e lkif -Ь (А,-—

iBi

t =

2

(AiCOS kj t -Ь B| sin kjt).

Вводя обозначения

Aj

— Csin9j,

Bi

= Ceos 9; и учитывая,

что С =

|/~ А? + В?2

, окончательно запишем

 

(A, cos k| t-f В! sin ki t)

= -j/"Af

+

Bf cos (kj t-f- 6 ,),

где tg 6 , = A .

Aj

 

 

кор­

Таким образом, в общем случае комплексно-сопряженных

ней переходная функция для движущего момента трактора

опи­

сывается выражением

 

 

 

 

2

 

hM (t) =

1 +

L,e-Ti‘cos(k,t + 0,),

(76)

где

L,- — Aj j/~Af -f- Bf.

 

88

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