畢業(yè)設(shè)計--單級齒輪減速器模型優(yōu)化設(shè)計

1、<p><b>　　畢業(yè)設(shè)計(論文)</b></p><p>　圖書分類號：</p><p>　密 級：</p><p><b>　　摘要</b></p><p>　　齒輪減速器是原動機和工作機之間的減速傳動裝置，它廣泛應用于日常生活和生產(chǎn)中，但傳統(tǒng)設(shè)計中仍存在著體積過大、效率過

2、低的問題，因此有必要對齒輪減速器進行優(yōu)化設(shè)計。針對一級齒輪減速器，選取合理的設(shè)計變量，以減速器體積最小為目標建立優(yōu)化設(shè)計的數(shù)學模型，并根據(jù)齒寬系數(shù)、模數(shù)、齒輪的應力和軸的彎曲強度等約束條件來確定約束函數(shù)。采用遺傳算法對模型進行優(yōu)化，經(jīng)過實數(shù)交叉、高斯變異和可行性規(guī)則選擇反復循環(huán)得出優(yōu)化后的參數(shù)，經(jīng)過機械校核得出優(yōu)化出的參數(shù)是滿足約束條件，實現(xiàn)了齒輪減速器體積最小的優(yōu)化目標。</p><p>　　關(guān)鍵詞 減速器；遺

3、傳算法；優(yōu)化設(shè)計；</p><p><b>　　Abstract</b></p><p>　　Gear reducer is the transmission device for speed deceleration, connecting prime motor and working machine，which has been widely used in o

4、ur daily life and production.But there still exist some problems in traditional design and manufacture, such as over-large volume and low efficiency. Therefore it is necessary to optimize the design of gear reducer. Afte

5、r selecting reasonable design variables to the gear reducer model, optimum mathematical model is set up aiming at the least volume. And the co</p><p>　　Keywords Gear Reducer Genetic Algorithm Optimal Desi

6、gn</p><p><b>　　目 錄</b></p><p><b>　　摘要I</b></p><p>　　AbstractII</p><p><b>　　1 緒論1</b></p><p><b>　　1.1減速器

7、1</b></p><p>　　1.1.1 減速器國內(nèi)外的現(xiàn)狀和發(fā)展1</p><p>　　1.1.2 齒輪減速器的介紹2</p><p>　　1.1.3 齒輪減速器的分類2</p><p>　　1.1.4 齒輪減速器的結(jié)構(gòu)2</p><p>　　1.1.5 齒輪減速器的工作原理3</p&

8、gt;<p><b>　　1.2優(yōu)化設(shè)計3</b></p><p>　　1.3本課題的主要內(nèi)容和解決問題4</p><p>　　2 齒輪減速器數(shù)學模型的建立5</p><p>　　2.1 齒輪減速器的數(shù)學模型5</p><p>　　2.1.1 建立一級齒輪減速器的數(shù)學模型5</p>

9、<p>　　2.1.2 確定約束函數(shù)6</p><p>　　3 遺傳算法程序設(shè)計12</p><p>　　3.1 智能算法選擇12</p><p>　　3.1.1 智能算法種類12</p><p>　　3.2 遺傳算法介紹12</p><p>　　3.3 遺傳算法程序?qū)崿F(xiàn)13</p>

10、;<p>　　3.3.1 遺傳算法中的參數(shù)13</p><p>　　3.3.2 遺傳算法的偽代碼14</p><p>　　3.3.3 編碼和解碼14</p><p>　　3.3.4 交叉15</p><p>　　3.3.5 變異16</p><p>　　3.3.6 選擇18</p>

11、<p>　　3.3.7 約束的處理19</p><p>　　4 機械設(shè)計與校核22</p><p>　　4.1 齒輪的設(shè)計與校核22</p><p>　　4.1.1 優(yōu)化結(jié)果22</p><p>　　4.1.2 大齒輪的強度的校核22</p><p>　　4.1.3 小齒輪的強度的校核24&

12、lt;/p><p>　　4.2 軸的設(shè)計和校核25</p><p>　　4.2.1 高速軸設(shè)計25</p><p>　　4.2.2 高速軸校核26</p><p>　　4.2.3 高速軸軸承壽命校核28</p><p>　　4.2.4 低速軸設(shè)計29</p><p>　　4.2.5 低速軸

13、校核30</p><p>　　4.2.6 低速軸軸承壽命校核32</p><p>　　4.3 齒輪減速器箱體的結(jié)構(gòu)尺寸33</p><p>　　4.4 減速器的潤滑和密封34</p><p><b>　　結(jié)論35</b></p><p><b>　　致謝36</b>

14、;</p><p><b>　　參考文獻37</b></p><p>　　謝謝朋友對我文章的賞識，充值后就可以下載此設(shè)計說明書（不包含CAD圖紙）。我這里還有一個壓縮包，里面有相應的word說明書（附帶：外文翻譯）和CAD圖紙。需要壓縮包的朋友聯(lián)系QQ客服1：1459919609或QQ客服2：1969043202。需要其他設(shè)計題目直接聯(lián)系?。?！ </p>

15、<p><b>　　1 緒論</b></p><p><b>　　1.1減速器</b></p><p>　　減速機是一種相對精密的機械，使用它的目的是降低轉(zhuǎn)速，增加轉(zhuǎn)矩。它的種類繁多，型號各異，不同種類有不同的用途。減速器的種類繁多[1]，按照傳動級數(shù)不同可分為單級和多級減速器；按照齒輪形狀可分為圓柱齒輪減速器、圓錐齒輪減速器和圓

16、錐－圓柱齒輪減速器；按照傳動的布置形式又可分為展開式、分流式和同軸式減速器，按照傳動和結(jié)構(gòu)特點來劃分，可分為下列六種減速器：齒輪減速器、蝸桿減速器、蝸桿齒輪減速器及齒輪-蝸桿減速器、擺線針輪減速器、諧波齒輪減速器和行星齒輪減速器。</p><p>　　對于上述六種減速器已經(jīng)有標準的系列產(chǎn)品，使用時只需要結(jié)合所需傳動功率、轉(zhuǎn)速、傳動比、工作條件和機械的總體布置等具體要求，從產(chǎn)品目錄和有關(guān)手冊中選取即可，只有在選不到

17、合適的產(chǎn)品時，才自行設(shè)計制造，此外還有一些專用的減速器（如在電梯、閥門開關(guān)等處用）。</p><p>　　1.1.1 減速器國內(nèi)外的現(xiàn)狀和發(fā)展</p><p>　　國外的減速器[2]，以德國、丹麥和日本的處于領(lǐng)先地位，它們特別是在制造工藝和材料的方面占據(jù)優(yōu)勢，它們的減速器工作具有可靠性好，使用壽命長的優(yōu)點。但其傳動形式仍是以定軸齒輪傳動為主。日本住友重工研制的FA型高精度減速器，美國Ala

18、n-Newton公司研制的X-Y式減速器在結(jié)構(gòu)上和傳動原理上都為當今世界先進的齒輪減速器。目前的減速器是向著小體積、大功率、高效率、大傳動比以及長壽命的趨勢發(fā)展。因此，除了不斷提高工藝水平、改進材料品質(zhì)外，還在傳動結(jié)構(gòu)和傳動原理上要深入探討和創(chuàng)新，平動齒輪傳動原理的出現(xiàn)就是一例。</p><p>　　國內(nèi)的減速器多以齒輪傳動、蝸桿傳動為主，但普遍存在著功率與重量比小，或者機械效率過低而傳動比大的問題。另外，工藝水

19、平和上材料品質(zhì)還有許多不足之處，特別是大型的減速器的不足之處更加的突出。60年代開始生產(chǎn)的少齒差傳動、擺線針輪傳動、諧波傳動等減速器具有傳動比大、體積小、機械效率高等優(yōu)點。但受其傳動的理論的限制，不能傳遞過大的功率，功率一般都是要小于40kw。由于在材料品質(zhì)方面、傳動的理論和工藝水平方面沒有突破，因此，國內(nèi)沒能從根本上解決傳遞功率大、機械效率高、傳動比大、重量輕、體積小等這些基本要求。改革開放以來,我國引進了一批先進的加工裝備。通過不斷

20、引進、消化和吸收國外先進技術(shù)以及科研攻關(guān),開始掌握了各種高速和低速重載齒輪裝置的設(shè)計制造技術(shù)。材料和熱處理質(zhì)量及齒輪加工精度都有較大的提高,通用圓柱齒輪的制造精度可從JB179—60的8～9級提高到GB10095—88的6級,高速齒輪的制造精度可穩(wěn)定在4～5 級。部分減速器采用硬齒面后,體積和重量明顯減小,承載能力、使用壽命、傳動效率有了大幅度的提高,對節(jié)能和提高主機的總體水平起到明顯的作用。從1988 年以來,我國相繼制定了50～60

21、 種齒</p><p>　　1.1.2 齒輪減速器的介紹</p><p>　　對于齒輪減速器，它是原動機和工作機之間的獨立減速傳動裝置[3]，由于其結(jié)構(gòu)緊湊、效率較高、傳遞運動準確可靠、使用維護方便，并可成批生產(chǎn)，所以在船舶汽車、機車、建筑用的重型機具中得到廣泛應用。并且從機械工業(yè)所用的加工機具以及自動化生產(chǎn)設(shè)備，到常見的電鐘表等，齒輪減速器都起到了重要的作用。盡管齒輪減速器廣泛的應用于日

22、常生活和生產(chǎn)中，但是它的設(shè)計和制造等方面仍存在著種種問題。</p><p>　　1.1.3 齒輪減速器的分類</p><p>　　齒輪減速器[3]按減速齒輪的級數(shù)可分為單級齒輪減速器、二級齒輪減速器、三級齒輪減速器和多級齒輪減速器幾種；按軸在空間的相互配置方式可分為立式和臥式減速器兩種；按運動簡圖的特點可分為展開式、同軸式和分流式減速器等。</p><p>　　單級

23、圓柱齒輪減速器的最大傳動比一般為8—10，作此限制主要為避免外廓尺寸過大。若要求i>10時，就應采用二級圓柱齒輪減速器。二級圓柱齒輪減速器應用于i為8—50及高、低速級的中心距總和為250—400mm的情況下。三級圓柱齒輪減速器，用于要求傳動比較大的場合。圓錐齒輪減速器和二級圓錐—圓柱齒輪減速器，用于需要輸入軸與輸出軸成90度配置的傳動中。因大尺寸的圓錐齒輪較難精確制造，所以圓錐—圓柱齒輪減速器的高速級總是采用圓錐齒輪傳動以減小其

24、尺寸，提高制造精度。</p><p>　　1.1.4 齒輪減速器的結(jié)構(gòu)</p><p>　　絕大多數(shù)減速器的箱體是用中等強度的鑄鐵鑄成[4]，重型減速器用高強度鑄鐵或鑄鋼。少量生產(chǎn)時也可以用焊接箱體。鑄造或焊接箱體都應進行時效或退火處理。箱體通常由箱座和箱蓋兩部分所組成，其剖分面則通過傳動的軸線。</p><p>　　對于兩軸系結(jié)構(gòu)[5]，由于采用直齒圓柱齒輪，不受

25、軸向力，因此兩軸均由滾動軸承支承。軸向位置由端蓋確定，而端蓋嵌入箱體上對應槽中。為了避免積累誤差過大，保證裝配要求，軸上各裝有一個調(diào)整環(huán)。</p><p>　　油面觀察結(jié)構(gòu)是通過油面指示片上透明玻璃的刻線，可看到油池中儲油的高度的結(jié)構(gòu)。當儲油不足時，應加油補足，保證齒輪的下部浸入油內(nèi)，從而滿足齒輪嚙合和軸承的潤滑。并且油標用于檢查油面高度，這樣可以保證有正常的油量，箱體上安裝油面指示片結(jié)構(gòu)的螺孔不能鉆通，避免機油

26、向外滲漏。</p><p>　　油封裝置是軸從透蓋孔中伸出，該孔與軸之間留有一定間隙的結(jié)構(gòu)。為了防止油向外滲漏和灰塵進入箱體內(nèi)，端蓋內(nèi)裝有毛氈密封圈，此圈應緊緊套在軸上。</p><p>　　對于透氣裝置，當減速器工作時，由于磨擦而產(chǎn)生熱，箱體內(nèi)溫度就會升高而引起揮發(fā)氣體熱膨脹，導致箱體內(nèi)壓力增高。因此，在頂部設(shè)計有透氣裝置，通過通氣塞的小孔使箱體內(nèi)的熱量能夠排出，從而避免箱體內(nèi)的壓力增高

27、。</p><p>　　軸套用于齒輪的軸向定位，它是空套在軸上的，因此內(nèi)孔應大于軸徑。齒輪端面必須超出軸肩，以確定齒輪與軸套接觸，從而保證齒輪軸向位置的固定。</p><p>　　放油孔和放油螺塞的作用是為了換油及清洗箱體時排除污油，其螺孔應低于油池底面，以便放盡機油，通常配有封油墊圈。</p><p>　　起吊裝置用于吊運箱蓋、箱座或整個減速器，包括吊環(huán)螺釘、吊耳

