CN114354015B - 热电偶校准系统测量不确定度的确定方法 - Google Patents
热电偶校准系统测量不确定度的确定方法 Download PDFInfo
- Publication number
- CN114354015B CN114354015B CN202111275518.2A CN202111275518A CN114354015B CN 114354015 B CN114354015 B CN 114354015B CN 202111275518 A CN202111275518 A CN 202111275518A CN 114354015 B CN114354015 B CN 114354015B
- Authority
- CN
- China
- Prior art keywords
- temperature
- thermocouple
- point
- uncertainty
- equalizing block
- 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.)
- Active
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000003780 insertion Methods 0.000 claims abstract description 51
- 230000037431 insertion Effects 0.000 claims abstract description 51
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 6
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 6
- 230000035945 sensitivity Effects 0.000 claims description 17
- 238000009529 body temperature measurement Methods 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 15
- 238000012795 verification Methods 0.000 abstract description 5
- 238000005245 sintering Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000008014 freezing Effects 0.000 description 2
- 238000007710 freezing Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000009827 uniform distribution Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明属于计量技术领域,涉及一种热电偶校准系统测量不确定度的确定方法,包括:1)选择热电偶Sa和Sb;2)将上述热电偶分别插入均温块内,反向串联形成热电偶校准系统;3)调整热电偶5个不同插入深度、均温块8个旋转角度,以及10个不同的测温点条件下,热电偶对应的温度差,并取算术平均值;4)测量均温块截面上插入热电偶的离轴孔的中心位置距均温块截面几何中心点的距离x;5)带入合成不确定度公式得到结果。本发明能实现特定条件下,热电偶检定结果的不确定度的量化计算,从而确定热电偶校准结果的准确程度,扩展了热电偶校准过程中的条件限制,提升测量结果的应用价值。
Description
技术领域
本发明属于计量技术领域,涉及一种热电偶校准系统测量不确定度的确定方法,用以解决热电偶检定系统核心区域温度控制不确定度的定量评估。
背景技术
热电偶是温度测量仪表中常用的测温元件,为了保证温度测量结果的准确性,需要对热电偶进行校准。目前,现有的热电偶校准系统通常采用大长径比的管状炉膛结构,核心控温区位于管状炉膛的中部位置,校准时,将被测电偶及标准电偶插入核心控温区的恒温块中,加热核心控温区至指定温度,再通过读取标准热电偶及被测热电偶的测温读数差值,确定被测热电偶在该温度点的误差。
而对于热电偶校准系统本身误差对测量结果的影响,现有技术主要是通过引入热电偶校准系统本身的控温误差参数,将该部分误差引入被测热电偶测量误差的不确定度中,但是采用此技术存在以下问题:(1)仅能粗略估算由热电偶校准系统引入的不确定度,无法更为精确的评估热电偶校准系统的温度测量不确定度;(2)对于被测电偶在径向或轴向偏离核心控温区几何中心的时对测量结果不确定度的影响,缺乏有效的量化计算手段,例如当热电偶由于机械结构限制,难以插入核心区域几何中心位置时,采用现有的计算手段,难以对测量结果的不确定度做出量化的评定,影响最终测量结果的测量不确定度的确定。
发明内容
为了解决现有背景技术中存在的技术问题,本发明提供了一种热电偶校准系统测量不确定度的确定方法,实现特定条件下,热电偶检定结果的不确定度的量化计算,从而确定热电偶校准结果的准确程度,扩展了热电偶校准过程中的条件限制,提升测量结果的应用价值。
为了实现上述目的,本发明提供的技术方案是:
一种热电偶校准系统测量不确定度的确定方法,包括以下步骤:
1)根据JJF1637、JJF1262标准选择符合I级的S型热电偶为Sa、Sb;
2)将上述Sa、Sb分别插入均温块的中心孔内,并反向串联Sa、Sb,形成热电偶校准系统;
3)调整Sa、Sb插入中心孔中多个不同的深度、旋转均温块角度在多个不同位置处,以及多个不同的测温点条件下,Sa及Sb分别对应的温度差,并取算术平均值;得到不同插入深度j、不同旋转角度k以及不同测温点i时,热电偶校准系统对应的温度差,记为Ti,j,k(i=0,j=0,k=0)至Ti,j,k(i=10,j=4,k=7),共获得多组数据;
4)使用卡尺或设备测量均温块截面上插入热电偶的离轴孔的中心位置距均温块截面几何中心点的距离x;
5)根据步骤3)和步骤4)得到的数据,带入下列公式中,得到合成不确定度,
式中:μc为合成不确定度;
μ(e1)被测电偶输入量引入的不确定度分量;
c1被电偶输入量引入不确定度的灵敏系数;
μ(e2)标准热电偶本身引入的不确定度;
c2被测热电偶与标准热电偶在被检点上证书标注的的微分热电动势之比引入的灵敏系数;
μ(e3)检测时标准热电偶测量重复性引入的不确定度分量;
c3被测热电偶与标准热电偶在被检点上测得的微分热电动势算术平均值之比引入的灵敏系数;
μ(e4)使用补偿导线引入的不确定度分量;
c4补偿导线引入不确定度分量的灵敏系数;
所述
式中,r实为测量均温块上离轴分布插入孔圆心所构成的圆的半径;h标实为标准热电偶偏离中心位置的长度;h被实为被检偶的深度偏移以一次测量过程中偏离中心位置最远的被测偶计算;
所述μ(e2)通过标准热电偶的证书获得;
所述
所述式中α为补偿导线在测试时环境温度点的示值误差;
所述c4=1;
所述
式中:e被(t)为被检电偶在某被校准温度点的热电动势;
被检电偶在上述温度点多次测量的算术平均值;
e标证标准热电偶在上述温度点时,证书给出的热电动势;
标准电偶在上述温度点多次测量的算术平均值;
e补补偿导线修正值;
S被被检热电偶在上述温度点的微分电动势;
S标标准热电偶在上述温度点的微分电动势。
进一步的,所述步骤2)中,均温块长度不小于20mm,半径不小于28mm.
进一步的,所述步骤2)中,均温块的径向半径中心位置及离中心轴x(mm)位置处,具有至少一组以截面圆心呈角对称分布的热电偶埋入孔,x≥10mm深度10mm以上,埋入孔直径与步骤1)中的S型热电偶相匹配。
进一步的,所述步骤2)中,反向串联时,均温块外部的热电偶引出线短接点,并接入补偿导线。
进一步的,所述步骤3)的具体过程包括:
3.1)在100℃至1100℃之间,每100℃设置一个待检测温度点,在各待检测温度点读取Sa及Sb的温度差值;
3.2)若Sa的温度差值小于0.1℃每分钟,或者在任意温度点Sa及Sb的温度差值小于0.1℃,则进入步骤3.5);反之,则进行步骤3.3)~步骤3.4);
3.3)将Sa插入均温块的中心孔底部,Sb插入均温块任一离轴孔底部,通过旋转均温块,使Sb插入孔位于Sa插入孔的正上方0点位置;并将均温块插入深度调整到被测系统温区中心点向内深度10cm处;
3.4)按照步骤3.1)进行相关检测,直至Sa的温度差值小于0.1℃每分钟,或者在任意温度点Sa及Sb的温度差值小于0.1℃;
3.5)按照步骤3.2)读取的Sa及Sb温度差,并取算术平均值,记做Ti,j,k(i=0,j=2,k=0);
3.6)将均温块顺时针旋转45°,插入深度不变,按照步骤3.1)再次进行测量,直至Sa的温度差值小于0.1℃每分钟,或者在任意温度点Sa及Sb的温度差值小于0.1℃;
3.7)在步骤3.6)位置处,读取Sa及Sb对应的温度差,并取算术平均值,记作Ti,j,k(i=0,j=2,k=1);
3.8)重复步骤3.6)和步骤3.7),读取在该插入深度上,Sb分别在Sa距0点位置45°、90°、135°、180°、225°、270°、315°这7个位置时,Sa及Sb的温度差,并分别取算术平均值,记为Ti,j,k(i=0,j=2,k=1)至Ti,j,k(i=0,j=2,k=7);
3.9)将均温块插入深度调整为被测系统温区中心点向内深度5cm处,重复步骤3.3)至步骤3.8),读取在该插入深度上,Sb分别在Sa距0点位置45°、90°、135°、180°、225°、270°、315°这7个位置上,Sa及Sb的温度差,并分别取算术平均值,Ti,j,k(i=0,j=1,k=1)至Ti,j,k(i=0,j=1,k=7);
3.10)将均温块插入深度向外抽出5cm,重复步骤3.3)至步骤3.9),直到获取插入位置分别为中心点向内10cm,5cm、中心点,中心点向外5cm,10cm共5个插入深度上,每个插入深度8个不同角度共40个温度点的数据Ti,j,k(i=0,j=0,k=0)至Ti,j,k(i=0,j=4,k=7);
3.11)按照步骤3.1)至步骤3.2)调整设定温度点,在100℃至1100℃各温度点重复步骤3.10),获得数据Ti,j,k(i=0,j=0,k=0)至Ti,j,k(i=10,j=4,k=7)共440个数据。
进一步的,所述步骤3.2)中,调整补偿导线的接线端进行数据读取,读取方式为采用不大于3秒的采样间隔,连续读取1分钟相应数据。
本发明的有益效果是:
本发明是利用现有的热电偶校准系统硬件设施,使用本方法所述步骤完成实验数据采集并计算相关参数后,可得到一组不确定度计算公式。通过该公式,一方面可以完成对特定条件下的热电偶检定结果的不确定度的量化的计算,从而确定在该种条件下的热电偶校准结果的准确程度,扩展了热电偶校准过程中的条件限制。另一方面,也可以使用本公式对热电偶校准系统本身的性能进行评估。为不同的热电偶校准系统的横向性能比较提供一个新的测评维度。
附图说明
图1为均温块结构示意图;
图2为两个热电偶的连接示意图。
具体实施方式
现结合附图以及实施例对本发明做详细的说明。
实施例
本实施例提供给的热电偶校准系统测量不确定度的确定方法,包括以下步骤:
1)根据JJF1637、JJF1262标准选择符合I级的S型热电偶为Sa、Sb;
具体的,被评定热电偶校准系统满足下列要求:
首先,被评定的热电偶校准系统满足JJF1637、JJF1262中对“恒温设备”的一般技术要求;
其次,被评定的热电偶校准系统至少进行过10次以上的“室温-极限工作温度-室温”的升降温循环;
2)将上述Sa、Sb分别插入均温块的中心孔内,并反向串联Sa、Sb,形成热电偶校准系统;
具体的,均温块满足下列要求:
要求1:在插入热电偶校准系统核心控温区后,能通过延伸部分读取其在核心控温区的插入深度和轴向转动角度。读数允许误差要求为深度1mm,角度1°。均温块长度不小于20mm,半径不小于28mm。
要求2:在径向半径中心位置及离中心轴x(mm),位置处,具有至少一组以截面圆心呈角对称分布的热电偶埋入孔,x≥10mm深度10mm以上,埋入孔直径与所需设备1中所述的S型热电偶相匹配;
要求3:所用材料最高工作温度能够达到1200℃;
参见图1,a:中心电偶插入孔;b:离轴的电偶插入孔;实施时,均温块可通过自行加工获得,也可通过对被评估热电偶校准系统本身配置的均温块进行一定的改造获得;
本实施例中,将2支符合I级的S型热电偶,分别定义为Sa、Sb,在保证绝缘的情况下尽可能将Sa及Sb顶端的熔结球靠近,并插入均温块中心孔底部;将均温块插入深度调整到被测系统温区中心点处;
参见图2,Sa及Sb的电偶连接方式采用“反向串联”的方式连接,需注意的是,其外部的电偶引出线短接点,仍需与其补偿导线与Sa及Sb的连接一样浸入冰点槽中并接入补偿导线,以消除电偶本身不均匀性引起的短接点电势波动对实验的影响。
图2中,a为均温块;b为冰点槽;c:Sa熔结球;d:Sb熔结球;e:补偿导线正级;f:补偿导线负级;h:Sa偶的SN级;i:Sa偶的SP级;j:Sb偶的SN级;k:Sb偶的SP级;g:Sa的SN级与Sb的SN级的短接点;m:由g引出的补偿导线负极;
3)调整Sa、Sb插入中心孔中多个不同的深度、旋转均温块角度在多个不同位置处,以及多个不同的测温点条件下,Sa及Sb分别对应的温度差,并取算术平均值;得到不同插入深度j、不同旋转角度k以及不同测温点i时,热电偶校准系统对应的温度差,记为Ti,j,k(i=0,j=0,k=0)至Ti,j,k(i=10,j=4,k=7),共获得多组数据;
本实施例中,调整Sa、Sb插入中心孔中5个不同的深度、旋转均温块角度在7个不同位置处,以及10个不同的测温点条件下,Sa及Sb分别对应的温度差,并取算术平均值;得到不同插入深入j、不同旋转角度k以及不同测温点i时,热电偶校准系统对应的温度差,记为Ti,j,k(i=0,j=0,k=0)至Ti,j,k(i=10,j=4,k=7),共获得440组数据;
具体的,步骤3)的具体过程包括:
3.1)在100℃至1100℃之间,每100℃设置一个待检测温度点,在各待检测温度点读取Sa及Sb的温度差值;
实施时,将按照前述热电偶安装和连接的电偶引出补偿导线的e,m端分别连接到高精度万用表正负上,通过Sa偶测量均温块温度,当到达设定温度并且由Sa测得的温度变化小于0.1℃每分钟时;
3.2)若Sa的温度差值小于0.1℃每分钟,或者在任意温度点Sa及Sb的温度差值小于0.1℃,则进入步骤3.5);反之,则进行步骤3.3)~步骤3.4);
实施时,将万用表正负极分别连接到图2的e、f端。读取Sa及Sb在该温度点的差值。在该项测试中,若由Sa测得的温度变化无法达到小于0.1℃每分钟,或者在任意温度点Sa及Sb的温度差值不小于0.1℃,则需对待评估系统进行相关调整以满足上述要求后,方可进行后续步骤。具体为:如每分钟温度变化值超差,则调整检定系统的控温电偶位置及PID参数等相关设置,若Sa及Sb的温度差值超差,则调整熔结球位置或更换其他符合要求的电偶;
3.3)将Sa插入均温块的中心孔底部,Sb插入均温块任一离轴孔底部,通过旋转均温块,使Sb插入孔位于Sa插入孔的正上方0点位置;并将均温块插入深度调整到被测系统温区中心点向内深度10cm处;
3.4)按照步骤3.1)进行相关检测,直至Sa的温度差值小于0.1℃每分钟,或者在任意温度点Sa及Sb的温度差值小于0.1℃;
实施时,在温度点100℃进行相关测量,当符合条件时继续;
3.5)按照步骤3.2)读取的Sa及Sb温度差,并取算术平均值,记做Ti,j,k(i=0,j=2,k=0);
实施时,数据的读取方式为采用不大于3秒的采样间隔,连续读取1分钟相应数据,
3.6)将均温块顺时针旋转45°,插入深度不变,按照步骤3.1)在进行测量,直至Sa的温度差值小于0.1℃每分钟,或者在任意温度点Sa及Sb的温度差值小于0.1℃;
实施时,插入深度不变,按照步骤3.1)在温度点100℃进行相关测量,当符合条件时继续;
3.7)在步骤3.6)位置处,读取Sa及Sb对应的温度差,并取算术平均值,记作Ti,j,k(i=0,j=2,k=1);
3.8)重复步骤3.6)和步骤3.7),读取在该插入深度上,Sb分别在Sa距0点位置45°、90°、135°、180°、225°、270°、315°这7个位置上,Sa及Sb的温度差,并分别取算术平均值,记为Ti,j,k(i=0,j=2,k=1)至Ti,j,k(i=0,j=2,k=7);
3.9)将均温块插入深度调整为被测系统温区中心点向内深度5cm处,重复步骤3.3)至步骤3.8),读取在该插入深度上,Sb分别在Sa距0点位置45°、90°、135°、180°、225°、270°、315°这7个位置上,Sa及Sb的温度差,并分别取算术平均值,Ti,j,k(i=0,j=1,k=1)至Ti,j,k(i=0,j=1,k=7);
3.10)将均温块插入深度向外抽出5cm,重复步骤3.3)至步骤3.9),直到获取插入位置分别为中心点向内10cm,5cm、中心点,中心点向外5cm,10cm共5个插入深度上,每个插入深度8个不同角度共40个温度点的数据Ti,j,k(i=0,j=0,k=0)至Ti,j,k(i=0,j=4,k=7);
3.11)按照步骤3.1)至步骤3.2)调整设定温度点,在100℃至1100℃各温度点重复步骤3.10),获得数据Ti,j,k(i=0,j=0,k=0)至Ti,j,k(i=10,j=4,k=7)共440个数据;
本实施例中,旋转角度在距离0点的45°、90°、135°、180°、225°、270°、315°几个位置上变换;均温块插入深度在中心点向内10cm,中心点向内5cm、中心点、中心点向外5cm以及中心点向外10cm共5个点上;温度是100℃至1100℃之间,每100℃设置一个待检测温度点;但是实际测量时,还可根据实际情况将旋转角度、插入深度以及检测温度点进行调整或设计。
4)使用卡尺或设备测量均温块截面上插入热电偶的离轴孔的中心位置距均温块截面几何中心点的距离x;
5)根据步骤3)和步骤4)得到的数据,带入下列公式中,得到合成不确定度;
式中:μc为合成不确定度;
μ(e1)被测电偶输入量引入的不确定度分量;
c1被电偶输入量引入不确定度的灵敏系数;
μ(e2)标准热电偶本身引入的不确定度;
c2被测热电偶与标准热电偶在被检点上证书标注的的微分热电动势之比引入的灵敏系数;
μ(e3)检测时标准热电偶测量重复性引入的不确定度分量;
c3被测热电偶与标准热电偶在被检点上测得的微分热电动势算术平均值之比引入的灵敏系数;
μ(e4)使用补偿导线引入的不确定度分量;
c4补偿导线引入不确定度分量的灵敏系数;
式中,r实为测量均温块上离轴分布插入孔圆心所构成的圆的半径;h标实为标准热电偶偏离中心位置的长度;h被实为被检偶的深度偏移以一次测量过程中偏离中心位置最远的被测偶计算;
μ(e2)通过标准热电偶的证书获得;
式中α为补偿导线在测试时环境温度点的示值误差;
c4=1;
式中:e被(t)为被检电偶在某被校准温度点的热电动势;
被检电偶在上述温度点多次测量的算术平均值;
e标证标准热电偶在上述温度点时,证书给出的热电动势;
标准电偶在上述温度点多次测量的算术平均值;
e补补偿导线修正值;
S被被检热电偶在上述温度点的微分电动势;
S标标准热电偶在上述温度点的微分电动势。
本实施例中,在步骤5)的数据计算阶段该阶段主要对数据采集阶段采集记录的数据进行处理,以获被测量系统的得最终测量不确定度数值。
第一、电偶检测原理简述
根据校准热电偶时常用的比较法,通过比较被测偶电动势和标准偶电动势的差值,并通过查阅标准偶及被测偶在各个校准点的微分电动势值,通过相关计算获取被测偶的测量温度误差值。其数学模型参见如下公式:
上式中,e被(t)为被检电偶在某被校准温度点的热电动势;
被检电偶在上述温度点多次测量的算术平均值;
e标证标准热电偶在上述温度点时,证书给出的热电动势;
标准电偶在上述温度点多次测量的算术平均值;
e补补偿导线修正值;
S被被检热电偶在上述温度点的微分电动势;
S标标准热电偶在上述温度点的微分电动势;
第二、不确定度的确定过程
本实施例提供的不确定度的确定方法,通过一系列的技术手段确定了上述测量结果的示值误差测量不确定度,其数学模型参见如下公式:
式中:μc为合成不确定度;
μ(e1)被测电偶输入量引入的不确定度分量;
c1被电偶输入量引入不确定度的灵敏系数;
μ(e2)标准热电偶本身引入的不确定度;
c2被测热电偶与标准热电偶在被检点上证书标注的的微分热电动势之比引入的灵敏系数;
μ(e3)检测时标准热电偶测量重复性引入的不确定度分量;
c3被测热电偶与标准热电偶在被检点上测得的微分热电动势算术平均值之比引入的灵敏系数;
μ(e4)使用补偿导线引入的不确定度分量;
c4补偿导线引入不确定度分量的灵敏系数;
第三、关于各分量灵敏系数的计算(c1~c4):
特别说明的是,由于测量过程中采用的电测设备测量的是电偶的热电势,实际使用时根据查表换算为温度值,故在某些计算过程中,温度的比值等运算可直接通过对应的电势之比表示。
同种被测电偶在某温度的微分电动势值是一致的。
根据实际测量时的被测电偶和标准电偶种类确定,标准电偶通常采用S型电偶;
本实施例中,当被测电偶为K型时,在500℃校准点时,通过查表得
当被检温度点或被检电偶类型变化时,需进行查表和相关计算以确定该参数具体值。
类似c2的,该值根据实际测量时的被测电偶和标准电偶种类确定,其绝对值与c2相同,仅符号取反。当被检温度点或被检电偶类型变化时,需进行查表和相关计算以确定该参数具体值。
c4:被测电偶与补偿导线在被测温度点微分电动势之比,通常这两者的值非常接近,实际使用中取1即可;
最后,各不确定度分量的具体计算(μ(e1)~μ(e4)):
μ(e1):由于通过前述数据采集过程所采集到得数据,包含了测量重复性,电测仪器误差,炉温均匀性误差,转换开关触电电势,参考端误差等输入量影响量。仅考虑由于温场均匀性引入的不确定度。
由于实际测量时,采用捆扎或分布式插入均温块的方式布置被测电偶及标准电偶。实际操作时,首先测量能够包容所有被测及标准热电偶结球点的圆柱空间的直径。采用捆扎法时测量捆扎后的热电偶集束体半径。采用分布式插入的,测量均温块上离轴分布插入孔圆心所构成的圆的半径。该值记作r实。
实际测量时的电偶插入深度偏移通过实际测量获得,分别测量稳定插入后,标准热电偶偏离中心位置的长度h标实,h被实被检偶的深度偏移以一次测量过程中偏离中心位置最远的被测偶计算。
考虑到该不确定度分量主要由温场均匀性因素影响,符合均匀分布,带入包含因子则
式中:x为均温块离轴孔离轴距离
式子中T(i,j,k)参数i取与实际测试温度点最接近的温度点对应的i值。
μ(e2)则是通过标准热电偶的证书获得;
μ(e3):该分量通过数据采集阶段获得的数据求出,根据实际测量时的条件,选取距被校准温度点最近的温度点,确定Ti,j,k中的参数i,根据实际测量时标准电偶插入深度位置选择参数j,之后通过下公式算出合成样本标准偏差。
并根据实际测量时的测量次数确定结果,其中n为实际测量时为获取算术平均值所进行的测量次数,通过查阅系统说明书获得。
μ(e4):当实际测量过程中使用补偿导线对被测热电偶冷端进行处理时引入该分量,否则该分量取0。当使用补偿导线时,需查询相关技术资料,确定在测试进行时的环境温度下,所使用的补偿导线的证书误差值。由于该误差主要由于材料均匀度,补偿导线材料与电偶材料微分热电势差异等因素引起,符合均匀分布,故引入时需考虑包含因子k。具体为
式中:α为补偿导线在测试时环境温度点的示值误差,通过查阅补偿导线的说明书或证书获得。
综上所述,本发明利用现有的热电偶校准系统硬件设施,使用本发明提供的方法完成实验数据采集并计算相关参数后,可得到一组不确定度计算公式。通过该公式,一方面可以完成对特定条件下的热电偶检定结果的不确定度的量化的计算,从而确定在该种条件下的热电偶校准结果的准确程度,扩展了热电偶校准过程中的条件限制。
另一方面,也可以使用本公式对热电偶校准系统本身的性能进行评估。为不同的热电偶校准系统的横向性能比较提供一个新的测评维度。
Claims (6)
1.一种热电偶校准系统测量不确定度的确定方法,其特征在于,包括以下步骤:
1)根据JJF1637、JJF1262标准选择符合I级的S型热电偶为Sa、Sb;
2)将上述Sa、Sb分别插入均温块的中心孔内,并反向串联Sa、Sb,形成热电偶校准系统;
3)调整Sa、Sb插入中心孔中多个不同的深度、旋转均温块角度在多个不同位置处、以及多个不同的测温点条件下,Sa及Sb分别对应的温度差,并取算术平均值;得到不同插入深度j、不同旋转角度k以及不同测温点i时,热电偶校准系统对应的温度差,记为Ti,j,k(i=0,j=0,k=0)至Ti,j,k(i=10,j=4,k=7),共获得多组数据;
4)使用卡尺或设备测量均温块截面上插入热电偶的离轴孔的中心位置距均温块截面几何中心点的距离x;
5)根据步骤3)和步骤4)得到的数据,带入下列公式中,得到合成不确定度,
式中:μc为合成不确定度;
μ(e1)被测电偶输入量引入的不确定度分量;
c1被电偶输入量引入不确定度的灵敏系数;
μ(e2)标准热电偶本身引入的不确定度;
c2被测热电偶与标准热电偶在被检点上证书标注的的微分热电动势之比引入的灵敏系数;
μ(e3)检测时标准热电偶测量重复性引入的不确定度分量;
c3被测热电偶与标准热电偶在被检点上测得的微分热电动势算术平均值之比引入的灵敏系数;
μ(e4)使用补偿导线引入的不确定度分量;
c4补偿导线引入不确定度分量的灵敏系数;
所述
式中,r实为测量均温块上离轴分布插入孔圆心所构成的圆的半径;h标实为标准热电偶偏离中心位置的长度;h被为被检偶的深度偏移以一次测量过程中偏离中心位置最远的被测偶计算;
所述μ(e2)通过标准热电偶的证书获得;
所述
所述式中α为补偿导线在测试时环境温度点的示值误差;
所述c4=1;
所述
式中:n为实际测量时为获取算术平均值所进行的测量次数,通过查阅系统说明书获得;
e被(t)为被检电偶在某被校准温度点的热电动势;
被检电偶在上述温度点多次测量的算术平均值;
e标证标准热电偶在上述温度点时,证书给出的热电动势;
标准电偶在上述温度点多次测量的算术平均值;
e补补偿导线修正值;
S被被检热电偶在上述温度点的微分电动势;
S标标准热电偶在上述温度点的微分电动势。
2.根据权利要求1所述的热电偶校准系统测量不确定度的确定方法,其特征在于,所述步骤2)中,均温块长度不小于20mm,半径不小于28mm。
3.根据权利要求2所述的热电偶校准系统测量不确定度的确定方法,其特征在于,所述步骤2)中,在均温块的径向半径中心位置及离中心轴x(mm)位置处,具有至少一组以截面圆心呈角对称分布的热电偶埋入孔,x≥10mm深度10mm以上,埋入孔直径与步骤1)中的S型热电偶相匹配。
4.根据权利要求3所述的热电偶校准系统测量不确定度的确定方法,其特征在于,所述步骤2)中,反向串联时,均温块外部的热电偶引出线短接点,并接入补偿导线。
5.根据权利要求4所述的热电偶校准系统测量不确定度的确定方法,其特征在于,所述步骤3)的具体过程包括:
3.1)在100℃至1100℃之间,每100℃设置一个待检测温度点,在各待检测温度点读取Sa及Sb的温度差值;
3.2)若Sa的温度差值小于0.1℃每分钟,或者在任意温度点Sa及Sb的温度差值小于0.1℃,则进入步骤3.5);反之,则进行步骤3.3)~步骤3.4);
3.3)将Sa插入均温块的中心孔底部,Sb插入均温块任一离轴孔底部,通过旋转均温块,使Sb插入孔位于Sa插入孔的正上方0点位置;并将均温块插入深度调整到被测系统温区中心点向内深度10cm处;
3.4)按照步骤3.1)进行相关检测,直至Sa的温度差值小于0.1℃每分钟,或者在任意温度点Sa及Sb的温度差值小于0.1℃;
3.5)按照步骤3.2)读取的Sa及Sb温度差,并取算术平均值,记做Ti,j,k(i=0,j=2,k=0);
3.6)将均温块顺时针旋转45°,插入深度不变,按照步骤3.1)再次进行测量,直至Sa的温度差值小于0.1℃每分钟,或者在任意温度点Sa及Sb的温度差值小于0.1℃;
3.7)在步骤3.6)位置处,读取Sa及Sb对应的温度差,并取算术平均值,记作Ti,j,k(i=0,j=2,k=1);
3.8)重复步骤3.6)和步骤3.7),读取在该插入深度上,Sb分别在Sa距0点位置45°、90°、135°、180°、225°、270°、315°这7个位置时,Sa及Sb的温度差,并分别取算术平均值,记为Ti,j,k(i=0,j=2,k=1)至Ti,j,k(i=0,j=2,k=7);
3.9)将均温块插入深度调整为被测系统温区中心点向内深度5cm处,重复步骤3.3)至步骤3.8),读取在该插入深度上,Sb分别在Sa距0点位置45°、90°、135°、180°、225°、270°、315°这7个位置上,Sa及Sb的温度差,并分别取算术平均值,Ti,j,k(i=0,j=1,k=1)至Ti,j,k(i=0,j=1,k=7);
3.10)将均温块插入深度向外抽出5cm,重复步骤3.3)至步骤3.9),直到获取插入位置分别为中心点向内10cm,5cm、中心点,中心点向外5cm,10cm共5个插入深度上,每个插入深度8个不同角度共40个温度点的数据Ti,j,k(i=0,j=0,k=0)至Ti,j,k(i=0,j=4,k=7);
3.11)按照步骤3.1)至步骤3.2)调整设定温度点,在100℃至1100℃各温度点重复步骤3.10),获得数据Ti,j,k(i=0,j=0,k=0)至Ti,j,k(i=10,j=4,k=7)共440个数据。
6.根据权利要求5所述的热电偶校准系统测量不确定度的确定方法,其特征在于,所述步骤3.2)中,调整补偿导线的接线端进行数据读取,读取方式为采用不大于3秒的采样间隔,连续读取1分钟相应数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111275518.2A CN114354015B (zh) | 2021-10-29 | 2021-10-29 | 热电偶校准系统测量不确定度的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111275518.2A CN114354015B (zh) | 2021-10-29 | 2021-10-29 | 热电偶校准系统测量不确定度的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114354015A CN114354015A (zh) | 2022-04-15 |
CN114354015B true CN114354015B (zh) | 2024-05-24 |
Family
ID=81096294
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111275518.2A Active CN114354015B (zh) | 2021-10-29 | 2021-10-29 | 热电偶校准系统测量不确定度的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114354015B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20090132921A (ko) * | 2008-06-23 | 2009-12-31 | 한국표준과학연구원 | 자동 열전대 온도 교정 시스템 및 방법 |
JP2014149252A (ja) * | 2013-02-04 | 2014-08-21 | Yoshiro Mizuno | 温度計の校正データの生成方法、校正データを格納した記憶装置、及びこの方法を採用した温度計 |
CN108680285A (zh) * | 2018-08-14 | 2018-10-19 | 上海市计量测试技术研究院 | 短支热电偶温度校验炉及短支热电偶的校准方法 |
CN109855583A (zh) * | 2018-11-16 | 2019-06-07 | 中国航发西安动力控制科技有限公司 | 三坐标测量机测量不确定度的确定方法 |
CN110595645A (zh) * | 2019-09-02 | 2019-12-20 | 天津市计量监督检测科学研究院 | 一种短型廉金属热电偶校准装置及其校准方法 |
CN111024265A (zh) * | 2019-11-26 | 2020-04-17 | 山东磐然仪器集团有限公司 | 一种管式炉进行热电偶标准的稳定性考核方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
MX340152B (es) * | 2013-12-19 | 2016-06-08 | Kaplun Mucharrafille Margarita | Sistema y metodo para calibracion y/o caracterizacion de instrumentos de medicion de temperatura por telemetria. |
-
2021
- 2021-10-29 CN CN202111275518.2A patent/CN114354015B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20090132921A (ko) * | 2008-06-23 | 2009-12-31 | 한국표준과학연구원 | 자동 열전대 온도 교정 시스템 및 방법 |
JP2014149252A (ja) * | 2013-02-04 | 2014-08-21 | Yoshiro Mizuno | 温度計の校正データの生成方法、校正データを格納した記憶装置、及びこの方法を採用した温度計 |
CN108680285A (zh) * | 2018-08-14 | 2018-10-19 | 上海市计量测试技术研究院 | 短支热电偶温度校验炉及短支热电偶的校准方法 |
CN109855583A (zh) * | 2018-11-16 | 2019-06-07 | 中国航发西安动力控制科技有限公司 | 三坐标测量机测量不确定度的确定方法 |
CN110595645A (zh) * | 2019-09-02 | 2019-12-20 | 天津市计量监督检测科学研究院 | 一种短型廉金属热电偶校准装置及其校准方法 |
CN111024265A (zh) * | 2019-11-26 | 2020-04-17 | 山东磐然仪器集团有限公司 | 一种管式炉进行热电偶标准的稳定性考核方法 |
Non-Patent Citations (5)
Title |
---|
Evaluating measurement uncertainty of thermocouples calibration;Sârbu, GC等;2018 INTERNATIONAL CONFERENCE AND EXPOSITION ON ELECTRICAL AND POWER ENGINEERING (EPE);20190301;第1-5页 * |
微波环境中热电偶的温度测量研究(英文);李琳 等;华南理工大学学报(自然科学版);20000328;第54-58页 * |
热电偶检定系统测量示值误差的不确定度评定;陆凤波;;中国高新技术企业;20080215(第04期);第99-100页 * |
热电偶自动检定系统测量结果的不确定评定;穆允宁;;山东冶金;20191020(第05期);第49-51页 * |
短型热电偶校准方法及测量不确定度分析探讨;杨金波 等;仪器仪表标准化与计量;20211025;第40-42页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114354015A (zh) | 2022-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105628732B (zh) | 一种测量Seebeck系数的装置及其方法 | |
CN102252770A (zh) | 一种测温补偿方法及修正型高精度温度计 | |
CN106092375B (zh) | 机载设备地面温度传感器的校验方法及校验仪器 | |
CN116678489B (zh) | 力传感器的质量控制方法、系统、设备及介质 | |
CN108152325B (zh) | 一种基于防护热板法校准热流计法导热仪的方法 | |
CN108680285A (zh) | 短支热电偶温度校验炉及短支热电偶的校准方法 | |
CN110595645B (zh) | 一种短型廉金属热电偶校准装置及其校准方法 | |
CN114487024B (zh) | 一种基于幂函数的钯合金氢气传感器的校准拟合方法 | |
CN114879775B (zh) | 一种温度控制方法、装置、系统、设备及存储介质 | |
CN111397772A (zh) | 异型热电偶的检定方法 | |
CN114354015B (zh) | 热电偶校准系统测量不确定度的确定方法 | |
CN111024265A (zh) | 一种管式炉进行热电偶标准的稳定性考核方法 | |
CN108519406B (zh) | 一种架空地线轴向的导体热阻和对流热阻的计算方法 | |
CN109596890A (zh) | 一种电阻温度系数测试仪 | |
CN113075571A (zh) | 一种锂离子电池ocv确定方法、装置及系统 | |
CN112098457B (zh) | 一种导热系数测量仪的多项式回归校准方法 | |
CN112903145B (zh) | 高低温探针台在片温度参数的校准方法及装置 | |
CN110715756A (zh) | 一种负载状况下管式炉温度场性能测试装置及其性能测试方法 | |
CN113960256A (zh) | 一种含水仪的温度补偿方法 | |
CN114002378A (zh) | 气体浓度传感器的浓度检测方法 | |
CN110118630B (zh) | 一种磁控真空计校准判断方法 | |
Cash et al. | Application of a digital computer to data acquisition and shield temperature control of a high‐temperature, adiabatic calorimeter | |
CN109060190A (zh) | 一种钢水测温定氧系统的校准方法 | |
Drnovsek et al. | Reduction of uncertainties in temperature calibrations by comparison | |
CN115541043B (zh) | 一种高温柜样品盒的测温装置及其测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |