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

книги из ГПНТБ / Барский И.Б. Динамика трактора

.pdf
Скачиваний:
21
Добавлен:
25.10.2023
Размер:
14.66 Mб
Скачать

где

 

дМ„

,

дМд

Л,

dio,

. А,

dh

 

А3 =

' д

 

 

dp

 

 

 

 

Уравнение турбокомпрессора. Работа турбокомпрессора в со­ стоянии устойчивого равновесия характеризуется равенством

Мт = Мк,

где Л4Т и Мк — моменты турбины и компрессора соответственно.. При введении малых возмущений на основании принципа

Даламбера можно записать

 

 

 

 

 

АМХ

К

= J4 ~~~~ >

(14)

 

 

 

 

 

 

 

 

at

 

где

/ 4

— момент инерции ротора

турбокомпрессора;

 

« 4

— угловая

скорость

ротора турбокомпрессора.

 

Найдем приращения ДМТ и АМК.

 

 

Из теории турбин известны следующие соотношения: момент,

развиваемый

турбиной,

 

 

 

 

 

 

 

 

 

Мт

= — QH т т ) т к ,

 

где

Q = [iF

1

Re Tl

у 2gHT

— расход газов через

турбину (для

 

 

 

0 P l

 

 

 

 

 

 

 

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

 

стационарного потока газов);

здесь

пропускной способности

турбины;

 

 

F — площадь выходного сечения турбины;

 

 

Pi — давление за турбиной;

 

 

 

Ti — температура газа за турбиной;

 

 

Я т — полный напор газов перед турбиной;

 

 

т)т к — полный к. п. д. турбины.

 

 

Коэффициент

напора

 

8gH?

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D 4

 

где

D — наружный диаметр турбинного колеса.

 

 

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

лей, можно без больших погрешностей принять

 

Г, =

(0,94-=-0,97)Г; р , =

1,03р0,

где Т — температура

газа перед турбиной;

Ро — давление окружающей среды.

 

Приняв такое допущение, выражения для момента М т

хода Q преобразуем следующим образом:

 

 

D2

О 2

Q£(D4;

М т = —

<2Ят]тк4

= —

ирас­

(15)

30

цУн = 1,84-

1(Г4

^ — Q T'

— = v,

(16)

 

 

p„DF

 

d)4

 

где

 

 

 

 

 

£ =

#r)T K ;

v =

iiVH.

 

 

