《ExcelR电子表格编程在化工计算中的应用.doc》由会员分享，可在线阅读，更多相关《ExcelR电子表格编程在化工计算中的应用.doc（12页珍藏版）》请在金锄头文库上搜索。
1、ExcelR电子表格编程在化工计算中的应用摘要 本文从化工计算中的气体状态方程、简单物料衡算、化学平衡过程、循环过程四个典型方面，举例阐明了ExcelR电子表格编程在化工计算中的应用方法、技巧和步骤。关键词 化工计算 物料衡算 ExcelR电子表格 化工计算是我校为高级技工大专班化工工艺专业开设的一门专业基础课程，该课程主要涉及了化工过程的物料衡算、能量衡算以及物料与能量联算等化工计算问题。同时书中还介绍了，应用计算机编程解题的方法，并在每章附有若干源程序。但由于本书采用ortran语言编程的解题方法，这对高职学生来说，理解Fortran算法语言还有一定的学习难度。而且若要实现上机教学，又因编
2、写程序的复杂性，需花费学生大量的时间和精力，这显然对高级技工大专班的教学是不现实的。我在实际教学中，为了解决该课程中存在的这方面不足，选用了普及率广、应用性强、通俗易懂的Microsoft ExcelR电子表格计算机编程算法，从而实现了把化工计算教学和计算机应用结合起来。ExcelR算法相对简单易学，只要向学生提供简单指导和关键步骤，学生能够不费太大时间和精力，掌握它的算法步骤、技巧和思维，从而实现利用计算机快速、准确地解决化工计算问题。本文就ExcelR电子表格编程在化工计算中的应用，从气体状态方程、简单物料衡算、化学平衡过程、循环过程四个方面，进行举例阐明。这四个方面在化工计算中具有典型性
3、，能够涵盖本专业教学中的化工物料衡算的全部内容。目的只在因起本专业的师生们对ExcelR算法有一个全新的认识。一、在实际气体状态方程（R-K方程）中的应用1、R-K方程Redlich和Kwong于1949年在范德华方程的基础上提出的一个两参数方程： (1)将(1)变型为：F(v)=P-RT/(Vm-b)+a/T1/2Vm(Vm+b)=0 (2)它是形式最简单的常用实际气体状态方程，一般情况下有一定的准确度。式中和b亦为物质的特征参数，可由临界温度Tc和临界压力pc按下式计算而得：0.4278R2Tc2.5/Pc，b0.0867RTc/Pc2、用RK方程求取500K和18atm下正丁烷的摩尔体积
4、。已知正丁烷的临界温度425.2；临界压力37.5atm。3、利用ExcelR解RK方程的算法步骤：ABCDEF1R-K方程2正丁烷3姓名 时间4参数单位数据单位5Tc425.2KT500K6Pc37.5atmP18atm7R0.08206L.atm/g.mol.K89公式结果10Tr=E5/B5Tr=B1011a=0.4278*B72*B52.5/B6a=B1112b=0.0867*B7*B5/B6b=B1213F(v)=E6-(B7*E5)/(E13-B12)+B11/(E50.5*E13*(E13+B12)v2.037514F(v)=B1315理想气体体积=B7*E5/E6v=B15步骤
5、1 编写如上图所示的电子表格。表格的上方应编写出标题、名称、姓名和日期等标记性内容，为今后的复习和查阅提供方便；步骤2 在参数与数据格中输入各参数。并将临界参数、温度、压力和状态参数的单元格命名为Tc、Pc、T、P和R；步骤3 在10至15行的单元格内列出如图中方程；步骤4 利用工具栏里的单变量求解或规划求解命令求解，即改变单元格v(单元格E13)使F(v)(单元格E14)等于零。操作时，会出现如下一个对话框，在空格内填写如下信息：目标单元格$E$14目标值0可变单元格E13点击确定按钮。答案出现在电子表格中：E13为：2.0375，E14为2.8494E-10，此值虽不为零，但相对目标值已经
6、接近零了；步骤5 编写时，应确保电子表格中输入的正确性和所有单位的一致性。为了检查结果是否正确，可作为参照，输入理想气体状态方程，察看两者的结果是否接近。如相对接近，就可认为正确；步骤6 通过工具栏的公式审核模式，进行公式与计算结果切换查询，查看到的计算结果如下图：ABCDEF1R-K方程2正丁烷3姓名 时间4参数单位数据单位5Tc425.2000KT500.0000K6Pc37.5000atmP18.0000atm7R0.0821L.atm/g.mol.K89公式结果10Tr1.1759Tr1.175911a286.3886a286.388612b0.0807b0.080713F(v)2.8
7、494E-10v2.037514F(v)0.000015理想气体体积2.2794v2.2794二、在简单过程物料衡算中的应用1、物料衡算 在化工生产中物料衡算是非常有用，利用物料衡算来确定一个过程或系统是否有成本效益的，也能估计该过程经济上是否具有可持续发展能力。一个化工过程的物料流程，包括输入物流和输出物流。2、简单过程 是指仅有一个设备或一个单元操作，且没有发生化学反应的过程。如混合、蒸馏、蒸发、干燥、吸收、结晶、萃取等，这些过程的物料衡算为总物料衡算式和各组成的衡算式，构成多元方程，用代数法求解。3、混和过程物料衡算简图 如一种废酸中加入浓硫酸和浓硝酸配成混合酸，它的物料流程如下图：混
8、合废酸 Xkg HNO3 0.23H2SO4 0.57 H2O 0.20浓硝酸Zkg HNO3 0.90H2SO4 0H2O 0.10浓硫酸Ykg HNO3 0H2SO4 0.93H2O 0.07混合酸HNO3 0.27 H2SO4 0.60H2O 0.134、利用ExcelR对混和过程进行物料衡算步骤：ABCDEF1混和过程物料衡算2姓名 时间3选择计算基准：100kg混合酸4废酸组成浓硫酸组成浓硝酸组成混合酸组成5HNO30.2300.900.276H2SO40.570.9300.607H2O0.200.070.100.138合计11119混合酸100Kg10输出水量=B9*F7输入水量=
9、C7*D11+D7*D12+E7*D13Kg11废酸量xKg12浓硫酸量yKg=(B9*F6-C6*$D$11)/D613浓硝酸量zKg=(B9*F5-C5*$D$11)/E514总物料量100Kg=SUM(D11:D13)步骤1 编写如上图所示的电子表格。在电子表格的上方编写计算基准，原则上任何一股物流都可选作为计算基准，但是计算基准选择的恰当，可使计算简化，避免错误。这里选择以混合酸的量100kg为基准；步骤2 在物料组成格中输入各已知变量（如流量、组成）及单位，对一些未知的变量，可用符号x、y、z表示出；步骤3 在10至14列的单元格内列出如图中方程；步骤4 利用工具栏中的单变量求解或规
10、划求解命令求解。操作时，会出现如下一个对话框，在空格内填写如下信息：目标单元格$D$14目标值100可变单元格D10点击确定按钮。答案出现在电子表格中：D11为41.6893、D12为38.9645、D13为19.3460、D14为：100；步骤5 通过对比输出水量（单元格B10）与输入水量（单元格D10）是否相等，来检查结果正确与否。最终结果如下图：ABCDEF1混和过程物料衡算2姓名 时间3计算基准：100kg混合酸4废酸组成浓硫酸组成浓硝酸组成混合酸5HNO30.200.90.276H2SO40.60.9300.67H2O0.20.070.10.138合计11119混合酸100Kg10输
11、出水量13输入水量13Kg11废酸量xKg41.689312浓硫酸量yKg38.964513浓硝酸量zKg19.346014总物料量100Kg100三、在有化学反应平衡过程的物料衡算中的应用工业上的化学反应的过程，各反应物的实际用量，并不等于化学反应式中的理论量。为了使所需的反应顺利进行，常常需要使一些反应物用量过量，或者产生副产物、或者存在不参加反应的组分等等，这些因素会使化工计算比无化学反应过程变得复杂。对于化学反应过程的物料衡算，除了需要建立物料衡算式以外，常常还需要利用反应的平衡关系来计算产物的平衡组成。1、化学平衡表达式化学平衡是指在一定的条件下,可逆反应中,当正、逆反应速率相等,反
12、应混合物各组成成分的含量不再改变，反应达最大限度时，所处的一种动态平衡状态。我们可利用化学平衡关系确定混合物的组成。而化学平衡关系是由平衡常数确定的。对气体的温度不低于0，压强不高于10atm时，实际气体通常可近视为理想气体。当它的反应式为：aA + bB cC + dD其反应平衡常数则由下式表示：2、化学反应平衡过程的物料衡算简图 在接触法硫酸生产中，SO2被氧化成SO3。反应式为： SO2 +1/2O2 SO3氧化过程的温度为570，压力为1.1atm，平衡常数K为14.9。其他已知条件见物料流程简图：反应过程 mol SO2 0.08O2 0.09 N2 0.83molSO2 xSO3 yO2 zN2 0.833、利用ExcelR对化学反应平衡过程进行物料衡算。ABCDE1反应平衡过程物料衡算2姓名 时间3反应方程式SO2 +1/2O2 SO34数据单位5T5706P1.1atm7K14.98计算基准：100mol输入气体9SO2平衡转化率x0.716810物料成分输入量输入组成输出量输出组成11SO280.08=B11*(1
