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

книги из ГПНТБ / Хачатрян Х.А. Стабильность работы почвообрабатывающих агрегатов

.pdf
Скачиваний:
2
Добавлен:
24.10.2023
Размер:
9.84 Mб
Скачать

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

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

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

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

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

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

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

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

80

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

14. АНАЛИТИЧЕСКОЕ ОПРЕДЕЛЕНИЕ УСТОЙЧИВОСТИ ПРОЦЕССА ПРИ ПОСТОЯННО Д Е Й С Т В У Ю Щ И Х ВОЗМУЩЕНИЯХ

Для различных процессов допустимое значение дисперсии и тем самым границы допустимого отклонения процесса нужно выбирать исходя из требований, предъявляемых к данной тех­ нологической операции. Естественно, эти требования должны быть обусловлены также конкретными условиями.

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

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

Динамические свойства агрегата характеризуются преобразо­ ваниями, происходящими в нем над входными воздействиями, в результате чего получаются выходные переменные. Наиболее простая одномерная динамическая система воздействию X (t) отвечает реакцией Y (t). Преобразования в системе могут быть любого вида (умножение, дифференцирование, интегрирование и т. д.). Воздействие же X (/) может быть как детерминированным, так и случайной функцией времени. Свойства системы, преобра­ зующие X {t) в Y {t), выражаются оператором данной одномерной динамической системы.

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

6 X . А. Хачатрян

81

являться линейной, т. е. описываться линейными дифферен­ циальными уравнениями.

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

Одномерная динамическая система может описываться линей­

ным дифференциальным

уравнением

 

 

 

 

 

^

+

a i

^

+

. . . + a

Y

{ t ) =

x m .

(57)

В общем

случае

 

 

 

 

 

 

 

 

 

 

 

v , , ,

,

dmZ(t)

 

,

, dni~lZ(t)

 

.

.

,

 

х

(0 = fro — ~

L

+

ьу

\

+

• • • + bmz

(t),

где Z (t) — воздействие, в общем случае являющееся функцией времени.

Решение уравнения (57) имеет вид:

пt

 

у (о = Ух (О+у2 (о = 2 с т

(о +

1 в е . ух

<58)

 

; = 1

о

 

 

где Y\

(t) — общее решение однородного

уравнения;

Y2 (t) —

частное

решение неоднородного уравнения;

г/( (t) •— система не­

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

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

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

[10].

 

 

 

 

При постоянных коэффициентах дифференциального уравне­

ния

(57)

 

 

 

 

 

 

g(t, t i )

= g ( t -

ti)

 

частные интегралы однородного уравнения получают

вид eV,

где

X/—корни

характеристического

уравнения.

 

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

производства большой

интерес

представляют процессы, носящие установившийся характер. Если

X (t) — стационарная

случайная функция, а система асимпто­

тически устойчива,

то

установившийся „.процесс

описывается

решением

 

 

 

 

 

t

 

y*{t)

= \g(t-h)X(h)dtl,

(59)

82

Стационарным случайным воздействиям устойчивая линейная система отвечает стационарным установившимся процессом. Его дисперсия определяется из выражения

t

t

 

 

 

De=]

IgV-tJgit-12)

Kx (ti-h)

dt,dt2

(60)

или

D„=\g(tJ

\g{k)Kx{t1-t2)dt2 dtlt

(61)

где Kx (tx — /2 ) — корреляционная

функция стационарной

слу­

чайной функции X (t).

 

 

Вывод и обоснование приведенных формул имеются в спе­

циальной литературе [5, 10, 11].

 

 

Определим дисперсию выходной

переменной одномерной дина­

мической системы, описываемой часто встречающимся уравнением

^

L + o l ^ ) . + e e y(0 = x W ,

(62)

где аг и а0 — постоянные коэффициенты.

 

 

Чтобы система, описываемая уравнением (62),.была устойчива,

необходимо аг > 0.

уравнения,

описыва­

Если известен

вид дифференциального

ющего поведение

системы, то при помощи

формул (60) или (61),

используя корреляционную функцию стационарного случайного

воздействия, можно определить дисперсию случайного

процесса

на

выходе

системы. Однако

вычислить Dy указанным

способом

довольно

сложно,

поэтому

иногда прибегают к

другим спо­

собам [5].

 

 

 

ах и а0

 

 

В зависимости от соотношения коэффициентов

урав­

нение (62) имеет три решения.

 

 

 

 

При а х

< 2 Yao

выражение (62) называется уравнением

колеба­

тельного звена

второго порядка. Весовая функция в этом случае

 

 

 

 

,

fli (t-ti)

 

 

R_.

 

 

g(t-t1)=±e

2 . s i n n ^ - ^ ) ,

 

 

( Ь 3 )

 

 

 

 

"1

 

 

 

 

 

2

 

#1

 

 

 

 

 

где

iii — a0

—.

 

 

 

 

 

Когда ax > 2 У~а~0, выражение (62) называется уравнением неколебательного звена второго порядка. Весовая функция при этом

 

g (t - tx) = е ~ "Г"( ' ~ h ) s k n2 { t - t j ,

(64)

2

а \

 

где п2 =

-~ — а0.

 

6*

83

Если

-j-

=

а0

(консервативное звено),

 

 

 

 

 

§

У - Ц

= е - ^ и

- ' 1 )

(t-h).

 

(65)

Допустим,

на

вход

системы

(62)

при

условии

at < 2 ] / а 0

подается

стационарная

случайная

функция

X (t), корреляцион­

ная функция

которой

 

 

 

 

 

 

.

/Сл- (т) = al-e ~ а 1 т 1 (cos

Рт +

-у sin Р | т |) .

(66)

Используя зависимости (60) и (63), можно записать

—00 —со

 

 

X sin пх

(t — tx) e

2

"

sin nx

(t — /2)

dtx

dt2.

 

 

 

После

подстановки

значений

и

интегрирования

получается

 

 

aal (сг +

р ? )

2 -

2/г2 )2 +

4 [а2 fa

-

2/i2) -

2

2

+

/г)]

L

 

«

a0h

 

[(а, — 2/г2 )2

+ 4 (а0 а2

— 2а2 Л2 +

2 )]2 +

а4

1

 

 

,

г ^ + 4 [ / г

г ( а Ч р 2 ) - а 2 К + а 2 ) ]

 

 

 

 

 

 

 

 

И + 4 ( Р 2 * 2 - Р 2 « 2 - * V ) j 2

+ «4 '

 

 

 

 

где

/г =

а2

== а0 а 2

— 02;

а3 = р2 — а 2

— а0 ;

а4

=

16а2 р2 х

Х ( а 3 + 2/г2).

Решение аналогичной задачи приводит А. А. Свешников [11 ] . Аналитическое выражение дисперсии выходной величины при подаче на вход системы случайной функции с корреляционной функцией (66) имеет сложный вид. Однако при других видах корреляционной функции случайной функции X (t), как увидим дальше, выражения для Dy получаются сравнительно простыми. Зависимость (67) и аналогичные ей другие уравнения выра­ жают связь между дисперсией выходной переменной и пара­ метрами системы при определенном виде возмущений. Имея эти выражения, часто удается параметры машины подобрать так, чтобы

получить желаемое значение дисперсии на выходе.

Дисперсию выходного переменного одномерной динамической системы иногда удается определить проще. Обозначим, как и раньше, выходную переменную Y (t). Пусть иа одномерную динамическую систему приложено только одно воздействие X (t).

84

Систему описывает дифференциальное уравнение, которое в сим­ волической записи имеет вид

Dn(p)Y(t)

= Qm(p)X(t),

(68)

где Dn и Qm — полиномы

степени

/г и т.

стационарности, то

Если решение (68) обладает

свойством

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

S . H = ^ S | S , H ,

(69)

где Sx (со) — спектральная

плотность входного

возмущения.

Отношение

 

 

 

Qm (tto)

Ф (tto)

(70)

Dn

(ico)

 

 

