化工原理课程设计-处理量70kmolh丙烯丙烷精馏筛板塔设计(全套图纸)(30页).doc

上传人:1595****071 文档编号:37040319 上传时间:2022-08-29 格式:DOC 页数:30 大小:911KB
返回 下载 相关 举报
化工原理课程设计-处理量70kmolh丙烯丙烷精馏筛板塔设计(全套图纸)(30页).doc_第1页
第1页 / 共30页
化工原理课程设计-处理量70kmolh丙烯丙烷精馏筛板塔设计(全套图纸)(30页).doc_第2页
第2页 / 共30页
点击查看更多>>
资源描述

《化工原理课程设计-处理量70kmolh丙烯丙烷精馏筛板塔设计(全套图纸)(30页).doc》由会员分享,可在线阅读,更多相关《化工原理课程设计-处理量70kmolh丙烯丙烷精馏筛板塔设计(全套图纸)(30页).doc(30页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、-化工原理课程设计-处理量70kmolh丙烯丙烷精馏筛板塔设计(全套图纸)-第 30 页化工原理课程设计处 理 量70kmol/h丙烯丙烷精馏筛板塔设计处 理 量:70kmol/h产品质量:(以丙稀摩尔百分数计)进 料:xf65塔顶产品:xD98塔底产品: xw2安装地点:大连总板效率:0.6塔板位置:塔底塔板形式:筛板回 流 比:1.6前 言 本设计说明书包括概述、流程简介、精馏塔、再沸器、辅助设备、管路设计和控制方案共七章。 说明中对精馏塔的设计计算做了详细的阐述,对于再沸器、辅助设备和管路的设计也做了正确的说明。 鉴于本人经验有限,本设计中还存在许多错误,希望各位老师给予指正 感谢老师的

2、指导和参阅!全套图纸加153893706目录第一章 概述4第二章 流程简介5第三章 精馏塔工艺设计5第四章 再沸器的设计15第五章 辅助设备的设计21第六章 管路设计26第七章 控制方案26附录一 主要符号说明27附录二 参考文献29第一章 概述 精馏是分离过程中的重要单元操作之一,所用设备主要包括精馏塔及再沸器和冷凝器。1 精馏塔精馏塔是一圆形筒体,塔内装有多层塔板或填料,塔中部适宜位置设有进料板。两相在塔板上相互接触时,液相被加热,液相中易挥发组分向气相中转移;气相被部分冷凝,气相中难挥发组分向液相中转移,从而使混合物中的组分得到高程度的分离。简单精馏中,只有一股进料,进料位置将塔分为精馏

3、段和提馏段,而在塔顶和塔底分别引出一股产品。精馏塔内,气、液两相的温度和压力自上而下逐渐增加,塔顶最低,塔底最高。本设计为筛板塔,筛板的突出优点是结构简单、造价低、塔板阻力小且效率高。过去认为它易漏液,易堵塞。然而经长期研究发现,只要设计合理和操作适当,晒斑人能够满足生厂上所需的操作弹性,而且效率比较高。小筛孔孔径一般36mm,大筛孔孔径一般1025mm。如果采用大筛板,堵塞问题即可解决。目前已发展成为日趋广泛的一种塔板。2 再沸器作用:用以将塔底液体部分汽化后送回精馏塔,使塔内气液两相间的接触传质得以进行。本设计采用立式热虹吸式再沸器,它是一垂直放置的管壳式换热器。液体在自下而上通过换热器管

4、程时部分汽化,由在壳程内的载热体供热。3 立式热虹吸特点:循环推动力:釜液和换热器传热管气液混合物的密度差。 结构紧凑、占地面积小、传热系数高。壳程不能机械清洗,不适宜高粘度、或脏的传热介质。塔釜提供气液分离空间和缓冲区。4 冷凝器 (设计从略) 用以将塔顶蒸气冷凝成液体,部分冷凝液作塔顶产品,其余作回流液返回塔顶,使塔内气液两相间的接触传质得以进行,最常用的冷凝器是管壳式换热器。第二章 方案流程简介1 精馏装置流程精馏就是通过多级蒸馏,使混合气液两相经多次混合接触和分离,并进行质量和热量的传递,使混合物中的组分达到高程度的分离,进而得到高纯度的产品。流程如下:原料(丙稀和丙烷的混合液体)经进

5、料管由精馏塔中的某一位置(进料板处)流入塔内,开始精馏操作;当釜中的料液建立起适当液位时,再沸器进行加热,使之部分汽化返回塔内。气相沿塔上升直至塔顶,由塔顶冷凝器将其进行全部或部分冷凝。将塔顶蒸气凝液部分作为塔顶产品取出,称为馏出物。另一部分凝液作为回流返回塔顶。回流液从塔顶沿塔流下,在下降过程中与来自塔底的上升蒸气多次逆向接触和分离。当流至塔底时,被再沸器加热部分汽化,其气相返回塔内作为气相回流,而其液相则作为塔底产品采出。2 工艺流程1)物料的储存和运输精馏过程必须在适当的位置设置一定数量不同容积的原料储罐、泵和各种换热器,以暂时储存,运输和预热(或冷却)所用原料,从而保证装置能连续稳定的

