零件参数设计的数学模型含matlab程序

1零件的参数设计的模型分析摘要本文以产品的成本和产品期望损失之和为目标函数，以标定值和容差为变量建立非线性优化模型。问题中产品的参数由7个零件参数确定。零件参数由标定值和容差共同确定。零件的容差越小，产品参数y偏离0y概率就越小，损失也就越小，但是成本变高。据此分析我们建立一个以损失和成本之和为目标函数的优化模型。我们通过模拟试验猜测题中零件的参数服从正态分布，其i和i分别是标定值和容差的三分之一。并采用2检验法进行了假设检验，确认其服从正态分布。在其基础上对经验公式进行简化处理，得到产品参数y关于零件参数x的线性函数：'y=24.58961x-5.99112x+14.66753x-4.02814x-1.15045x-0.05396x-1.15047x+3.4512y是ix的线性组合，再次证明产品参数服从正态分布，进而可求出损失期望和目标函数。目标函数为：Min=成本+损失期望。通过对题中所给数据的求解得到目标函数值为307.7万元。在对模型的求解中，零件的标定值是连续的，采用迭代法搜索；容差是离散的，采用穷举法搜索。通过matlab求得最优解为：Min=40.12725万元，此时x1=0.075，x2=0.225，x3=0.075，x4=0.075，x5=1.125，x6=18.0974，x7=0.8479。7个零件选取的容差等级依次为BBBCCBB。关键词：零件参数正态分布迭代法穷举法一、问题提出一件产品由多个零件组成，标志产品性能的参数取决于这些零件的参数。每个零件的参数是独立的，零件的参数是标定值和容差。假设每个零件不存在容差，则这件产品的参数是一个定值，但是这个假设不符合实际情况。实际生产过程中，零件的参数总是2出现在一个区间而不是一个点，即实际值总是偏离标定值的。当这些零件组装成产品时，产品的参数就不是一个定值，也将成为一个取值区间。如果产品的参数偏离原先设计值过多（y偏离0y0.3）这个产品就报废，带来9000元的损失；如果偏离的不太多（y偏离0y0.1）这个产品就成为次品，带来1000元的损失。产品的参数偏离设计值的多少是由多个零件参数的容差等级确定的。零件容差等级越高，产品参数的偏离值较小的概率就大，损失的费用也就越小，但是生产零件的成本就会变高。零件容差等级越低，产品参数的偏离值较大的概率就大，损失的费用也就愈大，但是生产零件的成本会降低。当批量生产时就会存在一个最优的零件容差等级组合，使期望成本与期望损失之和达到最小。本文就是要建立一个数学优化模型，来求解这个最优组合。二、基本假设1、假设1：7个零件的参数标定值均服从正态分布，且彼此独立。2、假设2：零件参数的容差为均方差的３倍。3、假设3：零件参数的目标值0y为1.50，当y偏离y00.1时，产品为次品，质量损失为1000元；当y偏离y00.3时，产品为废品，损失为9000元。4、假设4：产品的参数y只由七个零件标定值1x、2x、3x、4x、5x、6x、7x决定。三、符号说明ix零件参数的标定值y产品的参数'y优化设计时产品的参数1C单个产品的成本（元）32C单个产品的损失费用（元）产品参数的期望值产品参数的均方差ijD第i个参数第j容差等级零件的成本，j=1,2,3.假定C等为1，B等为2，A等为3ijR第j等容差与第i个参数的标定值的相对值fy产品参数y的概率密度函数()Fy产品参数y偏离目标值的概率函数1p产品为合格品的概率2p产品为次品的概率3p产品为废品的概率ijT容差等级选择矩阵（用1表示选择该容差，用0表示不选择该矩阵。且每一行有且只有一个值为1）.四、模型的分析建立与求解4.1模型的数据分析4.1.1对产品参数y随机分布的研究我们通过matlab软件对原设计中零件的7个参数分别随机选取了2万个值，将其代入经验公式，得到2万个y的值。我们用excel对y的分布进行了数据分析（如表1），平均1.7160654表一并得到了y值分布的直方图（如图1）随机y值分布直方图05001000150020002500300035001.41.461.521.581.641.71.761.821.881.9422.062.12y值频率频率图1根据直方图，我们不妨猜测y的随机分布函数服从正态分布。=x=1.7160，=S=0.1013。采用分布拟合检验的2检验法，根据如下的定理：标准误差0.010139标准差0.101393方差0.010281最小值1.531051最大值1.958957置信度(95.0%)0.0201195定理：若n充分大（n50）,0H：总体x的分布函数为()Fx，则当0H为真时（无论0H中的分布属何种分布），统计量总是近似地服从自由度为k-r-1的2分布；其中，r是被估计的参数个数。于是，若在假设0H下算得有22(1)kr其中221()kiiiifnpnp则在显著水平下拒绝0H，否则，就接受0H我们将y值分成26个组，用2检验法，代入，的估计值，得：2=12.3654因为220.050.05(1)(2611)36.415kr12.3654所以，将y的分布看成正态分布是可以接受的。4.1.2经验公式的简化由于y关于ix的函数是非线性的，由ix的均值求解y的均值，由ix的均方差求y的均方差都及其复杂，所以我们考虑通过多元函数的泰勒展开来实现，借助于matlab软件，我们得到了简化后的经验公式（程序见附录二），取常数项和一次项之后，结果如下：'y=24.58961x-5.99112x+14.66753x-4.02814x-1.15045x-0.05396x-1.15047x+3.45124.1.3原设计的总费用在原设计中，7个零件参数的标定值分别为：x1=0.1，x2=0.3，x3=0.1，x4=0.1，x5=1.5，x6=16，x7=0.75；容差均取最便宜的等级。而容差通常规定为均方差的3倍，易求得()ix3iijxR（其中ix为第i中零件参数的标定值，iR为第i中零件容差与标定6是的相对值）。根据随机变量线性组合数学期望及标准差的性质：11()(())nniiiiiiECXCEX，211(())(())nniiiiiiDCXCDX（iC为常数）得：()y=1.7256，()y0.1105所以假定y~N（1.7256，0.1105），则22()21()2yfye，()()yFyfydy由此可以得到：合格品的概率：1p=(1.6)(1.4)FF0.1262次品的概率：2p=(1.8)(1.6)(1.4)(1.2)FFFF=0.6234废品的概率：3p=0.2504原设计中容差均采用最便宜的等级，于是单个零件的成本费用为：C1=25+20+20+50+50+10+25=200（元）而单个零件的损失费用=1000次品的概率+9000废品的概率，即：C2=1000*2p+9000*3p=2877（元）总费用：C=1000（C1+C2）=307.7（万元）4.2模型的建立和求解74.2.1模型的建立对于单个产品，其费用由y偏离y0造成的损失和零件成本两部分构成，不妨设定目标函数：min1000(12)CC其中1C表示单个产品的成本，2C表示单个产品的由y偏离y0造成的损失费用。由题意，成本费用由7个零件参数所选择的容差等级决定，零件容差的大小决定了其制造成本，容差设计得越小，成本越高。所以成本费用可以表示为：73111ijijijCDT其中ijD表示第i个参数第j容差等级零件的成本，j=1,2,3（假定C等为1，B等为2，A等为3），根据题意，可以用矩阵的形式表示为：10000251000020501000020502005010050050100001000010251001000025100ijD（注：矩阵中所有等于10000的值都表示没有设定该等级容差，并不代表成本费用）而ijT是关于容差等级的一个选择矩阵，用1表示选择该容差，用0表示不选择该矩阵。且每一行有且只有一个值为1.如：8010100100100100100010ijT表示从x1到x7的容差等级分别为BCCCCCB.对于单个产品的损失费用，有如下设定：当产品参数y偏离y00.1时，产品为次品，质量损失为1,000元；当y偏离y00.3时，产品为废品，损失为9,000元，由于其偏离值的随机性，我们只能采用求数学期望的办法。首先，我们将7个零件参数的标定值代入经验公式，得到y的均值而()ix3iijxR，所以7221()()3iijiixRyw，ijR表示第j等容差与第i个参数的标定值的相对值，iw表示简化后经验公式中ix前面的系数。从而可以求到2'2()'21()2yfye，''''()()yFyfydy由此，合格品的概率：1p=(1.6)(1.4)FF次品的概率：2p=(1.8)(1.6)(1.4)(1.2)FFFF废品的概率：3p=1-(1.8)F+(1.2)F23C210009000pp综上所述，我们建立如下的数学模型：min1000(12)CC91234567310.0750.1250.2250.3750.0750.1250.0750.01251.1251.875.12200.56250.935011ijijjxxxxxstxxTT或4.2.2模型的求解我们通过matlab软件编程，得到目标函数的最优解为min=40.12725万元，比原设计减少了267.57275万元。此时的x1=0.075，x2=0.225，x3=0.075，x4=0.075，x5=1.125，x6=18.0974，x7=0.8479。7个零件选取的容差等级依次为BBBCCBB（程序见附录三）五、模型的评价与推广5.1模型的评价5.1.1对于泰勒展开式的误差分析我们将原设计中7个零件参数，代入原经验公式得到y=1.7256，而使用简化之后的经验公式得到'y=1.72597，相对误差'||*100%yyEy=0.021%，表明简化公式精确度很高。5.1.2模型的优点1.模型将原来复杂的经验公式简化为多元线性函数表达式，且误差极小，大大减少了计算的复杂度102.采用计算机模拟的方法，为前期确定y随x的随机分布情况提供了依据，对结果的验证也起到了很好的补充作用3.模型准确的解决了y偏离y0造成的损失和零件成本的最小值，并提供了7个参数标定值和容差的最佳选择5.1.3模型的缺点1．为了得到7个零件标定值的精确结果，势必将搜索步长减小，加大了算法的时间复杂度，使计算机运算时间更长。5.2模型的推广该模型很好解决了优化零件参数设计对产品参数合格率的提高，可以广泛应用于生产和生活，从而提高生产效益。七、参考文献【1】许伯强等,大学物理实验，江苏大学出版社，2001年2月【2】上海交通大学数学系，概率论与数理统计，科学出版社，2007年2月【3】韩中庚，数学建模方法及其应用，高等教育出版社，2009年6月【4】盛骤，谢式千等，概率论与数理统计，高等教育出版社，1989年八、附录8.1附录清单附录1：对经验公式做20000次模拟的matlab程序附录2：求解y关于ix的泰勒展开式的matlab程序附录3：求解模型的matlab程序8.2附录正文11附录1．symsx1x2x3x4x5x6x7Y1=174.42*(x1/x5)*((x3/(x2-x1))^0.85);Y2=(1-2.62*(1-0.36*((x4/x2)^-0.56))^(3/2)*(x4/x2)^1.16)/(x6*x7);Y=Y1*sqrt(Y2)mm=zeros(1,20000);fori=1:20000g1=normrnd(0