SPC·实战笔记
模块十 · 案例·误区·综合实战 · 高阶

案例:CUSUM 提前 3 天发现刀损

同一批缸孔数据,休哈特图第 54 班才报警,CUSUM 第 45 班就跳了。中间隔着 9 个班、整整 3 天。更妙的是,CUSUM 还顺手把"偏移是从第 34 班开始的"给算了出来,追溯范围一下子从两千件缩到七百件。

一、先讲个真事

发动机厂精镗缸孔,φ80.000±0.020 mm,CBN 镗刀,标称寿命 3000 件。检测方式是每班送一件到计量室做三坐标全尺寸,一天三个点。

这个工位一直挂着单值休哈特图,控制限按稳定期数据算的 79.988~80.012。三年下来图很好看,几乎不报警,大家挺满意。

后来我们在同一套数据上并行跑了一张 CUSUM。第 45 班,CUSUM 的下侧累积和冲破决策区间报警。班组长看休哈特图——点子好好地待在限内,离下限还差一截,规格更是宽得很,79.9944 距离下差 79.980 还有 0.0144 mm。谁也不愿意为一张"新图"停机换刀。

我们坚持要求做刀具检查。拆下来看,后刀面磨损带 VB 已经到 0.22 mm,刀尖有明显退让。换刀,孔径回到 80.000。

后续把数据补完整:如果不换刀,休哈特图要到第 54 班才会有点子跌破 LCL。9 个班,3 天。这 3 天里过程一直在往下走,虽然还没出废品,但已经在啃能力余量了——Cpk 从 1.67 掉到 1.20。

老师傅后来跟我说了句实在话:"休哈特图是等你摔下来才喊,CUSUM 是看你脚下打滑就喊。"

二、核心知识点拆解

2.1 CUSUM 到底累积了什么

休哈特图判异只看当前这一个点,前面 20 个点是什么样,它一概不记。这就注定它对小偏移迟钝:偏移 1σ 时,每个点超限的概率仍然很低,你得平均等 44 个点。

CUSUM(累积和)换了个思路:把每个点相对目标值的偏离量攒起来。单个偏离不起眼,攒够了就藏不住。

工程上用的是表格法双侧 CUSUM:

``` Cᵢ⁺ = max[ 0 , xᵢ − (μ₀ + K) + Cᵢ₋₁⁺ ] Cᵢ⁻ = max[ 0 , (μ₀ − K) − xᵢ + Cᵢ₋₁⁻ ]

初值 C₀⁺ = C₀⁻ = 0 当 Cᵢ⁺ > H 或 Cᵢ⁻ > H 时判异 ```

两个参数:

工程上最通用的一组是 k=0.5、h=5,对应 ARL₀ ≈ 465(受控时平均 465 点才误报一次),对 1σ 偏移的 ARL₁ ≈ 10.4。

CUSUM 表格法原理图。左侧一条带波动的数据序列,中心线 μ₀,上下各画一条虚线 μ₀±K(金橙),标注"落在这条带内的偏离被 K 吃掉";右侧对应画 C⁺、C⁻ 两条累积曲线(深蓝),水平红线 H,标出穿越点。图注:"K 管免赔额,H 管报警门槛"。
图 1 · CUSUM 表格法原理图。左侧一条带波动的数据序列,中心线 μ₀,上下各画一条虚线 μ₀±K(金橙),标注"落在这条带内的偏离被 K 吃掉";右侧对应画 C⁺、C⁻ 两条累积曲线(深蓝),水平红线 H,标出穿越点。图注:"K 管免赔额,H 管报警门槛"。

2.2 ARL 对照:什么时候该换武器

选图的唯一硬指标是 ARL(平均运行长度)。受控时 ARL₀ 越大越好(少误报),失控时 ARL₁ 越小越好(快检出)。

偏移 δ(σ 单位)休哈特 I 图(仅 3σ 准则)CUSUM(k=0.5, h=5)EWMA(λ=0.2, L=3)
0(受控)370465500
0.5σ15538.041
1.0σ43.910.411
1.5σ15.05.756.1
2.0σ6.304.014.2
3.0σ2.002.572.6

看最后一行:偏移到 3σ 这种大跳变,休哈特反而更快(2.00 < 2.57)。这不奇怪——CUSUM 要攒,攒需要时间;休哈特只看当下,大跳变一眼就见。

