试验设计课程案例分析资料.docx
《试验设计课程案例分析资料.docx》由会员分享,可在线阅读,更多相关《试验设计课程案例分析资料.docx(19页珍藏版)》请在冰点文库上搜索。
试验设计课程案例分析资料
实验设计课程案例分析
小组成员:
45120202白合兵
45120203田大越
45120212刘永勉
2.1正交实验设计概述
正交实验设计(Orthogonalexperimentaldesign)11是研究多因素多水平的一种设计方法,它是根据从全面实验中挑选出部分有代表性的点进行实验,正交实验设计又称正交设计或多因素优选设计,是一种合理安排、科学分析各实验因素的一种有效的数理统计方法。
它是在实践经验和理论认识的基础上,借助一种规格化的“正交表”,从众多的实验条件中确定出若干个代表性较强的实验条件,科学地安排实验,然后对实验结果进行综合比较,统计分析,探求各因素水平的最佳组合,从而得到最优或较优实验方案的一种实验设计方法。
正交实验设计的特点是用不太多的实验次数,找出实验因素的最佳水平组合,了解实验因素的重要性程度及交互作用情况,减少实验盲目性,避免资金浪费等。
它能以较少的实验次数找到较好的实验(生产)方案,由正交实验寻找出的优化参数(条件)与全面实验所找出的最优条件有一致的趋势。
正交实验设计具有正交性,使实验具备均衡分散和综合可比性。
此法应用方便,准确性高,在多因素条件下应用有很大的优越性,是一种高效率、快速、经济的实验设计方法。
日本著名的统计学家田口玄一将正交实验选择的水平组合列成表格,称为正交表。
例如作一个三因素三水平的实验,按全面实验要求,须进行33=27种组合的实验,且尚未考虑每一组合的重复数。
若按L9(3)3正交表按排实验,只需作9次,显然大大减少了工作量。
因而正交实验设计在很多领域的研究中已经得到广泛应用。
2.2正交实验设计基本程序
正交设计中常用的术语有:
指标、因子和水平。
正交设计把实验设计要考虑的结果和评价准则称为指标,一般以yi表示第i次实验的指标值;把对实验结果和对评价指标可能产生影响且在实验中明确了条件加以对比的因素称为因子,一般以大写字母表示;把每个因子在实验中的具体条件称为因子的水平,简称水平,一般以表示因子的大写字母加上脚标来表示。
对于多因素实验,正交设计是简单常用的一种设计方法,其设计程序12如图4所示。
Figure4Theprogrammingofpotassiumdiformate
2.3正交实验方案设计
2.3.1明确实验目的,确定实验指标
实验设计前必须明确实验目的,即本次实验要解决什么问题。
实验目的确定后,对实验结果如何衡量,即需要确定出实验指标。
实验指标可为定量指标,如强度、硬度、产量、出品率、成本等;也可为定性指标如颜色、口感、光泽等。
一般为了便于实验结果的分析,定性指标可按相关的标准打分或模糊数学处理进行数量化,将定性指标定量化。
2.3.2选择实验因素,确定实验水平,列出因素水平表
根据专业知识、以往的研究结论和经验,从影响实验指标的诸多因素中,通过因果分析筛选出需要考察的实验因素。
一般确定实验因素时,应以对实验指标影响大的因素、尚未考察过的因素、尚未完全掌握其规律的因素为先。
实验因素选定后,根据所掌握的信息资料和相关知识,确定每个因素的水平,一般以2-4个水平为宜。
对主要考察的实验因素,可以多取水平,但不宜过多(≤6),否则实验次数骤增。
因素的水平间距,应根据专业知识和已有的资料,尽可能把水平值取在理想区域。
2.3.3选择合适的正交表,进行表头设计
正交表的选择是实验设计的首要问题。
正交表选得太小,实验因素可能安排不下;正交表选得过大,实验次数增多,不经济。
正交表的选择原则是在能够安排下实验因素和交互作用的前提下,尽可能选用较小的正交表,以减少实验次数。
表头设计就是指将实验因素和交互作用合理地安排到所选正交表的各列中去的过程。
若实验因素间无交互作用,各因素可以任意安排;若要考察因素间的交互作用,各因素应按相对应的正交表的交互作用列表来进行安排,以防止设计“混杂”。
正交表是一整套规则的设计表格,表示方法记为Lm(rn),其中L为正交表的代号,m为实验的次数,r为水平数,n为列数,也就是可能安排最多的因素个数。
例如L9(33),(表1),它表示需作9次实验,最多可观察3个因素,每个因素均为3水平。
一个正交表中也可以各列的水平数不相等,我们称它为混合型正交表。
根据正交表的数据结构看出,正交表是一个m行n列的表,表11中,大写字母表示A、B、C三个因素,下标表示每个因素有1、2、3三个水平。
表1三因数三水平正交实验表
Table1Theorthogonaltesttableofthreefactorsandthreelevels
A
B
C
效益平均值Y
1
A1
B1
C1
2
A1
B2
C2
3
A1
B3
C3
4
A2
B1
C2
5
A2
B2
C3
6
A2
B3
C1
7
A3
B1
C3
8
A3
B2
C1
9
A3
B3
C2
正交表具有以下两项性质:
(1)每一列中,不同的数字出现的次数相等。
例如在三水平正交表中,任何一列都有下标“1”、“2”、“3”,且在任一列的出现次数均相等。
(2)任意两列中数字的排列方式齐全而且均衡。
例如在三水平情况下,任何两列(同一横行内)有序对共有9种,1.1、1.2、1.3、2.1、2.2、2.3、3.1、3.2、3.3,且每对出现数也均相等。
以上两点充分的体现了正交表的两大优越性,即“均匀分散性,整齐可比”。
通俗的说,每个因素的每个水平与另一个因素各水平各碰一次,这就是正交性。
2.3.4编制实验方案,按方案进行实验,记录实验结果
根据实验方案进行实验,记录实验结果。
2.4正交实验结果分析
极差分析法(R法)又称直观分析法,此法计算简便而直观,简单、易懂,是正交实验结果分析最常用的方法。
2.4.1确定实验因素的优水平和最优水平组合
最佳水平是指每个因子的各水平中使指标达最佳的水平。
为确定因子的最佳水平,必须确定该因子各水平对指标的影响。
分析A因素各水平对实验指标的影响。
A因素的i水平所对应的实验指标平均值用kAi表示(Kij—第j个因素第i个水平的所有实验结果指标值的均值)。
根据正交设计的特性,对A1、A2、A3来说,三组实验的实验条件是完全一样的(综合可比性),可进行直接比较。
如果因素A对实验指标无影响时,那么kA1、kA2、kA3应该相等,但kA1、kA2、kA3实际上不相等。
说明,A因素的水平变动对实验结果有影响。
因此,根据kA1、kA2、kA3的大小可以判断A1、A2、A3对实验指标的影响大小,k值大的对实验指标的影响大,为A因素的优水平。
同理,可以计算并确定B、C因素的优水平。
三个因素的优水平组合为实验的最优水平组合。
2.4.2确定因素的主次顺序
各因子对指标的影响是不同的,其重要性也各不相同。
为了评价各因子的重要性,需拟定一评价指标。
通常采用均方和S或极差R作为评价指标。
Sj—均方和
Sj=[(K1j-K2j)2+(K1j-K3j)2+(K2j-K3j)2]
MAX=MAX(K1j,,K2j,K3j)
MIN=MIN(K1j,,K2j,K3j)
极差=MAX-MIN
根据极差Rj的大小,可以判断各因素对实验指标的影响主次。
极差Rj各数据间的差距越大,说明该因子各水平相差悬殊,对指标的影响大,反之则小。
因此以均方和或极差可粗略的揭示出各因子的重要性。
2.4.3因子显著性的检验
因子的重要性只说明该因子相对其他因子的重要程度,而未说明该因子对指标影响的显著程度。
如果某因子对指标的作用不显著,则可排除该因子而使决策简化。
经显著性检验之后,可确定对指标有显著影响的因子、排除对指标影响不显著的因子。
在此基础之上可选择与确定最佳方案。
当计算出的F值大于临界值,k因子在α水平下作用显著,否则作用不显著。
最佳方案的确定方法是选择对指标有显著影响的因子中的最佳水平,对于对指标无显著影响的因子可不考虑,或根据实际情况决定。
2.4.4绘制因素与指标趋势图
以各因素水平为横坐标,实验指标的平均值(kjm)为纵坐标,绘制因素与指标趋势图。
由因素与指标趋势图可以更直观地看出实验指标随着因素水平的变化而变化的趋势,可为进一步实验指明方向。
2.5正交实验总结
1在实际生产中,最优条件的确定是灵活的。
对于主要的影响因素,一定要选最优水平,而对于次要因素,则应权衡利弊,综合考虑来选取优水平。
2极差分析得到的最优工艺条件并不一定在所实施的正交实验方案中。
为了考察最优条件的再现性,应追加验证性实验,从而进一步判断研究所找出的生产工艺条件是否最优。
3由极差分析得出的最优工艺条件,只有在实验所考察的范围内有意义。
4正交实验的极差分析简明、直观,是有效的数理统计方法,应提倡大力推广。
3.1材料与方法
3.1.1原料与仪器
甲酸(化学纯,85%)、无水碳酸钾(化学纯,99%)、蒸馏水、邻苯二甲酸氢钾(基准试剂,99.8%)、氢氧化钠(化学纯,96%)、酚酞;
三颈烧瓶,冷凝管,恒温水槽,磁力搅拌器,抽滤漏斗,真空干燥箱,电子天平。
3.1.2实验原理
甲酸与碳酸钾按4:
1的物质的量比在一定温度、时间条件下反应后,通过降温结晶,分离,干燥,得到的白色晶体二甲酸钾。
其化学反应
4HCOOH+K2CO3=2KH(HCOO)2+CO2+H3O
二甲酸钾的生产工艺流程图如下
3.1.3产品标定方法的确定
二甲酸钾纯度的测定有反相高效液相色谱法和酸碱滴定法。
反相高效液相色谱法是以二甲酸钾标准品的HPLC图谱和二甲酸钾产品HPLC图谱比较,运用峰高外标定量法即可计算出产品纯度;酸碱滴定法是用已知浓度的氢氧化钠溶液滴定二甲酸钾,其反应如下:
HCOOH-HCOOK+NaOH=NaCOOH-HCOOK+H2O
二甲酸钾与氢氧化钠是一比一反应,称取一定量二甲酸钾加水溶解,用已知浓度的氢氧化钠溶液滴定,滴入两滴酚酞做指示剂,到滴定终点记下所用氢氧化钠的体积,算出所用氢氧化钠的量就能算出二甲酸钾的量,从而可算出二甲酸钾的纯度。
综合考虑到实验条件和实验准确性,本实验用滴定法测定二甲酸钾纯度。
3.2实验步骤
3.2.1配制标准氢氧化钠溶液
用电子天平称取二甲酸钾产品,从准确性考虑:
量不能太少。
若定为取产品1.5g,二甲酸钾分子量为130.12,则滴定所用氢氧化钠的量为1.5g÷130.12g/mol=0.011528mol,用50ml碱式滴定管滴定则需用25.00ml,所配氢氧化钠溶液的浓度应该为0.01537mol÷0.025l=0.4611mol/l,九次实验,每次产品滴定3次至少要滴定27次,一次滴定要25.00ml,需用1000ml的容量瓶配制。
氢氧化钠纯度为96%,总共需称取粗氢氧化钠:
0.4611mol/L×1L×40g/mol÷0.96=19.21g.
3.2.1.1配制浓度为0.5mol/L的氢氧化钠溶液
把烧杯放在天平上,回零后称取氢氧化钠约20g,用蒸馏水溶解后,冷却到室温,用玻璃棒引流转移到1000ml容量瓶,玻璃棒和烧杯用蒸馏水洗三遍移入容量瓶,然后加蒸馏水至刻度线,塞紧瓶塞后摇匀即可。
3.2.1.2标定氢氧化钠溶液
取在105℃干燥至恒重的基准邻苯二甲酸氢钾约2g,精密称定质量为m=2.013克,加蒸馏水振摇,使其尽量溶解,加酚酞指示液2滴,用本液滴定,在接近终点时,溶液由无色显粉红色。
记录氢氧化钠溶液的体积v=22.30ml,计算氢氧化钠的准确浓度为0.44113mol/l,公式为CNaOH=m×99.8%÷M×1000÷V。
其中,m邻苯二甲酸氢钾质量,C氢氧化钠摩尔浓度,M氢氧化钠摩尔质量,V氢氧化钠溶液体积.
做三次平行实验,把实验数据及计算结果填入表2。
表2氢氧化钠标定
Table2Thecalibrationofsodiumhydroxide
邻苯二甲酸氢钾质量
m(g)
氢氧化钠溶液体积
V(ml)
氢氧化钠溶液浓度
C(mol/l)
2.013
22.30
0.44113
2.037
22.75
0.43756
2.026
22.70
0.43616
计算氢氧化钠溶液平均浓度:
CNaOH=(0.44113mol/l+0.43756mol/l+0.43616mol/l)÷3=0.4383mol/l
3.2.2合成二甲酸钾
本实验考察甲酸与碳酸钾合成二甲酸钾,考察指标为二甲酸钾的收率。
查二甲酸钾合成专利及各类文献,可知甲酸与碳酸钾合成二甲酸钾影响二甲酸钾收率的主要因数为反应摩尔比、反应温度、反应时间,其最优条件范围为反应摩尔比(甲酸:
碳酸钾=4.0~5.0),反应温度(50℃~60℃),反应时间(1.5h~3h)。
因此做三因数三水平正交实验,各因数的三个水平如下所示:
表3实验因数与水平
Table3experimentalfactorsandlevels
摩尔比(甲酸:
碳酸钾)
反应温度(℃)
反应时间(h)
1
4.0:
1.0
50
1.5
2
4.5:
1.0
55
2.0
3
5.0:
1.0
60
2.5
表4.摩尔比例与质量对应表
摩尔比
甲酸用量(85%)(g)
碳酸钾用量(99%)(g)
4.0:
1.0
64.9
41.882
4.5:
1.0
73.1
41.882
5.0:
1.0
81.2
41.882
具体实验正交表安排如下表:
表5实验正交表排列
Table5Thetableoforthogonaltest
A
摩尔比
B
反应时间(h)
C
反应温度(℃)
收率
Y
1
4.0:
1.0
1.5
50
2
4.0:
1.0
2.0
55
3
4.0:
1.0
2.5
60
4
4.5:
1.0
1.5
55
5
4.5:
1.0
2.0
60
6
4.5:
1.0
2.5
50
7
5.0:
1.0
1.5
60
8
5.0:
1.0
2.0
50
9
5.0:
1.0
2.5
55
为了更好更快完成实验,排除外界环境因数对实验产生误差,先把以上九次实验按反应温度分成三组:
表6.实验分组表
因子分组
A
摩尔比
B
反应时间(h)
C
反应温度(℃)
第一组
4.0:
1.0
1.5
50
4.5:
1.0
2.5
5.0:
1.0
2.0
第二组
4.0:
1.0
2.0
55
4.5:
1.0
1.5
5.0:
1.0
2.5
第三组
4.0:
1.0
2.5
60
4.5:
1.0
2.0
5.0:
1.0
1.5
3.2.2.1实验步骤
(一)接通电源,打开恒温水浴,把温度设定到50℃。
(二)分别取三个大烧杯三个小烧杯编号并称重,用小烧杯称取41.882克99%的碳酸钾试剂三份,用大烧杯称取64.9克、73.1克、81.2克85%的甲酸试剂各一份。
(用天平称取干燥时的烧杯重后,加入所需称量的试剂至所需量。
)
(三)向甲酸里缓慢加入碳酸钾,为减少甲酸挥发用玻片遮着点;反应过程中不断搅拌,一是可以使反应充分,二是可以使温度不太高。
(四)反应后,把三个烧杯一起放入恒温水浴,按表6所列反应时间依次取出各烧杯,在同一条件下冷却到室温结晶。
(五)安装好吸滤瓶与布氏漏斗,剪两张比漏斗内颈稍小的圆形滤纸。
室温为16℃时,启动循环水真空泵,滤纸用结晶产物中的少许液体润湿并将其贴在布氏漏斗内,放在吸滤瓶上,减压吸紧,然后将结晶产物倒入布氏漏斗中,十分钟后将吸滤瓶中的母液倒入原来装结晶产物的烧杯中清洗后倒入漏斗中再抽滤,母液循环清洗、抽滤至烧杯内壁所附产物很少时,抽滤完后停泵,将漏斗内的产物移至一干燥烧杯中。
3.2.2.2用氢氧化钠滴定产品
干燥大约24小时,称量产品重量直至重量不再减少,可认为产品已经干燥,记录产品重量;由反应式计算出二甲酸钾理论值,然后计算出二甲酸钾的收率,见表7。
收率=实际产品值÷理论产品值×100%
用电子天平准确称取制得的二甲酸钾晶体1.2-1.3g,溶解在80ml蒸馏水中,然后用已标定浓度的氢氧化钠溶液滴定,酚酞做指示剂,计算出二甲酸钾的纯度,见表8。
表7二甲酸钾收率计算结果
Table7Thereceivingrateofpotassiumdiformate’spurities
实验号
原烧杯重(g)
烘干后总重(g)
实际产品值(g)
理论产品值(g)
收率
1
165.201
226.018
60.817
78.730
0.772
2
89.333
144.901
55.568
78.722
0.706
3
151.875
203.675
51.800
78.730
0.658
4
190.037
237.971
47.934
78.730
0.609
5
73.881
120.289
46.408
78.704
0.590
6
136.303
182.389
46.086
78.780
0.585
7
179.215
213.613
34.398
78.707
0.437
8
81.389
114.931
33.542
78.737
0.426
9
194.110
226.480
32.370
78.760
0.411
表8二甲酸钾纯度的计算结果
Table8Thecalculationofpotassiumdiformate’spurities
实验号
二甲酸钾质量(g)
氢氧化钠体积(ml)
二甲酸钾纯度
收率
平均纯度
1
1.322
22.46
0.971
0.773
0.968
1.330
25.06
0.965
1.326
27.60
0.968
2
1.318
25.18
0.973
0.706
0.972
1.323
22.76
0.970
1.319
25.32
0.973
3
1.277
21.90
0.978
0.658
0.976
1.265
21.66
0.977
1.298
22.14
0.973
4
1.210
20.60
0.971
0.609
0.971
1.206
20.48
0.968
1.213
20.70
0.973
5
1.192
20.48
0.980
0.59
0.981
1.226
21.10
0.982
1.278
22.00
0.982
6
1.179
20.28
0.981
0.585
0.983
1.169
20.16
0.984
1.232
21.28
0.985
7
1.218
20.94
0.980
0.437
0.982
1.236
21.24
0.980
1.222
21.10
0.985
8
1.324
27.88
0.969
0.426
0.971
1.318
25.46
0.973
1.321
30.04
0.971
9
1.286
22.08
0.979
0.411
0.979
1.295
22.22
0.979
1.313
22.56
0.980
4实验结果分析
4.1对正交实验表分析
把实验结果填入正交实验表,用Excel计算各数据值。
表9正交实验分析结果
Table9Theanalysisoforthogonaltest
AHCOOH:
K2CO3
B
反应时间(h)
C
反应温度(℃)
D
收率(%)
1
4.0:
1.0
1.5
50
77.3
2
4.0:
1.0
2.0
55
70.6
3
4.0:
1.0
2.5
60
65.8
4
4.5:
1.0
1.5
55
60.9
5
4.5:
1.0
2.0
60
59.0
6
4.5:
1.0
2.5
50
58.5
7
5.0:
1.0
1.5
60
43.7
8
5.0:
1.0
2.0
50
42.6
9
5.0:
1.0
2.5
55
41.1
K1
71.23
60.63
59.47
K2
59.47
57.40
57.53
K3
42.47
55.13
56.17
S
1254.98
45.84
16.50
MAX
71.238
60.63
59.47
MIN
42.47
55.13
56.17
R
28.77
5.50
3.30
优水平
A1
B1
C1
表中K1,K2.K3,表示每个因素所进行的3次实验(3水平)所得二甲酸钾收率的平均值,即在代表因素A(摩尔比)的第一列中,将与水平1对应的l,2,3共3个二甲酸钾收率相加除3所得,其他类推。
Ri为平均数的极差(K中最大值与最小值之差)。
Ri愈大表明变化幅度大,说明该因素的水平变化对二甲酸钾收率的影响愈大,该因素即是最大因素。
注:
针对每个因素,其哪个水平的K值大,该水平即为该因素的最优水平。
其中,针对每个因素
Kij—第j个因素第i个水平的所有试验结果指标值的均值;
Sj—均方和
Sj=[(K1j-K2j)2+(K1j-K3j)2+(K2j-K3j)2]
MAX=MAX(K1j,,K2j,K3j)
MIN=MIN(K1j,,K2j,K3j)
极差=MAX-MIN
由上表4.1所列结果可知,就二甲酸钾的收率而言,A因素的第1个水平(物质量的比4.0:
1.0)、B因素的第1个水平(1.5h)、C因素的第1个水平(50℃)为最优工艺条件。
由表4.1中Rj的值还可以看出,影响二甲酸钾收率的因素主次依次排列为:
摩尔比>反应时间>反应温度,试验的最优水平组合为A1B1C1。
4.2绘制因素与指标趋势图
以各因素水平为横坐标,试验指标的平均值(kjm)为纵坐标,绘制因素与指标趋势图。
由因素与指标趋势图可以更直观地看出试验指标随着因素水平的变化而变化的趋势,可为进一步试验指明方向。
图5摩尔比因数的影响
Figure5Theinfluenceofthemonomermoleratio
随着甲酸与无水碳酸钾物质的量的比的增加,二甲酸钾的收率不断降低,这是由于增加甲酸过量会减少二甲酸钾从反应生成物中结晶析出数量。
其原因是没反应的甲酸对二甲酸钾的结晶析出起到阻碍作用。
摩尔比是影响二甲酸钾收率的主要因数。
图6反应时间因数的影响
Figure6Theinfluenceofreactiontime
反应时间对产品收率影响较小,主要是由于本反应属于简单的酸碱中和反应,反应较为迅速而剧烈,瞬间即可完成。
随着反应时间增加,产品收率反而减小,主