книги / Математическое моделирование и основы научных исследований в сварке. Статистическая обработка и планирование эксперимента-1
.pdfРис. 38. Листинг поиска оптимума: пятый и шестой циклы
171
Рис. 39. Листинг поиска оптимума: седьмой и восьмой циклы
172
Рис. 40. Листинг поиска оптимума: девятый цикл и окончательный расчет с определением оптимального варианта. Перевод кодированных значений факторов в реальные
173
Например, сначала изменяют в цикле от минимального до максимального значения факторы x1 и x2 , а остальные
фиксируют на основном уровне. По полученным данным строят двухмерные графики; на графике для толщины наплавленного слоя находят кривую, соответствующую заданной глубине проплавления, и мысленно переносят ее на второй график, на котором изображены кривые припуска на механическую обработку. Определяют точку, при которой пересекается кривая заданной толщины наплавленного слоя с кривой, соответствующей минимальному припуску. Находят примерное значение фактора x1 , фиксируют это значение, а в цикле изме-
няют уже факторы x2 и x3 . И так, последовательно изменяя в
цикле попарно все факторы, повторяют несколько раз, пока не получат оптимальный вариант.
Реализация предлагаемого варианта оптимизации представлена на рис. 36–40. Была задана толщина наплавленного слоя 2,5 мм. При реализации на компьютере в пакете Mathcad 9 циклов были определены оптимальные значения факторов, при которых механический припуск не превышал 0,1 мм. Перевод кодированных значений факторов в реальные производится по формуле
xр xk x0 ,
где хр – реальное значение фактора; хk – кодированное значение фактора; – интервал варьирования факторов.
В результате получим следующие оптимальные параметры наплавки:
•скорость подачи электродной проволоки x1 62,4 м/ч;
•скорость наплавки x2 34,8 м/ч;
•шаг наплавки x3 5 мм/об;
•рабочее напряжение x4 22 В;
•относительное смещение электрода от зенита x5 0,395.
174
5.3. Оптимизация состава электродного покрытия
Качество покрытых электродов определяется комплексом их сварочно-технологических, санитарно-гигиенических характеристик, а также механических свойств сварных соединений. Исследуем влияние компонентов покрытия электродов на такие важные характеристики, как стабильность горения сварочной дуги и критическая температура хрупкости металла шва. Решающее влияние на комплекс указанных выше свойств оказывает состав покрытия. В настоящих исследованиях за основу был принят состав покрытия электродов МР-3. Ферромарганец был заменен на ферросилиций, остальные компоненты варьировались. Из состава покрытия были выделены пять основных наиболее значимых компонентов, которые стали факторами оптимизации – это содержание в составе покрытия силикомарганца (Х1), слюды (Х2), ферротитана (Х3), графита (Х4) и
мрамора (Х5).
Общее число опытов Nоб определяют по формуле
Nоб = Nя + n + n0 ,
где Nя – количество опытов «ядра» плана; n – количество опытов в звездных точках; n0 – количество опытов на нулевом уровне.
При проектировании состава электродного покрытия необходимо учитывать максимальное количество характеристик электродов. Это приводит к невозможности выбора единственного параметра оптимизации электродного покрытия и к необходимости обобщения нескольких функций отклика, полученных в результате реализации пятифакторного плана экспериментов.
Уровни и интервалы варьирования факторов представлены в табл. 43.
175
|
|
|
|
Т а б л и ц а 4 3 |
||
|
|
|
|
|
|
|
|
Х1 |
Х2 |
Х3 |
Х4 |
Х5 |
|
Уровни и интервалы |
(мра- |
|||||
(SiMn) |
(слюда) |
(FeTi) |
(графит) |
|||
|
мор) |
|||||
|
|
|
|
|
||
Основной уровень |
14 |
12 |
2 |
1 |
12 |
|
Интервал варьирования |
4 |
3 |
1 |
1 |
3 |
|
Нижний уровень |
10 |
9 |
1 |
0 |
9 |
|
Верхний уровень |
18 |
15 |
3 |
2 |
15 |
|
Звездные точки |
22 |
18 |
4 |
3 |
18 |
|
6 |
6 |
0 |
0 |
6 |
||
|
Обобщение нескольких откликов в один количественный признак связано с рядом трудностей. Каждый отклик имеет свой физический смысл и свою размерность. Чтобы объединить их, необходимо ввести для каждого из откликов некоторую однотипную безразмерную шкалу – это делает их сравнимыми. Выбор шкалы – непростая задача, зависящая от априорных сведений об откликах.
Матрица 2-го порядка для 5-факторного плана представлена в табл. 44.
Адекватные математические модели были построены по двум параметрам оптимизации: критической температуре хрупкости и коэффициенту вариации по току, характеризующему стабильность горения дуги:
Kvar I ,
Iср
где I среднеквадратичное отклонение; Iср среднее значение сварочного тока.
Обе характеристики определялись с помощью статистической обработки осциллограмм сварочного тока. Сварка производилась электродами диаметром 4 мм при среднем значении сварочного тока 135 A.
176
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Т а б л и ц а |
4 4 |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
№ |
|
|
|
|
|
|
Х12 |
Х22 |
Х32 |
Х42 |
Х52 |
2 |
3 |
4 |
5 |
3 |
4 |
5 |
4 |
5 |
5 |
п/п |
Х0 |
Х1 |
Х2 |
Х3 |
Х4 |
Х5 |
1 Х |
1 Х |
1 Х |
1 Х |
2 Х |
2 Х |
2 Х |
3 Х |
3 Х |
4 Х |
|||||
|
|
|
|
|
|
|
|
|
|
|
Х |
Х |
Х |
Х |
Х |
Х |
Х |
Х |
Х |
Х |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
2 |
1 |
–1 |
1 |
1 |
1 |
–1 |
1 |
1 |
1 |
1 |
1 |
–1 |
–1 |
–1 |
1 |
1 |
1 |
–1 |
1 |
–1 |
–1 |
3 |
1 |
1 |
–1 |
1 |
1 |
–1 |
1 |
1 |
1 |
1 |
1 |
–1 |
1 |
1 |
–1 |
–1 |
–1 |
1 |
1 |
–1 |
–1 |
4 |
1 |
–1 |
–1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
–1 |
–1 |
–1 |
–1 |
–1 |
–1 |
1 |
1 |
1 |
5 |
1 |
1 |
1 |
–1 |
1 |
–1 |
1 |
1 |
1 |
1 |
1 |
1 |
–1 |
1 |
–1 |
–1 |
1 |
–1 |
–1 |
1 |
–1 |
6 |
1 |
–1 |
1 |
–1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
–1 |
1 |
–1 |
–1 |
–1 |
1 |
1 |
–1 |
–1 |
1 |
7 |
1 |
1 |
–1 |
–1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
–1 |
–1 |
1 |
1 |
1 |
–1 |
–1 |
–1 |
–1 |
1 |
8 |
1 |
–1 |
–1 |
–1 |
1 |
–1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
–1 |
1 |
1 |
–1 |
1 |
–1 |
1 |
–1 |
9 |
1 |
1 |
1 |
1 |
–1 |
–1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
–1 |
–1 |
1 |
–1 |
–1 |
–1 |
–1 |
1 |
10 |
1 |
–1 |
1 |
1 |
–1 |
1 |
1 |
1 |
1 |
1 |
1 |
–1 |
–1 |
1 |
–1 |
1 |
–1 |
1 |
–1 |
1 |
–1 |
11 |
1 |
1 |
–1 |
1 |
–1 |
1 |
1 |
1 |
1 |
1 |
1 |
–1 |
1 |
–1 |
1 |
–1 |
1 |
–1 |
–1 |
1 |
–1 |
12 |
1 |
–1 |
–1 |
1 |
–1 |
–1 |
1 |
1 |
1 |
1 |
1 |
1 |
–1 |
1 |
1 |
–1 |
1 |
1 |
–1 |
–1 |
1 |
13 |
1 |
1 |
1 |
–1 |
–1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
–1 |
–1 |
1 |
–1 |
–1 |
1 |
1 |
–1 |
–1 |
14 |
1 |
–1 |
1 |
–1 |
–1 |
–1 |
1 |
1 |
1 |
1 |
1 |
–1 |
1 |
1 |
1 |
–1 |
–1 |
–1 |
1 |
1 |
1 |
15 |
1 |
1 |
–1 |
–1 |
–1 |
–1 |
1 |
1 |
1 |
1 |
1 |
–1 |
–1 |
–1 |
–1 |
1 |
1 |
1 |
1 |
1 |
1 |
16 |
1 |
–1 |
–1 |
–1 |
–1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
–1 |
1 |
1 |
–1 |
1 |
–1 |
–1 |
17 |
1 |
2 |
0 |
0 |
0 |
0 |
4 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
18 |
1 |
–2 |
0 |
0 |
0 |
0 |
4 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
177
17 8
О к о н ч а н и е т а б л . 4 4
№ |
0Х |
1Х |
2Х |
3Х |
4Х |
5Х |
12Х |
22Х |
32Х |
42Х |
52Х |
Х1Х2 |
Х1Х3 |
Х1Х4 |
Х1Х5 |
Х2Х3 |
Х2Х4 |
Х2Х5 |
Х3Х4 |
Х3Х5 |
Х4Х5 |
п/п |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
19 |
1 |
0 |
2 |
0 |
0 |
0 |
0 |
4 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
20 |
1 |
0 |
–2 |
0 |
0 |
0 |
0 |
4 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
21 |
1 |
0 |
0 |
2 |
0 |
0 |
0 |
0 |
4 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
22 |
1 |
0 |
0 |
–2 |
0 |
0 |
0 |
0 |
4 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
23 |
1 |
0 |
0 |
0 |
2 |
0 |
0 |
0 |
0 |
4 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
24 |
1 |
0 |
0 |
0 |
–2 |
0 |
0 |
0 |
0 |
4 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
25 |
1 |
0 |
0 |
0 |
0 |
2 |
0 |
0 |
0 |
0 |
4 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
26 |
1 |
0 |
0 |
0 |
0 |
–2 |
0 |
0 |
0 |
0 |
4 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
27 |
1 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
28 |
1 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
29 |
1 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
30 |
1 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
31 |
1 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
32 |
1 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
Одним из наиболее удобных способов построения обобщенного отклика является обобщенная функция желательности Харрингтона.
Для построения обобщенной функции желательности натуральные значения частных откликов были преобразованы в безразмерную шкалу желательности, или предпочтительности. Ее назначение – установление зависимости соответствия между откликами и чисто субъективной оценкой экспериментатора предпочтительности (желательности) того или иного отклика.
Чтобы получить шкалу желательности, можно использовать готовые таблицы соответствий между отношениями предпочтения в эмпирической и числовой системах (см. табл. 19).
Эти отметки на шкале желательности соответствуют некоторым точкам кривой, описываемой уравнением
d exp[exp( y)] ,
где у – значение частного отклика в масштабе кодированной шкалы (рис. 41).
На оси ординат нанесены значения желательности, изменяющиеся от 0 до 1. По оси абсцисс указаны значения частного отклика, записанные в условном масштабе. На шкалах Ткр и Kvar, расположенных под кодированной шкалой частных откликов у, указаны натуральные значения частных откликов.
Распределение натуральных значений частных откликов по шкале не обязательно должно быть равномерным.
После преобразования частных откликов в частные желательности находят обобщенный показатель – обобщенную функцию желательности (D), которая задается как среднее геометрическое частных желательностей:
n
D n d i .
i 1
179
Рис. 41. Функция желательности
Способ задания обобщенной функции желательности таков, что если, с одной стороны, хотя бы одна частная желательность равна нулю, то и обобщенная функция тоже будет равна нулю, с другой стороны, обобщенная функция равна единице тогда и только тогда, когда все частные желательности равны единице. Обобщенная функция желательности весьма чувствительна к малым значениям частных желательностей.
|
|
|
|
|
|
Т а б л и ц а 4 5 |
|
|
|
|
|
|
|
|
|
Номер |
Натуральные |
Частные |
Обобщ. |
|
|||
значения |
желательно- |
Оценка по шкале |
|||||
элек- |
отклик |
||||||
откликов |
сти |
|
желательности |
||||
трода |
Kvar |
Tкр, С |
d1 |
d2 |
D |
|
|
|
|
||||||
1788 |
0,126 |
–20 |
0,72 |
0,368 |
0,515 |
Удовлетворительно |
|
1789 |
0,150 |
–55 |
0,63 |
0,840 |
0,727 |
Хорошо |
|
1790 |
0,141 |
–10 |
0,649 |
0,192 |
0,353 |
Плохо |
|
180 |
|
|
|
|
|
|