1、 1 C+课程设计报告 精馏塔理论塔板数的计算精馏塔理论塔板数的计算 化学化工学院化学化工学院 高级工程人才实验班（化学）高级工程人才实验班（化学） 赵静冰赵静冰 学号学号 15071101141507110114 2 目录 （1）设计报告* 3 （2）摘要* * 4 （3）程序设计思路*5 （4）流程框图* 6 （5）C 语言原始程序*15 （6）程序运行结果* 17 （7）后记*18 （8）参考文献*19 3 二、设计报告二、设计报告 1、撰写要求、撰写要求 （1）C 语言课程设计任务书； （2）目录； （3）摘要； （4）程序设计思路和流程框图； （5）C 语言原始程序； （6）程序运行
2、结果； （7）后记（程序调试过程出现问题的讨论） 。 2、印装和、印装和上交要求上交要求 （1）设计报告要求采用 A4 纸双面打印，订书针装订。 （2）将源程序按统一文件名 100a-bc 形式发给辅导老师 4 课程设计报告课程设计报告 摘要摘要：此次 C+程序课程设计，是在给定的设计条件下，利用梯 级图解法计算苯甲苯精馏塔理论塔板数。 第一个程序使用了直接计 算法解方程组，计算出塔顶 D 和塔底 W，使用最小二乘拟合法和高斯 消元法求得相平衡方程， 同时使用牛顿迭代的法求得 q 线方程与相平 衡方程的交点以及 q 线方程与精馏段操作线方程等的交点， 循环使用 牛顿迭代法最后求得精馏塔理论塔板
3、的数目和进料板的位置； 在设计 的过程中所需要的 C+语言方法有：循环、选择、数组、绘图等的综 合使用，最后完成整个设计过程。 5 程序设计思路程序设计思路：由于没有学习过化工原理方面的知识，所以初开 始，对于陌生的专有名词感觉这次的设计很不好做，不知从何下手。 后来老师给出了公式等推导原理， 并在老师讲解后我才大概懂得题目 的要求，参考了之前学过的课本和一些范例，再询问其他合作同学， 我终于有了设计思路。 我将整个程序设计过程分为两步， 计算部分和 作图部分，而计算又分为两个部分， 程序一过程使用了直接计算法解 方程组，计算出塔顶 D 和塔底 W，程序二过程先利用高斯消元法求得 相平衡方程，
4、再利用牛顿迭代法求得精馏塔理论塔板数的过程。 1 1、 高斯函数高斯函数：利用高斯主元素消去法可以求解平衡线方程中的系数 及常数 a0,a1,a2,a3。 2 2、 矩阵板块矩阵板块：我们在求解平衡线方程中的系数及常数 a0,a1,a2,a3 时，有一个问题，那就是矩阵的行数和列数不相同，不能完全 照搬书上所说的高斯主元素，而是需要利用数据的拟合。我们 可以用矩阵的转置与矩阵的乘法 （参考实际程序） ， 达到此目的。 3 3、 牛顿迭代法牛顿迭代法：通过翻阅资料，运用牛顿迭代法是求方程根，在所， 先定义一下牛顿迭代法，然后再在所需要的地方调用一下， 4 4、 循环板块循环板块：数组可以解决数据
5、庞大的问题，循环的过程中，输出 其角标即可。采用的是 do.while 语句。 6 C+设计流程图： 7 C C+原始程序：原始程序： 程序一 : /求解方程组 #include #define F 46.61 /F 为处理量 #define XF 0.45 /XF 为苯-甲苯混合液中苯的含量 #define XD 0.996 /XD 为塔顶产品摩尔分率 #define XW 0.0118 /XW 为塔底产品摩尔分率 int main() double D=0,W=0; /D 为塔顶，W 为塔底 D=(F*XF-F*XW)/(XD-XW); W=(F*XD-F*XF)/(XD-XW); cout=1e-6); return xn1; void main() loop: system(“cls“); couti;j-) s+=Hij*aj; ai=(Hi4-s)/Hii; coutxf; coutxd; coutxw; coutt; coutendl; R=t*Rmin; xq=xf; yq=R*xq/(R+1)+xd/(R+1); cout“xe=“xesetw(5)“ye=“yeendl; coutendl; cout“最小回流比为：“setw(5)