28、、吊鉤等。</p><p>　　1.1.5 齒輪減速器的工作原理</p><p>　　一級圓柱齒輪減速器是通過裝在箱體內(nèi)的一對嚙合齒輪的轉(zhuǎn)動[4]，動力從一軸傳至另一軸，實現(xiàn)減速的。動力由電動機通過皮帶輪傳送到齒輪軸，然后通過兩嚙合齒輪（小齒輪帶動大齒輪）傳送到軸，從而實現(xiàn)減速之目的。</p><p>　　減速器有兩條軸系——兩條配線，兩軸分別由滾動軸承支撐在箱體上

29、，采用過渡配合有較好的同軸度，從而保證齒輪嚙合的穩(wěn)定性。箱體采用分離式，沿軸線平面分為箱座和箱蓋，二者采用螺栓連接，這樣便于裝修。為了保證箱體上安裝軸承和端蓋的孔的正確形狀，兩個零件是在一起加工的。裝配時，他們之間采用兩銷定位，銷孔做成通孔。</p><p><b>　　1.2優(yōu)化設(shè)計</b></p><p>　　優(yōu)化設(shè)計是20世紀60年代初發(fā)展起來的一門新科學[5]

30、，它是將最優(yōu)化原理和計算技術(shù)應用于設(shè)計領(lǐng)域，能從眾多的設(shè)計方案中尋找盡可能完善的或最為適宜的設(shè)計方案的先進的設(shè)計方法，為工程設(shè)計提供一種重要的科學設(shè)計方案，從而大大提高設(shè)計效率和質(zhì)量。</p><p>　　一項機械產(chǎn)品的設(shè)計[6]一般是要經(jīng)過調(diào)查分析、設(shè)計方案擬定、技術(shù)設(shè)計、圖紙的繪制等環(huán)節(jié)。傳統(tǒng)設(shè)計方法通常是以調(diào)查分析為基礎(chǔ)的，參照相似或者同類的產(chǎn)品進行估算、經(jīng)驗的類比或?qū)嶒瀬泶_定初始設(shè)計方案。然后根據(jù)初始設(shè)計

31、方案的設(shè)計參數(shù)進行強度、剛度、穩(wěn)定性等性能的分析計算，檢查各性能是否滿足設(shè)計指標要求。很多的情況下，會出現(xiàn)工作人員需要對參數(shù)進行修改，往往是憑借經(jīng)驗或是直觀判斷，并不是根據(jù)某種理論精確計算出來的。所以在機械產(chǎn)品的設(shè)計的過程中仍有很大的改進提高的余地。</p><p>　　近年來發(fā)展起來的計算機輔助設(shè)計結(jié)合優(yōu)化設(shè)計的思想，使得能夠在優(yōu)化設(shè)計中不斷的選擇設(shè)計參數(shù)并評選出最優(yōu)的設(shè)計方案，相比于傳統(tǒng)的方法不僅可以加快設(shè)計

32、速度也能夠縮短設(shè)計周期。隨著科學技術(shù)的發(fā)展，也就相應的要求機械產(chǎn)品更新周期縮短，結(jié)合優(yōu)化設(shè)計方法與計算機輔助設(shè)計，使整個過程可以自動化，這也成為設(shè)計方法的一個重要發(fā)展趨勢。</p><p>　　減速器的優(yōu)化設(shè)計的方法主要包括以下幾種方法：第一、是利用PRO/E軟件中提供了強大的數(shù)控編程模塊Pro/NC進行優(yōu)化；第二、通過智能算法優(yōu)化，比如遺傳算法，粒子群算法，蟻群算法等；第三、通過MATLAB或者C語言進行優(yōu)化。

33、通過上面的三個方面可以進行優(yōu)化，能夠使實際設(shè)計效率得到明顯的提高，并獲得最優(yōu)解。</p><p>　　1) 參考文獻[7]和[8]應用Pro/E對一級圓柱齒輪減速器進行參數(shù)化設(shè)計，首先通過軟件建立數(shù)學模型，裝配減速器，模擬實際的減速器的運動狀態(tài)，對裝配之后的輸出軸進行有限元的分析，判斷軸與齒輪、軸承相互作用后的應力的分布情況，得出相應的參數(shù)和結(jié)論。這種方法不僅為減速器的數(shù)控加工的編程提供了很好的參數(shù)模型，也為其參

34、數(shù)的設(shè)計打下了基礎(chǔ)。</p><p>　　2) 參考文獻[9]和[10]是采用粒子群對一級齒輪減速器進行優(yōu)化設(shè)計，基于群體智能的全局隨機尋優(yōu)算法，它通過例子搜尋自身的個體最優(yōu)解和全局最優(yōu)解來對減速器模型進行優(yōu)化，此種方法類似于遺傳算法，最終實現(xiàn)了體積最小的目的。</p><p>　　3) 參考文獻[11]和[12]以一級齒輪減速器體積最小建立了優(yōu)化設(shè)計的數(shù)學模型，根據(jù)齒輪減速器本身的約

35、束條件，采用了Matlab軟件中的優(yōu)化工具箱進行優(yōu)化設(shè)計，整個過程初始參數(shù)輸入簡單，編程工作量小，提高了設(shè)計的質(zhì)量和和效率，同時也降低了減速器的制造成本。</p><p>　　4) 參考文獻[13]和[14]可知，遺傳算法可以對一級齒輪減速器進行優(yōu)化設(shè)計，建立數(shù)學模型，并結(jié)合齒寬系數(shù)，應力，彎曲強度等約束條件，經(jīng)過交叉、變異、選擇反復的進行迭代因而得出優(yōu)化后的解。</p><p>　　1

36、.3本課題的主要內(nèi)容和解決問題</p><p>　　針對傳統(tǒng)齒輪減速器出現(xiàn)的設(shè)計體積過大、成本增加、效率過低的問題，有必在此基礎(chǔ)上對齒輪減速器模型進行優(yōu)化設(shè)計。</p><p>　　基于優(yōu)化設(shè)計中的各種方法，對齒輪減速器進行優(yōu)化，首先針對一級齒輪減速器，建立數(shù)學模型，以減速器的軸和齒輪的體積最小為目標函數(shù)，結(jié)合優(yōu)化設(shè)計中的具體方法和實際模型的自身約束進行優(yōu)化，本課題主要是研究基于遺傳算法的

