零件参数设计的数学模型.doc

《零件参数设计的数学模型.doc》由会员分享，可在线阅读，更多相关《零件参数设计的数学模型.doc（17页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、. .零件参数设计的数学模型摘要：本文基于Y偏离Y0 造成的损失和零件本钱，根据原设计给定的标定值和容差，使用网格法和随机搜索法,利用计算机编程计算产品分别为正品、次品、废品时的概率，进而分析产品是正品、次品、废品的概率的稳定性，得到较为准确且合理的结果，最后求出原设计的总费用损失费+本钱费为313.4万元。本文通过分析参数x1,x2,x7对y的影响，在原设计的标定值附近找出一个使y在其附近的变化比拟稳定的点，并使y=1.5，再利用计算机仿真实验，综合判断容差等级方案，确定出比拟理想的标定值和容差等级方案：最后确定的方案比原设定节约费用271.2425万元。一、问题的重述一件产品由假设干零件组

2、装而成，标志产品性能的某个参数取决于这些零件的参数。零件参数包括标定值和容差两局部。进展批量生产时，标定值表示一批零件该参数的平均值，容差那么给出了参数偏离其标定值的容许X围。假设将零件参数视为随机变量，那么标定值代表期望值，在生产部门无特殊要求时，容差通常规定为均方差的3倍。在进展零件参数设计时，由于零件组装产品的参数偏离预先设定的目标值，所以造成质量损失，偏离越大，损失越大；且零件的容差大小决定了其制造本钱，容差设计的越小，本钱越高。有一种离子别离器某参数记作Y)由7个零件的参数(记作X1 ,X2 , X7)决定，经历公式为： Y=174.42Y的目标值(记作Y0)为1.50。假设Y偏离Y

3、00.1时，产品为次品，质量损失1000(元)；假设Y偏离Y00.3时，产品为废品，损失9000(元)。零件参数的标定值有一定的容许变化X围；容差分为A、B、C三个等级，用与标定值的相对值来表示，A等为1%，B等为5%，C等为10%。7个零件参数标定值的容许X围及不同容差等级零件的本钱(元)如下表(符号 / 表示无此等级零件):标定值容许X围 C等 B等 A等 X10.075,0.125 / 25 / X20.225,0.375 20 50 / X30.075,0.125 20 50 200 X40.075,0.125 50 100 500 X51.125,1.875 50 / / X612,

4、20 10 25 100 X70.5625,0.935 / 25 100现进展成批生产,每批产量1000个。在原设计中，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符合要求的零件只考虑自身本钱，而不再考虑其它因素的影响。符号约定M 表示成批生产时每批产量的个数，此题为1000个；a 表

5、示产品为次品时的质量损失为1000元；b 表示产品为废品时的质量损失为9000元；表示第i个零件参数对应的均方差；表示一批零件第i个零件参数的平均值，即期望值；表示第i个零件变量的新值；Ri表示变量Xi对的搜索区域；Kd表示区域缩减系数，其值正数；r 表示0，1之间服从均匀分布的伪随机数；k 表示随机概率的分布系数，是个正奇数；z y偏离的绝对值；P y偏离造成的损失；P 表示零件的本钱；Q y偏离造成的损失和零件本钱三、问题的分析由于标志产品性能的参数是由零件的参数所决定的。而零件的参数包括标定值和容差两局部。如果将零件参数视为随机变量，那么标定值代表期望值。那么，根据原理，在其中的概率为：

6、0.9974。显然，在此之外的概率为：0.0026。相比之下，在其之外的可以忽略不计。故此，在生产部门无特殊要求时，容差规定为均方差的3倍是合理的。由题意，我们还可以得到：容差与标定值的相对值可以判断容差的等级进而可以确定零件的本钱，即： A等： B等：0.10.3 C等：进展零件参数设计，就是要确定其标定值和容差。此时要考虑到产品的损失和零件本钱，而产品的损失和零件的本钱都是由零件参数决定。所以，我们就先从产品的零件参数着手，逐步求优。零件参数x1,x2,x7对y的影响由经历公式：来确定，因的目标值记作为1.50。且：当偏离时，产品为次品，质量损失为1,000元；当偏离时，产品为废品，损失为

