4结论 改进的SUMT方法，对用最小自由能原理求解化学平衡是行之有效的，计算精度已优于 国外同类程序；在微机上运算结果表明，30个变量以下的化学平衡问题完全可以在微机上运 行，10个变量在Pc/XT机上占机时5min,20个变量约100min,30个变量5h。 附：全文符号： G:体系总自由能；m:相数；：组元数，C:元素数；x:无相中j组元摩尔数：g:k相 中j组元的标准摩尔自由能；R:气体常数；T:绝对温度(K)；五(x):元素的等式约束向 量；E:k相中i元素在j组元中的系数；b,:i元素的总量；9：()：不等式约束的表达式向 量；T:惩罚因子；P(x,r):惩罚函数；f(2):目标函数；e1,2,e1:分别为收敛准则规定 的精度，一般为104~105。 致谢：忠谢何晓峰博士对本工作所提供的帮助！ 参考文献 .1 Smith W R.Theoretical Chemistr Advance and Perspectives.New York:Academic Press, 1980,5 2 Villars D S.J of Physical Chemistry,1959,63:521 3许志宏，王乐珊.无机热化学数据库.北京：科学出版社，1987 4 Gautam C W.AICHE J,1979,25:991 5 White C W.AICHE J,1981,27:466 6 Himmelblau D M.Applied Nonlinear Programming.New York:McGraw Hill Book Compa- ny,1972 7万耀青，梁庚荣.最优化计算方法常用程序汇编.北京：工人出版社，1983 8 Burin I,Knacke O.Thermochemical Properties of Inorganic Substances.USA:Springer Press 73),Supplement ('77 )1977 ·508·￾ 结论 改进的 ￾￾￾￾ 方法 ， 对用最小 自由能原理求解化学平衡是行之有效的 ， 计算精度已优于 国外 同类程序 ￾ 在微机上运算结果表 明 ， ￾ 个变量 以下的化学平衡问题完全可 以在微机上运 行 ， ￾￾个变量在 ￾￾￾￾ 机上 占机时 ￾￾￾￾ ， ￾￾个变量约 ￾￾￾￾￾￾ ， ￾￾个变量 ￾￾ 。 附 ￾ 全文符号 ￾ ￾ ￾ 体系总 自由能 ￾ ￾ ￾ 相数 ￾ ， ￾ 组元数 ￾ ￾ ￾ 元素数 ￾ 孙 ￾ ￾相 中 了组元摩尔数 ￾ 嵘 ￾ ￾ 相 中 ￾组元的标准摩尔 自由能 ￾ ￾￾ 气体常数 ￾ ￾ ￾ 绝对温度 ￾￾￾ ￾， ￾勃 ￾ ‘元素的等式约束 向 量 ￾ ￾， ￾ ￾相中 ￾元素在 ，组元 中的系数 ￾ ￾‘￾ ‘元素的总量 ￾ ￾ ， ￾幻 ￾ 不等式约束的表达式 向 量 ￾ ， ￾ 惩罚 因子 ￾ 尸 ‘ ， ， ￾ ￾ 惩罚 函数 ￾ ￾ ￾幻 ￾ 目标函数 ￾ 。￾ ， 。￾ ， 。￾ ￾ 分别为收敛准则规定 的精度 ， 一般为 一。 一 ‘一 ￾￾ 一 ￾。 致谢 ￾ 感谢何晓峰博士 对本工作所提供的帮助 ￾ 参 考 文 献 ￾￾￾￾￾ ￾ ￾ ￾ ￾￾￾￾￾￾￾￾￾】￾ ￾￾￾￾ ￾￾￾ ￾￾ ￾￾￾ ￾￾￾￾￾￾￾￾￾￾ ￾ ￾￾ ￾￾￾￾ ￾￾￾￾￾￾￾￾ ￾￾￾￾ ， ￾￾￾￾ ， ￾ ￾辽￾￾ ￾ ￾ ￾ ￾￾￾￾￾￾￾￾￾￾ ￾￾￾￾￾￾￾￾￾ ， ￾￾￾￾ ， ￾￾￾ ￾￾￾ 许志宏 ， 王乐珊 ￾ 无机热化学数据库 ￾ 北京 ￾ 科学 出版社 ， ￾￾￾￾ ￾￾￾饭￾ ￾ ￾ ￾ ￾￾ ￾￾ ￾ ， ￾￾￾￾ ， ￾￾￾ ￾￾￾ ￾￾￾￾ ￾ 锐 ￾ ￾￾￾￾￾ ￾， ￾￾￾￾ ， ￾￾￾ ￾￾￾ ￾如￾￾￾￾恤￾ ￾ ￾ ￾ ￾￾￾￾￾ ￾￾￾￾￾￾￾￾ ￾ ￾￾￾￾￾￾ ￾ ￾￾￾ ￾￾￾￾￾ ￾￾￾￾￾、￾ 凡￾ 致阅￾ ￾￾￾￾￾￾ ￾￾， ￾￾￾￾ 万耀青 ， 梁庚荣 ￾ 最优化计算方法常用程序汇编 ￾ 北京 ￾ 工人 出版社 ， ￾￾￾￾ ￾ ￾￾ ￾ ， ￾￾￾￾￾￾ ￾ ￾ ￾￾￾￾￾￾￾而份￾￾￾￾￾川 ￾￾ ￾ ￾￾￾￾￾￾￾￾ ￾￾比￾￾￾￾￾ ￾ ￾￾￾ ￾ ￾￾￾￾￾￾￾￾ ￾￾ ￾ ‘ ￾￾￾ ， ￾￾￾￾￾￾￾￾￾￾ ￾ ’ ￾￾ ￾ ， ￾￾￾￾ ￾￾￾