直读光谱法测定低合金钢中碳的测量不确定度评定报告

1、直读光谱法测定低合金钢中碳的测量不确定度一、 概 述1、目的用直读光谱法测定低合金钢中碳的测量不确定度。2、检测依据的标准GB/T4336-2002碳素钢和中低合金钢火花源原子发射光谱分析方法常规法3、检测使用的仪器设备和检测环境1直读光谱仪： Q8 Magellan德国利曼型多通道真空直读光谱仪。光栅刻线数高2400线/毫米、光栅焦距750毫米、波长范围110 nm 800 nm、真空度不高于0.03mmHg、色散一级谱线0.52 nm/mm 、二级谱线0.26 nm/mm、三级谱线0.13 nm/mm。工作环境: 相对湿度45%,温度23,t±2。2标准物质和试样GBW01328

2、-01333光谱标准物质(32×28),低合金钢试样(30mm×30mm×30mm),外表用60目的砂轮抛光。3校准曲线拟合每块标准物质激发3点(RSD3%),根据元素含量与元素的相对光强进行曲线拟合。4试验程序按GB/T4336-2002标准进行测定,分析线为C193.0nm,内标线为Fe 271.4nm。分析试样激发两次(两次结果超过本实验室内允许差时,须重新分析),结果取平均值。二、 不确定度来源 电火花直读光谱法测定低合金钢中碳的不确定度主要来自试样的不均匀性、电火花光源的不稳定性、光电倍增管负高压的不稳定性、光电流与光强的转换过程、温湿度变化引起的分光系

3、统的漂移、校准曲线的线性拟合、标准物质的定值。其中前4个因素对测定结果产生的不确定度集中表达在校准曲线回归拟合引入的不确定度；由于每次的实验环境是严格控制的,可以认为仪器标准化后的一定时间内温度和湿度几乎不发生变化,本文把温度和湿度变化引起的分光系统的漂移对测量结果所产生的不确定度予以忽略；因而电火花直读光谱法测定低合金钢中碳的不确定度主要来自校准曲线的回归拟合和标准物质的定值。三、 标准不确定度评定1、校准曲线非线性回归拟合对测定结果y0引入的标准不确定度1一元线性拟合对测定结果y0引入的标准不确定度一元线性回归的理论方程可以表示为Y=a+bx,a、b可由一组观测值确定。i=(yi-Yi)=

4、yi-(a+bxi)称为变差(又称残差),yi是实测值,Yi是按回归方程计算的预测值。变差平方和可用来表征实测值yi偏离回归方程预测值Yi的程度。一组实测值总的偏离程度可用变差平方和来定量描述。当回归方程拟合最正确时,Q必然最小。满足最小的条件是： 解方程可求得： 式中 在实际当中各实验点不一定都落在回归线上,围绕回归线有一定的离散性,其离散性的大小可用残差标准差S来表征。假设对y0进行p次重复测量,一元线性拟合对测量结果y0引入的标准不确定度为:2非线性回归拟合对测定结果y0引入的标准不确定度对于y与x之间呈非线性关系的情况下,需通过变量变换,将非线性回归变为线性回归,然后按照线性拟合的方法

5、处理数据。此时S的计算公式为：中yi是实测值,i是按回归方程计算的预测值。3二次校准曲线的回归拟合对C0(钢中碳的含量)引入的标准不确定度的评定a校准曲线拟合 采用电火花光源进行元素分析时,校准曲线一般是二次曲线,其方程为:y=a+bx+cx2,式中y是元素含量(%),x是相对光强IR(IR=105×I被测元素/I内标元素),a、b、c可通过一组观测值确定。校准曲线拟合所需的数据列于表1。计算a、b、c所得的数据列于表2。表1 校准曲线拟合数据y(%)0.03320.1880.2830.3920.5060.569备注x184992376431971403514717952742x28

6、2352326231778403994684052241x383172344732151401084691251707x()83502349131967402864697752230表2 计算a、b、c数据y(%)xx2x3x4xyx2y10.03328350697225005.82183E+114.86123E+15277.22231478720.188234915518270811.2963E+133.04513E+174416.3081037434930.2833196710218890893.26667E+131.04426E+189046.6612891946140.39240286

7、16229617966.53826E+132.634E+1815792.11263620212450.5064697722068385291.03671E+144.87014E+1823770.36211166602960.5695223027279729001.42482E+147.44184E+1829718.8715522165801.971220330182021118953.57747E+141.62996E+1983021.5333700330790根据最小二乘法进行拟合。a、b、c分别为：代入数据可得a=-0.03462；b=7.473E-6；c=8.022E-11。那么拟合的二

8、次校准曲线为：y(%)=-0.0346+7.473×10-6x+8.022×10-11x2b.碳的测定激发低合金钢Q345E试样两次,IR分别为28389、28486,代入校准曲线均值为C0=0.243%。c.校准曲线回归拟合对C0引入的标准不确定度进行变量置换,令那么二次校准曲线置换为：y(%)=cz+(a-b2/4c),代入a、b、c值,可得：y(%)=-0.2087+8.022×10-11z那么S、u(y0)的公式为： 计算所得的数据列于表3：y(%)0.03320.1880.2830.3920.5060.569y1(%)0.034700.18830.286

9、30.39760.49650.5827y2(%)0.032380.18260.28390.39820.49140.5747y3(%)0.033100.18470.28860.39420.49250.5663z13.0335E+94.9481E+96.1700E+97.5567E+98.7904E+99.8645E+9z23.0045E+94.8776E+96.1397E+97.5650E+98.7270E+99.7652E+9z33.0135E+94.9097E+96.1983E+97.5145E+98.7404E+99.6600E+9表3中 ，代入C0 z0=5.6372E+9 那么S=0.

10、00739% uC0=0.00557% =163、标准物质的定值对C0引入的标准不确定度的评定标准物质说明书提供了单次标准偏差S和n=8,通常以算术平均值的标准偏差作为标准物质的标准不确定度。标准物质的标准偏差列于表4。表4 标准物质的标准偏差标样GBW 01328GBW 01329GBW 01330GBW 01331GBW 01332GBW 01333y(%)0.03320.1880.2830.3920.5060.569S(%)0.00230.0060.0090.0050.0060.005(%)0.000810.002120.003180.001770.002120.00177选择最大的S作为标准物质的定值对C0引入的标准不确定度,那么u(RM)=0.0032% =64、合成标准不确定度的计算1合成标准不确定度的计算Ur(C0) =(u(C0)2+u(RM)2)1/2那么uc(C0)=0.0064%。2 有效自由度的计算 那么可近似认为：veff。2 扩展不确定度的评定a. t分布法