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

2328

.pdf
Скачиваний:
25
Добавлен:
07.01.2021
Размер:
5.63 Mб
Скачать

б)

Рис. 3.3

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

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

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

Основным документом, характеризующим механические и экономические (по расходу топлива) показатели двигателя внутреннего сгорания, является скоростная характеристика двигателя, дополненная перечнем значений его основных режимных параметров

[18, 27, 58, 71] (табл. 3.1).

Таблица 3.1

 

Основные режимы работы двигателя

 

 

 

Точка кривой

 

 

крутящего

Название режима работы

Режимные параметры

момента

двигателя

 

двигателя

 

 

О

Режим

максимальных

оборотов

n n

max

; Me 0

холостого хода

 

 

1

 

 

 

 

1

Номинальный

(максимальной

n nН ; Me МeН ;

мощности)

 

 

Ne

Nemax

 

 

 

2

Режим

среднего

 

крутящего

n n2; Me Мe2;

момента

 

 

 

Memax Me2 MeH

 

 

 

 

3

Максимального крутящего момента

n nМ max;

двигателя

 

 

Me Мemax

 

 

 

4

Режим минимальных

устойчивых

n nМ min;

оборотов

при

минимальном

Me Мemin

 

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

 

5

Режим

минимальных

оборотов

n n

min

; Me 0

холостого хода

 

 

 

 

 

 

 

На рис. 3.4 показана скоростная характеристика двигателя, снабженного всережимным регулятором, представляющая собой зависимость крутящего момента на коленчатом валу двигателя от частоты вращения, выраженной в оборотах в минуту: Me f (n). Такие характеристики получают в результате стендовых испытаний двигателей, регламентированных ГОСТ 14846-81 (СТ СЭВ 765-77). Одновременно в результате испытаний определяют часовой расход топлива как функцию частоты вращения коленчатого вала: GT f (n). В технической литературе и в рекламной информации обычно отсутствуют сведения о точке 4 характеристики двигателя. Точки 6 10 зависимости часового расхода топлива GT f (n) соответствуют рассмотренным режимам работы двигателя (см. табл. 3.1).

На рис. 3.5 показана скоростная характеристика двигателя, дополненная вычисляемыми величинами Ne Ne(n) и ge ge(n).

Рис. 3.4

Рис. 3.5

Эффективная мощность двигателя определяется по выражению Ne Me e , где e – угловая скорость вала двигателя. В мировой практике двигателестроения принято частоту n вращения вала двигателя измерять в об/мин. При вычислении мощности удобно угловую

скорость e выражать через частоту вращения n в об/мин по формуле

 

 

 

 

Ne

Me n

 

.

(3.16)

30

 

 

 

 

Удельный расход топлива вычисляют по формуле

 

ge

GT

1000,

(3.17)

 

 

Ne

 

где ge – удельный расход топлива, г/кВт ч; GT – часовой расход топлива на данном режиме работы внешней скоростной характеристики, кг/ч; Ne – мощность двигателя, кВт.

В настоящее время в связи с внедрением энергосберегающих приводов появилась необходимость выполнять расчеты топливноэкономической эффективности двигателей. В этих случаях оказывается недостаточной аппроксимация кривой крутящего момента тремя точками (1, 3, 4), как это принято в настоящее время. Для повышения точности расчетов, связанных с созданием энергосберегающих технических систем, необходимо между точками 1, 3 вводить несколько дополнительных точек 2 (21, 22 , 23 и т. п.). Рассмотрим аналитический метод аппроксимации характеристик двигателя любым количеством точек на примере использования четырех точек 1 4.

Запишем выражение зависимости мощности двигателя в функции частоты вращения в виде многочлена четвертой степени [18, 27]

 

 

 

 

 

 

n

 

n

 

 

 

2

 

 

 

 

n 3

n

4

 

 

 

 

 

Ne=Nemax a1

 

 

a2

 

 

 

 

 

 

 

 

a3

 

 

 

 

 

 

a4

 

 

 

 

