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

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

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

крутизне X. Эта величина зависит от /гд и называется коэффициен­ том статистической полезности одного наблюдения (статистиче­ ской полезностью):

Я = яд : яг.

(3.38)

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

Планы Г.2 и Г.З опираются на члены вариационного ряда, совпадающие с центром или близкие к нему. В отличие от них планы Г.1 (метод крайних значений) и, в известной мере, планы Г.4 (комбинация методов медианы и крайних значений) опи­ раются на первый и последний члены вариационного ряда, распре­ деления которых резко асимметричны и эксцессивпы. Аппрокси­ мация соответствующих частных оперативных характеристик функцией нормального распределения характеризуется более заметными погрешностями— до 0,02. Однако, как показали рас­ четы, для тестовых (рассчитанных на худшие результаты) ситуа­ ций влияние этих погрешностей на показатель эффективности СРК так мало, что не может повлиять на правильность выводов при выборе оптимального варианта.

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

Способ крайних значений Г.1 обсчитан применительно к слу­ чаю, когда измерение выполняется предельными калибрами, соответствующими допуску, т. е. так, как это обычно делает ра­ бочий. Таким образом, границы интервалов группировки в при­ мере соответствуют техническим границам 1~ и /+. От середины поля допуска, границы регулирования удалены на —0,5, 0,5

допуска

6 6 аА. При допуске б = 6 оА. границы интервалов груп­

пировки равны «гЛ = —3, «f2)i = 3.

Объем выборки п — 1 0

взят на

таком уровне, чтобы получить

параметры оперативной

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

Если сравнить способ крайних значений Г.1 со способом медианы Г.2 оказывается, что во втором случае можно обойтись вдвое меньшим объемом п выборки, получив даже более благо­ приятную оперативную характеристику с параметрами кру­ тизны X = 1,87 вместо X = 1,75 при планах Г. 1 (см. гр. 5). Ста­ тистическая полезность Н одного наблюдения при планах Г.1 равна 0,31, между тем, для планов Г,2 Н = 0,70 (см. гр.9). Иначе

89

говоря, при плане Г.1 надо выполнять в 2,3 раза (0,70:0,31) больше измерений, чем при плане Г.2. Сравнительно с планом А при плане Г.1 требуется в 3,4 раза (1,00 : 0,31) больше измере­ ний. Имея в виду, что план Г.1 является по сути отказом от ста­ тистической обработки данных, приведенные соотношения пред­ ставляют собой любопытную иллюстрацию полезности такой обработки. Надо добавить, что с возрастанием объема выборки

при планах Г.1 полезность

одного

наблюдения быстро

падает

и при значительных я = 2 0

способ

крайних значений

вообще

не рекомендуется. Кроме того, как указывалось, оперативная характеристика планов Г.1 заметно (на 0,02) и в неблагоприятную сторону отклоняется от оперативной характеристики планов А на периферии. Если это отклонение при я = 10 несущественно при выборе оптимальных вариантов СРК, то с возрастанием я оно дополнительно обесценивает планы Г.1.

Параметры у~ и у+ оперативной характеристики (см. гр. 6 и 7) при планах Г.1, Г.2 и А можно легко изменить в зависимости от экономической или технологической целесообразности. Тем не менее, полезно сравнить два часто встречающихся плана: а) план Г.1 с границами регулирования, совпадающими с границами поля допуска (см. табл. 3). При этом варианте можно ограничиться измерениями предельным инструментом — обычными калибрами; б) план А с английскими границами регулирования, удаленными

от границ поля допуска на ^ 3 ---- ^ ах.

В обоих случаях предположим одинаковую крутизну X = 1,75. Оказывается, что при плане А положение у+ оперативной харак­

теристики

справа равно: уд = 3 — ^ 3 ----р ^ -^#«1,7; при плане

Г.1— уф I

= 1,5.

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

