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

книги из ГПНТБ / Головинский В.В. Статистические методы регулирования и контроля качества. Расчет оптимальных вариантов

.pdf
Скачиваний:
7
Добавлен:
24.10.2023
Размер:
10.75 Mб
Скачать

среднему

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

отклонению

сг2б =

0,5 неискажен­

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

б)

уменьшением среднего

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

отклонения с сг2б =

0,5

до

ст2р = 0,4472.

 

 

Напомним, что для перехода от статистических к техническим

единицам измерения, надо помножить м.

о. z на ах. Например,

если <ух =

0,014 мм,

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

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

в сторону каждой из границ регулирования равняется 0,007 мм. Если ввести корректирующее правило, в соответствии с ко­ торым при нарушении правой границы регулирования оценка £ уровня настройки X уменьшается на м. о. zpax = 0,007 мм (и увеличивается соответственно при нарушении левой границы) — искажения информации в примере [и вообще при гауссовом рас­ пределении P+ (z) и Р“ (z) ] практически устраняется. В резуль­ тате расчет распределения ошибки настройки при настройке уточнениями выполняется так, как при независимой настройке. Практически, после проверки м. о. zB и ст2в указанным спо­ собом на обычных машиностроительных операциях либо обнару­ живается, что м. о. zB настолько мало, а ст2В так близко к сг2б, что искажением можно пренебречь, либо обнаруживается, что достаточно применить упомянутую коррекцию и принять агР равным о2б, переходя тем самым к условиям независимой на­ стройки (см. п. 4.2). В частности, во всех случаях, когда над интервалом ошибки выборочной оценки распределение F (v) близко к нормальному и этот интервал не достигает точки v = 0,

распределение P+ (z), в

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

а а2р < а2б. Приведем

некоторые данные для нормального рас­

пределения F (v), с г = 1 , на основании которых можно составить представление о порядке величин, характеризующих искажение

Р+ (z) сравнительно с б (г).

Параметры распределений вычислены

на основании соотношения

(4.9):

 

 

Y+

х

аг6

М. О. Z

azP+

3

1

1

1,75

0,73

2

2

0,5

0,49

0,46

1,2

2

0,5

0,35

0,46

0,8

3

0,33

0,14

0,32

Все распределения практически не отличаются от гауссова. Напомним, что единицей измерения является ах — среднее ква­ дратическое отклонение мгновенного распределения признака качества.

Еще раз напомним, что в случаях, когда, отличиями |3+ (z) и Р“ (z) от б (z) можно пренебречь, а также при коррекции вы­

борочных оценок и забракованных отклонений v на м. о. z, при вычислении ф (vHC) для настройки уточнениями можно пользо­ ваться соотношениями (4.4) и (4.6).

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

1С0

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

4.6. НАСТРОЙКА С ДОПОЛНИТЕЛЬНОЙ ПРОВЕРКОЙ

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

Настройки с дополнительной проверкой рассмотрены здесь применительно к случаю, когда при регулировках рабочий не пользуется контрольными данными (независимая настройка) или при настройке уточнениями перераспределением ошибки z вы­ борочной оценки, рассмотренным в п. 4.3, можно пренебречь. Распределение ошибки настройки с дополнительной проверкой, обозначаемое фдоп (инс), можно вычислить по той же схеме, что и распределение ф (инс) со следующими изменениями исходных данных в (4.4). Вместо р (ирг) вводится распределение ip (пнс) ошибки той настройки, которую выполняет рабочий, вычислен­

ное в соответствии

с (4.4). Вместо оперативной

характеристики

L (v)

в

формуле

(4.4) вводится оперативная

характеристика

Anon

(v)

плана дополнительной проверки. Таким

образом, осно­

ванием для вычисления ярдоп (унс) служит следующее

отношение:

Ф допКс) — mT f a c ) W » " .c)------.

(4.34)

J Ф (унс) Ьдоп (°нс) dvac

со

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

v2 = -р-;