7、9,000元。可见，选定的标定值x1,x2,x7使得y的值接近1.5，且在x1,x2,x7附近y的取值稳定在1.5附近。所以，我们所设计零件参数，就要尽可能使产品为正品的数量多，次品的数量少、尽量使废品不出现，从而使得总费用损失费+本钱费最小。四、模型的建立在原设计中，组成离子别离器的七个零件参数的标定值为： X1=0.1，X2=0.3，X3=0.1，X4=0.1，X5=1.5，X6=16，X7=0.75 将以上标定值代入公式：得出：显然：大于0.1且小于0.3由y的取值符合正态分布，可以看出在该标定值下，产品出现“次品和“废品的概率较大。由于零件的容差均取最廉价的等级，故此，可得出七个零件参

8、数可能的取值X围如下表： X1 X2 X3 X4 X5 X6 X7取值X围0.095,0.1050.27,0.330.09,0.110.09,0.111.35,1.6514.4,17.60.7125,0.7875为了计算造成的损失和零件本钱。我们给出了两个模型。模型一首先考虑造成的损失，由于给出的零件参数都有一定的容差，所以零件本钱即可确定。进一步，由零件的参数决定的产品参数也在一定的X围之内。而要确定损失，首要问题就是要确定生产一批产品中正品、次品、废品出现的数量。在此之前，我们先对一批产品中正品、次品、废品的概率做一计算；根据条件我们建立了以下的模型：其中为参数Xi标定值容许X围，为容差等

9、级。模型二利用随机搜索法，由于零件的参数是随机的参数且符合正态分布，所以，我们构造出另一模型：其中为随机概率的分布系数，是个正奇数，以保证值可正可负，其值通常取1，3，5，7等，其中K的值越大，那么所构成的函数就越窄，反之越缓。但是在K大于7时在多数情况下，对搜索不很有利，降低了收敛速度。所以，我们在对取值时应尽量避开大于7的数。由正态分布的特点可知：当=1时，显然是不可取的。但是，的取值有规律，即x的取值X围也就是零件的容差越小，就越大，反之越小。五、模型求解及结果分析模型一我们利用网格法亦称枚举法求解，把划定的区域分成假设干个“网格点，然后就各个网格所在的产品规格做一分析，得出正品、次品、

10、废品的概率，从而得出总的费用。于是得出求解方程所以从上式可以看出，求解需进展七次积分，如不利用计算机进展计算，显然很难得出结果，此时我们就编程利用计算机求解。在此，我们利用数学软件编程源程序及求解过程见附录1求解得： P=293.4万元由于零件的等级均取最廉价的，所以，零件的本钱为： P=20万元总的代价为： P总=P+P=293.4+20=313.4万元在此，我们为了使模型具有可靠性，还利用了数学软件在零件参数X围之内随机取值得出结果。当随机循环比拟小时，P总的变化比拟大，即P总的值不稳定，而当随机循环次数比拟大时，P总的值趋向一稳定值。我们把随机循环的次数为20万次与50万次的做一统计：

11、20万次时，P总=313.4万元； 50万次时，P总=314万元。由于在产品中只要出现一个废品，其费用就要增加9000元，而上面得出的结果只相差6000元。所以，可以验证以上得出的结果具有稳定性。模型二我们把模型二结合的数据，对模型中的参数做一分析：把记作零件参数的标定值，零件的容差决定了的取值X围。由于正是用来确定的取值。而是(-1,1)之间的值。所以，我们把记为。我们编程程序参见附录二：程序利用计算机求得P正、P次、P废的概率分别为0.09、0.695、0.215，求得在原设计中y 偏离y0造成的损失和零件本钱共283000元。在编程进展的随机搜索法中，我们发现和d的选择对算法效率有显著的

12、影响。当靠近最优点时，增大和减小d的值，可使P废的概率增大，经过一定次数的迭代，取d=1,K=3.这样我们的模型具有一定的稳定性和合理性。由于我们所建模型时伪随机数r的个数不同，导致在不同次数的计算中，r的值不能一一对应相等。r的个数越多，在我们所编程序中运行次数越多，即步长越小，搜索越细，相对来说计算结果就越准确，所以由于计算时间的限制我们的计算结果免不了会有误差存在。从以上两个模型结果可以看出，计算结果相差无几，这也许是由于随机误差的原因，因为只要在产品中增加一个废品，那么总费用将增加9000元，而两模型的结果相差不到两万元，故此，这点误差是可以容忍的。由于在模型二中，一些参数带有主观色彩

