JJF 1049-2024 温度传感器动态响应校准规范

中华人民共和国国家计量技术规范
JJF1049—2024
温度传感器动态响应校准规范
CalibrationSpecificationforTemperatureSensors’DynamicResponse
2024-09-18发布2025-03-18实施
国家市场监督管理总局 发布
温度传感器动态响应校准规范
CalibrationSpecificationforTemperature
Sensors’DynamicResponse
?????????????
?????
?????????????
??????
?
?
?
JJF1049—2024
代替JJF1049—1995
归口单位:全国温度计量技术委员会
主要起草单位:中国航空工业集团公司北京长城计量测试技术
研究所
上海仪器仪表自控系统检验测试所有限公司
参加起草单位:中国计量科学研究院
沈阳东大传感技术有限公司
本规范委托全国温度计量技术委员会负责解释
JJF1049—2024
本规范主要起草人:
赵 俭(中国航空工业集团公司北京长城计量测试技术
研究所)
胡林陶(中国航空工业集团公司北京长城计量测试技术
研究所)
郭爱华(上海仪器仪表自控系统检验测试所有限公司)
参加起草人:
廖 理(中国航空工业集团公司北京长城计量测试技术
研究所)
郑 玮(中国计量科学研究院)
王魁汉(沈阳东大传感技术有限公司)
JJF1049—2024
目 录
引言……………………………………………………………………………………… (Ⅱ)
1 范围…………………………………………………………………………………… (1)
2 引用文件……………………………………………………………………………… (1)
3 术语和定义…………………………………………………………………………… (1)
4 概述…………………………………………………………………………………… (1)
4.1 结构和原理………………………………………………………………………… (1)
4.2 用途………………………………………………………………………………… (2)
5 计量特性……………………………………………………………………………… (2)
6 校准条件……………………………………………………………………………… (2)
6.1 实验室环境条件…………………………………………………………………… (2)
6.2 计量标准及其他设备……………………………………………………………… (2)
7 校准项目和校准方法………………………………………………………………… (5)
7.1 校准项目…………………………………………………………………………… (5)
7.2 校准方法…………………………………………………………………………… (5)
8 校准结果表达………………………………………………………………………… (9)
9 复校时间间隔………………………………………………………………………… (10)
附录A 热响应时间作图法计算……………………………………………………… (11)
附录B 校准结果的测量不确定度评定……………………………………………… (13)
附录C 校准结果原始记录…………………………………………………………… (23)
附录D 校准证书格式………………………………………………………………… (24)
Ⅰ
JJF1049—2024
引 言
JJF1071 《国家计量校准规范编写规则》、JJF1001 《通用计量术语及定义》和
JJF1059.1《测量不确定度评定与表示》共同构成支撑本规范编写工作的基础性系列
规范。
本规范代替JJF1049—1995 《温度传感器动态响应校准规范》。与JJF1049—1995
相比,除编辑性修改外,本规范主要技术变化如下:
———增加了以激光器为动态温度激励源的校准方式;
———修改了水流环境为液流环境(见6.2.1);
———增加了测量不确定度评定内容(见附录B)。
本规范的历次版本发布情况为:
———JJF1049—1995。
Ⅱ
JJF1049—2024
温度传感器动态响应校准规范
1 范围
本规范适用于时间常数不小于1ms的接触式温度传感器在室温~1700 ℃、气流
马赫数0.95以下或液体流速1m/s以下的动态响应(热响应时间)校准,温度传感器
主要包括热电偶、热电阻式温度传感器。
2 引用文件
本规范无引用文件。
3 术语和定义
3.1 热响应时间 thermalresponsetime
当温度出现阶跃变化时,温度传感器的输出量变化到相当于该温度阶跃量的某个规
定百分数所需要的时间。变化到温度阶跃量的10%、50%、90%所需要的时间,分别
记为τ0.1、τ0.5、τ0.9。
3.2 时间常数 timeconstant
温度传感器的输出量变化到相当于温度阶跃量63.2%的热响应时间。符号为τ。
3.3 温度阶跃量 temperaturestep
温度阶跃变化之后与之前稳定状态温度的差值。
3.4 动态温度激励源 dynamictemperatureexcitationsource
用于产生温度变化的系统或装置。
3.5 弹射机构 ejectionmechanism
采用机械方式完成弹射,最终实现温度阶跃变化的机构。
3.6 激光激励源 laserexcitationmechanism
产生不同幅值和脉宽激光,实现温度阶跃(或脉冲)变化的装置。
3.7 气流总温 totaltemperatureofgasflow
气流绝热等熵滞止到速度为零时的温度。
3.8 气流总压 totalpressureofgasflow
气流绝热等熵滞止到速度为零时的压力。
4 概述
4.1 结构和原理
温度传感器通常由敏感元件、外壳、引线等组成,典型结构示意图如图1所示。
1
JJF1049—2024
图1 温度传感器典型结构示意图
由于温度传感器具有热惯性,在测量快速变化的温度时,不能立刻反映被测温度,
而是需要一定的时间才能达到热平衡状态。温度传感器的动态响应特性指温度传感器的
温度与被测介质温度变化量之间的关系,一般用微分方程或传递函数的形式表示。动态
响应特性的特征参数主要是热响应时间(含时间常数)。温度传感器的热响应时间不仅
与传感器本身的结构、材料有关,而且与流体介质、工况条件有关。
4.2 用途
温度传感器主要用于测量变化的气体、液体等流体温度。温度传感器的动态响应特
性,主要用于评价温度传感器的响应快慢。
5 计量特性
温度传感器动态响应计量特性见表1。
表1 温度传感器动态响应计量特性
序号计量特性名称说明
1 时间常数
2 热响应时间(τ0.1、τ0.5、τ0.9)
并非每类温度传感器均要求测量本表所列全部计量
特性,因此校准项目依据温度传感器实际情况及客
户要求确定
6 校准条件
6.1 实验室环境条件
实验室环境要求见表2。
表2 实验室环境要求
环境参数要求
温度15℃~25℃
相对湿度≤80%
6.2 计量标准及其他设备
校准所用的计量设备与仪器应经过计量技术机构检定(或校准),满足校准使用要
2
JJF1049—2024
求,并在有效期内。
6.2.1 计量标准
液流槽或校准风洞系统、动态温度激励源和高速采集设备共同构成计量标准,详见
表3。
表3 计量标准
序号计量标准提供的校准环境
1 液流槽,介质可以是水或油液流环境
2 高速校准风洞高速气流环境,气流马赫数0.1~0.95
3 低速校准风洞低速气流环境,气流速度5m/s~40m/s
4 动态温度激励源,可以
是弹射机构或激光激励源———
5 高速采集设备———
计量标准的技术要求见表4。
表4 计量标准的技术要求
序号计量标准技术要求用途
1 液流槽
液流温度
波动度在10min内优于温度阶跃量的±1%
均匀度
在被校传感器插入方向优于温度阶跃
量的±2%
液流速度
控制偏差在测量位置优于±0.05m/s
波动度在10min内优于液流速度的±5%
均匀度
在被校传感器插入方向优于液流速度
的±10%
液流流道
宽度不小于被校传感器直径的10倍
深度不小于100mm
液流环境
下校准的
计量标准
2 高速校
准风洞
试验段口径堵塞比不大于10%
气流马赫数
均匀度
在流道中心40%截面内优于气流马
赫数的±5%
控制偏差优于±0.005
气流方向角偏差在流道中心40%截面内优于±0.5°
气流温度波动度
在常温高速校准风洞上校准时,在
1min内优于±0.3 ℃;在热校准风
洞上校准,900 ℃以下时,在1min
内优于±1℃,900 ℃~1700 ℃时,
在1min内优于总温的±1%
高速(气流
马赫数0.1~
0.95)气流
环境下校
准的计量
标准
3
JJF1049—2024
表4 (续)
序号计量标准技术要求用途
2 高速校
准风洞
气流温度
均匀度
在常温高速校准风洞上校准时,在流
道中心40%截面内优于±0.5℃;在
热校准风洞上校准时,在流道中心
40%截面内优于总温的±5%,在流
道中心10%截面内优于总温的±1%
控制偏差900℃以下时,优于±10℃;900℃~
1700℃时,优于±20℃
高速(气流
马赫数0.1~
0.95)气流
环境下校
准的计量
标准
3 低速校
准风洞
气流速度
波动度在10min内优于气流速度的±5%
均匀度
在流道中心40%截面内优于气流速
度的±5%
气流温度
波动度在10min内优于温度阶跃量的±1%
均匀度
在流道中心40%截面内优于温度阶
跃量的±1%
气流流道宽度不小于被校传感器直径的10倍
低速(气流
速度5m/s~
40m/s)
气流环境
下校准的
计量标准
4 动态温度
激励源
执行时间
小于被校温度传感器时间常数
的10%
激励幅度不小于10℃
实现温度阶跃
5 高速采
集设备
电压
测量范围1mV~30V
最大允许
误差±0.5%
分辨力不低于16位
时间间隔
最大允许
误差
±(0.06×采样间隔+1×10-6的测量
间隔)
采样频率上限不低于1MHz
温度动态
信号采集,
建议使用
示波器
6.2.2 其他设备
其他设备与仪器及其技术要求见表5。
4
JJF1049—2024
表5 其他设备与仪器及其技术要求
序号设备或仪器测量范围技术要求用途
1 压力传感器0kPa~100kPa 最大允许误差:±0.15% 气流总压测量
2 总压探针0kPa~100kPa 总压系数:优于0.005 气流总压测量
3 压力传感器-20kPa~20kPa 最大允许误差:±0.15% 气流静压测量
4 大气压力计95kPa~105kPa 最大允许误差:±40Pa 大气压测量
5 风速计5m/s~40m/s 最大允许误差:
±(0.03m/s+5%×读数)
低速校准风洞的
气流速度测量
6 液体流速仪0.02m/s~2m/s 最大允许误差:
±(1%×读数+0.01m/s) 液流速度测量
7 冰点恒温器0℃ 温度:(0±0.1)℃ 热电偶冷端补偿
8 数据采集
处理系统
0V~1V,
0kΩ~5kΩ
对于电压测量,最大允许误差
为±(0.015%×读数+
0.005%×量程+3μV);
对于电阻测量,最大允许
误差为±(0.015%×读数+
0.0075%×量程+1.5mΩ)
采集温度传感器、
压力传感器的信号
7 校准项目和校准方法
7.1 校准项目
时间常数τ 和热响应时间τ0.1、τ0.5、τ0.9。
7.2 校准方法
温度传感器动态响应校准过程主要包括:产生稳定的校准工况(稳定的速度场、温
度场);使温度传感器接受动态温度激励,由高速采集设备(如示波器)采集温度传感
器的响应信号;计算出时间常数或热响应时间。
7.2.1 液流环境下的校准
7.2.1.1 校准装置准备
a)对于测量滑油、燃油、冷却液等液流温度的传感器,在液流槽中进行校准。
b)根据校准要求,选择水或油等合适的介质。
7.2.1.2 被校传感器的安装连接
a)将被校温度传感器装入弹射机构;
b)将被校温度传感器连接至高速采集设备,校准示意图如图2所示;
c)校准前,校准仪器充分预热。
5
JJF1049—2024
图2 温度传感器在液流环境下(正阶跃)的动态响应校准示意图
7.2.1.3 校准工况的选择
校准温度应选择在温度传感器的工作温度附近,初始时被校温度传感器在空气中,
阶跃量为10 ℃~270 ℃,建议至少在3个液流速度条件下进行校准,推荐液流速度
0.4m/s、0.7m/s、1.0m/s。也可根据用户要求,选择其他工况。
7.2.1.4 校准试验步骤
a)开启并调节液流槽,使装置内液流温度、液流速度达到要求值。
b)记录被校温度传感器的初始状态,包括液流温度、液流速度。
c)高速采集设备开始采集数据,启动弹射机构,使被校温度传感器迅速投入液流
槽介质中,记录被校温度传感器的温度与时间的变化关系,采集记录时间应在被校温度
传感器时间常数τ 预估值的10倍以上,如数据处理后发现采集记录时间达不到时间常
数τ 实际测量值的10倍以上,需要延长采集记录时间至时间常数τ 实际测量值的10倍
以上,重新校准。采样间隔应在时间常数τ 预估值的0.1%以下。
d)完成温度响应记录后,记录终止状态(完成温度阶跃后的状态),包括液流温
度、液流速度。
e)重复步骤a)~d),共测量3次。3次的测量值与其算术平均值的偏差均不得大
于10%,否则需要重复操作,直至满足要求为止。
7.2.1.5 数据处理
由采集的温度-时间响应曲线,经滤波后,直接用作图法算出时间常数或热响应时
间,详见附录A。
7.2.2 高速气流环境
7.2.2.1 校准装置准备
a)高速气流环境下的动态响应校准,在校准风洞上进行。
b)根据校准要求,选择合适的校准风洞。
c)将总压探针置于稳定段内,用于测量气流总压p0,用引压管连接总压探针和压
力传感器。
d)用引压管连接静压孔和压力传感器,以测量气流静压ps;对于开口试验段,也
可用大气压代表气流静压。
7.2.2.2 被校传感器的安装连接
6
JJF1049—2024
a)将被校温度传感器固定在试验段的安装座上。
b)将被校温度传感器连接至高速采集设备,校准示意图分别如图3、图4所示。
c)校准前,校准仪器和稳压电源充分预热。
图3 温度传感器在室温高速气流环境下(负阶跃)的动态响应校准示意图
图4 温度传感器在高温气流环境下(正阶跃)的动态响应校准示意图
7.2.2.3 校准工况的选择
通常采用温度传感器的使用工况或尽量接近其使用工况。为便于温度传感器之间相
互比较,推荐采用表6、表7中规定的工况,建议至少在3个流速(马赫数)条件下进
行校准。也可根据用户要求,选择其他工况。
表6 传感器使用温度范围与推荐的校准温度
传感器使用温度/℃ 室温~400 400~1000 1000~1700
校准温度/℃ 室温或200±10 600±10 1400±50
阶跃温度/℃ 100±20 200±20 200±20
7
JJF1049—2024
表7 传感器使用马赫数范围与推荐的校准马赫数
使用马赫数0.1~0.3 使用马赫数0.3~0.6 使用马赫数0.6~0.95
校准马赫数0.1、0.2、0.3 校准马赫数0.3、0.4、0.6 校准马赫数0.6、0.8、0.95
7.2.2.4 校准试验步骤
a)开启并调节风洞系统,使试验段的气流速度(马赫数)和气流温度达到要求值。
b)通过冷气流和热气流按不同比例掺混的方式,用动态温度激励源改变被校温度
传感器的温度,调整到所需的温度阶跃量,记录初始状态,包括气流总压、静压、大气
压、总温和被校温度传感器温度等。
c)高速采集设备开始采集数据,启动动态温度激励源,使被校温度传感器周围产
生温度突变,记录被校温度传感器的温度与时间的变化关系,采集记录时间应在被校温
度传感器时间常数τ 预估值的10倍以上。如数据处理后发现采集记录时间达不到时间
常数τ 实际测量值的10倍以上,需要延长采集记录时间至时间常数τ 实际测量值的
10倍以上,重新校准。采样间隔应在时间常数τ 预估值的0.1%以下。
d)完成温度响应记录后,记录终止状态,包括气流总压、静压、大气压、总温和
被校温度传感器温度等。
e)重复步骤a)~d),共测量3次。3次的测量值与其算术平均值的偏差均不得大
于10%,否则需要重复操作,直至满足要求为止。
7.2.2.5 数据处理
a)按式(1)计算气流马赫数Ma。
Ma= 2 κ-1
p0
ps κ-1 κ -1 ? ? ??
? ? ?? (
1)
式中:
p0———气流总压,Pa;
ps———气流静压,Pa;
κ ———绝热指数。
b)由采集的温度-时间响应曲线,经滤波后,直接用作图法算出时间常数或热响应
时间,详见附录A。
7.2.3 低速气流环境
7.2.3.1 校准装置准备
a)低速气流环境下的动态响应校准,在低速校准风洞上进行。
b)将风速计置于试验段内,用于测量气流速度。
7.2.3.2 被校传感器的安装连接
a)将被校温度传感器装入试验段内。
b)将被校温度传感器连接至高速采集设备,校准示意图如图5所示。
c)校准前,校准仪器充分预热。
8
JJF1049—2024
图5 温度传感器在低速气流环境下(负阶跃)的动态响应校准示意图
7.2.3.3 校准工况的选择
通常采用温度传感器的使用工况或尽量接近其使用工况。阶跃量为10 ℃~20 ℃,
建议至少在3个流速条件下进行校准,为便于温度传感器之间相互比较,推荐流速
5m/s、20m/s、30m/s。也可根据用户要求,选择其他工况。
7.2.3.4 校准试验步骤
a)开启并调节风洞系统,调整气流速度至校准点。
b)通过冷气流和热气流按不同比例掺混的方式,用动态温度激励源改变被校温度
传感器的温度,调整到所需的温度阶跃量,记录初始状态,包括气流速度和被校温度传
感器温度等。
c)高速采集设备开始采集数据,启动动态温度激励源,使被校温度传感器周围产
生温度突变,记录被校温度传感器的温度与时间的变化关系,采集记录时间应在被校温
度传感器时间常数τ 预估值的10倍以上。如数据处理后发现采集记录时间达不到时间
常数τ 实际测量值的10倍以上,需要延长采集记录时间至时间常数τ 实际测量值的
10倍以上,重新校准。采样间隔应在时间常数τ 预估值的0.1%以下。
d)完成温度响应记录后,再记录一次终止状态,包括气流速度和被校温度传感器
温度等。
e)重复步骤a)~d),共测量3次。3次的测量值与其算术平均值的偏差均不得大
于10%,否则需要重复操作,直至满足要求为止。
7.2.3.5 数据处理
由采集的温度-时间响应曲线,经滤波后,直接用作图法算出时间常数或热响应时
间,详见附录A。
8 校准结果表达
校准结果应在校准证书上反映。校准证书应至少包括以下信息:
a)标题:“校准证书”;
b)实验室名称和地址;
c)进行校准的地点(如果与实验室的地址不同);
d)证书的唯一性标识(如编号),每页及总页数的标识;
e)客户的名称和地址;
9
JJF1049—2024
f)被校对象的描述和明确标识;
g)进行校准的日期,如果与校准结果的有效性和应用有关时,应说明被校对象的
接收日期;
h)如果与校准结果的有效性和应用有关时,应对被校样品的抽样程序进行说明;
i)校准所依据的技术规范的标识,包括名称及代号;
j)本次校准所用计量标准的溯源性及有效性说明;
k)校准环境的描述;
l)校准结果及其测量不确定度的说明;
m)对校准规范的偏离的说明;
n)校准证书或校准报告签发人的签名、职务或等效标识;
o)校准结果仅对校准对象有效的声明;
p)未经实验室书面批准,不得部分复制证书的声明。
9 复校时间间隔
复校时间间隔建议为12个月。
由于复校时间间隔的长短是由仪器的使用情况、使用者、仪器本身质量等诸因素所
决定的,因此送校单位可根据实际使用情况自主决定复校时间间隔。
10
JJF1049—2024
附录A
热响应时间作图法计算
以某K型热电偶温度传感器为例,说明热响应时间的计算方法。校准工况:气流
总温为600℃,温度阶跃量为200℃,气流马赫数为0.2。
温度、时间、热电势和温度阶跃量分别用T 、t、E 和ΔT 表示,下标A、B、C均
表示曲线上的点位。
高速采集设备记录的波形是热电势随时间变化的曲线,如图A.1所示,其中A 为
阶跃前曲线水平段的最末端,选择A 作为阶跃的起始点,C为阶跃后曲线水平段的最
前端,选择C作为阶跃的终止点。
图A.1 热电势的时间响应曲线
温度阶跃前后的热电势见表A.1。
表A.1 阶跃前后的热电势
气流总温/℃ 气流马赫数
阶跃前热电势
EA/mV
阶跃后热电势
EC/mV
600 0.2 16.385 24.952
按照K型热电偶的参考函数,将热电势转换成温度,见表A.2。
表A.2 阶跃前后的温度
气流总温/℃ 气流马赫数
阶跃前温度
TA/℃
阶跃后温度
TC/℃
600 0.2 399.7 601.1
11
JJF1049—2024
相应的,温度随时间变化的曲线如图A.2所示。
图A.2 温度的时间响应曲线
温度阶跃量为:
ΔT =TC-TA=601.1℃-399.7℃=201.4℃
对应的曲线上B点温度为:
TB=TA+0.632ΔT =399.7℃+201.4℃×0.632=527.0℃
温度从399.7℃变化到527.0℃所用的时间,即B点的横坐标与A 点的横坐标之
差为:
tB-tA=2.36s-1.33s=1.03s
则该传感器的时间常数为:
τ=1.03s
其他热响应时间的计算方法与此类似,可参照进行。
激光脉冲激励的情况,在响应曲线的下降沿上分析,计算方法与此类似,可参照
进行。
12
JJF1049—2024
附录B
校准结果的测量不确定度评定
B.1 液流环境下时间常数τ 的测量不确定度评定
B.1.1 测量模型
τ=tB-tA (B.1)
式中:
tA———阶跃起点A对应的时间,s;
tB———响应曲线上B点对应的时间,s。
其中:
tB=t(TB)=t[xTC+(1-x)TA] (B.2)
式中:
TB———响应曲线上B点对应的温度,℃;
x ———温度变化到温度阶跃量的百分比,时间常数对应的x=63.2%;
TC———阶跃终点C对应的温度,℃;
TA———阶跃起点A对应的温度,℃。
液流速度v 为时间常数的强影响量,则有:
uc(τ)= u2p+u2(tB)+u2(tA)+ ∂τ
∂vu(v) ? ? ??
? ? ??2
(
B.3)
式中:
up ———测量重复性引入的不确定度;
∂τ
∂v ———时间常数对液流速度的偏导数;
v ———液流速度,m/s。
TA、TC由同一支传感器测得,因而可按完全正相关考虑,则有:
u(TB)=u[xTC+(1-x)TA]
=xu(TA)+(1-x)u(TA)=u(TA) (B.4)
因此,
u(tB)=∂tB
∂TB
u(TB)=∂tB
∂TB
u(TA) (B.5)
B.1.2 测量不确定度的来源
a)测量重复性引入的不确定度up。
b)u(TA)的来源如下:
1)温度传感器示值误差引入的不确定度;
2)阶跃起点A选取引入的不确定度;
3)数采系统测量误差引入的不确定度;
4)液流温度波动引入的不确定度。
13
JJF1049—2024
c)u(tA)的来源如下:
1)数据采样间隔引入的不确定度;
2)阶跃起点A选取引入的不确定度。
d)u(v)的来源如下:
1)流速测量仪的示值误差引入的不确定度;
2)液流速度波动引入的不确定度。
B.1.3 测量不确定度的评定
以某热电偶式温度传感器在速度为0.4 m/s的水流环境中、温度从室温阶跃到
80℃的时间常数校准为例,进行不确定度评定。
B.1.3.1 up的评定如下:
用极差法计算实验标准偏差s(τ)。在给定工况条件下共测量3次,3次时间常数τ
的测量结果分别为8.653s、8.524s、8.618s。从测量数据中找出最大值8.653s和最
小值8.524s,根据测量次数查表得到系数dn=1.693,按下式计算实验标准偏差s(τ):
s(τ)=τmax-τmin
dn =0.076s (B.6)
采用3次测量结果的平均值8.598s作为测量结果,测量重复性引入的不确定度up
按下式计算:
up=s(τ)
3 ≈0.044s (B.7)
B.1.3.2 u(TA)的评定如下:
a)温度传感器示值误差引入的不确定度
根据温度传感器的准确度等级,其最大允许误差为2.5℃,并经校准合格。由此引
入的不确定度按B类评定,假设测量值在允许误差范围内的概率分布为均匀分布,查
得包含因子k 为3,则其标准不确定度为:
u1=2.5℃
3 ≈1.4℃ (B.8)
b)阶跃起点A选取引入的不确定度
根据经验,由阶跃起点选取不准带来的温度误差为0.2℃。由此引入的不确定度按
B类评定,假设测量值在允许误差范围内的概率分布为均匀分布,查得包含因子k 为
3,则其标准不确定度为:
u2=0.2℃
3 ≈0.12℃ (B.9)
c)数采系统测量误差引入的不确定度
根据高速采集设备的准确度等级,其相对允许误差极限为±0.5%,并经校准合格。
由此引入的不确定度按B类评定,假设测量值在允许误差范围内的概率分布为均匀分
布,查得包含因子k 为3,则其标准不确定度为:
u3=0.5%×25℃
3 ≈0.07℃ (B.10)
14
JJF1049—2024
d)液流温度波动引入的不确定度
按B类评定,根据液流槽的校准结果,装置液流温度波动值不超过温度阶跃量的
±1%。假设测量值在允许误差范围内的概率分布为均匀分布,包含因子k 为3,则其
标准不确定度为:
u4=1%
3×55℃≈0.32℃ (B.11)
则有:
u(TA)= u2
1+u2
2+u2
3+u2
4 =1.4℃ (B.12)
B.1.3.3 u(tA)的评定如下:
a)数据采样间隔引入的不确定度
假定测量时设定的数据采样间隔为2ms。由此引入的不确定度按B类评定,假设
测量值在允许误差范围内的概率分布为均匀分布,查得包含因子k 为3,则其标准不
确定度为:
u5=0.001s
3 ≈0.0006s (B.13)
b)阶跃起点A选取引入的不确定度
根据经验,由阶跃起点选取不准带来的时间常数误差为±0.2s。由此引入的不确
定度按B类评定,假设测量值在允许误差范围内的概率分布为均匀分布,查得包含因
子k 为3,则其标准不确定度为:
u6=0.2s
3 ≈0.12s (B.14)
则有:
u(tA)= u2
5+u2
6 =0.12s (B.15)
B.1.3.4 u(v)的评定如下:
a)流速测量仪的示值误差引入的不确定度
根据流速测量仪的准确度等级,其最大允许误差为±(1%×读数+0.01m/s),并
经校准合格。由此引入的不确定度按B类评定,假设测量值在允许误差范围内的概率
分布为均匀分布,查得包含因子k 为3,则其标准不确定度为:
u7=1%×0.4m/s+0.01m/s
3 ≈0.008m/s (B.16)
b)液流速度波动引入的不确定度
按B类评定,根据液流槽的校准结果,装置的液流速度相对波动值不超过±5%。
假设测量值在允许误差范围内的概率分布为均匀分布,包含因子k 为3,则其标准不
确定度为:
u8=5%×0.4m/s
3 ≈0.011m/s (B.17)
c)计算u(v)
15
JJF1049—2024
u(v)= u27+u28 =0.014m/s (B.18)
B.1.4 相对合成标准不确定度ur(τ)
表B.1 不确定度分量一览表
不确定度分量不确定度来源评定方法分布k 值
标准不
确定度
up 测量的重复性A ——— ——— 0.044s
u1 温度传感器示值误差B 均匀3 1.4℃
u2 阶跃起点A 选取误差B 均匀3 0.12℃
u3 数采系统测量误差B 均匀3 0.07℃
u4 液流温度波动B 均匀3 0.32℃
u5 数据采样间隔B 均匀3 0.0006s
u6 阶跃起点A 选取误差B 均匀3 0.12s
u7 流速测量仪的示值误差B 均匀3 0.008m/s
u8 液流速度波动B 均匀3 0.011m/s
不确定度分量一览表见表B.1。其中各不确定度分量互不相关,则时间常数的相对
合成标准不确定度ur(τ)按下式计算:
ur(τ)=
u2p +u2(tB)+u2(tA)+ ∂τ
∂vu(v) ? ? ??
? ? ??
2
τ
=
u2p + ∂tB
∂TB
u(TA) ? ? ??
? ? ??
2
+ u(tA) 2+ ∂τ
∂vu(v) ? ? ??
? ? ??
2
τ
=16% (B.19)
其中,∂tB
∂TB
具体计算方法为:以温度为横坐标,时间为纵坐标,通过对TA、TB、
TC对应的时间tA、tB、tC进行多项式拟合,对拟合的多项式求导,求得B点处的切线
斜率,计算得到∂tB
∂TB=1。∂τ
∂v具体计算方法为:以流速为横坐标,时间常数为纵坐标,
通过对不同流速下对应时间常数点进行多项式拟合,对拟合的多项式求导,求得对应流
速点的切线斜率,计算得到∂τ
∂v=1。
B.1.5 相对扩展不确定度Ur(τ)
取k=2,则相对扩展不确定度Ur(τ)为
Ur(τ)=k·ur(τ)=32% (B.20)
16
JJF1049—2024
B.2 高速气流环境下时间常数τ 的测量不确定度评定
B.2.1 测量模型
τ=tB-tA (B.21)
式中:
tA———阶跃起点A对应的时间,s;
tB———响应曲线上B点对应的时间,s。
其中:
tB=t(TB)=t[xTC+(1-x)TA] (B.22)
式中:
TA———阶跃起点A对应的温度,℃;
TB———响应曲线上B点对应的温度,℃;
TC———阶跃终点C对应的温度,℃;
x ———温度变化到温度阶跃量的百分比。
则有:
uc(τ)= u2p+u2(tB)+u2(tA)+ ∂τ
∂Mau(Ma) ? ? ??
? ? ??
2 (B.23)
式中:
up ———测量重复性引入的不确定度;
∂τ
∂Ma ———时间常数对气流马赫数的偏导数;
Ma ———气流马赫数。
其中,u(Ma)按下式计算:
u(Ma)=2+(κ-1)Ma2
2κMa
u(p0)
p0
? ? ??
? ? ??
2+ u(ps)
ps
? ? ??
? ? ??
2 (B.24)
式中:
p0———气流总压;
ps———气流静压;
κ ———气体绝热指数,对于空气,通常取1.4,对于燃气,通常取1.33。
TA、TC由同一支传感器测得,因而可按完全正相关考虑,则有:
u(TB)=u[xTC+(1-x)TA]
=xu(TA)+(1-x)u(TA)=u(TA) (B.25)
因此,
u(tB)=∂tB
∂TB
u(TB)=∂tB
∂TB
u(TA) (B.26)
B.2.2 测量不确定度的来源
a)测量重复性引入的不确定度up。
b)u(TA)的来源如下:
1)温度传感器示值误差引入的不确定度;
17
JJF1049—2024
2)阶跃起点A选取引入的不确定度;
3)高速采集设备的示值误差引入的不确定度;
4)风洞气流的有效温度波动引入的不确定度。
c)u(tA)的来源如下:
1)数据采样间隔引入的不确定度;
2)阶跃起点A选取引入的不确定度。
d)u(p0)的来源如下:
1)总压传感器的示值误差引入的不确定度;
2)无汞气压计的示值误差引入的不确定度;
3)风洞气流的总压波动引入的不确定度;
4)数采系统引入的不确定度。
e)u(ps)的来源如下:
1)静压传感器的示值误差引入的不确定度;
2)无汞气压计的示值误差引入的不确定度;
3)风洞气流的静压波动引入的不确定度;
4)数采系统引入的不确定度。
B.2.3 测量不确定度的评定
以某热电偶式温度传感器在气流马赫数为0.4、气流温度从400℃阶跃到600℃的
时间常数校准为例,进行不确定度评定。
B.2.3.1 up的评定如下:
用极差法计算实验标准偏差s(τ)。在给定工况条件下共测量3次,3次时间常数τ
的测量结果分别为2.582s、2.569s、2.588s。从测量数据中找出最大值2.588s和最
小值2.569s,根据测量次数查表得到系数dn=1.693,按下式计算实验标准偏差s(τ):
s(τ)=τmax-τmin
dn =0.011s (B.27)
采用3次测量结果的平均值2.580s作为测量结果,测量重复性引入的不确定度up
按下式计算:
up=s(τ)
3 ≈0.006s (B.28)
B.2.3.2 u(TA)的评定如下:
a)温度传感器示值误差引入的不确定度
根据温度传感器的准确度等级,其最大允许误差为±0.4%,并经校准合格。由此
引入的不确定度按B类评定,假设测量值在允许误差范围内的概率分布为均匀分布,
查得包含因子k 为3,则其标准不确定度为:
u9=0.4%×400℃
3 ≈0.9℃ (B.29)
b)阶跃起点A选取引入的不确定度
根据经验,由阶跃起点选取不准带来的时间常数误差为±0.2℃。由此引入的不确
18
JJF1049—2024
定度按B类评定,假设测量值在允许误差范围内的概率分布为均匀分布,查得包含因
子k 为3,则其标准不确定度为:
u10=0.2℃
3 ≈0.12℃ (B.30)
c)高速采集设备的示值误差引入的不确定度
根据高速采集设备的准确度等级,其相对允许误差极限为±0.5%,并经校准合格。
由此引入的不确定度按B类评定,假设测量值在允许误差范围内的概率分布为均匀分
布,查得包含因子k 为3,则其标准不确定度u11为:
u11=0.5%×400℃
3 ≈1.2℃ (B.31)
d)风洞气流的有效温度波动引入的不确定度
按B类评定,估计风洞气流的总温波动值不超过±0.3℃。假设测量值在允许误差
范围内的概率分布为均匀分布,包含因子k 为3,则其标准不确定度u12为:
u12=0.3℃
3 ≈0.17℃ (B.32)
则有:
u(TA)= u2
9+u2
10+u2
11+u2
12 =1.5℃ (B.33)
B.2.3.3 u(tA)的评定如下:
a)数据采样间隔引入的不确定度
假定测量时设定的数据采样间隔为2ms。由此引入的不确定度按B类评定,假设
测量值在允许误差范围内的概率分布为均匀分布,查得包含因子k 为3,则其标准不
确定度u13为:
u13=0.001s
3 ≈0.0006s (B.34)
b)阶跃起点A选取引入的不确定度
根据经验,由阶跃起点选取不准带来的时间常数误差为±0.1s。由此引入的不确
定度按B类评定,假设测量值在允许误差范围内的概率分布为均匀分布,查得包含因
子k 为3,则其标准不确定度为:
u14=0.1s
3 ≈0.06s (B.35)
则有:
u(tA)= u2
13+u2
14 =0.06s (B.36)
B.2.3.4 u(p0)的评定如下:
a)总压传感器的示值误差引入的不确定度
根据压力传感器的准确度等级,其允许误差极限为±150Pa,并经校准合格。由此
引入的不确定度按B类评定,假设测量值在允许误差范围内的概率分布为均匀分布,
19
JJF1049—2024
查得包含因子k 为3,则其标准不确定度为:
u15=150Pa
3 ≈87Pa (B.37)
b)无汞气压计的示值误差引入的不确定度
根据无汞气压计的技术说明书规定,其允许误差极限为±40Pa,并经校准合格。
由此引入的不确定度按B类评定,假设测量值在允许误差范围内的概率分布为均匀分
布,查得包含因子k 为3,则其标准不确定度为:
u16=40Pa
3 ≈23Pa (B.38)
c)风洞气流的总压波动引入的不确定度
按B类评定,根据风洞校准结果,风洞气流的总压波动值不超过±10Pa。假设测
量值在允许误差范围内的概率分布为均匀分布,包含因子k 为3,则其标准不确定
度为:
u17=10Pa
3 ≈5.8Pa (B.39)
d)数采系统引入的不确定度
按B类评定,根据数采系统的技术说明书规定,其允许误差极限为±10Pa。按均
匀分布考虑,包含因子k 为3,则其标准不确定度为:
u18=10Pa
3 ≈5.8Pa (B.40)
e)计算u(p0)
u(p0)= u215+u216+u217+u218 =90Pa (B.41)
B.2.3.5 u(ps)的评定如下:
a)静压传感器的示值误差引入的不确定度
根据压力传感器的准确度等级,其允许误差极限为±6Pa,并经校准合格。由此引
入的不确定度按B类评定,假设测量值在允许误差范围内的概率分布为均匀分布,查
得包含因子k 为3,则其标准不确定度为:
u19=6Pa
3 ≈3.5Pa (B.42)
b)无汞气压计的示值误差引入的不确定度
同B.2.3.4b),u20=23Pa。
c)风洞气流的静压波动引入的不确定度
按B类评定,根据风洞校准结果,风洞气流的静压波动值不超过±10Pa。假设测
量值在允许误差范围内的概率分布为均匀分布,包含因子k 为3,则其标准不确定
度为:
u21=10Pa
3 ≈5.8Pa (B.43)
20
JJF1049—2024
d)数采系统引入的不确定度
同B.2.3.4d),u22=5.8Pa。
e)计算u(ps)
u(ps)= u219+u220+u221+u222 =25Pa (B.44)
f)计算u(Ma)
根据式(B.24)计算得到:
u(Ma)=0.0016 (B.45)
B.2.4 相对合成标准不确定度ur(τ)
表B.2 不确定度分量一览表
不确定
度分量
不确定度来源评定方法分布k 值
标准不
确定度
up 测量的重复性A ——— ——— 0.006s
u9 温度传感器的示值误差B 均匀3 0.9℃
u10 阶跃起点A 选取误差B 均匀3 0.12℃
u11 高速采集设备的示值误差B 均匀3 1.2℃
u12 风洞气流的有效温度波动B 均匀3 0.17℃
u13 数据采样间隔B 均匀3 0.0006s
u14 阶跃起点A 选取误差B 均匀3 0.06s
u15 总压传感器的示值误差B 均匀3 87Pa
u16 无汞气压计的示值误差B 均匀3 23Pa
u17 风洞气流的总压波动B 均匀3 5.8Pa
u18 数采系统的测量误差B 均匀3 5.8Pa
u19 静压传感器的示值误差B 均匀3 3.5Pa
u20 无汞气压计的示值误差B 均匀3 23Pa
u21 风洞气流的静压波动B 均匀3 5.8Pa
u22 数采系统的测量误差B 均匀3 5.8Pa
不确定度分量一览表见表B.2。其中各不确定度分量互不相关,则时间常数的相对
合成标准不确定度ur(τ)按下式计算:
ur(τ)=
u2p +u2(tB)+u2(tA)+ ∂τ
∂Mau(Ma) ? ? ??
? ? ??
2
τ
=
u2p + ∂tB
∂TB
u(TA) ? ? ??
? ? ??
2
+u2(tA)+ ∂τ
∂Mau(Ma) ? ? ??
? ? ??
2
τ
=3.7% (B.46)
21
JJF1049—2024
其中,∂tB
∂TB
具体计算方法为:以温度为横坐标,时间为纵坐标,通过对TA、TB、
TC对应的时间tA、tB、tC进行多项式拟合,对拟合的多项式求导,求得B点处的切线
斜率,计算得到∂tB
∂TB=0.5。∂τ
∂Ma具体计算方法为:以马赫数为横坐标,时间常数为纵
坐标,通过对不同马赫数下对应时间常数点进行多项式拟合,对拟合的多项式求导,求
得对应马赫数的切线斜率,计算得到∂τ
∂Ma=0.5。
B.2.5 相对扩展不确定度Ur(τ)
取k=2,则相对扩展不确定度Ur(τ)为:
Ur(τ)=k·ur(τ)=7.4% (B.47)
B.3 热响应时间的测量不确定度评定
热响应时间的测量不确定度评定方法与时间常数相同,可参照B.1~B.2进行。
22
JJF1049—2024
附录C
校准结果原始记录
委托方名称和地址
被校件描述
型号编号生产厂家
接收日期年 月 日校准日期年 月 日
校准用标准器描述
器具名称器具编号测量范围
准确度等级证书编号标准有效期
介质类型
校准依据文件
校准环境条件温度: ℃ 相对湿度: %
校准人员核验人员
编号次数
马赫数Ma
或速度
v/(m/s)
阶跃前
温度
T1/℃
阶跃后
温度
T2/℃
时间常数
τ/s
热响应
时间
τ0.1/s
热响应
时间
τ0.5/s
热响应
时间
τ0.9/s
123
平均
23
JJF1049—2024
附录D
校准证书格式
证书编号:
校准结果及数据
传感器型号传感器编号
马赫数Ma 或
速度v/(m/s)
阶跃前温度
T1/℃
阶跃后温度
T2/℃
温度阶跃量
ΔT/℃
时间常数
τ/s
相对扩展
不确定度
Ur(τ)
(k=2)
热响应
时间
τ0.1/s
相对扩展
不确定度
Ur(τ0.1)
(k=2)
热响应
时间
τ0.5/s
相对扩展
不确定度
Ur(τ0.5)
(k=2)
热响应
时间
τ0.9/ms
相对扩展
不确定度
Ur(τ0.9)
(k=2)
24
JJF1049—2024