37、齒輪減速器的優(yōu)化設(shè)計，在確定遺傳策略時，采用代數(shù)交叉、高斯變異和可行性規(guī)則選擇進行一代一代的優(yōu)化，最終得到最優(yōu)解，所要解決的主要問題是基于多約束函數(shù)的問題在遺傳算法中的處理，課題的最后是對遺傳算法得出的結(jié)果進行機械上的校核以滿足強度上的要求，主要涉及齒輪和軸的強度的校核，并繪制相應的CAD圖紙。</p><p>　　2 齒輪減速器數(shù)學模型的建立</p><p>　　在分析一個問題時通常并

38、不直接的分析這個系統(tǒng)[6]，往往要引入數(shù)學模型的概念，建立一個數(shù)學模型，便是找出問題中的相關(guān)的因素去定義和描述相應的問題。對于本課題中的齒輪減速器的優(yōu)化設(shè)計，首先就是要求將設(shè)計問題按規(guī)定的格式建立數(shù)學模型，選擇合適的優(yōu)化方法及程序，然后再通過計算機的計算，自動地獲得最優(yōu)設(shè)計方法，因此有必要建立齒輪減速器的數(shù)學模型來進行后續(xù)的程序優(yōu)化。</p><p>　　2.1 齒輪減速器的數(shù)學模型</p><

39、;p>　　隨著機械制造技術(shù)的不斷進步和提高，產(chǎn)品的輕型化、小型化道路已經(jīng)成為當今制造業(yè)的發(fā)展趨勢，對于齒輪減速器，在傳動比和功率相等的情況下，其體積小將會降低成本、節(jié)約資源等優(yōu)勢，因此在對齒輪減速器進行優(yōu)化設(shè)計時，可以按它們的體積之和為最小的原則來建立優(yōu)化設(shè)計的目標函數(shù)。</p><p>　　2.1.1 建立一級齒輪減速器的數(shù)學模型</p><p>　　如圖2.1所示是一級圓柱直齒輪

40、減速器的結(jié)構(gòu)簡圖，已知齒數(shù)比為u=3,輸入功率為P=100KW，主動齒輪的轉(zhuǎn)速為n1=1500 r/min，由于齒輪和軸的尺寸(即殼體內(nèi)的零件)是決定減速器體積大小的因素[12]，利用齒輪幾何尺寸及機構(gòu)尺寸的計算公式，殼體內(nèi)的齒輪和軸的體積可近似地表示為： </p><p>　　圖2-1一級圓柱齒輪減速器的結(jié)構(gòu)尺寸</p><p>　　小齒輪的體積：V1 =0.25π(d12-dz12)b

41、 </p><p>　　高速軸的體積：V2 =0.25πdz12 *(48+l)b</p><p>　　大齒輪的體積：V3 =0.25π(d22-dz22)- 0.25π(Dg22-dg22)(b2-c)-πd02c</p><p>　　低速軸的體積：V4 =0.25πdz22 *(96+l) </p><p>　　軸和齒輪的總體積：V =

42、V1+V2+V3+V4 V=0.25πb1(d12-dz12)+0.25πb2(d22-dz22)-0.25π(Dg22-dg22)(b2-c)-πd02c+0.25πl(wèi)(dz12+dz22)+12πdz12 +24πdz22</p><p>　　=0.25π(m2z12b-dz12b+m2z12u2b-dz22b-0.8b(mz1u-10m)2+2.05bdz12-0.05

43、b(mz1u-10m-1.6dz2)2+</p><p>　　48dz12 +96dz22 )；</p><p>　　上式中的各個符號的含義可由減速器的經(jīng)驗的結(jié)構(gòu)尺寸直接給出[15]，其計算公式為：</p><p>　　d1 = mz1 ； 式(2.1)

44、</p><p>　　d2 = mz2 ； 式(2.2)</p><p>　　Dg2 = umz1 -10m； 式(2.3)</p><p>　　dg2

45、 = 1.6dz2 ； 式(2.4)</p><p>　　d0 = 0.25(umz1 -10m-1.6dz2) ； 式(2.5)</p><p>　　c = 0.2b ；

46、 式(2.6)</p><p>　　由上式可知，當齒數(shù)比給定后，體積V取決于b, z1, m, l, dz1, dz2 這六個參數(shù)，則優(yōu)化設(shè)計的變量可取為：</p><p>　　x= [x1 x2 x3 x4 x5 x6]T = [ b z1 m l dz1 dz

47、2 ] T</p><p>　　將上面所設(shè)定的變量，從x1到x6 代入到上面所示的齒輪減速器的數(shù)學模型的表達式中，得出：</p><p>　　V=0.25*pi*(2.35*x(1)*x(2)2*x(3)2-85*x(1)*x(3)2+0.922*x(1)*x(6)2+51*x(1)*x(2)*x(3)2+x(4)*x(6)2+28*x(5)2+32*x(6)2-x(1)*x(5)2+0.

48、48*x(1)*x(2)*x(3)*x(6)-1.6*x(1)*x(3)*x(6))；</p><p>　　2.1.2 確定約束函數(shù)</p><p>　　為了避免產(chǎn)生根切現(xiàn)象，對于α=20o的標準直齒圓柱齒輪，應有齒數(shù)z1 應該大于zmin，即：</p><p>　　g1(x)= zmin - z1 ≤0</p><p>　　并且此閉式齒輪

49、減速器一般轉(zhuǎn)速較高，為了提高傳動的平穩(wěn)性，減小沖擊的振動，則齒數(shù)多一些為好。</p><p>　　齒輪的強度計算公式可知，齒輪越寬，承載能力也越高，因而齒輪不宜過載，但是增大齒寬又會使得齒面上的載荷分布更趨于不均勻，故齒寬系數(shù)應該取得適當，齒寬系數(shù)應滿足Фmax ≤b/d ≤Ф min，可參見下表3-3[1]。 </p><p>　　表3-3圓柱齒輪的齒寬系數(shù)Фd</p>

50、<p>　　故可以取Фmax =1.4, Ф min=0.9 ，即:</p><p>　　g2(x)= Фmin - b/mz1 ≤0</p><p>　　g3(x)= b/mz1- Фmax ≤0</p><p>　　對于傳遞動力的齒輪，模數(shù)不能過小，一般是m≥2，可見下表3-4所示的圓柱齒輪標準模數(shù)系列表[1]。</p><p&

51、gt;　　表3-4圓柱齒輪標準模數(shù)系列表（GB/T1357-1987）</p><p>　　即：g4(x)= 2 - m ≤0</p><p>　　根據(jù)工藝裝備條件，限制大齒輪的直徑d2不超過1500mm，故小齒輪的直徑d1 不應該超過500mm。</p><p>　　即：g5(x)= mz1 –500 ≤0</p><p>　　按照結(jié)構(gòu)的關(guān)

