自动化车床模型(数学建模).doc

自动化车床管理模型摘要 本文研究的是自动化车床管理中定期检查和预防性保全刀具问题。在现代技术下，被动地等待故障发生,然后投入较高资金处理出现的问题，这种传统的处理方法已经不符合工业生产和现代社会的发展要求。为解决此问题，我们共建立两个模型，使自动化车床管理方略更科学、更合理。对于问题一：我们通过一定的数学方法，巧妙地建立了生产每一个零件的平均损失费用（包括预防保全费用， 检查费用， 和故障造成的不合格品损失和修复费用，即)关于刀具定期更换间隔的单变量函数关系，并利用MATLAB等数学计算工具和多种方法，对进行逐个赋值,最终得到：当件时，取得最小值元，再根据与固定检查间隔之间的函数关系得到：件。对于问题二:此问中效益函数计算的费用与第一问相比，增加了错误判断带来的损失费用，我们将因误判带来的费用考虑到生产每一个零件的平均损失费用中，用与第一问类似的模型求解，得到当件时，取得最小值元，对应固定检查间隔件。对于问题三：保持问题二的情况，我们建立新的模型，并采取连续检查多个零件（最多3次）的方法,降低误判率，从而达到减少每个零件的平均管理费用，使模型更优化。最终得到在工序发生故障时误判率为0.208，比检查一次的误判率0.4减少0。192，误判率减小了50%；在工序正常时误判率为0。000792，比检查一次的误判率0。02减少0。019208，误判率降低了96.04，从而使模型得到优化。关键词:自动化车床管理 效益函数 正态分布 误判率1。问题重述一道工序用自动化车床连续加工某种零件，由于刀具损坏等原因该工序会出现故障,其中刀具损坏故障占95%, 其它故障仅占5。工序出现故障是完全随机的， 假定在生产任一零件时出现故障的机会均相同。工作人员通过检查零件来确定工序是否出现故障.现积累有100次刀具故障记录，故障出现时该刀具完成的零件数如附表。现计划在刀具加工一定件数后定期更换新刀具.已知生产工序的费用参数如下：故障时产出的零件损失费用 f=200元/件；进行检查的费用 t=10元/次；发现故障进行调节使恢复正常的平均费用 d=3000元/次(包括刀具费)；未发现故障时更换一把新刀具的费用 k=1000元/次。要求解的问题：1)假定工序故障时产出的零件均为不合格品，正常时产出的零件均为合格品， 试对该工序设计效益最好的检查间隔(生产多少零件检查一次)和刀具更换策略。2）如果该工序正常时产出的零件不全是合格品，有2为不合格品；而工序故障时产出的零件有40为合格品,60为不合格品。工序正常而误认有故障停机产生的损失费用为1500元/次。对该工序设计效益最好的检查间隔和刀具更换策略。3）在2)的情况， 可否改进检查方式获得更高的效益。附:100次刀具故障记录（完成的零件数) (见附录1）2。模型假设和符号说明2。1模型的假设1。换刀间隔和检查间隔很短，这段时间内产生的零件数可以忽略；2。检查时一旦发现不合格品生产立即停止；3。假设该工序只需用一个刀具；4。工序出现故障是随机的,且加工每个零件时出现故障的概率相同。2。2符号说明故障时产出的零件损失费用，本题=200元每次进行检查的费用,本题=10元/次发现故障进行调节使恢复正常的平均费用 (包括刀具费)，本题=3000元/次未发现故障时更换一把新刀具的费用，本题=1000元/次生产每个零件的平均管理成本（平均损失费用）零件生产定期检查间隔刀具定期更换间隔工序的平均故障间隔（包括刀具和故障和其他故障）相邻两次检查的后一次检查发现故障时， 件零件中不合格品的平均数平均故障率， 进行预防保全更换道具后,刀故障平均间隔非刀具故障（其他故障）平均间隔工序正常时不合格品率，本题第二问中=2%工序故障时合格品率，本题第二问中=40%工序正常而误认有故障停机产生的损失费用,本题=1500元/次因误判两次故障间隔产生的不合格品均值3. 问题分析本题研究的是自动化车床管理中定期检查和预防性保全刀具问题.由于发现故障进行调节使恢复正常的平均费用和故障时产出的零件损失费用较高，而检查费用和更换刀具费用较低，因此，在未进行定期检查和预防性更换刀具时损失费用明显要高,更新管理后，要确定最优的管理方案，得到合适的检查间隔和刀具更换策略,即刀具定期更换间隔和检查间隔，使损失费用最少。我们选择建立一个效益函数L(生产每个零件的平均损失费用）,通过求L最小值来确定刀具定期更换间隔和生产检查间隔.针对问题一：因为效益函数包括预防保全费用, 检查费用， 和故障造成的不合格品损失和修复费用(即)，我们按每个零件分摊，分别列式子算出、求和.在求解、过程中,根据样本数据确定刀具寿命的经验分布或拟合分布F （x ），并且得到相关参数，确定无预防性更换刀具时刀具故障的平均间隔和采取有预防性更换策略时刀具故障间隔的表达式,同时利用与之间的函数关系进行转换,最后得到关于刀具定期更换间隔的单变量函数,可通过取不同的步长逐个求出的值，找到的最小值和对应的值，并通过所得求出对应的生产检查间隔。针对问题二:第二问中条件变为该工序正常时产出的零件不全是合格品，有2%为不合格品;而工序故障时产出的零件有40为合格品，60%为不合格品且工序正常而误认有故障停机产生的损失费用为1500元/次。所以第二问中计算效益函数时要考虑两种误判情况带来的损失，一是工序正常时检查到不合格品,误判停机将使检查的费用增加；二是工序故障时检查到合格品，将继续生产直到下一次检查，使不合格数量增多造成的损失。效益函数的表达式因此有所改变，但建模方法基本与第一问相似.针对问题三：在第二问情况下只定期检查一个零件造成了许多误判的情况，从而使每个零件的平均管理费用较高。针对这种情况，我们采取适当连续检查多个零件降低误判率的方法达到减少每个零件的平均管理费用，是模型更优化.4。数据分析刀具故障时加工的零件服从正态分布根据所给的100个样本数据用Excel统计可得以下结果：表4.1 样本统计结果数据个数100平均值600标准偏差196。6291695最小值84最大值1153小数位数0区间个数16区间宽度66.875图4。1 样本分布直方图再利用6SQ统计软件进行卡方拟合优度检验得到如下结果：表4.2 刀具故障时加工的零件服从正态分布的检验结果假设检验零假设服从正态分布自由度9卡方统计量2。521839719p值0。980290368显著性水平0.05结果接受零假设5问题一解答5.1 模型一的建立首先，在预防保全措施下，定义平均每一个零件的损失费用(效益函数）为,表达式为： (1） 其中,为预防保全费用, 检查费用， 为故障造成的不合格品损失和修复费用。为工序平均故障间隔,为相邻两次检查的后一次检查发现故障时,件零件中不合格品的平均数。令在相邻两次检查的后一次发现故障的条件下，出现件不合格品的概率为：则得 ： (2）上式中， 上式经运算可得:由此得到与的函数关系,代入(1) 得: （3）然后给出的100 个数据分析算出无预防性更换时, 刀具故障平均间隔为=600 件，再根据题设刀具故障占95 , 非刀具故障占5 ， 故非刀具平均故障间隔为=11400件，其次由100 个数据确定刀具寿命的经验分布或拟合分布得密度函数： 为分布函数。其中,所以当进行预防保全定期更换刀具时, 刀故障的平均间隔为:工序的平均故障间隔由 和 决定, 满足 即： 由此得到和之间的函数关系.5。2模型一的求解对目标函数(3) 的参数进行如下优化: 1.给定， 计算出，代入（3)可得知是 的函数,易得当 时目标函数达到极小； 2. 按一定的步长取的值, 采取逐渐缩小范围并逐个求出 的极小值及相应的值， 其中使最小者所对应的和 即为所求。本步骤我们在MATLAB软件中计算（对应程序段见附录2）得到的结果如下： 元，此时对应的 件，件图5.1 取值与的变化关系图5。3结果分析: 通过逐步求解,我们得到当定期换刀间隔为342件,定期检查间隔为16件时,使得生产每个零件的平均管理成本（平均损失费用）达到极小值5。297元。由取值与的变化关系图我们发现，当步长的取值小于或大于342件时均会使平均损失费用增大。6。问题二的解答6.1确立新目标函数关系令,第二问的效益函数要考虑两种误判。一是工序正常时检查到不合格品误判停机, 将使损失的费用增加；二是工序故障时检查到合格品, 将继续生产直到下一次检查， 使不合格品数量增加，此时两次故障间由此产生的不合格品平均数为：上式中 为工序故障时的合格品率, 为工序在生产一零件时的平均故障率，故在第该问的条件下， 效益函数应为： (4)上式中是工序正常时零件的不合格品率， 元为第一种误判产生停机的损失费。6.2问题二的求解计算再应用与问题一相同的模型代入式（4）并用MATLAB软件求解（对应程序段件附录3）可得以下结果：元，此时对应的 件，件图6.1 取值与的变化关系图63结果分析： 在第二问的条件下运用和第一问相同的模型逐步求解，我们得到当定期换刀间隔为299件，定期检查间隔为18件时，使得生产每个零件的平均管理成本（平均损失费用）达到极小值7.381元。再由取值与的变化关系图我们发现，当步长的取值小于或大于299件时均会使平均损失费用增大.7.问题三的解答7。1模型二的建立由于在第二问情况下生产出的零件不完全是合格品或是不合格品，只定期检查一个零件时,不免会出现许多误判，造成不必要的损失.。今改进检查方式,采取适当连续检查多个零件降低误判率，尽量避免这样的损失，将获得更高的效益.具体方法如下：假设工序发生了故障:第一次检查第二次检查第三次检查结论结论正确性发生概率11未发生故障否0。16101未发生故障否0。096011未发生故障否0.096010发生故障是0。14400发生故障是0.036误判率=0.208(表中“1”表示检测到合格品，“0表示检查到不合格品)假设工序未发生了故障：第一次检查第二次检查第三次检查结论结论正确性发生概率11未发生故障是0.9604101未发生故障是0。019208011未发生故障是0。019208010发生故障否0.00039200发生故障否0.00004误判率=0。000792(表中“1”表示检测到合格品,“0表示检查到不合格品）7.2结论及分析在进行检查时，不能由一次检查结果妄下结论，这可能因误判造成大的损失.由以上模型可以得出：在工序发生故障时误判率为0.208，比检查一次的误判率0.4减少0。192；在工序正常时误判率为0.000792，比检查一次的误判率0。02减少0。019208。因为在工序故障时误判引起的损失为200元/件，工序正常时因误判停机造成的损失为1500/次，这远远大于检查费用10元/次。另外,此工序检查的错误率得到极大降低，对于其他工序的正常运行也有一定共贡献，带来了间接效益。所以采取适当连续多检查几次从而大幅度减小误判率的策略可以获得更高效益。8.模型的评价、改进和推广8.1模型的评价 优点：1.通过建立的模型确定了最优的检查间隔和道具更换间隔，从而将平均每一个零件的损失成本降低；2.优化后的模型二将检查的错误率极大地降低，提高了生产效率和效益;3.由于工序故障率较低,将平均每一个合格零件的损失费用近似等于平均每一个零件的损失费用,使模型得到简化. 缺点：1. 忽略了换刀间隔和检查间隔时间内产生的零件数;2.实际生产中还要考虑流水线上的其他工序的状态与此工序的相互影响，里没有考虑。8.2模型的改进由于流水线上的其他工序与此工序往往是关联的，尤其是故障率和误判率对其他工序有很大影响。因此各工序应看作整体,综合起来考虑,得到的评价参数对工业生产会更有实际指导意义。8.3模型的推广该模型可以广泛应用于汽车、飞机等机械零部件加工、电子仪器生产等生产车间的管理中,一方面有助于降低成本，另一方面提高机器的工作效率.参考文献：1 盛骤,谢式千,潘承毅，概率论与数理统计,北京 :高等教育出版社，2008。6。2 孙山泽，数学的实践与认识，北京：北京大学，2000。1（1）。附录附录1：100次刀具故障记录（完成的零件数）4593626245425095844337488155056124524349826407425657065936809266531644877346084281153593844527552513781474388824538862659775859755649697515628954771609402960885610292837473677358638699634555570844166061062484120447654564339280246687539790581621724531512577496468499544645764558378765666763217715310851附录2：问题一MATLAB求解程序：a=1：1500; l=1000./a+20./sqrt(1140*cxs(a)./(cxs（a)+11400）+31/114+3100./cxs(a）; plot(a，l) axis（0 900 0 30）ylabel(每个零件的平均管理成本/元)xlabel(刀具定期更换间隔/件）function f=cxe(x)f=normcdf（x,600,196。6);function u=cxs(m）k=0;for i=1：m1 k=i（cxe（i）cxe（i-1)+k;endu=(1。/cxe（m)）.*（k+m。*(1-cxe(m)）；最小值的判断：fun=inline(1000./x+20./sqrt（1140cxs（x）./（cxs（x）+11400)）+31/114+3100./cxs（x）)fun = Inline function: fun(x) = 1000./x+20。/sqrt（1140cxs（x）。/(cxs(x)+11400)）+31/114+3100。/cxs（x） x，feval，exitflag，output=fminbnd（fun，0，900)x = 342。0000feval = 5。2972exitflag = 1output = iterations： 31 funcCount: 32 algorithm: golden section search, parabolic interpolation message： 1x112 char附录3：问题二MATLAB求解程序:a=1:1500; l=1000./a+（1./sqrt(1140cxs（a）./（11400+cxs(a)）).*(10+30*(1-1./cxs(a）-1/11400）.sqrt（1140cxs（a)./(11400+cxs(a）)）)+。. （1140+cxs（a)）。/(114cxs(a)）。*（sqrt(1140cxs（a）。/（11400+cxs（a)+1）+。.。 400/3*sqrt(11400+cxs（a)./（114000*cxs(a)）+30*（11400+cxs（a）./(114cxs（a）; plot（a，l) axis(0 900 0 30) ylabel(每个零件的平均管理成本/元）xlabel（刀具定期更换间隔/件）function f=cxe(x） f=normcdf（x，600，196。6)；function u=cxs(m)k=0;for i=1:m1 k=i*（cxe(i）-cxe（i1）+k；endu=（1。/cxe(m）).*（k+m.*（1cxe(m)）;最小值的判断:fun=inline(1000./a+(1。/sqrt（1140*cxs(a）。/（11400+cxs(a）).*（10+30（1-1。/cxs(a）-1/11400）.sqrt（1140cxs(a)./(11400+cxs(a）)+(（1140+cxs（a）)。/(114cxs(a）).（sqrt(1140*cxs(a）./（11400+cxs（a）)+1)+400/3sqrt(（11400+cxs（a）)./（114000*cxs（a）)）+30（11400+cxs(a）./（114*cxs(a)）fun = Inline function: fun（x） = 1000./x+（1。/sqrt(1140cxs（x)./（11400+cxs(x）)）.(10+30(1-1./cxs（x)1/11400）.sqrt（1140*cxs（x)./(11400+cxs（x)）)+(（1140+cxs（x)。/（114*cxs（x）)）.(sqrt(1140cxs（x）。/(11400+cxs(x)）+1）+400/3*sqrt（(11400+cxs(x）。/（114000cxs(x）)+30（11400+cxs(x）./(114*cxs（x)） x，feval,exitflag,output=fminbnd(fun，0，900）x = 299.0001feval = 7。3814exitflag = 1output = iterations： 28 funcCount： 29 algorithm： golden section search, parabolic interpolation message： 1x112 char10