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

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

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

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

Нетрудно заметить, что, если корни характеристического по­ линома отличны от корней (74), передаточная функция опреде­ ляется подстановкой соответствующих их значений в выраже­ ние (76). В результате подстановки параметров исследуемой системы (методика определения последних представлена в гла­ ве IV) и вычисления значений коэффициентов получим выраже­ ние для передаточной функции ведущего момента на III пере­ даче трактора:

wu* (р) = _________Р2- ± 1-0-.102Р + 1 ' 3 2 : 103_________(77)

р4+ 1,1 • 102 р3+2,20-103 р®+1,54Ю4 р+6,30Ю4 * v '

Корни характеристического уравнения равны:

pi = —84,4; р2 = —16,6; р3,4= —3,0—Вь

Выражение переходной функции на III передаче трактораг

h”r(t)

= 1 +

0,017e-16’6t — 0,145-10 -4 e - 84’4t +

 

 

+

0,84е-з‘ cos (6 t +

1,49).

(78)

График функции приведен на рис. 11, а.

передаче трактора

Соответствующие

выражения

на

VII

имеют:

 

 

 

 

 

 

передаточная функция —

 

 

 

 

 

 

р 2 +

1,0-10* р +

1,32-Юз

 

W m ‘ ( p ) =

131,0 р4+1,02-102 р34-1,78■ 103 р 2+2,20-104 р + 1,76 -105

 

 

 

 

 

 

(79)

корни характеристического уравнения —

 

 

.р, = —84,2;

р2 = 12,6;

р3, 4 =

—2,89 ± 12,9 i;

 

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

h™(t) = 1- 0,89• 10~e~84’2t —0,123e-12’6t +

+ 0,92е-2-8 9 cos (12,9t + 1,52).

(80)

89

 

График

функции

приве­

 

ден на рис.

1 1 , 6 .

 

 

 

 

Для суждения о реакции

 

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

системы

на

 

возмущающее

 

воздейст­

 

вие

гармонического

 

типа

 

М с (t)

= M

c sincot

восполь­

 

зуемся

выражением

частот­

 

ных

характеристик.

Такие

 

характеристики

могут

быть

 

получены

непосредственно

 

из

передаточной

функции

 

(77)

путем

замены

в

них

 

оператора

дифференцирова­

 

ния

р мнимым

переменным

Рис. 11. Графики переходных функ­

ico (со — частота

гармониче­

ций на III и VII передачах

ской составляющей момента

 

сопротивления).

 

 

 

 

В результате такой подстановки получено комплексное чис­

ло, в общем виде описываемое ‘выражением

 

 

 

 

 

 

W m(!m) = и(ш) + iv( <в),

 

 

 

 

 

 

или, переходя к показательнымфункциям,

 

 

 

 

 

 

 

U* (to) +

V* (СВ)

е 1arc(g —

 

 

 

 

Величина ]/^u2 (to)-ф v2 (to) — есть модуль комплексного чис­

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

А(со) =]/~u2(co) +v2(co)

(82)

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

, ч

С 12 m6o>6 - fm4 ( й * + т 2 < о 2 + ш 0 .

/0 0 \

U ((В) =

------------------------------------ В “ 8 + 1 в

,

(8 8 )

 

“ ^-Нз “ 2+Iq

 

90

V(co) =

C12

ГП5 « 5+ П 1з «З+ПЦ «

(84)

Ij

«*+lg «5+ l4«4+lj <*>2 1q

 

 

В результате подстановки значений параметров системы и вычисления коэффициентов выражения для элементов амплитуд­

но-частотных

 

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

момента трактора

име­

ют вид:

 

 

 

 

 

 

на III передаче трактора —

 

 

 

и > ) =

 

 

« г>+7,3- Юз «<+9,24-105

8,4-1 Щ

(85)

—«*—7,3-103 «5—1,7- К 6 « 4+4,5•107 -со2—4,0-109

 

 

V(co)

=

_________ 2,6со5+2,09-10* соз+ l , b io ? «__________ _

 

-«3 -6,9 -10» ms+ 1,0- 10в «4+ 1 ,4-10* о)2—3,1 • 1010 ’

 

 

 

 

 

 

 

 