52、系，軸的支撐跨距應滿足：l ≥b+2△min +0.5 dz2 ，其中△為箱體內(nèi)壁到軸承中心線的距離，現(xiàn)去△min=20mm。</p><p>　　即：g6(x)= b+40 +0.5 dz2 –l ≤0</p><p>　　考慮齒輪的接觸疲勞強度條件進行考慮約束。</p><p>　　1)選擇齒輪材料和熱處理方法</p><p>　　對于正

53、常使用的一般的齒輪減速器，故大齒輪和小齒輪均選用45號鋼，為了制造的方便，可以采用軟齒面，小齒輪調(diào)質(zhì)處理，大齒輪正火處理，選用8級精度。</p><p>　　根據(jù)下面的表3-5常用齒輪材料及其力學特性可查出，小齒輪齒面的硬度為217~255HBW，可以取硬度值為240HBW進行計算；大齒輪的齒面硬度為162~217HBW，可以取硬度值為200HBW進行計算。</p><p>　　表3-5常

54、用齒輪材料及其力學特性</p><p>　　2）齒輪的許用接觸應力</p><p>　　齒輪的許用接觸應力見式(2.7)。</p><p><b>　　式(2.7)</b></p><p>　　式中 SH —— 表示疲勞強度安全系數(shù)；</p><p>　　Z——表示考慮應力循環(huán)次數(shù)影響的系數(shù)

55、，稱為壽命系數(shù)；</p><p>　　σH lim ——表示齒輪的疲勞極限；</p><p>　　查找《機械設(shè)計》圖10-21 齒輪的接觸疲勞強度的極限σHlim，可以查出：</p><p>　　σHlim1 = 580MPa ；</p><p>　　σHlim2 = 400MPa ；</p><p>　　通常情況

56、下[1]，可以取[σ]H=0.9σH lim，即有：</p><p>　　[σ]H1=0.9σH lim1=522MPa；</p><p>　　[σ]H2=0.9σHlim2=360MPa。</p><p>　　對于齒面接觸疲勞強度應該選取較小的值，即可以取[σ]H=360MPa。</p><p><b>　　3）載荷系數(shù)k<

57、/b></p><p>　　選用載荷系數(shù)k=1.4。</p><p>　　4)節(jié)點區(qū)域系數(shù)ZH </p><p>　　對于標準圓柱直齒輪，節(jié)點的區(qū)域系數(shù)ZH =2.5。</p><p><b>　　5）彈性系數(shù)ZE </b></p><p>　　根據(jù)下表3-6彈性影響系數(shù)所示，可以看出ZE

58、 =189.8MPa1/2。 。</p><p>　　表3-6 彈性影響系數(shù)ZE MPa1/2。</p><p>　　6）小齒輪軸上傳遞的扭矩T1</p><p><b>　　式(2.8)</b></p><p>　　式中 P ——傳遞功率；</p>

59、<p>　　n1 ——小齒輪的轉(zhuǎn)速；</p><p>　　又因為齒面接觸疲勞強度見式(2.9)。</p><p><b>　　式(2.9)</b></p><p>　　式中 T1 ——軸所受到的扭矩，T=T1；</p><p><b>　　k ——載荷系數(shù)；</b></p>

60、<p>　　ZH ——節(jié)點區(qū)域系數(shù)；</p><p>　　ZE ——彈性系數(shù)；</p><p><b>　　u ——齒數(shù)比；</b></p><p><b>　　b ——齒寬；</b></p><p><b>　　d——分度圓直徑；</b></p>

61、<p>　　根據(jù)上面所得出的數(shù)據(jù)代入上面的公式可得：</p><p>　　考慮齒輪的齒根彎曲疲勞強度條件進行考慮約束。</p><p>　　1）齒形系數(shù)YF和應力校正系數(shù)YS </p><p>　　查看表3-7和《齒輪減速器精解》圖8-8和圖8-9[4]</p><p>　　表3-7 齒形系數(shù)YF以及應力校正系數(shù)YS</p&g

62、t;<p>　　根據(jù)上表則可以得出：</p><p>　　YF1 =2.53 YF2 =2.22；</p><p>　　YS1 =1.62 YS2 =1.81 </p><p>　　2）齒輪的許用彎曲應力見式(2.10)。</p><p><b>　　式(2.10)</b&g

63、t;</p><p>　　根據(jù)《機械設(shè)計》查出彎曲疲勞極限應力為[1]：</p><p>　　σFlim1 = 220MPa</p><p>　　σFlim2 = 170MPa</p><p>　　由《減速器設(shè)計精解》圖8-11查出壽命系數(shù) YN1 = YN2 =1，由表8-20查出安全系數(shù)SF =1.25，故：</p>&

64、lt;p>　　所以可以根據(jù)齒根彎曲疲勞強度的公式，見式(2.11)。</p><p><b>　　式(2.11)</b></p><p><b>　　即得出：</b></p><p>　　按照軸的彎曲強度條件考慮，根據(jù)公式進行考察見式(2.12)。</p><p><b>　　式(2

65、.12)</b></p><p>　　式中 T ——軸所受到的扭矩，T=T1；</p><p>　　M——軸所受到的彎矩，單位是N.mm，又可根據(jù)公式M=Fn l /2，因為齒輪所受到的法向力Fn=2T1/d1cosα 其中對于標準的直齒圓柱齒輪，α=20o，即可以得出 T1l/(mz1cosα)= 677305 l/mz1；</p><p>　　

66、a——折合系數(shù)，是用來考慮彎矩產(chǎn)生的彎曲應力和由扭矩產(chǎn)生的扭轉(zhuǎn)剪切力的循環(huán)系數(shù)特性的不同。在這里可以取a=0.59；</p><p>　　[σ]b——軸的許用彎曲應力，根據(jù)《齒輪減速器精解》表8-32[4]查得軸的許用彎曲應力 [σ] b =60MPa。</p><p>　　表3-7 軸的許用彎曲應力 單位：MPa</p><p>　　W

