Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
2501.pdf
Скачиваний:
45
Добавлен:
07.01.2021
Размер:
13.38 Mб
Скачать

2. Бараз, В.Р. Корреляционно-регрессионный анализ связи показателей коммерческой деятельности с использованием программы Excel: учебное пособие / Екатеринбург: ГОУ ВПО «УГТУ – УПИ», 2005. – 102 с.

3. Чижонок В.Д. Теоретические основы и практические приложения логистики / В. Д. Чижонок. – М : Новое знание, 2015. – 320 с.

4.

Логистика: интеграция и оптимизация логистических бизнес–процессов в

целях поставок / В.В. Дыбская [и др.]. – Москва: Эксмо, 2014. – 939 с.

5.

Волгин, В. В. Логистика приемки и отгрузки товаров: практическое пособие /

В. В. Волгин. – Москва: Дашков и Кº, 2009. – 457 с.

6.

Герас мов, Б.И. Основы логистики / Б.И. Герасимов, В.В. Жариков, В.Д.

Жариков. - М.: ИНФРА-М, 2010. – 304 с.

7.

Голубч к, А.М. Транспортно-экспедиторский бизнес: создание, становление,

 

/ А.М. Голубч к. – Москва: ТрансЛит, 2011. – 317 с.

8.

Иванов, Д.А. Управление цепями поставок / Д.А. Иванов. – Санкт-Петербург:

С

Издательство Пол техн ческого университета, 2010. – 659 с.

9.

Интегр рованные логистические системы доставки ресурсов: (теория,

методолог я, орган зац я) / И.А. Еловой, И.А. Лебедева. – Минск: Право и

экономика, 2011. – 460 с.

управление

 

бА

УДК 629.3

ВЛИЯНИЕ ПОГРЕШНОСТЕЙ ИЗМЕРЕНИЙ НА ПАРАМЕТРЫ КОЛИЧЕСТВЕННЫХ СТРУКТУР

В.А. Щеглов, канд. техн. наук, доц., профессор кафедры АТиСА,

БГАРФ ФГБОУ ВО Д«КГТУ», г. Калининград

V.A. Shcheglov, Ph.D., associate Prof., Рrofessor of road transport and service vehicles, «BFFSA» FSEI HPE «Kaliningrad State Technical University», Kaliningrad

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

Abstract. Evaluated the impact of measurement errors on the parameters of the quantitative structures and the requirements to the instrument requirements. To obtain a stable solution, it is necessary to conduct the pilot measurement is not less than 20 modes.

Ключевые слова: погрешности, количественные структуры, входные параметры, выходные параметры, системы.

Keywords: error quantity structure, input parameters, output parameters, systems.

272

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

деление погрешностей на разновидности. СПогрешности могут быть разбиты на несколько групп [1]:

1. Погрешности измерения физических величин приборами.

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

3. Погрешности, связанные с постановкой математической задачи (погрешность задачи ли метода).

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

Разработка метод ки оценки погрешности измерений. Основные

способы установлен я

пределов допускаемых

погрешностей и

обозначен я классов точности измерения

установлены ГОСТ 8.401-80

(классы точности средств измерений). Погрешности измерения

формируются инструментальными и методическими погрешностями.

Основная предельная погрешность измерений δи входных и выходных

параметров рассчитана по формуле:

 

 

 

δи = δh + δn+ δmt ,

(1)

 

 

Д

где δh , δn , δt , δm - погрешности измерений соответственно хода рейки

ТНВД, частотыбАвращения вала ТНВД, времени замера навески и цикловой

подачи топлива (навески топлива).

 

 

 

Основная предельная погрешность измерений входных и выходных

параметров δи = 2,5 + 0,3 + 0,5+0,3 =3,6; δи = 3,6 %.

 

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

погрешности

И

δтп

минимизируются

множественностью реализаций и составляют 3 %.

 

Погрешность округления согласно данным работы [2] и при строгом

выполнении рекомендаций по округлению чисел, не превышает 0,1 %.

Оценка погрешностей, возникающих в результате флуктуации

теплофизических свойств топлива может быть оценена следующим

образом.

 

 

 

 

Физико-химические показатели топлив определяются динамической

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

топлива,

η = νγ.

 

(2)

 

 

Известно, что плотность жидкости γ изменяется нелинейно от температуры и для одного вида топлива при изменении температуры t в

273

пределах 10 - 50º С (рабочий предел изменения температуры топлива при работе ТА) плотность уменьшается не более, чем на 7 % [3].

Также указывается, что расход жидкости при турбулентном движении, которое, несомненно, имеет место в насосе и в форсунке,

практически не зависит от вязкости. При ламинарном движении, подчиняющемся закону Пуазейля [2] на объем протекающей жидкости непосредственно влияет динамическая вязкость.

Снесущественно. При проведении экспериментов температура топлива в баке установки поддерживалась за счет расположения установки в отапливаемом помещен в пределах 20 – 22º С, что позволяет сделать вывод с учетом рекомендаций ра оты [1] об уровне погрешности на

Таким образом, при существенном изменении кинематической вязкости в зав с мости от температуры плотность его изменяется

которыйПридает погрешность не олее ±0,5 % при применении электронных весов соответствующей точности, в то время как второй способ зависит от многих факторов (температура, плотность топлива и др.) и на разных стендах в зависимости от конструкции мерительного оборудования дает погрешность ±(0,5 - 3) %.

порядок н же, т.е. до 0,7 %.

ПогрешностьбАзадачи или метода δм. Начальная погрешность при принятой в работе методике построения количественных структур формируется за счет ограничений по аппроксимации экспериментальных данных полиномами второй степени. Это ограничение вызвано

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

необходимостью

ограничений

 

по

громоздкости

получаемых

зависимостей. Кроме того, полином второй степени еще способен

показать несистематические ошибки. Оценка точности аппроксимации

выполнена по коэффициенту корреляции

R2, который во всем поле

 

 

Д

режимов работы ТА не превышает 2 %.

 

 

С учетом потребностей эксперимента и методики его проведения

решена обратная задача определения суммарной погрешности при

установленной из условий решения общей задачи предельной

погрешности δu

функции u = f(x1,

x2, x3,… ., xn). Критерием успешного

 

 

 

 

И

проведения эксперимента и обработки его результатов является зависимость δΣ = δм < δu.

Решение обратной задачи выполнено по принципу равных влияний. Согласно этому принципу предполагается, что все частные дифференциалы одинаково влияют на образование общей абсолютной погрешности δu функции общего вида Di = f(h, n ):

