CN108549047A - 一种基于加速退化试验的电能表寿命评估方法 - Google Patents
一种基于加速退化试验的电能表寿命评估方法 Download PDFInfo
- Publication number
- CN108549047A CN108549047A CN201810521253.1A CN201810521253A CN108549047A CN 108549047 A CN108549047 A CN 108549047A CN 201810521253 A CN201810521253 A CN 201810521253A CN 108549047 A CN108549047 A CN 108549047A
- Authority
- CN
- China
- Prior art keywords
- electric energy
- energy meter
- accelerated
- life
- degradation
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R35/00—Testing or calibrating of apparatus covered by the other groups of this subclass
- G01R35/04—Testing or calibrating of apparatus covered by the other groups of this subclass of instruments for measuring time integral of power or current
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
Abstract
本发明涉及一种基于加速退化试验的电能表寿命评估方法,包括:1、采集电能表加速退化试验数据,估计得出电能表退化特征量在不同应力水平下退化轨迹模型的分布参数值;2、根据加速应力的类型选择加速退化模型,根据试验中退化特征量的失效阈值计算出不同应力水平下电能表的伪寿命值;3、根据加速退化模型计算出不同应力水平对应的加速因子,并利用加速因子把不同应力水平下的伪寿命值折算到正常工作应力水平下;4、选定电能表的寿命分布类型,并估计出寿命分布函数中的未知参数值;5、利用寿命分布函数评估得到电能表在正常工作应力水平下的寿命。与现有技术相比,本发明能够评估得到较准确的电能表寿命情况,更符合工程实际。
Description
技术领域
本发明涉及计量设备运行检定检测技术领域,尤其是涉及一种基于加速退化试验的电能表寿命评估方法。
背景技术
通过对智能电能表进行可靠性分析与寿命评估,电能表生产制造厂家能够提前发现产品设计的薄弱环节或元器件的质量缺陷,并针对性地进行改进,从而整体提高智能电能表的质量水平,保障智能电能表在现场能够长期可靠、稳定运行。
采用传统的可靠性评价方法对智能电能表进行可靠性评价需要花费很大的时间成本和人力成本,因此目前广泛采用加速寿命试验的方法对智能电能表的可靠性进行评价。加速寿命试验只记录试验样品的失效时间数据,而不考虑产品失效的具体过程,因此需要有一定数量的样品发生失效才能结束试验。如果试验中没有样品发生失效,或者失效样品的数量很少,那么就只能延长试验时间,否则无法利用获得的失效时间数据来对产品的可靠性进行评价,或者得出的评价结果与实际情况有很大的差距。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于加速退化试验的电能表寿命评估方法。
本发明的目的可以通过以下技术方案来实现:
一种基于加速退化试验的电能表寿命评估方法,包括以下步骤:
S1、采集电能表加速退化试验数据,估计得出电能表退化特征量在不同应力水平下退化轨迹模型的分布参数值;
S2、根据加速应力的类型选择加速退化模型,根据试验中退化特征量的失效阈值计算出不同应力水平下电能表的伪寿命值;
S3、根据加速退化模型计算出不同应力水平对应的加速因子,并利用加速因子把不同应力水平下的伪寿命值折算到正常工作应力水平下;
S4、以所有折算后的伪寿命值作为整体,选定电能表的寿命分布类型,并估计出寿命分布函数中的未知参数值;
S5、利用寿命分布函数评估得到电能表在正常工作应力水平下的寿命。
优选的,所述退化轨迹模型的分布参数值根据最小二乘法估计得出。
优选的,所述步骤S2中加速退化模型选择Hallberg-Peck模型。
优选的,所述退化特征量包括日计时误差和基本误差。
优选的,所述步骤S4中选定电能表的寿命分布类型的方法包括:
根据Anderson-Darling拟合优度检验方法比较各种分布类型的拟合情况,选择Anderson-Darling统计量最小的分布类型为寿命分布类型。
优选的,所述寿命分布类型选为两参数威布尔分布。
与现有技术相比,本发明利用性能退化数据得到电能表在正常工作时的寿命分布规律,进而对智能电能表的寿命进行评估,本方法能够准确地反映现场在运电能表的失效情况,评估得到较准确的电能表寿命情况,更符合工程实际。
附图说明
图1为本发明的方法流程图;
图2为基本误差伪寿命数据分布拟合图;
图3为日计时伪寿命数据分布拟合图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
实施例
本实施例中,电能表为智能电能表,对其采用加速退化试验,过程包括:
步骤(1)、确定引起现场电能表故障的原因。
步骤(2)、根据步骤(1)确定的故障原因结合摸底试验结果选择加速退化试验的加速应力类型,本实施例中选择温度和湿度作为加速退化试验的加速应力。
步骤(3)、通过摸底试验确定最大应力水平并设计加速应力组合。
加速退化试验最大加速应力的确定原则是,保证产品在最大应力下的性能退化机理不发生改变,并且与产品在现场的失效机理或故障机理一致。为确定加速退化试验的最大应力,本实施例中对某一批2级单相智能电能表进行摸底试验,试验中当施加的温-湿度应力为[85℃,95%]时,智能电能表的性能退化机理未改变,因此可选择温-湿度应力水平[85℃,95%]作为最大应力水平,即Tmax为85℃,RHmax为95%。加速退化模型选择Hallberg-Peck模型,可以用统一的线性形式表示为:
其中,ξi为产品的寿命特征;为应力Si的函数;a、b为常数。由上式可知智能电能表的寿命特性近似服从线性分布模型。为提高模型的“线性度”,在加速退化试验中,对于每一种应力(温度和湿度)类型,至少要使用3种不同的应力水平。定义一个中间的和低的温度应力水平,分别记作Tmed和Tmin;定义一个中间的和低的相对湿度应力水平,分别记作RHmed和RHmin。本实施例中将选择5组应力组合,如表1所示:
表1温-湿度加速应力组合
序号 | 应力S | T/℃ | RH/% |
S1 | Tmin,RHmax | 65 | 95 |
S2 | Tmed,RHmax | 75 | 95 |
S3 | Tmax,RHmin | 85 | 75 |
S4 | Tmax,RHmed | 85 | 85 |
S5 | Tmax,RHmax | 85 | 95 |
假设产品的使用应力水平为So,产品在应力So、Si下的寿命特征值为ξo、ξi,则温-湿度模型的加速因子为:
表1中不同应力组合的加速因子如表2所示:
表2不同应力组合的加速因子
步骤(4)、根据需要确定退化特征量及其失效阈值、试验样品数量、分配比例、试验结束时间和测量次数。
计量计费功能是智能电能表最基本也是最重要的功能,因此本方法选取基本误差和日计时误差作为加速退化试验主要测试的退化特征量。两个特征量的测量均是在被试样品施加额定电压UN、基本电流Ib、功率因数的条件下测量。
退化特征量的失效阈值D与JJG 596-2012《电子式交流电能表》中规定的电能表基本误差和日计时误差的误差限值有关。根据经验,一般取误差限值的55%~65%作为退化特征量的失效阈值,本实施例中,优选的选取60%,即
若被试样表的基本误差和日计时误差超过相应的失效阈值,即认为该样品失效。
本实施例中,试验样品取同一厂家同一生产批次的150只2级单相智能电能表进行加速退化试验,5种加速应力组合下的样品分配比例相同,均为0.2,即每组加速应力组合下加速退化试验的样品数量均为30只。
采取定时截尾的方法来结束试验,最大试验时间tMAX取1080小时,即45天。测量时间点等间隔,即每组加速应力下每隔24小时对退化特征量进行一次测量,均为45次。
步骤(5)、对各组试验样品分别持续施加相应的应力组合,每间隔设定时间卸载应力组合并测试需要的参数,直到达到试验结束时间,具体包括:
①分别对各组被试样表施加表1中的温-湿度应力组合,并保持恒定;
②每隔24h卸载所施加的温-湿度应力,恢复至常温常湿条件;
③执行测试,测试项目包括:基本误差、日计时误差、电池电压、起动试验、潜动试验、基本功能试验、显示、RS485通信,记录其性能退化数据;
④重复过程①、②、③,直至试验时间达到最大试验时间tMAX。
基本误差、日计时误差、起动试验、潜动试验等项目测试方法应符合最新版JJG307、JJG 596及国网智能电能表系列企业标准的要求;电池电压通过读取DL/T645-2007数据标识02800008H(时钟电池电压)和02800009H(停电抄表电池电压)获得;基本功能试验指电能表费控功能,包括报警和拉合闸操作;RS485通信测试指通过RS485方式抄读电表当前电量数据、瞬时电压电流功率数据、日冻结数据,若均抄读成功,则判定RS485通信正常,否则记录为异常。
如图1所示,一种基于上述加速退化试验的电能表寿命评估方法,包括以下步骤:
S1、采集上述电能表加速退化试验数据,估计得出电能表退化特征量在不同应力水平下退化轨迹模型的分布参数值。
退化轨迹模型的分布参数值根据最小二乘法估计得出。选用线性模型x=ωt+λ作为电能表基本误差和日计时误差的退化轨迹模型,ω和λ表示未知的分布参数,t表示时间,使用最小二乘法对加速退化试验得到的退化数据进行拟合,得到电能表基本误差和日计时误差在不同应力水平下退化轨迹。其中一个被试样品的退化轨迹模型为:
在应力水平S1下:
x1=0.0533t+0.0002
x2=0.0134t+0.0001;
在应力水平S2下:
x1=0.0845t+0.0001
x2=0.0211t+0.0001;
在应力水平S3下:
x1=0.1230t+0.0001
x2=0.0307t+0.0002;
在应力水平S4下:
x1=0.1788t+0.0002
x2=0.0447t+0.0001;
在应力水平S5下:
x1=0.2499t+0.0003
x2=0.0625t+0.0001。
S2、已选择加速退化模型为Hallberg-Peck模型,根据试验中退化特征量的失效阈值计算出不同应力水平下电能表的伪寿命值。
把电能表的基本误差和日计时误差的试验阈值代入不同应力下的退化轨迹模型,可得到该被试样品的伪寿命值t*,t1 *和t2 *分别表示不同应力下基本误差和日计时误差的伪寿命值,如表3所示。
表3不同应力水平下的伪寿命值
应力水平 | 伪寿命值t1 */天 | 伪寿命值t2 */天 |
S1 | 22.5103 | 21.3806 |
S2 | 14.2000 | 14.2133 |
S3 | 9.5553 | 9.7655 |
S4 | 6.7103 | 6.7092 |
S5 | 4.8007 | 4.7984 |
S3、根据加速退化模型计算出不同应力水平对应的加速因子,并利用加速因子把不同应力水平下的伪寿命值折算到正常工作应力水平下。
把不同应力水平下的伪寿命值t*折算到正常工作应力水平下得到电能表在正常工作应力水平下的伪寿命值,t1 *和t2 *分别表示正常应力水平下基本误差和日计时误差的伪寿命值,如表4所示。
表4正常工作应力水平下的伪寿命值
序号 | 伪寿命值t1 **/天 | 伪寿命值t2 **/天 |
1 | 3646.8937 | 3463.9710 |
2 | 3656.9260 | 3660.3511 |
3 | 3573.4911 | 3652.1017 |
4 | 3657.2477 | 3656.6482 |
5 | 3657.3653 | 3655.6131 |
S4、以所有折算后的伪寿命值作为整体,选定电能表的寿命分布类型,并估计出寿命分布函数中的未知参数值。
得到被试样品在正常工作应力水平下的伪寿命折算值后,可通过Anderson-Darling拟合优度检验方法比较各种分布类型的拟合情况,根据假设检验的P值确定样品的寿命分布类型并检验该寿命分布类型。
Anderson-Darling统计量具有良好的统计特性,不但可用于最优拟合检验,比较若干分布的拟合情况,以查看哪种分布是最佳分布类型,而且可以对特定分布类型进行拟合优度的假设检验。Anderson-Darling统计量可描述数据服从特定分布类型的程度,数据与分布拟合得越好,统计量的值越小。
Anderson-Darling的假设检验如下:
H0:寿命数据服从指定分布类型;
H1:寿命数据不服从指定分布类型。
假设显著性水平为0.05,如果Anderson-Darling检验的P值大于0.05,则可以得出寿命数据服从指定分布类型的结论。
经最优拟合检验,智能电能表的寿命分布模型采用两参数威布尔分布,根据现有方法可以得到其未知参数值,例如最大似然法。
S5、利用寿命分布函数评估得到电能表在正常工作应力水平下的寿命。
对电能表在正常工作应力水平下的伪寿命值进行威布尔分布拟合,基本误差和日计时误差的威布尔分布拟合结果分别如图2和图3所示。
Claims (6)
1.一种基于加速退化试验的电能表寿命评估方法,其特征在于,包括以下步骤:
S1、采集电能表加速退化试验数据,估计得出电能表退化特征量在不同应力水平下退化轨迹模型的分布参数值;
S2、根据加速应力的类型选择加速退化模型,根据试验中退化特征量的失效阈值计算出不同应力水平下电能表的伪寿命值;
S3、根据加速退化模型计算出不同应力水平对应的加速因子,并利用加速因子把不同应力水平下的伪寿命值折算到正常工作应力水平下;
S4、以所有折算后的伪寿命值作为整体,选定电能表的寿命分布类型,并估计出寿命分布函数中的未知参数值;
S5、利用寿命分布函数评估得到电能表在正常工作应力水平下的寿命。
2.根据权利要求1所述的一种基于加速退化试验的电能表寿命评估方法,其特征在于,所述退化轨迹模型的分布参数值根据最小二乘法估计得出。
3.根据权利要求1所述的一种基于加速退化试验的电能表寿命评估方法,其特征在于,所述步骤S2中加速退化模型选择Hallberg-Peck模型。
4.根据权利要求1所述的一种基于加速退化试验的电能表寿命评估方法,其特征在于,所述退化特征量包括日计时误差和基本误差。
5.根据权利要求1所述的一种基于加速退化试验的电能表寿命评估方法,其特征在于,所述步骤S4中选定电能表的寿命分布类型的方法包括:
根据Anderson-Darling拟合优度检验方法比较各种分布类型的拟合情况,选择Anderson-Darling统计量最小的分布类型为寿命分布类型。
6.根据权利要求1所述的一种基于加速退化试验的电能表寿命评估方法,其特征在于,所述寿命分布类型选为两参数威布尔分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810521253.1A CN108549047A (zh) | 2018-05-28 | 2018-05-28 | 一种基于加速退化试验的电能表寿命评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810521253.1A CN108549047A (zh) | 2018-05-28 | 2018-05-28 | 一种基于加速退化试验的电能表寿命评估方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108549047A true CN108549047A (zh) | 2018-09-18 |
Family
ID=63511301
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810521253.1A Pending CN108549047A (zh) | 2018-05-28 | 2018-05-28 | 一种基于加速退化试验的电能表寿命评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108549047A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110068789A (zh) * | 2019-05-23 | 2019-07-30 | 华立科技股份有限公司 | 电能表寿命的监测方法和电路 |
CN111025224A (zh) * | 2020-01-07 | 2020-04-17 | 中国电力科学研究院有限公司 | 一种智能电表的加速退化试验方法 |
CN111351697A (zh) * | 2020-03-24 | 2020-06-30 | 中国电子科技集团公司第二十四研究所 | 一种键合金丝可靠性评价方法 |
CN111707348A (zh) * | 2020-06-24 | 2020-09-25 | 中国电子产品可靠性与环境试验研究所((工业和信息化部电子第五研究所)(中国赛宝实验室)) | 光纤水听器寿命评价方法 |
CN112034282A (zh) * | 2020-08-19 | 2020-12-04 | 湖南海智机器人技术有限公司 | 一种基于加速因子的运动控制器加速退化试验方法 |
CN113032998A (zh) * | 2021-03-26 | 2021-06-25 | 北京中关村水木医疗科技有限公司 | 医疗器械寿命评估方法和装置 |
CN113032999A (zh) * | 2021-03-26 | 2021-06-25 | 北京中关村水木医疗科技有限公司 | 医疗设备使用寿命的预测方法及装置 |
CN113532497A (zh) * | 2021-07-12 | 2021-10-22 | 广州广电计量检测股份有限公司 | 一种传感器寿命评估方法及装置 |
CN114266006A (zh) * | 2021-12-17 | 2022-04-01 | 北京航空航天大学 | 一种加速退化试验测量不确定度的评定方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102708306A (zh) * | 2012-06-19 | 2012-10-03 | 华北电网有限公司计量中心 | 一种智能电表可靠寿命的预测方法 |
CN104076224A (zh) * | 2014-07-11 | 2014-10-01 | 国家电网公司 | 一种用电信息采集设备可靠性验证的试验方法 |
CN104181457B (zh) * | 2014-08-15 | 2017-01-18 | 中国电子科技集团公司第二十四研究所 | 一种半导体器件温湿度复合应力加速模型优选方法 |
-
2018
- 2018-05-28 CN CN201810521253.1A patent/CN108549047A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102708306A (zh) * | 2012-06-19 | 2012-10-03 | 华北电网有限公司计量中心 | 一种智能电表可靠寿命的预测方法 |
CN104076224A (zh) * | 2014-07-11 | 2014-10-01 | 国家电网公司 | 一种用电信息采集设备可靠性验证的试验方法 |
CN104181457B (zh) * | 2014-08-15 | 2017-01-18 | 中国电子科技集团公司第二十四研究所 | 一种半导体器件温湿度复合应力加速模型优选方法 |
Non-Patent Citations (5)
Title |
---|
杨舟: "基于回归分析的智能电表可靠寿命评估方法", 《工业技术创新》 * |
查国清 等: "基于多应力加速试验方法的智能电表寿命评估", 《北京航空航天大学学报》 * |
罗冉冉: "单向智能电能表可靠性加速寿命试验研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
贾占强 等: "基于步进加速退化试验的电子产品可靠性评估技术", 《系统工程理论与实践》 * |
陈金涛 等: "智能电能表的加速退化试验方案研究", 《电测与仪表》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110068789B (zh) * | 2019-05-23 | 2021-06-25 | 华立科技股份有限公司 | 电能表寿命的监测方法和电路 |
CN110068789A (zh) * | 2019-05-23 | 2019-07-30 | 华立科技股份有限公司 | 电能表寿命的监测方法和电路 |
CN111025224A (zh) * | 2020-01-07 | 2020-04-17 | 中国电力科学研究院有限公司 | 一种智能电表的加速退化试验方法 |
CN111351697A (zh) * | 2020-03-24 | 2020-06-30 | 中国电子科技集团公司第二十四研究所 | 一种键合金丝可靠性评价方法 |
CN111707348A (zh) * | 2020-06-24 | 2020-09-25 | 中国电子产品可靠性与环境试验研究所((工业和信息化部电子第五研究所)(中国赛宝实验室)) | 光纤水听器寿命评价方法 |
CN112034282A (zh) * | 2020-08-19 | 2020-12-04 | 湖南海智机器人技术有限公司 | 一种基于加速因子的运动控制器加速退化试验方法 |
CN113032998A (zh) * | 2021-03-26 | 2021-06-25 | 北京中关村水木医疗科技有限公司 | 医疗器械寿命评估方法和装置 |
CN113032999A (zh) * | 2021-03-26 | 2021-06-25 | 北京中关村水木医疗科技有限公司 | 医疗设备使用寿命的预测方法及装置 |
CN113032999B (zh) * | 2021-03-26 | 2024-05-17 | 北京中关村水木医疗科技有限公司 | 医疗设备使用寿命的预测方法及装置 |
CN113532497A (zh) * | 2021-07-12 | 2021-10-22 | 广州广电计量检测股份有限公司 | 一种传感器寿命评估方法及装置 |
CN113532497B (zh) * | 2021-07-12 | 2024-04-09 | 广电计量检测集团股份有限公司 | 一种速度传感器寿命评估方法及装置 |
CN114266006A (zh) * | 2021-12-17 | 2022-04-01 | 北京航空航天大学 | 一种加速退化试验测量不确定度的评定方法 |
CN114266006B (zh) * | 2021-12-17 | 2024-06-11 | 北京航空航天大学 | 一种加速退化试验测量不确定度的评定方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108549047A (zh) | 一种基于加速退化试验的电能表寿命评估方法 | |
CN109444791B (zh) | 一种电容式电压互感器误差状态评估方法和系统 | |
CN108037460B (zh) | 一种批产锂离子电池容量实时评估方法 | |
CN103048607B (zh) | 一种基于给定阈值的数控成品电路板性能退化测评方法 | |
CN108120949A (zh) | 一种智能电能表加速退化试验方法 | |
CN208459884U (zh) | 一种电能表的加速退化试验系统 | |
CN110286303A (zh) | 一种基于bp神经网络的同轴电缆绝缘老化状态评估方法 | |
CN103487781A (zh) | 一种基于加速老化的电子式互感器的可靠性评估方法 | |
CN115015623B (zh) | 一种基于大数据的矿山设备运行安全评价系统 | |
CN109284933A (zh) | 一种基于数理统计的电子式互感器状态评估系统及方法 | |
CN109375143B (zh) | 一种确定智能电能表剩余寿命的方法 | |
CN107907764B (zh) | 一种适用于智能仪表智能特性验证的检测方法和系统 | |
CN110764040B (zh) | 一种用于确定自动化检定系统误差测量功能的方法及系统 | |
CN117826058A (zh) | 一种关口电容式电压互感器运行状态在线监测方法及装置 | |
CN116773239A (zh) | 一种智能燃气表控制器可靠性寿命预测方法 | |
CN107561452A (zh) | 一种蓄电池组健康度的评估方法 | |
CN116699499A (zh) | 一种电表精度及可靠性自动检测系统 | |
CN116224208A (zh) | 电能表的误差检测方法、系统及存储介质 | |
CN111157938B (zh) | 一种评价自动化检定系统计量过程能力的方法和系统 | |
CN105548939B (zh) | 一种电压互感器状态检验方法 | |
CN114879071A (zh) | 一种基于平均电压的锂离子电池老化模式在线诊断方法 | |
CN112068067A (zh) | 流水线计量检定装置的表位隔离变压器带载误差检测系统 | |
CN109388829B (zh) | 一种电子产品寿命测算方法 | |
CN112147409A (zh) | 低压配电线路在线阻抗检测方法及装置 | |
CU et al. | Analysis of accelerated reliability testing data of electronic component in combat vehicles |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180918 |