欢迎订阅《FPGA学习入门100例教程》、《MATLAB学习入门100例教程》

一、理论基础

二、核心程序

三、测试结果

一、理论基础

原料物性表如下：

表3-1原料物理性质表
组分	质量流率	质量分数	摩尔流率	摩尔分数	分子量	常压沸点
丙烷	472.938	7.567	10.727	11.565	44.04	-41.85
异丁烷	213.938	3.423	3.681	3.969	58.12	-11.57
正丁烷	819.500	13.112	14.100	15.203	58.12	-0.35
异戊烷	1655.500	26.488	22.945	24.740	72.15	27.002
正戊烷	2420.188	38.723	33.544	36.167	72.15	36.224
异己烷	415.000	6.640	4.816	5.192	86.18	68.1
正己烷	252.938	4.047	2.935	3.165	86.18	68.7
合计	6250.000	100.000	92.747	100.000

表5-1 气相临界参数表
	Tc	Pc	Tcm	Pcm	Trm	Prm	Z
精馏段	369.82	4.25	138.94	1.60	0.83	0.22	0.85
	407.13	3.65	52.49	0.47
	425.16	3.80	208.91	1.87
	460.39	3.38	1.85	0.01
合计			402.19	3.95
提馏段	425.16	3.80	0.47	0.00	0.84	0.27	0.81
	460.39	3.38	163.73	1.20
	469.60	3.37	245.38	1.76
	507.35	2.97	38.06	0.22
	433.75	3.20	19.83	0.15
合计			467.47	3.34

第1步：确定回流罐压力Tb

塔顶冷凝器采用水作冷剂入口温度在25-30℃之间，可以估算出回流罐温度Tb为45℃。这个部分直接将Tb作为参数进行设置。

对应代码：

第2步：回流罐压力

由温度Tb查图可得各组分在该温度下的饱和蒸汽压，再根据泡点方程P=∑Pis*Xi可求得回流罐压力：

其中Pi和Xi的数据是已知数据，论文中数据有缺失，请自行补齐。

对应代码：

第3步：确定塔顶压力Pd、塔顶温度Td

Pd由公式Pd=Pb+△P计算得到，其中管线阻力△P=0.015MPa

Td则有安托万方程编程迭代：

对应代码：

第4步：确定塔釜压力Pw、塔釜温度Tw

填料塔取经验值精馏段压降0.007MPa，P进=Pd+△P =0.879MPa，Pw=P进+△P=0.886MPa.

即根据公式可得：

对应代码：

第5步：计算最小理论板数Nm

计算最小理论板数Nm

其中D的XL，XH数据为《表3-2T101清晰分割物料平衡表》中的塔顶馏出液的摩尔分数的最大值和最小值，W的XL，XH为《表3-2T101清晰分割物料平衡表》中的塔底馏出液的摩尔分数的最大值和最小值。

对应代码：

第6步：计算最小回流比

这里通过取q为计算得到，然后计算得到Rm。

对应代码：

第7步：实际回流比R和理论板数N，(R+1)N～R/Rm曲线

根据数据，绘制(R+1)N～R/Rm曲线，从而确定实际回流比R和理论板数N。

对应代码：

第8步：计算N精，N提

对应代码：

二、核心程序

clc;
clear;
close all;
warning off;%调用初始化参数
parameters;%model
[Pb,Pd,Pj,Pw,Tb,Td,Tj,Tw,Rm,Rs,Nm,Ns,Nj,Nt,Np]=func_models();fprintf('******************T-301******************\n\n');
% 回流罐P（MPa）
fprintf('回流罐P = %4.4f\n\n',Pb);
% 塔顶P（MPa）
fprintf('塔顶P   = %4.4f\n\n',Pd);
% 进料P（MPa）
fprintf('进料P   = %4.4f\n\n',Pj);
% 塔釜P（MPa）
fprintf('塔釜P   = %4.4f\n\n',Pw);% 回流罐T（℃）
fprintf('回流罐T = %4.4f\n\n',Tb);
% 塔顶T（℃）
fprintf('塔顶T   = %4.4f\n\n',Td);
% 进料T（℃）
fprintf('进料T   = %4.4f\n\n',Tj);
% 塔釜T（℃）
fprintf('塔釜T   = %4.4f\n\n',Tw);% Rmin
fprintf('Rmin    = %4.4f\n\n',Rm);
% R
fprintf('R       = %4.4f\n\n',Rs);
% Nmin
fprintf('Nmin    = %4.4f\n\n',Nm);
% N
fprintf('N       = %4.4f\n\n',Ns);
% N精
fprintf('N精     = %4.4f\n\n',Nj);
% N提
fprintf('N提     = %4.4f\n\n',Nt);
% N实际
fprintf('N实际   = %4.4f\n\n',Np);

