自动化车床管理问题详细解答含编程

自动化车床管理承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：三峡大学参赛队员(打印并签名)：1.2.3.指导教师或指导教师组负责人(打印并签名)：教师指导组日期：2010年7月18日赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：摘要本文要解决的是在一定的条件下如何安排刀具的更换策略和检查的次数与间隔使每个正品的平均消耗最低的问题，属于优化问题中的概率数理统计问题。通过对附加数据的数理统计分析，我们发现故障发生时所完成的零件数符合正态分布。每次换刀由两部分组成，第一部分为固定周期为T的换刀策略，第二部分为在进行零件的检查过程中发现故障后便立刻进行维修。检查零件我们采用非等间距方式，而后以生产每个正品所损耗的费用作为目标函数，最后利用计算机对第一部分的固定周期与检查间隔用穷举法一一比较各种情况并求解出最优解，最终算得的结果如下:第一问中平均每个正品消耗费用为4.1615，固有换刀周期为369，换刀次数为18.第二问中平均每个正品消耗费用为90268，固有换刀周期为306，换刀次数为28.【关键词】数理统计非等间距穷举法正态分布等比间距1、问题重述一道工序用自动化车床连续加工某种零件，由于刀具损坏等原因该工序会出现故障，其中刀具损坏故障占95%,其它故障仅占5%。工序出现故障是完全随机的,假定在生产任一零件时出现故障的机会均相同。工作人员通过检查零件来确定工序是否出现故障。现积累有100次刀具故障记录，故障出现时该刀具完成的零件数如附表。现计划在刀具加工一定件数后定期更换新刀具。已知生产工序的费用参数如下：已知生产工序的费用参数如下：1)故障时产出的零件损失费用f=200元/件；2)进行检查的费用t=10元/次；3)发现故障进行调节使恢复正常的平均费用d=3000元/次(包括刀具费)；4)未发现故障时更换一把新刀具的费用k=1000元/次；问题一：假定工序故障时产出的零件均为不合格品，正常时产出的零件均为合格品,试对该工序设计效益最好的检查间隔（生产多少零件检查一次）和刀具更换策略。问题二：如果该工序正常时产出的零件不全是合格品，有2%为不合格品；而工序故障时产出的零件有40%为合格品，60%为不合格品。工序正常而误认有故障仃机产生的损失费用为1500元/次。对该工序设计效益最好的检查间隔和刀具更换策略。问题三：在问题二的情况下,可否改进检查方式获得更高的效益。一百次故障记录数据如下：45936262454250958443374881550561245243498264074256570659368092665316448773460842811535938445275525137814743888245388626597758597556496975156289547716094029608856102928374736773586386996345555708441660610624841204476545643392802466875397905816217245315125774964684995446457645583787656667632177153108512、问题分析通过对车床故障数据的数理统计我们发现车窗出现故障时所生产的零件数符合正态分布函数，由于车床的故障由两方面造成，来自于95%刀具故障和来自与5%的其他因素，为了简化模型我们抓住主要因素忽略次要因素，认为5%的其他因素最终反映到对刀具的影响上，又由于车床的故障最终导致产品的质量，所以我们通过制定合适的检查方案可以及时的发现故障也可由检查的动态策略制定刀具的更换策略。对刀具的100次故障统计量如下：数据个数100平均值600样本标准偏差196.6291695最小值84最大值1153小数位数0区间个数20区间宽度53.5对于问题一：产品的合格与否直接来自与车床是否有故障，是一一对照关系，可以通过对零件的非等距检查来判断车床是否有故障，同时求出正品的期望值，产品的期望值，检查次数的期望值，废品的期望值(为产品的期望值-正品的期望值),提前更换刀具所花费用的平均期望值，及时更换刀具所花的费用期望，最后用平均花费量=生产所有产品所花的费用/正品的个数作为衡量方案好换的标准。对于问题二：与问题一比较我们发现其变化表现在产品的合格与否不和车床故障与否存在严格的一对一关系。也就导致了误判的发生，误判主要出现在两种情况，第一种当车床无故障时，存在2%为不合格品，我们采用一次检查一个零件的方法也就存在一定的误判可能性，把好的刀具当成了坏的刀具，将导致刀具的使用平均寿命减少；当车床有故障时，存在60%的不合格也存在40%的合格率，这是一个相对较相近的概率这也就大大增加了误判的可能性，使我们错误的把坏的刀具当成了好的刀具，这就加大了零件的损失率，当发现后更换时同样也是一笔客观的费用。我们同样可以求出正品的期望值，产品的期望值，检查次数的期望值，废品的期望值(为产品的期望值-正品的期望值),提前更换刀具所花费用的平均期望值，及时更换刀具所花的费用期望，误判的期望值，最后用平均花费量=生产所有产品所花的费用/正品的个数作为衡量方案好换的标准。3、模型的假设与说明(1)在生产过程中不考虑故障因素中的5%对刀具的影响；(2)每次发现故障都能够及时在生产下一个零件前对故障进行维修好；(3)发生误判时，只考虑停机生产的损失费用，1500元/次；(4)发生故障并进行调整好的平均费用是指所有可能情况的平均费用；(5)假设题目中的数据与前提都是正确的；(6)假设在进行检查时对生产速度.、质量无影响；(7)刀具的损坏时所生产的零件数所表现的频率表现出正态分布的特性；4、符号说明T相邻两次定期换刀之间生产的零件个数；0T首次检查时已经生产的零件个数；kT第k次检查所对应的零件数(k取1、2、3……)；t检查费用10元/次；d未发现故障换刀费用；()fx概率密度；fT分布密度；N总零件数；1U换刀前没有发生故障时总损失费用期望；2U换刀前已发生故障时总损失费用期望；U总的损失费用；1P在T时间段内因误判而导致的提前换刀的损失费用；2P在T时间段内因误判而导致的提后换刀的损失费用；n检查的次数；q等比数列的公比5、模型一的建立与分析正态性检验1、根据大样本场合(n50)的D检验验证刀具寿命记录概率的概率分布方式。1322112niiniinixDnxx0.282094790.02998598DnY正态分布的拒绝域为{Y=2Y或者Y=12Y}其中取=0.05，已知n=100则拒绝域，则拒绝域为{Y=-0.24或Y=1.31}其中Y=-1.2933，由于样本落入拒绝域，故在=0.05时可以认为刀具故障记录满足正态分布。2、概率密度函数求解3、于是失效概率密度函数为22(1)221()(,,)2xNftftxe其图像如下：其中11niixxN22211()1niisxxN由统计概率得x=600=196.6291695累计失效密度函数为20()(,,)()tNFtFtxfxdx其图像如下：相应的数据处理在附录中目标函数：合格零件单位期望损失：U=生产所有产品所花的费用U总/正品的个数N总假设自动化车床在连续运行中将发生N次更新刀具过程，可分为两种情况：1.换刀间隔T前尚未出现故障发生这种情况的更新间隔均为T，出现的次数等于刀具更新的总次数乘以以T为更新间隔情况下换刀前仍未出现故障的概率，即[1()]NFT,因此定期换刀前未出现故障的情况下的总损失等于这种情况下的道具更换次数[1()]NFT乘以单位更新过程中的损失费用1P,即11[1()]UNFTP2.换到间隔T前就出现故障，这是在故障发生后进行检查并进行维修换刀。这种情况下总的发生次数等于刀具更新的总次数乘以发生这种情况的概率即：()NFT，因此定期换刀前未出现故障的情况下的总损失等于这种情况下的道具更换次数()NFT乘以单位更新过程中的损失费用2P，即22()UNFTP而12UUU下面我们将对一个换刀间隔T的研究来求1P,2P1.求换刀间隔T前尚未出现故障时一次更新过程所消耗的费用1P(1)检查费用：等于检查的次数乘以单次检查所需的费用，即nt我们假设检查间隔成等比数列递减，则有0(1)1nTqTq其中0T为首次检查间隔，kT为第K次检查间隔，则有1kkTqT（2）.换刀费用k（3）不合格零件损失费用为0所以1PKnt于是，再换刀前出现故障的情况下总的损失费用1U为：一个换刀周期内的损失乘以这种情况下的可能发生的总的更新次数，即1[1()][0]UNFTKnt其中()FT为以T为更新周期情况下出现故障的概率，即为0到T内的累计失效密度函数0()()TFTfxdx2.求换刀间隔T前出现故障时一次更新过程所消耗的费用2P（1）换刀费用d（2）故障维修前所有损失费用：由于故障发生的随机性，因此可以发生在T内的任何位置，这部分的损失费用等于对于周期T内任意点x处发生故障所造成的损失与在x处可能发生故障概率的乘积进行积分的结果，即0()()()TfxWxdxFT其中()()fxFT表示在一个换刀周期T内任意的x处发生故障概率任意位置发生故障概率的损失费()Wx=检查费用+零件损失费1）检查费用同1P，等于nt2）零件损失费用等于从发生故障到维修检查之间产生的不合格零件数乘以单个零件损失费用，即()Txf所以20(){()}()TfxPntTxfdxdFT于是得到换刀前已出现故障的情况下的损失总费用2U为20()()({()})()TfxUNFTntTxfdxdFT因此工序总的期望损失为12UUU工序产生的合格零件总数为：换刀前没发生故障情况产生的合格零件总数加上换刀前发生故障情况下产生的合格零件总数，即0()[1()]()()TfxNFTTNFTxdxFT工序合格零件的单位期望损失为6、模型二的建立与分析目标函数：合格零件单位期望损失：cT=生产所有产品所花的费用U总/正品的个数N总假设自动化车床在连续运行中将发生N次更新刀具过程，可分为两种情况：1.换刀间隔T前尚未出现故障发生这种情况的更新间隔均为T，出现的次数等于刀具更新的总次数乘以以T为更新间隔情况下换刀前仍未出现故障的概率，即[1()]NFT,因此定期换刀前未出现故障的情况下的总损失等于这种情况下的道具更换次数[1()]NFT乘以单位更新过程中的损失费用1P,即11[1()]UNFTP2.换到间隔T前就出现故障，这是在故障发生后进行检查并进行维修换刀。这种情况下总的发生次数等于刀具更新的总次数乘以发生这种情况的概率即：()NFT，因此定期换刀前未出现故障的情况下的总损失等于这种情况下的道具更00()[1()][0](){{()}}()()[1()]()()TTfxNFTKntNFTntTxfdxdFTUfxNFTTNFTxdxFT换次数()NFT乘以单位更新过程中的损失费用2P，即22()UNFTP而12UUU下面我们将对一个换刀间隔T的研究来求1P