CN108984925A - 面向电子产品可靠性综合仿真分析的并行计算方法 - Google Patents
面向电子产品可靠性综合仿真分析的并行计算方法 Download PDFInfo
- Publication number
- CN108984925A CN108984925A CN201810822767.0A CN201810822767A CN108984925A CN 108984925 A CN108984925 A CN 108984925A CN 201810822767 A CN201810822767 A CN 201810822767A CN 108984925 A CN108984925 A CN 108984925A
- Authority
- CN
- China
- Prior art keywords
- analysis
- task
- sequence
- electronic product
- vibration
- 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.)
- Granted
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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
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)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种面向电子产品可靠性综合仿真分析的并行计算方法。步骤如下:1确定电子产品分散性,设定电子产品尺寸及影响电子产品寿命的热和振动的参数分散性。2构建载荷剖面,根据产品使用信息,建立电子产品进行热分析和振动分析的载荷剖面。3进行可靠性评估分析,对分散性参数进行随机抽样,并驱动相应的热分析和振动分析,进行可靠性评估分析。4建立计算任务执行序列,将所要执行的全部计算任务的序列随机重排,形成新的计算任务序列。5计算任务分配,按照新的计算任务序列将所有计算任务分配到计算节点中,并生成相应的数据文件F热,i,F振,i。6进行失效时间拟合,任务节点收集可靠性评估所得的失效时间并进行失效时间拟合。
Description
所属技术领域
本发明提供了一种面向电子产品可靠性综合仿真分析的并行计算方法,特别涉及到对电子产品进行热分析和振动分析后的可靠性综合仿真分析过程,同等条件下利用本发明所提供的方法可显著提高计算效率,减小计算时间。本发明属于可靠性与系统工程领域。
背景技术
在对电子产品进行可靠性综合仿真分析时,需要对电子产品进行产品建模,构建寿命剖面,对电子产品进行热分析、振动分析和可靠性评估,最终以报告形式呈现。可靠性综合仿真分析在短时间内需要处理大量数据,开发一种高效的并行计算方法就尤为重要。
并行计算技术是指把一台计算机或者服务器抽象成一个独立的计算节点,多个计算节点共同组成一个计算集群,节点与节点之间通过网络来通信。一次庞大的计算作业可以被拆分成一个个相对独立的计算任务,这些计算任务被均匀地分发到各个计算节点上,每个节点的计算任务完成后,管理节点会对任务重新调度分配,以进行结果的汇总。
本发明研究了面向电子产品可靠性综合仿真分析的计算过程,发明了一种新型高效的并行计算方法,即面向电子产品可靠性综合仿真分析的并行计算方法,可减小电子产品可靠性综合仿真分析的时间,提高计算效率。
发明内容
本发明的目的是为电子产品可靠性综合仿真分析提供一种并行计算方法,旨在减小对电子产品进行可靠性综合仿真分析时的任务时间,提高计算效率。
本发明的目的是提出了一种面向电子产品可靠性综合仿真分析的并行计算方法,主要包含以下步骤:
步骤一:确定电子产品分散性
根据电子产品设计情况收集相关数据信息,包括电子产品的设计信息、材料信息、封装信息、使用条件信息等相关数据,设定电子产品由于工艺或标准导致的分散性,包括材料参数及尺寸参数等,然后设定影响电子产品进行热分析和振动分析的参数分散性。
步骤二:构建载荷剖面
根据电子产品的使用信息,明确电子产品在全寿命周期中经历的各任务阶段时间及对应的环境条件,建立起电子产品进行热分析和振动分析的载荷剖面。
步骤三:进行可靠性评估分析
对电子产品在尺寸等具有分散性的属性值及影响电子产品寿命的热和振动等因素的参数进行随机抽样,并驱动在温度载荷剖面和振动载荷剖面下的热分析和振动分析,在分析结束后按照分析结果对电子产品进行不同失效物理模型下的可靠性评估分析,得到不同抽样结果及载荷剖面下的电子产品的失效时间t。
步骤四:建立计算任务执行序列
将对电子产品进行可靠性仿真分析的计算划分为热分析计算N热,i、振动分析计算N振,i及可靠性评估分析计算N评,i三部分,其中i=1,2…n。对于每一部分都需按照相应的算法对其进行计算。按照随机任务分配的方式,将所要执行的全部计算任务进行重新排序,形成新的任务计算序列。首先将所要执行的计算任务序列Q={N3n},对应生成随机数序列{R3n}。然后将序列{R3n}从小到大排序,排序后的序列为{S3n}。按照随机数序列{S3n}得到新的计算任务序列Q’={L3n}。
步骤五:计算任务分配
按照新的计算任务序列将所有计算任务分配到p个执行任务的计算节点中,计算节点执行完热分析和振动分析后生成相应的数据文件F热,i,F振,i,其中i=1,2…n。当计算节点所分得的可靠性评估任务N评,i无对应次数热分析和振动分析的数据文件F热,i,F振,i时,该计算节点暂停接受计算任务的分配,并执行尚未执行的F热,i,F振,i,且该计算任务从计算序列Q’={L3n}里消除,然后再执行所分得的可靠性评估任务N评,i,待此次可靠性评估任务执行完毕后再接受任务序列里任务的分配。
步骤六:进行失效时间拟合
当计算任务序列Q’={L3n}中的计算任务均执行完毕后,任务节点收集可靠性评估所得的失效时间并进行失效时间拟合,得出最终的可靠性分析结果,完成整个电子产品可靠性分析过程。
附图说明
图1为本发明中所述方法的整体架构框图
图2计算任务分配机制
具体实施方式
为使本发明的技术方案、特征及优点得到更清楚的了解,以下结合附图,作详细说明。
本发明给出了一种面向电子产品可靠性综合仿真分析的并行计算方法,可减小涉及到对电子产品进行热分析和振动分析后的可靠性综合仿真分析过程的计算时间。本发明的整体架构,见图1所示,下面以实例进一步说明本发明的实质内容,但本发明的内容并不限于此。
步骤一:确定电子产品分散性
根据电子产品设计情况收集相关数据信息,包括电子产品的设计信息、材料信息、封装信息、使用条件信息等相关数据,设定电子产品由于工艺或标准导致的分散性,包括材料参数及尺寸参数等,然后设定影响电子产品进行热分析和振动分析的参数分散性。
例1:电子产品由于工艺或标准导致某些尺寸参数的分布类型和分布参数不同,此外,电子产品上的某些参数的分散性会影响其进行热分析和振动分析。
电子产品具有分散性的属性参数及影响电子产品寿命的热和振动等因素的参数如表1所示:
表1具有分散性的参数
序号 | 子单元 | 属性 | 分布类型 |
1 | 焊点 | 高度 | 泊松分布 |
2 | 器件 | 长度/宽度 | 均匀分布 |
3 | 封装 | 引脚跨度 | 三角分布 |
4 | 基板 | 断裂韧性 | 威布尔分布 |
步骤二:构建载荷剖面
根据电子产品的使用信息,明确电子产品在全寿命周期中经历的各任务阶段时间及对应的环境条件,建立起电子产品进行热分析和振动分析的载荷剖面。
步骤三:进行可靠性评估分析
对电子产品在尺寸等具有分散性的属性值及影响电子产品寿命的热和振动等因素的参数进行随机抽样,并驱动在温度载荷剖面和振动载荷剖面下的热分析和振动分析,在分析结束后按照分析结果对电子产品进行不同失效物理模型下的可靠性评估分析,得到不同抽样结果及载荷剖面下的电子产品的失效时间t。
步骤四:建立计算任务执行序列
将对电子产品进行可靠性仿真分析的计算划分为热分析计算N热,i、振动分析计算N振,i及可靠性评估分析计算N评,i三部分,其中i=1,2…n。对于每一部分都需按照相应的算法对其进行计算。按照随机任务分配的方式,将所要执行的全部计算任务进行重新排序,形成新的任务计算序列。首先将所要执行的计算任务序列Q={N3n},对应生成随机数序列{R3n}。然后将序列{R3n}从小到大排序,排序后的序列为{S3n}。按照随机数序列{S3n}得到新的计算任务序列Q’={L3n}。
例2:随机任务分配方式,其关键在于将原来的计算任务序列Q={N3n}进行随机重排,产生新的计算任务序列Q’={L3n}。利用所提出的计算任务序列重新建立的方法能够实现计算任务对于各计算节点的随机分配,从而避免各计算节点同时计算的任务占用的峰值内存总和受到进程数占用峰值内存的最高值所限制。
步骤五:计算任务分配
按照新的计算任务序列将所有计算任务分配到p个执行任务的计算节点中,计算节点执行完热分析和振动分析后生成相应的数据文件F热,i,F振,i,其中i=1,2…n。当计算节点所分得的可靠性评估任务N评,i无对应次数热分析和振动分析的数据文件F热,i,F振,i时,该计算节点暂停任务分配的响应,并执行尚未执行的F热,i,F振,i,且该计算任务从计算序列Q’={L3n}里消除,然后再执行所分得的可靠性评估任务N评,i,待此次可靠性评估任务执行完毕后再接受任务序列里任务的分配。
例3:如图2所示,按照新的计算任务序列Q’={L3n}将所有计算任务分配到p个执行任务的计算节点中,当计算节点所分得的可靠性评估任务N评,i无对应次数热分析和振动分析的数据文件F热,i,F振,i时,该计算节点暂停任务分配的响应,并执行尚未执行的F热,i,F振,i,且该计算任务从计算序列Q’={L3n}里消除,然后再执行所分得的可靠性评估任务N评,i,待计算序列中的所有计算任务均被执行完毕时计算结束。
步骤六:进行失效时间拟合
当计算任务序列S={Ln}中的计算任务均执行完毕后,任务节点收集可靠性评估所得的失效时间并进行失效时间拟合,得出最终的可靠性分析结果,完成整个电子产品可靠性分析过程。
Claims (4)
1.面向电子产品可靠性综合仿真分析的并行计算方法,它包含以下步骤:
第一步:确定电子产品分散性:根据电子产品设计情况收集相关数据信息,包括电子产品的设计信息、材料信息、封装信息、使用条件信息等相关数据,设定电子产品由于工艺或标准导致的分散性,包括材料参数及尺寸参数等,然后设定影响电子产品进行热分析和振动分析的参数分散性。
第二步:构建载荷剖面:根据电子产品的使用信息,明确电子产品在全寿命周期中经历的各任务阶段时间及对应的环境条件,建立起电子产品进行热分析和振动分析的载荷剖面。
第三步:进行基于蒙特卡洛方法的可靠性评估:对影响电子产品寿命的热和振动等因素的参数进行随机抽样,并驱动在温度剖面和载荷剖面下的热分析和振动分析,在分析结束后按照分析结果对电子产品进行不同失效物理模型下的可靠性评估分析,得到电子产品的失效时间t。
第四步:建立计算任务执行序列:将对电子产品进行可靠性仿真分析的计算划分为热分析计算N热,i、振动分析计算N振,i及可靠性评估分析计算N评,i三部分,其中i=1,2…n。对于每一部分都需按照相应的算法对其进行计算。按照随机任务分配的方式,将所要执行的全部计算任务进行重新排序,形成新的任务计算顺序。首先将所要执行的计算任务序列Q={N3n},对应生成随机数序列{R3n}。然后将序列{R3n}从小到大排序,排序后的序列为{S3n}。按照随机数序列{S3n}得到新的计算任务序列Q’={L3n}。
第五步:计算任务分配:按照新的计算任务序列将所有计算任务分配到p个执行任务的计算节点中,计算节点执行完热分析和振动分析后生成相应的数据文件F热,i,F振,i,其中i=1,2…n。当计算节点所分得的可靠性评估任务N评,i无对应次数热分析和振动分析的数据文件F热,i,F振,i时,该计算节点暂停任务分配的响应,并执行尚未执行的F热,i,F振,i,且该计算任务从计算序列S={Ln}里消除,然后再执行所分得的可靠性评估任务N评,i,待此次可靠性评估任务执行完毕后再接受任务序列里任务的分配。
第六步:进行失效时间拟合:当计算任务序列S={Ln}中的计算任务均执行完毕后,任务节点收集可靠性评估所得的失效时间并进行失效时间拟合,得出最终的可靠性分析结果,完成整个电子产品可靠性分析过程。
2.根据权利要求1所述的面向电子产品可靠性综合仿真分析的并行计算方法,其特征在于:在第二步中所述的“构建载荷剖面”中,根据电子产品的使用信息,明确电子产品在全寿命周期中经历的各任务阶段时间及对应的环境条件,建立起电子产品进行热分析和振动分析的载荷剖面。
3.根据权利要求1所述的面向电子产品可靠性综合仿真分析的并行计算方法,其特征在于:在第四步中所述的“建立计算任务执行序列”中,按照随机任务分配的方式,将所要执行的全部计算任务进行重新排序,形成新的任务计算顺序。首先将所要执行的计算任务序列Q={N3n},对应生成随机数序列{R3n}。然后将序列{R3n}从小到大排序,排序后的序列为{S3n}。按照随机数序列{S3n}得到新的计算任务序列Q’={L3n}。
4.根据权利要求1所述的面向电子产品可靠性综合仿真分析的并行计算方法,其特征在于:在第五步中所述的“计算任务分配”中,按照新的计算任务序列将所有计算任务分配到p个执行任务的计算节点中,计算节点执行完热分析和振动分析后生成相应的数据文件F热,i,F振,i,其中i=1,2…n。当计算节点所分得的可靠性评估任务N评,i无对应次数热分析和振动分析的数据文件F热,i,F振,i时,该计算节点暂停任务分配的响应,并执行尚未执行的F热,i,F振,i,且该计算任务从计算序列S={Ln}里消除,然后再执行所分得的可靠性评估任务N评,i,待此次可靠性评估任务执行完毕后再接受任务序列里任务的分配。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810822767.0A CN108984925B (zh) | 2018-07-25 | 2018-07-25 | 面向电子产品可靠性综合仿真分析的并行计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810822767.0A CN108984925B (zh) | 2018-07-25 | 2018-07-25 | 面向电子产品可靠性综合仿真分析的并行计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108984925A true CN108984925A (zh) | 2018-12-11 |
CN108984925B CN108984925B (zh) | 2023-04-18 |
Family
ID=64550281
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810822767.0A Active CN108984925B (zh) | 2018-07-25 | 2018-07-25 | 面向电子产品可靠性综合仿真分析的并行计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108984925B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109783970A (zh) * | 2019-01-29 | 2019-05-21 | 北京航空航天大学 | 面向电子产品可靠性仿真分析的高效简易热分析方法 |
CN110222428A (zh) * | 2019-06-10 | 2019-09-10 | 哈尔滨工程大学 | 一种面向系统级封装sip器件的可靠性分析系统及方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005088483A1 (ja) * | 2004-03-10 | 2005-09-22 | Fujitsu Limited | 電子パッケージ信頼性予測装置および電子パッケージ信頼性予測プログラム |
CN103559418A (zh) * | 2013-11-20 | 2014-02-05 | 北京航空航天大学 | 一种基于故障物理的电子产品任务可靠度计算方法 |
CN106529090A (zh) * | 2016-12-12 | 2017-03-22 | 中国航天标准化研究所 | 一种航天电子类产品可靠性评估方法 |
-
2018
- 2018-07-25 CN CN201810822767.0A patent/CN108984925B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005088483A1 (ja) * | 2004-03-10 | 2005-09-22 | Fujitsu Limited | 電子パッケージ信頼性予測装置および電子パッケージ信頼性予測プログラム |
CN103559418A (zh) * | 2013-11-20 | 2014-02-05 | 北京航空航天大学 | 一种基于故障物理的电子产品任务可靠度计算方法 |
CN106529090A (zh) * | 2016-12-12 | 2017-03-22 | 中国航天标准化研究所 | 一种航天电子类产品可靠性评估方法 |
Non-Patent Citations (3)
Title |
---|
WUYANG PAN 等: "《Comparison of stochastic response surface method and Monte Carlo method for uncertainty analysis of electronics prognostics》", 《2015 IEEE CONFERENCE ON PROGNOSTICS AND HEALTH MANAGEMENT (PHM)》 * |
万博等: "航空电子产品可靠性仿真预计数据处理方法研究", 《电子产品可靠性与环境试验》 * |
骆明珠 等: "《基于PoF模型的电子产品可靠性参数计算方法》", 《系统工程与电子技术》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109783970A (zh) * | 2019-01-29 | 2019-05-21 | 北京航空航天大学 | 面向电子产品可靠性仿真分析的高效简易热分析方法 |
CN109783970B (zh) * | 2019-01-29 | 2021-01-15 | 北京航空航天大学 | 一种面向电子产品可靠性仿真分析的热分析方法 |
CN110222428A (zh) * | 2019-06-10 | 2019-09-10 | 哈尔滨工程大学 | 一种面向系统级封装sip器件的可靠性分析系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108984925B (zh) | 2023-04-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109559039A (zh) | 一种基于数字孪生技术的虚拟工厂业务协同系统 | |
Claes et al. | Merging event logs for process mining: A rule based merging method and rule suggestion algorithm | |
Israilevich et al. | The choice of an input-output table embedded in regional econometric input-output models | |
CN107358352A (zh) | 基于蒙特卡洛模拟的电力系统可靠性评估系统及方法 | |
CN109086546A (zh) | 信号链路信号质量评估方法、装置、设备及可读存储介质 | |
CN108984925A (zh) | 面向电子产品可靠性综合仿真分析的并行计算方法 | |
CN104182623A (zh) | 一种基于当量变化率计算的热工过程数据检测方法 | |
CN105677556A (zh) | 一种调度自动化主站系统软件测试方法及其平台 | |
CN109933488A (zh) | 芯片温度计算方法及芯片温度计算装置 | |
CN115328068A (zh) | 应用于工业生产的数字孪生系统 | |
CN110471857A (zh) | 人工智能模型性能容量的自动测试方法及装置 | |
CN105631612A (zh) | 基于大数据的公职人员个人绩效及能力评估系统及方法 | |
CN103699756B (zh) | 松耦合模拟系统及方法 | |
Srinivas et al. | Multi-site downscaling of maximum and minimum daily temperature using support vector machine | |
CN108369411A (zh) | 显示方法、显示程序以及显示控制装置 | |
CN106897505A (zh) | 一种考虑时‑空相关性的结构监测数据异常识别方法 | |
Lu et al. | Modeling and simulation of aircraft assembly line based on quest | |
CN102541811B (zh) | 基于按需计算的分析因子的数据分析装置及方法 | |
CN108564232A (zh) | 一种基于遗传算法的飞机部件装配车间总体布局优化方法 | |
CN111931323B (zh) | 存储器、加氢裂化设备故障预测方法、装置和设备 | |
CN104850923A (zh) | 半导体生产仿真系统 | |
US8180598B2 (en) | System and method for performing thermal analysis on a building through universal meshing | |
Roy et al. | Productivity Improvement of a fan manufacturing company by using DMAIC approach: A Six-Sigma practice | |
CN109388830A (zh) | 一种电路板寿命预测方法 | |
CN106934729A (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 |