实验设计DOE文档格式.docx

1、实验设计在质量控制的整个过程中扮演了非常重要的角色，它是我们产品质量提高，工艺流程改善的重要保证。实验设计已广泛运用了从航天业到一般生产制造业的产品质量改善、工艺流程优化甚至已运用到医学界。籍此课程，您将通过对产品质量，工艺参数的量化分析，寻找关键因素，控制与其相关的因素。根据实际需求，学习判别与选择不同的实验设计种类，设计你的实验步骤，发现如何控制各种影响因素，以最少的投入，换取最大的收益，从而使产品质量得以提升，工艺流程最优化。 在三十、四十年代，英、美、苏等国对实验设计法进行了进一步研究，并将其逐步推广到工业生产领域中，在冶金、建筑、纺织、机械、医药等行业都有所应用。二战期间，英美等国在

2、工业试验中采用实验设计法取得了显著效果，战后，日本将其作为管理技术之一从英美引进，对其的经济复苏起了促进作用。今天，实验设计已成为日本企业界人士、工程技术人员、研究人员和管理人员必备的一种通用技术。五十年代，田口玄一博士借鉴实验设计法提出了信噪比实验设计，并逐步发展为以质量损失函数、三次设计为基本思想的田口方法。八十年代，田口方法进入美国，得到了普遍关注。如今，实验设计技术的应用领域已经突破了传统的工业过程改进和产品设计范畴，广泛地渗透到商业布局、商品陈列、广告设计及产品包装的应用之中。我国在六十年代就曾对实验设计进行了研究和推广，八十年代又引入了田口方法，取得了一定成效。但实验设计作为一种质

3、量改进的有力武器，还尚未发挥它的全部威力。二、实验设计（DOE）方法实验设计（DOE）方法一直在发展、丰富和完善，美国人多里安谢宁（夏宁）（Dorian Shainin）总结了如下七种新的DOE应用技术：1、 多变图技术 （X-R图分析） 根据以往经验确定影响质量的可能要素，例如工人班次、机床、原材料、工艺变量等，每隔一段时间抽取连续生产的几件产品，按需观察的这些要素的几种情况分别测试质量特性，画成图表进行比较分析，以确定引起波动的原因。2、零件搜索技术 （零件交换实验） 根据以往经验确定影响产品质量的可疑零件，随机选取几个好的产品和坏的产品，将坏产品上的可疑零件与好产品上的对应零件进行交换，

4、重新装配后进行质量特性参数的测量、比较、分析，以找到影响产品质量的主要零件。3、成对比较技术 随机选取5对以上的好的和坏的产品，用各种方法测试其各种参数并比较其不同之处，以确定影响产品质量的主要原因。成对比较技术适用于不可拆卸的产品。4、变量搜索技术 与零件搜索有许多相似之处，都是逐个替换，进行测试比较，以搜索引起产品质量波动的主要原因，但变量搜索是针对变量的，而零件搜索则是针对零件的。5、完全析因技术 采用以上四项技术，寻找出4个以下的主要因素，这些因素按全部因素所有水平的一切组合逐次进行实验，研究这些因素的主效应和相互之间的交互作用，以确定最佳的因素水平匹配方案。6、 改进效果检查 令B为

5、改进后的工序，而C为改进前的工序。为了验证改进效果，可以设定风险率a（第类风险）和b（第类风险），随机抽取B和C两种产品（样本量由风险率决定），进行假设检验，以确定在规定的风险率下，B产品是否优于C产品。7、 实验设计的回归分析 对散布数据作出散布图，应用回归理论，诊断出对产品影响大和小的因素，从而找出影响产品质量波动的主要原因，并根据波动大小，确定各因素的目标值和容差。三、案例分析力克斯是某机械工业公司年轻的助理工程师，他接了老板交给的改进项目：用6西格玛方法（使用DOE手段）改善产品流水线合格率。 老板告诉他：“我们需要你改变生产流水线目前的混乱情况，自从我们更换了新的硬化剂供应商后，合格