А((в) = 11,4 У U * (u > )+

v2 (u);

 

на VII передаче трактора —

 

 

 

,

Л

=

«5+ 7,2 -103 «4+3,1-105

0)2-2,3-108

 

и(со)

---------------- :----------■----------------------------------- ;

 

 

 

 

—«8-6,9-103 «6+1Д). 106 0)4+1,4-108 « 2-3,1 •1010

 

/ ) _

___________ 2,6«5+2,09-lQ4«3 - 1 ,1 -Ю7« ___________

( 86)

^~ — «8-6,9-103 «5-1,0-106 «4+1,4-10* со»2— 3,1 •1010

 

Кривые

амплитудно-

 

частотных

 

характеристик

 

представлены на рис.

1 2 .

 

Вычислительные

работы

 

по определению переходных

 

процессов

частотных харак­

 

теристик

производились с

 

помощью

цифровой

элект­

 

ронно-вычислительной ма­

 

шины <«Промины1М».

 

 

Пользуясь

выражениями

 

для переходной

функции,

 

можно определить значения

 

момента

на

ведущей

оси

 

трактора

Mi2 (t) при

 

внеш­

 

нем нагружении, заданном в

Рис. 12. Амплитудно-частотные ха­

виде функции (65)

при раз­

личных вариантах

скорости

рактеристики на III и VII переда­

чах

нарастания нагрузки. В этом

91

случае удобно использование формулы интеграла Дюамеля (теоремы свертки) [126]:

M1 2 (t) = Мс (+ 0) hM(t) + I'

hM(t — т) d т,

(87)

J

UX

 

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

чае корней характеристического уравнения для интервала вре­ мени Т > t ]> 0 имеем:

 

2

 

 

 

 

 

M12(t) = —

+ ^ V

—-----(Tl cos 6

+

К, sin 0 )—

Т

I , T J

i f + k f Ul

1 1

 

1 = 1

 

 

 

 

 

e~Ti t [^. cos(k|t +

0 )

+ kj sin (M

+

0 )];

(8 8 )

В выражении

(8 8 )

А, и В/ — вещественные

части комплексных

чисел;L;= ]/A f

+

В? — модуль комплексного

числа; Ti и к;— ко­

эффициенты при вещественной и мнимой частях комплексно-со­ пряженных корней.

Для упрощения выражения [tj cos (kj t +0)4- kj sin (k} t

+

6 )]

введем обозначения:

 

 

 

 

 

