脱丙烷塔计算

新疆工程学院毕业论文（设计）2010 届题目五彩湾煤生烃潜力的研究专业应用化工技术学生姓名张营娣学号2010231422小组成员指导教师马燕老师完成日期2013-4-11新疆工业高等专科学校教务处印制新疆工程学院毕业论文（设计）任务书班级应化10-5（3）班专业应用化工技术姓名张营娣日期2013-4-111、论文（设计）题目：五彩湾煤生烃潜力的研究2、论文（设计）要求：（1）学生应在教师指导下按时完成所规定的内容和工作量，最好是独立完成。

（2）选题有一定的理论意义与实践价值，必须与所学专业相关。

（3）主题明确，思路清晰。

（4）文献工作扎实，能够较为全面地反映论文研究领域内的成果及其最新进展。

（5）格式规范，严格按系部制定的论文格式模板调整格式。

（6）所有学生必须在4月11日之前交论文初稿。

3、论文（设计）日期：任务下达日期 2013.3.5完成日期 2013.4.114、指导教师签字：新疆工程学院毕业论文（设计）成绩评定报告毕业论文答辩及综合成绩引言乙烯装置顺序分离流程中，最初均采用单塔脱丙烷，脱丙烷塔进料为碳三和碳四以上馏分，现在都用于乙烯装置双塔脱丙烷塔。

工艺对乙烯装置脱丙烷塔操作的基本要求是希望塔内能进行传质过程。

塔顶轻关键组分和塔底重能达到规定的分离纯度。

尽量提高产品的回收率，已获得较高的产量：尽量节约能源，使精馏过程中消耗的能源最少。

为此脱丙烷塔的自动控制也必须满足质量指标，物料平衡及余数条件等制要求。

脱丙烷塔的主要任务是切割C3和C4混合馏分，混合液进入精馏塔内（主要含C3和C4）进料为气态混合物.进料混合馏分经过脱丙烷塔切割分离，塔顶馏分被冷凝器冷凝后送至回流管中影响脱丙烷精馏操作因素有：进料量、成分、进料温度、再沸器加热量。

塔内蒸汽上升速度、回流量、塔顶底的采出量。

可操作变量有进料流量、塔底采出流量及再沸器加热脱丙烷塔所处环境为甲级防暴区域，工艺介质多为混合物、沸点低、易挥发、易爆生产装置处于露天低压导风向由西向东，冬夏季温差较大。

吹脱塔设计计算公式吹脱塔是工业用途中常用的应用设备，采用质疏机原理，可有效脱离污染物，去除气体、液体十分有效。

一、吹脱塔的设计原理：1. 质疏机原理：在吹脱塔内，通入的气体流动压力大，这使得污染物带有质量，无法直接结合气体而被带走，最终把污染物收集在吹脱塔内，达到净化气体的效果。

2. 内层受压力较高：在吹脱塔内，污染物被向上运动，靠近内壁会受到高压空气和抛物线型流动的作用，反复颠簸，影响污染物的同化，使其分布在塔内平衡状态，也被称为内层抛物线重洗。

3. 外层受重力较高：污染物被向上运动，向外部围墙时，受到重力层的作用，被吹脱到塔内上层被收集，也被称为外层重力重洗。

二、吹脱塔的计算公式：1. 吹脱塔流量计算公式：吹脱塔流量Q＝μ×V×δ×h2gt——μ是塔布赖特数，V是实际流动速度，δ是塔内温度的密度比，h是塔的液体高度，g是重力加速度。

2. 吹脱塔单位面积承受的动压和静压计算公式：动压Pd＝ρ×h2gt——ρ是介质的密度，h是塔的高度，g是塔的重力加速度。

静压Ps＝ρ×g×h，ρ为塔内介质的密度，g是塔的重力加速度，h为塔的液体高度。

3. 吹脱塔容量计算公式：吹脱塔容量V=Q/D——Q为吹脱塔每秒流失水量，D为塔节点之间乘积最大距离三、吹脱塔的选购要求：1. 首先要了解吹脱塔使用的介质的性能和分子式，以便选择合适的材料。

2. 根据实际使用的流量来确定塔的容量，以达到设备的理想功效。

3. 吹脱塔要采用工艺管线安装，设计用于外部环境配备；4. 工作条件下，考虑流量，吹脱塔的压力也需要考虑，并且要比所选择材料的载压等级更高；5. 吹脱塔使用环境：空气环境温度（℃）和相对湿度（％）。

6. 吹脱塔的安装方式和安装高度也要考虑。

7. 吹脱塔的运行功效以及排放的污染物的水平一定要达到规定的标准，以确保安全性。

