基于Matlab的精馏稳态模拟.pdfVIP

第45卷第l期 厦门大学学报(自然科学版) V01．45No．1 2006年1月 ofXiamen Journal Science) Jan．2006 University(Natural 基于Matlab的精馏稳态模拟 吴松涛，江青茵。，曹志凯，罗江水，余文武 (厦门大学化学化学化T学院化学工程与生物_72程系，福建厦门361005) 摘要：以Matlab为平台，开发面向多组分复杂精馏塔和吸收塔系统的精馏稳态模拟计算软件包XD-SOD．软件包包台了 泡点法、露点_I击、流率加和法和新松弛法等多种稳态精僖模拟计算算法．对4种不同的分离物累，采用不同的模拟算法进 行仿真验证．仿真结果表明，软件包能够有效地应用于精馏、吸收等分离过程和复杂的催化裂化分馏塔系统的稳态模拟。 并具有数值稳定性好；收敛快，适用面广，易于改进等特点． 关键词：Matlab,精馏-分馏塔}稳态仿真fXD-SOD 018 中图分类号：TQ015．9}TQ 文献标识码：A 精馏是化工生产中最主要的传质传热过程．精馏 复杂精馏塔或吸收塔系统．通用塔可具有任意板数，并 的稳态模拟，既是精馏塔设计和优化的重要手段，也是 可在任意层极上设置进料口：侧线气液相出料口以及 精馏过程动态模拟的基础．有关精馏稳态模拟的算法。 加热或冷却设备，因此灵活性很大，适用面广．本文假 迄今已发展得比较成熟，概括起来主要可以分成矩阵 定精馏塔内元化学反应，设塔内组分数为c，总级数为 法(如泡点法、露点法、流率加和法)、逐板计算法、不稳 N，且第一级为冷凝器，第N级为再沸器． 定方程法(如新松弛法)以及结合了前面几种算法长处 匿绕图1所示的通用模型塔任一平衡级j(1≤j 的联合算法等o]，表1给出了算法的优缺点比较．由于 描述精馏模型的MESH方程组庞大和高度非线性，很方程组中含有高度非线性函数．如相平衡常数H。一 难获得解析解，需用迭代计算逐步逼近求解，另外，由 于求解过程通常还包括各种气液平衡数据的计算，涉 及到多种矩阵运算和数值算法，如求解三对角矩阵的 往是一庞大的高度非线性方程组，要获得分析解是不 追赶法o]、求解热力学性质和物性时所要用到的插值 可能的，只能采用迭代的方法或其他数值法逐步逼近 和参数拟合等，因此，必须借助专业的计算机软件系统 求解．严格算法的核心就是联立求解这些非线性基本 才能实现精馏系统的严格计算．目前国外已开发了一 方程组，各种算法的不同之处仅在于联立求解这些基 些有关精馏计算的商业软件，比较著名的有ASPEN、本方程组所采用的方法和步骤互异[I]． Pro II等系列软件，但引进价格十分昂贵．而另一方 1．2算法库 面．以矩阵或数组为基本数据单元的高效的工程计算 语言——Matlab语言的出现，降低了在精馏系统模拟 序包(见表2)，各子程序包中还提供多种选择以保证 软件的开发难度D]．因此，本文采用Matlab开发精馏 算法的适用性和易于改进． 系统稳态计算软件包XD-SOD，并对不同分离物系进 主要包含的算法可分成以下三大类： 行稳态仿真． (1)矩阵法：泡点法(BP Method、SR—YMeth— Method)、流率加和法(SR—X I精馏稳态模拟软件包XD-SODod)； (2)不稳定方