(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111552692.7 (22)申请日 2021.12.17 (71)申请人 北京航空航天大 学 地址 100191 北京市海淀区学院路37号 (72)发明人 孙富强　郭宏轩　 (74)专利代理 机构 北京天汇航智知识产权代理 事务所(普通 合伙) 11987 代理人 黄川　史继颖 (51)Int.Cl. G06F 17/12(2006.01) G06F 17/13(2006.01) G06F 17/18(2006.01) G06F 30/20(2020.01) (54)发明名称 一种加速退化试验测量不确定度的评定方 法 (57)摘要 本发明涉及一种加速退化试验不确定性评 价方法， 基于现有加速退化模型， 以及模型参数 的估计过程， 得到加速退化模型参数的测量不确 定度、 可靠度的测量不确定度， 用于综合评价加 速退化试验的整体测量不确定度， 包括： 步骤一： 加速退化试验性能参数的测量不确定的评价， 包 括重复测量过程中的A类不确定度以及机器精度 影响的B类不确定度； 步骤二： 加速退化模型参数 的测量不确定度评定； 步骤三： 加速退化试验可 靠度的测量不确定度评定； 基于所述加速退化模 型参数的估计结果的测量不确定度， 以及正常应 力下试验的可靠度R(t)估计结果； 得到加速退化 试验可靠度估计结果的测量不确定度uc(R(t))。 本发明能够对加速退化试验的测量不确定度进 行综合评价， 更精准。 权利要求书2页 说明书12页 附图3页 CN 114266006 A 2022.04.01 CN 114266006 A 1.一种加速退化试验测量不确定度的评定方法， 其特征在于， 基于已知的加速退化模 型， 以及模型参数的估计过程， 结合测量不确定度的评价方法， 计算加速退化试验性能参数 的测量量的测 量不确定度、 加速退化模型参数 的测量不确定度、 加速退化试验可靠度的测 量不确定度， 实现对加速退化试验的测量 不确定度的综合评价， 具体步骤如下： 步骤一： 所述加速退化试验直接测量的性能参数的测量 不确定度的评价: (1)A类不确定度uA的计算： 性能参数测量量Y， 其满足重复性测量条件的n次重复测量中 的第i次测量值记为yi； ， n为自然数， i＝1, …,n， 重复测量所述测量量Y得到的性能值可被 用于A类评 定， 计算方法有 贝塞尔公式法、 极差法、 最大误差法； (2)B类不确定度uB的计算： B类不确定度指的是用不同于所述A类不确定度计算方法的 测量不确定度评定； 其来源有权威机构发布的量值、 有证标准物质的量值、 校准证书、 测量 仪器的测量精度、 经检定的测量仪器的准确度等级； (3)得到所述 性能参数测量 量Y的测量确定度u(yi)为： 在所述加速退化试验中， 所述性能参数测量量为已知信息， 根据上述公式可得到所述 性能参数测量 量的测量 不确定度； 步骤二： 所述加速退化模型参数的测量 不确定度评 定； 假设所述加速退化模型参数为θ＝{A,Ea,C, σ }， 其中模型参数A、 Ea、 C为常数， σ 是扩散系 数， 各参数均在建模过程中通过已知的模型参数估计过程计算得到； 基于对所述模型参数 A、 Ea、 C、 σ 的估计过程 得到各所述模型参数的测量不确定度u(A)、 u(Ea)、 u(C)以及u( σ )， 再根 据正常应力水平下的加速退化模型的退化率 μ0与所述模型参数之间的关系， 计算出正常应 力水平下的加速退化模型 退化率 μ0的测量不确定度u( μ0)； 步骤三： 所述加速退化试验可靠度的测量 不确定度评 定； 基于所述加速退化模型参数的估计过程， 以及正常应力水平下的所述加速退化试验的 可靠度R(t)估计过程； 得到所述加速退化试验的所述可靠度的测量不确定度uc(R(t))， 其 中， R(t)为所述加速退化试验在正常应力水平下的可靠度， 是时间t的函数。 2.根据权利要求1所述的加速退化试验测量不确定度的评定方法， 其特征在于， 在步骤 一中， 用贝塞尔公 式来计算测量数据的所述A类不确定度， 得到所述测量量Y的所述A类不确 定度如下： 其中， uA(yi)是单次测量结果的标准不确定度， 是n次测量平均值的标准不确定 度， 为n次重复测量 值的均值。 3.根据权利要求1所述的加速退化试验测量不确定度的评定方法， 其特征在于， 在步骤 一中， B类不确定度uB由所述测量仪器已知的所述测量精度计算得到， 具体通过： 得到被测 量值的可能值 区间半宽度aB， 再选取合适的因子kB， 则最终计算测定测量量的B类不确定度权　利　要　求　书 1/2 页 2 CN 114266006 A 2为： 4.根据权利要求1 ‑3中任一项所述的加速退化试验测量不确定度的评定方法， 其特征 在于， 在步骤二 或者步骤三中， 如果一个目标量的估计过程是通过多个相关量计算得到， 则 需要基于所述多个相关量的测量不确定度 计算所述目标量的合成测量不确定度uc； 假设所 述目标量f与n个相关量xi的关系式为f＝f(x1,x2,…,xn)； 则所述目标量f的合成测量不确 定度uc的一般性公式为： 其中， ρij表示不同相关量的测量不确定度u(xi)与u(xj)之间的相关系数， 通常为0 ‑1之 间的数值， 表示测量不确定度u(xi)的传递系数， 其中， n为大于等于2的自然数， i,j＝ 1,…,n， j为不同于i的自然数。 5.根据权利要求1所述的加速退化试验测量不确定度的评定方法， 其特征在于， 在步骤 二中， 所述加速退化模型为Peck模型时， 可以采用二元线性回归的不确定度评定方法对所 述各参数 的估计过程进行计算， 并结合所述合成测量不确定度的计算公式， 最终得到所述 参数A、 Ea、 C的所述测量 不确定度u(A)、 u(Ea)、 u(C)的数值。 6.根据权利要求4所述的加速退化试验测量不确定度的评定方法， 其特征在于， 在步骤 二中， 所述加速退化模型选择为Peck模型， 此时， 所述Peck模型的退化率为μ0＝exp(A‑Ea/ k0T0+C ln(H0))， Ea为激活能， H0为正常湿度水平的相对湿度值， 单位为％RH， T0为正常温度 水平的开氏温度， 单位为K， k0为波耳兹曼常数， 假设所述加速退化模型参数之间为强相关， 每两参数间的相关系数均为1， 由所述加速退化模 型退化率的公 式(11)， 以及所述合成测量 不确定度的计算公式(5)， 可 得正常应力下的加速退化模型 退化率 μ0的不确定度为： 7.根据权利要求5所述的加速退化试验测量不确定度的评定方法， 其特征在于， 在步骤 三中， 所述可靠度(R(t)是通过与所述退化率μ0以及的σ 是扩散系数估计过程计算得到， 所 述可靠度的不确定度为 合成的测量 不确定度， 具体为： 8.根据权利要求1所述的加速退化试验测量不确定度的评定方法， 其特征在于， 所述正 常应力水平选择为在所述正常湿度水平H0为50％RH， 所述正常温度水平T0为298.15K的情 况。权　利　要　求　书 2/2 页 3 CN 114266006 A 3