То обстоятельство, что изменение параметров у" и у+ для оперативных характеристик планов Г.1 связано с заменой обыч­ ного (по допуску) предельного инструмента инструментом, особо изготовленным, или связано с переходом к шкальному инстру­ менту, сильно обесценивает эти планы сравнительно с А, Г.2, Г.З в тех случаях, когда выбирается оптимальный вариант СРК.

Сравнивая план Г.2 с планом Г.З, видим, что в том случае, когда планы Г.З рассчитаны на фактическое соотношение между допуском б и ах средним квадратическим отклонением (см. гр. 8 ), эти разновидности экономически эквивалентны. Статистическая

полезность Я одного измерения в обоих

случаях равна

0,70,

и оба способа допускают применение как

шкального, так

и спе-

0 В- В. Головинский

 

81

циально изготовленного предельного измерительного инстру­ мента.

Для сравнения планов Г.З при б = бег* с английским вариан­ том вычислим удаление точки равновесия оперативной харак­ теристики от границы поля допуска: (3—yt. з) а* = (3 — 1,31) а* = = 1,69а*. Для английского варианта аналогичная величина

равна

За* — - ^ = За* ^ 1 —

1,5а*,

т.

е.

примерно

столько

же.

 

 

 

 

Все, что до сих пор сказано о способе калибров распределения,

относилось к операциям, на

которых б =

6 а*.

Но

практически

этот способ применяется с теми же интервалами группировки на всех операциях с б > 6 а*. Между тем, в случае, когда б = 6 а*, вследствие неустранимой неточности настройки, вероятность брака всегда превышает пренебрежимый уровень и технологически соотношение б = 6 а* недопустимо. В общем случае минимальный допуск, обеспечивающий пренебрежимо малую вероятность брака в машиностроении равен б = 9а*. Применительно к этому соот­ ношению вычислены параметры оперативной характеристики планов Г.З (вторая строка для плана Г.З в табл. 3). Оказывается,

что при соотношении

б =

9а* статистическая полезность одного

наблюдения при плане Г.З

падает сравнительно со случаем б =

= 6 а* незначительно.

Но удаление точки равновесия от границы

поля допуска равно (4,5—2,1) а* = 2,4а*, в то время как при

английском варианте

За* ^ 1 ----\~jw ) ~ ^ °х ^ — 0.57) =

1,29а*.

При планах Г.З и б =

9 риск излишней настройки ^при а =

-|----3

и q =

0,0013) равен 1

L (1,5) = 0 ,1 4 вместо 0,0013 при англий­

ских

границах.

 

 

Сопоставления с английским вариантом сделаны не потому, что он является эталоном оптимальности, а потому, что планы Г.З опираются на него в якобы нормальном случае, когда б = 6 а*- Возможны, например, операции с очень дешевой настройкой и дорогим браком, при которых в конечном счете выгодно согла­

ситься

с очень высокой вероятностью лишней настройки,

имея

в виду

одновременное снижение риска нарушить допуск

(как

это получается при планах Г.З и б = 9а*). Существенно то, что экономическая эффективность всех разновидностей планов класса Г, при которых границы интервалов группировки определяются не технико-экономическими особенностями операции, а зависят только от допуска, является делом случая, а не результатом пла­ номерного выбора оптимального варианта.

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

82

Сопоставление комбинированного способа медиан и крайних значений (планов Г.4) со способом медианы Г.2 привело к не­ сколько неожиданному результату. Казалось бы, усложнение плана добавлением пары границ должно было повысить его эффек­

тивность.

В действительности случалось наоборот— статисти­

ческая полезность

Я

одного наблюдения упала с Я = 0,70 до

Я = 0,58.

Впрочем,

в

этом нет ничего удивительного, если при­

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

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

теристиках,

заключается в

отыскании оптимальных параме­

тров X, у~ ,

у+.

Напомним, то, кроме этих величин, показатель

эффективности

зависит еще

от сроков выборочных проверок и

от заданного уровня настройки 86 и др.

 

 

3.7.

экономический смысл

 

 

 

