2022年自动化车床模型 .pdf

《2022年自动化车床模型 .pdf》由会员分享，可在线阅读，更多相关《2022年自动化车床模型 .pdf（15页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、1 自动化车床管理模型摘要本文研究的是自动化车床管理中定期检查和预防性保全刀具问题。在现代技术下，被动地等待故障发生， 然后投入较高资金处理出现的问题，这种传统的处理方法已经不符合工业生产和现代社会的发展要求。为解决此问题，我们共建立两个模型，使自动化车床管理方略更科学、更合理。对于问题一： 我们通过一定的数学方法， 巧妙地建立了生产每一个零件的平均损失费用L（包括预防保全费用1L, 检查费用2L, 和故障造成的不合格品损失和修复费用3L，即123LLLL）关于 刀具定 期更 换间 隔a的单变量 函数 关系 ，并利用MATLAB等数学计算工具和多种方法，对a进行逐个赋值，最终得到：当342a件

2、时，L取得最小值min5.297L元，再根据a与固定检查间隔n之间的函数关系得到：16n件。对于问题二： 此问中效益函数计算的费用与第一问相比，增加了错误判断带来的损失费用，我们将因误判带来的费用考虑到生产每一个零件的平均损失费用L中，用与第一问类似的模型求解， 得到当299a件时，L取得最小值min7.381L元，对应固定检查间隔18n件。对于问题三：保持问题二的情况， 我们建立新的模型， 并采取连续检查多个零件 （最多 3 次）的方法，降低误判率， 从而达到减少每个零件的平均管理费用，使模型更优化。最终得到在工序发生故障时误判率为0.208 ，比检查一次的误判率0.4 减少 0.192，误

3、判率减小了50% ；在工序正常时误判率为0.000792，比检查一次的误判率0.02 减少0.019208，误判率降低了 96.04%，从而使模型得到优化。关键词： 自动化车床管理效益函数正态分布误判率1. 问题重述名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 15 页 - - - - - - - - - 2 一道工序用自动化车床连续加工某种零件，由于刀具损坏等原因该工序会出现故障，其中刀具损坏故障占95%, 其它故障仅占 5% 。工序出现故障是完全随机的, 假定在生产

4、任一零件时出现故障的机会均相同。工作人员通过检查零件来确定工序是否出现故障。现积累有 100 次刀具故障记录，故障出现时该刀具完成的零件数如附表。现计划在刀具加工一定件数后定期更换新刀具。已知生产工序的费用参数如下：故障时产出的零件损失费用 f=200 元/ 件；进行检查的费用 t=10 元/ 次；发现故障进行调节使恢复正常的平均费用 d=3000 元/ 次( 包括刀具费 ) ；未发现故障时更换一把新刀具的费用 k=1000 元/ 次。要求解的问题：1)假定工序故障时产出的零件均为不合格品，正常时产出的零件均为合格品, 试对该工序设计效益最好的检查间隔（生产多少零件检查一次）和刀具更换策略。2

5、)如果该工序正常时产出的零件不全是合格品，有 2% 为不合格品； 而工序故障时产出的零件有 40% 为合格品，60% 为不合格品。 工序正常而误认有故障停机产生的损失费用为 1500 元/ 次。对该工序设计效益最好的检查间隔和刀具更换策略。3）在 2) 的情况 , 可否改进检查方式获得更高的效益。附:100 次刀具故障记录 (完成的零件数 ) （见附录 1）2. 模型假设和符号说明2.1 模型的假设1. 换刀间隔和检查间隔很短，这段时间内产生的零件数可以忽略；2. 检查时一旦发现不合格品生产立即停止；3. 假设该工序只需用一个刀具；4. 工序出现故障是随机的，且加工每个零件时出现故障的概率相同

6、。2.2 符号说明名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 15 页 - - - - - - - - - 3 f故障时产出的零件损失费用，本题f=200元t每次进行检查的费用，本题t=10 元/ 次d发现故障进行调节使恢复正常的平均费用 ( 包括刀具费 ) ，本题d=3000元/ 次k未发现故障时更换一把新刀具的费用，本题k=1000元/ 次L生产每个零件的平均管理成本（平均损失费用）n零件生产定期检查间隔a刀具定期更换间隔b工序的平均故障间隔（包括刀具和故障和其他

7、故障）c相邻两次检查的后一次检查发现故障时, n件零件中不合格品的平均数p平均故障率，1pba进行预防保全更换道具后，刀故障平均间隔m非刀具故障（其他故障）平均间隔p工序正常时不合格品率，本题第二问中p=2% q工序故障时合格品率，本题第二问中q=40% e工序正常而误认有故障停机产生的损失费用，本题e=1500元/ 次2因误判两次故障间隔产生的不合格品均值3. 问题分析本题研究的是自动化车床管理中定期检查和预防性保全刀具问题。由于发现故障进行调节使恢复正常的平均费用和故障时产出的零件损失费用较高，而检查费用和更换刀具费用较低，因此，在未进行定期检查和预防性更换刀具时损失费用明显要高，更新管名

