![]() Главная страница Случайная страница КАТЕГОРИИ: АвтомобилиАстрономияБиологияГеографияДом и садДругие языкиДругоеИнформатикаИсторияКультураЛитератураЛогикаМатематикаМедицинаМеталлургияМеханикаОбразованиеОхрана трудаПедагогикаПолитикаПравоПсихологияРелигияРиторикаСоциологияСпортСтроительствоТехнологияТуризмФизикаФилософияФинансыХимияЧерчениеЭкологияЭкономикаЭлектроника |
Оптимизация исследований технических свойств материалов
Из методов оптимизации исследований в настоящее время наиболее признаны методы планирования экспериментов и конечных элементов. В 50-е гг. XX в. Бокс и Уилсон предложили новое решение старой задачи – отыскание оптимальных условий протекания химических, физических и металлургических процессов. В ней рассматриваются процессы, зависящие от многих факторов, в условиях, когда механизм этих процессов неизвестен. В таком случае естественно прибегать к представлению результатов наблюдений полиномиальной моделью. Метод решения был назван планированием эксперимента и получил признание во всем мире. В России его активно развивал В.В. Налимов [18], большой вклад внес для практического использования метода при решении технических задач В.А. Вознесенский [5, 6], а в области исследований свойств бетонов – Е.Н. Львовский, В.Г. Зазимко и др. [11, 15]. В настоящее время сложилась стройная теория планирования эксперимента с достаточно сложным математическим аппаратом и терминологией. Основанное на построении экономичных планов (полные и дробные факторные планы, ортогональные латинские квадраты и сбалансированные блок-схемы) планирование эксперимента позволяет оптимизировать трудовые, временные и материальные ресурсы на проведение исследований и дать более точные оценки неизвестных параметров регрессии при равном числе измерений. Для освоения этого эффективного приема исследований рассмотрим правила формирования и реализации двухуровневых планов, получивших наибольшее распространение при факторном планировании эксперимента. Оптимальное планирование экспериментов. Целью экспериментальных исследований является поиск закономерности между оптимизируемыми величинами и варьируемыми факторами. При этом все эксперименты делятся на активные и пассивные. Пассивным экспериментом является такой, в котором варьируют одним фактором, сохраняя остальные на постоянном уровне. Результаты такого исследования анализируют путем построения графических зависимостей между оптимизируемыми величинами и варьируемым фактором. Активным экспериментом называют эксперимент, выполняемый по плану, предусматривающему все варианты сочетаний нескольких одновременно варьируемых факторов. Активный эксперимент называют планированным. Полная двухфакторная модель второго порядка имеет вид
и в зависимости от наличия коэффициентов
Рис. 6.1. Поверхности второго порядка и их изолинии
Если полином содержит только линейные члены ( Если полином содержит по одному фактору
Если полином содержит и линейные, и квадратичные эффекты (но Если полином содержит коэффициент
Абсолютное числовое значение коэффициента эффекта взаимодействия Проекция поверхности отклика Чтобы избежать необходимости использования полиномов высокого порядка, был предложен шаговый метод изучения поверхности отклика, напоминающий итерационный метод решения задач вычислительной математики.
Графически такое движение по градиенту в задаче с двумя незави-симыми переменными представлено на рис. 6.2. На рисунке нанесены кривые равного выхода для одного из технологических процессов. Эти кривые аналогичны кривым равной высоты на географических картах. Цикл попеременного движения повторяется много раз (получается своеобразное блуждание по лабиринту), он усложняется с ростом числа независимых переменных (рис. 6.3). Движение из точки О в направлении ОР соответствует наиболее крутому пути подъема по поверхности отклика (отсюда и название метода – метод крутого восхождения). В направлении ОР исследователь будет двигаться до тех пор, пока не попадет в точку Q. В окрестности точки Q нужно поставить вторую серию опытов и заново найти локальное линейное приближение поверхности отклика. Для сравнения на рис. 6.3 показано пунктиром движение по поверхности отклика при однофакторном эксперименте. В этом случае двигаются попеременно, сначала изменяя одну переменную и фиксируя другую, а затем, изменяя вторую переменную и фиксируя первую. Цикл попеременного движения повторяется много раз – получается своеобразное блуждание по лабиринту, который усложняется с ростом числа независимых переменных (рис. 6.3): М1; движение при х1 = сonst от точки М1 до частного максимума М2 и т. д.
Движение по градиенту давно известно в науке. Бокс и Уилсон использовали его в сочетании с дробным факторным экспериментом для локального описания поверхности отклика. Это и определило успех «крутого восхождения».
Рис. 6.4. Крутое восхождение
Для выбранной случайным образом достаточно малой области факторного пространства планируют дробный факторный эксперимент, проводят первую серию опытов (обычно из четырех) и строят линейную функцию отклика. Это еще не поиск экстремального значения функции, а предварительное отыскание направления дальнейшего поиска. Найдя градиенты уравнения (углы наклона поверхности в каждом направлении), повторяем операции и достигаем, наконец, вершины поверхности отклика. В этой области, как уже указывалось, проводят полный факторный эксперимент с определением не только линейных коэффициентов регрессии, но и всех учитываемых взаимодействий. Тот факт, что функция отклика в окрестности исследуемой точки почти не изменяется, еще не говорит, что исследователь находится вблизи точки максимума. Это может быть медленно поднимающийся гребень или гребень постоянной высоты; может быть седловидная точка, являющаяся максимальной по одному направлению и минимальной по другому; может быть точкой локального максимума. Именно поэтому в данном месте необходимо большое количество экспериментов и для решения задачи приходится переходить на рототабельные построения, познакомиться с которыми можно в [16]. Итак, для математического описания технологического явления требуется минимальное количество опытов, из которых информация извлекается с максимальной полнотой, и позволяет получить с заданной надежностью математическую модель процесса в аналитическом или графическом виде. При постановке задачи изучают опытные данные, полученные ранее. На основании этих данных выбирают как оптимизируемые величины, так и варьируемые факторы, которые должны быть независимыми друг от друга. Успех решения экспериментальной задачи в значительной мере зависит от того, насколько удачно выбран план эксперимента. Применяют ортогональные, рототабельные, Д-оптимальные и другие планы. Наиболее часто используют ортогональные симметричные планы, которым присущи минимальное число опытов, простота вычисления коэффициентов математической модели (уравнения). План обычно представляет матрицу, охватывающую изучаемую область. Если в эксперименте с двумя переменными х 1и х 2каждая переменная изменяется на двух уровнях (например при подборе состава бетона расход цемента х 1 изменяется на уровнях 280 и 320 кг/м3, а расход воды х 2 – на уровнях 190 и 210 л/м3), то все возможные комбинации варьируемых таким образом двух факторов будут исчерпаны перебором, приведенным в табл. 6.1.
Для упрощения расчетов от натуральных значений исследуемых величин переходят к кодовым переменным, принимающим на верхнем уровне значение +1, на нижнем –1. Этот переход выполняется по формуле
хi – кодовая переменная; Xi – натуральная переменная; X 0 – натуральное значение величины в центре эксперимента; D Xi – натуральное значение интервала варьирования. Полный факторный эксперимент для двух переменных, варьируемых на двух уровнях (планирование типа 22), приведен в табл. 6.1.
Таблица 6.1 План эксперимента
Табл. 6.1 называется матрицей планирования: каждая строка фиктивной переменной х 0 = + 1, в графах 3, 4 – значения х 1и х 2(собственно планирование); в графе 5 – значение х 1 х 2(эффект взаимодействия факторов в графе 6 – результаты наблюдений. Пользуясь планированием 22, можно определить коэффициенты регрессии неполного уравнения
В это уравнение переменные х 1 и х 2 входят не в абсолютном, а в относительном (кодированном) значении. Графически эту зависимость можно представить поверхностью отклика горизонтальной проекцией которой является квадрат (см. рис. 6.3). Число опытов (четыре) равно числу оцениваемых параметров ( Матрица планирования для 3 переменных – 23 получается из матрицы 22 при повторении её дважды: один раз при значении хi на нижнем уровне, второй раз – на верхнем. На рис. 6.5 показана схема планирования типа 23 с координатами точек куба в кодовом масштабе. Если приблизительно известна область, где находятся экспериментальные значения оптимизируемых величин, эксперименты проводят в этой почти стационарной области.
При изучении почти стационарной области возникает ряд сложных проблем. Если мы хотим описать эту часть поверхности отклика полиномом второго порядка, то переменные нужно варьировать уже на трех уровнях, но планы полного факторного эксперимента типа Зк здесь неприемлемы, так как они потребуют слишком большого числа опытов. Боксом и Уилсоном было предложено построение композиционных планов, ядром которых служат линейные ортогональные планы. Предполагается, что, попав в почти стационарную область, исследователь сначала ставит опыты, используя линейные планы. Затем, убедившись в том, что гипотеза линейности здесь не проходит, он достраивает линейный план до плана второго порядка; отсюда и название – композиционный план. Рассмотрим стратегию построения композиционного плана на примере задачи с тремя переменными (рис. 6.5). Сначала ставим опыты по линейному плану в точках, задаваемых вершинами правильного симплекса (тетраэдра). Эти точки являются подмножеством вершин того куба, которым задаются границы варьирования переменных в линейной задаче. На рис. 6.5 они обозначены зачерненными кружками. Далее ставим эксперимент в центре куба для проверки гипотезы адекватности. Если гипотеза адекватности не проходит, то достраиваются вершины куба (незачерненные точки), и затем добавляется еще часть так называемых «звездных» точек, образующих октаэдр. На рисунке эти точки обозначены звездочками. Здесь видно, что при переходе к плану второго порядка границы варьирования переменных расширяются. В результате получается композиционный план второго порядка, содержащий После постановки опыта вначале выполняется статистический, а затем технологический анализ полинома. Иногда полином исследуют на экстремум (тогда и задачу называют экстремальной). Но чаще, особенно при необходимости изучения нескольких оптимизирующих величин, определяют условный оптимум исходя из компромиссных требований, предъявляемых к оптимизируемым величинам. Такие задачи называют компромиссными. При решении компромиссных задач наиболее важным является технологический анализ полученной математической модели. Он позволяет установить влияние на оптимизируемую величину не только изучаемых факторов, но также эффектов их взаимодействия. В ряде случаев становится возможным достаточно точный прогноз улучшений условий изучаемого процесса. Воспользуемся примером оптимизации процесса тепловой обработки бетона, приведенной в [11], с помощью планированного эксперимента (табл. 6.2, 6.3). Его достоинство в том, что в нем четко прослежена вся цепочка статистического анализа результатов опыта. Таблица 6.2 Исходные данные
В качестве оптимизируемой величины принята прочность бетона в 28-суточном возрасте Rсж (или Варьируемыми факторами являются:
Выберем центр эксперимента Таблица 6.3 Планирование эксперимента
Исследования проводятся сразу почти в стационарной области, т. е. в области оптимума, так как из опыта (априори) эта область известна. Для описания процесса полиномом второго порядка выбираем ортогональный план полного факторного эксперимента 23. Кроме этого, в план введем фиктивную переменную х 0 = 1. Общее число опытов согласно плану N = 15 (табл. 6.3). После выполнения всех 15 экспериментов (каждый опыт имел 6 параллельных определений) были получены средние значения прочности у. Для получения ортогонального планирования второго порядка произведем преобразования квадратичных переменных:
Выберем величину интервала для звездных точек – для случая трех переменных a = 1, 215 (табл. 6.2, 6.3). Тогда
и для основного планирования
а для звездных точек
Эти преобразования необходимы для смещения при определении свободного члена полинома, а также при расчете коэффициентов квадратичных членов
Уравнение, которым описывается процесс, до преобразования имеет вид:
После преобразований
Все коэффициенты этого полинома при ортогональном планировании определяются и оцениваются независимо друг от друга. Это означает, что если тот или иной коэффициент окажется незначим, его можно отбросить, не пересчитывая все остальные. Коэффициенты полинома определяются по формуле
а для определения свободного члена уравнения воспользуемся выражением
Коэффициенты уравнений (6.6) и (6.7) имеют следующие значения
Знаменатель формулы (6.8) для определения коэффициентов
Знаменатель формулы (6.8) для определения коэффициента
+
Знаменатель формулы (6.8) для определения коэффициентов
После подстановки значений коэффициентов уравнение примет вид
Значимость коэффициентов уравнения определим по
где Дисперсию
Дисперсию воспроизводимости
Числитель формулы (6.12) определим из выражения
……………………………………………………………………………………… ……………………………………………………………………………………… Знаменатель формулы (6.12)
Вернемся к определению
Дисперсии для
определим следующим образом: Дисперсию свободного члена представим как
Теперь можно определить
Условие значимости коэффициентов уравнения имеет вид
Адекватность (тождественность) описания процесса термообработки полиномом второй степени проверим, используя
Дисперсия воспроизводимости определена ранее и равна Остаточную дисперсию определим так:
Число степеней свободы где Степень свободы ¦2 определена ранее (¦2 = 75). Определим дисперсию
|