ПАРАМЕТРОВ X, у~, у+

 

 

ОПЕРАТИВНОЙ ХАРАКТЕРИСТИКИ

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

Начнем с параметров положения у~ и у+. Оба эти параметра замечательны тем, что их уменьшение всегда приводит к увели­ чению числа настроек, что иногда оправдывается, а иногда не оправдывается уменьшением доли брака. В дальнейшем будем говорить только о положении оперативной характеристики справа у+, имея в виду, что все сказанное относится к положению оперативной характеристики слева у~ . Итак, число настроек при уменьшении у+ увеличивается, но при этом:

если у+ уменьшается до

----- 3 ^ 1 -----

то доля брака

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

6*

83

технологически оправданных настроек, т. е. таких, которые делаются в случае вероятности брака, превышающей пренебре­ жимо малый уровень, приравненный 0,0013. Будет ли увеличение числа настроек при этом оправдано экономически, зависит от потерь из-за брака и'затрат на настройки;

если у+

уменьшается

от у+ =

---- 3^1 -----^-) до у+ = —----

— 3 ^ 1

, то доля

брака

продолжает уменьшаться, что

по-прежнему сопровождается увеличением числа настроек. Однако число настроек растет не только за счет технологически оправданных, но и за счет технологически неоправданных на­ строек. Экономически уменьшение у+ на этом участке выгодно до тех пор, пока приращение экономии на браке — АУ превышает

приращение AR затрат н потерь на

настройках

и связанных

с ними измерениях А/С;

 

настроек особенно

число технологически неоправданных

быстро растет при у+< - | - , причем

при

таком

соотношении

возможны случаи, когда план проверки приводит к ошибке по знаку. Иначе говоря, иногда может возникнуть ложное пред­ ставление, что уровень настройки v надо уменьшить, хотя в дей­ ствительности его следовало бы увеличить. При всем том, если стоимость 1 ®6 брака относительно высока, а настройки дешевы,

значения у+< - ^ - бывают экономически оправданы.

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

Перейдем теперь к параметру крутизны X. Увеличение этого параметра всегда приводит к уменьшению доли брака и к умень­ шению числа технологически неоправданных настроек. Число технологически оправданных настроек, наоборот, возрастает. Увеличение X всегда приводит к возрастанию затрат на измерения и статистическую обработку данных (затраты /Сх) и, как уже сказано, затрат на настройки R. Увеличение X выгодно до тех пор, пока — АУ >• AR + АКг. Ясно, что при больших потерях в расчете на 1 % брака и небольших затратах на измерения кру­ тизна X должна быть относительно высокой. Более сложные

соотношения, связанные с

крутизной X, рассмотрены в связи

с вопросами оптимизации в

гл. 8 и 9.

Г л а в а 4

ВЕРОЯТНОСТНЫЕ СХЕМЫ И АЛГОРИТМЫ, СВЯЗАННЫЕ С СОБСТВЕННО НАСТРОЙКОЙ

4.1. ВСТУПИТЕЛЬНЫЕ ЗАМЕЧАНИЯ

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

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

рабочего

с целью обеспечить соответствие фактического значения

X мате­

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

шении повторения

попытки.

 

 

 

В дальнейшем

отклонения

у.

и. v, возникшие в

результате

настройки, обозначаются иис

и

именуются ошибкой

настройки

 

 

х т — ее

(4.1)

 

 

 

 

где Х нс •— м. о. х

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

 

Плотность распределения ошибки настройки оис обозначается ф (онс). Отклонения у. н. v, возникшие в результате отдельной регулировки, обозначаются vpr и именуются ошибкой регули­ ровки

V

=

(4-2)

где А рг — м. о. х после

регулировки.

обозначается

Плотность распределения

вероятностей ошибки и

Р К г )' Различаются независимые и уточняющие настройки. Первые

отличаются тем, что ошибка регулировки vpr при

них исчерпы­

вается технической погрешностью г/тех, в то

время

