4. Taguchi方法与正交试验设计
在注塑工艺参数优化中,传统做法往往是“把每个参数调到最好”,例如追求熔体温度、保压压力、冷却时间都落在某个理想点。但实际生产中,材料批次波动、环境温湿度变化、模具磨损、设备重复精度等噪声因素始终存在,单点最优解很容易在扰动下失效。Taguchi方法的核心思想不是消除噪声,而是让工艺参数组合对噪声不敏感,即实现鲁棒参数设计。本节围绕质量损失函数、信噪比、内表与外表设计以及正交表选择,结合注塑工艺给出可操作的讲解。
质量损失函数:把“偏离目标”变成可比较的代价
Taguchi提出,质量不是“是否落在规格内”的二元判断,而是产品特性偏离目标值即产生损失。对于望目特性,质量损失函数常写为:
L(y) = k · (y - m)²
k = A₀ / Δ₀²
其中,y为实际特性值,m为目标值,k为损失系数,A₀为在偏差Δ₀处产生的损失。以注塑件关键尺寸为例,若目标值m=50.00 mm,当偏差达到0.10 mm时客户返工损失为8元,则k=8/0.10²=800元/mm²。若某工艺组合下尺寸为50.06 mm,则单件质量损失L=800×0.06²=2.88元。这个数值比“合格/不合格”更细,能区分同样合格但波动不同的方案。
对于望小特性,如翘曲量、缩痕深度,损失函数可简化为L=k·y²;对于望大特性,如熔接痕强度,则用L=k/y²。注塑中多数尺寸和形位特性属于望目或望小,因此后续以望目特性为主展开。
信噪比:同时衡量均值与波动的综合指标
信噪比(SN比)把“均值靠近目标”和“波动小”合并成一个可比较的指标。望目特性的常用形式为:
η = 10·log₁₀( μ² / σ² ) (dB)
其中μ为特性均值,σ为标准差。η越大,说明信号(均值)相对噪声(波动)越强,鲁棒性越好。实际计算时,若每个试验条件下有n次重复,可先求均值与方差,再代入上式。例如某注塑试验条件下测得5个制件尺寸为50.04、50.06、50.05、50.07、50.03 mm,均值μ=50.05 mm,样本标准差s≈0.0158 mm,则η=10·log₁₀(50.05²/0.0158²)≈70.0 dB。若另一条件均值同为50.05 mm,但标准差为0.0316 mm,则η≈64.0 dB,说明后者虽然均值相同,但波动大一倍,鲁棒性明显更差。
望小特性常用η=-10·log₁₀(1/n·Σyᵢ²),望大特性常用η=-10·log₁₀(1/n·Σ1/yᵢ²)。在注塑中,若关注翘曲量,应选望小SN比;若关注熔接强度,应选望大SN比。
内表与外表:把可控因子和噪声因子分开安排
Taguchi方法最显著的结构是“内表—外表”双层设计。内表安排可控因子,如熔体温度、保压压力、保压时间、冷却时间、注射速度等;外表安排噪声因子,如材料批次、环境温度、模具温度波动、设备重复精度等。内表每个试验条件都要在外表若干噪声组合下运行,从而观察该条件在扰动下的表现。
以注塑为例,可控因子取4个,每个2水平:A熔体温度(230/245 ℃)、B保压压力(60/80 MPa)、C保压时间(8/12 s)、D冷却时间(15/20 s)。若用L8(2⁷)正交表安排内表,可得到8个工艺组合。噪声因子取材料批次(新料/回料)、环境温度(20/30 ℃),外表用L4(2³)安排4种噪声组合。这样每个内表条件需做4次试验,总试验数为8×4=32次。若不做内外表,直接全因子需要2⁴×2²=64次,且无法系统评估噪声影响。
正交表选择:用最少的试验看清主效应
正交表的选择取决于因子数、水平数和是否需要考察交互作用。注塑工艺中,若4个因子均为2水平,且暂不考虑交互,L8(2⁷)足够安排4个主效应,剩余列可作为误差估计。若需要考察A×B交互,则需把交互安排到指定列,L8仍可容纳部分交互,但需注意混杂。若因子为3水平,如熔体温度取230/240/250 ℃,则常用L9(3⁴)或L18(2¹×3⁷)。选择原则是:先确定因子与水平,再根据自由度选择正交表,保证试验次数大于总自由度,并优先选用标准表。
以L8(2⁷)为例,前4列安排A、B、C、D,第5列可作为误差列。试验顺序应随机化,避免机器预热、环境变化等系统误差混入因子效应。每个内表条件在外表4种噪声组合下各做一次,共32次试验。若资源有限,也可先做内表8次、外表仅取最恶劣2种组合,共16次,但噪声覆盖会变窄。
注塑案例:从数据到鲁棒参数
假设某注塑件关键尺寸目标为50.00 mm,按上述内表L8、外表L4进行试验。表1给出部分内表条件的汇总结果(示例数据,用于说明计算过程)。
| 内表条件 | A 熔体温度 | B 保压压力 | C 保压时间 | D 冷却时间 | 均值 μ / mm | 标准差 σ / mm | SN比 η / dB |
|---|---|---|---|---|---|---|---|
| 1 | 230 | 60 | 8 | 15 | 50.05 | 0.0158 | 70.0 |
| 2 | 230 | 60 | 12 | 20 | 50.03 | 0.0126 | 72.0 |
| 3 | 230 | 80 | 8 | 20 | 50.02 | 0.0100 | 74.0 |
| 4 | 230 | 80 | 12 | 15 | 50.06 | 0.0200 | 68.0 |
| 5 | 245 | 60 | 8 | 20 | 50.04 | 0.0250 | 66.0 |
| 6 | 245 | 60 | 12 | 15 | 50.07 | 0.0316 | 64.0 |
| 7 | 245 | 80 | 8 | 15 | 50.01 | 0.0080 | 75.9 |
| 8 | 245 | 80 | 12 | 20 | 50.02 | 0.0090 | 74.9 |
从SN比看,条件7的η=75.9 dB最大,条件8次之。条件7为A2B2C1D1,即熔体温度245 ℃、保压压力80 MPa、保压时间8 s、冷却时间15 s。但需注意,条件7的均值50.01 mm略低于目标50.00 mm,而条件3的均值50.02 mm、σ=0.0100 mm,η=74.0 dB,均值更接近目标。若同时考虑均值偏移,可对SN比最大的几个条件做调整:例如在条件7附近微调保压时间至9 s,观察均值是否回到50.00 mm附近,同时保持σ较小。
进一步做因子效应分析:对每个因子在各水平下的SN比求平均。以因子A为例,A1水平(条件1-4)的η均值为(70.0+72.0+74.0+68.0)/4=71.0 dB;A2水平(条件5-8)为(66.0+64.0+75.9+74.9)/4=70.2 dB。A1略优。因子B:B1为(70.0+72.0+66.0+64.0)/4=68.0 dB,B2为(74.0+68.0+75.9+74.9)/4=73.2 dB,B2明显优。因子C:C1为(70.0+74.0+66.0+75.9)/4=71.5 dB,C2为(72.0+68.0+64.0+74.9)/4=69.7 dB,C1略优。因子D:D1为(70.0+68.0+64.0+75.9)/4=69.5 dB,D2为(72.0+74.0+66.0+74.9)/4=71.7 dB,D2略优。综合主效应,较优组合为A1B2C1D2,即230 ℃、80 MPa、8 s、20 s。该组合未出现在原8次试验中,需做验证试验。
实施要点与常见误区
第一,噪声因子必须真实可发生。注塑中材料批次、环境温湿度、模具温度波动、设备液压油温等都可作为外表因子,但不要为了凑外表而引入不可能出现的极端条件。第二,SN比不是唯一指标,望目特性下应同时看均值是否靠近目标;若均值偏离过大,即使σ很小,损失函数仍可能偏高。第三,正交表选择后要检查自由度:4个2水平主效应占4个自由度,L8有7列,剩余3列可估计误差,若交互作用显著则需改用更大表。第四,试验顺序随机化,避免将“先做的条件”与“后做的条件”系统混淆。第五,验证试验不可省略,因为正交表分析得到的是主效应最优组合,未必等于全局最优,且可能存在未考察的交互。
通过质量损失函数把质量转化为经济代价,通过SN比把均值与波动合并评价,通过内表与外表分离可控因子和噪声因子,再借助正交表以较少试验获得主效应信息,注塑工艺参数鲁棒设计就形成了一条可计算、可验证、可复现的路径。最终目标不是找到“最好看”的单点参数,而是找到在材料、环境、设备波动下仍能稳定靠近目标值的参数组合。