稳健设计方法.docx

1、稳健设计方法机研142 孙利文 2100一、稳健设计方法在产品设计开发中的作用。稳健设计又称作鲁棒设计。是关于产品质量和成本的一种工程设计方 法。在产品或工艺系统设计中，正确的应用稳健设计的基本理论和方法可 以使产品在制造或使用中，或是在规定寿命内当结构或材料发生老化、变 质、工作环境发生微小的变化时，都能保证产品质量的稳定。通过稳健设计，可以使产品的性能对各种噪声因素的不可预测的变 化，拥有很强的抗干扰能力。产品性能将更加稳定、质量更加可靠。任何 一种产品影响其质量的因素有很多主要可分为两类：一类是在设计中 人们可以控制的因素如设计变量、变量的容差等；另一类是所谓的噪声因 素指由生产条件、使

2、用环境及时间等的变化而影响产品质量的因素如载 荷、儿何尺寸、工程材料特性的变异以及制造、安装误差等，其基本特点 是具有不确定性和随机性，是不可控制的因素。实际存在的不确定因素的 变化有可能导致产品的性能指标有较大的波动，使其功能劣化我至失效， 还有一些材料或元器件会随着时间的推移而发生失效等。对于这些因素有 两种处理方法:一是尽可能消除这些因素 这对可控因素是可以做到的，而 对噪声因素往往很难实现。即使能够消除也需要花费很大的代价；二是尽 量降低这些因素的影响。这是相对容易和低代价的方法，也就是使产品性 能对这些因素的变化不敏感，为了使所设计的产品在不确定因素的影响 下，其性能指标不仅能达到设

3、计要求，而且对各种不确定因素的变化不敏 感，就需要用稳健设计方法來实现。稳健设计就是使产品的性能对在制造 期间的变异或使用环境的变异不敏感，并使产品在其寿命周期内不管其参 数、结构发生漂移或老化在一定范围内都能持续满意地工作。二、 试验设计在稳健设计中的作用。试验设计就是运用正交试验法或优化方法确定零部件参数的最佳组 合，在系统内、外因素作用下，所产生的质量波动最小，即质量最稳定（健 壮）。试验设计的目的是根据系统设计中所确定的所有参数，通过多因素的 优选方法来考察三种干扰（内干扰，外干扰，产品间波动）对系统质量特 性的影响，寻求最佳的参数组合，以求得抗干扰性最佳的设计方案，使系 统质量特性波

4、动小，稳健性好，并且价格低廉。1） 通过对产品特性指标试验结构的分析，确定出在不可控因素影响 下的可控因素值的最佳组合或工作状态；2） 筛分出各可控因素对产品质量影响的重要程度，以利于有效控制;3） 可以实现原有材料代用、基本设计构型的比较、设计参数的选择 以保证产品质量在较宽的工作条件下具有稳定性；4） 用于建立经验模型、确定质量特性与设计变量（或工艺参数）之 间的函数关系，寻找最佳工作条件。三、 按照正交试验（直观分析法）的原理，解决你实际工作中的一个问题, 并总结成实验分析报告（以两个以上因素分析为例）。正交试验设计是用来科学地设计多种因素试验的一种方法，它利用一 套规格化的正交表安排试

5、验，得到的实验结果用数理统计进行处理，以得 到科学结论。正交表是试验设计的基本工具，它是根据均衡分布的思想， 运用组合数学理论构造的一种数学表格，均衡分布性是正交表的核心。正交试验设计（简称正交设计）的基本程序是设计试验方案和处理试 验结果两大部分。主要步骤可归纳如下：（1） 明确试验目的，确定考核指标（2） 挑因素，选水平。（3） 选择合适的正交表，进行表头设计。，（4） 确定试验方案。（5） 试验结果分析实例：某化工厂生产一种试剂，产率较低，希望通过试验探索好的生产工艺 以提高产率。考察的因子与水平如下表：表8-10考查试剂生产效率因素一览表因子ABC水平反应温度反应时间搅拌速度/摄氏度/

