热力学计算过程精.ppt

热力学计算过程第1页，本讲稿共11页授授 课内课内 容容1 CALPHAD技术2 热力学数据库TDB文件建立3 相图实验数据POP文件建立4 热力学参数优化:相互作用参数选择及par文件建立5 相图计算第2页，本讲稿共11页1 CALPHAD技术第3页，本讲稿共11页2 热力学数据库TDB文件的建立Element:(研究体系):ELEMENT CE FCC_A1 140.12 0 1.6600E+01!ELEMENT LA DOUBLE_HCP(ABAC)1.3891E+02 0 1.3600E+01!ELEMENT O 1/2_MOLE_O2(G)15.999 0 2.4502E+01!(元素符号 参考态结构,原子量,H298-H0K,S298K)Species:(没有参考态)SPECIES O2 O2!SPECIES CE+3 CE1/+3!SPECIES CEO2 CE1O2!SPECIES LA2O3 LA2O3!.第4页，本讲稿共11页2 热力学数据库TDB文件的建立Function:FUNCTION GHSERCE 298.-7160.519+84.23022*T-22.3664*T*LN(T)-.0067103*T*2-3.20773E-07*T*3-18117*T*(-1);1000.Y -79678.506+659.4604*T-101.32248*T*LN(T)+.026046487*T*2-1.930297E-06*T*3+11531707*T*(-1);2000.Y-14198.639+190.370192*T-37.6978*T*LN(T);4000.N!(SGTE unary database)第5页，本讲稿共11页2 热力学数据库TDB文件的建立Phase:(根据相图根据相图,每一相均包含有一系列的数据每一相均包含有一系列的数据)组员组员,结构结构,点阵点阵,表示相互作用的表达式，表示相互作用的表达式，Gibbs能，参数及其意义能，参数及其意义 PHASE LIQUID:L%1 1.0 !CONSTITUENT LIQUID:L:CE,LA:!PARAMETER G(LIQUID,CE;0)298.14+GLIQCE#;4000 N!PARAMETER G(LIQUID,LA;0)298.14 +GLALIQ;4000 N!PARAMETER L(LIQUID,CE,LA;0)298.14 V1+V2*T;6000 N!PARAMETER L(LIQUID,CE,LA;1)298.15 V3+V4*T;6000 N!(interactive parameter are added in par file)TDB 文件：第6页，本讲稿共11页3 TDB文件调用及Par文件建立调用TDB，并创建Par文件:S-S-F A-BA-PAR L(LIQUID,A,B;0)298V1+V2*T;6000NA-PAR.SAVESET-INTERACTIVE Go d Sw u A-B Def-ele A B GET GO PAR CREATE A-B根据体系特点增补相互作用参数第7页，本讲稿共11页4 实验数据收集及pop文件建立相图实验数据热力学实验数据 CREATE_NEW 1,1 CHANGE_STATUS PHASE LIQUID#2 BCC_A2 HCP_A3#2=FIX 0 SET_CONDITION P=P0 EXPERIMENT T=1713:20 EXPERIMENT Y(LIQUID#2,Y)=0.94:0.01EXPERIMENT X(BCC_A2,Y)=0.99:0.01 EXPERIMENT Y(HCP_A3#2,Y)=0.99:0.01第8页，本讲稿共11页5 热力学优化拟合实验数据（pop文件），优化相互作用参数 S-S-F A-BCOM A-BE-EREAD 1SASW-WE 0 1-LSE-WE 1C-EBS-O-V 1-10000E-EREAD 1C-E.第9页，本讲稿共11页6 热力学计算大部分实验数据可以在实验误差范围内被解释，参数合理S-S-F A-BCOM A-BE-EREAD 1SAGO P-3REIS-C P=101325 N=1 X(B)=0.1 T=800C-ES-A-V 1 X(B)0 1 0.01S-A-V 2 T 300 3000 10ADDMAPPOPL第10页，本讲稿共11页谢谢 谢谢 !第11页，本讲稿共11页