CC BY
23
Вычислительные технологии
Том 20, № 6, 2015
Моделирование процесса оптимального определения параметров состояния двигателя внутреннего сгорания измерительной экспертной системой
И. П. Добролюбов, О.Ф. Савченко*, В. В. Альт, С.Н. ОльшЕвский, Д. Н. Клименко
Сибирский физико-технический институт аграрных проблем, Новосибирск, Россия *Контактный e-mail: [email protected]
На основе анализа известных методов статистической обработки сигналов исследовано и обосновано применение оптимального метода обнаружения структурных параметров двигателя внутреннего сгорания, скрытых в высоком уровне других составляющих рабочих процессов. Рассмотрено использование измерительной экспертной системы для определения информативных параметров сигналов, отражающих эти структурные параметры при экспертизе технического состояния двигателей. Выполненное компьютерное моделирование оптимального обнаружения и измерения в условиях помех информативных параметров — амплитуды углового ускорения коленчатого вала и угла опережения подачи топлива — позволяют сформировать требования к наиболее эффективной структурной схеме измерительной экспертной системы.
Ключевые слова: двигатель внутреннего сгорания, ресурсные структурные параметры, сигнал, обнаружение, измерение, функция правдоподобия, пороговый уровень, оптимальное правило решения.
Введение
При экспертизе технического состояния двигателей внутреннего сгорания (ДВС) с помощью измерительной экспертной системы двигателей (ИЭСД) [1, 2] измеряются рабочие процессы (давления в камере сгорания, топливопроводе; линейные и угловые перемещения и др.), находятся диагностические признаки и по ним оценивается состояние ДВС. Предварительная обработка сигналов, поступающих с датчиков рабочих процессов, согласование с датчиками, усиление, фильтрация, дискретизация по времени, квантование по уровню осуществляются в ИЭСД устройством сопряжения с объектом. Затем сигналы обрабатываются процессорным устройством.
С целью принятия решений по оценке состояния ДВС обработка в ИЭСД проводится с применением методов анализа сигналов, полученных в результате проводимых в СибФТИ многолетних исследований. Среди них — статистические методы и частотный анализ рабочих процессов ДВС [3-7], выбор информативных признаков и параметров состояния [8], методы определения структурных ресурсных параметров ДВС и обобщенных показателей состояния [9], оценка погрешностей измерений, оптимизация
© ИВТ СО РАН, 2015
метрологического обеспечения ИЭСД [10-12], применение компьютерной настраиваемой модели ДВС [13, 14].
Проведенные исследования позволили обосновать требования к функциональным возможностям и структуре ИЭСД, разработать ее информационное, алгоритмическое и программное обеспечение и создать информационную технологию диагностирования ДВС [15-17]. Однако многообразие решаемых задач при расширении сферы применения ИЭС для оценки технического состояния отдельных агрегатов и систем двигателя, осуществлении экспертизы технического состояния для обеспечения работоспособности двигателя в процессе его эксплуатации требует совершенствования технического и программного обеспечения, решения оптимизационных задач по структуре, параметрам, методам анализа сигналов.
Как известно, работа ДВС характеризуется случайностью процессов подачи топлива и его сгорания от цикла к циклу, а также процессов трения в сопряжениях ДВС, вызывающих случайные отклонения сигналов, например угловых скорости и ускорения коленчатого вала u(t) и e(t) от их средних значений. Это приводит к необходимости рассматривать измеряемые процессы как аддитивную смесь полезного информативного (диагностического) случайного процесса и детерминированной составляющей, вызванной неуравновешенными составляющими крутящего момента. Кроме того, возникает необходимость оценить структурные параметры двигателя, диагностические признаки которых слабо проявляются на фоне других детерминированных и случайных процессов. Например, превышающие допустимый уровень такие ресурсные структурные параметры, как износы, жесткость рабочего процесса, люфты и др. (в значительной степени нелинейные элементы), характеризуются диагностическими сигналами, значение которых существенно ниже уровня измеряемых аддитивных и мультипликативных сигналов. Обнаружить и измерить такие признаки, даже при предварительном вычитании из результирующего сигнала зависимостей, известных для данного двигателя и отражающих другие структурные параметры, бывает достаточно сложно.
Цель настоящей работы — повышение эффективности идентификации технического состояния ДВС в эксплуатационных условиях с применением ИЭСД путем исследования и обоснования оптимального метода обнаружения структурных параметров ДВС, скрытых в высоком уровне других составляющих рабочих процессов с использованием известных вероятностных методов анализа сигналов, статистических методов оптимальной обработки сигналов на фоне помех [18-20].
1. Постановка задачи
Статистическая задача измерения параметров объекта экспертизы состоит в том, чтобы на основании измеренного значения вектор-столбца x(t) = {x\(t), x2(t),... ,Хм(¿)}, найденного при усреднении М реализаций (повторности) x(t) и априорной плотности распределения вероятностей (ПРВ) f (а) измеряемого вектора а = [а\, а\,..., ага} информативных диагностических признаков, получить оценку а*.
Сигнал, использующийся для обнаружения и измерения указанных признаков, описывается уравнением
x(t) = Auc (t, а, ß) + un (t, 7) ,
где А — дискретный случайный параметр, принимающий два значения: А0 = 0 и А\ = 1, что соответствует условиям отсутствия и наличия полезного сигнала uc (t, а, ß), отра-
жающего диагностируемые признаки; 0 = {^1, ... ,0т} — неинформативные параметры; ип (Ь, ^) — помеха, вызванная разбросом параметров рабочих процессов двигателя; ■у = {71,/у2,..., 11} — параметры помехи [19].
При дискретизации каждая реализация Хг (¿) [г = 1, М) состоит из Ь отсчетов. Общее число отсчетов (дискрет) т = ЬМ. Решение при обнаружении и измерении сигнала принимается по т-мерному вектору-столбцу х(Ь).
Оценка а* должна быть оптимальной, т. е. удовлетворять критерию минимума среднего риска (для квадратичной функции риска это эквивалентно требованию минимума средней квадратичной погрешности измерений [18, 20]):
+те +те
+те +те
г = М {г}
г(а*, а)/(а*, а) ¿а*¿а. = ¿х г(а*, а)/(х) / (а | х) ¿а
—те —те
— те —те
+те
г (а* | х) / (х) ¿х ^ ш1п,
+те
где г (а* | х) = / г (а*, а) / (а* | х) ¿а — условный средний риск.
—те
Для каждого измеренного значения х(£) дается вполне определенная оценка а* = а(ис). В выражении (1) ПРВ $ (х) процесса х(£) является заданной функцией и не может принимать отрицательных значений. Поэтому минимум подынтегрального выражения (1) достигается при г (а* | х) ^ шт.
Апостериорная ПРВ / (а* | х) = кх$ (а) / (х | а), где кх = 1// (х) — коэффициент,
+те
не зависящий от а. Из условия нормировки J / (а | х) йо. = 1 получаем
+те
1
/ («) / (х | ¿а
Кривая апостериорной ПРВ f (а | х) уже, чем априорной f (а). Условное отношение правдоподобия
I(х | а) = f (х | а) //га (х), где /п (х) — ПРВ помехи. Тогда имеем f (а | х) = ки£ (а) I (х | а), где Г —те I—1
Л (х)
к*.
/ (х)
/ (а) I (х | а) ¿а.
+те
1п/ (а | х) = 1п [/ (а)] + 1п [I (х | а)] + еопя^
При квадратичной функции риска для векторного и скалярного параметров условный средний риск:
г (а* | х)= [ (а* — а)Т \\В|| (а* - а) / (а | х) ¿а,
г (а* | х)= (а* — а)2! (а | х) ¿а.
к
X
Минимум (2) при симметричной неособенной матрице веса ||5||, характеризующей важность того или иного признака а, обеспечивается для векторного и скалярного параметров при следующих условиях:
те
a*opt = М (а | x} = af (а | x) da;
a*opt = М (а | x} = af (а | x) da f (а | x) da = af (a | x) da, (3)
где М {... } — математическое ожидание. При этом условный средний риск и средний риск погрешности измерения для скалярного параметра равны:
+те
г (а* | Х)=У (а - м {« I х})2 / {а I х} ^а = в {а I х} ;
—те
+те
г = ! И {а I х} f (х) вх,, (4)
—те
где Б {... } — дисперсия [18-20].
Зависимость г в (4) характеризует минимальное значение среднего квадрата погрешности измерения, которое обеспечивается условием для г (а* | х) в (4). Выполнение этого условия для любого значения х свидетельствует о том, что оценка параметра а* не имеет систематической погрешности: М {а* — а} = М {(а* — а) | х} = 0.
Допустимо полагать, что зависимость f (а | х) для процессов двигателя симметричная и унимодальная. Тогда оценки (3) и (4) совпадают с оценкой по критерию максимума апостериорной ПРВ.
Итак, ставится задача моделирования процесса оптимальной оценки информативных параметров, отражающих рабочий процесс ДВС, а именно амплитуды углового ускорения коленчатого вала, измеренной в угловом интервале работы данного цилиндра в цикле, сигнала, отражающего угол опережения подачи топлива, и совместно измеряемых параметров ДВС. Данную задачу необходимо рассмотреть также при наличии априорных данных.
2. Определение оптимальной оценки амплитуды углового ускорения коленчатого вала ДВС по цилиндрам
Одним из информативных параметров, отражающих рабочий процесс каждого цилиндра ДВС, является амплитуда углового ускорения коленчатого вала ис, измеренная в угловом интервале работы данного цилиндра в цикле. Априорную ПРВ этого скалярного параметра а = ис полагаем имеющей равномерное распределение, а помехи — нормальное (рис. 1, а):
f (ис) =-1-, иС1< ис < < иС2;
ис2 — ис1
f (ис) = 0, ис< иС1, ис> ис2.
Условная ПРВ одномерного сигнала f (х | ис) в функции неизвестного значения ис (рис. 1, а):
f (х | uc) = fn (х - uc) = (l/^2nan) exp [-(х - uc)2/2(гп2] ,
где ап — среднее квадратическое значение помехи.
На рис. 1, а представлена также кривая апостериорной ПРВ f (uc | х), полученная в результате перемножения функций f (uc) и f (х | ис) с учетом нормирующего множителя кх. Площадь под этой кривой равна единице. Функция f (uc | х) учитывает как результат измерения х, так и априорные данные о возможных значениях измеряемой величины сигнала ис и помехи. Существенное влияние на апостериорную ПРВ f (udx) оказывает уровень помех (рис. 1, б, в).
При экспертизе многоцилиндровых двигателей наблюдается значительное перекрытие измеренного углового ускорения контролируемого цилиндра угловыми ускорениями соседних с ним цилиндров, а также при идентификации нелинейностей: зазоров, сухого трения, зоны нечувствительности. Уровень этих сигналов существенно ниже уровня сигналов от работающих цилиндров. Измерения проводятся в условиях очень сильной помехи, при этом центр тяжести кривой апостериорной ПРВ и оптимальная оценка равны (uc\ + uc2)/2 независимо от измеренного значения х (рис. 1, в). Наоборот, для случая слабой гауссовой помехи (например, при экспертизе четырехцилиндровых двигателей) оптимальная оценка равна измеренному значению, так как х соответствует центру тяжести на кривой апостериорной ПРВ (рис. 1, б).
Для симметричных выпуклых кривых ПРВ центр тяжести всегда совпадает с максимумом кривой. Поэтому вместо оценки aopt = ucopt в качестве оптимальной оценки можно использовать абсциссу максимума кривой апостериорной ПРВ f (uc | ж), т. е. наи-
Рис. 1. Плотности распределения вероятностей для амплитуды ускорения коленчатого вала ис: а — формирование апостериорной ПРВ: кривая 1 — функция /(ис), 2 — /(х\ас), 3 — /(-ис|ж); б — ПРВ для слабой помехи (ап < ис2 — ис1): 1 — /(х^с) & /(-ис|ж), 2 — /(ис); в — ПРВ для сильной помехи (ап > ис2 — ис1): 1 — f (х^с) & (ис), 2 — f (ис)
X
Рис. 2. Интервал погрешностей измерения параметров ДВС: а — функция /(а|х), Ь = 2к<га при измерении одного параметра (одномерное измерение); б — функция /(а|х) = /(а1|х)(а2|х), 1 — 2ка2, 2 — 2ка1 для совместно измеряемых двух параметров (двумерное измерение)
более вероятную оценку. Для слабых гауссовых помех, так как ап не зависит от х, показатели качества измерения а* = и**:
г (а* | х) = И {а | х] = ап2, г = ап2.
Для сильных помех дисперсия апостериорной ПРВ совпадает с дисперсией априорной ПРВ, а ее значение также не зависит от х:
г = г | ж) = О {а | х] = И {а]
а 2
2
1 / /* 1 \ (^ ^ )2
И {а] = а2-йа — \ а-йа I =-,
] «2 — «1 \J «2 — 12
«1 \"1 /
где «2 — «1 = иС2 — иС1.
Оптимальная средняя квадратичная погрешность измерения при изменении помехи
«2 — «1 иС2 — исЛ
от очень сильной до слабой изменяется в -=-раз.
апл/12 апл/12
Если априорная ПРВ параметра а = ис неизвестна, то из-за множества влияющих параметров одномерную апостериорную ПРВ можно считать нормальной:
/ (а | х) = 1/^/2паа ехр [—(а — а*)2/2иа2] .
Дисперсия погрешностей точечной оценки скалярного параметра а = ис составляет
в2 -1 " 1п I (х | а*) .
с1а2
Вероятность попадания истинного значения а = ис при одномерном измерении (п = 1) в интервал I = 2каа (рис. 2, а) равна
к
Ро = I е-1'/2А1. \J2-k
о
Например, при Р0 = 0.9 значение к = 1/2оа = 1.64, I = 2каа = 3.28аа
3. Определение оптимальной оценки угла опережения подачи топлива
В качестве измеряемого параметра, отражающего угол опережения подачи топлива, можно рассматривать время задержки сигнала ¿з относительно верхней мертвой точки при известной частоте вращения коленчатого вала. Согласно (3) при измерении оптимальной оценкой является центр тяжести функции апостериорной ПРВ / {а* | х(£)} = / {¿* | х(^)} . Принимая в качестве оптимальной наиболее вероятную оценку = ¿* ^ определяем максимум функции / {а* | х(1)} = {1* | х(£)}:
в
— /{1з I х(1)} = 0, 1з = 1з*. &+
Так как на интервале изменения ¿з1 < Ьз < ¿з2 априорную ПРВ можно считать неизменной, достаточным является выполнение одного из двух условий:
в
—I {х(г) I г3} = 0, ¿з = ^*, в
— 1п I {х(1) I г3} = 0, г3 = г3*. и+
Поскольку значения I {х(^ | 1з} определяются для различных ожидаемых значений ¿з последовательно во времени (по множеству циклов ДВС), о них можно судить по огибающей напряжения отдельных импульсов на выходе оптимального фильтра. Такое же последовательное поступление данных во времени происходит и после некогерентного суммирования. Поэтому поиск максимума определяемой функции практически сводится к нахождению момента наступления максимума наибольшего импульса напряжения на выходе схемы некогерентного суммирования (накопления сигналов по множеству циклов ДВС [18-20]).
Влияние помехи (от слабой до сильной) на погрешность измерения параметра ¿з и вероятность попадания в интервал апостериорной ПРВ аналогичны рассмотренным в предыдущем пункте. Описанные ранее зависимости справедливы при измерении других скалярных структурных параметров ДВС (например, при обнаружении и измерении уровня зазоров в зоне перекладки поршня в цилиндре, уровня сухого трения при повышенной жесткости работы двигателя, зоны нечувствительности в сопряжениях ДВС и др.).
4. Определение оптимальной оценки совместно измеряемых параметров ДВС
При измерении векторного параметра a
a\ «2
(например, амплитуды индикатор-
ном составляющей углового ускорения отдельного цилиндра и смещения максимума этой составляющей относительно верхней мертвой точки по углу поворота коленчатого вала ДВС, характеризующего временное запаздывание или угол опережения впрыскивания топлива) сигнала uc(t) полагаем составляющие погрешности измерения а\ и а2 независимыми:
1/(7? 0 0 1/7??
Ñ, —
Апостериорная плотность вероятности
/ 1 х) = / 1 х) \ х)
1
2жа 1&2
ехр
!
(аг — а\*)2 (а2 — а2*У
о + о
о\
Оч
Оценки а* и а* соответствуют решениям уравнений правдоподобия d 1п 1/йа.1 = 0, А 1п1/ё,а2 = 0 при а\ = а* и а2 = о*2. Дисперсии погрешностей измерения а\ и о| определяются из условий
1
а
1,2
¿21п I
¿а2
а±=а*, а2 =а*
Полуоси эллипса погрешностей (рис. 2, б) при Р0 = 0.9 и п = 2 имеют размеры:
ках = 2.12аъ ка2 = 2.12а2.
2
5. Оптимальное обнаружение сигналов двигателя и измерение их параметров при наличии априорных данных
Априорные данные позволяют значительно уточнить оценку измерения параметров а. Допустимо полагать априорную ПРВ вектора параметров а нормальной:
/ (а) = (2п)-п/2\\Со\\1/2 ехр
где а** — его априорная оценка; \\С0\\ = измерения а.
Апостериорная ПРВ [18-20]:
f (а \ х) = (2п)-п/2\\Ср\\1/2 ехр
— \(а — а0о)Т ||Оо\ (а — а0)
\\С0\\Т — симметричная матрица точности
—2 (а—а*р)т \\сг\\ (а—
где а* = ||СР\\-1 (\\С0\\ а* + \|Сж\ а*) — оценка вектора параметров при измерении х; 11 ^^^ \ = \ \ М + \ \ \ \ — апостериорная матрица точности измерения х; | \ \ — матрица точности измерения х.
Оценка а*х является апостериорным математическим ожиданием вектора параметров а и соответствует максимуму апостериорной ПРВ / (а \ х), которая оптимальна по критерию минимума среднего риска как при квадратичной, так и при простой функции риска. Результирующая корреляционная матрица погрешностей \\СР\\ 1 получается путем обращения суммарной матрицы точности \\СР\\. С учетом того, что
\\сР\\-1 \\С0\\ + ||ср\\-1 \\a\l = \\ср\\-1 (МС0М + \\СХ
1,
оценка вектора параметров примет вид
а; = а*+\\ Ср\\ -1\\ С, \ \ (а: — о0).
В качестве примера рассмотрим измерение, когда вектор параметров сводится к скалярному параметру а (см. рис. 1). Априорная плотность вероятности / (а) имеет нормальное распределение с дисперсией аапостериорная плотность вероятности / (а \ х)
«2 2 имеет нормальное распределение с дисперсией а", а результирующая дисперсия — а°.
При этом ||Со|| = 1/сг°, ЦС^Ц = , Н^рН = 1/^2 = 1/ао + 1/аж. Оценка параметра а*
находится по функции
2 2 2 * ах * , *
~ .1 и * * , Р / * *\
= —--ап +—о-2 = ао +—2 (а% — ).
Р 2 I 2 и ' 2 I 2 ж и ^ 2
Результирующая апостериорная дисперсия не превышает как априорную а°, так и дисперсию измерения а'°: о^ < а°, о^ < Оценка а* представляет собой весовую сумму оценок ^ и а* с коэффициентами, составляющими в сумме единицу, поэтому а* располагается между этими оценками. В крайнем случае (а° ^ а"2) получаем значение «о ~ а1, т.е. априорная информация не улучшает результат измерения. В другом крайнем случае (а2 ^ а°) имеем значение а* ~ ^, т.е. априорная информация намного надежнее результата измерения. Весовой коэффициент невязки 0 < а'22/о2. < 1. Он
близок к единице при а"2 ^ а° и близок к нулю при о^ ^ а°
р / ^ х
и и близок к пулю при их ^^ ид .
Для определения результирующей оценки а* необязательно находить текущую оценку а* на каждом этапе измерения. Более удобно оценивать малую невязку а* — или величину, ей пропорциональную. При этом в качестве опорного значения можно в измерительном устройстве использовать номинальное или другое заданное значение (например, соответствующее определенному классу состояния ДВС). В качестве примера на рис. 3 представлено моделирование процесса оптимального обнаружения и измерения нелинейностей типа "люфт" в зоне перекладки поршня одного из цилиндров в режиме разгона ДВС на фоне помех (угловых ускорений от работы других цилиндров, разброса этих значений из-за случайности процессов сгорания топлива, трения в сопряжениях двигателя и др.) с применением статистических методов.
Заключение
В работе на основе методов статистической обработки сигналов исследовано и обосновано применение оптимального метода обнаружения структурных параметров ДВС, скрытых в высоком уровне других составляющих рабочих процессов, измерения информативных параметров сигналов, отражающих эти структурные параметры при экспертизе технического состояния двигателей измерительной экспертной системой.
Показано, что оценка измеряемых информативных параметров сигналов будет оптимальной, если выполняется требование минимума средней квадратичной погрешности измерений этих параметров, которое эквивалентно критерию минимума среднего риска.
В качестве примера рассмотрено оптимальное обнаружение и измерение в условиях слабой и сильной помех одного из информативных параметров, отражающих рабочий процесс каждого цилиндра ДВС, а именно амплитуды углового ускорения коленчатого вала, измеренной в угловом интервале работы данного цилиндра в цикле. Установлено, что для случая слабой гауссовой помехи (например, при экспертизе четырехцилиндровых двигателей) оптимальная оценка равна измеренному значению, так как в этом случае измеренный процесс соответствует центру тяжести кривой апостериорной плотности распределения вероятностей. Поэтому целесообразны применение накопителей сигналов и выделение сигналов на фоне помех с помощью фильтрации. Аналогичные исследования и выводы справедливы для сигнала, отражающего угол опережения подачи топлива, а также при совместной обработке этого сигнала с предыдущим сигналом.
150 Ф,град
150 ф>град
150 ф,град
150 Ф,град
Рис. 3. Моделирование статистическими методами процесса оптимального обнаружения и измерения нелинейностей типа "люфт" на фоне помех в зоне перекладки поршня (0-25° угла поворота коленчатого вала) одного из цилиндров в режиме разгона ДВС: а — структурная схема процесса моделирования; б — входное, максимальное, скользящее среднее и среднеквадратическое значения суммарного углового ускорения цилиндров (Results); в — импульс индикаторной составляющей углового ускорения одного цилиндра при появлении нелинейности "люфт" (Scope); г — импульс индикаторной составляющей суммарного углового ускорения цилиндров без нелинейности (Scope 2); д— импульсы индикаторной составляющей суммарного углового ускорения цилиндров с нелинейностью "люфт" (Scope 4); е — разностный сигнал между ускорением без нелинейностей и ускорением с наличием нелинейностей (Scope 3)
I
Ь
1
о to
CD
S
О -fc! CD О
0
go §
1
1
Й О
о
0
CD fcl CD
¡4
ы "а
go
1 £
о to
О
о
0 ч
1
¡4
со
При этом необходимо, по возможности, в качестве пороговых (опорных) значений использовать полученные ранее априорные статистические характеристики и параметры, отражающие различные классы состояний ДВС (номинальный, допустимый, предава-рийный, аварийный).
Компьютерное моделирование процесса оптимального обнаружения структурных параметров ДВС и их измерения на фоне помех статистическими методами позволяет сформировать наиболее эффективную структурную схему измерительного блока ИЭСД, обеспечивающую заданные нормативно-технические требования к диагностированию технического состояния ДВС.
Благодарности. Государственная работа. Фундаментальные и прикладные исследования в рамках реализации Программы фундаментальных научных исследований государственных академий наук на 2013-2020 годы.
Список литературы / References
[1] Савченко О.Ф., Добролюбов И.П., Альт В.В. Особенности и принципы построения экспертной системы оценки состояния ДВС // Современные методы, средства измерений и автоматизации для исследований физических процессов в сельском хозяйстве: сб. науч. тр. / РАСХН. Сиб. отд-ние. СибФТИ. Новосибирск, 1993. С. 59-73.
Savchenko, O.F., Dobrolyubov, I.P., Alt, V.V. Design features and principles for the expert system determining the condition of internal combustion engine // Modern Methods, Measurement and Automation Facilities for Researching the Physical Processes in Agriculture: Collection of Science Papers / RASHN. Sib. otdelenie. SibFTI. Novosibirsk, 1993. P. 59-73. (in Russ.)
[2] Добролюбов И.П., Савченко О.Ф., Альт В.В. Измерительная экспертная система для определения технического состояния двигателей внутреннего сгорания // Приборы и системы управления. 1998. № 12. С. 56-59.
Dobrolyubov, I.P., Savchenko, O.F., Alt, V.V. Measurement expert systems for determining the condition of internal combustion engine // Pribory i Sistemy Upravleniya. 1998. No. 12. P. 56-59. (in Russ.)
[3] Савченко О.Ф., Добролюбов И.П. Классификация состояния объектов статистическими методами с помощью измерительной экспертной системы // Методы и технические средства исследований физических процессов в сельском хозяйстве: c6. науч. тр. / РАСХН. Сиб. отд-ние. СибФТИ. Новосибирск, 2001. С. 70-81.
Savchenko, O.F., Dobrolyubov, I.P. Classification of objects condition by statistical methods with measurement expert system // Methods and Technical Means for Researching the Physical Processes in Agriculture: Collection of Science Papers / RASHN. Sib. otdelenie. SibFTI. Novosibirsk, 2001. P. 70-81. (in Russ.)
[4] Добролюбов И.П. Частотный анализ рабочих процессов ДВС с помощью измерительной экспертной системы // Тракторы и сельхозмашины. 2004. № 10. С. 17-19. Dobrolyubov, I.P. Frequency analysis of work processes in internal combustion engine with measurement expert system // Tractory i Sel'khozmashiny. 2004. No. 10. P. 17-19. (in Russ.)
[5] Савченко О.Ф. Экспертиза технического состояния ДВС на основе интеллектуального анализа данных // Труды ГОСНИТИ. 2008. Т. 101. С. 96-103.
Savchenko, O.F. Expertize of the technical condition of internal combustion engine based on data mining // Trudy GOSNITI. 2008. Vol. 101. P. 96-103. (in Russ.)
[6] Добролюбов И.П. Экспертиза состояния машин корреляционными методами // Тр. междунар. конф. "Энергообеспечение и энергосбережение в сельском хозяйстве". М.: ГНУ ВИЭСХ, 2010. Т. 1. С. 447-451.
Dobrolyubov, I.P. Expertize of the machinery's condition by correlation methods // Papers of Intern. Conf. "Power Supply And Power Economy in Agriculture". M.: GNU VIESKH, 2010. Vol. 1. P. 447-451. (in Russ.)
[7] Добролюбов И.П., Савченко О.Ф., Ольшевский С.Н. Экспертиза состояния дизель-генератора по статистическим характеристикам // Тр. междунар. конф. "Энергообеспечение и энергосбережение в сельском хозяйстве". М.: ГНУ ВИЭСХ, 2014. Т. 1. С. 185-190. Dobrolyubov, I.P., Savchenko, O.F., Olshevskiy, S.N. Expertize of the condition of a diesel generator by statistical characteristics // Papers of Intern. Conf. "Power Supply and Power Economy in Agriculture". M.: GNU VIESKH, 2010. Vol. 1. P. 185-190. (in Russ.)
[8] Добролюбов И.П., Савченко О.Ф. Выбор информативных признаков при использовании измерительной экспертной системы двигателя // Измерительная техника. 2005. № 2. С. 18-21.
Dobrolyubov, I.P., Savchenko, O.F. The choice of informative characteristics when using the expert measuring system of an engine // Measurement Techniques. 2005. Vol. 48, No. 2. P. 128-132.
[9] Добролюбов И.П., Савченко О.Ф. Выбор совокупности косвенных диагностических параметров для измерительной экспертной системы ДВС // Двигателестроение. 2012. № 2. С. 30-33.
Dobrolyubov, I.P., Savchenko, O.F. Selection of complex indirect diagnostics parameters for measurement expert system of an internal combustion engine // Dvigatelestroyeniye. 2012. No. 2. P. 30-33. (in Russ.)
[10] Добролюбов И.П., Савченко О.Ф. Оптимизация метрологического обеспечения измерительных экспертных систем ДВС // Тр. междунар. конф. "Энергообеспечение и энергосбережение в сельском хозяйстве". М.: ГНУ ВИЭСХ, 2010. Т. 1. С. 441-443. Dobrolyubov, I.P., Savchenko, O.F. Optimization of metrological supply for measurement expert systems of ICE // Papers of Intern. Conf. "Power Supply and Power Economy in Agriculture". M.: GNU VIESKH, 2010. Vol. 1. P. 441-443. (in Russ.)
[11] Добролюбов И.П. Повышение точности измерения характеристик двигателей внутреннего сгорания// Измерительная техника. 2007. № 2. С. 48-51.
Dobrolyubov, I.P. Increasing the precision of measurements of the characteristics of internal combustion engines // Measurement Techniques. 2007. Vol. 50. No. 2. P. 173-178.
[12] Добролюбов И.П., Савченко О.Ф., Ольшевский С.Н. Оптимизация обнаружения и измерения параметров ДВС измерительной экспертной системой // Ползуновский вестник. 2011. № 2/2. С. 275-279.
Dobrolyubov, I.P., Savchenko, O.F., Olshevskiy, S.N. Optimization of determining and measurement of ICE parameters by measurement expert system // Polzunovsky Vestnik. 2011. No. 2/2. P. 275-279. (in Russ.)
[13] Добролюбов И.П., Савченко О.Ф., Альт В.В., Ольшевский С.Н. Разработка компьютерной настраиваемой модели двигателя внутреннего сгорания // Вычисл. технологии. 2013. Т. 18. № 6. С. 54-61.
Dobrolyubov, I.P., Savchenko, O.F., Alt, V.V., Olshevskiy, S.N. Development of a modifiable computer model for internal combustion engine // Computational Technologies. 2013. Vol. 18. No. 6. P. 54-61. (in Russ.)
[14] Добролюбов И.П., Савченко О.Ф., Ольшевский С.Н. Принципы разработки компьютерной динамической модели автотракторных ДВС // Вестник НГАУ. 2014. № 2. С.141-146.
Dobrolyubov, I.P., Savchenko, O.F., Olshevskiy, S.N. Principles to design a computer dynamic model of autotractor ICE // Bulletin of NSAU. 2014. No. 2. P. 141-146. (in Russ.)
[15] Савченко О.Ф., Ольшевский С.Н., Добролюбов И.П. Информационная технология определения технического состояния тракторных двигателей // Механизация и электрификация сельского хозяйства. 2010. № 11. С. 27-30.
Savchenko, O.F., Olshevskiy, S.N., Dobrolyubov, I.P. Information technology for determining the technical condition of tractor's engine // Agriculture Mechanization and Electrification. 2010. No. 11. P. 27-30. (in Russ.)
[16] Техническое обеспечение измерительных экспертных систем машин и механизмов в АПК / В.В. Альт, И.П. Добролюбов, О.Ф. Савченко, С.Н. Ольшевский. Новосибирск: ГНУ СибФ-ТИ СО РАСХН, 2013. 523 с.
Technical support for measuring expert systems for machinery and devices in agriculture / Alt, V.V., Dobrolyubov, I.P., Savchenko, O.F., Olshevskiy, S.N. Novosibirsk: GNU SibFTI SO RASKHN, 2013. 523 p. (in Russ.)
[17] Патент RU 2428672 С2. МПК G01M15/04. Способ определения технического состояния двигателей внутреннего сгорания и экспертная система для его осуществления /И.П. Добролюбов, В.В. Альт, О.Ф. Савченко, С.Н. Ольшевский. Бюл. № 25. 10.09.2011.
Patent RU 2428672 С2. MPK G01M15/04. Method of determining the technical condition of ICE and expert system for its implementation / I.P. Dobrolyubov, V.V. Alt, O.F. Savchenko, S.N. Olshevskiy. Bulletin No. 25. 10.09.2011. (in Russ.)
[18] Фалькович С.Е., Хомяков Э.Н. Статистическая теория измерительных систем. М.: Радио и связь, 1981. 288 с.
Falkovich, S.E., Homyakov, E.N. Statistical theory of measurement systems. М.: Radio i Svyaz, 1981. 288 p. (in Russ.)
[19] Ширман Я.Д., Манжос В.Н. Теория и техника обработки радиолокационной информации на фоне помех. М.: Радио и связь, 1981. 416 с.
Shirman, Y.D., Manzhos, V.N. Theory and engineering for processing the radiolocation information with background noise. М.: Radio i Svyaz, 1981. 416 p. (in Russ.)
[20] Купер Дж., Макгиллем К. Вероятностные методы анализа сигналов и систем: Пер. с англ. М.: Мир, 1989. 376 с.
Cooper, G.R., McGillem, C.D. Probabilistic methods for signal and system analysis. New York: Oxford University Press, 1999. 496 p.
Поступила в 'редакцию 9 сентября 2015 г.
Modelling the process of optimal determination for the parameters of the conditions of internal combustion engine by measurement expert system
Dobroliubov, I.P., Savchenko, Oleg F.*, Olshevsky, S.N., Alt, V.V., Klimenko, D.N.
Siberian Physical-Technical Institute of Agrarian Problems, Novosibirsk, 630501, Russia Corresponding author: Savchenko, Oleg F., e-mail: [email protected]
Goals and methods: The problem of improvement of the efficiency for internal combustion engine (ICE) in operational conditions as well as the identification of its technical condition with the help of the measurement expert system is discussed. The
© ICT SB RAS, 2015
research and justification of the optimal method of detection of structural parameters of the engine, hidden in the high level of other components of the workflow was performed using known probabilistic methods of signal analysis, statistical methods of optimal signal processing on the background noise.
Results: The modelling of process for optimal estimate of informative parameters was conducted using a priori data. These parameters describe the operating process of ICE, which include: 1) the amplitude of angular acceleration of the crankshaft measured within the angular interval of operation of the specific cylinder during the cycle; 2) the signal that describes the advance angle of fuel injection; 3) jointly measured ICE parameters. It was revealed that the estimation of the measured information parameters of the signals will be optimal when the requirement of standard deviation minimum for these parameters is met, and that is equivalent for minimum average risk criteria.
Consider an optimal detection and measurement conditions for weak and strong interference amplitude of the angular acceleration of the crankshaft, measured in an angular range of operation of the cylinder in the cycle. It was found that in the case of a weak Gaussian noise (for example, the examination of four-cylinder engines) best estimate is equal to the measured value. In this case, a process corresponds to the measured center of gravity curve posterior probability density. Therefore it is advisable to use the signal accumulators and to distinguish the signals from the noise via filtration. Similar studies and conclusions are valid for the signal representing the advance fuel injection angle, as well as co-processing of the signal from the previous signal. It should, if possible, be used as a threshold (reference) values previously obtained from a priori statistical characteristics and parameters reflecting the different classes of states of the internal combustion engine (nominal, valid, pre-fault, fault).
Conclusion: The numerical modelling of the optimal detection and measurement of the background noise using statistical methods allows creating the most effective block diagram of the measuring unit. It provides the desired regulatory and technical requirements for diagnosing ICE's technical condition.
Keywords: ICE, resource structure parameters, signal, discovering, measurement, likelihood function, threshold level, optimal solution rule.
Acknowledgements. State engagement. Fundamental and applied research within the Fundamental Scientific Research Program for State Academies of Sciences for 2013-2020 years.
Received 9 September 2015
