<<
>>

4.1 Разработка алгоритма для реализации метрологического обеспечения контроля и диагностирования АТС.

Значения контролируемых показателей есть не что иное,.как случайный вектор, который характеризуется n-мерной плотностью вероятности fl(xi;...;x„), причем, вид функции f известен (определён из рассмотрения условий эксплуатации и производства АТС, либо путём обработки статистических данных), но каждый раз вид функции является априорной информацией о контролируемых показателях.

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

тов составляющих погрешности =4Ki4t + двлу +д^ . Данную процедуру можно представить в виде:

е А

K4j) =

0 « А

(4.1)

6 я'; otf™,« 0';

Области D и D' могут являться n-мерными параллелепипедами, т. е. А, < ?,• < Bj и А < А{ < $ < В/< В:

/•I

(4.2)

п

Здесь //, и Н, задаются подобно (4.2)

(4.3)

Учитывая введенные обозначения, имеем:

t' 2

Л=1-

(4.4)

Описание алгоритма (рис.4.1):

В первом блоке осуществляется запуск программы на выполнение

В блоке 2 обнуляются вспомогательные переменные 1|, 12, Ь, порядковой переменной i присваивается 1, где i=l,...,n - номер элемента в векторе измерений контролируемого параметра.

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

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

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

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

В блоке 10 осуществляется учет наработки при измерении текущего значения контролируемого параметра

В блоке 11 определяются величина контролируемого параметра с учетом эксплуатационной погрешности

9 В блоке 12 подсчитываются значения Н(х) и h(x)

10. В блоке 13 рассчитывается значение вспомогательной переменной

l=H*h

Рис. 4.1. Алгоритм расчета ошибок первого и второго рода при прямом многопараметрическом контроле.

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

В блоке 15 происходит прибавление к переменной I единице.

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

величин ошибок первого и второго рода в блоке 16 для данного контролируемого параметра.

1

C™ED

11-0

12=0 0*0 I -yoyimn [

еип:

I U. t»A<

¦Л

есрмироммм X Фсуциро» но о^.

11

¦у i

СбЬМЮ» Омам

t-Hh 13

IS |

kv. В,

21

R2»1-IW

17

19

w

Х*МодйП»

го :

\ I

2*_

СТОП » <^Нц|Ч<(>7+>

Ж | <Ч-<,1,<В|«0,*ЯЯч1п[П){<'—1 | 0»-0»*(В/0д)г ( 1

Рис. 4.2. Алгоритм расчета ошибок первого и второго рода при косвенном многопараметрическом контроле.

Описание алгоритма (рис.4.2):

1.В первом блоке осуществляется запуск программы на выполнение 22. В блоке 2 обнуляются вспомогательные переменные 12,1з, порядковой переменной i присваивается 1, где i=l,...,M - номер контролируемого показателя

В блоке 5 переменной] присваивается значение 1, где j=l, п - номер контролируемого параметра.

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

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

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

В блоке 17 вычисляется заносится в массив j-oe значение среднее квадратическое отклонение полной эксплуатационной погрешности для текущего контролируемого параметра

В блоке 19 осуществляется учет наработки для данного значения контролируемого параметра

В блоке 20 вычисляется и осуществляется ввод в массив j-ro измеренного значения контролируемого параметра с учетом абсолютной эксплуатационной погрешности

В блоке 22 осуществляется прибавление к переменной] единицы

11 В блоке 23 проверяется условие окончания вектора параметров и если оно выполняется, то управление передается блоку 25, где вычисляется значение контролируемого показателя из известной зависимости.

В блоке 26 переменной] присваивается значение единицы (управление передается первому значению в массиве параметров х)

В блоке 27 определяется процедура вызывается процедура расчета множественного коэффициента корреляции для данного параметра Xj

В блоке 28 проверяется условие коррелированности данного параметра с остальными, и если это так, то значение среднего квадратического отклонения контролируемого показателя рассчитывается в блоке 30 по формуле в блоке 30, иначе - по формуле в блоке 29.

