闪蒸过程的数学模型与控制过程仿真分析报告.pptVIP

） 题 目:闪蒸过程的数学模型与控制过程仿真 任务书 查找相关资料，熟悉Mathlab语言及Simulink工具箱的使用；熟悉动态模型建立方法，并建立闪蒸器的动态模型；用Simulink设计闪蒸器动态仿真模型，并进行调试；设计控制方案，并在Simulink中进行仿真调试；分析不同控制参数对控制过程的影响；整理资料，撰写论文，并提交论文打印文稿、仿真模型软件及符号说明清单；提交外文献翻译资料。 3.1 闪蒸简介 闪蒸是连续单级蒸馏过程。该过程使进料混合物部分汽化或部分冷凝得到含易挥发组分较多的蒸汽和含难挥发组分较多的液体。液体进料在一定压力下被加热，通过阀门绝热闪蒸到较低压力，在闪蒸罐内分离出气体。如果省略阀门，低压液体在加热器中被加热部分汽化后，在闪蒸罐内分两相。与之相反，如果设计合理，则离开闪蒸罐的汽液两相处于平衡状态。 3.2 平衡闪蒸过程方程 (4)归一化方程 变量数：3C+8个（F,TF ,pF ,T,p ,V,L,Q，zi ,yi ,xi） 方程总数：2C+3个 所以，模拟求解自由度为：(3C+8)-(2C+3)= C＋5 即，需规定变量数：C＋5个 其中进料变量数：C+3个（F,TF ,pF ,zi ） 表3-1 闪蒸计算类型 3.3 物性的计算方程 闪蒸模型的方程已全部列出，下面讨论物性的计算方程。一个组成为xi的多组分液体混合物，在一定压力P下将其加热至刚形成微小气泡时的温度称为该液体混合物在该压力下的泡点温度。将此混合物继续加热，升温到使其液相的最后一滴液体消失时的温度称为该混合物在该压力下的露点温度。当混合物处于泡点或露点温度下，其气液两相互成平衡状态。 平衡常数定义方程 3.4 闪蒸模型的解法 对于理想溶液或接近理想溶液，当进行平衡闪蒸时，根据已知条件的不同，闪蒸过程的操作型计算可分为多种不同的情况，在本课题中仅已一种情况为例。已知进料流量、进料组成、进料温度、闪蒸操作压力求闪蒸温度、气相组成、液相组成和流量。其计算方法是先假设闪蒸温度T，然后再根据进出物料焓相等的原则来校正T，直到T不再变化为止 。 汽化率目标式 迭代方程为 在进行闪蒸计算前，应首先判断进料混合物在指定的温度和压力下是否处于两相区，判据如下 温度计算式 绝热闪蒸的计算过程如图 3.5 闪蒸计算举例 已知液相闪蒸进料组成为：甲烷30%（摩尔分数，下同），正戊烷0.7%，；进料流率200kmol/h,进料温度22度。已知闪蒸罐操作压力110.84kPa,求闪蒸温度、汽相分率、汽液相组成和流率。 假设该物系为理想体系，各组分饱和蒸汽压的安托因方程为：甲烷 lnp0=15.2243-897.84/(T-7.16);正戊烷 ln p0=15.8333-2477.07/（T-39.49）；。p0的单位为mmHg,T单位为K。 液体比热容为： 甲烷CPV =34.33+0.05472T+3.66345*10-6T2-1.10113*10-8T3 ； 正戊烷CPV =114.93+0.34114T-1.8997*10-4T2+4.22867*10-8T3； 单位为 kJ/(mol·0C) 在Matlab中计算结果如下 闪蒸温度T= 5.28520C 汽化率 e=0.3050 液相流率：L= 139.0019 kmol/h 组成为：x(1)= 0.0064 x(2)= 0.9937 汽相流率：V=60.9981 组成为：y(1)= 0.9691 y(2)= 0.0307 通过计算发现，由于两组分物性参数相差较大，甲烷在常温下为气体，而正戊烷沸点达四十度以上，从而是分离比较完全。 4.1 动态平衡方程 4.2 仿真模型 状态空间法深入到动态空间内部，采用状态空间这种内部描述取代经典法的传递函数那种外部输入——输出系统描述。状态空间法可同时适应单输入单输出和多输入多输出。在本系统中有多输入多输出，因此选用状态方程来描述控制特性。 状态变量为 x1=p，x2=h，x3=T,x4=x;输入变量为u1=F,u2=L,u3=V,u4=TF。输出变量为x1=p,x2=h则状态方程为 : 输出方程描述系统的输出变量和状态变量、输入变量的变换关系。依工艺要求，确定闪蒸压力p与液位高度h为输出变量即:y1=p=x1y2=h=x2写成矩阵形式有 5.1 闪蒸过程的控制方案 在实际生产过程中，对闪蒸的控制方案一般采用如图所示的压力控制与液位控制即可。 控制系统方块图 5.2 建立S函数 本课题所建立的对象状态函数模型采用S-函数方法建立模块仿真过程。对于闪蒸过程的状态方程引入控制模块在此选用S- Function（system Function）,即S函数。S函数是system Function的简称，用它来写自己的Si