8、师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 15 页 - - - - - - - - - 4 理后，要确定最优的管理方案，得到合适的检查间隔和刀具更换策略，即刀具定期更换间隔和检查间隔，使损失费用最少。我们选择建立一个效益函数L（生产每个零件的平均损失费用），通过求L 最小值来确定刀具定期更换间隔a和生产检查间隔n。针对问题一：因为效益函数L包括预防保全费用1L, 检查费用2L, 和故障造成的不合格品损失和修复费用3L（即123LLLL） ，我们按每个零件分摊，分别列式

9、子算出1L、2L、3L求和。在求解1L、2L、3L过程中，根据样本数据确定刀具寿命的经验分布或拟合分布F ( x )，并且得到相关参数，确定无预防性更换刀具时刀具故障的平均间隔和采取有预防性更换策略时刀具故障间隔的表达式，同时利用n与a之间的函数关系进行转换，最后得到L关于刀具定期更换间隔a的单变量函数，可通过取不同的步长a逐个求出L的值，找到L的最小值和对应的a值，并通过所得a求出对应的生产检查间隔n。针对问题二： 第二问中条件变为该工序正常时产出的零件不全是合格品，有2% 为不合格品；而工序故障时产出的零件有40% 为合格品， 60% 为不合格品且工序正常而误认有故障停机产生的损失费用为1

10、500元/ 次。所以第二问中计算效益函数L时要考虑两种误判情况带来的损失，一是工序正常时检查到不合格品，误判停机将使检查的费用增加；二是工序故障时检查到合格品，将继续生产直到下一次检查，使不合格数量增多造成的损失。效益函数L的表达式因此有所改变，但建模方法基本与第一问相似。针对问题三：在第二问情况下只定期检查一个零件造成了许多误判的情况，从而使每个零件的平均管理费用较高。针对这种情况， 我们采取适当连续检查多个零件降低误判率的方法达到减少每个零件的平均管理费用，是模型更优化。4.数据分析刀具故障时加工的零件服从正态分布根据所给的 100个样本数据用 Excel 统计可得以下结果：名师资料总结

11、- - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 15 页 - - - - - - - - - 5 表4.1 样本统计结果数据个数100 平均值600 标准偏差196.6291695 最小值84 最大值1153 小数位数0 区间个数16 区间宽度66.875 22236101416149971311024681012141618117183.9250.8317.7384.6451.5518.4585.3652.2719.1786852.9919.8986.7105411206SQ.Net

12、00.020.040.060.080.10.120.140.16频率（ % ）正态概率图4.1 样本分布直方图再利用 6SQ 统计软件进行卡方拟合优度检验得到如下结果：表4.2 刀具故障时加工的零件服从正态分布的检验结果假设检验零假设服从正态分布自由度9 卡方统计量2.521839719 p 值0.980290368 显著性水平0.05 结果接受零假设5问题一解答名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 15 页 - - - - - - - - - 6 5.1 模型

13、一的建立首先，在预防保全措施下，定义平均每一个零件的损失费用（效益函数）为L，表达式为：123fcdktLLLLanb（1）其中，1L为预防保全费用 , 2L检查费用 , 3L为故障造成的不合格品损失和修复费用。b为工序平均故障间隔，c为相邻两次检查的后一次检查发现故障时，n件零件中不合格品的平均数。令在相邻两次检查的后一次发现故障的条件下，出现i件不合格品的概率为：1(1)1(1)ninpppp则得 ：111(1)1(1)nniiniiinpppcp（2）上式中1, 2, 3in，上式经运算可得：21112()2122nnnppco由此得到c与n的函数关系，代入 (1) 得：(1)2nfkt

14、dLanb（3）然后给出的 100 个数据分析算出无预防性更换时, 刀具故障平均间隔为=600 件，再根据题设刀具故障占95% , 非刀具故障占 5% , 故非刀具平均故障间隔为95%5%m=11400件，其次由 100 个数据确定刀具寿命的经验分布或拟合分布得密度函数：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 15 页 - - - - - - - - - 7 22()21( ),0,26xfxxxZe22()2001( )( )2xxxFxft dtdte为分布函