所以结论很明确:CUSUM 和休哈特不是替代关系,是分工关系。小漂移(刀具磨损、砂轮钝化、电极烧损、温度爬升)交给 CUSUM,大跳变(换错料、装夹失误、程序调错)交给休哈特,两张图挂在同一个屏上并行看。

ARL₁ 对比曲线图。横轴偏移量 0~3σ,纵轴 ARL₁(对数坐标),三条曲线:休哈特(灰)、CUSUM(深蓝)、EWMA(金橙)。在 δ≈2.5σ 处标出交叉点,注明"此处之后休哈特反超"。图注:"小漂移用 CUSUM,大跳变用休哈特"。
图 2 · ARL₁ 对比曲线图。横轴偏移量 0~3σ,纵轴 ARL₁(对数坐标),三条曲线:休哈特(灰)、CUSUM(深蓝)、EWMA(金橙)。在 δ≈2.5σ 处标出交叉点,注明"此处之后休哈特反超"。图注:"小漂移用 CUSUM,大跳变用休哈特"。

2.3 报警之后:偏移从哪一班开始的

这是 CUSUM 最被低估的一个功能。报警只告诉你"出事了",追溯范围却决定了要拦多少货。CUSUM 能直接算出偏移起点:

``` 偏移起始点 = i − N⁻ + 1 ```

N⁻ 是报警时刻 C⁻ 连续保持非零的点数(正侧同理用 N⁺)。偏移后的实际均值也能估:

``` μ̂ = μ₀ − K − C⁻ / N⁻ (下侧) μ̂ = μ₀ + K + C⁺ / N⁺ (上侧) ```

有了起点,追溯就从"全批扣留"变成"定点扣留",省下来的钱通常比 SPC 系统本身还贵。

三、车间真实案例 / 计算过程

(1)基线参数

稳定期 30 个班的数据算得:

``` μ₀ = 80.0000 mm ,σ̂ = MR̄/1.128 = 0.0040 mm 休哈特 I 图:UCL = 80.0120 ,CL = 80.0000 ,LCL = 79.9880 规格:79.980 ~ 80.020 ,Cpk 初始 = 0.020/(3×0.004) = 1.67 ```

CUSUM 参数取 k=0.5、h=5:

``` K = 0.5 × 0.0040 = 0.0020 mm → μ₀ − K = 79.9980 H = 5 × 0.0040 = 0.0200 mm ```

(2)第 31–45 班逐点手算(下侧 C⁻)

计算式:C⁻ᵢ = max[0, 79.9980 − xᵢ + C⁻ᵢ₋₁]

班次 i孔径 xᵢ (mm)79.9980 − xᵢC⁻ᵢ休哈特是否报警
3179.9992−0.00120
3279.9975+0.00050.0005
3379.9986−0.00060
3479.9968+0.00120.0012
3579.9973+0.00070.0019
3679.9981−0.00010.0018
3779.9962+0.00180.0036
3879.9970+0.00100.0046
3979.9958+0.00220.0068
4079.9966+0.00140.0082
4179.9955+0.00250.0107
4279.9963+0.00170.0124
4379.9948+0.00320.0156
4479.9957+0.00230.0179
4579.9944+0.00360.0215 > H=0.0200 → 报警

(3)定位偏移起点与偏移量

C⁻ 从第 34 班起连续非零,到第 45 班报警,所以:

``` N⁻ = 45 − 34 + 1 = 12 偏移起始班次 = 45 − 12 + 1 = 34 μ̂ = 80.0000 − 0.0020 − 0.0215/12 = 79.9980 − 0.00179 = 79.9962 mm ```

拿第 34~45 班的实测值直接求平均验证:

``` 实测均值 = 79.99621 mm ```

估计值 79.9962 和实测 79.99621 几乎重合,说明这个估计方法是靠谱的。偏移量 80.0000 − 79.9962 = 0.0038 mm,正好约 1σ。

(4)休哈特什么时候才报

不换刀继续跑,孔径按约 0.0008 mm/班 的速率下滑:

班次464748495051525354
孔径79.994079.993279.992579.991879.991079.990379.989679.988879.9879

第 54 班首次跌破 LCL=79.9880,休哈特报警。

``` CUSUM 报警:第 45 班 休哈特报警:第 54 班 提前量 = 9 个班 = 3 天(三班制) ```

