丙烯-丙烷体系对筛板塔顶精馏塔的设计（处理量：70回流比系数：1.6）

前言 本设计说明书包括概述、流程简介、精馏塔、再沸器、辅助设备、管路设计和控制方案共七章。 说明中对精馏塔的设计计算做了详细的阐述，对于再沸器、辅助设备和管路的设计也做了正确的说明。 鉴于设计者经验有限，本设计中还存在许多错误，希望各位老师给予指正 感谢老师的指导和参阅！目录第一章 概述1第二章 流程简介3第三章 精馏塔工艺设计5第四章 再沸器的设计19第五章 辅助设备的设计26第六章 管路设计34第七章 控制方案35设计心得及总结 36附件一 C程序 138附件二 matlab程序250附录一 主要符号说明51附录二 参考文献54第一章 概述 精馏是分离过程中的重要单元操作之一，所用设备主要包括精馏塔及再沸器和冷凝器。1 精馏塔精馏塔是一圆形筒体，塔内装有多层塔板或填料，塔中部适宜位置设有进料板。两相在塔板上相互接触时，液相被加热，液相中易挥发组分向气相中转移；气相被部分冷凝，气相中难挥发组分向液相中转移，从而使混合物中的组分得到高程度的分离。简单精馏中，只有一股进料，进料位置将塔分为精馏段和提馏段，而在塔顶和塔底分别引出一股产品。精馏塔内，气、液两相的温度和压力自上而下逐渐增加，塔顶最低，塔底最高。本设计为筛板塔，筛板的突出优点是结构简单、造价低、塔板阻力小且效率高。但易漏液，易堵塞。然而经长期研究发现其尚能满足生产要求，目前应用较为广泛。2 再沸器作用：用以将塔底液体部分汽化后送回精馏塔，使塔内气液两相间的接触传质得以进行。本设计采用立式热虹吸式再沸器，它是一垂直放置的管壳式换热器。液体在自下而上通过换热器管程时部分汽化，由在壳程内的载热体供热。立式热虹吸特点：循环推动力：釜液和换热器传热管气液混合物的密度差。 结构紧凑、占地面积小、传热系数高。壳程不能机械清洗，不适宜高粘度、或脏的传热介质。塔釜提供气液分离空间和缓冲区。3 冷凝器 （设计从略） 用以将塔顶蒸气冷凝成液体，部分冷凝液作塔顶产品，其余作回流液返回塔顶，使塔内气液两相间的接触传质得以进行，最常用的冷凝器是管壳式换热器。第二章 方案流程简介1 精馏装置流程 精馏就是通过多级蒸馏，使混合气液两相经多次混合接触和分离，并进行质量和热量的传递，使混合物中的组分达到高程度的分离，进而得到高纯度的产品。 流程如下： 原料（丙稀和丙烷的混合液体）经进料管由精馏塔中的某一位置（进料板处）流入塔内，开始精馏操作；当釜中的料液建立起适当液位时，再沸器进行加热，使之部分汽化返回塔内。气相沿塔上升直至塔顶，由塔顶冷凝器将其进行全部或部分冷凝。将塔顶蒸气凝液部分作为塔顶产品取出，称为馏出物。另一部分凝液作为回流返回塔顶。回流液从塔顶沿塔流下，在下降过程中与来自塔底的上升蒸气多次逆向接触和分离。当流至塔底时，被再沸器加热部分汽化，其气相返回塔内作为气相回流，而其液相则作为塔底产品采出。2 工艺流程1） 物料的储存和运输 精馏过程必须在适当的位置设置一定数量不同容积的原料储罐、泵和各种换热器，以暂时储存，运输和预热（或冷却）所用原料，从而保证装置能连续稳定的运行。2） 必要的检测手段 为了方便解决操作中的问题，需在流程中的适当位置设置必要的仪表，以及时获取压力、温度等各项参数。 另外，常在特定地方设置人孔和手孔，以便定期的检测维修。3） 调节装置由于实际生产中各状态参数都不是定值，应在适当的位置放置一定数量的阀门进行调节，以保证达到生产要求，可设双调节，即自动和手动两种调节方式并存，且随时进行切换。3 设备选用 精馏塔选用筛板塔，配以立式热虹吸式再沸器。4 处理能力及产品质量处理量： 70kmol/h产品质量：（以丙稀摩尔百分数计）进料：xf65塔顶产品：xD98塔底产品: xw2第三章 精馏塔工艺设计 第一节 设计条件1 工艺条件：饱和液体进料，进料丙稀含量xf65（摩尔百分数）塔顶丙稀含量 xD98，釜液丙稀含量 xw2，总板效率为0.6。 2操作条件：1）塔顶操作压力：P=1.62MPa（表压）2）加热剂及加热方法：加热剂水蒸气 加热方法间壁换热3）冷却剂：循环冷却水4）回流比系数：R/Rmin=1.6 3塔板形式：筛板 4处理量：qnfh=70kmol/h 5安装地点：大连 6塔板设计位置：塔顶 第二节 物料衡算及热量衡算一 物料衡算1换算将摩尔百分数换算成质量百分数：W=XMA/XMA+(1-X)MBxf65 wf63.93xD98 wD97.91xw2 wW1.91 将摩尔流量换算成质量流量：进料状态混合物平均摩尔质量：（MA为丙稀摩尔质量 MB为丙烷摩尔质量）M=xfMA+(1-xf)MB=0.65420.3544=42.7kg/kmol进料状态下的质量流量：qmfs=Mqnfh/3600=0.8303kg/s2求质量流量 qmDs + qmws = qmfsqmDswD + qmwswW = qmfswf解得： qmDs = 0.5364kg/s ; qmws= 0.2939kg/s 塔内气、液相流量：1）精馏段：L =RD; V =(R+1)D;2）提馏段：L=L+qF; V=V-(1-q)F; L=V+W; 二 热量衡算1）再沸器热流量：QR=Vr 再沸器加热蒸气的质量流量：GR= QR/rR2） 冷凝器热流量：QC=Vr冷凝器冷却剂的质量流量：GC= QC/(cl(t2-t1)第三节 塔板数的计算利用程序进行迭代计算：流程图如下：1泡点计算：泡点初值To总压PAntoinec常数outputy-1p0A p0BKA KB2塔板数计算：假设初值Tto、Tbo、2并输入Pt(绝)、xf 、xD 、xw 、qmfs 、qmDs 、qmws 、L计算1计算=(1+2)/2计算Rmin、R解得Nf、NPb=Pt0.98104N计算2计算并输出、Nf、N、R、qmLs、qmVs-(1+2)/2将2代入YN注：下标t、b分别表示塔顶、塔底参数。 计算过程包括：假设塔顶温度Tto=316K 经泡点迭代计算得塔顶温度Tt=316.13K塔顶压力Pt=1620+101.3=1721.3KPa 代入公式 计算并换算得：PAo=1726.955KPa ; PBo=1444.432KPa又 得：KA=1.003 ; KB=0.839BAKK=a1 1.1956设2 1.1888=(1+2)/2 1.1922计算过程包括：泡点进料：q=1 q线：x=xf 代入数据，解得 xe=0.65;ye=0.6889 =7.489 R=1.6Rmin=11.983为逐板计算过程：y1=xD=0.98ynynxn)1(-=aa直至xi xf 理论进料位置：第i块板进入提馏段：ynynxn)1(-=aa 直至xn xW 计算结束。理论板数：Nt=n（含釜）（具体程序见附件一）迭代结果：进料板Nf=i/0.6+1=52, 实际板数Np=(Nt-1)/0.6+1=112则塔底压力Pb=Pt+0.980.47Np= 1772.9KPa同可算得：塔底温度Tb=325.52K 2 1.1888 符合假设 所以假设成立，上述计算结果均为正确结果。塔内气、液相流量： 精馏段： qmLs=RqmDs=6.4275kg/s qmVs=(R+1)qmDs=6.9639kg/s 提馏段 ： q,mLs=qmLs+qmFs=7.2578kg/s q,mVs= qmVs =6.9639kg/s 第四节 精馏塔工艺设计1 物性数据常压43下，丙稀的物性数据：气相密度：V =28kg/ m3液相密度：L =470kg/ m3液相表面张力：=4.5mN/m2 初估塔径气相流量：qmVs=6.9639kg/s qVVs=qmVs/v=0.2487m3/s液相流量：qmLs=6.4275kg/s qVLs=qmLs/L=0.01368 m3/s两相流动参数： =0.2253初选塔板间距 HT=0.45m,查化工原理（下册）P107筛板塔泛点关联图，得：C20=0.061所以，气体负荷因子： =0.04526 液泛气速： 0.1798m/s 取泛点率0.7 操作气速：u = 泛点率 uf=0.1259 m/s 气体流道截面积： =1.9754 m2 选取单流型弓形降液管塔板，取Ad / AT=0.15; 则A / AT=1- Ad / AT =0.85 截面积: AT=A/0.85=2.324 m2 塔径： =1.72m 圆整后，取D=1.8m 符合化工原理书P108表6.10.1及P110表6.10.2的经验关联 实际面积： =2.5447 m2 降液管截面积：Ad=AT0.15= 0.3817 m2气体流道截面积：A=AT-Ad=2.1630 m2实际操作气速： = 0.1150 m/s 实际泛点率：u / uf =0.63943 塔高的估算 Np=112 有效高度：Z= HT Np=50.4m 釜液高度（略），进料处两板间距增大为0.7m设置8个人孔，每个人孔0.8m裙座取5m,塔顶空间高度1.5m,釜液上方气液分离高度取1.5m.设釜液停留时间为30min釜液高度：Z =0.44m 取其为0.5m 所以，总塔高h=Z+0.7+5+1.5+1.5+0.5=59.6m第五节 溢流装置的设计1 降液管（弓形）由上述计算可得：降液管截面积：Ad=AT0.15= 0.3817 m2由Ad/AT=0.15,查化工原理（下册）P113的图6.10.24可得：lw/D=0.82所以，堰长lw=0.82D=1.476m2 溢流堰 取E近似为1则堰上液头高： =0.029428m取堰高hw=0.04m,底隙hb=0.045m液体流经底隙的流速：ub =0.206m/sub0.5m/s 符合要求第六节 塔板布置和其余结构尺寸的选取取塔板厚度=3mm 进出口安全宽度bs=bs=80mm 边缘区宽度bc=50mm由Ad/AT=0.15,查化工原理（下册）P113的图6.10.24可得：bd/D=0.21所以降液管宽度：bd =0.21D=0.378m =8.542mr= =0.85m有效传质面积： = 1.4320 m2 取筛孔直径：do=7mm,取孔中心距：t=3.6do= 0.0252m开孔率： = =0.07 筛孔面积： = 0.1002 m2 筛孔气速： =2.4817m/s筛孔个数： =2604第七节 塔板流动性能校核1 液沫夹带量校核 Hf=2.5(hw+how)=0.1736m质量夹带率ev ： =0.000076kg液/kg气ev5s 满足要求 5 严重漏液校核 =0.0141m 满足稳定性要求 1.1257m/s第八节 负荷性能图1 过量液沫夹带线规定：ev = 0.1（ kg 液体 / kg气体） 为限制条件得： = 10671-167qVLh2/3 由上述关系可作得线2 液相上限线 整理出：qVLh=3.07lw=4.034 与y轴平行 由上述关系可作得线3 严重漏液线 将下式分别代入 近似取Co为前面计算的值 得：qVVh =a(b+cqVLh2/3)1/2其中：a= =4843.3 b=0.0056+0.13hw-h=0.01024 c= =0.000285得：qVVh =4843.3(0.01024+0.000285qVLh2/3)1/2 由上述关系可作得线4 液相上限线令 =5s 得： =90.693由上述关系可作得线5 浆液管液泛线Hd=HT+hW令 将 =0以及how与qVLh , hd 与qVLh ，hf 与qVVh , qVLh 的关系全部代入前式整理得： 式中：a= =42.6181109 b= =0.177 c= =267.4757108 d= =0.003724得： 上述关系可作得降液管液泛线上五条线联合构成负荷性能图作点为：qVLh =49.23m3/s qVVh =895.36 m3/s负荷性能图： 可见，线的位置偏上，所以它对操作的影响很小。放大后的负荷性能图：设计点位于四条线包围的区间中间稍偏下，操作弹性：qVVhmax / qVVhmin2.73所以基本满足要求（程序见附件二）第四章 再沸器的设计一 设计任务与设计条件 1选用立式热虹吸式再沸器 塔顶压力：1.7213MPa 压力降：Nphf=1120.1170.479.8103=0.0604MPa 塔底压力=1.7213+0.0604=1.7817MPa2再沸器壳程与管程的设计壳程管程温度（）100525压力（MPa绝压）0.10131.7817蒸发量：Db= q,mVs =6.9639kg/s3 物性数据1） 壳程凝液在温度（100）下的物性数据：潜热：rc=2257kj/kg热导率：c =0.683w/(m*K)粘度：c =0.283mPa*s密度：c =958.4kg/m32） 管程流体在（52.5 1.7817MPa）下的物性数据：潜热：rb=278.182kj/kg液相热导率：b =81.54mw/(m*K)液相粘度：b =0.071mPa*s液相密度：b =441.3kg/m3 液相定比压热容：Cpb= 3.090kj/(kg*k) 表面张力：b0.00377N/m气相粘度：v =0.071mPa*s气相密度：v =441.3kg/m3 蒸气压曲线斜率（t/P）=0.0000266 m2 K/kg 二 估算设备尺寸 热流量： = 1937235.522w 传热温差： =100-52.5=47.5K 假设传热系数：K=850W/( m2 K) 估算传热面积Ap =47.98 m2 拟用传热管规格为：382.5mm,管长L=3000mm 则传热管数： =134 若将传热管按正三角形排列，按式 NT =3a(a+1)+1;b=2a+1 得：b=13.552 管心距：t=0.048m 则 壳径： =0.7165m 取 D= L/D= 取 管程进口直径：Di=0.25m 管程出口直径：Do=0.3m 三 传热系数的校核1显热段传热系数K假设传热管出口汽化率 Xe=0.19则循环气量： =36.6522kg/s1) 计算显热段管内传热膜系数i 传热管内质量流速： di=38-22.5=33mm =0.118 = 310.687kg/( m2 s) 雷诺数： = 144403.9 普朗特数： =2.69058 显热段传热管内表面系数： = 1132.89w/( m2 K) 2) 壳程冷凝传热膜系数计算o 蒸气冷凝的质量流量： = 0.8583kg/s 传热管外单位润湿周边上凝液质量流量： =0.06kg/(m s) = 848.402122 管外冷凝表面传热系数： = 6547.988173w/ (m2 K) 3) 污垢热阻及管壁热阻 沸腾侧：Ri=0.000176 m2 K/w 冷凝侧：Ro=0.00021 m2 K/w 管壁热阻：Rw=b/w= 0.00005 m2 K/w 4）显热段传热系数 dm=(di+do)/2= 0.0355m = 611.4898w/( m2 K) 2. 蒸发段传热系数KE计算 传热管内釜液的质量流量：Gh=3600 qmws = 1118474.245kg/( m2 h) Lockhut-martinel参数： =1.1838 则1/Xtt=0.85 查设计书P96图329 得：E=0.2 在Xe=0.19 X3Xe=0.057的情况下 =0.2493 再查图329，=1 2）泡核沸腾压抑因数：=(E+)/2=0.6 泡核沸腾表面传热系数: =8104.605 w/( m2 K) 3）单独存在为基准的对流表面传热系数 ： = 1080.925w/( m2 K) 沸腾表面传热系数：KE 对流沸腾因子 ： = 1.7475227 两相对流表面传热系数: = 1888.94 w/( m2 K) 沸腾传热膜系数： = 6751.7 w/( m2 K) = 1266.69 w/( m2 K) 3.显热段及蒸发段长度 = 0.00319LBC = 0.00272L= 0.00957mLCD =L- LBC = 2.9904m4传热系数 = 1264.6 m2 实际需要传热面积： = 32.25 m25传热面积裕度： = 0.48780.150.2所以，传热面积裕度合适，满足要求四 循环流量校核1循环系统推动力：1）当X=Xe/3= 0.063333333时=3.6264 两相流的液相分率： = 0.3816 两相流平均密度： = 186.96kg/m3 2）当X=Xe=0.19 = 1.1838两相流的液相分率： = 0.2267 两相流平均密度： = 123.26kg/m3根据课程设计表319 得：L=1.02m, 则循环系统的推动力： = 6228.48pa 2循环阻力Pf： 管程进出口阻力P1 进口管内质量流速： = 747.05kg/(m2s)釜液进口管内流动雷诺数： = 2630462进口管内流体流动摩擦系数： = 0.0150进口管长度与局部阻力当量长度： =29.30m管程进出口阻力: =1112.343Pa 传热管显热段阻力P2 =310.687kg/(m2s) =144404 =0.0205 = 0.6512Pa 传热管蒸发段阻力P3 a. 气相流动阻力Pv3 =39.3537kg/(m2s) =18291.167 =0.030 =71.06Pab. 液相流动阻力PL3GL=G-Gv=271.33 kg/(m2s) = 126112 = 0.0210 = 158.464Pa = 1732.2382Pa 管内动能变化产生阻力P4 动量变化引起的阻力系数: = 2.58 = 564.423管程出口段阻力P5 a. 气相流动阻力Pv5 = 518.786kg/(m2s) = 98.569kg/(m2s) 管程出口长度与局部阻力的当量长度之和： = 35.0425m = 3333797.217 = 0.014775112 = 279.4712b. 液相流动阻力PL5 =212.118 kg/(m2s) = 896273 = 0.0164 = 97.6351Pa = 2735.64Pa 所以循环阻力：Pf=P1 + P2 + P3 + P4 + P5 = 6145.292 又因PD=6228.4845Pa 所以 =1.0135 第五章 辅助设备设计 一 辅助容器的设计 容器填充系数取：k=0.7 1进料罐（常温贮料） 20丙稀 L1 =526kg/m3 丙烷 L2 =500kg/m3 压力取1.818MPa 由上面的计算可知 进料 Xf=65% Wf=63.93% 则 =516.3 kg/m3 进料质量流量：qmfh=3600 qmfs=2989.08kg/h 取 停留时间：x为4天，即x=96h 进料罐容积： 793.978m3 圆整后 取V=794 m3 2回流罐（43）质量流量qmLh=3600RqmDs =23139kg/h设凝液在回流罐中停留时间为0.25h，填充系数=0.7则回流罐的容积 17.39 m3取V=18 m33塔顶产品罐质量流量qmDh=3600qmDs =1931.04 kg/h；产品在产品罐中停留时间为72h，填充系数=0.7则产品罐的容积 4225.98m3取V=4226m34 釜液罐取停留时间为5天，即x=120h质量流量qmWh=3600qmWs =1058.04 kg/h 则釜液罐的容积 404.86 m3取V=405m3二 传热设备 1进料预热器 用90水为热源，出口约为70走壳程 料液由20加热至45，走管程传热温差： 管程液体流率：qmfh=3600 qmfs=2989.08kg/h 管程液体焓变：H=401kj/kg 传热速率：Q= qmfsH=2989.08401/3600=332.99kw 壳程水焓变：H=125.6kj/kg 壳程水流率：q=9050kg/h 假设传热系数：K=650w/(m2K) 则传热面积： 圆整后取A=12m2 2 塔顶冷凝器拟用10水为冷却剂，出口温度为30。走壳程。管程温度为43管程流率：qmVs=6.9639kg/s取潜热r=302.54kj/kg传热速率：Q= qmVsr=2106.86kw壳程取焓变：H=125.5kj/kg则壳程流率：qc=Q/H=60435.82kg/h假设传热系数：K=650 w/(m2K)则传热面积： 圆整后 取A=152m23 塔顶产品冷却器拟用10水为冷却剂，出口温度为20。走壳程。管程温度由43降至25 管程流率：qmDs = 0.5364kg/s ; 取潜热：r=280kj/kg则传热速率：Q= qmDsr=150.192kw壳程焓变：H=84.0kj/kg则壳程流率：qc=Q/H=6437.14kg/h假设传热系数：K=650 w/(m2K)则传热面积 圆整后 取A=16 m24 釜液冷却器拟用10水为冷却剂，出口温度为20。走壳程。管程温度由52.5降到25管程流率：qmWs=0.2939kg/s丙烷液体焓变：H =282kj/kg传热速率：Q= qmVsH =82.88kw壳程取焓变：H=84.0kj/kg则壳程流率：qc=Q/H=3552kg/h假设传热系数：K=650 w/(m2K)则传热面积： 圆整后 取A=8 m2三 泵的设计1进料泵(两台，一用一备)取液体流速：u=0.5m/s液体密度： kg/ m3 qVfs = qmfs / =0.0016 m3/s 取d=65mm液体粘度 取=0.2相对粗糙度：/d=0.003查得：=0.025取管路长度：l=100m 取90度弯管4个，截止阀一个，文氏管流量计1个取则qVLh =5.97m3/h选取泵的型号：AY 扬程：3065m 流量：2.560m3 /s2回流泵（两台，一开一用）取液体流速：u=0.5m/s液体密度： kg/ m3 qVLs = qmLs / =0.0137 m3/s 液体粘度 取=0.2相对粗糙度：/d=0.00107查得：=0.02取管路长度：l=100m 取90度弯管4个，截止阀一个，文氏管流量计1个取则qVLh =49.44m3/h选取泵的型号：DSJH 扬程：38280m 流量：951740m3 /s3.釜液泵（两台，一开一用）取液体流速：u=0.4m/s液体密度： kg/ m3 qVWs = qmWs / =0.00089m/s 液体粘度 取=0.2相对粗糙度：/d=0.0045查得：=0.04取管路长度：l=40m 取90度弯管4个，截止阀一个，文氏管流量计1个取则qVLh =2.393m3/h该处泵扬程为负值，正常工作时不使用，但非正常工作或停止工作时，需要使用。选取泵的型号：GI 扬程：101510m 流量：0.190m3 /s第六章 管路设计进料管线取料液流速：u=0.5m/s则取管子规格685。其它各处管线类似求得如下：名称管内液体流速（m/s）管线规格（mm）进料管0.5703顶蒸气管1532510顶产品管0.5603回流管0.51803釜液流出管0.5324.5仪表接管/252.5塔底蒸气回流管151594第七章 控制方案 精馏塔的控制方案要求从质量指标、产品产量和能量消耗三个方面进行综合考虑。精馏塔最直接的质量指标是产品浓度。由于检测上的困难，难以直接按产品纯度进行控制。最常用的间接质量指标是温度。 将本设计的控制方案列于下表序号位置用途控制参数介质物性L(kg/m3)1FIC-01进料流量控制03000kg/h丙烷、丙稀L=516.32FIC-02回流定量控制01500kg/h丙稀L=469.83PIC-01塔压控制02MPa丙稀V=284HIC-02回流罐液面控制01m丙稀L=4705HIC-01釜液面控制03m丙烷L=4426TIC-01釜温控制4060丙烷L=442总结 在写这篇总结的时候我有了一种如释重负的感觉，说实话，两周的设计，其过程是痛苦和充满曲折的，在进行各种计算以及参数选择的时候，常会进入死胡同，看似无法解决，总会有一种濒临崩溃的感觉，但没有选择，只能硬着头皮做下去。问题是总会解决的，只要你付出努力，当你的迷茫达到一定的时候，就必然会看到一丝曙光。然后你可以沿着光的方向向前一步步艰难的迈进。也许每天，你的设计只会有一丝的进展，甚至你会发现某一天的努力都是徒劳的，因为你选择了一个错误的方向。面对这种情况，开始我自然是很困扰的，但后来我发现，即使缓慢，即使是错误，那些付出依然是有价值的。试了，错了，然后才会知道这样是错的，因此对于设计的每一步，我的映象都十分的深刻。 我可以自信的说，化原考试的话，我可以考到A，甚至A。然而，到了真正做设计的时候，才发现自己真的知之甚少，有时候甚至觉得无从下手。当设计终于做完的时候，我可以肯定，其中必定充满了很多很多的错误，但我完全可以坦然面对这些错误，因为进步正是在错了再改，一改再改的前提下产生的。经过这次课程设计，我深刻的体会到：从书本上的理论知识到真正的生产实践，期间的距离真是差了很远。现在我们是作设计，已经觉得很困难，到了下工厂操作的时候，必然又会遇上新的问题。但我们从来就是不惧怕困难的，在不断的征服困难的过程中，我们才能也必然会掌握这门技术。 这次课程设计完成后，我发现我对于化工原理知识的了解上升到了一个新的层面，对于设计过程中的每一步，我都能说出它的原理和具体做法。对于上课时涉及较少的工艺流程也熟悉了不少。此外，在做设计的过程中复习并掌握了许多计算机知识，例如C语言，EXCEL，MATLAB，AUTO-CAD等。总之，通过这次课程设计，丰富了我各个方面的知识，我受益匪浅。更希望各位老师能帮助指出我设计中的错误与不足之处，使我能不断提高进步。姓名：陈晓霞 班级：化机0104班 学号：200142044 指导教师：贺高红 韩志忠附件一 C程序 1#include /* 此头函数请不要删除 */#include #include int i,j,n,N,Nf;double T,pa,fn,fnd,x0,m,s,a,xe,ye,xd,xw,R,Rmin,pb,a1,a2,qmd,qmw,qmf,aa,dl,dv,Flv,Ht,Hf; static double p2=1,1;static double K2=1,1; static double X200; static double Y200;static double qml2;static double qmv2;static double A2=15.7027,15.7260;static double B2=1807.53,1872.46;static double C2=-26.15,-25.16;void pi(int i,double T) pi=exp(Ai-Bi/(T+Ci)*101.3/760; void Ki(int i,double pa) Ki=pi/pa;void sumy(double x0) s=K0*x0+K1*(1-x0); void test(double T,double pa,double x0) fn=0;fnd=1; do T=T-fn/fnd; for(i=0;i0.001); printf(T=%f ,T); printf(pa=%f ,pa);for(i=0;i2;i+) printf(p=%f ,pi);printf(K=%fn,Ki); a=K0/K1; printf(a=%f ,a); main()double f,c,c20,uf,u,A,At,Ad,b,Dc,D,hb,hw,lw,how,bs,bc,bd,Aa,x,r,fi,d0,Ao,uo,uo1,t,ev,co,g,ho,hl,be,ua,Fa,hs,hf,FI,hd,Hd,Hd1,to,ho1,k,d; test(316,1721.3,0.98);a1=a;doa2=a;a=(a1+a2)/2;xd=0.98;xw=0.02;qmd=0.5364;qmw=0.2939;qmf=0.8303;dl=470;dv=28,Ht=0.4