理想溶液二元相图计算.docVIP

实验四 理想溶液二元相图计算 一． 实验目的 1．了解相图在材料科学和工程中的意义及理想溶液二元匀晶相图的计算方法。 2．学会使用C语言或其它语言编程计算理想溶液的二元匀晶相图。 3．了解当前国内外相图计算软件的现状。 二．实验原理 1．理想溶液相图计算理论 理想液态混合物中任意一组分B的化学势为：=+RTln 对于1mol理想液态混合物来说，结合化学势的定义，上式可以写成： =+RTln（1） 设指定的温度、压力下，B物质有α、β两相存在，达到平衡时，则有： =（2） 在我们研究的系统中，有A、B两组分存在，对于组分A，将式（1）代入（2）式得： 多项整理得： （3） 同理可得： （4） 利用（3）和（4）式即可计算理想溶液平衡两相的组成。 2．相图计算软件简介 目前集成热化学数据库和相图计算软件的系统主要包括瑞典皇家工学院材料科学与工程系为主开发的 Thermo-Calc 系统和加拿大蒙特利尔多学科性工业大学计算热力学中心为主开发的 FACT（Facility for the Analysis of Chemical Thermodynamics）系统。这些软件的共同特点是集成了具有自洽性的热化学数据库和先进的计算软件。 Pandat合金相图软件与热力学计算软件是美国CompuThermLLC公司开发的用于计算多元合金相图和热力学性能的软件包。可用于计算多种合金的标准平衡相图和热力学性能，用户也可使用自己的热力学数据进行相图与热力学计算。 三．问题描述 NiO-MgO为液、固相连续互溶二元体系。已知 NiO和MgO的熔点分别为1960℃和800℃，熔化热分别为52.3kJ/mol和77.4kJ/mol，以纯液态NiO作为NiO的标准态，纯固态MgO作为MgO的标准态，用C语言或其它高级语言编程计算该二元体系相图。 四．问题分析 1．计算模型 在上述问题中分别以纯液态NiO作为NiO的标准态，纯固态MgO作为MgO的标准态，则和近似计算式为： 将（3）和（4）式用于NiO-MgO系，设液相（l）为β，固相（s）为α，则有： （5） 同理可得： （6） 又因为，，则（6）式可写成： （7） 联立（5）和式（7）得： （8） （9） 由（8）和式（9）即可计算得到NiO-MgO的完全固溶体相图。 五．实验内容 程序流程图 源程序代码 #includeprocess.h #includestdio.h #includestdlib.h #includeconio.h #includedos.h #includemath.h #includegraphics.h #define TN (1960+273) #define TM (2800+273) #define T_CH 2 #define NUM (int)((TM-TN)/T_CH)+5 float R=8.3145,H_chN=52300.0,H_chM=77400; int n,X1,X2,Y1,Y2; float e_x(float H_ch,int t,int t_fus) { if(TN=tt=TM)return exp(H_ch*(t_fus-t)/(t_fus*R*t)); else printf(The temperature is out of range\n); return 0; } float x_L(int t) { float exN,exM,x_L; exM=e_x(H_chM,t,TM); exN=e_x(H_chN,t,TN); x_L=(1-exN)/(exM-exN); return x_L; } float x_S(int t) { float x_S; x_S=x_L(t)*e_x(H_chM,t,TM); return x_S; } coord_sys() { int i,dx,dy; char y_interval[5],x_interval[5]; X1=getmaxx()/2-150;X2=getmaxx()/2+150; Y1=getmaxy()/2-200;Y2=getmaxy()/2+150; dx=(X2-X1)/10;dy=(Y2-Y1