脱丙烷塔校核计算单位尚川计算条件塔型板式设计压力MPa 1.1容器分段数（不包括裙座）1压力试验类型液压压力试验计入液柱高度H mm 28950试验压力（立试）MPa 1.375试验压力（卧试）MPa 1.659封头上封头下封头材料名称Q345R Q345R名义厚度mm 1212腐蚀裕量mm 22焊接接头系数11封头形状椭圆形椭圆形圆筒1 2 3 4 5 设计温度℃90圆筒长度mm 28500圆筒名义厚度mm 12圆筒内径mm 1600材料名称（即钢号）Q345R腐蚀裕量mm 2纵向焊接接头系数1环向焊接接头系数1圆筒外压计算长度mm 06 7 8 9 10 设计温度℃圆筒长度mm圆筒名义厚度mm圆筒内径mm材料名称（即钢号）腐蚀裕量mm纵向焊接接头系数环向焊接接头系数圆筒外压计算长度mm内件及偏心载荷介质密度kg/m3467塔釜液面离焊接接头的高度mm 1290塔板分段数 1 2 3 4 5 塔板型式浮阀塔板层数38每层塔板上积液厚度mm 60最高一层塔板高度mm 26500最低一层塔板高度mm 3200填料分段数 1 2 3 4 5 填料顶部高度mm填料底部高度mm填料密度kg/m3集中载荷数 1 2 3 4 5 集中载荷kg 5050集中载荷高度mm 245009500集中载荷中心至容器中mm 900900心线距离塔器附件及基础塔器附件质量计算系数 1.2基本风压N/m2450基础高度mm 220塔器保温层厚度mm 50保温层密度kg/m3145裙座防火层厚度mm 30防火层密度kg/m32000管线保温层厚度mm 30最大管线外径mm 950笼式扶梯与最大管线的相对90位置场地土类型I场地土粗糙度类别B地震烈度低于7度地震远近参数近震塔器上平台总个数2平台宽度mm 1200塔器上最高平台高度mm 2700塔器上最低平台高度mm 4500裙座裙座结构形式圆筒形裙座底部截面内径mm 1600裙座与壳体连接形式对接裙座高度mm 3200裙座材料名称Q245R裙座设计温度℃90裙座腐蚀裕量mm 2裙座名义厚度mm 20裙座材料许用应力MPa 113裙座上同一高度处较大孔个数2裙座较大孔中心高度mm 800裙座上较大孔引出管内径（或宽度）mm 800裙座上较大孔引出管厚度mm 16裙座上较大孔引出管长度mm 1000地脚螺栓及地脚螺栓座地脚螺栓材料名称地脚螺栓材料许用应力MPa 0注：以下设计参数均参照JB4710-92 表5-6 并计算确定地脚螺栓个数8地脚螺栓公称直径mm 64全部筋板块数16相邻筋板最大外侧间距mm 567.99筋板内侧间距mm 120筋板厚度mm 22筋板宽度mm 170盖板类型分块盖板上地脚螺栓孔直径mm 85盖板厚度mm 53盖板宽度mm 220垫板有垫板上地脚螺栓孔直径mm 67垫板厚度mm 24垫板宽度mm 120基础环板外径mm 1864基础环板内径mm 1344基础环板名义厚度mm 28计算结果容器壳体强度计算元件名称压力设计名义厚度(mm) 直立容器校核取用厚度(mm)许用内压(MPa) 许用外压(MPa)下封头 12 12 2.292 第1 段圆筒12 12 2.112 第1 段变径段第2 段圆筒第2 段变径段第3 段圆筒第3 段变径段第4 段圆筒第4 段变径段第5 段圆筒第5 段变径段第6 段圆筒第6 段变径段第7 段圆筒第7 段变径段第8 段圆筒第8 段变径段第9 段圆筒第9 段变径段第10 段圆筒上封头12 12 2.292裙座名义厚度(mm) 取用厚度(mm)20 20风载及地震载荷0－0 A－A 1－1(筒体) 1－1(下封头) 2－2 3－3 4－4 操作质量32816.8 31778.9 28332.6 28332.6最小质量24595.7 23557.8 20111.5 20111.5液压试验时质量87717.7 86679.8 24691 24691风弯矩9.9414e+08 9.53869e+088.35764e+088.35764e+08地震弯矩0 0 0 0偏心弯矩882900 882900 882900 882900最大弯矩9.95023e+08 9.54752e+088.36647e+088.36647e+08垂直地震力0 0 0 0应力计算σ110.00 0.00 44.00 44.00 σ12 3.76 2.58 5.53 5.53 σ1329.36 16.15 41.61 41.61 σ22 3.76 2.58 3.93 3.93 σ310.00 0.00 55.00 55.00 σ3210.05 7.03 4.82 4.82 σ338.83 4.85 12.51 12.51 [σ]t113.00 113.00 170.00 170.00 B 137.67 137.67 133.94 133.94组合应力校核σA180.08 80.08 许用值204.00 204.00 σA233.12 18.72 45.54 45.54 许用值135.60 135.60 160.72 160.72 σA362.70 62.70 许用值351.00 372.60σA418.88 11.88 17.33 17.33许用值165.20 165.20 160.72 160.72σ133.19 133.19许用值292.50 310.50校核结果合格合格合格合格注1: σi j中i 和j 的意义如下i=1 操作工况j=1 设计压力或试验压力下引起的轴向应力（拉）i=2 检修工况j=2 重力及垂直地震力引起的轴向应力（压）i=3 液压试验工况j=3 弯矩引起的轴向应力（拉或压）[σ]t设计温度下材料许用应力 B 设计温度下轴向稳定的应力许用值注2:σA1: 操作工况下轴向最大组合拉应力σA2: 操作工况下轴向最大组合压应力σA3: 液压试验时轴向最大组合拉应力σA4: 液压试验时轴向最大组合压应力σ: 试验压力引起的周向应力注3: 单位如下质量: kg 力:N 弯矩: N mm 应力: MPa计算结果地脚螺栓及地脚螺栓座基础环板抗弯断面模数mm3 4.63974e+08 基础环板面积mm2 1.31017e+06基础环板计算力矩N∙mm 16702.9 基础环板需要厚度mm 26.76 基础环板厚度厚度校核结果合格混凝土地基上最大压应力MPa 2.39地脚螺栓受风载时最大拉应力MPa 1.96 地脚螺栓受地震载荷时最大拉应力MPa 0.29地脚螺栓需要的螺纹小径mm 55.7337 地脚螺栓实际的螺纹小径mm 57.505 地脚螺栓校核结果合格筋板压应力MPa 60.81 筋板许用应力MPa 90.83 筋板校核结果合格盖板最大应力MPa 143.07 盖板许用应力MPa 140盖板校核结果¦不合格裙座与壳体的焊接接头校核焊接接头截面上的塔器操作质量kg28332.6 焊接接头截面上的最大弯矩N∙mm8.36647e+08对接接头校核搭接接头校核对接接头横截面mm285592 搭接接头横截面mm2对接接头抗弯断面模数mm3 3.38944e+07 搭接接头抗剪断面模数mm3对接焊接接头在操作工况下最大拉应力MPa21.44 搭接焊接接头在操作工况下最大剪应力MPa对接焊接接头拉应力许可值MPa81.36 搭接焊接接头在操作工况下的剪应力许可值MPa对接接头拉应力校核结果合格搭接焊接接头在试验工况下最大剪应力MPa搭接焊接接头在试验工况下的剪应力许可值MPa搭接接头拉应力校核结果主要尺寸设计及总体参数计算结果裙座设计名义厚度mm 20壳体和裙座质量kg 16682 附件质量kg 3336.39 内件质量kg 5724.42 保温层质量kg 2064.46 平台及扶梯质量kg 1268 操作时物料质量kg 3641.59 直立容器的操作质量kg 32816.8 直立容器的最小质量kg 24595.7 直立容器的最大质量kg 87717.7 液压试验时液体质量kg 58542.4 吊装时空塔质量kg 21163.2直立容器自振周期s 0.85 空塔重心至基础mm 15007.8 环板底截面上风弯矩N∙mm 9.9414e+08 环板底截面距离环板底截面上地震弯矩N∙mm 0 环板底截面上垂直地震力N 0操作时基础环板底截面的最大计算弯矩N∙mm 9.95023e+08风载对直立容器总的横推力N 50621.6 地震载荷对直立容器总的横推力N 0操作工况下容器顶部最大挠度mm 36.7446 容器许用外压MPa容器总容积mm3 5.85424e+10 直立容器总高mm 32162 第二振型自振周期s 0.14 第三振型自振周期s 0.05 注：内件质量指塔板质量,填料质量计入物料质量。

气分操作指标
一、脱丙烷塔工艺参数操作指标。

塔顶温度46℃左右
塔底灵敏板温度85℃左右
塔底温度95-100℃
塔顶压力 1.5-1.7MPa
二、脱丙烷他产品质量控制指标。

1、塔顶丙烷纯度：干气+丙烷+丙烯=99%以上为合格
2、塔底部碳四中丙烷含量小于1% 合格
三、紧急工况的处理。

（停水、电、汽、分）
1、遇到紧急工况需要停工时，首先停进料，关闭进料调节阀FV14501，关闭塔底温度控制调节阀TV14504，防止塔内超压，若压力高接近安全阀起跳压力时，打开回流罐顶部压力控制调节阀PV14503泄压。

2、紧急联系芳构化装置，切断C4采出至芳构化流量调节阀FV14502，并关闭后手阀。

防止串料。

3、关闭丙烷采出液位控制调节阀LV14502，并关闭后手阀。

防止不合格产品采出对产品罐的质量影响。

4、关闭回流泵P1405出口阀，及塔顶空冷。

（若停电时，只需关闭泵的出口阀。

）。

《过程控制工程》课程设计任务书一、设计题目：脱丙烷塔控制系统设计二、设计目的：1、掌握控制系统的基本构成、原理及设计的方法和步骤。

2、掌握控制方案的设计、仪表选型的方法及管道流程图、仪表接线图、仪表安装等图的绘制方法。

3、掌握节流装置和调节阀的计算。

4、了解信号报警及联锁系统的设计和顺序控制系统的设计。

5、了解过程控制设计的设计文件构成及编制。

6、通过理论联系实际，掌握必须的工程知识，加强对学生实践动手能力和协作完成工程设计任务能力的培养。

三、设计所需数据：1、主要工艺流程和环境特征概况脱丙烷塔的主要任务是切割C3和C4混合馏分，塔顶轻关键组分是丙烷，塔釜重关键是丁二烯。

主要工艺流程如附图1所示：第一脱乙烷塔塔釜来的釜液和第二蒸出塔的釜液混合后进入脱丙烷塔，进料为气液混合状态，液化率为0.28。

进料温度为32℃，塔顶温度为8.9℃，塔釜温度为72℃。

塔内操作压力基本恒定在0.75MPa(绝压)。

采用的回流比约为1.13。

冷凝器由0℃丙烯蒸发制冷,再沸器加热用的0.15 MPa(绝压)减压蒸汽由来自裂解炉的0.6 MPa(绝压)低压蒸汽与冷凝水混合制得的。

和其他精馏塔一样，脱丙烷塔也是一个高阶对象，具有对象通道多、内在机理复杂、变量间相互关联、动态响应慢、控制要求高等特点。

