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

Управление качеством

..pdf
Скачиваний:
18
Добавлен:
15.11.2022
Размер:
22.01 Mб
Скачать

Таблица 3 . 2

Значения максимальной погрешности ∆ S определения S

N, шт.

S, %

kσ

N, шт.

S, %

kσ

25

42,4

1,40

200

15,0

1,15

50

30,0

1,30

300

12,2

1,12

75

25,0

1,25

400

10,6

1,11

100

21,2

1,20

500

10,0

1,10

Вероятный процент брака вычисляется следующим образом. При рассеянии размеров, соответствующих закону нормального распределения Гаусса, принимается с погрешностью не более 0,27 %, что все детали партии имеют действительные размеры в пределах поля рассеяния:

6σ = xmax xmin,

где xmax, xmin – максимальное и минимальное значения параметра (размера). При этом площадь, ограниченная кривой нормального распределения и осью абсцисс, равна единице и определяет 100 % заготовок партии. Площадь заштрихованных участков (рис. 3.17)

 

представляет собой количест-

 

во деталей, выходящих по

 

своим

размерам за

пределы

 

допуска.

 

 

 

 

Для определения

коли-

 

чества годных деталей необхо-

 

димо найти площадь, ограни-

 

ченную кривой и осью абсцисс

 

на длине, равной допуску δ.

 

При симметричном

располо-

 

жении поля рассеяния относи-

 

тельно

поля

допуска

следует

 

найти значение интервала, оп-

Рис. 3.17. К определению количества ределяющего

половину

пло-

годных деталей

щади,

ограниченной

кривой

 

Гаусса иабсциссой х1 (х2).

211

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

x

 

1

x

2

 

 

 

F (x) = y dx =

 

ex

2σ

2

 

 

 

 

dx .

(3.29)

 

 

 

 

−∞

σ

−∞

 

 

 

 

Для случая, когда x = 0, распределение σ = 1 называют стандартным и функция распределения (3.29) имеет вид (рис. 3.18)

 

1

x

 

F (x) =

ex2 2dx .

(3.30)

 

 

−∞

 

Таким образом, если случайная величина Х следует закону нормального распределения, то вероятность появления случайной погрешности определяется площадью, ограниченной кривой f (x) и ее частью и осью абсцисс,

P{x1

 

} =

 

1

x2

 

 

 

 

< x < x2

 

ex

2

2

 

 

 

 

dx .

(3.31)

σ

 

 

 

 

 

 

x1

 

 

 

 

Подынтегральное значение есть элемент вероятности, равный площади прямоугольника с основанием dx и абсциссами x1 и x2, на-

зываемыми квантилями.

Рис. 3.18. Функция распределения F(x) и функция Лапласа Ф(x)

212

Произведем замену переменной: t = x / σ , dx = σ

dt,

P{x1 < x < x2 } =

1

t2

 

et2 2dt .

(3.32)

 

 

t1

 

Представим правую часть в виде суммы двух интегралов:

P{x1 < x < x2 } =

1

0

1

t2

et 2 2dt +

et 2 2dt .

 

 

t1

0

Интеграл вида

t

Ф(t ) = 1 et 2 2dt (3.33) 2π 0

носит название нормальной функции Лапласа. Значения этого интеграла сведены в табл. П. 3. Таким образом, указанная вероятность (3.31) сводится к разности нормальных функций Лапласа:

Р{ x1 < x < x2 } = Ф(t2) – Ф(t1).

(3.34)

Расчет количества годных деталей сводится к установлению величины t и определению Ф(t) по таблице с последующим пересчетом полученных величин в проценты или в число штук изделий.

В общем случае, когда x ≠ 0 , имеем следующую вероятность появления случайных погрешностей:

 

x2

 

x2

 

 

 

P{x1 xx2 }=

f (x ) dx= Ф

x

 

 

 

 

 

 

 

 

 

 

 

x1

 

 

σ

Фx1 x

. (3.35)

σ

Отметим свойства функции Лапласа: Ф(0) = 0; Ф(– х) = – Ф(х) (функция нечетная); Ф(∞) =1/2. Из рис. 3.18 видно, что кривые F(х) и Ф(x) эквидистантны.

Если в равенстве (3.354) положить х1 = – ∞, то

P (x x2 )=

1

+ Ф

x2

 

 

 

 

 

x

,

(3.36)

 

 