13、，使得计算结果就不能确定其完全可靠，但经过模型一及计算机随机发生器产生的结果检验。而且，当我们计算的循环次数越多，其结果越稳定。故此，模型二还是有一定的可信度。对于模型一，虽然比拟严密，但是计算量特别大，我们设计的程序运行将近两个小时，而模型二只需10分钟就可以得出结果。至于利用数学软件随机发生器计算结果，只是对模型进展验证的一种方法。六、重新设计零件参数由给定的值计算的结果：总费用的期望值313.4万元。可以看出，给定零件参数的标定值，其组成产品某参数在正品的X围之外，且总费用之大，简直不符合实际。对此，我们需重新设计零件参数，使得总费用的期望值降低。所以，我们需对原零件的参数做逐步微调。首

14、先，我们应分析各零件的敏感度零件参数对产品参数的影响程度。先把确定情况下产品参数对零件参数的偏导做一计算。显然，偏导越大，其敏感度就越大。也就是首先应调整的参数。 x1 x2 x3 x4 x5 x6 x7一阶偏导24.5896-5.9910614.6675-4.02809-1.15039-0.053925-1.15039如果对每一个符合条件的值都给予计算，其计算量之大是不可估量的，也是不可能的。故此，我们利用逐步规划，然后上机运行得出标定值比拟好的结果为：即新设计的标定值： x1 x2 x3 x4 x5 x6 x7 0.1 0.3 0.0988 0.1 1.72266 1.6 0.75 当标定

15、值一定的情况下，零件的等级组合有108种，下面我们就将一些组合列出，并计算其总费用值。为了使正品的概率增大来减小质量损失，从而使总的损失减小。首先我们取零件等级较高的情况，得出结果如下表： x1 x2 x3 x4 x5 x6 x7等级 B B A B C C B0.0050.0150.0009880.0050.1722661.60.0375 P正=0.824317 P次=0.164209 P废=0.0000651423 E(费用) =1000(25+50+200+100+50+10+25+10000.164209+90000.0000651423) =624795.278从上面计算的结果可以看

