CN103760814A - 基于特征的多品种小批量生产零件过程能力指数确定方法 - Google Patents

基于特征的多品种小批量生产零件过程能力指数确定方法 Download PDF

Info

Publication number
CN103760814A
CN103760814A CN201410043698.5A CN201410043698A CN103760814A CN 103760814 A CN103760814 A CN 103760814A CN 201410043698 A CN201410043698 A CN 201410043698A CN 103760814 A CN103760814 A CN 103760814A
Authority
CN
China
Prior art keywords
machining feature
process capability
sample
measure
tolerance
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
Application number
CN201410043698.5A
Other languages
English (en)
Other versions
CN103760814B (zh
Inventor
李迎光
王伟
郝小忠
刘长青
王慧洁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201410043698.5A priority Critical patent/CN103760814B/zh
Publication of CN103760814A publication Critical patent/CN103760814A/zh
Application granted granted Critical
Publication of CN103760814B publication Critical patent/CN103760814B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Landscapes

  • General Factory Administration (AREA)

Abstract

一种基于特征的多品种小批量生产零件过程能力指数确定方法,其特征是首先以加工特征为对象,把尺寸不同而加工工艺和公差等级相同且标准差与公差之比为常数的加工特征质量特性值定义为同一加工特征样本,其次将样本个体进行归一化,使样本个体均值相同即样本个体服从同一种分布,最后基于加工特征样本计算过程能力指数。本发明从根本上解决了多品种小批量生产的样本容量问题,实现了面向多品种小批量生产的过程能力指数由区间估计向点估计的转变。通过利用获得的过程能力指数能计算出加工特征质量特性值的不合格品率,根据不合格率对特征设计结构进行可加工性分析,选择合适的设计结构,避免100%检测,提高了检测的效率。

Description

基于特征的多品种小批量生产零件过程能力指数确定方法
技术领域
本发明涉及一种机械加工过程控制方法,尤其是一种能大幅度提高加工效率、减少检测次数而不影响总体产品质量的过程控制方法,具体地说是一种基于特征的多品种小批量生产零件过程能力指数的确认方法。
背景技术
过程能力指数(Process Capability Index,PCI)作为过程能力的定量评价,是统计过程控制(Statistical Process Control,SPC)理论的重要内容,现已被广泛应用于产品结构设计、公差优化、成本控制、工艺决策和质量管理等环节。传统的过程能力指数计算对统计样本容量有严格的要求,样本容量越大,获得的过程能力指数就越能反应实际的过程能力。多品种小批量生产由于样本容量较小而无法准确计算过程能力指数。
检索现有技术与文献发现,面向多品种小批量生产的过程能力指数计算方法目前主要有四类:
过程能力指数的区间估计方法。区间估计方法仍采用传统方法计算过程能力指数,只是所获得的过程能力指数具有一定的置信区间。样本容量越小,置信区间范围越大。若要求置信区间下限满足过程能力要求会造成实际的过程能力过于保守。因此,过程能力指数置信区间的实用性受到了很大限制。
基于成组技术的过程能力指数计算方法。为了增加样本容量,可利用成组技术将工艺过程相同或相似、形状相同而公差不同的零件分类成组,通过数据变换把分布类型相同的数据转换为服从同一分布的数据,进而用传统方法计算过程能力指数。成组的合理性依赖于相似度值的取值,相似度标准的建立是该类方法的难点,取值过低会导致样本纳入错误个体,取值过高,亦会造成样本容量不足。
基于Bootstrap的过程能力指数计算方法。Bootstrap方法实质上是一种模拟抽样统计推断方法。它通过重复采样扩充样本容量,进而得到样本统计量的经验分布。但是,基于Bootstrap的过程能力指数计算在本质上并没有增加样本容量,而是将已有的样本个体作为新的个体进行重复抽样。重复抽样的个体并不能完全反映过程能力的实际状况,只是通过重复抽样缩小了置信区间的范围。
基于Bayes估计的过程能力指数计算方法。该方法是根据当前的样本分布和未知分布参数θ的先验分布,利用条件概率分布求解方法估算未知参数θ的条件分布。利用Bayes估计计算过程能力指数的前提是已有先验分布,因此该方法只能计算先验分布已知的样本过程能力指数。
综上所述,针对多品种小批量生产,目前还没合适的过程能力指数计算方法,也没有合适的样本建立方法。
发明内容
本发明的目的是针对多品种小批量生产由于样本容量不足导致过程能力指数难以准确计算的问题,发明一种基于特征的多品种小批量生产零件过程能力指数的确认方法,同时利用获得的过程能力指数计算加工特征质量特性值的不合格品率,根据不合格率对特征设计结构进行可加工性分析,选择合适的设计结构;自动选择典型加工工艺方案,使得零件上每个质量特性值都达到所需的合格品率;对特征各加工状态进行检测必要性分析,避免100%检测,提高了检测的效率。
本发明的技术方案是:
一种基于特征的多品种小批量生产零件过程能力指数的确定方法,其特征是它包括以下步骤:
步骤1、建立加工特征样本,把尺寸不同而加工工艺和公差等级相同且标准差与公差之比为常数的加工特征的质量特性值定义为同一加工特征样本,从而扩充样本容量;
质量特性值是指:产品规格属性的固有要求。例如零件上筋特征的筋厚、孔特征的直径和深度等。
质量特性值的公差,例如筋厚的公差为
Figure BDA0000463806100000021
称为质量特性值的技术规格。
步骤2、将加工特征样本进行归一化处理,使样本个体均值相同即样本个体服从同一种分布;
步骤3、根据归一化后的加工特征样本计算出过程能力指数;
步骤4、根据计算所得的过程能力指数确定加工特征的不合格品率;
步骤5、根据加工特征的不合格品率对特征的设计结构进行可加工性分析,选择合格品率较高的设计结构;
步骤6、根据加工特征的不合格品率选择典型工艺方案,使得零件上所有加工特征的合格品率均符合要求;
步骤7、根据加工特征的不合格品率对加工特征进行检测必要性分析,避免100%检测,提高检测效率。
所述建立加工特征样本的方法为:
步骤1、选定加工特征的某个质量特性值,所述的质量特性值的公称尺寸不同,但公差等级和加工工艺相同;
步骤2、计算不同公称尺寸质量特性值的标准差、公差以及标准差与公差的比值;
步骤3、判定标准差与公差比值是否为常数,若为常数,则将这些质量特性值归为同一加工特征样本。
所述判定标准差与公差之比为常数的方法为:
步骤1、计算标准差在一定置信水平下的置信区间;
步骤2、计算标准差与公差比值在一定置信水平下的置信区间;
步骤3、落在置信区间内的质量特性值归为同一加工特征样本。
所述加工特征样本归一化处理包括以下步骤:
步骤1、进行公称尺寸转换;
步骤2、进行公差转换;
步骤3、进行方差一致性检验;
步骤4、均值一致性检验。
所述基于加工特征样本计算过程能力指数包括以下步骤:
步骤1、计算无偏移双侧规范情形下的过程能力指数;
步骤2、计算无偏移单侧规范情形下的过程能力指数;
步骤3、计算有偏移情形下的过程能力指数。
根据加工特征的不合格品率对特征设计结构进行可加工性分析的方法为:
步骤1、获取加工特征列表;
步骤2、获取加工特征的质量特性值列表;加工特征的质量特性值列表是指一个加工特征含有的多个质量特性值组成的列表。例如筋特征具有筋厚、筋高、筋长、表面粗糙度等多个质量特性值,每个质量特性值都有各自的技术规格(比如公差、粗糙度等级)。
步骤3、查询质量特性值在采用不同加工工艺时的过程能力指数Cp,判断是否大于1.33;
步骤4、若Cp大于等于1.33,则该质量特性值的不合格品率满足要求;
步骤5、若Cp小于1.33,则判定是否有可替代结构,若有给出可替代结构。
根据加工特征的不合格品率选择零件的典型工艺方案的方法为:
步骤1、获取加工特征列表;
步骤2、获取加工特征的质量特性值列表;
步骤3、查询质量特性值在采用不同加工工艺时的过程能力指数Cp
步骤4、选取质量特性值的过程能力指数Cp均大于等于1.33时的典型工艺方案作为实际加工方案。
根据加工特征的不合格品率对加工特征进行检测必要性分析的方法为:
步骤1、获取加工特征列表;
步骤2、获取加工特征的质量特性值列表;
步骤3、查询质量特性值在采用不同加工工艺时的过程能力指数Cp
步骤4、当过程能力指数Cp≥1.67时,表示过程能力过高,零件中间状态和最终状态均不需检测;当1.67>Cp≥1.33时,表示过程能力充分,零件中间状态不需检测,零件最终状态需检测;当1.33>Cp≥1.0时,表示过程能力一般,零件中间状态和最终状态均需检测;通过中间状态检测获得零件上关键尺寸的中间状态,确定判定是否超差,依据超差量调整刀轨,保证关键尺寸最终落在设计公差范围内;当Cp<1.0时,过表示程能力不足,零件最终状态需要检测,零件中间状态也需要检测;零件中间状态检测的目的是查找过程能力不足的原因,进而改进工艺或修改零件结构使得过程能力符合要求。
本发明的有益效果是:
1、提供了面向多品种小批量生产的加工特征样本建立方法。
2、提供了加工特征样本的归一化方法。
3、提供了基于加工特征样本计算过程能力指数的方法。
4、提供了基于过程能力指数对特征设计结构进行可加工性分析的方法。
5、提供了基于过程能力指数选择典型工艺方案的方法。
6、提供了基于过程能力指数进行检测必要性分析的方法。
附图说明
图1为验证本发明建立的筋特征样本。
图2为质量特性值的总体平均值μ与技术规范中心M不重合时不合格品率计算。
图3为基于过程能力指数的零件可加工性分析流程。
图4为本发明的槽特征典型工艺示意图。
具体实施方式
下面结合附图和实施例对本发明的技术方案进行详细说明。
如图1-3所示。
一种基于特征的多品种小批量生产零件过程能力指数的确定方法,它包括以下步骤:
步骤1、建立加工特征样本,把尺寸不同而加工工艺和公差等级相同且标准差与公差之比为常数的加工特征的质量特性值定义为同一加工特征样本,从而扩充样本容量;
加工特征样本建立步骤为:
1、选定加工特征的某个质量特性值,例如筋特征的筋厚、壁特征的壁厚、孔特征的孔径、孔深等。这些质量特性值的公称尺寸不同,但公差等级和加工工艺相同。
2、计算不同公称尺寸质量特性值的标准差、公差以及标准差与公差的比值。计算方法为:1)、计算标准差在一定置信水平下的置信区间;2)、计算标准差与公差比值在一定置信水平下的置信区间;3)、落在置信区间内的质量特性值归为同一加工特征样本。
3、判定标准差与公差比值是否为常数,若为常数,则将这些质量特性值归为同一加工特征样本。
为了验证本发明根据不同尺寸筋特征与加工工艺的对应关系,选取生产批量较大的民机零件建立了图1中的样本1至样本6。由于传统过程能力指数计算方法要求样本容量至少大于100,考虑卡方分布表可提供的数据,样本1至样本6的样本容量均取116。且样本1到样本6均服从正态分布。
由于样本服从N(u,σ),根据抽样分布定理有:
对于给定的置信水平r=1-α,样本总体方差σ2的置信区间为:
( ( n - 1 ) S 2 &chi; &alpha; / 2 2 ( n - 1 ) , ( n - 1 ) S 2 &chi; 1 - &alpha; / 2 2 ( n - 1 ) ) - - - ( 1 )
取α=0.05,查卡方分布表得:
&chi; &alpha; / 2 2 ( n - 1 ) = &chi; 0.05 / 2 2 ( 116 - 1 ) = &chi; 0.025 2 ( 115 ) = 146.57105 - - - ( 2 )
&chi; 1 - &alpha; / 2 2 ( n - 1 ) = &chi; 1 - 0.05 / 2 2 ( 116 - 1 ) = &chi; 0.975 2 ( 115 ) = 87.21279 - - - ( 3 )
将样本1和式(2)、(3)的数据代入式(1),得到方差σ2的置信区间为(0.0013,0.002185),标准差σ的置信区间为(0.036056,0.046744),σ/T的置信区间为(0.09014,0.11686)。
样本2、3、4、5和6的标准差与相应公差值的比值均落在该区间内。因此,在95%的置信水平下,样本1至样本6的标准差与公差之比为常数:
S 1 T 1 = S 2 T 2 = . . . = S i T i = . . . S n T n - - - ( 4 )
式(4)中,Si为区间内第i种尺寸筋厚的标准差,Ti为i种尺寸筋厚的公差,n为区间内筋厚尺寸的个数。
根据过程能力指数的定义Cpi=Ti/6Si,只要样本符合式(4),过程能力指数即相等。因此,尺寸不同而加工工艺和公差等级相同的加工特征质量特性值若存在标准差与公差之比为常数,就将其归为同一样本。
参照数理统计中总体、个体和样本的定义,加工特征总体、加工特征个体和加工特征样本定义如下:
加工特征总体、个体:加工特征某一质量特性值的全部可能观察值称为加工特征总体。该质量特性值的每一个观察值称为加工特征个体。以筋特征为例,取筋特征的厚度为质量特性值,则筋特征厚度的全部可能观察值为加工特征总体,每一个观察值为一个加工特征个体。
从加工特征总体、个体定义可以看出:若加工特征具有多个质量特性值,则该加工特征有多个加工特征总体。
加工特征样本:设加工特征某一质量特性值X是具有分布函数F的随机变量,若X1,X2,…,Xn是具有同一分布函数F的、相互独立的随机变量,则称X1,X2,…,Xn为分布函数F得到的容量为n的加工特征简单随机样本,简称加工特征样本。加工特征样本的观察值x1,x2,…,xn称为加工特征样本值,又称为该加工特征某一质量特性值X的n个独立的观察值。
步骤2、将加工特征样本进行归一化处理,使样本个体均值相同即样本个体服从同一种分布;
归一化处理的步骤为:
加工特征样本归一化方法的步骤为:
1、公称尺寸转换。转换公式为:
u ij = x ij - ( T U i + T L i ) 2 T U i - T L i = x ij - M i T i - - - ( 5 )
式(5)中,uij是加工特征样本中第i种尺寸的第j个观察值的归一化结果,xij是加工特征样本中第i种尺寸的第j个观察值,
Figure BDA0000463806100000071
是第i种尺寸的规格上限,是第i种尺寸的规格下限,Ti是第i种尺寸的技术规格。
2、公差转换。转换公式为:
T Ui u = T Ui - M i T i = T Ui - T Ui + T Li 2 T Ui - T Li = 1 2 ( T Ui - T Li T Ui - T Li ) = 1 2 - - - ( 6 )
T Li u = T Li - M i T i = T Li - T Ui + T Li 2 T Ui - T Li = - 1 2 ( T Ui - T Li T Ui - T Li ) = - 1 2 - - - ( 7 )
M i u = T Ui u + T Li u 2 = 1 2 + ( - 1 2 ) 2 = 0 - - - ( 8 )
T i u = T Ui u - T Li u = 1 2 - ( - 1 2 ) = 1 - - - ( 9 )
式(6)至(9)中,
Figure BDA0000463806100000078
分别是第i种尺寸样本公差上、下限的归一化值,
Figure BDA0000463806100000079
是第i种尺寸样本公差的归一化值。
3、方差一致性检验
以第一种尺寸样本和第二种尺寸样本为例:
H0为归一化后第一种尺寸和第二种尺寸方差没有显著差异,
H1为归一化后第一种尺寸和第二种尺寸方差有显著差异,
Figure BDA00004638061000000711
根据抽样分布定理有: ( n - 1 ) S 2 &sigma; 2 ~ &chi; 2 ( n - 1 ) .
用ni,μi和Si分别表示第i种尺寸样本的样本容量,样本平均值和样本标准差,则有:
( n 1 - 1 ) S 1 2 &sigma; 1 2 ~ &chi; 2 ( n 1 - 1 ) - - - ( 10 )
( n 2 - 1 ) S 2 2 &sigma; 2 2 ~ &chi; 2 ( n 2 - 1 ) - - - ( 11 )
式(10)和式(11)中,
S 1 2 = 1 n 1 - 1 &Sigma; j = 1 n 1 ( x 1 j - &mu; 1 ) 2 - - - ( 12 )
S 2 2 = 1 n 2 - 1 &Sigma; j = 1 n 2 ( x 2 j - &mu; 2 ) 2 - - - ( 13 )
又因为若U~χ2(n1),V~χ2(n2),且U和V相对独立,则随机变量
Figure BDA0000463806100000083
服从自由度为(n1,n2)的F分布,即F~F(n1,n2)。
因为
Figure BDA0000463806100000084
相互独立,所以有:
S 1 2 / &sigma; 1 2 S 2 2 / &sigma; 2 2 ~ F ( n 1 - 1 , n 2 - 1 ) - - - ( 14 )
p { F 1 - &alpha; 2 ( n 1 - 1 , n 2 - 1 ) &le; S 1 2 / &sigma; 1 2 S 2 2 / &sigma; 2 2 &le; F &alpha; 2 ( n 1 - 1 , n 2 - 1 ) } = 1 - &alpha; - - - ( 15 )
式(15)中,
Figure BDA0000463806100000089
是分布显著性水平为α时F分布的双侧分位数。
若H0成立,则 T 1 2 T 2 2 F 1 - &alpha; 2 ( n 1 - 1 , n 2 - 1 ) &le; S 1 2 S 2 2 &le; T 1 2 T 2 2 F &alpha; 2 ( n 1 - 1 , n 2 - 1 ) .
若H1成立,则 S 1 2 S 2 2 < T 1 2 T 2 2 F 1 - &alpha; 2 ( n 1 - 1 , n 2 - 1 ) S 1 2 S 2 2 > T 1 2 T 2 2 F &alpha; 2 ( n 1 - 1 , n 2 - 1 ) .
其余尺寸样本和第一种尺寸均行方差检验即完成加工特征样本的方差检验。
4、均值一致性检验
已知
Figure BDA00004638061000000813
归一化后有:
u ij = x ij - M i T i - - - ( 16 )
&mu; i u = 1 n i &Sigma; j = 1 n i u ij - - - ( 17 )
u i ~ N ( &mu; i u , &sigma; i 2 T i 2 ) - - - ( 18 )
式(17)和式(18)中,
Figure BDA00004638061000000817
是第i种尺寸样本的样本均值μi的归一化结果。
根据之前的证明结果令
Figure BDA00004638061000000818
则有:
u i ~ N ( &mu; i u , &sigma; 2 ) - - - ( 19 )
u &OverBar; i ~ N ( &mu; i u , &sigma; 2 n i ) - - - ( 20 )
u &OverBar; ~ N ( &mu; u , &sigma; 2 n ) - - - ( 21 )
&mu; u = 1 n &Sigma; i = 1 m &Sigma; j = 1 n i u ij = 1 n &Sigma; i = 1 m n i &mu; i u , n = &Sigma; i = 1 m n i - - - ( 22 )
式(21)和式(22)中,n是加工特征样本容量,m为尺寸类型的个数,μu是加工特征样本的样本均值。
H0为归一化后均值没有显著差异,
Figure BDA0000463806100000095
H1为归一化后均值有显著差异,不全相等。
可采用单因素方差分析法进行假设检验:
令SA表示尺寸不同引起的偏差,有:
令SB表示同尺寸样本内随机波动引起的偏差,有:
Figure BDA0000463806100000098
在给定显著性水平α下查F分布表,获得Fα(m-1,n-m)。
若H0成立,则 S A / ( m - 1 ) S B / n - m < F &alpha; ( m - 1 , n - m ) .
若H1成立,则 S A / ( m - 1 ) S B / n - m &GreaterEqual; F &alpha; ( m - 1 , n - m ) .
加工特征往往具有多种典型工艺方案,以槽特征为例,槽特征具有两种典型工艺方案,如图4所示,当采用典型工艺方案1时,槽特征转角处的表面粗糙度的过程能力指数仅为0.09,当采用典型工艺方案2时,槽特征转角处的表面粗糙度的过程能力指数为1.69。所以应该选用第二种典型工艺方案作为实际加工方案。
步骤3、根据归一化后的加工特征样本计算出过程能力指数;
过程能力指数计算常用的有无偏移双侧规范情形、无偏移单侧规范情形和有偏移情形的过程能力指数,其余类型的过程能力指数均可按照类似的方法计算。详述如下:
1:无偏移双侧规范情形过程能力指数计算
无偏移双侧规范情形的过程能力指数Cp计算如下:
C p = T 6 &sigma; = T U - T L 6 &sigma; - - - ( 23 )
式(23)中,T为样本的技术规格幅度,TU、TL分别为样本的上、下规格界限,σ为质量特性值分布的总体标准差。当样本分布服从正态分布时,可用样本的标准差S替换σ,得到Cp的估计值
Figure BDA0000463806100000102
S = 1 n - 1 &Sigma; i = 1 n ( X i - X &OverBar; ) 2 - - - ( 24 )
C ^ p = T 6 S = T U - T l 6 S - - - ( 25 )
式(24)中,Xi为第i个个体质量特性值,
Figure BDA0000463806100000105
是样本均值,n为样本容量。
加工特征样本的公差上下界限
Figure BDA0000463806100000106
公差中心Mu和公差带宽Tu
T U u = T Ui u = 1 2 - - - ( 26 )
T L u = T Li u = - 1 2 - - - ( 27 )
M u = T U u + T L u 2 = 1 2 + ( - 1 2 ) 2 = 0 - - - ( 28 )
T u = T U u - T L u = T i u = 1 - - - ( 29 )
加工特征样本的标准差Su
S u = 1 &Sigma; i m n i - 1 ( &Sigma; i = 1 m &Sigma; j = 1 n i ( u ij - &mu; u ) 2 ) - - - ( 30 )
加工特征样本的过程能力指数
Figure BDA00004638061000001012
C ^ p u = T u 6 S u = 1 6 S u - - - ( 31 )
2:无偏移单侧规范情形过程能力指数
当公差仅有上限要求而无下限要求时,则有上单侧过程能力指数CpU
C pU = T U - &mu; 3 &sigma; - - - ( 32 )
当样本分布服从正态分布时,可用样本的标准差S替换σ,用样本的平均值
Figure BDA00004638061000001015
替换μ得到CpU的估计值
Figure BDA00004638061000001016
C ^ pU = T U - X &OverBar; 3 S ( X &OverBar; < T U ) - - - ( 33 )
当公差仅有下限要求而无下限要求时,则有下单侧过程能力指数CpL
C pL = &mu; - T L 3 &sigma; - - - ( 34 )
当样本分布服从正态分布时,可用样本的标准差S替换σ,用样本的平均值
Figure BDA0000463806100000113
替换μ得到CpL的估计值
Figure BDA0000463806100000114
C ^ pL = X &OverBar; - T L 3 S ( X &OverBar; > T L ) - - - ( 35 )
将加工特征样本公差上限
Figure BDA0000463806100000116
公差下限均值μu和标准差Su代入式(33)和(35)得到加工特征样本无偏移单侧规范情形的过程能力指数:
C ^ pU u = T U u - &mu; u 3 S u = 1 2 - &mu; u 3 S u - - - ( 36 )
C ^ pL u = &mu; u - T L u 3 S u = &mu; u + 1 2 3 S u - - - ( 37 )
3:有偏移情形过程能力指数计算
当质量特性值的总体平均值μ与技术规范中心M不重合时,产品不合格率增大,Cp值降低,如图2所示,阴影部分表示产品不合格率。
定义质量特性值分布中心μ与技术规范中心M的偏差为ε,偏移度为K,则有:
ε=|M-μ|   (38)
K = &epsiv; T / 2 = 2 &epsiv; T - - - ( 39 )
过程能力指数Cp被修正为Cpk
C pk = ( 1 - K ) C p = ( 1 - K ) T 6 &sigma; - - - ( 40 )
当样本分布服从正态分布时,可用样本的标准差S替换σ,利用样本的平均值
Figure BDA00004638061000001112
替换μ得到Cpk的估计值
Figure BDA00004638061000001113
&epsiv; ^ = | M - X &OverBar; | - - - ( 41 )
K ^ = &epsiv; ^ T / 2 = 2 &epsiv; ^ T - - - ( 42 )
C ^ pk = ( 1 - K ^ ) T 6 S - - - ( 43 )
将加工特征样本公差中心Mu、样本均值μu、公差范围Tu和标准差Su代入式(41)至(43),得到多维特征的有偏移情形过程能力指数:
&epsiv; ^ u = | M u - &mu; u | - - - ( 44 )
K ^ u = &epsiv; ^ u T u / 2 = 2 &epsiv; ^ u T u = 2 &epsiv; ^ u - - - ( 45 )
C ^ pk u = ( 1 - K ^ u ) T u 6 S u - - - ( 46 )
步骤4、根据计算所得的过程能力指数确定加工特征的不合格品率;
步骤5、根据加工特征的不合格品率对特征的设计结构进行可加工性分析,选择合格品率较高的设计结构;具体步骤为:
1、获取加工特征列表;
2、获取加工特征的质量特性值列表(如筋厚、孔径等);
3、查询质量特性值在采用不同加工工艺时的过程能力指数Cp,判断是否大于1.33;
4、若Cp大于等于1.33,则该质量特性值的不合格品率满足要求;
5、若Cp小于1.33,则判定是否有可替代结构,若有给出可替代结构。
步骤6、根据加工特征的不合格品率选择典型工艺方案,使得零件上所有加工特征的合格品率均符合要求;选择方法为:
1、获取加工特征列表;
2、获取加工特征的质量特性值列表;
3、查询质量特性值在采用不同加工工艺时的过程能力指数Cp
4、选取质量特性值的过程能力指数Cp均大于等于1.33时的工艺方案作为典型工艺方案。
步骤7、根据加工特征的不合格品率对加工特征进行检测必要性分析,避免100%检测,提高检测效率。具体步骤为:
1、获取加工特征列表;
2、获取加工特征的质量特性值列表;
3、查询质量特性值在采用不同加工工艺时的过程能力指数Cp
4、当过程能力指数Cp≥1.67时,表示过程能力过高,零件中间状态和最终状态均不需检测;当1.67>Cp≥1.33时,表示过程能力充分,零件中间状态不需检测,零件最终状态需检测;当1.33>Cp≥1.0时,表示过程能力一般,零件中间状态和最终状态均需检测;通过中间状态检测获得零件上关键尺寸的中间状态,确定判定是否超差,依据超差量调整刀轨,保证关键尺寸最终落在设计公差范围内;当Cp<1.0时,过表示程能力不足,零件最终状态需要检测,零件中间状态也需要检测;零件中间状态检测的目的是查找过程能力不足的原因,进而改进工艺或修改零件结构使得过程能力符合要求。
本发明未涉及部分与现有技术相同或可采用现有技术加以实现。

