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

книги из ГПНТБ / Регулирование качества продукции средствами активного контроля

..pdf
Скачиваний:
6
Добавлен:
27.10.2023
Размер:
22.09 Mб
Скачать

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

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

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

функции. Корреляционная функция

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

зависимость между

сечениями

случайной функции,

относящимися

к различным

значениям

параметра

t.

 

 

 

 

Определим значения

функции X(t), которые она принимает в моменты

време­

ни tu

U,

• • .,

tm,

т. е. для

различных сечений случайной

функции. Каждому из зна­

чений

ti,

ti,

. . . , tm будет соответствовать п

значений

случайной

величины

(соб­

ственно случайных погрешностей). Величины интервалов между моментами време­

ни ti, h, . . . ,

tm обычно задают равноотстоящими друг от друга

и выбирают так,

чтобы была ясна основная тенденция изменения функции X{t).

Зарегистрирован­

ные значения

функции X (t) заносят в таблицу (табл. 1), к а ж д а я строка которой

соответствует определенной реализации, а число граф равно числу значений ар­

гумента (параметра t). Таким образом, в таблицу сведены результаты

п опытов

над системой

m случайных величин X(ti),

X(h),

. . . ,

X(tm).

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 1

 

Значения вероятностных характеристик

функции X

(t)

 

X(t)

 

h

 

h

tk

 

 

ti

 

tm

Xi(t)

A't

(*,)

Xi

(tt)

X,

(tk)

Xx(ti)

X,

(tm)

Xt(t)

X2

(h)

Хг

Ш

X2

(h)

Xt

