自动化车床管理

1、自动化车床管理摘要本文讨论了自动化车床连续加工零件的工序定期检查和刀具更换的最优策略，届丁优化问题。我们根据原始数据利用 EXCE软件进行统计分析，得出刀具 正常工作时间长的函数，建立了以期望损失费用为目标函数的数学模型。对丁问题一，我们假设所有的检查都为等间距检查， 以检查到的零件是否为 次品来判定工序是否正常，出现了故障进行维修并使其恢复正常， 若一直未出现 故障则当加工到定期换刀时刻就换刀， 利用概率论的相关知识，求出一个周期内 的期望损失费用E(L )，再求出期望零件个数E(T )，建立了以每个零件的期望损失费用L为目标函数的随机优化模型，利用MATLAB解得出检查间隔t。=27 ,换

2、刀间隔t =270，每个零件的期望损失费用L= 5.7742。对丁问题二，根据题目信息，不管工序是否正常都有可能出现正品和次品， 在问题一的基础上，利用概率知识调整了检查间隔中的不合格品所带来的损失费 用，同时加上了因误检停机而产生的费用，求出一个周期内的期望损失费用E (L ),再求出零件个数E (T ),建立了以每个零件的期望损失费用L为目标的随机优化模型，利用MATLAB解得出检查问to =46，换刀间隔t276 ,每个零件的期望损失费用L =10.3945。对丁问题三，在问题二的基础上，我们利用概率论的相关知识，将工序正常 工作的时间长由开始的近似等丁刀具无故障工作的时间长，改进为刀具

3、无故障工作时间长的95%其它的故障近似服从均匀分布，求出一个周期内的期望损失费 用E (L ),再求出零件个数E (T ),建立了以每个零件的期望损失费用L为目标的随机优化模型，利用MATLAB解得出检查问to =30，换刀间隔ti = 600 ,期望损失费用 L =10.4212。最后我们对模型进行了改进与推广，其中模型的推广对多道工序和多个零件 的复杂车床管理系统的生产有一定的指导意义。关键字：期望损失费用随机优化模型检查间隔换刀间隔一、问题重述1.1问题背景随着现代工业化的发展，自动化设备越来越多地在人们的生活、 生产中得到 应用。在工厂中，自动化车床的使用非常普遍，许多工序用自动化车床

4、连续加工 某种零件等。自动化车床生产有很多优点，只需要人工进行控制，定期检查进行 维修即可。生产的产品一般比人工生产的更快更精确。一般的车床生产都与刀具 有关，刀具的好坏直接影响到产品的优劣。由刀具引起的故障占机器故障的大部 分，及时地发现刀具故障对减少损失有重大的意义，同时如何制定有效的换刀策略是另外一个增加收益的方法。1.2该工序的相关信息现有一道工序用自动化车床连续加工某种零件，由于各种因素的影响该工序会出现故障，其中刀具损坏故障占95%,其它故障占5%假定在生产任一零件时 出现故障的机率均相同。工作人员通过检查零件来确定工序是否出现故障。现有100次刀具故障记录，故障出现时该刀具完成的

5、零件数如附表。现计划在刀具加 工一定件数后定期更换新刀具。生产过程中的部分费用如下：故障时产出的零件损失费用f=200元/件；进行检查的费用t=10元/次；发现故障进行调节使恢复正常的平均费用 d=3000元/次（包括刀具费）;未发现故障时更换一把新刀具的费用 k=1000元/次。附:100次刀具故障记录（完成的零件数）4593626245425095844337488155056124524349826407425657065936809266531644877346084281153593844527552513781474388824538862659775859755649697515

6、6289547716094029608856102928374736773586386996345555708441660610624841204476545643392802466875397905816217245315125774964684995446457645583787656667632177153108511.3需要解决的问题：题中附录给出了 100次刀具故障记录，根据生产过程中的部分费用等信息 我们需要采用数学建模的方法来帮助解决以下问题：问题一：假定工序故障后产出的零件均为不合格品，正常时产出的零件均 为合格品，根据相关条件，对该工序设计效益最好的检查间隔（生产多少零件检

