过程工艺与设备课程设计-丙烯---丙烷精馏装置设计70kmol 1.6.doc

- 49 -过程工艺与设备课程设计任务书丙烯-丙烷精馏装置设计学生姓名：班级：应化0612班学号：指导老师：完成时间：2008.7.518前言化工原理是化工及其相关专业学生的一门重要的技术基础课，其课程设计涉及多学科知识，包括化工，制图，控制，机械等各种学科，是一项综合性很强的工作；是锻炼工程观念和培养设计思维的好方法，是为以后的各种设计准备条件；是化工原理教学的关键环节，也是巩固和深化理论知识的重要环节。 本设计说明书包括概述、方案流程简介、精馏塔、再沸器、辅助设备、管路设计和控制方案共七章。 说明中对精馏塔的设计计算做了较为详细的阐述，对于再沸器、辅助设备和管路和控制方案的设计也做了简要的说明。 在设计过程中，得到了董宏光老师和肖武老师的指导，得到了同学们的帮助，同学们一起讨论更让我感受到设计工作是一种集体性的劳动，避免了不少错误。 鉴于本学生的经验和知识水平有限，设计中难免存在错误和不足之处，请老师给予指正 感谢老师的指导和参阅！目录 前言2任务书5第一章 概述111.1 精馏塔111.2 再沸器111.3 冷凝器12第二章 方案流程简介132.1 精馏装置流程132.2 工艺流程132.3 设备选用142.4 处理能力及产品质量14第三章 精馏过程系统设计153.1设计条件153.2物料衡算及热量衡算153.3塔板数的计算173.4精馏塔工艺设计193.5溢流装置的设计203.6塔板流动性能校核223.7负荷性能图25第四章 再沸器的设计284.1设计任务与设计条件284.2估算设备尺寸294.3传热系数的校核304.4循环流量校核33第五章 辅助设备的设计385.1冷凝器385.2其它换热设备385.3容器405.4 管路设计及泵的选择41第六章 管路设计45第七章 控制方案48设计心得及总结 49附件一 程序及运行结果51附录一主要符号说明54附录二 参考文献56 第一章 概 述 精馏是分离分离液体混合物最常用的一种单元操作，所用设备主体核心设备是精馏塔，辅助设备包括再沸器、冷凝器、储罐、预热器及冷却器。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/h5安装地点：大连 6塔板设计位置：塔顶第二节 物料衡算及热量衡算一 物料衡算1.单位换算（1）将摩尔百分数换算成质量百分数W=XMA/XMA+(1-X)MBxf65 wf63.93xD98 wD97.91xw2 wW1.91 （2）将摩尔流量换算成质量流量进料状态混合物平均摩尔质量：（MA为丙稀摩尔质量 MB为丙烷摩尔质量）M=xfMA+(1-xf)MB=0.65420.3544=42.7kg/kmol进料状态下的质量流量：qmfs=Mqnfh/3600=0.8303kg/s2.全塔物料衡算：= + =+=70 kmol/h，65%， 2 98解得：=45.9375koml/h ,= 24.0625koml/h.二 塔内气、液相流量塔内气、液相流量：1）精馏段：； 2）提馏段： 三 热量衡算1） 再沸器加热蒸气的质量流量：2） 冷凝器热流量： 冷凝器冷却剂的质量流量： 第三节 塔板数的计算相对挥发度的计算假设塔顶温度t=40.0 C 塔顶压力Pt=1.72MPa查P-K-T图得：kA=0.98 ；kB=0.90 则顶=kA/kB=0.98/0.90=1.089 ；假设精馏塔的塔板数是159块，每块板的压降为100mmH2O，液态丙烷密度近似为530kg/m3; 塔底压力为P=1.80Mpa;假设塔底温度为t=50C, kA=1.15 ；kB=1.10;则底=kA/kB=1.15/1.10=1.15 =（顶+底）/2=1.119；利用程序进行迭代计算：流程图如下：计算过程包括：给定平均相对挥发度：=1.119饱和液体即泡点状态进料：q=1 q线：x=xf 代入数据，解得 =12.13 R=1.6Rmin=19.41为逐板计算过程：y1=xD=0.98ynynxn)1(-=aa直至xi xf 理论进料位置：第i块板进入提馏段：ynynxn)1(-=aa 直至xn6mm取堰高hw=0.04m,底隙hb=0.040m液体流经底隙的流速：ub =0.39m/s0.5m/s 符合要求3.浮阀数及排列方式 1浮阀数 选取F1型浮阀，重型，阀孔直径d0=0.039 m初选阀动能因子F0=10，计算阀孔气速=1.87 浮阀个数 =176 2.浮阀排列方式 进出口安全宽度bs=bs=60mm边缘区宽度bc=50mm=0.313m = 0.627m = 0.95m有效传质面积：有效传质面积： =2.195选择等腰三角形排列，按t=125mm进行布孔，实排阀数n=168个重新计算塔板以下参数： 阀孔气速 =1.96m/s 动能因子 F1=2.078=10.5浮阀的开孔率 6.39%12%第六节 塔板流动性能校核1 液沫夹带量校核 =0.3040.8 由塔板上气相密度及塔板间距查化工单元过程及设备课程设计书图5-26得系数=0.134根据表5-11所提供的数据，K可取K=1。故不会产生过量的液沫夹带。Z=D-2=01.374m=2.512m2 =0.4220.82 塔板阻力hf的计算和核对塔板阻力hf= ho+hl+h（1）干板阻力ho临界气速 =1.469 Hd 所以不会发生液泛。4 液体在降液管中的停留时间液体在降液管中的停留时间应大于35s =8.129s5s 满足要求,则可避免严重的气泡夹带。 5 严重漏液校核当阀孔的动能因子低于5时将会发生严重漏夜，故漏液点的气速可取=5的相应孔流气速 =0.933 m/s =2.101.5 满足稳定性要求 第七节 负荷性能图以气相流量为纵坐标，液相流量为横作标1 过量液沫夹带线根据前面液沫夹带的较核选择表达式：由此可得液沫夹带线方程：=0.3941 -4.263 此线记作线（1）2 液相下限线对于平直堰，其堰上液头高度必须大于0.006m，取=0.006m ，即可确定液相流量的下限 取E=1，代入 lw，可求得lw的值，则 qnLh=3.07*lw=3.30m/h此线记作线（2）-与纵轴平行3. 严重漏液线当阀孔的动能因子低于5时将会发生严重漏夜，故取时，计算相应气相流量则/5 =275.2此线记作线（3） 与横轴平行4. 液相上限线 59.9由上述关系可作得线（4）5. 降液管液泛线Hd=HT+hW令 将 其中 =0 为避免降液管液泛的发生，应使 （*）。其中 =0.05 hf= ho+hl+h其中h可忽略不记 将各式代入（*）式可得液泛方程线：1.87*10-8 *=0.2087-2.701* qLh2/3-6.35* qLh 液相流量1020304050气相流量1023 998973945924第四章 再沸器的设计一 设计任务与设计条件1选用立式热虹吸式再沸器 其壳程以水蒸气为热源，管程为塔底的釜液。釜液的组成为（摩尔分数）丙烯=0.02，丙烷=0.98塔顶压力：1.72MPa塔底压力Pw=1.80 MPa2再沸器壳程与管程的设计壳程管程温度（）10050压力（MPa绝压）0.10131.803 物性数据1） 壳程凝液在温度（100）下的物性数据：潜热：rc=2258.4热导率：c =0.683w/(m*K)粘度：c =0.283mPas密度：c =958.4kg/m32） 管程流体在（50 1.70MPa）下的物性数据：潜热：rb=298液相热导率：b =0.086w/(mK)液相粘度：b =0.072mPas液相密度：b =446.1kg/m3 液相定比压热容：Cpb=3.2K 表面张力：b0.00301N/m气相粘度：v =0.008mPas气相密度：v =28.2kg/m3 蒸气压曲线斜率（t/P）=0.000244 m2 K/kg二 估算设备尺寸热流量： = MwVrb1000/3600= 3522956w传热温差： =50 假设传热系数：K=780W/( m2 K)估算传热面积Ap = 90.33 m2 拟用传热管规格为：252mm,管长L=3m则传热管数： =256 若将传热管按正三角形排列参考化工单元过程及设备课程设计P94表3-6得NT=301根b=19.08则 壳径： =653.7mm 取 D= 0.800m 取 管程进口直径：Di=0.50m 管程出口直径：Do=0.65m 三 传热系数的校核1显热段传热系数K假设传热管出口汽化率 Xe=0.23（对于丙烯丙烷取0.25左右）则循环气量： =51.4kg/s1) 计算显热段管内传热膜系数hi传热管内质量流速： di=25-22=21mm= 493.02kg/( m2 s)雷诺数： = 143798.710000普朗特数： =2.68 显热段传热管内表面系数： = 1868w/( m2 K)2) 壳程冷凝传热膜系数计算ho蒸气冷凝的质量流量： = 1.56kg/s传热管外单位润湿周边上凝液质量流量： =0.066 kg/(m s) = 932.6630%所以，传热面积裕度合适，满足要求四 循环流量校核1循环系统推动力：1）当X=Xe/3= 0.087时=2.94 两相流的液相分率： = 0.348两相流平均密度： = 173.6kg/m3 2）当X=Xe=0.23 = 0.929两相流的液相分率： = 0.201两相流平均密度： = 122.2g/m3根据课程设计表319 得：L=1.02m,则循环系统的推动力： =10657.8pa 2循环阻力Pf： 管程进出口阻力P1 进口管内质量流速： =261.78kg/(m2s)釜液进口管内流动雷诺数： = 1817903进口管内流体流动摩擦系数： =0.0154进口管长度与局部阻力当量长度：=58.02m管程进出口阻力: =137.48Pa 传热管显热段阻力P2 =493.02kg/(m2s) =143798.7=0.0205 = 24.323Pa 传热管蒸发段阻力P3 a. 气相流动阻力Pv3 G=493.02kg/(m2s) 取X=2/3Xe 则 =75.6kg/(m2s) =198442.2 =0.0196=416.7Pab. 液相流动阻力PL3GL=G-Gv=417.43kg/(m2s) =1095746.1=0.0161=659.8Pa = 8445.5Pa 管内动能变化产生阻力P4 动量变化引起的阻力系数:= 3.0 =1633.2管程出口段阻力P5 a. 气相流动阻力Pv5 = 154.9kg/(m2s) =35.6kg/(m2s) 管程出口长度与局部阻力的当量长度之和： = 75.26m =2894661.2 =0.0149=7.68pab. 液相流动阻力PL5=119.3 kg/(m2s) = 1076758.0=0.0161= 29.76Pa = 256.08Pa所以循环阻力：Pf=P1 + P2 + P3 + P4 + P5=10496.6pa 又因PD=10657.8Pa 所以 =1.015循环推动力略大于循环阻力，说明所设的出口汽化率Xe基本正确，因此所设计的再沸器可以满足传热过程对循环流量的要求。 第五章 其它辅助设备的选型第一节 冷凝器拟用10水为冷却剂，出口温度为30。走壳程。管程温度为40tm=18.2管程流率：qmVs=11.285kg/s取潜热r=302.54kJ/kg传热速率：Q= qmVsr=3414.18kw壳程取焓变：H=125.5kJ/kg则壳程流率：qc=Q/H=97936.6kg/h假设传热系数：K=1000 w/(m2K)则传热面积：A=Q/Ktm=187.63 m2圆整后 取A=188m2第二节 其它换热设备1进料预热器 用90水为热源，出口约为70走壳程 料液由20加热至50，走管程 传热温差：tm=44.81 管程液体流率：qmfh=3600 qmfs=2989kg/h 管程液体焓变：H=401kJ/kg 传热速率：Q= qmfsH=2349401/3600=261.65kw 壳程水焓变：H=125.6kJ/kg 壳程水流率：q=7496.9kg/h 假设传热系数：K=1000w/(m2K) 则传热面积：A=Q/Ktm=5.84 m2 圆整后 取A=6 m2 1 塔顶产品冷却器拟用10水为冷却剂，出口温度为20。走壳程。管程温度由40降至25 tm=13.95 管程流率：qmDs = 0.5359kg/s ; 取潜热：r=280kJ/kg则传热速率：Q= qmDsr=150.04kw壳程焓变：H=84.0kJ/kg则壳程流率：qc=Q/H=6430.2kg/h假设传热系数：K=1000 w/(m2K)则传热面积A=Q/Ktm=10.76 m2 圆整后 取A=11m22 釜液冷却器拟用10水为冷却剂，出口温度为20。走壳程。管程温度由50.0降到25tm=16.83 管程流率：qmWs=0.2941kg/s丙烷液体焓变：H =282kJ/kg传热速率：Q= qmWsH =82.9268kw壳程取焓变：H=84.0kJ/kg则壳程流率：qc=Q/H=3554.01kg/h假设传热系数：K=1000 w/(m2K)则传热面积：A=Q/Ktm=4.93 m2 圆整后 取A=5m2第三节 容器容器填充系数取：k=0.7 1进料罐（常温贮料） 20丙稀 L1 =526kg/m3 丙烷 L2 =500kg/m3 压力取1.8244MPa 由上面的计算可知 进料 Xf=65% Wf=63.93% 则 =516.3 kg/m3 进料质量流量：qmfh=3600 qmfs=2989kg/h 取 停留时间：x取为4天，即x=96h 进料罐容积： 793.96m3 圆整后 取V=794 m3 2回流罐（40）质量流量qmLh=3600RqmDs =38697kg/h设凝液在回流罐中停留时间为0.25h，填充系数取=0.7则回流罐的容积 26.27 m3取V=27 m33塔顶产品罐质量流量qmDh=3600qmDs =1929.06kg/h；产品在产品罐中停留时间为72h，填充系数取=0.7则产品罐的容积 377.22m3取V=378m34 釜液罐取停留时间为取5天，即x=120h质量流量qmWh=3600qmWs =1058.64 kg/h 则釜液罐的容积 362.96 m3取V=363m3第四节 管路设计及泵的选择一、泵的设计1进料泵(两台，一用一备)取液体流速：u=0.5m/s液体密度： kg/ m3 qVfs = qmfs / =0.00161 m3/s 0.064m取d=65mm液体粘度 取=0.2相对粗糙度：/d=0.003 246761查得：=0.025取管路长度：l=100m 取90度弯管4个，截止阀一个，文氏管流量计1个取则qVLh =5.9730m3/h选取泵的型号：AY 扬程：3065m 流量：2.560m3 /s2回流泵（两台，一开一用）取液体流速：u=0.5m/s液体密度： kg/ m3 qVLs = qmLs / =0.0229m3/s 0.241m液体粘度 取=0.2相对粗糙度：/d=0.00083 858106查得：=0.02取管路长度：l=100m 取90度弯管4个，截止阀一个，文氏管流量计1个取则qVLh =82.11m3/h选取泵的型号：DSJH 扬程：38280m 流量：951740m3 /s3.釜液泵（两台，一开一用）取液体流速：u=0.4m/s液体密度： kg/ m3 qVWs = qmWs / =0.000665m/s 0.046m液体粘度 取=0.2相对粗糙度：/d=0.00431066593查得：=0.04取管路长度：l=40m 取90度弯管4个，截止阀一个，文氏管流量计1个取则qVLh =2.393m3/h该处泵扬程为负值，正常工作时不使用，但非正常工作或停止工作时，需要使用。选取泵的型号：GI 扬程：101510m 流量：0.190m3 /s二、管路设计进料管线取料液流速：u=0.5m/s则0.064m取管子规格685。其它各处管线类似求得如下：名称管内液体流速（m/s）管线规格（mm）进料管0.5684塔顶蒸气管151218塔顶产品管0.5603回流管0.5684釜液流出管0.5455仪表接管/323塔底蒸气回流管151146第六章 管路设计1进料管线取料液流速：u=0.5m/s 体积流量V=0.00161则=0.064m取管子规格703的管材。其内径为0.064 m2塔顶蒸汽管：取原料流速：u=12m/s 体积流量：V=952则=0.134 m取管子规格1528.5 . 其内径为0.135m，其实际流速为u=4.3m/s3. 塔顶产品管取原料流速u=0.4m/s，其体积流量：V=4.07则=0.060m取管子规格684. 其内径为0.060 m，其实际流速为u=0.4m/s4. 回流管取原料流速：u=0.7m/s 体积流量：V=67则=0.135m取管子规格1528.5 . 其内径为0.135m，其实际流速为u=1.147m/s5釜液流出管取原料流速：u=0.3m/s 体积流量：V=2.387则=0.053 m取管子规格603.5. 其内径为0.053 m。6仪表接管选管规格：323 .7塔底蒸汽回流管取原料流速：u=20m/s 体积流量：qnvs=738.68则=0.035 m取管子规格374, 其内径为37mm，所求各管线的结果如下：名称管内液体流速（m/s）管线规格（mm）进料管0.5703顶蒸气管4.31528.5顶产品管0.4684回流管1.151528.5釜液流出管0.3603.5仪表接管/323塔底蒸气回流管101528.5第七章 控制方案 精馏塔的控制方案要求从质量指标、产品产量和能量消耗三个方面进行综合考虑。精馏塔最直接的质量指标是产品浓度。由于检测上的困难，难以直接按产品纯度进行控制。最常用的间接质量指标是温度。 将本设计的控制方案列于下表序号位置用途控制参数介质物性L(kg/m3)1FIC-01进料流量控制03000kg/h丙烷丙稀L=513.92FIC-02回流定量控制01500kg/h丙稀L=4703PIC-01塔压控制02MPa丙稀V=36.24HIC-02回流罐液面控制01m丙稀L=4705HIC-01釜液面控制03m丙烷L=4506TIC-01釜温控制4060丙烷L=44设计心得及总结这次课程设计的时间虽然不长，只有短短的两个星期，可以说是在这两个星期里的日子是痛苦、曲折又有着深刻的意义。在这次设计过程中我学到了在理论课上根本学不到的东西，同时也对自己所学的化工原理理论知识进行了巩固和深化，对于一些其他软件例如CAD制图软件和C语言编程有了进一步的学习。通过这次课程设计让我对自己是一名学化学的学生的身份进行了更深的思考。有些只是在书本上得理论知识如果运用到实际生产中时，会面临很多的问题，对于可行性还要进行深层次的讨论。我们应用化学是属于理科学生，对于一些工程问题要求得不如功课学生严格，但是我们也要严个要求自己，将自己所学的理论知识运用到实际生产过程中。这次我的任务是设计丙烯丙烷的精馏装置，涉及到精馏塔的工艺设计计算；塔底再沸器的设计计算； 管路尺寸的确定、管路阻力计算及泵的选择；其余辅助设备的计算及选型；控制仪表的选择参数；以及用3#图纸绘制带控制点的工艺流程图及主要设备(精馏塔或再沸器)的工艺条件图一张。这些知识在我们的化工原理理论课上根本就不要求的，考试也不会考的，但通过这次学习然我们对这方面的知识有了更深一层的了解。至少我们会运用自己所学的理论知识到工程操作中去。同时这次课程设计对于我们的自学水平也是一个挑战，有的知识我们需要自己去探索、思考、学习，对于以前没有掌握的软件也要学习，这无疑都是锻炼我们得好机会。记得在这次课程设计时，当老师讲课时就说过我们听完课以后刚开始设计时会感觉到无从下手，当我们真的实践时真的是无从下手，大脑一片空白，不知道从那个地方开始做，但是经过同学的讨论和认真学习，一步一步设计也就出来成果了。在各种计算和参数选择时，常常是得选择好多次，有的时候进退两难，有的时候是无从下手，当结果出来时发现不符合要求，那种滋味很痛苦，还得接着算下去。如果想得到正确结果，想学到真正的知识就得坚持。只要付出努力，我相信总有见到光明得一天。我们错了不要紧，关键是我们要从错误中吸取经验，以至于不再犯同样的错误。想想我们平时考试的题目真的是很简单，那些参数都是提前给好的，只要我们计算一些简单的计算而已，理论知识与实践真的是相差甚远，我们得不断的学习与思考。当我们完成设计时，心里无比的骄傲与自豪，有一种成就感。想想终于把自己学到的理论知识能应用到实际工程中了，将理论知识与实践初步的结合起来了。就像是前些日子我们刚刚实习完的时候，知道自己的专业知识是如何才能应用到世实际中一样，有种喜悦感。由于本人的专业知识和经验水平都有限，在设计中会存在一定的缺陷和不足，希望老师给与指正，在这里再一次谢谢老师的指导。通过这次的设计学习，让我更加热爱自己的专业。也感觉到我们身上的重任，作为21世纪得大学生，我们有责任担当起建设我国的重任，利用自己所学的专业知识建设我国。我们要抓住这个历练自己的好机会，为以后走上社会工作岗位打下坚实的基础。附件一 程序及运行结果# include # include main() double qnf=50, zf=0.65, q=1, xD=0.98, xw=0.02, a=1.138, qnD,qnw,xq,yq,hb=1.4,R,Rmin,x,y;int NT,NF,j;xq=zf,yq=a*xq/(1+(a-1)*xq);Rmin=(xD-yq)/(yq-xq);R=hb*Rmin;qnD=qnf*(zf-xw)/(xD-xw);qnw=qnf*(zf-xD)/(xw-xD);y=xD; x=y/(a-(a-1)*y);printf(y1=%fx1=%fn,y,x);for(j=2;xxw;j+)if(xxq) y=R*x/(R+1)+xD/(R+1);NF=j+1; Else y=(R*qnD+q*qnf)*x/(R*qnD+q*qnf-qnw)-qnw*xw/(R*qnD+q*qnf-qnw);x=y/(a-(a-1)*y);printf(y%d=%fx%d=%fn,j,y,j,x); NT=j-1;printf(NT=%dNF=%dqnD=%fqnw=%fn,NT,NF,qnD,qnw);printf(qnl=%fmol/hqnv=%fmol/hn,R*qnD,(R+1)*qnD