16、出，总的费用比给出的情况下减小了很多，我们为了进一步减小损失，把零件的容差调大，再计算其总费用如下表： x1 x2 x3x4 x5 x6 x7等级 B B B B C C B0.0050.0150.04940.0050.1722661.60.0375 P正=0.808299 P次=0.180162 P废=0.000130192 E(费用)=1000(25+50+50+100+50+10+25+10000.180162+90000.000130192 =491334上面计算结果说明：零件的等级降低后，其总费用显著减小，故此，我们再次把零件容差调大，再观察其总费用： x1 x2 x3 x4 x5

17、x6 x7等级 B B B C C C B0.0050.0150.04940.010.1722661.60.0375 P正=0.813755 P次=0.186095 P废=0.000150114 E(费用)=1000(25+50+50+50+50+10+25+10000.186095+90000.000150114) =447446此时，我们发现费用仍在减小，为了找到总损失最小的情况，继续调大零件容差，计算结果如下： x1 x2 x3 x4 x5 x6 x7等级 B C B C C C B0.0050.030.04940.010.1722661.60.0375此时，计算得到的总费用近似为49万

18、元。很明显，这时的总费用已经增加了，所以，在此之后的情况下，得出的总费用越来越高，故此，其调整方案也就越算越差，在此，我们就不一一列出。且对模型没有帮助。为了进一步寻求较优情况，我们再对上面的情况下，作进一步修改，由于x3的偏导较大即敏感度比拟大，使它的容差减小，再计算其总费用值。 x1 x2 x3 x4 x5 x6 x7等级 B B B C C B B0.0050.0150.04940.010.1722660.80.0375 P次=0.169104 P废=0.0000675851 P正=0.830828 E(费用)=1000(25+50+50+50+50+25+25+10000.169104

19、+90000.000067585) =444713很显然，以上计算即为上面求得的标定值下的最优情况。为了继续降低总费用，我们提出另一种计算零件标定值的方法。即：使各零件参数在标定值处偏导尽量小，且使偏导之和尽量小。再次计算得出的标定值如下： x1 x2 x3 x4 x5 x6 x70.0812460.3746150.1232920.1250001.25043512.001930.935利用上次标定值情况下零件的最优组合，求出此标定值下的最小费用，其计算如下： x1 x2 x3 x4 x5 x6 x7等级 B B B C C B B0.00406230.018730750.00616460.01

20、25000.12504350.60009650.04675 P次=0.146382 P废=0.0000214919 P正=0.853596 E(费用)=(25+50+50+50+50+25+25+10000.146382+90000.0000214919) =421575元经过上述一系列的计算，我们得出了一个比拟满意的结果，把总费用降低到421575元，比原设计的总费用降低了313.4-42.1575=271.2425万元七、模型推广对于任何一位设计工程师来说，总是愿意找出一个最优的设计方案，使所设计的工程设施或产品具有最好的使用性能和最底的材料消耗与制造本钱。而本模型的建立也正是为解决这种问

21、题，所以说本模型具有广泛的普遍性和适用性以及较高的推广价值。就此题来说：粒子别离器某参数由7个零件的参数决定，但在现实生活中其影响参数的因素是不定的，然而无论影响参数的因子有多少，我们都通过模型给出了一个比拟满意的方案。而且在现实的工厂生产中，影响整品的因素N会非常的多，这样如何利用比拟合理的方法解决这样的问题就显得尤为重要，这也正是本模型的意义所在。也就是说只要改变其中的局部系数本模型就可以适用机械化工业部门的生产。另外，这种随机搜索法没有固定的移动模式，而是在可行域内，适应目标函数的下降性质，向最优点作随机移动并靠近它。八、模型评价在本模型中我们首先用网格法，由于所取每个变量有不止一个离散

22、点，借用计算机编程进展计算，假设计算次数较少，那么在很短的时间内就可运行完毕。但无法满足拟合的精度要求假设计算次数较大，也就是说将其进展较细的细化，例如1010，据估计需要将近300小时，那么我们这三天时间是远远不够的。所以这种计算方案是不太合理的。而随机搜索法恰恰防止了网格法的运行时间长的缺点，并且它的合理性较大，总费用也较少。然而它有个缺点是K值较难准确确定，在模型里，我们是用试算法确定的，相对来说也有一定的误差，但误差较小，在这里可忽略不记。九、模型改良对于第一个模型，我们是非常易于理解的，它本质就是要计算满足要求的点落在正品、次品和废品的概率，从而确定费用的最小值。但是这种思想实现却非

23、常的麻烦，因为对于题目所给定的数据，我们要解决的是一个七维的函数，我们首先要将其细化，将其分成空间的个小的立方体，近似的依它中点落在的某个区间的概率来确定出现正品、次品和废品的概率，计算过程中我们要计算的是一个七重的积分，即使利用计算机编程也要花费大量的时间和精力。而对于模型二，却可以完全防止这种情况，因为我们构造的伪随机函数是非常简单的，即使通过改变步长，也不会带来太大的计算麻烦，虽然这个函数是我们随机构造的伪的随机函数，但通过我们的计算，发现它的计算结果完全呈正态分布，而且通过计算所得的结果与用原始方法计算的结果也大致一样，同时还可以比拟合理的检验某一组解的合理性。这充分说明这种算法是非常

24、合理的，但是由于它中间伴随有人为的模糊控制的因素，使它不可能十分的准确，所以我们认为可以通过第一个模型来确定解的大致位置以尽可能的缩小解的X围，再用第一个模型解出一个比拟准确的解，并代回第二个模型检验。这样可以到达减化计算的目的。十、模型优缺点模型二随机搜索法便于实现程序和实际使用，在较大的X围内模拟随机进展，构造简单，但准确性不高。模型一网格法通俗易懂，而且准确度高，但由于费时费工所以推广价值不大。参考文献 1、钢铁学院机械优化设计方法冶金工业。 2、中国数学会主办数学的实践与认识。 3、.周概容概率论与数理统计高等教育。 4、X惟信孟嗣宗机械最优化设计清化大学附录一.概率的计算程序清单输入

25、参数说明: dui-xi的标定值; dfi-xi对应的均方差; ni -区间的分割数; (i=1,2,.,7)源程序:s7du1_,df1_,n1_,du2_,df2_,n2_,du3_,df3_,n3_,du4_,df4_,n4_, du5_,df5_,n5_,du6_,df6_,n6_,du7_,df7_,n7_:=Module sc=0.0,xx1,s1=6.0*df1/n1,xx2,s2=6.0*df2/n2,xx3,s3=6.0*df3/n3, sf=0.0,xx4,s4=6.0*df4/n4,xx5,s5=6.0*df5/n5,xx6,s6=6.0*df6/n6, sz=0.0,i

26、1=0,i2=0,i3=0,xx7,s7=6.0*df7/n7,ddf,xc2,xc3,xc4,xc5, xc6,xc7,fyd, ddf=(df1*df2*df3*df4*df5*df6*df7); fyx1_,x2_,x3_,x4_,x5_,x6_,x7_:= 174.42*(x1/x5)*(x3/(x2-x1)0.85*Sqrt (1-2.62(1-0.36(x4/x2)(-0.56)1.5* (x4/x2)1.16)/(x6*x7); fsx1_,x2_,x3_,x4_,x5_,x6_,x7_:= 0.0016083*s1*s2*s3*s4*s5*s6*s7*Exp(-0.5)*( (

27、x1-du1)/df1)2+(x2-du2)/df2)2+(x3-du3)/df3)2+ (x4-du4)/df4)2+(x5-du5)/df5)2+(x6-du6)/df6)2+ (x7-du7)/df7)2)/df; xx1=du1-3*df1+s1/2; xc2=du2-3*df2+s2/2; xc3=du3-3*df3+s3/2; xc4=du4-3*df4+s4/2; xc5=du5-3*df5+s5/2; xc6=du6-3*df6+s6/2; xc7=du7-3*df7+s7/2; WhileAbsxx1-du13*df1, xx2=xc2; WhileAbsxx2-du23*d