假设该脱丙烷塔控制的主要目标是塔釜关键组分，可以再沸器的减压蒸汽流量为操纵变量构成控制系统，且此时再沸器的减压蒸汽流量是经常出现的扰动。

同时要保持塔进料稳定，以及塔釜液位与塔底A馏出物料均匀缓慢变化。

试设计自动控制，满足质量指标、物料指标、能量平衡及约束条件等要求。

脱丙烷塔所处的环境为甲级防爆区域,工艺介质为多种烃类混合物,沸点低、易挥发、易燃、易爆，生产装置处于露天，低压、低温。

主导风向由西向东。

2、仪表选型说明所选仪表应具有本质安全防爆性能等特点，电动Ⅲ型仪表在安全性、可靠性等方面已能满足要求。

电动仪表信号传送快且距离远，易与计算机配合使用，除控制阀外，可选用电动Ⅲ型仪表或采用数字式控制仪表。

脱丙烷塔原理
《脱丙烷塔原理》
嘿，今天咱就来讲讲脱丙烷塔原理哈。

话说有一次我去参观一个化工厂，那里面就有脱丙烷塔呢。

我就像个好奇宝宝似的，凑过去仔细瞧。

这个脱丙烷塔啊，就好像是一个超级大的魔法罐子。

它的作用呢，就是把混合物中的丙烷给分离出来。

你可以想象一下，就像我们在一堆糖果里，要把红色的糖果挑出来一样。

它是怎么做到的呢？其实就是利用不同物质沸点的差异啦。

就像我们煮东西，有的东西容易煮熟，有的就需要煮久一点。

在脱丙烷塔里，温度慢慢升高，那些沸点低的物质就先变成气体跑出来啦，而丙烷呢，就乖乖地留在里面，等着被分离出来。

我在那看着，工人们在旁边忙碌着，各种管子啊、仪表啊，感觉好神奇。

我都有点看入迷了，一直在那琢磨这脱丙烷塔的原理，想着这玩意儿可真厉害呀，能把那么复杂的混合物给分得清清楚楚。

总之呢，脱丙烷塔就是通过这种巧妙的方式，把我们需要的丙烷给分离出来啦，是不是很有意思呀！这就是我对脱丙烷塔原理的一点小观察和体验啦，哈哈。

设计一脱丙烷塔。

已知进料量h kmol /100，原料压力MPa 0.1，温度50℃,组成如下表。

塔操作压力0.817()MPa A ，塔顶设全凝器，塔底设再沸器。

分离要求：塔顶异丁烷含量为0.06，塔底丙烷含量为0.06解：（一）、用简捷法得到如下基本参数（二）LM 法1、初步确定理论级数1）设8=S 、2=n 、6=m （包括塔釜、进料板）、74.1=R 3252.75=D 6748.24=W 逐板计算，结果列表：2）设7=S 、2=n 、5=m （包括塔釜、进料板）、74.1=R 3252.75=D 6748.24=W 逐板计算3）比较进料板液摩尔分数已经接近，可进入第一次循环。

2、第一次循环 1）塔顶塔底量调整1585.035.165.684977.05652.0=+-==∆A d 4004.062.538.42206.01871.0-=+-=∆B d4509.07076.122107.02924.22716.02716.02107.0-=+-=∆C d 0073.09972.4037.00028.00102.00102.00370.0=+-=∆D d归零化，使得∑=∆0d ，i iw d∆-=∆-2）根据调整后的数据进行塔的逐板计算，结果列表，各板的汽液流率和摩尔分数列表 3）温度分布 4）计算各板气液流率 5）计算换热器热负荷 6）计算各板汽体液体流率 7）核算各板气液组成（1）各板汽液流率和温度确定相对挥发度 （2）逐板计算3、采用同样的方法，经过4次循环，结果如下：基本达到要求。

故理论板数为7.。

脱丙烷塔原理今天来聊聊脱丙烷塔的原理。

你知道吗，在化工生产的这个大工厂里啊，脱丙烷塔可是个很重要的角色呢。

就好比咱们家庭里把混合的杂粮分类一样，脱丙烷塔干的也是这么个分类的活儿，不过它分的是化工原料里的丙烷和其他物质。

这就像从一堆各种颜色、各种大小的豆子里把某一种特别的豆子拣出来，只不过脱丙烷塔使用的不是手，而是靠不同物质的沸点差异这个特性。

咱们先来说说这个沸点。

沸点啊，你可以理解为液体开始沸腾变成气体的温度。

比如说咱们烧水，到100℃水就开了，这个100℃就是水的沸点。

不同的物质沸点不一样，丙烷也有自己的沸点。

老实说，我一开始也不明白脱丙烷塔怎么能把丙烷和其他的物质分离开来。

这就要说到脱丙烷塔的工作过程了。

它就像一个超级精准的分拣机器，原料从塔底进入，然后加热，这时不同沸点的物质就开始了“变身秀”。

沸点低的丙烷就最先按耐不住啦，就像性子最急的人最先冲出去一样，它很容易就变成气体跑到塔的上部。

而那些沸点高的物质呢，可能还得在塔底附近多待会儿，它们变化相对比较慢，还保持着液态的状态。

说到这里，你可能会问，那怎么最终把丙烷和其他物质彻底分开呢？塔顶有专门的装置，可以把变成气体的丙烷收集起来，这样就把丙烷和其他的物质分离开来了。

打个比方，这就像我们把那些先冲出门去（指变成气体跑到塔顶）的“急性子”（丙烷）拦住，然后让其他没出门（仍然在塔底附近的物质）和先出去的这些分开，避免混在一起。

从更原理的层面来说呢，这是基于相平衡理论的。

物质在不同温度压力下会在气液两相之间进行分配，脱丙烷塔巧妙地利用了这个特性，根据其需要的分离目标调整温度和压力等条件。

实际应用的案例也很多啊。

比如说在石油炼制过程中，原油经过一系列的加工得到了很多混合的烃类物质，这时候就需要脱丙烷塔把丙烷这种比较轻的烃类分离开来。

如果不分开的话，这些不同的物质混在一起下一步的加工就不好进行啦，就像你做饭的时候一堆不同的调料混在一起没分开储存，取用的时候就很麻烦。

成绩：《过程控制工程》课程设计报告题目：脱丙烷塔控制系统设计学院：计算机与电子信息学院班级：自动化姓名：学号：指导教师：起止日期：2012年12月31日～2013年01月4日目录一、设计任务书 (2)二、设计说明书 (5)1、摘要2、基本控制方案的设计与分析3、节流装置的计算4、蒸汽流量控制阀口径的计算三、参考文献 (11)四、附图 (15)一、设计题目：《脱丙烷塔控制系统设计》二、设计目的：1、掌握控制系统的基本构成、原理及设计的方法和步骤。

2、掌握控制方案的设计、仪表选型的方法及管道流程图、仪表接线图、仪表安装等图的绘制方法。

3、掌握节流装置和调节阀的计算。

4、了解信号报警及联锁系统的设计和顺序控制系统的设计。

5、通过理论联系实际，掌握必须的工程知识，加强对学生实践动手能力和独立完成工程设计任务能力的培养。

三、设计所需数据：1、主要工艺流程和环境特征概况脱丙烷塔的主要任务是切割C3和C4混合馏分，塔顶轻关键组分是丙烷，塔釜重关键是组分丁二烯。

主要工艺流程如图1所示：第一脱乙烷塔塔釜来的釜液和第二蒸出塔的釜液混合后进入脱丙烷塔，进料为气液混合状态，液化率为0.28。

进料温度为32℃，塔顶温度为8.9℃，塔釜温度为72℃。

塔内操作压力为0.75MPa(绝压)。

采用的回流比约为1.13。

冷凝器由0℃丙烯蒸发制冷,再沸器加热用的0.15 MPa(绝压)减压蒸汽由来自裂解炉的0.6 MPa(绝压)低压蒸汽与冷凝水混合制得的。

和其他精馏塔一样，脱丙烷塔也是一个高阶对象，具有对象通道多、内在机理复杂、变量间相互关联、动态响应慢、控制要求高等特点。