Р = { ф Ы Едоп(унс) dvHc.

(4.35)

 

— СО

 

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

101

Г л а в а 5

ПЕРЕРАСПРЕДЕЛЕНИЕ ОТКЛОНЕНИЙ У. Н. v

ВРЕЗУЛЬТАТЕ КОНТРОЛЬНЫХ ПРОВЕРОК

ИВЕРОЯТНОСТЬ БРАКА

5.1.ПРЕДВАРИТЕЛЬНЫЕ ЗАМЕЧАНИЯ

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

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

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

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

после ^каждой настройки,

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

в этой

связи, используется схема

настройки с дополнительной

провер-

102

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

м.

о. Тг = м. о. Tj = ^ - \

Т2 = т 3= . . . = T j _ lt

где Тj, j = 1, 2, . . .,

/ — длительность межпроверочного про­

межутка. Число проверок в течение технологического промежутка равно / —1, иначе говоря, они выполняются только в начале межпроверочных промежутков, начиная со второго (случаи с отклонениями от такого порядка описаны в вариантах приме­ ров в гл. 7).

Вычисление распределения a (ивх) входных отклонений у. н. vBX

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

1) число

проверок

/ — 1 невелико, настройка

износостойка;

2)

то

же,

но износ

настройки неустраним; 3)

число проверок

/ —-1

на­

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

вероятности брака q удобней изложить совместно. Так это и сделано в данной главе.

5.2. ОПЕРАЦИИ С ИЗНОСОСТОЙКОЙ НАСТРОЙКОЙ

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

промежутка и единственное распределение а (ивх)

входного откло­

нения у. н.

ивх,

совпадающее с распределением

i|: (vnc) ошибки

настройки vllc.

 

 

 

 

 

 

 

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

от технических

границ

Г

и /+,

в которых должны

оставаться

отклонения

и =

х — X

 

 

 

 

 

—-----признака качества х во изоежание орака.

Технические

границы

выражены

в статистических

единицах:

 

 

 

-

X

/+ =

■'■6+' X

 

 

 

 

 

 

 

 

 

 

где А'б-, л,'б+ — нижняя и верхняя границы поля допуска. В этих обозначениях вероятность b (и) брака при отклонении у. н. v равна

b (и) = Ф (/ — v) -f Ф (v — /+),

(5.1)

103

где Ф (t) — функция гауссова распределения (см. табл. 1 при­ ложения).

Пример вычисления b (и) _в соответствии с (5.1) приведен

в табл. 7. Вероятность брака q за технологический промежуток в рассматриваемом случае равна

СО

Я= J Ш .сЬ) («нс) <fo|.c

СО

где ф (ц,1С) вычисляется тем или иным способом из числа изло­ женных в гл. 4; алгоритм практических вычислений и примеры содержатся в гл. 6 и 7.

Таблица 7

Вычисление вероятности брака b (и) 1~ = 3,5710; /+ = + 3,5710

i

vi

= l-- v t

 

4 =

 

P2=

b (v) =

 

 

 

= Vt - l+

 

= Ф('Г)

= Pi +Pi

I

2

3

 

4

5

6

7

 

 

 

 

 

 

 

—12

—4,2862

0,7142

<

— 4

0,7611

—о

0,7611

— 11

—3,9281

0,3571

<

4

0,6406

—0

0,6406

— 10

—3,5710

0,0000

< —4

0,5000

—о

0,5000

— 2

—6,7*142

—2,8568

< —- 4

0,0021

- o ' ’

ojo0 2 1

—I

—0,3571

—3,2139

 

—3,9281

0,0007

—о

0,0007

0

0,0000

—3,5710

 

—3,5710

0,0002

0,0002

0,0004

1

0,3571

—3,9281

 

—3,2139

~o

0,0007

0,0007

2

0,7142

< —4

 

—2,8568

~o

0,0021

0,0021

l'o"

3,5710

< —4

 

'0,0000 - o '

6,5000

0*5000

11

3,9281

< —4

 

0,3571

~0

0,6406

0,6406

12

4,2862

< —4

 

0,7142

~o

0,7611

0,7611

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

а(1)(«вх) = Ф(«нС)-

Кконцу первого межпроверочного промежутка распределе­ ние cot1) (ивых) выходного отклонения иВЬ1х остается таким же, каким было at1) (увх).

Вычислим элемент вероятности a<2) (ивх) dvBX для фиксиро­

ванного значения vlx входного отклонения ивх к началу второго промежутка для случая, когда настройка независима. Значе­

ние гфх может, во-первых, остаться от предыдущего межпрове-

104

рочного промежутка,

чему

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

элемент

вероятности

® (^вых) L (^вых) ^^вых1

где

Пвых = ^вх>

А (^вых)

оператив­