7、查一次）和刀具更换策略。问题二：如果该工序正常时产出的零件不全是合格品，有2利不合格品；而工序故障时产出的零件有40利合格品，60利不合格品。工序正常而误认有 故障停机产生的损失费用为1500元/次。根据相关条件，对该工序设计效益最好 的检查间隔和刀具更换策略.问题三：在问题二的基础上，从改进检查方式这个角度重新考虑，再建立 更优化模型，使该工序获得更高的效益。二、问题分析由题中信息可知，由于刀具损坏等原因会使工序出现故障，工序出现故障完 全是随机的，即在生产任意一个零件时都有可能发生故障。工作人员通过检查零件来确定工序是否出现故障,如果检查过于频繁,那 么工序就会经常处于正常状态而少生产出不

8、合格品，然而，这将使检查费用过 高；检查间隔过长，虽然可以减少检查费用，但由于不能及时发现故障而可能导 致大量不合格品出现，必将提高每个零件的平均损失费用。根据题目信息，刀具加工一定件数的零件后将定期更新刀具， 从而我们可以 通过确定最佳检查间隔和换刀问隔来减少损失。2.1对问题一的分析根据题目要求，我们假定所有的检查都为等间隔检查， 因为未发生故障时生 产的零件都是合格品，所以当发现零件不合格时就认为工序发生了故障， 从而停 机检查并使其恢复正常。若一直未发生故障，则当加工到定期更换刀具时刻，不 管是否发生了故障都进行换刀。计算平均费用可分为两种情况：（ 1）在换刀之 前未发生故障，记平均损

9、失费用为L1 , （2）在换刀之前发生了故障，记平均损失费用为Ln0然后以每个零件的期望损失费用为目标函数，运用 MATLAB等软件 进行编程求解使其最小。2.2对问题二的分析根据题目中所给的条件，我们还是假定所有的检查都为等间隔检查， 因为未 发生故障时次品率为2% ,发生故障时的正品率为40% ,所以不能单凭是否检查 到次品来判定工序是否正常，在工序正常时有可能误判，这样就会产生误检停机 费用，计算平均费用分为两种情况：（1）在换刀之前未发生故障，损失费用记为Pi , （2）在换刀之前发生了故障，损失费用记为 P2 ,然后以每个零件的期望损失费用为目标函数，运用MATLAB等进行编程求解使

10、其最小。2.3对问题三的分析由于刀具损坏故障占整个故障的95%，其他故障只占5% ,我们在问题二的 基础之上，可以认为整个工序正常工作的时间长近似等于刀具正常工作的时间长2的95%，即也服从正态分布N (0.95比(0.95b ),其它故障近似服从均匀分布。考虑到这点后，随着概率密度函数的改变，每个零件的期望损失费用将会改变， 最佳方案就会随之改变。三、符号说明f :每件不合格品的损失费用；t :每次检查的费用；d :发现故障进行调节使恢复正常的平均费用；k：未发现故障时更换一把新刀具的费用；t。：平均检查间隔；%:定期换刀间隔；n : 一个周期内的实际检查次数；h:工序正常而误认有故障停机产

11、生的损失费用；L :每个零件的期望损失费用；f (x ):刀具寿命的概率密度函数；x：出现故障时已经生产的零件个数；E(L): 一个周期内的期望损失总费用；E(T )：期望零件个数；s : 一个周期内的最多检查次数L1 :在定期换刀之前未发生故障的损失费用Ln :在定期换刀之前发生故障的损失费用四、模型假设1) 检查时间和换刀时间忽略不计；2) 所有的故障都为刀具故障；3) 刀具故障服从正态分布：4) 每次只抽查一个零件检查；5) s为整数,即 mod(t1,t0) = 06) 一道工序只需要一把刀具；五、模型的建立与求解5.1数据处理 5.1.1刀具正常工作的时间长的概率密度函数题中附录给出

