(68)

(69)

Число степеней свободы дисперсии адекватности равно разности чисел степеней свободы дисперсии остаточной и воспроизводимости.

1. ОПТИМИЗАЦИЯ ИССЛЕДУЕМЫХ ПРОЦЕССОВ

МЕТОДЫ ОПТИМИЗАЦИИ

Одним из первых этапов при оптимизации технологических процессов является определение критерия оптимизации – функции отклика Y, значение которой будет положено в основу оценки процесса при его оптимизации. Задача оптимизациии сводится к нахождению таких условий проведения технологического процесса, при которых критерий оптимизации достигает экстремума.

При исследовании технологических процессов аналитическая зависимость Y(X) (где Х — вектор управляемых, а следовательно, контролируемых факторов X=(X₁, Х₂.....Х_к) неизвестна, и исследователь не может найти экстремум путем решения системы дифференциальных уравнений

ду(Х)/дХ_j=0, j=l, 2...k, (72)

где k — число факторов.

Обычно исследователь может лишь наблюдать значения выходной величины при различных комбинациях варьируемых факторов. Принято называть геометрическое изображение функции отклика в факторном пространстве (X₁, Х₂.....Х_к) — поверхностью отклика.

При поиске экстремальной точки, в отличие от аналитического исследования, осуществляется локальное изучение поверхности отклика по результатам ряда опытов, специально поставленных около исходной точки. Движение к экстремуму в n-мерном пространстве независимых переменных осуществляется обычно не непрерывно, а шагами. Анализируя результаты экспериментов и сравнивая их с
результатами предыдущих, исследователь принимает решение о дальнейших действиях по поиску оптимума. Экстремальное значение отклика достигается с помощью многократного последовательного изучения поверхности отклика и продвижения в факторном пространстве. Существуют несколько экспериментальных методов оптимизации: Гаусса-Зайделя, случайного поиска, градиента, Кифера-Вольфовица, крутого восхождения, симплексный, различающихся способом определения направления движения и организацией самого движения.

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

Р ис.13. Поиск экстремумов функции отклика методами Гаусса-Зайделя (а), случайного поиска (б), градиента (в), Кифера-Вольфовица (г), крутого восхождения (д) и симплексным (е).

^{Таблица 14*}

^* - b_i – коэффициенты линейной аппроксимации поверхности отклика вблизи исходной точки, k - число факторов.

ОПТИМИЗАЦИЯ ПРИ МНОГОЭКСТРЕМАЛЬНОЙ ПОВЕРХНОСТИ ОТКЛИКА

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

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

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

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

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

ОБОБЩЕННЫЙ ПАРАМЕТР ОПТИМИЗАЦИИ

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

Для перехода к обобщенной функции желательности (D) необходимо преобразовать натуральные значения частных откликов в безразмерную шкалу желательности или предпочтительности (d).

Назначение шкалы желательности - установление соотношения между натуральным значением функции отклика (Y) и значением частной функции желательности (d). Наиболее часто используемый вариант шкалы желательности имеет интервал от нуля до единицы; значение d=0 соответствует абсолютно неприемлемому значению функции отклика, а d=1—самому лучшему его значению. Принятые отметки шкалы желательности приведены ниже (Табл. 15). Уровень отметки шкалы желательности объясняется удобством вычислений.

d =0,63 1— 1/е;

(73)

d =0,37 1/e.

Таблица 15

Преобразования с помощью шкалы позволяют привести любую функцию отклика (Y) к безразмерной величине. Построение шкалы желательности — это способ формализации (не лишенный субъективной оценки) представлений исследователя о важности частных откликов. После составления частных функций строят обобщенную функ
цию желательности:

D = (d₁d₂d₃...d_k)^{1 / k} (74)

Если хотя бы один из d_i=0, то какими бы ни были остальные d, обобщенная функция D=0, более того, D наиболее чувствительна к малым значениям частных функций желательности. Тем самым исключается принятие решений, при которых хотя бы один отклик принимает неудовлетворительное значение. С обобщенной функцией желательности можно выполнять все операции, как с любым откликом системы и, прежде всего, функцию D удобно использовать для оптимизации процесса.

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

1. Преобразование функций отклика в частные функции желательности;

2. Построение обобщенной функции желательности;

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

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

8. ИССЛЕДОВАНИЕ ПРОЦЕССОВ В ПРОИЗВОДСТВЕННЫХ УСЛОВИЯХ.

ТОЧНОСТЬ И СТАБИЛЬНОСТЬ ТЕХНОЛОГИЧЕСКИХ ПРОЦЕССОВ

По характеру вмешательства экспериментатора в исследуемый процесс эксперименты можно разделить на активные и пассивные.

Рассмотренный выше активный эксперимент предусматривает активное вмешательство в исследуемый процесс с изменением его по заранее намеченному плану. Достоинствами такого эксперимента являются:

• повышение эффективности исследования, т.к. происходит сокращение объема экспериментальных исследований;

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

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

• в условиях эксплуатации нет возможности широко варьировать факторы,

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

• нет возможности осуществлять произвольное варьирование факторов.

Достоинством пассивного эксперимента является отсутствие необходимости тратить время и средства на постановку опытов.

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

• неточность работы оборудования,

• погрешности инструмента,

• неточности приспособлений и технологической оснастки,

• неоднородность используемых материалов и сред,

• ошибки операторов при настройке и поддержании режимов работы,

• метрологические ошибки.

Эти погрешности могут быть охарактеризованы известными числовыми характеристикамии законами распределения. На рисунке 14 представлен пример распределения погрешностей выходного параметра х. Здесь x_max – x_min – поле рассеяния, характеризующее случайную составляющую производственных погрешностей,  - величина смещения центра группирования погрешностей х (среднего арифметического значения) относительно номинала М(х), характеризующая систематическую составляющую погрешностей.