CN114841573B - 一种指数寿命型产品的分组检验方案确定系统及方法 - Google Patents
一种指数寿命型产品的分组检验方案确定系统及方法 Download PDFInfo
- Publication number
- CN114841573B CN114841573B CN202210499694.2A CN202210499694A CN114841573B CN 114841573 B CN114841573 B CN 114841573B CN 202210499694 A CN202210499694 A CN 202210499694A CN 114841573 B CN114841573 B CN 114841573B
- Authority
- CN
- China
- Prior art keywords
- test
- samples
- group
- product
- risk
- 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
- 238000007689 inspection Methods 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000012360 testing method Methods 0.000 claims abstract description 260
- 238000004364 calculation method Methods 0.000 claims abstract description 7
- 238000009826 distribution Methods 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 3
- 238000003860 storage Methods 0.000 claims description 2
- 238000012795 verification Methods 0.000 claims 1
- 238000012544 monitoring process Methods 0.000 abstract description 7
- 238000010998 test method Methods 0.000 abstract description 7
- 238000005070 sampling Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 5
- 238000012806 monitoring device Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 210000001503 joint Anatomy 0.000 description 2
- 238000012827 research and development Methods 0.000 description 2
- 241000251468 Actinopterygii Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06395—Quality analysis or management
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Operations Research (AREA)
- Marketing (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Factory Administration (AREA)
- Testing Electric Properties And Detecting Electric Faults (AREA)
Abstract
本发明属于产品质量检验技术领域,公开了一种指数寿命型产品的分组检验方案确定系统及方法,包括:参数初始化模块获取生产方风险阈值、良品平均寿命、使用方风险阈值、差品平均寿命以及每组的试验计划耗时;样本数量计算模块基于初始化后的参数计算样本的总数量和接收数;分组试验模块分别进行第一组样本与第二组样本的检验试验;结果输出模块基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果。本发明不需配套产品状态在线监测设备,降低了相关的试验硬件条件,消除了试验完成时间的不确定性,能在计划的时间内准时完成试验,和相关标准的试验方法一样能把生产方风险和使用方风险都控制在阈值范围以内。
Description
技术领域
本发明属于产品质量检验技术领域,尤其涉及一种指数寿命型产品的分组检验方案确定系统及方法。
背景技术
目前,一般来说,正常使用的电子零部件寿命都服从指数分布,如:印制电路板插件、电子部件、电阻、电容、集成电路等。在理论上,当产品由很多部分组成,无论这些组成部分的寿命是什么分布,只要产品的任一部分出了故障,给予修复再投入使用,则较长时间之后,产品的寿命基本上即指数分布。
目前,在一系列的国际标准、国家标准和军队标准中,介绍了一些常用的指数寿命型产品的检验方法。例如,在GJB899A-2009中,给出了常用的指数寿命型序贯试验统计方案、指数寿命型定时截尾试验统计方案和指数寿命型全数试验统计方案等。这些方案中提及的“试验时间”实际上是所有样本的累计寿命时间,因此,使用样本的寿命数据是这些试验共同的特点。这意味着必须在试验中配置产品状态在线监测设备,用于捕捉产品样本发生故障的“瞬间”,从而得到参试产品样本的寿命数据。指数寿命型序贯试验可以“串行”方式让产品样本一个接一个地进行试验,因此只要有1套产品状态在线监测设备就能开展试验,但完成试验的时间就很漫长,往往是产品平均寿命的数倍甚至数十倍,试验的时间成本也就很高。指数寿命型定时试验能以“并行”方式让所有产品样本同时进行试验,虽然能快速完成试验,但这也意味着需要配置大量的产品状态在线监测设备(用于记录每个产品样本的寿命),大量的监测设备往往意味着较为严苛的试验硬件条件。
在实际工作中,除了试验结果的准确性,试验的时间成本、试验硬件条件等也是选择何种试验时需要考虑的重要因素。针对寿命服从指数分布的产品,目前这些标准的试验方式都需要配套一定数量的产品状态在线监测设备,完成试验的时间也是不确定的,且试验完成时间和以产品状态在线监测设备为代表的试验硬件条件往往是“鱼与熊掌,不可兼得”,这也限制了相关标准提供的指数寿命型试验方法在实际工作中的应用。
通过上述分析,现有技术存在的问题及缺陷为:现有试验方案的试验计划时间长且不确定,成本高,试验硬件条件苛刻,且不能应用于缺少产品状态监测设备等试验硬件条件不足的场景下。
发明内容
针对现有技术存在的问题,本发明提供了一种指数寿命型产品的分组检验方案确定系统及方法。
本发明是这样实现的,一种指数寿命型产品的分组检验方案确定系统包括:
参数初始化模块,用于获取生产方风险阈值、良品平均寿命、使用方风险阈值、差品平均寿命以及每组的试验计划耗时;同时当样本在每组的试验计划耗时内发生故障则为所述样本试验失败;
样本数量计算模块,用于基于初始化后的参数计算样本的总数量和接收数;
分组试验模块,用于分别进行第一组样本与第二组样本的检验试验;
结果输出模块,用于基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果。
本发明的另一目的在于提供一种应用于所述指数寿命型产品的分组检验方案确定系统的指数寿命型产品的分组检验方案确定方法,所述指数寿命型产品的分组检验方案确定方法包括:
以结合产品平均寿命和试验完成时间的可靠度为质量水平,利用基于贝塔分布的样本总数量和接收数计算方法,得到针对任意质量水平值、鉴定比和风险阈值的试验方案。
进一步,所述指数寿命型产品的分组检验方案确定方法包括以下步骤:
步骤一,获取生产方风险阈值、良品平均寿命、使用方风险阈值、差品平均寿命以及每组的试验计划耗时;同时当样本在每组的试验计划耗时内发生故障则为所述样本试验失败;
步骤二,基于初始化后的参数计算样本的总数量和接收数;分别进行第一组样本与第二组样本的检验试验;
步骤三,基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果。
进一步,所述基于初始化后的参数计算样本的总数量和接收数包括:
(1)进行参数初始化:
n=1
Ac=0
At=min(n-1,1.3n(1-q1));
其中,n表示样本总数量,Ac表示接收数;tw表示每组的试验计划耗时;At向上取整数;t0表示良品平均寿命,q0表示良品在试验期间正常工作的概率;t1表示差品平均寿命,q1表示差品在试验期间正常工作的概率;
(2)计算生产方风险概率、使用方风险概率:
其中,Pa表示样本总数量n、接收数Ac对应的生产方风险概率,Pb表示对应的使用方风险概率,B(·)为贝塔函数;
(3)判断Pa<α、Pb<β是否成立;若Pa<α且Pb<β都成立,则计算得到样本总数量n和接收数Ac计算完毕,进行第一组样本的检验试验,否则执行步骤(4);
(4)利用下式更新Ac:
Ac=Ac+1;
判断Ac与At的大小;若Ac>At,则n=n+1,Ac=0,并利用下式更新At:
At=min(n-1,1.3n(1-q1));
其中,At向上取整数;
(5)执行步骤(2)。
进一步,所述进行第一组样本的检验试验包括:
1)从批产品中随机抽取nz1个样本,所有样本同时开始试验,计划试验时间为tw;其中,nz1向上取整数;
2)在试验终点时刻tw,在第一组中统计因出现故障失败的样本数量r11;
3)判断r1是否满足预设关系式:
若r1>Ac成立,则试验完成时间Tw=tw,试验样本数量N=nz1,Nf=r1,则基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果;
若nz1-r1≥n-Ac成立,则试验完成时间Tw=tw,N=nz1,Nf=r1,则基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果;
若r1>Ac或nz1-r1≥n-Ac均不成立,则进行第二组样本的检验试验。
进一步,所述进行第二组样本的检验试验包括:
将第一组中未出现故障的样本共计nz1-r1个,视为新样本选入第二组;再从批产品中随机挑选nz2-nz1+r1个新样本,共同组成第二组样本共计nz2,同时开始试验,计划试验时间为tw;
在第二组的试验终点时刻tw,统计因出现故障失败的样本数量r2,则试验完成时间Tw=2tw,N=nz2+r1,Nf=r1+r2。
进一步,所述基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果包括:
若Nf≤Ac成立,则判断所处产品合格,否则判断所述产品不合格。
本发明的另一目的在于提供一种计算机设备,所述计算机设备包括存储器和处理器,所述存储器存储有计算机程序,所述计算机程序被所述处理器执行时,使得所述处理器执行所述指数寿命型产品的分组检验方案确定方法如下步骤:
步骤一,获取生产方风险阈值、良品平均寿命、使用方风险阈值、差品平均寿命以及每组的试验计划耗时;同时当样本在每组的试验计划耗时内发生故障则为所述样本试验失败;
步骤二,基于初始化后的参数计算样本的总数量和接收数;分别进行第一组样本与第二组样本的检验试验;
步骤三,基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果。
本发明的另一目的在于提供一种计算机可读存储介质,存储有计算机程序,所述计算机程序被处理器执行时,使得所述处理器执行所述指数寿命型产品的分组检验方案确定方法如下步骤:
步骤一,获取生产方风险阈值、良品平均寿命、使用方风险阈值、差品平均寿命以及每组的试验计划耗时;同时当样本在每组的试验计划耗时内发生故障则为所述样本试验失败;
步骤二,基于初始化后的参数计算样本的总数量和接收数;分别进行第一组样本与第二组样本的检验试验;
步骤三,基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果。
本发明的另一目的在于提供一种信息数据处理终端,其特征在于,所述信息数据处理终端用于实现所述指数寿命型产品的分组检验方案确定系统。
结合上述的技术方案和解决的技术问题,请从以下几方面分析本发明所要保护的技术方案所具备的优点及积极效果为:
第一、针对上述现有技术存在的技术问题以及解决该问题的难度,紧密结合本发明的所要保护的技术方案以及研发过程中结果和数据等,详细、深刻地分析本发明技术方案如何解决的技术问题,解决问题之后带来的一些具备创造性的技术效果。具体描述如下:
本发明把相关标准中以平均寿命作为指数寿命型产品质量水平,改为以结合了产品平均寿命和试验完成时间的可靠度为质量水平,消除了以往完成试验计划时间的不确定性,降低了试验的硬件条件。
本发明通过提供基于贝塔分布的样本总数量n和Ac接收数计算方法,补齐了现有标准中只提供典型良品可靠度值q0、鉴定比D和风险阈值α、β试验方案这一短板,实现了从指数寿命型试验到成败型试验的无缝对接,能设计出针对任意质量水平值q0、鉴定比D和风险阈值α、β情况下的试验方案,解决因缺少产品状态监测设备等试验硬件条件不足而无法采用相关标准的指数寿命型试验方法的难题,试验中在上一组中未发生故障的样本直接进入下一组继续试验,减少实际参试样本数量,可有效节省试验经费。
第二,把技术方案看做一个整体或者从产品的角度,本发明所要保护的技术方案具备的技术效果和优点,具体描述如下:
本发明不需配套产品状态在线监测设备,降低了相关的试验硬件条件,消除了试验完成时间的不确定性,能在计划的时间内准时完成试验,和相关标准的试验方法一样能把生产方风险和使用方风险都控制在阈值范围以内。
第三,作为本发明的权利要求的创造性辅助证据,还体现在本发明的技术方案是否克服了技术偏见:一直以来,指数寿命型产品采用指数寿命型试验进行检验,成败型产品采用成败型试验进行检验,指数寿命型试验和成败型试验被视为两种不同类型的抽样检验方法,在相关标准中都只是各自提供了有限数量的典型试验方案,但指数寿命型试验典型方案中的良品t0、差品t1与成败型试验典型方案中的良品q0、差品q1不能适配,即:当此良品t0是彼良品q0时,则此差品t1一定不是彼差品q1。本发明克服了该技术偏见,实现了从指数寿命型试验到成败型试验的有效转换,所提供的试验方案,在降低试验硬件条件,消除完成试验计划时间的不确定性的同时,达到了试验双方风险都低于阈值的要求。
附图说明
图1是本发明实施例提供的指数寿命型产品的分组检验方案确定系统结构示意图;
图2是本发明实施例提供的指数寿命型产品的分组检验方案确定方法流程图;
图3是本发明实施例提供的三种试验的产品接收概率仿真结果示意图;
图4是本发明实施例提供的三种试验的试验完成时间仿真结果示意图;
图5是本发明实施例提供的三种试验的样本总数量仿真结果示意图;
图中:1、参数初始化模块;2、样本数量计算模块;3、分组试验模块;4、结果输出模块。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
一、解释说明实施例。为了使本领域技术人员充分了解本发明如何具体实现,该部分是对权利要求技术方案进行展开说明的解释说明实施例。
如图1所示,本发明实施例提供的指数寿命型产品的分组检验方案确定系统包括:
参数初始化模块1,用于获取生产方风险阈值、良品平均寿命、使用方风险阈值、差品平均寿命以及每组的试验计划耗时;同时当样本在每组的试验计划耗时内发生故障则为所述样本试验失败;
样本数量计算模块2,用于基于初始化后的参数计算样本的总数量和接收数;
分组试验模块3,用于分别进行第一组样本与第二组样本的检验试验;
结果输出模块4,用于基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果。
本发明实施例提供的指数寿命型产品的分组检验方案确定方法包括:
以结合产品平均寿命和试验完成时间的可靠度为质量水平,利用基于贝塔分布的样本总数量和接收数计算方法,得到针对任意质量水平值、鉴定比和风险阈值的试验方案。
如图2所示,本发明实施例提供的指数寿命型产品的分组检验方案确定方法包括以下步骤:
S101,获取生产方风险阈值、良品平均寿命、使用方风险阈值、差品平均寿命以及每组的试验计划耗时;同时当样本在每组的试验计划耗时内发生故障则为所述样本试验失败;
S102,基于初始化后的参数计算样本的总数量和接收数;分别进行第一组样本与第二组样本的检验试验;基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果。
本发明实施例提供的基于初始化后的参数计算样本的总数量和接收数包括:
(1)进行参数初始化:
n=1
Ac=0
At=min(n-1,1.3n(1-q1));
其中,n表示样本总数量,Ac表示接收数;tw表示每组的试验计划耗时;At向上取整数;t0表示良品平均寿命,q0表示良品在试验期间正常工作的概率;t1表示差品平均寿命,q1表示差品在试验期间正常工作的概率;
(2)计算生产方风险概率、使用方风险概率:
其中,Pa表示样本总数量n、接收数Ac对应的生产方风险概率,Pb表示对应的使用方风险概率,B(·)为贝塔函数;
(3)判断Pa<α、Pb<β是否成立;若Pa<α且Pb<β都成立,则计算得到样本总数量n和接收数Ac计算完毕,进行第一组样本的检验试验,否则执行步骤(4);
(4)利用下式更新Ac:
Ac=Ac+1;
判断Ac与At的大小;若Ac>At,则n=n+1,Ac=0,并利用下式更新At:
At=min(n-1,1.3n(1-q1));
其中,At向上取整数;
(5)执行步骤(2)。
本发明实施例提供的进行第一组样本的检验试验包括:
1)从批产品中随机抽取nz1个样本,所有样本同时开始试验,计划试验时间为tw;其中,nz1向上取整数;
2)在试验终点时刻tw,在第一组中统计因出现故障失败的样本数量r1;
3)判断r1是否满足预设关系式:
若r1>Ac成立,则试验完成时间Tw=tw,N=nz1,Nf=r1,则基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果;
若nz1-r1≥n-Ac成立,则试验完成时间Tw=tw,N=nz1,Nf=r1,则基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果;
若r1>Ac或nz1-r1≥n-Ac均不成立,则进行第二组样本的检验试验。
本发明实施例提供的进行第二组样本的检验试验包括:
将第一组中未出现故障的样本共计nz1-r1个,视为新样本选入第二组;再从批产品中随机挑选nz2-nz1+r1个新样本,共同组成第二组样本共计nz2,同时开始试验,计划试验时间为tw;
在第二组的试验终点时刻tw,统计因出现故障失败的样本数量r2,则试验完成时间Tw=2tw,N=nz2+r1,Nf=r1+r2。
本发明实施例提供的基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果包括:
若Nf≤Ac成立,则判断所处产品合格,否则判断所述产品不合格。
下面结合具体实施例对本发明的技术方案作进一步说明。
实施例1:
本发明约定指数分布记为exp(μ),其概率密度函数伽马分布记为Ga(A,B),其概率密度函数/>其中Γ(A)为伽马函数,/>
本发明按惯例以平均寿命t来描述可靠性质量水平。对于试验双方(产品生产方和产品使用方),由于抽样试验始终存在着样本质量偏离产品总体质量的可能性,因而必然存在两种风险:弃真风险和取伪风险。
弃真风险由生产方承担,是指抽样检验结论为“该批产品质量不合格,被拒收”,但实际上该批产品的良品率达到或甚至超过生产方的要求,该批产品质量总体合格。良品的质量水平记为t0,弃真风险通常以质量为t0的产品被拒收的概率来描述,记为Pa,生产方可接受的Pa最大值记为α。弃真风险又称为生产方风险。
取伪风险由使用方承担,是指抽样检验结论为“该批产品质量合格,可接收”,但在实际使用中发现该批产品的质量低于使用方可接受的最低质量水平t1,t1对应的产品本发明称为差品(实际上是合格品的下限)。取伪风险通常以质量为t1的产品被接收的概率来描述,记为Pb,使用方可接受的Pb最大值记为β。取伪风险又称为使用方风险。
在开展指数寿命型试验前,由生产方和使用方约定各自的风险阈值α、β,各自关注的产品质量水平t0、t1(鉴别比),并以此确定其他试验方案的参数。
本发明的核心创新点是参考成败型试验模式,首先把相关标准中以平均寿命作为指数寿命型产品质量水平,改为以结合了产品平均寿命和试验完成时间的可靠度为质量水平,消除了以往试验时间的不确定性。由于此时良品/差品可靠度都不是相关标准中的典型值,本发明接着通过提供基于贝塔分布的样本总数量n和接收数Ac计算方法,补齐了现有标准中只提供典型良品可靠度值q0、鉴定比D和风险阈值α、β试验方案这一短板,实现了从指数寿命型试验到成败型试验的无缝对接,能设计出针对任意质量水平值q0、鉴定比D和风险阈值α、β情况下的试验方案,解决因缺少产品状态监测设备等试验硬件条件不足而无法采用相关标准的指数寿命型试验方法的难题,试验中在上一组中未发生故障的样本直接进入下一组继续试验,减少实际参试样本数量,可有效节省试验经费。
针对试验硬件条件有限(缺少产品状态在线监测设备),且想在计划时间期限内快速完成产品抽样检验试验的需求,本发明提供了一种指数寿命型产品的检验方法,能在把生产方风险和使用方风险稳定控制在预期的阈值范围以内的前提下,不需获取参试产品的寿命数据,可在计划时间内完成检验试验。为便于论述,本发明以分两组为例,建议每组的试验计划时间tw等于t1。本发明的具体步骤如下:
(1)初始化,输入以下参数。
生产方的相关参数:风险阈值α和良品平均寿命t0;
使用方的相关参数:风险阈值β和差品平均寿命t1。
每组的试验计划耗时tw。约定某样本若在tw内发生故障则视为该样本试验失败。
(2)计算样本总数量n和接收数Ac。
(2.1)令
令n=1,Ac=0;
令At等于n-1和1.3n(1-q1)中的最小数,并对At向上取整数。
(2.2)计算Pa、Pb。
令Pa描述了n、Ac对应的生产方风险概率,式中B(·)是贝塔函数;
令Pb描述了n、Ac对应的使用方风险概率。
(2.3)若Pa<α且Pb<β都成立,则样本总数量n和接收数Ac计算完毕,执行(3),否则执行(2.4)。
(2.4)更新Ac,令Ac=Ac+1后,若Ac>At,则令n=n+1,Ac=0,并更新At(令At等于n-1和1.3n(1-q1)中的最小数,并对At向上取整数)。
(2.5)执行(2.2)。
(3)对第一组样本开展检验试验。第一组的样本数量记为nz1,建议nz1等于n-Ac、1+Ac中的最大数,并对nz1向上取整数。
(3.1)从批产品中随机抽取nz1个样本,所有样本同时开始试验,计划试验时间为tw。
(3.2)在试验终点时刻tw,在第一组中统计因出现故障而失败的样本数量r1。
(3.3)若r1>Ac成立,则试验完成时间Tw=tw,N=nz1,Nf=r1,执行(5)。
若nz1-r1≥n-Ac成立,则试验完成时间Tw=tw,N=nz1,Nf=r1,执行(5)。
若r1>Ac或nz1-r1≥n-Ac都不成立,则执行(4)。
(4)对第二组样本开展检验试验。第二组的样本数量记为nz2,令nz2=n-nz1。
(4.1)第一组中未出现故障的样本,共有nz1-r1个,视为新样本选入第二组,然后再从批产品中随机挑选nz2-nz1+r1个新样本,共同组成第二组样本,同时开始试验,计划试验时间为tw。
(4.2)在第二组的试验终点时刻tw,统计因出现故障而失败的样本数量r2,则试验完成时间Tw=2tw,N=nz2+r1,Nf=r1+r2。
(5)输出试验结果。
若Nf≤Ac成立,则接收该批产品,否则拒收该批产品。
二、应用实施例。为了证明本发明的技术方案的创造性和技术价值,该部分是对权利要求技术方案进行具体产品上或相关技术上的应用实施例。
本发明的应用实施例如下:
现在需要对某批次产品进行抽样检验。已知该产品的寿命服从指数分布,生产方关注的良品平均寿命t0=300h,生产方风险阈值α=0.1;使用方关注的差品平均寿命t1=100h,使用方风险阈值β=0.1。现希望在200h内完成对该批次产品的抽样检验,应用本发明,开展抽样检验并给出是否接收的结论。
解:(1)由已知可得:α=0.1,t0=300h,β=0.1,t1=100h。令各组的试验终止时间tw=100h。
(2)计算样本总数量n和接收数Ac。
经计算,q0=0.7165,q1=0.3679。当n=14、Ac=6时,Pa=0.072、Pb=0.098,此时Pa<α且Pb<β都成立。
(3)令nz1=8,从该批产品中随机选择nz1个样本,对第一组样本开展检验试验。在试验终点时刻tw,在第一组中统计出现故障的样本数量r1。(假定)从试验现场得知r1=3。因此时r1>Ac或nz1-r1≥n-Ac都不成立,执行(4)。
(4)对第二组样本开展检验试验。第二组的样本数量nz2=6。第一组样本试验结束后,有5个样本未出现故障,这些样本加上再从该批产品中随机选择的1个样本,组成第二组样本,同时开展检验试验。在试验终点时刻tw,在第二组中统计因出现故障而失败的样本数量r2,(假定)从试验现场得知r2=1,Nf=r1+r2=4。
(5)因Nf≤Ac成立,因此可接收该批产品。
三、实施例相关效果的证据。本发明实施例在研发或者使用过程中取得了一些积极效果,和现有技术相比的确具备很大的优势,下面内容结合试验过程的数据、图表等进行描述。
可用仿真的方式验证本发明的有效性。以上述算例为例,对寿命服从指数分布、平均寿命在50h~400h范围内的产品,从产品接收概率、试验完成时间、试验样本总数量三个角度,分别模拟实现本发明的试验方案、GJB899A-2009的指数寿命型序贯试验方案和指数寿命型定时试验方案的抽样检验效果。其中指数寿命型序贯试验方案中的截尾时间为1169h、接收数为7。指数寿命型定时试验方案的截尾时间为927h、接收数为5,采用所有样本产品同时试验的模式,该方案的试验样本总数量为6(这也是该试验的最小样本总数量)。产品接收概率、试验完成时间、试验样本总数量仿真结果分别如图3、图4、图5.
图3表明,本发明的产品接收概率满足试验要求:对于平均寿命为小于100的差品,其接收概率未超过阈值α(α=0.1);对于平均寿命大于300的良品,其拒收概率未超过阈值β(β=0.1),达到了对试验方案的要求。在平均寿命50h~400h这个范围内,其产品接收概率不低于国标方案的接收概率。
图4表明,相比国标的试验方案,本发明的试验方案不仅消除了试验完成时间的不确定性,甚至还可能用时更短。
图5表明,由于本发明没有使用样本的寿命信息,因此试验样本总数量也明显大于这两种国标试验方案,这也是为降低试验硬件条件付出的代价。不过,上一组试验中未发生故障的样本参加下一组试验的做法,有可能把样本总数量的理论值14最多减小到实际值9,能在一定程度上抵消本发明所需样本总数量较大的不足。
应当注意,本发明的实施方式可以通过硬件、软件或者软件和硬件的结合来实现。硬件部分可以利用专用逻辑来实现;软件部分可以存储在存储器中,由适当的指令执行系统,例如微处理器或者专用设计硬件来执行。本领域的普通技术人员可以理解上述的设备和方法可以使用计算机可执行指令和/或包含在处理器控制代码中来实现,例如在诸如磁盘、CD或DVD-ROM的载体介质、诸如只读存储器(固件)的可编程的存储器或者诸如光学或电子信号载体的数据载体上提供了这样的代码。本发明的设备及其模块可以由诸如超大规模集成电路或门阵列、诸如逻辑芯片、晶体管等的半导体、或者诸如现场可编程门阵列、可编程逻辑设备等的可编程硬件设备的硬件电路实现,也可以用由各种类型的处理器执行的软件实现,也可以由上述硬件电路和软件的结合例如固件来实现。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,都应涵盖在本发明的保护范围之内。
Claims (9)
1.一种指数寿命型产品的分组检验方案确定系统,其特征在于,所述指数寿命型产品的分组检验方案确定系统包括:
参数初始化模块,用于获取生产方风险阈值、良品平均寿命、使用方风险阈值、差品平均寿命以及每组的试验计划耗时;同时当样本在每组的试验计划耗时内发生故障则为所述样本试验失败;
样本数量计算模块,用于基于初始化后的参数计算样本的总数量和接收数;
分组试验模块,用于分别进行第一组样本与第二组样本的检验试验;
结果输出模块,用于基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果;
所述基于初始化后的参数计算样本的总数量和接收数包括:
(1)进行参数初始化:
n=1,
Ac=0,
At=min(n-1,1.3n(1-q1));
其中,n表示样本总数量,Ac表示接收数;tw表示每组的试验计划耗时;At向上取整数;t0表示良品平均寿命,q0表示良品在试验期间正常工作的概率;t1表示差品平均寿命,q1表示差品在试验期间正常工作的概率;
(2)计算生产方风险概率、使用方风险概率:
其中,Pa表示样本总数量n、接收数Ac对应的生产方风险概率,Pb表示对应的使用方风险概率,B(·)为贝塔函数;
(3)判断Pa<α、Pb<β是否成立;若Pa<α且Pb<β都成立,则计算得到样本总数量n和接收数Ac计算完毕,进行第一组样本的检验试验,否则执行步骤(4);α表示生产方风险的阈值,β表示使用方风险的阈值;
(4)利用下式更新Ac:
Ac=Ac+1;
判断Ac与At的大小;若Ac>At,则n=n+1,Ac=0,并利用下式更新At:
At=min(n-1,1.3n(1-q1));
其中,At向上取整数;
(5)执行步骤(2)。
2.一种应用于如权利要求1所述指数寿命型产品的分组检验方案确定系统的指数寿命型产品的分组检验方案确定方法,其特征在于,所述指数寿命型产品的分组检验方案确定方法包括:
以结合产品平均寿命和试验完成时间的可靠度为质量水平,利用基于贝塔分布的样本总数量和接收数计算方法,得到针对任意质量水平值、鉴定比和风险阈值的试验方案。
3.如权利要求2所述指数寿命型产品的分组检验方案确定方法,其特征在于,所述指数寿命型产品的分组检验方案确定方法包括以下步骤:
步骤一,获取生产方风险阈值、良品平均寿命、使用方风险阈值、差品平均寿命以及每组的试验计划耗时;同时当样本在每组的试验计划耗时内发生故障则为所述样本试验失败;
步骤二,基于初始化后的参数计算样本的总数量和接收数;分别进行第一组样本与第二组样本的检验试验;
步骤三,基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果。
4.如权利要求3所述指数寿命型产品的分组检验方案确定方法,其特征在于,所述进行第一组样本的检验试验包括:
1)从批产品中随机抽取nz1个样本,所有样本同时开始试验,计划试验时间为tw;其中,nz1向上取整数;
2)在试验终点时刻tw,在第一组中统计因出现故障失败的样本数量r1;
3)判断r1是否满足预设关系式:
若r1>Ac成立,则试验完成时间Tw=tw,N=nz1,Nf=r1,则基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果;N表示被抽取并进行试验的样本数量,Nf表示在试验中出现故障失败的样本数量;Ac表示接收数;
若nz1-r1≥n-Ac成立,则试验完成时间Tw=tw,N=nz1,Nf=r1,则基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果;
若r1>Ac或nz1-r1≥n-Ac均不成立,则进行第二组样本的检验试验。
5.如权利要求4所述指数寿命型产品的分组检验方案确定方法,其特征在于,所述进行第二组样本的检验试验包括:
将第一组中未出现故障的样本共计nz1-r1个,视为新样本选入第二组;再从批产品中随机挑选nz2-nz1+r1个新样本,共同组成第二组样本共计nz2,同时开始试验,计划试验时间为tw;
在第二组的试验终点时刻tw,统计因出现故障失败的样本数量r2,则试验完成时间Tw=2tw,N=nz2+r1,Nf=r1+r2。
6.如权利要求3所述指数寿命型产品的分组检验方案确定方法,其特征在于,所述基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果包括:
若Nf≤Ac成立,则判断所处产品合格,否则判断所述产品不合格;Ac表示接收数。
7.一种计算机设备,其特征在于,所述计算机设备包括存储器和处理器,所述存储器存储有计算机程序,所述计算机程序被所述处理器执行时,使得所述处理器执行如权利要求2-6任意一项所述指数寿命型产品的分组检验方案确定方法如下步骤:
步骤一,获取生产方风险阈值、良品平均寿命、使用方风险阈值、差品平均寿命以及每组的试验计划耗时;同时当样本在每组的试验计划耗时内发生故障则为所述样本试验失败;
步骤二,基于初始化后的参数计算样本的总数量和接收数;分别进行第一组样本与第二组样本的检验试验;
步骤三,基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果。
8.一种计算机可读存储介质,存储有计算机程序,所述计算机程序被处理器执行时,使得所述处理器执行如权利要求2-6任意一项所述指数寿命型产品的分组检验方案确定方法如下步骤:
步骤一,获取生产方风险阈值、良品平均寿命、使用方风险阈值、差品平均寿命以及每组的试验计划耗时;同时当样本在每组的试验计划耗时内发生故障则为所述样本试验失败;
步骤二,基于初始化后的参数计算样本的总数量和接收数;分别进行第一组样本与第二组样本的检验试验;
步骤三,基于第一组样本与第二组样本的检验试验结果确定产品的总试验结果。
9.一种信息数据处理终端,其特征在于,所述信息数据处理终端用于实现如权利要求1所述指数寿命型产品的分组检验方案确定系统的功能。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210499694.2A CN114841573B (zh) | 2022-05-09 | 2022-05-09 | 一种指数寿命型产品的分组检验方案确定系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210499694.2A CN114841573B (zh) | 2022-05-09 | 2022-05-09 | 一种指数寿命型产品的分组检验方案确定系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114841573A CN114841573A (zh) | 2022-08-02 |
CN114841573B true CN114841573B (zh) | 2023-09-05 |
Family
ID=82571207
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210499694.2A Active CN114841573B (zh) | 2022-05-09 | 2022-05-09 | 一种指数寿命型产品的分组检验方案确定系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114841573B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116911703A (zh) * | 2023-06-29 | 2023-10-20 | 中国空间技术研究院 | 卫星体系协同观测能力的影响因素分析方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103646138A (zh) * | 2013-12-03 | 2014-03-19 | 北京航空航天大学 | 基于贝叶斯理论的定时截尾加速验收抽样试验优化设计方法 |
CN104614195A (zh) * | 2015-02-05 | 2015-05-13 | 中国人民解放军海军工程大学 | 一种纵向抽样型可靠性鉴定试验方法 |
CN104679983A (zh) * | 2014-12-30 | 2015-06-03 | 中国铁道科学研究院电子计算技术研究所 | 一种动车组零部件寿命分析算法 |
CN109388861A (zh) * | 2018-09-17 | 2019-02-26 | 中国人民解放军海军工程大学 | 一种指数型单元寿命分布参数估计方法 |
-
2022
- 2022-05-09 CN CN202210499694.2A patent/CN114841573B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103646138A (zh) * | 2013-12-03 | 2014-03-19 | 北京航空航天大学 | 基于贝叶斯理论的定时截尾加速验收抽样试验优化设计方法 |
CN104679983A (zh) * | 2014-12-30 | 2015-06-03 | 中国铁道科学研究院电子计算技术研究所 | 一种动车组零部件寿命分析算法 |
CN104614195A (zh) * | 2015-02-05 | 2015-05-13 | 中国人民解放军海军工程大学 | 一种纵向抽样型可靠性鉴定试验方法 |
CN109388861A (zh) * | 2018-09-17 | 2019-02-26 | 中国人民解放军海军工程大学 | 一种指数型单元寿命分布参数估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114841573A (zh) | 2022-08-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114839464B (zh) | 一种指数寿命型产品的快速检验方法、系统、设备及终端 | |
CN107478455B (zh) | 一种适用于威布尔分布型产品的定时截尾可靠性试验方法 | |
CN114841573B (zh) | 一种指数寿命型产品的分组检验方案确定系统及方法 | |
CN110634471B (zh) | 一种语音质检方法、装置、电子设备和存储介质 | |
TWI727323B (zh) | 返修板檢測裝置、方法及電腦可讀取存儲介質 | |
CN110969600A (zh) | 一种产品缺陷检测方法、装置、电子设备及存储介质 | |
CN111863649A (zh) | 芯片的成品测试方法、装置、终端设备和存储介质 | |
CN113011594A (zh) | 量子比特的校准方法及装置、量子计算机 | |
CN117233581B (zh) | 一种芯片测试方法、装置、设备及介质 | |
CN114971223B (zh) | 一种指数寿命型序贯试验方法、系统、设备、介质及终端 | |
CN111400644B (zh) | 一种用于实验室分析样品的计算处理方法 | |
CN114897349B (zh) | 一种成败型序贯抽样试验方案确定系统及方法 | |
CN114970120B (zh) | 一种指数寿命型电子零部件寿命试验方案确定系统及方法 | |
TWI490874B (zh) | 使用扇出/扇入矩陣之錯誤捕捉隨機存取記憶體支援技術 | |
CN115936266B (zh) | 轨道交通设备的可靠度预测方法、系统、设备和介质 | |
US20230086626A1 (en) | System and method for detection of anomalies in test and measurement results of a device under test (dut) | |
CN107562617B (zh) | 一种基于nc系统中的模块验证系统 | |
US20240193460A1 (en) | Data processing method and data processing apparatus | |
CN117435386A (zh) | 异常数据定位方法、装置、芯片测试机及可读存储介质 | |
CN116860652A (zh) | 一种软件质量的评价方法、装置、电子设备及存储介质 | |
KR20030087130A (ko) | 반도체 소자의 불량 분석 장치 및 방법 | |
CN115795772A (zh) | 测试数据的评估方法、系统、晶圆测试系统及存储介质 | |
JPH1152023A (ja) | Lsiの故障箇所推定方法及びlsiの故障箇所推定プログラムを記録した記録媒体 | |
CN114490361A (zh) | 测试脚本质量获取方法、装置、计算机设备和存储介质 | |
CN117472767A (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 |