12、了 100次刀具故障的记录，我们利用了 EXCELS件对这些数据 进行了相关的统计分析。我们采用了采用6SQ统计插件中的假设检验下的安德森 -达令正态性检验来对其进行正态分布的检验，在显著性水平a =0.1时，发现刀具故障服从正态分布N (H, cr2 ),其中 P =600, s =196.63。由此可知概率密度函数f x =.27：二刀具故障时生产的零件数和正态分布概率之间的关系如下图：其中横坐标为刀具故障时生产的零件数，纵坐标为正态分布概率|=j频数正态ta率柱状分布图卜面我们对正态分布进行检验： 利用72拟合检验法进行检验，我们用刀具寿命的最大值减去最小值，取 100为区间长度，将其分

13、成了 11个区间，分别算它们的频数，其中由丁最后两个区间的频数都为1,根据检验的原则，我们将它们合并为一个区间，再计算各数值在 区间出现的概率，其中n=100得到下表所示数据：区间频数团广吧fL明83.5 183.5391.287. 0313183.5 283.5393. 672. 4523283.5 383.56368.174.4064383.5 483.51214414.1310.1911483.5 583.52352918. 9827.8714583.5 683.52248419. 7924. 4568683. 5 783.51625616. 0215.98783, 5 883. 586

14、410. 076.3555883. 5 983. 55254. 915,0916983.5 1183.5242. 411.6598在显著性水平a =0.05下，结果如下：10 f 272 =£ 二 _n =105.4962 100 =5.4962 , Z2005 (7 ) = 14.067i x npi因为7720.05 (7 ),在可接受区间内，故服从正态分布。5.1.2刀具更换间隔在定期更换刀具之前，我们采用了等间距检查的方式对零件进行检查，若出现故障则进行调节使其恢复正常，若没有检查出故障，则到了定期更换刀具时刻进行换刀，为了简化模型，我们假定在正常换刀之前进行的是整数次检查，

