4. Оформление отчета

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

2. Привести симплекс для изучаемой системы с экспериментальными точками.

3 Привести значения рассчитанных коэффициентов и уравнения математических моделей зависимостей механических свойств от температуры и времени старения.

4. Привести рассчитанные значения для контрольных сплавов и результат оценки адекватности модели.

5. Привести график зависимости свойства от времени при постоянной температуре, построенный с помощью ЭВМ.

5. Литература

1. Ф.С. Новик. Планирование эксперимента на симплексе при изучении металлических систем. М.:Металлургия, 1985.

6. Контрольные вопросы

1. В каких случаях можно использовать симплекс для изучения зависимостей различных показателей от технологических параметров?

2. Как осуществляют перевод натуральных значений независимых

переменных в кодовые?

3. Как выбирают контрольные точки?

4. По каким критериям оценивают адекватность моделей?

5. В чем состоит целесообразность использования треугольных симплексов для оптимизация режима термообработки?

6. Как изменится предельное значение  при оценке адекватности по t-критерию, если дисперсия опыта S уменьшится?

ЛАБОРАТОРНАЯ РАБОТА N6

ПОСТРОЕНИЕ ЛИНЕЙНЫХ МОДЕЛЕЙ СОСТАВ - СВОЙСТВО

МЕТОДОМ ПАССИВНОГО ЭКСПЕРИМЕНТА

(2 ч)

1. Цель работы

Освоить методику регрессионного анализа для построения линейных моделей зависимости свойств от нескольких факторов по имеющимся экспериментальным данным.

Ознакомится с использованием функции ЛИНЕЙН (LINEST) программы EXCEL.

2. ТЕОРЕТИЧЕСКОЕ ВВЕДЕНИЕ

Эксперимент считают ПАССИВНЫМ в том случае, когда исследователь не проводит запланированные опыты, а только наблюдает за неуправляемым процессом, т.е. фиксирует и обрабатывает получаемые результаты. Типичным примером может служить математическая обработка экспериментальных данных, приведенных в различных литературных источниках. Другим примером является анализ определенных свойств сплавов разных заводских плавок, на основании которого можно построить зависимости этих свойств от состава в

пределах марки сплава.

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

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

На первом этапе многофакторного эксперимента часто ограничиваются линейными регрессионными моделями:

Y =b₀ + b_i X_i (6.1)

Для построения таких моделей необходимо знать значения всех k факторов (X₁₁ , X₁₂ , ..X_1n ; X₂₁ , X₂₂ , ..X_2n ;... X_k1 , X_k2 , ..X_kn) параметра оптимизации (Y₁ ,Y₂ , Y_n ) для всех n опытов.

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

f(X₁,X₂ ,...X_k) принимала минимальное значение, т.е.:

[Y -f?X ,X ,...X %] = min

В случае одной независимой переменной X₁, когда строят уравнения

Y = b₀ + X₁*b₁ (6.2)

коэффициенты b и b рассчитывают по следующим формулам:

b = [(Y_i )( X_i )) - (X_i )( X_i Y­_i )] / [n (X_i )) -(X_i ) ] (6.3)

