DOE案例(minitab实验设计).docVIP

DOE试验案列 （提高金属热处理后的强度） 按计划表完成试验并将试验结果填入表中。 利用Minitab软件，对结果做因子主效图、交互效应图和立方图。如图1-1(a)、(b)、(c)所示： 从(a)图可以看出：A(热处理温度)、B(升温时间)及D(恒温时间时间)主效应显著。 从(b)图可以看出B(升温时间)跟D(恒温时间)存在明显交互作用。 从(c)可以看出强度的最大值为574.5MPa，它对应的各因素水平分别是：热处理温度860℃、处理时间1.6h、升温时间3min、恒温时间60min；即：当选择热热处理温度860℃、处理时间1.6h、升温时间3min、恒温时间60min，可获得较好的强度结果。 作标准化效应的Pareto图和正态图，如图1-2(a)、(b)所示。 从上两图可以看出A、B、D显著，C不显著，BD交互作用处于临界点，做显著处理。 作残差图，如图1-3所示。 从上图可以看出：残差满足正态分布和随机波动的要求。无异常现象。 增加B*D项，对实验结果最方差和回归分析。如图1-4所示。 图1-4 从上图可以看出C(处理时间)不显著,需重新修订。 去掉C项，作再次回归分析。如图1-5所示。 图1-5 对修订回归方程再做残差诊断，残插图如图1-6所示。 从下图可以看出残差服从正态分布，无异状。 图1-6 修订后残插图 P值=0.9350.05，残插符合正态分布。 调优找出因子最佳方案。 当热处理温度=860℃,升温时间=3min，处理时间=1.6h，恒温时间=60min,强度最大值 Y=573.6.8 8 这个值最大 回归分析:强度MPa 与 A(热处理温度）℃, B(升温时间）/min, C（处理时间）/h, D（恒温时间）/min, B*D 回归方程为 强度MPa = 192 + 0.496 A(热处理温度）℃ - 63.6 B(升温时间）/min + 20.0 C（处理时间）/h - 2.52 D（恒温时间）/min + 1.46 B*D 自变量 系数 系数标准误 T P 常量 191.55 90.86 2.11 0.055 A(热处理温度）℃ 0.49625 0.06324 7.85 0.000 B(升温时间）/min -63.60 27.94 -2.28 0.040 C（处理时间）/h 20.00 12.65 1.58 0.138 D（恒温时间）/min -2.520 1.290 -1.95 0.073 B*D 1.4600 0.5059 2.89 0.013 S = 5.05902 R-Sq = 91.3% R-Sq（调整） = 87.9% 方差分析 来源 自由度 SS MS F P 回归 5 3479.57 695.91 27.19 0.000 残差误差 13 332.72 25.59 合计 18 3812.29 来源 自由度 Seq SS A(热处理温度）℃ 1 1576.09 B(升温时间）/min 1 1115.56 C（处理时间）/h 1 64.00 D（恒温时间）/min 1 510.76 B*D 1 213.16 p0.05不显著 回归分析:强度MPa 与 A(热处理温度）℃, B(升温时间）/min, D（恒温时间）/min, B*D 回归方程为 强度MPa = 222 + 0.496 A(热处理温度）℃ - 63.6 B(升温时间）/min - 2.52 D（恒温时间）/min + 1.46 B*D 自变量 系数 系数标准误 T P 常量 221.55 93.50 2.37 0.033 A(热处理温度）℃ 0.49625 0.06654 7.46 0.000 B(升温时间）/m