脱丙烷塔的自动控制应满足质量指标、物料指标、能量平衡及约束条件等要求。

脱丙烷塔所处的环境为甲级防爆区域,工艺介质为多种烃类混合物,沸点低、易挥发、易燃、易爆，生产装置处于露天，低压、低温。

主导风向由西向东。

2、仪表选型说明所选仪表应具有本质安全防爆性能等特点，电动Ⅲ型仪表在安全性、可靠性等方面已能满足要求。

on Base 1Const n As Integer = 3 '输出个数Const m As Integer = 3 '输入个数Const d As Integer = 100 ' 系统最大时滞Dim ts As Single '采样时间Dim tr(n) As Single '参考轨迹柔化系数Dim H(n) As Single '预测时域长度'标志位定义Dim APC_COM As IntegerDim APC_ON As IntegerDim APC_2_ON As IntegerDim T_APC_STA As IntegerDim MV_1_STA As IntegerDim MV_2_STA As IntegerDim MV_3_STA As IntegerDim DV_1_STA As IntegerDim DV_2_STA As Integer'模型参数Dim K11, K12, K13, K21, K22, K23, K31, K32, K33 As SingleDim T11, T12, T13, T21, T22, T23, T31, T32, T33 As SingleDim TAO11, TAO12, TAO13, TAO21, TAO22, TAO23, TAO31, TAO32, TAO33 As SingleDim KM11, KM12, KM13, KM21, KM22, KM23, KM31, KM32, KM33 As SingleDim TM11, TM12, TM13, TM21, TM22, TM23, TM31, TM32, TM33 As SingleDim TAOM11, TAOM12, TAOM13, TAOM21, TAOM22, TAOM23, TAOM31, TAOM32, TAOM33 As SingleDim K(n, m) As SingleDim t(n, m) As SingleDim tao(n, m) As SingleDim km(n, m) As SingleDim tm(n, m) As SingleDim taom(n, m) As Single'变量定义Dim I As IntegerDim J As IntegerDim g As IntegerDim AA(n, m) As SingleDim AAA() As SingleDim BB(n, 1) As SingleDim b(n) As SingleDim a(n, m) As SingleDim am(n, m) As SingleDim L(n, m) As SingleDim Lm(n, m) As SingleDim c(n) As SingleDim piduu1 As SingleDim piduu2 As SingleDim piduu3 As SingleDim ym(n, m) As Single '模型输出Dim ym1(n, m, d) As Single ' 模型上次输出Dim y0 As SingleDim yp0 As SingleDim yp(n) As Single ' 过程输出Dim y(n, m) As SingleDim ypav(n) As SingleDim ymm0 As SingleDim ym0(n) As SingleDim Y1(n, m, 1) As SingleDim u1(m, d) As SingleDim uu(m) As Single'力控数据通讯所需参数Dim read_datas(22) As Double '读数据数组Dim readtag As StringDim read_data As LongDim write_datas(6) As Double '写数据数组Dim writetag As StringDim write_data As LongPrivate Sub Form_Load()'过程参数K11 = 4.05: T11 = 50: TAO11 = 2K12 = 1.77: T12 = 60: TAO12 = 3K13 = 5.88: T13 = 50: TAO13 = 2K21 = 5.39: T21 = 50: TAO21 = 3K22 = 5.72: T22 = 60: TAO22 = 2K23 = 6.9: T23 = 40: TAO23 = 3K31 = 4.38: T31 = 33: TAO31 = 2K32 = 4.42: T32 = 44: TAO32 = 3K33 = 7.52: T33 = 19: TAO33 = 3K(1, 1) = K11: K(1, 2) = K12: K(1, 3) = K13: K(2, 1) = K21: K(2, 2) = K22: K(2, 3) = K23: K(3, 1) = K31: K(3, 2) = K32: K(3, 3) = K33t(1, 1) = T11: t(1, 2) = T12: t(1, 3) = T13: t(2, 1) = T21: t(2, 2) = T22: t(2, 3) = T23: t(3, 1) = T31: t(3, 2) = T32: t(3, 3) = T33tao(1, 1) = TAO11: tao(1, 2) = TAO12: tao(1, 3) = TAO13: tao(2, 1) = TAO21: tao(2, 2) = TAO22: tao(2, 3) = TAO23: tao(3, 1) = TAO31: tao(3, 2) = TAO32: tao(3, 3) = TAO33Text29.Text = Val(K(1, 1)): Text30.Text = Val(K(1, 2)): Text31.Text = Val(K(1, 3)): Text32.Text = Val(K(2, 1)): Text33.Text = Val(K(2, 2)): Text34.Text = Val(K(2, 3)): Text35.Text = Val(K(3, 1)): Text36.Text = Val(K(3, 2)): Text37.Text = Val(K(3, 3))Text38.Text = Val(t(1, 1)): Text39.Text = Val(t(1, 2)): Text40.Text = Val(t(1, 3)): Text41.Text = Val(t(2, 1)): Text42.Text = Val(t(2, 2)): Text43.Text = Val(t(2, 3)): Text44.Text = Val(t(3, 1)): Text45.Text = Val(t(3, 2)): Text46.Text = Val(t(3, 3))Text47.Text = Val(tao(1, 1)): Text48.Text = Val(tao(1, 2)): Text49.Text = Val(tao(1, 3)): Text50.Text = Val(tao(2, 1)): Text51.Text = Val(tao(2, 2)): Text52.Text = Val(tao(2, 3)): Text53.Text = Val(tao(3, 1)): Text54.Text = Val(tao(3, 2)): Text55.Text = V al(tao(3, 3))'模型参数KM11 = 4.05: TM11 = 50: TAOM11 = 2KM12 = 1.77: TM12 = 60: TAOM12 = 3KM13 = 5.88: TM13 = 50: TAOM13 = 2KM21 = 5.39: TM21 = 50: TAOM21 = 3KM22 = 5.72: TM22 = 60: TAOM22 = 2KM23 = 6.9: TM23 = 40: TAOM23 = 3KM31 = 4.38: TM31 = 33: TAOM31 = 2KM32 = 4.42: TM32 = 44: TAOM32 = 3KM33 = 7.52: TM33 = 19: TAOM33 = 3km(1, 1) = KM11: km(1, 2) = KM12: km(1, 3) = KM13: km(2, 1) = KM21: km(2, 2) = KM22: km(2, 3) = KM23: km(3, 1) = KM31: km(3, 2) = KM32: km(3, 3) = KM33tm(1, 1) = TM11: tm(1, 2) = TM12: tm(1, 3) = TM13: tm(2, 1) = TM21: tm(2, 2) = TM22: tm(2, 3) = TM23: tm(3, 1) = TM31: tm(3, 2) = TM32: tm(3, 3) = TM33taom(1, 1) = TAOM11: taom(1, 2) = TAOM12: taom(1, 3) = TAOM13: taom(2, 1) = TAOM21: taom(2, 2) = TAOM22: taom(2, 3) = TAOM23: taom(3, 1) = TAOM31: taom(3, 2) = TAOM32: taom(3, 3) = TA0M33'初始化变量For I = 1 To nFor J = 1 To mym(I, J) = 0y(I, J) = 0uu(J) = 0For g = 1 To dym1(I, J, g) = 0NextNextyp(I) = 0ypav(I) = 0ym0(I) = 0Nextpiduu1 = 0piduu2 = 0piduu3 = 0End SubPrivate Sub Timer1_Timer()'从力控读数据readtag = "T_APC_TS.pv,T_APC_H1.pv,T_APC_H2.pv,T_APC_H3.pv,T_APC_TR1.pv,T_APC_TR2.pv,T _APC_TR3.pv,TI_707_1_SP.pv,TIC_706_SP.pv,LIC_703_SP.pv,APC_COM.pv,APC_ON.pv,AP C_2_ON.pv,T_APC_STA.pv,FIC_730_MV.pv,TIC_706_MV.pv,FIC_704_MV.pv,FIC_702_DV.pv ,TIC_701_DV.pv,PID_U1.pv,PID_U2.pv,PID_U3.pv"read_data = Dbcom1.GetRealData(22, readtag, read_datas(1))Text1.Text = read_datas(1)Text2.Text = read_datas(2)Text3.Text = read_datas(3)Text4.Text = read_datas(4)Text5.Text = read_datas(5)Text6.Text = read_datas(6)Text7.Text = read_datas(7)Text14.Text = read_datas(8)Text15.Text = read_datas(9)Text16.Text = read_datas(10)Text20.Text = read_datas(11)Text21.Text = read_datas(12)Text22.Text = read_datas(13)Text23.Text = read_datas(14)Text24.Text = read_datas(15)Text25.Text = read_datas(16)Text26.Text = read_datas(17)Text27.Text = read_datas(18)Text28.Text = read_datas(19)Text17.Text = read_datas(20)Text18.Text = read_datas(21)Text19.Text = read_datas(22)'VB中变量赋值ts = Val(Text1.Text)H(1) = Val(Text2.Text)H(2) = Val(Text3.Text)H(3) = Val(Text4.Text)tr(1) = Val(Text5.Text)tr(2) = Val(Text6.Text)tr(3) = Val(Text7.Text)c(1) = Val(Text14.Text)c(2) = Val(Text15.Text)c(3) = Val(Text16.Text)APC_COM = Val(Text20.Text) APC_ON = Val(Text21.Text)APC_2_ON = Val(Text22.Text)T_APC_STA = Val(Text23.Text)MV_1_STA = Val(Text24.Text) MV_2_STA = Val(Text25.Text) MV_3_STA = Val(Text26.Text)DV_1_STA = Val(Text27.Text)DV_2_STA = Val(Text28.Text)piduu1 = Val(Text17.Text)piduu2 = Val(Text18.Text)piduu3 = Val(Text19.Text)'主程序部分For I = 1 To nb(I) = Exp(-ts / tr(I))For J = 1 To mam(I, J) = Exp(-ts / tm(I, J))Lm(I, J) = Round(taom(I, J) / ts) NextNextIf APC_COM = 1 ThenIf APC_ON = 1 ThenIf APC_2_ON = 1 Thenuu(1) = piduu1uu(2) = piduu2uu(3) = piduu3Text11.Text = uu(1)Text12.Text = uu(2)Text13.Text = uu(3)'模型计算子程序部分开始'回路一单独投用If MV_1_STA = 1 And MV_2_STA = 0 And MV_3_STA = 0 Then'********************************************ReDim AAA(1, 1) As Single'********************************************For I = 1 To 1For J = 1 To 1ym(I, J) = am(I, J) * ym1(I, J, 1) + km(I, J) * (1 - am(I, J)) * uu(J) 'ym1为y上一次的值For g = Lm(I, J) To 2 Step -1ym1(I, J, g) = ym1(I, J, g - 1)Nextym1(I, J, 1) = ym(I, J)NextNextFor I = 1 To 1y0 = 0yp0 = 0For J = 1 To 1yp(I) = yp0 + y(I, J)yp0 = yp(I)ypav(I) = y0 + y(I, J) + ym(I, J) - ym1(I, J, 1)y0 = ypav(I)NextNextFor I = 1 To 1ymm0 = 0For J = 1 To 1AA(I, J) = km(I, J) * (1 - am(I, J) ^ H(I))ym0(I) = ymm0 + (1 - am(I, J) ^ H(I)) * ym(I, J)ymm0 = ym0(I)NextBB(I, 1) = c(I) - b(I) ^ H(I) * (c(I) - yp(I)) - ypav(I) + ym0(I) Next'构造新的求逆矩阵AAA(1, 1) = AA(1, 1)Call Gauss_Inv(AAA(), 1)uu(1) = AAA(1, 1) * BB(1, 1)'If T_APC_STA = 0 Then '投用前赋值' uu(1) = piduu1'End IfText8.Text = yp(1)' Text11.Text = uu(1)End If'回路二单独投用If MV_1_STA = 0 And MV_2_STA = 1 And MV_3_STA = 0 Then'********************************************ReDim AAA(1, 1) As Single'********************************************For I = 2 To 2For J = 2 To 2ym(I, J) = am(I, J) * ym1(I, J, 1) + km(I, J) * (1 - am(I, J)) * uu(J) 'ym1为y上一次的值For g = Lm(I, J) To 2 Step -1ym1(I, J, g) = ym1(I, J, g - 1)Nextym1(I, J, 1) = ym(I, J)NextNextFor I = 2 To 2y0 = 0yp0 = 0For J = 2 To 2yp(I) = yp0 + y(I, J)yp0 = yp(I)ypav(I) = y0 + y(I, J) + ym(I, J) - ym1(I, J, 1)y0 = ypav(I)NextNextFor I = 2 To 2ymm0 = 0For J = 2 To 2AA(I, J) = km(I, J) * (1 - am(I, J) ^ H(I))ym0(I) = ymm0 + (1 - am(I, J) ^ H(I)) * ym(I, J)ymm0 = ym0(I)NextBB(I, 1) = c(I) - b(I) ^ H(I) * (c(I) - yp(I)) - ypav(I) + ym0(I)Next'构造新的求逆矩阵AAA(1, 1) = AA(2, 2)Call Gauss_Inv(AAA(), 1)uu(1) = AAA(1, 1) * BB(1, 1)' If T_APC_STA = 0 Then '投用前赋值' uu(2) = piduu2' End IfText9.Text = yp(2)' Text12.Text = uu(2)End If'回路三单独投用If MV_1_STA = 0 And MV_2_STA = 0 And MV_3_STA = 1 Then'********************************************ReDim AAA(1, 1) As Single'********************************************For I = 3 To 3For J = 3 To 3ym(I, J) = am(I, J) * ym1(I, J, 1) + km(I, J) * (1 - am(I, J)) * uu(J) 'ym1为y上一次的值For g = Lm(I, J) To 2 Step -1ym1(I, J, g) = ym1(I, J, g - 1)Nextym1(I, J, 1) = ym(I, J)NextNextFor I = 3 To 3y0 = 0yp0 = 0For J = 3 To 3yp(I) = yp0 + y(I, J)yp0 = yp(I)ypav(I) = y0 + y(I, J) + ym(I, J) - ym1(I, J, 1)y0 = ypav(I)NextNextFor I = 3 To 3ymm0 = 0For J = 3 To 3AA(I, J) = km(I, J) * (1 - am(I, J) ^ H(I))ym0(I) = ymm0 + (1 - am(I, J) ^ H(I)) * ym(I, J)ymm0 = ym0(I)NextBB(I, 1) = c(I) - b(I) ^ H(I) * (c(I) - yp(I)) - ypav(I) + ym0(I)Next'构造新的求逆矩阵AAA(1, 1) = AA(3, 3)Call Gauss_Inv(AAA(), 1)uu(3) = AAA(1, 1) * BB(3, 1)'If T_APC_STA = 0 Then '投用前赋值' uu(3) = piduu3' End IfText10.Text = yp(3)' Text13.Text = uu(3)End If'回路一、二投用If MV_1_STA = 1 And MV_2_STA = 1 And MV_3_STA = 0 Then'********************************************ReDim AAA(2, 2) As Single'********************************************For I = 1 To 2For J = 1 To 2ym(I, J) = am(I, J) * ym1(I, J, 1) + km(I, J) * (1 - am(I, J)) * uu(J) 'ym1为y上一次的值For g = Lm(I, J) To 2 Step -1ym1(I, J, g) = ym1(I, J, g - 1)Nextym1(I, J, 1) = ym(I, J)NextNextFor I = 1 To 2y0 = 0yp0 = 0For J = 1 To 2yp(I) = yp0 + y(I, J)yp0 = yp(I)ypav(I) = y0 + y(I, J) + ym(I, J) - ym1(I, J, 1)y0 = ypav(I)NextNextFor I = 1 To 2ymm0 = 0For J = 1 To 2AA(I, J) = km(I, J) * (1 - am(I, J) ^ H(I))ym0(I) = ymm0 + (1 - am(I, J) ^ H(I)) * ym(I, J)ymm0 = ym0(I)NextBB(I, 1) = c(I) - b(I) ^ H(I) * (c(I) - yp(I)) - ypav(I) + ym0(I)Next'构造新的求逆矩阵AAA(1, 1) = AA(1, 1)AAA(1, 2) = AA(1, 2)AAA(2, 1) = AA(2, 1)AAA(2, 2) = AA(2, 2)Call Gauss_Inv(AAA(), 2) '调用矩阵求逆函数'uu(1) = AAA(1, 1) * BB(1, 1) + AAA(1, 2) * BB(2, 1) '过程通道预测控制量uu(1)为回流量，uu(2)为蒸汽量uu(2) = AAA(2, 1) * BB(1, 1) + AAA(2, 2) * BB(2, 1)' If T_APC_STA = 0 Then '投用前赋值' uu(1) = piduu1' uu(2) = piduu2' End IfText8.Text = yp(1)Text9.Text = yp(2)' Text11.Text = uu(1)' Text12.Text = uu(2)End If'回路一、三投用If MV_1_STA = 1 And MV_2_STA = 0 And MV_3_STA = 1 Then'********************************************ReDim AAA(2, 2) As Single'********************************************ym(1, 1) = am(1, 1) * ym1(1, 1, 1) + km(1, 1) * (1 - am(1, 1)) * uu(1) ym(1, 3) = am(1, 3) * ym1(1, 3, 1) + km(1, 3) * (1 - am(1, 3)) * uu(3) ym(3, 1) = am(3, 1) * ym1(3, 1, 1) + km(3, 1) * (1 - am(3, 1)) * uu(1) ym(3, 3) = am(3, 3) * ym1(3, 3, 1) + km(3, 3) * (1 - am(3, 3)) * uu(3)For g = Lm(1, 1) To 2 Step -1ym1(1, 1, g) = ym1(1, 1, g - 1)Nextym1(1, 1, 1) = ym(1, 1)For g = Lm(1, 3) To 2 Step -1ym1(1, 3, g) = ym1(1, 3, g - 1)Nextym1(1, 3, 1) = ym(1, 3)For g = Lm(3, 1) To 2 Step -1ym1(3, 1, g) = ym1(3, 1, g - 1)Nextym1(3, 1, 1) = ym(3, 1)For g = Lm(3, 3) To 2 Step -1ym1(3, 3, g) = ym1(3, 3, g - 1)Nextym1(3, 3, 1) = ym(3, 3)For I = 1 To n Step 2y0 = 0yp0 = 0For J = 1 To m Step 2yp(I) = yp0 + y(I, J)yp0 = yp(I)ypav(I) = y0 + y(I, J) + ym(I, J) - ym1(I, J, 1)y0 = ypav(I)NextNextFor I = 1 To n Step 2ymm0 = 0For J = 1 To m Step 2ym0(I) = ymm0 + (1 - am(I, J) ^ H(I)) * ym(I, J)ymm0 = ym0(I)NextNextAA(1, 1) = km(1, 1) * (1 - am(1, 1) ^ H(1))AA(1, 3) = km(1, 3) * (1 - am(1, 3) ^ H(1))AA(3, 1) = km(3, 1) * (1 - am(3, 1) ^ H(3))AA(3, 3) = km(3, 3) * (1 - am(3, 3) ^ H(3))BB(1, 1) = c(1) - b(1) ^ H(1) * (c(1) - yp(1)) - ypav(1) + ym0(1) BB(3, 1) = c(3) - b(3) ^ H(3) * (c(3) - yp(3)) - ypav(3) + ym0(3)'构造新的求逆矩阵AAA(1, 1) = AA(1, 1)AAA(1, 2) = AA(1, 3)AAA(2, 1) = AA(3, 1)AAA(2, 2) = AA(3, 3)Call Gauss_Inv(AAA(), 2)uu(1) = AAA(1, 1) * BB(1, 1) + AAA(1, 2) * BB(3, 1)uu(3) = AAA(2, 1) * BB(1, 1) + AAA(2, 2) * BB(3, 1)' If T_APC_STA = 0 Then '投用前赋值' uu(1) = piduu1' uu(3) = piduu3' End IfText8.Text = yp(1)Text10.Text = yp(3)' Text11.Text = uu(1)' Text13.Text = uu(3)End If'回路二、三投用If MV_1_STA = 0 And MV_2_STA = 1 And MV_3_STA = 1 Then'********************************************ReDim AAA(2, 2) As Single'********************************************ym(2, 2) = am(2, 2) * ym1(2, 2, 1) + km(2, 2) * (1 - am(2, 2)) * uu(2) ym(2, 3) = am(2, 3) * ym1(2, 3, 1) + km(2, 3) * (1 - am(2, 3)) * uu(3) ym(3, 2) = am(3, 2) * ym1(3, 2, 1) + km(3, 2) * (1 - am(3, 2)) * uu(2) ym(3, 3) = am(3, 3) * ym1(3, 3, 1) + km(3, 3) * (1 - am(3, 3)) * uu(3)For g = Lm(2, 2) To 2 Step -1ym1(2, 2, g) = ym1(2, 2, g - 1)Nextym1(2, 2, 1) = ym(2, 2)For g = Lm(2, 3) To 2 Step -1ym1(2, 3, g) = ym1(2, 3, g - 1)Nextym1(2, 3, 1) = ym(2, 3)For g = Lm(3, 2) To 2 Step -1ym1(3, 2, g) = ym1(3, 2, g - 1)Nextym1(3, 2, 1) = ym(3, 2)For g = Lm(3, 3) To 2 Step -1ym1(3, 3, g) = ym1(3, 3, g - 1)Nextym1(3, 3, 1) = ym(3, 3)For I = 2 To ny0 = 0yp0 = 0For J = 2 To myp(I) = yp0 + y(I, J)yp0 = yp(I)ypav(I) = y0 + y(I, J) + ym(I, J) - ym1(I, J, 1)y0 = ypav(I)NextNextFor I = 2 To nymm0 = 0For J = 2 To mym0(I) = ymm0 + (1 - am(I, J) ^ H(I)) * ym(I, J)ymm0 = ym0(I)NextNextAA(2, 2) = km(2, 2) * (1 - am(2, 2) ^ H(2))AA(2, 3) = km(2, 3) * (1 - am(2, 3) ^ H(2))AA(3, 2) = km(3, 2) * (1 - am(3, 2) ^ H(3))AA(3, 3) = km(3, 3) * (1 - am(3, 3) ^ H(3))BB(2, 1) = c(2) - b(2) ^ H(2) * (c(2) - yp(2)) - ypav(2) + ym0(2) BB(3, 1) = c(3) - b(3) ^ H(3) * (c(3) - yp(3)) - ypav(3) + ym0(3)'构造新的求逆矩阵AAA(1, 1) = AA(2, 2)AAA(1, 2) = AA(2, 3)AAA(2, 1) = AA(3, 2)AAA(2, 2) = AA(3, 3)Call Gauss_Inv(AAA(), 2)uu(2) = AAA(1, 1) * BB(2, 1) + AAA(1, 2) * BB(3, 1)uu(3) = AAA(2, 1) * BB(2, 1) + AAA(2, 2) * BB(3, 1) ' If T_APC_STA = 0 Then '投用前赋值' uu(2) = piduu2' uu(3) = piduu3'' End IfText9.Text = yp(2)Text10.Text = yp(3)' Text12.Text = uu(2)' Text13.Text = uu(3)End If'所有回路全部投用If MV_1_STA = 1 And MV_2_STA = 1 And MV_3_STA = 1 ThenFor I = 1 To nFor J = 1 To mym(I, J) = am(I, J) * ym1(I, J, 1) + km(I, J) * (1 - am(I, J)) * uu(J) 'ym1为y上一次的值For g = Lm(I, J) To 2 Step -1ym1(I, J, g) = ym1(I, J, g - 1)Nextym1(I, J, 1) = ym(I, J)NextNextFor I = 1 To ny0 = 0yp0 = 0For J = 1 To myp(I) = yp0 + y(I, J)yp0 = yp(I)ypav(I) = y0 + y(I, J) + ym(I, J) - ym1(I, J, 1)y0 = ypav(I)NextNextFor I = 1 To nymm0 = 0For J = 1 To mAA(I, J) = km(I, J) * (1 - am(I, J) ^ H(I))ym0(I) = ymm0 + (1 - am(I, J) ^ H(I)) * ym(I, J)ymm0 = ym0(I)NextBB(I, 1) = c(I) - b(I) ^ H(I) * (c(I) - yp(I)) - ypav(I) + ym0(I)NextCall Gauss_Inv(AA(), n) '调用矩阵求逆函数'uu(1) = AA(1, 1) * BB(1, 1) + AA(1, 2) * BB(2, 1) + AA(1, 3) * BB(3, 1) '过程通道预测控制量uu(1)为回流量，uu(2)为蒸汽量，uu(3)为塔底采出量uu(2) = AA(2, 1) * BB(1, 1) + AA(2, 2) * BB(2, 1) + AA(2, 3) * BB(3, 1)uu(3) = AA(3, 1) * BB(1, 1) + AA(3, 2) * BB(2, 1) + AA(3, 3) * BB(3, 1)'If T_APC_STA = 0 Then '投用前赋值' End IfText8.Text = yp(1)Text9.Text = yp(2)Text10.Text = yp(3)End If'子程序部分结束End If '子控制器切除标志位End If '紧急切除标志位End If '通讯状态标志位End SubPrivate Sub Timer2_Timer()'向力控送数据writetag = "T_APC_YOUT1.pv,T_APC_YOUT2.pv,T_APC_YOUT3.pv,T_APC_U1.pv,T_APC_U2.pv,T_A PC_U3.pv"write_datas(1) = Val(Text8.Text)write_datas(2) = Val(Text9.Text)write_datas(3) = Val(Text10.Text)write_datas(4) = Val(Text11.Text)write_datas(5) = Val(Text12.Text)write_datas(6) = Val(Text13.Text)write_data = Dbcom1.SetRealData(6, writetag, write_datas(1))For I = 1 To nFor J = 1 To ma(I, J) = Exp(-ts / t(I, J))L(I, J) = Round(tao(I, J) / ts)NextNextIf APC_COM = 1 ThenIf APC_ON = 1 ThenIf APC_2_ON = 1 Then' 过程计算子程序部分开始'回路一单独投用If MV_1_STA = 1 And MV_2_STA = 0 And MV_3_STA = 0 Theny(1, 1) = a(1, 1) * Y1(1, 1, 1) + K(1, 1) * (1 - a(1, 1)) * u1(1, L(1, 1))Y1(1, 1, 1) = y(1, 1)For g = L(1, 1) To 2 Step -1u1(1, g) = u1(1, g - 1)Nextu1(1, 1) = uu(1)End If'回路二单独投用If MV_1_STA = 0 And MV_2_STA = 1 And MV_3_STA = 0 Theny(2, 2) = a(2, 2) * Y1(2, 2, 1) + K(2, 2) * (1 - a(2, 2)) * u1(2, L(2, 2))Y1(2, 2, 1) = y(2, 2)For g = L(2, 2) To 2 Step -1u1(2, g) = u1(2, g - 1)Nextu1(2, 2) = uu(1)End If'回路三单独投用If MV_1_STA = 0 And MV_2_STA = 0 And MV_3_STA = 1 Theny(3, 3) = a(3, 3) * Y1(3, 3, 1) + K(3, 3) * (1 - a(3, 3)) * u1(3, L(3, 3))Y1(3, 3, 1) = y(3, 3)For g = L(3, 3) To 2 Step -1u1(3, g) = u1(3, g - 1)Nextu1(3, 3) = uu(3)End If'回路一、二投用If MV_1_STA = 1 And MV_2_STA = 1 And MV_3_STA = 0 ThenFor I = 1 To 2For J = 1 To 2y(I, J) = a(I, J) * Y1(I, J, 1) + K(I, J) * (1 - a(I, J)) * u1(J, L(I, J))Y1(I, J, 1) = y(I, J)For g = L(I, J) To 2 Step -1u1(J, g) = u1(J, g - 1)Nextu1(J, 1) = uu(J)NextNextEnd If'回路一、三投用If MV_1_STA = 1 And MV_2_STA = 0 And MV_3_STA = 1 ThenFor I = 1 To 3 Step 2For J = 1 To 3 Step 2y(I, J) = a(I, J) * Y1(I, J, 1) + K(I, J) * (1 - a(I, J)) * u1(J, L(I, J))Y1(I, J, 1) = y(I, J)For g = L(I, J) To 2 Step -1u1(J, g) = u1(J, g - 1)Nextu1(J, 1) = uu(J)NextNextEnd If'回路二、三投用If MV_1_STA = 0 And MV_2_STA = 1 And MV_3_STA = 1 ThenFor I = 2 To 3For J = 2 To 3y(I, J) = a(I, J) * Y1(I, J, 1) + K(I, J) * (1 - a(I, J)) * u1(J, L(I, J))Y1(I, J, 1) = y(I, J)For g = L(I, J) To 2 Step -1u1(J, g) = u1(J, g - 1)Nextu1(J, 1) = uu(J)NextNextEnd If'所有回路全部投用If MV_1_STA = 1 And MV_2_STA = 1 And MV_3_STA = 1 ThenFor I = 1 To nFor J = 1 To my(I, J) = a(I, J) * Y1(I, J, 1) + K(I, J) * (1 - a(I, J)) * u1(J, L(I, J))Y1(I, J, 1) = y(I, J)For g = L(I, J) To 2 Step -1u1(J, g) = u1(J, g - 1)Nextu1(J, 1) = uu(J)NextNextEnd If'子程序部分结束End If '子控制器切除标志位End If '紧急切除标志位End If '通讯状态标志位。

脱丙烷塔高度
脱丙烷塔是石油化工中常用的设备之一，用于将原料中的丙烷与丁烷进行分离。

以下是有关脱丙烷塔高度的详细解释：
脱丙烷塔的高度取决于多个因素，其中最重要的是分离要求和操作条件。

一般来说，脱丙烷塔的高度在30-50米之间，具体高度需要根据工艺要求和设备规格来确定。

脱丙烷塔的操作原理是利用物质的沸点不同，通过加热和冷却的方法使不同沸点的物质得到分离。

在脱丙烷塔中，原料进入塔内后，经过加热和蒸发，丙烷和丁烷以气态形式上升，而重质烃则下沉。

为了使丙烷和丁烷更好地分离，塔内通常设置有多层塔板或填料，以增加气体的接触面积和分离效果。

脱丙烷塔的高度取决于多个因素，其中最重要的是分离要求和操作条件。

在确定脱丙烷塔的高度时，需要考虑以下因素：
1.分离要求：根据工艺要求，需要将原料中的丙烷和丁烷进行不同程度的分离。

分离要求
越高，所需的塔高越大。

2.操作条件：操作条件如进料流量、温度、压力等也会影响塔高。

在确定塔高时，需要考
虑这些因素对分离效果的影响。

3.设备规格：不同厂家和型号的脱丙烷塔设备规格不同，塔高也会有所不同。

