Математическое моделирование и основы научных исследований в сварке
..pdfНа рис. 45 приведен листинг моделирования 3-го и 4-го циклов и графики. В 3-м цикле (Т3) изменялось одновременно содержание графита (по оси ординат) и мрамора (по оси абсцисс). Содержание ферротитана в закодированном виде было принято равным 0,8, содержание остальных компонентов – на нулевом уровне.
Рис. 45. Листинг моделирования при поиске оптимального состава покрытия, при котором обеспечивается наименьшая критическая температура хрупкости металла шва. 3-й и 4-й циклы. Рассчитанные значения факторов
131
Стр. 131 |
ЭБ ПНИПУ (elib.pstu.ru) |
Из приведенного рисунка видно, что графит не влияет на Ткр, а увеличение содержания мрамора снижает критическую температуру хрупкости металла шва. Содержание мрамора в закодированном виде было принято равным 1,8. В 4-м цикле (Т4) снова изменяли содержание силикомарганца и слюды при найденных значениях остальных факторов: Х3 = 0,8; Х4 = 0; Х5 = 1,8. Видно, что при увеличении содержания мрамора для снижения критической температуры хрупкости количество слюды необходимо уменьшать, а содержание силико-марганца выдерживать на нулевом уровне.
Таким образом, чтобы получить металл шва с наименьшей критической температурой хрупкости, необходимо выдержать следующее содержание компонетов в составе покрытия: силикомарганец – 14 %; слюда – 6–7 %; ферротитан – 2,8 %; графит – 1 %; мрамор – 17,4 %.
Выше было показано, что для получения оптимального состава покрытия, обеспечивающего получение одновременно наименьшего коэффициента вариации (соответствующего высокой стабильности горения дуги) и наименьшей критической температуры хрупкости, можно использовать обобщенную функцию желательности.
5.3.3. Оптимизация состава покрытия по обобщенной функции желательности
Частные и обобщенная желательности приведены в табл. 37. Рассчитаны они по соответствующим экспериментальным данным и графику, приведенному на рис. 37. Оси желательности составляются исследователем на основе опыта. Для составления адекватной математической модели необходимо повторить расчеты в том же порядке, но за параметр оптимизации принять обобщенную желательность. На рис. 46 приведен листинг расчета коэффициентов нового уравнения регрессии и ихзначимости.
При сопоставлении двух приведенных матриц для доверительных интервалов и коэффициентов уравнения регрессии видно, что незначимыми являются коэффициенты при следующих членах уравнения:
X 4 , X 22 , X 42 , X 52 , X1 X 2 , X1 X 5 , X 2 X 3 , X 3 X 5 , X 4 X 5 .
На рис. 47 приведен листинг расчета адекватности модели и запись в отдельный файл значимых коэффициентов. Расчет показывает, что вычисленный критерий Фишера составляет 4,441 при табличном значениидля данных условий 4,619. Поскольку расчетный критерий Фишера
132
Стр. 132 |
ЭБ ПНИПУ (elib.pstu.ru) |
Рис. 46. Листинг расчета коэффициентов уравнения регрессии и их значимости для обобщенной желательности
меньше табличного, можно принять полученную модель адекватной. Тогда уравнение регрессии принимает следующий вид:
D = 0, 7 + 0, 036 X1 + 0, 038X 2 + 0, 055X 3 + 0, 049 X 5 − 0, 069 X12 − −0, 049 X 32 − 0, 025X1 X 3 − 0, 06 X1 X 4 − 0,1X 2 X 5 − 0, 0397 X 3 X 4 .
Поиск оптимального состава покрытия по наибольшей обобщенной желательности можно, как в предыдущих случаях, осуществить поочередным варьированием всех факторов от –2 до +2 (в закодированном
133
Стр. 133 |
ЭБ ПНИПУ (elib.pstu.ru) |
Рис. 47. Листинг расчета адекватности полученной модели
виде). На рис. 48 приведен листинг расчета для 1-го и 2-го циклов, при которых сначала изменялись в цикле силикомарганец и слюда, а затем слюда и ферротитан. По 1-му графику (D1) сначала выбираем количество силикомарганца – 0,15 (в закодированном виде).
Далее изменяем во 2-м цикле содержание слюды и ферротитана. По графику (D2) видно, что для более высокой обобщенной желательности содержание слюды (по оси ординат) необходимо выдерживать в пределах 1,8–2,0. Таким образом, после 1-го и 2-го циклов примерное содержание силикомарганца и слюды составляет соответственно 0,15 и 1,8.
На рис. 49 представлены листинг расчета и полученные результаты (в виде графиков изолиний) при 3-м и 4-м циклах. При 3-м цикле изменялось содержание ферротитана и графита, при 4-м – графита и мрамора.
134
Стр. 134 |
ЭБ ПНИПУ (elib.pstu.ru) |
Рис. 48. Листинг расчета при поиске оптимального состава покрытия по наибольшему значению обобщенной функции желательности.
1-й и 2-й циклы
По оси абсцисс на графике D3 можно выбрать оптимальное значение ферротитана (его значение составляет 1). По графику D4 видно, что содержание графита и мрамора должно находиться на минимуме. Примем содержание этих компонентов равным –1,8.
Для более точного определения оптимального содержания компонентов покрытия желательно повторить 1-й и 2-й циклы для силикомарганца, слюды и ферротитана. Связано это с тем, что при 1-м цикле (D1) ферротитан, графит и мрамор принимались на нулевом уровне. После расчета по четырем циклам было определено примерное содержание этих компонентов, поэтому повторный расчет позволит уточнить их содержание (рис. 50).
135
Стр. 135 |
ЭБ ПНИПУ (elib.pstu.ru) |
Рис. 49. Листинг расчета при поиске оптимального состава покрытия по наибольшему значению обобщенной функции желательности.
3-й и 4-й циклы
Таким образом, наибольшую обобщенную желательность можно получить при следующем содержании компонентов покрытия: силикомарганец – 16,0 %; слюда – 16,5 %; ферротитан – 2,7 %; графит – 0 %; мрамор – 7,5 %. По имеющимся моделям, характеризующим стабильность горения сварочной дуги (коэффициент Kv) и критическую температуру хрупкости (Ткр), можно рассчитать значения Kv и Ткр, соответствующие указанному выше содержанию компонентов. На рис. 51 приведен листинг расчета.
136
Стр. 136 |
ЭБ ПНИПУ (elib.pstu.ru) |
Рис. 50. Листинг повторного расчета по 1-му и 2-му циклу и расчет оптимального состава покрытия по наибольшему значению обобщенной функции желательности
137
Стр. 137 |
ЭБ ПНИПУ (elib.pstu.ru) |
Рис. 51. Листинг расчета величин Kv и Ткр при наибольшем значении обобщенной функции желательности
|
|
|
Таблица 38 |
|
|
|
|
|
Оптимальныйсоставдля получения |
Оптимальный состав |
|
|
наивысшей стабильности |
понаибольшему значению |
|
Компоненты |
инаименьшейкритической |
обобщенной функции |
|
покрытия, % |
температурыхрупкости |
желательности |
|
|
Kv = 0,062 |
Tкр = –94 °С |
D = 1; Kv = 0,119; |
|
Tкр = –70 °С |
||
|
|
|
|
Силикомарганец |
20–22 |
14,0 |
16,0 |
Слюда |
6–18 |
6,0–7,0 |
16,5 |
Ферротитан |
0 |
2,8 |
2,7 |
Графит |
0 |
1,0 |
0,0 |
Мрамор |
16–18 |
17,4 |
7,5 |
В табл. 38 обобщены полученные результаты моделирования. В реальности получить обобщенную желательность, равную единице, практически невозможно, как и Ткр = –94 °С, поэтому после моделирования необходимо провести дополнительные эксперименты при всех вариантах, указанных в табл. 38.
138
Стр. 138 |
ЭБ ПНИПУ (elib.pstu.ru) |
СПИСОК ЛИТЕРАТУРЫ
1.Романенко В.Н., Орлов А.Г., Никитина Г.В. Книга для начинающего исследователя-химика. – Л.: Химия, 1987. − 280 с.
2.Статистические методы в инженерных исследованиях (лабораторный практикум): учеб. пособие / В.П. Бородюк, А.П. Вощинин,
А.З. Иванов [и др.]. – М.: Высш. шк., 1983. – 216 с.
3.Павловский З. Введение в математическую статистику / под ред. Ф.Д. Лившица; пер. спольск. В.Д. Маникера. – М.: Статистика, 1967. – 286 с.
4.Поллард Дж. Справочник по вычислительным методам статистики / под ред. и с предисл. Е.М. Четыркина; пер. с англ. В.С. Занадворова. – М.: Финансы и статистика, 1982. – 344 с.
5.Спиридонов А.А., Васильев Н.Г. Планирование эксперимента при исследовании и оптимизации технологических процессов: учеб. пособие / Урал. политехн. ин-т. – Свердловск, 1975. – 140 с.
6.Спиридонов А.А. Планирование эксперимента при исследовании технологических процессов. – М.: Машиностроение, 1981. – 184 с.
7.Адлер Ю.П., Маркова Е.В., Грановский Ю.В. Планирование эксперимента при поиске оптимальных условий. – М.: Наука, 1976. – 254 с.
8.Горский В.Г., Адлер Ю.П. Планирование промышленных экспериментов. – М.: Металлургия, 1974. – 264 с.
9.Львовский Е.Н. Статистические методы построения эмпирических формул: учеб. пособие для втузов. – М.: Высш. шк., 1988. – 239 с.
10.Налимов В.В., Чернова Н.А. Статистические методы планирования экстремальных экспериментов. – М.: Наука, 1965. – 338 с.
11.Бондарь А.Г., Статюха Г.А. Планирование эксперимента в химической технологии. – Киев: Вища школа, 1976. – 183 с.
12.Тамразов А.М. Планирование и анализ регрессионных экспериментов в технологических исследованиях. – Киев: Наукова думка,
1987. – 175 с.
13.Рузинов Л.П. Статистические методы оптимизации химических процессов. – М.: Химия, 1980. – 199 с.
14.Сидняев Н.И. Теория планирования эксперимента и анализ статистических данных: учебное пособие для магистров. – М.: Юрайт,
2012. – 399 с.
15. Потапов Б.Ф., Бульбович Р.В., Крюков А.Ю. Начала инженерного творчества: учебное пособие. – Пермь: Изд-во Перм. гос. техн.
ун-та, 2010. – 189 с.
139
Стр. 139 |
ЭБ ПНИПУ (elib.pstu.ru) |
Учебное издание
Летягин Игорь Юрьевич
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
ИОСНОВЫ НАУЧНЫХ ИССЛЕДОВАНИЙ
ВСВАРКЕ
Ч. 1. Статистическая обработка и планирование эксперимента
Учебное пособие
Редактор и корректор Е.Б. Денисова
__________________________________________________________
Подписано в печать 18.04.2014. Формат 70×100/16. Усл. печ. л. 11,61. Тираж 15 экз. Заказ № 65/2014.
__________________________________________________________
Издательство Пермского национального исследовательского
политехнического университета.
Адрес: 614990, г. Пермь, Комсомольский проспект, 29, к. 113.
Тел. (342) 219-80-33.
Стр. 140 |
ЭБ ПНИПУ (elib.pstu.ru) |