67、——軸的抗彎截面系數(shù)，單位mm3 ，根據(jù)[1]《機械設(shè)計》可查出實心軸的彎曲截面系數(shù)W=0.1dz 3 。</p><p>　　因此可以根據(jù)上面所得到的公式和數(shù)據(jù)代入，即可得到：</p><p>　　綜上可知，一級圓柱直齒輪減速器以體積最小為優(yōu)化目標的優(yōu)化設(shè)計問題，是一個具有11個不等式約束的優(yōu)化問題，將x= [x1 x2 x3 x4 x5 x6]T = [ b

68、 z1 m l dz1 dz2 ] T代入其數(shù)學模型可以表示為：</p><p>　　Vmin=0.25*π*(2.35*x(1)*x(2)2*x(3)2-85*x(1)*x(3)2+0.922*x(1)*x(6)2+51*x(1)*x(2)*x(3)2+x(4)*x(6)2+28*x(5)2+32*x(6)2-x(1)*x(5)2+0.48*x(1)*x(2)*x(3)*x(6)-1.6*x(

69、1)*x(3)*x(6))。</p><p>　　化簡后的約束函數(shù)是：</p><p>　　g1(x)= 17- x(1) ≤0</p><p>　　g2(x)= 0.9 - x(1)/ ( x(3)x(2) )≤0</p><p>　　g3(x)= x(1)/ (x(3) x(2))- 1.4 ≤0</p><p>

70、　　g4(x)= 2 - x(3) ≤0</p><p>　　g5(x)= x(3) x(2) – 500 ≤0</p><p>　　g6(x)= x(1)+40 +0.5 x(6) –l ≤0</p><p>　　g7(x)= 731544.36/( x(2) x(3)x11/2) -360 ≤0</p><p>　　g8(x)= 7306

71、438/(x1 x(2) x(3)2) -176 ≤0</p><p>　　g9(x)= 2387703.86/( x1 x(2) x(3)2 ) -136 ≤0</p><p>　　g10(x)= [77305 l/ x(3) x(2)）2 + 375633.352 ] 2 / 0.1 x(5) 3 -60 ≤0</p><p>　　g11(x)= [225768

72、.3 l/ x(3) x(2)）2 + 375633.352 ] 2 / 0.1 x(6)3 -60 ≤0</p><p>　　3 遺傳算法程序設(shè)計</p><p>　　上一章已建立了數(shù)學模型，下面的步驟則是選用合適的優(yōu)化方法和計算機程序?qū)ζ溥M行自動的計算，對此所采用的優(yōu)化方法有MATLAB優(yōu)化、PRO/E軟件優(yōu)化、智能優(yōu)化方法等，本課題采用的是基于智能優(yōu)化中的遺傳算法的優(yōu)化設(shè)計。<

73、;/p><p>　　3.1 智能算法選擇</p><p>　　智能優(yōu)化方法中的智能計算也有人稱之為“軟計算”，是人們受自然（生物界）規(guī)律的啟迪[16]，根據(jù)其原理，模仿求解問題的算法。從自然界得到啟迪，模仿其結(jié)構(gòu)進行發(fā)明創(chuàng)造，這就是仿生學。這是我們向自然界學習的一個方面。另一方面，我們還可以利用仿生原理進行設(shè)計(包括設(shè)計算法)，這就是智能計算的思想。這方面的內(nèi)容很多，如人工神經(jīng)網(wǎng)絡(luò)技術(shù)、遺傳算

74、法、模擬退火算法和群集智能技術(shù)等。</p><p>　　3.1.1 智能算法種類</p><p>　　遺傳算法是屬于進化算法的一種,它通過模仿自然界的選擇與遺傳的機理來尋找最優(yōu)解[17]。遺傳算法有三個基本算子：選擇、交叉和變異。 但是遺傳算法的編程實現(xiàn)比較復雜，首先需要對問題進行編碼，找到最優(yōu)解之后還需要對問題進行解碼，另外三個算子的實現(xiàn)也有許多參數(shù)，如交叉概率和變異概率，并且這些參數(shù)的

75、選擇嚴重影響解的品質(zhì)，而目前這些參數(shù)的選擇大部分是依靠經(jīng)驗。</p><p>　　粒子群優(yōu)化算法是一種進化計算技術(shù)[10]，它是來自對鳥群捕食的行為的研究。粒子群算法和遺傳算法相似，也是一種基于迭代的算法。系統(tǒng)初始化一組隨機的解，通過反復的迭代搜尋最優(yōu)的值。相比于遺傳算法，它沒有交叉以及變異的過程，而是粒子在解空間追隨最優(yōu)的粒子進行搜索。目前粒子群的優(yōu)化算法已廣泛應用于函數(shù)優(yōu)化，神經(jīng)網(wǎng)絡(luò)訓練，模糊系統(tǒng)控制以及其他

76、遺傳算法的應用領(lǐng)域。</p><p>　　蟻群算法是一種用來在圖中尋找優(yōu)化路徑的機率型算法[5]。它來自螞蟻在尋找食物過程中發(fā)現(xiàn)路徑的行為。蟻群算法是一種模擬進化算法，通過將蟻群算法設(shè)計的結(jié)果與遺傳算法設(shè)計的結(jié)果進行比對發(fā)現(xiàn)，蟻群算法具有一種新的模擬進化優(yōu)化方法的應用價值。</p><p>　　對于本次的齒輪減速器優(yōu)化設(shè)計，采用遺傳算法進行優(yōu)化，遺傳算法還具有以下幾方面的優(yōu)越性[17]：首

77、先，它在搜索過程中不容易陷入局部最優(yōu)，即使在所定義的適應函數(shù)是不連續(xù)的、非規(guī)則的或有噪聲的情況下，它也能以很大的概率找到整體最優(yōu)解；其次，由于它固有的并行性，遺傳算法非常適用于大規(guī)模并行計算機。 </p><p>　　3.2 遺傳算法介紹</p><p>　　遺傳算法是模擬達爾文生物進化論的自然選擇和遺傳學機理的生物進化過程的計算模型[17]，是一種通過模擬自然進化過程搜索最優(yōu)解的方法，它

78、最初由美國Michigan大學J.Holland教授于1975年首先提出來的，并出版了頗有影響的專著 《Adaptation in Natural and Artificial Systems》，GA這個名稱才逐漸為人所知，J.Holland教授所提出的GA通常為簡單遺傳算法（SGA）。它是計算數(shù)學中用于解決最優(yōu)化的搜索算法，是進化算法的一種。進化算法最初是借鑒了進化生物學中的一些現(xiàn)象而發(fā)展起來的，這些現(xiàn)象包括遺傳、突變、自然選擇以及交