6、运行。2)必要的检测手段为了方便解决操作中的问题,需在流程中的适当位置设置必要的仪表,以及时获取压力、温度等各项参数。另外,常在特定地方设置人孔和手孔,以便定期的检测维修。3)调节装置由于实际生产中各状态参数都不是定值,应在适当的位置放置一定数量的阀门进行调节,以保证达到生产要求,可设双调节,即自动和手动两种调节方式并存,且随时进行切换。3 设备选用 精馏塔选用筛板塔,配以立式热虹吸式再沸器。4 处理能力及产品质量处理量: 70kmol/h产品质量:(以丙稀摩尔百分数计)进料:xf65塔顶产品:xD98塔底产品: xw2第三章 精馏塔工艺设计 第一节 设计条件1.工艺条件:饱和液体进料,进料丙

7、稀含量xf65(摩尔百分数)塔顶丙稀含量 xD98,釜液丙稀含量 xw2,总板效率为0.6。2操作条件:1) 塔顶操作压力:P=1.62MPa(表压)2) 加热剂及加热方法:加热剂水蒸气加热方法间壁换热3) 冷却剂:循环冷却水4) 回流比系数:R/Rmin=1.63塔板形式:筛板4处理量:qnfh=70kmol/h5安装地点:大连6塔板设计位置:塔底7. 物料密度:MA=42 MB=44 第二节 物料衡算及热量衡算一 物料衡算1换算将摩尔百分数换算成质量百分数:W=XMA/XMA+(1-X)MBxf65 wf63.93xD98 wD97.91xw2 wW1.91 将摩尔流量换算成质量流量:进料

8、状态混合物平均摩尔质量:(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 再沸器加热蒸气的质量流量:G

9、R= QR/rR2) 冷凝器热流量:QC=Vr冷凝器冷却剂的质量流量:GC= QC/(cl(t2-t1)假设初值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第三节 塔板数的计算利用程序进行迭代计算:流程图如下:泡点初值To总压PAntoinec常数outputy-1p0A p0BKA KB泡点计算2塔板数计算注:下标t、b分别表示塔顶、塔底参数。 计算过程包括:假设塔顶温度Tt0=31

10、6K 经泡点迭代计算得塔顶温度Tt=316.13K塔顶压力Pt=1620+101.3=1721.3KPa代入公式 计算并换算得:PAo=1726.955KPa ; PBo=1444.432KPa又 得:KA=1.003 ; KB=0.839BAKK=a1 1.1956设a2 1.1888=(a1+a2)/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(-

11、=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 物性数据常

12、压43下,丙烷的物性数据:气相密度:V =29.6kg/ m3液相密度:L =466kg/ m3液相表面张力:=4mN/m2 初估塔径气相流量:qmVs=6.9639kg/s qVVs=qmVs/v=0.2487m3/s液相流量:qmLs=6.4275kg/s qVLs=qmLs/L=0.01368 m3/s两相流动参数: =0.2190初选塔板间距 HT=0.45m,查化工原理(下册)P107筛板塔泛点关联图,得:C20=0.061所以,气体负荷因子: =0.04526 液泛气速: 0.1738m/s 取泛点率0.7 操作气速:u = 泛点率 uf=0.1216m/s 气体流道截面积: =2

13、.0452m2 选取单流型弓形降液管塔板,取Ad / AT=0.15; 则A / AT=1- Ad / AT =0.85 截面积: AT=A/0.85=2.406 m2 塔径: =1.75m 圆整后,取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.94573 塔高的估算 Np=112 有效高度:Z= HT Np=50.4m 釜液高度(略),进料处

14、两板间距增大为0.7m设置4个人孔,每个人孔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.04

15、5m液体流经底隙的流速: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第七节 塔板流

16、动性能校核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

17、.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

18、 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.67

19、25kg/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 =35kg/m3 蒸气压曲线斜率(t/P)=0

20、.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=800mm L/D=3.75 取 管程进口直径:Di=0.25m 管程出口直径:Do=0.3m 三 传热系数的校核1显热段传热系数K假设传热管

21、出口汽化率 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) 污垢热阻及管壁热阻 沸

22、腾侧: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)泡核沸腾压抑因数

23、:=(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

24、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)釜液进口管内流动雷诺数: = 2

25、630462进口管内流体流动摩擦系数: = 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 动量变化引起的阻力系数:

26、= 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第五章 辅助设备设计 一 辅

27、助容器的设计 容器填充系数取: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塔顶产品罐

28、质量流量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.99

29、kw 壳程水焓变: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管程流率:q

30、mDs = 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)则传热面积:

31、 圆整后 取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

32、查得:=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第七章 控制方案 精馏塔的控制方案要求从质量指标、产品产量和能量消耗三个方面进行综合考虑。精馏塔最直接的质量指标是产品浓度。由于检测上的困难,难以直接按产品纯度进行控制。最常用的间接质量指标是温度。 将本设计的

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 教育专区 > 小学资料

本站为文档C TO C交易模式,本站只提供存储空间、用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。本站仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知淘文阁网,我们立即给予删除!客服QQ:136780468 微信:18945177775 电话:18904686070

工信部备案号:黑ICP备15003705号© 2020-2023 www.taowenge.com 淘文阁