Методы оптимизации в металлургии: Методические указания к выполнению курсовой работы, страница 7

Кроме этого, критерием прекращения движение вдоль оси x1 в направлении возрастания (убывания) параметра оптимизации является не выполнения ограничений на значения параметра x1.

5.Достигнутая точка с максимальным значением параметра оптимизации является исходной для организации 2-го рабочего цикла, который осуществляется в направлении оси х2. Последний  организуется аналогичным образом и заканчивается достижением частного оптимума вдоль этой оси. После этого необходимо осуществить 3-й рабочий цикл вдоль оси х3 и т.д.

6.После завершения k-го рабочего цикла вдоль оси хk необходимо перейти к (k+1)-му циклу, который осуществляет движение вдоль оси хk‑1, затем к (k+2)-му циклу вдоль оси хk-2 и т. д. Поиск прекращается в точке, любое движение из которой приводит к уменьшению (или увеличению при поиске минимума) значения параметра оптимизации.

Иногда, стремясь ускорить сходимость метода, величину  подбирают так, чтобы при переходе от хk к хk+1 вдоль направления спуска (или подъема) обеспечивалось наибольшее (возможное) возрастание (убывание) целевой функции. Другими словами,  находят из условия минимума функции  одной переменной :

Характерной особенностью метода является необходимость продолжительной стабилизации всех изучаемых переменных, кроме одной по которой происходит движение. Таким образом, в методе покоординатного поиска, движение происходит по ломаной, состоящей из отрезков прямых, параллельных координатным осям (рис.4.).

Надпись: Рисунок 4 - Метод покоординатного поиска

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

Рассмотрим реализацию метода покоординатного поиска для оптимизации механических свойств алюминиевого сплава [3] с вышеизложенными ограничениями на химический состав (2) и на механические свойства сплава (5) .

Переведем уравнение прочности (3) из кодированного вида в натуральный вид с учетом зависимости (4):

В качестве начальной точки возьмем точку с координатами:

Х 1,0 = 5; Х 2,0 = 2, 5; Х 3,0 = 0, 5; Х 4,0 = 0, 5.

Шаг варьирования по каждой переменной выбираем по следующему выражению:

                                 (14)

где d - погрешность вычислений, равная 0, 25.

Таким образом, шаг варьирования по каждой переменной (кремнию, меди, магнию и марганцу) соответственно равен:

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

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

2. Выбираем движение по координатной оси х1 в сторону возрастание переменной х1, т.к. sв,1 > sв,2. Для проверки условий (5) по свойствам вычисляем значения параметров пластичности и твердости в точке Х1, 1. Полученные значения пластичности и твердости (d = 1, 813 %; НВ = 78, 95) удовлетворяют условия (5).

3. Находим следующее значение переменной х1,3:

Подставляем полученные значения в уравнение прочности (9) и получаем следующие значения sв,3 = 18, 502. Для проверки условий (5) по свойствам вычисляем значения параметров пластичности и твердости в точке Х1, 3. Полученные значения пластичности и твердости удовлетворяют условия (5).

Дальнейшее движение вдоль координатной оси х1 далее становится невозможным, учитывая ограничения на химический состав алюминиевого сплава (2).

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