σ

 

2

 

 

 

 

 

так как Ф(– ∞}= – Ф(∞) = –1/2. Положив в соотношении (3.35) х2 = ∞, находим:

213

P (x ≥ x1 )=

1

− Ф

x1

 

 

.

 

x

(3.37)

 

 

 

 

2

σ

 

 

 

Пример 1. При измерении сопротивлений делителя напряжения установлено, что среднее значение этого сопротивления R = 5,5 кОм, а среднее квадратическое отклонение σ= 1,5 кОм. Принимая нормальный закон распределения, найти вероятность появления сопротивлений свыше 10 кОм.

Решение. По равенству (3.37) и из таблиц (П3) находим:

Р(R > 10) = 1/2 – Ф[(10 – 5,5)/1,5] = 0,5 – 0,4986 = 0,0014.

Пример 2. Определить количество бракованных и годных деталей, если допуск на обработку δ = 0,10 мм. Среднее квадратическое отклонение S = 0,02 (получено по результатам замеров 75 штук). Общее количество обработанных деталей – 300 шт.

Решение.

1.Определяемрасчетноезначениеσ= kσ S = 1,25 0,02 = 0,025 мм.

2.Полефактическогорассеянияω= 6 σ= 6 0,025 = 0,15 ммпревосходит поле допуска δ = 0,1 мм; следовательно, условие обработки без брака не выполнено ипоявление брака возможно.

3.x = δ/2 = 0,1/2 = 0,05; t = x/σ = 0,05/0,025 = 2,0. Ф(t) = 0,4772,

что соответствует 47,72 % годных деталей для половины партии. Для

всей партии количество годных деталей – 95,44 %, или 286 шт., а бракованных 4,56 %, или 14 шт.

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

214

3.13.2. Экспоненциальное распределение

Экспоненциальное распределение характерно для внезапных отказов элементов и систем. Плотность вероятности экспоненциального распределения задается уравнением

f (x ) = λeλx , F (x) =1eλx , λ(x) = λ , x > 0,

(3.38)

где λ − параметр распределения, являющийся строго положительной константой.

Рис. 3.19. Графики плотности f(x), интенсивности отказов λ (x) (а) и функции F(x) экспоненциального распределения (б)

Среднее значение x и среднеквадратическое отклонение σ экспоненциального распределения совпадают и равны обратному значению параметра: x = σ = 1/λ . Графики функций F(х) и f (x) приведены на рис. 3.19. Отличительной особенностью экспоненциального распределения является то, что интенсивность отказов λ (х) постоянна, т.e. не зависит от аргумента (значения случайной величины).

3.13.3. Закон равной вероятности

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

215

Рис. 3.20. Распределение случайной величины по закону равной вероятности

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

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

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

При изменении случайной величины X в интервале от a до b плотность f (x) постоянна и равна С; вне этого интервала она равна нулю. Поскольку площадь, ограниченная кривой распределения, равна единице,

b

C (b a) =1 , отсюда C =

1

 

 

f (x)dx =1,

.

(3.39)

 

a

 

 

 

b a

 

Плотность распределения f (x) имеет вид:

 

 

 

f (x ) = (b a)1 ,

при

a xb

 

 

(3.40)

 

0,

при

x > b; x < a; x [a, b].

 

216

Функция распределения (рис. 3.21) и интенсивность отказов (см. рис. 3.20) имеют вид:

 

0,

 

 

при x < a ,

 

 

 

a

 

 

 

x

 

 

F (x) =

 

 

 

 

, при a

x

b ,

 

 

 

 

b

a

 

 

 

1,

 

 

при x > b ,

 

 

 

 

 

 

 

 

λ(x ) =

1

 

 

при a

x

b .

 

 

 

 

 

b x

 

 

Вычислим математическое ожидание и дисперсию:

(3.41)

(3.42)

 

 

 

 

+∞

 

 

b

 

 

1

 

 

 

 

 

1

 

 

x

2

 

 

b

 

a + b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M x =

x f (x)dx =x

 

 

 

dx =

 

 

 

 

 

 

 

 

=

,

 

 

 

 

 

 

 

b a

 

a 2

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

a

 

 

 

b

 

a

2

 

 

 

 

 

 

 

+∞

 

 

 

 

 

 

 

 

 

 

 

+∞

 

 

 

a + b

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

Dx = (x M x )

 