Далее в блоке 31 управление передается следующему элементу в массиве параметра

6.

В блоке 32 проверяется условие окончания массива параметров. Если это условие выполняется, то в блоке 3 вычисляется квадратный корень из среднего квадратического отклонения контролируемого показателя.

В блоке 4 вводится квантиль нормального распределения для данного контролируемого показателя

В блоке 6 вычисляется общая погрешность косвенного измерения для данного показателя

В блоке 8 определяется измеренное значение показателя

В блоке 10 подсчитываются значения Н(х) и h(x)

В блоке 12 рассчитывается значение вспомогательной переменной l=H*h

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

В блоке 16 пополняется на единицу счетчик номера контролируемого показателя

В блоке 18 происходит проверка условия достижения конца вектора контролируемых показателей и если оно выполняется, то в блоке 21 происходит вычисление величин ошибок первого и второго рода для данного контролируемого показателя. (?)

Видом* to МПОМК хммп) X, и X. Ргобт мрояжспЛ Р-.ИР* 18 •нпоро Л ОЧС1ИМ

пегрммосги ДХ 18 Г) Вдета

МГТС ¦MWK1

pax- 73 Рдеет нролхостей

Р. «Р, 1?

П9Л)чем« nrcpi iwm мА rac«verpoaX

Cnnuwnmw P.-.IX) Неиспрмностъ