274

 

df

.x

(i = 1, 2, …, n) .

(3)

 

 

 

 

i

 

 

dxi

 

 

В нашем случае в качестве таких функций принимаем

количественные структуры комплектов ТА.

 

Величина предельной абсолютной погрешности δu может быть задана

исходя из анализа работ [4] в которых в которых указывается на слабую

С

 

 

 

 

 

 

 

связь изменения функциональных и структурных параметров ТА в

пределах 0 - 5 % с выходными энергетическими и экономическими

показателями рабочего процесса дизеля. Принимаем относительную

предельно допуст мую погрешность δu

= 6 - 12 %. Конкретные значения

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

 

 

 

 

 

 

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

отклонен ю экспер ментальных точек от прямой. Таким образом, δu

определяет односторонний

 

доверительный

интервал

и

степень

 

каждого эталона.

 

 

 

 

 

На основан :

 

 

 

 

 

 

 

 

δu = Σ |du/dxi | δxi .

 

 

 

(4)

Полагая, что все слагаемые равны между собой, будем иметь:

 

Отсюда:

|du/dx1 |δx1 = |du/dx2|δx2 =…= |du/dxn |δxn = δu/n.

(5)

δxi = δu / (n |du/dxi|)

(i = 2).

 

 

(6)

 

 

 

Расчет допустимой погрешности измерения входных параметров (n,

h) позволил

получить значения предельной

погрешности

измерения:

h=1,0 %, n=2,2 %. Таким

о разом,

сформулированы

требования к

приборному обеспечению.

 

Д

 

Оценка

бАвлияния погрешностей измерения

информативных

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

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

Положим, что взаимосвязи между информативными параметрами заданы в виде суммы произвольных линейно зависимых операторов. На основе исходных данных из пространства, принятого за эталон, для определения неизвестных х1, х2, х3, …, хn образуем систему условных линейных однородных уравнений:

275

А х 0,

(7)

и эквивалентную ей систему нормальных уравнений:

 

АТ А х 0.

(8)

Из однородного характера системы (8) следует, что длина искомого

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

постоянного множителя. Это обстоятельство используем для

С

 

 

 

хn

= -1,0,

целесообразного нормирования вектора х. Например, полагая

запишем систему (8) в неоднородной форме

 

 

 

 

 

 

А

х f ,

 

(9)

где

A

- с мметр чная матрица порядка (n-1) (n-1).

 

 

рассмотрим

 

АТ А

 

Если элементы матр цы

А заданы точно и дефект матрицы

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

АТ А построена на основе параметров,

 

измеренных с погрешностями.

Поэтому оцен

х вл яние на устойчивость решения.

 

 

 

 

 

 

 

Наряду с (9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

возмущенное уравнение:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(А

 

 

)

 

 

f

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где А

и f

- соответственно возмущения матрицы

 

 

 

 

 

 

и вектора f.

 

 

 

A

 

 

 

 

 

 

 

б

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Предположим, что матрица A - невырожденная, поэтому существует

обратная матрица

 

 

1.

 

При

возмущении

 

 

 

 

 

на

 

 

 

 

 

 

 

 

 

обратная матрица

 

 

 

A

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

А

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

Д

 

 

 

получит возмущение (

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обозначим через:

 

А

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

) 1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

,

 

(11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

 

 

 

 

А

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

 

 

 

И

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

величины относительных возмущений матриц A и

 

 

 

A 1.

Как показано в

(15), справедлива оценка:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Н А

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 Н А

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Н

 

 

 

 

 

 

 

 

 

 

 

 

 

- число обусловленности матрицы

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

А

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Оценка

(12)

показывает,

что малое

 

относительное

возмущение

 

 

 

A

приводит к малому относительному возмущению A 1 только в том случае, кода число обусловленности матрицы A не слишком велико сравнительно с единицей.

276

Введем обозначения для относительных возмущений решений х и вектора f:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х

 

 

х х

 

;

f

 

 

 

f

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если величина

 

 

возмущения

 

 

 

 

 

 

такова,

 

что

матрица

 

 

 

 

f -

 

 

 

 

 

À

А

невырожденная, то возмущенное уравнение (10) и точное (9) будут иметь

единственные решения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

 

) 1 (f

 

);

 

 

x

 

f .

 

 

 

 

 

 

 

 

 

 

(14)

 

 

 

 

 

х

А

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

 

 

 

 

f

погрешности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из (14) наход м:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сx x [(E A 1

) 1 E]

A

1 f (E

A

1

 

) 1

A

1 .

(15)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

f

 

 

 

 

 

 

 

В номерах матр цы

 

 

 

A 1

 

и

вектора

 

 

f

 

величина абсолютной

 

преобразований

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

решен

я составит:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х

х

 

 

(Е А 1

) 1 Е

 

 

 

х

 

 

 

(Е

А

1

 

) 1

 

 

 

А

1

 

 

 

 

 

 

.

После

 

 

 

 

 

 

 

 

 

 

А

 

 

 

 

 

 

 

 

 

 

(16)

 

 

 

 

 

 

 

 

 

 

 

 

 

(16) для относительного возмущения решения

в (17) получена следующая оценка:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х

Н

 

 

 

(

 

f )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(17)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 Н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Д

 

 

 

 

 

 

 

Оценка (17) и устанавливает влияние погрешностей измерения

информативных параметров на устойчивость решений. С точки зрения

устойчивости решения важно,

 

чтобы число

обусловленности было не

слишком большим. Например, если использовать спектральную норму, то Н

равно отношениюмаксимального числа матрицы А кминимальному, т.е.:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

И

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

max

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Н

А

 

 

 

А

 

 

 

.

(18)

 

 

 

 

 

2

min

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из (17) и (18) следует, что малые относительные возмущения

матрицы

 

и вектора f вызывают

 

малое

относительное

возмущение

А

 

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

Следовательно, критической величиной, которая решает вопрос о физической надежности решения системы (9), является не определитель матрицы А, а отношение ее собственных чисел max и min. Квадратный корень этого отношения измеряет увеличение помех в направлении, соответствующем собственному числу min. Пока число Н не превышает

277

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

а1,а2 ,а3 ,...,аn решениенеявляетсякритическим.

Но если Н достигает величины 104 и более, увеличение помех в

направлении min

может привести к такому искажению решения, которое не

будетсоответствоватьфизическомусмыслузадачи[5].

 

 

Таким

образом, если

число

 

Н

принять в

качестве

меры

 

 

 

А

Т

 

 

 

 

 

обусловленности

матриц

А, то появляется

необходимость

А и

установить

пр знаки,

определяющие

обусловленность

матриц

относительно процесса вычисления собственных чисел. Важно отметить,

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

матриц относ тельно выч сления решений. Покажем это.

 

 

С

 

Т

 

Т

 

Т

 

Рассмотр м матр цу А А dA

А, близкую к матрицеА А. Пусть 1,

2, 3, …, n – собственные

 

 

и u1, u2, u3, …, un – собственные векторы

матрицы А

Т

А. Тогда:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т

 

 

Т

 

 

Т

 

 

 

числа

 

) (ui dui) ( i d i) (ui

dui).