называется частотной характеристикой линейной динамической системы и представляет собой отношение установившегося про­ цесса к гармоническому воздействию. В дальнейшем в качестве воздействия рассматриваются стационарные случайные функции X (t). Последние могут быть выражены интегральным канони­ ческим представлением [10]

со

X{t) = m.v + j V(w)t!aida,

(71)

— со

 

где тх — математическое ожидание; V (со) — «белый шум» пере­ менной со, интенсивность которого равна спектральной плот­ ности Sx (со) случайной функции X (t).

Координатными функциями (71) являются показательные функ­

ции

X (t, со) = е'"и<.

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

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

Если в формулах (60) и (61) дисперсия выходной переменной связана с корреляционной функцией входного возмущения слож­ ным интегральным выражением, то в зависимости (70) их спек­ тральные плотности связаны частотной характеристикой системы. Последняя имеет простой вид и может быть выражена параметрами системы. Это позволяет широко применять выражение (70) для определения дисперсии выходной переменной, а также при реше­ нии других практических задач.

85

Пользуясь выражением (70), можно вычислить частотную характеристику линейных систем. Например, для системы опи­ сываемой уравнением (62), частотная характеристика

 

 

Ф (ко)

=

 

1

 

 

 

 

 

а 0

+ aja

со2

 

 

 

Если на вход системы подается возмущение со спектральной

плотностью

 

 

 

 

 

 

 

 

 

 

 

 

 

5,(со)

=

Dx

 

 

 

 

 

то

дисперсия

выходной переменной

[10

 

 

 

 

D

_

a0at

А с

(и +

ar)

 

0 )

(72)

 

 

 

»

( a 2 + a x a + a

 

где

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

корреляционной

связи.

 

 

Когда этому же звену в качестве входа подается случайная

функция, обладающая

свойствами «белого

шума», т. е.

 

 

 

S ,

(со)

 

 

const,

 

 

то

дисперсия

выходной

переменной

на

некотором

диапазоне

частот [10]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(73)

 

Покажем

на примере,

 

как

можно

использовать

дисперсию

выходной переменной для выбора параметров динамической си­ стемы.

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

действием

случайного

возмущения

она отклоняется от

направления

движения

под углом

9 (рис. 17).

В центре тяжести машины дей­ ствует сила сопротивления F,

лв каждый момент времени направленная обратно абсолютной

 

 

 

скорости движения

этого

центра.

 

 

 

Пусть

центр

 

тяжести

 

машины

 

 

 

находится

в

точке

О.

Кроме

F

 

 

 

на

дисковые

батареи

действуют

 

 

 

осевые

силы

Rlt

R2, R3,

R4,

ко­

 

 

 

торые

помимо

глубины

обработки

п

-

.

почвы,

ее

свойств,

главным обра-

Рис.

17. Схема сил,

действующих

 

'

 

 

от значения

 

*

на дисковую почвообрабатывающую

3 0 м

зависят

углов

 

машину

 

атаки

дисков.

Но

после отклоне-

86

ния машины на угол 8 углы атаки дисков также изменяются. Поэтому можно принять

 

 

 

/?i

=

Ri = R sin (ф — 9);

 

 

 

 

 

#s

=

Яз = R sin Сф + 9),

 

 

где R — начальное значение осевой силы при

9 =

0; \\> — угол

атаки дисков.

 

 

 

 

Момент

сил RI, R2> R3< Rq относительно

точки

прицепа А

 

 

 

^MR

=

2R (2LC + d) cos2 tJ) sin

9,

 

где

L c

расстояние

центра тяжести машины до точки прицепа,

d =

00

v

 

 

 

 

 

Момент силы F относительно точки А, при небольших откло­ нениях машины,, когда можно принять sin 0 9, cos 9 1, на основании приведенного в § 12 примера,

 

 

 

 

FI? .

 

 

 

 

MF

= FLCB - j -

9,

 

где

vu — скорость

машины.

 

 

 

 

Следовательно, дифференциальное уравнение отклонения ма­

шины от направления движения

 

 

где

 

