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

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

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

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

X / / S Ui

fr-b-^-Ui)]

• h v t [ut I

Ui-,),

(209)

/ = і

 

і=і

 

 

где первый подналадочный импульс Ui = щ выбирается по априор­

ным характеристикам, а цо — начальное

значение

процесса

Подставляя выражение (209) в формулу

(201),

находим

следую­

щее выражение для

риска

 

 

 

 

 

А , = f х„*РЫ-

я я ( ^ |

ПР^)

/

n

^ - U

- b ^ - f

û (*л. Сл. !Ал. ^Л.)

 

 

+

2=1Я Г ;

({/, I Хі-г,

Ui-ùdQ.

 

(210)

 

Представим выражение

(210)

в следующем

виде

 

 

Dn

= j

"Л Гг (ut

1 (/,_,,

*, _ ,) {

J

Ы

U I

X

^

( * я _ і .

п

 

Q (х„,

и п ,

гп)

 

 

 

 

 

 

ХПР(ч)

ПИх,-

(/, +

vi +

L. + f/,)!

X

 

 

 

І=\

І=\

 

 

 

 

 

 

 

 

 

X r „ ( ( / j L . , t / J 4 ^ -

 

 

(211)

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

Перепишем плотность ô[x, — (Ц + и.* +

у + Ui)]

(i =

1, 2,

n)

в следующем

виде:

 

 

 

 

 

 

 

 

 

8 К , - ( я , - - £ / ,

) ] ( * =

1,

2,

п).

 

(212)

Подставляя выражение

(212)

в формулу

(211),

производя

ин­

тегрирование выражения в фигурных скобках по вектору

£ п и пола­

гая, что подналадчик обладает относительно управления

регуляр­

ной стратегией [146], находим, что оптимальный подналадочный

им­

пульс на п-ом такте определяется из условия минимума по Un

сле­

дующего выражения

 

 

 

 

 

 

 

«„=

ГѴ

^ Ы • П Р(ъ

I

П

P(Xi

I Uh

îx;-,

 

(213)

I J

^

\

 

 

 

 

 

 

 

О ir

il. 1

 

 

 

 

 

 

 

158

где Р (Xi\Ui,

\Xi, Ii) —-условная плотность

случайной величины л\-

 

при фиксированной совокупности (Ui, ц.,-, h)

 

параметров процесса,

относящихся к

тому

 

же і-му изделию.

 

 

Поясним

физический смысл величин, стоящих в формуле

(213)

под знаком интеграла. Функция Р 0 ) —априорная плотность на­ чального значения функциональной случайной составляющей;

функции Р (рі| \ii-i)—плотности вероятности перехода этой со­ ставляющей для двух соседних моментов времени. Эти функции от­

ражают априорные

сведения об изменениях

функциональной

слу-

 

п

 

 

 

чайной составляющей. Выражение вида UP

(xAUi,

іц, /,-) содержит

в себе информацию

об управляемом процессе.

Одновременно

это

выражение содержит информацию о собственно случайной состав­ ляющей (;„ процесса.

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

К изложенному тесно примыкают

некоторые результаты работы

[137]. В одном из разделов этой статьи рассматривается

задача оп­

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

 

хп = г(п, Ѳ) +

С„

( я = 1, 2, . . . , ) ,

(214)

где г (п, Ѳ) —систематическая

составляющая;

 

Ѳ — определяющий ее поведение неизвестный

параметр;

последовательность независимых случайных вели­ чин.

Рассматриваемый алгоритм управления использует в качестве оптимальной оценки Ѳ апостериорное математическое ожидание

этого параметра, найденное по результатам наблюдений хп

= (хі,

. . . , Хп) '•

 

 

Ѳ =

м ( Ѳ | х „ ) ,

(215)

где M — символ математического ожидания.

 

Как следует из результатов

работ [83, 137], оценку вида

(215)

используют и оптимальные по точности алгоритмы подналадки.

§ 18. О П Т И М А Л Ь Н Ы Е ПО ТОЧНОСТИ АЛГОРИТМЫ П О Д Н А Л А Д К И СТАЦИОНАРНЫХ ТЕХНОЛОГИЧЕСКИХ ПРОЦЕССОВ

Формула (195) дает общее представление о технологическом процессе. Случайная составляющая этого процесса образуется как сумма двух последовательностей {и.я } и (п=1, 2, . . . ) . В за­ висимости от статистических характеристик последовательности

159

{цп} технологические процессы разделим на стационарные и не­ стационарные. В данном параграфе построим оптимальные алго­ ритмы подналадки для двух классов стационарных технологических процессов, в следующем — для нестационарных.

I . Управление технологическим процессом с корреляционной

функцией Кц (т) =

ѳ | т 1

• Процесс этого типа является марков­

ским [145] и стационарным

[72].

Для

этого

процесса

переходная

плотность имеет вид

 

 

 

 

 

 

 

 

 

 

 

1-1

 

1

exp;

,

(216)

 

 

 

 

 

где r =

е~%.

Пусть также

 

 

 

 

 

 

 

 

 

 

 

- —

exp

 

 

(217)

 

 

 

 

 

 

2з.і

 

 

 

В силу

взаимной

независимости

случайных величин

(п

= 1,

2, . . . ) , произведение плотностей имеет вид

 

 

 

 

 

 

 

 

 

exp

(X.

mt

 

 

 

 

 

 

 

 

 

2s;

 

 

 

 

 

 

 

s? ( 2 . )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(218)

Подставляя выражения

(216)

— (218)

в формулу

(213),

полу­

чаем

 

 

 

 

00

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

2л-И

 

Г*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

з , ( 2 я )

(1

 

 

 

 

 

 

 

 

 

 

 

 

 

''n

 

 

 

 

У

exp

А

 

 

 

si=l-

(xi — ml u-i — £//)2

X

 

 

2а.2

(=1

 

 

2a?

 

 

 

 

 

^

'

 

 

 

 

 

 

 

 

X

^ „ t / u f t

. . .

<af(i,„.

 

 

 

(219)

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

Производя в формуле (219) интегрирование в соответствии с ме­ тодикой [102] и дифференцируя полученный результат по Un, нахо-

160

дим величину оптимального уровня настройки инструмента перед выпуском п-то изделия:

 

иѣа =

-mn

2s-,

 

 

(220)

где qn и p n — величины,

определяемые с помощью

рекуррентной

процедуры:

 

 

 

 

 

Р П =

А

 

 

 

 

(221)

 

г

р п _ 1 + ß. r + • 1

 

 

 

 

 

 

2 а;2С

 

 

 

 

 

 

1

\ ,

Х п - \

(222)

 

 

 

2з?/>„-і

I

£

 

 

2(1

г*\

<,

 

 

(223)

 

 

 

 

вспомогательная величина, зависящая от статистических характе­ ристик процесса. При п = 1, qx = 0

1

Pi = г -

2 ^

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

В соответствии с выражением (220), величина подналадочного импульса для гс-го изделия может быть выражена через соответст­ вующие уровни настройки:

« : = " ; -

к

- , =

- « -

i (f

-

y

= 4

(224)

Анализ выражения

 

(221)

показывает,

что

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

п} с ростом номера

n

монотонно

возрастает

до

какого-то

преде­

ла. Следовательно, при достаточно большом номере п величина р п

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

Очевидно,

при п-*-оо

 

 

рус =

lim ps =r- lim ps-\.

(225)

Подставляя выражение (225) в формулу (221) и решая квадрат­

ное уравнение относительно

р у с , находим

 

 

+ * +

( 2 2 6 )

1 1—2891

161

где

d=

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

^ - - ' « - о - і г -

( 2 2 7 )

В данной главе систематическая составляющая /„ считается де­ терминированной известной функцией времени. Поэтому ее компен­ сация может быть выполнена обычными методами [38]; в дальней­ шем рассмотрим лишь вопросы, связанные с компенсацией функ­ циональной случайной составляющей. Обозначим через Un* ту часть уровня настройки (227), которая учитывает случайную со­ ставляющую,

 

и'п = - г ^

(* = 1.

2.

• • • ) •

(228)

 

ЛРуС

 

 

 

 

Из формулы

(222) следует, что в установившемся режиме

спра­

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

 

 

 

 

qn

= а0<7„_і + 60 х„_і,

(я =

1,

2, • • • ),

(229)

где а0 и bo — некоторые коэффициенты,

зависящие

от статистиче­

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

управляемого процесса. Вид их

может быть найден из формулы (222).

 

Из выражения (228) следует,

что

 

 

q„-i = -

20'пРус.

 

(230)

Подставляя выражения (230)

и (228)

в формулу

(229), находим

связь между уровнями настройки (компенсирующими функциональ­ ную случайную составляющую) для двух соседних изделий:

 

и*п = a t / L i +

fccn-i,

(231)

где а = а0, о — —

—-.

 

 

 

Коэффициенты а и b связаны со статистическими

характерис­

тиками управляемого процесса соотношениями:

 

а

- ( т - Щ 1

+

і £ г }

< 2 3 2 )

 

' р

^ - — \.

(233)

 

?Рус

 

 

 

162

При пренебрежимо малых износах (maO) формулы (231) — (233) являются основой оптимального способа автоматической под­ наладки технологического процесса. В соответствии с этим спосо­ бом измеряют размеры каждого изделия, а величину и знак уров­

ня настройки инструмента

выбирают

совпадающими с

величиной

и знаком суммы предыдущего уровня

настройки и отклонения раз­

мера предыдущего изделия

от заданного уровня, взятых

с постоян-

 

 

г=0,1

 

 

 

 

0,5

 

 

 

 

0,9

 

 

1,0

5

10

15п

 

О

 

Рис. 63. Графики зависимостей ß „ от

п

 

ными коэффициентами, устанавливаемыми в зависимости от стати­ стических характеристик управляемого процесса 1 .

Конструктивные вопросы реализации данного способа подналад­ ки рассматриваются в гл. X I I .

Подставив соотношение (220) в выражение (219) и проинтегри­ ровав получившееся выражение по п Un-\), определим величи­ ну дисперсии управляемого процесса при оптимальном по точности способе подналадки:

+

( 2 3 4 )

где величина рп определяется по формуле (221). Физический смысл дроби 1/2 рп — вклад, вносимый функциональной составляющей u„. Графики зависимостей Dn* от п при а? —<т£ = ] и различных зна­ чениях г приведены на рис. 63.

В установившемся режиме

 

D* =

о? +

а2

 

~

.

(235)

 

ус

.

V-

 

 

 

 

 

 

 

 

(1 4- d)+

- f ( d - l )

 

Приведем

ряд частных случаев.

 

 

 

1. Пусть

функциональная

составляющая

образует

последова­

тельность независимых

случайных величин (г = 0). Тогда

 

 

ус

С 1

 

 

 

2. Функциональная составляющая цг а является случайной вели­ чиной (г — 1). Тогда

D* = а ? . ус С

Авторское свидетельство № 305342. Бюллетень изобретений, 1971, № 18.

11*

163

3. Дисперсия собственно случайной составляющей пренебрежи­ мо мала (0j = 0). Тогда

D*

- (1 — г2 )о2.

ус

V 1

' [Д.

Рассмотрим пример расчета точности ведения технологического процесса при оптимальном уровне настройки ІІ*п на каждом такте. Процесс обработки внут­ реннего диаметра колец карданных подшипников 804704 на автоматической ли­ нии ГПЗ-1 имеет следующие параметры:

г = е~ѳ

= 0,94; а* = 15 мкм2 ; а? = 1 мкм2 .

Суммарная дисперсия Оѵ

неуправляемого

процесса (при отсутствии трента)

равна

 

 

 

DL = cl+ а £ =

16 мкм2 .

Подставляя значения г, о\ , в выражение (235), находим, что при оп­ тимальном по точности алгоритме подналадки в установившемся режиме

D* c = 3,5 мкм2 .

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

Выражение (235) позволяет оценить «эффект управления» (при компенсированном тренте) при различных статистических харак-

 

 

_

-

 

і

 

r

-

 

- ,

 

 

 

 

 

 

нКѵ

"

i

 

 

i

f

пч

 

 

 

 

 

 

 

 

 

 

 

 

 

:

0.7

 

 

 

 

 

0,5

 

 

 

 

 

 

'

08 —

 

 

 

 

 

;

N

/ 4

^ ^

-

 

 

 

« £ _ J

 

 

 

 

 

 

 

 

0,9

I

 

 

 

 

 

 

 

 

 

 

 

 

 

0.9.5

j

 

 

 

 

 

{ .

 

 

I

 

~ ~

 

г=ьа.

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ß

 

 

5

 

W

tS

 

 

 

 

 

 

Рис. 64. Характеристика

«эффекта

управления»

процессом

 

 

теристиках

процесса

г,

о%,

. На рис. 64 изображены

графики

зависимости

отношения

дисперсий оптимального (по точности)

уп­

равляемого

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

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

дисперсий

сг2 / <j2

случайных

составляющих

при различных

значениях г =

е~в.

Здесь

Dz o? + a*

—- дисперсия

неуправляемого

процесса,

D*c = о? -j

 

 

—дисперсия

управляемого процесса

 

при

опти-

 

2/?у с

 

 

 

 

 

 

 

 

 

 

 

мальном управлении.

 

 

 

 

 

 

 

 

 

Из рис. 64 и формулы (235)

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

управ­

ления

возрастает

с увеличением

г = е ~ ѳ и становится

максималь­

ной при г — 1. Физически это соответствует вырождению

случайного

164

ц.п -2-

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

поддающуюся, естественно, компенсации.

 

 

П. Управление

технологическим процессом с

корреляционной

функцией К^{х) =

ст^е-9!-! cos сот. Корреляционную функцию

тако­

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

обработки

по на­

ружному диаметру колец роликовых подшипников на бесцентровошлифовальном автомате 6С133 [53, 54]. Процесс с корреляционной функцией такого вида является марковской стационарной последо­ вательностью второго порядка [145]. В дискретном случайном про­ цессе второго порядка стационарно связанными оказываются вели­ чины на трех соседних тактах: п, \in-u

Формула для риска Dn принимает в этом случае вид

X П

Ut, lL) • /7Г,<*2.

(236)

1 = 1

i-l

 

Оптимальный уровень настройки на n-ом такте определяется из условия минимума по Un выражения [6]

*n = J ХЫ Р (РіІИо) • ПР ЬЧI \ы-и I*-*) X ,

Щх„, ;

j . „ )

 

 

 

 

 

 

X П

Р (X; I U^ili)

dQ.

 

(237)

В работе [6] показано,

что оптимальный

уровень

настройки

на «-ом такте определяется

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

 

 

f r B = = _ / i m _ i »

( л = 1 , 2,

. . . ) ,

(238)

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

управления

 

 

 

 

Рп =

Pln—7Z

7Z

'

Pln-

 

 

= х ( і + ß) _

= _а_

°2

3 2

<?-* cos « ( 1 — <?_ 2 Ѳ cos 2«)

(239)

 

1 — e"2e cos2<u

 

165

1 — й 2 9 C O S - ы

а2 = а2(1 — ^ C O S 2 « / ) ;

 

 

 

0 2 =

а 2 ( l - ^ ( l - g ~ 2 8 C O s 2 u , )

 

 

 

 

 

 

 

° 2

^

 

 

1 - е - 2 8 COS« в»

 

 

 

 

 

выражаются

через параметры

априорной

и переходных плотностей

 

 

 

РЫ~=

 

 

^

е

х р і

- :

^ ]

;

 

 

 

 

 

 

 

— т = " е х

Р

(— (іА і — РИв)а/2о?};

 

 

 

 

 

 

ot

У 2 я

 

 

 

 

 

 

 

 

 

 

 

 

Ѵ-І-І) =

Ц=г ехр {(.^. ajx£_,

+

Pt*.,-_2)a/2af},

a g = e~e cos <o.

 

Oj i

2jc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К сожалению, автор работы [6] ограничился

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

тимального уровня настройки в форме

(238). Между тем, и для про­

цесса с корреляционной функцией К^(г) —о^е~®'

COSCOT

можно

построить рекуррентный алгоритм

(способ)

оптимальной по точно­

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

процесса.

 

 

 

р у с парамет­

Д л я этого надо найти

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

значение

ра р п

и выразить

через р у с

и уровень

настройки

Un

зеличину qn.

Оказывается,

что при отсутствии

трента

оптимальный уровень на-

 

л

 

 

 

 

 

 

 

 

л

 

 

 

л

 

 

стройки ип*

выражается

через

величины

U*n-u

хп-\,

U*n-2,

-<п

формулой

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и*п

= а£ / л _і +

Ьх„_і

+

cÜ*n-2

+ dxn-2,

(240)

где a,

Ь, с, d-— постоянные

коэффициенты, зависящие от статисти­

 

 

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

управляемого процесса.

Дисперсия отклонений

размеров

изделий

при

оптимальном

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

(238) определяется

равенством

 

 

 

 

 

 

 

 

D* =

а 2

 

+

 

 

 

 

 

 

 

При

П-+ со

 

 

 

 

 

_

 

J

 

 

 

 

 

 

 

 

 

 

 

 

ус

2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

2рус

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

§ 19. О П Т И М А Л Ь Н Ы Й ПО ТОЧНОСТИ АЛГОРИТМ П О Д Н А Л А Д К И НЕСТАЦИОНАРНОГО Т Е Х Н О Л О Г И Ч Е С К О Г О ПРОЦЕССА

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

+ £ A u - + С„+

І п = тп (л - 1, 2, . . . ) , (241)

<=і

 

166

где {ДЦІ}последовательность независимых нормальных величин со статистиками (0; о\).

Процессы с такой функциональной случайной составляющей описаны, например, в работах [12, 83]. Процесс с независимыми при­ ращениями является марковским [72] и нестационарным [145].

Пусть начальное значение р 0

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

по нормальному

закону

1

 

 

 

 

 

ехр

{ - Ж } .

 

(242)

 

 

 

Переходная плотность Р (р.;|ц,г-і)

= Р (ЦІ|Ц,,_І)

имеет

вид

Р{Ѵ, I !*,_,) =

1— е х р -

1 •(Pi—

\H-ï)2

(243)

 

% V2т.

\

 

 

 

Подставляя выражения (242) и (243) в формулу (213) и произ­ водя необходимые вычисления, находим оптимальный уровень на­

стройки

 

 

 

 

 

 

 

 

 

 

(J*

= - m n - ^ ,

 

 

 

(244)

где

 

 

 

 

 

 

 

 

 

Рі

= Р—

1

У = 2,

3,

• );

 

(245)

 

 

 

 

 

 

 

 

 

 

 

2*2

 

 

 

 

 

<7/

=

 

 

•и, 2-1

(t =

2,

3,

. . . ) ; (246)

 

2d?

 

+ ? + 2о2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2о?

 

 

 

 

 

При п =

1 : <7і• =

0, /?і = ß.

 

 

 

Dn

 

 

Подставляя выражение

(244) в формулу

для

и

интегрируя

полученное выражение по £/n -i, Xn-i, находим дисперсию при опти­ мальном по точности управлению

D* = a?-f-

(л =

1, 2, . . . ) .

(247)

При п->оо величина р„ стремится

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

значению

р.,с, которое получается из (245)

в результате

решения

квадратно­

го уравнения

 

 

 

 

 

1

1

=

0

 

Рус

 

 

и равно

ус

(248)

 

167