15、即mod (t1 ,t0 )=0 o5.2模型一的建立与求解5.2.1模型一的建立如果在换刀之前未发生故障，则损失包括两部分：(1) 检查费用sxt ;(2) 更换刀具费用k ；则此种情况下总的损失为L1 =s><t +k ;如果在换刀之前发生了故障，此时实际检查次数为n+1 ,假设前n次检查生产的都是正品，个数为x,则次品的个数为(n+1)xt0-x,此时损失包括三部分：(1) 检查费用为(n十1;(2) 发现故障进行调节使恢复正常的费用d ；(3) 损失费用 n 1 t0 -x f ；则此种情况下总的损失费用为Ln =(n +1 )心+d + (n +1四0 -xl f期望损失

16、为：E L = "1 f x dxi Ln f x dx's t0'nn zS期望零件个数： E(T)=sKt°x 广 f(xpx+£"°(n +1 阡 x f (x )dxst。n A n t°每个零件的期望损失费用：L =E T即min =L =，要使期望损失费用达到最低，则等价丁求最佳的to,ti，使L达到最小。5.2.2模型一的求解利用MATLAB对上述模型进行求解，可得到to =27 , ti = 270 , L =5.7742。即每生产27个零件检查一次，生产270个零件后进行定期换刀，每个零件的期 望损失

17、费用为5.7742。5.2.3模型一的评价此模型采用了等间隔检查的方式，简化了模型，其中对损失费用分两种情况 讨论，简单明了，易丁理解。将求解最佳检查间隔和换刀间隔转化求解最小期望 损失费用，使模型的目标性更强。但此模型为等间隔检查，在两次检查中次品率 可能会很高，这样次品损失费用就会增加。5.3模型二的建立与求解5.3.1模型二的建立如果在换刀之前未发生故障，则在换刀时刻已经生产的零件个数为ti ,根据题中的信息，这些零件中的废品率为 2%,则损失费用包括四部分：(1) 检查费用：sxt ;(2) 误检停机费用：hxsx2% ；(3) 正常工作时的次品损失费用：ti X f X2% ；(4)

18、 更换刀具费用：k ；则此种情况下总的损失费用为：L1 =sKt +k +t1 X f 乂2% +hKs2%如果在换刀之前已经发生故障，假设第n +1次检查出故障，则此时已经生产的零 件个数为：(n十1 )Mo ,前n次检查都是正常工作的，未发生故障时生产的零件个数为x (x =n勺。)，发生故障后生产的零件个数为：(n+1to -x ,根据题中信息，我们可知，正常工作时次品率为 2%,发生故障时次品率为60%,则损失费用包 括五部分：(1) 检查费用：(n +1 ><t ;(2) 误检停机费用：n x h x 2% ;(3) 正常工作时的次品损失费用：xx f X2% ；(4)

19、发现故障进行调节使恢复正常的费用d ;(5) 发生故障后次品的损失费用：n 1% -x I： f 60%则此种情况下总的损失费用为：L n = n 1 t h n 2% - x f 2% 牌n 1t0 x | ： f 60%期望损失为:E L = .st0L1 fn 1x dxi Ln f x dxn -0期望零件个数:E T =s t0COst0fs n 1 t0x dx I j n 1t0f x dxn 苴 n>t0每个零件的期望损失费用：L = E LE T即min =L =我，要使平均损失费用达到最低，则等价丁求最佳的t0,t1，使L达到最小。5.3.2模型二的求解利用MATLA

20、B对上述模型进行求解，可得到t0 =46 , t1 = 276 , L =10.3945即每生产46个零件检查一次，生产276个零件后进行定期换刀，每个零件的期 望损失费用为10.3945。5.3.3模型二的评价此模型考虑了工序正常工作和工序出现故障时产生的次品率，利用概率的相关知识，对模型简化，最后将求解最佳检查间隔和换刀间隔转化求解最小期望损 失费用，使模型的目标性更强。但模型考虑的是等间隔检查，当发现次品时就停 机检查，误检停机费用会增加，发现正品时就继续生产，次品损失费会增加。5.4模型三的建立与求解5.4.1 模型三的建立在前两个模型中，我们都认为所有的故障都为刀具故障，即整个工序正

21、常工 作的时间长近似等丁刀具正常工作的时间长，而事实上刀具故障只占整体故障的95%,在这个基础上，我们认为整个工序正常工作的时间长是刀具正常工作时间2长的95%，即它服从正态分布N (0.95卜(0.95。)。而其它故障近似服从均匀分 布。假设此时工序正常工作的时间长的概率密度函数为 P(x),则根据模型二中 的相关知识，建立的模型如下：如果在换刀之前未发生故障，则在换刀时刻已经 生产的零件个数为,根据题中的信息，这些零件中的废品率为 2%,则损失费 用包括四部分：(1) 检查费用：sxt ;(2) 误检停机费用：hxsx：2% ;(3) 正常工作时的次品损失费用：ti尺f 乂2% ；(4)

22、更换刀具费用：k ；则此种情况下总的损失费用为：L1 =st + k +t1 x f 乂2%十hXsX2%如果在换刀之前已经发生故障，假设第n +1次检查出故障，则此时已经生产 的零件个数为：(n +1)xt° ,前n次检查都是正常工作的，未发生故障时生产的零 件个数为x (x =n 7。)，发生故障后生产的零件个数为：(n+1户。x，根据题中信息，我们可知，正常工作时次品率为2%,发生故障时次品率为60%,则损失费用包括五部分：(1) 检查费用：(n +1 /t ;(2) 误检停机费用：n x h x 2% ；(3) 正常工作时的次品损失费用：x f 2% ;(4) 发现故障进行调

23、节使恢复正常的费用 d ;(5) 发生故障后次品的损失费用：n 1 t0 -x I： f 60%则此种情况下总的损失费用为：Ln =n 1 t h n 2% x f 2%，* n 1 t0 - x " f 60%期望损失为：E L ：= r L1 P x dxi Ln P x dxs.to- nn =0期望零件个数：E (T ) = sKt0P(x dx 十£ K* 先(n +1 V t0 x P (x )dxs t0n z0 n t0期望损失费用：l=ELE T即min =L =葺号，要使平均损失费用达到最低，则等价丁求最佳的t°,ti，使L达到最小。5.4.2

24、模型三的求解利用MATLAB对上述模型进行求解，可得到t0 =30 , t1 = 600 , L =10.4212即每生产30个零件检查一次，生产 600个零件后进行定期换刀，期望损失费 用为 10.4212。5.4.3模型三的评价此模型考虑到了工序出现故障并不是完全由刀具故障引起的，对概率密度函数进行了修正，最后将求解最佳检查间隔和换刀间隔转化求解每个零件的最小期 望损失费用，使模型的目标性更强。但此模型未将其它故障的概率算入工序发生 故障的概率中，实用价值较差。六、模型的优缺点6.1模型的优点1、本文建模思想易丁理解，模型操作性强；2、对零件的检查采取了等间隔抽查，简化了模型，使模型便丁建

25、立和求解；3、将每个零件的平均损失费用作为目标函数，建立了评估体系，既有利丁求出模型的最优解，乂比较符合实际生产中企业取舍方案的标准；6.2模型的缺点1、我们没有对模型进行模拟仿真：2、在模型一和模型二中，我们忽略了其他导致故障发生的原因，只考虑了刀具 故障。七、模型的改进与推广7.1模型的改进对丁问题二，由丁工序正常时产出的零件仍有2%为不合格品，而工序故障 时产生的灵件有40%为合格品，这样工作人员在通过定期检查单个零件来确定 工序是否出现故障的检查方式必然会导致两种误判(1)正常工序时因检查到不 合格零件而误认为出现故障；(2)工序发生故障后检查到的仍是合格品而认为 工序正常，这两种情况

26、都将造成很大损失.我们建议采取不等间距的检查方式， 分为以下几种情况：(1) 连续两次检查都为正品时，我们认为工序正常，继续生产；(2) 连续两次检查都为次品时，我们认为工序发生故障，进行维修使其恢复正(3) 常后再生产；(4) 连续两次检查中一次为正品，另一次为次品时，继续第三次检查，再进行判断；这样虽然会相应地增加检查费用，但大大降低了因误检而造成的损失，从 而使系统工序获得更高的效益.7.2 模型的推广本文所建立的模型针对的是单道工序加工单一零件的情况， 但可以扩展到多 道工序和多个零件的复杂车床管理系统。在多道工序中，我们可以通过统计分析 各道工序发生故障的概率，为工序的检查确定优先级

27、，有效的控制故障发生的次 数。当生产多个零件时，可以采用随机检查一组零件的方式， 把这组零件当做一 个整体来考虑，再综合利用以上模型，就可以求解出最优的检查间隔和换刀间隔。八、参考文献1 盛骤谢式千概率论与数理统计 高等教育出版社.2 蔡俊 可靠性工程学 黑龙7工科学技术出版社.3 沈玉波 冯敬海可修系统的最优检测更新模型数学的实践与认识4 朱道元 数学建模案例精选 科学出版社5 戴朝寿 孙世良数学建模简明教程 高等教育出版社6 楼顺天 陈生潭 雷虎明MATLAB5.X程序设计语言西安电子科技大学出版社7 宋来忠 王志明数学建模与实验 科学出版社九、附录模型一求解的MATLABS代码：k=1;for a=84:600for b=1:a-1if mod(a,b)=0p=normcdf(a,600,196.62917);c=1000+a/