CN106054601B - 确定防滑刹车控制装置低温故障分布的方法 - Google Patents
确定防滑刹车控制装置低温故障分布的方法 Download PDFInfo
- Publication number
- CN106054601B CN106054601B CN201610373928.3A CN201610373928A CN106054601B CN 106054601 B CN106054601 B CN 106054601B CN 201610373928 A CN201610373928 A CN 201610373928A CN 106054601 B CN106054601 B CN 106054601B
- Authority
- CN
- China
- Prior art keywords
- distribution
- low temperature
- triode
- failure
- formula
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
一种确定防滑刹车控制装置低温故障分布的方法,随机抽取防滑刹车控制装置中的元器件进行低温故障测试,根据测试数据进行计算,根据计算结果确定元器件的低温故障分布,作为工程中选择故障分布公式的依据,消除现有技术带来的误差和资源消耗。本发明确定采用威布尔分布评估防滑刹车控制装置的低温故障,证明防滑刹车控制装置的故障是耗损故障,纠正了现有技术中将耗损故障作为偶然故障处理的错误,以及在障原因、故障机理方面产生误差。
Description
技术领域
本发明涉及民用运输类飞机电子产品的低温故障分析领域,具体是一种确定防滑刹车控制装置低温故障分布的方法。
背景技术
现有技术进行电子产品的可靠性预计时,默认电子产品的故障服从指数分布,不管是什么故障模式,都这样处理。
防滑刹车控制装置是飞机电子防滑刹车系统中的附件,由飞机提供电源,接收机轮速度传感器感受的机轮转速变化电信号,防滑刹车控制装置根据该电信号进行着陆防滑刹车过程中的刹车压力控制。该防滑刹车控制装置能够完成正常着落防滑刹车控制,起飞线刹车控制,车轮间保护控制,接地保护控制,起落架收上后的机轮止转刹车功能。
国外现状:
国外一直采用指数分布进行电子产品的故障率预计,采用的标准有:美国标准MIL-HDBK-217F《电子设备可靠性预计手册》,该手册中有两种电子设备的故障率预计方法。
1)元器件计数法
元器件计数法采用电子元器件的基本故障率进行电子设备的故障率预计,在电子设备研制过程中采用,故障分布确定为指数分布。元器件计数法的数学模型为:
式中:
λEQUIP:电子设备的故障率;故障率×10-6 1/h;
λg:电子元器件的基本故障率;故障率×10-6 1/h;
πQ:第i种电子元器件的质量系数;
Ni:第i种电子元器件的数量;
n:设备所用电子元器件的种类数。
采用(1)式计算电子设备的故障率,由于按指数分布处理,电子设备的平均故障间隔时间MTBF计算方法为:
MTBF=1/λEQUIP (2)
采用(2)式计算电子设备的平均故障间隔时间MTBF。
在(1)式中电子设备的故障率仅考虑质量系数,未考虑环境和工作电应力的影响。
2)元器件应力法
与模型(1)不同,不直接计算电子设备的故障率,而是对每一种元器件建立具体的计算模型,不同种类的元器件建立不同种类的计算模型,在模型中综合考虑质量系数、电应力系数、环境系数,以CMOS电路为例的故障率计算模型:
λp=λBDπMFGπTπCD+λBPπEπQπPT+λEOS (3)
式(3)中:
λp:该CMOS电路的工作故障率;
λBD:该CMOS电路的基本故障率;
πMFG:该CMOS电路的制造工艺系数;
πT:该CMOS电路的温度系数;
πCD:该CMOS电路的复杂系数;
πE:该CMOS电路的使用环境系数;
πQ:该CMOS电路的质量系数;
πPT:该CMOS电路的封装类型系数;
λEOS:该CMOS电路的过载故障率。
在MIL-HDBK-217F标准中有不同元器件的故障率计算模型,将该设备中的所有电子元器件均按照MIL-HDBK-217F标准确定的模型计算在工作状态下和环境条件下的故障率,计算完成后,将所有电子元器件的故障率累加,就得到了该电子设备的故障率。
该故障率是基本故障率,未对故障原因进行区分。
国内现状:
国内可靠性工作起步较晚,一直参照国外的可靠性工作经验开展相应工作。国内也采用指数分布进行电子产品的故障分析、故障率预计,采用的标准有:GJB/Z299C《电子设备可靠性预计手册》,该手册中有两种电子设备的故障率预计方法:
1)元器件计数法
在GJB/Z299C中,元器件计数法采用电子元器件的基本故障率进行电子设备的故障率预计,这种方法在电子设备研制过程中采用,故障分布函数为指数分布。元器件计数法的数学模型为:
式中:
λGS:电子设备的故障率;故障率×10-6 1/h;
λGi:第i种电子元器件的通用故障率;故障率×10-6 1/h;
πQi:第i种电子元器件的通用质量系数;
Ni:第i种电子元器件的数量;
n:设备所用电子元器件的种类数。
在(4)式中电子设备的故障率仅考虑质量系数,未考虑环境和工作电应力的影响。
采用(4)式计算电子设备的故障率,由于按指数分布处理,电子设备的平均故障间隔时间MTBF计算方法为:
MTBF=1/λGS (5)
采用(5)式计算电子设备的平均故障间隔时间MTBF。
2)元器件应力法
在GJB/Z299C中,与模型(4)不同,不直接计算电子设备的故障率,而是对每一种元器件建立计算模型,不同种类的元器件建立不同种类的计算模型,在模型中综合考虑质量系数、电应力系数、环境系数,以单片数字电路为例的故障率计算模型为例:
λp=πQ[C1πTπV+(C2+C3)πE]πL (6)
式(3)中:
λp:该单片数字电路的工作故障率;
πQ:该单片数字电路的质量系数;
C1:单片数字电路复杂度故障率;
πT:单片数字电路温度应力系数;
πV:单片数字电路电压应力系数;
C2:单片数字电路电路复杂度故障率;
C3:单片数字电路封装复杂度故障率;
πE:该单片数字电路的使用环境系数;
πL:该单片数字电路的成熟系数;
在GJB/Z299C标准中有不同元器件的故障率计算模型,将该设备中的所有电子元器件均按照GJB/Z299C标准规定的模型计算在工作状态下和环境条件下的故障率,计算完成后,按照电子设备的故障服从指数分布,将所有电子元器件的故障率累加,就得到了该电子设备的故障率。
该故障率是基本故障率,未对故障原因进行区分。
国内外现有技术的共同特点是:
1)元器件计数法在研制初期缺少使用环境要求的条件下,粗略估计电子设备的故障率,按照电子设备的故障服从指数分布;
2)元器件应力法在已知使用环境要求和工作电应力的条件下计算电子设备的故障率,且假设电子设备的故障服从指数分布;
3)可靠性预计模型没有区分故障原因。
国内外现有技术的共同特点是将电子元器件的故障按照指数分布处理,且都颁布了标准。
现有技术的优点:国内外标准收集了大量的电子元器件故障率数据,利用手册数据和指数分布处理工程问题非常方便。
现有技术的缺点:
1)大量的试验和使用数据证明,电子元器件的故障并不一定服从指数分布,在不服从指数分布的情况下,用指数分步进行处理就带来很大误差;
2)在进行故障预计的过程中,未进行故障原因区分。
发明内容
为克服现技术采用指数分布分析电子产品的故障存在很大误差、且未区分故障原因的缺点,本发明提出了一种确定防滑刹车控制装置低温故障分布的方法。
本发明的具体过程是:
步骤1,确定电子元器件的低温故障分布测试方案
第一步,确定元器件低温抽样测试的种类
以半导体分立元器件作为低温测试样件,抽取三极管进行低温故障分布测试。
第二步,随机抽取同型号多个3AX31B三极管作为试验样件;各三极管的使用温度范围是:-55℃~125℃。
第三步,进行3AX31B三极管低温故障分布测试
首先确定3AX31B三极管低温故障原因;
其次,根据低温故障原因确定环境试验箱的温度为-55℃。将环境试验箱的温度降到-55℃并恒温,然后同时将各3AX31B三极管放进低温试验箱,通电工作条件下进行测试。测试条件为:施加共基极集电极-基极直流电压-6V,进行到各3AX31B三极管全部发生低温故障时测试结束。
第四步,确定累积失效概率计算公式
在样本数小于50个的条件下,通过近似中位秩公式(15)计算累积失效概率。
公式(15)中:Fn(ti)是试验中元器件的累积失效概率;i是第i个失效的元器件;0.32和0.36均为公式中的常数;n是试验元器件的总数。
第五步,对威布尔分布公式进行线性化处理
当该3AX31B三极管的低温故障分布测试数据回归分析相对于威布尔分布线性相关系数r小于0.9时,该3AX31B三极管的低温故障试验数据与威布尔分布拟合不好,对该3AX31B三极管的低温故障试验数据采用正态分布重新计算。
列出威布尔分布数学式:
式中:
t为时间;m为威布尔分布的形状参数,m值的大小表示故障分布的分散程度;t0为尺度参数,缩小和放大分布函数坐标尺度的作用;Fn(t)是威布尔分布的累积失效函数。
对式(7)等号两边取两次以e为底的自然对数线性化处理,得到(8)式:
列线性方程:y=bx+a (9)
(9)式是解析几何中常用的线性方程,y为bx+a的函数,a为x=0时方程(9)在y轴上的截距,b为系数。
依据对(9)式对(8)式进行变量转换:
ln是自然对数符号。
至此,完成了分布函数的线性化处理,具备在双对数坐标系中对三极管的低温测试数据进行计算的条件。
第六步,确定截距a、系数b和相关系数r:
采用最小二乘法确定截距a、系数b,定义a为x=0时方程(9)在y轴上的截距,b为系数。根据式(9)、(10),采用工程数学中的最小二乘法确定回归系数a、b的计算公式为:
工程中,的值等于b。为威布尔分布的形状参数m的估计值。
x=0时方程(9)在y轴上截距a的计算公式为:
威布尔分布的尺度参数t的计算公式为:
通过公式(14)确定相关系数r:
在威布尔分布双对数坐标系中计算三极管的低温故障数据,若服从威布尔分布,计算结束;若不服从威布尔分布,则重新在正态分布的双对数坐标系中计算三极管的低温故障数据;若不服从正态分布,则该三极管的研制质量存在隐患,应在进行质量改进后重新进行低温故障分布测试和计算。
步骤2,对三极管进行低温故障分布测试
按照步骤1确定的测试方案,进行各3AX31B三极管的低温故障分布测试,测试进行到各3AX31B三极管全部失效时结束。得到各3AX31B三极管的低温故障分布测试数据。
所述各3AX31B三极管的低温故障分布测试过程是:
1)将温度箱的温度降到-55℃;
2)将各3AX31B三极管放进温度箱并关闭箱门,施加共基极集电极-基极直流电压-6V,使各3AX31B三极管都通电工作;
3)当共发射极直流电流放大系数hFE小于等于10时确定为该三极管发生低温故障,测试结束。
步骤3,对各3AX31B三极管的低温故障分布分布测试数据进行计算
根据步骤2测试得到的各3AX31B三极管的低温故障数据,对这各测试数据进行计算,确定低温测试数据的故障分布。
根据各3AX31B三极管的低温故障分布测试数据列出最小二乘法计算表。
所述最小二乘法的回归计算表的填写过程为:
第一列为ti,具体为各三极管的低温故障时间数据,按照测试数据由小到大的原则在表1中排为同1列;第二列为通过近似中位秩公式得到的Fn(ti);第三列为lnti,即xi;第四列为yi,即lnln1/1-Fn(t i);第五列填写第六列填写第七列填写第八列填写第九列填写
其中:
lnti是对测试时间取对数,在线性方程(9)、(10)中用xi表示;xi是双对数坐标系中的横轴,i=1,2,3……n,表示第i个元器件发生故障的时间;
lnln1/1-Fn(ti)是对1/1-Fn(ti)取两次对数,在线性方程(12)中采用yi表示;其中Fn(ti)是故障时间的累积失效函数,取两次对数的工程含义是对威布尔分布公式进行线性化处理,i=1,2,3……n,表示第i个元器件发生故障的时间;
是(14)式中的计算部分,(14)式用于计算相关系数r,根据相关系数r确定各3AX31B三极管的低温故障数据是否服从威布尔分布。
各3AX31B三极管的测试数据按照最小二乘法计算方法,计算后全部填入所述最小二乘法的回归计算表。
步骤4,确定3AX31B三极管的低温故障分布
根据对3AX31B三极管低温故障测试数据的威布尔分布回归计算结果,若3AX31B三极管的低温故障服从威布尔分布,计算结束。
若经计算不服从威布尔分布,则应重新按照正态分布计算,若服从正态分布,计算结束。
若经计算不服从正态分布,则应对该3AX31B三极管进行质量改进,改进后重新进行计算,直至确定3AX31B三极管服从的概率分布,计算结束。
本发明随机抽取防滑刹车控制装置中的元器件进行低温故障测试,根据测试数据进行计算,根据计算结果确定元器件的低温故障分布,作为工程中选择故障分布公式的依据,消除现有技术带来的误差和资源消耗。
由于电子产品是在低温条件下发生的故障,因此制定低温测试方案,采用数理统计中的抽样方法,用样件的低温测试数据代表车防滑控制盒的低温测试数据,对取得的低温故障测试数据进行回归分析,得到电子产品的低温故障分布。
由于防滑刹车控制装置的低温故障是元器件发生的,本发明抽取该控制盒的元器件进行低温故障测试。
本发明确定了现有国内外标准采用指数分布进行电子产品故障率预计的方法与工程现实不符,当电子产品设计要求预计高精度的故障率时,应选择采用其他精度高的概率分布公式,使可靠性预计数据的精度达到设计要求。
为证明本发明的效果,本发明对一组分别采用威布尔分布和指数分布计算在-55℃条件下测试1000h的瞬时故障率数据,计算结果相差1.2257×109倍,证明本发明提出的根据使用要求确定故障分布的方法符合工程现实,避免由于计算误差带来的损失;证明本发明效果的计算过程为:
第一步,采用指数分布计算该3AX31B三极管在-55℃条件下测试1000h的故障率
当威布尔分布公式的位置参数为零,形状参数为1时就等于指数分布。3AX31B三极管使用到1000h的故障率为:
第二步,采用威布尔分布计算该三极管在-55℃试验1000h的故障率
本发明实施例通过低温测试已经证明,该三极管的形状参数m=3.6961,威布尔分布的位置参数为零,三极管使用到1000h的故障率为:
第三步,瞬时故障率计算对比
产品的故障率分有瞬时故障率和恒定故障率之分,在恒定故障率的条件下,故障率为常数。本发明故障率的大小是时间的函数,因此在进行故障率比较时,在相同的时间点上进行故障率大小的比较,是瞬时故障率。
λ(t威布尔)/λ(t指数)=0.1198×10-5/0.9774×10-14=1.2257×109
即在使用1000h的条件下,采用威布尔分布计算得到的故障率是指数分布的1.2257×109倍,采用指数分布得到的故障率与产品实际的故障率相差大,不宜作为精确计算方法。
与现有技术相比,本发明取得的有益效果为:
采用本发明技术进行防滑刹车控制装置的故障分布计算,得到了该防滑刹车控制装置低温故障的威布尔分布参数,形状参数的具体数值为3.6961,是大于1的数,不服从指数分布。确定了在对防滑刹车控制装置的低温故障分布参数进行分析时,采用威布尔分布预计防滑刹车控制装置的可靠性指标,能够得到平均故障前时间。
本发明和现有标准方法实施效果的综合比较:
1、现有国内外标准采用指数分布进行电子产品的可靠性预计,指数分布的故障率为常数,不适用于计算电子产品低温耗损故障的公式。
2、采用本发明确定的威布尔分布进行可靠性预计工作,能够提高电子产品的可靠性预计值与实际的可靠性指标的符合程度;
3、本发明确定采用威布尔分布计算防滑刹车控制装置的低温故障分布,测试数据证明防滑刹车控制装置的故障不是偶然故障,而是耗损故障。在可靠性工程领域产品的故障分为:早期故障、偶然故障、耗损故障,早期故障指使用初期出现的故障,偶然故障指偶然因素引起的故障。耗损故障指磨损、老化、腐蚀、疲劳等物理化学原因引起的故障。现有标准将电子产品的耗损故障作为偶然故障处理,在故障数据计算中产生误差,本发明的效果是纠正了这种差错。
具体实施方式
实施例1
本实施例以一种飞机防滑刹车控制装置为例,确定该防滑刹车控制装置的低温故障分布。
该防滑刹车控制装置的低温故障由电子元器件的低温故障引起,因此,确定该防滑刹车控制装置的低温故障分布就是确定组成该防滑刹车控制装置电子元器件的低温故障分布。确定所有电子元器件低温故障分布的测试时间长,费用高,在工程中也没有将全部电子元器件进行故障分布测试的先例。本实施例依据数理统计原理,采用抽样的方法确定防滑刹车控制装置的低温故障分布,就是用子样的故障分布表示母体的故障分布,子样指进行低温测试的元器件,母体指防滑刹车控制装置。
常用的故障分布有:
1)指数分布,故障率为恒定型,这是美国标准MIL-STD-217F、中国标准GJB299采用的故障分布;指数分布是概率论与数理统计中的一种常用分布。
2)正态分布,故障率为时间的增函数,适用于老化、腐蚀、磨损、疲劳耗损原因引起的故障,正态分布有方差和均值两个参数。
3)威布尔分布,故障率为时间的增函数,适用于老化、腐蚀、磨损、疲劳耗损原因引起的故障,威布尔分布有形状参数m、尺度参数t0、位置参数γ三个参数。
本实施例根据电子元器件的低温故障为机械故障,首先假设其故障服从威布尔分布,并进行电子元器件的低温故障分布测试,根据测试数据进行线性回归分析,若线性回归分析结果服从威布尔分布,线性回归分析结束。若线性回归分析不服从威布尔分布,则重新进行线性回归分析,直至得到所服从的故障分布。回归分析前应对威布尔分布公式进行线性化处理,线性化处理后的分析称为线性回归分析。
本实施例采用的“分布”是一个来自概率论与数理统计的常用术语,含义是任何一个随机变量都用概率分布的形式表示,说明示例:同一类产品早期发生的故障少,达到长寿命的也少,寿命集中在平均值附近的居多,使用时间是自变量,函数值的大小是因变量,统称为故障的密度函数。常用耗损故障的概率分布有正态分布、威布尔分布,分布的物理含义是产品的寿命不是确定值,而是随机变量,在该产品失效前只能预计该产品可能的失效时间范围,而不能确定该产品准确的失效时间。
步骤1,确定电子元器件的低温故障分布测试方案
第一步,确定元器件低温抽样测试的种类
现有元器件分为半导体分立元器件、半导体集成元器件两大类,分立、集成的制造方式根据设计用途确定,和元器件的低温故障原因无关,在低温条件作用下,分立、集成元器件都会发生故障,故障原因同为低温引起的材料性能飘移,严重时发生故障,这是半导体材料的物理属性。所以,抽样取分立元器件还是集成器件不影响故障分布的测试效果。本发明选择半导体分立元器件作为低温测试样件。
在防滑刹车控制装置的元器件中,在低温条件下出现故障的是三极管,当三极管出现低温故障时,防滑刹车控制装置在低温条件下不能正常工作。
在半导体分立元器件中,按照随机抽样原则抽取三极管进行低温故障分布测试。
第二步,随机抽取同型号4个3AX31B三极管作为试验样件
为了提高所述3AX31B三极管低温分布参数计算结果的精度,随机从同批次三极管中抽取四个作为测试样件,各三极管的使用温度范围是:-55℃~125℃,即所述各3AX31B三极管在-55℃~125℃的范围内都能够正常工作,是该种类3AX31B三极管出厂检验的合格判据之一。根据线性回归分析原理,在平面上二点确定一条直线,工程意义是最少抽取二个三极管,在原理上二个3AX31B三极管的低温故障数据就能够确定该3AX31B三极管的低温故障分布。但是,由于元器件制造过程中不可避免产生偏差,为了提高故障测试的精度,将抽取样件的数量增加一倍,共抽取四个三极管进行低温故障分布测试。
抽取三极管的型号是:3AX31B,其中3表示三极管;A表示PNP型锗材料低频小功率管;31、B:均为序号。该三极管发射极开路、集电极-基极反向截止电流ICBO小于等于10μA,基极开路、集电极-发射极反向截止电流ICEO小于等于750μA,VCE(sat):共发射极反向饱和电压小于等于0.65V;hFE:共发射极直流电流放大系数为50~150。
第三步,进行3AX31B三极管低温故障分布测试
本步骤的测试方案属于子样的测试方案,因为将所有元器件进行低温故障测试费用高,时间长。因此采用数理统计抽样测试的方案,用所抽取元器件的低温故障分布代表防滑刹车控制装置的低温故障分布,因为元器件的材料同为半导体,同为低温条件作用下发生的故障。在防滑刹车控制装置的元器件中,在低温条件下出现故障的是三极管,抽取三极管作为低温测试的元器件。
首先确定3AX31B三极管低温故障原因,包括:
1)低温条件下热载流子效应加剧,放大系数急剧下降,性能飘移,严重时发生低温故障;
2)低温焊点开裂;
3)低温凝露发生击穿或短路;
4)低温管脚脆断。
其次,根据低温故障原因确定环境试验箱的温度为-55℃。将环境试验箱的温度降到-55℃并恒温,然后同时将各3AX31B三极管放进低温试验箱,通电工作条件下进行测试。测试条件为:施加共基极集电极-基极直流电压-6V,进行到4个3AX31B三极管全部发生低温故障时测试结束。
第四步,确定累积失效概率计算公式
在样本数小于50个的条件下,通过近似中位秩公式计算累积失效概率。
近似中位秩公式:
在所述近似中位秩公式(15)中:
Fn(ti):试验中元器件的累积失效概率;
i:第i个失效的元器件;本实施例中,i=1~4;
0.32:公式中的常数;
0.36:公式中的常数。
n:试验元器件的总数,在本实施例中n=4。
第五步,对威布尔分布公式进行线性化处理
因为本发明列出的公式(7)为非线性,该公式是概率论与数理统计中的公式,不能直接进行回归分析,因此需要进行线性化处理。线性化处理的对象是威布尔公式。
经线性化后的威布尔公式,用于进行3AX31B三极管的低温故障测试数据的回归分析。3AX31B三级管是从该防滑刹车控制装置抽取进行低温故障分布测试的元器件。
本发明确定在分析过程中,当对3AX31B三极管的低温故障测试数据回归分析相对于威布尔分布的线性相关系数r大于等于0.9时,这些数据与威布尔分布的线性相关好,并且线性相关系数r越大测试数据与威布尔分布的拟合程度越好。
本发明确定,当该3AX31B三极管的低温故障分布测试数据回归分析相对于威布尔分布线性相关系数r小于0.9时,就认为该3AX31B三极管的低温故障试验数据与威布尔分布拟合不好,对该3AX31B三极管的低温故障试验数据采用正态分布重新计算。
若采用正态分布对该3AX31B三极管的低温故障测试数据回归分析结果,线性相关系数也小于0.9,则该批产品质量一致性不好,应进行质量改进,改进后重新进行该3AX31B三极管的低温测试和故障数据的计算工作。
列出威布尔分布数学式:
由
式中:
t为时间;m为威布尔分布的形状参数,m值的大小表示故障分布的分散程度;t0为尺度参数,缩小和放大分布函数坐标尺度的作用;Fn(t)是威布尔分布的累积失效函数。
根据式(7)等号两边取两次自然对数线性化处理,自然对数是以e为底的对数,得到(8)式:
列线性方程:y=bx+a (9)
(9)式是解析几何中常用的线性方程,y为bx+a的函数,a为x=0时方程在y轴上的截距,b为系数。
依据对(9)式对(8)式进行变量转换:
ln是自然对数符号,自然对数是以e为底的对数。
至此,完成了分布函数的线性化处理,具备在双对数坐标系中对三极管的低温测试数据进行计算的条件。
第六步,确定截距a、系数b和相关系数r。
在完成(10)式的变量转化后,线性方程(9)中a、b的物理意义是使三极管低温测试数据在威布尔分布双对数坐标系中的数据点与线性方程(9)所确定的回归直线距离最小,计算精度最高。
因为本发明是在假设三极管的低温测试数据服从威布尔分布的条件下完成计算过程的,最小二乘法具有判别测试数据是否服从威布尔分布的功能,当相关系数r大于等于0.9时,三极管的低温测试数据与回归直线的符合性好,三极管的低温测试数据服从威布尔分布;当相关系数r小于0.9时,三极管的低温测试数据与回归直线的符合性不好,三极管的低温测试数据不服从威布尔分布;在三极管的低温测试数据不服从威布尔分布的条件下,应采用正态分布重新计算三极管的低温测试数据。
根据双对数坐标系的原理,按照线性方程中{x,y}的对应关系,在双对数坐标系中相应列出试验数据根据式(10),采用工程数学的最小二乘法确定b:
工程中,的值等于b的值。为威布尔分布的形状参数m的估计值;
x=0时方程(9)在y轴上截距a的计算公式为:
尺度参数t的估计值为:
(13)式中为t的估计值,用表示针对尺度参数t进行计算后的估计值。e是自然对数中的e,在最小二乘法计算过程中,采用(7)式两边取自然对数的方法进行线性化处理。
本发明中相关系数r:
相关系数
至此,所确定的该防滑刹车控制装置电子元器件低温故障测试方案包括下列内容:
1)抽检测试的元器件种类为3AX31B三极管;
2)抽检测试的3AX31B三极管数量为4个;
3)测试结束的条件是3AX31B 4个三极管全部发生低温故障;
4)选取累积失效计算公式为近似中位秩公式;
5)在威布尔分布双对数坐标系中计算三极管的低温故障数据,若服从威布尔分布,计算结束;若不服从威布尔分布,则重新在正态分布的双对数坐标系中计算三极管的低温故障数据;若不服从正态分布,则该三极管的研制质量存在隐患,应在进行质量改进后重新进行低温故障分布测试和计算。
步骤2,对三极管进行低温故障分布测试
按照本发明步骤1确定的测试方案,进行4个3AX31B三极管的低温故障分布测试,测试进行到4个3AX31B三极管全部失效时结束。
4个3AX31B三极管的低温故障分布测试过程是:
1)将温度箱的温度降到-55℃;
2)将4个3AX31B三极管放进温度箱并关闭箱门,施加共基极集电极-基极直流电压-6V,使4个3AX31B三极管都同时通电工作;
3)当共发射极直流电流放大系数hFE小于等于10时确定为该三极管发生低温故障,测试结束。
经测试,4个三极管试验的试验发生故障的时间分别为:3840h,4320h,4560h,6480h。
步骤3,对4个3AX31B三极管的低温故障分布测试数据进行计算
根据本发明步骤2测试得到的4个3AX31B三极管的低温故障数据,本步骤对这4个测试数据进行计算,确定低温测试数据的故障分布。
根据4个3AX31B三极管的低温故障分布测试数据列出最小二乘法计算表。
为了方便计算,本发明列出最小二乘法的回归计算表,见表1。填写过程为:
第一列为ti,具体为4个三极管的低温故障时间数据,按照测试数据由小到大的原则在表1中排为同1列;第二列为通过近似中位秩公式计算得到的Fn(ti);第三列为lnti,即xi;第四列为yi,即lnln1/1-Fn(t i);第五列填写第六列填写第七列填写第八列填写第九列填写
其中:lnti是对测试时间取对数,在线性方程(9)、(10)中用xi表示;xi是双对数坐标系中的横轴;
lnln1/1-Fn(ti)是对1/1-Fn(ti)取两次对数,在线性方程(12)中采用yi表示;其中Fn(ti)是故障时间的累积失效函数,取两次对数的工程含义是对威布尔分布公式进行线性化处理;
是(14)式中的计算部分,(14)式用于计算相关系数r,根据相关系数r确定4个3AX31B三极管的低温故障数据是否服从威布尔分布。在概率论与数理统计中,测试数据与概率分布的符合性,常用术语“服从”表示。
4个3AX31B三极管的测试数据按照最小二乘法计算方法,计算后全部填入表1。
表1测试温度-55℃,共4件,试验数据威布尔分布回归分析,时间ti单位:h
(14)式是引用最小二乘法确定相关系数r的计算公式,属于工程数学中常用的方法。根据表1数据,由(14)式计算相关系数r的数值:
相关系数
相关系数r为0.9,3AX31B三极管的低温故障服从威布尔分布。继续评估威布尔分布参数。
由(11)式确定形状参数m的估计值:
公式(9)中,b为系数。
工程中,的值等于b的值。为威布尔分布的形状参数m的估计值。
由(12)式得x=0时方程(9)在y轴上的截距a的计算数值为::
由(13)式得尺度参数t的估计值为:
本步骤达到了下列目的:
1)利用3AX31B三极管低温故障测试数据,计算与威布尔分布的线性相关系数r的估计值为0.9,服从威布尔分布;
2)得到了3AX31B三极管低温故障的威布尔分布参数:形状参数m=3.6961,尺度参数t0=1.0231×1014。
步骤4,确定3AX31B三极管的低温故障分布
在可靠性工程中,常用概率分布表示故障发生的多少和时间的关系,分布是概率论与数理统计领域的常用语,概率论与数理统计是可靠性工程的基础之一。举例而言,在3AX31B三极管的低温故障服从威布尔分布的条件下,3AX31B三极管低温寿命短、低温寿命长的数量少,而在平均寿命时间的3AX31B三极管数量多,寿命的长短和时间之间存在随机函数关系,这种随机函数用概率分布的公式表示。
根据对3AX31B三极管低温故障测试数据的威布尔分布回归计算结果,若3AX31B三极管的低温故障服从威布尔分布,计算结束。
若经计算不服从威布尔分布,则应重新按照正态分布计算,若服从正态分布,计算结束。
若经计算不服从正态分布,则应对该3AX31B三极管进行质量改进,改进后重新进行计算,直至确定3AX31B三极管服从的概率分布,计算结束。
在电子元器件研制质量稳定的条件下,根据抽样测试的方法,用子样的故障分布表示母体的故障分布。不采用将所有元器件都进行低温故障测试的方法。
本发明确定,采用威布尔分布公式对3AX31B三极管低温故障测试数据进行回归分析的计算结果,3AX31B三极管低温故障与威布尔分布的相关系数为0.9,3AX31B三极管的低温故障服从威布尔分布。
本发明确定,3AX31B三极管低温故障的形状参数m=3.6961,故障率为递增型,而不是恒定型。所以,3AX31B三极管的低温故障不服从指数分布,指数分布的形状参数m=1,即故障率为恒定型。
本发明确定,对于由低温故障引起的电子元器件故障,确定为服从威布尔分布。
根据本发明结果,按照威布尔分布进行电子产品的故障分析。
Claims (4)
1.一种确定防滑刹车控制装置低温故障分布的方法,其特征在于,具体过程是:
步骤1,确定电子元器件的低温故障分布测试方案:
第一步,确定元器件低温抽样测试的种类,以半导体分立元器件作为低温测试样件;在半导体分立元器件中,抽取三极管进行低温故障分布测试;
第二步,随机抽取同型号多个三极管作为试验样件;各三极管的使用温度范围是:-55℃~125℃;
第三步,进行三极管低温故障分布测试,首先确定三极管低温故障原因;
其次,根据低温故障原因确定环境试验箱的温度为-55℃;将环境试验箱的温度降到-55℃并恒温;同时将各三极管放进低温试验箱,通电工作条件下进行测试;
测试条件为:施加共基极集电极-基极直流电压-6V,进行到各三极管全部发生低温故障时测试结束;
第四步,确定累积失效概率计算公式:
在样本数小于50个的条件下,通过近似中位秩公式计算累积失效概率;
公式(15)中:Fn(ti)是试验中元器件的累积失效概率;ti表示第i个试验中元器件的低温故障时间,i=1,2,3…n;0.32和0.36均为公式中的常数;n是试验元器件的总数,且n=4;
第五步,对威布尔分布公式进行线性化处理:
当该三极管的低温故障分布测试数据回归分析相对于威布尔分布线性相关系数r大于等于0.9时,认为服从威布尔分布;当线性相关系数r小于0.9时,该三极管的低温故障测试数据与威布尔分布拟合不好,对该三极管的低温故障测试数据采用正态分布重新计算;
列出威布尔分布数学式:
式中:t为尺度参数;m为威布尔分布的形状参数,m值的大小表示故障分布的分散程度;t0为参数,具有缩小和放大分布函数坐标尺度的作用;Fn(t)是威布尔分布的累积失效函数;
对式(7)等号两边取两次以e为底的自然对数线性化处理,得到:
列线性方程:y=bx+a (9)
y为bx+a的函数,a为x=0时方程(9)在y轴上的截距,b为系数;
依据(9)式对(8)式进行变量转换:
至此,完成了威布尔分布函数的线性化处理,具备在双对数坐标系中对三极管的低温测试数据进行计算的条件;
第六步,确定截距a、系数b和相关系数r:
按照双对数坐标系的原理,在双对数坐标系中相应列出测试数据;
采用最小二乘法确定截距a、系数b,通过公式(11)确定系数b:
工程中,的值等于b的值;为威布尔分布的形状参数m的估计值;
方程(9)在y轴上的截距a为:
威布尔分布尺度参数t估计值的计算公式为:
式中为t的估计值,用表示针对尺度参数t进行计算后的估计值;
通过公式(14)确定相关系数r:
在威布尔分布双对数坐标系中计算三极管的低温故障数据,若服从威布尔分布,计算结束;若不服从威布尔分布,则重新在正态分布的双对数坐标系中计算三极管的低温故障数据;若不服从正态分布,则该三极管的研制质量存在隐患,应在进行质量改进后重新进行低温故障分布测试和计算;
步骤2,对三极管进行低温故障分布测试:
按照步骤1确定的测试方案,进行各三极管的低温故障分布测试,测试进行到各三极管全部失效时结束;得到各三极管的低温故障分布测试数据;
步骤3,对各三极管的低温故障分布测试数据进行计算:
根据步骤2测试得到的各三极管的低温故障数据,对各测试数据进行计算,确定低温测试数据的故障分布;
根据各三极管的低温故障分布测试数据列出最小二乘法计算表;
步骤4,确定三极管的低温故障分布:
根据对三极管低温故障测试数据的威布尔分布回归计算结果,若三极管的低温故障服从威布尔分布,计算结束;
若经计算不服从威布尔分布,则应重新按照正态分布计算,若服从正态分布,计算结束;
若经计算不服从正态分布,则应对该三极管进行质量改进,改进后重新进行计算,直至确定三极管的低温测试数据服从的概率分布,计算结束。
2.如权利要求1所述确定防滑刹车控制装置低温故障分布的方法,其特征在于,所述各三极管的低温故障分布测试过程是:
Ⅰ将低温试验箱的温度降到-55℃;
Ⅱ将各三极管放进低温试验箱并关闭箱门,施加共基极集电极-基极直流电压-6V,使各三极管都同时通电工作;
Ⅲ当共发射极直流电流放大系数hFE小于等于10时确定为该三极管发生低温故障,测试结束。
3.如权利要求1所述确定防滑刹车控制装置低温故障分布的方法,其特征在于,所述最小二乘法的回归计算表的填写过程为:
第一列为ti,具体为各三极管的低温故障的时间数据,按照测试数据由小到大的原则在所述计算表中排为同1列;第二列为通过近似中位秩公式计算得到的Fn(ti);
第三列为lnti,即xi;第四列为yi,即lnln1/1-Fn(ti );第五列填写为第六列填写为第七列填写为第八列填写为第九列填写为
4.如权利要求3所述确定防滑刹车控制装置低温故障分布的方法,其特征在于,所述最小二乘法的回归计算表中:
lnti是对测试时间取对数,在线性方程(9)中用xi表示;xi是双对数坐标系中的横轴;
lnln1/1-Fn(ti)是对1/1-Fn(ti)取两次对数,在线性方程(9)中采用yi表示;
是(14)式中的计算部分,(14)式用于计算相关系数r,根据相关系数r确定各三极管的低温故障数据是否服从威布尔分布;
各三极管的测试数据按照最小二乘法计算,计算后全部填入所述最小二乘法的回归计算表。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610373928.3A CN106054601B (zh) | 2016-05-31 | 2016-05-31 | 确定防滑刹车控制装置低温故障分布的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610373928.3A CN106054601B (zh) | 2016-05-31 | 2016-05-31 | 确定防滑刹车控制装置低温故障分布的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106054601A CN106054601A (zh) | 2016-10-26 |
CN106054601B true CN106054601B (zh) | 2019-05-10 |
Family
ID=57172931
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610373928.3A Expired - Fee Related CN106054601B (zh) | 2016-05-31 | 2016-05-31 | 确定防滑刹车控制装置低温故障分布的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106054601B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106347706B (zh) * | 2016-10-27 | 2018-10-23 | 西安航空制动科技有限公司 | 消除防滑刹车控制装置低温故障的方法 |
CN111523185B (zh) * | 2020-04-22 | 2023-05-23 | 西安航空制动科技有限公司 | 一种飞机刹车控制装置的寿命评估方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5205619A (en) * | 1989-08-30 | 1993-04-27 | Wabco Westinghouse Fahrzeugbremsen Gmbh | Electronic circuit for surveillance of final amplifier and its load |
CN102096036A (zh) * | 2010-12-03 | 2011-06-15 | 华东光电集成器件研究所 | 一种集成三极管阵列电路测试装置 |
CN102928760A (zh) * | 2012-10-31 | 2013-02-13 | 许伟清 | Led老化测试系统及其测试方法 |
CN102981452A (zh) * | 2012-12-28 | 2013-03-20 | 吉林大学 | 数控机床三类功能部件的可靠性建模与可靠性评估方法 |
CN103294052A (zh) * | 2013-05-22 | 2013-09-11 | 西安航空制动科技有限公司 | 采用快速温度变化测试防滑刹车控制盒故障隐患的方法 |
CN103853148A (zh) * | 2013-05-09 | 2014-06-11 | 西安航空制动科技有限公司 | 防滑刹车控制盒低温步进试验的方法 |
-
2016
- 2016-05-31 CN CN201610373928.3A patent/CN106054601B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5205619A (en) * | 1989-08-30 | 1993-04-27 | Wabco Westinghouse Fahrzeugbremsen Gmbh | Electronic circuit for surveillance of final amplifier and its load |
CN102096036A (zh) * | 2010-12-03 | 2011-06-15 | 华东光电集成器件研究所 | 一种集成三极管阵列电路测试装置 |
CN102928760A (zh) * | 2012-10-31 | 2013-02-13 | 许伟清 | Led老化测试系统及其测试方法 |
CN102981452A (zh) * | 2012-12-28 | 2013-03-20 | 吉林大学 | 数控机床三类功能部件的可靠性建模与可靠性评估方法 |
CN103853148A (zh) * | 2013-05-09 | 2014-06-11 | 西安航空制动科技有限公司 | 防滑刹车控制盒低温步进试验的方法 |
CN103294052A (zh) * | 2013-05-22 | 2013-09-11 | 西安航空制动科技有限公司 | 采用快速温度变化测试防滑刹车控制盒故障隐患的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106054601A (zh) | 2016-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102789543B (zh) | 基于多试验信息的可靠性综合评估方法 | |
CN104677997B (zh) | 一种变压器油色谱在线监测差异化预警方法 | |
CN102654539B (zh) | 一种电子式互感器运行状态评价方法 | |
CN102449645B (zh) | 产品检查装置及产品检查方法 | |
CN115982602B (zh) | 一种光伏变压器电故障检测方法 | |
CN110293949A (zh) | 一种高速列车空气制动系统微小故障检测方法 | |
CN105138770B (zh) | 基于间接可靠性特征量的航天产品可靠性仿真评估方法 | |
CN104462808A (zh) | 安全水平位移与水位的滑动可变窗口动态数据拟合方法 | |
CN106054601B (zh) | 确定防滑刹车控制装置低温故障分布的方法 | |
CN106202776A (zh) | 一种运载火箭结构机构产品可靠性评估方法 | |
CN104573315A (zh) | 一种基于状态检修的输变电设备故障率的计算方法 | |
CN105891645B (zh) | 确定防滑刹车控制装置振动故障分布的方法 | |
CN106053984B (zh) | 确定防滑刹车控制装置高温故障分布的方法 | |
CN106855990B (zh) | 核电机组仪表通道测量误差论证方法 | |
CN103344937B (zh) | 智能电能表功耗检测设备及检测方法 | |
CN106442830B (zh) | 变压器油中气体含量告警值的检测方法和系统 | |
CN106644436B (zh) | 一种断路器机械特性的评定方法 | |
CN103217590A (zh) | 一种自动获取大气电场雷电预警特征参数阈值的方法 | |
CN103364669B (zh) | Gis设备运行状态在线检测方法与系统 | |
CN105425182A (zh) | 一种用于故障指示器生产环节的电流精度校准方法 | |
CN110865250A (zh) | 融合电流监测的配电设备状态监测装置及发热检测方法 | |
CN115792583A (zh) | 一种车规级芯片的测试方法、装置、设备及介质 | |
CN103426491B (zh) | 核电厂过程参数可靠性处理方法 | |
CN114001887B (zh) | 一种基于挠度监测的桥梁损伤评定方法 | |
CN115062026A (zh) | 一种用电数据的检测方法、装置、电子设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190510 Termination date: 20210531 |
|
CF01 | Termination of patent right due to non-payment of annual fee |