,

 

(3.18)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nH

 

nH

 

 

nH

nH

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где a1, a2 , a3, a4

– постоянные коэффициенты.

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая,

что

N

e

=M

e

 

e

;

 

 

 

N

emax

M

eN max

 

 

 

;

 

e

n

;

 

 

 

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

30

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

eH

nH

, после преобразований из (3.18) получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

30

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

n

2

n

 

3

 

 

 

 

 

 

 

Me =MeN max a1

a2

 

 

 

 

 

a3

 

 

 

 

 

a4

 

 

 

 

.

 

(3.19)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nH

 

 

nH

 

nH

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким

образом,

 

кривые

крутящего

 

момента Me

 

и мощности

двигателя Ne аппроксимируются с помощью единой системы аппроксимационных коэффициентов a1,…,a4 , для определения которых достаточно составить систему четырех уравнений типа (3.19) для точек 1 – 4 кривой крутящего момента двигателя (рис. 3.5).

Можно заметить,

что для точки 1 при n nH ; Мe =МeN max имеем

a1+a2 +a3+a4 =1.

Определять коэффициенты аппроксимации для

системы четырех уравнений типа (3.19) удобно с помощью Microsoft Excel.

В табл. 3.2 приведены результаты совместных решений систем уравнений для точек 1, 2, 3, 4 внешних ветвей скоростных характеристик двигателей.

Таблица 3.2

Коэффициенты аппроксимации крутящего момента

Марка

Номинальная

Коэффициенты аппроксимации кривых

крутящего момента и мощности двигателя

двигателя

мощность, кВт

a1

a2

a3

a4

 

 

А-01М

121,5

0,700885

0,153943

1,974576

–1,829404

ЯМЗ-236Н

175,0

0,967993

–0,348171

1,611048

–1,230870

SA6D155-4A

305,0

–0,260938

4,577909

–4,074223

0,757252

В-84М

647,2

0,008573

4,246418

–5,017605

1,762614

Для двигателя ЯМЗ-236Н уравнения крутящего момента и мощности двигателя имеют вид

 

 

 

 

n

 

 

 

 

n

2

n

3

Me =MeN max 0,9679 0,3481

 

 

 

 

1,6110

 

 

 

 

1,230

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

nH

 

nH

 

 

nH

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

n 2

 

n

 

3

n

 

 

4

Ne=Nemax 0,9679

 

 

0,3481

 

 

 

 

1,6110

 

 

 

1,230

 

 

 

.

 

 

 

 

 

 

 

 

 

nH

 

nH

 

nH

nH

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Функция часового расхода топлива двигателем является монотонной кривой второго порядка, поэтому ее достаточно аппроксимировать с помощью трех точек:

 

 

n

 

n

2

GT =GTN max b1

b2

 

 

b3

 

 

,

 

 

 

nH

nH

 

 

 

 

 

 

 

 

 

где b1, b2 , b3 – постоянные коэффициенты.

Для точки 7 (см. рис. 3.4), соответствующей номинальному режиму максимальной мощности, при n nH ; GT =GTN max имеем уравне-ние b1+b2 +b3=1.

В табл. 3.3 показаны значения коэффициентов bi для разных двигателей. Для двигателя ЯМЗ-236Н часовой расход топлива на внешней ветви скоростной характеристики определяется уравнением

 

 

 

 

n

 

n

2

GT =GTN

 

0,068579 1,48145

 

 

0,412871

 

 

.

 

 

 

 

max

 

nH

nH

 

 

 

 

 

 

 

 

 

 

 

Таблица 3.3

Коэффициенты аппроксимации часового расхода топлива

Марка

Номинальная

Коэффициенты аппроксимации кривых

часового расхода топлива двигателя

двигателя

мощность, кВт

b1

b2

b3

 

 

А-01М

121,5

–0,702022

3,511149

–1,809127

ЯМЗ-236Н

175,0

–0,068579

1,481450

–0,412871

SA6D155-4A

305,0

–0,193318