Yi — (osincp;

к,- =

шсо-ср.

 

 

Имея в виду соотношения

 

 

 

 

 

ш2 — yi + k'f; tg ср

=

JL ,

 

 

запишем

 

 

k|

 

 

 

 

 

 

 

[у,cos(kj t + 0 ) + kj sin (k, t +

0 )] =

(yf

+ kf)sin (kjt +

0 +

<p)

После соответствующих преобразований для интервала вре­ мени Т+> t > О

2

м 12 (f) = - f

+ 2- % У

R, + - ^ 4 ^ - sin (kit + а),

(89)

Т

I2 T

лг

 

 

i=l

У и

 

•92

где

* 1 = -

'Т + !<',В| :

г, - k f + T f ; «

= ? + ».

 

Т12+ к=

 

 

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

Т искомая функция

 

 

2

 

 

 

Mj2 (t) “

1 + 2

e_ V (А| sin k,t +

В,'cos kjt).

i = l

Для определения значений А,' и Bj воспользуемся очевид­ ными соотношениями:

|MJ2 (t)[t=T = 1 + 2

e~Tlт (-^i sin ki Т + В' cos к, Т) =

i=l

 

= I M l 2 (*) |t=Ti

2

M = 2 - Ti e~ Ti T(Ai s}n k>т + Bi cos ki? ) + i=i

kie~fiT (AJ cos k,-T — B'sink,-T)

dM j (t)

dt t-T

характеризующими краевые условия для функции в разных ин­ тервалах времени Т > t 0 и t > Т. Решение указанной системы уравнений дает значения AJ и В’, для данного интервала .време­

ни. Знание указанных величин позволяет представить искомую функцию в ‘виде

2

 

M;2 (t) = 1 + J ] R, е-т. tSin (k,t + а,).

(90)

1=1

Подставляя числовые значения входящих в соответствующие выражения величин, имеем для интервала времени T > t > 0 и УП передачи трактора (время нарастания нагрузки Т=0,25се/с)

93

M1 2 (t) = —

+ 0,C41 •10- 4 (e"64’2t — 1 ) + 0,036(e-12'6t

-1 )4 -

0,25

 

 

+

0,306 + 0,310 e-2-89f sin(]2,93t — 1,39).

(91)

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

M1 2 (t) = — + 0,306 + 0,310e“ 2’89t sin (12,93t — 1,39). (92}

0.25

Для интервала времени t > Т значения ведущего момента соответствуют выражению

М12 (t) 1 0,6le-2>89‘ sin(12,93t + 1,07). (93)

Графики переходных процессов для рассматриваемого случая нагружения на VII передаче трактора и темпах нарастания на­ грузки 0,25 и 0,5 сек представлены на рис. 13.

Работа трактора в неустановившемся нагрузочном режиме сопровождается постоянной подачей на крюк трактора возму­

щающих воздействий того

или иного типа. Для выяснения их

 

 

влияния

на

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

 

 

качества

колесной

машины'

 

 

необходимо произвести

ана­

 

 

лиз характера

и структуры

 

 

указанных возмущений.

 

 

Исходными данными для

 

 

такого

анализа

являются

 

 

сведения об изменении тя­

 

 

гового усилия

на различных

 

 

сельскохозяйственных

опе­

 

 

рациях,

полученные в

про­

 

 

цессе

динамометрирования

 

 

трактора. Проведенные из­

 

 

мерения показали,

что

коле­

 

 

бания

крюковой

нагрузки

 

 

носят незакономерный, слу­

 

 

чайный

характер

и пред­

Рис. 13. Переходные процессы

ставляют

собою стационар­

на VII передаче при темпах

нара­

ный

эртодический

случай­

стания нагрузки 0,25 и 0,5

сек

ный процесс.

Поэтому

для

94

характеристики источника возмущений следует применить методы теории случайных функций (87]. Случайная функция мо­ жет принимать различные конкретные значения. Всякая функ­ ция, описывающая исследуемый случайный процесс в резуль­ тате опыта, называется реализацией случайной функции (или возможным значением случайной функции).

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

Математическое ожидание случайной функции представляет собою среднее значение функции на заданном интервале

Т

т “

=

Т2

f x ( t ) d t ’

 

 

- Т

 

где 2Т-— интервал наблюдения при движении машины на уча­ стке.

Корреляционной функцией случайного процесса X(t) назы­ вается математическое ожидание произведения ординат случай­ ной величины для двух моментов времени К и t2 .

При практических расчетах .корреляционная функция рассчи­ тывается по формуле

 

 

 

N

 

 

 

R (t) = R ft* Д)

V х(чА)х

Ь + р)Л],

(95)

 

 

 

v=N

 

 

где

Л— длина

зачетного

интервала

разбиения

функции

 

(шаг);

 

 

 

 

 

N — число интервалов;

 

 

 

 

д

v 1 , 2 , . . . ; р — 4д-; р = 0 ,1 , ..;

 

t

и х — дискретные значения времени,

кратные Д.

 

Формула (95)

является основной для существующих методов

определения корреляционных функций.

 

 

Важнейшей характеристикой стационарной случайной функ­ ции, вытекающей из ее корреляционной функции, является спект­ ральная плотность

00

 

