車床自動化管理問題的進一步討論.doc

車床自動化管理問題的進一步討論 劉忠敏 程海川 王小勇 （重慶通信學院 400035）【摘要】本問題是隨機最優(yōu)化問題。通過分析單位正品零件的平均損耗費用，完美的建立了模型.用Mathematic求解，并用蒙特卡洛思想進行驗證，徹底解決此問題：【關鍵詞】 最優(yōu)化 、 隨機變量 、蒙特卡洛一、 符號的說明x：檢查間隔； y：刀具的更換周期； F：一個周期內所損耗的費用； t：刀具的壽命； g99(t)為密度函數(shù) H：一個周期內所生產的正品的零件數(shù) ； p=0.98； q=1-p； r=0.4； s=1-r；k=t/x； n=y/x； c：一個周期內所生產的每個正品零件所擔負的平均損耗費用；二、 模型假設只在x的倍數(shù)處檢查；5%非刀具故障符合 0 22800 上的均勻分布。其它略。三、 問題分析通過對99年數(shù)學建模A題的分析。知刀具出現(xiàn)故障符合正態(tài)分布tN（，2），其中=600、=196.629。為了使生產一定數(shù)量零件所損耗的費用盡量少，我們可以考察在一個周期內所生產的每個正品零件所負擔的費用，只要該費用最少，則生產一定數(shù)量零件所損耗的費用最少，因此需要先給周期下一個明確的定義?！局芷凇浚簱Q上新刀具開始生產至此刀具被更換，這之間生產的零件件數(shù)。首先研究不考慮5%非刀具故障的情況。四、 模型建立 由于給出的刀具壽命t是服從正態(tài)分布的，且在一個周期內可能出現(xiàn)三種情況：刀具壽命t大于刀具更換周期y；t落于y與nx之間，nx為離y最近的一個檢查點；t落于0到nx之間；因而在建立模型時應劃分為三段考慮。模型的建立： 由題目條件設t前生產的產品均為正品，其后為次品。當?shù)毒叩膲勖笥诟鼡Q周期時，則檢查費用為：y/x10；換刀費用為1000元。當?shù)毒叩膲勖∮诟鼡Q周期時，分為：nxty和tnx兩種情況。若tnx 則檢查費用為：（t/x+1）10；次品的損失費用為：200x（t/x+1-t/x）；更換刀具的費用為：3000；若nxty 則檢查費用為：y/x10；次品的損失費用為：（y-t）200；更換刀具的費用為：1000。所以一個周期內的損失費用為：則一個周期損耗費用的平均值為：一個周期所生產的正品零件的數(shù)目為:一個周期所生產的正品零件的平均數(shù)目為:目標函數(shù)為: 模型的建立：仿照模型的原理，建立本模型,費用函數(shù)：其中：當 0tnx 時工序正常,而誤認為有故障,停機產生的損失費用的為: A通項為工序有故障,由于r=40%而i+1次誤判為正常，從而認為正常換刀形成的費用。B通項為工序有故障,由于s=60%而i次誤判為正常，第i+1次認為工序故障換刀形成費用。一個周期內生產的正品數(shù)為：AA通項為工序有故障,由于r=40%而i+1次誤判為正常，從而認為正常換刀形成的正品數(shù)。BB通項為工序有故障,由于r=40%，s=60%而i次檢查認為工序正常,但第i+1次為故障，換刀形成正品數(shù)。模型的建立：檢查時發(fā)現(xiàn)次品則再查下一個產品，若為次品則認為工序故障。仿照模型(2)的原理，建立本模型，費用函數(shù)：當 0ty時的情形。上式的d為此情形的概率。雙方還應有3000與1000的費用差別(3000-1000/2),不過這對解（x，y）引響不大。五、 模型求解用Mathematica來求解模型，利用網格法對x,y進行搜索。模型(3)的程序為：Modulet,x,y,gg99,gg98,ee,m,ee1,ee2,ee3,een,n,k,i, p=0.98,q=0.02,r=0.4,s=0.6,x=10; Whilex30,x=x+1; y=300; Whiley350, n=Floory/x;k=Floort/x; gg99=(1/(196.629*(2*Pi)0.5)*E(-(t-600)2/(2*196.6292); gg98=1510*k*q2+200*t*q+10*k+s*200*(k+1)*x-t); ee1=NIntegrategg98*gg99,t,0,x*n +SumNIntegrate(1000+s*200*(y-n*x+i*x)+(i+1)*10)*(1-s2)(i+1)* gg99,t,x*(n-i-1),x*(n-i),i,0,9 +SumNIntegrate(3000+s*200*i*x+(i+1)*20)*(1-s2)i*s2* gg99,t,0,x*(n-i),i,0,9; ee2=NIntegrate(1510*n*q2+10*n+200*t*q+s*200*(y*TBG3333t)+1000)* gg99,t,x*n,y; ee3=NIntegrate(1510*n*q2+n*10+200*y*q+1000)*gg99,t,y,2000; ee=ee1+ee2+ee3; ee99=ee; een=NIntegrate(p*t+r*(k+1)*x-t)*gg99,t,0,n*x +SumNIntegrater*(y-n*x+i*x)*(1-s2)(i+1)*gg99,t,x*(n-i-1),x*(n-i),i,0,9 +SumNIntegrate(r*i*x)*(1-s2)i*s2*gg99,t,0,x*(n-i),i,1,9 +NIntegrate(p*t+r*(y-t)*gg99,t,n*x,y +NIntegratep*y*gg99,t,y,2000; a=NIntegrategg99,t,y,2000; ee=ee/een+a*s*200*x*(0.5+Sumi*(1-s2)i,i,1,10*s2)/11400 +a*10*s2*Sum(i+1)*(1-s2)i,i,1,5/11400+2500/11400; Print ee, ee99/een, ee99, een, y, x, ee1, ee2; y=y+1 模型（1） x=18 y=359 4.81元； 無5%時： x=18 y=359 4.47元。 比較 x=18 y=360 4.93元； x=18 y=360 4.60元；模型（2） x=45 y=314 10.12元； 無5%時： x=53 y=316 9.34元。 比較 x=45 y=315 10.38元； x=52 y=314 9.6元；模型（3） x=14 y=335 9.60元； 無5%時： x=23 y=321 9.05元。 比較 x=14 y=336 9.68元； x=23 y=322 9.15元。六、 結果分析與檢驗 模型（1）也可用如下方法計算: l=1000/y+10/x+200*(x/2)/cn+(3000-a1*(1000/y)/cn1/cn=1/a+1/b a=a1+a2 b=11400, a1為發(fā)生故障時所在周期平均生產的產品數(shù), x/2為平均次品數(shù).解為x=18,y=360 l=4.96 用蒙特卡洛思想來對求解的結果進行驗證，限于篇幅我們僅給出模型(3)TC仿真程序的核心部分。while(iy+1) hh=floor(i/x); if(tt=y) if(aa9i=0) if(i/x=hh) if(i=tt) ff=ff+10.0+200+10; if(aa9i+1=0&iy) ff=ff+3000+200; break; if(i=y) ff=ff+1000; break; if(hh-i/x!=0.0) ff=ff+200.0; if(aa9i=1) if(i=y) ff=ff+1010.0; in=in+1;break; if(i-tt0 & hh=i/x) ff=ff+10.0; in=in+1; if(i-tt0&(i!=y & hh-i/x=0.0) ff=ff+10.0; in=in+1; if(i-tt0 &(i!=y&hh-i/x!=0.0) in=in+1; if(tty) if(iy) if(aa9i=1 & hh-i/x!=0.0) in=in+1; if(aa9i=1 & hh-i/x=0.0) ff=ff+10.0; in=in+1; if(aa9i=0 & hh-i/x!=0) ff=ff+200.0; if(aa9i=0 & hh-i/x=0.0) ff=ff+10.0+200+10; if(aa9i+1=0 ) ff=ff+1500; if(i=y) if(aa9i=0) ff=ff+200+1000; break; if(aa9i=1) ff=ff+1000.0+10; in=in+1;break; i=i+1; 另一個難點在構造樣本上，此處省去。要求仿真法與math法所得結果要相同。 5%時 math法 仿真法 無5%時 math法 仿真法問題（2） （45 314） 10.12 10.09 (53 316) 9.341 9.332問題（3） (14 335) 9.60 9.62 (23 321) 9.051 9.064結果吻合很好，分析雙方中間結果也如此，從而說明了模型的完整性與準確性。順便說一下，問題（3）若采用不等距檢查改進不大?！緟⒖紩俊浚?）姜啟源 數(shù)學模型 北京 高等教育出版社 1987（2）孫山澤 車床自動化管理問題中國數(shù)學學會數(shù)學的實踐與認識2000年第一期（3）張韻華 符號計算系統(tǒng)實用教程 中國科學技術大學出版社 1998年9月A Further Study of Automatic Managing of Lathe LiuZhong_Min ChengHai_Chuan WangXiao_Yong (Journal of ChongQing Communication Institute )Abstract: In t