碳酸丙烯酯(PC)脱碳填料塔的设计

《碳酸丙烯酯(PC)脱碳填料塔的设计》由会员分享，可在线阅读，更多相关《碳酸丙烯酯(PC)脱碳填料塔的设计（28页珍藏版）》请在装配图网上搜索。

1、碳酸丙烯酯（PC）脱碳填料塔设计目 录化工原理设计任务书 3一、设计目的 3二、设计任务 3三、设计条件 3四、基础数据 4五、设计内容5一、计算前的准备61.CO2在PC中的溶解度关系62.PC密度与温度的关系73.PC蒸汽压的影响84.PC的粘度8二、物料衡算81.各组分在PC中的溶解量82.雾沫夹带量Nm3/m3PC93.溶液带出的气量Nm3/m3PC94.出脱碳塔净化气量105.计算PC循环量106.验算吸收液中CO2残量为0.15 Nm3/m3PC时净化气中CO2的含量107.出塔气体的组成11三、热量衡算121.混合气体的定压比热容122.液体的比热容133.CO2的溶解热144.

2、出塔溶液的温度145.最终的衡算结果汇总15四、设备的工艺与结构尺寸的设计计算161.确定塔径及相关参数15五、填料层高度的计算18六、填料层的压降26七、附属设备及主要附件的选型25八，塔总高计算 28九、设计概要表 29十、对本设计的评价与总结 28参考文献28+化工原理课程设计任务书一、化工原理课程设计目的、任务1. 培养学生查阅资料选用公式和搜索数据的能力2. 培养学生在填料吸收塔、精馏塔设计时，既考虑技术上的先进性和可行性，又考虑经济上的合理性并注意操作时的劳动条件和环境保护的正确设计思想。3. 培养学生能迅速准确的对填料塔进行工艺设计计算的能力4. 培养学生能用简洁的文字清晰的图表

3、来表达自己设计思想的能力二、设计任务碳酸丙烯脂（PC）脱出CO2气体填料吸收塔设计三、设计条件1、混合气（变换气）处理量： 40000 Nm3/h2、进塔混合气体成分： 原始数据表（均为体积%，下同)3、进塔吸收剂（碳酸丙烯酯PC）入塔浓度，自定；4、气液两相的入塔均选定为：305、出塔净化气中CO2浓度0.6%6、操作压力：1.6MPa进塔变换气CO2COH2N2合计体积百分数,%29.01.547.222.3100组分分压，MPa0.4640.0240.7550.3571.600组分分压，kgf/cm24.7350.2457.7013.63816.319原始数据表四、基础数据1.碳酸丙烯酯

4、（PC）的物理性质正常沸点，（）蒸汽压133.32-1Pa粘度，mPas分子量20430382050102.090.10.242.761.62温度，（）015254055（kg/m3）122412071198118411692.比热计算式 3.CO2在碳酸丙烯酯（PC）中的溶解度温度t，（）2526.737.84050亨利系数E101.3-1kPa81.1381.7101.7103.5120.84.CO2在碳酸丙烯酯（PC）中的溶解热可近似按下式计算（以表示）5.其他物性数据可查化工原理附录五、化工原理课程设计主要内容 1、工艺及设备设计（1）设计方案和工艺流程的说明（2）填料吸收塔的工艺计算

5、；吸收剂用量求取：最小吸收剂用量，吸收剂用量；操作线方程；填料塔径求取：选择填料，液泛速度，空塔气速，塔径及圆整，最小润湿速度求取及润湿速度的选取，塔径的校正；传质单元高度的求取；传质单元数的求取；填料层高度；单位填料层压降的求取；吸收塔高度计算；液体分布；再分布及分布器的选型；填料吸收塔的工艺流程图；（3）填料吸收塔设备设计 填料吸收塔附属结构的选型与设计； 全塔高度：包括上、下封头，裙座高度。 2、制图 包括工艺流程图、设备图。 3、编写设计说明书一：计算前的准备：1.CO2在PC中不同温度下亨利系数数据温度t，（）2526.737.84050亨利系数E101.3-1kPa81.1381.

