多效蒸发工艺设计计算

多效蒸发工艺设计计算多效蒸发工艺设计计算PAGE多效蒸发工艺设计计算目录第一章前言§1·1概述`第二章蒸发工艺设计计算§2·1蒸浓液浓度计算§2·2溶液沸点和有效温度差的确定§2·2·1各效由于溶液的蒸汽压下降所引起的温度差损失/§2·2·2各效由于溶液静压强所因引起的温度差损失§2·2·3由经验不计流体阻力产生压降所引起的温度差损失§2·3加热蒸汽消耗量和各效蒸发水量的计算§2·4蒸发器的传热面积和有效温度差在各效中的分布以及传热系数K的确定§2·5温差的重新分配与试差计算§2·5·1重新分配各效的有效温度差，§2·5·2重复上述计算步骤§2·6计算结果列表第三章NaOH溶液的多效蒸发优化程序部分§3·1具体的拉格朗日乘子法求解过程§3·2程序内部变量说明§3·3程序内容：§3·4程序优化计算结果§3·5优化前后费用比较第四章蒸发器工艺尺寸计算§4·1加热管的选择和管数的初步估计§4·1·1加热管的选择和管数的初步估计§4·1·2循环管的选择§4·1·3加热室直径及加热管数目的确定§4·1·4分离室直径与高度的确定§4·2接管尺寸的确定§4·2·1溶液进出§4·2·2加热蒸气进口与二次蒸汽出口§4·2·3冷凝水出口第五章、蒸发装置的辅助设备§5·1气液分离器§5·2蒸汽冷凝器§5·2·1冷却水量§5·2·2计算冷凝器的直径§5·2·3淋水板的设计§5·3泵选型计算§5·4预热器的选型第六章主要设备强度计算及校核§6·1蒸发分离室厚度设计§6·2加热室厚度校核第七章小结与参考文献： 符号说明希腊字母：c——比热容，KJ/　　　　　α――对流传热系数，Ｗ/m2.℃　　　d——管径，m　　　　　　　　　Δ――温度差损失，　℃　　　　　　　D——直径，m　　　　　　　　　η――误差，D——加热蒸汽消耗量，Kg/h　　η――热损失系数，f——校正系数，　　　　　　　　η――阻力系数，F——进料量，Kg/h　　　　　　　λ――导热系数，Ｗ/m2.℃g——重力加速度，s2　　　μ――粘度，h——高度，m　　　　　　　　　ρ――密度，Kg/m3H——高度，m　　　　　　　　　　　　k——杜林线斜率　　　　　　　　　　K——总传热系数，W/m2.℃　　　∑――加和L——液面高度，m　　　　　　　φ――系数L——加热管长度，mL——淋水板间距，m下标：n——效数　　　　　　　　　　　　1,2,3――效数的序号　　　　　　　　　　　　　　　　　　　n——第n效　　　　　　　　　　　0――进料的　　p——压强，Pa　　　　　　　　　　i――内侧q——热通量，W/m2　　　　　　　m――平均Q——传热速率，W　　　　　　　o――外侧r——汽化潜热，KJ/Kg　　　　　p――压强　　　R——热阻，m2.℃/W　　　　　　　s――污垢的S——传热面积，m2　　　　　　　w――水的　t——管心距，m　　　　　　　　　w――壁面的　T——蒸汽温度，℃u——流速，m/sU——蒸发强度，Kg/　　　　　上标：V——体积流量，m3/h　　　　　　　′：二次蒸汽的W——蒸发量，Kg/h　　　　　　′：因溶液蒸汽压而引起的W——质量流量，Kg/h〞：因液柱静压强而引起的x——溶剂的百分质量,％　　　　　：因流体阻力损失而引起的          第一章前言§1·1概述1蒸发及蒸发流程蒸发是采用加热的方法，使含有不挥发性杂质（如盐类）的溶液沸腾，除去其中被汽化单位部分杂质，使溶液得以浓缩的单元操作过程。蒸发操作广泛用于浓缩各种不挥发性物质的水溶液，是化工、医药、食品等工业中较为常见的单元操作。化工生产中蒸发主要用于以下几种目的：1获得浓缩的溶液产品；2、将溶液蒸发增浓后，冷却结晶，用以获得固体产品，如烧碱、抗生素、糖等产品；3、脱除杂质，获得纯净的溶剂或半成品，如海水淡化。进行蒸发操作的设备叫做蒸发器。蒸发器内要有足够的加热面积，使溶液受热沸腾。溶液在蒸发器内因各处密度的差异而形成某种循环流动，被浓缩到规定浓度后排出蒸发器外。蒸发器内备有足够的分离空间，以除去汽化的蒸汽夹带的雾沫和液滴，或装有适当形式的除沫器以除去液沫，排出的蒸汽如不再利用，应将其在冷凝器中加以冷凝。蒸发过程中经常采用饱和蒸汽间壁加热的方法，通常把作热源用的蒸汽称做一次蒸汽，从溶液蒸发出来的蒸汽叫做而次蒸汽。2蒸发操作的分类按操作的方式可以分为间歇式和连续式，工业上大多数蒸发过程为连续稳定操作的过程。按二次蒸汽的利用情况可以分为单效蒸发和多效蒸发，若产生的二次蒸汽不加利用，直接经冷凝器冷凝后排出，这种操作称为单效蒸发。若把二次蒸汽引至另一操作压力较低的蒸发器作为加热蒸气，并把若干个蒸发器串联组合使用，这种操作称为多效蒸发。多效蒸发中，二次蒸汽的潜热得到了较为充分的利用，提高了加热蒸汽的利用率。按操作压力可以分为常压、加压或减压蒸发。真空蒸发有许多优点：（1）、在低压下操作，溶液沸点较低，有利于提高蒸发的传热温度差，减小蒸发器的传热面积；（2）、可以利用低压蒸气作为加热剂；（3）、有利于对热敏性物料的蒸发；（4）、操作温度低，热损失较小。在加压蒸发中，所得到的二次蒸气温度较高，可作为下一效的加热蒸气加以利用。因此，单效蒸发多为真空蒸发；多效蒸发的前效为加压或常压操作，而后效则在真空下操作。3蒸发操作的特点从上述对蒸发过程的简单介绍可知，常见的蒸发时间壁两侧分别为蒸气冷凝和液体沸腾的传热过程，蒸发器也就是一种换热器。但和一般的传热过程相比，蒸发操作又有如下特点:(1)   沸点升高蒸发的溶液中含有不挥发性的溶质，在港台压力下溶液的蒸气压较同温度下纯溶剂的蒸气压低，使溶液的沸点高于纯溶液的沸点，这种现象称为溶液沸点的升高。在加热蒸气温度一定的情况下，蒸发溶液时的传热温差必定小于加热唇溶剂的纯热温差，而且溶液的浓度越高，这种影响也越显著。(2)   物料的工艺特性蒸发的溶液本身具有某些特性，例如有些物料在浓缩时可能析出晶体，或易于结垢；有些则具有较大的黏度或较强的腐蚀性等。如何根据物料的特性和工艺要求，选择适宜的蒸发流程和设备是蒸发操作彼此必须要考虑的问题。(3)   节约能源蒸发时汽化的溶剂量较大，需要消耗较大的加热蒸气。如何充分利用热量，提高加热蒸气的利用率是蒸发操作要考虑的另一个问题。4蒸发设备  蒸发设备的作用是使进入蒸发器的原料液被加热，部分气化，得到浓缩的完成液，同时需要排出二次蒸气，并使之与所夹带的液滴和雾沫相分离。  蒸发的主体设备是蒸发器，它主要由加热室和蒸发室组成。蒸发的辅助设备包括：使液沫进一步分离的除沫器，和使二次蒸气全部冷凝的冷凝器。减压操作时还需真空装置。兹分述如下： 由于生产要求的不同，蒸发设备有多种不同的结构型式。对常用的间壁传热式蒸发器，按溶液在蒸发器中的运动情况，大致可分为以下两大类：（1）循环型蒸发器  特点：溶液在蒸发器中做循环流动，蒸发器内溶液浓度基本相同，接近于完成液的浓度。操作稳定。此类蒸发器主要有  a.中央循环管式蒸发器，  b.悬筐式蒸发器c.外热式蒸发器，d.列文式蒸发器e.强制循环蒸发器。其中，前四种为自然循环蒸发器。（2）单程型蒸发器  特点：溶液以液膜的形式一次通过加热室，不进行循环。  优点：溶液停留时间短，故特别适用于热敏性物料的  蒸发；温度差损失较小，表面传热系数较大。  缺点：设计或操作不当时不易成膜，热流量将明显下降；不适用于易结晶、结垢物料的蒸发。  此类蒸发器主要有    a.升膜式蒸发器，    b.降膜式蒸发器，    c.刮板式蒸发器 本次设计采用的是中央循环管式蒸发器：结构和原理：其下部的加热室由垂直管束组成，中间由一根直径较大的中央循环管。当管内液体被加热沸腾时，中央循环管内气液混合物的平均密度较大；而其余加热管内气液混合物的平均密度较小。在密度差的作用下，溶液由中央循环管下降，而由加热管上升，做自然循环流动。溶液的循环流动提高了沸腾表面传热系数，强化了蒸发过程。这种蒸发器结构紧凑，制造方便，传热较好，操作可靠等优点，应用十分广泛，有"标准蒸发器"之称。为使溶液有良好的循环，中央循环管的截面积，一般为其余加热管总截面积的40%～100%；加热管的高度一般为1～2m；加热管径多为25～75mm之间。但实际上，由于结构上的限制，其循环速度一般在～s以下；蒸发器内溶液浓度始终接近完成液浓度；清洗和维修也不够方便。   第二章蒸发工艺设计计算§2·1蒸浓液浓度计算多效蒸发的工艺计算的主要依据是物料衡算和、热量衡算及传热速率方程。计算的主要项目有：加热蒸气（生蒸气）的消耗量、各效溶剂蒸发量以及各效的传热面积。计算的已知参数有：料液的流量、温度和浓度，最终完成液的浓度，加热蒸气的压强和冷凝器中的压强等。蒸发器的设计计算步骤多效蒸发的计算一般采用试算法。（1）      根据工艺要求及溶液的性质，确定蒸发的操作条件（如加热蒸气压强及冷凝器的压强），蒸发器的形式、流程和效数。（2）      根据生产经验数据，初步估计各效蒸发量和各效完成液的浓度。（3）      根据经验假设蒸气通过各效的压强降相等，估算个效溶液沸点和有效总温差。（4）      根据蒸发器的焓衡算，求各效的蒸发量和传热量。（5）      根据传热速率方程计算各效的传热面积。若求得的各效传热面积不相等，则应按下面介绍的方法重新分配有效温度差，重复步骤（3)至（5），直到所求得各效传热面积相等（或满足预先给出的精度要求）为止。F==h总蒸发量:W=F(1-)=(1-)=h并流加料蒸发中无额外蒸汽引出,可设:W2:W3=1::而W=W1+W2+W3=h由以上三式可得:W1=h;W2=h;W3=h;X1==;X2==;X3=§2·2溶液沸点和有效温度差的确定设各效间的压强降相等，则总压强差为：=P1-PK/=P=式中P---各效加热蒸汽压强与二次蒸气压强之差KPa，----第一次加热蒸气的压强KPa-----末效冷凝器中的二次蒸气的压强KPa各效间的压强差可求得各效蒸发室的压强即P1/=P1-Pi=P2/=P1-2Pi=P3/=Pk/= 由各效的二次蒸汽压强从手册中查得相应的二次蒸汽温度和汽化潜热列与下表中: 第一效第二效第三效二次蒸汽压强Pi/(KPa)   二次蒸汽温度Ti/(℃)(即下一效加热蒸汽温度)   二次蒸汽的汽化潜热(即下一效加热蒸汽的ri/)   多效蒸发中的有效传热总温度差可用下式计算： 有效总温度差式中-----有效总温度差，为各效有效温度差之和，℃。-----第一效加热蒸气的温度，℃。-----冷凝器操作压强下二次蒸气的饱和温度，℃。-------总的温度差损失，为各效温度差损失之和，℃，=/+在焓衡算式计入溶液的能缩热及蒸发器的热损失时，尚需考虑热利用系数一般溶液的蒸发，可取得（式中△x为溶液的浓度变化，以质量分率表示）。第i效的蒸发量的计算式为式中------第i效的加热蒸汽量，当无额外蒸汽抽出时=------第i效加热蒸气的汽化潜热------第i效二次蒸气的汽化潜热-----------原料液的比热---------水的比热，--------分别为第i效及第i-1效溶液的沸点-----------第i效的热利用系数无因次，对于加热蒸气消耗量，可列出各效焓衡算式并与式（3-2）联解而求得。第一效的焓衡量式为:W1=由相关手册查得cp0=(Kg.℃)，cpw=KJ/(Kg.℃)代入W1==(a)同理第二效的热衡算式为:W2=]所以W2=第三效的热衡算式为:W3=W3==+(c)又W1+W2+W3=(d)联立(a),(b),(c),(d)式,解得:W1=hW2=hW3=hD1=Kg/h§2·4蒸发器的传热面积和有效温度差在各效中的分布以及传热系数K的确定任意一效的传热速率方程为Si=式中---第i效的传热速率，W。----第i效的传热系数，W/（m2,℃）.---第i效的传热温度差，℃Si-------第i效的传热面积，m2在三效蒸发中，为了便于制造和安装，通常采用各效传热面积相等的蒸发器，即S1=S3=S3=S若求得的传热面积不相等，应根据各效面积相的原则重新分配各校的有效温度差。方法如下：总的传热系数K；称为公式（A）式中，---管外蒸汽冷凝传热系数与管内液体沸腾传热系数，----管外和管内的污垢热阻，-----管壁的导热系数，----管壁的厚度，，，-----加热管的外径、平均直径与内径，计算K0值主要在于求取关内溶液的沸腾传热系统。该值受溶液的性质、蒸发器的类型、沸腾传热的形式伊基蒸发操作的条件等许多因素的影响。因此，一般沸腾传热膜系数关联式的准确度较差。Q1=D1r1==h=T1-t1=查得RS0=×10-4m2.K/W,RSi=×10-4m2K/W,取规格为的加热管则d0=38mm,di==33mm,dm=,b=第一效：蒸汽在垂直管外的冷凝对流传热系数假设壁温tw1=150℃，特性温度t=（+150）/2=℃，==℃，由t值查得冷凝液的导热系数=，密度=m3，粘度=×10-3，又=Kg1===；验算:Re==<2000，所以1=由经验公式:(*)其中,：---沸腾液体与其蒸汽间的表面压力,N/m;----液体的导热系数,W/;P-----二次蒸汽压力,Pa;cL------液体的比热,J/;r-----液体的汽化潜热J/Kg;----液体的密度Kg/m3;-----蒸汽的密度Kg/m3-----液体的动力粘度Pas;-----水的动力粘度;-----管壁与沸腾液体间的传热温差K;C-----与表观液面高度百分比h有关的系数;h:-----加热管内含纯液柱高度占总管长的百分比;对一效查得：=N/m，P=Pa，C=，cL=;=;r=Kg;=m3;=m3;=×;=×10-3Pas;==;带入上式解得：1=带入公式（A）得K1=检验，带入数据解得tw=℃，与假设基本相同，故不必重算。所以K1=第二效：假设tw=135℃，特性温度t=（135+）/2=℃，==℃，由t值查得冷凝液的导热系数=，密度=m3，粘度=×10-3Pas，又=Kg由公式：2=得2==；验证：Re==<2000，所以2=由经验公式:因为t=（135+）/2=℃对二效查得：=N/m，P=KPa，C=，cL=;=;r=Kg;=m3;=m3;=×;=×10-3Pas;==K;带入上式解得：2=带入公式（A）得K2=检验，带入数据解得tw=℃，重新假设tw=℃，所以t=+/2=℃△    t=由t=℃查得=，密度=m3，粘度=×10-3Pas，又=Kg根据式2=得：2=；验证：验证：Re==<2000，所以2=因为t=（+）/2=℃对二效查得：=N/m，P=KPa，C=，cL=;=;r=Kg;=m3;=m3;=×;=×10-3Pas;=根据公式：得2=由公式（A）计算得到K2=检验带入数据解得tw=℃，与假设基本相同，故不必重算。所以K2=第三效根据式3=，Re<2000，验证：Re=假设tw3=110℃，特性温度t=（+110）/2=℃，==℃，由t值查得冷凝液的导热系数=，密度=m3，粘度=×10-3Pas，又=Kg由公式：3=得3==；验证：Re==<2000，所以2=由t=（110+）/2=℃对三效查得：=N/m，P=KPa，C=，cL=;=;r=Kg;=1275Kg/m3;=m3;=×;=×10-3Pas;==;经验公式:带入上式解得：3=带入公式（A）得K3=检验带入数据解得tw=℃，与假设基本相同，故不必重算。所以K3=即K1=,K2=,K3=。因为参数太多，值查的不准，计算的K值差别不大，优化效果不明显，故用经验值，K1=1860,K2=1280,K3=700分别带入公式：Si=得到S1==m2S2==S3==m2误差为：1--误差较大，故应调整各效的有效温度差，重复上述步骤。§2·5温差的重新分配与试差计算§2·5·1重新分配各效的有效温度差，平均传热面积:S==重新分配有效温度差:==℃==℃==℃§2·5·2重复上述计算步骤：§2·5·2·1由所求各效蒸汽量求各效溶液浓度X1=;X2=，X3=，§2·5·2·2计算各效溶液沸点末效溶液沸点和二次蒸汽压强保持不变,各种温度差损失可视为衡值，故末效溶液的沸点=.而=℃,则第三效加热蒸汽温度（即第二效二次蒸汽温度）为T3=T2/=+=℃,由于第二效二次蒸汽温度为T2/=℃,再由X2=查杜林曲线得:tA2=℃,所以t2=tA2+=+=℃同理：由t2=℃，=℃则T2=T1/=t2+=+=℃再由T1/=℃，X1=，查杜林曲线得:tA1=℃t1=+=℃；也可以由t1=T1=℃说明溶液的温度差损失变化不打，不必重新计算，故有效温度差为℃以上计算结果总结如下:效数第一效第二效第三效溶液沸点ti℃加热蒸汽的温度(℃)有效温度差(℃)§2·5·2·3各效焓衡算:二次蒸汽的汽化潜热和二次蒸汽温度Ti/(℃)(即下一效加热蒸汽)温度如下表：效数第一效第二效第三效二次蒸汽温度Ti/(℃)(即下一效加热蒸汽)温度二次蒸汽的汽化潜热(即下一效加热蒸汽的ri/)第一效：=－×（－）=W1=h第二效：W2=第三效：=－×－=W3=5[]Kg/hW3=－W1+07705W2因为：W1+W2+W3=h得：W1=h，W2=hW3=Kg/h，D=h与第一次热量恒算所得结果：W1=h，W2=h，W3=h比较并计算误差得=，=，=相对误差均在误差允许范围之内股计算得各效蒸发面积合理。其各效溶液浓度无明显变化不必再算。§2·5·2·4计算各效传热面积S1==m2S2==m2S3==m2误差计算得：1-=1－=<；所以误差允许取平均面积的S=§2·6计算结果列表效数123冷凝器加热蒸汽温度(℃)操作压强Pi/(KPa)溶液沸点ti℃ 完成液浓度(%) 蒸发水量WiKg/h 生蒸汽量DKg/h   传热面积Sim2 表中Pi/按T1/=T2=℃查得,P2/按T2/=T3=℃查得,Pk/≈P3/=第三章NaOH溶液的多效蒸发优化程序部分§3·1具体的拉格朗日乘子法求解过程具体的拉格朗日乘子法求解过程如下：（1）先输入笔算得相关数据（与优化计算的过程相比可视为常量的或者不变的量）a0=b0=,B=40000，c0=4180,c1=,c2=1,c=0，e=,f=,fc=，k1=1860,k2=1280,k3=700，o=7680,p=1000,R0=,r3=,y=,x0=,x3=,t0=20,t3=,T1=,Tk=,Tpm[3]={,,},Ta[3]={,,},T[3]={,,};（2）      根据公式w=f*(1-x0/x3);计算总的蒸发量（3）      根据Tpm[3]，Ta[3]，T[3]德具体数值计算相关量b视为不便的量。（4）      在输入笔算的第一效和第二效蒸气的汽化潜热为程序内部循环作准备（5）      根据公式w1=w/(1+r1/r2+r1/r3);w2=w/(1+r2/r1+r2/r3);w3=w-w1-w2;D1=w3*r3/R0;计算分别算出各效蒸发量并求出生蒸汽的量（6）      根据公式x1=f*x0/(f-w1);x2=f*x0/(f-w1-w2);进一步求出第一效和第二效的浓度（7）      根据根据拉格朗日乘子法，求出各效温度差损失dt1、dt2、dt3。(8)利用经验公式R1=R2=求变化后的第一效和第二效蒸气的汽化潜热。当满足条件时退出循环并开始计算面积，同时求出费用，如果不满足条件就自动传递变量第一效和第二效蒸气的汽化潜热，重新进入循环直至满足条件（9）分别计算年生蒸汽费用，蒸发器年折旧维修费用和整孔系统费用。并考虑节省的费用。§3·2程序内部变量说明程序内部变量说明：a0------计算年蒸汽费用时的常量,a--------双精度变量,A[3]----面积的一维变量,b------计算年蒸汽费用时的常量,b--------双精度变量,B--------蒸发器用不锈钢市场钢材价格,40000,c0--------水的比热容,c1---------电费,c2-------动力费附加系数,c----------双精度变量,,d---------双精度变量,D1-------生蒸汽-双精度变量dt1,dt2,dt3,------分别为温度差变量e------循环控制常量,E---------双精度变量,f----------溶液流量,fc-------蒸发器念折旧及维修率,h---------双精度变量,H[3]--------增资系数的一维变量,k1,k2,k3-------各效总的传热系数L---------双精度变量,M---------双精度变量,N---------双精度变量,m---------双精度变量,n---------双精度变量,o---------蒸发器年工作时间7680,p---------水的通常密度,R0---------生蒸汽的汽化潜热,r1,r2,r3--------各效汽化潜热的双精度变量,w-------总的蒸发量,w1,w2,w3-------各效蒸发量,vw--------双精度变量,y----------泵的效率,j,j1,j2,j3-------------费用的双精度变量，t0---------------冷凝水的进口温度20,t1,t2,t3-------------各效的溶液的沸点,T生蒸汽的温度,Tk---------------冷凝器的蒸汽温度,R1,R2-------第一效和第二效蒸气的汽化潜热双精度变量x0,x1,x2,x3---------原料液浓度和各效内部的浓度,Tpm[3]---------溶液静平均压强对应的饱和温度的一维变量,Ta[3]-----------由杜林曲线查出的各小溶液沸点的一维变量,T[3]-------------二次蒸汽温度的一维变量;§3·3程序内容：#include<>#include<>main(){inti;doublea0=,a,A[3],b0=,b,B=40000,c0=4180,c1=,c2=1,c=0,d=0,D1,dt1,dt2,dt3,e=,E,f=,fc=,h,H[3],k1=1860,k2=1280,k3=700,L,M=0,N,m,n,o=7680,p=1000,R0=,R1,R2,r1,r2,r3=,w,w1,w2,w3,vw,y=,j,j1,j2,j3,x0=,x1,x2,x3=,t0=20,t1,t2,t3=,T1=,Tk=,Tpm[3]={,,},Ta[3]={,,},T[3]={,,};for(i=0;i<3;i++){c+=Ta[i]-T[i];d+=Tpm[i]-T[i];}b=T1-Tk-c-d;scanf("%f,%f",&r1,&r2);w=f*(1-x0/x3);for(;;){w1=w/(1+r1/r2+r1/r3);w2=w/(1+r2/r1+r2/r3);w3=w-w1-w2;D1=w3*r3/R0;x1=f*x0/(f-w1);x2=f*x0/(f-w1-w2);m=D1*R0/k1;n=w1*r1/k2;h=w2*r2/k3;a=sqrt(m)+sqrt(n)+sqrt(h);dt1=b*sqrt(m)/a;dt2=b*sqrt(n)/a;dt3=b*sqrt(h)/a;R1=t1=t3+dt3+dt2;R2=if(fabs((R1-r1)/r1)<=e&&fabs((R2-r2)/r2)<=e)break;else{r1=R1;r2=R2;}}t1=T1-dt1;t2=t3+dt3;printf("xiangguanshuju:\n");printf("D1=%f,\nw=%f,w1=%f,w2=%f,\nR1=%f,R2=%f,\ndt1=%f,dt2=%f,dt3=%f,\nT1=%f,t1=%f,t2=%f,t3=%f,\nx0=%f,x0=%f,x1=%f,x2=%f,x3=%f\n",D1,w,w1,w2,w3,R1,R2,dt1,dt2,dt3,T1,t1,t2,t3,x0,x1,x2,x3);printf("xiamianjisuanmianji:\n");A[0]=m/dt1/;A[1]=n/dt2/;A[2]=h/dt3/;printf("mianjishi:\n");printf("A[0]=%f,\nA[1]=%f,\nA[2]=%f",A[0],A[1],A[2]);printf("\njisuanfeiyong:\n");j1=o*D1*(a0+b0*T1)/10000;for(i=0;i<3;i++){if(A[i]<=100)H[i]=1;elseif(A[i]<=200&&A[i]>100)H[i]=;elseH[i]=;}for(i=0;i<3;i++)printf("%f,%f\n",A[i],H[i]);for(i=0;i<3;i++)M+=(4400+B-620)**+*A[i])*H[i];j2=fc*M/10000;vw=w3*(r3*1000+c0*5)/p/c0/(Tk-5-t0)/3600;j3=c1*c2*o*(21+1450*vw*vw)*vw*p/102/y/10000;j=j1+j2+j3;printf("j=%f,\nj1=%f,\nj2=%f,\nj3=%f\n",j,j1,j2,j3);L=j1/j;E=j2/j;N=j3/j;printf("L=%f,E=%f,N=%f\n",L,E,N);return;}§3·4程序优化计算结果程序运行结果：D1=,w=,w1=,w2=,w3=,A[0]=,A[1]=,A[2]=,Dt1=；R2=；优化后的：总的费用：J=万元,各种费用所占的比例：年生蒸汽费用：J1=万元,L=,蒸发器年折旧维修费用：J2=万元,E=真空系统费用：J3=万元,N=§3·5优化前后费用比较等面积的各项费用：总的费用：J’=万元各种费用所占的比例：年生蒸汽费用：J1’=万元,L’=蒸发器年折旧维修费用：J2’=万元E’=真空系统费用：J3’=万元N=两者比较：总的节省费用：J0=J’-J=万元，节省费用比例:%年生蒸汽费用：J10=J1’-J1=万元,L0=,年折旧维修费用:J20=J2’-J2=万元,E0=,真空系统费用：J30=J3’-J3=万元,N0=,第四章蒸发器工艺尺寸计算蒸发器的主要结构尺寸（以下均以第一效为计算对象）我们选取的中央循环管式蒸发器的计算方法如下。§4·1加热管的选择和管数的初步估计§4·1·1加热管的选择和管数的初步估计蒸发器的加热管通常选用38*无缝钢管。加热管的长度一般为—2m，但也有选用2m以上的管子。管子长度的选择应根据溶液结垢后的难以程度、溶液的起泡性和厂房的高度等因素来考虑，易结垢和易起泡沫溶液的蒸发易选用短管。根据我们的设计任务和溶液性质，我们选用以下的管子。可根据经验我们选取：L=2M，38*可以根据加热管的规格与长度初步估计所需的管子数n’，=124（根）式中S=----蒸发器的传热面积，m2，由前面的工艺计算决定（优化后的面积）；d0----加热管外径，m；L---加热管长度，m；因加热管固定在管板上，考虑管板厚度所占据的传热面积，则计算n’时的管长应用（L—）m.§4·1·2循环管的选择循环管的截面积是根据使循环阻力尽量减小的原则考虑的。我们选用的中央循环管式蒸发器的循环管截面积可取加热管总截面积的40%--100%。加热管的总截面积可按n’计算。循环管内径以D1表示，则所以mm对于加热面积较小的蒸发器，应去较大的百分数。选取管子的直径为：循环管管长与加热管管长相同为2m。按上式计算出的D1后应从管规格表中选取的管径相近的标准管，只要n和n’相差不大。循环管的规格一次确定。循环管的管长与加热管相等，循环管的表面积不计入传热面积中。§4·1·3加热室直径及加热管数目的确定加热室的内径取决于加热管和循环管的规格、数目及在管板撒谎能够的排列方式。加热管在管板上的排列方式有三角形排列、正方形排列、同心圆排列。根据我们的数据表加以比较我们选用三角形排列式。管心距t为相邻两管中心线之间的距离，t一般为加热管外径的—倍，目前在换热器设计中，管心距的数据已经标准化，只要确定管子规格，相应的管心距则是定值。我们选用的设计管心距是：确定加热室内径和加热管数的具体做法是：先计算管束中心线上管数nc，管子安正三角形排列时，nc=*；其中n为总加热管数。初步估计加热室Di=t(nc-1)+2b’,式中b’=(1—d0.然后由容器公称直径系列，试选一个内径作为加热室内径并以该内径和循环管外景作同心圆，在同心圆的环隙中，按加热管的排列方式和管心距作图。所画的管数n必须大于初值n’，若不满足，应另选一设备内径，重新作图，直至合适。由于加热管的外径为38mm，可取管心距为48mm；以三角形排列计算，b’=(1—d0=*d0，Di=t(nc-1)+2b’=48*（13-1）+2*38*=690mm，选取加热室壳体内径为800mm鄙厚为12mm；§4·1·4分离室直径与高度的确定分离室的直径与高度取决于分离室的体积，而分离室的体积又与二次蒸汽的体积流量及蒸发体积强度有关。分离室体积V的计算式为：式中V-----分离室的体积，m3;W-----某效蒸发器的二次蒸汽量，kg/h;P-----某效蒸发器二次蒸汽量，Kg/m3，U-----蒸发体积强度，m3/(m3*s);即每立方米分离室体积每秒产生的二次蒸汽量。一般用允许值为U=~m3/(m3*s)根据由蒸发器工艺计算中得到的各效二次蒸汽量，再从蒸发体积强度U的数值范围内选取一个值，即可由上式算出分离室的体积。一般说来，各效的二次蒸汽量不相同，其密度也不相同，按上式计算得到的分离室体积也不会相同，通常末效体积最大。为方便起见，各效分离室的尺寸可取一致。分离室体积宜取其中较大者。确定了分离室的体积，其高度与直径符合关系，确定高度与直径应考虑一下原则：（1）分离室的高度与直径之比H/D=1~2。对于中央循环管式蒸发器，其分离室一般不能小于，以保证足够的雾沫分离高度。分离室的直径也不能太少，否则二次蒸汽流速过大，导致雾沫夹带现象严重。(2)在条件允许的情况下，分离室的直径尽量与加热室相同，这样可使结构简单制造方便。（3）高度和直径都适于施工现场的安放。现取分离室中U=（m3*s）；m3。H=，，D=§4·2接管尺寸的确定流体进出口的内径按下式计算式中-----流体的体积流量m3/s；U--------流体的适宜流速m/s，估算出内径后，应从管规格表格中选用相近的标准管。§4·2·1溶液进出口于并流加料的三效蒸发，第一效溶液流量最大，若各效设备尺寸一致的话，根据第一效溶液流量确定接管。取流体的流速为m/s；所以取ф规格管。§4·2·2加热蒸气进口与二次蒸汽出口各效结构尺寸一致二次蒸汽体积流量应取各效中较大者。所以取ф规格管。§4·2·3冷凝水出口冷凝水的排出一般属于液体自然流动，接管直径应由各效加热蒸气消耗量较大者确定。所以取ф规格管。 第五章、蒸发装置的辅助设备§5·1气液分离器蒸发操作时，二次蒸汽中夹带大量的液体，虽在分离室得到初步的分离，但是为了防止损失有用的产品或防止污染冷凝液，还需设置气液分离器，以使雾沫中的液体聚集并与二次蒸汽分离，故气液分离器或除沫器。其类型很多，我们选择惯性式除沫器，起工作原理是利用带有液滴的二次蒸汽在突然改变运动方向时，液滴因惯性作用而与蒸汽分离。在惯性式分离器的主要尺寸可按下列关系确定：D0=D1；D1：D2：D3=1：：2H=D3h=~D0--------二次蒸汽的管径，mD1--------除沫器内管的直径，mD2--------除沫器外管的直径，mD3--------除沫器外壳的直径，mH---------除沫器的总高度，mh---------除沫器内管顶部与器顶的距离，mD1=103mmD2=D3=206mmH=206mmh=选取二次蒸汽流出管：除雾器内管：除雾器外罩管：§5·2蒸汽冷凝器蒸汽冷凝器的作用是用冷却水将二次蒸汽冷凝。当二次蒸汽为有价值的产品需要回收或会严重地污染冷却水时，应采用间壁式冷却器。当二次蒸汽为水蒸气不需要回收时，可采用直接接触式冷凝器。二次蒸汽与冷凝水直接接触进行热交换，其冷凝效果好，被广乏采用。现采用多孔板式蒸汽冷凝器： §5·2·1冷却水量根据冷凝器入口蒸汽压强和冷却水进口温度可由图表可查得。-----冷却水量m3/h；-----所需冷凝的蒸汽量，Kg/h实际取§5·2·2由计算可知,进入冷凝器的二次蒸汽的体积流量可计算得到冷凝器的直径D取D=300mm§5·2·3淋水板的设计(1)淋水板数：D<500mm,淋水板取4块(2)淋水板间距：L末又L0=D+~m取L0=+=。L1=*L0=*=L2=*L1=L3=*L2=<所以取L0=。L1=L2=L3=（3）弓型淋水板的宽度：最上面一块=（~）D，其他各淋水板B=D/2+m所以=*~*，取=，其余取B=2++。(4)淋水板堰高h：D=300mm<500mm,取h=40mm(5)淋水孔径：冷却水循环使用,取d=10mm(6)淋水孔数：淋水孔冷却水流速,,考虑到长期操作易堵,则:N1=n*（1+12%）=134最上一板孔数为N1=n*（1+12%）=134取整为134个其他各板孔数为取整为125个§5·3泵的选型：(1)   真空泵真空度为,真空泵排气量---真空系统渗透空气量;---蒸发过程中料液释放的不凝性气体量,常可以忽略;---每小时冷却水能释放的空气量;---蒸发过程中流体的饱和蒸汽压的当量值;---不凝性气体机械夹带的生蒸汽量,常忽略;=h；=0==++)=+0+=h=0=++++=.根据排气量和所抽气体基本上是空气的前提条件,选用水喷射泵11/2BA-6系列功率极限真空度材料最大吸气量600mmHgA3钢h(2)离心泵q=F/ρ==h查IS型单级单吸离心泵性能表：选型号：IS50-32-125流量：h扬程20m转数：2900r/min轴功率：电机功率：效率60%必须汽蚀余量：§5·4预热器的选型：将h的10%NaOH溶液从20%加热到。加热介质采用500Kpa绝压的饱和水蒸气，冷凝液在饱和温度下排出：（1）      确定流体通入的空间：因加热介质是饱和蒸汽宜通入壳程NaOH溶液通入管内。（2）      确定流体的定位温度物性数据并选择列管式换热器的型式，饱和蒸汽压强为500Pa,绝压温度T=ºC。因是恒温，饱和蒸汽的定性温度下Tm=ºC，NaOH的定性温度tw=(20+=ºC.。两流体的温度差Tm-tm=两流体的温差大于50ºC，故选用浮式列管换热器。两流体在定性温度下的物性数据如下：物性/流体温度/ºC密度/Kg/m3μ/Cp/KJ/Kg.ºCλ/W/mºCΓKj/KgNaOH溶液 冷凝水     （3）      计算热负荷Q，按管内NaOH溶液计：即Q=Wc*Cp*（t2-t1）=3600**=607900w若忽略热损失，冷凝水又在饱和温度下排出，蒸汽的流量，可由热量恒算求得：Wh=Q/r=607900/2113200=s(4)计算平均温度差：∆tm=[(T-t1)-(T-t2)]/ln[(T-t1)/(T-t2)]=43ºC(5)选换热器规格：根据管外为水蒸汽冷凝来加热管内的NaOH溶液（μ<），总传热系数经验的范围为1160~4070。选Ko=1400W/（m2.ºC）So=Q/（Ko*∆tm）=10。9m2所以选FA系列浮头式换热器规格如下：壳径D---------325mm公称面积S---------10m2实际面积管长L------3m管程数N-----2管数n------〉6管子排列方法--------正三角形排列管子直径ф19*2mm第六章主要设备强度计算及校核筒体厚度计算公式设计厚度公式中,:设计压力;:焊缝系数::筒体内径;材料的许用应力,MPa压力实验选用水压试验,公式为---容器的有效厚度;---圆筒焊缝系数---试验时容器承受压力;---设计温度下屈服极限§6·1蒸发分离室厚度设计分离室壳体材料选用A3钢,查表得设计温度下,P==最小厚度查取:<3向上整和为,水压试验所以强度符合按标准管当Ri=1200mm时取10mm符合要求封头选用椭圆型封头,则厚度为焊接方便取封头厚度和筒体厚度相等即:10mm§6·2加热室厚度校核取进行校核PT=即所以水压实验强度符合要求查标准管，内径为900mm,最小壁厚为10mm,强度符合要求。  第七章课程设计心得时光荏苒，白驹过隙。转眼间，为期三周的化工原理课程设计就已经接近尾声了。回首望去，心情格外的开阔，感慨颇多。我忘不了和杨老师以及同学们一起度过的日日夜夜，忘不了我们组的几个人因为一个数据的来源而吵的面红耳赤，更忘不了看到一个个成果的喜悦………首先，要感谢杨老师能给我们提供一个进行实践锻炼的舞台。以前我们学过的知识只不过是纸上谈兵，而化工原理课程设计却是以门综合性课程，它不仅要求我们对化工设计有基础的了解，而且还要对化工原理、化工机械基础、化工热力学等一系列知识能够进行综合的运用，同时也对计算机软件和编程提出了较高的要求。，这次课程设计，我们小组共有5个人。作为组长，我首先对他们进行了分工：三个人负责计算，一各人负责编程，以各人负责画图。设计过程中，我们相互鼓励，相互帮助，发挥团队协作精神，齐心协力攻克了一个又一个难关，我主要负责画图。开始，我一边自学AutoCAD，一边进行计算。一段时间后，我开始着手进行流程图的画图。由于是自学，所以边学边画，速度比较慢，错误也比较多。但是在自己的摸索和同学们的相互帮助下，我逐渐掌握了AutoCAD的各种基本操作。特别感谢杨老师，在我画图的过程中给与了大量的指导和鼓励，使我的图更加完善。现在看到自己设计的流程图和设备图，心中无比充实。想起自己为画图和设计而熬过的几个不眠之夜，心中由衷地感激值得！这次设计，我不仅巩固了化工原理及相关知识，而且增强了团队的协作精神，同时也磨炼了意志。相信这次课程设计会让我们更加注意理论与实践的结合，成为我们人生中一笔宝贵的财富！参考文献：1、《化工原理课程设计》天津大学天津科学技术出版社柴城敬刘国维李阿娜编P7~~P74P1862、《化工原理》　　合肥工业大学　　　　　　　合肥工业大学出版社崔鹏　　　魏凤玉　编　P114~~P190P392~P395P412　3、《化工工艺设计手册》上册第一版（修订）国家医药管理局上海医药设计院编化学工业出版社P385~P391P671,P687,P654,P633,P667,4、《化学工程手册》第二版，上册第九篇蒸发，9-1~~9-20化学工业出版社，时均，汪家鼎，余国琮，陈敏恒主编5、《化工原理优化设计与解题指南》阮奇6、梁虎，王黎，朱平，多效蒸发系统优化设计研究。化学工程，1997，25（6）；48~517、华南工学院化工原理教研组，化工过程设备设计，广州；华南工学院出版社；1987。86~928、韦鹤平，最优化技术[m]，上海；同济大学出版社，1987，119~1089、阮奇，黄诗煌，叶长桑，陈文波，多效蒸发系统油画设计目标函数的建立与求解，计算机与应用化学，200110、时钧，汪家鼎，余国综，陈敏恒，化学工程，第二版，上卷，北京，化学工业出版社，1996，1`10811、张瑞生，沈才大，化学系统工程基础，上海，华东化工学院出版社，`23512、谭学富，李茂林，王红心等，氢氧化钠水溶液多效蒸发调优，沈阳化工学院学报，（1），25`30