Claims (8)

1.一种基于特征的多品种小批量生产零件过程能力指数的确定方法,其特征是它包括以下步骤:
步骤1、建立加工特征样本,把尺寸不同而加工工艺和公差等级相同且标准差与公差之比为常数的加工特征的质量特性值定义为同一加工特征样本,从而扩充样本容量;
步骤2、将加工特征样本进行归一化处理,使样本个体均值相同即样本个体服从同一种分布;
步骤3、根据归一化后的加工特征样本计算出过程能力指数;
步骤4、根据计算所得的过程能力指数确定加工特征的不合格品率;
步骤5、根据加工特征的不合格品率对特征的设计结构进行可加工性分析,选择合格品率较高的设计结构;
步骤6、根据加工特征的不合格品率选择典型工艺方案,使得零件上所有加工特征的合格品率均符合要求;
步骤7、根据加工特征的不合格品率对加工特征进行检测必要性分析,避免100%检测,提高检测效率。
2.如权利要求1所述的方法,其特征在于所述建立加工特征样本的方法为:
步骤1、选定加工特征的某个质量特性值,所述的质量特性值的公称尺寸不同,但公差等级和加工工艺相同;
步骤2、计算不同公称尺寸质量特性值的标准差、公差以及标准差与公差的比值;
步骤3、判定标准差与公差比值是否为常数,若为常数,则将这些质量特性值归为同一加工特征样本。
3.如权利要求1所述的方法,其特征在于所述加工特征样本归一化处理包括以下步骤:
步骤1、进行公称尺寸转换;
步骤2、进行公差转换;
步骤3、进行方差一致性检验;
步骤4、均值一致性检验。
4.如权利要求1所述的方法,其特征在于所述基于加工特征样本计算过程能力指数包括以下步骤:
步骤1、计算无偏移双侧规范情形下的过程能力指数;
步骤2、计算无偏移单侧规范情形下的过程能力指数;
步骤3、计算有偏移情形下的过程能力指数。
5.如权利要求1所述的方法,其特征在于根据加工特征的不合格品率对特征设计结构进行可加工性分析的方法为:
步骤1、获取加工特征列表;
步骤2、获取加工特征的质量特性值列表;
步骤3、查询质量特性值在采用不同加工工艺时的过程能力指数Cp,判断是否大于1.33;
步骤4、若Cp大于等于1.33,则该质量特性值的不合格品率满足要求;
步骤5、若Cp小于1.33,则判定是否有可替代结构,若有给出可替代结构。
6.如权利要求1所述的方法,其特征在于根据加工特征的不合格品率选择零件的典型工艺方案的方法为:
步骤1、获取加工特征列表;
步骤2、获取加工特征的质量特性值列表;
步骤3、查询质量特性值在采用不同加工工艺时的过程能力指数Cp
步骤4、选取质量特性值的过程能力指数Cp均大于等于1.33时的典型工艺方案作为实际加工方案。
7.如权利要求1所述的方法,其特征在于根据加工特征的不合格品率对加工特征进行检测必要性分析的方法为:
步骤1、获取加工特征列表;
步骤2、获取加工特征的质量特性值列表;
步骤3、查询质量特性值在采用不同加工工艺时的过程能力指数Cp
步骤4、当过程能力指数Cp≥1.67时,表示过程能力过高,零件中间状态和最终状态均不需检测;当1.67>Cp≥1.33时,表示过程能力充分,零件中间状态不需检测,零件最终状态需检测;当1.33>Cp≥1.0时,表示过程能力一般,零件中间状态和最终状态均需检测;通过中间状态检测获得零件上关键尺寸的中间状态,确定判定是否超差,依据超差量调整刀轨,保证关键尺寸最终落在设计公差范围内;当Cp<1.0时,过表示程能力不足,零件最终状态需要检测,零件中间状态也需要检测;零件中间状态检测的目的是查找过程能力不足的原因,进而改进工艺或修改零件结构使得过程能力符合要求。
8.如权利要求2所述的方法,其特征在于所述判定标准差与公差之比为常数的方法为:
步骤1、计算标准差在一定置信水平下的置信区间;
步骤2、计算标准差与公差比值在一定置信水平下的置信区间;
步骤3、落在置信区间内的质量特性值归为同一加工特征样本。
CN201410043698.5A 2014-01-29 2014-01-29 基于特征的多品种小批量生产零件过程能力指数确定方法 Active CN103760814B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410043698.5A CN103760814B (zh) 2014-01-29 2014-01-29 基于特征的多品种小批量生产零件过程能力指数确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410043698.5A CN103760814B (zh) 2014-01-29 2014-01-29 基于特征的多品种小批量生产零件过程能力指数确定方法

Publications (2)

Publication Number Publication Date
CN103760814A true CN103760814A (zh) 2014-04-30
CN103760814B CN103760814B (zh) 2016-04-06

Family

ID=50528069

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410043698.5A Active CN103760814B (zh) 2014-01-29 2014-01-29 基于特征的多品种小批量生产零件过程能力指数确定方法

Country Status (1)

Country Link
CN (1) CN103760814B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104133370A (zh) * 2014-07-02 2014-11-05 方志耕 一种基于过程能力的小样本最优质量控制方法
CN104809306A (zh) * 2015-05-13 2015-07-29 梧州学院 基于多质量指标的注塑成型工艺参数优化方法
CN105182903A (zh) * 2015-10-19 2015-12-23 南京航空航天大学 同一加工特征定义方法
CN107110661A (zh) * 2014-12-05 2017-08-29 赛峰飞机发动机公司 制造具有类似于公差叠加的尺寸要求的零部件的方法
CN107256001A (zh) * 2017-05-27 2017-10-17 四川用联信息技术有限公司 改进的衡量制造过程多元质量能力的算法
CN107256000A (zh) * 2017-05-27 2017-10-17 四川用联信息技术有限公司 改进的区域体积比衡量制造过程多元质量能力的算法
CN109545346A (zh) * 2018-11-13 2019-03-29 广州金域医学检验中心有限公司 检测系统的单侧能力评估方法及装置
CN109903822A (zh) * 2019-01-24 2019-06-18 广州金域医学检验中心有限公司 检测结果评估方法及装置、计算机可读存储介质
CN111177658A (zh) * 2020-01-06 2020-05-19 西安石油大学 一种规则几何形状零件尺寸波动的估计方法
CN111222080A (zh) * 2020-01-09 2020-06-02 西安石油大学 一种规则几何形状零件尺寸分布中心的估计方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000127023A (ja) * 1998-10-27 2000-05-09 Hitachi Ltd 半導体の研磨方法及びその製造方法
JP2003114713A (ja) * 2001-10-09 2003-04-18 Nippon Steel Corp 品質不良の原因解析方法
CN101561668A (zh) * 2009-05-15 2009-10-21 江苏网商软件有限责任公司 无油轴套烧结生产线在线监控系统及监控方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000127023A (ja) * 1998-10-27 2000-05-09 Hitachi Ltd 半導体の研磨方法及びその製造方法
JP2003114713A (ja) * 2001-10-09 2003-04-18 Nippon Steel Corp 品質不良の原因解析方法
CN101561668A (zh) * 2009-05-15 2009-10-21 江苏网商软件有限责任公司 无油轴套烧结生产线在线监控系统及监控方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
周祥才: "多品种小批量生产过程能力指数的计算与比较研究", 《常州工学院学报》 *
郑立伟等: "面向多品种小批量生产的工序能力指数计算", 《工业工程》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104133370A (zh) * 2014-07-02 2014-11-05 方志耕 一种基于过程能力的小样本最优质量控制方法
CN107110661B (zh) * 2014-12-05 2019-05-14 赛峰飞机发动机公司 制造具有类似于公差叠加的尺寸要求的零部件的方法
CN107110661A (zh) * 2014-12-05 2017-08-29 赛峰飞机发动机公司 制造具有类似于公差叠加的尺寸要求的零部件的方法
CN104809306A (zh) * 2015-05-13 2015-07-29 梧州学院 基于多质量指标的注塑成型工艺参数优化方法
CN105182903B (zh) * 2015-10-19 2018-03-09 南京航空航天大学 同一加工特征定义方法
CN105182903A (zh) * 2015-10-19 2015-12-23 南京航空航天大学 同一加工特征定义方法
CN107256000A (zh) * 2017-05-27 2017-10-17 四川用联信息技术有限公司 改进的区域体积比衡量制造过程多元质量能力的算法
CN107256001A (zh) * 2017-05-27 2017-10-17 四川用联信息技术有限公司 改进的衡量制造过程多元质量能力的算法
CN109545346A (zh) * 2018-11-13 2019-03-29 广州金域医学检验中心有限公司 检测系统的单侧能力评估方法及装置
CN109545346B (zh) * 2018-11-13 2021-10-19 广州金域医学检验中心有限公司 检测系统的单侧能力评估方法及装置
CN109903822A (zh) * 2019-01-24 2019-06-18 广州金域医学检验中心有限公司 检测结果评估方法及装置、计算机可读存储介质
CN111177658A (zh) * 2020-01-06 2020-05-19 西安石油大学 一种规则几何形状零件尺寸波动的估计方法
CN111222080A (zh) * 2020-01-09 2020-06-02 西安石油大学 一种规则几何形状零件尺寸分布中心的估计方法
CN111222080B (zh) * 2020-01-09 2023-11-28 西安石油大学 一种规则几何形状零件尺寸分布中心的估计方法

