1、内蒙古工业大学本科课程设计说明书第1章 蒸发方案设计1设计方案的确定 一般的加热蒸汽压强在0.50.8MPa范围内,此次选600 KPa。 加热蒸汽的确定需要考虑加热蒸汽温度的上限和下限。被蒸发的溶液有一个最高的蒸发温度,超过此温度蒸发就物料就会变质,破坏和分解,这是确定加热蒸气压强的一个依据通常所用饱和蒸汽的温度不超过180,超过时相应的压强就很高,这将增加加热的设备费和操作费。冷凝器的操作压强要80 KPa,(真空度),由于呼市当地大气压为100 KPa,,所以本次选取冷凝器的操作压强为20KPa 进入冷凝器冷凝需消耗大量冷却水,而且溶液粘度大,传热差。但对于那些热敏性物料的蒸发,为充分利
2、用热源还是经常采用。对混合式冷凝器,其最大的真空度取决于冷凝器内的水温和真空装置的性能。若第一效用较高压强的加热蒸汽,则末效可采用常压或真空蒸发,此时末效产生的二次蒸汽具有较高的温度,可以全部利用。而且各效操作温度高时,溶液黏度低,传热好。若一效加热蒸汽压强低,末效应采用真空操作。此时各效二次蒸汽温度低,所以此次选取冷凝器的操作压强为20KPa。蒸发的类型:本设计采用中央循环管式蒸发器 中央循环管式蒸发器结构紧凑,制造方便,操作可靠,故在工业上应用广泛,有所谓标准蒸发器。但设备的清洗和检修能保证较大的传热系数,满足生产工艺的要求;生产能力大,能完善分离液沫,尽量减慢传热面上垢层的生成;蒸发效数
3、的确定:本设计蒸发器效数采用3效.由于为充分利用热能,为提高热能利用效率,生产中一般采用多效蒸发,但并不是效数越多越好,效数的选择受经济上和技术上的限制。经济上的限制是指效数超过一定数时经济上不合算。多效蒸发中,随效数的增加,总蒸发量相同时所需蒸汽量减少,使蒸汽用量减少,使操作费用降低。但随效数增加,设备费成倍增长,而所节省的蒸汽量愈来愈少,所以无限制增加效数已无实际意义,最适宜的效数应使设备费和操作费二者之和为最小。技术上的限制是指效数过多,蒸发操作难于进行。一般工业秤中加热蒸汽压强和冷凝器操作压强都有一定限制,因此在一定操作条件下,蒸发器的理论总温度差为一定值。在效数增加时,由于各效温差损
4、失之和的增加,使总有效温差减小,分配到各效的有效温差小到无法保证各效发生正常的沸腾状态时,蒸发操将无法进行下去。综合考虑经济上和技术上的限制,选择三效。蒸发流程的确定:本设计蒸发流程采用并流模式。并流法变称顺流法,其料液和蒸汽呈并流。并流蒸发的优点有:(1)溶液从压力和温度较高的蒸发器流向压力和温度较低的蒸发器,料液能自动从前效进入后效,可省去输料泵;(2)前效的温度高于后效,料液从前效进入后效时呈过渡热状态,可以产生自蒸发;(3)操作简便,工艺条件稳定。缺点:随着溶液从前一效遂效流向后面各效,其组成增高,而温度反而降低,致使溶液的黏度增加,蒸发器的传热系数下降,所以对于随组成增加粘度变化很大
5、的料液不宜采用并流。逆流法溶液溶液的流向与蒸汽的流向相反,即加热蒸汽从第一效加入,料液从最后一效加入。逆流法的优点是:随溶液的组成沿着流动方向的增高,其温度也随之升高。因此因组成增高使黏度增大的影响大致与因温度升高使黏度降低的影响相抵,故各效溶液的黏度较为接近,各效的传热系数也大致相同。进料温度的确定:采用沸点进料。采用沸点进料,操作简单,计算简便。这样也会降低操作成本。多效蒸发的操作流程根据加热蒸汽与料液的流向不同,可分为并流、逆流、平流及错流四种。并流法变称顺流法,其料液和蒸汽呈并流。并流蒸发的优点有:(1)溶液从压力和温度较高的蒸发器流向压力和温度较低的蒸发器,料液能自动从前效进入后效,
6、可省去输料泵;(2)前效的温度高于后效,料液从前效进入后效时呈过渡热状态,可以产生自蒸发; 第2章 蒸发过程工艺计算计算的基本依据是物料衡算,热量衡算,传热速率方程三个基本关系。计算常采用逐次逼近的方法,其计算步骤如下:(1)假设冷凝水在饱和温度下排出,无额外蒸汽引出。(2)根据生产经验算据,初步估计各效蒸发量和各效完成液的浓度;(3)根据经验假设蒸汽通过各效的压强降相等,估算各效溶液沸点和有效总温差;(4)根据蒸发器的焓衡算,求各效的蒸发量和传热量;(5)根据传热速率方程计算各效的传热面积。若得的各效传热面积不相等则应按下面介绍的方法重新分配有效温度差,重复(3)(5),直到所求得的各效传热
7、面积满足要求为止。2.1 各效蒸发量和完成液组成的估算本设计的操作条件是:(1)NaOH水溶液处理量50kt/a,蒸发前NaOH溶液浓度为10%,温度30。完成液含NaOH浓度为32%;用三效并流蒸发装置。(2)加热蒸汽压强为600,冷凝器的压强为20。(3)各效蒸发器的总传热系数:=1350W/(),=1000 W/(),=400W/()。2.1.1 各效蒸发量和完成液浓度的估算:总蒸发量: (2-1)W:蒸发出的水量,F原料液流量,原料液初始浓度,完成液浓度。采用三效蒸发:n3 假设各效蒸发量相等,即:=1:1:1 4300/31430kg/h 因而初估各效完成液的浓度为: (2-2) 第
8、I效完成液的浓度:13% ; 第II效完成液的浓度:%第III效完成液的浓度:32%。2.1.2 估算各效二次蒸汽温度为求各效溶液沸点,需假定压强。一般加热蒸汽压强和冷凝器操作压强是已知的,其他各效二次蒸汽的压强可按各效间蒸汽压强降相等的假设来确定,即各效加热蒸汽压强与二次蒸汽压强之差为: (2-3)故第效二次蒸汽压强: (2-4)193;600193407;6002193214 ;600319021。将各效二次蒸汽的压力、附录中查到相应的二次蒸汽的温度和二次蒸汽的汽化潜热列于下表中:各效二次蒸汽的压力、二次蒸汽的温度和二次蒸汽的汽化潜热如下表效次二次蒸汽的压力, 40721421二次蒸汽的温
9、度,(即下一效蒸汽的温度)144122.1661.2二次蒸汽的汽化潜热,/(即下一效加热蒸汽的汽化潜热)2136.72199.22859.22.1.3 计算各效传热温度差各效传热温度差计算式为:=,=+ (2-5)式中为前一效二次蒸汽温度(即第效加热蒸汽温度),为第效溶液沸点,为第效二次蒸汽温度,为第效温差损失。各效温差损失: (2-6)(1)由于溶液中含有不挥发溶质所引起的沸点升高:=-式中:为与溶液压力相等时水的沸点,为第效溶液的沸点。(2)由于液柱静压头引起的沸点升高某些蒸发器操作时,器内需维持一定液位,因而溶液内部压强大于液面压强,致使实际沸点较液面为高,两者之差即为因溶液静压强一起的
10、温差损失。为简便计算溶液内部沸点升高按液面与底层的平均压强下水的沸点和二次蒸汽压强下水的沸点差估算。平均压强按静力学方程式计算:=+ ,则: = (2-7)式中:为蒸发器中液面与底层的平均压强,; 为对应下水的沸点,; 为二次蒸汽(即液面处)的压强,; 为对应p下水的沸点,; 为溶液的平均密度,; L为液层高度,。(3)由于管道流动阻力产生的沸点升高在多数蒸发中,各效二次蒸汽流到下一效加热室时,由于管道阻力使其压强降低,致使蒸汽的饱和温度相应降低,由此引起的温差损失即为。根据经验,一般取=1。 2.2 计算过程 2.2.1 各效二次蒸汽温度查NaOH水溶液浓度和密度表得各效蒸发器中溶液的平均密
11、度见下表:表 2-2 各效二次蒸发器中溶液的平均密度效次溶液的浓度13%18%32%溶液的密度103kg/m31.0141.0601.239对于NaOH水溶液可采用杜林经验公式计算:(1)+m0.142+150.752.71 (2-8)各效完成液浓度。所以得出:(k1)+m0.142+150.752.71=0.14213%144+150.7513%13%2.7113%4.85; 同理:= (k1)+m0.142+150.752.71= 0.142 18% 122.16 +150.75 18% 18% 2.7118% 7.5685 , (k1)+m0.142+150.752.71 = 0.142
12、 32% 122.16 +150.75 32% 32% 2.7132% =17.3236 。计算由于蒸发器中溶液静压强引起的温度差损失, (2-9)查表得出在各个状态下溶液的平均密度为,。得出:;。对应的,所以,。故;。2.2.2 计算传热温度差 各效传热温度差计算式为, (2-10)式中:第效二次蒸汽的温度;第效溶液的沸点。,;,;,将各效溶液的沸点列于下表中: 表2-3 各效溶液沸点效次溶液(温度)沸点,150.85132.7385.332.2.3 求各效蒸汽用量及各效蒸发量考虑溶液的稀释热,引入热利用系数 (2-11)第效:因沸点进料,故 =热利用系数0.980.7=0.980.7(13
13、%10%)=0.959查水蒸气表得:压力为600KPa的加热蒸汽的汽化热=2091.1,=150.85的二次蒸汽的汽化热=2136.7。则 = = =0.94 第效:引入热利用系数:热利用系数=0.980.7=0.980.7(18%13%)=0.945第效中溶液沸点=132.73的二次蒸汽的汽化热=2199.22,无额外蒸汽引出,=2136.7,已知/(kg),/(kg)。则:( =0.945(62503.7-4.187) =0.937+183 第效:热利用系数=0.980.7=0.980.7(32%18%)=0.882第效中溶液沸点=85.33的二次蒸汽的汽化热=2859.2,无额外蒸汽引出
14、,=2199.22,则()=0.882(62503.74.1874.187)=0.61690.06+338 +=4300,联立得出:=1553.8,=1460.6,=1551.58,=1207.535。则各效传热量分别为:=1553.82091.1=3.249=2136.71460.6=3.12=1551.582199.22=3.4122.2.4计算传热面积=()/ = 76.942.86/76.9=0.44,误差过大,所以调整各效的有效温度差,重复上述计算过程。2.3 重新计算过程2.3.1 有效温度差的重新分配=76.6457.8554.6=11=重复上述计算过程: 各效蒸发量: =146
15、0.6 =1551.6 =1207.5由于要完成指定的任务所以不变 ,所以有以下的计算过程:=32%2.3.2 重新求算各效溶液沸点由于第三效溶液的沸点:冷凝器的压力和完成液的组成均不变所以有,=17.3236+5.81+1=24.13=+=85.33+29=114.33第二效溶液的沸点:=114.33 查表可知=164 =2221.6= (k1) + mK=1+0.1420.18=1.02556M=150.750.180.182.710.18=4.3965=0.02556114.3+4.3965=7.32=1641.0609.812.52 = 177=116.06 =114.3=116.06
16、114.3=1.76=1=114.3+10.08=124.38第二效加热蒸汽的温度: =由=140.38 查表 (k1)+m=(0.1420.12)140.38150.750.120.122.710.12=4.2=374查表 =141141140.38=0.628=1则第一效液相沸点: +5.868=146.248核算热利用系数0.980.7=0.980.7(12%10%)=0.966 = = 0.9662091.12148.7 =0.94热利用系数=0.980.7=0.980.7(18%13%)=0.945( = 0.938(62503.7-4.187) =0.868+214 热利用系数=0
17、.980.7=0.980.7(32%18%)=0.952()=0.952(62503.74.1874.187)=0.6850.057+301 +=4300解得=1512.48 =1098.8 =967.46 =1609则各效传热量分别为:=2091.11609=3.36=1512.482148.7=3.25=2221.61098.8=2.42.3.4 计算各效传热面积=各项误差均小于,所以计算符合要求。=计算结果与初设值基本一致 传热面积加10%余量 S=56.47(1+0.1)=62.117将所有计算结果列入下表中:表2-6 蒸发器有关计算数据汇总效次冷凝器加热蒸汽温度,158.7138.6
18、261.2操作压力,60020溶液温度,146.24124.3885.33完成液浓度12%18%32%蒸发量,kg/h1512.481098.8967.46蒸汽消耗量,kg/h1609传热面积,m255.556.4257.5第3章 蒸发器工艺尺寸的设计中央循环管式蒸发器式蒸发器主体为加热室和分离室,加热室由直立的加热管束组成,管束中间为一根直径较大的中央循环管;分离室是汽液分离的空间,其主要结构尺寸包括:加热室和分离室的直径和高度,加热管和循环管的规格,长度及在花板上的排列方式等。这些尺寸取决于工艺计算结果,主要是传热面积。3.1 加热管的选择和管束的初步估计3.1.1 加热管的选择和管数的初
19、步估计加热管子长度的选择应根据溶液结垢的难易程度,溶液的起泡性和厂房的高度等因素来考虑。易结垢和起泡沫溶液的蒸发宜采用短管。本设计加热管选用382.5mm,长为3m。当加热管的规格与长度确定后,可由下式初估所需管数n: (3-1)式中: S为蒸发器的传热面积,m2 ,由前面的工艺计算决定; d0为加热管的外径,m; L为加热管长度,m。因加热管固定在管板上,考虑到管板厚度所占据的传热面积,计算n时管长用(L-0.1)m。为完成传热任务所需的最小实际管数n只有在管板上排列加热管后才确定。n=1833.1.2 循环管的选择 循环管的截面积是根据使循环阻力尽量减少的原则来考虑的。以截面积的60% 计
20、算,若以D1表示循环管内径,则: =(0.41) (3-2)因而: 对面积较小的蒸发器,应取较大的百分数。按上式计算出D1后,应从管规格表中选取管径相近的标准管,只要n与n相差不大,循环管的规格可一次确定。循环管的管长与加热管相等,循环管的表面积不计入传热面积。循环管选取3516mm的无缝不锈钢管。3.1.3 加热室直径及加热管数目的确定加热室的内径取决于加热管和循环管的规格,数目及在管板上的排列方式有正三角形,正方形,同心圆等,目前以三角形排列居多。管心距t为相邻两管中心线之间的距离,t一般为加热管外径的1.25-1.5倍。目前在换热器设计中,管心距的数值已经标准化,只要管子规格确定,相应的
21、管心距则为确定值。管心距值见表3-1。表3-1 不同加热管尺寸的管心距加热管外径d0,mm19253857管心距t,mm25324870加热室内径和加热管数采用做图法来确定,具体做法是;先计算管束中心线上管束nc:本设计采用正三角形排列加热管,即nc=1.1=1.1=12.7=138=38mm=t(-1)+2=48(131)+238=652mm根据初估加热室内径值和容器公称直径系列,试选一内径作为加热室内径,并以此内径和循环管外径作同心圆,在同心圆的环隙中,按加热管的排列方式和管心距作图。 由作图可知n=220n 壳体内径的标准尺寸列表3-2 壳体的尺寸标准壳体内径mm400-700800-1
22、0001100-15001600-2000最小壁厚mm8101214根据表3-2列表,本设计选取加热室壳体内径为900 mm,壁厚为10 mm。3.1.4 分离室直径与高度的确定分离室的直径和高度取决于分离室的体积,而分离室的体积又与二次蒸汽的体积流量及蒸发体积的强度有关。分离室体积V的计算公式为:m3 (3-3)式中:W某效蒸发的二次蒸汽量,kg/h; 某二次蒸发的二次蒸汽密度,kg/m3; U蒸发体积强度,m3/(m3s),一般允许值为1.1-1.5m/s。根据由蒸发工艺计算得到的各效二次蒸汽量,再选取适当的U值,即可得V。但各效二次蒸汽量,密度不同,按上式计算得到V值也不同,通常末效最大
23、。为方便计,各效分离室可取一致,分离室估积宜取其中较大者。各效蒸汽的密度见表3-3。表3-3 各效蒸汽的密度值效次第I效第II效第III效二次蒸汽温度,140.38114.361.2密度,kg/m31.9620.960.137将各值代入上述公式得: U取1.5 m/s。=1.3m3。分离室体积取其中较大者,V=1.3 m3分离室体积确定后,其高度H与直径D符合下列关系:V= (3-4)确定高度与直径时应考虑以下原则: (1) 分离室的高度与直径之比H:D =12。对中央循环管式蒸发器,其分离室高一般不小于1.8m,保证足够的雾沫分离高度。分离室的直径也不能太小,否则二次蒸汽流速过大将导致雾沫夹
24、带现象严重;(2) 在条件允许时,分离室直径应尽量与加热室相同,这样可使加热室结构简单,制造方便;(3) 高度和直径都适于施工现场的安放。本设计选取 H=1.8 分离室直径D=0.959mH:D =1.876 介于12之间而且大于加热室直径 满足要求。3.2接管尺寸的确定3.2.1溶液的进出口 对于并流加料的三效蒸发,第一效溶液的流量最大,若各效设备采用统一尺寸,应根据第一效溶液流量来确定接管。管径的计算公式为: d= (3-5)式中:为第一效溶液的体积流量; u为溶液的流速(u=13m/s)考虑。设计时进出口直径可取为一致。取u=2m/s,=926.1kg/m3 则d=34.6mm,查表得出
25、蒸发器溶液进出口管选取423mm的不锈钢无缝钢管。3.2.2加热蒸汽进口与二次蒸汽出口加热蒸汽进口: 若各效结构尺寸一致,则二次蒸汽体积流量应取各效中较大者,饱和蒸汽适宜流速u=70120m/s。计算公式同溶液进出口计算公式一样,加热蒸汽进口,取二次蒸汽流速为u=80m/s:=508/h=1575.5/h =/s加热蒸汽体积流量取其中较大者,V=8020.4/s则:d=188.35mm,加热蒸汽进口管取2036mm的不锈钢无缝钢管。二次蒸汽出口:=/s =0.316/s=1.96/s二次蒸汽出口:二次蒸汽体积流量取其中较大者,V=1.96/s,则:d=176.6mm二次蒸汽出口同样选取1945
26、mm不锈钢无缝钢管。3.2.3冷凝水出口由于冷凝水自然流出,所以冷凝水的排出属于自然流动, u=0.81.8m/s,接管直径由各效加热蒸汽消耗量较大者确定。取u=0.1m/s则: d=20.5mm冷凝水出口管取325mm的不锈钢无缝钢管。核算流速加热蒸汽进口:当d=42-6=36mm时 d= u=46250106/3600第4章 蒸发装置辅助设备的设计蒸发装置的辅助设备主要包括气液分离器,蒸汽冷凝器等辅助设备。4.2气液分离器气液分离器选惯性式除沫器。蒸汽操作时,二次蒸汽中夹带大量的液体,虽在分离室得到初步分离,但是为了防止有用的产品损失或污染冷凝液,还需设置气液分离器,以使雾沫中的液体凝聚并
27、与二次蒸汽分离,故气液分离器又称捕沫器或除沫器。本设计使用惯性式除沫器。惯性式除沫器的主要尺寸按下列关系确定:,=1:1.5:2H=,h =(0.40.5)式中:二次蒸汽的管径,m; 除沫器内管的直径,m; 除沫器外罩管的直径,m; -除沫器的外壳直径,m; H除沫器的总高度,m; h除沫器内管顶部与器顶的距离,m。因为=194mm=0.194m,=0.194m,=0.291m,=0.388m,H=0.338m,h=0.097m。最终选定标准为=21912mm不锈钢无缝钢管,=mm不锈钢无缝钢管,=3776.5mm不锈钢无缝钢管。 42蒸汽冷凝器 421冷凝器主要类型 蒸汽冷凝的作用是用冷凝水
28、将二次蒸汽冷凝。冷凝器分为直接冷凝器和 间接冷凝器。间接冷凝器价格昂贵,用水量大。当二次蒸汽为水蒸汽不需要加收时,可采用直接接触式冷凝器。直接接触式冷凝器的冷凝效果好,结构简单,操作方便,价格低廉,因此在本设计中选用直接接触式冷凝器。直接接触式冷凝器有多孔板式,水帘式,填充塔式及水喷射式等。水喷射式冷凝器的工作原理是冷凝水依靠泵加压后经喷嘴雾化使二次蒸汽冷凝,不凝汽也随冷凝水由排出管排出。此过程产生真空,则不需要真空泵就可造成和保持系统的真空度,但单位二次蒸汽所需的冷凝水量大,二次蒸汽量过大时不宜采用。各种形式蒸汽冷凝器的性能见表4-1。表4-1 各种形式冷凝器的性能表项目形式多层多孔板车单层
29、多孔板式水帘式填充塔式水喷射式水气接触面积大较小较大大最大压强降,Pa1067-2000小1333-3333较小大塔径范围大小均可不宜过大350mm100mm二次蒸汽量2t/h结构与要求较简单简单较简单,安装有一定要求简单不简单,加工有一定要求水量较大较大较大较大最大孔易堵塞适于腐蚀性蒸汽422冷凝器的设计与选用4221工作水量的计算 计算式 (4-1)式中: W冷凝水用量,kg/h; D冷凝蒸汽量,kg/h; I蒸汽的焓,kj/kg; ,冷却水进出口温度,; 冷却水平均比热,kj/(kg)。查得:61.2时蒸汽的焓值I=2606kj/kg,=4.187 kj/(kg),=12,=25,D=9
30、67.46 kg/h。=4.4kg/h。4.2.2.2 喷射器结构尺寸计算(1)喷嘴数n及喷嘴直径通过一个喷嘴的水的流速为:= (4-2)式中:工作水通过喷嘴的压强降,Pa ; 水的密度,kg/; 流量系数,可取0.930.96。喷嘴直径在水质清洁时可取58mm,一般为1222mm即可,选定直径后,喷嘴个数的确定: (4-3) 取=1000 kg/,=0.93,=(300-21)kPa, W=kg/h,=12mm=0.93=21.97m/s。,n=5 (2) 文氏管喉部直径kPa (3) 水喷射气其他各部尺寸的选取范围文氏管喉管长度:=(24)=322.6=67.8文氏管收缩口直径:=(1.5
31、81.78)=1.6022.6=36.16文氏管收缩段长度:=(1.75.15)()=4(36.1622.6)=54.24文氏管收缩角度:=左右文氏管扩散段直径:=1.78=1.7822.6=40.228文氏管扩散段长度:=(415)()=(40.22822.6)10=176.28文氏管扩散段角度:=3左右 取进水管直径与进蒸汽管直径可按一般原则计算,但管内水流速不宜太高,否则电耗增大,当真空度在7395kPa时,其最大流速为70120不等。所以最大流速为80m/s。详见后面表格。(4) 射流长度的决定喷射水的射流长度,是指喷射嘴出口处到聚焦点的水柱长度。该长度愈长,汽室就要愈高。喷射冷凝器的
32、射流长度与喷嘴直径关系大致如下:表4-2 喷射冷凝器的射流长度与喷嘴直径关系喷嘴直径,mm喷射水流长度,mm16201215581500200012001500200500所以喷嘴直径为12mm,喷射水流长度为1300mm。表5-3 除沫器尺寸项目尺寸文氏管喉管长度67.8mm文氏管收缩口直径36.16mm文氏管收缩段长度54mm文氏管收缩角度11文氏管扩散段直径40mm文氏管扩散段长度176mm文氏管扩散段角度36章 对本设计的评述通过此次对蒸发装置的设计我收获很多。首先是无论做什么事情都要动脑筋,就比如说用CAD画加热管的排列数目时我无从下手,同学给指点一下豁然开朗,其实很简单的事情只是不
33、去动脑,这样会浪费很多的时间,会影响设计进度。其次就是以后做事一定要提前做好准备,可以说是成败的关键,课设的第一天就遇到困难了,一些物性参数无法查阅,当时既着急又后悔,所以这也很关键。 再次就是做事一定要认真,我是一个很马虎的人,导致在大量的数据中弄得很乱,不得不重新进行计算,既伤神又费时,很不值得。此次设计的最大收获是学会了WPS的使用,编写一些公式非常方便,有利于节省时间。此次设计让我基本对蒸发装置有了很深的认识,丰富了我的视野,为以后工作奠定了一定的基础,个人认为这门课程非常有意义,最后要特别感谢这些天为我们服务的指导老师们,真的辛苦了,你们的辛苦换来了我如此大的收获,真心谢谢各位老师。如果此次设计中存在错误,请老师多多指点。参考文献1.柴诚敬,张国亮.化工流体流动与传热.北京:化学工业出版社.2007:298327(P337-339)2.朱孝钦.过程装备基础.化学工业出版社2008:(P200-216) 3.贾绍义,柴诚敬.化工原理课程设计.天津:天津大学出版社.2002:7310025