【文章內(nèi)容簡(jiǎn)介】 其誤差向量如下：s = R: [3x3 double] df: 2 normr: 由此可見(jiàn)，曲線1的嚙合誤差很小。對(duì)于曲線2，其誤差向量如下：s = R: [3x3 double] df: 2 normr: 由此可見(jiàn)，曲線2的嚙合誤差很小。對(duì)于曲線3，其誤差向量如下：s = R: [3x3 double] df: 2 normr: 由此可見(jiàn)，曲線3的嚙合誤差很小。對(duì)于曲線4，其誤差向量如下：s = R: [3x3 double] df: 2 normr: 由此可見(jiàn)，曲線4的嚙合誤差很小。對(duì)于曲線5，其誤差向量如下：s = R: [3x3 double] df: 2 normr: 由此可見(jiàn)，曲線5的嚙合誤差很小。當(dāng)HBS1和HBS2≤350或HBS1＞350,HBS2≤350時(shí)，在MATLAB中程序?yàn)閏lose allclcclearhold onx=::y1=[ ]y2=[ ]y3=[ ]y4=[ ]y5=[ ][p1,s]=polyfit(x,y1,2)[p2,s]=polyfit(x,y2,2)[p3,s]=polyfit(x,y3,2)[p4,s]=polyfit(x,y4,2)[p5,s]=polyfit(x,y5,2)khbeta11=polyval([ ],x)khbeta12=polyval([ ],x)khbeta13=polyval([ ],x)khbeta14=polyval([ ],x)khbeta15=polyval([ ],x)plot(x,y1,39。r39。,x,khbeta11,39。k39。)plot(x,y2,39。r39。,x,khbeta12,39。k39。)plot(x,y3,39。r39。,x,khbeta13,39。k39。)plot(x,y4,39。r39。,x,khbeta14,39。k39。)plot(x,y5,39。r39。,x,khbeta15,39。k39。)xlabel(39。φd39。)ylabel(39。khβ39。)title(39。齒輪傳動(dòng)的不均勻系數(shù)39。)gtext(39。139。)gtext(39。239。)gtext(39。339。)gtext(39。439。)gtext(39。539。)其運(yùn)行結(jié)果如圖3下面進(jìn)行誤差分析對(duì)于曲線1，其誤差向量如下：s = R: [3x3 double] df: 4 normr: 由此可見(jiàn)，曲線1的嚙合誤差很小。對(duì)于曲線2，其誤差向量如下：s = R: [3x3 double] df: 4 normr: 由此可見(jiàn)，曲線2的嚙合誤差很小。對(duì)于曲線3，其誤差向量如下：s = R: [3x3 double] df: 4 normr: 由此可見(jiàn)，曲線3的嚙合誤差很小。對(duì)于曲線4，其誤差向量如下：s = R: [3x3 double] df: 4 normr: 由此可見(jiàn)，曲線4的嚙合誤差很小。對(duì)于曲線5，其誤差向量如下：s = R: [3x3 double] df: 4 normr: 由此可見(jiàn)，曲線5的嚙合誤差很小。function [khb1,khb2]=gearkhb(yd1,yd2,yd3,yd4,zhch1,zhch2,fd)if yd1350amp。amp。yd2350 %計(jì)算載荷分布不均系數(shù)khb1 if zhch1==1 khb1=khb21(fd) elseif zhch1==2 khb1=khb22(fd) elseif zhch1==3 khb1=khb23(fd) elseif zhch1==4 khb1=khb24(fd) elseif zhch1==5 khb1=khb25(fd) endelseif (yd1=350amp。amp。yd2=350)||(yd1350amp。amp。yd2=350) if zhch1==1 khb1=khb11(fd) elseif zhch1==2 khb1=khb12(fd) elseif zhch1==3 khb1=khb13(fd) elseif zhch1==4 khb1=khb14(fd) elseif zhch1==5 khb1=khb15(fd) endendif yd3350amp。amp。yd4350 %計(jì)算載荷分布不均系數(shù)khb2 if zhch2==1 khb2=khb21(fd) elseif zhch2==2 khb2=khb22(fd) elseif zhch2==3 khb2=khb23(fd) elseif zhch2==4 khb2=khb24(fd) elseif zhch2==5 khb2=khb25(fd) endelseif (yd3=350amp。amp。yd2=350)||(yd4350amp。amp。yd2=350) if zhch2==1 khb2=khb11(fd) elseif zhch2==2 khb2=khb12(fd) elseif zhch2==3 khb2=khb13(fd) elseif zhch2==4 khb2=khb14(fd) elseif zhch2==5 khb2=khb15(fd) endend 動(dòng)載荷系數(shù)Kv2的計(jì)算在MATLAB中Kv程序?yàn)閒unction kv2=gearKV2(jd,v2)global jdglobal v2if jd=5 kv2=elseBB2=*()^AA2=50+56*()kv2=(AA2./(AA2+(200.*v2).^)).^(BB2)%kv2=fix(kv2*1000)/1000End 彈性影響系數(shù)Ze2的計(jì)算在MATLAB中，的程序?yàn)椋篺unction ze2=gearZE2(l1,l2)global l1global l2x=[ 。 0 。 0 0 。 0 0 0 ]。ze2=x(l1,l2)。5 VB程序?qū)崿F(xiàn) 優(yōu)化設(shè)計(jì)功能的實(shí)現(xiàn)通過(guò)點(diǎn)快捷方式按鈕“”，進(jìn)入減速器優(yōu)化設(shè)計(jì)主界面如圖4，在主界面點(diǎn)擊“確定”按鈕進(jìn)入優(yōu)化界面如圖5。輸入?yún)?shù)時(shí)需查表確定的點(diǎn)擊“幫助”按鈕如圖5。在圖中輸入已知參數(shù)和設(shè)計(jì)變量初始值后，點(diǎn)擊“運(yùn)行”按鈕，通過(guò)調(diào)用Matlab優(yōu)化設(shè)計(jì)模塊，自動(dòng)完成各零件的優(yōu)化設(shè)計(jì)計(jì)算，并將結(jié)果輸出到界面中如圖6。若計(jì)算結(jié)果不理想可以修改輸出結(jié)果中的個(gè)別參數(shù)，如果想再次優(yōu)化，點(diǎn)擊“重設(shè)”按鈕；點(diǎn)擊“建?！卑粹o，可以通過(guò)Solidworks生成優(yōu)化后的齒輪三維模型。 應(yīng)用實(shí)例已知功率P=5KW，高速級(jí)小齒輪轉(zhuǎn)速n=960，總傳動(dòng)比i=12,高速級(jí)傳動(dòng)比i1=4,載荷穩(wěn)定預(yù)期壽命10年，每年300個(gè)工作日，工作時(shí)間占20%，動(dòng)力機(jī)為電動(dòng)機(jī)，工作平穩(wěn)，傳動(dòng)不逆轉(zhuǎn)，齒輪對(duì)稱分布，人工選擇齒輪材料和熱處理，高速級(jí)小齒輪40Cr調(diào)質(zhì)，硬度240HB，大齒輪45鋼調(diào)質(zhì)，硬度200HB，許用應(yīng)力為 ；低速級(jí)小齒輪40Cr調(diào)質(zhì)，硬度260HB，大齒輪45鋼調(diào)質(zhì)，硬度230HB，許用應(yīng)力為 ；圖圖4. 優(yōu)化設(shè)計(jì)主界面在圖5界面中輸入已知參數(shù)和設(shè)計(jì)變量初始值，點(diǎn)擊“運(yùn)行”按鈕，經(jīng)過(guò)幾秒鐘等待，即可輸出優(yōu)化結(jié)果如圖6，選擇高速級(jí)齒輪副，點(diǎn)擊“建模按鈕”即可生成三維模型如圖7.表3優(yōu)化設(shè)計(jì)與傳統(tǒng)設(shè)計(jì)比較設(shè)計(jì)方法高速級(jí)齒輪模數(shù)m1低速級(jí)齒輪模數(shù)m2高速級(jí)小齒輪齒數(shù)z1低速級(jí)小齒輪齒數(shù)z3高速級(jí)齒輪傳動(dòng)比i1中心距A(mm)優(yōu)化率傳統(tǒng)設(shè)計(jì)3320364366優(yōu)化設(shè)計(jì)32033317結(jié)論通過(guò)本次畢業(yè)設(shè)計(jì)我完成了減速器的優(yōu)化設(shè)計(jì)，主要是二級(jí)減速器齒輪傳動(dòng)部分的優(yōu)化設(shè)計(jì)。通過(guò)以中心距最小為目標(biāo)函數(shù)建立數(shù)學(xué)模型，對(duì)與齒輪相關(guān)的系數(shù)進(jìn)行參數(shù)化處理來(lái)達(dá)到優(yōu)化的目的，最后通過(guò)Solidworks完成了三維造型。由于時(shí)間緊，任務(wù)重，而且是初次設(shè)計(jì)，所以難免存在一些問(wèn)題，特別是在一些參數(shù)的選定方面需要改進(jìn)，這樣可以大幅度的提高設(shè)計(jì)精度。例如：在利用MATLAB進(jìn)行曲線嚙合時(shí)還應(yīng)該多嘗試幾次其他的函數(shù)，然后進(jìn)行誤差分析與比較，經(jīng)比較后選擇誤差較小的函數(shù)來(lái)進(jìn)行曲線嚙合；在設(shè)計(jì)中個(gè)別的參數(shù)使用了定值(s=1,sf=)；在設(shè)計(jì)兩級(jí)齒輪減速器中沒(méi)有對(duì)軸、軸承等進(jìn)行設(shè)計(jì)。在以后盡量逐步完善使優(yōu)化更精確。圖5. 參數(shù)輸入界面圖6. 優(yōu)化結(jié)果界面圖7. 高速級(jí)齒輪副 參考文獻(xiàn)[1] [M].北京:機(jī)械工業(yè)出版社,1990.[2] 郇義鵬, 與 Matlab編 程接口應(yīng)用[J]，清華大學(xué)學(xué)報(bào)，2003,（5）：3438.[3] 朱志松,郭曉麗,[J].計(jì)算機(jī)應(yīng)用，2003,（9）：1819.[4] 劉長(zhǎng)榮,[M].北京：中國(guó)農(nóng)業(yè)科技出版社,2002.[5] [J].計(jì)算機(jī)應(yīng)用與軟件實(shí)現(xiàn)，2004,（12）：5152.[6] 蘇金明,黃國(guó)明,[M].北京：電子工業(yè)出版社，1999.[7] 和MATLAB 的混合編程優(yōu)化設(shè)計(jì)[M].北京：電子工業(yè)出版社，2007.[8] [M].北京:機(jī)械工業(yè)出版社,2000. [9] 馬興義, 應(yīng)用開(kāi)發(fā)指南[M].北京：機(jī)械工業(yè)出版社,2002. [10] 劉炳文,Visual Basic 程序設(shè)計(jì)基礎(chǔ)[M].北京：人民郵電出版社,1998111 [11] Farhi Marir，Ian Watson. Case一based reasoning:a catego2 rized bibliography[J]. The Knowledge Engineering Review. 1994，9 (4) :355381.[12] Nbolunt G,Clarke S. Artifical Intelligent and Design Automation System[J]. Journal of Engineering Design. 1994，5(4): 299314.[13] [D].蘭州：蘭州理工大學(xué),2006,110致謝四年大學(xué)生活即將結(jié)束，現(xiàn)在回想起來(lái)，在這四年當(dāng)中還是有很多學(xué)習(xí)時(shí)間沒(méi)有抓住，有很多知識(shí)沒(méi)有掌握。但總體來(lái)說(shuō)這四年大學(xué)生活是豐富多彩的，這要感謝大學(xué)生活中辛勤教育我的各位老師，無(wú)論在學(xué)習(xí)還是在生活以及做人方面各位老師都給予了指導(dǎo)，使我們受益匪淺。系里的領(lǐng)導(dǎo)和老師給我們提供了良好的環(huán)境，特別是指導(dǎo)教師張小芹，給予我極大的幫助，包括搜集資料，設(shè)計(jì)方案的提出，設(shè)計(jì)過(guò)程中的指導(dǎo)，設(shè)計(jì)結(jié)果的審核等，在此我對(duì)系里的領(lǐng)導(dǎo)以及老師們致以衷心的感謝。此外，各位同學(xué)也給了我很大的幫助，也同樣非常感謝他們。 通過(guò)這次設(shè)計(jì)，收獲很大，同時(shí)我也意識(shí)到了自身知識(shí)的貧乏，實(shí)踐能力的不足。在以后的學(xué)習(xí)和工作中，我會(huì)進(jìn)一步提高自身的能力，提高自己的素質(zhì)。附錄MATLAB總程序global i %總傳動(dòng)比global i1 %高速級(jí)傳動(dòng)比global m1 %模數(shù)global m2 %模數(shù)global p %功率global n %轉(zhuǎn)速global ze1 %彈性影響系數(shù)global ze2 %彈性影響系數(shù)global zh %區(qū)域系數(shù)global s %接觸疲勞安全系數(shù)global sf %彎曲疲勞安全系數(shù)global j1 %接觸疲勞極限global j2 %接觸疲勞極限global w1 %彎曲疲勞極限global w2 %彎曲疲勞極限global w3