过程参数检测及仪表课程设计-第二类标准节流装置设计命题

最新推荐文章于 2023-02-11 11:38:19 发布

ncepu_Chen

最新推荐文章于 2023-02-11 11:38:19 发布

阅读量4.2k

点赞数 11

分类专栏： # 复习笔记文章标签：华北电力大学过程参数检测及仪表课程设计测控技术与仪器段泉圣

本文链接：https://blog.csdn.net/ncepu_Chen/article/details/94203159

版权

复习笔记专栏收录该内容

12 篇文章 16 订阅

订阅专栏

pdf版本的下载地址: 过程参数检测及仪表课程设计-第二类标准节流装置设计命题（访问密码：3834）

设计任务

序号	项目	符号	单位	参数
1	被测介质名称			过热蒸汽
2	常用流量	$q_{m \, ch}$	t/h	600
3	最小流量	$q_{m\, min}$	t/h	400
4	最大流量	$q_{m \, max}$	t/h	800
5	工作压力	$P_1$	MPa	13.2
6	工作温度	$t_1$	°C	540
7	管道内径	$D_{20}$	mm	288
8	管道材料			20号钢
9	正常流量下的允许压力损失	$\delta p_y$	kPa	50

在这里插入图片描述

设计过程

求工作状态下各介质参数.介质的工作压力(绝对压力) $P_a = P + P_b = 13.3 MPa$ . 由介质的工作温度和工作压力查书后附录中相应表格得到工作状态下过热蒸汽的粘度 $\eta=30.83\times10^{-6} Pa \cdot s$ ,密度 $\rho_1 = 40.326 kg/m^3$ ,管道的线膨胀系数 $\lambda_D = 14.38 \times 10^{-6} / ^\circ C$ ,过热蒸汽的等熵指数取 $\kappa = 1.3$ .
求工作状态下的管道直径
$\begin{aligned} D_t & = D_{20} [1 + \lambda_D (t-20)] \\ & = 288 [1 + 14.38 \times 10^{-6} (540-20)] = 290.15 mm \end{aligned}$
计算常用流量下的雷诺数 $R_{eD}$
$R_{eD} = 0.354 \frac{q_m}{D_t \eta} = 0.354 \times \frac{600000}{290.15 \times 30.83 \times 10^{-6}} = 2.37 \times 10^7$
选取差压计,从用户对压力损失的要求出发,取
$\Delta P = 3 \delta_{p_y} = 3 \times 50 = 150 kPa$
选取最近的差压系列值,取 $\Delta P_{max} = 160kPa$ ,按照题目选择DP/HP1151电容式差压变送器.

本例中流体的最大流量为800t/h,所以流量计的流量刻度上限 $q_{mmax}^*$ 定为800t/h.
计算常用流量下的差压值
$\Delta P_{ch} = \left(\frac{q_{mch}}{q_{mmax}^{*}}\right)^2 \Delta P_{max} = \left(\frac{600}{800}\right)^2 160 = 90 kPa$

迭代计算 $\beta$ 值和 $d$ 值

求不变量A,得
$\frac{q_m}{0.004 D_t^2 \sqrt{\rho_1 \Delta P_{ch}}} = \frac{600000}{0.004 \times 290.15^2 \sqrt{40.326 \times 90000}} = 0.9352598$
设 $C_0=1$ , $\varepsilon_0=1$

根据公式
$\begin{aligned} & X_n = \frac{A}{C_{n-1} \varepsilon_{n-1}} \\ & \beta_n = (\frac{X_n^2}{1+X_n^2})^{\frac{1}{4}} \\ & \varepsilon = 1-(0.41+0.35 \beta^4) \frac{\Delta P}{\kappa P_1} \\ & C = 0.5959 + 0.0312 \beta^{2.1} - 0.1840 \beta^{8} +0.029 \beta^{2.5} (\frac{10^6}{R_{eD}})^{0.75} + 2.2860\frac{\beta^4}{D(1-\beta^4)} - 0.8560 \frac{\beta^3}{D} \\ & \delta_n = A - X_n C_n \varepsilon_n \end{aligned}$
进行迭代运算,从 $n = 3$ 起, $X_n$ 可用快速弦截法公式
$X_n= X_{n-1} - \delta_{n-1} \frac{X_{n-1} - X_{n-2}}{\delta_{n-1} - \delta_{n-2}}$
进行迭代运算,精度判别公式为 $E_n = \frac{\delta_n}{A}$ ,假设判断条件为 $|E_n| \le 5 \times 10^{-10}$