6、率就直线下降了，你的项目是使合格率重新回到原来92的历史水准。” 这个年轻助理工程师想，“我一定不能辜负公司的信任，可是要完成这个任务，不是件容易的事，听说已经有好几个前辈反复作试验，都没有成功呢”。 他想：我应该从一个SIPOC和绘制过程图开始，然后去找生产线的工程师菲利斯了解情况，寻找解决课题的方法。 菲利斯介绍说。“我们用一种硬化剂和其他成份混合组成一种粘合剂粘结工件。然后我们将批量的组装件装在盘里放进烤箱加热至155度100分钟，过去合格率始终在90以上。大约三个月前采购部指定了一个新的硬化剂供应商，使用后合格率立即下降到现在的水平，一天最好时也只有77。“所以你认为新的硬化剂有问题”

7、亚力克斯问道。“难道不是很明显的吗”菲利斯反问，“我们以前从来没有发生过这样的问题。“但是对于新的硬化剂可能有一套更好的烤箱温度和时间设置，”亚力克斯直率地说道，当他看到菲利斯睁大了眼睛，他后悔自己说话太冒失了。 “听着大学生男孩，当你还是高中生时我就在这个部门工作了，并且我们一直用155度100分钟。我们也试尝过用其它组合，但没有一个比155度100分钟更好。每次有人调整设置，使它过高或过低操作，合格率都会下降。除此以外，如果由于自动温度控制计失灵造成温度低于标准时也会发生同样的事情。目前我们的烤箱温控计没有问题，当合格率下降时，我做的第一件事就是用计量器检查烤箱温度，发现烤箱和温控计都正常

8、。而且我还要告诉你，对于新的硬化剂的最佳使用设置条件，我也验证过。我们上个月做过一些试验，发现目前的设置是最好的，合格率却只有75。采购部变更供应商，作出这么重要的决定之前应该事先与我们商量 ” “看来我需要联系一下供应商。”亚力克斯回答到，他暗自庆幸菲利斯将怒火发向采购部而并非他自己。“我可以拿到一份你手头数据的复印件吗”亚力克斯问。“小伙子，如果这能对产品合格率有帮助，当然没问题”菲利斯说。两天后，亚力克斯收到了菲利斯给的文件。包括试验数据记录和说明。亚力克斯： 这是你需要的新硬化剂的试验数据。首先我们在通常情况下生产了一批，合格率为76。然后我们降低了10度，在145度再生产一批，合格率

9、69。数据上来看，提高温度对合格率有好处，所以我们将烤箱温度调度至175度。没有获得预想的效果，我们只得到70合格率。由此判断，最佳温度仍然是155度一点也不意外，我已经说过，根据以往经验，这是最佳运行温度。 在我们证明了155度是最佳合格率点后两天，我们改变了烤制时间，烤制时间为60分钟时的一批合格率是52，烤制时间为120分钟时的一批合格率是73。 当你有了一个新硬化剂供应商来解决这问题时 或恢复用回原来的供应商 请告知我。“我得先看一下数据的状态”。亚力克斯一边想一边打开电脑并输入数据。首先亚力克斯使用三个数据观察合格率与温度之间的关系（见图1）：看来数据支持菲利斯的结论。目前设置的时间

10、和温度对应的最佳合格率点。如果改变时间或温度都会导致合格率下降。菲利斯还是对的。亚力克斯合上笔记本。明天我将和我的6西格玛指导讨论并听听他的意见。两天后，亚力克斯回到了车间。“小伙子，你不相信我的数据吗”菲利斯问道。 “你的数据是对的，菲利斯，我只是想再做一个你们通常生产情况下的DOE。你已经做了OFAT试验，一次改变一个因子。但是有可能存在时间与温度的交互作用，所以我们必须同时改变两个输入以便发现问题。我的指导告诉我你的数据可能错过了合格率更高点。 “好的，对我来说无所谓，按照你的想法去做吧。我下星期准备休假了。”菲利斯回答。他想，让他去烦助理员约翰吧，早晚亚力克斯会明白是这些垃圾硬化剂在作