(5)这 3 天值多少钱

时间点孔径Cpk_L =(x−79.980)/(3σ)状态
第 30 班(稳定)80.00001.67能力充足
第 45 班(CUSUM 报警)79.99441.20尚有余量,无废品
第 54 班(休哈特报警)79.98790.66余量耗尽,再跑 10 班必出废品

追溯这一块的收益更直观:没有起点定位,得从上一次刀具更换(第 20 班)开始全部扣留,约 2100 件;有了"第 34 班"这个起点,只扣第 34 班之后的约 750 件。少复检 1300 多件。

双联对比图。上图休哈特 I 图:横轴 30–56 班,CL/UCL/LCL 三线,点子从第 34 班起缓慢下滑,第 54 班第一个红点标"休哈特首次报警";下图 CUSUM C⁻ 图:横轴对齐,H=0.0200 红色水平线,C⁻ 曲线从第 34 班起爬升,第 45 班穿越 H 处标金橙星号"CUSUM 报警"。两图之间…
图 3 · 双联对比图。上图休哈特 I 图:横轴 30–56 班,CL/UCL/LCL 三线,点子从第 34 班起缓慢下滑,第 54 班第一个红点标"休哈特首次报警";下图 CUSUM C⁻ 图:横轴对齐,H=0.0200 红色水平线,C⁻ 曲线从第 34 班起爬升,第 45 班穿越 H 处标金橙星号"CUSUM 报警"。两图之间用金橙双向箭头连接第 45 与第 54 班,中间大字「提前 9 班 = 3 天」。

四、常见误区提醒

  1. k 和 h 拍脑袋乱设。 见过有人把 h 设成 2,图天天报警,一个月就没人看了;也见过 h 设 10,等于没装。k 必须按你要检出的偏移定(k = δ/2),h 按可接受的误报频率定(常规 h=4~5),设定依据要写进控制计划,不能凭手感。
  1. 上了 CUSUM 就把休哈特撤了。 从 ARL 表能看到,δ≥2.5σ 之后休哈特反而更快,而换错料、装夹失误这类大跳变恰恰就是大偏移。两张图并行,各管一段,撤掉任何一张都是给自己留缺口。
  1. μ₀ 和 σ 直接用全历史数据算。 这是最常见的技术性错误。如果历史数据里本来就含着漂移,算出来的 σ 会被撑大,K 和 H 跟着变大,CUSUM 的灵敏度当场归零。基线必须取自验证过受控的那一段,σ 用 MR̄/1.128 而不是整体标准差。

五、小结 & 互动提问

CUSUM 不玄,就两个公式、两个参数,Excel 拉一列就能算。它值钱的地方在两处:对 1σ 上下的小漂移比休哈特快四倍,报警后还能反推偏移起点,把追溯成本压下来。

咱做质量的,最怕的从来不是"突然坏了",而是"一直在慢慢坏,还没人喊"。CUSUM 就是专门对付后面这种的。

互动提问:你车间有哪些参数是"慢性病"型的——刀具磨损、砂轮钝化、镀液浓度衰减、模具型腔磨损?挑一个,把历史数据倒进 Excel 按 k=0.5、h=5 跑一遍,看看 CUSUM 会比现在的图提前多少个点报警?

如需配套互动练习、培训课件或软件试用,可搜索「东方大易」或访问大易官网 www.eastdayi.com(苏州东方大易管理咨询有限公司,电话 0512-65561989)联系客服获取。

【来源标注】本篇依据:GB/T 17989 系列《控制图》(累积和控制图的构造、参数选择与 ARL 评价);ISO 11462 系列(SPC 实施指南,含控制图选型与实施要素);能力指数计算参见 ISO 22514 系列。可配套互动页:iso-spc-tutorials/cusum_control_chart_guide.html(CUSUM 累积和控制图互动教程)、iso-spc-tutorials/ewma_control_chart_guide.html(EWMA 对照)、iso-spc-tutorials/special_control_chart_guide.html(特殊控制图);知识点单页见 iso-spc-kb/ 目录(KP101–KP310)。
苏州东方大易(大易培训中心)
专业汽车行业标准培训与咨询 · IATF 16949 / APQP / FMEA / SPC / CP / PPAP / VDA 6.3 内审员
官网 www.eastdayi.com | 咨询电话 0512-65561969 / 0512-65561989
相关文章推荐