15、数。其中，600，196.629所以当进行预防保全定期n更换刀具时 , 刀故障的平均间隔为：111( )(1)1()( )aaiiFiF iaFaFa工序的平均故障间隔b由a和m决定, 满足111abm即：111abm由此得到b和a之间的函数关系。5.2 模型一的求解对目标函数 (3) 的参数进行如下优化： 1. 给定a，计算出b, 代入(3) 可得知L是n的函数 , 易得当2btnf时目标函数L达到极小； 2. 按一定的步长取a的值, 采取逐渐缩小范围并逐个求出L的极小值及相应的n值, 其中使L最小者所对应的a和n即为所求。本步骤我们在MATLAB软件中计算（对应程序段见附录 2）得到的结果

16、如下：min5.297L元，此时对应的342a件，16n件名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 15 页 - - - - - - - - - 8 图 5.1 a取值与L的变化关系图5.3 结果分析：通过逐步求解，我们得到当定期换刀间隔为342 件，定期检查间隔为16 件时，使得生产每个零件的平均管理成本（平均损失费用）L达到极小值 5.297 元。由a取值与L的变化关系图我们发现，当步长a的取值小于或大于342 件时均会使平均损失费用L增大。6.问题二的解答6.

17、1 确立新目标函数关系令asn，第二问的效益函数要考虑两种误判。一是工序正常时检查到不合格品误判停机, 将使损失的费用增加； 二是工序故障时检查到合格品, 将继续生产直到下一次检查, 使不合格品数量增加 , 此时两次故障间由此产生的不合格品平均数为：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 15 页 - - - - - - - - - 9 211(1)(1)(1)1(1)1(1)12112jsjijnjniianppqqqunpqnnq上式中40%q为工序故障时的合

18、格品率 , p为工序在生产一零件时的平均故障率 , 故在第该问的条件下 , 效益函数L应为：2(1)1nfkdppteLanbb（4）上式中2%p是工序正常时零件的不合格品率,1500e元为第一种误判产生停机的损失费。6.2 问题二的求解计算再应用与问题一相同的模型代入式（4）并用 MATLAB 软件求解（对应程序段件附录3）可得以下结果：min7.381L元，此时对应的299a件，18n件名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 15 页 - - - - - -

19、- - - 10 图 6.1 a取值与L的变化关系图63 结果分析：在第二问的条件下运用和第一问相同的模型逐步求解，我们得到当定期换刀间隔为299 件，定期检查间隔为18件时，使得生产每个零件的平均管理成本（平均损失费用）L达到极小值 7.381 元。再由a取值与L的变化关系图我们发现， 当步长a的取值小于或大于 299 件时均会使平均损失费用L增大。7.问题三的解答7.1 模型二的建立由于在第二问情况下生产出的零件不完全是合格品或是不合格品, 只定期检查一个零件时 , 不免会出现许多误判 , 造成不必要的损失 . 。今改进检查方式 , 采取适当连续检查多个零件降低误判率，尽量避免这样的损失,

20、 将获得更高的效益。具体方法如下：假设工序发生了故障：第一次检查第二次检查第三次检查结论结论正确性发生概率1 1 未发生故障否0.16 1 0 1 未发生故障否0.096 0 1 1 未发生故障否0.096 0 1 0 发生故障是0.144 0 0 发生故障是0.036 误判率 =0.208（表中“ 1”表示检测到合格品，“ 0”表示检查到不合格品）假设工序未发生了故障：第一次检查第二次检查第三次检查结论结论正确性发生概率1 1 未发生故障是0.9604 1 0 1 未发生故障是0.019208 0 1 1 未发生故障是0.019208 0 1 0 发生故障否0.000392 0 0 发生故障

21、否0.00004 误判率 =0.000792名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 15 页 - - - - - - - - - 11 （表中“ 1”表示检测到合格品，“ 0”表示检查到不合格品）7.2 结论及分析在进行检查时，不能由一次检查结果妄下结论，这可能因误判造成大的损失。由以上模型可以得出：在工序发生故障时误判率为0.208 ，比检查一次的误判率0.4 减少0.192 ；在工序正常时误判率为0.000792 ，比检查一次的误判率0.02 减少 0.01

22、9208。因为在工序故障时误判引起的损失为200 元/ 件，工序正常时因误判停机造成的损失为1500/ 次，这远远大于检查费用10 元/ 次。另外，此工序检查的错误率得到极大降低，对于其他工序的正常运行也有一定共贡献，带来了间接效益。 所以采取适当连续多检查几次从而大幅度减小误判率的策略可以获得更高效益。8. 模型的评价、改进和推广8.1 模型的评价优点：1. 通过建立的模型确定了最优的检查间隔和道具更换间隔，从而将平均每一个零件的损失成本降低；2. 优化后的模型二将检查的错误率极大地降低，提高了生产效率和效益；3. 由于工序故障率较低，将平均每一个合格零件的损失费用近似等于平均每一个零件的损