(19)

А

А ui i ui;

(A

А dA

Пренебрегая в (54) величинами второго порядка малости, получим:

 

 

dA

Т

 

 

Т

dui idui d i ui.

 

(20)

 

 

А ui

A

 

 

 

б

 

 

Умножим (20) скалярно на ui

и учтем, что:

 

 

 

 

 

Т

 

dui;ui) i(dui;ui ).

 

(21)

 

 

 

(A

 

 

Тогда из (20) имеем:

 

 

 

 

 

 

 

 

 

 

 

Т

 

Т

 

(22)

 

 

 

d dA

; d dA .

 

 

 

 

Аi i

 

 

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

симметричных матриц.

И

 

 

 

 

 

 

 

 

 

 

 

 

 

Оценка степени линейной зависимости информативных параметров.

Теоретически

вопрос о

степени линейной зависимости

параметров

 

 

1 ,

 

2 ,

 

3 ,...,

 

n

решается

весьма просто: система (8) имеет

не равные

 

а

а

а

а

Т

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

278

насколько по своим свойствам построенная матрица близка к вырожденной.

Т

 

Действительно, определитель матрицы А А равен произведению ее

собственных чисел, т.е.:

 

Т

(23)

det A А 1 2 3... n .

Из приведенного тождества следует, что произведение собственных

чисел может быть весьма малым, если все множители i меньше единицы,

но больше погрешностей

змерения параметров а1 ,а2 ,а3 ,...,аn . Такое же

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

ближайшей

 

 

 

 

 

Т

 

 

множителей i

несколько весьма малых. Если учесть,

что расстояние до

Свырожденной матрицы определяется минимальным

собственным ч слом, то в первом случае условие линейной зависимости

 

 

 

 

 

 

 

Т

 

 

 

 

 

 

Т

 

 

больших

 

 

 

 

 

 

не выполняется, а во втором матрица А

 

А по своим свойствам близка к

вырожденной.

 

 

 

 

 

 

 

 

 

 

 

 

В

пр нц пе

спектр

со ственных

 

чисел

матрицы

А

А

может

соответствовать одному з следующих трех случаев [6]:

 

 

 

 

 

 

 

 

А

 

 

 

 

1. Среди чисел 1, 2, 3,

…, n имеются очень большие и малые.

Значения

 

 

чисел

заведомо

 

превышают

предполагаемый

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

уровень

погрешностей

 

элементов матрицы

А

А, а

значения малых чисел равны или меньше этого уровня.

 

 

 

 

 

 

 

 

 

 

Д

 

2. Нет четкого разрыва между собственными числами, т.е. каждое из

них сопоставимо с предполагаемым уровнем погрешностей

элементов

 

 

Т

 

 

 

 

 

 

 

 

 

 

 

 

матрицы А А.

 

 

 

 

 

 

 

 

 

 

 

 

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

числами, однако

значение min существенно

больше

предполагаемого

 

 

 

 

 

 

 

 

Т

И

уровня погрешностей элементов матрицы А А.

 

 

 

 

 

В

каждом

из названных

случаев

решение системы

нормальных

уравнений (8) имеет свои особенности.

Устойчивые решения системы (8) могут быть получены только в том

Т

случае, если в спектре собственных чисел 1, 2, 3, …, n матрицы А А (n-1) – чисел больше предполагаемого уровня погрешностей ее элементов и только число n равно или меньше этого уровня.

Случай, когда в спектре собственных чисел нет четкого разделения на большие и малые, соответствует условиям некорректной задачи. Здесь для определения устойчивого решения следует применять методы теории регуляризации [7].

279

Т

Если же все собственные числа матрицы А А существенно больше предполагаемого уровня погрешностей ее элементов, то это означает, что

параметры а1,а2 ,а3 ,...,аn не отвечают условию линейной зависимости и на их основе невозможно построить работоспособный диагностический эталон.

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

Действительно,

 

поскольку

параметры

 

а1,а2 ,а3 ,...,аn

заданы

с

погрешностями, то обычно полагают m n, что оправдывается

стремлен ем получ ть решение системы (8) с большей точностью. С

С

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

другой стороны,

m соответствует числу режимов, на которых следует

 

 

 

замеры параметров. Поэтому желательно, чтобы m было не

произвести

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ошибок

 

 

 

 

 

 

 

 

 

 

 

слишком больш м. Кроме этого, при увеличении m возрастает число

вычислительных операц й в процессе перехода от системы (7) к системе

(8) и, следовательно, накопление

 

 

 

 

 

 

 

 

 

 

в элементах матрицы А

А. В

конечном

 

счете

процесс

 

накопления

 

 

ошибок может

привести

к

 

 

 

 

 

 

 

 

 

А

 

 

 

существенному искажению решения.

 

 

 

 

 

 

 

 

 

 

 

 

 

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

обоснованно решить вопрос и о степени переопределенности условной

системы

(7).

 

Используем

 

для

 

 

 

этой

цели

известное

тождество,

в

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Д

Т

 

соответствии с которым сумма диагональных элементов матрицы

А

А

равна сумме ее собственных чисел, т.е.:

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

m

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

(ln a1i )2

(ln a 2i )

2 (ln a 3i

 

)

2 ... (ln a n 1,i )2

 

 

i 1

 

 

 

 

 

 

 

 

i 1

 

 

 

i 1

 

 

 

 

 

 

 

i 1

 

(24)

 

m

 

 

 

 

 

 

)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...

 

 

 

 

.

 

 

 

 

 

