CN1673693A - 一种简易的容器体积精确标定方法 - Google Patents
一种简易的容器体积精确标定方法 Download PDFInfo
- Publication number
- CN1673693A CN1673693A CN 200510046319 CN200510046319A CN1673693A CN 1673693 A CN1673693 A CN 1673693A CN 200510046319 CN200510046319 CN 200510046319 CN 200510046319 A CN200510046319 A CN 200510046319A CN 1673693 A CN1673693 A CN 1673693A
- Authority
- CN
- China
- Prior art keywords
- container
- pressure
- helium
- temperature
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 239000001307 helium Substances 0.000 claims abstract description 82
- 229910052734 helium Inorganic materials 0.000 claims abstract description 82
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 claims abstract description 82
- 229910000831 Steel Inorganic materials 0.000 claims abstract description 41
- 239000010959 steel Substances 0.000 claims abstract description 41
- 239000007789 gas Substances 0.000 claims abstract description 18
- 238000012360 testing method Methods 0.000 claims description 21
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 5
- 239000001257 hydrogen Substances 0.000 claims description 5
- 229910052739 hydrogen Inorganic materials 0.000 claims description 5
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 230000002045 lasting effect Effects 0.000 claims description 3
- 238000001179 sorption measurement Methods 0.000 abstract 1
- 238000012935 Averaging Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 239000011324 bead Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 102220008982 rs187686559 Human genes 0.000 description 1
- 102220223891 rs370721650 Human genes 0.000 description 1
Images
Landscapes
- Filling Or Discharging Of Gas Storage Vessels (AREA)
- Examining Or Testing Airtightness (AREA)
Abstract
本发明涉及一种简易的精确标定容器体积的方法。采用性质稳定、器壁吸附小的氦气作为标定气体,通过连续记录操作过程中各部分容器内的氦气压力及温度,根据氦气的MBWR方程和质量守恒定律建立以第一容器和第二容器的体积为未知数的二元一次方程。然后通过添加已知体积的钢球,改变其中一个容器体积,进行相同的操作。同样,根据氦气的MBWR方程和质量守恒定律建立以第一容器和添加钢球的第二容器的体积为未知数的二元一次方程。由于记录的数据有多组,所以可以建立多个二元一次方程组,解这些方程组可以获得多个容器体积值,计算出这些数值的平均值和中值,由操作者自行决定选用何值作为容器最终标定值。采用本发明可以精确标定容器体积。
Description
技术领域
本发明涉及标定容器体积的一种方法,特别提供了一种简易的精确标定容器体积的方法。
背景技术
当采用容量法测试材料的吸附气体曲线时必须首先标定测试设备各部分的体积。目前,基本通过量筒排水法标定各部分体积。其基本原理是将需要标定体积的容器充满气体,然后向装有水的量筒中排气,由于气体的充入使一部分体积的水被排除,通过容器中两次的水位差确定充入的气体的体积,同时移动外置胶管的位置使胶管内部液面和量筒内液面相平从而保证量筒内部所封闭的气体的压力为1个大气压,通过理想气体状态方程推算出需要标定部分的体积。该方法操作繁琐,通常需要两个人配合进行。另外由于量筒的读数精度有限,以及大气压力随地点,时间有变化,测出的容器体积有较高误差。
发明内容
本发明的目的在于提供一种简易的容器体积精确标定方法。
本发明的技术方案是:
一种简易的容器体积精确标定方法,包括如下步骤:
1)以氦气作为标定气体;
2)通过温度变送器测定所在的测试环境的温度,压力变送器测定通过阀门相连的第一容器和第二容器的实验压力;
3)首先通过真空泵对整个测试系统抽真空持续10~20分钟;停止抽真空,然后向系统内充入氦气使系统保持2~20kPa的压力,待压力变送器的读数稳定后,关闭阀门,记录此时的温度和压力分别记为Tc11、Py11;
4)向第一容器中充入压力为0.1~2MPa的氦气,测定并记录此时容器的温度Tc11、压力Pc11,打开阀门,第一容器中的氦气进入第二容器,当压力变送器读数稳定后,记录此时系统平衡压力Pb11以及此时第一容器的温度Tb11;
5)重复上述步骤4)过程多次,得到一系列向系统内充气过程的温度和压力数组,进行整个系统向外排气过程的测试;其操作原理与步骤4)相同,只是改不断向第一容器中充气为第二容器不断向大气中排气;记录测试过程的温度和压力值Pc1i,Py1i,Pb1i,Tc1i,Tb1i;
6)选用精度等级高的相同尺寸的微小钢珠,钢球的精度等级为≥G20,直径为3~15mm,放入第二容器中;
7)重复上述步骤4)、5)过程10次以上,获得一系列温度和压力值Pc2i,Py2i,Pb2i,Tc2i,Tb2i;
8)以MBWR方程计算相应温度及压力状态下的氦气压缩因子Z(P,T),然后根据公式
计算出各部分氦气的摩尔量,其中,R为气体常数,T为绝对温度,P为氢气的压力;
9)根据上述步骤5)、8)可以获得关于第一容器和第二容器体积的多个二元一次方程组,通过解这些方程组获得多个第一容器和第二容器体积值,计算出相应的平均值和中值,并做出体积值、平均值和中值的分布图;
10)整个计算及绘图过程由以下程序完成,整个程序包含5个函数:mbwr-comfactorHe函数、mbwr-eos函数、f函数、vcvy函数、pz函数;其中mbwr-comfactorHe函数用于计算氦气的压缩因子,输入变量为氦气的压力和温度,输出变量为氦气的压缩因子;mbwr-eos函数用于计算氦气的密度,输入变量为氦气的压力和温度,输出变量为氦气的密度,所用数值计算方法为割线法;f函数是由上述氦气的mbwr方程转换而来,其输入变量为氦气的密度、压力和温度,输出变量为压力,供mbwr-eos函数调用进行迭代运算;vcvy函数用于计算各次测得的容器体积值和对应的均值和中值,以及作出各次测得体积值、中值、平均值分布图,其输入变量为各次操作中容器的实测压力和温度,放入的钢球个数和钢球单个体积,输出变量为各次测得的容器体积值和对应的均值和中值;pz函数用于计算修正压力值,其输入变量为氦气的压力和温度,输出变量为氦气的修正压力。
所述步骤3)中抽真空,使真空度达到10-1-10-2Pa。
所述步骤5)重复步骤4)的次数为20-30次。
本发明的有益效果是:
本发明采用性质稳定、器壁吸附小的氦气作为标定气体,通过连续记录操作过程中各部分容器内的氦气压力及温度,根据氦气的MBWR方程和质量守恒定律建立以第一容器和第二容器的体积为未知数的二元一次方程。然后通过添加已知体积的钢球,改变其中一个容器体积,进行相同的操作。同样,根据氦气的MBWR方程和质量守恒定律建立以第一容器和添加钢球的第二容器的体积为未知数的二元一次方程。由于记录的数据有多组,所以可以建立多个二元一次方程组,解这些方程组可以获得多个容器体积值,计算出这些数值的平均值和中值,由操作者自行决定选用何值作为容器最终标定值。采用本发明可以精确标定容器体积。
附图说明
图1为本发明方法使用装置的简易示意图。图中,1第一容器;2温度变送器;3压力变送器;4阀门;5钢球;6第二容器。
图2为未添加钢珠及阀门4打开前,各部分内氦气的摩尔量。
图3为未添加钢珠及阀门4打开后,各部分内氦气的摩尔量。
图4为已添加钢珠及阀门4打开前,各部分内氦气的摩尔量。
图5为已添加钢珠及阀门4打开后,各部分内氦气的摩尔量。
图6为已测得第一容器1的多个体积值及其中值和平均值分布。
图7为已测得第二容器6的多个体积值及其中值和平均值分布。
图8a-c为MATLAB程序流程图。
具体实施方式
本发明方法使用装置如图1所示,该装置包括第一容器1、温度变送器2、压力变送器3、阀门4、钢球5、第二容器6,第一容器1上装有温度变送器2、压力变送器3,第二容器6内装有钢球5,第一容器1与第二容器6通过管道相连,管道上设有阀门4。
本发明提供了一种简易的容器体积精确标定方法,包括如下步骤:
1)以性质稳定,器壁吸附小的氦气作为标定气体。
2)通过温度变送器3测定所在的测试环境的温度,压力变送器2测定第一容器1和第二容器6的实验压力。
3)首先通过直联真空泵对整个测试系统抽真空持续10分钟,使真空度达到10-1-10-2Pa。停止抽真空,然后向系统内充入少量的氦气,使系统保持10kPa的压力,待压力变送器的读数稳定后,关闭阀门4,记录此时的温度和压力分别记为Tc11,Py11。
4)向第一容器1中充入一定压力的氦气(压力为0.1-2MPa),测定并记录此时第一容器1的温度Tc11、压力Pc11,打开阀门4,第一容器1中的氦气进入第二容器6,当压力变送器读数稳定后,记录此时系统平衡压力Pb11(该压力值既是第二容器6的压力也是第一容器1的压力)以及此时第一容器1的温度Tb11。
5)重复上述步骤4)过程多次,得到一系列向系统内充气过程的温度和压力数组。为使测出的体积更加反应真实情况,需进行整个系统向外排气过程的测试。其操作原理与步骤4)相同,只是改不断向第一容器1中充气为第一容器1不断向大气中排气。记录测试过程的温度和压力值Pc1i,Py1i,Pb1i,Tc1i,Tb1i。
6)选用表面光洁度及尺寸精度高的相同尺寸的微小钢珠若干,钢球的精度等级为≥G20,直径为3~15mm,放入第二容器6中。
7)重复上述步骤4)、5)过程多次,获得一系列温度和压力值Pc2i,Py2i,Pb2i,Tc2i,Tb2i。
8)以目前公认的描述氦气实际状态最佳的方程MBWR方程,计算相应于一定温度及压力状态下的氦气压缩因子Z(P,T),然后根据公式
计算出各部分氦气的摩尔量,由于氦气不会被吸附所以认为气态氦气的摩尔量在每次的测试过程中是保持不变的。
9)根据上述步骤5)、8)可以获得关于第一容器1和第二容器6体积的多个二元一次方程组,通过解这些方程组获得多个第一容器1和第二容器6体积值,计算出相应的平均值和中值,并做出体积值、平均值和中值的分布图。
10)整个计算及作图过程由自行编制MATLAB程序完成。
本发明采用性质稳定,器壁吸附小的氦气作为标定气体,通过连续记录操作过程中各部分容器内的氦气压力及温度,根据氦气的MBWR方程和质量守恒定律建立以第一容器1和第二容器6的体积为未知数的二元一次方程。然后通过添加已知体积的钢球5,改变其中一个容器体积(以下以改变第二容器6的体积进行操作),进行相同的操作。同样,根据氦气的MBWR方程和质量守恒定律建立以第一容器1和添加钢球5的第二容器6的体积为未知数的二元一次方程。由于记录的数据有多组,所以可以建立多个二元一次方程组,解这些方程组可以获得多个第一容器1和第二容器6体积值,将这些数值进行平均,取其平均值作为标定的体积。技术关键是如何通过计算机编程实现压缩因子的计算。氦气的MBWR方程一个非常复杂的方程,其形式为:
P=ρRT+ρ2(N0T+N1T1/2+N2+N3/T+N4/T2)
+ρ3(N5T+N6+N7/T+N8/T2)
+ρ4(N9T+N10+N11/T)+ρ5(N12)
+ρ6(N13/T+N14/T2)+ρ7(N15/T)
+ρ8(N16/T+N17/T2)+ρ9(N18/T2)
+ρ3(N19/T2+N20/T3)exp(γρ2)
+ρ5(N21/T2+N22/T4)exp(γρ2)
+ρ7(N23/T2+N24/T3)exp(γρ2)
+ρ9(N25/T2+N26/T4)exp(γρ2)
+ρ11(N27/T2+N28/T3)exp(γρ2)
+ρ13(N29/T2+N30/T3+N31/T4)exp(γρ2)
其中N0~N31,γ为描述氦气性质的参数,ρ为氦氢气密度,R为气体常数,T为绝对温度,P为氦气的压力,该方程描述的是氦气压力与其密度及温度的关系。通过压力变送器和温度变送器可以测出气体的压力及温度,通过编程实现计算机自动计算出氢气密度ρ,所采用的数值计算方法为割线法。然后,根据公式 计算出对应于该压力,温度下的氢气压缩因子Z(P,T)。由图2可得未添加钢珠及阀门4打开前第一容器1内氦气摩尔量为n1i,第二容器6内氦气摩尔量为n2i,由图3可得未添加钢珠及阀门4打开后第一容器1内氦气摩尔量为n3i,第二容器6内氦气摩尔量为n4i,,由于氦气不会被器壁吸附所以n1i+n2i=n3i+n4i。同样当添加已知体积的小球后建立可以建立相似的等式n1i’+n2i’=n3i’+n4i’。因为在实测过程中,由于每次操作过程时间很短,室温温度波动很小,所以可以得出tc1i=tb1i=t1o,tc2i=tb2i=t2o。通过上面的等式建立方程组,解出第一容器1和第二容器6的体积,由于建立了多个二元一次方程组所以解这些方程组会产生可以获得多个第一容器1和第二容器6的体积值,计算出相应的平均值和中值,并做出体积值、平均值和中值的分布图,供操作者判定实验测得的体积的可靠性,一般以中值(median)作为最终的标定值。
如图2所示,未添加钢珠及阀门4打开前,各部分内氦气的摩尔量;如图3所示,未添加钢珠及阀门4打开后,各部分内氦气的摩尔量;如图4所示,已添加钢珠及阀门4打开前,各部分内氦气的摩尔量;如图5所示,已添加钢珠及阀门4打开后,各部分内氦气的摩尔量;如图6所示,已测得第一容器1的多个体积值及其中值和平均值分布;如图7所示,已测得第二容器6的多个体积值及其中值和平均值分布。
如图8a-c所示MATLAB程序流程图,整个程序包含5个函数:mbwr-comfactorhe函数、mbwr-eos函数、f函数、vcvy函数、pz函数。其中,如图8b所示,mbwr-comfactorhe函数用于计算氦气的压缩因子,输入变量为氦气的压力和温度,输出变量为氦气的压缩因子;如图8c所示,mbwr-eos函数用于计算氦气的密度,输入变量为氦气的压力和温度,输出变量为氦气的密度,所用数值计算方法为割线法;如图8c所示,f函数是由上述氦气的mbwr方程转换而来,其输入变量为氦气的密度、压力和温度,输出变量为压力,此函数没有物理意义,目的是供mbwr-eos函数调用进行迭代运算;如图8a所示,vcvy函数用于计算各次测得的第一容器1和第二容器6体积值和对应的均值和中值,以及作出各次测得的体积值、中值、平均值分布图,其输入变量为各次操作中第一容器1和第二容器6的实测压力和温度,放入的钢球个数和钢球单个体积,输出变量为各次测得的第一容器1和第二容器6体积值和对应的均值和中值;如图8a所示,pz函数用于计算修正压力值,其输入变量为氦气的压力和温度,输出变量为氦气的修正压力。
以下是整个计算过程的MATLAB程序。
该部分程序用于计算氦气的压缩因子:
function z=mbwr_comfactorHe(P,T)
%z=mbwr_comfactorHe(P,T);
%The 32-term equation of state;
%Based on the Modified Bennedict-webb-rubin equation of state;
%P-kPa;
%T--k;
%在极端环境下有可能计算出的压缩因子不准;
%不过在常规温度273.15~573.15k,及常规压力0.001~100Mpa下,该程序计算出的氦气
%压缩因子是可靠的,而标定各部分体积时的测试环境包含于这个范围内;
%使用的数值计算方法为割线法。
R=8.314510;
ro=mbwr_eos(P,T);
z=P./ro./R./T;
function ro=mbwr_eos(P,T)
R=8.314510;
ro0=0.002;
ro1=P./R./T;
y=1;
while abs(y)>1e-7
ro2=(ro0.*feval(@f,ro1,P,T)-ro1.*feval(@f,ro0,P,T))./...
(feval(@f,ro1,P,T)-feval(@f,ro0,P,T));
ro0=ro1;
ro1=ro2;
y=feval(@f,ro1,P,T);
end
ro=ro1;
function y=f(ro,P,T)
rr=0.00831434;
gm=-0.33033259E-02;
g0=.4558980227431E-04;g1=.1260692007853E-02;g2=-.7139657549318E-02;
g3=.9728903861441E-02;g4=-.1589302471562E-01;g5=.1454229259623E-05;g6=-.4708238429298E-04;g7=.1132915223587E-02;g8=.2410763742104E-02;g9=-.5093547838381E-08;g10=.2699726927900E-05;g11=-.3954146691114E-04;
g12=-.1551961438127E-08;g13=.1050712335785E-07;g14=-.5501158366750E-07;
g15=-.1037673478521E-09;g16=.6446881346448E-12;g17=.3298960057071E-10;
g18=-.3555585738784E-12;g19=-.6885401367690E-02;g20=.9166109232806E-02;
g21=-.6544314242937E-05;g22=-.3315398880031E-04;g23=-.2067693644676E-07;
g24=.3850153114958E-07;g25=-.1399040626999E-10;g26=-.1888462892389E-11;
g27=-.4595138561035E-14;g28=.6872567403738E-14;g29=-.6097223119177E-18;
g30=-.7636186157005E-17;g31=.3848665703556E-17;
d0=ro;d1=ro.^2;d2=ro.^3;d3=ro.^4;d4=ro.^5;d5=ro.^6;d6=ro.^7;d7=ro.^8;d8=ro.^9;
d9=ro.^10;d10=ro.^11;d11=ro.^12;d12=ro.^13;
t0=sqrt(T);t1=T.^2;t2=T.^3;t3=T.^4;
f=exp(gm.*d1);
b0=d1.*T; b1=d1.*t0; b2=d1;
b3=d1./T; b4=d1./t1; b5=d2.*T;
b6=d2; b7=d2./T; b8=d2./t1;
b9=d3.*T; b10=d3; b11=d3./T;
b12=d4; b13=d5./T; b14=d5./t1;
b15=d6./T; b16=d7./T; b17=d7./t1;
b18=d8./t1;b19=d2.*f./t1;b20=d2.*f.t2;
b21=d4.*f./t1;b22=d4.*f./t3;b23=d6.*f./t1;
b24=d6.*f./t2;b25=d8.*f./t1;b26=d8.*f./t3;
b27=d10.*f./t1;b28=d10.*f./t2;b29=d12.*f./t1;
b30=d12.*f./t2;b31=d12.*f./t3;
p=(b0.*g0+b1.*g1+b2.*g2+b3.*g3+b4.*g4+b5.*g5+b6.*g6+b7.*g7+b8.*g8+b9.*g9+...b10.*g10+b11.*g11+b12.*g12+b13.*g13+b14.*g14+b15.*g15+b16.*g16+b17.*g17+...b18.*g18+b19.*g19+b20.*g20+b21.*g21+b22.*g22+b23.*g23+b24.*g24+b25.*g25+...b26.*g26+b27.*g27+b28.*g28+b29.*g29+b30.*g30+b31.*g31+rr.*d0.*T).*1000;
y=p-P;
该部分程序用于计算第一容器1和第二容器6的体积标定值和做出体积,中值,平均值分布图,用于判定测试值可靠性:
function[mean Vc,median Vc,mean Vy,median Vy,Vc,Vy]=vcvy(Pc1,Pc2,Py1,Py2,Pb1,Pb2,t10,t20,n,V0)
%[mean Vc,median Vc,mean Vy,median Vy,Vc,Vy]=vcvy(Pc1,Pc2,Py1,Py2,Pb1,Pb2,t10,t20,n,V0)
%Pc1--未添加钢珠及阀门4打开之前测得的第一容器1的压力,单位为kPa;
%Pc2--已添加钢珠及阀门4打开之前测得的第一容器1的压力,单位为kPa;
%Py1--未添加钢珠及阀门4打开之前测得的第二容器6的压力,单位为kPa;
%Py2--已添加钢珠及阀门4打开之前测得的第二容器6的压力,单位为kPa;
%Pb1--未添加钢珠及阀门4打开之后测得的系统的压力,单位为kPa;
%Pb2--已添加钢珠及阀门4打开之后测得的系统的压力,单位为kPa;
%t10--未添加钢珠的测试过程的温度,单位为℃;
%t20--已添加钢珠的测试过程的温度,单位为℃;
%n--添加钢珠个数;
%V0--单个钢珠的体积单位为ml;
%mean Vc--所计算出第一容器1的体积平均值;
%median Vc--所计算出第一容器1的体积中值;
%mean Vy--所计算出第二容器6的体积平均值;
%median Vy--所计算出第二容器6的体积中值;
%Vc--计算出第一容器1的体积数组;
%Vy--计算出第二容器6的体积数组;
T10=t10+273.15;
T20=t20+273.15;
n11=pz(Pc1,T10);
n12=pz(Pb1,T10);
n13=pz(Py1,T10);
c11=(n11-n12)./(n12-n13);
n21=pz(Pc2,T20);
n22=pz(Pb2,T20);
n23=pz(Py2,T20);
c21=(n21-n22)./(n22-n23);
Vy=n.*V0.*c21./(c21-c11);
mean Vy=mean(Vy);
median Vy=median(Vy);
a=mean Vy.*ones(length(Vy),1);
b=median Vy.*ones(length(Vy),1);
Vc=c1.*Vy;
mean Vc=mean(Vc);
median Vc=median(Vc);
c=mean Vc.*ones(length(Vy),1);
d=median Vc.*ones(length(Vy),1);
%分别做出体积,中值,平均值分布图,用于判定测试值可靠性。
plot(Vy,′rs′,′markerfacecolor′,′g′,′markersize′,10);
hold on
plot(a,′b-′);
hold on
plot(b,′r-′);
legend(′Vy′,′mean Vy′,′median Vy′);
figure
plot(Vc,′rs′,′markerfacecolor′,′g′,′markersize′,10);
hold on
plot(c,′b-′);
hold on
plot(d,′r-′);
legend(′Vc′,′mean Vc′,′median Vc′);
function y=pz(P,T)
y=P./mbwr_comfactorHe(P,T);
首先对本发明方法使用装置进行抽真空,持续10分钟。然后向系统充入氦气,10分钟后,待读数基本稳定后记录压力值和温度值,然后关闭阀门4,不断向其中的第一容器1充气,记录压力及温度值,然后打开阀门4,待读数稳定后记录压力值,获得25个相关的温度和压力数据。
然后添加7个直径为8mm的钢珠到第二容器6中,进行上述同样的试验,最终获得25个相关的温度和压力数据。两次测得的温度压力数据如表1所示。利用上述matlab程序计算得到第一容器1的体积中值为26.3655ml,体积平均值为24.4868ml,第二容器6的体积中值为22.5259ml,体积平均值为20.9895ml。对照图6和图7可以知道选择体积中值作为标定值更可靠。
表1
Pb1/kPa | Pb2/kPa | Pc1/kPa | Pc2/kPa | Py1/kPa | Py2/kPa | t10/℃ | t20/℃ |
314428642739229421031314128811646074392592393275922691465168417221151692508291760857638384 | 286533523696362932142667234220942028143512918821278202226121861995496252126502344189610519412351 | 2817262426321933173364312641057132295106426230751853884187217556653003501061163940451165 | 1863373639673577289022462089190119769731178563158926093079127632110641442752210615473948553837 | 3536314428642739229421031314128811646074392592393275921491465168417221151692508291760857638 | 416428653352369636293197266723422094202814351291882127820222612186199549625212651234418961049941 | 16171716.516.516.516.516.516.516.516.516.516.516.816.916.916.916.91716.516.516.816.816.516.5 | 16.516.516.516.516.51616161616161616161616161616161616161616 |
Claims (3)
1、一种简易的容器体积精确标定方法,其特征在于包括如下步骤:
1)以氦气作为标定气体;
2)通过温度变送器测定所在的测试环境的温度,压力变送器测定通过阀门相连的第一容器和第二容器的实验压力;
3)首先通过真空泵对整个测试系统抽真空持续10~20分钟;停止抽真空,然后向系统内充入氦气使系统保持2~20kPa的压力,待压力变送器的读数稳定后,关闭阀门,记录此时的温度和压力分别记为Tc11、Py11;
4)向第一容器中充入压力为0.1~2MPa的氦气,测定并记录此时容器的温度Tc11、压力Pc11,打开阀门,第一容器中的氦气进入第二容器,当压力变送器读数稳定后,记录此时系统平衡压力Pb11以及此时第一容器的温度Tb11;
5)重复上述步骤4)过程多次,得到一系列向系统内充气过程的温度和压力数组,进行整个系统向外排气过程的测试;其操作原理与步骤4)相同,只是改不断向第一容器中充气为第二容器不断向大气中排气;记录测试过程的温度和压力值Pc1i,Py1i,Pb1i,Tc1i,Tb1i;
6)选用精度等级高的相同尺寸的微小钢珠,钢球的精度等级为≥G20,直径为3~15mm,放入第二容器中;
7)重复上述步骤4)、5)过程10次以上,获得一系列温度和压力值Pc2i,Py2i,Pb2i,Tc2i,Tb2i;
8)以MBWR方程计算相应温度及压力状态下的氦气压缩因子Z(P,T),然后根据公式
计算出各部分氦气的摩尔量,其中,R为气体常数,T为绝对温度,P为氢气的压力;
9)根据上述步骤5)、8)可以获得关于第一容器和第二容器体积的多个二元一次方程组,通过解这些方程组获得多个第一容器和第二容器体积值,计算出相应的平均值和中值,并做出体积值、平均值和中值的分布图;
10)整个计算及绘图过程由以下程序完成,整个程序包含5个函数:mbwr-comfactorHe函数、mbwr-eos函数、f函数、vcvy函数、pz函数;其中mbwr-comfactorHe函数用于计算氦气的压缩因子,输入变量为氦气的压力和温度,输出变量为氦气的压缩因子;mbwr-eos函数用于计算氦气的密度,输入变量为氦气的压力和温度,输出变量为氦气的密度,所用数值计算方法为割线法;f函数是由上述氦气的mbwr方程转换而来,其输入变量为氦气的密度、压力和温度,输出变量为压力,供mbwr-eos函数调用进行迭代运算;vcvy函数用于计算各次测得的容器体积值和对应的均值和中值,以及作出各次测得体积值、中值、平均值分布图,其输入变量为各次操作中容器的实测压力和温度,放入的钢球个数和钢球单个体积,输出变量为各次测得的容器体积值和对应的均值和中值;pz函数用于计算修正压力值,其输入变量为氦气的压力和温度,输出变量为氦气的修正压力。
2、按照权利要求1所述简易的容器体积精确标定方法,其特征在于:所述步骤3)中抽真空,使真空度达到10-1-10-2Pa。
3、按照权利要求1所述简易的容器体积精确标定方法,其特征在于:所述步骤5)重复步骤4)的次数为20-30次。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100463199A CN1332185C (zh) | 2005-04-27 | 2005-04-27 | 一种简易的容器体积精确标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100463199A CN1332185C (zh) | 2005-04-27 | 2005-04-27 | 一种简易的容器体积精确标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1673693A true CN1673693A (zh) | 2005-09-28 |
CN1332185C CN1332185C (zh) | 2007-08-15 |
Family
ID=35046375
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2005100463199A Expired - Fee Related CN1332185C (zh) | 2005-04-27 | 2005-04-27 | 一种简易的容器体积精确标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1332185C (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101470016B (zh) * | 2007-12-28 | 2010-09-08 | 中国航天科技集团公司第五研究院第五一〇研究所 | 一种内置法测量容积的装置 |
CN103375164A (zh) * | 2012-04-13 | 2013-10-30 | 中国石油天然气股份有限公司 | 管线冻堵位置的判断方法及其使用的气体分压撬装设备 |
CN105115559A (zh) * | 2015-09-06 | 2015-12-02 | 中国人民解放军63653部队 | 基于二次微变容差压原理的容器容积测量装置及测量方法 |
CN105784054A (zh) * | 2016-03-15 | 2016-07-20 | 中国原子能科学研究院 | 一种适于标定不规则形状气腔容积的系统及方法 |
CN105980657A (zh) * | 2014-03-12 | 2016-09-28 | 界标制图有限公司 | 多种流体的共享状态方程表征 |
CN107490416A (zh) * | 2017-05-15 | 2017-12-19 | 宝沃汽车(中国)有限公司 | 测量车辆行李箱体积的方法 |
CN108759733A (zh) * | 2018-08-22 | 2018-11-06 | 上海力信能源科技有限责任公司 | 以气体测量不规则物体体积的装置及方法 |
CN109405919A (zh) * | 2018-11-27 | 2019-03-01 | 成都师范学院 | 白酒发酵糟醅体积测量装置及其使用方法 |
CN111912484A (zh) * | 2020-09-08 | 2020-11-10 | 东风汽车集团有限公司 | 快速测量汽车后备箱容积的装置及方法 |
CN116146445A (zh) * | 2022-07-12 | 2023-05-23 | 上海奔曜科技有限公司 | 机器人工作站中注射泵的标定方法和机器人工作站 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5484763A (en) * | 1977-12-19 | 1979-07-05 | Mitsubishi Electric Corp | Volume measuring method and apparatus |
JPS5920814A (ja) * | 1982-07-28 | 1984-02-02 | Shimadzu Corp | 体積測定方法および装置 |
DE10041051B4 (de) * | 2000-08-22 | 2006-08-10 | Fti Technologies Gmbh | Verfahren zur Volumenmessung durch Druckstoßbestimmung |
JP3756919B2 (ja) * | 2004-09-07 | 2006-03-22 | 日本ベル株式会社 | 死容積の変動量測定方法 |
-
2005
- 2005-04-27 CN CNB2005100463199A patent/CN1332185C/zh not_active Expired - Fee Related
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101470016B (zh) * | 2007-12-28 | 2010-09-08 | 中国航天科技集团公司第五研究院第五一〇研究所 | 一种内置法测量容积的装置 |
CN103375164A (zh) * | 2012-04-13 | 2013-10-30 | 中国石油天然气股份有限公司 | 管线冻堵位置的判断方法及其使用的气体分压撬装设备 |
CN103375164B (zh) * | 2012-04-13 | 2016-07-13 | 中国石油天然气股份有限公司 | 管线冻堵位置的判断方法 |
CN105980657A (zh) * | 2014-03-12 | 2016-09-28 | 界标制图有限公司 | 多种流体的共享状态方程表征 |
CN105115559A (zh) * | 2015-09-06 | 2015-12-02 | 中国人民解放军63653部队 | 基于二次微变容差压原理的容器容积测量装置及测量方法 |
CN105784054B (zh) * | 2016-03-15 | 2018-11-30 | 中国原子能科学研究院 | 一种适于标定不规则形状气腔容积的系统及方法 |
CN105784054A (zh) * | 2016-03-15 | 2016-07-20 | 中国原子能科学研究院 | 一种适于标定不规则形状气腔容积的系统及方法 |
CN107490416A (zh) * | 2017-05-15 | 2017-12-19 | 宝沃汽车(中国)有限公司 | 测量车辆行李箱体积的方法 |
CN108759733A (zh) * | 2018-08-22 | 2018-11-06 | 上海力信能源科技有限责任公司 | 以气体测量不规则物体体积的装置及方法 |
CN109405919A (zh) * | 2018-11-27 | 2019-03-01 | 成都师范学院 | 白酒发酵糟醅体积测量装置及其使用方法 |
CN111912484A (zh) * | 2020-09-08 | 2020-11-10 | 东风汽车集团有限公司 | 快速测量汽车后备箱容积的装置及方法 |
CN111912484B (zh) * | 2020-09-08 | 2021-05-18 | 东风汽车集团有限公司 | 快速测量汽车后备箱容积的装置及方法 |
CN116146445A (zh) * | 2022-07-12 | 2023-05-23 | 上海奔曜科技有限公司 | 机器人工作站中注射泵的标定方法和机器人工作站 |
CN116146445B (zh) * | 2022-07-12 | 2024-03-19 | 上海奔曜科技有限公司 | 机器人工作站中注射泵的标定方法和机器人工作站 |
Also Published As
Publication number | Publication date |
---|---|
CN1332185C (zh) | 2007-08-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1673693A (zh) | 一种简易的容器体积精确标定方法 | |
CN1687773A (zh) | 一种贮氢材料压力浓度温度性能精确测试方法 | |
CN1299108C (zh) | 光谱测定同位素气体的仪器 | |
CN1501440A (zh) | 可在连接状态下进行质量流控制器检查的半导体制造装置 | |
CN1623084A (zh) | 检漏仪 | |
CN1319185A (zh) | 气体类型识别系统 | |
CN1044978A (zh) | 具有测量流量功能的蝶形阀和用它测量流量的方法 | |
CN101047915A (zh) | 第三代时分同步码分多址移动终端自动校准的方法和装置 | |
CN1839358A (zh) | 流量控制装置及其调整方法 | |
CN1371470A (zh) | 直管科里奥利流量计 | |
CN1208614C (zh) | 透气性测定装置 | |
CN1466675A (zh) | 确定二氧化碳吸收剂吸收能力的同位素气体分析仪及方法 | |
CN106197896B (zh) | 一种轮胎内胎或内衬层气密性测定装置及气密性测定方法 | |
CN1926406A (zh) | 流量测定方法以及流量测定装置 | |
CN1779422A (zh) | 利用充放气过程测试容积的方法 | |
CN113340765B (zh) | 一种分子筛材料吸附性能检测装置及方法 | |
CN103344521B (zh) | 一种岩石体密度无损无污染的高精度测定方法 | |
CN1588002A (zh) | 降低实时荧光pcr仪器定量分析系统误差的分析方法及应用 | |
CN115128027B (zh) | 一种紫外差分光谱气体传感器环境压力补偿结构及方法 | |
TW202308002A (zh) | 多個製程腔室壓力的控制方法及半導體製程設備 | |
CN209783706U (zh) | 一种针对气体流量计的测量校正装置 | |
CN1685079A (zh) | 薄膜成膜方法、薄膜成膜装置和薄膜成膜过程的监视方法 | |
CN213749512U (zh) | 一种烟草专用通风率标准棒的测量装置 | |
JP3690923B2 (ja) | フレッシュコンクリートエアメータおよび空気量測定方法 | |
CN112697632A (zh) | 一种煤岩、页岩大样量重量法等温吸附测量装置及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20070815 Termination date: 20110427 |