книги из ГПНТБ / Коптев, В. В. Вопросы динамики сложных сельскохозяйственных агрегатов
.pdf1 ^, |
(к + 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