田口参数实验设计

《田口参数实验设计》由会员分享，可在线阅读，更多相关《田口参数实验设计（14页珍藏版）》请在装配图网上搜索。

1、教学案例一：田口参数实验设计1 田口方法源起实验设计是以概率论与数理统计为理论基础，经济地、科学地制定实验方案以便对实验 数据进行有效的统计分析的数学理论和方法。其基本思想是英国统计学家 R. A. Fisher 在进 行农田实验时提出的。他在实验中发现，环境条件难于严格控制，随机误差不可忽视，故提 出对实验方案必须作合理的安排，使实验数据有合适的数学模型，以减少随机误差的影响， 从而提高实验结果的精度和可靠度，这就是实验设计的基本思想。在三十、四十年代，英、美、苏等国对实验设计法进行了进一步研究，并将其逐步推广 到工业生产领域中，在冶金、建筑、纺织、机械、医药等行业都有所应用。二战期间，英美

2、 等国在工业试验中采用实验设计法取得了显著效果。战后，日本将其作为管理技术之一从英 美引进，对其经济复苏起了促进作用。今天，实验设计已成为日本企业界人士、工程技术人 员、研究人员和管理人员必备的一种通用技术。实验计划法最早是由日本田口玄一(G. Taguch)博士将其应用到工业界而一举成名的。 五十年代，田口玄一博士借鉴实验设计法提出了信噪比实验设计，并逐步发展为以质量损失 函数、三次设计为基本思想的田口方法。田口博士最早出书介绍他的理论时用的就是“实验 计划法一DOE ”，所以一般人惯以实验计划法或DOE来称之。但随着在日本产业界应用的普 及，案例与经验的累积，田口博士的理论和工具日渐完备，

3、整个田口的这套方法在日本产业 专家学者的努力之下，早已脱离其原始风貌，展现出更新更好的体系化内容。日本以质量工 程(Quality Enginerring)称之。但是，严格来讲，田口方法和DOE是不同的东西。田口方 法重视各产业的技术，着重快速找到在最低成本时的最佳质量。 DOE 则重视统计技术，着 重符合数学的严谨性。虽然学术界普遍认为田口方法缺少统计的严格性，但该方法还是以其 简单实用性广为工业界所应用和推广。先进国家对田口方法越来越重视，并且也已经取得了 很好的效果。该方法广泛应用于研发、技术改善、质量提升等部门。八十年代，田口方法进入美国，得到了普遍关注。如今，实验设计技术的应用领域已

4、经 突破了传统的工业过程改进和产品设计范畴，广泛地渗透到商业布局、商品陈列、广告设计 及产品包装的应用之中。我国在六十年代就曾对实验设计进行了研究和推广，八十年代又引 入了田口方法，取得了一定成效。但实验设计作为一种质量改进的有力武器，还尚未发挥它 的全部威力。2 田口方法基本思想和研究内容与传统的质量定义不同，田口玄一博士将产品的质量定义为：产品出厂后避免对社会造 成损失的特性，可用“质量损失”来对产品质量进行定量描述。质量损失是指产品出厂后“给 社会带来的损失”，包括：直接损失，如空气污染、噪声污染等；间接损失，如顾客对产品 的不满意以及由此而导致的市场损失、销售损失等。质量特性值偏离目标

5、值越大，损失越大， 即质量越差，反之，质量就越好。对待偏差问题，传统的方法是通过产品检测剔除超差部分 或严格控制材料、工艺以缩小偏差。这些方法一方面很不经济，而且有时技术上也难以实现。 田口方法是调整设计参数，使产品的功能、性能对偏差的起因不敏感，以提高产品自身的抗 干扰能力。为了定量描述，产品质量损失，田口提出了“质量损失函数”的概念，并以信噪 比来衡量设计参数的稳健程度。由此可见，田口方法是一种聚焦于最小化过程变异或使产品、过程对环境变异最不敏感 的实验设计方法。该方法是一种能设计出环境多变条件下能够稳定和优化操作的高效方法。一般而言，任何一个质量特性值Y在生产过程中均受很多因素的影响，田

6、口玄一博士将 影响质量特性的因素分为输入变量W、可控变量X和不可控变量Z，如图1所示。输入变 量非设计参数，可控变量是田口方法的设计对象，所谓可控变量，即可以调整和控制的参数， 这种变量通常称为信号因子。不可控变量，顾名思义，即不可控制的变量，也称为噪音因子 （Noise Factors），就是使质量特性偏离目标值的因素。田口玄一博士将噪音因子分为三类： 即外部噪音，如温度、湿度、灰尘等；内部噪音，如劣化等；产品间噪音，如制造缺失等。图 1 影响质量特性的关键因素解决的对策可由生产线外（Off li ne）质量控制与线上（On Line）质量控制两种。所谓 线外控制，即产品设计阶段和制造设计阶