(ln a

n,i

2

3

 

n

1

n

 

 

 

 

 

i 1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из (24) следует, что между степенью переопределенности матрицы А

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т

А имеется прямая взаимосвязь.

и свойствами эквивалентной матрицы А

 

Действительно,

если параметры

 

1,

 

2 ,

 

3 ,...,

 

n

линейно зависимы (ранг

а

а

а

а

матрицы А

Т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А = n-1), то увеличение суммы членовИлевой части тождества

должно компенсироваться суммой собственных чисел 1, 2, 3, …, n-1

и

каждое

из

 

этих

чисел

должно

 

 

быть больше

предполагаемого уровня

Т

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

Т

элементах матрицы А А.

280

Изложенными

соображениями и определяется необходимость

 

 

 

 

 

Т

 

 

 

 

 

 

А и А.

определения и анализа спектра собственных чисел матриц А

 

 

 

Т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отметим, что в силу построения матриц А А и А их собственные

числа всегда вещественные и положительные или равны нулю, а их

значения определяются в результате решения уравнений

 

 

 

 

 

С

Т

 

 

 

 

 

 

 

 

 

det( A

А Е ) 0;det( A E ) 0 .

(25)

Располагая спектрами собственных значений матриц

Т

 

 

А А и А,

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

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

нормальной с стемы уравнений (9), степени линейной зависимости информат вных параметров, т.е. получить систематический ответ на вопрос о пр годности ото ранных информативных параметров и

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

Выясн м мое количество уравнений условной системы для

Весь массив исходных данных (m=72) разделим на группы, в каждую из которых включим замеры параметров через каждые 10 режимов. На каждом шаге вычислений для количественнойДоценки свойств матриц

Т

C1 h2 n2 C2 h n2 C3

n2

C4

 

2

 

 

 

 

 

h

n

(26)

C

 

 

 

C

 

C

 

2

C

 

 

C

 

 

 

0.

h

n

n

h

h

g

ц

необход5 6 7

8

9

 

 

 

 

 

 

 

 

 

А

 

 

 

 

 

 

А А и А используем спектры их собственных чисел.

Результаты вычислений коэффициентов Сi в функции от степени переопределенности условной системы уравнений приведены на рис. 1.

Изложенный анализ показывает, что на основеИпараметров n, h и gц может быть построена количественная структура типа (25).

281

СР . 1. Изменен е коэффициентов количественных структур

 

 

в зав с мости от степени переопределенности

 

Из р с. 1 в дно, что при m>20 изменения значений Сi сопоставимы с

 

сходных данных, поэтому их можно считать достаточно

погрешностями

 

хорошим пр

жен ем к решению системы:

 

 

х f .

(27)

Степень

сходимости экспериментальных и расчетных

данных в

бл зависимости от степениАпереопределенности показана на рис. 2.

Д И

Рис. 2. Сходимость экспериментальных и расчетных данных в зависимости от степени переопределенности

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

282

Влияние масштабных факторов на устойчивость решений нормальной системы уравнений. Трудности получения устойчивых решений системы уравнений (8) порой без особой надобности усиливаются неподходящим выбором масштабов для измерения информативных параметров. Действительно, уравнение (9) можно трактовать как задачу разложения заданного вектора f по осям заданной системы координат, определяемым

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

столбцами

 

матрицы

А. Определитель

 

А

 

 

 

имеет следующую

геометрическую

интерпретацию.

Он

 

 

представляет

 

объем

параллелеп педа, образованного базисными векторами, т.е. столбцами

матрицы

 

. Этот определитель может

быть

 

очень малым по

двум

А

увеличить

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

причинам. Одна состо т в том, что векторы образуют весьма малые углы

Сдруг с другом. Другая причина объясняется весьма малыми длинами

некоторых

 

з векторов.

Поскольку параметры

а

1,

а

2 ,

а

3 ,...,

а

n

имеют

 

 

 

выбраны

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

могут быть

 

произвольно. Следовательно, неудачным выбором

масштабов

 

можно

 

 

 

косоугольность

выбранной

системы

координат

 

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

 

 

 

 

 

 

 

 

 

 

Если по как м-л

о соо ражениям при построении диагностической

модели

 

 

А

 

 

параметров,

то

используются

 

размерные

значения

 

 

 

неблагоприятное влияние масшта ных факторов частично может быть

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

матрицы условной системы уравнений [8].

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Д

 

 

При обработке условной системы уравнений методом наименьших

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

которая затем доводится до минимума. В этом процессе может оказаться,

что сами уравнения недостаточно сбалансированы. Тогда каждое

уравнение системы (7) перед возведением в квадрат полезно умножить на

 

 

 

 

 

 

 

 

И

произвольный множитель Кi.

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

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

283

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

Описанный процесс нормирования уравнений условной системы проводится поочередно по столбцам и строкам матрицы. Ясно, что Свыравнивание длины столбцов нарушит равновесие длины строк. Поэтому

требуется несколько повторных подгонок строк и столбцов, пока их суммы квадратов будут достаточно выровнены.

Если содержан е задачи не требует, чтобы исходные данные были

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

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

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

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

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

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

уравнений смещений, связанных с применением одних и тех же единиц измерения, но с различными началами отсчета.

являются наиболее подходящими исходными данными.

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

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

Проведена оценка влияние погрешностей измерений на параметры И

Библиографический список

1.Грехов, Л.В и др.. Конструкция, расчет и технический сервис топливоподающих систем дизелей: Учебное пособие. / Л.В. Грехов, И.И. Габитов, А.В. Неговора. - М: Легион-Автодата, 2013. – 292 с.

2.Грехов Л.В., Иващенко, Н.А., Марков, В.А Топливная аппаратура и системы управления дизелей: Учебник для вузов. -М.:Легион-Автодата, 2005. –344 с.

284

3. Загитов, И.И. Диагностирование геометрических параметров топливных струй многодырчатых форсунок [Текст] / И.И. Загитов, А.В. Неговора // Химия в сельском хозяйстве: материалы Всероссийской научно-практической конференции для студентов и аспирантов. – 2014. – С. 164-167.

4. Неговора, А.В и др.. Диагностирование топливной аппаратуры автотракторных дизелей [Текст] /А.В. Неговора., Л.В. Грехов, И.И. Габитов // Актуальные проблемы теории и практики современного двигателестроения сб. науч. тр. междунар. науч.-техн. конф. 100-лет Вибе. - Челябинск: ЮУрГУ, 2003. - 85 с.

5.

Култышев, С. Ю. Приближенная идентификация при измерениях с

погрешностями / . Ю. Култышев, Л. М. Култышева // Вестник ПГТУ, 2010. - № 15. -

. 53-61.

 

 

6.

Култышев, С.Ю., Култышева, Л.М. Имитационное моделирование реальных

объектов

параметр ческая идентификация функциональных систем // Вестн.

Пермск. гос. тех. ун-та. Матем. и прикл. матем.,1996. - № 1. - С. 17-22.

С

нестационарных объектов // Автоматика и

7.

Клейман, Е.Г. Идент

телемехан

ка, 1999. - № 10. - С. 3-45.

 

8. Мышенков, В.И., Мышенков, Е.В. Численные методы. Часть первая: Учебное

пособие для студентов спец альности 0101.07. – М.:МГУ,2001. – 120 с.: ил.

фикация

 

 

 

бА

УДК 629.08

ИССЛЕДОВАНИЕ НЕР ВНОМЕРНОСТИ РАБОТЫ ТОПЛИВНОЙ ПП Р ТУРЫ ИЗЕЛЕЙ

В.А. Щеглов, канд. техн. наук, доц., профессор кафедры АТиСА,

БГАРФ ФГБОУ ВО Д«КГТУ», г. Калининград

V.A. Shcheglov, Ph.D., associate Prof., Рrofessor of road transport and service vehicles, «BFFSA» FSEI HPE «Kaliningrad State Technical University», Kaliningrad

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

На этой основе открываются реальные возможности совершенствования методов комплектования элементов топливной аппаратуры (ТА), способствующих повышению надежности и экономичности дизелей.

Abstract. Developed schematic diagram of the experimental setup, which allows to study the cyclic fuel supply into the possible modes of sets of

285

fuel equipment. Experimental study of kits for the detection of cyclic patterns of change serves as a function of speed strips and speed and quantify the results of an experimental study.

On this basis, the real opportunities for improving the acquisition of elements of the fuel injection equipment help to establish the reliability and energy efficiency of diesel engines.

С

дизель, топливная аппаратура, цикловая подача

Ключевые слова:

топл ва, стенд, равномерность подачи топлива.

Keywords: diesel, fuel injection equipment, injection quantity, stand, доминирующuniformity of fuel supply.

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

ми.

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

Повышение неравномерности распределения нагрузки по цилиндрам

двигателя является одной

из непосредственных причин, влияющих на

 

 

 

Д

ухудшение заградительных характеристик, снижение экономичности и

эксплуатационногобАресурса и других нежелательных явлений.

полностью

Принято

считать,

что

 

нагрузка

цилиндров

обусловливается цикловыми подачами топлива [1].

 

Учитывая большое

влияние

технического состояния основных

 

 

 

 

 

И

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

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

286

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

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

Сосуществлять соответствующую подрегулировку.

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

дизелейования.

Требован я к стенду и методика исследования, оценка точности.

учетом поставленных задач, а также проведенного анализа существующ х спытательных стендов ТА дизелей отечественного и

1.Пр водбАстенда должен о еспечивать бесступенчатое устойчивого регулирование частоты вращения вала ТНВД, что может быть достигнуто применением в приводе стенда электромотора постоянного тока.

2.Регулирование частоты вращения вала ТНВД должно производится

вшироких пределах, при нижней границе от 30 мин-1 (работа на пусковых режимах) до 2000 мин-1 по верхнееДгранице. В соответствии с ГОСТ 8670

– 69 постоянство частоты вращения приводного вала должно быть в пределах ±5 об/мин в течении 5 мин. Большие частоты по верхней

границе, которые имеют место на современных испытательных стендах

(до 4300 мин-1) нецелесообразны по следующим соображениям:

- в прямой пропорции от частоты увеличиваютсяИтребования к мощности приводного электромотора;

- диапазон частот дизелей автотракторной техники находится в пределах по верхней границе 2600 – 3000 мин-1. С учетом возможного

перспективного форсирования дизелей объемного смесеобразования до 4000 мин-1 предельная частота вращения создаваемого стенда - 2000 мин-1;

3.Мощность приводного электродвигателя (КВт) определена по

зависимости [4].

 

Nн = 11· 10-5 рВср Vц nн iп ,

(1)

где рВср – среднее давление впрыска топлива (37 МПа); Vц

– цикловая

подача на номинальном режиме (0,110 г); nн – номинальная частота вращения вала ТНВД; iп – количество плунжерных пар ТНВД, работающих под нагрузкой.

287

4. Устройство для автоматического поддержания температуры топлива на входе в насос, предусмотренное ГОСТ 8670-82 для заводских испытательных стендов, работающих в длительных циклах, иметь для исследовательского стенда нецелесообразно по следующим причинам:

- длительность наибольшего цикла экспериментальных исследований по нашему опыту не превышает двух часов;

- температура отапливаемого помещения лаборатории, в которой расположен стенд, стабильно находится в пределах 18 - 21ºС;

- емкость топл вного бака стенда – 50 л, максимальный расход

топлива на одну форсунку - 12 кг/ч и при сливе нагретого после впрыска

топлива в запасной бак указанной емкости достаточно на 4 часа работы на

топливе, температура которого находится в пределах требований ГОСТ

С

 

 

8670-82 соответствует рекомендациям [5].

5. При рассмотрен

двух возможных способов измерения цикловой

подачи топл ва – весовой или о ъемный учтено, что первый является

более предпочт тельным по следующим причинам:

- возможность пр

 

электронных весов, обеспечивающих

менения

 

точность взвеш ван я навесок топлива ±0,5 г [6], что при средних

б

5.Рабочий о ъемАмерных мензурок должен быть таким, чтобы количество топлива, отмеренное дляДкаждой секции насоса, было не менее 85 см3. Времязамера должно бытьнеменее1 мин согласно ГОСТ 8670–58.

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

обеспечивать определение числа оборотов кулачкового вала насоса за время измерения производительности топливногоИнасоса с точностью ± 1 оборот.

7.Точность установки угла начала впрыска в соответствии с ГОСТ 10578 – 96 должна быть в пределах ±30`. Стенд должен быть оборудован оптическим датчиком с точностью ±15´ для определения действительного начала подачи топлива.

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

вимпульсности реактивного момента со стороны ТНВД. Наметилась тенденция упрощения привода стенда за счет использования мощного электродвигателя с регулированием электрических параметров и обеспечением равномерности вращения за счет обратной связи [6].

288

Для исследования неравномерности распределения цикловых подач по цилиндрам дизелей в широком диапазоне скоростных и нагрузочных режимов разработан и создан испытательный стенд на основе промышленного оборудования КИ-921М.

тенд предназначен для испытания топливных насосов высокого давления (ТНВД) с цикловой подачей до 250 мм³/цикл для дизелей, в т.ч. отвечающих нормам от Euro 3 до Euro 5.

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

управлен ем, насосов

форсунок систем «Common rаil».

На оборудован

КИ-921М выполнены следующие доработки:

- электродв гатель привода ТНВД переменного тока N=2 КВт

С

 

заменен на электродв гатель постоянного тока (ЭПТ) N=8 КВт;

пеногаситель

- в блок управлен я установлен реостат сопротивлением 3 Ом;

- установлено

мерительное устройство, включающее конус-

точностью змеренбя ±0,5 г; - выполнены стойки со сменными полками (для форсунок дизелей

(КПГ) с отверстием в нижней части диаметром 1,5 мм, мерную емкость (МЕ) о ъемом 500 мл, электронные весы (ВЭ) с

ЯМЗ и КамАЗ) для крепления форсунок;

- испытываемые Афорсунки (ФР) закреплены на специальной подставке, под которой размещены КПГ, МЕ, ВЭ. Конус-пеногаситель соединен с подставкой олтами и посажен на герметик, что позволяет исключить выход из КПГ топливаТНВДв мелкодисперсном состоянии и его потерю в окружающую среду;

- для фиксации и измерения хода рейки ТНВД в конструкцию регулятора частоты вращения насоса введены следующие изменения: грузы регулятора сведены друг к другу и скреплены скобой, рейка ТНВД отсоединена от главного рычага регулятораИи ее конец соединен со штангенциркулем, закрепленном на корпусе . Это позволяет

- установлен измерительный блок, включающий аналого-цифровой преобразователь (ИФ), ПЭВМ (М и СБ) с платой L-154, фотоэлектрический счетчик частоты вращения вала ТНВД;

- с использованием обтюратора в экспериментальном стенде использован фотодатчик ФПИ-2 для получения отметки верхней мертвой точки вала ТНВД и отсчета частоты вращения вала насоса (точность - 0,3 %).

Краткая характеристика стенда:

1. Мощность приводного электродвигателя - 8 кВт.

289

2.

Число одновременно испытываемых форсунок и плунжерных пар –

8 шт.

 

3.

Длины трубопроводов высокого давления - 415 мм.

4.

Рабочее тело - дизельное топливо по ГОСТ 305-82.

5.

Рабочее тело - дизельное топливо по ГОСТ 305-82.

6.

Максимальная средняя цикловая подача - не более 300 мм3/цикл.

С

7.

Метод замера - весовой.

8.

Тип и характеристика фотометрического отметчика ВМТ ФПИ-2:

 

- фотоэлектр ческ й, с фотоэлементом вакуумного типа;

- ш р на прорези обтюратора, мм – 1,5; - погрешность змерения, % - 0,3.

9. Д апазон частот вращения вала ТНВД – 300 – 1500 мин-1.

10. Д апазон хода

ТНВД – 0…15 мм.

11. Емкость топл вного

ака – 50 л.

рейки

бА

Рис. 1. Установка штангенциркуля дляДизмерения хода рейки топливного насоса

 

И

290

С

 

 

 

 

 

 

бак

 

 

 

 

 

 

. 2. Схема стенда:

 

 

РисБ – топл вный ; ТПН – топливоподкачивающий насос; ТНВД – топливный

 

насос высокого давления; ЭПТ – электродвигатель постоянного тока;

 

 

 

А

 

 

ФР – форсунка; Ф – фильтр; КПГ – конус-пеногаситель;

 

 

МЕ – мерная емкость; ВЭ – весы электронные

 

Стенд

может

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

для

насос-форсунок

с

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

электронным управлением, насосов и форсунок

систем «Common rаil».

 

Дизельное

 

 

 

 

 

 

Принцип работы стенда.

 

топливо (по ГОСТ 305-82)

с

кинематической вязкостью 6 - 8 мм2/с и температурой 25-30ºС забирается из топливного бака Б (рис. 2) топливоподкачивающим насосом ТПН через фильтр Ф и далее подается в магистраль низкого давления ТНВД. В

зависимости от установленного и зафиксированного хода рейки некоторый объем топлива подается плунжерной парой к форсунке ФР. Топливо из форсунки поступает в мензурку.

И

Время набора навески топлива на выбранном режиме замеряется секундомером. Прозрачный конус-пеногаситель позволяет визуально контролировать качество распыла топлива форсункой. Масса навески топлива в мерной емкости МЕ измеряется электронными весами ВЭ. Частота вращения вала ТНВД устанавливается при помощи реостата, включенного в цепь обмотки возбуждения электродвигателя ЭПТ.

Методика проведения эксперимента. После установки комплекта топливной аппаратуры необходимо согласно требованиям ГОСТ 10578-96 «Насосы топливные дизелей. Общие технические условия» проверить

291

гидроплотность плунжерных пар ТНВД, нагнетательных клапанов ТНВД, игл форсунок.

После включения стенда устанавливается необходимый скоростной режим и положение рейки ТНВД и добиться устойчивой работы. При работе стенда не должно быть подтеканий ТНВД, форсунок и

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

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

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

уровень [7].

Экспер ментальные сследования включают в себя два этапа. На

этапе доводки

отладки испытательного стенда

были выполнены

следующ е

 

:

 

 

 

 

 

равномерности

 

 

 

 

 

- проведена проверка гидравлической плотности испытуемых

форсунок по зап рающему конусу в соответствие с ГОСТ 10578-96 и

плунжерныхи клапанных пар (результаты проверкипредставлены втабл.1);

- проведенаработыпроверка форсунок на начало подачи топлива прибором

КИ-9917 и приведение всех испытуемых форсунок в соответствие с

требованиями нормативной заводской документации - 20 МПа;

- проведена проверка и регулировка испытуемых плунжерных пар

(ГОСТ 25708-83) на производительностьАи приведение ее в соответствие с

требованиями нормативной документации [8];

 

 

- проведена проверка пропускной способности форсунок (ГОСТ

10578-96) прокачкой топлива через форсунки (после проверки всех

элементов ТА на гидроплотность) соответствующей секцией топливного

 

 

 

Д

насоса при номинальных частоте вращения вала насоса и подаче топлива.

 

 

 

 

 

 

Таблица 1

 

Способы проверки элементов топливной аппаратуры

Проверяе-

 

Способ проверки и

 

При-

Контролируе-

Требования

 

мый элемент

исходные регулировки

 

бор

мый параметр

ГОСТ

 

1

 

2

 

3

4

5

 

 

 

 

 

 

И

 

Гидроплот-

Ручная прокачка,

 

КИ-

Падение

Время – 15 с,

 

ность иглы

давление начала впрыска

 

9917

давления от

соответствует

 

 

– 20 МПа

 

 

22 МПа до 21

ГОСТ 10579-88

 

 

 

 

 

 

МПа

 

 

Гидроплот-

Ручная прокачка,

 

КИ-

Отсутствие

Соответствует

 

ность иглы

давление начала впрыска

 

9917

потения на

ГОСТ 10579-88

 

по запираю-

– 20 МПа, создаваемое

 

 

носике

 

 

щему конусу

давление – 19 МПа

 

 

распылителя

 

 

292

Окончание таблицы 1

 

1

2

3

4

5

 

Качество

Ручная способом 60 – 80

КИ-

Качество

Туманообразное,

 

распыли-

впрысков в минуту

9917

распыленного

без сплошных

 

вания

 

 

топлива

струек?

 

форсункой

 

 

 

Соответствует

 

 

 

 

 

ГОСТ 9928 - 71

 

Герметич-

Опрессовка воздухом

 

Пропуск воз-

Время – 15 с,

 

ность уплот-

давлением не менее 0,4

 

духа в тече-

соответствует

 

нений низко-

МПа

 

нии 10 с не

ГОСТ 10579-88

 

го давлен я

 

 

допускается

 

 

Гидроплот-

Подача топл ва через

КИ-

Падение дав-

Время – 15 с,

 

ность плун-

штуцер насоса при

9917

ления от 20 до

соответствует

 

жерной пары

заторможенном приводе

 

19,5 МПа

ГОСТ 10578-96

С

 

 

 

 

 

плунжерной пары

 

 

 

 

-

Подача топл ва через

КИ-

Падение

Время – 15 с,

 

ность

л н ю подвода топлива

9917

давления от

соответствует

 

запорного

от топл воподкачи-

 

20 до 14 МПа

ГОСТ 10578-96

 

клапана

вающего насоса

 

 

 

 

Гидроплот

 

 

 

 

Подвиж-

Ручным спосо ом 30 – 40

КИ-

Число

Звук скрипа,

 

ность иглы

впрысков в м нуту,

9917

впрысков в

соответствует

 

 

давлен е начала впрыска

 

минуту

ГОСТ 25708-83

 

 

– 20 МПа

 

 

 

В результате ра оты стенда

ыли получены навески топлива G (г), за

время t (мин), снятые во всем поле возможных режимов работы ТА на

 

 

Д

 

строго фиксированных частоте n (350 – 1050 мин-1, с интервалом 100 мин-

1) и ходе рейкибАh (2 – 10 мм, с интервалом 1 мм).

 

 

Цикловая подачатопливаgц определялась по зависимости:

 

 

 

 

G

 

 

 

 

gц= t n .

 

 

(2)

Равномерность

подачи

топлива

оценивается

степенью

неравномерности:

δн = 2(gцmax - gцmin)/(gцmax - gцmin) 100 %,

 

(3)

 

 

где gцmax – весовое или объемное количество топлива, подаваемое секцией

с максимальной производительностью; gцmin – весовое или объемное

количество

топлива,

подаваемое

секциейИс минимальной

производительностью.

Неравномерность распределения цикловых подач по цилиндрам дизеля двух комплектов ТА представлены на рис. 3 и 4.

293

С

 

цилиндрам

 

Р с. 3. Неравномерность распределения цикловых подач

бА

по

дизеля комплекта 1

 

Д

 

И

Рис. 4. Неравномерность распределения цикловых подач

по цилиндрам дизеля комплекта 2

Из данных рисунков видно, что неравномерность распределения цикловых подач по цилиндрам дизеля для двух комплектов составляет: на режиме холостого хода 46 и 37 % соответственно, на номинальном режиме 8 и 6 % соответственно.

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

294

экспериментов, проведенных в идентичных условиях. Сначала вычисляется среднее арифметическое значение показателя по уравнению:

 

 

 

 

 

 

 

 

1

 

k

 

 

 

 

 

 

 

X

 

* x

 

 

 

(4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i ,

 

 

 

 

 

 

 

 

 

 

 

k i 1

 

 

 

 

где X - среднее арифметическое значение показателя; xi – значение

С

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получаемой величины; k – количество экспериментов (замеров).

 

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

дисперс , которая рассч тывается по уравнению:

 

 

 

 

 

 

2

 

1

 

 

 

k

 

2

 

 

вычисляется

 

 

 

 

i k 1

.

 

 

(6)

 

i 1

 

 

 

 

 

 

 

 

 

 

* (xi x) .

 

(5)

 

 

 

 

 

 

 

 

k 1 i 1

 

 

 

 

редняя

квадрат чная

 

 

погрешность

 

отдельного измерения

по уравнен ю:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

x *

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Оценкой

коэфф ц ента

 

 

вариации

 

служит

статистический

коэффиц ент вар ац , определяемый по формуле (%):

 

 

 

 

V 100*

.

 

 

 

(7)

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

Так как ни одна случайная погрешность (xi – X), с вероятностью

практически равной единице, не может выйти за пределы ±3

[7], то

 

 

 

 

 

 

 

 

Д

 

после оценки среднего квадратичного отклонения для каждой выборки

измеряемой

бАвеличины наблюдений, была произведена проверка на

наличие в выборках наблюдений грубых погрешностей. Наблюдения,

содержащие грубые погрешности были исключены из выборок

наблюдений, после чего Х и σ были вычислены заново.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

И

Проверка гипотезы о том, что результаты наблюдений принадлежат нормальному распределению, производилась с использованием критерия х2 - Пирсона. Доверительную вероятность Р=0,95 следует понимать так, если сделать большое число выборок, то в среднем до 95 % выборок получатся такие доверительные интервалы которые будут «накрывать» оцениваемое математическое ожидание [8].

Предельная абсолютная ошибка, характеризующая точность оценки математического ожидания с вероятностью γ, определяется по формуле:

;k

t ;k

 

 

 

 

 

,

(8)

 

 

 

 

 

 

 

n

 

295

где tγ;к – параметр распределения Стьюдента, как функция заданной вероятности γ и числа степеней свободы к; к – число степеней свободы (к = n-1).

Если в (9) выразим абсолютную ошибку δγ;к через относительную, то получим:

С

 

 

 

 

 

 

 

 

;k t ;k

 

 

 

,

 

(9)

 

 

X

 

 

 

 

 

 

n

 

 

где

γ;к – предельная относительная ошибка, гарантируемая

с

вероятностью γ при объеме выборки n=k+1.

 

 

 

 

 

 

 

 

При проведен

экспериментов выяснилось, что от цикла к циклу

наличием

 

 

 

 

 

 

 

имеет место некоторая нестабильность результатов, образующаяся в

результате нестандартности каждого отдельного цикла. Это объясняется

 

сложных ф з ческих явлений, обусловленных волновыми

явлениями в элементах ТА.

 

 

 

 

 

 

 

 

б

 

 

 

 

 

 

 

 

Для оценки степени неста ильности циклов и определения

необход мого кол чества реализаций циклов подачи топлива форсункой

на основе экспер ментальных данных, полученных в ходе исследований,

выполнены расчеты, основанные на методике, изложенной в работе [9].

 

 

 

А

 

 

 

Были получены по 14 реализаций статистической выборки цикловой

подачи топлива, снятых на номинальном режиме, режиме долевых

нагрузок и на холостом ходу в виде массива:

 

 

 

 

 

 

 

{ х1, х1, х1, …, хn},

 

 

 

(10)

 

Алгоритм расчетов с использованием Mathcad:

 

 

 

 

Д

 

 

 

- исключение из массива повторяющихся членов (корневой массив);

 

- сортировка массива, когда смежные члены комплектуются попарно

минимальными и

максимальными членами

выборки с пошаговым

уменьшением разницы между минимумом и максимумом, разбивка корневого массива на ряд массивов с пошаговым увеличением количества членов от 2 доi;

- определение статистического среднего m*x каждого массива;

- определение максимально возможной ошибки между смежными массивами i = (m*x i - m*x i+1)*100/ m*x I ;

- определение текущих значений верхней границы доверительного

интервала m1i

= + i ;

И

- определение текущих значений нижней границы доверительного

интервала m 2i = - ∆i ;

- определение доверительного интервала по условию |m1i| + |m2i| ≤ зад. При анализе полученных результатов, а также с учетом точности применяемой в ходе экспериментов аппаратуры заданная точность ∆зад = 0,03. При достижении этой точности определяется необходимое

количество реализаций.

296

Пример графического представления результатов расчета необходимого количества реализаций по наиболее важному параметру – цикловой подаче топлива, соответствующих номинальному, частичному режимам, режиму холостого хода приведены на рис. 5, 6, 7. Количество реализаций по соответствующим режимам работы ТА приведено в табл. 2.

Си

Рис. 5. ИнтервалыбАточности измерения цикловой подачи gц на номинальной частоте вращения вала ТНВД n=1050 мин-1 и ходе рейки h=10 мм

Д И

Рис. 6. Интервалы точности измерения цикловой подачи gц на средней частоте вращения вала ТНВД n=650 мин-1 и ходе рейки h=6 мм

297

С

 

 

 

 

 

 

реализаций

 

 

 

 

 

Рис. 7. Интервалы точности измерения цикловой подачи gц на режиме холостого

 

хода при частоте вращения вала ТНВД n=350 мин-1 и ходе рейки h=2 мм

 

 

 

Кол чество

 

для проведения исследований

Таблица 2

 

 

 

 

 

 

Диапазоны реж мов ра оты

 

 

Параметры

 

 

Частота,

 

Ход рейки,

Количество реализаций

Необходимое коли-

 

n, мин-1

 

h, мм

для параметра, gц, г/цикл

чество реализаций

 

1050

 

10

 

 

6

 

6

 

650

 

6

 

 

10

 

10

 

350

 

3

 

 

13

 

13

 

Выводы:

 

 

Д

 

1. РазработанабАпринципиальная схема и изготовлен стенд для

 

экспериментального исследования функционирования ТА дизелей в поле

 

возможных

режимов

работы.

Экспериментально

исследованы

закономерности изменения цикловой подачи топлива в зависимости от И

частоты вращения вала насоса и хода рейки. Получены экспериментальные данные по неравномерности распределения цикловых подач по цилиндрам дизелей для 2-х комплектов ТА.

2. Проведена оценка точности статистически среднего и выполнен расчет необходимого количества реализаций.

Библиографический список

1.Габитов, И.И., Неговора, А.В. Топливная аппаратура автотракторных дизелей/ И.И. Габитов, А.В. Неговора. – Уфа: Изд-во БашГАУ, 2004. –216 с.

2.Файнлейб, Б.Н. Топливная аппаратура автотракторных дизелей: Справочник.- 2-е изд., перераб. и доп. – Л.: Машиностроение, Ленингр отд., 1990. – 352 с.: ил.

3.Черноиванов, В.И. РТМ 10.457100 0001 - 01. Насосы топливные высокого давления автотракторных дизелей. Методы испытаний и регулирование Текст. / В.И. Черноиванов, А.Э. Северный, В.М. Корнеев [и др.]. - М.: ГОСНИТИ, 2001. -38 с.

298

4. Егошин, Е.А., Маков, И.П. Некоторые расчетные зависимости и статистические данные по подбору топливной аппаратуры: Учебное пособие. - Набережные Челны: КамПИ,1995. -104с.

5. Грехов, Л.В., Иващенко, Н.А., Марков, В.А. Топливная аппаратура и системы управления дизелей. Учебник для ВУЗов. – М.: Легион-Автодата, 2004. -344 с.

6. Корнеев, В.М. Особенности конструктивного исполнения стендов для испытания и регулирования топливных насосов высокого давления Текст. / B.М.

СКорнеев, М. Ю. Устинов // Международный научный журнал. 2010. - №3. -C. 77-81. 7. Шаламов, А.Н. Обработка результатов и оценка точности измерений при

многократных наблюден ях: учеб. пособие / А.Н. Шаламов, Б.А. Кудряшов, Т.М. Раковщик. – М.: МАДИ, 2016. -164 с.

8. Ж гадло, А.П., Иванов, А.Л., Саенко, М.М. Влияние различных факторов на испытанийвеличину равномерность подачи топлива в цилиндры дизеля. Вестник СибАДИ, выпуск 4 (32), 2013. – . 29 – 35.

9. Коньков, А.Ю., Лашко, В.А., Яранцев, М.В. Способ измерения цикловой подачи топл ва в ц л ндрах тепловозного дизеля в условиях эксплуатационных

. ВестнбАк ТОГУ №1 (24), 2012. –С. 119 – 124.

Д И

299

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