6、7101.7103.5120.8作图得：亨利系数与温度近似成直线，且kPa因为高浓度气体吸收，故吸收塔内CO2的溶解热不能被忽略。现假设出塔气体的温度为，出塔液体的温度为，并取吸收饱和度（定义为出塔溶液浓度对其平衡浓度的百分数）为70%，然后利用物料衡算结合热量衡算验证上述温度假设的正确性在40下，CO2在PC中的亨利系数E40=103.5101.3 kPa=10485 kPa1出塔溶液中CO2的浓度（假设其满足亨利定律）（摩尔分数）2根据吸收温度变化的假设，在塔内液相温度变化不大，可取平均温度35下的CO2在PC中溶解的亨利系数作为计算相平衡关系的依据。即： kPaCO2在PC中溶解的相平衡

7、关系，即：式中：为摩尔比，kmolCO2/kmolPC；为CO2的分压，kgf/cm2；T为热力学温度，K。用上述关联式计算出塔溶液中CO2的浓度有与前者结果相比要小，为安全起见，本设计取后者作为计算的依据。结论：出料（摩尔分数）2.PC密度与温度的关系利用题给数据作图，得密度与温度的关联表达式为（式中t为温度，；为密度，kg/m3）温度，（）015254055（kg/m3）122412071198118411693.PC蒸汽压的影响根据变换气组成及分压可知，PC蒸汽压与操作总压及CO2的气相分压相比均很小，故可忽略。4.PC的粘度 mPas（T为热力学温度，K）5.工艺流程确定：本次吸收采用

8、逆流吸收的方法。二、物料衡算1.各组分在PC中的溶解量查各组分在操作压力为1.6MPa、操作温度为40下在PC中的溶解度数据，并取其相对吸收饱和度均为70%，将计算所得结果列于下表（亦可将除CO2以外的组分视为惰气而忽略不计，而只考虑CO2的溶解）：CO2溶解量的计算如下：各个溶质溶解量的计算如下：（以CO2为例）通过第一部分已知CO2在40的平衡溶解度 Nm3/m3PC式中：1184为PC在40时的密度，102.09为PC的相对摩尔质量。CO2的溶解量为（10.859-0.15）0.7=7.496 Nm3/m3PC组分CO2COH2N2合计组分分压，MPa0.4640.0240.7550.3

9、571.60溶解度，Nm3/m3PC10.8590.0160.2230.22311.121溶解量，Nm3/m3PC7.4960.0110.1560.1567.819溶解气所占的百分数%95.880.141.991.99100.00说明：进塔吸收液中CO2的残值取0.15 Nm3/m3PC，故计算溶解量时应将其扣除。其他组分溶解度就微小，经解吸后的残值可被忽略。平均分子量：入塔混合气平均分子量：溶解气体的平均分子量：2.雾沫夹带量Nm3/m3PC以0.2 Nm3/m3PC计，各组分被夹带的量如下：CO2：0.20.29=0.058 Nm3/m3PCCO： 0.20.015=0.003 Nm3/m

10、3PCH2： 0.20.472=0.0944 Nm3/m3PCN2： 0.20.223=0.0446 Nm3/m3PC3.溶液带出的气量Nm3/m3PC各组分溶解量：CO2： 7.496 Nm3/m3PC 95.88%CO： 0.011 Nm3/m3PC 0.14%H2： 0.156 Nm3/m3PC 1.99%N2： 0.156Nm3/m3PC 1.99%7.819 Nm3/m3PC 100%夹带量与溶解量之和：CO2：0.028+7.203=7.554Nm3/m3PC 94.20%CO：0.003+0.011=0.014 Nm3/m3PC 0.175%H2：0.0944+0.156=0.2

11、50 Nm3/m3PC 3.12%N2：0.0446+0.156=0.201 Nm3/m3PC 2.51%8.019Nm3/m3PC 100%4.出脱碳塔净化气量以分别代表进塔、出塔及溶液带出的总气量，以分别代表CO2相应的体积分率，对CO2作物料衡算有：V1 =26058 Nm3/ h联立两式解之得V3=V1(y1-y2)/(y3-y2)=40000(0.290.006)/(0.94200.006)=12136.75Nm3/hV2 = V1 - V3 =27863.25 Nm3/ h5.计算PC循环量因每1 m3PC 带出CO2为7.554 Nm3 ，故有：L=V3y3/7.259=1213