1,8729045

–0,679586

В-84М

647,2

0,333129

0,666870

0,000001

Увеличение момента сопротивления Мс на валу двигателя, например, в точке 2 приводит к уменьшению числа оборотов и увеличению крутящего момента двигателя, т.е. в диапазоне изменения частоты вращения nM max n nH двигатель автоматически приспосабливается к изменению нагрузки, работает в устойчивом режиме.

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

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

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

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

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

задач и характера динамических процессов нагружения двигателя [70, 71].

На рис. 3.6 показана внешняя скоростная характеристика дизельного двигателя В2-550ТК-С5 фронтального погрузчика ТО-21 мощностью 406 кВт, которая устанавливает связь выходных координат Ме, GT, Ne, ge, ηe на установившихся режимах работы с частотой вращения n коленчатого вала.

Ме, Н·М

 

 

 

Ne, КВТ

Me

 

 

 

 

 

 

 

400

 

 

 

 

 

2500

 

 

 

380

2000

Ne

 

 

360

 

 

 

 

 

1500

 

 

 

340

 

1000

 

 

 

320

 

500

 

 

 

300

GТ, КГ/·Ч

 

GT

 

280

 

80

 

 

ge, г/КВТ·Ч

 

 

 

 

ηе

70

ge

 

 

240

0,5

60

 

 

200

 

ηe

 

0,4

50

 

 

160

 

 

 

0,3

40

1300

1500

1700

n, об/мин

 

1100

 

 

 

Рис. 3.6

 

 

Как показали результаты, изложенные в подразделах 3.1–3.4, дизельный двигатель воспринимает низкочастотные входные воздействия в диапазоне частот fп=0 5 Гц. Рабочие процессы фронтальных погрузчиков при наполнении ковша имеют тренды рабочих сопротивлений, которые характеризуются периодичностью выполнения

Тк =4 8 с [73, 77].

На рис. 3.7 показана примерная

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

погрузчика.

В разделе 11 будет рассмотрен

Рис. 3.7 аналитический метод расчета таких циклограмм, которые характеризуются периодом Тк процесса набора ковша: Тк=3 6 с. При этом полное время одного цикла погрузчика может составлять Тц 30 40 с и более. Технологические рабочие процессы фронтальных погрузчиков можно разделить на быстро протекающие и длительно протекающие динамические процессы. В соответствии с этим динамические процессы, выполняемые дизельными двигателями, тоже являются длительно протекающими или быстро протекающими динамическими процессами.

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

Пусть момент сопротивления на валу двигателя на некотором интервале времени изменяется по линейному закону

Mc At,

где А – интенсивность изменения момента сопротивления двигателя, Нм/с.

На рис. 3.8 показан переходный процесс крутящего момента двигателя при линейном законе нагружения с интенсивностью

A Me . Интеграл крутящего момента двигателя

t

 

 

1

t

(3.20)

 

 

Me A(t ) A e .

Выражение (3.20) описывает переходный процесс крутящего

момента двигателя Ме

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

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

увеличением времени

процесса второе

слагаемое уравнения (3.20) быстро уменьшается, поэтому момент на валу двигателя Ме превращается в линейную зависимость

Mc, Ме

 

Me A(t ),

 

 

 

 

 

 

 

 

 

 

 

 

 

Mc

где – время запаздывания двигателя.

 

 

 

 

Поскольку в

данной

работе

 

 

 

 

 

 

 

 

 

 

 

Ме

рассматриваются

длительно

 

 

 

 

 

 

 

 

 

At

 

протекающие динамические процессы,

 

 

 

 

 

 

то временем запаздывания двигателя

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

можно

пренебречь

как

малой

 

 

 

 

 

 

величиной по сравнению со временем

 

 

 

 

 

 

протекания

технологических

 

 

 

 

 

 

 

 

 

 

 

 

процессов подъема, опускания стрелы

 

 

 

 

 

 

и процессов копания.

 

 

 

 

 

Рис. 3.8

 

 

 

 

 

 

 

 

 

 