как при уточ­

няющих настройках погрешность ирг равна

 

 

ирг ^ У т е х + z

 

(4.3)

85

//тех — техническая погрешность регулировки, в результате ко­ торой вместо намеченного изменения d в отклонение у. н. v фак­ тически внесено изменение, равное d1\

_

х _gg

z — ошибка выборочной оценки и =

-------- отклонения у. н. и,

 

а х

полученной в ходе той проверки, в результате которой было забраковано и:

z = и V.

 

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

утех и z обозначаются

11 (Утех)> Р (^ )-

суммой: погрешности,

Техническая погрешность //тех является

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

Погрешность информации z совпадает с погрешностью выбо­

рочной оценки и, полученной при проверке забракованного откло­ нения у. и. V.

В зависимости от технологического содержания регулировок возможны разнообразные формы распределений р (vpr) ошибки регулировки и, в конечном итоге, распределений ошибки на­ стройки ф (инс). По этому поводу интересные данные приведены в книге Соломина [24]. Вероятностные схемы, рассмотренные в данной главе, пригодны для любых форм распределения р (ирг). Но в числовых примерах всюду для него предложена гауссова форма распределения, причем не только во избежание излишних (для иллюстративности) вычислительных осложнений, но и по­ тому, что она чаще других встречается при технологически пра­ вильных способах настройки.

4.2. НЕЗАВИСИМАЯ НАСТРОЙКА

Для вычисления показателя эффективности варианта СРК

5 (со) всегда необходимы данные о распределении ошибок на­ стройки ф (унс). При независимых настройках результат каждой очередной попытки не зависит от числа предшествующих неудач­ ных попыток и от значений допущенных ранее ошибок регули­ ровки. Для того чтобы ошибка настройки оказалась равной v'HC, необходимо, чтобы при условии удачной попытки реализовалось значение vpr = v'nc, вероятность чего равна р (ирг) L (ирг). По­

86

этому, воспользовавшись известным соотношением [6 , с. 54], опирающимся на формулу Байеса, можно записать

Ф (Рнс) =

Р(црг) L (^рг)

(4.4)

СО

 

J Р Фрг) L (ирг) & р г

 

здесь р (ирг) = т) (г/тех);

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

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

плана проверки крг, а знаменатель правой части равен вероят­

ности удачной попытки (независимо

от значения ипс).

Напомним,

что

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

вероятностей

технической

ошибки г/тех

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

мерностей, составляющих ос­

 

 

 

нову управления уровнем на­

 

 

 

стройки (см. гл. 1). Факти­

 

 

 

ческие данные о распределе­

 

 

 

нии 1] (г/тех) можно получить

 

 

 

расчетным путем или экспе­

 

 

 

риментально

на основании

 

 

 

обработки контрольных карт.

 

 

 

Соотношение (4.4) в неко­

 

 

 

тором смысле отражает влия­

 

 

 

ние на точность настройки

 

 

 

выборочных проверок и при­

 

 

 

нимаемых на этой основе ре­

 

 

 

шений. Технический

фактор

 

 

 

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

Рис. 4.

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

пределения

г]

(г/тех)

техни­

ностей т| ((/тех) и ф (кнс)

технической по­

ческой ошибки г/гех. Роль

грешности (/тех и ошибки

настройки инс

проверок и решений

харак­

 

при А, = У~7

теризуется

различием между

 

 

 

т] (г/тех) и ф (и,1С). На рис. 4 наглядно показано, насколько большое влияние на точность настройки и на вероятность брака могут ока­ зать проверки и решения в довольно обычном случае, когда ау = сх,.- Практически распределение ошибки настройки ф (н„с) вы­ числяется для последовательности дискретных значений * vit

* Для того чтобы впредь не возвращаться к этому вопросу, условимся, что далее в алгоритмах встречается следующая особенность обозначении, пока­ занная ниже применительно к ошибке настройки ннс, но относящаяся и к (/тех, ирг. ^вх. ивых и т. д.: ф(инс). ф(цНс)—плотность распределения вероятностей непре­

рывной случайной величины; ф (о,-); у; =

+ h— последовательность вероят­

ностей того,

что значение величины онс окажется в интервале округления ot- —

h

h

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

---- sg

Vi + -g-- Эта функция называется

меняет плотность распределения при дискретизации. При фиксированном / = С

 

2

(4.5)

Ф М

J h_

Ф Кс) dvuo ■

 

V..

г

 

Г

 

87

i

i'min, • ■

-I 2 , 1 ) 0 , 1 , 2 ,