12、6.750.9420/7.554=1513.479m3/h操作的气液比为V1/L=40000/1513.479=26.4296.验算吸收液中CO2残量为0.15 Nm3/m3PC时净化气中CO2的含量取脱碳塔阻力降为0.3kgf/cm2，则塔顶压强为16.319-0.3=16.019 kgf/cm2，此时CO2的分压为 kgf/cm2，与此分压呈平衡的CO2液相浓度为： 式中：1193为吸收液在塔顶30时的密度，近似取纯PC液体的密度值。计算结果表明，当出塔净化气中CO2的浓度不超过0.5%，那入塔吸收液中CO2的极限浓度不可超过0.216 Nm3/m3PC，本设计取值正好在其所要求的范围之内

13、，故选取值满足要求。入塔循环液相CO2：1513.4790.15=227.0227.出塔气体的组成出塔气体的体积流量应为入塔气体的体积流量与PC带走气体的体积流量之差。CO2：400000.29-7.5541513.479=167.179Nm3/h 0.60%CO： 400000.015-0.0141513.479=578.811Nm3/h 2.08%H2： 400000.472-0.2501513.479=18501.63Nm3/h 66.40%N2： 400000.223-0.2011513.479=8615.791Nm3/h 30.82%27863.411Nm3/h 100%计算数据总表