11、怪，他一定会在我休假回来之前走掉的。菲利斯度完假后回来。他端着咖啡来到助理员身边坐下。“那么，约翰，合格率怎么样了那孩子有没有用新的供应商来解决问题”菲利斯问。“是这样”，助理员回答，“我们为亚力克斯做了四批产品并且最少有一批合格率在80以上。有趣的是，恰恰是使用了新供应商提供的硬化剂。亚力克斯说他午饭后会来和你一起回顾试验结果。他在办公室留下了一张图表给你，真高兴我们最终认识到了问题所在” 约翰可以从菲利斯脸上看出他根本不相信。当菲利斯很快跳起来并冲向办公室去检查数据时，约翰笑了起来（见图5）。办公桌上有一封信。菲利斯： 看来可以通过提高温度并减少时间来提高合格率。当你回到办公室看到如果我们

12、能提高合格率到82时，我正在进行另外一次试验。我们试图通过减少炉中烤制时间来达到90以上的合格率。缩短时间将会大大提高你的生产效率。 祝你假期愉快 亚力克斯 “大家可能还记得我上次的发表，将两个因子放在一起做了合格率的DOE，结果我们发现用提高温度和缩短时间可以达到更高合格率。第一次DOE使合格率升高到82，而我们的目标是达到90合格率。我们发现时间和温度间有很强的交互作用，并决定进行另一次推论空间的DOE。 我们还与供应商一起讨论过，他们确认我们需要提高温度来提高合格率。在第一次DOE结论的基础上，我们又改变了推论空间，进行了第二次DOE。下面是第一和第二DOE的结果图。 “看来我们发现了一

13、个能获得95合格率的操作点，这个水平已经超过了使用老供应商硬化剂的合格率。工件只需在烤箱内烤制30分钟而不是原来的100分钟，因此节约了70的循环时间。工程师最初认为新供应商的硬化剂是合格率低的根本原因，我也猜想是对的。但是他漏掉了新硬化剂带给我们提高合格率缩短周期的机会，我们只需要探索设计空间并让数据说话。“接下来呢”指导者问。“我们进行了许多次验证，结果都是好的，他们将使用新设置点生产”亚力克斯肯定地说道。“听来很好，干得不错，亚力克斯”指导者夸奖道。总结 从以前的等位线图（见图6）可以看出，菲利斯的试验区域选在温度155度100分钟的范围内，而且只进行了OFAT试验，未能发现温度与时间的

14、交互作用，错过了最佳的设置区域。通过全因子试验和最优化设计，终于达到甚至超过了合格率的历史水平。Oliven:很好的案例，如果能将亚力克斯的DOE做法显示出来就完美了。我大概地理了一下思路：1.亚力克斯面对的问题已经有人作过了OFAT实验，其结果显示，单一得改变温度或时间的设置都会降低合格率。2.亚力克斯作了第一次的DOE全阶乘设计：温度HL=165,LL=145;时间HL=130,LL=70。实验结果显示，温度在165，时间在70时合格率达到了82%.3.亚力克斯并未停止前进，他又作了优化实验，也就是一般意义上更为复杂和精确的实验，如响应曲面法。最后根据等高线图，寻找到了合格率95%的鞍点。

15、4.最后他们对最佳的设定作了验证实验，从案例中看，似乎实验选取的最终设定为温度180，时间30分钟。但是我不太明白有几点：1. 亚力克斯似乎一直在拓宽他的因子水平的范围。温度：第一次145-165,第二次145-185,时间：第一次70-130,第二次30-150.这几乎是违背DOE实验原则的，我们分步实验会一步步缩小因子的有效设定范围，拓宽是不是有风险呢？当然也有可能因为优化实验是RSM，所以不怕非线性的影响。2.这个案例似乎想告诉我们实验因子的设定中，交互作用的影响是OFAT的一个软肋，以及实验因子水平设定如果太窄的话，可能会漏掉最佳的结果。 并不能把亚力克斯的实验作为一次正规的DOE实验

16、步骤，因为如果是一个黑带老手的话，OFAT就能看出交互作用来，DOE全阶乘实验就可以找到最佳结果。3.既然是交互作用影响显著，为何亚力克斯一口咬定提高温度缩短时间就可以提高合格率呢？他为何不降低温度延长时间呢？当然也许是工厂都需要快吧！四、实施方法 Experiment, DOE）是以概率论与数理统计为理论基础、合理安排试验的一种方法论,它研究如何高效而经济地获取数据信息、科学地分析处理、得出正确的结论。本文从案例入手,介绍DOE在各行各业的应用实例,然后介绍单因素试验设计方法和多因素试验设计的思路,进而讲解正交设计为主、简介回归正交设计、均匀设计和三次设计。本讲义假设学员已经掌握统计学关于数