23、失费用，使模型得到简化。缺点：1. 忽略了换刀间隔和检查间隔时间内产生的零件数；2. 实际生产中还要考虑流水线上的其他工序的状态与此工序的相互影响，里没有考虑。8.2 模型的改进由于流水线上的其他工序与此工序往往是关联的，尤其是故障率和误判率对其他工序有很大影响。因此各工序应看作整体，综合起来考虑，得到的评价参数对工业生产会更有实际指导意义。8.3 模型的推广名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 15 页 - - - - - - - - - 12 该模型可以广

24、泛应用于汽车、飞机等机械零部件加工、电子仪器生产等生产车间的管理中，一方面有助于降低成本，另一方面提高机器的工作效率。参考文献：1盛骤，谢式千 ,潘承毅，概率论与数理统计,北京：高等教育出版社，2008.6。2 孙山泽，数学的实践与认识，北京：北京大学，2000.1（1） 。附录附录 1：100 次刀具故障记录 (完成的零件数 ) 459 362 624 542 509 584 433 748 815 505 612 452 434 982 640 742 565 706 593 680 926 653 164 487 734 608 428 1153 593 844 527 552 513

25、781 474 388 824 538 862 659 775 859 755 649 697 515 628 954 771 609 402 960 885 610 292 837 473 677 358 638 699 634 555 570 84 416 606 1062 484 120 447 654 564 339 280 246 687 539 790 581 621 724 531 512 577 496 468 499 544 645 764 558 378 765 666 763 217 715 310 851 附录 2：问题一 MATLAB 求解程序：a=1:1500; l

26、=1000./a+20./sqrt(1140*cxs(a)./(cxs(a)+11400)+31/114+3100./cxs(a); plot(a,l) axis(0 900 0 30) ylabel( 每个零件的平均管理成本/ 元 )名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 15 页 - - - - - - - - - 13 xlabel( 刀具定期更换间隔/ 件 )function f=cxe(x)f=normcdf(x,600,196.6);function

27、 u=cxs(m)k=0;for i=1:m-1 k=i*(cxe(i)-cxe(i-1)+k;endu=(1./cxe(m).*(k+m.*(1-cxe(m); 最小值的判断：fun=inline(1000./x+20./sqrt(1140*cxs(x)./(cxs(x)+11400)+31/114+3100./cxs(x) fun = Inline function: fun(x) = 1000./x+20./sqrt(1140*cxs(x)./(cxs(x)+11400)+31/114+3100./cxs(x) x,feval,exitflag,output=fminbnd(fun,0,

28、900) x = 342.0000 feval = 5.2972 exitflag = 1 output = iterations: 31 funcCount: 32 algorithm: golden section search, parabolic interpolation名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 13 页，共 15 页 - - - - - - - - - 14 message: 1x112 char 附录 3：问题二 MATLAB 求解程序：a=1:1

29、500;l=1000./a+(1./sqrt(1140*cxs(a)./(11400+cxs(a).*(10+30*(1-1./cxs(a)-1/11400).sqrt(1140*cxs(a)./(11400+cxs(a)+. (1140+cxs(a)./(114*cxs(a).*(sqrt(1140*cxs(a)./(11400+cxs(a)+1)+.400/3*sqrt(11400+cxs(a)./(114000*cxs(a)+30*(11400+cxs(a)./(114*cxs(a); plot(a,l) axis(0 900 0 30) ylabel( 每个零件的平均管理成本/ 元 )

30、xlabel( 刀具定期更换间隔/ 件 )function f=cxe(x) f=normcdf(x,600,196.6);function u=cxs(m)k=0;for i=1:m-1 k=i*(cxe(i)-cxe(i-1)+k;endu=(1./cxe(m).*(k+m.*(1-cxe(m); 最小值的判断：fun=inline(1000./a+(1./sqrt(1140*cxs(a)./(11400+cxs(a).*(10+30*(1-1./cxs(a)-1/11400).sqrt(1140*cxs(a)./(11400+cxs(a)+(1140+cxs(a)./(114*cxs(a

31、).*(sqrt(1140*cxs(a)./(11400+cxs(a)+1)+400/3*sqrt(11400+cxs(a)./(114000*cxs(a)+30*(11400+cxs(a)./(114*cxs(a) fun = Inline function: fun(x) = 1000./x+(1./sqrt(1140*cxs(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(1140*cxs(x)./(11400+

32、cxs(x)+1)+400/3*sqrt(11400+cxs(名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 14 页，共 15 页 - - - - - - - - - 15 x)./(114000*cxs(x)+30*(11400+cxs(x)./(114*cxs(x) x,feval,exitflag,output=fminbnd(fun,0,900) x = 299.0001 feval = 7.3814 exitflag = 1 output = iterations: 28 funcCount: 29 algorithm: golden section search, parabolic interpolationmessage: 1x112 char 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 15 页，共 15 页 - - - - - - - - -