三、测试结果

代码仿真结果如下所示：

T201的模型数据反复验算了一下，你提供的数据出路太大，没法验算，对于T301，我们得到如下的仿真结果：

A16-11

基于matlab的精馏塔作业模拟仿真相关推荐

【基于MATLAB的火灾疏散模拟仿真】——安全隐患提前发现，疏散方案优化
[基于MATLAB的火灾疏散模拟仿真]--安全隐患提前发现,疏散方案优化随着城市化进程的不断加速,人口密度越来越大,特别是在高层建筑中.万一发生火灾,往往会对人的生命和财产造成严重损失.因此,火灾疏 ...
【数学建模】基于matlab GUI平行停车模拟仿真【含Matlab源码 1877期】
一.平行停车模拟仿真简介近年来, 随着我国经济的快速发展, 机动车保有量也迅速增长.截至2012年底, 全国机动车保有量已达2.4亿辆, 城市"停车难"的问题日趋严重.统计结果表 ...
基于matlab的杨氏干涉模拟仿真
目录 1.算法概述 2.仿真效果 3.MATLAB仿真源码 1.算法概述杨氏双缝干涉实验是物理学中最著名的实验之一.这个实验通过展示光从点光源到干涉图样的变化,揭示了光的波动特性.杨氏实验的结果可以 ...
基于matlab的杨氏双缝干涉模拟仿真+GUI界面
目录 1.算法仿真效果 2.MATLAB源码 3.算法概述 4.部分参考文献 1.算法仿真效果 matlab2022a仿真结果如下: 2.MATLAB源码 %********************* ...
基于matlab的航迹发生器模拟,仿真输出经纬度、高度、俯仰、航向角等
目录 1.算法仿真效果 2.MATLAB源码 3.算法概述 4.部分参考文献 1.算法仿真效果 matlab2022a仿真结果如下: <
matlab模拟Fraunhofer衍射,基于Matlab的夫琅禾费衍射光学仿真.doc
基于Matlab的夫琅禾费衍射光学仿真基于Matlab的夫琅禾费衍射光学仿真摘要计算机仿真技术是以多种学科和理论为基础,以计算机及其相应的软件为工具,通过虚拟试验的方法来分析和解决问题的一门综合 ...
matlab仿真炉温控制,基于MATLAB的炉温控制系统的仿真.doc
PAGE 11 PAGE 11 控制系统仿真课程大作业题目: 基于MATLAB的炉温控制系统的仿真院系名称: 电气工程学院专业班级: 自动F0904 学生姓名: 学号: 指导教师: 教师职称: ...
基于matlab的msk仿真,基于MATLAB的MSK系统的仿真
基于MATLAB的MSK系统的仿真Tag内容描述: 1.基于MATLAB的MSK系统的仿真设计1绪论11数字通信的发展通信按照传统的理解就是信息的传输与交换,为了传递消息,各种消息需要转换成电信号,消 ...
tcsc工作原理matlab仿真,基于Matlab的TCSC建模与仿真研究.doc
基于Matlab的TCSC建模与仿真研究基于Matlab的TCSC建模与仿真研究第17卷第5期 2006年1O月巾原T学院 JOURNALOFZHONGYUANINSTIT[ITEOFTECHN ...

基于matlab的精馏塔作业模拟仿真

一、理论基础

二、核心程序

三、测试结果

基于matlab的精馏塔作业模拟仿真相关推荐

最新文章

热门文章