CN105891645B - 确定防滑刹车控制装置振动故障分布的方法 - Google Patents
确定防滑刹车控制装置振动故障分布的方法 Download PDFInfo
- Publication number
- CN105891645B CN105891645B CN201610373390.6A CN201610373390A CN105891645B CN 105891645 B CN105891645 B CN 105891645B CN 201610373390 A CN201610373390 A CN 201610373390A CN 105891645 B CN105891645 B CN 105891645B
- Authority
- CN
- China
- Prior art keywords
- distribution
- solid
- css
- vibration fault
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
- G01M7/02—Vibration-testing by means of a shake table
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/327—Testing of circuit interrupters, switches or circuit-breakers
- G01R31/3277—Testing of circuit interrupters, switches or circuit-breakers of low voltage devices, e.g. domestic or industrial devices, such as motor protections, relays, rotation switches
- G01R31/3278—Testing of circuit interrupters, switches or circuit-breakers of low voltage devices, e.g. domestic or industrial devices, such as motor protections, relays, rotation switches of relays, solenoids or reed switches
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Electric Properties And Detecting Electric Faults (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)式计算电子设备的故障率,由于按指数分布处理,电子设备的平均故障间隔时间MTBF计算方法为:
MTBF=1/λGS (5)
采用(5)式计算电子设备的平均故障间隔时间MTBF。
在(4)式中电子设备的故障率仅考虑质量系数,未考虑环境和工作电应力的影响。
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,确定电子元器件的振动故障分布测试方案:
第一步,确定元器件振动抽样测试的种类:
以半导体分立元器件作为振动测试样件;在防滑刹车控制装置的元器件中,抽取固态继电器进行振动故障分布测试;
第二步,随机抽取同型号多个防滑刹车控制装置CSS-6固态继电器作为试验样件:
该固态继电器的使用温度范围是:-55℃~125℃;输入电压范围:4.0Vdc~7.0Vdc;输入电流范围:22mA;输出电流:200mA;功耗:260mW;重量为:5g;
第三步,进行CSS-6固态继电器振动故障分布测试:
首先确定CSS-6是该固态继电器振动故障原因;
其次,根据振动故障原因按照HB5830.5E随机振动谱输入三综合试验设备中振动台的控制计算机,将各CSS-6固态继电器按照防滑刹车控制装置的安装要求焊接在电路板上,将电路板固定在振动台的扩展头上;将三综合试验设备的温度箱温度设置为:-55℃恒温,在振动条件下更容易发生振动损伤;振动故障的测试条件为:输入电压5V,输入电流按照防滑刹车控制要求变化,变化范围为:0mA~2 0mA;在通电工作条件下进行测试;进行到各CSS-6固态继电器全部发生振动故障时测试结束;
第四步,确定累积失效概率计算公式:
根据抽试样本数量确定累积失效概率计算公式;
在样本数小于50个的条件下,采用近似中位秩公式计算累积失效概率:
近似中位秩公式:
在所述近似中位秩公式(15)中:
Fn(ti):试验中元器件的累积失效概率;
i:第i个失效的元器件;
0.32:公式中的常数;
0.36:公式中的常数;
n:试验元器件的总数;
第五步,对威布尔分布公式进行线性化处理:
列出威布尔分布数学式:
由得到:
式中:
t为测试时间;m为威布尔分布的形状参数,m的大小表示故障分布的分散程度;t0为尺度参数,缩小和放大分布函数坐标尺度的作用;
对式(7)等号两边取两次以e为底的自然对数线性化处理,得到:
列线性方程:y=bx+a (9)
(9)式是解析几何中常用的线性方程,y为bx+a的函数,a为x=0时方程(9)在y轴上的截距,b为系数;
依据对(9)式对(8)式进行变量转换:
ln是自然对数符号,自然对数是以e为底的对数;
至此,完成了威布尔分布函数的线性化处理,具备在双对数坐标系中对CSS-6固态继电器的振动测试数据进行计算的条件;
第六步,确定截距a、系数b和相关系数r:
按照双对数坐标系的原理,按照线性方程中{x,y}的对应关系,在双对数坐标系中相应列出测试数据i=1,2,3,4;采用最小二乘法确定截距a和系数b:
工程中,的值等于b的值;为威布尔分布的形状参数m的估计值;
x=0时方程(9)式在y轴上截距a的计算公式为:
威布尔分布尺度参数t估计值的计算公式为:
通过公式(14)确定相关系数r:
在威布尔分布双对数坐标系中计算CSS-6固态继电器的振动故障数据,若服从威布尔分布,计算结束;若不服从威布尔分布,则重新在正态分布的双对数坐标系中计算CSS-6固态继电器的振动故障数据;若不服从正态分布,则该CSS-6固态继电器的研制质量存在隐患,应在进行质量改进后重新进行振动故障分布测试和计算;
步骤2,对CSS-6固态继电器进行振动故障分布测试:
按照步骤1确定的测试方案,进行各CSS-6固态继电器的振动故障分布试验,且试验进行到各CSS-6固态继电器全部失效时结束;
所述各CSS-6固态继电器的振动故障分布测试过程是:
Ⅰ将HB5830.5E随机振动谱输入三综合试验设备中振动台的控制计算机;
Ⅱ将各CSS-6固态继电器按照防滑刹车控制装置的安装要求焊接在电路板上,将电路板固定在振动台的扩展头上;
Ⅲ将三综合试验设备的温度箱温度设置为:-55℃恒温,在振动和低温条件下更容易发生振动损伤;
Ⅳ振动故障的测试条件为:输入电压5V,输入电流按照防滑刹车控制要求变化,变化范围为:0mA~2 0mA;在各固态继电器同时通电工作条件下进行测试;直至各CSS-6固态继电器全部发生振动故障时测试结束。
步骤3,对各CSS-6固态继电器的振动故障分布测试数据进行计算
根据步骤2测试得到的各CSS-6固态继电器的振动故障数据,对各测试数据进行计算,确定振动测试数据的故障分布;
根据各CSS-6固态继电器的振动故障分布测试数据列出最小二乘法计算表;
所述最小二乘法的回归分析表中:
lnti是对测试时间取对数,在线性方程(12)中用xi表示;xi是双对数坐标系中的横轴;
lnln1/1-Fn(ti)是对1/1-Fn(ti)取两次对数,在线性方程(12)中采用yi表示;其中Fn(ti)是故障时间的累积失效函数,取两次对数的工程含义是对威布尔分布公式进行线性化处理;
是(14)式中的计算部分,(14)式用于计算相关系数r,根据相关系数r确定各CSS-6固态继电器的振动故障数据是否服从威布尔分布。
所述最小二乘法的回归分析表的填写过程为:
第一列填写ti,具体为各CSS-6固态继电器振动故障的时间数据,按照测试数据由小到大排列的原则在表1中排为同1列;第二列填写按照中位秩次法计算得到的Fn(ti);第三列填写lnti,即xi;第四列填写yi,即lnln1/1-Fn(ti);第五列填写第六列填写第七列填写第八列填写第九列填写
步骤4,确定CSS-6固态继电器的振动故障分布:
根据对CSS-6固态继电器的振动故障测试数据的威布尔分布回归计算结果,若CSS-6固态继电器的振动故障服从威布尔分布,计算结束;
若经计算不服从威布尔分布,则应重新按照正态分布计算,若服从正态分布,计算结束;
若经计算不服从正态分布,则应对该3AX31B固态继电器进行质量改进,改进后重新进行计算,直至确定3AX31B固态继电器服从的概率分布,计算结束。
本发明抽取防滑刹车控制装置中使用的元器件进行振动故障测试,根据测试数据确定元器件的振动故障分布,作为工程中选择故障分布公式的依据,消除现有技术带来的误差和资源消耗。
本发明是针对固态继电器管脚在振动条件下发生疲劳损伤的故障原因进行故障分布函数的研究工作。
所述管脚在振动条件下发生疲劳损伤、在振动和工作条件下发生故障和电子产品在长时间振动条件下发生累积损伤引起的故障均为机械损伤,机械损伤在质量受控的条件下,寿命散差小,国外对金属材料机械性能的控制要求是威布尔分布形状参数m≥8,这时电子产品振动损伤的寿命不服从指数分布,因为指数分布的形状参数m=1。指数分布是概率论与数理统计中的一种分布,其特征是故障率为恒定型。
由于指数分布计算简单,在工程中得到广泛应用,并被国外标准MIL-HDBK-217F、国内GJB299、GJB899采纳。
本发明技术方案的实施途径:由于振动原因发生故障,制定该防滑刹车控制装置振动测试方案进行测试,对测试取得的振动故障数据进行回归分析,得到该防滑刹车控制装置的振动故障分布。
为证明本发明的效果,本发明对一组固态继电器分别采用威布尔分布和指数分布计算在HB5830.5E振动谱规定的随机振动条件下测试1000h的瞬时故障率数据,计算结果相差5.0939×107倍,证明本发明提出的根据使用要求确定故障分布的方法,避免由于计算误差带来的损失。证明本发明效果的计算过程为::
第一步,采用指数分布计算该固态继电器测试1000h的故障率
当威布尔分布的位置参数为零,形状参数为1时就等于指数分布。固态继电器使用到1000h的故障率为:
第二步,采用威布尔分布计算该固态继电器测试1000h的故障率
本实施例通过试验已经证明,该固态继电器的形状参数m=3.569,威布尔分布的位置参数为零,固态继电器使用到1000h的故障率为:
第三步,瞬时故障率计算对比
产品的故障率分有瞬时故障率和恒定故障率之分,在恒定故障率的条件下,故障率为常数。本发明故障率的大小是时间的函数,因此在进行故障率比较时,在相同的时间点上进行故障率大小的比较,是瞬时故障率。
λ(t威布尔)/λ(t指数)=1.8985×10-9/0.3727×10-16=5.0939×107
即在使用1000h的条件下,采用威布尔分布计算得到的振动故障率是指数分布的5.0939×107倍,采用指数分布得到的故障率与产品实际的故障率相差大,不宜作为计算方法。
与现有技术相比,本发明取得的有益效果为:
采用本发明技术进行防滑刹车控制装置的振动故障分布参数计算,经计算形状参数等于3.569是大于1的数字,不服从指数分布。确定了在对防滑刹车控制装置的振动故障分布参数进行分析时,采用威布尔分布预计防滑刹车控制装置的可靠性指标,能够得到平均故障前时间。
本发明和现有标准实施效果的综合比较:
1、现有国内外标准采用指数分布进行电子产品的可靠性预计,指数分布的故障率为常数,不适作为计算电子产品振动故障的方法。
2、采用本发明确定的威布尔分布进行可靠性预计工作,能够提高电子产品的可靠性预计值与实际的可靠性指标的符合程度;
3、本发明确定采用威布尔分布计算防滑刹车控制装置的振动故障分布,证明防滑刹车控制装置的故障不是偶然故障,而是耗损故障。在可靠性工程领域产品的故障分为:早期故障、偶然故障、耗损故障,早期故障指使用初期出现的故障,偶然故障指偶然因素引起的故障,耗损故障指磨损、老化、腐蚀、疲劳等物理化学原因引起的故障。现有标准将耗损故障作为偶然故障处理,故障原因、故障机理分析不准确产生误差,本发明的效果是纠正了这种差错。
具体实施方式
实施例
本实施例是确定某型飞机防滑刹车控制装置振动故障分布的方法。
该防滑刹车控制装置的振动故障由电子元器件的振动故障引起,因此,确定该防滑刹车控制装置的振动故障分布就是确定组成该防滑刹车控制装置电子元器件的振动故障分布。确定所有电子元器件振动故障分布的测试时间长,费用高,在工程中也没有将全部电子元器件进行故障分布测试的先例。本实施例依据数理统计原理,采用抽样的方法确定防滑刹车控制装置的振动故障分布,就是用子样的故障分布表示母体的故障分布,子样指进行振动测试的元器件,母体指防滑刹车控制装置。
常用的故障分布有:
1)指数分布,故障率为恒定型,这是美国标准MIL-STD-217F、中国标准GJB299采用的故障分布;指数分布是概率论与数理统计中的一种常用分布。
2)正态分布,故障率为时间的增函数,适用于老化、腐蚀、磨损、疲劳耗损原因引起的故障,正态分布有方差和均值两个参数。
3)威布尔分布,故障率为时间的增函数,适用于老化、腐蚀、磨损、疲劳耗损原因引起的故障,威布尔分布有形状参数m、尺度参数t0、位置参数γ三个参数。
本实施例根据电子元器件的振动故障为机械故障,首先假设其故障服从威布尔分布,并进行电子元器件的振动故障分布测试,根据测试数据进行线性回归分析,若线性回归分析结果服从威布尔分布,线性回归分析结束。若线性回归分析不服从威布尔分布,则重新进行线性回归分析,直至得到所服从的故障分布。回归分析前应对威布尔分布公式进行线性化处理,线性化处理后的分析称为线性回归分析。
本实施例采用的“分布”是一个来自概率论与数理统计的常用术语,含义是任何一个随机变量都用概率分布的形式表示,说明示例:同一类产品早期发生的故障少,达到长寿命的也少,寿命集中在平均值附近的居多,使用时间是自变量,函数值的大小是因变量,统称为故障的密度函数。常用耗损故障的概率分布有正态分布、威布尔分布,分布的物理含义是产品的寿命不是确定值,而是随机变量,在该产品失效前只能预计该产品可能的失效时间范围,而不能确定该产品准确的失效时间。
失效是可靠性工程的翻译词,原文是:failure,在中国的工程领域,失效也叫故障,二个名词等同。
本实施例的具体过程是:
步骤1,确定电子元器件的振动故障分布测试方案
第一步,确定元器件振动抽样测试的种类
现有元器件分为半导体分立元器件、半导体集成元器件两大类,分立、集成的制造方式根据设计用途确定,和元器件的振动故障原因无关,在振动条件作用下,分立、集成元器件都会发生故障,故障原因同为振动引起的元器件管脚损伤、元器件性能飘移,严重时发生故障,这是半导体材料的物理属性。所以,抽样取分立元器件还是集成器件不影响故障分布的测试效果。
本发明选择半导体分立元器件作为振动测试样件。在防滑刹车控制装置的元器件中,在振动条件下出现故障的是固态继电器,当固态继电器出现振动故障时,防滑刹车控制装置不能正常工作,抽取固态继电器进行振动故障分布测试。
按照随机抽样原则抽取固态继电器进行振动故障分布测试。
第二步,随机抽取同型号4个CSS-6固态继电器作为测试样件
固态继电器由可控硅组成,是可控硅开关功能的一种应用。
本发明随机抽取4个CSS-6固态继电器作为振动故障分布的测试样件,该固态继电器的使用温度范围是:-55℃~125℃;输入电压范围:4.0Vdc~7.0Vdc;输入电流范围:22mA;输出电流:200mA;功耗:260mW;重量为:5g。根据线性回归分析原理,在平面上二点确定一条直线,工程意义是最少抽取二个CSS-6固态继电器,在原理上二个CSS-6固态继电器的振动故障数据就可以确定该CSS-6固态继电器的振动故障分布。但是,由于元器件制造过程中不可避免产生偏差,为了提高振动故障测试的精度,将抽取样件的数量增加一倍,共抽取四个CSS-6固态继电器进行振动故障分布测试。
CSS-6是该固态继电器的型号。
第三步,进行CSS-6固态继电器振动故障分布测试
本步骤的测试方案属于子样的测试方案,因为将所有元器件进行振动故障测试费用高,时间长。
首先确定CSS-6该固态继电器振动故障原因,包括:
1)管脚在振动条件下发生疲劳损伤;
2)在振动和工作条件下发生故障;
3)CSS-6固态继电器在长时间振动条件下发生累积损伤。
其次,根据振动故障原因按照HB5830.5E随机振动谱输入三综合试验设备中振动台的控制计算机,将4个CSS-6固态继电器按照防滑刹车控制装置的安装要求焊接在电路板上,将电路板固定在振动台的扩展头上。将三综合试验设备的温度箱温度设置为:-55℃恒温,在低温和振动条件下更容易发生振动损伤;振动故障的测试条件为:输入电压5V,输入电流按照防滑刹车控制要求变化,变化范围为:0mA~2 0mA;在通电工作条件下进行测试。进行到4个CSS-6固态继电器全部发生振动故障时测试结束。
第四步,确定累积失效概率计算公式
根据抽试样本数量确定累积失效概率计算公式。
在样本数小于50个的条件下,采用近似中位秩公式计算累积失效概率:
近似中位秩公式:
在所述近似中位秩公式(15)中:
Fn(ti):试验中元器件的累积失效概率;
i:第i个失效的元器件;本实施例中,i=1~4;
0.32:公式中的常数;
0.36:公式中的常数。
n:试验元器件的总数,在本实施例中n=4。
第五步,对威布尔分布公式进行线性化处理
因为本发明列出的公式(7)为非线性,该公式是概率论与数理统计中的公式,不能直接进行回归分析,因此需要进行线性化处理。线性化处理的对象是威布尔公式。
经线性化后的威布尔公式,用于进行CSS-6固态继电器的振动故障测试数据的回归分析。CSS-6固态继电器是从该防滑刹车控制装置抽取进行振动故障分布测试的元器件。
本发明确定在分析过程中,当对CSS-6固态继电器的振动故障测试数据回归分析相对于威布尔分布的线性相关系数r大于等于0.9时,这些数据与威布尔分布的线性相关好,并且线性相关系数r越大测试数据与威布尔分布的拟合程度越好。
本发明确定,当该CSS-6固态继电器的振动故障分布测试数据回归分析相对于威布尔分布线性相关系数r小于0.9时,就认为该CSS-6固态继电器的振动故障测试数据与威布尔分布拟合不好,对该CSS-6固态继电器的振动故障试验数据采用正态分布重新计算。
若采用正态分布对该CSS-6固态继电器的振动故障测试数据回归分析结果,线性相关系数也小于0.9,则该批产品质量一致性不好,应进行质量改进,改进后重新进行该CSS-6固态继电器的振动测试和故障数据的计算工作。
列出威布尔分布数学式:
由得到:
式中:
t为时间;m为威布尔分布的形状参数,m的大小表示故障分布的分散程度;t0为尺度参数,缩小和放大分布函数坐标尺度的作用;
对式(7)两边取两次自然对数进行线性化处理。两边取对数是工程数学中线性化处理的方法。自然对数是以e为底的对数,得到(8)式:
列线性方程:y=bx+a (9)
(9)式是解析几何中常用的线性方程,y为bx+a的函数,a为x=0时方程在y轴上的截距,b为系数。
依据(9)式对(8)式进行变量转换:
ln是自然对数符号,自然对数是以e为底的对数。
至此,完成了分布函数的线性化处理,具备在双对数坐标系中对CSS-6固态继电器的振动测试数据进行计算的条件。
第六步,确定截距a、系数b和相关系数r;
在完成(10)式的变量转化后,线性方程(9)中a、b的物理意义是使固态继电器振动测试数据在威布尔分布双对数坐标系中的数据点与线性方程(9)所确定的回归直线距离最小,计算精度最高。
因为本发明所用最小二乘法是在假设固态继电器的振动测试数据服从威布尔分布的条件下完成计算过程的,最小二乘法具有判别测试数据是否服从威布尔分布的功能,当相关系数r大于等于0.9时,固态继电器的振动测试数据与回归直线的符合性好,固态继电器的振动测试数据服从威布尔分布;当相关系数r小于0.9时,固态继电器的振动测试数据与回归直线的符合性不好,固态继电器的振动测试数据不服从威布尔分布;在固态继电器的振动测试数据不服从威布尔分布的条件下,应采用正态分布重新计算固态继电器的振动测试数据。
按照双对数坐标系的原理,按照线性方程中{x,y}的对应关系,在双对数坐标系中相应列出试验数据i=1,2,3,4。采用工程数学中的最小二乘法确定截距a和系数br,根据式(10),采用工程数学中的最小二乘法确定的计算公式为:
工程中,的值等于b的值。为威布尔分布的形状参数m的估计值。
x=0时方程(9)式在y轴上截距a的计算公式为:
威布尔分布的尺度参数t的计算公式为:
(13)式中为t的估计值,用表示针对尺度参数t进行计算后的估计值。e是自然对数中的e,在最小二乘法计算过程中,对(7)式两边取自然对数的方法进行线性化处理。
通过公式(14)确定相关系数r:
至此,所确定的该防滑刹车控制装置电子元器件振动故障测试方案包括下列内容:
1)抽检测试的元器件种类为CSS-6固态继电器;
2)抽检测试的CSS-6固态继电器数量为4个;
3)测试结束的条件是CSS-6固态继电器全部发生振动故障;
4)选取累积失效计算公式为近似中位秩公式;
5)在威布尔分布双对数坐标系中计算CSS-6固态继电器的振动故障数据,若服从威布尔分布,计算结束;若不服从威布尔分布,则重新在正态分布的双对数坐标系中计算CSS-6固态继电器的振动故障数据;若不服从正态分布,则该CSS-6固态继电器的研制质量存在隐患,应在进行质量改进后重新进行振动故障分布测试和计算。
步骤2,对CSS-6固态继电器进行振动故障分布测试
按照本发明步骤1确定的测试方案,进行4个CSS-6固态继电器的振动故障分布测试,且测试进行到4个CSS-6固态继电器全部失效时结束。
4个CSS-6固态继电器的振动故障分布测试过程是:
1)将HB5830.5E随机振动谱输入三综合试验设备中振动台的控制计算机;
2)将4个CSS-6固态继电器按照防滑刹车控制装置的安装要求焊接在电路板上,将电路板固定在振动台的扩展头上;
3)将三综合试验设备的温度箱温度设置为:-55℃恒温,在振动和低温条件下更容易发生振动损伤;
4)振动故障的测试条件为:输入电压5V,输入电流按照防滑刹车控制要求变化,变化范围为:0mA~20mA;在通电工作条件下进行测试;
5)进行到4个CSS-6固态继电器全部发生振动故障时测试结束。
经测试,4个CSS-6固态继电器振动测试发生故障的时间分别为:21373h,34196h,34196h,38471h,故障原因为管脚损伤。
步骤3,对4个CSS-6固态继电器的振动故障分布测试数据进行计算
根据本发明步骤2测试得到的4个CSS-6固态继电器的振动故障数据,本步骤对这4个测试数据进行计算,确定振动测试数据的故障分布。
根据4个CSS-6固态继电器的振动故障分布测试数据列出最小二乘法计算表。
为了方便计算,本发明列出最小二乘法的回归分析表见表1。填写过程为:
第一列填写ti,具体为4个CSS-6固态继电器的振动故障时间数据,按照测试数据由小到大的原则在表1中排为同1列;;第二列填写按照中位秩次法计算的Fn(ti);第三列填写lnti,即xi;第四列填写yi,即lnln1/1-Fn(ti);第五列填写第六列填写第七列填写第八列填写第九列填写
其中:
lnti是对测试时间取对数,在线性方程(12)中用xi表示;xi是双对数坐标系中的横轴;
lnln1/1-Fn(ti)是对1/1-Fn(ti)取两次对数,在线性方程(12)中采用yi表示;其中Fn(ti)是故障时间的累积失效函数,取两次对数的工程含义是对威布尔分布公式进行线性化处理;
是公式(14)中的计算部分,(14)式用于计算相关系数r,根据相关系数r确定4个CSS-6固态继电器的振动故障数据是否服从威布尔分布。在概率论与数理统计中,测试数据与概率分布的符合性,常用术语“服从”表示。
4个固态继电器的测试数据按照最小二乘法计算方法,计算后全部填入表1。
表1 HB5830.5E振动谱,-55℃,4个数据威布尔分布计算,时间ti单位:h
根据表1数据,由(14)式计算相关系数r的数值:
相关系数
相关系数r等于0.9426,CSS-6固态继电器的振动故障服从威布尔分布。继续评估威布尔分布参数。
由(11)式得形状参数m的估计值为:
b为公式(9)式中的系数。
工程中,的值等于b的值。为威布尔分布的形状参数m的估计值。
由(12)式得x=0时方程(9)在y轴上的截距a的计算数值为::
由(13)式得尺度参数t的估计值为:
本步骤达到了下列目的:
1)利用CSS-6固态继电器的振动故障测试数据,计算与威布尔分布的线性相关系数为0.9426,服从威布尔分布;
2)得到了CSS-6固态继电器振动故障的威布尔分布参数:形状参数m=3.569,尺度参数t0=2.6828×1016。
步骤4,确定CSS-6固态继电器的振动故障分布
在可靠性工程中,常用概率分布表示故障发生的多少和时间的关系,分布是概率论与数理统计领域的常用语,概率论与数理统计是可靠性工程的基础之一。举例而言,在CSS-6固态继电器的振动故障服从威布尔分布的条件下,CSS-6固态继电器的振动寿命短、振动寿命长的数量少,而在平均寿命时间的CSS-6固态继电器数量多,寿命的长短和时间之间存在随机函数关系,这种随机函数用概率分布的公式表示。
根据对CSS-6固态继电器的振动故障测试数据的威布尔分布回归计算结果,若CSS-6固态继电器的振动故障服从威布尔分布,计算结束。
若经计算不服从威布尔分布,则应重新按照正态分布计算,若服从正态分布,计算结束。
若经计算不服从正态分布,则应对该CSS-6固态继电器进行质量改进,改进后重新进行计算,直至确定CSS-6固态继电器服从的概率分布,计算结束。
在电子元器件研制质量稳定的条件下,根据抽样测试的方法,用子样的故障分布表示母体的故障分布。不采用将所有元器件都进行振动故障测试的方法。
本发明确定,采用威布尔分布公式对CSS-6固态继电器的振动故障测试数据进行回归分析的计算结果,CSS-6固态继电器的振动故障与威布尔分布的相关系数为0.9426,CSS-6固态继电器的振动故障服从威布尔分布。
本发明确定,CSS-6固态继电器的振动故障的形状参数m=3.569,故障率为递增型,而不是恒定型。所以,CSS-6固态继电器的振动故障不服从指数分布,指数分布的形状参数m=1,即故障率为恒定型。
本发明确定,对于由振动故障引起的电子元器件故障,确定为服从威布尔分布。根据本发明结果,按照威布尔分布进行电子产品的故障分析。
Claims (4)
1.一种确定防滑刹车控制装置振动故障分布的方法,其特征在于,具体过程是:
步骤1,确定电子元器件的振动故障分布测试方案:
第一步,确定元器件振动抽样测试的种类:
以半导体分立元器件作为振动测试样件;在防滑刹车控制装置的元器件中,抽取固态继电器进行振动故障分布测试;
第二步,随机抽取同型号多个防滑刹车控制装置CSS-6固态继电器作为试验样件:该固态继电器的使用温度范围是:-55℃~125℃;输入电压范围:4.0Vdc~7.0Vdc;
输入电流:22mA;输出电流:200mA;功耗:260mW;重量为:5g;
第三步,进行CSS-6固态继电器振动故障分布测试:
首先确定CSS-6固态继电器振动故障原因;
其次,根据振动故障原因将HB5830.5E随机振动谱输入三综合试验设备中振动台的控制计算机,将各CSS-6固态继电器按照防滑刹车控制装置的安装要求焊接在电路板上,将电路板固定在振动台的扩展头上;将三综合试验设备的温度箱温度设置为:-55℃恒温,在振动条件下更容易发生振动损伤;振动故障的测试条件为:输入电压5V,输入电流按照防滑刹车控制要求变化,变化范围为:0mA~2 0mA;在通电工作条件下进行测试;进行到各CSS-6固态继电器全部发生振动故障时测试结束;
第四步,确定累积失效概率计算公式:
根据抽试样本数量确定累积失效概率计算公式;
在样本数小于50个的条件下,采用近似中位秩公式计算累积失效概率:
近似中位秩公式:
在所述近似中位秩公式(15)中:
Fn(ti):试验中元器件的累积失效概率;
i:第i个失效的元器件;
0.32:公式中的常数;
0.36:公式中的常数;
n:试验元器件的总数;
ti:第i个失效的元器件的测试时间;
第五步,对威布尔分布公式进行线性化处理:
列出威布尔分布数学式:
由得到:
式中:
t为测试时间;m为威布尔分布的形状参数,m的大小表示故障分布的分散程度;
t0为尺度参数,缩小和放大分布函数坐标尺度的作用;
对式(7)等号两边取两次以e为底的自然对数线性化处理,得到:
列线性方程:y=bx+a (9)
(9)式是解析几何中常用的线性方程,y为bx+a的函数,a为x=0时方程(9)在y轴上的截距,b为系数;
依据对(9)式对(8)式进行变量转换:
ln是自然对数符号,自然对数是以e为底的对数;
至此,完成了威布尔分布函数的线性化处理,具备在双对数坐标系中对CSS-6固态继电器的振动测试数据进行计算的条件;
第六步,确定截距a、系数b和相关系数r:
按照双对数坐标系的原理,按照线性方程中{x,y}的对应关系,在双对数坐标系中相应列出测试数据采用最小二乘法确定截距a和系数b:
工程中,的值等于b的值;为威布尔分布的形状参数m的估计值;
x=0时方程(9)式在y轴上截距a的计算公式为:
威布尔分布尺度参数t0估计值的计算公式为:
通过公式(14)确定相关系数r:
在威布尔分布双对数坐标系中计算CSS-6固态继电器的振动故障数据,若服从威布尔分布,计算结束;若不服从威布尔分布,则重新在正态分布的双对数坐标系中计算CSS-6固态继电器的振动故障数据;若不服从正态分布,则该CSS-6固态继电器的研制质量存在隐患,应在进行质量改进后重新进行振动故障分布测试和计算;
步骤2,对CSS-6固态继电器进行振动故障分布测试:
按照步骤1确定的测试方案,进行各CSS-6固态继电器的振动故障分布试验,且试验进行到各CSS-6固态继电器全部失效时结束;
步骤3,对各CSS-6固态继电器的振动故障分布测试数据进行计算
根据步骤2测试得到的各CSS-6固态继电器的振动故障数据,对各测试数据进行计算,确定振动测试数据的故障分布;
根据各CSS-6固态继电器的振动故障分布测试数据列出最小二乘法计算表;
步骤4,确定CSS-6固态继电器的振动故障分布:
根据对CSS-6固态继电器的振动故障测试数据的威布尔分布回归计算结果,若CSS-6固态继电器的振动故障服从威布尔分布,计算结束;
若经计算不服从威布尔分布,则应重新按照正态分布计算,若服从正态分布,计算结束;
若经计算不服从正态分布,则应对该CSS-6固态继电器进行质量改进,改进后重新进行计算,直至确定CSS-6固态继电器服从的概率分布,计算结束。
2.如权利要求1所述一种确定防滑刹车控制装置振动故障分布的方法,其特征在于,所述各CSS-6固态继电器的振动故障分布测试过程是:
Ⅰ将HB5830.5E随机振动谱输入三综合试验设备中振动台的控制计算机;
Ⅱ将各CSS-6固态继电器按照防滑刹车控制装置的安装要求焊接在电路板上,将电路板固定在振动台的扩展头上;
Ⅲ将三综合试验设备的温度箱温度设置为:-55℃恒温,在振动和低温条件下更容易发生振动损伤;
Ⅳ振动故障的测试条件为:输入电压5V,输入电流按照防滑刹车控制要求变化,变化范围为:0mA~2 0mA;在各固态继电器同时通电工作条件下进行测试;直至各CSS-6固态继电器全部发生振动故障时测试结束。
3.如权利要求1所述一种确定防滑刹车控制装置振动故障分布的方法,其特征在于,所述最小二乘法的回归分析表的填写过程为:
第一列填写ti,具体为各CSS-6固态继电器振动故障的时间数据,按照测试数据由小到大排列的原则在表1中排为同1列;第二列填写按照中位秩次法计算得到的Fn(ti);第三列填写lnti,即xi;第四列填写yi,即lnln1/1-Fn(ti);第五列填写第六列填写第七列填写第八列填写第九列填写
4.如权利要求3所述一种确定防滑刹车控制装置振动故障分布的方法,其特征在于,所述最小二乘法的回归分析表中:
lnti是对测试时间取对数,在线性方程(12)中用xi表示;xi是双对数坐标系中的横轴;
lnln1/1-Fn(ti)是对1/1-Fn(ti)取两次对数,在线性方程(12)中采用yi表示;其中Fn(ti)是故障时间的累积失效函数,取两次对数的工程含义是对威布尔分布公式进行线性化处理;
是(14)式中的计算部分,(14)式用于计算相关系数r,根据相关系数r确定各CSS-6固态继电器的振动故障数据是否服从威布尔分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610373390.6A CN105891645B (zh) | 2016-05-31 | 2016-05-31 | 确定防滑刹车控制装置振动故障分布的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610373390.6A CN105891645B (zh) | 2016-05-31 | 2016-05-31 | 确定防滑刹车控制装置振动故障分布的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105891645A CN105891645A (zh) | 2016-08-24 |
CN105891645B true CN105891645B (zh) | 2018-10-09 |
Family
ID=56710567
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610373390.6A Expired - Fee Related CN105891645B (zh) | 2016-05-31 | 2016-05-31 | 确定防滑刹车控制装置振动故障分布的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105891645B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106483407B (zh) * | 2016-10-08 | 2019-03-19 | 西安航空制动科技有限公司 | 确定防滑刹车控制装置振动薄弱环节的方法 |
CN108459948B (zh) * | 2018-03-26 | 2021-03-09 | 华北电力大学(保定) | 系统可靠性评估中失效数据分布类型的确定方法 |
CN111413953A (zh) * | 2020-04-10 | 2020-07-14 | 西安航空制动科技有限公司 | 确定刹车控制装置振动测试参数的方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102426451A (zh) * | 2011-12-21 | 2012-04-25 | 西安航空制动科技有限公司 | 飞机防滑刹车控制盒的加速可靠性试验方法 |
CN102629300A (zh) * | 2012-03-15 | 2012-08-08 | 北京航空航天大学 | 一种基于灰色预测模型的步进应力加速退化数据评估方法 |
CN102945316A (zh) * | 2012-10-25 | 2013-02-27 | 华北电力大学 | 一种考虑协变量的继电保护设备失效率计算方法 |
CN103020438A (zh) * | 2012-11-29 | 2013-04-03 | 中国人民解放军军械工程学院 | 一种基于混合威布尔分布的航空发动机可靠性监测方法 |
CN103065052A (zh) * | 2013-01-07 | 2013-04-24 | 河南科技大学 | 一种测定机械产品理论寿命的方法 |
DE102012005955A1 (de) * | 2012-03-23 | 2013-09-26 | Robert Bosch Gmbh | Verfahren zur Bestimmung von Ausfallwahrscheinlichkeiten eines Bauteils |
CN103413048A (zh) * | 2013-08-15 | 2013-11-27 | 深圳供电局有限公司 | 基于三参数威布尔分布确定电网设备最佳退役时间的方法 |
CN103412986A (zh) * | 2013-07-23 | 2013-11-27 | 中国人民解放军海军航空工程学院 | 一种基于无故障数据分析的飞机定检内容优化方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012233814A (ja) * | 2011-05-06 | 2012-11-29 | Microsignal Kk | 振動試験制御装置 |
-
2016
- 2016-05-31 CN CN201610373390.6A patent/CN105891645B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102426451A (zh) * | 2011-12-21 | 2012-04-25 | 西安航空制动科技有限公司 | 飞机防滑刹车控制盒的加速可靠性试验方法 |
CN102629300A (zh) * | 2012-03-15 | 2012-08-08 | 北京航空航天大学 | 一种基于灰色预测模型的步进应力加速退化数据评估方法 |
DE102012005955A1 (de) * | 2012-03-23 | 2013-09-26 | Robert Bosch Gmbh | Verfahren zur Bestimmung von Ausfallwahrscheinlichkeiten eines Bauteils |
CN102945316A (zh) * | 2012-10-25 | 2013-02-27 | 华北电力大学 | 一种考虑协变量的继电保护设备失效率计算方法 |
CN103020438A (zh) * | 2012-11-29 | 2013-04-03 | 中国人民解放军军械工程学院 | 一种基于混合威布尔分布的航空发动机可靠性监测方法 |
CN103065052A (zh) * | 2013-01-07 | 2013-04-24 | 河南科技大学 | 一种测定机械产品理论寿命的方法 |
CN103412986A (zh) * | 2013-07-23 | 2013-11-27 | 中国人民解放军海军航空工程学院 | 一种基于无故障数据分析的飞机定检内容优化方法 |
CN103413048A (zh) * | 2013-08-15 | 2013-11-27 | 深圳供电局有限公司 | 基于三参数威布尔分布确定电网设备最佳退役时间的方法 |
Non-Patent Citations (2)
Title |
---|
基于Weibull分布的双馈风力发电机轴承寿命预测;王瑞军 等;《兰州交通大学学报》;20151231;第34卷(第6期);第117-121页 * |
基于设备状态振动特征的比例故障率模型可靠性评估;丁锋 等;《机型工程学报》;20091231;第45卷(第12期);第89-94页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105891645A (zh) | 2016-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106874582B (zh) | 一种电主轴加速寿命试验时间设计方法 | |
CN105891645B (zh) | 确定防滑刹车控制装置振动故障分布的方法 | |
CN109948860A (zh) | 一种机械系统剩余寿命预测方法及系统 | |
CN104677997B (zh) | 一种变压器油色谱在线监测差异化预警方法 | |
CN106919127B (zh) | 一种基于软件虚拟技术的物料料位检测方法 | |
CN102167259B (zh) | 起重机械载荷质量的测量方法 | |
CN116399402B (zh) | 一种用于生态环境监测的无线传感器的故障预警系统 | |
CN114295176A (zh) | 一种储油柜油位测算方法及装置 | |
CN115063017A (zh) | 一种中小跨径桥梁结构监测和评估系统及方法 | |
CN116934303A (zh) | 新能源汽车电池封装用耐温湿聚氨酯粘合剂性能检测系统 | |
CN115792364A (zh) | 一种基于待机状态分析的功耗计算方法 | |
CN110160615B (zh) | 车载称重方法、装置及系统 | |
CN106053984B (zh) | 确定防滑刹车控制装置高温故障分布的方法 | |
CN110488188B (zh) | 基于动态阈值的机组三维健康量化评价方法 | |
CN109613372B (zh) | 一种基于多元电网数据库的电网故障诊断方法 | |
CN102722636A (zh) | 一种针对航空电子系统面临电磁环境复杂度的量化评估方法 | |
CN116522096B (zh) | 基于动作捕捉的三维数字孪生内容智能制作方法 | |
CN117072520A (zh) | 一种基于plc控制器的液压元件检测系统 | |
CN106646326A (zh) | 一种电能计量装置智能化监测方法 | |
CN106598026A (zh) | 输变电设备状态检修时间检测方法和装置 | |
CN115792583A (zh) | 一种车规级芯片的测试方法、装置、设备及介质 | |
CN115062026A (zh) | 一种用电数据的检测方法、装置、电子设备及存储介质 | |
CN114858250A (zh) | 一种智慧水表的故障检测方法、系统及可读存储介质 | |
CN111523185B (zh) | 一种飞机刹车控制装置的寿命评估方法 | |
KR20110107481A (ko) | 피로도 센서기반 임베디드 정비관리시스템 |
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: 20181009 Termination date: 20210531 |
|
CF01 | Termination of patent right due to non-payment of annual fee |