Определение термодинамических активностей компонентов бронзы БрБ2
Рефераты >> Химия >> Определение термодинамических активностей компонентов бронзы БрБ2

Табл. 2.3. Аппроксимация зависимостей Q=f(T).

Линия

Полином

R2

1

Q = 33,285T - 17925

0,7556

Q = -0,1902T2 + 237T - 72123

0,7885

Q = 0,0128T3 - 20,674T2 + 11166T - 2E+06

0,9489

Q = -0,0001T4 + 0,2512T3 - 211,76T2 + 79045T - 1E+07

0,9586

Q = -1E-05T5 + 0,0294T4 - 31,331T3 + 16641T2 - 4E+06T + 5E+08

1

2

Q = -31,278T + 29506

0,9218

Q = 0,2428T2 - 291,32T + 98689

0,9959

Видно, что для линии 1 высоких значений R2 удаётся достичь только при больших степенях полинома. К сожалению, при этом не очень точно вычисляются их коэффициенты. К тому же, с такими зависимостями трудно работать. Всё это послужило причиной того, что от данного способа автор работы отказался.

Способ №2. Было принято решение разделить функции на три части соответствующие температурам для первой части, для второй и для третьей (на рис. 2.1 эти части разделены вертикальными прямыми). На каждом из этих отрезков зависимость можно аппроксимировать полиномом меньшей степени. Результаты приведены в таблице 2.4.

Табл. 2.4. Аппроксимация частей зависимости Q=Q(T).

Линия

Часть

Полином

R2

1

1

Q = 76,812T - 39259

0,9437

Q = -1,2995T2 + 1371,1T - 361006

1

2

Q = -46,012T + 24707

1

3

Q = 51,263T - 28567

0,9981

Q = -0,1545T2 + 228,27T - 79216

1

2

1

Q = -51,085T + 39360

0,9991

Q = -0,1052T2 + 53,71T + 13310

1

2

Q = -27,883T + 27204

1

3

Q = -13,086T + 19091

0,9994

Q = 0,0224T2 - 38,784T + 26444

1

Задав таким образом зависимости Q=f(T) как полиномы второй степени и зафиксировав один из параметров x, N, T, нужно решить систему (2.8). В этом случае система будет состоять из двух трансцендентных уравнений, и решить их совместно можно только численными методами. Автору работы не удалось этого сделать.

Поэтому было принято решение пожертвовать точностью аппроксимации функций Q=f(T) и определить их как линейные зависимости. В этом случае Q=aT+b и температура будет входить в уравнения системы (2.8) только в первой степени, что позволяет исключить её, как неизвестное.

Воспользуемся условными обозначениями, которые уже были использованы ранее.

Пусть , а . Тогда первое уравнение системы (2.9) запишется в виде:

(2.15)

Если перенести все слагаемые, содержащие Т, в левую часть, а все остальные – в правую часть уравнения, то получится:

(2.16)

Осталось только выразить температуру в явном виде:

(2.17)

Аналогично нужно выразить температуру и из второго уравнения системы (2.9):

(2.18)

(2.19)

(2.20)

Приравняв правые части равенств (2.17) и (2.20) и умножив их на -1, приведём уравнение к окончательному виду:

(2.21)

Параметра а и b определим из данных таблицы 2.4. Чтобы решить трансцендентное уравнение (2.21), нужно задаться одним из параметров x, или n и численными методами подобрать второй параметр, а затем определить и температуру по любому из уравнений (2.17) или (2.20).

Для решения была использована надстройка «поиск решения» пакета Microsoft Excel. Результаты решения представлены в таблице 2.5.

Табл. 2.5. Рассчитанный купол расслаивания твёрдого раствора при разных температурах

t, oC

Состав α-фазы (Cu)

Состав γ-фазы (Ni)

x1

x2

N1

N2

0

0,727

0,273

2,8E-06

0,999997

25

0,723

0,277

0,000014

0,999986

40

0,72

0,28

0,000035

0,999965

83

0,71

0,29

0,00027

0,99973

116

0,70

0,30

0,001

0,999

141

0,69

0,31

0,002

0,998

161

0,68

0,32

0,004

0,996

178

0,67

0,33

0,007

0,993

191

0,66

0,34

0,010

0,990

203

0,65

0,35

0,014

0,986

241

0,60

0,40

0,042

0,958

261

0,55

0,45

0,061

0,939

279

0,50

0,50

0,077

0,923

307

0,45

0,55

0,128

0,872

322

0,40

0,60

0,174

0,826

331

0,35

0,65

0,224

0,776

334

0,30

0,70

0,273

0,727

334

0,285

0,715

0,285

0,715


Страница: