直方图及工序能力分析

1、Nov, 20,2005zhong xin1直方图的概念直方图的概念直方图的作图步骤直方图的作图步骤直方图常见类型分析直方图常见类型分析直方图直方图 VS.VS.规格的比较分析规格的比较分析直方图使用要点直方图使用要点Nov, 20,2005zhong xin2直方图的概念直方图的概念直方图 (Histogram): 直方图: SPC七个常用品质管理工具之一. 是通过对数据加工整理，从而分析和掌握数据的分布状况和估算工序能力的一种方法。 从总体中随机抽取样本， 整理从样本中获得的数据，以矩形的形式得到的频率分布图形 矩形底是单元长度 矩形面积与频率成正比。 作用: 了解生产过程的状态及质量特性

2、值分布的情况 判断工序过程能否滿足质量要求 显示各种数值出现的相对频率，揭示数据的中心、散布及形状，快速阐明数据的潜在分布 为预测过程提供有用信息Nov, 20,2005zhong xin3 数据间隔 - 组距 h 组距内的工件数 - 频数 n n/N - 频率(N为样本容量), - 频率密度或分布密度 以组距 h为横坐标, 以 n 或 g 为纵坐标做图, 即可得到实际分布图直方图。 作直方图三步骤:确定样本容量N及实测数据确定组数 k、组距 h 和组界画直方图及计算分析hNng1直方图的作图步骤直方图的作图步骤Nov, 20,2005zhong xin4 一般样本容量N应在50个以上,较容易

3、显示出整体数据分布的情况。确定样本容量确定样本容量N N及实测数据及实测数据Nov, 20,2005zhong xin5确定组数确定组数 k k、组距、组距 h h 和组界和组界适当分组 ( k ) 组数太少, 会掩盖组內的变化情况组数太多, 会造影响数据分布的明显性, 难以看清分布的状况组数与样本容量 N 有关, 可按下表推荐值, 或按式 k=1+3.31*lgN 确定, 本例取k=7。组距 （h）最大值 Xmax=21um, 最小值 Xmin=11.5um组距 h=(Xmax-Xmin)/k =(21-11.5)/7 =1.36um圆整为 h=1.5um (按测量量具最小分辨值的整数倍进行

4、圆整,本例量具最小分辨值为0.5um)组界各组组界最好选在测量数据最后一位尾数的1/2处, 以免数据落在组界上, 如测量值尾数为0.5m, 组界应取在0.25um上。各组组界为 Xmin + (j-1)h h/2 各组的中值为 Xmin + (j-1)h (j = 1,2,3,k)统计频数、频率n/N和分布密度gNov, 20,2005zhong xin6根据表中有关数据画出直方图用分布密度作纵坐标有两个好处, 可避免因样本容量和组距不同而使分布图形状不同 每个小矩形面积恰好等于该组距内的工件出现的频率计算出样本均值Xbar和标准差s。Xbar表示样本中心,s反映样本分散的程度。分析直方图类型

5、,画规格线,分析工序状态.画直方图及计算分析画直方图及计算分析Minitab Graph HistogramNov, 20,2005zhong xin7直方图常见类型分析直方图常见类型分析11Nov, 20,2005zhong xin8直方图常见类型分析直方图常见类型分析22Nov, 20,2005zhong xin9直方图直方图 VS.VS.规格的比较分析规格的比较分析11Nov, 20,2005zhong xin10直方图直方图 VS.VS.规格的比较分析规格的比较分析22Nov, 20,2005zhong xin11有规格线的直方图可用来比较过程与要求。此时应确认直方图具有合适的比例。直

6、方图不应该单独使用。通常在它之前先构造一张链图或控制图。 因为直方图中的数据不是按时间顺序给出的, 经常掩盖了失控的事实评价直方图的模式以确定是否能够检测到任何形式的变化。比较不同时间段内的直方图。 观察直方图从一个时间段到下一个时间段的模式变化,对寻找过程改进的方法非常有用。根据数据来源的不同, 分别绘制直方图, 对数据分层。 例如, 对描述金属棒直径的直方图来说, 可单独作由不同供应商的原材料制造的直径直方图, 或由不同操作工或机器生产的棒径的直方图。 有时,这可能会揭示控制图都不能检测到的事情。直方图使用要点直方图使用要点Nov, 20,2005zhong xin12正态分布及其性质正态

7、分布及其性质非正态分布非正态分布正态性判定正态性判定非正态分布数据处理非正态分布数据处理Nov, 20,2005zhong xin13正态分布正态分布概率论己经证明,相互独立的大量微小随机变量,其总和的分布服从正态分布。大量的试验表明,在一次调整好的机床上,连续加工一批工件,若无变值系统误差的影响,误差是由一些相互独立的随机因素引起的,这些因素中又无明显优势者,其参数是服从正态分布的。正态分布的概率分布密度函数表达式为 式中 g(x) - 分布的概率密度 - 总体均值 - 总体标准差(均方差)正态分布曲线即高斯曲线如图所示exxg222)(21)(2)(1111NiiNiixNxNExcel

8、function=average () =stdevp ()Nov, 20,2005zhong xin14正态分布的性质正态分布的性质11正态分布曲线有以下几个主要特点 曲线呈钟形, 对称于平均值, 即g(+a) = g(-a) 在不变的情况下,变化只能影响分布曲线的位置而不影响分布曲线的形状和分散范围, 的变化是由常值系统误差引起的; 若不变而变化,曲线的位置不变, 但形状和分散的范围发生变化,值对正态分布的影响值对正态分布的影响Nov, 20,2005zhong xin15正态分布的性质正态分布的性质22正态分布曲线有以下几个主要特点 (续） 当x时, g(x) 0, 即g(x)以x轴为渐

9、近线; 在x=时, g(x)有最大值, 即 在x=处, 曲线有二拐点, 在二点之间曲线向上凸, 在二点之外曲线下凹 为 区间的面积, 等于该区间的概率, =0, =1 的正态分布为标准正态分布。任何不同的和的正态分布,都可以通过令 z= (x-)/进行坐标变换变成标准正态分布 标准正态分布从 0 到 z 区间的概率, 即该区间内曲线与横坐标所包含的面积。 不同 z 值的(z), 可查表求得21)(max ggaxxdxxg)()(axx1)(dxxg2221)(zezgdzedzzgzzzz020221)()(Nov, 20,2005zhong xin16正态分布的性质正态分布的性质33dze

10、dzzgzzzz020221)()(当 x=3( ) 时, 2(3)= 2*0.49865 = 99.73%, 即只有0.27% 的概率落在该范围之外, 可忽略不计, 因此一般取正态分布的分散范围为 3 (6)3xzNov, 20,2005zhong xin17正态分布的性质正态分布的性质44 如果加工后工件参数服从正态分布, 就可以利用正态分布的一些特点来分析加工误差。 由于正态分布的和是求不出来的,所以一般通过它的随机样本的均值 xbar 和标准差 s 来估计, 即2)(11111xxnsxnxniinii Excel functionxbar=average ()s=stdev ()No

11、v, 20,2005zhong xin18非正态分布非正态分布11工件加工后其参数服从正态分布是有条件的,并非所有情况都服从正态分布。 如在加工过程中存在较明显的变值系统误差 将两次调整下加工的工件混在一起,其参数分布将呈双峰形 (a) 如砂轮或刀具磨损显著和热平衡前的刀具热变形等,会形成平顶分布(b)和不对称分布 加工轴时偏向左,加工孔时偏向右(c) 用试切法加工轴和孔时,操作者主观上存在着宁可返修也不报废的倾向性,也往往出现不对称分布,加工轴时偏向右,加工孔时偏向左 工件的对称度、锥度、平行度、垂直度等误差是没有负值的。尽管其仍服从正态分布,由于其负值部分叠加到正值部分,就会出现(d)所示

12、的正值分布(差数模分布) 其他像径向跳动、端面跳动等误差,一般不考虑正负号其分布亦为正值分布(图(e)， (瑞利分布) 如寿命, 故障率等 正态分布的特性值也有可能因数据收集上的问题得到非正态分布结果。 因工序不稳定、组内变动混乱、工序外部原因引起组间变动、数据输入错误等Nov, 20,2005zhong xin19非正态分布非正态分布22Nov, 20,2005zhong xin20非正态分布非正态分布33非正态分布的分散范围, 不再是6,而是1/K*6 。K称为相对分布系数, 表示所研究的数据分布的不同分散性质, 即分布曲线的不同形状。如果以正态分布曲线作比较依据 (即正态曲线的K=l),

13、 各种不同分布的K值可参考下表。表中的为相对不对称系数Nov, 20,2005zhong xin21正态性判定正态性判定确定数据是否服从正态分布的过程 直方图: 绘制一个直方图。如果直方图和钟形相差很大,则拒绝正态性 离群数:找出离群数, 如果有离群数, 则可能拒绝正态性只有一个离群数,可能是一个错误, 或是偶然变异的结果。但是要仔细,因为即使只有一个离群数,也能对结果产生很大的影响.当识别其为一个错误时, 应对其进行修正或删除.但有些不是错误. 考察数据时, 可通过对包含离群数和不包含离群数的数据的分析来研究离群数的影响 正态分位数图:过程略1.31.21.11.00.90.80.70.61

14、050MFTAAFrequency1.41.31.21.11.00.90.80.70.6MFTAANov, 20,2005zhong xin22正态性判定正态性判定用 Minitab 的 Normality Test 功能 Stat Basic Statistics Normality Test In Variables, enter the columns containing the measurement data. If you like, use either of the dialog box options, then click OKP-Value: 0.005A-Square

15、d: 1.144Anderson-Darling Normality TestN: 50StDev: 0.173191Average: 0.8234331.381.281.181.080.980.880.780.680.58.999.99.95.80.50.20.05.01.001ProbabilityMFTAANormal Probability PlotP-Value: 0.459A-Squared: 0.351Anderson-Darling Normality TestN: 62StDev: 0.0085159Average: -0.00006450.010.00-0.01.999.9

16、9.95.80.50.20.05.01.001Probability(H0)ParallelismArmACAMn40Normal Probability PlotNov, 20,2005zhong xin23非正态分布非正态分布数据数据处理处理SPC, 6-Sigma 等使用的统计工具一般是以数据正态分布为基础的。若参数远离正态分布，则将上述计算方法所得结果用于工序分析将不能合理的反映工序的真实情况。若数据呈非正态分布，可先把数据转换为正态分布，或者用威布尔分布模型，可得更确切结果。威布尔分布包括指数分布和瑞利分布。使用 “Box-Cox transformation”,数据必须为正值。处理

17、方法： 收集数据后，用Normality Test检验是否正态分布 如果不是正态分布 数据转换为正态分布(Box cox变换等) 否则,利用其Data特性适合分布,统计性解释。Nov, 20,2005zhong xin24非正态分布非正态分布数据数据处理处理 若参数不符合正态分布，则将上述计算方法所得结果用于工序分析不是很合理。 作工序分析时，若数据呈非正态分布，可先把数据转换为正态分布、或者用威布尔分布模型，可得更确切结果。Nov, 20,2005zhong xin25判别加工误差的性质判别加工误差的性质确定工序能力确定工序能力估算不合格率估算不合格率Nov, 20,2005zhong xi

18、n26判别加工误差的性质判别加工误差的性质Nov, 20,2005zhong xin27判别加工误差的性质判别加工误差的性质 在工件的加工过程中, 如果没有变值系统误差, 只有随机误差时, 其参数分布应该服从正态分布。 如果实际分布与正态分布出入很大, 可根据其实际分布形状初步判断变值系统误差是什么类型。 常值系统误差可以根据分布的位置相对公差带位置来判断。=参照直方图及正态性判定部分。Nov, 20,2005zhong xin28评估工序能力评估工序能力Nov, 20,2005zhong xin29预测工序与规格要求的符合程度。帮助产品设计开发人员选择或修改工序，防止设计与制造脱节。辅助设立

19、工序控制的合适抽样区间。新机器或旧机器维修时的评估。不同供应商质量评比。当不同工序间相关时，可以提供工序规划的参考。降低工序的变异性。分辨操作人员的技能。等等 工序能力分析的用途工序能力分析的用途Nov, 20,2005zhong xin30工序能力工序能力工序能力：某工序处于控制状态(稳定状态)下所能达到的加工要求的能力。可用该工序参数分散范围来表示其工艺能力, 对于正态分布（或近似正态分布）, 其分散范围为6, 取6为其工序能力。产品质量受生产过程状态的影响, 而生产过程状态受到“5MlE”的影响，其综合效果反映了产品质量特性值的分布情况, 即和的状况。当“5MlE”受到完善的管理和控制时

20、, 通常已经消除了系统性因素的影响, 仅存在偶然性因素的影响。这时, 综合影响效果的质量特性值的概率分布, 反映了工序的实际加工能力。如前所述3,即 6范围内包含了99.73的质量特性值, 所以可以将工序能力定量表示为 6。在计算工序能力时，首先应对5M1E等条件充分标准化，并使作业活动处于受控状态。 6越大，則工序的实际精度越差，工序能力越小；6越小，則工序的实际精度越高，工序能力越大.Nov, 20,2005zhong xin31工序能力指数：就是表示工序能力满足产品质量标准的程度的评价指标。所谓产品质量标准，通常指产品规格、工艺规范、公差等。 只有在工序受控状态下获得的工序能力指数才对工

21、只有在工序受控状态下获得的工序能力指数才对工序控制有指导意义序控制有指导意义!工序能力指数一般用符号Cp 表示，则 Cp=T/6 T - 公差， - 总体标准差（或用样本标准差s)只有规格上限时，只有规格下限时，当有双规格限时，3USLCPU3LSLCPLTkMLSLUSLTTTkCCppk2,3262)1 (工序能力指数工序能力指数Nov, 20,2005zhong xin32Excel functionxbar= average (), s= stdev ()Cpk=IF(USL=,IF(LSL=,N/A,(mean-LSL)/3/s),IF(LSL=,(USL-mean)/3/s,MIN

22、(mean-LSL)/3/s,(USL-mean)/3/s)工序能力指数的计算工序能力指数的计算工序能力指数的计算在实际中对动态总体（生产过程）进行随机抽样，统计计算所收集的数据得到样本统计量，即样本的平均值 xbar和标准差s，用xbar和s去估计和。同时，可将 CPU 、 CPL 、Cpk 统一为：)3,3min() C , (Cmin CpkPLPUsLSLxsxUSLNov, 20,2005zhong xin33Cpk DPPM1.00.13 0.27 %1.10.05 0.10 %1.20.02 0.03 % 1.348.1 96.2 ppm1.413.4 26.7 ppm1.5 3

23、.4 6.8 ppm1.6 794 1589 ppb1.7 170 340 ppb1.8 33 67 ppb1.9 6 12 ppb2.0 1 2 ppb工序能力指数与工序能力指数与DPPMDPPMNov, 20,2005zhong xin34工序能力指数的评定工序能力指数的评定Nov, 20,2005zhong xin35sTCpk62提高工序能力指数的途径提高工序能力指数的途径提高工序能力指数的途径 影响工序能力指数的三個变量：T，s 调整工序加工的分布中心，减少偏移量 找出偏移规律; 调整设备; 刀具等加工定位装置; 改变操作倾向性方式; 精确量规 提高工序能力，减少分散程度 s 改进工

24、艺方法,优化工艺参数,补充中间工序,应用新材料、工艺、技术;环境改造;提高工具精度; 强化现场管理, 开展 QCC 改善活动; 缩短进货周期; 技术培训等 修订公差范围 T 当确信放宽公差范围不致影响产品质量时，有必要修订不切实际的现有公差范围Nov, 20,2005zhong xin36工序能力与工序性能工序能力与工序性能工序能力工序能力 工序性能工序性能 . WithinpLSLUSLC6OverallpLSLUSLP6),min(33puplpkWithinpuWithinplCCCUSLCLSLC),min(33puplpkOverallpuOverallplPPPUSLPLSLP主要

25、差异的在于 Within 和 Overall 的不同 (See appendix )虚线代表全部过程的变差:Overall每个小曲线代表组内变差, 或一个子组（某一时刻）的变差: WithinNov, 20,2005zhong xin37工序能力与工序性能工序能力与工序性能Nov, 20,2005zhong xin38工序能力指数的替代工序能力指数的替代 Cpm（k) 同时度量了相对于目标值 M 的工序分布和工序的中心趋势。工序平均值 = 目标值，则 Cpm =Cpk。工序平均值偏离目标值，则 Cpm Cpk。 Cpm 的偏差计算不是将数据与平均值相比，而是与目标值相比。因而偏离目标的测量值将

26、增加Cpm偏差Cpm与Cpk的差值代表了调整过程中心可带来的潜在收益。 targetSpec:1)()3,3min()3,3min(6)(661222222222MnMxLSLxxUSLLSLUSLCLSLUSLMxsLSLUSLLSLUSLCniiCpmCpmCpmpmkCpmpmNov, 20,2005zhong xin39工序能力指数的替代工序能力指数的替代0.22Nov, 20,2005zhong xin40工序能力指数总结工序能力指数总结Pp(k) vs Cp(k)，Ppm(k) vs Cpm(k)，计算公式相同;但Pp(k)和Ppm(k)反映长期工序能力，其标准偏差计算方法不同。主

27、要差异的在于 Within 和 Overall 的不同Nov, 20,2005zhong xin41估算不合格率估算不合格率Nov, 20,2005zhong xin42如果加工后, 产品参数服从正态分布, 且不存在常值系统误差, 只要满足T6,基本上就不会产生不合格品但实际上不存在常值系统误差似乎不可能。 如, 调整机床时不可能将尺寸分布中心调整的与公差带中心完全重合, 其偏离值设为。 此时不产生不合格品的条件为T6+2 如果6T 6+2时,可增加一个方向与相反的常值系统误差(如重新调整机床, 或重新选定刀具的尺寸), 来减小或消除,使分布中心接近或重合于公差带中心, 以保证不产生不合格品

28、如果T Quality Tools Capability Analysis (Normal)Nov, 20,2005zhong xin48综合应用综合应用实例实例 1-21-20.4450.4350.4250.4150.4050.395TargetUSLLSLProcess Capability Analysis for MN40 Bearing HeightPPM TotalPPM USLPPM USLPPM USLPPM Quality Tools Capability Analysis (Weibull)Nov, 20,2005zhong xin50综合应用综合应用实例实例 2-22-2

29、0.0160.0120.0080.0040.000USLUSLProcess Capability Analysis for MN40 Bearing ParallelismCalculations Based on Weibull Distribution ModelPPM TotalPPM USLPPM USLPPM Control Charts Box-Cox Transformation543210-1-2-3-4-5908070605040302095% Confidence IntervalStDevLambdaLast Iteration Info 20.668 20.662 2

30、0.663 -0.506 -0.562 -0.618StDevLambdaUpEstLowBox-Cox Plot for MFTAANov, 20,2005zhong xin52综合应用综合应用实例实例 3-23-2使用 Capability Analysis (Normal Distribution) with the optional Box-Cox power transformation 进行分析。Nov, 20,2005zhong xin53综合应用综合应用实例实例 3-33-30.150.140.130.120.110.100.090.080.07LSL*LSL*Process

31、Capability Analysis for MFTAABox-Cox Transformation, With Lambda = -0.5PPM TotalPPM LSL*PPM TotalPPM USLPPM LSLCpmCpkCPLCPUCpStDev* (Within)StDev (Within)Sample NMean*MeanLSL*LSLTarget*TargetUSL*USL2423.03 *2423.03 0.00 * 0.00 *0.940.94 * * 0.013423.159490 0.109886.8584 0.147446.0000 * * * *Expected PerformanceObserved PerformancePotential (Within) CapabilityProcess DataNov, 20,2005zhong xin54Nov, 20,2005zhong xin55Estimating Estimating WithinW