基于热电偶试验数据拟合的温度检测仪设计
在航空发动机排气温度检测、热处理、铸造等广泛领域都需要用到温度仪表来检测温度，在这些仪表中用得最多的就是热电偶仪表。然而由于热电偶的非线性特点，使得热电偶的测温电路非常复杂。例如为了实现非线性变化量的测量，传统测温仪表电路往往需要十余个集成运放和近百个外围电阻、电容等元件，不仅电路极为复杂，而且易受环境温度影响，测量精度也难以提高。另一方面，由于用以计算温度的热电偶标准热电势是在其冷端温度为0 ℃时测量的[3]，而实际应用中，热电偶冷端温度往往不是0 ℃，受环境影响在-50 ℃～+50 ℃范围内变化，这种大幅度的环境温度变化使热电偶的测量产生巨大误差，因此需对热电偶进行冷端补偿。同样，传统仪表的冷端补偿电路也存在电路复杂、调试困难等缺点。
针对上述传统测温电路的缺点，可对简单模拟电路的输入、输出数据进行拟合，将拟合结果交给数据处理功能强大的单片机处理，就可以使电路的设计结构大大简化，并获得更高的检测精度。本文在对大量试验结果的数据拟合和电路仿真的基础上，提出了新的热电偶测温仪表设计方法，以期为该类型测温仪表的生产提供参考。
1 热电偶的冷端温度补偿
在电极材料和冷端温度保持不变的情况下，热电偶的热电势为热端温度(热源温度)的单值函数。通常工作温度是由冷端为0 ℃时的标准热电势来反映，因此需要对非0 ℃的冷端进行温度补偿。其依据是中间温度定律[1]：
已知热电偶的热端温度为T(℃)，冷端温度为0 ℃时的热电势为E(标准电势)。如果：(1)该热端温度为Tn(环境温度)，冷端温度为0 ℃时的热电势为E1(环境电势); (2)该热电偶的热端温度为T，冷端温度为Tn时的热电势为E2(现场电势)。那么：
E=E1+E2 (1)
根据该定律，只要设法测出环境电势E1，并将其与热电偶测出的现场电势E2求和，即可计算出对应热源温度的标准热电势E，并由已知试验结果确定出热端温度T。
2 环境电势测量电路设计
热电偶的环境电势是环境温度的单值函数，所以只要利用热电阻等温度传感器测出环境温度就可以由已知的试验结果确定出环境电势。对K型热电偶，其环境电势随环境温度的变化关系可由试验数据的最小二乘拟合得到。其计算公式如下：
E1=0.0392×Tn+0.0340 (2)
式(2)中：E1为热电偶的环境电势(mV); Tn为环境温度(℃)。式(2)适用的温度范围为-50 ℃~+50 ℃。
式(2)中的环境温度Tn可由图1所示的电路实现。其中R22、R16、R17和R18组成测温电桥，其输出信号接差动集成运放UA741。R18为WZB型铂热电阻(0℃时，标称电阻值为100Ω)，R16为68Ω标准电阻。该电路可以获得近似恒流法的线性输出电压Vi2和Vi3。显然Vi3与热电阻R18存在线性关系，所以热电阻随温度的变化可以线性地反映到输出端Vo1。通过试验测出不同温度下的输出电压Vo1，将试验结果拟合可确定出环境温度Tn与输出电压Vo1的关系如下：
Tn=33.6080×Vo1-76.8451 (3)
式(3)中，Tn为环境温度(℃); Vo1为输出电压(V)。将模拟量Vo1经A/D转换芯片转换为数字量后送入单片机即可计算出环境温度Tn，再将Tn代入式(2)即可由单片机计算出环境电势E1。
3 现场电势测量电路设计
图2为测量现场电势E2的电路，热电偶所测出的现场电势由电阻R14输入到集成运算放大器UA741。由于集成运放并非完全理想，在精确的测量中，其输出电压Vo2与现场电势E2间的计算关系不能简单地采用理想集成运放的公式确定，而应依据试验建立。本文由数据拟合而确定的现场电势计算公式如下：
E2=16.3017×Vo2-1.2603 (4)
式中，E2为现场电势(mV); Vo2为检测电路的输出电压(V)。将Vo2经A/D转换后的数字量送入单片机，即可由式(4)计算出现场电势E2。
4 热电势测温信号处理电路设计
在确定了环境电势E1和现场电势E2后，即可由式(1)确定热电偶的标准电势E。根据试验结果[4]，可得出K型热电偶热端温度T和标准电势E存在如下拟合关系：
T＝24.5020×E-4.0446 (5)
式中，T为待检测的热端温度(℃)，E为热电偶冷端温度为0℃时的标准电势(mV)。图3为由式(5)拟合的曲线，可以看出，该拟合公式具有很高的计算精度。依据式(5)可由单片机程序计算出热端温度。
图4为热端温度信号处理电路图(及仿真效果)，主要由A/D转换器件ADC0832、单片机AT89S51、液晶显示器LCD1602、上拉电阻R2～R9等组成。
ADC0832主要完成模拟电压信号的采集与转换，该芯片有2个输入通道CH0和CH1，其中CH0输入现场电势信号，接图2的输出电压信号Vo2；CH1输入环境电势信号，接图2的输出电压信号Vo1，2个电势信号的采集次序受单片机AT89S51的指令控制。该芯片转换后的数字信号VAL与输入模拟电压Vi间的计算关系如下[5]:
VAL=51×Vi (6)
式中，VAL为输出的数字信号(0～255); Vi为输入的模拟电压(V)。当数字信号VAL输入单片机后，即可依据(6)式求出输入电压Vi，并进而由(2)～(3)式计算出环境电势E1和现场电势E2。将两者通过单片机求和即可算出标准电势E。