根据上述原理写出代码如下

% 不变量
A = 0.9352598;		% 设弦法不变量
Kappa = 1.3;		% 等熵指数
ReD = 2.37e7;		% 雷诺数
D = 290.15;			% 管道直径
P_1 = 13.2e6;		% 工作压力
Delta_P = 150e3;	% 最大压差
% 迭代变量
n = 1;				% 循环次数
X = zeros();		% 设弦法过程变量
Beta = zeros();		% 直径比 
C = zeros();		% 记录流出系数
Epsilon = zeros();	% 膨胀系数
Delta = zeros();	% 差值
E = zeros();		% 记录误差
% 设置初始流出系数和膨胀系数
C(1) = 1;
Epsilon(1) = 1;

% 循环体
while (n==1 || E(n)>5e-10)	% 误差小于5e-10时跳出循环
	% 迭代次数增加
	n = n+1;

	% 计算中间变量
	X(n) = A/(C(n-1)*Epsilon(n-1));						% 更新设弦法中间变量
	
	Beta(n) = (X(n)^2/(X(n)^2+1))^(0.25);						% 更新直径比

	% 标准喷嘴膨胀系数
	% quotient1 = Kappa*(Tau^(2/Kappa))/(Kappa-1);				% 根号下第一项
	% quotient2 = (1-Beta(n))/(1-(Beta(n)^4)*(Tau^(2/Kappa)));	% 根号下第二项
	% quotient3 = (1-Tau^((Kappa-1)/Kappa))/(1-Tau);			% 根号下第三项
	% Epsilon(n) = sqrt(quotient1*quotient2*quotient3);	% 更新膨胀系数
	% 标准孔板膨胀系数
	Epsilon(n) = 1 - (0.41+0.35*(Beta(n)^4))*(Delta_P)/(Kappa*P_1);% 更新膨胀系数

	% 标准孔板流出系数
	product1 = 0.5959;											% 加法第一项
	product2 = 0.0312*(Beta(n)^2.1);							% 加法第二项
	product3 = -0.1840*(Beta(n)^8);								% 加法第三项
	product4 = 0.029*(Beta(n)^2.5)*((1e6/ReD)^0.75);			% 加法第四项
	product5 = (Beta(n)^4)/(D*(1-(Beta(n)^4)));					% 加法第五项
	product6 = -0.8560*(Beta(n)^3)/D;							% 加法第六项
	C(n) = product1+product2+product3+product4+product5+product6;	% 更新流出系数

	Delta(n) = A - X(n)*C(n)*Epsilon(n);					% 更新差值

	E(n) = abs(Delta(n)/A);									% 计算误差 
    
end

result = [[0:(n-1)]; X; Beta; C; Epsilon; Delta; E];
%result = result(:,2:n);

得到结果如下

n	1	2	3	4	…	9	10
$X$	0.93526	1.611417	1.723145	1.734971	…	1.736195	1.736195
$\beta$	0.82648	0.921784	0.930003	0.930801	…	0.930882	0.930882
$C$	0.582146	0.544656	0.54097	0.540626	…	0.540591	0.540591
$\varepsilon$	0.996993	0.996524	0.996476	0.996472	…	0.996471	0.996471
$\delta$	0.392439	0.060642	0.006374	0.000599	…	4.08E-09	3.77E-10
$E$	0.419604	0.06484	0.006816	0.00064	…	4.36E-09	4.04E-10

当 $n = 8$ 时,求得得 $E$ 小于预定精度,因此得到
$\begin{aligned} & \beta = \beta_8 = 0.9309 \\& C = C_8 = 0.5406 \end{aligned}$

根据上述结果求出 $d_t$ ,得
$d_t = D_t \beta = 290.15 \times 0.9309 = 270.100635mm$

