CN108647458B - 一种结合制造工艺数据的机电产品退化建模方法 - Google Patents
一种结合制造工艺数据的机电产品退化建模方法 Download PDFInfo
- Publication number
- CN108647458B CN108647458B CN201810470417.2A CN201810470417A CN108647458B CN 108647458 B CN108647458 B CN 108647458B CN 201810470417 A CN201810470417 A CN 201810470417A CN 108647458 B CN108647458 B CN 108647458B
- Authority
- CN
- China
- Prior art keywords
- electronic product
- manufacturing process
- output characteristics
- degradation model
- effect coefficient
- 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
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种结合制造工艺数据的机电产品退化建模方法,所述方法首先根据机电产品的组成特点、制造工艺数据以及退化机理确定其输出特性退化模型函数形式;之后以该机电产品的制造工艺数据为基础,通过有限元仿真及近似建模方法确定退化模型中随机影响系数的分布情况;同时基于机电产品试验样本的退化数据,估计出退化模型中固定影响系数的值;最后,根据所确定退化模型函数形式以及随机影响系数与固定影响系数,给出结合制造工艺数据的机电产品退化模型。本发明解决了目前的退化建模方法由于无法对由制造工艺所决定的随机影响系数进行量化表述,而需要对其分布情况进行主观假设的问题。
Description
技术领域
本发明属于机电产品性能与可靠性分析技术领域,涉及一种机电产品退化建模方法。
背景技术
产品的可靠性是由设计、制造和管理共同决定的。对于机电产品而言,由于受到制造工艺水平的制约,其设计及制造过程所带来的参数分散性使得批次机电产品的个体之间不可避免地存在一定程度的差异性,而该差异性将直接影响机电产品的退化过程。制造工艺导致的个体差异性可视为不确定性信息。在考虑不确定性的退化建模与可靠性评估方面,Meeker和Escobar提出个体的退化轨迹模型中存在着固定影响系数以及随机影响系数。随机影响系数作为产品退化过程不确定性的具体体现,导致了产品退化轨迹的不唯一,而这种退化轨迹的区别正是由产品个体差异所引起的。
目前的研究中通常认为决定随机影响系数的产品个体差异由不可观测因素所导致,无法进行量化表述。因而,往往需要假定退化模型中的随机影响系数服从某种特定分布类型,从而根据退化数据进行分布检验及分布参数估计。但是,随着有限元仿真与近似建模技术的发展,对于许多产品来说,其个体间的差异性已经可以基于相关的设计及制造工艺数据通过仿真分析的方式具体给出。因而,结合机电产品自身的结构特点与失效机理,确定其退化模型函数形式,并以机电产品的制造工艺数据为基础,通过有限元仿真与近似建模技术实现对机电产品退化模型中随机影响系数分布情况的估计,则可以构建出能够反映批次产品退化过程的退化模型,并为进一步准确评价机电产品的可靠性水平奠定基础。
发明内容
为了解决在机电产品退化建模的相关研究中,未曾考虑制造工艺对产品退化过程的影响,以及无法对退化模型中的随机影响系数进行量化表述的问题,本发明提供了一种结合制造工艺数据的机电产品退化建模方法。
本发明的目的是通过以下技术方案实现的:
一种结合制造工艺数据的机电产品退化建模方法,包括以下步骤:
步骤一:根据机电产品的组成特点,构建其输出特性D与影响该输出特性D的n个底层性能参数pi(i∈1,…,n)的集合P={p1,p2,…,pn}之间的函数关系,所述函数关系中包含由制造工艺导致的退化模型随机影响系数;
步骤二:根据机电产品的退化机理,构建针对步骤一所述n个底层性能参数pi的n个退化物理模型;
步骤三:结合步骤一所述的函数关系与步骤二所述的退化物理模型,确定所述机电产品输出特性退化模型的函数形式,所述输出特性退化模型中包含体现所述机电产品退化过程的退化模型固定影响系数以及步骤一所述的退化模型随机影响系数;
步骤四:采用有限元仿真与机电产品输出特性近似建模方法,构建反映所述机电产品输出特性D与所述n个底层性能参数集合P以及m个内参数变化量ΔX={Δx1,Δx2,…,Δxm}之间函数关系的输出特性快速计算模型;
步骤五:在由制造工艺所导致的机电产品内参数ΔX波动范围内,应用蒙特卡洛抽样方法,随机抽样得到k个内参数组合,并将所得到的k个内参数组合分别代入步骤四所述的输出特性快速计算模型中,构建k个所述机电产品虚拟个体;
步骤六:在步骤一所述n个底层性能参数pi的变化范围内,分别等间隔抽取至少10个值,组成至少10个底层性能参数集合P的组合Pj∈(P1,P2,…,P10)(j∈1,…,10);
步骤七:将步骤六所述的Pj与步骤五所述的k个机电产品虚拟个体相结合,并代入步骤四所述的输出特性快速计算模型中,得到k组对应于所述Pj的所述输出特性Dlj∈(Dl1,Dl2,…,Dl10)的值,其中l∈(1,2,…,k);
步骤八:将步骤七中得到的k组Dlj与Pj代入步骤一所述的函数关系中,得到k个退化模型随机影响系数组合;
步骤九:根据步骤八得到的退化模型随机影响系数组合,估计步骤一所述各退化模型随机影响系数的分布情况;
步骤十:根据机电产品试验样本的退化数据,估计出步骤三所述输出特性退化模型中固定影响系数的值;
步骤十一:将步骤九所述的各退化模型随机影响系数分布情况与步骤十所述的固定影响系数值代入步骤三所述的机电产品输出特性退化模型中,完成结合制造工艺数据的机电产品退化建模。
相比于现有技术,本发明具有如下优点:
1、本发明采用有限元仿真与机电产品近似建模思想相结合的方式,以产品的具体制造工艺数据为基础,能够给出退化模型中随机影响系数的分布情况,解决了目前的退化建模方法由于无法对由制造工艺所决定的随机影响系数进行量化表述,而需要对其分布情况进行主观假设的问题。
2、本发明在对机电产品的退化建模过程中,充分考虑了制造工艺对产品退化过程的影响,并通过有限元仿真与机电产品近似建模相结合的方式将该影响量化表述为退化模型中随机影响系数的分布特性,能够在基于相同样本量试验数据的条件下更好地描述批次机电产品退化过程的统计特性,尤其对于因长寿命、高价值等原因而无法获取大量试验样本的机电产品而言,本发明为解决小子样条件下的批次产品退化建模以及可靠性评价问题提供了有效的手段。
附图说明
图1是本发明所述结合制造工艺数据的机电产品退化建模方法的流程图。
图2是本发明所述某型号电磁继电器的结构示意图。
图3是本发明所述某型号电磁继电器1000个Trls(0)的分布示意图。
图4是本发明所述某型号电磁继电器1000个参数c的分布示意图。
具体实施方式
下面结合附图对本发明的技术方案作进一步的说明,但并不局限于此,凡是对本发明技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,均应涵盖在本发明的保护范围中。
具体实施方式一:本实施方式提供了一种结合制造工艺数据的机电产品退化建模方法,所述方法首先根据机电产品的组成特点、制造工艺数据以及退化机理确定其输出特性退化模型函数形式,该模型中包括固定影响系数以及由制造工艺决定的随机影响系数;之后以该机电产品的制造工艺数据为基础,通过有限元仿真及近似建模方法确定退化模型中随机影响系数的分布情况;同时基于机电产品试验样本的退化数据,估计出退化模型中固定影响系数的值;最后,根据所确定退化模型函数形式以及随机影响系数与固定影响系数,给出结合制造工艺数据的机电产品退化模型。如图1所示,具体步骤如下:
步骤一:根据机电产品的组成特点,构建其输出特性D与影响该输出特性D的n个底层性能参数pi(i∈1,…,n)的集合P={p1,p2,…,pn}之间的函数关系,所述函数关系中包含由制造工艺导致的退化模型随机影响系数。本步骤中,所述的输出特性D与n个底层性能参数集合P之间的函数关系如公式(1)所示:
D=F(P|Θ1) (1)
式中,Θ1为由制造工艺导致的退化模型随机影响系数集合,表征批次产品间的个体差异性。
步骤二:根据机电产品的退化机理,构建针对步骤一所述n个底层性能参数pi的n个退化物理模型。本步骤中,所述n个底层性能参数pi的n个退化物理模型可表示为如公式(2)所示的关于时间t及应力水平S的函数关系:
式中,pi(t,S)表示产品在应力水平S下,在t时刻,其底层性能参数pi的值;为对应于所述退化物理模型的未知系数集合。
步骤一和步骤二中所述函数关系的常用数学表达式如表1所示:
表1
步骤三:结合步骤一所述的函数关系与步骤二所述的退化物理模型,确定所述机电产品输出特性退化模型的函数形式,所述输出特性退化模型中包含体现所述机电产品退化过程的退化模型固定影响系数以及步骤一所述的退化模型随机影响系数。本步骤中,所述输出特性退化模型的函数形式可通过公式(3)表示:
式中,Θ2统一表征P={p1,p2,…,pn}中各底层性能参数所对应的Θ2为所述的退化模型固定影响系数集合,Θ1为步骤一所述的退化模型随机影响系数集合。
步骤四:采用有限元仿真与机电产品输出特性近似建模方法,构建反映所述机电产品输出特性D与所述n个底层性能参数集合P以及m个内参数变化量ΔX={Δx1,Δx2,…,Δxm}之间函数关系的输出特性快速计算模型;
步骤五:在由制造工艺所导致的机电产品内参数ΔX波动范围内,应用蒙特卡洛抽样方法,随机抽样得到k个内参数组合,并将所得到的k个内参数组合分别代入步骤四所述的输出特性快速计算模型中,构建k个所述机电产品虚拟个体;
步骤六:在步骤一所述n个底层性能参数pi的变化范围内,分别等间隔抽取至少10个值,组成至少10个底层性能参数集合P的组合Pj∈(P1,P2,…,P10)(j∈1,…,10);
步骤七:将步骤六所述的Pj与步骤五所述的k个机电产品虚拟个体相结合,并代入步骤四所述的输出特性快速计算模型中,得到k组对应于所述Pj的所述输出特性Dlj∈(Dl1,Dl2,…,Dl10)的值,其中l∈(1,2,…,k);
步骤八:将步骤七中得到的k组Dlj与Pj代入步骤一所述的函数关系中,得到k个退化模型随机影响系数组合;
步骤九:根据步骤八得到的退化模型随机影响系数组合,估计步骤一所述各退化模型随机影响系数的分布情况;
步骤十:根据机电产品试验样本的退化数据,估计出步骤三所述输出特性退化模型中固定影响系数的值;
步骤十一:将步骤九所述的各退化模型随机影响系数分布情况与步骤十所述的固定影响系数值代入步骤三所述的机电产品输出特性退化模型中,完成结合制造工艺数据的机电产品退化建模。
具体实施方式二:下面结合图1至图4说明本实施方式,本实施方式对具体实施方式一的进一步说明。本实施方式以某型号电磁继电器为对象进行结合制造工艺数据的退化建模,具体步骤如下:
步骤一:根据公式(1)及所述电磁继电器的结构特点,构建其释放时间Trls(所述输出特性)关于簧片初力Fini(所述底层性能参数)的函数关系,如公式(4)所示:
Trls=F(Fini|Θ1)=kFini+c (4)
式中,Θ1={k,c}为所述退化模型随机影响系数集合。
步骤二:根据所述电磁继电器簧片的退化机理,构建簧片初力的退化物理模型,如公式(5)所示:
式中,Fini(t,T)为在温度应力T下贮存至t时刻的簧片初力值;Fini(0)为0时刻的簧片初力值;(a,b)为模型系数,对应于公式(2)中的C通常取为20。
步骤三:结合步骤一所述的函数关系与步骤二所述的退化模型,确定所述电磁继电器释放时间退化模型的函数形式,如公式(6)所示:
Trls(t,T)=kFini(t,T)+c
=kFini(0)(aθ+b)+c (6)
=[Trls(0)-c](aθ+b)+c
式中,Trls(t,T)为在温度应力T下贮存至t时刻的释放时间值;Trls(0)为0时刻的释放时间值;θ=(1.8T+491.67)(logt+C)。
步骤四:采用有限元仿真与机电产品输出特性近似建模方法,构建可反映所述电磁继电器释放时间与所述所述电磁继电器簧片初力值以及内参数变化量ΔX={Δx1,Δx2,…,Δxm}之间函数关系的释放时间快速计算模型。
步骤五:在由制造工艺所导致的机电产品内参数ΔX波动范围内,应用蒙特卡洛抽样方法,随机抽样得到1000个内参数组合,并将所得到的1000个内参数组合分别代入步骤四所述的释放时间快速计算模型中,构建1000个所述电磁继电器虚拟个体。
步骤六:在所述簧片初力的变化范围(0.05N-0.35N)内,等间隔抽取10个值,组成10个所述簧片初力值。
步骤七:将步骤六得到的10个簧片初力值与步骤五所述的1000个所述电磁继电器虚拟个体相结合,并代入步骤四所述的释放时间快速计算模型中,得到1000组对应于所述簧片初力值的所述释放时间的值。
步骤八:将步骤七中得到的1000组释放时间值与步骤六所述的10个簧片初力值共同代入公式(4)所述的函数关系中,得到如图3所示的1000个Trls(0)与如图4所示的1000个c。
步骤九:根据步骤八得到的1000组Trls(0)与c,得到Trls(0)的正态分布均值与标准差分别为1706.2277与141.5262,c的正态分布均值与标准差分别为3237.8428与299.6957。
步骤十:根据10个所述电磁继电器试验样本的退化数据,估计出公式(6)所述退化模型中a与b的值分别为-5.7060×10-5与2.0904。
步骤十一:将步骤九所述的Trls(0)与c分布情况与步骤十所述的a与b的值代入公式(6)所述的电磁继电器释放时间退化模型中,完成结合制造工艺数据的退化建模。
Claims (4)
1.一种结合制造工艺数据的机电产品退化建模方法,其特征在于所述方法步骤如下:
步骤一:根据机电产品的组成特点,构建其输出特性D与影响该输出特性D的n个底层性能参数pi的集合P={p1,p2,…,pn}之间的函数关系,所述函数关系中包含由制造工艺导致的退化模型随机影响系数,i∈1,…,n;
步骤二:根据机电产品的退化机理,构建针对步骤一所述n个底层性能参数pi的n个退化物理模型;
步骤三:结合步骤一所述函数关系与步骤二所述退化物理模型,确定所述机电产品输出特性退化模型的函数形式,所述输出特性退化模型中包含体现所述机电产品退化过程的退化模型固定影响系数以及步骤一所述退化模型随机影响系数;
步骤四:采用有限元仿真与机电产品输出特性近似建模方法,构建反映所述机电产品输出特性D与所述n个底层性能参数集合P以及m个内参数变化量ΔX={Δx1,Δx2,…,Δxm}之间函数关系的输出特性快速计算模型;
步骤五:在由制造工艺所导致的机电产品内参数变化量ΔX波动范围内,应用蒙特卡洛抽样方法,随机抽样得到k个内参数组合,并将所得到的k个内参数组合分别代入步骤四所述输出特性快速计算模型中,构建k个所述机电产品虚拟个体;
步骤六:在步骤一所述n个底层性能参数pi的变化范围内,分别等间隔抽取至少10个值,组成至少10个底层性能参数集合P的组合Pj∈(P1,P2,…,P10),j∈1,…,10;
步骤七:将步骤六所述Pj与步骤五所述k个机电产品虚拟个体相结合,并代入步骤四所述输出特性快速计算模型中,得到k组对应于所述Pj的所述输出特性Dlj∈(Dl1,Dl2,…,Dl10)的值,其中l∈(1,2,…,k);
步骤八:将步骤七中得到的k组Dlj与Pj代入步骤一所述函数关系中,得到k个退化模型随机影响系数组合;
步骤九:根据步骤八得到的退化模型随机影响系数组合,估计步骤一所述退化模型随机影响系数的分布情况;
步骤十:根据机电产品试验样本的退化数据,估计出步骤三所述输出特性退化模型中固定影响系数的值;
步骤十一:将步骤九所述退化模型随机影响系数分布情况与步骤十所述固定影响系数值代入步骤三所述机电产品输出特性退化模型中,完成结合制造工艺数据的机电产品退化建模。
2.根据权利要求1所述的结合制造工艺数据的机电产品退化建模方法,其特征在于所述步骤一中,输出特性D与n个底层性能参数集合P之间的函数关系如公式(1)所示:
D=F(P|Θ1) (1)
式中,Θ1为退化模型随机影响系数集合。
3.根据权利要求1所述的结合制造工艺数据的机电产品退化建模方法,其特征在于所述步骤二中,n个底层性能参数pi的n个退化物理模型表示为如公式(2)所示的关于时间t及应力水平S的函数关系:
式中,pi(t,S)表示产品在应力水平S下,在t时刻,其底层性能参数pi的值;为对应于所述退化物理模型的未知系数集合。
4.根据权利要求1所述的结合制造工艺数据的机电产品退化建模方法,其特征在于所述步骤三中,输出特性退化模型的函数形式通过公式(3)表示:
式中,Θ2统一表征P={p1,p2,…,pn}中各底层性能参数所对应的 为对应于所述退化物理模型的未知系数集合,Θ1为步骤一所述退化模型随机影响系数集合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810470417.2A CN108647458B (zh) | 2018-05-16 | 2018-05-16 | 一种结合制造工艺数据的机电产品退化建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810470417.2A CN108647458B (zh) | 2018-05-16 | 2018-05-16 | 一种结合制造工艺数据的机电产品退化建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108647458A CN108647458A (zh) | 2018-10-12 |
CN108647458B true CN108647458B (zh) | 2019-05-03 |
Family
ID=63756156
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810470417.2A Active CN108647458B (zh) | 2018-05-16 | 2018-05-16 | 一种结合制造工艺数据的机电产品退化建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108647458B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111506998B (zh) * | 2020-04-15 | 2022-02-08 | 哈尔滨工业大学 | 一种构建机电产品制造过程参数漂移故障特征样本库的方法 |
CN113255128A (zh) * | 2021-05-21 | 2021-08-13 | 哈尔滨工业大学 | 一种考虑产品分散性的退化机理辨识方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104573881A (zh) * | 2015-02-10 | 2015-04-29 | 广东石油化工学院 | 一种基于退化数据建模的服役设备剩余寿命自适应预测方法 |
CN106021685A (zh) * | 2016-05-16 | 2016-10-12 | 北京航空航天大学 | 一种考虑测量误差的退化可靠性分析方法 |
CN107515965A (zh) * | 2017-07-27 | 2017-12-26 | 北京航空航天大学 | 一种基于不确定过程的加速退化建模评估方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10289769B2 (en) * | 2013-08-26 | 2019-05-14 | The University Of Akron | Method for identifying stochastic information of heterogeneous materials |
-
2018
- 2018-05-16 CN CN201810470417.2A patent/CN108647458B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104573881A (zh) * | 2015-02-10 | 2015-04-29 | 广东石油化工学院 | 一种基于退化数据建模的服役设备剩余寿命自适应预测方法 |
CN106021685A (zh) * | 2016-05-16 | 2016-10-12 | 北京航空航天大学 | 一种考虑测量误差的退化可靠性分析方法 |
CN107515965A (zh) * | 2017-07-27 | 2017-12-26 | 北京航空航天大学 | 一种基于不确定过程的加速退化建模评估方法 |
Non-Patent Citations (3)
Title |
---|
仿真基混合效应模型加速退化试验;汪亚顺等;《机械工程学报》;20091231;第45卷(第12期);全文 |
基于性能退化和 Monte-Carlo 仿真的系统;周月阁等;《仪器仪表学报》;20140531;第35卷(第5期);全文 |
航天继电器贮存过程吸合时间退化机理研究;叶雪荣等;《电工技术学报》;20170630;第32卷(第11期);全文 |
Also Published As
Publication number | Publication date |
---|---|
CN108647458A (zh) | 2018-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107730127B (zh) | 基于输出特性初始分布的继电器贮存退化数据预测方法 | |
CN108647458B (zh) | 一种结合制造工艺数据的机电产品退化建模方法 | |
CN109657937A (zh) | 一种基于退化数据的产品可靠性评估与寿命预测方法 | |
CN114818348B (zh) | 考虑多应力耦合作用对产品退化影响的可靠性评估方法 | |
CN112883550A (zh) | 一种考虑多重不确定性的退化设备剩余寿命预测方法 | |
CN109598052B (zh) | 基于相关系数分析的智能电表寿命周期预测方法及装置 | |
CN110083910B (zh) | 一种基于nsga-ⅱ的混沌时间序列预测样本获取方法 | |
CN107479000A (zh) | 一种基于Box‑Cox变换与蒙特卡罗仿真的锂离子动力电池RUL预测方法 | |
CN106874634A (zh) | 基于逆高斯退化模型的剩余寿命贝叶斯预测方法 | |
CN107885930A (zh) | 继电器弹性金属材料加速贮存退化失效机理变化判别方法 | |
CN103324787A (zh) | 批产卫星热设计验证方法 | |
CN108710745B (zh) | 一种结合制造工艺数据的机电产品可靠性评估方法 | |
CN109191408A (zh) | 快速循环地面气象融合方法、装置及服务器 | |
CN109657260B (zh) | 一种考虑失效相关性的涡轮转子系统可靠性分配方法 | |
CN117077546B (zh) | 基于数据驱动的电力系统负荷建模方法及其系统 | |
CN107704691A (zh) | 一种加速应力可靠性鉴定试验优选方案设计方法 | |
CN107958129B (zh) | 模拟氧化锌压敏电阻微观电流分布的算法 | |
Li et al. | Storage life prediction method of the aerospace electromagnetic relays based on physics of failure and data-driven fusion | |
CN116307773A (zh) | 一种变电站二次设备可靠性估计方法 | |
CN111625995B (zh) | 一种集成遗忘机制和双超限学习机的在线时空建模方法 | |
CN113919162A (zh) | 基于仿真与多源实测数据融合的电压暂降风险预警方法 | |
CN111474905B (zh) | 一种机电产品制造过程参数漂移故障诊断方法 | |
CN109390946B (zh) | 一种基于多参数规划理论的最优概率潮流快速计算方法 | |
CN111506998B (zh) | 一种构建机电产品制造过程参数漂移故障特征样本库的方法 | |
CN104915890B (zh) | 一种适用于大区互联电网的负荷模型参数拟合方法 |
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 |