17、据处理的基础知识,但为便于学员统一基本概念,在开头穿插了关于数据变异、回归分析、方差分析的内容复习。1 引言试验设计在管理体系中的地位及其应用实例 过程性能的测量应兼顾各相关方的需求和期望,可包括:能力;反应时间;生产周期或生产能力;可信性的可测量因素;产量;组织的人员的有效性和效率;技术的应用;废物的减少;费用的分配和降低”。这里写出的9方面如何进行优化,无疑是组织经常面临的具体课题,设计开发（design and development）就是解决这些课题的重要环节。统计技术是解决上述课题的一种有效手段。这里的关键词是“安排试验”。它是专业技术和统计技术的结合。请看一些多因素DOE应用实例。

18、 【例1农作物栽培】某省供销社进行大棚丝瓜早熟栽培试验,选择4个播种期、2个品种、2种株行距进行条件组合试验,实现了提早丝瓜上市时间、提高产值的目的。 【例2材料替代】某市对豆腐制作中起凝固作用的GDL、盐卤、石膏3种物料配比进行条件组合试验,寻找出可代替原用的单一凝固剂而不影响产品质量、但避免了产生废渣的复合凝固剂,推出“环保型豆腐”新工艺。 【例3环保工程】某科研所对pH值、絮凝剂（聚合硅硫酸铝）投加量、其中的硅铝比等3种因素各4个水平进行条件组合试验,寻找处理印染废水时提高“色度去除率”的优化工艺条件。 【例4建筑工程】三峡大坝泄洪深孔抗冲耐磨混凝土性能对大坝安全至关重要,主要受水胶比、

19、用水量、粉煤灰品种、粉煤灰含量、砂率、引气剂、减水剂等影响。施工指挥部对该7种因素各两种选择进行条件组合试验,确定了最优的参数配合。 【例5质量检验】某市商检局研究在茶厂出厂检验、外贸收购验收及出口检验原用0324h恒重法和13027快速法检测茶叶中水分含量,为了改用更好的红外快速检测法,对温度、样品量、检测时间等3个影响检测误差的主要因素进行各3水平的条件组合试验,确定了利用红外水分分析仪以130称量约5g样品检测15min的最佳实验条件,在确保一致性的前提下相对于原法具有快速、简便、均匀、耗能少、占地小、结果准确等优点。 【例6节能降耗】指导实施的优化脱盐制水工艺、降低酸碱消耗的QCC活动

20、,通过两轮正交试验取得成功,相关成果获当年全国优秀QCC并经过全国“擂台赛”选拔作为我国代表出席国际QCC代表大会发表。 【例7焊接工艺1】现代通信电子电路的制造大量使用表面贴装元件以提高组装密度，然而许多电子产品应用中（如连接器）需要通孔插装式元件以获得高强度互连。近年来发展的通孔再流焊工艺,是用表贴元件的制造方法完成混装电路板的焊接,这种工艺取消了波峰焊工序,从而降低了生产成本。某省采用正交试验的方法分析多个因素对焊点强度及可靠性的影响规律,确定了通孔再流焊工艺优化的生产方案。 【例8焊接工艺2】波峰焊中无铅（Pb-free）焊接技术上必须解决一些现实问题,包括无挥发性有机化合物（VOC-

21、free）的助焊剂技术,以及修改工艺来接纳所要求的更高焊接温度。某作者经过正交试验确定了条件组合: 焊锡温度265270、接触时间1.8秒、预热时间（顶面）110、湿助焊剂量474 mg/dm2,取得成功。 【例9化工工艺】某化工产品期望指标堆密度（克/立方厘米）以0.50.6为佳。某厂为了寻求的最佳反应条件,从专业知识考虑选择了反应物、浓度1、2（摩尔/升）,反应温度3（）,加料方式4（取1,2）等4个可控因素,分别设定了取值范围和水平数值。进行试验取得产物、处理并测定其试验结果分析:用多元逐步回归方法建立数学模型:经逐步筛选,得回归方程为:= -0.004351 + 0.335441+0.