6、小时一水平301快二水平40lo 5中三水平502慢这是一个三水平的试验，我们可以在（罗）和岛仔）中选一张合适的 表。选择的原则是在试验因子能在正交表的列中安排得下的前提下，试验 次数越少越好。本例只有三个因子,故选用厶（扌）表，作9次试验就行了选 择了正交表后，将因子安排在矗厅）的表头上,我们将三个因子依次安排在1, 2, 3列，并且把表中各列的水平号用相应的实际因子水平写出来,就得到-张试验设计表。列123号反应温反应时间搅拌试验度(A)/摄氏(B)/小时速度(C)号度111 1130快212230lo 5中313 2330慢421 1240中522 lo3405慢6231402快7313

7、501慢8321表8-12计算表列号试验号1反应温度/摄氏度2反反应时间(B)/小时3搅拌速度(C)试验结 果产 率11(30)1 (1)1(快)82212 (lo5)2仲)81313376（慢）42(40)128052238562318273(50)1364832272933164I（一水平试验结果总和）239226236II（二水平试验结果总和）247238225III（三水平试验结果总和）2002222251/379o 775o 378o 7II/382o 379o 375o 0III/366o 774 o 075o 0极差R15o 65o 33o7按以上所设计的方案进行了 9次试验后,

8、将各次试验结果依次填入试 验计划表的最右边，并且在表上进行了 一系列的计算，形成了上表形式, 常称为计算表。现在根据这9次试验结果,来分析因素各水平对产率的影响。先看A因子（反应温度）。它的水平为30摄氏度的是第1, 2, 3号试验, 其总产率IA二82+81+76二239;它的水平是40摄氏度的是第4,5,6号试验, 其总产率IIA二80+85+82二247;它的水平是50摄氏度的是第7,8,9号试验, 其总产率 IIIA 二64+72+64二200。在A因子水平相同的三组试验中，不同水平的B因子（反应时间）和不 同水平的C因子（搅拌速度）都各出现一次。从整体上看，可以认为B,C两 因子对产

9、率的影响虽然在变动，但这种变动是均衡的。因此，比较这三个总 产率，就可以看出A因子各水平的差别对产率的影响。为便于说明,把上述 三个总产率都取平均值，分别得到IA/3二，IlA/3=, IIIA/3=,这是试剂的平均产率。显然A因子取40摄氏度最好50摄氏度最差。二者之差即极差J = 823 66.7二15.6 ,它表示反应温度40摄氏度与50摄氏度相比，试剂的 产率平均要提高。用同样的方法可以比较B因子和C因子的各水平的好与差。比较各因子极差的大小，就可以看出哪个因子对产率的影响大，哪个 因子影响小。反应温度的高低对试剂的平均产率的影响可以差到%，而搅拌 速度的快慢对试剂的平均产率的影响只差

10、到，显然反应温度是否合适要 比搅拌速度是否合适重要的多。根据这种比较，就可以回答本节开始提出 的三个问题了。(1)反应温度对产率影响最大，其次是反应时间，再其次是搅拌 速度。(2)反应温度是40摄氏度好，反应时间是小时好，搅拌速度是快 速好。(3)好的生产工艺是：&禺G?r卩反应温度40摄氏度；反应时间小时；搅拌速度快速。&色】这个条件在试验计划表中并没有出现，它是27次全面试验中的 一种。由此可见，用正交表安排试验确实具有很强的代表性。虽然只作了 9 次试验，但是通过对这9次试验结果的计算与分析，仍然不会漏掉最佳的 水平组合。以上利用比较各因子不同水平下试验结果平均值的方法就是直观分 析法，

11、也叫做综合比较法。显然，只有在均衡搭配的试验情况下，才能进行 综合分析，这也是正交表的一个特性,常称为“综合可比性”。四、举例说明稳健设计在机械工程中的应用(包括求解的过程)在汽车鼓式制动器设计中考虑不确定因素的影响，将可靠性优化理 论、可靠性灵敏度分析与稳健设计方法相结合，以制动效能因数为目标函 数建立制动器可靠性稳健优化数学模型。把制动力矩、摩擦衬片压力的可 靠性灵敏度溶入可靠性优化设计模型之中，将可靠性稳健优化设计转化为 满足可靠性要求的多目标优化问题。实例计算表明，稳健优化后的制动器 不仅有较高的制动效能和可靠性，还具有较低的可靠性灵敏度，取得了满 意的结果。1鼓式制动器的状态方程为避

12、免制动过程中车轮打滑，制动力矩不得超过车轮与地面的附着力矩, 一般希望车轮与地面的附着系数小于规定值。所以根据文献2可以推导 出制动力矩的状态方程为g： (x, y)=019Mn / mRl (1)Mx= (MfhA) / (R(cosB+fsinB)-fA)M：= (MfhA) / (R(cosB-fsinB)+fA)P =tan ( (cos2 a lcos2 a J / (23-sin2 a ：+sin2 u i)AMR (cos a i一cos a 2) / (cos2 a icos2 a 2)+ (2 a 3一sin2 a:+sin2a1)2)1/2式中,M为制动力矩,M-M1+M2

13、； r为制动鼓半径,R为蹄片支承销中心与制动鼓中心间的距离,h为制动蹄轴端至末端的距离，a - a 2分别是摩擦衬片的起始与终止点和鼓心连线的夹角，a 3= a 2- a 1为摩擦衬片的包角,b为制动鼓宽度,f为制动鼓与摩擦衬片间的摩擦系数，F为制动蹄促动力,m为汽车总质量，n为车轮数或制动器数,R1为午轮半径。摩擦衬片上承受的最大压力应小于规定值，根据文献2可以得到摩擦衬片的状态方程为3-4rfhsin ( a s/2) 2可靠性稳健优化设计由可靠性优化设计的基本理论可知，可靠性设计的目标是计算可靠度R=Qg (x) 0f x (X) dX(4)式中 fx(X)为基本随机参数向量 X=(x,

14、 y) = (xl, x2, xn; yl, y2, ym)T的联合概率密度，Xi 是设计变量(可控变量)，yi是设计参数(不可控变量)，X和y分别是设计变 量xi和设计参数yi构成的向量;g (X)为状态函数,可表示零部件的两种状态。g (x, y) F0为失效状态g(x, y) 0为安全状态(5)可靠性指标定义为3B=LgRg=Eg(X)Varg(X)(6)式中 Lg=E g (X) =g (X)和 Rg=Var g(X)=5g (X) 5XT2Var (X)分别是状态函数 g (X)的均值与标准差,Var (X)为基本随机参数的方差向量。这样一方面可以利用可靠性指标直接衡量构件的可靠性，

15、另一方面在基本随机参数向量X服从正态分布时，可以用失效点处状态表面的切平面近似地模拟极限状态表面，可以获得可靠度的一阶估计量R二5(B) (7)式中5(#)为标准正态分布函数。可靠度对基本随机参数向量X均值和方差的灵敏度为4 dRdXT=5R5B5B5Lg5Lg5X2(8)dRdVar(X)=5R5B5B5Rg5Rg5XT (9)将以上各式和已知条件代入式(6)和式(7),在基本随机参数服从正态分布的情况下,就可以求出可靠性指标B和可靠度R。可靠性稳健优化设计的基本思想是在可靠性优化设计模型的基础上，把可靠性灵敏度加到目标函数中，考虑约束的可行稳健性，将可靠性稳健优化设计归结为满足可靠性要求的

16、多目标优化问题。其数学模型可描述为minf 1 (?x, ?y)minf2 (?x, ?y)=Eni=l5R5xi2s. t. g(?x, ?y) -5-l(R0)Rg0qj(?x, ?y)+kE5qj5xi2R2xi+E5qj5yi2R2yi1/20(j=l, 2, m)xlixi3Rxixui (i=l, 2, n)(10)式中fl (?x, ?y)为原优化问题的目标函数,f2 (?x, ?y)为可靠度对设计变量x=(xl, x2, xn)T均值的灵敏度的平方和再开方,Rgl(x, y) = O19-Mn/niRl 0为设计所要求的可靠度,QJ (?x, ?y)为不等式约束函数,?x ?y

17、分别是设计变量和设计参数的均值,R2xi、R2yi分别是设计变量xi和设计参数yi的方差。3鼓式制动器的可靠性稳健优化设计311目标函数评价制动器性能好坏的最主要指标之一是制动器效能因数,它表征了制动器将一定大小的制动蹄促动力转化为制动器制动力矩并进而转化为地面制动力的能力。在相同的促动力下，效能因数越大,表明制动力矩越大,制动效果越好，工作效率越高。效能因素表示为K=M1+M2FR(11)因此，目标函数为fl(?x, ?y)=l/K (12)取设计变量 x= (r, R, h, Al, A2)T= (xl, x2, x3, x4, x5)T,设计参数 y二(m,Rl, f)T=(yl, y2

18、, y3)T。另外，要求制动器的可靠度对设计变量均值的灵敏度最小，则f2 (?x, ?y) =E5i=l5R5xi312约束条件(1)制动力矩的约束条件gl (?x,?y)-5-1(RO)Rgl0(14)(2)摩擦衬片压力的约束条件g2 (?x,?y)-5-1(RO)Rg20(15)g3 (?x,?y)-5-1(RO)Rg30(16)(3)自锁条件ql (x, y)=f-RcosBA-RsinB0(17)(4)压力分布均匀约束q2(x, y)=f-RA3hsin(A3/2) 0(19)(5)摩擦衬片的磨损特性约束q3(x, y) = 0(20)式中v为汽车制动初速度，t制动时间。(6)设计变量

19、边界约束综合以上目标函数和约束条件，根据式(10),可得制动器的可靠性稳健优 化模型。4实例计算某型汽车质量m为(Lm,Rm)二(3 500, 175) kg,轮胎半径为(LR1, RR1) = (01343, 01001 72)m,制动气压0153901588MPa,制动初速度v为(Lv, Rv) = (80, 4) km/h,制动鼓与摩擦衬 片间的摩擦系数f为(Lf,Rf) = (014, 0102)，设计可靠度要求R0为01999;设计变量的边界 为:0U6r0119 m,0U2R0118m, 0126h0138 m, O115A1 0180 rad, 1157A22109 rad,可认

20、为它们都服 从正态分布、相互独立的。根据加工公差和正态分如的3R法则，各设计变量的均方差可取相应均值的 015%3o根据建立的优化模型，按照相容决策支持问题法5,利用基于Matlab的序惯二次规划法(SQP)求解，求得可靠性稳健优化设计的最优解x*o可靠性稳健优化计算结果与常规优化及初始设计值见表1。表1优化结果xl/mx2/mx3/mx4/radx5/radfl(X)f2(X) /10-3初始设计值 0118 01128 0127 01262 2109 21151 41684 8常规优化 01188 01122 01302 01175 21007 31610 61168 4可靠性稳健优化 01187 01122 01314 01182 21059 31583 01954 3由表1可知，常规优化所得的制动器效能因数值要小一些,但其可靠度灵敏度明显偏大，而且有儿个约束处于边界上，一旦设计变量发生波动，就可能违反约束。而采用可靠性稳健优化方法虽然得到的制动器效能因数要稍大一些，但设计质量具有较好的不灵敏性，约束具有可行稳健性，确保了制动器的可靠性和稳健性。比较可知，可靠性稳健优化后的制动器结构参数比常规优化后 的更加合理。