Для

составления математических

 

 

 

 

 

 

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

В работах профессора Н.А. Ульянова [71] связь выходных параметров двигателя е и Ме устанавливается в аналитическом виде по линеаризованным статическим характеристикам. В последующем таким методом представления внешних характеристик двигателя с некоторыми усовершенствованиями пользовались М.В. Коваленко [39], А.Н. Подсвиров [48] и другие исследователи. Однако, как показал опыт использования этих моделей, они обладают недостатками, обусловленными громоздкостью записи соответствующих формул для всех линеаризованных участков и условий для их ограничения.

Перспективным для задания внешних характеристик дизельных двигателей является предложенный табличный метод [10].

Суть метода состоит в том, что любые функциональные зависимости, например Ме=Ме( е); GТ=GТ( е), задаются линеаризованными отрезками, концы которых являются узлами этих зависимостей и характеризуются соответствующими координатами, которые систематизируются в виде матрицы. Для получения текущих промежуточных величин при исследовании непрерывных процессов используют метод аппроксимации, основанный на полиномах Лагранжа.

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

многочленов Лагранжа.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

функции yi

Сначала

рассмотрим

 

простой

случай

 

выражения

 

одной переменной

семейством

точек

xi ,

принадлежащих кривой

y f (x).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общий вид формулы интерполяционного многочлена Лагранжа

степени n, принимающего значения yi

в точках xi

(i – целое число,

изменяющееся в пределах от 0 до n),

имеет вид [10]

 

 

 

 

 

 

 

 

 

 

n

(x x

0

)...(x x

 

)(x x

i

1

)...(x x

n

)

 

 

 

 

L

П

(x)

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

y

i

.

(x x )...(x

 

 

)(x

x

 

 

)...(x x

 

)

 

i 0

x

 

 

 

 

n

 

 

 

 

 

i

0

i

i 1

 

i

 

i 1

i

 

 

 

 

 

 

 

В зависимости от вида кривой используют интерполяционные многочлены Лагранжа первой или второй степени, т.е. кусочнолинейную или кусочно-квадратичную интерполяцию. В узлах интерполирования xi функция f (x) должна совпадать со своим интерполяционным многочленом. Если же x не является узлом интерполирования, то f (x) только приближенно равно LП (x).

Рассмотрим алгоритм, позволяющий запрограммировать

функцию

z двух переменных

 

 

z f (x,y)

дискретным

семейством

точек xi ,

y(K ). Для текущего значения y(K )

находим место этой точки

между кривыми K

и K 1.

Для этих кривых выполняем сечение по

параметру x

const

и находим

 

 

zK (x )

 

 

и

zK 1(x ).

Далее,

 

применяя

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

i

 

 

 

 

 

 

 

 

интерполяцию,

 

 

находим

 

 

 

искомое

 

 

значение

 

 

z f (x,y).

Математическая запись этого алгоритма имеет вид

 

 

 

 

 

 

 

 

zK (x)

(x xi )(x xi 1)

 

 

 

 

z

i 1

 

(x xi 1)(x xi 1)

z

i

 

 

 

 

 

 

 

 

(x

x )(x

x

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x x

)(x

i

x

)

 

 

 

 

 

i 1

 

i

i 1

 

 

i 1

 

 

 

 

 

 

 

 

 

i

i 1

 

 

 

i 1

 

 

 

 

 

 

 

 

(x xi 1)(x xi )

 

 

z

i 1

,

 

 

 

 

 

 

 

 

(3.21)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x

x

 

)(x

i 1

x )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

i 1

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где x [xi 1,

xi 1], i=1, 3, …, n–1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z (y,x)

 

 

(y y(K))

 

z(K 1)(x)

(y y(K 1))

 

z(K)(x), (3.22)

(y(K 1) y(K)

(y(K) y(K 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где y [y(K),

y(K 1)], K=1, m.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В данной работе выполнено обоснование правомерности данного подхода к решению рассматриваемой динамической задачи.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]