Also Published As

Publication number Publication date
CN103760814B (zh) 2016-04-06

Similar Documents

Publication Publication Date Title
CN103760814B (zh) 基于特征的多品种小批量生产零件过程能力指数确定方法
CN101620045B (zh) 基于时间序列的步进应力加速退化试验可靠性评估方法
CN109492683A (zh) 一种针对广域量测电力大数据数据质量的快速在线评估方法
CN104462184B (zh) 一种基于双向抽样组合的大规模数据异常识别方法
CN105676817B (zh) 不同大小样本均值-标准偏差控制图的统计过程控制方法
CN105956788A (zh) 一种输变电工程造价的动态管理控制方法
CN108520267B (zh) 一种基于时空特征的水文遥测数据异常检测方法
CN103954940B (zh) 雷达网基于交叉定位点聚类的集中式与分布式压制干扰鉴别方法
CN111046532B (zh) 一种基于肘形判据的同调机群聚类识别方法
CN104776827B (zh) Gps高程异常数据的粗差探测方法
CN105825288B (zh) 一种消除复杂系统中回归数据共线性问题的优化分析方法
Tian et al. An outliers detection method of time series data for soft sensor modeling
CN116520772A (zh) 基于多任务迁移学习的薄板零件数控加工质量预测方法
CN104376174A (zh) 一种基于线路阻抗比的交流线路参数辨识与修正方法
CN107808209B (zh) 基于加权kNN距离的风电场异常数据辨识方法
CN101780488B (zh) 内螺纹冷挤压加工质量在线预测方法
CN116662829B (zh) 一种场群风机标准功率曲线定义规则和偏差验证方法
CN116611743B (zh) 基于大数据的建筑工程建造质量评价方法
CN108052587B (zh) 基于决策树的大数据分析方法
CN114189313B (zh) 一种电表数据重构方法及装置
CN105183612A (zh) 服务器可用内存异常增长及运行状况的评估方法
CN112884348B (zh) 基于动态贝叶斯网络的航天起爆器生产偏差源诊断方法
Al-Sobky et al. A generalized piecewise regression for transportation models
CN105404736B (zh) 基于多源置信模糊信息的严酷度计算方法
Jie et al. The study for data mining of distribution network based on particle swarm optimization with clustering algorithm method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant