熱力學(xué)計算過程第1頁，共11頁，2023年，2月20日，星期六授課內(nèi)容1CALPHAD技術(shù)2熱力學(xué)數(shù)據(jù)庫TDB文件建立3相圖實驗數(shù)據(jù)POP文件建立4熱力學(xué)參數(shù)優(yōu)化:相互作用參數(shù)選擇及par文件建立5相圖計算第2頁，共11頁，2023年，2月20日，星期六1CALPHAD技術(shù)第3頁，共11頁，2023年，2月20日，星期六2熱力學(xué)數(shù)據(jù)庫TDB文件的建立Element

:(研究體系):ELEMENTCEFCC_A1140.1201.6600E+01!ELEMENTLADOUBLE_HCP(ABAC)1.3891E+0201.3600E+01!ELEMENTO1/2_MOLE_O2(G)15.99902.4502E+01!(元素符號參考態(tài)結(jié)構(gòu),原子量,H298-H0K,S298K)Species:(沒有參考態(tài))

SPECIESO2O2!SPECIESCE+3CE1/+3!SPECIESCEO2CE1O2!SPECIESLA2O3LA2O3!...........第4頁，共11頁，2023年，2月20日，星期六2熱力學(xué)數(shù)據(jù)庫TDB文件的建立Function:FUNCTIONGHSERCE298.-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!(SGTEunarydatabase)第5頁，共11頁，2023年，2月20日，星期六2熱力學(xué)數(shù)據(jù)庫TDB文件的建立Phase:(根據(jù)相圖,每一相均包含有一系列的數(shù)據(jù))組員,結(jié)構(gòu),點陣,表示相互作用的表達式，Gibbs能，參數(shù)及其意義PHASELIQUID:L%11.0!CONSTITUENTLIQUID:L:CE,LA:!PARAMETERG(LIQUID,CE;0)298.14+GLIQCE#;4000N!PARAMETERG(LIQUID,LA;0)298.14+GLALIQ;4000N!

PARAMETERL(LIQUID,CE,LA;0)298.14V1+V2*T;6000N!PARAMETERL(LIQUID,CE,LA;1)298.15V3+V4*T;6000N!(interactiveparameterareaddedinparfile)TDB文件：

第6頁，共11頁，2023年，2月20日，星期六3TDB文件調(diào)用及Par文件建立調(diào)用TDB，并創(chuàng)建Par文件:

S-S-FA-BA-PARL(LIQUID,A,B;0)298V1+V2*T;6000NA-PAR....SAVESET-INTERACTIVEGodSwuA-BDef-eleABGETGOPARCREATEA-B根據(jù)體系特點增補相互作用參數(shù)第7頁，共11頁，2023年，2月20日，星期六4實驗數(shù)據(jù)收集及pop文件建立相圖實驗數(shù)據(jù)熱力學(xué)實驗數(shù)據(jù)

CREATE_NEW1,1CHANGE_STATUSPHASELIQUID#2BCC_A2HCP_A3#2=FIX0SET_CONDITIONP=P0EXPERIMENTT=1713:20EXPERIMENTY(LIQUID#2,Y)=0.94:0.01EXPERIMENTX(BCC_A2,Y)=0.99:0.01EXPERIMENTY(HCP_A3#2,Y)=0.99:0.01第8頁，共11頁，2023年，2月20日，星期六5熱力學(xué)優(yōu)化擬合實驗數(shù)據(jù)（pop文件），優(yōu)化相互作用參數(shù)

S-S-FA-BCOMA-BE-EREAD1SASW-WE01-LSE-WE1C-EBS-O-V1-10000E-EREAD1C-E.....第9頁，共11頁，2023年，2月20日，星期六6熱力學(xué)計算大部分實驗數(shù)據(jù)可以在實驗誤差范圍內(nèi)被解釋，參數(shù)合理S-S-FA-BCOMA-BE-EREAD1SAGOP-3REIS-CP=101325N=1X(B)=0.1T=