丙烯-丙烷体系对筛板塔顶精馏塔的设计（处理量：140回流比系数：1.3）

课程设计 生化0203 屈明博 200241009过程工艺与设备课程设计（精馏塔及辅助设备设计） 班 级： 生化0201班 姓 名： 屈明博 指导老师： 目录概述3流程简介5精馏塔工艺设计7再沸器的设计21辅助设备的设计29管路设计36控制方案37附件一 C程序 138附件二 matlab程序244附录一 主要符号说明44附录二 参考文献47概述 精馏是分离过程中的重要单元操作之一，所用设备主要包括精馏塔及再沸器和冷凝器。精馏塔精馏塔是一圆形筒体，塔内装有多层塔板或填料，塔中部适宜位置设有进料板。两相在塔板上相互接触时，液相被加热，液相中易挥发组分向气相中转移；气相被部分冷凝，气相中难挥发组分向液相中转移，从而使混合物中的组分得到高程度的分离。简单精馏中，只有一股进料，进料位置将塔分为精馏段和提馏段，而在塔顶和塔底分别引出一股产品。精馏塔内，气、液两相的温度和压力自上而下逐渐增加，塔顶最低，塔底最高。本设计为筛板塔，筛板的突出优点是结构简单、造价低、塔板阻力小且效率高。但易漏液，易堵塞。然而经长期研究发现其尚能满足生产要求，目前应用较为广泛。再沸器作用：用以将塔底液体部分汽化后送回精馏塔，使塔内气液两相间的接触传质得以进行。本设计采用立式热虹吸式再沸器，它是一垂直放置的管壳式换热器。液体在自下而上通过换热器管程时部分汽化，由在壳程内的载热体供热。立式热虹吸特点：循环推动力：釜液和换热器传热管气液混合物的密度差。 结构紧凑、占地面积小、传热系数高。壳程不能机械清洗，不适宜高粘度、或脏的传热介质。塔釜提供气液分离空间和缓冲区。冷凝器 （设计从略） 用以将塔顶蒸气冷凝成液体，部分冷凝液作塔顶产品，其余作回流液返回塔顶，使塔内气液两相间的接触传质得以进行，最常用的冷凝器是管壳式换热器。第二章 方案流程简介精馏装置流程 精馏就是通过多级蒸馏，使混合气液两相经多次混合接触和分离，并进行质量和热量的传递，使混合物中的组分达到高程度的分离，进而得到高纯度的产品。 流程如下： 原料（丙稀和丙烷的混合液体）经进料管由精馏塔中的某一位置（进料板处）流入塔内，开始精馏操作；当釜中的料液建立起适当液位时，再沸器进行加热，使之部分汽化返回塔内。气相沿塔上升直至塔顶，由塔顶冷凝器将其进行全部或部分冷凝。将塔顶蒸气凝液部分作为塔顶产品取出，称为馏出物。另一部分凝液作为回流返回塔顶。回流液从塔顶沿塔流下，在下降过程中与来自塔底的上升蒸气多次逆向接触和分离。当流至塔底时，被再沸器加热部分汽化，其气相返回塔内作为气相回流，而其液相则作为塔底产品采出。工艺流程物料的储存和运输 精馏过程必须在适当的位置设置一定数量不同容积的原料储罐、泵和各种换热器，以暂时储存，运输和预热（或冷却）所用原料，从而保证装置能连续稳定的运行。必要的检测手段 为了方便解决操作中的问题，需在流程中的适当位置设置必要的仪表，以及时获取压力、温度等各项参数。 另外，常在特定地方设置人孔和手孔，以便定期的检测维修。3） 调节装置由于实际生产中各状态参数都不是定值，应在适当的位置放置一定数量的阀门进行调节，以保证达到生产要求，可设双调节，即自动和手动两种调节方式并存，且随时进行切换。设备选用 精馏塔选用筛板塔，配以立式热虹吸式再沸器。处理能力及产品质量处理量： 70kmol/h产品质量：（以丙稀摩尔百分数计）进料：xf65塔顶产品：xD98塔底产品: xw2第三章 精馏塔工艺设计 第一节 设计条件工艺条件：饱和液体进料，进料丙稀含量xf65（摩尔百分数）塔顶丙稀含量 xD99，釜液丙稀含量 xw1，总板效率为0.6。 2操作条件：1）塔顶操作压力：P=2.5MPa（表压）2）加热剂及加热方法：加热剂水蒸气 加热方法间壁换热3）冷却剂：循环冷却水(海水)4）回流比系数：R/Rmin=1.3 3塔板形式：筛板 4处理量：qnfh=140kmol/h 5安装地点：大连 6塔板设计位置：塔顶 第二节 物料衡算及热量衡算一 物料衡算1换算将摩尔百分数换算成质量百分数：W=XMA/XMA+(1-X)MBxf65 wf63.93xD98 wD97.91xw2 wW1.91 将摩尔流量换算成质量流量：进料状态混合物平均摩尔质量：（MA为丙稀摩尔质量 MB为丙烷摩尔质量）M=xfMA+(1-xf)MB=0.65420.3544=42.7kg/kmol进料状态下的质量流量：qmfs=Mqnfh/3600= 0.9489kg/s 2求质量流量 qmDs + qmws = qmfsqmDswD + qmwswW = qmfswf解得： qmDs = 0.613kg/s ; qmws= 0.3358kg/s 塔内气、液相流量：1）精馏段：L =RD; V =(R+1)D;2）提馏段：L=L+qF; V=V-(1-q)F; L=V+W; 二 热量衡算1）再沸器热流量：QR=Vr 再沸器加热蒸气的质量流量：GR= QR/rR冷凝器热流量：QC=Vr冷凝器冷却剂的质量流量：GC= QC/(cl(t2-t1)第三节 塔板数的计算利用程序进行迭代计算：流程图如下：1泡点计算：泡点初值To总压PAntoinec常数outputy-1p0A p0BKA KB2塔板数计算：假设初值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注：下标t、b分别表示塔顶、塔底参数。计算过程包括：假设塔顶温度Tto=316K 经泡点迭代计算得塔顶温度Tt=316.73K塔顶压力Pt=1620+101.3=1721.3KPa 代入公式 计算并换算得：PAo=1726.955KPa ; PBo=1444.432KPa又 得：KA=1.003 ; KB=0.839BAKK=a1 1.1956设2 1.1881=(1+2)/2 1.1922计算过程包括：泡点进料：q=1 q线：x=xf 代入数据，解得 xe=0.65;ye=0.6889 =7.489 R=1.4Rmin=10.45 为逐板计算过程：y1=xD=0.98ynynxn)1(-=aa直至xi xf 理论进料位置：第i块板进入提馏段：ynynxn)1(-=aa 直至xn xW 计算结束。理论板数：Nt=n（含釜）（具体程序见附件一）迭代结果： 进料板Nf=i/0.6+1=58, 实际板数Np= (Nt-1)/0.6 =125则塔底压力Pb=Pt+0.980.47Np= 1772.9KPa同可算得：塔底温度Tb=326.64K 2 1.1881 符合假设 所以假设成立，上述计算结果均为正确结果。塔内气、液相流量： 精馏段： qmLs=RqmDs=6.4058kg/s qmVs=(R+1)qmDs=7.01885kg/s 提馏段 ： q,mLs=qmLs+qmFs=7.3547kg/s q,mVs= qmVs =7.01885kg/s 第四节 精馏塔工艺设计物性数据常压43下，丙稀的物性数据：气相密度：V =28kg/ m3液相密度：L =470kg/ m3液相表面张力：=4.5mN/m初估塔径 气相流量：qmVs=6.9639kg/s qVVs=qmVs/v=0.2506737m3/s 液相流量：qmLs=6.4275kg/s qVLs=qmLs/L=0.01368 m3/s两相流动参数： =0.2227 初选塔板间距 HT=0.45m,查化工原理（下册）P107筛板塔泛点关联图，得：C20=0.061所以，气体负荷因子： =0.04526 液泛气速： 0.1798m/s 取泛点率0.7 操作气速：u = 泛点率 uf=0.1259 m/s 气体流道截面积： =1.9754 m2 选取单流型弓形降液管塔板，取Ad / AT=0.15; 则A / AT=1- Ad / AT =0.85 截面积: AT=A/0.85=2.324 m2 塔径： =1.72m 圆整后，取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.6394塔高的估算 Np=112 有效高度：Z= HT Np=56.25m 釜液高度（略），进料处两板间距增大为0.7m设置8个人孔，每个人孔0.8m裙座取5m,塔顶空间高度1.5m,釜液上方气液分离高度取1.5m.设釜液停留时间为30min釜液高度：Z =0.44m 取其为0.5m 所以，总塔高h=Z+0.7+5+1.5+1.5+0.5=65.45m 第五节 溢流装置的设计降液管（弓形）由上述计算可得：降液管截面积：Ad=AT0.15= 0.3817 m2由Ad/AT=0.15,查化工原理（下册）P113的图6.10.24可得：lw/D=0.82所以，堰长lw=0.82D=1.476m溢流堰 取E近似为1则堰上液头高： =0.029428m取堰高hw=0.04m,底隙hb=0.045m液体流经底隙的流速：ub =0.206m/s0.3m/sub0.45m/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第七节 塔板流动性能校核液沫夹带量校核 Hf=2.5(hw+how)=0.1736m质量夹带率ev=0.000076kg液/kg气ev5s 满足要求 严重漏液校核 =0.0141m 满足稳定性要求 1.1257m/s第八节 负荷性能图过量液沫夹带线规定：ev = 0.1（ kg 液体 / kg气体） 为限制条件得： = 10671-167qVLh2/3 由上述关系可作得线液相下限线 整理出：qVLh=3.07lw=4.034 与y轴平行 由上述关系可作得线严重漏液线 将下式分别代入 近似取Co为前面计算的值 得：qVVh =a(b+cqVLh2/3)1/2其中：a= =4843.3 b=0.0056+0.13hw-h=0.01024 c= =0.000285得：qVVh =4843.3(0.01024+0.000285qVLh2/3)1/2 由上述关系可作得线液相上限线令 =5s 得： =90.693由上述关系可作得线降液管液泛线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/h qVVh =895.36 m3/h负荷性能图： 可见，线的位置偏上，所以它对操作的影响很小。放大后的负荷性能图：设计点位于四条线包围的区间中间稍偏下，操作弹性：qVVhmax / qVVhmin2.73所以基本满足要求（程序见附件二）第四章 再沸器的设计一 设计任务与设计条件 1选用立式热虹吸式再沸器 塔顶压力：1.7213MPa 压力降：Nphf=1120.1170.479.8103=0.067363MPa 塔底压力=1.7213+0.0604=1.788663MPa 2再沸器壳程与管程的设计壳程管程温度（）7053.49压力（MPa绝压）0.10131.7817蒸发量：Db= q,mVs =7.01885kg/s 物性数据壳程凝液在温度（100）下的物性数据：潜热：rc=2257kj/kg热导率：c =0.683w/(m*K)粘度：c =0.283mPa*s密度：c =958.4kg/m3管程流体在（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 =441.3kg/m3 蒸气压曲线斜率（t/P）=0.0000266 m2 K/kg 二 估算设备尺寸 热流量： = 1959 143.28w 传热温差： =100-52.5=16.51K 假设传热系数：K=611.4898W/( m2 K) 估算传热面积Ap =194.05 m2 拟用传热管规格为：382.5mm,管长L=3000mm 则传热管数： =542 若将传热管按正三角形排列，按式 NT =3a(a+1)+1;b=2a+1 得：b=10.0498 管心距：t=0.048m 则 壳径： =0.548m 取 D= 800mm L/D=3.75 取 管程进口直径：Di=0.25m 管程出口直径：Do=0.3m 三 传热系数的校核1显热段传热系数K假设传热管出口汽化率 Xe=0.19则循环气量： =36.941kg/+s 计算显热段管内传热膜系数i 传热管内质量流速： di=38-22.5=33mm =0.13149 = 644.kg/( m2 s) 雷诺数： = 299629.5 普朗特数： =2.69058 显热段传热管内表面系数： = 2031.37w/( m2 K) 2) 壳程冷凝传热膜系数计算o 蒸气冷凝的质量流量： = 0.8393kg/s 传热管外单位润湿周边上凝液质量流量： =0.1049kg/(ms) = 1033.929 管外冷凝表面传热系数： = 3630.96w/ (m2 K) 3) 污垢热阻及管壁热阻 沸腾侧：Ri=0.0005117 m2 K/w 冷凝侧：Ro=0.00021 m2 K/w 管壁热阻：Rw=b/w= 0.00005117 m2 K/w 4）显热段传热系数 dm=(di+do)/2= 0.0355m = 589.49w/( m2 K) 2. 蒸发段传热系数KE计算 传热管内釜液的质量流量：Gh=3600 qmws = 2320767.243kg/( 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）泡核沸腾压抑因数：=(E+)/2=0.8 泡核沸腾表面传热系数: =8104.605 w/( m2 K) 3）单独存在为基准的对流表面传热系数 ： = 1080.925w/( m2 K) 1906.9 沸腾表面传热系数：KE 对流沸腾因子 ： = 1.211 两相对流表面传热系数: = 2309.34 w/( m2 K) 沸腾传热膜系数： =2309.34 w/( m2 K) = 614.20 w/( m2 K) 3.显热段及蒸发段长度 =0.0194 LBC = 0.00272L= 0.0583m LCD =L- LBC = 2.9417m 4传热系数 = 613.72 m2 实际需要传热面积： = 193.35 m25传热面积裕度： = 0.48780.150.2所以，传热面积裕度合适，满足要求四 循环流量校核1循环系统推动力：1）当X=Xe/3= 0.063333333时=9.9815 两相流的液相分率： = 0.567 两相流平均密度： = 441.3kg/m3 2）当X=Xe=0.19 = 3.19两相流的液相分率： = 0.3608 两相流平均密度： = 441.3kg/m3根据课程设计表319 得：L=1.02m, 则循环系统的推动力： = 3462.3pa 2循环阻力Pf： 管程进出口阻力P1 进口管内质量流速： = 752.6kg/(m2s) 釜液进口管内流动雷诺数： = 2 649 923 进口管内流体流动摩擦系数： = 0.0150 进口管长度与局部阻力当量长度： =29.30m管程进出口阻力: =11128.34Pa 传热管显热段阻力P2 =310.687kg/(m2s) =144404 =0.0205 = 1.102Pa传热管蒸发段阻力P3 气相流动阻力Pv3 =39.3537kg/(m2s) =18291.167 =0.030 =4.86Pa液相流动阻力PL3GL=G-Gv=271.33 kg/(m2s) = 126112 = 0.0210 = 158.464Pa = 636.412Pa 管内动能变化产生阻力P4 动量变化引起的阻力系数: = 2.58 = 190.48管程出口段阻力P5 气相流动阻力Pv5 = 518.786kg/(m2s) = 99.29kg/(m2s) 管程出口长度与局部阻力的当量长度之和： = 35.0425m = 419571.29 = 0.017776 = 149.47 液相流动阻力PL5 =423.33 kg/(m2s) = 1788698 = 0.01544 =366.28Pa = 23838.65Pa 所以循环阻力：Pf=P1 + P2 + P3 + P4 + P5 = 5794.97 又因PD=6228.4845Pa 所以 =1.0135 第五章 辅助设备设计 一 辅助容器的设计 容器填充系数取：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=3416kg/h 取 停留时间：x为4天，即x=96h 进料罐容积： 907.5m3 圆整后 取V=908 m3 2回流罐（43）质量流量qmLh=3600RqmDs =26444.6kg/h设凝液在回流罐中停留时间为0.25h，填充系数=0.7则回流罐的容积 19.87 m3取V=20 m33塔顶产品罐质量流量qmDh=3600qmDs =2206.903 kg/h；产品在产品罐中停留时间为72h，填充系数=0.7则产品罐的容积 4829.6m3取V=4830m3釜液罐取停留时间为5天，即x=120h质量流量qmWh=3600qmWs =1209.2 kg/h 则釜液罐的容积 462.70 m3取V=463m3二 传热设备 1进料预热器 用90水为热源，出口约为70走壳程 料液由20加热至45，走管程传热温差： 管程液体流率：qmfh=3600 qmfs=3904kg/h 管程液体焓变：H=401kj/kg 传热速率：Q= qmfsH=2989.08401/3600=434.3kw 壳程水焓变：H=125.6kj/kg 壳程水流率：q=11856kg/h 假设传热系数：K=650w/(m2K) 则传热面积： 圆整后取A=12m2 塔顶冷凝器拟用10水为冷却剂，出口温度为30。走壳程。管程温度为43管程流率：qmVs=7.954kg/s取潜热r=302.54kj/kg传热速率：Q= qmVsr=2407kw壳程取焓变：H=125.5kj/kg则壳程流率：qc=Q/H=69068.57kg/h假设传热系数：K=650 w/(m2K)则传热面积： 圆整后 取A=173m2塔顶产品冷却器拟用10水为冷却剂，出口温度为20。走壳程。管程温度由43降至25 管程流率：qmDs = 00.613029kg/s ; 取潜热：r=280kj/kg则传热速率：Q= qmDsr=171.65kw壳程焓变：H=84.0kj/kg则壳程流率：qc=Q/H=7356.73kg/h假设传热系数：K=650 w/(m2K)则传热面积 圆整后 取A=17 m2釜液冷却器拟用10水为冷却剂，出口温度为20。走壳程。管程温度由52.5降到25管程流率：qmWs=055kg/s丙烷液体焓变：H =282kj/kg传热速率：Q= qmVsH =94.72kw壳程取焓变：H=84.0kj/kg则壳程流率：qc=Q/H=4059kg/h假设传热系数：K=650 w/(m2K)则传热面积： 圆整后 取A=9 m2三 泵的设计1进料泵(两台，一用一备)取液体流速：u=0.5m/s液体密度： kg/ m3 qVfs = qmfs / =0.0018 m3/s 取d=65mm液体粘度 取=0.2相对粗糙度：/d=0.003查得：=0.025取管路长度：l=100m 取90度弯管4个，截止阀一个，文氏管流量计1个取则qVLh =6.823m3/h选取泵的型号：AY 扬程：3065m 流量：2.560m3 /s2回流泵（两台，一开一用）取液体流速：u=0.5m/s液体密度： kg/ m3 qVLs = qmLs / =0.0137 m3/s 液体粘度 取=0.2相对粗糙度：/d=0.00107查得：=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第七章 控制方案 精馏塔的控制方案要求从质量指标、产品产量和能量消耗三个方面进行综合考虑。精馏塔最直接的质量指标是产品浓度。由于检测上的困难，难以直接按产品纯度进行控制。最常用的间接质量指标是温度。 将本设计的控制方案列于下表序号位置用途控制参数介质物性L(kg/m3)1FIC-01进料流量控制03000kg/h丙烷、丙稀L=516.32FIC-02回流定量控制01500kg/h丙稀L=469.83PIC-01塔压控制02MPa丙稀V=284HIC-02回流罐液面控制01m丙稀L=4705HIC-01釜液面控制03m丙烷L=4426TIC-01釜温控制4060丙烷L=442附件一 C程序 1#include /* 此头函数请不要删除 */#include #include int i,j,n,N,Nf;double T,pa,fn,fnd,x0,m,s,a,xe,ye,xd,xw,R,Rmin,pb,a1,a2,qmd,qmw,qmf,aa,dl,dv,Flv,Ht,Hf; static double p2=1,1;static double K2=1,1; static double X200; static double Y200;static double qml2;static double qmv2;static double A2=15.7027,15.7260;static double B2=1807.53,1872.46;static double C2=-26.15,-25.16;void pi(int i,double T) pi=exp(Ai-Bi/(T+Ci)*101.3/760; void Ki(int i,double pa) Ki=pi/pa;void sumy(double x0) s=K0*x0+K1*(1-x0); void test(double T,double pa,double x0) fn=0;fnd=1; do T=T-fn/fnd; for(i=0;i0.001); printf(T=%f ,T); printf(pa=%f ,pa);for(i=0;i2;i+) printf(p=%f ,pi);printf(K=%fn,Ki); a=K0/K1; printf(a=%f ,a); main()double f,c,c20,uf,u,A,At,Ad,b,Dc,D,hb,hw,lw,how,bs,bc,bd,Aa,x,r,fi,d0,Ao,uo,uo1,t,ev,co,g,ho,hl,be,ua,Fa,hs,hf,FI,hd,Hd,Hd1,to,ho1,k,d; test(316,1721.3,0.98);a1=a;doa2=a;a=(a1+a2)/2;xd=0.98;xw=0.02;qmd=0.5364;qmw=0.2939;qmf=0.8303;dl=470;dv=28,Ht=0.45;xe=0.65;ye=a*xe/(1+(a-1)*xe);Rmin=(xd-ye)/(ye-xe);R=Rmin*1.6;Y0=xd;for(j=0;j+)Xj=Yj/(a+(1-a)*Yj);Yj+1=R*Xj/(R+1)+xd/(R+1);if(Xj0.65) Nf=(int)(j/0.6)+1;break;for(;j+) Xj=Yj/(a+(1-a)*Yj);Yj+1=(R+1.5238)*Xj/(R+1)-0.5238*xw/(R+1); if(Xj0.0001);a=aa;if(10*(j-1)%6!=0) N=(int)(j-1)/0.6)+1;else N=10*(j-1)/6;qml0=R*qmd;qmv0=(R+1)*qmd;qml1=qml0+qmf;qmv1=qmv0;printf(a=%f ,a);printf(R=%f ,R);printf(Nf=%d ,Nf);printf(N=%dn,N);for(i=0;i2;i+)printf(qml=%f ,qmli); printf(qmv=%f ,qmvi); printf(n);Flv=(qml0/qmv0)*sqrt(dv/dl);printf(Flv=%fn,Flv);f=4.5;c20=0.061;b=0.15; c=c20*pow(f/20),0.2);uf=c*sqrt(dl-dv)/dv);u=uf*0.7;A=qmv0/(dv*u);At=A/(1-b);Dc=sqrt(4*At/3.14159);printf(c=%f ,c);printf(uf=%f ,uf);printf(Dc=%fn,Dc);D=1.8;At=0.25*3.14159*D*D;A=(1-b)*At;Ad=At-A;u=qmv0/(A*dv);m=u/uf;b=Ad/At;printf(At=%f ,At);printf(A=%f ,A);pr