自动化车床模型

《自动化车床模型》由会员分享，可在线阅读，更多相关《自动化车床模型（16页珍藏版）》请在装配图网上搜索。

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

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

3、检查一次的误判率 0.4 减少 0.192，误 判率减小了 50%；在工序正常时误判率为 0.000792，比检查一次的误判率 0.02 减少 0.019208，误判率降低了 96.04%，从而使模型得到优化。关键词：自动化车床管理 效益函数 正态分布 误判率1.问题重述一道工序用自动化车床连续加工某种零件，由于刀具损坏等原因该工序会出现故 障，其中刀具损坏故障占 95%, 其它故障仅占 5%。工序出现故障是完全随机的, 假定在 生产任一零件时出现故障的机会均相同。工作人员通过检查零件来确定工序是否出现故 障。现积累有100 次刀具故障记录，故障出现时该刀具完成的零件数如附表。现计划在 刀具加

4、工一定件数后定期更换新刀具。已知生产工序的费用参数如下：故障时产出的零件损失费用 f=200 元/件；进行检查的费用 t=10 元/次； 发现故障进行调节使恢复正常的平均费用 d=3000 元/次（包括刀具费）； 未发现故障时更换一把新刀具的费用 k=1000 元/次。要求解的问题：1）假定工序故障时产出的零件均为不合格品，正常时产出的零件均为合格品, 试对 该工序设计效益最好的检查间隔（生产多少零件检查一次）和刀具更换策略。2）如果该工序正常时产出的零件不全是合格品，有 2%为不合格品；而工序故障时产 出的零件有 40%为合格品，60%为不合格品。工序正常而误认有故障停机产生的损 失费用为

5、1500 元/次。对该工序设计效益最好的检查间隔和刀具更换策略。3）在 2）的情况, 可否改进检查方式获得更高的效益。附:100 次刀具故障记录（完成的零件数） （见附录 1）2.模型假设和符号说明2.1 模型的假设1. 换刀间隔和检查间隔很短，这段时间内产生的零件数可以忽略；2. 检查时一旦发现不合格品生产立即停止；3. 假设该工序只需用一个刀具；4. 工序出现故障是随机的，且加工每个零件时出现故障的概率相同。2.2 符号说明f故障时产出的零件损失费用，本题f =200元t每次进行检查的费用，本题t =10元/次d发现故障进行调节使恢复正常的平均费用（包括刀具费），本题d =3000元/次k

6、未发现故障时更换一把新刀具的费用，本题k =1000元/次L生产每个零件的平均管理成本（平均损失费用）n零件生产定期检查间隔a刀具定期更换间隔b工序的平均故障间隔（包括刀具和故障和其他故障）c相邻两次检查的后一次检查发现故障时，n件零件中不合格品的平均数P平均故障率，P - 1 /b卩a进行预防保全更换道具后，刀故障平均间隔m非刀具故障（其他故障）平均间隔f p工序正常时不合格品率，本题第二问中P =2%q工序故障时合格品率，本题第二问中q =40%e工序正常而误认有故障停机产生的损失费用，本题e =1500元/次卩2因误判两次故障间隔产生的不合格品均值3. 问题分析本题研究的是自动化车床管理

7、中定期检查和预防性保全刀具问题。由于发现故障进 行调节使恢复正常的平均费用和故障时产出的零件损失费用较高，而检查费用和更换刀 具费用较低，因此，在未进行定期检查和预防性更换刀具时损失费用明显要高，更新管 理后，要确定最优的管理方案，得到合适的检查间隔和刀具更换策略，即刀具定期更换 间隔和检查间隔，使损失费用最少。我们选择建立一个效益函数L (生产每个零件的平 均损失费用)，通过求L最小值来确定刀具定期更换间隔a和生产检查间隔n。针对问题一：因为效益函数L包括预防保全费用L ,检查费用L ,和故障造成12的不合格品损失和修复费用L (即L = L + L + L)，我们按每个零件分摊，分别列3

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

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

10、以下结果：表4.1 样本统计结果0.120.060.080.04003 U5Q8 Mwaz9 C3-数据个数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 模型一的建立 首先，在预防保全措施下，定义平均每一个零件的损失费用(效益函数)为 L ， 表达式为：k

11、t c -f + dL = L + L + L = + +(1)1 2 3 a n b其中， L 为预防保全费用, L 检查费用, L 为故障造成的不合格品损失和修1 2 3复费用。b为工序平均故障间隔，C为相邻两次检查的后一次检查发现故障时，n件 零件中不合格品的平均数。令在相邻两次检查的后一次发现故障的条件下，出现件不合格品的概率为:(1 - P )n - 2 - PP =-11 - (1 - p )n则得 ：(1 - p )n-i p2)1 - (1 - p )n上式中 i = 1，2, 3 n上式经运算可得：2-112由此得到c与n的函数关系，代入(1)得:kL =-at (n + 1

12、) f3)然后给出的100个数据分析算出无预防性更换时，刀具故障平均间隔为卩=600 件，再根据题设刀具故障占95% , 非刀具故障占5% , 故非刀具平均故障间隔为 m = R 95% =11400件，其次由100个数据确定刀具寿命的经验分布或拟合分布得密度 5%函数：(x |L)2 e _2b 2(x _ 口)22b2 dt为分布函数。=F(x) = x f (t)dt = Je2兀b o其中，卩=600，b = 196.629所以当进行预防保全定期n更换刀具时, 刀故障的平均间隔为：F(a) i i F (i) F (i 1)+ a 1 F (a) I i=1工序的平均故障间隔b由卩和

