试验设计实验五.docVIP

实验5参数设计

一、实验目的

掌握参数设计的根本原理及实施步骤，能够通过信噪比等一些指标对试验结果进行分析，给出合理结论，能够利用SAS程序得到数据所用到的指标,学会对运行结果进行分析。

二、实验内容

1、〔可以考虑综合噪声因子〕某气动换向装置的设计中，关键是要使换向末速度V到达V0=960mm/s。据力学原理，换向末速度的运动方程式为

V=

其中A为换向行程〔待选参数〕，B为换向活塞直径〔待选参数〕，C为气缸内气压〔待选参数〕，D为换向阻力〔不可控，但可测，大约在750N±20N〕,E为系统重量(可在90kg±5kg取值),g=9800mm/s2为重力加速度.

在这个问题中可控因素、噪声因素及其水平见下表.

因素水平表

水平

1

2

3

可控因素

换向行程A/mm

52

56

60

换向活塞直径B/mm

22

24

26

气缸内气压C/(N/mm2)

2.2

2.6

3.0

噪声因素

换向行程A′/mm

A-0.2

A

A+0.2

换向活塞直径B′/mm

B-0.1

B

B+0.1

气缸内气压C′/(N/mm2)

C-0.2

C

C+0.2

换向阻力D

730

750

770

系统重量E/kg

85

90

95

请作参数设计，使V尽量接近目标值V0，而使波动小。

〔1〕内外表设计

把可控因素A、B、C放在正交表的第1，2，3列上,由于噪声因素比拟多，外表用综合噪声因子，由换向末速度的运动方程式可知，V是A、B、C的增函数，又是D、E的减函数，因此其综合噪声因子的两个水平分别为：

负侧最坏水平，

正侧最坏水平。

外表用此综合噪声因子，那么总试验次数只有18次，在内表的水平组合下把综合因素的两个水平带入公式，如内表第1号试验中可控因子A、B、C均取1水平,即,于是根据题目表立即可得到综合噪声因子的两个水平如下:

〔数据有误，设〕

在工作表中计算V和统计分析过程及结果如下：

在方差分析表中由于MSAMSE，故将MSA归入误差项得到新的MSE，根据方差分析的结果,因素B、C均不显著，但可确定因素的主次顺序为C、B高度显著的因子C应选使信噪比最大的,次显著因子B应选使信噪比最大的.这样一来,最正确水平组合应是,它是内表的第3号试验,该号的信噪比是9个试验中最大的.

灵敏度设计

第3号试验的可控因子水平组合使波动减少,但其均值与目标值还有一定的距离.这需要通过灵敏度分析来解决.

(2)在SAS中输入:

datacs;

inputobsabcdy1y2@@;

vm=(y1+y2)/2;

sn=10*log(2*y1*y2/(y2-y1)**2);

cards;

1 1 1 1 1 700 483.2873647

2 1 2 2 2 577.8193007 819.3768627

3 1 3 3 3 876.4763387 1097.790454

4 2 1 2 3 396.1407616 673.9892671

5 2 2 3 1 755.2203243 981.8290646

6 2 3 1 2 576.7601136 855.447042

7 3 1 3 2 598.1397521 838.9423648

8 3 2 1 3 399.5424241 718.2922821

9 3 3 2 1 788.4642126 1042.524478

;

run;

procprintdata=cs;

run;

procglmdata=cs;

classabcd;

modelvm=abc;

meansabc;

run;

结果如下：

输出结果1.SEQ输出结果1.\*ARABIC1方差分析表

结果分析:

由图4均值的方差分析表可知因素的主次顺序为B、C.结合稳健设计的结果可得可控因子分类表如表3,结果是B是比拟重要因子,C是稳健因子,于是可用力B这个因子使接近目标值V=960m/s

输出结果1.SEQ输出结果1.\*ARABIC2A、B、C对的影响

表1可控因子分类

类别

稳定性分析

灵敏度分析

因子

=1\*ROMANI

**

*

C重要因子)

=2\*ROMANII

*

B(调节因子)

=3\*ROMANIII

A(次要因子)

=4\*ROMANIV

无

表2信噪比的表

水平

A

B

C

1

30.41792

24.84277

23.13748

2

26.13041

26.2713

26.52736

3

26.06727

31.50152

32.95076

表3的表

水平