7、段的质量控制活动，通过实验设计，保证产品最佳 化和制造过程最佳化（主要是工艺参数的最佳化设计）。线外质量控制可以应用正交表、信 噪比（S/N）和损失函数来达成，强调有效率的实验和仿真，以减少变异。所谓线上质量控 制，是实际生产阶段的质量控制活动。田口式质量工程较关心线外质量控制，以降低成本、 提供最佳质量为目标；对于线上质量控制则以稳定制造过程为目标。田口方法的基本原理是通过控制可控因素的水平和配合，使产品和工艺对噪声因素的敏 感程度降低，从而使噪声因素对产品质量的影响作用减少和消除，以实现提高和稳定产品质 量的目的。田口玄一提出的“三次设计法”，即分三个阶段对产品质量进行优化：（ 1） 系统

8、设计：应用科学理论和工程知识对产品功能原型进行设计开发，这阶段完 成了产品的配置和功能属性；（2） 参数设计：在系统结构确定后进行参数设计。这一阶段以产品性能优化为目标 确定产品参数水平及配置，使工程设计对干扰源的敏感性最低；3） 容差设计：在参数确定的基础上，进一步确定这些参数的容差。系统设计、参数、容差设计等这三方面的内容构成田口方法的“线外质量控制”，田口 线外质量控制、质量损失函数和田口线上质量控制就构成了田口质量工程，如图 2 所示。一 般所讲田口参数设计是指田口线外质量控制，本章讨论的田口实验设计即为田口线外质量控 制。在田口线外质量控制中，参数设计是线外质量控制的核心，它通过实验

9、优化方法确定系 统各参数的最优组合，使产品对环境条件和其他噪声因素的敏感性降低。最终效果是在不提 高产品成本甚至降低成本的基础上使产品质量损失最小，可见，参数设计是获得高质量产品 的关键，也是田口方法的中心内容。系统设计是线外质量控制的基础和前提，容差设计是对 系统设计和参数设计的完善与提高，质量水平评价是对田口线外质量控制的效果评价与分 析。图 2 田口质量控制体系3 田口质量损失函数田口对产品质量提出了一个新概念，他认为：质量就是产品上市后给于社会的损失。一 般，一个产品的成本分为两个主要部分：销售前成本和出售后成本，前者是指制造成本，后 者是指产品销售给用户后由于产品质量的损失(质量特性

10、偏离目标值)所需的费用，这就是 上述产品质量定义中的“给予社会的损失”对此中损失，田口提出用质量损失函数来度量。 为了描述产品的质量损失，引入了以下几种类型质量特性的损失函数。1. 望目特性的质量损失函数望目特性质量损失函数适用于产品的输出特性y有一个确定的目标值y(通常不为零)，并且质量损失在目标值的两侧呈对称分布，如图3所示，这种质量特征称为望目特性。则质量损失函数为：L(y)二 K(y - y。)2(1)其中K是不依赖于y的常数，称为质量损失系数。若y离y0越近，则L(y)值越小，表明该项设计的质量损失小，功能质量好。图 3 功能质量损失函数式(1)说明，由于功能波动所造成的损失与偏离目

11、标值yO的偏差平方成正比。这也可以说 明，不仅不合格产品会造成损失，即使合格产品也会造成损失。输出特性值偏离目标值越远， 造成的损失越大。这就是田口玄一对产品质量概念的一个观点。由于产品的质量特性y表现为随机性，所以L(y)亦为随机变量，故有必要取L(y)的期望 值作为评定产品的质量水平。设有N件产品，若质量特性的N个测试值为儿，丁2,，yN， 则其质量损失可近似表示为：L(y)二 K -1 迓(y - y )2(2)N i 0L i=1称L(y)为这N件产品的平均质量损失。质量损失系数K的确定可以有两种方法确定，一种 是根据功能界限和相应的损失来确定；另一种是根据容差 y和相应的损失来确定。

12、2. 望小特性的质量损失函数有些产品的质量特征是：不取负值，越小越好，目标值为零；当其输出特性值增大时， 其性能逐渐变差，质量损失逐渐变大。这种质量特征称为望小特性。如计算机的响应时间、 汽车的污染、电子线路的电流损失、加工误差等，都属于这类的质量特性。这种情况下的质 量损失函数可由式1中令y0=0得：L(y)=Ky2(3)如图4所示，因为输出特性y只能取正值，故上式只取一侧的损失函数。另外，其质量损失 系数K也可以由功能界限和相应的损失，或者容差界限 y和相应的损失来确定。L(y)为随 机变量时，产品的功能质量水平用期望值反映。对于N件产品，若望小特性的测试值为yl, y2， ，yN，则平均

13、质量损失为：4)L( y) = ki=1图 4 望小特性的质量损失函数3. 望大特性的质量损失函数 有些产品的质量特性是：不取负值，越大越好，零值最差；当其输出特性值增大时，其 性能逐渐变好，质量损失逐渐变小，其理想的值是无穷大。这种质量特征称为望大特性，如 粘接强度等。很明显，望小特性的倒数与望大特性具有相同的性质，所以可以用y来代替 式子(3)中丁，即得望大特性得质量损失函数为：L(y) = K (1/y2)(5)如图5所示，其质量损失系数K可以由功能界限和相应的损失，或者容差界限和相应的 损失来确定。L(y)为随机变量时，则产品的质量水平用期望值反映，对于N件产品，若望大 特性的测试值为

14、y1，y2， ， yN，则平均质量损失为：图5望大特性的质量损失函数4 关于信噪比与正交表在通讯和电气工程中，为了对所选择设备的质量特征进行量化引入“信噪比”(输入信号强度与噪声强度之比) 这个概念。田口将这个概念引入到正交试验设计中，用它来模拟噪 声因素对质量特性的影响。1. 望目特性的信噪比设产品的望目特性值为y0,质量特性y服从正态分布yN(p y，o y),，信噪比计算公式为：SN = 101g丄(S -S 2)N m y6式中：b 2 二 S 2 仝(y - y)2 /(N -1)y yii=1;s = n( y)2m；y 为质量特性值，y为质量特性均值。i2. 望小特性的信噪比信噪

15、比计算公式为：SN = -10lg(l y2)Nii=13. 望大特性的信噪比信噪比计算公式为：SN = -10lg(l 丄)N =1 y 27)(8)(9)正交表是一些已经制作好的规格化的表，是正交试验设计的基本工具。正交表的每一列 等同于一个因素。每一列中，各水平重复出现的次数是相等的，并且任意两列中，各水平在 相同横向上的搭配也是均衡的。这些特定保证正交表安排的试验，具有均衡分散性、整齐可 比性。举例来说，对于3因素3水平的试验，若全部做全需要33=27次，而用正交表进行 的试验值需要9次，这9次在全体27次试验中是均衡分散的，具有很强的代表性。田口稳健设计中的参数设计一般都需要用到两个

16、正交表，即用于安排可控因素的正交 表，称为“内表”或“设计变量矩阵”；用于安排噪声因素的正交表称为“外表”或“不可 控因素矩阵”示例如表1所示。表中，A、B、C、D、E和F是可控因子，e是误差因子， 所谓误差因子，是田口方法提供的对未考虑到的系统可控因子的一个补充。若该误差因子对 系统性能影响较大，则说明还有为考虑到的其它的重要可控因子。U、V和W是噪音因子。表 1 因子正交试验设计的基本结构正父表类型内表L8（27）外表L4(23)噪声因子安排和行数试验因素可控因子安排及行数试验次序噪声因信1234素安排噪、列号12345671122U比X1212VSN试验次序、ABCDEFe1221W11

17、111111y11y12y13y14SN21112222y21y22y23y24SN?31221122y31y32y33y14SN341222211SN452121212sn562122121SN672211221sn782212112y81y82y83y84SN85 田口参数设计的流程参数设计是田口稳健设计的重要内容，它的工具是前述的正交试验设计和信噪比。参 数设计的工作步骤可按图6 的框图进行。图 6 田口参数设计流程图6 田口参数设计案例本案例选择了北美一家汽车零部件供应商生产的产品为，该零部件供应商主要向整车厂 供应汽车面板，如图7 所示，图中，实线部分表示铸型腔尺寸，虚线表示零件尺寸

18、，该汽车 面板采用模具注塑法工艺生产，待注塑并冷却后，存在大量零件收缩的质量问题，所谓零件 收缩，是指零件达到周围环境温度后的尺寸与铸模型腔尺寸之间的差异。该质量特性参数受 许多工艺参数和环境变量的影响。该注塑装置如图8 所示。对收缩率的控制目标或规格(%)：0.4+/- 0.1。对当前状态的评估结果为(样本容量 N=49)：Mean：0.45%；StDev：0.07%；Pp = 0.5Ppk = 0.24DPM = 250KCp=0.45Cpk=0.21由此可见，该注塑工序生产能力严重不足。图 7 汽车面板轮廓(1) 建立p图。通过分析，该注塑工序受下列因素的影响：不同的供应商、冷却时间、