0 +

сф + ао 0 = X (0,

(74)

F

..

 

2R(2Lc-Ld)cos^

+

F L c

 

п _

 

 

 

а ° -

 

щ

 

где

т — масса машины.

 

 

 

 

 

Таким образом,

получили

дифференциальное

уравнение вида

(62). Если при работе на машину действует возмущение, спек­ тральная плотность которого постоянна и имеет интенсивность с на определенном интервале частот («белый шум»), то согласно зависимости (73) для дисперсии отклонения центра тяжести

дисковой почвообрабатывающей

машины

 

У ~ F [2R (2LC +

d) cos21|> + FLC] '

(75)

Используя выражение (75), необходимо проверить выполняемость условия

a1<2\/raQ.

Из формулы (75) видно, что высокая интенсивность спектраль­ ной плотности входного возмущения и большая скорость движе­ ния увеличивают дисперсию Dy прямо пропорционально. При больших значениях угла атаки дисков равномерность хода ма­ шины уменьшается. Для выяснения влияния расстояния L c на Dy можно принять d = 0, тогда после сокращения в числителе

87

остается L c в степени единица. Это вполне объяснимо, так как при большом значении Ьс небольшим угловым отклонениям

машины соответствуют довольно заметные

перемещения точки

О

в поперечном направлении. Любопытно,

что с увеличением

d

устойчивость движения машины возрастает.

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

5 = D x а

хя а2 + со2

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

реляционной функцией вида (66), для определения Dy прихо­ дится использовать выражение (67). Такой анализ должен про­ изводиться для конкретных значений параметров, взятых при исследовании модельных образцов машин.

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

В § 12 рассматривалось влияние показателей переходных процессов на качество выполняемой технологической операции. Для оценки интенсивности затухания переходного процесса было использовано отношение

F - °

Для рассматриваемой дисковой почвообрабатывающей ма­ шины

 

Е

= - S - [2R (2LC +

d) cos2 ч|> +

F L J .

 

 

F Lc

 

 

 

Так

как при уменьшении Еот

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

затухания пере­

ходного

процесса

повышается,

то

время переходного периода

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

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

Изменения параметров системы могут влиять по-разному. Известно» что после переходного периода устойчивая система работает лишь под действием возмущения X (I). Следовательно, при выборе параметров машины необходимо останавливаться на результатах анализа установившегося режима работы. Но в сельском хозяй­ стве есть процессы, где такой подход может привести к ошибоч­ ным выводам.

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

жиме

под воздействием стационарного случайного возмущения

X{t).

При встрече с большим камнем, сопротивление которого

нельзя отнести к категории обычных возмущений, рыхлитель сильно отклонится от равновесного состояния. Начнется переход­ ный процесс, в котором большое значение будут иметь собственные колебания системы. Может случиться так, что при работе рых­ лителя установившийся режим будет незначительным. Остальное время будет приходиться на переходные процессы, чередующиеся случайным образом и иногда накладывающиеся один на другой. Фактически процесс превратится в нестационарный. Ввиду слож­ ности исследования такого процесса определить аналитически параметры системы невозможно. По-видимому, параметры ма­ шин, предназначенных для работы в таких условиях, можно выбрать исходя из условия интенсивного затухания переходных процессов. Поэтому при выборе параметров рыхлителей, рабо­ тающих на сильно каменистых почвах и в других аналогичных условиях, следует пользоваться выражением типа (51).

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

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

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

условиях.

 

 

 

 

 

 

 

Так, если динамическая система описывается линейным урав­

нением

 

 

 

 

 

 

 

 

 

 

 

 

^g!L

+

a0Y(t)

=

0,

( f l o > 0 ) ,

(76)

то

его

решение

имеет

вид

 

 

 

 

 

 

 

 

Y (t)

=

сг

cos

a0t

+

с 2 sin a0t.

 

 

Это решение

можно представить в форме

 

 

 

 

Y{t)

 

=

A sin (a0t

+ ср0),

 

где

А

и ср0 — постоянные.

 

 

 

 

 

89

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