13、m 决定 ,满足1b1+ma即：11 1+卩ma由此得到b和a之间的函数关系。5.2模型一的求解对目标函数(3) 的参数进行如下优化：1给定a，计算出b，代入(3)可得知L是的函数，易得当n -2 bt时目标函数L达到极小;2.按一定的步长取a的值，采取逐渐缩小范围并逐个求出L的极小值及相应 的n值，其中使L最小者所对应的a和n 即为所求。本步骤我们在MATLAB软件中计算(对应程序段见附录2)得到的结果如下：min L = 5.297 元，此时对应的 a = 342 件， n = 16 件3025205O100200300400500600700 SOO 900刀具定期更换间隔/件图5.1

14、a取值与L的变化关系图5.3 结果分析：通过逐步求解，我们得到当定期换刀间隔为 342 件，定期检查间隔为 16 件时，使 得生产每个零件的平均管理成本（平均损失费用）L达到极小值5.297元。由a取值 与L的变化关系图我们发现，当步长a的取值小于或大于342件时均会使平均损失费 用 L 增大。6.问题二的解答6.1 确立新目标函数关系令 a sn ，第二问的效益函数要考虑两种误判。一是工序正常时检查到不合格品误判停机, 将 使损失的费用增加；二是工序故障时检查到合格品, 将继续生产直到下一次检查, 使不 合格品数量增加,此时两次故障间由此产生的不合格品平均数为：-n卫Ti qT（iq）J_j

15、 二 1 i 二 1）（i n - （1 - p j1 - （1 - P ）a上式中q = 40%为工序故障时的合格品率，P为工序在生产一零件时的平均n +1 q 丁 + n 1-q故障率， 故在第该问的条件下， 效益函数 L 应为：2+ b4）k 1，L =+ t +（1 - p）n p e a n上式中P = 2%是工序正常时零件的不合格品率，e = 1500元为第一种误判产生停机的损失费。6.2 问题二的求解计算再应用与问题一相同的模型代入式（4）并用MATLAB软件求解（对应程序段件附录3）可得以下结果：min L = 7.381元，此时对应的a二299件，n = 18件3025205

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

17、误判率，尽量避免这样的损失,将获得更高的效益。具体方法如下：假设工序发生了故障：第一次检杳第二次检杳第三次检杳结论结论正确性发生概率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”表

18、示检测到合格品，“0”表示检查到不合格品）7.2 结论及分析在进行检查时，不能由一次检查结果妄下结论，这可能因误判造成大的损失。由以 上模型可以得出：在工序发生故障时误判率为 0.208，比检查一次的误判率 0.4 减少 0.192；在工序正常时误判率为 0.000792，比检查一次的误判率 0.02 减少 0.019208。 因为在工序故障时误判引起的损失为 200 元/件，工序正常时因误判停机造成的损失为 1500/次，这远远大于检查费用 10 元/次。另外，此工序检查的错误率得到极大降低， 对于其他工序的正常运行也有一定共贡献，带来了间接效益。所以采取适当连续多检查 几次从而大幅度减小误

19、判率的策略可以获得更高效益。8.模型的评价、改进和推广8.1 模型的评价优点：1.通过建立的模型确定了最优的检查间隔和道具更换间隔，从而将平均每一 个零件的损失成本降低；2.优化后的模型二将检查的错误率极大地降低，提高了生产效率和效益；3.由于工序故障率较低，将平均每一个合格零件的损失费用近似等于平均每 一个零件的损失费用，使模型得到简化。缺点：1. 忽略了换刀间隔和检查间隔时间内产生的零件数；2.实际生产中还要考虑流水线上的其他工序的状态与此工序的相互影响，里 没有考虑。8.2 模型的改进由于流水线上的其他工序与此工序往往是关联的，尤其是故障率和误判率对 其他工序有很大影响。因此各工序应看作

20、整体，综合起来考虑，得到的评价参数 对工业生产会更有实际指导意义。8.3 模型的推广该模型可以广泛应用于汽车、飞机等机械零部件加工、电子仪器生产等生产 车间的管理中，一方面有助于降低成本，另一方面提高机器的工作效率。参考文献：1 盛骤，谢式千,潘承毅，概率论与数理统计,北京 ：高等教育出版社，2008.6。2 孙山泽，数学的实践与认识，北京：北京大学，2000.1（1）。附录附录 1：100 次刀具故障记录（完成的零件数）459362624542509584433748815505612452434982640742565706593680926653164487734608428115359

21、3844527552513781474388824538862659775859755649697515628954771609402960885610292837473677358638699634555570844166061062484120447654564339280246687539790581621724531512577496468499544645764558378765666763217715310851附录 2：问题一 MATLAB 求解程序a=1:1500;l=1000./a+20./sqrt(1140*cxs(a)./(cxs(a)+11400)+31/114+310

22、0./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:m-1 k=i*(cxe(i)-cxe(i-1)+k;end u=(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

23、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,900)x =342.0000feval =5.2972exitflag =1output =iterations: 31funcCount: 32algorithm: golden section search, parabolic interpolationmessage: 1x112 char附录 3：问题二 MATLAB 求解程序：a=1:1500

24、;l=1000./a+(1./sqrt(1140*cxs(a)./(11400+cxs(a).*(10+30*(1-1./cxs(a)-1/114 00).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(每个零件的平均管理成本/元)xlabel(

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

26、40*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)4sqrt(1140*cxs(x)./(11400+cxs(x)+(1140+cxs(x)./(114*cxs(x).*(sqrt(1140*cxs(x)./(11400+cxs(x)+1)+400/3*sqrt(11400+cxs( x)./(114000*cxs(x)+30*(11400+cxs(x)./(114*cxs(x) x,feval,exitflag,output=fminbnd(fun,0,900)x =299.0001feval =7.3814exitflag =1output =iterations: 28funcCount: 29algorithm: golden section search, parabolic interpolation message: 1x112 char