Таким образом, для определения М т необходимо знать тем­ пературу газов на входе в турбину и расход газов через нее, час­ тоту вращения ротора турбокомпрессора и характеристику тур­ бины. Последняя, как правило, задается в виде эксперименталь­ ных зависимостей х\тк = г)(#) и р = ц ( # ) . Если допустить, что характеристики турбины в переходном и установившемся режи­ мах одинаковы, то т)(Я) и р ( # ) можно найти по результатам стендовых испытаний.

Из выражения (15)

MT = M(Q, I , со4).

Следовательно, приращение момента турбины можно опре­ делить из уравнения

ДЛ1Т = _™L_ AQ +

А£+™1-

А щ .

(17)

dQ

dt,

d(o4

 

Найдем значения частных производных, входящих в уравне­ ние (17):

дМг

D2

r

Qo

М-п

dQ

8g

Ь 0 « 4 -

Qo

~ Qo

аналогично

дМг

Мы

дМт

Мы

 

 

 

^40

Для определения приращений AQ и Д£ (от независимых па­ раметров системы) воспользуемся следующими функциональны­ ми зависимостями.

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

12600 v

где Vh — литраж двигателя.

В свою очередь, т)у зависит от угловой скорости ©i, плотности воздуха р и от количества поданного в цилиндр топлива, т. е. от

31

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

Q =

Q(©„ h, р).

 

Далее будет показано, что можно принять

 

Q =

Q(co„ р),

(18)

откуда

 

 

 

A Q = ^

X l

+ JS_u.

(19)

да>{

dp

 

Значение частных производных, входящих в уравнение (19), можно найти по результатам стендовых испытаний двигателя.

Приращение А£ (исходя из характеристики турбины) может быть выражено через Av:

 

=

A - Av.

 

(20)

 

 

av

 

 

Учитывая,

что v является

функцией

трех переменных:

v =

== v(Q, Т, 0 ) 4 ) ,

получим

 

 

 

 

Av = — - A Q + — Л Г + — х 4 .

21

 

3Q

дТ

dwt

 

Частные производные, входящие в выражение (21), опреде­ ляются следующим образом:

д\ _

у0

dv

_ v0

dv _

v0

dQ

Q0

дТ

TQ

дш4

w4 0

Температура T выпускных газов двигателя ношения количества воздуха и топлива в камере же от характера протекания рабочего процесса:

 

Г=Г((о,, h,

р).

Поэтому

 

 

 

a t

дТ

, дТ

дТ

А

Г = - —

х, + —-У

+ — и.

 

дш,

dh

dp

зависит от соот­ сгорания, а так­

(22)

/ 230 0 .

Частные производные, входящие в выражение (23), могут быть определены по результатам стендовых испытаний двига­ теля.

Таким образом, выражение (21) можно переписать в следу­ ющем виде:

Л

v0

/

dQ

, dQ

\ ,

v0

/

dT

, дТ

 

Q0

\

dati

dp

J

T0

\

dco,

dh

 

 

 

,

дТ

\ ,

v0

x 4 .

 

 

 

 

 

+ —— и

+

 

 

 

 

 

 

dp

J

co40

 

 

 

32

Подставив в выражение (17) зависимости

 

(19) — (21), полу­

чим

 

 

 

 

 

 

 

 

 

 

 

 

а

 

Л*™ /

dQ

 

,

dQ

\

,

 

 

Со

dv

Qo

V 5 ш 1

 

 

dp

 

 

/

fe0

Л

Г0

\

дях

 

dh

 

др /

 

Н ; — —г—-

 

 

Л 4 + Х4

 

 

 

Со

<?V

СО40

 

«40

 

 

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

 

 

 

 

 

 

АМТ

= аххх

+

а2«/ +

а3и + а4лг4,

 

(24)

где

л*,*

<JQ ,

м т о

 

ас

<э<э

 

 

 

 

 

 

 

 

 

Qo

да);

 

Co

 

dv

(Эй)!

 

 

 

 

a2--+

мто

 

dt,

v0

ar

 

 

 

 

 

Co

 

<?v

r 0

 

aw,

 

 

 

 

 

MT(>

 

dl

v0

(ЭГ .

 

 

 

 

=0

 

dv

T0

dh

 

 

 

 

 

Mro

dQ

,

M x o

_

aC

v0

_

dQ

_|_

 

 

 

 

_|_

---iu

dv

 

 

 

Qo

<?p

 

Co

 

Q0

 

dp

 

 

a+4 =

,

 

dl

dT

 

 

 

 

 

 

Co

 

dv

dp

 

T0

 

 

 

 

 

M T 0

dC

 

v0

MT

 

 

 

 

Co

dv

 

G)40

 

(o4 0

 

 

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

определить по его

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

которая

 

обычно задается

ввиде зависимостей

е= — = e(Q, (о4) и % = r\(Q, со4). Ро

Используя известные соотношения, получим

MK=J-QHK-±-

(25)

Г

(26)

Я к = - ^ - К е Г о ! . в * - 1

к1

 

где k — показатель адиабаты;

т]„ — адиабатический к. п. д. компрессора.

3 Зак . 830

33

Эти характеристики могут быть перестроены в виде функции M K =M(Q,co 4 ) .

Приращение момента может быть найдено из разложения

 

 

 

dQ

 

 

 

dat

 

 

 

 

дМк

dQ

 

д.Мк

 

dQ

 

дМк

 

х4-

 

dQ

• ^—•*! " т

 

— —

u "

i —

;

 

да>1

 

 

dQ

 

dp

 

d<»4

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

АМК

 

= а5х{

+ а6и + а7х4,

 

 

 

(27)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а5

=

dMK

 

dQ

 

 

 

 

 

 

 

 

dQ

 

доз,

'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dMK

 

dQ

 

 

 

 

 

 

 

 

 

 

dQ

 

dp

'

 

 

 

 

 

 

 

a7

 

 

dMK

 

 

 

 

 

 

 

 

 

 

d(i)t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

(24)

и

(27) в уравнение

(14), получим

 

dx

 

 

 

 

 

 

 

 

 

 

 

J—-i *-at

= clxl

 

+ a2y + c2u+csxi,

 

 

(28)

где

 

 

 

 

 

 

 

 

 

 

 

 

сх = ах — аъ; с2 = а3—а6; с3 = а4—а7.

 

 

Для определения плотности воздуха за компрессором поль­

зуются известным

выражением

 

 

 

 

 

 

 

 

I O V k

_

=

 

 

 

low

 

 

 

 

( 2 9 )

 

 

 

Rer.

 

 

 

 

 

 

 

Универсальную

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

компрессора

можно

пере­

строить в координатах:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p = p(Q,

со4).

 

 

 

 

(30)

Исходя из функциональной зависимости

(30),

 

 

 

 

 

dQ

 

 

 

<?со4

4

 

 

 

 

или, подставив AQ из выражения

(19),

после

некоторых

преоб­

разований имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u = Llxl

+ L2xi,

 

 

 

 

(31)

34

где

dp

dQ

dp

 

dQ

dw,

dco4

 

dp

dQ

dp

dQ

~~

dQ '

dp

dQ '

dp

Уравнения движения

валов

гидротрансформатора. Струк­

турная схема САР при наличии в трансмиссии

гидротрансформа­

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

 

 

Структурная схема составлена из условий, что динамическая схема МТА при установившемся режиме является двухмассовой.

 

Регулятор

и,

 

 

 

 

 

 

 

 

 

 

[Двигатель

Насосное

Турбинное Мс

 

 

колесо

 

колесо

 

 

 

 

УГоплибный.

 

 

Реакторное

1

насос

 

 

 

колесо

 

 

 

 

Рис. 12. Структурная схема САР МТА при гидро­ механической трансмиссии трактора

При наличии гидротрансформатора и жестком соединении вала двигателя и вала насосного гидротрансформатора уравне­ ние (5) двигателя можно записать так:

 

 

 

 

dt

 

 

 

 

 

(32)

 

 

 

 

д

 

 

 

 

 

где

/ \

— суммарный

момент

инерции

двигателя

и

насосного

 

Мн

колеса

гидротрансформатора;

 

 

 

 

 

— момент

на

валу насосного

колеса

гидротрансформа­

 

 

тора.

 

 

 

 

 

 

 

 

 

Исходя из того, что

 

 

 

 

 

 

 

 

 

 

— Мн

= М„ +

Мр,

 

 

 

уравнение (32) можно представить в следующем

виде:

 

 

 

 

 

do),

: Мд + М„

+ М р ,

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

где

Мгг — момент на

валу турбинного

 

колеса

гидротрансформа­

 

 

тора;

 

 

 

 

 

 

 

 

3*

Мр — момент на реакторе гидротрансформатора.

35

Перейдя от моментов к их отклонениям, можно записать

 

j\J^i-

= АМД

+ АМГ Т

+ АМР .

(33)

 

 

dt

 

 

 

 

Разложив зависимость (33) в ряд Тейлора и ограничившись

первыми степенями переменных, запишем

 

dt

(?С0[

dh

 

дют

дщ

 

,

дМр .

,

dMv

Л

 

 

-\

2- Ай)х Н

— Дсйь

 

или

 

<Зозх

 

(Эо>1

 

 

 

 

 

 

 

 

 

/ ;

^ = Л ^ 1

+

Л 2 (/+Л 4 д; 2 ,

(34)

где Х2 = Дсот приращение угловой скорости вала турбинного колеса гидротрансформатора или вала муфты сцепления;

д' _ сШд , дМгт , дМ0

д

 

_

сШг т

, дМс

4

Ц Л - 1 Г Т

_ ^

 

:

дшт

да>т

Таким образом, уравнение (34) описывает движение вала двигателя и вала насосного колеса гидротрансформатора. Этим же уравнением описывается движение вала насосного колеса

гидромуфты, но в коэффициентах Л' и Л 4 не будет членов ... М р

дМр

В соответствии со структурной схемой САР уравнение движе­ ния вала турбинного колеса гидротрансформатора

М Г Т - / Т ^ - = М С ,

 

 

 

 

dt

 

 

где Мс

— момент

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

действующий

на вал турбин­

 

ного колеса;

 

 

 

 

 

/ т

— момент

инерции вращающихся

частей

турбинного ко­

 

леса, включая момент инерции жидкости гидротранс­

 

форматора.

 

 

 

 

 

После разложения,в ряд Тейлора первого члена предыдуще­

го выражения получим

 

 

 

 

 

 

дЛ1г т

А

,

дМгт А

. dcoT

. . .

,,.

 

дшг

 

 

дш^

dt

 

 

или

 

 

 

 

 

 

 

 

-JT-2£-

dt

+ A6xl

+ Aex2=bMc{t),

(35)

 

 

 

 

 

 

 

36

где

. _ дМгг . , _ дМгг

Уравнение движения вала турбинного колеса гидромуфты не отличается от уравнения (35).

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

I

модель — трактор с механической трансмиссией,

двигатель

со свободным впуском;

 

 

 

 

I I

модель — трактор с механической трансмиссией,

двигатель

с газотурбинным

наддувом;

 

 

 

I I I модель — трактор

с

гидромеханической

трансмиссией,

двигатель со свободным

впуском.

 

 

I

модель включает уравнения (8) — (11). I I модель

включает

уравнения

(13а),

(9) — (11),

(28), (31). I I I модель

включает

уравнения

(34), (35), (9) —(11).

 

 

 

 

 

 

5. Методы исследования

 

 

 

 

 

 

математических моделей работы МТА

Возможности

электронного моделирования.

В зависимости

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

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

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

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

дования

систем следует пользоваться сочетанием этих

методов.

В этом

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

оптималь-

37

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

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

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

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

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

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

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

38

элементы системы, так и систему в целом, наблюдать на экране осциллографа или регистрировать колебания на входе в любой элемент системы и выходе из него.

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

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

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

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

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

6. Методы определения исходных данных для моделирования

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

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

39

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