ная характеристика. Во-вторых, может оказаться, что vDX = увых является результатом настройки, выполненной после того, как какое-либо из значений увых было забраковано при выборочной проверке. Этой возможности соответствует элемент вероятности

00

1J" ® (^вых) L (Увых) ^вых

00

Обозначив

 

 

 

со

 

 

 

 

 

 

 

 

 

 

Q (l)=

1

— СОJ

(»вых) L (у вых) Л'вых.

(5 -4 )

можно

записать *

при

 

vBX =

ивых = унс

 

 

а<2) (Овх) dvvx=

 

(“

(1) (»вых) L К ы х ) + Q (2)^ (Онс)) dv■

(5 -5 )

Для третьего и последующих межпроверочных промежутков

 

а 0) Кх) dvnx =

(®(/_1) (»вых) 1 Кых) + Q<;) Ф (о„с)1 dv.

(5.6)

Вероятность брака qU'* в течение межпроверочного проме­

жутка

/ равна

 

 

 

 

 

 

 

 

 

 

 

СО

 

 

 

 

д(п =

|

Ъ(и) gc(/) (v) cfo,

(5.7)

где b (v) определяется в соответствии с (5.2).

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

 

;<о

= j —1

;(/>

-jj )

 

+ 2

4 = J

 

(5.8)

 

 

 

 

/ = 2

 

 

Алгоритм и числовой пример вычисления q при независимой настройке для дискретной последовательности значений откло­ нения увх показаны в табл. 11.

Если настройка выполняется методом уточнений, распределе­

ние

ф (уис) зависит от распределения со (ивЫх) выходных отклоне­

ний

предыдущего

межпроверочного

промежутка и вычисляется

в соответствии с

изложенным в п.

4.3.

* Еще раз напомним, что подстрочные индексы upr, oHC, vBX, ивых указывают на ту схему перераспределения вероятностей, в результате которой получила очередное значение одна и та же переменная v непосредственно перед тем, как она становится входной величиной для какой-либо новой схемы. Эта условность существенно упрощает изложение, особенно в случаях, когда v становится аргу­ ментом функции, к которой распределение о не имеет отношения, но которая получает смысл (в рамках всей модели) только после определенного перераспре­ деления v, например L (овых), b (vBX) и др. Кроме того, мнемоническая симво­ лика облегчает ориентировку в многочисленных обозначениях распределений.

105

5.3. ОПЕРАЦИИ С НЕУСТРАНИМЫМ ИЗНОСОМ НАСТРОЙКИ

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

Схема вычисления распределений а </> (ивх), / = 1 , 2 , . . . , . / входных отклонений у. н. ицх в случае линейной динамики v в течение межпроверочного промежутка отличается от случая заведомо износостойкой настройки только тем, что в течение

промежутка Т отклонение у. н. v увеличивается на

w — 0х

где а г — угловой коэффициент уравнения линейной

статистиче­

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

Длительность первого межпроверочного промежутка в сред-

Т

нем равна у , следовательно, приращение отклонения у. н. v

соответственно равно у . С другой стороны, приращение w оди­

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

Вычисление распределения а (2> (vBX) входных отклонений при / = 2 выполняется в соответствии с формулой, которая фор­ мально не отличается от (5.5), но в ней распределение cot2) (уВЬ|Х) выходных отклонений к концу второго межпроверочного про­ межутка определяется соотношением

ш<2, (« вы х)== “

< 1,(0 вх — “ »);

 

в общем виде

 

 

ш</)( 0 = « (/

‘’ («вх — “»)•

(5.9)

Распределение сс</> (ипх) вычисляется в соответствии с форму­ лой, подобной (5.6).

Для вычисления вероятности брака b (v) в условиях линейной динамики отклонения у. н. v алгоритм (5.2) непригоден. Здесь

106

dbt/i

о

-4

чт

0,06

q 01- вез предваривче'лььой про­

 

 

вер/ш уровня нас троика

 

0.05

 

 

\

 

0,0k

1 ai=0.00tп

 

/

 

Г ц^бероятн ocmi брака при П/оеде 7Ь-

 

|

НЫХ COLтно цениях 0lim(l’вх)

 

а, =0,0001-

 

 

 

 

 

<а,=0,0001

 

 

 

P T R ?

 

0.01

 

qO^-npu дополнительной

 

проверке настройки

 

 

 

 

 

1 2 3

k £

6 7 8

9 10 j

 

Номера межпроверочных промежутков

Рис.

8. Изменения вероятности брака

в по­

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

при линейном

износе

настроенных

элементов

 

технологической системы

 

Рис. 7. Изменения плотности распределения ве­

роятностей aSfi (ивх) входных отклонений у. н. vax в последовательности межпроверочных про­ межутков при линейном износе настроенных эле­ ментов технологической системы (ах = 0,0001),

можно пользоваться очень удобным алгоритмом Джибра [41 ], переписав его для рассматриваемой задачи в следующем виде:

b (°) = 1 + 4 " — w) ф — w) + cp (/„ — w) /0Ф (/„) — ф (/„)]; (5.10)

10— Г v

при

0;

lv = l+ ■v

при

w > О,

где , /+ — технические (в статистических единицах) границы слева и справа; Ф (/) и ф (t) — функция гауссова распределения вероятностей и соответствующая ей плотность.

Вычисление вероятности брака qUi и q в течение межпрове­ рочного промежутка / и в среднем в течение технологического промежутка выполняется в соответствии с формулами'(5.7) и (5.8).

На рис. 7 показаны кривые а<0 (vBX) в последовательности межпроверочных промежутков, на рис. 8 — вероятности бра­

ка qh'K

 

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

характерны затухание колеба­

ний около постоянного уровня

и убывание промежутка между

точками его пересечения1. К этим важным особенностям мы вернемся в п 5.6.

5.4. ПРЕДЕЛЬНОЕ РАСПРЕДЕЛЕНИЕ ВХОДНЫХ ОТКЛОНЕНИЙ У. Н. цвх

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

Вэтих условиях для вычисления средней вероятности брака q

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

1 Фактически вычисленная

последовательность

приблизительно сов­

падает при целочисленных х с функцией типа у =

АеГ а

sin ф (х),

ф (х) —

монотонно возрастает.

В частности, на рис.

8 кривая

с постоянным

уровнем

0,026, начиная с / =

2, довольно точно соответствует функции

 

у = АеГ а {x~ b) sin 0—.V R W x2);

2 ^ x ^ R

при

А =

0,0125, а = 0,27,

 

6 = 2 ,

Д0 =

53,4;

Д = 49.

 

 

 

108

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

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

последовательностью

значений vc, i

= imln, . . ., —2,

—1, О,

1 , 2 , . . . , гтах (см. п.

3.2). Это значит,

что все состояния техноло-

 

 

(

h

гическои системы, которым соответствуют ивх в интервале (

vt— у ,

in -j- y j рассматриваются как одно и то же состояние, которому

соответствует vt.

Состояния системы, которым соответствуют дискретные зна­

чения vh обозначим через А\п , т = 1, 2, . . ., М, где / — номер межпроверочиого промежутка, т — номер в последовательности возрастания vt. Иначе говоря, взаимосоответствие состояний Ат и отклонений у. н. vt можно представить следующим образом:

vi : V w i -

*b

. . . .

V,

;

 

1 *т а х ’

Ат :Ai Ао,

Л_(<т1п.И)/1-1т1пЛ_/т1п+1,

■ ’

‘ т 1 п + ‘ тах=М -

 

 

В дальнейшем отклонение у. н. vBX = vt, соответствующее состоянию Ат , обозначено через vt (т) *. Рассмотрим условную вероятность того, что к началу /-го межпроверочного промежутка

возникает состояние системы Л5<;), если к началу (/—1)-го про­

межутка процесс был в состоянии

Здесь s и г для упроще­

ния дальнейших записей заменили

индекс т , причем г = т

в обозначении АЦ~1} состояния технологической системы к на­ чалу предыдущего межпроверочного промежутка, a s = т к на­ чалу данного промежутка.

Если рассматривается операция с заведомо износостойкой настройкой, вероятность prs перехода технологической системы

из

состояния А /( ~ 1) в состояние Л5(/)

определяется

следующими

соотношениями:

 

 

 

 

 

 

 

 

 

Р г г =

L (»* (г)) + [ 1

L (о, (г))] ф (vi (г));

|

 

{ j ^

 

Prs =

I1 — L (v (г))] ф (о, (s))

при /'

s,

(

 

 

где

L (vt (т)) — оперативная

характеристика

при

v =

vL (m);

О

 

 

 

настройки

при

 

vIIC =

vi (т).

ф (v{ (m)) — вероятность ошибки

 

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

109

Соседние файлы в папке книги из ГПНТБ