. . imax, у,-+1

vi

h, v0

0

с отсчетом

от дВ — заданного

уровня настройки.

Алгоритм

вы-

О

числения вероятностей ф (vt) в краткой записи сводится к сле­ дующему.

 

Алгоритм вычисления вероятностей ошибки Vi при независимой

 

настройке

I.

И с х о д н ы е д а н н ы е :

1.

Последовательность допустимых значений отклонения у. и. и,-.

О

2. Вероятность т| ((/;) технической ошибки f/;.

3. Оперативная характеристика L (о) плана выборочной проверки отклоне­ ний у. и. орг после каждой регулировки.

II. А л г о р и т м :

4 м -

(4.6)

 

2 j Пi‘Ji) 1 (vi)

Вразвернутом виде тот же алгоритм (для допустимых i):

1.а ,• = i°i (!/«) L (Vi).

2. Р = ]£а/. i

3. Ф (vi) = ^ .

Числовой пример реализации алгоритма приведен в табл. 5.

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

 

 

/“

00

 

vi

=

£

п р о - р у - 1) = 1 + р п _ ‘(7 - р )Т = ~ т *

(4-7)

где Р =

J

p(vpr)L (v pr)dvpr.

 

 

—оо

 

 

Расчеты этой величины приведены в гр. 5 табл. 5.

Формула (4.6) пригодна также и для уточняющих настроек, но при условии, что информация об отклонении у. и. v перед каждой регулировкой является результатом выборки, выполнен­ ной только для того, чтобы определить d — необходимое уточ­ нение отклонения у. н. v. В этом случае ошибка регулировки нрг рассматривается как сумма двух независимых случайных вели­ чин г/тех и z, причем распределение г совпадает с распределением ошибки выборочной оценки отклонения у. н. v. Например, если

88

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

 

 

 

 

 

о

i

vi

л (</,•)

L (ui)

о

=

 

 

 

=

M 0 £ ) 4 ui)

I

2

3

4

 

5

— 12

—4,2862

0,0010

0,2389

0,0002

— 11

—3,9281

0,0021

0,3594

0,0008

— 10

—3,5710

0,0042

0,5000

0,0021

—9

—3,2139

0,0072

0,6406

0,0046

—2

—0,Vi42

0,0885

0,9979

0,0883

— 1

—0,3571

0,0979

0,9993

0,0979

0

0,0000

0,1006

1,0000

0,1006

1

0,3571

0,0979

0,9993

0,0979

 

 

 

Р =

0,9596

 

 

 

 

 

1

 

 

 

Vi =

P- =

 

 

 

 

=

1,0421

Таблица 5

 

£

II

II

o le .

1

 

 

 

6

 

0,0002

0,0008

0,0022

0,0048

0,0920

0,1019

0,1048

0,1019

Утех н 2 распределены по нормальному закону (как это обычно бывает) и средние квадратические отклонения их соответственно

равны оу и az, то среднее квадратическое отклонение av

ошибки

регулировки

ирг равно

 

 

% . = V 4 + о!.

(4.8)

4.3.

ИСКАЖЕНИЕ ВЫБОРОЧНОЙ ИНФОРМАЦИИ

 

 

ОБ ОТКЛОНЕНИИ У. Н. v ПРИ

 

 

УТОЧНЯЮЩЕЙ НАСТРОЙКЕ

 

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

89

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