ьвгцвна (myuW 10 JMM1 У» га К»1

Рис. 4.3. Алгоритм расчета ошибок первого и второго рода при использовании алгоритма диагностирования.

Описание алгоритма расчета ошибок первого и второго рода при использовании алгоритма диагностирования (рис.4.3.):

Г ПУСК J

Ун—.

/ ВИА»«« /

/ aenvfii / / Р' /

В блоке 2 происходит введение вектора априорных вероятностей отказа того или иного узла АТС. В данном векторе априорные вероятности поставлены в вариационный ряд по убыванию.

В блоке 4 инициализируется цикл перебора всех узлов общим числом т.

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

В блоке 9 осуществляется генерация вектора значений параметров {X} характеризующих данный узел.

Затем в блоке 11 происходит выделение из вектора параметров {X} значений {Хг}, соответствующих условию работоспособности и {Хнг}, соответствующих условию неработоспособности.

Далее в блоке 13 производится расчет вероятностей методических ошибок первого и второго рода Pim и Р2м соответственно.

В блоке 16 происходит получение вектора значений погрешности измерения {ДХ} для выбранного узла.

В блоке 18 рассчитывается вектор измеренного значения параметра {X'}, из которого выделяются (блок 20) значения векторов {Х'г} и{Х;нг}, соответствующих условиям работоспособности и неработоспособности соответственно, применительно к измеренным значениям.

В блоке 23 вычисляются вероятности ошибок первого Р| и второго рода Р2 с учетом измерительной погрешности, а также эксплуатационных ошибок первого и второго рода Phl(Ct, и Рыхп (блок 3).

В блоке 5 проверяются условия верного вынесения диагноза при неисправности проверяемого узла. В случае положительного результата условия, делается вывод о неисправности текущего узла, что влечет за собой замену данного узла и проверку работоспособности АТС в блоках 10 и 14.

Если же результат (в блоке 5) отрицательный, то в блоках 8 и 12 проверяются условия неверности вынесения заключения об исправности или неисправности соответственно. В случае обнаружения неверного вынесения заключения, измерения данного узла повторяются с блока 6.

Затем осуществляется исключение максимального значения из вектора априорных вероятностей отказа и происходит переход на проверку следующего узла. Когда были измерены все узлы АТС и при этом неисправность АТС не выявлена, алгоритм прекращается.

Далее для осуществления расчета критерия (2.53) необходимо воспользоваться алгоритмом, показанным на рис.4.4.

6 сГ 9 -5 * 1 1 + ч г 1 3 1 0 1 1 2 .1 1 4

"т:

I 5

1 6 з:

2 I

22

/ ¦> У 2. 23

Рис.4.4. Алгоритм расчёта критерия оптимизации метрологического обеспечения

На рис.4.4 приведен алгоритм реализации методики технико-

экономической оценки выбора средств контроля на ЭВМ. Рассмотрим операции данного алгоритма.

1. Ввод исходных данных: -для АТС - число подсистем АТС - R; - для каждого обобщенного показателя - п^ -для варианта реализации системы контроля -

-для каждого контролируемого параметра Са/, Г, аЛх, а также данные о законах распределения контролируемого параметра и погрешности измере-

ния, т.е. вид закона f(x)j, (p(Ax)j и Проверка условия: f(x)j, (p(Ax)j заданы (да, нет).

Условие п.2 не выполняется. Принять, что (р(Ах) и f(x) - законы равной вероятности.

Условие п.2 выполняется. Проверка условия: axj задано (да, нет).

Условие П.4 выполняется. Расчет Рц и P2J в зависимости от вида реализуемой системы диагностирования.

Условие п.4 не выполняется. Нахождение максимумов функций Р} = f(Km, KJ, Р2 =f(Km, KJ.

Принять эти максимальные значения за Лу и P2j.

Вывод Pij и P2j.

Проверка условия: задано crdortJ, или АС - ЩАх) - Yfcr^) и Пдоп.

Задано стдщ. Приравнять сгл; = crdonj.

Задано АС = Ф(Ах) = и П^. Расчет выражения (1.35).

Решить уравнение (2.39) относительно о^.

Приравнять Gxj = (ТдопГ

Расчет Pldonj и P2donj по выбранным формулам.

Расчет PH.rdonj как сумму Р}доп] и P2dРасчет Рюупо формуле (2.40).

Расчет полноты контроля npi по формуле (2.38).

Проверка условия (2.37).

Условия (2.37) не выполняется, вывод сообщения об этом.

Расчет коэффициента К} по формуле (2.51).

Расчет значений критерия Q по формуле (2.53).

Оценка относительной погрешности расчета критерия Q - Sq, ПО формулам (3.15)-(3.16).

Выбор минимального Q.

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

<< | >>
Источник: Исакова Кира Сергеевна. ИССЛЕДОВАНИЕ И СОВЕРШЕНСТВОВАНИЕ МЕТРОЛОГИЧЕСКОГО ОБЕСПЕЧЕНИЯ ДИАГНОСТИРОВАНИЯ ПРИ ТЕХНИЧЕСКОМ ОБСЛУЖИВАНИИ И РЕМОНТЕ АВТОТРАНСПОРТНЫХ СРЕДСТВ (НА ПРИМЕРЕ ЭЛЕКТРООБОРУДОВАНИЯ). Диссертация на соискание ученой степени кандидата технических наук. Специальность 05.22.10-Эксплуатация автомобильного транспорта. Владимир - 2007. 2007

Еще по теме 4.1 Разработка алгоритма для реализации метрологического обеспечения контроля и диагностирования АТС.:

  1. ВВЕДЕНИЕ
  2. 1.1. Обзор способов и методов разработки метрологического обеспечения контроля и диагностирования технического состояния автотранспортных средств.
  3. 2.3. Разработка методики оценки характеристик достоверности прн использовании алгоритмов диагностирования с учетом методической составляющей погрешности, погрешности измерения н дополнительной погрешности.
  4. 4.1 Разработка алгоритма для реализации метрологического обеспечения контроля и диагностирования АТС.
  5. 4.3 Разработка метрологического обеспечения системы диагностирования технического состояния фар автотранспортных средств в режимах «ближний свет», «дальний свет» и «суммарный свет» при реализации стандартного метода измерений.