14、出脱碳塔净化气量进塔带出气量(V1)Nm3/h出塔气量(V2)Nm3/h溶液带出的总气量(V3)Nm3/h4000027863.2512136.75气液比26.429入塔气体平均分子量20.368溶解气体平均分子量42.823PC中的溶解量(溶解气量及其组成）40组分CO2COH2N2总量溶解度，Nm3/m3PC108590.02 0.22 0.22 11.32溶解量，Nm3/m3PC7.4960.01 0.16 0.16 7.819溶解体积流量Nm3/h11345.03916.648236.103236.10311833.892溶解气所占的百分数%95.880.14 1.99 1.99100

15、.00 出塔液相带出气量及其组成 40溶解量，Nm3/m3PC7.5540.014 0.25 0.20 8.019 体积流量Nm3/h11432.82121.189378.369 304.20912136.588溶解气所占的百分数%94.200.1753.122.51100.00 入塔气相及其组成 30体积流量Nm3/h11600600 188808920 40000溶解气所占的百分数%29.00 1.50 47.20 22.30 100.00 出塔气相的组成 35体积流量Nm3/h167.179578.811 18501.638615.79127863.411溶解气所占的百分数%0.602.

16、0866.4030.82100.00 入塔液相及其组成 30体积流量Nm3/h227.022227.022溶解气所占的百分数%100.00 1003、 热量衡算 在物料衡算中曾假设出塔溶液的温度为40，现通过热量衡算对出塔溶液的温度进行校核，看其是否在40之内。否则，应加大溶剂循环量以维持出塔溶液的温度不超过40。具体计算步骤如下：1.混合气体的定压比热容因未查到真实气体的定压比热容，故借助理想气体的定压比热容公式近似计算。理想气体的定压比热容：，其温度系数如下表：系数aBcdCp1（30）Cp2（32）CO24.7281.75410-2-1.33810-54.09710-98.929/37.

17、388.951/37.48CO7.373-0.30710-26.66210-6-3.03710-96.969/29.186.97/29.18H26.4832.21510-3-3.29810-61.82610-96.902/28.906.904/28.91N27.440-0.32410-26.410-6-2.7910-96.968/29.186.968/29.18表中Cp的单位为（kcal/kmol）/（kJ/kmol）进出塔气体的比热容Cpv2=Cpiyi =37.480.0060+29.180.0208+28.910.664+29.180.3082 =29.02 KJ/Kmol2.液体的比热

18、容溶解气体占溶液的质量分率可这样计算：质量分率为其量很少，因此可用纯PC的比热容代之。本设计题目中 kJ/kg文献查得 kJ/kg,据此算得： kJ/kg； kJ/kg本设计采用后者。3. CO2的溶解热kJ/kmolCO2文献查得 kJ/kmolCO2（实验测定值）本设计采用后者。CO2在PC中的溶解量为7.4961513.479=11345.039Nm3/h=506.475kmol/h故Qs=14654506.475=7421884.65kJ/h4.出塔溶液的温度设出塔气体温度为35，全塔热量衡算有：带入的热量（QV1+QL2）+ 溶解热量（Qs）= 带出的热量（QV2+QL1）Qv1=V

19、1Cpv1(Tv1T0)=4000031.4330/22.4=1683750 kJ/hQL2=L2CpL2(TL2T0)=1513.47911931.42630=20556533.39kJ/hQv2=V2Cpv2(Tv2T0)=27863.2529.0235/22.4=1263424.242kJ/hQL1=L1CpL1(TL1T0)=1828203.785 1.44TL1=1724717 TL1kJ/h式中：L1=1513.4791193+(12136.588-0.21513.479)42.823/22.4=1828203.785 kg/h1683750+20556533.39+7421884

20、.65=1263424.242+1724717TL1TL1=39.896.最终衡算结果：出塔气相及其组成（35）V2=27863.41Nm3/hCO2COH2N2167.2 578.81 18501.63 8615.79 Nm3/h0.60 2.0866.4030.82%QV2=1263424.24kJ/h入塔液相及其组成（30）L2=1513.479m3/hCO2COH2N2227.-Nm3/h-%QL2=20556533.39kJ/h入塔气相及其组成（30）V1=40000Nm3/hCO2COH2N22605811600 600 18880 8920Nm3/h29.01.547.222.3

21、%QV1=1093738kJ/h出塔液相带出气量及其组成（40）L1=12136.588kg/hCO2COH2N2694611433 21.19 378.37 304.21 Nm3/h94.200.183.122.51%脱碳塔溶解气量及其组成（40）L1=11833.892kg/hCO2COH2N27469.761134516.65 236.10 236.10 Nm3/h95.880.141.991.99%Qs=1683750kJ/h28四、设备的工艺与结构尺寸的设计计算1.确定塔径及相关参数 塔底气液负荷大，依塔底气液负荷条件求取塔径采用Eckert通用关联图法求取泛点气速，并确定操作气速。

22、入塔混合气体的质量流量:注：20.368为入塔混合气体的平均分子量10.653为出塔混合气体的平均分子量Mm2 = 440.006+280.0208+20.664+280.3082= 10.653kg/kmol塔底吸收液的质量流量:入塔混合气的密度（未考虑压缩因子）吸收液的密度（40）吸收液的粘度，依下式计算得到：mPas（平均温度35时的值）选mm塑料鲍尔环（米字筋），其湿填料因子，空隙率，比表面积，Bain-Hougen关联式常数。（1） 选用Eckert通用关联图法求解 关联图的横坐标：（v/l）0.5L/V=(12.83/1184)0.51197720/23508=5.408查Ecke

23、rt通用关联图得纵坐标值为0.0025，即：（2）选用Bain-Hougen关联式求解根据设计u=0.1m/s2求取塔径Vs=40000(0.1013/1.6)(303.15/273.15)=2810.644m3/h=0.781m3/sD=(40.781/3.140.1)0.5=3.154m本次设计取D=3200mm3核算操作气速u=4Vs/3.14D2=40.781/3.143.2=0.0972m/s则操作气体速度取u=0.10m/s合适4 核算径比D/d=3200/50=641015（满足鲍尔环的径比要求）5校核喷淋密度采用聚丙烯填料表面L喷，min=（MWR）at =0.08106.4=

24、8.512m3/(m2.h)（通常对于直径不超过75mm的散装填料，最小润湿速率可取为0.08m/（mh）;对于直径大于75mm的散装填料，可取为0.12m/（mh）L喷=（满足要求）五、填料层高度的计算塔截面积=0.785D2=0.7853.2=8.038因其他气体的溶解度很小，故将其他气体看作是惰气并视作为恒定不变，那么，惰气的摩尔流率又溶剂的蒸汽压很低，忽略蒸发与夹带损失，并视作为恒定不变，那么有L=1513.4791193/(102.0936008.038)=0.611kmol/(m2s)，吸收塔物料衡算的操作线方程为将上述已知数据代入操作线方程，整理得采用数值积分法求解，步骤如下：1

25、.将气相浓度在其操作范围内10等份，其等份间距为0.0284，并将各分点的y值代入式（1）计算出对应的x值，并列入后面表格中的第1、2列中。2.计算各分点截面处的气液相流率G=(1+Y)G L=(1+X)L 3.计算的传质系数 =1exp -1.45(33/39.1)0.75(225703/106.48.5248)0.1(2257032106.4/118421.27 108)-0.05(2257032/118439.1106.4)0.21由计算知awat=106.4式中：UL=kg/(m . h)、气体、液体的黏度，、气体、液体的密度，、溶质在气体、液体中的扩散系数， R通用气体常数，T系统温

26、度，K填料的总比表面积，填料的润湿比表面积，g重力加速度，1.27108m/h液体的表面张力，填料材质的临界表面张力，填料形状系数上述修正的恩田公式只适用于的情况，由计算得知u0.5uF 气膜吸收系数计算： 气体质量通量为 = 0.237106.510.8842.5510-4 =5.6910-3 = 0.23753.850.8842.5110-4 =2.8310-3液膜吸收系数计算： 液体质量通量为 = 0.009540.510.02497.06 = 0.896= 4.2610-3106.41.451.1 = 0.682= 0.896106.41.450.4 = 110.61故修正： =118

27、4/(102.09(1.620430+39.594)101.3=1.2910-3（稀溶液）计算准备：（1）两相摩尔流率与质量流率的转换气相平均分子量为：气相平均分子量为：33.35y+10.65VG=(33.35Y+10.65)G（稀溶液）（2）CO2在气相和液相中的扩散系数气相：分两步进行，定性温度取32.5。首先计算CO2在各组分中的扩散系数，然后再计算其在混合气体中的扩散系数。计算公式如下：DCO2-co=DCO2-H2=DCO2-N2= =（1-0.006）/（0.0208/8.6710-7+0.664/3.2810-6+0.3082/8.6210-7）=1.65410-6m2/s液相

28、：文献介绍了CO2在PC中扩散系数两个计算公式，定性温度取35。=1.1710-52/s （TK；mPas；Dcm2/s） =1.0110-5 2/s （TK；mPas；Dcm2/s）取大值（3）气液两相的粘度（纯组分的粘度）uG-CO2=1.3410-2(305.5/273.15)0.935=0.015mPas同理：uG-CO=0.018 mPas uG-H2=0.0093 mPas uG-N2=0.018 mPas为0、常压下纯气体组分的粘度，mPas 。m为关联指数（见下表）mmCO21.3410-20.935H20.8410-20.771CO1.6610-20.758N21.6610-

29、20.756气相：（气体混合物的粘度） =(0.280.015440.5+0.0250.018280.5+0.4720.009220.5+0.2230.018280.5)/(0.28440.5+0.025280.5+0.47220.5+0.223280.5)=0.0150 mPas液相： mPas=2.368 mPas（4）吸收液与填料的表面张力吸收液：=39.1 mPas填料：查教材，如聚乙烯塑料 mPas4.气相总传质单元数作CO2在PC中的相平衡曲线气相CO2的组成y（摩尔分率）0.0060 0.0500 0.1000 0.1500 0.2000 气相CO2的分压p（kgf/cm2）0.

30、0979 0.8158 1.6316 2.4473 3.2631 30对应的液相平衡组成x0.0010 0.0084 0.0168 0.0252 0.0336 35对应的液相平衡组成x1.0023 1.0196 1.0395 1.0598 1.0805 40对应的液相平衡组成x10.0537 10.4605 10.9518 11.4760 12.0369 将计算结果列表如下：因温度变化不大，故取平均温度下的数值作图得一直线，这说明CO在PC中的溶解情况满足亨利定律。但因操作关系不为直线，故仍需采用图解积分或数值积分。5.气相总传质单元数采用传质单元数的近似简化法计算图中数据源于下表数据，y、x

31、数据由操作线方程（1）计算而得。y*由y*=6.4283x-0.0002计算而得。y10-20.53.256.08.7511.5014.2517.0019.7522.5025.2528.00x10-20.0570.3120.5670.8211.0761.3311.5861.842.0952.3492.604y*10-20.3662.0033.6395.2696.9078.54410.1811.8113.4515.0816.72427.3569.5937.8625.9719.7915.9813.4111.5410.139.028.14做出图形得：现采用Smipson公式求区域的面积数值积分法（亦

32、可采用图解积分）6.气相总传质单元高度计算：由于对于PC，CO2为易吸收气体，为气膜控制10-7m填料层的有效传质高度=1.11910.828=12.117m设计高度H=1.2812.117=15.5m六、填料层的压降用Eckert通用关联图计算压降横坐标：（前已算出）纵坐标：0.00173查图得：50mmH2O/m七：附属设备及主要附件的选型1.塔壁厚 操作压力为1.6Mpa壁厚： 圆整后取27mm选用 27R钢板2液体分布器液体分布器是保持任一横截面上保证气液均匀分布。本次使用分布较好的槽盘式分布器。它具有集液、分液和分气三个功能，结构紧凑，操作弹性高，应用广泛。3除沫器除沫器用于分离塔顶

33、端中所夹带的液滴，以降低有价值的产品损失，改善塔后动力设备的操作。此次设计采用网丝除沫器。U=除沫器直径 4液体再分布器 液体向下流动时，有偏向塔壁流动现象，造成塔中心的填料不被润湿，故使用液体再分布器，对鲍尔环而言，不超过6m。故在填料3m处装一个再分布器。本次使用截锥式再分布器。5填料支撑板 填料织成板是用来支撑填料的重量，本次设计使用最为常用的栅板。本次塔径为26001400mm，使用四块栅板叠加，直径为850mm6塔的顶部空间高度 塔的顶部空间高度指顶第一层塔盘到塔顶封头的切线距离。为减少雾沫夹带的液体量，一般取1.21.5m，本次设计取1.2m八.塔附属高度及塔总高的计算塔上部空间高

34、度可取1.3m.液体分布器高度约为1.0m.考虑到气相接管所占空间高度，底部空间高度可取0.6m.若塔底液相停留时间按1min考虑，则塔釜液所占空间高度为：取h=3.2 m所以塔的附属高度为 得塔总高度: 即塔总高为21.38m九 设计概要表入塔混合气体的质量流量V36371.429 kg/h塔底吸收液的质量流量L1828203.785kg/h入塔混合气的密度p12.93Kg/m3吸收液的粘度2.368mPas填料因子120m-1空隙率0.9比表面积106.4m2/m3Bain-Hougen关联式常数A0.0942Bain-Hougen关联式常数K1.75uf0.14m/su0.1m/s塔径3

35、200mm喷淋密度L116.64m3/m2h塔截面积A5.307m2溶剂的摩尔流率L0.6068kmol/(m2s)惰气的摩尔流率G0.0438kmol/(m2s)CO2在PC中扩散系数1.1710-52/s气液两相的粘度0.015mPas吸收液与填料的表面张力39.1 mPas聚乙烯塑料的表面张力33.0 mPas气相总传质单元数NOG10.828气相总传质单元高度HOG1.119m设计高度H15.5m填料层的压降50mmH2O/m塔体壁厚27mm除沫器气速1.02m/s除沫器直径0.72m填料支撑板850mm顶部空间高度1.3m 十、设计评价及总结经过最近一周的时间的课程设计，让我受益匪浅

36、。我了解到化工单元操作的基本原理、典型设备的结构原理、操作性能和设计计算等。这次我的课程设计题目是“碳酸丙烯酯（PC）吸收过程填料吸收塔的设计”，填料塔是以塔内装有大量的填料为相接触构件的气液传质设备。填料塔具有结构较简单， 压降低，填料易用耐腐蚀材料制造等优点。本设计中，最初我采用 Eckert 通用关联图计算泛点气速，在填料的选择中，选择使用塑料鲍尔环（规格为 50mm50mm0.9mm）。在同类填料中，尺寸越小的，分离效率越高，但它的阻力将增加，通量减小，填料费用也增加很多。解决了上面的问题之后就是通过查找手册之类的书籍来确定辅助设备的选型，本设计我们所设计的填料塔持液量小，填料塔结构较

37、为简单，造价适合。不过，它的操作范围小，填料润湿效果差，当液体负荷过重时，易产生液泛，不宜处理易聚合或含有固体 悬浮物的物料等。在这次的化工原理课程设计的过程中，的确遇到了一些曾未想到的问题难题，但通过网络、书籍等，使得各项得以一一排除，这个遇到问题、查书籍、讨论、找思路、想方法的过是我们大学生的必须课。课程设计是我们在校大学生必须经过的一个过程，通过课程设计的锻炼，可以为我们即将来的毕业设计打下坚实的基础，是我们的明天更美好！参考文献：1 化工原理 上、下册 柴诚敬主编 2012.012 常用化工单元设备的设计 陈英南 刘玉兰主编 2005.043 化工AutoCAD制图应用及基础 周军，张秋利主编 2008.104 化工过程设计 王保国 王春艳 李会泉 石磊 译 2002.075 化工制图 赵慧清 蔡纪宁主编 2008.076 柴诚敬, 王军, 张缨化工原理课程设计M. 天津科学技术出版社2006.7 1化工原理课程设计 柴诚敬 刘国维 李阿娜编 2000.02