孔板流量计孔径的计算程序.doc

已知条件：质量流量Qmax,Qcom,Qmin（kg/h）；操作密度M （kg/m3）； 动力黏度N （mPa？s）； 20情况下工艺管道内径D20（m）； 正常操作压力P MPa（G）；正常操作温度T （）；管道材质的线膨胀系数rD； 孔板材质的线膨胀系数rd ；通过计算后所得的数据：正常操作温度T情况下工艺管道内径D（m）；雷诺数ReD；计算差压上限值CPcom 和CPmax；计算A22 ；计算孔板流量计的流出系数C2在知道了相应的工艺条件后，下面我就利用Turbo C 2.0编程软件编写孔板流量计孔径的计算程序。首先，在开始前，我还要先介绍一下Turbo C 2.0相关的几个基本函数，然后通过应用这些函数来编写孔板流量计孔径的计算程序。1 数学函数 包含在如表11所示：（注：使用时在源文件前加命令行：# include）、2 由上面的已知条件和Turbo C 2.0的数学函数，通过Turbo C 2.0编程软件编写孔板流量计孔径的计算程序如下：#includemain（） int n;const float pai=3.141592653289;const float b0=0.500000;const float C0=0.606000;float Qmax,Qcom,Qmin,M,N,D20,rD,rd,T=0.0,b,D,ReD,CPmax,CPcom,A2,L,C,C1,b1,b2,d,d20;printf（Qmax,Qcom,Qmin,M,N,D20,rD,rd,T=n）;scanf（%f%f%f%f%e%f%e%e%f,&Qmax,&Qcom,&Qmin,&M,&N,&D20,&rD,&rd,&T）;/输入已知的工艺参数/D=D20（1+rD（T-20）;ReD=4Qcom/（3600paiND）;CPcom=pow（4Qcomsqrt（1-pow（b0,4）/（3600paib0b0DDC0）,2）1/（2M）;CPmax=CPcompow（Qmax/Qcom）,2）;A2=4Qcom/（3600paiDDsqrt（225e3M）;printf（CPcom=%6.10fn,CPcom）;printf（CPmax=%6.10fn,CPmax）;if（ReD2e5） /利用Turbo C的if判断语句对雷诺数进行判断计算/b=pow（1+pow（0.60/A2）,2）,-0.25）;if（ReD1.0e-10;n+） /利用Turbo C的for循环语句进行迭代计算/b=b1;C=0.5961+0.0261pow（b,2）-0.216pow（b,8）+0.000521pow（pow（10,6）b）/ReD）,0.7）+（0.0188+0.0063pow（19000b/ReD）,0.8）pow（b,3.5）pow（pow（10,6）/ReD）,0.3）+（0.043+0.080exp（-10L）-0.123exp（-7L）（1-0.11pow（19000b/ReD）,0.8）（pow（b,4）/（1-pow（b,4）-0.031（2L/（1-b）-0.8pow（2L/（1-b）,1.1）pow（b,1.3）;b1=pow（1+pow（C/A2）,2）,-0.25）;C1=0.5961+0.0261pow（b1,2）-0.216pow（b1,8）+0.000521pow（pow（10,6）b1）/ReD）,0.7）+（0.0188+0.0063pow（19000b1/ReD）,0.8）pow（b1,3.5）pow（pow（10,6）/ReD）,0.3）+（0.043+0.080exp（-10L）-0.123exp（-7L）（1-0.11pow（19000b1/ReD）,0.8）（pow（b1,4）/（1-pow（b1,4）-0.031（2L/（1-b1）-0.8pow（2L/（1-b1）,1.1）pow（b1,1.3）;b2=pow（1+pow（C1/A2）,2）,-0.25）;printf（C1=%6.10fn,C1）; /输出计算孔板流量计的流出系数C1/printf（b2=%6.10fnn,b2）; /输出计算孔板流量计的直径比b2/d=Db2;d20=d/（1+rd（T-20）;printf（d20=%6.10fnn,d20）; /输出计算孔板流量计在20情况下孔板节流孔直径d20/以上就是计算孔板流量计孔径的计算程序。通过此程序可以很快的计算出孔板流量计在20情况下孔板节流孔孔径的大小。五 应 用实例一： 被测流体 ：水； 流体流量：最大流量Qmax500t/h;正常流量Qcom400t/h;最小流量Qmin200t/h 正常操作压力：P14.6 MPa（G）； 正常操作温度：T220 ； 20情况下工艺管道内径：D20233； 管道材质的线膨胀系数：rD11.16e-6/？； 孔板材质的线膨胀系数:rd=16.60e-6 /？； 操作密度：M850.9082 （kg/m3）； 动力黏度：N0.124 （mPa？s） ；按下Ctrl+F9运行计算程序：（注意单位换算）要求输入已知的工艺参数：Qmax, Qcom, Qmin, M, N, D20, rD, rd, T =500000 400000 200000 850.9082 124e-6 .233 11.16e-6 16.60e-6 233（回车）经过程序的运算得到如下的结果：（按下Alt+F5查看结果）CPcom=166259.4432591512CPmax=257980.3725675625C1=0.6023605334b2=0.5024122240d20=0.116966338实例二： 测流体 ：汽油； 流体流量：最大流量Qmax40000 kg/h;正常流量Qcom31250 kg/h;最小流量Qmin20000 kg/h 正常操作压力：P0.80 MPa（G）； 正常操作温度：T40 ； 20情况下工艺管道内径：D20102； 管道材质的线膨胀系数：rD11.16e-6/？； 孔板材质的线膨胀系数:rd=16.60e-6 /？； 操作密度：M720 （kg/m3）； 动力黏度：N0.3000 （mPa？s） ；按下Ctrl+F9运行计算程序：（注意单位换算）要求输入已知的工艺参数：Qmax, Qcom, Qmin, M, N, D20, rD, rd, T =40000 31250 20000 720 300e-6 .102 11.16e-6 16.60e-6 40（回车）经过程序的运算得到如下的结果：（按下Alt+F5查看结果）CPcom=31982.3496093750CPmax=52399.8828125000C1=0.6041094065b2=0.5007327199d20=0.051069240实例三： 测流体 ：碱液； 流体流量：最大流量Qmax6000 kg/h;正常流量Qcom5000 kg/h;最小流量Qmin3000 kg/h 正常操作压力：P0.65 MPa（G）； 正常操作温度：T40 ； 20情况下工艺管道内径：D2052； 管道材质的线膨胀系数：rD11.16e-6/？； 孔板材质的线膨胀系数:rd=16.60e-6 /？； 操作密度：M1112 （kg/m3）； 动力黏度：N2.0000 （mPa？s） ；按下Ctrl+F9运行计算程序：（注意单位换算）要求输入已知的工艺参数：Qmax, Qcom, Qmin, M, N, D20, rD, rd, T =6000 5000 23000 1112 2000e-6 .052 11.16e-6 16.60e-6 40（回车）经过程序的运算得到如下的结果：（按下Alt+F5查看结果）CPcom=7848.1137695312CPmax=11301.2841796875C1=0.6153961420b