验算压力损失
$\begin{aligned} \delta p_{ch} = \frac{\sqrt{1-\beta^4}-C\beta^2}{\sqrt{1-\beta^4}+C\beta^2} \Delta P_{ch} = \frac{\sqrt{1-0.9309^4}- 0.5406 \times 0.9309^2}{\sqrt{1-0.9309^4}+0.5406 \times 0.9309^2} \times 90 = 28.85 kPa < \delta P_y \end{aligned}$
压力损失验算合格.
验算流量
$\begin{aligned} q_m & = 0.004 \frac{C}{\sqrt{1-\beta^4}}\varepsilon d_t^2 \sqrt{\Delta P \rho_1} \\ & = 0.004 \times \frac{0.5406}{\sqrt{1-0.9309^4}}\times 0.9965 \times 270.10^2 \times \sqrt{90000 \times 40.326} \\ & = 600.11503t/h \approx q_{m \, ch} \end{aligned}$
流量验算合格
确定 $d_{20}$ .标准孔板材料为29号钢,查表得其线膨胀系数 $\lambda_d=14.38 \times 10^{-6} \, /^{\circ}C$ ,则
$d_{20} = \frac{d_t}{1+ \lambda_d (t-20)} = \frac{210.100635}{1+ 14.38 \times 10^{-6} \times (540-20)} = 208.54127mm$
确定安装位置.根据 $\beta=0.9309$ 和管路系统,查表得到
$\begin{aligned} & l_0 = 18D = 5.184m \\ & l_1 = 50D =14.4m \\ & l_2 = 8D = 2.304m\\ \end{aligned}$
计算流量测量的总不确定度.
- 流出系数 $C$ 的不确定度在上边程序中计算得到: $\frac{\delta C}{C} = \pm 0.540910\%$
- 可膨胀性系数的不确定度 $\frac{\delta \varepsilon}{\varepsilon} = \pm 4 \frac{\Delta \varepsilon}{\varepsilon} = \pm 0.027271 \%$
- 差压变送器不确定度 $\frac{\delta P}{P} = \frac{2}{3} \times \frac{1}{90000}\times (0.25\% \times 160000) = 0.2962\%$
- 另外,设温度,压力,密度测量的不确定度为 $\pm1\%$
则,总不确定度为
$\begin{aligned} \frac{\delta_{q_m}}{q_m} & = \pm \left[ (\frac{\delta_C}{C})^2 + (\frac{\delta_\varepsilon}{\varepsilon})^2 + (\frac{2\beta^4}{1-\beta^4})^2(\frac{\delta_D}{D})^2 + (\frac{2}{1-\beta^4})^2(\frac{\delta_d}{d})^2 + \frac{1}{4}(\frac{\delta_{\Delta P}}{{\Delta P}})^2 + (\frac{\delta_P}{P})^2 \right] \\ & = \pm \left\{0.540910^2 + 0.027271^2 + (\frac{0.9309^4}{1-0.9309^4})^2\times 0.1^2 + (\frac{2}{1-0.9309^4})^2\times0.05^2 +0.25\times(0.2962)^2 + 0.25 \right\}^{\frac{1}{2}} \% \\ & = \pm 0.94 \% \end{aligned}$

思考题: 讨论设计题目给出的条件变化对设计结果有何影响

直接探讨各个条件变化如何影响设计结果较困难,在这里,我们通过考察迭代过程的初始变量对迭代结果有何影响,再讨论题目给出的条件如何影响迭代的初始变量,从而得到设计题目给出的条件对设计结果的影响.

讨论初始参数对设计结果的影响

使用控制变量法,在保持其它变量处于本例数值上时,改变某一初始参数,从而探讨初始参数对设计结果的影响

不变量 $A$ 变化对设计结果的影响

令不变量 $A$ 在0~:1内变化,得到设计结果如下.

我们发现直径比 $\beta$ 与不变量 $A$ 呈正相关关系,膨胀系数 $\varepsilon$ 与不变量 $A$ 呈负相关关系,流出系数 $C$ 在 $A < = 0.352$ 之前,受 $A$ 影响不大,在 $A > 0.352$ 之后,流出系数 $C$ 与不变量 $A$ 呈负相关关系.
雷诺数 $R_{eD}$ 变化对设计结果的影响

令雷诺数 $R_{eD}$ 在1~1e7之间变化,得到设计结果如下.

我们发现在 $R_{eD} \in [1,1e4]$ 时,直径比 $\beta$ 与雷诺数 $R_{eD}$ 成正相关,流出系数 $C$ ,膨胀系数 $\varepsilon$ 与雷诺数 $R_{eD}$ 成负相关关系.在 $R_{eD} > 1e4$ 时,雷诺数 $R_{eD}$ 变化对设计结果几乎无影响. 可见在湍流状态下,雷诺数对设计结果的影响可以忽略.
管道内径 $D$ 变化对设计结果的影响

