零件的参数设计
摘要：本文通过研究成批生产粒子分离器某参数与7个零件参数的标定值和容
差的关系，建立了一个有约束条件的非线性规划模型。
首先，对参数y的经验公式进行线性化处理，确定了参数y近似服从正态分布，在此基础上, 以7个零件参数的标定值和容差等级为决策变量, 以生产产品的总费用为目标函数建立非线性规划模型,用matlab软件编程计算，对108种容差等级组合逐一进行求解，选取最优的参数标定值和容差等级组合，使得生产1000个产品的总费用最低，最后求得7个零件参数的标定值分别为0.0750，0.3750，0.1143，0.1200，1.2903，13.4292，0.5791，容差等级组合为B,B,B,C,C,B,B，总费用是42.2073万元，与原设计费相比，降低了265.4842万元。然后用蒙特卡洛随机数法对模型进行检验，模拟结果吻合很好。最后对模型进行了误差分析和灵敏度分析，对产品生产提出了合理化建议。
关键词： 非线性规划 标定值 容差 蒙特卡洛随机数法
一、 问题重述
一件产品由若干零件组装而成，标志产品性能的某个参数取决于这些零件的参数。零件参数包括标定值和容差两部分。进行成批生产时，标定值表示一批零件该参数的平均值，容差则给出了参数偏离其标定值的容许范围。若将零件参数视为随机变量，则标定值代表期望值，在生产部门无特殊要求时，容差通常规定为均方差的3倍。
进行零件参数设计，就是要确定其标定值和容差。这时要考虑两方面因素：一是当各零件组装成产品时，如果产品参数偏离预先设定的目标值，就会造成质量损失，偏离越大，损失越大；二是零件容差的大小决定了其制造成本，容差设计得越小，成本越高。
试通过如下的具体问题给出一般的零件参数设计方法。
粒子分离器某参数（记作y）由7个零件的参数（记作x1,x2,...,x7）决定，经验公式为：
??x4??1?2.62?1?0.36??x???2??x6x7?0.563?x1??x3??Y?174.42???x???x?x1?5??2????0.85?????2?x4????x??2?1.16
y的目标值（记作y0）为1.50。当y偏离y0?0.1时，产品为次品，质量损失为1,000元；当y偏离y0?0.3时，产品为废品，损失为9,000元。
零件参数的标定值有一定的容许范围；容差分为Ａ、Ｂ、Ｃ三个等级，用与标定值的相对值表示，Ａ等为?1%，Ｂ等为?5%，Ｃ等为?10%。7个零件参数标定值的容许范围，及不同容差等级零件的成本（元）如下表（符号／表示无此等级零件）： x1 x2 x3 x4 x5 x6 x7 标定值容许范围 [0.075,0.125] [0.225,0.375] [0.075,0.125] [0.075,0.125] [1.125,1.875] [12,20] [0.5625,0.935] Ｃ等 ／ 20 20 50 50 10 ／ Ｂ等 25 50 50 100 ／ 25 25 Ａ等 ／ ／ 200 500 ／ 100 100 现进行成批生产，每批产量1,000个。在原设计中，7个零件参数的标定值为：x1=0.1,x2=0.3,x3=0.1,x4=0.1,x5=1.5,x6=16,x7=0.75;容差均取最便宜的等级。
请你综合考虑y偏离y0造成的损失和零件成本，重新设计零件参数（包括标定值和容差），并与原设计比较，总费用降低了多少？
二、 模型的假设
（1）各个零件的参数服从正态分布且相互独立；
（2）各个零件参数的容差为均方差的３倍；
（3）经验公式有较高的精度，即不考虑经验公式的误差；
（4）当产品参数y?[1.4,1.6]时，产品无质量损失；当产品参数
y?[1.2,1.4)?(1.6,1.8]时，产品为次品，质量损失为
1000元；当产品参数
元。
y?(??,1.2)?(1.8,??)时，产品为废品，质量损失为9000
三、 符号说明
xi xiy*第i个零件参数的实际值（i?1,2,?,7）； 第i个零件参数的标定值（i?1,2,?,7）； 粒子分离器的某参数值； y0 粒子分离器参数y的目标值； 第i个零件标定值容许范围的下限（i?1,2,?,7）； 第i个零件标定值容许范围的上限（i?1,2,?,7）； n?1表示容差等级为An?2?i ?i Dn 时，最大偏移量与标定值的相对值，即D1?1ò?5%； ； ； 表示容差等级为B时，最大偏移量与标定值的相对值，即时，最大偏移量与标定值的相对值，即n?3表示容差等级为CZD3?10% 1000个产品因参数y偏离y0造成的损失和零件成本的总费用； 7个零件参数的容差等级组合在第k种情况下生产1000个产品的总费用（k?1,2,?,108）； 7个零件参数的容差等级组合在第k种情况下生产一件产品的平均费用（k?1,2,?,108）； 7个零件参数的容差等级组合在第k种情况下每个产品的零件成本（k?1,2,?,108）； Zk zk Ck cki 第i个零件在第k种容差等级组合中的成本（k?1,2,?,108，； i?1,2,?,7）7个零件参数的容差等级组合在第k种情况下平均每个产品质量损失的期望（k?1,2,?,108）； Ek 四、 问题分析
决定粒子分离器性能的某参数y由7个零件的参数决定，这些零件的参数包括标定值和容差两部分。在零件成批生产时，标定值表示一批零件该参数的平均
值，容差则给出了参数偏离其标定值的容许范围。零件的成本只由选择容差等级决定，容差设计得越小，成本越高。零件参数标定值的选择和容差等级的设计会影响到产品参数偏离目标值y0的大小，偏离的越大，损失越大。所以，为使生产1000个产品的总费用达到最小，应统筹规划，选取最优的零件参数标定值和容差。
假设每个零件的参数都符合正态分布，那么通过对参数y进行线性化，可以认为参数y也符合正态分布。根据经验公式，参数y的分布可由7个零件的参数分布求得，从而可以求出平均每个产品质量损失的期望。
现在把生产1000个产品的总费用设定为目标函数，则可通过求解生产每个产品的平均费用来确定总费用最小时的最优解。生产每个产品的平均费用由两部分组成：这个零件的成本和参数偏离目标值y0所造成质量损失的期望。
根据题目中所给的条件，我们知道对于每个零件，不同的容差等级对应不同的零件成本，所以一共有1?2?3?3?1?3?2?108种不同的容差等级组合，这108种容差等级组合对应了108种生产1000个产品的零件成本和产品质量损失的期望。可以先求出每种容差等级组合中7个零件标定值的最优解，再从108组最优解中选取总费用最低的解作为总体的最优解。
五、 模型的建立
5.1 确定参数y的统计特征
首先，为了简化模型，我们需把题目中所给出的关于产品参数y的经验公式进行线性化（线性化的可行性见附录一），将该公式在x?x*处进行泰勒展开，展开式如下：
7Y(x)?Y(x)?*??xi?1?Yxi?xii*7(xi?xi)??(*?(xi?1i?xi)),
*2又因为：
77?i?1(xi?xi)??Y(x)?[Y(x)?*2*??xi?1?Yxi?xii*(xi?xi)],
*所以可以略去佩亚诺余项，得：
7Y(x)?Y(x)?*??xi?1?Yxi?xii*(xi?xi),
*又整理得：
7Y(x)???xi?1?Yxi?xii*7xi?Y(x)??i?1*?Y?xixi?xi*xi ,
*记上式为：
7Y(x)??Mi?1ixi?M0 ,
其中：
Mi??Y?xixi?xi*7， M0?Y(x)??i?1*?Y?xixi?xi*xi.
*而由Xi~N(?i,?i2)，且Xi间相互独立，又由Y为Xi的线性组合，所以Y也服从正态分布，设Y~N(?',?'2)，则根据概率统计知识：
??'?Y(x*);?7 ?222??'??Mi?i.i?1?记Y的概率密度函数为f(Y),则：
12??'?(Y??')2?'22f(Y)?e.
这样就确定了参数y的统计特征。
5.2 构建目标函数
本问题是一个有约束条件的非线性规划问题。
我们已经知道有108种可能的7个零件参数容差等级组合的情况，如果求出这108种7个零件参数的容差等级组合所对应的108个总费用，比较这108个总费用的大小，取最小值，则对应的7个零件参数的标定值和容差组合即为题目所求。
现讨论7个零件参数的容差等级组合在第k种情况下，生产1000个产品的总费用Zk。 因为：