28、f2, xx3=xc3; WhileAbsxx3-du33*df3, xx4=xc4; WhileAbsxx4-du43*df4, xx5=xc5; WhileAbsxx5-du53*df5, xx6=xc6; WhileAbsxx6-du63*df6, xx7=xc7; WhileAbsxx7-du73*df7, fyd=Absfyxx1,xx2,xx3,xx4, xx5,xx6,xx7-1.5; If0.1fyd=0.3, sf=sf+fsxx1,xx2,xx3,xx4, xx5,xx6,xx7; Iffyd=0.1, sz=sz+fsxx1,xx2,xx3,xx4, xx5,xx6,xx

29、7; xx7=xx7+s7; xx6=xx6+s6; ; xx5=xx5+s5; ; xx4=xx4+s4; ; xx3=xx3+s3; i3=i3+1; Printi1, ,i2, ,i3 ; xx2=xx2+s2; i2=i2+1; xx1=xx1+s1; i1=i1+1; ; PrintSc= ,sc, ,Sf= ,sf, Sz= ,sz 运行结果:Mathematica 2.0 for MS-DOS 386/7Copyright 1988-91 Wolfram Research, Inc.In1:= s7.mIn2:= s70.1,0.05*0.1/3,5,0.3,0.05*0.3/3

30、,5,0.0998,0.01*0.0998/3,5, 0.1,0.05*0.1/3,5,1.72266,0.1*1.72266/3,5,16,0.1*16/3,20, 0.75,0.05*0.75/3,50 0 10 0 20 0 30 0 40 0 50 1 64 23 1204 24 1214 24 1224 24 1234 24 1244 24 125Sc= 0.164209 Sf= 0.0000651423 Sz= 0.824317In3:= 1000(310+0.164209*1000+0.0000651423*9000)Out3= 474795.附录二s1 = 0: s2 = 0

31、READ x11, d1, x22, d2, x33, d3, x44, d4, x55, d5, x66, d6, x77, d7 FOR r = 0 TO 1 STEP .01 FOR i = 1 TO 2 FOR j = 1 TO 2 FOR k = 1 TO 2 FOR l = 1 TO 2 FOR m = 1 TO 2 FOR n = 1 TO 2 FOR p = 1 TO 2 x1 = x11 + d1 * (-1) i * (2 * r - 1) 3 x2 = x22 + d2 * (-1) j * (2 * r - 1) 3 x3 = x33 + d3 * (-1) k * (

32、2 * r - 1) 3 x4 = x44 + d4 * (-1) l * (2 * r - 1) 3 x5 = x55 + d5 * (-1) m * (2 * r - 1) 3 x6 = x66 + d6 * (-1) n * (2 * r - 1) 3 x7 = x77 + d7 * (-1) p * (2 * r - 1) 3 NEXT p: NEXT n: NEXT m: NEXT l: NEXT k: NEXT j: NEXT i DATA 0.1,0.005,0.3,0.03,0.1,0.01,0.1,0.01,1.5,0.15,16,1.6,0.75,0.375 y = 174.42 * (x1 / x5) * (x3 / (x2 - x1) .85 * SQR(1 - 2.62 * (1 - .36 * (x2 /4) .56) 1.5 * (x4 / x2) 1.16) / (x6 * x7) s = ABS(y - 1.5) IF s .3 THEN IF s .1 THEN s1 = s1 + 1 ELSE s2 = s2 + 1 PRINT s1,s2=; s1, s2 NEXT r END. .word.