CN111695315B - 一种辐射环境下cpu失效仿真方法 - Google Patents
一种辐射环境下cpu失效仿真方法 Download PDFInfo
- Publication number
- CN111695315B CN111695315B CN202010454765.8A CN202010454765A CN111695315B CN 111695315 B CN111695315 B CN 111695315B CN 202010454765 A CN202010454765 A CN 202010454765A CN 111695315 B CN111695315 B CN 111695315B
- Authority
- CN
- China
- Prior art keywords
- output result
- breakdown voltage
- resistance
- reverse breakdown
- failure
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/32—Circuit design at the digital level
- G06F30/33—Design verification, e.g. functional simulation or model checking
- G06F30/3308—Design verification, e.g. functional simulation or model checking using simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2115/00—Details relating to the type of the circuit
- G06F2115/10—Processors
-
- 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/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- 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
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D10/00—Energy efficient computing, e.g. low power processors, power management or thermal management
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Tests Of Electronic Circuits (AREA)
Abstract
该发明公开了一种辐射环境下CPU失效仿真方法,属于可靠性实验、建模仿真与计算技术领域。解决了迫于试验周期和经济成本等因素的制约,整个评估体系将面临电源系统组件未能完全失效、试验对象数量有限、退化数据样本较少等一系列问题。本发明通过simulink仿真出驱动SIP电路中CPU板在辐射环境下的失效模型,通过误码率来度量其失效的过程,为解决在辐射环境这种苛刻条件下的实验时条件不足,数据缺乏等问题提供了更多的可行性和操作性,也为其他电路源件在特殊环境下的工作和失效机理的模型仿真研究提供了一些思路。
Description
技术领域
本发明属于可靠性实验、建模仿真与计算技术领域,特别涉及驱动SiP辐射可靠性仿真建模中CPU板在辐照强度下失效模型的仿真。
背景技术
电子产品是实现设备自动化、智能化的关键纽带,电源系统作为保障电子产品运行的直接动力来源,一旦退化失效,将导致整机瘫痪,造成巨大的经济损失。目前,国内外对电源系统的退化失效模式分析开展了广泛的研究,但是在某些特殊环境,特别是处于核辐射环境下的退化失效分析以及可靠性评估十分匮乏。因此,本文围绕多路输出电源系统开展了核辐射环境下的退化分析及可靠性评估工作,致力于构建一套系统硬件开发、试验方案设计、摸底试验探索以及理论分析的完整可靠性评价体系。基于前期对自主设计的电源系统在辐照环境下的3组摸底试验发现,该电源系统中有源器件对核辐射应力十分敏感,是导致有源器件退化失效的主要原因。然而,迫于试验周期和经济成本等因素的制约,整个评估体系将面临电源系统组件未能完全失效、试验对象数量有限、退化数据样本较少等一系列问题。
典型电路板是电子设备系统正常运行的必备条件,包含CPU板等在内的控制器件作为电子电路的核心,其工作的可靠性将会决定整个线性电源的可靠性性能。
发明内容
根据现有技术存在的问题,本发明提供了一种利用误码率这个数据来作为CPU板在辐射环境下失效的过程,解决了迫于试验周期和经济成本等因素的制约,整个评估体系将面临电源系统组件未能完全失效、试验对象数量有限、退化数据样本较少等一系列问题。
本发明技术方案为一种辐射环境下CPU失效仿真方法,该方法包括:
步骤1:采用如下函数作为误码率与辐照强度的关系;
wuma=0.00129173x3-0.465764727x2+55.9605636x-2239.4157835
其中,wuma为误码率,x为辐照强度;
步骤2:分别设置CPU中正向导通电阻阈值t1和反向击穿电压阈值t2;
步骤3:针对正向导通电阻和反向击穿电压各设置一个变化区间,正向导通电阻的变化区间为h1,反向击穿电压的变化区间为h2;
步骤4:变化区间随着误码率的变化会发生移动,采用如下公式作为变化区间与误码率之间的关系;
Zmin=t1-h1×wuma Zmax=h1×(1-wuma)+t1
Fmin=t2-h2×wuma Fmax=h2×(1-wuma)+t2
Zmin~Zmax为正向导通电阻的变化区间,Fmin~Fmax为反向击穿电压的变化区间;
步骤5:采用连续均匀分布的随机方法在正向导通电阻的变化区间内随机产生一个正向导通电阻的数值,然后与正向导通电阻阈值相比较,若产生的数值小于阈值则输出结果发生变化,否则不变;采用连续均匀分布的随机方法在反向击穿电压的变化区间内随机产生一个反向击穿电压的数值,然后与反向击穿电压阈值相比较,若产生数值小于阈值则输出结果发生变化,否则不变;
步骤6:采用步骤1到步骤5的方法在不同辐照强度下分别对正向导通电阻和反向击穿电压进行结果输出,得到一系列实际输出结果点;重复进行多次结果输出;得到同一辐照强度下的多个实际输出结果点;
步骤7:若实际输出结果点与额定输出结果点相同则认为可靠,不相同认为不可靠;将实际输出结果点进行重采样,重采样的采样点的可靠度不发生变化;
步骤8:对重采样后的数据根据同一辐照强度下正向导通电阻和反向击穿电压的实际输出结果点与额定输出结果点相比较,计算出该辐照强度下CPU的可靠度;进而得到不同辐照强度下的可靠度曲线。
进一步的,所述步骤7的方法为:
步骤7.1:输出结果点样本为x={x1,x2……xN};
步骤7.3:计算新样本的参数θ1,该参数为均值、分布中位数和标准差;
步骤7.4:采用步骤7.2和步骤7.3相同的方法获得更多新样本x2,x3…xn…,并相应的计算出参数θ2,θ3…θn…,得到参数集合θ={θ1,θ2,θ3…θn};
步骤7.5:计算参数集合θ的算术平均值和标准差,其算术平均值和标准差在实现设定的区间内认为重采样成功,否则放弃本次重采样结果,重新进行重采样。
进一步的,所述步骤8采用当前辐照强度下左右相邻的w个辐照强度点数据纳入一起,计算出的可靠度为当前辐照强度点的可靠度。
进一步的,所述步骤3中正向导通电阻阈值t1=350欧姆,反向击穿电压t2=0.4V;
所述步骤4中正向导通电阻的变化区间为h1=50,反向击穿电压的变化区间为h2=0.3;
所述步骤5中以步长0.00000001的方式产生随机数。
随着人造卫星和航天飞机等方面的发展,集成电路在核辐射环境中的可靠性问题越来越受重视,本发明通过simulink仿真出驱动SIP电路中CPU板在辐射环境下的失效模型,通过误码率来度量其失效的过程,为解决在辐射环境这种苛刻条件下的实验时条件不足,数据缺乏等问题提供了更多的可行性和操作性,也为其他电路源件在特殊环境下的工作和失效机理的模型仿真研究提供了一些思路。
附图说明
图1为实现辐照强度的不断上升以及计算误码率仿真结果的流程图。
图2为本发明仿真结果误码率变化曲线图。
图3为CPU在辐射环境下概率密度曲线仿真结果图。
图4为CPU在辐射环境下可靠度曲线仿真结果图。
具体实施方式
本发明的辐射环境下一种CPU板失效仿真模型包括如下步骤:
1.针对驱动SiP典型电路,根据其辐射可靠性建模仿真与计算数据需要,设计独特的辐射可靠性实验系统及方案,涉及的实验设备加工与制造计划通过外协完成,实验执行或操作计划派遣研究人员参与完成。
2.将误码率作为量化CPU板失效的实验数据,并对各个电路板的数据进行处理,分类。
根据驱动SiP典型电路中CPU板的组成电路及其功能,利用simulik设计出能够实现CPU板功能的仿真模型。
3.利用matlab和simulink的互通,编写matlab代码模拟辐射环境下辐照强度不断上升的过程,其步骤如下:
步骤1,先通过实验所得数据通过各种数学模型拟合得到最适合误码率与辐照强度函数关系的曲线(最后选择多项式数学模型):
步骤2,如图3可得到电路板输出变化的阈值为正向导通电阻350欧姆(低于阈值使本应该输出的“0”变为“1”)和反向击穿电压0.4V(低于阈值使本应输出的“1”变为“0”)。
编写代码设置一个变化区间,正向导通电阻的变化区间为50,反向击穿电压的变化区间为0.3;并且给定这个区间会随着误码率的变化而变化,而误码率又因为和辐照强度有函数关系,从而使得变化区间随着辐照强度的增加不断变化,其变化为正向导通电阻取值区间不断向左移动,反向击穿电压取值区间不断向左移动,最后使得输出的结果错误不断增加,得到误码率的增加。
4.增加管子个数,减少误差并且多次仿真得到多条误码率曲线,然后通过Bootstrap重抽样法扩大抽样,最后通过与实验数据得到的CPU电路板失效的阈值比较计算,得到误码率的退化分布族。
Bootstrap重抽样法
Bootstrap重抽样法是Efron提出的一种新的统计方法,该方法即不需要有足够的样本量,同时也不需要获取参数的总体分布,而是通过重采样的方法来估计参数的分布。Bootstrap方法的巨大优势在于它的简便性,它可以直接计算分布中复杂参数的标准差和置信区间等。并且,它还是一种控制和检验结果稳定性的合理方法,因为该方法会比使用样本方差和正态性假设得到的标准区间更加准确。同时Bootstrap方法相较于其它线性抽样方法,值得注意的是,Bootstrap估计方法可分为:参数Bootstrap和非参数Bootstrap。参数Bootstrap和非参数Bootstrap的区别在于参数Bootstrap已经从已有的先验知识当中对总体分布进行了假设,然后从该分布中有放回地采样得到新的样本;而非参数Bootstrap则是在总体分布完全未知的情况下,采用重抽样方法从原有的样本中抽样进行参数统计推断。
Bootstrap方法是基于对原样本进行有放回的重抽样,在Bootstrap方法中由于总体采样的局限性,统计方差的估计往往是由原始测量样本的重抽样来代替总体近似得到的。虽然该方法在发展过程中经过了一些改进,但是基本原理都是类似的,下面给出该算法的基本步骤:
步骤1:给定一组样本容量为N的原始测量样本;
步骤2:从给出的原始测量样本x中进行有放回地重抽样得到样本容量同样为N的Bootstrap样本。
步骤3:通过得到的Bootstrap样本计算需要求的参数,这些参数为均值、分布中位数或标准差;
步骤4:重复n次步骤2和步骤3,则可得到n组需要计算的参数;
步骤5:在得到参数后估计出比较接近总体参数的算术平均值;
步骤6:参数的标准差用于表征估计值的不确定度;
6.利用核密度估计法得到概率密度曲线。
7.利用蒙特卡罗法,在退化分布族中将每一个辐照强度下电压与失效阈值相比较,最终得出可靠度曲线。
从图2中我们可以得到一条辐照强度与误码率的曲线,从图我们可以看出随着辐照强度地不断增加,误码率从0开始逐渐增加,最后超过一定辐照强度后误码率变为1,即为此时已经全部失效。
首先我们多次仿真得到多条如图2所示的误码率曲线,然后根据CPU板失效的误码率的阈值(此处设定为5%),得到图3所示的概率密度曲线,从图3我们可以看出在110.5到111的区间概率密度最大,于是我们可以得到CPU板临近失效时其辐照强度为多少。
最后从图4中我们可以看到SIP驱动中的CPU板在辐射环境下的可靠度曲线。可以得出在辐照强度为110左右CPU板开始失效,到111.5左右完全失效。
表一,具体实施方式获得辐照强度和误码率之间的对应数据
Claims (4)
1.一种辐射环境下CPU失效仿真方法,该方法包括:
步骤1:采用如下函数作为误码率与辐照强度的关系;
wuma=0.00129173x3-0.465764727x2+55.9605636x-2239.4157835
其中,wuma为误码率,x为辐照强度;
步骤2:分别设置CPU中正向导通电阻阈值t1和反向击穿电压阈值t2;
步骤3:针对正向导通电阻和反向击穿电压各设置一个变化区间,正向导通电阻的变化区间为h1,反向击穿电压的变化区间为h2;
步骤4:变化区间随着误码率的变化会发生移动,采用如下公式作为变化区间与误码率之间的关系;
Zmin=t1-h1×wuma Zmax=h1×(1-wuma)+t1
Fmin=t2-h2×wuma Fmax=h2×(1-wuma)+t2
Zmin~Zmax为正向导通电阻的变化区间,Fmin~Fmax为反向击穿电压的变化区间;
步骤5:采用连续均匀分布的随机方法在正向导通电阻的变化区间内随机产生一个正向导通电阻的数值,然后与正向导通电阻阈值相比较,若产生的数值小于阈值则输出结果发生变化,否则不变;采用连续均匀分布的随机方法在反向击穿电压的变化区间内随机产生一个反向击穿电压的数值,然后与反向击穿电压阈值相比较,若产生数值小于阈值则输出结果发生变化,否则不变;
步骤6:采用步骤1到步骤5的方法在不同辐照强度下分别对正向导通电阻和反向击穿电压进行结果输出,得到一系列实际输出结果点;重复进行多次结果输出;得到同一辐照强度下的多个实际输出结果点;
步骤7:若实际输出结果点与额定输出结果点相同则认为可靠,不相同认为不可靠;将实际输出结果点进行重采样,重采样的采样点的可靠度不发生变化;
步骤8:对重采样后的数据根据同一辐照强度下正向导通电阻和反向击穿电压的实际输出结果点与额定输出结果点相比较,计算出该辐照强度下CPU的可靠度;进而得到不同辐照强度下的可靠度曲线。
2.如权利要求1所述的一种辐射环境下CPU失效仿真方法,其特征在于所述步骤7的方法为:
步骤7.1:输出结果点样本为x={x1,x2……xN};
步骤7.3:计算新样本的参数θ1,该参数为均值、分布中位数和标准差;
步骤7.4:采用步骤7.2和步骤7.3相同的方法获得更多新样本x2,x3…xn…,并相应的计算出参数θ2,θ3…θn…,得到参数集合θ={θ1,θ2,θ3…θn};
步骤7.5:计算参数集合θ的算术平均值和标准差,其算术平均值和标准差在实现设定的区间内认为重采样成功,否则放弃本次重采样结果,重新进行重采样。
3.如权利要求1所述的一种辐射环境下CPU失效仿真方法,其特征在于所述步骤8采用当前辐照强度下左右相邻的w个辐照强度点数据纳入一起,计算出的可靠度为当前辐照强度点的可靠度。
4.如权利要求1所述的一种辐射环境下CPU失效仿真方法,其特征在于所述步骤3中正向导通电阻阈值t1=350欧姆,反向击穿电压t2=0.4V;
所述步骤4中正向导通电阻的变化区间为h1=50,反向击穿电压的变化区间为h2=0.3;
所述步骤5中以步长0.00000001的方式产生随机数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010454765.8A CN111695315B (zh) | 2020-05-26 | 2020-05-26 | 一种辐射环境下cpu失效仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010454765.8A CN111695315B (zh) | 2020-05-26 | 2020-05-26 | 一种辐射环境下cpu失效仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111695315A CN111695315A (zh) | 2020-09-22 |
CN111695315B true CN111695315B (zh) | 2023-06-23 |
Family
ID=72478169
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010454765.8A Active CN111695315B (zh) | 2020-05-26 | 2020-05-26 | 一种辐射环境下cpu失效仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111695315B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0517238A1 (en) * | 1991-06-06 | 1992-12-09 | Hitachi, Ltd. | Method of and apparatus for estimating remaining service life of material being exposed to irradiation |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2979708B1 (fr) * | 2011-09-06 | 2014-06-06 | Eads Europ Aeronautic Defence | Procede de caracterisation de la sensibilite d'un composant electronique soumis a des conditions d'irradiation |
-
2020
- 2020-05-26 CN CN202010454765.8A patent/CN111695315B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0517238A1 (en) * | 1991-06-06 | 1992-12-09 | Hitachi, Ltd. | Method of and apparatus for estimating remaining service life of material being exposed to irradiation |
Non-Patent Citations (1)
Title |
---|
Radiation effect microscopy for failure analysis of microelectronic devices;G. Vizkelethy;Nuclear Instrument and Methods in Physics Research Section B: Beam Interaction with Materials and Atoms;第467-475页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111695315A (zh) | 2020-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100568249C (zh) | 用基本原理仿真辅助半导体制造过程的系统和方法 | |
US8806418B1 (en) | Scaled sigma sampling | |
KR20140024252A (ko) | 희귀 사건 실패율을 식별하기 위한 방법 및 시스템 | |
Yildiz et al. | Learning stochastic differential equations with gaussian processes without gradient matching | |
Wang et al. | Inference for Weibull competing risks model with partially observed failure causes under generalized progressive hybrid censoring | |
CN102663495B (zh) | 一种用于非线性器件建模的神经网络数据生成方法 | |
CN111695315B (zh) | 一种辐射环境下cpu失效仿真方法 | |
CN106649920A (zh) | 一种基于ibis的集成电路总剂量效应建模方法 | |
Martin et al. | Stochastic approximation and Newton's estimate of a mixing distribution | |
Solomentsev et al. | Diagnostics programs efficiency analysis in operation system of radioelectronic eguipment | |
Canelas et al. | Yield optimization using k-means clustering algorithm to reduce Monte Carlo simulations | |
CN108932197A (zh) | 基于参数Bootstrap重抽样的软件失效时间预测方法 | |
Dong et al. | Prognostics 102: efficient Bayesian-based prognostics algorithm in Matlab | |
CN106485392A (zh) | 互联发电系统可靠性评估方法 | |
Zaka et al. | Control charts for the shape parameter of power function distribution under different classical estimators | |
Krishnan et al. | Stochastic behavioral modeling of analog/mixed-signal circuits by maximizing entropy | |
CN114417653B (zh) | 基于不确定性的卫星组件电路板热可靠性分析方法 | |
CN113221460B (zh) | 一种基于神经网络回归的单粒子瞬态效应建模方法 | |
Qi et al. | The Device Compact Model Based on Multi-gradient Neural Network and Its Application on MoS 2 Field Effect Transistors | |
Xiong et al. | Analysis of a simple step-stress life test with a random stress-change time | |
CN111651337B (zh) | 一种sram存储器空间服役故障分类失效检测方法 | |
Peng et al. | Developing a reliability model of CNC system under limited sample data based on multisource information fusion | |
CN112667957A (zh) | 一种基于深度神经网络的智能电能表失效率预测方法 | |
CN108491606A (zh) | 一种材料强度分布获取方法 | |
Song et al. | A novel sequential sampling algorithm for reliability assessment of microgrids |
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 |