课程设计报告
( 2013—2014年度第一学期)
名 称： 过程参数检测及仪表 题 目： 第二类命题的计算 院 系： 控制与计算机工程学院 班 级： 学 号： 学生姓名： 指导教师： 段泉圣 设计周数： 一周 成 绩：
日期： 2014年 1 月14日
一、课程设计的目的与要求
1. 通过本次课程设计，加深自己对课程中理论基础的理解，掌握设计过程参数检测系统的基本思路、设计手段以及设计方法。
2. 掌握节流变压降流量测量系统的组成，各部分功能以及节流装置的计算。 3. 掌握标准节流装置第二类命题的设计与计算。 4. 能根据实际工程情况设计合理的流量检测系统。
二、设计正文
1. 设计题目
已知：被测流体（水蒸气）及工作压力p 1、工作温度t 1 、最大流量q m max 、常用流量q mch 、最小流量q m min 、管道材料、内径D 20、节流件前后管道及阻流件的情况。给出的限制条件有最小直管段L ≤ L x （实际直管段长度）；压力损失∆P ch ≤∆w y 要求完成以下工作：
（1）选择节流件形式（类型、取压方式）； （2）选择差压计；
（3）计算C 、ε 、β 、d 20； （4）计算所需最小直管段l 0、l 1 、l 2；
（5）计算最大压损∆P max ，并验算∆p max ≤∆p y ； （6）计算所设计的节流装置的测量不确定度；
本题的设计要求如下：
2. 设计中各参数定义
3. 具体设计过程：
(1)工作状态下各介质参数：
（1.1）根据介质工作状态下表压力P 0=0. 1MP ，绝对压力 1=2.1MP , 大气压力P
P a =P 1+P 0=2.1+0.1=2.2MP
再由介质工作温度t1=310℃，查表可得工作状态下水的密度ρ1=8.64118Kg /m ，动力粘度η=20.78⨯10Pa ⋅s ；
（1.2）根据管道材料为3号钢制作的新的无缝钢管和工作温度310℃，查表得管道材料的线膨胀系数λD =13.60⨯10
-6
3
-6
/ C 。
（1.3）过热蒸汽的等熵指数κ=1.3。 (2)计算工作状态下的管道直径： 根据如下公式计算管道直径：
D t = D 20 [1+λD (t 1-20)]
=500× [1+13.60×10-6×(310-20)] = 501.97mm
(3)计算常用流量和最小流量下的雷诺数：
根据如下公式计算
Re Dch =
4q mch 4⨯200000
==6.781⨯106 -3-6
πD t ηπ⨯501.97⨯10⨯20.78⨯10⨯3600
Re D min =
4q min 4⨯80000
==2.713⨯106 -3-6
πD t ηπ⨯501.97⨯10⨯20.78⨯10⨯3600
(4)一次装置的选定：
根据题中所给出要求，一次装置的形式为标准孔板，取压方式为法兰取压 (5)差压计的选择 满刻度流量值q m max ：
*q m max 的确定是根据用户给出的最大流量q m max 向上圆整到系列值中的临近值。
*
q k =a ⨯10n
a =1,1.25,1.6,2,2.5,3.2,4,5,6.3,8
本例中最大流量为250t/h,因此流量计的刻度上限q*m max定为250t/h；
型号：差压计的类型要根据使用条件与测量量精度进行选取。电容式差压变送器稳定性与准确性都较高，因此本题中选择DP1151电容式差压变送器检测差压信号。 对应的∆P max 的确定：节流件形式为标准孔板
∆P max =2.5∆ωy =2.5⨯56=140kP a
圆整要求：(1,1.6,2.5, 4,6) ⨯10
根据圆整要求，另a =1.6，则∆P max =160kPa 。
参考书本中附表的内容，可选用1151DP 电容式差压变送器，其量程范围为
n
（3175~19050）⨯9. 81Pa ，耐静压值为14MP ，精度为±0. 2，差压计的量程调在0~160kPa 。
验算：∆p max /p 1=160⨯10/2.1⨯10=0.07619
(7) 迭代计算β值和d 值
I. 求不变量A:
3
6
q mch 2200⨯1032
∆P ch =(*) ∆P max =() ⨯160=102.4 kPa
q mmax 250⨯103
将流量公式: q mch =
β Dt
中已知数移至等号一
π
22
4
边，形成一变量A ，即：
A =
4
q mch
Dt
=0.211023=
II. 假定初始值C 0=1; ①根据公式X n =
ε0=1，并根据公式逐个计算参数：
A
可得:
C n -1εn -1
X 1=
A C 0ε0
=
0.211023
= 0.211023
1⨯1
2X n 1/4
②根据公式: βn =(可得: ) 2
1+X n
X 121/40.21102321/4
β1=() =() = 0.454396 22
1+X 11+0.211023
3对于过热蒸汽, 可膨胀性系数 ○
ε=1-(0.41+0.35β4) ∆P /κρ1= 0.984061
4根据公式： ○
C =0.5959+0.0312β2.1-0.1840β8+0.0029β2.5(106/Re Dch ) 0.75+2.2860β4/D(1-β4)-0.8560β3/D = 0.601658
○5根据公式：δn =A -X n C n εn
就可通过前面确定的参数X n , C , n εn 计算出每一次的压力损失δn III. 按照上述公式进行迭代计算:
从n =3起进行迭代计算，X n 用快速弦截法公式:
X n =X n -1-δn -1
X n -1-X n -2
δn -1-δn -2
进行迭代计算，精度判别公式为E n =δn /A ，假设判别条件为
E n ≤1⨯10-7(n =1, 2, ⋅⋅⋅⋅⋅⋅) ，迭代计算结果如下表所示
表一:程序迭代结果
当n =4时，求得的E 小于预定精度，由此可得:
β=β4=0.57870465
C =C 4=0.60389936
IV. 求d t , 根据公式
d t =D t β=501.97⨯0.57870465=290.492047mm
(8)验算压力损失：
根据公式:
∆ωmax =(1-β1.9) ∆p ch =(1-0.578704651.9) ⨯102400=66.178429Kpa >∆ωy 压力
损失验算不合格;
差压法流量测量中不可避免地存在一定的压损，下图中的P2-P3为流体通过节流装置时产生的压损：
根据GB2624-93中标准节流装置的设计原理可知，管道流经标准孔板的压力损失为：
∆δ=(1-β1.9) ∆p ch
式中：β—节流件内径与管道内径之比d/D； △P —流量差压；
从上式中不难看出，在一定的流体情况下，差压△P 越大，压损越大; β值越小，压损越大。
由于差压△P ，直径比β与节流件的给定参数、约束条件相关，例如管道中的流量范围，工作压力与工作温度等。具体的分析过程如下：
○1从流量q 的角度考虑：管道内流束的流量q m 直接影响着压差∆p ch ，适当减小管道内正常流量q mch ，可减小∆p ch ，从而减小压损δp
2从工作压力P ○
1的角度考虑：若其他条件不变，提高管道工作压力
P 1，则过热蒸汽的密度ρ增大，
A 值减小，X 减小，直径比β
值减小，压损会相应的增加。因此为了减小压损以符合要求，可适当降低管道工作压力。
3从工作温度t 的角度考虑：其他条件不变，适当提高管道工作温度t ，流束密度ρ减○
小，A 值增加，X 增加，直径比β值增加，压损会相应的降低。因此为了减小压损以符合要求，可适当提高管道工作温度。 (9)验算流量：
q m ==
d t 0.983148⨯290.4920472=2000069.1303kg /h
(10)确定d 20：
选用的孔板材料为1C r 18N i 9T i ，温度t=310℃，查表可得节流件线性膨胀系数
λd =17.50⨯10-6/℃
根据公式：d 20=
d t 290.492047
==289.025244mm -6
1+λd (t -20) ⎡⎣1+17.50⨯10(310-20) ⎤⎦
(11)确定安装位置：
根据β=0.578704
l 0=18D t =9.03546m l 1=26D t =13.05122m l 2=7D t =3.51379m
标准孔板对管段粗糙度的限制可由表7-4查出，对于β=0.578704
Ks ≤D t ⨯4.2⨯10-4=501.97⨯4.2⨯10-4=0.211mm
查表可得新的无缝钢管的绝对粗糙度K s =0.05mm ~0.100mm
对于标准孔板，由β=0.578704
β≤0.6时，δC /C =±0.6%
2流束膨胀系数的不确定度： ○
δε∆P ch 102.4⨯103
=±4%=±4⨯%=±0.195048%，置信概率为95%。 6εp 12.1⨯10
3节流件孔径的不确定度： ○
δd
d
=±0. 05%，置信概率为95%
4管道直径不确定度： ○
δD
D
=±0. 1%
5密度的不确定度： ○
假设温度，压力的测量不确定度为±1%，则
δρ
1
ρ1
=±1%
6差压的不确定度： ○
选定的差压变送器为0.2级，则差压测量的不确定度：
δ∆p
21
=⨯⨯0.2%⨯160000=0.21% ∆p 3102400
7总的不确定度： ○
1
2
δq
⎡⎛δ⎫⎛δ⎫⎛2β4⎫⎛δ⎫⎛2⎫⎛δ⎫1⎛δ∆p ⎫1⎛δρ⎫⎤
d D
=±⎢ C ⎪+ ε⎪+ ⎪+ ⎪⎥⎪+ ⎪+ 4⎪ 4⎪ q m C ε1-βD 1-βd 4∆p ⎢⎭⎝⎭⎭⎝⎭⎝⎝⎭4⎝ρ⎭⎥⎣⎝⎭⎝⎭⎝⎦
2
2
2
2
m
2222
⎡⎛2⨯0.5787044⎫211⎤222⎛⎫22=±（0.6%）+（0.195048%）+0.1%+0.05%+0. 21%+%⎥) )()⎢ 4⎪(4⎪(44⎦⎝1-0.578704⎭⎝1-0.578704⎭⎣
=± 0.973255%
0.5
4. 结论总结
（1）选择节流件形式：标准孔板 法兰取压
（2）选择差压计：DP1151电容式差压变送器 精度0.2级 （3）计算C 、ε 、β 、d 20
C = 0.60389936
β= 0.57870465
d 20=289.025244mm
ε=0.98314882
（4）计算所需最小直管段l 0、l 1 、l 2 :
l 0=18D t =9.03546m l 1=26D t =13.05122m
l 2=7D t =3.51379m
（5）计算最大压损∆P max ，并验算∆P max ≤∆w y :
∆P max =160kPa 验算压力损失不合格 ，经分析通过以下方式改进:适当减小管道内正常流量q mch 、适当降低管道工作压力、适当提高管道工作温度
（6）计算所设计的节流装置的测量不确定度：
δq m
q m =± 0.973255%
5. 用Matlab 实现的源程序：
迭代：
clc;clear
format long
a=0.211023;c=1;e=1;re=6.7481*10^6;D=501.97;m=0.0487619;
x=a/c/e
b=((x*x)/(1+x*x))^(0.25)
c=0.5959+0.0312*b^2.1-0.1840*b^8+0.0029*b^2.5*(10^6/re)^0.75+2.2860*b^4/(D*(1-b^4))-0.8560*b^3/D
e=1-(0.41+0.35*b^4)*m/1.30
q=a-x*c*e
en=q/a
x1=x,q1=q;
x=a/c/e
b=((x*x)/(1+x*x))^(0.25)
c=0.5959+0.0312*b^2.1-0.1840*b^8+0.0029*b^2.5*(10^6/re)^0.75+2.2860*b^4/(D*(1-b^4))-0.8560*b^3/D
e=1-(0.41+0.35*b^4)*m/1.30
q=a-x*c*e
en=q/a
x2=x;q2=q;
n=2;
while (abs(en)>1*10^(-7))
x=x2-q2*((x2-x1)/(q2-q1))
x1=x2;
x2=x;
b=((x*x)/(1+x*x))^(0.25)
c=0.5959+0.0312*b^2.1-0.1840*b^8+0.0029*b^2.5*(10^6/re)^0.75+2.2860*b^4/(D*(1-b^4))-0.8560*b^3/D
e=1-(0.41+0.35*b^4)*m/1.30
q=a-x*c*e
q1=q2;
q2=q;
en=q/a
n=n+1;
end
n
不确定度：
format long
xl=(0.8)^2+(0.195048)^2+2*0.5787^4/(1-0.5787^4)*(0.1)^2+2/(1-0.5787^4)*(0.05)^2+0.25*(0.21)^2+0.25
x2=sqrt(xl)
6. 设计流程图
三. 课设心得
本次课程设计是本学期最后一门课程设计，经过对第二类命题设计的计算，我有以下体会：
在学习过程中要明确对变量的定义和把握。一个设计命题所涵盖的变量就需要前后翻找，精确地理解其意义；在计算过程中，需要反复查找附录中的表格，有的可以线性内插，有的不可以线性内插。在计算过程中，随时需要上网查找相关国标，以求结果准确性更高。
当计算的压力损失验算不合格时，通过思考找到了分析及改进的方式。在验算不合格的初期，总是在疑惑算错了，经过反复验算，终于体会到应该如何分析。
在计算过程中，尤其是迭代计算部分，用到Matlab 软件编程，在计算不确定部分，也可以通过软件计算。启发了我今后多应用相关软件减少计算量，获得更准确的结果。
这次课程设计让我对流量部分的基础知识有了一次全面的回顾，包括各种一次装置的取压方式及其对应各种物理量的经验公式都有所不同，需要我们设计的时候仔细比对，选择正确的公式进行计算。在学习理论课的过程中，对于流量的掌握还有很多不完美之处，很高兴通过这次课程设计有了一个查缺补漏的机会。
四. 参考文献
[1] 常太华，苏杰. 过程参数检测及仪表. 北京：中国电力出版社
课程设计报告
( 2013—2014年度第一学期)
名 称： 过程参数检测及仪表 题 目： 第二类命题的计算 院 系： 控制与计算机工程学院 班 级： 学 号： 学生姓名： 指导教师： 段泉圣 设计周数： 一周 成 绩：
日期： 2014年 1 月14日
一、课程设计的目的与要求
1. 通过本次课程设计，加深自己对课程中理论基础的理解，掌握设计过程参数检测系统的基本思路、设计手段以及设计方法。
2. 掌握节流变压降流量测量系统的组成，各部分功能以及节流装置的计算。 3. 掌握标准节流装置第二类命题的设计与计算。 4. 能根据实际工程情况设计合理的流量检测系统。
二、设计正文
1. 设计题目
已知：被测流体（水蒸气）及工作压力p 1、工作温度t 1 、最大流量q m max 、常用流量q mch 、最小流量q m min 、管道材料、内径D 20、节流件前后管道及阻流件的情况。给出的限制条件有最小直管段L ≤ L x （实际直管段长度）；压力损失∆P ch ≤∆w y 要求完成以下工作：
（1）选择节流件形式（类型、取压方式）； （2）选择差压计；
（3）计算C 、ε 、β 、d 20； （4）计算所需最小直管段l 0、l 1 、l 2；
（5）计算最大压损∆P max ，并验算∆p max ≤∆p y ； （6）计算所设计的节流装置的测量不确定度；
本题的设计要求如下：
2. 设计中各参数定义
3. 具体设计过程：
(1)工作状态下各介质参数：
（1.1）根据介质工作状态下表压力P 0=0. 1MP ，绝对压力 1=2.1MP , 大气压力P
P a =P 1+P 0=2.1+0.1=2.2MP
再由介质工作温度t1=310℃，查表可得工作状态下水的密度ρ1=8.64118Kg /m ，动力粘度η=20.78⨯10Pa ⋅s ；
（1.2）根据管道材料为3号钢制作的新的无缝钢管和工作温度310℃，查表得管道材料的线膨胀系数λD =13.60⨯10
-6
3
-6
/ C 。
（1.3）过热蒸汽的等熵指数κ=1.3。 (2)计算工作状态下的管道直径： 根据如下公式计算管道直径：
D t = D 20 [1+λD (t 1-20)]
=500× [1+13.60×10-6×(310-20)] = 501.97mm
(3)计算常用流量和最小流量下的雷诺数：
根据如下公式计算
Re Dch =
4q mch 4⨯200000
==6.781⨯106 -3-6
πD t ηπ⨯501.97⨯10⨯20.78⨯10⨯3600
Re D min =
4q min 4⨯80000
==2.713⨯106 -3-6
πD t ηπ⨯501.97⨯10⨯20.78⨯10⨯3600
(4)一次装置的选定：
根据题中所给出要求，一次装置的形式为标准孔板，取压方式为法兰取压 (5)差压计的选择 满刻度流量值q m max ：
*q m max 的确定是根据用户给出的最大流量q m max 向上圆整到系列值中的临近值。
*
q k =a ⨯10n
a =1,1.25,1.6,2,2.5,3.2,4,5,6.3,8
本例中最大流量为250t/h,因此流量计的刻度上限q*m max定为250t/h；
型号：差压计的类型要根据使用条件与测量量精度进行选取。电容式差压变送器稳定性与准确性都较高，因此本题中选择DP1151电容式差压变送器检测差压信号。 对应的∆P max 的确定：节流件形式为标准孔板
∆P max =2.5∆ωy =2.5⨯56=140kP a
圆整要求：(1,1.6,2.5, 4,6) ⨯10
根据圆整要求，另a =1.6，则∆P max =160kPa 。
参考书本中附表的内容，可选用1151DP 电容式差压变送器，其量程范围为
n
（3175~19050）⨯9. 81Pa ，耐静压值为14MP ，精度为±0. 2，差压计的量程调在0~160kPa 。
验算：∆p max /p 1=160⨯10/2.1⨯10=0.07619
(7) 迭代计算β值和d 值
I. 求不变量A:
3
6
q mch 2200⨯1032
∆P ch =(*) ∆P max =() ⨯160=102.4 kPa
q mmax 250⨯103
将流量公式: q mch =
β Dt
中已知数移至等号一
π
22
4
边，形成一变量A ，即：
A =
4
q mch
Dt
=0.211023=
II. 假定初始值C 0=1; ①根据公式X n =
ε0=1，并根据公式逐个计算参数：
A
可得:
C n -1εn -1
X 1=
A C 0ε0
=
0.211023
= 0.211023
1⨯1
2X n 1/4
②根据公式: βn =(可得: ) 2
1+X n
X 121/40.21102321/4
β1=() =() = 0.454396 22
1+X 11+0.211023
3对于过热蒸汽, 可膨胀性系数 ○
ε=1-(0.41+0.35β4) ∆P /κρ1= 0.984061
4根据公式： ○
C =0.5959+0.0312β2.1-0.1840β8+0.0029β2.5(106/Re Dch ) 0.75+2.2860β4/D(1-β4)-0.8560β3/D = 0.601658
○5根据公式：δn =A -X n C n εn
就可通过前面确定的参数X n , C , n εn 计算出每一次的压力损失δn III. 按照上述公式进行迭代计算:
从n =3起进行迭代计算，X n 用快速弦截法公式:
X n =X n -1-δn -1
X n -1-X n -2
δn -1-δn -2
进行迭代计算，精度判别公式为E n =δn /A ，假设判别条件为
E n ≤1⨯10-7(n =1, 2, ⋅⋅⋅⋅⋅⋅) ，迭代计算结果如下表所示
表一:程序迭代结果
当n =4时，求得的E 小于预定精度，由此可得:
β=β4=0.57870465
C =C 4=0.60389936
IV. 求d t , 根据公式
d t =D t β=501.97⨯0.57870465=290.492047mm
(8)验算压力损失：
根据公式:
∆ωmax =(1-β1.9) ∆p ch =(1-0.578704651.9) ⨯102400=66.178429Kpa >∆ωy 压力
损失验算不合格;
差压法流量测量中不可避免地存在一定的压损，下图中的P2-P3为流体通过节流装置时产生的压损：
根据GB2624-93中标准节流装置的设计原理可知，管道流经标准孔板的压力损失为：
∆δ=(1-β1.9) ∆p ch
式中：β—节流件内径与管道内径之比d/D； △P —流量差压；
从上式中不难看出，在一定的流体情况下，差压△P 越大，压损越大; β值越小，压损越大。
由于差压△P ，直径比β与节流件的给定参数、约束条件相关，例如管道中的流量范围，工作压力与工作温度等。具体的分析过程如下：
○1从流量q 的角度考虑：管道内流束的流量q m 直接影响着压差∆p ch ，适当减小管道内正常流量q mch ，可减小∆p ch ，从而减小压损δp
2从工作压力P ○
1的角度考虑：若其他条件不变，提高管道工作压力
P 1，则过热蒸汽的密度ρ增大，
A 值减小，X 减小，直径比β
值减小，压损会相应的增加。因此为了减小压损以符合要求，可适当降低管道工作压力。
3从工作温度t 的角度考虑：其他条件不变，适当提高管道工作温度t ，流束密度ρ减○
小，A 值增加，X 增加，直径比β值增加，压损会相应的降低。因此为了减小压损以符合要求，可适当提高管道工作温度。 (9)验算流量：
q m ==
d t 0.983148⨯290.4920472=2000069.1303kg /h
(10)确定d 20：
选用的孔板材料为1C r 18N i 9T i ，温度t=310℃，查表可得节流件线性膨胀系数
λd =17.50⨯10-6/℃
根据公式：d 20=
d t 290.492047
==289.025244mm -6
1+λd (t -20) ⎡⎣1+17.50⨯10(310-20) ⎤⎦
(11)确定安装位置：
根据β=0.578704
l 0=18D t =9.03546m l 1=26D t =13.05122m l 2=7D t =3.51379m
标准孔板对管段粗糙度的限制可由表7-4查出，对于β=0.578704
Ks ≤D t ⨯4.2⨯10-4=501.97⨯4.2⨯10-4=0.211mm
查表可得新的无缝钢管的绝对粗糙度K s =0.05mm ~0.100mm
对于标准孔板，由β=0.578704
β≤0.6时，δC /C =±0.6%
2流束膨胀系数的不确定度： ○
δε∆P ch 102.4⨯103
=±4%=±4⨯%=±0.195048%，置信概率为95%。 6εp 12.1⨯10
3节流件孔径的不确定度： ○
δd
d
=±0. 05%，置信概率为95%
4管道直径不确定度： ○
δD
D
=±0. 1%
5密度的不确定度： ○
假设温度，压力的测量不确定度为±1%，则
δρ
1
ρ1
=±1%
6差压的不确定度： ○
选定的差压变送器为0.2级，则差压测量的不确定度：
δ∆p
21
=⨯⨯0.2%⨯160000=0.21% ∆p 3102400
7总的不确定度： ○
1
2
δq
⎡⎛δ⎫⎛δ⎫⎛2β4⎫⎛δ⎫⎛2⎫⎛δ⎫1⎛δ∆p ⎫1⎛δρ⎫⎤
d D
=±⎢ C ⎪+ ε⎪+ ⎪+ ⎪⎥⎪+ ⎪+ 4⎪ 4⎪ q m C ε1-βD 1-βd 4∆p ⎢⎭⎝⎭⎭⎝⎭⎝⎝⎭4⎝ρ⎭⎥⎣⎝⎭⎝⎭⎝⎦
2
2
2
2
m
2222
⎡⎛2⨯0.5787044⎫211⎤222⎛⎫22=±（0.6%）+（0.195048%）+0.1%+0.05%+0. 21%+%⎥) )()⎢ 4⎪(4⎪(44⎦⎝1-0.578704⎭⎝1-0.578704⎭⎣
=± 0.973255%
0.5
4. 结论总结
（1）选择节流件形式：标准孔板 法兰取压
（2）选择差压计：DP1151电容式差压变送器 精度0.2级 （3）计算C 、ε 、β 、d 20
C = 0.60389936
β= 0.57870465
d 20=289.025244mm
ε=0.98314882
（4）计算所需最小直管段l 0、l 1 、l 2 :
l 0=18D t =9.03546m l 1=26D t =13.05122m
l 2=7D t =3.51379m
（5）计算最大压损∆P max ，并验算∆P max ≤∆w y :
∆P max =160kPa 验算压力损失不合格 ，经分析通过以下方式改进:适当减小管道内正常流量q mch 、适当降低管道工作压力、适当提高管道工作温度
（6）计算所设计的节流装置的测量不确定度：
δq m
q m =± 0.973255%
5. 用Matlab 实现的源程序：
迭代：
clc;clear
format long
a=0.211023;c=1;e=1;re=6.7481*10^6;D=501.97;m=0.0487619;
x=a/c/e
b=((x*x)/(1+x*x))^(0.25)
c=0.5959+0.0312*b^2.1-0.1840*b^8+0.0029*b^2.5*(10^6/re)^0.75+2.2860*b^4/(D*(1-b^4))-0.8560*b^3/D
e=1-(0.41+0.35*b^4)*m/1.30
q=a-x*c*e
en=q/a
x1=x,q1=q;
x=a/c/e
b=((x*x)/(1+x*x))^(0.25)
c=0.5959+0.0312*b^2.1-0.1840*b^8+0.0029*b^2.5*(10^6/re)^0.75+2.2860*b^4/(D*(1-b^4))-0.8560*b^3/D
e=1-(0.41+0.35*b^4)*m/1.30
q=a-x*c*e
en=q/a
x2=x;q2=q;
n=2;
while (abs(en)>1*10^(-7))
x=x2-q2*((x2-x1)/(q2-q1))
x1=x2;
x2=x;
b=((x*x)/(1+x*x))^(0.25)
c=0.5959+0.0312*b^2.1-0.1840*b^8+0.0029*b^2.5*(10^6/re)^0.75+2.2860*b^4/(D*(1-b^4))-0.8560*b^3/D
e=1-(0.41+0.35*b^4)*m/1.30
q=a-x*c*e
q1=q2;
q2=q;
en=q/a
n=n+1;
end
n
不确定度：
format long
xl=(0.8)^2+(0.195048)^2+2*0.5787^4/(1-0.5787^4)*(0.1)^2+2/(1-0.5787^4)*(0.05)^2+0.25*(0.21)^2+0.25
x2=sqrt(xl)
6. 设计流程图
三. 课设心得
本次课程设计是本学期最后一门课程设计，经过对第二类命题设计的计算，我有以下体会：
在学习过程中要明确对变量的定义和把握。一个设计命题所涵盖的变量就需要前后翻找，精确地理解其意义；在计算过程中，需要反复查找附录中的表格，有的可以线性内插，有的不可以线性内插。在计算过程中，随时需要上网查找相关国标，以求结果准确性更高。
当计算的压力损失验算不合格时，通过思考找到了分析及改进的方式。在验算不合格的初期，总是在疑惑算错了，经过反复验算，终于体会到应该如何分析。
在计算过程中，尤其是迭代计算部分，用到Matlab 软件编程，在计算不确定部分，也可以通过软件计算。启发了我今后多应用相关软件减少计算量，获得更准确的结果。
这次课程设计让我对流量部分的基础知识有了一次全面的回顾，包括各种一次装置的取压方式及其对应各种物理量的经验公式都有所不同，需要我们设计的时候仔细比对，选择正确的公式进行计算。在学习理论课的过程中，对于流量的掌握还有很多不完美之处，很高兴通过这次课程设计有了一个查缺补漏的机会。
四. 参考文献
[1] 常太华，苏杰. 过程参数检测及仪表. 北京：中国电力出版社