(td

Xt(tm)

 

Xi

(h)

Xi

(t2)

Xi

{tk)

Xi

(tù

Xi

(tm)

xn(t)

x n

Ci)

Xn

(h)

X„(tk)

Xn

(ti)

Xn

(tm)

По полученным значениям определяют основные вероятностные характеристи­ ки случайной функции:

математическое ожидание

('*) = —

;

(2)

дисперсию

Dx (tk) = —

;

(3)

28

среднее квадратическое отклонение

 

 

 

 

 

 

 

 

корреляционную

функцию

 

 

 

 

 

 

 

 

 

J [Xi (tk) - m ,

Ш

• \ХІ

(tt) -

mx

(*,)]

 

 

Kx (tk,

ti) -=

n — I

 

 

 

.

 

(5)

 

 

 

 

 

 

 

 

После этого, пользуясь рядом значений mx(ti),

mx(h),

 

 

mx(tm),

можно

построить график изменения величины mx(t)

(см. рис. 4).

 

Dx

(tu) и Кх

(tk,

ti).

Можно также графически выразить изменение

значений

При необходимости эти графики аппроксимируются

какими-либо аналитическими

выражениями, например, по способу наименьших квадратов. Функция mx(th)

ха­

рактеризует собой усредненную функциональную погрешность.

 

 

 

У функциональных случайных погрешностей

параметр

Dx(tn)

в общем

случае

является переменным, т. е. величина рассеивания собственно случайных погреш­

ностей изменяется во времени

(или в функции

какого-либо другого параметра).

На практике встречаются

процессы, протекающие при постоянных значениях Dx(t)

и mx(t).

При этом вероятностные

характеристики случайной функции не зависят

от значения t. Такие процессы

изменения функции X(t) называют стационарными

случайными

процессами. В данном

случае величина поля рассеивания собственно

случайных

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

является

постоянной.

 

Если

при каком-либо

процессе

условие mx(t)

= const не удовлетворяется, то

это не мешает изучать его как стационарный случайный процесс (разумеется, если удовлетворяется условие Dx(t) = const). Корреляционная функция стационарного случайного процесса является функцией не двух, а одного параметра.

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

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

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

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

Причинами, вызывающими появление собственно случайных по­

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

п р и п у с к о Е

на обработку, неоднородность материалов и термической

обработки

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

29

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

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

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

тический характер (в этом случае зазоры

обычно «выбираются» в

одну сторону). Местные технологические

погрешности для одного

прибора также носят систематический характер.

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

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

С и л о в ы е д е ф о р м а ц и и технологической

системы

зависят

от жесткости узлов станка, количества и жесткости

стыков,

глуби­

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

Т е п л о в ы е д е ф о р м а ц и и

при шлифовании зависят от кон­

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

(в частности, шпиндельных опор),

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

33

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

При других видах обработки (точение, фрезерование и др.) ко­ личество определяющих факторов и их характер примерно такие

же. Влияние всех этих факторов для массы

станков, инструментов

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

характер. Действие

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

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

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

Рассмотрим, например, случайные функциональные погрешно­ сти обработки, которые могут быть охарактеризованы графиком, приведенным на рис. 3,а. Функциональные усредненные погрешно­ сти характеризуют собой среднюю величину влияния на размеры деталей износа инструмента, тепловых и силовых деформаций тех­ нологической системы. Вместе с тем указанные факторы вызывают одновременно и собственно случайные погрешности размеров, кото­ рые характеризуются флюктуациями размеров отдельных деталей относительно линии /-/. Эти флюктуации являются следствием из­ менения от детали к детали величин силовых и тепловых деформа­ ций технологической системы и износа инструмента под влиянием неодинаковых величин припусков на обработку, неодинаковости ма­ териала и термической обработки заготовок, изменения режимов резания и т. д.

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

31

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

Собственно случайные погрешности размеров чащего всего под­ чиняются закону нормального распределения (закону Гаусса). Как известно, условием формирования закона Гаусса является наличие достаточно большого числа случайно действующих и незначитель­ ных по величине факторов. Если число определяющих факторов больше четырех и они примерно одинаковы по величине, то резуль­ тирующим практически будет являться закон Гаусса.

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

Условием возникновения закона Гаусса является отсутствие среди определяющих факторов доминирующих. Последние при-

Рис. 5. Искажение кривых Гаусса под влиянием сильно действую­ щих факторов

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

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

32

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

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

сматриваемого параметра

(под номиналом понимается значение па­

раметра,

от которого

отсчитываются

случайные

отклонения). Как:

следует

из

графика,

при

исход­

 

 

ном

гауссовском распределении

-^т?

 

кривая распределения

 

существен-

"~

 

но положительных

случайных по

 

 

грешностей

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

собой

 

 

одну

ветвь

кривой

Гаусса.

 

 

На рис. 6, б

также

сплошной

 

 

линией

изображена

 

кривая рас­

 

 

пределения

существенно

положи­

 

 

тельных

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

при усло­

 

/

вии, что их центр

группирования

 

 

 

не совпадает с номиналом. Сме­

 

 

щение центра

группирования по

. п

о

отношению

к

номиналу

может

Рис. 6. Кривые распределения сущест­

быть

вызвано

влиянием

какого-

венно положительных погрешностей

либо фактора

(например,

влияни­

 

 

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

ния размеров шага

резьбы окажется смещенным по отношению к

его номинальному

значению. В этом случае кривая распределения

приобретает асимметричный характер (погрешности шага являются

существенно

положительными

при

определении

диаметральной

компенсации погрешности

шага) .

 

 

 

 

Возникающий при этом закон распределения

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

зависимостью

 

 

 

(х-а)1

 

_

(х+а)

 

 

 

1

 

 

 

 

У =

_

2о«

+

е

3

(6)

 

 

 

V2г.

s

 

 

 

 

 

где а—среднее

квадратическое

отклонение

исходного гауссовско-

го распределения;

 

 

 

 

 

а — смещение центра

группирования

исходного

распределения.

3—2891

33

§ 7. ХАРАКТЕРИСТИКИ ЦЕНТРА Г Р У П П И Р О В А Н И Я СЛУЧАЙНЫХ

ПО Г Р Е Ш Н О С Т Е Й

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

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

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

 

 

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

ошибкой среднего арифметического

 

 

 

°, =

^ =

(7)

где

а — средняя квадратическая

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

 

N — число измерений

одной и той же величины (размера).

 

Понятие «средняя квадратическая ошибка среднего арифмети­

ческого» распространяется также и на партию

измеряемых дета­

лей. В этом случае

 

 

 

 

]/

N

VN

 

где

Ou — среднее квадратическое

отклонение

собственно случай­

 

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

 

 

о„ — средняя квадратическая

погрешность

измерения;

 

N — число деталей в измеряемой партии.

 

Параметр os представляет собой среднее квадратическое откло­ нение ряда значений средних арифметических.

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

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

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

34

щей медиане, достаточно, чтобы 50% всех отклонений располага­ лось по одну сторону от этой точки, а 50% — п о другую.

При дискретном законе распределения в качестве медианы мож­ но принять любое значение ряда, промежуточное между хт-\ и хт,. если удовлетворяются следующие условия:

 

 

m 1

я

 

 

 

m—1

п

m

 

п

 

2)

2 Я (*,•)<

2 У > Р

( х

{ ) >

2 Я ( 4

(Ю)'

 

/=1

i—m

i l

 

i = m + l

 

Имеется

в виду, что значения

хи

х2, ....

хп расположены

в по­

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

ложим, что имеется проба из пяти обработанных

на станке

деталей

с размерами 10,12;

10,01; 9,98; 10,07

и 10,03 мм.

Если

полученные

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

 

(9,98;

10,01;

10,03;

10,07; 10,12 мм)

и

считать, что Р (х{)

= Р (х2)

=

Р (х3)

=

Р

4 ) =

= Р (хъ) = — ,

то

согласно второму

условию

в

качестве

медианы

5

 

 

 

 

 

 

 

 

можно принять любое число, занимающее промежуточное положе­

ние между х2 = 10,01 мм и х3

=

10,03

мм.

 

 

 

принимают

При

нечетном числе

данных

п за

медиану обычно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п + 1

число, занимающее

среднее

положение

по

порядку,

т. е. —-— .

В данном

примере

(при

п =

5)

это

положение

третье:

медиана

те = хз =

10,03 мм. Если бы пришлось

оперировать лишь

четырьмя

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

бы

принять

среднее

арифметическое

из

двух

членов:

х2

=

10,01

мм

и

хг

ш л о

 

 

о

 

 

 

х,

+

х3

 

10,01

+

10,03

=

І П

А О

= 10,03

мм. В этом случае те=

 

—— =

 

 

 

10,02 мм.

Аналогично

поступают при любом

четном числе п. Среднее

ариф­

метическое

значение в данном

примере

равно х =

10,04 мм.

Таким

образом, для дискретных случайных величин медиана в общем слу­ чае не совпадает со средним значением. Одно из важных свойств медианы заключается в том, что грубые отклонения размеров (гру­ бые погрешности измерения и обработки) не оказывают существен­ ного влияния на ее значение. В то же время грубые погрешности оказывают значительное влияние на величину среднего арифмети­ ческого. Это объясняется тем, что при определении положения ме­ дианы вес грубой ошибки таков же, как и у других значений разме­ ров, в то время как при определении среднего арифметического вес этой ошибки может быть значительно больше, причем в первом случае вес грубой ошибки не зависит от ее величины, а во втором он изменяется пропорционально значению ошибки. Упрощенно это

можно показать на следующем примере.

Предположим, что есть

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

возрастания дискретных

3*

35

случайных величин х 2 ; Х3;

х 4 ; Х5 и Хв, среди

которых

величины

Х \ ~ х ъ распределены

нормально,

 

х6 представляет собой грубую по­

грешность измерения или обработки.

 

 

 

 

Значение среднего примерно

равно

 

 

 

 

 

+

• • • + * » ~ , у 3 .

 

 

( П )

 

 

5

 

 

 

 

 

 

Значение медианы всей выборки, включая и грубую

погреш­

ность Хв, равно

 

 

 

 

 

 

 

 

 

ха "і"

ХІ

^ Х 3 .

 

 

 

(12)

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

Среднее для всей выборки составляет

 

 

 

 

Х\ ~\~ х2

+ • • • ~\~ х6

*6

 

 

/ 1 0\

 

 

6

 

1

6 *

 

 

 

Таким образом, разность величины медианы и среднего в пер­

вом приближении равна

 

 

 

 

 

 

 

 

6

 

 

 

6

/

6

 

поскольку величина

3

 

 

невелика.

 

Медиана является более

грубой (менее

достоверной)

характе­

ристикой положения центра группирования по сравнению со сред­ ним арифметическим.

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

рическая квантиль порядка Р, где Р = F (£), асимптотически

нор­

мальна и имеет

параметры

 

 

 

 

 

 

 

 

 

 

5 >

- і

- т

/ ^

,

 

 

 

(15)

где Р = 1 —- q;

 

 

ta)

V

N*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N — число деталей в выборке.

 

 

(квантиль с Р =

 

 

Рассмотрим

распределение

медианы

Ѵг)

для

гауссовского закона с

параметрами

ц

и и.

Заметим,

что

при

F (Ï) = Ѵг величина g =

ц..

 

 

 

 

 

 

 

 

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

 

 

 

f(x)

=

б

^

,

 

06)

 

 

 

• j / 2 -

с

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

 

/(5) =

/

Ы

=

-

 

 

(17)

 

 

 

 

 

К г *

а

 

 

 

36

Таким образом, медиана эмпирической функции распределения асимптотически нормальна и имеет параметры

2 \

2

N

У 2л s

т. е.

(16)

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

2

(19)

УVN

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

? ( * ) = —1. < > а -

(20)

а

 

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

Определим среднее квадратическое отклонение медианы для закона равной вероятности, а также для композиции законов нор­ мального и равномерного распределения.

Из предыдущего следует, что эмпирическая медиана выборки для закона распределения с плотностью / (х) и интегральной функ­ цией F (х) асимптотически нормальна и имеет параметры

«, т ^ г • - Л г .

(21)

Асимптотическое распределение медианы для закона равной вероятности в интервале (—/, + / ) характеризуется следующими зависимостями:

37