79、叉等。</p><p>　　圖3-1 一般遺傳算法的流程圖</p><p>　　遺傳算法通常實現(xiàn)方式為一種計算機模擬。按照遺傳算法的工作流程，當用遺傳算法求解問題時，必須在目標問題實際表示與遺傳算法的染色體位串結(jié)構(gòu)之間建立聯(lián)系，也就是確定編碼和解碼的運算，本次的畢業(yè)設(shè)計采用的是實數(shù)編碼方式。然后定義適應度函數(shù)，確定遺傳策略，包括選擇群體的大小n、選擇、交叉、變異方法，以及確定交叉概率、變異

80、概率等遺傳參數(shù)。則可隨機初始化生成初始種群，計算群體中個體的適應度值，然后則運用交叉、變異、選擇遺傳算子進行作用于群體，產(chǎn)生下一代新的種群，并調(diào)用約束函數(shù)來判斷是否滿足來確定是否可以停止種群。</p><p>　　3.3 遺傳算法程序?qū)崿F(xiàn)</p><p>　　3.3.1 遺傳算法中的參數(shù)</p><p>　　染色體[18] ——可以叫做基因型個體(individua

81、ls)，一定數(shù)量的個體組成了群體(population)，群體中個體的數(shù)量叫做群體大小。</p><p>　　適應度函數(shù)(fitness function) ——就是問題中的全體對象與其適應度之間的一個對應關(guān)系，對象集合到適應度集合的一個映射。它一般是定義在論域空間上的一個實數(shù)值函數(shù)。 </p><p>　　種群(population) —— 是模擬生物種群而由若干個染色體組成的群體，它一

82、般是整個區(qū)域空間的一個很小的子集。遺傳算法就是通過在種群上實施所稱的遺傳操作，使其不斷更新?lián)Q代而實現(xiàn)對整個區(qū)域空間的搜索。</p><p>　　交叉概率(crossover rate) ——就是參加交叉運算的染色體個數(shù)占全體染色體總數(shù)的比例，記為Pc，取值范圍一般為0.4～0.99。由于生物繁殖時染色體的交叉是按一定的概率發(fā)生的，因此參加交叉操作的染色體也有一定的比例。本課題中交叉概率為0.6。</p>

83、;<p>　　變異概率(mutation rate) ——是指發(fā)生變異的基因位數(shù)所占全體染色體的基因總位數(shù)的比例，記為Pm，取值范圍一般為0.0001～0.1。由于在生物的繁衍進化過程中，變異也是按一定的概率發(fā)生的，而且發(fā)生概率一般很小。本課題中所取的變異概率為0.02。</p><p>　　3.3.2 遺傳算法的偽代碼</p><p>　　input: GA paramet

84、ers</p><p>　　output: best solution 輸出最優(yōu)解</p><p><b>　　begin</b></p><p>　　t 0; //進化種群代數(shù)</p><p>　　initialize P(t) by real number encoding; // 通過實數(shù)編碼進行初始化種群&

85、lt;/p><p>　　fitness eval(P);</p><p>　　while (not termination condition) do //不滿足終止條件時，循環(huán)</p><p>　　crossover P(t) to yield C(t) by arithmetic crossover; //交叉</p><p>　　mut

86、ation P(t) to yield C(t) by nonuniform mutation; // 變異</p><p>　　fitness eval(C); // </p><p>　　select P(t+1) from P(t) and C(t) by top popSize selection; //選擇</p><p>　　t t + 1; &

87、lt;/p><p><b>　　end</b></p><p>　　output best solution; //輸出最優(yōu)解</p><p><b>　　end // 結(jié)束</b></p><p>　　3.3.3 編碼和解碼</p><p>　　編碼就是把一個問題的可行解從其

88、解空間轉(zhuǎn)換到遺傳算法所能處理的搜索空間的轉(zhuǎn)換方法[17]。它將決定個體的基因排列形式，并決定個體從搜索空間的基因型變換到解空間的表現(xiàn)型時的解碼方法以及會影響交叉算子、變異算子等遺傳算子的運算方法。</p><p>　　二進制編碼方法是遺傳算法中最常用的一種編碼方法。使用的編碼符號集是由符號0和1所組成的，所構(gòu)成的個體基因型是一個二進制編碼符號串。對于本課題的遺傳算法的編碼采用浮點數(shù)編碼，個體的每個基因值用某一范圍

89、內(nèi)的一個浮點數(shù)來表示，個體的編碼長度等于其決策變量的個數(shù)。使用的是決策變量的真實值。 </p><p>　　采用浮點編碼應該要注意[17]，保證基因值在給定的區(qū)間限制范圍內(nèi)，使用遺傳操作所產(chǎn)生的新個體的基因值在給定的區(qū)間限制范圍內(nèi)，當用多個字節(jié)表示一個基因值時，交叉運算必須在兩個基因的分界字節(jié)處進行，而不能在某個基因的中間字節(jié)分隔處進行。</p><p>　　浮點編碼適合于在遺傳算法中表示

90、的范圍較大的情況，適合于精度要求較高的遺傳優(yōu)化問題，便于較大空間的遺傳搜索，它可以改善了遺傳算法的計算復雜性，提高了運算效率，便于遺傳算法與經(jīng)典優(yōu)化方法的混合使用，便于設(shè)計針對問題的專門知識的知識型遺傳操作，便于處理復雜的決策變量約束條件。</p><p><b>　　浮點編碼的程序：</b></p><p>　　function Ret=Code(LenChrom,

91、Opts,Bound)</p><p>　　case 'float' % float coding</p><p>　　Pick=rand(1,length(LenChrom));</p><p>　　Ret=[170 20 8 350 70 150]+5*Pick;</p><p>　　%Ret=Bound

92、(:,1)'+(Bound(:,2)-Bound(:,1))'.*Pick;</p><p><b>　　End</b></p><p><b>　　浮點解碼的程序：</b></p><p>　　function Ret=Decode(LenChrom,Bound,Code,Opts)</p>

93、<p><b>　　Ret=Code;</b></p><p><b>　　end</b></p><p><b>　　3.3.4 交叉</b></p><p>　　交叉就是互換兩個染色體某些位上的基因[17]。交叉算子的設(shè)計包括的內(nèi)容：確定交叉點的位置和確定部分基因的交換方式。</

94、p><p>　　交叉一般包括兩種交叉方式：單點交叉和算術(shù)交叉，其中算數(shù)交叉的操作對象是由浮點數(shù)編碼表示的個體，指由兩個個體的線性組合而產(chǎn)生出兩個新的個體。</p><p>　　執(zhí)行過程是：先確定兩個個體進行線性組合時的系數(shù) α，再按式(3.1)產(chǎn)生兩個新的子代個體。</p><p><b>　　式(3.1)</b></p><p