令雷诺数 $R_{eD}$ 在1mm~1000mm之间变化,得到设计结果如下.

我们发现管道内径 $D$ 变化对设计结果的影响趋势与管道内径大致相同.在 $\in [1,200]$ 时,直径比 $\beta$ 与管道内径 $D$ 成正相关,流出系数 $C$ ,膨胀系数 $\varepsilon$ 与管道内径 $D$ 成负相关关系.在 $D > 200$ 时,管道内径 $D$ 变化对设计结果影响不大.
工作压力 $P_1$ 变化对设计结果的影响:

令工作压力 $P_1$ 在1~1e8Pa之间变化,得到设计结果如下.

在工作压力 $P_1 \in [1,1e3]Pa$ 时,直径比 $\beta$ 与管道内径 $D$ 成正相关,流出系数 $C$ ,膨胀系数 $\varepsilon$ 与管道内径 $D$ 成负相关关系.在 $D > 1 e 3 P a$ 时,管道内径 $D$ 变化对设计结果影响不大.
差压变送器差压 $\Delta P$ 变化对设计结果的影响

令差压变送器差压 $\Delta P$ 在1~1e6Pa之间变化,得到设计结果如下.

我们发现,在差压变送器压差 $\Delta P < 0.1 MPa$ 时,差压变送器的最大差压随设计结果几乎不产生影响.

综上所铸,对当前模型来说,初始参数中不变量 $A$ 和管道内径 $D$ 能对设计结果造成较大影响,雷诺数 $R_{eD}$ ,工作压力 $P_1$ 和差压变送器选型对设计结果的影响不大.

讨论思考题中各个量的变化对设计结果的影响

在这里,我们通过分析思考题中给出的各个物理量的变化是如何影响迭代的初始参数来考察各个物理量是如何影响设计结果的.

工作压力

工作压力 $P_1$ 变化对设计结果的影响在上一节中已经讨论过,该节流装置工作压力 $P_1 = 13.2MPa$ ,根据我们画出的图像来看,在该数量级上工作压力 $P_1$ 的变化对设计结果几乎不产生影响. 这也很符合实际情况.
工作温度

工作温度 $t_1$ 的变化会影响管道内径 $D$ ,而管道内径 $D$ 的变化会引起 $A$ 的变化.曲线斜率可以看到,在当前情况下 $A$ 对设计结果的影响比 $D$ 的影响大. 因此当 $t_1$ 增大时 $A$ 减小引起的作用为主要作用,此时 $A$ 减小会导致直径比 $\beta$ 减小,膨胀系数 $\varepsilon$ 和流出系数 $C$ 增大
常用流量

常用流量 $q_{m \, ch}$ 的变化会同时影响 $A$ 和雷诺数 $R_{eD}$ .当前情况下雷诺数 $R_{eD}=2.37 \times 10^7$ 常用流量的改变难以改变当前处于湍流状态的事实,因此雷诺数的变化对设计结果的影响可以忽略.但常用流量的增大会引起 $A$ 的增大,从而导致直径比 $\beta$ 增大,膨胀系数 $\varepsilon$ 和流出系数 $C$ 减小.
最大流量

最大流量的变化会导致常用流量对应的差压 $\Delta P_{ch}$ 的变化,进而影响 $A$ .当最大流量增加, $\Delta P_{ch}$ 减小导致 $A$ 增大,从而导致直径比 $\beta$ 增大,膨胀系数 $\varepsilon$ 和流出系数 $C$ 减小.
最小流量

我们的设计过程中并没有涉及到最小流量,因此我认为最小流量的变化对设计结果几乎没有影响.但是要注意的是,当最小流量过小,会导致雷诺数 $R_{eD}$ 过小,使得管内流体的流动状态从湍流变为层流,测量效果受到影响.

pdf版本的下载地址: 过程参数检测及仪表课程设计-第二类标准节流装置设计命题（访问密码：3834）

ncepu_Chen

关注

11
点赞
踩
39

收藏

觉得还不错? 一键收藏
5
评论
过程参数检测及仪表课程设计-第二类标准节流装置设计命题

过程参数检测及仪表课程设计-第二类标准节流装置设计命题设计任务设计过程设计任务序号项目符号单位参数1被测介质名称过热蒸汽2常用流量qm&ThinSpace;chq_{m \, ch}qmcht/h6003最小流量qm&ThinSpace;minq_{m\, min}qmmint/h4004最大流量qm&...
复制链接

扫一扫