22、772154 （与2、3无显著相关关系）,且检验表明方程的可信度高;求回归方程的最优解:当四个依次为1.3、1.0、27、1时方程有最小值0.570;计算它们的全面组合共12种,得“最佳条件全排列”,从中得出结论:只要保持4取1方式 ,1在1.01.5范围内均能保证在0.50.6的所需范围内。 因此在进行DOE设计时，应确定因子及每个因子的取值范围。2 几个统计技术基本概念复习2.1数据变异的统计量一批参差不齐的数据,动态形成波动,静态形成分布。不同来源、不同性质的数据具有遵循不同的分布规律, 要考察的包括三个方面:数据的集中位置即平均水平、离散程度和分布形态。 平均水平用算术平均值代表,离散

23、程度用标准偏差代表,分布形态用某种分布规律的“适度”代表。设n个数据1、2、3in,其平均值和标准偏差分别用下面介绍的贝塞尔（Bessel）公式计算:平均值 变差”或偏差平方和 总体方差 总体标准（偏）差样本方差样本标准（偏）差设m组共n个数据为:1有1个,2有2个,j有j个m有m个（这里1,2,m称为频数）,则:总频数平均值“变差”或偏差平方和 总体方差、总体标准（偏）差、样本方差、样本标准（偏）差计算方法同上。另有一个显示离散程度的统计量“极差”R,等于各X中的最大值与最小值之差（取绝对值）。 【例题】n=5个数据:1,4,8,10,12;求以上各统计量: 总和=35,=7;平方和=325

24、,变差（325-35*35/5）/5=80;=80/5=16,= =4;=80/4=20,=20=4.47, 又R=12-1=11。本来,变差的定义是各（）的平方和,用本例数据不难验证。2.2 方差分析 （ANOVA, analysis of variance） 方差分析用于对误差因素（包括方法误差、模型误差、仪器误差、操作误差等等）而言各种因素对总变异的影响进行判断,回答这样一个问题:试验结果的变异究竟是误差所致,还是因素确实产生影响?【例1 单因素重复试验的方差分析】考察反应温度A（）对收率X（）的影响。对A因素选p=2个条件:A1=30、A2=40,各进行c=5次试验,结果如下表。粗看起

25、来,平均值A2A1,但各次试验单独看结论不同,例如就完全相反。要求判断结果的不同确实是因素A的不同水平起作用,还是仅仅由于试验误差所致。按下左表计算得SA与SE,然后作方差分析表如下右表,结论:因素A的水平差异所造成的影响并不显著。源数据汇总处理方 差 分 析 表因素水平试验结果 （重复c=5次）变差c2平均值平均值的变差SA=184.9方差来源偏差平方和S自由度fF比F检验临界值显著性=0.050.01A17578606183429.271.4因素Ap-1=11.335.3211.26A28962937185680.080.0误差E1109.2p*（c-1）=8全部pc个数据变差ST=129

26、4.1SE=1109.2总计T1294.1pc-1=9此行供参考,可不列出表中“F比”是要比较的因素的“均方和”对误差因素的“均方和”的比值,等于（SA/fA）（SE/fE）。从表中可以看出,操作过程实际上是把全部10个数据的变差ST分解为A和E两个因素所分别引起的两部分。本例中误差因素影响太大,表面上A的两个水平的平均值不同,实际上差异不算显著。如果在技术上觉得A因素应该产生显著影响,应进一步考虑查找和减少误差因素,以期得到符合实际的结论。【例2 双因素无重复试验的方差分析,p=5,q=4,c=1】（引自化工企业全面质量管理乙编上册用相同的原料按A1A5五种工艺本别用B1B4四个反应器合成某种化工产品各一锅,分别测得成品种某副产物含量Y（%）如下表中20个数据,问工艺A、反应器B两因素是否对该副产物含量有显著影响? 对表中横向计算变差与总和、纵向计算总和,再按i求SA、按j求SB,另求SE如表;然后作方差分析表如下右表。各因素自由度f计算为:fA= p-1=3,fB=q-1=2,fE=（p-1）（q-1）=6,fT= pq-