在选择脱丙
烷塔时，需要根据自己的工艺要求和设备规格进行选择。

总之，脱丙烷塔的高度是综合考虑多个因素的结果，需要根据具体情况进行确定。

本科毕业论文（设计）年产10万吨丙烯分离工段工艺设计姓名：指导教师：院系：化学化工学院专业：化学工程与工艺提交日期：2012年5月5日目录中文摘要 (1)外文摘要 (2)引言 (3)1.绪论 (3)1.1概述 (3)1.1.1简介 (3)1.1.2丙烯的性质 (3)1.1.3丙烯的用途 (3)1.2丙烯生产工艺选择及分离流程确定 (3)1.2.1生产工艺选择 (3)1.2.2分离流程确定 (4)1.3设计任务书 (5)2.工艺流程 (5)2.1工艺流程图 (5)2.2工艺流程简述 (6)3.物料衡算 (6)3.1设计依据 (6)3.2裂解气及各组分产量 (6)3.3各裂解产物的相对分子量 (7)3.4脱丙烷塔物料衡算 (7)3.5脱甲烷塔物料衡算 (10)3.6脱乙烷塔物料衡算 (12)3.7乙烯精馏塔物料衡算 (15)3.8丙烯精馏塔物料衡算 (16)4.热量衡算 (18)4.1乙烯精馏装置热量衡算 (18)4.2丙烯精馏装置热量衡算 (23)4.3脱甲烷精馏装置热量衡算 (26)4.4脱乙烷精馏装置热量衡算 (30)4.5脱丙烷精馏装置热量衡算 (33)5.设备选型 (37)5.1丙烯精馏塔 (37)5.1.1丙烯精馏塔操作压力及温度的确定 (37)5.1.2丙烯精馏塔密度、表面张力的计算 (39)5.1.3塔板数的确定 (42)5.1.4精馏塔主要尺寸计算 (44)5.1.5塔板流体力学验算 (50)5.1.6主要设备设计与选型 (53)5.1.7塔高的计算 (55)5.1.8浮阀塔设计一览表 (56)5.2换热器 (57)5.2.1试算和初选换热器规格 (57)5.2.2核算总传热系数 (58)6.生产安全及三废处理 (62)6.1生产安全 (62)6.2废气处理 (62)6.3废渣处理 (62)6.4废水处理 (62)结束语 (63)参考文献 (64)致谢 (65)附录 (66)年产10万吨丙烯分离工段工艺设计刘洋指导老师：崔秀云（黄山学院化学化工学院，黄山，安徽245041）摘要：本设计为年产10万吨丙烯分离工段工艺设计。