S((o) = |R(t) cos со t ch,

(96)

О

 

где со — текущая частота процесса.

95

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

С целью получения указанных статистических характеристик исследовались реализации случайного процесса при динамометрировании трактора МТЗ-бЛС на пахоте и культивации.

Обработка полученных данных производилась с помощью электронно-вычислительной ‘машины «Минск-14». Анализ резуль­ татов позволяет отметить следующее:

1.Корреляционная функция, характеризующая процесс соз­

дания

тягового усилия на пахоте, затухает за 0,5—0,7 сек

(рис.

14, а).

Рис. 14. Статистические характеристики тягового усилия на пахоте

икультивации2

2.Спектральная плотность, характеризующая спектр частот случайного процесса, имеет максимум дисперсии в зоне низко­ частотных составляющих порядка 0 — 6 рад/сек. Остальные ча-

96

стоты (от 7 до 50 рад!сек) представлены в спектре частот при­ близительно в равной мере (рис. 14, б).

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

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

170—220 кГ, для пахоты — 130— 140 кГ.

2.Спектральная плотность тягового усилия трактора на куль­ тивации на всех реализациях имеет ярко выраженный экстре­

мум. Максимум дисперсии соответствует частотам от 2,6 до 5 рад/сек. Указанные пределы частот являются преобладающи­ ми в спектре исследуемого случайного процесса (рис. 14, г).

3. Положение экстремумов кривых спектральной плотности, полученных для различных скоростей движения, показывает, что одновременно с возрастанием скорости движения агрегата часто­ ты, преобладающие в спектре, смещаются в сторону высокоча­ стотных составляющих. При работе на III передаче преобладаю­ щей является частота 2,6 рад]сек, на VI — 3,2, на VII — 4,1 рад!сек.

4. С повышением скорости движения агрегата возрастают показатели неравномерности тягового усилия. Среднеквадрати­ ческое отклонение значений крюковой нагрузки на III передаче составило 100 кГ, на VI — 186, на VII — 230 кГ.

5. При движении трактора на различных сельскохозяйствен­ ных операциях с различными скоростями движения корреляци­ онная функция затухает за 3,2—3,3 сек (рис. 14, г)..

Как известно [99], между спектральной плотностью входного сигнала Sxx (со) системы, имеющей частотную характеристику k(ico), и спектральной плотностью на выходе ее Syy (со) суще­

ствует зависимость

 

Syy(®) = fk M p S xx(o.),

(97)

где |k(io))( — модуль частотной характеристики.

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

Графики спектральных плотностей ведущего момента при., работе на пахоте и культивации представлены на рис. 15.

4. Зак. 64

97

 

Из

графиков видно, что прй

 

работе на пахоте и VII пере­

 

даче трактора в спектре частот

 

переменного

ведущего момента

 

преобладают частоты,

близкие

 

к 12 рад/'сек (кривая

1). При

 

работе^-на III передаче макси­

 

мальная энергия процесса со­

 

ответствует низкочастотным со­

 

ставляющим порядка

2 радfсек

 

(кривая 2 ).

на

культивации и

 

Работа

 

VII

передаче

трактора вызы­

 

вает колебания ведущего мо­

 

мента

с преобладающей часто­

 

той 4 рад/сек (кривая 3 ).Дан­

Рис. 15. Спектральные плотности ве­

ной

частоте

соответствует мак­

симум

дисперсии случайного

дущего момента на пахоте и куль­

тивации

процесса.

 

 

 

§ 3. Анализ переходных процессов в динамике колесного трактора при неустановившейся нагрузке

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

характера, принимались:

 

1.

Максимальные относительные отклонения ведущего мо­

мента трактора от установившихся значений

 

 

°м — ' м-М;"'Т 1М12макс — М12(со)].

(98)

 

М 12(оо)

 

2.

Коэффициент динамичности, характеризующий степень на­

пряженности привода

и ведущего ходового аппарата трактора,

 

 

М 1 2 макс

(99)

 

 

М 12(оо) ’

 

 

 

связанный с ам очевидной зависимостью р = 1 +

а. В выраже­

ниях

(98) и (99)

— установившееся после завершения пе­

98

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