f (x)dx =M x ( X

 

) M x

= x

 

 

f (x ) dx

 

 

 

 

=

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

b

2 1

 

 

(a + b)2

 

1 x3

 

b

a + b 2

 

(b a )2

 

 

 

 

 

 

 

 

 

 

 

 

 

= x

 

 

 

dx

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

a

 

b a

 

 

b a 3

 

a

 

 

 

 

 

 

 

 

 

 

 

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определяем среднее квад-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ратическое

 

отклонение

и

 

поле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

рассеяния:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ =

b a

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω= b a = 2σ 3 .

(3.43)

Рис. 3.21. График функции F(x)

 

 

 

 

 

 

 

Коэффициент асимметрии

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

 

 

 

 

 

Sk =

0

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

симмет-

рично). Для определения коэффициента эксцесса найдем четвертый центральный момент:

217

 

 

1

b

 

a + b

 

4

dx = (b a)

4

M4

=

x

 

 

.

 

2

 

 

 

b a a

 

 

 

80

 

Отсюда Ek = Mσ44 3 = −1, 2.

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

3.13.4. Закон Релея (эксцентриситета)

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

Распределение по закону Релея формируется тогда, когда случайная величина Rс представляет собой геометрическую сумму двух случайных величин Х и У (рис. 3.22), т.е.

Rс = X 2 +У2 ,

(3.44)

каждая из которых подчиняется закону Гаусса.

Закон Релея однопараметрический, и уравнение его кривой распределения имеет вид

y =

Rс

eRс2 2

,

(3.45)

σ2

 

 

 

 

где σ – среднее квадратическое отклонение значений координат х и у.

218

Рис. 3.22. Образование эксцентриситета

Рис. 3.23. Теоретическая

(радиуса-вектора R) втулки 1 при ее

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

обработке на цилиндрической оправке 2

закону Релея

при наличии зазора между оправкой

 

и отверстием втулки

 

Для теоретической кривой распределения по закону Релея (рис. 3.23) характерны крутой подъем восходящей ветви и более пологий спуск нисходящей ветви, вершина кривой более заострена, чем у кривой нормального распределения, и смещена от среднего значения переменной величины Rc в сторону начала координат. Из уравнения (3.45) следует, что при Rc = 0 y = 0, т.е. начало кривой распределения эксцентриситета совпадает с началом координат. Нисходящая ветвь этой кривой асимптотически приближается к оси абсцисс.

Основными параметрами закона Релея являются:

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

личины;

σRс среднее квадратическое отклонение Rc;

σ среднее квадратическое отклонение значений коорди-

нат х и у конца радиуса-вектора Rc.

Они связаны между собой следующими соотношениями:

σ = σRс / 0,655; Rср = 1,92 σRс = 1,257 σ.

(3.46)

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

219

ω = 5,252 σRс = 3,44 σ.

(3.47)

При распределении Релея, когда фактическое поле рассеяния превосходит поле допуска (ω > δ ), возможно появление брака

(см. рис. 3.23).

Общую площадь F(Rc), ограниченную кривой распределения, находят по интегральному закону распределения эксцентриситета:

F (Rc )

 

1

Rc

 

 

Rc

2

 

 

=

Rc exp

 

 

dRc ,

(3.48)

2

2

 

 

σ

0

 

 

 

 

 

который после подстановки величин t = Rc/σ принимает нормированный вид:

F ( Rc ) =1 et2 2

(3.49)

и табулируется аналогично функции Лапласа (табл. П. 4).

Пример. Рассчитать вероятный процент брака, если допуск на изготовление детали равен δ = 0,04 мм. В результате непосредственных измерений первых 25 деталей установлено среднее квадратическое отклонение S = 0,009 мм.

Решение. Расчетное значение среднего квадратического отклонения находим по формуле σ = kσ S и табл. 1.3: σRc = σ = kσ ·S =

= 1,4·0,009 = 0,0126 мм.

Фактическое поле рассеяния значений эксцентриситета находим по формуле (3.47):

ω = 5,252σRc = 5,252·0,0126 = 0,0662 мм.

При Rc= δ = 0,04 мм и t = 0,655δ Rc = 0,655·0,04/0,0126 = = 2,08. В соответствии с табл. П. 4 Ф(t) = 0,8851, т.е. количество годных деталей составляет 88,51 % и количество брака – остав-

шиеся 11,49 %.

220