典型的塔板效率经验数据石油化工塔名称塔板效率 (%)塔名称塔板效率 (%)脱乙烷塔 60～65二甲苯分离塔 90～95高压脱乙烷塔 50～60苯/甲苯/二甲苯分离塔 75～85脱丙烷塔 65～75苯/异丙苯分离塔 50～55脱丁烷塔 75～85吸收塔 20～35脱异丁烷/脱戊烷塔 80～90 解吸塔(再沸器供热) 40～50乙烷/乙烯分离塔 85～90 解吸塔(蒸汽汽提) 20～30丙烷/丙烯分离塔 90～95 气体汽提塔 7～10丁烷/丁烯分离塔 85～95 干燥塔 15戊烷/戊烯分离塔 85～95化工及其它塔名称塔板效率 (%) 塔名称塔板效率 (%) 醋酸乙烯装置聚乙烯醇装置醋酸精馏塔 60～70 聚醋酸乙烯分离塔(聚合一塔) 50～60醋酸乙烯精馏塔 55～65 醋酸乙烯/甲醇分离塔(聚合二塔) 50～55乙醛汽提塔20～25 醋酸乙烯分离塔(聚合三、四塔) 60～70丙酮萃取塔15～20 甲醇回收塔(回收一塔) 55～65洗涤塔30～40 甲醇回收塔(回收二、三塔) 50～60 在ASPENPLUS模拟软件中，关于塔板效率的问题，模拟软件普遍采用默认值。

