资源描述：
《线缆产品失效概率威布尔分布函数回归分析的计算机程序》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、线缆产品失效概率威布尔分布函数回归分析的计算机程序书京电线厂刘烨,、,摘要在电线电缆可靠性设计工作中为改变以往用人工速度慢精度低的问题提供了一种在,,已知失效数据的情况下建立失效函数分布形式的程序使线缆产品的失效分析能高、、。速精确多功能地进行主题词电线电缆产品可靠性程序设计电线电,缆产品在使用一段时间后不可,。基本思想避免发生失效例如电压击穿及机械损伤等虽然在某一产品上发生失效的事是不可预测:威布尔方程的形式为,F(=“的但是一批同型号产品在使用期内发生失`,`一p,卜(气与j效的情况却是有规律的这在概率学上称失。,式中F(分产品失
2、效率效函数分布通过大量失效数据的分析整理—,可,万形状参数衡量寿命离散程度以描绘产品失效规律从而准确判断该产—刀-,。尺度参数衡量寿命水平的单品在外界影响下的失效情况一在线,位长度缆产品失效分析中威布尔失效函,,数,位置参数即在此以前不会失分布用途最广而且线缆产品遵循威布尔—失效形式。效,:无论是电线电缆的平均寿命T与失效概此外尚有如下失效数据(1)在统计:率,云意义上表示平均失效程度的均值刀(O为F(T)电缆护套热寿命与失效概率=犷+”+,E“,r`F(约或是电缆绝缘层击穿电压V与失效概(矗),。率F(V)均属威布尔分布在线缆产品失
3、(2)表示各统计数据对均值偏离程度的方差,重要工作效分析中就是要建立失效函数的V(咨:)为形式及,各个参数值通过失效函数分布的建F(r`一·2立及,=)t各个参数的确定就可测定产品在各环叫介劲叶l(翻子。境应力情况下的工作寿命以前这种工作是(3)表示在外界应力影响下平均可靠度的可,:用手工将失效数据描在概率纸上然后凭目靠性刀(分为。,二(`)二。xp何力估算各参数值此种方法费时费力精度卜仁二牛、L、刀/,。低速度慢无法满足工作需要如果利用计,将威布尔方程变形并对其两端取二次,则可算机运行预定程序在短时间内完成上,自然对数则,述工作而且
4、可以精确完成失效数据分析处__1_「/云一犷、Ml。不一一下布万子=U工PI)l理的多项工作-J一刃气`夕L气、一,rIJ—/木文介绍以威布尔失效函数分布回归方,n,n=M`n(`一:(`万卜“。)。ha程建立的计算机程序一54一艺,。(`一刀1t之I二Jnl刀j一下万牙不户解网-L一卫又石)台司·““了二In〔z一咐N一、〔`n(`一r,〕B=叮In刀一〔名I}n(t一r)〕“N一’}一`,~则可将式(1)化为如下线性方程形式,二:一n`,“`丁丁刃万“Y二卫.X+B益{〔;」,~’二一,,一“,n““它在夕座标上是一条直线其斜率为
5、M!拭h-a」资}。截距为B·艺`n’“,,:亡万以此为例将式(1)变形得{[益〕`=。’`“·犷、、了一,一.叮、/乙n、了夕、.F一1飞’/M/1l3f可〕一`n闷孟,击L石丁二苏丈万」/了’/“二一=。二一`,pM·101·’卜下蛋丽{音青()以上变换后的方程均可表示为如下的线{「俞,,、,__:1~,性方程形式乙In(云一r)InIn、一下币不`_、`:=,,+,=,,一YAXB(I123……)刀〔In(卜犷)一〕“,:巡其中式(1)为_.、r__11,,_。一乙`n一-一`1(,犷)万`n`n四了,二``n了,二Inr)L厕
6、万」且(卜不万一艺Zn`一犷一1丁风l()NAl二万Bl=一In勺ǐJlesesese了!.,1艺In(卜r):式(2)为.,—一丽一—,/1、1,wJYZ=云2二上n而:人火丁二万可)以上三式均可表示为二lA:=州/、B:=r万f(灼犷二Z:f(万)式(3)为=3,r、l勺(M「1〕规f)一了“二’nL丁丁夕J匀,则万二,可如果数据为威布尔分布f(灼:二Z。:A3二丫/材B3=一甲甘/护及f(M)必有交点对交点讨论如下,、(l),:由于平面直线只有两个自由度欲求M如果数据服从威布尔分布f(灼与2,、,。;勺丫三个参数值必须将三式联立
7、了(万)有效且有一个交点(,,2)如果有三个交点数据必然有误程序原理,;(3)曲线必为单调上升或下降,r>=且X+B,(4)定义域为M>o叭对于回归方程Y其系数分。别为:(5)一般M值取小于7,。A=圣州丝对一般威布尔分布M值很少超过7对一里艺些二“一三(些)“/万于超过7的威布尔分布则分散过强,己经无勘=。,,B￡州N一A艺对N统计意义因此设M<7的范围程序、、,::二,:二,二:十Z,将M云刀分别代入可得其表达式为万O万7万(MM)/2然后取一一··。`,,`,、二:In(卜)`n`n:二。)犷值取M值并与M值相比如果河>Mf〔卷
8、〕,。:则增大M值此时程序得以循环当万与一:,65PRINTTAB(Jo)”pLEASEwAITING口值M之差大于迭代精度时说明此时数据不70l=0:Z:,+么1:吕,,MM二7M二(M)/服从威布尔分布应另选模型程序