b = [n-( X_iY_i)( X_i)( Y_i) / [n((X_i )-( X_i)] (6.4)

где X_i и Y_i - известные значений независимой и зависимой переменных в i-том измерении, а n- число измерений.

Процедуру расчета коэффициентов можно существенно сократить если использовать программы для ЭВМ, например EXCEL. Для работы с этой программой предварительно необходимо составить матрицу исходных значений неизвестных (X) и известных (Y ) параметров. Например, в случае нахождения регресионной зависимости какого-либо свойства (Y) от трех факторов (X₁ ,X₂ и X₃ ) по результатам n измерений исходная матрица имеет следующий вид:

X₁₁ X₂₁ X₃₁ Y₁

X₁₂ X₂₂ X₃₂ Y₂

...................

X_1n X_2n X_3n Y_n

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

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

3. ПОРЯДОК ВЫПОЛНЕНИЯ РАБОТЫ

Поставлена задача построить математическую модель влияния состава на заданные свойства алюминиевого сплава методом пассивного эксперимента, используя известные экспериментальные значения (см варианты задания).

1. Составить матрицу исходных значений заданных факторов и параметров оптимизации.

2. Используя программу "EXCEL" (функция "ЛИНЕЙН"), рассчитать коэффициенты регрессии модели и их доверительные интервалы.

3. По полученному регрессионному уравнению рассчитать значения свойства (Y) при заданном составе. Оценить адекватность расчетных значений по t-критерию.

4. Построить график зависимости свойства от содержания одного из компонентов при постоянном уровне остальных (см. вариант задания).

4. ОФОРМЛЕНИЕ ОТЧЕТА

1. Привести таблицу, включающую исходные данные, рассчитанные значения коэффициентов регрессии, их доверительные интервалы.

2. Привести уравнение регрессии.

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

4. Привести график влияния одного из компонентов на заданное свойство.

5. ЛИТЕРАТУРА

1. Румшиский Л.З. Математическая обработка результатов эксперимента. Наука, 1971, С.108-122, 124-125.

6. КОНТРОЛЬНЫЕ ВОПРОСЫ

1. Чем отличается активный эксперимент от пассивного?

2. Каковы преимущества пассивного эксперимента перед активным?

3. Какие недостатки пассивного эксперимента по сравнению с активным?

4. В чем заключается определение коэффициентов регрессионной модели методом наименьших квадратов?

5. Какие исходные данные необходимо иметь для построения регрессионной модели?

6. Как можно оценить адекватность регрессионной модели?

7. Как можно использовать регрессионные уравнения?

8. Можно ли рассчитать коэффициенты линейной регрессионной модели для 8 факторов по результатам 8 измерений?

9. Сколько степеней свободы имеет линейная регрессионная модель для 10 факторов при 15 измерениях?

РАСЧЕТ КОЭФФИЦИЕНТОВ ЛИНЕЙНОЙ РЕГРЕССИИ
СПЛАВ АК5М2		Т6			Предел прочности, МПа
N	Si,%	Cu,%	Mg,%	Fe,%	экспер.	расчет	разница
1	4.5	3.4	0.3	1.1	262	259.9517	2.048335
2	4.9	3.2	0.5	0.8	296
3	5.2	2.4	0.3	0.7	281
4	5.8	2.6	0.4	0.9	302
5	5.6	2.2	0.6	1	332
6	5.3	1.8	0.7	0.5	362	358.9322	3.067838
7	4.1	2.6	0.3	0.6	275
8	5.1	2.8	0.4	0.7	286
9	5.6	3.1	0.5	0.8	300
10	4.8	2.7	0.6	1	324
		b4(Fe)	b3(Mg)	b2(Cu)	b1(Si)	b0
	значен.	2.251487	158.1986	-21.6007	3.113536	269.4469
	довер.	8.557026	11.44866	3.908479	2.904088	18.6936
коррел.	Y-довер	0.989899	4.032747	#Н/Д	#Н/Д	#Н/Д
F-критер.	степ.св.	122.4968	5	#Н/Д	#Н/Д	#Н/Д
		7968.685	81.31525	#Н/Д	#Н/Д	#Н/Д
		#Н/Д	#Н/Д	#Н/Д	#Н/Д	#Н/Д
		% Si	% Cu	% Mg	% Fe	Gb, Mpa
		6.5	3.4	0.8	1.3	357.4361

Рис.6.1. Шаблон для расчета множественной регрессионных линейных моделей.

ЛАБОРАТОРНАЯ РАБОТА N7

ИСПОЛЬЗОВАНИЕ МАТРИЧНЫХ ФУНКЦИЙ ДЛЯ РАСЧЕТА ОБЪЕМНЫХ ДОЛЕЙ В МНОГОКОМПОНЕНТНЫХ СПЛАВАХ

(2 ч)

1. ЦЕЛЬ РАБОТЫ

Ознакомиться с принципом расчета объемных (Q_V) и массовых (Q_M) долей в многокомпонентных сплавах с использованием матриц.

Освоить расчет Q_V и Q_M для заданного четверного сплава в программе EXCEL.

2. ТЕОРЕТИЧЕСКОЕ ВВЕДЕНИЕ

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

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

1- низкая точность оценки Q_M, которая в значительной мере определяется качеством рисунка;

2- практическая невозможность использовать этот метод для сплавов с большим числом компонентов.

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

этом случае используется свойство многокомпонентных сплавов, из которого следует постоянство состава всех фаз независимо от состава сплава. Рассчитать относительные массовые количества фаз можно решением системы N линейных уравнений баланса, в которых концентрация каждого компонента в сплаве (X_i ) определяется как сумма произведений концентраций этого компонента в каждой фазе (Y_ij ) на массовую долю (Q_Mi ) данной фазы:

Y₁₁ Q_M1 + Y₁₂ Q_M2 + ...+ Y_1N Q_MN = X₁

Y₂₁ Q_M1 + Y₂₂ Q_M2 + ...+ Y_2N Q_MN = X₂

..................................

Y_N1 Q_M1 + Y_N2 Q_M2 + ...+ Y_NN Q_MN = X_N (7.2)

Поскольку из N концентраций независимыми являются только N-1, а сумма массовых долей всех фаз равна 1 (или 100%) для определения искомых значений QM можно воспользоваться более удобной системой N линейных уравнений по сравнению с (7.1):

Y₁₁ Q_M1 + Y₁₂ Q_M2 + ...+ Y_1N Q_MN = X₁

Y₂₁ Q_M2 + Y₂₂ Q_M2 + ...+ Y_2N Q_MN = X₂

..................................

Y_N-1;1 Q_M1 + Y_N-1;2 Q_M2 + ...+ Y_N-1;N Q_MN = X_N (7.2)

Q_M1 + Q_M2 + .....+ Q_MN = 1

В общем случае решение системы линейных уравнений с N неизвестными находится через определители N-го порядка. Для рассматриваемой системы (7.2) матрица главного определителя имеет следующий вид:

Y₁₁ Y₁₂ ......... Y_1N

Y₂₁ Y₂₂ ..........Y_2N

= .................................. (7.3)

Y_N1 Y_N2 ........ Y_NN

Частные определители D_1, D₂ ....D_N получаются из D заменой элементов соответствующего столбца свободными членами (X₁, X₂ ,....X_N ,1). Вычислив определители легко рассчитать искомые значения Q_M :

Q_M = D /D (7.4)

Для нахождения объемных долей необходимо знать кроме Q_M

плотности всех фаз ( _i):

Q_Vi =(Q_Mi /₁)/ (Q_M1 / ₁ + Q_M2 / ₂ + .....Q_MN / _N ) (7.5)

Поскольку ручное вычисление определителей высоких порядков слишком трудоемкая работа, целесообразно использовать специальные программы для ЭВМ, которые позволяют проводить основные операции с матрицами. В частности, программа EXCEL делает возможным сделать шаблоны, с помощью которых легко рассчитывать для любого сплава, вводя исходные значения концентраций (X_i и Y_ij ) и плотностей ( _i): в определенные ячейки и получая результат в других ячейках. Более того, по сделанному для какой-либо системы шаблону можно легко создать новый шаблон для другой системы.

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

3. ПОРЯДОК ВЫПОЛНЕНИЯ РАБОТЫ

Известно, что структу расплава четверной системы Al-Si-Cu-Mg в закаленном (с 500 ⁰С) состоянии содержит алюминиевый твердый раствора (Al)-1, кремниевый твердый раствор (Si)-2 и промежуточные фазы CuAl₂-3 и W(CuMgSiAl)-4. Соответствующие значения их массовых долей даны в табл.7.1.

1. Используя исходные данные, составить все необходимые уравнения и матрицы определителей и с помощью программы EXCEL рассчитать значения последних. Далее рассчитать значения массовых и объемных долей всех фаз в заданном сплаве.

2. Используя готовый шаблон для программы EXCEL построить график влияния массовой доли одного из компонентов на величину Q_Vдной из фаз в пределах 4 фазной области (см. вариант задания).

Примеры размещения исходных и искомых значений в ячейках рабочего листа EXCEL показан на рис.7.1.

4. ОФОРМЛЕНИЕ ОТЧЕТА

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

2. Привести систему линейных уравнений и матрицы всех определителей, которые необходимы для расчета искомых значений.

3. Привести график влияния массовой доли заданного компонента на величину QV заданной фазы.

Таблица 7.1

Исходные данные для определения массовых и объемных долей фаз в системе Al-Si-Mg-Cu (фазовая область (Al)-(Si)-CuAl2-W)

Компоненты	Концентрации в фазах и сплаве* мас.%
	Ф1	Ф2	Ф3	Ф4	Сплав**
Si,	1	0,05	0,05	0,05
Cu,	4	0,05	52	19,2
Mg,	0,4	0,05	0,05	33
плотность, г/см3	2,7	2,6	4,3	2,8
QM, мас.%	?	?	?	?
QV, об.%	?	?	?	?

^*стальное алюминий; ^** см. варианты заданий

5. КОНТРОЛЬНЫЕ ВОПРОСЫ

1. Какие недостатки расчетно-графического определения относительных массовых долей?

2. Что может означать отрицательное расчетное значение массовой

доли какой-либо фазы?

3. Как можно проверить правильность расчета объемных долей фаз?

4. Напишите систему уравнений, которую необходимо решить для расчета значений Q_M и Q_V трехфазного сплава трехкомпонентной системы.

5. То же что и в п.4 для пятифазного сплава пятикомпонентой системы.

6. Можно ли рассчитать значения Q_M и Q_V для 3-фазного 4- компонентного сплава ? Что для этого необходимо знать?

РАСЧЕТ МАССОВЫХ (QM) И ОБЪЕМНЫХ (QV) ДОЛЕЙ В ЧЕТВЕРНЫХ СПЛАВАХ
	A-B-C-D	(А-основа)
	ИСХОДНЫЕ ДАННЫЕ (B6:F8)			и (B17:B20)
КОМПО-	ФАЗА-1	ФАЗА-2	ФАЗА-3	ФАЗА-4	СПЛАВ
НЕНТЫ	F1	F2	F3	F4
B,%	0.8	99.5	0	32.1	6.8
C,%	4	0	52.5	19.2	6.3
D,%	0.3	0	0	33	2
	1	1	1	1	1
DD=	-157375		РЕЗУЛЬТАТ (E11:E15) и (G11:G15)
D1=	-136661		QM1,%=	86.83731	QV1,%=	87.98432
D2=	-6980.27		QM2,%=	4.435431	QV2,%=	4.666864
D3=	-5439.02		QM3,%=	3.456079	QV3,%=	2.19876
D4=	-8295.53		QM4,%=	5.271176	QV4,%=	5.150058
			SUM=	100	SUM=	100
Плотность фаз
R1=	2.7	g/cm*3
R2=	2.6	g/cm*3
R3=	4.3	g/cm*3
R4=	2.8	g/cm*3

Рис.7.1. Шаблон для расчета массовых и объемных долей фаз в четверных сплавах

ЛАБОРАТОРНАЯ РАБОТА N8

РАСЧЕТ ШИХТЫ ВТОРИЧНЫХ МНОГОКОМПОНЕНТНЫХ СПЛАВОВ

(2ч )

1. ЦЕЛЬ РАБОТЫ

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

Освоить расчет вторичной шихты для заданного пятикомпонентного сплава в программе EXCEL.

2. ТЕОРЕТИЧЕСКОЕ ВВЕДЕНИЕ

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

В общем случае для приготовления N-компонентного сплава требуется N исходных шихтовых сплавов, поэтому для расчета требуемых количеств последних ?M % можно использовать систему N уравнений баланса:

Y₁₁ Q_M1 + Y₁₂ Q_M2 + ...+ Y_1N Q_MN = X₁

Y₂₁ Q_M1 + Y₂₂ Q_M2 + ...+ Y_2N Q_MN = X₂

.................................. (8.1)

Y_N1 Q_M1 + Y_N2 Q_M2 + ...+ Y_NN Q_MN = X_N

или

Y₁₁ Q_M1 + Y₁₂ Q_M2 + ...+ Y_1N Q_MN = X₁

Y₂₁ Q_M1 + Y₂₂ Q_M2 + ...+ Y_2N Q_MN = X₂

.................................. (8.2)

Y_N-1;1 Q_M1 + Y_N-1;2 Q_M2 + ...+ Y_N-1;N Q_MN = X_N-1

Q_M1 + Q_M2 ......+ Q_MN = 1

где Y_iJ - концентрация i-компонента в j-шихтовом сплаве, мас.%

X_i - концентрация i-компонента в требуемом сплаве, мас.%

Q_Mj - доля j-сплава в шихте (Q_Mj =M_j /M_j )

M_j - количество требуемого сплава, кг

Решая систему уравнений (8.1) или (8.2) через определители, как это рассмотрено в лабораторной работе N7 можно получить искомые значения M_j (M_j =Q_Mj M_j ). При использовании программы EXCEL целесообразно создавать шаблоны для расчета шихты любых многокомпонентных сплавов. Для расчета по имеющемуся шаблону необходимо ввести в заданные ячейки концентрации компонентов в шихтовых материалах в требуемом сплаве X, а также количество последнего. После этого в других заданных ячейках моментально появятся искомые значения M_j. Очевидно, что наличие хотя бы одного отрицательного значения M_j говорит о невозможности приготовления заданного сплава X из данного набора шихтовых материалов. Для тройных и четверных сплавов это означает нахождение точки состава сплава X вне соответствующего треугольника или тетраэдра.

Пример шаблона для расчета вторичной шихты для приготовления 10-компонентого сплава показан на рис.8.1.

3. ПОРЯДОК ВЫПОЛНЕНИЯ РАБОТЫ

Известны составы десяти готовых сплавов в виде отходов, которые можно использовать в качестве шихты для приготовления пятикомпонентного сплава заданного состава. Варианты задания приведены в табл.8.1. Требуется рассчитать шихту для приготовления заданного сплава, используя готовый шаблон программы EXCEL.

4. ОФОРМЛЕНИЕ ОТЧЕТА

1. Привести таблицу, включающую исходные данные и рассчитанные значения количеств шихтовых сплавов.

2. Привести систему линейных уравнений и матрицы всех определителей, которые необходимы для расчета искомых значений.

5. КОНТРОЛЬНЫЕ ВОПРОСЫ

1. Что может означать отрицательное расчетное значение количества

какого-либо шихтового сплава?

2. Как можно проверить правильность расчета шихты?

3. Напишите систему уравнений, которую необходимо решить для расчета шихты шестикомпонентного сплава из 6 шестикомпонентных готовых сплавов?

4. Можно ли приготовить 3 компонентный сплав из 3 шихтовых 4 компонентных сплавов?

5. Можно ли приготовить 4-компонентный сплав из 3-х шихтовых 4-компонентных сплавов?

6. Можно ли приготовить 5 компонентный сплав из 4-х шихтовых 5- компонентных сплавов?

7. В чем состоит целесообразность использования вторичной шихты?

РАСЧЕТ ШИХТЫ ПЯТИКОМПОНЕНТНЫХ СПЛАВОВ
A-B-C-D-E		(A-основа)		мас.%
КОМПО-	ШИХ-1	ШИХ-2	ШИХ-3	ШИХ-4	ШИХ-5	СПЛАВ	МАССА
НЕНТЫ							грамм
B,%	3	5	0	0	0	2	80
C,%	3	0	20	0	0	2
D,%	0	4	0	0	0	1.5
E,%	0	0	0	20	0	1
A,%	94	91	80	80	100	93.5
				РЕЗУЛЬТАТ
D0=	-480000			M1=	3.333333	грамм
D1=	-20000			M2=	30	грамм
D2=	-180000			M3=	7.5	грамм
D3=	-45000			M4=	4	грамм
D4=	-24000			M5=	35.16667	грамм
D5=	-211000			СУММА=	80	грамм

Рис.7.1. Шаблон для расчета шихты пятикомпонентных вторичных сплавов

ДОМАШНЕЕ ЗАДАНИЕ

ПОСТРОЕНИЕ ФАКТОРНЫХ ПЛАНОВ ПЕРВОГО ПОРЯДКА И ИХ ПРИМЕНЕНИЕ В МЕТОДЕ КРУТОГО ВОСХОЖДЕНИЯ

(6 ч)