一般模拟软件在回归计算时候，在采用简捷计算和精确计算的模式时，塔板效率的预设就不同了，在精确计算模式，一般可以人工修改这个值。

塔板效率属于经验值，这和软件的内部默认还不一样。

例如当建立设计模式精算时，塔板效率是有一个范围的，可是实际中，往往可以突破这个范围。

例如，默认的塔板效率最高可能就是100％，而实际中，特定的塔系确有设计为100％以上的情况，主要是看实际的情况而定。

计算塔板效率可以有公式计算：有公式可算，一般可以使用D_B 法和O‘connel法塔板效率的计算Murphree效率（实验研究用的比较多，设计计算用的比较少）实际塔板上物料进出口浓度差与达到平衡时浓度差之比：全塔效率：在一定回流比和要求之下，所需理论板数与实际板数之比1、Drickaner-Bradford 法E0=0.17-0.616LOG（u）2、O‘connell方法（用的最多）E0= 49{ ua}-0.25u-塔顶底平均温度下进料的黏度a—轻重关键组分的平均相对挥发度3、其他方法AICHE方法（1958）Winckle方法（1963）吸收塔的效率最低，大约只有10-30%；解吸塔和吸收蒸出塔在30-45之间；一般的塔在60—90%之间。

脱丙烷精馏塔设计摘要化工生产中做处理的原料、中间产物有若干组分组成的混合物，在化工、炼油、医药、食品即环境保护等工业部门，精馏过程在能量计的驱动下，气液两相多次直接接触和分离，利用气液两相各相份挥发度不同时挥发组分由液相向气相转移，实现原料混合物中各组分同时进行传质传热过程。

塔设备是一种重要的单元操作设备。

它的应用面广、量大。

据统计，塔设备无论其投资费用还是所消耗的钢材重量，在整个过程设备中的比例都相当高。

例如，在化纤装置中，塔设备投资比例为44.9%；而在年产4.5 万吨丁二烯装置中，塔设备重量的比例高达54%之多。

随着石油、化工的迅速发展，塔设备的合理造型及设计将越来越受到关注和重视。

化工生产常需要进行液体混合物的分离已达到提纯或回收有用组分的目的。

互溶液体的分离有多种方法，精馏就是其中最常用的一种。

精馏塔是一种利用两组分的挥发度差异实现连续的高纯度分离的设备。

其中，回流是构成气、液两相接触传质的必要条件，也是精馏之区别于蒸馏所在。

本文设计的是脱丙烷精馏塔。

首先，根据已知的产品回收率进行了工艺计算，包括流程的确定、物料衡算、最小回流比的确定、最小理论塔板数的确定、塔板效率和实际塔板数的确定等。

然后对其结构进行了设计并得出流体力学计算结果。

综合以上设计及计算又得到塔的负荷性能图，以便对其性能有一个直观的了解。

本文最后对塔的附件进行设计并按照有关标准对其主要的部件进行强度和稳定性校核。

在完整地确定出结构和尺寸后，利用pore绘制了塔的主要零件图和塔的整体结构图。

关键词：精馏塔；工艺; 校核IAbstractTower is an important unit operation equipment in industries such as chemical engineering, o il refining, medicine, food and environmental protection. It is used widely. According to statistics, tower equipment, regardless of their investment costs or the amount of steel or the weight of equipment in the process, accounts for very . For example, in the fiber installations, the tower facility investment ratio is 44.9%. In an annual output of 45,000 tons of butadiene units, the ratio of the weight of tower equipment is as much as 54%.With the development of the petroleum, chemical industry is developing rapidly, reasonable design of the power will become more and more concerned.Chemical production often requires the separationof liquid mixtures that or recovery purposes. There are many ways of liquid separating. Distillation is one of the most commonly used. The use of distillation column is a two-point difference in the achievement of continuous volatility of the separation of constitutes a gas, liquid two-phase mass transfer contact with the necessary conditions for the distillation is distilled from theis designed in this article. First ofall, the basis of known products of the process is used to calculate the recovery rate. Including the identification process, material balance, the determination of the minimum reflux ratio, the minimum theoretical plate number of the identified tray efficiency and the actual determinationof the number plate. And then calculation of bear fruit designing carried out and reaching structure. The function designing and calculating the load getting a tower above synthesis is pursued. Finally, the tower accessories are designed and proofread according to carrying out the intensity and the stability on those main components in connection with the standard the main body of a book. Overall structural drawing in the picture and tower ascertaining out structure and the dimension queendrawn by making use of pore.Key Words：distillation; technology;check毕业设计（论文）原创性声明和使用授权说明原创性声明本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。