CN107478455A - 一种适用于威布尔分布型产品的定时截尾可靠性试验方法 - Google Patents
一种适用于威布尔分布型产品的定时截尾可靠性试验方法 Download PDFInfo
- Publication number
- CN107478455A CN107478455A CN201710777965.5A CN201710777965A CN107478455A CN 107478455 A CN107478455 A CN 107478455A CN 201710777965 A CN201710777965 A CN 201710777965A CN 107478455 A CN107478455 A CN 107478455A
- Authority
- CN
- China
- Prior art keywords
- weibull distribution
- test
- time
- obtaining
- parameter
- 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
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M99/00—Subject matter not provided for in other groups of this subclass
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种适用于威布尔分布型产品的定时截尾可靠性试验方法,其包括以下步骤:S1、根据历史数据获取威布尔分布的形状参数m;确定定时截尾试验方案的基本参数;S2、根据威布尔分布的形状参数m和定时截尾试验方案的基本参数确定试验方案的允许故障数c;S3、根据试验方案的允许故障数c确定试验方案的总试验时间T。本发明可以对威布尔分布型产品进行可靠性鉴定,更科学、合理地确定试验方案,大大改善产品可靠性评价水平,为产品的设计定型和批量生产把好关,从而提高产品的质量和可靠性。
Description
技术领域
本发明涉及产品可靠性试验技术领域,具体涉及一种适用于威布尔分布型产品的定时截尾可靠性试验方法。
背景技术
可靠性试验是为了了解、评价、分析和提高产品的可靠性而进行各种试验的总称,旨在暴露产品的缺陷,为提高产品的可靠性提供必要信息并最终验证产品的可靠性。换句话说,任何与产品故障或故障效应有关的试验都可以认为是可靠性试验。产品的可靠性是设计、制造出来的,是管理出来的,也是试验出来的,可靠性试验在产品的研制过程中,特别是保障产品达到设计的可靠性要求起到至关重要的作用,尤其是可靠性鉴定试验,它是验证产品的可靠性是否达到设计要求并作为产品设计定型的依据之一,是向使用方提供的一种合格证明。一般用于定型鉴定,是生产前的试验,为生产决策提供管理信息。
可靠性鉴定试验方法包括定数截尾试验、定时截尾试验和序贯截尾试验,其中定时截尾试验是可靠性鉴定试验中用得最多的试验方案,其优点是判决故障数及试验时间、费用在试验前已能确定,便于管理。
定时截尾试验抽验规则为从一批产品中随机抽取n个样品,对其试验到事先规定的截尾时间t时,停止试验。若在规定的(0,t]试验时间内,产品发生的总故障次数为r。将r与允许故障数c进行比较,若r≤c,说明产品的可靠性水平是合格的,接收产品;若r>c,则说明产品的可靠性水平是不合格的,拒收产品。即根据试验期内发生的总故障数r与允许故障数c的比较结果,来判定产品的可靠性水平是否满足要求。
定时截尾试验本质是一种抽样检验,即通过检验一部分样品来判定整批产品。理论上给定平均寿命θ0,当产品的平均寿命θ≥θ0时,产品合格,接收概率L(θ)=1;当产品的平均寿命θ<θ0时,产品不合格,接收概率L(θ)=0。
目前,国内外的大多数标准(如GJB 899A-2009、MIL-HDBK-781A)都以产品故障前工作时间服从指数分布这一假设制定。然而对于机械产品,其寿命或故障数据大多服从威布尔分布。此时,用指数分布的试验方案对故障数据服从威布尔分布的产品做可靠性鉴定试验,结果过于保守。
发明内容
针对现有技术中的上述不足,本发明提供的一种适用于威布尔分布型产品的定时截尾可靠性试验方法可以科学、合理地确定试验方案,大大提高工作效率和改善评价水平。
为了达到上述发明目的,本发明采用的技术方案为:
提供一种适用于威布尔分布型产品的定时截尾可靠性试验方法,其包括以下步骤:
S1、根据历史数据获取威布尔分布的形状参数m;确定定时截尾试验方案的基本参数;
S2、根据威布尔分布的形状参数m和定时截尾试验方案的基本参数确定试验方案的允许故障数c;
S3、根据试验方案的允许故障数c确定试验方案的总试验时间T。
进一步地,步骤S1中根据历史数据获取威布尔分布的形状参数m的方法为:
根据两参数威布尔分布概率密度函数:
得到故障间隔时间的似然函数:
和其余s个无失效数据的似然函数:
根据极大似然法原理,进而得到威布尔分布形状参数m的极大似然函数为:
根据威布尔分布形状参数m的极大似然函数和极大似然方程:
进一步得到两参数威布尔分布参数极大似然计算方程:
根据上述两参数威布尔分布参数极大似然计算方程得到威布尔分布形状参数m;
历史数据包括相似产品或同一产品历史使用数据中的样本数量s,尺度参数η,试验期内发生故障的次数k,对应的故障间隔时间t1≤t2≤…≤tk,各样本对应的截尾时间T1≤T2≤…≤Ts,设定截尾时间tτ,其中各样本对应的截尾时间和故障间隔时间均小于等于设定截尾时间tτ。
进一步地,步骤S1中定时截尾试验方案的基本参数包括:
待试验产品的样本数量n,生产方风险α,使用方风险β,平均故障间隔时间的检验下限θ1,平均故障间隔时间的检验上限θ0和鉴别比d,其中d=θ0/θ1。
进一步地,生产方风险α和使用方风险β的取值均包括10%、20%和30%;鉴别比d的取值包括1.5、2.0和3.0。
进一步地,步骤S2中根据威布尔分布的形状参数m和定时截尾试验方案的基本参数确定允许故障数c的方法为:
根据两参数威布尔分布模型及得到:
F(t)=1-exp(-λtm)
令上式中tm=x0,得到样本n在区间(0,x0]内发生i次故障的概率为:
进而得到在区间(0,x0]内发生总故障数r小于等于允许故障数c的接收概率L(θ):
根据威布尔分布的点估计公式及得到:
根据上式、接收概率L(θ)及产品抽样特征曲线L(θ0)=1-α、L(θ1)=β,得到:
根据泊松过程的统计推断,得:
进而得到试验方案的允许故障数c:
其中:n为样本数量;α为生产方风险;β为使用方风险;d为鉴别比;m为威布尔分布的形状参数;Γ(·)为伽马函数;和分别为自由度为2c+2的χ2分布的上1-α和β分位点。
进一步地,步骤S3中根据试验方案的允许故障数c确定试验方案的总试验时间T的方法为:
根据n个样本在试验截尾时间t内发生故障数r小于等于允许故障数c的接收概率L(θ):
得:
根据故障间隔时间θ服从威布尔分布模型,得到其累计失效次数的数学期望:
根据威布尔分布模型的期望公式和方差公式得到:
进而得到:
上式结合产品抽样曲线L(θ0)=1-α和L(θ1)=β的均值,得到单个产品的试验截尾时间t:
进而得到总试验时间:
T=nt
其中:n为样本数量;α为生产方风险;β为使用方风险;d为鉴别比;θ1为平均故障间隔时间的检验下限;m为威布尔分布的形状参数;c为允许故障数;Γ(·)为伽马函数;分别为自由度为2c+2的χ2分布的上1-α和β分位点。
本发明的有益效果为:本发明首先利用相似产品或历史使用数据获取威布尔分布的形状参数和确定定时截尾试验的基本参数,然后根据这些参数得到试验方案的允许故障数,最后确定试验方案的试验时间。利用本方法对威布尔分布型产品进行可靠性鉴定,可以更科学、合理地确定试验方案,大大改善产品可靠性评价水平,为产品的设计定型和批量生产把好关,从而提高产品的质量和可靠性。
附图说明
图1为本发明的流程图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
如图1所示,该适用于威布尔分布型产品的定时截尾可靠性试验方法包括以下步骤:
S1、根据历史数据获取威布尔分布的形状参数m;确定定时截尾试验方案的基本参数;
S2、根据威布尔分布的形状参数m和定时截尾试验方案的基本参数确定试验方案的允许故障数c;
S3、根据试验方案的允许故障数c确定试验方案的总试验时间T。
步骤S1中根据历史数据获取威布尔分布的形状参数m的方法为:
根据两参数威布尔分布概率密度函数(两参数威布尔属于威布尔中的一种类型):
得到故障间隔时间的似然函数:
和其余s个无失效数据的似然函数:
根据极大似然法原理,进而得到威布尔分布形状参数m的极大似然函数为:
根据威布尔分布形状参数m的极大似然函数和极大似然方程:
进一步得到两参数威布尔分布参数极大似然计算方程:
根据上述两参数威布尔分布参数极大似然计算方程得到威布尔分布形状参数m。
历史数据包括相似产品或同一产品历史使用数据中的样本数量s,尺度参数η,试验期内发生故障的次数k,对应的故障间隔时间t1≤t2≤…≤tk,各样本对应的截尾时间T1≤T2≤…≤Ts,设定截尾时间tτ,其中各样本对应的截尾时间和故障间隔时间均小于等于设定截尾时间tτ。
本发明示例中,直接假设威布尔分布的形状参数m已知,取m=1.2。
步骤S1中定时截尾试验方案的基本参数包括:
待试验产品的样本数量n,生产方风险α,使用方风险β,平均故障间隔时间的检验下限θ1,平均故障间隔时间的检验上限θ0和鉴别比d,其中d=θ0/θ1。
在实际过程中,由于是通过检查一部分样品来判定整批产品,所以可能会犯两类错误。第一类错误是将合格产品误判为不合格产品;第二类错误是将不合格产品误判为合格产品。犯第一类错误时会使生产方遭受损失,所以称犯这类错误的概率为生产方风险,一般用α表示;犯第二类错误时会使使用方遭受损失,所以称犯这类错误的概率为使用方风险,一般用β表示。而生产方和使用方都希望少承担风险,在这种情况下,给出了两个寿命界限θ0和θ1(θ0>θ1),当产品的平均寿命θ≥θ0时,以大概率接收整批产品,L(θ0)=1-α;当产品的平均寿命θ≤θ1时,以小概率接收整批产品,即L(θ1)=β。由以上条件就可确定抽样检验的试验方案。
(1)、对于样本数量n,《GJB 899A-2009》标准中指出:对于可靠性试验,若使用方无其他规定,每批产品至少应有2台接受试验。推荐的样本大小为每批产品的10%,但最多不超过20台,仅在特殊情况(如因安全或完成任务要求)下才采用全数试验(即100%抽样)。即试验方案所需样品数量n应按标准,或由生产方与使用方商定。
(2)、依据《GJB 899A-2009》,生产方风险α和使用方风险β通常取10%、20%和30%,且两类风险尽可能接近。
(3)、通常情况下,MTBF(平均故障间隔时间)的检验下限θ1应取为产品的最低可接受值,但MTBF的检验上限θ0不是取决于规定值,而应参照规定值选取。
(4)、鉴别比d通常选择的是1.5、2.0、3.0。在试验统计方案中,检验上限θ0、检验下限θ1和鉴别比d三个参数中只要确定其中两个,另一个也将随之确定。
本发明示例中,取产品样本数量n=8,决策风险α=β=20%,产品的MTBF检验下限为θ1=500h,鉴别比d=1.5。
步骤S2中根据威布尔分布的形状参数m和定时截尾试验方案的基本参数确定允许故障数c的方法为:
根据两参数威布尔分布模型及得到:
F(t)=1-exp(-λtm)
令上式中tm=x0,得到样本n在区间(0,x0]内发生i次故障的概率为:
进而得到在区间(0,x0]内发生总故障数r小于等于允许故障数c的接收概率L(θ):
上式根据威布尔分布的点估计公式及得到:
根据上式、接收概率L(θ)及产品抽样特征曲线L(θ0)=1-α、L(θ1)=β,得到:
上式根据泊松过程的统计推断,得:
进而得到试验方案的允许故障数c:
其中:n为样本数量;α为生产方风险;β为使用方风险;d为鉴别比;θ1为平均故障间隔时间的检验下限;m为威布尔分布的形状参数;c为允许故障数;Γ(·)是伽马函数;分别为自由度为2c+2的χ2分布的上1-α和β分位点。
本发明示例中,决策风险α=β=20%,鉴别比d=1.5,形状参数m=1.2。由于允许故障数c必须为整数,因此求解允许故障数公式的解法为计算c=0,1,2,……时的不同结果,将计算结果与dm值比较,选取最接近的c。在本发明示例中,c=11时,结果最为接近。故取试验方案的允许故障数c=11。
步骤S3中根据试验方案的允许故障数c确定试验方案的总试验时间T的方法为:
根据n个样本在试验截尾时间t内发生故障数r小于等于允许故障数c的接收概率L(θ):
得:
根据故障间隔时间θ服从威布尔分布模型,得到其累计失效次数的数学期望:
根据威布尔分布模型的期望公式和方差公式得到:
进而得到:
结合产品抽样曲线L(θ0)=1-α和L(θ1)=β的均值,得到单个产品的试验截尾时间:
进而得到总试验时间:
T=nt
其中:n为样本数量;α为生产方风险;β为使用方风险;d为鉴别比;θ1为平均故障间隔时间的检验下限;m为威布尔分布的形状参数;c为允许故障数;Γ(·)是伽马函数;分别为自由度为2c+2的χ2分布的上1-α和β分位点。
本发明示例中,产品样本数量n=8,决策风险α=β=20%,产品的MTBF检验下限为θ1=500h,鉴别比d=1.5,形状参数m=1.2,允许故障数c=11,因此可计算得到试验截尾时间t=979h,总试验时间T=7832h。
即本发明示例所确定的试验方案为:
对8个抽样产品进行试验时间为979h的定时截尾试验,如果故障数r≤c=11,那么认为该系列产品的可靠性水平合格,可接收该产品;如果r>c=11,则认为该系列产品的可靠性水平不合格,拒收该产品。
综上所述,本发明首先利用相似产品或历史使用数据获取威布尔分布的形状参数和确定定时截尾试验的基本参数,然后根据这些参数得到试验方案的允许故障数,最后确定试验方案的试验时间。利用本方法对威布尔分布型产品进行可靠性鉴定,可以更科学、合理地确定试验方案,大大改善产品可靠性评价水平,为产品的设计定型和批量生产把好关,从而提高产品的质量和可靠性。
Claims (6)
1.一种适用于威布尔分布型产品的定时截尾可靠性试验方法,其特征在于,包括以下步骤:
S1、根据历史数据获取威布尔分布的形状参数m;确定定时截尾试验方案的基本参数;
S2、根据威布尔分布的形状参数m和定时截尾试验方案的基本参数确定试验方案的允许故障数c;
S3、根据试验方案的允许故障数c确定试验方案的总试验时间T。
2.根据权利要求1所述的适用于威布尔分布型产品的定时截尾可靠性试验方法,其特征在于,步骤S1中根据历史数据获取威布尔分布的形状参数m的方法为:
根据两参数威布尔分布概率密度函数:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>m</mi>
<mi>&eta;</mi>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>t</mi>
<mi>&eta;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>t</mi>
<mi>&eta;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
<mo>)</mo>
</mrow>
</mrow>
得到故障间隔时间的似然函数:
<mrow>
<mi>l</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>,</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Pi;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mi>m</mi>
<mo>,</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Pi;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mfrac>
<mi>m</mi>
<mi>&eta;</mi>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mi>&eta;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mi>&eta;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
<mo>)</mo>
</mrow>
</mrow>
和其余s个无失效数据的似然函数:
<mrow>
<munderover>
<mo>&Pi;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</munderover>
<mi>R</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
根据极大似然法的原理,进而得到威布尔分布形状参数m的极大似然函数为:
<mrow>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>,</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Pi;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mfrac>
<mi>m</mi>
<mi>&eta;</mi>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mi>&eta;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mi>&eta;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
</mrow>
</msup>
<munderover>
<mo>&Pi;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</munderover>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mi>&eta;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
</mrow>
</msup>
</mrow>
根据威布尔分布形状参数m的极大似然函数和极大似然方程:
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>ln</mi>
<mi> </mi>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>,</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>m</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>ln</mi>
<mi> </mi>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>,</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>&eta;</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
</mrow>
进一步得到两参数威布尔分布参数极大似然计算方程:
<mrow>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msubsup>
<mi>t</mi>
<mi>i</mi>
<mi>m</mi>
</msubsup>
<mi>ln</mi>
<mi> </mi>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</munderover>
<msubsup>
<mi>T</mi>
<mi>j</mi>
<mi>m</mi>
</msubsup>
<mi>ln</mi>
<mi> </mi>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msubsup>
<mi>t</mi>
<mi>i</mi>
<mi>m</mi>
</msubsup>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</munderover>
<msubsup>
<mi>T</mi>
<mi>j</mi>
<mi>m</mi>
</msubsup>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>m</mi>
</mfrac>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>k</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mi>ln</mi>
<mi> </mi>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
</mrow>
根据上述两参数威布尔分布参数极大似然计算方程得到威布尔分布形状参数m;
历史数据包括相似产品或同一产品历史使用数据中的样本数量s,尺度参数η,试验期内发生故障的次数k,对应的故障间隔时间t1≤t2≤…≤tk,各样本对应的截尾时间T1≤T2≤…≤Ts,设定截尾时间tτ,其中各样本对应的截尾时间和故障间隔时间均小于等于设定截尾时间tτ。
3.根据权利要求2所述的适用于威布尔分布型产品的定时截尾可靠性试验方法,其特征在于:步骤S1中定时截尾试验方案的基本参数包括:
待试验产品的样本数量n,生产方风险α,使用方风险β,平均故障间隔时间的检验下限θ1,平均故障间隔时间的检验上限θ0和鉴别比d,其中d=θ0/θ1。
4.根据权利要求3所述的适用于威布尔分布型产品的定时截尾可靠性试验方法,其特征在于,所述生产方风险α和使用方风险β的取值均包括10%、20%和30%;所述鉴别比d的取值包括1.5、2.0和3.0。
5.根据权利要求4所述的适用于威布尔分布型产品的定时截尾可靠性试验方法,其特征在于,步骤S2中根据威布尔分布的形状参数m和定时截尾试验方案的基本参数确定允许故障数c的方法为:
根据两参数威布尔分布模型得到:
F(t)=1-exp(-λtm)
令上式中tm=x0,得到样本n在区间(0,x0]内发生i次故障的概率为:
<mrow>
<mi>Pr</mi>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
<mo>=</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>n&lambda;x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msup>
<mrow>
<mi>i</mi>
<mo>!</mo>
</mrow>
</mfrac>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>n&lambda;x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
</mrow>
进而得到在区间(0,x0]内发生总故障数r小于等于允许故障数c的接收概率L(θ):
<mrow>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>Pr</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>&le;</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>c</mi>
</munderover>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>n&lambda;x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msup>
<mrow>
<mi>i</mi>
<mo>!</mo>
</mrow>
</mfrac>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>n&lambda;x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
根据威布尔分布的点估计公式及得到:
<mrow>
<mi>&lambda;</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>&theta;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
</mrow>
根据上式、接收概率L(θ)及产品抽样特征曲线L(θ0)=1-α、L(θ1)=β,得到:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
<mo>-</mo>
<mi>&alpha;</mi>
<mo>=</mo>
<mstyle>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>c</mi>
</munderover>
</mstyle>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>nx</mi>
<mn>0</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msup>
<mrow>
<mi>i</mi>
<mo>!</mo>
</mrow>
</mfrac>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>nx</mi>
<mn>0</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
</mrow>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&beta;</mi>
<mo>=</mo>
<mstyle>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>c</mi>
</munderover>
</mstyle>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>nx</mi>
<mn>0</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msup>
<mrow>
<mi>i</mi>
<mo>!</mo>
</mrow>
</mfrac>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>nx</mi>
<mn>0</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
</mrow>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
上式根据泊松过程的统计推断,得:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>2</mn>
<msub>
<mi>nx</mi>
<mn>0</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
<mo>=</mo>
<msubsup>
<mi>&chi;</mi>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&alpha;</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>c</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>2</mn>
<msub>
<mi>nx</mi>
<mn>0</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
<mo>=</mo>
<msubsup>
<mi>&chi;</mi>
<mi>&beta;</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>c</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
进而得到试验方案的允许故障数c:
<mrow>
<msup>
<mi>d</mi>
<mi>m</mi>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&chi;</mi>
<mi>&beta;</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>c</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>&chi;</mi>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&alpha;</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>c</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中:n为样本数量;α为生产方风险;β为使用方风险;d为鉴别比;m为威布尔分布的形状参数;Γ(·)为伽马函数;分别为自由度为2c+2的χ2分布的上1-α和β分位点。
6.根据权利要求5所述的适用于威布尔分布型产品的定时截尾可靠性试验方法,其特征在于,步骤S3中根据试验方案的允许故障数c确定试验方案的总试验时间T的方法为:
根据n个样本在试验截尾时间t内发生故障数r小于等于允许故障数c的接收概率L(θ):
<mrow>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>Pr</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>&le;</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>c</mi>
</munderover>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mi>E</mi>
<mo>(</mo>
<mi>N</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msup>
<mrow>
<mi>i</mi>
<mo>!</mo>
</mrow>
</mfrac>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>n</mi>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
</msup>
</mrow>
得:
<mrow>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
</mfrac>
<msub>
<mi>&chi;</mi>
<mrow>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>c</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
根据故障间隔时间θ服从威布尔分布模型,得到其累计失效次数的数学期望:
<mrow>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>&ap;</mo>
<mfrac>
<mi>t</mi>
<mi>&theta;</mi>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>&theta;</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>&theta;</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
根据威布尔分布模型的期望公式和方差公式得到:
<mrow>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>t</mi>
<mi>&theta;</mi>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<mo>/</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>&Gamma;</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mn>1</mn>
</mrow>
进而得到:
<mrow>
<mi>t</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>&chi;</mi>
<mrow>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>c</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<mo>/</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>&Gamma;</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>&theta;</mi>
</mrow>
结合产品抽样曲线L(θ0)=1-α和L(θ1)=β的均值,得到单个产品的试验截尾时间t:
<mrow>
<mi>t</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mo>&times;</mo>
<msub>
<mi>&chi;</mi>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&alpha;</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>c</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&chi;</mi>
<mi>&beta;</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>c</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>4</mn>
<mi>n</mi>
</mrow>
</mfrac>
<mo>&times;</mo>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<mo>/</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>&Gamma;</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
</mrow>
进而得到总试验时间:
T=nt
其中:n为样本数量;α为生产方风险;β为使用方风险;d为鉴别比;θ1为平均故障间隔时间的检验下限;m为威布尔分布的形状参数;c为允许故障数;Γ(·)为伽马函数;分别为自由度为2c+2的χ2分布的上1-α和β分位点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710777965.5A CN107478455B (zh) | 2017-09-01 | 2017-09-01 | 一种适用于威布尔分布型产品的定时截尾可靠性试验方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710777965.5A CN107478455B (zh) | 2017-09-01 | 2017-09-01 | 一种适用于威布尔分布型产品的定时截尾可靠性试验方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107478455A true CN107478455A (zh) | 2017-12-15 |
CN107478455B CN107478455B (zh) | 2020-02-07 |
Family
ID=60603471
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710777965.5A Expired - Fee Related CN107478455B (zh) | 2017-09-01 | 2017-09-01 | 一种适用于威布尔分布型产品的定时截尾可靠性试验方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107478455B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108241909A (zh) * | 2018-01-24 | 2018-07-03 | 长安大学 | 一种机械设备再制造时机的预测方法 |
CN109344472A (zh) * | 2018-09-17 | 2019-02-15 | 中国人民解放军海军工程大学 | 一种估计机电件可靠性参数的方法 |
CN110069750A (zh) * | 2019-04-12 | 2019-07-30 | 电子科技大学 | 基于四参数威布尔分布的非精确分析方法 |
CN110096820A (zh) * | 2019-05-07 | 2019-08-06 | 中国人民解放军国防科技大学 | 单元寿命服从威布尔分布时k/n(G)***的剩余寿命预测方法 |
CN112052577A (zh) * | 2020-08-28 | 2020-12-08 | 中国航发贵阳发动机设计研究所 | 一种航空发动机可靠性指标规定值评价时机的确定方法 |
CN112485018A (zh) * | 2020-11-18 | 2021-03-12 | 西安电子科技大学 | 基于机理数据融合的机械设备能耗异常检测方法 |
CN112733088A (zh) * | 2021-03-30 | 2021-04-30 | 中国人民解放军国防科技大学 | 基于专家信息的威布尔型产品可靠性验收方案设计方法 |
CN112883590A (zh) * | 2021-03-22 | 2021-06-01 | 中国人民解放军国防科技大学 | ***可靠性验证试验方案确定方法 |
CN114580985A (zh) * | 2022-05-09 | 2022-06-03 | 中国人民解放军海军工程大学 | 一种电子产品的质量检测方法及*** |
CN114897349A (zh) * | 2022-05-09 | 2022-08-12 | 中国人民解放军海军工程大学 | 一种成败型序贯抽样试验方案确定***及方法 |
CN114970120A (zh) * | 2022-05-09 | 2022-08-30 | 中国人民解放军海军工程大学 | 一种指数寿命型电子零部件寿命试验方案确定***及方法 |
CN115825793A (zh) * | 2022-11-25 | 2023-03-21 | 楚能新能源股份有限公司 | 一种动力电池日历寿命鉴定试验方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002090267A (ja) * | 2000-09-21 | 2002-03-27 | Mitsui Eng & Shipbuild Co Ltd | 異常診断方法 |
US6996493B1 (en) * | 2004-02-11 | 2006-02-07 | The United States Of America As Represented By The Secretary Of The Navy | Programmed computation of predicted pressure loads on marine vessels |
CN102981452A (zh) * | 2012-12-28 | 2013-03-20 | 吉林大学 | 数控机床三类功能部件的可靠性建模与可靠性评估方法 |
CN104239694A (zh) * | 2014-08-28 | 2014-12-24 | 北京交通大学 | 一种城轨列车转向架的故障预测与视情维修方法 |
CN104462757A (zh) * | 2014-11-03 | 2015-03-25 | 南京航空航天大学 | 基于监测数据的Weibull分布可靠性序贯验证试验方法 |
CN105894133A (zh) * | 2016-05-12 | 2016-08-24 | 华中科技大学 | 一种风电机组部件维修及备品备件需求预测方法 |
CN106443221A (zh) * | 2015-08-04 | 2017-02-22 | 北京电子工程总体研究所 | 一种用于电子产品的可靠性鉴定方法 |
-
2017
- 2017-09-01 CN CN201710777965.5A patent/CN107478455B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002090267A (ja) * | 2000-09-21 | 2002-03-27 | Mitsui Eng & Shipbuild Co Ltd | 異常診断方法 |
US6996493B1 (en) * | 2004-02-11 | 2006-02-07 | The United States Of America As Represented By The Secretary Of The Navy | Programmed computation of predicted pressure loads on marine vessels |
CN102981452A (zh) * | 2012-12-28 | 2013-03-20 | 吉林大学 | 数控机床三类功能部件的可靠性建模与可靠性评估方法 |
CN104239694A (zh) * | 2014-08-28 | 2014-12-24 | 北京交通大学 | 一种城轨列车转向架的故障预测与视情维修方法 |
CN104462757A (zh) * | 2014-11-03 | 2015-03-25 | 南京航空航天大学 | 基于监测数据的Weibull分布可靠性序贯验证试验方法 |
CN106443221A (zh) * | 2015-08-04 | 2017-02-22 | 北京电子工程总体研究所 | 一种用于电子产品的可靠性鉴定方法 |
CN105894133A (zh) * | 2016-05-12 | 2016-08-24 | 华中科技大学 | 一种风电机组部件维修及备品备件需求预测方法 |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108241909B (zh) * | 2018-01-24 | 2022-04-26 | 长安大学 | 一种机械设备再制造时机的预测方法 |
CN108241909A (zh) * | 2018-01-24 | 2018-07-03 | 长安大学 | 一种机械设备再制造时机的预测方法 |
CN109344472A (zh) * | 2018-09-17 | 2019-02-15 | 中国人民解放军海军工程大学 | 一种估计机电件可靠性参数的方法 |
CN109344472B (zh) * | 2018-09-17 | 2023-02-07 | 中国人民解放军海军工程大学 | 一种估计机电件可靠性参数的方法 |
CN110069750A (zh) * | 2019-04-12 | 2019-07-30 | 电子科技大学 | 基于四参数威布尔分布的非精确分析方法 |
CN110096820A (zh) * | 2019-05-07 | 2019-08-06 | 中国人民解放军国防科技大学 | 单元寿命服从威布尔分布时k/n(G)***的剩余寿命预测方法 |
CN110096820B (zh) * | 2019-05-07 | 2020-08-07 | 中国人民解放军国防科技大学 | 单元寿命服从威布尔分布时k/n(G)***的剩余寿命预测方法 |
CN112052577B (zh) * | 2020-08-28 | 2022-09-06 | 中国航发贵阳发动机设计研究所 | 一种航空发动机可靠性指标规定值评价时机的确定方法 |
CN112052577A (zh) * | 2020-08-28 | 2020-12-08 | 中国航发贵阳发动机设计研究所 | 一种航空发动机可靠性指标规定值评价时机的确定方法 |
CN112485018A (zh) * | 2020-11-18 | 2021-03-12 | 西安电子科技大学 | 基于机理数据融合的机械设备能耗异常检测方法 |
CN112485018B (zh) * | 2020-11-18 | 2021-09-10 | 西安电子科技大学 | 基于机理数据融合的机械设备能耗异常检测方法 |
CN112883590A (zh) * | 2021-03-22 | 2021-06-01 | 中国人民解放军国防科技大学 | ***可靠性验证试验方案确定方法 |
CN112883590B (zh) * | 2021-03-22 | 2022-10-18 | 中国人民解放军国防科技大学 | ***可靠性验证试验方案确定方法 |
CN112733088A (zh) * | 2021-03-30 | 2021-04-30 | 中国人民解放军国防科技大学 | 基于专家信息的威布尔型产品可靠性验收方案设计方法 |
CN112733088B (zh) * | 2021-03-30 | 2021-06-18 | 中国人民解放军国防科技大学 | 基于专家信息的威布尔型产品可靠性验收方案设计方法 |
CN114970120A (zh) * | 2022-05-09 | 2022-08-30 | 中国人民解放军海军工程大学 | 一种指数寿命型电子零部件寿命试验方案确定***及方法 |
CN114897349A (zh) * | 2022-05-09 | 2022-08-12 | 中国人民解放军海军工程大学 | 一种成败型序贯抽样试验方案确定***及方法 |
CN114580985A (zh) * | 2022-05-09 | 2022-06-03 | 中国人民解放军海军工程大学 | 一种电子产品的质量检测方法及*** |
CN114897349B (zh) * | 2022-05-09 | 2023-09-05 | 中国人民解放军海军工程大学 | 一种成败型序贯抽样试验方案确定***及方法 |
CN114970120B (zh) * | 2022-05-09 | 2023-10-24 | 中国人民解放军海军工程大学 | 一种指数寿命型电子零部件寿命试验方案确定***及方法 |
CN115825793A (zh) * | 2022-11-25 | 2023-03-21 | 楚能新能源股份有限公司 | 一种动力电池日历寿命鉴定试验方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107478455B (zh) | 2020-02-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107478455B (zh) | 一种适用于威布尔分布型产品的定时截尾可靠性试验方法 | |
CN114935703B (zh) | 一种变频组件自动测试方法、装置及*** | |
CN107391373B (zh) | 基于AutoIT的性能自动化测试方法 | |
CN105702595B (zh) | 晶圆的良率判断方法以及晶圆合格测试的多变量检测方法 | |
CN112733088B (zh) | 基于专家信息的威布尔型产品可靠性验收方案设计方法 | |
CN114720899A (zh) | 退役电池梯次利用分选方法及***、电子设备、存储介质 | |
CN114839464A (zh) | 一种指数寿命型产品的快速检验方法、***、设备及终端 | |
CN105426647B (zh) | 基于可靠度先验信息融合的冷备***可靠度估计方法 | |
CN106886620B (zh) | 航天器测试资源优化配置方法 | |
CN116954624B (zh) | 基于软件开发包的编译方法、软件开发***及服务器 | |
CN117706363A (zh) | 一种光伏储能电池健康状态监控***及方法 | |
CN111124816B (zh) | 一种基于hdp算法的服务器日志分析方法及*** | |
CN114841573B (zh) | 一种指数寿命型产品的分组检验方案确定***及方法 | |
CN115792583A (zh) | 一种车规级芯片的测试方法、装置、设备及介质 | |
CN114971223A (zh) | 一种指数寿命型序贯试验方法、***、设备、介质及终端 | |
CN115407915A (zh) | 眼图工具实现方法、***、电子设备及存储介质 | |
CN111859296B (zh) | 一种基于装备使用期间的测试性指标评估方法及*** | |
CN113962558A (zh) | 一种基于生产数据管理的工业互联网平台评价方法及*** | |
CN114002574A (zh) | 半导体结构的测试方法 | |
CN107491576B (zh) | 基于性能退化数据的导弹部件可靠性分析方法 | |
CN106776236A (zh) | 监测程序的执行的方法和装置 | |
CN108269004B (zh) | 产品寿命分析方法及终端设备 | |
CN112445632A (zh) | 基于故障数据建模的hpc可靠性评估方法 | |
CN112633093B (zh) | 一种装备状态检测方法及*** | |
CN114580985B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200207 Termination date: 20200901 |
|
CF01 | Termination of patent right due to non-payment of annual fee |