2 方案论证…………………………………………….………………… (4)
2.1 精馏塔类型……………………………………………………………（4）
2.2 精馏压力………………………………………………………………（4）
2.3 进料方式………………………………………………………………（4）
2.4 填料类型………………………………………………………………（5）
2.5 加热方式………………………………………………………………（5）
2.6 回流比的选择…………………………………………………………（6）
3 数学模型的建立…… …………………………………………...……（6）
3.1 精馏塔塔体年投资折旧费及维修费用………………………………（6）
3.2 冷凝器年运转费用……………………………………………………（8）
3.3 再沸器年运转费用…………………………………………………（8）
3.4 填料年折旧费…………………………………………………………（8）
4 数学模型的求解………………………………………………...……（8）
4.1 数学模型决策变量分析………………………………………………（8）
4.2 主要工艺参数的求解…………………………………………………（9）
4.2.1 塔径的计算……………………………………………………….…（9）
4.2.2 理论塔板数的计算………………………………………….……（10）
4.2.3 冷却水用量及冷凝器传热面积的计算………………..….....(11)
4.2.4 冷凝器冷却水最佳出口温度的确定……………………….…（11）
4.2.5 直接蒸汽加热蒸汽流量的计算……………………………...…（12）
4.3 数学模型的求解………………………………………………………（12）
4.3.1 单变量最优化方法………………………………………………….（13）
4.3.2 优化设计程序框图………………………………………………….（13）
4.3.2.1 求解过程………………………………………………………..…（13）
4.3.2.2 函数调用关系…………………………………………………..…（14）
5 优化设计计算……………………………………………………...……（14）
5.1 数据预处理……………………………………………………………（14）
5.1.1进塔物料的计算………………………………………………………（14）
5.1.2塔顶蒸气温度的计算………………………………………..……（15）
5.1.3等板高度HETP的计算………………………………………………（15）
5.1.4产品汽化潜热的计算…………………………………………...…（15）
5.1.5最小回流比的确定………………………………………………（16）
5.1.6理论板数N的求取…………………………………………………...（16）
5.1.7空塔气速的计算………………………………………………………（16）
5.2 塔径的计算……………………………………………………………（17）
5.3 填料层高度的计算……………………………………………………（17）
5.4精馏塔塔体年投资折旧费及维修费用……………………….…….（17）
5.5冷凝器年运转费用…………………………………………….…….（18）
5.5.1冷凝器冷却水最佳出口温度的确定………………………...;…（18）
5.5.2冷却水用量及冷凝器传热面积的计算……………….......…...（18）
5.5.3的计算………………………………………….……….……….…（18）
5.6再沸器年运转费用或加热蒸汽费用……………….…….…..…….（18）
5.7填料年折旧费用…………………………………….…….…..…….（19）
5.8 汽液负荷……………………………………………….….….…….….（19）
5.9 年总费用与回流比的关系…………………………….….….………（19）
6 填料塔水力学性能校核………………………………….….…………（20）
6.1 泛点率的校核………………………………………….….………….（21）
6.2核算径比………………………………………………….….…………（20）
6.3核算喷淋密度…………………………………………….….…………（21）
6.4 填料塔压降…………………………………………….….……….…（21）
7 附属设备的设计与选型………………………………….….…………（21）
7.1 塔顶冷凝器…………………………………………….…..……….…（21）
7.2 接管选型……………………………………………….….……….…（28）
7.3冷却水输送泵………………………………………….….……………（32）
7.4 填料支承结构………………………………………….….……….…（34）
7.5 填料压紧装置………………………………………….….……….…（34）
7.6 液体分布装置………………………………………….….……….…（34）
7.7 液体收集再分布装置………………………………….….……….…（34）
7.8 气体分布装置………………………………………….….……….…（34）
8 设计结果汇总（以三线表分类汇总）………………….………….…（34）
设计心得……………………………………………………….………...…（35）
参考文献……………………………………………………….……………（36）
附录一 甲醇—水系统的主要物理性质
附录二 甲醇—水汽液平衡数据
附录三 优化设计程序源代码
附录四 填料精馏塔装配示意图
摘要：本说明书介绍了甲醇—水溶液填料精馏塔的优化设计，主要内容包括了此次课程设计的计算机编程—最佳回流比的优化计算、塔的主要尺寸设计（包括塔板的板面设计、阻力损失等）、辅助设备选型、填料精馏塔图纸的绘制等若干重要环节。本文详细阐述了设计的思路，计算贯穿在整个设计中，最后得出一定条件下的最优化设计方案，并在附录中填加了优化设计的程序清单。
关键词：甲醇精馏；填料塔；优化设计
1 前言
本次课程设计任务为设计一甲醇—水溶液填料精馏塔，要求处理量： 20000（吨/年）、料液浓度：15%（wt%）、产品浓度：99.5%（wt%）、回收率： 99.9%、填料类型：鲍尔环、每年实际生产时间：7200小时/年。通过对甲醇—水填料精馏塔的优化设计，提出对于一定工艺要求的最优化方案，从而达到节能和节省费用的目的。
在化工生产中，精馏是最常用的单元操作,，是分离均相液体混合物的最有效方法之一。在化学工业中，总能耗的40%用于分离过程，而其中的95%是精馏过程消耗的因此，有必要开辟多种途径来回收利用余热，降低再沸器能耗，实现精馏节能。同时，精馏所需费用在生产装置的总投资及操作费中占了相当大的比例。当今世界对甲醇的需求量极大，而甲醇的精馏也越来越受到重视，因此甲醇的精馏的研究也越来越重要。甲醇精馏塔的优化设计无论是对节省投资，还是降低能耗，都具有非常重要的意义。
为了使填料塔的设计获得满足分离要求的最佳设计参数(如理论板数、热负荷等) 和最优操作工况(如进料位置、回流比等) ,准确地计算出全塔各处的组分浓度分布(尤其是腐蚀性组分) 、温度分布、汽液流率分布等,常采用高效填料塔成套分离技术。而且,20 世纪80 年代以来,以“高效填料及塔内件”为主要技术代表的新型填料塔成套分离工程技术在国内受到普遍重视。由于其具有高效、低阻、大通量等优点,广泛应用于化工、石化、炼油及其它工业部门的各类物系分离。进入20 世纪90 年代,高效填料塔成套分离工程技术开始向行业化、复合化、节能化、大型化方向发展,如复合塔。所谓复合塔(Compound Tray) 是指人们将塔板与填料有机地结合起来而形成的一种新型塔板。其目的在于将塔板的优点和填料的优势加以互补。此种复合塔具有效率高、通量大及压降小的性能。在国内,复合塔板已在溶剂回收、酒精、丙酮和甲醇精馏中成功应用。
实践证明，对于甲醇精馏分离能力，不能仅依靠高效能的塔内件应用于塔内以达到降低能耗，而甲醇精馏工艺流程的改进也可以明显降低能耗．
2 方案论证
2.1 精馏塔类型
本次设计采用填料塔精馏。精馏是气液两相之间的传质过程，而传质过程是由能提供气液两相充分接触的塔设备完成，并要求达到较高的传质效率。根据塔内气液接触部件的结构型式，可分为板式塔与填料塔两大类。板式塔内设置一定数量塔板，气体以鼓泡或喷射形式穿过板上液层进行质量、热量传递，气液相组成呈阶梯变化，属于逐级接触逆流操作过程。填料塔内装有一定高度的填料层，液体自塔顶填料表面下流，气体逆流而上，（也有并流向下者）与液相接触进行质量、热量传递,气液相组成沿塔高连续变化,属于微分接触操作过程[4].
工业对塔设备的主要要求：①生产能力大；②传质、传热效率高；③气流的摩擦阻力小；④操作稳定，适应性强，操作弹性大；⑤结构简单、材料耗用量少；制造安装容易，操作维修方便。此外还要求不易堵塞、耐腐蚀等。实际上，任何塔设备都难以满足上述要求，因此，设计者应根据塔型特点、物系性质、生产工艺条件、操作方式、设备投资、操作与维修费用等技术经济评价以及设计经验等因素，依矛盾的主次，综合考虑，选择适宜的塔型。
2.2 精馏压力
塔的操作压力的选择实际上是塔顶和培底温度的选取问题[5]。在塔顶产品的组成被决定以后，塔顶的温度和压力只能选定一项。精馏操作通常可在常压、减压和加压下进行。确定操作压力时，必须根据所处理的物料性质，兼顾技术上的可行性和经济上的合理性进行全面考虑。操作压力常取决于冷凝温度。一般除热敏性物料外，凡通过常压蒸馏不难实现分离要求，并能用江河水或循环水将馏出物冷凝下来的系统，都应采用常压蒸馏；对热敏性物料或混合液沸点过高的系统则宜采用减压蒸馏；对常压下馏出物的冷凝温度过低的系统，需提高塔压或采用深井水、冷冻盐水作为冷却剂；而常压下呈气态的物料必须采用加压蒸馏。例如苯乙烯常压沸点为145.2℃,而将其加热到102℃以上就会发生聚合，故苯乙烯应采用减压蒸馏；脱丙烷、丙烯塔操作压力提高到1765kPa时，冷凝温度约50℃，便可用江河水或循环水进行冷凝冷却，则运转费用减少[6]. 本次设计的对像为甲醇—水溶液物系，故采用常压精馏。
2.3 进料方式
进料可以是过冷液体(＞1)、饱和液体(＝1)、气液两相(1＞＞o)、饱和蒸气(＝o)或过热蒸汽(＜o)。不同的进料状态对塔的热负荷、塔径和所需的培板数都有影响，但进料状态主要取决于系统的前一工序的物料状态。从设计的角度来看，饱和液体进料时，精馏段和提馏段的气液流率基本相近，两段塔径可以相同以便于设汁和制造，操作上也比较容易控制。因此，如果原料为过冷液体，可考虑加设原料预热器，将料液预热至泡点。从操作费用的角度来看，对于高温精馏(即塔顶用冷却水冷凝，塔釜用低压蒸汽加热)，当／值大时宜采用较小的值；当／值较小时宜采用值较大的气液相混合进料；对于低温精馏(如塔顶用173K液态乙烯蒸发制冷，塔釜用273K丙烯蒸汽加热)，不论／值如何，均以采用较高的值为经济[6]。对于具体情况，需要综合考虑设备费用和操作费用两方面的因素。因此，本次设计采用泡点进料方式。
2.4 填料类型
填料的种类很多，根据装填方式的不同，可分为散装填料和规整填料两大类。散装填料是一个个具有一定几何形状和尺寸的颗粒体，一般以随机的方式堆积在塔内，又称为乱堆填料或颗粒填料。散装填料根据结构特点不同，又可分为环形填料、鞍形填料、环鞍形填料及球形填料等。较典型的散装填料有拉西环，鲍尔环，阶梯环，矩鞍形、弧鞍形，环矩鞍形等。规整填料是按一定的几何图形排列，整齐堆砌的填科。规整填料种类很多，根据其几何结构可分为格栅填料、波纹填料、脉冲填料等，工业上应用的规整填料绝大部分为波纹填料。波纹填料按结构分为网波纹填料和板波纹填料两大类，可用陶瓷、塑料、金同等材质制造。规整填料其特点是压降低、分离效率高，特别适用于精密精馏及真空精馏装置，为难分离物系、热敏性物系的精馏提供了有效的手段。尽管其造价高，但因性能优良仍得到了广泛的应用。
本次设计填料采用鲍尔环，方式为乱堆。鲍尔环它的结构是在普通拉西环的壁上开一层或两层长方形孔。开孔时，孔材不完全从环上断开，而是断开四边形的三条边，保留另一边，并使其开出舌状弯向环的中心，几乎在环中心对接起来。上下两层孔的位置是错开的。—般孔的面积为整个环壁的35％左右。这样．气、液体便可从一个个孔中流过，流通性能改善，对于同样的孔隙率流动阻力大为降低，环的内表面积得以充分利用。只外，由于开孔后保留的舌片向中心弯去，所以液体的分布较为均匀，改进了拉西环使液休向壁偏流的缺点。正固如此，与拉西环相比．鲍尔环具有生产能力大、阻力小、效率局、操作弹性大等优点。
2.5 加热方式
本次设计选用直接蒸加热。塔釜料液的加热方式可以是间接加热或直接加热。通常情况下，如果要求加热的温度低于180℃，一般都采用饱和水蒸气作为加热剂(也可以使用其他工序的热载体)。为了避免加热蒸汽冷凝产生的冷凝水与塔内的物料混合，塔底可设再沸器采用间接方式加热，同时管线上应安装冷凝水排除器(疏水器)，以便在防止蒸汽逸出的情况下及时地排除冷凝水。如果精馏的物料是水溶液，水为难挥发组分，则可采用直接蒸汽加热的方式，塔底只需设置简单的蒸汽分布管，不需要设置再沸器。如果塔底要求加热的温度超过了180℃左右，则应考虑采用其他的高温热源，如烟道气等。关于加热蒸汽温度的选择，应考虑经济效益，传热温差不宜选取得过大，以能够使跺脚传热维持在核状沸腾阶段为宜。过高的蒸汽温度或压力不仅不利于传热，而且还将导致设备费用和操作费用大大增加。
2.6回流比的选择
适宜的回流比是指精馏过程中设备费用与操作费用两方面之和为最低时的回流比。精馏过程的主要设备费用有精馏塔，再沸器和冷凝器，当回流比最小时，塔板数为无穷大，故设备费用最大，当回流比略大于最小回流比时，塔板数便从无穷多锐减到某一值，塔的设备费用随之锐减，当回流比继续增大时，塔板数仍随之减少，但已较缓慢。但是，由于回流比的增加，导致上升蒸气量随之增加，从而使塔径、再沸器、冷凝器等尺寸相应增大，设备费用随之上升。对特殊物系的场合，则应根据实际需要选定回流比。在进行课程设计时，也可参考同类生产的经验值选定，必要时可选若干个R值，利用吉利兰图（简捷法）[5]求出对应理论塔板，曲线，从中找出适宜操作回流比。
3 数学模型的建立
优化设计模型以填料精馏塔年总费用最低为优化目标，用菲波那契法求解单变量优化问题，回流比R为变量建立数学模型如下：
3-1
| 式中 | 填料精馏塔年总费用，元/年； | ||
| 填料精馏塔塔体年投资折旧费及维修费用，元/年； | |||
| 冷凝器年运转费用，元/年； | |||
| 再沸器年运转费用或加热蒸汽费用，元/年； | |||
| 填料年折旧费用，元/年。 |
3-2
| 式中 | 0.06 | 精馏塔年维修费用比率； | |
| 年折旧率；一般可取折旧年限5~8， =0.125~0.2； | |||
| 塔体费用，元； | |||
| 通货膨胀系数，缺乏数据时按年增5%计算。 |
3-3
其中： ，
| 式中 | 塔质量，kg；（化工系统工程基础） | ||
| 美元与人民币兑换率； | |||
| 塔高，m； | |||
| 塔壁厚，m；常压3mm，可取5mm； | |||
| 碳钢密度，； | |||
| 塔径，m； | |||
| 理论塔板数； | |||
| 等板高度，m； | |||
| 塔两端高度，m；包括塔顶空间、塔底空间、裙座（常用化工单元设备的设计）、及塔内件、进料等空间的总和。 |
本次设计采用下式计算：
| 式中 | 液体表面张力，N/m； | ||
| 液体粘度，Pa.s。 |
| 填料种类 | 填料尺寸 | ||
| #25（1英寸） | #40（1.5英寸） | #50（2英寸） | |
| 金属鲍尔环 | 1.3104 | 1.5747 | 1.8380 |
| 金属矩鞍环 | 1.3104 | 1.5499 | 1.7482 |
| 陶瓷矩鞍环 | 1.3104 | 1.5953 | 1.9029 |
3.2 冷凝器年运转费用
3-4
| 式中 | 冷却水价格，元/kg；可取0.00012~0.0004； | ||
| 冷却水用量，kg/h； | |||
| 年工作时间，h/a； | |||
| 冷凝器价格回归系数； | |||
| 冷凝器价格回归指数； | |||
| 冷凝器压力校正系数； | |||
| 冷凝器材质校正系数； | |||
| 冷凝器传热面积，。 |
直接蒸汽加热
3-5
| 式中 | 蒸汽价格，元/kg；可取0.03~0.05，常压取0.03~0.035。 |
吸收塔填料费与填料填充体积和填料类型有关，即
3-6
式中
| 填料层高度，m； | ||
| 填料单价，元/m3。 |
4.1数学模型决策变量分析
在完成同样任务要求的情况下，填料精馏塔系统的年总费用越小越好，这也是优化设计的最终目的。所以本次设计中以精馏塔系统的年总费用为目标函数，建立数学模型，精馏塔系统的年总费用包括填料精馏塔塔体年投资折旧费及维修费用、冷凝器年运转费用、再沸器年运转费用或加热蒸汽费用、填料年折旧费用。其中因变量是回流比，通过确定一最适合回流比使年总费用最小，从而达到求解。
4.2 主要工艺参数的求解
4.2.1 塔径的计算
塔径可按照下式计算
3-7
式中
| 汽相体积流量，m3/h； | ||
| 空塔气速，m/s，其中； | ||
| TD | 塔顶蒸汽的露点温度，℃； | |
| 液泛分率，无因次，可取 0.5~0.8，一般在0.7左右较好。 |
3-8
式中
| 与填料类型有关的常数，B-H常数； | ||
| 液相质量流率，kg/m2.s； | ||
| 液相粘度，Pa.s； | ||
| 干填料因子，1/m。 |
3-9
式中
| 填料因子； | ||
| 液相密度校正因子（非水系统需校正）； | ||
| 汽相质量流率，kg/m2.s。 |
（1）在填料塔中不同位置的泛点气速不同，故要对塔顶和塔底进行比较；
（2）表面张力的混合规则为：（、为两组分）；
（3）由于优化设计结果还未得到，未得到；所以在计算以上参数时，回流比可以在区间内选择一适宜值，它对优化设计结果没有影响，因为的范围较大。
4.2.2 理论塔板数的计算
4.2.2.1 相平衡关系的表示
对理想溶液，其相平衡关系为：；
对非理想溶液，其相平衡关系可以从实验数据中利用三次样条插值得到或通过回归实验数据得到相平衡关系表达式，其形式主要有：
（余国琮式），（阮奇式）
| 式中 | 、 | 分别为液相和汽相的摩尔组成； | |
| 、 | 相平衡关系回归系数。 |
4.2.2.2 的计算
（全凝器）……<
注意点：
（1）<时，更换操作线方程，注意区别直接与间接蒸汽加热时的提馏段操作线方程的不同；
（2）若相平衡关系是分段表示的，则必须判断汽相组成所在的区间来选择相平衡关系式；
（3）总理论塔板数为，不对理论板数进行圆整是为了保证回流比和塔板数一一对应的关系。
4.2.3 冷却水用量及冷凝器传热面积的计算
3-10
| 式中 | 冷凝器的热负荷，kJ/h； | ||
| 冷却水比热，； | |||
| 、 | 冷却水进出口温度，； | ||
| 冷凝器总传热系数，；1400~3000 | |||
| 冷凝器传热推动力，； | |||
| 塔顶蒸汽温度(露点)，； | |||
| 馏出物汽化潜热，kJ/kmol。 |
3-11
3-12
将式3-11和3-12代入式3-4得：
3-13
∴
4.2.4 冷凝器冷却水最佳出口温度的确定
其中，
为的解，因为为非线性方程，无法直接求解，故应采用数值方法，如牛顿迭代法（需对求导）。
∴
4.2.5 直接蒸汽加热蒸汽流量的计算
式3-5中的为加热蒸汽的流量，所以：
，kmol/h
3-14
4.3数学模型的求解
4.3.1单变量最优化方法
求单变量函数的极小值（或极大值）的方法有两类。一类是不需要计算目标函数的导数的，通常称为直接法；另一类是需要计算导数的，通常称为微分法或间接法。由于工程问题的复杂性，在进行化工单元过程最优化时所建立的目标函数往往是无法求导的，只能用不需要求导数的各种直接法搜索求解。本次优化设计采用菲波拿契（Filonacci）法.菲波拿契法是基于Filonacci数列的一种搜索法，在相同的计算次数下，它比其它方法得到的精度高，或者说在同样的精度要求下，它所需要的计算次数最少。所以是一个效率高的方法。用Filonacci法求一元函数极小值点的程序框图如图1所示。
图1 Fibonacci法的框图
4.3.2 优化设计程序框图
4.3.2.1 求解过程
4.3.2.2 函数调用关系
4.3.3 优化设计程序
本次优化设计程序采用C语言编写，具体设计程序见附录三，具体输出结果如下表．
表１ C程序输出结果
| 项目 | 结果 | 单位 |
| 最适合回流比 | 2.602938 | - |
| 最小回流比 | 2.479000 | - |
| 理论塔板数 | 24.999993 | 块 |
| 塔高 | 17.424996 | m |
| 最佳进料位置 | 9.069996 | m |
| 塔径 | 0.721729 | m |
| 塔体费用 | 169103.943408 | 元／年 |
| 冷凝器冷却水最佳出口温度 | 49.256676 | ℃ |
| 传热面积 | 43.502518 | m2 |
| 冷却水用量 | 13663.301269 | kg/h |
| 加热蒸汽用量 | 850.221372 | kg/h |
| 精馏塔塔体年投资折旧费及维修费用 | 122492.132447 | 元／年 |
| 冷凝器年运转费用 | 54358.991466 | 元／年 |
| 再沸器年运转费用，直接蒸汽加热 | 214255.785772 | 元／年 |
| 填料年折旧费 | 7270.030779 | 元／年 |
| 填料精馏塔年总费用 | 398376.940463 | 元／年 |
| 塔底釜液摩尔分数 | 0.000071 |
5.1 数据预处理
5.1.1进塔物料的计算
m2
5.1.2塔顶蒸气温度的计算
塔顶蒸气温度可由内插法求得，kg/h
由附录一：甲醇—水溶液汽液相平衡数据（摩尔组成）可查得：
t1=65.00℃ x1=0.95 ; t2=.50℃ x2=1.00
代入 得=.58℃
5.1.3等板高度HETP的计算
近似取塔平均操作压强为l01．3kPa，填料为#25金属鲍尔环，n=1.3104 故：
塔顶：由手册可得=.58℃下
=16.80mN/m =65.36mN/m =0.327 cP =0.441cP
=16.91mN/m
=0.328cP
=0.52521m
塔底：按100℃水算 由手册查得 =58.8 mN/m =0.283cP
=0.40335m
所以 =0.43m
取安全系数为1.2
所以=0.557
5.1.4产品汽化潜热的计算
由手册查得水的汽化潜热为2285kJ/kg 甲醇的汽化潜热为1101kJ/kg
所以35324.3kJ/kg
5.1.5最小回流比的确定
由常压下甲醇—水溶液的汽液相平衡数据可回归成
代入计算可得=0.3492
由得=2.479
5.1.6理论板数N的求取
由甲醇—水溶液的平衡方程与操作线方程（包括提馏段与精馏段操作线方程）进行逐板计算得出N=24.99（具体过程见附录：优化设计程序编程部分）。
5.1.7空塔气速的计算
由手册可得=.58℃下
=753.0
令
由得=0.9876
取
=0.0286
由手册查得#25金属鲍尔环的填料因子 液相密度校正因子
代入式：
可得=1.2695
取液泛分率=0.7
所以空塔气速
5.2 塔径的计算
塔径可按照式3-7计算
其中
所以=0.734m
圆整塔径取=0.8m
5.3 填料层高度的计算
填料层高度=
由表１知道最佳进料位置＝9.070m，假设的＝3.5m
所以提馏段高度＝－＝9.070－3.5＝5.57m
精馏段高度＝－5.57=8.355m
5.4精馏塔塔体年投资折旧费及维修费用
塔体材料为碳钢，取=3.5m, =0.005m, =7860
=17.20m
=1762 kg
根据式3-3得塔体费用=170741元/年
取=4.187， =0.2
代入式3-2得=185872元/年
5.5冷凝器年运转费用
5.5.1冷凝器冷却水最佳出口温度的确定
其中 ，
为的解，因为为非线性方程，无法直接求解，故应采用数值方法，如牛顿迭代法（需对求导）。
取=7200h, =487, =0.72, =1, =6.5, =.58℃, =20℃
=4.174,=1400
代入各项数据可得=49.26℃
5.5.2冷却水用量及冷凝器传热面积的计算
根据式3-10得=1725980 kJ/h
根据式3-11,3-12得=14132kg /h, =45.04
5.5.3 的计算
取=0.0002元/kg,
由式3-13得=55878元/年；
5.6再沸器年运转费用或加热蒸汽费用
进料为泡点进料，,
根据式3-14得加热蒸汽的流量=879.5 kg/h
再由式3-5得=215486元/年
5.7填料年折旧费用
取=6384元/m3 ,根据式3-6得=8788元/年
5.8 汽液负荷
5.8.1汽相负荷
下，塔顶蒸汽密度为1.15kg/m3，塔顶=31.92g/mol
5.8.2液相负荷
，回流液密度
5.9年总费用与回流比的关系
由5.4、5.5、5.6、5.7可以得出常规设计年总费用为＝4034元/年
表1 常规设计与优化设计的比较
| 回流比 | 1 | 1.00001 | 1.001 | 1.005 | 1.01 | 1.015 | 1.02 |
| 年总费用 | 9853 | 493707 | 434275 | 416856 | 409990 | 406321 | 403737 |
| 增加费用% | 63.13 | 23.93 | 9.01 | 4. | 2.92 | 1.99 | 1.35 |
| 回流比 | 1.025 | 1.03 | 1.035 | 1.038 | 1.04 | 1.05 | 1.06 |
| 年总费用 | 402071 | 401658 | 400776 | 399781 | 399799 | 398376 | 399212 |
| 增加费用% | 0.93 | 0.82 | 0.6 | 0.35 | 0.35 | 0 | 0.21 |
| 回流比 | 1.07 | 1.08 | 1.09 | 1.1 | 1.2 | 1.3 | 1.4 |
| 年总费用 | 399220 | 399921 | 400328 | 400188 | 411424 | 426002 | 441259 |
| 增加费用% | 0.21 | 0.39 | 0.49 | 0.45 | 3.28 | 6.93 | 10.76 |
| 回流比 | 1.5 | 1.6 | 1.7 | 1.8 | 1.9 | 2 | |
| 年总费用 | 457279 | 474383 | 491301 | 508465 | 525403 | 543550 | |
| 增加费用% | 14.79 | 19.08 | 23.33 | 27.63 | 31. | 36.44 |
6 填料塔水力学性能校核
6.1 泛点率校核
由结果得
所以
代入式
得=3.1146 m/s
所以泛点率=28.53%
6.2核算径比
6.3核算喷淋密度
由手册查得最小润湿率=0.08 m3/(m.h) 填料比表面积=209m2/m3
所以m3/(m.h)
实际操作中
回流液
回流液压体积
喷淋密度m3/(m.h) >
6.4 填料塔压降
横坐标
纵坐标
由文献[7]查Ekert图可得压力降为8.8mmH2O/m
7 附属设备的设计与选型
7.1 塔顶冷凝器
7.1.1 初估冷凝器传热面积
7.1.1.1 冷凝器传热量
7-1
| 式中 | 冷凝器传热量，； | ||
| 精馏段汽相流量，； | |||
| 冷凝器中汽相冷凝潜热，； | |||
| 塔顶产品流量，； | |||
| 、 | 回流比和最佳回流比。 |
7.1.1.2 冷凝器传热推动力
在冷凝器中塔顶产品在露点温度下由汽相冷凝成液相产品，故冷凝器传热属于一侧有相变化的情况，不论冷凝器是单管程还是多管程，其对数平均温度差均以逆流时的对数平均温度差计算，不必进行校正。
7-2
| 式中 | 冷凝器传热推动力，； | ||
| 冷凝器汽相（第一块塔板汽相）露点温度，； | |||
| 冷凝器中冷却水进口温度，； | |||
| 冷凝器中冷却水最佳出口温度，。 |
若，则冷凝器应装有温度补偿装置或采用浮头式换热器。
7.1.1.3 初估冷凝器传热面积
7-3
| 式中 | 冷凝器传热面积，； | ||
| 冷凝器总传热系数，。 |
其中传热系数是在优化设计中初估的，它不是冷凝器中传热系数的实验值，所以计算的传热面积也是初估值，根据初估的面积，可以从换热器系列表中选择一台适宜的换热器。
选择一台固定管板式列管换热器，其尺寸为：
公称直径D=500mm 公称压力P=1MPa
管程数N=2 管子根数nt=1
换热面积A=56.6 管长L=4.5m
管子成正三角形排列 管子规格mm
7.1.2 冷凝器选型
7.1.2.1 冷凝器传热系数的校核
由于传热系数在优化设计前预先初估，它与冷凝器传热系数的实际值必须吻合，否则，优化设计就失去了实际意义。
7-4
| 式中 | 、 | 冷凝器管外、内对流给热系数，； | |
| 、 | 冷凝器管外、内污垢热阻，； | ||
| 冷凝器管壁导热系数，； | |||
| 、 | 冷凝器管外、内径，m； | ||
| 冷凝器管壁平均直径，m； | |||
| 冷凝器管壁厚度，m。 | |||
| 冷凝器总传热系数，；。 |
、、、,由手册查得。
管外为有机物甲醇，管内为未经处理的井水，针对不同的体系可查得、；
（1）管外的计算
对水平管外冷凝的对流给热系数可用下式计算
7-5
其中
| 式中 | 冷凝液质量流量，；； | ||
| 冷凝液平均分子量； | |||
| 冷凝器管长，m； | |||
| 冷凝器总管数； | |||
| 冷凝液粘度，；（注意：粘度、密度、导热系数的混合规则） | |||
| 冷凝液密度，； | |||
| 冷凝液导热系数，。 |
假设冷凝器壁温，可求得，查得此温度下：
、、
把以上数据代入式7-5可得
求得后根据牛顿粘性定律
7-6
| 式中 | 冷凝器管外壁温度，； | ||
| 冷凝器管外表面积，。 |
代入式7-6计算得管外侧壁温=62.153，的计算值与初估值62.15接近，迭代计算成功。所以
（2）管内的计算
当、、时，管内对流给热系数可用下式计算
7-7
该式适用于低粘度流体（），本次设计对像为甲醇—水溶液，为低粘度流体，故该式可用。
当雷诺数小于10000时，应乘以校正因子。
| 式中 | 冷却剂（水）导热系数，； | ||
| 冷凝器管内径，m；即； | |||
| 冷凝器管内流速，m/s； | |||
| 冷凝器管内流体密度，； | |||
| 冷凝器管内流体粘度，； | |||
| 冷凝器管内流体的比热，； |
①的单位为，但在式（7）中应采用单位；
② 物性参数、、、等应由平均温度查表或计算，；
③流速的计算：
管内流体（冷却水）的用量：，；则流速为
7-8
式中为单管程的管子数，若对双管程共有200根管，则。
计算如下：
由手册查得此此温度下管内液体水的物性参数如下：
、、
、；
管内径
管内流体（冷却水）的用量： =13661.98
代入式7-8得=0.15 m/s
雷诺数=4024.43<10000,要进行校核。
校正因子=0.8052
把以上相关数据代入式7-7得
把， 代入7-4得
7.1.2.2 冷凝器传热面积的校核
由1.2.1可求得，从而实际所需的冷凝器传热面积<56.6，
其裕度：
所选则的换热器可用。
7.1.2.3 冷凝器管程、壳程流动阻力
（1）管程阻力损失
7-9
其中： ，
| 式中 | 每程直管压降，； | ||
| 每程局部阻力，； | |||
| 壳程数； | |||
| 每壳程的管程数； | |||
| 冷凝器管内流体密度，； | |||
| 冷凝器管内流体流动阻力系数。 |
由7.1.2.1得=0.15 m/s, ,、；
由手册查得所选换热器的=1， =2；
管内雷诺数=4024.43；
由手册查得碳钢的粗糙度,
所以相对粗糙度为
查图[7]得摩擦因素
=128.33,=33.55
把以上数据代入式7-9得
（2）壳程阻力损失
7-10
其中： ，
| 式中 | 换热器壳程内径，m； | ||
| 折流板数目；蒸汽冷凝时折流挡板间距只有480、600mm两种类型； | |||
| 冷凝器壳程当量直径，m； | |||
| 冷凝器壳程流速，m/s；；设有折流挡板时，；无折流挡板时，； | |||
| 冷凝器壳程折流挡板间距，m； | |||
| 管子外径，m； | |||
| 管子中心距，m。 |
计算如下：
所选换热器=0.025m, 管子中心距=0.032m, =500mm；
当量直径=0.020m
蒸汽冷凝时无折流挡板， =0;
=0.1158
=3.15 m/s
塔顶
=68.20 , =0.323
代入式7-10得=45.
7.2 接管选型
针对不同的流体，选择适宜的流速，由求得接管的直径，依据它选择合适的接管型号，最后校核接管中的实际流速，本次设计所选择的钢管材料均为无缝钢管。
需选型的接管主要有：
进料接管：0.4~0.8m/s，泵送1.5~2.5m/s
冷却水接管：1.0~2.5m/s
塔顶蒸汽接管：12~20m/s
塔顶产品接管：0.5~1.0m/s
塔底产品接管：0.5~1.0m/s
塔顶产品回流接管：0.2~0.5m/s，泵送1.0~2.5m/s
塔底加热蒸汽接管：20~40m/s(<295kPa)，40~60m/s(<785kPa)，80m/s(>2950kPa)
注意：计算中为体积流量，必须针对不同流体、不同状态（气体），不同组成（气体、液体）计算。计算结果汇总如表6。
7.2.1进料接管
由设计条件知=2777.78kg/m3 , =15%
由附录一运用内插法得进料温度=88.48℃,
由手册查得该温度下=727.52kg/m3, =966.29kg/ m3
计算该温度下进料混合液体的密度为
设=0.4m/s,则=51.38mm
设=0.8m/s,则=36.33mm
选取，则=0.66m/s
7.2.2冷却水接管
由前可知冷却水用量为13662kg/h,密度为998.2kg/m3
所以求得=13.68 m3/h
设=1.0m/s,则=69.59mm
设=2.5m/s,则=44.01mm
选取，则=1.94m/s
7.2.3塔顶蒸汽接管
下，塔顶蒸汽密度为1.15kg/m3
质量流量=0.4188kg/s,所以=1311.11 m3/h
设=12m/s,则=196.63mm
设=20m/s,则=152.31mm
选取，则=18.58m/s
7.2.4塔顶产品接管
塔顶产品为甲醇质量分数99.5%的甲醇――水溶液
产品，可求得=0.56 m3/h
设=0.5m/s,则=19.84mm
设=1.0m/s,则=14.03mm
选取，则=0.50m/s
7.2.5塔底产品接管
所选加热方式为直接蒸汽加热，所排出废水（含微量甲醇）的量与相等[7]，
密度按纯水算可得=3.35 m3/h
设=0.5m/s,则=48.69mm
设=1.0m/s,则=34.43mm
选取，则=0.74m/s
7.2.6塔顶产品回流接管
根据前述可得塔顶产品回流液的=1.45m3/h
设=0.2m/s,则=50.6mm
设=0.5m/s,则=32mm
选取，则=0.32m/s
7.2.7塔底加热蒸汽接管
塔底加热蒸汽为100℃水蒸汽，蒸汽用量为=1425.76m3/h
设=20m/s,则=158.83mm
设=40m/s,则=112.31mm
选取，则=32.29m/s
7.3 冷却水输送泵
7.3.1 塔高计算
7-11
| 式中 | 塔高，m； | ||
| 塔顶空间高度，m； | |||
| 理论塔板数； | |||
| 等板高度，m； | |||
| 塔内件及人孔、手孔、进料位置等空间的总高度； | |||
| 塔釜空间高度，m；保证釜液10~15min的储量； | |||
| 裙座高度，m。 |
根条据条件取=1.2m, =1.8m, , =1.3m;
塔釜空间高度要保证釜液10~15min的储量，假设15min的储量，由7.2.5有=3.35 m3/h＝0.056 m3/min
由，可得＝1.667m,取＝1.8m
代入式7-11得=19.8m
7.3.2 冷却水输送泵选型
输水泵的选型主要计算体积流量和扬程。
7-12
7-13
其中： ，，
| 式中 | 体积流量，； | ||
| 扬程，m； | |||
| 冷却水流量，kg/h； | |||
| 冷却水密度，。 |
计算如下：
冷却水温度为20℃，由手册查得此温度下水的， =1.005cP
由7.2可知冷却水接管为， =13.68 m3/h， 1.94m/s
所以由式7-12得=13.69 m3/h
雷诺数=96344；
由手册查得碳钢的粗糙度,
所以相对粗糙度为
查图[7]得摩擦因素
取管长
则
由7.3.1知＝19.8m
代入式7-13得=22.784
查手册可选泵的型号IS65-50-160，具体数据如下表：
表2 泵的选型
转速
| r/min | 流量/m3 | 扬程m | 效率 | 轴功率 | 电机功率 | NPSH |
| 2900 | 15 | 35 | 54% | 2.65 | 5.5 | 2.0 |
本次填料塔设计为散装填料，由于塔径为0.8m,选用分块栅板型支承装置，塔内设有支承圈，支承圈为的扁钢。
7.5 填料压紧装置
每层填料设置一块压紧装置，选用压紧栅板，在其下方，根据填料的规格敷设一层金属网，并将其与压紧栅极固定。
7.6 液体分布装置
在塔内回流液接管处设置一喷头式液体分布器，喷头标准为H45-0010-2．.
7.7 液体收集再分布装置
为减小壁流现象，当填料层较高时需进行分段，故需设置液体收集及再分布装置。本次填料塔设计在精馏段分段填料层之间设置一分盘式液体收集分布器，在进料位置设置一遮板式液体收集再分布器．
7.8 气体分布装置
本次填料塔设计在塔釜内设置一孔管式气体分布装置，其能有效的对气体进行分布。
填料精馏塔装配示意图见附录４．
8 设计结果汇总
表3 工艺参数表
| 参数 | 数值 | 单位 | 参数 | 数值 | 单位 |
| 处理量 | 144.03 | Kmol/h | 产品汽化潜热 | 35324.3 | kJ/kmol |
| 进料浓度 | 0.0903 | 摩尔分率 | 塔顶蒸汽温度 | ．58 | ℃ |
| 产品浓度 | 0．9911 | 摩尔分率 | 塔顶产品流率 | 13．11 | Kmol/h |
| 回收率 | 99．9% | 等板高度HETP | 0．557 | m | |
| 实际生产时间 | 7200 | 小时/年 | 空塔气速 | 0．8887 | m/s |
| 项目 | 费用 / 元/年 | 项目 | 费用/ 元/年 |
| 精馏塔塔体设备费 | 169104 | 冷凝器设备费 | 54358 |
| 精馏塔塔体年投资折旧费用及维修费 | 122492 | 冷却水费用 | 19675 |
| 填料年折旧费用 | 7270 | ||
| 加热蒸汽费用 | 214256 | ||
| 参数 | 参数值 | 单位 |
| 塔材料 | 碳钢 | - |
| 塔材料密度 | 7800 | kg/m3 |
| 塔壁厚度 | 5 | mm |
| 塔径 | 0.8 | m |
| 塔高 | 20234 | m |
| 填料类型 | 鲍尔环# | - |
| 填料比表面 | 209 | m 2/ m 3 |
| 填料层高度 | 13．925 | m |
| 精馏段填料层高 | 8.355 | m |
| 精馏段填料层分层数 | 2 | - |
| 精馏段填料层第一层高度 | 4 | m |
| 精馏段填料层第二层高度 | 4.355 | m |
| 提馏段填料层高 | 5.570 | m |
| 提馏段填料层分层数 | 1 | - |
| 填料压降 | 8.8 | mmH2O/m |
| 接管 | 物流 | 型号 | 流量 m3/h | 流速 m/s | 适宜流速范围m/s |
| 进料管 | 15%wt甲醇—水溶液 | 2.99 | 0.66 | 0.5~0.8 | |
| 塔顶液相回流管 | 1.45 | 0.32 | 0.2~0.5 | ||
| 塔顶蒸汽管 | 99.5%甲醇蒸汽 | 1311.09 | 18.58 | 12~20 | |
| 塔顶产品管 | 0.56 | 0.50 | 0.5~1.0 | ||
| 冷却水输送管 | 冷却水 | 13.68 | 1.94 | 1.0~2.5 | |
| 塔底残液管 | 3.35 | 0.74 | 0.5~1.0 | ||
| 塔底蒸汽管 | 1425.76 | 32.29 | 20~40 |
本次课程设计为期两周,虽然时间有点紧,做的也很辛苦，不过我的收获确实不小．
首先，我对填料精馏有了一个系统的认识，同时也对以前学过的化工原理相关知识做了一个系统的复习．虽然以前在化工原理课程中对填料塔有一个初步的认识，但所了解的内容真的很少．在这次课程设计查找资料的过程中，我清楚的认识了填料塔的设计包括了节能效果最好的塔体的设计，塔内件、冷凝器、离心泵的选择等．对于一些填料的类型，规格，效果等方面也有了一定的了解．同时了解了国内外对于填料塔的一些最新进展．
其次，我对于精馏的整个优化设计流程有了大概的了解，明白了计算机在化工过程设计的应用前景不可限量．本次设计要求用C语言编写程序，由于有一年多没有看过C语言，而且一开始对于整个优化设计的流程以及一些参数的算法没有了解清楚，以至编程开始的第一天速度很慢．通过在宿舍一晚上时间的加班，第二天终于对整个程序的设计过程有了一个比较明确的了解．在老师和同学的帮助下，程序编辑运行成功，所得出的结果与老师提供的查询结果较吻合．我也彻底明白，掌握并精通一门计算机编程语言对以后步入社会是多么的重要．
再次，这次课程设计中要求要作填料塔的装配示意图．最后一周的时间白天基本就是泡在教室画图，虽然以前学过工程制图而且自以为学得还不错，但是在画图的过程中还是觉得很痛苦．对于一些画图的基本知识还算印象深刻，但是有一些细节东西还是没有一下子注意到，最后还是通过老师、同学和自己查的一些有关化工设备制图方法的资料才画好了图．
最后，谢谢老师无私的指导和同学们的帮助，使我在设计过程中的问题得到解答．
参考文献：
[1]赵静妮.填料塔技术的现状与发展趋势.兰州石化职业技术学院学报.2002,2(3) : 25-28
[2]王梦华.精馏过程节能技术探讨.齐鲁石油化工，2003,31(4):324-326
[3]陈 强,王树楹.塔器技术的发展现状与展望. 现代化工,1997 ,17 (11) :19.
[4]王明辉.化工单元过程课程设计. 北京：化学工业出版社,2002
[5]刘雪缓,汤景凝.化工原理课程设计.
[6] 柴诚敬，刘国维，李阿娜．化工原理课程设计．天津：天津科学技术出版社，1994
[7]谭天恩,麦本熙,丁惠华.化工原理下册.北京：化学工业出版社,1998
[8]陈敏恒,从德滋,方图南.化工原理.下册.北京：化学工业出版社,1993
[9]谭天恩,麦本熙,丁惠华.化工原理上册.北京：化学工业出版社,1990.6
[10]王志魁.化工原理.北京：化学工为出版社,2004
[11]乔敬玲.甲醇精馏塔的设计.太原科技.1995,4:16-17
附录一 甲醇—水系统的主要物理性质
附录二 甲醇—水汽液平衡数据(摩尔组成)
| t | x | y | t | x | y |
| 100.00 | 0.00 | 0.000 | 75.30 | 0.40 | 0.729 |
| 96.40 | 0.02 | 0.134 | 73.10 | 0.50 | 0.779 |
| 93.50 | 0.04 | 0.234 | 71.20 | 0.60 | 0.825 |
| 91.20 | 0.06 | 0.304 | 69.30 | 0.70 | 0.870 |
| .30 | 0.08 | 0.365 | 67.60 | 0.80 | 0.915 |
| 87.70 | 0.10 | 0.418 | 66.00 | 0.90 | 0.958 |
| 84.40 | 0.15 | 0.517 | 65.00 | 0.95 | 0.979 |
| 81.70 | 0.20 | 0.579 | .50 | 1.00 | 1.000 |
| 78.00 | 0.30 | 0.665 |
附录三 优化设计程序源代码
*填料精馏塔优化设计程序*/
#include #include #define XF 0.0903 /*进料组成*/ #define F 144.03 /*进料量 Kmol/h*/ #define Q 1 /*进料热状态*/ #define D 13.11 /*塔顶产品流率 Kmol/h*/ #define XD 0.9911 /*塔顶产品组成 摩尔分率*/ #define TD .58 /*塔顶蒸汽温度 ℃*/ #define RD 35324.3 /*产品汽化潜热 kJ/kmol*/ #define Cpa 6384 /*填料单价 元/m3*/ #define HETP 0.557 /*等板高度HETP m*/ #define U 0.8887 /*空塔气速 m/s*/ #define HA 3.5 /*塔体两端空间高度 m*/ #define Ps 7860 /*塔体材料密度 kg/m3 */ #define B 0.005 /*塔体壁厚 mm*/ #define F1 1 /*冷凝器压力校正系数*/ #define F2 6.5 /*冷凝器材质校正系数*/ #define AA 487 /*冷凝器回归系数*/ #define BB 0.72 /*冷凝器回归指数*/ #define KD 1400 /*冷凝器传热系数 kJ/(m2. ℃.h)*/ #define T1 20 /*冷却水进口温度 ℃*/ #define CW 0.0002 /*冷却水价格 元/kg*/ #define CP 4.174 /*冷却水比热 kJ/(kg. ℃)*/ #define Cz 0.035 /*蒸汽价格 元/kg*/ #define Time 7200 /*年处理时间 h/年*/ #define C 7.6 /*美元汇率*/ #define Fc 0.2 /*年折旧率*/ #define Rmin 2.479 /*最小回流比*/ #define MSI 2.786 /*通货膨胀系数*/ /*函数列表*/ /*塔径DT*/ double dt(double r) {double DT; DT=sqrt(D*22.4*(273.15+TD)*(r+1)/(3600*0.785*U*273.15)); return(DT); } /*平衡关系曲线函数*/ double p(double y) {double x,a,b; b=y/(1-y); if(y>0 && y<1.00 ) a=pow(b/3.3874,1/0.7977); x=a/(1+a); return(x); } /*塔底釜液摩尔分数*/ double xw(double r) {double XW; XW=(F*XF-D*XD)/(r*D+F); return(XW); } /*精馏段操作线*/ double y1(double x,double r) {double Y1; Y1=r/(r+1)*x+XD/(r+1); return(Y1); } /*提馏段操作线*/ double y2(double x,double r) {double Y2,w,v0,XW; w=r*D+F; v0=(r+1)*D; XW=xw(r); Y2=w*x/v0-w*XW/v0; return(Y2); } /*理论塔板数*/ double n(double r) {double y[100],x[100],N,i=1,XW; y[1]=XD;XW=xw(r); x[1]=p(y[1]); while(x[i]>XW) {if(x[i]>XF) y[i+1]=y1(x[i],r); else y[i+1]=y2(x[i],r); i++; x[i]=p(y[i]); } N=i-1+(x[i-1]-XW)/(x[i-1]-x[i]); return(N); } /*总塔高H*/ double h(double r) {double H,N=n(r); H=N*HETP+HA; return(H); } /*塔质量Ws*/ double ws(double r) {double H,N,DT,ws; DT=dt(r); N=n(r); H=N*HETP+HA; ws=3.14*DT*(H+0.8116*DT)*B*Ps; return(ws); } /*塔体费用CH*/ double ch(double r) {double CH,H,WS,DT; DT=dt(r); H=h(r); WS=ws(r); CH=C*exp(6.95+0.1808*log(WS)+0.02468*log(WS)*log(WS)+0.0158*H/DT); return(CH); } /* 精馏塔塔体年投资折旧费及维修费用J1*/ double j1(double r) {double CH,J1; CH=ch(r); J1=MSI*(Fc+0.06)*CH; return(J1); } /*冷凝器冷却水最佳出口温度t2opt 的确定*/ double t2(double r) {double T2,cd,tm=T1,tm0,f,f1; cd=1.3*MSI*AA*BB*F1*F2*Fc*pow((r+1)*D*RD,BB-1)/(pow(KD,BB)*pow(TD- T1,BB-1)); do {tm0=tm; f=-CW*Time/CP+cd*pow((tm0-1)/(tm0*log(tm0)),(1-BB))*(tm0-1-log(tm0)); f1=cd*pow((tm0-1)/(tm0*log(tm0)),2-BB)*((BB-1)*pow(tm0-1-log(tm0),2)/pow(tm0-1,2)+log(tm0)); tm=tm0-f/f1; } while(fabs(tm-tm0)>=1e-6); T2=(TD-T1-TD*tm)/(-tm); return(T2); } /*冷却水用量W*/ double w(double r) {double W,T2; T2=t2(r); W=(r+1)*D*RD/(CP*(T2-T1)); return(W); } /*传热面积AD*/ double ad(double r) {double T2,AD; T2=t2(r); AD=(r+1)*D*RD*log((TD-T1)/(TD-T2))/(KD*(T2-T1)); return(AD); } /*冷凝器年运转费用J2*/ double j2(double r) {double W,AD,J2; W=w(r); AD=ad(r); J2=CW*W*Time+1.3*MSI*AA*F1*F2*pow(AD,BB)*Fc; return(J2); } /*加热蒸汽的流量Z*/ double z(double r) {double Z; Z=((r+1)*D-(1-Q)*F)*18; return(Z); } /*再沸器年运转费用，直接蒸汽加热J3*/ double j3(double r) {double J3,Z; Z=z(r); J3=Z*Cz*Time; return(J3); } /*填料年折旧费J4*/ double j4(double r) {double J4,N,DT; N=n(r); DT=dt(r); J4=(3.14/4)*HETP*Cpa*Fc*DT*DT*N; return(J4); } /*填料精馏塔年总费用J*/ double j(double r) {double J1,J2,J3,J4,J; J1=j1(r); J2=j2(r); J3=j3(r); J4=j4(r); J=J1+J2+J3+J4; return(J); } /*进料位置高度Hf*/ double hf(double r) {double HF,y[100],x[100],N1,i=1,H; H=h(r); y[1]=XD; x[1]=p(y[1]); while(x[i]>XF) {y[i+1]=y1(x[i],r); i++; x[i]=p(y[i]); } N1=i;HF=H-N1*HETP; return(HF) ; } /*菲波拿契法求解最佳回流比及其它相关数据*/ main() { double a=1.001*Rmin,b=2*Rmin,L[100],e=1e-6,u1,u2,f1,f2,m,i=1,q,k,r,N,H,HF,DT,CH,T2,AD,W,Z,J1,J2,J3,J4,J,XW; L[1]=1;L[2]=1;L[3]=2; while(L[i+2]<(b-a)/e) {i++;L[i+2]=L[i+1]+L[i];} u1=a+(b-a)*L[i]/L[i+2]; f1=j(u1); q=i+2; k=1; m=0; do {if(m==0){u2=a+(b-a)*L[q-k]/L[q-k+1];f2=j(u2);} else {u1=a+(b-a)*L[q-k-1]/L[q-k+1];f1=j(u1);} if(f1 k=k+1; } while(k!=q-1); r=(a+b)/2; N=n(r); H=h(r); HF=hf(r); DT=dt(r); CH=ch(r); T2=t2(r); AD=ad(r); W=w(r); Z=z(r); J1=j1(r); J2=j2(r); J3=j3(r); J4=j4(r); J=j(r); XW=xw(r); printf("最适合回流比R=%f\\n",r); printf("最小回流比Rmin=%f\\n",Rmin); printf("理论塔板数N=%f\\n",N); printf("塔高H=%f\\n",H); printf("最佳进料位置hf=%f\\n",HF); printf("塔径DT=%f\\n",DT); printf("塔体费用CH=%f\\n",CH); printf("冷凝器冷却水最佳出口温度t2opt=%f\\n",T2); printf("传热面积AD=%f\\n",AD); printf("冷却水用量W=%f\\n",W); printf("加热蒸汽用量Z=%f\\n",Z); printf("精馏塔塔体年投资折旧费及维修费用J1opt=%f\\n",J1); printf("冷凝器年运转费用J2=%f\\n",J2); printf("再沸器年运转费用，直接蒸汽加热J3=%f\\n",J3); printf("填料年折旧费J4=%f\\n",J4); printf("填料精馏塔年总费用J=%f\\n",J); printf("塔底釜液摩尔分数XW=%f\\n",XW); }
Copyright © 2019-2025 51dongshi.net 版权所有
违法及侵权请联系：TEL:177 7030 7066 E-MAIL:11247931@qq.com 本站由北京市万商天勤律师事务所王兴未律师提供法律服务