19、融化温度、螺杆转速、填充时间、填充压力、喷嘴直径、模子壁的温度、小球尺寸的内在变 异、小球的研磨度。可以将上述工艺参数进行分类，根据田口参数设计基本思想，可以得到 如下的分类：输入变量( W ) ：小球供应商( Supplier ) ( A ) ；信号因子(X)：冷却时间(Cool time)(B)融化温度( Melt Temp.)( C)； 螺杆转速( Screwspeed)( D)； 填充时间(Fill Time)(E)； 填充压力( Fill Press)( F)； 喷嘴直径( Nozzle Dia.)( G)； 模子壁的温度(H)；噪音因子(Z)：小球的再研磨度(Z1)； 小球尺寸的内

20、在变异(Z2)；结果如图9 所示。目前的工艺参数设置为：A：小球供应商(II)；B：冷却时间：29sec；C：融化温度:250F；D：螺杆转速：150rpm；E：填充时间：4sec；F：填充压力:50Mpa；G：喷嘴直径：6.5cm；H：模子壁的温度(MoldWTemp.)： 75OF；久小球供应商+漏斗 G:喷嘴直径(66.5, 7 cm) (If n) Hopptar图8注塑装置喷嘴NonleMatcirialHeaters射腔加热器铸型腔尺寸-赛欧零件尺寸B:冷却时间(30, 29, 28 sec)C：融化温度(240, 250, 260旧E:填充时间(3,4, 5 sec)F：填充压力

21、(45, 50, 55畑)H：模子壁的温度( 70, 65旧信号冋控)变量(Xrs)B:冷却时间(30, 29, 28 sec)C：融化温度(246 250, 260旧D:螺杆转速(125,150,175 rpm)E:填充时间(3, 4, 5 sec)F：填充压力(45, 50, 55册)G:喷嘴直径(6县7 cm)H:模子壁时温度(75, 70, 65販)水目前设置：粗体 且下划线输入变量(Wfs)A:小球供应商(1,11)输出(Yrs) *零件收缩率()不可控变量(Zfs)Z1:小球研磨度(5%, 10%)Z2：小球尺寸的忻在交异(小-中)图9构建p图(2) 创建内外表。为进行信噪比的计算

22、分析，首先需要创建内表和外表。内表是输入变 量(W)和过程参数(X)的实验设计(如，部分设计因子设计)。外表则是噪音变量(Z) 的实验设计，注塑成型内表：DOE因子和水平如下：A：小球供应商(I, II)B:冷却时间(30, 29, 28 sec)C：融化温度(240, 250, 260 F)D:螺杆转速(125, 150, 175 rpm)E:填充时间(3, 4, 5 sec)F:填充压力(45, 50, 55 MPa)G:喷嘴直径(6, 65, 7 cm)H:模子壁的温度(75, 70, 65 F)由此可见，共8个因子，其中第1个因子2个水平；后7个因子3个水平，则全部可 能实验组合为：

23、37 * 21 = 4,374，如此之多实验一为实验资源所不容许，二是实验时间不允 许，所以有必要借助实验设计进行分析。如果用Mi ni tab软件进行分析，则该实验是8个因 子的混合水平设计，分别为2层和3层。步骤为：(1)选择“Taguchi Design”出现如图 10 (a)所示对话框，选择“Mixed Level Design” 和 “Number of factors： 8”； (2)选择“Display Available Desig ns”，确定可能的设计方案，如图 10 (b)所示，“ Mixed 2-3 Level Design”中的L18符合本例中的已知数据；(3)在创建

24、设计后，点击“factors”进一步确定变量、层次和交互作用，如图10 (c)所示；(4)确定因子及水平，Minitab会给出所有可能的交互作用，由分析人员进行选择如图10 (d)所示。(a)(b)(c)(d)图10 Mi ni tab创建内表这样，Minitab将创建一个实验(非随机化)，该实验设计表显示了田口 18次实验组合, 如表2所示。表 2 内表Run SupplierBCcclTimeCMeltTempDScrewSpeedEFillTimeFFillPressGNczzleDiaHMcldW/Temp1111111112112222223113333334121122335122

25、23311612331122713121323813232131913313212102113322111212113321221322113132212313214222312131522312321162313231217232131231823321231然后建立外表，由前面可知，实验噪音变量有小球再研磨度（低水平 5%和高水平 10%）， 小球内在的尺寸变异（低水平为小尺寸，高水平为大尺寸），2 个噪音变量分别设置为2 个 水平，则有4 种组合，如表3 所示。表 3 外表再研磨度（ Regin d）低水平（L）咼咼水平（H）小球尺寸（Pellet size）低水平（L）1112高水平（

26、H）2122综合表（2）和表（3）可得到综合的内外表，如表4 所示。表 4 综合的内外表RunASupplieBQccITim/MeltTemipD：ScrewSpe、E FillTimpRFillPresrB：CoolTimfCMeltTempDScrawSpa、E FillTim=FFillPresSG：NozzleDaH MoldWTenS11S12S21S221111111110 530 580 60 612112222220 430 430 530 513113333330 270 280 350 394121122330 250 290 310 375122233110 430 4

27、50 460 496123311220 280 30 380 427131213230 140 19UbwO0 250 288132321310 20 240 270 349133132120 340 350 360 3810211332210 440 490 560 6411212113320 450 480 60 6412213221130 40 460 5505813221231320 410 50 57v w0 6514222312130 180 260 290 3515223123210 360 370 490 5616231323120 150 230 280 3417232131

28、230 30 370 360 43-18.23321231027035049053（3）计算信噪比。该信噪比为望目特性，信噪比计算公式如式（7）所示。表6 是信噪比 的计算结果。通过Minitab,可以得到其均值响应表和S/N响应表，分别如表7和表8所示，主效应图如 图 11 所示。表 6 信噪比计算结果RunS11S12S21S22AvgSS/N10.530.580.60.610.580.03562.4220.430.430.530.510.4750.05261.9130.270.280.350.390.32250.05741.5040.250.290.310.370.3050.051.57

29、50.430.450.460.490.45750.0252.5260.280.30.380.420.3450.06611.4470.140.190.250.280.2150.06241.0780.20.240.270.340.26250.05911.3090.340.350.360.380.35750.01712.64100.440.490.560.640.53250.08691.57110.450.480.60.640.54250.09181.54120.40.460.550.580.49750.08261.56130.410.50.570.650.53250.10211.43140.180

30、.260.290.350.270.07071.16150.360.370.490.560.4450.09681.33160.150.230.280.340.250.08040.99170.30.370.360.430.3650.05321.67180.270.350.490.530.410.12111.06表 7 均值响应表LevelA:BCDEFGH10.3688890.4916670.4025000.4325000.3937500.4304170.4020830.44791720.4272220.3925000.3954170.4312500.3725000.3916670.396250.

31、41708330.3100000.3962500.3304170.4279170.3720830.3958330.329167Delta0.0583330.1816670.0070830.1020830.0554170.0583330.3958330.329167Rank4.51.07.03.06.04.58.02.0表8 S/N响应表LevelA:BCDEFGH118.196717.520315.104418.628614.499516.368718.832017.0048213.685715.757516.851415.939014.411916.534914.987316.5860314

32、.545715.867813.255918.912114.919914.004314.2327Delta4.51102.97461.74705.37274.50021.61514.82772.7721Rank35714826Main Effects Plot for MeansA:Supplier C:MeltTemp E:Fiirime G:NozzleDiaB:CoorTime D: Screws peed F:FillPress H:MoldWTenA B C D E F G H图 11 均值响应的主效应图(4)根据S/N值和均值响应确定因子设置。遵循以下三个步骤：第1 步：寻找那些因子

33、对应的均值delta排序较大，S/N排序较小； 选择最优设置以满足目标：均值越接近于目标值越好；第2步：在余下的因子中，寻找较高S/N排序的因子； 确定使 S/N 最大化；第3步：对于那些具有较低均值delta排序和较低S/N排序的因子，所作选择也许比较容易 实现健壮；由第1步可以得到推荐的设置结果为：B-2, H-2, F-2；由第2步可以得到推荐的设置 结果为：D-1, F-1, A-1, E-3, G-1；由第3步可以得到推荐的设置结果为：C-2。用Minitab 可以在此基础上进行预测，可以得到如下结果：S/N=28.711, Mean=0.441667,此结果的均值响应偏离较大，考虑到 E 因素所选水平的均值偏离较大，故将 E 由水平3 调整为水平1， 得到新的结果为：S/N=24.984, Mean=0.4075, Stdev=0.221924,新的设置改善较为明显, 通过这样的分析，就可以最终将可控因子的水平组合确定下来，即：A-1, B-2, C-2, D-1, E-1, F-2, G-1, H-2