开始 开始 读入优化参数 由终端输入初值 /读入优选实验数据 相平衡方 热力学函数表 程的构筑 达式的构筑 G、G2EG、△H 计算平衡成 实验数据变换， 实验数据变换， 计算不同温度时 分与温度 构筑目标函数 二分法计算四、” 的热力学函数值 线性回归拟合 计算值与实验值的比较 EG、 EG.△H Y 绘图 <达精度否下 输出优化参数 N N 将所有数据与图形输入数据库 <发散否？> 输人数据库、结束 结束 图1 QCVFIT程序框图 图2 QCOMP程序框图 Fig.1 The skeleton chart of QCVFIT program Fig.2 The skeleton chart of QCOMP program 3.2优化结果 将上述实验信息代入QCVFIT程序，优化得到ω的表达式(n=0)为： w=-49369.66-201884.3Y+867799.9y2-1320276r3+625541.6r4 (18) 中间化合物的生成自由能为： 3KCL()+YCl3(I)=K3YCle(s) △G=-30315.06+5.439T (J/mol) (19) KCI(I)+2YCls()=KY2Cl(s) △G=-36675.96十23.304T (J/mol) (20) 将优化结果代入QCOMP程序，计算该体系的相图和热力学性质，并与实验值比较（图3和图 4)。 3.3讨论 (1)图3是优化分析后给出的“推荐”相图，从图中可看出，除YC一侧的液相线与 实验值有差别外，其余部分的计算值与实验值吻合很好。按照文献〔3)的方法对相图的端部 进行检验后认为，计算的YC端的液相线走向比实验值的走向更趋合理。图3与文献〔3)给 出的优化相图计算基本上没有差别。 ·615·匕卜 实验数据变换 ， ‘ 实验数据变 换 ， 构筑 目标函数 二分法计算 。 、 叮 相 平衡方 热力学 函数表 程的构筑 达式的构筑 计算平衡成 计算不同温度时 分与温度 的热力学 函数值 图 程序框图 址甘 仁 图 地 程序框图 州 优化结果 将上述实验信息代入 程序 ， 优化得到 。 的表达式 月一 为 。 一 一 ’ 厂 一 十 刁 中间化合物的生成 自由能为 。 吕 △ 一 十 刚 二 △ 一 十 珊 将优化结果代入 材， 程序 ， 计算该体系的相图和热力学性质 ， 并与实验值 比较 图 和 图 。 讨论 图 是优化分析后 给 出的 “ 推荐 ” 相 图 ， 从 图 中可看出 ， 除 一侧的液相 线与 实验值有差别外 ， 其余部分的计算值与实验值吻合很好 。 按照文献 〔 〕 的方法对相 图的端部 进行检验后认为 ， 计算的 端的液相线走 向比实验值的走 向更趋合理 。 图 与文献 〔 〕 给 出的优化相 图计算基本上没有差别 。 · ·