95、>　　交叉的程序的代碼如下：</p><p>　　function Ret=Cross(PCross,LenChrom,Individuals,SizePop,Opts,Pop)</p><p>　　for i=1:SizePop</p><p>　　% select two children at random</p><p>　　P

96、ick=rand(1,2);</p><p>　　while prod(Pick)==0</p><p>　　Pick=rand(1,2);</p><p><b>　　end</b></p><p>　　Index=ceil(Pick.*SizePop);</p><p>　　% random

97、position of crossover</p><p>　　Pick=rand;</p><p>　　while Pick==0</p><p>　　Pick=rand;</p><p><b>　　end</b></p><p>　　if Pick>PCross</p>

98、<p><b>　　continue;</b></p><p><b>　　end</b></p><p>　　% random position of crossover</p><p>　　Pick=rand;</p><p>　　while Pick==0</p>&l

99、t;p>　　Pick=rand;</p><p><b>　　end</b></p><p>　　Pos=ceil(Pick.*sum(LenChrom));</p><p>　　Pick=rand;</p><p>　　V1=Individuals.Chrom(Index(1),Pos);</p>

100、<p>　　V2=Individuals.Chrom(Index(2),Pos);</p><p>　　Individuals.Chrom(Index(1),Pos)=Pick*V2+(1-Pick)*V1;</p><p>　　Individuals.Chrom(Index(2),Pos)=Pick*V1+(1-Pick)*V2;</p><p><

101、;b>　　end</b></p><p>　　Ret=Individuals.Chrom;</p><p><b>　　end</b></p><p><b>　　3.3.5 變異</b></p><p>　　遺傳算法過程中的變異是將個體染色體編碼串中的某些基因座上的基因值用該基因

102、座的其他等位基因來替代[17]，從而形成一個新的個體。其目的是為了改善遺傳算法的局部搜索能力，維持群體的多樣性，防止出現(xiàn)早熟現(xiàn)象。</p><p>　　變異算子的設(shè)計包括的內(nèi)容是：確定變異點的位置和確定基因值替換方法。變異的方法有兩種：單點變異和高斯變異。</p><p>　　單點變異是針對二進制編碼，是等概率的隨機確定一個基因座，以概率改變個體分量的取值。</p><

103、p>　　高斯變異是指進行變異操作時，用符合均值為u，方差為 σ2的正態(tài)分布的一個隨機數(shù)來替換原有基因值從而產(chǎn)生一個新的個體。高斯變異的操作對象是由浮點數(shù)編碼表示的個體。</p><p>　　執(zhí)行過程是：確定進行變異的個體分量，然后則按式產(chǎn)生一個新的子代個體分量。</p><p>　　本次遺傳算法的變異采用高斯變異，變異程序的代碼如下：</p><p>　　Fu

104、nction Ret=Mutation(PMutation,LenChrom,Individuals,SizePop,Opts,Bound,Pop)</p><p>　　for i=1:SizePop</p><p>　　% select child at random</p><p>　　Pick=rand;</p><p>　　while

105、 Pick==0</p><p>　　Pick=rand;</p><p><b>　　end</b></p><p>　　Index=ceil(Pick*SizePop);</p><p>　　Pick=rand;</p><p>　　if Pick>PMutation</p>

106、<p><b>　　continue;</b></p><p><b>　　end</b></p><p>　　% mutation position</p><p>　　Pick=rand;</p><p>　　while Pick==0</p><p>　　

107、Pick=rand;</p><p><b>　　end</b></p><p>　　Pos=ceil(Pick*sum(LenChrom));</p><p>　　v=Individuals.Chrom(i,Pos);</p><p>　　v1=v-Bound(Pos,1);</p><p>　

108、　v2=Bound(Pos,2)-v;</p><p>　　Pick=rand;</p><p>　　if Pick>0.5</p><p>　　Delta=v2*(1-Pick^((1-Pop(1)/Pop(2))^2));</p><p>　　Individuals.Chrom(i,Pos)=v+Delta;</p>

109、<p><b>　　else</b></p><p>　　Delta=v1*(1-Pick^((1-Pop(1)/Pop(2))^2));</p><p>　　Individuals.Chrom(i,Pos)=v-Delta;</p><p><b>　　end</b></p><p>&

110、lt;b>　　end</b></p><p>　　Ret=Individuals.Chrom;</p><p><b>　　End</b></p><p><b>　　3.3.6 選擇</b></p><p>　　遺傳算法中使用選擇算子來對群體中的個體進行優(yōu)勝劣汰操作，即有下面的原

111、則：適應度較高的個體被遺傳到下一代群體中的概率較大；適應度較低的個體被遺傳到下一代群體中的概率較小。</p><p>　　對于本課題區(qū)別與一般的遺傳算法的地方主要在于要處理多個約束，而處理約束優(yōu)化的關(guān)鍵是如何對可行域邊界進行有效的，適度的，可控的搜索。本次遺傳算法處理約束優(yōu)化的方法是可行性規(guī)則法?；谌齻€準則[9]：1.可行解總優(yōu)于不可行解；2.兩個解都可行，適應度小的為優(yōu)；3.兩個解都不可行，違反約束小的為優(yōu)；

112、</p><p>　　本次選擇的偽代碼如下：</p><p>　　if X1，X2 均滿足約束；</p><p>　　比較X1，X2所對應的體積，越小的越好；</p><p><b>　　結(jié)束 </b></p><p>　　if X1滿足約束條件，X2不滿足約束條件；</p>

113、;<p><b>　　選擇X1；</b></p><p><b>　　結(jié)束</b></p><p>　　If X2滿足約束條件，X1不滿足約束條件；</p><p><b>　　選擇X2 ；</b></p><p><b>　　結(jié)束</b&g

114、t;</p><p>　　if X1，X2 均不滿足約束條件；</p><p>　　比較X1，X2 違反約束的程度，違反約束程度越小越好。</p><p><b>　　結(jié)束</b></p><p>　　在選擇過程中將會調(diào)用名為Panduan的子程序，判斷是否滿足約束條件，并遵循處理約束的三個法則，則會有如下的程序代

115、碼：</p><p>　　function Ret=panduan(x)</p><p>　　Ret=0; %滿足條件</p><p>　　if 17-x(2)>0</p><p>　　Ret=Ret-abs(17-x(2));</p><p><b>　　end</b></p&g

116、t;<p>　　if 0.9-x(1)/(x(2)*x(3))>0</p><p>　　Ret=Ret-abs(0.9-x(1)/(x(2)*x(3)));</p><p><b>　　end</b></p><p>　　if x(1)/(x(2)*x(3))-1.4>0</p><p>　　R