CN108399278B - 一种电子装备多因素加速因子计算方法 - Google Patents
一种电子装备多因素加速因子计算方法 Download PDFInfo
- Publication number
- CN108399278B CN108399278B CN201810069541.8A CN201810069541A CN108399278B CN 108399278 B CN108399278 B CN 108399278B CN 201810069541 A CN201810069541 A CN 201810069541A CN 108399278 B CN108399278 B CN 108399278B
- Authority
- CN
- China
- Prior art keywords
- stress
- failure
- acceleration
- electronic equipment
- damage
- 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.)
- Active
Links
- 230000001133 acceleration Effects 0.000 title claims abstract description 167
- 238000004364 calculation method Methods 0.000 title claims description 33
- 230000006378 damage Effects 0.000 claims abstract description 146
- 230000007246 mechanism Effects 0.000 claims abstract description 72
- 238000000034 method Methods 0.000 claims abstract description 30
- 230000002860 competitive effect Effects 0.000 claims abstract description 9
- 230000008878 coupling Effects 0.000 claims description 41
- 238000010168 coupling process Methods 0.000 claims description 41
- 238000005859 coupling reaction Methods 0.000 claims description 41
- 230000009471 action Effects 0.000 claims description 24
- 230000007797 corrosion Effects 0.000 claims description 23
- 238000005260 corrosion Methods 0.000 claims description 23
- 239000000956 alloy Substances 0.000 claims description 13
- 229910045601 alloy Inorganic materials 0.000 claims description 12
- 230000001186 cumulative effect Effects 0.000 claims description 12
- 238000009825 accumulation Methods 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 8
- 230000004913 activation Effects 0.000 claims description 6
- 230000001808 coupling effect Effects 0.000 claims description 6
- 230000014509 gene expression Effects 0.000 claims description 6
- 238000002474 experimental method Methods 0.000 claims description 3
- 125000004122 cyclic group Chemical group 0.000 claims description 2
- 238000012360 testing method Methods 0.000 abstract description 34
- 230000035882 stress Effects 0.000 description 119
- 230000000694 effects Effects 0.000 description 8
- 239000000463 material Substances 0.000 description 8
- 238000006243 chemical reaction Methods 0.000 description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 230000007613 environmental effect Effects 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 230000018109 developmental process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 239000002184 metal Substances 0.000 description 4
- 229910052751 metal Inorganic materials 0.000 description 4
- 239000007769 metal material Substances 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000006872 improvement Effects 0.000 description 3
- 238000010998 test method Methods 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- 208000027418 Wounds and injury Diseases 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000001351 cycling effect Effects 0.000 description 2
- 230000006353 environmental stress Effects 0.000 description 2
- 208000014674 injury Diseases 0.000 description 2
- 150000002500 ions Chemical class 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000002861 polymer material Substances 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000012356 Product development Methods 0.000 description 1
- 229910052581 Si3N4 Inorganic materials 0.000 description 1
- CZMRCDWAGMRECN-UGDNZRGBSA-N Sucrose Chemical compound O[C@H]1[C@H](O)[C@@H](CO)O[C@@]1(CO)O[C@@H]1[C@H](O)[C@@H](O)[C@H](O)[C@@H](CO)O1 CZMRCDWAGMRECN-UGDNZRGBSA-N 0.000 description 1
- 229930006000 Sucrose Natural products 0.000 description 1
- 239000003929 acidic solution Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 239000012670 alkaline solution Substances 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000001311 chemical methods and process Methods 0.000 description 1
- 229910052681 coesite Inorganic materials 0.000 description 1
- 229910052906 cristobalite Inorganic materials 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000003487 electrochemical reaction Methods 0.000 description 1
- 238000005530 etching Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 230000009477 glass transition Effects 0.000 description 1
- RVRCFVVLDHTFFA-UHFFFAOYSA-N heptasodium;tungsten;nonatriacontahydrate Chemical compound O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.[Na+].[Na+].[Na+].[Na+].[Na+].[Na+].[Na+].[W].[W].[W].[W].[W].[W].[W].[W].[W].[W].[W] RVRCFVVLDHTFFA-UHFFFAOYSA-N 0.000 description 1
- 230000003301 hydrolyzing effect Effects 0.000 description 1
- 238000013095 identification testing Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 150000002739 metals Chemical class 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000000377 silicon dioxide Substances 0.000 description 1
- 235000012239 silicon dioxide Nutrition 0.000 description 1
- 229910000679 solder Inorganic materials 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 229910052682 stishovite Inorganic materials 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000005720 sucrose Substances 0.000 description 1
- 238000002076 thermal analysis method Methods 0.000 description 1
- 230000000930 thermomechanical effect Effects 0.000 description 1
- 229910052905 tridymite Inorganic materials 0.000 description 1
- 238000003466 welding Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种电子装备多因素加速因子计算方法,包括:针对影响电子设备寿命的不同因素,分别构建单点应力损伤模型并且计算得到对应的加速因子;基于电子装备的特点以及单点应力损伤模型,构建整个电子装备基于竞争失效的故障模型;基于构建的竞争失效故障模型,根据统计等效原则计算得到电子装备的加速因子。该方法针对于当前设备结构及失效模式复杂的电子装备,通过在温度、湿度、振动及电应力及其组合等各种应力条件下的失效机理和失效形式进行分析并构建相应的应力失效模型,能够根据不同因素及故障模型,相应根据系统的寿命模型得到电子装备的加速因子,完全适用于整机装备的加速试验,进而使得电子装备的加速试验更加准确、可靠。
Description
技术领域
本发明涉及设备加速试验相关技术领域,特别是指一种电子装备多因素加速因子计算方法。
背景技术
现代科学技术的发展和工业水平的提高,加快了电子装备的材料、元器件和工艺质量以及可靠性的改进速度,并且对可靠性试验技术也提出了新要求。因此,如何评估高可靠、长寿命产品的可靠性是目前可靠性工程领域亟待解决的重要问题。对于平均故障间隔时间(MTBF,Mean Time Between Failure)长达数千甚至上万小时的产品,进行传统的可靠性试验无论是从周期上,还是经费上都是无法承受的。所以,为了探索更加高效的方法以获取高可靠电子装备,迫切需要研究加速可靠性试验方法。
目前的加速试验主要为基于模型的试验方法,所用的模型包括物理模型和统计模型。但该方法的应用受到加速模型的限制,只适用于元器件及材料级产品。而对于电子装备,由于其包含大量的器件和材料,结构较为复杂。并且在其全寿命周期内经历的任务阶段及面临的环境因素多样,包括自然环境因素、诱发环境因素以及通电工作应力的共同作用,并且不同部件的环境敏感应力不同,导致多因素综合作用对电子装备的失效模式和失效机理的影响更为复杂。为满足装备的研制要求,开展多因素加速可靠性试验是一个必然发展趋势,其中的一个关键技术是解决电子装备在温度、湿度、振动、电应力及其组合模型等多应力条件下的加速因子计算难题。
因此,在实现本申请的过程中,发明人发现现有技术至少存在以下问题:当前的加速寿命试验针对结构以及寿命影响因素复杂的电子设备难以获得准确有效的加速因子。
发明内容
有鉴于此,本发明的目的在于提出一种电子装备多因素加速因子计算方法,能够根据失效机理以及应力作用方式的不同,提出不同的应力损伤模型,进而计算得到准确有效的加速因子,为后续相关的加速试验提供可靠的依据和前提。
基于上述目的本发明提供的一种电子装备多因素加速因子计算方法,包括:
针对影响电子设备寿命的不同因素,分别构建单点应力损伤模型并且计算得到对应的加速因子;
基于电子装备的特点以及单点应力损伤模型,构建整个电子装备基于竞争失效的故障模型;
基于构建的竞争失效故障模型,根据统计等效原则计算得到电子装备的加速因子。
可选的,所述构建单点应力损伤模型并且计算得到对应的加速因子的步骤包括:
若失效是由于单因素单失效机理的应力导致的,则需要计算单应力条件下的加速因子;其中,
对于温度应力,通过阿伦尼斯模型得到温度应力对应的加速模型;
对于温度循环应力,利用Manson-Coffin模型计算得到加速因子;
对于振动应力,通过逆幂律模型计算得到加速因子;
对于腐蚀模型,通过法拉第定理及腐蚀损伤相等原则,计算得到加速因子。
可选的,所述构建单点应力损伤模型并且计算得到对应的加速因子的步骤包括:
若失效是由单点的不同失效机理在独立条件下导致的,则所述构建单点多应力失效机理独立损伤模型并且计算得到对应的加速因子的步骤包括:
采用miner损伤累积法则构建不同失效机理对该单点的累积损伤模型;计算表达式如下:
在多应力独立作用下,计算单点多应力独立累积损伤下的第一单点寿命值;
在多应力独立加速条件下,计算单点多应力独立累积损伤在加速条件下的第二单点寿命值;
将第一单点寿命值除以第二单点寿命值,计算得到单点多应力失效机理独立损伤加速因子。
可选的,当失效形式为电子装备互联部位的热疲劳失效和振动疲劳失效,得到累积损伤模型,计算表达式为:
D=Dc+Df=1;
其中,Tc为热循环应力c单独作用下的应力损伤寿命;Tf为振动应力f单独作用下的应力损伤寿命;
在双应力独立作用下,互联部位的应力损伤寿命t为:
在多应力加速条件下,热循环应力c及振动应力f对应的加速因子分别为Ac,Af,则互联部位等效加速寿命tA为:
对于互联部位热循环应力和振动应力累积损伤失效情况,其加速因子计算式为:
可选的,所述构建单点应力损伤模型并且计算得到对应的加速因子的步骤包括:
若失效是由单点的不同失效机理相互耦合导致的,则所述构建单点多应力失效机理耦合损伤模型并且计算得到对应的加速因子的步骤包括:
根据单点多应力失效机理耦合条件下的应力损伤机理,构建耦合条件下的累积损伤模型;
在多应力耦合作用下,计算单点多应力耦合累积损伤下的第三单点寿命值;
在多应力耦合加速条件下,计算单点多应力耦合累积损伤在加速条件下的第四单点寿命值;
将第三单点寿命值除以第四单点寿命值,计算得到单点多应力失效机理耦合损伤加速因子。
可选的,当失效形式为电子装备合金结构的疲劳-蠕变失效时,耦合条件下的累积损伤模型为:
D=Dc+Df=Dcr;
其中:Tc为蠕变应力c作用下的应力损伤寿命;Tf为疲劳应力f作用下的应力损伤寿命;Dcr表示临界总损伤,当D>Dcr时发生失效;当D<Dcr结构处于安全状态;基于最简单的累积损伤理论认为总损伤达到1时结构失效,但实验发现临界损伤并不等于1,其数值与损伤应力的施加顺序有关;
在双应力耦合作用下,电子装备合金结构的应力损伤寿命t为:
并设在多应力耦合加速条件下,蠕变应力c及疲劳应力f对应的加速因子分别为Ac,Af,则合金结构等效加速寿命tA为:
对于电子装备合金结构疲劳-蠕变应力耦合累积损伤失效情况,其加速因子计算式为:
可选的,所述构建基于竞争失效的故障模型的步骤还包括:
基于竞争失效模型,在已知组成电子装备部件、器件等单点寿命分布及其参数的情况下,综合底层信息计算电子装备整机的寿命分布及其参数。
可选的,所述构建基于竞争失效的故障模型的步骤还包括:
通过构建的竞争失效故障模型,基于电子装备在加速状态下的失效与自然贮存失效的等效关系,获取其的加速因子;
对于电子装备,其组成部件、器件的寿命分布为逆高斯分布,其平均无故障工作时间为:Ti=μi;系统的平均寿命计算式为:
在加速条件下,系统的等效加速寿命为:
其中,μA为电子装备在多应力加速条件下的等效加速寿命;A为电子装备的加速因子;μS表示电子装备实际使用条件下平均寿命;μi表示电子装备所包含失效单点i下的平均寿命;Ai为失效单点i多应力加速条件加速因子,n为电子装备所包含的失效单点数量;
计算得出逆高斯(IG)寿命分布电子装备的加速因子为:
从上面所述可以看出,本发明提供的电子装备多因素加速因子计算方法,针对于当前结构及失效模式复杂的电子装备,通过在温度、湿度、振动以及电应力及其组合等各种应力条件下的失效机理和形式进行分析和构建相应的应力失效模型,进而能够根据不同的因素以及相应的故障模型,相应的根据系统的寿命模型计算得到电子装备的加速因子,使得计算得到的加速因子完全适用于整机装备的加速试验,从而克服不同部件之间的失效差异,进而使得电子装备的加速试验更加准确、可靠。
附图说明
图1为本发明提供的电子装备失效机理及失效模式的相互关系的示意图;
图2为本发明提供的电子装备多因素加速因子计算方法的一个实施例的流程图;
图3为本发明提供的单点多应力失效机理独立条件下累积损伤加速因子计算的流程图;
图4为本发明提供的单点多应力失效机理耦合条件下累积损伤加速因子计算的流程图;
图5为本发明提供的电子装备多应力加速因子计算综合方案示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
需要说明的是,本发明实施例中所有使用“第一”、“第二”、“第三”、“第四”的表述均是为了区分两个相同名称非相同的实体或者非相同的参量,可见“第一”“第二”“第三”“第四”仅为了表述的方便,不应理解为对本发明实施例的限定,后续实施例对此不再一一说明。
本申请针对当前加速试验中整机设备例如整机电子设备难以准确计算得到加速因子的问题,提出一种多因素情况下计算加速因子的方法。而为了更加清楚的理解本申请方案的特点及效果,首先需要分析设备失效机理以及相应损伤模型。参照图1所示,为本发明提供的电子装备失效机理及失效模式的相互关系的示意图。由图可知,电子装备失效的应力、失效模式及失效机理之间呈现如下关联映射关系:应力→失效机理→失效模式。这种关联映射关系不是单纯的直线关系,例如:敏感应力2经由失效机理M导致失效模式2发生;敏感应力1经由失效机理1导致失效模式1发生;而应力2和应力M共同作用经由失效机理2同样也可导致失效模式1发生。所以电子设备最终的失效是由多种可能的因素导致的。
而为进行电子装备多应力加速因子计算,首选需要建立多因素应力损伤模型。申请人发现,典型电子装备多应力损伤可分为单点应力损伤模型和系统应力损伤模型。对于单点损伤模型,需要分为两种情况:1)失效点的失效是由于单因素单失效机理的应力损伤作用;2)失效点的失效是多因素多失效机理的应力损伤作用。对于单点多因素应力损伤,又可分为两种情况:1)如果假设各种失效机理之间互相独立,不发生相互影响,需要采用损伤累计方法来处理;2)如果考虑不同失效机理的耦合作用,需要采用考虑耦合项的多因素应力累积损伤加速因子综合方法。对于电子装备系统,由于考虑的是多个失效点的应力损伤的综合,需要采用故障竞争方法来处理。
参照图2所示,为本发明提供的电子装备多因素加速因子计算方法的一个实施例的流程图。所述电子装备多因素加速因子计算方法包括:
步骤101,针对影响电子设备寿命的不同因素,分别构建单点应力损伤模型并且计算得到对应的加速因子;也即任何系统首先是由单点构成,所以先构建并计算单点的加速因子,而且需要针对不同的影响因素分析分析计算。
步骤102,基于电子装备的特点以及单点应力损伤模型,构建整个电子装备基于竞争失效的故障模型;通过竞争失效的模式使得整个电子装备中的单点相互关联起来,然后得到适用于整个电子装备的故障模型。
步骤103,基于构建的竞争失效故障模型,根据统计等效原则计算得到电子装备的加速因子。
由上述实施例可知,所述电子装备多因素加速因子计算方法针对于当前结构及失效模式复杂的电子装备,通过在温度、湿度、振动以及电应力及其组合等各种应力条件下的失效机理和形式进行分析和构建相应的应力失效模型,进而能够根据不同的因素以及相应的故障模型,相应的根据系统的寿命模型计算得到电子装备的加速因子,使得计算得到的加速因子完全适用于整机装备的加速试验,从而克服不同部件之间的失效差异,进而使得电子装备的加速试验更加准确、可靠。
在本申请一些可选的实施例中,若失效点的失效是由于单因素单失效机理的应力损伤作用导致的,则只需开展各单应力条件下的加速因子测定即可。其中,典型的单因素应力损伤机理包括高分子材料的老化、焊点的热疲劳、结构的振动疲劳及金属材料的腐蚀等。因此,本申请以温度应力、温度循环应力、振动应力以及腐蚀应力为例进行计算说明。
其中,温度应力对应加速因子的测定包括:瑞典物理学家Svante Arrhenius研究了温度激发类化学过程,在研究温度对酸催化蔗糖水解转化反应的基础上总结出:某产品的性能退化速率与激活能的指数成反比,与温度倒数的指数成反比,在大量数据的基础上总结出了著名的经验公式,即阿伦尼斯模型:
ξ=AeE/KT (1)
其中,ξ为产品寿命特征,如平均寿命、中位寿命、特征寿命等;A是未知常数;E表示产品激活能,与材料有关;K为波尔兹曼常数,值为8.617×10-5ev/K;T为绝对温度,单位为K。
阿伦尼斯模型主要用于以温度为加速应力的加速寿命试验的数据分析,可以将模型两边取对数,进一步得到加速模型
lnξ=a+bT (2)
其中,a=lnA,b=E/K都是未知参数,需要进行估计或者根据经验设置。
阿伦尼斯模型表明寿命特征的对数是温度倒数的线性函数。该模型在电子产品的加速寿命试验中得到了广泛的应用。但是,这种模型只考虑了温度应力对物质的化学及物理性质变化的影响,因此仅适用于单应力加速。
可选的,除了通过加速寿命试验获得产品的激活能外,还可以采用基于热分析技术的测试得到产品的激活能。例如,利用动态热机械分析仪得到的高分子材料试件的测试频率f和相应的玻璃化转变温度Tg,两者的关系遵循WLF方程:
其中,R为理想气体常数,取值为8.345J/mol/K;c为常数。将f、Tg等试验数据带入进行回归分析,便可计算得到激活能E。
在另一些可选的实施例中,针对于温度循环应力,其加速因子的测定方法包括:
温度循环最常用的加速模型是Coffin-Manson模型,其加速因子是
其中,ΔT1,ΔT2为不同的温度变化幅值,q是与失效机理有关的常数,常见的q值见表1所示。
表1材料加速系数表
材料 | q |
柔性材料、焊锡 | 1~3 |
硬的合金材料 | 3~5 |
易碎的材料,Si,SiO2,Si3N4 | 6~9 |
进一步,由于Coffin-Manson的温度循环加速因子只考虑了温度变化幅值,而未考虑温度变化速率,因此,可以对Coffin-Manson模型进行修正。
其中,T为温度;F为使用频率;n一般取1.9;m为0.3333,Tmax1、Tmax2为温度循环中最高温度。
在一些可选的实施例中,振动应力加速因子的测定方法包括:
加速振动应力的选取原则是在加严应力下产品出现的故障模式和机理与正常应力下一致,即不会因为振动应力量级的加大而导致新的失效模式和机理产生。在产品研制中,环境鉴定试验即产品的功能和耐久振动试验是产品经受的最大振动量级,在加速试验方案中振动应力的选取应保证不产生故障模式和机理的变化以及过应力试验。
振动应力的加严是通过统一提高功率谱密度曲线来实现的,即不改变振动能量的频率分布。
MIL-STD-810G中的振动疲劳等价折算公式与GJB1032中的振动等价折算公式的表达式是:
式中,W和W'分别是折算前和折算后的振动功率谱密度,T和T'分别是折算前和折算后的振动时间。
在一些可选的实施例中,腐蚀应力加速因子的测定方法包括:
基于绝大多数金属腐蚀是由电化学反应引起的,金属在水中阳极反应的通式为:
与此同时,在金属上留下等量的电子从阳极流向阴极,中性或碱性溶液中生成
在酸性溶液中反应为:
有上式可知,腐蚀总是伴随着离子迁移。电化学作用时间越长,离子迁移越多,腐蚀越严重。金属电化腐蚀过程需要适当厚度的水膜。水膜只有几个分子厚时,腐蚀速度极小;当水膜厚度增加为100A~1μm时,腐蚀速度逐渐增加;当水膜厚度为1μm~1mm时腐蚀速度增至最大;当水膜厚增至几百微米时,由于氧扩散变得困难,腐蚀速度下降。
利用电化学原理如Faraday定律,按腐蚀损伤相等原则,将不同温度、湿度、介质浓度,折算到某一试验室便于实现的条件上,起到加速腐蚀作用。对应的计算公式为:
式中,v表示腐蚀速率,um/a;ik表示腐蚀电流密度,单位为uA/cm2;M为被氧化元素的摩尔质量,单位为g/mol;d表示被氧化元素的密度。
通过测定不同温度、湿度下典型金属材料腐蚀电流来建立不同组合对应的当量折算系数。最终计算得到加速腐蚀速率,也即加速因子。
参照图3所示,为本发明提供的单点多应力失效机理独立条件下累积损伤加速因子计算的流程图。若失效是由单点的不同失效机理在独立条件下导致的,则所述构建单点多应力失效机理独立损伤模型并且计算得到对应的加速因子的步骤包括:
步骤201,采用miner损伤累积法则构建不同失效机理对该单点的累积损伤模型;
步骤202,在多应力独立作用下,计算单点多应力独立累积损伤下的第一单点寿命值;
步骤203,在多应力独立加速条件下,计算单点多应力独立累积损伤在加速条件下的第二单点寿命值;
步骤204,将第一单点寿命值除以第二单点寿命值,计算得到单点多应力失效机理独立损伤加速因子。
具体的,若失效点的失效是多因素多失效机理的应力损伤作用,并且各种失效机理之间互相独立,不发生相互影响,可以采用损伤累计方法来计算失效点的加速因子。典型的单点多应力失效机理独立条件下的应力损伤机理包括电子装备互联部位的热疲劳失效和振动疲劳失效两种机理造成的损伤等。
在多应力损伤作用下,电子装备的失效过程可以看作是损伤程度趋于临界损伤值的累积过程,不同环境应力都会对电子装备产生影响,都会对产品的最终失效做出“贡献”,而且这一点点的损伤不断累积起来,直至电子产品最终失效。如果认为每一环境应力对电子产品的损伤量只与最终损伤量的大小有关,而与应力本身无关,则最终仅考虑损伤量的叠加效应,这就是损伤累积假设。
对于同一失效点发生的不同失效机理,需要考虑各种应力对该位置点造成的综合损伤,进行损伤累积,获得该潜在失效点的总损伤量。具有代表性的是Miner的线性累积损伤法则,参见式(21)。损伤累积分析的结果是产品在给定时间内,各潜在点的累积损伤量。当累积损伤量达到1时,失效就会发生。因此,可以进一步计算得到该产品中各潜在失效点的失效时间。
其中,D为多应力损伤量,l为应力的种类数,ti为第i种应力损耗的寿命,Ti为第i种应力的寿命。
以电子装备互联部位的热疲劳失效和振动疲劳失效为例,其累积损伤模型表述为:
D=Dc+Df=1 (22)
其中:Tc为热循环应力c单独作用下的应力损伤寿命;Tf为振动应力f单独作用下的应力损伤寿命。因此,在双应力独立作用下,互联部位的应力损伤寿命t为:
在多应力加速条件下,热循环应力c及振动应力f对应的加速因子分别为Ac,Af,则互联部位等效加速寿命tA为:
因此,对于互联部位热循环应力和振动应力累积损伤失效情况,其加速因子计算式为:
参照图4所示,为本发明提供的单点多应力失效机理耦合条件下累积损伤加速因子计算的流程图。若失效是由单点的不同失效机理相互耦合导致的,则所述构建单点多应力失效机理耦合损伤模型并且计算得到对应的加速因子的步骤包括:
步骤301,根据单点多应力失效机理耦合条件下的应力损伤机理,构建耦合条件下的累积损伤模型;
步骤302,在多应力耦合作用下,计算单点多应力耦合累积损伤下的第三单点寿命值;
步骤303,在多应力耦合加速条件下,计算单点多应力耦合累积损伤在加速条件下的第四单点寿命值;
步骤304,将第三单点寿命值除以第四单点寿命值,计算得到单点多应力失效机理耦合损伤加速因子。
具体的,若失效点的失效是多因素多失效机理的应力损伤作用,并且各种失效机理具有相互耦合的作用,则不能简单用各个应力造成的损伤简单叠加来描述,进而需要研究耦合失效的问题,计算多应力失效机理耦合条件下单点的加速因子。典型的单点多应力失效机理耦合条件下的应力损伤机理包括电子装备金属材料的腐蚀-疲劳失效,以及疲劳-蠕变失效等。
以电子装备合金结构的疲劳-蠕变失效为例,可以将耦合条件下的累积损伤模型表述为:
D=Dc+Df=Dcr (27)
其中:Tc为蠕变应力c作用下的应力损伤寿命;Tf为疲劳应力f作用下的应力损伤寿命。
其中,Dcr表示临界总损伤,当D>Dcr时发生失效;D<Dcr结构处于安全状态。最简单的累积损伤理论认为总损伤达到1时结构失效,但是实验发现临界损伤并不等于1,其数值与损伤应力的施加顺序有关。
因此,在双应力耦合作用下,电子装备合金结构的应力损伤寿命t为:
在多应力耦合加速条件下,蠕变应力c及疲劳应力f对应的加速因子分别为Ac,Af,则合金结构等效加速寿命tA为:
因此,对于电子装备合金结构疲劳-蠕变应力耦合累积损伤失效情况,其加速因子计算式为:
可以看出,单点多应力失效机理耦合条件下累积损伤加速因子计算式,与单点多应力失效机理耦合条件下累积损伤加速因子综合计算式相同。
在本申请一些可选的实施例中,所述构建基于竞争失效的故障模型的步骤还包括:
基于竞争失效模型,在已知组成电子装备部件、器件等单点寿命分布及其参数的情况下,综合底层信息计算电子装备整机的寿命分布及其参数。
其中,对于电子装备,由于其包含大量的非金属材料,电子元器件,结构较为复杂。因此,电子装备可以认为是由多个单点组成的一个系统。对于包含多个失效单点的电子装备系统,可以认为各个单点的失效之间是竞争关系,在相同的工作环境条件作用下,具有最短失效时间的失效单点的失效时间就是电子装备的寿命。
竞争性故障模型可定义为,如果系统有K种失效方式,而每一种失效方式都独立的作用于系统,且都对应一定的失效时间,其中任何一种失效都会引起系统失效,在所有的失效中,最早产生的那种失效出现时,将导致系统失效,即失效因素T:
T=min{T1,T2,...,Tn} (31)
任一个失效因素起作用时,对应的可靠度为
其中,λi(t)是对应第i个失效因素的失效率。
当n个因素同时起作用时,系统的可靠度为:
整机的总失效率将是对应时刻t的n个独立的失效率之和,即
λ(t)=λ1(t)+λ2(t)+...+λn(t) (34)
对于电子装备,其组成部件、器件的寿命分布为逆高斯分布,其平均无故障工作时间为:Ti=μi;系统的平均寿命计算式为:
在加速条件下,系统的等效加速寿命为:
其中,μA为电子装备在多应力加速条件下的等效加速寿命;A为电子装备的加速因子;μS表示电子装备实际使用条件下平均寿命;μi表示电子装备所包含失效单点i下的平均寿命;Ai为失效单点i多应力加速条件加速因子,n为电子装备所包含的失效单点数量;
由上式,计算得出逆高斯寿命分布电子装备的加速因子为:
也即最终计算得到电子设备的加速因子。
参照图5所示,为本发明提供的电子装备多应力加速因子计算综合方案示意图。本申请针对于电子装备在其全寿命周期内经历的任务阶段及面临的环境因素多样的特点,在单点单因素应力损伤加速因子计算方法的技术上,首次提出了单点多应力失效机理独立条件下累积损伤加速因子计算方法及单点多应力失效机理耦合条件下累积损伤加速因子计算方法;同时,针对目前加速试验由于受到加速模型的限制,只适用于元器件及材料级产品的问题,本申请提出了基于竞争失效的系统多应力加速因子综合计算方法,为开展结构复杂电子装备整机加速试验奠定基础。因此,本申请至少包括以下优点:(1)对于MTBF长达数千小时的产品,进行传统的可靠性试验无论是从周期上,还是经费上都是无法承受的。本申请提出的一种电子装备多因素加速因子综合计算方法,为开展电子装备的可靠性加速试验奠定了基础,在提高可靠性试验效果的同时短试验时间、降低试验费用,具有很高的工程应用价值。(2)通过计算电子装备在温度、湿度、振动、电应力及其组合模型等多应力条件下的加速因子,可以为开展电子装备整机加速试验奠定基础。克服了目前基于模型的加速寿命试验方法需要产品结构简单、具有单一失效模式及失效机理,只适用于元器件及材料级产品的限制。电子装备整机加速试验中各器件及材料的边界条件更为真实,因而试验结果更为真实可信。
所属领域的普通技术人员应当理解:以上任何实施例的讨论仅为示例性的,并非旨在暗示本公开的范围(包括权利要求)被限于这些例子;在本发明的思路下,以上实施例或者不同实施例中的技术特征之间也可以进行组合,步骤可以以任意顺序实现,并存在如上所述的本发明的不同方面的许多其它变化,为了简明它们没有在细节中提供。
本发明的实施例旨在涵盖落入所附权利要求的宽泛范围之内的所有这样的替换、修改和变型。因此,凡在本发明的精神和原则之内,所做的任何省略、修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种电子装备多因素加速因子计算方法,其特征在于,包括:
针对影响电子设备寿命的不同因素,分别构建单点应力损伤模型并且计算得到对应的加速因子;
所述构建单点应力损伤模型并且计算得到对应的加速因子的步骤包括:
若失效是由于单因素单失效机理的应力导致的,则需要计算单应力条件下的加速因子;其中,
对于温度应力,通过阿伦尼斯模型得到温度应力对应的加速模型;所述加速模型表示为lnξ=a+bT,其中,a=lnA,b=E/K,ξ为产品寿命特征,A是常数,E表示产品激活能,K为波尔兹曼常数,值为8.617×10-5ev/K,T为绝对温度;
对于温度循环应力,利用Coffin-Manson模型计算得到加速因子;对所述Coffin-Manson模型进行修正,修正后模型表示为 其中,T为温度,ΔT1为第一温度变化幅值,ΔT2为第二温度变化幅值,F为使用频率,F1为电子装备的第一使用频率,F2为电子装备的第二使用频率,n和m为常数,Tmax 1和Tmax 2为温度循环中的两个最高温度;
对于振动应力,通过逆幂律模型计算得到加速因子;
若失效是由单点的不同失效机理在独立条件下导致的,则构建单点多应力失效机理独立损伤模型并且计算得到对应的加速因子,包括:
采用miner损伤累积法则构建不同失效机理对该单点的累积损伤模型;计算表达式如下:
在多应力独立作用下,计算单点多应力独立累积损伤下的第一单点寿命值;
在多应力独立加速条件下,计算单点多应力独立累积损伤在加速条件下的第二单点寿命值;
将第一单点寿命值除以第二单点寿命值,计算得到单点多应力失效机理独立损伤加速因子;
当失效形式为电子装备互联部位的热疲劳失效和振动疲劳失效,得到累积损伤模型,计算表达式为:
D=Dc+Df=1;
其中,Tc为热循环应力c单独作用下的应力损伤寿命;Tf为振动应力f单独作用下的应力损伤寿命;
在双应力独立作用下,互联部位的应力损伤寿命t为:
在多应力加速条件下,热循环应力c及振动应力f对应的加速因子分别为Ac,Af,则互联部位等效加速寿命tA为:
对于互联部位热循环应力和振动应力累积损伤失效情况,其加速因子计算式为:
若失效是由单点的不同失效机理相互耦合导致的,则构建单点多应力失效机理耦合损伤模型并且计算得到对应的加速因子,包括:
根据单点多应力失效机理耦合条件下的应力损伤机理,构建耦合条件下的累积损伤模型;
在多应力耦合作用下,计算单点多应力耦合累积损伤下的第三单点寿命值;
在多应力耦合加速条件下,计算单点多应力耦合累积损伤在加速条件下的第四单点寿命值;
将第三单点寿命值除以第四单点寿命值,计算得到单点多应力失效机理耦合损伤加速因子;
当失效形式为电子装备合金结构的疲劳-蠕变失效时,耦合条件下的累积损伤模型为:
D=Dc+Df=Dcr;
其中:Tc为蠕变应力c作用下的应力损伤寿命;Tf为疲劳应力f作用下的应力损伤寿命;Dcr表示临界总损伤,当D>Dcr时发生失效;当D<Dcr结构处于安全状态;基于最简单的累积损伤理论认为总损伤达到1时结构失效,但实验发现临界损伤并不等于1,其数值与损伤应力的施加顺序有关;
在双应力耦合作用下,电子装备合金结构的应力损伤寿命t为:
并设在多应力耦合加速条件下,蠕变应力c及疲劳应力f对应的加速因子分别为Ac,Af,则合金结构等效加速寿命tA为:
对于电子装备合金结构疲劳-蠕变应力耦合累积损伤失效情况,其加速因子计算式为:
基于电子装备的特点以及单点应力损伤模型,构建整个电子装备基于竞争失效的故障模型;
基于构建的竞争失效故障模型,根据统计等效原则计算得到电子装备的加速因子。
2.根据权利要求1所述的方法,其特征在于,所述构建基于竞争失效的故障模型的步骤还包括:
通过构建的竞争失效故障模型,基于电子装备在加速状态下的失效与自然贮存失效的等效关系,获取电子装备的加速因子;
对于电子装备,其组成部件、器件的寿命分布为逆高斯分布,电子装备平均无故障工作时间为:Ti=μi;电子装备的平均寿命计算式为:
在加速条件下,电子装备的等效加速寿命为:
其中,μA为电子装备在多应力加速条件下的等效加速寿命;A为电子装备的加速因子;μS表示电子装备实际使用条件下平均寿命;μi表示电子装备所包含失效单点i下的平均寿命;Ai为失效单点i多应力加速条件加速因子,n为电子装备所包含的失效单点数量;
计算得出逆高斯寿命分布电子装备的加速因子为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810069541.8A CN108399278B (zh) | 2018-01-24 | 2018-01-24 | 一种电子装备多因素加速因子计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810069541.8A CN108399278B (zh) | 2018-01-24 | 2018-01-24 | 一种电子装备多因素加速因子计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108399278A CN108399278A (zh) | 2018-08-14 |
CN108399278B true CN108399278B (zh) | 2021-11-30 |
Family
ID=63094206
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810069541.8A Active CN108399278B (zh) | 2018-01-24 | 2018-01-24 | 一种电子装备多因素加速因子计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108399278B (zh) |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109655683B (zh) * | 2018-12-13 | 2021-07-20 | 广州广电计量检测股份有限公司 | 一种汽车电子产品热疲劳寿命加速测试方法 |
CN111380771B (zh) * | 2019-04-30 | 2022-11-08 | 中国人民解放军63961部队 | 具有两阶段失效机制的长期贮存产品贮存寿命的测定方法 |
EP3757703A1 (en) * | 2019-06-27 | 2020-12-30 | Siemens Aktiengesellschaft | System, apparatus and method of condition based management of one or more electro-mechanical systems |
CN110399658B (zh) * | 2019-07-09 | 2021-01-22 | 湖北文理学院 | 电池的加速因子值计算方法、装置、设备及存储介质 |
CN110414086B (zh) * | 2019-07-10 | 2023-01-17 | 北京华安中泰检测技术有限公司 | 一种基于灵敏度的综合应力加速因子计算方法 |
CN110414087B (zh) * | 2019-07-10 | 2023-01-24 | 北京华安中泰检测技术有限公司 | 一种多应力耦合作用下的电子产品寿命评估方法 |
CN111125868B (zh) * | 2019-11-08 | 2021-06-15 | 航天科工防御技术研究试验中心 | 一种电子产品的多应力寿命评估方法及装置 |
CN111079250B (zh) * | 2019-11-08 | 2021-06-08 | 航天科工防御技术研究试验中心 | 一种电子产品疲劳寿命评估、评估模型建立方法及装置 |
CN111079253B (zh) * | 2019-11-12 | 2023-12-12 | 中国电子产品可靠性与环境试验研究所((工业和信息化部电子第五研究所)(中国赛宝实验室)) | 电力电子器件寿命评估方法、装置和计算机设备 |
CN111079255A (zh) * | 2019-11-15 | 2020-04-28 | 湖南海智机器人技术有限公司 | 一种基于加速因子的电子调速器加速寿命试验方法 |
CN111475932B (zh) * | 2020-03-26 | 2023-05-16 | 青岛海尔空调器有限总公司 | 压缩机的测试方法及装置 |
CN111898236B (zh) * | 2020-05-25 | 2024-01-09 | 中国航天标准化研究所 | 基于失效大数据的加速贮存试验加速因子分析方法 |
CN111881539A (zh) * | 2020-05-25 | 2020-11-03 | 中国航天标准化研究所 | 基于失效大数据的电子整机加速贮存试验加速因子风险率分析方法 |
CN112464441B (zh) * | 2020-11-04 | 2023-06-30 | 北京强度环境研究所 | 一种电子产品多维向量加速因子表征方法 |
CN112784414B (zh) * | 2021-01-22 | 2022-06-10 | 北京航空航天大学 | 一种多部件整机贮存寿命置信下限评估方法 |
CN113378368B (zh) * | 2021-06-03 | 2022-07-22 | 中国人民解放军32181部队 | 一种基于非线性退化轨迹模型的加速因子评估方法 |
CN113609760B (zh) * | 2021-07-19 | 2023-12-05 | 北京圣涛平试验工程技术研究院有限责任公司 | 集成电路寿命预估方法及装置、电子设备和存储介质 |
CN114088117B (zh) * | 2021-11-30 | 2024-06-25 | 中国兵器工业集团第二一四研究所苏州研发中心 | 一种复杂工况下mems惯性器件可靠性评测方法 |
CN114239326B (zh) * | 2022-02-28 | 2022-08-02 | 中国电子产品可靠性与环境试验研究所((工业和信息化部电子第五研究所)(中国赛宝实验室)) | 产品可靠性加速系数评估方法、装置和计算机设备 |
CN115081200B (zh) * | 2022-06-13 | 2024-05-28 | 北京理工大学 | 复杂设备的加速因子及失效边界域分析方法 |
CN115468759B (zh) * | 2022-09-22 | 2024-01-23 | 中国核动力研究设计院 | 一种电动阀门热备寿命加速试验方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103336901A (zh) * | 2013-06-27 | 2013-10-02 | 中国人民解放军国防科学技术大学 | 一种竞争失效相关加速寿命试验统计分析方法 |
CN104298853A (zh) * | 2014-08-26 | 2015-01-21 | 北京航空航天大学 | 一种多应力多机理作用下柱塞泵产品综合加速因子优化方法 |
CN104392073A (zh) * | 2014-12-12 | 2015-03-04 | 中国航空综合技术研究所 | 一种基于故障物理的电子产品可靠性加速试验方法 |
CN105868543A (zh) * | 2016-03-25 | 2016-08-17 | 航天科工防御技术研究试验中心 | 基于逆高斯寿命分布的贮存寿命试验加速因子评估方法 |
CN106650204A (zh) * | 2016-09-27 | 2017-05-10 | 北京航空航天大学 | 一种产品故障行为耦合建模和可靠性评估方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090096405A1 (en) * | 2007-10-15 | 2009-04-16 | General Electric Company | Method and system for remotely predicting the remaining life of an ac motor system |
CN101710368B (zh) * | 2009-12-21 | 2011-09-28 | 北京航空航天大学 | 基于多源退化数据的贝叶斯可靠性综合评估方法 |
CN104142293B (zh) * | 2014-07-14 | 2016-08-17 | 山东科技大学 | 一种镁合金表面薄液膜电化学实验装置及其实验方法 |
CN106872575B (zh) * | 2017-01-05 | 2020-01-14 | 航天科工防御技术研究试验中心 | 一种塑封器件分层缺陷的分级风险评价方法 |
-
2018
- 2018-01-24 CN CN201810069541.8A patent/CN108399278B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103336901A (zh) * | 2013-06-27 | 2013-10-02 | 中国人民解放军国防科学技术大学 | 一种竞争失效相关加速寿命试验统计分析方法 |
CN104298853A (zh) * | 2014-08-26 | 2015-01-21 | 北京航空航天大学 | 一种多应力多机理作用下柱塞泵产品综合加速因子优化方法 |
CN104392073A (zh) * | 2014-12-12 | 2015-03-04 | 中国航空综合技术研究所 | 一种基于故障物理的电子产品可靠性加速试验方法 |
CN105868543A (zh) * | 2016-03-25 | 2016-08-17 | 航天科工防御技术研究试验中心 | 基于逆高斯寿命分布的贮存寿命试验加速因子评估方法 |
CN106650204A (zh) * | 2016-09-27 | 2017-05-10 | 北京航空航天大学 | 一种产品故障行为耦合建模和可靠性评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108399278A (zh) | 2018-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108399278B (zh) | 一种电子装备多因素加速因子计算方法 | |
Uckun et al. | Standardizing research methods for prognostics | |
Kumaraswamidhas et al. | Bearing degradation assessment and remaining useful life estimation based on Kullback-Leibler divergence and Gaussian processes regression | |
CN111859658B (zh) | 一种产品贮存寿命与可靠性评估方法 | |
Elwany et al. | Real-time estimation of mean remaining life using sensor-based degradation models | |
Vichare et al. | Environment and usage monitoringof electronic products for health assessment and product design | |
CN109829137A (zh) | 一种周期应力下非线性退化设备的寿命预测方法及系统 | |
CN106169124B (zh) | 系统级产品可靠性综合评估置信推断方法 | |
Cruz-Manzo et al. | Inductive effect on the fuel cell cathode impedance spectrum at high frequencies | |
Ismail et al. | Gaussian process regression remaining useful lifetime prediction of thermally aged power IGBT | |
Ismail et al. | Remaining useful life estimation for thermally aged power insulated gate bipolar transistors based on a modified maximum likelihood estimator | |
Ahmad et al. | Estimating the remaining useful life of bearings using a neuro-local linear estimator-based method | |
CN108846239A (zh) | 一种基于温湿度的弹性环氧树脂加速贮存试验及评估方法 | |
CN114388855B (zh) | 一种燃料电池寿命预测方法 | |
US20240110861A1 (en) | Multifunctional corrosion probe system | |
Quan et al. | Time-resolved water measurement in a PEM fuel cell using high-resolution neutron imaging technique | |
CN111079270B (zh) | 一种基于二元混合随机过程的轴承剩余寿命预测方法 | |
CN115099044A (zh) | 一种涂层寿命预测方法、系统、计算机设备和存储介质 | |
Park et al. | Evaluation of silver-coated magnesium bipolar plate for lightweight PEM fuel cell stack | |
Duan et al. | Generalized exponential‐dispersion process model for degradation analysis under nonlinear condition | |
Zhao et al. | Fuzzy scalar quantisation based on hidden Markov model and application in fault diagnosis of wind turbine | |
WO2019216225A1 (ja) | 水素侵入挙動推定方法、水素侵入挙動推定装置及び水素侵入挙動推定プログラム | |
Wang et al. | A cumulative exposure model for remaining useful life prediction of pitot tube based on dynamic operating data | |
JP6571552B2 (ja) | 予測装置および予測方法 | |
Hou et al. | Accelerated degradation Wiener model for lithiumion battery considering individual difference |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |